海底滑坡海啸的颗粒流耦合模型

海底滑坡的运动可能引发海啸,破坏离岸设施,威胁海岸带安全。国内外关于海底滑坡引发海啸的研究方兴未艾。采用Mih颗粒流模型控制具弱黏聚力的砂土滑坡运动,利用两相流模型计算岩土体-水体相互作用及RNG湍流模型控制水体运动,构建了基于颗粒流模型的海底滑坡海啸全耦合数值分析方法。通过简单水槽水下滑坡案例进行了海底滑坡海啸全过程研究。数值分析再现了变形滑体的不均一运动、密度分异流动、水滑机制和以波谷为典型特征的涌浪波等典型海底滑坡及海啸现象,这表明数值模型具有有效性。许多海域（包括中国南海北部）都存在弱黏聚力和无黏聚力的水下滑坡,该数值方法值得推广和进一步研究完善。     

A rigorous finite volume model to simulate subaerial and submarine landslide-generated waves

This paper presents a new landslide-generated wave (LGW) model based on incompressible Euler equations with Savage-Hutter assumptions. A two-layer model is developed including a layer of granular-type flow beneath a layer of an inviscid fluid. Landslide is modeled as a two-phase Coulomb mixture. A well-balanced second-order finite volume formulation is applied to solve the model equations. Wet/dry transitions are treated properly using a modified non-linear method. The numerical model is validated using two sets of experimental data on subaerial and submarine LGWs. Impulsive wave characteristics and landslide deformations are estimated with a computational error less than 5 %. Then, the model is applied to investigate the effects of landslide deformations on water surface fluctuations in comparison with a simpler model considering a rigid landslide. The model results confirm the importance of both rheological behavior and two-phase nature of landslide in proper estimation of generated wave properties and formation patterns. Rigid slide modeling often overestimates the characteristics of induced waves. With a proper rheological model for landslide, the numerical prediction of LGWs gets more than 30 % closer to experimental measurements. Single-phase landslide results in relative errors up to about 30 % for maximum positive and about 70 % for maximum negative wave amplitudes. Two-phase constitutive structure of landslide has also strong effects on landslide deformations, velocities, elongations, and traveling distances. The complex behaviors of landslide and LGW of the experimental data are analyzed and described with the aid of the robust and accurate finite volume model. This can provide benchmark data for testing other numerical methods and models.     

Risk management study on impulse waves generated by Hongyanzi landslide in Three Gorges Reservoir of China on June 24, 2015

On June 24, 2015, Hongyanzi slope located in Wushan County of the Three Gorges Reservoir collapsed, generating 5–6-m-high impulse waves, which overturned 13 boats, killed 2 persons, and injured 4 persons. It is the second incident of landslide-generated impulse waves since the 175-m experimental impoundment in 2008. The emergency investigation shows that Hongyanzi landslide is a bedding soil landslide with a volume of 23?×?10⁴ m³ induced by a series of triggering factors such as rainfall, flooding upstream, and reservoir drawdown. The nonlinear Boussinesq water wave model is used to reproduce the impulse waves generated by the landslide of June 24th. The numerical simulation results suggest that the wave propagation process was influenced by the T-shaped geomorphic conditions of river valley, and the coastal areas in the county seat were the major wave-affected areas, which is opposite to the landslide. The numerical wave process accord well with the observed incident, and the investigation values were in good agreement with the calculated values. Moreover, the worst-case scenario of the 7?×?10⁴ m³ deformation mass beside Hongyanzi landslide is potential to generate impulse waves, which was predicted with the same numerical model. This adjacent deformation mass will probably generate impulse waves with maximum height and run-up of 2.2 and 2.0 m, respectively, and only a very few areas in the water course had waves rising to a height of 1 m or above. The research results provide a technical basis for emergency disposal to Hongyanzi landslide and navigation restriction in Wushan waterway. More importantly, it pushes the risk management of the navigation based on the impulse wave generated by landslide. It is advised that the Three Gorges Reservoir and other reservoirs around the world should put more efforts in performing special surveys and studies on the potential hazards associated with landslide-generated impulse waves.     

水库滑坡涌浪传播有限元数值模拟

库岸滑坡涌浪危害巨大,正确预测库区可能的滑坡涌浪非常重要,是工程可行性论证的重要内容之一。将浅水控制方程应用于滑坡涌浪数值模拟,控制方程采用简便且具有较高精度的两步Taylor-Galerkin方法进行求解。通过算例对数值模型的应用进行验证,结果表明,涌浪产生后将以入水点为源点迅速向四周推进并不断衰减,且随传播距离的增加浪高降幅逐渐减小。计算的涌浪高度及水位变化规律与实测资料吻合得很好,并将滑坡涌浪的沿程传播过程可视化。研究结果表明,文中方法模拟滑坡涌浪传播是有效可行的,可用于滑坡涌浪灾害的预测和防治。     

Modeling the interaction between fluid flow and particle movement with discontinuous deformation analysis (DDA) method

Coupled hydro‐mechanical processes in granular media represent the interaction between solid particles movement and fluid flow during external and/or internal loading. The processes attract attention in geotechnical engineering since they cause many disasters such as landslide, slope collapse, boiling/quicksand and soil liquefaction, etc. This study presents a numerical method combining discontinuous deformation analysis (DDA) for mechanical calculation and finite element method for fluid flow simulation to model the interaction between solid particles' movement and fluid flow from microscopic point of view. The term of hydraulic pressure was newly formulated and introduced to the original DDA. The simplified examples are given to verify the new method, and the computational results correlate well with the theoretical calculations, further development is also considered. Copyright © 2010 John Wiley & Sons, Ltd.     

Analysis of waves generated by Gongjiafang landslide in Wu Gorge, three Gorges reservoir, on November 23, 2008

At 4:40p.m. on November 23, 2008, the Gongjiafang slope collapsed on the north bank of Yangtze River in Wu Gorge of Three Gorges Reservoir. The 380,000-m³ sliding mass consisted mainly of cataclastic rock. A video record of the major sliding incident was analyzed using the general laws of physical motion. The analysis indicated that the maximum speed and maximum acceleration of the sliding mass were 11.65?m/s and 2.23?m/s², respectively, and that the maximum amplitude and the propagation velocity of the water wave near the landslide were 31.8?m and 18.36?m/s, respectively. Wave run-up investigation indicated that the maximum run-up on shore was 13.1?m, which declined to 1.1?m at Wushan dock 4?km away. The incident causes no casualties, but did result in economic losses of RMB five million. The numerical simulation model GEOWAVE was used to simulate and reproduced the impulse wave generated by the landslide; the results were in good agreement with the observed incident. The numerical simulation data were then applied to analyze the decay and amplification effects of the landslide wave in the river course. The field investigations and witness information provide valuable materials for the studies of landslide kinematics and impulse waves generated by landslides. In addition, the research results provide a useful reference for future similar waves generated by landslides in reservoirs.     

基于偏振特性的下孔法剪切波到时判别问题探讨

基于偏振特征的剪切波到时判别在下孔法波速测试中被广泛应用，然而实际测试中时有发现压缩波初动出现反相的现象，因此针对压缩波不偏振而剪切波偏振的理论基础进行论证很有必要。通过三维有限元分析模型模拟了地表敲击下孔法波速测试，并基于时域逐步积分动力有限元方法求解了地表正反向水平敲击下不同深度观测点振动响应。结果表明：地表水平正反向激励下，各深度观测点中与激励方向平行的水平分量中压缩波也出现反相。为进一步探讨这与实际工程的认知相悖的现象，数值模拟了荷载倾斜、探头旋转、探头倾斜等实际测试可能出现的情形。结果显示：荷载倾斜与探头倾斜的综合影响会导致正反向敲击激励下压缩波的初动相位不反相。进一步，提出了一种剪切波旅时计算新方法以减少每次敲击之间敲击角度、探头偏转和倾斜变化带来的影响，数值计算得到基于新方法求解的剪切波速接近模型介质剪切波速值，其优于利用峰－峰法和互相关法确定旅时而得到的波速结果。     

非结构网格上的三维浅水流动数值模型

针对当前复杂环境水流模拟的需求,建立了新型的基于特征型高分辨率数值算法的三维非结构网格浅水动力模型。模型采用有限体积法离散sigma坐标下的三维浅水方程,运用Roe黎曼近似解评估水平界面通量。模型网格拟合边界能力强,可根据需要局部加密;格式数值性能优良,具有守恒性、单调迎风性、高数值分辨率等特性。同时,应用干湿判别法处理动边界,以适应浅滩地形漫/露过程模拟的需要。封闭水池内部风生环流、干河床上溃坝过程和长江口实际潮流场的模拟从不同侧面展示了模型的特点,结果表明它能够准确地预测水流的三维流动结构,而且计算简单高效,具有良好的数值稳定性。     

电阻率层析成像法探测滑坡面正演模拟与反演成像研究

滑坡灾害已成为危及人民生命财产安全的重大问题,确定滑坡体和滑坡面的分布是滑坡治理的关键因素,电阻率层析成像法因对物性差异的界面响应敏感,被广泛应用到滑坡面的探测中。首先,对现有的滑坡地电模型进行概化,得到了4种常见的滑坡的简化模型;其次,针对不同类型的滑坡,采用有限元方法进行了系统的数值正演研究,得到了层析成像法探测滑坡面的响应特征,提出了不同类型滑坡面的快速识别方法。然后,采用最小二乘法开展了反演研究,总结了层析成像法探测滑坡面的反演成像规律,结合数值正演的响应特征,提出了层析成像法探测滑坡面的异常识别特征。最后,在山东沂源西里镇滑坡进行了层析成像法探测,成功探明了滑坡面的分布,表明了层析成像法探测滑坡面是行之有效的,为滑坡的治理方案设计奠定了基础     

A pore-scale numerical model for non-Darcy fluid flow through rough-walled fractures

This paper aims at developing a pore-scale numerical model for non-Darcy fluid flow through rough-walled fractures. A simple general relationship between the local hydraulic conductivity and the flow velocity is proposed. A new governing equation for non-Darcy fluid flow through rough-walled fractures is then derived by introducing this relationship into the Reynolds equation. Based on the non-linear finite element method, a self-developed code is used to simulate the non-Darcy fluid flow through fractures. It is found that the macroscopic results obtained by the numerical simulation agree well with the experimental results. Furthermore, some interesting experimental observations can be reproduced.     

煤矿水平孔复杂地质模型方位电磁波响应数值模拟

煤岩界面的预先、精准识别是实现煤矿巷道自动化掘进和煤矿智能化的关键技术之一,方位电磁波仪器探测距离较大,能够分辨岩性界面和界面方位,在油田测井中已取得成功应用,但在煤矿领域应用较少。为了研究方位电磁波测井在煤矿复杂三维(3D)模型中的探测性能,使用有限元数值模拟方法,考察了方位电磁波测井在煤矿测量环境中,不同发射频率和源距时受钻孔内流体的影响,以及采空区和起伏地层边界对方位电磁波测量的影响。结果表明:对于煤矿测井中的常见情况——孔中流体为空气时,方位电磁波响应受钻孔的影响较小;方位电磁波测量信号能够反映地层界面的起伏变化,其探测能力与发射频率、源距有关;方位电磁波响应对低阻采空区较敏感,但在高频发射时可以利用相位差信号实现对高阻采空区的探测。方位电磁波方法为煤矿井下水平孔煤岩界面和采空区测量提供了一种新技术手段,具有较大的应用潜力。      

A 2D integrated FEM model for surface water–groundwater flow of slopes under rainfall condition

Numerical modeling of water infiltration in slopes under rainfall conditions, especially under rainstorm conditions, is a fundamental problem for slope stability assessment. To obtain representative results, surface water–groundwater flow models are incorporated in the simulation. Based on finite element representation of Richards’ equation and of kinematic wave equations, an integrated 2D numerical model (IMCR2D) of the surface water–groundwater system was established. The model has a symmetrical matrix that modifies the flux boundary according to the runoff solution on the slope. IMCR2D was verified using two laboratory experiments, and it showed good agreement with numerical and experimental results. Additional numerical examples were used to study the effect of flux supply from runoff on infiltration. In comparison with SimMd (an existing method), IMCR2D displayed advantages in cases where surface runoff develops in an upper low-permeability section of the slope and flows down into a high-permeability section of the slope. To illustrate the advantages of the new method, the seepage field and stability condition of a case study in the Three Gorges Hydroelectric Reservoir were analyzed using IMCR2D and SimMd. The deformation of a landslide in part reflects its stability, and therefore, we also used displacement monitoring data to estimate the variation of stability conditions from that aspect. Comparison of the two numerical models indicated that flux supply greatly affects the seepage field, and that rainfall plays an important role in landslide stability evaluation, but only when considering flux supply from upper slope surface runoff.     

Computational fluid dynamics modelling of landslide generated water waves

This paper describes the application of detailed computational fluid dynamics (CFD) to simulate the formation and propagation of waves generated by the impact of landslide material with water. The problem is schematised as a multiphase–multicomponent fluid flow: compressible air, water and transported alluvial material. The landslide simulation is performed by means of a hybrid approach: as a rigid solid body slipping down along an inclined slope until it starts penetrating the water body. The CFD model solves the Navier–Stokes equations with the RNG k-ɛ turbulence closure scheme and the volume of fluid multiphase method, which maintains the interface as a sharp front. The governing equations are solved using the commercial CFD code, FLUENT. The computed results are compared with experimental data reported in the literature. The model is then applied to simulate the 1958 Lituya bay Tsunami event with a 2D a simplified geometry and the results are compared to others found in literature.     

福建德化霞碧滑坡渗流稳定性模拟分析

降雨入渗导致土体孔隙水压力的增大造成土体有效抗剪强度的降低,进而发生滑坡。本文建立滑坡饱和-非饱和渗流有限元模型,考虑孔隙水压力随时间和空间任意不规则分布,计算特定降雨量、降雨强度、降雨历时不同滞后时段基于渗流场的滑坡稳定性。     

Comprehensive analyses of the initiation and landslide-generated wave processes of the 24 June 2015 Hongyanzi landslide at the Three Gorges Reservoir,China

Reservoir landslides pose a great threat to shipping safety, human lives and properties, and the operation of the hydropower station. In this paper, the 24 June 2015 Hongyanzi landslide at the Three Gorges Reservoir is considered as an example to study the initiation mechanism and landslide-generated wave process of a reservoir landslide. The finite difference method and limit equilibrium analysis are used to analyze the deformation and failure characteristics of the Hongyanzi slope. Simulation results show that a large deformation (about 358 mm) happens in the shallow deposits under intermittent rainfall condition, and the slope is in a limit state. At the same time, continuous rapid drawdown of the water level (about ?0.55 m/day during 8–24 June 2015) reduced the support and accelerated the drainage of the water for the bank slope. A coupling effect of intermittent rainfall and rapid drawdown of the water level was the triggering factor of the 24 June Hongyanzi landslide. Landslide-generated wave process was simulated using a fluid–solid coupling method by integrating the general moving object collision model. Simulation results show that the landslide-generated wave is dominated by the impulse wave, which is generated by sliding masses entering the river with high speed. The maximum wave height is about 5.90 m, and the wave would decay gradually as it spreads because of friction and energy dissipation. To prevent reservoir landslides, the speed for the rising or drawdown of the water level should be controlled, and most importantly, rapid drawdown should be avoided.     

基于颗粒流方法的堆积层滑坡运动过程模拟

以颗粒流离散元为研究方法对勉县杨家湾十组堆积层滑坡破坏方式与运动过程进行数值模拟研究。通过PFC2D双轴模拟试验所标定的岩土体宏观模拟参数与室内试验所获取的宏观实测参数进行对比,确定堆积层滑坡所需的颗粒细观参数,然后将标定的细观参数代入堆积层滑坡模型,对滑坡破坏方式及运动过程进行模拟研究。结果表明:滑坡破坏在初始阶段蠕滑变形累积,滑坡体挤压坡脚,直至坡脚产生剪切破坏,并向上牵引发展,使得滑坡整体顺接触面破坏下滑并堆积于坡脚,表现为典型的牵引式渐进破坏,结果与实际情况基本吻合。研究认为,采用颗粒流方法对堆积层滑坡破坏与运动过程的模拟研究具有较高的适用性,对该类滑坡防治具有一定的参考意义。      

黄岗滑坡稳定性评价及数值分析

根据对黄岗滑坡的分析,其属于蠕滑-拉裂型滑坡。本文阐述了滑坡的基本特征,并研究分析了滑坡的成因机制及演化过程;在成因机制分析基础上,采用极限平衡法对该滑坡天然状态及暴雨情况下的稳定性进行了计算和评价;另外,还采用有限元数值分析软件Midas对滑坡进行数值模拟,并对数值模拟成果和传统计算的结果进行对比分析,为滑坡治理设计提供了重要的依据。     

The August 27, 2014, rock avalanche and related impulse water waves in Fuquan,Guizhou, China

At about 8:30 p.m. on 27 August 2014, a catastrophic rock avalanche suddenly occurred in Fuquan, Yunnan, southwestern China. This landslide and related impulse water waves destroyed two villages and killed 23 persons. The impulse waves occurred after initiation of the landslide, caused by the main part of the slide mass rapidly plunging into a water-filled quarry below the source area. The wave, comprising muddy water and rock debris, impacted the opposite slope of the quarry on the western side of the runout path and washed away three homes in Xinwan village. Part of the displaced material traveled a horizontal distance of about 40 m from its source and destroyed the village of Xiaoba. To provide information for potential landslide hazard zonation in this area, a combined landslide–wave simulation was undertaken. A dynamic landslide analysis (DAN-W) model is used to simulate the landslide propagation before entering the quarry, while Fluent (Ansys Inc., USA) is used to simulate the impulse wave generation and propagation. Output data from the DAN-W simulation are used as input parameters for wave modeling, and there is good agreement between the observed and simulated results of the landslide propagation. Notably, the locations affected by recordable waves according to the simulation correspond to those recorded by field investigation.     

三峡树坪滑坡动力学的有限元模拟

三峡工程是迄今为止最大的水利工程,对库区滑坡灾害的监测和机制研究一直是重要的研究课题。本文利用Terra SAR-X的强度图进行相关计算,求解出2009年5月20日至8月5日期间三峡树坪滑坡的形变场。该形变场特征和树坪滑坡体的地形特征吻合甚好,位移大小、方向和三峡大学对滑坡体的野外观测结果基本吻合。以此高精度位移场为外部约束,结合野外观测资料对滑坡体介质力学性质进行分类并选取边界条件,利用有限元方法对滑坡活动进行动力学计算模拟。计算过程中对滑坡体的滑动面形状、因降雨引起材料参数变化和三峡水库水位等因素分别反演和调整,得出符合其变形和发展过程的滑坡动力学特征。发现软弱带的物性参数决定滑坡体总体滑动量,滑坡体的物性参数决定位移分布的峰值位置。在确定了滑坡动力学特征之后,进一步讨论降雨和库区水位下降对滑坡产生的贡献权重,得出降雨是树坪滑坡的决定因素。