[發明專利]一種水文集合預報后處理方法

在審

• 申請號: 201810974789.9
• 申請日: 2018-08-24
• 公開（公告）號: CN109325209A
• 公開（公告）日: 2019-02-12
• 發明（設計）人: 王國強;王溥澤;姚繼平;韓子叻;朱虹霖
• 申請（專利權）人: 北京師范大學;山東省環境規劃研究院
• 主分類號: G06F17/18
• 分類號: G06F17/18;G06F17/16
• 代理公司: 北京市廣友專利事務所有限責任公司 11237
• 代理人: 張仲波

1.一種水文集合預報后處理方法，其特征在于，包括：

步驟1、獲取概率密度函數

式中，Q_o是觀測流量，Q_s是模擬流量，Q_f是預報流量，d為預設參數；

步驟2、獲取歷史觀測流量Q_o,a、歷史模型模擬流量Q_s,a、模型預報流量Q_s,f；當概率密度函數服從高斯正態分布時，則求解上式以獲取總線性后處理模型GLMPP：

步驟3、利用總線型模型求出預報時期的未來觀測流量Q_o,f。

2.根據權利要求1所述的水文集合預報后處理方法，其特征在于，所述步驟2還包括：

當流量不服從正態分布時，通過正態轉換將流量矩陣轉換成正態分布：

令Z_1,2＝Z₁|Z₂為由預報矢量得到的預報量，GLMPP模型可用下式表達：

Z_1,2＝A·Z₂+B·E E＝N(0,1)

式中，A，B是參數矩陣，E是標準正態分布；

令則Z的協方差矩陣為：

則求解出參數矩陣

3.根據權利要求2所述的水文集合預報后處理方法，其特征在于，所述步驟3具體包括：

步驟31、構建待預測時間段的觀測矩陣和預報矩陣；

Q₁＝[Q_o,f(N_y×N_b,N_f)]

Q₂＝[Q_s,f(N_y×N_b,N_f)；Q_o,a(N_y×N_b,N_a)；Q_s,a(N_y×N_b,N_a)]；

其中N_y為預報與觀測數據的年數；N_b為待預測的時間段；N_a為N_b之前的預設數量個時間段；N_f為N_b之后的預設數量個時間段；

步驟32、計算預報矢量得到的預報量Z_1,2，然后通過投影轉換將預報量Z_1,2轉換成未來觀測流量矩陣[Q_o,f]。