王雷 解明禮 黃會寶 柯虎 高強(qiáng)

摘要：針對不同規(guī)?；氯胨a(chǎn)生的涌浪開展三維數(shù)值分析，分析不同規(guī)?；氯胨T發(fā)涌浪災(zāi)害特征，如涌浪高度、涌浪速度、對岸爬高等，探討不同規(guī)?；氯胨l(fā)涌浪對大壩的影響。利用FLOW-3D數(shù)值模擬方法對滑坡失穩(wěn)過程、涌浪形成及傳播、涌浪爬升、涌浪回流的全過程進(jìn)行模擬分析。結(jié)果顯示：310萬m³滑坡入水產(chǎn)生涌浪在對岸最大爬高為54.5 m，壩前涌浪高度為6.69 m，涌浪在右岸壩肩處有小范圍漫壩；80萬m³滑坡入水產(chǎn)生涌浪在對岸最大爬高為26.00 m，壩前涌浪高度為5.38 m，涌浪對大壩安全無影響。結(jié)果表明：310萬m³滑坡入水誘發(fā)涌浪與80萬m³滑坡入水誘發(fā)涌浪相比，致災(zāi)性較強(qiáng)。

關(guān)鍵詞：不同規(guī)模滑坡；涌浪災(zāi)害；FLOW-3D數(shù)值模擬；致災(zāi)性

中圖分類號：P642.22;TV139.23文獻(xiàn)標(biāo)識碼：B文章編號：1001-9235（2024）02-0018-11

Difference Analysis of Wave Disaster Characteristics Induced by Landslides of Different Water Entry Scales

WANG Lei¹，XIE Mingli^1*，HUANG Huibao²，KE Hu²，GAO Qiang²

（1.State Key Laboratory of Geohazard Prevention and Geoenvironmental Protection，Chengdu University of Technology，Chengdu 610059，China;2.CHN Energy Dadu River Hydropower Development Co.，Ltd.，Chengdu 610041，China）

Abstract：This paper conducts a three-dimensional numerical analysis on the surges generated by landslides of different water entry scales，and analyzes the characteristics of surge disasters induced by landslides of different water entry scales， such as surge height， surge speed，and bank climbing height.Meanwhile，the impact of surges caused by landslides of different water entry scales on the dam is explored.The FLOW-3D numerical simulation method is employed to simulate and analyze the entire process of landslide instability，surge formation and propagation，surge climbing，and surge backflow.The results show that the maximum climbing height of the surge generated by the 3.1 million m³landslide of water entry is 54.5 m on the opposite bank，and the surge height in front of the dam is 6.69 m.The surge has a small area of overflow at the right bank dam shoulder.The surge generated by the 0.8 million m³landslide of water entry has a maximum climbing height of 26.00 m on the opposite bank，and the surge height in front of the dam is 5.38 m，without influence exerted by the surge on the dam safety.The results indicate that the induced surge caused by 3.1×10⁶m³landslide of water entry is more catastrophic than that brought by 0.8×10⁶m³landslide of water entry.

Keywords：landslides of different scales;surge disaster;FLOW-3D numerical simulation;catastrophability

近幾年來，由于能源需求迅速增長，中國水電工程的建設(shè)數(shù)量明顯增加，水力發(fā)電已成為中國的主要電力來源^［1］。大渡河流域蘊(yùn)藏著豐富的潛在水資源，但該區(qū)域經(jīng)常面臨滑坡的風(fēng)險(xiǎn)^［2］。

全世界報(bào)道的300多座大壩倒塌事件中，約有10%是由山谷斜坡上形成的山體滑坡引起的^［3］，85%的滑坡發(fā)生在大壩蓄水后2 a內(nèi)^［4］。例如，1963年Vajont滑坡災(zāi)害就是由于蓄水和相關(guān)的地下靜水壓力突然變化造成的。雖然大壩沒有倒塌，但250 m高的沖擊波沖過大壩，導(dǎo)致數(shù)個村鎮(zhèn)被毀，并造成近3 000人死亡^［5］。水庫滑坡失穩(wěn)可能會引起巨大的涌浪，涌浪可能會危及到人民的生命財(cái)產(chǎn)安全，還可能對大壩和水電站的穩(wěn)定運(yùn)行帶來威脅，因此對滑坡涌浪開展研究十分關(guān)鍵。

楊渠鋒等^［6］研究了三峽庫區(qū)陡巖滑坡的初始涌浪高度沿河道傳播的衰減規(guī)律，提出波高衰減方程；曹婷等^［7］建立三維地質(zhì)模型，研究不同滑坡體的自身特性與涌浪爬高的關(guān)系，并提出涌浪爬高計(jì)算公式；肖莉麗等^［8］建立 三峽庫區(qū)河道模型，利用方差法和回歸法推導(dǎo)出水舌高度與長度計(jì)算公式；Koo等^［9］運(yùn)用非線性水波動力學(xué)，分析了不同水位工況下涌浪的傳播過程及其特征；彭輝等^［10］使用彎曲河道型水庫物理模型分析不同水深時首浪高度的變化；Li等^［11］采用TsunamiSquares方法模擬了水位動態(tài)變化下滑坡涌浪的傳播特征；謝海清等^［12］運(yùn)用無量綱分析方法和非線性回歸方法對涌浪災(zāi)害的形成過程進(jìn)行模擬；黃筱云等^［13］利用流體動力學(xué)軟件FlOW-3D模擬滑坡引發(fā)的涌浪在V形河道中的形成和傳播過程；鄧成進(jìn)等^［14］利用三維數(shù)值模擬方法對滑坡涌浪對壩體的影響進(jìn)行研究，分析了壩前最大浪高、壩面水壓力等特性；徐文杰等^［15］基于SPH-DEM流固耦合算法模擬了滑坡涌浪的傳播過程，通過模擬得到漫頂時的洪峰流量、涌浪對大壩的最大水壓力。

上述研究都是針對同一規(guī)?；氯胨T發(fā)的涌浪災(zāi)害特征進(jìn)行分析，而鮮有對不同規(guī)模滑坡入水誘發(fā)涌浪災(zāi)害特征差異性的研究。

本文的研究重點(diǎn)如下：①確定瀑布溝庫首右岸拉裂變形體支護(hù)加固前潛在不同規(guī)模滑坡；②模擬分析滑坡-涌浪災(zāi)害時空演化過程；③對比分析不同規(guī)?；氯胨T發(fā)涌浪災(zāi)害特征（涌浪高度、涌浪速度、對岸涌浪爬高）；④對比分析不同規(guī)模滑坡入水引發(fā)涌浪對大壩的影響。

1 拉裂變形體概述

瀑布溝水庫位于岷江水系的大渡河中游，地處四川省西部雅安市漢源縣和甘洛縣交界地帶，經(jīng)度102°23′～102°50′E，緯度29°8′～29°27′N。研究區(qū)距離成都市270 km，距上游漢源縣城約28 km，距石棉縣城約80 km。水庫的水域面積約84 km²，是西南地區(qū)最大的人工湖，因?yàn)槠鋷熘芪挥跐h源縣境內(nèi)，又被當(dāng)?shù)厝朔Q作漢源湖。區(qū)內(nèi)交通便利，以鐵路、公路為主，京昆高速G5、國道G108、省道S306橫穿境內(nèi)，是川西交通次樞紐，縣道和鄉(xiāng)道縱橫交錯，四通八達(dá)。隨著大渡河流域水利水電工程的不斷開發(fā)，導(dǎo)致該地區(qū)的滑坡、泥石流等各種地質(zhì)災(zāi)害頻繁爆發(fā)，對當(dāng)?shù)鼐用裆?cái)產(chǎn)安全構(gòu)成極大威脅。

受河谷下切的影響，谷坡巖體快速卸荷拉裂，在瀑布溝水電站右岸壩軸線上游約798 m處形成了拉裂變形體，拉裂變形體位置及全貌見圖1。其前、后緣高程分別為730、1 187 m，高差約為457 m，兩側(cè)有深10～35 m沖溝切割，邊坡坡度約為40°～50°。

拉裂變形體地層以前震旦系淺變質(zhì)玄武巖、震旦系下統(tǒng)蘇雄組凝灰?guī)r為主，第四系崩坡積物基本上分布于高程980 m以上，結(jié)構(gòu)松散至稍密。拉裂變形體的工程地質(zhì)剖面見圖2。現(xiàn)場地質(zhì)調(diào)查及勘察資料顯示，760～820 m高程巖體陡傾角裂隙發(fā)育，巖體多出現(xiàn)傾倒拉裂；820～970 m高程巖體發(fā)育有結(jié)構(gòu)面EW/N∠25°～35°、N26W/NE∠35°～40°，其已基本成為控制性結(jié)構(gòu)面；970～1 180 m高程巖體卸荷特征明顯，強(qiáng)卸荷段巖體呈碎裂-散體結(jié)構(gòu)，巖塊之間發(fā)生了相互錯動。巖體的風(fēng)化程度相對較低，卸荷和松弛現(xiàn)象比較明顯。巖體強(qiáng)卸荷深一般為30～60 m；弱卸荷深一般為65～90 m。巖體卸荷主要表現(xiàn)為沿中陡傾坡內(nèi)的傾倒拉裂及其與順坡結(jié)構(gòu)面組合形成滑移拉裂。

岸坡卸荷作用強(qiáng)度與地貌密切相關(guān)，谷坡越高越陡，卸荷作用越強(qiáng)。巖體卸荷與巖體聲波縱波速度V_p有較好的相關(guān)關(guān)系，強(qiáng)卸荷帶各種巖體V_p=3.2～4.2 km/s；弱卸荷帶各種巖體V_p＜3.0 km/s。

由于低線公路（大致高程為875 m）施工（采用硐室大爆破方式開挖），同時路塹邊坡嚴(yán)重超挖形成岸坡切腳，導(dǎo)致坡體松弛變形和表層覆蓋層的失穩(wěn)。硐室大爆破造成了坡體淺表部碎裂結(jié)構(gòu)巖體結(jié)構(gòu)面的張開和巖塊間的相互位移，并在拉裂變形體后緣覆蓋層與基巖接觸面間形成從上游往下游逐漸尖滅的拉裂縫。路塹邊坡的切腳引起岸坡的應(yīng)力重分布，在路塹邊坡頂部形成應(yīng)力集中，造成變形破壞。

2 潛在不同規(guī)模滑坡

2.1 淺表層巖土體變形破壞

支護(hù)加固前拉裂變形體的淺表層巖土體穩(wěn)定性問題可能出現(xiàn)在970～1 180 m高程硐7號勘探平硐附近（岸坡上部）處于強(qiáng)風(fēng)化帶的碎裂-散體結(jié)構(gòu)淺表層巖體，該部分巖體的變形破壞主要表現(xiàn)為傾倒拉裂、滑移拉裂或滑移壓致拉裂，一旦該部分巖體發(fā)生滑移，則可能形成滑坡。根據(jù)滑坡周界和淺表層滑動的滑動面，確定淺表層滑坡一次下滑平面面積約為3.63萬m²，滑體平均垂直厚度約為22 m，因此推測淺表層滑坡方量為80萬m³，見圖3a。

2.2 深部巖土體變形破壞

在拉裂變形體支護(hù)加固前，強(qiáng)烈的拉裂變形可能導(dǎo)致深部巖土體穩(wěn)定性問題，主要出現(xiàn)在表層滑塌區(qū)到強(qiáng)卸荷底界范圍內(nèi)的巖體。該區(qū)域的巖體出現(xiàn)了不同類型的變形，包括傾倒拉裂、滑移拉裂以及滑移壓致拉裂等。這一部位巖體一旦發(fā)生滑移，則可能形成滑坡。根據(jù)滑坡周界和深部滑動的滑動面，確定深部滑坡一次下滑平面面積約為8.85萬m²，滑體平均垂直厚度約為35 m，因此推測深部滑坡方量為310萬m³，見圖3b。

3 涌浪分析計(jì)算模型

FLOW-3D是基于有限差分法對控制方程進(jìn)行求解的數(shù)值計(jì)算軟件，在模擬過程中將計(jì)算區(qū)域設(shè)置為三維網(wǎng)格，每個網(wǎng)格在計(jì)算過程中都有多個相關(guān)變量的值在每個網(wǎng)格上體現(xiàn)，其中網(wǎng)格表面包括3個方向的速度u、v、w，網(wǎng)格中心包括壓強(qiáng)p、流體體積分?jǐn)?shù)F、密度ρ、湍流動能q、紊流動能消耗率以及黏度μ等變量，這能夠更真實(shí)地還原滑坡涌浪的實(shí)際運(yùn)動過程。

本文利用顆粒流模型、漂移-通量模型、重力和非慣性框架模型以及黏流和湍流模型相互組合分析滑坡涌浪形成、傳播過程。其中控制方程包括連續(xù)方程和動量方程，見式（1）—（4）：

連續(xù)方程：

動量方程：

式中 u、v、w——x、y、z 3個方向上的速度；A_x、A_y、A_z——流體所流經(jīng)單元格x、y、z 3個方向上的面積比例；G_x、G_y、G_z——流體在網(wǎng)格x、y、z 3個方向上的加速度；f_x、f_y、f_z——流體在網(wǎng)格x、y、z 3個方向上的阻力加速度；V_F——流體體積分?jǐn)?shù)；ρ——流體密度；P——作用在流體上的壓強(qiáng)。

在模擬中采用計(jì)算流體力學(xué)方法（CFD）求解水體流動。對滑坡可采用的計(jì)算簡化模型有多種，例如剛性體、可變形體和顆粒體模型等。按照瀑布溝庫首右岸拉裂變形體可能發(fā)生的失穩(wěn)破壞類型，本文采用顆粒流模型進(jìn)行滑坡運(yùn)動過程模擬能比較真實(shí)的反映其實(shí)際情況。固體顆粒和流體形成混合物的兩相效應(yīng)采用漂移-通量模型近似分析，在漂移通量模型中，可以來定義可能發(fā)生的最大固體分?jǐn)?shù)，一旦達(dá)到極限體積分?jǐn)?shù)，漂移速度即為零，充分利用不同計(jì)算模型的優(yōu)點(diǎn)，對滑坡及涌浪演進(jìn)過程進(jìn)行數(shù)值模擬分析。

由滑坡產(chǎn)生的沖擊波形成的大壩水庫中的水流通常是湍流。水庫中由沖擊波引起的湍流的特征是“疊加的高頻隨機(jī)波動”。因此，本文中計(jì)算模型采用RNGk- ε湍流模型，它可以更好地表述涌浪的復(fù)雜運(yùn)動和能量耗散。

構(gòu)建出的滑坡誘發(fā)涌浪的計(jì)算模型見圖4，每個模型均采用結(jié)構(gòu)化正交網(wǎng)格的非均勻網(wǎng)格方法劃分計(jì)算區(qū)域，并對網(wǎng)格進(jìn)行漸近加密以優(yōu)化計(jì)算區(qū)域，總單元格數(shù)約為9.9×10⁶個。模型底部為無滑移邊界條件的實(shí)心墻，周圍側(cè)壁為自由流出邊界，頂部為對稱邊界。模型計(jì)算的整個過程一直處于重力作用下，重力加速度為9.8 m/s²。

4 滑坡-涌浪災(zāi)害時空演化過程分析

基于前文中的涌浪分析計(jì)算模型，模擬了支護(hù)加固前正常蓄水位工況下（850 m）310萬m³的滑坡入水形成涌浪傳播動態(tài)演進(jìn)過程。由圖5可知，t=6 s左右時，由于滑坡重力的作用，涌浪逐漸形成；t=12 s，隨著滑坡滑落，涌浪高度越來越大，呈圓弧狀向外側(cè)傳播；t=15 s，滑坡體全部進(jìn)入水中，涌浪在對岸開始爬升；t=23 s，滑坡堆積于河道底部，此時對岸涌浪爬升高度最大，然后涌浪開始回流，回流水體與新形成波峰相互碰撞造成了涌浪傳播圈，涌浪繼續(xù)向上下游傳播；t=36 s，涌浪傳播至壩址附近，然后壩前水面呈持續(xù)上升態(tài)勢；t=52 s，在右岸壩肩處有小范圍漫壩，然后壩前涌浪開始回落衰退。

由圖6可知，在滑坡運(yùn)動方向上的涌浪傳播過程大致可分為3個階段。①初始涌浪形成及傳播階段?；率Х€(wěn)后在重力作用下進(jìn)入水中，涌浪高度迅速增加同時向相鄰水域擴(kuò)散。在滑坡高速沖擊水面后，水面向前翻卷，形成了初始浪高。②涌浪爬升階段。由于河谷較窄小，涌浪傳播到對岸的時間為15 s，由于慣性力的推動，涌浪沿坡面爬升，23 s左右涌浪爬高最大，達(dá)到54.5 m。③回流階段。28 s后，滑坡已不在運(yùn)動，對岸爬坡涌浪逐漸回流。

為研究滑坡涌浪在河道方向的傳播規(guī)律，選取見圖7的A、B、C、D、E 5個監(jiān)測點(diǎn)進(jìn)行涌浪高度監(jiān)測。

由圖8可知，拉裂變形體失穩(wěn)破壞進(jìn)入水中后，涌浪波峰隨即擴(kuò)散到A、B、C、D、E各監(jiān)測點(diǎn)的時間各為52、23、25、44、59 s。涌浪到達(dá)A點(diǎn)用時比涌浪到達(dá)幾近等距離的D點(diǎn)用時更長，這證明了在涌浪的傳播過程中，水庫地形條件對涌浪的傳播物理特性有較強(qiáng)的干擾。各監(jiān)測點(diǎn)最大浪高分別為8.7、12.0、6.4、4.3、1.6 m，且發(fā)生在第1個波峰，后續(xù)最大浪高出現(xiàn)衰減。

5 不同規(guī)?；氯胨l(fā)涌浪災(zāi)害特征研究

5.1 沿滑坡運(yùn)動方向不同規(guī)模滑坡入水引發(fā)涌浪高度研究

為了將涌浪高度在滑坡運(yùn)動方向的波動精準(zhǔn)地表達(dá)出來，在滑坡運(yùn)動方向設(shè)置3個監(jiān)測點(diǎn)：滑坡入水處a點(diǎn)、河道中央b點(diǎn)和對岸岸邊c點(diǎn)，見圖9。

由圖10可知，310萬m³滑坡入水產(chǎn)生的初始涌浪高度為8.8 m。最大浪高用時13 s傳播至b點(diǎn)，此時b點(diǎn)的最大涌浪高度為17 m。對岸處c點(diǎn)涌浪高度在18 s時達(dá)到22.1 m；80萬m³滑坡入水產(chǎn)生的初始涌浪高度為6.7 m，最大浪高用時10 s傳播至b點(diǎn)，此時b點(diǎn)的最大涌浪高度為11.3 m。對岸處c點(diǎn)涌浪高度在15 s時達(dá)到12.6 m。

5.2 不同規(guī)模滑坡入水引發(fā)涌浪對岸爬高研究

為研究不同規(guī)?；氯胨笮纬傻挠坷嗽趯Π兜呐栏撸瑢⒂坷嗽诨逻\(yùn)動方向的最大高度隨時間變化過程進(jìn)行記錄，并繪制了對應(yīng)的曲線，見圖11。

不同規(guī)模滑坡入水引發(fā)涌浪沿滑坡運(yùn)動方向的傳播高度變化見圖11，滑坡入水方量為310萬m³時，對岸爬坡浪高度為54.5 m，當(dāng)滑坡入水方量為80萬m³時，涌浪的傳播高度顯著降低，對岸爬坡浪高度為26 m。

5.3 沿河道方向不同規(guī)?；氯胨l(fā)涌浪高度與速度研究

為研究不同規(guī)模滑坡入水后形成的涌浪在河道方向的運(yùn)動特征，以滑坡入水點(diǎn)為界，按上下游分組，將河道中最大浪高隨涌浪傳播距離變化過程以及河道中涌浪速度隨涌浪傳播距離變化過程進(jìn)行記錄，并繪制了對應(yīng)的涌浪速度-距離曲線，見圖12，最大浪高-距離曲線，見圖13。

不同規(guī)模滑坡入水引發(fā)涌浪在河道中的傳播速度見圖12，涌浪在向上游的傳播過程中，不同規(guī)模滑坡入水產(chǎn)生涌浪的速度都在某段距離上呈現(xiàn)出上升態(tài)勢，然后速度變小。

最后階段，涌浪速度與傳播距離成反比。涌浪向下游的傳播過程中，不同規(guī)?；氯胨?，涌浪傳播速度都會出現(xiàn)“下降—上升—下降”的階段，同時速度波動較大。

在滑坡涌浪速度處于上升階段時，310萬m³滑坡入水引發(fā)涌浪傳播速度的增幅大于80萬m³滑坡入水引發(fā)涌浪傳播速度的增幅。在整個傳播過程中，310萬m³的滑坡入水引發(fā)涌浪傳播速度更大，即310萬m³的滑坡入水引發(fā)涌浪傳播至上下游的時間更短。

不同規(guī)?；氯胨l(fā)涌浪在河道中的最大涌浪高度見圖13，310萬m³滑坡入水產(chǎn)生涌浪沿河道方向的最大浪高約為17 m。由于310萬m³滑坡入水產(chǎn)生的動能更大，導(dǎo)致其沿河道方向的最大浪高均大于80萬m³滑坡入水產(chǎn)生涌浪沿河道方向的最大浪高。

5.4 不同規(guī)?；氯胨l(fā)涌浪對大壩的影響分析

大壩在滑坡入水點(diǎn)下游約798 m處，大壩為礫石土心墻堆石壩，壩頂高程為856 m，根據(jù)圖13b可知，310萬m³滑坡體入水產(chǎn)生涌浪在壩前的最大高度為6.69 m，在右岸壩肩處出現(xiàn)漫壩，見圖14。漫壩持續(xù)時間短，過流量有限，所以涌浪對大壩和下游的危害極??；80萬m³滑坡體入水誘發(fā)涌浪在壩前的高度為5.38 m，未超過大壩留有的安全裕度，所以無漫壩風(fēng)險(xiǎn)。

6 數(shù)值模擬結(jié)果合理性分析

6.1 最大滑速計(jì)算

根據(jù)滑坡涌浪研究現(xiàn)狀、基本計(jì)算方法及其優(yōu)缺點(diǎn)、使用普遍性、合理性等，結(jié)合研究區(qū)的工程地質(zhì)環(huán)境、基本特征、形成機(jī)制、巖土特性，采用潘家錚法^［16］來計(jì)算滑坡滑速。

滑坡速度采用條分方法，沿滑面方向和垂直滑面方向建立平衡方程，假定當(dāng)前條塊所受的前后塊體作用力合力沿滑面方向，并認(rèn)為垂直滑面方向的加速度為零，對塊體加速度進(jìn)行求解，然后用基本運(yùn)動學(xué)原理，并考慮水阻力，求出塊體運(yùn)動速度。

滑坡體在滑動過程的某一時刻，可以認(rèn)為所有條塊的水平加速度相等。所以可以解出為：

設(shè)條塊i初速度、末速度分別為v_i1、v_i2，滑動距離為L_i，滑動時間為T_i，進(jìn)一步可得：

式中 W_i——第i塊體的重量，kN；M_i——第i塊體質(zhì)量，kg；c_i——第i塊體滑帶的動內(nèi)聚力，kPa；φ_i——第i塊體滑帶的動摩擦角，（°）；L_i——第i塊體底面的長度，m；α_i——第i條塊滑面傾角，（°）；U_i、U_i+1——第i塊體水下部分所受靜水壓力，kN；U_bi——第i塊體所受浮托力，kN。

根據(jù)前面所述的計(jì)算過程進(jìn)行計(jì)算，得到滑坡在水庫正常蓄水位條件下的最大滑速，310萬m³滑坡最大速度為28 m/s，80萬m³滑坡最大速度為15 m/s。

6.2 涌浪傳播規(guī)律分析

傳統(tǒng)的潘家錚法計(jì)算涌浪傳播時單純考慮滑速大小及庫水深度計(jì)算滑體入水產(chǎn)生的初始涌浪高度顯然是不合理的，影響浪高的主要因素還取決于滑坡規(guī)模、入水角度、河道特征等，殷坤龍等^［17］在潘家錚法的基礎(chǔ)上，參考物理模擬實(shí)驗(yàn)的觀測數(shù)據(jù)，考慮滑坡相對單寬體積lt/bh，滑坡相對寬度w/b，總結(jié)出具有參考意義的涌浪計(jì)算公式，由于試驗(yàn)中滑體較小，且所有滑體全部入水，因此在實(shí)際計(jì)算中需要乘以相關(guān)系數(shù)修正。

初始浪高公式：

式中 相關(guān)系數(shù)為0.76;H_max——最大首浪高度；h——水庫平均水深；l——滑體入水時的平均長度；v——最大滑速；t——滑體入水時平均的厚度；b——河道寬度；w——滑體平均寬度。

涌浪向?qū)Π兜膫鞑ネ嬖谙魅酰?/p>

距滑體不同距離的涌浪高度為：

式中 相對系數(shù)為0.93；H_x——距離滑坡x處的涌浪波高。

80萬m³滑坡的計(jì)算參數(shù)取值見表1，310萬m³滑坡的計(jì)算參數(shù)取值見表2?，F(xiàn)對80萬m³滑坡、310萬m³滑坡滑落引發(fā)的涌浪高度進(jìn)行計(jì)算，得到的初始浪高與對岸涌浪爬高見表3，得到的涌浪沿河傳播規(guī)律見表4、5。

大壩在滑坡下游約798 m處，計(jì)算得到310萬m³滑坡體入水產(chǎn)生涌浪在壩址處的最大高度為6.06 m，80萬m³滑坡體入水產(chǎn)生涌浪在壩址處的最大高度為4.96 m（表6）。

由表6可知，根據(jù)理論公式法得到的涌浪高度結(jié)果和使用FlOW-3D軟件進(jìn)行數(shù)值模擬計(jì)算得到的涌浪高度結(jié)果較相近，所以FlOW-3D數(shù)值模擬方法得到的結(jié)果較可靠。

7 結(jié)論

本文依托于瀑布溝庫首右岸拉裂變形體，模擬分析了滑坡-涌浪災(zāi)害時空演化過程，研究其支護(hù)加固前潛在不同規(guī)?；氯胨T發(fā)涌浪災(zāi)害的差異性。結(jié)論如下。

a）拉裂變形體在支護(hù)加固前若發(fā)生淺表層巖土體變形破壞則可能形成80萬m³的滑坡，若發(fā)生深部巖土體變形破壞則可能形成310萬m³的滑坡。

b）310萬m³滑坡入水引發(fā)的初始涌浪高度為8.8 m、對岸涌浪高度為22.1 m、對岸涌浪爬高為54.5 m和壩址處涌浪高度為6.69 m；80萬m³滑坡入水引發(fā)的初始涌浪高度為6.7 m、對岸涌浪高度為12.6 m、對岸涌浪爬高為26 m和壩址處涌浪高度為5.38 m。

c）310萬m³滑坡入水誘發(fā)涌浪在右岸壩肩處有小范圍漫壩，漫壩持續(xù)時間較短，過流量有限，涌浪不會危及大壩及下游安全；80萬m³滑坡入水誘發(fā)涌浪對大壩安全無影響。310萬m³滑坡入水誘發(fā)涌浪與80萬m³滑坡入水誘發(fā)涌浪相比，易形成漫堤以及產(chǎn)生沿河道傳播的沖擊性涌浪，可能對沿岸基礎(chǔ)設(shè)施造成破壞，威脅庫區(qū)安全。

d）采用FlOW-3D數(shù)值模擬的方法，模擬了310萬m³滑坡-涌浪災(zāi)害時空演化過程，根據(jù)滑坡涌浪運(yùn)動特征將滑坡-涌浪災(zāi)害的演進(jìn)過程劃分為4個階段：滑坡失穩(wěn)過程、涌浪形成及傳播、涌浪爬升、涌浪回流。

e）運(yùn)用理論公式法得到的涌浪高度結(jié)果和使用FlOW-3D軟件進(jìn)行數(shù)值模擬計(jì)算得到的涌浪高度結(jié)果較相近，所以FlOW-3D數(shù)值模擬方法得到的結(jié)果較可靠。

本文僅討論了瀑布溝庫首右岸拉裂變形體支護(hù)加固前可能形成的不同規(guī)?；?，并進(jìn)行滑坡涌浪模擬，對于支護(hù)加固后拉裂變形體的穩(wěn)定性有待后期作進(jìn)一步研究。

參考文獻(xiàn)：

［1］夏婷，鄭聲安，任艷，等.全球重點(diǎn)國家和區(qū)域水電開發(fā)程度對標(biāo)分析［J］.水力發(fā)電學(xué)報(bào)，2022，41（5）：1-11.

［2］羅涵，巨莉，羅茂盛.四川省水電開發(fā)基地水電開發(fā)程度評估［J］.水利水電科技進(jìn)展，2023，43（1）：71-80.

［3］ERTUN? A.The geological problems of the large dams constructed on the Euphrates River［J］.Engineering Geology，1999，51（3）：167-182.

［4］RIEMER W.Landslides and reservoirs［C］//Proceedings 6^thInternational Symposium on Landslides.Rotterdam：Balkema Publishers，1995：213-224.

［5］竇思軍，徐興倩，李繼國，等.降雨條件下水庫滑坡變形特征及形成機(jī)制分析［J］.人民珠江，2020，41（9）：32-38.

［6］楊渠鋒，王平義，喻濤，等.三峽庫區(qū)陡巖滑坡涌浪爬高試驗(yàn)分析［J］.中國地質(zhì)災(zāi)害與防治學(xué)報(bào)，2014，25（3）：43-48.

［7］曹婷，王平義，胡杰龍.基于物理模型試驗(yàn)的庫區(qū)巖質(zhì)滑坡涌浪爬高研究［J］.南水北調(diào)與水利科技，2018，16（5）：159-164.

［8］肖莉麗，殷坤龍，王佳佳，等.基于物理模擬試驗(yàn)的庫岸滑坡沖擊涌浪［J］.中南大學(xué)學(xué)報(bào)（自然科學(xué)版），2014，45（5）：1618-1626.

［9］KOO W，KIM M H.Numerical modeling and analysis of waves induced by submerged and aerial/sub-aerial landslides［J］.KSCE Journal of Civil Engineering，2008，12（2）：77-83.

［10］彭輝，吳凡，金科，等.庫岸滑坡涌浪首浪高度試驗(yàn)研究［J］.水利水電技術(shù)，2017，48（12）：95-100.

［11］LI Y，CHEN L X，YIN K L，et al.Quantitative risk analysis of the hazard chain triggered by a landslide and the generated tsunami in the Three Gorges Reservoir area［J］.Landslides，2021，18（2）：667-680.

［12］謝海清，蔣昌波，鄧斌，等.狹窄型庫區(qū)河道滑坡涌浪的形成及其傳播規(guī)律［J］.交通科學(xué)與工程，2017，33（4）：45-50.

［13］黃筱云，劉燦，程永舟，等.V型河道下滑坡涌浪的傳播與爬高［J］.長沙理工大學(xué)學(xué)報(bào)（自然科學(xué)版），2017，14（1）：70-77.

［14］鄧成進(jìn)，黨發(fā)寧，陳興周.庫區(qū)滑坡涌浪傳播及其與大壩相互作用機(jī)理研究［J］.水利學(xué)報(bào)，2019，50（7）：815-823.

［15］徐文杰.滑坡涌浪流-固耦合分析方法與應(yīng)用［J］.巖石力學(xué)與工程學(xué)報(bào)，2020，39（7）：1420-1433.

［16］潘家錚.建筑物的抗滑穩(wěn)定和滑坡分析［M］.北京：水利出版社，1980．

［17］殷坤龍，劉藝梁，汪洋，等.三峽水庫庫岸滑坡涌浪物理模型試驗(yàn)［J］.地球科學(xué)（中國地質(zhì)大學(xué)學(xué)報(bào)），2012，37（5）：1067-1074.

（責(zé)任編輯：高天揚(yáng)）