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

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

0 引言

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

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

1 主要技術方法

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

1．1 浙江突發(fā)地質災害特點及易發(fā)區(qū)分布

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

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

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

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

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

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

伴隨動態(tài)定量估測降水精確度分析，用評估因子(E)隨時進行反饋訂正：

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

式中，QR—雷達估測雨量(插值到站點)，QG—對應自動站雨量。

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

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

1．2．2 3 h 定量降水預報

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

1．3 滑坡泥石流預警模型

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

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

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

圖3 2007年10月7日21：00短時臨近預報(1，3 h)與實況對比

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

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

由此得出，當日某時刻未來3 h后激發(fā)雨量R：

R=Y3(3 h雨量預報)+R21(21 h實際累計降水量)。

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

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

W(i，j)為預警指數(shù)，N(i，j)為地質災害易發(fā)分級指數(shù)，P為發(fā)生災害概率，A，B為權重系數(shù)(對新出現(xiàn)的災害點通過調整該點的權重系數(shù)也可使預報模型適應其改變)。

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

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

2 個例應用

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

圖4 2009年8月13日15：00地質災害預報產(chǎn)品(黑三角符號上方格點為災害實況位置)

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

3 結語

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

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

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

［1］ 葉 勇．浙江省小流域山洪災害與預警技術研究［D］．浙江大學碩士學位論文．2008，5．

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

［3］ 殷坤龍，張桂榮，龔日祥，等．基于Web2GIS的浙江省地質災害實時預警預報系統(tǒng)設計［J］．水文地質工程地質，2003，30(3)：19-22．

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

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

［6］ 馬 力，廖代強．重慶市山體滑坡氣象條件等級預報業(yè)務系統(tǒng)［J］．應用氣象學報，2003，14(1)：122-124．

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

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

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

［10］ 趙 放，冀春曉，任鴻翔，等．臺風螺旋帶回波識別追蹤及結構分析技術應用［J］．浙江氣象，2008，29(1)：15-21．

［11］ 趙 放，冀春曉，任鴻翔，等．應用多普勒雷達制作近海臺風臨近預報技術研究［J］．氣象，2008，34(5)：64-74．

［12］ 李鐵鋒，叢威青．基于Logistic回歸及前期有效雨量的降雨誘發(fā)型滑坡預測方法［J］．中國地質災害與防治學報，2006，17(1)：33-35．