1. 引言
干旱或半干旱地区的河流一般为季节性河流，汛期河道洪水往往具有暴涨暴落的特点，行洪期内的河道渗漏补给是本区域地下水资源的主要来源，也是当地人民生活和生产的主要水源。确定河道渗漏量及渗漏规律，不仅对于确定地下水补给量，而且对于季节性河道的洪水演进计算也具有重要意义[1] -[4] 。
土壤入渗主要包括垂向入渗和侧向入渗。聂卫波等在非饱和土壤水分运动理论基础上，利用Richards模型模拟不同试验处理下的一维垂向入渗过程，并用SWMS-1D软件对模型求解，通过试验得出，土壤质地、容重对一维垂向入渗累计入渗量、入渗率和湿润锋运移距离影响较大，土壤初始含水率和入渗水头对其影响较小[5] 。马东豪等利用一维垂向入渗试验确定了Brooks-Corey模型参数，并对其方法进行检验。结果表明，只要入渗历时足够长，利用预测的参数，数值解和提出的近似解均能很好的模拟土壤水分入渗过程[6] 。范严伟等通过对不同初始含水率条件下多种典型的土壤垂直一维入渗特性进行模拟，发现土壤累积入渗量变化过程符合Philip模型，提出了包括初始含水率、饱和含水率和饱和导水率的Philip模型参数的确定方法，通过验证表明，所建模型可反映土壤水分垂直一维入渗特点[7] 。
本文采用季节性河床砂样，进行了垂向土柱非饱和入渗试验，研究季节性河道行洪前期影响垂向入渗的主要因素，为后续开展现场入渗试验、建立季节性河道入渗模型以及研究季节性河道渗漏对洪水演进影响提供科学依据。
2. 材料与方法
2.1. 供试样品
试验砂样采自青岛大沽河干流河床，选取河道表面砂土以下10 cm的砂样。采样点为有代表性的大沽河上游、中游和下游，上游取样点位于莱西市的孙受镇，中游取样点位于即墨市的李家庄村，下游取样点位于胶州市的李哥庄村。分别命名为粗砂、中砂和细砂。砂样的颗粒级配(质量百分数)和水理性质分别参见表1和表2。
2.2.试验装置
试验用土柱底部直径为15 cm，柱高100 cm，材质为有机玻璃，土柱底部密布小孔，可以防止砂土流失，土柱下部加有由有机玻璃做的存水室，土柱中的水可以从小孔流到下面，存水柱侧壁设有出水口；土柱侧壁上每隔10 cm埋设1个张力计，用以测定张力变化情况，共设有8根张力计。土柱的供水装置为马氏瓶，材质为有机玻璃，直径25 cm，高90 cm。马氏瓶顶部中间插有一根有机玻璃管，为马氏瓶提供定水头。试验装置见图1。
2.3.试验方法
试验时间为2015年06月~08月，试验过程中实验室温度变化范围为21.7℃~30.8℃，湿度为25%~82%。试验前，要首先去掉砂样中的杂质，然后加入一定量水并加以搅拌，使含水量均匀。加水时注意加水量，水量太多则张力计水银柱无法显示读数。装砂样前，在土柱底部铺设一层纱网以防止砂粒下漏，堵塞底部小孔。然后将砂样装入试验装置中。每装一定的高度时，要均匀夯实，夯实砂样时要注意力度，太大则有可能压坏埋在砂里的陶土头，太小则不能保证密实度。保证砂样与土柱的边壁接触良好，水流不会从边壁渗漏。为了保证两层
Table 1. Grain-size distribution of sand samples
表1. 砂样的颗粒级配(%)
Table 2. Hydraulic properties of sand samples
表2. 砂样的水理性质
Figure 1. Device of vertical column, 1-Markov bottle, 2-Soil column, 3-Tensiometer, 4-Outlet, 5-Water storage chamber
图1. 垂向土柱装置，1-马氏瓶，2-土柱，3-张力计，4-出水口，5-存水室
砂样间的均匀性，在下一层砂样装填之前，要把砂土表面划毛。用一定量的砂样来测定初始含水率。砂样装完后，用薄膜将土柱口封住，防止砂样水分蒸发，影响实验结果。在马氏瓶侧壁加上直尺，记下初始读数，通过水量下降可以求得入渗量。
待到张力计水银柱读数上升到最大值时，实验开始。试验前，记下张力计水银柱初始读数和马氏瓶中液面高度。先在土柱砂样上方铺设塑料膜以防止水分下渗，在塑料膜上加水至与马氏瓶中有机玻璃管底部相同的高度，此时记为0时，打开马氏瓶开关并迅速抽出塑料膜开始入渗试验。实验过程中要保证马氏瓶气密性良好。然后每隔单位时间记录下马氏瓶侧壁直尺读数和张力计读数，得到实验过程中入渗量变化和土柱内水分变化情况。
3. 试验结果与分析
对三组砂样分别做两组垂向入渗实验，控制不同初始含水率，其中，细砂的两组初始含水率分别为1.02%和5.11%，中砂初始含水率分别为1.38%和3.51%，粗砂初始含水率分别为1.83%和5.92%。
(1) 垂向水分分布
试验得到垂向入渗中张力计水银柱读数随时间的变化过程，根据试验结果，以垂向土柱的深度为纵轴，土柱表面为纵轴0点，向下方向为正；根据不同砂样的水分特征曲线和张力计读数得到不同深度在一定时间的体积含水率，以体积含水率为横轴，绘制不同砂样垂向入渗时，在不同入渗时间和不同深度砂样的体积含水率。本文分别绘制了细砂、中砂和粗砂在较小含水率时入渗4 min时垂向不同深度的水分分布，分别绘制了细砂、中砂和粗砂在较大含水率时入渗4 min和入渗8 min时垂向不同深度的水分分布。如图2所示。
由图2可知，同种砂样在不同初始含水率时，入渗速度不同，不同砂样在初始含水率相近时，入渗速度也不同。对于同种砂样，初始含水率不同，垂向入渗速度相差比较大，初始含水率越小，入渗速度越快，初始含水率越大，入渗速度越慢。细砂样在初始含水率等于1.02%时，土柱内距离进水口最近的砂样只用4.5 min即可渗完，而当细砂初始含水率为5.11%时，同样位置的砂样则需要6 min。同样对于中砂，在初始含水率等于1.38%的时候，土柱只用5 min即可入渗完成，而当初始含水率为3.51%时，同样位置的砂样则需要7 min。对于粗砂，当初始含水率为1.83%时，入渗在5 min即可完成，而当初始含水率为5.92%时，同样深度的砂样则需要6 min。由土壤入渗理论解释如下。
垂向入渗前期入渗速度的快慢，也即湿润峰运移的速度可通过入渗水流湿润区的土水势梯度的概念进行解释。非饱和土壤水的运动是由土水势差引起的，水势梯度越大，运动越快，反之亦然。土壤初始含水率越高，水分入渗锋面的土水势越高，则水分入渗锋面与表层之间的平均势梯度越小，土壤的入渗速度也就越小。因此，同一砂样初始含水率小的入渗速度较快。粗砂、细砂入渗速度不同主要是因为砂质结构不同，细砂相比于粗砂，干容重大，结构更密实，大空隙数量少，小空隙多，比表面积大，相同含水率下，细砂的土水势高(基质势大)，水分入渗锋面与表层之间的平均势梯度低，因此细砂的入渗速度就慢些。
(2) 入渗率变化曲线
马氏瓶中液面高度随时间逐渐降低，根据液面降低的高度可以求得入渗量，入渗量除以单位时间内土柱的断面面积即为随时间变化的垂向入渗率。根据试验结果，得到入渗率随时间的变化过程。以入渗率的大小为纵轴，入渗时间为横轴，绘制不同砂样在不同含水率时入渗率变化曲线，结果见图3。
由图3可知，对于不同的砂样，入渗率的变化趋势大致相同，入渗开始时，由于含水率较小，入渗速度快，入渗率大，随着试验的进行，入渗速度逐渐降低，入渗率逐渐减小，最后趋于一个稳定值。对于同一种砂样，当初始含水率不同的时候，砂样入渗率的稳定值也不相同，若初始含水率较小，则入渗率的稳定值较大，且从入渗开始到入渗率减小到稳定值的时间较长。当初始含水率分别为1.02%和5.11%时，细砂的稳定入渗率约为0.027 mL/(s·cm2)和0.021 mL/(s·cm2)；当初始含水率分别为1.38%和3.51%时，中砂的稳定入渗率约为0.028 mL/(s·cm2)
(a) (b) 
(c) (d)
(e) (f)
Figure 2. Curves for vertical moisture distribution of sand
图2. 砂样垂向水分分布曲线
(a) (b) 
(c) (d)
(e) (f)
Figure 3. Curves for vertical-infiltration rate of sand varying with time
图3. 砂样垂向入渗率随时间变化曲线
和0.026 mL/(s·cm2)；当初始含水率分别为1.83%和5.92%时，粗砂的稳定入渗率约为0.030 mL/(s·cm2)和0.029 mL/(s·cm2)。对于不同的砂样，当初始含水率相近时，粗砂的稳定入渗率值比细砂要大，这是由于粗砂颗粒间的空隙大，相对于细砂更容易入渗。
4. 结论
垂向土柱非饱和入渗试验结果表明：
(1) 不同的砂样，干容重越大，入渗率越小。细砂的最大干容重最大，入渗率最小，粗砂的最大干容重最小，入渗率最大；
(2) 同种砂样的不同初始含水率，入渗率也不相同。当初始含水率小时，入渗率大，当初始含水率大时，入渗率小；
(3) 垂向入渗初始阶段，由于含水率小，所以入渗率大。随着时间的进行，入渗率逐渐减小，当砂样逐渐饱和，入渗率逐渐趋于稳定值。
(4) 当初始含水率分别为1.02%和5.11%时，细砂的稳定入渗率约为0.027 mL/(s·cm2)和0.021 mL/(s·cm2)；当初始含水率分别为1.38%和3.51%时，中砂的稳定入渗率约为0.028 mL/(s·cm2)和0.026 mL/(s·cm2)；当初始含水率分别为1.83%和5.92%时，粗砂的稳定入渗率约为0.030 mL/(s·cm2)和0.029 mL/(s·cm2)。粗砂的稳定入渗率大于细砂。
基金项目
国家自然科学基金资助项目(51309137, 51379097)。