Probabilistic modeling of seismically triggered landslides using Monte Carlo simulations

The 2004 Chuetsu earthquakes of Niigata (Japan) triggered numerous landslides, and the most widespread types of landslides were highly disrupted, relatively shallow slides and soil (debris) flows. This paper presented a method to evaluate slope instability using Newmark displacement on a pixel-by-pixel basis in a given area. The proposed method was able to integrate Newmark displacement modeling and Monte Carlo simulations within geographical information systems. In the modeling, an empirical attenuation relationship was utilized to calculate Arias intensity over this study area, and the variability of geotechnical parameters was taken into account to calculate coseismic landslide displacement. Before deriving the displacement from related inputs, the Monte Carlo simulations ran 1,500 times and generated 1,500 displacement values for each grid cell, and then means and standard deviations of displacement were calculated and probabilistic distributions can be obtained. Finally, given 10 cm as a threshold value of displacement, estimated probabilities of displacement exceeding 10 cm were shown as a map of seismic landslide hazards. The resulting hazard map was classified into four categories from very low to high level.     

基于贡献权重迭加法的滑坡风险区划

贡献权重迭加法是对滑坡本底因子和承灾体因子在滑坡发育中的贡献率进行统计后,通过贡献率均值化、规一化处理,利用权重转换模型计算出每一个因子内部的权重——自权重W和因子相互之间的权重——互权重W＇。将滑坡因子、承灾体因子和两者的自权重、互权重分别相乘迭加,得到滑坡危险度和易损度区划结果,再将2项结果相乘得出滑坡风险区划。     

基于蒙特卡罗法的多级黄土滑坡可靠性分析

滑坡稳定性分析目前工程中一般采用极限平衡算法。为了确定滑坡的稳定性系数,在计算过程中要将岩土体、计算模型简化,忽略一些影响滑坡稳定的次要因素,把影响滑坡稳定的各种主要因素作为确定参数。由于滑坡岩土体本身的非均质性及其参数确定具有随机性和变异性,用唯一的稳定系数不能客观地反映整个滑坡的真实安全程度。以陕西省夏呀河滑坡为例,通过野外详细调查、工程勘探和室内土工试验,初步确定该滑坡存在四级滑动面,难以用简单的单一滑动面计算其稳定性。因此,先采用推力传递系数法对滑坡体上的四级滑坡分别进行滑坡稳定性计算,得到该滑坡体上的四级滑坡稳定系数。再采用蒙特卡罗法对该四级滑坡进行可靠性分析,得到其可靠度,定量地表达夏呀河四级滑坡的安全程度。同时对比分析夏呀河四级滑坡的稳定系数和失稳概率,综合评价该滑坡体上的四级滑坡的稳定性及其风险概率,为滑坡的工程治理及预测预警、灾害危险性评价提供依据。     

分布式光纤传感技术在滑坡监测中的应用

将分布式光纤传感技术引入滑坡监测,既可得到整个滑坡体的概要特征,又能提高监测效率。FBG与BOTDR是两种最具代表性的分布式光纤传感技术,FBG通过测量其反射光波长的变化获得应变或温度值。BOTDR通过测定后向布里渊散射光的频移实现分布式温度、应变测量。FBG传感器灵敏度高,但只能实现离散点的准分布式测量,BOTDR可实现分布式、长距离、不间断测量,但其空间分辨率不高。笔者提出FBG与BOTDR联合监测滑坡的方案。在巫山残联滑坡,在整个滑坡体上铺设监测光纤,利用BOTDR获得整个滑坡体的概要信息;在滑坡体变形的关键部位———变形缝安装FBG传感器,获得某些关键部位的应变信息,从而实现由点到线再到面的监测,获得滑坡体较完整的应变信息。     

地震学定量预报指标的一种提取方法——对比筛选法

鉴于目前地震学综合定量预报指标的缺乏和预报工作的急需,尝试使用“对比筛选法(简称CSM方法)”进行地震学定量预报指标的提取试验。较之以往作法的进展在于:1.同时使用“有震”和“无震”两类样本对比筛选;2.对不同地区的地震学参数进行了归一化处理。这样做的显著优点是:1.可以较有效地提取“有震异常”和“正常变化”指标;2.提取的异常和预报指标具有定量化和普适性特点。 试验研究使用大华北地震区资料,研究对象取中强地震。经内符和外推检验,证明该方法提取的异常和预报指标有效性和实用性较高。     

地震综合效应场函数及其在地震预报中的应用──（一）地震综合效应场函数及其计算方法

在系统地分析了目前各种测震学地震预报方法科学思路的基础上,认为测震学地震预报方法基本上可以分为两大类。一类是以已经发生的一些地震作为未来可能发生的地震的“因”,即由于已经发生的地震对区域应力场的影响,导致未来发生较强地震。这一类包括的预报方法较多,如空区、条带、ｂ值、地震迁移、相关地震等等及其由此衍生出来的各种方法。另一类是把已经发生的一些地震作为区域应力场增强的“果”,即已经发生的地震是区域应力场增强过程中的一种反映,而未来地震不一定是已经发生的地震所导致的结果。这一类包括“地震窗口”、小震群活动等方法。针对第一类方法,各种预报方法都是力图从地震三要素中提取未来地震的信息,而具体作法又都是利用地震三要素这个多维空间的某个剖面。为了从地震活动诸要素的多维空间提取综合信息,我们对每个地震加入了破裂面方位,构成了地震第四要素,并依据地震４要素建立了地震综合效应场函数。地震综合效应场函数概括了多种测震学地震预报方法的科学思路和预报经验,从而可以形成测震学的综合预报方法。     

Prediction of the occurrence of Related Strong Earthquakes in Italy

In the seismic flow it is often observed that a Strong Earthquake (SE), is followed by Related Strong Earthquakes (RSEs), which occur near the epicenter of the SE with origin time rather close to the origin time of the SE. The algorithm for the prediction of the occurrence of a RSE has been developed and applied for the first time to the seismicity data of the California-Nevada region and has been successfully tested in several regions of the world, the statistical significance of the result being 97%. To date it has been possible to make five successful forward predictions, with no false alarms or failures to predict.The algorithm is applied here to the Italian territory, where the occurrence of RSEs is a particularly rare phenomenon. Our results show that the standard algorithm is successfully directly applicable without adjustment of the parameters. Eleven SEs are considered. Of them, three are followed by a RSE, as predicted by the algorithm, eight SEs are not followed by a RSE, and the algorithm predicts this behaviour for seven of them, giving rise to only one false alarm. Since, in Italy, often the series of strong earthquakes is relatively short, the algorithm has been extended to handle such a situation. The result of this experiment indicates that it is possible to attempt to test a SE, for the occurrence of a RSE, soon after the occurrence of the SE itself, performing timely preliminary recognition on reduced data sets. This fact, the high confidence level of the retrospective analysis, and the first successful forward predictions, made in different parts of the World, indicates that, even if additional tests are desirable, the algorithm can already be considered for routine application to Civil Defence.     

首都减灾圈的思考与建议

首都减灾圈,系首都减轻自然灾害预测防治圈。1991年12月20—21日在北京召开了《首都圈自然灾害及其减灾对策研讨会》。本文根据这次会议所提供的材料,在从整体上实现减灾的思想指导下,就首都减灾圈的“成灾背景”、“首都减灾圈的组成”、“首都减灾圈的灾害预测与防治状况”、“首都减灾圈的灾害关联性分析”、“首都减灾圈的减灾实效预估”和“首都减灾圈的减灾对策与实施”等六个带有共同性的问题,进行了讨论。供制定首都圈减灾方案时参考。     

新疆地震预报研究

新疆的地震预报是1970年开始的。二十多年来在“边观测、边研究、边预报”、“多路探索、多兵种联合作战”、“走综合预报之路”的方针指导下取得了长足的进展。建设了43个地震台站,投入189套仪器,建成了遍布全疆的地震无线通讯网,开展了历史地震调查和地震烈度区划工作,建立健全了地震会商预报制度,进行了多方面的地震预报研究。实际地震预报统计分析表明,扣除自然发震概率之后,趋势预报的成功率约0.3,短临预报的成功率约0.1。在前兆台网控制范围内取得了一些震例,说明地震确实是有前兆的,但又是非常复杂的。地震预报作为科学难题还有漫长的路要走。作者简要地讨论了地震观测的间接性和地震异常的离散性、难以区别的地震异常和地壳变动异常、建立在复杂现象基础上的地震前兆以及地震的混沌性对地震预报的影响等科学问题。     

A Block-Based Theoretical Model Suited to Gravitational Sliding

An original theoretical model has been devised to simulate mass flow over hill slopes due to gravitational sliding. The sliding mass is discretized into a sequence of contiguous blocks which are subjected to gravitational forces, to bottom friction and to surface resistance stresses that are generally negligible for subaerial flows, but are relevant for submarine slides. The blocks interact with each other while sliding down the hill flanks because of internal forces that dissipate mechanical energy and produce a momentum exchange between the individual blocks, yet conserving the total momentum of the mass. Internal forces are expressed in terms of interaction coefficients depending on the instantaneous distance between the block centers of mass, which is a measure of the deformation experienced by the blocks: the functional dependence includes three parameters, namely the interaction intensity ¯, the deformability parameter and the shape parameter , by means of which a wide range of interaction types can be fully accounted for. The time integration is performed numerically by solving the equations for the block velocities and positions at any time ti by means of the block accelerations at the previous time t_i-1, and by subsequently updating the block accelerations, which allows to proceed iteratively to the following times. The model has been tested against laboratory results available from literature and by means of several numerical experiments involving a simplified geometry both for the sliding body and the basal surface, with the purpose of clarifying the influence of the model parameters on the slide dynamics. The model improves the performance of the existing kinematic models for slides, moreover preserving an equivalent numerical simplicity. Future applications and possible improvements of this model are suggested.