[发明专利]小麦自然水分亏缺干旱评估模型的GIS建模方法

说明书

技术领域

本发明涉及一种小麦自然水分亏缺干旱评估模型的GIS建模方法，属于农业灾害评估的技术领域。

背景技术

小麦是我国重要的粮食作物之一，小麦干旱灾害严重影响我国粮食安全，对小麦干旱的评估显得尤为重要。

自然水分亏缺率模型为作物自然供水量与需水量的差值占需水量的百分比，较好地反映了土壤、植物和气象三方面因素的综合影响，比较真实地反映出作物水分亏缺状况，是常用的作物干旱诊断方法之一。计算过程如下：

式中，E为小麦全生育期的需水量，W为小麦全生育期的自然供水量；单位均为mm。

E＝Kc×ET0(2)

式中，Kc为作物系数，小麦全生育期Kc取1.04；ET0为可能蒸(mm),采用FAO(1998)推荐的Penman-Monteith公式求得。

式中：ET0为参考作物蒸散量(mm/d)；Rn为净辐射[M J/(m2·d)]；G为土壤热通量[M J/(m2·d)]，当计算时段在10～30d时，日平均土壤热通量值很小，可忽略不计；而在逐日或长期估算时，G值通常很重要；U 2为2m高处风速(m·s-1)；Δ为饱和水气压温度关系曲线上在Ta处的切线斜率(kPa/℃)；γ为干湿表常数(kPa/℃)；Ta为日平均气温(℃)；Ea为饱和水气压(kPa)；Eb为实际水气压(kPa)。

小麦自然供水量(W)包括三部分：①小麦的土壤有效底墒(W1)；②小麦全生育期内的有效降水量(W2)；③地下水供水量(W3)。计算公式为：

W＝W1+W2+W3(4)

W1＝(Wt-Wd)×ρ×h×0.1(5)

式中，W1为小麦播种前的土壤有效底墒量，单位为mm；Wt为小麦播种前的实际土壤湿度，单位为％；Wd为凋萎湿度，取值6.5，单位为％；ρ为土壤容重1.44，单位为g/cm3；h为土层厚度，单位为cm；0.1为单位换算系数。

W2＝P(6)

其中P为实际降水量，单位为毫米(mm)。

在实际计算小麦供水量中，地下水供水量W3可忽略不计。

模型构造器(ModelBuilder)是ArcGIS提供的构造地理分析和处理工作流和脚本的图形化数据建模工具。ModelBuilder由输入数据、空间处理工具和输出数据等3个基本结构组成。

空间插值方法，即通过已获取的样点数据的规律，外推或内插到整个研究区域为面数据的方法。

栅格计算器是建立复杂的应用数学模型的基本模块，基于数学运算符以及数学函数将数值连续的空间栅格数据进行计算，得到所需要的空间栅格数据。

自然水分亏缺率模型涉及较多的数据和大量复杂的公式，在现实应用中存在易出错、效率低等弊端。本文将该模型与GIS相结合，逐步分块建立基于GIS的小麦自然水分亏缺干旱评估模型，从而改善原数学模型的缺点，提高模型的使用效率。

发明内容

本发明的目的在于克服现有技术存在的不足，提供一种便于使用者理解，提高模型使用效率的小麦自然水分亏缺干旱评估模型的GIS建模方法。

本发明采用的技术方案为：一种小麦自然水分亏缺干旱评估模型的GIS建模方法，包括以下步骤：

步骤一：首先将自然水分亏缺率模型分成2个部分：小麦全生育期的需水量E和小麦全生育期的自然供水量W。再将小麦全生育期的需水量E公式中的小麦可能蒸散量ET0分为8个部分。

步骤二：选取研究区一段时间的气象站点数据与土壤墒情报表，在ModelBuilder建模平台上采用反距离加权插值工具将点数据处理为数值连续的空间栅格数据。

步骤三：利用栅格计算器对小麦可能蒸散量ET0中的变量进行叠加分析，分别建模，包括饱和水气压温度关系曲线上在T_a处的切线斜率Δ、湿表常数γ、土壤热通量G、2m高处风速U2、饱和水气压es、日平均气温Ta、实际水气压ea以及净辐射R_n的建模。

步骤四：根据步骤三中已经建立的模型，同样利用ModelBuilder中的栅格计算器对小麦可能蒸散量ET0进行建模，并利用小麦可能蒸散量ET0模型对小麦全生育期需水量E进行建模。

步骤五：采用ModelBuilder中的栅格计算器分别对小麦的土壤有效底墒W1、小麦全生育期内的有效降水量W2及地下水供水量W3进行建模，同时对麦全生育期的自然供水量W进行建模。