張?zhí)m新，賈超

(1.山東省地礦工程勘察院，山東 濟(jì)南 250014；2.山東大學(xué)海洋研究院，山東 青島 266237)

0 引言

巖體被大量結(jié)構(gòu)面切割成各種大小不一的塊，當(dāng)塊體空間狀態(tài)滿足某些條件時，就容易滑落而失穩(wěn)，形成關(guān)鍵塊體[1]。通常情況下，巖體強(qiáng)度遠(yuǎn)大于結(jié)構(gòu)面強(qiáng)度，20世紀(jì)80年代中期，石根華和R.E.Goodman提出了一種分析非連續(xù)塊體穩(wěn)定性的理論，即關(guān)鍵塊體理論[2～4]。將巖塊視為剛體，進(jìn)行力學(xué)和幾何學(xué)分析，并依據(jù)“有限性定理”和“可動性定理”分析塊體穩(wěn)定性[5～8]。該方法計算簡單高效，成果可靠，已被工程界高度認(rèn)可，普遍應(yīng)用于各種地下洞室和隧道開挖工程的圍巖穩(wěn)定分析。之后，國內(nèi)外學(xué)者對其進(jìn)行了發(fā)展和完善，應(yīng)用范圍廣泛。然而，傳統(tǒng)的關(guān)鍵塊體分析方法大都假設(shè)塊體只受重力的作用，而忽略了內(nèi)部塊體對臨空面塊體的作用力[6]?？紤]巖體漸進(jìn)破壞及塊體間力傳遞的分析方法在國內(nèi)外尚不多見，此外，被錨固在一起的塊體可能形成新的不穩(wěn)定的整體，因此還要對錨固后的塊體進(jìn)行二次評估。

本文應(yīng)用基于力傳遞算法的改進(jìn)的關(guān)鍵塊體理論并結(jié)合矢量分析方法，對泰山抽水蓄能電站地下洞室隨機(jī)塊體進(jìn)行了穩(wěn)定性分析和預(yù)測。對比分析了力傳遞前后塊體的穩(wěn)定性，對已錨固的塊體進(jìn)行了二次評估，驗算了廠房原有錨桿支護(hù)系統(tǒng)的支護(hù)效果，推薦不穩(wěn)定塊體的支護(hù)措施和加固方案。

1 計算理論與方法

基于力傳遞算法的關(guān)鍵塊體分析方法是在關(guān)鍵塊體理論的基礎(chǔ)上針對堅硬或半堅硬地下洞室圍巖穩(wěn)定性評估，將塊體分層，考慮巖體漸進(jìn)破壞應(yīng)用力傳遞算法計算不同批次塊體對臨空面塊體的作用力，分析臨空面塊體的穩(wěn)定性，最終確定合理的圍巖支護(hù)方案。在力傳遞分析之后，可以把之前視為固定的相鄰塊移除而單獨(dú)考慮要研究的塊體。考慮塊體失穩(wěn)的自由滑落、單平面滑動和雙平面滑動3種形式。只需分析計算臨空面上塊體的穩(wěn)定即可，因為只要臨空面上塊體安全穩(wěn)定就能保證地下廠房的穩(wěn)定。

1.1 單塊體穩(wěn)定性分析方法

巖體中的結(jié)構(gòu)面的走向、傾向和傾角等產(chǎn)狀發(fā)育具隨機(jī)的特點(diǎn)，因此產(chǎn)生的隨機(jī)塊體也具有隨機(jī)的特點(diǎn)[7～11]。在地質(zhì)實測資料的基礎(chǔ)上，統(tǒng)計分析研究區(qū)節(jié)理裂隙產(chǎn)狀信息，模擬生成結(jié)構(gòu)面網(wǎng)絡(luò)，搜索關(guān)鍵塊體，計算分析其穩(wěn)定性，是傳統(tǒng)關(guān)鍵塊體理論分析的一般步驟[12～14]。單塊體穩(wěn)定性分析是后續(xù)分析的基礎(chǔ)。錯綜復(fù)雜的節(jié)理裂隙將邊坡或隧道切割如圖1所示，形成的塊體如圖2所示。

圖1 被結(jié)構(gòu)面切割的巖體

圖2 塊體滑動矢量分析示意圖

三種不同的失效模式下，對塊體滑動進(jìn)行向量分析，可移動性及滑動方向的確定方法如下:首先，定義點(diǎn)積

(1)若?i/Ti=-1，則塊體為穩(wěn)定塊體；

(2)若?i/Ti＜0，則為自由滑落塊，此時滑動方向sd=r；

(3)若 ?i/-1＜Ti=0且 ?k/Tk≥0,k≠i，則為單平面滑落，此時滑動方向sd=(n1×r)×ni或sd=r-(r?n1)n1；

(4)若

?i/-1＜Ti=0&?j/-1＜Ti=0,

?k/Tk≥0,k≠i&k≠j，則為雙平面滑動，此時滑動方向為 sd=sign[(n1× n2)?r)](n1× n2)。

上式中r為驅(qū)動力向量R的單位向量；ni為塊體上第i個切割面的單位法向量；sd為失穩(wěn)塊體的滑動方向；ni為第i 塊塊體節(jié)理裂隙面上的單位法向量[15～17]。需要注意的是力傳遞后的穩(wěn)定分析中，塊體上最終的驅(qū)動力已經(jīng)包括了后面批次塊體對臨空面批次關(guān)鍵塊體的作用力。

1.2 塊體間力傳遞過程

1.2.1 力傳遞算法

在傳統(tǒng)的關(guān)鍵塊體分析中，研究某特定塊體時假定與其相鄰的所有塊體都固定不動(Warburton，1987)。事實上，最終塊體上的驅(qū)動力應(yīng)該包含內(nèi)部塊體間的作用力。但是，目前的關(guān)鍵塊體分析方法并沒有考慮這樣不同批次關(guān)鍵塊體之間力的傳遞，這有可能導(dǎo)致我們對臨空面批次關(guān)鍵塊體上滑動力估計不足，由此得出的支護(hù)方案有可能不能保證其安全穩(wěn)定。為了考慮不同批次關(guān)鍵塊體之間的影響，引入一種力傳遞算法。對于圍巖支護(hù)設(shè)計而言只有滑移面指向洞室內(nèi)部的塊體才是要研究的關(guān)鍵塊體。如果這些關(guān)鍵塊體被支護(hù)固定在其原來的位置而不滑動就能保證地下洞室整體的安全穩(wěn)定性。內(nèi)部塊體的作用力逐層傳遞到臨空面批次關(guān)鍵塊體，由此計算的支護(hù)力將會更加準(zhǔn)確。

在力的傳遞過程中，塊體的內(nèi)部作用可以分為穩(wěn)定影響和不穩(wěn)定影響兩類，如圖3所示。

圖3 塊體間的兩類作用

力傳遞分析后，可以把之前視為固定的相鄰塊移除而單獨(dú)考慮要研究的塊體。塊體間可傳遞力的接觸面定義為有效接觸面，后面批次關(guān)鍵塊體產(chǎn)生的不平衡力，通過有效接觸面?zhèn)鬟f到前面批次關(guān)鍵塊體上。力傳遞算法的基本步驟如下:

(1)判定接觸面是否為有效接觸面。關(guān)鍵塊體之間的接觸面是否為有效接觸面取決于塊體的滑動方向向量n(n1,n2,n3)與接觸面的法向量m(m1,m2,m3)的相對關(guān)系。當(dāng)n?m＞0時，接觸面即為能傳遞力的有效接觸面；

(2)計算有效接觸面?zhèn)€數(shù)及面積。對于每一個位于后面第i+1批次的關(guān)鍵塊體Bm和與其相鄰的位于前面第i批次的關(guān)鍵塊體Bn，統(tǒng)計并識別它們二者之間的接觸面的個數(shù)。若某個塊體產(chǎn)生的滑動力傳遞到多個塊體，有效接觸面的面積為接觸面在滑動方向法線方向的投影，且傳遞的滑動力的大小與投影面積成正比；

(3)當(dāng)滑動力由后面第i+1批次塊體傳遞到前面第i批次關(guān)鍵塊體后，應(yīng)利用基于向量分析法的關(guān)鍵塊體理論重新評定前面第i批次關(guān)鍵塊體的安全穩(wěn)定性；

(4)重復(fù)過程(1)～(3)，直到后面所有批次關(guān)鍵塊體上的不平衡滑動力全部傳遞到第一批次也就是臨空面批次關(guān)鍵塊體。至此，我們根據(jù)臨空面關(guān)鍵塊體上這個最終滑動力來確定支護(hù)力的大小及方式[12]。

假設(shè)后面批次關(guān)鍵塊體在力傳遞前后沒有位移，據(jù)此力傳遞模型傳遞力可通過式(2)計算得到:

式中，F(xiàn)i為傳遞到第i塊關(guān)鍵塊體上的不平衡力；F為前面批次關(guān)鍵塊體產(chǎn)生的總不平衡力；Si為塊體間有效接觸面積為前面批次塊體與后面批次所有塊體接觸的總面積。通過漸近破壞分析方法，內(nèi)部所有塊體的不平衡力將會逐次傳遞到臨空面塊體。

1.2.2 穩(wěn)定評價

傳統(tǒng)的關(guān)鍵塊體理論在進(jìn)行塊體穩(wěn)定性評價計算安全系數(shù)時，只考慮了塊體的自重而忽略了后面批次關(guān)鍵塊體對前面批次關(guān)鍵塊體的作用力，這樣計算得到的塊體安全系數(shù)有可能與實際情況偏差過大。本文基于力傳遞算法的關(guān)鍵塊體理論穩(wěn)定評價，是在傳統(tǒng)安全系數(shù)計算的基礎(chǔ)上考慮了后面批次關(guān)鍵塊體的不平衡力后計算的到的，更加符合實際情況。安全系數(shù)計算的概化公式如下:

已被錨桿錨固的關(guān)鍵塊安全系數(shù)定義如下:

(b)單平面滑動:

(c)雙平面滑動:

式(3)～(5)中，n 為錨桿數(shù)量；G 為塊體自重；GNi和Gsi分別為塊體自重在滑動面i上垂向和切向分量；Gsij為平面i和j交線方向上塊體重力的分量；RBg為錨桿錨固力在重力方向上的分量；RBni和RBsi分別為錨桿錨固力在平面i垂向和切向方向上的分量；RBsij為錨桿錨固力在平面i和j交線方向上的分量；φj，ci和Ai分別是摩擦角、粘聚力及平面i的面積[18～22]。

2 工程算例分析

2.1 工程概況

泰山抽水蓄能電站樞紐由上水庫、下水庫、輸水系統(tǒng)、地下廠房洞室群及地面開關(guān)站等組成，裝機(jī)容量4×250 MW。地下廠房及輸水系統(tǒng)布置在上水庫右岸橫嶺南坡呈北東向展布的山梁內(nèi)，上覆巖體約210～240 m，平均埋深為224.7 m，頂拱高程為140.275 m，底板高程為86.60 m。主廠房尺寸為180.0×24.5×52.3 m(長×寬×高)，走向為 N40°W。地下洞室群巖層主要為交代式花崗巖、斑紋狀混合巖、混合花崗巖、黑云斜長片麻巖夾斜長角閃巖等，巖層主要產(chǎn)狀為 N60～80°E SE(NW)∠70～85°。本區(qū)層理發(fā)育，為ⅠⅠⅠ類巖石，巖體飽和抗拉強(qiáng)度為0.6～0.7 Mpa。工程區(qū)最大水平主應(yīng)力方向為N60～80°E，水平主應(yīng)力介于7～12 Mpa。通過廠房的斷層有 f9、f17、f18、f25、f37、fc1 等 6 條，f25 斷層規(guī)模稍大，其余均較小。NNE及NEE向陡傾角節(jié)理較發(fā)育?？紤]到地下主廠房的穩(wěn)定問題，尤其頂拱隨機(jī)塊體的安全穩(wěn)定尤其重要，因此利用改進(jìn)的關(guān)鍵塊體理論對地下主廠房隨機(jī)塊體進(jìn)行穩(wěn)定性計算，對錨固的關(guān)鍵塊體進(jìn)行二次評估，復(fù)核錨桿支護(hù)方案，保證地下廠房安全穩(wěn)定運(yùn)行具有重要意義。

2.2 參數(shù)選取及模型建立

對工程區(qū)巖體節(jié)理裂隙進(jìn)行統(tǒng)計分析，重點(diǎn)研究對廠房安全穩(wěn)定影響較大及發(fā)育程度高的結(jié)構(gòu)面，統(tǒng)計結(jié)果見表1，根據(jù)實際情況，本文對表1中的結(jié)構(gòu)面按(1，3，4)，(1，3，6)，(2，3，4)，(3，5，7)，(1，5，6)，(1，6，7)，(2，3，6)，(2，5，6)，(1，2，6)進(jìn)行面組合。

繪制節(jié)理裂隙極點(diǎn)散點(diǎn)及等值線圖，并結(jié)合廠房軸線走向繪制極點(diǎn)等值線見圖4。

對傾向分布進(jìn)行擬合，結(jié)果如圖5所示。

根據(jù)工程區(qū)節(jié)理和巖石精細(xì)化采集描述的結(jié)果及地下洞室設(shè)計、地質(zhì)勘查實測資料，確定隨機(jī)塊體穩(wěn)定分析相關(guān)物理力學(xué)參數(shù)如表2。

建立地下廠房三維數(shù)值模型，廠房橫面尺寸為180m×24.5m×52.3m，模型尺寸 180m×125m×150m，錯綜復(fù)雜的結(jié)構(gòu)面將水電站地下洞室切割，如圖6所示。

表1 廠房區(qū)巖體結(jié)構(gòu)面產(chǎn)狀統(tǒng)計成果

表2 地下廠房穩(wěn)定分析參數(shù)確定

圖4 節(jié)理極點(diǎn)散點(diǎn)分布統(tǒng)計圖

圖5 優(yōu)勢節(jié)理傾向分布擬合圖

圖6 地下廠房三維模型

2.3 隨機(jī)塊體穩(wěn)定計算及結(jié)果分析

通過對該水電站地下廠房地質(zhì)條件的統(tǒng)計分析知，頂拱及上下游邊墻處存在塊體滑動的可能較大，頂拱處形成的各層次塊體二維圖形見圖7。

圖7 頂拱各層塊體二維圖形

經(jīng)節(jié)理裂隙切割后塊體的典型空間形式，典型塊體主要有4種，見圖8。

圖8 塊體形態(tài)

只要臨空面關(guān)鍵塊體保持穩(wěn)定就能保證洞室穩(wěn)定，考慮力傳遞計算臨空面關(guān)鍵塊體的安全系數(shù)來評定洞室穩(wěn)定性，逐一分析計算各不利組合形成的臨空面塊體穩(wěn)定性，見圖9。

圖9 臨空面塊體搜索與穩(wěn)定計算

為保守起見，切割面按貫通率為百分之百計算。分別計算力傳遞前后塊體驅(qū)動力結(jié)果，見圖10。

對結(jié)果進(jìn)行分析，頂拱及上下游邊墻塊體驅(qū)動力并不大一般在180 kN范圍內(nèi)。統(tǒng)計知，力傳遞后塊體總驅(qū)動力普遍增大，只有4塊有降低趨勢，編號為 3592、7564、1876 和 2364，3 塊在側(cè)壁，1 塊在頂拱，占總塊體的15%。一塊不變，占總塊體的3%。分析原因，是由于內(nèi)部塊體傳遞來的力是重力形成的，主要為向下的分量，而向下的力對頂拱處塊體最為不利，側(cè)壁處塊體，有不利的下滑力和抗滑力兩個影響。工程區(qū)節(jié)理裂隙發(fā)育，巖體被切割成的塊體體積較小，對某個特定的塊體而言抑制其下滑的塊體少，且按結(jié)構(gòu)面完全貫通巖體考慮。因此需要的支護(hù)力也會相應(yīng)的改變。

為確定塊體最終所需的支護(hù)力強(qiáng)度，確定合適的錨桿支護(hù)系統(tǒng)，須進(jìn)一步計算塊體的安全系數(shù)，計算結(jié)果見圖11。

圖10 力傳遞前后驅(qū)動力計算結(jié)果

圖11 力傳遞前后安全系數(shù)計算結(jié)果

對于力傳遞前后的塊體安全系數(shù)可知，由于力傳遞的效應(yīng)造成的塊體安全性的改變是比較明顯的。由塊體安全系數(shù)計算結(jié)果知，僅3塊關(guān)鍵塊體的安全系數(shù)有所增加，占11%，其他塊體的安全系數(shù)普遍下降，而且下降幅度各不相同，4塊塊體由穩(wěn)定塊變成了不穩(wěn)定塊，2塊由原來的不穩(wěn)定變成穩(wěn)定塊。對比圖11和圖12知，驅(qū)動力增大的塊體，安全系數(shù)不一定降低，而且安全系數(shù)降低的幅度與驅(qū)動力改變的大小也無類似關(guān)系，這是由于力傳遞后在臨空面上塊體的總驅(qū)動力變化的同時總抗滑力也在改變，二者的變化無內(nèi)在聯(lián)系，它是結(jié)構(gòu)面大小、產(chǎn)狀、貫穿度等多種因素共同決定的。經(jīng)進(jìn)一步計算，對安全系數(shù)小于1.5的塊體進(jìn)行加錨支護(hù)，推薦采用如下的錨桿支護(hù)系統(tǒng)f25@1.5×1.5 m，L=8.4 m，錨桿交錯間隔布置形式，支護(hù)力為58.24 kN/m2，結(jié)構(gòu)面摩擦系數(shù)取0.54，錨固后二次評估計算得到所有關(guān)鍵塊安全系數(shù)均大于1.5，滿足安全要求，因此原有的錨桿支護(hù)系統(tǒng)由于其原有安全余量，力傳遞計算后仍能滿足塊體穩(wěn)定的要求，說明錨桿在該洞室安全性上發(fā)揮了重要作用。

3 結(jié)語

本文針對巖體地下洞室穩(wěn)定性問題擴(kuò)展了傳統(tǒng)的關(guān)鍵塊體理論分析方法，充分考慮了巖體內(nèi)部結(jié)構(gòu)面空間分布的不確定性。采用基于力傳遞算法的擴(kuò)展的關(guān)鍵塊體理論對洞室圍巖穩(wěn)定性進(jìn)行分析，該分析方法更為科學(xué)合理。通過本文研究可以得到以下結(jié)論:

在巖質(zhì)地下洞室開挖工程中，隨機(jī)塊體直接滑動、單平面滑動和雙平面滑動是常見的破壞模式，對于產(chǎn)狀復(fù)雜的各向異性結(jié)構(gòu)面，結(jié)合力小，節(jié)理裂隙發(fā)育，進(jìn)行隨機(jī)塊體穩(wěn)定分析十分必要。

考慮洞室內(nèi)部關(guān)鍵塊體的先后作用，采用漸進(jìn)破壞分析方法，應(yīng)用力傳遞算法，研究關(guān)鍵塊體不連續(xù)巖體力學(xué)特性，并對相應(yīng)的塊體進(jìn)行穩(wěn)定分析，是評價地下洞室穩(wěn)定性的有效方法之一。

泰山抽水蓄能電站地下洞室塊體穩(wěn)定性的改變是各種因素共同作用的結(jié)果，具有一定的隨機(jī)性，對隨機(jī)塊體穩(wěn)定性及錨桿系統(tǒng)計算復(fù)核表明塊體安全系數(shù)隨支護(hù)過程發(fā)生改變，選定的系統(tǒng)錨桿支護(hù)方案能夠滿足地下洞室塊體穩(wěn)定性要求。

本文分析是基于一個確定性的巖體開挖模型計算的。考慮實際上難以準(zhǔn)確測量和統(tǒng)計出全部不連續(xù)節(jié)理裂隙的空間分布特征，及完整巖體和節(jié)理裂隙的力學(xué)參數(shù)，所以將來采用概率方法來充分考慮這些節(jié)理裂隙巖體各參數(shù)及其空間分布特征的不確定性將更具科學(xué)意義。