江金濤，鐵永波

（1.中國地質(zhì)大學（北京），北京100083；2.中國地質(zhì)科學院，北京100037；3.中國地質(zhì)調(diào)查局成都地質(zhì)調(diào)查中心，四川成都610081）

四川安寧河流域姑巴溝泥石流動力學特征與減災對策*

江金濤1，2，3，鐵永波3

（1.中國地質(zhì)大學（北京），北京100083；2.中國地質(zhì)科學院，北京100037；3.中國地質(zhì)調(diào)查局成都地質(zhì)調(diào)查中心，四川成都610081）

姑巴溝是四川安寧河流域內(nèi)最具代表性的典型低頻泥石流之一?；谝巴庹{(diào)查的基礎上，對姑巴溝泥石流的動力學特征進行了計算，并對2011年泥石流的成因機制進行了分析。研究認為，姑巴溝泥石流的觸發(fā)原因是由持續(xù)強降雨引發(fā)的溝道徑流與上游冰湖溢水的綜合作用，水動力補給過程具有復合特征。針對姑巴溝泥石流爆發(fā)突然、流體流量大流速高及其威脅對象的分布特征，提出了針對性的防災減災建議，對揭示姑巴溝泥石流的形成過程及動力學特征具有實際意義，并可為安寧河流域內(nèi)其它具有相似孕災背景的低頻泥石流形成機理的認識提供參考。

安寧河流域；低頻泥石流；動力學特征；冕寧縣；姑巴溝

泥石流是山區(qū)常見的一種突發(fā)性地質(zhì)災害，按照暴發(fā)頻率可分為高頻、中頻和低頻泥石流。低頻泥石流暴發(fā)周期長（30年以上），規(guī)模大，危害重，具有強隱蔽性［1］。溝口附近居民往往將此類泥石流平緩的堆積扇作為安全地帶，修建人口密集的村舍、工廠和學校等［2］。較高頻泥石流而言，低頻泥石流致災性更強。1999年12月在委內(nèi)瑞拉北部阿維拉山區(qū)加勒比海沿岸的8個州暴發(fā)的約300年一遇的群發(fā)性泥石流，造成3萬余人死亡，直接經(jīng)濟損失達100億美元［3］；2010年8月在甘南藏族自治州舟曲縣城發(fā)生的400年一遇的低頻泥石流導致該縣一半被淹，造成1 744人死亡和失蹤，直接經(jīng)濟損失數(shù)億元［4］；2012年6月在四川省寧南縣矮子溝發(fā)生的特大泥石流，導致下游基礎設施損失慘重，造成40人死亡和失蹤，產(chǎn)生了極大的社會影響［5］。

研究表明，四川安寧河流域5個縣市（喜德、冕寧、西昌、會理、德昌）均有大量泥石流溝分布，僅冕寧縣內(nèi)就有泥石流溝143條，泥石流溝密度高達32.1條／103km2［6-7］。泥石流活動成為該區(qū)最常見的地質(zhì)災害之一，常常造成家毀人亡等重大生命和財產(chǎn)損失，極其不利于當?shù)氐姆€(wěn)定與發(fā)展。本次研究以具有典型代表性的姑巴溝低頻泥石流為例，通過查明其動力學特征及分析形成機制，對認識安寧河流域內(nèi)其他泥石流災害的成因及防災減災提供科學依據(jù)。

1 流域基本概況

姑巴溝位于四川省涼山彝族自治州北部的冕寧縣彝海鄉(xiāng)境內(nèi)，安寧河東，距冕寧縣城區(qū)NNE方向約16 km。流域面積約66.7 km2，主溝長10.5 km，縱比降202‰。主溝北側發(fā)育有一條較大的支溝（小溝），其余支溝的流域面積均相對較小。流域內(nèi)最高海拔4 630 m，為安寧河東側分水嶺—小相嶺，常年積雪；溝口最低海拔2 102 m，相對高差2 528 m。根據(jù)其地貌及泥石流物源補給特征，將姑巴溝劃分為清水區(qū)、形成區(qū)、流通區(qū)和堆積區(qū)四個部分。清水區(qū)海拔3 249～4 630 m，呈喇叭狀，利于匯水形成表面徑流；形成區(qū)海拔2 533～3 249 m，V型谷，溝岸坡度約40°，坡面泥石流廣泛發(fā)育；流通區(qū)海拔2 288～2 533 m，呈V型谷，溝岸坡度35°，溝道平直，利于流泥石流快速沖出；堆積區(qū)海拔2 102～2 288 m，縱比降46‰，呈U型谷，溝床寬數(shù)百米，泥石流流體物質(zhì)到達堆積區(qū)流速逐減，開始堆積。如圖1所示。

2011年6月16日20時左右，四川省冕寧縣彝海鄉(xiāng)遭受3 h強降雨襲擊，總降雨量101 mm，并引發(fā)泥石流［8］。暴發(fā)泥石流之前，小溝先漲水，大溝緊隨其后。主溝兩岸受災情況迥異：左岸布落2人死亡，15人失蹤；右岸布里僅受洪水淹沒，漫水高出地面1 m。此次泥石流造成當?shù)亟煌ê屯ㄐ胖袛啵? 890人受災。泥石流發(fā)生后，當?shù)卣嚓P部門開展了溝道清淤及簡易的河道治理工程，但這些工程難以對其潛在風險進行有效控制，為此，查明該溝泥石流的形成過程及機制，并提出有針對性的防災減災措施，最大程度降低泥石流帶來的潛在風險極為迫切。

2 泥石流成因及流動力學特征

2.1 姑巴溝泥石流成因分析

以2011年姑巴溝泥石流過程為例，主要從泥石流形成的三個條件：地形、水動力及物源，對該溝泥石流的成因進行分析。

2.1.1 地形特征

姑巴溝支溝發(fā)育，相對高差大，平面上呈圈椅狀，清水區(qū)和形成區(qū)縱寬達8 km，溝口僅數(shù)百米，這一形態(tài)極利于地表徑流匯集。溝道較平直，上游呈“V”型谷，跌水和陡坎發(fā)育，下游呈“U”型谷。在流域的頂端，位于主溝和小溝的源頭各有面積約0.1 km2的冰磧湖。

圖1 姑巴溝流域形態(tài)圖

2.1.2 水動力補給特征

流域內(nèi)具有冕寧縣所處區(qū)域上的典型亞熱帶季風氣候，降雨集中在6-9月份，常年平均降雨量1 094 mm。16日20時冕寧縣突發(fā)強降雨，3 h降雨量達101 mm。降雨過程中上游冰磧湖水體溢出，具有相對較大的勢能，對溝道內(nèi)松散土體的沖刷作用較強。降雨形成的表面徑流填充土體孔隙，使之浸潤、飽和，并逐漸加大飽和土厚度，土體隨后起動奔流。強降雨是激發(fā)此次泥石流的主要因素，冰湖溢水在泥石流的形成中有重要作用。

2.1.3 物源補給特征

姑巴溝以及整個安寧河流域受安寧河深斷裂帶影響明顯。斷裂帶兩側節(jié)理裂隙密度達30～60條／m2，巖體破碎，極易風化。姑巴溝內(nèi)巖性復雜，涵蓋了從震旦紀以來的3大巖類和第四紀冰磧物，強風化帶厚2～25 m［7，9］。調(diào)查發(fā)現(xiàn)海拔范圍2 533～3 249 m溝道坡面第四紀殘坡積物較厚，是本次泥石流災害的啟動物源；2 288～2 533 m之間的溝床內(nèi)第四紀沖洪積物和冰磧物受到了溝道徑流揭底侵蝕后，溝谷中堆積物在極短的時間內(nèi)被席卷。流體攜帶的固體物質(zhì)如同“滾雪球”一般快速增大，形成這次大規(guī)模泥石流。

按動力條件劃分，2011年姑巴溝泥石流屬于水力類泥石流［10］。

2.2 動力學參數(shù)特征

2.2.1 流體流速

泥石流流速是泥石流的重要屬性，在評價泥石流危險性過程中以及防治工程中有重大意義。姑巴溝泥石流土體粘粒含量約5%，產(chǎn)流過程中，與入滲為主的水結合，屬粘性泥石流。按照王繼康等在1982年提出的粘性泥石流流速經(jīng)驗公式，推導姑巴溝泥石流的流速［11］：

式中：Vc為粘性泥石流流速，m／s；Hc為泥石流泥深（或泥石流水力半徑Rc），m；Ic為泥石流流面（或溝底）坡度，‰；Kc為粘性泥石流流速系數(shù)。

本次調(diào)查中，一共實測姑巴溝中下游不同部位的7個剖面（P1-P7）。在這些剖面的溝岸，能夠清楚地分辨出泥石流發(fā)生時留下的泥位。通過實測剖面數(shù)據(jù)，可以獲得泥石流泥深Hc和剖面處泥石流流面坡度Ic。

各剖面處泥石流流速見表1。剖面P1～P3、P5～P6位于形成區(qū)末端，P7位于流通區(qū)頂端，P4位于堆積區(qū)頂端。各剖面流體流速差異較大，堆積區(qū)流體流速降低明顯，此時固體物質(zhì)以淤積為主。受水電站的影響，位于剖面P5下游約500 m的剖面P7，流速相對降低3.85 m／s。

表1 剖面流速

2.2.2 流量

泥石流流量計算方法：可分為雨洪修正法與形態(tài)調(diào)查法兩種。

（1）雨洪修正法

雨洪修正法認為泥石流流量為雨洪流量（Qp）和固體流量（QH）之和，自1959年以來，在驗證各家泥石流流量計算公式的基礎上，考慮粘性泥石流有堵塞因素的實際，將公式改為［11］：

式中：Qc為泥石流流量，m3／s；Qp為清水流量，m3／s。用泥石流溝所在地區(qū)的小徑流公式計算清水流量；Dc為泥石流堵塞系數(shù)。粘性泥石流有陣流堵塞，而稀性泥石流沒有；φ為泥石流修正系數(shù)。計算法φ＝（ρc-1）／（ρH-ρc）（ρH由查表獲得［11］）。雨洪法計算泥石流流量見表2。

表2 雨洪法計算泥石流流量

（2）形態(tài)調(diào)查法

調(diào)查泥石流泥痕和溝槽形態(tài)斷面特征，根據(jù)流速（νC）的計算結果，與泥石流通過斷面面積（WC）相乘。按下式計算泥石流流量Qc：

在計算過程中，分別選取了小溝和主溝上的泥石流剖面P1、P7（表3）。

表3 形態(tài)調(diào)查法計算泥石流流量

流域內(nèi)泥石流峰值流量Qc應P1和P7剖面流量之和，最后計算出泥石流流量Qc＝2 929.75 m3／s。計算結果與雨洪法在頻率p＝1%時得出的計算結果2 736.16 m3／s基本吻合。暴雨發(fā)生頻率與對當?shù)卮迕癫稍L結果相仿。故認為2011年6月姑巴溝發(fā)生流量2 929.75 m3／s的特大型泥石流為100年一遇。

2.2.3 泥石流沖壓力

泥石流沖擊力是泥石流防治工程設計的重要參數(shù)。

（1）泥石流整體沖壓力

采用推薦的泥石流體整體沖壓力計算公式［12］：

式中：δ為泥石流整體沖壓力，kPa；γc為泥石流重度，姑巴溝泥石流取1.544（kN／m3）；Vc為泥石流流速，取舍去剖面P5數(shù)據(jù)的平均值7.368 m／s；α為建筑物受力面與泥石流沖擊壓力方向的夾角，取90°；g為重力加速度，取9.8 m／s2；λ為建筑物形狀系數(shù)，圓形建筑物取1.0，矩形建筑物取1.33，方形建筑物取1.47。本次布落、布里建筑物主要以矩形為主，λ取1.33。

通過計算，姑巴溝泥石流整體沖壓力為111.48 kPa。

（2）泥石流中大石塊沖擊力

泥石流中石塊沖擊力的計算參照以下公式［13］：

式中：Pd為為泥石流中石塊的沖擊力，Pa；γ為動能折減系數(shù)，對于圓端正面撞擊，采用γ＝0.3；Vc為泥石流大石塊流速，為4.56 m／s；Q為石塊質(zhì)量，kg，按離溝口最大粒徑1.3 m來算；α為受力面與泥石流沖擊力方向的夾角，取90°；C1、C2為巨石及攔擋圬工的彈性變形系數(shù)，C1＋C2＝0.005。

姑巴溝泥石流石塊比重為2.59 t／m3，近似認為石塊與被撞擊物之間的接觸面積為圓面，最后得到姑巴溝泥石流受災體附近大石塊的沖擊力為1.13 kPa。

2.2.4 泥石流爬高

泥石流在運動過程中，若突然受阻會產(chǎn)生沖起爬高，如橋墩、擋壩、陡壁等。由以下公式可以求得沖起爬高值［12］：

式中：ΔZ為泥石流沖起爬高值，m；Vc為泥石流流速，m／s；g為自由落體加速度，m／s2；α為動能修正系數(shù)，為迎面坡度的函數(shù)（迎面坡度90°，α取1.0）。α取1.0作為近似計算。

小溝在布里村東并流進主溝，流向迅速南偏約30°。位于凹岸的布里村受泥石流固體物質(zhì)沖擊較大，村民的生命財產(chǎn)安全受到嚴重威脅。選擇距布里村最近的P4剖面流速5.68 m／s，計算出泥石流最大沖爬高為2.63 m，高于布里村地面線約1 m。

3 泥石流的防災減災建議

針對姑巴溝泥石流爆發(fā)突然、流體流量大流速高和基礎設施分布特殊的特點，對該溝提出可行性防災減災建議?？梢宰鳛榱饔騼?nèi)其他具有相似特征的泥石流的防災參考。

3.1 排導工程措施建議

小溝在匯流進主溝之前，溝口段1 km向左偏轉(zhuǎn)約30°形成彎曲溝道，布里村臨近被侵蝕的凹岸。如果再暴發(fā)同等規(guī)模泥石流，必定會危及布里村。因此建議在彎道處布置長1 km、高1.7 m的V型全襯砌槽，使進入排導槽流體保持密度保持不變，利用水動力增大水深，重組動力束流，減小阻力，加大流速，防止漫流改道。并在布里村溝岸一側輔以擋墻，以抵擋流體對凹岸的侵蝕。依照泥石流最大沖爬高計算結果，擋墻應該高出地面至少2 m。

3.2 導流工程措施建議

布落村地勢較低，正好處于小溝溝口軸線上，是本次泥石流災害的重災區(qū)。該村靠近溝槽一側有前后兩次修葺的導流堤，所用材料均為溝道內(nèi)泥石流橢圓形漂石。舊導流堤高約0.6 m，堤角近70°，與溝床流向平行，絕大部分在這次特大型泥石流中被沖毀或掩埋，本次受災慘重與舊導流堤高度不夠有關?，F(xiàn)有導流堤于泥石流發(fā)生之后修建，高約1.5 m，堤邊坡坡度50°，軸線基本與主溝溝槽走向一致，距布落村直線距離10 m。經(jīng)計算，在保持導流堤邊坡坡度50°情況下，為了有效防治同等規(guī)模的泥石流災害，導流堤坡高應達到2.48 m?，F(xiàn)有導流堤的高度和走向?qū)τ诜罏臏p災均不理想。而在增加導流堤邊坡坡度到70°情況下，2 m高導流堤能防止流體爬越。因此應該在現(xiàn)有導流堤外圍，沿布落村修建一條向溝口發(fā)散長1.8 km的導流堤，迎水面邊坡坡度70°，堤高2 m。

3.3 監(jiān)測預警措施建議

建議姑巴溝建立高程控制的雨量計監(jiān)測預報和專業(yè)泥位監(jiān)測預警系統(tǒng)。在海拔2 500 m和3 500 m溝谷處建立雨量計，用確定的降雨臨界值指標，預報未來一段時間可能發(fā)生泥石流災害；在海拔2 800 m處安裝超聲波探頭實時監(jiān)控泥位高度，當超過危險值時，就可以進行警報。

4 結論

姑巴溝是安寧河流域內(nèi)典型的低頻泥石流溝，其致災性極強。2011年6月暴發(fā)的泥石流為100年一遇，粘性泥石流，泥石流的水動力來自3 h強降雨形成的地表徑流與流域上游冰磧湖水體混合形成溝道徑流。泥石流的物源主要來自溝床內(nèi)堆積的冰磧物，在被徑流沖刷后啟動形成泥石流。

姑巴溝泥石流嚴重威脅到下游及溝口的小型電站、數(shù)十戶房屋和高速公路等基礎設施，對當?shù)厝嗣袢罕娚敭a(chǎn)及基礎設施構成極大威脅，其潛在風險較高，需要對該溝泥石流開展工程治理，并輔以監(jiān)測預警設施，最大程度降低泥石流的潛在風險。

四川省安寧河流域地質(zhì)構造活動強烈，巖體破碎，加之該地區(qū)早期冰川活動遺留的大量松散冰磧物在流域上游廣泛分布，在河流侵蝕作用下穩(wěn)定性較差，為這一地區(qū)泥石流的形成提供了必要的地形和物源條件，使得這一區(qū)域內(nèi)泥石流極為發(fā)育。探索這一區(qū)域內(nèi)泥石流的成因及動力學特征，對該地區(qū)泥石流的防災減災具有防災減災意義。

［1］ 唐邦興，周必凡，吳積善，等.中國泥石流［M］.北京：商務印書館，2000：60-71.

［2］ 陳曉清，崔鵬，韋方強.良好植被區(qū)泥石流防治初探［J］.山地學報，2006，24（3）：333-339.

［3］ 韋方強，謝洪，Jose L Lopez，等.委內(nèi)瑞拉1999年特大泥石流災害［J］.山地學報，2000，18（6）：580-582.

［4］ 余斌，楊永紅，蘇永超，等.甘肅省舟曲“8.7”特大泥石流調(diào)查研究［J］.工程地質(zhì)學報，2010，18（4）：437-444.

［5］ 鐵永波，胡凱衡.基于遙感解譯的典型低頻泥石流形成機制研究—以四川省寧南縣矮子溝泥石流為例［J］.災害學，2014，29（3）：77-80.

［6］ 李秀珍，劉希林，蘇鵬程.四川涼山州安寧河流域泥石流危險性評價［J］.防災減災工程學報，2005，25（4）：426-430.

［7］ 嚴光輝，李健生，羅世清.四川省冕寧縣安寧河流域泥石流［J］.山地研究，1991，9（2）：123-130.

［8］ 李志杰，蘇冠群，朱光榮，等.四川山洪災害防治非工程措施建設顯作用［EB／OL］.［2015-03-27］.http：／／www.chinawater.com.cn／newscenter／slyw／201106／t20110623＿152680.html.

［9］ 趙良政.川西-冕寧瀘沽地區(qū)早更新世雙橋冰期之證據(jù)［J］.地球科學-中國地質(zhì)大學學報，1990，15（5）：495-500.

［10］唐川，章書成.水力類泥石流起動機理與預報研究進展與方向［J］.地球科學進展，2008，23（8）：787-793.

［11］王繼康，黃榮鑒，丁秀燕.泥石流防治工程技術［M］.北京：中國鐵道出版社，1996：37-84.

［12］國土資源部.DZ／T0239-2004泥石流災害防治工程勘察規(guī)范［S］.北京：中國標準出版社，2006.

［13］吳積善.泥石流及其綜合治理［M］.北京：科學出版社，1993：164-191.

Dynam ic Characteristics and M itigation of Guba Debris Flow in Anning Valley in Sichuan Province

Jiang Jintao1，2，3and Tie Yongbo3
（1.China University of Geosciences（Beijing），Beijing 100083，China；2.Chinese Academy of Geological Sciences，Beijing 100037，China；3.Chengdu Center of China Geological Survey，Chengdu 610081，China）

Guba debris flow is one of themost typical debris flows with low return period in Anning Valley，Sichuan Province.Based on the field survey，we calculate the dynamic properties and analyze the genetic mechanism of Guba debris flow occurred in 2011.Results show that，combined effectof gully runoff from sustained heavy rainfall and overflow from glacial lakes upstream triggered the disaster.The water power supply processes show complex features.According to the characteristics of sudden outbreak，heavy flow，high speed and particular position of the affected targets in Guba debris flow，proposals are put forward for prevention and reduction.Ithas practical significance to reveal forming process and dynamic properties of Guba debris flow，which could also be a reference to understanding the formation mechanism of low frequency debris flow with similar backgrounds in Anning Valley.

Anning Valley；debris flow with low return period；dynamic properties；Mianning County；Guba gully

X43；P642.23

A

1000-811X（2015）04-0216-04

10.3969／j.issn.1000-811X.2015.04.037

江金濤，鐵永波.四川安寧河流域姑巴溝泥石流動力學特征與減災對策［J］.災害學，2015，30（4）：216-219.［Jiang Jintao and Tie Yongbo.Dynamic characteristics and mitigation of Guba debris flow in Anning Valley，Sichuan Province［J］.Journal of Catastrophology，2015，30（4）：216-219.］

2015-03-27

2015-05-22

科技部基礎性工作專項（2011FY110100-5）；中國地質(zhì)調(diào)查局公益性地調(diào)項目（12120113010200）

江金濤（1991-），男，湖北新洲人，碩士研究生，主要從事環(huán)境地質(zhì)和災害地貌方向研究.

E-mail：jiang564487079＠126.com