基于Newmark模型的概率地震滑坡危险性分析方法研究

本文以天水地区为研究区,结合地震潜在震源区模型和Newmark位移预测方程,采用概率地震危险性分析方法,计算了该地区50年超越概率10%水平下的Newmark位移。同时,根据天水地区50年超越概率10%下的阿里亚斯烈度,并结合Newmark位移与阿里亚斯烈度的关系式,计算了天水地区在遭受50年超越概率10%下的阿里亚斯烈度影响时,潜在滑坡体产生的Newmark位移分布。通过比较上述两种方法得到的天水地区不同Newmark位移的分布特征,本文认为二者虽然存在较大差异,但其空间分布特征均能反映天水地区每个场点处的相对滑坡危险性。对滑坡危险性水平进行分区的结果显示,天水地区60%以上的区域具有高地震滑坡危险性,50%以上的区域具有甚高地震滑坡危险性。本文的研究结果可以作为天水地区地震危险性及风险评估的参考资料,也可以作为天水地区城市规划、土地使用规划、地震应急准备以及其它公共政策制定的参考资料。      

APPLICATION OF TWO SIMPLIFIED NEWMARK MODELS TO THE ASSESSMENT OF LANDSLIDES TRIGGERED BY THE 2008 WENCHUAN EARTHQUAKE

At present, with the wide application of the Newmark method, various Newmark empirical formulas with different ground motion parameters have been fitted by many researchers based on global strong-motion records. However, the existing study about the Wenchuan earthquake does not quantitatively evaluate the applicability of different Newmark models based on the actual landslides distribution. The aim of this paper is to present a comparison between observed landslides from the 2008 Wenchuan earthquake and predicted landslides using Newmark displacement method based on different ground motion parameters. The factor-of-safety map and critical acceleration(a_c)map in the study area are obtained by using the terrain data and geological data. The distribution of Arias intensity(I_a)and PGA in the study area is obtained by using the attenuation formulas of Arias intensity(I_a)and PGA, which is regressed by Wenchuan ground motion records. Based on the distribution of Arias intensity(I_a)and PGA parameters, we obtained the predicted locations of landslide using Newmark regression equations which are generated using global strong-motion records. The results shows that the assessment results can better reflect the macroscopic distribution characteristics of co-seismic landslides, most predicted landslide cells are distributed on the two sides of the Beichuan-Yingxiu Fault, especially the Pengguan complex rock mass in the hanging wall. The abilities to predict landslide occurrence of the two Newmark simplified models are different. On the whole, the evaluated result of simplified model based on parameter I_a is better than that based on PGA parameter. The GFC values obtained by the Newmark model of I_a and PGA parameters are 65.7% and 34.9%respectively. The evaluated result based on I_a can better reflect the macro distribution of coseismic landslides. The Ls_Pred value based on the Newmark model of parameter I_a is 26.5%, and the Ls_Pred value based on the Newmark model of PGA parameter is 10.3%. However the total area of predicted landslides accounts for 2.4% of the study area, which indicates that the predicted landslide cells are greater than the observed landslide cells. This reminds us that depending on the current input of shear strength and ground-motion parameters, we can only conduct landslide hazard assessment in macro areas, the ability to predict landslide can be improved using more accurate topographic data and input parameters.     

Real-time mapping of earthquake-induced landslides

Rigid sliding block analysis is a common analytical procedure used to predict the potential for earthquake-induced landslides for natural slopes. Currently, predictive models provide the expected level of displacement as a function of the characteristics of the slope (e.g., geometry, strength, yield acceleration) and the characteristics of earthquake shaking (e.g., peak ground acceleration, peak ground velocity). These predictive models are used for developing seismic landslide hazard maps which identify zones with risk of earthquake-induced landslides. Alternatively, these models can be combined with Shakemaps to generate “near-real-time” Slidemaps which could be used, among others, as a tool in disaster management. Shakemaps (a publicly available free service of the United States Geological Survey, USGS) provide near-real-time ground motion conditions during the time of an earthquake event. The ground motion parameters provided by a Shakemap are very useful for the development of Slidemaps. By providing ground motion parameters from an actual earthquake event, Shakemaps also serve as a tool to decouple the uncertainty of the ground motion in sliding displacements prediction. Campania region in Italy is studied for assessing the applicability of using Shakemaps for regional landslide-risk assessment. This region is selected based on the availability of soil shear strength parameters and the proximity to the 1980 Irpina (M _w  = 6.9) Earthquake.     

潜在地震滑坡危险区区划方法

不同地区地震活动的强度和频率是不同的.基于地震危险性分析的地震滑坡危险研究在综合了地震烈度、位置、复发时间等因素的基础上,考虑了地震动峰值加速度时空分布的特点,可以有效地应用于潜在地震滑坡危险区区划.以汶川地震灾区为研究对象,根据研究区的地质构造、地震活动特点等划分出灾区的潜在震源区,对该区进行地震危险性分析,并在此基础上采用综合指标法做出基于地震危险性分析的地震滑坡危险性区划.所得地震滑坡危险性区划按照滑坡危险程度分为高危险、较高危险、较低危险和低危险四级,表示未来一段时间内研究区在遭受一定超越概率水平的地震动作用下,不同地区地震滑坡发生的可能程度. 本文给出的地震滑坡危险性区划结果中,汶川地震滑坡崩塌较发育的汶川、北川、茂县等部分区域均处于高危险或较高危险区域;在对具有较高DEM精度的北川擂鼓镇地区所作的地震滑坡危险性区划中,汶川地震中实际发生的地震滑坡灾害与地震滑坡危险区划结果表现出较好的一致性.对区域范围而言,基于地震危险性分析的地震滑坡区划,可为初期阶段的土地规划使用及重大工程选址提供参考.     

基于修正岩土体强度参数的简化纽马克位移法地震滑坡危险性快速评估技术

地震滑坡危险性评估可为震后应急响应等提供科学的决策依据。纽马克位移法可不依赖同震滑坡编目快速评估同震滑坡危险性。工程岩体物理力学参数是该方法的核心参数之一,但其赋值过于单一,难以反映复杂地质背景下岩体强度的空间差异性。针对上述问题,本文在分析地震滑坡影响因子的基础上,选择距断层距离、高程和距水系距离作为影响岩体强度的评价指标并建立岩体强度评价模型,获得区域岩体强度修正系数,进而修正传统方法的临界加速度。结合震后的即时地震动峰值加速度,采用简化纽马克位移法计算边坡累积位移,开展地震滑坡危险性快速评估,并以汶川M_W7.9地震的地震滑坡危险性评估为例验证本文方法。结果表明,相对于传统方法,本文方法划分的地震滑坡危险区与同震滑坡分布更加一致。     

考虑动态临界加速度的地震边坡永久位移预测模型研究

在区域边坡地震危险性评价中主要采用永久位移预测模型进行地震边坡永久位移计算.永久位移预测模型以Newmark滑块理论为基础,通过大量实测地震时程记录统计拟合得出.针对Newmark理论中滑动面抗剪强度参数保持不变和已有位移预测模型的计算位移小于实测位移的问题,利用动态临界加速度理论,分别构建含有峰值加速度和阿里亚斯强度...     

Earthquake-reactivated landslide scenarios in Southern Italy based on spectral-matching input analysis

The Tyrrhenian portion of the Calabria region (southern Italy) is particularly prone to landslides as a consequence of intense morphodynamic processes. These processes affect the slopes that are composed of highly jointed metamorphic rock masses. Moreover, the frequent intense rainfalls and the up to Mw 7.0 regional earthquakes represent the main landslide triggering factors. An area of approximately $45\,\hbox {km}^{2}$ was selected as a test site in the context of a regional project aimed at reconstructing possible earthquake-reactivated landslide scenarios (i.e., referred to already existing landslide masses). An inventory map led to the identification of 175 landslides, including rock slides, earth slides and rock falls. Ground-motion scenarios based on a spectral-matching method were derived to evaluate the expected earthquake-induced displacements of the existing landslides. Naturally recorded acceleration time histories were selected from international ground-motion databases based on a similarity index and considered representative of the seismological features of the considered seismic sources (i.e., epicentral distance, magnitude, focal mechanism). Spectral attenuation was considered, according to well-established attenuation laws, to define the expected response spectrum at the outcropping bedrock corresponding to each existing landslide. Subsequently, the selected natural records were modified to guarantee spectral matching with the attenuated response spectra at each landslide site. The derived time histories were used to compute co-seismic displacements via the classic Newmark’s sliding-block method. Different scenarios of co-seismic landslide displacements or collapse were generated for different pore-water pressure hypotheses. The strongest $\hbox {Mw}>6$ seismic scenario (Messina Straits seismogenic source) indicated an exceedance probability of earthquake-induced co-seismic landslide collapse varying from 20 to 55 % with the increasing severity of the pore-water pressures. This probability corresponds to a percentage of co-seismic landslide displacements up to 40 % of the total inventoried landslides. The exceedance probability indicated that co-seismic landslide collapse drops below 20 % for $\hbox {Mw}<6$ seismic scenarios. In contrast, if a uniform probability is assumed for the seismic action occurrence, i.e., return periods of 475 and 2,475 years, the total percentage of landslide co-seismic displacements could be as high as 70 and 90 %, respectively, for the considered pore-water pressures.     

Rigorous landslide hazard zonation using Newmark's method and stochastic ground motion simulation

Researchers and practitioners in earthquake engineering have recognized geographic information systems (GIS) to be a significant tool in modeling spatial phenomenon related to hazard and risk. GIS, as an engineering tool, has been primarily used for its spatial data storing and presentation features. Models are often simplified to be more compatible with the light computational capabilities of many GIS. If not simplified, heavy computations are generally performed external to the GIS. A prototype vector-based GIS was developed that employs a rigorous approach to Newmark's displacement method for assessing earthquake triggered landslide hazards. The rigorous Newmark's analysis provides desirable flexibility by allowing input of actual ground motions. The prototype hazard GIS incorporates a popular shot filtered noise technique for generating artificial ground motions. The rigorous approach was compared to a popular simplified approach for computing Newmark displacements. Distribution of regional displacements was found to be similar with the simplified approach giving more and larger extreme displacements. The rigorous approach is suitable for large scales to model various seismic scenarios and their effect on seismically induced landslide potential.     

基于阿里亚斯烈度估值的概率性地震危险性分析——以四川丹棱县及其周缘为例

以编制《中国地震动参数区划图》时所构建的地震潜在震源区模型和地震活动性模型为基本输入,在四川省丹棱县及其周缘地区开展基于阿里亚斯烈度的概率地震危险性分析,计算50年超越概率10%的阿里亚斯烈度(Arias Intensity,Ia)。结果表明:丹棱县及其周缘绝大部分地区的Ia值都在0.11m/s以上,部分地区在0.32m/s甚至0.54m/s以上,具有较高的地震诱发滑坡风险,应当加强人员密集区的地震诱发滑坡危险性评估。根据不同Ia预测方程计算得到的Ia分布有较大差异,因此在计算Ia时应考虑多个Ia预测方程,对最终结果进行加权平均,以减小Ia结果的不确定性。同时还发现Ia值与1.0s的谱加速度具有很好的相关性,这也印证了Ia和1.0s谱加速度与砂土液化的相关性。     

GIS支持下的地震诱发滑坡危险区预测研究

为了满足对地震诱发滑坡危险区预测的不断增长的迫切要求，灾害评价成为帮助决策过程重要的基础工具之一。即使地震滑坡危险性各组份的评价很困难，但地理信息可辅助提出这种灾害制图的有关方法。描述了用于地理信息系统识别和定量计算不同地震滑坡危险区的技术方法，确定了地震烈度、地形坡度、岩土体类型和现存滑坡密度共4个因子参与的地震诱发滑坡危险性分析。在ARC／INFO DRID支持下，进行叠合分析，由此编制了云南省地震诱发滑坡危险区预测图。由地貌学家提出的地震诱发滑坡预测为规划和工程师提供了对区域规划和建筑工程有价值的技术方法。     

浙江省文泰震区地震滑坡危险性评估

地震滑坡灾害的震前预测与震后快速评估已成为减轻地震次生灾害的重要手段之一。本文使用简化Newmark模型,设定地震震级(M_S5.0),利用区域地质图、数字高程模型等基础数据,考虑地形对地震动的放大效应,对文泰震区潜在同震滑坡区域开展评估工作。研究表明,干燥与饱和状态下,设定地震作用下研究区内地震滑坡高危险区均主要分布在距设定震中15km以内的范围内,其分布与区内岩土体处于临界稳定状态的分布趋势相同。区内水库坝址与水库库体未受到潜在同震滑坡的影响,划定的重点关注区内位于潜在滑坡体下方的千秋门村、驮加村、高西村、杜山村、南峤村、包坑村、龙前村、新厂村以及各级公路易受到同震滑坡的影响,应提升重点关注区内承灾体的风险防范能力,尽可能减少潜在同震滑坡对区内生命财产安全造成的威胁。     

区域边坡地震危险性评价理论研究进展

地震滑坡会对自然环境和人民生活带来极大破坏,在大区域范围内对边坡遭受地震影响的程度进行评价判断,主要采用具有预测性质的潜在地震诱发滑坡危险性评估的方法。因此,深入研究边坡地震危险性分析的基本理论并作出符合实际的危险性分布图对工程建设和灾害防治具有较大意义。梳理并阐述具有预测性质的边坡地震危险性分析所需数据资料及各类数据的研究进展,将利用永久位移法进行边坡地震危险性评价所需数据归纳为三种：(1)边坡在地震影响下破坏程度的判定依据;(2)区域地震动参数如峰值加速度、阿里亚斯强度;(3)边坡坡体基本参数如黏聚力、摩擦角、重度、滑块厚度、坡角等。边坡地震危险性评价的准确程度与这三类数据的研究程度与准确性有关。文章对三类数据分别详细阐述各自的研究现状与成果应用,最后分析理论存在的部分问题以及以后的研究方向。     

Reliability analysis of three-dimensional dynamic slope stability and earthquake-induced permanent displacement

Different models were developed for evaluating the probabilistic three-dimensional (3D) stability analysis of earth slopes and embankments under earthquake loading using both the safety factor and the displacement criteria of slope failure. In the 3D analysis, the critical and total slope widths become two new and important parameters.The probabilistic models evaluate the probability of failure under seismic loading considering the different sources of uncertainties involved in the problem, i.e. uncertainties stemming from the discrepancies between laboratory-measured and in-situ values of shear strength parameters, randomness of earthquake occurrence, and earthquake-induced acceleration. The models also takes into consideration the spatial variabilities and correlations of soil properties.Five probabilistic models of earthquake-induced displacement were developed based on the non-exceedance of a limited value criterion. Moreover, a probabilistic model for dynamic slope stability analysis was developed based on 3D dynamic safety factor.These models are formulated and incorporated within a computer program (PTDDSSA).A sensitivity analysis was conducted on the different parameters involved in the developed models by applying those models to a well-known landslides (Selset landslide) under different levels of seismic hazard.The parametric study was conducted to evaluate the effect of different input parameters on the resulting critical failure width, 3D dynamic safety factor, earthquake-induced displacement and the probability of failure. Input parameters include: average values and coefficients of variations of water table, cohesion and angle of friction for effective stress analysis, scales of fluctuations in both distance and time, hypocentral distance, earthquake magnitude, earthquake strong shaking period, etc.The hypocentral distance and earthquake magnitude were found to have major influence on the earthquake-induced displacement, probability of failure (i.e. probability of allowable displacement exceedance), and dynamic 2D and 3D safety factors.     

Factors Controlling Seismic Susceptibility of the Sele Valley Slopes: The Case of the 1980 Irpinia Earthquake Re-Examined

Most of the documented slope failures triggered by the 1980 Irpinia earthquake (M_s 6.9) occurred in the upper Sele valley epicentral area (southern Italy). The early investigations revealed some puzzling characteristics of the slope failure distribution, i.e., (i) the higher landslide concentration on the valley slopes located farther away from the earthquake fault; (ii) the predominance of re-activations over first-time movements. The analyses of factors controlling the landslide concentrations indicates that the differences in hydrological setting and in slope were the two main causal factors whereas the seismic shaking, according to the radiation pattern modelling, could have been characterised by a relatively low rate of decrease across the valley. The aspect of the slopes did not play a significant role. The differences in groundwater conditions between the western and eastern valley sides were probably enhanced by the earthquake. In addition to the probable pore-water pressure rise, the seismic shaking caused large increases in the flow of springs draining the western aquifer, and this made the adjacent flysch slopes more prone to landsliding. Data from the available literature suggest that the effects of earthquake-induced groundwater release on seismic landslide distribution is especially important for normal-fault events. The Sele valley case also indicates that the slope of the pre-existing landslides is an important factor controlling their susceptibility to seismic re-activations.     

Slope instabilities triggered by the 11th May 2011 Lorca earthquake (Murcia,Spain): comparison to previous hazard assessments and proposition of a new hazard map and probability of failure equation

The Lorca Basin has been the object of recent research aimed at studying the phenomena of earthquake-induced landslides and its assessment in the frame of different seismic scenarios. However, it has not been until the 11th May 2011 Lorca earthquakes when it has been possible to conduct a systematic approach to the problem. In this paper we present an inventory of slope instabilities triggered by the Lorca earthquakes which comprises more than 100 cases, mainly rock and soil falls of small size (1–100  \(\hbox {m}^{3}\) ). The distribution of these instabilities is here compared to two different earthquake-triggered landslide hazard maps: one considering the occurrence of the most probable earthquake for a 475-years return period in the Lorca Basin \((\hbox {M}_{\mathrm{w}}=5.0)\) based on both low- and high-resolution digital elevation model (DEM); and a second one matching the occurrence of the \(\hbox {M}_{\mathrm{w}}=5.2\) 2011 Lorca earthquake, which was performed using the higher resolution DEM. The most frequent Newmark displacements related to the slope failures triggered by the 2011 Lorca earthquakes are lower than 2 cm in both the hazard scenarios considered. Additionally, the predicted Newmark displacements were correlated to the inventory of slope instabilities to develop a probability of failure equation. The fit seems to be very good since most of the mapped slope failures are located on the higher probability areas. The probability of slope failure in the Lorca Basin for a seismic event similar to the \(\hbox {M}_{\mathrm{w}}\) 5.2 2011 Lorca earthquake can be considered as very low (0–4 %).     

Time probabilistic evaluation of seismically induced landslide hazard in Irpinia (Southern Italy)

A recently proposed method, which incorporates the Newmark model to evaluate the earthquake-induced landslide hazard at regional scale, was applied to Irpinia, one of the most seismically active regions of Italy. The method adopts a probabilistic approach to calculate values of critical acceleration a_c representing the minimum strength required for a slope not to fail at a fixed probability level in a given time interval. Regional probabilistic hazard maps were generated for the two failure types most common in Irpinia (slump–earthflows and rock falls). The results suggest that quite moderate critical acceleration (0.05–0.08 g) could suffice to keep the slope failure probability low. However, the available data indicate that potential slide surfaces with a_c below these values could be common in Irpinia, where, perhaps in relation to particular geo-environmental conditions, a relative large number of marginally stable slopes might survive other destabilising actions and fail even on occasion of not particularly strong earthquake shaking.     

Assessing the rate of seismogravitational denudation of the relief of southeastern Altai: The Chagan-Uzun R. Basin

This paper considers the influence of earthquake-induced slope processes on relief changes. We used the parameters of Holocene maximum earthquake-induced landslides resulting from past great earthquakes in order to use these data for the Chagan-Uzun R. basin for deriving the rate of Holocene seismogravitational denudation for the highest southeastern part of the Russian Altai, which shows the highest seismic activity. It is shown that the influence of moderate earthquakes on this relief is negligibly small compared with that of larger earthquakes, while the landslides due to aftershocks provide a significant contribution in seismic regions where the aftershock process does not obey the Bath-Omori laws. Our numerical estimate for the rate of seismogravitational denudation obtained by statistical methods was corroborated by measuring the volumes of Holocene seismogravitational movements using the method of detailed profiling. It is shown that the volume of an earthquake-induced landslide is largely controlled by the size of the associated earthquake and by climate.     

Probabilistic PGA and Arias Intensity maps of Kyrgyzstan (Central Asia)

New probabilistic seismic hazard and Arias Intensity maps have beendeveloped for the territory of the Kyrgyz Republic and bordering regions.Data were mainly taken from the seismic catalogue of Kyrgyzstan and partlyfrom the world seismic catalogue. On the base of seismicity and activetectonics, seismic zones were outlined over the area. For these,Gutenberg-Richter laws were defined using mainly instrumental data, butregarding also historical events. Attenuation of acceleration inside the targetarea could not be determined experimentally since existing strong motiondata are insufficient. Therefore, empirical laws defined for other territories,principally Europe and China, were applied to the present hazardcomputations. Final maps were calculated with the SEISRISKIII programaccording to EUROCODE8 criteria, i.e. for a period of 50 years with90% probability of non-exceedance. For long-term prediction, 100 yearsmaps with 90% probability of non-exceedance have been developed. Theprocedure used for seismic hazard prediction in terms of PGA (PeakGround Acceleration) was also applied to Arias intensities in order to beable to define regional seismogenic landslide hazard maps.     

STUDY ON THE DISTRIBUTION PATTERN OF EARTHQUAKE- TRIGGERED LANDSLIDES BASED ON SEISMIC LANDSLIDE SUSCEPTIBILITY ANALYSIS: A CASE STUDY OF LANDSLIDES TRIGGERED BY THE MS6.5 LUDIAN EARTHQUAKE IN 2014

On August 3, 2014, an M_W6.5 earthquake occurred in Ludian County, Yunnan Province, which triggered significant landslides and caused serious ground damages and casualties. Compared with the existing events of earthquake-triggered landslides, the spatial distribution of co-seismic landslides during the Ludian earthquake showed a special pattern. The relationship between the co-seismic landslides and the epicenter or the known faults is not obvious, and the maximum landslide density doesn't appear in the area near the epicenter. Peak ground acceleration (PGA), which usually is used to judge the limit boundary of co-seismic landslide distribution, cannot explain this distribution pattern. Instead of correlating geological and topographic factors with the co-seismic landslide distribution pattern, this study focuses on analyzing the influence of seismic landslide susceptibility on the co-seismic distribution. Seismic landslide susceptibility comes from a calculation of critical acceleration values using a simplified Newmark block model analysis and represents slope stability under seismic loading. Both DEM (SRTM 90m)and geological map (1 ︰ 200^￣￣000)are used as inputs to calculate critical acceleration values. Results show that the most susceptible slopes with the smallest critical accelerations are generally concentrated along the banks of rivers. The stable slopes, which have the larger critical accelerations and are comparably stable, are in the places adjacent to the epicenter. Comparison of the distribution of slope stability and the real landslides triggered by the 2014 M_W6.1 Ludian earthquake shows a good spatial correlation, meaning seismic landslide susceptibility controls the co-seismic landslide distributions to a certain degree. Moreover, our study provides a plausible explanation on the special distribution pattern of Ludian earthquake triggered landslides. Also the paper discusses the advantages of using the seismic landslide susceptibility as a basic map, which will offer an additional tool that can be used to assist in post-disaster response activities as well as seismic landslides hazards zonation.     

基于高分卫星遥感影像的地震应急滑坡编目与分布特征探讨——以2017年8月8日九寨沟7.0级地震为例

地震应急是减轻地震灾害的重要途径之一。地震应急工作具有时间紧迫、事关重大的特点。2017年8月8日四川九寨沟M_S7.0级地震发生后,为快速、准确地提供地震引发的滑坡灾害分布,本研究基于震后第一天获取到的高分辨率遥感影像（高分二号卫星影像、北京二号卫星影像）,通过人工目视解译的方法初步建立了四川九寨沟地震滑坡编目。结果表明,该地震至少触发了622处同震滑坡,分布在沿使用影像边界框定的面积为3919km2的区域内。本研究还利用这个地震滑坡编目,统计了九寨沟地震滑坡数量和滑坡点密度（LND）与地形（坡度、坡向）、地震（地震烈度、震中距）等因素的关系。结果表明九寨沟地震滑坡多发生在坡度为20°—50°的区域内,滑坡的易发性随着坡度的增加而增加。受地震波传播方向的影响,E、SE向是地震滑坡较易发生的坡向。滑坡的易发程度和地震烈度呈正相关,即随着烈度的增大,滑坡易发性增大。滑坡易发性还随着震中距增加而降低,这是由于地震波能量随震中距的增加而衰减导致的。