复合非饱和土中水-热-盐多场耦合分析

高国耀; 郭伟; 周凤玺

doi:10.11779/CJGE2022S1010

复合非饱和土中水-热-盐多场耦合分析 English Version

高国耀^1,,
郭伟^1, ,,
周凤玺²

1.
天津大学建筑工程学院, 天津 300000
2.
兰州理工大学土木工程学院, 甘肃兰州 730000

基金项目:

国家自然科学基金项目 51878446

国家自然科学基金项目 52171273

国家自然科学基金项目 52108335

详细信息

作者简介:
高国耀(1994—)，男，博士研究生，主要从事非饱和土、冻土的多场耦合研究。E-mail: Guoyao-Gao@163.com

通讯作者:
郭伟, E-mail: guow@tju.edu.cn

中图分类号: TU443
计量
- 文章访问数: 211
- HTML全文浏览量: 34
- PDF下载量: 39
出版历程
- 收稿日期: 2022-09-27
- 网络出版日期: 2023-02-06
- 刊出日期: 2022-11-30

Multi-field coupling analysis of water-heat-salt in composite unsaturated soils

1.
School of Civil Engineering, Tanjin University, Tanjin 300000, China
2.
School of Civil Engineering, Lanzhou University of Technology, Lanzhou 730000, China

摘要

摘要: 研究了非饱和黏土上覆非饱和砂土组成的复合非饱和土中的水-热-盐多场耦合问题。基于多孔介质的质量守恒和能量守恒方程，选取温度、孔隙压力和含盐率以及它们的梯度作为状态变量，建立了非稳态条件下水热盐多场耦合问题的状态方程组。将状态方程组进行Laplace变换后，进而求解该强耦合的非线性变系数微分方程组的频域解，并通过反演法获得时域解。并将数值解与自主设计试验的试验结果进行了对比，验证了所建立数学模型的准确性。分析了复合非饱和土体在温度梯度下，底部非饱和砂土的厚度和饱和渗透系数等对整个非饱和土体中水分场和盐分场分布以及迁移规律的影响，继而对盐渍土的改良工作和工程土障的建设提供试验和理论依据。
- 复合非饱和土 /
- 传热-传质 /
- Laplace变换 /
- 反演法
Abstract: The multi-field coupling problem of water, heat and salt in the composite unsaturated soils composed of unsaturated clay and unsaturated sand is studied. Based on the mass and energy conservation equations for porous media, the temperature, pore pressure, salt content and their gradients are selected as the state variables, and the state equations for the multi-field coupling problem of water, heat and salt under unsteady conditions are established. After the Laplace transformation of the state equations, the frequency-domain solution for the strongly coupled nonlinear variable coefficient differential equations is solved, and the time-domain solution is obtained by the inversion method. The numerical solution is compared with the experimental results of the self-designed tests, and the accuracy of the established mathematical model is verified. The influences of the thickness and saturated permeability coefficient of the unsaturated sand at the bottom of the composite unsaturated soils on the distribution and migration law of water and salt fields in the whole unsaturated soils under the temperature gradient are analyzed, and then the test and theoretical basis are provided for the improvement of saline soils and the construction of engineering soil barriers.
- composite unsaturated soil /
- heat and mass transfer /
- Laplace transform /
- inversion method

HTML全文

0. 引言

随着中国基础设施建设的快速发展，在非饱和盐渍土地区进行盐渍土的改良以及工程土障建设中，均涉及到水-热-气-传质等多物理场多相的耦合作用。Reshetin等^[1]基于扩散定律和温度、湿度和蒸汽含量之间的函数关系，建立了关于温度场与湿度场的动力学数学模型，但数学模型中没有考虑非饱和土中水力梯度变化对耦合作用的影响。Jia等^[2]在建立传热传质耦合模型时，考虑到水蒸汽对传热传质效应的影响。基于压力梯度和温度梯度等多种驱动因素建立了多种有关非饱和多孔介质的水-热耦合的数学模型，但忽略了非饱和多孔介质中水力梯度随时间和空间的变化规律且对气相的组分划分较为笼统^[3-6]，Mohammad等^[7]考虑干燥气体与水蒸气的迁移对能量守恒的影响，但在在能量迁移的方式中忽略了对流作用影响。He等^[8]在非饱和土水-热-力多场耦合的数学模型，没有考虑多场多相耦合作用中溶质的存在对多孔介质中水分场、温度场迁移规律的影响。

多孔介质中水-热-盐迁移现象是各影响因素相互耦合的过程，且迁移规律随时间和空间不断变化。Cleall等^[9]，HERNÁNDEZ-LÓPEZ等^[10]通过试验孔隙的曲率因子和土体温度梯度的增强因子对水蒸汽通量有较大影响，Chen等^[11]发现随着盐分在一端的大量聚集，溶液中的盐分会在盐分梯度作用下继续向该处迁移。周凤玺等^[12-13]建立了完善热-湿-盐多物理场多相耦合的数学模型。但模型仅考虑均质非饱和多孔介质中的热湿盐多场耦合，并没有对复合非饱和土中的水-热-盐迁移规律进行进一步的研究且缺乏试验验证。

本文全面研究了复合非饱和土中的液态水、液相中的干燥气体和盐分以及孔隙中水蒸气等对水分场、盐分场和温度场随时间变化在空间上的分布规律的影响，并在能量守恒中考虑盐分、干燥气体以及各组分导热系数等分项对温度场的影响，建立了非稳态条件下复合非饱和土中多场多相耦合的状态方程。通过试验分单层验证了数学模型的准确性而后对复合非饱和土中水-盐-热的多场耦合问题进行参数分析。

1. 数学模型

1.1 质量守恒方程

组分α在π相中的质量守恒方程表示为^[12-13]

$\frac{\partial ({\phi }^{\pi }{\rho }_{\alpha }^{\pi })}{\partial t}+\nabla \cdot ({\phi }^{\pi }{\rho }_{\alpha }^{\pi }{\boldsymbol{v}}^{\alpha })+\nabla \cdot {\boldsymbol{j}}_{\alpha }^{\pi }={\dot{m}}_{\alpha }^{\pi }\text{，} \text{，}$

(1)

式中， ${\boldsymbol{j}}_\alpha ^\pi$ 为α组分在π相的非对流通量矢量， ${\phi ^\pi }$ 为π相的体积分数， $\rho _\alpha ^\pi$ 为α组分在π相的质量分数， ${{\boldsymbol{v}}^\alpha }$ 为α组分的速度矢量^[12-13]。

由式（1）可得水分平衡方程为

$\frac{{\partial (n{S^{\text{l}}}\rho _{\text{w}}^{\text{l}})}}{{\partial t}} + \frac{{\partial (n{S^{\text{g}}}\rho _{\text{w}}^{\text{g}})}}{{\partial t}} + \nabla\cdot {\text{(}}n{S^{\text{g}}}\rho _{\text{w}}^{\text{g}}{{\boldsymbol{v}}^{\text{g}}}{\text{)}} + \nabla\cdot {\boldsymbol{j}}_{\text{w}}^{\text{g}} + \\ \;\;\;\;\;\nabla \cdot (n{S}^{\text{l}}{\rho }_{\text{w}}^{\text{l}}{\boldsymbol{v}}^{\text{l}})+\nabla \cdot ({\rho }_{\text{w}}^{\text{l}}{q}^{\text{l}})+\nabla \cdot ({\rho }_{\text{w}}^{\text{g}}{q}^{\text{g}})=0\text{ }\text{，}$

(2)

式中，饱和度 ${S^{\text{l}}}$ 为孔隙水压力P_l与孔隙气压力P_g的函数，存在表达式P_l + P_g = 1。液相、气相通量^[12-13]存在表达式：

${q^{\bf{l}}} =- {K^{\bf{l}}}(\nabla {P_{\bf{l}}} + \rho _{\text{w}}^{\bf{l}}g) \text{，}$

(3)

${q^\text{g}} =- {K^\text{g}}(\nabla {P_\text{g}} + {\rho ^\text{g}}g) \text{，}$

(4)

式中，K^l，K^g为水动力参数^[12-13]，g为重力加速度矢量。

非对流水通量遵循Philip等^[14]所建立的模型，气相中水蒸汽非对流通量的变化可以归因于系统中水分的变化和系统中温度的变化：

${\boldsymbol{j}}_{\text{w}}^{\text{g}} = {\boldsymbol{j}}_{{\text{vw}}}^{\text{g}} + {\boldsymbol{j}}_{{\text{vT}}}^{\text{g}} =- {D_{{\text{vw}}}}\nabla {P_{\text{l}}} + {D_{{\text{vw}}}}\nabla {P_{\text{g}}} - {D_{{\text{vT}}}}\nabla T 。$

(5)

其中由温度和水分引起的分子扩散系数存在表达式^[12-13]：

${D_{{\text{vT}}}} = n{S^{\text{g}}}\rho _{\text{w}}^{\text{g}}{D_{{\text{atm}}}}v\tau {f_{{\text{Tv}}}}\left[ {\frac{{4974.0}}{{{T^2}}} + \frac{{{M_{\text{w}}}\psi }}{{R{T^2}\rho _{\text{w}}^{\text{l}}}}} \right] \text{，}$

(6)

${D_{{\text{vw}}}} = n{S^{\text{g}}}\rho _{\text{w}}^{\text{g}}{D_{{\text{atm}}}}v\tau \frac{{{M_{\text{w}}}}}{{RT\rho _{\text{w}}^{\text{l}}}} \text{，}$

(7)

式中，D_atm为水蒸汽分子的扩散系数，D_atm = 2.16×10^-5(T/273.15)^1.8，τ为曲率因子^[12-13]。

结合式（1）可得干燥气体的质量平衡方程为

$\frac{{\partial (nH{S^{\text{l}}}\rho _{\text{a}}^{\text{l}})}}{{\partial t}} + \nabla\cdot ({\phi ^{\text{l}}}\rho _{\text{a}}^{\text{l}}{{\boldsymbol{v}}^{\text{l}}}) + \nabla\cdot {\boldsymbol{j}}_{\text{a}}^{\text{l}} +\\ \;\;\;\;\;\; \frac{\partial (n{S}^{\text{g}}{\rho }_{\text{a}}^{\text{g}})}{\partial t}+\nabla \cdot ({\phi }^{\text{g}}{\rho }_{\text{a}}^{\text{g}}{v}^{\text{g}})+\nabla \cdot {j}_{\text{a}}^{\text{g}}=0\text{ }。$

(8)

干燥气体在液相中不发生扩散，即 ${\boldsymbol{j}}_{\text{a}}^{\text{l}} = 0$ 。非对流项中的干燥气体的非对流通量与水蒸气的非对流通量存在表达式 ${\boldsymbol{j}}_{\text{a}}^{\text{g}} =- {\boldsymbol{j}}_{\text{w}}^{\text{g}}$ ，水中气体的密度存在^[12-13]：

$\rho _{\text{a}}^{\text{l}} \approx \rho _{\text{a}}^{\text{g}} = {\rho _{\text{a}}} = \frac{{{M_{\text{a}}}}}{{RT}}({P_{\text{g}}} - {P_{\text{v}}}) \text{，}$

(9)

式中，P_v为水蒸气压力，M_a为干气摩尔质量，R为气体常数。

结合式（1）可以得出液相中盐分的质量守恒方程为

$\frac{{\partial ({\phi ^{\text{l}}}\rho _{\text{p}}^{\text{l}})}}{{\partial t}} + \nabla\cdot (\rho _{\text{w}}^{\text{l}}\omega {q^{\text{l}}}) + \nabla\cdot {\boldsymbol{j}}_{\text{p}}^{\text{l}} = 0 。$

(10)

由于研究土壤液相中盐分得浓度较低，不考虑溶质的溶解与析出。液相中溶质非对流通量 ${\boldsymbol{j}}_{\text{p}}^{\text{l}}$ 可以表示为

${\boldsymbol{j}}_{\text{p}}^{\text{l}} =- n{S^{\text{l}}}\rho _{\text{w}}^{\text{l}}\tau {D_{\text{p}}}\nabla \omega \text{，}$

(11)

式中，D_p为溶质的扩散系数^[12-13]。

1.2 能量平衡方程

系统中能量平衡的一般形式为^[15]

$\frac{{\partial {\mathit{\Phi} _\text{h}}}}{{\partial t}} + Q + \nabla\cdot {q_\text{h}} = 0 \text{，}$

(12)

式中，Φ_h为内能，Q为潜热， $\nabla\cdot {q_\text{h}}$ 表示热对流和传导。

考虑局部平衡以及不同组分温度相等的条件下，土壤的内能变化可表示为

$\frac{{\partial {\mathit{\Phi} _\text{h}}}}{{\partial t}} = \left[ {{\text{(}}1 - n{\text{)}}{\rho ^\text{s}}{c_\text{s}} + nH{S^\text{l}}{\rho _\text{a}}{c_\text{a}} + } \right. \\ \;\;\;\;\;\;\;\; n{S}^{\text{g}}{\rho }_{\text{a}}{c}_{\text{a}}+n{S}^{\text{l}}{\rho }_{\text{w}}^{\text{l}}{c}_{\text{w}}^{\text{l}}+n{S}^{\text{g}}{\rho }_{\text{w}}^{\text{g}}{c}_{\text{w}}^{\text{g}}]\frac{\partial T}{\partial t}\text{ }。$

(13)

在非饱和土中由汽化和凝结引起的能量变化可表示为

$Q = \dot m_\text{w}^\text{g}{H_{\text{gw}}} \text{，}$

(14)

式中， ${H_\text{gw}}$ 为液相中水分的汽化焓^[16]， ${H_{\text{gw}}} = {L_{\text{gw}}} -$ $(c_\text{w}^\text{l} - c_\text{w}^\text{g})$ $(T - {T_0})$ ， $\dot m_\text{w}^\text{g}$ 为蒸发率。

热流方程包括固相热传导、液相对流、气相对流，存在表达式^[14]：

$\begin{aligned} q_{\mathrm{h}}= & q^{\mathrm{T}}+\rho_\alpha^\pi E_{\mathrm{T} \alpha}^\pi q^\alpha+E_{\mathrm{T} \alpha}^\pi {\boldsymbol{j}}_\alpha^\pi \\ = & -\left(\lambda_{\mathrm{s}}^{1-n} \lambda_{\mathrm{w}}^{n S^l} \lambda_{\mathrm{p}}^{n S^{\mathrm{p}}}\right) \nabla T+\left(\rho_a c_a+\rho_{\mathrm{w}}^{\mathrm{g}} c_{\mathrm{w}}^{\mathrm{g}}\right)\left(T-T_0\right) q^{\mathrm{g}}+ \\ & \left(c_{\mathrm{w}}^{\mathrm{g}}-c_a\right)\left(T-T_0\right) {\boldsymbol{j}}_{\mathrm{w}}^{\mathrm{g}}+c_{\mathrm{p}}^1 \boldsymbol{j}_{\mathrm{p}}^{\mathrm{1}}\left(T-T_0\right)+ \\ & \left(\rho_a c_a H+\rho_{\mathrm{w}}^1 c_{\mathrm{w}}^1\right)\left(T-T_0\right) q^1。 \end{aligned}$

(15)

以往的研究忽略温度对溶液热容与导热系数的影响，溶液热容和导热系数与溶质的浓度有关^[17]

$\left.\begin{array}{l} c_{\mathrm{w}}^{\prime}=\beta_2 c_{\mathrm{s}}+c_{\mathrm{w}}, \\ \lambda_{\mathrm{w}}=\beta_3 c_{\mathrm{s}}+\lambda_{\mathrm{l}}, \end{array}\right\}$

(16)

式中， ${\beta _2}$ = -167， ${\beta _3}$ = 0.00493。

1.3 复合非饱和土非稳态问题求解

方程（2），（8），（10），（12）给出了复合非饱和土中水-气-盐-热多场耦合过程的控制方程，耦合方程是将孔隙水压力、孔隙气压力、温度、含盐率及其一阶偏导为未知量的封闭方程。将控制方程进行拉普拉斯变换后通过数值方法计算可得控制方程的频域解，而后基于Hausdoff矩问题将频域解 $u(x)$ 转化为时域解 ${q_N}(t)$ 。

Hausdoff矩问题存在积分关系式

${\displaystyle {\int }_{0}^{1}{x}^{n}u(x)\text{d}x={\mu }_{n}}\text{ (}n=0,1,2,\cdot \cdot \cdot \text{) }。$

(17)

Hausdorff矩问题是一个典型的不适定问题， $u(x)$ 的正则化近似解为

${p_N}(x) = \sum\limits_{i = 0}^N {{\lambda _i}{L_i}(x)} 。$

(18)

最终Laplace逆变换的近似解为

${q_N}{\text{(}}t{\text{)}} \equiv {p_N}({{\text{e}}^{ - t}}) = \sum\limits_{i = 0}^N {{\lambda _i}{L_i}({{\text{e}}^{ - t}})} 。$

(19)

基于式（19）Hausdoff矩问题的稳定化算法获得Laplace逆变换的解。

2. 试验验证

试验采用的土样是取自甘肃省兰州市的粉质黏土，该土样的颗粒相对质量密度为2.70 g/cm³，渗透系数约为2.89×10^-6 m/s。将土样压制成高为280 mm、直径为200 mm，试验土柱的孔隙率n为0.3。初始时刻试验土柱的温度为25℃，顶板与底板的温度分别设为30，15℃。两端的孔隙气压力为101 kPa、初始质量含水率为10 %、初始含盐率为7%。

表 1，2分别给出了试验粉质黏土土样各物理成分的热力学参数和Van Genuchten（VG）模型^[13，15]参数。经试验和计算可得粉质黏土土样在VG模型中所对应的拟合参数分别为n_vg=1.69，m_vg=0.0147。图 1~3分别给出了试样中水分，盐分和温度的分布规律，并与本文所提出的教学模型进行验证，两者吻合较好。

表 1 多孔介质中各物理成分的热力学参数

Table 1. Thermodynamic parameters of physical components in porous media

参数	数值	参数	数值
${c_{\text{a}}}$ /(J·kg^-1·K^-1)	1000	${L_{{\text{gw}}}}$ /(J·kg^-1)	2413000
$c_{\text{w}}^{\text{l}}$ /(J·kg^-1·K^-1)	4180	${k_0}$ /(m·s^-1)	2.89×10^-6
$c_{\text{w}}^{\text{g}}$ /(J·kg^-1·K^-1)	1900	g/(m·s^-2)	9.8
$c_{\text{p}}^{\text{l}}$ /(J·kg^-1·K^-1)	3750	$\rho _{\text{w}}^{\text{g}}$ /(kg^-1·m^-3)	1000
${\lambda _{\text{g}}}$ /(W·m^-1·K^-1)	0.024	${M_{\text{a}}}$ /(kg^-1·mol^-1）	0.0288
${\lambda _{\text{s}}}$ /(W·m^-1·K^-1)	2.93	${M_{\text{w}}}$ /(kg^-1·mol^-1)	0.018016
${\lambda _{\text{l}}}$ /(W·m^-1·K^-1)	0.56	$R$ /(J·mol^-1·K^-1)	8.3144

下载: 导出CSV

| 显示表格

表 2 粉质黏土VG模型参数

Table 2. VG model parameters of silty clay

参数	数值	参数	数值
$n$	0.3	${n_{{\text{vg}}}}$	1.698
$S_{{\text{res}}}^{\text{l}}$	0.15	v	1.4
$S_{{\text{sat}}}^{\text{l}}$	0.85	${f_{\tau {\text{v}}}}$	1.2
${m_{{\text{vg}}}}$	0.411	${P_{\text{s}}}$ /kPa	68

下载: 导出CSV

| 显示表格

图 1 体积含水率的对比验证

Figure 1. Comparative verification of volume water content

下载: 全尺寸图片幻灯片

图 2 含盐率的对比验证

Figure 2. Comparative verification of salt content

下载: 全尺寸图片幻灯片

图 3 温度的对比验证

Figure 3. Comparative verification of temperature

下载: 全尺寸图片幻灯片

3. 参数分析

表 2，3^[18]分别给出了粉质黏土与砂土的VG模型参数。非饱和黏土上覆非饱和砂土的高度为300 mm保持不变，对砂土不同厚度、渗透系数等条件下复合非饱和土中的水分、盐分随时间变化的分布情况以及变化规律进行分析。砂土的孔隙率为0.5，上覆粉质黏土的孔隙率为0.3、土体上下边界的温度分别为27℃，3℃，上下边界的孔隙气压力均为102 kPa。复合非饱和土初始时刻的温度为15℃、初始时刻的孔隙气压力为102 kPa，为了让迁移规律更加明显砂土与黏土中初始的孔隙水压力设为-800 kPa、含盐率为0.12%。

表 3 砂土的VG模型参数

Table 3. VG model parameters of sand soil

参数	数值	参数	数值
n	0.5	${n_{{\text{vg}}}}$	2.45
$S_{{\text{res}}}^{\text{l}}$	0.07	v	2
$S_{{\text{sat}}}^{\text{l}}$	0.93	${f_{\tau {\text{v}}}}$	1.4
${m_{{\text{vg}}}}$	0.59	${P_{\text{s}}}$ /kPa	32

下载: 导出CSV

| 显示表格

图 4，5分别给出了在砂土厚度分别为50，100，150 mm，且复合非饱和土厚度始终为300 mm条件下，复合土体中孔压和溶质的分布情况。从图 4中可以明显地看出随着砂土厚度的增加，复合土体中砂土中的盐分是降低的。图 5给出了距热端25，75，125，75，225，275 mm位置处复合非饱和土孔隙水压力随时间的变化情况。随着砂土厚度增大，相同位置处复合非饱和土热端的孔隙水压力越小，继而含水率越小。当孔隙率为0.5的砂土厚度越大，砂土对复合非饱和土中盐分迁移和水分迁移阻碍能力就越大。

图 4 不同砂土厚度下含盐率的变化

Figure 4. Variation of salt content under different sand thicknesses

下载: 全尺寸图片幻灯片

图 5 不同砂土厚度下不同位置处孔隙水压力随时间变化

Figure 5. Variation of pore water pressure with time at different locations under different sand thicknesses

下载: 全尺寸图片幻灯片

为了分析砂土渗透系数对复合非饱和土中水-热-盐耦合行为的影响，当砂土与粉质黏土的孔隙率以及复合非饱和土体的边界和初始条件保持不变。图 6绘出了砂土厚度为100 mm，复合非饱和土厚度为300 mm，砂土的渗透系数分别为1.2×10^-4，4.2×10^-4，8.2×10^-4 m/s条件下，复合非饱和土体中盐分的变化情况。由图中曲线可以看出：随着砂土渗透系数的增大，复合非饱和土热端处溶液中的盐分浓度是增加。图 7给出了距低温端25，75，125，175，225，275 mm位置处孔隙水压力随时间的变化情况，随着砂土渗透系数增大，相同位置处复合非饱和土热端处的孔隙水压力越小，继而含水率越小。

图 6 不同渗透系数下含盐率的变化

Figure 6. Change of salt content under different permeability coefficients

下载: 全尺寸图片幻灯片

图 7 不同渗透系数下不同位置处孔隙水压力随时间变化

Figure 7. Change of pore water pressure with time at different positions under different permeability coefficients

下载: 全尺寸图片幻灯片

4. 结论

文章建立了完善的双层非饱和多孔介质中水-热-盐耦合运移的数学模型，并进行相应的参数分析和试验研究。

（1）复合非饱和土中，盐分在温度梯度的作用下逐渐向热端聚集，水分在温度梯度的作用下逐渐向冷端聚集。

（2）随着砂土厚度的增大，复合非饱和土低温端聚集的水分是逐渐增加的，高温端的水分是逐渐减少的。高温端砂土中聚集的盐分会随着砂土厚度的增大而逐渐降低，低温端粉质黏土中的盐分随着砂土厚度的增大而增加。

（3）随着砂土渗透系数的增大，复合非饱和土低温端的水分降低。高温端砂土中聚集的盐分会随着砂土渗透系数的增大而逐渐增加，低温端粉质黏土中的盐分随着砂土渗透系数的增大而减小。

图 1 体积含水率的对比验证

Figure 1. Comparative verification of volume water content

下载: 全尺寸图片幻灯片

图 2 含盐率的对比验证

Figure 2. Comparative verification of salt content

下载: 全尺寸图片幻灯片

图 3 温度的对比验证

Figure 3. Comparative verification of temperature

下载: 全尺寸图片幻灯片

图 4 不同砂土厚度下含盐率的变化

Figure 4. Variation of salt content under different sand thicknesses

下载: 全尺寸图片幻灯片

图 5 不同砂土厚度下不同位置处孔隙水压力随时间变化

Figure 5. Variation of pore water pressure with time at different locations under different sand thicknesses

下载: 全尺寸图片幻灯片

图 6 不同渗透系数下含盐率的变化

Figure 6. Change of salt content under different permeability coefficients

下载: 全尺寸图片幻灯片

图 7 不同渗透系数下不同位置处孔隙水压力随时间变化

Figure 7. Change of pore water pressure with time at different positions under different permeability coefficients

下载: 全尺寸图片幻灯片

表 1 多孔介质中各物理成分的热力学参数

Table 1 Thermodynamic parameters of physical components in porous media

参数	数值	参数	数值
${c_{\text{a}}}$ /(J·kg^-1·K^-1)	1000	${L_{{\text{gw}}}}$ /(J·kg^-1)	2413000
$c_{\text{w}}^{\text{l}}$ /(J·kg^-1·K^-1)	4180	${k_0}$ /(m·s^-1)	2.89×10^-6
$c_{\text{w}}^{\text{g}}$ /(J·kg^-1·K^-1)	1900	g/(m·s^-2)	9.8
$c_{\text{p}}^{\text{l}}$ /(J·kg^-1·K^-1)	3750	$\rho _{\text{w}}^{\text{g}}$ /(kg^-1·m^-3)	1000
${\lambda _{\text{g}}}$ /(W·m^-1·K^-1)	0.024	${M_{\text{a}}}$ /(kg^-1·mol^-1）	0.0288
${\lambda _{\text{s}}}$ /(W·m^-1·K^-1)	2.93	${M_{\text{w}}}$ /(kg^-1·mol^-1)	0.018016
${\lambda _{\text{l}}}$ /(W·m^-1·K^-1)	0.56	$R$ /(J·mol^-1·K^-1)	8.3144

下载: 导出CSV

表 2 粉质黏土VG模型参数

Table 2 VG model parameters of silty clay

参数	数值	参数	数值
$n$	0.3	${n_{{\text{vg}}}}$	1.698
$S_{{\text{res}}}^{\text{l}}$	0.15	v	1.4
$S_{{\text{sat}}}^{\text{l}}$	0.85	${f_{\tau {\text{v}}}}$	1.2
${m_{{\text{vg}}}}$	0.411	${P_{\text{s}}}$ /kPa	68

下载: 导出CSV

表 3 砂土的VG模型参数

Table 3 VG model parameters of sand soil

参数	数值	参数	数值
n	0.5	${n_{{\text{vg}}}}$	2.45
$S_{{\text{res}}}^{\text{l}}$	0.07	v	2
$S_{{\text{sat}}}^{\text{l}}$	0.93	${f_{\tau {\text{v}}}}$	1.4
${m_{{\text{vg}}}}$	0.59	${P_{\text{s}}}$ /kPa	32

下载: 导出CSV

参考文献(18)

[1]	RESHETIN O L, ORLOV S Y. Theory of heat and moisture transfer in a capillary-porous body[J]. Technical Physics, 1998, 43(2): 263–264. doi: 10.1134/1.1258982
[2]	JIA Yong ying, WANG Zhi guo. The study of coupled transfer of heat, moisture and air in reservoir porous media[J]. Applied Mechanics and Materials, 2012, 271/272: 1195–1200. doi: 10.4028/www.scientific.net/AMM.271-272.1195
[3]	任荣. 非等温条件下土壤水热耦合迁移数值模拟研究[D]. 太原: 太原理工大学, 2018. REN Rong. Numerical Simulation of Coupled Soil Water and Heat Transfer under Non-Isothermal Conditions[D]. Taiyuan: Taiyuan University of Technology, 2018. (in Chinese)
[4]	DEB S K, SHUKLA M K, SHARMA P, et al. Coupled liquid water, water vapor, and heat transport simulations in an unsaturated zone of a sandy loam field[J]. Soil Science, 2011, 176(8): 387–398. doi: 10.1097/SS.0b013e318221f132
[5]	安然. 地下水流与气流、热流的典型耦合问题及解耦条件研究[D]. 北京: 中国地质大学(北京), 2015. AN Ran. Research on Typical Coupling Problems and the Decoupling Conditions for Groundwater Flow Linked with Air Flow and Thermal Flow[D]. Beijing: China University of Geosciences, 2015. (in Chinese)
[6]	姜建梅. 基于滨海平原区浅层地下水对土壤水汽热耦合运移规律的影响研究[D]. 天津: 天津大学, 2015. JIANG Jian-mei. The Effect of Shallow Groundwater on Coupled Soil Water Vapor and Heat Transport in Coastal Plain[D]. Tianjin: Tianjin University, 2015. (in Chinese)
[7]	JAHANGIR M H, SADRNEJAD S A. A new coupled heat, moisture and air transfer model in unsaturated soil[J]. Journal of Mechanical Science and Technology, 2012, 26(11): 3661–3672. doi: 10.1007/s12206-012-0839-z
[8]	HE Z Y, ZHANG S, TENG J D, et al. A coupled model for liquid water-vapor-heat migration in freezing soils[J]. Cold Regions Science and Technology, 2018, 148: 22–28. doi: 10.1016/j.coldregions.2018.01.003
[9]	CLEALL P J, SINGH R M, THOMAS H R. Vapour transfer in unsaturated compacted bentonite[J]. Géotechnique, 2013, 63(11): 957–964. doi: 10.1680/geot.12.P.147
[10]	HERNÁNDEZ-LÓPEZ M F, GIRONÁS J, BRAUD I, et al. Assessment of evaporation and water fluxes in a column of dry saline soil subject to different water table levels[J]. Hydrological Processes, 2014, 28(10): 3655–3669. doi: 10.1002/hyp.9912
[11]	CHEN Y P, SHI M H, LI X C. Experimental investigation on heat, moisture and salt transfer in soil[J]. International Communications in Heat and Mass Transfer, 2006, 33(9): 1122–1129. doi: 10.1016/j.icheatmasstransfer.2006.06.013
[12]	周凤玺, 高国耀. 非饱和土中热-湿-盐耦合作用的稳态分析[J]. 岩土力学, 2019, 40(6): 2050–2058. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201906004.htm ZHOU Feng-xi, GAO Guo-yao. Steady-state analysis of the heat-moisture-salt coupling for unsaturated soil[J]. Rock and Soil Mechanics, 2019, 40(6): 2050–2058. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201906004.htm
[13]	周凤玺, 高国耀. 非饱和土中热湿盐多场耦合过程分析[J]. 岩土工程学报, 2019, 41(5): 813–820. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201905004.htm ZHOU Feng-xi, GAO Guo-yao. Multi-field coupling process of heat- moisture-salt in unsaturated soil[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(5): 813–820. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201905004.htm
[14]	PHILIP J R, DE VRIES D A. Moisture movement in porous materials under temperature gradients[J]. Transactions, American Geophysical Union, 1957, 38(2): 222. doi: 10.1029/TR038i002p00222
[15]	WANG W Q, RUTQVIST J, GÖRKE U J, et al. Non-isothermal flow in low permeable porous media: a comparison of Richards' and two-phase flow approaches[J]. Environmental Earth Sciences, 2011, 62(6): 1197–1207. doi: 10.1007/s12665-010-0608-1
[16]	NASSAR I, HORTON R. Heat, water, and solution transfer in unsaturated porous media: I: theory development and transport coefficient evaluation[J]. Transport in Porous Media, 1997, 27: 17–38. doi: 10.1023/A:1006583918576
[17]	肖泽岸, 赖远明, 尤哲敏. 单向冻结过程中NaCl盐渍土水盐运移及变形机理研究[J]. 岩土工程学报, 2017, 39(11): 1992–2001. doi: 10.11779/CJGE201711006 XIAO Ze-an, LAI Yuan-ming, YOU Zhe-min. Water and salt migration and deformation mechanism of sodium chloridesoil during unidirectional freezing process[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(11): 1992–2001. (in Chinese) doi: 10.11779/CJGE201711006
[18]	陈世平, 李毅, 高金芳. 覆膜开孔入渗-蒸发条件下夹砂层土壤水、盐、热变化规律[J]. 中国农村水利水电, 2011(11): 47-51. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNSD201111014.htm CHEN Shi-ping, LI Yi, GAO Jin-fang. The movement of soil water, solute and heat for saline-alkali soil with sand layers under infiltration and plastic film hole evaporation[J]. China Rural Water and Hydropower, 2011(11): 47–51. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZNSD201111014.htm

施引文献

资源附件(0)

图(7) / 表(3)

计量

文章访问数: 211
HTML全文浏览量: 34
PDF下载量: 39
被引次数: 0

0. 引言
1. 数学模型
1.1 质量守恒方程
1.2 能量平衡方程
1.3 复合非饱和土非稳态问题求解
2. 试验验证
3. 参数分析
4. 结论

复合非饱和土中水-热-盐多场耦合分析 English Version

作者简介: 高国耀(1994—)，男，博士研究生，主要从事非饱和土、冻土的多场耦合研究。E-mail: Guoyao-Gao@163.com

通讯作者: 郭伟, E-mail: guow@tju.edu.cn

计量

出版历程