趙 放 王東法 楊 軍 陳 列 徐月飛

(1．浙江省氣象臺(tái)，浙江 杭州310017;2．臺(tái)州市氣象局，浙江椒江331002;3．衢州市氣象局，浙江 衢州324000)

0 引言

地質(zhì)災(zāi)害調(diào)查表明，滑坡在我國(guó)地質(zhì)災(zāi)害中所占比例高達(dá)51% ，而暴雨誘發(fā)的滑坡比重竟占滑坡總數(shù)的90%［1］。浙江省發(fā)生的滑坡、泥石流等斜坡地質(zhì)災(zāi)害以突發(fā)暴雨型為主。1996年8月1日，受9608號(hào)臺(tái)風(fēng)影響，青田尾礦庫(kù)大壩沖毀，形成巨大的泥石流災(zāi)害，造成近百人傷亡，并沖毀廠房、民房，損毀大量農(nóng)田，直接經(jīng)濟(jì)損失達(dá)3000萬(wàn)元。2005年9月3日19：00—23：00，受臺(tái)風(fēng)后部氣流影響臨安市西部昌化地區(qū)遭遇特大強(qiáng)降雨，引發(fā)多處泥石流災(zāi)害，沖毀民房132間，死亡9人。2009年8月13日23：00，暴雨引發(fā)清涼峰鎮(zhèn)發(fā)生重大山體滑坡，造成一間3層民房被掩埋，11人死亡2人失蹤。因此，重視和加強(qiáng)由暴雨等氣象因素誘發(fā)的突發(fā)性地質(zhì)災(zāi)害的研究意義重大。造成地質(zhì)災(zāi)害的因素有很多，持續(xù)或急劇的降水是最主要的因素之一。目前，尚未研究出對(duì)降雨誘發(fā)滑坡泥石流進(jìn)行定量預(yù)測(cè)預(yù)報(bào)的成熟方法。另外，加強(qiáng)災(zāi)害預(yù)警信息產(chǎn)品開(kāi)發(fā)并使之與城市網(wǎng)格化的管理理念相一致，也是突發(fā)公共事件預(yù)警處理網(wǎng)格化管理的迫切需求。

以氣象方法進(jìn)行地質(zhì)災(zāi)害的預(yù)報(bào)國(guó)內(nèi)外作了不少研究［2］，地質(zhì)災(zāi)害氣象等級(jí)預(yù)報(bào)是當(dāng)前常用的一種方法，以短期預(yù)報(bào)為主。主要是采用基于統(tǒng)計(jì)資料的臨界雨量閾值判別法，來(lái)獲取滑坡、泥石流等發(fā)生的指標(biāo)。如殷坤龍等［3］指出降雨量的閾值受地區(qū)、滑坡類型等影響十分明顯，滑坡發(fā)生的概率和滑坡數(shù)量不僅與降雨量的大小成正相關(guān)關(guān)系，而且與滑坡發(fā)生的當(dāng)天降雨及前期降雨特征關(guān)系密切。根據(jù)單九生等［4］對(duì)誘發(fā)江西滑坡的降水特征的統(tǒng)計(jì)分析，王錫穩(wěn)等［5］對(duì)隴南“5·31”特大泥石流災(zāi)害的成因分析，馬力等［6］對(duì)重慶市山體滑坡氣象條件等級(jí)預(yù)報(bào)的研究，以及杜惠良等［7］針對(duì)滑坡、泥石流應(yīng)用MM5中尺度預(yù)報(bào)模式的應(yīng)用等研究，都表明暴雨或持續(xù)性降雨是誘發(fā)山體滑坡的最主要因素，并與連續(xù)降水的累計(jì)值有關(guān)。然而，由于形成滑坡災(zāi)害的因素很多，地質(zhì)背景條件差異很大，它和降雨的關(guān)系十分復(fù)雜，傳統(tǒng)的預(yù)報(bào)方法仍然帶有很強(qiáng)的主觀隨意性和非客觀化的分析和預(yù)測(cè)因素。一是由于泥石流災(zāi)害的發(fā)生在時(shí)空上存在明顯的不確定性，預(yù)報(bào)因素包含的地質(zhì)背景條件及其機(jī)理的信息相對(duì)不足;二是常規(guī)降水的定量化預(yù)報(bào)、降水落區(qū)的精度和定位等仍有很大的偏差;三是由于地形復(fù)雜、強(qiáng)降水分布不均，用空間尺度較大的資料作滑坡泥石流災(zāi)害的監(jiān)測(cè)和預(yù)警還存在著許多困難和不足。隨著現(xiàn)代化氣象裝備的發(fā)展，加強(qiáng)以雷達(dá)、衛(wèi)星、自動(dòng)站等新的監(jiān)測(cè)手段的應(yīng)用，并結(jié)合數(shù)值預(yù)報(bào)、高分辨地理信息等技術(shù)，在傳統(tǒng)的氣象地質(zhì)災(zāi)害預(yù)測(cè)方法的基礎(chǔ)上，發(fā)展和完善該領(lǐng)域的研究與業(yè)務(wù)，建立時(shí)空更加精細(xì)的突發(fā)性的氣象災(zāi)害預(yù)警方法，具有十分重要的理論意義與實(shí)用價(jià)值。

1 主要技術(shù)方法

新一代多普勒天氣雷達(dá)可以獲得高時(shí)空分辨率的降水強(qiáng)度信息，是監(jiān)測(cè)臺(tái)風(fēng)、強(qiáng)對(duì)流等天氣，并進(jìn)行臨近預(yù)報(bào)的重要手段之一。利用雷達(dá)數(shù)字組網(wǎng)、自動(dòng)站等監(jiān)測(cè)手段，結(jié)合數(shù)值預(yù)報(bào)產(chǎn)品、地理信息等技術(shù)，形成連續(xù)滾動(dòng)的高分辨(0．01經(jīng)緯度格距)0～3 h定量降水預(yù)報(bào)，通過(guò)數(shù)字高程模型、遙感數(shù)據(jù)(地質(zhì)圖、土壤圖、土地利用現(xiàn)狀)獲取浙江省境內(nèi)地質(zhì)災(zāi)害易發(fā)區(qū)數(shù)據(jù)分布，在Logistic回歸模型基礎(chǔ)上對(duì)傳統(tǒng)的災(zāi)害預(yù)警模型進(jìn)行改進(jìn)，以動(dòng)態(tài)方式逐時(shí)遞進(jìn)的引入前期有效累積降雨量，以網(wǎng)格化的處理方式，建立具有高分辨(約1 km×1 km)的強(qiáng)降雨誘發(fā)型滑坡、泥石流的定點(diǎn)定量臨近預(yù)警客觀分析方法。

1．1 浙江突發(fā)地質(zhì)災(zāi)害特點(diǎn)及易發(fā)區(qū)分布

浙江省地形復(fù)雜，島嶼眾多，地勢(shì)由西南向東北呈階梯狀傾斜。浙西南山地多千米以上山嶺，最高山峰黃茅尖海拔1929 m;中部多500 m左右的山地丘陵;東北部平原區(qū)屬長(zhǎng)江三角洲一部分。當(dāng)前，斜坡地質(zhì)災(zāi)害是浙江省危害最嚴(yán)重的地質(zhì)災(zāi)害災(zāi)種，它主要發(fā)生于山地小流域中。山地斜坡地質(zhì)災(zāi)害是水流和巖土地質(zhì)體相互作用下致使巖土體發(fā)生變形破壞、滑移、墜落或流滾而釀成的災(zāi)害，崩塌、滑坡、泥石流的比例各占17%，74%，9%［1］。根據(jù)近年來(lái)的災(zāi)害統(tǒng)計(jì)數(shù)據(jù)和浙江地質(zhì)背景條件及其機(jī)理的信息［8］，按照地質(zhì)災(zāi)害預(yù)報(bào)5個(gè)等級(jí)的劃分標(biāo)準(zhǔn)(1級(jí)為可能性很小，2級(jí)為可能性較小，3級(jí)為可能性較大(注意級(jí))，4級(jí)為可能性大(預(yù)警級(jí))，5級(jí)為可能性很大(警報(bào)級(jí)))，將滑坡(泥石流)災(zāi)害分布特點(diǎn)及對(duì)觸發(fā)因素的分析(地形(坡度和形態(tài)等))、斷裂構(gòu)造及地層巖性信息，土壤濕度、土壤含水量、植被/裸露地作為輸入變量，把當(dāng)前已知的并有具體發(fā)生日期和位置記載的滑坡災(zāi)害點(diǎn)作為輸出變量，在GIS平臺(tái)上從數(shù)字高程模型(DEM)、遙感數(shù)據(jù)、地質(zhì)圖、土壤圖、土地利用現(xiàn)狀圖數(shù)據(jù)庫(kù)進(jìn)行調(diào)用分析，形成浙江省境內(nèi)高分辨格點(diǎn)(1 km×1 km)的地質(zhì)災(zāi)害易發(fā)區(qū)分布數(shù)據(jù)(圖1)，并將其作為影響該格點(diǎn)上地質(zhì)災(zāi)害發(fā)生概率大小的因子之一。

圖1 浙江地質(zhì)災(zāi)害易發(fā)區(qū)分布(分辨率1 km×1 km)

1．2 資料處理與網(wǎng)格選取

1．2．1 雷達(dá)動(dòng)態(tài)降水反演

降水因子至關(guān)重要，主要依靠雷達(dá)、自動(dòng)站等資料分析獲得。根據(jù)浙江地理情況，選25°N～33°N，116°E ～124°E 所圍區(qū)域，進(jìn)行資料時(shí)空同化和網(wǎng)格量化。網(wǎng)格點(diǎn)水平方向的經(jīng)緯度分辨率為0．01°×0．01°，通過(guò) GIS 平臺(tái)與地質(zhì)災(zāi)害易發(fā)區(qū)分布數(shù)據(jù)相結(jié)合，形成高密度的預(yù)警網(wǎng)格單元元素(約1 km ×1 km)。雷達(dá)時(shí)空和網(wǎng)格量化處理，以某一部雷達(dá)體掃時(shí)間所獲得的基數(shù)據(jù)為基準(zhǔn)，截取其余雷達(dá)在6 min時(shí)間段內(nèi)所生成的數(shù)據(jù)，集成為區(qū)域范圍的雷達(dá)體掃數(shù)據(jù)(理論上時(shí)間誤差小于3 min)。根據(jù)雷達(dá)的仰角、方位和斜距在雷達(dá)球坐標(biāo)系中的位置，用內(nèi)插方法給該網(wǎng)格點(diǎn)回波強(qiáng)度賦值，得到該網(wǎng)格點(diǎn)上的分析值。取笛卡爾坐標(biāo)網(wǎng)格在垂直方向分辨率為1 km，水平方向的經(jīng)緯度分辨率為0．01°×0．01°。多個(gè)雷達(dá)的資料重疊區(qū)網(wǎng)格單元i的反射率值的公式為：

雷達(dá)降水反演，在聚類分析的基礎(chǔ)上使用改進(jìn)的最佳窗概率配對(duì)法計(jì)算Z-I關(guān)系中的系數(shù)A和b，實(shí)時(shí)計(jì)算基本反射率因子Z和自動(dòng)站降水量I的動(dòng)態(tài)關(guān)系。采用變分技術(shù)對(duì)估算的小時(shí)降水強(qiáng)度進(jìn)行校準(zhǔn)，并進(jìn)行誤差定量分析。該方法將雨量計(jì)單點(diǎn)測(cè)量精度較高和雷達(dá)能測(cè)量降水時(shí)空分布的優(yōu)點(diǎn)結(jié)合起來(lái)，可獲得比單純使用雷達(dá)作定量降水估測(cè)更為精確的結(jié)果［9］。

伴隨動(dòng)態(tài)定量估測(cè)降水精確度分析，用評(píng)估因子(E)隨時(shí)進(jìn)行反饋訂正：

評(píng)估因子(E)：E=QR/QG

式中，QR—雷達(dá)估測(cè)雨量(插值到站點(diǎn))，QG—對(duì)應(yīng)自動(dòng)站雨量。

如以 2007年 10月 07日 11：00—12：00(UTC)降水過(guò)程雷達(dá)定量估測(cè)降水精確度評(píng)估因子(E)為例(圖2)。該降水過(guò)程(全省自動(dòng)站≥1．0 mm降水743個(gè))的評(píng)估因子檢驗(yàn)結(jié)果分別為 0．98(≥1．0 mm)和 0．95(≥10．0 mm)。

圖2 2007年10月07日12：00(UTC)降水過(guò)程定量估測(cè)降水精確度評(píng)估因子(E)

1．2．2 3 h 定量降水預(yù)報(bào)

臨近預(yù)報(bào)正在向自動(dòng)化、快速化、高時(shí)空分辨率方向發(fā)展。近年來(lái)，通過(guò)采用和引進(jìn)國(guó)內(nèi)外先進(jìn)的技術(shù)方法和手段，浙江省的定量降水預(yù)報(bào)系統(tǒng)得以不斷改進(jìn)，主要的方法是：1 h的降水臨近預(yù)報(bào)，用改進(jìn)的交叉相關(guān)法識(shí)別和矩心跟蹤法進(jìn)行回波團(tuán)跟蹤外推，再根據(jù)最新時(shí)刻的雨量數(shù)據(jù)，作1 h降水估計(jì)線性外推;2～3h的降水預(yù)報(bào)：考慮到降水率變化及簡(jiǎn)單外推會(huì)造成較大的累計(jì)誤差，采用數(shù)值模式的輸出結(jié)果對(duì)外推進(jìn)行匹配和適當(dāng)控制。首先，以1 h降水預(yù)報(bào)與模式預(yù)報(bào)對(duì)應(yīng)的網(wǎng)格點(diǎn)的預(yù)報(bào)進(jìn)行陣列式匹配分析，調(diào)整模式在降水落區(qū)上的時(shí)空差;其次，將TREC風(fēng)場(chǎng)與相應(yīng)時(shí)刻各高度層模式輸出的引導(dǎo)氣流進(jìn)行比較，搜索出最具相關(guān)的引導(dǎo)氣流的空間分布，并以此引導(dǎo)場(chǎng)對(duì)1 h回波作2～3 h的外推，同時(shí)以模式物理量的輸出結(jié)果對(duì)雷達(dá)外推回波團(tuán)進(jìn)行強(qiáng)弱變化(生消)控制，進(jìn)而獲得2～3 h降水預(yù)報(bào)?；夭ㄍ馔浦?，考慮到浙江經(jīng)常受到臺(tái)風(fēng)、東向暴雨等特有的天氣系統(tǒng)的影響，針對(duì)螺旋帶結(jié)構(gòu)中尺度對(duì)流群既具有系統(tǒng)性的移向趨勢(shì)又有旋轉(zhuǎn)性的運(yùn)動(dòng)特點(diǎn)，增加了半拉格朗日訂正的算法。此外，外推中引入方向變量因子，通過(guò)判斷回波連續(xù)兩個(gè)時(shí)段內(nèi)移動(dòng)方向上的變化，在外推時(shí)加上旋轉(zhuǎn)修正，從而提高了帶狀結(jié)構(gòu)回波移動(dòng)的預(yù)報(bào)時(shí)效。經(jīng)檢驗(yàn)，3 h 強(qiáng)回波外推預(yù)報(bào)相關(guān)系數(shù)在0．4 ～0．85［10，11］。實(shí)際應(yīng)用表明，上述方法極大地提高了3h定量降水的客觀預(yù)報(bào)能力，圖3為2007年10月7日21：00 1 h和3 h的降水預(yù)報(bào)與實(shí)況降水對(duì)比，降水預(yù)報(bào)與實(shí)況的疊合比較一致。

1．3 滑坡泥石流預(yù)警模型

地質(zhì)災(zāi)害的預(yù)報(bào)實(shí)質(zhì)是對(duì)觸發(fā)條件的預(yù)報(bào)，從統(tǒng)計(jì)學(xué)角度出發(fā)，對(duì)其產(chǎn)生影響的各因子數(shù)據(jù)可以作為自變量。而災(zāi)害的發(fā)生與不發(fā)生可以作為分類因變量。由于不是連續(xù)變量，線性回歸將不適用于此類自變量和因變量之間的推導(dǎo)關(guān)系，通常采用對(duì)數(shù)線性模型。

Logistic回歸模型是對(duì)數(shù)線性模型的一種特殊形式［12］，設(shè)P為某事發(fā)生的概率，取值范圍為［0，1］，xm為自變量，建立回歸方程：

α 為常數(shù);βi(i=1，2，…，m) 為邏輯回歸系數(shù)。

圖3 2007年10月7日21：00短時(shí)臨近預(yù)報(bào)(1，3 h)與實(shí)況對(duì)比

強(qiáng)降水與地質(zhì)災(zāi)害等級(jí)關(guān)系，通常以“前期有效降雨量”的概念，取自當(dāng)天及前幾天的逐日降雨量用于滑坡災(zāi)害分析。所謂“前期有效降雨量”是指前期降雨進(jìn)入巖土體并一直滯留，直至發(fā)災(zāi)當(dāng)天的雨量。國(guó)內(nèi)外學(xué)者對(duì)此已作過(guò)相應(yīng)的研究，并提出了計(jì)算進(jìn)入巖土體降雨量的經(jīng)驗(yàn)公式：

式中：ra，當(dāng)天的前期有效降雨量;k，有效雨量系數(shù)(一般取0．84);rn，前第n天的降雨量;

由此得出，當(dāng)日某時(shí)刻未來(lái)3 h后激發(fā)雨量R：

R=Y3(3 h雨量預(yù)報(bào))+R21(21 h實(shí)際累計(jì)降水量)。

各網(wǎng)格點(diǎn)上，某時(shí)刻未來(lái)3 h后實(shí)際有效累積降雨量 Y(i，j)：

通過(guò)統(tǒng)計(jì)浙江省1950年以來(lái)地質(zhì)災(zāi)害的歷史資料，確定前期有效累積降雨量和當(dāng)日降雨量引發(fā)地質(zhì)災(zāi)害的最小降雨量界限(該界限之下未發(fā)生災(zāi)害)和最大降雨量界限(該界限之上必定發(fā)生災(zāi)害)，從而得到災(zāi)害概率發(fā)生的回歸方程和預(yù)警指數(shù)：

W(i，j)為預(yù)警指數(shù)，N(i，j)為地質(zhì)災(zāi)害易發(fā)分級(jí)指數(shù)，P為發(fā)生災(zāi)害概率，A，B為權(quán)重系數(shù)(對(duì)新出現(xiàn)的災(zāi)害點(diǎn)通過(guò)調(diào)整該點(diǎn)的權(quán)重系數(shù)也可使預(yù)報(bào)模型適應(yīng)其改變)。

依據(jù)上面的公式即可得到各個(gè)網(wǎng)格點(diǎn)上發(fā)生災(zāi)害概率的界限和地質(zhì)災(zāi)害臨近預(yù)報(bào)預(yù)警指數(shù)W(i，j)，根據(jù)實(shí)時(shí)監(jiān)測(cè)覆蓋全省地區(qū)水平分辨率為0．01°×0．01°經(jīng)緯度格點(diǎn)的地質(zhì)災(zāi)害臨近預(yù)報(bào)預(yù)警指數(shù)W(i，j)，實(shí)現(xiàn)定量化分析和氣象災(zāi)害等級(jí)判斷。若預(yù)測(cè)未來(lái)3 h達(dá)到臨界降雨量及出現(xiàn)較高的預(yù)警，則可發(fā)布格點(diǎn)所對(duì)應(yīng)的地理位置3 h時(shí)效的地質(zhì)災(zāi)害的預(yù)警。

對(duì)2005—2008年期間浙江發(fā)生的降雨誘發(fā)的109個(gè)滑坡泥石流災(zāi)害調(diào)查實(shí)例，應(yīng)用該預(yù)報(bào)預(yù)警模型作驗(yàn)證，在災(zāi)害發(fā)生的地點(diǎn)均給出了高級(jí)別(3級(jí)以上)的預(yù)警指數(shù)，無(wú)漏報(bào)現(xiàn)象，3級(jí)以下的預(yù)警指數(shù)存在一定范圍的虛報(bào)(也和災(zāi)情的收集不完整有關(guān))，但從實(shí)用的角度看仍然可以起到較好的警示意義。事實(shí)上，在一般情況下由暴雨所引發(fā)的地質(zhì)災(zāi)害，從泥石流發(fā)生到成災(zāi)過(guò)程之間還有一定的時(shí)間，最大雨強(qiáng)的出現(xiàn)到泥石流暴發(fā)有0．5～2 h的時(shí)間間隔，如避難信息能在這段時(shí)間間隔內(nèi)發(fā)布，就有可能提前3～5 h采取措施，避免重大的人員傷亡和經(jīng)濟(jì)損失。

2 個(gè)例應(yīng)用

2009年8月13日下午，浙江省臨安清涼峰鎮(zhèn)出現(xiàn)暴雨天氣，16：00—18：00雨量達(dá)到88 mm。23：30，暴雨引發(fā)該鎮(zhèn)林竹村梓棚塢發(fā)生重大山體滑坡，造成了嚴(yán)重災(zāi)害。圖4所示為8月13日15：00 3 h實(shí)效降雨量預(yù)報(bào)(圖4a)和網(wǎng)格點(diǎn)上地質(zhì)災(zāi)害預(yù)警指數(shù)(圖4b)。經(jīng)檢驗(yàn)，預(yù)報(bào)預(yù)警模型在災(zāi)害發(fā)生的地點(diǎn)給出了高級(jí)別的預(yù)警，預(yù)報(bào)與實(shí)況基本吻合。

圖4 2009年8月13日15：00地質(zhì)災(zāi)害預(yù)報(bào)產(chǎn)品(黑三角符號(hào)上方格點(diǎn)為災(zāi)害實(shí)況位置)

2010年受熱帶風(fēng)暴“莫蘭蒂”以及其減弱后的低壓環(huán)流與冷空氣相結(jié)合的影響，浙江全省大部分地區(qū)普降中到大雨 。9月11日，杭州市東部早晨(05：00—07：00)出現(xiàn)強(qiáng)降雨天氣過(guò)程，主要出現(xiàn)在主城區(qū)、蕭山和富陽(yáng)東部，全市累計(jì)面雨量35．2 mm，其中主城區(qū)81．5 mm、蕭山80．9 mm、富陽(yáng)30．9 mm、桐廬19．3 mm、余杭15．3 mm。杭州市蜀山街道金西村、聞堰鎮(zhèn)、湘湖旅游度假區(qū)、杭州市濱江區(qū)長(zhǎng)河街道湯家井村美女山公墓等16個(gè)地點(diǎn)，相繼發(fā)生了泥石流、滑坡等小面積地質(zhì)災(zāi)害。預(yù)報(bào)預(yù)警模型業(yè)務(wù)試運(yùn)行系統(tǒng)，從05：00起對(duì)災(zāi)害發(fā)生的地點(diǎn)連續(xù)給出了高級(jí)別的預(yù)警(圖略)，為氣象服務(wù)及時(shí)提供了準(zhǔn)確的預(yù)報(bào)預(yù)警信息。

3 結(jié)語(yǔ)

(1)利用雷達(dá)、自動(dòng)站等新的監(jiān)測(cè)手段，并結(jié)合數(shù)值預(yù)報(bào)產(chǎn)品和地理信息等技術(shù)并結(jié)合地質(zhì)背景條件，通過(guò)實(shí)現(xiàn)降水及實(shí)效降雨量的格點(diǎn)化預(yù)報(bào)，從時(shí)間和空間尺度上加密了對(duì)滑坡泥石流災(zāi)害的監(jiān)測(cè)和預(yù)警，使精細(xì)化到鄉(xiāng)、鎮(zhèn)的泥石流滑坡災(zāi)害預(yù)警預(yù)報(bào)成為可能。

(2)最大雨強(qiáng)出現(xiàn)到泥石流暴發(fā)一般有0．5～2 h的時(shí)間間隔，從泥石流發(fā)生到成災(zāi)過(guò)程之間還有一定的時(shí)間。如避難信息能在這段時(shí)間間隔內(nèi)發(fā)布，就有可能提前3～5 h采取措施，避免重大的人員傷亡和經(jīng)濟(jì)損失。另外，將數(shù)值預(yù)報(bào)結(jié)果與預(yù)警模型連接，即可獲得更長(zhǎng)時(shí)效及相應(yīng)分辨率格點(diǎn)的災(zāi)害預(yù)報(bào)。

(3)雖然根據(jù)歷史資料發(fā)生災(zāi)害概率的界限，實(shí)現(xiàn)了格點(diǎn)化的定量化分析，但地質(zhì)災(zāi)害的發(fā)生存在著許多不確定的因素，滑坡泥石流預(yù)警模型需要在今后的實(shí)際應(yīng)用中不斷檢驗(yàn)和優(yōu)化。

［1］ 葉 勇．浙江省小流域山洪災(zāi)害與預(yù)警技術(shù)研究［D］．浙江大學(xué)碩士學(xué)位論文．2008，5．

［2］ 周創(chuàng)兵，李典慶．暴雨誘發(fā)滑坡致災(zāi)機(jī)理與減災(zāi)方法研究進(jìn)展［J］．地球科學(xué)進(jìn)展．2009，24(5)：477-487．

［3］ 殷坤龍，張桂榮，龔日祥，等．基于Web2GIS的浙江省地質(zhì)災(zāi)害實(shí)時(shí)預(yù)警預(yù)報(bào)系統(tǒng)設(shè)計(jì)［J］．水文地質(zhì)工程地質(zhì)，2003，30(3)：19-22．

［4］ 單九生，魏 麗，劉修奮，等．誘發(fā)江西滑坡的降水特征分析［J］．氣象，2004，30(1)：13-15．

［5］ 王錫穩(wěn)，陶建紅，馮 軍，等．隴南“5·31”特大泥石流災(zāi)害成因分析［J］．氣象，2004，30(10)：43-46．

［6］ 馬 力，廖代強(qiáng)．重慶市山體滑坡氣象條件等級(jí)預(yù)報(bào)業(yè)務(wù)系統(tǒng)［J］．應(yīng)用氣象學(xué)報(bào)，2003，14(1)：122-124．

［7］ 杜惠良，鈕學(xué)新，殷坤龍．浙江省滑波、泥石流氣象條件分析及其預(yù)報(bào)研究［J］．熱帶氣象學(xué)報(bào)，2005，21(6)：642-650．

［8］ 陳海燕，潘小凡，吳利紅．浙江泥石流氣象特征分析［J］．災(zāi)害學(xué)，2005，20(1)：61-64．

［9］ 冀春曉，陳聯(lián)壽，徐祥德，等．多普勒雷達(dá)資料動(dòng)態(tài)定量估測(cè)臺(tái)風(fēng)小時(shí)降水量的研究［J］．熱帶氣象學(xué)報(bào)，2008，24(2)：147-155．

［10］ 趙 放，冀春曉，任鴻翔，等．臺(tái)風(fēng)螺旋帶回波識(shí)別追蹤及結(jié)構(gòu)分析技術(shù)應(yīng)用［J］．浙江氣象，2008，29(1)：15-21．

［11］ 趙 放，冀春曉，任鴻翔，等．應(yīng)用多普勒雷達(dá)制作近海臺(tái)風(fēng)臨近預(yù)報(bào)技術(shù)研究［J］．氣象，2008，34(5)：64-74．

［12］ 李鐵鋒，叢威青．基于Logistic回歸及前期有效雨量的降雨誘發(fā)型滑坡預(yù)測(cè)方法［J］．中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào)，2006，17(1)：33-35．