Influence mechanism of rock joints on blasting fragmentation
-
摘要: 结构面对爆破块度形成的影响一直是爆破工程中的热点问题，与工程的长期运行安全和经济效益关系密切。首先基于玉龙喀什水利枢纽开展级配料爆破原位试验，通过对比岩体天然块度与爆破块度的对应关系研究了结构面对爆破块度形成的影响特性，结果表明对于5 mm以下的小粒径块度，爆破荷载与岩石本身的动抗压强度起主要作用；对于中大粒径块度，应力波的扰动和结构面对原始岩体的切割是最重要的影响因素，试验结果表明结构面对级配形成影响具有明显的分区特性；进而采用LS-DYNA的二次开发技术开展爆破块度数值仿真，结果表明由于粉碎区的冲击波峰值大，结构面对小粒径块度的形成影响并不明显，但结构面的存在改变了应力波传播过程，对于中大粒径块度则有明显影响；数值仿真与原位试验相互印证。最后在传统的KUZ-RAM模型基础上，引入修正的小粒径级配预测方法，采用自学习的方式修正岩体参数项，提出了基于分区机制的岩体爆破块度预测模型，结合工程数据与已有的爆破块度预测模型进行验证，结果表明本文所提的爆破级配分区预测模型在预测精度方面有所提高。Abstract: The influences of rock joints on the formation of blasting fragmentation have always been a hot issue in blasting engineering, which is one of the most important factors affecting the construction efficiency and long-term operation safety. Based on the explosion tests on the graded materials of Yulongkashi Water Conservancy Project, the influences of rock joints on the formation of blasting fragmentation are studied. The test results show that the fragmentation with rock particle size below 5 mm is mainly formed by the blasting loads and dynamic compressive strength but rarely by the rock joints. However, the influences of the rock joints on the fragmentation rock particle size above 5 mm are significant, which is reflected in the partition characteristics for the formation of rock blasting fragmentation. Further, the numerical simulation of rock blasting fragmentation is carried out by using the LS-DYNA secondary development technology, and the mechanical mechanism of the rock joints on rock blasting fragmentation is analyzed. The results show that the rock joints greatly affect the stress wave propagation in the fracture zone but produce a small effect on the formation of the fragmentation wiht small rock particle size, which is verified by the test results. Based on the dual influences of rock blasting fragmentation and rock joints, a partition prediction model for the rock blasting fragmentation is proposed. The rock blast fragmentation for many test samples is predicted by the proposed model and the existing models. By comparing the predicted results, it is found that the proposed model has a higher accuracy in predicting rock blasting fragmentation.
-
0. 引言
降雨入渗是造成滑坡最为重要的诱因之一[1]，降雨入渗会降低土体抗剪强度，增加土体重度，进而导致边坡下滑力增大，诱发边坡失稳。发展一种有效的降雨入渗模型进而在最大程度上降低降雨对边坡的危害，是降雨作用下边坡稳定性分析的关键。目前常用的边坡降雨入渗分析方法主要有[2]：①基于Richards方程的数值分析方法；②基于物理入渗模型的解析计算方法，如Green-Ampt（GA）入渗模型、Smith模型和Philip模型等。数值分析方法需要不断迭代求解降雨入渗控制方程，计算过程较为繁琐，对于复杂问题将十分耗时，并且若是网格尺寸和时间步长选择不当，将会导致计算不收敛。相较之下，以GA入渗模型为代表的解析计算方法，模型参数具有明确的物理意义，计算过程便捷，在降雨入渗下边坡稳定性研究中得到了广泛应用[3-14]。如Chen等[3]基于GA入渗模型研究了边坡降雨入渗过程，推导了边坡降雨入渗率计算表达式。张杰等[4]考虑气阻作用建立了分层假定的雨水入渗计算分析的GA模型。Zhang等[5]基于GA模型建立了可以计算时变破坏概率和失效时间的可靠度分析框架，并开发了基于电子表格的降雨型边坡破坏概率预测方法。Dou等[6]考虑土体饱和渗透系数变异性，基于GA入渗模型建立了雨水入渗-重分布分析模型。潘永亮等[7]建立了可综合考虑地下水位、初始含水率及不同降雨工况的改进GA模型。蒋水华等[8]提出了一种可考虑土体参数空间变异性的修正GA模型，并探讨了降雨入渗与多参数空间变异性作用下边坡失稳机理。
虽然上述研究大大地拓展了GA入渗模型的应用范围，但是忽略了土壤中客观存在的湿润区非饱和过渡层。大量试验发现湿润锋至土体入渗表面之间的土体并非处于完全饱和状态，在饱和区和天然区之间至少存在一个过渡层，并且过渡层内不同深度处土体体积含水量也不尽相同[9-10]。为此，张杰等[11]基于土体分层假定，推导了GA模型降雨历时与湿润锋深度之间的函数关系。胡海军等[12]在分层GA模型基础上，合理考虑湿润锋处吸力作用和水分剖面形状，分析了积水和降雨作用下非饱和重塑黄土的降雨入渗过程。Yao等[13]将饱和层与非饱和层厚度设定为入渗区深度的50%，分别对无限长边坡和马家沟滑坡进行了稳定性分析，并发现分层GA模型对研究强降雨事件下边坡稳定性问题十分有效。温馨等[14]针对黄土入渗中存在的非饱和区，提出了一种修正的GA模型来表征土体实际渗透系数随埋深的变化规律。
另一方面，受沉积、后沉积等地质作用的影响，岩土体参数存在固有的空间变异性。现有研究[8, 15-16]表明土体水力参数的空间变异性是影响降雨入渗下边坡渗流计算的关键因素。为客观评价降雨入渗下边坡失稳机理及稳定性，需要合理表征岩土体参数空间变异性的影响。然而，目前缺乏适用于考虑土体渗透系数空间变异性多层非均质边坡的降雨入渗模型。为此，提出了一种改进的Green-Ampt模型，通过4种常见的边坡工况验证该模型的有效性，进而探讨降雨入渗下多层非均质边坡含水率分布及稳定性问题。
1. 改进Green-Ampt入渗模型
1.1 Green-Ampt入渗模型
Green-Ampt（GA）入渗模型是1911年由Green等[17]提出的，最初用于解决初始干燥土体在薄层积水时的入渗问题。该模型假定雨水在入渗时存在可将湿润区和非湿润区明确分开的水平湿润锋面（图 1）。GA模型假设土体体积含水量呈矩形分布，湿润区为土体饱和层，其体积含水量为饱和体积含水量θs，非湿润区为天然层，其体积含水量为初始体积含水量θi。
![]()
图 1 经典GA入渗模型的含水率分布Figure 1. Distribution of water content of classic Green-Ampt model据此，土体入渗率f可表示为
f=−ksdhdz=−ks(hi−zf)−h0zf−0=ks(1+h0−hizf)=ks(1+h0+Sfzf)。 (1) 式中：ks为土体饱和渗透系数；h0为地面处积水水头；zf为湿润锋深度；hi为湿润锋下方的水头；Sf为湿润锋处土体概化基质吸力水头，是一个关于含水量的函数，可通过相关试验或下面的推导方程进行求解[18]：
Sf=∫1Kr(θi)ψ(Kr)dKr1−Kr(θi)=3λ+23λ+1ψb[1−Kr(θi)3λ+13λ+2]1−Kr(θi)≈3λ+23λ+1ψb。 (2) 式中：ψb为土体进气值；Kr为湿润锋处相对饱和渗透系数；λ为表征土体孔径分布特征的拟合系数。进而，可获得累积降雨入渗量I为
I = (θs−θi)zf。 (3) 可见，该模型计算简单，且模型参数具有明确的物理意义。然而，针对边坡降雨入渗问题，该模型存在如下不足：①在模拟降雨入渗时，该模型没有考虑到降雨强度R对入渗率f的影响，将会高估入渗量I，导致计算的湿润锋深度偏大。②该模型没有考虑客观存在的过渡层，大量试验研究[19-21]表明，湿润锋至土体入渗表面之间的土体并非完全处于饱和状态，即在饱和区和天然区之间应至少存在一个过渡层。如不考虑客观存在的过渡层，将会高估湿润锋以上的土重，低估边坡安全系数。③当只考虑参数（如饱和渗透系数）变异性时，传统的GA模型可以采用临界渗透系数划分降雨入渗边界条件，但是一旦考虑土体参数空间变异性，整个土层剖面的渗透系数呈空间非均匀分布，不同位置处数值大小不同，此时采用临界渗透系数划分降雨入渗边界条件不再适用[8]。
1.2 改进Green-Ampt入渗模型
为克服上述不足，十分必要对经典GA入渗模型进行改进，使其能够处理考虑土体参数空间变异性的多层非均质边坡降雨入渗问题。实际工程边坡或天然边坡通常存在一定的坡度，因而积水易形成坡面径流，其中的关键难题是如何确定不同时刻下边坡降雨入渗边界条件。鉴于采用临界渗透系数划分降雨入渗边界条件不适用考虑参数空间变异性的多层非均质边坡，为此蒋水华等[8]提出采用开始积水时间tp划分入渗边界，本文将该GA模型称为“修正GA模型”。由GA入渗模型[11, 13]可知，土体入渗率f随着累计入渗量I的增加而逐渐减小。对于倾角为α的边坡而言，当土体的累计入渗量达到某一累积入渗量Ip值时，边坡坡面开始积水[5, 7]，此时f = Rcosα，Ip的计算表达式为
Ip=(θs−θi)Sfcosα(R/ks−1)。 (4) 式中：Ip为积水时刻tp所对应的累积入渗量；α为边坡倾角；R为降雨强度。据此可得边坡开始积水的时间tp为
tp=Ip/(Rcosα)。 (5) 根据开始积水时间tp可划分入渗边界条件如下：
（1）当降雨历时t < tp时，入渗区土体处于非饱和状态，雨水将全部入渗至土体中，此时的入渗边界由流量控制，可得到入渗控制方程为
f=Rcosα=k(θ)+ks[ψr(θ)−ψr(θi)]zf。 (6) 式中：θ为土体体积含水量；ψr(θ)为湿润锋处土体相对吸力水头[8]，
ψr(θ)=∫−∞ψk(θ)ksdψ=∫−∞ψ(ψbψ)3λ + 2dψ=(ψbψ)3λ + 1ψb3λ + 1=Se3+1/λψb3λ + 1。 (7) 式中：θr为土体残余体积含水量；Se为有效饱和度，Se=(θ−θr)/(θs−θr)。
随后，依据GA模型的矩形入渗原理可推导相应的湿润锋深度为
zf=Iθ−θi = Rtcosαθ−θi。 (8) 将式（8）代入式（6）中，可得到非线性方程：
f=Rcosα=k(θ)+ks[ψr(θ)−ψr(θi)](θ−θi)Rtcosα。 (9) 求解式（9）所示的非线性方程，可得到不同降雨历时下土体体积含水量θ分布，进而结合式（8）可以计算得到对应的湿润锋深度。
（2）当降雨历时t ≥ tp时，此时边坡表面开始积水，因降雨强度超过土体允许入渗率，故土体表面入渗区处于饱和状态，此时的入渗边界由水头控制[2, 6, 8]。降雨入渗过程可分为自由入渗和积水入渗两个阶段。其入渗率f和累计入渗量I可为[11, 13]
f(t)={Rcosα(t⩽tp)ks(cosα+Sf/zf)(t>tp), (10) I=∫t0f(t)dt。 (11) 按照式（11）进行积分得到累积入渗量，进而得到入渗饱和区湿润锋深度zf为
zf=Iθs−θi。 (12) 因此，一旦确定了土体饱和渗透系数ks和相应的降雨历时t，便可针对上述两种降雨入渗边界，分别求解对应的土体湿润锋深度和含水量分布。
另一方面，大量相关试验研究[9-10, 19]表明：湿润锋至土体入渗表面之间的土体真实体积含水量并非完全处于饱和状态，即在饱和层和天然层之间应至少存在一个过渡层。根据文献[9~13]，将GA入渗模型进行分层处理，增加客观存在的过渡层，如图 2（b）所示。从图 2中可以看出，改进后的土体体积含水量分布包含饱和层、过渡层和天然层3部分。其中饱和层深度为zs，过渡层深度为zt。
![]()
图 2 改进的GA入渗模型的含水量分布Figure 2. Distribution of water content of improved Green-Ampt model文献[9，10]的相关试验数据表明，采用椭圆曲线分布来表征过渡区含水量分布具有较高的计算精度，故本文也采用椭圆曲线来描述过渡层含水量分布，由此可得修正的土体含水量分布为
θ(z)={θs (饱和层, 0≤z≤zs)θi+(θs−θi)√1−((z−zs)2zt2) (过渡层, zs≤z≤zh)θi (天然层, z>zh)。 (13) 式中：zh为整个入渗区深度，
zh=zs+zt。 (14) 随着降雨入渗的推进，饱和层、过渡层和天然层的深度都会发生变化。彭振阳等[10]通过试验发现过渡层深度与总入渗区深度zh存在一个比例关系，并且该比例与整个入渗区深度之间具有良好的线性关系：
zt=ηzh, η=azh+b。} (15) 式中：η为过渡层占整个入渗区深度的比例；a和b为拟合系数。
在过渡层修正过程中，应保持累计入渗量不变，即I=IM。I可通过式（8）和（11）计算得到，IM为增加过渡层后的累计入渗量，其值应等于饱和层入渗量Is和过渡层入渗量It之和。其中对于饱和层，入渗量Is =(θs−θi)zs；对于过渡层，由于其为椭圆曲线分布，故其入渗量等于1/4的椭圆面积It = 0.25π(θs−θi)zt，据此可得总的累计入渗量为
I=IM=Is+It=(θs−θi)zs+0.25π(θs−θi)zt。 (16) 联立式（14）~（16），便可计算得到饱和层深度zs和过渡层深度zt，进而可求得整个边坡土层的入渗区深度zh及体积含水量分布。本文将所提模型称为“改进GA模型”
2. 多层非均质边坡稳定性分析
降雨诱发滑坡通常是沿边坡浅层面发生，在剪切力作用下趋于破坏，滑动面深度远小于滑动面长度，故常采用图 3所示的无限长边坡模型分析浅层滑坡稳定性。根据莫尔-库仑破坏准则，基于极限平衡分析的无限长边坡模型安全系数计算公式为[15, 22]
FS = τfτm = c′+[(σn−ua)−σs]tanφ′Wsinαcosα。 (17) 式中：τf为土的抗剪强度；τm为沿潜在破坏面任意一点的剪应力；W为单位土体切片的重量；c′，φ′分别为土体有效黏聚力和有效内摩擦角；(σn−ua)为单位土条底部净法向力，ua为孔隙气压力。根据文献[2，6，8，15]，取ua = 0，σn=γtzcos2α，其中γt为土体有效重度，z为滑动面深度，α为无限长边坡倾角；σs为吸应力，
σs=−Seψ(θ)。 (18) 式中：ψ(θ)为基质吸力，是关于土体体积含水量θ的函数。
为了考虑降雨入渗过程中因体积含水量变化而引起的土体单位重量变化，对滑动面以上土体部分进行积分得到土体重量为
W=∫z0γtdz。 (19) 式中：土体有效重度γt[15]为
γt=γd+θ(z)γw。 (20) 式中：γd和γw分别为土体干重度和水重度。
综上，基于本文提出的改进GA模型进行降雨入渗下非均质边坡稳定性分析的步骤如下：
（1）收集并确定土体参数统计特征（包括均值、标准差、自相关函数、波动范围及概率分布等），沿垂直方向将边坡划分为H/d个等厚的均质土层（H为坡高，d为随机场单元尺寸）。
（2）基于土体参数统计特征，采用蒙特卡洛模拟（MCS）方法产生N组独立标准正态随机样本，随后基于Karhunen-Loève（KL）级数展开方法[15-16, 23]生成N组土体参数随机场实现值，并将其一一赋值给对应的土层。
（3）采用式（5）计算开始积水时间tp，并与降雨历时t进行比较，划分降雨入渗边界条件，再根据式（8），（12），（13）~（16），求解考虑客观存在过渡层的降雨入渗下边坡体积含水量分布。
（4）采用式（17）计算边坡每个土层及湿润锋处安全系数，从中可获得（H/d+1）个安全系数，并将其中最小的安全系数视作该边坡安全系数，对应的滑动面则为该边坡最危险滑动面。
3. 算例分析
下面以图 3所示的非饱和无限长边坡为例[6, 24]，验证提出的改进GA模型的有效性。边坡坡高H = 3 m，倾角α = 50°，下部为不透水基岩，降雨强度R恒定为5 mm/h。采用Brook等[25]提出的土水特征曲线模型描述土体渗透系数和基质吸力与体积含水量之间的函数关系为
k(θ)=ksSe3+2/λ, (21) ψ(θ)=ψbSe−1/λ。 (22) 式中：k(θ)为土体渗透系数。根据文献[6，24]，得到土体物理力学参数取值如表 1所示。另外根据文献[10]的试验结果，取式（15）中的系数a和b分别为-0.003，0.8712。
表 1 土体物理力学参数取值Table 1. Values of physical parameters of soil计算参数 取值 计算参数 取值 饱和渗透系数ks/(mm·h-1) 3 有效黏聚力c′/kPa 5 有效内摩擦角φ′/(°) 28 进气值ψb/kPa 2.752 土体干重度γ /(kN·m-3) 16.217 湿润锋处土体概化基质吸力Sf /mm 424.3 饱和含水量θs 0.335 残余含水量θr 0.068 初始含水量θi 0.148 土体孔隙分布特征参数λ 0.319 通过以下4种常见的边坡工况，验证采用本文提出的改进GA模型进行降雨入渗下边坡体积含水量计算及稳定性分析的有效性：①均质边坡；②两层边坡（上层渗透性大、下层渗透性小）；③两层边坡（上层渗透性小、下层渗透性大）；④考虑任意一次参数随机场典型实现的非均质边坡。同时，选用Richards方程数值解和修正GA入渗模型[8]计算结果进行对比分析。其中，Richards方程数值解采用一维流动模型来模拟边坡垂直入渗过程，描述一维水流在边坡上的入渗过程控制方程为[26]
∂θ∂t=∂∂z[k(θ)(∂h∂z−cosα)]。 (23) 式中：θ∈[θr,θs]；h为压力水头。
本文借助Hydrus-1D软件[27]对式（23）所示Richards方程进行数值求解，可得到不同降雨历时t下不同深度z处压力水头h(z, t)和土体体积含水量θ(z, t)分布。本文将Richards方程数值解视作精确解进行判断对比分析。
3.1 工况1：均质边坡
基于表 1中的土体参数均值，对图 3所示均质边坡模型进行降雨入渗分析。图 4比较了采用修正GA模型和提出的改进GA模型计算的体积含水量分布和安全系数及Richards方程数值解。
![]()
图 4 不同降雨历时下边坡体积含水量和安全系数分布的比较（工况1）Figure 4. Comparison of distribution of water content and factors of safety of slope under different rainfall durations (Scenario Ⅰ)由图 4可知，采用修正GA模型计算的体积含水量分布与提出的改进GA模型计算结果和Richards方程数值解均存在较大差别，这是由于修正GA模型假定雨水呈矩形入渗，没有考虑客观存在的过渡层。在此基础上进行边坡稳定性分析，其边坡FS沿坡体深度的变化如图 4（b）所示。由图 4可知，随着降雨的进行，雨水入渗区也不断下移，饱和层内土体基质吸力基本丧失，导致入渗区内边坡FS随着埋深的增大而逐渐减小。到达过渡层后，由于土体并非完全处于饱和状态，基质吸力的作用未完全丧失，边坡FS先下降而随后又呈增加趋势。接着，虽然土体天然层内土体含水量和基质吸力均保持初始状态，但由于土体深度不断增加，导致土条自重逐渐增加，坡体滑动力增大，此时的边坡FS表现为随着埋深的增大而逐渐降低，并在埋深z = 3.0 m的不透水层处达到最小值。
表 2比较了不同方法计算的降雨历时20，36，60 h下边坡安全系数和最危险滑动面深度。表中，入渗区是指雨水入渗至坡体内的区域，包括饱和层和过渡层（例如降雨历时60 h，图 4（a）中入渗区可视为距离坡面0~1.25 m区域）；整个边坡则包含入渗区和天然区。将其分开进行处理，可分析不同降雨历时下边坡最危险滑动面位置和确定边坡失稳是否因降雨所诱发。若边坡最危险滑动面位于入渗区，可认为边坡失稳因降雨造成。由表 2可知，边坡安全系数随着降雨历时的增加而减小。相较于修正GA模型，提出的改进GA模型与基于Richards方程数值方法计算的安全系数吻合。如降雨历时60 h，提出的改进GA模型计算的边坡安全系数为1.22，与基于Richards方程数值方法计算的1.19差别很小。表 2中采用修正GA模型计算的边坡安全系数偏小的原因是，该模型假定雨水呈矩形入渗，导致计算的体积含水量θ偏大，进而高估了土条重量W，造成边坡安全系数偏小。相较之下，提出的改进GA模型因考虑了客观存在的过渡层，计算的边坡体积含水量、安全系数以及最危险滑动面位置均与Richards方程数值解吻合，这表明提出的改进GA模型能够有效分析均质边坡降雨入渗问题。此外，降雨历时20，36 h时，3种方法均识别出边坡最危险滑动面位于z = 3.0 m不透水层处。降雨历时60 h，相较于修正GA模型计算结果（1.03 m），本文提出的改进GA模型结果（0.99 m）与精确解（0.99 m）更吻合，说明提出方法能够准确识别边坡最危险滑面位置。
表 2 边坡安全系数和最危险滑动面深度的比较（工况1）Table 2. Comparison of factors of safety and depths of critical slip surface of slope (Scenario Ⅰ)类别 降雨历时/h Richards方程数值解 修正GA模型 提出的改进GA模型 入渗区安全系数 20 2.79 2.35 2.78 36 1.72 1.53 1.74 60 1.19 1.10 1.22 整个边坡安全系数 20 1.36 1.36 1.36 36 1.34 1.34 1.34 60 1.19 1.10 1.22 最危险滑动面深度 20 3.00 m 3.00 m 3.00 m 36 3.00 m 3.00 m 3.00 m 60 0.99 m 1.03 m 0.99 m 3.2 工况2：两层边坡（上层渗透性大、下层渗透性小）
为了进一步说明提出的改进GA模型的有效性，本节对两层边坡（上层渗透性大、下层渗透性小）进行降雨入渗分析。第一层土饱和渗透系数为3.5 mm/h，土层厚度为0.5 m；第二层土饱和渗透系数为3 mm/h，土层厚度为2.5 m。两个土层的其他参数均相同，见表 1。图 5比较了工况2下采用上述2种模型和Richards方程数值方法计算的土体体积含水量分布。由图 5可知，提出的改进GA模型计算结果同样与Richards方程数值解吻合。表 3比较了采用不同方法计算的降雨历时20，36，60 h安全系数及最危险滑动面深度。可见，提出的改进GA模型与Richards方程数值方法计算的边坡安全系数依然吻合。如降雨历时20 h，Richards方程数值方法计算的入渗区安全系数为2.83，提出的改进GA模型结果为2.81，而修正GA模型结果为2.39，进一步体现了提出的改进GA模型的有效性及考虑客观存在的过渡层的必要性。
![]()
图 5 不同降雨历时下边坡体积含水量分布的比较（工况2）Figure 5. Comparison of distribution of water content of slope under different rainfall durations (Scenario Ⅱ)表 3 边坡安全系数和最危险滑动面深度的比较（工况2）Table 3. Comparison of factors of safety and depths of critical slip surface of slope (Scenario Ⅱ)类别 降雨历时/h Richards方程数值解 修正GA模型 提出的改进GA模型 入渗区安全系数 20 2.83 2.39 2.81 36 1.71 1.55 1.76 60 1.19 1.10 1.22 整个边坡安全系数 20 1.36 1.36 1.36 36 1.34 1.34 1.34 60 1.19 1.10 1.22 最危险滑动面深度 20 3.00 m 3.00 m 3.00 m 36 3.00 m 3.00 m 3.00 m 60 0.98 m 1.03 m 0.99 m 3.3 工况3：两层边坡（上层渗透性小、下层渗透性大）
本节进一步分析两层边坡（上层渗透性小、下层渗透性大）来验证改进GA模型的有效性。第一层土饱和渗透系数为3 mm/h，土层厚度为0.5 m；第2层土饱和渗透系数为3.5 mm/h，土层厚度为2.5 m。这两个土层的其他参数均相同，见表 1。图 6比较了采用上述2种模型和Richards方程数值方法计算的边坡体积含水量分布。表 4比较了基于不同方法计算的边坡安全系数。与工况2类似，提出的改进GA模型与Richards方程数值方法计算的安全系数一致，说明提出的改进GA模型可有效求解多层边坡降雨入渗问题。
![]()
图 6 不同降雨历时下边坡体积含水量分布的比较（工况3）Figure 6. Comparison of distribution of water content of slope under different rainfall durations (Scenario Ⅲ)表 4 边坡安全系数和最危险滑动面深度的比较（工况3）Table 4. Comparison of factors of safety and depths of critical slip surface of slope (Scenario Ⅲ)类别 降雨历时/h Richards方程数值解 修正GA模型 提出的改进GA模型 入渗区安全系数 20 2.79 2.35 2.78 36 1.72 1.53 1.74 60 1.19 1.10 1.22 整个边坡安全系数 20 1.36 1.36 1.36 36 1.34 1.34 1.34 60 1.19 1.10 1.22 最危险滑动面深度 20 3.00 m 3.00 m 3.00 m 36 3.00 m 3.00 m 3.00 m 60 1.01 m 1.04 m 1.00 m 3.4 工况4：考虑一次参数随机场典型实现的非均质边坡
工程实际中，土体饱和渗透系数ks呈空间非均质分布。本文将边坡沿垂直方向上划分为60个均质土层，厚度为d = 0.05 m[6]。基于KL级数展开方法来离散ks随机场，利用高斯型自相关函数模拟土体ks的空间自相关性，垂直自相关距离lv取0.5 m。当KL级数展开截断项数取6时，计算的期望能比率因子ε为95.67%，可以满足随机场离散精度要求（大于95%）。图 7给出了土体饱和渗透系数随机场的一次典型实现，图中颜色较深部分表示ks值较大区域，颜色较浅部分表示ks值较小区域。
![]()
图 7 饱和渗透系数随机场一次典型实现分布图Figure 7. Typical realization of random filed of hydraulic conductivity图 8比较了不同方法计算的降雨历时8，36，60 h边坡体积含水量分布。由图 8可知，在3个不同降雨历时下，提出的改进GA模型计算的体积含水量分布和Richards方程数值解同样吻合，入渗区和整个边坡的安全系数和最危险滑动面深度几乎相同（见表 5），同样说明了基于提出的改进GA模型进行多层非均质边坡降雨入渗分析的有效性。由表 5可知，降雨历时8，36 h，由这2种模型和Richards方程数值方法计算的整个边坡安全系数皆小于入渗区安全系数，此时边坡最危险滑动面均位于z = 3 m的不透水层处。然而，在降雨历时60 h时，提出的改进GA模型和Richards方程数值方法均表明边坡最危险滑动面位于z = 3 m不透水层处。而基于修正GA模型计算的整个边坡安全系数等于入渗区安全系数，即修正GA模型认为此时最危险滑动面处于入渗区（距离坡面约0.85 m），显然错误评估了边坡失稳模式。这是由于该典型实现下土层表面渗透性低，湿润锋推进缓慢，虽然湿润锋处安全系数在不断减小，但是始终大于不透水层处安全系数，故而边坡仍然沿着不透水层处发生失稳破坏。综上，采用修正GA模型会造成对边坡失稳机理的错误认识，相比之下采用提出的改进GA模型能够准确识别最危险滑动面位置，获得更合理的边坡安全系数。另外需要说明的是，提出的改进Green-Ampt模型和Richards方程数值解得到的过渡层分布存在误差，并且该误差随着降雨历时的增加而增大，见图 4（a），5，6，这是因为本文没有考虑式（15）中系数a和b随着降雨历时的增加有可能发生变化。尽管如此，两者计算的湿润锋深度随着降雨历时的增加均保持一致。
![]()
图 8 工况4不同降雨历时下边坡体积含水量分布比较Figure 8. Comparison of distribution of water content of slope under different rainfall durations (Scenario Ⅳ)表 5 边坡安全系数和最危险滑动面深度的比较（工况4）Table 5. Comparison of factors of safety and depths of critical slip surface of slope (Scenario Ⅳ)类别 降雨历时/ h Richards方程数值解 修正GA模型 提出的改进GA模型 入渗区安全系数 8 6.58 5.62 6.43 36 1.95 1.72 1.98 60 1.40 1.24 1.39 整个边坡安全系数 8 1.37 1.37 1.37 36 1.35 1.35 1.35 60 1.33 1.24 1.33 最危险滑动面深度 20 3.00 m 3.00 m 3.00 m 36 3.00 m 3.00 m 3.00 m 60 h 3.00 m 0.85 m 3.00 m 为了说明这4种工况下边坡体积含水量分布之间的差别，图 9比较了这4种工况降雨历时36，60 h的边坡含水量分布曲线。由图 9（a）可知，工况1，2，3的边坡含水量分布曲线重合，而工况4入渗区较小。这是由于工况4边坡表面渗透性小，雨水入渗率较低，进而导致入渗区小。工况1，3由于上层饱和渗透系数相同，故其在降雨历时36 h下含水量分布曲线重合。工况2与工况1，3的边坡含水量分布曲线重合的原因是，在降雨初期，降雨历时未超过积水时刻tp，雨水入渗率等于Rcosα，此时的入渗量仅与降雨历时相关，与土体饱和渗透系数大小无关。因此，工况1，2，3在降雨历时36 h入渗量相差不大，对应的边坡含水量分布曲线基本重合。相比之下，降雨历时60 h，工况4由于其自身的渗透性较低，故入渗区小；工况3下由于下层土体渗透性大，导致入渗量增大，进而入渗深度更大；工况1，2由于下层土体的饱和渗透系数一致，故这两种工况不仅入渗率和入渗量相等，而且边坡体积含水量分布也相同，见图 9（b）。
![]()
图 9 不同工况下边坡体积含水量分布的比较Figure 9. Comparison of distribution of water content of slope for different scenarios4. 结论
考虑降雨入渗作用下边坡土体中客观存在的过渡层，提出了一种适用于多层非均质边坡降雨入渗分析的改进GA模型，为合理确定不同降雨历时下多层非均质边坡的入渗区深度和土体体积含水量分布奠定了基础。同时，以无限长边坡模型为例，进行降雨入渗下边坡稳定性分析，并与采用Richards方程数值解和修正GA模型计算的边坡安全系数和最危险滑动面深度进行了对比分析，验证了提出的改进GA模型的有效性。
（1）修正GA模型局限于雨水呈矩形入渗，没有考虑雨水入渗过程中土体客观存在的过渡层，导致降雨入渗过程中单位土体重度被高估，计算的边坡安全系数偏小，造成对边坡失稳机理的错误认识。
（2）提出的改进GA模型因考虑了边坡土体入渗区内客观存在的过渡层，故计算的土体含水量分布、安全系数和滑动面分布不仅贴近客观实际，而且与Richards方程数值解吻合，能够较好地解决考虑土体饱和渗透系数空间变异性的多层非均质边坡降雨入渗问题。
（3）相比于修正GA模型，提出的改进GA模型能够准确识别非均质边坡最危险滑动面位置，为降雨型滑坡加固设计及风险防控提供了一条有效途径。如工况4中，降雨历时60 h时，修正GA模型认为最危险滑动面位于入渗区，而改进GA模型和Richards方程计算结果均证明边坡最危险滑动面位于不透水层处。
（4）工程实际中存在多种环境因素对边坡降雨入渗过程具有重要的影响，包括边坡初始状态（如含水量，土壤类型等）的不均匀分布、边坡倾角大小、覆土厚度和降雨雨型等，融入多种环境因素影响的边坡降雨入渗模型还有待进一步研究。
-
![]()
图 4 原生结构面切割岩体形成天然块度示意图
Figure 4. Schematic diagram of natural fragmentation shaped by original joints cutting through rock mass
![]()
图 6 不同试验区的岩体爆破级配筛分曲线
Figure 6. Screening curves of blasting gradation of rock mass in different test regions
![]()
图 7 小粒径块度对比分析示意图
Figure 7. Illustration of comparative analysis of blasting fragmentation with small particle sizes
![]()
图 10 特定爆破参数花岗岩中应力波的衰减规律
Figure 10. Attenuation laws of stress wave in granite with specific blasting parameters
![]()
图 13 不同区域等效应力变化对比曲线
Figure 13. Comparison curves of variation of equivalent stress in different regions
![]()
图 16 有无结构面条件下的块度分布情况对比
Figure 16. Comparison of fragamentation distribution with and without joints
![]()
图 18 平均块度与单耗的拟合关系
Figure 18. Regression relationship between average fragmentation size and specific explosive consumption
表 1 现场爆破专项试验的部分参数
Table 1 Parameters of blasting fragmentation tests
试验序号 岩体体积节理数 孔距/m 排距/m 孔深/m 堵塞长度/m 钻孔角度/(°) 筛分方法 1 < 3 4 3 10 2.5 90 人工 2 3~10 4 3 10 2.5 90 人工 3 10~20 4 3 10 2.5 90 人工 4 3~10 4 3 10 2.5 90 人工 5 10~20 4 3 10 2.5 90 人工 
下载: 导出CSV
表 2 岩体结构面参数
Table 2 Parameters of rock joints
统计区段 分组 倾向/(°) 倾角/(°) 间距/m 迹长/m ① NEE/NW 332 75 1.24 1.57 NNE/NW 290 70 1.91 2.63 NNE/SE 116 62 1.79 2.98 NNW/SW 256 70 2.45 3.65 ② NWW/NE 19 67 1.91 2.93 NEE/NW 345 74 1.58 3.66 NNW/NE 71 63 0.93 4.18 NNE/NW 290 68 2.50 5.04
下载: 导出CSV
表 3 计算参数
Table 3 Parameters for calculation
密度/(kg·m-3) 动弹性模量/GPa 泊松比 动抗拉强度/MPa 损伤参数 KIC/(MN·m-3/2) 损伤系数λ/(kg·J-1) k m 2700 45 0.16 5 2.33×1024 7 1.13 0.0001
下载: 导出CSV
-
[1] 高启栋, 卢文波, 冷振东, 等. 岩石爆破中孔内起爆位置对爆炸能量传输的调控作用研究[J]. 岩土工程学报, 2020, 42(11): 2050-2058. doi: 10.11779/CJGE202011010 GAO Qidong, LU Wenbo, LENG Zhendong, et al. Regulating effect of detonator location in blast-holes on transmission of explosion energy in rock blasting[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(11): 2050-2058. (in Chinese) doi: 10.11779/CJGE202011010
[2] 吴新霞, 彭朝晖, 张正宇. 面板堆石坝级配料开采爆破块度预报模型及爆破设计参数优化研究[J]. 工程爆破, 1996, 2(4): 95-100. https://www.cnki.com.cn/Article/CJFDTOTAL-GCBP604.020.htm WU Xinxia, PENG Zhaohui, ZHANG Zhengyu. Study on optimization of blasting fragmentation forecasting model and blasting design parameters of the gradating material mining of the rockfill dam with concrete facing[J]. Engineering Blasting, 1996, 2(4): 95-100. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCBP604.020.htm
[3] 张有才, 朱传云. 堆石坝级配料爆破块度分布模型的研究[J]. 爆破, 2005, 22(1): 44-47. https://www.cnki.com.cn/Article/CJFDTOTAL-BOPO20050100A.htm ZHANG Youcai, ZHU Chuanyun. Blasting fragmentation forecasting model for concrete panel rock-fill dam[J]. Blasting, 2005, 22(1): 44-47. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BOPO20050100A.htm
[4] 蔡建德, 郑炳旭, 汪旭光, 等. 多种规格石料开采块度预测与爆破控制技术研究[J]. 岩石力学与工程学报, 2012, 31(7): 1462-1468. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201207020.htm CAI Jiande, ZHENG Bingxu, WANG Xuguang, et al. Research on blasting control technique and block size prediction of different dimensions stones[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(7): 1462-1468. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201207020.htm
[5] 吴新霞, 胡英国, 刘美山, 等. 水利水电工程爆破技术研究进展[J]. 长江科学院院报, 2021, 38(10): 112-120, 147. https://www.cnki.com.cn/Article/CJFDTOTAL-CJKB202110023.htm WU Xinxia, HU Yingguo, LIU Meishan, et al. Research progress of blasting technology in hydropower engineering[J]. Journal of Yangtze River Scientific Research Institute, 2021, 38(10): 112-120, 147. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-CJKB202110023.htm
[6] GRADY D E, KIPP M E. Dynamic rock fragmentation[J]. Fracture Mechanics of Rock, 1987, 10: 429-475.
[7] CHUNG S H, KATSABANIS P D. Fragmentation prediction using improved engineering formulae[J]. Fragblast, 2000, 4(3/4): 198-207.
[8] CUNNINGHAM C V B. The Kuz-Ram fragmentation model–20 years on[C]// Brighton Conference Proceedings. Brighton, UK: European Federation of Explosives Engineer, 2005.
[9] KUZNETSOV V M. The mean diameter of the fragments formed by blasting rock[J]. Soviet Mining Science, 1973, 9(2): 144-148. doi: 10.1007/BF02506177
[10] ROSIN P A U L. Laws governing the fineness of powdered coal[J]. Journal of Institute of Fuel, 1933, 7: 29-36.
[11] 武仁杰, 李海波, 于崇, 等. 基于统计分级判别的爆破块度预测模型[J]. 岩石力学与工程学报, 2018, 37(1): 141-147. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201801014.htm WU Renjie, LI Haibo, YU Chong, et al. Model for blasting fragmentation prediction based on statistical classification[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(1): 141-147. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201801014.htm
[12] 吴发名, 刘勇林, 李洪涛, 等. 基于原生节理统计和爆破裂纹模拟的堆石料块度分布预测[J]. 岩石力学与工程学报, 2017, 36(6): 1341-1352. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201706005.htm WU Faming, LIU Yonglin, LI Hongtao, et al. Fragmentation distribution prediction of rockfill materials based on statistical results of primary joints and simulation of blasting cracks[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(6): 1341-1352. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201706005.htm
[13] 王仁超, 连嘉欣, 邸阔. 结合深度学习和NCFS算法的堆石料粒度分布智能检测方法[J]. 水利学报, 2021, 52(9): 1103-1115. https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB202109010.htm WANG Renchao, LIAN Jiaxin, DI Kuo. Intelligent detection method of rockfill particle size distribution based on deep-learning and NCFS algorithm[J]. Journal of Hydraulic Engineering, 2021, 52(9): 1103-1115. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLXB202109010.htm
[14] 周继国, 张智宇, 王建国, 等. 最小抵抗线对岩石破碎块度影响的能量分析[J]. 地下空间与工程学报, 2021, 17(1): 135-142. https://www.cnki.com.cn/Article/CJFDTOTAL-BASE202101015.htm ZHOU Jiguo, ZHANG Zhiyu, WANG Jianguo, et al. Energy analysis on the effect of minimum burden on rock fragmentation[J]. Chinese Journal of Underground Space and Engineering, 2021, 17(1): 135-142. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BASE202101015.htm
[15] PAINE A S, PLEASE C P. An improved model of fracture propagation by gas during rock blasting: some analytical results[J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1994, 31(6): 699-706.
[16] SHOCKEY D A, CURRAN D R, SEAMEN L. Fragmentation of rock under dynamic loads[J]. Int J Rock Mech Sci & Geomech Abstr, 1974, 11: 303-317.
[17] ROSSMANITH H P, FOURNEY W L. Fracture initiation and stress wave diffraction at cracked interfaces in layered media: I brittle/brittle transition[J]. Rock Mechanics, 1982, 14(4): 209-233. doi: 10.1007/BF01596617
[18] 张继春. 节理岩体爆破的损伤机理及其块度模型[J]. 中国有色金属学报, 1999, 9(3): 21-27. https://www.cnki.com.cn/Article/CJFDTOTAL-ZYXZ199903043.htm ZHANG Jichun. Damage mechanism of blasting in jointed rock masses and its fragment size model[J]. The Chinses Journal of Nonferrous Metals, 1999, 9(3): 21-27. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZYXZ199903043.htm
[19] 郭文章, 王树仁. 节理岩体爆破块度预测的力学模型探讨[J]. 爆破, 1997, 14(3): 31-34. https://www.cnki.com.cn/Article/CJFDTOTAL-BOPO199703006.htm GUO Wenzhang, WNAG Shuren. Discussion on the mechanical model of jointed rock mass fragmentation forecast by blasting[J]. Blasting, 1997, 14(3): 31-34. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BOPO199703006.htm
[20] ALER J, DU MOUZA J, ARNOULD M. Evaluation of blast fragmentation efficiency and its prediction by multivariate analysis procedures[J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts, 1996, 33(2): 189-196.
[21] 石根华. 一般自由面上多面节理生成、节理块切割与关键块搜寻方法[J]. 岩石力学与工程学报, 2006, 25(11): 2161-2170. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200611002.htm SHI Genhua. Producing joint polygons, cutting joint blocks and finding key blocks for general free surfaces[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(11): 2161-2170. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200611002.htm
[22] 邬爱清. 基于关键块体理论的岩体稳定性分析方法及其在三峡工程中的应用[J]. 长江科学院院报, 2019, 36(2): 1-7. https://www.cnki.com.cn/Article/CJFDTOTAL-CJKB201902004.htm WU Aiqing. Series methods of analyzing rock mass stability based on key block theory and their applications to Three Gorges Project[J]. Journal of Yangtze River Scientific Research Institute, 2019, 36(2): 1-7. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-CJKB201902004.htm
-
期刊类型引用(5)
1. 郭旭辉，朱鸿鹄，吴冰，高宇新，胡乐乐，曹鼎峰. 基于人工神经网络的黄土含水率光纤被动感测技术研究. 岩土力学. 2025(02): 653-664 . 
百度学术
2. 高磊，袁泽，王勤，高明军. 基于OFDR技术的土体含水率模型研究. 工程地质学报. 2025(01): 96-105 . 百度学术
3. 严辉，林沛元. 深圳市岩溶地层标准贯入击数神经网络模型. 地质科技通报. 2025(02): 305-321 . 百度学术
4. 李杰，朱鸿鹄，吴冰，刘喜凤，王家琛，曹鼎峰，施斌. 下蜀土降雨入渗光纤感测及渗透系数估算研究. 工程地质学报. 2024(02): 601-611 . 百度学术
5. 朱鸿鹄. 工程地质界面：从多元表征到演化机理. 地质科技通报. 2023(01): 1-19 . 百度学术
其他类型引用(0)
计量
- 文章访问数: 464
- HTML全文浏览量: 61
- PDF下载量: 82
- 被引次数: 5
