(遼寧省遼陽水文局，遼寧 遼陽 111000)

遼河流域位于我國東北部，多年平均年降水量在520～860mm之間，主要集中在每年的6—9月份并且占全年的65%以上，流經(jīng)量在時(shí)空上的分布以西部區(qū)域?yàn)橹?；不同區(qū)域的地勢結(jié)構(gòu)存在一定差異，整體呈現(xiàn)出東南高、西北低的變化趨勢，高差處于216～840m范圍。流域范圍內(nèi)大小支流縱橫，降雨在時(shí)空分布上極不均衡，暴雨或強(qiáng)降雨較多，并且主要集中在汛期。據(jù)不完全統(tǒng)計(jì)，該流域每3年左右就會(huì)發(fā)生一次一般洪水，每隔7～8年的時(shí)間就會(huì)出現(xiàn)一次較大洪水。河道行洪排澇受上下游水位及主要支流狀況影響較為顯著，水位變化關(guān)系復(fù)雜。對(duì)河道匯流參數(shù)利用歷史水文資料進(jìn)行率定的傳統(tǒng)水文學(xué)法，往往難以客觀反映河道匯流變化規(guī)律和斷面形狀，在洪水預(yù)報(bào)中容易造成洪峰時(shí)間與洪峰量等要素的模擬與實(shí)際狀況存在較大的偏差，并對(duì)防汛工作產(chǎn)生不利影響[1-6]。

據(jù)此，本文對(duì)遼河流域中下游洪水演進(jìn)利用EasyRiv1D模型進(jìn)行模擬分析。該模型能夠綜合考慮河道斷面形成、水工構(gòu)筑物調(diào)度規(guī)則以及河流汊點(diǎn)等因素，系統(tǒng)、定量地認(rèn)識(shí)與把握河網(wǎng)洪水演進(jìn)過程，揭示河道匯流與河道斷面變化規(guī)律。本文成果可為政府部門防洪減災(zāi)提供可靠的科學(xué)性依據(jù)和技術(shù)理論數(shù)據(jù)支撐。

1 構(gòu)建EasyRiv1D洪水演進(jìn)模型

1.1 模型簡介

EasyRiv1D模型是由我國水利部門自主研發(fā)的一種能夠支持不同時(shí)空尺度快速模擬、DEM建模和參數(shù)自動(dòng)識(shí)別并具有自主知識(shí)產(chǎn)權(quán)的水力學(xué)模型。該模型是以河道水流不可壓縮與水流均質(zhì)為基本條件，并且假定河道斷面、坡降無較大變化幅度，然后利用圣維南方程組進(jìn)行模擬分析[7]。其方程組如下：

?A/?t+?Q/?x=q

(1)

(2)

以上式中 A——水流斷面的面積，m2；

Q——水流斷面的流量，m3/s；

q——基本單位河流長度上支流的流量，m3/s；

x——河流長度，m；

t——匯流時(shí)間，h；

Z——水位,m;

C——泄洪系數(shù);

R——水力半徑,m；

α——?jiǎng)恿康男Ｕ禂?shù)。

1.2 河道下滲參數(shù)耦合系統(tǒng)

遼河流域洪水演進(jìn)受河道下滲參數(shù)的影響作用較為明顯，地下水與地表水之間相互作用、相互影響，其相互關(guān)聯(lián)的變量是決定二者定耦合系統(tǒng)的關(guān)鍵，因此，地下水流可用如下方程表示：

(3)

(4)

(5)

以上式中Bgw——潛水含水層寬度，m；

K——滲透系數(shù);

M——含水層厚度,m;

S——給水度;

P——河道濕周,m；

b——河道防滲層厚度，m；

h——地下水水位，m；

Z,Z0——河道水位,河底高程或者河道防滲層底面的高程，m；

qg——交換量，正數(shù)為流入,負(fù)數(shù)為流出。

基于三級(jí)解法與加權(quán)離散對(duì)河網(wǎng)水流方程進(jìn)行求解，同時(shí)利用追趕分析法與邊權(quán)離散計(jì)算地下水流方程。其求解的主要流程如下：對(duì)初始側(cè)流項(xiàng)qg利用初始地下、地表水位和初始條件進(jìn)行求解，并假定初始值為河道下一刻側(cè)向流入量，即qg0=qg；然后對(duì)地下水位與河道內(nèi)水位進(jìn)行計(jì)算，并根據(jù)新的計(jì)算結(jié)果對(duì)比分析qg0、qgn，如果滿足收斂精度則繼續(xù)計(jì)算，否則重新調(diào)整qg0=qgn，收斂時(shí)，停止計(jì)算。

2 遼河流域河道洪水演進(jìn)模擬

2.1 EasyRiv1D模型概化

EasyRiv1D模型計(jì)算的重要基礎(chǔ)和前提條件是對(duì)其前處理，主要有如下內(nèi)容：河網(wǎng)拓?fù)浣Y(jié)構(gòu)輸入與徑流計(jì)算過程的處理；關(guān)鍵參數(shù)與河段關(guān)鍵斷面的輸入處理；模型運(yùn)行規(guī)則與水工建筑物尺寸的處理[8]，初始邊界條件的界定等。本文以遼河流域中下游為研究對(duì)象，整體連續(xù)模擬中下游河道洪水演進(jìn)流量過程和沿程水位，從而依據(jù)分析結(jié)果為洪水風(fēng)險(xiǎn)圖的編制和制定防洪減災(zāi)措施提供數(shù)據(jù)支持。模型上下邊界條件分別為福德店水文站出庫流量與遼源水文站水文過程。

2.2 模型參數(shù)

以河道實(shí)測斷面資料作為模型所需數(shù)據(jù)的主要來源，各斷面之間的距離平均為2～3km，本文共設(shè)計(jì)斷面260個(gè)，主要資料信息見表1。

表1 遼河中下游河道斷面信息

構(gòu)建模型的關(guān)鍵是對(duì)不同河段的滲透系數(shù)進(jìn)行科學(xué)合理的確定，遼河流域內(nèi)上下游與各河段之間存在著地下水埋藏深度、地層地質(zhì)巖石的特性等下墊面條件的很大差異，故此可以依照河道下滲其他項(xiàng)目檢測相關(guān)結(jié)論、河床巖性等資料，確定河道下滲在洪水演進(jìn)過程中的相關(guān)參數(shù)，結(jié)果見表2。

表2 河段滲透系數(shù)與地質(zhì)巖性

度量河道水流阻力的重要參數(shù)為河道糙率，該參數(shù)是影響河道泄洪能力預(yù)測、水庫回水末端位置和蓄水高度模擬，各種泥沙沖淤和水力計(jì)算等方面的關(guān)鍵參數(shù)。遼河中下游河段的河道水動(dòng)力初始模擬的糙率值見表2，并通過試算確定適合的糙率值。利用曼寧公式與謝才公式相結(jié)合計(jì)算糙率，公式如下：

(6)

式中Q——斷面流量，m3/s；

A——斷面面積，m2；

R——水力半徑，m；

J——水力梯度。

2.3 水力學(xué)模型預(yù)報(bào)分析

對(duì)遼河中下游河道15場洪水的徑流與洪峰進(jìn)行模擬，其中12個(gè)場次作為率定期樣本，3個(gè)場次作為驗(yàn)證期樣本。結(jié)果顯示，模擬相對(duì)誤差在各個(gè)水文站存在較大差異，其原因可能與流量實(shí)測資料較差相關(guān)。在大水時(shí)水位流量關(guān)系曲線為遼河水系的流量檢驗(yàn)主要依據(jù)，隨著人類活動(dòng)與經(jīng)濟(jì)社會(huì)對(duì)河道斷面影響的不斷增大，水文資料已不能準(zhǔn)確、客觀地反映河流斷面實(shí)際情況。并且，在汛期強(qiáng)降雨時(shí)遼河流域修建的橡膠壩，可出現(xiàn)人造洪峰的現(xiàn)象并使得河道驗(yàn)證期與率定期存在較大的誤差[9-10]。各站點(diǎn)模擬成果特征值見表3。

表3 各站點(diǎn)模擬成果特征值

表中H、T、R分別為渾河、太子河、饒陽河支流場洪，從驗(yàn)證期3個(gè)場次洪水特征值可以看出，實(shí)測值與H、T、R水文站的數(shù)值模擬具有良好的一致性，模擬結(jié)果基本能夠反映河流水文斷面流量變化特征，并且實(shí)測洪峰時(shí)間與模擬結(jié)果基本保持一致，EasyRiv1D模型具有良好的擬合效果，在我國東北區(qū)域的洪水預(yù)報(bào)中具有良好的實(shí)際應(yīng)用價(jià)值。另外，各水文站的確定性系數(shù)均高于0.6，由此表明在洪水預(yù)報(bào)過程中模型表現(xiàn)出較強(qiáng)的適用性與較高的精度。

2.4 定量預(yù)報(bào)的不確定性分析

為了避免對(duì)預(yù)報(bào)結(jié)果不確定性的影響，可利用不確定處理器HUP定量評(píng)估洪水預(yù)報(bào)結(jié)果[11-12]。該模型首先假定實(shí)測流量H0，利用Hn、Sn;(n=1,2,…,N)分別表示實(shí)測流量、模型確定性預(yù)測流量，N為預(yù)見期長度。假定hn、sn分別為Hn、Sn的取值，模型起始預(yù)測時(shí)間H0=h0。因此可利用貝葉斯公式計(jì)算概率密度函數(shù)φn，結(jié)合似然函數(shù)和先驗(yàn)密度函數(shù)進(jìn)行估計(jì)，計(jì)算公式如下：

(7)

確定性預(yù)報(bào)結(jié)果sn是模型耦合結(jié)果，利用式(7)可得到在hn、sn條件下預(yù)報(bào)量hn的概率分布φn(hn|sn,h0)，對(duì)其預(yù)報(bào)結(jié)果的不確定性和相似性加以分析和評(píng)估。本文以渾河典型水文站的5場次洪水為案例進(jìn)行了分析，并對(duì)其利用式(7)和求算方法獲取90%預(yù)報(bào)區(qū)間等數(shù)據(jù)，計(jì)算結(jié)果見表4。

表4 渾河洪水預(yù)報(bào)不確定性定量分析結(jié)果

由表4計(jì)算結(jié)果可知，預(yù)報(bào)區(qū)間為90%的離散度處于0.25～0.28區(qū)間，并且百分比均大于95%。通過上述計(jì)算和分析可知，耦合模型預(yù)測結(jié)果值的不確定性相對(duì)比較低，對(duì)于確定性預(yù)報(bào)結(jié)果Q50%的預(yù)報(bào)結(jié)果來說，洪峰流量預(yù)測結(jié)果的序列可靠性與精確性相當(dāng)高，具有一定的應(yīng)用價(jià)值。

3 結(jié) 論

本文基于遼河流域?qū)崪y水文資料和水量平衡原理構(gòu)建了一維水動(dòng)力學(xué)模型，預(yù)測分析了各代表性河段的徑流與洪峰誤差，以期為該流域洪水治理和防汛工作的開展提供一定參考，得出的主要結(jié)論如下：

a.模擬相對(duì)誤差在各個(gè)水文站存在較大差異，其原因可能與流量實(shí)測資料較差相關(guān)。隨著人類活動(dòng)與經(jīng)濟(jì)社會(huì)對(duì)河道斷面影響的不斷增大，水文資料已不能準(zhǔn)確、客觀地反映河流斷面實(shí)際情況，在汛期強(qiáng)降雨時(shí)遼河流域修建的橡膠壩，可出現(xiàn)人造洪峰的現(xiàn)象并使得河道驗(yàn)證期與率定期存在較大的誤差。

b.實(shí)測值與H、T、R水文站的數(shù)值模擬具有良好的一致性，模擬結(jié)果基本能夠反映河流水文實(shí)際狀況，具有良好的擬合效果，并且各水文站的確定性系數(shù)均高于0.6。EasyRiv1D模型在洪水預(yù)報(bào)過程中表現(xiàn)出較強(qiáng)的適用性與較高的精度。

c. EasyRiv1D模型考慮了河流下滲損失的影響作用，通過模型的預(yù)測可以為北方干旱半干旱地區(qū)的洪水預(yù)測的實(shí)際應(yīng)用提供理論依據(jù)和數(shù)據(jù)支撐，同時(shí)，為遼河區(qū)域水資源保護(hù)和開發(fā)利用提供一定參考和借鑒。