土层结构对非饱和毛细水盐运移的影响

2021-08-26吕擎峰单小康赵彦旭贾博博

中国地质灾害与防治学报 2021年4期

吕擎峰，单小康，赵彦旭，贾博博，杲斐

（1.兰州大学土木工程与力学学院，甘肃兰州 730000；2.中国电建集团中南勘测设计研究院有限公司，湖南长沙 410000；3.中国铁建二十一局集团有限公司，甘肃兰州 730050）

0 引言

土体沉积受搬运介质和搬运方式的影响，往往形成不同土层结构的层状土体。层状土层之间孔隙结构和水力学性质的不连续性[1−3]，导致层状土体和均质土体的水盐运移规律存在显著差异。

层状土体的水盐运移规律在土壤学领域已经取得了一些研究成果。李韵珠等[4]指出壤土中的黏土夹层对水盐运移具有阻滞作用。史文娟等[5]通过室内土柱试验发现不同层位的夹砂层对水盐运移存在抑制或促进作用。HUANG 等[6]发现砂土和砂黄土互层结构可以限制水分向上运移，抑制土体蒸发。李毅等[7]认为砂土夹层和黏土夹层均有阻水作用。黏土的阻水能力源于其低渗透性，砂土层的阻水能力则是因为孔隙水由细粒土层向砂土层运移时存在滞后效应，因此在层状土层界面处水分运移均会出现停滞[8−9]。张莉等[10]研究发现在壤土中设置砂层能够加速上覆土层脱盐和抑制下伏土层返盐，该方法较依赖于降雨的淋滤作用，对干旱半干旱地区的适用性还有待研究。ALIMI 等[11]设计了冲积土和黏质砂土以及冲积土和黏土的双层土柱试验，结果表明黏土层的阻渗能力更强，而黏土层厚是延迟溶质运移的关键因素。LIU 等[12]对不同土体的水分和溶质含量进行了长期监测，发现均质土中水盐入渗较快，而层状土中水流受阻，盐分在局部累积。

层状土层水盐运移规律研究不仅与土壤水盐管理相关，同时涉及地基承载力计算、隐伏岩溶塌陷和路基盐胀等工程问题[13−16]。目前层状土层水盐运移研究侧重于讨论入渗过程的水盐动态变化，而针对层状土层毛细水盐运移研究尚且不足。本文基于某高铁路基上拱变形处的现场调查，设计了两种不同粒径土层二元结构组合（黄土-砂质粉土和黄土-粉质黏土）的室内土柱模型试验，讨论了不同地下水补给条件下土层结构对土体毛细水分布和盐分累积的影响。研究结果旨在为层状土区硫酸盐病害防治提供试验依据。

1 试验背景

西北地区某高铁路基上拱变形处的地层剖面和水盐分布如图1 所示，表层覆盖薄层砾石土，上覆粉质砂土层，层厚约80 cm，下伏厚层圆砾石土层。粉质砂土层的含水率和离子含量均高于下伏圆砾土层，40 cm 层位处盐分聚积，SO24−和 Cl−含量接近2%，下伏圆砾石土层离子含量则非常低。由于区域内地下水埋深较大（约9 m），而且粗粒土中毛细上升高度较小，距地表400 cm 土层中基本没有地下毛细水补给，驱动盐分运移的水分主要来源于降雨入渗和地表径流。

图1 上拱变形处土层结构与水盐垂直分布图Fig.1 The soil layer structure and water-salt distribution in heave site

由图2 可知，未上供变形处剖面上覆砂砾土层，下伏粉质黏土层，含水率和 Cl−含量在界面处呈断层式分布。上覆粗粒土层含水率约0.5%～2%， Cl−含量均小于0.08%；下伏细粒土层含水率在9%～13%之间， Cl−含量均大于0.57%。 S O24−在深度120～400 cm 土层中含量较高，深度0～80 cm 土层中含量较低。和 Cl−最大含量均小于1%，盐分未发生局部聚积，且表层含盐量较低。下伏粉质黏土层的毛细上升高度较大，有利于地下毛细水补给，通过下伏土层含水率较大推测距地表400 cm土层内有一定的毛细水补给。

图2 未上拱变形处土层结构与水盐垂直分布图Fig.2 The soil layer structure and water-salt distribution in contrast site

上拱变形处土层呈上细下粗型结构，盐分在上覆细粒土层中聚积，而未上拱变形处土层呈上粗下细型结构，盐分主要储存在下伏细粒土层中，离子分布较为均一，表层土体未发生盐分聚积现象。上述现象表明土层结构对水盐分布存在显著影响，而水盐分布的差异性对硫酸盐病害发育存在潜在影响。

2 试验设计与方法

2.1 试验设计

根据现场调查情况设计了室内土柱试验，其装置如图3 所示，由多层有机玻璃筒和补水瓶组成。有机玻璃筒内径20 cm，高度分为40 cm 和10 cm 两种，中间用法兰盘连接。模拟试验土柱上覆土层厚30 cm，下伏土层厚10 cm，现场调查发现区域内地下水含盐量较低，盐分主要分布在土层中，因此将下伏土层设置为质量分数5%的含盐土层。含盐土层下方根据不同的地下水条件设置不同隔断层，标准砂层模拟有毛细水补给的情况，卵砾石层则模拟无毛细水补给的情况，隔离层的另一作用是防止含盐层的盐分扩散至蒸馏水中。

图3 试验装置Fig.3 Test apparatus

试验组设计如表1 所示。上覆土层均为洗盐后的黄土，下伏土层分别为砂质粉土和粉质黏土，黄土-砂质粉土构成上细下粗型土层结构，黄土-粉质黏土构成上粗下细型土层结构，土柱ΙA 和ΙB 中存在毛细水补给，ⅡA 和ⅡB 中则无毛细水补给。

表1 试验组设计Table 1 Test group design

试验黄土取自甘肃省兰州市东岗镇，砂质粉土与粉质黏土由洗盐后的黄土掺石英砂和黄土掺高岭土得到。三种土样的基本参数见表2，颗分曲线见图4。

表2 试验土样的物理性质Table 2 Physical properties of soil samples

图4 试验土样的粒径分布曲线Fig.4 Grain-size distribution curves of soil samples

土样以分层夯实的方式填入有机玻璃筒中，上覆黄土层采用15%最优含水率和90%的压实度，即干密度1.62 g/cm3，下伏土层的含水率与上覆土层保持一致，考虑到不同土样的天然孔隙比存在显著差异，因此干密度通过统一的压实度进行控制，压实度均设置为90%，砂质粉土层夯实干密度为1.8 g/cm3，粉质黏土层夯实干密度为1.56 g/cm3。

2.2 试验过程

试样制作完成后在室内进行蒸发试验，分别在10，20，40，60 d 取样，前三次通过筒壁上的预留取样孔进行取样，取样深度为0，7，14，21，28，32，36，40 cm；60 d 时拆除土柱取样，并增加2，4 cm 取样深度，土样采用烘干法进行含水率测定。取10 g 干土样按照土水比1∶5 配成土水混合溶液，在 SHZ–88 型振荡器上振荡 30 min，将震荡后的土悬液放入TDL-SA 型离心机中离心30 min，用针管抽出离心管上部清液，经过 0.45 µm 孔径微孔膜过滤，得到土样浸出液，采用DIONEX 公司生产的ICS–2500 型离子色谱仪测定60 d 土样浸出液中的易溶盐离子含量。

3 试验结果与分析

3.1 毛细水补给土柱的水分运移

毛细水补给土柱ΙA 和ΙB 含水率沿深度分布结果见图5，两组试验土柱含水率在10 d 时已经基本稳定，试验周期内同一层位的含水率略有波动，整体分布趋势没有显著变化。土柱ΙA 和ΙB 上覆黄土层的含水率随深度增大整体呈线性增大，拟合线性曲线的斜率分别为4.94 和4.93，含水率梯度基本一致。ΙA 土柱上覆土层中的各层位含水率均大于ΙB，非饱和导水率K(θ)和扩散率D(θ)均为含水率的连续增函数，因此ΙA 土柱上覆土层中的导水率和扩散率均大于IB，IA 土柱毛细水运移速率较大。含水率较大时，砂质粉土层的导水率远大于粉质黏土，毛细水通过砂质粉土层向黄土土层运移相对较为容易；而粉质黏土层孔隙小，导水率低，水分运移速率缓慢，对土柱整体的水分运移存在阻滞作用。

图5 毛细水补给土柱含水率沿深度分布Fig.5 Distribution of water content along depthwith capillary water supply

IA 和IB 土柱含水率在土层界面处发生突变，ΙA 土柱界面处32 cm 与28 cm 含水率沿深度减小方向增加8.1%～9.1%，ΙB 土柱含水率减少3.6%～4.2%，界面含水率变化存在差异。含水率的突变是土层孔隙结构的不连续性导致的，从吸力角度考虑，当土柱的含水率达到平衡状态时，土体吸力与孔隙水势能具有连续性，30 cm 处上下界面的吸力大小一致；由于不同土层的颗粒大小、孔隙分布存在差异，同一土体吸力对应的不同土层的体积含水率存在一特定差值，这一差值在土体中表现为土层界面含水率突增或突降。

为了更好的分析和解释界面含水率的变化差异，通过Arya-Paris 模型[17]和Fredlund-Xing 模型[18]拟合了试验中三种土层的土水特征曲线。将土层的颗粒分布、干密度和颗粒密度代入AP 模型预测出多个吸力与体积含水率对应点，再由FX 模型拟合得到完整的土水特征曲线（图6）。

图6 三种土层的土水特征曲线Fig.6 Soil-water characteristic curves of three kinds of soil layers

同一吸力对应的三种土层体积含水率不同，细粒土的体积含水率均大于粗粒土的体积含水率。土层之间的含水率差值除了与土层性质有关，还随吸力变化。低吸力段，黄土-砂质粉土的含水率差值较大；高吸力段，黄土-粉质黏土的含水率差值较大。IA 和IB 土柱土层界面处的含水率较大，吸力较小，因此IA 土柱的界面含水率差值大于IB。

3.2 毛细水补给土柱的盐分运移

60 d 时土柱主要离子垂直分布情况如图7，ΙA 和ΙB 离子整体分布均呈 “Γ”型，盐分在表层聚积。土柱中的初始值较大，表层含量大于 N a+，ΙA 土柱表层含量高达15.64%，ΙB 土柱表层含量约7.69%，仅为ΙA 土柱的1/2，IB 土柱其他层位的含量均高于ΙA。 N a+分布情况与一致，ΙA 土柱表层Na+含量大于ΙB，其他层位则小于IB。两组土柱和 N a+的分布差异表明在相同的时间内，ΙA 土柱向上的盐分迁移更为剧烈，表层积盐更多，上细下粗型土层结构相对于上粗下细型结构更有利于土体盐分运移和表层盐分积累。

图7 毛细水补给土柱60 d 离子分布情况Fig.7 Ions distribution on the 60th day with capillary water supply

土层结构对水分运移的影响直接导致了盐分运移的差异，下伏砂质粉土层的存在加速了土柱的毛细水盐运移，导水率较低的粉质黏土层则阻碍了毛细水盐运移。除了对水分运移的影响，黏土颗粒较强的吸附性也延缓了土体中的溶质运移，黏粘土颗粒比表面积较大，结合水膜较厚，结合水膜中的离子需要先扩散至毛细水中，后随着毛细水迁移，这一过程降低了含盐层的脱盐速率，表层盐分积累速率相应减小。通过60 d 表层的泛盐现象可以直观的看出两组土柱的积盐差异（图8），ΙA 土柱表层的盐壳厚度较大，表层积盐较多，IB 土柱表层盐壳则较薄，盐结晶较少。由于硫酸盐结晶膨胀，土柱表层土体隆起，表面发生盐胀破坏。

图8 毛细水补给土柱60 d 表层积盐情况Fig.8 Surface salt accumulation on the 60th day with capillary water supply

3.3 无毛细水补给土柱的水分运移

土柱ⅡA 和ⅡB 中设置了卵砾石层隔断毛细水补给，模拟了地下水埋深较大或毛细水上升高度较小条件下的水盐运移，其含水率结果见图9。ⅡA 和ⅡB 土柱随着试验进行，ⅡA 和ⅡB 土柱整体含水率逐渐降低，在40～60 d 时，含水率变化非常小，表明土柱中的水分散失速率减小。一方面是因为随着表层一定深度土层含水率的减小，蒸发速率降低；另一方面则由于孔隙水含量的减小，毛细作用减弱，孔隙水主要受短程吸附作用，以薄膜水的形式附着在土颗粒表面，水分运移速率大幅降低，因此土柱内部的水分散失速率减小。

图9 无毛细水补给土柱含水率沿深度分布曲线Fig.9 Distribution of water content along depth without capillary water supply

土柱ⅡA 和ⅡB 深度30 cm 土层界面处的含水率变化趋势与毛细水补给土柱一致，但含水率差值大小存在显著差异。ⅡA 土柱界面含水率差值保持在1.68%～2.01%之间，ⅡB 土柱保持在2.78%～3.14% 之间。由图6 已知在高吸力状态下，黄土-砂质粉土土层界面含水率差值小于黄土-粉质黏土，因此ⅡB 土柱界面含水率差值大于ⅡA。砂质粉土层的水分向上覆黄土层迁移相对较多，在10～40 d 内，ⅡA 土柱上覆黄土层的整体含水率以及表层含水率明显大于ⅡB，在40～60 d 时，ⅡA 土柱下伏土层的含水率变化极小，上下土层之间的水分运移基本停止，而ⅡB 土柱下伏土层的含水率在该时间段内仍在减小，说明ⅡA 土柱水分运移速率明显大于ⅡB，主要原因是孔隙结构变化导致土水性质的差异，其次则因为黏土颗粒的短程吸附作用较强，结合水膜较厚，粉质黏土层中的薄膜水含量相对于砂质粉土层较多，而毛细水含量较少。

3.4 无毛细水补给土柱的盐分运移

无毛细水补给土柱60 d 离子分布如图10 所示。土柱表层的和 N a+的含量小幅增大，但并未发生表层积聚，无毛细水补给条件下，土柱中的水分运移速率逐渐降低直至接近停止，因此盐分的上升高度有限，并未全部迁移至表层。ⅡA 土柱上覆黄土层中各层位的含量均大于ⅡB，下伏含盐层的含量则小于ⅡB，Na+分布趋势与一致。ⅡA 土柱下伏含盐层60 d 的含量相比于初始值降低了64%，ⅡB 降低了45%。试验结果表明在无毛细水补给条件下，上细下粗型土柱下伏含盐层水盐向上的迁移总量和迁移速率均大于上粗下细型土柱。