賈澄澄，卓 莉，肖明礫，謝紅強，何江達

(四川大學 水力學與山區(qū)河流開發(fā)保護國家重點實驗室，水利水電學院，四川 成都 610065)

某水電站位于甘孜州丹巴縣境內，正常蓄水位2 020.00 m，死水位2 018.50 m，正常蓄水位以下庫容9.58 萬m3，調節(jié)庫容4.5 萬m3，水電站裝機容量90 MW，額定水頭157 m，設計引用流量66 m3/s。水電站右岸有壓隧洞引水，隧洞長11.173 km，廠址位于丹巴大橋下游約2.7 km的大渡河右岸漫灘上。丹巴滑坡位于大金川右岸丹巴縣城白呷山高陡斜坡下部，是在古滑坡的基礎上逐漸發(fā)育形成并向古滑坡后部左側擴張。其高程為1 887～2 199 m，前后緣高程差約為223 m，滑坡體寬約為250～280 m，縱向長度約為290 m，滑坡體平均厚約30 m，為一特大型堆積體滑坡。

鑒于水電站引水隧洞距離上覆丹巴滑坡體距離較近，爆破施工過程中不可避免對周圍建筑的安全造成不利的影響，會產(chǎn)生較大爆破振動效應，對上覆滑坡體造成動力擾動[1,2]，危及上覆滑坡體的穩(wěn)定性，誘發(fā)滑坡失穩(wěn)。因此，本文結合水電站區(qū)域地形地質條件、隧洞設計參數(shù)及施工爆破參數(shù)等，建立有限元模型，計算施工爆破振動對滑坡體的動力擾動效應，分析爆破作用下振動位移、速度、加速度特性以及振動衰減規(guī)律，基于GeoStudio有限元軟件，采用剛體極限平衡法和有限元法兩種方法計算滑坡穩(wěn)定安全系數(shù)，對上覆滑坡體穩(wěn)定性進行評價。

1 引水隧洞-滑坡體有限元計算

1.1 有限元計算模型

滑坡與隧洞整體有限元模型見圖1，隧洞局部模型見圖2。有限元計算選取引水隧洞距丹巴滑坡垂直距離最近剖面為研究對象。計算模型洞中心坐標為(2 205.7，1 938)，底面取至高程1 499 m，向上取至坡頂，鉛直向高度1 344.6 m，模型左側邊界至大渡河河床中心，右側邊界距離隧洞中心線331.3 m，水平長度共計2 085.3 m。有限元網(wǎng)格采用4節(jié)點等參單元對滑坡進行離散，單元總數(shù)共1 935 個，節(jié)點總數(shù)1 930 個。二維有限元的計算坐標系選定為：X軸為水平方向，由滑坡指向洞；Y軸為垂直于X軸且鉛直向上。

圖1 水電站丹巴滑坡和隧洞整體有限元網(wǎng)格Fig.1 Finite element mesh of Danba landslide and tunnel

圖2 水電站引水隧洞局部有限元網(wǎng)格Fig.2 Local finite element mesh of diversion tunnel of hydropower station

1.2 有限元計算參數(shù)及邊界條件

有限元模型左右及底邊界面采取黏彈性人工邊界，以便計算區(qū)域內的外行波在截斷邊界處不會反射回來。黏彈性邊界[3,4]推導采用坐標系見圖3。

平面柱面波的運動方程為：

(1)

cs=(G/ρ)1/2

(2)

式中：v為介質的切向位移；cs為剪切波在連續(xù)介質中的傳播速度；G為連續(xù)介質的剪切模量；ρ為介質的質量密度。

圖3 黏彈性邊界坐標系示意圖Fig.3 Diagram of viscoelastic boundary coordinate system

根據(jù)類似工程的原理，在邊界上簡單地施加黏性邊界，將其中的正應力σ和正應變ε相應地換為剪應力τ和剪應變γ，位于邊界rb上的剪應力為：

(3)

任一半徑rb處的應力與該處的速度和位移的關系可以表示如下：

(4)

式(4)可以寫成如下式：

(5)

式(5)等價于一阻尼系數(shù)為ρcs的阻尼器并聯(lián)上一個剛度系數(shù)為G/2rb的線性彈簧，可以精確地模擬滿足式(3)的剪切波由有限域向無限域的傳播。

各滲透分區(qū)巖體的各向同性滲透系數(shù)見表1。

表1 巖土體及襯砌物理力學參數(shù)表Tab.1 Physical and mechanical parameters of rock and soil

2 爆破方案

有限元計算中施工爆破參數(shù)見表2，布孔及網(wǎng)路圖見圖4。

根據(jù)水電站引水隧洞沿線的地質條件，引水隧洞距丹巴滑坡垂直距離最近洞段圍巖為Ⅲ～Ⅳ類，設計開挖斷面為馬蹄形，面積約41.58 m2，開挖爆破每次循環(huán)進尺最大允許值為3 m。

表2 引水隧洞圍巖爆破設計參數(shù)統(tǒng)計表Tab.2 Statistics of blasting design parameters ofsurrounding rock of diversion tunnel

圖4 引水隧洞圍巖爆破網(wǎng)絡圖Fig.4 Blasting network diagram of surrounding rock of diversion tunnel

表3 引水隧洞圍巖爆破荷載參數(shù)Tab.3 Blasting parameters of surrounding rock of diversion tunnel

引水隧洞圍巖爆破荷載參數(shù)見表3，等效動荷載時程三角形分布見圖5。炸藥采用2號巖石乳化炸藥，爆速為3 500 m/s，炸藥密度為1 000 kg/m3，根據(jù)爆破參數(shù)和爆破網(wǎng)絡圖，采用施加等效荷載[5]的方法，得到不同雷管段別的爆破峰值壓力 、等效動荷載 及其時程三角形分布。巖體的阻尼比取0.03，覆蓋層巖石的阻尼比取0.05，爆破荷載壓力上升時間為100 μs，正壓作用時間為100 μs，主掏槽孔起爆不考慮。

圖5 微差爆破等效動荷載時程三角形分布圖Fig.5 Triangular distribution of equivalent dynamic load time history of millisecond blasting

3 引水隧洞施工爆破振動效應

選擇邊坡距爆源同高程不同水平距離的7個特征點(P1～P7)、爆源正上方不同高程的5個特征點(S1～S5)研究爆破振動特性，特征點分布位置示意圖見圖6。

圖6 特征點分布位置示意圖Fig.6 Map of feature points distribution

表4給出了各特征點的振動位移、速度、加速度峰值結果。限于篇幅僅給出了部分特征點的振動位移、速度、加速度時程曲線圖7-10以及同高程特征點振速峰值與爆心距的衰減曲線圖11。

在隧洞爆破荷載作用下，特征點都出現(xiàn)了不同頻率、振幅的振動現(xiàn)象。離爆源的距離越近，質點的波動特性越強，質點的振動頻率和振幅越大。當爆心距增大到一定程度后，質點的位移、速度及加速度的變化小到可以忽略。從表4中可以看出，P1測點處出現(xiàn)水平向振動最大位移為0.009 cm，水平向最大振動速度為10.8 cm/s，水平向最大振動加速度為44 000 cm/s2；S1測點處出現(xiàn)鉛直向振動最大位移為0.012 cm，鉛直向最大振動速度為15 cm/s；鉛直向最大振動加速度為58 000 cm/s2。

表4 各特征點的振動位移、速度、加速度峰值統(tǒng)計表Tab.4 Statistics of vibration displacement, velocity and acceleration

從圖7、8可以看出，隧洞圍巖巖體對施工爆破振動應力波具有明顯的衰減作用，質點振動隨著爆源距的增大，起振時間越晚，振幅衰減率越慢。P1特征點幾乎在爆破同時開始振動，而且位移振幅在很短的時間內衰減至0附近，P5特征點振動時間較起爆時間延遲約0.15 s，振動位移在1 s內衰減了極值的90%；從圖9、10來看，S5特征點振動時間較起爆時間延遲約0.12 s，在振動2 s后，其豎直向速度和加速度均衰減至極值的50%。

圖7 特征點P1的振動位移時程曲線圖Fig.7 Time history curve of vibration displacement of point P1

圖8 特征點P5的振動位移時程曲線圖Fig.8 Time history curve of vibration displacement of point P5

圖9 特征點S5的振速時程曲線圖Fig.9 Time history curve of vibration velocity of point S5

圖10 特征點S5的振動加速度時程曲線圖Fig.10 Time history curve of vibration acceleration of point S5

從特征點振速峰值衰減曲線可以看出，振速峰值在距爆心距15 m以內衰減速率最快，振速峰值與爆心距之間近似呈線性關系，在爆心距15 m以后，振速峰值與爆心距之間呈非線性關系，衰減速率逐漸變緩。在爆心距100 m范圍外，振速峰值小于0.2 cm/s，表明在距爆源100 m范圍外基本不受爆破影響。

圖11 同高程特征點振速峰值與爆心距的衰減曲線圖Fig.11 Attenuation curve of the peak velocity of the same point and the distance of the center of detonation

根據(jù)國家標準[6]規(guī)定的爆破振動安全允許標準，選擇最不利安全允許質點振速0.15 cm/s作為丹巴滑坡體的爆破振動控制標準。從圖11可以看出：當爆心距為120 m時，質點的振速峰值為0.15 cm/s，離引水隧洞最近的滑坡體質點(距離1 420 m)振速極值僅為0.024 cm/s，遠小于標準下限值0.15 cm/s，表明引水隧洞的施工爆破不會對丹巴滑坡體安全造成影響。

4 滑坡穩(wěn)定性分析

基于有限元軟件GeoStudio，分別采用剛體極限平衡法和有限元法，計算施工爆破前后丹巴滑坡的穩(wěn)定安全系數(shù)，滑坡的施工爆破擾動作用效應采用擬靜力法計算，且考慮最不利情況，以附加水平最大加速度1.35 cm/s2作為振動計算加速度，同時根據(jù)抗震設計規(guī)范[7]考慮振動效應的折減系數(shù)和質點的動態(tài)分布系數(shù)計算滑坡體質點的動荷載慣性力。兩種計算方法所得爆破前后的滑坡體穩(wěn)定性計算成果見表5。

表5 爆破前后滑坡體穩(wěn)定性計算成果表Tab.5 Calculation results of landslide stability beforeand after blasting

可以看出：采用剛體極限平衡法時，在無爆破振動作用下，滑坡體的穩(wěn)定安全系數(shù)為1.358；在爆破作用下，滑坡受到微弱的振動作用，其安全系數(shù)為1.336，較無爆破作用時略有降低，但均滿足規(guī)范要求，滑坡處于穩(wěn)定狀態(tài)。有限元法計算結果與剛體極限平衡法相似，在無爆破振動作用下，滑坡體的穩(wěn)定系數(shù)為1.381；在爆破作用下，其安全系數(shù)為1.353，均滿足規(guī)范要求。有限元法結果均大于剛體極限平衡法結果，這主要是由于兩種計算方法的差異，以及滑弧上條分位置與應力積分點位置不同造成滑弧位置上的差異所致。在引水隧洞爆破作用下，丹巴滑坡雖然受到微弱的振動作用，但不會造成丹巴滑坡體的失穩(wěn)。

5 結 論

(1)離爆源的距離越近，質點的振動頻率和振幅越大，且振動衰減也較快，當距離爆心距足夠大時，水平向的振動幅值和鉛直向值差不多，振動位移、速度、加速度的時程曲線圖相似，此時爆破產(chǎn)生的振動影響很小。

(2)振速峰值在距爆心距15 m以內衰減速率最快，在爆心距15 m以后衰減速率逐漸變緩，在距爆源100 m范圍外基本不受爆破影響。

(3)剛體極限平衡法和有限元法計算得到的滑坡穩(wěn)定安全系數(shù)均滿足規(guī)范要求，表明引水隧洞爆破不會影響丹巴滑坡體的穩(wěn)定和安全。

□

[1] 謝承煜,羅周全,賈楠,熊立新,程貴海. 露天爆破振動對臨近建筑的動力響應及降振措施研究[J]. 振動與沖擊,2013, 2013,32(13):187-193.

[2] 劉超英,葛雙成. 水庫輸水隧洞施工爆破振動對大壩安全影響監(jiān)測及分析[J]. 地震工程學報,2015,37(2):367-371.

[3] 朱彬. 城市地下工程結構抗震分析研究[D]. 西安：西安科技大學,2005.

[4] DEEKS A J, RANDOLPH M F. Axisymmetric time-domain transmitting boundaries[J]. Journal of Engineering Mechanics, 1994,120(1):25-42.

[5] 張玉成,楊光華,劉鵬,姚捷,姜燕. 爆破荷載在數(shù)值計算中的等效施加方法研究[J]. 地下空間與工程學報,2012,8(1):56-64.

[6] GB6722-2003，爆破安全規(guī)程[S].

[7] SL 203-1997，水工建筑物抗震設計規(guī)范[S].