Study on the energy transfer efficiency of explosive blasting with different coupling medium
-
摘要: 不耦合装药爆破可有效降低孔壁峰值压力，改善爆破效果。针对不同耦合介质爆破时爆炸能量的传递问题，考虑岩体应变率效应，理论分析了爆炸作用下岩体变形及破坏特征，得到不同耦合介质爆破时理论爆炸能量的传递效率，并结合数值模拟研究了岩体性质、炸药类别及不耦合装药系数对不同耦合介质爆破时爆炸能量传递效率差异的影响。研究结果表明，爆炸能量传递效率与耦合介质相关，装药结构相同时，水耦合爆破比空气耦合爆破爆炸能量传递效率高；不同耦合介质爆破爆炸能量传递效率的差异受爆破岩体、炸药性质及不耦合装药系数等因素影响；装药系数相同耦合介质不同的爆破，岩体强度越高，不同耦合介质爆炸能量传递效率差别越大；岩体性质相同时，不同耦合介质爆破间能量传递效率差异随不耦合装药系数的增大而增大，对于乳化炸药在粉砂岩中的爆破，不耦合装药系数由1.28增至3.44时，水耦合爆破传递入周围岩体的能量由空气耦合爆破的1.45倍增至6.52倍。研究结果可对优化爆破设计、改善爆炸能量分布提供参考。Abstract: Decoupled charge blasting can effectively reduce the peak pressure of hole wall and improve the blasting effect. Aiming at the issues associated with the explosive energy transferred into rock mass with different coupling medium, the deformation and failure characteristics of rock mass under explosion were analyzed theoretically with the consideration of strain rate effect of the rock mass, and the theoretical energy transfer efficiency of blasting with different coupling medium was obtained. Combined with the numerical simulations, the effects of rock mass properties, explosive categories and decoupled charge coefficient on the energy transfer efficiency of blasting with different coupling medium were studied. The results show that the energy transferred into rock mass from explosive of decoupling charge blasting is related to the coupling medium, the energy transfer efficiency of water coupling blasting is higher than that of air coupling blasting for the same charge structure and same blasting medium. When the charge structure and blasting rock mass are the same but the coupling medium is different, the energy transferred into the rock mass from explosive will be different, depending on blasting rock mass, explosive categories and decoupling charge coefficient. When blasting with the same decoupling charge coefficient but different coupling medium, the higher the rock mass strength, the greater the difference of energy transfer efficiency between different coupling medium. For blasting with the same rock mass properties and same charge structure, the difference of energy transfer efficiency between air coupling blasting and water coupling blasting increases with the increase of decoupling charge coefficient. For emulsion explosive exploding in siltstone, when the decoupling charge coefficient increases from 1.28 to 3.44, the energy transferred from water coupling blasting to surrounding rock mass increases from 1.45 to 6.52 times of air coupling blasting. The research results are of great reference significance for optimizing blasting design, improving explosion energy distribution and increasing the explosion energy utilization rate.
-
Key words:
- decoupled charge /
- coupling medium /
- explosion energy /
- transfer efficiency /
- rock mass
-
钻爆法是岩体开挖的主要手段，如何实现爆炸能量的高效利用和爆破危害效应的有效控制是工程爆破技术的研究热点[1]。工程实践和研究表明[2-3]，径向不耦合装药爆破能有效降低孔壁峰值压力，减少粉碎区范围，改善爆破效果，提高爆炸能量利用率，较耦合装药爆破优势显著。而空气和水是工程爆破中炮孔壁和炸药间常见的两种耦合介质，因其物理性质及动力学特性不同，对爆炸能量传递及爆破效果的改善也存在差异。
目前，爆炸应力波动态作用以及高温高压爆生气体准静态作用联合破岩是岩体爆破破碎主流观点，即爆炸能量分为冲击波能及爆生气体能，炸药能量通过做功传递给岩体，冲击波能主要消耗在粉碎岩体、形成初始裂隙以及诱发岩体振动方面，爆生气体能则耗散于爆腔及初始裂隙的扩张、破碎岩体的抛掷[4]等方面。Brinkman[5]通过套管试验，观察了模型在爆炸冲击波或爆生气体单独作用下的破坏特征，证实了爆破岩体破裂是由爆破冲击波能和爆生气体能综合作用产生的。Livingston[6]通过爆破漏斗试验，总结提出了炸药在岩体中爆破后能量的几种存在形式，主要包括弹性变形能、冲击破裂能、破碎能、抛掷能及空气冲击波能，定性分析了爆炸能量的分布方向。Sanchidrián等[7]通过单孔爆破试验计算了岩体内炸药爆破后的能量分布，计算结果表明，地震波能量占炸药总能量的1%～3%，破碎能占炸药总能量的2%～6%，动能占炸药总能量的3%～21%。Hong等[8]通过霍普金森压杆试验得出：爆炸产生的冲击波能量撞击炮孔壁后分散为反射能、透射能及岩体破碎能，且岩体破碎能的分配随总入射能的增加近似线性提高。
炸药爆炸传递至岩体中能量的多少与炸药性能、岩体性质和装药结构有关，为提升炸药的能量利用率，学者们对炸药（耦合装药和不耦合装药）与矿岩的匹配进行了广泛的研究，提出了波阻抗匹配[9]、全过程匹配[10]和能量匹配[11]等观点。采用不耦合装药爆破时，耦合介质不同，爆破效果和传至岩体能量的多少也有显著差异。陈世海等[12]从理论上对水介质耦合装药爆破与炸药耦合装药的破岩效能进行了讨论，认为水介质耦合装药爆破更能提高炸药的能量利用率；宗琦等[13]推导了水介质耦合条件下的孔壁峰值压力和破岩范围，认为炮孔水介质耦合装药比空气介质耦合装药更能提高爆炸能量利用率，增强破岩能力。Jang[14]通过数值模拟和现场试验分析了水垫层对岩体爆破的影响，发现爆炸能量与水相互作用在孔底产生了均匀连续的压力转换，使爆块尺寸分布更加均匀，没有巨石产生。顾文彬等[15]、Xia等[16]分析了装药结构对能量传递的影响，发现爆炸能量的传递受到炸药与岩石的波阻抗、不耦合装药系数及炮孔内耦合介质的压缩程度影响，不同装药结构能量利用存在差别。但是目前的分析多停留在试验及定性分析阶段，缺乏定量化指标。
本文中将基于波动方程、率相关的岩体破坏分区模型，求解不同耦合介质爆破时爆炸能量的理论传递效率，并通过数值模拟分析耦合介质对爆炸能量传递效率的影响。
1. 不耦合装药爆破孔壁爆炸荷载峰值计算
不耦合装药爆破时，炸药爆炸产生冲击波，冲击波在传播过程中先后与耦合介质及岩体发生碰撞，并在介质交界面发生透反射。在不同介质的交界面上，界面两侧应存在应力和位移连续，同时任意冲击波均满足质量、动量及能量守恒方程[17]：
ρ(D−u)=ρ0(D−u0) (1) p−p0=ρ0(D−u0)(u−u0) (2) e−e0=(p+p0)2(1ρ0−1ρ)+Q (3) 式中：下标0表示波前参数，其余为波后参数；Q为介质的比能，若介质不释放能量，则Q=0；p、ρ、u、D分别为冲击波的压力、密度、质点速度及传播速度，e为比内能。
1.1 空气耦合爆破
空气耦合爆破时，炮孔壁入射压力与爆生气体的膨胀过程有关，当不耦合系数较小时，爆生气体在炮孔内仅经历等熵膨胀，此时孔壁入射压力为[9]：
pi=ρeD22(k+1)(dcdb)2k (4) 当不耦合系数较大时，爆生气体在炮孔内经历等熵膨胀及绝热膨胀，此时孔壁入射压力为[9]：
pi=(pwpk)γ/kpk(dcdb)2γ (5) 式中：pi为炮孔壁入射压力，pw为平均爆轰压力，pk为临界压力，ρe、D分别为炸药密度及爆速，dc、db分别为装药直径及炮孔直径，k、γ为绝热指数。
根据爆生气体与岩体交界面上的位移连续，可得[3]：
ut=ui+ur (6) 式中：ui为入射波孔壁质点速度，ut和ur分别为透射波及反射波孔壁质点速度。
根据波的质量和动量守恒可得任意介质中波的质点速度，即：
ui=√pi(1ρi0−1ρi) (7) ur=ui−√(pr−pi)(1ρi−1ρr) (8) ut=√pt(1ρt0−1ρt) (9) 式中：ρi0为入射波波前密度，即空气初始密度；ρt0为透射波波前密度，即原岩密度。
多方气体的状态方程为[6]：
e=pτγ−1 (10) 由气体状态方程结合反射波守恒方程可得：
ρrρi=(γ+1)pr+(γ−1)pi(γ−1)pr+(γ+1)pi (11) 由交界面上应力连续条件可得：
pt=pi+pr (12) 取岩石的状态方程为：
pt=A′[(ρtρt0)k′−1] (13) 式中：A′为常数，与岩体性质有关；k′为岩石的等熵指数。
联立式（6）～（13）可得空气耦合装药孔壁透射峰值压力。
1.2 水耦合爆破
水耦合爆破时，爆生产物首先与水发生碰撞，在水中激起冲击波，在爆生产物与水的交界面上同样满足位移及应力连续条件。由瞬态爆轰的爆轰波波阵面参数基本公式可得[3]：
ui=Dk+1 (14) 根据爆轰产物多方曲线状态方程及反射波动量守恒方程可得：
ur=2kDk2−1[1−(ptpi)(k−1)/2k] (15) 将式（14）～（15）代入式（6）可得：
ut=Dk+1{1+2kk−1[1−(ptpi)(k−1)/2k]} (16) 水的状态方程可表示为：
pt=A″[(ρtρw0)k″−1] (17) 式中：A′′为常数，k′′为水的等熵指数，ρw0为水密度。
联立式（9）、（16）～（17）可得水耦合爆破水中冲击波初始峰值压力。
冲击波在水中传播伴随着能量耗散及几何衰减，延长药包装药水中超压随相对距离的衰减关系为[18]：
Δp={A1ˉr2.491≤ˉr≤2A2ˉr1.452＜ˉr≤5A3ˉr0.635＜ˉr≤240 (18) 式中：A1为水中冲击波初始峰值压力，
ˉr 为相对距离，ˉr =r/rc，其中r为距炮孔中心距离，rc为药卷半径。由式（18）可得水岩界面入射压力，水耦合爆破孔壁透射峰值压力推导过程与空气耦合装药相近，联立式（6）、（12）～（13）、（17）可得水耦合爆破孔壁透射峰值压力。
2. 基于率相关的岩体爆破破坏范围计算
在爆炸荷载作用下，炮孔周围岩体发生破坏，炮孔近区岩体由近及远分为粉碎区、裂隙区及弹性振动区,但这种岩体破坏分区模型忽略了裂隙区的环向承载力，与实际分布存在较大偏差。冷振东等[19]根据现场试验观察提出，在粉碎区及裂隙区间，存在以受压及剪切破坏为主的破碎区，破碎区内冲击波衰减较缓。岩体在爆炸荷载作用下处于强动力响应状态，炮孔近区岩体应变率为：
˙ε=urbα(rrb)−(α+1)/rb (19) 式中：urb为孔壁峰值质点速度；rb为炮孔半径；α为压力衰减指数，α=2±μ/(1−μ)，在不同破坏区域α的取值不同。
岩体中任意一点的有效应力为[20]：
σi=Bpt(rrb)−α (20) 式中：
B=1√2[(1+λ)2+(1+λ2)−2μ(1−μ)(1−λ)2]1/2 ，μ为岩体泊松比，λ为侧向压力系数。根据有效应力强度准则，当σi满足下式，岩体将会发生破坏：
σi={κσcd˙ε1/3˙ε≥˙εc0σtd(1+(˙ε˙εt0)γ)˙ε≥˙εt0 (21) 式中：
κ 为静单轴抗压强度缩小系数，σcd为岩石静态抗压强度，σtd为静态抗拉强度，˙εc0 、˙εt0 分别为抗压强度及抗拉强度临界应变率。联立式（19）～（21），可得不耦合装药爆破岩体破坏范围，粉碎区范围：
r1=(κσcdBpt)31−2α1(urbα1rb)11−2α1rb (22) 破碎区范围：
r2=(κσcdBσr1)31−2α2(ur1α2r1)11−2α2r1 (23) 裂隙区范围：
r3=(σtdBpt(r1/rb)−α1(1+(˙εθ˙εt0)γ))1/−α3r1 (24) 式中：σr1、ur1为粉碎区边界位置质点应力及振动速度，
˙εθ=ur1α3(r3/r1)−(α3+1)/r1 ； α1、 α2、 α3分别为粉碎区、破碎区及裂隙区应力波衰减指数。3. 不耦合装药爆破传递入岩体内能量计算
3.1 冲击波能量
3.1.1 冲击波爆腔膨胀
在爆炸冲击荷载作用下，孔壁岩石将会压缩，爆腔发生膨胀。爆腔膨胀在岩石压碎过程中持续存在，在破碎区边缘，冲击波引起的爆腔膨胀同步结束，粉碎区及破碎区均消耗了大量冲击波能量。冲击波扩胀爆腔的过程中，扩胀范围内岩体质量守恒，因此有：
(r22−r2b)ρm=∫r2rk2ρrdr (25) 式中：rk为冲击波扩胀后爆腔的半径， r2为岩体破碎区范围，ρm为原岩密度，ρ为冲击波波阵面上的岩石瞬态密度。
爆腔扩胀过程中，其波阵面后岩石密度变化很小，因此冲击波阵面上岩石瞬态密度ρ可用压缩后孔壁处岩石密度ρr代替，根据波阵面上岩石质量守恒以及岩体Hugoniot方程可得：
ρr=a+burba+(b−1)urbρm (26) 式中：a、b为与岩性相关的参数。
联立式（25）～（26）可得，冲击波爆腔扩胀半径为：
rk=[r22−(r22−r2b)ρm/ρr]1/2 (27) 由此可得冲击波消耗在岩体爆腔扩胀上的能量为：
w1=∫rkrb2πrptˉr−αdr (28) 3.1.2 冲击波裂隙扩张
冲击波在岩体内逐渐衰减，在破碎区边缘基本衰减为应力波，此时径向应力已经小于岩体抗压强度，但环向应力仍大于岩体抗拉强度，会引起岩体开裂。岩体中r处的应力波切向应力为：
σθ=λptˉr−α (29) 此时裂隙长度为a，若σθ>[σtd]，裂隙将进一步扩展，假定裂隙扩展到a+δa时，切向应力衰减到抗拉强度，裂隙停止扩展；此时裂隙在切向应力作用下产生了切向位移，切向力对岩体做功，在此过程中切向力做功近似为：
δw=2∫δa012σθuθda (30) 由断裂力学可得[21]：
σθ=k2√2πa (31) uθ=k1Em√δa−a2π(1+μ)(f+1) (32) 式中：f=(3−4μ)，k1、k2为长度不同的裂隙端部的应力强度因子，为简化计算可近似认为裂隙扩展过程中强度因子保持不变，Em为岩石弹性模量。
联立式（30）～（32）可得裂隙扩展过程中切向应力所做的功为：
w2=n∫r3rbδwdr=∫r3rbnk22Em(1−μ2)dr (33) 式中：n为径向主裂隙条数，根据相关试验文献，n=4～12。
3.1.3 应力波引起的弹性变形
随着应力波的进一步衰减，裂隙区外应力波不再造成岩体破坏，仅引起弹性振动，柱状坐标系下单位体积岩体发生弹性变形需要的能量为[22]：
ΔE=12(σrεr+σθεθ) (34) 由变形区r→∞得到，应力波引起岩石弹性变形所做的功为：
w3=∫+∞r32πrΔEdr=π(1+λ2)Emp2tr2b2(α−1) (35) 3.2 准静态爆生气体能量
3.2.1 准静态荷载扩腔
在冲击荷载作用之后，爆生气体膨胀压缩耦合介质，以准静态荷载的形式作用在孔壁上，爆腔进一步扩大，当爆腔内准静态压力p等于围岩压力ps时，准静态爆腔扩胀结束。不考虑大气压力及围岩自重影响，岩体中围岩压力近似为三轴应力状态下岩石的屈服强度σs，其中σs=(ρm
c2p /σcd)1/4σcd，其中cp为岩体纵波波速。考虑到爆生气体在空气与水中的膨胀差异，分别对空气耦合及水耦合爆破准静态扩腔过程进行分析。
（1）空气耦合装药
空气耦合爆破时，爆生气体的膨胀过程与距起爆中心的距离有关，根据式（4）～（5）可得炮孔内任一点爆生气体的膨胀压力为：
p={ρeD22(k+1)(rrb)2kp≥pk(pwpk)γkpk(rrb)2γp<pk (36) 根据准静态爆腔扩胀条件及炮孔内爆生气体膨胀压力（式（36））可得准静态荷载爆腔扩张半径
r′k ，进而可得爆生气体膨胀准静态做功为：w4=∫r′kr22πrpdr (37) （2）水耦合装药
冲击波扩腔结束后，爆生气体膨胀压缩水至充满爆腔，设此时膨胀压缩过程结束后水的径向压缩量为δ，由水的质量守恒可得压缩后水的密度为：
ρ′w=(r2b−r2c)ρwr2b−(rc+δ)2 (38) 根据流体力学理论，水在压缩过程中满足：
dp=−2kwrr2b−r2dr (39) 联立式（36）、（38）～（39）可得冲击波扩腔结束后爆生气体半径
r′2 。准静态爆生气体扩腔完成后，由准静态爆腔扩胀条件可得炮孔内压力平衡方程：
p=σs=ph (40) 式中：ph为不耦合介质水中压力。
联立式（36）～（40）可得准静态爆腔扩张结束后爆生气体半径
r′k ，此过程爆生气体膨胀准静态做功为：w4=∫r′kr′22πrpdr (41) 3.2.2 准静态荷载作用下裂纹扩展
研究表明[23]，岩体在应力波作用下产生初始裂隙，随后裂隙在爆生气体准静态荷载作用下进一步扩展。考虑到裂隙长度通常远大于炮孔尺寸，因此可将炮孔尺寸近似为裂隙长度的一部分，此时裂纹尖端的应力强度因子为[24]：
K1=21−N[L(t)+rbπ]1/2∫L(t)+rb0p(x,t)−σ{[L(t)+rb]2−x2}1/2dx (42) 式中：L(t)为准静态扩展裂纹长度，N为裂纹尖端初始损伤，p(x,t)为裂隙内准静态荷载分布，σ为远场应力。
当满足下式时裂纹扩展终止：
K1=Ka (43) 假定最终裂纹长度为L，对应炮孔内压力为p，联立式（36）、（42）～（43）可求得考虑初始损伤情况下最终裂纹长度L，进而求得爆生气体准静态作用裂纹扩展做功为：
w5=n∫Lr3(1−μ2)EK21dr (44) 3.3 爆炸传递入岩体总能量
炸药爆炸传递入岩体内的能量最终消耗于粉碎区、破碎区、裂隙区及弹性振动区，因此各部分能量总和即为爆炸传递入岩体内能量，它们与炸药能量的比值即为不耦合装药爆破爆炸能量传递效率，即：
η=w1+w2+w3+w4+w5E (45) 式中：E为单位长度炸药的总能量。
以乳化炸药、花岗岩为例，计算典型装药结构的不同耦合介质爆破时爆炸能量传递的效率。乳化炸药（ρe=1 300 kg/m3，D=4 000 m/s）的体积爆热为4.192 GJ/m3，花岗岩力学参数如表1所示[25]。
表 1 花岗岩力学参数Table 1. Mechanical parameters of granite岩石类型 密度/（kg·m−3） 纵波波速/（m·s−1） 抗压强度/MPa 抗拉强度/MPa 弹性模量/GPa 泊松比 花岗岩 2700 5500 150 15 68 0.24 不同耦合介质爆炸能量传递效率计算结果如表2所示。根据计算结果可知，水介质相较于空气介质具有更好的能量传递效果。同时，相较于空气，水的流动黏度大、压缩性低、冲击波衰减缓慢，造成的粉碎区范围大于空气耦合，耗散于粉碎区的能量更多，计算结果与相关试验结论接近。
表 2 不同耦合介质爆炸能量传递效率Table 2. Energy transfer efficiency of blasting explosion with different coupling mediumK 耦合介质 w1/E0 w2/E0 w3/E0 w4/E0 w5/E0 η 1.28 空气 0.123 0.036 0.041 0.004 0.077 0.281 水 0.249 0.083 0.075 0.079 0.145 0.631 上述研究建立在爆炸荷载垂直作用孔壁的假定基础上，未能完全真实考虑爆炸荷载斜入射以及反射带来的影响，同时爆炸荷载在炮孔耦合介质中的衰减也较难得出理论解，需要采用部分经验公式，这些假定及经验公式为理论上定量分析不同耦合介质爆破的具体能量差异带来困难，可能造成分析结果出现一定误差。下面结合高精度数值模拟方法，分析不同耦合介质爆破时爆炸能量传递效率的差异及其影响因素。
4. 不同耦合介质爆破爆炸能量传递效率数值模拟
4.1 计算模型及参数
为进一步研究不同耦合介质爆破时爆炸能量传递效率的差异以及爆破介质、不耦合系数及炸药种类对能量传递效率差异的影响，基于实际开挖过程中岩体种类的分布、常用爆破参数及工业炸药类型，选择粉砂岩、石灰岩及花岗岩3类性质差异明显的岩体，分别代表软、硬及坚硬类岩体介质，同时选用6种常用不耦合装药系数，以及乳化炸药和多孔粒状铵油炸药2种典型工业炸药开展相关研究。
数值计算采用非线性动力有限元软件LS-DYNA，计算模型如图1所示，为1/4炮孔模型，尺寸为1.5 m×1.5 m×3 m，为确保能较准确地模拟出真实的爆炸效应，将计算模型中炸药和耦合介质单元的网格尺寸控制在2 mm左右，邻近岩石单元尺寸也与炸药/耦合介质单元近似相同，模型网格数量约38万，不同工况时保证模型尺寸及网格划分基本一致，计算时间步与模型最小网格匹配；模型边界包括无反射边界、对称边界及自由边界；采用流固耦合算法模拟爆炸荷载的冲击作用。
炸药采用JWL状态方程进行模拟，其爆炸过程中的压力和内能及相对体积之间的关系为：
p=A(1−ωR1V)e−R1V+B(1−ωR2V)e−R2V+ωE0V (46) 式中：p为爆轰压力，V为爆轰产物的相对体积，其余为方程常数。计算炸药相关参数取值见表3。
表 3 计算炸药参数Table 3. Table of charge parameters岩体材料采用双线性随动硬化塑性模型，该模型可以较好地考虑岩体破坏过程中的应变硬化，应变率敏感性以及毁伤失效等特征，与理论分析中基于应变率的岩体破坏分区模型相对应，数值模拟中采用的岩石力学参数如表4所示，与理论计算参数对应，均取自岩石力学参数手册推荐范围[25]。
表 4 计算岩体力学参数表Table 4. Table of mechanical parameters of rock mass used in the simulations岩石类型 密度/
（kg·m−3）弹性模量/
GPa泊松比 屈服应力/
MPa切线模量/
GPa粉砂岩 2170 6.7 0.25 39.2 0.6 石灰岩 2600 32.5 0.25 72.9 3.0 花岗岩 2700 68.0 0.24 150.0 7.0 LS-DYNA中空气采用mat_null材料模型结合*EOS_Linear_Polynomial状态方程描述；水采用mat_null材料模型结合*EOS_Grüneisen状态方程描述，相关参数取值见文献[19]。
4.2 计算结果
图2中给出了部分典型工况下的爆炸能量时程曲线。由图2可知，相同装药结构下，相较于空气介质，水介质具有更好的传能效果，且爆炸能量增长速率显著快于空气介质，这与水的压缩性显著低于空气、黏滞度大于空气有关，与耦合介质的物理性质差异相匹配；随着时间的推移，爆炸传至岩体内的能量趋于稳定，水耦合爆炸能量传递效率显著高于空气耦合装药，数值模拟结果与前述理论分析结果基本一致。
图 2 典型工况下爆炸能量时程曲线Figure 2. Time history curve of explosion energy under typical working conditions不同工况下不同耦合介质爆破时爆炸能量传递效率比值如图3所示。
图 3 不同装药结构下不同耦合介质能量传递效率比值Figure 3. The ratio of energy transfer efficiency for different coupling medium under different charge structures由计算结果可知，不同耦合介质间能量传递效率的差值并非定值，受到爆破介质、炸药性质以及不耦合装药系数的影响。在相同装药条件下，花岗岩水耦合爆破比空气耦合爆破时爆炸能量传递效率的提升更加明显，石灰岩、粉砂岩次之，即岩体强度越高能量传递效率差别越大；炸药性质对不同耦合介质爆炸能量传递效率的差值也存在一定的影响，采用乳化炸药爆破时不同耦合介质能量传递效率差异比采用铵油炸药爆破小；当采用同种炸药爆破相同岩体介质时，能量传递效率的差异随不耦合系数的增大逐渐增大，在不耦合系数较小时，不同耦合介质能量传递的效率基本相同，极端情况接近耦合装药情况下，两者能量传递效率基本相当；当不耦合系数较大时，水耦合爆破比空气耦合爆破具有更好的能量传递效果；以乳化炸药在粉砂岩中起爆为例，不耦合系数由1.28增至3.44时，水耦合爆破传递入周围岩体的能量由空气耦合爆破的1.45倍增至6.52倍。
5. 结 论
通过理论及数值模拟分析，研究了不同耦合介质爆破时爆炸能量传递的效率及其影响因素，得到以下主要结论：
（1）基于波动方程及界面连续条件，求解了不耦合装药爆破孔壁爆炸荷载峰值，进而得到了基于率相关的岩体破坏区范围，并进行了分区能量计算，得到不同耦合介质爆破时爆炸能量的理论传递效率；
（2）装药结构相同时，水耦合爆破比空气耦合爆破减少了爆炸荷载在耦合介质中的衰减，爆炸荷载峰值压力高、作用时间长，提高了爆炸能量传递的效率；
（3）不同耦合介质爆破时爆炸能量传递效率的差值并非定值，受到爆破介质、炸药性质以及不耦合装药系数的影响，相同装药条件下，硬岩（花岗岩、石灰岩）比软岩（粉砂岩）不同耦合介质爆破时能量传递效率的差别更大；
（4）爆破工况相同时，不同耦合介质爆破时爆炸能量传递效率的差异随不耦合系数的增大逐渐增大，以乳化炸药在粉砂岩中起爆为例，不耦合系数由1.28增至3.44时，水耦合爆破传递入周围岩体的能量由空气耦合爆破的1.45倍增至6.52倍。
-
图 2 典型工况下爆炸能量时程曲线
Figure 2. Time history curve of explosion energy under typical working conditions
图 3 不同装药结构下不同耦合介质能量传递效率比值
Figure 3. The ratio of energy transfer efficiency for different coupling medium under different charge structures
表 1 花岗岩力学参数
Table 1. Mechanical parameters of granite
岩石类型 密度/（kg·m−3） 纵波波速/（m·s−1） 抗压强度/MPa 抗拉强度/MPa 弹性模量/GPa 泊松比 花岗岩 2700 5500 150 15 68 0.24 
下载: 导出CSV
表 2 不同耦合介质爆炸能量传递效率
Table 2. Energy transfer efficiency of blasting explosion with different coupling medium
K 耦合介质 w1/E0 w2/E0 w3/E0 w4/E0 w5/E0 η 1.28 空气 0.123 0.036 0.041 0.004 0.077 0.281 水 0.249 0.083 0.075 0.079 0.145 0.631
下载: 导出CSV
表 4 计算岩体力学参数表
Table 4. Table of mechanical parameters of rock mass used in the simulations
岩石类型 密度/
（kg·m−3）弹性模量/
GPa泊松比 屈服应力/
MPa切线模量/
GPa粉砂岩 2170 6.7 0.25 39.2 0.6 石灰岩 2600 32.5 0.25 72.9 3.0 花岗岩 2700 68.0 0.24 150.0 7.0
下载: 导出CSV
-
[1] 中国爆破行业协会. 中国爆破行业中长期科学与技术发展规划(2016～2025年) [M]. 北京: 中国爆破行业协会, 2015: 2−5. [2] 王志亮, 李永池. 工程爆破中径向水不耦合系数效应数值仿真 [J]. 岩土力学, 2005, 26(12): 1926–1930. DOI: 10.3969/j.issn.1000-7598.2005.12.012.WANG Z L, LI Y C. Numerical simulation on effects of radial water-decoupling coefficient in engineering blast [J]. Rock and Soil Mechanics, 2005, 26(12): 1926–1930. DOI: 10.3969/j.issn.1000-7598.2005.12.012. [3] 王伟, 李小春. 不耦合装药下爆炸应力波传播规律的试验研究 [J]. 岩土力学, 2010, 31(6): 1723–1728. DOI: 10.3969/j.issn.1000-7598.2010.06.008.WANG W, LI X C. Experimental study of propagation law of explosive stress wave under condition of decouple charge [J]. Rock and Soil Mechanics, 2010, 31(6): 1723–1728. DOI: 10.3969/j.issn.1000-7598.2010.06.008. [4] 王文龙. 钻眼爆破[M]. 北京: 煤炭工业出版社, 1984: 196−209. [5] BRINKMANN J R. Separating shock wave and gas expansion breakage mechanisms [C] // Proceedings of the Second International Symposium on Rock Fragmentation by Blasting. Keystone, Colorado, International Organising Committee of Rock Fragmentation by Blasting, 1989: 6−15. [6] LIVINGSTON C W. Fundamental concepts of rock failure [J]. Quarterly of the Colorado school of mines, 1956, 51(3): 1–11. [7] SANCHIDRIÁN J A, SEGARRA P, LÓPEZ L M. Energy components in rock blasting [J]. International Journal of Rock Mechanics and Mining Sciences, 2007, 44(1): 130–147. DOI: 10.1016/j.ijrmms.2006.05.002. [8] HONG L, ZHOU Z L, YIN T B, et al. Energy consumption in rock fragmentation at intermediate strain rate [J]. Journal of Central South University of Technology, 2009, 16(4): 677–682. DOI: 10.1007/s11771-009-0112-5. [9] 杨善元. 岩石爆破动力学基础[M]. 北京: 煤炭工业出版社, 1993: 143−144. [10] 郭子庭, 吴从师. 炸药与岩石的全过程匹配 [J]. 矿冶工程, 1993, 13(3): 11–15.GUO Z T, WU C S. The exeplosive-rock matching throughout the blasting process [J]. Mining and Metallurgical Engineering, 1993, 13(3): 11–15. [11] 赖应得. 论炸药和岩石的能量匹配 [J]. 工程爆破, 1995, 1(2): 22–26.LAI Y D. On energy match between explosive and rock [J]. Engineering Blasting, 1995, 1(2): 22–26. [12] 陈士海, 崔新壮. 水偶合装药与全偶合装药爆破效能讨论 [J]. 爆破, 1998, 15(2): 10–13.CHEN S H, CUI X Z. Study on blasting efficacy of water coupling and explosive coupling charges [J]. Blasting, 1998, 15(2): 10–13. [13] 宗琦, 田立, 汪海波. 水介质不耦合装药爆破岩石破坏范围的研究和应用 [J]. 爆破, 1998, 29(2): 42–46. DOI: 10.3963/j.issn.1001-487X.2012.02.011.ZONG Q, TIAN L, WANG H B. Study and application on rock damage range by blasting with water-decoupled charge [J]. Blasting, 1998, 29(2): 42–46. DOI: 10.3963/j.issn.1001-487X.2012.02.011. [14] JANG H, HANDEL D, KO Y, et al. Effects of water deck on rock blasting performance [J]. International Journal of Rock Mechanics and Mining Sciences, 2018, 112: 77–83. DOI: 10.1016/j.ijrmms.2018.09.006. [15] 顾文彬, 王振雄, 陈江海, 等. 装药结构对爆破震动能量传递及爆破效果影响研究 [J]. 振动与冲击, 2016, 35(2): 207–211. DOI: 10.13465/j.cnki.jvs.2016.02.035.GU W B, WANG Z X, CHEN J H, et al. Influence of charge structure on the energy transfer of blasting vibration and explosive effect [J]. Journal of Vibration and Shock, 2016, 35(2): 207–211. DOI: 10.13465/j.cnki.jvs.2016.02.035. [16] XIA W J, LU W B, LI R Z, et al. Effect of water-decked blasting on rock fragmentation energy [J]. Shock and Vibration, 2020, 2020: 8194801. DOI: 10.1155/2020/8194801. [17] 王礼立. 应力波基础M]. 2版. 北京: 国防工业出版社, 2005: 25−29. [18] 王伟, 李小春, 石露, 等. 深层岩体松动爆破中不耦合装药效应的探讨 [J]. 岩土力学, 2008, 29(10): 2837–2842. DOI: 10.16285/j.rsm.2008.10.009.WANG W, LI X C, SHI L, et al. Discussion on decoupled charge loosening blasting in deep rock mass [J]. Rock and Soil Mechanics, 2008, 29(10): 2837–2842. DOI: 10.16285/j.rsm.2008.10.009. [19] 冷振东. 岩石爆破中爆炸能量的释放与传输机制[D]. 武汉: 武汉大学, 2017: 53−60. [20] 戴俊. 柱状装药爆破的岩石压碎圈与裂隙圈计算 [J]. 辽宁工程技术大学学报(自然科学版), 2001, 20(2): 144–147. DOI: 10.3969/j.issn.1008-0562.2001.02.005.DAI J. Calculation of radii of the broken and cracked areas in rock by a long charge explosion [J]. Journal of Liaoning Technical University (Natural Science), 2001, 20(2): 144–147. DOI: 10.3969/j.issn.1008-0562.2001.02.005. [21] 布洛克. 工程断裂力学基础[M]. 王克仁, 何明元, 高桦, 译. 北京: 科学出版社, 1980: 6−22. [22] 宗琦, 杨吕俊. 岩石中爆炸冲击波能量分布规律初探 [J]. 爆破, 1999, 16(2): 1–6.ZONG Q, YANG L J. Shock energy distribution of column charge in rock [J]. Blasting, 1999, 16(2): 1–6. [23] 卢文波, 陶振宁. 爆生气体驱动的裂纹扩展速度研究 [J]. 爆炸与冲击, 1994, 14(3): 264–268.LU W B, TAO Z Y. A study of fracture propagation velocity driven by gases of explosion products [J]. Explosion and Shock Waves, 1994, 14(3): 264–268. [24] 徐颖, 孟益平, 宗琦, 等. 断层带爆炸裂隙区范围及裂纹扩展长度的研究 [J]. 岩土力学, 2002, 23(1): 81–84. DOI: 10.3969/j.issn.1000-7598.2002.01.018.XU Y, MENG Y P, ZONG Q, et al. Study on range of cranny and length of fissure expansion in fault zone [J]. Rock and Soil Mechanics, 2002, 23(1): 81–84. DOI: 10.3969/j.issn.1000-7598.2002.01.018. [25] 水利水电科学研究院. 岩石力学参数手册[M]. 北京: 水利电力出版社, 1991: 429−434. [26] 夏祥, 李海波, 李俊如, 等. 岭澳核电站二期工程基岩爆破安全阈值分析 [J]. 岩土力学, 2008, 29(11): 2945–2951, 2956. DOI: 10.3969/j.issn.1000-7598.2008.11.010.XIA X, LI H B, LI J R, et al. Research on vibration safety threshold for rock under blasting excavation [J]. Rock and Soil Mechanics, 2008, 29(11): 2945–2951, 2956. DOI: 10.3969/j.issn.1000-7598.2008.11.010. [27] 刘军. 岩体在冲击载荷作用下的各向异性损伤模型及其应用 [J]. 岩石力学与工程学报, 2004, 23(4): 635–640. DOI: 10.3321/j.issn:1000-6915.2004.04.020.LIU J. Anisotropic damage model and its application to rock materials under impact load [J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(4): 635–640. DOI: 10.3321/j.issn:1000-6915.2004.04.020. 期刊类型引用(20)
1. 马衍坤，张正锐，张纪辉. 水耦合作用下高压空气冲击煤岩破裂规律与声发射响应试验研究. 振动与冲击. 2025(02): 1-10+29 . 
百度学术2. 谭铭. 炸药性能对爆破地震波传播与能量衰减规律影响研究. 中国矿业. 2024(01): 193-199 . 百度学术3. 马泗洲，刘科伟，杨家彩，李旭东. 不耦合装药下岩石爆破块体尺寸的分布特征. 爆炸与冲击. 2024(04): 122-140 . 
本站查看4. 杨建华，彭超，叶志伟，冷振东，魏彬. 深部岩体爆破冲击波能量分布特征. 兵工学报. 2024(06): 1735-1746 . 百度学术5. 李昊明，李祥龙，王建国，胡启文，武永博. 不同耦合介质下岩石损伤特性数值模拟. 有色金属工程. 2024(07): 121-129 . 百度学术6. 叶志伟，陈明，杨建华，姚池，张小波，周创兵. 隧道爆破炮孔堵塞结构运动规律与长度优化研究. 中国公路学报. 2024(08): 204-215 . 百度学术7. 余诚. 环保严控区高陡灰岩路堑抑尘控爆施工技术研究. 工程与建设. 2024(05): 1138-1140+1149 . 百度学术8. 杨凯. 双向荷载下不耦合装药对石灰岩预裂爆破的影响研究. 煤. 2023(02): 104-108 . 百度学术9. 周贤舜，张学民，武朝光，胡涛，陈鑫磊，段亚. 基于低碳减排的隧道水封爆破优化效果研究. 铁道科学与工程学报. 2023(03): 996-1007 . 百度学术10. 付晓强，麻岩，俞缙，戴良玉，黄凌君. 隧道爆破振动信号时频谱增强优化分析. 矿业科学学报. 2023(03): 348-356 . 百度学术11. 冷振东，高启栋，卢文波，陈明，周桂松，范勇. 岩石钻孔爆破能量调控理论与应用技术研究进展. 金属矿山. 2023(05): 64-76 . 百度学术12. 张学富，李宁杰，周杰，韩风雷，张建. 异形孔水包爆破围岩粉碎范围研究. 振动与冲击. 2023(11): 131-138 . 百度学术13. 张宪堂，董国庆，余辉，张经双，周红敏，朱宝合，殷振. 围压下空孔直径对直眼掏槽爆破振动的影响. 山东科技大学学报(自然科学版). 2023(03): 44-52 . 百度学术14. 滕树强，吴亮，代颖，刘汉龙. 岷江龙溪口水下岩体钻爆的爆破能量特征研究. 水运工程. 2023(10): 104-109 . 百度学术15. 刘万通，徐振洋，张久洋，王雪松，刘鑫. 孔底间隔介质对岩体损伤规律研究. 有色金属工程. 2023(10): 82-94 . 百度学术16. 孙得志，郭俊庆，杨永康，赵国飞，谢二伟，宋选民，刘亚龙. 水介质预裂爆破顶板裂纹的扩展. 工程爆破. 2023(06): 100-109 . 百度学术17. 李康贵，夏鹏，米元桃，陈明，李桐. 水耦合强抛掷台阶爆破技术研究. 爆破. 2022(02): 23-29 . 百度学术18. 周桂松，钟冬望. 绿色爆破的爆炸能量转化机制. 金属矿山. 2022(07): 35-41 . 百度学术19. 谢先启，黄小武，姚颖康，何理，伍岳. 露天深孔台阶精细爆破技术研究进展. 金属矿山. 2022(07): 7-18 . 百度学术20. 杨赛群，李洪伟，吴立辉，刘伟. 影响光面爆破效果和空孔处能量分布的试验研究. 工程爆破. 2022(06): 58-65 . 百度学术其他类型引用(20)
-
下载:
百度学术
计量
- 文章访问数: 895
- HTML全文浏览量: 433
- PDF下载量: 188
- 被引次数: 40