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.

Paleostress Reconstruction from 3D Seismic Data and Slip Tendency in the Northern Slope Area of the Bongor Basin, Southwestern Chad

Paleostress plays a significant role in controlling the formation, accumulation, and distribution of reservoirs, and this could be an important factor in controlling the production of hydrocarbons from the unconventional reservoirs. In this study, we will use 3D seismic reflection data to perform the slip-tendency-based stress inversion to determine the stress field in the basement of the northern slope area in the Bongor Basin. The dataset for this technique is easily available in the oil and gas companies. The stress inversion results from the basement of the northern slope area of Bongor basin show that the maximum principal stress axis (σ1) is oriented vertically, the intermediate principal stress axis (σ2) is oriented N18° and the minimum principal stress axis (σ3) is oriented N105°, and σ2/σ1 = 0.60 and σ3/σ1 = 0.29. The findings of this paper provide significant information to understand the fault reactivation at the critical stage of hydrocarbon accumulation and the regional tectonic evolution.

Efficiency and effectiveness of systems for the treatment of domestic wastewater based on subsurface flow constructed wetlands in Jarabacoa, Dominican Republic

Constructed wetlands(CwW)are well known nature-based systems for water treatment.This study evaluated the efficiency and effectiveness of seven domestic wastewater treatment systems based on horizontal flow CWs in Jarabacoa,the Dominican Republic.The results showed that the CWs were efficient in reducing the degree of contamination of wastewater to levels below the Dominican wastewater discharge standards for parameters such as the 5-day biochemical oxygen demand(BOD5)and chemical oxygen demand,but not for the removal of phosphorus and fecal coliforms.In addition,a horizontal flow subsurface wetland in the peri-urban area El Dorado was evaluated in terms of the performance of wastewater treatment in tropical climatic conditions.The concentrations of heavy metals,such as zinc,copper,chromium,and iron,were found to decrease in the effluent of the wetland,and the concentrations for nickel and manganese tended to increase.The levels of heavy metals in the effluent were lower than the limit values of the Dominican wastewater discharge standards.The construction cost of these facilities was around 200 USD per population equivalent,similar to the cost in other countries in the same region.This study suggested some solutions to the improved performance of CWs:selection of a microbial flora that guarantees the reduction of nitrates and nitrites to molecular nitrogen,use of endemic plants that bioaccumulate heavy metals,combination of constructed wetlands with filtration on activated carbon,and inclusion of water purification processes that allow to evaluate the reuse of treated water.

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.

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.

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.

基于Construct 3的射击类游戏设计与开发——以飞机大战为例

Construct 3作为小而精的H5跨平台游戏引擎,强大的交互性和用户友好性使其在2D游戏开发方面独具优势。借助Construct 3重构经典射击类游戏“飞机大战”,通过Event Sheet实现游戏中背景循环滚动、动画导入、角色创建与随机运动、游戏难度控制、游戏变量控制等主要功能。经过测试,游戏可打包发布成Web、移动设备和PC端.exe格式,游戏交互、特效运行流畅,游戏变量与难度控制处于合理范围。

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.

Proper orthogonal decomposition based seismic source wavefield reconstruction for finite element reverse time migration

The large storage requirement is a critical issue in cross-correlation imaging-condition based reverse time migration(RTM),because it requires the operation of the source and receiver wavefields at the same time.The boundary value method(BVM),based on the finite difference method(FDM),can be used to reconstruct the source wavefield in the reverse time propagation in the same way as the receiver wavefield,which can reduce the storage burden of the RTM data.Considering that the FDM cannot well handle models with discontinuous material properties and rough interfaces,we develop a source wavefield reconstruction strategy based on the finite element method(FEM),using proper orthogonal decomposition(POD)to enhance computational efficiency.In this method,we divide the whole time period into several segments,and construct the POD basis functions to get a reduced order model(ROM)for the source wavefield reconstruction in each segment.We show the corresponding quantitative analysis of the storage requirement of the POD-FEM.Numerical tests on the homogeneous model show the effectiveness of the proposed method,while the layered model and part of the Marmousi model tests indicate that the POD-FEM can keep an excellent balance between computational efficiency and memory usage compared with the full-stored method(FSM)and the BVM,and can be effectively applied in imaging.

Sparse Seismic Data Reconstruction Based on a Convolutional Neural Network Algorithm

At present,the acquisition of seismic data is developing toward high-precision and high-density methods.However,complex natural environments and cultural factors in many exploration areas cause difficulties in achieving uniform and intensive acquisition,which makes complete seismic data collection impossible.Therefore,data reconstruction is required in the processing link to ensure imaging accuracy.Deep learning,as a new field in rapid development,presents clear advantages in feature extraction and modeling.In this study,the convolutional neural network deep learning algorithm is applied to seismic data reconstruction.Based on the convolutional neural network algorithm and combined with the characteristics of seismic data acquisition,two training strategies of supervised and unsupervised learning are designed to reconstruct sparse acquisition seismic records.First,a supervised learning strategy is proposed for labeled data,wherein the complete seismic data are segmented as the input of the training set and are randomly sampled before each training,thereby increasing the number of samples and the richness of features.Second,an unsupervised learning strategy based on large samples is proposed for unlabeled data,and the rolling segmentation method is used to update(pseudo)labels and training parameters in the training process.Through the reconstruction test of simulated and actual data,the deep learning algorithm based on a convolutional neural network shows better reconstruction quality and higher accuracy than compressed sensing based on Curvelet transform.

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.

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.

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.