蒸发条件下宁夏银北地区盐渍土水分运移数值模拟.docVIP

蒸发条件下宁夏银北地区盐渍土水分运移数值模拟 　　摘要：对盐渍土中的水盐动态进行定量分析是提出适宜于宁夏银北地区盐渍土调控措施的基础工作。蒸发条件下土壤水分运动模拟初步研究表明，将有限体积法应用到求解蒸发条件下的土壤水分运动模拟计算中，模拟值与实测值吻合较好；春季土壤含水率模拟曲线在纵断面较为一致，地表蒸发量对0.2～0.6 m的土壤含水率有较大的影响。 　　关键词：土壤水分；蒸发；数值模拟；盐渍土 　　中图分类号：S152.7 文献标识码：A 文章编号：0439-8114（2015）15-3638-04 　　DOI：10.14088/j.cnki.issn0439-8114.2015.15.013 　　Abstract： It’s the basis of the work to use the numerical simulation method for finding appropriate control measures in Yinbei region of Ningxia province. The results showed that there was a good agreement between simulated and measured values when the finite volume method was selected to simulate the movement of soil moisture evaporation. The simulated curve of soil moisture in the horizontal section was consistent with that of the real values. Surface evaporation had a greater impact on the soil moisture of 0.2～0.6 m. 　　Key words： soil moisture； evaporation； numerical simulation； saline soil 　　宁夏回族自治区属于黄河中上游半干旱-半漠境盐渍区，盐渍化土地面积为20 718 km2，占自治区总面积的40%。其中尤以银北地区分布最广，盐渍化面积为17 070 km2，约占全区盐渍化土地总面积的81%左右，该地区含盐量0.6%～1.0%的重度盐渍土就占到本区盐渍化土壤面积25.66%[1]。由于银北地区地形低洼，排水不畅，地下水位高致使水中的盐分在强烈的蒸发作用下堆积于土壤内，这是形成土壤盐渍化的主要原因。通过对该地区的盐化土地实施水利工程改造和农业生物改良调控措施，盐渍化得到了整体控制[2]。但近几年来该地区局部的盐渍化现象明显加重，主要是由于打浅井使用高矿化度地下水进行灌溉及淋洗，加剧了土壤的盐害[3]。这就需要提出更为适宜银北地区土壤盐渍化调控措施，作为该地区农业可持续发展和生态环境保护的基础。 　　本研究是通过建立能够描述宁夏银北地区盐渍土中水盐运移的数学模型，采用数值计算方法进行模拟计算，对盐渍土中的水盐动态进行定量分析。分析模拟结果，掌握水分、盐分在盐渍土体中的基本运动规律，为确定合理的土壤盐渍化调控模式提供科学依据。 　　1 材料与方法 　　1.1 土壤水分运动数学模型 　　土壤水分运移模型是基于对流扩散方程和初边值条件共同构建的模型，模型中的各类土壤水动力参数均需要通过试验来确定。该模型属于非线性偏微分方程，当前最有效的求解方法是数值计算方法。有限差体积法[4]是当前求解流动问题数值计算方法中最成功的方法，该方法是将方程中的每一项在所划分的控制体积进行积分，各项积分均有明确的物理意义，构成线性方程组进行求解计算。 　　1.2.1 有限体积法离散过程 采用有限体积方法对二维非饱和土壤水分运动方程进行线性离散。 　　其中，？准0为初值含水率（cm3/cm3），D（？准）为非饱和土壤扩散率（cm2/min），为非饱和土壤导水率（cm/min），？着为水量交换强度（cm/d）。 　　对方程在图1所示的控制体内积分可得： 　　由于全隐式格式可以在较大时间步长保持计算结果的稳定性[5]，所以采用全隐格式对扩散项进行时间上的积分，设？驻x=？驻z=？驻d则？驻V=A?？驻x=A?？驻z，得： 　　1.2.2 边界条件处理 上边界条件为地表，与当地的气象因素相关联，即与垂直向的水量交换强度？着（cm/d）有关。当？着为负数时为该时期的蒸发强度。根据监测蒸发量数据即可得到整个监测区域内的水量交换强度值。 　　采用上述相同方法给出边界条件在计算区域内离散得到各个端点的离散系数，见表1。以上系数通过编写M