楊尚遠，趙 懷

(1.中廣核工程有限公司，廣東深圳518124；2.甘肅地質工程勘察院，甘肅張掖734000)

寺溝泥石流沖擊力對影響因素的敏感性分析

楊尚遠*1，趙 懷2

(1.中廣核工程有限公司，廣東深圳518124；2.甘肅地質工程勘察院，甘肅張掖734000)

泥石流整體沖擊力和石塊沖擊力對不同因素的敏感性存在明顯差異，以寺溝泥石流為例，在取得水文與環(huán)境地質條件資料、設計頻率下的各因素基準值和各因素變動范圍基礎上，最終分析比較寺溝泥石流整體沖擊力和石塊沖擊力對10個因素的敏感系數。據此做出相應的治理措施建議，確定沖擊力計算中需提高精度水平的參數。

寺溝；泥石流沖擊力；影響因素；敏感系數

1 概述

中國是泥石流多發(fā)的國家，泥石流暴發(fā)造成的損失是巨大的。在泥石流危險評估和防治工程設計中泥石流沖擊力是最重要的力學參數[1]，在泥石流綜合治理中如何有效地減小泥石流沖擊力變得尤為重要，比較泥石流沖擊力對各影響因素的敏感性，對在泥石流綜合治理中提出科學合理的治理措施、在防治工程設計中提高泥石流沖擊力的計算精度具有重要意義。

唐川等根據巖石的不同種類及軟硬程度對巖土體進行分類，對巖性進行泥石流敏感性計算[2-3]等等。本文抽取泥石流沖擊力的10個影響因素：流域面積、溝槽縱坡度、年最大24h降雨量均值、匯流參數、降雨歷時等于24h的徑流系數、泥石流密度、堵塞系數、建筑物形狀系數、受力面與沖擊力夾角、石塊質量，引入敏感系數概念，定義敏感系數計算公式，比較各因素對泥石流整體沖擊力、石塊沖擊力的敏感系數，以寺溝泥石流為例進行敏感性比較、做出寺溝泥石流綜合治理措施建議和確定沖擊力計算中需提高精度水平的參數。

2 敏感性分析方法

泥石流沖擊力分為泥石流整體沖擊力和石塊沖擊力，需分別計算敏感性，如圖1所示。步聚如下：

圖1 泥石流沖擊力敏感性分析技術路線

（1）泥石流整體沖擊力和石塊沖擊力受各種內在因素和外部因素變動的影響，諸如氣溫、降雨、巖土類型與性質、植被、流域、泥石流溝形態(tài)、人為活動等，抽取泥石流整體沖擊力和石塊沖擊力的10個影響因素：流域面積、溝槽縱坡度、年最大24h降雨量均值、匯流參數、降雨歷時等于24h的徑流系數、泥石流密度、堵塞系數、建筑物形狀系數、受力面與沖擊力夾角、石塊質量。

（2）選擇10個因素的基準參數，計算基準條件下泥石流整體沖擊力和石塊沖擊力，在此基礎上變動其中的某一參量，其它參數固定不變，計算出該參量在其變化范圍內變動時沖擊力隨之變化的結果。

（3）根據上述計算結果，引入敏感系數S：

式中，對于泥石流整體沖擊力，η1=， η2=；對于石塊沖擊力，η1=， η2=，ΔX為某因素的變化量，Xmax-Xmin為某因素最大變化量，Δδ、ΔG為δ、G對應ΔX的變化量，δO、GO為δ、G的基準值。

（4）根據計算結果繪制各因素 X的Sδ-X、SG-X曲線，比較各因素基準值處的敏感系數大小，分析總結，做出泥石流綜合治理措施的建議和確定沖擊力計算中需精度水平的參數。

其中，傳統(tǒng)泥石流流速是按

計算，在實際應用中式（2）的泥石流流速計算精確度不高[4]，所以本文應用推理公式計算清水流理[5]（計算方法在此不贅述），據此計算泥石流流量，由于應用推理公式計算泥石流流量計算精度較式（2）計算精度高，試采用公式（3）來計算泥石流流速：

式中：QC——泥石流流量；

A——過泥斷面面積。

3 沖擊力敏感性分析

3.1 水文與環(huán)境地質條件

寺溝泥石流位于甘肅省隴南市武都縣玉皇鄉(xiāng)。東靠黃坪鄉(xiāng)，西鄰三河鄉(xiāng)，北連郭河鄉(xiāng)，南與琵琶鄉(xiāng)接壤，寺溝溝口位于福津河右側石家村。溝口地理坐標東經105°12′12″，北緯33°18′45″。其氣象條件見圖2。流域平面形態(tài)近似“殘月”形，流域內溝谷分布形態(tài)呈“桃葉”形，流域平均縱向長度約6.47km，平均寬度2.08km，流域面積13.5km2。流域最高點位于山尖上，海拔高度2500m，最低點位于寺溝與福津河交匯處，海拔高度1180m，相對高差1320m。全流域溝谷為深切割“V”型谷，具有岸坡陡峻、切割深度較大等特點。主溝縱長3.05km，主溝形成點即兩支溝交匯處高程為1800m，與溝口相對高差620m，主溝平均縱坡降203.5‰，其中上游溝段縱坡陡峻，平均縱坡360.9‰，而下游至岷江入口段縱坡略緩，縱坡150‰～180‰，且呈現陡緩相間的空間變化特征。溝腦兩條支溝主溝道長分別為3.58km、1.64km，左側支溝溝谷比降125‰，右側支溝溝谷比降183‰。溝槽總體上較為狹窄，形態(tài)呈“U”形，溝槽寬度20～80m，但局部地段也較為狹窄，形成寬窄相間的變化特點，局部形成卡口。

圖2 寺溝氣象要素圖

按照不同分類方法，寺溝為暴雨型、坡面侵蝕型、溝谷型、高頻、稀性、衰退期泥石流。

3.2 各因素參數選取與敏感性計算

寺溝泥石流治理工程按2%的設計頻率設計，10個影響因素的基準值見表1，各因素最大值、最小值確定如下：

表1 泥石流整體沖擊力、石塊沖擊力在各因素基準值處對各因素的敏感系數

（1）流域面積。根據遙感資料，我們發(fā)現上游左岸部分流域有可能通過工程措施使其減小，面積為Fmax=13.5km2，Fmin=12.6km2。

（2）主溝縱坡度。對寺溝調查發(fā)現寺溝各段溝槽縱坡度不同，上游溝段縱坡陡峻，平均縱坡度20°，下游至岷江入口段縱坡低緩，平均縱坡度8.5°,則選取βmax=20°，βmin=8.5°。

（3）年最大24h降雨量均值。對寺溝多年的降雨量統(tǒng)計發(fā)現在不同年代計算出的年最大24h降雨量均值在86～93mm內變動，確定-H24max=93mm，-H24min= 86mm。

（4）匯流參數。因寺溝流域系數為2.71，查《匯流參數值表》可知，此流域系數下的匯流參數的變化范圍為0.2～1.3，選取mmax=1.3，mmin=0.2。

（5）降雨歷時等于24h的徑流系數。查《降雨歷時等于24h的徑流系數表》可知，據寺溝最大24h降雨量、流域內土體類型可能變化情況，確定αmax=0.75，αmin= 0.6。

（6）泥石流密度。泥石流密度的下限為清水密度，上限為 1900kg/m3，則 γCmax=1900kg/m3，γCmin= 1000kg/m3。

（7）堵塞系數。查《堵塞系數表》可知堵塞系數在1～13內取值，則DCmax=13，DCmin=1。

（8）建筑物形狀系數。圓形建筑物為1.0，矩形建筑物為1.33，方形建筑物為1.47，則λmax=1.47，λmin= 1.0。

（9）受力面與沖擊力夾角。取θmax=90°，θmin=0°。

（10）石塊質量。此設計頻率下泥石流能夠搬運的最大石塊質量為Wmax=3073kg，Wmin=0kg。

按照敏感性計算步驟得出各因素不同值對應的敏感系數，繪制敏感系數—影響因素回歸曲線，并確定曲線函數表達式，見圖3～圖12。

圖3 泥石流整體沖擊力、石塊沖擊力敏感系數—流域面積關系曲線

圖4 泥石流整體沖擊力、石塊沖擊力敏感系數—主溝縱坡度關系曲線

圖5 泥石流整體沖擊力、石塊沖擊力敏感系數—年最大24h降雨量均值關系曲線

圖6 泥石流整體沖擊力、石塊沖擊力敏感系數—匯流參數關系曲線

圖7 泥石流整體沖擊力、石塊沖擊力敏感系數—降雨歷時為24h徑流系數關系曲線

圖8 泥石流整體沖擊力、石塊沖擊力敏感系數—泥石流密度關系曲線

圖9 泥石流整體沖擊力、石塊沖擊力敏感系數—堵塞系數關系曲線

圖10 泥石流整體沖擊力、石塊沖擊力敏感系數—建筑物形狀系數關系曲線

圖11 泥石流整體沖擊力、石塊沖擊力敏感系數—受力面與沖擊力夾角關系曲線

圖12 泥石流整體沖擊力、石塊沖擊力敏感系數—石塊質量關系曲線

圖3～圖12說明，在設計頻率下，泥石流整體沖擊力和石塊整體沖擊力對10個因素的敏感性大小有較大差異。

4 結論與建議

寺溝泥石流整體沖擊力對堵塞系數、匯流參數、泥石流密度3個因素的敏感性大，對主溝縱坡度、建筑物形狀系數、降雨歷時等于24h的徑流系數、年最大24h降雨量均值、流域面積5個因素的敏感性中等，對受力面與沖擊力夾角的敏感性小，對石塊質量不具敏感性。

寺溝泥石流石塊沖擊力對堵塞系數、匯流參數、泥石流密度、石塊質量4個因素的敏感性大，對主溝縱坡度、降雨歷時等于24h的徑流系數、年最大24h降雨量均值、流域面積4個因素的敏感性中等，對受力面與沖擊力夾角的敏感性小，對建筑物形狀系數不具敏感性。

綜述，對于寺溝泥石流整體沖擊力和石塊沖擊力，堵塞系數、匯流參數、泥石流密度、石塊質量為敏感性大的因素，提出以下建議：

（1）由于4個因素的變化對泥石流整體沖擊力和石塊沖擊力影響較大，在寺溝泥石流沖擊力計算中須提高4者的精度水平；

（2）在寺溝泥石流綜合治理過程中要控制寺溝泥石流的整體沖擊力和石塊沖擊力，關鍵在減小寺溝的堵塞系數、匯流參數、泥石流密度、石塊質量。開展以下措施：

①堵塞系數由流通區(qū)溝槽彎曲度、寬窄不均勻性、泥石流粘性稠度決定，治理中需排除流通區(qū)中的大石塊，平直順暢流通區(qū)，減少生活垃圾向泥石流溝的排放，減小徑流對物源區(qū)粘性土的沖蝕，在物源區(qū)積極植樹造林；

②匯流參數由雨洪特征、清水匯集溝槽特性、土壤植被、流域特征系數決定，泥石流治理時需增加清水匯集區(qū)的植被覆蓋率特別是灌土林的覆蓋率、增加形成區(qū)溝槽的大卵石、大礫石的數量、杜絕開墾梯田；

③泥石流密度由降雨量、物源區(qū)巖土類型與性質、植被覆蓋率所影響，泥石流治理中應增加植被覆蓋率、治理物源區(qū)的不良地質災害（如崩塌、滑坡、危巖體）、增加出露巖土體的強度減小其風化速度；

④石塊質量受石塊重度、石塊體積、石塊強度的影響，在物源區(qū)應防止大石塊進入流通區(qū)，排除流通溝槽的大石塊，同時對物源區(qū)的不良地質災害進行治理。

[1]費祥俊，舒安平.泥石流運動機理與災害防治[M].北京：清華大學出版社，2004.

[2]唐川.金沙江流域(云南境內)山地災害危險性評價[J].山地學報,2004,22(4):451-460.

[3]唐川.云南怒江流域泥石流敏感性空間分析[J].地理研究, 2005,24(2):178-185.

[4]中國科學院成都山地災害與環(huán)境研究所.泥石流研究與防治[M].四川：四川科學技術出版社，1989.

[5]葉鎮(zhèn)國.土木工程水文學[M].北京：人民交通出版社，2000.

P642.23

A

1004-5716(2015)03-0019-05

2013-06-05

2014-03-31

楊尚遠（1975-），男（漢族），廣西靈山人，工程師，現從事地質環(huán)境災害預測與防治工作。