CN110580657B - 农业灌溉需水量预测方法 - Google Patents

Abstract

本发明公开了一种农业灌溉需水量预测方法，其包括选择基准期，选取播种面积大于设定比例的作物作为典型作物；根据气象信息计算典型作物在生育期内的灌溉用水定额；获取节水区域内各气象站点的旬降雨量，并基于旬降雨量计算典型作物的有效降雨量；计算典型作物单位面积的灌溉需水量；判断是否存在至少一种典型作物单位面积的灌溉需水量大于连续实测年单位面积的实际灌溉用水量的平均值，若是，对所有典型作物单位面积的灌溉需水量进行修正；否则采用每种典型作物单位面积的灌溉需水量和下一年度种植面积，计算下一年度的总需水量。本方案通过实际灌溉用水量对计算得到的灌溉需水量进行修正，可以提高灌溉需水量预测的准确性和科学评估节水效果。

Description

农业灌溉需水量预测方法

技术领域

本发明涉及农业灌溉技术，具体涉及一种农业灌溉需水量预测方法。

背景技术

用水总量控制和定额管理是全面落实“节水优先”方针政策的重要抓手。农业作为用水大户，由于其用水过程不仅与农业种植结构有关，而且与气候和节水技术措施的推广及节水灌溉制度的实施情况等众多因素密切相关，复杂的作用因素使得目前难以合理确定有关农业用水量，通常采用历史种植结构和相应条件的灌溉定额水平的整体平均核定，缺乏必要的校验，难以真实的体现出农业节水措施带来的效用，从而导致目前绝大多数地区的农业用水基准偏高，增加农业用水成本，难以实现农业灌溉节水，水资源不能合理高效利用，使得种植时成本投入升高，农业产值降低。

发明内容

针对现有技术中的上述不足，本发明提供的农业灌溉需水量预测方法能够准确地预测下一年度总需水量。

为了达到上述发明目的，本发明提供的技术方案为：

提供一种农业灌溉需水量预测方法，其包括：

S1、获取节水区域的农业统计资料，根据农业统计资料选取播种面积大于设定比例的作物作为典型作物；

S2、获取节水区域内气象站在基准期内每旬的气象信息，并根据气象信息计算基准期内典型作物在生育期内的灌溉用水定额；

S3获取节水区域内各气象站点的旬降雨量并基于旬降雨量计算基准期典型作物的有效降雨量；

S4、根据基准期及同种典型作物的灌溉用水定额和有效降雨量，计算典型作物单位面积的灌溉需水量；

S5、判断是否存在至少一种典型作物单位面积的灌溉需水量大于连续实测年单位面积的实际灌溉用水量的平均值，若是，进入步骤S6，否则进入步骤S7；

S6、根据典型作物在基准期内所有的单位面积灌溉需水量和连续实测年的单位面积的实际灌溉用水量，对所有典型作物单位面积的灌溉需水量进行修正；

S7、采用每种典型作物单位面积的灌溉需水量和下一年度种植面积，计算下一年度的总需水量。

本发明的有益效果为：本方案通过基准期内的气象信息计算的灌溉用水定额和有效降雨量确定作物单位面积的灌溉需水量，并在灌溉需水量小于实际灌溉用水量时，基于实际灌溉用水量对灌溉需水量进行修正，以使得预测的灌溉需水量无限接近真实值，以提高预测的灌溉需水量的准确性，进而保证得到的总需水量预测的准确性。通过对灌溉需水量的逐年修正实现逼近真实值的过程，可以减少灌溉用水量，提高农业灌溉用水效益。

节水区域的管理者通过准确性较高的总需水量，能够在满足农业部规划用水的情况下，尽量保证较大的种植面积，从而保证农作物最终的经济收益。

附图说明

图1为农业灌溉需水量预测方法的流程图。

具体实施方式

下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。

参考图1，图1示出了农业灌溉需水量预测方法的流程图，如图1所示，该方法S包括步骤S1至步骤S7。

在步骤S1中，获取节水区域的农业统计资料，根据农业统计资料选取播种面积大于设定比例的作物作为典型作物；

在本方案中，将节水区域农作物总播种面积之比超过10％的农作物确定为典型作物，其对应的播种面积记为A_i,(i＝1,2,3，......n)，其中i为节水区域典型农作物种类，n为典型作物总数量。

在步骤S2中，获取节水区域内气象站在基准期内每旬的气象信息，并根据气象信息计算基准期内典型作物在生育期内的灌溉用水定额；

在本发明的一个实施例中，步骤2进一步包括：

S21、根据气象信息，采用彭曼公式计算典型作物在基准期的逐旬潜在蒸散量：

其中，ET₀为典型作物基准期时段内逐旬潜在蒸散量；Δ为饱和水汽压-温度曲线斜率；R_n为典型作物冠层的净辐射；G为增热土壤消耗的能量；γ为湿度计常数；T为平均气温；U₂为2m高风速；e_a为饱和水汽压；e_d为实际观测水汽压；

S22、根据同种典型作物的逐旬作物潜在蒸散量及生育期内的逐旬作物系数，计算典型作物的旬需水量：

ET_i,j,k＝ET₀×Kc_i,k

其中，ET_i,j,k为典型作物i在第j年生育期的第k旬的旬需水量，单位为mm；Kc_i,k是典型作物i在生育期的第k旬的逐旬作物系数；

S23、累加生育期内所有旬需水量作为典型作物的灌溉用水定额ET_i，j：

在步骤S3中，获取节水区域内各气象站点的旬降雨量，并基于旬降雨量计算典型作物的有效降雨量；

实施时，本方案优选基于旬降雨量和旬需水量计算典型作物的有效降雨量进一步包括：

S31、根据旬降雨量，采用空间插值方法计算典型作物在基准期的旬面雨量：

当节水区域的地形小于设定高程时，采用泰森多边形法计算典型作物在基准期的旬面雨量；

当节水区域的地形大于等于设定高程时，采用距离加权反比方法计算典型作物在基准期的旬面雨量。

S32、根据同种典型作物同一旬的旬面雨量和旬需水量得到旬降雨量PE_i,j,k：

PE_i,j,k＝min(P_i,j,k,ET_i,j,k)

其中，P_i,k,j为典型作物i在第j年生育期的第k旬的旬面雨量，单位为mm；

S33、累加生育期内所有旬降雨量作为典型作物的有效降雨量PE_i,j：

在步骤S4中，根据基准期及同种典型作物的灌溉用水定额和有效降雨量，计算典型作物单位面积的灌溉需水量；典型作物i单位面积的灌溉需水量q_i的计算方法包括：

S41、计算典型作物i在第j年的单位面积的灌溉需水量IR_i,j：

IR_i,j＝ET_i,j-PE_i,j，(i＝1,2,3，......n,j＝1,2,3，...m)；

S42、根据灌溉需水量IR_i,j，计算典型作物i单位面积的灌溉需水量q_i：

在步骤S5中，判断是否存在至少一种典型作物单位面积的灌溉需水量大于连续实测年单位面积的实际灌溉用水量的平均值，若是，进入步骤S6，否则进入步骤S7；

在步骤S6中，根据典型作物在基准期内所有的单位面积灌溉需水量和连续实测年的单位面积的实际灌溉用水量，对所有典型作物单位面积的灌溉需水量进行修正；

其中，对典型作物单位面积的灌溉需水量进行修正的计算公式为：

其中，q_i为典型作物i单位面积的灌溉需水量，单位为m³；q_ij为典型作物i在第j年的单位面积灌溉需水量，单位为m³；m为基准期；q_i,g为典型作物i在第g年的实际灌溉用水量，该数据可以通过各级政府部门发表的用水情况及对应典型作物种植面积，直接计算而得，单位：m³；n为典型作物的总数量；k'为连续实测年。

其中基准期的大小为5年，其为下一年度的上5年，比如下一年度为2020，那么基准期指的是2015、2016、2017、2018、2019。连续实测年可以选取邻近预测年的前面连续几年数据。

在步骤S7中，采用每种典型作物单位面积的灌溉需水量和下一年度种植面积，计算下一年度的总需水量；步骤S7进一步包括：

S71、采用每种典型作物单位面积的灌溉需水量和下一年度种植面积，计算每种典型作物下一年度的理论需水量：

Q_i＝q_i×A_i,(i＝1,2,3，......n)

其中，Q_i为典型作物i下一年度的理论需水量；A_i为典型作物i下一年度的种植面积；

S72、累加所有典型作物的理论需水量得到下一年总需水量Q：

在本发明的一个实施例中，农业灌溉需水量预测方法还包括修正下一年度的总需水量：

获取下一年度的规划用水量(其可以通过各级政府部门发表的规划进行确定)，并比较规划用水量和总需水量；

当总需水量大于下一年度的规划用水量时，根据调整策略调整典型作物的种植面积，并返回步骤S7；以及

当总需水量小于等于下一年度的规划用水量时，输出下一年度的总需水量和所有典型作物的当前种植面积。

其中，调整策略为：选取经济价值最小的典型作物，并按设定阈值缩小种植面积；或者选取灌溉需水量最大的典型作物，并按设定阈值缩小种植面积；或者将所有典型作物的种植面积按预设比例缩小；设定阈值大于预设比例。

下面结合具体的实例对本方案提供的方法进行说明：

本实施例以河北省沧州节水示范试验区为例：以2010-2014年为基准期，2015年作为下一年度，对其灌溉需水量进行预测。

典型作物确定

以《河北省统计年鉴》、《石家庄统计年鉴》和《沧州市社会经济统计信息》以依据，获得基准期(2010-2014年)沧州节水灌溉试验区农业种植结构和相应作物的播种面积，见表1，表中其他包括高粱和豆类；计算每种作物播种面积占区域作物总播种面积占比，超过10％确定区域典型农作物，确定典型作物主要为冬小麦、玉米。

表1新华区作物种植面积统计

气象信息的获取

通过国家气象站网和沧州地区气象监测数据，为避免雨量站稀疏影响气象信息面值确定的准确性，选择沧州气象站，选择距离新华区位置相对近的气象站——沧州、泊头、沧县三个气象站，选择基准期气象站逐日典型气象要素信息，主要气象要素包括降水、气温、风速、湿度、水汽压等日数据，汇总获得旬气象信息。

基准期内典型作物在生育期内的灌溉用水定额确定

采用PM公式，利用所获取的气象信息，计算典型作物基准期逐旬潜在蒸散量ET0，新华区基准期逐旬潜在蒸散量参见表2；然后乘以典型作物逐旬作物系数Kc，得到基准期典型农作物逐旬用水量。

表2新华区潜在蒸散量

基准年	2010	2011	2012	2013	2014
						潜在蒸散量(mm)	1021.33	1015.73	1071.50	956.13	898.97
1月上旬	8.04	6.77	6.67	7.30	8.33
						1月中旬	8.12	7.31	6.92	7.03	8.38
1月下旬	10.25	9.70	8.09	11.24	9.46
						2月上旬	10.80	10.90	9.84	8.88	11.32
2月中旬	11.91	13.00	11.20	11.36	12.80
						2月下旬	12.94	11.94	11.13	11.13	10.32
3月上旬	21.24	17.57	22.59	16.49	19.16
						3月中旬	26.92	21.26	23.76	23.87	22.62
3月下旬	33.26	32.15	31.89	28.59	32.59
						4月上旬	35.38	34.22	35.87	35.52	25.85
4月中旬	41.83	37.64	39.69	29.22	34.10
						4月下旬	40.96	39.42	44.44	35.92	40.56
5月上旬	40.30	45.44	47.19	44.83	37.94
						5月中旬	50.20	46.58	53.22	44.70	49.04
5月下旬	57.24	54.96	54.59	58.10	54.23
						6月上旬	52.78	51.63	50.34	48.20	51.26
6月中旬	53.48	58.92	56.13	43.53	49.94
						6月下旬	49.12	45.61	55.48	48.59	41.15
7月上旬	48.11	46.14	50.29	38.75	38.58
						7月中旬	44.23	47.78	40.54	47.11	36.99
7月下旬	48.67	44.20	48.93	44.11	36.82
						8月上旬	39.61	36.23	39.47	35.82	32.25
8月中旬	36.06	35.84	44.88	37.94	30.87
						8月下旬	40.51	41.77	45.38	38.44	32.95
9月上旬	31.39	32.10	36.17	32.89	27.58
						9月中旬	26.30	33.28	32.52	24.80	25.62
9月下旬	28.76	27.13	28.42	22.56	23.26
						10月上旬	23.65	26.21	28.53	21.86	19.47
10月中旬	20.71	22.13	21.32	22.27	16.12
						10月下旬	17.12	20.67	20.26	16.62	14.71
11月上旬	12.03	13.43	15.63	13.68	10.65
						11月中旬	10.47	11.28	13.46	11.79	8.76
11月下旬	8.69	10.85	10.77	7.69	7.38
						12月上旬	6.81	8.41	9.38	7.71	6.29
12月中旬	6.56	6.18	7.53	8.38	5.76
						12月下旬	6.85	7.10	8.97	9.21	5.90

其中，作物系数的确定结合河北省和沧州市农作物生长发育试验数据获得，新华区冬小麦的生育期从9月下旬到翌年6月中旬，夏玉米从6月中下旬到9月中下旬，冬小麦和夏玉米整个生育期作物系数见表3；累加生育期内所有旬用水量作为典型作物冬小麦和夏玉米的灌溉用水定额，具体见表4。

表3基准期主要作物系数(Kc)表

表4基准期逐年典型作物灌溉需水定额和平均灌溉需水量

典型作物生育期有效降水量确定。

选址的区域内以气象站点的旬降水量为基础，本次用泰森多边形法计算旬面雨量，结合典型作物旬需水量，计算典型作物生育期内的有效降雨量。典型作物灌溉基准期用水定额确定为：冬小麦230mm，夏玉米85.6mm。

采用本方案的步骤S4可以计算得到冬小麦在基准期单位面积平均灌溉用水量229.9mm(153.3m³/亩)；夏玉米在基准期单位面积平均灌溉用水量85.6mm(57.1m³/亩)。

实际中，根据2010-2014年政府部门下达的用水计量，一方面要压减冬小麦种植面积，同时要求区域内实现节水规划项目，即实施结构调整和冬小麦夏玉米水肥一体化节水技术措施。通过以上政策的实施，在2010-2014的连续实测年得到的试验区域的单位面积的实际灌溉用水量的平均值夏小麦仅为156.4m³/亩，夏玉米仅为55.4m³/亩。

比较基准期典型作物单位面积上灌溉需水量和实测典型作物单位面积灌溉后，发现夏玉米的实测单位面积灌溉需水量小于基准期理论值，需要进行修正。基于此采用步骤S6方式对夏小麦和夏玉米单位面积的灌溉需水量进行修正，修正后的典型作物冬小麦单位面积需水量为138.9m³/亩，较单位面积理论用水量偏少9.3％；夏玉米单位面积需水量为53.4m³/亩，较单位面积理论用水量偏少6.5％。修正前后各计算指标的对比见表5。

表5实测值修正前后典型作物基准期单位面积需水量变化

综上所述，通过本方案提供的预测方法得到的单位面积的灌溉需水量比较接近实际测量值，相对而言能够比较准确地预测下一年单位面积的灌溉需水量，从而保证了下一年度总需水量预测的准确性。

Claims (8)

1.农业灌溉需水量预测方法，其特征在于，包括：

S1、获取节水区域的农业统计资料，根据农业统计资料选取播种面积大于设定比例的作物作为典型作物；

S2、获取节水区域内气象站在基准期内每旬的气象信息，并根据气象信息计算基准期内典型作物在生育期内的灌溉用水定额；

S3、获取节水区域内各气象站点的旬降雨量，并基于旬降雨量计算基准期典型作物的有效降雨量；

S4、根据基准期及同种典型作物的灌溉用水定额和有效降雨量，计算典型作物单位面积的灌溉需水量；

S5、判断是否存在至少一种典型作物单位面积的灌溉需水量大于连续实测年单位面积的实际灌溉用水量的平均值，若是，进入步骤S6，否则进入步骤S7；

S6、根据典型作物在基准期内所有的单位面积灌溉需水量和连续实测年的单位面积的实际灌溉用水量，对所有典型作物单位面积的灌溉需水量进行修正：

其中，q_i为典型作物i单位面积的灌溉需水量；q_ij为典型作物i在第j年的单位面积灌溉需水量；m为基准期，以年计；q_i,g为典型作物i在第g年的实际灌溉用水量；n为典型作物的总数量；k'为连续实测年；

S7、采用每种典型作物单位面积的灌溉需水量和下一年度种植面积，计算下一年度的总需水量。

2.根据权利要求1所述的农业灌溉需水量预测方法，其特征在于，还包括修正下一年度的总需水量：

获取下一年度的规划用水量，并比较规划用水量和总需水量；

当总需水量大于下一年度的规划用水量时，根据调整策略调整典型作物的种植面积，并返回步骤S7；以及

当总需水量小于等于下一年度的规划用水量时，输出下一年度的总需水量和所有典型作物的当前种植面积。

3.根据权利要求2所述的农业灌溉需水量预测方法，其特征在于，所述调整策略为：选取经济价值最小的典型作物，并按设定阈值缩小种植面积；或者选取灌溉需水量最大的典型作物，并按设定阈值缩小种植面积；或者将所有典型作物的种植面积按预设比例缩小；所述设定阈值大于预设比例。

4.根据权利要求1所述的农业灌溉需水量预测方法，其特征在于，所述步骤2进一步包括：

S21、根据气象信息，采用彭曼公式计算典型作物在基准期的逐旬潜在蒸散量：

其中，ET₀为典型作物基准期时段内逐旬潜在蒸散量；Δ为饱和水汽压-温度曲线斜率；R_n为典型作物冠层的净辐射；G为增热土壤消耗的能量；γ为湿度计常数；T为平均气温；U₂为2m高风速；e_a为饱和水汽压；e_d为实际观测水汽压；

S22、根据同种典型作物的逐旬作物蒸散量及生育期内的逐旬作物系数，计算典型作物的旬需水量：

ET_i,j,k＝ET₀×Kc_i,k

其中，ET_i,j,k为典型作物i在第j年生育期的第k旬的旬需水量；Kc_i,k是典型作物i在生育期的第k旬的逐旬作物系数；

S23、累加生育期内所有旬需水量作为典型作物的灌溉用水定额ET_i，j：

其中，N为基准期作物生长期的旬总数。

5.根据权利要求4所述的农业灌溉需水量预测方法，其特征在于，所述基于旬降雨量计算基准期典型作物的有效降雨量进一步包括：

S31、根据旬降雨量，采用空间插值方法计算典型作物在基准期的旬面雨量；

S32、根据同种典型作物同一旬的旬面雨量和旬需水量得到旬降雨量PE_i,j,k：

PE_i,j,k＝min(P_i,j,k,ET_i,j,k)

其中，P_i,k,j为典型作物i在第j年生育期的第k旬的旬面雨量；

S33、累加生育期内所有旬降雨量作为典型作物的有效降雨量PE_i,j：

6.根据权利要求5所述的农业灌溉需水量预测方法，其特征在于，所述典型作物i单位面积的灌溉需水量q_i的计算方法包括：

S41、计算典型作物i在第j年的单位面积的灌溉需水量IR_i,j：

IR_i,j＝ET_i,j-PE_i,j，(i＝1,2,3，......n,j＝1,2,3，...m)；

S42、根据灌溉需水量IR_i,j，计算典型作物i单位面积的灌溉需水量q_i：

7.根据权利要求5所述的农业灌溉需水量预测方法，其特征在于，采用空间插值方法计算典型作物在基准期的旬面雨量进一步包括：

当节水区域的地形高程和采用的气象站点数与计算区域内的距离小于设定值时，采用泰森多边形法计算典型作物在基准期的旬面雨量；

当节水区域的地形高程和采用的气象站点数与计算区域内的距离大于等于设定值时，采用距离加权反比方法计算典型作物在基准期的旬面雨量。

8.根据权利要求1所述的农业灌溉需水量预测方法，其特征在于，所述步骤S7进一步包括：

S71、采用每种典型作物单位面积的灌溉需水量和下一年度种植面积，计算每种典型作物下一年度的理论需水量：

Q_i＝q_i×A_i,(i＝1,2,3，......n)

其中，Q_i为典型作物i下一年度的理论需水量；A_i为典型作物i下一年度的种植面积；

S72、累加所有典型作物的理论需水量得到下一年总需水量Q：

Images