大洋河洪水預報技術研究

張 靜

（盤錦市水利服務中心,遼寧 盤錦 124010）

1 概述

1.1 流域概況

大洋河流域位于遼寧省東部，降水受西太平洋暖濕空氣及地形的影響，時空分布不均。降水量分布自西向東增加（多年平均降水量800 mm～1200 mm），年際變化大，年內集中在6 月～9 月（約占全年的76%～79%）。

全長230.2 km，是遼東半島地區最大的獨流入海河流。上游分為東、西兩支，西支為東洋河，發源于岫巖縣西北部偏嶺鄉的唐望山；東支為哨子河，發源于岫巖縣北部的古洞嶺。一般河寬300 m～400 m，河道平均坡度為1.1‰。水系發育良好，因河水流量大如汪洋之意而得名。

流域呈扇形形狀，流域面積6554 km2，地勢西高東低，多山峰丘陵。土壤多為山石土和淤土，有輕度侵蝕，植被覆蓋度40%～60%。大洋河流域屬扇形山區型流域，加之降水時空分布不均，洪水具有集流時間短，來勢迅猛的特點，洪峰暴漲暴落。

河流流經岫巖縣、鳳城市、東港市和大孤山經濟區等重要縣區，一旦發生洪澇災害必然帶來較大損失，因此在現有防洪工程的基礎上研究該流域洪水預報技術，對構建大洋河洪水防治體系具有重大意義。

1.2 工程設施

大洋河工程設施主要包括水文站網、水庫、防洪工程等。現有岫巖、沙里寨和龍王廟3 處水文站。基本情況見表1。

表1 大洋河水文站基本情況

干、支流上建有土門子大型水庫；黑山、羅圈背、刁家壩、廉家壩4 座中型水庫；前營、自由、江山頭、小甸子、臥龍趙、唐家隈子、躍進7 座小（一）型水庫。

防洪工程措施主要包括：在上游山區封山育林，涵養水源；在中游支流修建水庫，攔洪削峰，調蓄水量；上下游修建護岸、防洪堤、防潮堤，保護農田，已經初步形成以堤防為主的防洪體系。

1.3 研究內容

研究大洋河岫巖站、沙里寨站、龍王廟站洪水預報模型。依據大洋河流域預報方案和洪水預報系統要求，建立大洋河流域洪水預報水文模型，進行基于河、庫、站的流域聯合預報。

2 洪水預報

采用基于水文模型和網絡地理信息系統的結構，通過算法和邊界條件的人機交互，實現洪水預報。

2.1 水文模型構建

2.1.1 雙層指數模型

根據流域內各點包氣帶的蓄水容量不同的特點，繪成流域蓄水容量曲線。在雙層指數模型中，采用指數形式的蓄水容量曲線，見式（1）。

式中：F 為流域面積；F0為產流面積；a 為流域參數，反映流域蓄水容量分布特征；P 為降水量。

本模型所用飽和產流計算公式為：

式中：R 為徑流深；P 為降雨量；Pa為前期影響雨量；a 為流域參數。

匯流過程為流域上的地面凈雨、表層凈雨和地下凈雨轉化為流域出口斷面流量的過程。本模型的匯流采用單位線法。

2.1.2 三水源新安江模型

模型包括蒸散發計算、產流量計算和分水源計算三部分。本研究采用分布式三水源新安江模型，將大洋河流域劃分為多塊單元，通過產匯流計算得出單元出流過程；進行出口以下的河道洪水演算，獲得流域出流過程；最后將各單元的出流過程相加，得到流域總出流過程。

模型蒸散發計算采用三層蒸發模式，產流計算采用蓄滿產流模型。蓄水容量線型為拋物線型。

式中：f 為產流面積；F 為總流域面積；W 為已產流面積的蓄水容量；WM 為流域平均蓄水容量；b 為參數。

根據式（3）可推出流域最大張力水容量為：

流域張力水蓄水容量為：

流域產流量為：當P-K×EM+A

當P-K×EM+A≥WM 時（全流域蓄滿），

式中：R 為產流量；P 為降水量；K 為蒸散發折算系數；EM 為流域蒸發能力。

水源劃分通過自由水蓄水庫進行，匯流模型的三水源采用線性水庫法。

2.2 流域洪水模擬預報

在前述水文模型的基礎上進行流域水系河庫聯合洪水預報。

2.2.1 邊界條件

流域洪水模擬預報一般涉及降雨量及其時空分布、下墊面條件（流域蓄水量）、水庫防洪調度規程、河道行洪能力等諸多因素。為便于實際應用，本研究對這些邊界條件及因素進行概化，經過分析最終確定洪水模擬預報條件的邊界條件。

（1）流域蓄水量情況分為：飽和、半飽和。

（2）24 h 面降雨量設置分為：P=50 mm、100 mm、150 mm、200 mm、250 mm、300 mm，時程分布見表2。

表2 雨量時程分布(Δt=3 h)

（3）大型水庫按“汛限水位、死水位”與“汛限水位平均值”調洪；中小型水庫按“汛限水位、死水位”進行調洪。

2.2.2 技術方法

本研究主要使用經驗預報模型和數學預報模型。經驗預報模型涉及降雨徑流相關、上下游水位（流量）相關、河道經驗演算系數等方法。數學預報模型使用雙層指數模型和新安江模型。

3 應用實例

3.1 資料準備

降雨徑流資料來源于1966 年～2012 年系列水文年鑒。岫巖站流域平均雨量采用石灰窯、王家堡子和岫巖三站算術平均計算；沙里寨站流域平均雨量選取黃花甸子、谷家店、王家堡子、石灰窯、白家、岫巖、文家街、沙里寨8 處雨量站用算術平均法計算。各站蒸發計算均同岫巖站。徑流深使用實測洪水過程割去基流和前期退水過程，再逐時段累計。流域日蒸散發能力、蒸散發折算系數從《遼河流域實用洪水預報方案》中獲取。

3.2 水文模型構建

本研究使用的水文模型包括雙層指數模型和三水源新安江模型。按2.1 節構建流域水文模型，并進行參數率定。本文以大洋河沙里寨站為例，闡述各類模型產匯流模塊中參數的率定方法。

3.2.1 雙層指數模型

大洋河流域屬于濕潤地區，下墊面上層為枯枝落葉層，易蓄滿，與指數曲線不完全相符。鑒于此，將下墊面分為兩層，上層按蓄滿產流計算，下層按指數模型計算。雙層指數模型產流參數包括流域蓄水容量Im，上層蓄水容量Sup，上層系數b。經率定，參數選為Im=110 mm，Sup=30 mm，b=0.25。

依據2.1 節方法進行匯流計算，共分析了5 條單位線（如圖1 所示）。在進行洪水預報時可根據實時降雨的雨型、降雨強度和暴雨中心位置等情況綜合考慮，選擇適當的單位線進行計算。

圖1 大洋河流域匯流單位線

3.2.2 三水源新安江模型

大洋河流域最大蓄水量Wm采用110 mm，流域蓄水量W0采用單層計算的方法，計算公式為：

式中：K 為蒸散發折算系數；Pt為凈雨量。W0的計算日期是從每年的6 月1 日開始，以5 月上旬雨量的三分之一和五月中下旬雨量的三分之二之和作為W0的起始值。W0的上限以Wm控制。

流域的蓄水容量曲線采用b 次拋物線型。本流域b 值采用0.3。

根據上述參數計算并點繪P-W0-R 降水徑流相關圖，模型關系線與實際點據配合較好。用1966 年～2012 年中16 次洪水資料進行校核，徑流深方案合格率為81.2%。

匯流采用三水源滯后演算模型。坡地匯流階段，地面徑流直接進入河網；壤中流經過水庫消退作用入流河網；地下徑流經地下蓄水庫的消退入流河網。各種水源的河網匯流特性均主要受制于河槽水力學條件，三者之和經滯后演算法匯流至單元出口斷面。河道匯流采用分段馬斯京根連續演算法，將各單元出口流量演算至流域出口，再疊加。

三水源新安江模型率定后的產匯流參數見表3（以沙里寨站為例）。

表3 大洋河沙里寨站三水源新安江模型產匯流參數

3.3 流域洪水預報

按2.2 節進行“以流域干流為主線、以水庫和河道水文站為控制節點”的流域水系河庫聯合洪水模擬預報。以土門子水庫（大型）和大洋河沙里寨水文站為例，模擬預報成果見表4、表5。

表4 土門子水庫洪水預報成果

表5 大洋河沙里寨站洪水模擬預報成果

4 結論

本文使用雙層指數模型、三水源新安江模型構建大洋河各站洪水預報方案。依據大洋河流域預報方案和洪水預報系統要求，以大洋河干流為主、以干、支流水文站點為節點，建立大洋河流域洪水預報水文模型，為大洋河流域的河、庫、站聯合預報、全流域水系河庫洪水模擬預報提供參考。未來可據此進一步進行實時預報與洪水預警相關研究，最終構建洪水預報實時系統。