劉 玉 馬 宇

(1.廣東省水文局汕頭水文分局，廣東 汕頭 515041；2.黃河水利委員會黃河水資源中心，河南 鄭州 450003)

螺河發(fā)源于蓮花山脈的陸河縣南萬鎮(zhèn)境內(nèi)三神凸東坡，全長102km，流域面積1356km2。上游山區(qū)植被情況較好，大部分為林木荒草覆蓋；中、下游草木稀少，水土流失嚴重；下游地勢低洼，受潮水頂托，澇漬災(zāi)害較嚴重。因此，進行蕉坑站洪水預(yù)報研究對陸河縣、陸豐市防洪減災(zāi)工作具有積極意義。

目前國內(nèi)外開發(fā)的水文模型眾多，根據(jù)模型構(gòu)建基礎(chǔ)可以分為基于物理的模型、概念性模型及黑箱模型。我國自主研發(fā)的概念性新安江模型得到了國際上的廣泛認可，并且在我國濕潤及半濕潤地區(qū)得到了廣泛的推廣應(yīng)用[1]。螺河流域地處粵東沿海暴雨高區(qū)，平均年徑流量為19.02億m3，平均年降雨量為2140mm。多年的研究和實踐表明，此類氣候地區(qū)的降雨產(chǎn)流機制主要是蓄滿產(chǎn)流，江河常年水量豐沛，具備新安江流域模型和馬斯京根河道模型的應(yīng)用條件。本文基于廣東省建立的中小河流洪水預(yù)報預(yù)警平臺，應(yīng)用新安江模型編制了螺河蕉坑站洪水預(yù)報方案，并對方案精度的影響進行了分析研究。

1 方法與數(shù)據(jù)

1.1 新安江模型

新安江模型由河海大學(xué)趙人俊教授設(shè)計,是國內(nèi)第一個完整的流域概念性水文模型。20世紀80年代中期，借鑒山坡水文學(xué)概念及國內(nèi)外產(chǎn)匯流理論提出了三水源新安江模型。三水源新安江模型將凈雨劃分成地面徑流、壤中流以及地下徑流，其中地面徑流采用單位線進行匯流計算，壤中流和地下徑流經(jīng)過線性水庫的調(diào)蓄分別作為壤中流出流和地下水出流[2]。

1.2 研究區(qū)概況

蕉坑水文站于1955年3月9日由廣東省水利電力廳在廣東省陸豐縣(現(xiàn)陸豐市)河?xùn)|鎮(zhèn)蕉坑村設(shè)立，位于螺河中下游。1980年，因站網(wǎng)改造，蕉坑站基本水尺斷面遷往上游410m處，并改名為蕉坑(二)水文站。蕉坑(二)水文站的集水面積為1104km2，螺河干流的彎曲指數(shù)為2.0，河床比降為2.69‰。該站為收集螺河水文資料，結(jié)合螺河流域開發(fā)、水質(zhì)監(jiān)測及水情需要、水資源服務(wù)、水利工程服務(wù)而設(shè)，屬于二類精度站，主要任務(wù)為水位、流量、降雨量、蒸發(fā)、墑情、水質(zhì)監(jiān)測等。蕉坑(二)水文站流域概況見圖1。

圖1 蕉坑(二)站以上流域概況

1.3 數(shù)據(jù)選擇及處理

蕉坑(二)站以上閉合流域共有多個雨量站和水文站，經(jīng)對已有資料進行分析，河口、羅經(jīng)壩、萬全、丁洋、南告、各安、麥湖、青年、屯埔、下葫10 站的雨量數(shù)據(jù)較為可靠，平均控制面積110km2，且在流域內(nèi)的空間分布較為均勻[3]。采用泰森多邊形法得到各雨量站權(quán)重，進而求取流域平均雨量。雨量站權(quán)重見表1。

表1 蕉坑(二)站流域雨量站網(wǎng)權(quán)重

選用蕉坑(二)水文站1999—2014年平均逐月蒸發(fā)的均值作為模型計算資料，見表2。

表2 蕉坑(二)水文站1999—2014年平均逐月蒸發(fā)量

《水文情報預(yù)報規(guī)范》(GB/T 22482—2008)要求預(yù)報方案編制使用不少于10年的水文氣象資料，應(yīng)選取大、中、小洪水各代表性年份，且濕潤地區(qū)洪水場次不少于50次。因此選取1971—2005年資料作為率定期，考慮優(yōu)先選擇雨洪關(guān)系較好場次洪水，從中選取大、中、小共50場次洪水過程樣本供方案編制，選取2006—2007年2場洪水進行檢驗。

2 洪水預(yù)報方案構(gòu)建

2.1 新安江模型建模

本文采用三水源新安江模型構(gòu)建洪水預(yù)報方案，蒸散發(fā)采用三層蒸發(fā)模式計算，將土壤層劃分為上層、下層和深層，產(chǎn)流采用蓄滿產(chǎn)流模型，匯流采用滯后演算法。依托廣東省水文局組織建立的中小河流洪水預(yù)報系統(tǒng)平臺，根據(jù)流域特性確定模型的輸入、輸出，在此基礎(chǔ)上確定模型結(jié)構(gòu)，進行模型率定、校核，增加模型方法、作業(yè)預(yù)報和方案管理等[4]。

蕉坑(二)水文站以上為閉合流域，產(chǎn)流采用蓄滿產(chǎn)流模型，匯流采用三水源滯后演算法。方案計算步長為3h，預(yù)見期為24h，流域圈畫面積為1104km2。多年綜合水位流量關(guān)系見圖2，蕉坑(二)水文站模型預(yù)報結(jié)構(gòu)見圖3。

圖2 蕉坑(二)站多年綜合水位流量線

圖3 蕉坑(二)站模型預(yù)報結(jié)構(gòu)

2.2 模型參數(shù)率定及結(jié)果

采用自動優(yōu)選與人工調(diào)參相結(jié)合的方式率定參數(shù)，部分率定過程見圖4、圖5和表3，圖中藍色實線為實測流量，紅色實線為人工率定流量，率定參數(shù)見表4。

圖4 1973090220編號洪水率定過程

圖5 41990091105編號洪水率定過程

表3 蕉坑(二)站新安江模型方案評定

表4 蕉坑(二)站新安江模型參數(shù)

3 洪水預(yù)報方案分析

根據(jù)新安江模型參數(shù)模擬計算各次洪水樣本流量過程，統(tǒng)計洪峰流量、峰現(xiàn)時間、過程水量的相對誤差。由表4可知，蕉坑(二)水文站50場次洪水樣本中，洪峰平均相對誤差為22.11%，次洪產(chǎn)流量平均誤差為20.25%,平均過程確定性系數(shù)為0.76，根據(jù)《水文情報預(yù)報規(guī)范》(GB/T 22482—2008)，為乙等方案。

在近期的實際作業(yè)預(yù)報中，嘗試用此方案進行作業(yè)預(yù)報，選擇2018年1場洪水樣本對方案進行驗證。用該方案模擬20180830洪水，計算結(jié)果為8月31日12時出現(xiàn)洪峰流量3196m3/s，實測8月31日12時出現(xiàn)洪峰流量2810m3/s，預(yù)報峰現(xiàn)時間與實測時間一致，基本可以滿足作業(yè)預(yù)報使用。實時作業(yè)預(yù)報結(jié)果見圖6。

圖6 蕉坑(二)站20180830場次洪水預(yù)報結(jié)果

本次預(yù)報精度為乙等，經(jīng)過分析有以下幾點原因：一是由于選取歷史資料時間跨度較大，流域內(nèi)下墊面條件變化及河道斷面變化對水文過程都有較大的影響，城鎮(zhèn)化對流域下墊面的改變,使洪水特性較天然流域發(fā)生較大變化[5],因此通過一套模型參數(shù)來代表長時間多方面的流域特性，存在一定的不合理性。二是本站歷史洪水資料的洪峰流量大多相對較小，在相對誤差計算時對方案整體誤差評定影響較大。三是多年來站點上下游水庫、水閘的建設(shè)已破壞天然河道產(chǎn)匯流環(huán)境，部分歷史場次洪水資料受水利工程人為調(diào)控影響，出現(xiàn)了回水頂托[6]等情況，因此進一步降低了方案的確定性系數(shù)。在實際工作中，應(yīng)當(dāng)根據(jù)流域洪水特點建立不同預(yù)報方案，結(jié)合預(yù)報員長期經(jīng)驗，針對洪水實時特點，交互調(diào)整不同方案及參數(shù)進行作業(yè)預(yù)報，才能更好地提高洪水預(yù)報效果。

4 結(jié) 語

本文以螺河蕉坑(二)水文站為例，研究了洪水預(yù)報方案編制中新安江模型的應(yīng)用過程，并且在實際洪水作業(yè)預(yù)報中對本方案進行了進一步檢驗。洪水預(yù)報方案的建立其實是從眾多歷史資料中尋找預(yù)報流域的水文特性的工作，模型參數(shù)就是流域水文特性的概化，然后通過這些模型參數(shù)預(yù)報流域水文過程。但是，實際流域水文過程復(fù)雜多變，受諸多條件影響，傳統(tǒng)的由一組模型參數(shù)組成的預(yù)報方案很難準確地預(yù)報水文過程。諶潔等[7]指出近年水文學(xué)研究中均假定水循環(huán)的大氣過程、產(chǎn)流過程和匯流過程基本是漸變和緩變的，即服從水文一致性假定。但隨著全球氣候條件的不斷變化、人類活動疊加的影響，越來越多的研究成果表明，水文現(xiàn)象不確定性增加，一致性遭到破壞。建立在一致性假定基礎(chǔ)上的傳統(tǒng)徑流預(yù)報方法，預(yù)報精準度受到重大影響。

目前，在蕉坑(二)站實際工作中采用了多種預(yù)報方案，并結(jié)合歷史經(jīng)驗和雨洪條件進行作業(yè)預(yù)報?；诒敬畏治黾叭粘９ぷ鹘?jīng)驗，可以通過對歷史洪水的暴雨中心、降雨過程、上下游洪水特征等條件進行分類，來建立不同預(yù)報方案，組成預(yù)報方案集，在進行實時洪水預(yù)報中根據(jù)當(dāng)前雨洪特點選用不同預(yù)報方案集，得到集合預(yù)報結(jié)果。并且可以將新場次洪水進行分類，不斷更新方案集，進而降低流域特性變化、洪水資料對方案預(yù)報精度的影響。集合洪水預(yù)報，可作為蕉坑(二)水文站洪水預(yù)報發(fā)展研究的方向之一。