鄧 鵬李致家郭 元孫善磊

(1.南京信息工程大學(xué)水文氣象學(xué)院,江蘇南京 210044;2.河海大學(xué)水文水資源學(xué)院,江蘇南京 210098)

淮河上游地區(qū)降雨徑流分析及模擬

鄧 鵬1,李致家2,郭 元2,孫善磊1

(1.南京信息工程大學(xué)水文氣象學(xué)院,江蘇南京 210044;2.河海大學(xué)水文水資源學(xué)院,江蘇南京 210098)

為了掌握淮河流域的水情信息,對淮河干流王家壩以上流域的降雨徑流進行分析和模擬。采用流域內(nèi)1985—2007年的站點歷史資料對降雨徑流進行分析,并且應(yīng)用分布式SWAT模型,采用SCE-UA算法進行參數(shù)自動校準(zhǔn)并進行檢驗,建立降雨徑流月尺度模擬模型。分析發(fā)現(xiàn)流域年降水量50%集中在6—8月,并且淮南山區(qū)息縣及潢川子流域的月平均徑流和徑流系數(shù)都大于屬于淮北平原區(qū)的班臺子流域。息縣、潢川、班臺及王家壩月徑流的模擬結(jié)果確定性系數(shù)分別達到0.77、0.81、0.69和0.80?；春由嫌蔚貐^(qū)的降雨存在明顯的季節(jié)集中性,并且淮南山區(qū)和淮北平原區(qū)因地勢不同使產(chǎn)流特征不同。流域中水利工程設(shè)施的應(yīng)用對天然降雨徑流過程的改變產(chǎn)生很大影響,使班臺子流域的模擬結(jié)果差于其他子流域。研究區(qū)屬于徑流資料短缺地區(qū),分布式水文模型的建立充分反映了下墊面的已知信息,彌補了徑流資料短缺的不足之處。

淮河上游地區(qū);降雨徑流;水文模擬;分布式水文模型

淮河處在我國南北氣候過渡帶,汛期雨量豐沛但分布不均[1]。每當(dāng)汛期大暴雨時,淮河上游山地洪水迅速匯集,河道水位迅速抬高。加上淮河中游地勢平坦,河道比降小造成洪水排泄不暢,使淮河中上游地區(qū)防洪任務(wù)艱巨。為了防御洪水,淮河上游地區(qū)建立了一定數(shù)量的水庫、行蓄洪區(qū)和分洪道,這些工程使用頻繁,使淮河洪水的傳播過程更加復(fù)雜[2-5]。

淮河流域不利的氣候和地形組合,使流域內(nèi)洪澇災(zāi)害頻繁發(fā)生。最近10年來,淮河流域降水量年際變化很大,旱澇交替出現(xiàn)。2001年,淮河流域發(fā)生四季連續(xù)干旱的大旱;2003年、2007年淮河發(fā)生流域性大洪水;2005年淮河干流上游及部分支流出現(xiàn)較大洪水[6-7]。關(guān)于淮河流域降雨方面,大多數(shù)是關(guān)于極值降雨的變化及分布研究[6]。有關(guān)淮河徑流的研究,主要包括氣候變化情景下徑流的預(yù)測[8-10],以及淮河流域近年來大洪水事件的氣象成因[11-12]等,有關(guān)淮河流域歷史徑流的統(tǒng)計分析及變化趨勢研究較少見。

由于人為因素及下墊面條件不均的影響,淮河上游地區(qū)降雨徑流規(guī)律變得更加復(fù)雜。本文將通過歷史資料對淮河上游地區(qū)降雨徑流情況進行分析,并對降雨徑流過程進行分布式月時段水文模擬。月模型適合長時間的水文模擬,更注重水平衡過程,并且對資料要求不高,易于應(yīng)用[13]。該研究對淮河上游地區(qū)區(qū)域水文情勢的分析具有重要作用。

1 研究區(qū)域及數(shù)據(jù)

以淮河干流王家壩以上流域作為研究區(qū)域,應(yīng)用ArcGIS處理90 m分辨率SRTM高程數(shù)據(jù)[14],提取流域信息,如圖1所示。王家壩以上流域總面積30672 km2。流域上有多座水庫,其中直接對王家壩徑流產(chǎn)生影響的是圖1所示的6座水庫。其中宿鴨湖水庫包括桂莊和夏屯2個出流點,故圖1中有7個水庫出流點。以上游水文站息縣、潢川、班臺和水庫出流點劃分子流域如圖1所示。流域的南部地區(qū)多為山地,匯流迅速,為山溪性河流。流域北部為平原性河流,與淮河干流息縣以上和淮南山區(qū)顯著不同。息縣、潢川、班臺以上集水面積分別為10232 km2、2050 km2、11280 km2。

應(yīng)用1km分辨率2000年土地利用圖[15]對流域情況分析如表1所示。表中息縣、潢川、班臺指這些站點的上游集水面積流域,王家壩以上流域即是整個流域。從表中可以統(tǒng)計出息縣、潢川、班臺、王家壩以上流域的森林地覆蓋率分別為11.77%、19.31%、1.69%、5.90%,說明森林植被主要覆蓋于息縣和潢川以上的山區(qū)性地區(qū)。各流域的主要土地利用為農(nóng)田。采用10 km分辨率1999年土壤類型圖[16]對流域土壤情況進行分析,結(jié)果列于表1中。從表1可以看出息縣、潢川以上流域多為黏壤土,而班臺以上流域多為沙壤土。搜集了流域內(nèi)1985—2007年40個雨量站點的時段雨量資料,以及息縣、潢川、班臺、王家壩和各水庫出流點在汛期的徑流資料。搜集了流域內(nèi)1961—2007年寶豐、駐馬店、信陽、阜陽4個氣象站點資料,包括日最高、最低及平均氣溫,風(fēng)速,相對濕度及日照時數(shù)等。資料站點的分布如圖1所示。

2 模擬方法

采用SWAT模型對流域的降雨徑流過程進行月時段模擬。SWAT模型是由美國農(nóng)業(yè)部農(nóng)業(yè)研究中心開發(fā)的分布式流域水文模型,該模型可以在一個大型復(fù)雜的流域內(nèi),模擬或預(yù)測長期的降雨、土壤、土地利用和管理措施條件下對流域產(chǎn)流、產(chǎn)沙和化學(xué)污染物負荷的影響。SWAT模型中以水量平衡為基礎(chǔ),模擬包括降雨、徑流、下滲、蒸發(fā)蒸騰、基流、壤中流等水文過程。模型中采用SCS徑流曲線數(shù)法對流域地表徑流量進行模擬,采用動力貯水方法計算壤中流,降水中超過土壤蓄水容量的水分補給地下含水層。蒸散發(fā)量的計算采用Penman-Monteith公式計算蒸發(fā)能力,在此基礎(chǔ)上計算實際蒸散發(fā),包括植物截留層植物散發(fā)以及土壤水蒸發(fā)。模型針對水文響應(yīng)單元(HRU)計算匯流時間,再采用水力學(xué)公式分別計算坡面匯流和河道匯流,水庫匯流采用馬斯京根法[17]。

隨著計算機技術(shù)的迅猛發(fā)展,參數(shù)自動優(yōu)選方法得到了普遍應(yīng)用,大大加快了參數(shù)優(yōu)選的速度,增加了模擬結(jié)果的客觀性和可信度。水文模型參數(shù)優(yōu)選一直是水文預(yù)報的重要內(nèi)容,國內(nèi)外研究較多的水文模型參數(shù)優(yōu)選方法主要有單純形法、SCE-UA、遺傳算法等。SWAT模型中參數(shù)自動校準(zhǔn)采用SCE-UA數(shù)學(xué)算法[18]。參數(shù)校準(zhǔn)有一個原則:先校準(zhǔn)上游子流域的參數(shù),當(dāng)上游子流域的參數(shù)校準(zhǔn)好后便不再改變,接著校準(zhǔn)下游流域參數(shù)。因而這里先分別校準(zhǔn)上游息縣、潢川、班臺子流域的參數(shù),最后校準(zhǔn)這些子流域至王家壩區(qū)間的參數(shù)。另外研究區(qū)的徑流資料有短缺,這里先用模型對整個徑流過程進行連續(xù)模擬,再對有實測資料的徑流值和模擬值之間采用上述方法進行模型校準(zhǔn),以觀測值和模擬值差值的平方和最小作為自動校準(zhǔn)的目標(biāo)函數(shù)。采用相對誤差δ、復(fù)相關(guān)系數(shù)R2和確定性系數(shù)DC作為模型適用性評價指標(biāo)。

3 結(jié)果與討論

采用泰森多邊形插值法計算流域面雨量,得到1985—2007年各子流域及整個流域的月平均雨量變化過程如圖2所示。從圖2可以看出流域降水量從1月到7月迅速增加,然后從7月到12月迅速減小。流域年降水量的50%集中在6—8月。1985—2007年4—10月月平均徑流深變化過程如圖3所示。從圖2、圖3可以看出,潢川以上流域面雨量最大,徑流量也最大。潢川以上流域面積較小,且為山地地形流域,雨量較大的一個原因是山地對上升氣流的阻礙使雨量加大,而徑流量大的原因是山區(qū)地形使產(chǎn)流過程陡漲陡落,土壤對雨量的再分配作用比其他流域小。班臺流域徑流量最小,這是由于班臺以上流域地勢平緩,為平原區(qū)流域,洪水漲落都比較慢,土壤及河網(wǎng)對降水量的再分配作用比較明顯。計算得到息縣、潢川、班臺、王家壩以上流域4—10月的徑流系數(shù)分別為0.41、0.45、0.28、0.35,同樣說明息縣、潢川這樣的山地性流域產(chǎn)流效率要大于班臺流域及整個王家壩流域。各流域的降雨徑流年際變化線性趨勢如表2所示。其中班臺流域汛期徑流量的增加達到了0.05的信度檢驗。潢川流域雖然降雨量較其他流域大,但是其降雨量年變化反而呈減少趨勢。

降雨徑流的模擬計算中,根據(jù)敏感性分析結(jié)果,選定8個參數(shù)作為模型敏感參數(shù),即河道有效導(dǎo)水率CH-K2、SCS徑流曲線系數(shù)CN2、土壤蒸發(fā)補償系數(shù)ESCO、地下水再蒸發(fā)系數(shù)GW-REVAP、回歸流產(chǎn)生時需要淺層蓄水中的初始水深GWQMN、再蒸發(fā)發(fā)生時淺層蓄水層水深REVAPM、土壤可利用水量比SOL-AWC、土壤飽和導(dǎo)水率SOL-K。這8個參數(shù)的取值范圍如下:CH-K2=0～150 mm/h,CN2=35～98,ESCO=0～1,GW-REVAP=0.02～0.2,GWQMN=0～5000 mm,REVAPM=0～500 mm,SOL-AWC=0～1,SOL-K=0～2000 mm/h。利用息縣、潢川、班臺、王家壩4個水文站1985—1998年月徑流實測數(shù)據(jù),進行月尺度徑流校準(zhǔn)。校準(zhǔn)的模型參數(shù)值結(jié)果見表3。校準(zhǔn)完成后,使用1999—2007年的徑流觀測數(shù)據(jù)進行驗證。圖4分別是各個測站觀測徑流與模擬徑流的比較。表4為各個測站徑流校準(zhǔn)與驗證期間的評價系數(shù)值。

從表4中可以看出,息縣和潢川流域的模擬結(jié)果好于班臺流域,__________這是因為班臺以上流域有眾多水利工程設(shè)施,如板橋水庫、薄山水庫、宿鴨湖水庫、楊莊滯洪區(qū),以及很多行蓄洪區(qū)、閘門設(shè)施等。本模擬考慮到水庫放水對徑流的影響,而各種蓄、滯洪區(qū)設(shè)施的應(yīng)用對徑流的影響無法考慮,故模擬效果較差。比較而言,息縣和潢川流域為山區(qū)地形流域,除了已經(jīng)考慮的上游水庫,其他水利工程設(shè)施相對較少,降雨徑流關(guān)系比較有規(guī)律,故模擬結(jié)果較好。淮河上游地區(qū)人類活動對徑流過程的改變產(chǎn)生很大影響,特別是對于水利工程設(shè)施頻繁使用的流域,無法反映出天然降雨徑流的相關(guān)性過程。從各指標(biāo)中發(fā)現(xiàn)整個王家壩徑流模擬過程達到一般標(biāo)準(zhǔn)水平。

4 結(jié) 語

通過分析淮河上游地區(qū)的降雨徑流規(guī)律,發(fā)現(xiàn)淮河上游50%的年降雨量集中在6—8月,故產(chǎn)生季節(jié)性洪水的風(fēng)險較大?；茨仙絽^(qū)徑流系數(shù)大于淮北平原區(qū)徑流系數(shù),故淮南地區(qū)降雨轉(zhuǎn)換為河道洪水的幾率較高,防范突發(fā)性大洪水的任務(wù)更加艱巨。采用分布式水文模型進行月徑流模擬,因受水利工程設(shè)施使用的影響,可以得到一般標(biāo)準(zhǔn)水平的模擬結(jié)果。本次所選擇的研究區(qū)屬于徑流資料短缺地區(qū),所以建立網(wǎng)格分布式模型充分反映下墊面的已知信息,彌補徑流資料短缺的不足之處,并且由降雨等資料驅(qū)動模型得到連續(xù)的徑流模擬過程,然后根據(jù)實測徑流部分進行模型參數(shù)的率定與檢驗。該模擬計算可為無資料地區(qū)或資料短缺地區(qū)的水文模擬提供一定的參考。

[1]鄧鵬,李致家.近50年淮河流域年降水量變化分析[J].水電能源科學(xué),2013,31(4):14-17.(DENG Peng,LI Zhijia.Variation analysis of annual precipitation events in Huaihe River Basin in past 50 years[J].Water Resources and Power,2013, 31(4):14-17.(in Chinese))

[2]LI Zhijia,BAO Hongjun,XUE Cangsheng,et al.Real-time flood forecasting of Huai River with flood diversion and retarding areas[J].Water Science and Engineering,2008,1(2):10-24.

[3]包紅軍,趙琳娜,李致家.淮河具有行蓄洪區(qū)河系洪水預(yù)報水力學(xué)模型研究[J].湖泊科學(xué),2011,23(4):635-641.(BAO Hongjun,ZHAO Linna,LI Zhijia.Hydraulic model for flood forecast of river basin with flood diversion and flood retarding areas of Huaihe River[J].Journal of Lake Sciences,2011,23(4):635-641.(in Chinese))

[4]李致家,于莎莎,李巧玲,等.降雨-徑流關(guān)系的區(qū)域規(guī)律[J].河海大學(xué)學(xué)報:自然科學(xué)版,2012,40(6):597-604.(LI Zhijia,YU Shasha,LI Qiaoling,et al.Regional pattern of rainfall-runoff relationship[J].Journal of Hohai University:Natural Sciences,2012,40(6):597-604.(in Chinese))

[5]鄧鵬,李致家.3種水文模型在淮河息縣流域洪水模擬中的比較[J].河海大學(xué)學(xué)報:自然科學(xué)版,2013,41(5):377-382.(DENG Peng,LI Zhijia.Comparison of three hydrological models in flood simulation for Xixian Basin of Huaihe River[J].Journal of Hohai University:Natural Sciences,2013,41(5):377-382.(in Chinese))

[6]王嘉濤,梁樹獻,徐慧.淮河流域2000—2009年水勢分析[J].中國防汛抗旱,2011(1):21-24.(WANG Jiatao,LIANG Shuxian,XU Hui.Analysis of water resources in Huai River Basin between 2000 and 2009[J].China Flood&Drought Management,2011(1):21-24.(in Chinese))

[7]LI Zhijia,WANG Lili,BAO Hongjun,et al.Rainfall-runoff simulation and flood forecasting for Huaihe Basin[J].Water Science and Engineering,2008,1(3):24-35.

[8]佘敦先,夏軍,張永勇,等.近50年來淮河流域極端降水的時空變化及統(tǒng)計特征[J].地理學(xué)報,2011,66(9),1200-1210.(SHE Dunxian,XIA Jun,ZHANG Yongyong,et al.The trend analysis and statistical distribution of extreme rainfall events in the Huaihe River Basin in the past 50 years[J].Acta Geographica Sinica,2011,66(9):1200-1210.(in Chinese))

[9]高歌,陳德亮,徐影.未來氣候變化對淮河流域徑流的可能影響[J].應(yīng)用氣象學(xué)報,2008,19(6):741-748.(GAO Ge, CHEN Deliang,XU Ying.Impact of climate change on runoff in the Huaihe River Basin[J].Journal of Applied Meteorological Science,2008,19(6):741-748.(in Chinese))

[10]郝振純,鞠琴,王璐,等.氣候變化下淮河流域極端洪水情景預(yù)估[J].水科學(xué)進展,2011,22(5):605-614.(HAO Zhenchun,JU Qin,WANG Lu,et al.Prediction of extreme floods in Huaihe River Basin under climate change[J].Advances in Water Science,2011,22(5):605-614.(in Chinese))

[11]王黎娟,管兆勇,何金海.2003年淮河流域致洪暴雨的環(huán)流背景及其與大氣熱源的關(guān)系[J].氣象科學(xué),2008,28(1):1-7.(WANG Lijuan,GUAN Zhaoyong,HE Jinha.The circulation background of the extremely heavy rain causing severe floods in Huaihe River valley in 2003 and its relationships to the apparent heating[J].Scientia Meteorologica Sinica,2008,28(1):1-7.(in Chinese))

[12]矯梅燕,金榮花,齊丹.2007年淮河暴雨洪澇的氣象水文特征[J].應(yīng)用氣象學(xué)報,2008,19(3):257-264.(JIAO Meiyan, JIN Ronghua,QI Dan.Meteorological and hydrological characteristics of flood related torrential rain over Huaihe River Basin in 2007[J].Journal of Applied Meteorological Science,2008,19(3):257-264.(in Chinese))

[13]XU C Y,SINGH V P.A review on monthly water balance models for water resources investigations[J].Water Resources Management,1998(12):31-50.

[14]SRTM90.International Scientific Data Service Platform[EB/OL].[2012-09-19].http://datamirror.csdb.cn/dem/files/ys.jsp.

[15]UMD.1km Global land cover[EB/OL].[2010-05-04].http://www.geog.umd.edu/landcover/1km-map/meta-data.html.

[16]NASA.GLDAS global soils dataset of Reynolds[EB/OL].[2010-05-04].http://www.ngdc.noaa.gov.

[17]肖倉軍,周文斌,羅定貴,等.非點源污染模型:SWAT用戶應(yīng)用指南[M].北京:地質(zhì)出版社,2010:5-10.

[18]DUAN Q,GUPTA V K,SOROOSHIAN S.Effective and efficient global optimization for conceptual rainfall runoff models[J].Water Resources Research,1992,28(4):1015-1031.

Rainfall-runoff analysis and simulation for upstream area of Huaihe River

DENG Peng1,LI Zhijia2,GUO Yuan2,SUN Shanlei1
(1.College of Hydrometeorology,Nanjing University of Information Science and Technology,Nanjing 210044,China; 2.College of Hydrology and Water Resources,Hohai University,Nanjing 210098,China)

In order to obtain information on the water regime of the Huaihe River Basin,the rainfall and runoff in the catchment upstream of the Wangjiaba sub-basin,in the main stream of the Huaihe River,were analyzed and simulated.Historical data from 1985 to 2007 in the basin were used for rainfall-runoff analysis.The distributed SWAT model was used for monthly rainfall-runoff simulation,and the SCE-UA algorithm was used for automatic calibration of parameters.The results show that 50%of annual rainfall was concentrated in the period from June to August.The mean monthly runoff and runoff coefficients of the Xixian and Huangchuan sub-basins in the southern mountainous area were higher than those of the Bantai sub-basin in the northern plain area.The determination coefficients of the simulation results of monthly runoff in the Xixian,Huangchuan,Bantai,and Wangjiaba subbasins were 0.77,0.81,0.69,and 0.80,respectively.The rainfall in the upstream area of the Huaihe River areas shows a significantly seasonal centrality.The characteristics of runoff of the southern mountainous area and northern plain area are different due to terrain difference.The water conservancy facilities have significant influence on the change of rainfall-runoff processes,resulting in the simulated results of the Bantai sub-basin being inferior to those of other sub-basins.The study area is a runoff data-deficient area,and the distributed hydrological model adequately reflects the given information of the underlying surface,which can offset this deficiency.

upstream area of Huaihe River;rainfall-runoff;hydrological simulation;distributed hydrological model

P339

:A

:1000-1980(2014)05-0377-05

10.3876/j.issn.1000-1980.2014.05.001

2013-07 30

南京信息工程大學(xué)?？蒲袉淤M項目(20100407);江蘇高校優(yōu)勢學(xué)科建設(shè)工程項目(PAPD);江蘇省普通高校研究生科研創(chuàng)新計劃(CXZZ12-0241)

鄧鵬(1983—),男,江蘇興化人,講師,博士,主要從事水文模型和水文氣象研究。E-mail:dp@nuist.edu.cn