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.     

A non-hydrostatic model for the numerical study of landslide-generated waves

This paper proposes and demonstrates a two-layer depth-averaged model with non-hydrostatic pressure correction to simulate landslide-generated waves. Landslide (lower layer) and water (upper layer) motions are governed by the general shallow water equations derived from mass and momentum conservation laws. The landslide motion and wave generation/propagation are separately formulated, but they form a coupled system. Our model combines some features of the landslide analysis model DAN3D and the tsunami analysis model COMCOT and adds a non-hydrostatic pressure correction. We use the new model to simulate a 2007 rock avalanche-generated wave event at Chehalis Lake, British Columbia, Canada. The model results match both the observed distribution of the rock avalanche deposit in the lake and the wave run-up trimline along the shoreline. Sensitivity analyses demonstrate the importance of accounting for the non-hydrostatic dynamic pressure at the landslide-water interface, as well as the influence of the internal strength of the landslide on the size of the generated waves. Finally, we compare the numerical results of landslide-generated waves simulated with frictional and Voellmy rheologies. Similar maximum wave run-ups can be obtained using the two different rheologies, but the frictional model better reproduces the known limit of the rock avalanche deposit and is thus considered to yield the best overall results in this particular case.     

基于二维浅水方程的滑坡体兴波数值模型

基于有限体积方法和结构化网格,建立了海底滑坡引起的波浪传播数值模型。模型控制方程为考虑了海床随时间变化的二维浅水方程。采用中心迎风格式计算控制体界面数值通量,采用线性重构技术、局部海床高程处理技术和全隐式离散底摩阻项,保证了格式的和谐性、守恒性和水深非负性,有效处理了海岸动边界问题。时间积分采用具有强稳定性质的二阶龙格-库塔方法（ssp RK）。针对滑坡体兴波经典算例开展数值模拟,将计算结果与解析解、实验结果及其他模型计算结果进行比较和分析。结果表明,对于所考虑的计算工况,模型能较合理地模拟滑坡兴波的产生、传播和爬高过程。     

A new equation for estimating the maximum surface settlement above tunnels excavated in soft ground

Nowadays, due to urbanization and population increase, need for metro tunnels, has been considerably increased in urban areas. Common characterization of urban area tunnels is that they are excavated in very shallow depths and soft ground. In such excavations, main challenge for tunneling is low bearing capacity and easy deformation characteristic of the ground. Tunnel face instability and the potential surface settlement are the most hazardous factors that should be considered in all tunneling methods applied in urban areas. Incorrect estimation of the maximum surface settlement value can lead to irreparable damages to the buildings and other nearby structures. There are several published relationships concerned with field measurements and analytical solutions to estimate the amount of the maximum surface settlement value due to tunneling. These relationships are not precise for calculating the aimed values. Therefore, providing accurate equations for estimation of these values is certainly useful. First purpose of this study is to determine the effective parameters such as geotechnical factors (cohesion, internal friction angle, density, Young’s modulus and Poisson’s ratio), and engineering factors (tunnel depth, tunnel diameter and face support pressure) on the maximum surface settlement value. In this study, three metro project constructions namely Istanbul, Tehran, Mashhad in the Middle East were chosen. FLAC3D (Itasca Consulting Group 2002) was used for detailed numerical analysis. The second aim is to present better equations in estimating the maximum surface settlement-based actual data set from several tunnel projects and numerical modeling. The results from the new estimation equation are compared with results of empirical and field observations. The maximum surface settlement values obtained from the new equation have good agreement with the actual results for three different metro case studies.     

约束条件对滑坡涌浪特征影响的室内试验研究

水库滑坡约束条件影响其运动过程的几何形态,是滑坡涌浪预测的重要参数之一.为了探究约束条件对滑坡涌浪特征(波高、波幅与周期)的影响,采用正交试验设计法开展了 54组滑坡涌浪室内模型试验,并基于统计学理论对约束散体和半约束散体的涌浪特征进行了分析.结果表明:涌浪波周期基本不受滑体约束条件的影响;而半约束散体模型的波高和波幅...     

水库滑坡涌浪风险研究现状和灾害链风险管控实践

滑坡作为水库库区主要地质灾害类型之一,其风险研究一直是近年来的研究热点.水库滑坡涌浪的产生使滑坡灾害的影响范围由滑坡源本身扩散到上下游数千米,极大地扩大了滑坡风险的承灾体类型与数量以及灾害损失程度.因涉及交叉学科领域,滑坡涌浪风险评估是滑坡风险灾害链评价的难点与前沿课题.本文综合了前人近几十年来的研究成果,首先从危险性...     

Extension of the empirical Alyavdin equation for representing batch grinding data

The Alyavdin equation for batch grinding data is:

$\text{1 ? P(χ, t) = [1 ? P(χ, 0)]}\text{exp}\text{?}\text{c(x)t}{}^{\mathrm{p}}\text{]}$

where P(χ,t) is the weight fraction less than size χ after grinding time t, c (χ) is constant with t and p is a constant close to one. It is shown that this equation is illogical (except for a single size of feed) unless c (χ) varies with P(χ,0), which makes the equation of little utility. A new empirical equation is developed for finite size intervals:

which reduces to the Alyavdin equation for a single size of feed, and which gives consistent computations for any feed size distribution. Techniques are given for determining K_i, γ values from sets of batch grinding data. The values are then used to predict size distributions for other times and other feed size distributions. The equation was quite successful in predicting size distributions in batch milling: (a) providing the feed size distribution was not un-natural, that is, not truncated or (b) if a truncated feed was used, the values of K_i and γ are determined from size distributions of grinding of the same type of feed. Thus, K_i, γ are not, unfortunately, completely independent of the starting feed size distribution.     

A numerical model in predicting the initial karst development in porous limestone

The numerical model proposed in this paper represents a new method in predicting initial karst spatial distribution of carbonate aquifers. Karstification Index, \( \varPhi_{\text{K}} \), is a measure of the possible largest amount of calcite that can be dissolved from a unit area of rock and transported by groundwater flow during a unit length of time. Mathematically, it is defined as \( \varPhi_{\text{K}} = q \cdot s \), where q is the flux of groundwater and s is the rock solubility; and Karstification Index, \( \varPhi_{\text{K}} \), is the function of flow, heat, and calcium concentration. The numerical model simulates the time- and space-dependent susceptibility to, and process of, karstification by iterative calculations between \( \varPhi_{\text{K}} \) and hydraulic conductivity. A sample calculation is presented on a hypothetical limestone plateau with steep escarpment, a spring, and a river on one side, a vertical no flow boundary on the other side, and a horizontal impermeable base. The rock has an equivalent porous-medium hydraulic conductivity of K?=?5E?8 ms^?1, which is characteristic to non-karstified limestone. Initially, the spring is located at the intersection of the water table and the escarpment’s face some distance above the river. The simulation result showed that the water table is lowered in 2325 years across the plateau, but the spring’s drop is stopped after 65 years. Karstification is most intensive along the water table due to the higher value of Karstification Index which resulted from the intensive flux of the water. The continual decline of the water table results in a laterally elongated and progressively expanding body of karstified rock and finally as a horizontal conduit at river level. The concentration of Ca²⁺ and conduit diameter decrease from the spring toward the impermeable vertical wall on the opposite side and downward to the base of the plateau. The results are consistent with general theories and empirical observations of karst and cave evolution.     

跨孔电波衰减成像初始振幅估算方法比较

跨孔电磁法衰减层析成像是一种利用电磁波振幅信息的方法,通过发射端到接收端电磁波的振幅变化来反演介质衰减常数分布。发射端振幅也就是初始振幅,一般情况下是未知的,它的精度很大程度上影响到层析成像结果,需要在反演前得到或者通过特殊反演方法来处理。本文总结了4种初始振幅处理方法——线性拟合法、矩阵反演法、双频电磁波法以及相邻道比值法,通过合成数据验证了这四种方法的可行性,并且指出了每种方法的优缺点:线性拟合法适合物性变化不大的情况;矩阵反演法对物性情况要求不高,但计算量较大;双频电磁波法能直接得到电导率分布,但只适合良导体情况;相邻道比值法适用情况最广,但容易受干扰影响。     

Evaluating accuracy for two empirical methods in predicting settlement of drilled shafts

Several theoretical, empirical and semi-empirical methods are available in the literature to predict settlement of drilled shafts in sandy soils. In the Arabian Gulf countries, specifically in the United Arab Emirates, equations and procedure from the rest of the world are being used in analysis and design of drilled shafts without proper validation. It is the aim of this study to assess the applicability and evaluate the accuracy of two well known, and commonly used methods for pile prediction in the United Arab Emirates (UAE), namely Vesic (1977) and Poulos (1979), via comparison with data from field pile load tests conducted on shafts drilled in the region. Some of these tests were conducted for the purpose of this study, while others were made available through the courtesy of International Piling Contractors who are active in the region (e.g. Bauer International and Swiss Borings). Pile load test data were analyzed to back-calculate the model parameters related to settlement under different loading stages. Geological data and soil properties were obtained from studies conducted at the relevant sites. An effort is made to correlate soil properties with the prediction models. Statistical analysis is conducted to assess the accuracy of the results obtained from the two methods at different stages of loading via those obtained from pile load tests. Moreover, a detailed parametric study is conducted to assess the effect of the related parameters on the predicted pile settlement and the estimated settlement at different stages of loading. The study concluded with a recommendation of the most appropriate models and procedures to be followed for predicting the settlement of drilled shafts in the UAE, together with useful charts and correlation relations. Results showed that settlement values predicted by Vesic (1977) and Poulos (1979) overestimates the true values. This revised version was published online in July 2006 with corrections to the Cover Date.     

A new X-ray investigation of maximum microcline crystalline solutions

Refined unit-cell dimensions and molar volumes of a mircocline-low albite ion-exchange series previously studied by Waldbaum and Robie (1971) have been redetermined. Members in the sodic half of the series were found to have significantly larger b dimensions than those previously reported, resolving the major difference between the above study and that of Orville (1967). W _v, the Margules parameter related to the volume of mixing, is significantly lower (0.078 cal/bar mol) than those determined in either of the investigations above. However, it is close to values of W _v determined for several alkali feldspar series with topochemically monoclinic Al-Si distributions. Thus, the dependence of K-Na volume of mixing properties on Al-Si distribution in alkali feldspars is doubtful.     

一种新的钻井液静失水量方程

实验室测试钻井液的API失水量结束后,仪器内剩余钻井液的粘度比测试前有所增大。为了解释这一现象,对原钻井液静失水量方程推导中所采用的渗滤模式提出了不同看法,以一种不同的渗滤模式为基础提出了新的钻井液静失水量方程,此方程可以解释这一现象。依照此方程,实验室测试API失水量时仍然可用7.5 min失水量的测试结果乘以2得到的数值代替30 min失水量的测试结果。     

A new exact solution of Kompaneets equation for a shock front

An analytical solution of the Kompaneets equation is found for a medium whose density varies according to a hyperbolic tangent law, from a more diffuse (interstellar) medium to a denser medium (a cloud). The law of motion of the leading points of the shock front is discussed. Intermediate asymptotics describing the acceleration and change of the shape of the shock front are analyzed in detail. The obtained solution can be used to investigate analytically the evolution of the shock-front shape at a boundary with a molecular cloud.     

An empirical equation of effective shaking duration for moderate to large earthquakes

Natural Hazards - The duration of strong shaking is particularly important for assessing building performance, potential landslides and liquefaction hazards. The results of this investigation can...     

考虑滑脱效应的低阶煤动态渗透率预测新模型

为准确预测低阶煤动态渗透率变化规律,在煤岩立方体模型基础上,考虑基质孔隙和滑脱效应对渗透率的影响,建立低阶煤动态渗透率预测新模型,并对影响绝对渗透率和滑脱系数的因素进行敏感性分析,讨论了甲烷和氮气对基质收缩与滑脱效应的影响。研究表明：基于“火柴棍”假设建立的模型是新模型不考虑基质孔隙时的一个特例,P-M、S-D模型与新模型相比基质收缩作用更加明显,考虑基质收缩与滑脱效应的新模型更具实用性。气体郎格缪尔应变是影响基质收缩的关键,煤岩绝对渗透率能否反弹是割理压缩、基质孔隙膨胀、基质弹性形变和基质收缩4个因素共同作用的结果。相同条件下,甲烷的基质收缩强于氮气,氮气的滑脱效应强于甲烷,影响滑脱系数的因素包括内因和外因,滑脱系数与割理宽度随孔隙压力变化时呈现相反规律。滑脱效应和基质收缩效应共同提升气测渗透率,煤岩孔隙压力越低,二者对渗透率的提升作用越明显。     

Physical model tests to determine the mechanism of submarine landslides under the effect of sea waves

Submarine landslides are a common type of disaster which threaten property and the safety of human life. To effectively prevent and control such disasters, we conduct a series of large-scale physical model tests to determine the mechanism of submarine landslides. First, a large-scale physical model test system is designed and developed, including flume test frame, wave-making system, wave-absorbing system, and data monitoring system. In the tests, we investigate the effect of different sea waves by changing the parameters of the wave-making system and the influence of the slope inclination by constructing different models. Data regarding the wave pressure acting on the slope surface, seepage pressure, and displacement are monitored during the test procedure. The test results show that the seepage pressure in the faults varies cyclically with the sea waves and is lower at internal points than at outcrops. If the wave loading time is sufficiently long, the seepage pressure and displacement deformation in the fault zone will gradually increase. In other words, failures in fault zones precede submarine landslides. The weak fault zone provides the preferred sliding surface, and the sea waves supply the external dynamic energy for submarine landslides. The conclusions provide guidelines for similar engineering and research.

     

The southern Rhine graben: A new view of the initial phase

One of the puzzling features of the southern end of the Rhine graben is the Dinkelberg-Tabular Jura block on the eastern shoulder of the graben. It is dissected by a large number of faults, the most notable ones forming a field of narrow little grabens and half-grabens whose bordering faults converge at the level of the Middle Triassic evaporites, which points to décollement at that horizon. The little grabens were traditionally considered to be of Oligocene age, coeval with the main taphrogenesis of the Rhine graben. Two hypotheses were offered for their formation, one ascribing them to extension on the extrados of large basement folds, the other to gravity sliding on paleoslopes. Recent field work uncovered overwhelming evidence for an Eocene age of the little grabens, the time of the initial phase of Rhine graben formation. At that time there were neither large basement folds nor paleoslopes of any significance, and therefore the two hypotheses offered until now do not work. However, the map-view pattern of the field of faults offers a somewhat unusual way out of the dilemma. This pattern is most prominently displayed in the Dinkelberg area north of the Rhine. There a lane of narrow décollement grabens with a mean NNE strike is confined within the NW- striking Dinkelberg graben, which is much wider and rooted in the basement. It is also very shallow, with a subsidence on the order of 100 m. The lane of décollement grabens forms a dextral en-échelon pattern with respect to the Dinkelberg graben, suggesting stretching of the post-evaporite sequence above a basement essentially extended by strike slip. This model, though not as clearly expressed, is also compatible with the data in the rest of the Dinkelberg-Tabular Jura block. It also fits surprisingly well a theoretical model by Withjack and Scheiner (1982) that predicts a dominance of strike-slip in the marginal area of a system consisting of extension superimposed on doming.     

A new estimate of global methane flux from onshore and shallow submarine mud volcanoes to the atmosphere

A new estimate of global methane emission into the atmosphere from mud volcanoes (MVs) on land and shallow seafloor is presented. The estimate, considered a lower limit, is based on 1) new direct measurements of flux, including both venting of methane and diffuse microseepage around craters and vents, and 2) a classification of MV sizes in terms of area (km²) based on a compilation of data from 120 MVs. The methane flux to the atmosphere is conservatively estimated between 6 and 9 Mt y^–1. This emission from MVs is 3–6% of the natural methane sources and is comparable with ocean and hydrate sources, officially considered in the atmospheric methane budget. The total geologic source, including MVs, seepage from seafloor, microseepage in hydrocarbon-prone areas and geothermal sources, would amount to 35–45 Mt y^–1. The authors believe it is time to add this parameter in the Intergovernmental Panel on Climate Change official tables of atmospheric methane sources.GEM     

A new empirical estimator of coseismic landslide displacement for Zagros Mountain region (Iran)

Earthquake-induced landslides are responsible worldwide for significant socioeconomic losses and historically have a prominent position in the list of natural hazards affecting the Iran plateau. As a step toward the development of tools for the assessment and the management of this kind of hazard at regional scale, an empirical estimator of coseismic displacements along potential sliding surfaces was obtained through a regression analysis for the Zagros region, a mountainous Iranian region subjected to earthquake-induced landslides. This estimator, based on the Newmark’s model, allows to evaluate the expected permanent displacement (named “Newmark displacement”) induced by seismic shaking of defined energy on potential sliding surface characterized by a given critical acceleration. To produce regression models for Newmark displacement estimators, a data set was constructed for different critical acceleration values on the basis of 108 accelerometric recordings from 80 Iranian earthquakes with moment magnitudes between 3.6 and 7. The empirical estimator has a general form, proposed by Jibson (Eng Geol 91:209–218, 2007), relating Newmark displacement to Arias intensity (as parameter representing the energy of the seismic forces) and to critical acceleration (as parameter representing the dynamic shear resistance of the sliding mass). As an example of application, this relation was employed to provide a basic document for earthquake-induced landslide hazard assessment at regional scale, according to a method proposed by Del Gaudio et al. (Bull Seismol Soc Am 93:557–569, 2003), applied to the whole Iranian territory, including Zagros region. This method consists in evaluating the shear resistance required to slopes to limit the occurrence of seismically induced failures, on the basis of the Newmark’s model. The obtained results show that the exposure to landslide seismic induction is maximum in the Alborz Mountains region, where critical accelerations up to ~0.1 g are required to limit the probability of seismic triggering of coherent type landslides within 10% in 50 years.     

Development of an empirical model for predicting peak breach flow of landslide dams considering material composition

The existing empirical models do not consider the influence of material composition of landslide deposits on the peak breach flow due to the uncertainty in the material composition and the randomness of its distribution. In this study, based on the statistical analyses and case comparison, the factors influencing the peak breach flow were comprehensively investigated. The highlight is the material composition-based classification of landslide deposits of 86 landslide cases with detailed grain-size distribution information. In order to consider the geometric morphology of landslide dams and the potential energy of dammed lakes, as well as the material composition of landslide deposits in an empirical model, a multiple regression method was applied on a database, which comprises of 44 documented landslide dam breach cases. A new empirical model for predicting the peak breach flow of landslide dams was developed. Furthermore, for the same 44 documented landslide dam failures, the predicted peak breach flow obtained by using the existing empirical models for embankment and landslide dams and that obtained by using the newly developed model were compared. The comparison of the root mean square error (E_rms) and the multiple coefficient of determination (R²) for each empirical model verifies the accuracy and rationality of the new empirical model. Furthermore, for fair validation, several landslide dam breach cases that occurred in recent years in China and have reliable measured data were also used in another comparison. The results show that the new empirical model can reasonably predict the peak breach flow, and exhibits the best performance among all the existing empirical models for embankment and landslide dam breaching.