楊繼華，盛 謙，朱澤奇，冷先倫

（1.中國科學院武漢巖土力學研究所 巖土力學與工程國家重點實驗室，武漢 430071；2.黃河勘測規(guī)劃設計有限公司，鄭州 450003）

1 引 言

我國西南地區(qū)興建了大批的水電工程，由于地形條件的限制，大都采用了地下式廠房，并形成了大跨度、高邊墻的地下洞室群。同時，該地區(qū)處于高地震烈度區(qū)，如果發(fā)生強震，地下洞室群有發(fā)生破壞的可能，如在2008年的汶川地震中，多座山嶺隧道和水工地下廠房發(fā)生破壞[1－2]。因此，研究大型地下巖體洞室群的抗震安全性十分重要。許多學者對此進行研究，并得出一些有益的結論。陳健云等[3]采用阻尼影響抽取法研究了圍巖動剛度的動力特性，提出了地下結構抗震分析的實用算法。李海波等[4]運用動應力集中因子代表值的概念，對溪洛渡水電站超大型地下洞室群安全性進行了評價。隋斌等[5]對地震荷載作用下大型地下洞室群進行動態(tài)模擬，采用新的劈裂判據對震后可能出現(xiàn)的劈裂損傷范圍進行了預測。王如賓等[6]研究了自然地震波作用下，有無襯砌工況下廠房洞室相對位移、點安全系數變化趨勢。

巖體在地震力作用下的響應是由靜態(tài)變成動態(tài)的力學過程，由于動力響應是一個復雜的非線性系統(tǒng)，考慮到地下巖體的失穩(wěn)是由于地震力的加、卸載所致，并且其演化過程具有明顯的非線性特征，因此，可借鑒地震預報領域的加、卸載響應比理論進行描述[7－9]。石崇等[10]建立了地震波穿越地下巖體洞室模型，采用動力時程分析法，驗證了加、卸載響應比法的正確性；姜彤等[11－12]將加、卸載比應用于巖體邊坡動力響應分析并取得了較好的效果。

本文應用加、卸載響應比理論，將作用在地下巖體洞群圍巖的地震力時程視為加、卸載過程，通過動力時程分析，對地下巖體洞室群圍巖的地震動力響應進行了一些有益的探索。

2 巖體地震響應的加、卸載響應比理論

2.1 加、卸載響應比理論

加、卸載響應比理論是由我國學者尹祥礎提出的[13－14]，最初應用于地震預測中。其主要思路是：非線性系統(tǒng)的失穩(wěn)前兆是系統(tǒng)在極其微弱的加載下產生巨大的響應，因此，可以利用系統(tǒng)的加荷響應率與卸荷響應率的比值定量地描述非線性系統(tǒng)偏離穩(wěn)定狀態(tài)（或接近失穩(wěn)狀態(tài)）的程度。

從巖體受力破壞的角度分析，加、卸載響應比理論就是將巖體受力過程中的荷載P和相應的某響應參數R的變化之比值ΔR/ΔP定義為響應率。由于ΔR/ΔP是一個有量綱的量，單位取的不同，其數值也會不同，這將引起許多不便，為此進一步引入無量綱化量——加、卸載響應比Y：

式中：Y為響應比；P為荷載；R為響應；下標+、-分別表示加載和卸載。

當巖體在彈性變形階段時，Y≈1；進入非線性變形階段后，Y＞1；巖體失穩(wěn)破壞時，Y→∞。根據這一規(guī)律，可以用Y這一參量來描述巖體趨向破壞失穩(wěn)的程度。

2.2 地下洞室群圍巖加、卸載響應參數的選取

在地震力作用下，地下洞室群圍巖的加、卸載過程是一個時間序列，同樣，圍巖的響應過程也是一個時間序列，通過對這兩個時間序列進行計算，即可獲得圍巖加、卸載響應比隨時間的變化規(guī)律。在巖體地震動力時程的加、卸載響應比分析中，一般將輸入的加速度時程作為加、卸載參數，將輸出的加速度時程、速度時程或位移時程作為響應參數[10－12]。由動力分析的原理可知[15]，輸出的速度時程和位移時程通常是由加速度時程通過1次和2次積分得到的，盡管地震波在輸入之前進行了修正，但一些周期較長的分量仍然會殘留下來，這對輸出的加速度時程沒有影響，但可能造成經過1次和2次積分得到的速度時程和位移時程的失真。如果以速度時程和位移時程作為加、卸載響應參數同樣會造成加、卸載響應比的失真，因此，采用加速度作為響應參數進行加、卸載響應比計算會更加合理。

2.3 加、卸載響應比計算方法

加、卸載過程一般采用區(qū)間直接平均法[10－11]，即將地震荷載周期性變化作為對地下洞室群圍巖的加、卸載手段，按地震持時T劃分n個加、卸載區(qū)間，將加、卸載區(qū)間內的地震荷載時間序列pi（加載時pi為正，卸載時pi為負）將洞室群圍巖在加、卸載區(qū)間內的加速度ai作為地下洞室圍巖對地震荷載的響應，應用加、卸載響應比理論確定圍巖的加、卸載響應比Y。

第i個區(qū)間的地震加、卸載和響應平均值為

根據加、卸載響應比理論，可得地震過程中的加、卸載響應比序列為

2.4 加、卸載響應區(qū)段的確定

一般來說，加、卸載區(qū)間的劃分應既能反映整個地震過程響應的特征，又能有效減少計算量，加、卸載區(qū)間劃分過少，將使計算結果趨于均化，且會漏掉峰值信息，不能真實反映地震過程響應變化過程；區(qū)間劃分過多，將極大增大計算量，不利于分析計算成果，另外還會出現(xiàn)因某一區(qū)間全為加載，計算的響應比為無窮大的情況。因此，區(qū)間劃分的確定實際上是一個均衡這兩方面因素的過程。本文經過反復試算，將加、卸載區(qū)間取為0.5 s。

3 工程背景

3.1 工程概況

白鶴灘水電站位于四川省寧南縣和云南省巧家縣交界的金沙江下游，其引水發(fā)電建筑采用地下式，其中廠區(qū)三大洞室主副廠房洞、主變洞、尾水調壓室平行布置，洞室尺寸均以長×寬×高給出，主副廠房洞尺寸為 439 m×32.2 m×78.5 m、439 m×29 m×78.5 m，主變洞尺寸為400 m×20.5 m×33.2 m，尾水調壓室尺寸為321.6 m×27.6 m×103.5 m。白鶴灘水電站地下廠房洞室群處在地震活動強烈的高山峽谷地區(qū)，地震基本烈度為Ⅷ度，工程區(qū)域地形地貌復雜，區(qū)域地應力水平較高，地質條件復雜，發(fā)育有多條層間錯動帶，其中層間錯動帶 C3、C4和C5對地下洞室群穩(wěn)定性影響較大。

3.2 數值模型

通過有限差分程序FLAC3D建立白鶴灘水電站13號機組剖面的數值模型。X軸與廠房軸線垂直，以指向下游方向為正；Y軸為廠房軸線方向；Z軸為鉛直方向，以向上為正。區(qū)域范圍為X向由主廠房中心線向上游方向約延伸450 m，向下游方向約延伸450 m，共900 m；由于采用假三維模擬二維的方法，因此，模型Y軸取2 m；Z向由高程300 m延伸至地表。在研究區(qū)域內發(fā)育的巖層主要有新鮮狀隱晶質玄武巖、斑狀玄武巖夾杏仁玄武巖、玄武質角礫熔巖等，巖質堅硬，巖體新鮮較完整，巖塊嵌合緊密，無區(qū)域斷裂切割，構造型式以斷層、節(jié)理裂隙、層間錯動帶為主要特征，壩區(qū)巖體可分為Ⅱ、Ⅲ、Ⅳ、Ⅴ共4大類。模型中僅考慮對地下洞室群圍巖穩(wěn)定影響較大的Ⅱ、Ⅲ類圍巖及層間錯動帶C3、C4和C5，忽略位于坡面附近影響較小的全、強、弱風化層。層間錯動帶C3、C4和C5采用軟弱夾層模擬。在模型中沒有考慮襯砌的影響，為滿足動力時程計算對單元尺寸的要求，將Ⅱ、Ⅲ類圍巖處最大單元尺寸取為10 m，將層間錯動帶處最大網格尺寸取為2 m。模型網格見圖1。巖體及層間錯動帶的物理力學參數見表1。

3.3 地震動輸入

在地震動力時程分析中，輸入的地震波最好是本場區(qū)的實測地震記錄，但白鶴灘水電站場地沒有實測的地震記錄可以選用，因此，輸入的地震波采用與白鶴灘水電站場區(qū)斷裂帶和地質背景相似的2008年汶川地震中綿竹清平臺站處的實測地震記錄。根據地震危險性分析，白鶴灘水電站場區(qū) 100年超越概率2%水平向峰值加速度為3.4 m/s2，經過峰值加速度調整后的汶川地震波時程曲線如圖2所示。根據地殼基本結構及固體介質波的折射規(guī)律，在地震響應分析中通常認為從震源發(fā)出的地震波到達地殼表層時入射方向垂直于地表，故在計算分析時在模型邊界上輸入地震波的入射方向垂直于底邊表面，同時考慮到地震波的剪切作用對地下洞室群破壞作用最大，將入射地震波設置為水平向振動的S波。邊界條件設置為黏性邊界。選用瑞利阻尼，通過對巖體動態(tài)響應的速度譜分析，在數值分析時阻尼比取為0.05，中心頻率取為1.0。

圖1 13號機組剖面數值模型Fig.1 Numerical model of cross-section No.13

表1 巖體物理力學參數Table1 Physico-mechanical parameters of rock masses

圖2 汶川地震波加速度時程曲線Fig.2 Acceleration time-history curve of Wenchuan earthquake

4 計算結果分析

洞室在地震作用下的破壞多表現(xiàn)為相互錯動、擠壓、彎曲，因此，可以采用洞室內壁質點的響應特征來衡量洞室受地震影響程度[16]。通過動力時程計算，可以得出洞室各關鍵部位的加速度、位移和主應力差時程曲線，通過式（2）可計算出白鶴灘水電站右岸地下廠房洞室群 13號機組剖面洞室群圍巖內壁關鍵部位的加速度加、卸載響應比，如圖 3所示。

圖3 13號機組剖面洞室群圍巖加速度加、卸載響應比峰值Fig.3 Acceleration loading/unloading response ratio peak value of the underground cavern groups of cross-section No.13

由加、卸載響應比理論可知，響應比值表示巖體非線性變形程度，由圖3可以看出，三大洞室關鍵部位加速度加、卸載響應比峰值變化范圍在4.05～11.52之間，其中層間錯動帶 C4在尾調室出露部位響應比峰值較大，達到11.52，這表明此部位進入非線性變形程度較大。為研究洞周圍巖的穩(wěn)定狀態(tài)，現(xiàn)選取2個典型響應比值所在部位，即主廠房拱頂和層間錯動帶 C4在尾調室上游邊墻出露部位，結合加速度、位移及應力進行綜合分析。

4.1 主廠房頂拱

圖4分別給出了主廠房拱頂部位的加速度加、卸載響應比時程曲線，加速度時程曲線、位移時程曲線及主應力差時程曲線。

由圖4可知，主廠房拱頂部位的加速度加、卸載響應比在波形上表現(xiàn)為多脈沖的形式。在地震剛發(fā)生時，加速度、位移和主應力差均未達到峰值，響應比出現(xiàn)了極值，這表明圍巖系統(tǒng)由靜止進入振動狀態(tài)后產生了巨大的響應。在隨后的振動過程中，當加速度出現(xiàn)較大幅值時（如在8、15、17.5 s），響應比均出現(xiàn)極值，同時位移及主應力差也出現(xiàn)較大幅度的波動。在地震最后6 s，加速度幅值逐漸衰減至0，位移值亦趨于0，主應力差趨于穩(wěn)定，響應比在1上下小幅波動。按照加、卸載響應比理論，當巖體失穩(wěn)破壞時，加、卸載響應比趨于無窮大，在整個地震過程中，拱頂部響應比沒有出現(xiàn)這種情況，這表明拱頂部位巖體雖發(fā)生了非線性變形，但整體是穩(wěn)定的，沒有發(fā)生失穩(wěn)破壞。另外，巖體響應位移時程曲線與加速度時程曲線有明顯的區(qū)別，可以看出，位移時程的波動頻率明顯低于加速度時程。

圖4 主廠房拱頂部位加速度加、卸載響應比Fig.4 Acceleration loading/unloading response ratio of main powerhouse top arch

4.2 尾調室上游邊墻層間錯動帶C4出露部位

圖5分別給出了尾調室上游邊墻層間錯動帶C4出露部位加速度加、卸載響應比時程曲線、加速度時程曲線、位移及主應力差時程曲線。

由圖5可以看出，尾調室上游邊墻層間錯動帶C4出露部位加速度加、卸載響應比表現(xiàn)為多脈沖形式。地震剛發(fā)生時，響應比出現(xiàn)極值，這表明此時是圍巖處于最危險狀態(tài)的時刻之一。與拱頂響應比相比，尾調室上游邊墻層間錯動帶C4出露部位響應比非線性信息明顯強烈，極值個數也明顯增加，在加速度峰值時刻出現(xiàn)最大響應。震后響應比雖趨于1，加速度幅值歸零，但產生了110 mm的永久位移，圍巖應力得到一定程度的釋放，主應力差減小了4.04 MPa，這說明錯動帶部位圍巖經受了較大的變形與應力狀態(tài)調整，但并沒有發(fā)生整體失穩(wěn)破壞。

圖5 層間錯動帶C4與尾調室上游邊墻切割部位加速度加、卸載響應比Fig.5 Acceleration loading/unloading response ratio of intersection of bedding fault zone No.C4and upstream sidewall of tailrace surge chamber

當錯動帶在洞室出露時，洞室的穩(wěn)定性受錯動帶兩側圍巖的相對變形影響較大。圖6給出了錯動帶C4與尾調室上游邊墻上、下盤的相對位移曲線。可以看出，錯動帶上、下盤巖體相互錯動產生了37 mm的相對位移，這對尾調室的邊墻的穩(wěn)定性會產生不利影響。

圖6 層間錯動帶C4與尾調室上游邊墻出露部位相對位移Fig.6 Relative displacement of intersection of bedding fault zone No.C4and upstream sidewall of tailrace surge chamber

5 結 論

（1）在本文的計算條件下，三大洞室圍巖關鍵部位加速度加、卸載響應比峰值變化范圍在4.05～11.52之間，各關鍵部位的響應比均在震后不久出現(xiàn)極值，這從理論上解釋了洞室圍巖由靜態(tài)進入動態(tài)所發(fā)生的巨大系統(tǒng)變化。

（2）主廠房拱頂部位的峰值響應比為6.12，在地震過程中，響應比均未出現(xiàn)趨于無窮大的趨勢，綜合分析應力及位移特征，主廠房拱頂雖進入非線性變形狀態(tài)，但并沒有失穩(wěn)破壞。

（3）層間錯動帶C4在尾調室上游邊墻出露部位巖體峰值響應比達到11.52，巖體非線性變形強烈，震后產生了110 mm的絕對位移，且主應力降低了4.04 MPa，表明錯動帶部位巖體經受了較大變形與應力調整，從響應比趨于1可以看出，圍巖并沒有發(fā)生整體破壞，但錯動帶上、下盤產生了37 mm的相對位移，這對尾調室的邊墻的穩(wěn)定性會產生不利影響。

本文研究方法可以應用于一般地下巖體結構地震動力穩(wěn)定分析中。如何通過數值計算方法確定地下洞室圍巖系統(tǒng)失穩(wěn)破壞的臨界加、卸載響應還有待于進一步的研究。

[1]李天斌. 汶川特大地震中山嶺隧道變形破壞特征及影響因素分析[J]. 工程地質學報,2009,16(6): 742－750.LI Tian-bin. Failure characteristics and influence factor analysis of mountain tunnels at epicenter zones of great Wenchuan earthquake[J]. Journal of Engineering Geology,2009,16(6): 742－750.

[2]左雙英,肖明. 映秀灣水電站大型地下洞室群三維非線性損傷地震響應數值分析[J]. 水力發(fā)電學報,2009,28(5): 127－133.ZUO Shuang-ying,XIAO Ming. 3D numerical analysis of nonlinear damage seismic reponse of large underground caverns of Yingxiuwan hydropower station[J]. Journal of Hydroelectric Engineering,2009,28(5): 127－133.

[3]陳健云,胡志強,林皋. 超大型地下洞室群的三維地震響應分析[J]. 巖土工程學報,2001,23(4): 494－498.CHEN Jian-yun,HU Zhi-qiang,LIN Gao. 3-D seismic response study of large scale undergroud group caverns[J].Chinese Journal of Geotechnical Engineering,2001,23(4): 494－498.

[4]李海波,劉博,呂濤,等. 一種簡單的巖體地下洞室地震安全性評價方法[J]. 巖土力學,2009,30(7): 1873－1882.LI Hai-bo,LIU Bo,Lü Tao,et al. A simple safety assessment method for underground engineering in rock under seismic loading[J]. Rock and Soil Mechanics,2009,30(7): 1873－1882.

[5]隋斌,朱維申,李曉靜. 地震荷載作用下大型地下洞室群的動態(tài)響應模擬[J]. 巖土工程學報,2008,30(12):1877－1882.SUI Bin,ZHU Wei-shen,LI Xiao-jing. Simulation on dynamic response of large underground opening complex under seismic loads[J]. Chinese Journal of Geotechnical Engineering,2008,30(12): 1877－1882.

[6]王如賓,徐衛(wèi)亞,石崇,等. 高地震烈度區(qū)巖體地下洞室動力響應分析[J]. 巖石力學與工程學報,2009,28(3):568－575.WANG Ru-bin,XU Wei-ya,SHI Chong,et al. Dynamic response analysis of rock underground caverns in highly seismic region[J]. Chinese Journal of Rock Mechanics and Engineering,2009,28(3): 568－575.

[7]苗勝軍,樊少武,蔡美峰,等. 基于加卸載響應比的載荷巖石動力學特征試驗研究[J]. 煤炭學報,2009,34(3):329－333.MIAO Sheng-jun,FAN Shao-wu,CAI Mei-feng,et al.Dynamic characteristics experiment study of loading rock based on theory of load/unload response ratio[J]. Journal of China Coal Society,2009,34(3): 329－333.

[8]王裕倉,尹祥礎,彭克銀,等. 加卸載響應比的數值模擬[J]. 地球物理學報,1999,42(5): 669－676.WANG Yu-cang,YIN Xiang-chu,PENG Ke-yin,et al.Numerical simulation on load/unload response ration theory[J]. Chinese Journal of Geophysics,1999,42(5):669－676.

[9]張暉輝,尹祥礎,梁乃剛,等. 模擬引潮力作用下巖石破壞前兆的實驗研究——加卸載響應比(LURR)理論和能量加速釋放(AER)[J]. 巖石力學與工程學報,2005,24(17): 3172－3179.ZHANG Hui-hui,YIN Xiang-chu,LIANG Nai-gang,et al.Experimental study of rock failure precursor using simulated tide stress—load/unload response ratio (LURR)and accelerating energy release(AER)[J]. Chinese Journal of Rock Mechanics and Engineering,2005,24(17): 3172－3179.

[10]石崇,徐衛(wèi)亞,周家文,等. 動力時程響應的加、卸載響應比描述研究[J]. 巖土力學,2007,28(增刊): 743－747.SHI Chong,XU Wei-ya,ZHOU Jia-wen,et al. Dynamic stability research by using load/unload ratio method[J].Rock and Soil Mechanics,2007,28(Supp.): 743－747.

[11]姜彤,馬莎,許兵,等. 邊坡在地震作用下的加、卸載響應規(guī)律研究[J]. 巖石力學與工程學報,2004,23(22):3803－3807.JIANG Tong,MA Sha,XU Bing,et al. Study of laws of load-unload response of slope under seismic load[J].Chinese Journal of Rock Mechanics and Engineering,2004,23(22): 3803－3807.

[12]姜彤,馬瑾,許兵. 基于加卸載響應比理論的邊坡動力穩(wěn)定分析方法[J]. 巖石力學與工程學報,2007,26(3):626－631.JIANG Tong,MA Jin,XU Bing. Dynamic analysis method of slope stability based on LURR theory[J].Chinese Journal of Rock Mechanics and Engineering,2007,26(3): 626－631.

[13]尹祥礎. 地震預測新途徑探索[J]. 中國地震,1987,3(1):1－7.YIN Xiang-chu. Study for a new approach to earthquake prediction[J]. Earthquake Research in China,1987,3(1): 1－7.

[14]尹祥礎,尹燦. 非線性系統(tǒng)的失穩(wěn)與前兆——響應比理論及其應用[J]. 中國科學(B輯),1991,21(5): 512－518.YIN Xiang-chu,YIN Can. Instability and precursor of nonlinear system—Response ratio theory and its application[J]. Science in China(Series B),1991,21(5):512－518.

[15]Itasca Consulting Group,Inc. FLAC(Fast Lagrangian Analysis of Continua) User's Guide,Version 5.0[M].Mineapolis: Itasca Consulting Group,Inc.,2005.

[16]張玉敏,盛謙,朱澤奇,等. 深度衰減效應對大型地下洞室群強震響應影響分析[J]. 巖土力學,2010,31(10):3197－3203.ZHANG Yu-min,SHENG Qian,ZHU Ze-qi,et al.Influence of depth attenuation on seismic response of large underground cavern group[J]. Rock and Soil Mechanics,2010,31(10): 3197－3203.