吳 亮李 鳳盧文波陳 明許 鋒

(1.武漢科技大學理學院,湖北武漢430065; 2.中鐵港航-武漢科技大學爆破技術研究中心,湖北武漢430065; 3.武漢大學水資源與水電工程科學國家重點實驗室,湖北武漢430072; 4.貴州新聯(lián)爆破工程集團有限公司,貴州貴陽550002)

爆破擾動下鄰近層狀圍巖隧道的穩(wěn)定性與振速閾值*

吳 亮1,2,李 鳳1,2,4,盧文波3,陳 明3,許 鋒1,2

(1.武漢科技大學理學院,湖北武漢430065; 2.中鐵港航-武漢科技大學爆破技術研究中心,湖北武漢430065; 3.武漢大學水資源與水電工程科學國家重點實驗室,湖北武漢430072; 4.貴州新聯(lián)爆破工程集團有限公司,貴州貴陽550002)

爆炸荷載下鄰近既有層狀圍巖隧道的迎爆側易發(fā)生開裂、破壞,確保既有隧道圍巖在爆破擾動下的安全是隧道近接施工中的關鍵問題?；趶椥粤W和結構力學原理,建立了既有層狀圍巖隧道迎爆側最危險點的簡化力學模型,得到了最危險點的應力計算表達式,分析了隧道半徑與巖層厚度的比值、巖層的傾角、巖層所對應的圓心角和隧道埋深等因素對既有隧道圍巖穩(wěn)定性的影響規(guī)律,并依據(jù)巖體最大拉應力判據(jù)確定了既有層狀隧道圍巖穩(wěn)定的質點臨界振動速度,結合地震波傳播時的衰減規(guī)律,可以得出爆破施工時的最大單響藥量,從而為隧道近接爆破動力擾動施工提供理論依據(jù)。

圍巖穩(wěn)定;爆破;隧道;層狀圍巖;振動速度

隨著我國國民經(jīng)濟的迅速發(fā)展,大量的鐵路單線需要修建復線,如建設中的寶成復線、株六復線等。復線隧道施工爆破時產(chǎn)生的振動對既有隧道的影響一般都比較大,特別是對于復雜地質(如層狀地質)的情況,有可能使鄰近既有隧道損傷甚至坍塌。另外,地下開采、公路隧道以及水電站硐室群等工程都涉及到爆破對既有層狀隧道擾動的問題。自然界中具有層狀結構的沉積巖占陸地面積的三分之二,許多變質巖也具有層狀構造特征,因而層狀巖體是地下工程中經(jīng)常遇到的一種巖體,研究復雜地質條件下既有隧道圍巖損傷破壞的影響因素意義重大。

長期以來,隧道圍巖穩(wěn)定性判定沒有科學合理的方法,常以洞周位移或塑性區(qū)經(jīng)驗值作為穩(wěn)定判定指標[1]。由于洞周位移與隧道的振動存在一定的聯(lián)系,故現(xiàn)實工程中一般通過監(jiān)測洞周壁的振動速度監(jiān)控隧道圍巖的穩(wěn)定性。許多學者對鄰近均質隧道圍巖的穩(wěn)定性開展了研究,胡文清等[2]最早將有限元強度折減法引入隧道圍巖穩(wěn)定性分析中,并對木寨嶺隧道進行了分析;文獻[3-7]中對均質隧道進行了系統(tǒng)研究,提出了隧道的破壞機理及剪切安全系數(shù)和抗拉安全系數(shù),并建立了相應的計算方法。邊克信[8]建議采用一維應力波理論,根據(jù)巖石抗拉強度或巖石拉伸極限應變推算臨界振動速度;朱瑞賡等[9]提出以巖體動力強度與隧洞所受動、靜應力之和相平衡來確定爆破振動影響下無襯砌隧洞出現(xiàn)不同破壞時隧洞圍巖臨界振動速度的計算方法;P.K.Singh[10]對爆破沖擊荷載作用下臨近煤礦巷道的動力損傷問題進行了數(shù)值分析;J.Torano等[11]綜合考慮了不同介質、裝藥量、起爆順序等因素,通過建立三維有限元模型研究了隧道爆破振動速度的分布規(guī)律;劉國華等[12]基于加權雙剪強度準則,提出了鄰近爆破的安全評估方法。另外,國內學者對近距離爆破下隧道周邊振動場分布規(guī)律也進行了大量的研究[13-19],分析了不同爆心距、地質條件、間距、布置方位等情況下既有隧道圍巖與襯砌的響應規(guī)律。

目前,國內外學者已經(jīng)對新建隧道爆破擾動下既有隧道的圍巖穩(wěn)定性進行了大量的分析,但是定量的理論研究相對而言還比較少。本文中,基于彈性力學和結構力學原理建立既有層狀圍巖隧道迎爆側的簡化模型,研究新建隧道爆破擾動下既有層狀圍巖隧道的穩(wěn)定性影響因素及變化規(guī)律,得出既有層狀圍巖隧道的臨界振動速度,以期為工程施工提供理論依據(jù)。

1 計算模型及過程

1.1 計算模型及假設

在層狀地質條件下鄰近隧道開挖過程中,爆破擾動可能會引起鄰近既有層狀圍巖隧道的振動,繼而危及既有隧道的安全,特別是既有隧道迎爆側的安全。隧道層狀圍巖受力情況如圖1所示,其中巖層與水平方向的夾角為θ,與隧道相交部分在隧道中可見長度為2L,對應的圓周角為δ,巖層厚度為h,隧道半徑為R,鄰近隧道洞壁的層狀圍巖受到覆蓋圍巖的重力為G2,層狀巖體自重為G1,鄰近隧道爆破開挖而產(chǎn)生的爆破擾動荷載為Fd。

根據(jù)層狀圍巖的實際受力及位移情況,對其作如下假設:(1)該層狀圍巖可以等效為兩端固定的梁;(2)巖體為均質各向同性材質;(3)梁的破壞形式為拉伸破壞;(4)爆破荷載以水平入射;(5)梁的重力均勻分布在梁上,為常體力,作用于梁上的力可分解為平行于梁的力和垂直于梁的力。

1.2 計算過程

梁(層狀圍巖)的受力情況如圖2和圖3所示,其中x、y分別為平行于梁和垂直于梁的方向;σd和τd分別為爆破荷載與層狀巖體上覆蓋圍巖重力的合力作用于梁垂直方向的正應力和平行方向的切應力;fx和fy為梁自重在x、y方向上的分力,即

圖1 隧道層狀圍巖受力示意圖Fig.1 Force diagram of layered rock surrounding the tunnel

圖2 梁的受力情況Fig.2 Stress analysis of beam

圖3 解除右側約束后梁的受力Fig.3 Stress after removing constraint of right side

由圖2可知,當梁僅受水平分力作用時,運用半逆解法[20]可假設梁沿y方向的應力為:

引入應力函數(shù)Φ,有:

式中:σx為梁沿x方向的應力,τxy為沿y軸方向作用于梁x面上的切應力。將式(1)代入式(2)得:

式中:f1(y)和f2(y)表示任意關于y的一次函數(shù)。將式(3)代入相容方程得:

由于式(4)對于任意的x都成立,故解得:

式中:A1、B1、C1、D1和E1為常數(shù)。將其代入式(2)得:

將結構一端的約束解除代之以相應的力,得如圖3所示的結構受力。

根據(jù)力法原理[21]解得:

式中:FN為拉力,FS為剪切力,M為力矩。引入邊界條件:解得:

當梁僅受垂直方向的力時,則由圖2中梁的受力情況可知:

重復上述方法,引入應力函數(shù)后求解得到:

式中:A2、B2、C2、D2、E2、F2、G、H2和K2為常數(shù)。同理,將結構一端的約束解除代之以相應的力,得到梁的受力情況如圖3所示。

利用力法原理解得:

因為yz平面是梁和荷載的對稱面,所以應力分布應當對稱于yz面;σx和σy是x的偶函數(shù),τxy是x的奇函數(shù),綜合以上條件可知:E2=F2=G=0。引入邊界條件:(σy)y=h/2=-σd,(σy)y=-h/2=0,解得:

綜上所述,梁的中間截面上最危險點x方向拉應力的表達式為:

利用圖1中參數(shù)可以將σd、τd、fx、fy表示如下:σd=Fdsinθ+G2cosθ;τd=G2sinθ-Fdcosθ;fx=ρgsinθ;fy=ρgcosθ;G2=ρgH,其中ρ為巖體密度,H為隧道埋深。另外,層狀巖層與隧道邊界相交的兩點間的圓弧對應的圓心角為δ,故L=Rsin(δ/2)。將以上表達式代入式(13)可得:

2 計算結果分析

爆破擾動載荷Fd=ρcv,其中ρ=2.3 g/cm3,縱波波速c=3.6 km/s,v為爆破振動速度。運用控制變量法,除以某參數(shù)作為變量的情況外,各參數(shù)取值為:δ=10°,H=20 m,R=3 m,h=0.4 m,θ=45°。將數(shù)據(jù)代入式(14),通過計算得出,不同爆破振動速度下巖梁正應力(σx)隨各影響因素的變化規(guī)律如圖4～圖9所示。需要說明的是,計算過程中直接取覆蓋層重力作為巖層的覆蓋壓力。由于實際施工中巖層存在應力擴散作用,本文計算所得結果偏大,如需精確計算,應結合地應力的實測數(shù)據(jù)。

圖4為不同爆破質點振速條件下巖層厚度與巖梁正應力(σx)的關系曲線。圖4顯示:σx隨著巖層厚度的增大而減小,兩者近似呈反比例關系。當隧道圍巖僅受靜力作用(v=0 cm/s)時,σx的變化規(guī)律與文獻[22]中武當山層狀隧道監(jiān)測得出的結論一致,表明本計算模型是合理的。隨著爆破質點振速的增加,巖體應力也不斷增大;但隨著巖層厚度的增加,σx的增幅逐漸減小。

圖5為巖層傾角與巖梁正應力的關系曲線。當巖層僅受靜力作用時,σx隨著巖層傾角的增加而減小,并且均小于1 MPa,因此既有層狀隧道在不受爆破擾動時是穩(wěn)定的。當有爆破擾動時,σx隨巖層傾角的增加而增加,但增加速率逐漸減小,在85°時達到最大值;σx的增長速率隨著爆破擾動的增加而逐漸增加;當巖層處于水平時,爆破動力擾動對巖梁正應力沒有貢獻,此時巖梁正應力主要是由靜壓力引起的。由于爆破動擾動易使巖層沿結構面因發(fā)生相互錯動而破壞,故新建隧道施工產(chǎn)生的爆破擾動易使既有層狀圍巖隧道的迎爆側巖體發(fā)生節(jié)理面間的剪切破壞。

圖4 巖層厚度與巖梁正應力關系曲線Fig.4 Stress vs.the thickness of rock layer

圖5 巖層傾角與巖梁正應力關系曲線Fig.5 Stress vs.the angle of rock layer

圖6 隧道半徑與巖梁正應力關系曲線Fig.6 Stress vs.tunnel radius

圖7 覆蓋層厚度與巖梁正應力關系曲線Fig.7 Stress vs.the thickness of covering layer

圖6為隧道半徑與巖梁正應力的關系曲線。當圍巖僅受靜力作用時,σx隨著隧道半徑的增加而增加,但是增長較慢,其值較小。當有爆破擾動時,σx與隧道半徑成正比,且增幅速率呈逐漸增加趨勢;σx的增長速率與爆破擾動強度成正比,這是因為當隧道半徑增加時,梁的長度增加,其剛度減小,故圍巖質點的動應力增大。圖7為覆蓋層厚度(H)與巖梁正應力關系曲線。由圖7可知:σx隨著覆蓋層厚度的增大而增大,且基本呈現(xiàn)線性關系,其斜率與爆破擾動無關,因此,當既有隧道僅受靜力作用且埋深達到一定深度時,要注意隧道的安全支護。

圖8為層狀圍巖所對應夾角與應力關系曲線。由圖8可知:σx及其增幅均隨著層狀圍巖所對應夾角的增大而增大,且爆破擾動下其增幅更加明顯。圖9為隧道半徑和層狀厚度比與巖梁正應力關系曲線。圖9顯示:σx隨著半徑與厚度之比的增加而增加,且呈冪指數(shù)規(guī)律變化;層狀圍巖所對應夾角越大,σx的增長速率越大;并且R/h與δ乘積相等時,σx取值相同。

圖8 層狀圍巖所對應夾角與應力關系曲線Fig.8 Stress vs.the corresponding angle of rock layer

圖9 隧道半徑和層狀厚度比與巖梁正應力關系曲線Fig.9 Stress vs.the ratio of tunnel radius to thickness

3 安全振速

根據(jù)最大主拉應力理論判據(jù),為了使隧道周壁不出現(xiàn)裂紋,則

式中:[σ]t為最大安全應力。由式(15)可得,爆破施工時,鄰近層狀隧道最危險點出現(xiàn)裂紋的臨界爆破振動速度(即最危險點所能承受的最大爆破質點振動速度)為

由式(16)可知,當埋深及層狀的傾角一定時,R/h與δ乘積越大,最大振動速度越小。根據(jù)既有隧道的半徑和埋深、層狀圍巖的厚度、層狀圍巖相對隧道的位置等可以計算得到既有隧道的安全質點振動峰值速度;再結合爆破質點振動衰減公式,即可確定出爆破施工時的最大單響藥量。

4 結 論

基于彈性力學和結構力學原理,建立了既有層狀圍巖隧道迎爆側最危險點的簡化力學模型,得到了最危險點的應力計算表達式,分析了各相關因素對既有隧道圍巖穩(wěn)定性的影響規(guī)律,得到以下結論。

(1)σx隨著爆破質點振速的增加而增大;但是隨著巖層厚度的增加,σx的增幅逐漸減小。

(2)σx隨著層狀傾角的增大而增大,且增大速率逐漸減小,在85°時達到最大值;σx增長速率隨著爆破擾動的增大而逐漸增大,當巖層處于水平時,σx主要是由靜壓力引起,爆破擾動沒有貢獻。

(3)σx及其增幅均隨著隧道半徑的增大而增大,并且σx的增幅與爆破擾動強度成正比。

(4)σx隨著覆蓋層厚度的增大而增大,兩者基本呈線性關系,其比例系數(shù)與爆破擾動無關。

(5)σx及其增幅均隨著層狀圍巖所對應夾角的增大而增大,且爆破擾動下其增幅更加明顯。

(6)σx隨著半徑與厚度之比的增大而增大,且兩者呈冪指數(shù)關系。

(7)利用計算所得最危險點最大振速的計算公式,結合爆破質點振動衰減公式,可以確定出爆破施工時的最大單響藥量。

參考文獻:

[1] 張志強,何本國,關寶樹.節(jié)理巖體隧道圍巖穩(wěn)定性判定指標合理性研究[J].現(xiàn)代隧道技術,2012,49(1):12-19. Zhang Zhiqiang,He Benguo,Guan Baoshu.A study of the criterion for the stability of a tunnel in a jointed rock mass[J].Modern Tunnelling Technology,2012,49(1):12-19.

[2] 胡文清,鄭穎人,鐘昌云.木寨嶺隧道軟弱圍巖段施工方法及數(shù)值分析[J].地下空間,2004,24(2):194-197. Hu Wenqing,Zheng Yingren,Zhong Changyun.Construction method and numerical analysis on Muzhailing tunnel [J].Underground Space,2004,24(2):194-197.

[3] 張黎明,鄭穎人,王在泉,等.有限元強度折減法在公路隧道中的應用探討[J].巖土力學,2007,28(1):97-101. Zhang Liming,Zheng Yingren,Wang Zaiquan,et al.Application of strength reduction finite element method to road tunnels[J].Rock and Soild Mechanics,2007,28(1):97-101.

[4] 鄭穎人,趙尚毅,鄧楚鍵,等.有限元極限分析法發(fā)展及其在巖土工程中的應用[J].中國工程科學,2006,8(12):39-61. Zheng Yingren,Zhao Shangyi,Deng Chujian,et al.Development of finite element limit analysis method and its application to geotechnical engineering[J].Engineering Science,2006,8(12):39-61.

[5] 鄭穎人,邱陳瑜,張紅,等.關于土體隧洞圍巖穩(wěn)定性分析方法的探索[J].巖石力學與工程學報,2008,27(10):1968-1980. Zheng Yingren,Qiu Chenyu,Zhang Hong,et al.Exploration of stability analysis methods for surrounding rocks of soil tunnel[J].Chinese Journal of Rock Mechanics and Engineering,2008,27(10):1968-1980.

[6] 鄭穎人,邱陳瑜,宋雅坤,等.土質隧洞圍巖穩(wěn)定性分析與設計計算方法探討[J].后勤工程學院學報,2009,25(3): 1-9. Zheng Yingren,Qiu Chenyu,Song Yakun,et al.Exploration of stability analysia and design calculation methods of surrounding rocks in soil tunnel[J].Journal of Logistical Engineering University,2009,25(10):1-9.

[7] 鄭穎人,徐浩,王成,等.隧洞破壞機理及深淺埋分界標準[J].浙江大學學報,2010,44(10):1851-1856. Zheng Yingren,Xu Hao,Wang Cheng,et al.Research on the failure mechanism of tunnel and dividing line standard between shallow and deep bury[J].Journal of Zhejiang University,2010,44(10):1851-1856.

[8] 邊克信.爆破地震對地下構筑物的影響[J].中國礦山工程,1981,1(6):25-36. Bian Kexin.Influence of blasting seismic on underground structure[J].China Mine Engineering,1981,1(6):25-36.

[9] 朱瑞賡,李錚.爆炸地震波的現(xiàn)場測量及其安全距離[J].爆炸與沖擊,1982,2(1):60-67. Zhu Ruigeng,Li Zheng.In-situmeasurements of blasting seismic waves and their safety distance[J].Explosion and Shock Waves,1982,2(1):60-67.

[10] Singh P K.Blast vibration damage to underground coal mines from adjacent open-pit blasting[J].International Journal of Rock Mechanics and Mining Sciences,2002,39(8):959-973.

[11] Torano J,Rodríguez R,Diego I,et al.FEM models including randomness and its application to the blasting vibrations prediction[J].Computers and Geotechnics,2006,33(1):15-28.

[12] 劉國華,王振宇.爆破荷載作用下隧道的動態(tài)響應與抗爆分析[J].浙江大學學報,2004,38(2):204-209. Liu Guohua,Wang Zhenyu.Dynamic response and blast-resistance analysis of a tunnel subjected to blast loading [J].Journal of Zhejiang University,2004,38(2):204-209.

[13] 陽生權.小線間距施工隧道爆破地震影響下既有隧道圍巖線性動力分析[J].工程爆破,1998,4(1):1-6. Yang Shengquan.Linear dynamic analysis on rock-media of old tunnel affected by blasting vibration of new tunnel with small line-clearance[J].Engineering Blasting,1998,4(1):1-6.

[14] 張慶松,李利平,李術才,等.小間距隧道爆破動力特性試驗研究[J].巖土力學,2008,29(10):2655-2660. Zhang Qingsong,Li Liping,Li Shucai,et al.Experimental study of blasting dynamic vibration of closely adjacent tunnel[J].Rock and Soild Mechanics,2008,29(10):2655-2660.

[15] 王明年,潘曉馬,張成滿,等.鄰近隧道爆破振動響應研究[J].巖土力學,2004,25(3):412-414. Wang Mingnian,Pan Xiaoma,Zhang Chengmam,et al.Study of blasting vibration influence on close-spaced tunnel[J].Rock and Soild Mechanics,2004,25(3):412-414.

[16] 畢繼紅,鐘建輝.鄰近隧道爆破震動對既有隧道影響的研究[J].工程爆破,2004,10(4):69-73. Bi Jihong,Zhong Jianhui.Study on influence of blasting vibration from excavation of a new tunnel on existed tunnel[J].Engineering Blasting,2004,10(4):69-73.

[17] 吳亮,鐘冬望.不同布置條件下鄰近隧道掘進爆破對既有隧道的影響[J].煤炭學報,2009,34(10):1339-1343. Wu Liang,Zhong Dongwang.Effect of tunneling blasting on the existing adjacent tunnel under different conditions [J].Journal of China Coal Society,2009,34(10):1339-1343.

[18] 王新宇,邵珠山,喬汝佳.小凈距下穿鐵路隧道爆破震動的響應研究[J].應用力學學報,2013,30(4):527-532. Wang Xinyu,Shao Zhushan,Qiao Rujia.Study on the dynamic responses of blasting vibration on undercrossing railway tunnel with small distance[J].Chinese Journal of Applied Mechanics,2013,30(4):527-532.

[19] 蘇永貴,趙躍堂.爆炸長波作用拱形坑道在不同側壓力系數(shù)下的應力與破壞分析[J].國防交通工程與技術,2015, 13(1):50-53. Su Yonggui,Zhao Yuetang.A stress-and-failure analysis of an arched tunnel under the action of the long-wave explosion with different lateral-pressure coefficients[J].Traffic Engineering and Technology for National Defence, 2015,13(1):50-53.

[20] 徐芝綸.彈性力學[M].北京:高等教育出版社,2006:41-45.

[21] 李廉錕.結構力學[M].北京:高等教育出版社,2004:130-132.

[22] 徐穎,陳建平,左昌群.武當群片巖隧道剪切型縱向裂縫破壞機理[J].中南大學學報(自然科學版),2013,44(7): 2950-2957. Xu Ying,Chen Jianpping,Zuo Cangqun.Shear type longitudinal cracks failure mechanism in Wudang group of schist tunnel[J].Journal of Central South University(Science Edition),2013,44(7):2950-2957.

Vibration velocity threshold of a tunnel adjacent to surrounding layered rocks under blasting load

Wu Liang1,2,Li Feng1,2,4,Lu Wenbo3,Chen Ming3,Xu Feng1,2
(1.College of Science,Wuhan University of Science and Technology,Wuhan430065,Hubei,China; 2.Blasting Technology Research Center,CRPCE-WUST,Wuhan430065,Hubei,China; 3.State Key Laboratory of Water Resource and Hydropower Engineering Science,Wuhan University,Wuhan430072,Hubei,China; 4.Guizhou Xinlian Blasting Engineering Group Co Ltd,Guiyang550002,Guizhou,China)

Under a blasting load,the front side adjacent to blasting in an existing tunnel surrounded by layered rocks is apt to be destroyed,and therefore it is of critical importance to keep the tunnel safe from the destructive effects of the surrounding rocks’cracks and damages that occur in a blasting disturbance.In this work,based on elastic mechanics and structural mechanics,we established a simplified mechanical model for the most dangerous point of the tunnel’s blasting side,and obtained the stress formula that can calculate the stress at this point.We also analyzed the regular influence of such different factors as the ratio of the tunnel radius to the thickness of layer,the inclination angle, the central angle and depth of the surrounding layered rocks on the stress.In addition,using the theory of the maximum tensile stress,we calculated the surrounding rocks’critical vibration velocity that ensures the stability of the tunnel and obtained,according to the attenuation patterns of seismic wave propagation,the maximum charge of detonation for the construction of a new tunnel.

stability of surrounding rock;blasting;tunnel;layered rock;vibration velocity

O382.2國標學科代碼:13015

:A

10.11883/1001-1455(2017)02-0208-07

(責任編輯 王玉鋒)

2015-09-07;

:2015-11-18

國家自然科學基金項目(51004079;51479147);湖北省自然科學基金項目(2014CFB822);水利部科技推廣計劃項目(TG1522)

吳 亮(1980- ),男,博士,副教授,wuliangwust@sina.com。