冷漬溝泥石流成因分析及危險性評價

［關鍵詞］泥石流；成災因素；數值模擬；危險性分區

泥石流成因分析和危險性評價是泥石流防治工作中的重要組成部分[1-3]。黃潤秋等通過分析汶川地震地質災害的發育特征及演化分布規律發現，震后不良地質災害的發育為暴雨型泥石流的暴發提供了豐富的固體松散物質[4]；華曙光研究基底侵蝕作用下泥石流動力過程，構建了數學模型和流變本構的數值[5]；楊濤等利用數值模擬，研究泥石流溝道的危險性分區，通過危險性分區，為防災減災提供參考[6]。

文章在前人研究的基礎上，通過現場調查結合室內分析，研究冷漬溝泥石流的致災因素、力學參數、溝道危險性分區等特征，為泥石流工程治理提供依據。

1. 泥石流基本特征

冷漬溝為常年流水溪溝，溝域呈口袋狀（溝口地理坐標：102°11'17″E，28°36'33″N），流域匯水面積0.69 km²，主溝長2.0 km，相對高差592 m，平均縱坡比296.0 ‰，溝道總體上緩下陡，中下游溝道狹窄，局部陡坎發育，溝道斷面呈“V”型。溝域內發育三疊系上統白果灣組二段（T₃bg²）的砂巖、泥巖，震旦系蘇雄組一段（Z₁s¹）灰巖，以及第三系昔格達組（NQx）等地層。

2. 泥石流致災因素分析

2.1 地形地貌

冷漬溝位于安寧河左岸，流域地處深切割高中山、高山寬脊溝谷地貌區。兩岸斜坡坡度較大，中、下部30 ～ 50°，局部有大于70°的陡崖。溝域內松散固體物源豐富，主要發育照壁山滑坡，滑坡區坡面總體較緩，結構松散，在降雨驅動下易發生局部滑塌堵潰溝道。

2.2 物源

溝域內松散固體物源主要集中于上游形成區段的照壁山滑坡及流通區（圖1）。通過調查，溝域共發現4處崩滑物源、4處溝道物源、1處坡面物源。物源總儲量與動儲量分別采用弓形均高法、溝谷下切侵蝕堵潰型模式進行計算[7]，計算物源總儲量為357.77×104 m3，物源啟動方式以岸坡側蝕滑塌、揭底沖刷為主，能參與泥石流的動儲量約75.85×104 m3。

2.3 水源

溝域內降水充沛，年均降雨量1222 mm，日最大降雨量218.6 mm，小時最大降雨量150 mm。據《四川省中小流域暴雨洪水計算手冊》，冷漬溝流域的1/6 h、1 h、6 h、24 h 降雨量平均值分別為14 mm、40 mm、70 mm、90 mm。參照泥石流激發雨量研究結果[8]，冷漬溝流域降雨量超過激發降雨量，具備良好的水源條件。

3. Massflow的泥石流危險性評價

3.1 軟件介紹

Massflow軟件目前在泥石流研究方面應用較多，以速度快、規模大、可拓展性強等特點被廣泛使用。

3.2 模型原理

泥石流流體運動過程要遵循質量和動量守恒定律。Massflow軟件將復雜的地質模型簡化為流動層、穩定層兩個界面，通過深度積分建立連續介質力學方程，計算公式見公式（1）～（4），流體的上下界面分別用Z_s和Z_b表示，流體厚度h=Z_s-Z_b，用u、v、w 代表流體速度在x、y、z軸上的分量。

式中：ρ 為流體的密度，kg/m³；t代表時間，s；τ_ij 代表不同方向上的應力分量；g_x、g_y、g_z 分別代表重力加速度g 在x、y、z 軸上的分量，表達成矩陣形式如下：

3.3 基礎數據及參數取值

（1）模型建立將冷漬溝的地形數據DEM轉換成Massflow軟件識別的ASCII碼，建立計算模型。

（2）容重本次泥石流容重的確定采用室內顆分試驗。通過現場調查，冷漬溝泥石流顆粒分選性較差，與黏性泥石流特征相符，根據余斌容重計算公式[9]，冷漬溝泥石流容重為2.11 g/cm³（圖2）。

（3）流量

假定泥石流暴發頻率與降雨頻率相互吻合，根據雨量信息結合溝道調查分析，冷漬溝泥石流是在小時降雨量150 mm雨量下激發（20年一遇P=5%）。根據水利工程設計和鐵一院清水流量計算公式，泥沙修正系數φ=0.31，堵塞系數Dc取值1.2，計算參數如表1：

（4）流量控制線

流量過程線用概化五邊形控制，即以單次泥石流時間的1/3段作為分界點，然后以峰值流量的1/4和1/3作為兩個分界點下的流量，繪制出泥石流暴發的流量過程線；（5）摩擦模型：選用Voellmy模型，主要考慮了紊動系數，可以限制顆粒流的運動速度。

4. 模型結果及分析

通過對冷漬溝泥石流的多次反演，模擬得出冷漬溝泥石流的泥深及流速，根據常鳴提出的泥石流強度影響及危險性分區標準[10]，將模擬得出的泥深及泥深與速度乘積關系進行強度分級，然后將泥石流強度與重現周期相結合，得到冷漬溝流域危險性分區概況（圖2）。

分析可得：冷漬溝泥石流泥深最大為7.50 m，主要堆積于溝口與主河安寧河，地形寬緩變化的地方，泥石流堆積物較厚，因為地形由陡變緩，泥石流流速變小，利于流體物質堆積；冷漬溝泥石流最大流速8.55 m/s，主要分布于冷漬溝溝道中游及中下游，地形越陡，泥石流流速越大，高陡的地形條件為泥石流流動創造了勢能空間。結合實地調查發現，流速的大小還與植被的稀疏及溝道寬緩條件有關；冷漬溝泥石流溝道高危險區占77%，中危險區占13%，低危險區占10%，高危險區主要分布于冷漬溝中下游及溝道與主河安寧河交匯處，且危險性區域較大。

5. 結論

（1）通過對冷漬溝溝道的調查分析，夏季充沛的降雨、特殊的地形地貌及溝道上游滑坡堆積體提供的豐富物源是泥石流主要的致災因素。

（2）經現場試驗及室內分析，冷漬溝泥石流容重為2.11 g/cm³，以黏性為主；通過Massflow數值模擬，再現冷漬溝泥石流運動狀態，模擬得出冷漬溝泥石流最大泥深7.50 m，最大流速8.55 m/s，其中，溝道高危險區占77%，主要分布于冷漬溝中下游及溝道與主河安寧河交匯處，且危險性區域較大。

（3）根據分析泥石流物源分布與數量計算，結合流域內降雨數據、地形地貌及溝道堆積物特征，分析得出：冷漬溝再次暴發更大規模泥石流的可能性極大。