摘"要：探討極端強降雨對于泥石流易發(fā)性的影響，為泥石流災(zāi)害預(yù)警與防治提供支持。以北京市曹家坊泥石流溝為研究區(qū)，選取泥砂沿程補給長度比、沿溝松散物儲量、產(chǎn)沙區(qū)松散物平均厚度等因子，結(jié)合遙感手段參照《地質(zhì)災(zāi)害危險性評估技術(shù)規(guī)范》，對研究區(qū)在“23·7”強降雨前后的泥石流易發(fā)性進行打分并評價。結(jié)果表明：泥沙沿程補給長度比、沿溝松散物儲量、產(chǎn)沙區(qū)松散物平均厚度、河溝堵塞程度、河溝近期一次變幅等5項因子更易受到極端強降雨的影響；在沖於作用影響下，曹家坊泥石流溝的易發(fā)性由“低易發(fā)”升級為“中易發(fā)”，建議對淤積物源進行清理。

關(guān)鍵詞：泥石流；易發(fā)性評價；“23·7”強降雨；曹家坊泥石流溝

Influences of \"23·7\" heavy rainfall on the hidden danger of debris flow at Caojiafang， Beijing

ZHAO Danning¹， JIAO Runcheng¹， YANG Chun²

（1.Beijing Institute of Geological Hazard Prevention， Beijing 100120， China；

2.Twenty-First Century Aerospace Technology Co.， Ltd.， Beijing 100096， China）

Abstract： To explore the influence of extremely heavy rainfall on the susceptibility of mudslides， and provide support for early warning and prevention of mudslide disasters， this paper takes Caojiafang mudslide ditch in Beijing as the study area and selects 15 factors such as sediment recharge length ratio， loose material storage along the ditch， and the average thickness of loose material in the sand-producing area. The standardized scoring method is applied in combination with remote sensing to score the susceptibility and make comparisons of the mudslides before and after the heavy rainfall of \"23·7\" in the study area. The results show that 5 factors， including the length ratio of along-channel recharge， loose material storage along the gully， the average thickness of loose material in the sand-producing area， the degree of blockage of the gully， and the recent variation of the gully， are more susceptible to the influence of extreme rainfall. Under the decisive influence of the flushing effect， the susceptibility of Caojiafang mudslide gully is upgraded from \"low susceptibility\" to \"medium susceptibility\". The susceptibility of Caojiafang mudslide ditch is also upgraded from \"low susceptibility\" to \"medium susceptibility\" under the decisive influence of alluvial action.

Keywords： mudslide; vulnerability assessment; \"23·7\" heavy rainfall; Caojiafang mudslide ditch

泥石流是一種山區(qū)常見的地質(zhì)災(zāi)害，它的發(fā)生嚴重影響了社會經(jīng)濟的發(fā)展，給人民生命財產(chǎn)的安全帶來了巨大威脅（李寧等，2022），因此，泥石流易發(fā)性的分析與研究，對于泥石流災(zāi)害的監(jiān)測預(yù)警以及防治有著重要的意義（南赟等，2020；鮑進進，2019）。

目前泥石流易發(fā)性評價的方法主要分為兩類：定性評價和定量評價。定性評價主要基于個人的經(jīng)驗和認知，主觀性強，人為因素干擾較大；定量評價主要基于統(tǒng)計學(xué)，并與經(jīng)驗相結(jié)合，如多因子評價法、GIS信息模型法、層次分析法、灰色關(guān)聯(lián)度法、人工神經(jīng)網(wǎng)絡(luò)評價法等，客觀性較強（陳攀等，2022；崔志超等，2020；李坤等，2022；羅守敬等，2021；張群等，2023）。在進行泥石流易發(fā)性評價時，確定易發(fā)性主控因子及其貢獻率十分重要，這其中還有諸多問題需要解決（劉永垚等，2018）。

在易發(fā)性評價單元上，現(xiàn)有的評價方法中評價單元一般分為３種，即柵格單元、自然流域（王峰等，2023；王向前等，2023）和行政區(qū)域。柵格單元的優(yōu)勢是可以利用GIS處理矩陣數(shù)據(jù)，十分便捷，缺點是破壞了流域單元的完整性，不利于結(jié)合地形地貌等環(huán)境因素；以行政區(qū)域作為評價單元有利于為各級政府和管理部門進行防災(zāi)救災(zāi)決策提供科學(xué)依據(jù)和技術(shù)支持；以自然流域作為評價單元，可以保證流域的完整性，利于易發(fā)性判別因子的獲取，還可以充分結(jié)合周邊地質(zhì)環(huán)境要素（黃來源，2022）。隨著GIS技術(shù)的發(fā)展，基于流域單元的泥石流易發(fā)性評價已成為主流。

北京西山地區(qū)富含大量煤炭資源，曾經(jīng)存在大量采煤區(qū)，截止2020年煤礦都已關(guān)閉。煤礦開采初期占用大量土地，造成地表植被破壞，導(dǎo)致土壤持水能力減弱，開采過程中會在露天堆積大量煤矸石，一旦有暴雨發(fā)生，流域內(nèi)容易發(fā)生泥石流（王海芝等，2022；王文等，2022）。針對此類泥石流隱患，北京市政府先后投入大量人力物力對其進行了治理。

在“23·7”強降雨過程中，北京市很多泥石流隱患發(fā)生了溝道侵蝕和沖淤，在一定程度上改變了泥石流的性質(zhì)和特征，對它的易發(fā)性產(chǎn)生了影響。這些泥石流隱患迫切需要對其易發(fā)性進行重新評估。

本文以曹家坊泥石流溝為研究對象，對其災(zāi)前和災(zāi)后的物源分布情況、溝道動力條件、植被覆蓋等情況進行提取和對比，探索強降雨過程對泥石流溝的易發(fā)性造成的影響。

1 "研究區(qū)概況和數(shù)據(jù)源

1.1 "研究區(qū)概況

1）地質(zhì)背景

曹家坊泥石流溝位于北京市房山區(qū)史家營鄉(xiāng)，處于中山區(qū)域，此區(qū)域一般山體陡峭，溝道落差和縱比降較大，較利于泥石流的形成。

根據(jù)《北京市區(qū)域地質(zhì)志》（北京市地質(zhì)礦產(chǎn)局，1991），該區(qū)第四紀時期的上升速度明顯要高于第三紀，且從中更新世晚期開始表現(xiàn)為快速上升，最快抬升速度接近1 mm·a^-1，依據(jù)北京市地質(zhì)志里關(guān)于新構(gòu)造運動的論述，推測今后山區(qū)仍將有加速上升的發(fā)展趨勢。因此，該泥石流溝目前處于抬升區(qū)。

曹家坊溝內(nèi)主要發(fā)育侏羅紀地層，大部分地區(qū)巖性以較硬的砂礫巖和火山巖為主；中上游地區(qū)巖性軟硬相間，其局部發(fā)育窯坡組含煤地層。綜合來看：曹家坊溝內(nèi)巖石風(fēng)化程度較低，巖性較硬；斷裂構(gòu)造較為發(fā)育，包含有多處斷層破碎帶（圖1）。

曹家坊泥石流溝內(nèi)主要物源包括煤矸石堆（未治理和治理后2種類型）、自然松散堆積物、危巖體、道路切坡、旱地覆土、災(zāi)害體、局部水體以及人類工程形成的堆積物等。

曹家坊溝歷史沒有發(fā)生泥石流災(zāi)害的記載。

2）工程治理情況

曹家坊礦區(qū)由于開采歷史較長，區(qū)內(nèi)森林植被損毀，水土流失、采空塌陷問題突出，山體崩塌、泥石流等地質(zhì)災(zāi)害易發(fā)，自2010年起政府相關(guān)部門對曹家坊礦區(qū)開展了全面生態(tài)修復(fù)，修復(fù)面積160多公頃，礦山修復(fù)治理面積約69.84 hm2，修建行洪渠4000余米，對溝道內(nèi)物源進行了攔擋、穩(wěn)固以及植被恢復(fù)，植被覆蓋率大幅度提升（圖2）。

1.2 "數(shù)據(jù)源

采用時相分別為2023年5月份和2023年10月份的數(shù)字正射影像圖（Digital Orthophoto Map，DOM）數(shù)據(jù)，分辨率都是0.1 m。獲取的數(shù)據(jù)無霧、霾覆蓋，色調(diào)均勻，反差適中，云量覆蓋率小于"0.5%，且不壓蓋目標地物，滿足本次工作的要求。同時收集了研究區(qū)的大比例尺航攝數(shù)據(jù)用來制作高精度數(shù)字高程模型（Digital Elevation Model，DEM）。

2 "研究方法

在研究區(qū)“23·7”強降雨前后兩期高精度DOM基礎(chǔ)上，通過建立遙感解譯標志，提取地質(zhì)災(zāi)害隱患背景信息和泥石流溝物源分布，崩塌隱患以及受威脅對象分布信息，結(jié)合高精度DEM提取溝道動力綜合條件因子，并分別開展災(zāi)害隱患易發(fā)性評價，通過對比分析，探索強降雨過程對泥石流溝的易發(fā)性造成的影響。

2.1 "源數(shù)據(jù)處理

1）幾何校正與配準

采用特征點匹配法對兩期影像進行配準，提取兩期影像中相似的特征點，如角點、邊緣，然后通過相似點間的相似性度量。本次共選擇10個特征點，采用一階多項式變換（仿射），殘差值控制在1個像素以下（表1）。

2）DEM制作

利用高精度定位定姿系統(tǒng)（position and orientation system，POS）數(shù)據(jù)進行生產(chǎn)，自動DEM匹配采用DEM制作模塊可以完成0.5 m格網(wǎng)間距DEM自動匹配，自動匹配DEM的效果如圖3所示。

DEM編輯在二三維一體化的環(huán)境下完成，通過點線面編輯等工具以及立體檢查完成DEM的局部編輯，提高了DEM精度。

2.2 "遙感解譯

利用大比例尺航攝數(shù)據(jù)，綜合運用遙感地質(zhì)解譯和遙感信息提取技術(shù)開展地質(zhì)災(zāi)害遙感綜合解譯。

1）解譯內(nèi)容

地形背景解譯：在充分利用研究區(qū)已有研究成果和資料的基礎(chǔ)上，通過遙感解譯補充解譯與地質(zhì)災(zāi)害發(fā)育有關(guān)的地形信息。

泥石流物源信息提?。豪萌藱C交互信息提取方法，通過二維與三維相結(jié)合的方式獲取每條泥石流隱患溝內(nèi)物源類型、分布情況。

溝道動力綜合條件信息提?。夯诤綌z數(shù)據(jù)，利用三維空間分析手段，通過高精度DEM數(shù)據(jù)進行匯水和徑流模擬計算，對泥石流溝流域范圍、流域形態(tài)、溝系發(fā)育、溝床比降、流域山坡坡度及溝道形態(tài)等情況進行自動提取。

崩塌信息提?。夯诟叻謹?shù)據(jù)，借助數(shù)字三維模型，提取崩塌隱患大型裂縫和危巖分布情況。

2）解譯方法

本次研究中，結(jié)合不同提取目標物的特點，綜合運用人機交互解譯和計算機自動信息提取兩種方法，遵循“從已知到未知、先易后難、逐步解譯”的原則，充分利用各種分析推理方法進行解譯。

建立研究區(qū)各類信息提取目標物解譯標志，進行人機交互解譯（表2）。

水體分為河流、坑塘水面及水庫等，是泥石流隱患溝中發(fā)育的物源類型之一。

坑塘水面一般為淡藍至深藍色調(diào)，隨水體的加深，色調(diào)逐漸加深，內(nèi)部色調(diào)較均勻，與周邊地物邊界清晰，易于判讀

結(jié)合高精度DEM能夠精確獲取泥石流水動力綜合條件因子，如河流縱坡、河岸山坡坡度、產(chǎn)沙區(qū)橫斷面、流域面積、流域相對高差等（圖4）。

3 "強降雨泥石流易發(fā)性影響分析

參照DB 11/T 893－2021《地質(zhì)災(zāi)害危險性評估技術(shù)規(guī)范》G.4“泥石流易發(fā)程度等級評判表”對曹家坊溝泥石流易發(fā)性進行評價，評價過程選取15個因子，按照暴雨工況對評價因子的可影響程度對因子進行分類并打分，最終按照總分“gt;115”“86～115”“l(fā)t;86”將評價結(jié)果分為“高易發(fā)”“中易發(fā)”“低易發(fā)”3個等級。

3.1 "不受暴雨影響的因子

1）不良地質(zhì)體嚴重程度

通過解譯，暴雨前共在泥石流溝內(nèi)發(fā)現(xiàn)崩塌、滑坡、危巖體等不良地質(zhì)體5處，自然成因和人工切坡成因均有。在“23·7”暴雨后未發(fā)現(xiàn)新增的地質(zhì)災(zāi)害。

依據(jù)區(qū)域地質(zhì)背景調(diào)查，研究區(qū)內(nèi)的砂礫巖和火山巖質(zhì)地偏硬，并且風(fēng)化不強烈，在暴雨的影響下，不容易產(chǎn)生新的不良地質(zhì)體。而區(qū)域內(nèi)工程治理比較完善，也極大降低了不良地質(zhì)體產(chǎn)生的幾率。

2）溝口泥石流堆積活動程度

在溝道治理工程中，曹家坊溝溝口的堆積物已經(jīng)被改造，建設(shè)成為居民地，河道受堆積物和人類工程活動影響向西南側(cè)偏移（圖5），在暴雨后主河道沒有受到泥石流的影響，這次暴雨形成的沖淤物也沒有沖出溝口對主河道造成影響，此項因子沒有變化。

3）巖性影響

曹家坊溝內(nèi)發(fā)育巖性較硬，風(fēng)化程度不高，此項因子不受強降雨的影響。

4）泥石流溝水動力條件因子

河溝縱坡、溝岸山坡坡度、產(chǎn)沙區(qū)橫斷面、流域面積、流域相對高差這5項屬于泥石流溝水動力條件因子，是泥石流形成的關(guān)鍵因素之一。

利用DEM結(jié)合高精度DOM分析計算可得：曹家坊溝流域面積3.57 km²，主溝流動方向約為140°；溝岸山坡平均坡度約為34.43°；主溝長約3.75 km，彎曲系數(shù)為1.10，溝系密度6.7 km·km^-2；最大高差為1 050 m，平均坡降約為28.6%；產(chǎn)沙區(qū)溝槽橫斷面呈V字型（圖6）。

這5項因子在“23·7”強降雨后都不發(fā)生變化。

5）區(qū)域構(gòu)造影響程度

由地質(zhì)背景可知調(diào)查區(qū)域范圍位于抬升區(qū)，7級地震烈度帶內(nèi)，斷裂構(gòu)造發(fā)育。此因子暴雨后不發(fā)生變化。

6）植被覆蓋率

由于暴雨沖刷導(dǎo)致泥土稀松，使得部分植被被損毀，通過土地利用類型自動分類和植被指數(shù)兩種計算方法，分別獲得流域的植被覆蓋率，再計算二者平均值，得出“23·7”強降雨前植被覆蓋率為81.52%，強降雨后植被覆蓋率為75.65%。

3.2 "受暴雨影響的因子

1）動力綜合條件因子

溝道堵塞程度和河溝最近一次沖淤變幅都屬于動力綜合條件。

原溝道經(jīng)過治理后，主河道流通區(qū)已經(jīng)十分通暢，“23·7”強降雨后溝道內(nèi)中上游煤矸石大量沖淤，造成溝道內(nèi)多處發(fā)生嚴重堵塞（圖7、圖8）。

在“23·7”強降雨之前，曹家坊泥石流溝最近一次較明顯沖淤發(fā)生于2012年，沖淤變幅程度相對較低，普遍不超過0.1 m；“23·7”強降雨后，溝道內(nèi)物源淤積明顯（圖9），有些位置物源已經(jīng)漫過人工泄洪渠，平均厚度達到4 m。因此，河溝最近一次沖淤變幅接近4 m。

2）物源變化相關(guān)因子

泥砂沿程補給長度比、沿溝松散物儲量、產(chǎn)沙區(qū)松散物平均厚度等3項因子都與物源變化相關(guān)。

泥砂沿程補給長度比是指泥沙沿程補給累計長度與主溝長度比。根據(jù)泥石流溝遙感解譯成果，通過圖上量測可得曹家坊溝主溝長度，泥沙沿程補給長度則從流域西北側(cè)分水嶺開始測量，分水嶺開始至溝谷直到溝口，對其中可能對泥石流產(chǎn)生補給的區(qū)段（如第四系松散沉積物分布區(qū)段、基巖風(fēng)化強烈的區(qū)段、礦渣與生活固廢等人工堆積物分布區(qū)段等，但無補給源的裸露基巖分布區(qū)段應(yīng)去除）分別量測其長度，然后累計計算總長度。

曹家坊泥石流溝經(jīng)測量泥石流溝主溝長度約為3"740 m，在“23·7”強降雨前綜合量測泥沙補給長度約為2 169 m，“23·7”強降雨后，因大量第四系松散沉積物剝落，溝道兩側(cè)植被覆蓋層被雨水沖刷侵蝕，新增物源主要為廢棄的煤矸石，災(zāi)后新增物源范圍位于流域中段（圖10），沿主溝長度約460 m，泥石流補給區(qū)段明顯加長，綜合量測泥沙補給長度約為2 618 m。

經(jīng)計算補給長度比在“23·7”強降雨前泥砂沿程補給長度比約為58%，強降雨后可達到約70%。

“23·7”強降雨后，溝道內(nèi)沖出大量物源，其中以原治理后的煤矸石堆為主，輔以少量旱地覆土和自然松散堆積物，通過遙感解譯，對曹家坊溝內(nèi)各類物源儲量進行計算，估算結(jié)果顯示該溝內(nèi)主要物源儲量總計可達19.69×10⁴ m³。

“23·7”強降雨前，主要產(chǎn)沙區(qū)位于泥石流溝中上游，物源為煤矸石堆和松散堆積物，通過影像分析和野外調(diào)查，各處煤矸石堆、松散物的堆積厚度多數(shù)在1 m以下；“23·7”強降雨后，大量煤矸石堆及少量松散物被沖出，產(chǎn)沙區(qū)堆積物厚度最厚處可達4 m，最薄處也在1 m以上。

4 "泥石流溝易發(fā)程度評價

由前文可得，在對曹家坊泥石流溝進行易發(fā)性評價的15個因子中，有5項在“23·7”強降雨后發(fā)生了變化，它們分別是泥沙沿程補給長度比、沿溝松散物儲量（×10⁴"m³·km^-2）、產(chǎn)沙區(qū)松散物平均厚度、河溝堵塞程度、植被覆蓋率。

泥沙沿程補給長度比由58%變?yōu)?0%；沿溝松散物儲量（×10⁴"m³·km^-2）由0.95變?yōu)?.5；產(chǎn)沙區(qū)松散物平均厚度由1"m變?yōu)? m；河溝堵塞程度由沒有堵塞變?yōu)閲乐囟氯?；植被覆蓋率由81.52%變?yōu)?5.65%（表3）。

根據(jù)泥石流溝的易發(fā)程度進行評價，結(jié)合泥石流溝基本特征調(diào)查成果，可知極端降雨可以影響，泥沙沿程補給長度比、沿溝松散物儲量、產(chǎn)沙區(qū)松散物平均厚度、河溝堵塞程度、河道近期一次變幅等因子，導(dǎo)致泥石流易發(fā)程度增大，如曹家坊溝泥石流在“23·7”強降雨前易發(fā)程度評價為“低易發(fā)”，在“23·7”強降雨后易發(fā)性評價為“中易發(fā)”（表3），因此建議在暴雨過后，對存在溝道兩側(cè)植被覆蓋層被雨水沖刷侵蝕明顯、溝道內(nèi)沖出大量物源等現(xiàn)象的泥石流溝，及時開展易發(fā)性的重新評估。

5 "結(jié)論

1）經(jīng)研究可得，對泥石流易發(fā)性進行評價時需要提取的因子中，泥沙沿程補給長度比、沿溝松散物儲量、產(chǎn)沙區(qū)松散物平均厚度、河溝堵塞程度、河溝近期一次變幅等5項因子更易受到強降雨的影響。

2）曹家坊溝在“23·7”強降雨之后，其易發(fā)性由“低易發(fā)”升級為“中易發(fā)”。其中沖淤作用起到了決定性作用，建議對其淤積物源進行清理或固定。

3）由于泥石流溝下游有酒店、景區(qū)和村莊，且溝道內(nèi)已經(jīng)發(fā)生沖淤，說明目前治理工程效果不足以應(yīng)對“23·7”這種極端強降雨情況，建議在經(jīng)過充分勘查研究后，進一步提升工程防治的等級。

參考文獻

鮑進進，2019. 北京西山地區(qū)泥石流危險性評價方法體系研究［D］.中國地質(zhì)大學(xué)（北京）.

北京市地質(zhì)礦產(chǎn)局，1991.北京市區(qū)域地質(zhì)志［M］.北京：地質(zhì)出版社

陳攀，葛永剛，孫慶敏，梁馨月，2022.基于小流域單元的泥石流易發(fā)性評價［J］.科學(xué)技術(shù)與工程，22（29）：12 764-12 771.

崔志超，王俊豪，崔傳峰，潘國營，2020.基于層次分析法和模糊數(shù)學(xué)相結(jié)合的甘肅東鄉(xiāng)八丹溝泥石流易發(fā)性評價［J］.中國地質(zhì)災(zāi)害與防治學(xué)報，31（1）：44-50.

黃來源，2022.基于流域形態(tài)特征分析的北京北山泥石流易發(fā)性評價因子研究［J］.城市地質(zhì)，17（4）：385-393.

李坤，趙俊三，林伊琳，周豹，2022.基于不同斜坡單元劃分方法和BP神經(jīng)網(wǎng)絡(luò)的泥石流易發(fā)性評價［J］.測繪通報"（8）：68-74.

李寧，唐川，史青云，龔凌楓，亓星，2022.九寨溝震區(qū)“6·21”泥石流成因與致災(zāi)機制研究［J］.工程地質(zhì)學(xué)報，30（3）：740-750.

劉永垚，第寶鋒，詹宇，"STAMATOPOULOS C A，"2018.基于隨機森林模型的泥石流易發(fā)性評價：以汶川地震重災(zāi)區(qū)為例［J］.山地學(xué)報，36（5）：765-773.

羅守敬，王珊珊，付德荃，2021.北京山區(qū)突發(fā)性地質(zhì)災(zāi)害易發(fā)性評價［J］.中國地質(zhì)災(zāi)害與防治學(xué)報，32（4）：126-133.

南赟，曹穎，李巖，"2020. 新形勢下北京市突發(fā)地質(zhì)災(zāi)害防治工作思路探析［J］. 城市地質(zhì)，"15（3）："233-238.

王峰，楊帆，江忠榮，吳鄂，汪冠，2023.基于溝域單元的康定市泥石流易發(fā)性評價［J］.中國地質(zhì)災(zāi)害與防治學(xué)報，34（3）：145-156.

王海芝，曾慶利，許冰，胡福根，于淼，2022.北京“7·21”特大暴雨誘發(fā)的地質(zhì)災(zāi)害類型及其特征分析［J］.中國地質(zhì)災(zāi)害與防治學(xué)報，33（2）：125-132.

王文，馬超，張志，張巖，路璐，2022.北京云蒙山片區(qū)“7·16”群發(fā)性暴雨泥石流特征研究［J］.山地學(xué)報，40（4）：626-636.

王向前，丁明濤，何雨楓，2023.基于RF模型的甘孜藏族自治州鄉(xiāng)城縣泥石流易發(fā)性評價［J］.災(zāi)害學(xué)，38（3）：222-226.

張群，馮輝，賈三滿，張沁瑞，賈磊，2023.基于CF與Logistic回歸模型耦合的地質(zhì)災(zāi)害易發(fā)性評價［J］.城市地質(zhì)，18（1）：17-25.