巖體灌漿的數(shù)值模擬

1、巖土工程學(xué)報Chinese Journal of Geotechnical Engineering第27卷笫8期2005年 8月Vol. 27 No. 8Aug. 2005巖體灌漿的數(shù)值模擬羅平平I,朱岳明' ,趙詠梅何山'(1.河海大學(xué) 水利水電工程學(xué)!江蘇南京210098; 2.劇0大學(xué)iilft工程學(xué)BL江西 南呂330047；3.泰安市誠信I程趣設(shè)監(jiān)理有限公司，山東泰安271000 >摘 要：認(rèn)為漿液只在連通性和滲透性好的裂隙網(wǎng)絡(luò)中流動，在模擬巖體裂隙網(wǎng)絡(luò)慕礎(chǔ)|二建立J'裂隙網(wǎng)絡(luò)賓漢姆漿 液的滲透?？?，在此模樂上對灌漿進行實時模擬.另外還研究了涎漿過程中裂

2、隙變形及對灌漿的彫響。計知結(jié)果農(nóng)明 通過該模型可預(yù)測漿液在巖體中滲透狀態(tài)，改進灌漿參數(shù).關(guān)鍵詞：裂隙巖體：灌漿：數(shù)值模擬：節(jié)理：賓漢姆漿液中圖分類號：TU472文獻標(biāo)識碼：A文章編號x 1000 - 4548(2005)080918-04作者簡介：羅平平(1976)男.江西上高人，博士主要從事巖體灌漿理論的研究。Numerical simulation of grouting in rock massLUO Pmg-pmg1,1 ZHU Yue-ming ZHAO Yong-mei HE Shan1(1. Hohai University. Nanjing 21009S. China: 2 N

3、anchang Uaivercity. Nanchang 330047 Chiu; 3. Tai'An Chenfxin ConztructionSuperintendence Ltd., Tai An 271000. China)Abstract: It is considered that giouts only flow in fractuie network with good connectivity and peneuability. Based on the simulation of fiacture network of rock mass, the seepage

4、model of the Binghanuan grouts flowing in fracture network was developed. On the basis of the model, grouting was simulated In addition, the deformation of fractures and its effects upon grouting were studied duiing the process of grouting. The results indicate that the grouting effect can be predic

5、ted and improred by the model.Key words: fracmred rock mass; grouting; numerical simulation; joint; Binghamian grouts© 1994-2008 China Academic Journal Electronic Publishing House. All rights reserved, 巖土工程學(xué)報Chinese Journal of Geotechnical Engineering第27卷笫8期2005年 8月Vol. 27 No. 8Aug. 20050引 言研究賓

6、漢姆漿液在巖體裂隙內(nèi)的滲透規(guī)律，對J- 改進灌漿工藝、優(yōu)化灌漿參數(shù)等仃匝咚意義。目前 的巖體灌漿理論只限J：單一裂隙或裂隙均勻分布的巖 體難以反映現(xiàn)場復(fù)雜的裂隙分布情況。Hassler和郝 桝對裂隙巖體灌漿進行了有益嘗試，不過他們研究 的是T頓漿液，而水灰比W/Cvi的水泥漿屬安漢姆 流體叫 且他們沒仃右慮裂隙變形對灌漿的影響。楊 米加建J'裂隙巖體灌漿模型，研究J'賓漢姆漿液 的滲透規(guī)律，但他只建立了等寬裂隙網(wǎng)絡(luò)，II其模型 中的流速為冗漢姆漿液在水平裂隙中的平均流速,但 巖體中的裂隙分布錯綜復(fù)雜。本文將在梵模型上進行 有益改進。1巖體灌漿裂隙網(wǎng)絡(luò)模型人然巖體裂隙分布狀況及兀

7、復(fù)雜，13前難以確定 裂隙精確位置和產(chǎn)狀。忖前主耍采用Monte-Carlo法 來模擬巖體裂隙網(wǎng)絡(luò).Samaiueg151 > Bmgham61 Shanley1 Priest和Hudson18 9.潘別腳等對此有深 入研究.Monte-Carlo法是先用實測統(tǒng)計建立的裂隙幾 何參數(shù)(包括巖體裂隙組數(shù)，每組的產(chǎn)狀、仙距、跡 長)概率模熨，然厲利用隨機數(shù)產(chǎn)生服從這種概率模 型的隨機變最，最終生成巖體裂隙網(wǎng)絡(luò)圖。圖1巖體裂隙網(wǎng)絡(luò)圖Fig. 1 Fracture net of rock inass衣1是某壩址區(qū)巖體裂隙幾何參數(shù)統(tǒng)計，裂隙網(wǎng)牧稿日期：2001- 1008© 1994-2

8、008 China Academic Journal Electronic Publishing House. All rights reserved, 第8期羅平平.等巖體灌漿的數(shù)值模擬921裂做 分組傾角（止態(tài)m（負(fù)指數(shù)m斷沖（均勻）An跡長（負(fù)指數(shù)）m隙如D均值方差均伉方差均值方差均值方差T1709.00806191.52.50.50 001021207.21.00.43.01.42.4080 00123109.30.720.32.52.01.80.40.000S絡(luò)模擬結(jié)果見圖1 （圖中粗線段肋為灌漿孔/M9Jsl表1某壩址區(qū)巖體裂隙幾何參數(shù)統(tǒng)計Table 1 Statistics of

9、 geometric paiameteiof rock on dam site2裂隙巖體灌漿數(shù)值模型首先假定：漿液在裂隙中為層流；巖體為完 全彈性；漿液為賓漢姆漿液；裂隙是光滑的，11 裂隙在灌漿壓力卜不擴展延伸。二維巖體裂隙網(wǎng)絡(luò)苴構(gòu)成耍素主要足裂隙交點及 交點間的裂隙段，根據(jù)這一兒何特征，町把裂隙交點 作為節(jié)點，"點間的裂隊段作為線單尤。在此裂隙網(wǎng) 絡(luò)中取-由節(jié)點i和該節(jié)點相交的M'條線單元組成 的均衡域，根據(jù)質(zhì)鼠守恒方程，可得節(jié)點i處的漿液 流動方程（非穩(wěn)定流）：（M' 曲dr工乞-工® +0=-么Vs1 ）i Vs1 Ji（i = l,2,N）,其屮竹

10、為與節(jié)點i相連的/線元內(nèi)的流眉；M'為與 節(jié)點j相連的線元數(shù)；©為i點的垂苴補給彊；Qj為S M匯源項；h為節(jié)點i上的灌漿壓.力;d嚴(yán)十工勺1廠2 7-15.為裂隙以i點為中心的表征單元域內(nèi)彈性儲（釋） 水系數(shù)；勺和-分別為與，點柑連的_/線尤寬度和長 度；/為灌漿時間。賓漢姆漿液流經(jīng)單-傾斜裂隙的平均流速為u =竺-（/ ± sin0） +, r° -f,12v3vegp （I ± sm 0） 4pv 其屮&為裂隙傾角，當(dāng)漿液從裂隙上部滲流到卜部時 取正值，否則取負(fù)值。若忽略裂隙傾斜時的貞力影響，則n為國ge2°M=/+,12

11、v 3wgp rpv此時流經(jīng)該裂隙的流鼠為-g，,可roe2q = ue = -/ + 12v3vg-pI4pvtt屮g為車力加速度：。対裂隙寬度："為漿液運動 粘度；/為與i點相連的某銜接線尤的滲透梯度；q為 漿液剪切強度；Q為漿液密度。(4)把式（3）代入（1）式可得_ 工®j +Q = _4普 /=1 Jiu,(心1,2,N), Z1Zo(5)式屮 0 =w'12v； 4pvl）若區(qū)域內(nèi)有N個節(jié)點，M個線尤，則式（4）可 用張量表示為c 小s < f、dhaj,j - aucoi + Q = 一扌 G/J 二(i = 1,2,N; J = 1,2,M)

12、,(6)i為自由標(biāo)，J為啞標(biāo)；A = «.vxw為裂隙網(wǎng)式屮絡(luò)的NxM階銜接矩陣，該矩陣屮尤素®可表述為j線元銜接于i節(jié)點.II指向i節(jié)點方向,0.J線元銜接丁i節(jié)點,-Lj線元銜接“節(jié)點, 且指向i節(jié)點方向o灌漿床力導(dǎo)致裂隙的擴張和巖體的變形.裂隙的 張開引起注漿磺的增丿川C巖體裂隙邊界受灌漿床力的 作用引起的變形可用半無限邊界上作用仃均布力的彈 性變形Boussinesq公式求解（如圖2）,求得裂隙變形 最為=)ap厶(i_in厶),7rE= 2(1/r)APL(l - In L),7rE分別為裂隙中點和端點由灌漿床力P引起的變形最；E、"分別為巖體彈模和泊松

13、系 數(shù)：厶為裂隙長度。式（6）, （7）構(gòu)成了裂隙巖體灌漿控制方程。根據(jù)以上結(jié)論，作者編制了相應(yīng)程序，對賓漢姆 漿液（主耍是水泥漿）在裂隙網(wǎng)絡(luò)屮的滲透規(guī)律進行 了研究.3漿液在裂隙網(wǎng)絡(luò)中的滲流特性分析某壩址區(qū)巖體裂隙網(wǎng)絡(luò)和計算域見圖1.其中頂部居中位置一深10 m灌漿孔，假定左右兩側(cè)節(jié)點為 出漿點，壓力為0NlPa°主耍灌漿參數(shù)為：整個灌漿 孔段灌漿壓力為2 MPa ,漿液密度 p = 1.5xlO3kg/m5 ,動丿J粘度 =35 MPa s ,儲水 系數(shù)5,-03,漿液初始剪切r0-14 Pa.整個灌漿時 間為600 So圖2灌漿壓力引起裂隙面變形Fig 2 The deform

14、ation of fissures caused by grouting pressure 計算分兩種丁況：第一種匸況是將整個灌漿時間 600 s平均分成10個時步，不考股裂隙面的變形；第 二種工況也是將整個灌漿時間600 s平均分成10個時 步，但考慮裂隙變形和灌漿的相耳影響。限J：篇幅， 只給出這兩個工況第一和址后一個時步分別不考渥和 考電裂隙變形情況卜，節(jié)點灌漿壓力等值線圖（分別 見圖36）。圖3不考慮裂隙面變形/-60s節(jié)點灌漿壓力等值線(MPa)Fig 3 Contour of the nodal pressure when r 60 s and the deformation of

15、 fracture is ignored (MPa)圖4不考慮裂隙面變形r-600s節(jié)點灌漿壓力等值線(MPa)Fig. 4 Contoui of the nodal piessuie when t 600 s and the deformation of fracture is ignoied (MPa)圖6考慮裂隙面變形600 s節(jié)點灌漿壓力等值線（MPa）Fig. 6 Contour of the nodal pressure when f 60 s and thedeformation of fracture is considered （MPa）從圖4, 6可以對比看出,若考慮由丁灌

16、漿壓力而 導(dǎo)致裂隙而變形，漿液滲透速度加快，摻透范用也更 廣。為直觀反映裂隙面變形對漿液流動影響，以"24 和“110節(jié)點為例（節(jié)點位置見圖1標(biāo)識）,筆者作了 芳電裂隙面變形前后灌漿壓力隨灌漿時間變化的對比 圖（圖7）。從圖7可以看出：離灌漿孔越近的節(jié)點， 節(jié)點斥力越人；考慮裂隙面變形時的節(jié)點灌漿斥力要 人不考偲裂隙面變形的節(jié)點灌漿爪力：巧必灌漿壓 力逐漸趨穩(wěn)定。浹289節(jié)點足孤點，不論是否考農(nóng)裂 隙變形，梵圧力總是0 MPa。這些都與實際巖體灌從圖3, 5可以看出,灌漿60s時，漿液只分布 在灌漿孔兩側(cè)附近。但是隨看灌漿的進行，漿液便沿 若裂隙網(wǎng)絡(luò)交叉滲透，逐漸深入。裂隙隨機分布、

17、相 互交義的這種網(wǎng)絡(luò)特征極人影響了漿液滲透特性，從 圖4, 6可以看出，盡管已經(jīng)灌漿了 600 s,但局部區(qū) 域節(jié)點灌漿壓力還是0 MPa。這由這些節(jié)點所處 裂隙都足孤芒的，不ri接與灌漿孔相連，或沒仃間接 通過其它裂隙與灌漿孔相連，不是漿液滲透路徑。圖5考慮裂隙面變形/-60s節(jié)點灌漿壓力等值線(lPa)Fig. 5 Contoui of the nodal pressure when f 60 s and thedeformation of fracture is considered (MPa)漿特征相符合。-*24 I.況w- »24 |.況二 -*-1101.況一 +011

18、OT 況二圖7 “24和“110的節(jié)點壓力時程對比Fig. 7 Comparison between the nodal pressures of Node 24 and Node 110從圖8可以看出若考世裂隙變形，灌漿初始階段 灌漿就也越多，但衰減快，意味著將更快結(jié)束灌漿， 另外還可看出灌漿吊隨時間是遞減的，根據(jù)此特征可 估計出，若不考慮裂隙變形和灌漿的相互影響，完全 灌漿時間大致要900 s,否則灌漿時間只需720 s。計 算表明灌漿到600 s時，不考慮裂隙變形影響時，總 灌漿駅為76.21,否則總灌漿最為791。號慮了衆(zhòng)隙變形° 60 IM) 30()420 M0 600圖

19、8計算時步的濯漿呈隨時間變化Fig.8 Vaiiation of grout quantity with time4結(jié)語（1）該裂隙巖體灌漿數(shù)值模型較好，能實時反映 灌漿進程以及漿液在裂隙網(wǎng)絡(luò)屮的流動規(guī)律。不足Z 處在r模型只考慮了漿液的單相流，事實上流體是以 兩相（漿液和水）艮至三相（漿液、水和氣）的形式 在巖體裂隙內(nèi)流動。另外模型沒令考慮裂隙粗糙度對 漿液流動的影響，這些都是今后的研究方向。(2) 巖體裂隙隨機分布、相互交叉的網(wǎng)絡(luò)特征， 影響著漿液滲透特性。若不考虎裂隙擴展延伸，漿液 總是沿著相互連接并與灌漿孔相連的裂隙流動。(3) 該算例中計算時步的灌漿屆隨灌漿時間是遞 減的，若考慮裂隙

20、變形與灌漿相互影響，則衰減速度 更快，根據(jù)此特征可估計出結(jié)束灌漿的時間。(4) 節(jié)點灌漿壓力隨時間遞增,并逐步趨J：穩(wěn)定。參考文獻：1 Hassler Lars, Hakansson Ulf. Hakan Stille Computer- simulated flow of grouts in jointed rockJ. Tunneling and Undergiound Space Tech. 1992, 7(4):441 - 446.2 郝 竹,王介強，何修仁.巖體裂隙注漿的計克機模擬J%七工程學(xué)報.1999. 21(6):727 - 7303 Ski Caijun, Day Robert L. Acceleration of the leactivity of flyash by chemical activationJ. Cement and Concrete Reseaich. 1995.25 15 -21.