农田土壤水分渗深度模型.docVIP

农田土壤水分渗透深度模型 　　摘要：为了准确预测1次降水过程的土壤水分渗透深度，选择有代表性的站点，采用人工试验、土壤水分自动站观测相结合的方法，在分析河南省2010―2013年118个土壤水分自动观测资料和对应台站降水资料的基础上，进行回归分析，建立降水渗透深度模型。结果表明，各方程R值均大于0.80，F值均通过了P 中国论文网 /1/viewhtm　　关键词：降水量；土壤质量含水率；渗透深度；模型；麦田 　　中图分类号： S152.7+2 文献标志码： A 文章编号：1002-1302（2016）07-0419-03 　　土壤水分渗透是一个非常复杂的动态过程，不仅受到土壤性质、土壤初始含水量、降水量等因素的影响，而且与下垫面性质（植被类型、坡度、盖度等）、降水强度以及水温等因素有显著关系[1-2]。国内外已有大量关于土壤水分渗透规律的研究，取得了一些重要成果。国外常用的有Green-Ampt、Philip渗透模型，这些模型有明确的物理学意义，便于建立其特征参数与土壤物理特征间的关系，已经得到了广泛的应用[3-5]。原鹏飞等通过室内和野外模拟降水入渗试验，揭示了降水量对渗透深度、渗透速率的影响[6]。赵同应等分别利用不同的方法计算了降水渗透深度，对墒情、雨情评价有重要作用[7-9]。虽然前人对土壤渗透模型的研究较多，但是多偏重于机制方面，结构复杂，参数难于确定，给实际生产应用带来不便[10]。本研究探讨以大量的自动化观测资料和人工试验资料为依据，构建土壤水分渗透深度模型，以期为农田管理、灌溉决策和防旱抗旱措施的制定提供科学依据，为农业气象服务提供技术支撑，具有一定的科学意义和现实意义。 　　1 材料与方法 　　1.1 资料来源 　　收集了河南省118个土壤水分自动观测站2010―2013年土壤水分观测资料，各站点资料主要包括：土壤质量含水率、土壤水文物理常数（土壤容重、田间持水量和凋萎湿度）。土壤湿度采用河南省气象科学研究所和中国电子科技集团公司第27研究所共同研制的Gstar-1土壤水分自动监测仪进行数据测定，对自动土壤水分数据进行数据质量分析，剔除异常的土壤水分数据。降水量资料分别来自对应地区的气象台站。 　　根据河南省不同土壤类型的分布特点，分别选择有代表性的站点（郑州、鹤壁和黄泛区）进行人工试验观测。并利用郑州市农业气象试验站大型水分观测试验场进行人工模拟降水、作物截流量、墒情和渗透深度观测试验。人工降水机器装置采用“QYJY-501”式便携式野外模拟降水机，降水高度为4 m，降水强度为15 mm/h，模拟面积为4×8.25 m2，降水均匀系数在80%以上，模拟降水选择在无风天气条件下进行。 　　1.2 研究方法 　　农田水分平衡方程： 　　式中：Wt为雨后土壤含水量，mm；W0为初始土壤含水量，mm；P为过程降水量，mm；Ic为作物截留量，mm；ET为农田蒸散量，mm；r为地表径流量，mm；D为渗漏量，mm；G为地下水补给量，mm。 　　在一次降水过程中，由于河南省大部分地区为平原的特点，农田地势较平坦，一般情况下降水产生的径流较小，r=0；降水期间地表水汽含量接近饱和，同时降水历时一般较短，因此农田蒸散量ET=0；河南地区农田地下水位常年高于2 m，地下水对土壤水分的影响较小，G=0；渗漏量与降水量、地下水位土质、土壤水分状况有关[11-12]，本研究不予考虑；作物截留量与作物叶面积指数相关，根据前人研究结果，按照发育期进行不同级别划分[13-15]。方程（1）简化为： 　　以悬着水渗透模型为基础，选择对降水渗透深度影响较大的初始土壤湿度、过程降水量为分析因子[16-17]。应用2010―2013年各台站数据，首先根据降水过程前后自动站土壤含水量变化情况，计算各台站渗透深度值，统计各台站降水量和对应层次的平均初始土壤含水量，建立与渗透深度的相关关系；然后，对模型进行回代检验，并进行预报值、实测值的相关分析，利用绝对误差对模型精度进行统计分析。其中，数据统计分析使用SPSS19.0软件完成，最终通过C#和ArcGIS Engine8.3实现模型的可视化。 　　2 结果与分析 　　2.1 渗透深度的观测和计算 　　降水渗透深度的观测采用实测和自动站数据分析相结合的方法进行。2013年3月4日进行降水量为20 mm的模拟降水过程，并进行降水前、降水后墒情和土壤渗透深度观测。由图1可见，10、20 cm土壤水分在降水后有明显的增量，与人工观测土壤渗透深度为20 cm相同。同样，其他台站根据各层次土壤水分变化情况，可以推算出各自降水过程的土壤渗透深度值。 　　2.2 冬小麦不同生育期降水截留量 　　研究表明，作物冠层截留量除与作物本身特征因素（叶面积指数、植株地面生物量、株高及作物干燥度等）有关外，还