王洪心　黃新平　王秋碩

【摘要】河溶站是湖北省重要防汛站，洪峰水位預(yù)報(bào)是防汛關(guān)鍵要素，由于其位于平原地區(qū)，坡降平緩，水位流量關(guān)系在中、高水呈連時(shí)序（主要以順時(shí)針繩套線為主），在利用降雨徑流關(guān)系預(yù)報(bào)出洪峰流量后，預(yù)報(bào)洪峰水位存在瓶頸，有時(shí)誤差較大，本文從實(shí)時(shí)校正出發(fā)，研究洪峰水位預(yù)報(bào)方法，提高預(yù)報(bào)精度，為防汛當(dāng)好參謀作用。

【關(guān)鍵詞】最大流量；最高水位；漲水面；實(shí)時(shí)訂正；預(yù)報(bào)水位

1、基本概況

1.1流域概況

沮漳河是長江中游上段干流左岸一級(jí)支流，分東、西兩支，東支為漳河，西支為沮河，西支較長為干流，發(fā)源于保康縣歐家店大灣，自北向南流經(jīng)保康、南漳、遠(yuǎn)安、當(dāng)陽、枝江等縣市，于荊州市李埠鎮(zhèn)臨江寺入長江，1993年，實(shí)施沮漳河下游出口段改道，河口上移15.5km，減少了荊江大堤李埠以下14.5km與沮漳河共堤，在枝江市七星臺(tái)鎮(zhèn)鴨子口村匯入長江。流域全長344km，集水面積7305km2。在當(dāng)陽市玉泉鎮(zhèn)北門有百里長渠分流灌溉，于枝江市問安鎮(zhèn)鳳臺(tái)回歸沮漳河。下游兩岸有堤防，左岸堤自棗林崗與荊江大堤連，因其泄洪能力較小，若與長江洪水遭遇，可威脅荊江大堤，歷史上曾多次自然潰口和破垸行洪。

西支沮河于當(dāng)陽市河溶鎮(zhèn)兩河口與漳河匯合，2015年，實(shí)施河溶鎮(zhèn)附近漳河改道，兩河交匯口由原河溶大橋上移約9km 至張家灣，全長230km，集水面積3353 km2，東支漳河源出?？悼h龍坪鄉(xiāng)黃龍洞溝，自北向南流經(jīng)漳河水庫，途經(jīng)南漳、荊門、遠(yuǎn)安縣，于當(dāng)陽市河溶鎮(zhèn)兩河口匯入沮漳河，全長190km，集水面積2968 km2。

1.2測站特性

河溶水文站位于沮漳河原兩河口下游800m處，集水面積6338 km2（含漳河水庫面積2212 km2），主要觀測項(xiàng)目是水位、流量?；舅哂^測采用自記水位計(jì)，測驗(yàn)河段兩岸為堤防控制，順直長度約500m，主槽偏于左岸，最大河寬450m，當(dāng)水位超過47.70m時(shí)，右岸開始漫灘，水位達(dá)48.80m時(shí)，右岸全部漫灘直達(dá)堤腳。河溶站屬平原河流站，坡降平緩，流速較小，Z～Q關(guān)系受洪水漲落、沖淤變化和長江回水頂托等影響，水位流量關(guān)系復(fù)雜，洪水期線型多為逆時(shí)針繩套狀或其變形，歷年采用連時(shí)序法推流，該站設(shè)防水位為47.50m，警戒水位為48.50m，保證水位為50.00m，斷面附近堤頂高程為52.80m。

河溶站多年平均降水量995.0mm，汛期（4～9月）降水量約占全年的75%，降水空間分布均勻，過程較為集中。年內(nèi)暴雨洪水主要出現(xiàn)在汛期，尤以5-7月為最多；流域呈平行狀，坡降小，匯流慢。

河溶站歷年實(shí)測最高水位50.49m，歷年最大流量2360 m3.s-1，歷年實(shí)測最低水位39.86m，歷年最小流量-5.00 m3.s-1，實(shí)測最大流速2.98 m.s-1。年水位變幅最大可達(dá)10 m。

2、洪峰流量預(yù)報(bào)洪峰水位模型建立

該站洪水由暴雨形成，其洪峰流量采用降雨徑流預(yù)報(bào)，精度能滿足要求。由于水位流量關(guān)系為繩套或其變形，通過洪峰流量確定洪峰水位存在較大誤差，而洪峰水位是防汛指揮的關(guān)鍵因素。本文通過實(shí)時(shí)校正法，對(duì)提高水位預(yù)報(bào)精度進(jìn)行探討。

2.1相同流量不同水位分析

選取洪峰流量接近，但洪峰水位相差較大的典型洪水場次，如表1。

2.2由洪峰流量預(yù)報(bào)洪峰水位

2.2.1繪制最大流量～最高水位綜合線

逆時(shí)針繩套曲線各特征值出現(xiàn)順序?yàn)樽畲蟊冉?、最大流速、最大流量、最高水位。選取2002-2016年洪峰水位高于44.00m以上洪水場次，由于最大流量與最高水位不是同時(shí)出現(xiàn)，現(xiàn)繪制洪峰最大流量與最高水位綜合關(guān)系線，此線用于解決水位流量在同場洪水過程中不是單一對(duì)應(yīng)關(guān)系，見圖3。

2.2.2繪制洪峰流量～相應(yīng)水位綜合線（漲水面）

選取2002-2016年典型洪水場次，繪制最大流量與其相應(yīng)水位綜合關(guān)系線（漲水面），此線用于對(duì)洪峰預(yù)報(bào)水位的實(shí)時(shí)修正，解決相同的洪峰流量，在不同的場次時(shí)，其洪峰水位差異大的問題，見圖3。

2.2.3洪峰流量預(yù)報(bào)洪峰水位步驟

由預(yù)報(bào)洪峰流量查最大流量～最高水位綜合線得Z初估作為最高水位預(yù)報(bào)初估值；由預(yù)報(bào)流量查最大流量～相應(yīng)水位綜合線得Z查線作為漲水面最大流量對(duì)應(yīng)的標(biāo)準(zhǔn)水位；以實(shí)時(shí)最大流量相應(yīng)水位Z實(shí)時(shí)為依據(jù)，計(jì)算實(shí)時(shí)預(yù)報(bào)水位訂正值Z訂正=（Z實(shí)時(shí)- Z查線）/2；最后確定最高水位預(yù)報(bào)值Z修正= Z初估+ Z訂正。

3、實(shí)例運(yùn)用

3.1由最大流量預(yù)報(bào)最高水位方法

以2007年8月31日?qǐng)龃魏樗疄槔?，說明由最大流量（實(shí)例以實(shí)測值代替預(yù)報(bào)洪峰流量）預(yù)報(bào)最高水位的方法如下，從表1可知：本場洪水最大流量1190 m3.s-1，最高水位Zmax=47.77m。

（1）由最大流量查圖3中1線，初估最高水位預(yù)報(bào)值Z初估=48.25m。

（2）由最大流量查圖3中2線，得綜合線上最大流量對(duì)應(yīng)水位（漲水面）Z查線=47.53m。

（3）在洪水水文要素摘錄表中，查得洪峰流量出現(xiàn)時(shí)，實(shí)際相應(yīng)水位Z實(shí)時(shí)=47.00m。

（4）計(jì)算最高水位預(yù)報(bào)訂正值，Z訂正=（Z實(shí)時(shí)- Z查線）/2=（47.00-47.53）/2=-0.265m。

（5）實(shí)時(shí)修正最高水位預(yù)報(bào)值Z修正= Z初估+ Z訂正=48.25-0.265=47.98m。

（6）最高水位預(yù)報(bào)誤差δz= Z修正- Zmax=47.98-47.77=0.21m。

3.2檢驗(yàn)

1973年以來，洪峰水位高于46.50m的洪水共計(jì)46場，按照上述方法進(jìn)行洪峰水位預(yù)報(bào)，根據(jù)《水文情報(bào)預(yù)報(bào)規(guī)范》規(guī)定，以0.30m作為水位預(yù)報(bào)許可誤差的上限，直接采用圖3中1線（最大流量最高水位關(guān)系線）預(yù)報(bào)水位合格率為54.3%；經(jīng)實(shí)時(shí)校正后，有42場洪峰水位預(yù)報(bào)合格，合格率為91.3%，達(dá)甲級(jí)精度等級(jí)，由于篇幅限制，水位預(yù)報(bào)誤差統(tǒng)計(jì)表略。

4、結(jié)論

根據(jù)實(shí)時(shí)最大流量、相應(yīng)水位對(duì)最高預(yù)報(bào)水位進(jìn)行修正，能顯著提高最高水位預(yù)報(bào)精度，河溶站防汛關(guān)鍵是需要最高水位值，想辦法提高水位預(yù)報(bào)精度尤其重要。經(jīng)對(duì)洪峰過程統(tǒng)計(jì)，該站最大流量與最高水位出現(xiàn)時(shí)間間隔在4～12h之間，若要獲得更長預(yù)見期，可采用最大流量出現(xiàn)前的實(shí)時(shí)流量值（小于最大流量10%左右）進(jìn)行修正。根據(jù)上述方法，實(shí)時(shí)流量需要測量、分析計(jì)算、查線獲取；另外，最大流量預(yù)報(bào)精度直接影響最高水位的預(yù)報(bào)精度；在受長江高水回水頂托影響期間，需按回水頂托規(guī)律進(jìn)行修正，不適用上述方法分析。