GPU-acceleration 3D rotated-staggered-grid solutions to microseismic anisotropic wave equation with moment tensor implementation

To improve the accuracy of microseismic inversion,seismic anisotropy and moment tensor source should be carefully considered in the forward modelling stage.In this study,3D microseismic anisotropy wave forward modelling with a moment tensor source was proposed.The modelling was carried out based on a rotated-staggered-grid(RSG)scheme.In contrast to staggered-grids,the RSG scheme defines the velocity components and densities at the same grid,as do the stress components and elastic parameters.Therefore,the elastic moduli do not need to be interpolated.In addition,the detailed formulation and implementation of moment-tensor source loaded on the RSG was presented by equating the source to the stress increments.Meanwhile,the RSG-based 3D wave equation forward modelling was performed in parallel using compute unified device architecture(CUDA)programming on a graphics processing unit(GPU)to improve its efficiency.Numerical simulations including homogeneous and anisotropic models were carried out using the method proposed in this paper,and compared with other methods to prove the reliability of this method.Furthermore,the high efficiency of the proposed approach was evaluated.The results show that the computational efficiency of proposed method can be improved by about two orders of magnitude compared with traditional central processing unit(CPU)computing methods.It could not only help the analysis of microseismic full wavefield records,but also provide support for passive source inversion,including location and focal mechanism inversion,and velocities inversion.

Seismic safety assessment with non-Gaussian random processes for train-bridge coupled systems

Extensive high-speed railway(HSR)network resembled the intricate vascular system of the human body,crisscrossing mainlands.Seismic events,known for their unpredictability,pose a significant threat to both trains and bridges,given the HSR’s extended operational duration.Therefore,ensuring the running safety of train-bridge coupled(TBC)system,primarily composed of simply supported beam bridges,is paramount.Traditional methods like the Monte Carlo method fall short in analyzing this intricate system efficiently.Instead,efficient algorithm like the new point estimate method combined with moment expansion approximation(NPEM-MEA)is applied to study random responses of numerical simulation TBC systems.Validation of the NPEM-MEA’s feasibility is conducted using the Monte Carlo method.Comparative analysis confirms the accuracy and efficiency of the method,with a recommended truncation order of four to six for the NPEM-MEA.Additionally,the influences of seismic magnitude and epicentral distance are discussed based on the random dynamic responses in the TBC system.This methodology not only facilitates seismic safety assessments for TBC systems but also contributes to standard-setting for these systems under earthquake conditions.

RESEARCH ON THE DEFORMATIONS IN QINGHAI_TIBET PLATEAU AND ITS MARGINS BY INVERTING SEISMIC MOMENT TENSORS AND GPS VELOCITIES

We have determined approximate average rates of deformation in the Qinghai_Tibet plateau and its margins from the GPS data for last 10 years and the moment tensors from earthquakes between 1900 and 1999.We also determined the strain rate (seismic strain rate) associated with the seismic deformation using 254 M w ≥5.0 earthquakes,and estimated the shortening and extension rates for every block in the area as well.We also estimated the strain rate (geodetic strain rate)by 80 GPS sites’ velocity vectors and analyzed characteristic of kinematics by two kinds of strain rates and discussed earthquake potential in the area.As a result,the deformation rates from seismic moment tensors and from GPS velocities are basically agreed with each other.It is feasible to analyze seismic risk by comparing geodetic strain rate with seismic strain rate based on the opinion that strain energy will be released through earthquake.It is concluded that there is no strong earthquake potential (>M7) in the Qinghai_Tibet plateau and its margins,but there is earthquake potential (>M5) in middle Tibet in a few years.

Effectiveness of two conventional methods for seismic retrofit of steel and RC moment resisting frames based on damage control criteria

This study investigates the efficiency of two types of rehabilitation methods based on economic justification that can lead to logical decision making between the retrofitting schemes. Among various rehabilitation methods, concentric chevron bracing（CCB） and cylindrical friction damper（CFD） were selected. The performance assessment procedure of the frames is divided into two distinct phases. First, the limit state probabilities of the structures before and after rehabilitation are investigated. In the second phase, the seismic risk of structures in terms of life safety and financial losses（decision variables） using the recently published FEMA P58 methodology is evaluated. The results show that the proposed retrofitting methods improve the serviceability and life safety performance levels of steel and RC structures at different rates when subjected to earthquake loads. Moreover, these procedures reveal that financial losses are greatly decreased, and were more tangible by the application of CFD rather than using CCB. Although using both retrofitting methods reduced damage state probabilities, incorporation of a site-specific seismic hazard curve to evaluate mean annual occurrence frequency at the collapse prevention limit state caused unexpected results to be obtained. Contrary to CFD, the collapse probability of the structures retrofitted with CCB increased when compared with the primary structures.

Applications of seismic moment tensor inversion in fast response to earthquakes

Using the technique of seismic moment tensor inversion, the source mechanisms of 10 earthquakes with Ms5.2that occurred in China from November 1996 to January 1998 were determined rapidly. The determined resultswere sent as 'Bulletins of Source Mechanism Parameters of Earthquakes' to the Seismic Regime Guards' Office,China Seismological Bureau, and the relevant provincial seismological bureaus. These bulletins have played rolein the fast response to large earthquakes.

Optimal seismic retrofit model for steel moment resisting frames with brittle connections

Based on performance-based seismic engineering, this paper proposes an optimal seismic retrofit model for steel moment resisting frames（SMRFs） to generate a retrofit scheme at minimal cost. To satisfy the acceptance criteria for the Basic Safety Objective（BSO） specified in FEMA 356, the minimum number of upgraded connections and their locations in an SMRF with brittle connections are determined by evolutionary computation. The performance of the proposed optimal retrofitting model is evaluated on the basis of the energy dissipation capacities, peak roof drift ratios, and maximum interstory drift ratios of structures before and after retrofitting. In addition, a retrofit efficiency index, which is defined as the ratio of the increment in seismic performance to the required retrofitting cost, is proposed to examine the efficiencies of the retrofit schemes derived from the model. The optimal seismic retrofit model is applied to the SAC benchmark examples for threestory and nine-story SMRFs with brittle connections. Using the retrofit efficiency index proposed in this study, the optimal retrofit schemes obtained from the model are found to be efficient for both examples in terms of energy dissipation capacity, roof drift ratio, and maximum inter-story drift ratio.

A new methodology for energy-based seismic design of steel moment frames

A procedure is proposed whereby input and hysteretic energy spectra developed for single-degree-of-freedom （SDOF） systems are applied to multi-degree-of-freedom （MDOF） steel moment resisting frames. The proposed procedure is verified using four frames, viz., frame with three-, five-, seven- and nine-stories, each of which is subjected to the fault- normal and fault-parallel components of three actual earthquakes. A very good estimate for the three- and five-story frames, and a reasonably acceptable estimate for the seven-, and nine-story frames, have been obtained. A method for distributing the hysteretic energy over the frame height is also proposed. This distribution scheme allows for the determination of the energy demand component of a proposed energy-based seismic design （EBSD） procedure for each story. To address the capacity component of EBSD, a story-wise optimization design procedure is developed by utilizing the energy dissipating capacity from plastic hinge formation/rotation for these moment frames. The proposed EBSD procedure is demonstrated in the design of a three-story one-bay steel moment frame.

Evaluation of the seismic moment of the April 20,2013 Lushan earthquake

The April 20, 2013 Lushan earthquake which occurred in Sichuan, China had only moderate thrust. However, the computed seismic moments （M0） for the Lushan earthquake calculated by several institutions differ significantly from 0.4 × 1019 to 1.69 ×019 Nm, up to four times difference. We evaluate ten computed Mos by using normal mode observations from superconducting gravimeters in China's Mainland. We compute synthetic normal modes on the basis of moment tensor solutions and fit them to the observed normal modes. Comparison of our results indicates that Mo is the main cause for some large differences between observations and synthetics. We sug- gest that a moment magnitude of Mw6.6, corresponding to a Mo of 0.97-1.08 × 1019 Nm, characterizes the size and strength of the seismic source of the Lushan earthquake.

Estimating seismic moments and Lg Q using Lg spectra

A multi-event and multi-station inverse method is presented in the paper to simultaneously estimate the seismic moments (M0) and source corner frequencies (fc) of several Jiashi (Xinjiang, China) earthquakes, as well as the apparent Lg Q models for the paths from Jiashi to eight seismic stations (WMQ, AAK, TLG, MAKZ, KUR, VOS, ZRN and CHK) in Central Asia. The resultant seismic moments correlate well with the M0 values obtained by Harvard University using the centroid moment tensor (CMT) inversion and the surface-wave magnitudes as well. After the correction by a typical value of average radiation coefficient for regional SV waves, the M0 values from Lg spectral inversion are still close to the corresponding values obtained from CMT inversion. The obtained ap- parent Q0Lg values (Lg Q at 1 Hz) are consistent with the tectonic features of corresponding propagation paths. The Q0Lg values are 351±87, 349±86 and 300±27 for the paths from Jiashi to AAK, TLG and MAKZ, respectively. They are smaller than Q0Lg values for the paths to KUR, VOS, ZRN and CHK, which are 553±72, 569±58, 550±57 and 603±65, respectively. These results agree with the condition that the paths to AAK, TLG and MAKZ mainly propagate through the mountainous Tianshan area where relatively strong seismic activities and large variations of topography are exhibited, while the paths to KUR, VOS, ZRN and CHK mainly propagate through the stable area of Kazak platform. The Q0Lg value for the path to WMQ is 462±56. This is also in agreement with the condition that the path to WMQ is basically along the border area between Tianshan Mountain and Tarim Basin, and along this path the variations of topography and crustal thickness are moderate in comparison with that along the path to MAKZ.

Seismic moment tensor representations and radiation patterns in unbounded media with elipsoidal cavities driven by low frequency pressure

Seismic waves generated in elastic media by ellipsoidal cavities driven by low frequency pressure can be solved by Eshelby′s method. In this paper we prove that this method, originally presented for elastostatic ellipsoidal inclusion problems, can also be applied to dynamic problems on the assumption that wavelengths are much longer than the dimensions of cavities. Comparison between the approximate solution for spherical cavities and its relevant precise solution shows that this method may be used for radii of cavities smaller than 100 m and for frequencies below 10 Hz. By this approach we show that, for low frequency seismic waves in the far field an explosive load on an ellipsoidal cavity can be equivalent to three dipoles oriented along the principal axes of that cavity. Seismic radiation patterns are given analytically and results thus obtained show that an explosion detonated in an ellipsoidal cavity may radiate significant S waves.

Seismic design factors for RC special moment resisting frames in Dubai,UAE

This study investigates the seismic design factors for three reinforced concrete （RC） framed buildings with 4, 16 and 32-stories in Dubai, UAE utilizing nonlinear analysis. The buildings are designed according to the response spectrum procedure defined in the 2009 International Building Code （IBC＇09）. Two ensembles of ground motion records with 10% and 2% probability of exceedance in 50 years （10/50 and 2/50, respectively） are used. The nonlinear dynamic resPonses to the earthquake records are computed using IDARC-2D. Key seismic design parameters are evaluated; namely, response modification factor （R）, deflection amplification factor （Cd）, system overstrength factor （Ωo）, and response modification factor for ductility （Rd） in addition to inelastic interstory drift. The evaluated seismic design factors are found to significantly depend on the considered ground motion （10/50 versus 2/50）. Consequently, resolution to the controversy of Dubai seismicity is urged. The seismic design factors for the 2/50 records show an increase over their counterparts for the 10/50 records in the range of 200%-400%, except for the D~ factor, which shows a mere 30% increase. Based on the observed trends, perioddependent R and Cd factors are recommended if consistent collapse probability （or collapse prevention performance） in moment frames with varying heights is to be expected.

Full moment tensor inversion constrained by doublecouple focal mechanism for induced seismicity

In this study,we propose a new method to determine full moment tensor solution for induced seismicity.This method generalizes the full waveform matching algorithm we have developed to determine the double-couple(DC)focal mechanism based on the neighbourhood algorithm.One major difference between the new method and the former one is that we adopt a new misfit function to constrain the candidate moment tensor solutions with respect to a reference DC solution in addition to other misfit terms characterizing the waveform matching.Through synthetic tests using a real passive seismic survey geometry,the results show the new constraint can help better recover the DC components of inverted moment tensors.We further investigate how errors in the velocity model and source location affect the moment tensor solution.The synthetic test results indicate that the constrained inversion is robust in recovering both the DC and non-DC components.We also test the proposed method on several real induced events in an oil/gas field in Oman using the same observation system as synthetic tests.While it is found that the full moment tensor solutions without using the DC constraints have much larger non-DC components than solutions with the DC constraints,both solutions are able to fit the observed waveforms at similar levels.The synthetic and real test results suggest the proposed DC constrained inversion method can reliably retrieve full moment tensor solutions for the induced seismicity.

Control Parameters of Magnitude—Seismic Moment Correlation for the Crustal Earthquakes

In connection with conversion from energy class KR (KR = log10E R, where ER — seismic energy, J) to the universal magnitude estimation of the Tien Shan crustal earthquakes the development of the self-coordinated correlation of the magnitudes (mb , ML, Ms ) and KR with the seismic moment M0 as the base scale became necessary. To this purpose, the first attempt to develop functional correlations in the magnitude—seismic moment system subject to the previous studies has been done. It is assumed that in the expression M (mb , ML , Ms) = Ki + zi log10M0 , the coefficients ki? and zi? are controlled by the parameters of ratio?(where;f0 —corner frequency, Brune, 1970, 1971;M0, N×m). According to the new theoretical predictions common functional correlation of the advanced magnitudes Mm (mbm = mb , MLm = ML , MSm = MS ) from log10M0,? log10t0? and the elastic properties (Ci) can be presented as , where , and , for the averaged elastic properties of the Earth’s crust for thembmthe coefficients Ci= –11.30 and di = 1.0, for MLm: Ci = –14.12, di = 7/6;for MSm : Ci = –16.95 and di = 4/3. For theTien Shan earthquakes (1960-2012 years) it was obtained that , and on the basis of the above expressions we received that MSm = 1.59mbm – 3.06. According to the instrumental data the correlation Ms = 1.57mb – 3.05 was determined. Some other examples of comparison of the calculated and observed magnitude - seismic moment ratios for earthquakes of California, the Kuril Islands, Japan, Sumatra and South America are presented.

Seismic Evaluation of Steel Moment Resisting Frames (MRFs)—Supported by Loose Granular Soil

Soil underneath a structure might affect the behavior and the overall response of the structure in seismic events. The role of loose soil conditions and the inclusion of soil-structure interaction (SSI) in the analysis are important issues that need to be addressed. Since steel structures are light, two configurations designed as spatial and perimeter are considered to study the effect of soil on the steel structural frames for the same building. The paper provides a parametric analysis on the influence of SSI on the overall performance of MRFs (Moment Resisting Frames) according to the provisions of Saudi Building Code (SBC) [1]. A case study has been developed in which spatial and perimeter moment resisting frames of 12, 6 and 3 stories residential buildings are designed using Saudi Building Code (SBC) prescriptions. A modal response spectrum analysis has been carried out to see the influence of SSI on the fundamental period of vibration, top story displacement and inter-story drift limitations. Moreover, a static non-linear analysis has been performed to investigate the performance of frames, thus allowing to identify the influence of SSI on the structural design of steel MRFs.

A Catalogue of Earthquakes With Moment Magnitude of the Center-Southern Asia Test Area of the Global Seismic Hazard Assessment Program

We compile a regional catalogue of earthquakes with moment magnitude of the Center-Southern Asia test area(20°-35°N,85°-105°E)for the Global Seismic Hazard Assessment Program(GSHAP).There are significant inhomogeneous and uncompleted data,and n uniformity of earthquake magnitudes in this test area because this region is situated on the boundaries of many countries,such as China,India,Nepal,Vietnam,etc.We establish a relationship between Gutenberg surface-wave magnitude and IASPEI surface-wave magnitude,which can be used for conversion of different magnitude scales into moment magnitude for this catalogue.A catalogue of events with Mw≥6.0 of this test area is given at the end of this paper.

Moment tensor and stress inversion solutions of acoustic emissions during compression and tensile fracturing in crystalline rocks

We investigate the accuracy and robustness of moment tensor(MT)and stress inversion solutions derived from acoustic emissions(AEs)during the laboratory fracturing of prismatic Barre granite specimens.Pre-cut flaws in the specimens introduce a complex stress field,resulting in a spatial and temporal variation of focal mechanisms.Specifically,we consider two experimental setups:(1)where the rock is loaded in compression to generate primarily shear-type fractures and(2)where the material is loaded in indirect tension to generate predominantly tensile-type fractures.In each test,we first decompose AE moment tensors into double-couple(DC)and non-DC terms and then derive unambiguous normal and slip vectors using k-means clustering and an unstructured damped stress inversion algorithm.We explore temporal and spatial distributions of DC and non-DC events at different loading levels.The majority of the DC and the tensile non-DC events cluster around the pre-cut flaws,where macro-cracks later develop.Results of stress inversion are verified against the stress field from finite element(FE)modeling.A good agreement is found between the experimentally derived and numerically simulated stress orientations.To the best of the authors’knowledge,this work presents the first case where stress inversion methodologies are validated by numerical simulations at laboratory scale and under highly heterogeneous stress distributions.

2023年12月18日甘肃积石山M_(S)6.2地震多点源地震矩张量反演研究

基于区域宽频带数字地震台网观测到的地震记录,采用多点源地震矩张量反演方法对2023年甘肃积石山M_(S)6.2地震进行了地震矩张量反演研究。结果显示,该地震震源机制解的NE倾的节面(走向312°,倾角48°,滑动角76°)为发震断裂,主体地震矩从震源向地表沿滑动矢量的方向(方位角236°,倾伏角40.5°)释放,主破裂持续时间约为5.8 s,与此同时,地震破裂区断层面上发生了次级的双侧破裂。

中国青藏高原M≥6.0地震的辐射能量测定

地震辐射能量主要由高频体波携带,与震源动态特征直接相关,能够有效地弥补地震矩和矩震级对震源动态过程和高频成分描述的不足,为地震应急和灾害评估提供更为全面的参考。本文使用自主研发的地震辐射能量测定软件测定了中国青藏高原自1990年以来M≥6.0浅源地震的地震辐射能量,并初步分析了中国青藏高原地区浅源地震的能量释放规律。研究结果表明:(1)利用该软件测定的中国青藏高原地区的34例浅源地震所得的测定结果稳定可靠。在一般情况下,地震的矩震级M_(W)与能量震级M_(e)并不相等,为了更全面地评估地震动效应,不仅需要考虑地震断层错动的静态特征,还需要考虑震源的动态特性。因此,在表示震源特性时,矩震级M_(W)和能量震级M_(e)各有特点,优势互补。(2)改则地震和尼玛地震的能量震级分别为6.47和6.00,地震辐射能量分别为1.26×10^(14) J和0.25×10^(14) J。两次震源位置相近、矩震级相同的地震,所释放的能量相差5倍之多。为了探究改则地震与尼玛地震能量释放差异的原因,我们通过S变换对两次地震的波形进行了时频分析。结果表明,震源动态过程的不同导致了地震辐射能量释放过程的差异。因此能量震级和地震辐射能量等动态震源参数比矩震级更适合描述震源的动态过程和地震的潜在破坏性。(3)结合地震矩资料得到中国青藏高原地区浅源地震的平均能矩比为1.9×10^(-5),是全球浅源地震平均能矩比的1.6倍。能矩比的大小与震源机制有关,走滑型地震的平均能矩比高于倾滑型地震。(4)中国青藏高原浅源地震的能矩比分布存在区域特征,东部和西部平均能矩比分别为2.25×10^(-5)和1.62×10^(-5),东部明显高于西部。能矩比的区域性差异与地质构造背景有关,地震辐射能量能够反映该地区的地质构造状态。同一区域内地震的能矩比也存在差别,中国青藏高原地区浅源地震的能矩比范围为5.03×10^(-6)—4.80×10^(-5),反映了不同断层上发生的地震能量释放过程存在差异。因此,为全面反映地震的震源特性有必要进行能量震级和矩震级的联合测定。

大震应急产品产出与服务平台建设

地震时空强参数已不能满足日益增长的最大限度减轻地震灾害的需求。作为国家级地震监测中心,中国地震台网中心联合各研究所和试点省局,建设了大震应急产品产出与服务平台。中强震发生后,该平台准实时产出除地震三要素之外的多类地震参数和图件,为及时准确评估受灾区域和灾损程度、揭示地震成因和致灾机理提供了丰富的基础资料,第一时间服务于救援决策和震情会商,显著提升了国家级地震监测中心的应急处置和践行减轻地震灾害的能力。

W-LCFST柱抗震性能及弯矩-曲率关系研究

目的研究L形宽肢钢管混凝土组合柱(W-LCFST柱)的抗震性能及弯矩-曲率计算方法,为工程应用提供设计依据。方法应用ABAQUS有限元分析软件建立考虑混凝土塑性损伤和延性金属损伤的有限元模型,在验证模型合理准确的基础上,对比分析等肢截面、不等肢截面、轴压比及钢管厚度等参数对抗震性能的影响;基于W-LCFST柱的参数化分析结果和破坏现象,提出柱截面屈服、峰值和极限状态弯矩-曲率的简化计算方法。结果等肢W-LCFST柱截面高宽比从4.0增至5.0时,初始刚度和正向峰值承载力分别提高了16.9%和35.0%;X向截面高宽比从4.0增至5.0时,初始刚度和正向峰值承载力分别提高了15.3%和31.8%;Y向截面高宽比从4.0增至5.0时,初始刚度提高了9.4%,正向峰值承载力提高了18.0%;轴压比从0.2增至0.6时,正向峰值承载力和延性分别降低了13.8%和6.8%;钢板厚度从6 mm增至10 mm后,峰值承载力和延性分别提升了25.7%和4.9%。结论有限元模拟、试验和理论计算的结果吻合度较好,误差在10%以内,验证了弯矩-曲率计算方法的准确性;等肢截面高宽比增加,增大X向肢长和增加钢管厚度使柱的承载力、初始刚度和延性有明显提高。