確定延時爆破間隔時間的新方法

張小軍 汪旭光 于亞倫 楊德強

(1.北京科技大學土木與資源工程學院，北京 100083；2.北京礦冶研究總院，北京 100160)

隨著城市的發展，現有的城市空間尤為不足，人們開始進行城市改造和向地下空間探索。但是城市拆除工程環境復雜、城市隧道埋深比較淺，并且毗鄰建筑物和各類地下管線，所以對爆破振動控制的要求非常苛刻。目前工程中常用的減振技術有分散裝藥、延時爆破、預裂爆破、減小斷面或進尺等，甚至采用靜態破碎的方法等，以保證爆破施工的順利進行[1-2]。

延時爆破是目前應用最廣泛的減振手段之一。

廖先葵等在礦山大規模爆破中，采用對稱分區，延時爆破等方法，實現了爆破地震波干擾減振的目的[3]；田振農等通過時頻分析方法對隧道爆破振動的一般特征進行了分析，結果表明爆破地震波是多列爆炸波疊加的結果[11]；傅洪賢等結合蘭渝鐵路的建設，在人和場隧道進行電子雷管降振試驗、電子雷管與非電雷管聯合降振試驗，得出隧道電子雷管單孔連續起爆降振技術，與隧道非電雷管爆破相比，爆破振動降低80%以上[8]；根據Basebi等的工程試驗研究[10]，采用單孔單響的電子雷管起爆，爆破振動可以降低約三分之一；高文學等結合基礎開挖工程，探討了電子雷管起爆技術及其在復雜環境下深孔爆破工程的應用，研究表明電子雷管可以精確且合理地延時時間，可以達到良好的降振效果，結合預裂控制爆破，降振效果更加顯著[12]。上述研究基本都證明了延時爆破可以達到減振的目的，并且效果顯著，但是都沒有明確提出確定延時爆破間隔時間的方法。

本研究主要依據波的疊加原理和波的干涉，對爆破地震波的疊加效應進行了理論分析，得出延時爆破的最佳間隔時間。同時以理論分析為基礎，利用波形峰值包絡線的思想與爆破地震波疊加的波動方程，推導出了爆破地震波疊加的簡析式。最后通過案例，利用單炮次爆破振動波形圖，摘取其峰值數據，擬合出了測點處質點振動速度的表達式，根據不同延時間隔時間爆破振動的理論波形圖，確定出最佳的延時爆破間隔時間。

1 爆破地震波的疊加規律

爆破地震波的傳播是極為復雜的，質點振動強度的衰減規律受到諸多因素的影響。在延時爆破中，某質點的振動受到若干炮孔炸藥爆炸的影響。為了實現減振的目的，需要對雷管的延時間隔時間進行確定。以下是通過建立模型，理論計算雷管的精確延時間隔時間。

設相鄰的2個炮孔A和炮孔B距測點C的距離近似相等，爆破引起的振動近似為簡諧振動。炮孔和測點示意圖如圖1所示。

圖1 炮孔與測點

圖1中，A、B為相鄰的2個炮孔，C為振動測點；A、B到C點的距離近似相等，為X。

由炮孔A爆炸產生的振動波形圖如圖2。其波動方程為

y=asinwt,

(1)

式中，a為振幅；W為圓頻率；t為時間；y為A點處質點在時刻t離開其平衡位置的位移。

圖2 炮孔A爆炸產生的振動波形

由于振動是從A點到C點傳播的，所以C點處質點的振動落后于A點，其落后時間為X/c，c指波的傳播速度，所以C點在時刻t的位移等于A點處質點在時刻(t-X/c)的位移。所以C點處質點振動的波動方程為

(2)

因為

(3)

cy=λ,

(4)

式中，f為頻率；λ為波長；T為周期。

即

(5)

為了實現炮孔A和炮孔B是延時爆破，設炮孔B延時Δt爆炸，其波動方程為

y=asinw(t-Δt).

(6)

由于振動是從炮孔B到測點C傳播的，所以C點處質點的振動落后于B點,其落后時間為X/c，所以C點在時刻t的位移等于B點質點在時刻(t-X/c)的位移，所以C點處質點振動的波動方程為

(7)

根據式(3)、式(4)，式(7)變換為

(8)

爆破地震波是由不同頻率、不同幅值的波在一個有限時間范圍內組合的隨機過程。地震波的最大幅值、頻率和持續時間被稱為地震波的3要素。而最大振幅又與速度、加速度密切相關，故地震波的速度、頻率和持續時間是表征地震波強度的3個必不可少的參量，其中振速是最常用的，那么式(5)對時間t進行求導,即為A炮孔爆炸C處質點振動速度：

(9)

式(8)對時間t進行求導,即為B炮孔爆炸C處質點振動速度：

(10)

炮孔A和炮孔B的振動傳播到測點C時，C點的振動就是2列振動波的疊加，那么測點C處質點振動的合速度為

vc=vac+vbc=

(11)

根據振動的合成公式，式(11)合成為

(12)

其中，

(13)

(14)

從式(11)、式(13)中可以發現，這2列波在空間任一點的位相差為

Δφ=2πcΔt/λ,

(15)

而Δφ的取值影響測點C振動合成式(13)中振速峰值b的大小。所以當：

(16)

根據式(13)，測點C處質點的振動速度峰值最大。

當：

(17)

根據式(13)，測點C處質點的振動速度峰值最小，那么延時爆破的最佳間隔時間為

k=0,1,2,…

(18)

式中，T為質點振動周期。

所以延時爆破的最佳間隔時間，即半個周期或半個周期的奇數倍。

2 振速包絡線以及合成振速近似簡析式

事實上，爆破振動并不是簡單的簡諧振動，質點在振動過程中總會受到一些阻力，振動的能量要不斷轉化為其他形式的能量。根據振動物體的動能公式式(19)，振動物體的動能跟振速的平方成正比，動能隨時間減少，振速也隨時間減小，如果不給予能量補充，那么經過一段時間后，振動就要停止，如圖3所示。

(19)

式中，Ek為振動物體動能；m為質量；v為物體振動速度。

為了更好地描述爆破振動的特性，引入爆破振動速度波形峰值包絡線的概念，通過爆破振動速度波形峰值包絡線的簡析式與第1節中依據爆破地震波疊加理論推導的公式相結合，得出測點處質點振動速度的近似解析式。

2.1 振速包絡線

不管在爆破振動的監測過程中還是對爆破振動效應的預防中，往往比較關心的是爆破振動的峰值速度。那么通過對爆破振動速度波形峰值摘取，然后對這些摘取的數據進行擬合，即可得到光滑的爆炸地震波波形峰值包絡線，如圖3所示。

圖3 爆破振動速度

從圖3中以及其他的爆破振動速度波形圖中，可以得到爆破振動波形峰值包絡線的近似解析式，其形式如下：

V=k′e-βt,

(20)

式中，V為介質質點振動速度峰值，cm/s；k′為與炮孔炸藥量、爆心距等有關的系數；β為與爆破場地環境有關的指數；t為時間，s。

其中k′跟炮孔的裝藥量、測點距爆源的距離、爆破場地環境等有關系，可以近似看作：

(21)

則式(20)變換為

(22)

2.2 合成振速近似解析式

為了研究爆破振動的質點振動速度，需要知道質點振動速度的簡析式，根據第1節的理論推導，將式(9)作變量代換，令：

(23)

則式(9)變換為

(24)

同理將式(10)作變量代換，令：

(25)

則式(10)變換為

(26)

所以式(25)、式(26)合成為

(27)

即：

(28)

式中，

(29)

根據式(22)則：

(30)

(31)

(32)

根據實際情況，測點處的質點振動要歷經3個階段，第1個階段是傳播階段，即從炮孔爆炸到爆破振動傳遞到測點處，需要一定的時間，在這個階段測點處的質點振動速度為零；第2個階段是單振動階段，即最先傳遞到測點處的爆破振動引起測點處質點振動；第3個階段是雙振動階段，即2炮孔爆炸引起的振動都傳遞到測點處引起測點處質點振動,并且各個階段之間的振動波形曲線會不連續，會發生跳躍。所以測點處質點振動速度的表達式為

(33)

3 案 例

3.1 工程概況

某居住小區Ⅱ標土石方爆破工程東面緊鄰公路，西面160 m為一些土窯洞，這些土窯洞將要作為歷史文物，需要保護。北面為已開挖完成的施工場地，南面為山體。周邊環境如圖4所示。為了防止爆破對土窯洞造成影響，爆破工程采取延時爆破減振措施。

圖4 周邊環境

3.2 數據擬合

為了確定延時爆破間隔時間，首先進行單炮次試爆，根據單炮次的爆破振動速度波形圖，進行數據擬合。單炮次藥量15 kg，試爆地點選在施工區域爆破點附近，監測點距離160 m。其測點處質點振動速度波形見圖5。

根據單炮次爆破振動速度波形圖，選取圖中峰值的8個點，見圖5中的標記，將其整理到表1。

根據表1中的8組數據，按照式(22)進行非線性曲線擬合，得到單炮次爆破振動速度峰值包絡線，見圖6。

圖5 單炮次爆破振動波形

項 目時間/s振速/(cm/s)數據10．0150．140數據20．0620．070數據30．1200．060數據40．2000．020數據50．2300．021數據60．2750．024數據70．3000．010數據80．3320．010

圖6 單炮次爆破振動速度峰值包絡線

其擬合公式為

(34)

3.3 不同延時間隔時間測點處質點的振動速度

根據段如泰在《用地震定位方法研究礦山爆破中的地震波傳播速度》一文中的敘述，爆破地震波的傳播速度為v=4 498.85 m/s(或王進強在《地震波傳播速度原位試驗及計算》一文中敘述爆破地震波的平均傳播速度為4 810 m/s)，其地質條件與該爆破工程類似，所以在該工程中爆破地震波的傳播速度取 4 498.85 m/s[14]；爆破地震頻率在10～100 Hz，這里取頻率25 Hz，周期為0.04 s，波長為179.95 m。不同炮孔距測點的距離看作近似相等?，F對不同延時間隔時間的2炮孔爆炸引起測點處質點振動速度進行討論，每個炮孔取藥量為15 kg。

(1)Δt=0 s。根據式(29)、式(32)得：

A=2×0.135 8e-7.932(t-0.035 5),

(35)

tanφ=0.842 6,φ=0.699,

(36)

則測點處質點振動速度方程為

(37)

測點處質點振動速度波形圖見圖7所示。

圖7 Δt=0 s測點處質點振動速度波形

(2)Δt=0.02 s。根據式(29)、式(32)得：

A=0.030 9e-7.932t

(38)

tanφ=0.859 0,φ=0.709 7

(39)

則測點處質點振動速度方程為

(40)

測點處質點振動速度波形圖見圖8所示。

圖8 Δt=0.02 s質點振動速度波形

(3)Δt=0.025 s。根據式(29)、式(32)得

A=0.156 8e-7.932t，

(41)

tanφ=-1.658 5,φ=-1.028 2，

(42)

則測點處質點振動速度方程為

(43)

測點處質點振動速度波形圖見圖9所示。

圖9 Δt=0.025 s質點振動速度波形

(4)Δt=0.05 s。根據式(29)、式(32)得：

A=0.323 0e-7.932t，

(44)

tanφ=-0.283 849,φ=-0.276 6，

(45)

則測點處質點振動速度方程為

(46)

測點處質點振動速度波形圖見圖10所示。

圖10 Δt=0.05 s 質點振動速度波形

從圖7、圖8、圖9、圖10中可以看出炮孔爆炸的瞬間，測點處的質點振動速度為零，也就是地震波還沒有傳播到測點處，其傳播需要一定的時間；圖7中測點處質點振動速度峰值最大。圖8中測點處質點振動速度峰值最小，但是剛開始的時候還是比較大的，不可忽視。

根據圖8，測點處質點振動速度峰值最小，那么該爆破工程選擇延時間隔時間為0.02 s的雷管是最佳的，但是現場選用的是奧瑞凱高精度導爆管雷管，其段別有限，只有延時25 ms、50 ms等的雷管供選用，那么通過比較圖9和圖10，間隔0.025 s的振速峰值要比間隔0.05 s的振速峰值小，同時圖10中可以看出會出現反復的情況，即振速將低后有驟升的現象，所以該爆破工程建議選用延時間隔時間為25 ms的高精度導爆管雷管來實施爆破。

3.4 現場爆破

根據上述的理論分析，實際單孔裝藥量15 kg，最大段藥量15 kg，總裝藥量420 kg，孔內500 ms高精度導爆管雷管，孔間25 ms的高精度導爆管雷管實施爆破，在測點處監測到的波形圖如圖11所示。

圖11 測點處爆破振動速度波形

從圖11可以看出選用延時間隔25 ms的高精度雷管，測點處質點的振速大部分都在0.1 cm/s之內，只有個別超出0.1 cm/s，但小于0.15 cm/s，與圖9的波形圖相吻合，基本沒有對窯洞以及其他設施造成破壞，爆破取得了成功。

4 結 論

(1)通過建立數學模型、理論分析，得出了延時爆破的最佳間隔時間，即半個周期或半個周期的奇數倍。

(2)利用單炮次的爆破試驗，推導出了2孔延時爆破振動速度的理論簡析式，從而可以畫出測點處質點振動的速度波形圖，使研究問題更加形象化。

(3)2孔延時爆破測點處的質點振動要歷經3個階段，分別是傳播階段、單振動階段、雙振動階段。

(4)根據測點處質點振動速度表達式畫出的理論波形圖，可以確定最佳的延時爆破間隔時間。如果現有爆破器材實現不了最佳間隔時間，可以根據理論波形圖，選擇較優的延時爆破間隔時間。

(5)本研究對延時爆破間隔時間的確定提供了新的思路與方法。

[1] 孫冰峰,傅洪賢,管曉明，等.基于建筑物模態分析的隧道爆破減振技術[J].中國公路學報,2016,29(10):86-94.

Sun Bingfeng,Fu Hongxian,Guan Xiaoming,et al.Decreasing vibration technology for tunnel blasting based on building modal analysis[J].China Journal of Highway and Transport,2016,29(10):86-94.

[2] 李鎮明,張 瑀,周廣春.單自由度體系有阻尼自由振動振幅包絡線的性質研究與數學證明[J].數學的實踐與認識，2015,45(21):227-232.

Li Zhenming,Zhang Yu,Zhou Guangchun.Characteristics of the envelope of damped free vibrations of SDOF systems and its proof[J].Mathematics in Practice and Theory,2015,45(21):227-232.

[3] 廖先葵.深孔爆破大規模干擾減震機理及其應用[J].爆 破,1993(4):1-6.

Liao Xiankui.Large scale interference shock absorption mechanism and its application in deep hole blasting[J].Blasting,1993(4):1-6.

[4] 劉清泉,李玉民.城市地下工程的減震爆破技術[J].煤炭科學技術,1991(10):13-15.

Liu Qingquan,Li Yumin.Cushion blasting technique for city underground project[J].Coal Science and Technology,1991(10):13-15.

[5] 薛孔寬,唐光榮,曹恒安，等.分段微差爆破地震效應的疊加分析[J].爆破,1991(3):67-71.

Xue Likuan,Tang Guangrong,Cao Heng′an,et al.Overlapping analysis of seismic effects of particle borehole blasting[J].Blasting,1991(3):67-71.

[6] 魏曉林,劉 翼.爆破拆除樓房時塌落振動的預測[J].工程爆破,2016,22(2):13-18.

Wei Xiaolin,Liu Yi.Prediction of collapse vibration of building demolished by blasting[J].Engineering Blasting,2016,22(2):13-18.

[7] 張 迪.地鐵隧道精確延時爆破振動傳播規律與控制試驗研究[D].北京:中國礦業大學(北京),2016.

Zhang Di.Vibration Propagation Regularity and Control Experiment Research of Metro Tunnel with Precise Delay[D].Beijing:China University of Mining and Technology (Beijing),2016.

[8] 陳繼強,劉為洲.多孔爆破振動強度的單孔波形疊加計算[J].金屬礦山,2000(8):23-25.

Chen Jiqiang,Liu Weizhou.Single hole waveform superimposed calculation of the vibration strength of multihole blasting[J].Metal Mine,2000(8):23-25.

[9] 陳士海,燕永峰,戚桂峰，等.微差爆破將震效果影響因素分析[J].巖土力學,2011,32(10):303-308.

Chen Shihai,Yan Yongfeng,Qi Guifeng,et al.Analysis of influence factors of interence vibration reduction of millisecond blasting[J].Rock and Soil Mechanics,2011,32(10):303-308.

[10] Basebi T,Thomas R M.A study of moving mesh methods applied to a thin flame propagating in a detonator delay element[J].Computers and Mathematics with Applications,2003,45(1):131-163.

[11] 田振農,孟祥棟,王國欣.城區隧道電子雷管起爆錯相減震機理分析[J].振動與沖擊,2012,31(21):108-111.

Tian Zhennong,Meng Xiangdong,Wang Guoxin.Mechanism analysis of fault-phase vibration reduction for tunnel blasting initiated by electronic detonators in city area[J].Journal of Vibration and Shock,2012,31(21):108-111.

[12] 胡建華,尚俊龍,羅先偉,等.單孔爆破振動監測與衰減規律多元線性化回歸[J].振動與沖擊,2013,32(16):49-53.

Hu Jianhua,Shang Junlong,Luo Xianwei,et al.Monitoring of single-hole blasting vibration and detection of its attenuation law by using multiple linear regression[J].Journal of Vibration and Shock,2013,32(16):49-53.

[13] 喻 軍,劉松玉,童立元.淺埋隧道爆破振動空洞效應[J].東南大學學報:自然科學版,2010,40(1):176-179.

Yu Jun,Liu Songyu,Tong Liyuan.Hollow effect induced by blasting vibration in shallow tunnels[J].Journal of Southeast University:Natural Science Edition,2010,40(1):176-179.

[14] 段如泰,金振奎,索重輝,等.用地震定位方法研究礦山爆破中的地震波傳播速度[J].科學技術與工程，2011,11(36):8967-8977.

Duan Rutai,Jin Zhenkui,Suo Chonghui,et al.Earthquake location method used in velocity study of mine explosion stimulating seismic wave[J].Science Technology and Engineering,2011,11(36):8967-8977.

[15] 胡鶴飛.MATLAB及應用[M].北京:北京郵電大學出版社,2012.

Hu Hefei.MATLAB and Its Application[M]:Beijing University of Posts and Telecommunications Press,2012.

[16] 李洪濤,盧文波,舒大強,等.爆破地震波的能量衰減規律研究[J].巖石力學與工程學報,2010,29(1):3364- 3369.

Li Hongtao,Lu Wenbo,Shu Daqiang,et al.Study of energy attenuation law of blast-induced seismic wave[J].Chinese Journal of Rock Mechanics and Engineering,2010,29(1):3364-3369.

[17] 龔 敏,趙振振,吳昊駿,等.針對建構筑物不同振速要求的隧道爆破藥量計算方法[J].煤炭學報,2016,41(7):1747-1753.

Gong Min,Zhao Zhenzhen,Wu Haojun,et al.Computational method of blasting charges under different vibration velocity limits of buildings and structures around tunnel[J].Journal of China Coal Society,2016,41(7):1747-1753.