matlab模拟三体运动_[2020.Vol 126]实验室条件下模拟渠灌土壤水分运动的新入渗模型...

A new infiltration model for simulating soil water movement in canal irrigation under laboratory conditions

实验室条件下模拟渠灌土壤水分运动的新入渗模型

1

介绍

合理设计灌溉系统可以提高灌溉用水效率。分析灌溉过程中不同时期的入渗趋势、累积入渗速率和湿润锋的迁移距离，可以促进灌溉系统的合理设计。通过田间试验分析土壤水分入渗趋势需要大量的人力物力，并且田间工作受季节条件等因素的制约。通过室内入渗试验可以相对简单有效地获得入渗趋势，并可以直接观察入渗过程。因此，这类实验受到了广泛的关注。为了确定影响渠道二维入渗过程的主要因素，建立了一个准确有效的土壤入渗模型，动态模拟了农田灌溉渠道的入渗过程。在室内试验的基础上，建立了由供水单元和长方体有机玻璃单元组成的渗滤系统。建立了五个不同的初始条件，分别为土壤容重，初始含水量、河道底部宽度、入渗头和坡度系数。在此基础上，研究了人工渠道断面的累积入渗和湿润锋过程，建立了试验设置条件下的二维土壤水分运动方程。采用HYDRUS-2D软件对模型进行求解，并对模拟值与实测值进行了误差分析。

2

材料方法

(1)室内试验

设备包括长方体有机玻璃容器和供水系统(马里奥特烧瓶)(Fig 1)。透明有机玻璃容器主要用于储存土壤样品和观察入渗湿润锋的迁移，供水系统则为入渗过程提供不同水头压力情景下的水分。

试验土样于2017年5月在哈尔滨市东北农业大学水利与土木工程学院实验田采集。在10-30cm的深度采集土壤样本。根据美国农业部(USDA)制定的土壤分类标准，本试验采用的土壤类型为砂壤土。室内实验于2017年春季在东北农业大学水利与土木工程学院进行。实验步骤见Table 1。

(2)数值模拟

设计试验条件下的渠道入渗过程可简化为二维非饱和土壤水分运动问题，用Richards方程描述：

用van Genuchten-Mualem模型描述了Richards方程中非饱和导水率、土壤基质势与含水量的关系：

渗水量计算公式如下：

渗透系数计算公式如下：