劉志春， 高新強， 朱永全

(石家莊鐵道大學(xué)土木工程學(xué)院，河北 石家莊 050043)

0 引言

隨著可持續(xù)發(fā)展戰(zhàn)略的深化，隧道工程所引起的地下水環(huán)境問題逐漸引起人們的高度重視，隧道工程建設(shè)與運營常導(dǎo)致地下水位下降、地下水資源流失及水文循環(huán)失衡，從而引起生態(tài)及生活環(huán)境惡化［1-3］。對于隧道工程與地下水的相互作用問題，以往研究大多集中在地下水對隧道工程的影響及對策，如防排水原則、結(jié)構(gòu)形式、外水壓力處理及人工降水引起的地面沉降等研究［2-5］，而隧道工程對地下水環(huán)境的影響方面研究較少，且大多限于隧道涌水量的預(yù)測研究方面，個別開展的隧道工程對地下水環(huán)境影響研究也僅處于定性及半定量階段［1，6，7］。

在如何有效地評價隧道對地下水環(huán)境的影響問題上，建立切實可行、行之有效的裂隙巖體滲流模型顯得尤為重要。裂隙巖體的滲流模型，根據(jù)巖體裂隙特征和裂隙性質(zhì)的把握準(zhǔn)確程度，主要為等效孔隙(連續(xù))介質(zhì)模型、雙重介質(zhì)模型和網(wǎng)絡(luò)裂隙介質(zhì)模型等［8］。然而，裂隙巖體本身具有復(fù)雜性，其滲流特征具有多樣性，計算模型如要完全反映裂隙巖體的滲流特征基本上是不可能的。因此采用等效連續(xù)介質(zhì)模型，旨在研究如何進(jìn)行模型假定、邊界條件選取、模型材料選擇、參數(shù)選取等，才能較好地反映裂隙巖體隧道的滲流特征。

1 裂隙巖體的滲流特征

(1)巖體的裂隙性。裂隙介質(zhì)是相對于孔隙介質(zhì)和溶穴介質(zhì)，是指含裂隙水的巖體。裂隙巖體由于裂隙發(fā)育、分布的方向性和不均勻性，其滲透具有明顯的各向異性和不均一性特點。

(2)巖體的滲透性。水只能沿連通孔隙或裂隙滲透，巖體滲透性大小主要取決于巖石孔隙的大小、方向及其相互連通情況，用滲透系數(shù)來衡量。

(3)巖體滲流的不均勻性。是指在巖體空間內(nèi)，不同位置上滲透系數(shù)大小不同。宏觀巖體滲流不均勻性主要是由于不連續(xù)面的發(fā)育程度的差異，通常表現(xiàn)為分帶性和成層性。

(4)巖體的各向異性。由于巖體介質(zhì)具有非均質(zhì)各向異性，要用滲透張量來描述巖體介質(zhì)各個方向上的不同滲透性能。通常在節(jié)理裂隙密集展布的方向上，巖體的滲透性占主導(dǎo)優(yōu)勢。一般經(jīng)過坐標(biāo)變換的各向異性為非軸向各向異性。

除此，裂隙巖體還具有非飽和性及耦合特性等［8-9］。

2 裂隙巖體滲流模型

(1)等效連續(xù)介質(zhì)模型是把巖塊-裂隙系統(tǒng)等效成連續(xù)介質(zhì)，用連續(xù)介質(zhì)理論描述滲流方程，認(rèn)為水流服從Darcy定律，可直接應(yīng)用孔隙介質(zhì)滲流理論，操作性好。

(2)雙重介質(zhì)模型也是一種連續(xù)介質(zhì)模型，認(rèn)為巖塊孔隙系統(tǒng)(包括微裂隙)和裂隙系統(tǒng)連續(xù)地充滿整個區(qū)域，即把裂隙巖體看作是具有不同水力參數(shù)的兩種連續(xù)介質(zhì)的疊加體，兩種連續(xù)介質(zhì)中的滲流場均建立在Darcy定律的基礎(chǔ)上，并依據(jù)兩種介質(zhì)間的水交換項來聯(lián)立求解各自的滲流場。

(3)網(wǎng)絡(luò)裂隙介質(zhì)模型認(rèn)為巖塊本身不透水，整個地下水運動是通過裂隙網(wǎng)絡(luò)來進(jìn)行的。

由于裂隙本身的復(fù)雜性，后兩種模型在解決實際工程問題時均有其局限性。

3 流體偏微分方程

模型采用GeoStudio軟件SEEP/W模塊，該模塊用于分析土體或巖體地下水滲流和超孔隙水壓力消散問題［10-11］。二維滲流偏微分方程通?？杀硎鋈缦?/p>

式中，H為總水頭;Kx為x方向的水力滲透系數(shù);Ky為y方向的水力滲透系數(shù);Q為應(yīng)用邊界通量;θ為單位土體中含水量的變化;t為時間。

該方程說明了在某一時刻點流進(jìn)和流出一個單元體流量通量的差等于土體中含水量的變化。更具體的說，單位時間內(nèi)，在x和y方向及外表面方向通量總的變化速率等于單位土體中含水量的變化率。

體積水含量的變化取決于土體應(yīng)力狀態(tài)和特性的變化。在飽和土和非飽和土的條件下，應(yīng)力狀態(tài)用兩種狀態(tài)變量來描述。這些應(yīng)力狀態(tài)變量分別是(σ－ua)和(σ－uw)，這里σ是總應(yīng)力，ua空隙氣壓力，uw空隙水壓力。

SEEP/W程序假設(shè)總應(yīng)力恒定，也就是說，對土體既沒有加載也沒有卸載。SEEP/W還假定:在瞬態(tài)處理過程中，空隙氣壓力為定值且等于大氣壓力。這意味著(σ－ua)保持為定值，對體積含水量的變化沒有影響。含水量變化相應(yīng)地僅僅取決于(σ－uw)狀態(tài)變量的變化，其中，ua為定值，含水量的變化僅是于空隙水壓力的函數(shù)。結(jié)果是，含水量的變化僅與空隙水壓變化有關(guān)，如式(2)所示

式中，mw為含水量曲線的斜率?？偹^H可以這樣定義

式中，uw為空隙水壓;γw為水容重;y為位置水頭。

4 裂隙巖體隧道非均質(zhì)、各向異性等效滲流模型

采用SEEP/W模擬軟件建立等效連續(xù)介質(zhì)水文地質(zhì)模型，為增加模型擬真性，在模型中考慮巖層的非均質(zhì)性、非軸向各向異性［12］。

4.1 模型材料

模型材料中按實際地層厚度取值，考慮材料的非均質(zhì)性。根據(jù)結(jié)構(gòu)面分布情況，考慮平行于結(jié)構(gòu)面方向的滲透系數(shù)為Ky'，垂直于結(jié)構(gòu)面方向的滲透系數(shù)為Kx'，定義各向異性特征系數(shù)Kr=Ky'/Kx'。Kr的取值與節(jié)理發(fā)育程度緊密相關(guān)，結(jié)合以往經(jīng)驗，Ⅱ、Ⅲ、Ⅳ、Ⅴ級圍巖Kr分別取值為10、50、100、500。

計算斷面選取為垂直于隧道軸向橫斷面，巖層結(jié)構(gòu)面真傾角轉(zhuǎn)化成結(jié)構(gòu)面在計算斷面處的視傾角γ=arctan(tan αsin(90°－β))，其中，α為巖層結(jié)構(gòu)面(真)傾角;β為巖層結(jié)構(gòu)面傾向。

4.2 材料參數(shù)

對于穩(wěn)態(tài)分析，材料參數(shù)主要為包括滲透系數(shù)K、結(jié)構(gòu)面視傾角γ、各向異性特征系數(shù)Kr。圍巖滲透系數(shù)K取值根據(jù)施工期間實測涌水量，以及水井水位監(jiān)測結(jié)果通過經(jīng)驗公式及數(shù)值模型反演確定。帷幕注漿及徑向注漿所形成的注漿圈滲透系數(shù)按圍巖水泥注漿效果試驗經(jīng)驗統(tǒng)計取值。噴射混凝土滲透系數(shù)按設(shè)計抗?jié)B等級、實際施工效果試驗和經(jīng)驗統(tǒng)計按P2抗?jié)B等級取值。二次襯砌模筑混凝土滲透系數(shù)按設(shè)計抗?jié)B等級P8取值。盲管滲透系數(shù)按砂礫石滲透系數(shù)取值。結(jié)構(gòu)面視傾角γ根據(jù)實際結(jié)構(gòu)面傾向及傾角計算確定。Kr按圍巖分級確定。

4.3 模型邊界取值范圍

計算模型上邊界為隧道施工前初始地下水位線，下邊界根據(jù)隧道埋深情況取隧道底面以下50～100 m。左、右邊界根據(jù)地形情況、地質(zhì)情況選擇確定，其選取原則為:

(1)分水嶺位置作為模型左、右邊界;

(2)斷層結(jié)構(gòu)面位置做為模型的左、右邊界;

(3)河流及水庫位置做為模型的左、右邊界;

(4)如地面較平坦，無上述各種邊界，可根據(jù)影響半徑的經(jīng)驗公式［13］取值，即R=215.5+510.5K，其中，R為隧道涌水影響半徑;K為巖體滲透系數(shù)。

4.4 邊界條件

穩(wěn)態(tài)分析中，左、右邊界采用定水頭邊界條件，其水頭等于邊界處水位高度;下邊界為不透水邊界。

上邊界，即隧道施工之前的初始地下水線，隧道施工排水后，其水頭(負(fù)壓力水頭)是未知的，其流量在地形及地下水位起伏較大時也是未知的，只有在地形及地下水位水平時上邊界流量均為0(不考慮降雨入滲因素)。此時，對于上邊界可采用q=0(或Q=0，q為單元流量，Q為節(jié)點流量)潛在浸潤面檢驗邊界，也即通過迭代求解找出浸潤面位置。當(dāng)考慮降雨入滲時，可按當(dāng)?shù)啬杲涤炅亢徒涤耆霛B系數(shù)，計算出每天的降雨入滲量，然后可設(shè)q=定值的潛在浸潤面檢驗邊界。對于隧道洞壁邊界條件可設(shè)為壓力水頭為0，或同樣采用檢驗的潛在浸潤面q=0邊界，二者是等效的。當(dāng)然也可以根據(jù)需要設(shè)為固定流量邊界。

5 應(yīng)用實例驗證

石板山隧道位于石太客運專線河北省井陘縣之內(nèi)，設(shè)計為單洞雙線隧道，最大埋深322 m，全長7505 m，隧道開挖最大跨度近15 m，起訖里程DIK49+585～DIK57+090，是全線的重點控制工程。隧道工程區(qū)屬于太行山東麓基巖低山區(qū)，屬于北方干旱、半干旱區(qū)，降水量不大，地表植被較發(fā)育。隧道范圍地層主要巖性為第四系上更新統(tǒng)坡洪積層(Q3dl+pl)新黃土，碎石類土，沖洪積層(Q3al+pl)粉質(zhì)黏土;寒武系中統(tǒng)徐莊組(ε2x)泥質(zhì)條帶灰?guī)r、頁巖、石灰?guī)r，寒武系下統(tǒng)毛莊組(ε1mz)白云巖、白云質(zhì)泥巖;元古界長城系高于莊組(Chg)白云巖、長石石英砂巖，大紅峪組(Chd)石英巖狀砂巖、長石石英砂巖夾粉砂巖，甘陶河群變長石石英砂巖、云母片巖、綠泥片巖;太古界南營巖組(Arn)花崗片麻巖;局部為燕山期(δu)輝綠巖侵入體。地下水分為變質(zhì)巖風(fēng)化裂隙水、碳酸鹽巖巖溶裂隙水、松散層孔隙水，以及碎屑巖夾碳酸鹽巖孔隙裂隙巖溶水。地下水主靠大氣降水直接滲入及基巖山區(qū)地下水徑流的側(cè)向補給。

隧道區(qū)范圍內(nèi)分布有多處泉、井和水庫。隧道施工過程中，曾造成井泉干涸，水庫水位下降，影響到居民的生產(chǎn)生活用水。施工中采取帷幕注漿、補充徑向注漿、全包式柔性防水層、提高噴射混凝土防水性能、及時支護封閉等系列防排水技術(shù)體系。合理評價隧道對地下水資源的影響顯得尤為重要。

沿隧道縱向選取14個代表性斷面進(jìn)行滲流分析，建立計算模型分析隧道施工引起的地下水滲流規(guī)律，并得出地下水位變化及地下水資源流失量，如圖1為DIK52+300斷面計算模型。在未考慮降雨季節(jié)變化情況下，得隧道運營階段分區(qū)段排水量分布如圖2所示。

圖1 DIK52+300計算模型

圖2 成洞運營階段分區(qū)段隧道排水量

由圖2，DIK51+800、52+300代表區(qū)段按實際施作帷幕注漿考慮，其余區(qū)段采用常規(guī)復(fù)合式襯砌，得運營階段隧道總排水量為1669.74 m3/d。從圖中區(qū)段排水量曲線可見，施作帷幕注漿區(qū)段排水量明顯下降，說明帷幕注漿是防止地下水資源流失的有效手段。2008年8月28日和10月9日雨季兩次測試的隧道總排水量分別為2082.24 m3/d和1607.04 m3/d。計算結(jié)果與實測結(jié)果較為吻合，但計算結(jié)果未考慮降雨季節(jié)變化，兩次實測為雨季測試，故計算結(jié)果實際上比實測結(jié)果略偏大。

6 結(jié)語

(1)采用非均質(zhì)、各向異性等效連續(xù)介質(zhì)模型模擬裂隙巖體滲流規(guī)律及地下水資源流失量，并應(yīng)用于實際工程中，計算結(jié)果較為符合實際。

(2)帷幕注漿區(qū)段排水量明顯下降，是防止地下水資源流失的有效手段。

(3)計算石板山隧道運營期間隧道排水量1669.74 m3/d，說明隧道建設(shè)對區(qū)域地下水資源有一定影響，運營期間應(yīng)采用有效措施進(jìn)行地下水資源的合理再利用。

［1］劉丹，楊立中，于蘇俊.華鎣山隧道排水的生態(tài)環(huán)境問題及效應(yīng)［J］.西南交通大學(xué)學(xué)報，2001，36(3):308-313.

［2］劉志春.地鐵施工與地下水滲流的環(huán)境效應(yīng)分析及應(yīng)用［D］.石家莊:石家莊鐵道學(xué)院交通工程系，2003.

［3］高新強.高水壓山嶺隧道襯砌水壓力分布規(guī)律研究［D］.成都:西南交通大學(xué)土木工程學(xué)院，2005.

［4］黃濤，楊立中.滲流與應(yīng)力耦合環(huán)境下裂隙圍巖隧道涌水量的預(yù)測研究［J］.鐵道學(xué)報，1999，21(6):75-80.

［5］王秀英.巖溶隧道堵水限排襯砌外水壓力及結(jié)構(gòu)設(shè)計研究［D］.北京:北京交通大學(xué)土木建筑工程學(xué)院，2005.

［6］蔣忠信.隧道工程與水環(huán)境的相互作用［J］.巖石力學(xué)與工程學(xué)報，2005，24(1):121-127.

［7］肖了林.隧道修長對區(qū)域地下水位的影響分析［J］.公路隧道，2000(1):43-52.

［8］周志芳.裂隙介質(zhì)水動力學(xué)原理［M］.北京:高等教育出版社，2007.

［9］林學(xué)鈺，廖資生，趙勇勝，等.現(xiàn)代水文地質(zhì)學(xué)［M］.北京:地質(zhì)出版社，2005.

［10］中仿科技.專業(yè)高效的巖土工程設(shè)計分析軟件 GeoStudio［EB/OL］.［2011-2-26］.http://geostudio.cntech.com.cn/intro.html.

［11］JOHN KRAHN.Seepage Modeling with SEEP/W［EB/OL］.［2011-2-26］.http://www.geo-slope.com.

［12］劉志春.干旱山區(qū)隧道施工地下水環(huán)境保護技術(shù)研究［R］.石家莊:石家莊鐵道大學(xué)，2010.

［13］朱大力，李秋楓.預(yù)測隧道涌水量的地下水動力方法［J］.鐵路地質(zhì)與路基，1995(4):32-37.