大伙房流域降雨徑流模型

1、大伙房流域降雨徑流模型大伙房流域降雨徑流預報模型又簡稱“DHF”模型，該模型于1973年由遼寧省大伙房水庫管理局劉愛杰、王本德等人提出，至今已使用30余年，為水庫洪水調度做出了很大貢獻?！癉HF”模型是適用于我國濕潤地區(qū)的超滲產(chǎn)流模型，目前已在遼寧省多個水庫的水情自動測報系統(tǒng)中使用，效果較好。建立在“DHF”模型基礎上的降雨徑流預報方案，在大伙房流域經(jīng)過調試和精度驗證后進行使用，在使用中平均精度令人滿意。尤其在“957”特大洪水調度中，發(fā)揮了顯著作用，準確預報出了第一非常溢洪道溢流時間和水庫最高庫水位，為省防制定調度決策提供了科學依據(jù)，使水庫工程發(fā)揮了強大的調蓄作用，最大限度地配合了下游搶險，

2、共減免下游直接經(jīng)濟損失74.89億元。“DHF”模型由兩部分組成，一是八參數(shù)超滲產(chǎn)流計算模型，引用雙層入滲曲線進行扣損計算，并以拋物線描述表層蓄水量和下層滲率的分布狀況；二是八參數(shù)變強度、變速度的經(jīng)驗單位線匯流計算模型，參數(shù)隨降雨分布而變，采用“前期影響凈雨”描述匯流速度的變化。這是一個集總的概念模型，模型的參數(shù)多半在滿足其物理意義的前提下確定，只有6個需要優(yōu)選法選定或試錯法確定。1 大伙房模型概化流程流域下墊面分為表層、下層和地下水蓄存三部分，計算流程如圖10-1所示。2 大伙房模型產(chǎn)流計算產(chǎn)流模型將下墊面分為表層，下層和深層三部分。表層土壤中的張力水蓄量與植物截流、填洼儲存合稱表層蓄水量，

3、其極值為表層蓄水容量；下層土壤中的張力水蓄量稱為下層蓄水量，其極值為下層蓄水容量；地下水儲水層的蓄水量以表示，其極值為地下水庫蓄水容量。圖10-1 大伙房產(chǎn)流模型概化流程圖其中：S0-表層蓄水容量；U0-下層蓄水容量；D0-下層蓄水容量（U0）與地下水庫蓄水容量（V0）之和；g-不透水面積占全流域面積之比值；Kc-流域蒸散發(fā)能力與大水體蒸發(fā)量的比值；A-表層蓄水容量曲線形狀參數(shù)；B-下滲率拋物線分配曲線形狀參數(shù)；K-下層下滲曲線曲率；2.1 不透水面積上的徑流計算 根據(jù)上面的流程圖，不透水面積上直接產(chǎn)生徑流，通過下式來計算： （1） （2）式中P為時段降雨量，KC是流域蒸散發(fā)能力與ER大水體蒸

4、發(fā)量的比值，C為日降雨觀測的時段數(shù)，g為流域不透水面積占總流域面積的比值。2.2 地面壤中流的計算地面壤中流通過下式計算 （3）2.2.1 表層扣損R的計算當時， （4）當時， （5）式中Sa是表層流域平均蓄水量，Sm是相應于Sa的點蓄水量；a為曲線形狀參數(shù)。Sm和Sa由下式計算得到： （6） （7）式中為表層點最大蓄水量，當時，則，并有 （8）其中通過下式來計算： （9）2.2.2 下滲計算時段平均總下滲率與的關系: (10)時段平均下層下滲率與的關系: (11)其中，和的關系如下: (12) 和分別是反映總下滲和下層下滲曲線曲率系數(shù)；是的函數(shù)，二者的關系式如下: (13)是下層蓄水容量與地

5、下水庫蓄水容量之和。2.2.3 徑流深計算公式總徑流y的計算公式： (14)當時，得 (15)地面壤中流的計算公式：當時， (16)當時， (17)和通過下列兩式計算 (18) (19)地下徑流的計算： (20)其中為與之比值，在干旱地區(qū)全部滲入深層，不形成地下徑流，則；在濕潤地區(qū)全部形成地下徑流流入河道，所以。為地下水庫下滲強度，為總徑流與壤中流之差，即： (21)2.3 大伙房模型蒸散發(fā)計算（1）當降雨量超過流域蒸發(fā)能力時，即凈滲雨強，產(chǎn)生下滲，土壤蓄水量增加，此時：表層蓄水量計算公式 (22)式中t為單位日，以8:00為日分界；以式（6）計算當時，上式可寫成 (23)依照流程圖的下滲程序

6、，逐時段計算式， (24)式中i為時段序；、分別用（9）、（4）、（5）式進行計算。下層土壤蓄水計算公式 據(jù)流程圖，逐時段計算的公式 (25)若逐日計算，可用下式計算若時 (26)若時 (27)式中的與分別用式（19）及求出。（2）當日流域平均降雨量小于流域蒸發(fā)能力，即時，土壤蓄水量將減少進行蒸散發(fā)計算表層蒸散發(fā)公式假定：a.表層蒸發(fā)等于流域平均蒸散發(fā)能力；b.每次下滲后的初次蒸散發(fā)皆從地表開始。為表層可蒸發(fā)量，計算公式為： (28)式中： (29)因此，日表層蓄水量計算式為 (30)下層蓄水量計算公式假定：a.表層蓄水量蒸發(fā)盡后，該點的下層開始蒸發(fā)；b.下層蒸發(fā)量與下層蓄水量成比例關系。因此

7、， (31)當時，則 (32)這樣，日下層土壤蓄水量以下式計算 (33)2.4 大伙房模型匯流計算大伙房匯流模型是根據(jù)多年預報經(jīng)驗，總結出的8參數(shù)變強度、變匯流速度的經(jīng)驗單位線，其線型公式為式中 匯流曲線，是與t的函數(shù)； 匯流曲線底寬；t 匯流時間（時刻）；K3 匯流曲線的比例系數(shù)；n 地下徑流與地面壤中流匯流曲線底寬的比例系數(shù)；AA、DD、CC反映匯流曲線形狀的參數(shù)；地表徑流與地下徑流具有不同的匯流曲線形狀，可用形狀參數(shù)區(qū)分；不同的降雨分布地表徑流也具有不同的匯流曲線，其參數(shù)也不同。是時段凈雨（徑流）與前期影響凈雨（已產(chǎn)生的徑流）的函數(shù)，模型中假定：式中 特征河長 ；流域最大河長；擬合線的斜

8、率；前期影響凈雨量，表征河槽蓄水狀況，制約匯流速度等；前期凈雨影響程度衰減系數(shù)，反映已產(chǎn)生的徑流在河槽中消退速度；2.5 大伙房模型產(chǎn)流參數(shù)優(yōu)選大伙房模型產(chǎn)流參數(shù)的最優(yōu)化問題是多目標的有約束的非線性規(guī)劃問題。目標函數(shù)取為N次洪水產(chǎn)流的擬合精度和N次洪水計算產(chǎn)流與實際產(chǎn)流誤差總和的線性加權，各取0.5的權重進行求解。將問題變?yōu)椋菏街?ACU1-N次洪水產(chǎn)流的擬和精度 TYTC-N次洪水計算產(chǎn)流與實際產(chǎn)流誤差總和 其中N為洪水的次數(shù)，TOTAL為N次洪水的實際凈雨的總和，TYTC為N次洪水的預報凈雨的總和， N合格為擬合的合格洪水的次數(shù)。若某次洪水滿足以下條件之一即為合格：計算產(chǎn)流與實際產(chǎn)流的絕對

9、誤差小于3mm；計算產(chǎn)流與實際產(chǎn)流的相對誤差小于20%且絕對誤差小于20mm；計算時判斷合格的絕對誤差和相對誤差的標準可以根據(jù)具體情況來設定；本程序采用可變容差法求解問題，把多個約束求極小值問題變?yōu)橐粋€單約束求極小值問題。設所求解的非線性規(guī)劃問題為：本算法將上述兩問題化為單約束問題進行求解:式中T(x)是約束破壞的估計量，表達式為：其中，Uj為Heaviside算子。是第K步搜索中給出的關于可行性的可變容差準則值，表達式為：其中 t-初始多面體的大小; m-等式約束數(shù)目;r=n-m，的自由度;-多面體中第K次搜索時的第i個頂點值;-多面體（除）外的形心頂點;由于各參數(shù)的初值相差較大，而且變化范

10、圍也各不相同，所以將參數(shù)歸一化，即初值都取為1.0。根據(jù)物理意義將E6、E7、E8、E9的變化范圍取為0.5-10.0，其它參數(shù)取為0.5-5.0。2.6 大伙房模型匯流參數(shù)優(yōu)選大伙房匯流模型需要優(yōu)選的參數(shù)有B0、K0、K、n、DD、CC，其中DD和CC根據(jù)降雨的分布分為上游、中游和下游三套參數(shù)，所以需要優(yōu)選的參數(shù)一共有10個。K3和AA可以通過優(yōu)選出的參數(shù)帶入下列兩式計算求得：大伙房匯流參數(shù)優(yōu)選的基礎是多次歷史洪水，其最優(yōu)化問題是多目標的有約束的非線性規(guī)劃問題。目標函數(shù)取為N次洪水過程的擬合精度、N次洪水計算洪峰與實際洪峰誤差總和以及N次洪水的計算峰現(xiàn)時間和實際的峰現(xiàn)時間誤差總和的線性加權，

11、各取0.4、0.3、0.3的權重進行求解。將問題變?yōu)椋菏街?QE1-N次洪水匯流的擬和精度; QE2-N次洪水計算峰值與實際峰值誤差總和; QE3-N次洪水計算峰現(xiàn)時間與實際峰現(xiàn)時間誤差總和;本程序采用可變容差法求解問題，把多個約束求極小值問題變?yōu)橐粋€單約束求極小值問題.設所求解的非線性規(guī)劃問題為:本算法將上述兩問題化為下式進行求解:式中T(x)是約束破壞的估計量，表達式為:其中，Uj為Heaviside算子。是第K步搜索中給出的關于可行性的可變容差準則值，表達式為：其中 t-初始多面體的大小; m-等式約束數(shù)目;r=n-m，的自由度;-多面體中第K次搜索時的第i個頂點值;-多面體（除）外的形

12、心頂點;2.7 大伙房產(chǎn)流模型參數(shù)為提高產(chǎn)流預報精度，按降雨開始時的起始流量大小進行分類。分類預報參數(shù)如表10-1表10-2所示。表10-1 降雨起始流量³100m3/s的大伙房產(chǎn)流模型參數(shù)gABKS0U0D00.054.1291.1810.2845100130續(xù)表10-1 流域蒸發(fā)能力Ed月份6789日雨量 < 5mm3.74.34.33.75mm日雨量 < 30mm1.91.91.91.9日雨量30mm1.01.01.01.0表10-2 降雨起始流量<100m3/s的大伙房產(chǎn)流模型參數(shù)gABKS0U0D00.054.1291.3710.46980100140續(xù)表10-2 流域蒸發(fā)能力Ed月份6789日雨量 < 5mm3