基于三維隨機(jī)裂隙網(wǎng)絡(luò)的壩基多孔分序灌漿數(shù)值模擬

李瑞金,趙夢琦,王曉玲,祝玉珊

(天津大學(xué)水利工程仿真與安全國家重點(diǎn)實(shí)驗(yàn)室,天津300072)

基于三維隨機(jī)裂隙網(wǎng)絡(luò)的壩基多孔分序灌漿數(shù)值模擬

李瑞金,趙夢琦,王曉玲,祝玉珊

(天津大學(xué)水利工程仿真與安全國家重點(diǎn)實(shí)驗(yàn)室,天津300072)

大壩基巖內(nèi)節(jié)理裂隙數(shù)量眾多、分布隨機(jī)性強(qiáng)。針對裂隙灌漿數(shù)值模擬研究現(xiàn)狀,在耦合多孔的三維隨機(jī)裂隙網(wǎng)絡(luò)巖體灌漿模型的基礎(chǔ)上,提出了多孔分序灌漿實(shí)現(xiàn)方法,并結(jié)合實(shí)際工程,進(jìn)行了考慮隙寬隨機(jī)分布的裂隙灌漿數(shù)值模擬研究,分析了灌漿孔序?qū){液擴(kuò)散壓力場和單寬流量的影響規(guī)律。通過模擬灌漿量和灌漿時間與實(shí)際值的對比分析,驗(yàn)證了該方法的可靠性。

灌漿工程；數(shù)值模擬；三維隨機(jī)裂隙網(wǎng)絡(luò)；隙寬隨機(jī)分布；多孔分序

0 引 言

在水利水電工程中[1],大壩基巖內(nèi)通常含有數(shù)量眾多、隨機(jī)分布的節(jié)理裂隙,構(gòu)成了地下滲流通道,威脅水電站、水庫以及下游地區(qū)人民的生命和財(cái)產(chǎn)安全。為防止壩基滲漏,通常采用帷幕灌漿,構(gòu)建防滲帷幕。由于灌漿工程的隱蔽性和巖體裂隙分布的復(fù)雜性,裂隙巖體灌漿理論的研究目前仍滯后于工程實(shí)踐的進(jìn)展[2- 4]。隨著計(jì)算機(jī)技術(shù)的發(fā)展,數(shù)值模擬技術(shù)逐步成為研究裂隙巖體灌漿問題的主要方法。開展裂隙巖體灌漿的數(shù)值模擬研究,對指導(dǎo)水電工程大壩基巖裂隙巖體灌漿工程具有重要理論與現(xiàn)實(shí)意義。

在裂隙灌漿數(shù)值模擬方面,國內(nèi)外學(xué)者都進(jìn)行了一定的研究?；诙S離散元軟件UDEC,Omid等[5]模擬了漿液在兩組規(guī)則相交裂隙內(nèi)的擴(kuò)散過程和灌漿參數(shù)對漿液擴(kuò)散半徑的影響分析,并建議采用3DEC進(jìn)行空間裂隙巖體灌漿數(shù)值模擬；Xu等[6]模擬了多個灌漿孔在平面裂隙網(wǎng)絡(luò)內(nèi)的擴(kuò)散規(guī)律；范佳俊[7]采用蒙特卡洛方法構(gòu)建了二維離散裂隙網(wǎng)絡(luò),并根據(jù)模擬結(jié)果提出了最優(yōu)灌漿參數(shù)；基于FEM耦合VOF方法,Chen等[8]利用立方定律模擬了漿液在劈裂縫中的擴(kuò)散過程；基于Comsol軟件,劉健等[4]進(jìn)行了單一水平裂隙的灌漿數(shù)值模擬和試驗(yàn)的對比研究；趙鵬[9]建立了三維交叉裂隙模型,研究了漿液在交叉裂隙中的擴(kuò)散形態(tài)和壓力場分布規(guī)律；張慶松[10]、張超[11]、俞文生[12]等同樣針對單一裂隙進(jìn)行了數(shù)值模擬研究。上述裂隙灌漿數(shù)值模擬主要針對單裂隙和平面裂隙網(wǎng)絡(luò),且多以單孔灌漿為研究對象,缺少考慮多孔分序、隙寬隨機(jī)分布的三維隨機(jī)裂隙巖體灌漿數(shù)值模擬研究。

本文在耦合多孔的三維隨機(jī)裂隙網(wǎng)絡(luò)巖體灌漿模型的基礎(chǔ)上,提出多孔分序灌漿實(shí)現(xiàn)方法,結(jié)合實(shí)際工程,進(jìn)行了考慮隙寬隨機(jī)分布的裂隙灌漿數(shù)值模擬研究,重點(diǎn)分析了灌漿孔序?qū){液擴(kuò)散壓力場和單寬流量的影響規(guī)律,通過模擬灌漿量和灌漿時間與實(shí)際值的對比分析,驗(yàn)證了該方法的可靠性。

表1 基巖裂隙統(tǒng)計(jì)結(jié)果

1 數(shù)學(xué)模型

1.1 漿液擴(kuò)散模型

水泥基漿液流變性質(zhì)復(fù)雜,根據(jù)水灰比(W/C)一般可分為牛頓型(W/C>1.0)、賓漢姆型(0.8

(1)

基于光滑平行板裂隙模型,對牛頓型漿液,當(dāng)漿液在裂隙內(nèi)的流動為層流狀態(tài)時,裂隙中的漿液運(yùn)動滿足Navier-Stokes方程[14-15],考慮無滑移邊界條件,光滑裂隙的單寬流量可采用立方定律描述。實(shí)際情況下,裂隙面具有一定的粗糙度,通常采用流量等效原理,用等效水力隙寬uh代替粗糙裂隙物理隙寬um,其關(guān)系可表示為

uh=fum

(2)

式中,f為粗糙度修正因子,反映了裂隙面粗糙度對漿液流動的影響。因此,漿液擴(kuò)散模型為

(3)

式中,q為單寬流量；g為重力加速度；ρ為漿液密度；J為水力梯度。

1.2 漿液量統(tǒng)計(jì)

根據(jù)模擬灌漿流量與灌漿時間的關(guān)系曲線,采用積分原理統(tǒng)計(jì)灌漿量。計(jì)算公式為

(4)

式中,Va為多孔模型模擬灌漿總量；Vij為第i排、第j孔的模擬灌漿量,i=1,2,3；j=1,2,…,m。

1.3 多孔分序灌漿實(shí)現(xiàn)方法

實(shí)際灌漿工程通常采用分序加密的原則進(jìn)行灌漿,灌漿排序?yàn)橄掠闻拧嫌闻拧虚g排；灌漿孔序?yàn)棰裥蚩住蛐蚩住笮蚩住?/p>

為了真實(shí)反映分序灌漿過程,同時簡化求解流程,保證計(jì)算效率,在進(jìn)行具體的模擬計(jì)算時,同排同序孔采用相同的灌漿壓力和灌漿時間。各序孔都將經(jīng)歷3種狀態(tài),即待灌、在灌、已灌。每種狀態(tài)設(shè)置如下：①待灌孔設(shè)置為不透水邊界；②在灌孔設(shè)置為壓力入口邊界,壓力為模擬灌漿壓力；③已灌孔保持原模擬灌漿壓力邊界條件設(shè)置狀態(tài)。

裂隙巖體的灌漿過程本質(zhì)上是漿液在灌漿壓力作用下,由灌漿孔壁至出滲邊界形成穩(wěn)定滲流場的過程,滲流沿相互交叉貫通的裂隙網(wǎng)絡(luò)進(jìn)行。已有研究認(rèn)為[16-18],當(dāng)漿液擴(kuò)散的前端壓力接近地下靜水壓力時,漿液將停止流動。因此,灌漿壓力的分布在一定程度上反映了漿液分布狀態(tài)。從該角度分析,保留已灌孔的壓力場分布狀態(tài)近似等效于保留已灌孔附近漿液的賦存狀態(tài),限制相鄰在灌孔漿液的擴(kuò)散范圍及灌漿量,進(jìn)而體現(xiàn)孔序的影響。

2 耦合多孔的三維隨機(jī)裂隙巖體灌漿模型

本研究以某水電站廠①壩段帷幕灌漿工程為依托,開展耦合多孔的三維隨機(jī)裂隙巖體灌漿數(shù)值模擬的建模和數(shù)值模擬分析。根據(jù)廠①壩段地質(zhì)編錄資料,得到與該壩段實(shí)測裂隙幾何參數(shù)符合的分布類型。基巖裂隙統(tǒng)計(jì)結(jié)果見表1。

基于Monte-Carlo方法,獲得廠①壩段的三維隨機(jī)裂隙網(wǎng)絡(luò)模擬參數(shù),通過引入灌漿孔,得到耦合多孔的三維隨機(jī)裂隙灌漿模型。由于直接構(gòu)建基于整個廠①壩段的灌漿模型(最大孔深130 m)數(shù)據(jù)量過于龐大,難以進(jìn)行有限差分網(wǎng)格的離散和模擬計(jì)算，本文擇取該壩段中間的2組灌漿孔段為研究對象(埋深58～70 m,灌漿孔段12～13),介紹灌漿模擬的完整實(shí)施流程。

圖1為廠①壩段第12～13灌漿孔段耦合多孔的三維隨機(jī)裂隙巖體灌漿模型。模型尺寸為18 m×7.5 m×12 m。灌漿孔采用梅花樁形布置,以正二十棱柱表示。灌漿孔布置見圖2?；谟邢薏罘志W(wǎng)格劃分技術(shù),對灌漿模型進(jìn)行離散,得到網(wǎng)格離散模型(見圖3)。網(wǎng)格劃分尺寸為0.2 m,共劃分網(wǎng)格數(shù)1 580 337 個。

圖1 多孔隨機(jī)裂隙灌漿模型(單位：m)

圖2 灌漿孔布置(單位：m)

3 灌漿數(shù)值模擬分析

3.1 隙寬分布

通過參數(shù)賦值的方法為每條裂隙賦予特定的等效水力隙寬值,該隙寬值服從均勻分布。所得三維隨機(jī)裂隙網(wǎng)絡(luò)模型等效水力隙寬分布結(jié)果見圖4。

圖3 網(wǎng)格劃分結(jié)果

圖4 多孔灌漿模型等效水力隙寬分布(單位：m)

3.2 邊界條件

邊界條件設(shè)置見圖5。即頂面為不透水邊界,底面和四周為初始靜水壓力邊界,作為漿液出滲面,灌漿孔內(nèi)壁為壓力入口邊界。初始靜水壓力由地下水位線決定。計(jì)算得到,四周靜水壓力為0.628 MPa,底面為0.687 MPa。模擬灌漿壓力見表2。

圖5 灌漿模型邊界條件

表2 模擬灌漿壓力 MPa

3.3 漿液參數(shù)

實(shí)際灌漿工程中,本段模型所用漿液以水灰比為3∶1的單一純水泥漿為主,密度1 292 kg/m3,動力粘度2.50 mPa·s。

3.4 結(jié)果分析

3.4.1 壓力場

以下游排分序灌漿模擬過程為例進(jìn)行詳細(xì)分析,并對其他排完成模擬灌漿時的結(jié)果進(jìn)行對比分析。圖6為下游排各序孔分序灌漿數(shù)值模擬壓力場變化過程。由圖6可知,對于單一灌漿孔,壓力場整體上呈現(xiàn)由灌漿孔壁沿徑向到裂隙面內(nèi)遞減。相鄰Ⅰ序孔由于間距過大,壓力分布相互之間無影響；Ⅱ序孔與Ⅰ序孔之間壓力場存在部分重疊區(qū),說明Ⅱ序孔的擴(kuò)散已受到Ⅰ序孔的限制；Ⅲ序孔充填了Ⅰ序孔和Ⅱ序孔之間的裂隙,使上游排連接成一個整體,構(gòu)成密閉的防滲帷幕。上游排模擬灌漿壓力大小關(guān)系為：Ⅰ序孔<Ⅱ序孔<Ⅲ序孔,這反映了先灌孔對后灌孔的影響,即先灌孔封堵了后灌孔的一部分?jǐn)U散路徑,使后灌孔需要更大的壓力才能擴(kuò)散到預(yù)期的范圍,達(dá)到設(shè)計(jì)的灌漿效果。

圖6 下游排各序孔灌漿完成時壓力場變化(單位：Pa)

圖7為上游排和中間排模擬灌漿完成時的壓力場分布狀態(tài)。由圖7可知,當(dāng)上游排完成灌漿時,整體構(gòu)成了一個高壓力區(qū),說明漿液已完全充填上游排貫通的裂隙網(wǎng)絡(luò),形成防滲帷幕。同時,上、下游排的壓力場之間出現(xiàn)了重疊區(qū),說明已有漿液擴(kuò)散到了中間排區(qū)域的裂隙網(wǎng)絡(luò)內(nèi),但由于重疊區(qū)的壓力相對較小,漿液充填效果難以保證。中間排各序孔的灌漿壓力普遍大于上、下游排,這是由于上、下游排已完成了灌注,顯著封堵了漿液的有效擴(kuò)散路徑,導(dǎo)致中間排需要更高的壓力才能保證漿液對裂隙網(wǎng)絡(luò)的充填效果；當(dāng)中間排灌漿完成時,各序孔連接成一個高壓力區(qū),和上、下游排連成一體,共同構(gòu)成一道封閉的防滲帷幕。

從整體上分析,各序孔灌漿壓力基本符合以下規(guī)律：對于同序孔,下游排<上游排<中間排；對于同排孔,Ⅰ序孔<Ⅱ序孔<Ⅲ序孔。這反映了分序灌漿的特點(diǎn),即后灌孔漿液的擴(kuò)散會受到已灌孔漿液分布的阻礙,需要提高灌漿壓力,保證灌漿效果。

3.4.2 流量

圖8為上游排Ⅲ序孔灌注過程中典型單寬流量場的分布情況。由圖8可知,由于保留了已灌孔的壓力入口條件,已經(jīng)灌注結(jié)束的Ⅰ、Ⅱ序孔附近依然存在部分注入流量,但是相對于Ⅲ序孔的注入流量可以忽略不計(jì),在統(tǒng)計(jì)灌漿量時,由這部分流量產(chǎn)生的灌漿量并未做統(tǒng)計(jì)。從局部放大圖可以看出,每個灌漿孔段附近的單寬流量場分布特征與壓力場分布規(guī)律相同。灌孔的漿液擴(kuò)散受到已灌孔漿液分布的限制時,漿液的擴(kuò)散方向被迫改變,進(jìn)而使?jié){液擴(kuò)散到其他未被充填或充填程度低的裂隙網(wǎng)絡(luò)內(nèi),從而提高了防滲帷幕整體質(zhì)量。

圖9為上游排各序孔模擬漿液流量與灌漿時間的關(guān)系。由圖9可知,上游排各序孔模擬灌漿流量隨著時間的增長而逐漸下降并最終趨于穩(wěn)定,各序孔的最大灌漿流量和穩(wěn)定注入流量滿足Ⅰ序孔>Ⅱ序孔>Ⅲ序孔的關(guān)系,這說明隨著灌漿的進(jìn)行,漿液擴(kuò)散難度逐步增大,這與實(shí)際情況相符。此外,各序孔的穩(wěn)定注入流量都在1 L/min以上,這是因?yàn)閮H采用了水灰比為3∶1的單一漿液灌注,未考慮變漿過程。

圖7 分序灌漿數(shù)值模擬完成時的壓力場(單位：Pa)

圖8 典型單寬流量場分布特征(1 232.0 min)

圖9 各序孔模擬灌漿流量與灌漿時間的關(guān)系

3.4.3 灌漿量

各序孔灌漿量模擬結(jié)果與實(shí)際值對比見表3及圖10。從表3可知,該段多孔模型模擬灌漿量為40 169.85 L,實(shí)際灌漿量為45 747.50 L,相對誤差為-12.2%,說明灌漿總量的模擬結(jié)果與實(shí)際基本吻合。各序孔的模擬灌漿量與實(shí)際值基本相符,但部分存在較大誤差,如上游排Ⅰ序孔和中間排Ⅱ序孔,這體現(xiàn)了隨機(jī)裂隙網(wǎng)絡(luò)灌漿模型的不確定性。從圖10可知,Ⅰ序孔的總模擬灌漿量偏低,比Ⅰ序孔實(shí)際總灌漿量少32.86%；上游排的總模擬灌漿量偏低,比實(shí)際上游排總灌漿量少26.41%。分析可能與灌漿孔相交裂隙的數(shù)目和幾何特征等有關(guān),再次反映了隨機(jī)裂隙網(wǎng)絡(luò)灌漿模型的不確定性。

表3 各序孔模擬灌漿量與實(shí)際灌漿量結(jié)果對比

圖10 模擬灌漿量與實(shí)際值對比

3.4.4 灌漿時間

各序孔模擬單序灌漿時間為224.00 min,實(shí)際單序灌漿時間均值為221.32 min,相對誤差1.211%,模擬灌漿時間與實(shí)際情況基本相符。

4 結(jié) 語

本文對耦合多孔的三維隨機(jī)裂隙巖體灌漿數(shù)值模擬研究進(jìn)行了有益的探索,得出如下結(jié)論：

(1)由于已灌孔的漿液分布有效封堵了部分裂隙網(wǎng)絡(luò),增大了后序灌漿孔漿液擴(kuò)散的難度,各序孔的灌漿壓力基本符合分序灌漿的特點(diǎn)。同排各序孔的漿液流量基本符合Ⅰ序孔>Ⅱ序孔>Ⅲ序孔的關(guān)系。

(2)所選壩段模擬灌漿總量與實(shí)際灌漿總量的相對誤差為-12.2%,說明所構(gòu)建的灌漿模型在反映整體灌漿結(jié)果方面具有一定的可靠性。

(3)模擬單序灌漿時間和實(shí)際單序灌漿時間均值相對誤差為1.211%,驗(yàn)證了模擬灌漿時間和實(shí)際灌漿時間的一致性,保證了模擬灌漿過程與實(shí)際灌漿過程的相似性。

(4)由于隨機(jī)裂隙模型本身具有的不確定性,該模型在描述具體某一孔段灌漿量時的精度有待提高。建議繼續(xù)開展耦合單灌漿孔的三維隨機(jī)裂隙巖體灌漿數(shù)值模擬精細(xì)研究。

[1]閆福根. 水利水電工程壩基灌漿統(tǒng)一模型與分析理論及應(yīng)用[D]. 天津： 天津大學(xué), 2014．

[2]蘇培莉. 煤巖體裂隙網(wǎng)絡(luò)注漿的計(jì)算機(jī)模擬[J]. 煤炭科學(xué)技術(shù), 2011, 39(10)： 29- 33．

[3]許萬忠, 潘進(jìn)兵, 周治平, 等. 節(jié)理裂縫巖體注漿滲透模型分析[J]. 中國鐵道科學(xué), 2010, 31(3)： 47- 51．

[4]劉健, 劉人太, 張霄, 等. 水泥漿液裂隙注漿擴(kuò)散規(guī)律模型試驗(yàn)與數(shù)值模擬[J]. 巖土力學(xué)與工程學(xué)報(bào), 2012, 31(12)： 2445- 2452．

[5]OMID S, H?KAN S, SEYED R T. Numerical and analytical analyses of the effects of different joint and grout properties on the rock mass groutability[J]. Tunnelling and Underground Space Technology, 2013, 38(9)： 11- 25．

[6]XU Yong, PENG Weihong. Research on multiple holes grouting of fractured rock mass[J]. Applied Mechanics and Materials, 2012, 256- 259(12)： 547- 551．

[7]范佳俊. 高放廢物處置庫北山預(yù)選區(qū)裂隙巖體注漿特性研究[D]. 徐州： 中國礦業(yè)大學(xué), 2015．

[8]CHEN Tielin, ZHANG Liangyi, ZHANG Dingli. An FEM/VOF hybrid formulation for fracture grouting modelling[J]. Computers and Geotechnics, 2014, 58(9)： 14- 27．

[9]趙鵬. 基于漿液粘度時變特性的交叉裂隙注漿擴(kuò)散[D]. 濟(jì)南： 山東大學(xué), 2016．

[10]張慶松, 張連震, 張霄, 等. 基于漿液黏度時空變化的水平裂隙巖體注漿擴(kuò)散機(jī)制[J]. 巖石力學(xué)與工程學(xué)報(bào), 2015, 34(6)： 1198- 1210．

[11]張超, 李樹剛, 曾強(qiáng), 等. 漿液在煤巖體裂隙中滲透規(guī)律的數(shù)值模擬[J]. 煤礦安全, 2015, 46(12)： 19- 22, 26．

[12]俞文生, 李鵬, 張霄, 等. 可變傾角單裂隙動水注漿模型試驗(yàn)研究[J]. 巖土力學(xué), 2014, 35(8)： 2137- 2143, 2149．

[13]阮文軍. 注漿擴(kuò)散與漿液若干基本性能研究[J]. 巖土工程學(xué)報(bào), 2005, 27(1)： 69- 73．

[14]張有天. 巖石水力學(xué)與工程[M]. 北京： 中國水利水電出版社, 2005．

[15]倪紹虎, 何世海, 汪小剛, 等. 裂隙巖體水力學(xué)特性研究[J]. 巖石力學(xué)與工程學(xué)報(bào), 2012, 31(3)： 488- 498．

[16]劉嘉材. 裂縫灌漿擴(kuò)散半徑研究[C]∥中國水利水電科學(xué)院科學(xué)研究論文集. 北京： 水利出版社, 1982: 186- 195．

[17]朱俞仿. 裂隙傾角對灌漿擴(kuò)散半徑的影響分析[J]. 中國礦山工程, 2014, 43(2)： 55- 62．

[18]MOHAJERANI S, BAGHBANAN A, BAGHERPOUR R, et al. Grout penetration in fractured rock mass using a new developed explicit algorithm[J]. International Journal of Rock Mechanics & Mining Sciences, 2015, 80(12)： 412- 417．

(責(zé)任編輯 楊 健)

Numerical Simulation of Multi-hole Sequence Grouting of Dam Bedrock Based on Three-Dimensional Stochastic Fracture Network

LI Ruijin, ZHAO Mengqi, WANG Xiaoling, ZHU Yushan
(State Key Laboratory of Hydraulic Engineering Simulation and Safety, Tianjin University, Tianjin 300072, China)

The number of joint fractures in dam bedrock is large and its distribution is wildly random. Aiming at the research status of grouting numerical simulation in fractured rock mass, the method of sequence grouting is proposed based on three-dimensional stochastic fractured network grouting model coupled with multi-hole, and the numerical simulation of fractured rock mass grouting with the consideration of fracture aperture random distributions is carried out for practical engineering. The influences of grouting sequence on slurry diffusion pressure field and discharge per unit width are studied. The reliability of the method is verified by comparing the simulation values of grouting quantity and time with actual values．

grouting engineering; numerical simulation; three-dimensional stochastic fracture network; aperture random distribution; multi-hole sequence

2016- 11- 28

國家自然科學(xué)基金項(xiàng)目(51409186)；天津市應(yīng)用基礎(chǔ)與前沿技術(shù)研究計(jì)劃項(xiàng)目(14JCQNJC09100)

李瑞金(1990—),男,山東菏澤人,碩士研究生,主要從事水利工程模擬與優(yōu)化研究．

TV543

A

0559- 9342(2017)05- 0064- 06