[發(fā)明專利]一種基于主控因子修正的參照作物騰發(fā)量預(yù)測方法

在審

• 申請?zhí)枺?/td>: 202011473915.6
• 申請日: 2020-12-15
• 公開（公告）號: CN112507558A
• 公開（公告）日: 2021-03-16
• 發(fā)明（設(shè)計）人: 彭致功;宮儷芹;張寶忠;魏征;蔡甲冰;陳鶴;董潔;張倩;王雅琦;韓聰穎;謝天慧
• 申請（專利權(quán)）人: 中國水利水電科學(xué)研究院
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06F119/08
• 代理公司: 北京正華智誠專利代理事務(wù)所(普通合伙) 11870
• 代理人: 李林合

1.一種基于主控因子修正的參照作物騰發(fā)量預(yù)測方法，其特征在于，包括以下步驟：

S1、收集預(yù)測地PMT模型計算所需歷史實測氣象數(shù)據(jù)；

S2、將所述步驟S1收集的預(yù)測地歷史實測氣象數(shù)據(jù)代入PMT模型中計算出歷年參照作物騰發(fā)量；

S3、根據(jù)所述步驟S2得到的歷年參照作物騰發(fā)量采用基于多參數(shù)法確定PMT模型中的主控因子；

S4、收集預(yù)測地PMT模型計算所需天氣預(yù)報值；

S5、采用基于貢獻率的卡爾曼濾波法對所述步驟S4中收集的天氣預(yù)報值中的PMT模型主控因子預(yù)報值進行修正；

S6、將所述步驟S5中修正完成的PMT模型主控因子預(yù)報值和所述步驟S4收集的其它預(yù)測地天氣預(yù)報值代入PMT模型中進行計算，得到基于主控因子修正的參照作物騰發(fā)量預(yù)測結(jié)果。

2.如權(quán)利要求1所述的一種基于主控因子修正的參照作物騰發(fā)量預(yù)測方法，其特征在于：所述步驟S1具體包括收集過去30年以上的長系列預(yù)測地日最高氣溫、日最低氣溫和風(fēng)速的實測值。

3.如權(quán)利要求2所述的一種基于主控因子修正的參照作物騰發(fā)量預(yù)測方法，其特征在于，所述步驟S2包括以下分步驟：

S21、根據(jù)日最高氣溫和日最低氣溫計算出平均溫度T，計算公式為：

T＝(T_max+T_min)/2

其中，T_max為日最高氣溫，T_min為日最低氣溫；

S22、根據(jù)日最高氣溫和日最低氣溫計算出飽和水汽壓e_s，計算公式為：

S23、根據(jù)平均溫度T計算出飽和水汽壓曲線斜率Δ，計算公式為：

S24、根據(jù)日最低氣溫計算出實際水汽壓e_a，計算公式為：

S25、根據(jù)待測地緯度值和當(dāng)天序數(shù)J計算當(dāng)天天文輻射R_a，計算公式為：

其中，G_sc為太陽常數(shù)，δ為太陽傾角，d_r為日地相對距離，ω_s為日落時角；

根據(jù)日最高氣溫，日最低氣溫和天文輻射R_a計算出太陽短波輻射R_s，計算公式為：

R_s＝K_r(T_max-T_min)^0.5R_a

其中，K_r為地區(qū)調(diào)整系數(shù)；

根據(jù)日最高氣溫，日最低氣溫、天文輻射R_a和太陽短波輻射R_s計算出作物表面上凈輻射量R_n，計算公式為：

其中，Z為預(yù)測地高程；

S26、根據(jù)不同高程的風(fēng)速值計算出2m高處風(fēng)速值u₂，計算公式為：

其中，u_z為z米高處風(fēng)速值。

S27、根據(jù)所述步驟S21到S26的計算結(jié)果通過PMT模型計算歷年參照作物騰發(fā)量，計算公式為：

其中，ET₀為參照作物騰發(fā)量，G為土壤熱通量，γ為濕度計常數(shù)。