盧 洪 沛， 胡 卸 文

(1.甘孜藏族自治州自然資源和規(guī)劃局，四川 康定 626000；2.西南交通大學(xué)地球科學(xué)與環(huán)境工程學(xué)院，四川 成都 610031)

0 引 言

隨著人類在林區(qū)活動不斷加劇，森林火災(zāi)逐漸演變?yōu)橐环N常見的林區(qū)災(zāi)害，據(jù)不完全統(tǒng)計，近年來中國平均每年發(fā)生森林火災(zāi)約1萬多次。例如2020年3月30日四川省西昌市發(fā)生森林火災(zāi)，火災(zāi)燒毀林地約1 000公頃，2016年3月15日，四川省九龍縣突發(fā)森林火災(zāi)，火災(zāi)燒毀林地約160公頃[2]，森林火災(zāi)造成林區(qū)內(nèi)大量動植物死亡，對林區(qū)生態(tài)系統(tǒng)造成嚴(yán)重破壞，產(chǎn)生大量的煙塵和溫室氣體污染大氣環(huán)境[3]。此外，火災(zāi)后引發(fā)水土流失、崩塌滑坡、泥石流等不良次生地質(zhì)災(zāi)害，造成巨大的人員傷亡和財產(chǎn)損失[4]。

火災(zāi)后泥石流在成因上與普通泥石流有較大差異，其松散固體物質(zhì)主要是在火災(zāi)后產(chǎn)生，按照啟動條件的不同，可分為表面徑流型和淺層滑坡型[5]。目前，國內(nèi)外對火災(zāi)后泥石流的研究較少，對泥石流的研究多集中在震后型泥石流和普通型泥石流[6-8]，火災(zāi)后泥石流的研究主要集中在火災(zāi)改變土壤的水力性質(zhì)以及增加物源的動態(tài)儲量[9-13]；ColeR[14]通過對比試驗發(fā)現(xiàn)，在火災(zāi)區(qū)域內(nèi)土壤含水量和淺表土的吸附性顯著降低；McGuire[15](2019)指出火災(zāi)后流域爆發(fā)泥石流的臨界閾值降低約35%；雷鳴宇[16]等(2021)研究發(fā)現(xiàn)火災(zāi)后植被根系死亡導(dǎo)致土體抗剪強度下降，從而解釋了火災(zāi)后泥石流的啟動機理。

2018年2月16日，雅江縣惡古鄉(xiāng)馬益西村后山突發(fā)森林火災(zāi)，過火面積20.59 km2，其中嚴(yán)重?zé)龤娣e6 km2，涉及惡古村、下馬益西村。森林火災(zāi)過后，同年6月暴發(fā)泥石流災(zāi)害。本文以洞擁過溝為典型研究對象，探究火災(zāi)后泥石流的運動過程以及形成機制。

1 研究區(qū)概況

洞擁過溝位于甘孜州雅江縣惡古村，屬雅礱江水系一級支流惡古河左岸支流， 距雅江縣城58 km。出露地層由老到新依次為三疊系上統(tǒng)(侏倭組T3zh)地層、第四系(Q)堆積層。地理位置處于川西高原與四川盆地過渡帶，地跨揚子陸塊區(qū)與羌塘-三江造山系兩大一級構(gòu)造單元，新構(gòu)造活動異常強烈，使區(qū)域內(nèi)地勢不斷抬升，兩側(cè)風(fēng)化剝蝕搬運作用顯著，屬構(gòu)造剝蝕侵蝕高山地貌。研究區(qū)位于活動頻繁的鮮水河地震帶附近，地震頻發(fā)，地震基本烈度為7度。

洞擁過溝位于青藏高原東緣，屬于大陸性季風(fēng)高原型氣候。年均氣溫一般在4～11 ℃，最低氣溫-16.0 ℃，最高氣溫為36.1 ℃。年平均降雨量705.2 mm，年降水量最大值為1 001.3 mm，最小值為497.0 mm，年平均蒸發(fā)量1 090 mm。但區(qū)內(nèi)降雨量時間分布不均，主要集中于5～9月。

2 洞擁過溝泥石流基本特征

2.1 泥石流災(zāi)害的歷史及特征

洞擁過溝歷史上多次發(fā)生泥石流。1987年，洞擁過溝因森林火災(zāi)曾暴發(fā)大規(guī)模泥石流，導(dǎo)致溝口11戶民房被埋。1997年7月，泥石流攜帶的大石塊堆積于溝道下游和溝口，堆積平均長度約150 m，平均寬度為100 m，平均厚度約為2 m，初步估算該次泥石流固體沖出量約3×104m3。2018年2月惡古鄉(xiāng)流域發(fā)生大面積深林火災(zāi)，發(fā)生火災(zāi)之前洞擁過溝屬于低頻泥石流溝，但是火災(zāi)發(fā)生后，同年的雨季(5～9月)，該流域多次暴發(fā)暴雨型泥石流災(zāi)害，最為嚴(yán)重的一次發(fā)生在6月中旬，在強降雨作用下(6 h降雨量33.6 mm)泥石流持續(xù)時間約30 min，泥石流攜裹溝口前緣滑坡堆積體沿溝道沖刷，溝口最高處泥石流泥位高度3.5 m以上，造成兩岸農(nóng)戶房屋損毀，護(hù)堤破壞。該次泥石流沖出固體物質(zhì)約1.05×104m3，堆積物沖至惡古河，造成惡古河淤塞，上游形成堰塞湖，嚴(yán)重威脅了溝道下游村鎮(zhèn)及惡古村農(nóng)戶的生命財產(chǎn)安全。

洞擁過溝與主河惡古河呈90°左右的銳角相交，不利于泥石流物質(zhì)排泄；利用形態(tài)調(diào)查法計算出泥石流發(fā)生時惡古河洪峰流量約為100 m3/s，洞擁過溝泥石流在P=5%下泥石流峰值流量為32.27 m3/s，惡古河流量為P=5%泥石流峰值流量的3倍，洞擁過溝泥石流屬高重度的粘性泥石流，泥石流沖出物易淤積，難易被主河河水稀釋并搬運，因此，暴發(fā)泥石流災(zāi)害時，其造成主河惡古河堵塞和潰決的可能性較大。

2.2 泥石流類型特征

根據(jù)現(xiàn)場調(diào)查訪問，1997、2018年洞擁過溝泥石流漿體濃稠，按泥石流災(zāi)害防治工程勘察規(guī)范中泥石流漿體稠度特征表的標(biāo)準(zhǔn)，其容重應(yīng)在1.8～2.3 g/cm3之間，結(jié)合現(xiàn)場配漿試驗，實際容重在1.89 g/cm3，為偏粘性泥石流。根據(jù)記錄，洞擁過溝泥石流暴發(fā)的頻率為低頻泥石流。值得注意的是，2018年森林火災(zāi)后已連續(xù)兩年發(fā)生暴雨泥石流災(zāi)害，而且具有發(fā)育為中頻泥石流的可能，因其發(fā)育于溝谷地帶，屬于溝谷型泥石流。泥石流發(fā)生所需要的松散固體物質(zhì)來源是主溝兩側(cè)的崩滑和坡面侵蝕以及運動過程中溝床底部的溝道物源，屬于滑坡及坡面侵蝕型泥石流。泥石流的激發(fā)條件是短歷時強降雨，屬于暴雨型泥石流。因此，洞擁過溝泥石流屬于粘性—低頻—暴雨—溝谷型—滑坡及坡面侵蝕性類泥石流。

2.2 泥石流堆積特征

對洞擁過溝泥石流堆積區(qū)不同位置的泥石流堆積物進(jìn)行采樣，在室內(nèi)進(jìn)行顆粒級配試驗并繪制顆粒級配曲線，曲線說明洞擁過溝內(nèi)不同位置的水動力條件和沖淤特征。溝道內(nèi)(1號，2號)土石比均較小，縱坡比大，水動力條件強，強大的水動力條件使溝道內(nèi)的細(xì)小顆粒難以停留在流通區(qū)溝道內(nèi)，被泥石流裹挾帶走，僅留下粗顆粒物質(zhì)在溝內(nèi)。流通區(qū)(3號)由于森林火災(zāi)后坡面侵蝕物源、溝道岸坡滑坡物源等聚集，而地表水勢能尚未在向下徑流中經(jīng)巨大的溝道比降迅速增大，其裹挾泥砂的能力相對較弱，溝內(nèi)出現(xiàn)輕微淤積。在下游流通區(qū)(4號)土石比自上游向下游總體呈現(xiàn)逐漸增大的趨勢特征，表明堆積區(qū)水動力條件迅速變?nèi)酰沽黧w中細(xì)顆粒物質(zhì)淤積。

2.3 泥石流運動參數(shù)特征

2.3.1 泥石流流體重度

洞擁過溝火災(zāi)后泥石流容重的計算標(biāo)準(zhǔn)，是根據(jù)2018年6月暴發(fā)的泥石流災(zāi)害現(xiàn)場調(diào)查測定的。測定方法是在洞擁過溝泥石流主要流通堆積區(qū)段，采集泥石流堆積物攪拌配制成泥石流漿體，配置好的漿體要與泥石流災(zāi)害后溝道殘留的泥痕的顏色、形狀和厚度等特征相符，并找到目睹過泥石流暴發(fā)的當(dāng)?shù)卮迕?，將配置的樣品讓其確認(rèn)，然后測量配置漿體的容重，其計算公式為：

(1)

式中γc為泥石流重度(g/cm3)；Gc為配制泥漿重量(g)；V為配制泥漿體積(cm3)。

計算結(jié)果(表1)表明洞擁過溝泥石流重度為1.76～2.02 g/cm3，平均值為1.89 g/cm3。

表1 洞擁過溝泥石流容重計算表

2.3.2 泥石流的流速

泥石流流速是重要的泥石流動力學(xué)性質(zhì)之一，洞擁過溝泥石流流速采用東川公式計算：

Vc=KHc2/3Ic1/5

(2)

式中:Vc為泥石流斷面平均流速(m/s)；H為泥位高度(m);K為泥位縱坡率；Ic為清水河床糙率。

計算結(jié)果表明(表2)，洞擁過溝泥石流平均流速為6.63 m/s。

表2 洞擁過溝泥石流流速計算表

2.3.3 泥石流峰值流量

洞擁過溝泥石流峰值流量按《四川省水文手冊》提供的公式計算，且先計算暴雨洪峰流量，公式如下：

(3)

式中：Qp為暴雨洪峰流量(m3/s)；ψ為洪峰徑流系數(shù)；F為匯水面積(km2)；S為暴雨雨強(mm/h)，P=2%，32.73 mm/h，P=5%，27.6 mm/h ；i為最大暴雨強度(mm/h)；n為暴雨公式指數(shù)；τ為流域匯水時間(h)。

泥石流峰值流量計算：

Qc=(1+φ)Qp·Dc

(4)

φ=(γc-γw)/(γH-γc)

(5)

式中Qc為頻率為P的泥石流洪峰值流量(m3/s)；Qp為頻率為P的暴雨洪水設(shè)計流量(m3/s)：φ為泥石流泥沙修正系數(shù)；γc為泥石流重度(t/m3)；γw為清水重度(t/m3)；γH為泥石流中固體物質(zhì)比重(t/m3)；Dc為泥石流堵塞系數(shù)。

洞擁過溝暴雨洪峰流量計算表見表3，泥石流峰值流量計算結(jié)果見表4。計算結(jié)果表明在溝口位置，P=5%的暴雨洪峰流量為9.02 m3/s，對應(yīng)的泥石流峰值流量為32.27 m3/s。

表3 洞擁過溝暴雨洪峰流量計算表

表4 洞擁過溝泥石流峰值流量計算表

2.3.4 一次泥石流過流總量

一次泥石流過流總量按照蔣忠信編寫的《震后山地地質(zhì)災(zāi)害治理工程設(shè)計概要》[17]提供的計算公式進(jìn)行計算：

(6)

式中Q為泥石流持續(xù)時間(s)，根據(jù)調(diào)查訪問確定;T為一次泥石流過程總量(m3)；QC為泥石流最大流量(m3/s)。

計算結(jié)果(表5)表明，洞擁過溝泥石流溝口處在P=5%時，一次泥石流過流總量為20 446.48 m3。

表5 洞擁過溝一次泥石流過流總量計算表

3 泥石流的成因與機制

泥石流災(zāi)害的形成主要取決于流域內(nèi)物源的數(shù)量。流域內(nèi)動物源的數(shù)量豐富時，在一定的短歷時強降雨條件下，發(fā)生泥石流的可能性較大[18]。內(nèi)外動力地質(zhì)作用為洞擁過溝泥石流的暴發(fā)提供了豐富的松散固體物質(zhì)，如火災(zāi)使洞擁過溝上游形成區(qū)面蝕加劇坡面物源增加；洞擁過溝道上游及其支溝地形陡峻，溝谷縱坡大，是泥石流暴發(fā)的動力勢能條件；溝谷降雨豐富匯水條件好，是泥石流暴發(fā)的水源條件。

3.1 地形地貌條件

洞擁過溝地勢北西高南東底，溝域形態(tài)近似帶狀，溝道匯水面積2.26 km2，溝域縱向長度3.45 km，平均寬度0.6 km，平均縱坡比降為457‰。研究區(qū)內(nèi)高程最高點位于西北側(cè)分水嶺處海拔高度為4 210 m，最低海拔位于洞擁過溝匯入惡古河口處為2 631 m，相對高差1 579 m。研究區(qū)內(nèi)地勢險峻，溝谷縱坡比降大，泥石流形成區(qū)的溝道縱坡比降在400‰以上。因此，為泥石流的暴發(fā)和運動提供了良好的地勢條件。

根據(jù)洞擁過溝的溝谷變化、物源數(shù)量等將其劃分為：清水區(qū)、補給區(qū)、堆積區(qū)(表6)。

表6 洞擁過溝流域分區(qū)表

清水區(qū)高程4 210～3 440 m，區(qū)內(nèi)基巖大多裸露，局部分布有少量冰蝕堆積物，溝源山頂基本無植被覆蓋；溝道二側(cè)溝谷陡峭，崩滑堆積體厚度大，該區(qū)段為泥石流的形成提供了一定量的物源和較為豐富的水源和水動力條件。補給區(qū)高程分布為2 705～3 440 m， 溝谷岸坡陡峭，松散堆積體在區(qū)內(nèi)分布范圍廣，崩塌滑坡等不良地質(zhì)現(xiàn)象在區(qū)內(nèi)零星分布，尤其是在溝道中下段發(fā)育有大型滑坡，為泥石流的發(fā)生和運動提供了豐富的物源基礎(chǔ)。堆積區(qū)高程為2 705～2 631 m，平均縱坡220‰，地形平緩，不利于泥石流堆積物流通。

3.2 物源條件

松散固體物質(zhì)是泥石流暴發(fā)的物源基礎(chǔ)[19-22]。洞擁過溝泥石流松散固體物質(zhì)儲量較豐富，物源類型可分為集中物源、溝道物源、坡面物源。根據(jù)物源的潛在滑動面、膠結(jié)程度、密實度等條件將其劃分為不穩(wěn)定和穩(wěn)定物源。根據(jù)現(xiàn)場勘察，溝域內(nèi)物源總量194.3×104m3，動物源總量67.16×104m3。

此外值得注意的是，由于2018年2月大面積、高強度持續(xù)火燒，林火導(dǎo)致坡面上植被大面積死亡，固坡能力下降，林火形成大量的坡面細(xì)小堆積物，例如灰燼和泥沙，灰燼是火災(zāi)后泥石流的一種特殊物源。林火發(fā)生后，大量的植被及枯枝落葉被燒毀炭化，形成了灰燼層?；覡a層疏松多孔，質(zhì)地較輕，在前期降雨中容易飽和，在暴雨時很容易起動，作為泥石流的一種物源，直接參與泥石流活動。由于高溫烘烤作用，坡表土壤的物理力學(xué)性質(zhì)發(fā)生了較大改變，形成了較松散的泥沙層，這是坡面物源的重要來源，它們也是構(gòu)成火災(zāi)后泥石流直接的物源。此外，火災(zāi)后大量的樹枝和樹根等木材也成為火燒跡地主要的碎屑物質(zhì)。在坡面流水的沖刷作用下，部分細(xì)顆粒在雨水的沖刷作用下容易順斜坡向坡腳搬運成為泥石流坡面物源。洞擁過溝物源類型見表7。

表7 洞擁過溝物源類型表

3.3 泥石流的形成機制

在短歷時強降雨作用下洞擁過溝暴發(fā)泥石流，火災(zāi)后的淺表物源被泥石流夾帶。流域內(nèi)相對高差大，地形陡峭動力勢能強，使洞擁過溝泥石流規(guī)模不斷增大，流速不斷加快。在強大的沖擊侵蝕作用下，泥石流沿途掏蝕流通堆積區(qū)溝岸坡腳，下切侵蝕作用不斷增強，崩塌滑坡物源和溝道內(nèi)松散堆積物不斷參與泥石流的運動過程。

洞擁過溝在2018年深林火災(zāi)未發(fā)生以前，流域植被平均覆蓋率約80%，火災(zāi)發(fā)生后，溝道中上游火災(zāi)區(qū)域基本無植被覆蓋，導(dǎo)致崩滑坡面侵蝕等不良地質(zhì)災(zāi)害強烈發(fā)育，松散固體物質(zhì)增加。根據(jù)現(xiàn)場勘察流域內(nèi)松散固體物源總量為159.48×104m3，可能參與泥石流運動的物源53.26×104m3，溝道內(nèi)沿途泥沙等物源補給充分，加之泥石流形成區(qū)的溝道縱坡比降在400‰以上，相對高差大，溝谷形態(tài)為“V”，因此，內(nèi)外地質(zhì)作用強烈，地勢不斷抬升，兩側(cè)風(fēng)化剝蝕搬運作用顯著。以上因素均為泥石流暴發(fā)的有力條件。按照《泥石流災(zāi)害防治工程勘查規(guī)范》[23]，對影響洞擁過溝泥石流的因素綜合評價，發(fā)生程度屬易發(fā)。

4 結(jié) 論

(1)洞擁過溝泥石流屬于粘性—低頻—暴雨—溝谷型—滑坡及坡面侵蝕性類泥石流，泥石流平均重度為1.89 g/cm3，平均流速為6.63 m/s，P=5%時泥石流一次固體物質(zhì)沖出量為20 446.48 m3；洞擁過溝泥石流P=5%時的泥石流峰值流量為32.27 m3/s ，P=2%時泥石流峰值流量為42.01 m3/s，惡古河流量為P=5%泥石流峰值流量的3倍。

(2)火災(zāi)導(dǎo)致流域內(nèi)崩滑等不良地質(zhì)現(xiàn)象增多，侵蝕加劇，泥石流暴發(fā)所需要的松散固體物質(zhì)增加，根據(jù)現(xiàn)場勘察流域內(nèi)松散固體物源總量為159.48×104m3，可為泥石流運動提供物質(zhì)補給的物源約53.26×104m3。

(3)在短歷時降雨條件激發(fā)下，洞擁過溝暴發(fā)泥石流的可能性較大，惡古河輸沙能力一般，不利于泥石流物質(zhì)的排泄，泥石流沖出物易淤積，暴發(fā)泥石流災(zāi)害時其造成主河惡古河堵塞和潰決的可能性較大。