四川省典型區域滑坡泥石流與降水的關系

劉海知， 馬振峰， 范廣洲

(1.成都信息工程大學 大氣科學學院, 四川 成都 610225 2.四川省氣候中心， 四川 成都610072)

四川省典型區域滑坡泥石流與降水的關系

劉海知1， 馬振峰2， 范廣洲1

(1.成都信息工程大學 大氣科學學院, 四川 成都 610225 2.四川省氣候中心， 四川 成都610072)

[目的] 對四川典型區域滑坡泥石流與降水的關系進行研究，為進一步分析四川地區氣候變化與地質災害分布規律之間的關系提供基礎。 [方法] 四川滑坡、泥石流災害多由降水引發，在收集整理近15 a滑坡泥石流災害事件后，確定出四川滑坡、泥石流的典型區域，進而分析典型區域滑坡、泥石流與降水氣候特征的關系。 [結果] 四川省滑坡典型區域在川東地區， 泥石流典型區域在阿壩州地區。近年典型區域滑坡、泥石流近年頻次明顯增加，典型區域滑坡屬于較長日數降水誘發型，與15 d有效降水密切相關，前期降水中，強降水所占比例較大；典型區域泥石流屬于較短日數降水誘發型，與3 d有效降水密切相關，前期降水中，中小型降水所占比例較大。 [結論] 研究區15 d有效降水量和3 d有效降水與降水氣候特征間的關系密切，可分別作為典型區域滑坡、泥石流預報因子。

滑坡; 泥石流; 降水; 氣候特征

地處我國西南的四川地區，位于地貌第一階梯向第二階梯的過度帶上，南北大斷裂帶的鮮水河斷裂(龍門山斷裂)—安寧河斷裂—小江斷裂—紅河斷裂縱觀其間，巨大的地形起伏，復雜的地質構造，加上西南季風影響形成的獨特氣候，導致四川地區歷來是我國滑坡、泥石流的高發區。災害威脅著山區的交通干線、農田水利等基礎設施安全，成(都)昆(明)鐵路沿線就分布著3005條泥石流溝，占全國鐵路沿線的30%，為全國之最。如1997年6月5日，四川省美姑縣則租滑坡轉化為特大泥石流，造成四個村莊1527人受傷，151人死亡。許多學者展開了短期和短臨地質災害預報系統研究，相關科研機構和業務單位也發展了業務化運行系統。譚炳炎[1-3]根據成昆鐵路暴雨泥石流歷史數據，大致給出了泥石流爆發的雨量等級。李軍、周成虎[4]在分析香港地區滑坡大小與降水關系時，進行了多種周期的前期降水的定量化研究，發現不同類型的滑坡出現與前期降水有不同的關聯。國外研究中，Berruti在通過研究布雷西亞6—18世紀滑坡與洪水的關系時，建立了該省地質災害的歷史分布圖。Keefer等[5]依據土體滲透能力、含水量變化，采用靜力學理論推導降雨強度及延時與泥石流發生的關系。而目前針對7 d以上的較長時間地質災害預報研究比較缺乏，因而制約了地質災害中期以上預報研究和業務發展，難以滿足防災減災業務需求。因此，收集四川地區歷史上發生的地質災害案例記錄，統計四川地區滑坡、泥石流災害的年際分布規律和年內分布規律，為進一步分析四川地區氣候變化與地質災害分布規律間的關系提供基礎，對四川地區減災體系建設，謀求社會經濟的持續發展具有現實意義。

1 資料與方法

1.1 資料處理

(1) 災情資料。1999—2013年四川省滑坡、泥石流災害歷史資料，考慮到災情資料的完整性和客觀性，把1999年作為起始年份，15 a災情資料主要來自于中國科學院水利部成都山地災害與環境研

究所、四川省國土廳，少部分由四川省民政廳和四川省救災年鑒補充。統計時，以“縣”和“天”為單位，即：同一天某縣多地發生滑坡、泥石流災害記為該縣當天發生災情一次[6]。最終得到四川省1999—2013年時間精確到日，地點精確到縣的滑坡災害個例2746個，泥石流災害個例960個。

(2) 降水資料。四川省156個觀測站的逐日降水資料；來自于四川省氣候中心，剔除峨眉、新都、郫縣3個在1999—2013年間無滑坡、泥石流發生的縣的測站，最終得到153個站點的降水資料。災害點密度遠遠大于氣象站個數，確定災害點最近的氣象站，其降水作為該災害點的雨量。

1.2 研究方法

主要采用相關分析、線性回歸的方法[7]，確定四川滑坡泥石流典型區域，研究典型地區地質災害發生與降水氣候特征的關系，探求時空分布規律，分析典型區域災害—氣候要素耦合關系。

2 典型區域滑坡、泥石流災害分析

2.1 滑坡、泥石流典型區域

圖1分別為四川滑坡、泥石流多年平均分布圖，空心點表示氣象站位置。由圖1可以看出，滑坡災害多出現在四川東部，高頻區在巴山地區，巴中、達州部分地區年均滑坡發生頻次在4次以上；泥石流災害主要發生于四川西部，高頻區在川西以北的阿壩地區，阿壩州部分縣年均泥石流爆發頻次達到3次以上。因此將滑坡年均發生頻次在3次及其以上的廣元(市區)、南江縣、萬源縣、儀隴縣、宣漢縣、達州(市區)、南充(市區)7個區縣作為四川滑坡典型區域，將泥石流年均發生頻次在2次及其以上的壤塘縣、金川縣、馬爾康縣、汶川縣、理縣、黑水縣6個縣作為四川泥石流典型區域。

圖1 四川省地區滑坡、泥石流多年平均分布圖

2.2 滑坡、泥石流災害時空變化

圖2為典型區域滑坡、泥石流災害時空分布，滑坡頻次在2003年出現極大值，泥石流頻次在2010年出現極大值。兩類災害數量近年都明顯增加，一是因為該地區人類活動的增強，二是由于地質災害統計工作的更加客觀全面。兩類災害在月際分布上大致呈現正態分布，都主要集中在夏季爆發，夏季發生的次數明顯多于其他季節，最高值月份出現在7月，與降雨在月際分布上吻合。在空間分布上，滑坡災害頻次在前三的區縣依次是達州(市區)、宣漢縣、南江縣，其中達州(市區)滑坡爆發頻次遠遠多于其他代表區縣。泥石流發生頻次在前三的區縣依次為理縣、汶川縣和金川縣，各個區縣的災害頻次相差不大。

圖2 典型區域滑坡泥石流災害時空分布

3 降水誘發滑坡、泥石流災害分析

3.1 滑坡、泥石流災害與前期降水的關系

為探究前期累計雨量與滑坡、泥石流發生的關系，以典型區域近6年中(2008—2013年)發生滑坡、泥石流災害次數與其前期降雨量因子做相關性分析，結合災害滯后于降雨的時間，選定前期累計雨量有效時段為15 d。各站點所在區域氣候條件不同；地表徑流的產生、水分的蒸發等過程，使得進入巖土體的雨量小于實際記錄雨量；且一次降雨不一定會導致地質災害發生，每次降雨中也只有部分降雨量對災害的發生起作用，因此用有效降雨量來表示前期累計雨[8-9]，前期有效降雨量是指災害發生前的降雨過程中對巖土體穩定性產生作用的雨量[9]，地形陡峭，土壤吸水能力較好，有效雨量則較大，反之則較小。以災害典型區域各區縣的逐日降水量計算其各次災害前期的有效降水量﹝公式(1)﹞，計算得出滑坡、泥石流災害與不同累計降雨量的相關系數(表1)。因為滑坡、泥石流災害典型區域地形地貌條件復雜，前期有效降水量對各區縣地質災害的影響程度不完全一致[11]，需要分別確定各區縣的最佳有效降水系數﹝如公式(2)﹞。

(1)

式中：α——有效降水系數(0<α<1)，表示巖土體對雨水滯留能力的大小，它由區域內巖土體總體性質決定[10]；Rc——有效降雨量(mm)；R0——災害發生當天的降雨量(mm)；Rn——災害發生前n天的降水量(mm)；α——典型區縣的最佳有效降水系數。

(2)

式中：σ——該處有效降雨量標準差。

由表1可知，典型區域滑坡與前15d有效降雨量相關系數最大，后續依次是前10d有效雨量和前1d有效雨量；可見誘發滑坡的動態因子多為長時間降水；典型區域泥石流與前3d有效降雨量相關系數最大，后續依次是前2d有效雨量和前1d有效雨量，可以看出泥石流由相對較短時間降水激發。

3.2 不同級別降水與滑坡、泥石流的關系

分別選取與滑坡相關性最高的前15d有效雨量和與泥石流相關性最高的前3d有效雨量作為分析對象，按照每25mm一個量級對有滑坡前15d的進行分組，每5mm一個量級對有泥石流前3d的進行分組，得到兩類災害頻次在各個量級有效雨量的分布(圖3)。

由圖3可以看出，在兩類災害在各個有效降雨量級的分布上都比較離散，沒有明顯規律性；即在很高雨量和很低雨量的情況下均有滑坡、泥石流災害發生；這符合實際，因為降雨量很少的條件下，若地表水沿優勢通道滲入地下，也能誘發滑坡、泥石流[12]。

表1 典型區域滑坡、泥石流發生前15 d內 降雨量因子與災害次數的相關性分析

滑坡次數泥石流次數X44MM?時間/d 相關系數 排序X40MM?相關系數 排序 0123456789101112131415074807820751075607630768076807720773077507830778077907810780078816315141311111098276451084908530854088908400837083008270824082108230822082608250824082343215678111615149101113

為避免造成有效雨量越大，災害爆發率越低的直觀誤解，圖3分別給出了滑坡、泥石流頻次以及頻次累計百分比與前期有效降降雨的關系；典型區域滑坡頻次累計百分比25%,50%,75%對應的前15 d有效降雨量級分別是[25 mm，50 mm),[100 mm，125 mm),[150 mm，175 mm)，即前15 d有效降雨達

到[25 mm，50 mm)時已有25%的滑坡發生；達到[100 mm，125 mm)時已有50%的滑坡發生；達到[150 mm，175 mm)時已有75%的滑坡發生。典型區域泥石流頻次累計百分比25%，50%，75%對應的前3 d有效降雨量級分別是[10 mm，15 mm)，[20 mm，25 mm)，[30 mm，35 mm)，即前3 d有效雨量達到[10 mm，15 mm)時已有25%的泥石流發生；達到[20 mm，25 mm)時已有50%的泥石流發生；達到[30 mm，35 mm)時已有75%的泥石流發生。

將滑坡、泥石流前期降雨中有1 d及其以上日降水量≥50 mm(暴雨)，25 mm<日降水量≤50 mm(大雨)和日降水量<25 mm(中、小雨)的災情個例分為3組，[25 mm,50 mm)個例組中不包含≥50 mm的情況，<25 mm個例組中不包含≥50 mm和[25 mm,50 mm)的情況。從典型區域滑坡前15 d降水情況來看(表2)，日降水量≥50 mm條件下發生災害的次數遠遠多于其他雨量下的災害次數，前15 d平均有效降水為226.7 mm；典型區域泥石流前3 d降水中日降水量<25 mm條件下誘發災害的次數明顯大于其他雨量下的災害次數，日降水量≥50 mm條件下誘發泥石流個例數為0，前3 d平均有效降水為26.7 mm。這與實況相符，阿壩州大部分地區極少出現日降雨量≥50 mm的強降水天氣，滑坡典型區川東多低山丘陵，水汽條件好，地形抬升作用對強降水天氣系統(如西南渦)有加強作用，易造成強降水。

表2 滑坡泥石流發生前期降水情況

4 滑坡、泥石流與降水氣候特征的關系

滑坡、泥石流災害多由異常降水引發，為探究異常降水與多年平均降水的關系，利用典型區域269個滑坡個例和175個泥石流個例的降水特征數據，分析滑坡發生前15 d有效降水量、泥石流發生前3 d有效降水量與多年平均降水量、雨季降水量間的關系。分析結果顯示，滑坡前15 d有效降雨量與年降雨量、雨季降雨量線性相關；泥石流發生前3 d有效雨量與年降雨量無明顯相關性，與雨季降水量線性相關，相關性用函數關系式描述詳見表3—4。在對滑坡的前15 d有效降水和泥石流的前3 d有效降水分別與年暴雨雨量、年大雨雨量間關系進行分析后得出，15 d有效雨量與年暴雨雨量線性相關，與大雨以及以下的雨量不存在相關性；誘發泥石流前3 d雨量與大雨雨量存在一定線性關系，與其他雨量不具有明顯關系，相關性用函數關系式描述詳見表3—4。根據以上分析結果可見，滑坡前15 d降水與區域降水氣候背景間表現出較強規律性，可作為典型區域滑坡預報因子，泥石流前3 d降水同區域降水氣候背景間有較強相關，可作為典型區域泥石流預報因子。

表3 滑坡前15 d雨量擬合結果

圖3 典型區域滑坡泥石流與前期降水關系

降水函數類型R2常數項一次項雨季降水量線性031033570035大雨雨量線性0340194250126

5 結 論

(1) 四川滑坡典型區域在川東地區；廣元(市區)、南江縣、萬源縣、儀隴縣、宣漢縣、達州(市區)、南充(市區)7個區縣年均滑坡爆發頻次在3次及其以上；四川泥石流典型區域在阿壩州地區，壤塘縣、金川縣、馬爾康縣、汶川縣、理縣、黑水縣年均泥石流爆發頻次達到2次及其以上。

(2) 典型區域滑坡、泥石流近年頻次明顯增加，主要在夏季爆發。

(3) 典型區域滑坡與前15 d有效雨量相關性最大，屬于長時間降水誘發型，典型區域滑坡發生的前期降水中強降水所占比例較大；典型區域泥石流與前3 d有效雨量相關性最大，屬于較短時間降水誘發型，典型區域泥石流發生的前期降水當中中小型降水所占比例較大。

(4) 與滑坡發生密切相關的15 d降水量與年降水量、雨季降水量、暴雨雨量線性相關，與小于50 mm雨量線性無關，與泥石流發生密切相關的3 d降水量與雨季降水量、大雨雨量線性相關，與年降水量和暴雨雨量線性無關。

[1] 譚炳炎,楊大文,石勝國.暴雨泥石流預報的研究[J].鐵道學報,1992(3):92-101.

[2] 譚炳炎,段愛英.山區鐵路沿線暴雨泥石流預報的研究[J].中國鐵道科學,1994(4):67-78.

[3] 譚炳炎.二十世紀中國鐵路沿線泥石流防治理論與實踐[J].鐵道工程學報,2005(S1):369-372.

[4] 李軍,周成虎.香港地區滑坡體積與前期降水關系分析[J].自然災害學報,2002,11(2):37-45.

[5] Keefer K D, Wilson C R, Mark K R, et al. Real-time landslide warning during heavy rainfall [J]. Science,1987,238(4829)：921-924.

[6] 唐川,朱靜.云南滑坡泥石流研究[M].北京:商務印書館,2003.

[7] 黃嘉佑.氣象統計分析與預報方法[M].北京:氣象出版社,2000.

[8] 謝劍明,劉禮領,殷坤龍,等.浙江省滑坡災害預警預報的降雨閥值研究[J].地質科技情報,2003,22(4):101-105.

[9] 張國平.有效雨量和滑坡泥石流災害概率模型[J].氣象,2014，40(7),886-890.

[10] 叢威青,潘懋,李鐵鋒,等.降雨型泥石流臨界雨量定量分析[J].巖石力學與工程學報,2006,25(S1):2808-2812.

[11] 胡娟,閔穎,李華宏,等.云南省山洪地質災害氣象預報預警方法研究[J].災害學,2014,29(1):62-66.

[12] 惠鑫.陜北黃土高原地質災害與降雨關系研究[D].陜西 西安：長安大學,2010.

Relationship Between Landslide/Debris Flow and Rainfall in Typical Region of Sichuan Province

LIU Haizhi1， MA Zhenfeng2， FAN Guangzhou1

(1.SchoolofAtmosphericSciences,ChengduUniversityofInformationTechnology,Chengdu,Sichuan610225,China; 2.SichuanClimatCenter,Chengdu,Sichuan610072,China)

[Objective] We studied the relationship between climatic characteristics of rainfalls and geological disasters in typical regions of Sichuan Province to provide basis for analyzing the relation between climatic variation and geological disasters distribution in Sichuan Province. [Methods] Disaster cases were collected in recent 15 years and the typical regions of landslide and debris flow in Sichuan Province were marked out to analyze the relationship between geological disasters and the climatic characteristics of rainfalls. [Results] Landslides mainly occurred in the east of Sichuan Province, while debris flows mainly occurred in the west of Sichuan Province. The frequencies of the two kinds of disasters increased in recent years. Landslides in typical regions were caused by days precipitation, and specifically 15 days effective precipitation was responsible for it. Beside that, heavy rains were ordinarily found in early stage of the collected landslide cases. Debris flows in typical regions were caused by few-day rainfall, and was closely related to 3 days effective precipitation. Moderate rain and light rain were frequently observed in 3 days effective precipitation before debris flow outbreaks. [Conclusion] 15 days effective precipitation and 3 days effective precipitation were both associated closely to climate characteristics of rainfalls: the former can be the factor for landslide forecast, and the latter can be the factor for debris flow forecast.

landslide; debris flow; rainfall; climatic characteristics

2016-05-05

2016-05-25

國家自然科學基金項目“華西秋雨年際變化規律及其形成機理研究”(41275097)； 四川省科技支撐計劃項目(15C0098)

劉海知(1991—)，男(漢族)，四川省榮縣人，碩士研究生，研究方向為短期氣候診斷。E-mail:lhz1012@aliyun.com。

馬振峰(1961—)，男(漢族)，內蒙古自治區四子王旗人，碩士，研究員，碩士生導師，主要從事氣候預測，氣候變化及生態環境研究。E-mail:mzf616@aliyun.com。

10.13961/j.cnki.stbctb.2016.06.013

A

1000-288X(2016)06-0073-05

X43, P64

文獻參數： 劉海知， 馬振峰， 范廣洲, 等.四川省典型區域滑坡泥石流與降水的關系[J].水土保持通報，2016,36(6)：073-077.