■李歡仲佳鑫劉曉莉穆增光
(1寧夏回族自治區(qū)遙感測繪勘查院(寧夏回族自治區(qū)遙感中心)寧夏銀川750021;2寧夏國土資源調(diào)查監(jiān)測院寧夏銀川750021;3寧夏地質(zhì)博物館寧夏銀川750021;4寧夏國土資源地理信息中心寧夏銀川750021)
寧夏降水型地質(zhì)災害潛勢預報模型的創(chuàng)建研究
■李歡1仲佳鑫2劉曉莉3穆增光4
(1寧夏回族自治區(qū)遙感測繪勘查院(寧夏回族自治區(qū)遙感中心)寧夏銀川750021;2寧夏國土資源調(diào)查監(jiān)測院寧夏銀川750021;3寧夏地質(zhì)博物館寧夏銀川750021;4寧夏國土資源地理信息中心寧夏銀川750021)
本文基于對1982~2004年寧夏地質(zhì)災害與降水資料進行分析,并在引發(fā)寧夏干旱區(qū)主要地質(zhì)災害相應氣象條件的框架下,采用統(tǒng)計學方法實施建模操作,最終就寧夏降水型地質(zhì)災害潛勢預報模型予以建立。
寧夏降水型地質(zhì)災害潛勢預報模型
寧夏作為地質(zhì)災害高發(fā)區(qū),其不僅在災害類型上較多,且在分布方面也較廣。近些年來由于地質(zhì)災害所造成的經(jīng)濟損失高達上百億,且還造成了重大的人員傷亡。至此,寧夏回族自治區(qū)遙感中心借助自治區(qū)科技廳自然科學基金開展本次地質(zhì)災害預報研究,且通過相關部門提供的資料,并就氣象部門所提供的降水資料予以結合,最后就寧夏地質(zhì)災害在具體形成方面的相應規(guī)律及特點進行分析,創(chuàng)建了寧夏降水型地質(zhì)災害潛勢預報模型。
由于泥石流及滑坡等地質(zhì)災害較多發(fā)生于局地,且實質(zhì)為小概率事件,故此具有較大的預報難度,然而經(jīng)較多資料統(tǒng)計分析可知,寧夏降水方面和泥石流及滑坡等災害的發(fā)生具有緊密關聯(lián),并且在時間方面則具有通常的一致性。此種關系首先在年際變化方面予以表現(xiàn),通過對22年以來相應資料進行分析可知,寧夏地質(zhì)災害發(fā)生方面密切相關與大氣降水周期,也可說未相關于地殼活動的影響。雖然在1996年~2003年在具體的年均降水量方面并不算大,然而發(fā)生強降雨及暴雨方面的相應高發(fā)頻次,則與地質(zhì)災害易發(fā)區(qū)內(nèi)正好重疊,所以在此區(qū)域內(nèi)有較多地質(zhì)災害發(fā)生,比如在1984年一級1985年期間,在地質(zhì)災害易發(fā)區(qū)內(nèi),由于未出現(xiàn)相應強降雨及暴雨,因此,此地在所造成的地質(zhì)災害方面也相對較少。另外,在月際變化方面也有表現(xiàn),寧夏地區(qū)所發(fā)生的全部地質(zhì)災害均在相應4~10月,此時間段也是個月在具體的平均降水量方面較多的季節(jié),特別是在本地區(qū)的7~8月,其在具體的地質(zhì)災害發(fā)生頻次方面占據(jù)總體比重的80%,此狀況同時也是所有季節(jié)當中暴雨發(fā)生次數(shù)最多的,而二者之間也呈現(xiàn)出較好的對應關系。
通過對寧夏多年來地質(zhì)災害和氣象資料實施綜合統(tǒng)計分析,初步將寧夏地質(zhì)災害和降水之間相應關系予以建立。依據(jù)寧夏地質(zhì)及氣候等所存在的差異,此外,還根據(jù)地質(zhì)災害發(fā)生的具體頻次狀況,將寧夏在地質(zhì)災害易發(fā)區(qū)域方面實施多個劃分操作,依據(jù)各個區(qū)域相應雨量站各次地質(zhì)災害發(fā)生時間,就對應的雨量資料進行統(tǒng)計;并對各個區(qū)域內(nèi),與各次地質(zhì)災害對應時間段,在具體的最大平均雨量方面進行計算,然后將N次地質(zhì)災害方面各時段最大雨量平均值方面的最小值進行統(tǒng)計,其也就是各時間段地質(zhì)災害臨界雨量初值。
2.1 預報因子的選取
依據(jù)上述當中的降雨分型,選取其中的5個和降水有關的預報因子,依據(jù)已經(jīng)存在的諸如蔣家溝模型等研究成果,首先將有效雨量設定為Ra,將其定義為:
在公式當中,Ra表示為有效雨量;RL則表示為自災害發(fā)生之日起相應前1天算起的相應雨量,其中災害發(fā)生的當天為l=0,而在災害發(fā)生前1天則為l=1,而前2天則表示為L=2;K則表示為遞減系數(shù)。將公式(1)經(jīng)轉(zhuǎn)換可得出:
從公式(2)中可知,其實際計算的雨量為10天之后的,其在有效雨量的貢獻方面則不大,如若將14天為例,則0.814則是一個較小的數(shù)。所以,在結果方面和采用20天雨量計算方面所得出的有效雨量方面所存在的差別并不大。
2.2 地質(zhì)災害潛勢預報建模
根據(jù)以往方法,在實際建立模型過程中,可將地質(zhì)災害發(fā)生概率為10%定義為可能性很小,25%定義為可能性較小,50%定義為可能性大,而95%則定義為可能性很大,而這些均為臨界概率。而將與這些概率相對應的因子取值,則可當做寧夏各地區(qū)地質(zhì)災害發(fā)生概率方面相應臨界值。
2.2.1 將因子的災害潛勢等級a(1)予以確定
現(xiàn)以寧夏中北部地區(qū)為例,針對第一個因子,就當天降雨量f(i)=f(1)為例:如若f(1)<1mm,則寧夏中北部地質(zhì)災害發(fā)生數(shù)在整個寧夏災害發(fā)生數(shù)的比重為10%之下,則取a(1)=1;如若1mm≤f(1)<5mm,則寧夏中北部地質(zhì)災害發(fā)生數(shù)在整個寧夏災害發(fā)生數(shù)的比重為10%之上,則可取a(1)=2;如若5mm≤f(1)<26mm,則寧夏中北部地質(zhì)災害發(fā)生數(shù)在整個寧夏災害發(fā)生數(shù)的比重為25%之上,則可取a(1)=3;如若26mm≤f(1)<65mm,則寧夏中北部地質(zhì)災害發(fā)生數(shù)在整個寧夏災害發(fā)生數(shù)的比重為50%之上,則可取a(1)=4;如若65mm≤f(1)<154mm,則寧夏中北部地質(zhì)災害發(fā)生數(shù)在整個寧夏災害發(fā)生數(shù)的比重為75%之上,則可取a(1)=5;……,如若f(1)≥250.0mm,則寧夏中北部地質(zhì)災害發(fā)生數(shù)95%之上,也就是地質(zhì)災害一定發(fā)生,可將a取值為5。
2.2.2 將總的地質(zhì)災害潛勢等級值a給與確定
將統(tǒng)計公式予以運用,即
在公式當中,a(i)=1(當f(i)<時);a(i)=2(當f1(i)≤f(i)<f2(i)時);a(i)=3(當f2(i)≤f(i)<f3(i)時);a(i)=4(當f3(i)≤f(i)<f4(i)時);a(i)=5(當f3(i)≥f4(i)時)(4)。如若因子值相比于地質(zhì)災害發(fā)生的臨界值相比,出現(xiàn)超過后者狀況后,由于本因子在地質(zhì)災害發(fā)生的條件方面已經(jīng)達到,所以地質(zhì)災害預報等級直接劃定為最后,即a=5,也就是在f(i)>f4(i) f4(i)時,將其中的a取為5。
2.2.3 對每個分區(qū)的災害潛勢等級預報予以確定
依據(jù)公式(3)相應計算結果,可將預報結果進行等級劃分,分別劃分為5級,即當a=1時,其地質(zhì)災害發(fā)生可能性很??;當a=2時,其地質(zhì)災害發(fā)生可能性較??;當a=3時,其地質(zhì)災害發(fā)生可能性較大;當a=4時,其地質(zhì)災害發(fā)生可能性大;當a=5時,其地質(zhì)災害發(fā)生可能性很大。
[1]陳豫英,趙光平,王紅英.寧夏汛期地質(zhì)災害氣象預報預警模型研究 [J].災害學,2005,20(3):25-28.
P694[文獻碼]B
1000-405X(2016)-6-414-1