李　明

(陜西省氣象臺(tái), 陜西 西安710015)

?

陜南中小河流洪水氣象統(tǒng)計(jì)預(yù)報(bào)模型業(yè)務(wù)化研究與應(yīng)用*

李明

(陜西省氣象臺(tái), 陜西 西安710015)

摘要:在確定中小河流選取標(biāo)準(zhǔn)及其洪水氣象預(yù)報(bào)等級(jí)的基礎(chǔ)上，面向中小河流洪水氣象預(yù)報(bào)業(yè)務(wù)的需求提出一種統(tǒng)計(jì)預(yù)報(bào)方法，即根據(jù)雷達(dá)定量降水估測的流域面雨量、前期影響雨量方法確定流域的徑流深，綜合考慮中小河流的前1 d水位、徑流深、當(dāng)日水位，利用2013-2014年的6-9月陜南中小河流的水文觀測站逐日水位、流域內(nèi)自動(dòng)氣象站降雨量、雷達(dá)回波資料，建立以水位預(yù)報(bào)為目標(biāo)的二元回歸方程。以子午河為例建立的洪水氣象預(yù)報(bào)方程和其他中小河流在2015年汛期的業(yè)務(wù)應(yīng)用結(jié)果表明：該統(tǒng)計(jì)方法建立的洪水氣象預(yù)報(bào)方程和實(shí)際情況較為吻合，能較為準(zhǔn)確、簡潔地預(yù)報(bào)中小河流發(fā)生洪水的氣象可能性等級(jí)。

關(guān)鍵詞:陜南地區(qū)；中小河流；洪水；氣象預(yù)報(bào)；二元回歸

我國的大江大河經(jīng)過幾十年的基礎(chǔ)設(shè)施建設(shè)已具備了較強(qiáng)的防汛抗洪能力和一定的預(yù)報(bào)水平，取得較好社會(huì)效益。但是針對(duì)中小河流洪水的氣象預(yù)報(bào)預(yù)警技術(shù)研究和業(yè)務(wù)應(yīng)用很少，防洪建設(shè)亦較薄弱、防洪標(biāo)準(zhǔn)較低，而且其洪水產(chǎn)生突發(fā)性強(qiáng)、歷時(shí)短，造成的損失常常比較嚴(yán)重，所以中小河流洪水的預(yù)報(bào)和防御已成為當(dāng)前防洪減災(zāi)的重要工作。為了提高中小河流洪水氣象預(yù)報(bào)預(yù)警水平，氣象部門從2011年起開展針對(duì)強(qiáng)降水誘發(fā)中小河流洪水的氣象預(yù)警技術(shù)的研究和業(yè)務(wù)應(yīng)用工作，將傳統(tǒng)的天氣預(yù)報(bào)向氣象災(zāi)害預(yù)報(bào)預(yù)警這一領(lǐng)域推進(jìn)，以期有效提升防災(zāi)減災(zāi)科學(xué)的決策水平。強(qiáng)降水誘發(fā)的中小河流洪水的氣象預(yù)報(bào)技術(shù)應(yīng)用研究，其核心是流域定量降水的確定以及定量降水和洪水之間關(guān)系的確定，氣象預(yù)報(bào)工作者深入探討了定量降水估測和預(yù)報(bào)技術(shù)在水文預(yù)報(bào)中的應(yīng)用[1-7]，水文模型是確定定量降水和洪水關(guān)系的重要手段之一[8-9]，近年來眾多水文、氣象工作者在水文模型應(yīng)用方面進(jìn)行了大量的研究[10-18]，但是由于水文模型在結(jié)構(gòu)上與實(shí)際水文空間分散性和不均勻性輸入不匹配，分布式水文模型水文單元的劃分、空間參數(shù)的確定、產(chǎn)匯流機(jī)制的確定和模型算法的實(shí)現(xiàn)等較為復(fù)雜，因此各種水文模型具有不確定性，效果并不令人滿意，因此對(duì)該方面進(jìn)行的理論研究比較多，而實(shí)際業(yè)務(wù)應(yīng)用比較少。因此，一些學(xué)者探討了統(tǒng)計(jì)方法在洪水預(yù)報(bào)中的應(yīng)用[19-26]，統(tǒng)計(jì)方法是黑箱問題，雖然不著重考慮洪水產(chǎn)生的具體物理過程，但是也能夠得到較為理想的預(yù)報(bào)結(jié)果,且操作簡單快捷,因此是實(shí)際業(yè)務(wù)工作中開展洪水預(yù)報(bào)的有效途徑之一。本文從業(yè)務(wù)需求和可操作性出發(fā)，利用多年的氣象、水文樣本資料，提出了一種簡潔的中小河流洪水氣象預(yù)報(bào)統(tǒng)計(jì)方法，經(jīng)檢驗(yàn)，此方法對(duì)中小河流的洪水預(yù)報(bào)能夠取得較為理想的預(yù)報(bào)結(jié)果。

1資料與方法

1.1資料

研究選取的資料為陜南秦巴山區(qū)中小河流水文站2013-2014年汛期6-9月逐日的水位、流域內(nèi)自動(dòng)氣象站的降水觀測資料、雷達(dá)回波資料。文中時(shí)間均為北京時(shí)間。

1.2方法

首先根據(jù)中小河流的警戒、保證、漫堤(壩)水位和普查情況，以水位為標(biāo)準(zhǔn)，將中小河流洪水的氣象預(yù)報(bào)等級(jí)分為4級(jí)；根據(jù)中小河流的當(dāng)日水位、前1 d水位、徑流深建立以水位預(yù)報(bào)為目標(biāo)的二元回歸方程；利用二元回歸方程，根據(jù)前1 d水位、預(yù)報(bào)徑流深實(shí)現(xiàn)未來24 h的預(yù)報(bào)水位，比較預(yù)報(bào)水位和各中小河流洪水的4個(gè)等級(jí)標(biāo)準(zhǔn)水位，確定洪水發(fā)生可能性的氣象預(yù)報(bào)等級(jí)。

表1　選取的陜南秦巴山區(qū)中小河流基本信息

表2　陜南中小河流洪水四級(jí)氣象預(yù)報(bào)標(biāo)準(zhǔn)

2中小河流洪水氣象預(yù)報(bào)技術(shù)方法

2.1確定中小河流

中小河流的選取標(biāo)準(zhǔn)：根據(jù)文獻(xiàn)[26]和中國氣象局編寫的《中小河流洪水、山洪、泥石流、滑坡災(zāi)害致災(zāi)臨界(面)雨量確定技術(shù)指南》，將中小河流定義為流域面積大于200 km2小于3 000 km2；此外本文還考慮沒有匯流，有實(shí)時(shí)監(jiān)測的水位、流量信息、警戒和保證水位。根據(jù)該標(biāo)準(zhǔn)，選取了陜南秦巴山區(qū)14條中小河流(表1)，它們均屬于長江流域漢江水系和丹江水系。

2.2中小河流洪水氣象預(yù)報(bào)等級(jí)劃分

水文部門一般從低到高分為警戒、保證和漫堤(壩)水位三個(gè)標(biāo)準(zhǔn)。由于陜南秦巴山區(qū)中小河流河道附近常有農(nóng)田等，考慮對(duì)其影響，結(jié)合陜西氣象部門的普查情況，將低于警戒水位1.3 m作為陜南秦巴山區(qū)中小河流發(fā)生洪水氣象可能性的最低等級(jí)，其他等級(jí)按水文部門的標(biāo)準(zhǔn)確定，陜南中小河流洪水的四級(jí)氣象預(yù)報(bào)標(biāo)準(zhǔn)如表2所示。

2.3中小河流洪水氣象預(yù)報(bào)方程建立

中小河流未來水位的變化和該河流當(dāng)前的水位、流域內(nèi)未來可能產(chǎn)生的徑流深有密切關(guān)系，因此本文采用流域當(dāng)日徑流深、河流水文站前1 d水位作為自變量，水文站當(dāng)日水位作為因變量，采用回歸方法，建立以水位預(yù)報(bào)為目標(biāo)的二元回歸方程。

2.4徑流深因子的確定

陜南秦巴山區(qū)年降水量在800～1 200 mm之間，屬于北亞熱帶濕潤氣候區(qū)，因此該區(qū)域的徑流深計(jì)算采用蓄滿產(chǎn)流方法。

徑流深計(jì)算公式[25]為：

Rs=Ap-(100-Pa)。

(1)

式中：Rs是徑流深，Ap是流域面雨量，Pa是流域前期影響雨量。即徑流深為流域面雨量和土壤蓄水量的折減量之差。若當(dāng)日面雨量Rs≤Ap-(100-Pa)，則當(dāng)日徑流深為0。

2.5流域面雨量(Ap)確定

流域面雨量(簡稱面雨量)是指一定流域范圍內(nèi)的平均降水量。面雨量估算的準(zhǔn)確性直接關(guān)系到洪水預(yù)報(bào)的準(zhǔn)確性。面雨量采用兩種方法：一種是針對(duì)流域已經(jīng)產(chǎn)生降雨的面雨量計(jì)算方法，即采用雷達(dá)回波和自動(dòng)氣象站降水觀測資料相結(jié)合的雷達(dá)定量降水估測(Quantity Precipitation Estimation，QPE)方法，稱之為流域?qū)崨r面雨量，再建立回歸方程計(jì)算歷史逐日徑流深；另一種是針對(duì)流域預(yù)報(bào)降雨的面雨量計(jì)算方法，采用泰森多邊形方法，稱之為流域預(yù)報(bào)面雨量，用于實(shí)時(shí)預(yù)報(bào)時(shí)洪水氣象預(yù)報(bào)方程中流域內(nèi)未來徑流深的計(jì)算。

2.5.1流域?qū)崨r面雨量方法

陜南秦巴山區(qū)地形地貌復(fù)雜，自動(dòng)氣象站點(diǎn)密度有限且布點(diǎn)多以行政區(qū)劃為基礎(chǔ)，不能完全代表流域的降水情況，因此流域的實(shí)況面雨量采用雷達(dá)定量降水估測(QPE)方法，即雷達(dá)回波和自動(dòng)氣象站降水觀測資料融合技術(shù)，空間分辨率達(dá)1 km×1 km。對(duì)流域內(nèi)的所有格點(diǎn)的QPE進(jìn)行算數(shù)平均即得到流域?qū)崨r面雨量。目前全國氣象部門使用的雷達(dá)定量降水估測(QPE)方法是以長江中下游降水情況研發(fā)的，并不完全適用于陜南秦巴山區(qū)，因此需要針對(duì)陜南的降水情況進(jìn)行本地化，技術(shù)方法如下所示。

雷達(dá)定量估測降水的方法有多種，最常用的是Z-I關(guān)系法，即應(yīng)用雷達(dá)氣象方程測得降水回波的反射率因子，然后根據(jù)反射率因子和降水強(qiáng)度之間的關(guān)系來推算降水強(qiáng)度。

(1)采用最優(yōu)化處理法建立Z-I關(guān)系(Z為雷達(dá)反射率因子，I為降水強(qiáng)度)

最優(yōu)化處理法：用雨量計(jì)測得每小時(shí)雨量Gi, 同時(shí)由雷達(dá)測得時(shí)空相對(duì)應(yīng)的dBz值, 目前常用的Z-I關(guān)系式為Z=AIb，先假定上式中A和b的值, 把dBz轉(zhuǎn)換成雨強(qiáng), 再對(duì)雨強(qiáng)進(jìn)行時(shí)間累積得到雷達(dá)估算的每小時(shí)雨量Hi，選取判別函數(shù)為

(2)

式中：i為樣本的序號(hào),是不同時(shí)間的取樣值；N為樣本總數(shù)。不斷調(diào)整A和b 的值, 使CTF2為最小。使CTF2為最小的A和b就是這個(gè)統(tǒng)計(jì)樣本總體的最優(yōu)參數(shù)，即CTF2值為最小的參數(shù)A和b，使得Hi最逼近Gi。

(2)分組統(tǒng)計(jì)求取Z-I關(guān)系

從雷達(dá)探測角度, 直接得到的就是dBz的值,對(duì)dBz值分組求取Z-I關(guān)系, 以及用分組Z-I關(guān)系轉(zhuǎn)換雨強(qiáng),計(jì)算機(jī)處理簡單易行，且不同強(qiáng)度的雷達(dá)回波在一定程度上代表了不同類型的降水。本文按照dBz值大小把樣本分為8組, 即5～20dBz，21～25dBz，26～30dBz，31～35dBz，36～40dBz，41～45dBz，46～50dBz，50dBz以上，分別用最優(yōu)化處理方法得出Z-I關(guān)系。CTF2場分析,取A 參數(shù)為100～400, 步長為1，b參數(shù)為0.5～4.5, 步長為0.01，不斷地改變A和b的值,得到CTF2隨A、b的兩維分布, 使CTF2場最小，則得到了最優(yōu)化Z-I關(guān)系的A、b值。

(3)資料處理

利用2009-2012年汛期降水資料和多普勒天氣雷達(dá)相應(yīng)的6min一次的體掃資料，把自動(dòng)站雨量計(jì)觀測到的1h降水定為一個(gè)樣本。I為1h降雨量，dBz為1h的平均回波強(qiáng)度。

將雷達(dá)回波掃描面上1km×1km像素點(diǎn)位置的極坐標(biāo)變換到直角坐標(biāo)，挑選地面雨量站對(duì)應(yīng)掃描面上像素點(diǎn)的雷達(dá)回波強(qiáng)度。雷達(dá)每6min進(jìn)行一次完整的體掃，1h有10次體掃，即每個(gè)點(diǎn)有10次掃描回波資料，將各像素點(diǎn)10次的強(qiáng)度求平均，作為統(tǒng)計(jì)用的回波強(qiáng)度資料。

根據(jù)陜南安康、漢中多普勒雷達(dá)60次降水過程，按照最優(yōu)化方法定義,當(dāng)判別函數(shù)、平均絕對(duì)誤差、平均相對(duì)誤差最小時(shí),雷達(dá)定量測量降水關(guān)系中的A、b為最佳系數(shù)。得出適合陜南降水的最優(yōu)Z-I關(guān)系(表3)。經(jīng)檢驗(yàn)，該系數(shù)對(duì)雷達(dá)估測陜南降水有較好的改善。

(4)流域?qū)崨r面雨量計(jì)算

根據(jù)本地化的雷達(dá)定量降水估測(QPE)產(chǎn)品，得到流域內(nèi)每個(gè)格點(diǎn)逐小時(shí)的定量降水估測值(QPE)，將流域內(nèi)所有格點(diǎn)的QPE累加，并進(jìn)行算術(shù)平均，得到逐小時(shí)的流域面雨量，然后將逐小時(shí)的面雨量相加得到流域內(nèi)24h實(shí)況面雨量。

2.5.2流域預(yù)報(bào)面雨量方法

流域的降水預(yù)報(bào)是業(yè)務(wù)人員在客觀降水預(yù)報(bào)產(chǎn)品基礎(chǔ)上主觀訂正的鄉(xiāng)鎮(zhèn)站點(diǎn)降水預(yù)報(bào)值，由于站點(diǎn)的不均勻性，面雨量估算需要利用離散的點(diǎn)雨量推算流域內(nèi)的平均降水量，本文流域預(yù)報(bào)面雨量采用泰森多邊形方法。

泰森多邊形方法：根據(jù)流域的數(shù)字高程地圖(DEM),利用ArcGis軟件生成河流水系數(shù)據(jù)文件，截取所需河流的流域，再根據(jù)流域范圍內(nèi)及周圍的雨量站點(diǎn)，利用ArcGis制作流域所在區(qū)域的泰森多邊形，最后采用權(quán)重法計(jì)算流域面雨量，即：

表3　本地化的陜南區(qū)域雷達(dá)估測降水Z-I關(guān)系系數(shù)

(3)

式中:Ap為流域面雨量，Ri為站點(diǎn)i的雨量，Ai為站點(diǎn)i所代表的面積，A為流域總面積，n為泰森多邊形個(gè)數(shù)。

2.6前期影響雨量計(jì)算

2.6.1前期影響雨量公式

降雨開始時(shí)，流域的濕潤程度對(duì)降雨產(chǎn)生徑流量影響很大。因此，在產(chǎn)流計(jì)算中需要考慮這一因素。本文采用前期影響雨量Pa作為衡量流域干濕程度的指標(biāo)，以反映流域蓄水量的大小。前期影響雨量大，則凈雨多，徑流大；反之，則凈雨少，徑流也小。

Pa的計(jì)算公式有多種，本文采用文獻(xiàn)[25]中的算式如下：

Pa,t+1=Ka(Pa,t+Pt)，

(4)

且控制Pa,t+1≤wm。

(5)

式中：Pa,t，Pa,t+1分別為第t天和第t+1天開始時(shí)的前期影響雨量(mm)。Pt為第t天的流域降雨量(mm)，根據(jù)流域?qū)崨r面雨量方法即雷達(dá)定量降水估測(QPE)方法計(jì)算確定。Ka為流域蓄水量的折減系數(shù)，各月近似取一個(gè)平均值，等于1-Em/Wm，Em為流域平均日蒸散量，Wm為流域最大蓄水量。

2.6.2前期影響雨量公式中折減系數(shù)(Ka)的確定

陜西水文部門一般把陜南秦巴山區(qū)土壤蓄水量Wm取為100mm，利用衛(wèi)星遙感技術(shù)確定流域蒸散量是有效手段之一[27-28]，根據(jù)文獻(xiàn)[28]中對(duì)陜南秦巴山區(qū)蒸散量的研究成果，利用折減系數(shù)公式

Ka=1-Em/Wm，

(6)

計(jì)算出各月的日平均折減系數(shù)，如表4所示。

表4　陜南秦巴山區(qū)汛期各月多年平均日蒸散量和折減系數(shù)

2.6.3前期影響雨量計(jì)算

計(jì)算前期影響雨量Pa，首先確定流域內(nèi)土壤蓄水量達(dá)到100 mm的具體日期：流域內(nèi)第一次出現(xiàn)暴雨及以上量級(jí)天氣且中小河流水位上漲日即為流域內(nèi)土壤蓄水量達(dá)到100 mm的時(shí)間，即Pa=100 mm。然后從土壤蓄水量達(dá)100 mm的日期起，根據(jù)每日的實(shí)況面雨量(Pt)、前一日的前期影響雨量(Pa,t)，根據(jù)前期影響雨量公式，計(jì)算出每日的前期影響雨量(Pa,t+1)。

3中小河流洪水氣象預(yù)報(bào)方程-以子午河為例

3.1子午河流域

子午河(圖1)屬于漢江的支流，水文站為兩河口，警戒水位、保證水位、漫壩水位分別為546.8 m，548.5 m，551.8 m，兩河口站以上子午河流域總面積2 818 km2，主要包括佛坪和寧陜2個(gè)縣級(jí)行政區(qū)，流域大部分處于寧陜縣，區(qū)域內(nèi)地形復(fù)雜，自然植被良好，流域內(nèi)多年平均降水量900 mm左右，是典型的亞熱帶濕潤氣候區(qū)。

圖1　子午河流域(藍(lán)色區(qū)域)在陜南的位置示意圖

統(tǒng)計(jì)資料樣本為2013-2014年2年內(nèi)的6-9月的子午河逐日水位，以及對(duì)應(yīng)的流域內(nèi)5個(gè)自動(dòng)氣象觀測站逐小時(shí)降雨、漢中多普勒天氣雷達(dá)每6 min一次的回波資料，樣本共計(jì)78個(gè)。

3.2預(yù)報(bào)方程的建立

該流域土壤蓄水量達(dá)100 mm的日期：2013年為5月29日，2014年為5月26日。根據(jù)前期影響雨量公式，逐日計(jì)算2013年5月29日、2014年5月26日蓄水量為100 mm后各日的前期影響雨量。建立流域洪水氣象預(yù)報(bào)方程所需的徑流深(x1)和前一日水位(x2)、當(dāng)日水位(Y)數(shù)據(jù)，如表5所示。

根據(jù)二元回歸方法，子午河水位預(yù)報(bào)方程為：

Y=260.10+0.08x1+0.52x2。

(7)

式中：Y為預(yù)報(bào)水位；x1為徑流深；x2為前1 d水位。統(tǒng)計(jì)分析表明相關(guān)系數(shù)達(dá)0.94。

3.3預(yù)報(bào)方程檢驗(yàn)

以2015年6月29日子午河洪水為例。2015年6月28日08:00-29日08:00，子午河流域普降暴雨、大暴雨，子午河流域29日08:00，兩河口水位站流量達(dá)3 960 m3/s，水位達(dá)550.71 m，超保證水位(548.5 m)2.21 m。本次強(qiáng)降水過程主要是本地生成的一個(gè)中β尺度對(duì)流云團(tuán)在佛坪、寧陜?cè)斐傻拇蟊┯?，在子午河產(chǎn)生洪水，該對(duì)流云團(tuán)(圖2a)孤立存在，周圍都是弱層狀云區(qū)，強(qiáng)降水區(qū)正好覆蓋整個(gè)子午河流域，降水持續(xù)近12 h，從28日08:00-29日08:00子午河流域降水量分布圖(圖2b圓圈處)，可以明顯看出在子午河流域大部分區(qū)域降雨達(dá)100 mm以上，和對(duì)流云團(tuán)的位置基本一致，降水強(qiáng)度大。

子午河流域面雨量確定：根據(jù)泰森多邊形方法計(jì)算確定28日08:00-29日08:00子午河流域的面雨量為107.3 mm。

表5　建立子午河流域洪水氣象預(yù)報(bào)方程所需樣本主要因子的數(shù)據(jù)

圖2　造成子午河流域強(qiáng)降水的對(duì)流云團(tuán)和6月28日08:00-29日08:00陜南累積降水分布(單位：mm)

子午河流域蓄水量為100 mm的日期確定：2015年6月24日08:00-25日08:00子午河流域出現(xiàn)大到暴雨，部分區(qū)域出現(xiàn)暴雨，河流水位明顯上漲，因此確定25日子午河流域土壤蓄水量達(dá)最大100 mm；根據(jù)前期影響雨量公式，29日流域的前期影響雨量為100 mm。

徑流深：根據(jù)徑流深公式(1)計(jì)算出6月29日徑流深為107.3 mm。

子午河29日08:00水位計(jì)算：前1 d 28日08:00水位為543.82 m，29日徑流深為107.3 mm。因此，根據(jù)3.2節(jié)確定的子午河流域洪水氣象預(yù)報(bào)方程，29日08時(shí)預(yù)報(bào)水位Y=260.10+0.08×107.3+0.52×543.82，Y為551.45 m。

29日08:00觀測水位為550.71 m，計(jì)算的29日08時(shí)預(yù)報(bào)水位高出觀測水位0.74 m，根據(jù)2.2節(jié)中小河流洪水等級(jí)劃分標(biāo)準(zhǔn)，子午河洪水的四級(jí)氣象預(yù)報(bào)標(biāo)準(zhǔn)如表6所示。根據(jù)表6可知，本次洪水氣象預(yù)報(bào)等級(jí)為Ⅱ級(jí)，預(yù)報(bào)方程的計(jì)算等級(jí)和實(shí)況等級(jí)一致。

4其他中小河流預(yù)報(bào)效果

利用該方法，在2015年汛期對(duì)陜南中小河流洪水氣象預(yù)報(bào)業(yè)務(wù)試用，結(jié)果表明共有18條次中小河流發(fā)生洪峰，預(yù)報(bào)出現(xiàn)洪水12條次，空?qǐng)?bào)2條次，漏報(bào)4條次，準(zhǔn)確率66.7%。從預(yù)報(bào)等級(jí)來看，12次預(yù)報(bào)出現(xiàn)洪水，8次預(yù)報(bào)等級(jí)和實(shí)況等級(jí)一致，3次等級(jí)偏低一個(gè)等級(jí)，1次偏高一個(gè)等級(jí)。因此從2015年汛期應(yīng)用來看，該方法建立的中小河流洪水氣象預(yù)報(bào)方程具有較好的可應(yīng)用性。

表6　子午河洪水四級(jí)等級(jí)標(biāo)準(zhǔn)

5結(jié)論與討論

面向當(dāng)前業(yè)務(wù)需求，根據(jù)標(biāo)準(zhǔn)確定中小河流及其洪水四級(jí)氣象預(yù)報(bào)等級(jí)。綜合考慮中小河流的前1 d水位、徑流深和當(dāng)日水位，建立以水位預(yù)報(bào)為目標(biāo)的中小河流洪水二元回歸方程，并探討了建立方程時(shí)所需的面雨量、前期影響雨量、徑流深等確定方法。在2015年汛期，把文中所確定的這套方法應(yīng)用到了陜南14條中小河流的洪水氣象預(yù)報(bào)業(yè)務(wù)，結(jié)果表明該方法對(duì)無匯流的中小河流洪水氣象預(yù)報(bào)具有一定的可用性，在實(shí)際應(yīng)用中方便、快捷。同時(shí)該方法的漏報(bào)率高于空?qǐng)?bào)率，這需要在今后進(jìn)行兩方面的改進(jìn)：一是降水預(yù)報(bào)的準(zhǔn)確性特別是落區(qū)預(yù)報(bào)需要更為準(zhǔn)確，這直接影響中小河流洪水的氣象預(yù)報(bào)準(zhǔn)確率；二是由于建立洪水氣象預(yù)報(bào)方程所采用的資料時(shí)間段為2013-2014年，這兩年強(qiáng)降水多為系統(tǒng)性、區(qū)域性過程，而2015年多為短時(shí)性、局地性強(qiáng)降水過程，因此造成24 h內(nèi)流域面上降水量較小，導(dǎo)致漏報(bào)，所以今后需要進(jìn)一步根據(jù)降水性質(zhì)分類別建立方程，提高預(yù)報(bào)準(zhǔn)確率。

致謝：感謝陜西省農(nóng)業(yè)遙感信息中心董金芳提供中小河流示意圖。

參考文獻(xiàn)：

[1]包紅軍,趙琳娜.基于集合預(yù)報(bào)的淮河流域洪水預(yù)報(bào)研究[J]. 水利學(xué)報(bào),2012,43(2):215-224.

[2]李俊,廖移山,張兵,等. 集合數(shù)值預(yù)報(bào)方法在山洪預(yù)報(bào)中的初步應(yīng)用[J].高原氣象,2007,26(4):854-861.

[3]王莉莉,陳德輝,趙琳娜.GRAPES氣象-水文模式在一次洪水預(yù)報(bào)中的應(yīng)用[J].應(yīng)用氣象學(xué)報(bào),2012,23(3):274-284.

[4]彭濤,宋星原,殷志遠(yuǎn),等. 雷達(dá)定量降水估算在水文模式汛期洪水預(yù)報(bào)中的應(yīng)用試驗(yàn)[J].氣象,2010,36(12):50-55.

[5]張亞萍,程明虎,徐慧,等.雷達(dá)定量測量降水在佛子嶺流域徑流模擬中的應(yīng)用[J].應(yīng)用氣象學(xué)報(bào),2007,18(3):295-305.

[6]崔春光,彭濤,沈鐵元,等.定量降水預(yù)報(bào)與水文模型耦合的中小河流域汛期洪水預(yù)報(bào)試驗(yàn)[J].氣象,2010,36(12):56-61.

[7]趙坤,葛文忠,劉國慶,等. 在雷達(dá)測雨和洪水預(yù)報(bào)中自適應(yīng)卡爾曼濾波法的應(yīng)用[J].高原氣象,2005,24(6):956-965

[8]金鑫,郝振純,張金良.水文模型研究進(jìn)展及發(fā)展方向[J].水土保持研究,2006,13(4):197-202.

[9]崔春光, 彭濤, 殷志遠(yuǎn),等. 暴雨洪澇預(yù)報(bào)研究的若干進(jìn)展[J].氣象科技進(jìn)展, 2011,1(2):32-37.

[10]劉志雨, 謝正輝. TOPKAPI模型的改進(jìn)及其在淮河流域洪水模擬中的應(yīng)用研究[J].水文, 2003,23(6):1-7.

[11]石嵐,馮震,徐麗娜,等.SWAT模型在黃河河萬區(qū)間入庫徑流模擬中的應(yīng)用[J].高原氣象,2012,31(5):1446-1453.

[12]方國華,金德鋼,劉俊,等.甬江流域河網(wǎng)水利計(jì)算模型研究及應(yīng)用[J].災(zāi)害學(xué), 2009,24(2):41-45.

[13]董先勇,樊明蘭,繆韌. 流域水文模型在川中丘陵區(qū)小流域水沙過程研究中的應(yīng)用[J].水土保持研究, 2007,14(3):300-305.

[14]Zeng Xinmin, Zhao Ming,Su Bingkai,et al. Simulation of hydrological model as coupled to a regional climate model[J]. Advances in Atmospheric Sciences,2003,20(2):227-236.

[15]湯成友,項(xiàng)祖?zhèn)?繆韌,等.水箱模型在大尺度流域?qū)崟r(shí)洪水預(yù)報(bào)模型研制中的應(yīng)用[J].水文,2007,27(5):36-51

[16]劉惠民,鄧慧平,孫菽芬,等.陸面模式SSIB耦合TOMMODEL對(duì)流域水文模擬影響的數(shù)值試驗(yàn)[J].高原氣象,2013,32(3):829-838.

[17]許繼軍,楊大文,蔡治國,等. 基于分布式水文模擬的三峽區(qū)間洪水預(yù)報(bào)—模型構(gòu)建及驗(yàn)證[J].水文,2008,28(1):32-37.

[18]劉志雨,侯愛中,王秀慶. 基于分布式水文模型的中小河流洪水預(yù)報(bào)技術(shù)[J].水文,2015,35(1):1-6.

[19]王曉東,董煒華,褚麗娟, 等. 吉林西部旱澇災(zāi)害與降水對(duì)應(yīng)關(guān)系的分析[J].災(zāi)害學(xué),2015,30(4):53-58.

[20]肖瑛,李振興,董玉華.偏差遞歸神經(jīng)網(wǎng)絡(luò)在區(qū)間洪水預(yù)報(bào)中的應(yīng)用研究[J].洛陽大學(xué)學(xué)報(bào),2007,22(2):68-71.

[21]金琪,王繼竹,高琦.基于致洪暴雨預(yù)報(bào)的三峽水庫洪水預(yù)估系統(tǒng)及試驗(yàn)評(píng)估[J].災(zāi)害學(xué),2012,27(3):54-58.

[22]錢鏡林, 張曄,劉國華.基于模糊聚類算法的洪水預(yù)報(bào)研究[J].自然災(zāi)害學(xué)報(bào),2010,19(4):129-134.

[23]李鴻雁, 劉寒冰, 苑希民,等.提高人工神經(jīng)網(wǎng)絡(luò)洪水預(yù)報(bào)精度的研究[J]. 自然災(zāi)害學(xué)報(bào), 2002, 11(1):57-61.

[24]李常斌,李文艷,楊文瑾,等.黃土高原祖厲河流域日降雨徑流過程模擬[J].蘭州大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,48(4):20-26.

[25]雒文生, 宋星原. 洪水預(yù)報(bào)與調(diào)度[M]. 武漢: 湖北科學(xué)技術(shù)出版社, 2000: 19-61.

[26]柳林,安會(huì)靜.中小河流洪水預(yù)報(bào)的難點(diǎn)與解決方案探討[J]. 海河水利, 2013(5):51-53.

[27]楊紅娟,叢振濤,趙巖,等.葉爾羌河流域綠洲蒸散量的遙感估算[J].干旱區(qū)研究,2012，29(3):479-486.

[28]范建忠,李登科,高茂盛. 基于MOD16的陜西省蒸散量時(shí)空分布特征[J].生態(tài)環(huán)境學(xué)報(bào),2014,23(9): 1536-1543.

*收稿日期：2015-12-28修回日期：2016-02-23

作者簡介：李明(1974-)，男，漢族，新疆昌吉人，高級(jí)工程師，從事氣象災(zāi)害預(yù)警研究.E-mail：lm0223@sina.cn

中圖分類號(hào)：X43

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1000-811X(2016)03-0054-07

doi:10.3969/j.issn.1000-811X.2016.03.009

Statistical Methods of Meteorological Forecasting for Small-and Medium-sized Rivers Flood and Their Applications in South Part of Shaanxi Province

LI Ming

(ShaanxiMeteorologicalObservatory,Xi’an710015,China)

Abstract:Based on determining standard of the small- and medium-sized rivers and its flood levels, It is put forward that a statistical forecasting method for small- and medium-sized rivers flood can meet the operational demand in this paper. Runoff depth is calculated on the base of area rainfall of river basin from radar quantity precipitation estimation and previous influence rainfall. Comprehensively considering the water level of the last day, runoff depth and water level of the same day, the binaryregression equation aimed at water level prediction is constructed by using the data of daily water level and the rainfall amount of automatic meteorological station, radar echo over river basin from June to September during 2013-2014.The example of Ziwu river and the other rivers are operated during the flood season in 2015, the results show that the flood forecasting equations constructed by this method can forecast accurately the class of meteorological possibility for the small- and medium-sized river flood, and it is very simple.

Key words:south part of Shaanxi province; small- and medium-sized river; flood; meteorological forecast; binary equation

李明.陜南中小河流洪水氣象統(tǒng)計(jì)預(yù)報(bào)模型業(yè)務(wù)化研究與應(yīng)用[J].災(zāi)害學(xué), 2016,31(3):54-59,65.[LI Ming. Statistical Methods of Meteorological Forecasting for Small-and Medium-sized Rivers Flood and Their Applications in South Part of Shaanxi Province[J].Journal of Catastrophology, 2016,31(3)：54-59,65.]

金項(xiàng)目：公益性行業(yè)(氣象)科研專項(xiàng)(GYHY201206029)；中國氣象局2014年氣象關(guān)鍵技術(shù)集成與應(yīng)用項(xiàng)目(CMAGJ2014M52)；陜西省氣象局2015年重點(diǎn)科研項(xiàng)目(2015Z-2)