蒙特卡罗模拟在区域地震滑坡灾害评价中应用

2004年10月23日,日本新泻中越地区发生史上最强地震之一,震级达到Mw6.6。这次地震诱发上千个规模大小不一、形式各样的滑坡,造成一定的灾害损失和人员伤亡。因此,为了有效地避免和减轻这方面的灾害,有必要有效地、科学地预测和评价地震诱发滑坡的空间发生规律。为此以一定的地区作为研究对象,提出了蒙特卡罗随机模拟和Newmark滑动体位移分析法在区域地震斜坡失稳概率分析中的应用。该法充分考虑到岩土参数和相关地震参数空间变异性,结合地理信息系统空间分析的优势,以网格模型为基础,分析该区斜坡滑动体位移分布规律,并以2cm为失稳临界值探讨网格斜坡破坏的概率。实例证明,所提出的模型有效地预测了滑坡发生的空间分布规律。     

Assessment of shallow landslide susceptibility using the transient infiltration flow model and GIS-based probabilistic approach

This study proposes a probabilistic analysis method for modeling rainfall-induced shallow landslide susceptibility by combining a transient infiltration flow model and Monte Carlo simulations. The spatiotemporal change in pore water pressure over time caused by rainfall infiltration is one of the most important factors causing landslides. Therefore, the transient infiltration hydrogeological model was adopted to estimate the pore water pressure within the hill slope and to analyze landslide susceptibility. In addition, because of the inherent uncertainty and variability caused by complex geological conditions and the limited number of available soil samples over a large area, this study utilized probabilistic analysis based on Monte Carlo simulations to account for the variability in the input parameters. The analysis was performed in a geographic information system (GIS) environment because GIS can deal efficiently with a large volume of spatial data. To evaluate its effectiveness, the proposed analysis method was applied to a study area that had experienced a large number of landslides in July 2006. For the susceptibility analysis, a spatial database of input parameters and a landslide inventory map were constructed in a GIS environment. The results of the landslide susceptibility assessment were compared with the landslide inventory, and the proposed approach demonstrated good predictive performance. In addition, the probabilistic method exhibited better performance than the deterministic alternative. Thus, analysis methods that account for uncertainties in input parameters are more appropriate for analysis of an extensive area, for which uncertainties may significantly affect the predictions because of the large area and limited data.     

利用模糊数学方法建立汶川地震滑坡灾害评判标准

利用汶川地震主震在龙门山震区的强震资料,建立了利用地震动参数估计斜坡位移的Newmark累积位移模型。借助于模糊数学中的隶属函数方法,实现了对模型估计累积位移的隶属度定量计算,进而用以判定滑坡危险性分级。为使计算结果定量化,进一步利用模糊转换方法,将定性的滑坡危险性判定结果转换为定量的概率表达,使滑坡判定结果可以继续应用于其他环节的计算分析。以汶川地震在龙门山震区的滑坡分布为例,对模型的计算结果及其定量转换结果与实际情况进行了分析,表明模型结果与实际情况非常吻合。分析认为,以概率方式可以使滑坡灾害的空间分布描述更为准确,文中的方法使之得以实现。     

A Hazard Assessment Method for Potential Earthquake-Induced Landslides – A Case Study in Huaxian County, Shaanxi Province

The hazard assessment of potential earthquake-induced landslides is an important aspect of the study of earthquake-induced landslides. In this study, we assessed the hazard of potential earthquake-induced landslides in Huaxian County with a new hazard assessment method. This method is based on probabilistic seismic hazard analysis and the Newmark cumulative displacement assessment model. The model considers a comprehensive suite of information, including the seismic activities and engineering geological conditions in the study area, and simulates the uncertainty of the intensity parameters of the engineering geological rock groups using the Monte Carlo method. Unlike previous assessment studies on ground motions with a given exceedance probability level, the hazard of earthquake-induced landslides obtained by the method presented in this study allows for the possibility of earthquake-induced landslides in different parts of the study area in the future. The assessment of the hazard of earthquake-induced landslides in this study showed good agreement with the historical distribution of earthquake-induced landslides. This indicates that the assessment properly reflects the macroscopic rules for the development of earthquake-induced landslides in the study area, and can provide a reference framework for the management of the risk of earthquake-induced landslides and land planning.     

Seismic-event-based methodology to obtain earthquake-induced translational landslide regional hazard maps

A seismic-event-based methodology to generate earthquake-induced translational landslide maps using Newmark method is proposed. The steps are: (1) to construct a GIS-based geotechnical database; (2) to identify those areas that are susceptible to the occurrence of translational landslides based on available geological information; (3) to compute a static safety factor; (4) to compute the critical acceleration that defines the threshold acceleration required to cause a displacement; (5) to characterize the seismic hazard as a set of stochastic events, collectively exhaustive and mutually exclusive, that fully describes the hazard spatial distribution and annual frequency of occurrence (in accordance with the earthquake location, depth and magnitude) with the appropriate ground-motion prediction equations; (6) to compute the Newmark displacement; and finally, (7) to carry out a probabilistic translational landslide hazard analysis to estimate an exceedance rate of a given displacement. This methodology is applied to Mexico, and maps for return periods of 150 and 500 years are presented. Results shown in maps are estimations of where translational landslides may occur and should be useful to carry out local studies to elaborate recommendations of site specific hazard reduction plans as well as to calculate insurance rates. In addition, these results are useful to identify civil protection actions, risk management at regional and local level, and land use planning, as well as for promoting more detailed vulnerability and risk studies at different scales.     

基于简化Newmark模型的长白山天池火山诱发崩塌滑坡危险性评价

火山喷发过程所伴生的地震活动会诱发大量的崩塌滑坡次生灾害,其所造成的人员财产损失甚至超过火山活动本身。2002年以来长白山天池火山区地震活动的异常,表明火山深部的岩浆正在发生变化,天池火山存在喷发的危险。地震崩塌滑坡的危险性区划是降低生命财产损失的有效手段。将火山伴生地震作为崩塌滑坡灾害的诱发因素并据此设置地震参数,利用简化的Newmark累积位移模型,考虑地形因素对地震的放大效应,对长白山地区天池火山喷发下次生崩塌、滑坡灾害的危险性进行评价。通过探讨不同地震震级下的危险性分区结果,认为不同地震参数的设置对危险性分区结果没有影响。将研究区划分为极高、高、中等、低、极低等5个危险等级,其中,极高危险区主要分布在3个区域:以天池口为中心,40km为半径的范围内;沿江乡—两江镇—松江镇条带区域;长白县境内鸭绿江沿岸区域。     

Recent landslide activity in Manaihan,East Belgium

Past landslides have been recognized in the Battice area in E-Belgium. In contrast to the other inactive landslides, the Manaihan landslide responded immediately to heavy rainfall events in the last two decades. This study aims to map its spatial extent and the dominant surface features; to measure surface displacement using GPS; to investigate subsurface structure with Cone penetration test (CPT) and corings; and to determine the depth of the shear surface by inclinometers. Results show a partial landslide reactivation. Surface velocities range between 20 and 40 cm/year and are strongly dependent on winter rainfall. CPT results give clear boundaries between the landslide mass and the undisturbed bedrock in the head scarp. Distinct shear surfaces have been determined with displacement rates up to 15.8 mm in 21 days. Further research should apply geophysical methods for two-dimensional information on the ground, investigate geotechnical properties of the landslide mass, model slope instability, and determine the influence of a sewage pipe crossing the central landslide mass as a potential cause for landslide activity.     

Landslide susceptibility mapping using Bayesian approach in the Sultan Mountains (Akşehir,Turkey)

Landslides cause heavy damage to property and infrastructure, in addition to being responsible for the loss of human lives in many parts of the Turkey. The paper presents GIS-based spatial data analysis for landslide susceptibility mapping in the regions of the Sultan Mountains, West of Akşehir, and central part of Turkey. Landslides occur frequently in the area and seriously affect local living conditions. Therefore, spatial analysis of landslide susceptibility in the Sultan Mountains is important. The relationships between landslide distributions with the 19 landslide affecting parameters were analysed using a Bayesian model. In the study area, 90 landslides were observed. The landslides were randomly subdivided into 80 training landslides and 10 test landslides. A landslide susceptibility map was produced by using the training landslides. The test landslides were used in the accuracy control of the produced landslide susceptibility map. Approximately 9% of the study area was classified as high susceptibility zone. Medium, low and very low susceptibility zones covered 8, 23 and 60% of the study area, respectively. Most of the locations of the observed landslides actually fall into moderate (17.78%) and high (77.78. %) susceptibility zones of the produced landslide susceptibility map. This validates the applicability of proposed methods, approaches and the classification scheme. The high susceptibility zone is along both sides of the Akşehir Fault and at the north-eastern slope of the Sultan Mountains. It was determined that the surface area of the Harlak and Deresenek formations, which have attained lithological characteristics of clayey limestone with a broken and separated base, and where area landslides occur, possesses an elevation of 1,100–1,600 m, a slope gradient of 25°–35° and a slope aspect of 22.5°–157.5° facing slopes.     

Influence of strong motion characteristics on permanent displacement of slopes

Landslides triggered by moderate to major earthquakes are a recognized seismic hazard. Arias Intensity (I _a) is a key intensity measure of the ground motion, but significant duration is widely used to define strong motion duration. We calculate Newmark’s displacements using earthquake records bracketing a broad range of Arias Intensity and significant duration employing Newmark’s rigid block method and a number of yield accelerations. Total landslide displacement increases with the increase in the energy content of the ground motion (I _a) above a threshold. Such threshold may be expressed as a function of yield acceleration of the slope regardless of the ground motion characteristics. Newmark displacement decreases with increasing duration for earthquakes with similar energy content. The wide scatter in the results converges when using formal dimensional analysis. Self-similar symmetry may facilitate the assessment of the performance of slopes during earthquakes. The mathematical framework for probabilistic determination of landslides displacement may be a useful aid to estimate the likelihood of landslide hazards provided that the geotechnical properties of the slopes are known.     

基于不同位移预测模型的地震滑坡危险性评估研究-以天水地区为例

近年来,Newmark累积位移分析方法经过不断的改进和应用成为国际主流的地震滑坡危险性评估方法之一,众多学者基于位移预测模型开展区域地震滑坡危险性评估,然而鲜有针对不同位移模型对评估结果影响的定量研究。以天水地区为例,基于不同的位移预测模型开展地震滑坡危险性评估,对比位移模型对地震滑坡危险性评估的影响,探讨建立适用于我国的Newmark位移预测模型。结果表明:基于不同位移预测模型评估所得的地震滑坡危险性结果整体趋势一致,均能区分区域地震滑坡危险性等级的相对差异,但在同样的危险性分级标准下,所得中、高危险区的分布范围有较大差异。这与位移模型的函数形式及其区域相关性有关,在引入Newmark累积位移分析方法开展地震滑坡危险性评估的同时,应尽快建立考虑地震动衰减特征和工程地质背景的Newmark位移预测模型,为中国潜在地震滑坡危险性预测评估、震后滑坡快速评估等提供技术支撑。      

The 2010 Yushu earthquake triggered landslide hazard mapping using GIS and weight of evidence modeling

The Yushu County, Qinghai Province, China, April 14, 2010, earthquake triggered thousands of landslides in a zone between 96°20′32.9″E and 97°10′8.9″E, and 32°52′6.7″N and 33°19′47.9″N. This study examines the use of geographic information system (GIS) technology and Bayesian statistics in creating a suitable landslide hazard-zone map of good predictive power. A total of 2,036 landslides were interpreted from high-resolution aerial photographs and multi-source satellite images pre- and post-earthquake, and verified by selected field checking before a final landslide-inventory map of the study area could be established using GIS software. The 2,036 landslides were randomly partitioned into two subsets: a training dataset, which contains 80 % (1,628 landslides), for training the model; and a testing dataset 20 % (408 landslides). Twelve earthquake triggered landslide associated controlling parameters, such as elevation, slope gradient, slope aspect, slope curvature, topographic position, distance from main surface ruptures, peak ground acceleration, distance from roads, normalized difference vegetation index, distance from drainages, lithology, and distance from all faults were obtained from variety of data sources. Landslide hazard indices were calculated using the weight of evidence model. The landslide hazard map was compared with training data and testing data to obtain the success rate and predictive rate of the model, respectively. The validation results showed satisfactory agreement between the hazard map and the existing landslide distribution data. The success rate is 80.607 %, and the predictive rate is 78.855 %. The resulting landslide hazard map showed five classes of landslide hazard, i.e., very high, high, moderate, low and very low. The landslide hazard evaluation map should be useful for environmental recovery planning and reconstruction work.     

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

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

Centrifugal model test on a riverine landslide in the Three Gorges Reservoir induced by rainfall and water level fluctuation

Frequent soil landslide events are recorded in the Three Gorges Reservoir area, China, making it necessary to investigate the failure mode of such riverside landslides. Geotechnical centrifugal test is considered to be the most realistic laboratory model, which can reconstruct the required geo-stress. In this study, the Liangshuijing landslide in the Three Gorgers Reservoir area is selected for a scaled centrifugal model experiment, and a water pump system is employed to retain the rainfall condition. Using the techniques of digital photography and pore water pressure transducers, water level fluctuation is controlled, and multi-physical data are thus obtained, including the pore water pressure, earth pressure, surface displacement and deep displacement. The analysis results indicate that: Three stages were set in the test (waterflooding stage, rainfall stage and drainage stage). Seven transverse cracks with wide of 1–5 mm appeared during the model test, of which 3 cracks at the toe landslide were caused by reservoir water fluctuation, and the cracks at the middle and rear part were caused by rainfall. During rainfall process, the maximum displacement of landslide model reaches 3 cm. And the maximum deformation of the model exceeds 12 cm at the drainage stage. The failure process of the slope model can be divided into four stages: microcracks appearance and propagation stage, thrust-type failure stage, retrogressive failure stage, and holistic failure stage. When the thrust-type zone caused by rainfall was connected or even overlapped with the retrogressive failure zone caused by the drainage, the landslide would start, which displayed a typical composite failure pattern. The failure mode and deformation mechanism under the coupling actions of water level fluctuation and rainfall are revealed in the model test, which could appropriately guide for the analysis and evaluation of riverside landslides.     

四川九寨沟Ms7.0级地震滑坡应急快速评估

基于地震滑坡危险性评估的Newmark累积位移模型,利用震前获取的震区地形数据、区域地质资料,结合地震动近实时获取技术,开展了四川九寨沟M_s7.0级地震诱发滑坡的应急快速评估。地震滑坡位移分析结果表明,同震滑坡活动的中—高强度区分布在断层两侧宽约4 km的带状区域内,整体沿北西方向延伸。其中,极震区的丰雪塘、日则和干海子等城镇驻地及附近道路的滑坡强度相对较高;震前、震后影像对比表明九寨沟地震诱发的滑坡类型以浅表型碎屑流及小规模崩塌为主,且同震碎屑流多是在震前已有碎屑流的基础上进一步活动扩展而来,震后汛期泥石流隐患也不容忽视;通过典型地区滑坡位移分析结果与震前、震后影像对比,表明滑坡位移分析结果能够较好的反映同震滑坡的宏观分布特征,但在场地尺度上吻合程度欠佳,后续将通过提升岩性和地形等数据质量进行改进。研究结果可为灾情研判提供宝贵信息,对提高灾害应急救援效率具有重要意义。     

The contribution of PSInSAR interferometry to landslide hazard in weak rock-dominated areas

In the Grande da Pipa river basin, north of Lisbon, 64 % of the total number of landslides inventoried is totally or partially included in a lithological unit composed by marl, clay, and sandstone intercalation complex that is present in 58 % of the study area. The Persistent Scatterer synthetic aperture radar interferometry technique is applied to a data set of TerraSAR-X SAR images, from April of 2010 to March of 2011, firstly to the Laje-Salema test site and further exported to the Grande da Pipa river basin. This work’s specific objectives are the following: (i) to assess the potential of the Persistent Scatterer displacement maps to the identification of new landslides/unstable areas and in the redefinition of landslide limits, (ii) to update the landslide state of activity, and (iii) to evaluate the capacity of the Persistent Scatterer deformation maps in assessing landslide susceptibility at the regional scale. Based on this approach, it was possible to increment the number of landslides and to redefine the landslide limits in the test site in 3.8 %. For 39 landslides, it was possible to update the landslide state of activity, in particular from dormant to reactivated or dormant-reactivated (23 landslides) or from stabilized to reactivated (5 landslides). Landslide susceptibility map based in Persistent Scatterer deformation rates, independently validated with a deep rotational slide map, obtained the best value of area under the curve (0.668).     

Regional landslide susceptibility zoning with considering the aggregation of landslide points and the weights of factors

In this paper, we propose a methodology for landslide susceptibility assessment at a regional scale in Yunnan, southwestern province of China. A landslide inventory map including 3,242 landslide points was prepared for the study area. Five factors recognized as correlated to landslide (namely, lithology, relative relief, tectonic fault density, rainfall, and road density) were analyzed and mapped in geographic information system. An index expressing the correlation between each factor and landslides [called class landslide susceptibility index (CLSI)] was proposed in the study. While analyzing landslide distribution in a large area, point aggregation might be expected. To quantify the uncertainty caused by aggregation, class landslide aggregation index was proposed. To account for the importance of each of the factors in the landslide susceptibility assessment, some weights were calculated by means of analytic hierarchy process. We propose a weighted class landslide susceptibility model (WCLSM), obtained by the combination of CLSI values of each factor with the correspondent weight. WCLSM performance in the study area was evaluated comparing the results obtained by first modeling all landslides and then by performing a time partition. The model was run including only landslides that occurred before 2009 and then validated with respect to landslides that occurred after 2009. The prediction–rate curve shows that the WCLSM model provides a good prediction for the study area. Of the study area, 21.4 % shows very high and high susceptibility and includes the 87.7 % of the number of landslides that occurred after 2009.     

Diagnosis of landslide risk for individual buildings: insights from Prahova Subcarpathians,Romania

The aim of this study is to quantify the landslide risk for individual buildings using spatial data in a GIS environment. A landslide-prone area from Prahova Rivers’ Subcarpathian Valley was chosen because of its associated landslide hazards and its impact upon human settlements and activities. The bivariate landslide susceptibility index (LSI) was applied to calculate the spatial probability of landslides occurrence. The Landslide Susceptibility Index map was produced by numerically adding the weighted thematic maps for slope gradient and aspect, water table, soil texture, lithology, built environment and land use. Validation curves were obtained using the random-split strategy for two combinations of variables: (a) all seven variables and (b) three variables which showed highest individual success rates with respect to landslides occurrences (slope gradient, water table and land use). The principal pre-disposing factors were found to be slope steepness and groundwater table. Vulnerability was established as the degree of loss to individual buildings resulting from a potential damaging landslide with a given return period in an area. Risk was calculated by multiplying the spatial probability of landslides by the vulnerability for each building and summing up the losses for the selected return period.     

Estimation of Seismic Landslide Hazard in the Eastern Himalayan Syntaxis Region of Tibetan Plateau

The eastern Himalayan syntaxis is one of the most tectonically active and earthquake-prone regions on Earth where earthquake-induced geological disasters occur frequently and caused great damages. With the planning and construction of Sichuan-Tibet highway, Sichuan-Tibet railway and hydropower development on the Yarlung Zangbo River etc. in recent years, it is very important to evaluate the seismic landslide hazard of this region. In this paper, a seismic landslide hazard map is produced based on seismic geological background analysis and field investigation using Newmark method with 10% PGA exceedance probabilities in future 50 years by considering the influence of river erosion, active faults and seismic amplification for the first time. The results show that the areas prone to seismic landslides are distributed on steep slopes along the drainages and the glacier horns as well as ridges on the mountains. The seismic landslide hazard map produced in this study not only predicts the most prone seismic landslide areas in the future 50 years but also provides a reference for mitigation strategies to reduce the exposure of the new building and planning projects to seismic landslides.     

GIS-based statistical analysis of the spatial distribution of earthquake-induced landslides in the island of Lefkada, Ionian Islands, Greece

This is the first landslide inventory map in the island of Lefkada integrating satellite imagery and reports from field surveys. In particular, satellite imagery acquired before and after the 2003 earthquake were collected and interpreted with the results of the field survey that took place 1 week after this strong (Mw?=?6.3) event. The developed inventory map indicates that the density of landslides decreases from west to east. Furthermore, the spatial distribution of landslides was statistically analyzed in relation to the geology and topography for investigating their influence to landsliding. This was accomplished by overlaying these causal factors as thematic layers with landslide distribution data. Afterwards, weight values of each factor were calculated using the landslide index method and a landslide susceptibility map was developed. The susceptibility map indicates that the highest susceptibility class accounts for 38 % of the total landslide activity, while the three highest classes that cover the 10 % of the surface area, accounting for almost the 85 % of the active landslides. Our model was validated by applying the approaches of success and prediction rate to the dataset of landslides that was previously divided into two groups based on temporal criteria, estimation and validation group. The outcome of the validation dataset was that the highest susceptibility class concentrates 18 % of the total landslide activity. However, taking into account the frequency of landslides within the three highest susceptibility classes, more than 85 %, the model is characterized as reliable for a regional assessment of earthquake-induced landslides hazard.