[發明專利]一種中小流域產匯流參數自動率定方法

在審

• 申請號: 201810206472.0
• 申請日: 2018-03-13
• 公開（公告）號: CN108491605A
• 公開（公告）日: 2018-09-04
• 發明（設計）人: 趙玲玲;張豫;蔡德所;陳彩霞;張瑞虎;溫美麗;何奕霏
• 申請（專利權）人: 廣州地理研究所
• 主分類號: G06F17/50
• 分類號: G06F17/50;G06Q10/06
• 代理公司: 南寧東智知識產權代理事務所(特殊普通合伙) 45117
• 代理人: 巢雄輝;汪治興

本發明屬于產匯流參數率定的技術領域，尤其涉及一種中小流域產匯流參數自動率定方法，針對直接應用CN值查算表進行產匯流計算誤差大的問題，現提出以下方案，包括以下步驟：選取待測流域，通過水文手冊查算統計該流域歷時的面雨量、徑流量和前期影響雨量，對該流域面雨量、徑流量和前期影響雨量的歷時數據進行篩選采樣，選取統計某一年中每個月的面雨量、徑流量和前期影響雨量。本發明選取典型的面雨量、徑流量和前期影響雨量的數據，進行分析評估，將采樣的面雨量、徑流量和前期影響雨量的數據輸入到編程中，編程按照特定的公式進行優選，自動優選出最優的產匯流參數，并選擇某一組數據對參數進行驗證，保證產匯流參數率定的準確性。

技術領域

本發明涉及產匯流參數率定的技術領域，尤其涉及一種中小流域產匯流參數自動率定方法。

背景技術

流域的產匯流計算是研制水文預報方法、解決許多水文水資源實際問題的重要依據，所謂產流，是指流域中各種徑流成分的生成過程，流域的產流計算主要是根據產流理論通過不同的方法構建數學模型來計算流域的產流量，主要的計算方法有降雨徑流相關法，濕潤地區的蓄滿產流模型和干旱地區的超滲產流模型，流域匯流是降落在流域上的雨水，從流域各處向流域出口斷面匯集的過程，它包括地表水、地下水的水流匯集，在對中小流域的產匯流計算時，常需要用到產匯流參數進行計算，二產匯流參數隨著不同地區的自然地理情況和水文條件的差異而變化，傳統計算過程中直接應用CN值查算表進行產匯流計算，計算誤差大。

發明內容

基于背景技術存在的技術問題，本發明提出了一種中小流域產匯流參數自動率定方法。

本發明提出的一種中小流域產匯流參數自動率定方法，包括以下步驟：

S1：選取待測流域，通過水文手冊查算統計該流域歷時的面雨量、徑流量和前期影響雨量；

S2：對該流域面雨量、徑流量和前期影響雨量的歷時數據進行篩選采樣，選取統計某一年中每個月的面雨量、徑流量和前期影響雨量；

S3：判斷采樣數據是否是典型情況，通過該流域的歷時文獻和流域資料對面雨量、徑流量和前期影響雨量的采樣數據進行分析評估；

S4：將該流域統計的面雨量、徑流量和前期影響雨量的數據曲線圖繪制在紙上，進行分析；

S5：將面雨量、徑流量和前期影響雨量的采樣數據通過編程按照特定的公式自動優選出最優的參數，完成產匯流參數的率定；

S6：選取歷史數據中的一組面雨量、徑流量和前期影響雨量的數據，將優選出的產匯流參數代入相應的關系函數檢驗，進行分析評估。

優選地，所述S1中面雨量是由泰森多邊形加權平均求該流域上的面平均雨量，徑流量是由徑流總量除以流域面積。

優選地，所述S1中前期影響雨量的計算公式為：

其中P為面雨量，R為徑流量，n為降雨前天數取10，K為土壤含水率的日消退系數，EM為流域日蒸發能力，按照不同的月份選用相應的平均值，WM為流域最大需水量，為前期影響雨量。

優選地，所述S4中面雨量、徑流量和前期影響雨量的數據曲線圖繪制時，應注意三條曲線之間不能相交，且三條曲線繪制的距離要適中。

優選地，所述S5中的特定公式如下：

其中為采樣數據中徑流量系列的平均值，是參數m的函數，越大，計算結果越優，約束條件一般取0.10<m<0.25，通過編程對目標函數進行優化計算可以確定m。

優選地，所述S6中關系函數的公式如下：

其中S為待測流域當時最大缺水量。