魏恒志，陳洋波，劉永強，董禮明，徐章耀，王幻宇

(1. 河南省白龜山水庫管理局，河南 平頂山 467031；2. 中山大學(xué)地理科學(xué)與規(guī)劃學(xué)院，廣州 510275)

白龜山水庫位于沙河上游，河南省平頂山市區(qū)西南部，自20世紀60年代建成以來，在防治沙河洪水中發(fā)揮了關(guān)鍵作用。20世紀70年代上游的昭平臺水庫建成后，2庫聯(lián)合運用，進一步提升了沙河流域的洪水防治能力，有效減輕了沙潁河及淮河的防洪壓力。白龜山水庫已成為沙河洪水防治的骨干性水庫[1]。圖1為沙河白龜山水庫以上流域簡圖。

昭平臺水庫下距白龜山水庫51 km，控制沙河流域面積1 430 km2，昭平臺水庫建成后，白龜山水庫實際控制流域面積為1 310 km2，該區(qū)域一般稱為昭-白區(qū)間流域，本文稱白龜山水庫控制流域，簡稱白龜山流域。白龜山流域水系發(fā)達，支流眾多，流域面積較大的支流有9條，見圖1。

白龜山流域位于河南省常遇暴雨中心，洪水發(fā)生頻繁，空間分布極不均勻，洪水預(yù)報的難度大。白龜山流域?qū)偕絽^(qū)性河流，下墊面條件復(fù)雜，流域內(nèi)人類活動劇烈，加大了洪水預(yù)報的難度。昭平臺水庫的建成運用，改變了白龜山流域的洪水形成規(guī)律，進一步加大了白龜山水庫入庫洪水預(yù)報的難度。采用常規(guī)的集總式洪水預(yù)報模型無法滿足白龜山水庫防洪調(diào)度對白龜山水庫入庫洪水預(yù)報的精度要求，影響了水庫防洪調(diào)度效益的充分發(fā)揮。

圖1 沙河流域白龜山水庫以上簡圖Fig.1 Sketch map of Baiguishan Reservoir catchment

分布式物理水文模型是流域洪水預(yù)報模型的最新發(fā)展，由于它將流域劃分成精細化的單元，可以充分考慮流域下墊面的異質(zhì)性及降雨在流域空間分布上的不均勻性，從而提高流域洪水預(yù)報的精度，被稱為新一代流域洪水預(yù)報模型。近年來精細化的流域下墊面特性數(shù)據(jù)的大量涌現(xiàn)，以及模型參數(shù)優(yōu)選方法的進展，使得分布式物理水文模型進入了實用階段。自世界上第1個分布式物理水文模型SHE模型[2]提出以來，國內(nèi)外已提出了1批分布式物理水文模型，代表性的如VIC模型[3]、WetSpa模型[4]、Vflo模型[5]、WEP模型[6]、LL模型[7]和流溪河模型[8，9]等。

本文采用分布式物理水文模型——流溪河模型，針對白龜山水庫入庫洪水開展了研究，提出了白龜山水庫入庫洪水預(yù)報模型，并采用PSO算法優(yōu)選了模型參數(shù)。對歷史洪水的模擬發(fā)現(xiàn)，該模型具有較好的入庫洪水預(yù)報精度，洪水模擬效果明顯優(yōu)于NAM模型的模擬效果，可應(yīng)用于白龜山水庫入庫洪水預(yù)報。

1 數(shù)據(jù)收集與流溪河模型構(gòu)建

1.1 流域特性數(shù)據(jù)挖掘

流溪河模型建模所需的流域物理特性數(shù)據(jù)包括DEM、土地利用類型和土壤類型。本文研究采用的DEM數(shù)據(jù)來自于美國航天飛機雷達地形測繪計劃公共數(shù)據(jù)庫的DEM數(shù)據(jù)[10]，空間分辨率為90 m×90 m，見圖2(a)。土地利用類型數(shù)據(jù)采自于美國地質(zhì)調(diào)查局(USGS)30″×30″全球土地覆蓋數(shù)據(jù)庫[11]，空間分辨率為1 000 m×1 000 m，經(jīng)過重采樣處理得到空間分辨率為90 m×90 m的白龜山流域土地利用類型數(shù)據(jù)，見圖2 (b)。土壤類型數(shù)據(jù)取自于國際糧農(nóng)組織(FAO) 于2008年發(fā)布的中國土壤分布數(shù)據(jù)(http:∥www.isric.org/)，空間分辨率為1 000 m×1 000 m，經(jīng)過重采樣獲得空間分辨率為90 m×90 m的白龜山流域土壤類型數(shù)據(jù)，見圖2 (c)。

圖2 白龜山流域物理特性數(shù)據(jù)Fig.2 The terrain data of Baiguishan Reservoir catchment

1.2 實測洪水過程資料

白龜山水庫水情自動測報系統(tǒng)現(xiàn)有12個自動雨量站，雨量站在流域上的空間位置見圖1。本文研究過程中，收集了11場白龜山水庫入庫洪水資料，12個雨量站的降雨量和入庫流量均以小時為時段進行了整理。對雨量站的實測降雨，通過泰森多邊形法進行空間插值，得到各網(wǎng)格單元上的面雨量。

2 流溪河模型構(gòu)建與參數(shù)優(yōu)選

2.1 流溪河模型及其參數(shù)優(yōu)選方法簡介

流溪河模型[8,9]是由陳洋波等提出的流域洪水預(yù)報分布式物理水文模型。流溪河模型采用DEM將流域從水平方向劃分成單元，包括邊坡單元、河道單元和水庫單元，每個單元有獨立的流域物理特性和降雨量；在單元流域上進行蒸散發(fā)量及產(chǎn)流量的計算，各單元上產(chǎn)生的徑流量通過匯流網(wǎng)絡(luò)進行逐單元匯流至流域出口。邊坡匯流采用運動波法進行計算，河道匯流采用擴散波法進行計算，水庫匯流采用平移法進行計算，壤中流匯流采用Capmell公式進行計算，地下徑流匯流則采用線性水庫法計算。流溪河模型每個單元上共有14個參數(shù)，分成4種類型，包括氣象類參數(shù)、地形類參數(shù)、土壤類參數(shù)和植被類參數(shù)。早期的流溪河模型參數(shù)優(yōu)選是一個半自動化的方法[12,13]，需要人工干預(yù)，過程繁復(fù)，工作量大，不易獲取全局最優(yōu)解。流溪河模型現(xiàn)在采用優(yōu)化算法進行參數(shù)自動優(yōu)選[14,15]，既提高了模型參數(shù)優(yōu)選的效率，也提高了模型的性能。

2.2 流溪河模型構(gòu)建

構(gòu)建流溪河模型就是以DEM為依據(jù)，對流域進行單元劃分，提取匯流網(wǎng)絡(luò)，對單元進行分類，在此基礎(chǔ)上，對河道斷面尺寸進行估算。采用空間分辨率為90 m×90 m的DEM對流域進行劃分，按照流溪河模型中的單元分類方法對單元進行分類，將單元分成了7 251個水庫單元、2 210個河道單元和359 991個邊坡單元。河道劃分為2級河網(wǎng)，參照Google Earth遙感影像，設(shè)置了22個河道結(jié)點，將河道分成了33個虛擬河段，并估算了各個虛擬河道的斷面寬度、側(cè)坡及底坡。單元分類及河道虛擬結(jié)點和虛擬河段的劃分結(jié)果見圖 3。

圖3 模型結(jié)構(gòu)示意圖Fig.3 Schematic diagram of model structure

2.3 模型初始參數(shù)推求

流溪河模型是一個分布式物理水文模型，每個單元上均采用不同的模型參數(shù)。流向和坡度為不可調(diào)參數(shù)，直接根據(jù)DEM計算確定?？烧{(diào)參數(shù)需先確定一個參數(shù)初值，再對參數(shù)進行優(yōu)選。其中，蒸發(fā)系數(shù)v是個非常不敏感的參數(shù)，本文采用推薦值0.7作為其初值。邊坡糙率與河道糙率為土地利用類型的可調(diào)參數(shù)，本文根據(jù)文獻[16]的推薦值確定初值。土壤類型可調(diào)參數(shù)包括土壤厚度Zs、土壤飽和含水率Csat、田間持水率Cfc、飽和水力傳導(dǎo)率Ks、土壤凋萎含水率Wl和土壤特性b共6個。土壤特性b本文采用流溪河模型的推薦值2.5作為初值，飽和含水率、田間持水率和飽和水力傳導(dǎo)率本文采用由Arya等人提出的土壤水力特性計算器[17]來計算確定初值。限于篇幅，這些值在本文未列出。

3 可調(diào)參數(shù)自動優(yōu)選

選擇2010071809號洪水，采用粒子群算法對白龜山流域流溪河模型可調(diào)參數(shù)進行自動優(yōu)選，圖4列出了參數(shù)優(yōu)選計算過程中的部分結(jié)果。

圖4(a)為參數(shù)優(yōu)選過程中適應(yīng)值(目標函數(shù)值)的變化過程，圖4(b)為參數(shù)的變化過程，圖4(c)為模擬的洪水過程。從圖4中可以看出，隨著尋優(yōu)進程的推進，模型參數(shù)值和目標函數(shù)值趨近于穩(wěn)定。表 1列出了參數(shù)優(yōu)選結(jié)果。

統(tǒng)計該場洪水模擬效果的評價指標，確定性系數(shù)為0.974，相關(guān)系數(shù)為0.994，水量平衡系數(shù)為0.93，過程相對誤差為30.2%，洪峰相對誤差為0.5%，洪水模擬效果優(yōu)良。

4 模型驗證

4.1 洪水模擬

采用上述優(yōu)選的模型參數(shù)，對收集的其余10場洪水進行了模擬，統(tǒng)計了4個評價指標，見圖5、表2。

圖4 參數(shù)優(yōu)選過程部分結(jié)果Fig.4 Evolution process of parameter optimization

符號KsnManZsbBsBwCsatCfcvCwSs參數(shù)值0.5711.4951.4991.4781.4110.5231.4711.0010.5530.5810.9450.504

從上述的模擬結(jié)果來看，本文建立的白龜山水庫入庫洪水預(yù)報流溪河模型對實測的白龜山水庫入庫洪水過程有較好的模擬效果，其中，確定性系數(shù)平均達到0.857，相關(guān)系數(shù)達0.957，洪峰誤差為10.8%。根據(jù)《水文情報預(yù)報規(guī)范》 ( SL250-2000)，場次洪水模擬預(yù)報的峰值流量誤差在實測值的20%以內(nèi)的為合格，則本論文中建立的白龜山水庫入庫洪水預(yù)報模型的場次洪水模擬預(yù)報僅有1場不合格，預(yù)報方案可用于白龜山水庫入庫洪水預(yù)報。

圖5 2種模型模擬的洪水過程對比Fig.5 Comparison of hydrological processes simulated by two models

4.2 與NAM模型的比較

為了與集總式模型的效果進行對比，以檢驗分布式模型對白龜山水庫入庫洪水預(yù)報是否具有優(yōu)越性，采用NAM模型開展對比研究。NAM模型[18-20]是Nielsen和Hansen于1973年提出的集總式概念性降雨徑流模型，在世界一些不同氣候類型地區(qū)得到不同程度的應(yīng)用[21-23]，目前在國內(nèi)流域洪水預(yù)報中應(yīng)用不多，應(yīng)用于水庫入庫洪水預(yù)報的則還沒有。白龜山水庫入庫洪水預(yù)報NAM模型采用與本文相同的洪水過程建立，使用其中的3場洪水進行參數(shù)率定，對其他8場洪水進行模擬，表2列出了2種模型對10場洪水(其中1995072402與2000070300這2場洪水為NAM模型用于參數(shù)率定的洪水)模擬效果的4個統(tǒng)計指標，圖5繪出了2種模型模擬的洪水過程圖。

從上述的模擬結(jié)果來看，本文建立的白龜山水庫入庫洪水預(yù)報流溪河模型對實測洪水的模擬效果明顯優(yōu)于NAM模型，流溪河模型模擬的洪水過程明顯好于NAM，這說明流溪河模型在應(yīng)用于白龜山水庫入庫洪水預(yù)報時的效果優(yōu)于NAM模型，也可以說明分布式模型在應(yīng)用于白龜山水庫入庫洪水預(yù)報時的效果優(yōu)于NAM模型。

5 結(jié) 語

本文采用分布式物理水文模型-流溪河模型建立了白龜山水庫入庫洪水預(yù)報模型，并對模型參數(shù)進行了自動優(yōu)選。對實測洪水過程進行模擬的結(jié)果表明，模型的確定性系數(shù)平均達到0.857，相關(guān)系數(shù)達0.957，洪峰誤差為10.8%，模擬效果明顯優(yōu)于NAM模型，可應(yīng)用于白龜山水庫入庫洪水預(yù)報模型。

表2 2種模型模擬的洪水過程統(tǒng)計指標對比Tab.2 Comparison of statistic index of the flood processes simulated by two models

□

[1] 袁晶瑄, 王本德，田 力. 白龜山水庫防洪預(yù)報調(diào)度方式研究及風(fēng)險分析[J]. 水力發(fā)電學(xué)報, 2010,(2):132-138.

[2] Abbott M B. An introduction to the European hydrologic system-system hydrologue European, 'SHE', 2: structure of a physically based [J]. Distributed Modeling System, Journal of Hydrology, 1986:87.

[3] Liang X, D P Lettenmaier, E F Wood, et al. A simple hydrologically based model of land surface water and energy fluxes for general circulation models [J]. Geophys. Res., 1994,99(D7):14 415-14 428.

[4] Wang Z, Batelaan O, De Smedt F. A distributed model for water and energy transfer between soil, plants and atmosphere (WetSpa) [J]. Phys. Chem. Earth, 1996,21:189-193.

[5] Vieux B E, J E Vieux. VfloTM: a real-time distributed hydrologic model[C]∥ Proceedings of the 2nd Federal Interagency Hydrologic Modeling Conference. July 28-August 1, 2002, Las Vegas,Nevada. Abstract and paper on CD-ROM, 2002.

[6] Jia Y, Ni G, Kawahara Y, et al. Development of WEP model and its application to an urban watershed [J]. Hydrological Processes, 2001,15(11):2 175- 2 194.

[7] 李 蘭,鐘名軍. 基于GIS 的LL-Ⅱ分布式降雨徑流模型的結(jié)構(gòu)[J]. 水電能源科學(xué),2003,21 (4).

[8] Yangbo Chen, Qiwei Ren, Fenghua Huang, et al. Liuxihe model and its modeling to river basin flood [J]. Journal of Hydrologic Engineering,2011,16(1):33-50.

[9] 陳洋波,任啟偉,徐會軍,等. 流溪河模型I:原理與方法[J]. 中山大學(xué)學(xué)報(自然科學(xué)版),2010,49(1):105-112.

[10] Sharma A, Tiwari K N. A comparative appraisal of hydrological behavior of SRTM DEM at catchment level[J]. J. Hydrol., 2014,519:1 394-1 404.

[11] Loveland T R, Reed B C, Brown J F, et al. Development of a global land cover characteristics database and IGBP discover from 1-km AVHRR data[J]. Int. J. Remote Sens., 2000,21:1 303-1 330.

[12] 陳洋波,任啟偉,徐會軍,等. 流溪河模型II:參數(shù)確定[J]. 中山大學(xué)學(xué)報(自然科學(xué)版),2010,49(2),105-112.

[13] 陳洋波,流溪河模型[M]. 北京：科學(xué)出版社,2009.

[14] 徐會軍,陳洋波,曾碧球，等.SCE-UA算法在流溪河模型參數(shù)優(yōu)選中的應(yīng)用[J].熱帶地理,2012,32(1):32-37.

[15] Yangbo Chen, Ji Li, Huijun Xu. Improving flood forecasting capability of physically based distributed hydrological model by parameter optimization[J]. Hydrology & Earth System Sciences, 2016,20:375-392.

[16] Sugawara M.Tank model[M]. Computer Models of Watershed Hydrology,1995:165-214.

[17] Arya L M, J F Paris. A physioempirical model to predict the soil moisture characteristic from particle-size distribution and bulk density data[J].Soil Sci. Soc. Am. J, 1981,45:1 023-1 030.

[18] Shamsudin S, Hashim N. Rainfall-runoff simulation using MIKE 11 NAM[M]. Journal of Civil Engineering,2002,15(2):1-13.

[19] Danish Hydraulic Institute(DHI).MIKE11: a modeling system for rivers and channels user-guide manual[R].DHI,2004.

[20] 林 波,劉琪璟,尚 鶴,等. MIKE 11/NAM模型在撓力河流域的應(yīng)用[J]. 北京林業(yè)大學(xué)學(xué)報,2014,(5):99-108.

[21] Loliyana, V D, Patel, P L. Lumped conceptual hydrological model for Puma river basin[J]. Sadhana- Academy Proceedings in Engineering Sciences,2015,40(8):2 411-2 428.

[22] Tawatchai Tingsanchali,Mahesh Raj Gautam. Application of tank, NAM, ARMA and neural network models to flood forecasting[J]. Hydrol. Process.,2000:1 414.

[23] Yee Cai Goh,Zurina Zainol,Mohd Zaki Mat Amin. Assessment of future water availability under the changing climate: case study of Klang River Basin, Malaysia[J]. International Journal of River Basin Management,2016:141.

[24] Sherman L K.Stream flow from rainfall by the Unit-Graph Method[R].Eng. News-Rec., 1932,108:501-505.

[25] Horton R E. The role of infiltration in the hydrologic cycle[J]. Transactions of the American Geophysical Union,1933,14:446-460.

[26] Bicknell B R, Imhoff J C, Kittle Jr J L, et al.Hydrological simulation rrogram-fortran(HSPF): uesr's manual for release 12.2[M]. US Environmental Protection Agency,Athens,Georgia,2005.

[27] Hydrologic Engineering Center.HEC-1, flood hydrograph package: uesr's manual[R]. Davis,California USA:US Army Corps of Engineers, Hydrologic Engineering Center,1981.

[28] Bergstiom S,Singh V P. The HBV model[M]. Computer Models of Watershed Hydrology,1995:443-476.

[29] 趙人俊.流域水文模型——新安江模型與陜北模型[M].水利出版社,1984.

[30] Freeze R A, Harlan R L. Blueprint for a physically-based, digitally simulated [J]. Hydrologic Response Model, Journal of Hydrology, 1969,(9):237-258.

[31] Kouwen N. WATFLOOD: a micro-computer based flood forecasting system based on real-time weather radar [J]. Canadian Water Resources Journal, 1988,13(1):62-77.

[32] Wigmosta M S, L W Vail, D P Lettenmaier. A distributed hydrology-vegetation model for complex terrain [J]. Water Resources Research, 1994,30(6):1 665-1 669.

[33] Julien P Y, Saghafian B, Ogden F L. Raster-based hydrologic modeling of spatially-varied surface runoff [J]. Water Resources Bulletin, 1995,31(3):523-536.

[34] Yang D, Herath S, Musiake K. Development of a geomorphologic properties extracted from DEMs for hydrologic modeling[J]. Annual Journal of Hydraulic Engineering JSCE, 1997,(47):49-65.