水利工程论文华北地区夏玉米田间水分转化规律的研究.docVIP

水利工程论文-华北地区夏玉米田间水分转化规律研究摘要：本文利用一维垂直非饱和土壤水运动计算模型，通过模拟计算夏玉米生育期各生长阶段不同水文年、不同灌水定额下田间水分转化规律，对棵间蒸发，深层渗漏与灌水定额之间的数量关系进行了定量研究，并建议华北地区夏玉米节水灌溉定额为40m/亩. 关键词：夏玉米 田间水分转化 灌水定额 ?1? 田间水分转化的动力学模式　　田间水分的转化过程主要包括降雨或灌水的入渗，径流，根区土壤水分运移，根系吸水，蒸发蒸腾，根区以下土壤水运动等.描述土壤水运动的动力学方程是根据达西定律和连续方程结合推导出来的，由于方程是非线性的，所以只有在特定的初始和边界条件下才能用解析法求解.自60年代以来，随着电子计算机技术的发展，人们能够借助计算机运用数值模拟的方法，对一般条件下的土壤水分运动问题求解，并得到了满意的结果.本文采用HYDRUS模型 软件 进行数值模拟.　　1.1? HYDRUS模型? HYDRUS模型 软件 ［8］是美国盐碱实验室在Worm模型的基础上的改进版，用于模拟计算一维垂直非饱和流和溶质运移，考虑了作物根系吸水和土壤持水能力的滞后影响，适应于恒定或非恒定的边界条件，具有灵活的输入输出功能，可用来模拟非匀质土壤，最多可模拟5种20层土壤.模型中方程解法采用Gakerin线性有限元法.　　1.2? 方程和边界条件? 在忽略土壤侧向水流运动，仅考虑一维垂向运移时，有根系吸水项的土壤水分运动方程为［1］：c(h)(h)/(t)=()/(Z)k(h)(h)/(Z)-k(h)-s(z,t)　　(1)式中c(h)——比水容重［L-1］,c(h)=dθ/dh,h——土壤压力水头［L］，θ——体积含水量［L3L-3］,k(h)——水力传导度［LT-1］，s(z,t)——单位体积根系吸水率［T-1］,z——土壤深度［L］、向下为正，t——时间［T］.　　初始条件：h(z,t)=h0(z), t=0.　　上边界条件：h(0,t)=h0(t),??? (积水条件)[-k(h)(h)/(z)+k(h)]/z=0=q0(t)　　(2)式中h0(t)——已知压力水头［L］，q0(t)——净通量，正通量表示下渗，负通量表示蒸发(棵间蒸发).　　下边界条件：h(l,t)=hi(t),[-k(h)(h)/(z)+k(h)]/z=l=q1(t)或(h)/(z)/zl=0,? (3)2? 作物根系吸水模式　　本文采用VanGenuchten模型［8］，即s(z,t)=Ep(t)ξ′,(z)σ(h,h0)　　(4)式中S(z,t)——单位根据吸水率，Ep(t)——作物最大蒸腾率，σ(h,h0)——盐分应力函数，反映土壤盐分对田间根系吸水的影响.σ(h,h0)=(1)/(1+(h+h0)/(h50))p　　　(5)式中h——压力水头，h0——渗透压，与溶液浓度C有关，h0=a1c,a1简单换算系数，其值依赖于压力水头和浓度表达的单位.h50——作物潜在蒸腾率减少50%的土水势，对玉米而言，h50为-0.25—-0.65MPa，一般取-0.43MPa(Ehler,1983),P——经验常数，P≈3；ξ′(z)——根系密度分布函数相对值.? (6)其中 Lr——根层深度，ξ(z)——根系密度分布函数，ξ(z)是根深Lr的函数，实际应用中将Lr分为若干层且认为每层内根系分布是均匀的.根系密度采用层内干根重占根区总干根重的比值表示，由田间实测获得.3? 非饱和导水率的确定　　非饱和导水率K是土壤水分运动的重要参数，在对非饱和土壤水运动基本方程进行数学分析时，无论用解析解或数值解的方法，都要用到它.获得K的方法有两种：一是实验方法，如瞬时剖面法，垂直下渗通量法，垂直土壤稳定蒸发法，结壳法等；另一种是间接法，即根据已知的水分特征曲线和其他条件推导出的函数形式，如Van Genunchten(1980)将土壤水分特征曲线的函数形式与Mualem(1976)导出的用来预测非饱和导水率的函数形式相结合得到如下函数关系［5］：　　　　　　　　　　　　θ(h)=θr+(θs-θr)/(［1+(αh)n］m),(7)　　　　　　　　　　　　k(l)=kss1/21[1-(1-s1/ml)m]2,(8)　　　　　　　　　　　　　　　sl=(θ-θr)/(θs-θr),(9)式中θs——饱和含水量；θr——残余含水量；Ks——饱和导水率；Sl——相对饱和度，α,n,m——拟合参数，通过非线性最小二乘法对室内试验获得的土壤持水数据(θ,h)进行拟合求模糊.将式(7)代入式(8)得：k(h)=ks［1+(αh)n］-2/m｛1-［1-(1+αh)n］-m｝2　　(10)表1 VG非饱和导水率函数有关参数土层深度/cm饱和