何夏裕 劉蕾 王勇智

收稿日期：2023-12-10

作者簡介：何夏裕（1999—），男，浙江紹興人，主要從事地質(zhì)工程研究。#通信作者：劉蕾 ，E-mail：769245865@qq.com。

摘 要：泥石流運動參數(shù)特征及危險性評價是對米林比定村段內(nèi)活動人員和公路等設施安全的重要依據(jù)。為了研究比定溝對區(qū)域內(nèi)的影響，采用實地勘察、室內(nèi)試驗、遙感分析等相結(jié)合的方法，分析了泥石流溝運動參數(shù)和危害特點，對溝道進行泥石流危險性評價，為米林的災害防治提供科學依據(jù)。

關鍵詞：泥石流；運動參數(shù)；危險性

中圖分類號：P642.23 文獻標志碼：B文章編號：2095–3305（2024）03–0-03

米林位于西藏東南地區(qū)，是我國泥石流災害最為嚴重的地區(qū)之一，這與其特殊的地理、氣候和地質(zhì)條件密不可分[1]。位于雅魯藏布江中段、念青唐古拉山與喜馬拉雅山交匯處，由于地震活動頻繁，導致該區(qū)域地質(zhì)巖層較為破碎。同時，因為該地區(qū)降雨量較為豐富，而泥石流作為一種具有突發(fā)性、高度破壞性、淤積性等特性的自然災害，給國家?guī)砭薮蟮膿p失，甚至威脅人民的生命安全[2]。目前，許多專家先后提出過一些泥石流危險評價的方法，如模糊數(shù)學法、投影追蹤法

等[3-4]?；谏鲜龇椒ǖ难芯?，以米林市比定溝作為研究對象，借用熵理論并選取一定的評價因子作為評價指標進行危險性評價，為比定溝后期治理提供依據(jù)[5]。

1 研究區(qū)概況

1.1 位置與地形地貌

比定溝流域位于西藏自治區(qū)米林市比定村向西1.8 km處，溝口坐標為94°5′41.140"E，29°11′38.806"N

（圖1）。該流域整體由西南—北東向展布，山體被深而長的縱谷所切割，地勢陡峭，地形崎嶇，屬于復合山地地貌和河谷地貌，溝谷形態(tài)呈“I”字形。該流域面積0.635 km2，主溝長度為1.978 km，流域內(nèi)高差為1 155 m，

地勢南高北低，縱坡降為584‰，平均坡度為38°，低坡度主要集中在溝槽和分水嶺地帶，而高坡度主要集中分布在山坡兩岸不同的高度帶上。

1.2 地質(zhì)及氣候水文條件

在構造條件上該區(qū)域位于南迦巴瓦楔形構造的西側(cè)，并從附近露出的地層可以得到，該區(qū)域地層巖性為新元古界—中元古界念青唐古拉巖群八拉巖組（Pt2-3b）的片巖，及第四系全新統(tǒng)風積層（Q4heol）。片巖主要分布于流域的主支溝溝道內(nèi)；風積層主要分布于流域溝口山體斜坡地帶，多為塊石與顆粒物。

在氣候上屬于高原溫帶半濕潤季風氣候區(qū)，氣候特點為降水集中，雨熱同季，蒸發(fā)量大。因溝道與米林市距離較近，因此引用米林氣象站1987—2020年的降雨數(shù)據(jù)統(tǒng)計，降雨多集中在6—9月，年平均氣溫8.6 ℃，平均年降水量為702 mm，年降雨量最高為891.9 mm。

該流域內(nèi)蒸發(fā)量大、水汽充沛，夏季雨水集中，在干濕交替的作用下，容易使得土體失去穩(wěn)定性，從而變得松散，加大泥石流的物源量。

2 泥石流運動參數(shù)的計算

2.1 泥石流容重

泥石流容重的大小直接反映了泥石流對周圍環(huán)境的影響程度，同時也反映了泥石流混合物中含有固體物質(zhì)的數(shù)量。鑒于比定溝沒有相關的歷史紀錄，且時間跨度較長，根據(jù)黏粒含量計算泥石流容重和泥石流溝的易發(fā)程度進行數(shù)值評估，以便更精確地測定泥石流的容重[6]。通過現(xiàn)場采集到的樣品，并進行篩分，并利用馬爾文激光粒度儀器測量樣品中黏粒含量，結(jié)果見圖2。

rc=-1.32×103x7-5.13×102x6+8.91×102x5-55x4+34.6x3-67x2+12.5x+1.55（1）

式（1）中，rc為泥石流容重g/cm3；x為泥石流土樣中黏粒占的含量。

比定溝根據(jù)計算得出：黏粒含量為83.99%，對應容重為1.651 g/cm3；綜合得分為98分，對應容重為1.676 g/cm3。因土樣在現(xiàn)場長期受到雨水的影響，容易導致樣品中黏粒含量的減少，與實際容重存在偏差。故選取較大的值作為泥石流容重的最終值，泥石流容重為1.676 g/cm3。

2.2 暴雨洪峰流量

根據(jù)條件分析，該溝道為暴雨型泥石流溝。由于比定溝泥石流發(fā)生時間過于久遠，已經(jīng)無法準確確定。因選用“小流域暴雨洪峰流量計算方法”算出暴雨洪峰流量，通過米林氣象站1987—2020年的數(shù)據(jù)分析所得年最大日降雨量均值（H24）和變差系數(shù)（H6），再利用P-Ⅲ曲線理論的相關方法可確定工作區(qū)1 h和6 h最大降雨量均值（H1和H6）和變差系數(shù)（CV1和CV6）。計算結(jié)果見表1。

QP=0.278ψF（2）

式（2）中，QP為暴雨洪水設計流量（m3/s）；ψ為洪峰徑流系數(shù)；s為暴雨雨力；τ為匯流時間（h）；n為暴雨指數(shù)；F為流域面積。

2.3 泥石流峰值流量

該溝采用雨拱法，根據(jù)不同頻率下的暴雨洪峰流量，選用合適的溝道堵塞系數(shù)，計算出泥石流峰值流量?？紤]到溝槽彎道發(fā)育但曲率不大，溝道寬度有一定變化，局部有陡坎、卡口分布，還聚集著較豐富的松散堆積物源，并對溝道造成一定程度的堵塞，因此確定溝床堵塞情況為中等，取值2.0～2.5。

QC=（1+φ）QP·DC（3）

φ=（γC -γW）/（γH -γC）（4）

利用式（3）和式（4），計算比定溝100年一遇的峰值流量，φ為泥沙修正系數(shù)取0.51；γH為固體物質(zhì)比重取26.5 kN/m3；γW為清水重度9.8 kN/m3；DC為流堵塞系數(shù)取2.41，經(jīng)計算得出比定溝百年一遇的泥石流峰值流量為16.93 m3/s。

2.4 一次泥石流總量和沖出固體物質(zhì)

基于固體物質(zhì)總量和泥石流峰值流量統(tǒng)計關系的經(jīng)驗模型、五邊形法和依據(jù)流域面積增加系數(shù)的修正五邊形法估算，公式如下：

WC=0.264TQC（5）

WS=（γC -γW）WC/（γH -γW）（6）

式（5）中，WC為一次泥石流總量；T為泥石流歷時取1 200 s。經(jīng)計算，比定溝泥石流一次總量為0.54萬m3，

固體物質(zhì)總量為0.22萬m3。

2.5 泥石流流速

由以上數(shù)據(jù)可知，該流域為稀性泥石流溝，道內(nèi)較為粗糙，現(xiàn)場可見多個由于大顆粒聚集而形成的陡坎，通用計算公式如下。

VC=R2/3IC1/2（7）

式（7）中，VC為斷面平均流速（m/s）；R為水力半徑；IC為流域比降為0.095；泥沙深度取3 m；φC為泥沙修正系數(shù)0.65；n為粗糙系數(shù)7.8，根據(jù)式（7）可得泥石流斷面流速為3.03 m/s。

2.6 泥石流最大沖起高度與沖壓力

鑒于比定溝沒有相關的歷史記錄，且時間跨度較長無法核實，故按照《泥石流災害防治工程勘查規(guī)范》提供的計算公式計算泥石流動力學特征，得到最大沖起高度和整體沖壓力。

?HC=（8）

δ=λV 2C sinα（9）

式（8）和式（9）中，?HC為最大沖起高度（m）；δ為整體沖壓力（Pa）；α為建筑物受力面與泥石流沖壓力方向的夾角，取90°；λ為建筑物形狀系數(shù)，矩形取1.33。 計算結(jié)果表明，比定溝的最大沖起高度0.47 m、沖擊力為20.46 kPa。

3 泥石流危險性評價

泥石流的危險評價是對其發(fā)展趨勢的預測，而在泥石流的形成過程中，其影響因素較多，這些因素在不同的時空條件下，它的影響程度又大不相同。將危險度作為一個定量，系統(tǒng)地借助了熵理論用以表達各因子的效用價值，以此準確評價泥石流的危險性。

在劉希林等[7]提出的泥石流危險評價方法的基礎

上，選用容重、流域面積、主溝長度、單位面積物源量、泥石流整體沖壓力、流域相對高差、流域切割密度、不穩(wěn)定溝床比例8個因素作為危險度評價因子，其中，泥石流規(guī)模按照100年一遇取值（表2）。

采用熵值法對泥石流溝的危險性評價，因熵值法的系數(shù)K需要2個以上對象時才能計算，故將比定溝相鄰的雙溪溝和牌林溝參與運算[8-10]。

3.1 評價指數(shù)標準

通過對原始數(shù)據(jù)的處理，達到消除量綱的目的。本文以大者為優(yōu)的指標。

γij=（10）

式（10）中，γij為量綱化處理后的數(shù)據(jù)；min（ j）為某項因子中最小值；max（ j）為某項因子中最大值；xij為原始數(shù)據(jù)。

3.2 熵值計算公式

Ej=KPijlnPij（11）

（12）

式（11）和式（12）中，Ej為第j項因子熵值；Pij為第i項指標在泥石流災害中出現(xiàn)的頻率；m為樣本的個數(shù)；n表示因子的個數(shù)[11-13]。

3.3 計算評價指標權重值

（13）

3.4 計算危險性值

Fj=Wjrij（14）

確定各因子權重值：U={0.15，0.08，0.19，0.08，0.15，0.13，0.13，0.10}，最終計算得到泥石流危險值分別為0.40，按照泥石流危險性分級所示：0.4～0.6屬于中度危險，可知比定溝屬于中度危險。

4 結(jié)論

基于上述泥石流運動參數(shù)計算與熵值法對比定溝8個影響因子進行分析，得出以下結(jié)論：

（1）該流域的地形地貌、地質(zhì)情況及水源條件皆有利于泥石流的形成。根據(jù)比定溝泥石流發(fā)育特征與參數(shù)計算得知，在規(guī)模上屬于小型泥石流。

（2）基于熵值法對該流域的危險度評價結(jié)果，其處于中度危險之中，會對區(qū)域內(nèi)活動人員和公路等設施安全產(chǎn)生較大影響。

參考文獻

[1] 鐘祥浩，劉淑珍，王小丹，等.西藏高原生態(tài)安全研究[J].山地學報，2010，28（1）：1-10.

[2] 巫建暉，張正波，田金昌，等.川藏公路山地災害特征及對西藏可持續(xù)發(fā)展的影響[J].水土保持研究，2008（4）：142-144.

[3] 陳偉，任光明，左三勝.泥石流危險度的模糊綜合評判[J].水土保持研究，2006（2）：138-139，171.

[4] 王念秦，姚勇.基于模糊數(shù)學和權的最小平方法的泥石流易發(fā)性評價方法[J].災害學，2008（2）：5-9.

[5] 王英杰，王磊，榮起國，等.基于熵值理論的泥石流評價因子選取[J].交通運輸工程學報，2014，14（2）：28-33.

[6] 陳寧生，楊成林，李歡.基于漿體的泥石流容重計算[J].成都理工大學學報（自然科學版），2010，37（2）：168-173.

[7] 劉希林，莫多聞.泥石流風險及溝谷泥石流風險度評價[J].工程地質(zhì)學報，2002（3）：266-273.

[8] 杜國梁，高金川，李瑞冬.基于模糊綜合評價法和熵值法的泥石流危險度評價[J].安全與環(huán)境工程，2013，20（5）：15-17，32.

[9] 孫健，劉海，劉欽，等.小容河流域溝谷泥石流地質(zhì)災害形成條件分析及風險評價[J].安徽地質(zhì)，2020，30（3）：208-213.

[10] 陳果，賈蒼琴.溝谷泥石流災害風險評價[J].內(nèi)蒙古大學學報（自然科學版），2009，40（2）：233-238.

[11] 孫麗蘋，王瑞龍.河南淇縣某溝谷泥石流災害現(xiàn)狀及防治措施分析[J].地下水，2017，39（6）：123-124.

[12] 韓俊，王保云.基于原型網(wǎng)絡的云南怒江州泥石流災害易發(fā)性評價與區(qū)劃[J].中國地質(zhì)災害與防治學報，2023，34（5）：117-129.

[13] 羅健，郭寧，楊崢.安寧河流域西昌北段東岸溝谷泥石流危險度評價[J].中國地質(zhì)災害與防治學報，2018，29（4）：40-49.