減振孔減振效應的數值模擬研究

趙蕊　李雅詩　王建新　彭琦　唐超華　朱建波

摘 要：利用顯式動力有限元軟件AUTODYN研究了減振孔對巖石地下洞室爆破產生的爆炸應力波傳播與衰減以及對地表和鄰近地下洞室動態響應的影響與規律。通過與爆破模型試驗結果進行對比，驗證了AUTODYN模擬地下洞室爆炸過程和應力波傳播的有效性和準確性。通過數值模擬分析了不同減振孔參數對鄰近洞室和地表的質點峰值振動速度（PPV）和減振孔屏障前后測點的隔振率的影響，揭示了減振孔參數變化對減振效果的影響規律。數值模擬結果表明：減振孔直徑、孔間距、孔距爆源距離及孔排數均會對鄰近地下洞室和地表的PPV及減振孔屏障前后測點的隔振率產生影響。隨著減振孔直徑、減振孔排數的增大，或者減振孔間距、減振孔距爆源距離的減小，鄰近地下洞室和地表測點的PPV呈現減小趨勢，減振孔屏障前后測點的隔振率呈現顯著增大趨勢，即減振孔減振效果增強。此外，隨減振孔參數變化，減振孔前后和減振孔間巖體前后測點的隔振率呈現顯著的差異性，前者明顯大于后者，兩者差值隨孔直徑的增加先增大后減小，隨爆心距的增加迅速減小，隨孔排數的增加沒有明顯變化。減振孔間距變化會影響測點與減振孔的相對位置，從而影響測點的隔振率大小。

關鍵詞：爆破;減振孔;質點峰值振動速度;隔振率

中圖分類號：TU94? ?文獻標志碼：A ??文章編號：2096-6717（2021）02-0049-11

Abstract： The explicit dynamic FEM-based code AUTODYN was used to study the effect and law of barrier holes on blast-induced stress wave propagation and attenuation， as well as the dynamic responses of adjacent underground opening and ground. First of all，? compared with the results of blasting model test， the effectiveness and accuracy of AUTODYN? to simulate the process of underground opening explosion and stress wave propagation are verified. Then， the effects of different barrier hole parameters on peak particle velocity （PPV） of the adjacent chamber and ground， and the vibration-isolation rates of measuring points before and after barrier hole screen were analyzed through the numerical simulation method， and the pattem of varying barrier hole parameters on vibration reduction effect is revealed. The simulated results indicated that barrier hole diameter， hole spacing， the distance between hole and explosion center， and hole row have effects on PPVs of adjacent underground opening and ground and the vibration-isolation rates of measuring points before and after barrier hole screen. With increasing barrier hole diameter and barrier hole row， or decreasing barrier hole spacing and distance between barrier hole and explosion center， PPVs of the adjacent underground opening and ground show a decreasing trend and vibration-isolation rates of measuring points before and after barrier hole screen show a great increasing trend， that is， the vibration reduction effect of barrier hole? is enhanced. In addition， with change of barrier hole parameters， the vibration-isolation rates of the measuring points before and after barrier hole and the measuring points before and after rock mass between barrier holes show significant differences. The former is obviously larger than the latter， and the difference between the two increases first and then decreases with increase of hole diameter， and decreases rapidly with increase of the distance between hole and explosion center. But it doesnt change significantly with the increase of barrier hole row. The change of barrier hole spacing could affect the relative position of the measuring point and the barrier hole， thus determining the vibration-isolation rate of the measuring point.

Keywords：blasting; barrier hole; peak particle velocity （PPV）; vibration-isolation rate

工程爆破作為主要的施工手段被廣泛應用于礦山開采、隧道開挖、基礎設施建設等工程領域，與此同時，也帶來了諸多負面效應。炸藥在巖石中爆炸產生的能量主要用于對周圍介質做功，另有部分能量以應力波的形式向外傳播[1-2]。爆炸應力波的傳播對爆破施工場地附近的建（構）筑物、邊坡、鄰近地下洞室等會產生嚴重危害[3-6]。因此，如何在不影響工程質量、工期的前提下利用減振技術降低爆破振動對周邊的影響，一直是工程爆破亟待解決的關鍵問題。

天然巖體中存在如斷層、節理、層理等不連續面。爆炸應力波在巖體中傳播時，在不連續面處會發生反射、折射，導致波的能量減少、幅值減小、振動速度變小[7-9]。通過改變傳播介質的物理性質，影響爆炸應力波在巖體中的傳播與衰減，是目前地下工程爆破施工中行之有效的減振手段之一。依據此減振機理，在震源和被保護體之間設置減振孔，可以加速爆炸應力波的衰減，減小爆破振動對周圍構筑物造成的影響。減振孔鉆孔既不會產生額外振動，也不受開挖巖體條件的限制，具有施工工藝簡單、布置方式靈活、對環境影響破壞小等優點，現已廣泛應用于基坑、礦山、邊坡、隧道爆破開挖、人防工程等領域[10-12]。

許多學者通過試驗分析和數值模擬對減振孔的減振效果及影響因素進行研究。試驗分析方面，徐亮[13]通過現場試驗研究減振孔布置對地表振動的影響。結果發現，布置減振孔可明顯減小水平切向爆破能量，而對水平徑向和垂直向的能量分布無明顯影響。Uysal等[14]進行了3排減振孔的現場爆破試驗，與孔前測點PPV相比，孔后測點PPV降低了18%。鄭文富[15]通過模型試驗發現，增加減振孔排數可導致徑向主頻增大，而對切向和垂直向主頻影響不大。

數值模擬研究方面，Lei等[16]通過LS-DYNA數值模擬發現，增大減振孔直徑和減小孔間距，可以明顯提高減振率，擴大振動屏蔽范圍。楊典光[17]通過數值模擬研究發現，增加減振孔直徑和孔排數可以有效降低振動速度，但并未進一步探究機理。Bian等[18]采用離散元UDEC模擬研究減振孔對爆破地表振動的影響，通過對比水平振速和垂直振速發現，減振孔參數變化對垂直方向振速影響更大。Dohyun[19]通過AUTODYN進行隧道爆破施工時減振孔對地表振動的減振研究，對減振率隨孔參數變化進行曲線擬合，提出安全系數1.2。惠峰等[20]通過數值模擬研究了小凈距隧道爆破中減振孔排數和爆心距對爆破應力波衰減的影響，但參數設置較少。

目前，減振孔減振效果的研究大多局限于爆破施工過程中減振孔對地表振動的影響，而對于巖石地下洞室爆破中減振孔對鄰近洞室的動態響應以及對減振孔屏障前后測點的PPV影響規律的認識不足。鑒于此，通過對減振孔的地下洞室爆破進行數值模擬，系統研究了不同孔參數對地下鄰近洞室和地表振動響應的影響規律，探究了減振孔減振機理及減振孔減振效果機制，深入分析了減振孔屏障前后不同位置測點隔振率隨孔參數變化的差異性。

1 數值模擬方法及其試驗驗證

AUTODYN是一種顯式動力有限元分析軟件，用來解決固體、流體、氣體及其相互作用的高度非線性動力學問題[21]。由于其具有界面友好、集成方便、材料庫完善、求解復雜問題快速、求解精度高等特點，現已廣泛應用于軍工、航空、巖土、化學、汽車等領域。

1.1 試塊制作

模型試驗中長方體水泥砂漿試塊模型如圖1（a）、（b）所示，試塊尺寸為800 mm×600 mm×400 mm，模型相似比例為1∶20。兩個洞室直徑為200 mm，凈間距為200 mm，洞室頂部距離試塊上表面為200 mm，長400 mm，貫通試塊。左邊為爆炸洞室，右邊為鄰近洞室。在兩個洞室之間設置1排減振孔，減振孔間距為50 mm，減振孔距爆源距離為150 mm，減振孔直徑分別為5、10、15 mm，每種直徑減振孔模型設置兩次平行試驗。

砂漿試塊采用普通硅酸鹽水泥（PC32.5）和細黃砂配制而成。在制備砂漿試塊的同時，用相同的混合物澆筑5塊100 mm×100 mm×100 mm的標準試樣。用CTS-25非金屬超聲波探測儀測得標準試樣的縱波速度。用TJW-1000型電液伺服巖石多功能試驗機對其進行單軸壓縮試驗，測得其抗壓強度和彈性模量。在計算出試塊的密度、縱波速度、彈性模量之后，由式（1）求出試塊的泊松比。取5個試塊所測數據的平均值，最終得到砂漿試塊的物理力學參數，如表1所示。

試驗所用炸藥組成為含5%石墨的2 g黑索金 （RDX）粉末和含1 g RDX的雷管，見圖1（c）。玻璃管內徑為15 mm，RDX粉末裝入玻璃管后高13.5 mm，RDX裝藥密度為1 260 kg/m3。用紙板將RDX藥柱卡于爆破洞室中央，洞室兩端用黃泥填堵密實，使剩余洞室空腔長10 cm，則RDX裝載密度為0.95 kg/m3。裝載密度為炸藥總質量與爆破洞室體積之比[2]。

1.2 試驗測試系統

由于爆破測振傳感器難以預埋于試塊內部，許多學者[22-23]在進行爆破模型試驗時通常采用預埋應變片的方式來監測試塊的爆破振動情況。基于此，對水泥砂漿試塊進行超動態應變測試。測試系統主要由應變片、橋盒、LK2107A型超動態應變儀、TST3406動態測試分析儀（含采集卡）、特制靜電磁屏蔽線組成，用于監測試塊中由爆炸應力波引起的應變變化情況，該測試系統示意圖如圖1（d）、（e）所示。

1）試驗中采用型號為BE120-3AA（11）、120 Ω的箔式環氧酚醛類電阻應變片。其敏感柵尺寸為2.8 mm×2 mm，基底尺寸為6.4 mm×3.5 mm，靈敏系數為2.11%±1%，引線為4 cm鍍銀銅線。每塊模型預埋3個應變磚，如圖1（b）所示，分別位于鄰近洞室左側壁1、減振孔前方2和減振孔后方3。應變磚尺寸為20 mm×20 mm×10 mm，應變磚埋深為300 mm，每個應變磚上粘貼1個徑向應變片，如圖1（f）所示，用來記錄測點沿x方向的應變（圖1（b））。應變磚與水泥砂漿試塊的材料配比相同，以保持兩者的彈性模量、波阻抗等一致，在應變片表面涂抹環氧樹脂進行防潮絕緣處理后，澆筑于水泥砂漿試塊中。

2）LK2107A型超動態應變儀能夠測量材料及結構發生任意變形時產生的動態應變，可進行動力強度研究及試驗應力分析。

在測量中，為防止周圍電磁干擾，影響監測結果，將測量導線與電源線絞扭，測量導線方向垂直于電源線方向，使得每一絞的感應電流與下一絞的感應電流相反，在絞扭線的外面采用較厚的銅金屬屏蔽套包裹，并將屏蔽套兩端接地。對信號電路采用一點接地方式，即儀器的外殼和銅金屬屏蔽套在同一點接地[23]。

1.3 AUTODYN數值模型

將爆破模型試驗簡化為一個2D平面應變問題，基于AUTODYN程序，建立減振孔直徑變化時RDX炸藥爆炸的2D仿真模型，模型尺寸及測點布置和試驗相同，如圖1（b）所示。模型長800 mm，高600 mm。在模型水平中軸線處布置3個測點，測點1位于鄰近洞室左側壁，測點2位于減振孔屏障前方，測點3位于減振孔屏障后方。模型邊界設為自由邊界。建立的AUTODYN模型材料包括巖石、炸藥（RDX）和空氣。巖石簡化為均勻、各向同性的彈性材料，材料參數選取爆破模型試驗中水泥砂漿試塊的物理力學參數，如表1所示。空氣材料從軟件自帶的材料庫中選擇。RDX能量方程滿足JWL（Jones–Wilkins–Lee）狀態方程[21]

空氣和RDX炸藥選用多材料Euler算法模擬，巖石材料用Lagrange算法模擬。巖石材料和空氣材料通過Lagrange-Euler耦合實現接觸，Lagrange-Euler耦合是解決流固氣耦合建模問題常用的有效方法。模型網格采用漸變網格，在炸藥和減振孔處網格劃分較密，遠離此處，網格逐漸增大。最小網格為0.3 mm，最大網格為3 mm，大約25萬個網格單元，提高計算速度的同時也保證了計算精度。模型網格如圖2所示。在AUTODYN-2D模型中，認為炸藥沿爆破洞室空腔長度方向均勻分布，則裝載密度可以簡化為式（3）[2]。

1.4 對比驗證

通過對比分析地下洞室爆炸中不同減振孔直徑下巖石測點的隔振率和鄰近洞室的峰值應變模型試驗結果和數值模擬結果（表3），驗證AUTODYN軟件模擬本研究的有效性和準確性。圖1（b）中，減振孔屏障前后的兩個對稱測點2和3組成一組，測得其峰值應變，通過式（4）求得該組測點的隔振率，作為表征減振孔減振效果的一種方法。

當孔徑為5 mm時，試驗和模擬所得測點1、2、3的應變時程曲線如圖3所示。測點的峰值應變取曲線的首波峰值（第一峰值應變）。圖中顯示爆炸產生的應力波在砂漿介質中傳播時，在測點位置首先產生壓縮應變;隨著應力波的傳播，遇到孔隙或塊體邊界等自由面時，生成反射拉伸波，產生拉伸應變;隨后壓應變和拉應變交替出現，形成小的波形波動。模擬和試驗得到的首波峰值接近，但試驗測得的首個反射拉伸波峰值大于數值模擬的結果。這可能與試驗中水泥砂漿試塊的非均質性和多相多孔性有關，其存在導致應力波多次反射后疊加產生較高幅值的反射拉伸波，從而生成較大的拉伸應變。

圖4和圖5分別為隨減振孔直徑變化時，鄰近洞室左壁測點1水平方向第一峰值應變和隔振率的試驗與模擬結果對比圖。可以看出，隨著減振孔直徑的增加，鄰近洞室左壁測點1的第一峰值應變逐漸減小，減振孔屏障前后測點的隔振率逐漸增大。這表明減振孔對爆炸應力波有屏蔽作用，且當減振孔直徑增大時，爆炸應力波通過減振孔時能量消耗增大，隔振效果增強。

對比數值結果和試驗結果可以看出，兩者變化趨勢一致，大小相近。但鄰近洞室左壁測點1的第一峰值應變的數值模擬結果大于試驗結果，隔振率的數值模擬結果小于試驗結果。主要原因可能是：1）試樣存在不均勻分布的問題，且由于攪拌不均而含有大量孔隙，從而導致爆炸應力波衰減;2）在有限元軟件模擬時，忽略了巖石材料（水泥砂漿介質）本身阻尼對應力波的衰減作用。因此，數值計算得出的鄰近洞室的應變峰值大于試驗結果，相應的隔振率小于試驗結果。

表4為模擬與試驗結果對比的誤差分析。由表4可以看出，相對誤差大多在20%以內，屬于工程允許范圍。當孔徑為10 mm時，隔振率的相對誤差超過了20%，這可能是由于該次模型試驗測量儀器誤差或模型試塊非均勻性較高，從而導致爆炸應力波通過減振孔時衰減較大。

綜上所述，有限元軟件AUTODYN可有效、準確地模擬地下洞室爆炸、爆炸應力波傳播以及減振孔減振效果等問題。

2 減振孔減振效應的數值模擬

利用AUTODYN建立含減振孔的2D爆破隔振效應模型，分別探討了減振孔直徑、孔間距、孔距爆源距離和孔排數對爆炸應力波的傳播和衰減的影響。

2.1 數值模型

采用的數值模型如圖6所示，模型尺寸為30 m×18 m，兩個圓形地下洞室直徑為6 m，洞室頂部距離模型上表面為6 m，兩個洞室凈間距為6 m，左側洞室為爆炸洞室，右側為鄰近洞室。在兩個洞室之間布設減振孔。共布置16個測點，其中測點1和測點2位于右側鄰近洞室左側壁和頂部，用于監測鄰近洞室的動態響應過程。測點3至測點6位于鄰近洞室上方的地表，間距為2 m，測點3位于鄰近洞室左側壁正上方，用于監測地表的動態響應過程。在減振孔屏障前后各布置了5個測點，測點7-8、11-12、15-16位于減振孔前后，測點9-10、13-14位于減振孔間巖體前后，用于監測爆炸應力波通過減振孔前后的能量變化。

為了避免模型四周反射的應力波影響數值結果的準確性，對模型左右及下邊界進行人工無反射處理，模型上邊界是地表，為自由邊界。巖石參數見表5[8]，炸藥采用模型試驗中的RDX參數，見表2。炸藥裝載密度取15 kg/m3，由式（3）算得，RDX炸藥半徑為327 mm。模型網格采用漸變網格，由炸藥和減振孔向外，網格由密變疏，最小網格為1 cm，最大網格為17.5 cm，大約18萬個網格單元。模型網格如圖7所示。

2.2 減振孔直徑對減振效果的影響

數值模型如圖6（a）所示，S1為0.5 m，S2為4 m，n為1排。

圖8（a）、（b）、（c）分別為對應不同減振孔直徑下地下洞室測點1-2處PPV、地表測點3-6處PPV、減振孔屏障前后測點7-16處隔振率變化曲線。從數值模擬結果可以看出，測點1-2（鄰近洞室左側壁和洞頂）和測點3-6（地表）處PPV均隨著減振孔直徑的增加呈明顯的減小趨勢，說明大直徑減振孔可更有效衰減爆炸應力波，與試驗結論相同。對此分析認為，爆炸應力波傳播至減振孔時，會產生反射、透射和繞射現象。由于減振孔中空氣的存在，部分應力波能量發生耗散，其余波能以反射波和繞射波的形式繼續在塊體中傳播。當減振孔直徑增大時，一方面，入射波波能耗散增加，使得繞射進入減振孔后方的應力波減少，并且由于繞射路徑的增長，繞射傳播過程中，應力波的衰減增加;另一方面，減振孔間的凈間距減小，使得直接通過孔間介質傳播的應力波減少，從而亦導致應力波衰減增大。

隨著減振孔直徑的增加，隔振率呈明顯的增大趨勢，與試驗結果一致。當減振孔直徑從0.1 m增大到0.4 m時，隔振率從大約10%增大到50%，說明減振孔直徑變化對減振孔減振效果影響顯著。從圖8（c）中可以看出，隔振率大小呈現明顯的兩種分布：位于減振孔前后的測點7-8、11-12、15-16隔振率大小接近，而位于減振孔間巖體前后的測點9-10、13-14隔振率大小接近，前者明顯大于后者。兩者差值隨減振孔直徑增加逐漸增大，在孔徑為0.3 m時達到最大，之后逐漸縮小。這是因為，相比于孔間巖體，減振孔對爆炸應力波的削弱作用更強。并且，與位于孔間巖體前后的測點相比，位于孔前后的測點距減振孔更近，因此，可以獲得更好的隔振效果。但當孔徑增大至0.4 m時，孔間隙變得很小，減振孔可近似看作連續屏障，因此，兩者差值縮小。

2.3 減振孔間距對減振效果的影響

數值模型如圖6（a）所示，Φ為0.1 m，S2為4 m，n為1排。

圖9（a）、（b）、（c）分別為減振孔間距變化時，對應的地下洞室壁測點1-2處PPV、地表測點3-6處PPV、減振孔屏障前后測點7-16處隔振率的變化曲線。從數值模擬結果可以看出，測點1-2處PPV隨著減振孔間距的增加呈緩慢增加趨勢，測點3-6處PPV隨減振孔間距的增加呈明顯增加趨勢，說明隨減振孔間距的增大，減振效果減弱。這是由于減振孔間距變大時，一方面，更多的爆炸應力波可直接從減振孔間的介質傳播到減振孔后方，另一方面，減振孔離散性變大，使得減振孔對爆炸應力波的整體屏蔽作用變差。因此，在實際工程中，建議在爆源和保護區之間布設較密集的減振孔，可以得到更好的減振效果，減小爆炸應力波對保護區的破壞作用。

圖9（c）中，測點7-16處隔振率隨減振孔間距的增加呈明顯的減小趨勢，當減振孔間距由0.2 m增大到0.5 m時，隔振率由大約22%減小到8%。另外，5組測點與減振孔的相對位置會隨減振孔間距的變化而變化。當孔間距較小時（0.2 m或0.3 m），減振孔整體更接近連續屏障，5組測點的隔振率大小相近;當孔間距較大時（0.4 m或0.5 m），減振孔之間離散性增加，5組測點的隔振率大小呈現差異性。

2.4 減振孔距爆源距離對減振效果的影響

數值模型如圖6（a）所示，Φ為0.1 m，S1為0.3 m，n為1排。

減振孔距爆源距離變化對減振孔減振效果影響變化規律見圖10所示。由圖10（a）可以看到，隨著減振孔與爆源之間距離的增大，測點1-2（鄰近洞室左側壁和洞頂）處PPV呈現緩慢增大的變化規律。由圖10（b）可以看到，隨著減振孔與爆源之間距離的增大，測點3-6（地表）處PPV呈現明顯增大的變化規律，說明減振孔減振效果隨爆心距的增加而減弱。分析認為，應力波隨著傳播距離的增加而逐漸衰減。當減振孔距離爆源較近時，波的能量較大，減振孔的存在可以有效阻隔大部分爆炸應力波的傳播;而當減振孔距爆源較遠時，相較于前者，更多的波能在經過減振孔之前，發生傳播耗散，減振孔對爆炸應力波的屏蔽作用變小，從而減振效果減弱。因此，在實際工程中，為了增強減振效果，更好地屏蔽爆炸應力波的傳播，建議在距爆源較近的位置布設減振孔。

由圖10（c）可以看到，隨減振孔與爆源之間距離的增大，隔振率呈現明顯減小的變化規律，且當減振孔距離爆源較近時（爆心距為4～5 m），隔振率減幅大，爆心距的變化對減振孔的減振效果影響程度大;而當減振孔距離爆源較遠時（爆心距為5～7 m），隔振率減幅變小，爆心距的變化對減振孔的減振效果影響程度變小。另外，位于減振孔前后的3組測點7-8、11-12、15-16的隔振率大于減振孔間巖體前后的2組測點9-10、13-14的隔振率，且差距隨爆心距的增加而減小，當爆心距大于5 m時，5組測點的隔振率大小基本相同。

2.5 減振孔排數對減振效果的影響

數值模型如圖6（d）所示，Φ為0.1 m，S1為0.3 m，S2為4 m，S3為0.3 m。

減振孔排數變化對減振孔減振效果影響的數值模擬結果見圖11所示。由圖11（a）、（b）可以看出，測點1-2（鄰近洞室左側壁和洞頂）處PPV和測點3-6（地表）處PPV均隨減振孔排數的增大明顯減小。這主要是由于爆炸應力波通過更多排數的減振孔時會發生多次透、反射和繞射，傳播方向被干擾，傳播路徑增長，使得爆炸應力波能量發生更多的耗散，從而導致通過減振孔的應力波幅值減弱。因此，在實際工程中，在一定范圍內合理布設多排減振孔可以起到更好的減振效果。

由圖11（c）可以看出，隨著減振孔排數的增加，這5組測點的隔振率均明顯增大。孔前后和孔間巖體前后得到的結果依然表現出較明顯的差異，總體上前者的隔振率更大，隔振效果更好。隨孔排數增加，兩者變化幅度一致，差值沒有明顯變化。

此外，對比圖8（a）、圖9（a）、圖10（a）、圖11（a）可以看出，當減振孔直徑變化時，鄰近洞室PPV變化比較明顯;而當減振孔間距、孔距爆源距離及孔排數變化時，鄰近洞室PPV雖有增大或減小的趨勢，但變化不明顯。而對比圖8（c）、圖9（c）、圖10（c）、圖11（c）可以看出，減振孔參數變化時，減振孔屏障前后測點的隔振率大小均有明顯變化。其中，減振孔直徑變化時，隔振率的變化幅度最大。鄰近洞室PPV變化不明顯的原因可能是由于鄰近洞室離爆源較遠，爆炸應力波到達此處時，大部分能量已經耗散，所以，鄰近洞室的振動絕對值變化不明顯。相比而言，減振孔屏障前后測點的隔振率變化比較明顯，主要是由于減振孔參數變化對靠近減振孔測點的PPV影響較大。

3 結論

針對地下巖體工程中的爆破減振問題，通過有限元數值模擬軟件AUTODYN建立了含減振孔的數值模型，經爆破試驗驗證后，系統地分析了在爆炸作用下地下洞室不同參數減振孔對鄰近洞室和地表的動態響應特性的影響，探究了減振孔參數對減振孔減振效果的影響規律。主要結論如下：

1）減振孔尺寸、位置和排列的變化均會對減振孔減振效果產生影響。隨著減振孔直徑和孔排數的增加，或者隨著孔間距和孔距爆源距離的減小，地下鄰近洞室和地表的PPV減小，減振孔屏障前后測點的隔振率明顯增大，減振效果更顯著。

2）減振孔直徑變化時，不同位置測點的隔振率大小呈現如下兩種分布：減振孔前后測點的隔振率大小接近，而減振孔間巖體前后測點的隔振率大小接近。且前者明顯大于后者，兩者差值隨孔直徑的增大先增大后縮小。

3）減振孔距爆源距離變化時，減振孔前后的3組測點隔振率大于減振孔間巖體前后的2組測點隔振率，兩者差值隨爆心距的增加迅速縮小。減振孔排數變化時，隔振率大小呈現同樣的分布，但兩者差值隨孔排數增加沒有明顯變化。

4）減振孔間距變化會影響測點與減振孔的相對位置。當孔間距較小時，5組測點的隔振率大小相近;而當孔間距較大時，5組測點的隔振率大小呈現差異性。

5）減振孔直徑對鄰近洞室振動和隔振率的影響程度最大，而孔間距、減振孔距爆源距離和孔排數等的影響程度相對較小。

參考文獻：

[1] 楊仁樹， 丁晨曦， 王雁冰， 等. 爆炸應力波與爆生氣體對被爆介質作用效應研究[J]. 巖石力學與工程學報， 2016， 35（Sup2）： 3501-3506.

YANG R S， DING C X， WANG Y B， et al. Action-effect study of medium under loading of explosion stress wave and explosion gas [J]. Chinese Journal of Rock Mechanics and Engineering， 2016， 35（Sup2）： 3501-3506.（in Chinese）

[2] DENG X F， CHEN S G， ZHU J B， et al. UDEC-AUTODYN hybrid modeling of a large-scale underground explosion test [J]. Rock Mechanics and Rock Engineering， 2015， 48（2）： 737-747.

[3] ZHU J B， LI Y S， WU S Y， et al. Decoupled explosion in an underground opening and dynamic responses of surrounding rock masses and structures and induced ground motions：A FEM-DEM numerical study [J]. Tunnelling and Underground Space Technology， 2018， 82： 442-454.

[4] DENG X F， ZHU J B， CHEN S G， et al. Numerical study on tunnel damage subject to blast-induced shock wave in jointed rock masses [J].Tunnelling and Underground Space Technology， 2014， 43： 88-100.

[5] LI H B， XIA X， LI J C， et al. Rock damage control in bedrock blasting excavation for a nuclear power plant [J]. International Journal of Rock Mechanics and Mining Sciences， 2011， 48（2）： 210-218.

[6] LI X B， LI C J， CAO W Z， et al. Dynamic stress concentration and energy evolution of deep-buried tunnels under blasting loads [J]. International Journal of Rock Mechanics and Mining Sciences， 2018， 104： 131-146.

[7] DENG X F， ZHU J B， CHEN S G， et al. Some fundamental issues and verification of 3DEC in modeling wave propagation in jointed rock masses [J]. Rock Mechanics and Rock Engineering， 2012， 45（5）： 943-951.

[8] ZHU J B， DENG X F， ZHAO X B， et al. A numerical study on wave transmission across multiple intersecting joint sets in rock masses with UDEC [J]. Rock Mechanics and Rock Engineering， 2013， 46（6）： 1429-1442.

[9] LI J C， LI N N， CHAI S B， et al. Analytical study of ground motion caused by seismic wave propagation across faulted rock masses [J]. International Journal for Numerical and Analytical Methods in Geomechanics， 2018， 42（1）： 95-109.

[10] 孫崔源， 薛里， 劉世波， 等. 深基坑減震孔降震效果試驗和數值模擬研究[J]. 兵工學報， 2017， 38（Sup1）： 237-242.

SUN C Y， XUE L， LIU S B， et al. Experimental and numerical simulation on vibration reduction effect of damping holes for deep foundation pits [J]. Acta Armamentarii， 2017， 38（Sup1）： 237-242.（in Chinese）

[11] 孫振， 謝事臻， 張憲堂， 等. 減震孔在小間距隧道爆破開挖中的降振效應[J]. 山東科技大學學報（自然科學版）， 2019， 38（1）： 25-31， 39.

SUN Z， XIE S Z， ZHANG X T， et al. Vibration reduction effect of damping holes in blasting excavation of tunnels with small spacing [J]. Journal of Shandong University of Science and Technology （Natural Science）， 2019， 38（1）： 25-31， 39.（in Chinese）

[12] 黃富強， 張大巖， 夏祥. 核電站核島基礎爆破開挖減震控制[J]. 土工基礎， 2012， 26（1）： 30-32.

HUANG F Q， ZHANG D Y， XIA X. Blasting vibration control of foundation pit excavation in nuclear unit at a nuclear power plant [J]. Soil Engineering and Foundation， 2012， 26（1）： 30-32.（in Chinese）

[13] 徐亮. 多排密集鉆孔爆破減振理論及設計方法研究[D]. 福建 泉州： 華僑大學， 2014.

XU L.Study of the blasting damping theory and the design method of multi-row drilling holes [D]. Quanzhou， Fujian： Huaqiao University， 2014.（in Chinese）

[14] UYSAL O， ERARSLAN K， CEBI M A， et al. Effect of barrier holes on blast induced vibration [J]. International Journal of Rock Mechanics and Mining Sciences， 2008， 45（5）： 712-719.

[15] 鄭文富. 減振孔參數對降振效率影響研究[D]. 湖南 湘潭： 湖南科技大學， 2016.

ZHENG W F.Study of damping hole parameters influence on blasting vibration reduction [D]. Xiangtan， Hunan： Hunan University of Science and Technology， 2016.（in Chinese）

[16] LEI Z， KANG Q， ZHAO M S， et al. Numerical simulation and experimental study on vibration-decreasing function of the barrier holes [J]. Applied Mechanics and Materials， 2014， 602/603/604/605： 53-59.

[17] 楊典光. 減震孔對爆破地震的隔振效果研究[D]. 長沙： 長沙理工大學， 2015.

YANG D G.The study of vibration isolation mechanism and effect of damping hole in blasting seismic [D]. Changsha： Changsha University of Science & Technology， 2015.（in Chinese）

[18] BIAN K F， LI H B， LIU Y Q. Numerical simulation for the damping effect of empty holes on blasting vibration [J]. Advanced Materials Research， 2014， 936： 1490-1495.

[19] PARK D. Reduction of blast-induced vibration intunnelling using barrier holes and air-deck [C]// Proceedings of the 12th ISRM Congress， October 18-21， 2011， Beijing， International Society for Rock Mechanics and Rock Engineering， 2011： 17-26.

[20] 惠峰， 李志龍， 徐全軍， 等. 空孔對爆破地震波減震作用的數值模擬分析[J]. 爆破， 2012， 29（4）： 58-61， 116.

HUI F， LI Z L， XU Q J， et al. Numerical simulation on effect of empty hole on damping of blasting vibration [J]. Blasting， 2012， 29（4）： 58-61， 116.（in Chinese）

[21] AUTODYN User Manual [M]. Version 6.1， Century dynamics， 2005.

[22] 高祥濤， 解北京. 巖石中爆炸應變波信號測試與分析[J]. 爆破， 2013， 30（2）： 41-46.

GAO X T， XIE B J. Analysis and measurement of blasting strain wave signal in rock [J]. Blasting， 2013， 30（2）： 41-46.（in Chinese）

[23] LIU F， GUO Z R， LV H， et al. Test and analysis of blast wave in mortar test block [J]. International Journal of Rock Mechanics and Mining Sciences， 2018， 108： 80-85.

[24] 項大林， 榮吉利， 李健， 等. 基于KHT程序的RDX基含鋁炸藥JWL狀態方程參數預測研究[J]. 北京理工大學學報， 2013， 33（3）： 239-243.

XIANG D L， RONG J L， LI J， et al. JWL equation of state parameters prediction of RDX-based aluminized explosive based on KHT code [J]. Transactions of Beijing Institute of Technology， 2013， 33（3）： 239-243.（in Chinese）

（編輯 王秀玲）