基于三維離散元法的地下洞室定位塊體地震響應

朱付廣,王 禹,宋康林,劉繼國,朱澤奇

(1.湖北省電力勘測設計院,湖北武漢 430040;2.中交第二公路勘察設計研究院有限公司,湖北武漢 430056; 3.中國科學院武漢巖土力學研究所,湖北武漢 430071)

基于三維離散元法的地下洞室定位塊體地震響應

朱付廣1,王 禹1,宋康林2,劉繼國2,朱澤奇3

(1.湖北省電力勘測設計院,湖北武漢 430040;2.中交第二公路勘察設計研究院有限公司,湖北武漢 430056; 3.中國科學院武漢巖土力學研究所,湖北武漢 430071)

為分析定位塊體在不同支護工況下的地震響應特征與動態穩定性響應,以大崗山水電站地下廠房為工程實例,采用塊體理論和三維離散元軟件(3DEC)建立了安裝間某典型定位塊體數值模型,并將圍巖模擬開挖位移與實測變形監測數據進行對比,論證塊體與結構面力學參數及相關模型建立的合理性。對比該定位塊體在不同支護工況下的地震響應規律與動態穩定性,結果表明,應用離散單元方法分析塊體地震響應與穩定性是可行的,錨桿支護對定位塊體具有較好的抗震作用。

地下洞室;塊體理論;三維離散元法;地震響應;大崗山水電站

在不連續數值分析方法中,離散元與DDA、流形元等方法一起,是被認為適用于不連續巖體變形與穩定分析的有力工具[1]。離散元法是針對節理巖體提出的一種適用于模擬巖土體大位移的數值方法,與工程中常用的有限元方法相比,離散元法在三大定律的滿足上有所不同:①從平衡方程上看,離散元法采用牛頓第二定律控制,按圍繞各剛體塊體形心的力平衡和力矩平衡來滿足;②從變形協調方程來看,各剛體塊體間不再位移連續,而是允許大變形和斷裂分開,可以模擬巖體不連續結構面的滑移與開裂;③從本構方程上看,離散元法避開了復雜的本構關系推導,采用在剛性塊體間設置不同種類彈簧和阻尼(法向剛度和阻尼、切向剛度和阻尼)來反映材料的應力 位移關系[2]。離散元法避開大規模剛度矩陣的形成,采用動態松弛法求解,可以高效地求解非線性大變形問題。

在地震動載荷作用下,地震波將在斷層、巖脈等結構面上產生復雜的透射、折射、反射等傳播現象,個別部位的巖體就可能產生松動滑移甚至失穩垮落。節理巖體地下洞室開挖的穩定性問題一直是巖土工程界比較關注的問題,但對于地下廠房中由結構面和高邊墻、大跨度頂拱切割形成的可動塊體,對其動力特性和穩定性目前研究較少。周劍等[3]運用離散元軟件研究了結構面產狀與節理數量影響下應力波傳播規律,結果表明結構面數量越多,應力波穿過后造成動力響應越弱,驗證了斷裂的隔震效應與上盤效應;楊志華等[4]將強震作用下斷層破裂及其引起的隧道圍巖破壞過程定性地劃分為5個主要階段:彈性應力集中階段、破裂起始階段、破裂加速階段、穩定破裂階段和破裂逐漸停止階段;王光綸等[5-7]對剛性塊體進行了動力試驗和離散元法動力分析參數選擇的研究,得到了地震影響下剛性塊體穩定性的一些非常有價值的結論;石崇等[8-10]運用二維離散元與三維離散元軟件模擬了巖質邊坡失穩過程,結果表明高烈度地震作用下巖質邊坡破壞是從表層局部塊體失穩開始的;李海波等[11-12]針對地震荷載作用下順層巖質邊坡安全系數進行了分析研究,給出了在離散元中判斷塊體滑動變形最終穩定與否的方法。

本文以大崗山水電站地下廠房施工期揭露地質資料為基礎,運用三維離散元軟件(3DEC)建立安裝間某典型定位塊體的數值模型,研究定位塊體在不同支護工況下的地震響應特征與動態穩定性。

1 工程概況及模型建立

大崗山水電站壩址位于四川省大渡河中游上段雅安市石棉縣境內。電站裝機容量為2600MW。引水、發電建筑物由電站進水口、4條壓力管道、地下廠房、主變室、尾水調壓室、2條尾水洞等組成。地下廠房系統采用主副廠房、主變室、尾水調壓室3大洞室平行布置,主廠房開挖尺寸為206.00 m×30.80 m× 73.78m(長×寬×高),主變室開挖尺寸為144.00 m× 18.80 m×25.10 m(長×寬×高),尾水調壓室長為130.00m,凈跨度為20.50～24.00m,室高為75.08m。廠區巖體新鮮較完整,呈塊狀～次塊狀結構,巖塊嵌合緊密,有斷層f57、f58、f59和f60等較大規模的軟弱結構面和延伸較長的β6、β80、β81等輝綠巖脈穿過廠房洞群區。廠區裂隙(Ⅴ級結構面)優勢產狀大致有5組,總體上以近SN/E向、近EW/N(或S)向陡傾角結構面最為發育。

工程區位為SN向、NW向和NE向多組斷裂構造交匯部位,區域地質構造及地震地質背景復雜,新構造活動較強烈,其區域構造穩定性對工程建設影響較大。水電站壩址50年超越概率10%基巖水平向峰值加速度為0.251g,相應地震基本烈度為8度。安裝間頂拱施工期間揭露較大型巖脈如圖1所示,由圖1可直觀看出,巖脈βc3、βc4、βc5相互切割,形成了一體積達46.87 m3、埋深12.83 m的定位塊體,其施工期和運營期穩定性引起了設計和施工方的極大關注,有鑒于此,本文利用3DEC建立安裝間定位塊體模型,模型尺寸330 m×330 m×60 m(分別為x、y、z向尺寸,如圖2所示),通過模擬開挖和施加支護與地震荷載,研究該定位塊體開挖穩定性以及地震響應規律,為地下廠房施工期和運營期穩定性評價提供分析資料和依據。

圖1 安裝間巖脈地質示意圖

圖2 三維離散元模型

計算模型中巖體材料選用各向同性彈性本構模型,結構面采用面接觸的Coulomb滑動模型,巖體為II類圍巖,物理力學指標如下:密度2 650 kg/m3,彈性模量20 GPa,泊松比0.25;結構面物理力學指標如下:法向剛度20 GPa/m,切向剛度10 GPa/m,內摩擦角35°,黏聚力0.5 MPa。

2 模型驗證

地下洞室開挖順序見圖3。每開挖一步,迭代計算直至整個計算模型穩定為止,然后再進行下一步的開挖。洞室開挖后發生指向臨空面的位移,洞周位移云圖如圖4所示?？梢钥闯?拱頂開挖完成后最大位移值為2.276 cm。

圖3 地下廠房開挖順序示意圖(單位:m)

在上層排水廊道具備鉆孔施工條件后,由排水廊道向主廠房鉆孔,施工時在監測孔CSP2-1中埋設了滑動測微計監測管,此時主廠房僅進行了中導洞部分的先期開挖,如圖5所示。

圖4 完全開挖后洞周位移云圖

圖5 監測孔布置示意圖

將該洞段附近安裝的滑動測微計位移監測結果與該計算剖面相應部位的位移進行對比,圖6給出了主廠房上游邊墻部位CSP2-1監測孔利用滑動測微計實測位移和無支護條件下的計算位移過程線。從安裝間開挖完成后洞周圍巖變形的數量級來看,可以認為圍巖在洞室開挖后未發生大的變形。位移實測值與計算值大小在同一個量級,相差不大,且隨開挖施工過程的變化趨勢也基本一致,說明本文所建模型和相關力學參數取值合理,可以作為下一步地震響應研究的基礎。

圖6 監測與計算位移過程線

3 地震響應分析

3.1 無支護塊體地震穩定性計算與分析

由于工程區域缺乏現有的實測地震波數據,采用與大崗山水電站斷裂帶和地質背景相似的汶川地震中郫縣走石山地震臺記錄的地震波來研究該電站地下廠房塊體的地震響應。利用低頻濾波器將輸入地震波時程中高頻部分濾去,并經基線校正后得到如圖7所示的地震速度時程曲線。地下洞室開挖模擬達到穩定狀態后,將模型四周(水平向)設置為黏滯邊界,同樣將計算模型底部邊界(豎直向)設置為黏滯邊界,并輸入校正后的汶川地震中郫縣走石山3個方向上的地震動作用,開展定位塊體的地震響應計算,如圖8所示。

圖7 郫縣走石山3個方向的地震速度波

圖8 地震響應分析計算模型

在離散元法地震反應分析中,可通過監測該定位塊體質心部位的絕對位移響應時程和該定位塊體質心與組成該定位塊體四周結構面的相對滑動位移規律,來分析研究定位塊體的地震響應特征以及地震穩定性。圖9給出了定位塊體在郫縣走石山地震波作用下的3個方向上的絕對位移響應時程曲線。可以看出該定位塊體在位移地震響應過程中整體上呈與圍巖共同振動的趨勢,但最終3個方向上都產生了一定的殘余位移。

圖9 定位塊體3個方向的位移響應時程

為了進一步分析該殘余位移產生的過程,圖10給出了該定位塊體沿結構面滑動的位移時程曲線(位移為塊體對應于圍巖的位移,下同)。通過對塊體地震工況下產生的永久位移分段分析可知:由于地震波起始階段振幅較小,因而塊體AB段位移較小;地震波輸入約20 s時3個方向上都幾乎達到峰值,導致塊體在BC段的滑動位移產生了躍遷;地震波在25～45 s時3個方向都處于劇烈的振動階段,由此導致塊體CD段的位移也產生較大的波動;地震波在DE段趨于平緩,塊體的位移相應也趨于平穩,并最終產生5cm左右永久位移。3DEC軟件的計算未能考慮巖土力學參數的退化效應,故塊體滑動了5 cm仍能在兩個滑動面阻滑力作用下保持穩定,但廠房拱頂定位塊體5 cm的變形已經影響正常使用,屬于破壞狀態。

圖10 定位塊體的滑動位移時程

3.2 支護條件下塊體地震穩定性計算與分析

依據洞室支護設計資料實際情況(圖11),對該定位塊體施加2根長9m的預應力錨桿和2根長6 m的普通砂漿錨桿,布錨間距都為1.2 m,進行支護后的地震穩定性計算,其中錨桿是通過3DEC軟件中的命令“struct”與關鍵詞“cable”和“prop”施加的;圖12為支護后該定位塊體的地震滑動位移響應時程曲線。可將施加支護后的塊體位移時程劃分為3個階段:AB段時地震作用較小,幾乎沒有產生位移;BC段隨著地震作用增大,該定位塊體相對位移在2.5 mm左右震蕩,且位移波動幅度較大,但由于錨桿力的約束,該定位塊體相對位移沒有產生圖11所示的躍遷;CD段地震波逐漸減弱,塊體的位移也趨于平穩,最終在錨桿力、塊體下滑力與兩個滑面阻滑力綜合平衡下產生了1 mm左右的殘余位移。對比施加支護前后的相對位移響應規律,說明錨桿對于定位塊體具有較好的抗震作用。

圖11 定位塊體錨桿支護

圖12 支護條件下定位塊體的滑動位移時程

4 結 語

本文以大崗山水電站地下廠房為工程實例,以在施工期揭露地質資料為基礎,運用3DEC建立安裝間某典型定位塊體的數值模型,開展了安裝間的開挖穩定性分析以及定位塊體的地震響應分析,比較了定位塊體在支護條件下的位移響應特點和規律,結果表明應用離散單元方法分析定位塊體地震響應與穩定性是可行的,錨桿對于定位塊體具有較好的抗震作用。

[1]周家文,徐衛亞,石崇.基于3DEC的節理巖體邊坡地震影響下的楔體穩定性分析[J].巖石力學與工程學報, 2007,26(增刊1):3402-3409.(ZHOU Jiawen,XU Weiya, SHI Chong.Wedge stability analysis of jointed rockmass slope considering seismic influences based on 3DEC[J]. Chinese Journal of Rock Mechanics and Engineering,2007, 26(Sup1):3402-3409.(in Chinese))

[2]王泳嘉,邢紀波.離散單元法及其在巖土力學中的應用[M].沈陽:東北工學院出版社,1991.

[3]周劍,張路青,胡瑞林,等.大型結構面產狀影響下應力波傳播規律研究[J].巖土力學與工程學報,2011, 30(4):769-780.(ZHOU Jian,ZHANG Luqing,HU Ruilin,et al.Study of rules of stress waves propagation under various attitudes of large-scale fractures.[J]. Chinese Journal of Rock Mechanics and Engineering. 2011,30(4):769-780.(in Chinese))

[4]楊志華,蘭恒星,張永雙,等.強震作用下穿越斷層隧道圍巖力學響應研究[J].工程地質學報,2013,21(2): 171-181.(YANG Zhihua,LAN Hengxing,ZHANG Yongshuang,et al.Rock dynamic response of railway tunnel traversing through fault under strong earthquake. [J].Journal of Engineering Geology,2013,21(2):171-181.(in Chinese)) [5]王光綸,張楚漢,彭岡,等.剛塊體動力試驗與離散元法動力分析參數選擇的研究[J].巖石力學與工程學報1994,13(2):124-133.(WANG Guanglun,ZHANG Chuhan,PENG Gang,et al.Dynamic test of rigid blocks and parameter study in dynamic analysis of distinct elements.[J].Chinese Journal of Rock Mechanics and Engineering,1994,13(2):124-133.(in Chinese))

[6]李世海,汪元年.三維離散元計算參數選取方法研究[J].巖石力學與工程學報2004,23(21):3642-3651. (LIShihai, WANG Yuannian.Selection study of computational parameters for DEM in geomechanics.[J]. Chinese Journal of Rock Mechanics and Engineering, 2004,23(21):3642-3651.(in Chinese))

[7]王衛華,李夕兵,胡盛斌.模型參數對3DEC動態建模的影響[J].巖石力學與工程學報,2005,24(增刊1): 4790-4797.(WANG weihua,LI xibing,HU Shengbin. Effect of model parameters on 3DEC dynamic modeling. [J].Chinese Journal of Rock Mechanics and Engineering, 2005,24(Sup1):4790-4797.(in Chinese))

[8]石崇,王盛年,劉琳.地震作用下陡巖崩塌顆粒離散元數值模擬研究[J].巖石力學與工程學報,2013,32(增刊1):2798-2805.(SHI Chong,WANG Shengnian,LIU Lin.Research of avalanche disaster numerical simulation based on granular discrete element method of high-steep slope under seismic loads[J].Chinese Journal of Rock Mechanics and Engineering,2013,32(Sup1):2798-2805.(in Chinese))

[9]賀續文.基于離散單元法的節理巖體邊坡穩定性分析[D].湘潭:湘潭大學,2010.

[10]李新坡,何思明.節理巖質邊坡破壞過程的PFC2D數值模擬分析[J].四川大學學報:工程科學版,2010,42 (增刊 1)):70-75.(LI Xinpo,HE Siming.Numerical analysis of the failure of heavily jointed rock slopes using PFC2D[J].Journal of Sichuan University:Engineering Science,2010,42(Sup1):70-75.(in Chinese))

[11]李海波,肖克強,劉亞群.地震荷載作用下順層巖質邊坡安全系數分析[J].巖石力學與工程學報,2007,26 (12):2385-2393.(LI Haibo,XIAO Keqing,LIU Yaqun.Factor of safety analysis of bedding rocking slope under seismicload.[J].Chinese JournalofRock Mechanics and Engineering,2007,26(12):2385-2393. (in Chinese))

[12]李海波,蔣會軍,趙堅,等.動荷載作用下巖體工程安全的幾個問題[J].巖石力學與工程學報,2003,22(11): 1887-1891.(LI Haibo,JIANG Huijun,ZHAO jian,et al. Some problems about safety analysis of rock engineering under dynamic load.[J].Chinese Journal of Rock Mechanics and Engineering,2003,22(11):1887-1891. (in Chinese))

Seismic response assessment of positioning block in underground cavern based on 3D DEM//

ZHU Fuguang1, WANG Yu1,SONG Kanglin2,LIU Jiguo2,ZHU Zeqi3
(1.Hubei Electric Power Survey&Design Institute,Wuhan 430040; 2.China Second Highway Survey,Design and Research Institute of Communications,Wuhan 430056,China;3.Institute of Rock and Soil Mechanics,Chinese Academy of Sciences,Wuhan 430071,China)

In this paper,we analyze the earthquake response characteristics and dynamic stability of positioning block under different support conditions by taking the underground plant of Dagangshan Hydropower Station as an example.In doing it so,a typical positioning block numerical model of erection bay based on the block theory and three-dimensional distinct element code(3DEC)have been built.Additionally,we perform a comparison between simulative dig displacement of surrounding rock and real deformation monitoring data,then demonstrates the parameter of block and structural planes and the feasibility of building relative models.Overall,the comparisons of earthquake response characteristics and dynamic stability of positioning block under different support conditions show that it is feasible to analyze block earthquake and stability with discrete element method.Additionally,the bolting support presents a better antiearthquake effect for positioning block.

underground cavern;key block;3D DEM;seismic response;Dagangshan Hydropower Station

TV731.6;P315.9

:A

:1006-7647(2014)04-0020-04

10.3880/j.issn.1006-7647.2014.04.005

2013-0717 編輯:熊水斌)

國家自然科學基金(51009130,51279202)

朱付廣(1983—),男,江蘇鹽城人,工程師,碩士,主要從事巖土力學參數和圍巖穩定性研究。E-mail:331584689@qq.com