A new approach for estimating seismic damage of buried water supply pipelines

Conventional damage prediction methods for lifeline structures are primarily based on peak ground motion measurements. However, line structures such as lifelines suffer damage that is mainly induced by the strain of the ground and therefore are likely to be vulnerable to sharp spatial changes in the ground motion. In this study, we propose a measure for evaluating the damage incurred by underground water supply pipelines based on the spatial gradient of the peak ground velocity (PGV), in an attempt to quantify the effects of the geospatial variabilities in the ground motion on pipeline damage. We investigated the spatial distribution of the damage caused to water pipelines during the Niigata‐ken Chuetsu earthquake on October 10, 2004 (Japan Meteorological Agency magnitude (M_JMA) of 6.8) and the Kobe earthquake on January 17, 1995 (M_JMA7.3) and compared the surveyed damage with the PGV distribution as well as with the gradients of the PGV calculated around the damage areas. For the Kobe earthquake, we used the PGV distribution obtained by the strong‐motion simulation performed by Matsushima and Kawase 1 . In case of the Chuetsu earthquake, we estimated the ground motion using a broadband‐frequency‐based strong‐ground‐motion simulation method based on a multiasperity source model. In both cases, we calculated the gradients of the PGV along the geographical coordinates, with the amplitude of the PGV gradient vector being employed as the damage estimator. Our results show that the distribution of damage to underground water supply pipelines exhibits a greater correlation with the gradients of the PGV than with the PGV itself. Thus, the gradient of the PGV is a useful index for preparing initial‐screening hazard maps of underground facilities. Copyright © 2017 John Wiley & Sons, Ltd.     

Seismic response analysis on the stability of running vehicles

The seismometer network of the Japanese expressway system has been enhanced since the 1995 Kobe earthquake. Using earthquake information from the instruments, the expressways are closed if the peak ground acceleration (PGA) is larger than or equal to 80cm/s². The aim of this regulation is to avoid secondary disasters, e.g. cars running into the collapsed sections. However, recent studies on earthquake damage have revealed that expressway structures are not seriously damaged under such‐level of earthquake motion. Hence, we may think of relaxing the regulation of expressway closure. But before doing this, it is necessary to examine the effects of shaking to automobiles since the drivers may encounter difficulties in controlling their vehicles and traffic accidents may occur. In this study, a vehicle was modelled with a six‐degree‐of‐freedom system and its responses were investigated with respect to PGA, peak ground velocity (PGV) and Japan Meteorological Agency (JMA) seismic intensity using five ground motion records. It was observed that the response of the vehicle shows a larger amplitude for the record that has larger response spectrum in the long period range compared to other records. However, similar response amplitudes of the vehicle were observed for all the records with respect to the JMA seismic intensity. The response characteristics of the vehicle model may be very useful for decision‐making regarding the relaxation of the expressway closure under seismic motion. Copyright © 2002 John Wiley & Sons, Ltd.     

Effect of seismic retrofit of bridges on transportation networks

The objective of this research is to determine the effect earthquakes have on the performance of transportation network systems. To do this, bridge fragility curves, expressed as a function of peak ground acceleration (PGA) and peak ground velocity (PGV), were developed. Network damage was evaluated under the 1994 Northridge earthquake and scenario earthquakes. A probabilistic model was developed to determine the effect of repair of bridge damage on the improvement of the network performance as days passed after the event. As an example, the system performance degradation measured in terms of an index, “Drivers Delay,“ is calculated for the Los Angeles area transportation system, and losses due to Drivers Delay with and without retrofit were estimated.     

Effect of earthquake ground motions on fragility curves of highway bridge piers based on numerical simulation

Fragility curves express the probability of structural damage due to earthquakes as a function of ground motion indices, e.g., PGA, PGV. Based on the actual damage data of highway bridges from the 1995 Hyogoken‐Nanbu (Kobe) earthquake, a set of empirical fragility curves was constructed. However, the type of structure, structural performance (static and dynamic) and variation of input ground motion were not considered to construct the empirical fragility curves. In this study, an analytical approach was adopted to construct fragility curves for highway bridge piers of specific bridges. A typical bridge structure was considered and its piers were designed according to the seismic design codes in Japan. Using the strong motion records from Japan and the United States, non‐linear dynamic response analyses were performed, and the damage indices for the bridge piers were obtained. Using the damage indices and ground motion indices, fragility curves for the bridge piers were constructed assuming a lognormal distribution. The analytical fragility curves were compared with the empirical ones. The proposed approach may be used in constructing the fragility curves for highway bridge structures. Copyright © 2001 John Wiley & Sons, Ltd.     

井下、地表地震计记录波形对比——以松原地震台为例

分析了相邻井下、地表地震计记录的背景噪声及近震、远震波形特征,以及相关的振幅谱及近震、远震接收函数波形特征。结果表明,井下地震计记录的背景噪声比地表地震计小1个数量级,其近震、远震波形记录较清晰、简洁;井下、地表地震计记录的远震振幅谱和径向接收函数具有较好的一致性,但背景噪声、近震振幅谱、近震接收函数、远震切向接收函数波形存在一定差异;由地表地震计得到的近震震级比井下地震计的大0.12,地表地震计到井下地震计之间的P波速度约为2 km/s。     

椭圆空洞对盾构隧道地震易损性曲线影响研究

采用增量动力分析方法,探究水平向地震下地层空洞对盾构隧道地震响应特征的影响规律。针对管片损伤及周边地层应力,选择弯矩比作为性能评价指标,峰值加速度(PGA)及峰值速度(PGV)作为衡量地震强度指标(IM),阐明椭圆形空洞对管片抗震性能的影响,得到隧道结构的地震易损性曲线。研究表明：椭圆空洞加大了浅埋盾构隧道的地震破坏概率;PGA与PGV均可作为IM并获得相应的隧道易损性曲线;使用弯矩比作为破坏指标,PGV作为地震动指标,其对应的易损性曲线对地层变异性更敏感。研究结论可为潜在空洞发育区防震方案的制定提供参考。     

Numerical simulation of strong ground motion for the <Emphasis Type="Italic">M</Emphasis><Subscript>s</Subscript>8.0 Wenchuan earthquake of 12 May 2008

The Wenchuan earthquake of 12 May 2008 is the most destructive earthquake in China in the past 30 years in terms of property damage and human losses. In order to understand the earthquake process and the geo-morphological factors affecting the seismic hazard, we simulated the strong ground motion caused by the earthquake, incorporating three-dimensional (3D) earth structure, finite-fault rupture, and realistic surface topography. The simulated ground motions reveal that the fault rupture and basin structure control the overall pattern of the peak ground shaking. Large peak ground velocity (PGV) is distributed in two narrow areas: one with the largest PGV values is above the hanging wall of the fault and attributed to the locations of fault asperities and rupture directivity; the other is along the northwestern margin of the Sichuan Basin and caused by both the directivity of fault rupture and the amplification in the thick sediment basin. Rough topography above the rupture fault causes wave scattering, resulting in significantly larger peak ground motion on the apex of topographic relief than in the valley. Topography and scattering also reduce the wave energy in the forward direction of fault rupture but increase the PGV in other parts of the basin. These results suggest the need for a localized hazard assessment in places of rough topography that takes the topographic effects into account. Finally, had the earthquake started at the northeast end of the fault zone and ruptured to the southwest, Chengdu would have suffered a much stronger shaking than it experienced on 12 May, 2008. Supported by the U.S. National Science Foundation (Grant Nos. EAR 0738779 and OCE 0727919), the National Basic Research Program of China (Grant No. 2004CB418404), and partially by the National Nature Science Foundation of China (Grant No. 40521002)     

Models for reproducing the damage scenario of the Lorca earthquake

A damage scenario modelling is developed and compared with the damage distribution observed after the 2011 Lorca earthquake. The strong ground motion models considered include five modern ground motion prediction equations (GMPEs) amply used worldwide. Capacity and fragility curves from the Risk-UE project are utilized to model building vulnerability and expected damage. Damage estimates resulting from different combinations of GMPE and capacity/fragility curves are compared with the actual damage scenario, establishing the combination that best explains the observed damage distribution. In addition, some recommendations are proposed, including correction factors in fragility curves in order to reproduce in a better way the observed damage in masonry and reinforce concrete buildings. The lessons learned would contribute to improve the simulation of expected damages due to future earthquakes in Lorca or other regions in Spain with similar characteristics regarding attenuation and vulnerability.     

Observational fragility functions for residential stone masonry buildings in Nepal

This paper outlines the seismic vulnerability of rural stone masonry buildings affected by the 2015 Gorkha earthquake sequence. Summary of field observation is presented first and empirical fragility curves are developed from the detailed damage assessment data from 603 villages in central, eastern and western Nepal. Fragility curves are developed on the basis of 665,515 building damage cases collected during the post-earthquake detailed damage assessment campaign conducted by Government of Nepal. Two sets of fragility functions are derived using peak ground acceleration and spectral acceleration at 0.3 s as the intensity measures. The sum of the results highlights that stone masonry buildings in Nepal are highly vulnerable even in the case of low to moderate ground shaking. The results further indicate that in the case of strong to major earthquakes, most of the stone masonry buildings in Nepal would sustain severe damage or collapse.     

Analytical seismic fragility functions for highway and railway embankments and cuts

A fundamental tool in seismic risk assessment of transportation systems is the fragility curve, which describes the probability that a structure will reach or exceed a certain damage state for a given ground motion intensity. Fragility curves are usually represented by two‐parameter (median and log‐standard deviation) cumulative lognormal distributions. In this paper, a numerical approach, in the spirit of the IDA, is applied for the development of fragility curves for highways and railways on embankments and in cuts due to seismic shaking. The response of the geo‐construction to increasing levels of seismic intensity is evaluated using a 2D nonlinear finite element model, with an elasto‐plastic criterion to simulate the soil behavior. A calibration procedure is followed in order to account for the dependency of both the stiffness and the damping to the soil strain level. The effect of soil conditions and ground motion characteristics on the response of the embankment and cut is taken into account considering different typical soil profiles and seismic input motions. This study will provide input for the assessment of the vulnerability of the road/railway network regarding the performance of the embankments and cuts; therefore, the level of damage is described in terms of the permanent ground displacement in these structures. The fragility curves are estimated based on the evolution of damage with increasing earthquake intensity, which is described by PGA. The proposed approach allows the evaluation of new fragility curves considering the distinctive features of the element's geometry, the input motion, and the soil properties as well as the associated uncertainties. A relationship between the computed permanent ground displacement on the surface of the embankment and the PGA in the free field is also suggested based on the results of the numerical analyses. Finally, the proposed fragility curves are compared with existing empirical data and the limitations of their applicability are outlined. Copyright © 2015 John Wiley & Sons, Ltd.     

地震预警PGV-Pd关系参数的距离分段特征

地震预警系统需要在破坏性的地震波到来前快速估算地震参数和地震动参数,以对可能出现的地震灾害进行预测,对重要工程、人员密集区域发布警报信息.以Pd估测PGV的方法是地震预警研究涉及的一种重要问题,该方法利用初至P波触发后前几秒的峰值位移(Pd)对最终地震动峰值速度(PGV)进行估算,以满足预警的需要.本文对2016年在日本发生的熊本地震及其前震、余震的震中距100km以内、矩震级大于4级、井下基岩PGA5cm·s~(-2)和地表PGA20cm·s~(-2)的Kik-net强地震动记录进行处理分析,用于研究以Pd估测PGV的方法.将获得的强地震动数据按震中距从0~100km平均划分为5个区段,在记录时间3~10s范围内将Pd的计算取8个时间窗,分别对每个震中距区段、每一个Pd的时间窗下的PGV-Pd数据进行线性拟合,最终提出了一套应对不同震中距对位移幅值连续追踪测定PGV的算法.对每一个震中距区段的研究表明,震中距会对PGV-Pd关系产生影响.对5次地震进行验证分析,认为基于基岩记录估测基岩PGV的准确度高于基于地表记录估测地表PGV的准确度;对震中距进行分段的PGV估测方法准确度高于不考虑震中距因子的估测方法.最后拓展了将井下基岩Pd估测井下基岩PGV这种原地地震预警方法,使其能够为异地P波预警方法服务.     

Numerical simulation of strong ground motion for the M_s8.0 Wenchuan earthquake of 12 May 2008

The Wenchuan earthquake of 12 May 2008 is the most destructive earthquake in China in the past 30 years in terms of property damage and human losses. In order to understand the earthquake process and the geo-morphological factors affecting the seismic hazard, we simulated the strong ground mo-tion caused by the earthquake, incorporating three-dimensional (3D) earth structure, finite-fault rupture, and realistic surface topography. The simulated ground motions reveal that the fault rupture and basin structure control the overall pattern of the peak ground shaking. Large peak ground velocity (PGV) is distributed in two narrow areas: one with the largest PGV values is above the hanging wall of the fault and attributed to the locations of fault asperities and rupture directivity; the other is along the north-western margin of the Sichuan Basin and caused by both the directivity of fault rupture and the ampli-fication in the thick sediment basin. Rough topography above the rupture fault causes wave scattering, resulting in significantly larger peak ground motion on the apex of topographic relief than in the valley. Topography and scattering also reduce the wave energy in the forward direction of fault rupture but increase the PGV in other parts of the basin. These results suggest the need for a localized hazard as-sessment in places of rough topography that takes the topographic effects into account. Finally, had the earthquake started at the northeast end of the fault zone and ruptured to the southwest, Chengdu would have suffered a much stronger shaking than it experienced on 12 May, 2008.     

Ground-motion prediction equations for the intermediate depth Vrancea (Romania) earthquakes

We present the regional ground-motion prediction equations for peak ground acceleration (PGA), peak ground velocity (PGV), pseudo-spectral acceleration (PSA), and seismic intensity (MSK scale) for the Vrancea intermediate depth earthquakes (SE-Carpathians) and territory of Romania. The prediction equations were constructed using the stochastic technique on the basis of the regional Fourier amplitude spectrum (FAS) source scaling and attenuation models and the generalised site amplification functions. Values of considered ground motion parameters are given as the functions of earthquake magnitude, depth and epicentral distance. The developed ground-motion models were tested and calibrated using the available data from the large Vrancea earthquakes. We suggest to use the presented equations for the rapid estimation of seismic effect after strong earthquakes (Shakemap generation) and seismic hazard assessment, both deterministic and probabilistic approaches.     

Long-period ground motions from the 2003 Tokachi-oki earthquake

We study characteristics of long-period ground motions from the 2003 Tokachi-oki earthquake (Mj 8.0), a large interplate earthquake, based on spatial distribution maps and attenuation relationships for four kinds of peak ground velocity (PGV) value. The first kind (PGV(WB)) is obtained from a maximal value of vector sum of the three-component, wide-band velocity seismograms, and the other three kinds (PGV(BP10), PGV(BP20), and PGV(BP30)) are obtained from a maximal value of vector sum of the three-component, narrow band-pass filtered velocity seismograms (the central periods are 10, 20, and 30 s). The spatial distribution maps for all kinds of PGV value show azimuth dependence; the PGV values in Hokkaido, northern side of the epicenter are larger than those in Tohoku, southwestern side of the epicenter, when compared at a comparable distance. We find that the features result from the radiation pattern of long-period surface waves, that is, the source effect. The attenuation relationships show the following trends: The PGV(WB) values are larger than the sum of the PGV(BP10), PGV(BP20), and PGV(BP30) at distances (D) less than 200 km, while the PGV(WB) values are comparable to the sum of the PGV(BP20) and PGV(BP30) at D > 200 km. This indicates that the PGV(WB) values at D < 200 km are affected by ground motions with periods less than 10 s, while long-period surface waves mainly contribute to the PGV(WB) values at D > 200 km. The basin site effects generate a patchy pattern in the spatial distribution maps and a large scattering in the attenuation relationships for the PGV(WB) and PGV(BP10) values. Finally, we conclude that the PGV(WB) values from the 2003 Tokachi-oki earthquake are controlled by the radiation pattern of long-period S and surface waves and various basin site effects.     

不同地震动作用下隔震结构减震影响分析

基于地震动加速度峰值和速度峰值比值（PGA／PGV）选用了13条不同频谱的地震波,对不同参数条件下多层基础隔震结构地震响应进行非线性时程分析,得到了不同地震动作用下隔震结构模型的地震反应。计算结果表明：随PGA／PGV增大,隔震结构减震率整体上呈增大趋势,但是PGA／PGV对减震率的影响程度还受结构自振周期,上部结构侧移刚度的影响,应综合考虑。     

Characteristics of the strong ground motions from the 6 April 2009 L’Aquila earthquake,Italy

An M_w 6.25 earthquake occurred on April 6, 2009 at 03:33 a.m. local time, in the Abruzzo region (Central Italy), close to the city of L’Aquila. The earthquake ruptured a North-West (NW)–South-East (SE) oriented normal fault dipping toward the South-West (SW), with the city of L’Aquila lying a few kilometers away on the hanging wall.The main shock has been recorded by fifty-eight accelerometric stations: the highest number of digital recordings ever obtained in Italy for a single earthquake, one of the best-recorded earthquakes with a normal fault mechanism. Very high values of peak ground acceleration (0.3–0.65 g) were observed close to the center of L’Aquila (6 stations at zero JB distance from the fault). The earthquake caused severe loss of lives (299 victims and 1500 injured) and damage (about 18000 unusable buildings) in the epicentral area.In this study we analyze the ground motion characteristics of both the main shock in terms of peak ground acceleration (PGA), peak ground velocity (PGV), and pseudo-acceleration response spectra (5% of damping ratio). In particular, we compare the pseudo-acceleration response spectra for horizontal directions with the EC8 design spectrum and the new Italian building code (NTC08). In order to understand the characteristics of the ground motions induced by L’Aquila earthquake, we also study the source-related effects and site response of the strong motion stations that recorded the seismic sequence. A novel method is used for the analysis of inter-station and site-specific H/V spectral ratios for the main event and for 12 aftershocks.     

Real-time prediction of earthquake potential damage: A case study for the January 8, 2022 MS 6.9 Menyuan earthquake in Qinghai,China

It is critical to determine whether a site has potential damage in real-time after an earthquake occurs, which is a challenge in earthquake disaster reduction. Here, we propose a real-time Earthquake Potential Damage predictor (EPDor) based on predicting peak ground velocities (PGVs) of sites. The EPDor is composed of three parts: (1) predicting the magnitude of an earthquake and PGVs of triggered stations based on the machine learning prediction models; (2) predicting the PGVs at distant sites based on the empirical ground motion prediction equation; (3) generating the PGV map through predicting the PGV of each grid point based on an interpolation process of weighted average based on the predicted values in (1) and (2). We apply the EPDor to the 2022 M_S 6.9 Menyuan earthquake in Qinghai Province, China to predict its potential damage. Within the initial few seconds after the first station is triggered, the EPDor can determine directly whether there is potential damage for some sites to a certain degree. Hence, we infer that the EPDor has potential application for future earthquakes. Meanwhile, it also has potential in Chinese earthquake early warning system.     

Simulation of accelerograms from simplified deterministic approach for the 23rd October 2004 Niigata-ken Chuetsu,Japan earthquake

A simple hybrid approach for the simulation of strong ground motion is presented in this paper. This approach is based on the deterministic modelling of rupture plane initially started by Midorikawa, Tectonophysics 218:287–295, (1993) and further modified by Joshi, Pure Appl Geophys (PAGEOPH) 8:161, (2004). In this technique, the finite rupture plane of the target event is divided into several subfaults, which satisfy scaling relationship. In this paper, simulation of strong ground motion due to a rupture buried in a earth medium consisting of several layers of different velocities and thicknesses is made by considering (1) transmission of energy at each layer; (2) frequency filtering properties of medium and earthquake source; (3) correction factor for slip of large and small magnitude earthquakes and (4) site amplification ratio at various stations. To test the efficacy of the developed technique, strong motion records were simulated at different stations that have recorded the 2004 Niigata-ken Chuetsu, Japan earthquake (M _s 7.0). Comparison is made between the simulated and observed velocity and acceleration records and their response spectra. Distribution of peak ground acceleration, velocity and displacement surrounding the rupture plane is prepared from simulated and observed records and are compared with each other. The comparison of synthetic with the observed records over wide range of frequencies shows that the present technique is effective to predict various strong motion parameters from simple deterministic model which is based on simple regression relations and modelling parameters.     

Analysis of Applicability of Seismic Intensity Estimation Using Fourier Spectrum of Ground Motions for Iranian Earthquakes

The intensity scales in different forms provide valuable information on regional earthquake effects. In this paper, a theoretical model which has been developed recently for seismic intensity estimation is re-examined by employing strong motion records from ten Iranian earthquakes. The analysis results confirmed the capability of the implemented method to estimate the seismic intensity in terms of the MMI scale based on a Fourier spectrum in the study area. The predicted intensity values were compared with another technique utilizing peak ground velocity (PGV) as a predictor. To reveal the high potential of the adopted approach, the theoretical isoseismal map was developed for the 1978 Tabas, Iran earthquake (Mw = 7.4) based on a stochastic finite-fault modelling of ground motions. Results showed good compatibility of predicted intensity values while the historical earthquake records are not enough for a given site.     

Seismic fragility of base-isolated water storage tanks under non-stationary earthquakes

Seismic fragility curves for fixed-base and base-isolated liquid storage tanks are developed under non-stationary earthquakes, and their seismic performance are compared. The correlation between different earthquake intensity measure (IM) parameters and peak response quantities of the base-isolated liquid storage tanks are investigated. The failure criteria are chosen based on (1) the elastic buckling strength of the tank wall, which is defined in terms of critical base shear and critical overturning moment, and (2) in terms of the critical isolation displacement. The uncertainty involved is considered in the earthquake characteristics. Non-stationary earthquake ground motions are generated using Monte Carlo (MC) simulation. Influence of the isolator characteristic parameters and modeling approaches on the seismic fragility of the base-isolated liquid storage tanks is also investigated. Peak ground acceleration is found to be the well correlated IM parameter with the peak response quantities of the base-isolated liquid storage tanks. Substantial decrease in the seismic fragility of the base-isolated liquid storage tanks is observed as compared to the fixed-base tanks. Significant influence of the isolator characteristic parameters on the seismic fragility of the base-isolated liquid storage tanks are reported in the present study.