中国地震烈度衰减的概率模型

本文在考虑地震烈度评定和预测的不确定性的基础上,提出建立烈度衰减概率模型的方法。根据我国东部和西部(大致以东经104为界)地区已发生的破坏性地震的等震线图,提供了两个地区地震烈度衰减的概率模型。最后,简要讨论了如何将这一模型应用于地震危险度评定的基本思想。      

Inversion-based data-driven time-space domain random noise attenuation method

Conventional time-space domain and frequency-space domain prediction filtering methods assume that seismic data consists of two parts, signal and random noise. That is, the so-called additive noise model. However, when estimating random noise, it is assumed that random noise can be predicted from the seismic data by convolving with a prediction error filter. That is, the source-noise model. Model inconsistencies, before and after denoising, compromise the noise attenuation and signal-preservation performances of prediction filtering methods. Therefore, this study presents an inversion-based time-space domain random noise attenuation method to overcome the model inconsistencies. In this method, a prediction error filter (PEF), is first estimated from seismic data; the filter characterizes the predictability of the seismic data and adaptively describes the seismic data’s space structure. After calculating PEF, it can be applied as a regularized constraint in the inversion process for seismic signal from noisy data. Unlike conventional random noise attenuation methods, the proposed method solves a seismic data inversion problem using regularization constraint; this overcomes the model inconsistency of the prediction filtering method. The proposed method was tested on both synthetic and real seismic data, and results from the prediction filtering method and the proposed method are compared. The testing demonstrated that the proposed method suppresses noise effectively and provides better signal-preservation performance.     

青海玛多7.4级地震烈度快速评估

利用青海地震烈度衰减模型和基于断层最短距离的地震动参数衰减模型对青海玛多7.4级地震的烈度分布进行快速评估,将评估结果与实际烈度图对比分析.研究结果表明:两种模型可以在地震正式速报后1分钟内得到地震烈度分布的快速评估结果,时效性强,可以在最短的时间内为应急指挥和救援工作提供参考;青海地震烈度衰减模型计算结果比较理想化,...     

Earthquake-induced landslide displacement attenuation models and application in probabilistic seismic landslide displacement analysis

A landslide displacement (DLL) attenuation model has been developed using spectral intensity and a ratio of critical acceleration coefficient to ground acceleration coefficient. In the development of the model,a New Zealand earthquake record data set with magnitudes ranging from 5.0 to 7.2 within a source distance of 175 km is used. The model can be used to carry out deterministic landslide displacement analysis,and readily extended to carry out probabilistic seismic landslide displacement analysis. DLL attenuation models have also been developed by using earthquake source terms,such as magnitude and source distance,that account for the effects of earthquake faulttype,source type,and site conditions. Sensitivity analyses show that the predicted DLL values from the new models are close to those from the Romeo model that was developed from an Italian earthquake record data set. The proposed models are also applied to an analysis of landslide displacements in the Wenchuan earthquake,and a comparison between the predicted and the observed results shows that the proposed models are reliable,and can be confidently used in mapping landslide potential.     

2021年5月21日云南漾濞6.4级地震灾害损失快速评估结果分析

本文针对2021年5月21日云南漾濞6.4级地震,选取不同的地震烈度衰减关系模型,对各模型地震影响场评估结果与发布的地震烈度图进行对比分析,并对地震影响范围不确定性进行研究。选取多种死亡人数评估模型,分别计算各模型在不同地震影响场下的死亡人数、人口分布数量,探讨各地震影响场模型下的人口分布特征及影响人员死亡的主要因素。通过对比分析可知,导致此次地震灾害损失评估结果与真实地震现场结果不同的主要原因是地震影响场分布、人口分布、房屋建筑抗震能力偏差、地形地貌、次生灾害等多种因素不同。研究结果表明,有效提高地震灾害损失快速评估精确性的途径为提高地震影响场评估精度,提高人口分布、房屋建筑等数据空间分布评估精度,后期专家检验等。     

地震动衰减模型在玛多Mw7.3级地震影响场快速评估中的应用

2021年5月22日青海省玛多县发生M_w7.3级地震。震后,根据初步估计的断层走向和破裂长度,基于YU15地震动衰减模型和三种NGA-West2（Next Generation Attenuation-West2）地震动衰减模型快速产出地震区震动图及理论烈度图。在获得强震记录和地表破裂长度信息后,对预测结果进行修正。通过比较理论烈度与调查烈度,并结合震动图分布形态以及衰减模型在2016年新疆呼图壁M_w6.0地震中的应用情况对四种地震动衰减模型的适用性进行了分析。结果表明:在台网稀疏地区,基于地震动衰减模型可在震后快速获得地震动分布,并产出具有应用价值的地震影响场;NGA-West2模型在断层破裂较长的大震中表现优于YU15模型,而在中强地震中后者适用性更强;近实时强震动记录可用来检验模型的适用性并对预测结果进行修正;断层破裂尺度、震源机制和破裂过程等信息的准确估计可有效提高地震影响场预测精度。     

关于地震烈度衰减模型的系统偏差

本文讨论在地震危险性分析中目前使用广泛的地震烈度衰减模型,即Cornell建立的点源模型和洪华生(A.H-S.Ang)等建立的断层-破裂模型。与实际烈度资料的比较表明,在等震线狭长的区域,点源模型有系统偏差,通常这种偏差不算严重。为了克服上述偏差,断层-破裂模型引进了破裂长度的概念,但其衰减与实际资料出入较大。与实际地震等震线的比较和计算表明,对大地震沿断层或等震线长轴方向,断层破裂模型总是高估高烈度区而低估低烈度区。此外,本文分析了断层破裂模型产生系统偏差的原因,讨论了点椭圆模型的合理性,并用地震危险性分析实例做了说明。      

Macroseismic attenuation in the Campanian area, southern Italy

The main objective of the present study is to evaluate seismic attenuation relationships for the Campanian area (southern Italy) using the felt intensity report data obtained from comprehensive historical databases (DOM 4.1). We focused our attention on the Campania region because it is characterized by a high seismic hazard and risk, particularly in the Naples area and its suburbs. In order to derive an attenuation relationship for the area, we fitted the observed data by using several functions. We found that a linear plus logarithmic model gives the best fits for the data in the Campanian region. Most of the attenuation relationships proposed up to now for the Italian Peninsula have an isotropic behavior and do not always properly describe the macroseismic attenuation. Therefore, in order to check the possible dependence of the attenuation on the azimuth of the seismic rays, we divided our data set in octants and performed for each of them the same analysis we carried out for the whole data set. The obtained results differ from octant to octant and the differences, besides being associated with the source effects, could be interpreted as probably due to the existence of lateral variations in the lithological and physical features of the crust at different depths, which could affect the patterns of attenuation.     

Attenuation relation of Arias intensity for Zagros Mountains region (Iran)

Arias intensity is considered as a shaking parameter suitable for characterizing earthquake impact on ground stability. Within the framework of a study aimed at providing tools for the assessment of hazards related to earthquake-induced slope failures, Arias intensity attenuation relations were determined for the Zagros Mountains region, an active tectonic belt elongated NW–SE in the western and south-western part of Iran. The calculation of relation coefficients was based on strong-motion data of earthquakes located in the Zagros area and recorded by Iranian stations managed by the Building and Housing Research Center of Iran (BHRC). Five models of attenuation relation were considered and their coefficients were estimated through a least-square regression analysis. The relations obtained were then applied to a data sample different from that used for regression and the root mean square (RMS) of residuals was examined in order to compare effectiveness of different relations in probabilistic estimates. Furthermore a comparison made with attenuation relations obtained for Alborz and the central part of Iran showed significant differences possibly related to structural differences.     

New Seismic Zoning Map of Chin a and Its Application in Earthquake Disaster Mitigation

The New Seismic Zoning Map of China was prepared from 1987 to 1990 and officially promulgated in 1991.In comparison with the previous two seismic zoning maps prepared in 1957 and 1977,some new methods were applied to upgrade the method currently used for seismic hazard analysis.First,a probabilistic method was used instead of the deterministic analysis was used for previous mapping.Second,by taking advantages of the long history of historical seismic data in China,the nonhomogeneity of seismicity both in space and in time has been fully considered and hence the over-and/or underestimation of seismic hazard could be avoided.Third,the results of middle-and long-term earthquake prediction based on tectonic evidence have been incorporated into seismic hazard analysis.In addition,the attenuation laws for both intensity and peak acceleration of strong motion as the mapping parameters are also presented.Finally,an evaluation of the New Seismic Map and its effect on engineering application,such as aseismic desi     

An empirical spectral ground-motion model for Iran

A new ground-motion prediction equation for 5%-damped horizontal spectral acceleration applicable to Iran is presented. On the basis of analysis of variance (ANOVA), selected West-Eurasian records are added to an existing dataset of Iranian accelerometric data to yield a ground-motion prediction equation applicable at wider ranges of magnitude and distance. The advantages of using this model rather than those proposed previously for Iran are discussed by considering the distribution of residuals against the explanatory variables, magnitude and distance. The applicability of the proposed model, as well as those of several other models developed for shallow crustal environments, is also investigated by means of statistical tools. The results reveal the overall suitability of the new model as well as the validity of models developed using mainly Eurasian strong-motion records. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users.     

南北地震带地区地震烈度衰减关系研究

收集了我国南北地震带地区1970年~2012年的85次地震事件的烈度等震线资料,建立了震中烈度与震级之间的经验关系,利用长轴和短轴椭圆模型拟合得到了该地区的地震烈度衰减关系。同时还搜集整理了近年来基于我国各地区的地震资料,研究得到的地震烈度衰减关系,并将其中与南北地震带存在地域重合的研究结果与本文结果进行对比分析。该研究结果能够较好地反映南北地震带地区的地震烈度衰减规律,对该地区的地震灾害快速评估具有一定的参考意义。     

预测多次波的逆散射级数方法与SRME方法及比较

逆散射级数方法和SRME方法是典型的两类基于波动方程的无需地下结构信息的自由表面多次波预测方法,本文详细讨论了这两类方法的基本思想、方法原理和实现思路,并对2维逆散射级数方法进行了降维简化,推导给出了1.5维逆散射级数衰减自由表面多次波的方法,减少了计算量,并降低了对3维规则观测系统的要求.利用经典SMAART模型数据测试比较了2维逆散射级数方法、1.5维逆散射级数方法及2维SRME方法预测自由表面多次波的效果及优缺点,比较了三者在方法实现、并行计算效率及对数据要求上的异同,并结合分析比较结果,给出了实际应用中方法选择的建议.     

Strong ground-motion relations for Mexican interplate earthquakes

We derive strong ground-motion relations for horizontal components of pseudo-acceleration response spectra from Mexican interplate earthquakes at rock sites (NEHRP B class) in the forearc region. The functional form is obtained from the analytical solution of a circular finite-source model. For the regression analysis we use a recently proposed multivariate Bayesian technique. The resulting model has similar accuracy as those models derived from regional and worldwide subduction-zone databases. However, there are significant differences in the estimations computed from our model and other models. First, our results reveal that attenuation in Mexico tends to be stronger than that of worldwide relations, especially for large events. Second, our model predicts ground motions for large earthquakes at close distances to the source that are considerably larger than the estimations of global models. Lack of data in this range makes it difficult to identify the most appropriate model for this scenario. Nevertheless, according to the available data at the city of Acapulco, our model seems to estimate seismic hazard more adequately than the other models. These new relations may be useful in computing seismic hazard for the Mexican forearc region, where no similar equations had been yet proposed.     

河西重点监视区烈度预测的研究

在前人研究的基础上，建立了适合于重点监视区烈度预测的概率方法的数学模型，并对该模型中影响预测结果的因素进行了灵敏度分析。发现预测的不确定性主要来自烈度衰减规律。根据概率分析中对正态分布函数不确定性校正的原理，提出了适于重点监视区内烈度预测中不确定性的校正方法，并将该方法应用于河西重点监视区。在确定震级－频度关系的过程中，考虑了地震活动性周期及震级间隔两方面的因素，使得求出的震级－频度关系更为精确     

Spatially constrained attenuation compensation in the mixed domain

Seismic attenuation compensation is a spectrum-broadening technique for enhancing the resolution of non-stationary seismic data. The single-trace attenuation compensation algorithms ignore the prior information that the seismic reflection events are generally continuous along seismic traces, thus, the compensated result may have poor spatial continuity and low signal-to-noise ratio. To address this problem, we extend the single-trace approaches to the multi-trace algorithms and furthermore propose a multi-trace attenuation compensation with a spatial constraint. The frequency-space prediction filters are the key to construct this spatial regularization. We test the effectiveness of the proposed spatially constrained attenuation compensation algorithm by applying both synthetic and field data examples. Synthetic data tests indicate that the proposed multi-trace attenuation compensation approach can provide a better compensated result than single-trace attenuation compensation algorithm in terms of suppressing noise amplification and guaranteeing structural continuities. Field data applications further confirm its stability and practicality to improve seismic resolution.     

Estimates of site seismicity rates using ill-defined macroseismic data

A new approach to the problem of site seismic hazard analysis is proposed, based on intensity data affected by uncertainties. This approach takes into account the ordinal and discrete character of intensities, trying to avoid misleading results due to the assumption that intensity can be treated as a real number (continuous distribution estimators, attenuation relationships, etc.). The proposed formulation is based on the use of a distribution function describing, for each earthquake, the probability that site seismic effects can be described by each possible intensity value. In order to obtain site hazard estimates where local data are lacking, the dependence of this distribution function with the distance from the macroseismic epicenter and with epicentral intensity is examined. A methodology has been developed for the purpose of combining such probabilities and estimating site seismicity rates which takes into account the effect of uncertainties involved in this kind of analysis. An application of this approach is described and discussed.     

A PRELIMINARY STUDY ON THE METHOD OF SEISMIC INTENSITY ASSESSMENT BASED ON RESIDENTIAL BUILDING DATA AND HIGH RESOLUTION REMOTE SENSING IMAGES

After destructive earthquakes, the assessment result of seismic intensity is an important decision-making basis for emergency rescue, recovery and reconstruction. This job requires higher timeliness by government and society. Because remote sensing technology is not affected by the terrible traffic conditions on the ground after the earthquake, large-scale seismic damage information in the earthquake area can be collected in a short time by the remote sensing image. The remote sensing technique plays a more and more important role in rapid acquisition of seismic damage information, emergency rescue decision-making, seismic intensity assessment and other work. On the basis of previous studies, this paper proposes a new method to assess seismic intensity by using remote sensing image, i.e. to interpret the building collapse rate of a residential quarter after an earthquake by high-resolution remote sensing images. If there already are detailed building data and building structure vulnerability matrix data of a residential area, we can calculate the building collapse rate under any intensity values in this residential area by using the theory of earthquake damage prediction. Assuming that the building collapse rate interpreted by remote sensing is equal to the building collapse rate predicted by using the existing data, it will be easy to calculate the actual seismic intensity of the residential area in this earthquake event. Based on this idea, according to the relevant standard specifications issued by China Earthquake Administration, this paper puts forward some functional models, such as the calculation model of building collapse rate based on remote sensing, the data matrix model of residential building structure, the prediction function matrix model of residential building collapse rate and the prediction model of residential building collapse rate. A formula for calculating seismic intensity by using remote sensing interpretation of collapse rate is also proposed. To test and verify the proposed method, this paper takes two neighboring blocks of Jiegu Town after the Yushu M7.1 earthquake in Qinghai Province as an example. The building structure matrix of the study block was constructed by using pre-earthquake 0.6m resolution satellite remote sensing image(QuickBird, acquired on November 6, 2004), post-earthquake 0.2m aerial remote sensing image(acquired by National Bureau of Surveying and Mapping, April 15, 2010) and some field investigation data. The building collapse rate in the two blocks was calculated by using the interpretation results of seismic damage from the Remote Sensing Technology Coordinating Group of China Seismological Bureau. The seismic damage matrix of building structures in Yushu area is constructed by using the abundant scientific data of the scientific investigation team of the project “Comprehensive Scientific Investigation of the Yushu M7.1 Earthquake in Qinghai Province” of China Seismological Bureau. On this basis, the collapse rate prediction function of different structures in Yushu area is constructed. According to the prediction function of collapse rate and the building structure matrix of the two blocks, the building collapse rate under different intensity values is predicted, and the curve of intensity-collapse rate function is drawn. By comparing the building collapse rate interpreted by remote sensing and the intensity-collapse rate function curve of this two blocks, the seismic intensity of both blocks are calculated to be the same value: Ⅸ degree, which is consistent with the results of the field scientific investigation of the earthquake. The validation shows that the method proposed in this paper can effectively avoid the influence caused by the difference of seismic performance of buildings and accurately evaluate seismic intensity when using remote sensing technique. The method has certain application value for earthquake emergency work.     

Seismic Intensity Zoning Map of China (1990) and Its Explanation

The Seismic Intensity Zoning Map of China(1990)was based on the probabilistic method of seismic hazard analysis.In compiling the map,the characteristics of inhomogeneity of earthquake distribution both in space and time in China are considered sufficiently,and some necessary modifications in the model of seismic hazard analysis are carried out.Based on the analysis of the seismic activity and seismotectonic environment,26 seismic provinces are divided first as the statistical elements of the seismicity analysis; the seismic potential source areas are then divided in the seismic provinces.The 733 potential source areas with various upper limit magnitudes have been divided in the country.According to the reliable time domain of earthquake data with various magnitude intervals,the b values in magnitude-frequency relationship are calculated in the seismic provinces.According to the analysis of the inhomogeneity of seismicity distribution both in space and time,the annual average occurrence rates of the eart     

Regional adjustment factors for three NGA-West2 ground-motion prediction equations to be applicable in northern Iran

Ground-motion prediction equations (GMPEs) are essential tools in seismic hazard studies to estimate ground motions generated by potential seismic sources. Global GMPEs which are based on well-compiled global strong-motion databanks, have certain advantages over local GMPEs, including more sophisticated parameters in terms of distance, faulting style, and site classification but cannot guarantee the local/region-specific propagation characteristics of shear wave (e.g., geometric spreading behavior, quality factor) for different seismic regions at larger distances (beyond about 80 km). Here, strong-motion records of northern Iran have been used to estimate the propagation characteristics of shear wave and determine the region-specific adjustment parameters for three of the NGA-West2 GMPEs to be applicable in northern Iran. The dataset consists of 260 three-component records from 28 earthquakes, recorded at 139 stations, with moment magnitudes between 4.9 and 7.4, horizontal distance to the surface projection of the rupture (R _JB) less than 200 km, and average shear-wave velocity over the top 30 m of the subsurface (V _S30) between 155 and 1500 m/s. The paper also presents the ranking results for three of the NGA-West2 GMPEs against strong motions recorded in northern Iran, before and after adjustment for region-dependent attenuation characteristics. The ranking is based on the likelihood and log-likelihood methods (LH and LLH) proposed by Scherbaum et al. (Bull Seismol Soc Am 94: 2164–2185, 2004, Bull Seismol Soc Am 99, 3234–3247, 2009, respectively), the Nash–Sutcliffe model efficiency coefficient (Nash and Sutcliffe, J Hydrol 10:282–290, 1970), and the EDR method of Kale and Akkar (Bull Seismol Soc Am 103:1069–1084, 2012). The best-fitting models over the whole frequency range are the ASK14 and BSSA14 models. Taking into account that the models’ performances were boosted after applying the adjustment factors, at least moderate regional variation of ground motions is highlighted. The regional adjustment based on the Iranian database reveals an upward trend (indicated as high Q factor) for the selected database. Further investigation to determine adjustment factors based on a much richer database of the Iranian strong-motion records is of utmost important for seismic hazard and risk analysis studies in northern Iran, containing major cities including the capital city of Tehran.