于竹娟，羅貴東

(達(dá)州市氣象局，四川 達(dá)州635000)

達(dá)州市地處大巴山南麓，是四川省的人口大市、工農(nóng)業(yè)重鎮(zhèn)，是四川省、重慶市和陜西省三省市結(jié)合部的商貿(mào)中心和交通樞紐，素有“川東明珠”之美譽(yù).但是市域地質(zhì)環(huán)境條件復(fù)雜，降雨充沛，人類工程活動(dòng)強(qiáng)烈，因此地質(zhì)災(zāi)害災(zāi)情重、種類全、頻率高，是四川省地質(zhì)災(zāi)害易發(fā)區(qū)之一.

地質(zhì)災(zāi)害是指包括自然因素或者人為活動(dòng)引發(fā)的危害人民財(cái)產(chǎn)和生命安全與地質(zhì)作用有關(guān)的災(zāi)害.達(dá)州市地質(zhì)災(zāi)害以滑坡、崩塌、不穩(wěn)定斜坡、地面塌陷為主.[1]達(dá)州地質(zhì)災(zāi)害多發(fā)易發(fā)，除了與本身地質(zhì)構(gòu)造特點(diǎn)有關(guān)外，還與降水、人類工程活動(dòng)等自然條件有關(guān).[2-11]大量研究表明，[12-17]地質(zhì)災(zāi)害與降水的時(shí)空分布有著密切聯(lián)系，降雨強(qiáng)度越大，地質(zhì)災(zāi)害越容易發(fā)生.

2017年我們采用基于前期降水量和當(dāng)日降水量的降水閾值方法，[18]研發(fā)了基于概率分布的達(dá)州市地質(zhì)災(zāi)害與降水關(guān)系的地質(zhì)災(zāi)害預(yù)報(bào)預(yù)警平臺(tái).前期降水通過氣象站實(shí)際觀測(cè)得到，未來降水則使用了通過我們本地訂正的四川省氣象臺(tái)精細(xì)化平臺(tái)的格點(diǎn)定量降水預(yù)報(bào).四川省氣象臺(tái)2017年發(fā)布的水平分辨率10KM的格點(diǎn)化定量降水預(yù)報(bào)產(chǎn)品，通過本地預(yù)報(bào)員人工訂正后具有較高的分辨率和準(zhǔn)確性.現(xiàn)在我們對(duì)這個(gè)產(chǎn)品的實(shí)際應(yīng)用進(jìn)行定量研究，對(duì)進(jìn)一步改進(jìn)和完善兩種產(chǎn)品，提供預(yù)警效果具有很重要的意義.

持續(xù)性強(qiáng)降雨由于持續(xù)時(shí)間長，累計(jì)雨量大，巖土?xí)浞诌_(dá)到飽和，容易引發(fā)大規(guī)模的滑坡、泥石流.2019年“6.20-6.22”持續(xù)性強(qiáng)降雨造成全市9.1萬人不同程度受災(zāi)，直接經(jīng)濟(jì)損失達(dá)3.58億元，其中，引發(fā)萬源、大竹多起地質(zhì)災(zāi)害.針對(duì)本次過程，達(dá)州市自然資源和規(guī)劃局、達(dá)州市氣象局聯(lián)合發(fā)布了三期地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)黃色(Ⅲ級(jí))或橙色(Ⅱ級(jí))預(yù)警.本文嘗試以“6.20-6.22”持續(xù)性強(qiáng)降雨引發(fā)的地質(zhì)災(zāi)害為對(duì)象，定量檢驗(yàn)達(dá)州市地質(zhì)災(zāi)害預(yù)報(bào)預(yù)警平臺(tái)發(fā)布的地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)警產(chǎn)品，總結(jié)經(jīng)驗(yàn)，以期更好地使用地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)模型和格點(diǎn)化定量降水預(yù)報(bào)，提高預(yù)警效果.

1 過程概述

15日開始，我市出現(xiàn)連陰雨天氣，20日到22日全市陸續(xù)出現(xiàn)強(qiáng)降雨.20日萬源北部暴雨；21日大竹、渠縣部分地方出現(xiàn)暴雨；22日宣漢大部地方大到暴雨，其余地方小到中雨，局部暴雨.此次連續(xù)性降雨天氣過程，全市共有53個(gè)站點(diǎn)累積雨量超過100mm，最大大竹永勝302mm.特別是在20到22日的強(qiáng)降雨過程中，有11站次小時(shí)雨強(qiáng)超過50 mm·h-1，萬源竹峪和大竹永勝兩小時(shí)達(dá)大暴雨.20到22日的強(qiáng)降雨過程，一是雨強(qiáng)大，二是時(shí)間集中，造成大竹、萬源等多地出現(xiàn)滑坡、泥石流和崩塌.

2 資料來源與檢驗(yàn)方法

2.1 資料來源及說明

地質(zhì)災(zāi)害災(zāi)情險(xiǎn)情資料來源于達(dá)州市地質(zhì)環(huán)境監(jiān)測(cè)站，包含了地質(zhì)災(zāi)害發(fā)生時(shí)間、地點(diǎn)、規(guī)模、受災(zāi)人口、原因、財(cái)產(chǎn)損失等信息.降水?dāng)?shù)據(jù)，包含國家氣象站和區(qū)域自動(dòng)氣象站觀測(cè)數(shù)據(jù)、達(dá)州市氣象臺(tái)24h格點(diǎn)化定量降水預(yù)報(bào)產(chǎn)品.達(dá)州市地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)警產(chǎn)品預(yù)報(bào)時(shí)效為24h.產(chǎn)品預(yù)報(bào)時(shí)段從當(dāng)日20時(shí)至次日20時(shí).

2.2 預(yù)警檢驗(yàn)方法

預(yù)警檢驗(yàn)方法我們采用了兩種評(píng)價(jià)標(biāo)準(zhǔn).

(1)按照中國地質(zhì)災(zāi)害防治工程行業(yè)協(xié)會(huì)2018年發(fā)布的《地質(zhì)災(zāi)害區(qū)域氣象風(fēng)險(xiǎn)預(yù)警標(biāo)準(zhǔn)(試行)》(下面簡稱《標(biāo)準(zhǔn)》)中關(guān)于預(yù)警準(zhǔn)確率的評(píng)價(jià)方法規(guī)定，如果有地質(zhì)災(zāi)害點(diǎn)落入預(yù)報(bào)區(qū)范圍內(nèi)，則表示此次預(yù)報(bào)準(zhǔn)確，如果無地質(zhì)災(zāi)害點(diǎn)落入預(yù)報(bào)區(qū)范圍內(nèi)，則表示此次預(yù)報(bào)不準(zhǔn)確.將預(yù)報(bào)準(zhǔn)確的次數(shù)除以總的預(yù)報(bào)次數(shù)即為預(yù)報(bào)準(zhǔn)確率，計(jì)算公式如下：

式中：p—預(yù)報(bào)準(zhǔn)確率；m—預(yù)報(bào)準(zhǔn)確的次數(shù)；n—預(yù)報(bào)次數(shù).

(2)2016年10月，中國氣象局下發(fā)了《暴雨誘發(fā)地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)警業(yè)務(wù)規(guī)范》(以下簡稱《規(guī)范》).《規(guī)范》規(guī)定，檢驗(yàn)預(yù)報(bào)準(zhǔn)確率時(shí)以縣級(jí)區(qū)域和單點(diǎn)(滑坡、泥石流)為檢驗(yàn)對(duì)象.當(dāng)以縣級(jí)區(qū)域?yàn)闄z驗(yàn)對(duì)象時(shí)，在該預(yù)報(bào)時(shí)段內(nèi)任意一點(diǎn)出現(xiàn)了地質(zhì)災(zāi)害，則記錄為一次預(yù)報(bào)正確.如果在上述區(qū)域，該預(yù)報(bào)時(shí)段所有的地質(zhì)災(zāi)害點(diǎn)都未出現(xiàn)地質(zhì)災(zāi)害，則記錄為一次空?qǐng)?bào).

具體方案如下：

命中率：

(2)

(3)

(4)

其中，NA為預(yù)警服務(wù)產(chǎn)品發(fā)布正確次數(shù)，NB為空?qǐng)?bào)次數(shù)，NC為漏報(bào)次數(shù)，ND為無預(yù)警服務(wù)產(chǎn)品預(yù)報(bào)正確次數(shù).

3 “6.20-6.22”持續(xù)性強(qiáng)降水及其引發(fā)的地質(zhì)災(zāi)害

3.1 “6.20-6.22”持續(xù)性強(qiáng)降水分析

圖1是2019年6月19日20時(shí)至22日20時(shí)的過程累積雨量分布，我市北部和南部普遍出現(xiàn)了大到暴雨，部分地方降雨量超過了100mm，且均伴有小時(shí)雨強(qiáng)超過20mm的短時(shí)強(qiáng)降水.累積降雨量超過100毫米的地方主要出現(xiàn)在萬源西北部、渠縣南部和大竹部分鄉(xiāng)鎮(zhèn).

圖1 2019年6月19日20時(shí)-22日20時(shí)達(dá)州市累積降雨量

3.2 地質(zhì)災(zāi)害概況

根據(jù)達(dá)州市地質(zhì)環(huán)境監(jiān)測(cè)站提供的災(zāi)情和險(xiǎn)情信息，本次持續(xù)性降水過程中，地質(zhì)災(zāi)害主要隨雨帶的移動(dòng)而發(fā)生變化，共兩個(gè)縣出現(xiàn)了地質(zhì)災(zāi)害(見表1).20日地質(zhì)災(zāi)害主要出現(xiàn)在萬源市竹峪鎮(zhèn)一帶，21日地質(zhì)災(zāi)害主要出現(xiàn)在大竹縣蒲包鄉(xiāng)和萬源市皮窩廟坡一帶，22日大竹縣蒲包鄉(xiāng)一帶、23日萬源皮窩鄉(xiāng)一帶均再次出現(xiàn)地質(zhì)災(zāi)害.地質(zhì)災(zāi)害類型包括滑坡、泥石流和崩塌.此次地質(zhì)災(zāi)害共威脅77戶居民，291人；轉(zhuǎn)移受災(zāi)群眾24人，20戶86人成功避險(xiǎn)，直接經(jīng)濟(jì)損失176萬元，潛在經(jīng)濟(jì)損失1450萬元.針對(duì)上述地質(zhì)災(zāi)害，選取了6月19日20時(shí)至21日20時(shí)的地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)報(bào)產(chǎn)品進(jìn)行檢驗(yàn).

表1 2019年6月21日-22日地質(zhì)災(zāi)害

3.3 致災(zāi)分析

地質(zhì)災(zāi)害受降水、巖組特征和人類經(jīng)濟(jì)活動(dòng)等多種因素影響.達(dá)州地勢(shì)由北向南傾斜，大巴山橫直在萬源、宣漢北部，明月山、銅鑼山、華鎣山由北而南縱臥其間，屬“川東褶皺剝蝕—侵蝕低山丘陵嶺谷區(qū)”.區(qū)內(nèi)地質(zhì)災(zāi)害主要發(fā)育于低山、深丘部位，中丘、淺丘部位次之.根據(jù)劉德兵[1]等對(duì)達(dá)州地質(zhì)災(zāi)害區(qū)劃的研究，達(dá)州市地質(zhì)災(zāi)害易發(fā)區(qū)的主要分布區(qū)域包括：東北部大巴山中山峽谷易發(fā)區(qū)、西南部低山深丘易發(fā)區(qū)、東南部華鎣山低山易發(fā)區(qū)、人類工程活動(dòng)影響區(qū)和重點(diǎn)防治隱患點(diǎn)區(qū)域.進(jìn)一步分析地質(zhì)災(zāi)害發(fā)生地，此次過程發(fā)生地質(zhì)災(zāi)害的地方均位于高易發(fā)區(qū).

3.4 地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)警產(chǎn)品

根據(jù)前期降水與未來24小時(shí)格點(diǎn)預(yù)報(bào)，達(dá)州市地質(zhì)災(zāi)害預(yù)報(bào)預(yù)警平臺(tái)從20日到22日連續(xù)三天預(yù)報(bào)了地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)Ⅲ級(jí)以上等級(jí).20日萬源和宣漢北部有Ⅲ級(jí)以上風(fēng)險(xiǎn)，其中，萬源西北部有Ⅱ級(jí)風(fēng)險(xiǎn)；21日萬源北部有Ⅱ級(jí)到Ⅲ級(jí)風(fēng)險(xiǎn)；22日全市大部地方有Ⅱ級(jí)到Ⅲ級(jí)風(fēng)險(xiǎn).

4 預(yù)警檢驗(yàn)結(jié)果與分析

4.1 預(yù)警檢驗(yàn)

目前地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)警業(yè)務(wù)只發(fā)布Ⅲ級(jí)及以上的地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)等級(jí)，下面對(duì)Ⅲ級(jí)及以上氣象風(fēng)險(xiǎn)等級(jí)預(yù)報(bào)進(jìn)行不分級(jí)別的檢驗(yàn).

(1)按照《標(biāo)準(zhǔn)》檢驗(yàn)，20日準(zhǔn)確率100%，21日50%，22日100%，三天地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)警準(zhǔn)確率達(dá)到83.3%.

(2)按照《規(guī)范》進(jìn)行檢驗(yàn)，結(jié)果如表2所示：在“6.20-6.22”持續(xù)性強(qiáng)降水引發(fā)的地質(zhì)災(zāi)害中，全市6個(gè)縣在6月20日至22日期間，地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)警產(chǎn)品預(yù)報(bào)正確3個(gè)、漏報(bào)1個(gè)、空?qǐng)?bào)6個(gè)、無預(yù)警服務(wù)產(chǎn)品預(yù)報(bào)正確8個(gè)，預(yù)報(bào)命中率61.1%，漏報(bào)率5.6%，空?qǐng)?bào)率33.3%.根據(jù)這個(gè)定量檢驗(yàn)結(jié)果，在業(yè)務(wù)使用時(shí)需要注意消除空?qǐng)?bào)率，特別需要考慮地質(zhì)環(huán)境和前期累積雨量以及后期短時(shí)強(qiáng)降水等因素的影響.

4.2 萬源地質(zhì)災(zāi)害個(gè)例與降水關(guān)系分析

6月15日開始，萬源出現(xiàn)持續(xù)降雨，19日20時(shí)至20日20時(shí)，萬源北部普遍出現(xiàn)大到暴雨，萬源國家氣象觀測(cè)站59.6mm， 21日普遍小到中雨.受此影響，6月20日到21日萬源市竹峪、虹橋及廟坡一帶多處出現(xiàn)滑坡、泥石流等地質(zhì)災(zāi)害.

地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)警需綜合考慮降水、地質(zhì)條件等因素.降水不但要考慮持續(xù)降水影響，還與雨強(qiáng)有關(guān).出現(xiàn)地質(zhì)災(zāi)害的竹峪、虹橋一帶，前期不但持續(xù)降水，19日夜間到20日白天還出現(xiàn)短時(shí)強(qiáng)降水.以竹峪鎮(zhèn)為例，從圖2可看出，20日07時(shí)開始雨強(qiáng)加大，08時(shí)小時(shí)雨強(qiáng)21.1 mm·h-1，09時(shí)小時(shí)雨強(qiáng)達(dá)86 mm·h-1，同一時(shí)段虹橋小時(shí)雨強(qiáng)也達(dá)42.6 mm·h-1.

圖2 2019年6月19日20時(shí)-20日20時(shí)萬源市竹峪鎮(zhèn)逐小時(shí)雨量

廟坡鄉(xiāng)在20日、21日均出現(xiàn)了滑坡，僅從19日20時(shí)至20日20時(shí)雨量看，廟坡區(qū)域站12.7mm，但出現(xiàn)地質(zhì)災(zāi)害的點(diǎn)靠近花萼山和皮窩一帶，皮窩、花萼山雨量分別為47mm、82.5mm，且20日10時(shí)至12時(shí)分別出現(xiàn)了33.4 mm·h-1、59.8 mm·h-1的短時(shí)強(qiáng)降水(見圖3).

圖3 2019年6月19日20時(shí)-21日20時(shí)萬源市皮窩鄉(xiāng)逐小時(shí)雨量

從地質(zhì)條件方面分析，萬源北部山區(qū)屬于地質(zhì)災(zāi)害重點(diǎn)防治區(qū)，20日主要發(fā)生的滑坡區(qū)屬于竹峪-虹橋滑坡高易發(fā)亞區(qū)，位于萬源市西北部，地貌以構(gòu)造侵蝕中山嶺脊-峽谷地形和構(gòu)造侵蝕中山單斜山-不對(duì)稱“V”型谷為主，海拔520～1600m，地形起伏較大，切割深，出露三疊和侏羅系地層，巖性以砂巖、泥巖、灰?guī)r、礫巖、頁巖為主，中間夾雜煤層，地質(zhì)環(huán)境脆弱，現(xiàn)狀地質(zhì)災(zāi)害發(fā)育.

21日主要發(fā)生的廟坡一帶滑坡區(qū)屬于鐘?！舷赂咭装l(fā)亞區(qū)，位于萬源市東北部，屬大竹河流域，地貌以侵蝕溶蝕高中山巖溶褶皺山和巖溶斷塊山為主，海拔520～2200m，最大高差近1700m，區(qū)域北西-南東向緊密褶皺和斷層均較為發(fā)育，出露地層有震旦、寒武、奧陶和志留系地層，巖性以灰?guī)r、白云巖、砂巖、硅質(zhì)巖、頁巖為主，局部有采煤等活動(dòng)，現(xiàn)狀地質(zhì)災(zāi)害發(fā)育.

因此，對(duì)于地質(zhì)災(zāi)害預(yù)報(bào)高易發(fā)區(qū)，在地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)警時(shí)，空?qǐng)?bào)率較小，但在分析風(fēng)險(xiǎn)等級(jí)時(shí)，不僅要考慮累積降雨量，小時(shí)雨強(qiáng)也應(yīng)該特別關(guān)注.

4.3 大竹地質(zhì)災(zāi)害個(gè)例與降水關(guān)系分析

6月21日、22日大竹縣蒲包鄉(xiāng)一帶出現(xiàn)了山體滑坡，塌方量1.065*105m3的,21日平臺(tái)漏報(bào)了該起地質(zhì)災(zāi)害.

6月20日20時(shí)至21日20時(shí)，大竹東部和北部出現(xiàn)大到暴雨，有5個(gè)站點(diǎn)雨量超過100毫米，國家觀測(cè)站降雨量121.3毫米.蒲包鄉(xiāng)42.2毫米，蒲包鄉(xiāng)附近的永勝鄉(xiāng)降雨量達(dá)218.9毫米.19日ECMWF、西南區(qū)域模式等客觀模式均漏報(bào)了這次過程；格點(diǎn)化定量降水產(chǎn)品也只考慮了小到中雨，主觀降水預(yù)報(bào)把大到暴雨時(shí)間推遲到了21日.因此，平臺(tái)預(yù)報(bào)地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)Ⅳ級(jí)，認(rèn)為發(fā)生地質(zhì)災(zāi)害的氣象風(fēng)險(xiǎn)等級(jí)較低.

在發(fā)生地質(zhì)災(zāi)害的蒲包鄉(xiāng)和永勝鄉(xiāng)附近，實(shí)況降水比預(yù)報(bào)偏大明顯.由圖4和圖5看出，地質(zhì)災(zāi)害發(fā)生地有連續(xù)降水，最大小時(shí)雨強(qiáng)分別為14.6mm·h-1、71.4 mm·h-1.雨量大、雨勢(shì)強(qiáng)、降雨時(shí)間長等均是引發(fā)蒲包鄉(xiāng)山體滑坡地質(zhì)災(zāi)害的有利氣象因素.由此可見，雨量預(yù)報(bào)偏小會(huì)造成地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)報(bào)等級(jí)偏低，導(dǎo)致預(yù)警漏報(bào).

圖5 2019年6月20日20時(shí)-22日20時(shí)大竹縣永勝鄉(xiāng)逐小時(shí)雨量

5 結(jié) 論

(1)持續(xù)性降水和短時(shí)強(qiáng)降水是造成此次地質(zhì)災(zāi)害主要原因.準(zhǔn)確的降水量級(jí)預(yù)報(bào)是地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)警的基礎(chǔ)，降水量級(jí)預(yù)報(bào)偏小、雨量預(yù)報(bào)偏小會(huì)造成地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)報(bào)等級(jí)偏低，導(dǎo)致地質(zhì)災(zāi)害漏報(bào)；反之則可能引起空?qǐng)?bào).

(2)地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)警平臺(tái)成功預(yù)報(bào)了“6.20-6.22”期間的地質(zhì)災(zāi)害，按照地質(zhì)行業(yè)的《標(biāo)準(zhǔn)》預(yù)報(bào)準(zhǔn)確率達(dá)到83.3%；按照中國氣象局《規(guī)范》命中率達(dá)到61.1%，漏報(bào)率5.6%，空?qǐng)?bào)率33.3%.各項(xiàng)指標(biāo)均表明，預(yù)警平臺(tái)有很好的指導(dǎo)作用.

(3)按照中國氣象局暴雨誘發(fā)地質(zhì)災(zāi)害氣象風(fēng)險(xiǎn)預(yù)警業(yè)務(wù)規(guī)范，達(dá)州市地質(zhì)災(zāi)害風(fēng)險(xiǎn)預(yù)警客觀預(yù)報(bào)方法空?qǐng)?bào)率較漏報(bào)率高，實(shí)際業(yè)務(wù)預(yù)警發(fā)布還需綜合預(yù)警發(fā)布前降水特點(diǎn)和降水量級(jí)、主觀降水量級(jí)訂正、以及達(dá)州地質(zhì)災(zāi)害發(fā)生特點(diǎn)等進(jìn)行訂正應(yīng)用.