高压实膨润土与孔隙溶液物理作用机制研究进展

廖饶平; 陈永贵; 刘聪; 叶为民; 乌东北; 王琼

doi:10.11779/CJGE20230811

高压实膨润土与孔隙溶液物理作用机制研究进展 English Version

廖饶平^1,,
陈永贵^{1, 2, ,},
刘聪¹,
叶为民^{1, 2},
乌东北³,
王琼^{1, 2}

1.
同济大学土木工程学院地下建筑与工程系, 上海 200092
2.
土木工程防灾减灾全国重点实验室(同济大学), 上海 200092
3.
同济大学化学科学与工程学院, 上海 200092

基金项目:

国家自然科学基金项目 42125701

国家自然科学基金项目 41977232

国家自然科学基金项目 42030714

上海市教委科研创新计划项目 2023ZKZD26

中央高校基本科研业务费

土木工程Ⅰ类高峰学科建设经费 2022-3-ZD-08

详细信息

作者简介:
廖饶平(1993—)，男，博士研究生，主要从事环境地质与非饱和土力学研究。E-mail: lrp_liao@tongji.edu.cn

通讯作者:
陈永贵, E-mail: cyg@tongji.edu.cn

中图分类号: TU413
计量
- 文章访问数: 605
- HTML全文浏览量: 66
- PDF下载量: 234
出版历程
- 收稿日期: 2023-08-21
- 网络出版日期: 2024-03-24
- 刊出日期: 2024-11-30

Advances in physical interaction mechanism between highly compacted bentonite and pore solution

LIAO Raoping^1,,
CHEN Yonggui^{1, 2, ,},
LIU Cong¹,
YE Weimin^{1, 2},
WU Dongbei³,
WANG Qiong^{1, 2}

1.
Department of Geotechnical Engineering, College of Civil Engineering, Tongji University, Shanghai 200092, China
2.
State Key Laboratory of Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China
3.
School of Chemical Science and Engineering, Tongji University, Shanghai 200092, China

摘要

摘要: 作为高放废物处置库的工程屏障，高压实膨润土在长期服役过程中将受到围岩地下水及其化学成分的渗入作用，膨胀性能不断衰减，最终威胁处置安全。在阐述孔隙溶液对高压实膨润土水化膨胀过程影响规律的基础上，总结了高压实膨润土与孔隙溶液物理作用机制的最新研究成果。结果表明，孔隙溶液对高压实膨润土的物理作用机制包括晶层膨胀、扩散双电层膨胀和吸附作用等3种。其中，孔隙溶液对晶层膨胀的作用与浓度有关，低浓度时表现为促进作用，高浓度时表现为抑制作用，取决于孔隙溶质吸力与临界吸力的差；孔隙溶液对扩散双电层的抑制作用是导致有效孔隙通道扩大、渗透和扩散系数增大的主要原因；pH对膨润土表面活性位点和核素水解的影响是引起吸附特性变化的主要原因，背景离子的竞争吸附作用致使膨润土对核素离子的吸附量显著减少。目前，孔隙溶液作用的参数概化、有效孔隙量化和吸附化学模型的研究仍有不足。因此，进一步优化本构模型中孔隙溶液作用的化学参数，明确不同尺度孔隙的等效量化研究，构建约束条件下压实膨润土多组分竞争吸附模型仍是今后需要深入研究的重点方向。
- 深地质处置库 /
- 高放射性废物 /
- 高压实膨润土 /
- 孔隙溶液 /
- 作用机制 /
- 进展
Abstract: The swelling performance of highly compacted bentonite deteriorates due to infiltrating rock groundwater and chemical components during its long-term operation, ultimately posing a threat to disposal safety. The recent researches on the physical interaction between compacted bentonite and pore solution are reviewed based on the phenomena related to the influences of the pore solution on buffering properties. The results show that the physical mechanism of the pore solution on highly compacted bentonite includes swelling of crystal layer, swelling of diffusion double layers and adsorption effects. The effects of the pore solution on the swelling of crystal layer are related to its concentration, and they are promoted at low concentration and inhibited at high concentration, which depends on the difference between the suction and the critical suction of the poresolution. The inhibition of the pore solution on the swelling of double layers is the main factor for the enlarged pore channels, higher permeability and diffusion coefficients. The pH change and nuclide hydrolysis alter the adsorption characteristics, while the competitive background ion adsorption reduces the nuclide capacity of bentonite. The current shortcomings include parameter generalization for pore solution effects, effective porosity quantification and adsorption models. Therefore, the further optimization of chemical parameters of pore solution in the constitutive model, the clarification of the equivalent quantification of the pores at different scales, and the establishment of a multi-component competitive adsorption model under the constraint of compacted bentonite are still the key directions for the further researches in the future.
- deep geological repository /
- high-level radioactive waste /
- highly compacted bentonite /
- pore solution /
- physical interaction mechanism /
- advance

HTML全文

0. 引言

城市排水管网肩负着城市污水、雨水排放的重任,是现代城市重要的基础设施。但中国大部分市政排水管道铺设于20世纪80年代以前,随着使用年限的临近,管道老化严重,管道下部脱空、沉降、管道腐蚀、渗漏、断裂等灾害不断发生,导致管道无法正常使用^[1],更严重的是因此产生的地面塌陷,地下水污染等事故频发。排水管道的传统修复方式主要是现场整体开挖进行局部修复或管道整体更换,该方式成本昂贵且对城市道路及周围居民生活造成严重影响。高聚物注浆修复技术是一种高效的非开挖修复技术,通过向管道结构脱空渗漏处注射高聚物注浆材料,利用非水反应类高聚物材料本身的膨胀性及固化后的高黏结性,达到填充脱空、封堵渗漏和抬升沉降管道的目的（如图1,2所示）,目前该技术已成功应用于多种地下管道修复项目中^[2-3]。

图 1 高聚物注浆填充管道脱空

Figure 1. Polymer grouting fills the void outside pipe

下载: 全尺寸图片幻灯片

图 2 高聚物注浆抬升管道沉降

Figure 2. Polymer grouting lifts the settlement pipe

下载: 全尺寸图片幻灯片

对于高聚物修复管道的力学性能,目前多利用有限元软件进行模拟或进行试验探究,对修复后管道计算理论和方法的研究相对不足。如徐建国等^[4]采用ABAQUS软件计算了地下管道高聚物注浆修复后的受力情况,结果表明高聚物注浆后,管道内最大拉压应力和最大竖向位移均明显减小；李斌等^[5]对修复后管道在车辆荷载下的力学性能进行了有限元分析,结果显示管道经高聚物修复后,其受力变形恢复到了正常管道水平；张倩倩^[6]利用有限元软件分析了高聚物浆液膨胀扩散对管道应力应变的影响；Wang等^[7]进行了高聚物注浆修复钢筋混凝土排水管道在冲击荷载与车辆荷载下的足尺试验,结果表明高聚物注浆可以有效地修复脱空管道。

通常对埋地管道进行计算分析时,需要考虑管与周围介质的相互作用。国内外学者对此进行了大量研究,如Kjartanson等^[8]通过将管道模化为一系列梁单元,管周的填充土模化为一个非线性材料,从而提出了不开槽埋设管道的管-土相互作用有限元分析法；Lee等^[9]利用有限元模型来反映土壤特性、回填和原位条件,提出了考虑管-土相互作用的地下管道水平挠度方程；Jung等^[10]利用Mohr-Coulomb屈服准则建立了粒状土壤中管-土相互作用的有限元模型；可文海等^[11]利用Pasternak模型模拟管-土相互作用,计算了盾构隧道下穿管道施工引起的管道处土体竖向位移；王小龙等^[12]基于Winkler地基梁理论,建立了埋地钢管在局部悬空时管-土相互作用的力学模型；刘全林等^[13]基于地埋管道受力特征的实测结果,建立了地埋管道与土的相互作用分析组合模型,并给出了模型参数的确定方法。上述文献研究多集中于管-土相互作用,而对于高聚物注浆修复管道则需要考虑管-土-高聚物之间的相互作用。因此本文基于弹塑性理论建立考虑管-土-高聚物相互作用的Winkler模型和Vlazov模型,采用传递矩阵法,编制了MATLAB计算程序,经对地下管道高聚物注浆修复前后的受力变形特性的对比计算,并通过参照现场试验结果,探讨了Winkler模型和Vlazov模型的准确性和合理性。最后采用Vlazov模型分析计算了多工况下管道脱空高聚物注浆修复后的受力变形特性。

1. 管-土-高聚物相互作用分析模型

1.1 管-土-高聚物相互作用Winkler模型

Winker所提出的土介质理想化弾性模型是假设土介质表面任一点处的位移 $ω$ 与作用在该点的应力 $q$ 成正比而与作用在其他各点的应力无关^[14],即

$q (x, y) = k ω (x, y)$ ,

(1)

式中, $q (x, y)$ 为作用在土介质表面的压力, $ω (x, y)$ 为土介质表面的沉降, $k$ 为基底反力系数。

方程（1）是Winker模型的响应函数。该模型将土介质理想化为一系列各自独立的弹簧单元,其弹簧常数为 $k$ 。将土与高聚物都简化为弹性材料即可得出管-土-高聚物相互作用Winkler模型,如图3所示。

图 3 管-土-高聚物相互作用Winkler模型

Figure 3. Winkler model for pipe-soil-polymer interaction

下载: 全尺寸图片幻灯片

对于其中的模型参数 $k_{1}$ ,可按地基承载板试验结果确定,或根据加州承载比试验结果确定。而模型参数 $k_{2}$ 可通过求解无限弹性介质中圆形边界上均匀径向位移,从而求得基底反力系数 $k_{2}$ 。

对于平面应变的轴对称问题,取应力函数^[15]为 $Φ = A \ln r + B r^{2} \ln r + C r^{2} + D$ ,当不计体力时则有

$\begin{array}{l} σ_{r} = \frac{1}{r} \frac{\partial Φ}{\partial r} = \frac{A}{r^{2}} + B (1 + 2 \ln r) + 2 C, \\ σ_{θ} = \frac{\partial^{2} Φ}{\partial r^{2}} = - \frac{A}{r^{2}} + B (3 + 2 \ln r) + 2 C, \\ τ_{r θ} = 0 。 \end{array}}$

(2)

式中, $Φ$ 为应力系数, $A, B, C, D$ 为系数, $r$ 为半径, $σ_{r}$ , $σ_{θ}$ 分别为径向应力与切向应力。

由位移单值性可得 $B = 0$ ,其次在距离原点无穷远处,由圣维南原理可知 $σ_{r} = σ_{θ} \overset{}{\to} 0$ 从而可得 $C = 0$ 。将 $B = 0$ , $C = 0$ 代入式（2）化简可得

$\begin{array}{l} σ_{r} = \frac{1}{r} \frac{\partial Φ}{\partial r} = \frac{A}{r^{2}}, \\ σ_{θ} = \frac{\partial^{2} Φ}{\partial r^{2}} = - \frac{A}{r^{2}} 。 \end{array}}$

(3)

由平面应变的胡克定律,并将式（3）代入可得

$\begin{array}{l} ε_{r 1} = \frac{1}{E_{g}} (σ_{r} - ν_{g} σ_{θ}) = \frac{A (1 + ν_{g})}{E_{g} r^{2}}, \\ ε_{r 2} = \frac{1}{E_{s}} (σ_{r} - ν_{s} σ_{θ}) = \frac{A (1 + ν_{s})}{E_{s} r^{2}} 。 \end{array}}$

(4)

式中, $ε_{r 1}$ , $ε_{r 2}$ 分别为高聚物层与土层的应变, $E_{g}$ , $ν_{g}$ 分别为高聚物的弹性模量与泊松比, $E_{s}$ , $ν_{s}$ 分别为土的弹性模量与泊松比。

圆形边界上均匀径向位移可为

$u = \int_{R}^{R + D} ε_{r 1} d r + \int_{R + D}^{\infty} ε_{r 2} d r$ ,

(5)

式中,D为高聚物层厚度。

将式（4）代入式（5）从而得到径向位移：

$u = \frac{A}{R + D} (\frac{1 + ν_{s}}{E_{s}} - \frac{1 + ν_{g}}{E_{g}}) + \frac{A (1 + ν_{g})}{E_{g} R}$ 。

(6)

按Winker模型中弹性参数 $k$ 的定义,将式（3）,（6）代入,则得到高聚物层的基底反力系数为

$k_{2} = \frac{σ_{r}}{u} = \frac{1}{\frac{R^{2}}{R + D} (\frac{1 + ν_{s}}{E_{s}} - \frac{1 + ν_{g}}{E_{g}}) + \frac{R (1 + ν_{g})}{E_{g}}}$ 。

(7)

1.2 管-土-高聚物相互作用Vlazov模型

Winkler模型由于在描述自然土体的连续性态中所固有的缺陷以及Winkler模型只适用于地基土（压缩层）很薄的情况,而当压缩层较厚时就会产生较大误差。为了改进Winkler模型的上述缺陷,通过在弹性连续介质模型中引入约束或简化位移分布与应力的某些假设,从而得到一种新的模型,称为双参数模型^[14]。Vlazov模型是一种双参数模型,它所列出的模型公式是以变分法的应用为基础,通过引进一些能简化各向同性线弹性连续介质基本方程的位移约束而导出的。其响应函数为

$q (x, y) = K ω (x, y) - T \nabla^{2} ω (x, y)$ ,

(8)

式中, $q (x, y)$ 为作用在土介质表面的压力, $ω (x, y)$ 为土介质表面的沉降, $\nabla^{2}$ 为拉普拉斯算子, $K$ 为压应力作用下对土介质变形的度量,T为对邻近土单元体传递性的度量。

该模型也适用于非均匀的弹性层,将土与高聚物都简化为弹性材料即可得出管-土-高聚物相互作用分析Vlazov模型,如图4所示。

图 4 管-土-高聚物相互作用Vlazov模型

Figure 4. Vlazov model for pipe-soil-polymer interaction

下载: 全尺寸图片幻灯片

可通过求解双参数模型的地基控制方程来求得模型参数。弹性层在平面应变状态下,在 $x$ - $z$ 平面内假设土层中任一点M的位移满足

$\bar{u} (x, z) = 0, \bar{ω} (x, z) = ω (x) φ (z), [φ (0) = 1, φ (H) = 0]$ ,

(9)

式中, $φ (z)$ 为描述 $z$ 方向位移 $ω (x, z)$ 的变化函数。

利用平面应变条件下的应力与应变关系可得

$\begin{array}{l} σ_{z 1} = \frac{E_{1}}{1 - ν_{1}^{2}} ω (x) φ^{'} (z), \\ τ_{x z 1} = \frac{E_{1}}{2 (1 + ν_{1})} ω^{'} (x) φ (z), \\ σ_{z 2} = \frac{E_{2}}{1 - ν_{2}^{2}} ω (x) φ^{'} (z), \\ τ_{x z 2} = \frac{E_{2}}{2 (1 + ν_{2})} ω^{'} (x) φ (z) 。 \end{array}}$

(10)

式中 $E_{1} = E_{g} / (1 - ν_{g}^{2})$ , $ν_{1} = ν_{g} / (1 - ν_{g})$ ； $E_{2} = E_{s} /$ $(1 - ν_{s}^{2})$ , $ν_{2} = ν_{s} / (1 - ν_{s})$ ； $E_{g}$ , $ν_{g}$ 为高聚物的弹性模量与泊松比, $E_{s}$ , $ν_{s}$ 为土的弹性模量与泊松比。

利用拉格朗日虚功原理,取虚位移 $δ \bar{u} = 0$ , $δ \bar{ω} =$ $φ (z) δ ω$ ,得到外力所作的虚功为

$\begin{matrix} U_{e} = b q (x) φ (0) δ ω d x + b \int_{0}^{H_{1}} \frac{\partial τ_{x z 1}}{\partial x} φ (z) δ ω d x d x + \\ b \int_{H_{1}}^{H} \frac{\partial τ_{x z 2}}{\partial x} φ (z) δ ω d x d z \end{matrix}$ 。

(11)

内力所作的虚功为

$U_{i} = - b \int_{0}^{H_{1}} σ_{z 1} φ^{'} (z) δ ω d x d z - b \int_{H_{1}}^{H} σ_{z 2} φ^{'} (z) δ ω d x d z$ 。

(12)

由 $U_{e} + U_{i} = 0$ ,代入式（10）并消去 $δ ω d x$ 可得

$q (x) = K ω (x) - 2 T \frac{d^{2} ω (x)}{d x^{2}}$ 。

(13)

式中, $K = \frac{E_{1}}{1 - ν_{1}^{2}} \int_{0}^{H_{1}} {[φ^{'} (z)]}^{2} d z + \frac{E_{2}}{1 - ν_{2}^{2}} \int_{H_{1}}^{H} {[φ^{'} (z)]}^{2} d z$ , $T =$ $\frac{E_{1}}{4 (1 + ν_{1} ）} \int_{0}^{H_{1}} φ^{2} （ z) d z + \frac{E_{2}}{4 (1 + ν_{2} ）} \int_{H_{1}}^{H} φ^{2} （ z) d z$ 。

对于函数 $φ (z)$ ,Vlazov提出了线性与非线性表达式分别为

$φ_{1} (z) = 1 - z / H, φ_{2} （ z) = \frac{sh [γ (H - z) / L]}{sh [γ h / L]} 。$

式中 H为计算弹性层厚度,L为地埋管道计算长度, $γ$ 为反映基础梁特性的常数,通常取1.5。本文中对于高聚物层,取线性表达式,对于土层,采用非线性表达式。将其代入K,T计算表达式可得模型参数如下：

$K_{1} = \frac{E_{2}}{H {(1 - ν_{2})}^{2}} \cdot \frac{γ H}{2 L} (\frac{sh (γ H / L) ch (γ H / L) + γ H / L}{{sh}^{2} (γ H / L)}),$

$T_{1} = \frac{E_{2} H}{12 (1 + ν_{2})} \cdot \frac{3 L}{2 γ H} (\frac{sh (γ H / L) ch (γ H / L) - γ H / L}{{sh}^{2} (γ H / L)}),$

$K_{2} = \frac{E_{1}}{H_{1} {(1 - ν_{1})}^{2}} + \frac{E_{2}}{H {(1 - ν_{2})}^{2}} \cdot \frac{γ H}{2 L} \cdot [sh (γ H / L - γ H_{1} / L) \cdot$

$ch (γ H / L - γ H_{1} / L) + γ H / L - γ H_{1} / L] / {sh}^{2} (γ H / L),$

$T_{2} = E_{1} H_{1} / 12 (1 + ν_{1}) + E_{2} H / 12 (1 + ν_{2}) \cdot 3 L / 2 γ H \cdot [sh (γ H / L -$

$γ H_{1} / L) \cdot ch(γ H / L - γ H_{1} / L) - γ H / L + γ H_{1} / L)] / {sh}^{2} (γ H / L) 。$

2. 基于传递矩阵法的管道平面应变问题分析

2.1 传递矩阵法

传递矩阵法是以受力杆件的基本微分方程为基础,采用矩阵相乘来求出杆件各截面的内力和位移值,从而对杆件进行力学分析的一种方法^[16]。其计算过程为

$S_{i} = U_{i} U_{i - 1} \dots U_{2} U_{1} S_{0}$ ,

(14)

式中, $S_{0}$ 为初始状态矢量, $S_{i}$ 为末端状态矢量, $U_{1}$ , $U_{2}$ , $\dots$ , $U_{i}$ , $U_{i - 1}$ 为各段的传递矩阵。

2.2 管道计算模型描述

为计算简便,将管道简化为平面应变问题,忽略管道横截面的不圆度,将其假设为圆环形,并且考虑作用在管道上的荷载为对称荷载,因此只需对管道的一半进行计算分析。如图5所示,取单位长度管道,管道壁厚为t,半径为R,将管道的一半划分为n个单元,同时由于管道的不同位置,管土的变形趋势不同并且不同位置土体的压实程度也会有区别,所以与管道相接触的高聚物和土体也相应的被划分为三部分,分别为：底部支承区、胸腔区、管顶区。

图 5 管道计算模型

Figure 5. Computational model for buried pipe

下载: 全尺寸图片幻灯片

将划分的n个单元都假设为梁单元,并按顺序刚性连接,取第 $i$ 单元作为研究对象,其受力简图如图6所示。

图 6 单元受力简图

Figure 6. Forces acting on elements

下载: 全尺寸图片幻灯片

图6中, $p_{i}$ , $q_{i}$ 为 $i$ 节点处的径向荷载和切向荷载, $p_{i + 1}$ , $q_{i + 1}$ 为 $i + 1$ 节点处的径向荷载和切向荷载, $l$ 为单元 $i$ 的长度, $θ$ 为 $i$ 单元与总体坐标系中 $X$ 轴的夹角。

设在外力作用下,管道 $i$ 单元中 $i$ 节点处的切向位移为 $u_{i}$ ,径向位移为 $v_{i}$ ,转角为 $ω_{i}$ ,弯矩为 $M_{i}$ ,轴力为 $N_{i}$ ,剪力为 $T_{i}$ 。在局部坐标系中,由平衡条件和物理方程可得

$\begin{array}{l} \frac{d T}{d x} = - K^{'} y + 2 T \frac{d^{2} y (x)}{d x^{2}} - p_{i} (1 - \frac{x}{l}) - p_{i + 1} \frac{x}{l}, \\ \frac{d N}{d x} = q_{i} (1 - \frac{x}{l}) + q_{i + 1} \frac{x}{l}, \\ \begin{matrix} \frac{d M}{d x} = T & , & \frac{d u}{d x} = - \frac{N}{E S}, \end{matrix} \\ \begin{matrix} \frac{d^{2} v}{d x^{2}} = \frac{M}{E I} & , & \frac{d v}{d x} = ϖ 。 \end{matrix} \end{array}}$

(15)

求解式（15）,并满足边界条件：当 $x = 0$ 时, $u = u_{i}$ , $ν = ν_{i}$ , $ω = ω_{i}$ , $M = M_{i}$ , $N = N_{i}$ , $T = T_{i}$ ；当 $x = l$ 时, $u = u_{i + 1}$ , $ν = ν_{i + 1}$ , $ω = ω_{i + 1}$ , $M = M_{i + 1}$ , $N = N_{i + 1}$ , $T = T_{i + 1}$ 。可得

${S_{i + 1}} = [U_{i, i + 1}] {S_{i}}$ ,

(16)

式中, ${S_{i}} = {u_{i}, ν_{i}, ω_{i}, M_{i}, N_{i}, T_{i}, 1}$ 为 $i$ 节点处的内力和位移列阵, $[U_{i, i + 1}]$ 为 $7 \times 7$ 的传递矩阵, ${S_{i + 1}} =$ ${u_{i + 1}, ν_{i + 1}, ω_{i + 1}, M_{i + 1}, N_{i + 1}, T_{i + 1}, 1}$ 为 $i + 1$ 节点处的内力和位移列阵。

通过坐标转换矩阵 $R$ 可将式（16）变为整体坐标系下的传递关系：

${{\bar{S}}_{i + 1}} = [{\bar{U}}_{i, i + 1}] {{\bar{S}}_{i}} 。$

(17)

式中 $[{\bar{U}}_{i, i + 1}] = R [U_{i, i + 1}] R^{- 1}$ ,为 $i$ 单元的传递矩阵, ${{\bar{S}}_{i}}$ 、 ${{\bar{S}}_{i + 1}}$ 分别为整体坐标系下 $i$ 节点与 $i + 1$ 节点处的内力和位移列阵, $R$ 为 $i$ 单元的坐标转换矩阵。已知各单元的传递矩阵后,将其代入式（14）可得

${{\bar{S}}_{n}} = [{\bar{U}}_{n, n + 1}] [{\bar{U}}_{n - 1, n}] \dots [{\bar{U}}_{2, 3}] [{\bar{U}}_{1, 2}] {{\bar{S}}_{1}}$ 。

(18)

由 $u_{1}$ = $ω_{1}$ = $T_{1}$ , $u_{n}$ = $ω_{n}$ = $T_{n}$ ,代入式（18）即可求得 $ν_{1}, M_{1}, N_{1}$ 和 $ν_{n}, M_{n}, N_{n}$ 。

再以 $1$ 节点为已知量,代入

${{\bar{S}}_{i}} = [{\bar{U}}_{i, i + 1}] [{\bar{U}}_{i - 1, i}] \dots [{\bar{U}}_{2, 3}] [{\bar{U}}_{1, 2}] {{\bar{S}}_{1}}$ ,

(19)

即可求得任意节点的内力和位移。求出节点的位移后,可以确定各单元相应的接触压力。结合土体的弹-塑性变形法则和管道本身的变形量,部分接触压力将被限制以满足土的塑性条件,然后经迭代计算使之满足所有的边界条件和变形准则。

2.3 程序实现

根据上述计算过程,利用MATLAB编制了管道应力计算程序,其计算流程图如图7所示。

图 7 管道应力计算程序流程图

Figure 7. Flow chart of stress program for pipe

下载: 全尺寸图片幻灯片

3. 试验概述及数据处理

管道试验在室外现场区域进行,开挖长、宽、深为8 m×2.3 m×2.3 m的试验坑。选择3个完整的管段和两个半管段,有效长度为8 m。管道放置到沟槽中后,开始回填,每层厚约0.2 m,用小夯夯实土壤。然后对同一根管道先后进行正常、脱空、注浆修复等工况的试验。通过在管内外壁黏贴应变片来测量荷载作用下管道的内外应变。试验所用管道、土壤和高聚物性能参数如表1所示,试验管道剖面图及现场试验图如图8所示。

表 1 试验所用管道、土壤和高聚物性能参数表

Table 1. Performance parameters of pipes, soil and polymer

材料	密度/(kg·m^-3)	弹性模量/MPa	泊松比	黏聚力/kPa	内摩擦角/(°)
管道	2400	30000	0.3	—	—
填土	1980	60	0.3	40	30
高聚物	240	20	0.2	—	—

下载: 导出CSV

| 显示表格

图 8 试验管道剖面图及现场试验图

Figure 8. Sectional drawing of test pipes and pictures of field tests

下载: 全尺寸图片幻灯片

为了将试验值与模型计算结果进行更好的对比,需要对试验数据进行一定处理。根据文献[17]可知管道弯矩与内外管壁应变之间的关系如下：

$M = E I (\frac{ε_{in} - ε_{out}}{h}),$

(20)

式中,M为管道弯矩,E为管道弹性模量,I为截面惯性矩, $ε_{in}$ , $ε_{out}$ 分别为内外管壁应变,h为管壁厚度。

在管道设计中,管道自身承受外部荷载的能力通过管道刚度来定义,其与管壁的几何形状和管道材料类型有关。同时在ASTM D 2412—11^[18]中,管道刚度还被定义为：使管道产生一定百分比变形（即5%）的荷载与沿荷载施加方向管道内径变形量的比值,因此已知管道刚度和荷载大小可以根据

$Δ y = \frac{F}{P S}$

(21)

近似求出相应的管道径向变形。式中, $Δ y$ 为管道径向变形, $F$ 为管壁上的径向应力, $P S$ 为管道刚度。

4. 模型计算及与试验结果对比

依据试验情况,取静荷载作用管段为计算管道,采用不同的管-土-高聚物相互作用模型,其计算结果如图9,10所示。图9（a）,（b）分别表示正常、脱空和高聚物修复情况下,采用Winkler模型计算所得的管道径向变形和环向轴力对比图,同样的图10（a）,（b）表示采用Vlazov模型计算所得的结果图。从图9,10可知,修复后管道的径向变形会略大于正常值,而轴力与正常值基本一致,正常与高聚物修复后的径向变形最大差值都出现在管道底部,分别为0.4 mm（采用Winkler模型）与1 mm（采用Vlazov模型）。由此可见,高聚物注浆后管道的变形与力学性能可以恢复到正常水平。在脱空条件下,管道底部的径向变形与管顶的轴力会出现明显的增大现象,与正常条件相比,底部的径向变形增大了3.9倍（采用Winkler模型）与5.3倍（采用Vlazov模型）,同样的,顶部的环向轴力分别增大了1.6倍和2.0倍。说明底部脱空对于管道的力学特性会有很大的影响。

图 9 Winkler模型计算结果

Figure 9. Calculated results of Winkler model

下载: 全尺寸图片幻灯片

图 10 Vlazov模型计算结果

Figure 10. Calculated results of Vlazov model

下载: 全尺寸图片幻灯片

表2为正常与修复条件下管-土-高聚物相互作用分析模型参数的比较。由表2可知,Winkler模型计算出的修复后地基反力系数为正常值的66%,而Vlazov模型达到了正常值的91%,说明Vlazov模型计算出的修复后结果会与正常值更为接近。图11给出了采用不同分析模型得到的径向变形和环向弯矩与试验值的对比。由图11（a）,（b）可以看出,模型计算得出的径向变形和环向弯矩变化趋势图与试验值基本一致,与试验值相比最大径向变形差值都为2 mm,误差较小,说明建立的计算方法是可靠的。并且Vlazov模型计算结果明显比Winkler模型计算结果大,这是由于Vlazov模型增加了土体的侧向传力作用,考虑了土介质的荷载传递性,比Winkler模型合理。由图11（a）可知,采用Vlazov模型计算的径向变形与试验值较接近,除管顶处,计算误差在2%～25%,在45°管腰处,模型计算结果与试验值最接近,Winkler模型计算结果与试验值误差较大,达9%～55%。由图11（b）可知,在环向弯矩计算中,Winkler模型计算结果与试验值的误差在管底、管侧、管顶分别为35%,21%,12%,同理Vlazov模型与试验值的误差分别为18%,16%,11%,因此,采用Vlazov模型计算的环向弯矩与试验值较接近。综上所述,相比于Winkler模型,Vlazov模型与管道实际受力情况符合得更好。

表 2 相互作用分析模型参数

Table 2. Parameters of interaction analysis model

模型参数/(N·m^-3)	底部支承区（修复）	底部支承区（正常）	$\frac{修复}{正常} / %$
Winkler模型K	68945	105230	66
Vlazov模型K	82842	91434	91
Vlazov模型T	66	91	73

下载: 导出CSV

| 显示表格

图 11 模型计算结果与试验对比

Figure 11. Comparison between model and test results

下载: 全尺寸图片幻灯片

5. 不同因素对修复后管道内力和变形的影响

埋地管道的内力和变形受到管道埋深、管径、管道材料、地面荷载及周围填土的性质等诸多因素的影响,同时由于利用Vlazov模型来模拟管-土-高聚物的相互作用更加符合实际情况,因此下文将采用Vlazov模型对管道脱空高聚物注浆修复后不同工况下的管道进行计算。

5.1 不同埋深的影响

图12分别为覆土深度是1.0,1.5,2 m时修复后管道的径向变形和环向轴力对比图。从图中可以看出,不同埋深时,径向变形和轴力的变化趋势是一致的,随着埋深增大,其数值会明显增大。在图12（a）中,由于选用坐标系的不同,虽然管道底部变形为正值,顶部为负值,但二者都是向下的,顶部变形量会大于底部变形量。而对于管道的环向弯矩和轴力,从图12（b）可看出,轴力分布在管顶处呈压缩状态,在管腰下部呈拉伸状态,且轴向最大拉应力大于最大压应力。

图 12 不同埋深下管道的计算值

Figure 12. Calculated values of pipes at different buried depths

下载: 全尺寸图片幻灯片

5.2 不同附加荷载的影响

图13为管顶施加不同大小荷载时,修复后管道的径向变形和环向轴力对比图。随着管道附加荷载的增加,管顶处的径向变形以及管腰处的轴力都出现了明显的增大。从图13（a）可知,荷载为0.1,0.15 MPa时管底的径向变形分别为0.05 MPa时管底径向位移的1.1倍和1.2倍；其管顶的径向变形分别为0.05 MPa时管顶径向位移的1.3倍和1.5倍,由此可见,同样荷载变化下,管顶的变形变化率要大于管底,这是因为相对于管顶,管底的埋深更大,而随着土壤厚度的增加,与土压力相比,附加载荷对管道的影响会变小,同时管土相互作用加强,这也在一定程度上限制了管道的变形。

图 13 不同荷载下管道的计算值

Figure 13. Calculated values of pipes under different loads

下载: 全尺寸图片幻灯片

5.3 不同土壤弹性模量的影响

图14为改变回填土的弹性模量时,修复后管道的径向变形和环向弯矩对比图。从图14（a）可知,在管腰以下,高聚物注浆修复区的径向变形随回填土弹性模量的增大而增大,在非脱空区则呈现相反趋势,图14（b）中的弯矩也有相同的变化趋势。由此可见,非脱空区的管土相互作用随着土体弹性模量的增大而加强,但在高聚物注浆修复区,由于高聚物的加入,与原来的土层一起形成了复合土层,管道与复合土层之间的相互作用则随着土体弹性模量的增大而减弱。从图14（b）可以看出,在非脱空区,管腰下部和管顶处的弯矩呈相反的变化趋势,这是因为,管道内力随回填土弹性模量的变化而变化是两个方面效应共同作用的结果^[19],增大管周土体弹性模量一方面增加了土体对管道的弹性抗力,增强了管土相互作用,使管壁内力分布更加均匀,从而降低了管道内力,但另一方面也增大了管道基础的刚度导致了更大的管顶竖向应力,这会使管道内力随之增大。在管顶区,由于埋深较浅,管顶竖向应力起主导作用,因此图14（b）中管顶区的弯矩随回填土弹性模量的增大而增大,而在管腰以下,埋深相对较深,管土相互作用起主导作用,弯矩也就随之减小。

图 14 不同土壤弹性模量下管道的计算值

Figure 14. Calculated values of pipes under different elastic moduli of soil

下载: 全尺寸图片幻灯片

5.4 不同管壁厚度的影响

图15为不同管壁厚度下,修复后管道的径向变形和环向弯矩对比图。从图中可以看出,管道的环向弯矩随着管壁厚度的增大而增大,径向变形则呈相反趋势。这是由于,同一直径管道,管壁厚度不同,管道的抗弯刚度也就不同,增加管壁厚度实际就是增加了管道的截面刚度,从而导致管道径向变形的减小,但刚度的增大也会使管土之间的变形协调能力降低,使得管顶竖向应力增大,从而增大了管周的应力分布,因此管道的弯矩也随之增大。

图 15 不同管壁厚下管道的计算值

Figure 15. Calculated values of pipes with different pipe thicknesses

下载: 全尺寸图片幻灯片

6. 结论

为了对高聚物注浆修复脱空渗漏的地下管道开展受力变形特性分析,本文基于弹塑性理论分别建立了考虑管-土-高聚物相互作用的Winkler模型和Vlazov模型,采用传递矩阵法,编制了MATLAB计算程序,经对地下管道高聚物注浆修复前后的受力变形特性的对比计算,并通过参照现场试验结果,探讨了Winkler模型和Vlazov模型的准确性和合理性。最后采用Vlazov模型开展了管道脱空高聚物注浆修复后受力变形特性研究,工况包括管道埋深、荷载集度、回填土弹性模量、管道壁厚等因素,得到以下6点结论。

（1）考虑管-土-高聚物相互作用的Winkler模型和Vlazov模型计算得出的径向变形和环向弯矩计算值与试验值基本一致,验证了模型的准确性和合理性,Vlazov模型计算结果与试验值更为接近。

（2）管道底部的脱空对管道的整体受力变形力学特性有很大影响,底部脱空区使管道的径向变形和轴向内力明显增大。但经高聚物注浆修复后,管道的径向变形和轴向内力接近正常管道,说明高聚物注浆具有较好的修复管道功能。

（3）随着管道埋深和附加荷载的增大,管道的径向变形、环向弯矩和环向轴力都有不同程度的增大。

（4）管道的径向变形随回填土弹性模量的增大而减小,随着管道周围填充土弹性模量的增大,管土相互作用也随之加强,更有利于发挥土介质对外荷载的支承作用,从而对管道的整体受力有利。

（5）管道的环向弯矩和轴力随着管壁厚度的增大而增大,径向变形则呈相反趋势,管壁厚度的增加使得管道截面刚度增大,从而导致管道径向变形的减小。

（6）高聚物注浆修复技术是一种高效的非开挖修复技术,通过向管道结构脱空渗漏处注射高聚物注浆材料,利用非水反应类高聚物材料本身的膨胀性及固化后的高黏结性,达到填充脱空、封堵渗漏和抬升沉降管道的目的,目前该技术已成功应用于地下管道的渗漏脱空修复工程。本研究为管道高聚物注浆修复设计施工提供了一定的理论参考。

图 1 高放废物深地质处置库及多重屏障缓冲示意图

Figure 1. Schematic diagram of deep geological repository and multiple barrier buffer for high-level waste

下载: 全尺寸图片幻灯片

图 2 高压实膨润土微观结构组成^[24]

Figure 2. Microstructural composition of highly compacted bentonite^[24]

下载: 全尺寸图片幻灯片

图 3 不同吸力下压实膨润土膨胀力时程曲线^[31]

Figure 3. Time-swelling pressure curves with increase of suction^[31]

下载: 全尺寸图片幻灯片

图 4 吸力对蒙脱石水化层间距和层叠体厚度的影响^[25-26]

Figure 4. Effects of suction on interlayer distance and laminate thickness of montmorillonite^[25-26]

下载: 全尺寸图片幻灯片

图 5 膨润土水化膨胀机制

Figure 5. Hydration swelling mechanism of bentonite

下载: 全尺寸图片幻灯片

图 6 溶液作用下膨胀力三阶段演化特征^[30]

Figure 6. Three stages of swelling pressure hydrated with salt solutions^[30]

下载: 全尺寸图片幻灯片

图 7 溶液对膨润土最终膨胀力的影响^[35-39]

Figure 7. Effects of ion species and concentration on final swelling pressure of compacted bentonite^[35-39]

下载: 全尺寸图片幻灯片

图 8 层叠体颗粒受力平衡示意图

Figure 8. Schematic diagram of stress balance of layered particles

下载: 全尺寸图片幻灯片

图 9 离子类型和浓度对膨润土膨胀指数的影响^{[15, 43]}

Figure 9. Effects of ion species and concentration on free swell index of bentonite^{[15, 44]}

下载: 全尺寸图片幻灯片

图 10 溶液浓度与渗透吸力的关系^[45]

Figure 10. Relationship between concentration of different solutions and osmotic suction^[45]

下载: 全尺寸图片幻灯片

图 11 离子类型对高压实膨润土膨胀力的影响^[16]

Figure 11. Effects of ion species on final swelling pressure of compacted bentonite^[16]

下载: 全尺寸图片幻灯片

图 12 离子类型对高压实膨润土膨胀变形的影响^[15]

Figure 12. Effects of ion species on final swelling deformation of compacted bentonite^[15]

下载: 全尺寸图片幻灯片

图 13 Stern理论扩散双电层构造^[45]

Figure 13. Diffusion double-layer structure of Stern theory^[45]

下载: 全尺寸图片幻灯片

图 14 离子类型及浓度对膨润土饱和渗透系数的影响^[46]

Figure 14. Effects of ion species and concentrations on saturation permeability coefficient of compacted bentonite^[46]

下载: 全尺寸图片幻灯片

图 15 离子强度对膨润土中核素离子扩散的影响^[47]

Figure 15. Effects of ionic strength on diffusion of nuclides^[47]

下载: 全尺寸图片幻灯片

图 16 蒙脱石表面吸附作用示意图^[49]

Figure 16. Schematic diagram of surface adsorption of montmorillonite^[49]

下载: 全尺寸图片幻灯片

图 17 pH对GMZ-Na基膨润土吸附核素的影响^[52-56]

Figure 17. Effects of pH on adsorption of nuclides by GMZ-Na bentonite^[52-56]

下载: 全尺寸图片幻灯片

图 18 不同温度中U(Ⅵ)离子形态分布随pH的变化^[49]

Figure 18. Species distribution of U(Ⅵ) ions with pH change under different temperatures^[49]

下载: 全尺寸图片幻灯片

图 19 离子强度对GMZ-Na基膨润土吸附U(Ⅵ)的影响^[54]

Figure 19. Effects of ionic strength on adsorption of U(Ⅵ) ions by GMZ-Na bentonite^[54]

下载: 全尺寸图片幻灯片

图 20 阴阳离子类型对GMZ-Na基膨润土吸附U(Ⅵ)的影响^[54]

Figure 20. Effects of cationic and anion species on adsorption of U(Ⅵ) by GMZ-Na bentonite^[54]

下载: 全尺寸图片幻灯片

图 21 背景离子类型对U(Ⅵ)吸附分配系数的影响^[60]

Figure 21. Effects of ion species on equilibrium constants of U(Ⅵ) ^[60]

下载: 全尺寸图片幻灯片

表 1 Na基膨润土表面活性位点的反应机制^[51]

Table 1 Reactions of active sites on Na-bentonite^[51]

表面活性位点反应	双层模型参数(lgK)	非静电模型参数(lgK)
$\equiv X{\text{Na}} \rightleftharpoons \equiv {X^ - } + {\text{N}}{{\text{a}}^ + }$	-1.58	-1.58
$\equiv X{\text{Na}} + {{\text{H}}^ + }f \equiv XH + {\text{N}}{{\text{a}}^ + }$	1.50	2.10
$\equiv {\text{AlOH}} + {{\text{H}}^ + }f \equiv {\text{AlOH}}_{\text{2}}^{\text{ + }}$	6.15	5.83
$\equiv {\text{AlOH}} + {\text{O}}{{\text{H}}^ - }f \equiv {\text{AlO + H}}_{\text{2}}^{}{\text{O}}$	-9.27	-7.02
$\equiv {\text{SiOH}} + {\text{O}}{{\text{H}}^ - }f \equiv {\text{Si}}{{\text{O}}^ - }{\text{ + H}}_{\text{2}}^{}{\text{O}}$	-9.06	-8.75

下载: 导出CSV

表 2 GMZ-Na基膨润土表面活性位点的分布密度^[50]

Table 2 Distribution density of active sites on surface of GMZ-Na bentonite ^[50]

活性位点类型	密度
蒙脱石层间阳离子交换位点，（ $\equiv {\rm X^ - }$ ）	1.35×10^-5 mol/m², (66.5±7 cmol/kg)
铝氧八面体边缘羟基位点，（ $\rm \equiv AlOH$ ）	9.39×10^-7 mol/m², (4.61 cmol/kg)
硅氧四面体边缘羟基位点，（ $\rm \equiv SiOH$ ）	1.88×10^-6 mol/m², (9.23 cmol/kg)

下载: 导出CSV

参考文献(60)

[1]	TAKAYAMA Y, TACHIBANA S, IIZUKA A, et al. Constitutive modeling for compacted bentonite buffer materials as unsaturated and saturated porous media[J]. Soils and Foundations, 2017, 57(1): 80-91. doi: 10.1016/j.sandf.2017.01.006
[2]	LU P H, YE W M, HE Y. A constitutive model of compacted bentonite under coupled chemo-hydro-mechanical conditions based on the framework of the BExM[J]. Computers and Geotechnics, 2023, 158: 105360. doi: 10.1016/j.compgeo.2023.105360
[3]	BOURG I C, BOURG A C M, SPOSITO G. Modeling diffusion and adsorption in compacted bentonite: a critical review[J]. Journal of Contaminant Hydrology, 2003, 61(1/2/3/4): 293-302.
[4]	GUO G L, FALL M. Advances in modelling of hydro-mechanical processes in gas migration within saturated bentonite: a state-of-art review[J]. Engineering Geology, 2021, 287: 106123. doi: 10.1016/j.enggeo.2021.106123
[5]	GUIMARãES L D, GENS A, SÁNCHEZ M, et al. A chemo-mechanical constitutive model accounting for cation exchange in expansive clays[J]. Géotechnique, 2013, 63(3): 221-234. doi: 10.1680/geot.SIP13.P.012
[6]	ALONSO E E, VAUNAT J, GENS A. Modelling the mechanical behaviour of expansive clays[J]. Engineering Geology, 1999, 54(1/2): 173-183.
[7]	DOMINIJANNI A, MANASSERO M, PUMA S. Coupled chemical-hydraulic-mechanical behaviour of bentonites[J]. Géotechnique, 2013, 63(3): 191-205. doi: 10.1680/geot.SIP13.P.010
[8]	BENNETHUM L S, MURAD M A, CUSHMAN J H. Macroscale thermodynamics and the chemical potential for swelling porous media[J]. Transport in Porous Media, 2000, 39(2): 187-225. doi: 10.1023/A:1006661330427
[9]	徐永福. 膨胀土的水力作用机理及膨胀变形理论[J]. 岩土工程学报, 2020, 42(11): 1979-1987. doi: 10.11779/CJGE202011002 XU Yongfu. Hydraulic mechanism and swelling deformation theory of expansive soils[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(11): 1979-1987. (in Chinese) doi: 10.11779/CJGE202011002
[10]	RUAN K L, KOMINE H, ITO D, et al. Hydraulic conductivity and X-ray diffraction tests of unsaturated bentonites with a multi-ring and their predictions by pores distributions[J]. Engineering Geology, 2022, 306: 106738. doi: 10.1016/j.enggeo.2022.106738
[11]	KOMINE H. Theoretical equations on hydraulic conductivities of bentonite-based buffer and backfill for underground disposal of radioactive wastes[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2008, 134(4): 497-508. doi: 10.1061/(ASCE)1090-0241(2008)134:4(497)
[12]	TORSTENFELT B, ALLARD B, KIPATSI H. Measurements of ion mobilities in clay[J]. Soil Science, 1985, 139(6): 512-516. doi: 10.1097/00010694-198506000-00006
[13]	MUURINEN A, PEMTILÄ-HILTUNEN P, UUSHEIMO K. Diffusion of chloride and uranium in compacted sodium bentonite[J]. MRS Online Proceedings Library, 1988, 127(1): 743-748.
[14]	TRIPATHY S, SRIDHARAN A, SCHANZ T. Swelling pressures of compacted bentonites from diffuse double layer theory[J]. Canadian Geotechnical Journal, 2004, 41(3): 437-450. doi: 10.1139/t03-096
[15]	陈龙. 复杂化学环境下膨润土膨胀变形研究[D]. 绵阳: 西南科技大学, 2021. CHEN Long. Study on the Swelling Deformation of Bentonite in Complex Chemical Environment[D]. Mianyang: Southwest University of Science and Technology, 2021. (in Chinese)
[16]	张艺东. 碱溶液及离子溶液作用下高庙子膨润土膨胀性及渗透性研究[D]. 绵阳: 西南科技大学, 2017. ZHANG Yidong. Study on the Permeability and Expansibility of Alkali Solution and Ion Solution of Gaomiaozi Bentonite[D]. Mianyang: Southwest University of Science and Technology, 2017. (in Chinese)
[17]	VEGA F A, COVELO E F, ANDRADE M L. Competitive sorption and desorption of heavy metals in mine soils: influence of mine soil characteristics[J]. Journal of Colloid and Interface Science, 2006, 298(2): 582-592. doi: 10.1016/j.jcis.2006.01.012
[18]	POLCARO A M, MASCIA M, PALMAS S, et al. Competitive sorption of heavy metal ions by soils[J]. Environmental Engineering Science, 2003, 20(6): 607-616. doi: 10.1089/109287503770736122
[19]	高子瑞, 陈涛, 徐永福. 盐溶液对膨润土膨胀性的影响[J]. 岩土力学, 2018, 39(1): 249-253. GAO Zirui, CHEN Tao, XU Yongfu. Effect of salt solution on swelling characteristics of bentonite[J]. Rock and Soil Mechanics, 2018, 39(1): 249-253. (in Chinese)
[20]	陈永贵, 李昆鹏, 马婧, 等. 化学作用下高庙子膨润土屏障性能演化行为[J]. 工程地质学报, 2022, 30(1): 71-82. CHEN Yonggui, LI Kunpeng, MA Jing, et al. Evolution of barrier properties for gaomiaozi bentonite under chemical effects[J]. Journal of Engineering Geology, 2022, 30(1): 71-82. (in Chinese)
[21]	马婧, 陈永贵, 刘聪, 等. 化学作用下压实膨润土膨胀力响应机制研究进展[J]. 岩土工程学报, 2023, 45(10): 2042-2051. doi: 10.11779/CJGE20220911 MA Jing, CHEN Yonggui, LIU Cong, et al. Research progress on the swelling pressures mechanisms of compacted bentonite under chemical conditions[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(10): 2042-2051. (in Chinese) doi: 10.11779/CJGE20220911
[22]	MITCHELL J K. Fundamentals of Soil Behavior[M]. 2nd ed. New York: Wiley, 1993.
[23]	GARCÍA-ROMERO E, LORENZO A, GARCÍA-ⅥCENTE A, et al. On the structural formula of smectites: a review and new data on the influence of exchangeable cations[J]. Journal of Applied Crystallography, 2021, 54: 251-262. doi: 10.1107/S1600576720016040
[24]	李昆鹏, 陈永贵, 叶为民, 等. 高压实膨润土孔隙结构特征研究进展[J]. 岩土工程学报, 2022, 44(3): 399-408. doi: 10.11779/CJGE202203001 LI Kunpeng, CHEN Yonggui, YE Weimin, et al. Advances in studies on pore structure of highly compacted bentonite[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(3): 399-408. (in Chinese) doi: 10.11779/CJGE202203001
[25]	SAIYOURI N, HICHER P Y, TESSIER D. Microstructural approach and transfer water modelling in highly compacted unsaturated swelling clays[J]. Mechanics of Cohesive-Frictional Materials, 2000, 5(1): 41-60. doi: 10.1002/(SICI)1099-1484(200001)5:1<41::AID-CFM75>3.0.CO;2-N
[26]	SAIYOURI N, TESSIER D, HICHER P Y. Experimental study of swelling in unsaturated compacted clays[J]. Clay Minerals, 2004, 39(4): 469-479. doi: 10.1180/0009855043940148
[27]	WANG Q. Hydro-mechanical behaviour of bentonite-basedmaterials used for high-level radioactive waste disposal[D]. Paris: Ecole Des Ponts Paris Tech, 2012.
[28]	KOMINE H, YASUHARA K, MURAKAMI S. Swelling characteristics of bentonites in artificial seawater[J]. Canadian Geotechnical Journal, 2009, 46(2): 177-189. doi: 10.1139/T08-120
[29]	李亚楠, 苏锐, 周志超, 等. 北山新场BS34钻孔岩样在不同温度下的水岩作用[J]. 核化学与放射化学, 2022, 44(3): 386-392. LI Yanan, SU Rui, ZHOU Zhichao, et al. Water-rock interaction of granite from borehole BS34 at different temperatures[J]. Journal of Nuclear and Radiochemistry, 2022, 44(3): 386-392. (in Chinese)
[30]	LIU L N, CHEN Y G, YE W M, et al. Effects of hyperalkaline solutions on the swelling pressure of compacted Gaomiaozi (GMZ) bentonite from the viewpoint of Na+ cations and OH–anions[J]. Applied Clay Science, 2018, 161: 334-342. doi: 10.1016/j.clay.2018.04.023
[31]	梁维云, 韦昌富, 张芹, 等. 膨润土吸湿过程中膨胀力演化及水分分布特征[J]. 岩土工程学报, 2023, 45(2): 283-291. doi: 10.11779/CJGE20211496 LIANG Weiyun, WEI Changfu, ZHANG Qin, et al. Swelling pressure evolution and water distribution characteristics of bentonite during wetting process[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(2): 283-291. (in Chinese) doi: 10.11779/CJGE20211496
[32]	叶为民, 刘樟荣, 崔玉军, 等. 膨润土膨胀力时程曲线的形态特征及其模拟[J]. 岩土工程学报, 2020, 42(1): 29-36. doi: 10.11779/CJGE202001003 YE Weimin, LIU Zhangrong, CUI Yujun, et al. Features and modelling of time-evolution curves of swelling pressure of bentonite[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(1): 29-36. (in Chinese) doi: 10.11779/CJGE202001003
[33]	PONS C H, ROUSSEAUX F, TCHOUBAR D. Utilisation du rayonnement synchrotron en diffusion Aux petits angles pour l'étude du gonflement des smectites: Ⅰ. Etude du systeme eau-montmorillonite-Na en fonction de la temperature[J]. Clay Minerals, 1981, 16(1): 23-42. doi: 10.1180/claymin.1981.016.1.02
[34]	SUZUKI S, PRAYONGPHAN S, ICHIKAWA Y, et al. In situ observations of the swelling of bentonite aggregates in NaCl solution[J]. Applied Clay Science, 2005, 29(2): 89-98. doi: 10.1016/j.clay.2004.11.001
[35]	CHEN Y G, DONG X X, ZHANG X D, et al. Combined thermal and saline effects on the swelling pressure of densely compacted GMZ bentonite[J]. Applied Clay Science, 2018, 166: 318-326. doi: 10.1016/j.clay.2018.10.001
[36]	SUN D A, ZHANG L, LI J, et al. Evaluation and prediction of the swelling pressures of GMZ bentonites saturated with saline solution[J]. Applied Clay Science, 2015, 105/106: 207-216.
[37]	KARNLAND O, OLSSON S, NILSSON U. Mineralogy and Sealing Properties of Various Bentonites and Smectite-Rich Clay Materials[R]. Stockholm: Svensk Kärnbränslehantering Ab, 2006.
[38]	LEE J O, LIM J G, KANG I M, et al. Swelling pressures of compacted Ca-bentonite[J]. Engineering Geology, 2012, 129/130: 20-26.
[39]	JADDA K, BAG R. Variation of swelling pressure, consolidation characteristics and hydraulic conductivity of two Indian bentonites due to electrolyte concentration[J]. Engineering Geology, 2020, 272: 105637.
[40]	THYAGARAJ T, RAO S M. Osmotic swelling and osmotic consolidation behaviour of compacted expansive clay[J]. Geotechnical and Geological Engineering, 2013, 31(2): 435-445.
[41]	DI MAIO C, SANTOLI L, SCHIAVONE P. Volume change behaviour of clays: the influence of mineral composition, pore fluid composition and stress state[J]. Mechanics of Materials, 2004, 36(5/6): 435-451.
[42]	项国圣, 徐永福, 陈涛, 等. 盐溶液中膨润土膨胀变形的分形模型[J]. 岩土力学, 2017, 38(1): 75-80. XIANG Guosheng, XU Yongfu, CHEN Tao, et al. Fractal model for swelling deformation of bentonite in salt solution[J]. Rock and Soil Mechanics, 2017, 38(1): 75-80. (in Chinese)
[43]	李彩霞, 李俊, 徐猛, 等. 氯盐溶液对钠基膨润土垫层膨胀性能的影响[J]. 土木与环境工程学报(中英文), 2023, 45(1): 97-104. LI Caixia, LI Jun, XU Meng, et al. Influence of chlorine salt solution on swelling properties of sodium bentonite cushion[J]. Journal of Civil and Environmental Engineering, 2023, 45(1): 97-104. (in Chinese)
[44]	李晓月, 徐永福. 盐溶液中膨润土膨胀变形的计算方法[J]. 岩土工程学报, 2019, 41(12): 2353-2359. doi: 10.11779/CJGE201912022 LI Xiaoyue, XU Yongfu. Method for calculating swelling deformation of bentonite in salt solution[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(12): 2353-2359. (in Chinese) doi: 10.11779/CJGE201912022
[45]	TRIPATHY S, BAG R, THOMAS H R. Effect of Stern-layer on the compressibility behaviour of bentonites[J]. Acta Geotechnica, 2014, 9(6): 1097-1109.
[46]	ZHU C M, YE W M, CHEN Y G, et al. Influence of salt solutions on the swelling pressure and hydraulic conductivity of compacted GMZ01 bentonite[J]. Engineering Geology, 2013, 166: 74-80.
[47]	WU T, WANG Z F, WANG H, et al. Salt effects on Re(Ⅶ) and Se(Ⅳ) diffusion in bentonite[J]. Applied Clay Science, 2017, 141: 104-110.
[48]	JIN Q, SU L, MONTAVON G, et al. Surface complexation modeling of U(Ⅵ) adsorption on granite at ambient/elevated temperature: experimental and XPS study[J]. Chemical Geology, 2016, 433: 81-91.
[49]	LI S C, WANG X L, HUANG Z Y, et al. Sorption–desorption hysteresis of uranium(Ⅵ) on/from GMZ bentonite[J]. Journal of Radioanalytical and Nuclear Chemistry, 2016, 310(2): 671-678.
[50]	刘福强, 叶远虑, 郭宁, 等. Eu(Ⅲ)在Na基高庙子膨润土上的吸附作用: 实验和构模研究[J]. 中国科学: 化学, 2013, 43(2): 242-252. LIU Fuqiang, YE Yuanlv, GUO Ning, et al. The adsorption of Eu(Ⅲ) on Gaomiaozi Na-bentonite: experimental and modeling study[J]. Scientia Sinica (Chimica), 2013, 43(2): 242-252. (in Chinese)
[51]	LU S S, XU H, WANG M M, et al. Sorption of Eu(Ⅲ) onto Gaomiaozi bentonite by batch technique as a function of pH, ionic strength, and humic acid[J]. Journal of Radioanalytical and Nuclear Chemistry, 2012, 292(2): 889-895.
[52]	CHEN Y G, ZHU C M, SUN Y H, et al. Adsorption of La(Ⅲ) onto GMZ bentonite: effect of contact time, bentonite content, pH value and ionic strength[J]. Journal of Radioanalytical and Nuclear Chemistry, 2012, 292(3): 1339-1347.
[53]	SUN Z, CHEN Y G, MU X, et al. Graphene oxide-modified organic Gaomiaozi bentonite for Yb(Ⅲ) adsorption from aqueous solutions[J]. Materials Chemistry and Physics, 2021, 274: 125176.
[54]	LI S C, WANG X L, HUANG Z Y, et al. Sorption and desorption of uranium(Ⅵ) on GMZ bentonite: effect of pH, ionic strength, foreign ions and humic substances[J]. Journal of Radioanalytical and Nuclear Chemistry, 2016, 308(3): 877-886.
[55]	DONG Y H, LIU Z J, LI Y Y. Effect of pH, ionic strength, foreign ions and humic substances on Th(Ⅳ) sorption to GMZ bentonite studied by batch experiments[J]. Journal of Radioanalytical and Nuclear Chemistry, 2011, 289(1): 257-265.
[56]	BRADBURY M H, BAEYENS B. Sorption of Eu on Na- and Ca-montmorillonites: experimental investigations and modelling with cation exchange and surface complexation[J]. Geochimica et Cosmochimica Acta, 2002, 66(13): 2325-2334.
[57]	COVELO E F, VEGA F A, ANDRADE M L. Competitive sorption and desorption of heavy metals by individual soil components[J]. Journal of Hazardous Materials, 2007, 140(1/2): 308-315.
[58]	GOMES P C, FONTES M P F, DA SILVA A G, et al. Selectivity sequence and competitive adsorption of heavy metals by Brazilian soils[J]. Soil Science Society of America Journal, 2001, 65(4): 1115-1121.
[59]	LIU C, XU Q W, XU Y W, et al. Characterization of adsorption behaviors of U(Ⅵ) on bentonite colloids: batch experiments, kinetic evaluation and thermodynamic analysis[J]. Journal of Radioanalytical and Nuclear Chemistry, 2022, 331(1): 597-607.
[60]	杜作勇, 王彦惠, 李东瑞, 等. 膨润土对U(Ⅵ)的吸附机理研究[J]. 核技术, 2019, 42(2): 22-29. DU Zuoyong, WANG Yanhui, LI Dongrui, et al. Adsorption mechanism of U(Ⅵ) by bentonite[J]. Nuclear Techniques, 2019, 42(2): 22-29. (in Chinese)

施引文献

资源附件(1)

其他相关附件
- PDF格式
  01-202412-G23-0811一文审稿意见与作者答复点击下载(1820KB)

图(21) / 表(2)

计量

文章访问数: 605
HTML全文浏览量: 66
PDF下载量: 234
被引次数: 0

0. 引言
1. 管-土-高聚物相互作用分析模型
1.1 管-土-高聚物相互作用Winkler模型
1.2 管-土-高聚物相互作用Vlazov模型
2. 基于传递矩阵法的管道平面应变问题分析
2.1 传递矩阵法
2.2 管道计算模型描述
2.3 程序实现
3. 试验概述及数据处理
4. 模型计算及与试验结果对比
5. 不同因素对修复后管道内力和变形的影响
5.1 不同埋深的影响
5.2 不同附加荷载的影响
5.3 不同土壤弹性模量的影响
5.4 不同管壁厚度的影响
6. 结论

0. 引言
1. 管-土-高聚物相互作用分析模型
1.1 管-土-高聚物相互作用Winkler模型
1.2 管-土-高聚物相互作用Vlazov模型
2. 基于传递矩阵法的管道平面应变问题分析
2.1 传递矩阵法
2.2 管道计算模型描述
2.3 程序实现
3. 试验概述及数据处理
4. 模型计算及与试验结果对比
5. 不同因素对修复后管道内力和变形的影响
5.1 不同埋深的影响
5.2 不同附加荷载的影响
5.3 不同土壤弹性模量的影响
5.4 不同管壁厚度的影响
6. 结论

参考文献(60)

施引文献

资源附件(1)

高压实膨润土与孔隙溶液物理作用机制研究进展 English Version

作者简介: 廖饶平(1993—)，男，博士研究生，主要从事环境地质与非饱和土力学研究。E-mail: lrp_liao@tongji.edu.cn

通讯作者: 陈永贵, E-mail: cyg@tongji.edu.cn

计量

出版历程

Advances in physical interaction mechanism between highly compacted bentonite and pore solution

0. 引言

1. 管-土-高聚物相互作用分析模型

1.1 管-土-高聚物相互作用Winkler模型

1.2 管-土-高聚物相互作用Vlazov模型

2. 基于传递矩阵法的管道平面应变问题分析

2.1 传递矩阵法

2.2 管道计算模型描述

2.3 程序实现

3. 试验概述及数据处理

4. 模型计算及与试验结果对比

5. 不同因素对修复后管道内力和变形的影响

5.1 不同埋深的影响

5.2 不同附加荷载的影响

5.3 不同土壤弹性模量的影响

5.4 不同管壁厚度的影响

6. 结论

其他相关附件

PDF格式

计量

出版历程

目录

0. 引言

1. 管-土-高聚物相互作用分析模型

1.1 管-土-高聚物相互作用Winkler模型

1.2 管-土-高聚物相互作用Vlazov模型

2. 基于传递矩阵法的管道平面应变问题分析

2.1 传递矩阵法

2.2 管道计算模型描述

2.3 程序实现

3. 试验概述及数据处理

4. 模型计算及与试验结果对比

5. 不同因素对修复后管道内力和变形的影响

5.1 不同埋深的影响

5.2 不同附加荷载的影响

5.3 不同土壤弹性模量的影响

5.4 不同管壁厚度的影响

6. 结论

作者简介:
廖饶平(1993—)，男，博士研究生，主要从事环境地质与非饱和土力学研究。E-mail: lrp_liao@tongji.edu.cn

通讯作者:
陈永贵, E-mail: cyg@tongji.edu.cn