原状黄土的结构性与渗气特性的关联性初探

张登飞; 陈存礼; 舒迎涛; 庞腾腾

doi:10.11779/CJGE202107021

原状黄土的结构性与渗气特性的关联性初探 English Version

1.
西北大学地质学系大陆动力学国家重点实验室,陕西　西安　710069
2.
西安理工大学岩土工程研究所,陕西　西安　710048
3.
哈尔滨电气国际工程有限责任公司,黑龙江　哈尔滨　150028
4.
西安市建筑设计研究院有限公司,陕西　西安　710054

基金项目:

国家自然科学基金项目 41907233

支撑服务川藏铁路规划建设地质调查 20190505

川西—藏东地区交通廊道活动构造与地质调查 DD20201123

详细信息

作者简介:
张登飞(1987— ),男,博士,主要从事黄土灾害与黄土力学方面的研究工作。E-mail：dfzhang87@nwu.edu.cn

中图分类号: TU411
计量
- 文章访问数: 230
- HTML全文浏览量: 30
- PDF下载量: 115
出版历程
- 收稿日期: 2020-10-25
- 网络出版日期: 2022-12-02
- 刊出日期: 2021-06-30

Correlation between structure of intact loess and gas permeability

1.
State Key Laboratory of Continental Dynamics, Department of Geology, Northwest University, Xi'an 710069, China
2.
Institute of Geotechnical Engineering, Xi'an University of Technology, Shaanxi Province, Xi'an 710048, China
3.
Harbin Electric International Company Limited, Heilongjiang Province, Harbin 150028, China
4.
Xi'an Architectural Design and Research Institute Co., Ltd., Xi'an 710054, China

摘要

摘要: 对3个地区的不同含水率原状黄土进行了等向压缩与渗气性同步测定试验,分析了原状黄土等向压缩过程中变形特性与渗气特性,探讨了压缩特性与渗气性之间的关联性通过引入用渗气性定义的结构性参数,提出了描述湿载耦合作用时原状黄土压缩变形的表达式。研究结果表明：原状黄土的渗气性变化与压缩特性的变化过程呈良好的关联性,渗气性可反映原状黄土组构与湿度不同而引起结构势差异及其应力作用引起黄土结构势的损伤过程；把结构屈服时某一含水率与参考含水率土的渗气系数之比定义为结构比参数,通过将不同含水率下结构比参数引入到等向应力变量中而得到了具有归一性的压缩曲线；据此,提出了以参考含水率的原状黄土的压缩曲线为基础,统一、连续描述非饱和原状黄土屈服前后压缩变形模式,其计算与试验结果皆吻合较好。
- 土的结构性 /
- 等向压缩特性 /
- 渗气特性 /
- 原状黄土
Abstract: The isotropic compression tests with simultaneous gas permeability measurements are performed for intact loess with various water contents in three regions. The behaviours of compression deformation and change of induced gas permeability of intact loess are analyzed. The correlation between gas permeability and compression deformation is discussed. The formula of compression curve of intact loess is proposed to characterize the compressive deformation under the coupling of moisture and stress, in which the structural parameter defined by gas permeability is included as a known parameter. The results suggest that changes of gas permeability of intact loess are found to have a good correlation with volumetric deformation during isotropic compression tests. The structural potential difference of intact loess caused by the fabric and moisture and its damage induced by the stress can be reflected by the change in gas permeability. The gas permeability ratio corresponding to the yield stress at the ratio of a certain water content to the reference water content is defined as the structureal ratio parameter, and it is introduced into the isotropic stress variable to lead to a normalized compression curve of intact loess with various water contents. Regarding the compression curve at the reference water content as the reference state, a unified and continued description of the compression deformation before and after yielding for unsaturated intact loess is proposed, and the calculated results agree well with the test ones.
- soil structure /
- compression /
- gas permeability /
- intact loess

HTML全文

0. 引言

地下水的原位修复是环境岩土领域的重要研究课题^[1]。近年来,以纳米材料作为修复剂的地下水原位修复技术得到了广泛的关注^[2-3],目前常用于污染物处理的纳米材料包括纳米零价铁、铁氧化物纳米颗粒、铁的双金属颗粒等^[4-5],其中纳米零价铁（nano zero-valent iron,nZVI）是最常见的一种。纳米零价铁在电镜下的单个颗粒粒径约为40～200 nm,颗粒内部为铁单质内核、表面为氧化物外壳,具有特殊的核壳结构^[6],对不同类型的污染物均具有修复能力^[7-8]。纳米零价铁的比表面积大、吸附能力强^[9],国内外的研究表明纳米零价铁对于水中的氯代烃、磷、铅、汞、砷、硝酸盐等多种污染物均有较好的去除效果^[9]。

纳米零价铁在多孔介质中的运移能力是该技术应用的关键问题^[7]。目前,通常将纳米零价铁注入被污染的地下含水层进行原位修复,因此纳米零价铁在多孔介质中的运移能力对其修复效果有显著的影响^[10]。之前的研究表明,未修饰的纳米零价铁在地下水中有较强的团聚趋势,从而抑制了其在地下系统中的运移,表面修饰等方法常被用来提高纳米零价铁的运移能力^[11]。地下水中的磷和砷等阴离子主要通过吸附在修复材料表面而去除,但负载在纳米零价铁表面的污染物仍然存在于地下系统之中,若其进一步在地下系统中运移,则有可能引发环境风险,因此有必要研究负载污染物的纳米零价铁的运移能力^[12-13]。地下水的环境化学条件对地下系统中颗粒的运移有重要的影响,如地下水的流速、离子强度等条件的改变将影响纳米零价铁的运移能力^[14-16]。此外,多孔介质特性也是影响颗粒运移的重要因素^[17],如多孔介质的粒径级配将影响孔隙结构及孔隙流场分布,进而影响纳米零零价铁的运移行为。因此,研究磷负载纳米零价铁在不同环境化学和多孔介质条件下的运移能力对于其风险分析和防控有着重要的意义。

纳米零价铁在水中的粒径随着其不同修饰方法而变化,其粒径通常处于10～10 μm^[18-19],属于胶体颗粒的粒径范围^[19]。研究表明,胶体在地下系统中通过对流和水动力弥散过程而运移,由于表面沉积（surface deposition）、颗粒阻塞（straining）以及粒桥（hydrodynamic bridging）等物理作用而滞留在多孔介质中^[20-21],如图1所示。表面沉积是指胶体颗粒由于附着在孔隙结构的表面而留存在多孔介质中,当颗粒靠近多孔介质表面时颗粒与多孔介质表面之间的存在范德华力和双电层力等力的作用,这些作用产生了使颗粒停留的力矩（adhesive torque）,同时孔隙流体作用在胶体颗粒上的拖曳力和升力则产生了使其运移的力矩（applied torque）,两个力矩的相对大小决定了该颗粒是否会发生表面沉积^[22]。颗粒阻塞是指当颗粒的粒径较大无法通过多孔介质中的狭小通道因而留存在多孔介质中,文献[23]称其为筛滤作用,Bradford等^[24]在研究中指出当胶体粒径与多孔介质粒径之比大于0.003时会发生颗粒阻塞。粒桥是指颗粒粒径小于孔隙通道尺寸时由于多个颗粒同时到达通道入口处而形成的类似拱的结构,该过程与颗粒粒径、团聚特性、流速等多种因素有关^[25]。此外,由于纳米零价铁较强的团聚趋势^[26],导致运移过程中其团粒粒径可能发生变化,使得其运移行为更为复杂^[18]。

图 1 多孔介质中胶体颗粒滞留机理

Figure 1. Retention mechanisms of colloids in porous media

下载: 全尺寸图片幻灯片

本文的研究目标是探讨磷负载对于纳米零价铁运移特性的影响,分析磷负载纳米零价铁在不同地下水及多孔介质条件下的运移能力。首先在实验室中合成了聚苯胺修饰的纳米零价铁（记为nZVI）和磷负载的聚苯胺修饰的纳米零价铁（记为PS-nZVI）,并对两种纳米零价铁的基本特性进行了分析,进而通过沉降试验研究了两种纳米零价铁在水中的胶体稳定性,最后通过柱试验对比研究了磷负载对nZVI运移能力的影响,并分析了流速、离子浓度、多孔介质粒径级配等因素对PS-nZVI运移能力的影响。

1. 材料

本文在实验室中合成了nZVI和PS-nZVI,通过电镜图片和平均粒径来表征其结构特性,测量zeta电位来反映其表面带电特性；用于模拟多孔介质的材料包括3种粒径的玻璃珠、天然砂以及高岭土。

1.1 纳米零价铁

采用液相还原法在实验室中合成纳米零价铁,并采用聚苯胺对纳米零价铁进行表面修饰^[27]。制备过程中,将6 g氯化铁溶于80 mL去离子水中,加入0.8 mL苯胺,以600 rpm的转速搅拌5 min,然后静置2 d,使其中的苯胺充分聚合形成聚苯胺。在氮气的保护下,将100 mL 1 mol/L的硼氢化钠溶液滴定到上述的聚苯胺-氯化铁混合液中,充分反应30 min后将反应容器中的黑色固体离心分离,冷冻干燥后得到黑色的nZVI粉末。

在去离子水中配置100 mL 150 mg/L的nZVI溶液,超声处理2 h后加入1 mL磷元素质量浓度为2 g/L的磷酸二氢钾溶液,在摇床上震荡60 min以上保证达到吸附平衡。之后用450 nm滤膜过滤,取滤膜上的固体部分重新溶于100 mL去离子水中,获得PS-nZVI,试验测得PS-nZVI的磷负载量为78.1 mg/g。

利用扫描电镜（SEM,Merlin,德国蔡司）和透射电镜（TEM,JEM-2010F,日本电子光学实验室）获取nZVI的形貌特征。分别在去离子水中配置150 mg/L的nZVI和PS-nZVI溶液,用激光粒度分析仪（Mastersizer 2000,英国马尔文公司）测定其平均粒径,并用表面电位仪（NanoBrook 90Plus Zeta,美国布鲁克海文公司）测定二者表面的zeta电位。

1.2 多孔介质材料

本文用于模拟多孔介质的玻璃珠、天然砂及高岭土的粒径级配曲线如图2所示,细玻璃珠、中玻璃珠、粗玻璃珠的平均粒径分别为0.11,0.57,1.50 mm,所用天然砂（ISO标准砂,厦门艾思欧标准砂有限公司）平均粒径为0.83 mm,高岭土的平均粒径为0.026 mm。

图 2 试验所用玻璃珠和天然砂的粒径级配曲线

Figure 2. Grain-size distribution curves of glass beads and sand in column tests

下载: 全尺寸图片幻灯片

2. 试验方法

本文采用沉降试验和柱试验来探究nZVI和PS-nZVI的团聚趋势以及运移特性。胶体在不同化学条件下的团聚趋势可以利用DLVO理论进行分析。

2.1 沉降试验

研究表明遮光度与颗粒悬浮在水中部分的浓度成线性关系,因此遮光度随时间的变化可以反应颗粒在水中的沉降过程^[28]。分别在去离子水中配置150 mg/L nZVI、150 mg/L PS-nZVI、75 mg/L PS-nZVI与75 mg/L高岭土混合溶液,超声处理2 h后利用分光光度计（DR5000,加拿大哈希公司）在508 nm波长下测定上述体系的遮光度随时间的变化,以研究nZVI和PS-nZVI在水中的稳定性,以及高岭土对于PS-nZVI稳定性的影响。

2.2 柱试验

一维柱试验装置由蠕动泵、柱和自动收集装置组成,其示意图如图3所示。柱内径3 cm,高15 cm,在上下各0.5 cm填充粒径为2.5 mm的大玻璃珠以确保形成一维条件,中间填充不同的多孔介质。试验过程中,将0.8个孔隙体积（柱中孔隙水总体积）的待研究液体分别从柱的顶部注入,随后持续注入5个孔隙体积的背景液,并收集柱底部的出流液体,用硝酸酸化后测量铁和磷的出流浓度,绘制出流曲线。试验结束后,将柱拆开后分层取出多孔介质,用硝酸酸化后测量沉积在多孔介质中的铁的含量,得到铁的滞留曲线。出流曲线和滞留曲线均除以注入浓度以获得相对液相浓度和相对固相浓度,本文中元素浓度采用电感耦合等离子体发射光谱法（ICP-OES,Prodigy7,美国黎曼实验室）测量。

图 3 柱试验装置示意图

Figure 3. Column test setup

下载: 全尺寸图片幻灯片

为研究磷负载对nZVI以及不同条件下PS-nZVI的运移特性,本文共开展10组柱试验,其中每组试验均进行了2次以上的重复试验,试验序号及条件如表1所示。表1中所述溶液均在去离子水中配置,并利用氯化钠调节其离子强度,填充过程中始终保持液面在多孔介质的上方以防止气泡留存,通过控制装填多孔介质的质量将柱的孔隙率控制在0.30～0.32,在填充完成后,通过10个孔隙体积的背景液以去除多孔介质及孔隙液中环境化学条件的影响。对于含有高岭土的试验组,将高岭土按照质量比5%与天然砂充分混合后进行填充,之后通过10个孔隙体积的背景液,在此过程中部分高岭土将随水流出,并最终达到稳定。

表 1 柱试验条件

Table 1. Conditions of column tests

序号	纳米零价铁	质量浓度/(mg·L^-1)	离子强度/(mmol·L^-1)	流速/(m·d^-1)	多孔介质
1	nZVI	150	0.1	16.3	中玻璃珠
2	PS-nZVI	150	0.1	16.3	中玻璃珠
3	PS-nZVI	150	0.1	4.1	中玻璃珠
4	PS-nZVI	150	1.0	4.1	中玻璃珠
5	PS-nZVI	150	10.0	4.1	中玻璃珠
6	PS-nZVI	150	0.1	4.1	细玻璃珠
7	PS-nZVI	150	0.1	4.1	粗玻璃珠
8	PS-nZVI	150	0.1	4.1	天然砂
9	PS-nZVI	150	0.1	16.3	天然砂
10	PS-nZVI	150	0.1	4.1	天然砂+5%高岭土

下载: 导出CSV

| 显示表格

2.3 DLVO理论

DLVO理论常用来评价胶体的稳定性^[29]。将胶体颗粒假设为理想球体,通过绘制两个球体之间的能量随距离的变化关系来分析胶体的稳定性,该能量随距离的变化曲线称为能量曲线,曲线中极大值点的高度称为势垒。通常势垒越高代表颗粒之间的斥力越大,胶体体系越稳定,而势垒越低则表明胶体有较强的团聚趋势。DLVO理论中总能量V_t计算的能量包括范德华力V_vdw、双电层力V_ES和波恩力V_B所带来的能量,其计算表达式如下^[30-33]：

$V_{t} = V_{vdw} + V_{ES} + V_{B}$ ,

(1)

$V_{vdw} = - \frac{A a}{12 s (1 + \frac{14 s}{λ})}$ ,

(2)

$V_{ES} = 2π ε_{0} ε_{r} a ψ_{1}^{2} \ln (1 + e^{- κ s})$ ,

(3)

$V_{B} = \frac{A σ_{c}^{6}}{7560} [\frac{8 a + h}{{(2 a + 7)}^{7}} + \frac{6 a - h}{h^{7}}]$ 。

(4)

式中,A为纳米零价铁的Hamaker常数（J）,取10^-19 J ^[28],a为纳米零价铁颗粒的半径（nm）,s为两个颗粒之间的距离（nm）,ε₀为真空介电常数（C²/J·m）,取8.85×10^-12 C²/J·m,ε_r为水的相对介电常数,大小为78.4,ψ₁纳米零价铁的表面电位（mV）,这里取为zeta电位, $σ_{c}$ 为碰撞直径（nm）,通常取为0.5 nm^[34],κ为Debye长度的倒数（m^-1）,按如下公式计算：

$κ = \sqrt{\frac{2 N_{A} e^{2} I}{ε_{r} ε_{0} k_{B} T}}$ ,

(5)

式中,N_A为Avogadro常数,取6.02×10²³ mol^-1,e为元电荷（C）,取1.6×10^-19 C,I为离子强度（mmol/L）,k_B为Boltzmann常数（J/K）,取1.38×10^-23 J/K,T为绝对温度（K）,取298 K。

3. 结果与讨论

3.1 nZVI团粒结构和zeta电位

nZVI的SEM和TEM图像如图4所示,可以看出nZVI主要以链状结构形式存在,单个粒径范围为50～150 nm。测得nZVI和PS-nZVI在水中的平均粒径为4.255,3.814 μm,这表明nZVI在水中团聚成直径在微米级的团粒,并且吸附磷对其平均粒径的影响不大。nZVI和PS-nZVI的zeta电位值分别为-10.5,-28.1 mV,这是由于纳米零价铁表面吸附的磷酸根离子带负电荷,使其表面电位显著降低。

图 4 nZVI团粒电镜图像

Figure 4. SEM and TEM images of nZVI aggregate

下载: 全尺寸图片幻灯片

3.2 沉降试验

nZVI、PS-nZVI以及PS-nZVI与高岭土的混合溶液相对遮光度随时间的变化如图5所示,结果表明nZVI在试验的前45 min快速聚沉,而后进入缓慢聚沉的过程,在8 h后相对遮光度降低了92%,而PS-nZVI相比于nZVI则具有较高的稳定性,在试验过程中始终保持着缓慢聚沉的特性,8 h后相对遮光度仅降低了37%。nZVI和PS-nZVI的稳定性可以用DLVO理论的能量曲线来解释,其计算结果如图6所示。从图6可以看出,吸附磷之后PS-nZVI的能量曲线的势垒显著提高,这是由于表面负电荷的增加引起的更高静电斥力所致,高势垒对应着高胶体稳定性,这与沉降试验中所得到的结果一致。

图 5 溶液相对遮光度随时间的变化

Figure 5. Relative optimal densities of suspensions over time

下载: 全尺寸图片幻灯片

图 6 nZVI与PS-nZVI的DLVO理论能量曲线

Figure 6. Interaction energy curves by DLVO theory for nZVI and PS-nZVI

下载: 全尺寸图片幻灯片

加入高岭土后,8 h后相对遮光度降低幅度达到84%,这表明在富含高岭土的环境下,PS-nZVI的聚沉将加快。Kim等^[35]的研究中曾提到,带电不均匀的高岭土颗粒可以在有多聚物修饰的纳米零价铁之间起到连接作用,从而促进纳米零价铁的团聚,这与本文针对PS-nZVI所开展的试验结论一致。

3.3 柱试验

（1）磷负载对nZVI运移的影响

nZVI和PS-nZVI的柱试验结果如图7所示,其中图7（a）,（b）分别展示了相对出流曲线和相对滞留曲线的结果。结果表明PS-nZVI的相对出流曲线显著高于nZVI,二者的液相回收率分别为69.6%和2.7%。二者在入口处沉积浓度均较高,这一现象对于nZVI来说尤为显著,这种颗粒在柱的入口处大量沉积的现象通常被称为超指数型滞留曲线（hyper- exponential retention profile）现象,研究者对该现象出现的原因提出了不同的解释^[36-39],如颗粒的粒桥以及颗粒阻塞作用,这两种机理均与颗粒的粒径密切相关。根据沉降试验的结论可知,nZVI在试验过程中伴随着剧烈的团聚和沉降,因此在颗粒穿过多孔介质的过程中伴随着粒径的增加,颗粒的粒桥和颗粒阻塞更为显著；而PS-nZVI的粒径分布则相对较为稳定,能够在多孔介质中运移更长的距离。应用纳米零价铁进行污染场地原位修复时,有必要了解负载污染物的纳米零价铁的运移能力以及环境因素的影响,从而合理评估和控制环境风险。

图 7 nZVI与PS-nZVI柱试验结果

Figure 7. Results of column tests of nZVI and PS-nZVI

下载: 全尺寸图片幻灯片

（2）地下水条件对PS-nZVI运移的影响

图8（a）,（b）分别展示了流速对PS-nZVI的相对出流曲线和相对滞留曲线的影响。结果表明随着流速的降低,PS-nZVI的运移能力减弱,其液相出流浓度降低至50.3%,滞留曲线的浓度也有所上升。流速对颗粒出流能力的影响主要包括多个方面：一方面是由于流速的提升,颗粒随水通过孔隙的速度加快,从而使得对于相同体积的注入溶液,颗粒在多孔介质中存在的时间缩短,颗粒粒径的变化也较小,因此抑制了粒桥和颗粒阻塞的发生；另一方面高流速对于靠近多孔介质表面的颗粒有更强的水动力作用^[40]。陈星欣等^[23]的试验研究表明在颗粒粒径一定的条件下,增加流速会导致回收率的增加；Bradford等^[41]的研究表明,随着流速的增加,水流作用在颗粒上的力矩增加,该力矩有助于颗粒离开多孔介质表面,因此能够减少颗粒的表面沉积；Kermani等^[42]在孔隙尺度进行了数值模拟,结果表明高流速条件下颗粒的团聚作用和表面沉积作用都有所降低。

图 8 PS-nZVI在不同流速条件下的柱试验结果

Figure 8. Results of column tests of PS-nZVI under different velocity conditions

下载: 全尺寸图片幻灯片

图9（a）展示了不同离子强度条件下PS-nZVI的相对出流曲线,对应情况的能量曲线如图9（b）所示。结果表明随着离子强度的提高,PS-nZVI的运移能力减弱,随着离子强度从低到高试验,液相回收率分别为50.3%,22.2%,17.2%。离子强度可以影响颗粒与颗粒之间的团聚特性,随着离子强度的升高,颗粒表面的双电层厚度降低,从而导致了势垒的降低。因此,高离子强度条件更有利于颗粒的团聚,从而使得颗粒的粒径增加^[43],有利于粒桥和颗粒阻塞的产生,nZVI滞留量提高。同时,高离子强度条件下,颗粒与多孔介质表面的相互作用更强,更有利于通过表面沉积的方式滞留在多孔介质中^[16]。

图 9 PS-nZVI在不同离子强度条件下的柱试验结果

Figure 9. Results of column tests of PS-nZVI under different ionic strength conditions

下载: 全尺寸图片幻灯片

（3）多孔介质粒径级配对PS-nZVI运移的影响

多孔介质粒径级配对于PS-nZVI运移能力的影响如图10所示。图10（a）表明PS-nZVI在3种不同粒径玻璃珠中的运移能力呈现出中玻璃珠>粗玻璃珠>细玻璃珠的趋势,三者中的液相回收率分别为50.3%,41.0%,4.1%。多孔介质的孔隙尺寸随着粒径的增加而增大,粒桥和颗粒阻塞随着孔隙的增大而减弱,因此随着多孔介质粒径的增加,粒桥和颗粒阻塞所导致的颗粒滞留量将减少。然而,Bradford等^[41]的数值模拟工作表明,当多孔介质的粒径增加且表观流速一致时,胶体颗粒在多孔介质表面所受到的力矩随着多孔介质粒径的增加而降低,这是由于流体作用在胶体颗粒表面的拖曳力与流速梯度成正比。当表观流速一致时,大颗粒中的流速梯度较低,因此在大颗粒表面更容易发生表面沉积。这一现象表明,颗粒运移能力随着多孔介质粒径的变化将呈现出不同的趋势,当以粒桥或颗粒阻塞等与颗粒粒径相关的滞留作用为主时,运移能力随着多孔介质粒径的增加而增大；当以表面沉积的滞留作用为主时,运移能力则随着多孔介质粒径的增加而降低。本研究中PS-nZVI在细玻璃珠中出流能力很低是由于粒桥和颗粒阻塞所致,而在粗玻璃珠中出流能力低于中玻璃珠则是主要是由于表面沉积作用所致。

图 10 PS-nZVI在不同多孔介质中的相对出流曲线

Figure 10. Relative breakthrough curves of PS-nZVI in different porous media

下载: 全尺寸图片幻灯片

PS-nZVI在天然砂中的出流曲线如图10（b）所示。天然砂的平均粒径接近中玻璃珠（介于中玻璃珠和粗玻璃珠之间）,PS-nZVI在天然砂中的运移能力则低于玻璃珠,这是由于天然砂的粒径范围较大,含有的小颗粒起到了填充大孔隙的作用,导致更多的纳米零价铁滞留在小孔隙中^[44]。同时,天然砂中PS-nZVI的出流早于玻璃珠,这反映了PS-nZVI在天然砂中的水动力弥散作用高于玻璃珠,主要原因是玻璃珠粒径相对均匀,而天然砂粒径级配分布较广,孔隙结构更复杂。含高岭土的天然砂中,PS-nZVI的运移能力很弱,结合图5所示的沉降试验结果可知,高岭土能够促进PS-nZVI的团聚,在实际工程中,PS-nZVI可能会富集在含有高岭土的土层中。

4. 结论

本文在实验室中合成了nZVI和PS-nZVI,通过沉降试验和一维柱试验研究了磷负载对于nZVI运移能力的影响,并探究了PS-nZVI在不同流速、离子强度、多孔介质粒径级配等条件下的运移能力,得到以下4点结论。

（1）表面电位对于纳米零价铁的运移能力有显著影响。磷负载后nZVI的表面负电荷增加进而导致颗粒团聚趋势降低,团聚趋势的变化可以用DLVO理论解释。PS-nZVI在多孔介质中的运移能力显著高于nZVI,这可能引发一定的环境风险。

（2）地下水条件对PS-nZVI的运移有显著的影响。随着离子浓度的升高,PS-nZVI的运移能力减弱,这是由于高离子强度有利于颗粒的团聚,同时增加颗粒与多孔介质表面的相互作用；出流能力随着流速的提高而增加,这是由于流速的提升缩短了颗粒团聚的时间,并提高了颗粒在流场中所受到的力矩。

（3）多孔介质的粒径级配对纳米零价铁的运移能力有显著影响。PS-nZVI在不同粒径玻璃珠中的运移能力有显著区别（中玻璃珠>粗玻璃珠>细玻璃珠）,这是由于PS-nZVI的滞留过程包括多种作用机理,其中粒桥、颗粒阻塞随着多孔介质粒径的增加而减弱,表面沉积随着多孔介质粒径的增加而增强。

（4）多孔介质中细颗粒将影响纳米零价铁的运移行为。天然砂中的细颗粒起到了填充孔隙的作用,PS-nZVI在天然砂中的运移能力低于在平均粒度接近的玻璃珠中的运移能力；天然砂中复杂孔隙结构增强了弥散作用,PS-nZVI的出流时间早于玻璃珠；高岭土的存在能够促进PS-nZVI的团聚和沉积,说明PS-nZVI可能富集在含有高岭土的土层中。

图 1 原状黄土的颗粒级配曲线

Figure 1. Grain-size distribution curves of intact loess

下载: 全尺寸图片幻灯片

图 2 三轴渗气仪

Figure 2. Triaxial equipment for air permeability

下载: 全尺寸图片幻灯片

图 3 原状黄土的渗气速率与气压梯度的典型关系曲线

Figure 3. Relationship between rate of gas flow and gas-pressure gradient at various given isotropic stresses for intact loess

下载: 全尺寸图片幻灯片

图 4 不同含水率下原状黄土的e–lnp和k_a–lnp关系曲线

Figure 4. Curves of e-lnp and k_a-lnp of intact loess with various water contents

下载: 全尺寸图片幻灯片

图 5 不同含水率原状黄土结构屈服应力比 $χ_{0}$ 与结构比参数 ${χ^{'}}_{0}$ 之间关系

Figure 5. Relationship between ratio of yield stress and structural ratio parameter of intact loess with various water contents

下载: 全尺寸图片幻灯片

图 6 结构比参数对不同区域不同含水率原状黄土压缩变形的归一化曲线

Figure 6. Normalized curves of compressive deformation of intact loess with various water contents in different areas by using structural ratio parameter

下载: 全尺寸图片幻灯片

图 7 非饱和原状黄土的e–lnp压缩试验与计算结果

Figure 7. Comparison between calculated and measured results from isotropic compression tests at various water contents for intact loess

下载: 全尺寸图片幻灯片

表 1 原状黄土的物理性质

Table 1 Physical properties of intact loess

地区	含水率w_n/%	干密度ρ_d/(g·cm^-3)	液限w_L/%	塑限w_P/%	塑性指数I_p	土样定名
西安	15.2	1.30	30.9	19.8	11.1	粉质黏土
天水	11.0	1.45	27.5	17.9	9.6	粉土
兰州	8.0	1.30	24.9	16.5	8.4	粉土

下载: 导出CSV

表 2 渗气试验的控制条件与试验结果

Table 2 Control conditions and results of air permeability tests

地区	初始含水率w₀/%	等向应力p/kPa						屈服应力p_y/kPa	渗气系数k_ay/(cm·s^-1)
西安	1.5	0	50	100	200	300	400	265	6.68×10^-2
	8.5	0	50	100	200	300	400	210	4.61×10^-2
	13.5	0	50	100	200	300	400	155	2.47×10^-2
	15.2	0	50	100	200	300	400	135	2.14×10^-2
	16.6	0	50	100	200	300	400	115	1.12×10^-2
	21.8	0	50	100	200	300	400	100	8.40×10^-3
天水	11.0	－	50	100	200	－	400	240	9.50×10^-4
	13.0	－	50	100	200	－	400	180	7.44×10^-4
	15.0	－	50	100	200	－	400	145	5.34×10^-4
	19.0	－	50	100	200	－	400	105	3.51×10^-4
兰州	8.0	0	50	120	200	300	－	110	1.02×10^-2
	10.5	0	50	120	200	300	－	105	8.20×10^-3
	12.5	0	50	120	200	300	－	92	5.60×10^-3
	14.5	0	50	120	200	300	－	75	2.56×10^-3
	18.7	0	50	120	200	300	－	68	7.19×10^-4
注：“－”表示未在对应的等向应力下进行渗气试验。

下载: 导出CSV

参考文献(30)

[1]	沈珠江. 土体结构性的数学模型——21 世纪土力学的核心问题[J]. 岩土工程学报, 1996, 18(1): 95-97. doi: 10.3321/j.issn:1000-4548.1996.01.015 SHEN Zhu-jiang. Mathematical model for soil structure—The core topic of soil mechanics in the 21st century[J]. Chinese Journal of Geotechnical Engineering, 1996, 18(1): 95-97. (in Chinese) doi: 10.3321/j.issn:1000-4548.1996.01.015
[2]	谢定义, 齐吉琳. 土结构性及其定量化参数研究的新途径[J]. 岩土工程学报, 1999, 21(6): 651-656. doi: 10.3321/j.issn:1000-4548.1999.06.003 XIE Ding-yi, Qi Ji-lin. Soil structure characteristics and new approach in research on its quantitative parameter[J]. Chinese Journal of Geotechnical Engineering, 1999, 21(6): 651-656. (in Chinese) doi: 10.3321/j.issn:1000-4548.1999.06.003
[3]	JIANG M J, ZHANG F F, HU H J, et al. Structural characterization of natural loess and remolded loess under triaxial tests[J]. Engineering Geology, 2014, 181: 249-260. doi: 10.1016/j.enggeo.2014.07.021
[4]	高国瑞. 黄土湿陷变形的结构理论[J]. 岩土工程学报, 1990, 12(4): 1-10. doi: 10.3321/j.issn:1000-4548.1990.04.001 GAO Guo-rui. A structure theory for collapsing deformation loess soils[J]. Chinese Journal of Geotechnical Engineering, 1990, 12(4): 1-10. (in Chinese) doi: 10.3321/j.issn:1000-4548.1990.04.001
[5]	雷祥义. 中国黄土的孔隙类型与湿陷性[J]. 中国科学(B辑), 1987, 17(12): 1309-1318. https://www.cnki.com.cn/Article/CJFDTOTAL-JBXK198712008.htm LEI Xiang-yi. Classification of micro-structure and collapsibility of the loess[J]. Scientia Sinica, 1980, 17(12): 1309-1318. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JBXK198712008.htm
[6]	WANG J D, LI P, MA Y, et al. Change in pore-size distribution of collapsible loess due to loading and inundating[J]. Acta Geotechnica, 2020, 15(5): 1081-1094. doi: 10.1007/s11440-019-00815-9
[7]	沈珠江, 胡再强. 黄土的二元介质模型[J]. 水利学报, 2003, 34(7): 1-6. doi: 10.3321/j.issn:0559-9350.2003.07.001 SHEN Zhu-jiang, HU Zai-qiang. Binary medium model for loess[J]. Journal of Hydraulic Engineering, 2003, 34(7): 1-6. (in Chinese) doi: 10.3321/j.issn:0559-9350.2003.07.001
[8]	李宏儒, 胡再强, 赵凯, 等. 结构性土二元介质本构模型及破损率影响因素的研究[J]. 岩土力学, 2012, 33(增刊1): 67-72. doi: 10.16285/j.rsm.2012.s1.014 LI Hong-ru, HU Zai-qiang, ZHAO Kai, et al. Structural soil binary-medium constitutive model and factor of breakage ratio influence[J]. Rock and Soil Mechanics, 2012, 33(S1): 67-72. (in Chinese) doi: 10.16285/j.rsm.2012.s1.014
[9]	陈正汉. 非饱和土与特殊土力学的基本理论研究[J]. 岩土工程学报, 2014, 36(2): 201-272. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201402002.htm CHEN Zheng-han. On basic theories of unsaturated soils and special soils[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(2): 201-272. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201402002.htm
[10]	邵帅, 褚峰, 邵生俊. Q₃结构性黄土的扰动状态本构模型试验研究[J]. 岩石力学与工程学报, 2016, 35(7): 1494-1500. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201607020.htm SHAO Shuai, CHU Feng, SHAO Sheng-jun. Experimental study on constitutive model of structural Q₃ loess based on disturbed state concept[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(7): 1494-1500. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201607020.htm
[11]	姚志华, 连杰, 陈正汉, 等. 考虑细观结构演化的非饱和Q₃原状黄土弹塑性本构模型[J]. 岩土力学, 2018, 39(5): 1553-1563. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201805002.htm YAO Zhi-hua, LIAN Jie, CHEN Zheng-han, et al. An elastic plastic constitutive model for unsaturated Q₃ undisturbed loess considering meso-structured evolution[J]. Rock and Soil Mechanics, 2018, 39(5): 1553-1563. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201805002.htm
[12]	邵生俊, 周飞飞, 龙吉勇. 原状黄土结构性及其定量化参数研究[J]. 岩土工程学报, 2004, 26(4): 531-536. doi: 10.3321/j.issn:1000-4548.2004.04.021 SHAO Sheng-jun, ZHOU Fei-fei, LONG Ji-yong. Structural properties of loess and its quantitative parameter[J]. Chinese Journal of Geotechnical Engineering, 2004, 26(4): 531-536. (in Chinese) doi: 10.3321/j.issn:1000-4548.2004.04.021
[13]	骆亚生, 谢定义, 邵生俊, 等. 复杂应力状态下的土结构性参数[J]. 岩石力学与工程学报, 2004, 23(24): 4248-4251. doi: 10.3321/j.issn:1000-6915.2004.24.027 LUO Ya-sheng, XIE Ding-yi, SHAO Sheng-jun, et al. Structural parameter of soil under complex stress conditions[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(24): 4248-4251. (in Chinese) doi: 10.3321/j.issn:1000-6915.2004.24.027
[14]	陈存礼, 高鹏, 胡再强. 黄土的增湿变形特性及其与结构性的关系[J]. 岩石力学与工程学报, 2006, 25(7): 1352-1360. doi: 10.3321/j.issn:1000-6915.2006.07.009 CHEN Cun-li, GAO Peng, HU Zai-qiang. Moistening deformation characteristic of loess and its relation to structure[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(7): 1352-1360. (in Chinese) doi: 10.3321/j.issn:1000-6915.2006.07.009
[15]	冯志焱. 非饱和黄土的结构性定量化参数与结构性本构关系研究[D]. 西安: 西安理工大学, 2008. FENG Zhi-yan. Research on Soil Structural Parameters and Structural Constitutive Model of Unsaturated Loess[D]. Xi'an: Xi'an University of Technology, 2008. (in Chinese)
[16]	刘奉银, 张昭, 周冬. 湿度和密度双变化条件下的非饱和黄土渗气渗水函数[J]. 岩石力学与工程学报, 2010, 29(9): 1907-1914. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201009022.htm LIU Feng-yin, ZHANG Zhao, ZHOU Dong. Density-saturation-dependent air-water permeability function of unsaturated loess[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(9): 1907-1914. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201009022.htm
[17]	CHEN C L, ZHANG D F, ZHANG J. Influence of stress and water content on air permeability of intact loess[J]. Canadian Geotechnical Journal, 2017, 54(9): 1221-1230. doi: 10.1139/cgj-2016-0186
[18]	PICCOLI I, SCHJØNNING P, LAMANDÉ M, et al. Coupling gas transport measurements and X-ray tomography scans for multiscale analysis in silty soils[J]. Geoderma, 2019, 338(3): 576-584.
[19]	ZHAI X F, HORN R. Dynamics of pore functions and gas transport parameters in artificially ameliorated soils due to static and cyclic loading[J]. Geoderma, 2019, 337(3): 300-310.
[20]	姚志华, 陈正汉, 黄雪峰, 等. 非饱和Q₃黄土渗气特性试验研究[J]. 岩石力学与工程学报, 2012, 31(6): 1264-1273. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201206024.htm YAO Zhi-hua, CHEN Zheng-han, HUANG Xue-feng, 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) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201206024.htm
[21]	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
[22]	BERISSO F E, SCHJONNING P, KELLER T, et al. Gas transport and subsoil pore characteristics: anisotropy and long-term effects of compaction[J]. Geoderma, 2013, 195/196: 184-191.
[23]	袁俊平, 张锋, 王启贵, 等. 裂隙对压实膨胀土渗气性影响试验[J]. 水利水电科技进展, 2014, 34(3): 34-38. https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD201403009.htm YUAN Jun-ping, ZHANG Feng, WANG Qi-gui, et al. The influences of fissure on gas permeability of the compacted expansive soil[J]. Advances in Science and Technology of Water Resources, 2014, 34(3): 34-38. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SLSD201403009.htm
[24]	NGUYEN V, PINEDA J A, ROMERO E, et al. Influence of soil microstructure on air permeability in compacted clay[J]. Géotechnique, 2021, 71(5): 373-391.
[25]	李喜安, 刘锦阳, 郭泽泽, 等. 马兰黄土孔隙结构参数与渗透性关系研究[J]. 工程地质学报, 2018, 26(6): 1415-1423. https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201806002.htm LI Xi-an, LIU Jin-yang, GUO Ze-ze, et al. Study on relationship between pore structure parameters and permeability of malan loess[J]. Journal of Engineering Geology, 2018, 26(6): 1415-1423. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ201806002.htm
[26]	陈存礼, 张登飞, 张洁, 等. 等向应力条件下原状Q₃黄土的渗气特性研究[J]. 岩土工程学报, 2017, 39(2): 287-294. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201702015.htm CHEN Cun-li, ZHANG Deng-fei, 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) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201702015.htm
[27]	陈存礼, 曹程明, 王晋婷, 等. 湿载耦合条件下结构性黄土的压缩变形模式研究[J]. 岩土力学, 2010, 31(1): 39-45. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201001008.htm CHEN Cun-li, CAO Cheng-ming, WANG Jin-ting, et al. Modeling compression of structural loess under coupling of stress and moisture[J]. Rock and Soil Mechanics, 2010, 31(1): 39-45. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201001008.htm
[28]	LIU Z B, XIE S Y, SHAO J F, et al. Multi-step triaxial compressive creep behaviour and induced gas permeability change of clay-rich rock[J]. Géotechnique, 2018, 68(4): 281-289.
[29]	尹振宇. 天然软黏土的弹黏塑性本构模型：进展及发展[J]. 岩土工程学报, 2011, 33(9): 1357-1369. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201109011.htm YIN Zhen-yu. Elastic viscoplastic models for natural soft clay: review and development, 2011, 33(9): 1357-1369. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201109011.htm
[30]	姚仰平, 刘林, 罗汀. 砂土的UH模型[J]. 岩土工程学报, 2016, 38(12): 2147-2153. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201612004.htm YAO Yang-ping, LIU Lin, LUO Ting. UH model for sands[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(12): 2147-2153. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201612004.htm

施引文献(9)

期刊类型引用(4)

1.	邹新军，杨紫健，吴文兵. 非饱和土地基中端承桩对SH波的水平地震响应. 岩土工程学报. 2024(01): 72-80 . 本站查看
2.	林浩，郑长杰，丁选明. 层状地基中海洋大直径管桩水平动力响应分析. 岩土力学. 2024(06): 1873-1883 . 百度学术
3.	胡枫，宋春雨，时昆鹏，李得睿，陈龙珠. 砂土中桩基动力响应振动台模型试验. 实验室研究与探索. 2023(01): 79-84 . 百度学术
4.	黄娟，胡钟伟，余俊，李东凯. 考虑黏性的液化土中水平振动桩基桩顶阻抗研究. 岩土力学. 2023(05): 1445-1456 . 百度学术