非饱和盐渍黏土物理化学作用的应力依赖特性

王立业; 周凤玺; 周立增; 梁玉旺

doi:10.11779/CJGE20221559

非饱和盐渍黏土物理化学作用的应力依赖特性 English Version

王立业^1,,
周凤玺^{1, 2, ,},
周立增³,
梁玉旺¹

1.
兰州理工大学土木工程学院, 甘肃兰州 730050
2.
西部土木工程防灾减灾教育部工程研究中心, 甘肃兰州 730050
3.
西北工业大学力学与土木建筑学院, 陕西西安 710072

基金项目:

国家自然科学基金项目 11962016

国家自然科学基金项目 51978320

甘肃省基础研究创新群体项目 20JR5RA478

甘肃省教育厅：优秀研究生“创新之星”项目 2022CXZX-451

详细信息

作者简介:
王立业(1993—)，男，博士研究生，主要从事盐渍土试验和理论模型方面的研究。E-mail: gwly1024@163.com

通讯作者:
周凤玺, E-mail: geolut@163.com

中图分类号: TU43
计量
- 文章访问数: 347
- HTML全文浏览量: 39
- PDF下载量: 83
出版历程
- 收稿日期: 2022-12-27
- 网络出版日期: 2024-04-09
- 刊出日期: 2024-03-31

Stress-dependent properties of physicochemical interaction of unsaturated saline clay

1.
School of Civil Engineering, Lanzhou University of Technology, Lanzhou 730050, China
2.
Engineering Research Center of Disaster Mitigation in Civil Engineering of Ministry of education, Lanzhou 730050, China
3.
School of Mechanics, Civil Engineering and Architecture, Northwestern Polytechnical University, Xi'an 710072, China

摘要

摘要: 为了揭示物理化学作用对含盐溶液非饱和黏土压缩行为的影响及其与应力水平的依赖特性，对孔隙含蒸馏水、氯化钠和硫酸钠溶液及控制基质吸力条件下的试样进行了一维压缩试验。然后，根据试验结果测定了不同条件下黏土的压缩指数、次固结系数和屈服应力，并标定了它们随基质吸力和渗透吸力的变化规律。此外，对不同基质吸力和渗透吸力下非饱和盐渍黏土主、次固结行为展开了深入分析，明确了非饱和盐渍黏土物理化学作用的应力依赖特征，并对非饱和盐渍黏土的LC屈服行为进行了探索。结果表明：不同物理化学力下非饱和盐渍黏土的次固结系数与压缩指数的比值 ${C_{\text{a}}}/{C_{\text{c}}}$ 以及屈服应力可以使用基质吸力和渗透吸力来统一描述；由塑性加载区压缩曲线斜率以及特征参量、 ${C_{\text{a}}}/{C_{\text{c}}}$ 与竖向应力的相关性指出非饱和盐渍黏土物理化学作用与应力水平息息相关；此外，化学-水力-力学耦合作用下非饱和盐渍黏土的LC屈服曲线是由MLC屈服曲线和OLC屈服曲线组成的一条光滑曲线。
- 非饱和黏土 /
- 孔隙盐溶液 /
- 物理化学作用 /
- 应力依赖性 /
- LC屈服行为
Abstract: To reveal the impact of physicochemical effects on the compressive behaviors of unsaturated clay containing salt solution and its dependence properties on stress level, one-dimensional compression tests are performed on the specimens with pores containing distilled water, sodium chloride solution, sodium sulfate solution and controlled matric suction conditions. Then, the compression index, secondary compression coefficient and yield stress of clay under different conditions are measured according to the test results, and their variation laws with the matric suction and osmotic suction are calibrated. Furthermore, the stress-dependent characteristics of physicochemical action are clarified through an in-depth analysis of the primary and secondary consolidation behaviors of unsaturated saline clay under different matric suctions and osmotic suctions, and the LC yielding behaviors of unsaturated saline clay are explored. The results show that the ratio of the secondary compression coefficient to the compression index ${C_{\text{a}}}/{C_{\text{c}}}$ and the yield stress of unsaturated saline clay at different physicochemical forces can be described uniformly using the osmotic suction and matric suction. From the slope of compression curve in the plastic loading zone and the correlation between characteristic parameter, ${C_{\text{a}}}/{C_{\text{c}}}$ and the vertical stress, it is noted that the physicochemical action of unsaturated saline clay is closely related to the stress level. Moreover, the LC yield curve of unsaturated saline clay is a smooth curve composed of MLC yield curve and OLC yield curve under chemo- hydro-mechanical coupling.
- unsaturated clay /
- pore salt solution /
- physicochemical action /
- stress dependence /
- LC yield behavior

HTML全文

0. 引言

根据国家“十四五”规划中对能源体系的发展计划，至2025年非化石能源使用占比提升至20%，需建设一批风电、光电、水电等清洁能源基地，同时加快抽水蓄能电站和新型储能技术规模化应用^[1]。压缩空气储能技术（compress air energy storage，简称CAES）是新型的大型能源存储技术，采用压缩空气作为工质对能量进行存储和释放，具备改善电力品质、提高电力系统稳定可靠性以及效率高，布置灵活、投资和运行费用较少等优点，一直是国内外研究热点^[2]。

1969年世界上第一个商业大型CAES电站在德国Huntorf诞生^[3]，20世纪70年代美国能源部开始对CAES电站进行示范性工程计划，该计划关注的两个主要问题为：①CAES长期运行过程储层稳定性标准；②尽量少用或不用化石燃料补充燃烧的CAES系统可行性^[4]。直至目前CAES的研究的出发点依旧是紧密围绕CAES电站的安全可靠性和经济性这两个关键问题，针对经济性的研究主要围绕压缩热量处理^[5-6]，压缩空气热力学过程^[7-10]、能量计算^[11-12]，系统设计及能量利用效率计算^[13-14]等。针对安全可靠性的研究主要围绕地下储气库稳定性这一问题展开，其研究成果远少于经济性研究成果。

CAES电站地下储气库深处地下，本身受力影响因素众多，受力特性复杂，再加上空气在储气库中频繁循环压缩和释放，造成压力和温度载荷的频繁循环变化，使长期运行过程中地下储气库可靠性和稳定性的精准预测和评估成为难点。目前多数学者的研究成果均存在较大局限性，如Carranza-torres等^[15]采用解析解和有限差分两种方法研究了浅埋圆形洞室的地质稳定性；Kim等^[16]采用极限平衡法探究了储气库设计参数对地面隆起的影响；夏才初等^[17]采用数值模拟方法分析了衬砌围岩受力特性。然而以上研究过程中并未对温度这一关键物理因素对储气库受力影响进行考虑。蒋中明等^[18]、周舒威等^[19-20]、夏才初等^[21]对温度的影响进行考虑，采用热力耦合理论分别采用FLAC3D、解析解和COMSOL研究了储气库衬砌围岩的受力特性，但忽略了运行过程中气体渗漏对稳定性的影响。周瑜等^[22]则考虑了储气库的渗漏问题，但忽略了温度的作用。目前仅有Kim等^[23]、Rutqvist等^[24]采用TOUGH-FLAC对储气库衬砌围岩进行了热流固耦合（thermo-hydro-mechanical，简称THM）研究。

以上研究成果虽有进步性，然而对多物理场耦合过程考虑的还不够全面，尤其针对有无地下水两种情况下储气库气密性和稳定性的研究仍需进一步深入。为探究地下水对CAES硬岩地下储气库运行过程中的气密性和稳定性影响，本文基于压缩空气热力学、热流固耦合理论和两相渗流理论，分别建立了非稳态边界条件下衬砌围岩单相和两相THM耦合数值模型，对比分析了地下水对CAES电站运行过程中储气库稳定性的影响，为此类工程的选址、设计、建设和运行提供科学依据。

1. 孔隙介质多物理场耦合理论和数值实现方法

当CAES电站选址位置无地下水或者地下水位较低时，运行过程中储气库不受地下水影响，考虑气体渗漏过程和压缩空气的温度变化和洞壁传热过程时，储气库衬砌围岩的受力特性可采用单相THM耦合理论进行描述。

1.1 单相THM耦合理论

基于线性热弹性假设，在考虑渗流耦合作用和温度应力作用下时，总应变为应力导致的应变、渗流压力导致的应变和热应变之和，固体总应变控制方程^[25]表示如下：

。

${\sigma _{ij}}{\text{ = }}2G{\varepsilon _{ij}} + \lambda {\delta _{ij}}{\varepsilon _{\text{V}}} - {\beta _{\text{T}}}{\delta _{ij}}({T_{\text{s}}} - {T_{{\text{s0}}}}) + \alpha {\delta _{ij}}(p - {p_0})。$

(1)

式中： $G{\text{ = }}E/2(1 + \nu )$ ； $\lambda {\text{ = }}2\nu G/(1 - 2\nu )$ ； ${\varepsilon _{ij}} = ({u_{i,j}} + {u_{j,i}}{\text{)/2}}$ ； ${\beta _{\text{T}}} = {\alpha _{\text{T}}}E/(1 - 2\nu )$ 。其中： ${\sigma _{ij}}$ 为应力二阶张量分量（Pa）；G为岩石剪切模量（Pa）； ${\varepsilon _{ij}}$ 为应变张量的分量；E为岩石弹性模量（Pa）； $\nu$ 为泊松比； $\lambda$ 为拉梅系数（Pa）； ${\delta _{ij}}$ 为Kronecker符号，当i=j时，δ_ij=1，当i≠j时， ${\delta _{ij}}$ =0； ${\varepsilon _{\text{v}}}$ 为体积应变； ${\beta _{\text{T}}}$ 为热膨胀因子； ${\alpha _{\text{T}}}$ 为热膨胀系数（K^-1）；T_s为固体温度（K）；T_s0为固体初始温度（K）；α为有效应力系数，无量纲；p为孔隙渗流压力（Pa）； ${p_0}$ 为初始孔隙压力（Pa）；u_i和u_j为分别为位移在i和j方向上的分量。

不考虑衬砌结构的配筋和损伤破坏，将其视作理想弹性材料，围岩屈服准则采用莫尔库仑模型：

。

${\sigma _1} - {\sigma _3} = 2c\cos \varphi - ({\sigma _1} + {\sigma _3})\sin \varphi 。$

(2)

式中： ${\sigma _1}$ ， ${\sigma _3}$ 分别为第一、三主应力（N/m²）；c为黏聚力（MPa）；φ为内摩擦角（°）。

孔隙内流体流动满足达西定理，同时考虑固体流体应力变形和温度膨胀过程，热流固耦合作用下渗流方程如下：

。

$\frac{\theta }{{{\beta _{\text{l}}}}}\frac{{\partial p}}{{\partial t}} = \frac{K}{{{\mu _{\text{g}}}}}(\nabla p + {\rho _{\text{f}}}g) - \alpha \frac{{\partial {\varepsilon _{\text{V}}}}}{{\partial t}} + (\theta {\alpha _{\text{l}}} + (1 - \theta ){\alpha _{\text{s}}})\frac{{\partial T}}{{\partial t}} 。$

(3)

式中：θ为孔隙介质孔隙率； ${\beta _{\text{l}}}$ 为流体体积模量（Pa）；K为多孔基质渗透率（m²）； ${\mu _{\text{g}}}$ 为气体动力黏度（Pa·s）； ${\rho _{\text{f}}}$ 为流体密度（kg/m³）；g为重力加速度（m/s²）； ${\alpha _{\text{l}}}$ 和 ${\alpha _{\text{s}}}$ 分为流体和固体组分热膨胀系数（K^-1）；T为温度（K）。

其中考虑体积变形和热膨胀作用时多孔介质孔隙率和渗透率的变化如下所示^[26]：

$\theta {\text{ = }}1 - \frac{{(1 - {\theta _0})}}{{1 + {\varepsilon _{\text{v}}}}}(1 - \frac{{\Delta p}}{{{K_{\text{s}}}}} + {\alpha _{\text{s}}}\Delta T)\text{，}$

(4)

$K{\text{ = }}{K_0}{{\text{(}}{\theta \mathord{\left/ {\vphantom {\theta {{\theta _0}}}} \right. } {{\theta _0}}}{\text{)}}^3}。$

(5)

气体在多孔介质中的流动，在温度差异下产生对流传热和传导传热过程，此处采用局部热平衡假设，假设多孔基质和孔隙内气体温度在极短时间内达到热平衡状态，此时的多孔基质内传热控制方程如下所示：

$\begin{array}{l} ({\theta _{\text{s}}}{\rho _{\text{s}}}{C_{{\text{sp}}}} + (1 - {\theta _{\text{s}}})\rho {C_{\text{p}}})\frac{{\partial T}}{{\partial t}} + {\theta _{\text{s}}}K'{\alpha _{\text{T}}}{T_0}\frac{{\partial {\varepsilon _{\text{V}}}}}{{\partial t}} + \\ \rho {C_{\text{p}}}{u_{\text{g}}} \cdot \nabla T - \nabla \cdot ({\theta _{\text{s}}}{k_{\text{s}}} + (1 - {\theta _{\text{s}}})k)\nabla T = Q。 \end{array}$

(6)

式中： ${\theta _{\text{s}}}$ 为多孔介质内基质体积分数； ${\rho _{\text{s}}}$ 为基质材料密度（kg/m³）；C_sp为基质材料恒压热容（J/(kg·K)）；C_p为气体恒压热容（J/(kg·K)）； $\nabla T$ 为温度梯度（K/m）；k_s为固体基质导热率（W/(m·K)）；k为气体导热率（W/(m·K)）；Q为热源（J/(m³ s)）。

1.2 两相THM耦合理论

当储气库处于地下水影响范围内时，整个物理过程为两相热流固耦合过程，其中固体力学理论和耦合理论如上一节基本相同，但式（1）中的孔隙压力为气相和液相的平均压力；气体渗漏过程由单相渗流变为气液两相渗流，根据达西定律和毛管力驱动可得其控制方程如下：

$\left.\begin{array}{l}\frac{\theta }{{\beta }_{\text{l}}}\frac{\partial p}{\partial t}+\nabla \cdot \left[-K\left({s}_{1}\frac{{k}_{\text{r}}{}_{1}}{{\mu }_{1}}+{s}_{2}\frac{{k}_{\text{r}}{}_{2}}{{\mu }_{2}}\right)(\nabla p+({s}_{1}{\rho }_{1}+{s}_{2}{\rho }_{2})g)\right]\\ \text{ }=-\alpha \frac{\partial {\varepsilon }_{\mathrm{V}}}{\partial t}+(\theta {\alpha }_{\text{l}}+(1-\theta ){\alpha }_{\text{s}})\frac{\partial T}{\partial t}\text{ }\text{，}\\ \frac{\partial (\theta {c}_{1})}{\partial t}+\nabla \cdot \left[-K\left({s}_{1}\frac{{k}_{\text{r}}{}_{1}}{{\mu }_{1}}+{s}_{2}\frac{{k}_{\text{r}}{}_{2}}{{\mu }_{2}}\right)\nabla p\cdot {c}_{1}\right]\\ \text{ }=\nabla \cdot \left[\frac{{k}_{\text{r1}}}{{\mu }_{1}}k({s}_{1}-1)\frac{\partial {p}_{\text{ec}}}{\partial {s}_{1}}\cdot \nabla {c}_{1}\right]\text{ }。\end{array}\right\}$

(7)

其中， ${\beta _{\text{l}}} = {s_1} \times {\beta _{\text{w}}} + {s_2} \times {\beta _{\text{g}}}$ ， ${c_1} = {s_1}{\rho _1}{k_{{\text{r1}}}} = {s_1}^{{l_{{\text{vg}}}}}(1 -$ ${(1 - {s_1}^{{1 \mathord{\left/ {\vphantom {1 {{m_{{\text{vg}}}}}}} \right. } {{m_{{\text{vg}}}}}}})^{{m_{{\text{vg}}}}}}{)^2}$ ， ${k_{{\text{r2}}}} = {(1 - {s_1})^{{l_{{\text{vg}}}}}}{(1 - {s_1}^{1/{m_{{\text{vg}}}}})^{2{m_{{\text{vg}}}}}}$ 。

式中：s₁和s₂分别为湿润相和非湿润相体积分数； ${\beta _{\text{w}}}$ 和β_g分别为液体和气体体积模量（Pa）； ${\rho _1}$ 和 ${\rho _2}$ 分别为湿润相和非湿润相流体密度（kg/m³）； ${k_{{\text{r1}}}}$ 和 ${k_{{\text{r2}}}}$ 分别为湿润相和非湿润相流体相对渗透率；l_vg和m_vg分别为模型参数；c₁为湿润相流体在孔隙中含量（kg/m³）。

采用局部热平衡假设，且根据体积平均理论，多孔基质内两相流体传热控制方程如下所示：

$\begin{array}{l} ({\theta _{\text{s}}}{\rho _{\text{s}}}{C_{{\text{sp}}}} + (1 - {\theta _{\text{s}}})\rho C{p_{{\text{ef}}}})\frac{{\partial T}}{{\partial t}} + {\theta _{\text{s}}}K'{\alpha _{\text{T}}}{T_0}\frac{{\partial {\varepsilon _{\text{V}}}}}{{\partial t}} + \\ \rho {C_{{\text{pf}}}}u\nabla T - \nabla ({\theta _{\text{s}}}{k_{\text{s}}} + (1 - {\theta _{\text{s}}}){k_{{\text{ef}}}})\nabla T = Q。 \end{array}$

(8)

式中： ${\theta _{\text{s}}}$ 为多孔介质内基质体积分数； $C{p_{{\text{ef}}}}$ 为流体有效比热容（J/(kg·K)）， $C{p_{{\text{ef}}}}$ =s₁·Cp₁+ s₂·Cp₂，Cp₁和Cp₁分为流体相1和2的比热容；k_ef为流体有效导热率（W/(m·K)），k_ef= s₁·k₁+ s₂·k₂，k₁和k₂分别为流体相1和2的导热系数。

1.3 储气库THM耦合动态边界模型

关于CAES电站运行过程中储气库内压缩空气热力学模型，He等^[27]做出了相当有益的工作，本部分内容在其研究结果上进行了细化推导和演化，增加了渗漏项。考虑洞壁传热和气体渗漏过程的压缩空气热力学模型如下：

$\left.\begin{array}{l}\text{d}m/\text{d}t={m}_{\text{in}}-{m}_{\text{out}}-{m}_{\text{e}}\text{ }\text{，}\\ {m}_{\text{e}}={\displaystyle \oint u{\rho }_{\text{g}}D\text{d}}s\text{ }\text{，}\\ \text{d}T/\text{d}t=(\dot{Q}+{m}_{\text{in}}Cp({T}_{\text{in}}-T)+V\text{d}p/\text{d}t)/mCp\text{ }\text{，}\\ \text{d}p/\text{d}t=1/V[\kappa R{m}_{\text{in}}{T}_{\text{in}}-\kappa \text{R}{m}_{\text{out}}T-\kappa \text{R}{m}_{\text{e}}T+\\ (\kappa -1){h}_{\text{c}}A({T}_{\text{w}}-T)]\text{ }。\end{array}\right\}$

(9)

式中：m为储气库内压缩空气质量（kg）； ${m_{{\text{in}}}}$ 和 ${m_{{\text{out}}}}$ 分别为充放气速率（kg/s）； ${m_{\text{e}}}$ 为储气库整体渗流渗漏速率（kg/s）；D为圆柱形储气库长度（m）；s为圆柱形储气库横断面周长（m）；V为储气库容积（m³）； ${T_{{\text{in}}}}$ 为充气温度（K）；R为气体常数，取值287.06 J/(K kg)；к为绝热指数； ${h_{\text{c}}}$ 为储气库洞壁对流换热系数，（W/(m²·K)）；A为洞壁的对流换热面积（m²）； ${T_{\text{w}}}$ 为储气库洞壁温度（K）。

数值实现时，储气库压缩空气实时变化的压力作为应力和压力边界作用在硐室内壁边界上；压缩空气实时变化的温度则作为热传导分析的边界作用在洞壁表面。

1.4 THM耦合数值分析流程

由于渗流速度u和储气库洞壁温度 ${T_{\text{w}}}$ 受到压缩空气压力p和温度T影响，以COMSOL数学全局ODE（ordinary differential equation）来建立压缩空气热力学模型，将耦合时变渗漏和洞壁传热的空气温度T压力p作为储气库洞室内壁非稳态边界。流场采用弱形式PDE（partial differential equation）在达西定律与两相达西模块进行二次开发建模。耦合数值分析模型计算流程如图 1所示。

图 1 储气库衬砌围岩两相热流固耦合模型计算流程

Figure 1. Calculation flow of two-phase THM coupling model for lining and surrounding rock of gas reservoir

下载: 全尺寸图片幻灯片

1.5 模型及程序验证

为验证热流固模型的正确性，本文采用经典的一维饱和土柱弹性热固结问题进行验证，该问题前后经过Noorisha等^[28]、孙致学等^[25]学者的验证。利用COMSOL建立该热流固耦合数值模型并求解，模型高为7 m，宽1 m。模型初始孔隙水压力为10 kPa，模型初始温度为10℃，土柱顶部向下施加10 kPa均布载荷，同时受60℃温度作用，且为自由渗流边界。模型左右边界和底部的力学边界为法向约束，热学边界为绝热边界，渗流边界为不透水边界。模型及网格如图 2所示，模型参数如文献[25]。

图 2 热固结验证模型

Figure 2. Verification model of thermal consolidation

下载: 全尺寸图片幻灯片

热固结数值模型计算所得到的沉降位移、渗流压力、温度验证结果分别如图 3~5所示，由图可知，本模型与数值解所得出结果在数值和规律上一致性较好，位移曲线先下降，最后开始上升，其过程代表含义为孔隙水排除后，孔隙压力消散，有效应力增加，位移增加，当孔隙水排净后，温度传递导致的膨胀作用导致位移向上。从位移、压力和温度对比可表明该数值模型的正确性。

图 3 不同高度沉降位移验证结果

Figure 3. Verification results of displacement at different heights

下载: 全尺寸图片幻灯片

图 4 不同高度渗流压力验证结果

Figure 4. Verification results of seepage pressure at different heights

下载: 全尺寸图片幻灯片

图 5 不同高度温度验证结果

Figure 5. Verification results of temperature at different heights

下载: 全尺寸图片幻灯片

2. 算例研究

以某拟建地下硬岩CAES储气库为例，建立其横断面二维数值模型，图 6为地下储气库模型及边界设置，储气库为圆柱形横向布置，半径7.5 m，洞室长度300 m，衬砌密封层采用低渗混凝土，厚度为0.5 m，洞室埋深250 m，围岩为花岗岩。假定地下水位在地平面位置。模型采用二维轴对称模型进行建模，为尽量减弱边界尺寸对计算结果的影响，围岩研究范围定为500 m。

图 6 CAES储气库数值模型及边界示意图

Figure 6. Numerical model and boundary diagram of CAES gas reservoir

下载: 全尺寸图片幻灯片

电站运行以1 d为一个周期，分别为充气（8 h）—高压储气（4 h）—放气（4 h）—低压储气（8 h）。入气温度为15℃，储气库的运行压力上下限为7~10 MPa。充放气速率调整结果为初始阶段以50.3 kg/s的充气速率充气24 h，正式运行时充气速率为47.5 kg/s，放气速率为95 kg/s。

依照Kim等^[23]、Rutqvist等^[24]提出的地下水两相THM模型计算流程，模型计算流程分4步：①初始自重应力平衡过程计算，清除初始自重位移。②地下储气库的开挖过程计算，计算开挖排水之后7 d的储气库位移变形和围岩内孔压变化。③计算衬砌施工后1 d内的衬砌围岩的变形位移和孔压变化。④计算包括初始充气阶段在内的60个运行循环过程中的压缩空气热力学过程和储气库两相THM耦合过程。

模型左侧为对称边界，对力学模型底部为固定约束、右边为法向约束、顶部为自由面；对渗流模型底部、右侧和顶部均为随埋深变化的常压和水相饱和度边界，对传热模型底部、右侧和顶部为常温边界。计算相关参数如表 1所示。

表 1 计算参数取值表

Table 1. Values of calculation parameters

参数	值/单位	参数	值/单位	参数	值/单位	参数	值/单位
储气库半径r	7.5 m	两相流参数l_vg	0.5	围岩初始孔隙率θ_r0	0.02	衬砌恒压热容Cp_c	880 J/(kgK)
储气库长度D	300 m	两相流参数m_vg	0.595	围岩恒压热容Cp_r	930 J/(kg K)	衬砌导热率λ_c	2 W/(m·K)
衬砌密封层厚度d	0.5 m	围岩弹性模量E_r	20 GPa	围岩导热率λ_r	1.8 W/(m·K)	围岩热膨胀系数	5.0×10^-51/K
初始压力P₀	1 atm	衬砌弹性模量E_c	23 GPa	衬砌热膨胀系数	1.2×10^-5 1/K	围岩内摩擦角	50 °
初始温度T₀	15 ℃	围岩密度ρ_r	2650 kg/m³	衬砌初始孔隙率θ_c0	0.01	围岩泊松比ν_r	0.18
充气温度T__in	15 ℃	衬砌密度ρ_c	2300 kg/m³	衬砌初始渗透率K_c0	1×10^-20 m²	衬砌泊松比ν_c	0.2
毛细入口压力p_ec	1.47 MPa			围岩初始渗透率K_r0	2×10^-16 m²	围岩黏聚力c	2 MPa

下载: 导出CSV

| 显示表格

3. 计算成果及分析

3.1 储气库压缩空气热力学过程变化

压缩空气热力学特性变化过程如图 7所示，由图可知，第一个充气阶段，压缩空气温度和压力随时间迅速升高，空气压力随时间线性升高，在24 h迅速上升至10 MPa，温度先迅速升高后缓慢升高至415 K。压缩空气的温度和压力在放气过程中迅速降低，随后的运行循环中，空气温度会波动降低，稳定时的温度波动上下限279~314 K，波动范围35 K。

图 7 压缩空气热力学特性变化过程

Figure 7. Variation of thermodynamic characteristics of compressed air

下载: 全尺寸图片幻灯片

3.2 地下水对储气库密封性影响

密封性以储气库内气体质量为对象，储气库渗漏率以储气库气体减少为负，增加为正。两种模型储气库60个运行循环整体渗漏速率变化过程如图 8所示，有无地下水两种地质条件下的储气库渗漏速率对运行过程周期性波动变化。由于空气压力的动态变化，初始时衬砌和围岩内渗流压力较低，但是储气库内气体压力迅速升高，因此初始循环时气体渗漏最快。随着渗流压力扩散，衬砌和围岩内渗流压力升高，此时渗漏速率逐步降低，最后形成稳定波动。单相THM模型渗漏速率稳定时波动范围为-0.034~0.0165 kg/s，两相THM模型渗漏速率稳定波动范围为-0.0135~0.011 kg/s。值得注意的是，在放气阶段，气体渗漏为正，即气体在放气阶段由衬砌围岩向储气库内渗流，原因为放气阶段由于气体压力迅速降低3 MPa，衬砌结构内渗流压力的变化存在滞后性，其压力高于储气库内气体压力，此时渗流方向向内。

图 8 两种模型储气库整体渗漏速率变化规律

Figure 8. Variation of overall leakage rate of two models for gas reservoir

下载: 全尺寸图片幻灯片

储气库每个循环气体渗漏质量变化规律如图 9所示，每个循环气体渗漏质量变化规律随循环过程总体成减小并稳定的趋势。单相THM模型在第4个循环是渗漏质量达到稳定，稳定时每个循环气体质量损失为440 kg，两相THM模型在10个循环时达到稳定，稳定时每个循环气体损失为95 kg。该储气库稳定运行时充气速率47.5 kg/s，充气时间8 h，每个循环充气质量为1.37×10⁶ kg，单相THM模型每个循环气体渗漏质量占充气质量0.032%，两相THM模型每个循环气体渗漏质量占充气质量0.007%。两个模型稳定运行时每个循环气体损失均小于1%，表明当衬砌密封层渗透率为1×10^-20 m²，考虑其受力过程产生的应变而导致渗透率变化，该衬砌密封层密封效果满足要求。且地下水存在时，会导致气体质量损失降低78%，地下水可有效辅助密封。

图 9 储气库每个循环气体渗漏质量随运行过程变化规律

Figure 9. Variation of gas leakage quality in each cycle of gas reservoir with operation process

下载: 全尺寸图片幻灯片

3.3 地下水对储气库结构受力特性影响

储气库衬砌围岩结构的受力状态决定了储气库的安全稳定性，两种模型在稳定运行过程充气段末衬砌和围岩有效应力分布如图 10所示，由图可知，当衬砌结构强度高于围岩时，衬砌结构受力高于围岩，是主要受力结构，衬砌顶部和底部会出现应力集中现象，单相THM模型充气段末衬砌最大有效应力为27.3 MPa，而两相THM模型衬砌最大有效应力为15.7 MPa。当储气库地层中有地下水作用时，衬砌结构顶部和顶部的集中有效应力会有较大幅度降低，充气段末的最大应力降低幅度为42.5%。

图 10 充气段末衬砌和围岩受力分布云图

Figure 10. Stress distribution of lining and surrounding rock at end of charging stage

下载: 全尺寸图片幻灯片

分析衬砌密封层结构的主要受力形态和可能发生的破坏形式，有助于对储气库结构的安全稳定性做出准确评估。此处第三、第一主应力按照岩土力学规定进行划分。储气库衬砌结构第60循环充气和放气段末第三、第一主应力分布如图 11所示，由图 11可知，①放气段末衬砌结构第三主应力高于充气段末，即放气段结构失稳危险程度高于充气段；②衬砌结构在充气和放气阶段第三主应力大范围表现为正，即衬砌结构大范围承受拉应力，拉应力较高位置为顶部和底部，侧壁位置第三主应力为负，表现为压应力；③两相THM模型衬砌结构充放气段第三、第一主应力均小于单相THM模型，即有地下水时，衬砌结构最大拉应力可以降低52%。

图 11 典型运行阶段末两种模型衬砌第三和第一主应力分布

Figure 11. Distribution of first and third principal stresses of lining structures of two models at end of typical operation stage

下载: 全尺寸图片幻灯片

为消除开挖和施工变形位移差异对衬砌结构最终变形位移结果的影响，以运行初期为参照对两种模型衬砌结构典型位置测点位移变化进行研究，其规律如图 12所示。衬砌径向位移随空气压力升高迅速变大，然后随充放气状态起伏波动，稳定运行时径向位移波动范围为1.5 mm。洞室侧壁水平位移大于顶部竖向位移，即运行过程中储气库横向膨胀变形大于纵向，整体变形成椭圆形膨胀。单相THM模型侧壁水平位移和顶部竖向位移分别比两相THM模型大0.8 mm和0.4 mm。

图 12 衬砌测点径向位移变化规律

Figure 12. Variation of radial displacement of monitoring points for lining

下载: 全尺寸图片幻灯片

衬砌典型测点环向应变变化规律如图 13所示，衬砌结构顶部应变大于侧壁，该规律与Kim等^[29]研究成果相一致。衬砌顶部环向应变为拉应变，侧壁环向应变状态随充放气变化，充气时为拉应变，放气时为压应变。单相THM模型衬砌测点环向应变均大于两相THM模型，单相THM模型衬砌最大环向应变为0.0011，该应变结果与相似工况下的夏才初等^[17]的研究成果相类似，而两相THM模型最大为0.0006。可知地下水可有效降低衬砌顶部拉应变。

图 13 衬砌测点环向应变变化规律

Figure 13. Variation of circumferential strain at different measuring points for lining

下载: 全尺寸图片幻灯片

储气库运行过程中围岩塑性区面积变化和分布规律如图 14所示。储气库围岩塑性区为剪切塑性区，稳定运行时两相THM模型围岩塑性区成环形分布于洞室周围，而单相THM模型围岩塑性区成矩形分布于洞室周围，在洞室顶部和底部向围岩中延伸，最远可达围岩内部20 m。开挖围岩塑性区集中于洞室腰部围岩内，随着运行过程增加，洞室顶部和底部围岩有效应力增加，塑性区逐渐向顶部和底部围岩中发展。初始运行时单相THM模型塑性区面积为332 m²，随后的两个循环中急速增长，之后缓慢增长至第10个循环时达到稳定，稳定时围岩塑性区面积为1337 m²。初始运行时两相THM模型塑性区面积为248 m²，在第3个循环时达到稳定，稳定时塑性区面积为356 m²。地下水可以有效降低运行过程中围岩塑性区分布。根据此项结果，围岩在施工过程中需要进行加固处理，在没有地下水影响时需对洞室顶部、底部深部区域进行加固处理。

图 14 围岩塑性区面积随时间变化及分布规律

Figure 14. Variation and distribution laws of plastic zone area of surrounding rock with time

下载: 全尺寸图片幻灯片

3.4 讨论

储气库THM耦合模型结果和两相THM耦合模型结果对比如表 2所示。相同条件下当有地下水时，储气库渗漏速率从-0.034~0.016 kg/s降低至-0.013~0.011 kg/s，每个循环渗漏百分比从0.035%降低至0.008%，地下水的存在可有效辅助气体密封。

表 2 有无地下水对储气库渗漏和受力特性影响

Table 2. Influences of groundwater on leakage and mechanical.characteristics of gas reservoir

工况	渗漏速率/(kg·s^-1)	渗漏百分比%	衬砌最大有效应力/MPa	衬砌最大第三主应力/MPa	围岩最大有效应力/MPa	衬砌最大拉应变	围岩塑性区面积/m²
无地下水	-0.034 ~ 0.016	0.032	27.3	29.6	7.63	0.0011	1337
有地下水	-0.013 ~ 0.011	0.007	15.7	13.3	7.57	0.0006	356
注：表中应力符号规定为拉正压负，渗漏速率以气体从储气库向外为负

下载: 导出CSV

| 显示表格

衬砌出现应力集中，当有地下水时，衬砌最大有效应力从27.3 MPa降低至15.7 MPa，最大第三主应力也从29.6 MPa降低至13.3 MPa，最大剪应力从31 MPa降低至17.9 MPa，衬砌最大拉应变降低0.0005。围岩剪应力大小和分布均出现较大范围降低，塑性区面积从1337 m²降低至356 m²。

地下水存在时，衬砌外侧和围岩内孔隙压力增高，有效应力降低，尤其是衬砌结构内，衬砌最大有效应力、最大第三主应力和剪应力均降低1/3。地下水的影响会使围岩塑性区面积也出现了大幅度降低，降低幅度为92%，其原因为有地下水时，洞周围岩内孔压增加，有效应力降低、同时剪应力降低，由于正应力降低导致的抗剪强度也随之降低，但是在本文模型的洞室形状、洞室埋深、地下水位和运行条件下，围岩的剪应力降低幅度高于抗剪强度的降低幅度，因此围岩剪切塑性区分布区域减小。有效应力和塑性区的降低幅度主要由孔隙压力增高幅度决定，孔隙压力的增高幅度主要由埋深与地下水位相对关系控制。由此可知，地下水的存在会大幅度降低衬砌有效应力和减小围岩塑性区面积。

4. 结论

对CAES储气库密封可靠性和安全稳定性的研究是保证CAES电站可靠运行的基础，地下储气库的密封效果和安全稳定性的影响因素众多，耦合过程复杂，同时耦合变量之间高度非线性增加了对储气库安全性和可靠性评价的难度。本文采用热流固耦合理论、两相渗流理论、压缩空气热力学理论对储气库运行过程中涉及到的流体状态变化、空气渗漏、衬砌围岩单相THM和两相THM耦合过程进行了建模。该模型计算了有无地下水两种情况下储气库从开挖至60 d运行过程，主要得出以下3点结论。

（1）在该储气库以本文模型设定方案运行时，衬砌渗透率为1×10^-20 m²时，储气库渗漏率低于1%，可满足密封要求。地下水的存在可使储气库每个循环气体损失量降低78%。

（2）该模型储气库运行时衬砌为受力主体，衬砌顶部和底部出现应力集中，且表现为受拉应力作用，因此衬砌顶部和底部需注意配加密筋以满足受力要求。运行时储气库总体成椭圆形膨胀变形，横向变形大于竖向。

（3）地下水存在时可大幅度降低衬砌结构有效应力和围岩塑性区面积。

图 1 非饱和盐渍黏土的一维压缩试验结果

Figure 1. One-dimensional compression test results of unsaturated saline clay

下载: 全尺寸图片幻灯片

图 2 非饱和盐渍黏土在压缩过程中体积应变的变化情况

Figure 2. Change in volume strain of unsaturated saline clay during compression process

下载: 全尺寸图片幻灯片

图 3 不同含水和含盐条件下试样的吸力

Figure 3. Suctions of samples under different water and salt contents

下载: 全尺寸图片幻灯片

图 4 非饱和盐渍黏土物理化学作用应力依赖性的示意图

Figure 4. Schematic of stress dependence of physicochemical interaction of unsaturated saline clay

下载: 全尺寸图片幻灯片

图 5 压缩指数随基质吸力和渗透吸力的变化规律

Figure 5. Variation of compression index with matric suction and osmotic suction

下载: 全尺寸图片幻灯片

图 6 非饱和盐渍黏土物理化学作用的应力依赖性分析

Figure 6. Analysis of stress dependence of physicochemical interaction of unsaturated saline clay

下载: 全尺寸图片幻灯片

不同应力水平下 ${C_{\text{a}}}/{C_{\text{c}}}$ 与渗透吸力和基质吸力的关系

Relationship among ${C_{\text{a}}}/{C_{\text{c}}}$ , osmotic suction and matric suction at different stress levels

下载: 全尺寸图片幻灯片

图 8 拟合参数随竖向应力的变化规律

Figure 8. Variation of fitting parameters with vertical stress

下载: 全尺寸图片幻灯片

图 9 屈服应力随基质吸力和渗透吸力的变化规律

Figure 9. Variation of yield stress with matric suction and osmotic suction

下载: 全尺寸图片幻灯片

图 10 非饱和盐渍黏土的LC屈服曲线

Figure 10. LC yield curves of unsaturated saline clay

下载: 全尺寸图片幻灯片

表 1 试验用土中易溶盐的离子种类和含量

Table 1 Ion species and contents of soluble salts in experimental soil

阳离子/%					阴离子/%			总含盐量/%	pH值
${{\text{K}}^ + }$	${\text{N}}{{\text{a}}^ + }$	${\text{M}}{{\text{g}}^{2 + }}$	${\text{G}}{{\text{a}}^{2 + }}$	${\text{CO}}_3^{2 - }$	${\text{HCO}}_3^ -$	${\text{SO}}_4^{2 - }$	${\text{C}}{{\text{l}}^ - }$	总含盐量/%	pH值
0.0113	0.4470	0.0481	0.0546	—	0.0235	0.3222	0.0912	1.0	7.47

下载: 导出CSV

表 2 不同基质吸力和渗透吸力下的次固结系数

Table 2 Secondary consolidation coefficients under different conditions of matric suction and osmotic suction

s = 0 kPa		s = 50 kPa						s = 100 kPa		s = 200 kPa
π = 457 kPa		π = 457 kPa		π = 5630 kPa		π = 5902 kPa		π = 457 kPa		π = 457 kPa		π = 5630 kPa		π = 5902 kPa
σ/kPa	C_a×10²	σ/kPa	C_a×10²	σ/kPa	C_a×10²	σ/kPa	C_a×10²	σ/kPa	C_a×10²	σ/kPa	C_a×10²	σ/kPa	C_a×10²	σ/kPa	C_a×10²
100	0.1165	100	0.1165	100	0.0764	97	0.0872	101	0.1156	97	0.1161	99	0.0792	97	0.1126
200	0.1685	200	0.1222	200	0.1650	201	0.1676	194	0.1160	194	0.1141	200	0.1451	201	0.1540
300	0.2186	300	0.1318	300	0.1555	291	0.1657	291	0.1220	291	0.1140	300	0.1669	291	0.1673
400	0.2184	401	0.1755	400	0.2141	388	0.2339	387	0.1625	388	0.1159	400	0.1256	387	0.1675
600	0.2672	601	0.1959	600	0.2644	581	0.2626	581	0.1739	581	0.1645	599	0.2098	581	0.2201
800	0.2911	802	0.2645	800	0.3076	775	0.3529	775	0.2040	775	0.1744	800	0.2511	774	0.2722
1000	0.3623	1003	0.3075	1001	0.3952	969	0.4083	969	0.2915	969	0.2510	1000	0.3028	968	0.3217

下载: 导出CSV

参考文献(34)

[1]	张莎莎, 谢山杰, 杨晓华, 等. 火山灰改良粗粒硫酸盐渍土路基填料及其作用机理研究[J]. 岩土工程学报, 2019, 41(3): 588-594. doi: 10.11779/CJGE201903023 ZHANG Shasha, XIE Shanjie, YANG Xiaohua, et al. Action mechanism of coarse particle sulfate soil subgrade modified by volcanic ash[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(3): 588-594. (in Chinese) doi: 10.11779/CJGE201903023
[2]	ZHANG Y, YE W M, CHEN B, et al. Desiccation of NaCl-contaminated soil of earthen heritages in the Site of Yar City, northwest China[J]. Applied Clay Science, 2016, 124-125: 1-10. doi: 10.1016/j.clay.2016.01.047
[3]	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.
[4]	YAO CQ, CHEN P, MA TT, et al. Physicochemical effect on shear strength characteristics of clayey soils based on ring-shear experiment[J]. Canadian Geotechnical Journal, 2020, 57(12): 1820-1831. doi: 10.1139/cgj-2019-0513
[5]	陈云敏. 环境土工基本理论及工程应用[J]. 岩土工程学报, 2014, 36(1): 1-46. doi: 10.11779/CJGE201401001 CHEN Yunmin. A fundamental theory of environmental geotechnics and its application[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(1): 1-46. (in Chinese) doi: 10.11779/CJGE201401001
[6]	THYAGARAJ T, DAS A P. Physico-chemical effects on collapse behaviour of compacted red soil[J]. Géotechnique, 2017, 67(7): 559-571. doi: 10.1680/jgeot.15.P.240
[7]	孙重初. 酸液对红粘土物理力学性质的影响[J]. 岩土工程学报, 1989, 11(4): 89-93. http://cge.nhri.cn/cn/article/id/9296 SUN Zhongchu. The effect of acidizing fluid on the physico mechanical properties of red clay[J]. Chinese Journal of Geotechnical Engineering, 1989, 11(4): 89-93. (in Chinese) http://cge.nhri.cn/cn/article/id/9296
[8]	陈声凯, 凌建明, 罗志刚. 路基土回弹模量应力依赖性分析及预估模型[J]. 土木工程学报, 2007, 40(6): 95-99, 104. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC200706018.htm CHEN Shengkai, LING Jianming, LUO Zhigang. Stress-dependent characteristics and prediction model of the resilient modulus of subgrade soils[J]. China Civil Engineering Journal, 2007, 40(6): 95-99, 104. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC200706018.htm
[9]	弋晓明, 李术才, 王松根, 等. 非饱和粉土回弹模量的应力依赖性与水敏感性耦合分析[J]. 山东大学学报(工学版), 2013, 43(2): 84-88, 110. https://www.cnki.com.cn/Article/CJFDTOTAL-SDGY201302015.htm YI Xiaoming, LI Shucai, WANG Songgen, et al. Coupling analysis of stress dependence and water sensitivity for the resilient modulus of unsaturated silt soil[J]. Journal of Shandong University (Engineering Science), 2013, 43(2): 84-88, 110. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SDGY201302015.htm
[10]	SUN D A, SHENG D C, CUI H B, et al. A density-dependent elastoplastic hydro-mechanical model for unsaturated compacted soils[J]. International Journal for Numerical & Analytical Methods in Geomechanics, 2010, 31(11): 1257-1279.
[11]	XU X T, JIAN W B, WU N S, et al. Void ratio-dependent water retention model for a deformable residual clay[J]. International Journal of Geomechanics, 2020, 20(8): 04020131-1. doi: 10.1061/(ASCE)GM.1943-5622.0001773
[12]	CUI Y J, DELAGE P. Yielding and plastic behaviour of an unsaturated compacted silt[J]. Géotechnique, 1996, 46(2): 291-311. doi: 10.1680/geot.1996.46.2.291
[13]	RAMPINO C, MANCUSO C, VINALE F. Experimental behaviour and modelling of an unsaturated compacted soil[J]. Canadian Geotechnical Journal, 2000, 37(4): 748-763. doi: 10.1139/t00-004
[14]	高登辉, 陈正汉, 郭楠, 等. 干密度和基质吸力对重塑非饱和黄土变形与强度特性的影响[J]. 岩石力学与工程学报, 2017, 36(3): 736-744. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201703023.htm GUO Denghui, CHEN Zhnghan, HUO Nan, et al. The influence of dry density and matric suction on the deformation and the strength characteristics of the remolded unsaturated loess soils[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(3): 736-744. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201703023.htm
[15]	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.
[16]	张彤炜, 邓永锋, 刘松玉, 等. 渗透吸力对重塑黏土的压缩和渗透特性影响的试验研究[J]. 岩土工程学报, 2014, 36(12): 2260-2266. doi: 10.11779/CJGE201412014 ZHANG Tongwei, DENG Yongfeng, LIU Songyu, et al. Experimental investigation of osmotic suction effect on hydro-mechanical behaviour of remolded clay[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(12): 2260-2266. (in Chinese) doi: 10.11779/CJGE201412014
[17]	WITTEVEEN P, FERRARI A, LALOUI L. An experimental and constitutive investigation on the chemo-mechanical behaviour of a clay[J]. Géotechnique, 2013, 63(3): 244-255. doi: 10.1680/geot.SIP13.P.027
[18]	颜荣涛, 赵续月, 于明波, 等. 盐溶液饱和黏土的等向压缩特性[J]. 岩土力学, 2018, 39(1): 129-138. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201801017.htm YAN Rongtao, ZHAO Xuyue, YU Mingbo, et al. Isotropic compression characteristics of clayey soil saturated by salty solution[J]. Rock and Soil Mechanics, 2018, 39(1): 129-138. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201801017.htm
[19]	周凤玺, 王立业, 赖远明. 饱和盐渍土的一维蠕变试验与模型研究[J]. 岩土工程学报, 2020, 42(1): 142-149. doi: 10.11779/CJGE202001016 ZHOU Fengxi, WANG Liye, LAI Yuanming. One-dimensional creep tests and model studies on saturated saline soil[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(1): 142-149. (in Chinese) doi: 10.11779/CJGE202001016
[20]	王立业, 周凤玺, 秦虎. 饱和盐渍土分数阶蠕变模型及试验研究[J]. 岩土力学, 2020, 41(2): 543-551. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202002026.htm WANG Liye, ZHOU Fengxi, QIN Hu. Fractional creep model and experimental study of saturated saline soil[J]. Rock and Soil Mechanics, 2020, 41(2): 543-551. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202002026.htm
[21]	ALONSO E E, GENS A, JOSA A. A constitutive model for partially saturated soils[J]. Géotechnique, 1990, 40(3): 405-430. doi: 10.1680/geot.1990.40.3.405
[22]	陈正汉, 扈胜霞, 孙树国, 等. 非饱和土固结仪和直剪仪的研制及应用[J]. 岩土工程学报, 2004, 26(2): 161-166. http://cge.nhri.cn/cn/article/id/11362 CHEN Zhenghan, HU Shengxia, SUN Shuguo, et al. Development & application of consolidation apparatus and direct shear apparatus for unsaturated soils[J]. Chinese Journal of Geotechnical Engineering, 2004, 26(2): 161-166. (in Chinese) http://cge.nhri.cn/cn/article/id/11362
[23]	陈正汉. 重塑非饱和黄土的变形、强度、屈服和水量变化特性[J]. 岩土工程学报, 1999, 21(1): 82-90. http://cge.nhri.cn/cn/article/id/10260 CHEN Zhenghan. Deformation, strength, yield and moisture change of a remolded unsaturated loess[J]. Chinese Journal of Geotechnical Engineering, 1999, 21(1): 82-90. (in Chinese) http://cge.nhri.cn/cn/article/id/10260
[24]	THYAGARAJ T, THOMAS S R, DAS A P. Physico -chemical effects on shrinkage behavior of compacted expansive clay[J]. International Journal of Geomechanics, 2017, 17(2): 06016013. doi: 10.1061/(ASCE)GM.1943-5622.0000698
[25]	LEONG E, RAHARDJO H, HE L. Factors affecting the filter paper method for total and matric suction measurements[J]. Geotechnical Testing Journal, 2002, 25(3): 8198.
[26]	SHENG D C. Review of fundamental principles in modelling unsaturated soil behaviour[J]. Computers and Geotechnics, 2011, 38(6): 757-776. doi: 10.1016/j.compgeo.2011.05.002
[27]	RAMPINO C, MANCUSO C, VINALE F. Laboratory testing on an unsaturated soil: equipment, procedures, and first experimental results[J]. Canadian Geotechnical Journal, 1999, 36(1): 1-12. doi: 10.1139/t98-093
[28]	LU N, LIKOS W J. Suction stress characteristic curve for unsaturated soil[J]. Journal of Geotechnical & Geoenvironmental Engineering, 2006, 132(2): 131-142.
[29]	周凤玺, 王立业, 赖远明. 饱和盐渍土渗透吸力的回顾及研究[J]. 岩土工程学报, 2020, 42(7): 1199-1210. doi: 10.11779/CJGE202007003 ZHOU Fengxi, WANG Liye, LAI Yuanming. Review and research on osmotic suction of saturated saline soils[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(7): 1199-1210. (in Chinese) doi: 10.11779/CJGE202007003
[30]	PHILIP J R. Unitary approach to capillary condensation and adsorption[J]. Journal of Chemical Physics, 1977, 66(11): 5069-5075. doi: 10.1063/1.433814
[31]	LU N, ZHANG C. Soil Sorptive Potential: Concept, Theory, and Verification[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2018, 145(4): 04019006.
[32]	赵建斌, 申俊敏, 张艳聪. 考虑土体应力依赖性的水泥混凝土路面基层力学响应研究[J]. 公路, 2014, 59(12): 54-58. https://www.cnki.com.cn/Article/CJFDTOTAL-GLGL201412012.htm ZHAO Jianbin, SHEN Junmin, ZHANG Yancong. Research on mechanical responses of the base for cement concrete pavement considering stress dependence of subgrade[J]. Highway, 2014, 59(12): 54-58. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GLGL201412012.htm
[33]	KONG L M, YAO Y P, QI J L. Modeling the combined effect of time and temperature on normally consolidated and overconsolidated clays[J]. Acta Geotechnica, 2020, 15(3): 2451-2471.
[34]	MESRI G, GODLEWSKI P M. Time and stress-compressibility interrelationship[J]. Journal of the Geotechnical Engineering Division, 1977, 103(5): 417-430. doi: 10.1061/AJGEB6.0000421