考虑过渡层多层空间变异性的边坡降雨入渗模型及稳定性分析应用

蒋水华; 元志镕; 刘贤; 黄劲松; 周创兵

doi:10.11779/CJGE20230996

考虑过渡层多层空间变异性的边坡降雨入渗模型及稳定性分析应用 English Version

蒋水华^{1, 2,},
元志镕^{1, 2},
刘贤^1, ,,
黄劲松¹,
周创兵^{1, 2}

1.
南昌大学工程建设学院, 江西南昌 330031
2.
南昌大学流域碳中和教育部工程研究中心, 江西南昌 330031

基金项目:

国家自然科学基金优秀青年基金项目 52222905

国家自然科学基金项目 52179103

国家自然科学基金项目 42272326

江西省自然科学基金项目 20242BAB24001

江西省自然科学基金项目 20232ACB204031

详细信息

作者简介:
蒋水华(1987—)，男，江西九江人，教授，博士生导师，主要从事水工岩土工程可靠度和灾害风险防控方面的研究工作。E-mail: sjiangaa@ncu.edu.cn

通讯作者:
刘贤, E-mail：liux597@mail2.sysu.edu.cn

中图分类号: TU47
计量
- 文章访问数: 384
- HTML全文浏览量: 54
- PDF下载量: 120
出版历程
- 收稿日期: 2023-10-10
- 网络出版日期: 2024-06-10
- 刊出日期: 2025-01-31

Rainfall infiltration model considering spatial variability of multiple layers in transition layer and its application in slope stability analysis

JIANG Shuihua^{1, 2,},
YUAN Zhirong^{1, 2},
LIU Xian^1, ,,
HUANG Jinsong¹,
ZHOU Chuangbing^{1, 2}

1.
School of Infrastructure Engineering, Nanchang University, Nanchang 330031, China
2.
Engineering Research Center of Watershed Carbon Neutralization (Ministry of Education), Nanchang University, Nanchang 330031, China

摘要

摘要: 建立合理的降雨入渗模型是揭示降雨诱发边坡失稳机制及灾害防控的重要前提。传统的Green-Ampt模型未考虑土壤分层及雨水入渗形成的过渡层分布，导致计算的入渗率存在较大偏差，难以适用于空间变异性边坡。提出土层入渗率的分层求解方法，可由不同土层间入渗率大小关系确定过渡层厚度，据此提出考虑土体饱和渗透系数空间变异性边坡降雨入渗分析的改进Green-Ampt模型。进而应用改进模型进行降雨入渗下均质和非均质无限长边坡渗流及稳定性分析，并与传统Green-Ampt模型计算结果和Richards方程数值解加以对比。分析结果表明：相比于传统的Green-Ampt模型，利用改进模型计算的边坡体积含水量分布和稳定性系数与Richards方程数值解更为吻合，可以更好地为非均质边坡降雨入渗分析及降雨型滑坡灾害防控提供理论依据。此外，发现过渡层厚度与土体饱和渗透系数、过渡层顶部入渗率和体积含水量之间存在依赖关系，而与雨水入渗的总深度无直接关系。
- 边坡稳定性分析 /
- 空间变异性 /
- 降雨入渗 /
- 改进Green-Ampt模型 /
- 过渡层
Abstract: Establishing a reasonable rainfall infiltration model is an important prerequisite for revealing the rainfall-induced slope failure mechanism and disaster prevention and control. The traditional Green-Ampt model does not consider the distribution of soil stratification and transition layer formed by rainwater infiltration. This results in a large deviation on the calculated infiltration rate. Thus, the traditional model is difficult to apply to the spatially varying slopes. A method is proposed for calculating the infiltration rate of an arbitrary soil layer. The thickness of transition layer is estimated based on the relationship among the infiltration rates underlying different soil layers. Based on this, an improved Green-Ampt model is proposed to analyze the rainfall infiltration process in the slope considering the spatial variability of saturated hydraulic conductivity of soil. The improved Green-Ampt model is further applied to an infinite slope to analyze its seepage and stability for both homogeneous and heterogeneous soils under the rainfall infiltration. The results obtained from the improved model are systematically compared with those obtained from a traditional Green-Ampt model and the numerical solutions of Richards equation. The results indicate that the distribution of water content and factor of stability calculated from the proposed improved model are more consistent with the numerical solutions of Richards equation than those of the traditional Green-Ampt model. The proposed improved model can lay a solid theoretical foundation for analyzing the rainfall infiltration processes in the heterogeneous slopes, and formulating effective measures for the prevention and control of rainfall-induced landslide disasters. Additionally, it is found that there is a dependence between the thickness of transition layer and the saturated hydraulic conductivity as well as the infiltration rate and volumetric water content at its top, while it is not directly related to the total depth of rainwater infiltration.
- slope stability analysis /
- spatial variability /
- rainfall infiltration /
- improved Green-Ampt model /
- transition layer

HTML全文

0. 引言

注浆是地下工程的灵魂，在地下工程的各类主工法（如明挖、暗挖和盾构）施工中的地层加固与止水两个方面发挥越来越重要的作用^[1-3]。为确保注浆质量，目前工程上注浆多数是靠分段注浆工艺（前进式或后退式）来实现的，在分段注浆设计中，P-Q-T 3个关键技术参数（分别为某段注浆的注浆压力、注浆速率和注浆持续时间等3个设计参数，均为一具体的值，下同）对注浆效果、成本和环境扰动程度有重要影响。但一般情况下在注浆过程中注浆压力和注浆速率（分别以p和q表示）是随时间t变化的，这种情况下如何用一个具体的P和Q值来指导注浆施工？目前的注浆设计计算理论尚不完善，现场施工多依赖于经验，致使设计与施工实际严重脱节，最终造成注浆施工中工程风险和环境风险均很大。这是目前注浆设计施工中普遍存在的问题，尤其是分段注浆压力P的确定上表现更明显，例如：P过低，无法得到满足工程需要的注浆量或扩散半径；有时P过大，造成冒浆、围岩与邻近结构位移超限等问题。

按照不同的介质和工艺，注浆可分为渗透注浆、劈裂注浆、圧密注浆和充填注浆等4类。对于一般的水泥-水玻璃浆液，在中砂粒径以上的地层中以渗透注浆为主。针对渗透注浆，1938年法国学者Maag提出了注浆史上第一个计算公式，即著名的Maag公式。之后，渗透注浆成为整个注浆理论的基础和研究热点，但仍然存在诸多不足，其中之一就是在对浆液黏时变的考虑上。

0.1 传统渗透注浆理论存在的问题

为方便后续讨论，将传统的Magg公式与柱面扩散公式整合在一起表达^[4]（注浆压力统一采用水头高度表达，下同），即

${h_{\text{g}}} = \frac{{{\mu _{{\text{g0}}}}{q_{\text{g}}}}}{{{\mu _{\text{W}}}K}}\int_{{r_{\text{g}}}}^{R(T)} {\frac{{{\text{d}}r}}{{A(r)}}} + {h_0} \text{，}$

(1)

其中， $A(r)=\{\begin{array}{l}4\text{π}({r}^{2}-{r}_{g}^{2})\begin{array}{cc}& \text{(}球面扩散\text{)}\end{array}\\ 2\text{π}L(r-{r}_{\text{g}})\begin{array}{cc}& \text{(}柱面扩散\text{)}\end{array}\end{array}$ ， ${q_{\text{g}}}T =$ $\int_{{r_{\text{g}}}}^{R(T)} {A(r){\text{d}}r \cdot n}$ 。

式中：h_g为注浆压力水头[L]；q_g为注浆速率[L³T^-1]； ${\mu _{{\text{g0}}}}$ 为浆液初始黏滞系数[ML^-1T^-1]，为时间常量； ${\mu _{\text{W}}}$ 为水的黏滞系数[ML^-1T^-1]；K为地层渗透系数[LT^-1]；r_g为注浆孔半径[L]；T为注浆持续时间[T]；R(T)为达到T时的浆液扩散半径[L]；h₀为注浆影响范围（近似为扩散半径R(T)附近）以外的地层中孔隙水压力水头[L]；A(r)为浆液扩散到r时的面积[L²]；r为离注浆孔中心的距离变量[L]；L为柱面扩散段的长度[L]；n为地层孔隙率。

式（1）定量地描述了渗透注浆p和q两个过程量之间的关系（分别对应式中的h_g和q_g），同时，由于推导过程中，黏滞系数按照时间常量考虑的（即整个过程中均等于初始黏滞系数u_g0），因此当q为时间的恒量时，p也可为常量^[4]，这样式（1）中h_g可直接用于P的压力水头取值，为注浆设计提供一定依据^[4]。

由于式（1）中黏滞系数按时间常量考虑，在理论上与在含水层中注水并无本质区别。但实际上，出于注浆加固和止水目的，与化学性质稳定的水体存在着本质区别，浆液的黏滞系数注浆过程中是随时间快速增大的（即所谓“黏时变”，time-dependent viscosity）^[5-13]，即

${\mu _{\text{g}}}^\prime (t) > 0\begin{array}{*{20}{c}} {}&{(t \in [0,{T_{\text{g}}}]} \end{array}) 。$

(2)

式中：μ_g(t)为任意时刻t的黏滞系数，即黏时变函数[ML^-1T^-1]；t为注浆时间[T]；T_g为浆液的凝胶时间[T]；其他符号同前。

浆液黏时变作用势必会引起p-q两个新的时间变量。大量的工程实践也表明，对于渗透注浆这种不改变地层结构的注浆类型，在同一个注浆段内，p一般为t的增函数，而q为一般t的减函数^[14]，即

${h_{\text{g}}}^\prime (t) \geqslant 0{\text{ (}}t \in [0,T]) \text{，}$

(3a)

$q'(t) \leqslant 0{\text{ (}}t \in [0,T]) 。$

(3b)

式中：h_g(t)为任意时刻t的注浆压力p(t)对应的水头高度[L]；q(t)为任意时刻t的注浆速率[L³T^-1]；其他符号同前。

式（1）不能很好地反映式（2），（3）中的实际情况，这是造成传统渗透注浆理论分析的P和实际相差较大的主要原因之一。

0.2 黏时变渗透注浆理论研究现状

随着注浆实践中上述问题越来越突出，相关学者开展了较多的理论研究，这对改进实际渗透注浆设计计算水平起到积极促进作用。目前的研究成果主要是基于黏时变函数某种数学处理而达到对传统注浆理论式（1）一定改进，总体上可以分为如下两大类。

（1）黏时变函数平均值法

该类方法比较直观，直接利用黏时变函数μ_g(t)的平均值直接置换式（1）中的μ_g0，从而得到其注浆压力^[8-10]，即

${h_{{\text{gq}}}} = \frac{{{q_{\text{g}}}\int_0^T {{\mu _{\text{g}}}(t){\text{d}}t} }}{{{\mu _{\text{W}}}KT}}\int_{{r_{\text{g}}}}^{R(T)} {\frac{{{\text{d}}r}}{{A(r)}}} + {h_0} 。$

(4)

式中：h_gq为用黏时变函数平均值法计算的某段注浆压力水头高度[L]。

（2）黏时变函数调和平均值法

该类方法是将式（1）中的μ_g0置换成黏时变函数调和平均值，得到其注浆压力^[11-13]，即

${h_{{\text{gp}}}} = \frac{{{q_{\text{g}}}T}}{{{\mu _{\text{W}}}K\int_0^T {\frac{{{\text{d}}t}}{{{\mu _{\text{g}}}(t)}}} }}\int_{{r_{\text{g}}}}^{R(T)} {\frac{{{\text{d}}r}}{{A(r)}}} + {h_0} 。$

(5)

式中：h_gp为用黏时变函数调和平均值法计算的某段注浆压力水头高度[L]。

在式（4），（5）中，当注浆速率q_g按时间常量考虑时，与之对应的注浆压力h_gq和h_gp也是时间的常量，类似于式（1）的思路，h_gq和h_gp直接作为分段注浆压力设计值。

但是，式（4），（5）在数学推导上存在如下不足：仅考虑到式（2）中的黏时变函数这一个时间变量，而并未充分地考虑到式（3）所表征的p-q双时间变量，同时在概念上也未充分厘清过程量（p，q）与对应具体技术参数（P，Q）的区别及联系。因此，这种对于黏时变的考虑是“不彻底的”，实际上还是一种“常量方法”，无非是将传统注浆理论式（1）中的常量 ${\mu _{{\text{g0}}}}$ 置换成了另一个常量（黏时变函数的平均值或黏时变函数调和平均值），存在理论上的不足。由此得到的P用于设计是否合理，以及上述两种计算方法之间又存在何种关系，也缺乏足够的科学论证和实践的检验。

为系统地揭示对黏时变渗透注浆过程和提高注浆设计计算水平，开展了如下研究工作：①基于定积分原理，提出了p-q双时间变量这一普遍情况下其所对应的P-Q工程定义式，并分析了其物理意义、误差可控性和问题可解性。②基于渗流理论和黏时变理论，建立了黏时变渗透注浆的解析模型，系统地揭示了黏时变渗透注浆复杂过程，探讨了过程量h_g(t)的客观存在性和唯一性。③为克服P显式解推导的困难，基于积分不等式性质和系统的理论推导，得到了P的上、下限通解，讨论了其科学性和普适性。④结合实际工程需要，提出了指数黏时变函数下P的上、下限特解，为黏时变渗透注浆设计实践提供了科学依据和方法基础。

1. 技术参数Q和P的工程定义

为充分考虑p，q双时间变量这一普遍客观实际，同时又便于工程设计应用，取p和q在各自注浆持续时间T内的积分平均值分别定义为该注浆段的P与Q，即

$Q{\text{ = }}\frac{{\int_0^T {q(t){\text{d}}t} }}{T}\begin{array}{*{20}{c}} {}&{{\text{(}}t \in [0,T]} \end{array}) 。$

(6)

式中：Q为分段注浆速率[L³T^-1]，其他符号同前。

${H_{\text{g}}}{\text{ = }}\frac{{\int_0^T {{h_{\text{g}}}(t){\text{d}}t} }}{T}\begin{array}{*{20}{c}} {}&{(t \in [0,T]} \end{array}) 。$

(7)

式中：H_g为分段注浆压力水头高度[L]，其他符号同前。

（1）物理意义

上述工程定义不是单纯数学意义上的平均，而是具有明确物理意义的：式（6）表征了单位时间内该段总注浆量，即该段平均注浆速率；同样，式（7）该段注浆压力的平均冲量。同时，式（6），（7）充分厘清了某注浆段的具体技术参数（P，Q）与相应过程量（p，q）的区别及联系，且反映在数值上，后者是确定前者依据的这一逻辑关系。

（2）误差的可控性

当T较小时，式（6），（7）具有较高精度；当T较大时，可以根据需要将T分为若干个足够小的时间段，然后分别采用类似如式（6），（7）方式表达。

（3）问题的可解性

T和Q可以根据总注浆量要求和浆液性质综合确定，对于P，可以通过建立H_g(Q, T)关系式来求解P。

2. 黏时变注浆过程的解析模型

由于Q与H_g中分别隐含过程量q(t)和h_g(t)，因此，建立H_g(Q, T)的首要工作就是建立q(t)和h_g(t)之间的关系，即黏时变渗透注浆过程解析模型。

2.1 黏时变效应对浆液渗透系数的影响

Bear^[17]在Nutting^[18-19]的工作基础上，提出了任意流体在多孔介质中的渗透系数表达式，利用这一思路可以表达浆液在地层中渗透系数，即

${K_{\text{g}}}(t) = \frac{{{\rho _{\text{g}}}g}}{{{\mu _{\text{g}}}(t)}}k\begin{array}{*{20}{c}} {}&{(t \in [0,{T_{\text{g}}}]} \end{array}) 。$

(8)

式中：K_g(t)为注浆后任意t时刻浆液在地层中的渗透系数[LT^-1]；ρ_g为浆液的密度[ML^-3]；g为重力加速度[LT^-2]；k为地层固有渗透率[L²]；其他符号同前。

类似于式（8），可以得到水在地层中的渗透系数：

$K = \frac{{{\rho _{\text{w}}}g}}{{{\mu _{\text{w}}}}}k 。$

(9)

式中： ${\rho _{\text{w}}}$ 为水的密度[ML^-3]；其他符号同前。

由式（8）和（9）消去渗透率k这一不易确定的参数，K_g(t)可以进一步表达为

${K_{\text{g}}}(t) = \frac{{{\rho _{\text{g}}}{\mu _{\text{W}}}}}{{{\rho _{\text{W}}}{\mu _{\text{g}}}(t)}}K\begin{array}{*{20}{c}} {}&{(t \in [0,{T_{\text{g}}}]} \end{array}) 。$

(10)

2.2 解析模型

（1）物理方程

考虑单孔一般扩散形式下注浆情况，顺着浆液流向取一截平面，如图 1所示。注浆孔半径为r_g，注浆压力为p_g，扩散半径为R(t)（图 1中简写为“R”），地层中孔隙水压力为p₀（相应水头高度为h₀），距离注浆孔中心任意距离r处t时刻的浆液压力为p(r, t)（图 1中简写为“p”，相应水头高度为h(r, t)）。

图 1 单孔注浆浆液流向截面示意图

Figure 1. Section of flow direction of slurry under single hole-grouting

下载: 全尺寸图片幻灯片

根据Darcy定律，可以得到注浆过程中任意断面的流量与压力水头的关系，即黏时变渗透注浆过程的物理方程：

$q(t) = {K_{\text{g}}}(t)\left( { - \frac{{\partial h(r,t)}}{{\partial r}}} \right)A(r)\begin{array}{*{20}{c}} {}&{(t \in [0,T]} \end{array}) 。$

(11)

式中：h(r, t)为注浆扩散半径内任意位置r任意时刻t的浆液压力p(r, t)对应的水头高度[L]，其他符号同前。

将式（10）代入式（11），可以进一步得到

$q(t) = \frac{{{\mu _{\text{w}}}K}}{{{\mu _{\text{g}}}(t)}}\left( { - \frac{{\partial h(r,t)}}{{\partial r}}} \right)A(r)\begin{array}{*{20}{c}} {}&{(t \in [0,T]} \end{array}) 。$

(12)

（2）几何方程

根据任意时刻t注入的浆液量和地层中填充的浆液量相等的原则，结合多孔介质流体力学中有效孔隙率的概念^[17]，可以得到如下浆液扩散几何方程，即

$\int_0^t {q(t)} {\text{d}}t = \int_{{r_{\text{g}}}}^{R(t)} {A(r){\text{d}}r} .{n_{\text{e}}} 。$

(13)

式中：n_e为地层的有效孔隙率，即能与外界发生流体交换且内部互相连通的那部分孔隙体积与地层总体积之比；其他符号同前。

（3）边界条件

根据当前注浆技术领域普遍认识^[4]，注浆扩散半径处的浆液压力为零，即等于该处的地下水水头，可以得到渗流边界条件，即

${\left. h \right|_{r = R(t)}} = {h_0} 。$

(14)

从式（14）可以看出，黏时变渗透注浆的渗流边界为随时间t而发生位置变化的动边界。

式（12）~（14）综合反映了h(r, t)，q(t)，u_g(t)和R(t) 4个时间变量之间的相互制约关系，系统地揭示了黏时变渗透注浆复杂的过程。

2.3 模型解h_g(t)的存在性和唯一性

联立方程（12）和（14），并利用分离变量法，可得到地层中浆液压力表达式，即

$h(r,t){\text{ = }}\frac{{q(t){\mu _{\text{g}}}(t)}}{{{\mu _{\text{W}}}K}}\int_r^{R(t)} {\frac{{{\text{d}}r}}{{A(r)}}} + {h_0}\begin{array}{*{20}{c}} {}&{(t \in [0,T]} \end{array}) 。$

(15)

令式（15）中r=r_g，可以得到模型的解：

${h_{\text{g}}}(t){\text{ = }}\frac{{q(t){\mu _{\text{g}}}(t)}}{{{\mu _{\text{W}}}K}}\int_{{r_{\text{g}}}}^{R(t)} {\frac{{{\text{d}}r}}{{A(r)}}} + {h_0}\begin{array}{*{20}{c}} {}&{(t \in [0,T]} \end{array}) 。$

(16)

式中：h_g(t)为任意时刻注浆压力水头[L]；其他符号同前。

当q(t)和μ_g(t)这两个时间变量已知的条件下，理论上可由式（13）和式（16）两个方程唯一确定注浆压力时间函数h_g(t)，同时并行得到扩散半径时间函数R(t)。

3. P的上、下限通解

鉴于h_g(t)解的存在性和唯一性，理论上可以由此联立式（6），（7）唯一确定关系式H_g(Q，T)。但由于h_g(t)表达式（16）右边除了q(t)这个时间变量外，还隐含另一个时间变量μ_g(t)，根据分部积分原理，H_g(Q，T)很难像h_g(q(t), t)那样显式表述出来。为克服这一推导上的困难，充分利用积分不等式相关性质（详见附录A），进行H_g的上、下限解研究。

3.1 下限解H_g1

对式（16）进行变形，得到

$\frac{{{h_{\text{g}}}(t) - {h_0}}}{{{\mu _{\text{g}}}(t)}}{\text{ = }}\frac{{q(t)}}{{{\mu _{\text{W}}}K}}\int_{{r_{\text{g}}}}^{R{\text{(}}t)} {\frac{{{\text{d}}r}}{{A(r)}}} \begin{array}{*{20}{c}} {}&{(t \in [0,T]} \end{array}) 。$

(17)

对式（17）两边同时在[0，T]上做定积分后除以T，得到

$\frac{{\int_0^T {\frac{{[{h_{\text{g}}}(t) - {h_0}]}}{{{\mu _{\text{g}}}(t)}}} {\text{d}}t}}{T}{\text{ = }}\frac{{\int_0^T {\frac{{q(t)}}{{{\mu _{\text{W}}}K}}\int_{{r_{\text{g}}}}^{R(t)} {\frac{{{\text{d}}r}}{{A(r)}}{\text{d}}t} } }}{T}\begin{array}{*{20}{c}} {}&{(t \in [0,T]} \end{array}) 。$

(18)

由式（2）知，μ_g(t)是t的增函数，因此1/μ_g(t)是t的减函数，由式（3a）知，h_g(t)是t的增函数。因此式（18）左端为单调性相反的两个同域函数积的平均值，根据“单调性相反的两个函数积的平均值不大于二者平均值的积”这一积分不等式性质（证明过程详见附录A），并经过整理，可以得到

$\frac{{\int_0^T {\frac{{[{h_{\text{g}}}(t) - {h_0}]}}{{{\mu _{\text{g}}}(t)}}} {\text{d}}t}}{T} \leqslant \frac{{\int_0^T {[{h_{\text{g}}}(t) - {h_0}]} {\text{d}}t\int_0^T {\frac{{{\text{d}}t}}{{{\mu _{\text{g}}}(t)}}} }}{{{T^2}}} 。$

(19)

式中： $“=”仅当{{h}^{\prime }}_{\text{g}}(t)\equiv 0时成立$ 。

将式（7）代入式（19），并联立式（18），经过整理，得到如下不等式，即

${H_{\text{g}}} \geqslant \frac{{\frac{T}{{\int_0^T {\frac{{{\mu _{\text{w}}}}}{{{\mu _{\text{g}}}(t)}}{\text{d}}t} }}}}{K}\int_0^T {\frac{{q(t)}}{T}} \int_{{r_{\text{g}}}}^{R(t)} {\frac{{{\text{d}}r}}{{A(r)}}{\text{d}}t} + {h_0} 。$

(20)

式中： $“=”仅当{{h}^{\prime }}_{\text{g}}(t)\equiv 0时成立$ 。

显然，式（20）的右端为 ${h'_{\text{g}}}(t) \equiv 0$ （即恒压注浆）下成立的H_g下限解，为方便后续讨论，利用分部积分原理，并结合几何方程式（13），对式（20）右端做进一步简化表达，即

${H_{{\text{g1}}}} = \frac{1}{K}[{F_1}(T) - {G_1}(T)] + {h_0}\begin{array}{*{20}{c}} {}&{(t \in [0,T]} \end{array}) 。$

(21)

其中：

${F_1}(Q,T) = {\alpha _1}Q\int_{{r_{\text{g}}}}^{R(T)} {\frac{{{\text{d}}r}}{{A(r)}}} \text{，}$

${G_1}(Q,T) = \frac{{{\alpha _1}{n_{\text{e}}}}}{T}\int_0^T {[\int_{{r_{\text{g}}}}^{R(t)} {A(r){\text{d}}r]} \frac{{R'(t)}}{{A[R(t)]}}{\text{d}}t} \text{，}$

${\alpha _1} = \frac{T}{{\int_0^T {\frac{{{\mu _{\text{w}}}}}{{{\mu _{\text{g}}}(t)}}{\text{d}}t} }} 。$

式中：H_g1为H_g的下限解，在恒压注浆条件下成立[L]；F₁(Q, T)为恒压注浆下定边界影响因子[L²T^-1]；G₁(Q, T)为恒压注浆条件下动边界影响因子[L²T^-1]；α₁为恒压注浆条件下的黏时变影响系数，无量纲量；其他符号同前。

3.2 上限解H_g2

对式（17）两边同时在[0, T]上做定积分后除以T，得到

$\begin{gathered} \frac{\int_0^T h_{\mathrm{g}}(t) \mathrm{d} t}{T}=\frac{\int_0^T q(t) \mu_{\mathrm{g}}(t) \int_{r_{\mathrm{g}}}^{R(t)} \frac{\mathrm{d} r}{A(r)} \mathrm{d} t}{\mu_{\mathrm{W}} K T}+h_0 \\ t \in[0, T] \quad \text { 。 } \end{gathered}$

(22)

由式（3b）知，q(t)是t的减函数，由式（2）知μ_g(t)是t的增函数，而 $\int_{{r_{\text{g}}}}^{R(t)} {\frac{{{\text{d}}r}}{{A(r)}}}$ 是t的增函数，因此式（22）右端第一项中的 $\frac{{\int_0^T {q(t)[{\mu _{\text{g}}}(t)\int_{{r_{\text{g}}}}^{R(t)} {\frac{{{\text{d}}r}}{{A(r)}}]} {\text{d}}t} }}{T}$ 为单调性相反的两个函数积的平均值，根据类似于式（19）的方法，可以得到如下积分不等式，即

$\begin{aligned} & \frac{\int_0^T q(t)\left[\mu_{\mathrm{g}}(t) \int_{r_{\mathrm{g}}}^{R(t)} \frac{\mathrm{d} r}{A(r)}\right] \mathrm{d} t}{T} \\ & \leqslant \frac{\int_0^T q(t) \mathrm{d} t \int_0^T\left[\mu_{\mathrm{g}}(t) \int_{r_{\mathrm{g}}}^{R(t)} \frac{\mathrm{d} r}{A(r)}\right] \mathrm{d} t}{T^2} 。 \end{aligned}$

(23)

式中： $"="仅当{{h}^{\prime }}_{\text{g}}(t)\equiv 0时成立$ ， $"="仅当{q}^{\prime }(t)\equiv$ $0时$ $成立$ 。

将式（22）代入式（23）中，并联立式（6），（7），经过整理，得到

$\begin{gathered} H_{\mathrm{g}} \leqslant \frac{Q}{\mu_{\mathrm{W}} K T} \int_0^T\left[\mu_{\mathrm{g}}(t) \int_{r_{\mathrm{g}}}^{R(t)} \frac{\mathrm{d} r}{A(r)}\right] \mathrm{d} t+h_0 \\ t \in[0, T] 。 \end{gathered}$

(24)

式中： $“=”仅当{q}^{\prime }(t)\equiv 0时成立$ 。

显然，式（24）右端为 $q'(t) \equiv 0$ 条件下（即恒速注浆）H_g的上限解，为方便后续讨论，对式（24）右端做进一步简化表达，即

${H_{{\text{g2}}}} = \frac{1}{K}[{F_2}(T) - {G_2}(T)] + {h_0} \text{，}$

(25)

其中， ${F_2}(Q,T) = {\alpha _2}Q\int_{{r_{\text{g}}}}^{R(T)} {\frac{{{\text{d}}r}}{{A(r)}}}$ ，

${G_2}(Q,T) = \frac{Q}{{{\mu _{\text{W}}}T}}\int_0^T {[\int_0^t {{\mu _{\text{g}}}{\text{(}}t{\text{)d}}t]} \frac{{R'(t){\text{d}}t}}{{A[R(t)]}}} \text{，}$

${\alpha _2} = \frac{{\int_0^T {\frac{{{\mu _{\text{g}}}(t)}}{{{\mu _{\text{w}}}}}{\text{d}}t} }}{T} 。$

式中：H_g2为H_g的上限解，在恒速注浆条件下成立[L]； ${\alpha _2}$ 为恒速注浆条件下黏时变影响系数，无量纲量；F₂(Q, T)为恒速注浆下定边界影响因子[L²T^-1]；G₂(Q, T)为恒速注浆条件下动边界影响因子[L²T^-1]；其他符号同前。

上述下限解H_g1和上限解H_g2的推导分别基于恒压和恒速的各种理想假设条件，而实际注浆中压力p和速率q是同时变化的，因此下限解H_g1和上限解H_g2主要工程意义是为满足一定Q和T所需要H_g的最小值与最大值。

3.3 关于解的讨论

式（21），（25）在推导过程中出于数学表达的便利，引入了定边界影响因子和动边界影响因子。下面从二者表达式上，进一步分析各自物理意义。

（1）定边界影响因子

定边界影响因子未含扩散半径的时间变化过程量R(t)，而仅含有最终的扩散半径R(T)（参见F₁₍Q，T)和F₂(Q，T)的表达式）。相当于注浆过程中地层中浆液压力边界自始至终都在R(T)的位置，这与含水层中注水过程中的水力传导边界位置类似^[19]。但实际上，浆液扩散范围R(t)以外是不存在浆液压力的（式（14）），因此在R(T)域上并非任意时刻注浆压力传导都是连续过程，若整个T内统一用定边界R(T)来计算分段注浆压力上、下限解，显然均是偏大的。

（2）动边界影响因子

动边界影响因子含有扩散半径的时间变化过程量R(t)（参见G₁(Q, T)和G₂(Q, T)表达式），反映了注浆过程中由于浆液扩散过程而对定边界下的分段注浆压力的一种“消散作用”。相同条件下 $R'(t)$ 越大，G₁(Q, T)和G₂(Q, T)也越大，反映了浆液扩散速率越大，对分段注浆压力消散作用越明显。

（3）关于 $R'(t)$ =0的一个特例

假设如下一个极端的特例，即

$R'(t) = 0{\text{ (}}t \in [0,7]) 。$

(26)

则很容易得到上、下限解分别蜕化为如下两个表达式，即

$H_{{\text{g1}}}^* = \frac{{{\alpha _1}}}{K}\int_{{r_{\text{g}}}}^{R(T)} {\frac{{{\text{d}}r}}{{A(r)}}} + {h_0}{\text{ (}}t \in [0,T]) \text{，}$

(27)

$H_{{\text{g2}}}^* = \frac{{{\alpha _2}}}{K}\int_{{r_{\text{g}}}}^{R(T)} {\frac{{{\text{d}}r}}{{A(r)}}} + {h_0}{\text{ (}}t \in [0,T]) 。$

(28)

对比既有黏时变研究成果，容易发现，在数学形式上，下限解 $H_{{\text{g1}}}^*$ 与黏时变函数调和平均值法（式（5））完全一致^[11-13]，而上限解 $H_{{\text{g2}}}^*$ 与黏时变函数平均值法（式（4））完全一致^[8-10]。因此，既有研究成果可以分别理解为上、下限通解忽略各自动边界影响因子下的一种近似简化形式。但式（26）仅仅是一种纯数学上假设，而缺乏足够的科学依据和工程背景。因此，在进行一般情况下P的上、下限解分析时，需要同时考虑上述两类边界影响因子。

4. 指数黏时变函数下的上、下限特解

在前述通解基础上，针对目前常用的指数黏时变函数条件下的特解做进一步研究。

根据Honma等^[5]研究成果，指数黏时变函数可表达为

${\mu _{\text{g}}}(t) = {\mu _{{\text{g0}}}} \cdot {{e} ^{at}}\begin{array}{*{20}{c}} {}&{(t \in [0,T]} \end{array}) 。$

(29)

式中：a为凝胶常数[T^-1]，反映浆液黏滞系数随注浆时间t变化的快慢程度。

根据式（29），结合式（21），（25），可以得到指数黏时变函数下的恒压、恒速两种条件下黏时变影响系数，即

$\alpha _{\text{1}}^{\text{e}} = \frac{{{\mu _{{\text{g0}}}}aT\exp (aT)}}{{{\mu _{\text{w}}}[\exp (aT) - 1]}} \text{，}$

(30a)

$\alpha _{\text{2}}^{\text{e}} = \frac{{{\mu _{{\text{g0}}}}}}{{{\mu _{\text{w}}}aT}}[\exp (aT) - 1 。$

(30b)

式中： $\alpha _{\text{1}}^{\text{e}}$ 为指数黏时变函数下恒压注浆黏时变影响系数； $\alpha _{\text{2}}^{\text{e}}$ 为指数黏时变函数下恒速注浆黏时变影响系数；其他符号同前。

以式（30）所示的两类指数函数黏时变影响系数为基础，根据前述上下限通解研究成果，结合式（1），进行进一步推导，得到球面扩散和柱面扩散这两个常见浆液扩散模式下的下、上限特解如表 1和表 2。

表 1 球面扩散模式下P的上、下限特解

Table 1. Special limit solution of P under spherial diffusion mode

解的类型	解的表达式	符号说明
下限解	$\begin{gathered} H_{{\text{g1}}}^{\text{s}} = \frac{{{\mu _{{\text{g0}}}}aT\exp (aT)}}{{{\mu _{\text{w}}}[\exp (aT) - 1]K}} \cdot \hfill \\ \left\{ {\frac{Q}{{4{\text{π }}}}\left[ {\frac{1}{{{r_g}}} - \frac{1}{{{R^s}(T)}}} \right]} \right. - \hfill \\ \left. {\frac{{{n_{\text{e}}}{{[{R^{\text{s}}}(T) - {r_{\text{g}}}]}^2}[{R^s}(T) + 2{r_{\text{g}}}]}}{{6{R^{\text{s}}}(T)T}}} \right\} + {h_0} \hfill \\ \end{gathered}$	$\begin{aligned} &R^{\mathrm{s}}(T)=\\ &\sqrt[3]{\frac{Q T}{4 \pi n_{\mathrm{e}}}+r_{\mathrm{g}}^3} \end{aligned}$
上限解	$H_{{\text{g2}}}^{\text{s}} = \frac{{{\mu _{{\text{g0}}}}Q}}{{4{\text{π }}{\mu _{\text{W}}}TK}} \cdot$ $\left\{ {\frac{{[\exp (aT) - 1]}}{{a{r_{\text{g}}}}} - \int_0^T {\frac{{{\text{exp(}}a{\text{t)}}}}{{{{\bar R}^{\text{s}}}{\text{(}}t{\text{)}}}}{\text{d}}t} } \right\}$	$\begin{aligned} &\bar{R}^{\mathrm{s}}(t)=\\ &\sqrt[3]{\frac{Q t}{4 \pi n_{\mathrm{e}}}+r_{\mathrm{g}}^3} \end{aligned}$

下载: 导出CSV

| 显示表格

表 2 柱状扩散模式下P的上、下限特解

Table 2. Special limit solution of P under cylindrical diffusion mode

解的类型	解的表达式	符号说明
下限解	$\begin{gathered} H_{{\text{g1}}}^{\text{c}} = \frac{{{\mu _{{\text{g0}}}}aT\exp (aT)}}{{2{\mu _{\text{w}}}K[\exp (aT) - 1]}}\left\{ {\left( {\frac{Q}{{{\text{π }}L}} + \frac{{{n_{\text{e}}}}}{T}r_{\text{g}}^{\text{2}}} \right) \cdot } \right. \hfill \\ \ln \frac{{{R^{\text{c}}}(T)}}{{{r_{\text{g}}}}} - \left. {\frac{{{n_{\text{e}}}}}{{2T}}[{R^{\text{c}}}{{(T)}^2} - r_{\text{g}}^{\text{2}}]} \right\} \hfill \\ \end{gathered}$	$\begin{gathered} {R^{\text{c}}}(T) = \hfill \\ \sqrt {\frac{{QT}}{{{\text{π }}L{n_{\text{e}}}}} + r_{\text{g}}^{\text{2}}} \hfill \\ \end{gathered}$
上限解	$\begin{gathered} H_{{\text{g2}}}^{\text{c}} = \frac{{{\mu _{{\text{g0}}}}Q}}{{2{\text{π }}KL{\mu _{\text{w}}}T}}\left[ {\int_0^T {\exp (at){\text{ln}}{{\bar R}^{\text{c}}}(t)} {\text{d}}t - \begin{array}{{20}{c}} {} \\ {} \end{array}} \right. \hfill \\ \left. {\begin{array}{{20}{c}} {} \\ {} \end{array}\frac{{\exp (aT) - 1}}{a}{\text{ln}}{r_{\text{g}}}} \right] \hfill \\ \end{gathered}$	$\begin{gathered} {{\bar R}^{\text{c}}}(t) = \hfill \\ \sqrt {\frac{{Qt}}{{{\text{π }}L{n_{\text{e}}}}} + r_{\text{g}}^{\text{2}}} \hfill \\ \end{gathered}$

下载: 导出CSV

| 显示表格

上述两类扩散模式的特解表明，下限解均为显式解，工程应用比较便利；而上限解略为复杂，主要是均含有特殊积分，但可以借助数值积分手段予以实现。

5. 结论

浆液的黏时变效应对渗透注浆有重要影响，传统的注浆理论没有考虑到这一点，用于设计中在注浆压力分析存在较大误差。而目前黏时变注浆理论更多的是一种半经验数学处理，缺乏足够的科学内涵和工程验证。基于系统理论研究，得到以下4点结论。

（1）基于定积分原理的P-Q工程定义从概念上厘清了注浆设计中的具体技术参数值与注浆施工过程量p-q之间的区别与联系，同时具有较好的误差可控性和问题可解性。

（2）黏时变渗透注浆解析模型系统揭示了h(r, t)、q(t)、u_g(t)和R(t)等4个时间变量之间的相互制约关系，较全面而精确地描述了黏时变渗透注浆复杂的过程，形成了黏时变注浆理论研究的重要基础。

（3）基于P-Q工程定义、渗透注浆过程的解析模型和积分不等式等科学理论方法的严格推导，得到了一般黏时变渗透注浆条件下的分段注浆压力P的下限通解（恒压条件下）与上限通解（恒速条件下）这一普遍意义成果，解的讨论进一步说明了其科学性与普适性。

（4）指数黏时变函数下的分段注浆压力的上、下限特解表达式较为简单实用，为该类型浆液黏时变函数下渗透注浆工程实践提供了科学依据和技术支撑。

附录A单调性相反的两个同域函数积分不等式性质及其证明

积分不等式性质：函数f(x)与g(x)在 $x \in [a,b]$ 上二阶连续可导，且f(x)与g(x)的单调性相反，则该两个函数在[a, b]上积的平均值不大于二者的平均值的积，即

$\begin{aligned} &\frac{\int_a^b f(t) g(t) \mathrm{d} t}{b-a} \leqslant \frac{\int_a^b f(t) \mathrm{d} t}{b-a} \cdot \frac{\int_a^b g(t) \mathrm{d} t}{b-a} \\ "= &\text { "当且仅当 } f^{\prime}(t) \equiv 0 \text { 或 } g^{\prime}(t) \equiv 0 \text { 时成立， } \end{aligned}$

(A)

证明：构造函数

$\begin{gathered} F(x)=(x-a) \int_a^x f(t) g(t) \mathrm{d} t- \\ \int_a^x f(t) \mathrm{d} t \int_a^x g(t) \mathrm{d} t \quad x \in[a, b], \end{gathered}$

(A-1)

对式(A-1)两边对x求导数，即

$$ \begin{aligned} &F^{\prime}(x)=(x-a) \cdot f(x) g(x)+ \\ &\quad \int_a^x f(t) g(t) \mathrm{d} t-f(x) \int_a^x g(t) \mathrm{d} t-g(x) \int_a^x f(t) \mathrm{d} t \quad \end{aligned} $ 。$

(A-2)

对式(A-2)两边继续对x求导数并经过整理，可以得到函数F(x)的二阶导数，即

$\begin{aligned} &F^{\prime \prime}(x)=f^{\prime}(x)[g(x)(x-a)- \\ &\left.\quad \int_a^x g(t) \mathrm{d} t\right]+g^{\prime}(x)\left[(x-a) f(x)-\int_a^x f(t) \mathrm{d} t\right] 。 \end{aligned}$

(A-3)

根据已知条件f(t)与g(t)单调性相异，考虑问题对称性和方便讨论，不妨令

$f'(t) \geqslant 0 \text{，}$

(A-4a)

$g'(t) \leqslant 0 。$

(A-4b)

由式(A-4b)知，g(t)为时间变量t单调减函数，根据积分中值定理，则有

$\begin{aligned} &g(x)(x-a)-\int_a^x g(t) \mathrm{d} t \leqslant 0 \\ &"=\text { "当且仅当 } g^{\prime}(t) \equiv 0 \text { 时成立。 } \end{aligned}$

(A-5)

将式(A-5)和式(A-4a)相乘，可以得到式(A-3)右端第一项非正，即

$\begin{aligned} &f^{\prime}(x)\left[g(x)(x-a)-\int_a^x g(t) \mathrm{d} t\right] \leqslant 0 \\ &"=\text { "当且仅当 } f^{\prime}(t) \equiv 0 \text { 或 } g^{\prime}(t) \equiv 0 \text { 时成立。 } \end{aligned}$

(A-6)

由式(A-4a)知，f(t)为时间变量t单调增函数，根据积分中值定理，则有

$\begin{aligned} &(x-a) f(x)-\int_a^x f(t) \mathrm{d} t \geqslant 0 \\ &\quad"=\text { "当且仅当 } f^{\prime}(t) \equiv 0 \text { 时成立。} \end{aligned}$

(A-7)

由式(A-7)与式(A-4b)相乘，得到式(A-3)右端第二项非正，即

$$ \begin{aligned} &g^{\prime}(x)\left[(x-a) f(x)-\int_a^x f(t) \mathrm{d} t\right] \leqslant 0 \\ &\quad"=\text { "当且仅当 } f^{\prime}(t) \equiv 0 \text { 或 } g^{\prime}(t) \equiv 0 \text { 时成立 } 。 \end{aligned} $ 。$

(A-8)

由式(A-6)与式(A-8)相加，可以得到式(A-3)非正，即

$\begin{array}{cc}{F}^{″}(x)\le 0& "="当且仅当{f}^{\prime }(t)\equiv 0或{g}^{\prime }(t)\equiv 0时成立。\end{array}$

(A-9)

由式(A-9)并根据二阶导数的定义可知， $F'(x)$ 为变量x的单调减函数，则有

$\begin{aligned} &F^{\prime}(x) \leqslant F^{\prime}(a)=0 \\ &\quad " \leqslant \text { "中的" }=\text { "当且仅当 } f^{\prime}(t) \equiv 0 \text { 或 } g^{\prime}(t) \equiv 0 \text { 时成立。 } \end{aligned}$

(A-10)

由式(A-10)并根据一阶导数性质可知，F(x)为变量x的单调减函数，则有

$\begin{aligned} &F(x) \leqslant F(a)=0 \\ &\quad " \leqslant \text { "中的" = "当且仅当 } f^{\prime}(t) \equiv 0 \text { 或 } g^{\prime}(t) \equiv 0 \text { 时成立。 } \end{aligned}$

(A-11)

联立式(A-11)和式(A-1)，并令x=b，通过整理可以得到下式：

$\begin{aligned} &\frac{\int_a^b f(t) g(t) \mathrm{d} t}{b-a} \leqslant \frac{\int_a^b f(t) \mathrm{d} t}{b-a} \cdot \frac{\int_a^b g(t) \mathrm{d} t}{b-a} \\ &"=\text { "当且仅当 } f^{\prime}(t) \equiv 0 \text { 或 } g^{\prime}(t) \equiv 0 \text { 时成立 } \end{aligned}$

(A-12)

式(A)得证。

图 1 Brook-Corey模型描述的土-水特征曲线

Figure 1. Soil-water characteristic curve described using Brook-Corey model

下载: 全尺寸图片幻灯片

图 2 提出模型的体积含水量分布图

Figure 2. Distribution of volumetric water content by proposed model

下载: 全尺寸图片幻灯片

图 3 无限长边坡计算模型

Figure 3. Diagram of an infinite slope mode

下载: 全尺寸图片幻灯片

图 4 不同降雨历时下边坡体积含水量分布的比较

Figure 4. Comparison of distributions of volumetric water content in slope under different rainfall durations

下载: 全尺寸图片幻灯片

图 5 不同降雨历时下过渡层体积含水量和入渗率关系的比较

Figure 5. Comparison of relationships between volumetric water content and infiltration rate of transition layer for different rainfall durations

下载: 全尺寸图片幻灯片

图 6 过渡层厚度逐层推进计算示意图

Figure 6. Diagram for calculating thickness of transitional soil layer through layer by layer

下载: 全尺寸图片幻灯片

图 7 饱和层基质吸力变化规律示意图

Figure 7. Diagram for variation of matric suction in saturated soil layer

下载: 全尺寸图片幻灯片

图 8 不同降雨历时下边坡体积含水量分布的比较

Figure 8. Comparison of distributions of volumetric water content of slope under different rainfall durations

下载: 全尺寸图片幻灯片

图 9 不同降雨历时下边坡稳定性系数分布的比较

Figure 9. Comparison of factors of stability of slope under different rainfall durations

下载: 全尺寸图片幻灯片

图 10 不同降雨历时下边坡体积含水量分布的比较

Figure 10. Comparison of distributions of volumetric water content of slope under different rainfall durations

下载: 全尺寸图片幻灯片

图 11 不同降雨历时下边坡稳定性系数分布的比较

Figure 11. Comparison of factors of stability of slope under different rainfall durations

下载: 全尺寸图片幻灯片

表 1 土体物理力学参数取值

Table 1 Values of physical mechanical parameters of soil

计算参数	单位	量纲	取值
饱和渗透系数k_s	cm/h	LT^-1	0.3
有效内摩擦角 $\varphi '$	°	1	28
土体干重度 ${\gamma _{\text{t}}}$	kN/m³	ML^-2T^-2	16.217
饱和体积含水量 ${\theta _{\text{s}}}$		1	0.335
初始体积含水量 ${\theta _{\text{i}}}$		1	0.148
残余体积含水量 ${\theta _{\text{r}}}$		1	0.068
有效黏聚力 $c'$	kPa	L^-1MT^-2	5
进气值 ${\psi _{\text{b}}}$	kPa	L^-1MT^-2	2.752
初始基质吸力 $\psi ({\theta _{\text{i}}})$	kPa	L^-1MT^-2	120.356
土体孔隙分布特征参数λ		1	0.319
注： $\psi ({\theta _{\text{i}}})$ 是将初始体积含水量代入式（4）土水特征曲线函数中计算得到。

下载: 导出CSV

参考文献(21)

[1]	FROUDE M J, PETLEY D N. Global fatal landslide occurrence from 2004 to 2016[J]. Natural Hazards and Earth System Sciences, 2018, 18(8): 2161-2181. doi: 10.5194/nhess-18-2161-2018
[2]	蒋水华, 刘贤, 黄劲松, 等. 多层非均质边坡降雨入渗分析的改进Green-Ampt模型[J]. 岩土工程学报, 2024, 46(6): 1177-1186. doi: 10.11779/CJGE20230205 JIANG Shuihua, LIU Xian, HUANG Jinsong, et al. An improved Green-Ampt model for rainfall infiltration analysis of multi-layered heterogeneous soil slopes[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(6): 1177-1186. (in Chinese) doi: 10.11779/CJGE20230205
[3]	王文焰, 汪志荣, 王全九, 等. 黄土中Green-Ampt入渗模型的改进与验证[J]. 水利学报, 2003, 34(5): 30-34. WANG Wenyan, WANG Zhirong, WANG Quanjiu, et al. Improvement and evaluation of the Green-Ampt model in loess soil[J]. Journal of Hydraulic Engineering, 2003, 34(5): 30-34. (in Chinese)
[4]	彭振阳, 黄介生, 伍靖伟, 等. 基于分层假设的Green-Ampt模型改进[J]. 水科学进展, 2012, 23(1): 59-66. PENG Zhenyang, HUANG Jiesheng, WU Jingwei, et al. Modification of Green-Ampt model based on the stratification hypothesis[J]. Advances in Water Science, 2012, 23(1): 59-66. (in Chinese)
[5]	王雪冬, 李世宇, 孙延峰, 等. 考虑非饱和浸润层厚度和累积入渗量的改进Green-Ampt模型[J]. 水文地质工程地质, 2021, 48(6): 64-71. WANG Xuedong, LI Shiyu, SUN Yanfeng, et al. An improved Green-Ampt model for rainfall infiltration in the inner dumping site of an open pit coal mine[J]. Hydrogeology & Engineering Geology, 2021, 48(6): 64-71. (in Chinese)
[6]	黄良誉, 何廷全, 周成, 等. 边坡植被恢复中植被水泥土Green-Ampt入渗模型的改进与应用[J]. 岩土工程学报, 2022, 44(增刊1): 183-188. doi: 10.11779/CJGE2022S1033 HUANG Liangyu, HE Tingquan, ZHOU Cheng, et al. Improvement and application of Green-Ampt infiltration model of vegetation cement soil in slope vegetation restoration[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(S1): 183-188. (in Chinese) doi: 10.11779/CJGE2022S1033
[7]	ZHANG L L, ZHANG L M, TANG W H. Rainfall-induced slope failure considering variability of soil properties[J]. Géotechnique, 2005, 55(2): 183-188. doi: 10.1680/geot.2005.55.2.183
[8]	DOU H Q, HAN T C, GONG X N, et al. Probabilistic slope stability analysis considering the variability of hydraulic conductivity under rainfall infiltration-redistribution conditions[J]. Engineering Geology, 2014, 183: 1-13. doi: 10.1016/j.enggeo.2014.09.005
[9]	蒋水华, 刘贤, 黄发明, 等. 考虑多参数空间变异性的降雨入渗边坡失稳机理及可靠度分析[J]. 岩土工程学报, 2020, 42(5): 900-907. doi: 10.11779/CJGE202005012 JIANG Shuihua, LIU Xian, HUANG Faming, et al. Failure mechanism and reliability analysis of soil slopes under rainfall infiltration considering spatial variability of multiple soil parameters[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(5): 900-907. (in Chinese) doi: 10.11779/CJGE202005012
[10]	马世国. 强降雨条件下基于Green-Ampt入渗模型的无限边坡稳定性研究[D]. 杭州: 浙江大学, 2014. MA Shiguo. Study on Infinite Slope Stability under Heavy Rainfall Based on Green-Ampt Infiltration Model[D]. Hangzhou: Zhejiang University, 2014. (in Chinese)
[11]	GREEN W H, AMPT G A. Studies on soil physics: flow of air and water through soils[J]. Journal of Agricultural Science, 1911, 4(1): 1-24.
[12]	雷志栋. 土壤水动力学[M]. 北京: 清华大学出版社, 1988. LEI Zhidong. Soil Hydrodynamics[M]. Beijing: Tsinghua University Press, 1988. (in Chinese)
[13]	FREDLUND D G, RAHARDJO H. Soil Mechanics for Unsaturated Soils[M]. New York: Wiley, 1993.
[14]	BROOKS R H, COREY A T. Hydraulic Properties of Porous Media[M]. Fort Collins: Colorado State University, 1964.
[15]	李强, 贾森, 李鑫, 等. 考虑非饱和浸润区的改进Green-Ampt模型[J]. 岩土力学, 2022, 43(12): 3484-3492, 3502. LI Qiang, JIA Sen, LI Xin, et al. An improved Green-Ampt model considering unsaturated infiltration zone[J]. Rock and Soil Mechanics, 2022, 43(12): 3484-3492, 3502. (in Chinese)
[16]	RICHARDS L A. Capillary conduction of liquids through porous mediums[J]. Physics, 1931, 1(5): 318-333.
[17]	SIMUNEK J, VAN GENUCHTEN M T, SEJNA M. The Hydrus-1D software package for simulating the movement of water, heat, and multiple solutes in variably saturated media, Version 4.16, HYDRUS Software Series 3[R]. California: Department of Environmental Sciences, University of California Riverside, 2013.
[18]	CHO S E. Failure distribution analysis of shallow landslides under rainfall infiltration based on fragility curves[J]. Landslides, 2020, 17(1): 79-91.
[19]	LU N, GODT J W, WU D T. A closed-form equation for effective stress in unsaturated soil[J]. Water Resources Research, 2010, 46(5): W05515.
[20]	陈正汉, 郭楠. 非饱和土与特殊土力学及工程应用研究的新进展[J]. 岩土力学, 2019, 40(1): 1-54. CHEN Zhenghan, GUO Nan. New developments of mechanics and application for unsaturated soils and special soils[J]. Rock and Soil Mechanics, 2019, 40(1): 1-54. (in Chinese)
[21]	赵振兴, 何建京. 水力学[M]. 2版. 北京: 清华大学出版社, 2010. ZHAO Zhenxing, HE Jianjing. Hydraulics[M]. 2nd ed. Beijing: Tsinghua University Press, 2010. (in Chinese)

施引文献(7)

期刊类型引用(2)

1.	刘要来，王堡生，周红波，赵二峰，李章寅. 基于凸集比例因子和WOA-Kriging模型的重力坝非概率可靠性分析. 长江科学院院报. 2025(03): 164-170 . 百度学术
2.	魏博文，张升，袁冬阳，徐富刚. 基于概率—模糊—区间混合模型和改进分枝限界法的重力坝可靠性分析方法. 水利学报. 2022(12): 1476-1489 . 百度学术