基于時間序列的中壩隧道涌水分析

張 吉力

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

基于時間序列的中壩隧道涌水分析

張 吉力

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

隧道涌水是四川巖溶地區隧道開挖過程的重大影響因素之一，決定隧道的順利施工與否.以中壩隧道為例，通過工程地質勘查和滲流管道調查，結合降雨量與涌水量的時間序列分析，更加準確地判斷富水溶腔的連通性，闡述該類隧道涌水形成機理，有效地為預測和解決隧道涌水問題提供借鑒.

隧道涌水;工程地質;巖溶;時間序列分析

0 引言

巖溶地區公路、鐵路隧道的施工過程中，涌水、突水等問題對隧道施工安全和進度有著重大影響，往往由于對涌水補給模式和富水溶腔連通性的判斷錯誤，或處置措施的不當，導致重大的人員傷亡和財產損失.此外，隧道長期疏排水無疑會導致隧道地區的水環境發生改變，甚至對居民飲水造成影響.所以，在巖溶地區隧道施工過程中要重點針對涌水、突水等工程地質問題做詳細的分析和判斷.而在隧道的掘進過程中遇到高壓涌水的不良地質災害時，正確的分析和判斷隧道涌水、突水現象形成的原因和機理，可以給災害處置方案提供有力的參考，幫助解決和處理涌水、突水問題.時間序列分析方法最早由Jenkins等于1968年提出，之后被廣泛用于水文研究領域，尤其是用于研究巖溶地區的水文地質過程［1－3］，本研究嘗試通過一定時間監測隧道涌水量和降雨量，分析涌水對降雨的響應特征，定性定量地判定涌水的水源補給形式，以達到了解涌水水源補給情況的目的.

1 研究區概況及工程地質條件

1.1 研究區概況

2012年7月14日上午11時，中壩隧道(見圖1)進口施工段里程D9K55+221掌子面發生局部高壓涌水現象.涌水點出現在掌子面左側拱腰位置的三處垂直分布的炮眼附近，水柱沿掌子面方向噴出，首日噴距最大達到8 m，且呈現出逐漸增大的趨勢.涌水由清澈漸變為渾濁，同時攜帶大量泥沙.通過現場勘察，可以確定隧道涌水上方有暗河通過，但是僅僅通過掌子面涌水情況很難判斷該涌水的水源補給機制和富水溶腔的連通性.

圖1 隧道區平面示意圖

1.2 工程地質條件

研究區內構造條件較為簡單，區域內所涉及的地層單元主要為第四系坡洪積()層、坡殘積)層，三疊系雷口坡組(T2l)、嘉陵江組(T1j)、飛仙關組(T1f)，二疊系龍潭、長興組(P2l+c)，茅口棲霞組(P1m+q).為一單斜構造，未見斷裂構造體通過，局部發育小型褶皺.層面在進口處多變，其他大部分以22°緩傾角為主，圍巖以Ⅲ、Ⅳ級圍巖為主，節理裂隙成組發育.

2 涌水機制分析

受到飛仙關組(T1f)非可溶巖地層的隔擋，以P1地層為主的水文地質單元對隧道無影響.因此，隧址區內隧道所穿越的三疊系雷口坡、嘉陵江組灰巖含水層是隧道開挖過程中的主要涌水源.

該區域內巖溶發育的特點是地表水和地下暗河管道相連通，巖溶以水平溶蝕與垂直溶蝕相互交替呈現，形成規模較大的地下巖溶通道，尤其是嘉陵江組巖溶十分發育，主要表現為地表大量的溶蝕洼地、落水洞和巖溶槽谷，并且順層發育有地下暗河管道系統［4］，如圖 2 所示.

圖2 隧道區匯水情況平面示意圖

隧道D9K54+100～D9K56+230段切割穿越三疊系下統嘉陵江組(T1j)灰巖，可溶巖段長2 130 m，隧道在D9K55+221掌子面處發生涌突水災害.水平距離隧道軸線該300 m處地表發育有一個大型的溝槽——紙槽溝，槽谷長900 m，呈 N38°W 向展布，寬220 m，中間洼地底部最低高程為999 m，具規模大，匯水條件好的特點.同時，隧道頂部地表巖溶十分發育，多見溶洞、溶蝕槽谷、洼地、落水洞，具備很好的地下水賦存空間和運移通道，有利于地表水的徑流及排泄.隧址區所調查巖溶負地形大多發育在800～1 000 m高程左右，現場調查顯示，隧址區洼地、落水洞、巖溶槽谷等巖溶形態星羅棋布，均受區內的構造線控制，順層沿軸向呈串珠狀展布，尤其位于隧道東側的溶蝕現象，以紙槽溝為起始點順巖層走向延伸地表發育的各個巖溶形態的底部高程具有階梯狀的下降趨勢，并且呈現一定的方向指示性，如圖3所示.

圖3 洼地、落水洞高程分布圖

研究區暗河管道沿線地表發育的落水洞、洼地和溶蝕槽谷均是很好的匯水區域，為地下暗河的補給提供了良好的導水通道.經詳細勘察，暗河位置如圖4所示.

圖4 地下暗河管道平面示意圖

圖5 地下暗河管道與隧道縱剖面示意圖

圖6 地下暗河管道與隧道橫斷面示意圖

同時，地表溶蝕的發育高程均位于隧道頂板之上，為此暗河管道也在隧道上部，隧道位于地下水水位以下的水平徑流帶中，如圖5、6所示地下暗河管道的發育對于隧道的開挖存在較大的安全隱患，如2012年7月14日，在隧道開挖至里程D9K55+221處曾發生突然性的大型涌突水災害.

3 降雨和涌水量監測分析

中壩隧道發生高壓涌水情況之后，本研究通過地方氣象局的相關降雨量信息和現場涌水量實際監測情況做了實際的統計和比較，詳見表1、圖7.

表1 中壩隧道降雨量和涌水量數據比較

圖7 涌水量、涌水量監測跟蹤圖

根據以上數據信息，可以發現隧道高壓涌水的最大涌水量出現在7月30日，實際降雨的最大降雨量值出現在7月23日，由此可以發現該涌水量與地表降雨量有明顯的滯后現象.在此基礎上，本研究采用基于時間序列及功率譜的方法分析了隧道涌水對于地表降水的響應關系，結果如圖8、9所示.

由圖8可知，降雨時間序列的自相關函數在k=0～5附近下降很快，并且當k≥5時的所有r(k)均小于0.2.此表明，降雨時間序列無記憶效應，為純隨機序列，這與圖9的自功率譜分析結果(無明顯峰值)相一致.然而，隧道涌水時間序列卻表現出完全不同的變化特征，首先是其自相關函數出現了明顯的拖尾，在0～20 d滯后期內R(k)均超過0.2，表明在經歷一次輸入脈沖后系統的記憶長度可達20 d;其次，從自功率譜看，S(f)在整個頻域范圍內未出現明顯的峰值，由此表明隧道涌水時間序列不存在小于40 d的短周期，但不排除有大于本次觀測時段長度長周期的可能，這從另一個側面反映出隧址區巖溶地下水系統的調節能力較強，致使隧道涌水與降雨時間序列的結構特征出現了明顯的差異.

圖8 涌水量、降雨量相關函數圖

圖9 涌水量、降雨量自功率譜密度函數圖

此外，降雨(P)與隧道涌水(TQ)時間序列的互相關功率譜函數、互功率相位譜、譜相關函數如圖10～12所示.

圖10 涌水量、降雨量互功率譜密度函數圖

圖11 涌水量、降雨量功率譜相位函數圖

圖12 涌水量、降雨量譜相干函數圖

由圖10可知，在整個時域范圍內，降雨和隧道涌水的互相關函數有只有一處比較明顯的峰值，出現在k為9時.但總體而言，二者的相關性并不很強，超過95%置信水平的個數僅約1/5.相位譜用于估計一個時間序列對另一時間序列的主導作用，如超前或滯后.由于四象限反正切變換得到的相位在［－π，π］之間存在相位卷繞現象，故需展開得到實際的相位變化［5］.由φ12(f)=2πfd可知，φ12(f)為f的線性函數，斜率為K=2πd.據此，通過擬合圖11便可求得隧道涌水對降雨響應的平均滯后時間為，

中壩隧道涌水處暗河高度約為200 m，可大致估算出隧址區含水層的平均垂直滲透系數為16.82 m/d.相干譜可以理解為2個時間序列在不同頻率時的關聯程度，用于表征輸入和輸出時間序列在某個頻率處是否存在同類型的振蕩，因此可用來度量輸入和輸出之間的線性相干程度.若相干譜值接近1，則說明輸入信號的某個變化在輸出中也得到了成比例的反映;相反，若相干譜值遠小于1，則說明系統中存在很強的非線性特征.由圖12可知，隧道涌水和降雨的相干系數介于0～0.5，平均0.176，表明二者之間的相干程度降低.這與互相關分析和相位譜反映出的信息一致.

結合涌水量和降雨量的時間序列分析，隧道涌水量與降雨情況有著明顯的時域上的差異，涌水量作為降雨量的響應序列，涌水量和降雨量的關系并不密切，主要受滲流通道的影響，涌水的水源補給主要是含水層滲流裂隙補給，而非與暗河通道聯通直接補給.

4 結論

巖溶地區的隧道建設過程中，經常遇到隧道涌水情況，中壩隧道涌水，通過對隧道施工區區域地質情況和涌水點區域工程地質勘查分析、對涌水量和降雨量的量化分析可以準確的判斷隧道涌水的性質和形式，從而判定隧道涌水的補給形式的滲流通道，然后決定隧道施工的處理辦法.

中壩隧道涌水的發生機制屬于復合型涌水，由于隧道上方約200 m處有暗河經過，通過層面和裂隙滲流在隧道前方形成富水溶腔，而非地表水直接補給，隧道開挖至此處時，發生高壓涌水.經過一段時間的涌水泄壓，隨著富水溶腔中的水涌出，涌水量、水壓也隨之減小.

根據后期的超前鉆孔揭露的情況看，與判斷的情況相吻合，時間序列分析與地質勘察相結合的方法對于探究巖溶地區的涌水成因分析有較好的效果.

［1］劉建.巖溶隧道涌水對降雨的響應特征［J］.煤田地質與勘探，2010，38(2):32 －35.

［2］徐則民，黃潤秋，羅杏春.特長巖溶隧道涌水預測的系統辨識方法［J］.水文地質工程地質，2002，29(4):50 －54.

［3］李蒼松，王石春，陳成宗.巖溶地區鐵路工程地質研究［J］.鐵道工程學報，2005，22(S1):357 －363.

［4］劉建.銅鑼山隧道涌水對降水的響應特征研究［J］.工程勘察，2010，38(2):37 －41.

［5］劉招偉，張民慶.某鐵路隧道巖溶突水的連通性特征［J］.隧道建設，2006，26(1):6 －9.

Analysis of Zhongba Tunnel Gushing Based on Time Series

ZHANG Jie
(The State Key Laboratory of Geohazard Prevention and Geoenvironment Protection，Chengdu University of Technology，Chengdu 610059，China)

The tunnel gushing is one of the major influencing factors of tunnel excavation in karst areas of Sichuan，which decides whether the construction of the tunnel is successful or not.This paper takes Zhongba tunnel from Xuyong to Dacun railway as an example.By conducting engineering geology exploration and seepage pipes survey，and combining time series analysis of the rainfall amount and water inflow amount，we can judge more accurately the connectivity of rich aqueous chambers.This paper discusses such kind of tunnel gushing formation mechanism，which effectively provides a reference for predicting and solving the problem of tunnel gushing.

tunnel gushing;engineering geology;karst;time series analysis

P641

A

1004－5422(2014)01－0096－05

2013－11－05.

張 吉力(1987—)，男，碩士研究生，從事巖土工程技術研究.