基于HBV模型的開(kāi)都河致災(zāi)洪水臨界雨量分析

樊靜，江遠(yuǎn)安，陳穎，陳鵬翔，白素琴，余行杰

（新疆氣候中心，新疆烏魯木齊 830002）

基于HBV模型的開(kāi)都河致災(zāi)洪水臨界雨量分析

樊靜，江遠(yuǎn)安，陳穎，陳鵬翔，白素琴，余行杰

（新疆氣候中心，新疆烏魯木齊 830002）

應(yīng)用HBV水文模型對(duì)開(kāi)都河流域焉耆站以上段進(jìn)行了模擬，通過(guò)對(duì)模型敏感參數(shù)的率定（1982—1985年），Nash效率系數(shù)達(dá)到0.73，說(shuō)明模型能夠很好地反映開(kāi)都河流域的徑流過(guò)程；對(duì)模型進(jìn)行驗(yàn)證的結(jié)果表明在1986—1989年模擬結(jié)果的Nash效率系數(shù)達(dá)到0.75，比率定期的Nash效率系數(shù)高，模型在開(kāi)都河流域有較高的適用性，能夠反映出流域年、月、日徑流量的變化趨勢(shì)。通過(guò)率定后的HBV模型，計(jì)算了不同流量和水位所對(duì)應(yīng)的臨界致災(zāi)雨量，當(dāng)基礎(chǔ)水位為6.27m，流量為100m3/s時(shí)，臨界一級(jí)（保證流量800m3/s)和二級(jí)(警戒流量500m3/s）致災(zāi)面雨量分別為50mm和30mm，隨著基礎(chǔ)水位和流量的增加，臨界致災(zāi)雨量的值在減小。將通過(guò)HBV水文模型計(jì)算的臨界雨量與歷史上發(fā)生洪水的24 h和120 h的實(shí)際降水量進(jìn)行對(duì)比分析，發(fā)現(xiàn)應(yīng)用HBV水文模型推算的臨界面雨量大于實(shí)際24 h降水量，小于前5 d的降水合計(jì)量，因此在實(shí)際的業(yè)務(wù)應(yīng)用中還需要對(duì)結(jié)果進(jìn)行進(jìn)一步的分析和驗(yàn)證。

HBV模型；開(kāi)都河；臨界雨量

山洪災(zāi)害形成的機(jī)理非常復(fù)雜，諸多因素中降雨是引發(fā)山洪災(zāi)害最直接的外在動(dòng)力因素[1]。當(dāng)一個(gè)小流域某時(shí)段內(nèi)降雨量達(dá)到或者超過(guò)某一量級(jí)和強(qiáng)度時(shí)，形成的洪水流量剛好為河道的安全泄洪能力，大于該降雨量將可能引發(fā)山洪災(zāi)害，把此時(shí)的降雨量稱為臨界雨量[2]。臨界雨量是目前山洪災(zāi)害預(yù)測(cè)預(yù)報(bào)的重要指標(biāo)，對(duì)山洪災(zāi)害防治有著重要的意義[3]。對(duì)于臨界雨量的確定，根據(jù)有無(wú)實(shí)測(cè)水文和降水資料，可采用不同的計(jì)算方法。山洪災(zāi)害臨界雨量分析方法主要采用統(tǒng)計(jì)歸納法[4]、暴雨等值線分析法和地質(zhì)與地貌條件綜合分析法[5]、水位反推法[6]、山洪災(zāi)害實(shí)例調(diào)查法[7]、災(zāi)害與降雨頻率分析法[8]等。近期張亞萍等[9]提出一種基于水文模擬建立中小河流洪水氣象風(fēng)險(xiǎn)等級(jí)臨界累積面雨量指標(biāo)的方法，對(duì)于有水位、流量等監(jiān)測(cè)資料的中小河流，可通過(guò)水文模擬方法，在進(jìn)行洪水氣象風(fēng)險(xiǎn)等級(jí)指標(biāo)確定時(shí)考慮流域的水文響應(yīng)。目前，新疆區(qū)域關(guān)于臨界雨量的研究較少。本文以開(kāi)都河流域?yàn)槔诹私馄渲聻?zāi)降雨特征的基礎(chǔ)上，利用HBV水文模型法計(jì)算流域的臨界致災(zāi)雨量并選取1960年以來(lái)的多次洪水過(guò)程對(duì)確定的臨界雨量進(jìn)行驗(yàn)證，探討計(jì)算得出的臨界雨量的合理性。通過(guò)此項(xiàng)研究，以期為開(kāi)展洪水氣象預(yù)警業(yè)務(wù)提供參考，也為新疆區(qū)域臨界雨量的確定探索新方法。同時(shí)，很多中小河流并無(wú)水文監(jiān)測(cè)資料，對(duì)有水文監(jiān)測(cè)資料的流域進(jìn)行臨界累積面雨量指標(biāo)的研究并與地理信息結(jié)合，是進(jìn)一步研究無(wú)水文監(jiān)測(cè)資料流域洪水氣象風(fēng)險(xiǎn)指標(biāo)的重要基礎(chǔ)[10-11]。

1 研究區(qū)概況

開(kāi)都河流域是塔里木河流域的主要源流之一，位于新疆天山南麓，塔里木盆地北緣，屬于雪冰融水和雨水混合補(bǔ)給的河流。開(kāi)都河發(fā)源于天山中部山區(qū)，上游流經(jīng)巴音布魯克草原，下游流經(jīng)焉耆盆地，最后注入中國(guó)最大的內(nèi)陸淡水湖——博斯騰湖，河流全長(zhǎng)560 km。深居歐亞大陸腹地，屬大陸性溫帶干旱氣候，日照充足，太陽(yáng)輻射強(qiáng)。夏季炎熱，冬季寒冷，晝夜溫差及年較差大，降水少，蒸發(fā)能力強(qiáng)烈，空氣干燥。開(kāi)都河流域自北至南高差較大，大致可劃分為北部天山高寒半干旱半濕潤(rùn)區(qū)，海拔2 500～ 5 000m，年平均氣溫-3.5～6.5℃，年降水量300～600mm；中部是低山丘陵及沖積扇戈壁干旱區(qū)，海拔1 500～2 500m，年平均氣溫7.3～9.2℃，年降水量100～300mm；南部戈壁平原溫帶干旱區(qū)，海拔1 000～ 1 500m，年平均氣溫9.2～11.5℃，年降水量70～100mm。

2 資料和方法

2.1 資料

氣象資料：分布在開(kāi)都河流域內(nèi)和周邊的6個(gè)氣象站（巴音布魯克、大西溝、巴輪臺(tái)、焉耆、和碩、和靜站）1960年1月1日—2010年12月31日的降水量、日平均氣溫、日最高氣溫和日最低氣溫等氣象資料。

水文資料：焉耆水文站1982年1月1日—1989年12月31日的日平均流量和水位數(shù)據(jù)，焉耆水文站的警戒流量和保證流量等指標(biāo)數(shù)據(jù)。1960—1989年焉耆站超過(guò)警戒流量（500m3/s）和保證流量（800m3/s）的典型洪水過(guò)程的日均流量和水位資料。

其他資料：全疆1:5萬(wàn)DEM數(shù)據(jù)；全國(guó)土壤持水力和土地利用數(shù)據(jù)；全疆暴雨洪澇相關(guān)歷史災(zāi)情、實(shí)時(shí)災(zāi)情。

2.2 方法

采用HBV水文模型對(duì)開(kāi)都河流域焉耆站以上段的水文過(guò)程進(jìn)行模擬，對(duì)水文模型進(jìn)行率定和驗(yàn)證，通過(guò)水文模型驗(yàn)證的雨—洪關(guān)系求得開(kāi)都河流域的臨界面雨量。

2.2.1 HBV水文模型

HBV模型是20世紀(jì)70年代瑞典氣象水文研究院（SMHI）開(kāi)發(fā)研制的水文預(yù)報(bào)模型，是一種基于DEM劃分子流域模擬積雪、融雪、實(shí)際蒸散量、土壤水分儲(chǔ)存、地下水埋深和徑流等機(jī)制的半分布式的降雨徑流模型[11]，經(jīng)過(guò)幾十年的發(fā)展,模型的不同版本已經(jīng)在全世界100多個(gè)國(guó)家得到成功應(yīng)用[12]。

HBV模型具有靈活易用的特點(diǎn)，模型的參數(shù)數(shù)據(jù)（如降雨、溫度、流量等）是以二進(jìn)制文件的形式存放于數(shù)據(jù)庫(kù)中，數(shù)據(jù)的轉(zhuǎn)換和共享方便。HBV模型預(yù)報(bào)的預(yù)見(jiàn)期可以是1～24 h，可利用模型計(jì)算降水量、降雪量、土壤的含水量及其變化過(guò)程、流域或水庫(kù)的入流和出流等。實(shí)際應(yīng)用證明HBV水文模型是一個(gè)經(jīng)過(guò)良好測(cè)試的、易操作的洪水預(yù)報(bào)工具[13]。

2.2.2 臨界面雨量的確定方法

基于研究流域的水文和氣象資料的情況，整理適用的流量和降水量序列，建立流量與降水量的定量關(guān)系，通過(guò)HBV水文模型模擬的雨—洪關(guān)系并根據(jù)警戒流量、保證流量等指標(biāo)確定一、二級(jí)洪水臨界致災(zāi)雨量，詳細(xì)的計(jì)算過(guò)程見(jiàn)圖1[14]。

圖1 HBV水文模型計(jì)算臨界雨量的流程

3 結(jié)果分析

3.1 開(kāi)都河流域降水的氣候特征分析

選擇巴音布魯克、大西溝、巴輪臺(tái)、焉耆、和碩、和靜站的逐日降水資料作為開(kāi)都河流域的降水進(jìn)行計(jì)算分析，巴音布魯克、大西溝、巴輪臺(tái)、焉耆、和碩、和靜站的年降水量分別為268.6mm、477.5mm、206mm、75.2mm、86.8mm和62.6mm，年際降水變率分別為0.58mm/a、1.87mm/a、1.26mm/a、0.45mm/a、0.71、0.53mm/a；最大日降水量分別為46.3mm、40.3mm、79.7mm、39.7mm、57.3mm、49.5mm。

3.2 利用HBV水文模型計(jì)算致災(zāi)洪水臨界雨量

3.2.1 數(shù)據(jù)的整理和準(zhǔn)備

（1）使用ARCGIS軟件的水文分析模塊從全疆1:50 000 DEM數(shù)據(jù)中提取開(kāi)都河流域并對(duì)子流域進(jìn)行劃分，根據(jù)開(kāi)都河的實(shí)際情況來(lái)調(diào)整子流域面積閾值的大小從而確定合理的子流域個(gè)數(shù)，并對(duì)劃分好的子流域編號(hào)。

（2）統(tǒng)計(jì)氣象站點(diǎn)的日最高氣溫、平均氣溫、日最低氣溫和降水?dāng)?shù)據(jù)，并按照模型輸入?yún)?shù)的要求整理。

（3）雨量的插值使用國(guó)家氣候中心推薦的R軟件，該插值方法為氣象常用普通克里格，通過(guò)雨量插值得到每個(gè)子流域中心點(diǎn)位置的降水量。

（4）計(jì)算子流域的匯流時(shí)間，使用GIS軟件提取每個(gè)子流域的河流長(zhǎng)度，并使用公式flow length*0.012/1 000計(jì)算每個(gè)子流域的匯流時(shí)間。

（5）流域內(nèi)土地利用和土壤持水力使用國(guó)家氣候中心提供的全國(guó)土地利用和土壤持水力柵格數(shù)據(jù)，在GIS軟件中進(jìn)行數(shù)據(jù)重采樣，將土地利用、土壤持水力、子流域和高程?hào)鸥駭?shù)據(jù)的投影、柵格大小以及范圍統(tǒng)一，存儲(chǔ)為二進(jìn)制格式。

3.2.2 模型的率定和驗(yàn)證

整理好的數(shù)據(jù)通過(guò)HBV水文模型的預(yù)處理程序計(jì)算各子流域的集水面積、匯流時(shí)間等參數(shù)，最后運(yùn)行HBV水文模型主程序中對(duì)流域徑流過(guò)程進(jìn)行模擬。當(dāng)模型的各類(lèi)數(shù)據(jù)輸入并運(yùn)行成功后,為了使模型更符合研究區(qū)實(shí)際情況,需要對(duì)一些非物理參數(shù)和難以確定的物理參數(shù)進(jìn)行調(diào)整,即通過(guò)實(shí)測(cè)數(shù)據(jù)對(duì)模型進(jìn)行率定和驗(yàn)證[15]。本文在模擬過(guò)程中，1982—1985年作為模型的率定期，1986—1989年作為模型的驗(yàn)證期。在校準(zhǔn)和驗(yàn)證的過(guò)程中,通常采用Nash系數(shù)（ME）[16]考查模擬結(jié)果與實(shí)測(cè)值的接近程度,評(píng)估模型的適用性。

式中：Qobs為實(shí)測(cè)值，Qsim為模擬值，Qobs實(shí)測(cè)平均值，n為實(shí)測(cè)數(shù)據(jù)個(gè)數(shù)。當(dāng)Qobs=Qsim時(shí)，ME=1。ME越接近于1，表明模型效率越高，適用性越好。若ME為負(fù)值，表明模型模擬平均值比實(shí)測(cè)平均值的可信度低。

通過(guò)反復(fù)調(diào)整模型中敏感參數(shù)的值（表1），對(duì)開(kāi)都河流域參數(shù)濾定NASH系數(shù)（ME）達(dá)到的最高值為0.73。

表1 HBV模型的敏感參數(shù)及參數(shù)率定結(jié)果

完成模型的率定后需要對(duì)模型進(jìn)行驗(yàn)證（圖2），結(jié)果表明在1986—1989年期間，模型模擬結(jié)果的Nash效率系數(shù)達(dá)到0.75，比率定期Nash效率系數(shù)高，說(shuō)明模型在開(kāi)都河流域有較好的適用性，可以用來(lái)對(duì)流域的徑流過(guò)程進(jìn)行模擬。通過(guò)對(duì)率定期和驗(yàn)證期模擬結(jié)果進(jìn)行對(duì)比分析，可以看出模擬結(jié)果整體的趨勢(shì)與實(shí)測(cè)值較吻合，基本能夠反映出年、月、日徑流量的變化，但是對(duì)峰值的模擬存在誤差。

圖2 HBV模型率定結(jié)果（1982—1985年）和驗(yàn)證結(jié)果（1986—1989年）

3.3 臨界雨量計(jì)算

率定后的HBV模型可以反映流域降水—流量的定量關(guān)系，進(jìn)一步結(jié)合焉耆水文站洪水流量與水位的關(guān)系（圖3），建立降水—流量—水位之間的關(guān)系，從而推算河水達(dá)到警戒流量、保證流量時(shí)的臨界面雨量。

圖3 焉耆站流量—水位概化關(guān)系圖

通過(guò)對(duì)HBV水文模型的驗(yàn)證，當(dāng)模擬的流量達(dá)到警戒、保證流量時(shí)，即認(rèn)定此時(shí)的面雨量為所對(duì)應(yīng)洪水等級(jí)的臨界（面）雨量。由于開(kāi)都河歷次洪水過(guò)程的起漲水位一般都在6.0m以上，通過(guò)對(duì)收集到的洪水資料進(jìn)行統(tǒng)計(jì)后最終以6.27m作為起始水位，將不同前期水位與雨量代入到模型中，再以焉耆站警戒水位7.67m、保證水位8.93m作為臨界判別條件，分別得出不同前期水位下對(duì)應(yīng)的各級(jí)臨界（面）雨量值（表2），需要注意的是這里的雨量值是焉耆站以上流域面雨量。從表2可以看出，臨界致災(zāi)雨量與流域前期的流量直接相關(guān)，前期流量越大，臨界致災(zāi)雨量就越小。

表2 開(kāi)都河流域臨界（面）雨量表

3.41960年以來(lái)流域洪水過(guò)程臨界雨量驗(yàn)證

查閱歷史上焉耆水文站逐日日均流量超警戒（保證）流量的洪水次數(shù)，焉耆水文站共出現(xiàn)9次（表3），均超過(guò)警戒流量但未達(dá)到保證流量。

表3 焉耆水文站出現(xiàn)洪峰流量超警戒（保證）流量的洪水過(guò)程

根據(jù)陳穎等[17]對(duì)烏魯木齊河流域的研究結(jié)果，在天山山區(qū)五天累計(jì)降水量與洪水的相關(guān)性最高。因此對(duì)超過(guò)警戒流量但未達(dá)到保證流量的洪水過(guò)程進(jìn)行分析，分別計(jì)算24 h降水量和前5 d的合計(jì)降水量，并將其與模擬計(jì)算降水量進(jìn)行對(duì)比分析（表3）。從表3可以看出，模擬計(jì)算的降水量較24 h實(shí)際降水值偏大，同時(shí)小于前5 d的降水合計(jì)，說(shuō)明模型考慮了流域徑流的匯流過(guò)程，但是在匯流時(shí)間的時(shí)效上還需進(jìn)一步的驗(yàn)證。需要說(shuō)明的是表3中的9次洪水日數(shù)反映并非累計(jì)降水量最大對(duì)應(yīng)流量最大，其中最大24 h降水量為38.39mm，對(duì)應(yīng)的日均流量為518m3/s；最小24 h降水量不足2mm，但是由于前期累計(jì)降水量大，5 d合計(jì)降水量達(dá)到47.25mm，對(duì)應(yīng)的日均流量為511m3/s，說(shuō)明前期的河流流量和水位與致災(zāi)臨界雨量直接相關(guān)，臨界面雨量值應(yīng)該是動(dòng)態(tài)調(diào)整的。

4 結(jié)論

通過(guò)HBV水文模型在開(kāi)都河流域的應(yīng)用，表明模型能夠很好地模擬開(kāi)都河流域徑流過(guò)程，可以反映出流域年、月、日徑流量的變化趨勢(shì)，說(shuō)明模型在西北區(qū)域具有較高的適用性。應(yīng)用HBV水文模型，計(jì)算了不同流量和水位所對(duì)應(yīng)的臨界致災(zāi)面雨量，當(dāng)基礎(chǔ)水位為6.27m，流量為100m3/s時(shí)，臨界一級(jí)（保證流量800m3/s）和二級(jí)（警戒流量500m3/s）致災(zāi)雨量分別為50mm和30mm，隨著基礎(chǔ)水位和流量的增加，臨界致災(zāi)雨量的值在減少，前期的河流流量和水位與致災(zāi)臨界雨量直接相關(guān)，致災(zāi)臨界雨量值應(yīng)該是動(dòng)態(tài)調(diào)整的。將通過(guò)HBV水文模型計(jì)算的臨界雨量與歷史上發(fā)生洪水的24 h和120 h實(shí)際降水量進(jìn)行對(duì)比分析，發(fā)現(xiàn)應(yīng)用HBV水文模型推算的臨界面雨量大于實(shí)際24 h降水量，同時(shí)小于前5 d的降水合計(jì)量，因此在實(shí)際的業(yè)務(wù)應(yīng)用中還需要對(duì)結(jié)果進(jìn)行進(jìn)一步的分析和驗(yàn)證。

HBV水文模型集降雨徑流和融雪徑流模擬于一體,通過(guò)概化考慮流域高程、植被等下墊面的差異，從而開(kāi)展降雨徑流模擬。HBV模型在開(kāi)都河流域的成功應(yīng)用，對(duì)于西北地區(qū)有融雪性洪水補(bǔ)給的河流徑流模擬來(lái)說(shuō)具有借鑒意義。但是，從模型模擬的結(jié)果來(lái)看，夏季的模擬效果好于冬季，可能與開(kāi)都河流域復(fù)雜的徑流過(guò)程有關(guān)，因?yàn)樵谙募緦?duì)開(kāi)都河流域徑流影響最大的因素是降水，而在冬、春季對(duì)開(kāi)都河流域徑流影響的主要因素是氣溫和積雪面積[18]。此外，在模型的應(yīng)用過(guò)程中，主要用收集到的日均流量資料來(lái)對(duì)模型進(jìn)行率定和校正，缺乏洪峰資料，影響了模型對(duì)洪水峰值的模擬精度，這可能是造成致災(zāi)臨界雨量值偏大的主要原因之一。因此，未來(lái)需要收集典型的洪水過(guò)程資料，進(jìn)一步提高模型在夏季對(duì)洪水過(guò)程的模擬精度。

對(duì)于新疆區(qū)域而言，河流流域面積較大，使用水文模型模擬最大的不足在于流域內(nèi)氣象和水文資料的缺乏，區(qū)域自動(dòng)氣象觀測(cè)站觀測(cè)時(shí)間較短且數(shù)據(jù)質(zhì)量沒(méi)有保證，國(guó)家一般氣象站數(shù)量少、分布稀疏，而且多數(shù)氣象站點(diǎn)分布在平原區(qū)，對(duì)于流域上游和山區(qū)氣象狀況的代表性不強(qiáng)，因此后期通過(guò)提升區(qū)域自動(dòng)觀測(cè)站數(shù)據(jù)質(zhì)量并將其引入模型，會(huì)進(jìn)一步提升模擬效果。

[1]曲曉波，張濤，劉鑫華，等.舟曲“8·8”特大山洪泥石流災(zāi)害氣象成因分析[J].氣象，2010，36（10）：102-105.

[2]全國(guó)山洪災(zāi)害防治領(lǐng)導(dǎo)小組.山洪災(zāi)害臨界雨量分析計(jì)算細(xì)則[S].2003.

[3]樊建勇，單九生，管珉，等.江西省小流域山洪災(zāi)害臨界雨量計(jì)算分析[J].氣象，2012，38（9）：1110-1114.

[4]陳桂亞，袁雅鳴.山洪災(zāi)害臨界雨量分析計(jì)算方法研究[J].水資源研究，2004，25（4）：36-40.

[5]張玉龍，王龍，李靖，等.云南省山洪災(zāi)害臨界雨量空間插值分析方法研究[J].云南農(nóng)業(yè)大學(xué)學(xué)報(bào)，2007，22（4）：570-573.

[6]王仁喬，周月華，王麗，等.湖北省山洪災(zāi)害臨界雨量及降雨區(qū)劃研究[J].高原氣象，2006：25（2）：330-334.

[7]葉勇，王振宇，范波芹.浙江省小流域山洪災(zāi)害臨界雨量確定方法分析[J].水文，2008：28（1）.

[8]譚萬(wàn)沛，韓慶玉.四川省泥石流預(yù)報(bào)的區(qū)域臨界雨量指標(biāo)研究[J].災(zāi)害學(xué)，1992，（2）：37-42.

[9]張亞萍，劉德，廖峻，等.一種基于水文模擬建立中小河流洪水氣象風(fēng)險(xiǎn)等級(jí)指標(biāo)的方法[J].暴雨災(zāi)害，2012，31（4）：351-357.

[10]張建云，何惠.應(yīng)用地理信息進(jìn)行無(wú)資料地區(qū)流域水文模擬研究[J].水科學(xué)進(jìn)展，1998，9（4）：345-350.

[11]Bergstrom S.The HBVmodel[M].In：Singh V P. ComputerModelsofWatershedHydrology.Water Resources Publications，1995：443-476.

[12]李霖.HBV水文預(yù)報(bào)模型及與之集成的水文模型系統(tǒng)介紹[J].水利水文自動(dòng)化，2004：6.

[13]高超，翟建青，陶輝，等.巢湖流域土地利用/覆被變化的水文效應(yīng)研究[J].自然資源學(xué)報(bào)，2009，24（10）1794-1802.

[14]劉蘇峽，夏軍，莫興國(guó).無(wú)資料流域水文預(yù)報(bào)(PUB計(jì)劃)研究進(jìn)展[J].水利水電技術(shù)，2005，36（2）：9-12.

[15]金婧靚，王飛兒.SWAT模型及其應(yīng)用與改進(jìn)的研究進(jìn)展[J].東北林業(yè)大學(xué)學(xué)報(bào)，2010，38（12）：111-114.

[16]Nash J E，Sutcliffe J V.River flow forecasting through conceptualmodels part I:A discussion of principles[J]. Journal of Hydrology，1970，10（3）：282-290.

[17]陳穎，陳鵬翔，江遠(yuǎn)安，等.HBV水文模型與統(tǒng)計(jì)模型計(jì)算烏魯木齊河臨界致災(zāi)雨量的對(duì)比分析[J].暴雨洪澇災(zāi)害風(fēng)險(xiǎn)普查論文集（待出版）.

[18]李倩，李蘭海，包安明.開(kāi)都河流域積雪特征變化及其與徑流的關(guān)系[J].2012，34（1）：91-97.

The Critical Rainfall Calculation in Kaidu River Based on HBV Hydrologicalmodel

FAN Jing，JIANG Yuanan，CHEN Ying，CHEN Pengxiang，BAI Suqin，YU Xingjie
（Xinjiang Climate Center，Urumqi 830002，China）

The HBV hydrologicalmodel was used to simulate the runoff of Kaidu Riverabove the Yanqi hydrologic station in this paper.Based on HBVmodel calibrationand correction,the results show that the runoff process was basicallyaccord with themeasured which can reflect changes of runoff trend in the year,monthand daily.But thereare certain errors on the peak runoff of the simulation results.Through the HBVmodel simulation,calculation the different flowand water level correspond to the critical rainfall.Meanwhile,the paper compared HBV calculation results to 24 hour/120 hour rainfall when the river floods in history.The results show that rainfall by HBV hydrologicalmodel calculation is greater than theactual 24 hour precipitation,maybeassociated with the simulation of flood peak flow.

HBV hydrologicalmodel;Kaidu River;critical Rainfall

S421

B

1002-0799（2014）06-0031-05

10.3969/j.issn.1002-0799.2014.06.005

2014-09-06；

2014-10-14

新疆氣象局面上項(xiàng)目（MS201402）資助。

樊靜（1980-），女，工程師，現(xiàn)從事氣候變化分析及其影響評(píng)估工作。E-mail：fanjing365@126.com

樊靜，江遠(yuǎn)安，陳穎，等.基于HBV模型的開(kāi)都河致災(zāi)洪水臨界雨量分析[J].沙漠與綠洲氣象，2014,8（6）：31-35.