郭佳奇，錢 源，黃 猛

(河南理工大學(xué) 土木工程學(xué)院，河南 焦作 454003)

因周邊隱伏溶洞的隱蔽性，使隧道施工中易于突然遭受到突水涌泥等巖溶災(zāi)害，危害嚴(yán)重。為此眾多學(xué)者開展了大量研究工作。任美鍔等[1]研究發(fā)現(xiàn)巖溶引發(fā)的隧道工程災(zāi)害主要表現(xiàn)為隧道突然涌水、隧道頂部溶洞充填物垮塌等問題。汪從錦[2]結(jié)合魯布革電站引水隧洞討論了不同洞徑圓形隱伏溶洞對隧道圍巖位移的影響；吳夢軍等[3-4]通過現(xiàn)場試驗(yàn)和有限元方法分析了溶洞規(guī)模和位置對公路隧道圍巖穩(wěn)定性的影響。李科[5]采用數(shù)值模擬和模型試驗(yàn)方法研究了隧道拱部不同大小、距隧道不同距離的溶洞對隧道圍巖變形的影響。劉招偉等[6]采用有限差分軟件分析了溶洞位于隧道上側(cè)方和下側(cè)方時(shí)，防突巖層安全厚度的取值范圍。莫陽春等[7]研究了隧道與其底部或側(cè)部的隱伏溶腔間圍巖向2個相反方向變形的特征。張鋒等[8]運(yùn)用有限元方法分析了隧道側(cè)部溶洞對圍巖穩(wěn)定性的影響。謝舉[9]利用數(shù)值方法研究了不同影響因素下巖溶隧道圍巖的受力及穩(wěn)定性。郭明[10]采用COMSOL Multiphysics軟件研究了在不同類型和分布形態(tài)的隱伏溶洞影響下隧道圍巖滲流場、應(yīng)力場和位移場的變化規(guī)律。蘇海健[11]采用有限差分軟件研究了承壓水壓力、保護(hù)層厚度等對巖溶隧道圍巖位移場、塑性區(qū)及主應(yīng)力的影響規(guī)律。朱國偉[12]研究了長昆客運(yùn)專線隧道周邊非充填型溶洞對隧道施工過程中圍巖應(yīng)力、位移及塑性區(qū)的影響。

上述文獻(xiàn)的研究對象為隧道周邊的隱伏干溶洞，且多采用有限元方法，較難準(zhǔn)確刻畫隧道與高壓富水隱伏溶腔間防突巖層突水過程。因此，本文首先從微觀上分析巖溶隧道的突水機(jī)理，然后選取工程中較為常見且對隧道危害最大的隧道邊墻外賦存高壓富水溶腔這一類型的新建山嶺隧道作為研究對象，采用離散元方法分析巖溶隧道防突巖層漸進(jìn)破壞過程中滲流場、位移場等的變化規(guī)律，研究成果對提高隧道巖溶災(zāi)害的可預(yù)見性具有意義。

1 巖溶隧道突水機(jī)理及類型

1.1 巖溶隧道突水機(jī)理

1)巖溶水的軟化和溶蝕作用

隧道與周邊富水隱伏溶腔間防突巖層往往處于飽水狀態(tài)，會受到巖溶水軟化和溶蝕的雙重劣化作用。巖溶水軟化作用降低了防突巖層強(qiáng)度；而軟化的防突巖層在應(yīng)力、化學(xué)溶蝕、巖溶水壓力季節(jié)性變化等耦合作用下，巖溶裂隙增加，防突能力進(jìn)一步下降，甚至喪失。

2)巖溶水壓力對裂隙巖體的有效應(yīng)力作用

巖溶隧道防突巖層中普遍存在著孔隙、裂隙、層理等不同大小的非連續(xù)結(jié)構(gòu)面，甚至還有規(guī)模更大的斷層面，它們構(gòu)成了巖溶隧道防突巖層突水的潛在通道。巖體孔隙、裂隙中地下水的力學(xué)效應(yīng)表現(xiàn)為孔隙、裂隙中的水壓力作用和軟化作用，2種作用存在聯(lián)系，其綜合的力學(xué)效應(yīng)可采用下式表達(dá)[13]。

Δτ=σ(tanφ-tanφw)+ptanφw+c-cw

式中：Δτ為巖溶水引起的防突巖層抗剪強(qiáng)度降低值；σ為正應(yīng)力；p為巖溶水壓力；c，φ分別為浸水前巖體的黏聚力和內(nèi)摩擦角；cw，φw分別為浸水后巖體的黏聚力和內(nèi)摩擦角。

巖溶水壓力對裂隙巖體的力學(xué)作用主要體現(xiàn)在巖溶水對裂隙巖體有效應(yīng)力的控制上。巖溶隧道開挖過程中，開挖面內(nèi)巖體從母巖上脫離，包括防突巖層在內(nèi)的隧道圍巖因失去支撐作用而卸荷，深埋高地應(yīng)力巖溶隧道這一現(xiàn)象更為明顯。由上式可知，因卸荷而導(dǎo)致的正應(yīng)力σ下降，會進(jìn)一步增強(qiáng)巖溶水壓力p對防突巖層巖體強(qiáng)度的影響。

3)巖溶水壓力的水楔效應(yīng)

圖1為巖溶隧道開挖工程中常見的一種工況。在隧道掌子面前方有一小傾角的逆斷層F，隧道底板承受巖溶水壓力p的作用。開挖前隧道底板上作用著原始垂向地應(yīng)力σ0，此時(shí)斷層處于閉合狀態(tài)。施工時(shí)開挖范圍內(nèi)σ0被解除，而在開挖區(qū)四周分布著的支承壓力σ(x)部分作用于隧道底板上。隨著σ0的解除隧道底板產(chǎn)生一定程度的彈性恢復(fù)變形，即產(chǎn)生垂直方向的上升位移，斷層面之間的距離會因此而改變，這樣以來原本閉合的斷層將會張開，有巖溶承壓水趁機(jī)滲入，斷層中的水分別對斷層上下兩盤產(chǎn)生法向壓應(yīng)力，進(jìn)而引起兩斷層面的法向位移，促使承壓水進(jìn)一步向上導(dǎo)升。即整個過程中在壓力作用下巖溶水像楔子一樣沿著斷層面逐漸向上“擠”，這就是巖溶水壓力的水楔效應(yīng)。

圖1 巖溶水壓力在導(dǎo)通斷層中的水楔效應(yīng)示意

4)巖溶水流的沖刷擴(kuò)徑作用

水流的沖刷擴(kuò)徑作用是指在巖溶水流的沖刷下突水通道橫斷面尺寸不斷擴(kuò)大。該作用使防突巖層突水口不斷擴(kuò)大，是突水規(guī)模持續(xù)增大的內(nèi)因。

1.2 巖溶隧道突水類型

根據(jù)突水成因?qū)r溶隧道突水分為直接揭露型、純劈裂型和綜合破壞型3種類型[14]。

1)直接揭露型突水。直接揭露型突水主要是指隧道開挖直接揭露富水溶腔、地下暗河等引起的突水。該類突水源于隧道掘進(jìn)時(shí)對巖層的主動破壞。

2)純劈裂型突水。當(dāng)巖溶水壓力較高，其影響遠(yuǎn)大于圍巖應(yīng)力重分布的影響時(shí)，可忽略圍巖應(yīng)力重分布的影響。在巖溶水高壓作用下圍巖劈裂導(dǎo)致突水事故發(fā)生，即發(fā)生純劈裂型突水[15]。

3)綜合破壞型突水。綜合破壞型突水是指在巖溶水壓力與圍巖應(yīng)力重分布的共同作用下致使隧道與隱伏巖溶構(gòu)造間防突巖層破壞，進(jìn)而喪失隔水能力而發(fā)生的突水。根據(jù)水文地質(zhì)條件、防突巖層介質(zhì)類型及防突巖層裂隙發(fā)育情況，綜合破壞型突水又可細(xì)分為整體失穩(wěn)破壞、自然營力下的水力劈裂破壞、關(guān)鍵塊體失穩(wěn)破壞和滲透失穩(wěn)破壞。

2 防突巖層突水離散元分析

巖溶隧道突水機(jī)理分析屬于定性層面的研究，無法反映巖溶隧道防突巖層突水過程及該過程中圍巖變形場、滲流場的演化特征。本節(jié)以邊墻外賦存高壓富水溶腔這一類型的新建山嶺隧道為研究對象，通過離散元方法對此進(jìn)行分析。

2.1 計(jì)算模型及模擬過程

計(jì)算模型見圖2。左側(cè)隧道直墻段高6 m，頂部半圓拱直徑8 m，跨度8 m；右側(cè)富水隱伏溶洞直徑 10 m，橫斷面圓形；中間為3 m厚的防突巖層。取隧道底板中心為計(jì)算模型的坐標(biāo)原點(diǎn)，模型頂部距隧道底板45 m，底部距隧道底板30 m。模型水平方向長65 m，為溶洞直徑的6.5倍。

圖2 計(jì)算模型

根據(jù)隧道埋深情況將模型上部巖體重量換算成均布荷載施加在模型頂面。頂面為自由面，模型底面施加位移約束，模型左右側(cè)面施加應(yīng)力約束，側(cè)壓力系數(shù)取1.2；模型底面及左右側(cè)面均不透水。圍巖中分布2組相互正交的節(jié)理，1組與水平面(x軸正方向)夾角為15°，節(jié)理間距1 m；另1組與水平面(x軸正方向)夾角為105°，節(jié)理間距2 m。

表1 圍巖力學(xué)參數(shù)

表2 節(jié)理力學(xué)參數(shù)

2.2 模擬結(jié)果與分析

1)位移場變化時(shí)防突巖層的破壞情況

巖溶隧道防突巖層突水過程中位移場的變化見圖3?？梢姡核淼篱_挖后在圍巖應(yīng)力釋放初期防突巖層整體發(fā)生了較大位移，尤其是隧道側(cè)的防突巖層，由于位移約束消失向隧道內(nèi)部變形明顯，隧道周邊其他位置位移較小(step=700)；圍巖應(yīng)力繼續(xù)釋放，隧道側(cè)的防突巖層開始垮落，但整體尚處于穩(wěn)定狀態(tài)(step=2 200)； step=4 200 時(shí)防突巖層進(jìn)一步破壞，隧道左側(cè)拱頂和底板也發(fā)生明顯位移，溶腔左側(cè)頂部亦出現(xiàn)了不穩(wěn)定跡象，故隧道開挖后應(yīng)及時(shí)支護(hù)；應(yīng)力釋放到一定值時(shí)防突巖層破壞，突水通道基本形成，此時(shí)其他位置圍巖也出現(xiàn)了滑動跡象(step=6 200)，裂隙增大。

圖3 巖溶隧道防突巖層突水過程中位移場的變化

2)滲流場變化時(shí)巖溶水壓力的分布情況

巖溶隧道防突巖層突水過程中巖溶水壓力的分布見圖4?？梢姡涸趪鷰r應(yīng)力釋放初期巖溶水壓力主要分布在溶腔頂部和臨溶腔側(cè)的防突巖層中(step=700～2 200)，但step=2 200 時(shí)巖溶水壓力在防突巖層中的分布范圍明顯增加，且?guī)r溶水壓力在隧道底板中也有分布，反映出巖溶水已通過防突巖層滲透到隧道中，并且隧道底部巖層已和溶腔建立了水力聯(lián)系。隨著隧道圍巖應(yīng)力的繼續(xù)釋放(step=4 200)，巖溶水壓力分布范圍繼續(xù)增大，溶腔周圍特別是頂部已分布高壓巖溶水，隧道頂部巖層也和溶腔建立了水力聯(lián)系，此時(shí)防突巖層中已遍布巖溶水壓力，隧道底板及左墻腳巖層內(nèi)亦有巖溶水壓力分布，但壓力較小。step=6 200 時(shí)溶腔和隧道頂部水壓力分布范圍有所減小，且隧道頂部巖層與溶腔之間的水力聯(lián)系已中斷，隧道底板及左墻腳巖層中巖溶水壓力已和溶腔中水壓力相當(dāng)，突水通道已經(jīng)基本形成，防突巖層處于無阻水能力狀態(tài)，富水隱伏溶腔水位下降。

3 結(jié)論

1)巖溶隧道突水機(jī)理主要為巖溶水的軟化和溶蝕作用、巖溶水壓力對裂隙巖體的有效應(yīng)力作用、巖溶水壓力的水楔效應(yīng)和巖溶水流的沖刷擴(kuò)徑作用。

2)巖溶隧道突水可分為直接揭露型、純劈裂型和綜合破壞型。綜合破壞型又可進(jìn)一步細(xì)分為整體失穩(wěn)破壞、自然營力下的水力劈裂破壞、關(guān)鍵塊體失穩(wěn)破壞及滲透失穩(wěn)破壞。

3)隧道開挖后隨圍巖應(yīng)力逐漸釋放，防突巖層經(jīng)歷顯著位移、局部垮落和整體破壞3個階段，最終形成突水通道。在隧道開挖初期巖溶水壓力分布范圍逐漸增大，隧道周圍的巖體逐步和高壓富水溶腔建立了水力聯(lián)系；開挖后期溶腔水位下降，巖溶水壓力分布范圍減小。