基于天气预报信息的参考作物需水量预报研究

白依文, 鲁梦格, 程浩楠, 严冬, 孙怀卫

PDF(2820 KB)
中国农村水利水电 ›› 2021 ›› (10) : 175-179,184.
农田水利

基于天气预报信息的参考作物需水量预报研究

作者信息 +

Research on the Reference Crop Evapotranspiration Forecast Based on Weather Forecast Information

Author information +
稿件信息 +

摘要

获得作物需水信息是农业灌溉管理与水资源优化配置的基础。为探索利用公共天气预报信息进行作物需水预测,研究尝试以PM公式和HS公式为基础,将天气预报信息进行定量转化后预测参考作物需水量,以期解决典型地区参考作物需水量预报难题。通过以武汉市为典型研究区,收集了武汉站2020年3月16日至2020年4月15日的每日对未来七天的天气预报数据,使用天气预报数据与PM公式和HS公式得出的预报值与使用历史监测数据和PM公式采用实测数据得出的标准值进行对比并对本文使用的方法进行了验证。结果表明,不同预报期内,预报值与标准值有一致的变化趋势,两种方法在1~2 d预报期的MAE平均值分别为0.69和0.675 mm/d;RMSE平均值分别为0.91和0.905 mm/d,基本满足需水量精度要求,而随着预报期的增加准确性在降低;两种方法的对比显示,HS法的相关系数总是相对较高且预报后期比PM法稳定,PM主要误差是由实际日照时数导致。当实际日照时数预报精度较差时推荐使用HS法,此方法具有良好的物理基础并且数据容易获取。

Abstract

Obtaining crop water requirement information is the basis of agricultural irrigation management and optimal allocation of water resources. In order to explore the use of public weather forecast information to predict reference crop evapotranspiration, this paper is an attempt to use PM model and HS model as the basis for quantitatively transforming weather forecast information to predict reference crop evapotranspiration so as to solve the problem of water demand forecast in typical areas. By using Wuhan as the typical study area, the daily forecast data of Wuhan Station from March 16, 2020, to April 15, 2020, daily for seven days in the future are collected. The results show that the predicted value and the standard value have the same trend in different forecast period, MAE of the two methods in 1~2 days forecast periods are 0.69 and 0.76 mm/d, RMSE are 0.91 and 0.905 mm/d respectively, which basically meet the precision requirements of water demand, while with the increase of forecast period, the accuracy is decreasing. The comparison between the two methods shows that the correlation coefficient of HS is always higher than that of PM and the prediction during later period is more stable. The main PM error is caused by the actual sunshine hours. HS method is recommended when the prediction accuracy of actual sunshine hours is poor. This method has good physical basis and data is easy to obtain.

关键词

参考作物腾发量 / 公共天气预报 / Penman-Monteith法 / Hargreaves-Samani法

Key words

reference crop evapotranspiration / public weather forecast / Hargreaves-Samani / Penman-Monteith

基金

国家重点研发计划项目(2018YFC1508301)
国家自然科学基金项目(51879110)
湖北省水利重点科研项目(HBSLKY201818)

引用本文

导出引用
白依文 , 鲁梦格 , 程浩楠 , 严冬 , 孙怀卫. 基于天气预报信息的参考作物需水量预报研究[J].中国农村水利水电, 2021(10): 175-179,184
Yi-wen BAI , Meng-ge LU , Hao-nan CHENG , Dong YAN , Huan-wei SUN. Research on the Reference Crop Evapotranspiration Forecast Based on Weather Forecast Information[J].China Rural Water and Hydropower, 2021(10): 175-179,184
当前,我国农业用水占据了总用水量的60%左右,这其中的大部分是灌溉用水。获得作物需水信息,是农业灌溉管理与水资源优化配置的基础。作物需水量可以通过参考作物腾发量(ET 0)获得1ET 0的准确预报对提高灌溉系统管理水平、农业水资源利用率以及缓解水资源短缺具有十分重要的意义。
随着数值天气预报精度的提高,基于气象数据和物理过程的间接法受到广泛关注,并认为具有更好的准确性。Duce等在美国加利福尼亚和意大利的站点,运用数值天气预报获得的气象数据,评估并比较了几种常用的计算方法,得到的研究结果证明利用数值天气预报能够很好地进行预测2。因此,许多学者使用天气预报数据进行了预测研究34,以便为当地的灌溉计划制定者提供可行的操作方案。如Cai等将日常免费的天气预报信息转换成相应的计算数据,从而利用Penman-Monteith(PM)方法构建对北京等站点的实时灌溉预报模型5;罗玉峰等6、晏成明等7将气温数据输入Hargreaves-Samani(HS)公式中分别对西藏林芝站、广东青年运河灌区湛江站进行预报。
目前已提出的ET 0间接预测方法或模型超过50种,其中综合法的PM公式因其较高的准确性和良好的适用性被广泛使用89。温度法的HS公式10因为使用的气象数据较少、可操作性强,应用也很广泛。也有一部分研究人员利用掌握的大量历史数据构建GARCH类模型11、回归型支持向量机(SVR)预报模型12、BP人工神经网络模型13-15或多元线性回归模型等,都有较高的准确度但计算相对复杂。为此,探讨适合湖北省武汉的参考作物腾发量预报方法,验证典型地区的ET 0预测模型准确性以解决当地选择作物模型难题。
本文旨在运用中国气象局地面资料中武汉站的气象预报数据以及历史实测数据,将PM公式和HS公式计算得出的预报值与运用实测数据与PM公式得出的计算值进行比较,从而评价使用公共天气预报进行预测的效果,为当地农业经营者提供可操作性强、准确性高的灌溉调度计划。

1 数据和方法

1.1 数 据

武汉站位于中国的中部地区,地理坐标位置为东经114.03°,北纬30.36°,站点海拔高程为23.6 m。从中国天气网(http://www.weather.com.cn/)收集了武汉站从2020年3月16日至2020年4月15日的逐日对未来7天的气象预报数据。从中国气象数据网(http://data.cma.cn/)收集了武汉站2020年3月16日至2020年4月21日的逐日历史实测数据。天气预报数据主要包括天气类型、最高气温、最低气温、平均气温、风力等级;历史实测数据主要包括最高气温、最低气温、平均气温、平均风速、相对湿度、实际日照时间。

1.2 PM预测ET 0的预报方法

采用联合国粮农组织(FAO)推荐的计算参考作物腾发量的Penman-Monteith方法1
ET0,PM=0.408ΔRn-G+γ900Tmean+273U2es-eaΔ+γ1+0.34U2
式中:ET 0,PM是PM公式计算得到的参考作物腾发量,mm/d;R n为作物表面净辐射,MJ/( m2·d);T mean代表平均气温,单位为℃;U 2为2 m高处的平均日风速,m/s;γ是湿度常数,kPa/℃;e se a分别代表饱和蒸气压和实际水汽压,kPa;Δ是饱和水汽压—气温曲线的斜率,kPa/℃。
利用FAO-56 Penman-Monteith(PM)公式计算ET 0,不仅需要当日的气温、平均风速,还需要太阳辐射以及实际水汽压等参数。关于天气预报中变量转换相关内容见表1
表1 天气预报转换参数

Tab.1 Parameters transformed form weather forecast

PM变量 天气预报 转换变量
R n 天气类型 实际日照时数
U 2 风力等级 10m高度风速
e a 最低温度 假设为露点温度
Cai等5提出了一种量化每日天气预报数据来间接计算PM公式所需气象参数的方法。先将不同的天气类型转化为日照时数(表2所示),其中最大可能日照时数由日期与站点位置信息获取。再将日照时数转化为短波辐射Rs,太阳净辐射(R n)根据净短波辐射减去净长波辐射得到,即完成了天气类型估算太阳辐射。
表2 天气类型量化指标表

Tab.2 Quantitative indicators of weather type

天气类型 晴天 晴转多云 多云 阴天 雨天
n/N 0.9 0.7 0.5 0.3 0.1
注: n为实际日照小时数,h;N为最大可能日照小时数,h。
公共天气预报中的风力等级同样是定性预报。根据中国气象总局发布的地面气象观测标准16,风力等级对应于距地面10 m高处的风速如表3所示,而后根据Allen提出的风速转化公式,可将z m处的风速对应到2 m处的风速1
  u2=uz4.87ln (67.8z-5.42)
式中:u 2为离地面2 m高处的风速,m/s;uz 离地面z m高处的风速,m/s。
表3 风力等级量化指标表

Tab.3 Quantitative indicators of wind level

风力 0 1 2 3 4 5 6
u 10/(m·s-1 0.1 1.0 2.0 4.0 7.0 9.0 12.0
u 2/(m·s-1 0.1 0.7 1.5 3.0 5.2 6.7 9.0
注:武汉站预报风力等级整体较小,截取至6级风力;表中风速为平均值。
z取10 m时,得到的转化系数为0.747 95,从而得到风力等级与2 m高处的风速的对应关系。值得说明的是,尽管FAO对参考作物的定义为高度0.12 m,参考作物处的风速相对较低,而PM公式对风速相对不敏感。Allen et al.计算全球两千个气象站点的平均风速得到的均值为2 m/s,因此在风速参数缺失的情况下,可假定风速为2 m/s的定值进行计算1。事实上,杨洋对全国4个气候区总计61个研究站点进行风速预报的评价,对比了2 m/s定值、预报值、多年逐日均值和年均值,发现75%的站点中2 m/s定值的预测效果是最好的17
实际水汽压无法直接获得,通常是由相对湿度或露点温度转化而来。由于公共天气预报缺乏露点温度的预测,所以假设了日最低气温即为露点温度。事实上,夜间的气温下降到最低时,空气湿度接近饱和,因此最低气温是对露点温度的良好估计值1。实际水汽压的计算公式表示如下:
ea=e0(Tmin)=0.611exp 17.27TminTmin+237.3
式中:ea 为实际水汽压,kPa;t 0为饱和水滤质,kPa;T min为日最低气温,℃。

1.3 HS预测ET 0的预报方法

Hargreaves等18提出的当太阳辐射、相对湿度和风速数据难以获取时仅仅基于日最高气温和日最低气温计算的预报方法,是典型的基于温度的简单计算方法。其公式为:
ET0,HS=0.408CRaTmax-TminE Tmax+Tmin2+T
式中:ET 0,HS为预报值,mm/d;Ra 为太阳辐射(地球外辐射),可由日序数以及站点的地理纬度计算得到,MJ/(m2·d);T max为日最高气温,℃;T min为日最低气温,℃;CET为公式的3个参数,分别取0.002 3、0.5、17.8。

1.4 ET0预报精度统计指标

本研究中,用实测数据与FAO-PM公式计算出的结果作为标准值与运用预报数据计算出的预报值进行比较。为了更好地利用公共天气预报进行ET 0预测的精度检验,运用3个统计评价指标对标准值与预报值进行比较。统计指标分别为均方根误差(RMSE)、平均绝对误差(MAE)和相关系数(r)。

2 结果与讨论

2.1 公共天气预报的准确性分析

主要评价的气象因子包括最高气温、最低气温、风速。由于无法将预报过程中的最低温度和标准值计算过程中的相对湿度进行直接对比,因此选取它们共同的解析后的实际水汽压来进行对比,从而实现天气预报信息的评价。此外在ET 0的预测中,由于实际日照小时数是通过天气类型转换得到的,该转化方式以及天气类型预报的准确性应一起评价(表4)。
表4 武汉站参数预报统计评价指标结果

Tab.4 Statistical evaluation index of parameter forecast at Wuhan Station

预报期/d 1 2 3 4 5 6 7
最高温度 r 0.86 0.87 0.87 0.86 0.87 0.86 0.75
RMSE/℃ 2.74 2.66 2.78 2.71 2.78 3.05 3.51
MAE/℃ 1.54 1.91 2.20 2.06 2.19 2.43 3.02
最低温度 R 0.72 0.69 0.68 0.65 0.73 0.75 0.73
RMSE/℃ 2.64 2.73 3.02 3.22 2.86 2.77 2.81
MAE/ 2.02 2.14 2.37 2.40 2.16 2.21 2.15
风速 R 0.75 0.77 0.77 0.67 0.61 0.59 0.46
RMSE/(m·s-1 1.79 1.88 1.95 1.69 1.73 1.64 1.17
MAE/(m·s-1 1.45 1.53 1.57 1.38 1.38 1.38 0.87
日照实数 R 0.26 0.17 0.25 0.25 0.26 0.19 0.06
RMSE/h 4.25 4.26 4.10 4.31 4.23 4.77 4.64
MAE/h 3.48 3.51 3.35 3.58 3.44 3.93 3.40
实际水汽压 R 0.68 0.67 0.63 0.53 0.58 0.65 0.70
RMSE/kPa 0.27 0.27 0.28 0.34 0.30 0.27 0.26
MAE/kPa 0.19 0.19 0.20 0.22 0.20 0.18 0.18
最高空气温度(T max)和最低气温(T min)用来计算饱和水汽压、实际水汽压以及净长波辐射。表3显示,在研究期内从预测期第一天至第七天,总体上预报的最高气温和最低气温的MAERMSE都在增加,而r值都在减小,这说明了随着预报期的增加,武汉站各气温预报的准确度在下降。此外,将最高气温与最低气温之间的统计指标分别对比可以看出,预报期小于五天时,最高气温的预报精度要优于最低气温的预报精度。另外值得注意的是,最低气温的预报精度随预报期增加呈现先下降后上升的趋势。可能的原因是本研究的时间跨度小,且为春季,武汉属于亚热带季风气候,其日最高气温的变化在该季节波动不大,相对稳定,应当考虑到季节变化趋势对最高气温和最低气温的影响。总体来看,气温预报的精度并不算太高。
随着预报期的增加,武汉站风速预报统计指标r值在下降,风速预报的相关性在逐渐降低。MAERMSE,二者并没有明显的随预报期变化的趋势。直接对比预测值与历史数据可以发现,绝大多数的风速预报都是小于三级(取为1.5 m/s)或是三到四级(取4.1 m/s),而实际测得的日平均风速往往比预报的风力等级转化得到的风速还要低,可见该季节武汉风力较弱,从而使得MAERMSE的值较小也较为稳定。此外对比以往研究,如Yang等19将国内6个站点的风力预报与实测数据进行对比发现相关系数仅为0.06~0.27,而本次研究中武汉站风力预报的r值相对较高,其变化范围为0.46~0.77。本研究中较高的风速预报精度对提高ET 0预报的精度有积极的影响。
准确地计算短波辐射是利用PM公式进行高精度预报的关键,而天气类型的定量转化是整个短波辐射计算的重要环节。随预报期增加,日照时数的预报精度在下降,与气温预报出现相似的结果。且日照时数的r值显著低于其他各类气象因素的r值,因此日照实数可能是造成ET 0预报产生误差的主要原因,这个结论与Yang等提到的相同。
预报过程中,实际水汽压是在假定最低气温近似露点温度的情况下得到的[式(3)],而标准值计算过程中,选用了相对湿度而没有使用露点温度,二者不同的计算原理导致实际水汽压计算结果的差异。预报期从1~7 d,实际水汽压的相关系数r先减小后增大,而MAERMSE先增大后减小,总的变化趋势与最低气温随预报期的变化趋势相同。同时通过相关系数的对比可以看出,实际水汽压的相关系数与最低气温的相关系数近似。因此,可以认为最低气温预测的准确性对实际水汽压的准确性有一定影响,使用最低气温估算露点温度的方法是可靠的,同时Allen通过不同物理参数计算实际水汽压的公式是比较精确的,具有相当的等价性。
公共天气预报的气象数据被应用于PM和HS公式的计算,在探讨ET 0预报精度前,说明各气象数据的预报值与实测值间的误差,以及随着预测期不同,误差变化的趋势等尤为重要。从而可以分析在实际运用过程中,参数的气象预报对ET 0预报准确性的影响。

2.2 PM公式与HS公式计算值对比

为验证利用公共天气预报进行ET 0预报的准确性,用HS公式和PM公式进行同时期的预测,并将两种方法的预测结果与以实测数据通过FAO56-PM公式所得出的标准值进行比较。本研究由于观测数据量小时间跨度短,两种预报方法未经过率定,直接选用公式推荐值得出预报值。
预报期分别为1、4、7 d时,预报值和计算值的散点图如图1所示。在不同预报期,HS预报结果的相关系数要比PM预报高,尤其在预报期为7 d时,PM公式预测值的相关系数极低。但两个公式的预测值与标准值相比都有一定的偏差,HS与PM总是高估了实际的ET 0。造成这种现象可能的原因除了气象数据预报存在误差外,HS未考虑日照实数和湿度对腾发量的影响。
图1 预报期为1、4、7天的PM与HS预报值和标准值散点图

Fig.1 Scatter diagram of PM and HS forecast value and standard value with forecast period of 1, 4 and 7 days

Full size|PPT slide

预报期从一天到七天,两种方式预报结果的统计指标变化趋势如图2所示。总的来说,随预报期增加,两种方式的MAERMSE增加且r值下降,有相同的变化趋势。HS方法的MAE式中要比PM方法大,RMSE在预报期4天后超过了PM方法,而相关系数总是保持在高于PM的水平,总的来说HS方法的误差较为稳定,如果能够对武汉站的参数进行率定,应该能将误差消除一部分。二者的相关系数有高度相似的变化趋势,进一步印证了温度是预报ET 0比较重要的气象参数。
图2 预报期从1~7 d的PM与HS预报的统计指标

Fig.2 Statistical index of PM and HS forecast with forecast period from 1 to 7 days

Full size|PPT slide

2.3 PM公式与HS公式ET 0预报的计算准确性分析验证

不同预报期计算结果的折线图如图3所示,由于研究期处于春季,ET 0的值主要分布在2~4 mm/d。结果显示不同预报期的预报值与标准值之间有着一致的变化趋势。对比不同的预报期可以发现,预报期越短,预报值与标准值的波动越一致。反之,差异越大。
图3 预报期为1~7 d的预报值与实测值的折线图

Fig.3 Line chart of forecast value and measured value with forecast period of 1~7 days

Full size|PPT slide

随着预报期的增加,ET 0预报的r下降,RMSEMAE增大(如表5所示),说明预报的精度随预报期增加而下降。这与各气象因素预报值的准确性随预报期增加而下降有关。PM法与HS法7天预报值的MAE均值分别为0.82和0.90 mm/d,RMSE的均值分别为1.06和1.05 mm/d,误差还是不容忽视的。而两种方法在1~2 d预报期内的预报值与实测值吻合相对较好[如图3(a)~(b)所示],具有一定的准确性。PM法与HS法1~2 d预报期内预报值的MAE均值分别为0.69和0.76 mm/d;RMSE均值分别为0.91和0.905 mm/d;r均值分别为0.675和0.715。单独观察PM法的相关系数,当仅考虑1~6 d的预报期时,均值为0.53。但是预报期为7 d时,相关系数变化较大,仅为0.26,并使均值下降到0.49,将这一现象与天气预报评价和实际水汽压评价进行对比可以发现,预报期为7 d时的最低气温以及实际水汽压的相关系数与前几个预报期相差不大,而最高气温、实际日照小时数和风速变化较大,分别由0.86下降到0.75、0.19下降到0.06、0.59下降到0.46。说明此时预报的误差与这3个因素的误差有密切关系。此外,从5 d预报期到6 d预报期,预报的r实现了一定幅度的增长,对比同一时期天气预报统计评价指标可以看到,该阶段内只有最低气气温预报的相关系数存在增长现象,并带动了实际水汽压相关系数从0.58到0.65的增长。说明了最低气温的准确预报对提高预报的准确性有较大的影响。温度法HS在同一预报期r也实现了较大的增长,进而说明了温度预报准确性的重要。
表5 武汉站ET 0预报统计评价指标结果

Tab.5 Statistical evaluation index of ET 0 at Wuhan Station

预报期/d 1 2 3 4 5 6 7 M
PM r 0.70 0.65 0.45 0.45 0.43 0.49 0.26 0.49
RMSE/℃ 0.96 0.86 1.06 1.10 1.08 1.11 1.23 1.06
MAE/℃ 0.73 0.65 0.79 0.91 0.86 0.92 0.91 0.82
HS r 0.72 0.71 0.53 0.48 0.50 0.63 0.48 0.58
RMSE/℃ 0.95 0.86 0.96 1.10 1.09 1.11 1.25 1.05
MAE/℃ 0.79 0.73 0.80 0.95 0.97 0.98 1.10 0.90
注: M为平均值。

3 结 论

本文以武汉站为典型地区,实现了利用公共天气预报信息运用PM和HS公式进行短期ET 0预测。为了评价该方法,收集了历史实测数据得到ET 0的实测值,将其作为标准值与预报结果进行对比。得到的主要结论如下:
(1)ET 0预报的精度在预报期1~2 d尚可,PM与HS方法预报期从1~2 d的MAE均值分别为0.69 mm/d和0.76 mm/d;RMSE均值分别为0.91 mm/d和0.905 mm/d;r均值分别为0.675和0.715。
(2)随预报期增加,天气预报准确性下降,从而导致两种方法的ET 0预报的准确性下降。预报期第七天准确性下降程度较大,预报期在两天内的预报值较为理想。
(3)关于武汉的短期公共天气预报对气温预测的准确性尚可接受,所以PM其误差主要是由实际日照小时数引起的。
(4)当气温精度较高时,不需要收集及转换大量输入气象数据,可利用系数校正后的HS公式对武汉地区进行需水量预报;若日照实数预报精度较高时,可采取率定后的PM法进行预报。
本文所选用的短期ET 0预报方法与时间序列方法相比更具有物理基础;与数值天气预报(NWP)方法相比较而言其天气预报的数据更容易获取,因此其未来广泛使用的可能性更大。

参考文献

1
ALLEN R G PEREIRA L S RAES D, et al. Crop evapotranspiration Guidelines for Crop water requirements[R]. FAO Irrigation and Drainage Paper No.56, 1998.
2
DUCE P SNYDER R L SOONG S T, et al. Forecasting reference evapotranspiration. Third International Symposium on Irrigation of Horticultural Crops1999537(1):135-141.
3
蔡甲冰,刘钰,雷廷武,等. 根据天气预报估算参照腾发量[J]. 农业工程学报2005(11): 19-23.
4
张倩,段爱旺,王广帅,等. 基于天气预报的参照作物腾发量中短期预报模型研究[J]. 农业机械学报201546(5):107-114.
5
CAI J LIU Y LEI T, et al. Estimating reference evapotranspiration with the FAO Penman-Monteith equation using daily weather forecast message [J]. Agricultural and Forest Meteorology2007145:22-35
6
罗玉峰,李思,彭世彰,等. 基于气温预报和HS公式的参考作物腾发量预报[J]. 排灌机械工程学报201331(11):987-992.
7
晏成明,刘梦,罗玉峰 .基于公共天气预报的广东青年运河灌区参考作物腾发量预报方法比较[J].排灌机械工程学报201836( 8):38-743.
8
刘钰,蔡林根. 参照腾发量的新定义及计算方法对比[J]. 水利学报1997(6):28-34.
9
赵琪,罗玉峰, 彭世彰, 等. 基于天气预报和Penman-Monteith公式的短期逐日参考作物腾发量预报[J]. 节水灌溉2014(1):1-4,8.
10
刘 梦,罗玉峰,汪文超,等. 基于天气预报的漳河灌区参考作物腾发量预报方法比较[J]. 农业工程学报201733(19):156-162.
11
孙怀卫,严冬,陈皓锐,等. 参考作物腾发量的GARCH类模型模拟与比较[J]. 农业工程学报201531(7):131-136.
12
江显群,陈武奋,邵金龙 .基于公共天气预报的参考作物腾发量预报[J].排灌机械工程学报201937(12):1 077-1 081.
13
SLAVISA T BRANIMIR T MIOMIR S. Fore-casting of Reference Evapotranspiration by Artificial Neural Networks. 2003129(6):454-457.
14
崔远来,马承新,沈细中, 等. 基于进化神经网络的参考作物腾发量预测[J]. 水科学进展2005(1):76-81.
15
徐俊增,彭世彰,张瑞美, 等. 基于气象预报的参考作物蒸发蒸腾量的神经网络预测模型[J]. 水利学报2006(3):376-379.
16
中国气象总局. 地面气象观测标准[M]. 北京:气象出版社,2003.
17
杨洋. 基于公共天气预报的参考作物腾发量预报模型比较[D]. 武汉:武汉大学,2018.
18
HARGREAVESGEORGE H ALLENRICHARD G. History and Evaluation of Hargreaves Evapotranspiration Equation. 2003129(1):53-63.
19
YANG Y CUI Y L LUO Y F, et al. Short-term forecasting of daily reference evapotranspiration using the Penman-Monteith model and public weather forecasts[J]. Agricultural Water Management2016,177.
PDF(2820 KB)

740

访问

0

引用

详细情况

段落导航
相关文章

/