摘 要：在對庫區(qū)滑坡涌浪資料的整理分析基礎(chǔ)上，概述典型庫區(qū)滑坡涌浪實(shí)例及其災(zāi)害。滑坡涌浪會給庫區(qū)周圍居民的生命財(cái)產(chǎn)安全帶來極大威脅，為了更好地防治滑坡涌浪災(zāi)害，該文從滑坡涌浪爬高方面出發(fā)，對近年來滑坡涌浪在物理模型實(shí)驗(yàn)和數(shù)值模擬研究進(jìn)展進(jìn)行綜述，在物理模型實(shí)驗(yàn)方面，研究者們通過調(diào)整不同的參數(shù)，如滑坡形狀、斜度、體積等，模擬滑坡涌浪的爬高過程，并得出一些重要的結(jié)論。在數(shù)值模擬方面，研究者們通過建立數(shù)學(xué)模型，采用數(shù)值方法計(jì)算滑坡涌浪的運(yùn)動和爬升過程。這種方法可以很好地理解滑坡涌浪的動力學(xué)行為，并預(yù)測其爬高過程。未來應(yīng)該繼續(xù)深入研究影響滑坡涌浪爬高行為的因素，以期取得更多的突破和進(jìn)展。

關(guān)鍵詞：滑坡涌浪；爬高；物理模型實(shí)驗(yàn)；數(shù)值模擬；綜述

中圖分類號：P642.22 文獻(xiàn)標(biāo)志碼：A 文章編號：2095-2945（2024）27-0115-04

Abstract： Based on the collation and analysis of the data of surge landslides in the reservoir area， this paper summarizes the typical examples of landslide surges and their disasters in the reservoir area. The landslide surge will bring great threat to the life and property safety of the residents around the reservoir area. In order to better prevent the landslide surge disaster， this paper starts from the aspect of landslide surge climbing， this paper summarizes the research progress of landslide surge in physical model experiment and numerical simulation in recent years. In physical model experiment， researchers simulate the climbing process of landslide surge by adjusting different parameters， such as landslide shape， slope， volume and so on， and some important conclusions are drawn. In the aspect of numerical simulation， researchers establish mathematical models and use numerical methods to calculate the movement and climbing process of landslide surges. This method can well understand the dynamic behavior of landslide surge and predict its climbing process. In the future， we should continue to study the factors that affect the climbing behavior of landslide surge， in order to make more breakthroughs and progress.

Keywords： landslide surge; climbing; physical model experiment; numerical simulation; summary

隨著經(jīng)濟(jì)的快速發(fā)展，中國在水利水電工程領(lǐng)域?qū)崿F(xiàn)了迅猛發(fā)展。特別是高壩建設(shè)技術(shù)，已經(jīng)逐步趨于成熟，并開始在全球范圍內(nèi)占據(jù)領(lǐng)先地位。目前，我國已經(jīng)成功建造了超過五百座中大型水庫，其中不乏世界級的重要水利工程，例如三峽、錦屏、白鶴灘和烏東德等。這些高壩建成后形成的大型水庫，在降雨、地震及水庫水位變化等多種因素的影響下會產(chǎn)生庫岸滑坡，從而導(dǎo)致涌浪的產(chǎn)生?；掠坷俗鳛橐环N由滑坡體滑入水中引發(fā)的次生災(zāi)害，其破壞力往往遠(yuǎn)超滑坡本身，經(jīng)常會對周圍居民、沿岸建筑物以及設(shè)備設(shè)施造成極大損失，甚至是慘痛的災(zāi)難性后果。每年全世界有數(shù)百人因?yàn)榛掠坷硕鴨适?963年10月9日，意大利的瓦伊昂水庫因?yàn)樾钏僮鲗?dǎo)致左岸約2.7億 m3的巖土體失穩(wěn)并滑入庫區(qū)，由此產(chǎn)生的涌浪在對岸攀升了大約250 m的高度，漫過壩頂超100 m，摧毀大壩壩體內(nèi)所有設(shè)備設(shè)施，在1分鐘左右的時(shí)間造成約3 000人死亡[1]；1971年3月，秘魯?shù)膴W德斯由于巖石崩塌而形成的滑坡涌浪，在數(shù)秒鐘內(nèi)吞沒了丘加礦區(qū)400～600人，同時(shí)損壞了大量設(shè)施設(shè)備[2]；1985年6月，長江三峽庫區(qū)新灘鎮(zhèn)北岸由于巖崩而造成約3 000萬m3的滑坡體沖入長江，形成的滑坡涌浪在對岸的爬高為49 m，導(dǎo)致2人失蹤，8人受傷，10人死亡[3]；在2009年7月，由于持續(xù)降雨的影響，瀾滄江上的小灣水電站遭遇了一次嚴(yán)重的滑坡涌浪事故，其中約有100萬m3的土體滑入了瀾滄江中。這次滑坡產(chǎn)生了超過30 m高的涌浪，不幸造成了14人失蹤[4]。表1中詳細(xì)列出了我國歷史上的典型滑坡涌浪災(zāi)害案例。

表1 我國典型滑坡涌浪災(zāi)害實(shí)例

產(chǎn)生滑坡涌浪的過程是復(fù)雜的，按照其發(fā)展變化，主要分為3個階段：產(chǎn)生、傳播、爬高。在第一階段，即涌浪的起始階段，滑坡體在眾多因素的共同影響下開始滑動，其原本具有的較高勢能轉(zhuǎn)變?yōu)閯幽?，?dāng)這些動能沖擊水面時(shí)，便形成了能量巨大的涌浪。這些涌浪能夠迅速破壞附近的建筑和設(shè)施設(shè)備，給社會造成巨大損失[5]。第二個階段是涌浪的傳播階段，相較于涌浪的產(chǎn)生階段，其破壞力更大，影響的范圍更廣，由于能量傳遞作用，導(dǎo)致河道中的庫水位會瞬間上升，產(chǎn)生水位差而形成壓力差，同時(shí)還會混雜水體回流以及波浪疊加等現(xiàn)象，導(dǎo)致涌浪流動不均勻而增大破壞力。涌浪的傳播不僅會對沿岸的建筑物和船舶造成威脅，還會對大壩壩面和庫區(qū)岸坡產(chǎn)生極強(qiáng)的沖刷，使其穩(wěn)定性降低。第三個階段是涌浪的爬高階段，對于狹窄的河道，涌浪爬高階段的破壞力高于前兩個階段[6]。涌浪爬高一旦超過大壩，會給沿岸居民的生命安全帶來極大的威脅，更嚴(yán)重是會沖刷大壩，影響壩體結(jié)構(gòu)穩(wěn)定安全，給沿岸的居民和建筑物帶來極大的危害。目前，對于滑坡涌浪爬高研究的系統(tǒng)梳理較少，本文基于以往物理模型實(shí)驗(yàn)和數(shù)值模擬，對滑坡涌浪爬高的研究成果進(jìn)行了總結(jié)并提出展望。研究成果將豐富滑坡涌浪研究體系，提高后續(xù)研究的有效性與針對性。

1 物理模型實(shí)驗(yàn)

Mueller等[7]基于物理模型實(shí)驗(yàn)造出了脈沖波列，并研究了在爬坡角度β為18.4、45和90°時(shí)，滑坡涌浪首波的爬高和漫壩過程。Synolakis[8]對孤立波爬高進(jìn)行了研究，提出了非破碎波理論的近似解，并做了大量實(shí)驗(yàn)來進(jìn)行比較，結(jié)果表明非線性理論能夠較好地預(yù)測孤立波的爬高。然而，Synolakis提出的近似解僅適用于爬坡角度β為2.9°。Evers等[9]研究了在爬坡角度β為10～90°范圍內(nèi)脈沖波的爬高，發(fā)現(xiàn)現(xiàn)有的預(yù)測方程得出的理論值與實(shí)際結(jié)果相比偏小。Lo等[10]針對斜坡上的單、雙孤立波的爬高和反射過程，發(fā)現(xiàn)第一波產(chǎn)生的反射強(qiáng)度會影響第二波的爬高高度，第二波的爬高高度隨第一波的分離時(shí)間而變化。彭勇等[11]研究了在斜坡上多個孤立波的爬高和相互作用規(guī)律，實(shí)驗(yàn)結(jié)果表明第二個孤立波與第三個孤立波的爬高放大系數(shù)與相鄰的波峰間距有關(guān)，同時(shí)根據(jù)實(shí)驗(yàn)結(jié)果提出了等波幅三孤立波在緩坡上的爬坡經(jīng)驗(yàn)公式。汪洋等[12]根據(jù)流體力學(xué)方程、非恒定流的連續(xù)性方程和沿程水頭損失理論，以新灘滑坡為例對滑坡涌浪的爬坡高度進(jìn)行計(jì)算。殷坤龍團(tuán)隊(duì)[13]運(yùn)用正交試驗(yàn)法，綜合考量了滑坡體量、滑動面角度、岸邊坡度、水深和入水速度等多個因素對滑坡涌浪產(chǎn)生的影響。他們以三峽庫區(qū)的白水河滑坡事件為原型，構(gòu)建了一個按1∶200縮小比例的河道物理模型，通過實(shí)驗(yàn)研究，他們成功推導(dǎo)出了適用于三峽庫區(qū)滑坡涌浪爬高的計(jì)算公式[13]。S？覸levik等[14]研究者基于透水性泥沙斜坡的分析，根據(jù)實(shí)驗(yàn)結(jié)果建立了孤立波爬高的經(jīng)驗(yàn)公式。Fuchs團(tuán)隊(duì)[15]利用二維水槽的實(shí)體模型實(shí)驗(yàn)，深入探究了孤立波在傳播過程中的特性及其爬坡行為。McFall等[16]則利用三維物理實(shí)驗(yàn)手段，構(gòu)建了以鵝卵石作為滑體材料的物理實(shí)驗(yàn)?zāi)Ｐ?，以此來研究孤立波的爬高現(xiàn)象，其物理模型如圖1所示，通過此裝置來研究孤島以及水平滑道下的滑體運(yùn)動和滑坡涌浪爬坡規(guī)律。Liyanage等[17]研究人員探討了在不同的波高、水深和斜坡角度條件下滑坡涌浪的爬坡特性，通過進(jìn)行實(shí)驗(yàn)并結(jié)合現(xiàn)有的實(shí)驗(yàn)數(shù)據(jù)與模型分析，構(gòu)建了一個用于預(yù)測孤立波爬升高度的模型。

2 數(shù)值模擬研究

Titov研究團(tuán)隊(duì)[18]采用了有限差分法來求解非線性方程組，以此構(gòu)建了孤立波在斜坡上破碎并爬升的數(shù)值模型，為了解決數(shù)值穩(wěn)定性問題，他們對邊界條件進(jìn)行了調(diào)整，盡管如此，波浪破碎位置的數(shù)值模擬結(jié)果仍顯示出振蕩現(xiàn)象。Grilli等[19]運(yùn)用邊界元方法對孤立波的傳播行為及其爬升高度進(jìn)行了研究，該模型能夠模擬水體粒子運(yùn)動的垂直分布，但是在處理波浪破碎時(shí)產(chǎn)生的射流與自由表面的碰撞時(shí)，計(jì)算過程會中斷，因此該模型無法計(jì)算孤立波在岸邊的最大爬升高度。黃波林團(tuán)隊(duì)[20]運(yùn)用水動力學(xué)模型，并融合地理信息系統(tǒng)（GIS）技術(shù)，共同研發(fā)了滑坡涌浪災(zāi)害的快速評估系統(tǒng)軟件，簡稱FAST，該系統(tǒng)能夠根據(jù)不同的工況迅速評估滑坡涌浪的傳播特性和爬升高度，以三峽地區(qū)茅草坡滑坡為例，如圖2所示，通過對比系統(tǒng)計(jì)算結(jié)果與實(shí)際滑坡事件的數(shù)據(jù)，對FAST系統(tǒng)的實(shí)用性進(jìn)行了驗(yàn)證。Barranco等[21]對潰壩生成的鉆孔進(jìn)行了數(shù)值模擬，提出了潰壩滑坡涌浪爬高公式。楊渠鋒等[22]以三峽斗巖滑坡作為研究對象，模擬了在不同的堆放組合情況下矩形滑坡體涌浪的產(chǎn)生和傳播過程，得出了適合陡巖滑坡的涌浪爬高公式。Viroulet等[23]采用了SPHysics和Gerris兩種數(shù)值模型來研究滑坡涌浪問題，SPHysics是基于拉格朗日無網(wǎng)格法來精確描述撞擊階段，Gerris是基于兩相有限體積法來研究滑坡涌浪波的傳播和爬高特性。首先模擬了2個經(jīng)典算例：二維楔形塊體沿斜面下滑入水試驗(yàn)和滑塊垂直入水試驗(yàn)，試驗(yàn)數(shù)據(jù)與兩個模型計(jì)算結(jié)果對比吻合較好。

圖1 鵝卵石滑坡物理模型

圖2 茅草坡位置圖

3 展望

滑坡涌浪作為一種由滑坡體入水引發(fā)的次生災(zāi)害，其爬高特性對于評估和預(yù)防水庫區(qū)災(zāi)害風(fēng)險(xiǎn)具有重要意義。本文通過系統(tǒng)梳理滑坡涌浪爬坡特性的研究現(xiàn)狀，并分析影響其特性的多種因素。由于滑坡涌浪的爬坡特性受到多因素的影響，目前的理論分析、物理實(shí)驗(yàn)和數(shù)值模擬方法均無法全面涵蓋所有相關(guān)因子。因此，對庫區(qū)滑坡涌浪的爬高特性研究必須采取具體問題具體分析的方法。研究人員可以在以下幾個方面對滑坡涌浪的爬高特性進(jìn)行深入研究。

1）大部分物理模型試驗(yàn)?zāi)壳岸际窃诙S情況下進(jìn)行的。在這種二維條件下，地形邊界對滑坡涌浪產(chǎn)生的反射、衍射和折射效應(yīng)很難被充分考慮進(jìn)去，導(dǎo)致涌浪爬高預(yù)測的結(jié)果較為單一，并且無法充分考慮到沿岸地區(qū)的影響。三維涌浪模型因其能夠提供更為全面和精確的模擬結(jié)果，正逐漸成為該領(lǐng)域研究發(fā)展的重要趨勢。

2）依托于工程實(shí)例，探究滑坡涌浪爬高特性，對滑坡涌浪的爬高進(jìn)行預(yù)測并給出指導(dǎo)性意見。

3）滑坡體的物理特性是影響涌浪爬坡特性的關(guān)鍵因素之一?；麦w的體積、質(zhì)量、形狀以及滑坡物質(zhì)的顆粒大小和密度都會對涌浪的形成和爬高產(chǎn)生顯著影響。研究人員可以通過收集和分析不同類型滑坡事件的數(shù)據(jù)，建立滑坡體物理特性與涌浪爬坡特性之間的關(guān)聯(lián)模型。

4）滑坡體入水時(shí)的初始條件，如入水速度、入水角度和入水深度等也是影響涌浪爬高特性的重要因素。這些初始條件決定了涌浪的起始動能和動量，從而影響其在斜坡上的爬升能力。為了更好地理解這些初始條件對涌浪爬高特性的影響，對于不同入水條件下的涌浪行為還需要深入展開研究。

5）地形條件，尤其是滑坡發(fā)生區(qū)域的岸邊地形，對涌浪的爬坡特性有著直接的影響。岸邊地形的坡度、長度和材料特性都會改變涌浪的傳播路徑和爬升能力。通過對不同地形條件下的涌浪爬坡特性進(jìn)行研究，可以為水庫區(qū)的災(zāi)害預(yù)防和應(yīng)急響應(yīng)提供科學(xué)依據(jù)。

6）環(huán)境因素，如風(fēng)、波浪和降雨，也會對滑坡涌浪的爬坡特性產(chǎn)生影響。這些環(huán)境因素可能會改變水體的表面條件，從而影響涌浪的傳播和爬升。因此，考慮這些環(huán)境因素在涌浪爬坡過程中的作用并將其納入模型中，對于準(zhǔn)確預(yù)測涌浪爬高具有重要意義。

7）數(shù)值模擬方法的發(fā)展為滑坡涌浪爬坡特性的研究提供了新的工具。通過采用先進(jìn)的數(shù)值模擬技術(shù)，如有限元分析、有限差分法和邊界元法等，研究人員可以在計(jì)算機(jī)中模擬復(fù)雜的涌浪行為，從而更準(zhǔn)確地預(yù)測涌浪的爬坡高度。然而，數(shù)值模擬的準(zhǔn)確性依賴于輸入?yún)?shù)的精確性和模型的合理性，因此，對模型參數(shù)的校準(zhǔn)和驗(yàn)證是未來研究滑坡涌浪的一個重點(diǎn)方向。

滑坡涌浪的爬高特性研究是一個多學(xué)科、多因素的復(fù)雜問題。研究人員需要綜合考慮滑坡體物理特性、滑坡體入水條件、環(huán)境因素以及數(shù)值模擬方法等多個方面，才能更全面地理解和預(yù)測滑坡涌浪的爬高特性。未來的研究應(yīng)該著重于提高模型的準(zhǔn)確性，完善參數(shù)的測量和校準(zhǔn)方法，以及發(fā)展更為精確和實(shí)用的數(shù)值模擬技術(shù)，以便更好地服務(wù)于水庫區(qū)滑坡涌浪災(zāi)害的預(yù)防和應(yīng)對工作。

參考文獻(xiàn)：

[1] PANIZZO A， DE G P， DI R M， et al. Great landslide events in Italian artificial reservoirs[J]. Natural Hazards and Earth System Sciences， 2005，5（5）：733-740.

[2] COUSTON L A， MEI C C， ALAM M R. Landslide tsunamis in lakes[J].Journal of Fluid Mechanics，2015，772：784-804.

[3] 殷坤龍，姜清輝，汪洋.新灘滑坡運(yùn)動全過程的非連續(xù)變形分析與仿真模擬[J].巖石力學(xué)與工程學(xué)報(bào)，2002，21（7）：959-962.

[4] 石傳奇，安翼，楊家修.滑坡涌浪的三維SPH方法模擬及其工程應(yīng)用[J].中國科學(xué)：物理學(xué) 力學(xué) 天文學(xué)，2015（10）：52-60.

[5] LEVIN B W， NOSOV M. Physics of tsunamis[M].Dordrecht： Springer，2009.

[6] 李榮輝，蔣昌波，鄧斌，等.近壩庫區(qū)狹窄河道滑坡涌浪高度及其傳播規(guī)律[J].交通科學(xué)與工程，2016，32（2）：79-84.

[7] MUELLER D R. Wave run-up and overtopping on dams; Auflaufen und Ueberschwappen von Impulswellen an Talsperren[J].1995.

[8] SYNOLAKIS C E. The runup of solitary waves[J].Journal of Fluid Mechanics，1987，185：523-545.

[9] EVERS F M， BOES R M. Impulse wave runup on steep to vertical slopes[J]. Journal of Marine Science and Engineering，2019，7（1）：8.

[10] LO H Y， PARK Y S， LIU P L F. On the run-up and back-wash processes of single and double solitary waves-An experimental study[J]. Coastal engineering，2013，80：1-14.

[11] 彭勇，吳衛(wèi)，劉樺.三孤立波在緩坡上爬高及相互作用的實(shí)驗(yàn)研究[J].水動力學(xué)研究與進(jìn)展：A輯，2018，33（4）：445-454.

[12] 汪洋，殷坤龍.水庫庫岸滑坡涌浪的傳播與爬高研究[J].巖土力學(xué)，2008，29（4）：1031-1034.

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

[14] S？覸LEVIK G， JENSEN A， PEDERSEN G. Experimental investigation of impact generated tsunami; related to a potential rock slide， Western Norway[J].Coastal Engineering， 2009，56（9）：897-906.

[15] FUCHS H. Solitary impulse wave run-up and overland flow[R]. Swiss Federal Institute of Technology ETH，2013.

[16] MCFALL B C， FRITZ H M. Runup of granular landslide-generated tsunamis on planar coasts and conical islands[J]. Journal of Geophysical Research： Oceans， 2017，122（8）：6901-6922.

[17] LIYANAGE CVT， TANAKA N. A unified runup formula for solitary waves on a plane beach[J]. Ocean Engineering， 2020， 216： 108038.

[18] TITOV V V， SYNOLAKIS C E. Modeling of breaking and nonbreaking long-wave evolution and runup using VTCS-2[J]. Journal of Waterway， Port， Coastal， and Ocean Engineering，1995，121（6）：308-316.

[19] GRILLI S T， SVENDSEN I A， SUBRAMANYA R. Breaking criterion and characteristics for solitary waves on slopes[J]. Journal of waterway， port， coastal， and ocean engineering， 1997，123（3）：102-112.

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

[21] BARRANCO I， LIU P L F. Run-up and inundation generated by non-decaying dam-break bores on a planar beach[J]. Journal of Fluid Mechanics， 2021，915：A81.

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

[23] VIROULET S， CéBRON D， KIMMOUN O， et al. Shallow water waves generated by subaerial solid landslides[J]. Geophysical Journal International， 2013，193（2）：747-762.