賈金曉, 李天斌, 孟陸波, 劉陽飛, 曹海洋

(地質災害防治與地質壞境保護國家重點實驗室(成都理工大學)，成都 610059)

?

千枚巖隧道破碎帶TSP超前預報解譯標志

賈金曉, 李天斌, 孟陸波, 劉陽飛, 曹海洋

(地質災害防治與地質壞境保護國家重點實驗室(成都理工大學)，成都 610059)

千枚巖等軟巖隧道施工中易出現垮塌等地質災害。以典型千枚巖隧道——汶馬高速公路鷓鴣山隧道為例，以地震波傳播理論為依據，結合工程地質分析及地震波正演模擬，探討TSP預報中針對千枚巖構造破碎帶在不同含水率狀況下的地震波反射特性及解譯標志，結果表明：縱波及橫波遇千枚巖構造破碎帶反射現象明顯，橫波較縱波反射更強烈；破碎帶富水時，千枚巖泥化變軟，縱波反射加強；從深度偏移角度分析，破碎帶內正負反射頻繁變化，單個反射層延伸性差。研究結果可以提高TSP在軟巖隧道中的預報準確性。

千枚巖；TSP；破碎帶；正演模擬；解譯標志

以千枚巖地層為代表的典型軟巖，因其強度低、易變形、遇水易軟化的工程特性，在隧道施工過程中易出現大變形、垮塌等災害。軟巖隧道中的構造破碎帶、泥化帶等軟弱地質體更是引發(fā)隧道垮塌、大變形的危險地帶。目前，超前地質預報作為隧道信息化動態(tài)設計和施工的重要組成部分，及時預測掌子面前方不良地質情況，采取有效預防措施，確保隧道安全施工，已在工程界得到廣泛認可。

TSP超前預報技術被認為是國內外最先進的隧道長距離預報方法之一。不少專家學者針對隧道不良地質體的TSP超前地質預報進行了研究，取得了一定的研究成果。其中，李術才等[1]、李利平等[2]、周宗青等[3]研究了特長深埋隧道基巖裂隙水的綜合預測方法。許振浩等[4]分析了巖溶隧道斷層破碎帶、地下水充填型、泥加石充填型、軟弱夾泥充填型溶洞及富水巖層的地震波反射特性。王延壽等[5]運用TSP超前預報法對秦安高風險隧道斷層(碎裂帶)、巖性接觸帶和富水區(qū)等不良地質現象進行預報研究。張劍等[6]將超前預報技術用于軟巖隧道進行綜合判別圍巖狀況。自TSP超前地質預報方法在國內推廣以來，在很多工程項目中有所運用。在理論方面，許多專家學者也做了研究，而多數專家學者關于TSP地震波反射特性及解譯標志的研究多基于硬質巖體。而對于軟巖隧道，TSP超前地質預報的解譯標志研究甚少。

鷓鴣山隧道位于四川省阿壩州理縣境內，起訖地面海拔高度為3 050～4 623.8 m，全長約8.8 km，屬典型高原特長隧道。本文研究區(qū)域主要為隧道進口段千枚巖地層，該地層受米亞羅斷裂影響，局部巖體軟弱破碎，多呈薄層狀構造，節(jié)理、裂隙極為發(fā)育，整體性差，強度低，遇水泥化變軟，在開挖過程中，多次發(fā)生掌子面坍塌，嚴重影響施工安全。本文在運用TSP法對鷓鴣山隧道進行超前預報過程中，結合地震波正演模擬，分析千枚巖破碎帶在干燥和富水情況下TSP方法的解譯標志，為后期隧道超前地質預報及類似工程提供借鑒。

1 TSP超前地質預報解譯準則

TSP系統(tǒng)的數據處理由其自帶的TSPwin數據處理軟件完成，輸入相關處理參數可得到掌子面前方預報段范圍內的深度偏移圖，各里程段圍巖的縱波速度(vP)、橫波速度(vS)、泊松比(μ)、密度(ρ)、彈性模量(E)等參數，以及波阻抗反射界面位置。通過相應的分析預報掌子面前方的圍巖地質情況。

根據 Dickmann及國內學者對TSP解譯方法研究，總結了TSP解譯準則[7-14]：(1)深度偏移圖中，弧線表示波阻抗反射界面，反射能量的強弱根據偏移大小確定，正反射振幅即紅色區(qū)域，表示巖質較硬，負反射振幅即藍色區(qū)域，表示巖質較軟；(2)縱波速度與巖體硬度呈一定正相關性，若vP下降，則表明裂隙密度或孔隙度增加；(3)若橫波反射較強，則表明巖層含水或流體；(4)vP/vS或泊松比有較大的增加，常常因流體的存在而引起。

2 典型例子

例一隧道K181+265～K181+376里程段TSP預報成果參數曲線(圖1)，根據巖體物性參數曲線和P波深度偏移圖分析，A區(qū)巖體物性參數曲線(圖1)中，縱橫波速度均有明顯降低，橫波速度降低10.91%，縱波速度降低2.3%，縱橫波速比升高，同時動態(tài)彈性模量明顯降低；同時P波深度偏移(圖2)中，在區(qū)域前端界面，以強烈的負反射開始，在區(qū)域后端界面，以強烈的正反射結束，區(qū)域內巖質較均一，正負反射層少，表現為前端界面負反射的延續(xù)，單個反射條帶寬。預測該段圍巖為軟弱夾層，巖質較均一，力學強度低，地下水較發(fā)育。

開挖揭露情況顯示A區(qū)域(K181+306～K181+334)為千枚巖富水破碎帶，鑲嵌結構，巖體濕潤，局部線狀出水，千枚巖崩解為細顆粒，局部夾塊狀千枚巖，泥化染手現象嚴重，有明顯流體特性，在K181+322掌子面發(fā)生坍塌(圖3)。預報結果與開挖情況對比見表1。

圖1 預報結果2D視圖Fig.1 2D view of prediction results位置：K181+265～K181+376

圖2 P波深度偏移圖Fig.2 Diagram showing P-wave depth shift位置：K181+265～K181+376

例二

ZK181+038～ZK181+150里程段TSP預報成果參數曲線(圖4)，從巖體物性參數曲線分析，B區(qū)域(ZK181+038～ZK181+060)縱橫波速均起伏變化，橫波速度降低10.95%，縱波速度降幅1.8%，縱橫波速比升高，動態(tài)彈性模量降低；從P波深度偏移角度分析(圖5)，在區(qū)域前端界面，以強烈的負反射開始，在區(qū)域后端界面，以強烈的正反射結束，局部段正負反射層頻繁變化。預測該段圍巖軟硬巖交替出現，預報段內除中間局部巖體外，前后兩段巖體極破碎，力學強度較低，地下水較發(fā)育，可表現為滲滴水-線狀出水。

圖3 K181+322掌子面坍塌Fig.3 Photograph showing tunnel face collapse in K181+322

圖4 預報結果2D視圖Fig.4 2D view prediction results位置：ZK181+038～ZK181+150

開挖揭露情況顯示，B區(qū)域為千枚巖破碎帶，巖體松散破碎，呈薄片狀，地下水弱發(fā)育，局部滲滴水，在ZK181+054掌子面發(fā)生坍塌，坍塌段長4 m，坍塌現場情況見圖6。預報結果與開挖情況對比見表1。

圖5 P波深度偏移圖Fig.5 Diagram showing P-wave depth shift位置：ZK181+038～ZK181+150

Table 1 Comparison of prediction results with excavation situation for K181+306～K181+334 and ZK181+038～ZK181+060

位置預報結果開挖情況(A)K181+306～K181+334圍巖強度低，局部含軟弱夾層，注意防止局部坍塌。滲滴水?線狀出水K181+322～K181+332段，千枚巖破碎帶，鑲嵌結構，局部為千枚巖細顆粒，無層理，巖體泥化染手現象嚴重，地下水局部線狀滲出(B)ZK181+038～ZK181+060預報段內除中間局部巖體外，前后兩段巖體極破碎，力學強度較低，地下水發(fā)育，可表現為滲滴水?線狀出水ZK181+042～ZK181+046和ZK181+052～ZK181+058里程段，巖體松散破碎，受構造應力影響，巖體呈薄片狀，地下水弱發(fā)育，滲滴狀出水

圖6 ZK181+054掌子面坍塌Fig.6 Photograph showing tunnel face collapse in ZK181+054

對比分析例一和例二，實際開挖情況顯示：2個特殊段揭露的巖體均為千枚巖破碎帶。不同的是，例一中千枚巖破碎帶地下水較為發(fā)育，表現為線狀出水，局部巖體泥化變軟，圍巖松散呈顆粒狀；例二中千枚巖破碎帶地下水發(fā)育較弱，千枚巖呈薄片狀。

TSP預測結果顯示：地震波遇2種千枚巖破碎帶時，橫波降幅較為一致，約為10.9%；例一中，破碎帶富水時，縱波降幅更大。

3 不良地質體地震波響應正演模擬

通過地震正演模擬技術，查明地震波在復雜地質體中的傳播特征，是研究地震波傳播規(guī)律的有效途徑。地震波場正演方法主要有2種：幾何射線法和波動方程法。幾何射線法可以精確地映射出地震波的射線路徑和旅行時間等運動學特征，但對于一些復雜的巖性信息和地質構造容易產生盲區(qū)；而波動方程法以地震波動力學性質為依據，能夠更好地模擬復雜地層的地震波場特征。

Tesseral-2D是一個基于波動方程的地震波場數值模擬軟件包。通過Tesseral-2D高階差分網格，來模擬地震波在介質中的傳播，可以考慮速度和密度的任意變化，模擬復雜地質剖面的數字模型，快速而精確地計算在介質中傳播的P波和SV波，幫助研究者較好地了解地下介質中的地震波動過程。

結合上述典型例子，采用Tesseral-2D正演模擬方法，分析千枚巖破碎帶在地下水弱發(fā)育和富水2種情況時的地震波傳播規(guī)律，可以為隧道超前預報提供解釋依據，對TSP解譯成果進行驗證，提高解釋精度。

3.1 模型設置

本文采用Tesseral-2D軟件進行正演分析，設定2種地質預報模型如圖7所示，分別為千枚巖破碎帶(干燥)及富水破碎帶(含水率10%)。考慮到地震波反射干擾，設置模型大小為400 m×600 m，破碎帶寬度均為10 m，設定20炮逐點放炮激發(fā)地震波，炮孔間隔1.5 m，最前炮孔距離破碎帶20 m；在距離最后一個炮孔15 m位置，每間隔1 m設置1道檢波器，共4道接收并記錄地震波；震源峰值頻率設為1 000 Hz，子波類型為對稱單峰(Symmetrical)，地質體物性參數見表2。

3.2 正演模擬結果分析

千枚巖干燥破碎帶及富水破碎帶(t=20 ms)波場如圖8所示，模擬結果顯示地震波在破碎帶前后界面均發(fā)生強烈的反射現象；觀測不同時刻波場快照可發(fā)現：縱波速度較橫波速度快，縱波提前于橫波到達；根據反射振幅分析，橫波反射振幅更大，反射能量更強。對比2種地質模型可以看出，破碎帶富水時的縱波反射能量強于地下水不發(fā)育的破碎帶。

圖7 兩種地質模型Fig.7 Two geological models

Table 2 Physical and mechanical parameters

vP/m·s-1vS/m·s-1ρ/g·cm-3千枚巖390020002．67水140001．00

圖8 波場快照(t=20 ms)Fig.8 The snapshot of the wave field

通過千枚巖破碎帶的TSP測試成果及正演模擬結論，千枚巖因其特有的質軟、遇水泥化崩解的特性，使其對橫波的反射效果增強；而硬質巖體破碎帶內“縱波遇斷層破碎帶反射較強，若巖層富水橫波也反射較強”[4]，破碎帶對縱波的響應更加明顯。

4 結 論

綜合分析TSP現場測試成果及Tesseral-2D軟件對千枚巖破碎帶2種含水條件下的地質模型正演模擬結果，可以得出不同含水條件下千枚巖破碎帶的TSP解譯準則。共同點：(1)縱波和橫波均對千枚巖破碎帶有較明顯反射現象，橫波降幅約為10.9%，縱波降幅約為2%；(2)橫波反射較縱波強；(3)vP/vS增大，動態(tài)彈性模量降低。不同點：當千枚巖破碎帶富水時，縱波反射現象增強。

千枚巖破碎帶與硬質巖體破碎帶相比，千枚巖對橫波的反射效果增強。

[1] 李術才，李樹忱，張慶松，等．巖溶裂隙水與不良地質情況超前預報研究[J].巖石力學與工程學報，2007, 26(2):217-225. Li S C, Li S C, Zhang Q S,etal. Forecast of karst-fractured groundwater and defective geological conditions[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26 (2): 217-225. (in Chinese)

[2] 李利平，石少帥，李術才，等．特長深埋隧道裂隙水綜合預測方法與應用[J].地下空間與工程學報，2013, 9(3): 603-609. Li L P, Shi S S, Li S C,etal. Comprehensive prediction of fissure water in extra-long deep tunnel and its application[J]. Chinese Journal of Underground Space and Engineering, 2008, 9(3): 603-609. (in Chinese)

[3] 周宗青，李術才，李利平，等．特長深埋隧道基巖裂隙水探測與應用研究[J].地下空間與工程學報，2012, 8(1):99-104. Zhou Z Q, Li S C, Li L P,etal. Detection of bedrock fissure water and its application in extra long and deep tunnel engineering[J]. Chinese Journal of Underground Space and Engineering, 2008, 8(1): 99-104. (in Chinese)

[4] 許振浩，李術才，張慶松，等．TSP超前地質預報地震波反射特性研究[J].地下空間與工程學報，2008,4(4):640-644. Xu Z H, Li S C, Zhang Q S,etal. Reflection characteristic of seismic wave in TSP advance geological prediction[J]. Chinese Journal of Underground Space and Engineering, 2008, 4(4): 640- 644. (in Chinese)

[5] 王延壽，謝婉麗，葛瑞華，等．TSP超前預報法在秦安高風險隧道中的應用[J].科學技術與工程，2014, 14(34):97-103. Wang Y S, Xie W L, Ge R H,etal. Application of TSP advanced prediction in Qinan high-risk tunnel[J]. Science Technology and Engineering, 2014, 14(34): 97-103. (in Chinese)

[6] 張劍．軟巖隧道采用超前地質預報及巖性分析綜合判析圍巖狀況[J].公路，2014(9):193-196. Zhang J. Identification of the wall rock statues in soft rock tunnel by means of geologic forecast and rock properties analysis[J]. Highway, 2014(9): 193-196. (in Chinese)

[7] 戴前偉，何剛，馮德山．TSP-20在隧道超前預報中的應用[J].地球物理學進展，2005, 20(2):460-464. Dai Q W, He G, Feng D S. Application of the TSP-20 System in geological advanced prediction of tunnel[J], Progress in Geophysics, 2005, 20(2): 460-464. (in Chinese)

[8] 張前進．基于TSP多波關系的圍巖穩(wěn)定性等級判定方法研究[D].北京:中國地質大學檔案館，2013. Zhang Q J. Research on the Method of Crade Determination for Surrounding Rock Stability Based on TSP Multiwave Relation [D]. Beijing: The Archive of China University of Geosciences, 2013. (in Chinese)

[9] 周慶國．TSP-203系統(tǒng)與ProEx地質雷達在超前地質預報中的應用[D]. 北京:中國地質大學檔案館，2010. Zhou Q G. Application of TSP-203 and ProEx Ground Penetrating Radar in Geological Prediction [D]. Beijing: The Archive of China University of Geosciences, 2010. (in Chinese)

[10] 何剛．TSP-203系統(tǒng)在隧道超前地質預報中的應用研究[D].長沙:中南大學檔案館，2005. He G. Study of TSP-203 System for Advanced Geological Prediction in Tunnels [D]. Changsha: The Archive of Central South University, 2005. (in Chinese)

[11] 仵軍勝．地震波超前預報系統(tǒng)解譯技術應用研究[J].國防交通工程與技術，2005(3):65-69. Wu J S．The study of the interpretation technology of TSP geologic forecast[J]. Traffic Engineering and Defense Technology, 2005(3): 65- 69. (in Chinese)

[12] 齊傳生．TSP-202 隧道地震波超前地質預報系統(tǒng)的應用[J]. 世界隧道，1999(1):36-40. Qi C S．Application of TSP-202 tunnel seismic wave prediction system on forecast[J]. World Tunnel, 1999(1): 36-40. (in Chinese)

[13] 雷啟云，諶文武，張景科．巖體波速在TSP-203解譯中的應用[J].西部探礦工程，2005(3):122-124. Lei Q Y, Chen W W, Zhang J K. The application of rock mass velocity in TSP-203’s interpreting[J]. West-China Exploration Engineering, 2005(3):122-124. (in Chinese)

[14] 趙宇．GPR及TSP在隧道超前地質預報中的解譯標志研究[D].成都:成都理工大學檔案館，2009. Zhao Y. Study on the Interpretation of GPR and TSP in Tunnel Geological Forecast[D]. Chengdu: The Archive of Chengdu University of Technology, 2009. (in Chinese)

Study of seismic wave reflection characteristics of unfavorable geological bodies in phyllite tunnel

JIA Jin-xiao, LI Tian-bin, MENG Lu-bo, LIU Yang-fei, CAO Hai-yang

StateKeyLaboratoryofGeohazardPreventionandGeoenvironmentProtection,ChengduUniversityofTechnology,Chengdu610059,China

Collapse and other geological disasters often occur in the process of tunnel construction in soft rock like phyllite. Accurate prediction of geologic bodies in front of tunnel face is particularly important. Taking the typical phyllite tunnel, the Zegu mountain tunnel in Wenchuan-Maerkang Highway of Sichuan as example, on the basis of seismic wave propagation theory, combined with engineering geological analysis and forward modeling, the reflectivity and interpretation signs of seismic waves under different water conditions for phyllite fracture zone in the TSP prediction are discussed. It shows that reflection phenomenon is obvious when longitudinal wave and transverse wave meet fracture zone and transverse wave reflection phenomenon is more intense than the longitudinal wave. When fracture zone is rich in water, phyllite becomes softer and transverse wave reflection intensified. Analysis of the perspective of depth migration reveals that normal reflection and negative reflection shift frequently and the extensibility of single reflecting layer are poor. The result of this research improves the accuracy of TSP prediction in the construction of soft rock tunnel.

phyllite; TSP; structural crushing zone; forward modelling; interpretation signs

10.3969/j.issn.1671-9727.2016.06.13

1671-9727(2016)06-0751-06

2015-11-04。 [基金項目] 國家自然科學基金項目(41102189); 地質災害防治與地質環(huán)境保護國家重點實驗室團隊重點項目(SKLGP2009Z002); 成都理工大學優(yōu)秀創(chuàng)新團隊培育計劃項目(HY0084) 。

賈金曉(1990-),男,碩士研究生,研究方向：地質災害研究與防治, E-mail:754206159@qq.com。

U456.33; P642

A