胡小昕，朱澤奇

（1.中南建筑設(shè)計(jì)院股份有限公司，武漢 430000；2.中國(guó)科學(xué)院武漢巖土力學(xué)研究所巖土力學(xué)與工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，武漢 430071）

0 引 言

在地下巖土工程中，地下洞室有著非常重要的地位［1］，其在水電站、公路鐵路隧道、地下民防空間、地下核電廠等大型建筑中得到廣泛的應(yīng)用，因此地下洞室結(jié)構(gòu)安全的重要性不言而喻，而圍巖穩(wěn)定性一直是地下洞室安全性保證的重要一環(huán)。很多學(xué)者利用數(shù)值分析方法研究圍巖的變形破壞與穩(wěn)定性，但由于巖體力學(xué)參數(shù)的不確定性，近年來(lái)位移反演分析方法在巖體地下洞室工程中得到廣泛應(yīng)用［2-10］。位移反演分析是在地下工程中利用現(xiàn)場(chǎng)得到的位移推算圍巖力學(xué)性質(zhì)參數(shù)的一種方法［3］，具有良好的理論性與實(shí)用性。馮夏庭等［4］很早就將位移反演分析方法與神經(jīng)網(wǎng)絡(luò)方法結(jié)合，提出了一種用于位移反分析的進(jìn)化神經(jīng)網(wǎng)絡(luò)方法。高延法［5］利用動(dòng)態(tài)位移反演分析方法對(duì)巖移“四帶”模型進(jìn)行了擬合研究，認(rèn)為位移反演分析方法在覆巖結(jié)構(gòu)中非常有效。楊志法等［6］考慮松動(dòng)圈的影響，提出反演正算綜合預(yù)測(cè)法以提高位移反演分析方法的精度。王登剛等［7］通過(guò)將遺傳算法加入位移反演分析方法中，在考慮先驗(yàn)信息及抗噪音能力的情況下，實(shí)現(xiàn)同時(shí)反算彈性模量與泊松比。田茂霖等［8］為解決巖體非線性破壞問(wèn)題，基于Hoek-Brown 準(zhǔn)則與差分進(jìn)化算法，提出了更適用于現(xiàn)場(chǎng)的位移反演分析方法。王少杰等［9］在橫觀各向同性巖體中將差分進(jìn)化法與位移反演分析方法結(jié)合，得到了精度更高的參數(shù)。黃初濤等［10］將位移反演分析法應(yīng)用到公路隧道圍巖參數(shù)的確定中，比較了不同位移反演分析方法的優(yōu)劣。但是，基于連續(xù)介質(zhì)的位移反演分析不能夠很好地解決工程情況極其復(fù)雜的圍巖穩(wěn)定問(wèn)題，尤其是當(dāng)圍巖包含大量的裂隙時(shí)，目前常見(jiàn)的位移反演分析方法沒(méi)有較好地加以考慮，有關(guān)這方面的文獻(xiàn)報(bào)導(dǎo)也較少。

由于地質(zhì)作用的多樣性，巖體被不同尺度和方向的裂隙切割，形成空間變異性和各向異性強(qiáng)的裂隙介質(zhì)，構(gòu)成了復(fù)雜的裂隙網(wǎng)絡(luò)系統(tǒng)［11］。對(duì)于巖體裂隙網(wǎng)絡(luò)模型的應(yīng)用，目前主要將其與有限元、離散元等數(shù)值方法相結(jié)合進(jìn)行巖體裂隙滲流計(jì)算。于青春等［12］基于逆方法提出了巖體三維裂隙網(wǎng)絡(luò)的三維面狀滲流模型。劉曉麗等［13］利用RFNM 等裂隙網(wǎng)絡(luò)生成程序，編制了基于滲流力學(xué)理論的軟件GeoCAAS，較好的模擬了巖石裂隙的隨機(jī)性。相關(guān)裂隙網(wǎng)絡(luò)模型能夠?qū)α严稁r體的流體力學(xué)性質(zhì)進(jìn)行較好的模擬，但對(duì)于裂隙本身的力學(xué)性質(zhì)研究較少，一個(gè)重要的原因是目前對(duì)于裂隙結(jié)構(gòu)面的力學(xué)參數(shù)獲取還存在較多的問(wèn)題，相關(guān)理論與測(cè)試方法均需開(kāi)展深入研究。

基于目前現(xiàn)狀，本文依托某水電站地下廠房開(kāi)展裂隙網(wǎng)絡(luò)有限元反演分析，重點(diǎn)研究裂隙圍巖的建模方法與隨機(jī)模型的反演方法，得到了一些有益的認(rèn)識(shí)。

1 裂隙網(wǎng)絡(luò)模型

1.1 工程背景

某水電站地下廠房系統(tǒng)中，主廠房開(kāi)挖跨度30.80 m。主變室開(kāi)挖跨度為18.80 m。尾水調(diào)壓室凈跨度20.50～24.00 m。副廠房下部廠（橫）0+130.10～廠（橫）0+145.50 區(qū)段為III 類圍巖，在開(kāi)挖過(guò)程中出現(xiàn)裂隙發(fā)育導(dǎo)致的邊墻局部失穩(wěn)現(xiàn)象，如圖1照片所示。經(jīng)地質(zhì)調(diào)查，副廠房廠（橫）0+130.10～廠（橫）0+145.50區(qū)段主要發(fā)育四組裂隙，如表1。

圖1 副廠房上游邊墻下部局部失穩(wěn)Fig.1 Local instability of the lower part of the upstream side wall of the auxiliary power house

表1 主要發(fā)育裂隙產(chǎn)狀Tab.1 Occurrence of main developed fissures

圖2 為副廠房廠（橫）0+130.10～廠（橫）0+145.50 區(qū)段頂拱部位的地質(zhì)素描圖，在統(tǒng)計(jì)結(jié)構(gòu)面時(shí)可通過(guò)布置若干測(cè)線計(jì)算出結(jié)構(gòu)面線密度和間距。圖3基于結(jié)構(gòu)面圓盤(pán)模型給出了副廠房廠（橫）0+130.10～廠（橫）0+145.50 區(qū)段的隨機(jī)裂隙跡線圖的一個(gè)樣本。

圖2 副廠房頂拱開(kāi)挖面地質(zhì)素描圖Fig.2 Geological sketch of top arch excavation face of auxiliary plant

圖3 副廠房隨機(jī)裂隙跡線圖Fig.3 Trace map of random cracks in auxiliary power house

1.2 裂隙網(wǎng)絡(luò)模型

針對(duì)某水電站副廠房圍巖裂隙結(jié)構(gòu)統(tǒng)計(jì)結(jié)果，采用有限元軟件RS2建立接近于真實(shí)狀態(tài)的裂隙巖體模型。RS2提供了強(qiáng)大的裂隙網(wǎng)絡(luò)生成器，結(jié)構(gòu)面采用Goodman 單元模擬，其剪切行為可服從Mohr-Coulomb 滑動(dòng)準(zhǔn)則或Barton-Bandis 準(zhǔn)則。同時(shí)RS2 裂隙網(wǎng)絡(luò)生成器中還提供了多種概率分布函數(shù)，總結(jié)前人的研究，本文結(jié)構(gòu)面間距服從負(fù)指數(shù)分布，結(jié)構(gòu)面產(chǎn)狀服從Fisher 分布［14］，結(jié)構(gòu)面采用Baecher 圓盤(pán)模型，其跡長(zhǎng)服從對(duì)數(shù)正態(tài)分布［15］。

為了簡(jiǎn)化問(wèn)題，計(jì)算模型中僅考慮副廠房，不考慮主變室和尾調(diào)室的影響，建立副廠房廠（橫）0+140.00 剖面的裂隙網(wǎng)絡(luò)數(shù)值模型如圖4所示。將副廠房圍巖視為“巖石”與結(jié)構(gòu)面的組合，其屈服破壞均按Mohr-Coulomb強(qiáng)度準(zhǔn)則考慮。根據(jù)參考文獻(xiàn)［16，17］，其力學(xué)參數(shù)取值見(jiàn)表2。

圖4 副廠房裂隙網(wǎng)絡(luò)有限元計(jì)算模型（局部）Fig.4 Finite element calculation model of fracture network of auxiliary powerhouse（local）

表2 巖石與結(jié)構(gòu)面力學(xué)參數(shù)Tab.2 Parameters of rock and structural plane

裂隙網(wǎng)絡(luò)生成具有隨機(jī)性，需要用到服從不同分布形式的隨機(jī)數(shù)。以圖4 生成的裂隙網(wǎng)絡(luò)模型為考察對(duì)象，針對(duì)副廠房圍巖中較為發(fā)育的一組裂隙進(jìn)行間距統(tǒng)計(jì)，利用裂隙網(wǎng)絡(luò)生成器的裂隙信息輸出功能得到圖4計(jì)算模型中裂隙網(wǎng)絡(luò)的所有裂隙參數(shù)，圖5 給出了其中一組裂隙的間距的統(tǒng)計(jì)分布圖。結(jié)合圖5 和其他結(jié)構(gòu)面參數(shù)統(tǒng)計(jì)結(jié)果，可以認(rèn)為生成的裂隙網(wǎng)絡(luò)基本符合相關(guān)隨機(jī)性的設(shè)定，滿足表1 中a 組結(jié)構(gòu)面參數(shù)的隨機(jī)性特征。

圖5 裂隙間距的負(fù)指數(shù)分布檢驗(yàn)Fig.5 Negative exponential distribution test of fracture spacing

為了考察圖4 中裂隙網(wǎng)絡(luò)模型的工程適應(yīng)性，對(duì)副廠房開(kāi)挖進(jìn)行數(shù)值模擬，將該洞段附近安裝的多點(diǎn)位移計(jì)位移監(jiān)測(cè)結(jié)果，與該計(jì)算剖面相應(yīng)部位的位移進(jìn)行對(duì)比。圖6 給出了頂拱部位的M4-4CFX 孔口實(shí)測(cè)位移和裂隙網(wǎng)絡(luò)模型計(jì)算位移隨施工步的關(guān)系曲線。位移實(shí)測(cè)值與計(jì)算值大小在同一個(gè)量級(jí)，量值相差不大，且隨開(kāi)挖施工過(guò)程的變化趨勢(shì)也基本一致，說(shuō)明本文建立的裂隙網(wǎng)絡(luò)模型以及相應(yīng)的材料力學(xué)參數(shù)取值合理，可以作為下一步位移反演分析的基礎(chǔ)。

圖6 M4-4CFX監(jiān)測(cè)與計(jì)算位移過(guò)程線Fig.6 Displacement curves of monitoring point M4-4CFX

2 裂隙圍巖位移反演

2.1 隨機(jī)性問(wèn)題

由于結(jié)構(gòu)面參數(shù)的隨機(jī)性使得建立的裂隙網(wǎng)絡(luò)模型也具有較大隨機(jī)性。在正分析時(shí)，材料力學(xué)參數(shù)作為輸入，圍巖變形作為輸出，對(duì)于隨機(jī)裂隙網(wǎng)絡(luò)模型，則意味著一個(gè)輸入對(duì)應(yīng)理論上的無(wú)數(shù)個(gè)輸出，這一情況對(duì)于反演分析是不利的，會(huì)讓輸入?yún)?shù)的影響變得不確定。因此本文針對(duì)隨機(jī)裂隙網(wǎng)絡(luò)對(duì)圍巖變形破壞的影響進(jìn)行了研究，以表2的力學(xué)參數(shù)作為輸入，進(jìn)行20個(gè)隨機(jī)裂隙網(wǎng)絡(luò)模型開(kāi)挖試算，將計(jì)算剖面內(nèi)副廠房頂拱和邊墻位移作為考察對(duì)象，圖7 給出了副廠房頂拱位移的概率分布曲線。

圖7說(shuō)明基于隨機(jī)裂隙網(wǎng)絡(luò)模型得到的圍巖變形具有正態(tài)分布的特點(diǎn)，因此將圍巖變形的數(shù)學(xué)期望值作為輸出量是合理且符合數(shù)學(xué)規(guī)律的。

圖7 計(jì)算位移概率分布圖Fig.7 Probability distribution of calculated displacement

2.2 反演方法

結(jié)合裂隙網(wǎng)絡(luò)巖體力學(xué)特性，本文選取了對(duì)地下廠房副廠房圍巖變形影響較大的5個(gè)參數(shù)進(jìn)行參數(shù)反演。巖石變形模量E（20～30 GPa）、垂直廠房軸線方向的側(cè)壓力系數(shù)Kx（0.5～0.9）、廠房軸線方向的側(cè)壓力系數(shù)Kz（1.1～1.9）、結(jié)構(gòu)面法向剛度Kn和切向剛度Ks。

參數(shù)反演的目標(biāo)是副廠房所在計(jì)算斷面附近布設(shè)的多點(diǎn)位移計(jì)監(jiān)測(cè)位移。按照正交設(shè)計(jì)原則構(gòu)造計(jì)算樣本，對(duì)每個(gè)參數(shù)取其可能的取值區(qū)間內(nèi)的5 個(gè)水平，得到25 組計(jì)算方案。進(jìn)而采用裂隙網(wǎng)絡(luò)有限元進(jìn)行正分析計(jì)算，得到監(jiān)測(cè)點(diǎn)處的位移計(jì)算值。將得到的位移計(jì)算值與對(duì)應(yīng)的輸入?yún)?shù)組合在一起，作為神經(jīng)網(wǎng)絡(luò)位移反演分析模型的一個(gè)學(xué)習(xí)樣本，一共得到了25組樣本，用于反演分析模型的訓(xùn)練和檢測(cè)。

2.3 反演分析結(jié)果

根據(jù)本文介紹的反演分析方法，對(duì)某水電站副廠房廠（橫）0+140.00 剖面巖體力學(xué)參數(shù)進(jìn)行了反演。綜合力學(xué)參數(shù)反演結(jié)果與現(xiàn)場(chǎng)巖體力學(xué)試驗(yàn)結(jié)果，可以認(rèn)為：巖石變形模量大約在19～22 GPa 左右，垂直廠房軸線方向的側(cè)壓力系數(shù)大約為0.6～0.7 左右、廠房軸線方向的側(cè)壓力系數(shù)大約為1.6～1.8 左右，結(jié)構(gòu)面法向剛度2 800 MPa/m，切向剛度約為1 200 MPa/m。通過(guò)對(duì)比不同計(jì)算方案的圍巖變形，可以認(rèn)為裂隙結(jié)構(gòu)面法向與切向剛度對(duì)于圍巖變形的影響是較大的，本文反演方法得到的結(jié)構(gòu)面剛度是基本合理的，也是目前結(jié)構(gòu)面相關(guān)參數(shù)獲取的一種有效的方法。

針對(duì)副廠房廠（橫）0+140.00 計(jì)算剖面，用反演獲得的參數(shù)經(jīng)RS2 正向計(jì)算得出的相應(yīng)開(kāi)挖步的位移增量和監(jiān)測(cè)位移增量在量值上基本相當(dāng)，如表3 所示。其中，M42CFS 測(cè)點(diǎn)的位移實(shí)測(cè)值雖未參與反演，但位移計(jì)算值與實(shí)測(cè)值較為吻合，這進(jìn)一步表明本文建立的裂隙網(wǎng)絡(luò)模型與所采用反演分析方法在工程上是可行的，且效果較為顯著。

表3 多點(diǎn)位移計(jì)監(jiān)測(cè)結(jié)果與計(jì)算結(jié)果對(duì)比 mmTab.3 Comparison of monitoring results and analysis results

圖8為采用反演力學(xué)參數(shù)的正向計(jì)算得到的副廠房開(kāi)挖完成后裂隙破壞與圍巖變形圖?？梢钥闯?，開(kāi)挖完成后圍巖向洞內(nèi)變形，其變形形態(tài)受裂隙切割影響較大，圍巖變形最大約為30 mm。總的來(lái)說(shuō)，對(duì)于與墻面成小角度相交的裂隙，由于應(yīng)力釋放的影響不同，在裂隙面上易發(fā)生上下層面的剪切錯(cuò)動(dòng)變形；對(duì)于陡傾角結(jié)構(gòu)面，在下游邊墻表現(xiàn)為反傾切割模式，開(kāi)挖后在反傾結(jié)構(gòu)面處發(fā)生張開(kāi)變形，且與該處緩傾結(jié)構(gòu)面組合切割形成了可動(dòng)塊體，是邊墻施工穩(wěn)定的隱患；陡傾角結(jié)構(gòu)面與上游邊墻切割時(shí)則易形成多層薄層下滑體，給墻體的穩(wěn)定性帶來(lái)不利影響，在副廠房下游邊墻主要有零星掉塊的可能。

為方便比較，圖8 中用紅色表示發(fā)生破壞的裂隙，可以看出，上游邊墻多為陡傾角結(jié)構(gòu)面發(fā)生破壞，而下游邊墻則主要是緩傾裂隙發(fā)生破壞。開(kāi)挖計(jì)算結(jié)果與施工期地質(zhì)調(diào)查及圍巖穩(wěn)定評(píng)價(jià)基本相符，可以反映該局部洞段圍巖的主要變形特征與失穩(wěn)機(jī)制。

圖8 副廠房裂隙破壞與圍巖變形Fig.8 Fracture failure and surrounding rock deformation of auxiliary powerhouse

圖9 為副廠房開(kāi)挖后的最大主應(yīng)力等色區(qū)圖?？梢园l(fā)現(xiàn)，由于裂隙結(jié)構(gòu)面在開(kāi)挖卸荷過(guò)程中發(fā)生較大的剪切張開(kāi)變形，加劇圍巖應(yīng)力釋放，使得副廠房裂隙圍巖主要表現(xiàn)為應(yīng)力松弛。其中緩傾結(jié)構(gòu)面對(duì)于圍巖的應(yīng)力松弛影響較大，應(yīng)力松弛區(qū)與圍巖發(fā)生松動(dòng)變形的區(qū)域基本一致，說(shuō)明裂隙圍巖在開(kāi)挖卸荷過(guò)程中的剪切張開(kāi)變形是圍巖應(yīng)力松弛的主要原因。

圖9 副廠房最大主應(yīng)力分布Fig.9 Maximum principal stress distribution of auxiliary powerhouse

3 結(jié) 論

依托某水電站地下廠房開(kāi)展裂隙網(wǎng)絡(luò)有限元反演分析，研究裂隙圍巖建模方法與隨機(jī)模型的反演方法，得到了以下認(rèn)識(shí)：

（1）根據(jù)某水電站地下洞室施工期開(kāi)挖揭示地質(zhì)資料，以前人研究為基礎(chǔ)，建立副廠房局部范圍巖體的裂隙網(wǎng)絡(luò)模型，并通過(guò)抽樣與試算檢驗(yàn)所建立的裂隙網(wǎng)絡(luò)模型及其力學(xué)參數(shù)的合理性。

（2）考慮裂隙圍巖變形的隨機(jī)性和正態(tài)分布特征，形成基于裂隙網(wǎng)絡(luò)有限元的位移反演分析方法。基于反演力學(xué)參數(shù)的正向計(jì)算結(jié)果與施工期地質(zhì)調(diào)查及圍巖穩(wěn)定評(píng)價(jià)基本相符，可以反映裂隙圍巖主要變形特征與失穩(wěn)機(jī)制，表明本文建立的裂隙網(wǎng)絡(luò)模型與所采用反演分析方法在工程上的可行性。

（3）裂隙網(wǎng)絡(luò)有限元計(jì)算成果表明裂隙結(jié)構(gòu)面在開(kāi)挖卸荷過(guò)程中主要發(fā)生剪切與張開(kāi)變形，加劇圍巖應(yīng)力釋放，其中緩傾結(jié)構(gòu)面對(duì)于圍巖的應(yīng)力松弛影響較大。應(yīng)力松弛區(qū)與圍巖發(fā)生松動(dòng)變形的區(qū)域基本一致，說(shuō)明裂隙圍巖在開(kāi)挖卸荷過(guò)程中的剪切張開(kāi)變形是圍巖應(yīng)力松弛的主要原因。

針對(duì)裂隙網(wǎng)絡(luò)有限元的工程應(yīng)用以及裂隙結(jié)構(gòu)面力學(xué)參數(shù)反演進(jìn)行了嘗試性的研究，取得了一些初步的認(rèn)識(shí)，對(duì)于具體工程問(wèn)題的解決具有較好的參考價(jià)值。為了方便研究，本文假設(shè)所有裂隙力學(xué)性質(zhì)相同，而這與裂隙網(wǎng)絡(luò)的復(fù)雜性是不相符的，作者將在下一步的工作中結(jié)合相關(guān)反演方法重點(diǎn)研究不同裂隙組的力學(xué)參數(shù)的隨機(jī)性特點(diǎn)。