TransQuake: A transformer-based deep learning approach for seismic P-wave detection

Recent years,we have witnessed the increasing research interest in developing machine learning,especially deep learning which provides approaches for enhancing the performance of microearthquake detection.While considerable research efforts have been made in this direction,most of the state-of-the-art solutions are based on Convolutional Neural Network(CNN)structure,due to its remarkable capability of modeling local and static features.Indeed,the globally dynamic characteristics contained within time series data(i.e.,seismic waves),which cannot be fully captured by CNN-based models,have been largely ignored in previous studies.In this paper,we propose a novel deep learning approach,TransQuake,for seismic P-wave detection.The approach is based on the most advanced sequential model,namely Transformer.To be specific,TransQuake can exploit the STA/LTA algorithm for adapting the three-component structure of seismic waves as input,and take advantage of the multi-head attention mechanism for conducting explainable model learning.Extensive evaluations of the aftershocks following the 2008 Wenchuan MW 7.9 earthquake clearly demonstrates that TransQuake is able to achieve the best detection performance which excels the results obtained using other baselines.Meanwhile,experimental results also validate the interpretability of the results obtained by TransQuake,such as the attention distribution of seismic waves in different positions,and the analysis of the optimal relationship between coda wave and P-wave for noise identification.

Introduction to a recently released dataset entitled CSNCD:A Comprehensive Dataset of Chinese Seismic Network

China and its adjacent regions are heavily affected by earthquake disasters,which has led to strong demand and vigorous development of the Chinese Seismic Network(CSN).Since 2004,China Earthquake Networks Center(CENC)has been authorized by China Earthquake Administration(CEA)to take charge of CSN management,mainly observing product output and quality control.

Probabilistic analysis of tunnel face seismic stability in layered rock masses using Polynomial Chaos Kriging metamodel

Face stability is an essential issue in tunnel design and construction.Layered rock masses are typical and ubiquitous;uncertainties in rock properties always exist.In view of this,a comprehensive method,which combines the Upper bound Limit analysis of Tunnel face stability,the Polynomial Chaos Kriging,the Monte-Carlo Simulation and Analysis of Covariance method(ULT-PCK-MA),is proposed to investigate the seismic stability of tunnel faces.A two-dimensional analytical model of ULT is developed to evaluate the virtual support force based on the upper bound limit analysis.An efficient probabilistic analysis method PCK-MA based on the adaptive Polynomial Chaos Kriging metamodel is then implemented to investigate the parameter uncertainty effects.Ten input parameters,including geological strength indices,uniaxial compressive strengths and constants for three rock formations,and the horizontal seismic coefficients,are treated as random variables.The effects of these parameter uncertainties on the failure probability and sensitivity indices are discussed.In addition,the effects of weak layer position,the middle layer thickness and quality,the tunnel diameter,the parameters correlation,and the seismic loadings are investigated,respectively.The results show that the layer distributions significantly influence the tunnel face probabilistic stability,particularly when the weak rock is present in the bottom layer.The efficiency of the proposed ULT-PCK-MA is validated,which is expected to facilitate the engineering design and construction.

An adaptive physics-informed deep learning method for pore pressure prediction using seismic data

Accurate prediction of formation pore pressure is essential to predict fluid flow and manage hydrocarbon production in petroleum engineering.Recent deep learning technique has been receiving more interest due to the great potential to deal with pore pressure prediction.However,most of the traditional deep learning models are less efficient to address generalization problems.To fill this technical gap,in this work,we developed a new adaptive physics-informed deep learning model with high generalization capability to predict pore pressure values directly from seismic data.Specifically,the new model,named CGP-NN,consists of a novel parametric features extraction approach(1DCPP),a stacked multilayer gated recurrent model(multilayer GRU),and an adaptive physics-informed loss function.Through machine training,the developed model can automatically select the optimal physical model to constrain the results for each pore pressure prediction.The CGP-NN model has the best generalization when the physicsrelated metricλ=0.5.A hybrid approach combining Eaton and Bowers methods is also proposed to build machine-learnable labels for solving the problem of few labels.To validate the developed model and methodology,a case study on a complex reservoir in Tarim Basin was further performed to demonstrate the high accuracy on the pore pressure prediction of new wells along with the strong generalization ability.The adaptive physics-informed deep learning approach presented here has potential application in the prediction of pore pressures coupled with multiple genesis mechanisms using seismic data.

Seismic performance evaluation of water supply pipes installed in a full-scale RC frame structure based on a shaking table test

As an important part of nonstructural components,the seismic response of indoor water supply pipes deserves much attention.This paper presents shaking table test research on water supply pipes installed in a full-scale reinforced concrete(RC)frame structure.Different material pipes and different methods for penetrating the reinforced concrete floors are combined to evaluate the difference in seismic performance.Floor response spectra and pipe acceleration amplification factors based on test data are discussed and compared with code provisions.A seismic fragility study of displacement demand is conducted based on numerical simulation.The acceleration response and displacement response of different combinations are compared.The results show that the combination of different pipe materials and different passing-through methods can cause obvious differences in the seismic response of indoor riser pipes.

Time-domain dynamic constitutive model suitable for mucky soil site seismic response

Soil nonlinear behavior displays noticeable effects on the site seismic response.This study proposes a new functional expression of the skeleton curve to replace the hyperbolic skeleton curve.By integrating shear modulus and combining the dynamic skeleton curve and the damping degradation coefficient,the constitutive equation of the logarithmic dynamic skeleton can be obtained,which considers the damping effect in a soil dynamics problem.Based on the finite difference method and the multi-transmitting boundary condition,a 1D site seismic response analysis program called Soilresp1D has been developed herein and used to analyze the time-domain seismic response in three types of sites.At the same time,this study also provides numerical simulation results based on the hyperbolic constitutive model and the equivalent linear method.The results verify the rationality of the new soil dynamic constitutive model.It can analyze the mucky soil site nonlinear seismic response,reflecting the deformation characteristics and damping effect of the silty soil.The hysteresis loop area is more extensive,and the residual strain is evident.

Monitoring seismicity in the southern Sichuan Basin using a machine learning workflow

Monitoring seismicity in real time provides significant benefits for timely earthquake warning and analyses.In this study,we propose an automatic workflow based on machine learning(ML)to monitor seismicity in the southern Sichuan Basin of China.This workflow includes coherent event detection,phase picking,and earthquake location using three-component data from a seismic network.By combining Phase Net,we develop an ML-based earthquake location model called Phase Loc,to conduct real-time monitoring of the local seismicity.The approach allows us to use synthetic samples covering the entire study area to train Phase Loc,addressing the problems of insufficient data samples,imbalanced data distribution,and unreliable labels when training with observed data.We apply the trained model to observed data recorded in the southern Sichuan Basin,China,between September 2018 and March 2019.The results show that the average differences in latitude,longitude,and depth are 5.7 km,6.1 km,and 2 km,respectively,compared to the reference catalog.Phase Loc combines all available phase information to make fast and reliable predictions,even if only a few phases are detected and picked.The proposed workflow may help real-time seismic monitoring in other regions as well.

Major methods of seismic anisotropy

Seismic anisotropy reveals that seismic wave velocity, amplitude, and other physical properties show variations in different directions, which can be divided into lattice-preferred orientation(LPO) and shape-preferred orientation(SPO) according to its physical mechanisms. The main methods for studying seismic anisotropy include shearwave splitting analysis, P-wave travel time inversion and surface-wave tomography, etc. There are some differences and correlations among these methods. Seismic anisotropy is an important way to reveal the dynamic processes of crust-mantle evolution, and it is significant for monitoring crustal stress changes and improve seismic exploration studies. With the help of long-term observation, the application of machine learning techniques and combining inversion based on multiple phases would become potential developments in seismic anisotropy studies. This may improve the understanding of complex seismic anisotropic models, such as multiple layers anisotropy with an oblique axis of symmetry.

Dynamically triggered seismicity on a tectonic scale:A review

Earthquakes triggered by dynamic disturbances have been confirmed by numerous observations and experiments.In the past several decades,earthquake triggering has attracted increasing attention of scholars in relation to exploring the mechanism of earthquake triggering,earthquake prediction,and the desire to use the mechanism of earthquake triggering to reduce,prevent,or trigger earthquakes.Natural earthquakes and large‐scale explosions are the most common sources of dynamic disturbances that trigger earthquakes.In the past several decades,some models have been developed,including static,dynamic,quasi‐static,and other models.Some reviews have been published,but explosiontriggered seismicity was not included.In recent years,some new results on earthquake triggering have emerged.Therefore,this paper presents a new review to reflect the new results and include the content of explosion‐triggered earthquakes for the reference of scholars in this area.Instead of a complete review of the relevant literature,this paper primarily focuses on the main aspects of dynamic earthquake triggering on a tectonic scale and makes some suggestions on issues that need to be resolved in this area in the future.

Remnant Paleoproterozoic Subduction or Lithospheric Drip Initiation at the Yilgarn Craton Margin:Constraints from P-wave Tomography

The Paleoproterozoic was a critical time in whether modern-style plate tectonics had become globally dominant(e.g.,Wan et al.,2020).The Capricorn Orogen witnessed the assembly of the Pilbara and Yilgarn Cratons and an exotic microcontinent,the Glenburgh Terrane,to form the West Australia Craton(WAC)through two collisional orogenic events,the 2215–2145 Ma Ophthalmian and 2005–1950 Ma Glenburgh Orogenies(Johnson et al.,2013;Fig.1).Compared to other Proterozoic orogenic belts in Australia,the Capricorn Orogen preserves‘complete'opposing continental margin successions,together with intervening arc fragments associated with oceanic closure and foreland basins associated with collisional loading(Cawood et al.,2009).

Potential of the Arkhangelsk seismic network for European Arctic monitoring

The Arkhangelsk Seismic Network(ASN)of the N.Laverov Federal Center for Integrated Arctic Research of the Ural Branch of the Russian Academy of Sciences,founded in 2003,includes 10 permanent seismic stations located on the coasts of the White,Barents,and Kara Seas and on the Arctic archipelagos of Novaya Zemlya,Franz Josef Land,and Severnaya Zemlya.The network is registered with the International Federation of Digital Seismograph Networks and the International Seismological Center.We used not only ASN data to process earthquakes but also the waveforms of various international seismic stations.The 13,000 seismic events were registered using ASN data for 2012-2022,and for 5,500 of them,we determined the parameters of the earthquake epicenters from the European Arctic.The spatial distribution of epicenters shows that the ASN monitors not only the main seismically active zones but also weak seismicity on the shelf of the Barents and Kara Seas.The representative magnitude of ASN was ML,rep=3.5.The level of microseismic noise has seasonal variations that affect the registration capabilities of each station included in the ASN and the overall sensitivity of the network as a whole.In summer,the sensitivity of the ASN decreased owing to the increasing microseismic and ambient noises,whereas in winter,the sensitivity of the ASN increased significantly because of the decrease.

The 2023 Turkey earthquake doublet: Earthquake relocation, seismic tomography, and stress field inversion

On February 6,2023,two earthquakes with magnitudes of M_(W) 7.8 and M_(W) 7.5 struck southeastern Turkey,causing significant casualties and economic losses.These seismic events occurred along the East Anatolian Fault Zone,a convergent boundary between the Arabian Plate and the Anatolian Subplate.In this study,we analyze the M_(W) 7.8 and M_(W) 7.5 earthquakes by comparing their aftershock relocations,tomographic images,and stress field inversions.The earthquakes were localized in the upper crust and exhibited steep dip angles.Furthermore,the aftershocks occurred either close to the boundaries of low and high P-wave velocity anomaly zones or within the low P-wave velocity anomaly zones.The East Anatolia Fault,associated with the M_(W) 7.8 earthquake,and the SürgüFault,related to the M_(W) 7.5 earthquake,predominantly experienced shear stress.However,their western sections experienced a combination of strike-slip and tensile stresses in addition to shear stress.The ruptures of the M_(W) 7.8 and M_(W) 7.5 earthquakes appear to have bridged a seismic gap that had seen sparse seismicity over the past 200 years prior to the 2023 Turkey earthquake sequence.

Resilience-incorporated seismic risk assessment of precast concrete frames with“dry”connections

A resilience-incorporated risk assessment framework is proposed and demonstrated in this study to manifest the advantageous seismic resilience of precast concrete frame(PCF)structures with“dry”connections in terms of their low damage and rapid recovery.The framework integrates various uncertainties in the seismic hazard,fragility,capacity,demand,loss functions,and post-earthquake recovery.In this study,the PCF structures are distinguished from ordinary reinforced concrete frame(RCF)structures by characterizing multiple limit states for the PCF based on its unique damage mechanisms.Accordingly,probabilistic story-wise pushover analyses are performed to yield story-wise capacities for the predefined limit states.In the seismic resilience analysis,a step-wise recovery model is proposed to idealize the functionality recovery process,with separate considerations of the repair and non-repair events.The recovery model leverages the economic loss and downtime to delineate the stochastic post-earthquake recovery curves for the resilience loss estimation.As such,contingencies in the probabilistic post-earthquake repairs are incorporated and the empirical judgments on the recovery parameters are largely circumvented.The proposed framework is demonstrated through a comparative study between two“dry”connected PCFs and one RCF designed as alternative structural systems for a prototype building.The results from the risk quantification indicate that the PCFs show reduced loss hazards and lower expected losses relative to the RCF.Particularly,the PCF equipped with energy dissipation devices at the“dry”connections largely reduces the expected economic loss,downtime,and resilience loss by 29%,56%,and 60%,respectively,compared to the RCF.

Fracture network characterisation of the Balmuccia peridotite using drone-based photogrammetry,implications for active-seismic site survey for scientific drilling

The presence of discontinuities(e.g.faults,fractures,veins,layering)in crystalline rocks can be challenging for seismic interpretations because the wide range of their size,orientation,and intensity,which controls the mechanical properties of the rock and elastic wave propagation,resulting in equally varying seismic responses at different scales.The geometrical characterisation of adjacent outcrop discontinuity networks allows a better understanding of the nature of the subsurface rocks and aids seismic interpretation.In this study,we characterise the discontinuity network of the Balmuccia peridotite(BP)in the IvreaeVerbano Zone(IVZ),northwestern Italy.This geological body is the focus of the Drilling the Ivrea eVerbano zonE(DIVE),an international continental scientific drilling project,and two active seismic surveys,SEismic imaging of the Ivrea ZonE(SEIZE)and high-resolution SEIZE(Hi-SEIZE),which aim to resolve the subsurface structure of the DIVE drilling target through high-resolution seismic imaging.For fracture characterisation,we developed two drone-based digital outcrop models(DOMs)at two different resolutions(10^(-3)-10 m and 10^(-1)-10^(3)m),which allowed us to quantitatively characterise the orientation,size,and intensity of the main rock discontinuities.These properties affect the seismic velocity and consequently the interpretation of the seismic data.We found that(i)the outcropping BP discontinuity network is represented by three more sets of fractures with respect to those reported in the literature;(ii)the discontinuity sizes follow a power-law distribution,indicating similarity across scales,and(iii)discontinuity intensity is not uniformly distributed along the outcrop.Our results help to explain the seismic behaviour of the BP detected by the SEIZE survey,suggesting that the low P-wave velocities observed can be related to the discontinuity network,and provide the basic topological parameters(orientation,density,distribution,and aperture)of the fracture network unique to the BP.These,in turn,can be used for interpretation of the Hi-SEIZE seismic survey and forward modelling of the seismic response.

A reweighted damped singular spectrum analysis method for robust seismic noise suppression

(Multichannel)Singular spectrum analysis is considered as one of the most effective methods for seismic incoherent noise suppression.It utilizes the low-rank feature of seismic signal and regards the noise suppression as a low-rank reconstruction problem.However,in some cases the seismic geophones receive some erratic disturbances and the amplitudes are dramatically larger than other receivers.The presence of this kind of noise,called erratic noise,makes singular spectrum analysis(SSA)reconstruction unstable and has undesirable effects on the final results.We robustify the low-rank reconstruction of seismic data by a reweighted damped SSA(RD-SSA)method.It incorporates the damped SSA,an improved version of SSA,into a reweighted framework.The damping operator is used to weaken the artificial disturbance introduced by the low-rank projection of both erratic and random noise.The central idea of the RD-SSA method is to iteratively approximate the observed data with the quadratic norm for the first iteration and the Tukeys bisquare norm for the rest iterations.The RD-SSA method can suppress seismic incoherent noise and keep the reconstruction process robust to the erratic disturbance.The feasibility of RD-SSA is validated via both synthetic and field data examples.

Detailed seismic zoning of the East Kazakhstan region in the Republic of Kazakhstan

Kazakhstan is currently drafting new construction regulations that comply with the major provisions of the Eurocodes.Such regulations are created on the basis of seismic zoning maps of various degrees of detail,developed by our Institute of Seismology using a new methodological approach for Kazakhstan.The article is about creating the first normative map of the Detailed Seismic Zoning on a probabilistic foundation for the Republic of Kazakhstan’s East Kazakhstan region.We carried out the probabilistic assessment of seismic hazard using a methodology consistent with the main provisions of Eurocode 8and updated compared with that used in developing maps of Kazakhstan’s General Seismic Zoning and seismic microzoning of Almaty.The most thorough and current data accessible for the area under consideration were combined with contemporary analytical techniques.Updates have been done to not only the databases being used but also the way seismic sources were shown,including active faults now.On a scale of 1:1000000,precise seismic zoning maps of the East Kazakhstan region were created for two probabilities of exceedance:10%and 2%in 50 years in terms of peak ground accelerations and macroseismic intensities.The obtained seismic hazard distribution is generally consistent with the General Seismic Zoning of Kazakhstan’s previous findings.However,because active faults were included and a thoroughly revised catalog was used,there are more pronounced zones of increased danger along the fault in the western part of the region.In the west of the territory,acceleration values also increased due to a more accurate consideration of seismotectonic conditions.Zoning maps are the basis for developing new state building regulations of the Republic of Kazakhstan.

Hcable for Time-Lapse Seismic Monitoring of Marine Carbon Capture and Storage

To ensure project safety and secure public support, an integrated and comprehensive monitoring program is needed within a carbon capture and storage(CCS) project. Monitoring can be done using many well-established techniques from various fields, and the seismic method proves to be the crucial one. This method is widely used to determine the CO_(2) distribution, image the plume development, and quantitatively estimate the concentration. Because both the CO_(2) distribution and the potential migration pathway can be spatially small scale, high resolution for seismic imaging is demanded. However, obtaining a high-resolution image of a subsurface structure in marine settings is difficult. Herein, we introduce the novel Hcable(Harrow-like cable system) technique, which may be applied to offshore CCS monitoring. This technique uses a highfrequency source(the dominant frequency>100 Hz) to generate seismic waves and a combination of a long cable and several short streamers to receive seismic waves. Ultrahigh-frequency seismic images are achieved through the processing of Hcable seismic data. Hcable is then applied in a case study to demonstrate its detailed characterization for small-scale structures. This work reveals that Hcable is a promising tool for timelapse seismic monitoring of oceanic CCS.

On electro-acoustic characteristics of a marine broadband sparker for seismic exploration

The construction of major marine infrastructure projects and the exploration and development of deep-sea mineral resources require fine imaging of seabed strata and structures.The highresolution marine seismic exploration based on a high broadband sparker source is an important approach to reveal seabed stratum and reservoir structure,and identify geohazard.To optimize the performance of sparker seismic source,we investigated the electro-acoustic characteristics of spark discharge under conditions of different charging voltages and electrode numbers.Results show that the sound source level increased with the increase of the charging voltage,whereas the main frequency decreased when the charging voltage increases.In addition,it was found that the charging capacitance had more obvious influence on the main frequency than the sound source level did.Although the load energy decreased with increasing electrode number,the sound source level still increased but the main frequency decreased.Meanwhile,the primary to bubble(P/B)ratio increased with the increase of the electrode number.To gain a deeper insight into the electro-acoustic characteristics,we investigate the relationship between sound source level and power peak,from which a good correlation was observed.A more practical statistical analysis on the rise rate of current was processed,and a perfect logarithmic function was derived.Furthermore,we found that the main frequency was most possibly subjected to the electrical energy,especially the charging energy per electrode.The results indicate that the charging energy per electrode less than 10 J could increase the main frequency to above 300 Hz.At last,the main frequency could be reduced to 20 Hz when the charging energy of a single-electrode discharge was enhanced to over 4 kJ.This study shall be helpful in developing a sparker seismic source and improving the performance for marine engineering exploration and geohazard assessment.

The seismicity in the middle section of the Altyn Tagh Fault system revealed by a dense nodal seismic array

The left-lateral Altyn Tagh Fault(ATF) system is the northern boundary of the Qinghai-Xizang Plateau, separating the Tarim Basin and the Qaidam Basin. The middle section of ATF has not recorded any large earthquakes since1598 AD, so the potential seismic hazard is unclear. We develope an earthquake catalog using continuous waveform data recorded by the Tarim-Altyn-Qaidam dense nodal seismic array from September 17 to November23, 2021 in the middle section of ATF. With the machine learning-based picker, phase association, location, match and locate workflow, we detecte 233 earthquakes with M_L-1–3, far more than 6 earthquakes in the routine catalog. Combining with focal mechanism solutions and the local fault structure, we find that seismic events are clustered along the ATF with strike-slip focal mechanisms and on the southern secondary faults with thrusting focal mechanisms. This overall seismic activity in the middle section of the ATF might be due to the northeastward transpressional motion of the Qinghai-Xizang Plateau block at the western margin of the Qaidam Basin.

Erratum to:A review of the wave gradiometry method for seismic imaging

Earthq Sci 2023(36):254-281 Doi:10.1016/j.eqs.2023.04.002 In the original version of this article,the important funding source was inadvertently omitted in Acknowledgement section.The corrected one is as follows.