高金燕，史 輝，梁 博

（河北省桃林口水庫(kù)管理局，河北 秦皇島066400）

1 基本情況

桃林口水庫(kù)位于灤河主要支流青龍河上，是一座集城市供水、農(nóng)業(yè)灌溉為主，兼防洪和水力發(fā)電等綜合功能的大型水利樞紐工程，擔(dān)負(fù)著秦皇島市城市生活、工業(yè)、港口用水和灤河中下游地區(qū)農(nóng)田的灌溉任務(wù)。壩址至河源長(zhǎng)210km，壩址以上控制流域面積5060km2，占灤河流域面積的20.4%，多年平均降雨量500～700mm，多年平均徑流量9.60億m3，總庫(kù)容8.59億m3，以多年調(diào)節(jié)方式運(yùn)行。流域地處燕山山脈東段，屬東亞季風(fēng)氣候區(qū)，降雨量年內(nèi)及年際變化較大，徑流年內(nèi)集中，年際變化懸殊。形成暴雨的天氣系統(tǒng)主要為西風(fēng)槽和切變兩種，暴雨一般發(fā)生在每年的7月下旬至8月上旬，降雨量占全年的80%，是形成大洪水的主要原因，由于暴雨歷時(shí)短、強(qiáng)度大及地面坡陡流急，一次洪水過(guò)程多表現(xiàn)為峰高量大，陡漲陡落的形式。

根據(jù)水庫(kù)下游防洪需要，2013年對(duì)水庫(kù)調(diào)度運(yùn)用方式進(jìn)行了調(diào)整，由不承擔(dān)防洪任務(wù)調(diào)整為承擔(dān)20年一遇及其以下標(biāo)準(zhǔn)洪水的防洪任務(wù)，主汛期將汛限水位143.4m調(diào)整為138.5m控制運(yùn)用，防洪庫(kù)容2.04億m3，調(diào)洪庫(kù)容2.40億m3。

2 精細(xì)化降雨預(yù)報(bào)系統(tǒng)

為了進(jìn)步延長(zhǎng)洪水預(yù)見(jiàn)期，為優(yōu)化水庫(kù)調(diào)度創(chuàng)造條件，桃林口水庫(kù)引入了針對(duì)上游流域的精細(xì)化降雨預(yù)報(bào)系統(tǒng)。精細(xì)化降雨預(yù)報(bào)系統(tǒng)，利用中尺度WRF模式及快速循環(huán)同化技術(shù)，在對(duì)比試驗(yàn)的基礎(chǔ)上，選取了適合本地的積云過(guò)程、微物理過(guò)程和輻射過(guò)程的方案組合，以美國(guó)NECP（美國(guó)國(guó)家環(huán)境預(yù)報(bào)中心）的GFS全球預(yù)報(bào)資料作為背景場(chǎng)，提供預(yù)報(bào)時(shí)段0～48h的天氣形勢(shì)及多種物理量的精細(xì)預(yù)報(bào)。近些年來(lái)，隨著水庫(kù)水情信息采集站點(diǎn)設(shè)備的更新，精細(xì)化降雨預(yù)報(bào)系統(tǒng)充分利用了水情遙測(cè)信息，并在此基礎(chǔ)上分析開(kāi)發(fā)，采用先進(jìn)的地圖發(fā)布平臺(tái)、XML技術(shù)、Flex開(kāi)發(fā)技術(shù)和GIS圖表化表達(dá)技術(shù)，實(shí)現(xiàn)了實(shí)時(shí)水情信息與空間數(shù)據(jù)一體化存儲(chǔ)，并通過(guò)水情站與預(yù)警區(qū)域的關(guān)聯(lián)，實(shí)現(xiàn)了水雨情信息的地理網(wǎng)格化區(qū)域預(yù)警，有效提高了水庫(kù)應(yīng)對(duì)暴雨洪水的能力，提高了洪水的預(yù)見(jiàn)期。

3 精細(xì)化系統(tǒng)的組成與功能

桃林口水庫(kù)精細(xì)化降雨預(yù)報(bào)系統(tǒng)是在GIS系統(tǒng)的基礎(chǔ)上開(kāi)發(fā)的，系統(tǒng)由實(shí)時(shí)雨量數(shù)據(jù)提取、雨量預(yù)報(bào)數(shù)據(jù)傳輸、 面雨量水量計(jì)算和WEBGIS平臺(tái)顯示4部分組成。

3.1 實(shí)時(shí)雨量數(shù)據(jù)

用于提取各數(shù)據(jù)源的實(shí)時(shí)雨量，包括遙測(cè)雨量、氣象雨量等，將各數(shù)據(jù)源雨量進(jìn)行插補(bǔ)，以最大限度地增加觀測(cè)密度。雨量預(yù)報(bào)數(shù)據(jù)傳輸系統(tǒng)，將降雨預(yù)報(bào)數(shù)據(jù)（0～48h數(shù)值預(yù)報(bào)）從氣象部門(mén)傳輸?shù)剿畮?kù)管理局的服務(wù)器中。

3.2 面雨量水量計(jì)算系統(tǒng)

能夠計(jì)算水庫(kù)流域整個(gè)地區(qū)或任意區(qū)域的實(shí)時(shí)面雨量和預(yù)報(bào)面雨量，并提供預(yù)警報(bào)警功能，當(dāng)某區(qū)域的實(shí)際降雨或預(yù)報(bào)降雨超過(guò)一定閾值，或者未來(lái)0～48h預(yù)測(cè)降雨超過(guò)一定閾值，均進(jìn)行閃爍報(bào)警。

3.3 實(shí)時(shí)降雨預(yù)報(bào)

數(shù)據(jù)分為等值線圖、等值面圖、降水熱圖、面雨量分析列表4種表現(xiàn)方式，并對(duì)降水熱圖以動(dòng)畫(huà)方式演示，體現(xiàn)降水過(guò)程的動(dòng)態(tài)變化。

3.4 WEBGIS平臺(tái)顯示系統(tǒng)

用于在網(wǎng)頁(yè)上展示桃林口水庫(kù)流域電子地圖，并提供各種地圖操作和計(jì)算操作。整個(gè)網(wǎng)頁(yè)系統(tǒng)采用一個(gè)整體架構(gòu)模式，全部?jī)?nèi)容放置在菜單面板上，每打開(kāi)一個(gè)網(wǎng)頁(yè)，均會(huì)以標(biāo)簽的方式顯示在構(gòu)架之中，各個(gè)頁(yè)面可以進(jìn)行快速切換。

4 系統(tǒng)運(yùn)用

4.1 2018年“安比”“摩羯”臺(tái)風(fēng)暴雨洪水概況

2018年第10號(hào)臺(tái)風(fēng)“安比”和14號(hào)臺(tái)風(fēng)“摩羯”，以其異常罕見(jiàn)的路徑襲擊江浙滬地區(qū)之后，一路北上抵達(dá)山東，進(jìn)而影響北方地區(qū)。受臺(tái)風(fēng)“安比”“摩羯”及北部高氣壓影響，桃林口水庫(kù)所在青龍河流域先后發(fā)生了2018年“7·24”“8·13”兩場(chǎng)強(qiáng)降雨，累計(jì)平均降雨量200.2mm，最大降雨量站龍王廟，降雨量338mm。對(duì)應(yīng)兩場(chǎng)強(qiáng)降雨，桃林口水庫(kù)相繼出現(xiàn)2018年1號(hào)洪水、2號(hào)洪水。降雨過(guò)程及特征如表1。

表1 臺(tái)風(fēng)降雨特征

4.2 2018年降水預(yù)報(bào)與實(shí)際降水對(duì)比分析

2018年7月22日，臺(tái)風(fēng)“安比”在上海市崇明島登陸后，根據(jù)氣象部門(mén)預(yù)報(bào)，河北省將出現(xiàn)一次大范圍強(qiáng)降水過(guò)程，中東部降水量可達(dá)暴雨到大暴雨，局地可能超過(guò)200mm。

2018年7月22日18：00，河北省防指啟動(dòng)了防汛Ⅳ級(jí)應(yīng)急響應(yīng)，水庫(kù)利用衛(wèi)星云圖、雷達(dá)等多種網(wǎng)絡(luò)氣象信息手段，密切跟蹤臺(tái)風(fēng)路徑，研判降雨發(fā)展和落雨區(qū)域。2018年7月23日16：00，臺(tái)風(fēng)“安比”位于山東沂水縣境內(nèi)，中心氣壓990hPa，風(fēng)力8級(jí)，移速移向20km/h，北北西，水庫(kù)及時(shí)啟用精細(xì)化降雨預(yù)報(bào)系統(tǒng)，預(yù)報(bào)顯示7月23日20：00至25日20：00，水庫(kù)流域普面降暴雨到大暴雨，預(yù)報(bào)流域平均降雨量136.5mm，預(yù)報(bào)結(jié)果如圖1（a），水庫(kù)流域降雨實(shí)際監(jiān)測(cè)平均降雨量133.4mm，監(jiān)測(cè)結(jié)果如圖1（b），對(duì)比看出，預(yù)報(bào)和實(shí)際在降雨量和降雨分布上接近，流域平均降雨量預(yù)報(bào)誤差為2.30%。

圖1 “7·23”預(yù)報(bào)降雨量與實(shí)際降雨量

2018年8月12日，臺(tái)風(fēng)“摩羯”在浙江溫嶺沿岸登陸后，根據(jù)氣象部門(mén)預(yù)報(bào)，8月13日到15日，承德南部、唐山、秦皇島和石家莊、滄州以南地區(qū)有大雨到暴雨，局地有大暴雨。13日16：00，臺(tái)風(fēng)“摩羯”位于安徽省壽縣境內(nèi)，中心氣壓995hPa，風(fēng)力8級(jí)，移速移向30km/h，北西，水庫(kù)運(yùn)用精細(xì)化降雨預(yù)報(bào)系統(tǒng)，預(yù)報(bào)顯示8月13日20：00至15日20：00，流域普面降大到暴雨，南部邊緣及起河支流有大暴雨，平均降雨量達(dá)100～150mm，預(yù)報(bào)結(jié)果如圖2（a），預(yù)報(bào)流域平均降雨量70.6mm，監(jiān)測(cè)結(jié)果如圖2（b），實(shí)際監(jiān)測(cè)流域平均降雨量62.8mm，對(duì)比看出，預(yù)報(bào)和實(shí)際在降雨量和降雨分布上較接近，流域平均降雨量預(yù)報(bào)誤差為12.8%。

圖2 “8·13”預(yù)報(bào)降雨量與實(shí)際降雨量

4.3 為水庫(kù)科學(xué)調(diào)度提供重要依據(jù)

“安比”強(qiáng)降雨導(dǎo)致桃林口水庫(kù)與下游分水樞紐區(qū)間來(lái)水?dāng)y帶泥沙，致使分水樞紐水體濁度較大，對(duì)城市供水產(chǎn)生影響，7月26日18：00，根據(jù)預(yù)報(bào)后期降雨情況，水庫(kù)緊急啟動(dòng)一臺(tái)機(jī)組發(fā)電，流量29m3/s，提前預(yù)泄，起到對(duì)城市引水進(jìn)行清水置換處理、提前騰出庫(kù)容、 為洪水消峰創(chuàng)造條件及增加發(fā)電效益等綜合作用。而后，根據(jù)“摩羯”殘余環(huán)流的發(fā)展變化，水庫(kù)分別于8月14日14：00、8月15日11：30將泄水流量加大到180，260m3/s；8月20日11：00將泄水流量減少到130m3/s，8月21日11：00庫(kù)水位穩(wěn)定在137.60m，入庫(kù)流量與出庫(kù)流量保持平衡，水庫(kù)關(guān)閉泄洪閘門(mén)，轉(zhuǎn)入正常供水調(diào)度。

本次洪水調(diào)度調(diào)洪最高蓄水位138.26m，距控制目標(biāo)水位138.50m低24cm，最大限度地發(fā)揮了水庫(kù)消峰攔洪的作用，1號(hào)洪水全部攔蓄庫(kù)內(nèi)，削峰100%；2號(hào)洪水入庫(kù)洪峰流量455m3/s削減至104.5m3/s，削峰77.03%。

精準(zhǔn)的預(yù)測(cè)預(yù)報(bào)，為水庫(kù)科學(xué)研判提供了有利條件，贏得了寶貴的時(shí)間，確保了水庫(kù)上、下游人民群眾的生命安全，將洪水災(zāi)害損失降至最低。

4.4 為水庫(kù)洪水災(zāi)害預(yù)警提供技術(shù)支持

（1）桃林口水庫(kù)精細(xì)化降雨預(yù)報(bào)系統(tǒng)不僅可作為獨(dú)立系統(tǒng)使用，還可為洪水預(yù)報(bào)系統(tǒng)提供準(zhǔn)確的雨量輸入值，為水庫(kù)攔蓄與下泄提供及時(shí)準(zhǔn)確的科學(xué)依據(jù)。

（2）通過(guò)防汛調(diào)度中的實(shí)際應(yīng)用，逐步完善本流域小尺度降雨數(shù)值預(yù)報(bào)的實(shí)用系統(tǒng)，對(duì)流域內(nèi)不同單位面積的降雨進(jìn)行分布式的預(yù)報(bào)及監(jiān)測(cè)，從而滿足分布式水文模型建模應(yīng)用的數(shù)據(jù)要求。

5 結(jié)語(yǔ)

根據(jù)洪水預(yù)報(bào)結(jié)果，結(jié)合水庫(kù)蓄水、來(lái)水情況、退水預(yù)測(cè)及供水水質(zhì)等因素，提出下泄水量的調(diào)度建議，為河北省防指科學(xué)調(diào)度指揮提供科學(xué)依據(jù)。