红黏土覆盖层水力参数演变规律与防渗设计讨论

焦卫国; 刘振男; 季永新; 张玥; 贺明卫; 廖化荣

doi:10.11779/CJGE202201003

红黏土覆盖层水力参数演变规律与防渗设计讨论 English Version

焦卫国^{1, 3,},
刘振男^{1, 3},
季永新²,
张玥^{1, 3},
贺明卫^{1, 3},
廖化荣^{1, 3}

1.
贵州理工学院土木工程学院, 贵州贵阳 550003
2.
中国建筑第四工程局有限公司工程技术研究院(技术中心), 广东广州 510665
3.
贵州理工学院乡村建设工程技术研究中心, 贵州贵阳 550003

基金项目:

国家自然科学基金项目 52069005

黔科合基础项目 [2019]1143，ZK[2021]295

中建四局科技项目 1700202100019003

贵州理工高层次人才项目 XJGC20190912

详细信息

作者简介:
焦卫国(1983—)，男，博士，从事非饱和土力学、环境土工方面的研究和教学工作。E-mail: 805810460@qq.com

中图分类号: TU43
计量
- 文章访问数: 298
- HTML全文浏览量: 23
- PDF下载量: 131
出版历程
- 收稿日期: 2020-07-06
- 网络出版日期: 2022-09-22
- 刊出日期: 2021-12-31

Evolution laws of hydraulic parameters of red clay covers and design of seepage prevention

JIAO Wei-guo^{1, 3,},
LIU Zhen-nan^{1, 3},
JI Yong-xin²,
ZHANG Yue^{1, 3},
HE Ming-wei^{1, 3},
LIAO Hua-rong^{1, 3}

1.
School of Civil Engineering, Guizhou Institute of Technology, Guiyang 550003, China
2.
Engineering Technology Research Institute (Technology Center) of Chinaconstruction Fourth Engineering Division Co., Ltd., Guangzhou 510665, China
3.
Rural construction Engineering Technology Research Center of Guizhou Institute of Technology, Guiyang 550003, China

摘要

摘要: 以红黏土为对象，通过室内单元体并结合填埋场现场原位试验，从建设施工和建成后长期服役两个不同时间尺度对土质覆盖层水力参数开展了5 a的跟踪监测；分析对比了从实验室到现场、从建设施工到建成长期服役水力参数的劣化衰减规律。结果表明：①5 a长期服役中，无植被红黏土覆盖层入渗系数从10^-7 cm/s增大到10^-3 cm/s，增大4个数量级；有植被条件入渗系数从10^-7 cm/s增大到10^-6 cm/s，增大1个数量级。②自然气候长期服役覆盖层中大孔隙逐渐增多，降雨雨强过大（如暴雨）导致土质覆盖层防渗能力和储水能力下降。工程设计中应根据当地降雨参数进行适当放大调整。③自然气候下经历反复吸、脱湿滞回循环，土质覆盖层的储水能力与降雨前土层吸湿起点（或初始吸力）和吸湿路径有关。采用室内吸湿起点初始吸力-1500 kPa的主吸湿曲线设计偏于保守，结果相对安全。
- 尾矿堆场 /
- 填埋场 /
- 封场覆盖 /
- 植被土 /
- 长期服役 /
- 水力参数 /
- 植被根系 /
- 防渗设计
Abstract: Through the laboratory unit and in-situ tests in landfill site, the hydraulic parameters of red clay soil covers have been monitored for more than five years on two different time scales of construction and long-term service. The degradation laws of hydraulic parameters from laboratory to field as well as from construction to long-term service are analyzed and compared. The results show that: (1) During the five years of long-term service, the infiltration coefficient of red clay covers without vegetation changes from 10^-7 cm/s to 10^-3 cm/s, with the variation fluctuated by 5 orders of magnitude, and the infiltration coefficient with vegetation is 2 orders of magnitude from 10^-7 cm/s to 10^-6 cm/s. (2) The macropores in the covers gradually increase in the long-term service of natural climate, and the excessive rainfall (such as rainstorm) leads to the decrease of water storage capacity of soil overburden. In engineering design, it should be adjusted according to the local rainfall parameters. (3) With the repeated cycles of moisture absorption and desorption, the water storage capacity is related to the initial matrix suction and path of moisture absorption before rainfall. The design of the main moisture absorption curve with the starting point of laboratory moisture absorption (initial matric suction) of -1500 kPa is conservative, and the results are relatively safe.
- tailings stockyard /
- landfill /
- soil cover /
- vegetation soil /
- long-term service /
- hydraulic parameter /
- vegetation root system /
- design of seepage prevention

HTML全文

0. 引言

管棚超前预支护技术能在一定范围内较大幅度提高天然土体强度、刚度，有效减小地表沉降、拱顶沉降，确保隧道掘进支护安全，且因具有施工简便、经济、支护距离长和安全效益高的特点而成为经常采用的超前预支护方法。

目前国内外学者关于管棚预支护技术的作用机理和加固效果开展了大量研究。Hisatake等^[1]、Juneja等^[2]通过离心机模型试验，考虑了开挖效果的影响，研究了管棚支护的作用效果。郭璇等^[3]采用设计正交模型试验的方法，研究了管棚预支护的作用效果以及管棚对土压力分布的影响。周顺华等^[4]在原位观测和室内模型离心模型试验的基础上，提出了管棚作用的棚架原理。董新平^[5]则基于棚架体系建立了管棚的空间分析模型，研究了管棚支护作用的主要特征。Chungsik等^[6]、Ali等^[7]通过建立三维数值模型研究了管棚的受力和变形特征以及预支护效果。在理论分析方面，由于基于弹性地基梁理论的管棚分析模型具有简单实用的特点，故而具有很高的进一步研究完善的价值。李忠凯^[8]通过建立基于Winkler地基模型的管棚分析模型，推导出一个开挖循环过程中管棚内力和变形的计算表达式。郑俊杰等^[9]考虑了开挖对基床系数的影响，提出了变基床系数下管棚的地基梁模型。丁祖德等^[10]则对弹性地基梁模型中基床系数的取值展开详细研究。王炳智等^[11]通过改进弹性地基约束形式，采用Pasternak弹性地基理论建立管棚的分析模型并进行了求解。王道远等^[12]通过对比两个地基模型的计算结果与现场试验值的比较，得到基于Pasternak弹性地基理论的管棚模型更具合理性，并通过讨论地层反力系数以及管棚直径对管棚作用效果的影响，得出通过调整注浆改善地层反力系数或直径的方法来控制管棚变形效果甚微的结论。武松等^[13]根据弹性理论的刚度等效法，将管棚的一维弹性地基梁模型拓展至二维弹性地基各向异性板模型，并采用有限元法对模型进行了求解。

综上所述，针对基于弹性地基梁理论的管棚分析模型的研究已较为成熟，但在模型的建立过程中，对荷载的分布范围和管棚约束形式的认识上还存在矛盾，并且通过对模型的求解，并没有得到关于管棚设计和施工参数对管棚支护效果影响的一般性结论。预支护区段开挖支护完成后管棚全长范围内力和变形的最终分布并不明确。本文结合管棚超前预支护技术在实际运用过程中的问题以及国内外学者对管棚预支护技术研究的基础上，考虑掌子面失稳段对管棚约束的弱化以及隧道开挖支护动态施工过程，通过建立基于Pasternak弹性地基梁理论的管棚分析模型，采用理论推导的方法研究管棚预支护技术的支护效果以及受力特点，并系统地研究各设计施工参数对管棚力学响应的影响，从而为管棚的设计施工提供参考和借鉴。

1. 管棚分析模型的建立

根据管棚作用下隧道开挖过程的施工经验以及管棚既有设计方法，可采用Euler-Bernoulli梁来对管棚的受力和变形进行分析计算。根据管棚的受力特点，可将管棚沿纵向分为初期支护段、开挖段和未开挖段3部分来进行分析（见图 1）。位于洞口的导向墙由于具有独立的基础和仅承受管棚作用的特点，可认为导向墙对管棚的约束形式为固定端约束，且无初始位移。在隧道开挖过程中，开挖段在围岩压力的作用下将力传递给导向墙或初期支护以及掌子面前方围岩，管棚发生变形。将该开挖段支护后，管棚则以该变形状态固结在初期支护中，初期支护前端对管棚的约束可认为具有初始位移和初始转角的固定端约束，而此处管棚的变形值即为下个循环开挖分析时管棚起固定端处的初始变形值。在掌子面前方，可认为由于隧道的开挖，掌子面在地应力作用下发生失稳，则该段围岩将无法承受管棚变形引起的附加应力，则可认为掌子面失稳段与开挖段具有相同的性质。失稳段前方围岩承受因管棚变形而产生的附加力，该段管棚受围岩的弹性地基约束。

图 1 管棚预支护段施工示意图

Figure 1. Construction diagram of pre-support section of pipe roof

下载: 全尺寸图片幻灯片

在单个循环进尺的分析中，以初期支护前端作为分析模型的起点，为具有初始竖向位移 ${\omega _{\text{0}}}$ 和初始转角 ${\theta _{\text{0}}}$ 的固定端。管棚承受开挖段及掌子面失稳段上方的围岩压力。根据莫尔–库仑强度准则，可假设掌子面沿开挖方向失稳的破裂面与铅垂线的夹角为45° $- \varphi {\text{/2}}$ ，若单个循环开挖进尺长度为 $a$ ，则围岩压力作用范围为 $s = a + H\tan ({45^ \circ } - \varphi /2)$ ，其中 $H$ 为上台阶高度。由于一个循环开挖进尺长度较短，可认为作用于管棚上的围岩压力值在该分析段保持不变。而围岩压力在纵向方向上的连续变化可离散为每个开挖区段上的常围岩压力。破裂面前方围岩以弹性地基约束形式作用于管棚，从而建立单个循环进尺的管棚分析模型（见图 2）。在管棚超前支护作用下隧道的循环开挖，其力学行为等价于单个循环进尺的管棚分析模型的不断移动，通过迭代计算可得到管棚全长范围最终的内力和变形分布。

图 2 单个循环进尺的管棚分析模型

Figure 2. Analysis model for pipe roof with single cycle footage

下载: 全尺寸图片幻灯片

通过对管棚受力过程的分析，针对管棚的分析模型提出以下假设：①导向墙或初期支护对管棚的约束假设为固定端约束；②导向墙处管棚无初始挠度和转角；③假设掌子面前方失稳段围岩对管棚不具有约束作用；④掌子面前方岩体对管棚的约束反力采用Pasternak弹性地基理论来确定。

围岩对管棚的弹性地基约束根据管棚长度可分为两种类型。类型Ⅰ，埋入围岩中的管棚长度足够长，该条件下可认为管棚B端的边界条件为：当 $x \to \infty$ 时， $\omega (x) = 0$ ， $\theta (x) = 0$ ，分析图示如图 2（a）所示。类型Ⅱ，以类型Ⅰ为基本分析模型，随着在管棚预支护范围内开挖的不断推进，在对预支护段的最后一个循环进尺进行分析时，B端将发生变形。此时B端的边界条件为：当 $x = s + l$ 时， $M(x) = 0$ ， $Q(x) = 0$ ，其中 $l$ 为搭接长度，分析如图 2（b）所示。

管棚作为松散破碎地层中开挖隧道所采用的预支护技术，其预支护地层往往是破碎的散体介质，因此在进行管棚结构的分析计算时，采用松散围岩压力计算方法确定作用于管棚结构上的荷载比较符合实际情况。本模型浅埋隧道围岩压力的计算采用全土柱理论，深埋隧道围岩压力的计算采用普氏理论，则围岩压力可表示为统一的表达式：

。

${q^ * } = \gamma \cdot h。$

(1)

式中 $\gamma$ 为围岩重度（kN/m³）； $h$ 为上覆松动围岩的高度（m），浅埋隧道为取隧道埋深，深埋隧道取压力拱高度。

Pasternak弹性地基模型是在Winkler弹性地基模型的基础上假定弹簧之间存在剪切作用，引入地基剪切模量 ${G_{\text{p}}}$ ，则弹性地基作用于梁上的约束反力为

。

$p(x) = k{b^ * }\omega (x) - {G_{\text{p}}}{b^ * }\frac{{{{\text{d}}^2}\omega (x)}}{{{\text{d}}{x^2}}}。$

(2)

式中 $k$ 为地基基床系数（kN/m³）； $\omega (x)$ 为 $x$ 点处梁的挠度（m）； ${G_{\text{p}}}$ 为地基剪切模量（kN/m）； ${b^ * }$ 为地基剪切层的等效宽度（m），满足 ${b^ * } = b[1 + \sqrt {{G_{\text{p}}}{\text{/}}k} /b]$ ， $b$ 为管棚钢管的管径。Pasternak弹性地基模型相比Winkler地基模型通过引入地基剪切模量，弥补了Winkler模型地基不连续的缺陷，因此更符合实际情况，并且相比弹性半空间模型，Pasternak模型大大降低了计算复杂度。

2. 管棚分析模型的求解

分析Euler-Bernoulli梁在受外荷载作用发生变形时，在梁上任选微段作为分析单元，其中，剪力、弯矩及荷载的正方向均如图 3所示。

图 3 梁单元分析图示

Figure 3. Diagram of beam element analysis

下载: 全尺寸图片幻灯片

根据平衡条件可得

，

$\frac{{{\text{d}}Q(x)}}{{{\text{d}}x}} = - q(x)\text{，}$

(3)

。

$\frac{{{\text{d}}M(x)}}{{{\text{d}}x}} = Q(x)。$

(4)

再根据梁的挠曲线近似微分方程：

。

$EI\frac{{{\text{d}}{\omega ^2}(x)}}{{{\text{d}}{x^2}}} = - M(x)。$

(5)

将式（4），（5）代入式（3），即可得到分布荷载作用下，梁的变形控制微分方程为

。

$EI\frac{{{{\text{d}}^4}\omega (x)}}{{{\text{d}}{x^4}}} = q(x)。$

(6)

单个循环进尺的管棚分析模型如图 2所示，根据荷载和约束的分布形式，可分段建立管棚的控制微分方程并求解，得到OA段控制微分方程：

，

$EI\frac{{{{\text{d}}^4}\omega (x)}}{{{\text{d}}{x^4}}} = q\text{，}$

(7)

式中， $q = j{q^ * }$ ， $j$ 为管棚环向间距（m），即作用于单根梁上围岩压力的宽度。其通解为

。

${\omega _{\text{1}}}(x) = \frac{q}{{24EI}}{x^4} + {B_1}{x^3} + {B_2}{x^2} + {B_3}x + {B_4}。$

(8)

根据式（5）~（7）表示的梁的弯矩、剪力和荷载与梁挠度之间的微分关系，将挠度表达式进行求导运算，即可得到管棚转角、弯矩和剪力的分段表达式。其中OA段转角、弯矩和剪力的解析表达式分别为

，

${\theta _1}{\text{(}}x{\text{)}} = \frac{q}{{6EI}}{x^3} + 3{B_1}{x^2} + 2{B_2}x + {B_3}\text{，}$

(9)

${M_{\text{1}}}(x) = - EI\left( {\frac{q}{{2EI}}{x^2} + 6{B_1}x + 2{B_2}} \right) \text{，}$

(10)

${Q_1}(x) = - EI\left( {\frac{q}{{EI}}x + 6{B_1}} \right)。$

(11)

根据AB段的荷载分布形式，建立控制微分方程：

$EI\frac{{{{\text{d}}^4}\omega (x)}}{{{\text{d}}{x^4}}} - {G_{\text{p}}}{b^ * }\frac{{{{\text{d}}^2}\omega (x)}}{{{\text{d}}{x^2}}} + k{b^ * }\omega (x) = 0\text{，}$

(12)

令 ${\lambda ^{\text{4}}}{\text{ = }}\frac{{k{b^ * }}}{{4EI}}$ ， $\alpha = \sqrt {1 + \frac{{{G_{\text{p}}}{\lambda ^2}}}{k}}$ ， $\beta = \sqrt {1 - \frac{{{G_{\text{p}}}{\lambda ^2}}}{k}}$ ，求得其通解为

${\omega _2}(x) = {{\text{e}}^{ - \alpha \lambda x}}({B_5}\cos \beta \lambda x + {B_6}\sin \beta \lambda x) + \\ \;\;\;\;\;{{\text{e}}^{\alpha \lambda x}}({B_7}\cos \beta \lambda x + {B_8}\sin \beta \lambda x)。$

(13)

同理求得AB段转角、弯矩和剪力的解析表达式分别为

${\theta _2}(x) = - {B_5}{{\text{e}}^{ - \alpha \lambda x}}(\beta \lambda \sin \beta \lambda x + \alpha \lambda \cos \beta \lambda x) + \\ \;\;\;\;\;\;\;\;\;{B_6}{{\text{e}}^{ - \alpha \lambda x}}(\beta \lambda \cos \beta \lambda x - \alpha \lambda \sin \beta \lambda x) + \\ \;\;\;\;\;\;\;\;\;{B_7}{{\text{e}}^{\alpha \lambda x}}(\alpha \lambda \cos \beta \lambda x - \beta \lambda \sin \beta \lambda x) + \\ \;\;\;\;\;\;\;\;\; {B_8}{{\text{e}}^{\alpha \lambda x}}(\beta \lambda \cos \beta \lambda x + \alpha \lambda \sin \beta \lambda x)\text{，}$

(14)

${M_2}(x) = - EI\{ {B_5}{{\text{e}}^{ - \alpha \lambda x}}[{\text{(}}{\alpha ^2} - {\beta ^2}{\text{)}}{\lambda ^2}\cos \beta \lambda x + 2\alpha \beta {\lambda ^2}\sin \beta \lambda x] + \\ \;\;\; {B_6}{{\text{e}}^{ - \alpha \lambda x}}[({\alpha ^2} - {\beta ^2}){\lambda ^2}\sin \beta \lambda x - 2\alpha \beta {\lambda ^2}\cos \beta \lambda x] + \\ \;\;\; {B_7}{{\text{e}}^{\alpha \lambda x}}[({\alpha ^2} - {\beta ^2}){\lambda ^2}\cos \beta \lambda x - 2\alpha \beta {\lambda ^2}\sin \beta \lambda x] + \\ \;\;\; {B_8}{{\text{e}}^{\alpha \lambda x}}[({\alpha ^2} - {\beta ^2}){\lambda ^2}\sin \beta \lambda x + 2\alpha \beta {\lambda ^2}\cos \beta \lambda x]\} \text{，}$

(15)

${Q_2}(x) = - EI\{ {B_5}{{\text{e}}^{ - \alpha \lambda x}}[(3\alpha {\beta ^2} - {\alpha ^3}){\lambda ^3}\cos \beta \lambda x + ({\beta ^3} - 3{\alpha ^2}\beta ){\lambda ^3} \cdot \\ \sin \beta \lambda x] + {B_6}{{\text{e}}^{ - \alpha \lambda x}}[(3\alpha {\beta ^2} - {\alpha ^3}){\lambda ^3}\sin \beta \lambda x - ({\beta ^3} - 3{\alpha ^2}\beta ){\lambda ^3} \cdot \\ \cos \beta \lambda x] + {B_7}{{\text{e}}^{\alpha \lambda x}}[({\alpha ^3} - 3\alpha {\beta ^2}){\lambda ^3}\cos \beta \lambda x + ({\beta ^3} - 3{\alpha ^2}\beta ){\lambda ^3} \cdot \\ \sin \beta \lambda x] + {B_8}{{\text{e}}^{\alpha \lambda x}}[({\alpha ^3} - 3\alpha {\beta ^2}){\lambda ^3}\sin \beta \lambda x + (3{\alpha ^2}\beta - {\beta ^3}){\lambda ^3} \cdot \\ \;\;\;\;\;\; \cos \beta \lambda x]\} \text{，}$

(16)

式中， ${B_1}$ ， ${B_2}$ ， ${B_3}$ ， ${B_4}$ ， ${B_5}$ ， ${B_6}$ ， ${B_7}$ ， ${B_8}$ 均为待求积分常数。

2.1 Ⅰ类管棚分析模型的求解

类型Ⅰ的模型示意图如图 2（a）所示，对于AB段梁，由于B端满足 ${\omega _2}(x)\left| {_{x \to \infty } = 0} \right.$ ， ${\theta _2}(x)\left| {_{x \to \infty } = 0} \right.$ 的边界条件，则代入式（13）可求得 ${B_7}{\text{ = }}{B_{\text{8}}}{\text{ = 0}}$ 。由边界条件 ${\omega _1}(x)\left| {_{x = 0} = {\omega _0}} \right.$ ， ${\theta _{\text{1}}}(x)\left| {_{x = 0} = {\theta _0}} \right.$ ， ${\omega _1}(x)\left| {_{x = s} = {\omega _2}(x)} \right.\left| {_{x = s}} \right.$ ， ${\theta _1}(x)\left| {_{x = s} = {\theta _2}(x)} \right.\left| {_{x = s}} \right.$ ， ${M_1}(x)\left| {_{x = s} = {M_2}(x)} \right.\left| {_{x = s}} \right.$ ， ${Q_1}(x)\left| {_{x = s}} \right. =$ ${Q_2}(x)\left| {_{x = s}} \right.$ ，得方程组：

$\left[ {\begin{array}{*{20}{c}} {\text{0}}&{\text{0}}&{\text{0}}&{\text{1}}&{\text{0}}&{\text{0}} \\ {\text{0}}&{\text{0}}&{\text{1}}&{\text{0}}&{\text{0}}&{\text{0}} \\ {{s^3}}&{{s^2}}&s&1&{{\varphi _{35}}}&{{\varphi _{36}}} \\ {3{s^2}}&{2s}&1&0&{{\varphi _{45}}}&{{\varphi _{46}}} \\ {6s}&2&0&0&{{\varphi _{55}}}&{{\varphi _{56}}} \\ 6&0&0&0&{{\varphi _{65}}}&{{\varphi _{57}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{B_1}} \\ {{B_2}} \\ {{B_3}} \\ {{B_4}} \\ {{B_5}} \\ {{B_6}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{\omega _0}} \\ {{\theta _0}} \\ {{\psi _3}} \\ {{\psi _4}} \\ {{\psi _5}} \\ {{\psi _6}} \end{array}} \right]\text{，}$

(17)

式中， ${\varphi _{35}} = - {{\text{e}}^{ - \alpha \lambda s}}\cos \beta \lambda s$ ， ${\varphi _{3{\text{6}}}} = - {{\text{e}}^{ - \alpha \lambda s}}\sin \beta \lambda s$ ， ${\varphi _{{\text{4}}5}} =$ ${\text{e}}^{-\alpha \lambda s}(\beta \lambda \mathrm{sin}\beta \lambda s+\alpha \lambda \mathrm{cos}\beta \lambda s)，{\phi }_{\text{46}}=-{\text{e}}^{-\alpha \lambda s}(\beta \lambda \mathrm{cos}\beta \lambda s-$ $\alpha \lambda \mathrm{sin}\beta \lambda s)，{\phi }_{\text{5}5}=-{\text{e}}^{-\alpha \lambda s}[2\alpha \beta {\lambda }^{2}\mathrm{sin}\beta \lambda s+({\alpha }^{2}-{\beta }^{2}){\lambda }^{2}\cdot$ $\mathrm{cos}\beta \lambda s]，{\phi }_{\text{56}}=-{\text{e}}^{-\alpha \lambda s}[-2\alpha \beta {\lambda }^{2}\mathrm{cos}\beta \lambda s+({\alpha }^{2}-{\beta }^{2}){\lambda }^{2}\cdot$ $\mathrm{sin}\beta \lambda s]，{\phi }_{\text{6}5}=-{\text{e}}^{-\alpha \lambda s}[({\beta }^{3}-3{\alpha }^{2}\beta ){\lambda }^{3}\mathrm{sin}\beta \lambda s-({\alpha }^{3}-3\alpha {\beta }^{2})\cdot$ ${\lambda }^{3}\mathrm{cos}\beta \lambda s]，{\phi }_{\text{66}}={\text{e}}^{-\alpha \lambda s}[({\beta }^{3}-3{\alpha }^{2}\beta ){\lambda }^{3}\mathrm{cos}\beta \lambda s+({\alpha }^{3}-$ $3\alpha {\beta ^2}){\lambda ^3}\sin \beta \lambda s]$ ， ${\psi _3} = - q{s^4}/(24EI)$ ， ${\psi }_{\text{4}}=-q{s}^{\text{3}}/(\text{6}EI)，$ ${\psi _{\text{5}}} = - q{s^{\text{2}}}/(2EI)$ ， ${\psi _{\text{6}}} = - qs/(EI)$ 。

通过式（17）即可求得积分常数 ${B_1}$ ， ${B_2}$ ， ${B_3}$ ， ${B_4}$ ， ${B_5}$ ， ${B_6}$ ，将其回代到挠度、转角、弯矩和剪力表达式中，即可求得满足类型Ⅰ分析模型边界条件下管棚的变形和受力分布。

2.2 Ⅱ类管棚分析模型的求解

类型Ⅱ的模型示意图如所示，由边界条件 ${\omega _1}(x)\left| {_{x = 0} = {\omega _0}} \right.$ ， ${\theta _{\text{1}}}(x)\left| {_{x = 0} = {\theta _0}} \right.$ ， ${\omega _1}(x)\left| {_{x = s} = {\omega _2}(x)} \right.\left| {_{x = s}} \right.$ ， ${\theta _1}(x)\left| {_{x = s} = {\theta _2}(x)} \right.\left| {_{x = s}} \right.$ ， ${M_1}(x)\left| {_{x = s} = {M_2}(x)} \right.\left| {_{x = s}} \right.$ ， ${Q_1}(x)\left| {_{x = s} = } \right.$ ${Q_2}(x)\left| {_{x = s}} \right.$ ， ${M_{\text{2}}}(x)\left| {_{x = s + l} = } \right.0$ ， ${Q_2}(x)\left| {_{x = s + l}} \right. = 0$ ，得方程组：

$\left[ {\begin{array}{*{20}{c}} {\text{0}}&{\text{0}}&{\text{0}}&{\text{1}}&{\text{0}}&{\text{0}}&{\text{0}}&{\text{0}} \\ {\text{0}}&{\text{0}}&{\text{1}}&{\text{0}}&{\text{0}}&{\text{0}}&{\text{0}}&{\text{0}} \\ {{s^3}}&{{s^2}}&s&1&{{\varphi _{35}}}&{{\varphi _{36}}}&{{\varphi _{37}}}&{{\varphi _{38}}} \\ {3{s^2}}&{2s}&1&0&{{\varphi _{45}}}&{{\varphi _{46}}}&{{\varphi _{47}}}&{{\varphi _{48}}} \\ {6s}&2&0&0&{{\varphi _{55}}}&{{\varphi _{56}}}&{{\varphi _{57}}}&{{\varphi _{58}}} \\ 6&0&0&0&{{\varphi _{65}}}&{{\varphi _{66}}}&{{\varphi _{67}}}&{{\varphi _{68}}} \\ 0&0&0&0&{{\varphi _{75}}}&{{\varphi _{76}}}&{{\varphi _{77}}}&{{\varphi _{78}}} \\ 0&0&0&0&{{\varphi _{85}}}&{{\varphi _{86}}}&{{\varphi _{87}}}&{{\varphi _{88}}} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{B_1}} \\ {{B_2}} \\ {{B_3}} \\ {{B_4}} \\ {{B_5}} \\ {{B_6}} \\ {{B_7}} \\ {{B_8}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{\omega _0}} \\ {{\theta _0}} \\ {{\psi _3}} \\ {{\psi _4}} \\ {{\psi _5}} \\ {{\psi _6}} \\ 0 \\ 0 \end{array}} \right]\text{，}$

(18)

式中， ${\varphi _{35}} = - {{\text{e}}^{ - \alpha \lambda s}}\cos \beta \lambda s$ ， ${\phi }_{3\text{6}}=-{\text{e}}^{-\alpha \lambda s}\mathrm{sin}\beta \lambda s，{\phi }_{37}=$ $-{\mathrm{e}}^{\alpha \lambda s}\mathrm{cos}\beta \lambda s，$ ${\phi }_{3\text{8}}=-{e}^{\alpha \lambda s}\mathrm{sin}\beta \lambda s，{\phi }_{\text{4}5}={\text{e}}^{-\alpha \lambda s}(\beta \lambda \mathrm{sin}\beta \cdot$ $\lambda s + \alpha \lambda \cos \beta \lambda s)$ ， ${\varphi _{{\text{46}}}} = - {{\text{e}}^{ - \alpha \lambda s}}(\beta \lambda \cos \beta \lambda s - \alpha \lambda \sin \beta \lambda s)$ ， ${\varphi _{{\text{47}}}} = - {{\text{e}}^{\alpha \lambda s}}(\alpha \lambda \cos \beta \lambda s - \beta \lambda \sin \beta \lambda s)$ ， ${\varphi _{{\text{48}}}} = - {{\text{e}}^{\alpha \lambda s}} \cdot (\beta \lambda$ $\mathrm{cos}\beta \lambda s\text{+}\alpha \lambda \mathrm{sin}\beta \lambda s)，{\phi }_{\text{5}5}=-{\text{e}}^{-\alpha \lambda s}[2\alpha \beta {\lambda }^{2}\mathrm{sin}\beta \lambda s+({\alpha }^{2}-$ ${\beta }^{2}){\lambda }^{2}\mathrm{cos}\beta \lambda s]，{\phi }_{\text{56}}=-{\text{e}}^{-\alpha \lambda s}[-2\alpha \beta {\lambda }^{2}\mathrm{cos}\beta \lambda s+({\alpha }^{2}-{\beta }^{2})$ ${\lambda ^2}\sin \beta \lambda s]$ ， ${\varphi _{{\text{57}}}} = - {{\text{e}}^{\alpha \lambda s}}[ - 2\alpha \beta {\lambda ^2}\sin \beta \lambda s + ({\alpha ^2} - {\beta ^2}){\lambda ^2}$ $\mathrm{cos}\beta \lambda s]，$ ${\varphi _{{\text{58}}}} = - {{\text{e}}^{\alpha \lambda s}}[2\alpha \beta {\lambda ^2}\cos \beta \lambda s + ({\alpha ^2} - {\beta ^2}){\lambda ^2} \cdot \sin \beta$ $\lambda s]，{\phi }_{\text{6}5}=-{\text{e}}^{-\alpha \lambda s}[({\beta }^{3}-3{\alpha }^{2}\beta ){\lambda }^{3}\mathrm{sin}\beta \lambda s-({\alpha }^{3}-3\alpha {\beta }^{2})\cdot {\lambda }^{3}$ $\mathrm{cos}\beta \lambda s]，{\phi }_{\text{66}}={\text{e}}^{-\alpha \lambda s}[({\beta }^{3}-3{\alpha }^{2}\beta ){\lambda }^{3}\mathrm{cos}\beta \lambda s+({\alpha }^{3}-3\alpha {\beta }^{2})\cdot {\lambda }^{3}$ $\mathrm{sin}\beta \lambda s]，{\phi }_{\text{67}}=-{\text{e}}^{\alpha \lambda s}[({\alpha }^{3}-3\alpha {\beta }^{2}){\lambda }^{3}\mathrm{cos}\beta \lambda s\text{+}({\beta }^{3}-3{\alpha }^{2}\beta )\cdot$ ${\lambda }^{3}\mathrm{sin}\beta \lambda s]，{\phi }_{\text{68}}=-{\text{e}}^{\alpha \lambda s}[({\alpha }^{3}-3\alpha {\beta }^{2}){\lambda }^{3}\mathrm{sin}\beta \lambda s-({\beta }^{3}-3{\alpha }^{2}\beta )\cdot$ ${\lambda ^3}\cos \beta \lambda s]$ ， ${\varphi _{{\text{7}}5}} = {{\text{e}}^{ - \alpha \lambda (s + l)}}[2\alpha \beta {\lambda ^2}\sin \beta \lambda (s + l) + ({\alpha ^2} -$ ${\beta ^2})$ ${\lambda ^2}\cos \beta \lambda (s + l)]$ ， ${\varphi _{{\text{76}}}} = {{\text{e}}^{ - \alpha \lambda (s + l)}}[ - 2\alpha \beta {\lambda ^2}\cos \beta \lambda (s +$ $l)+({\alpha }^{2}-{\beta }^{2}){\lambda }^{2}\mathrm{sin}\beta \lambda (s+l)]，{\phi }_{\text{77}}={\text{e}}^{\alpha \lambda (s+l)}[-2\alpha \beta {\lambda }^{2}\cdot$ $\sin \beta \lambda (s + l) + ({\alpha ^2} - {\beta ^2}){\lambda ^2}\cos \beta \lambda (s + l)]$ ， ${\varphi _{{\text{78}}}} = {{\text{e}}^{\alpha \lambda (s + l)}} \cdot$ $[2\alpha \beta {\lambda }^{2}\mathrm{cos}\beta \lambda (s+l)+({\alpha }^{2}-{\beta }^{2}){\lambda }^{2}\mathrm{sin}\beta \lambda (s+l)]，{\phi }_{\text{8}5}={\text{e}}^{-\alpha \lambda (s+l)}\cdot$ $[({\beta ^3} - 3{\alpha ^2}\beta ){\lambda ^3}\sin \beta \lambda (s + l) - ({\alpha ^3} - 3\alpha {\beta ^2}) \cdot {\lambda ^3}\cos \beta \lambda (s +$ $l)]，{\phi }_{\text{86}}={\text{e}}^{-\alpha \lambda (s+l)}[(3{\alpha }^{2}\beta -{\beta }^{3}){\lambda }^{3}\mathrm{cos}\beta \lambda (s+l)-({\alpha }^{3}-3\alpha {\beta }^{2})\cdot$ ${\lambda }^{3}\mathrm{sin}\beta \lambda (s+l)]，{\phi }_{\text{87}}={\text{e}}^{\alpha \lambda (s+l)}[({\beta }^{3}-3{\alpha }^{2}\beta )\cdot {\lambda }^{3}\mathrm{sin}\beta \lambda (s+$ $l)+({\alpha }^{3}-3\alpha {\beta }^{2}){\lambda }^{3}\mathrm{cos}\beta \lambda (s+l)]，{\phi }_{\text{88}}={\text{e}}^{\alpha \lambda (s+l)}\cdot [(3{\alpha }^{2}\beta -$ ${\beta }^{3}){\lambda }^{3}\mathrm{cos}\beta \lambda (s+l)+({\alpha }^{3}-3\alpha {\beta }^{2}){\lambda }^{3}\mathrm{sin}\beta \lambda (s+l)]，{\psi }_{3}=$ $- q{s^4}/(24EI)$ ， ${\psi _{\text{4}}} = - q{s^{\text{3}}}/({\text{6}}EI)$ ， ${\psi _{\text{5}}} = - q{s^{\text{2}}}/(2EI)$ ， ${\psi _{\text{6}}} =$ $- qs{\text{/}}(EI)$ 。

同理，通过式（18）即可求得积分常数 ${B_1}$ ， ${B_2}$ ， ${B_3}$ ， ${B_4}$ ， ${B_5}$ ， ${B_6}$ ， ${B_7}$ ， ${B_8}$ ，将其回代到挠度、转角、弯矩和剪力表达式中，即可求得满足类型Ⅱ分析模型边界条件下管棚的变形和受力分布。

3. 案例计算

南宁至崇左铁路引入机场隧道起讫里程DK25+615—DK28+935，全长3320 m。其中DK25+940—DK27+234段长度1294 m为暗挖区间，埋深3.4~15.2 m，施工工法为CRD法（350 m）和三台阶临时仰拱（944 m）；该隧道暗挖区间存在浅埋、开挖下穿既有道路、临近既有构（建）筑物。隧址内岩层上部为杂填土，下部为粉质黏土、泥质砂岩夹砂质泥岩、砾岩，弱风化至全风化，岩溶发育，层面起伏变化大。本隧道进口明暗分界处拱部设置一环35 m长Φ108 mm×6 mm管棚加强支护，管棚环向间距40 cm，具体参数见表 1。

表 1 设计参数取值

Table 1. Values of design parameters

基床系数/(kN·m^-3)	地基剪切模量/(kN·m^-1)	隧道埋深/m	内摩擦角/(°)	围岩重度/(kN·m^-3)	等效惯性矩/m⁴	等效弹性模量/kPa	开挖高度/m	管棚间距/m	进尺深度/m
30000	2800	6	45.6	24	6.68×10^-6	7.89×10⁷	2.4	0.4	0.6

下载: 导出CSV

| 显示表格

本文以隧道进口明暗分界处拱顶单根钢管为研究对象，则该钢管的挠曲变形与隧道因开挖引起的拱顶沉降值相同。对管棚分析模型中的参数分别确定取值，并利用MATLAB编程进行应用分析。第一个循环进尺下管棚挠度、转角、弯矩和剪力在开挖方向上的分布曲线如图 4所示。

图 4 管棚的变形及内力分布

Figure 4. Distribution of deformation and internal force of pipe roof

下载: 全尺寸图片幻灯片

图 4表明，管棚挠度在OA段逐渐增大，且在A点处达到最大值7.2 mm。AB段挠度逐渐减小，当计算深度为4 m时，该循环开挖对该处管棚的变形无影响。随计算深度的加深，弯矩呈震荡形式趋近于零，其中最大值位于固定端处为-23.5 kN·m。剪力的分布同样随计算深度的增大而收敛于零，最大值同样位于固定端处，大小为62.6 kN。当计算深度达到4 m时，管棚的弯矩和剪力将趋近于零，由此可得一个循环开挖区段对管棚的影响范围长度为4 m。

通过以上单个开挖区段的分析，得到管棚在该分析区段的变形和受力特征。该区段前0.6 m的变形与受力将随初期支护的完成而与初期支护固结在一起，0.6 m处的变形值将成为下一个分析区段的位移初值。通过迭代计算，可得到管棚预支护范围内30 m全部开挖支护完成后隧道拱顶下沉分布如图 5所示，管棚全长范围内内力的分布情况如图 6所示。

图 5 管棚预支护范围内拱顶下沉纵向分布

Figure 5. Longitudinal distribution of vault settlement within scope of pre-support of pipe roof

下载: 全尺寸图片幻灯片

图 6 管棚全长范围的内力分布

Figure 6. Distribution of internal force in whole length of pipe roof

下载: 全尺寸图片幻灯片

图 5表明，管棚的挠曲变形随掌子面的不断推进而逐渐增大。当拱顶沉降达到22.8 mm时，随掌子面的进一步推进，拱顶沉降将基本保持不变。

王道远等^[12]通过现场试验得出，管棚的变形整体呈凹形分布，在开挖初期，管棚的变形速率较快。而图中管棚变形分布特征符合这一规律，且管棚变形的数量级基本一致，从而证明了该分析模型的合理性和实用性。

图 6表明，管棚的弯矩和剪力的最不利截面均发生在导向墙的约束处，其中弯矩最大值为23.5 kN·m，剪力最大值为62.6 kN。管棚全长范围内弯矩基本保持在12 kN·m，剪力基本保持在52 kN。

综上所述，利用该分析模型可以系统地对管棚的支护效果和受力特征进行预测。通过对算例的数值计算，得到管棚预结构的受力和变形分布特征，从而体现了管棚预支护体系的主要特征。

4. 参数分析

本节在前述理论研究的基础上，通过调整管棚的环向间距、管径、开挖进尺长度以及搭接长度的设计参数，研究各设计参数对管棚受力和变形的影响形式和程度，从而指导管棚的设计和施工。

4.1 环向间距对管棚受力和变形的影响分析

取其他参数不变，当管棚环向间距j为0.40，0.60，0.80，1.00 m时，分别计算并绘出管棚的变形分布曲线（见图 7）和内力分布曲线（见图 8），研究环向间距对管棚变形和受力的影响。

图 7 环向间距对管棚变形的影响

Figure 7. Influences of circumferential spacing on deformation of pipe roof

下载: 全尺寸图片幻灯片

图 8 环向间距对管棚受力的影响

Figure 8. Influence of circumferential spacing on stress of pipe roof

下载: 全尺寸图片幻灯片

表明，管棚环向间距的变化不改变其变形分布形式，最大挠度 ${\omega _{\max }}{\text{(mm)}}$ 和最大转角 ${\theta _{\max }}$ (°)发生的位置不变，其中最大挠度和最大转角均与环向间距之间满足线性增长关系，关系式分别为

$\left. \begin{array}{l}{\omega }_{\mathrm{max}}=18.159j\text{ }\text{，}\\ {\theta }_{\mathrm{max}}=1.3\text{8}63j\text{ }。\end{array} \right\}$

(19)

表明，管棚环向间距同样不改变管棚的内力分布形式，最不利截面不随环向间距的变化而变化，管棚的最大弯矩 ${M_{\max }}{\text{(kN}} \cdot {\text{m)}}$ 和最大剪力 ${Q_{\max }}{\text{(kN)}}$ 均与环向间距之间满足线性增长关系，关系式分别为

$\left. \begin{array}{l}{M}_{\mathrm{max}}=58.657j\text{ }\text{，}\\ {Q}_{\mathrm{max}}=156.61j\text{ }。\end{array} \right\}$

(20)

4.2 管径对管棚受力和变形的影响分析

取其他参数不变，当管径 $b$ 分别为0.089，0.108，0.152，0.180 m时，分别计算并绘制出管棚的挠度分布曲线（见图 9）和弯矩分布曲线（见图 10），研究管径对管棚变形和受力的影响。当管棚的管径发生变化时，管棚的等效截面惯性矩和等效弹性模量均发生变化，根据刚度等效法^[14]求得不同管径下的等效截面惯性矩和等效弹性模量见表 2。

图 9 管径对管棚挠度的影响

Figure 9. Influence of pipe diameter on deflection of pipe roof

下载: 全尺寸图片幻灯片

图 10 管径对管棚弯矩的影响

Figure 10. Influences of pipe diameter on bending moment of pipe roof

下载: 全尺寸图片幻灯片

表 2 不同管径的截面计算参数

Table 2. Cross section calculation parameters of different pipe diameters

管径/m	等效惯性矩/(10^-6 m⁴)	等效弹性模量/(10⁷ kPa)
0.089	3.08	9.24
0.108	6.68	7.89
0.152	26.00	5.89
0.180	51.50	5.07

下载: 导出CSV

| 显示表格

从可以看出，管径越粗，管棚的挠曲变形越平缓，随着管径的增大，最大挠度 ${\omega _{\max }}$ 逐渐减小，并且随着管径的不断增大对管棚挠曲变形的影响将越来越小。最大挠度 ${\omega _{\max }}$ 与管径b(m)之间满足

${\omega _{\max }} = 0.0397{b^{ - 2.346}}。$

(21)

从可以看出，随着管径的增大，弯矩的危险截面位置不变，管棚所受的最大弯矩 ${M_{\max }}$ （kN·m）也随着增大。最大弯矩与管径之间满足

${M_{\max }} = - 115.67{b^2} + 100.87b + 13.917。$

(22)

4.3 开挖进尺长度对管棚受力和变形的影响分析

取其他参数不变，当开挖进尺长度 $a$ 分别为0.5，0.75，1，1.25，1.5 m时，分别绘制出管棚变形分布曲线（图 11）和内力分布曲线（图 12），研究开挖进尺长度对管棚变形和受力的影响。

图 11 开挖进尺长度对管棚变形的影响

Figure 11. Influences of excavation footage length on deformation of pipe roof

下载: 全尺寸图片幻灯片

图 12 开挖进尺长度对管棚受力的影响

Figure 12. Influences of excavation footage length on stress of pipe roof

下载: 全尺寸图片幻灯片

表明，随着开挖进尺长度的增大，管棚最大变形位置不断前移，且最大值不断增大。通过开挖进尺长度与最大变形关系曲线图可得到，开挖进尺长度与最大挠度 ${\omega _{\max }}$ 和最大转角 ${\theta _{\max }}$ 之间呈非线性增长关系，进尺越长，变形增长速度越快，关系式分别为

$\left. \begin{array}{l}{\omega }_{\mathrm{max}}=9.6241{a}^{2}+1.9131a+2.6503\text{ }\text{，}\\ {\theta }_{\mathrm{max}}=0.2622{a}^{2}+0.487a+0.1679\text{ }。\end{array} \right\}$

(23)

说明开挖进尺长度是控制管棚变形的关键施工技术指标。

表明，开挖进尺长度的增大同样会使管棚的内力分布前移，但最不利截面总位于固定端处。最大弯矩 ${M_{\max }}$ 与开挖进尺长度之间满足非线性关系式：

${M_{\max }} = 4.5114{a^2} + 17.093a + 11.583\text{，}$

(24)

式中，开挖进尺长度 $a$ （m）， ${M_{\max }}$ 最大弯矩（kN·m）。而最大剪力 ${Q_{\max }}$ （kN）与开挖进尺长度之间满足线性关系式：

${Q_{\max }} = 30.841x + 44.13。$

(25)

4.4 搭接长度对管棚受力和变形的影响分析

由类型Ⅱ分析模型可计算管棚搭接长度对管棚作用效果的影响，搭接长度 $D = H\tan ({45^ \circ } - \varphi /$ $2) + l$ ，分别取搭接长度为1.1，1.3，1.5，5 m，分别绘制出管棚挠度分布曲线以及搭接长度与最大挠度关系曲线（图 13）和弯矩分布曲线以及搭接长度与最大弯矩关系曲线（图 14），研究搭接长度对管棚变形和受力的影响。

图 13 搭接长度对管棚挠度的影响

Figure 13. Influences of lap length on deflection of pipe roof

下载: 全尺寸图片幻灯片

图 14 搭接长度对管棚弯矩的影响

Figure 14. Influences of lap length on bending moment of pipe roof

下载: 全尺寸图片幻灯片

图 13表明，随着搭接长度的增大，管棚发生最大挠度的位置和挠度分布形式基本不变。当搭接长度大于2.5 m时，最大挠度将趋近于7.3 mm而保持不变。

图 14表明，随着搭接长度的增大，管棚的最不利截面位置和弯矩的分布形式基本不变。当搭接长度大于2.5 m时，最大弯矩将趋近于23.5 kN·m而保持不变。综上所述，在该工程条件下，当管棚搭接长度大于2.5 m时，搭接长度的变化将对管棚的变形和受力基本没有影响。

5. 结论

（1）通过建立和求解基于Pasternak弹性地基理论的单个循环进尺下管棚分析模型可得管棚的受力和变形的分布特征。考虑隧道开挖支护动态施工过程，通过迭代计算可有效预测管棚预支护段开挖支护完成后管棚整体的受力和变形分布，且其理论值与已有文献实测结果具有相同的分布特征，从而证明了模型的合理性。

（2）管棚环向间距的变化不改变管棚内力和变形的分布形式，且最不利截面位置保持不变，管棚的受力和变形与环向间距之间呈线性增长关系，随环向间距的增大而增大。

（3）选用的管径越粗，管棚的挠曲变形越平缓，挠曲变形越小，并且随着管径的不断增大，管径对管棚挠曲变形的影响将越来越小。随着管径的增大，弯矩的危险截面位置不变，管棚的最大弯矩与管径之间呈非线性增长关系。

（4）随着开挖进尺长度的增大，管棚最大变形位置不断前移，最大变形值与开挖进尺长度之间呈非线性增长关系，开挖进尺越长，变形增长速度越快。说明开挖进尺长度是控制管棚变形的关键施工技术指标。开挖进尺长度越大，最大弯矩和最大剪力越大。

（5）管棚的最大内力及变形随搭接长度的增长而迅速趋于收敛。在本文案例工况下，当管棚搭接长度大于2.5 m时，搭接长度的变化将对管棚的变形和受力基本没有影响。

图 1 土质封场覆盖层水分存储-释放循环防渗机理物理模型

Figure 1. Model for seepage prevention mechanism of water storage-release cycle in soil cover

下载: 全尺寸图片幻灯片

图 2 现场^#1、^#2试验区张力计和TDR埋设

Figure 2. Embedment of tensiometer and TDR in test areas No. 1 and 2

下载: 全尺寸图片幻灯片

图 3 现场双套环原位入渗试验简图

Figure 3. Principle of field double-ring in infiltration tests

下载: 全尺寸图片幻灯片

图 4 现场覆盖层红黏土原状土取样

Figure 4. Sampling of field undisturbed red clay soil

下载: 全尺寸图片幻灯片

图 5 现场无植被2015—2020年服役期红黏土入渗曲线

Figure 5. Infiltration curves of red clay without vegetation from 2015 to 2020

下载: 全尺寸图片幻灯片

图 6 现场有植被2015—2020年服役期红黏土入渗曲线

Figure 6. Infiltration curves of red clay with vegetation from 2015 to 2020

下载: 全尺寸图片幻灯片

图 7 2015—2020年各年有植被与无植被条件渗透系数对比

Figure 7. Comparison of permeability coefficients with and without vegetation from 2015 to 2020

下载: 全尺寸图片幻灯片

图 8 自然气候下2015—2020年服役期红黏土土水特征曲线（脱湿过程）

Figure 8. Soil-water characteristic curve of red clay in natural climate from 2015 to 2020 (moisture desorption)

下载: 全尺寸图片幻灯片

图 9 自然气候下2015—2020年服役期红黏土土水特征曲线（吸湿过程）

Figure 9. Soil-water characteristic curves of red clay in natural climate from 2016 to 2020 (moisture absorption process)

下载: 全尺寸图片幻灯片

图 10 2015—2020年间各年土水特征曲线变化分布范围图

Figure 10. Envelope map of distribution range of soil-water characteristic curves from 2015 to 2020

下载: 全尺寸图片幻灯片

图 11 体积压力板仪测试吸湿土水特征曲线与大孔隙“瓶颈（或墨水瓶）效应”^[16]

Figure 11. Soil-water characteristic curve of hygroscopic soil measured by volume pressure plate apparatus and "bottleneck effect" of macropores

下载: 全尺寸图片幻灯片

图 12 长期服役原状红黏土有效储水率θ_a

Figure 12. Effective water storage rate (θ_a) of undisturbed red clay in long-term service

下载: 全尺寸图片幻灯片

图 13 现场自然降雨气候下红黏土土水特征点（第一场雨）

Figure 13. Soil-water characteristic points of red clay in natural rainfall climate (the 1st rain)

下载: 全尺寸图片幻灯片

图 14 现场自然降雨气候下红黏土土水特征点（第二场雨）

Figure 14. Soil-water characteristic points of red clay under natural rainfall climate (the 2nd rain)

下载: 全尺寸图片幻灯片

图 15 土水特征曲线滞回圈内各曲线定义和示意图^[17]

Figure 15. Definition and schematic diagram of each curve in hysteresis loop of soil-water characteristic curves

下载: 全尺寸图片幻灯片

图 16 现场红黏土在两场自然降雨气候条件下的土水特征响应

Figure 16. Soil-water characteristic responses of red clay in two natural rainfalls

下载: 全尺寸图片幻灯片

表 1 现场长期监测试验工况设计

Table 1 Design of field long-term monitoring tests

测试区	位置	植被条件	有/无仪器埋设	埋设深度
^#1	4级平台	无	张力计、水分TDR	15，45 cm（各2支）
^#2	5级平台	有	张力计、水分TDR	15，45 cm（各2支）
^#3	6级平台	有、无（约一半）	无	无

下载: 导出CSV

表 2 现场和室内（实验室）开展的测试项目与测试方法

Table 2 Items and methods of field and indoor (laboratory) tests

现场	仪器和方法	室内	仪器和方法
现场覆盖层红黏土体积含水率和基质吸力监测	TDR/张力计联合监测	重塑红黏土基本物理性质和参数	室内土工常规测试
现场覆盖层红黏土双套环原位入渗	下文详述	重塑土室内变水头饱和渗透系数	室内变水头渗透试验
现场降雨等气候条件	翻斗式雨量计和气象信息网	室内重塑土的土水特征曲线	体积压力板仪（0~-1500 kPa）
现场覆盖层红黏土裂缝观测	详见文献[8]	现场取回的覆盖层原状土土水特征曲线	体积压力板仪（0~-1500 kPa）

下载: 导出CSV

表 3 红黏土基本物理力学参数

Table 3 Basic physical and mechanical parameters of red clay

天然含水率/%	天然密度/(g·cm^-3)	相对质量密度	液限/%	塑限/%	最优含水率/%	最大干密度/(g·cm^-3)	压缩系数a_1-2/MPa^-1	内摩擦角/(°)	黏聚力/kPa	孔隙比
31.5	1.67	2.71	70.65	40.32	23.8	1.61	0.21	14.5	18.7	1.13

下载: 导出CSV

表 4 自然气候下2015—2020年服役期红黏土土水特征曲线V-G模型拟合参数（脱湿过程）

Table 4 Fitting parameters of V-G model of soil-water characteristic curves of red clay in natural climate from 2015 to 2020 (moisture desorption process)

服役时间	饱和含水率θ_s/%	残余含水率θ_r/%	与进气值有关的倒数a	脱水速率n
建造之初(2015年)	49.05	18.78	2.59	1.01
第1年 (2016年)	50.27	20.12	6.31	1.16
第2年 (2017年)	50.31	19.74	8.52	1.19
第3年 (2018年)	52.41	22.09	13.33	1.58
第4年 (2019年)	53.95	19.47	10.86	1.55
第5年 (2020年)	51.67	20.97	11.07	1.60

下载: 导出CSV

表 5 自然气候下2015—2020年服役期红黏土土水特征曲线V-G模型拟合参数（吸湿过程）

Table 5 Fitting parameters of V-G model of soil-water characteristic curves of red clay in natural climate from 2015 to 2020 (moisture absorption process)

服役时间	饱和含水率θ_s/%	残余含水率θ_r/%	与进气值有关的倒数a	脱水速率n
建造之初(2015年)	34.10	18.78	43.30	1.55
第1年 (2016年)	37.10	20.12	27.16	1.51
第2年 (2017年)	36.08	19.74	30.63	1.54
第3年 (2018年)	32.94	19.74	32.15	1.75
第4年 (2019年)	36.10	22.09	42.47	2.11
第5年 (2020年)	32.10	20.97	28.98	1.90

下载: 导出CSV

表 6 长期服役原状红黏土有效储水率θ_a

Table 6 Effective water storage rate (θ_a) of undisturbed red clay in long-term service

水力条件	时间/a	田间持水率θ_c/%	枯萎持水率θ_r/%	有效储水率θ_a/%
脱湿	0	36.08	18.78	17.30
	1.07	37.15	20.12	17.03
	2.03	37.85	19.74	18.11
	3.0	39.05	22.09	16.96
	4.1	38.35	19.47	18.88
	5.0	37.65	20.97	16.68
吸湿	0	31.00	18.78	12.22
	1.07	31.75	20.12	11.63
	2.03	37.15	19.74	17.41
	3.0	32.94	22.09	10.85
	4.1	29.06	19.47	9.59
	5.0	28.43	20.97	7.46
*考虑测试过程中天平称量最小精度、体积压力板仪滴定管肉眼读数最小精度和测试加压压力表精度等，有效储水率θ_a的误差率为1.2%。

下载: 导出CSV

参考文献(17)

[1]	CHEN C. Meteorological Conditions for Design of Monolithic Alternative Earthen Covers(AEFCs)[D]. Madison: University of Wisconsin, 1999.
[2]	BENSON C H, ALBRIGHT W H, ROESLER A C, et al. Evaluation of final cover performance: Field data from the Alternative Cover Assessment Program (ACAP)[C]// Proceedings of Waste Management. 2002.
[3]	ZHANG Wen-jie, GENG Xiao, Performance and mechanism of capillary-barrier evaportranspiration cover of landfills[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(3): 454–459. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202102002.htm
[4]	张文杰, 林午, 董林兵. 垃圾填埋场毛细阻滞型腾发封顶模型试验[J]. 岩土力学, 2014, 35(5): 1263–1268. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201405007.htm ZHANG Wen-jie, LIN Wu, DONG Lin-bing. Model test of a capillary barrier evapotranspiration cover for landfills[J]. Rock and Soil Mechanics, 2014, 35(5): 1263–1268. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201405007.htm
[5]	ALBRIGHT W H, BENSON C H, GEE G W, et al. Field water balance of landfill final covers[J]. Journal of Environmental Quality, 2004, 33(6): 2317–2332. doi: 10.2134/jeq2004.2317
[6]	ALBRIGHT W H, BENSON C H, GEE G W, et al. Field performance of a compacted clay landfill final cover at a humid site[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2006, 132(11): 1393–1403. doi: 10.1061/(ASCE)1090-0241(2006)132:11(1393)
[7]	NI J J, LEUNG A K, NG C W W, et al. Investigation of plant growth and transpiration-induced matric suction under mixed grass–tree conditions[J]. Canadian Geotechnical Journal, 2017, 54(4): 561–573. doi: 10.1139/cgj-2016-0226
[8]	BENSON C H, SAWANGSURIYA A, TRZEBIATOWSKI B, et al. Postconstruction changes in the hydraulic properties of water balance cover soils[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2007, 133(4): 349–359. doi: 10.1061/(ASCE)1090-0241(2007)133:4(349)
[9]	VAN GENUCHTEN M T. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils[J]. Soil Science Society of America Journal, 1980, 44(5): 892–898. doi: 10.2136/sssaj1980.03615995004400050002x
[10]	焦卫国, 季永新, 张玥, 等. 红黏土土质覆盖层水力参数随服役时间演变规律分析[J]. 岩土工程学报, 2021, 43(6): 1059–1068. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202106013.htm JIAO Wei-guo, JI Yong-xin, ZHANG Yue, et al. Evolution of hydraulic parameters of red clay cover with service time [J]. Chinese Journal of Geotechnical Engineering, 2021, 43(6): 1059–1068. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202106013.htm
[11]	谭邦宏. 红黏土作为垃圾填埋场防渗垫层的研究进展[J]. 四川建筑, 2015, 35(6): 107–108, 111. doi: 10.3969/j.issn.1007-8983.2015.06.039 TAN Bang-hong. Research progress of red clay as impermeable cushion in landfill[J]. Sichuan Architecture, 2015, 35(6): 107–108, 111. (in Chinese) doi: 10.3969/j.issn.1007-8983.2015.06.039
[12]	彭玉林, 龚爱民, 孙海燕, 等. 垃圾填埋场中改性红黏土防渗料的性能研究[J]. 人民长江, 2011, 42(增刊2): 163–165, 169. https://www.cnki.com.cn/Article/CJFDTOTAL-RIVE2011S2057.htm PENG Yu-lin, GONG Ai-min, SUN Hai-yan, et al. Study on properties of modified red clay anti-seepage materials in landfill sites[J]. Yangtze River, 2011, 42(S2): 163–165, 169. (in Chinese https://www.cnki.com.cn/Article/CJFDTOTAL-RIVE2011S2057.htm
[13]	王云, 关爱军, 沈峰, 等. 红黏土路堤边坡降水影响深度模拟分析[J]. 公路工程, 2015, 40(5): 50–55. doi: 10.3969/j.issn.1674-0610.2015.05.012 WANG Yun, GUAN Ai-jun, SHEN Feng. Simulation analysis on the precipitation influence depth of laterite soil embankment slope[J]. Highway Engineering, 2015, 40(5): 50–55. (in Chinese) doi: 10.3969/j.issn.1674-0610.2015.05.012
[14]	张文杰, 董林兵, 王顺玉. 山谷型填埋场汇水区域径流补给规律研究[J]. 土木工程学报, 2011, 44(增刊2): 173–176. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC2011S2042.htm ZHANG Wen-jie, DONG Lin-bing, WANG Shun-yu. Investigation on water inflow from the catchment area of a valley-type landfill[J]. China Civil Engineering Journal, 2011, 44(S2): 173–176. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC2011S2042.htm
[15]	冶运涛, 伍靖伟, 王兴奎. 双套环测定土壤渗透系数数值模拟分析[J]. 灌溉排水学报, 2007, 26(3): 14–18. https://www.cnki.com.cn/Article/CJFDTOTAL-GGPS200703003.htm YE Yun-tao, WU Jing-wei, WANG Xing-kui. Numerical simulation of double-ring measuration on hydraulic conductivity[J]. Journal of Irrigation and Drainage, 2007, 26(3): 14–18. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GGPS200703003.htm
[16]	LU N, W LIKOS J W. Unsaturated Soil Mechanics[M]. Beijing: Higher Education Press, 2012.
[17]	张雪东, 赵成刚, 刘艳, 等. 土水特征曲线(SWCC)的滞回特性模拟研究[J]. 工程地质学报, 2010, 18(6): 920–925. doi: 10.3969/j.issn.1004-9665.2010.06.017 ZHANG Xue-dong, ZHAO Cheng-gang, LIU Yan, et al. Simulation and hysteresis model for soil-water characteristic curves[J]. Journal of Engineering Geology, 2010, 18(6): 920–925. (in Chinese) doi: 10.3969/j.issn.1004-9665.2010.06.017

施引文献(20)

期刊类型引用(8)

1.	邹增富，唐世禄，廖云朋，代荣，谢昌志. 直墙拱型引水隧洞洞口段管棚预加固技术分析. 四川建筑. 2024(02): 132-134+138 . 百度学术
2.	张小波，荣耀，张雕，孙洋. 考虑初支滞后性的软弱地层隧道进洞段管棚注浆加固效果研究. 公路. 2024(08): 425-433 . 百度学术
3.	王晓东. 基于最大开挖进尺的管棚设计参数研究. 铁道建筑技术. 2023(03): 40-43 . 百度学术
4.	张小波，张雕，姚池，杨建华，蒋水华，荣耀. 考虑微拱效应的隧道洞口段管棚合理间距确定方法及参数分析. 岩土力学. 2023(06): 1625-1635 . 百度学术
5.	姬云鹏，房灵国，唐昊天，张兴丽，王祥金，白云天，李健，赵红华. 隧道开挖全过程管棚支护受力研究. 现代隧道技术. 2023(03): 123-138 . 百度学术
6.	石钰锋，张涛，曹成威，阳军生，陈祥胜. 基于双参数地基的隧道预支护拱棚壳体力学模型. 工程科学与技术. 2023(04): 142-152 . 百度学术
7.	李奎双，王峥，尹学鑫，程星源. 浅埋软弱地层中超大跨度地铁隧道管棚支护效果研究. 四川建筑. 2023(04): 133-137 . 百度学术
8.	施英，罗春. 隧道管棚超前支护作用机理及其影响因素研究. 人民长江. 2022(11): 130-135 . 百度学术