爆破振動下偏壓隧道洞口段邊坡穩定性分析

鄭明新，伍明文

爆破振動下偏壓隧道洞口段邊坡穩定性分析

鄭明新1, 2，伍明文1, 2

(1. 華東交通大學 土木建筑學院，江西 南昌 330013；2. 江西省巖土工程基礎設施安全與控制重點實驗室，江西 南昌 330013)

通過有限元軟件建立三維雙側壁偏壓隧道洞口段邊坡數值模型，提出3種爆破工況；分析3種爆破工況下邊坡的振速響應規律，通過邊坡的位移分析坡體的變形特性，依據最小二乘法對監測數據進行回歸分析得出爆破振動衰減公式，制定每個工況單段最大起爆藥量安全閾值。由解析法求出不同工況在靜、動態下坡體的穩定性變化規律。研究結果表明：隨著爆心距的增大，坡面振速衰減呈現出近快、遠慢的傳播規律；工況1相比其他工況對邊坡整體及坡腳的穩定性影響最大；通過對坡體測點各向振速監測數據進行回歸分析，得出在安全控制標準內工況1?工況3的單段最大起爆藥量分別控制在3.4，3.4及3.5 kg之內效果最好；由解析法得出，靜力狀態下開挖隧道時邊坡均處于穩定狀態。左、中導坑在爆破開挖時邊坡安全系數分別為1.13和1.18均小于1.20，極有可能誘發局部坡體滑移失穩。

爆破振動；偏壓隧道洞口段；邊坡穩定性；安全閾值；最小二乘法

爆破振動作用下隧道洞口段邊坡的穩定性及控制是隧道初始開挖過程中面臨的一個急需解決的問題。通常認為，爆破振動對隧道洞口段邊坡穩定性的影響主要體現在：爆破動荷載作用下產生的慣性力提高了坡體的下滑力，誘發其慣性失穩；同時，過大的動應力導致巖體結構面的擴展和增加，爆破損傷弱化了邊坡巖體的力學參數，使得隧道洞口段邊坡整體抗剪能力降低，導致邊坡滑移失穩[1]。因此在施工中如何控制爆破盡最大可能減輕對隧道洞口段邊坡穩定性的影響，以保證隧道安全施工順利進行就顯得極為迫切和重要。國內外相關學者對此進行了大量的研究，并取得了一系列成果。馬沖等[2]針對含軟弱夾層邊坡在爆破振動作用下坡體的振速響應規律，同時根據邊坡穩定性判據,制定出爆破振速安全閾值，提升了隧道洞口段坡體在爆破振動下的穩定性。吳燕開等[3]針對爆破荷載引起的巖石振動效應進行模擬研究，得到了巖質高邊坡在爆破振動作用下的質點振動衰減及位移變化規律，用于今后指導施工。賈黨育[4]采用數值法模擬隧道爆破開挖施工過程，分析了爆破施工對洞口段邊坡穩定性的影響規律。曹蘭柱等[5]為提高礦山邊坡安全性，通過礦山邊坡作為研究對象，量化爆破荷載對邊坡損傷程度，建立邊坡爆破振動損傷函數關系式,計算邊坡破壞時間，確定最危險區域。何理等[6]基于量綱分析理論，分析影響邊坡爆破振動速度相關物理量,推導并改進高程因素影響的爆破振動速度預測公式。王智德等[7]針對順層巖質邊坡爆破開挖，研究坡體振速的傳播規律以及的損傷程度，研究取得了爆破區任意點的振動速度與損傷深度存在對應關系。以上研究主要針對礦區爆破振動對邊坡的影響，然而在隧道爆破開挖對洞口段邊坡的動力響應分析研究還較少，同時對地處偏壓地段的隧道洞口坡體在爆破振動下的穩定性研究也較少。本文以福建漳州新建花都Ⅱ號隧道工程進口右洞作為研究對象，通過有限元數值模擬分析偏壓隧道洞口段爆破施工對坡體的影響，提出相應的控制措施及邊坡爆破振速安全閾值，旨在為相似工程的分析提供方法，這既是提升坡體穩定性的前提，也是優化爆破設計的基礎。

1 工程地質概況

某隧道工程位于福建省漳州市，隧道右洞全長 985 m，隧道最大埋深 102 m。場區表層為少量的殘坡積層(Qdl-el)，下伏基巖主要為侏羅系南園組凝灰熔巖(J3n)和燕山早期侵入的花崗巖(γ52)。進口右洞邊坡傾角為30°~35°呈現左底右高的形態，進口右洞隧道采用雙側壁導坑方式開挖，由于不規范的爆破開挖，導致隧道上方坡體發生局部滑移、開裂。

通過對現場地質的調查分析可知，影響坡體穩定性的內在因素主要為偏壓、坡度大及豐富的地下水，同時隧道走向與邊坡走向相互垂直(如圖1(a))；爆破振動為影響邊坡穩定性的外在因素，頻繁的爆破導致邊坡巖層發生錯動，產生爆破裂隙(圖１(b))。在內外因素的反復劣化作用下，導致坡體發生局部滑移、開裂。

(a) 坡體；(b) 邊坡裂縫

2 動力分析模型建立與參數選取

2.1 計算模型

根據現場實際的工程地質情況，利用MIDAS- GTS數值軟件建立數值計算模型，如圖2(a)所示。洞口段采用雙側壁方式開挖，斷面凈空尺寸為 15.27 m(寬)和10.27 m(高)。隧道支護和襯砌結構參數如圖2(b)所示；計算模型橫截面尺寸及測點分布圖見圖2(c)，測點以隧道拱頂為中點向兩側布設，布設在=10 m的斷面上；圖2(d)為坡體支護結構體系。

隧道、圍巖采用 M-C 屈服準則，圍巖及注漿體采用實體單元模擬，初支采用板單元模擬，錨桿與超前小導管采用植入式桁架模擬。模型計算范圍橫向為取20 m；向取92 m，向64 m，如圖2(a)所示。模型邊界若采用彈性邊界條件在進行動力分析時會由于波的反射作用給計算結果帶來較大誤差，所以采用1972 年Lysmer和Wass 提議的黏性邊界條件。巖土體及相關支護結構物理力學參數見表1。

(a)計算模型；(b) 隧道支護和襯砌結構參數；(c) 坡體橫截面尺寸；(d) 坡體支護結構

2.2 爆破荷載的確定

在數值計算分析中，單響藥量、加載邊界、荷載加卸時間和峰值應力等是確定爆破荷載值的主要參數。通常將爆破荷載簡化為一條具有加載與卸載過程的三角波型荷載。爆破荷載確定后，將爆破荷載以垂直壓力方式作用在開挖隧道洞周壁上[8]。通過現場施工爆破每個導坑最大的單響藥量，其爆破荷載峰值見表2。

爆破荷載峰值max的計算公式為[9]

爆破荷載達到峰值壓力所需的升壓時間為8至12 ms，卸荷時間為 40至120 ms[10]。本次加壓時間取 10 ms，卸荷時間取 100 ms，為了最佳地分析完整的爆破過程，總計算時間取 1 000 ms。

表1 物理力學參數

表2 隧道單響藥量的計算爆破荷載峰值

2.3 數值模擬方案設計

本文數值模擬主要分析各導坑在最大單響藥量工況下，爆破振動對已支護偏壓隧道邊坡在動力響應作用下對其穩定性的影響規律進行研究，工況詳情如下：

共分3種工況，工況1~工況3依次為左導坑爆破、中導坑爆破及右導坑爆破。每次爆破開挖進尺為2 m ，鑒于施工現場中導坑開挖進尺累積達到10 m(第5個爆破開挖循環)時，拱頂附近土體發生如圖1所示的坡體局部開裂，所以每個工況均以第5個爆破開挖循環作為本次研究對象。

3 爆破振動作用下坡體動力響應 規律

偏壓隧道洞口段邊坡在爆破振動作用下，其變形特性可以通過各測點振動速度、各向位移值及應力應變增量來顯示坡體在爆破振動下的變形及穩定情況。

3.1 測點振度響應規律

由圖2(a)及圖3可知：1) 隨著爆心距的增大，測點速度峰值逐漸減小，呈現出薩道夫斯基傳播規律、衰減規律遠慢近快，其中和向振速峰值變化趨勢相比向更大；2) 在工況2中，測點5主振速峰值在所有工況中最大，為22.16 cm/s小于一般巖質邊坡質點振動安全閾值速度25 cm/s[11]，右導坑巖性較佳且覆土相比其他工況更厚，所以工況3所有測點各個方向振速均小于其他2種工況在各向的振速值；3) 在所有工況中，每個工況的和向峰值振速總體上大于向峰值振動速度，所以在同類工程中可只監測向與向振動速度達到簡化爆破振動監測的目的；4) 工況1在坡腳(測點1)處的合成峰值振速均大于其他工況在坡腳的合成峰值振速，可見左導坑爆破開挖對坡腳穩定性影響較大，需密切關注坡體變形情況，必要時可降低爆破藥量，分多次弱爆方可提高邊坡的穩定性。

3.2 邊坡測點位移變化規律

爆破振動結束后，所有工況在爆源處周圍巖體產生整體向坡外的變形，如圖4所示。

由圖5可知：1) 工況1最大總位移為5.83 mm，是工況2的1.25倍和工況3的7.29倍，在所有工況中最大；2) 工況1 在坡腳處測點1的總位移最大為2.27 mm。左導坑爆破對整個偏壓隧道邊坡的穩定性影響最大，對坡腳的擾動是所有工況中之最；坡腳擾動后在偏壓作用下，將可能誘發整個坡體下滑最終導致邊坡失穩。3) 由圖5(a)及5(b)可知：中導坑爆破振動對拱頂上方土體擾動最強，導致測點7發生下沉并右移，所以中導坑爆破振動極有可能誘發坡頂發生局部滑移、開裂；4)鑒于工況3的坡體位移最小，所以其對坡體的影響最小。

綜上可得，左、中導坑爆破開挖對坡體的穩定性影響較大，必要時需控制爆破藥量、分多次弱爆，方可達到減振護坡的效果。

3.3 邊坡爆破振動衰減公式的回歸分析

隧道爆破后會急劇釋放能量，并以波動的形式向外傳播，引起坡體內巖土體介質的快速振動，然而爆破振動衰減快。依據現有的研究成果，當前國內外大多采用薩道夫斯基經驗衰減公式對爆破振動數據進行回歸分析。

(a) x向振速；(b) y向振速；(c) z向振速；(d) 合成振速

(a) 工況1；(b) 工況2；(c) 工況3

(a) y向位移；(b) z向位移；(c) 總位移

式中：為爆破產生的質點振動速度，cm/s；為爆源與測點的直線距離，m；為爆破單段最大起爆藥量，kg；為與地質條件、圍巖特性有關的系數；為爆破衰減系數。

對式(1)兩邊取對數得

令ln=，ln=，ln(1/3/)=，=，則將式(2)變換得到一元線性回歸公式如下：

利用式(1)~式(3)對表3數據進行最小二乘法回歸，得薩道夫斯基經驗衰減公式(表(4))。(注：由前面研究結果可知，工況1和2對坡體的穩定性影響較大，同時在2種工況下向與向振速比向更大，所以表3僅考慮工況1和2下向和向及合成振速值。)

表3 隧道邊坡爆破振動測點數據

由表4中的薩道夫斯基經驗衰減公式可以推算，不同導坑在爆破中的振速峰值在安全控制標準以內的單段最大起爆藥量。

表4 薩道夫斯基經驗衰減公式

注：1Y向指工況1方向振速，其余等同。

由工程實際及安全控制標準可得知：min=7.8 m，max=25 cm/s，將上述公式代入表4進行計算，得工況1及工況2單段最大起爆藥量均為3.4 kg，鑒于工況3采用現場設計爆破藥量對坡體穩定性影響小同時也能滿足施工作業需要，所以其最小爆破藥量可為3.5 kg。

3.4 邊坡穩定性分析

通過對坡體在不同工況下所對應靜動態下的穩定性安全系數分析，可得爆破振動下的邊坡相比靜態下的邊坡在安全性變化情況。

3.4.1 靜力狀態下邊坡穩定性分析

通過現場局部坡體滑移動態可知；滑移體主要為表層坡積土的滑動為主，所以圖6中，邊坡土體的物理力學參數以坡積土為主。其他幾何參數為，臺階坡腳=30°，臺階高度=29 m。

圖6 邊坡穩定性計算簡圖

依據邊坡穩定性系數的定義與力的平衡條件可得靜力邊坡穩定性系數為：

若直接利用式(4)求的最小值相對困難，為了求解的便利，將三角函數用坡高與水平距離來表示，由圖6的幾何關系可得。

將式(5)代入式(4)，得：

為求穩定性系數的最小值，對式(6)導數等于0，得：

化簡可得：

從圖6可以得出，當≦cot時，滑面不在坡內，因此，通過方程(5)可得解為：

由式(9)中的代入式(5)和式(6)可分別求得邊坡在靜力狀態下的最小穩定系數及對應穩定性系數下的最小滑移面及水平方向的夾角分別為min=1.25，min=21.74°。

為了驗證理論計算的準確性，在同等條件的基礎上建立三維隧道邊坡數值模型，通過強度折減法分別模擬左、中及右導坑在靜力開挖下的坡體穩定性安全系數分別為：1.22，1.24及1.33，其與理論計算值相近。

3.4.2 爆破動力狀態下邊坡穩定性分析

上一節分析了靜力作用下邊坡的穩定性，其最小穩定性系數=1.22>1.20[12](一級公路邊坡穩定性安全系數為1.20)，所以在靜力狀態下邊坡是穩定的。但邊坡在承受爆破振動力作用時，坡體中的潛在滑體會附加一慣性力，使潛在滑體的下滑力增大，從而降低了邊坡的穩定性，提高了邊坡失穩坍塌的幾率。

爆破振動引起的作用于邊坡潛在滑體并指向爆破方向的等效靜力，為：

式中：k為等效靜力系數；為滑體自重；為爆破動力系數，其值通常在0.1~0.2之間，本次取0.15；為地震系數是巖體質點振動加速度與重力加速度(9.81 m/s2)的比值。

將3種工況坡腳處的加速度代入式(11)可得工況1~工況3因爆破振動引起的等效水平靜力分別為滑體重力的6.92%，4.60%及2.73%，且令其為。

所以，在爆破振動下邊坡的穩定性系數如下：

將不同工況的值與式(5)分別代入式(12)解得不同工況下最小穩定性系數分別為：1.13，1.18及1.26。

由表5與圖7可知：工況1和2的邊坡穩定性安全系數均小于1.20，所以左、中導坑爆破開挖極易誘發局部坡體滑移；其中工況1的安全系數相比工況2小4.24%；在邊坡穩定性損失系數中，工況1穩定性損失值在所有工況中最高，分別是工況2的1.52倍和工況3的1.40倍。綜上可得，左導坑爆破開挖對坡體穩定性影響最大，在必要時可降低爆破藥量、分多次弱爆方可提高坡體的穩定性。

表5 邊坡靜/動態安全系數值

圖7 不同工況系下邊坡安全系數曲線

3.5 邊坡振速監測數據分析

采用爆破測振儀對各導坑爆破振動對上方邊坡振速進行監測，每個測點安裝一個振速傳感器，測點分布見圖2(c)所示。各導坑爆破振動時對上方邊坡測點的各向振速如表6所示。

由表3與表6可知：邊坡測點的實測振速略小于數值計算的振速，出現偏差說明現場地質及巖層分布復雜性等相關因素所致，但二者在數值上較為相近，并且各向振速的大小及分布位置相接近；最后對不同工況的各向振速監測數據進行了回歸分析可知：工況1和2的單段最大起爆藥量均為3.4 kg，其結果與數值仿真結果相近。

表6 邊坡振速監測數據

注：為便于與表3對比，因此表6中的監測數據為工況1和2下的向和向及合成振速。

4 結論

1) 隧道邊坡爆破振速衰減規律為，隨著爆心距增大，振速隨之減小，越靠近爆心衰減幅度越大；左導坑爆破對坡腳的動力響應最大，中導坑合成振速最大為22.18 cm/s。

2) 工況1對偏壓隧道洞口段邊坡位移變形影響最大同時其也對坡腳的穩定性影響最強，總位移值達2.27 mm；工況2對拱頂上方土體擾動最大，致使測點7發生下沉并右移，極有可能誘發坡頂處發生局部滑移失穩。

3) 根據最小二乘法對坡體測點各向振速數據與進行回歸分析，得出向、向及矢量合成的薩道夫斯基經驗衰減公式，得出在安全控制標準內工況1(左導坑)、工況2(中導坑)及工況3(右導坑)的單段最大起爆藥量分別控制在 3.4，3.4及3.5 kg 之內效果較好，為后期爆破開挖提供了參考與借鑒。

4) 由解析法分析可知，在靜力作用下開挖左、中及右導坑時所對應的邊坡穩定性系數均大于1.20，故隧道在靜力開挖狀態下坡體是穩定的；但在爆破振動作用下開挖左及中導坑時，邊坡穩定性系數分別為1.13和1.18均小于1.20極易誘發邊坡局部滑移失穩，同時左導坑在爆破作用下穩定性損失系數最大為7.37%。

5) 通過對坡體測點振速監測值與數值仿真值對比與回歸分析可知：兩者數值結果、振速大小與分布位置、各導坑單段最大起爆藥量大體相近可為后期同類工程的分析提供參考。

[1] 鄒玉君, 嚴鵬, 劉琳, 等. 白鶴灘水電站壩肩邊坡爆破振動對周邊民房影響評價及控制[J]. 振動與沖擊, 2018, 37(1): 248?258. ZOU Yujun, YAN Peng, LIU Lin, et al. Assessment and control of the blasting vibration effect on local residences during the excavation of abutment slope in the Baihetan project[J]. Journal of Vibration and Shock, 2018, 37(1): 248?258.

[2] 馬沖, 詹紅兵, 姚文敏, 等. 爆破振動作用下含軟弱夾層邊坡穩定性及安全判據[J]. 爆炸與沖擊, 2018, 38(3): 563?571. MA Chong, ZHAN Hongbing, YAO Wenmin, et al. Stability and safety criterion of a slope with weak interlayer under blasting vibration[J]. Explosion and Shock Wave, 2018, 38(3): 563?571.

[3] 吳燕開, 胡曉士, 石玉斌, 等. 爆破開挖對巖質高邊坡穩定性的影響分析[J]. 爆破, 2017, 34(4): 66?72. WU Yankai, HU Xiaoshi, SHI Yubin, et al. Influence of blasting excavation on high rock slope stability[J]. Blasting, 2017, 34(4): 66?72.

[4] 賈黨育. 隧道爆破開挖對洞口段邊坡穩定性影響的數值分析[J]. 礦冶工程, 2017, 37(1): 25?28. JIA Dangyu, Numerical analysis for effect of tunnel excavation by blasting on slope stability of portal section [J]. Ming and Metallurgical Engineering, 2017, 37(1): 25?28.

[5] 曹蘭柱, 王珍, 王東. 邊坡受爆破振動損傷及其穩定性分析[J]. 中國安全科學學報, 2018, 28(2): 75?80. CAO Lanzhu, WANG Zhen, WANG Dong. Law of blasting vibration damage and stability analysis to slope in open-pit mine[J]. China Safety Science Journal, 2018, 28(2): 75?80.

[6] 何理, 鐘冬望, 陳晨, 等. 巖質高邊坡開挖施工的爆破振動監測與分析[J]. 金屬礦山, 2017(1): 6?10. HE Li, ZHONG Dongwang, CHEN Chen, et al. Monitoring and analysis of blasting vibration in high rocky slope excavation[J]. Metal Mine, 2017(1): 6?10.

[7] 王智德, 夏元友, 周雄, 等. 順層巖質邊坡爆破的振動控制及損傷特性[J]. 爆炸與沖擊, 2017, 37(1): 27?36. WANG Zhide, XIA Yuanyou, ZHOU Xiong, et al. Blasting vibration control and damage characteristics rock slope[J]. Explosion and Shock Wave, 2017, 37(1): 27?36.

[8] 夏一鳴. 鄰近隧道爆破施工對既有高鐵隧道的振動安全分析[D]. 南昌: 華東交通大學, 2017. XIA Yiming. Analysis of vibration impact caused by adjacent tunnel blasting on existing high-speed railway tunnel[D]. Nanchang: East China Jiaotong University, 2017.

[9] 賈磊, 解詠平, 李慎奎. 爆破振動對鄰近隧道襯砌安全的數值模擬分析[J]. 振動與沖擊, 2015, 34(11): 173? 177, 211. JIA Lei, XIE Yongping, LI Shenkui. Numerical simulation for impact of blasting vibration on nearby tunnel lining safety[J]. Journal of Vibration and Shock, 2015, 34(11): 173?177, 211.

[10] 包德勇.小間距隧道爆破振動的影響分析[J].鐵道標準設計, 2008(11): 89?92.BAO Deyong. Analysis on the influence of blasting vibration in small space tunnel[J]. Railway Standard Design, 2008(11): 89?92.

[11] 歐陽建華, 李雷斌, 吳亮, 等. 巖質凸形邊坡體爆破質點振動監測與分析[J]. 爆破, 2015, 32(1): 54?56, 60. OUYANG Jianhua, LI Leibin, WU Liang, et al. Monitoring and analysis of particle vibration of convex rock slope under blasting loading[J]. Blasting, 2015, 32(1): 54?56, 60.

[12] 熊文林. 山區高速公路隧道洞口邊坡穩定性的數值模擬[J]. 安全與環境工程, 2013, 20(5): 108?112. XIONG Wenlin. Numerical simulation analysis of the slope stability of highway tunnel in mountainous area[J]. Safety and Environmental Engineering, 2013, 20(5): 108?112.

Stability analysis of slope at biased tunnel portal section under blasting vibration

ZHENG Mingxin1, 2, WU Mingwen1, 2

(1. School of Civil and Architectural Engineering, East China Jiaotong University, Nanchang 330013, China; 2. Key Laboratory of Geotechnical Engineering Infrastructure and Security Control, East China Jiaotong University, Nanchang 330013, China)

A 3D numerical model of the slope at biased tunnel portal Section with double side drift method was established by using finite element software. Firstly, the vibration velocity response law of slope under three blasting conditions was analyzed, secondly the deformation characteristics were analyzed by the displacement of the slope, Then based on the least squares method analysis was carried out on the monitoring data, it was concluded that blasting vibration attenuation formula for working condition of each single dose period of maximum detonating safety threshold. Finally, the stability change law of the static and dynamic slope body under different working conditions was obtained by analytical method. the study indicates that with the increase of blaster center distance value, velocity attenuation slope presents a fast, nearly far slow propagation law, elevation amplification effect; Compared with other working conditions, condition 1 has the greatest influence on the stability of the whole slope and the slope foot. By means of the regression analysis of the monitoring data of the velocity of each direction vibration of the slope to get attenuation formula under different working conditions, the formula is concluded thatin safety control standard condition 1-conditions 3 single segment of initiation largest doses respectively control within 3.4 kg, 3.4 kg and 3.5 kg effect is best, Finally, it can be obtained by analytic method that the slope is in a stable state when the tunnel is excavated in the static state, while the safety coefficient of the slope stability is 1.13 and 1.18, respectively, which are less than 1.20 in the blasting excavation of the left and middle pilot tunnel, which may induce local slope sliding instability.

blasting vibration; biased tunnel portal section; slope stability; safety threshold; least squares method

U455

A

1672 ? 7029(2019)08? 2018 ? 10

10.19713/j.cnki.43?1423/u.2019.08.019

2018?10?22

國家自然科學基金資助項目(51568022)；江西省自然科學基金資助項目(20171BAB206056)

鄭明新(1966?)，男，江西南昌人，教授，博士，從事滑坡與路基、隧道工程的教學和研究工作；E?mail：zhmx@ecjtu.jx.cn

(編輯 蔣學東)