Effects of site geology on seismic ground motion: early history

This is the 6^th contribution in the series of Historical Notes on seminal concepts in earthquake engineering and structural dynamics. It documents the origins and early developments (from the 1880s through 1992) of the effects of site geology on seismic ground motion. Copyright © 2014 John Wiley & Sons, Ltd.     

利用区域地震记录测定地震辐射能量

地震辐射能量作为描述地震大小的物理量,可为地震应急和灾害评估提供重要参考.现有快速测定地震辐射能量的方法多使用远震记录（35°≤Δ < 80°）,受到几何扩展和频率相关衰减校正方法的限制,利用区域地震记录（5°≤Δ < 35°）测定地震辐射能量的难度较大.因此,本文发展了一种利用区域地震记录快速测定地震辐射能量的方法.研究结果表明：（1）利用该方法可以测定5.0级以上地震的辐射能量,弥补了由于低信噪比和台站分布影响,导致利用远震记录只能稳定测定6.0级以上地震辐射能量的不足;（2）将该方法应用于2009—2021年发生在中国大陆的66次M_W>5.0地震,结果显示74%的单台区域能量震级与远震能量震级的偏差在±0.3以内;对于44次M_W≥5.5地震,区域结果与远震结果基本一致,86%的事件区域能量震级与远震能量震级的差在±0.2之间;（3）结合地震矩资料,得到中国大陆地区地震的能矩比范围为5.2×10^-6~8.1×10^-5,平均能矩比为2.4×10^-5;走滑型地震的平均能矩比高于倾滑型地震,能矩比整体分布存在区域特征,西域地块地震的平均能矩比高于青藏地块东部地震,青藏地块西部地震的平均能矩比最低,同一区域内的地震能矩比也存在差别,这表明不同断层上发生的地震能量释放过程存在明显差异;（4）该方法具有广泛的应用前景,可以在地震台网的日常工作中测定地震辐射能量,为科学研究、地震应急和地震灾害评估提供新的地震参数.

     

大地震与全球地震台网

引言 地震学家经常开玩笑说,＂阻止地震的最好办法就是布设地震台站。＂为了记录地震信号或开展有针对性的有关地球结构的研究就需要布设地震仪,而地震的发生往往是变幻莫测的,在时间和空间上都绝对是非均匀的,所以安装地震仪的辛苦有时会徒劳无功。     

中国历史地震烈度表研究

在比较分析以往烈度表的基础上,着重增加了社会反响标志;对Ⅹ-Ⅻ度的房屋建筑物和地表现象标志进行了调整与补充,完善了作为12阶烈度表相应的《中国历史地震烈度表》.文中对烈度表的各项标志作了简要说明,并列举了国内外10次历史地震事件的评定实例.本文提出的历史地震烈度表,保持了以往烈度表的适用性与一致性.     

抗震设计中多遇地震的问题和讨论

指出作为截面抗震验算标准的“小震”混淆了两种明显不同的概念:其一是对应于50年超越概率63％的地震动参数(以下简称“多遇小震”);其二是50年超越概率10％的地震动参数经过折算得到的参数(以下简称“折算小震”).上述混淆导致抗震设计中实际上执行着不同的标准.文中讨论了中、小地震动参数之间关系的复杂性,包括从地震危险性分析得到的基岩面中、小地震峰值加速度的离散性,一致概率谱的差异,以及中、小地震由于上层反应非线性影响导致的场地相关谱的系统性差异.指出重大工程的抗震设计以“多遇小震”取代“折算小震”进行截面验算可能导致反应谱特征周期Tg明显低估的问题,并提出了相应建议.     

A ground motion prediction equation for JMA instrumental seismic intensity for shallow crustal earthquakes in active tectonic regimes

The JMA (Japan Meteorological Agency) seismic intensity scale has been used in Japan as a measure of earthquake ground shaking effects since 1949. It has traditionally been assessed after an earthquake based on the judgment of JMA officials. In 1996 the scale was revised as an instrumental seismic intensity measure (I_JMA) that could be used to rapidly assess the expected damage after an earthquake without having to conduct a survey. Since its revision, Japanese researchers have developed several ground motion prediction equations (GMPEs) for I_JMA using Japanese ground motion data. In this paper, we develop a new empirical GMPE for I_JMA based on the strong motion database and functional forms used to develop similar GMPEs for peak response parameters as part of the PEER (Pacific Earthquake Engineering Research Center) Next Generation Attenuation (NGA) project. We consider this relationship to be valid for shallow crustal earthquakes in active tectonic regimes for moment magnitudes ( M ) ranging from 5.0 up to 7.5–8.5 (depending on fault mechanism) and rupture distances ranging from 0 to 200 km. A comparison of this GMPE with relationships developed by Japanese researchers for crustal and shallow subduction earthquakes shows relatively good agreement among all of the relationships at M 7.0 but relatively poor agreement at small magnitudes. Our GMPE predicts the highest intensities at small magnitudes, which together with research on other ground motion parameters, indicates that it provides conservative or upwardly biased estimates of I_JMA for M <5.5. Copyright © 2010 John Wiley & Sons, Ltd.     

Precursors of large Kamchatka earthquakes based on monitoring of seismic noise

This paper presents results from a monitoring study in high frequency seismic noise in Kamchatka during 1992–2006 and reports their use for predicting large regional earthquakes (M ≥ 6.0) in an epicentral distance range Δ within 400 km. The prediction is based on an original method using earth tides as the standard external excitation to study patterns of seismic emission. The method as used at present is described in detail. Guidelines are delineated for future refinement of the prediction method.     

Calculating the seismic energy of earthquakes and explosions

An analysis of the Zoeppritz-Wiechert and Golitsyn equations commonly used to calculate the energy of an earthquake showed that these equations are at variance with the physical assumptions made before deriving them. The value of seismic energy as found from the equations currently used for the purpose is overestimated by 1–2 orders of magnitude. It is suggested that the seismic energy of earthquakes can be found from its value over a half-period of the maximum swing.     

天然地震与人工爆破的波形小波特征研究div

研究了如何从天然地震和人工爆破事件的波形记录中提取出有效、适用的波形特征,以用于对爆破事件的识别．首先对波形记录进行了4层小波包变换;然后对变换得到的最后一层小波包系数提取3种波形特征:能量比特征、香农熵特征及对数能量熵特征;最后利用upsilon;-SVC支持向量分类机对这3种特征的分类能力进行了外推检验．通过选用不同地区、不同台站、不同震级的天然地震与人工爆破的波形记录,力求提取的特征量能尽可能地反映天然地震与人工爆破波形的本质区别,尽量弱化震中距、震级等因素对识别效果的影响．结果表明,上述３种特征中以香农熵特征的识别效果最好,能反映天然地震与人工爆破的本质区别,可作为识别天然地震与人工爆破的一个有效判据．      

Effects of aftershocks on peak ductility demand due to strong ground motion records from shallow crustal earthquakes

Aftershocks induced by a large mainshock can cause additional damage to structures and infrastructure, hampering building reoccupation and restoration activities in a post‐disaster situation. To assess the nonlinear damage potential due to aftershocks, this study investigates the effects of aftershocks by using real as well as artificially generated mainshock–aftershock sequences. The real mainshock–aftershock sequences are constructed from the Pacific Earthquake Engineering Research Center—Next Generation Attenuation database for worldwide shallow crustal earthquakes; however, they are deemed to be incomplete because of missing records. To supplement incomplete real dataset, artificial sequences are generated on the basis of the generalized Omori's law, and a suitable aftershock record selection procedure is then devised to simulate time‐series data for mainshock–aftershock sequences. The results from nonlinear dynamic analysis of inelastic single‐degree‐of‐freedom systems using real and artificial sequences indicate that the incremental effects of aftershocks on peak ductility demand using the real sequences are relatively minor and that peak ductility demand estimates based on the generalized Omori's law are greater, particularly in the upper tail, than those for the real sequences. The results based on the generalized Omori's law also suggest that the aftershock effects based on the real sequences might underestimate the aftershock impact because of the incompleteness of the real dataset. Copyright © 2012 John Wiley & Sons, Ltd.     

Response of offshore towers to strong motion earthquakes

Presented is a stochastic method of analysis of offshore towers subjected to strong motion earthquakes. A zero mean ergodic Gaussian process of finite duration is used to characterize horizontal ground acceleration and full fluid structure interaction effects are included. Numerical results for four representative deep water towers having heights of 475, 675, 875 and 1075 ft are compared with corresponding results obtained by the response spectrum method of analysis. Particular emphasis is placed on the maximum or extreme values of total transverse shear and total overturning moment. Comparisons are made with code values and the role of ductility is briefly discussed.     

Predicting earthquakes by analyzing accelerating precursory seismic activity

During 11 sequences of earthquakes that in retrospect can be classed as foreshocks, the accelerating rate at which seismic moment is released follows, at least in part, a simple equation. This equation (1) is ,where is the cumulative sum until time,t, of the square roots of seismic moments of individual foreshocks computed from reported magnitudes;C andn are constants; andt _fis a limiting time at which the rate of seismic moment accumulation becomes infinite. The possible time of a major foreshock or main shock,t _f,is found by the best fit of equation (1), or its integral, to step-like plots of versus time using successive estimates oft _fin linearized regressions until the maximum coefficient of determination,r ²,is obtained. Analyzed examples include sequences preceding earthquakes at Cremasta, Greece, 2/5/66; Haicheng, China 2/4/75; Oaxaca, Mexico, 11/29/78; Petatlan, Mexico, 3/14/79; and Central Chile, 3/3/85. In 29 estimates of main-shock time, made as the sequences developed, the errors in 20 were less than one-half and in 9 less than one tenth the time remaining between the time of the last data used and the main shock. Some precursory sequences, or parts of them, yield no solution. Two sequences appear to include in their first parts the aftershocks of a previous event; plots using the integral of equation (1) show that the sequences are easily separable into aftershock and foreshock segments. Synthetic seismic sequences of shocks at equal time intervals were constructed to follow equation (1), using four values ofn. In each series the resulting distributions of magnitudes closely follow the linear Gutenberg-Richter relation logN=a–bM, and the productn timesb for each series is the same constant. In various forms and for decades, equation (1) has been used successfully to predict failure times of stressed metals and ceramics, landslides in soil and rock slopes, and volcanic eruptions. Results of more recent experiments and theoretical studies on crack propagation, fault mechanics, and acoustic emission can be closely reproduced by equation (1). Rate-process theory and continuum damage mechanics offer leads toward understanding the physical processes.     

Orientation effect on ground motion measurement for Mexican subduction earthquakes

The existence of the principal directions of the ground motion based on Arias intensity is well-known. These principal directions do not necessarily coincide with the orientations of recording sensors or with the orientations along which the ground motion parameters such as the peak ground acceleration and the pseudo-spectral acceleration (PSA) are maximum. This is evidenced by the fact that the maximum PSA at different natural vibration periods for horizontal excitations do not correspond to the same orientation. A recent analysis carried out for California earthquake records suggests that an orientation-dependent ground motion measurement for horizontal excitations can be developed. The main objective of this study is to investigate and provide seismic ground motion measurements in the horizontal plane, including bidirectional horizontal ground motions, for Mexican interplate and inslab earthquake records. Extensive statistical analyses of PSA are conducted for the assessment, The analysis results suggest that similar to the case of California records, the average behavior of the ratio of the PSA to the maximum resulting PSA can be approximated by a quarter of an ellipse in one quadrant; and that the ratio can be considered to be independent of the value of the maximum resulting PSA, earthquake magnitude, earthquake distance and the focal depth. Sets of response ratios and attenuation relationships that can be used to represent a bidirectional horizontal ground motion measurement for Mexican interplate and inslab earthquakes were also developed.     

Variation of human perceptiveness of earthquakes during seismic sequences

This paper describes the variation of the human perception of earthquakes under the peculiar condition of seismic sequence occurrence. To this aim, we chose and have analyzed four seismic sequences that reflect the most common models of seismicity occurrence in the Italian territory. Our data always refer to the epicentral area, so that the contribution of the epicentral distance to the earthquake sensitiveness is considered constant. To search for recurrent behavior of people in the perception of earthquakes, we crossed seismological data to sensitivity data coming from the archives of the Italian Civil Protection. In each sequence, we individuate a number of “indicators”, such as the rate between felt and not-felt events and the minimum magnitude for which 100% of the events are felt. Such indicators are proportional to the energy released by the sequence and can be particularly affected by the time protraction of the sequence, and/or by the different familiarity with the seismic phenomenon, and/or by the intensity of the emotional impact of a seismic crisis. The observed trends of sensitiveness show a significant variation of the people perceptivity within the sequence development. From the analysis of the dataset, one can distinguish the variation of perceptivity due to emotional factors from those due to physical factors. An important observation, which is confirmed by a statistical analysis, is that the source depth does not affect the perceptiveness of earthquakes. Finally, we find important differences between on-sequence and off-sequence human perceptiveness that are connected with the different levels of alert.     

地震活动关联图与新疆的强震

对新疆1970～2002年的地震资料根据不同范围不同震级进行全时空扫描,认为新疆台网监测能力较好地区在绝大多数6级以上地震前,地震活动关联图较为清晰。提出了预警地震的概念,认为预警地震在短临预报方面有一定的前景。     

软土隧道基于地震响应的输入地震动排序

为探讨“最不利地震动”概念在地下结构抗震设计中的适用性,以软土地铁区间隧道为对象建立相应的地层-结构动力分析模型。以直径变形率为分析指标,基于动力时程方法研究18条不同输入地震动作用下隧道结构动力响应的分布及差异性,得出基于隧道地震响应的输入地震动排序,并通过调幅手段对比分析了地面峰值加速度（PGA）和隧道埋深变化对隧道结构地震动响应排序的影响规律。最后,评价了不同输入地震动参数,包括峰值加速度、峰值速度、峰值位移、绝对累积速度（CAV）和阿里亚斯（Arias）强度（I_A）与隧道地震响应之间的相关性。分析结果表明:① 随着PGA从0.5 m/s2增加到2 m/s2,地震动排序发生明显变化,并且不同输入地震动引起的隧道地震响应差异显著提高,最不利地震动引起的直径变形率与平均值的比值从1.1增加到1.9;② 隧道从浅埋到深埋的过程中,地震动排序结果基本保持不变;③ PGA为2 m/s2时,隧道地震响应与基岩面峰值速度（PBV）的相关性最好,相关系数达到0.94,其次是与基岩面峰值位移（PBD）和I_A,相关系数分别为0.62和0.48,相关性最差的是基岩面峰值加速度（PBA）和CAV,相关系数仅为0.37和0.22。研究结论可为今后软土隧道的输入地震动选择提供科学依据。      

Ground motion duration effects on nonlinear seismic response

The study presented in this paper addresses the question of which nonlinear demand measures are sensitive to ground motion duration by statistical analyses of several case studies. A number of single degree of freedom (SDOF) structures were selected considering: (1) four oscillation periods; (2) three evolutionary and non‐evolutionary hysteretic behaviours; (3) two target ductility levels. Effects of duration are investigated, by nonlinear dynamic analysis, with respect to six different demand indices ranging from displacement ductility ratio to equivalent number of cycles. Input is made of six real accelerogram sets representing three specific duration scenarios (small, moderate and large duration). For all considered demand quantities time‐history results are formally compared by statistical hypothesis test to asses the difference, if any, in the demand concerning different scenarios. Incremental dynamic analysis curves are used to evaluate duration effect as function of ground motion intensity (e.g. spectral acceleration corresponding to the SDOF's oscillation period). Duration impact on structural failure probability is evaluated by fragility curves. The results lead to the conclusion that duration content of ground motion is statistically insignificant to displacement ductility and cyclic ductility demand. The conclusions hold regardless of SDOF's period and hysteretic relationship investigated. Copyright © 2005 John Wiley & Sons, Ltd.     

剪力墙结构对等效地震动输入的影响

采用ANSYS软件建立了考虑土-剪力墙结构相互作用的二维时域动力计算模型,分别考虑了不同土体、上部结构层数对等效输入的影响.根据计算结果,上述因素对等效输入的平动影响不大,但对基础的转动影响很大.特别地,桩基对基础的转动角位移影响很大.因此,分析土-结构相互作用的等效输入,应着重研究基础的转动输入.     

Assessment of the vertical distribution on seismic ground motion

It is very important for the facilities such as nuclear power plants to infer seismic force loading on the earthquake stability assessment of the building foundation and the surrounding slope. The purpose of this paper was to propose a method to evaluate underground seismic coefficients, taking into account dynamic response along the depth in horizontally multi-layered ground. The dynamic property of the seismic coefficient was analyzed on the basis of earthquake records observed at hard and soft rock sites mostly found in Tertiary deposits and sedimentary ground sites of the Pleistocene and Holocene epoch. The evaluation methods of a vertical distribution on underground seismic coefficients were proposed for a few calculation methods on the classified layered ground. Extended evaluation for underground seismic coefficients was confirmed with respect to some multi-layered ground during strong motion.