我國是世界上地質(zhì)災(zāi)害最為嚴(yán)重的國家之一，其中鏈生型災(zāi)害造成的損失超過50%?；露陆且活惖湫玩溕鸀?zāi)害，在我國西部高山峽谷地區(qū)的歷史強(qiáng)震中頻繁發(fā)生，如唐家山堰塞壩（2008汶川地震）、紅石巖堰塞壩（2014魯?shù)榈卣穑?，?yán)重影響水電、交通等重大工程安全，強(qiáng)震區(qū)滑坡堵江災(zāi)害風(fēng)險防控技術(shù)研究需求迫切。然而，地震—滑坡—堵江災(zāi)害涵蓋強(qiáng)震下巖土失穩(wěn)導(dǎo)致滑坡啟動、滑體沿斜坡運(yùn)動、滑體入江堆積成壩多個物理環(huán)節(jié)，致災(zāi)因子眾多，鏈生過程復(fù)雜，發(fā)展高效、準(zhǔn)確的災(zāi)害鏈推演技術(shù)對于風(fēng)險防控至關(guān)重要。

　　針對這一問題，中國科學(xué)院武漢巖土力學(xué)研究所巖土工程抗震安全組研究人員通過深入分析地震—滑坡—堵江災(zāi)害鏈的能量演化特征，建立了各個物理環(huán)節(jié)的動量守恒方程，充分利用物質(zhì)點法（MPM）兼具拉格朗日和歐拉算法的優(yōu)勢，發(fā)展了巖土體與水體相互作用的兩相雙質(zhì)點MPM。該方法在一套背景網(wǎng)格下，分別采用一套固相質(zhì)點（如M-C模型）表征巖土體、一套液相質(zhì)點（如人工狀態(tài)方程）描述水體；固、液相質(zhì)點同時與背景網(wǎng)格的計算節(jié)點相互映射，映射過程中兩套質(zhì)點之間發(fā)生動量交換，從而實現(xiàn)水力強(qiáng)耦合模擬。首先，通過模擬美國Lituya海灣滑坡巖土體與水體的相互作用過程，兩相雙質(zhì)點MPM的有效性得到了驗證；在此基礎(chǔ)上，推演了2011年湖北十堰市二荒村滑坡堵江災(zāi)害鏈生全過程，通過分析速度場與動能，揭示了高速滑動、涌浪制動、堆積成壩三個環(huán)節(jié)呈現(xiàn)的滑體運(yùn)移模式、滑體—水體作用機(jī)制、堰塞壩堵江特征；最后，建立了強(qiáng)震區(qū)U型河谷滑坡涌浪概化模型，通過定義系統(tǒng)能量轉(zhuǎn)換比，從能量角度系統(tǒng)揭示了災(zāi)害鏈生過程中滑坡與水體的動力相互作用機(jī)制。研究成果為強(qiáng)震區(qū)滑坡堵江災(zāi)害鏈推演與風(fēng)險評估提供了理論方法與技術(shù)支持。

　　研究成果得到了中國科學(xué)院青年創(chuàng)新促進(jìn)會（2021325）、面上項目（51779250）、中國科學(xué)院國際合作局國際伙伴計劃項目（131551KYSB20180042）的資助，相關(guān)研究發(fā)表在Journal of Rock Mechanics and Geotechnical Engineering、Engineering Analysis with Boundary Elements、工程科學(xué)與技術(shù)（原四川大學(xué)學(xué)報（工程科學(xué)版））等學(xué)術(shù)期刊。

　　論文鏈接：

　　https://www.sciencedirect.com/science/article/pii/S1674775522000762

　　https://www.sciencedirect.com/science/article/abs/pii/S0955799721002770

　　https://jsuese.scu.edu.cn/html/2022/3/202101133.html

圖1 強(qiáng)震區(qū)滑坡堵江災(zāi)害鏈生過程的物理分析模型

圖2 兩相雙質(zhì)點MPM的求解流程示意圖

圖3 Lituya滑坡誘發(fā)涌浪過程數(shù)值與實驗結(jié)果對比

圖4 湖北二荒村滑坡堵江災(zāi)害鏈全過程推演

圖5 災(zāi)害鏈生過程的能量演化示意圖