TSP斷層模型數值模擬

林 義，劉爭平，王朝令，肖 締

1.西南交通大學地球科學與環境工程學院，成都 610031 2.中鐵西北院科學研究院有限公司，蘭州 730000 3.四川農業大學土木工程學院，成都 611830

?

TSP斷層模型數值模擬

林 義1,2，劉爭平1，王朝令3，肖 締3

1.西南交通大學地球科學與環境工程學院，成都 610031 2.中鐵西北院科學研究院有限公司，蘭州 730000 3.四川農業大學土木工程學院，成都 611830

隧道在施工開挖中會遇到各種地質問題，其中以斷層和軟弱帶居多，目前隧道地質預報主要采用TSP(tunnel seismic prediction)系統進行。雖然TSP技術應用廣泛，但目前對它的研究工作主要集中于工程應用實例，采用正演模擬方法進行的研究較少。筆者采用有限元方法模擬隧道地震波場，采用波場快照與時間記錄相結合的方法研究斷層對隧道地震波場傳播的影響，并對含斷層模型的時間記錄進行了反演處理，得到了數值模型的速度云圖和反射層位圖。數據處理結果表明：采用TSP Win軟件默認值處理得到的速度云圖與模型設定的斷層位置一致；根據反射層位圖，對異常速度帶的層狀模型來說，P波預報的準確性更高。研究表明，TSP系統具有良好的抗噪性能。通過對工程實例的處理，驗證了數值模擬所得結論。

TSP；TSPwin；波場快照；隧道地震預報；斷層

0 引言

TSP(tunnel seismic prediction)系統是瑞士Amberg公司開發的、進行隧道地質超前預報的地球物理勘探設備。研究原理是在隧道內[1]采用小藥量的炸藥爆破作為震源，按照地震勘探中震源與檢波器的互換原理，在隧道的邊壁進行一系列激發，只有兩個檢波器接收，這樣可以實現單炮激發、多道接收的效果。由震源激發產生的地震波遇到地層界面、節理面，特別是斷層破碎帶、溶洞、暗河、巖溶陷落柱、淤泥帶等不良地質界面時，將產生反射波?；诘卣鹂碧皆恚瑱z波器所接收到的反射波傳播速度、延遲時間、波形、強度和方向等均與界面的性質和產狀密切相關[2]，并通過不同振幅、同相軸等波形特征表現出來。TSP主要用于超前預報隧道掌子面前方、上方和下方不良地質體的性質、位置和規模，最大預報距離300～500 m，最高可以分辨直徑為1 m的地質體[3]。

作為工程建設中的一種實用技術，TSP具有追求短平快的特點，前人已從工程應用的角度對TSP進行了研究和分析。劉志剛等[1]對TSP基本原理、技術狀況及其在隧道隧洞超前預報中能夠解決的技術問題、應用效果、能達到的技術指標以及發展進行研究，論述了其存在的問題和解決的途徑。何剛[4]分析了分辨率、頻帶寬度及激發炸藥藥量之間的相互關系，提出了對TSP接收系統的改進意見，并深入剖析系統測量的誤差因素。李堅[5]應用TSP對地層分界面、破碎帶、巖溶等進行預報，預報結果得到了后續施工開挖的驗證。此外:Ashida[6]以隧道爆破和TBM(tunnel boring machine)震動作為震源，對隧道掌子面前方的地質情況做了成像研究，成像結果與前方實際地質情況一致；Bohlen等[7]研究了隧道掌子面上RSSR現象，即隧道地震波場中的Rayleigh面波在掌子面上轉換為S波,遇到反射界面后反射到掌子面上再轉換為Rayleigh面波;王朝令等[8]基于有限元法對隧道地震波場進行了模擬，在其中討論了頻散、邊界條件等環節。

筆者應用數值模擬方法模擬存在斷層的TSP地震波場，提取其時間記錄和波場快照，分析存在斷層的情況下的地震波場特征;再應用TSP后處理軟件TSPwin對得到的時間記錄進行處理，對比模型圖與提取到反射界面的差異，并研究TSP系統的抗噪性能;最后通過工程實例驗證了研究結果,以期為TSP隧道地震預報提供研究依據。

1 TSP系統

TSPwin是TSP系統的數據處理軟件。處理的結果是形成反映相關界面或地質體反射能量的影像點圖和隧道平面、剖面圖，以供工程技術人員解譯。數據處理流程如圖1所示，主要處理步驟包括帶通濾波、初至拾取、Radon變換、Q值濾波、速度分析、深度偏移和反射界面提取。

圖1 TSP數據處理流程圖Fig.1 Flow chart of TSP data processing

①帶通濾波：通過Fourier變換得到地震數據的頻譜，然后根據有效波的頻率范圍設置濾波窗口，將頻率在窗口之外的干擾波濾除。

②初至拾?。禾崛≈边_波到達檢波器的起跳時刻點。

③Radon變換：Radon變換是將t-x域地震記錄轉變為τ-p/q域數據的一種數學變換[9]。τ-p/q變換后各種波的運動學特征和動力學特征會發生變化，使得某些問題在τ-p/q域中解決更為方便[10]。

④Q值濾波[11]：Q為衡量地震波在傳播過程中能量衰減的品質因子，基于它可以分析傳播介質的非彈性情況。

⑤速度分析：速度是地震勘探的關鍵參數，速度分析的主要目的是為水平疊加、偏移等處理步驟提供處理參數。

⑥深度偏移[12]：深度偏移是將檢波器接收到的時間剖面轉換成以空間坐標表示的地質構造剖面。

⑦反射界面提取：提取到反射界面的橢圓形曲線代表界面的反射振幅。其明顯程度表示物理反射界面的波阻抗強度對比，反射振幅越明顯，則反射系數越大，巖石介質的彈性波阻抗分界面屬性差異也就越大。

2 TSP地震波場數值模擬

為了模擬松動圈存在時的TSP地震波場，建立如圖2所示的模型圖。震源采用300 Hz的Ricker子波，這個頻率與TSP實際工程數據采集的震源主頻相近。通常地震子波是由2至3個或多個相位組成的地震脈沖，最常用的Ricker子波由1個波峰和2個波谷組成[10]。在本文的數值模擬研究中，采用0相位的Ricker子波，其表示公式如下：

(1)

式中：u表示位移；f0表示主頻；t代表時刻；t0代表中心時刻。

TSP數據采集時，采用單個三分量接收，24炮激發；根據炮檢互換原理，相當于單炮激發，24道三分量檢波器接收。為了方便計算，在設置模型時，即采用單炮激發多道接收，使得計算模型的觀測系統與TSP實際采集時相同。需要說明的是，在模型數據處理過程中，速度云圖和反射層圖都是從震源位置起算。

模型(圖2)周邊加載了黏彈性邊界，它可以吸收邊界反射波[13]。模型的長度為135 m，高度為75 m，隧道高度6 m，有限元網格劃分為0.5 m×0.5 m，偏移距15 m，道間距1 m，采樣間隔62.5 μs，采樣點數7 218，采樣時窗451.125 ms。由于隧道開挖過程中存在松動圈[14]，因此模擬時在隧道邊壁中設置了斷層帶，其厚度為1 m。工作面前面斷層帶總共設置3層，層厚3 m，層間距15 m，從震源到第一、二、三層異常帶的水平距離分別為65、80、95 m，每個區域內的縱波速度(vP)、橫波速度(vS)以及密度(ρ)如圖2所示。

圖2 模型圖Fig.2 Model diagram

作為地震勘探方法在隧道預報中的一種應用，TSP的分辨率遵循地震勘探分辨率的Knapp準則，TSP能區分的厚度是波長的1/2。本文模型中縱波的波長為4 000/300≈13.3，故TSP能區分6.65 m的厚度。模型中雖然斷層的厚度為3 m，但斷層間距為15 m，故可以分辨。模型中斷層只是作為一個反射層進行勘探，這是因為TSP探測的目的主要是勘探斷層的位置和數量，對斷層厚度可以不予區分。

a.x分量;b.y分量。圖3 數值模擬的時間記錄Fig.3 Time recording of numerical simulation

由于模擬是在二維中進行的，故只能取到x、y兩個分量。原始時間剖面如圖3所示。圖3a中的同相軸依次為直達縱波、直達橫波、左邊界角點繞射橫波和反射橫波。由于反射縱波能量較小,在圖中沒有顯示。時間剖面x分量(圖3a)各同相軸起跳點明顯，y分量(圖3b)各同相軸一致性好且較清晰，信噪比和分辨率都比較高。需要指出的是RSSR現象，即隧道內沿邊壁傳播的Rayleigh面波在工作面角點位置轉換為S波，以體波的形式向前傳播，遇到界面反射回到角點后轉換為Rayleigh面波的現象[7]。圖3a中前方界面的反射橫波經RSSR轉換在隧道內產生Rayleigh面波，由于其與橫波速度相差很小，使得60～80 ms時間段內x分量剖面的反射波同相軸混亂，這影響了反射橫波的分辨率，也降低了地震解釋的可靠性。為了降低這種干擾，提高資料解釋的可靠性，觀測系統應盡量采用大偏移距，使得RSSR現象在隧道內產生Rayleigh面波與前方界面的反射橫波分離開來。

圖4是模型激發y分量的波場快照。t=1 ms時，震源開始激發；t=5 ms時觀察到S波傳播；t=8、10 ms時波傳播到左邊界；t=12 ms時可以看到波被左邊界吸收；t=15、18 ms時繼續在接收排列上傳播；t=20、22 ms時，波傳播到底邊界，底邊界將其吸收；t=25、28、30 ms時，波傳播到工作面角點位置；t=32、35 ms時，RSSR轉換出的S波向前方傳播，以角點為震源產生新的波陣面。此過程是波到達斷層之前的。

圖4 模擬y分量的波場快照Fige.4 The y-component snapshot of simulation

t=38 ms時波到達第一個反射層位置；t=40 ms時S波開始反射回工作面；t=42 ms時波通過第一個反射層，到達第一個反射層與第二個反射層之間的空間位置；t=45 ms時，波傳播到達第二個反射層；t=48、50 ms時，波穿過第二個反射層，此時從第一個反射層返回的波又傳播到排列上；t=52、55 ms時，波傳播到第三個反射層；t=58、60 ms時，穿過第三個反射層，傳播到右邊界上；t=62、65 ms時，從第一個反射層返回的波已經通過接收排列；t=68、70 ms，不斷有從第二、三反射層上的波向排列上傳播；t=72～80 ms，反射波完全穿過排列，此時仍有大量的雜波在快照上存在。

3 模擬數據處理

按照如圖1所示的處理流程應用TSPwin對得到的時間記錄進行處理，各處理步驟均采用TSPwin默認值進行，可以得到如圖5所示速度云圖：橫坐標表示預報距離，計算以震源位置起算；縱坐標表示以隧道向兩邊展布半徑。圖5a中P波速度云圖的低速帶存在位置為75～110 m，模型(圖2)設定的斷層帶位置為65～100 m，兩者非常接近。圖5b中S波速度云圖低速帶在80～120 m，與模型(圖2)設置的斷層帶也比較接近；但所呈現的形態一致性不高，分析是受到RSSR轉換的干擾，使得速度分布云圖出現偏差。綜合兩張圖可知，P波異常帶的速度分布無論是位置還是形態都非常接近模型設定的斷層帶區域，表明采用TSPwin默認值處理所得到的P波速度分布云圖與模型設定非常接近。

圖5 速度云圖Fig.5 Velocity Contours

圖6是提取得到的反射層圖。P波反射界面所處位置與模型設置非常接近，且形態也一致，預報的反射界面正好處于隧道軸線上(圖6a)，這正是隧道超前預報所需要解決的。S波提取反射界面(圖6b)與模型設定位置比較接近，但形態接近程度相對P波較差。對隧道地震超前預報來說，用P波進行預報的準確性比較高，而S波提取到的反射界面準確性相對稍差。

對于如本模型這樣的多層異常帶模型來說，S波沒有P波準確，P波的速度分布和反射界面都與模型設定的位置非常吻合；對隧道地震超前預報來說，P波預報的準確性更高。

圖6 提取反射層Fig.6 Layers of extracting reflection

4 TSP抗噪性分析

為研究TSP的抗干擾能力，對數值模擬所提取的時間記錄加入隨機噪聲，得到如圖7所示的低信噪比信號。從圖7中可以看到：信號湮沒在隨機噪聲中，整個剖面不容易辨識，x分量中除了直達波比較明顯外，反射波幾乎被噪聲所掩蓋；相比于x分量，y分量信噪比稍高。

對加噪聲信號采用TSPwin進行處理，經過帶通濾波、初至拾取、反射波提取、P-S波分離、深度偏移和反射層提取，得到圖8所示的速度云圖和圖9的反射層提取圖。

圖8所示為加噪聲的速度云圖，對比模型(圖2)可以發現，在65 ～100 m區間內，P波的速度云圖存在低速區，這與模型設定相近；相比之下，S波的速度云圖在60～80 m范圍內也存在低速區，但與模型的對應程度較差。

圖9為加噪聲數據所提取的反射層，對比模型(圖2)可以看到：在70～100 m區間內， P波反射界面所處位置與模型設置相近，且形態也較一致，但是相比于圖8未加噪聲的P波提取反射層，對應性有所降低，且反射層強度也有所下降；S波在80～100 m存在反射層，這接近模型的設定位置，與圖8未加噪聲的S波提取反射層相比，一致性和對應性相近。

加噪聲的數據處理表明TSP具有良好的抗噪性能。這是因為TSP在數據采集時會在隧道邊壁打孔，并采用環氧樹脂將套管固定在接收孔中，然后將傳感器與套管緊密耦合，以此接收地震信號，使得原始時間記錄具有很高的信噪比，因此TSP系統所具有的抗噪性能可以滿足實際工程需要。

a.x分量;b.y分量。圖7 加噪聲的時間記錄Fig.7 Time recording with additive noise

圖8 加噪聲的速度云圖Fig.8 Velocity contours with additive noise

圖9 加噪聲的提取反射層Fig.9 Layers of extracting reflection with additive noise

5 實例

考慮到隧道超前預報方法的施工特點，若要驗證預報結果的準確性，只能選取已經開挖完的隧道，比較開挖之前的預報結果和開挖后揭示的地質情況，對比兩者的差異。依據常規TSP觀測排列的布置方法，在隧道邊墻布置地震波信息接收孔，孔深2 m；按約1.5 m的間距布置24個激發孔，孔深1.5 m左右，向下傾斜約10°；每個激發孔裝填藥量為100 g，激發孔和接收孔基本保持在同一個高度。

圖10 實例P波速度云圖(a)和反射層圖(b)Fig.10 Real case of P-wave velocity contours(a) and layers(b) of extracting reflection

圖10是實例P波的速度云圖和反射層圖。圖10a中，110～140 m存在異常低速區域。由于是實例預報，地震波場比較復雜，因此提取的反射層比較多，分布位置比較分散。其中，在110～140 m，存在比較多的反射層，形態以豎向反射層為主(圖10b)。它表示本段圍巖節理裂隙發育，圍巖破碎，推測該段局部存在節理密集帶或破碎帶。開挖后，120～135 m段中間偏右存在黃黏土夾層，在拱腰至拱頂向內延伸約1.5～3.0 m，是一處破碎帶構造，與預報結果相吻合。

6 結論與建議

筆者建立斷層模型進行數值模擬，通過波場快照研究隧道地震波場的傳播特征，對得到的時間剖面按照TSP的常規流程進行處理。得到如下結論：

1)在速度分布云圖中，斷層的速度呈現比周圍速度低的特征，在提取的反射層圖中斷層呈現豎向條狀的特征。

2)對數值模擬所得的時間記錄添加噪聲，結果表明TSP具有良好的抗噪性能。

3)工程實例表明P波預報是比較準確的，為工程中應用TSP進行隧道地震預報的數據處理和資料解釋提供理論依據。

由于隧道施工對TSP時效性的要求較高，采集環境復雜，且地質體多樣，文中的結論是基于垂直多斷層的情況，因此需要研究不同情況下的TSP預報效果，并開展三維隧道地震波場數值模擬下的TSP預報效果研究，這是需要進一步研究的工作。

[1] 劉志剛,劉秀峰.TSP (隧道地震勘探) 在隧道隧洞超前預報中的應用與發展[J].巖石力學與工程學報,2003,22(8):1399-1402. Liu Zhigang，Liu Xiufeng. TSP Application and Deve-lopment in Tunnel Lead Forecast[J]. Chinese Journal of Rock Mechanics and Engineering,2003，22(8):1399-1402.

[2] Alimoradi A, Moradzadeh A, Naderi R, et al.Prediction of Geological Hazardous Zones in Front of a Tunnel Face Using TSP-203 and Artificial Neural Networks[J]. Tunnelling and Underground Space Technology, 2008, 23(6): 711-717.

[3] 楊峰.TSP隧道超前預報中典型地質異常體的射線追蹤數值模擬研究[D].成都：西南交通大學，2007. Yang Feng. The Ray Tracing Numerical Simulation Study of Typical Geological Anomalous Body in Tunnel Advance Prediction[D]. Chengdu: Southwest Jiaotong University, 2007.

[4] 何剛.TSP-203 系統在隧道超前地質預報中的應用研究[D].長沙:中南大學，2005. He Gang. The Research of Application About TSP-203 in Tunnel Geological Predicatio[D]. Changsha:Central South University, 2005.

[5] 李堅.TSP法在鐵路客運專線隧道超前地質預報工作中的應用前景[J].鐵道勘察，2006，31(6):45-49. Li Jian. Application Prospect of the TSP Method in the Lead Geological Predication of the Dedicated-Passenger Railway Tunnels[J]. Railway Survey, 2006, 31(6): 45-49.

[6] Ashida Y. Seismic Imaging Ahead of a Tunnel Face with Three-Component Geophones[J]. International Journal of Rock Mechanics and Mining Sciences, 2001,38(6): 823-831.

[7] Bohlen T, Ulrich L,Wolfgang R, et al. Rayleigh-to-Shear Wave Conversion at the Tunnel Face-From 3D-FD Modeling to Ahead-of-Drill Exploration[J].Geophysics, 2007, 72(6): 67-79.

[8] 王朝令, 劉爭平, 黃云艷,等.隧道地震預報波場的有限元數值模擬[J].吉林大學學報：地球科學版，2014，44(4): 1369-1381. Wang Zhaoling，Liu Zhengping, Huang Yunyan, et al. Wavefield Modeling Based on the Finite Element Method for the Tunnel Seismic Prediction[J]. Journal of Jilin University: Earth Science Editon, 2014,44(4):1369-1381.

[9] 王朝令，劉爭平.τ-p變換在隧道反射地震超前預報波場分離中應用的數值模擬研究[J].地球物理學進展，2012，27(5): 2216-2225. Wang Zhaoling, Liu Zhengping. The Modelling Study ofτ-pTransform Applied in the Tunnel Seismic Prediction[J]. Progress in Geophysics, 2012, 27(5):2216-2225.

[10] 王朝令.隧道地震超前預報中波場分離與反演方法的數值模擬研究[J].成都:西南交通大學，2012. Wang Zhaoling. The Numerical Simulation Research on Separation of Wave Field and Inversion Method in Tunnel Seismic Prediction[D]. Chengdu: Southwest Jiaotong University, 2012.

[11] 徐彥, 蘇有錦,秦嘉政.Q值研究動態[J].地震研究，2004，27(7):385-389. Xu Yan, Su Youjin, Qin Jiazheng. Development ofQ-Value Studies[J]. Journal of Seismological Research, 2004, 27(7): 385-389.

[12] 王有新, 史政軍,李少英.深度偏移的發展與應用概況[J].石油地球物理勘探，1994，29(4): 526-531,534. Wang Youxin, Shi Zhengjun, Li Shaoying. Development of Depth Migration and Its Application[J]. Oil Geophysical Prospecting, 1994.29(4):526-531,534.

[13] 王朝令,劉爭平.黏彈性邊界條件在ANSYS有限元波場模擬中的實現[J].大地測量與地球動力學，2012，32(2): 28-31. Wang Zhaoling, Liu Zhengping. Realization of Vis-coelastic Boundary Condition in Wave Field Simulation with Ansys Finite Element Software[J]. Journal of Geodesy and Geodynamics, 2012, 32(2): 28-31.

[14] 宋宏偉, 王闖,賈穎絢.用地質雷達測試圍巖松動圈的原理與實踐[J].中國礦業大學學報，2002，31(4):43-46. Song Hongwei, Wang Chuang, Jia Yingxuan. Priciple of Measuring Broken Rock Zone Around Underground Roadway with GPR and Its Application[J]. Journal of China University of Mining& Technology, 2002, 31(4): 43-46.

Numerical Simulation of TSP Fault Model

Lin Yi1,2, Liu Zhengping1, Wang Zhaoling3，Xiao Di3

1.FacultyofGeosciencesandEnvironmentalEngineering,SouthwestJiaotongUniversity,Chengdu610031,China2.NorthwestResearchInstituteCo,.LtdofC.R.E.C,Lanzhou730000,China3.CollegeofCivilEngineering,SichuanAgriculturalUniversity,Chengdu611830,China

During tunnel excavation, a variety of geological disasters might be encountered, such as faults, caves, et.al. Tunnel seismic prediction (TSP) is adopted to mitigate the possible damages. Although TSP technology is used widely, the research about TSP is currently focused on its engineering application cases. We use the finite element method to simulate the tunnel seismic wave field, employ wave field snapshots and time recording method on the impact of faults on the characteristics of the propagation of tunnel seismic wave field, and inversely process the time record of model containing the fault. The digital model of the velocity scattered image and the reflection interface position are obtained, and the fault position from velocity scattered image processed with the default values set by using TSPwin is agreed to the one from the model. In respect to the layered model for an abnormal velocity zone, P-wave is more precise. The system of TSP is strong for its feature of anti-noise. The numerical simulation is verified finally through the process and analysis of the engineering cases.

TSP；TSPwin；snapshot；tunnel seismic prediction；fault

10.13278/j.cnki.jjuese.201506305.

2015-02-28

國家自然科學基金項目(40874051，51508358)

林義(1980--)，男，博士研究生，高級工程師，主要從事隧道地震超前預報研究，E-mail:766708597@qq.com

王朝令(1980--)，男，副教授，博士，主要從事隧道地震波場模擬研究，E-mail:wong8010@gmail.com。

10.13278/j.cnki.jjuese.201506305

P631.4

A

林義，劉爭平，王朝令,等. TSP斷層模型數值模擬.吉林大學學報:地球科學版,2015,45(6):1870-1878.

Lin Yi, Liu Zhengping, Wang Zhaoling, et al. Numerical Simulation of TSP Fault Model.Journal of Jilin University:Earth Science Edition,2015,45(6):1870-1878.doi:10.13278/j.cnki.jjuese.201506305.