撫順渾河流域洪水分析及水位短期預(yù)測

楊 驍

（遼寧省撫順?biāo)木?遼寧 撫順 113006）

1 引言

我國在過去幾十年中,由于河岸泛濫以及導(dǎo)致人口過剩的城市化,洪水的脆弱性預(yù)計會上升。在過去的半個世紀(jì)里,國內(nèi)與洪水有關(guān)的死亡人數(shù)以及相關(guān)的經(jīng)濟(jì)損失急劇增加。這可能會在極端洪水事件發(fā)生后產(chǎn)生嚴(yán)重后果,人口和經(jīng)濟(jì)發(fā)展不斷增加,造成持續(xù)壓力[1-2]。此外,由于流域內(nèi)的測量站有限且管理不善導(dǎo)致對流域內(nèi)洪水防治不及時,洪水發(fā)生時造成嚴(yán)重?fù)p失[3-4]。

撫順渾河流域易發(fā)生洪水,在過去十年中,由于暴雨頻發(fā),該地區(qū)的洪水,盡管做出應(yīng)對策略,但洪水預(yù)報服務(wù)和洪水警報方案仍然不夠完善,這些方案旨在提前了解水體狀況,并在某種情況下提供不同的警報,以減輕洪水影響。本文以渾河有研究對象,利用撫順氣象站1974 年～2017 年的月降水觀測資料（來自中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)）。同時,從撫順市水文局收集了1964 年～2017 年期間渾河撫順?biāo)恼緶喓铀缓团潘康臄?shù)據(jù)。評估了集水區(qū)對降雨和洪水模式的響應(yīng),通過提前預(yù)報降雨情形預(yù)判洪水發(fā)生時間,從而提前做出策略,減輕洪水造成的損失。

2 洪水模式分析

對長期降雨分布的分析表明, 多年平均降雨量為776.1 mm, 月平均降水量最低為8 mm , 發(fā)生在1 月份, 最高為196.3 mm,發(fā)生在7 月份。雨季主要分布在6、7、8 月,月均降雨量均超過100 mm,月最大降雨量分別為223.5 mm、436.1 mm 和346.5 mm。降雨量高于平均值的年份是更容易發(fā)生洪水的年份,其中1995 年、2010 年、2013 年和2016 年降水量均高于1000 mm,與渾河下游洪水息息相關(guān)。

對撫順氣象站的降雨量和渾河撫順?biāo)恼镜乃唬ㄋ恼緮嗝婧哟驳撞康剿矫婢嚯x）進(jìn)行聯(lián)合分析有助于確定洪水模式。尤其是渾河流域經(jīng)歷了不同的降雨季節(jié),形成了不同洪水類型?？煞譃槿悾▓D1）：分別為單洪峰的洪水（一個長波）、多峰洪水、漫延洪水（一波無明顯峰值）。

圖1 洪水類別

表1 撫順氣象站長期（1974 年～2016 年）月降雨量分布

圖1（a）表示三個水文年（1968 年、1978 年和1994 年）的最高洪峰的洪水。研究顯示,從4 月開始到9 月中旬,出現(xiàn)了單次洪水峰值。從雨季開始,水位緩慢上升,形成洪峰流量。通常,這類洪水的水位最高。洪峰過后,洪水的水位緩慢下降,直到旱季開始圖1 （a）。

具有多個峰值的洪水如圖1（ b）所示,有多個波峰；通常,第一次波峰在3 月底開始,5 月底開始衰退。第二次波峰始于6 月,7 月開始衰退,而第三次波峰始于8 月,9 月結(jié)束。1965 水文年共出現(xiàn)三次峰值其中4 月觀測到第一次波峰的高水位為5.1 m,7 月中旬觀測到第二次波峰的最高水位為4.3 m,8 月中旬觀測到第三次波峰的最高水位為4.8 m；1984 年水文年出現(xiàn)多次峰值；5 月上旬觀測到兩次水位峰值,分別為4.9 m 和4.6 m,7 月上旬觀測到兩次水位峰值,分別為5.2 m 和4.8 m,8 月初觀測到一次水位峰值,為5.2 m；1991 水文年觀測到的水位峰值較低最大為8 月中旬的4 m。

圖1（c）表示水文年份（1966 年、1971 年和1972 年）的漫延洪水。通常,這類洪水的水位較高,在整個雨季有幾個中等峰值。峰值流量與其他峰值沒有顯著差異。這些洪水表明,水位從3 月開始上升,到9 月中旬結(jié)束。與其他水文年相比,1972 水文年的水位最高。這類洪水的最高洪峰出現(xiàn)在1972 水文年的7 月份,水位為4.1 m。

此后,對52 個水文年的洪水進(jìn)行了分析。還計算了平均值、最小/最早和最大/最晚水位及其日期（表2）。

表2 1965 年～2016 年渾河洪水參數(shù)匯總

表2顯示,在有記錄的52 個水文年中,洪水平均從3 月24 日開始,到9 月21 日結(jié)束。6 月29 日平均出現(xiàn)洪峰流量,最高水位為6.00 m。9 月21 日平均洪水結(jié)束,渾河洪水的平均持續(xù)時間約為180 天。雨季開始前的最低水位和峰值流量之間的水位上升約為3.10 m。

3 多元回歸法預(yù)測水位

為了預(yù)測渾河撫順?biāo)恼镜娜账?使用了渾河撫順?biāo)恼镜臄?shù)據(jù)和撫順氣象站的日降雨量總和。預(yù)測前兩天的水位（Hn-2）對應(yīng)于河床內(nèi)一定的積水水平。預(yù)測前一天的水位變化（ΔHn-1）對應(yīng)于特定集水區(qū)內(nèi)的當(dāng)前趨勢,作為對有無降雨的反應(yīng)。在預(yù)測前一天（Xn-1）使用降雨量,以比較收到的降雨量和水位上升的反應(yīng)。此類相關(guān)性的公式如下所示：

式中：Hn為渾河撫順?biāo)恼镜乃?cm；Hn-2為預(yù)測前兩天渾河撫順?biāo)恼镜乃?cm；ΔHn-1為預(yù)測前一天渾河撫順?biāo)恼镜乃蛔兓?cm；a0、a1、a2為通過普通最小二乘法計算的系數(shù)；X為撫順氣象站的日降雨量,mm；n為預(yù)測發(fā)布日期。

由于某些天的降雨量值等于“0”,因此對兩類方程進(jìn)行了多元回歸分析：

1）降雨量大于“0”的天數(shù)-水位上升的潛在天數(shù)。

2）對于降雨量等于“0”的天數(shù)-由于集水區(qū)內(nèi)的現(xiàn)有水量而具有靜態(tài)水流的天數(shù)。

3.1 降雨量大于“0”的水位預(yù)測

收到的多元回歸系數(shù)為0.87,表明所建立的回歸（2）適用于短期預(yù)測交付。渾河撫順?biāo)恼窘涤炅看笥凇?”天的水位預(yù)測公式如下：

降雨量等于“0”時的水位預(yù)測

多元回歸系數(shù)為0.90,表明所建立的回歸（3）適用于短期預(yù)測交付。因此,系數(shù)0.90 表明,在無降雨的天數(shù)內(nèi),所選變量與渾河撫順?biāo)恼镜乃恢g存在很強(qiáng)的關(guān)系。在降雨量小于“0”的天數(shù)內(nèi),渾河撫順?biāo)恼镜乃活A(yù)測公式如下：

3.2 多元回歸質(zhì)量評價

圖2顯示了2015 水文年的預(yù)測和觀測水位,并以最小流量天數(shù)進(jìn)行了相同的預(yù)測,通過方程2 和3 得出的多元回歸公式與觀測數(shù)據(jù)相比提供了良好的結(jié)果。建立的方程式2 和方程式3 的質(zhì)量評估基于允許誤差和成功預(yù)測的百分比。方程2 和3 的結(jié)果表明,91%和96%以上的預(yù)測是成功的。因此,建議在交付周期為1 天的渾河水位預(yù)測中使用式（2）和式（3）。

4 結(jié)論

渾河流域經(jīng)歷了不同的降雨季節(jié),下游形成了不同洪水類型：單峰洪水、多峰洪水和蔓延洪水。一般來說,洪水每年約有4 個峰值,但在某些年份可能達(dá)到6 個,而在其他年份,洪水可能只有1 個明顯的峰值。為了預(yù)測渾河撫順?biāo)恼镜娜账?使用了渾河撫順?biāo)恼镜臄?shù)據(jù)和撫順氣象站的降雨量數(shù)據(jù)。這些數(shù)據(jù)通過兩個多元回歸分別對有降雨和無降雨的天數(shù)進(jìn)行分析。對建立的降雨天數(shù)方程質(zhì)量的評估表明,預(yù)測的成功率達(dá)到91%～96%。為提前做好洪水防御措施提供了理論依據(jù)。