Constraints on triggered seismicity and its control on permeability evolution

Triggered seismicity is a key hazard where fluids are injected or withdrawn from the subsurface and may impact permeability. Understanding the mechanisms that control fluid injection-triggered seismicity allows its mitigation. Key controls on seismicity are defined in terms of fault and fracture strength, second-order frictional response and stability, and competing fluid-driven mechanisms for arrest. We desire to constrain maximum event magnitudes in triggered earthquakes by relating pre-existing critical stresses to fluid injection volume to explain why some recorded events are significantly larger than anticipated seismic moment thresholds. This formalism is consistent with several uncharacteristically large fluid injection-triggered earthquakes. Such methods of reactivating fractures and faults by hydraulic stimulation in shear or tensile fracturing are routinely used to create permeability in the subsurface. Microearthquakes (MEQs) generated by such stimulations can be used to diagnose permeability evolution. Although high-fidelity data sets are scarce, the EGS-Collab and Utah FORGE hydraulic stimulation field demonstration projects provide high-fidelity data sets that concurrently track permeability evolution and triggered seismicity. Machine learning deciphers the principal features of MEQs and the resulting permeability evolution that best track permeability changes – with transfer learning methods allowing robust predictions across multiple eological settings. Changes in permeability at reactivated fractures in both shear and extensional modes suggest that permeability change (Δk) scales with the seismic moment (M) of individual MEQs as Δk∝M. This scaling relation is exact at early times but degrades with successive MEQs, but provides a method for characterizing crustal permeability evolution using MEQs, alone. Importantly, we quantify for the first time the role of prestress in defining the elevated magnitude and seismic moment of fluid injection-triggered events, and demonstrate that such MEQs can also be used as diagnostic in quantifying permeability evolution in the crust.

From earthquake resistance structure to earthquake resilience city–urban seismic resilience assessment

Earthquakes pose significant perils to the built environment in urban areas.To avert the calamitous aftermath of earthquakes,it is imperative to construct seismic resilient cities.Due to the intricacy of the concept of urban seismic resilience(USR),its assessment is a large-scale system engineering issue.The assessment of USR should be based on the notion of urban seismic capacity(USC)assessment,which includes casualties,economic loss,and recovery time as criteria.Functionality loss is also included in the assessment of USR in addition to these criteria.The assessment indicator system comprising five dimensions(building and lifeline infrastructure,environment,society,economy,and institution)and 20 indicators has been devised to quantify USR.The analytical hierarchy process(AHP)is utilized to compute the weights of the criteria,dimensions,and indicators in the urban seismic resilience assessment(USRA)indicator system.When the necessary data for a city are obtainable,the seismic resilience of that city can be assessed using this framework.To illustrate the proposed methodology,a moderate-sized city in China was selected as a case study.The assessment results indicate a high level of USR,suggesting that the city possesses strong capabilities to withstand and recover from potential future earthquakes.

Seismic design of variable cross-section damped steel support frame for post equipment

High-voltage electrical post equipment is generally installed on steel supports,which amplifies the seismic inputs and degrades the seismic performance of equipment.This study proposed a variable cross-section damped steel support frame(VCDFS)with viscous dampers to reduce seismic responses of both tall and low-rise electrical equipment.The VCDFS contains a trapezoidal damper layer to generate rocking motions,enabling the diagonal viscous dampers to dissipate seismic inputs.A theoretical model of post equipment with VCDFS is established,and an optimal design procedure is proposed.The analysis shows that the remaining static stiffness ratio λ_(k) is the key parameter that determines the effectiveness of the VCDFS.The VCDFS reduces the average displacement and stress response of a post insulator by 39.4%and 44.6%,respectively,together with a significant decrease in the dynamic amplification factor.Therefore,it is recommended to use the VCDFS instead of the conventional latticed-steel frame in earthquake zones.

Influence of a non-free field subsurface seismic history on the liquefaction resistance of sandy soils

Seismic-induced liquefaction of sandy soils can fail foundations in the vicinity of buildings.To investigate the effect of a non-free field subsurface seismic history on the ability of saturated sandy soils to resist liquefaction,four shaking events with different accelerations were input to the sandy soils in the non-free-field.The results of the study revealed that:(1)Shallow soils that are not free-field undergo acceleration amplification effects after being subjected to seismic loading.(2)Building overburden pressure reduces the sensitivity of the shallow soils directly below in small and moderate earthquakes,which are more prone to rearranging and forming unstable structures under strong seismic effects.The excess pore pressure response on the load side resembles that of a free site,with the depth range of the liquefaction strength of soils affected by the seismic history,increasing progressively as input seismic intensity increases.(3)After experiencing earthquakes of different intensities,the excess pore pressure directly below the building overburden pressure at 0.1 m and 0.2 m is greater than that at the side.At the same time,the side of the building structure is more prone to liquefaction than the soil directly below it.

Seismic isolation design and resilience improvement of railway station considering the influence of near-fault pulse-like ground motions

To improve the resilience of railway stations,a typical station was selected as the research object,and an isolation design was introduced.Twenty-four groups of near-fault pulse-like ground motions were selected.The seismic resilience of the no-isolation railway stations(NIRS)and the isolation railway stations(IRS)were compared to provide a numerical result of the improvement in resilience.The results show that in the station isolation design,the station's functional requirements and structural characteristics should be considered and the appropriate placement of isolation bearings is under the waiting room.Under the action of a rare earthquake,the repair cost,repair time,rate of harm and death of the IRS were decreased by 8.04 million,18.30 days,6.93×10^(-3)and 1.21×10^(-3),respectively,when compared to the NIRS.The IRS received a seismic resilience grade of three-stars and the NIRS only one-star,indicating that rational isolation design improves the seismic resilience of stations.Thus,for the design of stations close to earthquake faults,it is suggested to utilize appropriate isolation techniques to improve their seismic resilience.

A logistic-Lasso-regression-based seismic fragility analysis method for electrical equipment considering structural and seismic parameter uncertainty

Damage to electrical equipment in an earthquake can lead to power outage of power systems.Seismic fragility analysis is a common method to assess the seismic reliability of electrical equipment.To further guarantee the efficiency of analysis,multi-source uncertainties including the structure itself and seismic excitation need to be considered.A method for seismic fragility analysis that reflects structural and seismic parameter uncertainty was developed in this study.The proposed method used a random sampling method based on Latin hypercube sampling(LHS)to account for the structure parameter uncertainty and the group structure characteristics of electrical equipment.Then,logistic Lasso regression(LLR)was used to find the seismic fragility surface based on double ground motion intensity measures(IM).The seismic fragility based on the finite element model of an±1000 kV main transformer(UHVMT)was analyzed using the proposed method.The results show that the seismic fragility function obtained by this method can be used to construct the relationship between the uncertainty parameters and the failure probability.The seismic fragility surface did not only provide the probabilities of seismic damage states under different IMs,but also had better stability than the fragility curve.Furthermore,the sensitivity analysis of the structural parameters revealed that the elastic module of the bushing and the height of the high-voltage bushing may have a greater influence.

Cretaceous Sandbody Characters at Shallow-Water Lake Delta Front and the Sedimentary Dynamic Process Analysis in Songliao Basin,China

Based on the abundant information from drilling, cores, and logging, the influence of topography, size of rivers and lakes, climate changes and the lake level＇s fluctuation on the sandbodies at shallow-water delta front are systematically summarized and the sedimentary dynamic processes are analyzed. The interwell communication among the sandbodies and their planar distribution revealed from the hydrodynamic features of the development wells are integrated during the analysis. The fundamental requirements for the development of the shallow-water delta included flat topography and uniform subsiding rate. The delta plain was connected smoothly with the wide delta front and predelta, without the three-fold structure of topset, foreset, and bottomset as defined in the Gilbert Delta Model. Because of the weak fluvial effect and the lake energy is strong, the small and scattered shallow-water delta is destroyed by the scouring-backwashing, coastal current, and lake wave, resulting in the coastal sheet deposition. As the fluvial effect became stronger and the lake energy became weaker, the shape of the shallow-water deltas transferred from sheets to lumps and then branches.

Solution of shallow-water equations using least-squares finite-element method

A least-squares finite-element method （LSFEM） for the non-conservative shallow-water equations is presented. The model is capable of handling complex topography, steady and unsteady flows, subcritical and supercritical flows, and flows with smooth and sharp gradient changes. Advantages of the model include： （1） sources terms, such as the bottom slope, surface stresses and bed frictions, can be treated easily without any special treatment; （2） upwind scheme is no needed; （3） a single approximating space can be used for all variables, and its choice of approximating space is not subject to the Ladyzhenskaya-Babuska-Brezzi （LBB） condition; and （4） the resulting system of equations is symmetric and positive-definite （SPD） which can be solved efficiently with the preconditioned conjugate gradient method. The model is verified with flow over a bump, tide induced flow, and dam-break. Computed results are compared with analytic solutions or other numerical results, and show the model is conservative and accurate. The model is then used to simulate flow past a circular cylinder. Important flow charac-teristics, such as variation of water surface around the cylinder and vortex shedding behind the cylinder are investigated. Computed results compare well with experiment data and other numerical results.

Least-squares finite-element method for shallow-water equations with source terms

Numerical solution of shallow-water equations （SWE） has been a challenging task because of its nonlinear hyperbolic nature, admitting discontinuous solution, and the need to satisfy the C-property. The presence of source terms in momentum equations, such as the bottom slope and friction of bed, compounds the difficulties further. In this paper, a least-squares finite-element method for the space discretization and θ-method for the time integration is developed for the 2D non-conservative SWE including the source terms. Advantages of the method include： the source terms can be approximated easily with interpolation functions, no upwind scheme is needed, as well as the resulting system equations is symmetric and positive-definite, therefore, can be solved efficiently with the conjugate gradient method. The method is applied to steady and unsteady flows, subcritical and transcritical flow over a bump, 1D and 2D circular dam-break, wave past a circular cylinder, as well as wave past a hump. Computed results show good C-property, conservation property and compare well with exact solutions and other numerical results for flows with weak and mild gradient changes, but lead to inaccurate predictions for flows with strong gradient changes and discontinuities.

An Anti-multipath Frequency Hopped Communication Technique in Shallow-water Acoustic Channels

This paper introduces a frequency-hopped (FH) communication system to anti-intersymbol interferences (ISI) caused by the multipath propagation in shallow-water acoustic channels, and uses high-speed digital signal processor (DSP) and serial ADC (MAX121) chip to demodulate received signal efficiently based Fast Fourier Transform (FFT) algorithm. The field experimental results show: a data rate of 1Kbit/s with the bit error rates on the order of 10 -4 is demonstrated at 2000 m in the shallow-water acoustic channel of Xiamen harbor, and the key techniques of the system is analyzed in the paper.

Parallel Computation of Shallow-water Model on Workstations Cluster

ＰａｒａｌｌｅｌＣｏｍｐｕｔａｔｉｏｎｏｆＳｈａｌｌｏｗ－ｗａｔｅｒＭｏｄｅｌｏｎＷｏｒｋｓｔａｔｉｏｎｓＣｌｕｓｔｅｒＳｏｎｇＪｕｎｑｉａｎｇ；ＳｕｎＡｎ－ｃｌａｎｇ，；ＬｉＸｉａｏｍｅｉ（ｅｐａｒｔｍｅｎｔＯｆＣＯ，ｍｐ＇ｕｔｅｒＳｃｉｅｎｃｅ，...

Stratal Slice Recognition of Thin Shallow-Water Delta Sandbodies in the Songliao Basin

Objective China＇s petroleum exploration has entered a new stage of finding deeply buried thin sandbodies lbr the abundant oil resources they contain. Here thin sandbodies refer to those less than 10 m in thickness, or even less than 1-2 m. It is difficult to depict thin-layer sandbodies of different genetic types using conventional core, well logging and seismic data due to their limited vertical resolution in petroliferous basins. However, seismic sedimentology provides a new research method especially tbr thin sandbody interpretation, i.e., validating interpreted sedimentary sandbodies from 3D seismic data based on horizontal resolution, stratal slice and seismic geomorphology interpretation. At present, a series of studies on seismic sedimentology in North America marine basins and elsewhere have been completed successfully and are relevant to the exploration and development of oil and gas fields.

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.

A Similarity Technique for Solving Two-Layer Shallow-Water Equations

This paper is devoted to the analysis of the two-layer shallow-water equations representing gravity currents. A similarity technique which is the characteristic function method is applied for this study. The application of the characteristic function method makes it possible to obtain the similarity forms depending on a group of infinitesimal transformations. Thus, the number of independent variables is reduced by one and the governing partial differential equations with the auxiliary conditions reduce to a system of ordinary differential equations with the appropriate auxiliary conditions. Numeric solutions are presented and discussed.

A Method for Reducing Ocean Wave-Induced Magnetic Noises in Shallow-Water MT Data Using a Complex Adaptive Filter

In shallow-water areas,the marine magnetotelluric(MT)method faces a challenge in the investigation of seabed conductivity structures due to electrical and magnetic noises induced by ocean waves,which seriously contaminate MT data.Ocean waves can affect electric and magnetic fields to different extents.In general,their influence on magnetic fields is considerably greater than that on electric fields.In this paper,a complex adaptive filter is adopted to reduce wave-induced magnetic noises in the frequency domain.The processing results of synthetic and measured MT data indicate that the proposed method can effectively reduce wave-induced magnetic noises and provide reliable apparent resistivity and phase 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.

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.

Transformation of Nonlinear Surface Gravity Waves under Shallow-Water Conditions

This article describes transformation of nonlinear surface gravity waves under shallow-water conditions with the aid of the suggested semigraphical method. There are given profiles of surface gravity waves on the crests steepening stages, their leading edges steepening. There are discussed the spectral component influence on the transformation of surface wave profile.

Numerical Simulation of Nonlinear Surface Gravity Waves Transformation under Shallow-Water Conditions

This work considers the problems of numerical simulation of non-linear surface gravity waves transformation under shallow bay conditions. The discrete model is built from non-linear shallow-water equations. Are resulted boundary and initial conditions. The method of splitting into physical processes receives system from three equations. Then we define the approximation order and investigate stability conditions of the discrete model. The sweep method was used to calculate the system of equations. This work presents surface gravity wave profiles for different propagation phases.

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.