四川省涼山州美姑縣瓦扎乃曲泥石流特征值分析與防治措施

劉盛達，張光輝，李彥彬，李成志

（四川省煤田地質局一三五隊 四川 瀘州 646000）

瓦扎乃曲溝泥石流位于四川省涼山州美姑縣巴古鄉境內，屬于美姑河左岸一級支流，為暴雨溝谷型泥石流。目前，瓦扎乃曲溝域內存在多處崩塌、滑塌以及人工采石棄渣堆現象，為瓦扎乃曲泥石流的暴發提供了較好的物源條件。在雨季期間容易誘發該泥石流的暴發，瓦扎乃曲泥石流一旦再次暴發將直接威脅溝口堆積區巴古鄉牛洛村三組23戶127人及牛洛小學師生120人的生命財產安全，以及區內S103省道的正常運營。因此，研究瓦扎乃曲溝泥石流的基本特征、形成條件、動力學參數，預測暴發規模，提出相應的防治方案十分必要。

1 瓦扎乃曲泥石流流域基本特征

瓦扎乃曲溝屬于美姑河左岸一級支流，流域面積約2.36km2，整體呈南-北向展布，縱向溝長約3.4km，平均縱向坡降約 259.35‰，流域最高點位于后緣山頂分水嶺，高程約 2600m，最低高程為瓦扎乃曲溝溝口與美姑河階地交匯處，高程約1670m，相對高程約930m。區內山高坡陡，具有一條主溝和三條支溝，1#支溝長約800m，與瓦扎乃曲溝中游左岸交匯，溝道較狹窄，平均縱坡降約447.84‰，溝道平均寬約2.5m，水流量較小，支溝流域內植被稀疏，物源豐富，頂部呈漏斗狀，匯水條件較好。2#支溝位于瓦扎乃曲溝上游右岸形成流通區與清水區交界，長約735m溝道平均寬約3.5m，平均縱坡降約417.16‰，頂部植被茂盛，匯水條件好。3#支溝位于瓦扎乃曲溝上游左岸形成流通區與清水區交界，長約985m溝道平均寬約4m，平均縱坡降約321.6‰，頂部植被茂盛，匯水條件好。

瓦扎乃曲溝整個流域形態分布較為對稱。兩側岸坡植被較茂盛，岸坡縱向坡度較陡一般 30°～50°，最大可達80°，這種地形條件有利于匯集雨水，能夠為泥石流提供較好的水源條件，且溝道內具有常年的流水，水流流量約3～6L/s，匯水及水源條件均較好。

1.1 泥石流匯水區（清水區）條件

清水區位于高程2100～2600m左右，相對高差約500m，平均縱坡降約321.6‰，清水區內山高坡陡，匯水面積2.26km2，區內坡面植被覆蓋率較高，多以喬、灌木為主。區內地形地貌條件能夠起到匯集大氣降水及阻止地表水快速下滲的作用，并能保證地表水迅速匯集至流域溝道中，為泥石流的形成提供良好的水源條件。

1.2 泥石流形成流通區（物源區）條件

形成流通區位于高程1700～2100m左右，平均縱坡降約255.35‰，為明顯的中山峽谷地貌，兩岸岸坡高度一般100～600m，坡度一般30°～50°，最大坡度約80°，區內溝道彎曲、時寬時窄，溝道寬度一般在 3～6m左右，溝谷多呈“V”字型，存在多處跌水，形成流通區內由于坡度較陡，地震后岸坡巖體松動，表層松散物質大量存在，區內存在2處人工棄渣和多處崩滑塌現象，加之豐富的溝道堆積物源，是瓦扎乃曲溝泥石流爆發形成的主要物源來源。

1.3 泥石流溝口堆積區條件

堆積區位于泥石流溝溝口，高程1670～1690m之間，相對高差約20m，縱坡比降122.88%，主要為牛洛村三組居民的耕地及居住地，S103省道瓦扎乃曲溝橋函橫跨溝道，堆積區長度約 230m，面積約0.07km2，堆積體厚度約 5～8m。堆積物巖性為塊碎石土，大小混雜堆積，固結度較高，顆粒成分以白云巖及含礫砂巖為主，粒徑以10～50cm為主、最大可達3.0m。

堆積區溝道相對較寬闊，溝道寬7～22m，侵蝕、切割較嚴重，溝道內目前基本無水流通過，由于該區主要為牛洛村三組居民居住生活用地，故區內受人類工程活動影響較大，一般表現為耕地和建房等，地形地貌受到較大破壞。

2 瓦扎乃曲泥石流形成條件

2.1 水源條件

2.1.1 降雨條件

根據《四川省中小流域暴雨洪水計算手冊》所附暴雨量等值圖，美姑縣巴古鄉地區的1/6h、1h、6h、24h多年最大暴雨量平均值分別為 12.5mm、35mm、60mm、70mm。變異系數分別為 0.35、0.40、0.40、0.35。通過暴雨手冊可以了解到瓦扎乃曲泥石流不同頻率下的雨強。從雨強信息可以看出美姑降雨量變化較大，其中1/6h可達12.55mm。暴雨為泥石流形成的基本影響因素。

2.1.2 流域匯流條件

1）水系和地形因素對匯水的影響。

瓦扎乃曲溝流域主要為深切割中山侵蝕地形，總體上地形陡峻，有利于降雨的匯聚。流域內主支溝溝谷縱坡較大，且表現出上游陡，向下游逐漸縱坡變緩的特點，谷寬則表現為從上游至下游逐漸變大的特點，這種特點控制下，流域水動力條件強大，輸沙能力較強。

此外，流域內具有三條泥石流支溝，各支溝長度、縱坡、溝道形態特征表現出一定的差異性，因而各支溝匯水形成洪峰達到溝口的時間也呈現錯峰到達的特點，這種匯水特征也決定了泥石流運動過程中主支溝峰值流量的疊加特點，當各支溝群發泥石流時，可能出現多峰和陣流的現象。

2）地層巖性分布對匯水的影響

從瓦扎乃曲溝流域地層巖性的分布情況看，上游各支溝溝源部位以基巖為主，覆蓋層較薄，各支溝主要為崩坡積物分布，溝道局部平緩地段以及溝口牛洛村三組區主要為泥石流堆積物分布。

基巖出露區總體上有利于降雨形成地表徑流，而松散堆積層分布區則滯水作用較強，部分降水補充地下水，因而相對不利于形成地表徑流。

3）植被分布對匯水的影響

瓦扎乃曲流域內植被分布特征表現為：上游地區植被覆蓋率較高，僅局部滑塌及坡面侵蝕對植被造成了一定程度的破壞；支溝一中下游地區植被覆蓋率較低，局部受人類農業活動、放牧影響。

植被發育地區滯水作用較強烈，不利于降雨形成地表徑流，而植被覆蓋較差的地區則更加有利于地表徑流的形成。從泥石流松散固體物源的啟動情況和分布特征來看，森林植被發育的地區啟動物源量相對較少，而森林植被相對稀少的地區啟動和參與泥石流活動的物源量則較多，這種特征除與植被發育根固作用阻止起動外，也與植被發育區對降雨的滯水作用具有很大的關系。

2.2 物源條件

泥石流的形成與發生與充足的物源補給密切相關，瓦扎乃曲溝泥石流物源根據形成條件不同主要可分為：崩滑堆積物源（人工棄渣）、溝道堆積物源和坡面侵蝕物源三種。

根據勘查統計，瓦扎乃曲溝流域內崩滑類物源主要分布有14處（包括人工棄渣2處），坡面侵蝕物源分布3處。崩滑總量20.66×104m3，可能參與泥石流活動的動儲量約3.25×104m3；溝道堆積固體物源總量為8.5×104m3，可能參與泥石流活動的動儲量為2.7×104m3；坡面侵蝕固體物源總量為2.86×104m3，可能參與泥石流活動的動儲量為 0.64×104m3。共計有松散固體物源量 32.02×104m3，可能參與泥石流活動的動儲量為6.59×104m3。

3 瓦扎乃曲泥石流特征值參數

3.1 泥石流流體重度

3.1.1 現場配方法

本次勘查共在現場做了 6組現場稱重試驗，在現場由目擊者配制和認定泥石流樣品，最后算出泥石流容重。根據式：

表1 泥石流流體現場稱重表

γc—泥石流容重（t/m3）；P1—桶與泥石流體總重（g）；P2—空桶重量（g）；V—桶內所盛泥石流體體積（cm3）；通過現場試驗各數據如表1。

1）查表法。根據《泥石流災害防治工程勘查規范》（DZ/T0220-2006）附錄 G 泥石流溝易發程度數量化評分標準對瓦扎乃曲溝泥石流進行了評分，評分結果為115分，查表G.2得出泥石流流體重度為1.793t/m3，泥砂修正系數為1.932。

2）經驗公式計算法。

根據蔣忠信教授《震后山地地質災害治理工程設計概要》對泥石流重度以及泥沙修正系數進行估算；對泥石流重度大于 1.5t/m3的泥石流可根據黏粒含量 P05(粒徑＜0.05mm，小數表示)和粗粒含量 P2（粒徑＞2mm，小數表示）按余斌公式估算：

據試驗數據分析，P05=0.053，P2=0.736，代入余斌公式得：

其中：γc—泥石流重度；γw—清水的重度；γH—泥石流中固體物質重度，取 2.4t/m3；經計算，泥石流泥沙修正系數：φc=1.69。

3）綜合取值。由于瓦扎乃曲溝泥石流沒有監測資料，只能通過配方法、查表法和經驗公式計算三種方法來確定，而這三種方法又各有其特點。配方法只能對已發生過，且有人目擊的泥石流進行測定，且測定結果只能代表當時的一次泥石流發生的結果，如本次勘查測定的瓦扎乃曲溝泥石流重度是參照2013年泥石流的情況測定的，由于此次泥石流發生時水量較少，屬黏性泥石流，因而測定結果相對吻合。相對而言，查表法是在現狀調查基礎上帶預測性的重度值結果，但查表法取值相對保守，經驗公式計算法相對前兩種方法有一定的可靠性，因經驗公式是經過成千上萬現場實驗推導而出，具有一定的可靠性，且三種方法結果基本相同，因此本次泥石流重度的綜合取值采用查表法、配方法和經驗公式計算法取平均值作為本次泥石流流體重度的最終結果，即γc=1.88t/m ，φc=1.69

3.2 泥石流流量

3.2.1 雨洪法

1）清水流量計算。用該方法計算是選用《泥石流災害防治工程勘查規范》（DZ/T0220-2006）中雨洪法計算泥石流流量，首先查《四川省中小流域暴雨洪水計算手冊》附圖2-1～2-8，附表6得到美姑縣巴古鄉地區10min、1h、6h、24h相應暴雨計算參數，并按下列公式進行計算：

首先計算頻率為1%、2%、5%、10%的暴雨洪水設計流量Qp，按推理公式：

式中，Sp--某頻率的雨力，mm/h。其計算式為：Sp=Htp·tn-1

Htp--設計頻率最大t小時暴雨量；n--暴雨參數；τ--匯流時間,h。L --主溝長度，km；J --主溝底坡降，‰；F--流域面積，km2；μ--產流參數,mm/h；

計算過程如下：

匯流參數；m=0.221θ0.204=0.296

通過計算得出的瓦扎乃曲溝各斷面部位暴雨洪峰流量值詳見表2。

2）泥石流峰值流量計算。

式中，Qc—泥石流斷面峰值流量（m3/s）；φ—泥沙修正系數，根據計算泥沙修正系數為1.69；Qp—暴雨洪峰流量（計算結果見表 2）；Dc—堵塞系數，按勘查規范表 I.1查表確定堵塞程度為中等，其值為1.8。

據此，采用雨洪法求得泥石流峰值流量如表3。

表2 瓦扎乃曲溝暴雨洪峰流量計算表

表3 瓦扎乃曲溝泥石流峰值流量計算表

3.2.2 形態調查法

對已發生的泥石流流量可采用形態調查法進行計算，根據調查得到的泥石流泥位及據此求得的泥石流流速，再乘上調查求得的過流斷面求取泥石流流量。

形態調查法計算公式為：Qc=WcVc

式中Qc—泥石流斷面峰值流量（m3/s）；Wc—泥石流過流斷面面積（m2）；Vc—泥石流斷面平均流速（m/s）。

由于泥石流重度為1.88t/m3，為黏性泥石流，因此其流速采用鐵二院黏性泥石流計算公式：

式中：Vc—泥石流流速（m/s）；Hc—平均泥深（m）；Ic—泥位縱坡率，以溝道縱坡率代替。

nc_——黏性泥石流河床糙率，根據《泥石流災害防治工程勘查規范》（DZ/T0220-2006）表I.4。

并根據蔣忠信教授《震后山地地質災害治理工程設計概要》中采用彎道超高差對流速反算進行復核；V=｛Rg{[Δh/B-tanφ-c/(Hγc o s2θ)]｝據此求得各斷面位置泥石流流速特征值詳見表4，泥石流流量計算結果詳見下表5。

表4 瓦扎乃曲溝泥石流流速計算表

表5 瓦扎乃曲溝泥石流流量形態調查法計算表

3.2.3 綜合取值

通過表3和5的計算結果，采用形態調查法求得的泥石流峰值流量計算結果P=5%的雨洪法求得的結果相差不大，但采用形態調查法計算結果沒有暴雨頻率的概念[1]，僅能代表當次泥石流的特征值，并且由于年代久遠當次泥石流痕跡調查誤差很可能較大，而雨洪法則根據現有溝域面積、溝域植被發育分布情況和徑流系數進行計算，具有預測性質，因此，泥石流峰值流量采用表3雨洪法計算求得的結果。

3.3 一次泥石流過流總量

一次泥石流過流總量按照《泥石流災害防治工程勘查規范》（DT/T0220-2006）附錄I提供的計算公式進行計算：

式中，Q—一次泥石流過程總量（m3）；T—泥石流歷時（s）；Qc—泥石流最大流量（m3/s）；

據上節計算結果，2006年瓦扎乃曲溝泥石流為近期發生的規模相對較大的一次泥石流，對比斷面法計算結果，其流量相當于P=5%的雨洪法計算的峰值流量，綜合確定泥石流峰值流量為86.38m3/s；據訪問，泥石流集中歷時約15分鐘，即T900s；按上式計算的本次泥石流沖出量為Q=2.05×104m3。

3.4 一次泥石流固體沖出物

一次泥石流固體沖出物按照《泥石流災害防治工程勘查規范》（DT/T0220-2006）附錄I提供的計算公式進行計算：

式中QH—一次泥石流沖出固體物質總量（m3）；Q—一次泥石流過程總量（m3）γc—泥石流重度（t/m3）；γw—水的重度（t/m3）；γH—泥石流固體物質的重度（t/m3）。

據前面的計算結果和分析測試結果，瓦扎乃曲溝 2006年泥石流暴發時泥石流沖出量為 Q=2.05×104m3，泥石流重度γC=1.88t/m3，水的重度γW=1t/m3，泥石流固體物質重度γH=2.4 t/m3，據此計算出一次性固體物質沖出總量為QH=1.29×104m3。

4 瓦扎乃曲泥石流防治措施

瓦扎乃曲溝泥石流治理思路主要以攔排結合，通過攔擋壩穩固溝內堆積物源，并有一定庫容量攔截泥石流固體物源，并調節下游洪峰流量，同時在溝口段修建排導槽，確保洪峰安全通過，不至對溝口民房再次造成危害。

治理方案：攔砂壩+排導槽+單邊防護堤。

因瓦扎乃曲溝流域內物源儲量豐富，其主要來源為溝道兩側松散崩、滑塌物源及溝道堆積物源，針對這一特征可在形成流通區下游地形寬緩地段設置 1～2級，攔蓄泥石流固體沖出物，通過回淤壓腳降低溝道坡降及固體沖出物；以及針對堆積區溝道沖刷、切割嚴重，可修建排導槽（鋪底），為保護右側居民在堆積區右岸修建單邊防護堤與排導槽相連，將攔擋后剩余泥石流流體直接引向美姑河排泄。

5 結論與建議

1）瓦扎乃曲溝泥石流為黏性泥石流，屬暴雨型泥石流，暴發頻率為中頻，暴發規模多為小型。

2）根據瓦扎乃曲溝泥石流形成條件、運動和堆積特征結合泥石流危害對象的分布，本次瓦扎乃曲溝泥石流防治工程采用攔排結合的治理思路，采用“排導槽+單邊防護堤+攔砂壩”治理方式。

3）建議地方政府制定專門的規劃，加強溝域植被保護和和水土保持，通過長期的生物防治方案，減少泥石流活動物源，加強溝域地質生態環境保護，避免因不合理人類工程活動加劇不良地質現象，增加泥石流物源。

4）建議治理工程完成前進一步加強地質災害群測群防工作，加強泥石流活動的監測預警，避免因泥石流災害造成人員傷亡和大的經濟財產損失。

[1]劉德兵.茂縣瓦窯溝泥石流動力學特征分析與防治對策研究.川煤地勘,21（2013）1-0096-06

[2]黃潤秋,李為樂.“5·12”汶川大地震觸發地質災害的發育分布規律研究[J].巖石力學與工程學報,2008,27(12)： 2585-2592.

[3]黃潤秋.“5·12”汶川大地震地質災害的基本特征及其對災害后重建影響的建議[J].中國地質教育,2008, (2)：21-24.

[4]蔣忠信.震后山地地質災害治理工程設計綱要.西南交大出版社,2013.7.

[5]四川省一三五勘察設計有限公司.涼山州美姑縣巴古鄉瓦扎乃曲泥石流勘查報告.2014.05