汪文強(qiáng)，張新平，宋振東

(中鐵二院工程集團(tuán)有限責(zé)任公司，四川 成都 610031)

1 引 言

大瑞鐵路高黎貢山隧道是亞洲在建鐵路第一長(zhǎng)隧道，工程地質(zhì)條件復(fù)雜，施工難度大，風(fēng)險(xiǎn)高[1]。其出口段采用TBM施工，施工過程中在斷層破碎帶附近遭遇多次卡機(jī)，嚴(yán)重影響施工進(jìn)度[2]。因此，在TBM隧道施工過程中開展超前地質(zhì)預(yù)報(bào)是必要的，可以為潛在的卡機(jī)段提供依據(jù)進(jìn)行預(yù)處理，從而降低卡機(jī)風(fēng)險(xiǎn)。目前國(guó)內(nèi)外在TBM隧道施工時(shí)采用的超前地質(zhì)預(yù)報(bào)物探方法主要有地震負(fù)視速度法及TSP法(瑞士)、TGP法(國(guó)產(chǎn))、TRT法[3](美國(guó))、TST法(國(guó)產(chǎn))、HSP法[4](國(guó)產(chǎn))等[5]。高黎貢山隧道采用TSP法、李術(shù)才等[6]提出的三維激發(fā)極化法主動(dòng)源、被動(dòng)源三維地震進(jìn)行綜合物探法預(yù)報(bào)。為了研究卡機(jī)段圍巖的TSP法地震波場(chǎng)傳播規(guī)律及對(duì)卡機(jī)段的影響，本文基于Tesseral 2D對(duì)斷層破碎帶附近卡機(jī)段TSP法進(jìn)行數(shù)值模擬[7]，并與實(shí)際卡機(jī)段TSP法預(yù)報(bào)結(jié)果相對(duì)比，旨在觀察地震波場(chǎng)在斷層破碎帶內(nèi)的傳播特征，分析TSP物性參數(shù)[8]，提高預(yù)報(bào)準(zhǔn)確率，為TBM施工預(yù)報(bào)卡機(jī)可能段落提供依據(jù)。

Tesseral 2D是一個(gè)商業(yè)化的全波場(chǎng)模擬軟件包，可以構(gòu)建復(fù)雜地質(zhì)模型。Tesseral 2D基于FDTD有限差分法[9]，這樣可以很快地模擬非常復(fù)雜的地質(zhì)介質(zhì)，包括固體和液體聯(lián)合的介質(zhì)和各向異性介質(zhì)。本文采用Tesseral 2D建立在高黎貢山隧道TBM段采用的地震波反射法TSP模型，模擬TSP法地震波在造成TBM卡機(jī)的斷層破碎帶中的波場(chǎng)傳播特征。

2 卡機(jī)段TSP正演模擬

2.1 TBM卡機(jī)機(jī)理

TBM卡機(jī)包括刀盤被卡和護(hù)盾被卡。通常情況下刀盤被卡主要發(fā)生在斷層破碎帶內(nèi)，采用刀盤瞬時(shí)脫困扭矩進(jìn)行脫困處理。護(hù)盾被卡的原因一般是由于圍巖大變形，導(dǎo)致TBM護(hù)盾被卡，且作用護(hù)盾上的圍巖應(yīng)力引起的阻力超過了TBM脫困推力[2]。高黎貢山隧道TBM掘進(jìn)過程中遭遇了多次刀盤和護(hù)盾卡機(jī)。截至2019年底，累計(jì)發(fā)生17次不同規(guī)模的卡機(jī)事件，卡機(jī)造成的停工時(shí)間短至8天，長(zhǎng)達(dá)數(shù)月，同時(shí)卡機(jī)段落也從10～43 m不等，給施工造成了嚴(yán)重的損失。因此，研究TBM卡機(jī)段TSP法地震波的傳播特征，提高預(yù)報(bào)準(zhǔn)確率，是非常必要的。

2.2 模型建立及參數(shù)設(shè)置

TSP法觀測(cè)系統(tǒng)布置為24個(gè)激發(fā)孔和2個(gè)接受孔，為“多發(fā)一收”模式，本次模擬為減少運(yùn)算時(shí)長(zhǎng)，根據(jù)“炮檢互換”理論，采用“一發(fā)多收”模式，建立24個(gè)接收器和2個(gè)激發(fā)孔。

為了模擬TBM施工過程中遭遇斷層破碎帶卡機(jī)的TSP地震波場(chǎng)傳播特征，本次正演模擬建立長(zhǎng)1 000 m×寬1 000 m的均勻全空間模型模擬隧道開挖所在的圍巖介質(zhì)，巖性為花崗巖。隧道開挖從模型高度500 m處開始，掌子面位于模型566 m 處，隧道高度10 m。造成卡機(jī)的斷層破碎帶位于隧道掌子面前方20 m處，模擬高度100 m，寬度20 m，垂直于隧道走向。模型詳細(xì)物性參數(shù)見表1，軟件建立的模型見圖1。

表1 模型物性參數(shù)

圖1 模型示意圖(橫、縱坐標(biāo)單位：m)Fig.1 Model diagram

圖2 不同時(shí)刻波場(chǎng)快照Fig.2 Snapshot of wave field at different moments

2.3 卡機(jī)段斷層地震波響應(yīng)

本次正演采用聲波方程模擬[10]。圖2為不同時(shí)刻截取的地震波場(chǎng)快照[10]，圖3為24個(gè)檢波器的道集記錄。

選取第24個(gè)檢波器進(jìn)行波場(chǎng)分析。讀取波場(chǎng)快照和道集記錄，結(jié)合模型參數(shù)：山體圍巖縱波波速VP=4 000 m/s，斷層破碎帶縱波波速VP=1 500 m/s。

圖2(a)為地震波第一次到達(dá)24#檢波器，從軟件波場(chǎng)快照和道集記錄上直接讀取時(shí)間值，T1=0.017 s，于是計(jì)算可得到接收器與激發(fā)孔距離為68 m，實(shí)際模擬距離為66 m，誤差為2 m。

圖2(b)為地震波第一次到達(dá)斷層前界面，從軟件波場(chǎng)快照和道集記錄上直接讀取時(shí)間值，T1=0.027 5 s，于是計(jì)算可得到接收器與斷層前界面距離為110 m，實(shí)際模擬距離為106 m，誤差為4 m。

圖2(c)為地震波傳播到斷層前界面后發(fā)生反射，再次傳播到24#接收器的波場(chǎng)快照，從軟件波場(chǎng)快照和道集記錄上直接讀取時(shí)間值，T1=0.037 s，于是計(jì)算可得到地震波的傳播總距離為148 m，實(shí)際模擬地震波的傳播總距離為146 m，誤差為2 m。

圖2(d)為為地震波傳播到斷層后界面的波場(chǎng)快照，從軟件波場(chǎng)快照和道集記錄上直接讀取時(shí)間值，T1=0.058 5 s，亦可得知地震波在斷層破碎帶內(nèi)傳播時(shí)間T2=0.031 s，于是計(jì)算可得到地震波的傳播總距離為156.5 m，實(shí)際模擬地震波的傳播總距離為146 m，誤差為10.5 m，計(jì)算得到斷層破碎帶寬度為46.5 m，實(shí)際模擬斷層破碎帶寬度為40 m，誤差為6.5 m。

圖2(e)為斷層后界面地震波反射傳播至24#接收器的波場(chǎng)快照，從軟件波場(chǎng)快照和道集記錄上直接讀取時(shí)間值，T1=0.092 s，于是計(jì)算可得到地震波的傳播總距離為213 m，實(shí)際模擬地震波的傳播總距離為226 m，誤差為13 m。

圖3 道集記錄Fig.3 gather records

2.4 正演結(jié)果誤差分析

根據(jù)本文第2.2節(jié)數(shù)值模擬結(jié)果，得出模擬計(jì)算誤差與傳播時(shí)間(亦可為傳播距離)的關(guān)系，見圖4。從圖4中可看出，隨著地震波場(chǎng)傳播時(shí)間(傳播距離)增大，數(shù)值模擬的計(jì)算值與實(shí)際建模值誤差有增大趨勢(shì)。分析誤差原因，主要有以下三點(diǎn)：

1)在建模計(jì)算波場(chǎng)快照的時(shí)間步長(zhǎng)值(Δt)偏大，導(dǎo)致時(shí)間讀取不精確。

2)斷層破碎帶模型建立較為簡(jiǎn)單，山體圍巖與斷層過渡不夠線性，以致于產(chǎn)生局部邊界問題，導(dǎo)致地震波場(chǎng)傳播局部出現(xiàn)畸變。同時(shí)斷層破碎帶內(nèi)不同邊界反射波疊加，導(dǎo)致需要觀察的反射波相位讀取變困難。

圖4 誤差分析Fig.4 Error analysis

3)數(shù)值模擬中，山體圍巖和斷層破碎帶內(nèi)縱波速度分別為4 000 m/s、1 500 m/s，速度相差較大，導(dǎo)致斷層計(jì)算寬度與實(shí)際模型誤差較大。

2.5 Tesseral 2D數(shù)值模擬結(jié)果分析

圖5 TSP道集記錄Fig.5 TSP gather records

根據(jù)數(shù)值模擬正演模擬結(jié)果可知，波長(zhǎng)快照可以清楚的反應(yīng)地震波場(chǎng)在斷層破碎帶內(nèi)的傳播規(guī)律，24個(gè)接收器采集的地震道集記錄可以輕易讀取各模型界面的反射信號(hào)，從而計(jì)算得知斷層破碎帶的位置及寬度。因此，利用Tesseral 2D軟件對(duì)TSP法在斷層破碎帶內(nèi)的數(shù)值模擬是可行的，其結(jié)果可以為實(shí)際TSP數(shù)據(jù)處理中反射層提取和計(jì)算提供理論參考。

數(shù)值模擬計(jì)算斷層寬度即傳播距離時(shí)，與實(shí)際建模存在誤差，從本文第2.3節(jié)分析可知，需要調(diào)整波長(zhǎng)快照拍攝時(shí)間步長(zhǎng)，優(yōu)化各模型之間的邊界條件等，從而得出更精確的模擬結(jié)果。不過，有限差分及時(shí)間步長(zhǎng)減小，隨之而來的是幾何上升的數(shù)據(jù)計(jì)算量。

在TBM隧道實(shí)際進(jìn)行TSP超前預(yù)報(bào)的情況極為復(fù)雜，數(shù)值模擬參數(shù)設(shè)置始終無法做到與現(xiàn)場(chǎng)情況完全吻合，這就需要在做數(shù)值模擬時(shí)，盡可能選取與實(shí)際相接近的模型參數(shù)，從而得到更準(zhǔn)確的數(shù)值模擬結(jié)果，為實(shí)際TSP數(shù)據(jù)分析提供參考。

3 實(shí)例

3.1 卡機(jī)段TSP實(shí)施概況

本次TBM卡機(jī)段為燕山期花崗巖地層，巖性為花崗巖，整體呈塊石狀構(gòu)造，弱風(fēng)化。TSP實(shí)施掌子面里程為PDK221+958,TBM卡機(jī)發(fā)生里程為PDK221+918，卡機(jī)影響帶寬度為20 m。該段后期推測(cè)為斷層影響帶，圍巖破碎軟弱。

3.2 TSP預(yù)報(bào)結(jié)果分析

TSP一般采用24個(gè)激發(fā)孔激發(fā)，2個(gè)接收器接收。為了更直觀地對(duì)比實(shí)際采集的TSP道集記錄與數(shù)值模擬道集記錄的差別，對(duì)TSP采集的道集進(jìn)行偏移距計(jì)算、數(shù)據(jù)旋轉(zhuǎn)、時(shí)域?yàn)V波、帶通濾波、初至拾取處理、炮能量均衡，得到圖5處理后的道集記錄。從道集記錄中，可以讀取清晰的直達(dá)波相位，通過增益控制，亦能讀取出反射波相位(圖5)。實(shí)際采集的TSP數(shù)據(jù)由于現(xiàn)場(chǎng)條件復(fù)雜，施工機(jī)械震動(dòng)、激發(fā)孔聲波干擾等，導(dǎo)致信噪比較低，直接從道集記錄上讀取反射波組較為困難。因此，運(yùn)用TSPwin軟件對(duì)數(shù)據(jù)進(jìn)行處理，從2D物性參數(shù)圖和3D云圖分析數(shù)據(jù)，更具有指導(dǎo)意義。

圖6和圖7為TSP數(shù)據(jù)處理后的速度3D云圖和2D物性參數(shù)圖，從圖中可以明顯看到，在PDZK221+926～+900段縱波速度、密度和楊氏模量均下降，推斷該段圍巖破碎，為斷層破碎帶，易造成TBM卡機(jī)。TSP預(yù)報(bào)結(jié)果與實(shí)際揭示及TBM運(yùn)行狀況一致。

圖6 2D物性參數(shù)Fig.6 2D physical property parameter map

4 結(jié) 論

1)利用Tesseral 2D對(duì)TSP法進(jìn)行數(shù)值模擬，觀察地震波在斷層卡機(jī)段傳播規(guī)律，證明了TSP法在TBM卡機(jī)段超前預(yù)報(bào)的實(shí)用性，對(duì)實(shí)際數(shù)據(jù)采集有一定的參考意義；

2)正演模擬中對(duì)斷層寬度計(jì)算與實(shí)際模型存在一定偏差，分析為波場(chǎng)快照提取精度、斷層模型邊界過渡非線性和波阻抗差異過大導(dǎo)致，實(shí)際TSP反演計(jì)算應(yīng)適當(dāng)修正。