吳季寰，張春山，楊為民，孟華君，郭 涵，萬飛鵬

（1.中國地質(zhì)科學(xué)院地質(zhì)力學(xué)研究所，北京 100081；2.自然資源部活動構(gòu)造與地質(zhì)安全重點實驗室，北京 100081；3.自然資源部陜西寶雞地質(zhì)災(zāi)害野外科學(xué)觀測研究站，北京100081；4.中國地質(zhì)大學(xué)（北京），北京 100083）

引言

青藏高原東北緣地處西秦嶺構(gòu)造帶西延［1］，位于南北向地震帶與祁連山地震帶的交匯區(qū)，活動斷裂發(fā)育，地震活動性強，山地災(zāi)害頻發(fā)。坐落在此的山地城鎮(zhèn)，如甘肅岷縣、武都等地區(qū)，多屬泥石流高發(fā)地［2，3］，泥石流災(zāi)害已嚴(yán)重制約區(qū)域經(jīng)濟社會發(fā)展。

甘肅岷縣朱家溝是青藏高原東北緣典型的小流域泥石流溝，歷史上多次發(fā)生泥石流災(zāi)害，如1982年“8·8”朱家溝泥石流堵斷洮河中上游，形成長達(dá)15 km的回水區(qū)，河水越過防洪堤侵入岷縣城區(qū)，淹沒建筑物1 700余棟，毀壞房屋600余間，造成100余人傷亡，并導(dǎo)致下游約148.67 hm2農(nóng)田絕收減產(chǎn)。近年來，受2008年汶川地震和2013年“7·22”岷縣-漳縣Ms6.6地震影響，朱家溝內(nèi)物源劇增，導(dǎo)致該小流域泥石流危險性急劇升高，極有可能再次爆發(fā)大規(guī)模泥石流災(zāi)害事件，極大威脅流域內(nèi)及城區(qū)居民生命財產(chǎn)安全，因此亟需開展相應(yīng)的泥石流災(zāi)害危險性評估，為防災(zāi)減災(zāi)提供科學(xué)依據(jù)。

泥石流危險評價評價方法大致分為多因子疊加法［4-8］和數(shù)值模擬辦法［9-12］兩類，二者區(qū)別在于：1）評價方式不同：多因子疊加法以泥石流影響因素為數(shù)據(jù)基礎(chǔ)，對流域進行綜合評價分區(qū)，直接體現(xiàn)流域整體危險度，數(shù)值模擬法基于對流域內(nèi)所爆發(fā)泥石流的泥深與流速的模擬，著重于對泥石流動態(tài)過程的反演，并以此為依據(jù)進行危險性分區(qū)；2）適用范圍不同：多因子疊加法適用于大范圍區(qū)域，而數(shù)值模擬辦法更適用于小流域或具體溝谷；3）評價目的不同：多因子疊加法旨在評價各流域于未來階段暴發(fā)泥石流的可能性，數(shù)值模擬法注重將泥石流強度與重現(xiàn)周期相結(jié)合［11］，模擬不同頻率下泥石流的運動過程，預(yù)測泥石流的運動特征及堆積區(qū)域。

本文結(jié)合以上兩類方法，通過詳細(xì)現(xiàn)場調(diào)查和資料分析，采用熵權(quán)與變異系數(shù)融合的多因子疊加法開展朱家溝泥石流易發(fā)性評價；在此基礎(chǔ)上，進行強降雨條件（P=1.6%，P=1%）下朱家溝及其高易發(fā)支溝的泥石流運動數(shù)值模擬，并與野外科考資料對比驗證，最后以流速、泥深、堆積范圍為基礎(chǔ)進行危險性區(qū)劃，該研究可為當(dāng)?shù)啬嗍鳛?zāi)害防治提供依據(jù)。

1 流域泥石流概況

朱家溝坐落于甘肅岷縣城區(qū)北部，是洮河一級支溝，流域面積約14 km2，主溝長度約6.7 km，與縣城一河之隔。流域位于青藏高原東北緣，東昆侖斷裂帶和西秦嶺北緣斷裂帶的應(yīng)變傳遞和構(gòu)造轉(zhuǎn)換的中間過渡區(qū)，屬中強震多發(fā)區(qū)，基本地震動峰值加速度為0.15 g，溝內(nèi)哈樹—麻路灘斷裂為高角度逆沖兼右旋走滑斷裂，傾向北東，傾角約55～70°，斷層上盤以二疊系炭板巖夾石英砂巖、粉砂巖為主，表層中風(fēng)化，風(fēng)化殼厚度2～3 m，屬中等硬度塊狀碎屑巖組，下盤以三疊系板巖、砂質(zhì)板巖、粉砂巖為主，表層強風(fēng)化，風(fēng)化殼厚度5～8 m，屬軟弱層狀碎裂巖。強表層巖土體受地震擾動，結(jié)構(gòu)疏松，崩塌滑坡的觸發(fā)閾值降低，在地震觸發(fā)效應(yīng)下，崩滑發(fā)育，溝道中上游淤積嚴(yán)重，具備豐富的泥石流固體物源條件。

朱家溝地勢北高南低，屬西秦嶺北部構(gòu)造侵蝕剝蝕中低山地貌，海拔介于2 276～3 017 m，相對高差741 m，主溝縱比降約72～110‰，溝道上窄下寬，自北到南曲形延伸，溝道斷面自上而下由“V”型向“U”型過渡，兩岸分布沖溝33條，支溝17條，其中郎家溝、松樹溝、朱林溝和哈固溝是朱家溝的主要支溝。流域地處溫帶半濕潤與高寒濕潤氣候過渡帶，多年平均降水量560.8 mm，降水集中于5～9月，占全年降水總量的78%以上，且多為夏季短時強降雨，受其激發(fā)泥石流活動頻繁，嚴(yán)重威脅城區(qū)及流域內(nèi)居民生命財產(chǎn)安全。

綜上，朱家溝一直是岷縣實施泥石流防災(zāi)減災(zāi)的重要地區(qū)。在岷縣特大冰雹山洪泥石流災(zāi)后重建和岷縣-彰縣地震災(zāi)后重建時，分別在朱家溝、松樹支溝、郎家支溝、王家支溝溝口修建排導(dǎo)槽，在松樹溝修建攔擋壩3道以及在王家支溝修建攔擋壩1道。

2 基于熵權(quán)與變異系數(shù)融合算法的泥石流易發(fā)性評價

2.1 融合算法原理

信息熵是用于度量事物的不確定性［13］，在地質(zhì)災(zāi)害領(lǐng)域表達(dá)為影響因子攜帶信息熵的大小與災(zāi)害發(fā)生的概率成反比。通過信息熵計算熵權(quán)的過程如下：

（1）建立標(biāo)準(zhǔn)化評價矩陣：

式中m、n、x ij分別代表評價對象數(shù)、評價因子數(shù)以及第i個評價對象的第j個因子的標(biāo)準(zhǔn)值。

（2）確定第j個因子的信息熵權(quán)：

式中e j為第j個因子的信息熵，w1j為第j個因子的熵權(quán)。

變異系數(shù)是一個表達(dá)概率分布離散程度的歸一化量度，在評價體系中，某因子的變異系數(shù)愈大，則該因子愈能代表體系內(nèi)部的差異性，對評價結(jié)果影響愈大。第j個因子變異系數(shù)權(quán)重如下：

式中Cv j為第j個因子的變異系數(shù)，w2j為第j個因子的變異系數(shù)權(quán)。

相對熵（式4），可用于衡量兩組離散向量之間差異。

根據(jù)相對熵的定義，引入融合向量W=(w1,w2,…,w n)與待評價向量組Wp×n，使該向量與熵權(quán)法權(quán)重向量以及變異系數(shù)法權(quán)重向量距離之和最?。?3］。建立距離規(guī)劃模型以求解融合權(quán)重：

式中p、n的物理意義分別為災(zāi)害評價模型數(shù)以及災(zāi)害評價因子數(shù)，w i為第j個因子的融合權(quán)重。

2.2 泥石流易發(fā)程度的影響因子

從泥石流形成條件：降水、溝谷、物源三方面出發(fā)，選取影響泥石流發(fā)育的主要因素：形狀系數(shù)、流域面積、溝谷形心至斷層距離、溝谷密度、溝道長度、主溝縱比降、溝口單位寬度清水流量峰值、物源分布、相對高差、流域平均坡度、植被覆蓋率，淺表層工程巖組共計12個評價因子，各因子實際取值見表1。

表1 研究區(qū)各因子實際取值Table 1 The actual value of each factor in the study area

由于研究區(qū)范圍較小，同期次降雨下各支溝雨量接近，故降雨量不作為評價因子。但因一次性短歷時、強降雨形成泥石流過程中，流域內(nèi)雨水短期入滲、蒸發(fā)較少，多以地表徑流形式激發(fā)物源啟動，并最終從溝口流出，所以溝口處單位寬度清水流量峰值可以直觀體現(xiàn)該流域產(chǎn)匯流能力。流域附近歷史最大日降雨量為61.5 mm（2010年5月10日），據(jù)岷縣氣象局2020年6～8月逐日降雨量資料，7月18日（以下簡稱7·18）59.3 mm的降雨量為2020年度當(dāng)?shù)赜昙咀畲蠼邓?。本文基于FLO-2D的Rain模塊，降雨分布如圖2，選擇適用于西部山區(qū)小流域的SKK模型［14］，進行朱家溝“7·18”強降雨-徑流模擬，以提取各支溝溝口清水流量過程峰值，見表1。

圖2 朱家溝“7.18”降雨曲線圖Fig.2 Zhujia gully'7.18'rainfall curve

依據(jù)行業(yè)規(guī)范［15］及相關(guān)研究成果［6.16］，將各因子分為：強（IV）、中（III）、弱（II）及不發(fā)育（I）四級，并將分級結(jié)果代入融合算法計算權(quán)重。分級標(biāo)準(zhǔn)與權(quán)重計算結(jié)果見表2，分級結(jié)果見圖3。

圖3 因子易發(fā)程度量化分級Fig.3 Quantitative classification of factors’susceptibility

表2 評級因子發(fā)育等級及權(quán)重Table 2 Development grade and weight of rating factors

2.3 泥石流易發(fā)程度評價結(jié)果

經(jīng)多因子疊加計算，朱家溝內(nèi)微流域泥石流易發(fā)程度評價結(jié)果介于［0.424，0.624］，利用統(tǒng)計學(xué)中的自然間斷點法進行分類（圖4）：高易發(fā)泥石流溝（0.526，0.624］、中易發(fā)泥石流溝（0.464，0.526］、低易發(fā)泥石流溝［0.424，0.464］。結(jié)合2020年8月野外實地調(diào)查分析（見圖1），流域內(nèi)主要支溝：哈固溝、朱林溝、松樹溝和郎家溝的均為泥石流高易發(fā)溝。

圖1 朱家溝流域地質(zhì)簡圖Fig.1 Geological schematic map of Zhujiagully watershed

圖4 朱家溝泥石流易發(fā)程度評價結(jié)果圖Fig.4 Assessment results of debris flow susceptibility in Zhujia gully

3 基于FLO-2D的朱家溝泥石流危險性評價

綜合考慮朱家溝內(nèi)各支溝泥石流易發(fā)程度和獲取的ALSO12.5 m數(shù)字高程數(shù)據(jù)（DEM）的精度，選取哈固溝、朱林溝、松樹溝、郎家溝等高易發(fā)泥石流支溝以及主溝開展評價。

3.1 工況設(shè)置及計算參數(shù)

FLO?2D以非牛頓流體和中央有限差分法為基礎(chǔ)，迭代計算流體運動控制方程，模擬洪水、泥石流運動過程［12］。計算條件設(shè)置如下：

（1）降雨條件：根據(jù)國家氣象局提供的1981-2010年岷縣逐日降水?dāng)?shù)據(jù)，利用指數(shù)外推法對岷縣地區(qū)年降雨量頻率曲線（圖5）進行擬合（式6），可知雨量59.3 mm的2020年度朱家溝流域最大單日降雨：“7·18”強降雨，年降雨頻率為1.6%，屬50年一遇暴雨事件，設(shè)置此次降雨作為模擬條件。此外，根據(jù)式6選取單日雨量99.38 mm作為百年一遇極端降水條件進行模擬。

圖5 岷縣地區(qū)年降雨量頻率擬合曲線Fig.5 Fitting curve of annual frequency of rainfall in Minxian County

式中：R為降雨量，P為年降雨頻率，擬合公式置信概率為95%。

（2）有無防治工程條件：歷經(jīng)“5·10”特大冰雹和“7·22”岷縣-漳縣地震后，朱家溝泥石流危險性大幅增加，且緊鄰縣城風(fēng)險較大，定西國土局為有效防治地質(zhì)災(zāi)害，在模擬區(qū)內(nèi)修建工程措施見表3。

表3 朱家溝泥石流治理工程參數(shù)Table 3 Parameters of debris flow control project in Zhujia gully

組合降雨、防治工程條件，模擬3種工況下高易發(fā)溝道泥石流的危險性，分別為：①無防治工程時“7·18”強降雨（50年一遇）泥石流，②有防治工程時“7·18”強降雨（50年一遇）泥石流以及③有防治工程時百年一遇強降雨泥石流爆發(fā)。以工況①、②作為對照組評判工程措施效用，以工況②、③作為對照組分析極端氣象下泥石流運動趨勢，以極端工況③作為目標(biāo)組進行危險性區(qū)劃。

計算所需參數(shù)設(shè)置如下：1）重度及體積濃度：據(jù)李國營等在朱家溝的工作成果，主溝、郎家溝、松樹溝、朱林溝和哈固溝的泥石流容重γ（kN/m3）分別為［17.46，17.59，17.53，16.16，17.16］，流域平均容重為17.2 kN/m3，平均比重Gs=2.6 g/cm3，經(jīng)計算流域內(nèi)平均孔隙比為e=1.228，體積濃度C V=0.463，對應(yīng)泥砂比Rns=0.70；2）層流阻滯系數(shù)：本文參考侯圣山等［10］研究，用工程地質(zhì)類比法，取層流阻滯系數(shù)K=2 280；3）屈服應(yīng)力及粘滯系數(shù)：根據(jù)1）中參數(shù)，聯(lián)立粘滯系數(shù)η（Pas）和屈服應(yīng)力τy（kpa）計算公式以及王裕宜等［17］統(tǒng)一的泥砂比—體積濃度—流變參數(shù)關(guān)系式（式7），可得［α1，α2，β1，β2］為［0.012 2，0.015 6，17.41，18.28］，以及［η，τy］為［30.146，56.992］。

式中：η為粘滯系數(shù)，τy為屈服應(yīng)力，CV為體積濃度，Rn s為泥砂比，α1、α2、β1、β2為實驗系數(shù)。

（3）曼寧粗糙系數(shù)：利用王裕宜等［17］統(tǒng)一阻力糙率系數(shù)表征公式，結(jié)合實際調(diào)查主溝、郎家溝、松樹溝、朱林溝和哈固溝的泥深h（m）為［1.7，1.3，1.3，1.5，1.2］，可得各溝曼寧系數(shù)nc為［0.29，0.14，0.14，0.22，0.10］；5）集水點及泥石流流量過程曲線：集水點取自各溝的泥石流形成區(qū)底部流通區(qū)頂部，根據(jù)“7·18”降雨-徑流模擬提取的集水點處清水流量線和放大系數(shù)B F=1/（1-CV）=1.862可得點位處泥石流流量過程曲線（圖6）。

圖6 溝谷集水點分布及泥石流流量過程曲線Fig.6 The distribution of catchment points of gullies and flow process curves of debris flow

3.2 模擬結(jié)果

圖7在實際降雨條件下朱家溝泥石流運動過程的模擬結(jié)果。由圖7可知：（1）“7·18”強降雨時（圖7（a）、（b）、（d）、（e）），無工程條件下最大泥速6.03 m/s，最大泥深4.88～8.07 m，位于溝口堆積扇中部以及溝道匯集處；有工程條件下最大泥速4.54 m/s，三道攔擋壩前最大泥深為7.42～8.58 m，郎家溝至主溝排導(dǎo)槽段受后方補給與前方運動受阻影響，淤積較深可達(dá)4.2～5.75 m，溝口堆積區(qū)最大深度為3.23 m，攔擋壩與排導(dǎo)槽起到攔截回淤、減速、導(dǎo)流控界的作用。（2）百年一遇強降雨時（圖7（c）、（f）），水動力條件增強，雨水潤濕溝道，朱家溝內(nèi)泥速提升、泥深增加，泥石流物質(zhì)沖出溝口，堵塞洮河，直接威脅岷縣縣城。

圖7 多工況下泥深、泥速分布Fig.7 Distribution of mud depth and mud velocity under different conditions

經(jīng)驗證（表4），模擬精度高于74.12%，能反映“7·18”降雨條件下區(qū)內(nèi)泥石流運動特征與堆積范圍，且說明百年一遇降雨條件下泥石流預(yù)測較為可信。

表4 模擬結(jié)果驗證表Table 4 Validation table of simulation results

3.3 泥石流危險性評價

實地調(diào)查表明，朱家溝泥石流直進性明顯，破壞形式以沖毀、淤埋農(nóng)田、民居、公路和擠壓河道，泥速、泥深與危險性呈正相關(guān)關(guān)系，參考侯圣山等［10］在岷縣耳陽河流域的研究（表5），在百年一遇強降雨條件和2020年“7·18”強降雨條件下，對朱家溝泥石流進行危險性分級（圖8）。

表5 泥石流危險性分區(qū)標(biāo)準(zhǔn)Table 5 Criteria for hazard zoning of debris flow

泥石流危險性評價結(jié)果表明：“7·18”泥石流總危險區(qū)面積為371 875 m2，高危險區(qū)占危險區(qū)總面積的45.72%，主溝中下游均處于高危險區(qū)，溝口排導(dǎo)槽將高危險性區(qū)限制于槽中，抑制堆積扇成形，此外松樹溝內(nèi)三道攔擋壩處由于起到攔擋回淤作用表現(xiàn)為高危險性；中危險區(qū)占危險區(qū)總面積的30.70%，分布在松樹溝、朱林溝以及主溝上游區(qū)域、溝口排導(dǎo)槽出口周圍；低危險區(qū)占危險區(qū)總面積的23.58%，分布于哈固溝和郎家溝。P=1%概率下，泥石流總危險區(qū)面積為500 625 m2，與“7·18”降雨條件相比，危險區(qū)面積增大34.62%。其中高危險區(qū)占總面積71.62%，增加了188 525 m2，具體表現(xiàn)為主溝溝頭處危險性上升，朱家溝主溝沿岸房屋建筑易受泥石流威脅，以及溝口沖積扇成型，主要威脅洮河左岸溝口處與右岸30 m范圍內(nèi)城區(qū)建筑；中危險區(qū)面積約占18.34%，包括朱林溝、哈固溝和松樹溝，但由于松樹溝內(nèi)攔擋工程淤滿，超過極限狀態(tài)，使其周邊危險性較高；低危險區(qū)面積約占10.04%，郎家溝整體危險性較低，溝口排導(dǎo)設(shè)施仍可持續(xù)倒流。

總體而言，朱家溝泥石流危險性大。其每年雨季向洮河瀉流，泥石流物質(zhì)在溝河交匯處淤積，抬高河床，在極端的氣象條件下（圖8（b））有阻斷洮河，形成堰塞壩的風(fēng)險，進而會迫使上游水位增漲，淹沒岷縣城區(qū)。并且朱家溝溝口下游2.5 km處為龍王臺水電站，溝內(nèi)流體瀉入洮河后，存在部分流體順流向下匯入水電站庫區(qū)，擠占庫容的可能性。

圖8 朱家溝泥石流危險性分區(qū)圖Fig.8 Zhujia gully debris flow risk zoning map

4 結(jié)論與討論

本文在青藏高原東北部山地泥石流災(zāi)害調(diào)查的基礎(chǔ)上，選取具有區(qū)域代表性的泥石流孕災(zāi)小流域——朱家溝為研究區(qū)，用于探索適合青藏高原東北緣小流域泥石流溝的危險性評價方法。所得方法將依托熵權(quán)-變異系數(shù)融合算法的多因子疊加模型和數(shù)值模擬技術(shù)相結(jié)合，前者通過衡量泥石流形成的地形、地質(zhì)條件評價小流域內(nèi)各溝道泥石流發(fā)育程度，判別在降雨誘發(fā)下，區(qū)內(nèi)易于暴發(fā)泥石流并匯入主溝，參與主溝泥石流流體物質(zhì)補給，最終共同致災(zāi)的支溝；后者則結(jié)合前者結(jié)果，引入泥石流觸發(fā)因素，通過模擬極端降雨條件下流域內(nèi)泥石流的爆發(fā)過程和運動特征，進行泥石流危險性區(qū)劃，二者相互補充，充分考慮青藏高原東北緣小流域的泥石流孕災(zāi)條件，用于泥石流危險性評價具有充分的合理性。經(jīng)2020年朱家溝“7·18”泥石流實例驗證，該法具有較高的準(zhǔn)確度，用其進行朱家溝泥石流危險性評價得到以下結(jié)論和認(rèn)識：

（1）朱家溝內(nèi)各支溝泥石流易發(fā)程度主要受溝口單寬清水流量峰值、主溝長度、淺表層巖組等條件控制，區(qū)內(nèi)17條支溝中高易發(fā)溝6條，中易發(fā)溝5條，低易發(fā)溝6條，其中哈固溝、朱林溝、松樹溝和郎家溝等主要支溝均為高易發(fā)程度泥石流溝。

（2）泥石流流速和溝谷淤積數(shù)值模擬結(jié)果表明，在50年一遇降雨強度以下，朱家溝溝口的排導(dǎo)槽以及松樹溝內(nèi)三座攔擋壩對溝內(nèi)泥石流起到了導(dǎo)流控界和攔截減速的作用，但基本已達(dá)到其能力極限。在百年一遇的極端降雨條件下，對溝口居民和房屋、岷縣縣城與下游龍王臺水電站威脅較大，建議在主溝修建攔擋壩等防治工程，并及時清淤。

（3）“7·18”朱家溝泥石流危險區(qū)集中分布于主溝、哈固溝、朱林溝、松樹溝和郎家溝溝道，高危險區(qū)分布于主溝中下游區(qū)域。當(dāng)遭遇百年一遇強降雨條件（P=1%），朱家溝泥石流危險區(qū)面積約500 625 m2，較“7·18”泥石流危險性區(qū)面積增加34%，其中高危險區(qū)分布于朱家溝主溝、溝口處與洮河右岸30 m范圍，對岷縣縣城構(gòu)成重大威脅。評價結(jié)果可為岷縣縣城的地質(zhì)災(zāi)害防治、國土空間規(guī)劃提供支撐。

小流域泥石流危險性區(qū)劃的影響因素眾多，從泥石流高易發(fā)支溝的選取到模擬時所選定的極端降雨條件均直接影響評價結(jié)果。本文雖建立起包含了12個地形地質(zhì)方面致災(zāi)因子的泥石流易發(fā)程度評價體系，但忽略了因子相關(guān)性，未能將高度相關(guān)因子剔除或整合，一定程度上影響了高易發(fā)支溝的選擇，同時，僅重現(xiàn)了百年一遇降雨條件下的泥石流過程，難以全面反映爆發(fā)不同規(guī)模泥石流時流域的危險性分布。因此，為得到更準(zhǔn)確的評價結(jié)果，在下一步的工作中應(yīng)重點提煉泥石流致災(zāi)因子，凝練評價體系，以及開展不同雨強條件下小流域泥石流危險性綜合區(qū)劃。

致謝

本文在撰寫過程中使用了岷縣氣象局提供的朱家溝流域降雨資料以及岷縣岷陽鎮(zhèn)朱家溝泥石流治理工程勘察資料，在此表示誠摯的感謝。