岩石钻孔爆破粉碎区计算模型的改进

冷振东; 卢文波; 陈明; 严鹏; 胡英国

doi:10.11883/1001-1455(2015)01-0101-07

岩石钻孔爆破粉碎区计算模型的改进

doi: 10.11883/1001-1455(2015)01-0101-07

冷振东^{1, 2},
卢文波^{1, 2, ,},
陈明^{1, 2},
严鹏^{1, 2},
胡英国^{1, 2}

1.
武汉大学水资源与水电工程科学国家重点实验室，湖北武汉 430072
2.
武汉大学水工岩石力学教育部重点实验室，湖北武汉 430072

基金项目: 国家自然科学基金项目(51125037，51279135)；国家重点基础研究发展计划(973计划)项目(2011CB013501)；高等学校博士学科点专项科研基金项目(20110141110026)

详细信息

作者简介:
冷振东(1989—), 男, 博士研究生

通讯作者:
卢文波, wblu@whu.edu.cn

中图分类号: O383.1;TV554
计量
- 文章访问数: 3275
- HTML全文浏览量: 339
- PDF下载量: 489
- 被引次数: 94
出版历程
- 收稿日期: 2013-05-03
- 修回日期: 2013-07-25
- 刊出日期: 2015-01-25

Improved calculation model for the size of crushed zone around blasthole

Leng Zhen-dong^{1, 2},
Lu Wen-bo^{1, 2
, ,},
Chen Ming^{1, 2},
Yan Peng^{1, 2},
Hu Ying-guo^{1, 2}

1.
State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, Hubei, China
2.
Key Laboratory of Rock Mechanics in Hydraulic Structural Engineering Ministry of Education, Wuhan University, Wuhan 430072, Hubei, China

摘要

摘要: 为了研究炮孔周围岩石的破坏机理，准确预测粉碎区的范围，提出了一种计算钻孔爆破粉碎区范围的改进模型。该四分区模型考虑了破裂区内侧的环向压应力和炮孔空腔膨胀的影响，假定粉碎区为丧失了内聚力但仍具有内摩擦力的散体介质。采用弹塑性力学理论推导了柱状装药起爆条件下的岩石钻孔爆破粉碎区半径公式。计算结果表明，岩石钻孔爆破粉碎区范围通常为1.2~5.0倍炮孔半径，不同种类岩石的粉碎区范围差别很大。与其他计算模型相比，本模型的计算结果与实验数据更吻合。
- 爆炸力学 /
- 计算模型 /
- 环向应力 /
- 粉碎区 /
- 改进 /
- 空腔膨胀
Abstract: To explore the breakage mechanism of rocks around the blasthole and to accurately predict the size of the crushed zone in drilling and blasting, an improved calculation model for calculation the size of the crushed zone was presented. The four-region model was established with hoop compressive stress in the inner part of fractured zone and cavity expansion effect taken into account. The material in the crushed zone is assumed to be granular medium without cohesion but with internal friction. On this basis, the formula for the crushed zone radius in drilling blasting is derived by elastic-plastic mechanics theory. The analysis shows that the size of the crushed zone ranges from 1.2 to 5.0 times of the blasthole radius, and there are notable discrepancies among rock types. Compared with other models, the improved model is in better agreement with the experiment data.
- mechanics of explosion /
- calculation model /
- hoop stress /
- crushed zone /
- improvement /
- cavity expansion

HTML全文

近十年来, 随着现代化矿山的建设和西部水电的开发, 钻孔爆破作为石方开挖的最常用技术得到了快速发展。大量的研究资料表明, 炸药破碎岩石的能量仅占炸药总能量的20%~30%。N.V.Melnikov^[1]指出, 在传统装药结构的爆破中, 有不少于50%的爆炸能量浪费在粉碎圈和破裂圈内侧部分的过度破碎上, 而且极度粉碎的岩石很容易堵塞破裂区形成的裂纹通道, 阻碍爆炸气体向裂纹中的进一步扩散, 影响了气体的“气楔”作用^[2], 减少了岩石的破坏范围。同时, 在粉碎区内比表面积大的颗粒会吸收爆生气体的大量热能, 降低了炸药能量的有效利用率。因此, 如何控制岩石的爆破粉碎区范围对于提高炸药能量的有效利用率、降低炸药单耗和工程成本具有重要意义。

对岩石在炸药作用下的破坏范围已有较深入的研究^[3~6]。目前粉碎区计算模型的结果差异很大, 主要有以下几种:G.Szuladzinski^[7]模型假设炮孔周围的岩石为弹性体, 粗略地假定炸药的有效能量为炸药爆热的三分之二, 不能反映装药结构对粉碎区范围的影响; N.Djordjevic^[8]基于Griffith强度准则的计算模型只适用于脆性岩石; A.A.Il’yushin^[9]基于Mohr-Coulumb准则的计算模型认为粉碎区岩石仍具有黏聚力, 但是A.Vovk等^[10]在石灰岩和混凝土中的爆破实验结果表明, Il’yushin模型的结果偏大。另一方面Il’yushin公式的推导过程中, 将爆腔膨胀过程中的气体绝热指数取为常数, 因此该公式不适用于不耦合系数较大的情况, 从而使公式具有很大的局限性。S.S.Kanchibotla等^[11]和S.Esen等^[12]基于工程经验和实验统计的计算公式简单, 但是结果离散性较大; 国内学者戴俊^[13]基于Mises屈服准则的计算模型等。由于炮孔空腔的膨胀, 作用在炮孔上的压力会降低^[14], 上述模型除Il’yushin外均未考虑炮孔空腔膨胀对炮孔压力的影响。

本文中提出一种计算钻孔爆破粉碎区范围的改进模型, 由于破裂区是连接粉碎区和弹性区的约束, 此区域不可能被径向裂纹完全破坏, 因此将破裂区划分为两部分, 内侧部分(破裂Ⅰ区)的介质为塑性破坏, 环向压应力不为零; 外侧部分(破裂Ⅱ区)介质受到径向裂缝破坏, 丧失了环向承载力。假定粉碎区为丧失了黏聚力的散体介质, 但是破碎颗粒之间仍然具有内摩擦力。在此基础上, 考虑炮孔空腔膨胀对炮孔压力降低的影响, 推导柱状装药起爆条件下的岩石钻孔爆破粉碎区半径公式。

1. 爆破破坏分布及计算条件

1.1 爆破破坏分布

炸药爆炸后, 炮孔由里往外岩石依次承受剧烈的爆炸冲击波、应力波和地震波作用, 岩石介质的连续性发生改变, 呈现出不同的破碎和损伤状态。炮孔壁在向外扩展一定距离后逐渐稳定下来, 最终形成膨胀空腔。根据周围岩石的破坏程度, 炮孔周围的岩石可以划分为不同的区域, 不同的学者对破坏分区的定义各不相同^[15]。现有的计算模型通常把爆破作用的最终影响范围划分为粉碎区、径向开裂区和弹性变形区3部分, 如图 1(a)所示。

图 1 传统模型和改进模型的岩石钻孔爆破破坏分区示意图

Figure 1. Image of the damage zones surrounding a blasthole under blasting in traditional and improved calculation model

下载: 全尺寸图片幻灯片

现有的模型^[6-7,9]认为粉碎区和弹性区之间的整个区域被径向裂纹完全破坏, 因而岩石只能传递径向应力, 没有环向承载力, 即σ_θ=0, 从而简化物理过程。然而在实际过程中, 破裂区是连接粉碎区和弹性区的约束, 不可能造成径向完全的破坏, 特别是破裂区的内侧部分岩石受到极高的径向压应力, 由于泊松效应, 必然会受到周围岩石的约束, 因此环向压应力作用明显, 不能忽略其作用。因此, 我们把破裂区划分为2个部分, 内侧部分(破裂Ⅰ区)的介质为塑性破坏, 需要考虑环向应力的影响:σ_θ≠0;外侧部分(破裂Ⅱ区)介质受到径向裂缝破坏, 丧失了环向承载力:σ_θ=0, 改进模型能更好地反映炮孔周围岩石的实际破坏情况。

图 1(b)为改进模型中的钻孔爆破破坏分区示意图, 在该四分区模型中, 各破坏分区的边界定义如下:粉碎区:a(t)≤r ≤b_*(t); 破裂Ⅰ区:b_*(t) < r ≤b_Ⅰ(t); 破裂Ⅱ区:b_Ⅰ(t) < r ≤b_Ⅱ(t); 弹性变形区:b_Ⅱ(t) < r ≤∞。以上各式中, a(t)为膨胀空腔的半径, b_*(t)为粉碎区的半径, b_Ⅰ(t)为破裂Ⅰ区半径, b_Ⅱ(t)为破裂Ⅱ区的半径。

1.2 计算条件

在岩石介质中有一无限长的圆柱形空腔, 受到一个内部沿轴向均布的爆炸荷载的作用, 作如下假设:(1)圆柱形空腔沿轴向无限延伸, 可将问题视为轴对称平面应变问题; (2)粉碎区岩石为各向同性、不可压缩且丧失黏聚力的散体介质, 但是颗粒之间仍然具有内摩擦力; (3)爆生气体的膨胀过程为绝热膨胀, 忽略进入岩石裂隙的爆生气体的体积。

2. 公式推导

2.1 弹性变形区

利用柱坐标系描述问题, 在弹性区的应力分布为:

$\sigma_{r}=\sigma_{0}\left[1-\left(\frac{b_{\mathrm{II}}}{r}\right)^{2}\right]+\sigma_{r=b \mathbb{I}}\left(\frac{b_{\mathrm{II}}}{r}\right)^{2}, \quad \sigma_{\theta}=\sigma_{0}\left[1-\left(\frac{b_{\mathrm{II}}}{r}\right)^{2}\right]-\sigma_{r=b}\left(\frac{b_{\mathrm{II}}}{r}\right)^{2}$

(1)

式中:σ_{r=b_Ⅱ}为弹性区的内边界(r=b_Ⅱ)处的径向应力, σ₀为岩石初始应力。在弹性变形区和破裂Ⅱ区的交界面上, 环向应力达到岩石的抗拉强度:σ_θ=-[σ_t], 故由式(1)得σ_{r=b_Ⅱ}=[σ_t]+2σ₀。

弹性区的位移

$u=\frac{1+\nu}{E} \frac{b_{\mathbb{I I}}^{2}}{r}\left(\left[\sigma_{\mathrm{t}}\right]+\sigma_{0}\right)$

(2)

2.2 破裂Ⅱ区

破裂Ⅱ区的特征是介质受到裂缝破坏, 丧失了环向承载力, 但是径向仍为弹性, 类似于径向柱杆, 主要把破裂Ⅰ区传来的压力过渡到弹性区介质中去。破裂Ⅱ区满足σ_θ=0, 平衡微分方程可以简化为:

$\frac{\mathrm{d} \sigma_{r}}{\mathrm{~d} r}+\frac{\sigma_{r}}{r}=0$

(3)

在破裂Ⅱ区外边界上有σ_r=[σ_t]+2σ₀, 内边界上有σ_r=[σ_c], 其中[σ_c]为岩石单轴抗压强度。由此得到破裂Ⅱ区的径向应力

$\sigma_{r}=\frac{\left(\left[\sigma_{\mathrm{t}}\right]+2 \sigma_{0}\right) b_{\mathrm{II}}}{r}=\frac{\left[\sigma_{\mathrm{c}}\right] b_{\mathrm{I}}}{r}$

(4)

由式(2)知r=b_Ⅱ时, 则破裂Ⅱ区内边界r=b_Ⅰ处的位移

$u_{b_{\mathrm{I}}}(t)=\frac{1+\nu}{E} b_{\mathrm{I}}\left[\sigma_{\mathrm{c}}\right]\left(\frac{\left[\sigma_{\mathrm{t}}\right]+2 \sigma_{0}}{\left[\sigma_{\mathrm{t}}\right]+\sigma_{0}}-(1-v) \ln \frac{\left[\sigma_{\mathrm{t}}\right]+2 \sigma_{0}}{\left[\sigma_{\mathrm{c}}\right]}\right)$

(5)

2.3 破裂Ⅰ区

破裂Ⅰ区为塑性破坏区, 在破裂Ⅰ区产生大量裂缝, 导致介质体积的膨胀, 因此需要考虑该区域岩石的剪胀作用。采用非关联流动法则来描述破裂Ⅰ区岩石的剪胀特性:

$h \varepsilon_{r}^{\mathrm{p}}+\varepsilon_{\theta}^{\mathrm{p}}=0$

(6)

式中:h为破裂Ⅰ区岩石的剪胀率。

可以求得破裂Ⅰ区的位移

$u(t)=\frac{1+\nu}{E}\left[\sigma_{\mathrm{c}}\right]\left[\frac{(1-h)(1-\nu)}{1+h} r+b_{\mathrm{I}}^{1+1 / h} L r^{-1 / h}\right]$

(7)

为了表述简便, 记

把Mohr-Coulomb屈服准则代入平衡微分方程得到破裂Ⅰ区径向应力的分布:

$\sigma_{r}=\frac{1+\sin \theta}{2 \sin \theta}\left[\sigma_{\mathrm{c}}\right]\left(\frac{b_{\mathrm{I}}}{r}\right)^{\frac{2 \sin \theta}{1+\sin \theta}}-\frac{1-\sin \theta}{2 \sin \theta}\left[\sigma_{\mathrm{c}}\right]$

(8)

2.4 粉碎区

在高温高压的爆炸气体的作用下, 装药附近的岩石受到强烈的压缩剪切作用, 结构被完全破坏, 形成紧挨着空腔壁的粉末化区域。该区域的岩石可以视为各向同性、不可压缩并且没有粘聚力的散体介质, 但是破碎颗粒之间仍然具有内摩擦力。在粉碎区采用没有凝聚内成分的Mohr-Coulomb强度准则:

$\left(\sigma_{r}-\sigma_{\theta}\right) / 2=\left[-\left(\sigma_{r}+\sigma_{\theta}\right) \sin \theta\right] / 2$

(9)

在粉碎区外边界满足σ_r=σ_s, σ_s为多向应力条件下的岩体动抗压强度, $\sigma_{\mathrm{s}}=\left[\sigma_{\mathrm{c}}\right]\left(\frac{\rho_{\mathrm{m}} c_{\mathrm{p}}^{2}}{\left[\sigma_{\mathrm{c}}\right]}\right)^{\frac{1}{4}}$ ^[16]; ρ_m为岩石密度; c_p是岩石介质中的纵波速度。

将(9)式代入平衡微分方程中解得粉碎区的径向应力

$\sigma_{r}=\sigma_{\mathrm{s}}\left(b_{*} / r\right)^{\frac{2 \sin \theta}{1+\sin \theta}}$

(10)

由界面上径向应力的连续性条件, b_Ⅰ=ξ b_*, 其中ξ为比例系数, 。

由粉碎区的不可压缩条件, 解得粉碎区的位移

$u(t)=\left(\xi^{1+1 / h} L+\frac{1-h}{1+h}\right) \frac{1+\nu}{E}\left[\sigma_{\mathrm{c}}\right] b_{*}^{2} r^{-1}$

(11)

u(t)对b_*(t)求导得

$\frac{\partial u}{\partial b_{*}}=2\left(\xi^{1+1 / h} L+\frac{1-h}{1+h}\right) \frac{1+\nu}{E}\left[\sigma_{\mathrm{c}}\right] \frac{b_{*}}{r}$

(12)

当≪1时, 式子成立, 其中v(r)为粉碎区内某一点的质点速度。

则在膨胀空腔的壁面(r=a(t))上有

$a \mathrm{~d} a=2 \frac{1+\nu}{E}\left(\xi^{1+1 / h} L+\frac{1-h}{1+h}\right)\left[\sigma_{\mathrm{c}}\right] b_{*} \mathrm{~d} b_{*}$

(13)

在初始时刻(t=0时), 粉碎区由膨胀空腔壁面开始产生, 此时有a=b_*=r_b, 其中r_b为炮孔半径。

对方程(13)两边积分, 整理得

$a_{\mathrm{m}}^{2}=2 \frac{1+\nu}{E}\left(\xi^{1+1 / h} L+\frac{1-h}{1+h}\right)\left[\sigma_{\mathrm{c}}\right] b_{* \mathrm{~m}}^{2}+\left\{1-2 \frac{1+\nu}{E}\left(\xi^{1+1 / h} L+\frac{1-h}{1+h}\right)\left[\sigma_{\mathrm{c}}\right]\right\} r_{\mathrm{b}}^{2}$

(14)

式中:b_*m为粉碎区的最大半径, a_m为膨胀空腔的最大半径。

令, 式(14)可以简化成

$\frac{b_{* \mathrm{~m}}}{a_{\mathrm{m}}}=\sqrt{1 / K+(1-1 / K)\left(\frac{a_{\mathrm{m}}}{r_{\mathrm{b}}}\right)^{-2}}$

(15)

在膨胀空腔的后续扩展中, 炮孔压力可以由两阶段的Jonse-Miller绝热方程来确定^[18]:

$p_{m}=\left\{\begin{array}{ll} p_{\mathrm{b}}\left(\frac{a_{\mathrm{m}}}{r_{\mathrm{b}}}\right)^{-2 \gamma_{1}} \\ p_{\mathrm{b}}\left(\frac{a_{\mathrm{m}}}{r_{\mathrm{b}}}\right)^{-2 \gamma_{2}}\left(\frac{r_{\mathrm{b}}}{r_{\mathrm{k}}}\right)^{2\left(\gamma_{1}-\gamma_{2}\right)} & a_{\mathrm{m}} \gt r_{\mathrm{k}} \end{array}\right.$

(16)

式中:p_b为作用在炮孔壁上的初始压力; p_m为膨胀空腔半径达到最大时作用在孔壁上的压力; 两个阶段的绝热指数分别取γ₁=3, γ₂=1.27;r_k为与临界爆腔压力p_k对应的临界爆腔半径, 。

在r=a_m处, 有σ_r=p_m, 即

$p_{\mathrm{m}}=\sigma_{\mathrm{s}}\left(\frac{b_{* \mathrm{~m}}}{a_{\mathrm{m}}}\right)^{\frac{2 \sin \theta}{1+\sin \theta}}$

(17)

联立(15)~(17)得

$\sigma_{\mathrm{s}}\left[1 / K+(1-1 / K)\left(\frac{a_{\mathrm{m}}}{r_{\mathrm{b}}}\right)^{-2}\right]^{\frac{\sin \theta}{1+\sin \theta}}=\left\{\begin{array}{ll} p_{\mathrm{b}}\left(\frac{a_{\mathrm{m}}}{r_{\mathrm{b}}}\right)^{-2 \gamma_{1}} \\ p_{\mathrm{b}}\left(\frac{a_{\mathrm{m}}}{r_{\mathrm{b}}}\right)^{-2 \gamma_{2}}\left(\frac{r_{\mathrm{b}}}{r_{\mathrm{k}}}\right)^{2\left(\gamma_{1}-\gamma_{2}\right)} \end{array} \quad \begin{array}{l} a_{\mathrm{m}} \leqslant r_{\mathrm{k}} \\ a_{\mathrm{m}} \gt r_{\mathrm{k}} \end{array}\right.$

(18)

由式(18)即可求得最大膨胀空腔半径与炮孔半径的比值(a_m/r_b), 带入式(15)可求得粉碎区的范围(b_*m/r_b)。

对于空腔膨胀比较明显的情况, (1-1/K) (a_m/r_b)^-2≈0, 则式(15)可以简化为

$b_{* \mathrm{~m}}=a_{\mathrm{m}} K^{-1 / 2}$

(19)

联立方程(16)~(17)、(19)可得柱状装药起爆条件下的粉碎区半径公式:

$b_{* \mathrm{~m}}=\left\{\begin{array}{ll}r_{\mathrm{b}}\left(\frac{p_{\mathrm{b}}}{\sigma_{\mathrm{s}}} K^{\frac{\sin \theta}{1+\sin \theta}}\right)^{\frac{1}{2 \gamma_{1}}} K^{-1 / 2} & a_{\mathrm{m}} \leqslant r_{\mathrm{k}} \\ r_{\mathrm{b}}\left(\frac{p_{\mathrm{b}}}{\sigma_{\mathrm{s}}} K^{\frac{\sin \theta}{1+\sin \theta}}\right)^{\frac{1}{2 \gamma_{2}}}\left(\frac{r_{\mathrm{b}}}{r_{\mathrm{k}}}\right)^{\frac{\gamma_{1}-\gamma_{2}}{\gamma_{2}}} \mathrm{~K}^{-1 / 2} & a_{\mathrm{m}}>r_{\mathrm{k}}\end{array}\right.$

(20)

由公式(20)可以看出, 粉碎区半径主要受以下因素的影响:(1)岩石特性, 包括动抗压强度[σ_c]、动抗拉强度[σ_t]、弹性模量E、泊松比ν、剪胀特性h和初始应力σ₀; (2)炮孔压力p_b, 炮孔压力p_b又与炸药的密度ρ₀、爆轰速度D以及装药结构有关; (3)炮孔半径r_b。

3. 不同计算模型对比

表 1给出了计算中采用的4种岩石的物理力学参数, 炸药分别采用铵油炸药(ρ₀=0.9 g/cm³, D=3 600 m/s)、乳化炸药(ρ₀=1.05 g/cm³, D=4 100 m/s)和Gurit炸药(ρ₀=1.0 g/cm³, D=2 200 m/s)。

表 1 各种岩石的物理力学参数

Table 1. Physico-mechnical properties of different rocks

岩石类型	[σ_c]/MPa	[σ_t]/MPa	E/GPa	θ/(°)	ρ_m/(kg·m^-3)	ν
石英岩	250	20	70	48	2 670	0.26
花岗岩	150	15	50	50	2 670	0.24
石灰岩	60	6	35	43	2 500	0.23
粉砂岩	40	3	30	41	2 600	0.22

下载: 导出CSV

| 显示表格

表 2中列出了各种计算模型在不同岩石不同炸药类型条件下的粉碎区范围。同种岩石采用相同的装药结构时, 铵油炸药产生的粉碎区比乳化炸药的小。当炮孔压力值一定时, 不同种类的岩石的粉碎区范围差别很大, 对于强度较大的硬岩, 粉碎区半径不到2倍炮孔半径; 而对于强度较小的软岩, 粉碎区半径可达3~5倍甚至更大的炮孔半径。绝大多数学者认为在工程爆破中粉碎区的半径不会超过3~5倍炮孔半径^[18], 这与本文的计算相符合。本文中列举的4种岩石的计算结果表明, 各种计算模型得到的粉碎区的范围差别较大, 本文计算结果与S.Esen和戴俊的结果较为接近, 介于两者之间。Il’yushin以及Kanchibotla的计算结果最大, 引起这种显著差异的原因之一是各种模型对于粉碎区的定义不同, Il’yushin和Kanchibotla模型认为发生过塑性破坏、裂隙网状分布的破裂Ⅰ区也属于粉碎区。

表 2 不同粉碎区范围计算模型的计算结果对比

Table 2. Comparison of crushing zone models in full scale blasting conditions

岩石类型	装药结构	炸药类型	b_*m/r_b
岩石类型	装药结构	炸药类型	Il’yushin 模型	Szuladzinski 模型	Djordjevic 模型	Kanchibotla 模型	Esen 模型	戴俊模型	改进模型
石英岩	耦合	铵油	5.08	1.68	1.74	3.42	1.00	1.54	1.16
石英岩	耦合	乳化	5.45	1.81	2.14	4.20	1.04	1.79	1.31
花岗岩	耦合	铵油	5.81	2.17	2.01	4.41	1.06	1.79	1.36
	耦合	乳化	6.22	2.34	2.48	5.42	1.39	2.08	1.55
	K=3	Gurit	2.91	1.79	1.00	1.76	1.00	1.00	1.00
石灰岩	耦合	铵油	9.04	3.43	3.18	6.97	1.71	2.71	2.18
石灰岩	耦合	乳化	9.69	3.70	3.91	8.58	2.25	3.12	2.52
粉砂岩	耦合	铵油	10.76	4.20	4.50	8.54	2.11	3.19	2.76
粉砂岩	耦合	乳化	11.53	4.53	5.54	10.50	2.77	3.67	2.95
注：K为径向不耦合系数，K=炮孔直径\装药直径。

下载: 导出CSV

| 显示表格

4. 与实验结果的对比

由于爆破后, 岩石被破碎和移动, 钻孔爆破的粉碎区半径很难被直接测得, 因此相关的实验数据十分缺乏。本文中收集了A.Vovk等^[10]在石灰岩和混凝土中的4组爆破实验数据, M.Olsson等^[19]在花岗岩中的1组不耦合装药钻孔爆破数据, S.Slaughter^[20]在煤矿中的3组爆破数据。各种计算模型在相应的爆破条件下的计算结果(b_*m)和实测数据(b'_*m)对比如图 2所示。可见, 与其他模型相比, 改进的计算模型和实验结果能够更好地吻合。

图 2 不同计算模型的粉碎区计算结果b_*m和实验数据

*m对比

Figure 2. Comparison of different models with full scale blasting experiments

下载: 全尺寸图片幻灯片

M.Olsson等^[19]通过在花岗岩中的不耦合装药钻孔爆破实验来测量粉碎区的范围, 采用直径64 mm的炮孔和直径22 mm的Gurit炸药。实验发现炮孔周围几乎观察不到粉碎区。由表 2中不耦合装药一栏的计算结果可知, Il’yushin、Szuladzinski和Kanchibotla这3种模型的粉碎区半径和炮孔半径的比值b_*m/r_b分别为2.91、1.79、1.76, 而本文计算结果为1.00, 即几乎不会产生粉碎区, 与实验结果很吻合。其中, 由于Szuladzinski公式中的粉碎区范围只与炸药密度、单位质量炸药释放的热量及岩石抗压强度有关, 不适用于不耦合装药结构的情况; Il’yushin和Kanchibotla这2种模型在不耦合装药结构时的计算结果误差较大, Djordjevic、S.Esen、戴俊和本文的改进模型可以有效反映不耦合装药对粉碎区半径的影响。

5. 结论与讨论

提出了一种计算钻孔爆破粉碎区范围的改进模型, 认为破裂区内侧部分环向应力不为零, 考虑了炮孔空腔膨胀对炮孔压力降低的影响, 推导了柱状装药起爆条件下的岩石钻孔爆破粉碎区半径公式, 并与实验数据进行了对比, 主要有以下结论:(1)破裂Ⅰ区环向应力的非零假定(σ_θ≠0)使改进模型更能反映炮孔周围岩石的实际受力情况。与其他计算模型相比, 改进模型能与实验数据更好地吻合。(2)基于改进模型推导的粉碎区半径公式可以反映岩石特性、炸药特性、装药结构及炮孔半径对粉碎区范围的综合影响。计算表明, 不同种类岩石的粉碎区范围差别很大, 通常粉碎区半径为1.2~5.0倍炮孔半径。

在实际岩石钻孔爆破中, 岩石的破坏机理非常复杂, 本文仅仅是在基于若干假设的前提下对该现象的讨论, 对于岩石钻孔爆破的破坏分区机制还需要更进一步研究。

图 1 传统模型和改进模型的岩石钻孔爆破破坏分区示意图

Figure 1. Image of the damage zones surrounding a blasthole under blasting in traditional and improved calculation model

下载: 全尺寸图片幻灯片

图 2 不同计算模型的粉碎区计算结果b_*m和实验数据*m对比

Figure 2. Comparison of different models with full scale blasting experiments

下载: 全尺寸图片幻灯片

表 1 各种岩石的物理力学参数

Table 1. Physico-mechnical properties of different rocks

岩石类型	[σ_c]/MPa	[σ_t]/MPa	E/GPa	θ/(°)	ρ_m/(kg·m^-3)	ν
石英岩	250	20	70	48	2 670	0.26
花岗岩	150	15	50	50	2 670	0.24
石灰岩	60	6	35	43	2 500	0.23
粉砂岩	40	3	30	41	2 600	0.22

下载: 导出CSV

表 2 不同粉碎区范围计算模型的计算结果对比

Table 2. Comparison of crushing zone models in full scale blasting conditions

岩石类型	装药结构	炸药类型	b_*m/r_b
岩石类型	装药结构	炸药类型	Il’yushin 模型	Szuladzinski 模型	Djordjevic 模型	Kanchibotla 模型	Esen 模型	戴俊模型	改进模型
石英岩	耦合	铵油	5.08	1.68	1.74	3.42	1.00	1.54	1.16
石英岩	耦合	乳化	5.45	1.81	2.14	4.20	1.04	1.79	1.31
花岗岩	耦合	铵油	5.81	2.17	2.01	4.41	1.06	1.79	1.36
	耦合	乳化	6.22	2.34	2.48	5.42	1.39	2.08	1.55
	K=3	Gurit	2.91	1.79	1.00	1.76	1.00	1.00	1.00
石灰岩	耦合	铵油	9.04	3.43	3.18	6.97	1.71	2.71	2.18
石灰岩	耦合	乳化	9.69	3.70	3.91	8.58	2.25	3.12	2.52
粉砂岩	耦合	铵油	10.76	4.20	4.50	8.54	2.11	3.19	2.76
粉砂岩	耦合	乳化	11.53	4.53	5.54	10.50	2.77	3.67	2.95
注：K为径向不耦合系数，K=炮孔直径\装药直径。

下载: 导出CSV

参考文献(20)

[1]	Melnikov N V. Influence of explosive charge design on results of blasting[C]//International Symposium on Mining Research. Clark G B. London, Pergamon, 1962: 147-155.
[2]	朱红兵, 卢文波, 吴亮.空气间隔装药爆破机理研究[J].岩土力学, 2007, 28(5): 986-989. http://d.wanfangdata.com.cn/Periodical/ytlx200705025 Zhu Hong-bing, Lu Wen-bo, Wu Liang. Research on mechanism of air-decking technique in bench blasting[J]. Chinese Journal of Rock and Soil Mechanics, 2007, 28(5): 986-989. http://d.wanfangdata.com.cn/Periodical/ytlx200705025
[3]	王明洋, 邓宏见, 钱七虎.岩石中侵彻与爆炸作用的近区问题研究[J].岩石力学与工程学报, 2005, 24(16): 2859-2863. http://www.cnki.com.cn/Article/CJFDTotal-YSLX200516009.htm Wang Ming-yang, Deng Hong-jian, Qian Qi-hu. Study on problems of near cavity of penetration and explosion in rock[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(16): 2859-2863. http://www.cnki.com.cn/Article/CJFDTotal-YSLX200516009.htm
[4]	Jimeno E L, Carcedo F J A. Drilling and blasting of rocks[M]. Taylor & Francis Group, 1995.
[5]	Ouchterlony F, Moser P. On the branching-merging mechanism during dynamic crack growth as a major source of fines in rock blasting[C]//The 10th International Symposium on Rock Fragmentation by Blasting. New Delhi, India, 2012: 65-75.
[6]	钱七虎.岩石爆炸动力学的若干进展[J].岩石力学与工程学报, 2009, 28(10): 1945-1968. http://d.wanfangdata.com.cn/Periodical/yslxygcxb200910001 Qian Qi-hu. Some advances in rock blasting dynamics[J]. Chinese Journal of Rock Mechanics and Engineering, 2009, 28(10): 1945-1968. http://d.wanfangdata.com.cn/Periodical/yslxygcxb200910001
[7]	Szuladzinski G. Response of rock medium to explosive borehole pressure[C]//Proceedings of the Fourth International Symposium on Rock Fragmentation by Blasting-Fragblast-4. Vienna, Austria, 1993: 17-23.
[8]	Djordjevic N. A two-component of blast fragmentation[C]//Proceedings of the Sixth International Symposium on Rock Fragmentation by Blasting-Fragblast. Johannesburg, South Africa, 1999: 213-222.
[9]	Hustrulid W A. Blasting principles for open pit mining-Theoretical foundations[M]. Rotterdam: Balkema, 1999.
[10]	Vovk A A, Mikhalyuk A V, Belinskii I V. Development of fracture zones in rocks during camouflet blasting[J]. Journal of Mining Science, 1973, 9(4): 383-387. doi: 10.1007/BF02501169
[11]	Kanchibotla S S, Valery W, Morrell S. Modelling fines in blast fragmentation and its impact on crushing and grinding[C]//Proceedings of Explo'99-A Conference on Rock Breaking. Kalgoorlie, Australia, 1999: 137-181.
[12]	Esen S, Onederra I, Bilgin H A. Modelling the size of the crushed zone around a blast hole[J]. International Journal of Rock Mechanics and Mining Sciences, 2003, 40(4): 485-495. http://www.sciencedirect.com/science/article/pii/S1365160903000182
[13]	戴俊.岩石动力学特性与爆破理论[M].北京: 冶金工业出版社, 2002.
[14]	Hagan T N, Gibson I M. Lower blasthole pressures: A means of reducing costs when blasting rocks of low to moderate strength[J]. International Journal of Mining and Geological Engineering, 1988, 6(1): 1-13. doi: 10.1007/BF00881023
[15]	Ghosh A. Fractal and numerical models of explosive rock fragmentation[D]. Tucson: The University of Arizona, 1990.
[16]	Rakishev B, Rakisheva Z B. Basic characteristics of the stages of rock massif destruction by explosive crushing[C]//Asia-Pacific Symposium on Blasting Techniques. Kunming, China, 2011: 65-69.
[17]	亨利奇.爆炸动力学及其应用[M].熊建国, 等, 译.北京: 科学出版社, 1987.
[18]	杨善元.岩石爆破动力学基础[M].北京: 煤炭工业出版社, 1991.
[19]	Olsson M, Bergqvist I. Crack lengths from explosives in multiple hole blasting[C]//The Fifth International Symposium on Rock Fragmentation by Blasting. Montreal, Quebec, Canada, 1996: 187-191.
[20]	Slaughter S. Investigation of coal fines[R]. Internal Report, Julius Kruttschnitt Mineral Research Centre(JKMRC), Australia, 1991.

施引文献

期刊类型引用(49)

1.	雷涛，康普林，叶海旺，李宁，王其洲. 柱状药包爆破过程中应力波叠加与岩体裂隙分布的方向效应研究. 岩石力学与工程学报. 2024(02): 399-411 . 百度学术
2.	孙鹏昌，杨广栋，卢文波，范勇，孟海利，薛里. 考虑岩体破坏分区的岩石爆破爆炸荷载历程研究. 爆炸与冲击. 2024(03): 171-186 . 本站查看
3.	王洪刚，贾永胜，余浩天，罗鹏，黄炳林，周俊汝. 基于根底控制的深孔台阶爆破超深优选方法. 爆破. 2024(01): 44-50 . 百度学术
4.	叶志伟，陈明，杨建华，姚池，张小波，周创兵. 隧道爆破炮孔堵塞结构运动规律与长度优化研究. 中国公路学报. 2024(08): 204-215 . 百度学术
5.	崔雪姣，李启月，冷振东，姚颖康，周建敏，赵明生. 基于能量传输效率控制的现场混装炸药与岩石匹配方法. 爆破. 2024(03): 9-15 . 百度学术
6.	杨国梁，邹泽华，张赫，李峰. 径向不耦合装药爆破下页岩的动态应变分布及损伤分形特征. 爆破. 2024(03): 26-32+42 . 百度学术
7.	李亮，陈嘉骏，赵炼恒，何可培，胡世红，李华隆. 基于统一强度理论的深埋隧道围岩爆破致裂范围评估（英文）. Journal of Central South University. 2024(07): 2341-2364 . 百度学术
8.	阿比尔的，胡尊镕，傅林，韩亚峰，刘明维. CO_2相变致裂岩石的损伤范围的计算及影响因素研究. 爆破器材. 2024(06): 41-49 . 百度学术
9.	刘树新，庄宇，戴谦君，安帅，曹飞. 盲天井中深孔掏槽爆破布孔优化及围岩损伤数值分析. 金属矿山. 2024(12): 44-51 . 百度学术
10.	蓝亦辉. 临近爆破处既有建筑物振动响应分析. 福建建材. 2023(01): 65-70 . 百度学术
11.	冷振东，高启栋，卢文波，陈明，周桂松，范勇. 岩石钻孔爆破能量调控理论与应用技术研究进展. 金属矿山. 2023(05): 64-76 . 百度学术
12.	楼晓明，陈诗伟，李广斌，牛明远，林日宗，姚炳金. 耦合装药条件下不同孔径孔壁冲击压力的阶段特征. 爆炸与冲击. 2023(08): 168-181 . 本站查看
13.	蒋宏杰，卢文波，王高辉，刘义佳，王洋. 大体积混凝土水下接触爆炸破坏分区特征分析. 爆炸与冲击. 2023(10): 15-29 . 本站查看
14.	章彬彬，程金明，占汪妹，赵东波，廖述能，郑中华. 不耦合装药下炮孔—空孔距离对预裂爆破效果的影响探究. 工程爆破. 2023(05): 79-85 . 百度学术
15.	聂根良，程晓红，周文海，胡才智，凌晓. 冲击载荷作用下炮孔近区节理岩体的动力响应. 矿业研究与开发. 2023(11): 72-80 . 百度学术
16.	冷振东，周桂松，刘令，曹进军，侯国荣. 高原大断面隧道现场混装爆破关键技术研究. 地下空间与工程学报. 2023(S2): 890-900 . 百度学术
17.	李芳涛，胡志平，陈南南，张永辉，安学旭. 爆破荷载作用下隧道围岩裂隙范围计算方法研究. 振动与冲击. 2022(08): 260-269 . 百度学术
18.	周桂松，钟冬望. 绿色爆破的爆炸能量转化机制. 金属矿山. 2022(07): 35-41 . 百度学术
19.	魏东，陈明，卢文波，李康贵，王高辉. 岩体性能变化条件下台阶爆破根底的产生机制研究. 岩土力学. 2022(S1): 490-500 . 百度学术
20.	马力，徐甜新，李克民，孙进步，刘宇，张奇峰，薛飞，靳新宇，李瑞行，王恒荣. 露天煤矿抛掷爆破技术研究现状及发展趋势. 煤炭工程. 2022(09): 134-141 . 百度学术
21.	周文海，胡才智，包娟，郑俊杰，梁瑞. 含节理岩体爆破过程中应力波传播与裂纹扩展的数值研究. 力学学报. 2022(09): 2501-2512 . 百度学术
22.	李桐，陈明，叶志伟，卢文波，魏东，郑祥. 混凝土含水裂隙中爆炸压力传播的模型试验研究. 岩土力学. 2022(S2): 205-213 . 百度学术
23.	孙鹏昌，卢文波，杨招伟，孟海利，薛里. 白鹤滩坝肩边坡开挖爆破损伤预测研究. 水力发电学报. 2022(10): 30-41 . 百度学术
24.	李永祺，梁正召，钱希坤，刘红波. 应力波形对岩石爆生裂纹扩展机制影响的数值模拟. 工程科学学报. 2022(12): 2057-2068 . 百度学术
25.	梁瑞，李生荣，包娟，周文海. 高地应力下岩体的爆破损伤及能量特性. 高压物理学报. 2022(06): 127-138 . 百度学术
26.	周春国，王荣富，刘兴堂，封磊，刘刚. 地下厂房岩壁梁爆破松弛试验及稳定性影响研究. 岩土工程技术. 2022(06): 502-506 . 百度学术
27.	魏东，陈明，叶志伟，卢文波，李桐. 基于应变率相关动力特性的岩体爆破破坏区范围研究. 工程科学与技术. 2021(01): 67-74 . 百度学术
28.	尹岳降，夏文俊，卢文波，刘建成，陈明，于永军. 灰岩矿爆破开采岩粉成因及其控制研究. 工程爆破. 2021(06): 44-51 . 百度学术
29.	梁瑞，周文海，余建平，李珍宝，杜超飞，王敦繁. 冲击载荷作用下岩体拉-压损伤破坏的边坡抛掷爆破模拟. 高压物理学报. 2019(01): 82-91 . 百度学术
30.	费鸿禄，苏强，蒋安俊，洪陈超. 爆破载荷下隧道围岩破坏裂隙范围研究. 爆破器材. 2019(02): 51-56 . 百度学术
31.	梁瑞，俞瑞利，周文海，陈鹏辉，陈必港. 基于LS-DYNA模拟的准直眼掏槽爆破技术研究. 有色金属工程. 2019(08): 101-107 . 百度学术
32.	岳志坤，夏文俊，李福千，陈明，尹岳降，聂攀. 大型停车场岩石基坑快速开挖技术. 建筑施工. 2019(12): 2114-2116 . 百度学术
33.	李俊平，叶浩然，侯先芹. 高应力下硬岩巷道掘进端面钻孔爆破卸压动态模拟. 安全与环境学报. 2018(03): 962-967 . 百度学术
34.	冷振东，刘亮，周旺潇，周桂松，杜华善. 起爆位置对台阶爆破爆堆形态影响的离散元分析. 爆破. 2018(02): 50-55+100 . 百度学术
35.	钟权，冷振东，彭峥，刘放. 节理岩体隧洞开挖爆破损伤特性及爆破方案研究. 长江科学院院报. 2018(02): 89-94 . 百度学术
36.	涂书芳，冷振东. 空气间隔装药爆破损伤特性数值模拟研究. 水电与新能源. 2018(02): 16-20+44 . 百度学术
37.	李启月，吴正宇，黄武林. 直眼掏槽空孔效应的计算模型改进与分析. 采矿与安全工程学报. 2018(05): 925-930 . 百度学术
38.	杨跃宗，邵珠山，熊小锋，米俊峰. 岩石爆破中径向和轴向不耦合装药的对比分析. 爆破. 2018(04): 26-33+146 . 百度学术
39.	费鸿禄，洪陈超. 应力波和爆生气体共同作用下裂隙区范围研究. 爆破. 2017(01): 33-36+107 . 百度学术
40.	冷振东，卢文波，范勇，陈明，严鹏. 侧向起爆条件下的爆炸能量分布及其对破岩效果的影响. 爆炸与冲击. 2017(04): 661-669 . 本站查看
41.	张帆，林从谋，肖绍清，温智捷，殷榕鹏，杨宾. 地铁盾构孤石预爆破破碎范围计算方法. 工程爆破. 2017(01): 34-38 . 百度学术
42.	吴发名，刘勇林，李洪涛，姚强. 基于原生节理统计和爆破裂纹模拟的堆石料块度分布预测. 岩石力学与工程学报. 2017(06): 1341-1352 . 百度学术
43.	曲艳东，孙从煌，章文姣，孔祥清，马石磊. 深孔间隔装药爆破对不同孔壁介质的影响. 工程爆破. 2016(03): 6-14 . 百度学术
44.	吴亮，鲁帅，许锋，蔡路军，曾国伟. 矿岩爆破破碎机理、块度分布与测量技术研究动态. 金属矿山. 2016(07): 47-53 . 百度学术
45.	陈鹏辉，雷涛. 基于ALE算法的单孔台阶爆破数值模拟研究. 有色金属(矿山部分). 2016(01): 72-76 . 百度学术
46.	刘亮，卢文波，陈明，严鹏，王高辉. 钻爆开挖条件下岩体临界破碎状态的损伤阈值统计研究. 岩石力学与工程学报. 2016(06): 1133-1140 . 百度学术
47.	刘亮，郑炳旭，陈明，宋锦泉，王高辉，严鹏. 起爆方式对台阶爆破根底影响的数值模拟分析. 爆破. 2015(03): 49-54+78 . 百度学术
48.	楼晓明，陈鹏辉，周文海. 基于岩粉分段装药的爆破技术研究. 有色金属(矿山部分). 2015(06): 65-69 . 百度学术
49.	Leng Zhendong，Lu Wenbo，Chen Ming，Yan Peng，Hu Yingguo. A new theory of rock-explosive matching based on the reasonable control of crushed zone. Engineering. 2014(06): 32-38 . 必应学术

其他类型引用(45)

资源附件(0)

访问统计

图(2) / 表(2)

计量

文章访问数: 3275
HTML全文浏览量: 339
PDF下载量: 489
被引次数: 94

1. 爆破破坏分布及计算条件
1.1 爆破破坏分布
1.2 计算条件
2. 公式推导
2.1 弹性变形区
2.2 破裂Ⅱ区
2.3 破裂Ⅰ区
2.4 粉碎区
3. 不同计算模型对比
4. 与实验结果的对比
5. 结论与讨论

岩石钻孔爆破粉碎区计算模型的改进

doi: 10.11883/1001-1455(2015)01-0101-07

作者简介: 冷振东(1989—), 男, 博士研究生

通讯作者: 卢文波, wblu@whu.edu.cn

计量

出版历程