干濕循環(huán)作用下裂隙巖體滲流特性演化規(guī)律研究

摘要：

為了研究水庫(kù)消落帶巖體的滲流特性，以烏東德庫(kù)區(qū)白云巖為研究對(duì)象，開展干濕循環(huán)水-巖作用試驗(yàn)。基于單裂隙滲流試驗(yàn)研究了干濕循環(huán)作用下裂隙巖體滲流特性演化規(guī)律，采用理論分析方法研究了裂隙巖體各向異性滲透特性。研究結(jié)果表明：當(dāng)圍壓恒定時(shí)，單裂隙滲流量隨滲透水壓力的增加線性增加；當(dāng)滲透水壓力恒定時(shí)，隨著圍壓的增加滲流量呈指數(shù)函數(shù)關(guān)系遞減；隨著水-巖作用周期的增加，裂隙滲流量呈先驟減小后緩慢增大的變化規(guī)律。裂隙巖體滲透特性存在明顯的各向異性特征，并且水-巖作用后主滲透系數(shù)和滲透張量減小，滲透橢圓半徑增加，滲透主方向保持不變。

關(guān)" 鍵" 詞：

裂隙巖體； 滲流特性； 滲流量； 各向異性； 滲透橢圓； 干濕循環(huán)； 烏東德庫(kù)區(qū)

中圖法分類號(hào)： TU45

文獻(xiàn)標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2024.07.026

收稿日期：

2023-08-16

；接受日期：

2023-11-03

基金項(xiàng)目：

中國(guó)三峽建工（集團(tuán)）有限公司科研項(xiàng)目（WDD/0569）;國(guó)家自然科學(xué)基金聯(lián)合基金項(xiàng)目（U2034203）;國(guó)家自然科學(xué)基金青年基金項(xiàng)目（52009067）

作者簡(jiǎn)介：

施" 炎，男，工程師，碩士，主要從事水利水電工程地質(zhì)勘查研究。E-mail：451435467@qq.com

通信作者：

孫旭曙，男，副教授，博士，主要從事水工巖石力學(xué)研究。E-mail：sunxs@ctgu.edu.cn

Editorial Office of Yangtze River. This is an open access article under the CC BY-NC-ND 4.0 license.

文章編號(hào)：1001-4179（2024） 07-0198-06

引用本文：

施炎，王團(tuán)樂，宛良朋，等.干濕循環(huán)作用下裂隙巖體滲流特性演化規(guī)律研究

［J］.人民長(zhǎng)江，2024，55（7）：198-203.

0" 引 言

水利水電工程在運(yùn)行期間庫(kù)水位頻繁升降，尤其是日調(diào)節(jié)抽水蓄能電站的上下水庫(kù)，水位升降高度可達(dá)20～40 m，由此導(dǎo)致庫(kù)岸邊坡滲流場(chǎng)的頻繁變化。與此同時(shí)，庫(kù)水呈弱酸或堿性，并含有多種化學(xué)物質(zhì)，可使巖石中的礦物成分與其發(fā)生以物理和化學(xué)過程為主的溶解或沉淀作用，這一過程也稱為水-巖作用；該作用進(jìn)而改變裂隙面的形貌特征和滲流特性，直接影響著庫(kù)岸邊坡滲流場(chǎng)的變化。岸坡滲流場(chǎng)的變化又將引起巖體應(yīng)力狀態(tài)的改變，從而進(jìn)一步影響庫(kù)岸邊坡的穩(wěn)定性。因庫(kù)水位周期性變化引起庫(kù)岸邊坡失穩(wěn)現(xiàn)象日益增多。因此，研究干濕循環(huán)作用下裂隙巖體滲流場(chǎng)的滲透特性對(duì)水利水電工程安全運(yùn)營(yíng)具有重要的指導(dǎo)作用。

近年來，庫(kù)岸邊坡水-巖作用引起了眾多學(xué)者的關(guān)注。如王士天等［1］總結(jié)了庫(kù)區(qū)水-巖作用的類型和特征，闡述了由水-巖作用導(dǎo)致的庫(kù)區(qū)地質(zhì)災(zāi)害問題。李克鋼［2］通過室內(nèi)試驗(yàn)研究水-巖作用下白云巖變形、強(qiáng)度、破壞特征等力學(xué)特性的變化規(guī)律。Hurowitz［3］、Alt-Epping［4］、鄧華鋒［5］等通過花崗巖、砂巖、頁(yè)巖及紅層軟巖等多種巖石的水-巖作用試驗(yàn)，發(fā)現(xiàn)水-巖作用導(dǎo)致巖石彈性模量、抗壓強(qiáng)度、黏聚力及內(nèi)摩擦角等發(fā)生不同程度的劣化，其力學(xué)性能的劣化直接影響庫(kù)岸邊坡的穩(wěn)定性。文獻(xiàn)［2-5］主要分析了水-巖作用對(duì)小尺寸巖塊力學(xué)性能的影響，而天然巖體中含有數(shù)量眾多的節(jié)理、裂隙等不連續(xù)面，地下水流經(jīng)該類不連續(xù)面將產(chǎn)生水-巖作用，從而進(jìn)一步影響巖體的變形、應(yīng)力狀態(tài)、滲流特性，且其宏觀效應(yīng)變化更顯著?；诖耍糠謱W(xué)者從應(yīng)力場(chǎng)和滲流場(chǎng)耦合角度研究了裂隙的滲流特性。如Walsh［6-7］、Tsang［8］、速寶玉［9］、王媛［10］等通過裂隙巖體滲流和應(yīng)力耦合作用建立了裂隙滲流量和應(yīng)力的關(guān)系。此外，趙宗勇等［11］通過數(shù)值模擬方法分析了裂隙過流面積對(duì)滲流量的影響，張修碩等［12］基于立方定律得到不同形式裂隙流量計(jì)算公式。段玲玲［13］、王珂［14］等學(xué)者通過水-巖作用下裂隙巖體滲流試驗(yàn)分析了水-巖作用對(duì)裂隙滲流特性的影響，但其研究過程中均以單一裂隙為研究對(duì)象，未進(jìn)一步研究大尺寸裂隙網(wǎng)絡(luò)滲流特性。因此，黃庭威［15］、青松勇［16］、王振偉［17］等學(xué)者開展了巖體裂隙網(wǎng)絡(luò)滲流各向異性特征研究，發(fā)現(xiàn)裂隙幾何參數(shù)是裂隙滲透張量和滲透主方向產(chǎn)生變化的主要因素，但其研究過程中缺少考慮庫(kù)岸邊坡水-巖作用條件對(duì)巖體裂隙網(wǎng)絡(luò)滲流特性的影響。

因此，本文以烏東德庫(kù)區(qū)三臺(tái)危巖體白云巖為研究對(duì)象，設(shè)計(jì)了干濕循環(huán)試驗(yàn)?zāi)M消落帶巖體“浸泡+自然風(fēng)干”的循環(huán)作用環(huán)境。通過單裂隙滲流試驗(yàn)獲得裂隙滲流特性變化規(guī)律，基于蒙特卡洛方法建立三臺(tái)危巖體裂隙網(wǎng)絡(luò)模型，并通過理論分析開展巖體裂隙網(wǎng)絡(luò)滲流各向異性特征研究。

1" 試樣制備及試驗(yàn)方案

1.1" 試樣制備

本文取烏東德庫(kù)區(qū)三臺(tái)危巖體白云巖為研究對(duì)象，通過X射線衍射儀獲得白云巖的主要礦物為白云石（MgCa（CO3）2）。根據(jù)規(guī)范和單裂隙水巖耦合儀對(duì)巖樣尺寸的要求，首先制備成50 mm（直徑）×100 mm（高度）的圓柱形標(biāo)準(zhǔn)試樣。經(jīng)巴西劈裂試驗(yàn)獲得具有完整裂隙面的單裂隙試樣，見圖1，通過三維形貌掃描獲得裂隙面分形維數(shù)為2.211，算術(shù)平均高度Sa為0.801 mm。為了避免裂隙面微細(xì)觀結(jié)構(gòu)及形貌特征的離散造成裂隙滲流量差異較大，故取一組試樣開展水-巖作用條件下單裂隙滲流試驗(yàn)。

1.2" 試驗(yàn)方案

1.2.1" 水-巖作用試驗(yàn)方案

三臺(tái)危巖體位于烏東德峽谷右岸。由三臺(tái)壩址預(yù)可行性階段地下水質(zhì)資料可知：工作區(qū)地下水類型為HCO3-Mg·Ca型，礦化度469 mg/L，pH為8.4，江水為HCO3-Mg·Ca型，礦化度250 mg/L，pH為8.3，具有微腐蝕性。因此，配制pH為8.5的飽和碳酸氫鈉溶液作為浸泡溶液，模擬巖樣所賦存的水環(huán)境條件。消落帶巖體長(zhǎng)期處于“浸泡+自然風(fēng)干”循環(huán)狀態(tài)，共設(shè)計(jì)了12期干濕循環(huán)試驗(yàn)過程，每期試驗(yàn)靜水浸泡20 d，自然風(fēng)干2 d。在初始狀態(tài)和第1、3、5、8、12周期開展單裂隙滲流試驗(yàn)。

1.2.2" 單裂隙滲流試驗(yàn)方案

借助自主研發(fā)的HYS-4-2型單裂隙水巖耦合儀開展單裂隙白云巖滲流試驗(yàn)，如圖2所示。試驗(yàn)設(shè)計(jì)了3個(gè)圍壓（1，2，3 MPa）和3個(gè)滲透水壓力（0.2，0.5，0.8 MPa），分別由液壓泵和氣控閥控制圍壓和滲透水壓力，控制精度分別不大于0.1 MPa和0.02 MPa，在初始狀態(tài)和第1，3，5，8，12周期記錄滲流量并分析其變化規(guī)律。水-巖作用試驗(yàn)與單裂隙滲流試驗(yàn)交替進(jìn)行，為了避免人工操作對(duì)裂隙面貼合度和裂隙錯(cuò)位的影響，裂隙面閉合后觀察試樣表面紋理連接情況如圖1所示，并在滲流試驗(yàn)中通過有機(jī)硅膠密封固定。

2" 單裂隙滲流特性分析

2.1" 滲透水壓力對(duì)滲流量的影響

干濕循環(huán)作用下不同試驗(yàn)周期單裂隙滲流量見表1。由表1可知，不同周期、不同圍壓條件下滲流量與滲透水壓力變化關(guān)系基本相似。因此，取第12周期數(shù)據(jù)繪制滲流量與滲透水壓力變化曲線，如圖3所示。結(jié)果表明：當(dāng)圍壓恒定時(shí)，隨著滲透水壓力增加單裂隙滲流量線性增加，用線性函數(shù)y=Ax+B能較好擬合，其中，系數(shù)A隨著圍壓的增大而減小。當(dāng)滲透水壓由0.2 MPa增加至0.8 MPa時(shí)，滲流量增加3.1～7.4倍。

2.2" 圍壓對(duì)滲流量的影響

從表1取第12周期數(shù)據(jù)繪制圍壓與滲流量的變化曲線見圖4。結(jié)果表明：當(dāng)滲透水壓力恒定時(shí)，隨著圍壓增加單裂隙滲流量減小，其變化趨勢(shì)可用指數(shù)函數(shù)y=Ce-Dx較好擬合，且隨著滲透水壓力的增加系數(shù)C增大。當(dāng)圍壓由1 MPa增加到3 MPa時(shí)，滲流量減小為初始流量的2%～31%，與滲透水壓力相比，滲流量隨圍壓變化最為顯著。

2.3" 試驗(yàn)周期對(duì)滲流量的影響

水-巖作用過程中，滲流量隨試驗(yàn)周期變化較為明顯，以3 MPa圍壓作用下單裂隙滲流量變化為例，繪制滲流量與試驗(yàn)周期的變化曲線見圖5。結(jié)果表明：水巖作用初期，滲流量出現(xiàn)驟減；隨后，滲流量隨試驗(yàn)周期的增加呈線性增加，與第一周期相比，第12周期的滲流量分別增加了48.92%～72.69%。

對(duì)初始狀態(tài)及第12周期水-巖作用后裂隙表面開展SEM掃描電鏡測(cè)試，獲得微細(xì)觀結(jié)構(gòu)圖像如圖6所示。由測(cè)試結(jié)果可知：經(jīng)多次干濕循環(huán)作用可使得裂隙表面礦物軟化溶解，微細(xì)觀孔隙增加，導(dǎo)致裂隙表面形貌特征向復(fù)雜多樣化變化；在圍壓作用下進(jìn)一步改變了裂隙中水流通道，最終導(dǎo)致水巖作用后單裂隙滲流量逐漸增加。

3" 干濕循環(huán)作用下裂隙巖體各向異性滲透特性

3.1" 計(jì)算模型及邊界條件

從烏東德庫(kù)區(qū)三臺(tái)危巖體露頭、平硐獲取兩組主要裂隙參數(shù)見表2。其中間距和斷距服從正態(tài)分布?；贛onte-Carlo方法生成裂隙網(wǎng)絡(luò)模型并對(duì)其進(jìn)行連通處理，如圖7所示，模型尺寸為20 m×20 m。為研究裂隙巖體各向異性滲透特性，在裂隙網(wǎng)絡(luò)模型中心取10 m×10 m的研究區(qū)域，并以中心點(diǎn)為圓心，依次旋轉(zhuǎn)15°、30°、45°、60°、75°獲得不同方向的裂隙網(wǎng)絡(luò)模型。

在裂隙網(wǎng)絡(luò)滲流計(jì)算中假定地下水只能在裂隙中流動(dòng)，文中以初始狀態(tài)和第12周期（圍壓1 MPa、滲透水壓力0.8 MPa）測(cè)量的滲流量為基礎(chǔ)，通過立方定律計(jì)算獲得水-巖作用前后單裂隙的水力傳導(dǎo)系數(shù)，滲流計(jì)算邊界條件及水力傳導(dǎo)系數(shù)如下：邊界條件為入口壓力水頭80 m，出口壓力水頭0 m，不透水邊界流量0，初始狀態(tài)裂隙水力傳導(dǎo)系數(shù)1.83×10-3 m/s，水巖作用后裂隙水力傳導(dǎo)系數(shù)5.69×10-4 m/s。借助Geo-studio軟件開展?jié)B流計(jì)算，假定研究區(qū)域內(nèi)裂隙網(wǎng)絡(luò)滲流量與同尺寸連續(xù)介質(zhì)滲流量相同，即可通過連續(xù)介質(zhì)模型計(jì)算等效滲透系數(shù)：

Q=kJfA

（1）

式中：k為等效滲透系數(shù)，A為截面面積，Jf為水力坡度。

3.2" 裂隙滲透主方向及各向異性特征研究

假設(shè)坐標(biāo)系xoy下滲透系數(shù)為k，坐標(biāo)系x′oy′下滲透系數(shù)為k′，θ為坐標(biāo)軸的夾角（圖8），則坐標(biāo)軸夾角的正弦和余弦值分別為：lx′x=cosθ，ly′y=cosθ，lx′y=sinθ，ly′x=-sinθ，滲透系數(shù)張量K′和K的分量之間存在如下關(guān)系［12］：

ki′j′=kijli′ilj′j

（2）

設(shè)主滲透系數(shù)k1、k2的方向與坐標(biāo)系x軸和y軸方向一致，則θ方向的滲透系數(shù)為

kθ=k1cos2θ+k2sin2θ

（3）

用橢圓方程表示為

cos2θa2+sin2θb2=1r2

（4）

設(shè)1a2=k1，1b2=k2，1r2=kθ，即旋轉(zhuǎn)角θ與kθ滿足橢圓曲線關(guān)系。1k1、1k2分別為橢圓的主軸長(zhǎng)度［12］，橢圓長(zhǎng)軸與x軸正向的夾角為滲透主方向。由公式（5）得滲透張量的各個(gè)分量：

Kxx=12k1+k2+12k1-k2cos2θKyy=12k1+k2-12k1-k2cos2θKxy=Kyx=12k1-k2sin2θ

（5）

式中：θ為滲透主方向。綜上可得裂隙巖體二階滲透張量為

K=KxxKxy

KyxKyy

（6）

根據(jù)不同旋轉(zhuǎn)角模型對(duì)應(yīng)的滲透系數(shù)及橢圓半徑關(guān)系，在極坐標(biāo)和直角坐標(biāo)系中繪制滲透系數(shù)和橢圓半徑的變化曲線如圖9所示。由橢圓方程得主滲透系數(shù)

k1=1.01×10-6 m/s，k2=1.02×10-5 m/s，滲透主方向?yàn)?0.2°，滲透張量為K=6.46-4.54-4.544.79×10-6 m/s。

同理，得水-巖作用后滲透系數(shù)及橢圓半徑的變化曲線如圖10所示。由橢圓方程得主滲透系數(shù)k1=3.13×10-7 m/s，k2=3.16×10-6 m/s，滲透主方向?yàn)?0.2°，滲透張量為K=2.00-1.40-1.401.48×10-6 m/s。綜上可知，裂隙巖體具有明顯的各向異性特征，并且水-巖作用后滲透橢圓半徑增加，主滲透系數(shù)和滲透張量減小，而滲透主方向保持不變。

4" 結(jié) 論

（1） 當(dāng)圍壓恒定時(shí)，單裂隙滲流量隨滲透水壓的增加線性增加；當(dāng)滲透水壓恒定時(shí)，隨著圍壓增加滲流量呈指數(shù)函數(shù)關(guān)系遞減；隨著試驗(yàn)周期的增加，裂隙滲流量呈先驟減后緩慢增加并趨于穩(wěn)定的變化規(guī)律。

（2） 無論水-巖作用前后，裂隙網(wǎng)絡(luò)滲透特性均表現(xiàn)出明顯的各向異性特征，相比于初始狀態(tài)而言，水-巖作用后裂隙巖體主滲透系數(shù)和滲透張量減小，滲透橢圓半徑增加，滲透主方向保持不變。

參考文獻(xiàn)：

［1］" 王士天，劉漢超，張倬元，等.大型水域水巖相互作用及其環(huán)境效應(yīng)研究［J］.地質(zhì)災(zāi)害與環(huán)境保護(hù)，1997（1）：70-90.

［2］" 李克鋼.水巖物理作用下巖石力學(xué)特性研究［M］.北京：冶金工業(yè)出版社，2016.

［3］" HUROWITZ J A，F(xiàn)ISCHER W W，Contrasting styles of water-rock interaction at the mars exploration rover landing sites［J］.Geochimica Et Cosmochimica Acta，2014，127：25-38.

［4］" ALT-EPPING，DIAMOND L W，HARING，M O，et al.Prediction of water-rock interaction and porosity evolution in a granitoid-hosted enhanced geothermal system，using constraints from the 5 km Basel-1 well ［J］.Applied Geochemistry，2013，38：121-133.

［5］" 鄧華鋒，支永艷，段玲玲，等.水-巖作用下砂巖力學(xué)特性及微細(xì)觀結(jié)構(gòu)損傷演化［J］.巖土力學(xué)，2019，40（9）：3447-3456.

［6］" WALSH J B.Effect of pore pressure and confining pressure on fracture permeability［J］.Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstract，1981，18（5）：429-435.

［7］" WALSH J B，GROSENBAUGH M A.A new model for analyzing the effect of fractures on compressibility［J］.Journal of Geophysical Research：Solid Earth，1979，84（B7）：3532-3536.

［8］" TSANG Y W，WITHERSPOON P A.Hydromechanical behavior of a deformable rock fracture subject to normal stress［J］.Journal of Geophysical Research：Solid Earth，1981，86（B10）：9287-9298.

［9］" 速寶玉，張文捷，盛金昌，等.滲流-化學(xué)溶解耦合作用下巖石單裂隙滲透特性研究［J］.巖土力學(xué)，2010，31（11）：3361-3366.

［10］王媛.單裂隙面滲流與應(yīng)力的耦合特性［J］.巖石力學(xué)與工程學(xué)報(bào)，2002（1）：83-87.

［11］趙宗勇，宛良朋，朱思軍.劈裂砂巖過流面積對(duì)滲流的影響研究［J］.人民長(zhǎng)江，2012，43（3）：88-91.

［12］張修碩，王志良，張瑞，等.巖體中不同形式裂隙對(duì)滲流壓力差損耗的影響［J］.人民長(zhǎng)江，2017，48（5）：71-78.

［13］段玲玲，鄧華鋒，齊豫，等.水-巖作用下單裂隙灰?guī)r滲流特性演化規(guī)律研究［J］.巖土力學(xué)，2020，41（11）：3671-3679.

［14］王珂，盛金昌，郜會(huì)彩，等.應(yīng)力-滲流侵蝕耦合作用下粗糙裂隙滲流特性研究［J］.巖土力學(xué)，2020，41（增1）：30-40.

［15］黃庭威.裂隙巖體REV尺寸及等效滲透系數(shù)影響因素研究［D］.合肥：合肥工業(yè)大學(xué)，2021.

［16］青松勇.隧道裂隙圍巖各向異性滲流場(chǎng)特性研究［D］.成都：西南交通大學(xué)，2021.

［17］王振偉，馬克，田洪圓，等.基于RFPA2D-Flow軟件對(duì)裂隙巖體滲透特性表征單元體的研究［J］.煤炭學(xué)報(bào)，2019，44（10）：3012-3021.

（編輯：鄭 毅）

Seepage characteristics of fractured rock mass under dry-wet cycling action

SHI Yan1，WANG Tuanle1，WAN Liangpeng2，SUN Xushu3，ZUO Xiaopeng3，GUO Xiaoping3

（1.Three Gorges Survey and Research Institute Co.，Ltd.，Wuhan 430074，China；" 2 China Three Gorges Construction Group Co.， Ltd.， Chengdu 610041，China；" 3.Key Laboratory of Geological Hazards on Three Gorges Reservoir Area of MoE，China Three Gorges University，Yichang 443002，China）

Abstract：

In order to study the seepage characteristics of fractured rock mass in the water-level fluctuation zone of reservoirs，the dry-wet cycling water-rock interaction test was carried out on dolomite in Wudongde Reservoir area.Based on the single fissure seepage test，the evolution law of seepage characteristics of fractured rock mass under the action of dry-wet cycling was studied，and the anisotropic seepage characteristics of fractured rock mass were studied by theoretical analysis method.The results showed that when the confining pressure was constant，the seepage flow of single crack increased linearly with the increase of seepage water pressure.When the seepage water pressure was constant，the seepage flow decreased exponentially with the increase of confining pressure.With the increasing of water-rock interaction period，the fracture seepage flow decreased suddenly and then slowly increased.The permeability characteristics of fractured rock mass had obvious anisotropy，the principal permeability coefficient and permeability tensor decreased，the radius of permeability ellipse increased，and the principal permeability direction remains unchanged after water-rock interaction.

Key words：

fractured rock； seepage characteristics； seepage flow； anisotropy； permeability ellipse； dry-wet cycling； wudongde Reservoir area