三峽庫(kù)區(qū)白家包滑坡涌浪計(jì)算分析

盧書強(qiáng) 黃波林 梁衛(wèi) 張巷生

摘要：

對(duì)白家包滑坡穩(wěn)定性及其涌浪災(zāi)害進(jìn)行研究對(duì)三峽庫(kù)區(qū)人民生命財(cái)產(chǎn)安全保護(hù)具有重要意義，將風(fēng)險(xiǎn)量化也可以有效反映滑坡災(zāi)害影響程度。通過(guò)美國(guó)土木工程師協(xié)會(huì)推薦的公式求解滑坡速度，建立數(shù)學(xué)模型求解滑坡涌浪傳播浪高與傳播距離，然后與水波動(dòng)力學(xué)理論改進(jìn)的FAST數(shù)值模擬軟件結(jié)果進(jìn)行對(duì)比分析。分析結(jié)果表明：庫(kù)水位175 m工況下，滑坡速度為2.19 m/s，最大首浪高度3.30 m，對(duì)岸爬高3.35 m，傳播至長(zhǎng)江主航道及其支流最大浪高沒(méi)有超過(guò)1.50 m;庫(kù)水位145 m工況下，滑坡滑速為3.44 m/s，最大首浪高度5.73 m，對(duì)岸爬高10.65 m，傳播至長(zhǎng)江主航道及其支流最大浪高沒(méi)有超過(guò)2.0 m。數(shù)學(xué)模型模擬結(jié)果較FAST數(shù)值模擬結(jié)果偏小。

關(guān) 鍵 詞：

滑坡涌浪; 滑坡速度;最大首浪高度; 改進(jìn)FAST模型; 白家包滑坡

中圖法分類號(hào)： P642

文獻(xiàn)標(biāo)志碼： A

DOI：10.16232/j.cnki.1001-4179.2021.07.017

0 引 言

自三峽水庫(kù)2003年蓄水以來(lái)，白家包滑坡頗受國(guó)內(nèi)外專家學(xué)者的關(guān)注，其形成機(jī)理[1-3]及階躍型變形特征[4-5]是專家學(xué)者重點(diǎn)研究的方向。其中，關(guān)于涌浪研究，以滑坡速度、最大浪高、沿程浪高、對(duì)岸爬高的研究較為普遍，而滑坡機(jī)理研究較少。盧書強(qiáng)等[1]根據(jù)現(xiàn)場(chǎng)地質(zhì)調(diào)查及勘查資料，結(jié)合10多年的監(jiān)測(cè)數(shù)據(jù)，分析得出滑坡體滲透系數(shù)較低、周期性蓄水導(dǎo)致滑坡穩(wěn)定性發(fā)生周期性變化、庫(kù)水位升降速率是影響滑坡穩(wěn)定性的重要因素。馬倩等[6]采用物理模型試驗(yàn)對(duì)滑坡涌浪展開研究，認(rèn)為影響涌浪爬坡高度的因素主要有6種，并以新灘滑坡為例推出了滑坡涌浪爬高計(jì)算公式，通過(guò)與傳統(tǒng)經(jīng)驗(yàn)公式對(duì)比，驗(yàn)證了推斷公式的可靠性。魯芃等[7]利用Geo-Studio和FLAC3D對(duì)白家包滑坡展開不同庫(kù)水位日降幅對(duì)滑坡穩(wěn)定性影響的數(shù)值模擬分析，認(rèn)為該滑坡目前處于欠穩(wěn)定和不穩(wěn)定狀態(tài)，增大庫(kù)水位日降幅至1.2 m/d，滑坡將處于不穩(wěn)定狀態(tài)。黃波林等[8-9]基于地質(zhì)災(zāi)害涌浪計(jì)算公式和局部水頭損失理論，建立了涌浪公式計(jì)算體系，并對(duì)龔家方崩滑體涌浪進(jìn)行了分析，計(jì)算出滑體下滑速度、最大浪高及最大爬高?；糁緷萚10-11]基于水波動(dòng)力理論，利用FAST模擬軟件對(duì)黑石板滑坡進(jìn)行涌浪計(jì)算分析及預(yù)測(cè)，得知滑坡入江速度為5.56 m/s，最大涌浪高度37.2 m，對(duì)岸最大爬高36.0 m。趙永波等[12]通過(guò)改進(jìn)的FAST模擬軟件對(duì)千將坪滑坡進(jìn)行涌浪分析，得到最大浪高38.8 m，最大爬高36.7 m。黃波林等[13-14]引入水波動(dòng)力學(xué)模型，結(jié)合GIS技術(shù)開發(fā)出FAST軟件，以茅草坡滑坡為例，計(jì)算出滑坡最大涌浪高度為25.0 m，最大涌浪爬高為12.5 m。王偉等[15]通過(guò)分析滑坡產(chǎn)生的脈沖波，提出了耦合的DDA-SPH方法來(lái)解決固液相互作用問(wèn)題，評(píng)價(jià)了耦合方法在滑坡運(yùn)動(dòng)與波廓建模中的準(zhǔn)確性，證明了耦合方法能夠準(zhǔn)確地捕捉滑坡運(yùn)動(dòng)與波廓。本次研究假設(shè)白家包滑坡成災(zāi)，在此基礎(chǔ)上進(jìn)行涌浪風(fēng)險(xiǎn)等級(jí)評(píng)估及風(fēng)險(xiǎn)區(qū)域劃分。

1 白家包滑坡概況

白家包滑坡位于秭歸縣歸州鎮(zhèn)向家店村，距香溪河入江口2.5 km，衛(wèi)星圖見圖1。

大地坐標(biāo)：經(jīng)度110°45′33.4″，緯度 30°58′59.9″，白家包滑坡展布于香溪右岸，前緣直抵香溪河?；录舫隹谖挥诟叱?25～135 m之間;滑坡后緣以基巖為界，高程265 m;左側(cè)以山脊下部基巖為界，右側(cè)以山梁為界。前緣寬500 m，后緣寬300 m，均寬400 m，縱長(zhǎng)約550 m，滑坡面積22萬(wàn) m2?；缕旅嫫露?0°～15°，滑體前緣臨江段坡度為20°，中部平緩，坡度為10°～12°，滑坡平均坡度約15°。滑坡平面形態(tài)呈短舌狀，深層滑體前緣厚20～30 m，中部厚47 m，后緣厚10～40 m，滑體平均厚度為45 m，滑體體積990萬(wàn)m3。淺層滑體前緣厚10～20 m，中部厚35 m，后緣厚10～40 m，滑體平均厚度30 m，滑體體積660萬(wàn)m3。滑坡全貌如圖2所示。

2 滑坡滑速計(jì)算

滑坡滑速計(jì)算采用美國(guó)土木工程師協(xié)會(huì)推薦的公式（1）計(jì)算，計(jì)算剖面如圖3所示。

滑坡速度計(jì)算公式：

式中：α為滑面傾角;W 為滑體單寬重量;f、c為滑動(dòng)時(shí)滑面抗剪強(qiáng)度參數(shù);H為滑體重心距離水面的位置;l為滑塊與滑面接觸面長(zhǎng)（沿滑動(dòng)方向）。

3 白家包滑坡涌浪數(shù)學(xué)模型

3.1 數(shù)學(xué)模型建立

白家包滑坡為一典型土質(zhì)滑坡，滑坡坡面平均坡度約15°，采用緩傾角分別建立滑坡涌浪最大首浪高度、傳播浪高數(shù)學(xué)計(jì)算模型。

3.1.1 最大首浪高度計(jì)算

緩傾角滑坡涌浪最大首浪高度數(shù)學(xué)計(jì)算模型為

式中：Hmax為最大首浪高度;g為重力加速度;l為滑坡體入水長(zhǎng)度，本文取550 m;w為滑坡體入水寬度，本文取400 m;t為滑坡體入水厚度，取30 m;v為滑坡體入水最大速度，本文取3.44 m/s（庫(kù)水位145 m）和2.19 m/s（庫(kù)水位175 m）;h為水深，分別取65 m（庫(kù)水位145 m）和98 m（庫(kù)水位175 m）;α為滑動(dòng)面傾角，取15°;b為滑坡入水?dāng)嗝婧拥缹挾?，分別為440 m（庫(kù)水位145 m）和650 m（庫(kù)水位175 m）。

3.1.2 涌浪傳播浪高計(jì)算

緩傾角滑坡涌浪傳播浪高計(jì)算數(shù)學(xué)模型推導(dǎo)公式為

式中：Hp為涌浪沿滑坡橫斷面?zhèn)鞑ブ聊程幍膫鞑ダ烁?，m;x為橫斷面某處至滑坡點(diǎn)的距離，m。

3.1.3 白家包滑坡涌浪爬高計(jì)算

滑坡涌浪緩傾角爬高數(shù)學(xué)模型公式為

式中：HR為正對(duì)岸涌浪爬高，m;α0為對(duì)岸爬坡坡角，（°）。

3.2 數(shù)學(xué)模型計(jì)算結(jié)果

將相關(guān)參數(shù)代入公式（2），可以得到在145 m和175 m庫(kù)水位工況下，滑坡首浪最大高度分別為5.73 m和3.30 m，代入式（3）可以得到傳播浪高與傳播距離關(guān)系式：

可以看出：當(dāng)庫(kù)水位為145 m時(shí)，距滑坡點(diǎn)440 m的對(duì)岸涌浪爬高為10.65 m，距離滑坡2 km長(zhǎng)度范圍內(nèi)涌浪高度衰減約1.52 m;當(dāng)水位為175 m時(shí)，距滑坡點(diǎn)650 m的對(duì)岸涌浪爬高3.35 m，距離滑坡2 km長(zhǎng)度范圍內(nèi)滑坡涌浪高度衰減約0.58 m。

4 白家包滑坡涌浪數(shù)值模型驗(yàn)證

本文采用FAST中的水庫(kù)陸地滑坡涌浪源模型進(jìn)行計(jì)算分析。白家包滑坡涌浪計(jì)算區(qū)域的地形數(shù)據(jù)通過(guò)NASA開源30 m矢量化的地形圖得到，同時(shí)參考了一些文獻(xiàn)中所述水深數(shù)據(jù)，綜合以上地形數(shù)據(jù)形成計(jì)算地形。計(jì)算區(qū)域東西長(zhǎng)約18 720 m，南北長(zhǎng)約30 680 m。根據(jù)庫(kù)水位及輸入?yún)?shù)的不同，選用適合的柵格進(jìn)行模擬計(jì)算。采用65 m×65 m的柵格將計(jì)算域劃分為472行、288列。計(jì)算區(qū)域是從香溪河鄭家坪至河口（見圖4）。FAST軟件應(yīng)用了美國(guó)特拉華大學(xué)開發(fā)的四階包辛奈斯克模型開源程序FUNWAVE作為傳播及爬高計(jì)算模型。FUNWAVE是基于Wei的Boussinesq精確解的帶散射的非線性波浪模型，是完全非線性的，能夠模擬各種波而不僅限于長(zhǎng)波。FUNWAVE的主要功能是產(chǎn)生波浪源和邊界條件，模擬波浪傳播以及波浪傳播過(guò)程中波浪破碎和爬高。

計(jì)算每時(shí)步代表0.558 s，設(shè)置計(jì)算3 000時(shí)步，約1 674 s（約28 min）。按照波速24 m/s計(jì)算，涌浪上下游均可以傳播約40.1 km。傳播距離超過(guò)計(jì)算區(qū)域;設(shè)置的時(shí)步數(shù)可以滿足河道涌浪的分析。

利用FAST的水庫(kù)陸地滑坡涌浪源進(jìn)行涌浪源計(jì)算，計(jì)算工況有2種：庫(kù)水位175 m工況和145 m工況，156 m水位由于介于175～145 m之間，其涌浪大小也將介于其中。

4.1 庫(kù)水位175 m工況下涌浪災(zāi)害分析

白家包滑坡涌浪數(shù)值模擬的結(jié)果如圖5所示，滑體入水后產(chǎn)生涌浪，最大涌浪高度約為4.20 m。當(dāng)T=77 s時(shí)，形成最大爬高約4.50 m。此后，波的傳播方向由最開始的滑動(dòng)方向?yàn)橹鬓D(zhuǎn)化為沿河道方向?yàn)橹?，加大了河道的縱向流動(dòng)性。

由于涌浪傳播區(qū)河道蜿蜒延伸，涌浪在河道中傳播時(shí)推進(jìn)方向多變，傳播浪衰減強(qiáng)烈。當(dāng)T為135 s時(shí)，最大涌浪傳播至香溪河河口馬槽嶺，最大浪高約1.40 m，最大爬高約2.00 m。T為180 s時(shí)，最大涌浪傳播至長(zhǎng)江對(duì)岸郭家壩鎮(zhèn)處，最大浪高約0.50 m。

當(dāng)涌浪朝香溪河上游傳播，T為149 s時(shí)，最大涌浪傳播至喬家壩處，最大浪高約0.95 m。T為182 s時(shí)，最大涌浪傳播至萬(wàn)古寺處，最大浪高約0.70 m。從萬(wàn)古寺往上游，河道變窄，水體有所變淺。當(dāng)T為279 s時(shí)，最大涌浪傳播至鹽關(guān)，最大浪高約0.60 m。T為463 s時(shí)，最大涌浪至游家河時(shí)，最大浪高約為0.20 m。

當(dāng)涌浪在長(zhǎng)江干流上傳播，T為667 s時(shí)，最大涌浪傳播至卜莊河溝頭，最大浪高低于0.10 m;當(dāng)T為488 s時(shí)，最大涌浪傳播至吒溪河河口，最大浪高低于0.10 m;當(dāng)T為385 s時(shí)，最大涌浪傳播至屈原鎮(zhèn)場(chǎng)址，最大浪高低于0.10 m;當(dāng)T為573 s時(shí)，最大涌浪傳播至九畹溪口，最大浪高低于0.10 m;當(dāng)T為485 s時(shí)，最大涌浪傳播至泄灘，最大浪高約0.1 m。當(dāng)T為478 s時(shí)，最大涌浪傳播至蘇溪溝溝源頭時(shí)，最大浪高約0.10 m。

可以看出，滑坡涌浪的主要影響區(qū)在香溪河內(nèi)，這一點(diǎn)與其他支流內(nèi)發(fā)生的滑坡涌浪類似，例如千將坪滑坡。當(dāng)蓄水至175 m后，水體變淺引起的涌浪雍高效應(yīng)消失，涌浪消減迅速。

T為120 s后，河道區(qū)域內(nèi)皆為低于1 m的波浪傳播，水面基本進(jìn)入小幅震蕩階段。此時(shí)船舶行駛形成的浪仍可能與涌浪疊加形成較大的波浪，特別是小型船只通行仍具有一定風(fēng)險(xiǎn)性。

急劇衰減區(qū)長(zhǎng)約2 km，浪高從3.10 m下降至0.80 m，這一急劇衰減區(qū)以滑坡入水段為主，是涌浪危害航道的重點(diǎn)區(qū)域（見圖6）。由于河流蜿蜒，衰減機(jī)制不一，平緩衰減區(qū)僅近似滿足緩斜線形式下降規(guī)律，平均100 m內(nèi)涌浪下降高度為0.01～0.02 m不等，上下游的下降率相近。在河道上游，水變淺且河道變狹窄，在上游區(qū)域涌浪又有局部少量雍高現(xiàn)象。水波的折射、反射和疊加作用使得沿程河道中的波高并非呈簡(jiǎn)單單一下降趨勢(shì)，而是一個(gè)復(fù)雜的波變化衰減過(guò)程。支流內(nèi)傳播區(qū)河流兩岸的浪高也大多不一致。

圖7為支流涌浪雍高的情況和河流兩岸浪高不對(duì)稱情況。從圖7可見，河流的中線-水最深處的涌浪一般較低，兩側(cè)淺水區(qū)涌浪雍高。支流內(nèi)一般溝頭涌浪小，而入口處涌浪大一些。在河道剖面方向上，除急劇衰減區(qū)外，涌浪傳播過(guò)程中深水區(qū)的浪高明顯低于淺水區(qū)的浪高，說(shuō)明淺水區(qū)岸坡加劇了波浪壅高和爬高。在溝谷內(nèi)和地形急劇變窄區(qū)域，涌浪高度明顯升高，出現(xiàn)放大效應(yīng)。在河道由寬突然變窄的峽口區(qū)域，如在支流上游峽谷區(qū)和凹槽地形中，涌浪曲線顯現(xiàn)出雍高爬升的現(xiàn)象。而在河道快速變寬緩的區(qū)域，涌浪則進(jìn)一步加速衰減，例如下游河道和長(zhǎng)江河道。

4.2 庫(kù)水位145 m工況下涌浪災(zāi)害分析

由圖8可見：庫(kù)水位145 m時(shí)白家包滑坡形成的涌浪位于河道中心線附近，涌浪高度約6.10 m，抵達(dá)對(duì)岸時(shí)（T=72 s），形成最大爬高約14.10 m。此后，波的傳播方向由最開始的滑動(dòng)方向?yàn)橹鬓D(zhuǎn)化為沿河道方向?yàn)橹?，加大了河道的縱向流動(dòng)性。

由于涌浪傳播區(qū)河道蜿蜒延伸，涌浪在河道中傳播時(shí)推進(jìn)方向多變，傳播浪衰減強(qiáng)烈。當(dāng)T為162 s時(shí)，最大涌浪傳播至香溪河河口馬槽嶺，最大浪高約2.20 m，最大爬高約2.00 m。

T為210 s時(shí)，最大涌浪傳播至長(zhǎng)江對(duì)岸郭家壩鎮(zhèn)處，最大浪高約1.20 m。當(dāng)涌浪朝香溪河上游傳播，T為155 s時(shí)，最大涌浪傳播至喬家壩處，最大浪高約2.20 m;T為296 s時(shí)，最大涌浪傳播至萬(wàn)古寺處，最大浪高約1.20 m。從萬(wàn)古寺往上游，河道變窄，水開始明顯變淺，涌浪又有所雍高。T為522 s時(shí)，最大涌浪傳播至鹽關(guān)，最大浪高約1.40 m;T為659 s時(shí)，最大涌浪至游家河時(shí)，雍高至4.70 m。游家河上游河道水位低于145 m，只有非常淺的自然河道，涌浪傳播基本消失。但是當(dāng)T為598 s時(shí)，最大涌浪傳播至蘇溪溝溝源頭時(shí)，最大浪高約達(dá)1.10 m。

當(dāng)涌浪在長(zhǎng)江干流上傳播，T為810 s時(shí)，最大涌浪傳播至卜莊河溝頭，最大浪高約0.60 m;T為587 s時(shí)，最大涌浪傳播至吒溪河河口，最大浪高約0.30 m;T為433 s時(shí)，最大涌浪傳播至屈原鎮(zhèn)場(chǎng)址，最大浪高約0.20 m;T為642 s時(shí)，最大涌浪傳播至九畹溪口，最大浪高約0.10 m;T為547 s時(shí)，最大涌浪傳播至泄灘，最大浪高約0.30 m?？梢钥闯觯夯掠坷说闹饕绊憛^(qū)在香溪河內(nèi)，這一點(diǎn)與其他支流內(nèi)發(fā)生的涌浪類似，例如千將坪滑坡。在香溪河上游，由于河道變窄，水體變淺。涌浪在向上游傳播過(guò)程中，出現(xiàn)先衰減后雍高的現(xiàn)象。T為120 s后河道區(qū)域內(nèi)皆為低于1 m的波浪傳播，水面基本進(jìn)入小幅震蕩階段。

通過(guò)涌浪最大浪高圖和剖面圖可知，涌浪可分為急劇衰減區(qū)和平緩衰減區(qū)。由于白家包滑坡處于支流的平直段內(nèi)，在急劇衰減區(qū)內(nèi)上下游基本對(duì)稱。急劇衰減區(qū)長(zhǎng)約2 km，從5.70 m下降至1.70 m，平均100 m內(nèi)涌浪高度下降約0.20 m，這一急劇衰減區(qū)以滑坡入水段為主，是涌浪危害航道的重點(diǎn)區(qū)域（見圖9）。由于河流蜿蜒，衰減機(jī)制不一，平緩衰減區(qū)僅近似滿足緩斜線形式下降規(guī)律，平均100 m內(nèi)涌浪下降高度為0.01～0.20 m不等，上下游下降率相近。

在河道上游，由于水變淺且河道變狹窄，在上游區(qū)域涌浪又有局部少量雍高現(xiàn)象，加之水波的折射、反射和疊加作用，使得沿程河道中的波高并非呈簡(jiǎn)單單一下降趨勢(shì)，而是一個(gè)復(fù)雜的波變化衰減過(guò)程。支流內(nèi)傳播區(qū)河流兩岸的浪高也大多不一致。因此，平緩衰減區(qū)的浪高一般為起伏形下降。該區(qū)域是滑坡涌浪危害的拓展區(qū)域，長(zhǎng)度非常長(zhǎng)。在河道剖面方向上，除急劇衰減區(qū)外，涌浪傳播過(guò)程中深水區(qū)的浪高明顯低于淺水區(qū)的浪高，表明淺水區(qū)岸坡加劇了波浪壅高和爬高。在溝谷內(nèi)和地形急劇變窄區(qū)域，涌浪高度明顯升高，出現(xiàn)放大效應(yīng)。在河道由寬突然變窄的峽口區(qū)域，如在支流上游峽谷區(qū)和凹槽地形中，涌浪曲線顯現(xiàn)出雍高爬升的現(xiàn)象。在河道快速變寬緩的區(qū)域，涌浪則進(jìn)一步加速衰減，例如下游河道和長(zhǎng)江河道。

從圖10可以看出：各支流內(nèi)滑坡涌浪情況復(fù)雜多變且浪高較大，譬如卜莊河內(nèi)涌浪波幅大于0.50 m，蘇溪溝支流內(nèi)的涌浪波幅甚至大于1.00 m。這些支流雖遠(yuǎn)離滑坡體，但仍受涌浪影響。

5 結(jié) 論

（1） 采用美國(guó)土木工程師協(xié)會(huì)推薦法計(jì)算出在庫(kù)水位175 m和145 m時(shí)，白家包滑坡滑動(dòng)速度分別為2.19 m/s和3.44 m/s。

（2） 庫(kù)水位175 m工況下，數(shù)學(xué)模型計(jì)算得到的滑坡涌浪首浪高度為3.30 m，而改進(jìn)FAST涌浪模擬軟件模擬的結(jié)果為4.20 m，數(shù)值模擬結(jié)果略大于數(shù)學(xué)模型結(jié)果。

庫(kù)水位145 m工況下，數(shù)學(xué)模型計(jì)算得到的滑坡涌浪首浪高度為5.73 m，改進(jìn)FAST涌浪模擬軟件模擬結(jié)果為6.10 m，數(shù)值模擬結(jié)果同樣略大于數(shù)學(xué)模型結(jié)果。

（3） 庫(kù)水位175 m工況下，數(shù)學(xué)模型計(jì)算得到涌浪對(duì)岸爬高約3.35 m，改進(jìn)FAST涌浪模擬軟件得出的滑坡涌浪對(duì)岸爬高為4.20 m，數(shù)值模擬結(jié)果偏大。

庫(kù)水位145 m工況下，數(shù)學(xué)模型計(jì)算得到涌浪對(duì)岸爬高約10.65 m，改進(jìn)FAST涌浪模擬軟件得出的滑坡涌浪對(duì)岸爬高為13.00 m，數(shù)值模擬結(jié)果偏大。

（4） 兩種工況下，F(xiàn)AST數(shù)值模擬涌浪傳播衰減幅度略大于涌浪傳播數(shù)學(xué)模型衰減幅度。

參考文獻(xiàn)：

[1] 盧書強(qiáng)，張國(guó)棟，易慶林，等.三峽庫(kù)區(qū)白家包階躍型滑坡動(dòng)態(tài)變形特征與機(jī)理[J].南水北調(diào)與水利科技，2016，14（3）：144-149.

[2] 劉伊凡，盧書強(qiáng)，管琪.模糊數(shù)學(xué)理論和數(shù)值分析在滑坡評(píng)價(jià)中的應(yīng)用[J].中國(guó)水土保持，2019（5）：46-49，69.

[3] YAO W M，LI C G，ZUO Q J，et al.Spatiotemporal deformation characteristics and triggering factors of Baijiabao landslide in Three Gorges Reservoir region，China[J].Geomorphology，2019，343：1712-1723.

[4] ZHAO N G，HU B，YI Q L，et al.The coupling effect of rainfall and reservoir water level decline on the Baijiabao landslide in the Three Gorges Reservoir Area，China[J].Geofluids，2017，132：71-82.

[5] XU S L，NIU R Q.Displacement prediction of Baijiabao landslide based on empirical mode decomposition and long short-term memory neural network in Three Gorges area，China[J].Computers and Geosciences，2018，111：312-316.

[6] 馬倩，薛宏程，任昱，等.水庫(kù)滑坡涌浪爬坡特征試驗(yàn)研究[J].人民長(zhǎng)江，2020，51（2）：195-199，218.

[7] 魯芃，王世梅，郭振，等.庫(kù)水位日降幅對(duì)白家包滑坡穩(wěn)定性影響分析[J].人民長(zhǎng)江，2017，48（9）：50-53，69.

[8] 黃波林，陳小婷，殷躍平，等.滑坡崩塌涌浪計(jì)算方法研究[J].工程地質(zhì)學(xué)報(bào)，2012，20（6）：909-915.

[9] 黃波林.水庫(kù)滑坡涌浪災(zāi)害水波動(dòng)力學(xué)分析方法研究[D].武漢：中國(guó)地質(zhì)大學(xué)，2014.

[10] 霍志濤，黃波林，張全，等.三峽庫(kù)區(qū)黑石板滑坡涌浪分析[J].水利水電技術(shù)，2020，51（1）：115-122.

[11] 霍志濤，黃波林，王魯琦，等.三峽庫(kù)區(qū)三門洞滑坡潛在涌浪風(fēng)險(xiǎn)研究[J].災(zāi)害學(xué)，2019，34（4）：107-112.

[12] 趙永波，王健，黃波林，等.基于改進(jìn)的FAST模型的千將坪滑坡涌浪研究[J].華南地質(zhì)與礦產(chǎn)，2017，33（4）：411-418.

[13] 黃波林，殷躍平，王世昌，等.GIS技術(shù)支持下的滑坡涌浪災(zāi)害分析研究[J].巖石力學(xué)與工程學(xué)報(bào)，2013，32（增2）：3844-3851.

[14] 黃波林，殷躍平.水庫(kù)區(qū)滑坡涌浪風(fēng)險(xiǎn)評(píng)估技術(shù)研究[J].巖石力學(xué)與工程學(xué)報(bào)，2018，37（3）：621-629.

[15] WANG W，CHEN G Q，ZHANG H，et al.Analysis of landslide-generated impulsive waves using a coupled DDA-SPH method[J].Engineering Analysis with Boundary Elements，2016，64：71-79.

（編輯：胡旭東）

引用本文：

盧書強(qiáng)，黃波林，梁衛(wèi)，等.

三峽庫(kù)區(qū)白家包滑坡涌浪計(jì)算分析

[J].人民長(zhǎng)江，2021，52（7）：101-107.

Calculation and analysis of surge wave caused by Baijiabao landslide

in Three Gorges Reservoir area

LU Shuqiang1，2，3，HUANG Bolin1，2，3，LIANG Wei1，2，3，ZHANG Xiangsheng1，2，3

（1.National Field Observation & Research Station of Landslide in Three Gorges Reservoir Area of Yangtze River，China Three Gorges University，Yichang 443002，China; 2.Engineering Technology Research Center for Geological Hazard Prevention and Control of Hubei Province，Yichang 443002，China; 3.Key Laboratory of Disaster Prevention & Mitigation of Hubei Province，China Three Gorges University，Yichang 443002，China）

Abstract：

The study on the stability of Baijiabao landslide and its surge disaster is of great significance to the safety protection of people′s lives and property in the Three Gorges Reservoir area.At the same time，the risk quantification can effectively reflect the impact of landslide disasters.Firstly，the landslide velocity was solved by the formula recommended by the American Civil Engineers Association，and then the mathematical model was established to solve the landslide surge height and propagation distance.The calculated results were compared with the simulation results from improved FAST numerical simulation software based on the hydrodynamic theory.The analysis results show that when the reservoir water level is 175 m，the landslide velocity is 2.19 m/s，the maximum first wave height is 3.30 m，the climbing height to opposite shore is 3.35 m，and the maximum wave height propagating to the main channel of the Changjiang River and its tributaries does not exceed 1.5 m.When the reservoir water level is 145 m，the landslide speed is 3.44 m/s，the maximum first wave height is 5.73 m，the climbing height to opposite shore is 10.65 m，and the maximum wave height propagating to the main channel of the Changjiang River and its tributaries does not exceed 2.0 m.The mathematical model results are smaller than FAST numerical simulation results.

Key words：

landslide surge;landslide velocity;maximum first wave height;improved FAST model;Baijiabao landslide