周 亮，孫茂前

(國(guó)核電力規(guī)劃設(shè)計(jì)研究院，北京 100095)

1 工程概述

某核電站規(guī)劃分三期建成，其中三期工程進(jìn)行廠區(qū)場(chǎng)地開(kāi)挖后，在場(chǎng)地南部燈塔山與廠區(qū)連接處的燈塔嶺東北側(cè)形成140余米長(zhǎng)的人工邊坡。其距擬建核島廠房最近距離190m，邊坡走向137°，坡腳開(kāi)挖起自廠坪標(biāo)高14.5m，坡高11m～20m。

2 巖體結(jié)構(gòu)面網(wǎng)絡(luò)模擬

在現(xiàn)場(chǎng)勘測(cè)大量結(jié)構(gòu)面調(diào)查統(tǒng)計(jì)資料的基礎(chǔ)上，應(yīng)用巖體結(jié)構(gòu)面網(wǎng)絡(luò)模擬理論、方法和程序，對(duì)該核電站人工邊坡進(jìn)行結(jié)構(gòu)面幾何形態(tài)要素的統(tǒng)計(jì)分析和巖體結(jié)構(gòu)面網(wǎng)絡(luò)模擬的應(yīng)用。為了能夠全面反映人工邊坡巖體結(jié)構(gòu)面的分布和組合狀況，采用測(cè)線量測(cè)法，在測(cè)區(qū)內(nèi)選擇有代表性、結(jié)構(gòu)面出露條件較好的典型地段進(jìn)行了結(jié)構(gòu)面系統(tǒng)量測(cè)。根據(jù)統(tǒng)計(jì)的結(jié)構(gòu)面樣本179 條進(jìn)行結(jié)構(gòu)面分組，其中圖1為傳統(tǒng)極點(diǎn)投影圖，可以看出測(cè)區(qū)內(nèi)巖體結(jié)構(gòu)面發(fā)育較密集的有3組，傾向分別為NEE、SE和NW，表1為現(xiàn)場(chǎng)量測(cè)的結(jié)構(gòu)面產(chǎn)狀(同時(shí)分解成三組形式表示結(jié)構(gòu)面的分組結(jié)果)；而采用動(dòng)態(tài)聚類法分析結(jié)構(gòu)面分組時(shí)，大致以此3組作為初始聚類參考，圖2為動(dòng)態(tài)聚類結(jié)構(gòu)面分組，優(yōu)勢(shì)產(chǎn)狀見(jiàn)表2。

表1 結(jié)構(gòu)面產(chǎn)狀 單位：°

圖1 極點(diǎn)投影圖

圖2 動(dòng)態(tài)聚類結(jié)構(gòu)面分組

表2 優(yōu)勢(shì)產(chǎn)狀

結(jié)構(gòu)面分組完成后，在分析各組結(jié)構(gòu)面幾何形態(tài)要素概率分布形式的基礎(chǔ)上，建立該組結(jié)構(gòu)面幾何形態(tài)要素概率模型。現(xiàn)以現(xiàn)場(chǎng)量測(cè)的第3組結(jié)構(gòu)面樣本幾何形態(tài)要素?cái)?shù)據(jù)為例，建立巖體結(jié)構(gòu)面概率模型。在此，簡(jiǎn)單的以圖3和表3表達(dá)概率模型實(shí)現(xiàn)計(jì)算效果及相應(yīng)參數(shù)。

圖3 各幾何形態(tài)要素分布直方圖

表3 各幾何形態(tài)要素特征值

可以看出，結(jié)構(gòu)面傾角基本上服從正態(tài)分布，傾向接近于正態(tài)分布，跡長(zhǎng)基本服從負(fù)指數(shù)分布，間距也大致接近于負(fù)指數(shù)分布。

使用巖體節(jié)理分布Monte-Carlo計(jì)算機(jī)輔助模擬程序，根據(jù)各組結(jié)構(gòu)面幾何形態(tài)要素概率模型，生成與實(shí)測(cè)結(jié)構(gòu)面具有相同幾何形態(tài)要素及分布形式的切面結(jié)構(gòu)面網(wǎng)絡(luò)圖。圖4為巖體結(jié)構(gòu)面二維網(wǎng)絡(luò)模擬圖，模擬結(jié)果可以很好的顯示三組結(jié)構(gòu)面的分布與交錯(cuò)關(guān)系。

3 邊坡巖體模型

圖4 結(jié)構(gòu)面二維網(wǎng)絡(luò)模擬圖

坡體巖性主要為侏羅系上侏羅統(tǒng)東興組棕黃色長(zhǎng)石石英砂巖和紫紅色泥巖。其中，地表淺層巖體基本上全風(fēng)化為土層，屬于地表植被根系的作用層，其厚度小，對(duì)邊坡整體穩(wěn)定性不構(gòu)成影響；強(qiáng)風(fēng)化巖石也為散體狀，其性質(zhì)不受結(jié)構(gòu)面的控制，因此全風(fēng)化和其下部的強(qiáng)風(fēng)化帶的物理力學(xué)性質(zhì)按土層來(lái)確定，可不考慮結(jié)構(gòu)面的影響。由于邊坡體主體為軟弱～較軟巖體，屬?gòu)椝苄越橘|(zhì)，計(jì)算采用理想彈塑性本構(gòu)模型，只施加重力荷載，不考慮構(gòu)造應(yīng)力，屈服準(zhǔn)則采用Mohr-Coulomb準(zhǔn)則。計(jì)算邊界取到足夠遠(yuǎn)處，即人工邊坡開(kāi)挖后，引起邊坡應(yīng)力、應(yīng)變影響范圍以外。垂直邊界面水平方向加約束，水平邊界面垂直方向加約束，為固定支座。廠址區(qū)地震烈度為Ⅵ度，邊坡上沒(méi)有構(gòu)筑物，因此不考慮地震的影響，按靜力情況模擬。根據(jù)實(shí)測(cè)的典型工程地質(zhì)剖面，以人工開(kāi)挖形成1:1坡比進(jìn)行邊坡穩(wěn)定性計(jì)算。圖5核電站人工邊坡剖面巖體含三組結(jié)構(gòu)面網(wǎng)絡(luò)圖。

圖5 邊坡剖面巖體結(jié)構(gòu)面網(wǎng)絡(luò)圖

以有限元法邊坡穩(wěn)定性分析，用節(jié)理單元模型處理時(shí)，將巖體結(jié)構(gòu)面幾何模型轉(zhuǎn)化為計(jì)算模型，具體實(shí)現(xiàn)過(guò)程在此不再贅言。其中，節(jié)理單元接觸面法向剛度Kn取巖塊的抗剪強(qiáng)度、切向剛度Ks取巖塊彈性模量。

在邊坡區(qū)選取中風(fēng)化砂巖和泥巖兩類巖性發(fā)育三組結(jié)構(gòu)面的代表性統(tǒng)計(jì)窗1處，統(tǒng)計(jì)結(jié)果見(jiàn)表4。

表4 結(jié)構(gòu)面統(tǒng)計(jì)窗數(shù)據(jù)統(tǒng)計(jì)

同樣，用節(jié)理巖體損傷模型處理時(shí)，也可直接應(yīng)用有限元法中巖體損傷模型，對(duì)燈塔嶺人工邊坡概化模型處理。將現(xiàn)場(chǎng)量測(cè)的結(jié)構(gòu)面抗壓強(qiáng)度JCS、結(jié)構(gòu)面粗糙度系數(shù)JRC代入式Kni= -7.15+1.75JRC+(JCS/JRC)，可求得結(jié)構(gòu)面法向剛度Kni，再和幾何形態(tài)要素中節(jié)理盤半徑統(tǒng)計(jì)平均值α一同代入式即可得結(jié)構(gòu)面?zhèn)鲏合禂?shù)Cv，可見(jiàn)表5。

表5 結(jié)構(gòu)面統(tǒng)計(jì)窗數(shù)據(jù)統(tǒng)計(jì)

在此基礎(chǔ)上，進(jìn)行損傷張量計(jì)算可得到考慮結(jié)構(gòu)面影響傳壓和傳剪系數(shù)修正后的整體損傷張量，同時(shí)也可得到三個(gè)特征值及與之對(duì)應(yīng)的相互垂直的三個(gè)特征向量，依此相應(yīng)折減便可得到穩(wěn)定性計(jì)算時(shí)變形參數(shù)，結(jié)果見(jiàn)表6。

表6 巖體變形參數(shù)選取

采用RMR巖體分類對(duì)工程地質(zhì)巖組評(píng)分，然后按Hoek-Brown經(jīng)驗(yàn)方程在給定巖石的單軸抗壓強(qiáng)度σc條件下，以不同的σ3計(jì)算相應(yīng)的σm，用一組σ3和σ1m繪制摩爾圓及其包絡(luò)線，在包絡(luò)線上截取cm值，并量取φm。此方法即為對(duì)巖體抗剪強(qiáng)度的弱化處理，用弱化處理后的巖體抗剪強(qiáng)度和工程經(jīng)驗(yàn)綜合確定穩(wěn)定性計(jì)算時(shí)的強(qiáng)度參數(shù)。各類巖組計(jì)算參數(shù)取值見(jiàn)表7。

表7 巖體強(qiáng)度參數(shù)選取

4 人工邊坡穩(wěn)定性評(píng)價(jià)

在某核電站人工邊坡穩(wěn)定性評(píng)價(jià)中，以有限元法兩種不同的模型處理，圖6和圖7分別為兩種模型穩(wěn)定性分析圖。無(wú)論是節(jié)理單元模型，還是損傷模型都清楚地顯示了邊坡滑動(dòng)帶。同時(shí)，自坡腳起大致圓弧形的一定范圍內(nèi)變形很大，所以節(jié)理單元模型位移等值線圖在此區(qū)域很密集，而損傷模型應(yīng)力矢量圖表示的應(yīng)力矢量同樣也在近圓弧形范圍內(nèi)自坡腳起向外逐漸減小。在巖體邊坡穩(wěn)定性評(píng)價(jià)中，兩種數(shù)值計(jì)算模型不同之處的最直接體現(xiàn)就是穩(wěn)定系數(shù)大小的差別，這種穩(wěn)定系數(shù)之間大小差別反映到最危險(xiǎn)滑動(dòng)面搜索結(jié)果上就是滑動(dòng)面位置和范圍的不同。對(duì)于不同滑動(dòng)面的差異可以解釋為是由兩種不同數(shù)值計(jì)算模型本身方法的差異造成，以節(jié)理單元模型考慮時(shí)只是對(duì)整個(gè)坡體范圍內(nèi)的結(jié)構(gòu)面進(jìn)行特殊的單元剖分處理，有限元法計(jì)算時(shí)可以看作每條結(jié)構(gòu)面和周圍巖體建立彈塑性力學(xué)方程進(jìn)行力與位移的計(jì)算；而損傷模型對(duì)結(jié)構(gòu)面的折減則是建立模型時(shí)過(guò)渡為對(duì)整個(gè)坡體范圍內(nèi)巖體強(qiáng)度和變形的折減，所以巖體強(qiáng)度和變形的降低也要在有限元法計(jì)算時(shí)相應(yīng)的引起潛在滑動(dòng)面的“擴(kuò)展”。

圖6 節(jié)理單元模型位移分布圖

圖7 節(jié)理巖體損傷模型應(yīng)力矢量圖

事實(shí)上，對(duì)于考慮受結(jié)構(gòu)面影響的巖體邊坡穩(wěn)定性評(píng)價(jià)采用有限元法計(jì)算時(shí)，有限元計(jì)算受數(shù)值模擬方法限制，文中模擬進(jìn)行了一定的理想化處理，所以效果相對(duì)較好。除了在模擬中參數(shù)還有待進(jìn)一步的探討外，穩(wěn)定性評(píng)價(jià)時(shí)的有限元模擬最重要的還是對(duì)于當(dāng)前使用的有限元軟件其內(nèi)部計(jì)算核心信息一無(wú)所知，正是這種知識(shí)上的缺失很大程度上限制了研究的深度，因此，自用有限元軟件的開(kāi)發(fā)至關(guān)重要，從這一意義上說(shuō)，對(duì)于節(jié)理化巖體邊坡穩(wěn)定性研究無(wú)論空間還是內(nèi)容都很廣闊，還有很多其涉及到的有意義的工作需要進(jìn)行。

[1]孫衛(wèi)軍，周維垣.裂隙巖體彈塑性—損傷本構(gòu)模型[J].巖石力學(xué)與工程學(xué)報(bào)，1990，9(2).

[2]范雷，王亮清，唐輝明.節(jié)理巖體結(jié)構(gòu)面產(chǎn)狀的動(dòng)態(tài)聚類分析[J].巖土力學(xué)，2007，28(11).

[3]賈洪彪，等.巖體結(jié)構(gòu)面三維網(wǎng)絡(luò)模擬理論與工程應(yīng)用[M].北京：科學(xué)出版社，2008.

[4]凌建明.節(jié)理裂隙巖體損傷力學(xué)研究中的若干問(wèn)題[J].力學(xué)進(jìn)展，1994，25(2).

[5]楊延毅.節(jié)理裂隙巖體損傷—斷裂力學(xué)模型及其在巖體工程中的應(yīng)用[D].北京：清華大學(xué)，1990.

[6]周維垣，楊延毅.節(jié)理巖體力學(xué)參數(shù)取值研究[J].巖土工程學(xué)報(bào)，1992，14(5).

[7]李同錄，等.節(jié)理巖體力學(xué)參數(shù)的選取與應(yīng)用[J].巖石力學(xué)與工程學(xué)報(bào)，2004，23(13).

[8]陳征宙，等.巖體節(jié)理網(wǎng)絡(luò)模擬技術(shù)研究[J].巖土工程學(xué)報(bào)，1998，20(1).