趙考生

（天津市水文水資源勘測(cè)管理中心，天津 300061）

網(wǎng)格新安江模型在于橋水庫(kù)洪水預(yù)報(bào)中的應(yīng)用

趙考生

（天津市水文水資源勘測(cè)管理中心，天津 300061）

根據(jù)于橋水庫(kù)流域的地理概況和水文氣象特性，利用網(wǎng)格新安江模型GRID-XAJ對(duì)流域歷史資料進(jìn)行降雨徑流模擬及參數(shù)率定，研制洪水預(yù)報(bào)方案。從模擬結(jié)果上看，網(wǎng)格新安江模型能有效反映流域下墊面變化及降水空間分布的離散對(duì)降水徑流的影響。

于橋水庫(kù)；洪水預(yù)報(bào)；網(wǎng)格新安江模型；參數(shù)率定

根據(jù)于橋水庫(kù)流域的地理概況和水文氣象特性，考慮人類活動(dòng)對(duì)下墊面的影響，應(yīng)用ArcView對(duì)于橋流域的數(shù)字高程模型DEM進(jìn)行識(shí)別和相應(yīng)的處理，選擇網(wǎng)格新安江模型GRID-XAJ對(duì)流域歷史資料進(jìn)行降雨徑流模擬及參數(shù)率定，模擬各子流域降雨、洪水過(guò)程，預(yù)報(bào)全流域暴雨之后的入庫(kù)水量、水位、流量及峰現(xiàn)時(shí)間等，建立洪水預(yù)報(bào)方案，提高洪水預(yù)報(bào)精度。

1 水文資料選取和子流域劃分

水文資料主要包括降水量、徑流量、水庫(kù)水位及流域地下水位等，資料年限為1960～2012年。本研究采用于橋水庫(kù)的庫(kù)面蒸發(fā)資料作為于橋水庫(kù)流域的蒸發(fā)值處理，以流域的DEM數(shù)據(jù)為基礎(chǔ)，根據(jù)生成的集水面積和水系流向，客觀充分利用地形資料，并結(jié)合人工經(jīng)驗(yàn)、流域內(nèi)雨量站、水文站分布情況等劃分自然子流域。

于橋水庫(kù)流域的控制面積2060km2，分為前毛莊、水平口、淋河橋、水庫(kù)區(qū)間4個(gè)一級(jí)子流域，其中每個(gè)一級(jí)子流域再用自然子流域方法來(lái)劃分二級(jí)子流域，劃分原則是盡可能利用水文站作為子流域的控制斷面，根據(jù)流域內(nèi)雨量站和水文站的布設(shè)情況及自然流域的邊界，用Arcgis軟件把流域（包括于橋水庫(kù)庫(kù)區(qū)水域）劃分為17個(gè)二級(jí)子流域。

2 模型參數(shù)

2.1 年降水徑流相關(guān)關(guān)系

采用1960～2012年各子流域內(nèi)平均降水量、實(shí)測(cè)徑流量（或反推入庫(kù)流量）、天然徑流量和用水量、流域面積等數(shù)據(jù)系列，做降雨徑流相關(guān)圖。其中天然徑流量為實(shí)測(cè)徑流量與用水量的總和，用水量為地表水的消耗量，該數(shù)據(jù)通過(guò)實(shí)際調(diào)查和統(tǒng)計(jì)分析得到。分別對(duì)黎河前毛莊以上、沙河水平口以上、淋河龍門口以上、庫(kù)區(qū)間4個(gè)一級(jí)子流域，自1960～2012年降雨徑流資料進(jìn)行相關(guān)法分析，得到降雨徑流關(guān)系。經(jīng)統(tǒng)計(jì)分析，1960～1980年、1981～1999年、2000～2012年這3個(gè)時(shí)段的降雨徑流關(guān)系線存在明顯偏差。不同時(shí)期、相同降雨產(chǎn)生的徑流量明顯變化，說(shuō)明不同時(shí)段內(nèi)下墊面條件發(fā)生了變化，本流域內(nèi)對(duì)徑流有明顯影響的人類活動(dòng)，主要是水土流失、城市發(fā)展和地下水過(guò)度超采引起的地面下沉。地下漏斗的形成，使包氣帶增大，徑流量減小。

本研究中，前毛莊子流域分為4個(gè)時(shí)期進(jìn)行模型參數(shù)率定，分別為1960～1980、1981～1984、1985～1999、2000年以后；其他3個(gè)子流域分為3個(gè)時(shí)期進(jìn)行模型參數(shù)的率定，分別為1960～1980，1980～1999，2000年以后。由資料的完整度和預(yù)見期的要求，各子流域均選用1h作為次洪水的預(yù)報(bào)時(shí)段。

2.2 張力水容量Wm與地下水埋深的關(guān)系

選取于橋水庫(kù)流域內(nèi)泉水頭、大柳樹、南新城、新店子4處地下水位觀測(cè)站，分別統(tǒng)計(jì)其歷年的年平均埋深、年最大埋深、年最小埋深的3種特征值。對(duì)系列分析發(fā)現(xiàn)：4處地下水位站的年平均埋深、年最大埋深、年最小埋深3個(gè)統(tǒng)計(jì)特征值全部呈增大趨勢(shì)。從各站地下水埋深的多年變化情況看，地下水位的下降，增加了包氣帶厚度，提高了田間持水能力，從而導(dǎo)致匯出徑流量的減少。

Wm表示流域張力水蓄水容量，計(jì)算公式：

式中 La為包氣帶厚度（mm）；均可以根據(jù)流域土壤類型通過(guò)查土壤參數(shù)統(tǒng)計(jì)表獲取，因此只要知道每個(gè)柵格單元的La，即可獲得Wm在流域的空間分布。

在自然界中，影響包氣帶厚度的因素較多，很難進(jìn)行直接推求。La可通過(guò)與地形指數(shù)及土壤類型對(duì)應(yīng)的土壤水分常數(shù)進(jìn)行估算，可假定地形指數(shù)大的地方包氣帶較薄而地形指數(shù)小的地方包氣帶較厚。這與實(shí)際情況也基本相符，一般而言，流域內(nèi)地形指數(shù)大的地方大多位于河道附近，而這些區(qū)域的地下水埋深較淺，包氣帶相對(duì)較薄；相反，地形指數(shù)小的地方基本位于流域的上游山坡，遠(yuǎn)離河道，包氣帶相對(duì)較厚。因此，可以假設(shè)流域上地形指數(shù)最大的柵格單元對(duì)應(yīng)的張力水蓄水容量最小，而地形指數(shù)最小的柵格單元對(duì)應(yīng)的張力水蓄水容量最大。

同時(shí)Wm表示土層最大可能缺水量，與包氣帶厚度有關(guān)。地下水埋深大，包氣帶厚度大，蓄滿包氣帶所需水量增大，即Wm加大。反之亦然。將地形指數(shù)與地下水埋深相結(jié)合，改進(jìn)式（1）的Wm計(jì)算方法，可得：

式中 Sup表示流域最大缺水量；TO表示地形指數(shù)；TOmax表示最大地形指數(shù)值，TOPO表示地形指數(shù)閾值；VE表示地下水埋深，VEmin指地下水最小埋深值；VERT表示地下水埋深閾值，由上式可知，Wm值隨著地下水埋深的增加而增加。

分別選取泉水頭、大柳樹站代表前毛莊子流域，新店子站代表水平口、水庫(kù)區(qū)間子流域，南新城代表淋河橋子流域。利用代表站地下水埋深與張力水容量進(jìn)行定量分析，按照降雨徑流相關(guān)年代的劃分，隨著年代的推移，地下水埋深增加，相應(yīng)的張力水容量也在增大。

3 模型參數(shù)率定及歷史洪水模擬

采用網(wǎng)格新安江模型對(duì)各子流域進(jìn)行歷史洪水模擬計(jì)算，根據(jù)于橋水庫(kù)流域年降雨徑流關(guān)系分析和張力水容量Wm與地下水埋深相關(guān)分析，各子流域不同時(shí)期取不同參數(shù)值。網(wǎng)格新安江模型各子流域不同時(shí)期相同參數(shù)如表1。

表1 網(wǎng)格新安江模型各子流域不同時(shí)期相同參數(shù)

3.1 前毛莊子流域

前毛莊子流域采用4個(gè)時(shí)期內(nèi)的7場(chǎng)典型洪水資料系列進(jìn)行洪水模擬。模擬的洪水徑流量相對(duì)誤差合格率為71%，洪峰合格率為100%，確定性系數(shù)為0.83，模擬結(jié)果良好。網(wǎng)格新安江分布式模型采用三水源劃分徑流，雖然洪量模擬有的場(chǎng)次偏大，但洪峰模擬值較好。

自20世紀(jì)80年代以來(lái)，隨著地下水埋深不斷下降，張力水容量、自由水容量不斷增大，產(chǎn)生這種現(xiàn)象的主要原因：①社會(huì)經(jīng)濟(jì)發(fā)展，人口增加，生活水平提高等，使用水量大大增加；②“引灤入津”工程的建立，使下墊面發(fā)生了變化，同時(shí)流域內(nèi)近年來(lái)多處修建水庫(kù)塘壩，也攔蓄了一部分水量，致使模擬的徑流量大于實(shí)測(cè)值。

3.2 水平口子流域

水平口子流域采用3個(gè)時(shí)期內(nèi)的12場(chǎng)典型洪水資料系列進(jìn)行洪水模擬。模擬的洪水徑流量相對(duì)誤差合格率為87.5%，洪峰合格率為100%，確定性系數(shù)為0.84。流域內(nèi)水庫(kù)、塘壩的建成使用，使原來(lái)的產(chǎn)匯流發(fā)生了變化，其中般若院水庫(kù)和上關(guān)水庫(kù)以上來(lái)水由這2座中型水庫(kù)控制，水庫(kù)塘壩在汛期末蓄水造成預(yù)報(bào)徑流值偏大，同時(shí)受人類活動(dòng)，社會(huì)經(jīng)濟(jì)發(fā)展的影響，流域內(nèi)用水量也大大增加，因此在模型參數(shù)率定時(shí)，為了達(dá)到較好的模擬效果，隨著年達(dá)增加，逐步加大參數(shù)Wm和SM的值。

3.3 淋河橋子流域

淋河橋子流域位于于橋水庫(kù)流域的中下游，2006年3月龍門口水庫(kù)建成，2012年5月原龍門口水文站下移設(shè)置為淋河橋站，控制面積由原來(lái)的126km2增加到210km2。龍門口以上區(qū)間的降雨徑流首先匯入龍門口水庫(kù)，經(jīng)調(diào)蓄后下泄流量經(jīng)淋河橋站匯入于橋水庫(kù)。本次參數(shù)率定以2012年以后實(shí)際情況為準(zhǔn)，選取2場(chǎng)洪水。模擬的洪水徑流量相對(duì)誤差合格率為100%，洪峰合格率為100%，確定性系數(shù)為0.75，模擬結(jié)果良好。

3.4 水庫(kù)區(qū)間子流域

水庫(kù)區(qū)間子流域采用3個(gè)時(shí)期內(nèi)的11場(chǎng)典型洪水資料系列進(jìn)行洪水模擬。模擬次洪水徑流量相對(duì)誤差合格率為73%，洪峰合格率為73%，確定性系數(shù)為0.52，洪峰、洪量模擬結(jié)果良好。由于于橋水庫(kù)沒有入庫(kù)水文控制站，實(shí)測(cè)入庫(kù)流量過(guò)程是利用庫(kù)水位變化反推而來(lái)的時(shí)段概化過(guò)程，而預(yù)報(bào)入庫(kù)洪峰是瞬時(shí)流量過(guò)程，故在確定性系數(shù)上與其他3個(gè)有瞬時(shí)實(shí)測(cè)流量過(guò)程的子流域略微偏低。

由模型特征值統(tǒng)計(jì)結(jié)果可以看出，網(wǎng)格新安江模型在洪量、洪峰與確定性系數(shù)的精度與合格率達(dá)到了較好效果，且網(wǎng)格新安江模型具有三水源的產(chǎn)匯流結(jié)構(gòu)，對(duì)于缺乏資料的地區(qū)也有較好的模擬精度，適合于橋水庫(kù)實(shí)時(shí)洪水預(yù)報(bào)作業(yè)。

4 洪水預(yù)報(bào)方案建立

將前毛莊、水平口以及淋河橋站（原龍門口站）的實(shí)測(cè)流量資料作為于橋水庫(kù)區(qū)間的入流，對(duì)于橋水庫(kù)庫(kù)區(qū)間進(jìn)行網(wǎng)格新安江模型模擬，產(chǎn)流匯流采用網(wǎng)格新安江分布式模型計(jì)算，選取11場(chǎng)典型洪水，計(jì)算出的匯流量經(jīng)過(guò)水庫(kù)庫(kù)面產(chǎn)流，最終得到預(yù)報(bào)流量；將其與于橋水庫(kù)計(jì)算得到的反推入庫(kù)流量進(jìn)行對(duì)比，率定分析參數(shù)，預(yù)報(bào)方案結(jié)構(gòu)如圖1，模擬特征值結(jié)果如表2。

表2 流域網(wǎng)格新安江模型模擬成果統(tǒng)計(jì)

圖1 預(yù)報(bào)方案結(jié)構(gòu)

由表2可以看出，本方案的徑流深合格率、洪峰合格率相等，為73%，且模型的確定性系數(shù)為0.51。實(shí)際預(yù)報(bào)工作中，可以選擇各子流域2000年以后率定出的模型參數(shù)進(jìn)行預(yù)報(bào)作業(yè)。

5 結(jié)語(yǔ)

通過(guò)對(duì)各子流域采用網(wǎng)格新安江模型（GRIDXAJ），進(jìn)行歷史降雨徑流過(guò)程模擬和模型參數(shù)率定。從模擬結(jié)果上看，模型效果較好，并在此基礎(chǔ)上提出了于橋水庫(kù)全流域入庫(kù)洪水過(guò)程的預(yù)報(bào)方案。

［1］魏兆珍.海河流域下墊面要素變化及其對(duì)洪水的影響研究［D］.天津：天津大學(xué)，2013.

［2］劉文具.流域下墊面變化對(duì)流域洪水的影響分析［J］.地下水，2013，35（4）：142-143.

［3］芮孝芳，黃國(guó)如.分布式水文模型的現(xiàn)狀與未來(lái)［J］.水利水電科技進(jìn)展，2004，24（2）：55-58.

［4］李致家，張珂，姚成.基于 GIS的 DEM和分布式水文模型的應(yīng)用比較［J］.水利學(xué)報(bào)，2006（8）：1022-1028.

［5］姚成.基于柵格的分布式新安江模型構(gòu)建與分析［D］.南京：河海大學(xué)，2007.

［6］李致家，胡偉升，丁杰，等.基于物理基礎(chǔ)與基于柵格的分布式水文模型研究［J］.水力發(fā)電學(xué)報(bào)，2012（2）：5-13，32.

［7］謝永玉.淮北平原降雨入滲補(bǔ)給系數(shù)隨地下水埋深變化特征［J］.地下水，2012（1）：9-11.

［8］李世峰，馬小雷，付丹平.地下水埋深的影響因素分析及模型研究［J］.地下水，2014（3）：21-23.

［9］孫仕軍，丁躍元，馬樹文.地下水埋深較大條件下井灌區(qū)土壤水分動(dòng)態(tài)變化特征［J］.農(nóng)業(yè)工程學(xué)報(bào)，2003，19（2）：70-74.

［10］陳建峰.降雨入滲補(bǔ)給規(guī)律的分析研究［J］.地下水，2010，32（2）：30-31.

［11］夏自強(qiáng)，李瓊芳.土壤水資源的變化和補(bǔ)給特征研究［J］.水文，2001，21（5）：1-5.

［12］Dawes W.Modelling the effects of climate and land cover change on roundwater recharge insouth-west Western Australia［J］.Hydrology and Earth System Sciences，2012，16（8）：2709-2722.

Application of GRID-XAJ model in Yuqiao Reservoir flood forecasting

ZHAO Kao-sheng
（Tianjin Hydrology and Water Resource and Survey Management Center，Tianjin 300061，China）

Based on the geography and hydrological and meteorological characteristics of Yuqiao reservoir watershed，the rainwater runoff simulation and parameter rate were simulated by GRID-XAJ model，and the flood forecasting scheme was developed.The results showed that the model can effectively reflect the influence of the variation of the underlying surface and the spatial distribution of the precipitation on the precipitation runoff.

Yuqiao Reservoir；hydrological forecasting；GRID-XAJ model；parameter calibration

TV122

B

1672-9900（2017）05-0030-04

2017-08-23

趙考生（1963-），男（漢族），河北任縣人，高級(jí)工程師，主要從事水利水電工程技術(shù)管理工作，（Tel）13502109302。

（責(zé)任編輯：王艷肖）