藏東南直白溝冰崩型泥石流孕災(zāi)條件與動(dòng)力特征

龍 飛龔 誠黃 海焦成才

〔1.西藏自治區(qū)地質(zhì)環(huán)境監(jiān)測總站,西藏 拉薩850000;2.中國地質(zhì)科學(xué)院探礦工藝研究所,四川 成都611734;3.中國地質(zhì)大學(xué)(武漢)工程學(xué)院,湖北 武漢430074;4.北京歐遠(yuǎn)致科技有限公司,北京100085〕

地震和全球氣候變暖是青藏高原地質(zhì)災(zāi)害的兩個(gè)關(guān)鍵驅(qū)動(dòng)力[1-2]。地震在誘發(fā)同震災(zāi)害的同時(shí)[3],還對(duì)周邊巖體具有震裂效應(yīng)[4-5],而氣溫升高則對(duì)冰川消退和高山區(qū)的凍融作用具有明顯加劇作用[6],兩者的疊加導(dǎo)致冰凍圈的地震影響區(qū)泥石流非?；钴S。如1950年墨脫地震影響下[7],直白溝和古鄉(xiāng)溝相繼暴發(fā)特大堵江型冰川泥石流災(zāi)害[8-9],2017年米林地震后,色東普在4個(gè)水文年暴發(fā)了7次泥石流[10-11]。前人對(duì)于周邊區(qū)域不同類型冰川泥石流開展了大量研究,如冰崩型泥石流(天摩溝)[12]、冰水復(fù)合型泥石流(古鄉(xiāng)溝)[13]、冰湖潰決型泥石流(光謝錯(cuò))等[14],探索了冰川泥石流災(zāi)害演化、觸發(fā)條件、啟動(dòng)機(jī)制、動(dòng)力過程及風(fēng)險(xiǎn)評(píng)價(jià),為深入認(rèn)識(shí)冰川泥石流災(zāi)害提供了重要理論基礎(chǔ)。在取得豐碩成果的同時(shí),由于冰川泥石流多發(fā)育于人跡罕至的極高山區(qū),現(xiàn)場觀測與監(jiān)測工作開展較少,使得我們對(duì)于冰川泥石流的觸發(fā)條件、動(dòng)力參數(shù)等特征還存在較多不足。

2017年11月18日,藏東南南迦巴瓦地區(qū)發(fā)生Ms 6.9級(jí)地震,震中位于色東普溝源頭(29.89°N,95.04°E)[15]。地震在雅魯藏布江大峽谷附近誘發(fā)大量同震地質(zhì)災(zāi)害[16],大部分為冰凍圈的巖土體失穩(wěn),這些災(zāi)害多數(shù)都與冰川或者高山積雪相關(guān)。同時(shí),地震作用還破壞了周邊巖體和冰川的穩(wěn)定性,震后高位冰崩、巖崩高發(fā)(以色東普溝為代表),這些高位地質(zhì)體失穩(wěn)后,沿溝道鏟刮侵蝕,形成大規(guī)模泥石流,并多次堵塞雅魯藏布江形成堰塞湖,回淹范圍達(dá)上百公里[17]。2020年9月10日下午5:00左右,直白溝暴發(fā)大規(guī)模泥石流,造成直白大橋損毀,公路路面淤埋破壞,中斷交通近半年之久,所幸未造成人員傷亡。同時(shí),大量泥石流固體物質(zhì)沖入雅魯藏布江,堵塞主河行洪斷面50%?！?·10”冰川泥石流災(zāi)害暴發(fā)時(shí)溫度與氣溫均無明顯異常,其觸發(fā)機(jī)制有別于其他震區(qū)的泥石流事件[3]。本文通過野外調(diào)查、遙感解譯、現(xiàn)場觀測和災(zāi)害監(jiān)測等研究手段,分析該次泥石流災(zāi)害事件的觸發(fā)條件、演進(jìn)過程和動(dòng)力參數(shù),為進(jìn)一步認(rèn)識(shí)冰崩型泥石流的觸發(fā)與運(yùn)動(dòng)具有重要現(xiàn)實(shí)意義和理論意義。

1 流域概況

直白溝位于西藏自治區(qū)林芝市米林縣派鎮(zhèn),為雅魯藏布江右岸一級(jí)支流。流域地處青藏高原喜馬拉雅東構(gòu)造和雅魯藏布江水汽通道交匯處[18-19],鄰近雅魯藏布江大拐彎。研究區(qū)地質(zhì)背景復(fù)雜,活動(dòng)構(gòu)造發(fā)育,地震頻發(fā),地貌類型復(fù)雜多樣。直白溝流域面積57.50 km2,主溝溝長11.50 km,地勢總體東高西低,海拔最高點(diǎn)為著名的南迦巴瓦峰,海拔高程為7 527 m,最低海拔為主河交匯口,海拔2 953 m,相對(duì)高差達(dá)4 574 m。溝道上游發(fā)育著名的則隆弄冰川,冰川的進(jìn)退變化為泥石流孕災(zāi)提供了豐富的冰磧物物源。受冰川作用,流域水系復(fù)雜,一直延伸到出山口的冰川側(cè)磧壟將匯流通道分割成3個(gè)水系,為泥石流的運(yùn)動(dòng)提供了有利的溝道地形條件。

直白溝流域復(fù)雜的地質(zhì)構(gòu)造和充沛水汽補(bǔ)給為泥石流的啟動(dòng)提供了驅(qū)動(dòng)力。在地震和氣候變暖的內(nèi)外動(dòng)力耦合作用下,研究區(qū)歷史上多次暴發(fā)泥石流災(zāi)害,對(duì)溝口直白村居民區(qū)、交通設(shè)施以及河流等危害嚴(yán)重(表1)?？傮w上,直白溝泥石流具有冰川失穩(wěn)和巖土體失穩(wěn)兩種激發(fā)模式,1950年和1968年暴發(fā)的3次泥石流災(zāi)害均由冰川躍動(dòng)誘發(fā);2020年泥石流災(zāi)害則為高位冰巖崩誘發(fā)。

表1 直白溝泥石流災(zāi)害事件統(tǒng)計(jì)

如圖1所示,2020年9月10日直白溝暴發(fā)的大規(guī)模泥石流損毀了溝口直白大橋和公路,并堵塞主河。

圖1 直白溝2020年9月10日泥石流災(zāi)害災(zāi)損情況

泥石流中的巨塊石撞擊在凈空約25 m高的直白大橋,造成主橋橋面下沉15 cm,水平錯(cuò)位移動(dòng)30 cm,泥石流從溝道右岸平臺(tái)漫流至公路,淤埋公路約150 m,所幸未造成人員傷亡,大量泥石流固體物質(zhì)沖入雅魯藏布江,堵塞主河行洪斷面50%。

2 泥石流形成條件

2.1 流域地貌特征

直白溝發(fā)育典型冰川地貌,主溝溝谷呈現(xiàn)U形地貌,溝谷寬度達(dá)800～1 200 m,在強(qiáng)烈構(gòu)造運(yùn)動(dòng)作用下,形成陡峻溝床縱比降,主溝平均縱比降203.73‰。根據(jù)泥石流形成過程,可將直白溝自上游而下劃分為形成區(qū)、流通補(bǔ)給區(qū)和堆積區(qū)(圖2)。

圖2 直白溝溝道縱剖面圖

海拔3 750 m以上主溝匯水區(qū)內(nèi)冰川發(fā)育,基巖陡立,發(fā)育零星高位崩滑體,冰川和巖體的失穩(wěn)常誘發(fā)大規(guī)模泥石流,劃分為泥石流形成區(qū);出山口以上段溝道兩岸分布冰磧物臺(tái)地,溝道下切嚴(yán)重,形成陡峻的岸坡,尤其是冰磧物岸坡,形成近90°的直立陡崖,在溝道流水側(cè)蝕作用下,極易形成補(bǔ)給泥石流的松散物源,劃分為泥石流流通補(bǔ)給區(qū);出山口至交匯口段為泥石流堆積區(qū)。因此形成區(qū)陡峻的溝道和斜坡,有利于泥石流物源的啟動(dòng),流通區(qū)順直深切的水力斷面在泥石流輸移過程中產(chǎn)生“束水攻沙”效應(yīng),有利于物質(zhì)輸移和泥石流龍頭的形成[20]。

2.2 物源特征

通過多期遙感解譯和野外調(diào)查,查出直白溝物源主要有3類:高位崩滑物源、冰磧物源和溝床物源。流域物源儲(chǔ)量豐富,總量估算為6.45×106m3,其中高位崩滑物源量約7.50×105m3,主要分布于海拔5 000 m以上的基巖裸露區(qū),其形成與凍融循環(huán)作用和地震活動(dòng)密切相關(guān);冰磧物源量約4.70×106m3,主要分布于冰川下游段主溝溝道兩側(cè)和支溝源頭;溝床物源量約1.00×106m3,主要分布于主溝溝道中,為泥石流、碎屑流及斜坡侵蝕等搬運(yùn)作用在溝道中堆積的松散物質(zhì)。

2.3 冰川發(fā)育特征

由歷史災(zāi)害事件資料分析[8],冰川是直白溝泥石流重要孕災(zāi)體。受氣候變暖影響,直白溝冰川退縮嚴(yán)重,主溝冰川長度由1950年的7.4 km縮短至現(xiàn)在的5.9 km。冰川變化速度總體呈現(xiàn)慢—快—慢的趨勢,近20 a來冰川變化較慢,末端高程變化約10 m[21]。冰川長度的縮小和厚度的減薄均有利于冰川整體穩(wěn)定性,其賦存狀態(tài)的變化對(duì)直白溝泥石流的啟動(dòng)機(jī)制具有顯著影響,相較于20世紀(jì)50年代,現(xiàn)狀由冰川運(yùn)動(dòng)直接誘發(fā)的泥石流的可能性已經(jīng)大大降低。

2.4 氣象條件

水熱組合和長期氣候變化是冰川泥石流的重要誘因,冰川泥石流觸發(fā)時(shí)的氣象條件不僅僅由降雨決定,如色東普、天摩溝等,均在少雨的情況下暴發(fā)泥石流[22]。研究區(qū)對(duì)全球氣候變暖的響應(yīng)非常敏感,1981年以來,青藏高原地區(qū)氣溫升溫0.60℃/10 a。米林縣汛期平均氣溫與1—10月平均氣溫均為線性增長趨勢,汛期增速0.35℃/10 a,1—10月平均氣溫增速為0.39℃/10 a。在快速增溫的背景下,冰川不穩(wěn)定性增加,消融加劇,冰崩誘發(fā)的泥石流還將持續(xù)甚至加強(qiáng),危險(xiǎn)性增加。綜上所述,直白溝溝域內(nèi)賦存豐富的松散物源,在氣候變化下的冰雪運(yùn)動(dòng)強(qiáng)烈,溝道堵塞嚴(yán)重且縱比降陡峻,具備泥石流形成的三大主要條件。因此在水熱條件的有利組合誘發(fā)下,直白溝于2020年9月10日暴發(fā)特大規(guī)模泥石流災(zāi)害。

3 泥石流災(zāi)害特征

3.1 觸發(fā)條件

根據(jù)安裝于直白溝溝口(海拔2 900 m)的降雨量和附近白努弄巴(海拔3 000 m)的氣溫監(jiān)測數(shù)據(jù),“9·10”泥石流暴發(fā)當(dāng)日未出現(xiàn)集中降雨,溝口24 h降雨量僅有1.00 mm,氣溫在一周內(nèi)也無明顯強(qiáng)烈變化(圖3)。從全年降雨量分布來看,直白溝溝口的降雨量監(jiān)測顯示9月的降雨量僅為23.50 mm,僅占全年降雨總量的6.80%,屬于少雨季節(jié)(表2)。由此可見,直白溝的“9·10”泥石流并非由于短周期內(nèi)的氣候突變誘發(fā),而是長期以來的氣候變化疊加2017年的米林地震的影響,高寒山區(qū)強(qiáng)烈的凍融作用驅(qū)使高位巖體失穩(wěn)而形成的泥石流。對(duì)比1968年暴發(fā)的冰川躍動(dòng)泥石流事件,同樣為無降雨的情況下誘發(fā)災(zāi)害,由此推測直白溝冰川泥石流與周邊的冰水復(fù)合型泥石流災(zāi)害在啟動(dòng)機(jī)制上存在差異[23],其觸發(fā)條件中水熱優(yōu)勢組合并非唯一條件。

圖3 “9·10”泥石流暴發(fā)時(shí)氣溫與降雨條件

表2 直白溝溝口2020年降雨量監(jiān)測值分月統(tǒng)計(jì)

3.2 泥石流災(zāi)害過程

通過泥石流災(zāi)害暴發(fā)前后的遙感影像分析發(fā)現(xiàn),“9·10”泥石流主要是在2#支溝活動(dòng)。受2017年米林地震影響,2#支溝源頭發(fā)育了一處同震崩塌,形成冰巖崩碎屑流堆積于山坡冰磧臺(tái)地上(圖4a)。2020年9月10日,上述崩塌所在海拔5 500 m的分水嶺再次發(fā)生高位巖崩,大量碎屑物沿斜坡傾瀉而下,并撞擊下部平臺(tái)上的冰川,致使冰川前緣發(fā)生冰崩(圖4b)。冰塊混合巖土體一起沖擊在冰磧平臺(tái)上,鏟刮2017年形成的冰巖崩堆積體和冰磧物,快速運(yùn)動(dòng)的冰巖碎屑流進(jìn)入狹窄的溝道中,在冰雪的相變作用和溝道徑流補(bǔ)給作用下,碎屑流逐步轉(zhuǎn)化為泥石流。高速運(yùn)動(dòng)的流體沿途鏟刮溝道兩側(cè)冰磧壟和溝床物質(zhì),不斷補(bǔ)給松散固體物質(zhì)并放大災(zāi)害規(guī)模。在泥石流沖出出山口后,受溝谷寬度擴(kuò)張影響,龍頭撞擊岸坡沖起部分流體沿著溝道右岸平臺(tái)漫流,大量松散物質(zhì)堆積于公路路面。泥石流主流部分從直白大橋橋下排導(dǎo)進(jìn)入雅魯藏布江,將原160 m寬的主河擠占至80 m。

圖4 “9·10”泥石流啟動(dòng)區(qū)遙感影像

從現(xiàn)場拍攝的觀測視頻上看,“9·10”冰川泥石流具有典型黏性泥石流特征,流體結(jié)構(gòu)整體性強(qiáng),巨塊石含量高,滾動(dòng)的巨塊石對(duì)流速和峰值流量均具有放大效應(yīng)(如圖5所示)。

圖5 “9·10”冰川泥石流觀測視頻截圖(受監(jiān)測設(shè)備影響圖像精度較低)

3.3 泥石流溝道沖淤特征

利用泥石流前后遙感影像,基于溝道微地貌變化,“9·10”泥石流過程中的溝道沖淤特征上可劃分為4個(gè)區(qū)段(如圖6所示)。2#支溝匯入主溝之前主要表現(xiàn)為沖刷作用為主,但在不同區(qū)段,由于災(zāi)害體物理性質(zhì)差異,在不同時(shí)空尺度上表現(xiàn)出差異化的沖淤特征。在源頭斜坡區(qū)域主要表現(xiàn)為以沖刷為主,部分由陡變緩的地形突變段存在明顯淤積,顯示出物質(zhì)交換現(xiàn)象。當(dāng)災(zāi)害體為冰巖碎屑流時(shí),其運(yùn)動(dòng)性遠(yuǎn)低于固液兩相的泥石流,但是這種高勢能驅(qū)動(dòng)下的高流速碎屑流對(duì)斜坡依然具有強(qiáng)烈的沖刷作用,災(zāi)害體前端快速改變斜坡地貌,侵蝕沖刷形成溝槽(圖4b),但碎屑流末端的能量相對(duì)較弱,在地形平緩段出現(xiàn)局部的堆積。當(dāng)碎屑流高速進(jìn)入2#支溝后,碎屑流中的冰和雪產(chǎn)生相變,不斷補(bǔ)給水體進(jìn)入到流體中,碎屑流逐步轉(zhuǎn)化為泥石流,并在沿溝道運(yùn)動(dòng)過程中,侵蝕岸坡和溝床,以滾雪球方式不斷放大災(zāi)害規(guī)模,如圖6d所示。泥石流進(jìn)入主溝后,得到主溝充沛的徑流補(bǔ)給和稀釋,同時(shí)主溝溝道遠(yuǎn)寬于支溝溝道,泥石流流速逐步減少,但依然具有一定的沖刷能力,在侵蝕岸坡和溝床底部的同時(shí)也沿彎道的寬緩平臺(tái)產(chǎn)生少量堆積(圖6c),總體表現(xiàn)為沖淤平衡狀態(tài)(圖6b)。泥石流進(jìn)入堆積區(qū)后,災(zāi)害體沿堆積扇擴(kuò)散堆積,塑造新的堆積扇形態(tài)(圖6a)。泥石流在溝口產(chǎn)生大量固體物質(zhì)堆積,堆積范圍面積達(dá)6.50×104m2,平均堆積厚度約4 m,根據(jù)堆積厚度特征,可劃分為4個(gè)堆積區(qū)(見圖1),其中Ⅰ區(qū)堆積厚度0.2～2 m,Ⅱ區(qū)堆積厚度1～3 m,Ⅲ區(qū)堆積厚度2～5 m,Ⅳ區(qū)堆積厚度3～10 m,估算總計(jì)沖出固體物質(zhì)2.5×105m3,其中進(jìn)入主河約2.00×105m3。

圖6 “9·10”2#支溝泥石流沿溝道的沖淤特征

3.4 泥石流動(dòng)力參數(shù)計(jì)算

“9·10”泥石流暴發(fā)前,西藏自治區(qū)地質(zhì)災(zāi)害主管部門在直白溝開展了地質(zhì)災(zāi)害監(jiān)測工作,布設(shè)了次聲儀、泥位計(jì)、降雨溫度站以及視頻等監(jiān)測設(shè)備。在該次泥石流暴發(fā)時(shí),捕獲到泥石流過流視頻,同時(shí)當(dāng)?shù)鼐用襁€通過手機(jī)現(xiàn)場拍攝了泥石流視頻。根據(jù)這些視頻資料,結(jié)合地形測繪資料,可計(jì)算出該次冰川泥石流的流速和峰值流量兩個(gè)關(guān)鍵動(dòng)力參數(shù)。流速計(jì)算公式如公式(1),計(jì)算結(jié)果詳見表3。

表3 出山口處泥石流流速觀測結(jié)果

式中:v為斷面流速(m/s);L為運(yùn)動(dòng)距離(m);T為運(yùn)動(dòng)時(shí)間(s)。

峰值流量基于流速和過流斷面進(jìn)行計(jì)算,計(jì)算公式如公式(2),計(jì)算結(jié)果詳見表4。

表4 出山口處泥石流峰值流量計(jì)算結(jié)果

式中:Q為峰值流量(m3/s);v為斷面流速(m/s);S為過流斷面(m2)。

以觀測所得流速對(duì)泥石流災(zāi)害防治工程勘查規(guī)范所推薦的黏性泥石流流速計(jì)算公式進(jìn)行對(duì)比分析分析結(jié)果見表5。

表5 流速經(jīng)驗(yàn)公式檢算

式中:v為斷面流速(m/s);n為泥石流粗糙系數(shù)公式(3)取0.45,公式(4)取0.27;H為泥石流流深(m);J為溝床縱比降。

計(jì)算結(jié)果分析表明,直白溝“9·10”冰川泥石流與暴雨型泥石流在動(dòng)力參數(shù)上具有顯著的差距,對(duì)比規(guī)范推薦流速計(jì)算公式計(jì)算,理論計(jì)算值與觀測值差距50%以上,由此按照泥痕法計(jì)算峰值流量具有同樣的誤差率。結(jié)合觸發(fā)的氣象條件,進(jìn)一步證明了冰川泥石流成因機(jī)制與動(dòng)力特征的差異性,傳統(tǒng)的泥石流動(dòng)力計(jì)算分析方法無法有效刻畫冰川泥石流運(yùn)動(dòng)特征。

4 泥石流形成機(jī)制討論

(1)泥石流驅(qū)動(dòng)力。該次泥石流事件暴發(fā)前7 d內(nèi)并無明顯降雨過程,氣溫也沒有明顯異常,由此可知泥石流中的主要水源并非來自降雨或者融水,這與傳統(tǒng)水力驅(qū)動(dòng)型泥石流存在巨大差異。結(jié)合前后遙感影像分析發(fā)現(xiàn),該次泥石流是由2#支溝出現(xiàn)的高位崩塌并鏟刮了斜坡上的冰川和冰磧臺(tái)地而演化形成。高寒山區(qū)的巖體穩(wěn)定性與凍融作用密切相關(guān)。直白溝泥石流暴發(fā)前7 d的平均最高溫度和最低溫度分別為21.5℃和13.1℃,根據(jù)海拔對(duì)氣溫影響規(guī)律,研究區(qū)相應(yīng)的0℃等溫線高程范圍為5 100～6 550 m。激發(fā)泥石流的高位巖崩海拔正處于0℃等溫線波動(dòng)范圍,這也從側(cè)面進(jìn)一步說明泥石流災(zāi)害是凍融循環(huán)作用的驅(qū)動(dòng)。同時(shí)2017年地震對(duì)周邊巖土體的震裂效應(yīng)也加劇了凍融循環(huán)作用的強(qiáng)度,未來周邊山體可能還存在類似的巖體失穩(wěn)隱患將是冰川泥石流的重要激發(fā)源。

(2)高流速特征。前述可知該次冰川泥石流出山口位置的龍頭最大流速達(dá)17.72 m/s,且遠(yuǎn)高于經(jīng)驗(yàn)公式計(jì)算的理論值,這說明該次泥石流災(zāi)害的能量轉(zhuǎn)化不同于傳統(tǒng)意義上的降雨型泥石流,其運(yùn)動(dòng)過程中的能量受到初始激發(fā)點(diǎn)物源的高勢能控制。啟動(dòng)該次泥石流的高位巖崩(海拔高度5 500 m)距離撞擊的斜坡相對(duì)高差約600 m,距離溝床底部(海拔高度4 225 m)達(dá)1 275 m。高勢能巖體失穩(wěn)形成的碎屑流具有極快的運(yùn)動(dòng)速度,同時(shí)冰雪的相變作用對(duì)流體輸移具有減阻等作用更進(jìn)一步減小了勢能向動(dòng)能轉(zhuǎn)化過程中的能量損耗,經(jīng)驗(yàn)公式中的阻力系數(shù)可能遠(yuǎn)低于規(guī)范中所界定的范圍。

(3)水文過程特征。在流速之外該次泥石流災(zāi)害的峰值流量同樣極高達(dá)到9 923.20 m3/s。按照流域匯流的水文理論,僅57.50 km2匯水面積的直白溝在無明顯降雨情況下,基本不可能形成滿足該次泥石流所需的水量,因此冰川泥石流的水源來源明顯不僅僅是降雨和冰雪融水。從災(zāi)害過程分析,該次泥石流的水體主要來源于巖崩撞擊斜坡上的冰川和冰磧平臺(tái)上裹挾的冰和雪,固態(tài)水體以集中補(bǔ)給方式參與泥石流活動(dòng)。泥石流龍頭峰值流量與龍身流量比值達(dá)15.67,流量過程呈現(xiàn)明顯三角形而非五邊形,這也說明存在集中水體補(bǔ)給的泥石流的沖出量計(jì)算不能按照五邊形概化模型進(jìn)行分析。

5 結(jié)論

(1)近70 a來直白溝冰川泥石流的驅(qū)動(dòng)力有冰川躍動(dòng)和高位冰崩兩種,均受到地震和氣候變暖的影響,尤其是地震誘發(fā)的同震災(zāi)害和震裂山體的演化。泥石流災(zāi)害的啟動(dòng)過程中的水熱優(yōu)勢條件組合并非必要條件。

(2)2020年“9·10”冰川泥石流災(zāi)害的演進(jìn)過程可分為巖崩—冰崩—冰巖碎屑流—泥石流—堵江,泥石流水源主要來自于冰雪為代表的固態(tài)水,相對(duì)于匯流形成的泥石流,水土耦合過程顯著異于降雨補(bǔ)給型泥石流。

(3)冰崩型泥石流災(zāi)害具有高初始勢能和集中補(bǔ)給水源的特點(diǎn),進(jìn)而導(dǎo)致其高流速和高峰值流量的運(yùn)動(dòng)學(xué)特點(diǎn),沿溝道運(yùn)動(dòng)過程中的沖刷和淤積特征主要受固體物質(zhì)控制,而非水動(dòng)力條件控制。

(4)冰崩型泥石流災(zāi)害的動(dòng)力學(xué)分析表明,冰水相變作用對(duì)流速和峰值流量均具有顯著放大效應(yīng)。現(xiàn)有規(guī)范推薦的流速、峰值流量和一次性沖出固體物質(zhì)量計(jì)算公式不適宜于此類災(zāi)害的計(jì)算,誤差在50%以上。