一种原位测量非饱和土的气体扩散系数和渗透系数的方法与装置

刘红位; 王孟奇; 詹良通; 冯嵩; 吴涛

doi:10.11779/CJGE20221228

一种原位测量非饱和土的气体扩散系数和渗透系数的方法与装置 English Version

刘红位^1,,
王孟奇¹,
詹良通^{2, 3},
冯嵩^{4, 3, ,},
吴涛^{5, 3}

1.
福州大学紫金地质与矿业学院, 福建福州 350108
2.
浙江大学岩土工程研究所, 浙江杭州 310058
3.
浙江大学软弱土与环境土工教育部重点实验室, 浙江杭州 310058
4.
福州大学土木工程学院, 福建福州 350108
5.
江苏科技大学土木工程与建筑学院, 江苏镇江 212000

基金项目:

国家自然科学基金项目 52178320

国家自然科学基金项目 42177120

国家重点研发计划 2019YFC1806003

详细信息

作者简介:
刘红位(1987— )，女，博士，副教授，主要从事环境岩土工程方面的教学和科研工作。E-mail: hliuan@connect.ust.hk

通讯作者:
冯嵩, E-mail: sfengaa@connect.ust.hk

中图分类号: TU43
计量
- 文章访问数: 253
- HTML全文浏览量: 30
- PDF下载量: 80
出版历程
- 收稿日期: 2022-10-06
- 网络出版日期: 2024-05-14
- 刊出日期: 2024-04-30

Method and apparatus for measuring in-situ gas diffusion coefficient and permeability coefficient of unsaturated soils

LIU Hongwei^1,,
WANG Mengqi¹,
ZHAN Liangtong^{2, 3},
FENG Song^{4, 3, ,},
WU Tao^{5, 3}

1.
Zijin School of Geology and Mining, Fuzhou University, Fuzhou 350108, China
2.
Institute of Geotechnical Engineering, Zhejiang University, Hangzhou 310058, China
3.
Key Laboratory of Soft Soils and Geoenvironmental Engineering of Ministry of Education, Zhejiang University, Hangzhou 310058, China
4.
College of Civil Engineering, Fuzhou University, Fuzhou 350108, China
5.
College of Civil Engineering and Architecture, Jiangsu University of Science and Technology, Zhenjiang 212000, China

摘要

摘要: 气体扩散系数D_g和渗透系数K_g是研究非饱和土中气体运移的两个重要参数。推导了稳态条件下，非饱和土中气体运移参数的计算理论。该理论基于示踪气体浓度C与气体流量q₁之间及压力梯度P与气体流量q₂之间的变化规律，分别确定非饱和土的D_g，K_g。依据推导的理论公式，设计并研制了一种原位测量非饱和土的D_g，K_g的装置。开展土柱试验，利用新装置测量非饱和土的D_g，K_g，并采用单元试验和数值模拟验证新方法和装置的有效性。研究表明新方法测量的非饱和土的D_g，K_g与单元试验及数值模拟结果吻合良好。基于数值模拟开展的敏感性分析结果揭示：D_g测量准确度随土体的质量含水率w减少而降低，而K_g的测量准确度受含水量的影响不显著；非饱和土的D_g，K_g的测量准确度均随着注气球半径增加而降低，但是注气球半径对K_g的测量准确度的影响更为显著，当注气球半径增加到0.05 m时，其准确度约为55%；土体各向异性降低测量准确度，当竖向与水平向的气体运移参数的比值D_gv/D_gh和K_gv/K_gh均为0.1时，D_g，K_g的测量准确度仍可达到40%以上。
- 非饱和土 /
- 气体扩散系数 /
- 气体渗透系数 /
- 原位测量装置 /
- 数值模拟
Abstract: The gas diffusion (D_g) and permeability (K_g) coefficients are the two important parameters for the study on gas transport in unsaturated soils. In the present study, the theory of calculating D_g and K_g for unsaturated soils under steady state is derived. D_g and K_g of unsaturated soils are determined based on the relationship between tracer gas concentration (C) and gas flow rate (q₁) and that between pressure gradient (p) and gas flow rate (q₂), respectively. A device for in-situ measurement of D_g and K_g in unsaturated soils is designed and developed. The soil column tests are conducted to measure D_g and K_g of unsaturated soils with different water contents. Thereafter, the element tests and numerical simulation analysis are used to verify the effectiveness of the new method and device. It is found that D_g and K_g obtained by the new method are in good agreement with the results of the element tests and numerical simulations. The results of the sensitivity analysis show that the measurement accuracy of D_g reduces as the soil water content decreases, while the measurement accuracy of K_g is independent of water content of soil. The measurement accuracies of both D_g and K_g decrease with the increase of the radius of the aeration bulb, but the measurement of K_g is more sensitive to the change of the radius. The measurement accuracy of K_g is about 55% of the true value, when the radius of the aeration bulb is 0.05 m. The soil anisotropy reduces the measurement accuracy; when the values of D_gv/D_gh and K_gv/K_gh is 0.1, the measurement accuracies of D_g and K_g can still reach more than 40%.
- unsaturated soil /
- gas diffusion coefficient /
- gas permeability coefficient /
- in-situ measurement device /
- numerical simulation

HTML全文

0. 引言

潜水面上升可能诱发边坡失稳，艾伯凡矿灾^[1]（Aberfan Disaster）、沃楚西特水坝^[2]（Wachusett Dam）等就属于这一破坏模式。Brand^[3]建议研究边坡的破坏机理时，应该考虑边坡内土体单元的应力应变响应。边坡内潜水面上升的应力路径可以简化为等剪应力排水路径^[3-6]。针对等剪应力排水路径下土体的力学响应，研究者们进行了一系列工作^{[2, 6-7]}。Chu等^[2]通过试验证明了等剪应力排水路径下土体的非稳定态与边坡失稳存在相关性。Lashkari等^[7]利用离散元数值模拟等剪应力排水路径，发现松砂和密砂的非稳定态存在区别。此外，研究者还发现各向异性对等剪应力排水路径的非稳定态存在显著影响。Fanni等^[8]通过试验发现大主应力方向从垂直于沉积面旋转到45°时，非稳定态的应力比不断减小，指出在评估边坡稳定时不能忽略各向异性的影响。

一些研究者也在利用本构模型模拟该应力路径。Lashkari等^[7]采用了各向异性临界状态塑性本构模型^[9]模拟了不同孔隙比和剪应力值下的等剪应力排水试验，表明了模型具有模拟该应力路径的能力。但本构模型能否模拟各向异性条件下的等剪应力排水试验，目前仍然未知。

各向异性对砂土的力学响应有显著影响^[10-16]。如何用本构模型描述饱和砂土的力学响应是目前的研究热点之一。为总结各向异性对于砂土力学响应影响的规律，研究者们进行了大量物理试验和离散元数值试验。物理试验方面如Yoshimine等^[17]和Nakata等^[18]研究了不同加载方向和中主应力系数对丰浦砂应力应变关系的影响，发现当大主应力平行于沉积面、中主应力系数为1时，土样表现出更强的剪缩特性，更易发生静态液化。离散元数值试验方面如Fu等^[19]再现了直剪试验，研究各向异性砂土的峰值强度。这些试验为本构模型的建立和验证提供了基础和依据，各类模型通过模拟试验证明了自身的有效性^[20-23]。如Wang等^[22]通过离散元数值试验总结各向异性对塑性模量和剪胀的影响，在CycLiq本构模型^[24]的基础上加入了组构各向异性，并通过离散元试验和物理试验对新模型进行了验证。

本文利用Wang等^[22]提出的各向异性CycLiq本构模型，对Fanni等^[8]进行的一系列定轴不排水试验与定轴等剪应力排水试验进行了模拟，验证了该模型对于各向异性砂土在不同主应力方向下、不同应力路径力学响应的模拟能力，并研究了初始组构强度对于模型模拟定轴等剪应力排水试验结果的影响。

1. 各向异性CycLiq模型

1.1 本构模型

基于排水与不排水循环扭剪试验现象研究，张建民等^[25-26]将砂土的体应变 ${\varepsilon _\text{v}}$ 分解为有效球应力变化引起的弹性体应变 ${\varepsilon _\text{vc}}$ 和剪切引起的塑性体应变 ${\varepsilon _\text{vd}}$ ，又将剪切引起的塑性体应变 ${\varepsilon _\text{vd}}$ 进一步划分为可逆性剪切体应变 ${\varepsilon _\text{vd, re}}$ 与不可逆性剪切体应变 ${\varepsilon _\text{v, ir}}$ 。在此基础上，王刚等^[27]和Wang等^[24]建立了一个能够模拟饱和砂土液化前后应力应变响应的弹塑性循环本构模型，并在模型中引入了状态变量，使得模型能够统一描述砂土在不同状态下的静动力响应。为考虑各向异性，Wang等^[22]基于试验和离散元数值模拟总结的规律，对模型进行了改进，引入了组构张量，在塑性模量与剪胀率的计算中引入了各向异性的影响，使得模型能够描述不同主应力方向下的响应。

CycLiq模型的应力应变关系可以用下式来描述：

$\boldsymbol{\dot \varepsilon} = \frac{{p \boldsymbol{\dot r}}}{{2G}} + \left( {\frac{\boldsymbol{r}}{{2G}} + \frac{\boldsymbol{I}}{{3K}}} \right)\dot p + \left( {\boldsymbol{n} + \frac{D}{3}\boldsymbol{I}} \right)\left\langle L \right\rangle 。$

(1)

式中：p为有效球应力；r为偏应力张量s与球应力p之比；G和K为弹性剪切模量和体积模量；I为单位张量；n为塑性偏应变流动方向与加载方向；D为剪胀率；L为塑性加载指数。参数上加(·)表示增量；〈〉为Macaulay计算符，当x > 0时，〈x〉= x，当x ≤ 0时，〈x〉= 0。

弹性剪切模量和体积模量采用Richart等^[28]提出的线性弹性模量公式，在此不再赘述。塑性加载指数L通过下式计算：

$L = \frac{{p \boldsymbol{\dot r}:\boldsymbol{n}}}{H} 。$

(2)

式中：H为塑性模量。

砂土临界状态采用Been等^[29]提出的状态变量Ψ = $e - {e_\text{c}}$ ，但 ${e_\text{c}}$ 与p的关系采用了与Wang等^[22]所使用的指数关系^[30]不同的对数关系：

${e_\text{c}} = \varGamma - {\lambda _\text{e}}\ln (p) 。$

(3)

这一改动是因为Fanni等^[8]给出了本文模拟试样临界状态下试样的e-p曲线，并采用式（3）的对数关系进行了拟合，且e-p曲线的对数关系被广泛采用，许多模型采用了这种关系^[31-32]，因此对模型中的公式进行了替换。不过需要指出的是，在液化状态时，对数关系不再适用，因此砂土的动力本构模型大都采用指数关系^[33-34]。本文中不涉及液化内容，因此可以选取对数关系进行模拟。

为了考虑各向异性的影响，Wang等^[22]在模型中引入了组构张量 ${\boldsymbol{F}_\text{n}}$ ，并考虑了组构张量对于塑性模量与剪胀率的影响。组构张量 ${\boldsymbol{F}_\text{n}}$ 的演化采用了如下公式：

$\text{d}{\boldsymbol{F}_\text{n}} = \left\langle L \right\rangle c\left[ {\boldsymbol{n} - (1 + D){\boldsymbol{F}_\text{n}}} \right] 。$

(4)

式中：c为模量参数，代表组构张量演化的速率。通过总结离散元数值试验的规律，Wang等^[22]发现组构张量的模|| ${\boldsymbol{F}_\text{n}}$ ||与组构张量和加载方向的联合不变量A_n = ${\boldsymbol{F}_\text{n}}$ : n之差|| ${\boldsymbol{F}_\text{n}}$ ||-A_n可以作为一个良好的指标，以反映塑性模量和剪胀率随加载角度的变化。以水平沉积的试样在三轴压缩中的加载角度为0°，随着加载角度的增加，指标|| ${\boldsymbol{F}_\text{n}}$ ||-A_n不断增大，同时塑性模量不断减小，剪胀率从不断增加。因此，将|| ${\boldsymbol{F}_\text{n}}$ ||-A_n引入到了塑性模量与剪胀率的计算公式中：

$C = \exp \left[ { - {\Delta _1}{{(\left\| {{\boldsymbol{F}_\text{n}}} \right\| - {A_\text{n}})}^{{1 \mathord{\left/ {\vphantom {1 2}} \right. } 2}}}} \right] \text{，}$

(5)

$\begin{align}H = &\frac{2}{3}hg\left( \theta \right)G\exp ( - {n^p}\Psi ) \cdot \\& \left( {\frac{{M\exp ( - {n^{\rm{p}}}\Psi )}}{{{M_\text{m}}}}\frac{{\bar \rho }}{\rho } - \frac{1}{C}} \right) \text{，} \end{align}$

(6)

${D_\text{ir}} = \frac{{{d_\text{ir}} + {\Delta _2}(\left\| {{\boldsymbol{F}_\text{n}}} \right\| - {A_\text{n}})}}{{{d_\text{ir}}}}{D_\text{ir, iso}} 。$

(7)

式中：h，Δ₁，Δ₂，M，d_ir， ${n^{\rm{p}}}$ 为模型参数；g(θ)为形状函数，这里采用的是Zhang等^[35]建议的函数； ${M_\text{m}}$ 为历史最大应力比；D_{ir, iso}为各向同性条件下的不可逆剪胀率，详细公式可见Wang等^[24]。

1.2 参数标定

Fanni等^[8]试验中采用的砂为一种细硅砂，颗粒相对质量密度为2.64，平均粒径D₅₀为0.21 mm，均匀系数C_u为2.3，颗粒呈圆形至亚圆形。因缺乏相关试验，大部分参数参照了Wang等^[22]标定的丰浦砂参数，仅对部分可确定的临界状态参数（n^p，Γ， ${\lambda _\text{e}}$ ）进行了重新标定。标定结果如表 1所示。

其中，Γ和 ${\lambda _\text{e}}$ 由Fanni等^[8]根据试验拟合得到，在更换了Wang等^[22]模型中e_c与p的关系式后可直接使用。展示了 ${n^{\rm{p}}}$ 的标定方法，以Fanni等^[8]进行的一组定轴等剪应力排水试验（中主应力系数b = 0.2，加载角度α = 30°）结果参照进行标定。图中点代表试验所得试样到达非稳定态时的有效球应力与孔隙比，虚线和实线分别为n^p = 6和12时的模拟结果。可以看到当其他参数不变时， ${n^{\rm{p}}}$ 值越大，模拟得到的非稳定态时有效球应力随孔隙比变化越快，曲线更加陡峭。这是由于等剪应力排水路径中的非稳定态代表试样处于峰值应力比状态，而峰值应力比可由Mexp(- ${n^{\rm{p}}}$ Ψ)决定。

图 1 模型参数n^p的标定方法演示

Figure 1. Calibration method for model parameter n^p

下载: 全尺寸图片幻灯片

除了表 1中列出的模型参数外，初始组构强度，即初始组构张量的模||F_nin||，也是一个重要的参数。它代表了试样在试验开始时的各向异性程度，||F_nin||越大，试样的各向异性越强。不同的制样方式制得的试样往往有着不同的各向异性程度^[36]，因此需要取不同的初始组构强度。本次模拟中涉及两种不同的制样方式，湿装夯实法（moist tamped, MT）和干砂雨法（dry pluviated, DP），其中干砂雨法得到的试样孔隙比要小于湿装夯实法，结合Yang等^[37]、叶斌等^[38]通过物理试验与离散元模拟均发现干砂雨法的各向异性强于湿装夯实法，可以推测干砂雨法在制样中，由于受到重力影响，土颗粒的排列具有定向性，而湿装夯实法在制样中，土颗粒更倾向于随机排列，导致其各向异性相对较弱。最终，湿装夯实法制得的试样初始组构强度||F_nin||取为0.5，干砂雨法制得的试样初始组构强度||F_nin||取为0.7。

2. 不同主应力方向下模型的响应

2.1 定轴不排水剪切试验

Fanni等^[8]进行了一系列空心圆柱定轴不排水剪切试验。试验中，初始有效球应力均为171 kPa，大小主应力之差均为64 kPa，中主应力系数b分别为0，0.2，0.5，大主应力与沉积方向的夹角α分别为0°，22.5°，30°和45°，两种制样方法，即湿装夯实法和干砂雨法，各5组试验，如图 2，3所示。其中，湿装夯实法所得试样孔隙比在0.761~0.785，干砂雨法所得试样孔隙比在0.674~0.697。采用表 1中的模型参数对试验进行模拟，模拟结果如图 2，3所示。图 2（a），（b）中横纵坐标为Fanni等^[8]给出的八面体剪应变和八面体偏应力，图 3同理。

图 2 不同中主应力系数与主应力方向下湿装夯实法试样定轴不排水试验结果与模拟

Figure 2. Undrained torsional shear hollow cylinder tests and simulated results at different α and b of moist tamped loose sand

下载: 全尺寸图片幻灯片

图 3 不同中主应力系数与主应力方向下干砂雨法试样定轴不排水试验结果与模拟

Figure 3. Undrained torsional shear hollow cylinder tests and simulated results at different α and b of dry pluviated loose sand

下载: 全尺寸图片幻灯片

表 1 各向异性CycLiq模型参数

Table 1. Model parameters of anisotropic CycLiq model

G_o	κ	h	d_{re, 1}	d_{re, 2}	d_ir	α	γ_{d, r}	n^p	n^d	M	λ_e	Γ	Δ₁	Δ₂	c
100	0.008	4.0	0.7	30	0.3	30	0.05	12	2	1.2	0.008	0.761	1.35	0.8	5

下载: 导出CSV

| 显示表格

图 2展示了湿装夯实法试样定轴不排水试验结果与模拟。湿装夯实法得到的试样均为松砂，各试验和模拟均出现静态液化。总体上来看，随着b和α的增加，试验与模拟中试样的峰值偏应力均不断减小。从图 2（a），（c）中可以看出，试验中b = 0.2的3组曲线极为接近，甚至α = 30°的偏应力要高于α = 22.5°，这主要是试验中不同试样的孔隙比相差较大造成的，在模拟中该现象并不明显。此外，b = 0.5的模拟结果在加载初期剪应力出现了明显的减小，这主要是由于本构模型中组构对塑性模量的影响描述不够准确，在大孔隙比下放大了组构对塑性模量的影响。

图 3展示了干砂雨法试样定轴不排水试验结果与模拟。相比于湿装夯实法，干砂雨法得到的试样孔隙比更小，可认为密砂。各组模拟的曲线与试验结果均有较好的吻合。对于b = 0.2的3组试验和模拟，α = 22.5°时均表现出明显的剪胀，而在α = 45°时则表现出强烈的剪缩。对于b = 0，α = 0°和b = 0.2，α = 22.5°两组试验，试验和模拟的曲线在图 3（a），（b）均出现了交叉，这是由于两者的孔隙比相差较大，在孔隙比、中主应力系数、主应力方向三者的共同作用下出现了图中的效果。

2.2 定轴等剪应力排水试验

在定轴不排水剪切试验的基础上，Fanni等^[8]进行了一系列定轴等剪应力排水试验，以探究各向异性对等剪应力排水应力路径的影响。试验中采用了湿装夯实法（MT）和干砂雨法（DP）两种制样方法，初始有效球应力均为171 kPa，大小主应力之差均为64 kPa，中主应力系数b分别为0，0.2，主应力与沉积方向的夹角α分别为0°，22.5°和45°。每个夹角均用不同的孔隙比进行试验，对于湿装夯实法，b = 0，α = 0°组采用了6个不同的孔隙比，b = 0.2，α = 22.5°组采用了4个，b = 0.2，α = 45°组采用了10个；对于干砂雨法， b = 0.2，α = 22.5°组1个，b = 0.2，α = 45°组2个。采用表 1中的模型参数进行模拟，模拟中保持偏应力不变，当有效球应力无法再减小时，即试样达到峰值应力比时，认为此时试样处于非稳定态，模拟结束。

图 4给出了湿装夯实法试样在不同b和α下的定轴等剪应力排水试验结果与模拟。图中圆代表初始状态，三角形代表试样达到非稳定态，即处于峰值应力比状态。当相同b和α时，非稳定态的有效球应力p随孔隙比e的增加而增加，3组中模拟与试验均吻合较好。当孔隙比e不变、α增加时，模拟与试验得到的有效球应力均增加。

图 4 不同中主应力系数与主应力方向下湿装夯实法试样定轴等剪应力排水试验结果与模拟

Figure 4. Constant shear drained tests and simulated results at different α and b of moist tamped loose sand

下载: 全尺寸图片幻灯片

图 5给出了干砂雨法试样在不同b和α下的定轴等剪应力排水试验结果与模拟。b = 0.2，α = 45°组试验与模拟结果相近，规律也与湿装夯实法中一致；模拟结果在定量上与试验结果存在差距，推测这是由于模型参数或者初始组构强度不准确导致的。当α增加时，模拟与试验得到的有效球应力均增加，这与湿装夯实法中的规律是一致的。

图 5 不同主应力方向下干砂雨法试样定轴等剪应力排水试验结果与模拟

Figure 5. Constant shear drained tests and simulated results at different α and b of dry pluviated loose sand

下载: 全尺寸图片幻灯片

图 6给出了试验和模拟中非稳定态的结果与趋势线。可以看到试验得到的非稳定态在e-p空间内具有良好的线性关系，对模拟得到的非稳定态进行线性回归，得到R²均在0.99以上，表明模型对等剪应力比试验具有良好的模拟能力。对比湿装夯实法的3组试验，可以发现模拟得到的3条趋势线之间的间距要略小于试验所得，这主要受初始组构强度的影响，具体情况在下文中阐释。

图 6 不同中主应力系数、主应力方向下定轴等剪应力排水试验中非稳定态试验与模拟结果及各自的趋势线

Figure 6. Instability states of constant shear drained tests and simulations at different α and b

下载: 全尺寸图片幻灯片

为了说明初始组构强度对于定轴等剪应力排水条件下土体响应的影响，在数值模拟中对不同初始组构强度条件进行了分析。中给出了不同初始组构强度对非稳定态模拟结果的影响。从图中可以发现初始组构越强，不同主应力方向下的非稳定态差异越大，即峰值应力比差异越大。式（5），（6）中可以解释这一现象。试样达到峰值应力比时，有塑性模量H = 0。当初始组构强度增加，即|| ${\boldsymbol{F}_\text{n}}$ ||增加时，指标|| ${\boldsymbol{F}_\text{n}}$ ||-A_n = || ${\boldsymbol{F}_\text{n}}$ ||(1 – n_F: n)（其中n_F指组构张量的方向）的值域[0, 2|| ${\boldsymbol{F}_\text{n}}$ ||]随之扩大，式（5）所得C的值域也相应扩大。对于定轴试验，初始组构强度增加，导致了不同的b和α计算得到C值之间的差异增加，因H = 0，进而使得式（6）中历史最大峰值应力比M_m（即非稳定态时的应力比）的差异增加，在剪应力q不变的情况下，非稳定态时的球应力p之间的差异也就增加了。这些反映在图 7中即为趋势线之间的间距增加。

图 7 初始组构强度对定轴等剪应力排水试验非稳定态e-p关系的影响

Figure 7. Effects of initial fabric norm on instability state of constant shear drained loading at different α and b

下载: 全尺寸图片幻灯片

3. 结论

本文利用各向异性CycLiq本构模型，模拟了Fanni等^[8]进行的定轴不排水试验与定轴等剪应力排水试验，对比试验结果与模拟结果，可以得到以下3点结论。

（1）对于定轴不排水剪切试验，该模型与试验均表现为随着中主应力系数和加载角度的增加，土体剪缩增强，更易发生静态液化。在定量上，模型模拟的应力路径和应力应变曲线与试验也有较好的吻合。

（2）对于定轴等剪应力排水试验，模型与试验均表现为非稳定态时的平均有效应力随孔隙比和加载角度的增加而增加，且相同加载角度的非稳定态在e-p空间内呈线性分布。总体来说，CycLiq模型对各向异性下砂土等剪应力排水响应的模拟效果较好，具有运用到实际工程中的潜力。

（3）不同初始组构强度对于等剪应力排水试验的结果有显著的影响。初始组构强度越高，不同主应力方向作用下的非稳定态时平均有效应力相差越大。

图 1 非饱和土中球形气流示意图

Figure 1. Schematic diagram of spherical gas flow in unsaturated soil

下载: 全尺寸图片幻灯片

图 2 非饱和土的气体扩散系数和气体渗透系数原位测试装置示意图

Figure 2. Schematic diagram of in-situ instruments of gas diffusion coefficient and gas coefficient permeability of unsaturated soils

下载: 全尺寸图片幻灯片

图 3 数值模型尺寸和网格剖分

Figure 3. Sizes and meshing of numerical model

下载: 全尺寸图片幻灯片

图 4 土柱试验的原始测量记录

Figure 4. Raw measurement records of soil column tests

下载: 全尺寸图片幻灯片

图 5 土柱试验和数值模拟值

Figure 5. Soil column tests and numerical simulations: (a) relationship between CO₂ concentration and time, (b) relationship between gas pressure and influx of gas

下载: 全尺寸图片幻灯片

图 6 单元试验，土柱试验和数值模拟值的结果对比

Figure 6. Comparison of measured results among element tests, soil column tests and numerical simulations

下载: 全尺寸图片幻灯片

图 7 注气球半径对测量准确度的影响

Figure 7. Effects of aeration bulb radius (r₀) on measurement accuracy

下载: 全尺寸图片幻灯片

图 8 不同深度对测量准确度的影响

Figure 8. Effects of different depths on measurement accuracy

下载: 全尺寸图片幻灯片

图 9 土体各向异性对测量准确度的影响

Figure 9. Effects of soil anisotropy on measurement accuracy

下载: 全尺寸图片幻灯片

表 1 土柱试验用土的基本性质

Table 1 Index properties of soils used in column tests

基本性质	土柱试验	参考文献
土的分类	SM	ASTM D 2487
相对质量密度	2.66	ASTM D 854
颗粒组成		ASTM D 422
粗砂（2~4.75 mm）/%	0
中砂（0.425~2 mm）/%	53.60
细砂（0.075~0.425 mm）/%	26.13
粉粒（0.002~0.075 mm）/%	15.87
黏粒（＜0.002 mm）/%	4.40
液塑限		ASTM D 4318
液限/%	37.8
塑限/%	28.2
塑性指数	9.6
最优含水率/%	17.12	ASTM D 698
最大干密度/(g·cm^-3)	1.84	ASTM D 698

下载: 导出CSV

表 2 敏感性分析模拟工况

Table 2 Summary of sensitivity analysis

系列	r₀/ m	注气球埋置深度/m	K_gv/ (10^-7 m·s)	D_gv/ (10^-7 m·s)	土体各向异性（D_gv/D_gh，K_gv/K_gh）^①
1	0.01~0.05	0.35	6.43, 4^② 4.54, 1.43	31.0, 7.46, 2.71	1
2	0.02	0.05~0.45			1
3	0.01	0.35			0.1~1
注：①D_gv，D_gh分别为竖向与水平向的土体的气体扩散系数，K_gv，K_gh分别为竖向与水平向的土体的气体渗透系数。②w=12.5%时，K_gv=6.43×10^-7 m/s，D_gv=31.0×10^-7 m²/s；w=18.0%时，K_gv=4.54×10^-7 m/s，D_gv=7.46×10^-7 m²/s；w=21.0%时，K_gv=1.43×10^-7 m/s，D_gv=2.71×10^-7 m²/s。

下载: 导出CSV

参考文献(43)

[1]	ALLAIRE S E, LAFOND J A, CABRAL A R, et al. Measurement of gas diffusion through soils: comparison of laboratory methods[J]. Journal of Environmental Monitoring, 2008, 10(11): 1326-1336. doi: 10.1039/b809461f
[2]	LIPIEC J, HORN R, PIETRUSIEWICZ J, et al. Effects of soil compaction on root elongation and anatomy of different cereal plant species[J]. Soil and Tillage Research, 2012, 121: 74-81. doi: 10.1016/j.still.2012.01.013
[3]	CHEN F Q, ZHAO N K, FENG S, et al. Effects of biochar content on gas diffusion coefficient of soil with different compactness and air contents[J]. Environmental Science and Pollution Research, 2020, 27(17): 21497-21505. doi: 10.1007/s11356-020-08594-7
[4]	FREDLUND D G, RAHARDJO H. Soil Mechanics for Unsaturated Soils[M]. New York: Wiley, 1993.
[5]	KUNCORO P H, KOGA K, SATTA N, et al. A study on the effect of compaction on transport properties of soil gas and water: I Relative gas diffusivity, air permeability, and saturated hydraulic conductivity[J]. Soil and Tillage Research, 2014, 143: 172-179. doi: 10.1016/j.still.2014.02.006
[6]	MOLDRUP P, OLESEN T, KOMATSU T, et al. Tortuosity, diffusivity, and permeability in the soil liquid and gaseous phases[J]. Soil Science Society of America Journal, 2001, 65(3): 613-623. doi: 10.2136/sssaj2001.653613x
[7]	DESUTTER T M, SAUER T J, PARKIN T B, et al. A subsurface, closed-loop system for soil carbon dioxide and its application to the gradient efflux approach[J]. Soil Science Society of America Journal, 2008, 72(1): 126-134. doi: 10.2136/sssaj2006.0101
[8]	PINGINTHA N, LECLERC M Y, BEASLEY J P, et al. Assessment of the soil CO₂ gradient method for soil CO₂ efflux measurements: comparison of six models in the calculation of the relative gas diffusion coefficient[J]. Tellus B: Chemical and Physical Meteorology, 2010, 62(1): 47. doi: 10.1111/j.1600-0889.2009.00445.x
[9]	ZHAN T L T, YANG Y B, CHEN R, et al. Influence of clod size and water content on gas permeability of a compacted loess[J]. Canadian Geotechnical Journal, 2014, 51(12): 1468-1474. doi: 10.1139/cgj-2014-0126
[10]	姚志华, 陈正汉, 黄雪峰, 等. 非饱和Q₃黄土渗气特性试验研究[J]. 岩石力学与工程学报, 2012, 31(6): 1264-1273. doi: 10.3969/j.issn.1000-6915.2012.06.023 YAO Zhihua, CHEN Zhenghan, HUANG Xuefeng, et al. Experimental research on gas permeability of unsaturated Q₃ loess[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(6): 1264-1273. (in Chinese) doi: 10.3969/j.issn.1000-6915.2012.06.023
[11]	陈云敏. 环境土工基本理论及工程应用[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
[12]	MCINTYRE D S, PHILIP J R. A field method for measurement of gas diffusion into soils[J]. Soil Research, 1964, 2(2): 133. doi: 10.1071/SR9640133
[13]	LAI S H, TIEDJE J M, ERICKSON A E. In situ measurement of gas diffusion coefficient in soils[J]. Soil Science Society of America Journal, 1976, 40(1): 3-6. doi: 10.2136/sssaj1976.03615995004000010006x
[14]	NICOT J P, BENNETT P C. Shallow subsurface characterization of gas transport in a playa wetland[J]. Journal of Environmental Engineering, 2015, 124(11): 1038-1046.
[15]	WEEKS E P, EARP D E, THOMPSON G M. Use of atmospheric fluorocarbons F-11 and F-12 to determine the diffusion parameters of the unsaturated zone in the Southern High Plains of Texas[J]. Water Resources Research, 1982, 18(5): 1365-1378. doi: 10.1029/WR018i005p01365
[16]	秦冰, 陆飏, 张发忠, 等. 考虑Klinkenberg效应的压实膨润土渗气特性研究[J]. 岩土工程学报, 2016, 38(12): 2194-2202. doi: 10.11779/CJGE201612007 QIN Bing, LU Yang, ZHANG Fazhong, et al. Gas permeability of compacted bentonite considering Klinkenberg effect[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(12): 2194-2202. (in Chinese) doi: 10.11779/CJGE201612007
[17]	陈存礼, 张登飞, 张洁, 等. 等向应力条件下原状Q₃黄土的渗气特性研究[J]. 岩土工程学报, 2017, 39(2): 287-294. doi: 10.11779/CJGE201702012 CHEN Cunli, ZHANG Dengfei, ZHANG Jie, et al. Gas permeability of intact Q₃ loess under isotropic stresses[J]. Chinese Journal of Geotechnical Engineering, 2017, 39(2): 287-294. (in Chinese) doi: 10.11779/CJGE201702012
[18]	苗强强, 陈正汉, 张磊, 等. 非饱和黏土质砂的渗气规律试验研究[J]. 岩土力学, 2010, 31(12): 3746-3750, 3757. doi: 10.3969/j.issn.1000-7598.2010.12.009 MIAO Qiangqiang, CHEN Zhenghan, ZHANG Lei, et al. Experimental study of gas permeability of unsaturated clayey sand[J]. Rock and Soil Mechanics, 2010, 31(12): 3746-3750, 3757. (in Chinese) doi: 10.3969/j.issn.1000-7598.2010.12.009
[19]	陈正汉. 非饱和土与特殊土力学[M]. 北京: 中国建筑工业出版社, 2022. CHEN Zhenghan. Mechanics for Unsaturated and Special Soils[M]. Beijing: China Architecture & Building Press, 2022. (in Chinese)
[20]	CHEN R, HUANG J W, ZHOU C, et al. A new simple and low-cost air permeameter for unsaturated soils[J]. Soil and Tillage Research, 2021, 213(12): 105083.
[21]	MOHAMMADI M H, VANCLOOSTER M. A simple device for field and laboratory measurements of soil air permeability[J]. Soil Science Society of America Journal, 2019, 83(1): 58-63. doi: 10.2136/sssaj2018.03.0114
[22]	JALBERT M, DANE J H. A handheld device for intrusive and nonintrusive field measurements of air permeability[J]. Vadose Zone Journal, 2003, 2(4): 611-617. doi: 10.2136/vzj2003.6110
[23]	IVERSEN B V, SCHJøNNING P, POULSEN T G, et al. In situ, on-site and laboratory measurements of soil air permeability: boundary conditions and measurement scale [J]. Soil Science, 2001, 166(2): 97-106. doi: 10.1097/00010694-200102000-00003
[24]	SHAN C, FALTA R W, JAVANDEL I. Analytical solutions for steady state gas flow to a soil vapor extraction well[J]. Water Resources Research, 1992, 28(4): 1105-1120. doi: 10.1029/91WR02986
[25]	OLSON M S, TILLMAN F D, CHOI J W, et al. Comparison of three techniques to measure unsaturated-zone air permeability at Picatinny Arsenal, NJ[J]. Journal of Contaminant Hydrology, 2001, 53(1/2): 1-19.
[26]	ZHAN L T, QIU Q W, XU W J, et al. Field measurement of gas permeability of compacted loess used as an earthen final cover for a municipal solid waste landfill[J]. Journal of Zhejiang University-Science A (Applied Physics & Engineering), 2016, 17(7): 541-552.
[27]	陈正汉, 谢定义, 王永胜. 非饱和土的水气运动规律及其工程性质研究[J]. 岩土工程学报, 1993, 15(3): 9-20. doi: 10.3321/j.issn:1000-4548.1993.03.002 CHEN Zhenghan, XIE Dingyi, WANG Yongsheng. Experimental studies of laws of fluid motion, suction and pore pressures in unsaturated soil[J]. Chinese Journal of Geotechnical Engineering, 1993, 15(3): 9-20. (in Chinese) doi: 10.3321/j.issn:1000-4548.1993.03.002
[28]	ROUF M A, SINGH R M, BOUAZZA A, et al. Gas permeability of partially hydrated geosynthetic clay liner under two stress conditions[J]. Environmental Geotechnics, 2016, 3(5): 325-333. doi: 10.1680/envgeo.14.00009
[29]	QIU Q W, ZHAN L T, LEUNG A K, et al. A new method and apparatus for measuring in situ air permeability of unsaturated soil[J]. Canadian Geotechnical Journal, 2021, 58(4): 514-530. doi: 10.1139/cgj-2019-0733
[30]	SCHJØNNING P, EDEN M, MOLDRUP P, et al. Two-chamber, two-gas and one-chamber, one-gas methods for measuring the soil-gas diffusion coefficient: validation and inter-calibration[J]. Soil Science Society of America Journal, 2013, 77(3): 729-740. doi: 10.2136/sssaj2012.0379
[31]	赵能凯. 根系特征对非饱和土体气体渗透与扩散的影响规律[D]. 福州: 福州大学, 2021. ZHAO Nengkai. Influence of Root Characteristics on Gas Permeation and Diffusion in Unsaturated Soils[D]. Fuzhou: Fuzhou University, 2021. (in Chinese)
[32]	MOLINS S, MAYER K U. Coupling between geochemical reactions and multicomponent gas and solute transport in unsaturated media: a reactive transport modeling study[J]. Water Resources Research, 2007, 43(5): 687-696.
[33]	NG C W W, FENG S, LIU H W. A fully coupled model for water-gas-heat reactive transport with methane oxidation in landfill covers[J]. Science of the Total Environment, 2015, 508: 307-319. doi: 10.1016/j.scitotenv.2014.11.037
[34]	REID R C, PRAUSNITZ J M, POLING B E. The Properties of Gases and Liquids[M]. 4th ed. New York: McGraw-Hill, 1987.
[35]	ALBERTY R A, SILBEY R J. Physical Chemistry[M]. New York: Wiley, 1992.
[36]	PARKER J C. Multiphase flow and transport in porous media[J]. Reviews of Geophysics, 1989, 27(3): 311-328. doi: 10.1029/RG027i003p00311
[37]	LAEMMEL T, MAIER M, SCHACK‐KIRCHNER H, et al. An in situ method for real‐time measurement of gas transport in soil[J]. European Journal of Soil Science, 2017, 68(2): 156-166. doi: 10.1111/ejss.12412
[38]	BARRAL C, OXARANGO L, PIERSON P. Characterizing the gas permeability of natural and synthetic materials[J]. Transport in Porous Media, 2010, 81(2): 277-293. doi: 10.1007/s11242-009-9398-x
[39]	AUBERTIN M, AACHIB M, AUTHIER K. Evaluation of diffusive gas flux through covers with a GCL[J]. Geotextiles and Geomembranes, 2000, 18(2/3/4): 215-233.
[40]	MOLDRUP P, OLESEN T, SCHJØNNING P, et al. Predicting the gas diffusion coefficient in undisturbed soil from soil water characteristics[J]. Soil Science Society of America Journal, 2000, 64(1): 94-100. doi: 10.2136/sssaj2000.64194x
[41]	HAMAMOTO S, MOLDRUP P, KAWAMOTO K, et al. Effect of particle size and soil compaction on gas transport parameters in variably saturated, sandy soils[J]. Vadose Zone Journal, 2009, 8(4): 986-995. doi: 10.2136/vzj2008.0157
[42]	施建勇, 赵义. 气体压力和孔隙对垃圾土体气体渗透系数影响的研究[J]. 岩土工程学报, 2015, 37(4): 586-593. doi: 10.11779/CJGE201504002 SHI Jianyong, ZHAO Yi. Influence of air pressure and void on permeability coefficient of air in municipal solid waste(MSW)[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(4): 586-593. (in Chinese) doi: 10.11779/CJGE201504002
[43]	FENG S, SUN J X, ZHAN L T, et al. A new method and instrument for measuring in situ gas diffusion coefficient and gas coefficient of permeability of unsaturated soil[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2023, 149(7): 04023041. doi: 10.1061/JGGEFK.GTENG-10871