安陽市地質(zhì)災(zāi)害分級(jí)預(yù)報(bào)系統(tǒng)

劉勤娣

摘要 通過分析安陽市地質(zhì)災(zāi)害時(shí)空分布特征，利用安陽市自動(dòng)雨量站資料得出降水與安陽市地質(zhì)災(zāi)害的關(guān)系判別式，運(yùn)用上級(jí)氣象臺(tái)的指導(dǎo)預(yù)報(bào)產(chǎn)品和前期降水，采取計(jì)算機(jī)程序自動(dòng)讀取數(shù)據(jù)制作未來安陽市地質(zhì)災(zāi)害分級(jí)預(yù)報(bào)。

關(guān)鍵詞 地質(zhì)災(zāi)害；時(shí)空分布；分級(jí)預(yù)報(bào)；河南安陽

中圖分類號(hào) P694 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 1007-5739（2016）05-0253-02

Abstract Through the analysis of space-time distribution characteristics of geological disasters in Anyang City，using the data of Anyang City automatic precipitation station，the relationship between precipitation and Anyang City geological disasters discriminant was obtained，forecast products and prior precipitation were forcasted guiding by superior meteorological station，geological hazard classification forecast in Anyang City were made taking the computer program automatically read data.

Key words geological disasters；space-time distribution；classification prediction；Anyang Henan

安陽市位于河南省的最北部地區(qū)，東部為平原地帶，西部為太行山山系或丘陵地區(qū)，東西落差較大，地形和地質(zhì)構(gòu)造復(fù)雜，其中的林州市和安陽縣西部地區(qū)的山丘地區(qū)面積3 424 km2（占全市總面積的46%），是地質(zhì)災(zāi)害的多發(fā)區(qū)，每年5—9月因強(qiáng)降水引發(fā)的山體滑坡、泥石流、塌陷等都有發(fā)生，直接威脅著人民群眾的生命財(cái)產(chǎn)安全，同時(shí)也給當(dāng)?shù)亟?jīng)濟(jì)發(fā)展造成很大損失[1-3]。

能夠造成地質(zhì)災(zāi)害的因素很多，包括人為的亂砍濫伐，不規(guī)范的采礦，巖石的分化演變等，但主要誘發(fā)因素是降水，特別是強(qiáng)降水造成的滲透徑流，使土石抗剪強(qiáng)度降低，極容易造成地質(zhì)災(zāi)害[4-6]。

本文主要是通過分析安陽市近20年的地質(zhì)災(zāi)害發(fā)生情況，找出地質(zhì)災(zāi)害發(fā)生的多發(fā)區(qū)和易發(fā)區(qū)，對(duì)誘發(fā)地質(zhì)災(zāi)害的強(qiáng)降水進(jìn)行分型，在此基礎(chǔ)上，利用縣局或鄉(xiāng)鎮(zhèn)自動(dòng)雨量站過去72 h的降水資料及未來24 h的降水預(yù)報(bào)，對(duì)可能發(fā)生的地質(zhì)災(zāi)害進(jìn)行定量化預(yù)報(bào)，并通過《安陽市新一代預(yù)報(bào)業(yè)務(wù)流程系統(tǒng)》實(shí)現(xiàn)預(yù)報(bào)結(jié)果的自動(dòng)化輸出。

1 地質(zhì)災(zāi)害的分布特征

1.1 空間分布

按照地質(zhì)勘查部門的調(diào)查，將其分為：一是特大型。災(zāi)害很嚴(yán)重，危害程度很大，受威脅人數(shù)大于1 000人，預(yù)測(cè)經(jīng)濟(jì)損失1 000萬元以上。二是大型。災(zāi)害嚴(yán)重，危害程度大，受威脅人數(shù)100～1 000人，預(yù)測(cè)經(jīng)濟(jì)損失500萬～1 000萬元。三是中型。災(zāi)害較嚴(yán)重，危害程度較大，受威脅人數(shù)10～100人，預(yù)測(cè)經(jīng)濟(jì)損失100萬～500萬元。四是小型。災(zāi)害較輕，危害程度不嚴(yán)重，受威脅人數(shù)小于10人，預(yù)測(cè)經(jīng)濟(jì)損失小于100萬元。

統(tǒng)計(jì)表明，近20年共發(fā)生特大型地質(zhì)災(zāi)害2次，占1.6%；大型25次，占20.3%；中型55次，占44.7%；小型39次，占31.7%（圖1）。

大型地質(zhì)災(zāi)害分布情況：在發(fā)生的25次中，以林州的任村、橫水、河順3個(gè)鄉(xiāng)鎮(zhèn)最為顯著，分別為7、5、4次（表1）。

中型地質(zhì)災(zāi)害分布情況：中型地質(zhì)災(zāi)害發(fā)生的次數(shù)最多，共發(fā)生55次中，占總次數(shù)的45%，以林州的任村、橫水、河順、紅旗渠、石板巖、東姚6個(gè)鄉(xiāng)鎮(zhèn)最為顯著，分別為10、9、9、7、5、5次（表2），6個(gè)鄉(xiāng)鎮(zhèn)合計(jì)占83%。

小型地質(zhì)災(zāi)害分布情況：在發(fā)生的39次中，除林州的原康、河順2個(gè)鄉(xiāng)鎮(zhèn)（分別為9、6次）外，其他地區(qū)分布非常分散（表3）。

1.2 時(shí)間分布

安陽市1991—2010年共發(fā)生山體滑坡、泥石流、塌陷等地質(zhì)災(zāi)害123次，平均每年發(fā)生6.2次，統(tǒng)計(jì)結(jié)果表明，地質(zhì)災(zāi)害的多發(fā)時(shí)段為每年的5—9月，集中時(shí)段為7—8月，這和強(qiáng)降水的分布呈一致性。

2 地質(zhì)災(zāi)害與降水

正如前文所述，地質(zhì)災(zāi)害的發(fā)生有2個(gè)最決定性因素，一個(gè)是當(dāng)?shù)氐牡匦魏偷刭|(zhì)構(gòu)造，一個(gè)是降水。對(duì)于一個(gè)地區(qū)而言，地形、地質(zhì)構(gòu)造和地表情況無法改變，因此地質(zhì)災(zāi)害的預(yù)報(bào)，就是對(duì)觸發(fā)條件的預(yù)報(bào)，也就是說對(duì)一個(gè)地區(qū)（鄉(xiāng)鎮(zhèn)）而言，前期累計(jì)降水的正確評(píng)價(jià)和對(duì)未來強(qiáng)降水的正確預(yù)報(bào)是地質(zhì)災(zāi)害預(yù)報(bào)的關(guān)鍵。從1991—2010年林州、安陽縣日降水大于等于30 mm的詳細(xì)情況來看，對(duì)照地質(zhì)災(zāi)害的發(fā)生情況，可以發(fā)現(xiàn)，大于30 mm的降水日數(shù)遠(yuǎn)小于地質(zhì)災(zāi)害的發(fā)生次數(shù)，因此說引發(fā)地質(zhì)災(zāi)害的主要是強(qiáng)降水，一次強(qiáng)降水，能誘發(fā)多處（次）地質(zhì)災(zāi)害，而同對(duì)于一次降水，不同地區(qū)發(fā)生地質(zhì)災(zāi)害的程度也不相同，也就是說不同地區(qū)發(fā)生地質(zhì)災(zāi)害的臨界值不同。例如：1996年8月3—4日，安陽市中西部普降大暴雨，林州市8月4日降水量192.5 mm，林州的紅旗渠等地先后12處發(fā)生了山體滑坡，紅旗渠部分渠段被山洪沖毀或被泥石流淤滿，林州市后峪村被山洪沖的面目全非。

把引發(fā)地質(zhì)災(zāi)害的降水分為當(dāng)日強(qiáng)降水型（前72 h降水不足10 mm）和累計(jì)降水型，表4給出了兩型降水與地質(zhì)災(zāi)害的統(tǒng)計(jì)關(guān)系。

3 地質(zhì)災(zāi)害等級(jí)劃分及預(yù)報(bào)

3.1 地質(zhì)災(zāi)害的等級(jí)劃分

按照國(guó)家氣象中心關(guān)于地質(zhì)災(zāi)害預(yù)報(bào)預(yù)報(bào)等級(jí)的劃分標(biāo)準(zhǔn)，將安陽地質(zhì)災(zāi)害分為5個(gè)等級(jí)：1級(jí)為可能性很?。?級(jí)為可能性較?。?級(jí)為可能性較大；4級(jí)為可能性大；5級(jí)為可能性很大。同時(shí)考慮到前文中對(duì)地質(zhì)災(zāi)害的分型，特大型在20年中僅發(fā)生了2次，當(dāng)作小概率事件，在預(yù)報(bào)中不予考慮；大型對(duì)應(yīng)5級(jí)；中型對(duì)應(yīng)4級(jí)；小型對(duì)應(yīng)3級(jí)；而3級(jí)以上為對(duì)外發(fā)布級(jí)別，1級(jí)和2級(jí)只在預(yù)報(bào)時(shí)作參考。也就是說，在預(yù)報(bào)5級(jí)時(shí)，要考慮大型地質(zhì)災(zāi)害出現(xiàn)的可能性較大，預(yù)報(bào)4級(jí)時(shí)中型地質(zhì)災(zāi)害出現(xiàn)的可能性較大，預(yù)報(bào)3級(jí)時(shí)小型地質(zhì)災(zāi)害出現(xiàn)的可能性較大。

3.2 地質(zhì)災(zāi)害的預(yù)報(bào)判別和流程

引發(fā)地質(zhì)災(zāi)害需要充足的水源動(dòng)力將泥土運(yùn)走，同時(shí)考慮到前期降水量達(dá)到一定值，水分的滲透使地表變得疏松，一次再發(fā)性中等強(qiáng)度的降水也會(huì)引發(fā)地質(zhì)災(zāi)害。因此，將地質(zhì)災(zāi)害分為當(dāng)日強(qiáng)降水型和累計(jì)降水型。建立地質(zhì)災(zāi)害預(yù)報(bào)判別式：

z=∑R+（R12+R12-24）

其中，z為地質(zhì)災(zāi)害預(yù)報(bào)等級(jí)；∑R為前72 h累積雨量（從前期站點(diǎn)報(bào)文中提取）；R12為未來12 h降水預(yù)報(bào)編碼（從當(dāng)日本站對(duì)省臺(tái)所發(fā)指導(dǎo)報(bào)文中提取，下同）；R12-24為未來12～24 h降水預(yù)報(bào)編碼（提取方法同上）。判別式中各變量由程序自動(dòng)識(shí)別，運(yùn)算結(jié)果直接輸出為地質(zhì)災(zāi)害等級(jí)。并規(guī)定，一是未來無降水不報(bào)；二是未來有降水時(shí)只考慮報(bào)文中的04、08、09、10、11、12（說明：03和07為穩(wěn)定的陣雨和小雨，不考慮，04為雷陣雨，時(shí)空分布不均勻，有可能產(chǎn)生局地的強(qiáng)降水，雖然范圍可能不大，但考慮到服務(wù)效果應(yīng)盡量不漏報(bào)，所以適當(dāng)考慮）；三是前72 h累計(jì)雨量不足10 mm（當(dāng)日強(qiáng)降水型）時(shí)前期雨量不考慮，只考慮未來的預(yù)報(bào)；四是前72 h累計(jì)雨量大于等于10 mm時(shí)（累計(jì)降水型），2項(xiàng)綜合相加判別z=∑R+（R12+R12-24），其中，z為地質(zhì)災(zāi)害等級(jí)；∑R為前72 h累計(jì)雨量（從前期站點(diǎn)報(bào)文中提取）；R12為未來12 h降水預(yù)報(bào)量，R12-24為未來12～24 h降水量。

4 程序設(shè)計(jì)

實(shí)現(xiàn)在《安陽新一代天氣預(yù)報(bào)業(yè)務(wù)流程系統(tǒng)》下資料自動(dòng)提取，自動(dòng)運(yùn)行，預(yù)報(bào)結(jié)果自動(dòng)化化輸出。一是從前期各站報(bào)文中自動(dòng)提取前54 h雨量（或手工輸入鄉(xiāng)鎮(zhèn)自動(dòng)雨量站點(diǎn)前72 h雨量資料）并求合計(jì)。二是校對(duì)前期雨量，缺報(bào)或網(wǎng)絡(luò)不通時(shí)人為的手工補(bǔ)齊資料。三是從當(dāng)日對(duì)省臺(tái)的指導(dǎo)報(bào)報(bào)文中提取未來12、12～24 h降水預(yù)報(bào)報(bào)文，提取遇故障時(shí)手工輸入。四是綜合判別運(yùn)算。五是輸出地質(zhì)災(zāi)害預(yù)報(bào)結(jié)果。六是把預(yù)報(bào)結(jié)果在相應(yīng)的地理位置上閃動(dòng)警示（圖2）。

說明：如果輸入的過去雨量和預(yù)報(bào)是縣局的資料，程序輸出的結(jié)果為大型地質(zhì)災(zāi)害多發(fā)區(qū)的預(yù)報(bào)級(jí)別，對(duì)外發(fā)布時(shí)對(duì)中型和小型地質(zhì)災(zāi)害多發(fā)區(qū)的預(yù)報(bào)應(yīng)逐級(jí)遞減一個(gè)級(jí)別。如果輸入的為鄉(xiāng)鎮(zhèn)自動(dòng)雨量站的資料，程序輸出的結(jié)果即為那個(gè)鄉(xiāng)鎮(zhèn)的預(yù)報(bào)結(jié)果。

5 預(yù)報(bào)檢驗(yàn)

2011年5—9月，預(yù)報(bào)無3級(jí)以上地質(zhì)災(zāi)害，實(shí)況安陽區(qū)域內(nèi)未出現(xiàn)3級(jí)以上地質(zhì)災(zāi)害。

2012年7月29日發(fā)布林州和安陽縣西部4級(jí)地質(zhì)災(zāi)害預(yù)警。實(shí)況是7月30日至8月1日，林州地區(qū)出現(xiàn)暴雨，局部大暴雨天氣過程。大部分鄉(xiāng)鎮(zhèn)雨量都超過了100 mm，達(dá)到大暴雨量級(jí)，最大降雨中心出現(xiàn)在臨淇、茶店、五龍3個(gè)鄉(xiāng)鎮(zhèn)，雨量分別為236.6、231.7、212.0 mm。實(shí)況是林州臨期、五龍和安陽縣銅業(yè)發(fā)生不同程度山體滑坡災(zāi)害。

2013年7月7日預(yù)報(bào)8—14日林州和安陽縣西部3級(jí)地質(zhì)災(zāi)害預(yù)警，實(shí)況是從8日20：00至9日6：00，林州的五龍鎮(zhèn)和馬安腦降水分別達(dá)到74.0、62.3 mm，林州市馬安腦發(fā)生山體滑坡，部分路段被毀壞。

2014年未預(yù)報(bào)3級(jí)以上地質(zhì)災(zāi)害，實(shí)況是安陽區(qū)域內(nèi)未發(fā)生3級(jí)以上地質(zhì)災(zāi)害。

6 結(jié)語

系統(tǒng)應(yīng)用了鄉(xiāng)鎮(zhèn)自動(dòng)雨量站的降水資料，使預(yù)報(bào)更精細(xì)化；運(yùn)行后能夠直接輸出地質(zhì)災(zāi)害的分級(jí)預(yù)報(bào)結(jié)果，使預(yù)報(bào)定量化，減少了每個(gè)預(yù)報(bào)員之間的人為預(yù)報(bào)偏差；系統(tǒng)運(yùn)行方便快捷，為預(yù)報(bào)員節(jié)約了時(shí)間。

7 參考文獻(xiàn)

[1] 林孝松，郭躍.滑坡與降雨的耦合關(guān)系研究[J].災(zāi)害學(xué)，2001，16（2）：87-92.

[2] 陳百煉.降水誘發(fā)地質(zhì)災(zāi)害的氣象預(yù)警方法研究[J].貴州氣象，2002，26（4）：4-7.

[3] 劉傳正.區(qū)域地質(zhì)災(zāi)害預(yù)警的理論與方法研究[J].水文地質(zhì)工程地質(zhì)，2004（3）：1-6.

[4] 靳冰凌.鶴壁市地質(zhì)災(zāi)害預(yù)報(bào)預(yù)警系統(tǒng)[J].河南氣象，2005（2）：23-24.

[5] 賀為民.地質(zhì)災(zāi)害危險(xiǎn)性評(píng)估分級(jí)因素的探討[J].災(zāi)害學(xué)，2013，28（3）：112-116.

[6] 張國(guó)平，宋建洋，邵小路，等.近期我國(guó)大型地質(zhì)災(zāi)害預(yù)報(bào)與避險(xiǎn)案例分析[J].自然災(zāi)害學(xué)報(bào)，2015（3）：20-27.