基于块系构造的大规模地下爆炸不可逆位移计算方法*

李 杰，陈 伟,2，施存程，王明洋

(1.陆军工程大学爆炸冲击防灾减灾国家重点实验室,江苏 南京 210007；2.军委后勤保障部工程质量监督中心,北京 100850)

大量地下爆炸实验资料表明，基于连续介质力学概念的分析模型不能准确描述实际岩体的变形破坏过程[1-4]：(1)实际的大当量地下爆炸实验数据表明，在岩体的破裂区之外，距离爆心很远的距离上，发生了大的断层面的不可逆位移[5-9]，不可逆位移的范围要远远超出按连续介质模型的计算结果(要超出一个数量级)；(2)大当量地下爆炸时，岩体变形呈现强烈的变形局部化现象，在地下核爆炸实验中监测到的岩体变形数据，因测点位于岩块内或岩块边缘而导致其应变值相差2～3个量级[2-4]；(3)在岩体中进行地下核爆炸时，岩体裂隙、节理两侧的位移相差很大，甚至符号(方向)相反[10-13]，表明在岩体中由断裂、裂隙、节理所分割的岩块发生了旋转和相互间的滑移，这些现象和规律利用传统的连续介质力学模型无法得到，是目前岩体动力学研究领域所面临的非常重要而实际的研究课题。

钱七虎等[2,12]曾指出，地下爆炸时连续介质力学模型不适用性产生的物理本质主要在于没有考虑岩体作为地质体的非连续构造的块体性，以及没有考虑由岩块块体性所导致的附加自由度。均匀连续介质模型相对于块体集合来说是增加了“约束”，其计算结果较之实测数据必然低估了地运动参数。

因此，对于大规模地下爆炸不可逆位移的估算，必须考虑岩体的块体构造特性。实际上，岩体不是连续介质，在长期的地质构造运动作用下，大量的被低强度介质所充填的裂隙将完整的岩体切割成尺度各异的岩块。在地下爆炸作用远区，冲击波的作用不足以破坏强度较高的岩块，局部不可逆现象的产生是块体受限转动以及裂隙间填充物被破坏的结果，岩体的变形主要特征是块体沿结构破裂面的相对位移。正是基于上述思想，Kocharyan等[4]从力学角度对爆炸扰动条件下构成岩体的岩块受限运动进行受力分析，得到了地下大规模爆炸时激活块体的尺度和不可逆位移范围。本文中在文献[4]的基础上，基于岩体的块系构造理论对Kocharyan计算模型进行完善，补充其在激活块体尺度计算式中忽略的条件，重新推导并改进了深部大规模爆炸的不可逆范围表达式。

1 岩体的块系构造理论

众所周知，岩体是由不同层次的构造单元所组成的集合体，单元之间存在裂隙并由比岩石强度弱的材料填充，按照Sadovsky[14]块系构造等级学说，以及Kurlenya等[15]、戚承志等[16]的观点，从大到大陆板块或区域性岩体，小到细砂的结构组成甚或晶体的结构单元，岩体的结构具有嵌入特性或者层次重复性，在大的部分中嵌入小的部分，而后者又有更小的部分嵌入，如图1所示。根据岩体中已有的构造破坏对介质构造单元按大小进行排序，每一等级的岩块都由其对应等级的断层或者裂隙所构成的网所划分，最大的等级i=0对应于我们所需要研究的爆炸扰动范围，1级岩块由0级岩块范围内最大的断层或者裂缝来区分，逐次类推，这样弱扰动范围内的每一构造单元都可以被赋予等级i=1,2,3，…，并且相同等级的岩块具有相同的特征尺寸Li。第i级构造单元的变形特性由其组成部分(块体及其间隙)的变形及强度参数决定[10-11]，当由于爆炸所产生的应力幅值低于块体之间填充物的强度时，整个岩体结构并不会发生破坏，此时可将整个岩体视为均质连续体并且其运动及变形行为可由经典的弹性理论描述。当爆炸应力提高但又不超过岩块强度时，构成岩体的岩块将在爆炸和高地应力的双重作用下发生受限平动或转动，从而导致岩块间填充物发生形变，在此过程中块体的形状和体积可能并不发生明显的形变，完整岩体分解为大量块体的原因可认为是一个或者若干个软弱面在某一空间区域达到填充体的极限应力状态，致使岩体发生分裂进而形成块体单元，此时所形成的构造单元的尺度应当由施加于介质的作用的尺度及强度所决定，并随时间和空间分布而发生变化。

假设在爆炸扰动波的传播过程中，尺度为L*的块体能够被激活，L*与自然岩块系的块体尺度Li相对应，关于Li尺度的确定，需要提及的是Sadovsky块系构造等级学说的一个重要发现，即构成岩体的岩块的尺寸存在着成簇效应。Sadovsky等[17]通过对于泥炭的破碎、石英的侵蚀、土壤颗粒的分析，以及地下爆炸中岩石的破碎分析、露天采石场的岩石破碎分析、在建设水电站时利用地震地质方法对于岩石非均匀性的分析、利用地质方法对于地壳非均匀性的分析、地壳板块的分析、天体尺寸的分析等，发现了块体数随块体尺寸的分布中存在着优势尺寸，而且相邻的2个优势尺寸之比λ=Li/Li+1相对稳定，主要的变化范围为λ=2.0～5.5，λ被称为岩石的嵌入系数[16-17]。

同时实验还发现，裂纹的张开尺寸ai与被其所切割的相应等级的岩块尺度Li之间存在着统计的岩石力学不变量[16-17]μΔ=ai/Li=(0.5～2)×10-2，对于理解岩体在爆炸作用下的不可逆变形来说，岩石力学不变量μΔ可能具有重要的意义，μΔ的发现不仅证明了广泛应用于岩体应力应变理论计算中圣文南变形协调原理的不再适用[15]，同时利用μΔ可以给出实际岩体中岩块可能“自由度”的量值。

2 地下大规模爆炸诱发岩块不可逆位移的推导

在上述岩体的块系构造理论基础上，采用如图1所示的理想化计算模型，则可从力学的角度建立块体的运动学方程，从而针对具体的爆炸波分析构造介质的变形。

假设在给定爆炸当量Q作用下，在距离爆炸中心r处，尺度为L*的块体能够被激活，激活块体的尺度L*与岩体中的自然岩块尺寸Li对应，如图2所示。在计算过程中首先做如下假设：

(1)由于岩块的刚度远大于块系间填充物的刚度，因此在岩块的受限运动过程中，假设岩块为刚体；

(2)对于远场弱扰动环境，应力波的作用不足以破坏强度较高的岩块，局部不可逆现象的产生主要是块体受限运动导致块体间裂隙充填物被破坏的结果；

(3)由于爆炸扰动所激活的块体尺寸是均匀的，并且是标准的刚性立方体，所形成的块体单元的运动是独立的；

(4)相邻块体尺度之间满足关系：

Li/Li+1=k2k2>1

(1)

(5)块体之间的间隙与块体尺度之间满足关系：

ai/Li=k1

(2)

式中：k1、k2为岩体不变量值，按照Sadovsky[14]块系构造理论：

(3)

假设爆炸前，岩块间的初始相对位移、速度为零。在岩体中发生地下爆炸时，实际测试表明，在岩体中由断裂、裂隙、节理所分割的岩块发生了旋转和相互间的滑移(见图3)，从而使岩块与岩块间的软弱夹层产生剪切应变：

ε=Ui/ai

(4)

由于假设岩块是刚性的，即岩块不可压缩，因此第i级岩块的转动(或者平动)极限Ui*取决于刚性岩块接触面之间的间隙宽度ai及其充填层的极限应变εmax，即：

Ui*=ζai

(5)

式中：ζ为充填层物质的极限压缩比

ζ=εmax<1

(6)

假设岩体的构造符合如图1、4所示的构造层次及分布形式，在大规模爆炸扰动的作用下L*=Li层次以上的岩块(Li+1，Li+2，…)均发生转动或者平动。

首先,对于平动模型(见图4(a))，Li,Li+1,Li+2,…之间的层次结构决定了这n者有一个公共边界xjxj+1，那么该边界的宏观位移将是以上各尺度块体的位移之和：

U=Ui+Ui+1+Ui+2+…

(7)

假设各尺度岩块均达到运动极限状态，可得到受限平动条件下激活岩块的宏观极限位移：

(8)

对于受限转动模型，由于块体的位移Ui与转角ωi之间存在下列关系:

|Ui|=(Li/2)|ωi|

(9)

因此极限转角：

Δωi*=Ui*/(Li/2)=2ζai/Li

(10)

同样，Li,Li+1,Li+2,…之间的层次结构决定了这n者有一个公共端点xj，研究xj端点的位移，那么该端点的宏观位移将是以上各尺度块体的宏观位移之和：

U=Ui+Ui+1+Ui+2+…

(11)

最不利的情况为所有层次的岩块均达到了转动极限，则根据公式(10)～(11)可得到受限平动条件下激活岩块的宏观位移为：

(12)

式中:Δωi*(xj),Δωi+1*(xj),…分别表示端点xj(图4(b)中两种尺度块体的公共端点)因为尺度分别为Li,Li+1,…的块体的转动而引起的转角位移。

由式(8)、(12)可知，对于受限平动和转动，所引起的宏观极限位移表达式相同。

考虑到岩体的块系构造等级关系，即将公式(1)～(2)代入公式(12)可以得到：

(13)

因此从公式(13)可以看出，尺寸为L*的岩块的运动对岩体的变形具有极大的影响，而更高层次等级岩块(较小岩块)的运动则具有较小数量级的影响，如果忽略更高层次等级岩块运动的影响，公式(13)可以写成：

U*≈ζ(k1L*+k1(L*/k2))

(14)

对上式进行计算可得到激活块体的尺度计算式为：

L*≈k2U*/ζ(1+k2)k1

(15)

现根据上面所述的概念，针对具体的爆炸波形式评价激活块体的尺度和构造介质的变形，对于给定爆炸当量的地下大规模爆炸，实验测得的爆炸扰动波的形式为：

(16)

通过速度对时间积分得到爆炸扰动条件下岩体运动的位移表达式：

(17)

当t=τ时,上述位移达到最大值：

Umax=2v0τ/π

(18)

式中：Umax为通过实验测得的距离爆心r处的岩体最大位移，而U*为尺度为L*的最大激活岩块发生受限运动所能够获得的最大位移，如果令U*=Umax，并将公式(18)代入(15)，得到激活块体的尺度与爆炸当量及爆心距之间的关系：

(19)

上述公式还可以写成：

(20)

从公式(19)～(20)中可以看出，激活块体的尺度由爆炸当量和距离爆心的比例距离决定，爆炸当量越大，则激活的块体尺度越大。

对公式(19)进行变换得到地下爆炸不可逆范围的计算公式：

(21)

根据公式(14)，ε*与ζ存在下列关系：

ε*=(U*/L*)≈k1ζ(1+(1/k2))

(22)

对于具体的软弱充填物质，ζ为定值，同时考虑到k1、k2的取值，则从公式(22)可以看出ε*的取值主要取决于k1，即岩石力学不变量μΔ对岩体的局部不可逆变形做出了积极贡献。

将式(22)代入式(21)可得到，以ε*评价的不可逆位移范围计算公式：

(23)

3 讨 论

利用式(23)评估地下爆炸不可逆位移范围，需要确定方程中的4个未知参数A、B、n、m。其中A、n可以利用地下爆炸实验的监测数据得到，Radionov等[18]统计了涵盖采矿爆破和地下核爆炸的大量地下爆炸实验，给出了A和n的统计值，如表1所示。

表1 不同岩石的常数值[18]Table 1 Parameters of different rocks

B和m的直接监测成果较少，可通过其他实验数据间接获得。对地下大当量爆炸的实验处理结果表明[18]：当r=Rcavity=(7～12)Q1/3时，τ≈10-3Q1/3；r=Rcrush=(30～40)Q1/3时，τ≈20×10-3Q1/3,这里Rcavity和Rcrush分别为爆炸的空腔半径和破碎区半径。通过上述数据得到：Rcrush≈(3～5)Rcavity，从而有：m≈1.86～2.72。如果令m≈2.3，Rcrack≈7Q1/3，可以得到：B≈1.14×10-5。

以花岗岩中爆炸为例进行计算，令纵波速度cP≈5 km/s，则对于给定的爆炸当量(计算取单位爆炸当量Q≈1 kt)，则利用公式可以得到不可逆位移的范围，计算结果列于表2中,表中Rd为计算所得的不可逆位移范围，v*为相应的位移速度。

表2 利用公式(23)计算得到的不可逆位移范围Table 2 Irreversible displacement region radius calculated by Eq.(23)

通过计算可以看出，计算结果与俄罗斯学者的实测实验数据[2,13,18]基本吻合，即：Rd/Q1/3≈(650～1 400) m/kt1/3,v*=0.05～0.15 m/s。这一计算结果远远超出按连续介质力学计算所得到的结果[12,19-21]：Rd/Q1/3≈(80～120) m/kt1/3。

4 结 论

主要基于岩体块系构造理论，对Kocharyan等[4]的工作进行完善和改进，推得了大规模地下爆炸激活块体尺寸和诱发远区局部不可逆位移计算公式，主要结论如下：(1)本文工作与Kocharyan等[4]的推导过程的主要区别在于，Kocharyan的计算公式只考虑了岩块在爆炸扰动作用下的受限转动，而本文中则考虑了平动和转动；此外，本文中引入了充填层物质的极限压缩应变ζ，重新定义了岩体的等效特征应变ε*。(2)由ε*的定义公式可以看出，岩石力学不变量μΔ(或者k1)对于岩体的局部不可逆变形做出了积极贡献，它给出了实际岩体中岩块运动的可能“自由度”的量值。(3)由于岩体中节理裂隙(或者说多余自由度)的存在，使得实际岩体中不可逆位移的范围远远大于按连续介质力学所得的岩体不可逆位移范围。(4)本文推导公式的计算结果与实测结果基本一致。