薛海波
(遼寧省柴河水庫管理局有限責(zé)任公司,遼寧 鐵嶺 112000)
水文模型是為模擬水循環(huán)過程而構(gòu)建,是水文科學(xué)研究的重要手段與方法之一,流域水文模擬經(jīng)歷了由系統(tǒng)模型到概念性模型再到物理模型、由集總式模擬到分布式模擬的發(fā)展歷程,每種模型都有各自的理論與應(yīng)用特點(diǎn)。對于一個(gè)特定流域,如何選擇最適合的水文模型是水文過程模擬的關(guān)鍵。為此,本文選定遼寧省柴河水庫以上流域?yàn)檠芯繉ο?,?yīng)用目前在遼寧地區(qū)具有代表性的大伙房(DHF)模型和遼寧指數(shù)模型水進(jìn)行洪水預(yù)報(bào),并對2種模型進(jìn)行了分析比較。
柴河起源于遼寧省清原縣北樂山岔,是遼河干流左側(cè)的一個(gè)較大型支流,柴河水庫控制面積1355km2,水庫以上河流長度123.5km。柴河流域降雨年際、年內(nèi)、場次分配不均。流域多年平均降雨量為737mm。系列最大值1185.2mm(2010年),最小值為493.6mm(1997年),其比值為2.4。全年降雨主要集中在6—9月份,占到全年降雨總量的75.0%.
大伙房洪水預(yù)報(bào)模型(簡稱為DHF模型)于1973年由遼寧省大伙房水庫管理局提出。在東北地區(qū)得到了廣泛應(yīng)用,并取得良好的效果。DHF模型分為產(chǎn)、匯流2部分,產(chǎn)流部分是根據(jù)超滲產(chǎn)流原理,采用雙層入滲曲線進(jìn)行扣損計(jì)算;匯流部分是變強(qiáng)度、變速度的經(jīng)驗(yàn)單位線匯流計(jì)算模型。
遼寧指數(shù)模型是由遼寧省水文總局提出的一種基于蓄滿產(chǎn)流模式的水文預(yù)報(bào)模型,參數(shù)簡單、效果較好,在遼寧地區(qū)大范圍使用。遼寧指數(shù)模型共由蓄滿產(chǎn)流和經(jīng)驗(yàn)單位線匯流2個(gè)部分組成。其中有4個(gè)蓄滿產(chǎn)流模型參數(shù),利用柴河水庫降雨徑流相關(guān)曲線計(jì)算產(chǎn)流;使用綜合單位線進(jìn)行匯流模擬結(jié)果的計(jì)算。
DHF模型在超滲產(chǎn)流預(yù)報(bào)基礎(chǔ)之上,建立起經(jīng)驗(yàn)單位線匯流模型計(jì)算公式如下:
(1)
由于我國北方流域具有典型的特點(diǎn),參數(shù)率定過程需在DHF模型中引入蒸發(fā)計(jì)算,選擇合適的取值區(qū)域,以保證參數(shù)真實(shí)可靠。利用人工優(yōu)選與遺傳算法相互結(jié)合的方法來進(jìn)行反復(fù)迭代,最終選取出最優(yōu)參數(shù)值,具體見表1。相應(yīng)的蒸發(fā)參數(shù)見表2。
表1 大伙房模型產(chǎn)流參數(shù)率定結(jié)果
表2 大伙房模型蒸發(fā)參數(shù)
遼寧指數(shù)模型以蓄滿產(chǎn)流為機(jī)理,產(chǎn)流的計(jì)算公式如下:
(2)
(1)降雨量、流量資料采用1975—2013年系列,其中1981年以后入庫流量資料采用水量平衡法反推求得。
(2)流域平均雨量采用各雨量站算術(shù)平均值或加權(quán)平均值。
(3)徑流深計(jì)算為入庫洪水過程割去基流和前期退水過程,用累積法計(jì)算。
(4)流域內(nèi)蒸發(fā)量采用清原、鐵嶺站資料。
(5)流域土壤蓄水量,流域蓄水量采用雙層扣損法計(jì)算。從5月1日起算,初始值為60mm。
(6)降雨徑流相關(guān)圖。鑒于原參數(shù)的準(zhǔn)確性與穩(wěn)定性,仍采用柴河水庫原參數(shù)值。
利用矩法和極大似然法,并參考柴河水庫的經(jīng)驗(yàn)來進(jìn)行最終的參數(shù)優(yōu)選,經(jīng)過反復(fù)迭代得到率定結(jié)果見表3—4。結(jié)合柴河水庫在實(shí)際運(yùn)用中的經(jīng)驗(yàn)總結(jié),所選參數(shù)取值均在在合理范圍內(nèi),認(rèn)為參數(shù)值合理可靠。
表3 遼寧指數(shù)模型產(chǎn)流參數(shù)率定結(jié)果
表4 遼寧指數(shù)模型蒸發(fā)參數(shù)
依據(jù)GB/T 22482—2008《水文情報(bào)預(yù)報(bào)規(guī)范》的要求,判斷水庫洪水預(yù)報(bào)精度。
根據(jù)DHF模型和遼寧指數(shù)模型參數(shù)列表,采用選定的25場歷史洪水?dāng)?shù)據(jù)模擬柴河水庫的產(chǎn)流情況,模擬結(jié)果匯總見表5。
表5 DHF模型和遼寧指數(shù)模型產(chǎn)流結(jié)果
根據(jù)表格中的計(jì)算結(jié)果,DHF模型洪水產(chǎn)流計(jì)算合格率僅為64%,未達(dá)到甲級(jí)預(yù)報(bào)水平,模擬結(jié)果不理想。預(yù)報(bào)凈雨值比實(shí)測凈雨值平均偏大約6.13mm,其中存在的誤差并非是偶然誤差,而是系統(tǒng)誤差。
根據(jù)圖表中的數(shù)據(jù),遼寧指數(shù)模型產(chǎn)流預(yù)報(bào)結(jié)果合格率為92%;預(yù)報(bào)凈雨值與實(shí)測凈雨值較為接近,平均偏大約1.82mm,模擬預(yù)報(bào)效果較好。
基于DHF模型和遼寧指數(shù)模型模擬柴河水庫洪水實(shí)例,產(chǎn)流計(jì)算結(jié)果分析如下。
(1)DHF產(chǎn)流模型的預(yù)報(bào)凈雨大于實(shí)測凈雨,產(chǎn)流合格率較低僅為64%,說明存在系統(tǒng)誤差;不合格的9場洪水降雨前期土壤含水量均不高,其中6場洪水模擬凈雨值集中在20~50mm之間,實(shí)測洪峰流量和洪水總量均未超過5年一遇標(biāo)準(zhǔn)(19940807場次洪水的洪峰流量略小于5年一遇標(biāo)準(zhǔn)),說明大伙房模型不適用于柴河流域的小洪水產(chǎn)流預(yù)報(bào)。
(2)遼寧指數(shù)模型的預(yù)報(bào)凈雨與實(shí)測凈雨較為接近,預(yù)報(bào)合格率為92%,模擬預(yù)報(bào)效果較好。大洪水19950725產(chǎn)流誤差為1.5mm,19940816洪水產(chǎn)流誤差0.9mm,誤差均較小。預(yù)報(bào)不合格的場次洪水均為當(dāng)年第一場較大降雨,可能前期土壤缺水量大,實(shí)際產(chǎn)流量較小,導(dǎo)致預(yù)報(bào)值偏大較多。
基于以上產(chǎn)流計(jì)算結(jié)果,選擇使用遼寧指數(shù)模型進(jìn)行下一步的匯流模擬計(jì)算。匯流計(jì)算使用到的計(jì)算方法為謝爾曼時(shí)段單位線,按照雨強(qiáng)、雨量、時(shí)空分布等不同的分類指標(biāo)進(jìn)行匯流的模擬預(yù)報(bào)。計(jì)算結(jié)果見表6。
表6 遼寧指數(shù)模型匯流結(jié)果
上述匯流結(jié)果表明。
(1)對于頻率大于5年一遇標(biāo)準(zhǔn)的洪水,8場次中有7場預(yù)報(bào)結(jié)果合格,合格率為87.5%;19950725場次洪水降雨較大,達(dá)到了100年一遇洪水標(biāo)準(zhǔn),洪水產(chǎn)匯流預(yù)報(bào)誤差較小,洪峰流量預(yù)報(bào)值3332m3/s接近實(shí)際值3310m3/s,峰現(xiàn)時(shí)間無誤差,使用遼寧指數(shù)模型模擬柴河流域大洪水取得良好效果。
(2)對于頻率小于5年一遇標(biāo)準(zhǔn)的洪水,17場次中有13場預(yù)報(bào)結(jié)果合格,合格率為76.5%,雖然滿足預(yù)報(bào)要求但沒有達(dá)到理想的效果。產(chǎn)生這種結(jié)果的原因主要是柴河水庫作為多年調(diào)節(jié)的大型水庫,調(diào)蓄洪水能力比較強(qiáng),小洪水對流域整體以及水庫安全產(chǎn)生的影響較小。
(1)基于超滲產(chǎn)流的DHF模型不適用于柴河流域水庫洪水預(yù)報(bào)。原因?yàn)槟P痛嬖谙到y(tǒng)誤差,可能是由于原始資料誤差引起,也可能是在進(jìn)行計(jì)算的過程中選取的模型與流域擬合度不高而產(chǎn)生。
(2)基于蓄滿產(chǎn)流的遼寧指數(shù)模型適合模擬柴河流域水庫洪水預(yù)報(bào)。產(chǎn)流預(yù)報(bào)合格率達(dá)到92%,且降雨量級(jí)越大,預(yù)報(bào)精度越高;匯流預(yù)報(bào)合格率達(dá)87.55%,產(chǎn)匯流預(yù)報(bào)均為甲級(jí)預(yù)報(bào)水平,較大洪水的洪量預(yù)報(bào)誤差均較小,滿足預(yù)報(bào)調(diào)度方式的要求。
(3)遼寧指數(shù)模型在小洪水的應(yīng)用中存在不足。小于5年一遇標(biāo)準(zhǔn)的小洪水預(yù)報(bào)中,雖然滿足預(yù)報(bào)要求但沒有達(dá)到理想的效果,需要在以后的研究中進(jìn)一步優(yōu)化模型參數(shù)。
(4)本文屬于對柴河水庫汛期洪水預(yù)報(bào)的探索性研究,后續(xù)可以在本文基礎(chǔ)上,就汛限水位動(dòng)態(tài)控制在水庫實(shí)時(shí)調(diào)度中的應(yīng)用問題進(jìn)行更加深入的研究。