某水電站邊坡卸荷巖體宏觀力學(xué)參數(shù)反演對比研究

蔡 健 ，劉 杰

（1.三峽大學(xué)土木與建筑學(xué)院，湖北 宜昌 443002；2.湖北省水利水電規(guī)劃勘測設(shè)計院，湖北 武漢 430064）

0 引 言

巖體宏觀力學(xué)參數(shù)在有限元計算中占有十分重要的地位，沒有正確的力學(xué)參數(shù)，數(shù)值模擬計算結(jié)果不能令人信服。由于巖體試驗(yàn)尺寸及現(xiàn)場選點(diǎn)條件的限制，復(fù)雜巖體的力學(xué)參數(shù)很難通過室內(nèi)外試驗(yàn)直接測得。本文在總結(jié)前人研究成果的基礎(chǔ)上，基于巖體卸荷原理和分級分區(qū)方法的基礎(chǔ)上，對某水電站邊坡區(qū)域進(jìn)行了合理的劃分，并選取具有代表性的區(qū)域進(jìn)行二維與三維參數(shù)反演對比分析，最終得到卸荷巖體宏觀力學(xué)參數(shù)，并進(jìn)一步驗(yàn)證了本文研究方法的精確性和合理性。

1 基本方法

1.1 反演方法

以現(xiàn)場地質(zhì)調(diào)查、材料試驗(yàn)為基礎(chǔ)，與數(shù)值計算相結(jié)合，根據(jù)巖體結(jié)構(gòu)面的分布狀況、級別及巖體的不同性質(zhì)，對巖體力學(xué)參數(shù)采用分區(qū)分級計算，最終對整個工程巖體力學(xué)參數(shù)進(jìn)行評價。具體計算步驟如下:

（1）根據(jù)工程地質(zhì)勘探成果，對給定的整個工程巖體各個區(qū)域分區(qū)分塊。分區(qū)分塊主要考慮巖石的類型及強(qiáng)度、風(fēng)化程度、地應(yīng)力分布、結(jié)構(gòu)面的產(chǎn)狀、間距及連通率、力學(xué)參數(shù)、地下水情況等。

（2）根據(jù)各個區(qū)塊內(nèi)所含節(jié)理的尺寸，對各區(qū)塊劃分不同尺寸級別，數(shù)值模擬計算試塊全面、真實(shí)地反映其對應(yīng)巖體的巖性組合和結(jié)構(gòu)特征。

（3）根據(jù)每個區(qū)塊的邊界和外荷條件，對各區(qū)塊內(nèi)的分級計算試塊逐級進(jìn)行大三軸試驗(yàn)數(shù)值模擬，對每個區(qū)塊的巖體參數(shù)進(jìn)行評定。計算中應(yīng)模擬區(qū)塊內(nèi)每條裂隙、每個層面。對同一區(qū)塊不同級別計算試塊，計算順序是先小尺寸后大尺寸，先小裂隙后大裂隙。第1級試塊內(nèi)為巖石和小節(jié)理裂隙材料，計算參數(shù)和本構(gòu)關(guān)系直接采用室內(nèi)試驗(yàn)成果，通過第1級計算得到試塊材料的力學(xué)參數(shù)。第2級計算試塊內(nèi)含第1級計算試塊材料和中等裂隙材料。根據(jù)以上方法，依次計算各區(qū)塊內(nèi)不同級別的試塊，直至計算到該區(qū)塊最大級別試塊，最大試塊計算所得的參數(shù)就為該區(qū)塊的力學(xué)參數(shù)值。

（4）把各區(qū)塊和大斷層及軟弱厚層組合在一起，建立整體計算試塊，計算得到工程整體的力學(xué)參數(shù)，對整個工程的巖體進(jìn)行評價。

1.2 計算公式

在主應(yīng)力空間中，彈性狀態(tài)應(yīng)力、應(yīng)變關(guān)系服從廣義虎克定律。在平面應(yīng)變狀態(tài)下，應(yīng)變增量Δε3=0，計算模型在應(yīng)力增量Δσ1、Δσ2的作用下，通過有限元計算可得到Δε1和Δε2，則可求出泊松比μ、變形模量E。由所求的變形參數(shù)，可得到單向應(yīng)力狀態(tài)下的應(yīng)變增量和應(yīng)力增量之間的關(guān)系式。三維有限元分區(qū)分級參數(shù)模擬同二維有限元，根據(jù)廣義虎克定律推導(dǎo)出泊松比μ′和變形模量E′計算公式為

2 反演分析

2.1 區(qū)塊劃分

根據(jù)邊坡巖體的類型及強(qiáng)度、風(fēng)化程度、剪切破碎帶、裂隙及斷層的走向和傾向、結(jié)構(gòu)面的間距及連通率、力學(xué)參數(shù)、地下水情況等因素將邊坡劃分為不同的區(qū)塊，選取其中3個區(qū)塊 （2、3號和4號）進(jìn)行二維、三維參數(shù)反演對比分析，其區(qū)塊劃分見圖1。

圖1 區(qū)塊劃分

（1）第2區(qū)塊:位于左岸弱下風(fēng)化層內(nèi)邊坡上，巖性主要為灰白色英安巖，巖體呈鑲嵌狀、塊裂狀結(jié)構(gòu)，結(jié)構(gòu)面發(fā)育，貫穿性結(jié)構(gòu)面不多見，結(jié)構(gòu)面延展差，多閉合，巖塊間嵌合力較好。

（2）第3區(qū)塊:位于左岸弱上風(fēng)化層內(nèi)邊坡上，巖性主要為灰白色英安巖，巖體呈鑲嵌狀、塊裂狀結(jié)構(gòu)，結(jié)構(gòu)面較發(fā)育。

（3）第4區(qū)塊:位于右岸弱上風(fēng)化層內(nèi)邊坡上，巖性主要為灰黑色英安巖，巖體呈鑲嵌狀結(jié)構(gòu)，結(jié)構(gòu)面較發(fā)育。

2.2 區(qū)塊二維參數(shù)反演

根據(jù)裂隙大小、性質(zhì)、產(chǎn)狀等因素，將第2區(qū)塊分4級計算。其概化模型試件見圖2，初始參數(shù)見表1。

在ADINA中建立好宏觀力學(xué)參數(shù)分析模型后，設(shè)定材料參數(shù)進(jìn)行計算，運(yùn)行后處理文件，提取出節(jié)點(diǎn)數(shù)據(jù)，對節(jié)點(diǎn)應(yīng)力、應(yīng)變進(jìn)行統(tǒng)計。最后，用巖體初始計算參數(shù)計算得到模擬試塊的宏觀變形參數(shù)。第2區(qū)塊巖體二維宏觀參數(shù)計算結(jié)果見表2。同理，第3、第4區(qū)塊巖體二維宏觀參數(shù)計算結(jié)果見表3、4。

2.3 區(qū)塊三維參數(shù)反演

三維有限元參數(shù)反分析也采用二維參數(shù)分析中各區(qū)域巖體與破碎帶的力學(xué)參數(shù)作為數(shù)值模擬的基本參數(shù)，將幾種不同的巖體和剪切破碎帶斷層等組合在一起構(gòu)成計算試塊，并對其進(jìn)行地質(zhì)力學(xué)參數(shù)反分析，計算試塊的選定也與三維有限元模型結(jié)合起來考慮。進(jìn)行三維模擬區(qū)塊斷層和結(jié)構(gòu)面過程中，各級區(qū)塊大小以及結(jié)構(gòu)面的分布和發(fā)育狀況均與二維模型保持一致，故此處不再贅述。

表1 第2區(qū)塊巖體的初始參數(shù)

圖2 第2區(qū)塊計算概化二維模型

表2 第2區(qū)塊巖體的宏觀參數(shù)計算成果

表3 第3區(qū)塊巖體的宏觀參數(shù)計算成果

表4 第4區(qū)塊巖體的宏觀參數(shù)計算成果

根據(jù)裂隙大小，性質(zhì)，產(chǎn)狀等因素，將第2區(qū)塊分為4級計算。其分級概化模型試件見圖3，初始參數(shù)見表1。在ADINA中建立好宏觀力學(xué)參數(shù)分析模型后，設(shè)定材料參數(shù)進(jìn)行計算，運(yùn)行后處理文件，提取出節(jié)點(diǎn)數(shù)據(jù)，對節(jié)點(diǎn)應(yīng)力、應(yīng)變進(jìn)行統(tǒng)計，最后將統(tǒng)計得到的數(shù)據(jù)代入公式1、2，計算得到模擬試塊的宏觀變形參數(shù)，計算結(jié)果見表5。同理，第3、第4區(qū)塊巖體三維宏觀參數(shù)計算結(jié)果見表6、7。

圖3 第2區(qū)塊分級計算概化三維模型

表5 第2區(qū)塊巖體三維宏觀參數(shù)計算成果

表6 第3區(qū)塊巖體三維宏觀參數(shù)計算成果

表7 第4區(qū)塊巖體三維宏觀參數(shù)計算成果

2.4 對比分析

3個區(qū)塊二維與三維模擬計算變形模量、泊松比對比見圖4。從圖4可以看出，二維和三維模擬裂隙對計算結(jié)果的影響程度幾乎一致。進(jìn)行三維模擬時，變形模量的折減幅度比二維模擬略大，趨勢基本一致。同樣，泊松比的變化趨勢基本一致，均隨著計算級數(shù)的增大而增大，部分存在跳躍現(xiàn)象，其原因可能是:在三維模型中，由于結(jié)構(gòu)面的存在，隨著頂面荷載的施加，另外2個側(cè)向臨空面中間會出現(xiàn)應(yīng)變集中凸出，使該面的平均位移增加，即側(cè)向應(yīng)變增大，故反求得出的泊松比隨之減小。

圖4 二維與三維計算變形模量、泊松比對比

3 結(jié)論

（1）本文采用分區(qū)分級計算試塊法解決了室內(nèi)試驗(yàn)對巖體參數(shù)評價的局限性和以經(jīng)驗(yàn)為主的方法提供巖體的力學(xué)參數(shù)對巖體參數(shù)選取的隨機(jī)性。

（2）綜合考慮巖體的結(jié)構(gòu)類型與特點(diǎn)，將室內(nèi)試驗(yàn)成果運(yùn)用到數(shù)值仿真分級模擬巖體的參數(shù)，由此得出的反演結(jié)果更具科學(xué)性。

（3）從二維到三維參數(shù)反演，精確地模擬巖體內(nèi)部真實(shí)應(yīng)力、應(yīng)變狀態(tài)，大大提高了巖體參數(shù)分析的精確性。

（4）二維與三維的計算結(jié)果比較吻合，更加驗(yàn)證了分區(qū)分級方法計算巖體力學(xué)參數(shù)的準(zhǔn)確性和合理性。

［1］李建林，哈秋齡.三峽卸荷巖體宏觀力學(xué)參數(shù)三維數(shù)值模擬［J］.武漢水利電力大學(xué) （宜昌）學(xué)報，1997，9（3）:1-6.

［2］楊學(xué)堂，哈秋船，張永興，等.裂隙巖體宏觀力學(xué)參數(shù)數(shù)值仿真模擬研究［J］.水力發(fā)電，2004，30（7）:14－16.

［3］周濟(jì)芳，李建林，劉杰，等.隔河巖水電站邊坡巖體宏觀力學(xué)參數(shù)的有限元反演分析［J］.巖石力學(xué)與工程學(xué)報，2004，23（增1）:4506-4508.

［4］帥春，何江達(dá)，謝紅強(qiáng)，等.某水電站樞紐區(qū)邊坡天然穩(wěn)定性及巖體力學(xué)參數(shù)反演［J］.成都理工大學(xué)學(xué)報:自然科學(xué)版，2009，36（5）:511-515.

［5］李建林，王樂華，楊學(xué)堂，等.錦屏水電站壩廠區(qū)卸荷巖體力學(xué)參數(shù)研究［J］.巖土力學(xué)，2004，25（增）:17-20.

［6］駱世威，劉杰，肖蕾，等.西藏如美水電站裂隙巖體力學(xué)參數(shù)模擬反分析［J］.人民長江，2012，43（10）:84-87.

［7］李建林.巖石邊坡卸荷巖體宏觀參數(shù)研究［M］.北京:中國建筑工業(yè)出版社，1996.

［8］張志強(qiáng).非貫通裂隙巖體破壞細(xì)觀特征及其宏觀力學(xué)參數(shù)確定方法［D］.西安:西安理工大學(xué)，2009.

［9］李建林.三峽工程永久船間陡高邊坡巖體宏觀力學(xué)參數(shù)研究［D］.重慶:重慶建筑大學(xué)，1996.