斥水性土壤中水分运动规律的数值模拟

发布时间：2020-05-16 21:26

【摘要】：土壤斥水性(Soil water repellency,简称SWR)是一种广泛存在的土壤物理性质,对土壤水分运动产生较大影响。随着近年来数值模拟被广泛应用于土壤水分运动,研究斥水性土壤中水分运动的数值模拟对全面了解斥水性土壤中水分运动规律具有指导意义。本文结合不同情况下的试验数据,结合HYDRUS-1D模型,确定其初始和边界条件,进行了均质和层状斥水性土壤入渗、蒸发及遮雨棚下斥水土壤中夏玉米生长过程中根系吸水规律的数值模拟。模型模拟的效果评价采用相关系数(R~2)、相对均方根误差(RRMSE)和纳什效率系数(NSE)表达。对于室内入渗和蒸发试验,通过应用RETC软件对van-Genutchen土壤水分特征曲线进行土壤水力参数拟合,结合实测土壤水力参数残余含水率(θ_r)、饱和含水率(θ_s)和饱和导水率(K_s),对土壤水力参数α和n经过率定和验证,进而应用HYDRUS-1D模型对均质和层状斥水性土壤入渗、蒸发过程水分运动进行了模拟。遮雨棚下斥水土壤中夏玉米种植试验共5个处理(CK、WR1、WR2、WR3和WR4),斥水级别为亲水和轻微斥水,生育期内LAI与积温的关系通过高斯三参数方程进行拟合,根据LAI将Penman-Monteith计算的腾发量划分为蒸发量和蒸腾量。应用HYDRUS-1D根据实测的2016年含水率对土壤水力参数进行反算,得到K_s、α和n,结合实测的θ_r和θ_s,应用2017年的实测含水率数据进行参数的验证。最后对试验中较难直接观测的不同剖面的逐日含水率变化值以及根系吸水量进行了模拟。论文得出了以下主要结论:(1)应用HYDRUS-1D模拟亲水、轻微斥水和强度斥水级别的不同质地土壤水平吸渗和垂直入渗过程的水分运动过程模拟效果整体较好。对于亲水土壤,RRMSE、R~2和NSE的变化范围分别为2.2%-16.6%、0.93-0.999和0.783-0.998;对于轻微斥水土壤,RRMSE、R~2和NSE的变化范围分别为1.5%-11.7%、0.942-0.996和0.796-0.996;对于强斥水土壤,RRMSE、R~2和NSE的变化范围分别为8.4%-21.1%、0.876-0.99和0.774-0.989。对于同一土壤,随着斥水性的增加,模拟的效果有所降低,但是整体模拟效果较好,肯定了HYDRUS-1D在斥水性土壤水分运动模拟中的适用性。积水高度分别为4cm、6cm、8cm和10cm条件下,强度斥水土壤的预测模拟结果清晰地揭示了累积入渗量、湿润锋和剖面含水率随入渗时间的变化规律,为进一步研究不同质地、不同斥水性土壤中水分运动以及溶质运移等提供依据。(2)基于严谨的土壤水力参数率定和验证,应用HYDRUS-1D对层状斥水性土壤积水入渗过程中的土壤水分运动规律进行了数值模拟。层状土壤分别为亲水、轻微斥水、强度斥水和严重斥水X土以及亲水砂土,有X土夹砂土和砂土夹X土两个层状土壤试验处理。亲水土壤的夹层位置分别为5cm、10cm和20cm,斥水性土壤的夹层位置均为5cm。对于X夹砂情况,模拟的累积入渗量、湿润锋和剖面含水率的R~2变化范围分别为0.912-0.999、0.951-0.999和0.940-0.993;对于砂夹X情况,累积入渗量、湿润锋和剖面含水率的实测值与模拟值R~2变化范围分别为0.957-0.999、0.976-0.999和0.952-0.999。对于斥水性X土夹砂土以及砂土夹X土中夹层位置为10和20cm的情景,进行水分入渗过程的预测。不同斥水级别层状土壤水分运动过程的模拟结果显示,斥水性夹层土壤滞水造成水分回填从而降低水分入渗速率,很大程度影响水分运动过程,尤其是对X土夹砂土情景,当夹层位置为10cm时,土壤斥水性对水分运动过程的影响最为显著。对于砂土夹亲水X土处理,夹层位置为10cm对水分运动影响最大;对于砂土夹斥水(轻微、强度和严重斥水)X土,夹层位置为5cm对水分运动影响最大。土壤斥水性对水分运动的影响比夹层位置更为明显,尤其是砂土夹X土情景的试验处理,当夹层X土的斥水性由亲水增加到严重斥水,各个处理上层和夹层的水分变化不大,但是下层土壤的含水率显著降低。综上,HYDRUS-1D为研究斥水性层状土壤水分运动规律提供了有力的帮助。(3)对均质斥水性X土、砂土和盐碱土蒸发过程中的水分运动过程进行了模拟,斥水级别分别为亲水、轻微斥水和强度斥水。对于X土和盐碱土,累积蒸发量随着斥水级别的增加呈现明显减小的规律,砂土无明显规律。R~2、RRMSE和NSE值一致表明均质斥水性土壤蒸发率的模拟效果与累积蒸发量和含水率的模拟效果相比稍差。层状斥水性土壤蒸发试验处理与层状斥水性土壤入渗试验处理相同。对于亲水性层状土蒸发试验,累积蒸发量随着夹层位置的增加呈现显著上升的趋势。不同斥水级别的X土夹砂土的累积蒸发量随斥水级别的增加呈现下降的趋势。X土夹砂土的试验处理中,蒸发率的模拟效果较好。对于亲水性砂土夹X土,累积蒸发量随着夹层位置的增加而增加,但是在砂土夹斥水性X土试验中,累积蒸发量无明显规律,模拟效果较其他蒸发试验的稍差。通常HYDRUS-1D在斥水性土壤蒸发过程水分运动规律的模拟中具有一定的适用性。(4)遮雨棚下斥水性土壤中夏玉米生长情况下,在参数验证阶段,模拟5个试验处理含水率的R~2大于0.91,RRMSE小于19.4%,NSE大于0.95。对于2016年试验,模拟CK、WR1、WR2、WR3和WR4试验处理的蒸散量(ET_a)的R~2分别是0.82、0.86、0.87、0.74和0.73;对于2017年试验,CK、WR1、WR2、WR3和WR4试验处理的ET_a的R~2分别是0.81、0.74、0.71、0.72和0.66,土壤表层0.5cm处逐日含水率的模拟R~2大于0.65。因此,HYDRUS-1D的模拟效果整体较好。两年度5个试验处理的根系吸水模拟结果表明,逐日根系吸水在生长旺盛时期即播后50天左右达到峰值,尤其是2016年。根系吸水用来维持夏玉米生长,斥水性土壤中的根系吸水值要比亲水土壤中小。随着初始WDPT值的增加,根系吸水呈现明显下降的趋势。较弱的根系吸水对夏玉米的长势以及产量产生较大的影响,两年度4个斥水土壤处理的根系吸水较弱,土壤蒸发较强,夏玉米可吸收的土壤水分很少,因此斥水土壤中夏玉米的生长受到抑制。
【图文】：

技术路线图,斥水性,土壤水分入渗,夏玉米

蒸发运动的规律；4）遮雨棚控水条件下玉米生长过程中土壤水分运动过程的数值模拟。斥夏玉米试验分 2016 和 2017 两年度，试验土样为X土，斥水等级为亲水根据初始 WDPT 值分为 5 个试验处理。选用 2016 年试验数据进行夏玉系吸水模拟土壤水力参数的率定，，选用 2017 年试验数据进行土壤水力参出两年度夏玉米生长过程中根系吸水量，揭示斥水性土壤中夏玉米产量因。技术路线文旨在研究斥水性土壤中水分运动规律的数值模拟。首先应用 HYDRU平和垂直方向均质斥水性土壤中水分入渗运动进行模拟，以确定 HYDR性土壤的适用性；然后对层状斥水性土壤水分入渗过程、均质斥水性土层状斥水性土壤蒸发过程进行模拟，全面揭示室内斥水性土壤水分入渗律，进而进行遮雨棚下斥水性土壤中夏玉米种植试验根系吸水模拟，揭造成作物减产可能的原因。本论文的研究思路如图 1-1 所示。

曲线,水分特性,级别,盐碱土

土样（g kg-1）i（s）（g kg-1）i（s）（g kg-1）（s）水平方向X土 0.00 2.2 0.40 24 0.67 33砂土 0.00 0.9 0.40 7 0.67 11砂浆黑土 0.00 1.5 0.40 14 0.67 36盐碱土 0.00 3.2 0.40 41 0.67 85垂直方向X土 0.00 3.1 0.40 32 0.67 31砂土 0.00 0.5 0.40 9 0.67 13砂浆黑土 0.00 2.2 0.40 18 0.67 35盐碱土 0.00 4.0 0.40 27 0.67 10采用离心机测定亲水、轻微斥水和强度斥水的X土、砂土、砂浆黑土和盐碱土的水分特性曲线，测定结果见图 2-1。
【学位授予单位】：西北农林科技大学
【学位级别】：硕士
【学位授予年份】：2019
【分类号】：S152.7

【相似文献】