查元源，刘洁颖，伍靖伟，刘 松

（1.武汉大学水资源与水电工程科学国家重点实验室，武汉 430072；2.长江水利委员会长江勘测规划设计研究院规划处，武汉 430010）

亭子口水利枢纽低地垫高农田防护工程保水性能研究

查元源1，刘洁颖2，伍靖伟1，刘 松1

（1.武汉大学水资源与水电工程科学国家重点实验室，武汉 430072；2.长江水利委员会长江勘测规划设计研究院规划处，武汉 430010）

针对亭子口水利枢纽农田防护工程低地垫高方案的保水性能，在室内大型土柱试验的基础上，建立基于Richards方程的一维饱和-非饱和水流运动数学模型，模拟分析不同防护方案条件下土壤水分的运动过程和保水性能。结果表明：目前拟采用的70 cm耕作层防护方案通过适度压实保水性能良好，所建模型和所得相关成果可以为类似工程提供参考。

亭子口水利枢纽；低地垫高；保水性能；数值模拟

为了保护亭子口水利枢纽库区耕地免受库区尾水淹没，解决移民生产生活问题，农田防护工程拟采用低地垫高方案，即先将耕作层剥离，然后用砂卵石混合土层垫高至安全高度，再将耕作层回填。低地垫高方案实施后，由于土层结构和土壤水运动条件的改变，可能会带来一些问题：下垫层的材料为沙砾石混合土材料，透水性较强，保水性能较差，可能形成俗称的“倒蒙金”土，会导致渗漏、土肥流失等问题。为了验证工程的可行性，同时取得工程施工时的必要参数，作者在武汉大学水资源与水电工程国家重点实验室进行了大型土柱保水性能试验，并在此基础上建立了一维饱和-非饱和水流运动数学模型，进而对不同防护方案条件下土壤水分运动过程和保水性能进行了数值模拟分析。

1 试验介绍

1.1 试验装置及土样材料

试验在大型土柱中进行。土柱由多层40 cm× 40 cm×30 cm的方形钢制容器拼接而成；每层土柱壁以10 cm为间距安装TDR土壤含水率探头，实时监测土壤含水率变化；最下一层土柱底部留有排水口，监测排水情况。

试验耕作层土样为在库区罗家桥和安全坝两地挖取的土样，分别代表保护区的砂土和粉砂土类型。在取土前进行了土壤原位测试，获取了耕作层土样的容重、孔隙度、田间持水量等基本参数。垫层材料来自库区河床附近的天然砂卵石料混合土。试验装置及填充土料如图1所示。

图1 大型土柱试验装置Fig.1 Large-scale soil column test device

1.2 试验方法和方案

试验方法主要是在土柱上部人工灌水，观测水层的深度变化；通过TDR观测不同深度土壤含水率随入渗时间的变化；在土柱底部收集排水，观测排水时间和排水量。通过上述3种数据定量分析该工程方案的保水性能。

各试验方案主要改变灌溉水量、土壤初始含水率、耕作层土壤质地、耕作层厚度、耕作层干密度，以对比不同试验方案下土柱保水性能的差异。

2 数学模型

2.1 控制方程

一维垂向土壤水流运动用Richards方程描述，即

式中：θ为土壤体积含水率；h为土壤水分负压（cm）；垂直坐标z向下为正（cm）；t为时间（h）；S为根系吸水项或其它源汇项（h-1）；K为土壤非饱和水力传导度（cm·h-1）。

2.2 定解条件

边界条件：模拟中，下边界排水，采用渗出面边界。其表达式为，当h＜0时，q＝0；否则h＝0。其中：h为边界处土壤负压；q为边界流量通量。上边界灌溉，采用大气边界，其表达式为，h＞0时，θ＝θs；h≤0时，q＝0。其中θ为边界处土壤含水率。本文中上边界均为初始时刻建立20 cm水层，自由入渗。

初始条件：h（z，0）＝h0（z），即初始时刻的土壤剖面负压已知。

2.3 土壤水分运动参数

2.3.1 土壤水分运动参数模型

土壤水分运动参数主要包括土壤水分特征曲线参数和非饱和土壤水力传导度，可应用van Genuchten［1］模型来表示：

式（2）至式（6）共有5个未知参数，分别为θr，θs，α，n，Ks。其中θr为残余体积含水量；θs为饱和体积含水量；Ks为饱和水力传导度（cm·h-1）；Se为饱和度。α（cm-1），n，m和l是反映土壤水分曲线特征的经验常数，本文中l取为0.5，m由式（4）求得。2.3.2 土壤水分运动参数选择

根据实测的土壤颗分曲线和天然干密度，利用RETC软件拟合，得到天然干密度下试验土料的土壤水分运动参数，如表1所示。

不同干密度下，土壤上述各参数值的变化将在方案分析中具体说明。

表1 天然干密度下的土壤水分运动参数Table1 Soil hydraulic properties in natural density

3 部分试验数据及模型验证

多组试验取得了各时间含水率剖面、灌水水层厚度、累积排水量3类数据。数据显示本工程采用低地垫高方案是合理的。限于篇幅，本文对试验数据不作详细说明，仅列出单组试验部分数据，以确定模型和参数的选择是否合适，数值模拟是否合理。

该单组试验耕作层为70 cm罗家桥砂土，下垫层为30 cm砂砾石混合土，土壤水分运动参数见表1，初始含水量θ0＝0.17，并且是按天然干密度均匀填装的。用HYDRUS-1D软件模拟和试验中TDR监测分别得到的耕作层土壤剖面体积含水率结果对比如图2所示。

图2 土柱耕作层不同深度含水率试验值与模拟值比较Fig.2 Com parisons ofmoisture content at different depths from experiment and simulation at arable layer in soil column

从图2可以发现，试验得出的土壤入渗速度与模拟值是一致的，而且各个时间段的剖面含水率曲线也符合很好。但图中也反映出存在一些差异，产生这些差异的一个可能原因是土柱试验的土层是分层手工填装的，土壤干密度难以达到完全均匀，加上土壤自身的自然沉降作用，导致下部土壤干密度有所增高，因此模拟结果有细微差别，但整体模拟结果与实测一致。

图3是耕作层下边界累计排水量的数值模拟和实测结果，两者非常吻合。这进一步说明，所建立的模型和所使用的参数能够模拟农田防护工程低地垫高方案的土壤水分运动过程。

图3 土柱下边界累计排水量试验值与模拟值比较Fig.3 Com pare of cumulative drainage volumes from experiment and simu lated data at lower boundary in soil column

从图3还可以看出，灌水后2 d共排水不到4 cm，即计算土层（包含30 cm的下垫层）能蓄住80%的灌溉水量。而且试验和模拟结果均表明，充分灌溉2 d后，罗家桥砂土的平均含水率为0.22～0.29；安全坝粉砂土为0.30～0.38，均超过天然田间持水量，能满足一次灌溉维持长时间作物需水要求。

4 方案分析

4.1 不同耕作层土壤质地方案

耕作层土壤质地不同，其保水性能可能不同。为此，再取安全坝的粉砂土作为耕作层土壤，耕作层仍为70 cm，下垫层仍为30 cm，进行模拟计算。两种情况下设置的观测点（15 cm深和45 cm深处）体积含水率变化对比如图4所示。

图4 不同耕作层土壤质地下的土壤含水率变化Fig.4 Variety ofmoisture content of arable layers with different soils

模拟结果表明，粉砂土的平均含水率高于砂土。这说明，采取不同的耕作土壤，保水性能会有较大差异，工程中应尽可能选择保水性能较好的土壤。

4.2 不同耕作层土壤深度方案

根据渗流垂直于土壤分层面时渗透系数主要决定于渗透系数值较小的土壤［2］，可知低地垫高方案中起主要保水作用的是渗透系数相对小的耕作层土壤。将耕作层调整为100 cm、下垫层仍为30 cm，模拟结果与耕作层70 cm、下垫层30 cm的方案对比如图5所示。

图5 不同耕作层厚度下的土壤含水率变化Fig.5 Variety ofmoisture content of arab le layers w ith different thicknesses

从图中可以看出，耕作层厚度分别为70 cm和100 cm方案的入渗速度是同步的。而由于罗家桥砂土的n值较大，退水较快，导致同一时间同一深度下是100 cm耕作层厚度的土壤含水率较小。这说明，在耕作层厚度达到70 cm后，继续增加耕作层土壤厚度不能明显增强保水作用。考虑到植物根系生长一般需求耕作层大于等于70 cm，工程施工时可以取土壤耕作层厚度为70 cm。

4.3 不同土壤干密度方案

不同压实程度，亦即不同干密度，必然带来土壤的渗水性能和保水性能发生变化。通过不同干密度下土壤吸力试验，得知当某种土壤的干密度ρd变大时，参数θs，Ks会减小；而θr，α，n的变化幅度非常小。因此这里假设密度改变时，同种土壤的参数θr，α，n保持不变。

理论上，土壤饱和时，水完全充满土壤的孔隙，那么土壤饱和含水率与土壤干密度存在以下理论关系，即

式中：p为土壤孔隙度；ρd为土壤干密度（g·cm-3）；ρs为土粒密度（g·cm-3）。

但实际中土壤水分往往难以完全充满孔隙，相关试验也发现干密度改变时，土壤饱和含水率只能近似按式（7）变化，有时需要对其进行修正［3］。本次试验通过测量罗家桥砂土不同压实程度时的干密度和相应饱和含水量，得到了二者在一定范围内的相关关系为

对于罗家桥砂土，ρs＝2.69。由式（8）知，当罗家桥砂土ρd分别为1.37 g／cm3（0.9倍原干密度）和1.66 g／cm3（1.1倍原干密度）时，对应的θs分别为0.481和0.375。

对于饱和渗透系数与土壤干密度的关系，本文中采用刘杰公式［4］，

式中：Ks单位为cm／s；d20为等效粒径；对于罗家桥砂土取为0.008 1cm。

根据土壤密度和孔隙度的关系，可将式（9）改写为

当罗家桥砂土的干密度变为1.37 g／cm3和1.66 g／cm3时，代入式（10）计算得饱和渗透系数从4.37 cm／h，分别变为6.70 cm／h和3.09 cm／h。该数据得到了试验的验证。

干密度采取4种不同方案：①田间干密度，即1.51 g／cm3；②田间干密度的1.1倍；③田间干密度的0.9倍；④耕作层上层40 cm采用0.9倍田间干密度，下层30 cm采用1.1倍田间干密度。模拟结果对比如图6所示。

图6 不同耕作层土壤干密度方案下的土壤含水率剖面Fig.6 M oisture content of arable layersw ith different soil densities

从图6可以看出，采用干密度为天然干密度的0.9～1.1倍时，土壤灌水48 h后的含水量（亦即田间持水量）均大于天然田间持水量0.23。

此外，还可以看出方案4为最佳方案。方案4的初始入渗速度与方案3同步，但退水曲线比方案3滞后，因此，在灌水后同一时间上方案4在土壤40 cm深度以内保持的水分更多。实际上，这样做就使得土壤结构有了“蒙金土”的特征。采取“上松下密”的原则，在不显著增加入渗速度的情况下，提高了计划湿润层土壤的田间持水量，减少了深层渗漏量。同样可知增大垫层干密度对保水持水都非常有利。

通过试验也得到了相同的结论。同时试验还表明：当土壤过于松散，在灌水等条件下，土壤会在重力作用下自然沉降、压实；而使土壤非常密实既增加成本也不利于作物生长。因此，干密度变动应有一个范围，可根据实际情况参照方案4进行施工。

5 结 论

本文根据水库淹没区低地垫高农田防护工程的实际需要，在室内进行大型土柱试验，同时建立低地垫高双层土壤结构的非稳定土壤水分运动数值模型，对不同防护方案进行了模拟分析，得到如下结论：

（1）通过试验检验了所建立的低地垫高土壤水分运动模型的合理性，模型可适用于其他低地垫高农田防护工程的保水性能分析；

（2）室内试验和数值模拟结果表明，亭子口水利枢纽库区农田防护工程采用低地垫高方案保水性能良好；

（3）耕作层土壤厚度取70 cm即可。施工时，耕作层采用上松下密的原则，并宜在天然干密度的0.9～1.1倍变动比较合适；垫层尽量压实，保证大于其天然干密度。

［1］ Van GENEUCHTEN，M.Th.A closed-form equation for predicting the hydraulic conductivity of unsaturated soils［J］.Soil Science Society of America Journal.1980，44：892-898.

［2］ 薛禹群.地下水动力学［M］.北京：地质出版社.1997.（XUE Yu-qun.Groundwater dynamics［M］.Beijing：Geological Publishing House.1997.（in Chinese））.

［3］ 吕殿青，邵明安，刘春平.容重对土壤饱和水分运动参数的影响［J］.水土保持学报，2006，20（3）：154-157.（LU Dian-qing，SHAOMing-an，LIUChun-ping.Effect of bulk density on soil saturated watermovement parameters［J］.Journal of Soil and Water Conservation，2006，20（3）：154-157.（in Chinese））

［4］ 沙金煊.农田不稳定排水理论与计算［M］.北京：中国水利水电出版社，2004：35-38.（SHA Jin-xuan.Theory and Calculation of Farmland Unsteady Drainage［M］.Beijing：China Waterpower Press，2004：35-38.（in Chinese） ）

（编辑：周晓雁）

Study on W ater Retention Performance of Low land Heightening Farm land Protection Project in Tingzikou Water Control Pivot

ZHA Yuan-yuan1，LIU Jie-ying2，WU Jing-wei1，LIU Song1
（1.State Key Laboratory ofWater Resources and Hydropower Engineering Sciences，Wuhan University，Wuhan 430072，China；2.Changjiang Institute of Survey Planning Design and Research，Wuhan 430010，China）

To research thewater retention performance of lowland heightening farmland protection project in Tingzikou Water Control Pivot，large-scale soil column experiments and one-dimensional saturated-unsaturated flow simulations have been performed.The results show：（1）70 cm thick arable soil layer with appropriate compaction can ensure satisfying water retention performance；（2）the established model and the results obtained in the paper can provide a reference for similar projects.

Tingzikou Water Control Pivot；lowland heightening；water retention performance；numerical simulation

TV93；S277.7

A

1001-5485（2010）04-0016-04

2009-05-07

查元源（1987-），男，安徽怀宁人，硕士研究生，主要从事地下水资源与环境方面的研究，（电话）13100616800（电子信箱）zhayuan87＠126.com。