The 2025 M_(w)7.7 Myanmar earthquake highlighted the challenge of near-fault seismic intensity field reconstruction due to sparse seismic networks.To address this limitation,a framework was proposed integrating seismi...The 2025 M_(w)7.7 Myanmar earthquake highlighted the challenge of near-fault seismic intensity field reconstruction due to sparse seismic networks.To address this limitation,a framework was proposed integrating seismic wave simulation with a data-constrained finite-fault rupture model.The constraint is implemented by identifying the optimal ground motion models(GMMs)through a scoring system that selects the best-fit GMMs to mid-and far-field China Earthquake Networks Center(CENC)seismic network data;and applying the optimal GMMs to refine the rupture model parameters for near-fault intensity field simulation.The simulated near-fault seismic intensity field reproduces seismic intensities collected from Myanmar’s sparse seismic network and concentrated in≥Ⅷintensity zones within 50 km of the projected fault plane;and identifies abnormal intensity regions exhibiting≥Ⅹintensity along the Meiktila-Naypyidaw corridor and near Shwebo that are attributed to soft soil amplification effects and near-fault directivity.This framework can also be applied to post-earthquake assessments in other similar regions.展开更多
We develop and implement a Stochastic Discrete Event Simulation(SDES)algorithm to model the housing re-covery trajectory after an extreme event.The algorithm models discrete events and their underlying uncertainties i...We develop and implement a Stochastic Discrete Event Simulation(SDES)algorithm to model the housing re-covery trajectory after an extreme event.The algorithm models discrete events and their underlying uncertainties in each construction phase.Specifically,the algorithm is developed for the Government Assisted Owner Driven(GAOD)reconstruction system to simulate long-term recovery trajectory.SDES,as a flexible modeling approach,can simulate any housing recovery scenario that follows phased reconstruction.The 2015 M 7.8 Gorkha earthquake sequence in Nepal is considered the extreme event,with 796,245 buildings requiring reconstruction.We present some recovery trajectories from severely hit,crisis hit,and earthquake hit parishes,comparing them with the actual reconstruction progress.We also assess quality and improvement of reconstructed buildings using seismic fragility functions,compared to pre-earthquake constructions.Housing recovery uncertainties are dissected in relation to reconstruction pace.We conclude that the vast majority of the reconstructed buildings followed the Build Back Better(BBB)approach and missed the opportunity to pursue the Build Back Resilient(BBR)approach due to multifaceted challenges ranging from unclear policies to economic constraints.We critically assess the GAOD vs Owner Driven(OD)recovery framework and conclude that insurance-supported and technically assisted OD approach could be the most suitable model for post extreme event housing recovery.展开更多
Gas extraction practice has been proven for the clear majority of coal mines in China to be unfavorable using drill holes in the coal seam. Rather, mining-induced fractures in the goaf should be utilized for gas extra...Gas extraction practice has been proven for the clear majority of coal mines in China to be unfavorable using drill holes in the coal seam. Rather, mining-induced fractures in the goaf should be utilized for gas extraction. To study gas migration in mining-induced fractures, one mining face of 10 th Mine in Pingdingshan Coalmine Group in Henan, China, has been selected as the case study for this work. By establishing the mathematical model of gas migration under the influence of coal seam mining, discrete element software UDEC and Multiphysics software COMSOL are employed to model gas migration in mining-induced fractures above the goaf. The results show that as the working face advances, the goaf overburden gradually forms a mining-induced fracture network in the shape of a trapezoid, the size of which increases with the distance of coal face advance. Compared with gas migration in the overburden matrix, the gas flow in the fracture network due to mining is far greater. The largest mining-induced fracture is located at the upper end of the trapezoidal zone, which results in the largest gas flux in the network. When drilling for gas extraction in a mining-induced fracture field, the gas concentration is reduced in the whole region during the process of gas drainage, and the rate of gas concentration drops faster in the fractured zone. It is shown that with gas drainage, the gas flow velocity in the mininginduced fracture network is faster.展开更多
To reduce the spatial simulation error generated by the finite difference method,previous researchers compute the optimal finite-difference weights always by minimizing the error of spatial dispersion relation.However...To reduce the spatial simulation error generated by the finite difference method,previous researchers compute the optimal finite-difference weights always by minimizing the error of spatial dispersion relation.However,we prove that the spatial simulation error of the finite difference method is associated with the dot product of the spatial dispersion relation of the finite-difference weights and the spectrum of the seismic wavefield.Based on the dot product relation,we construct a L_(2) norm cost function to minimize spatial simulation error.For solving this optimization problem,the seismic wavefield infor-mation in wavenumber region is necessary.Nevertheless,the seismic wavefield is generally obtained by costly forward modeling techniques.To reduce the computational cost,we substitute the spectrum of the seismic wavelet for the spectrum of the seismic wavefield,as the seismic wavelet plays a key role in determining the seismic wavefield.In solving the optimization problem,we design an exhaustive search method to obtain the solution of the L_(2) norm optimization problem.After solving the optimization problem,we are able to achieve the finite-difference weights that minimize spatial simulation error.In theoretical error analyses,the finite-difference weights from the proposed method can output more accurate simulation results compared to those from previous optimization algorithms.Furthermore,we validate our method through numerical tests with synthetic models,which encompass homogenous/inhomogeneous media as well as isotropic and anisotropic media.展开更多
The March 28,2025 Myanmar earthquake generated ground shaking that was perceptible throughout Myanmar and adjacent regions.This study simulated three-component ground motions across the affected region using an improv...The March 28,2025 Myanmar earthquake generated ground shaking that was perceptible throughout Myanmar and adjacent regions.This study simulated three-component ground motions across the affected region using an improved stochastic finite-fault method to systematically assess seismic impacts.Observed near-field recordings at MM.NGU station was used to determine the reliability of the theoretically derived stress drop as input for simulation.Far-field recordings constrained the frequency-dependent S-wave quality factors(Q(f)=283.305f^(0.588))for anelastic attenuation modeling.Comparisons of peak accelerations between simulation and empirical ground-motion models showed good agreement at moderate-to-large distances.However,lower near-fault simulations indicate a weaker-than-average source effect.Analysis of simulated instrumental seismic intensity revealed key patterns.Maximum intensity(Ⅹ)occurred in isolated patches within the ruptured fault projection,correlating with shallow high-slip areas.TheⅨ-intensity zone formed a north-south elongated band centered on fault projection.Significant asymmetry inⅧ-intensity distribution perpendicular to the fault strike was observed,with a wider western extension attributed to lower shear-wave velocities west of the fault.Supershear rupture behavior enhanced ground motions,expanding intensity ranges by~20%compared to sub-shear rupture.This study reveals the integrated effects of fault geometry,slip spatial distribution,rupture velocity,and site condition in governing ground motion patterns.展开更多
In this study, the spatial distributions of stress and fracture fields for three typical underground coal mining layouts, Le, non-pillar mining （NM）, top-coal caving mining （TCM） and protective coal-seam mining （...In this study, the spatial distributions of stress and fracture fields for three typical underground coal mining layouts, Le, non-pillar mining （NM）, top-coal caving mining （TCM） and protective coal-seam mining （PCM）, are modeled using discrete element software UDEC, The numerical results show that different mining layouts can lead to different mining-induced stress fields, resulting in diverse fracture fields, For the PCM, the mining influenced area in front of the mining faces is the largest, and the stress concentration factor in front of the mining faces is the lowest, The spatial shapes of the mining-induced fracture fields under NM, TCM and PCM differ, and they are characterized by trapezoidal, triangular and tower shapes, respectively, The fractal dimensions of mining-induced fractures of the three mining layouts decrease in the order of PCM, TCM and NM, It is also shown that the PCM can result in a better gas control effect in coal mines with high outburst potential, The numerical results are expected to provide a basis for understanding of mining-induced gas seepage fields and provide a reference for high- efficiency coal mining,展开更多
Seismic quantitative reservoir simulations and characterizations have played a vital role in exploring stratigraphic traps,such as lateaggradational lowstands prograding wedge systems(LPWS)within lowstands systems tra...Seismic quantitative reservoir simulations and characterizations have played a vital role in exploring stratigraphic traps,such as lateaggradational lowstands prograding wedge systems(LPWS)within lowstands systems tracts(LST).However,seismic data acquisition operations are always dominated by exceptional seismic coherent noise events,e.g.,multiples,which reduce the signal strengths of the sourcegenerated incident seismic waves within vertically and laterally heterogeneous earth systems.Hence,these noise events create hurdles in predicting paleo-depositional impedance(PDI),paleo-thickness(PTS),paleo-dense fractured networks,erosional and depositional zones,faultcontrolled migrations,and types of seismic reflection configurations(SRFC),which are key elements in developing stratigraphic pinch-out traps.This research utilizes the state-of-the-art technologies of spectral wavelet-based instantaneous time-frequency analysis and seismic waveform frequency-controlled porosity-constrained static reservoir simulation(FDPVS)tools to quantify the LPWS inside the Onshore Basin,Pakistan.The use of conventional amplitude-based seismic attributes,such as the average energy,remained a better tool for deciphering the overall geological architecture of the LPWS.Conventional FDPVS realizations resolved a PDI of−1.391 gm./c.c.^(*)m/s to−0.97 gm./c.c.^(*)m/s for LPWS with PTS of 12 and 20 m,respectively.A 0.9 km lateral extent of paleo-dense fractured networks(PDFN)with a strong linear regression R^(2)=0.93 was also resolved.Average energy attribute-based instantaneous frequency FDPVS realizations enabled the imaging of parallel-toprograding SRFC with resolved magnitudes of−0.259 gm./c.c.^(*)m/s for PDI,20 m for PTS,and 0.73 km for PDFN with linear regression transforms at R^(2)=0.92,which indicates the deposition of onlap fill facies inside the LPWS during extensive sea-level fall.These realizations have also resolved frequency-controlled fault migrations on 27-Hz spectral waveform-based amplitude plots with 2.174 gm./c.c.^(*)m/s PDI for conduit fault systems and 27-Hz with 0.585 gm./c.c.^(*)m/s PDI for sealing fault systems.All these structural configurations are completely sealed up by transgressive seals of transgressive systems tracts and,hence,developed into pure stratigraphic-based oil and gas plays.This research has strong implications for side-tracking drilling locations and provides an analogue for basins with similar geology and stratigraphy worldwide.展开更多
Strong low-frequency energy beneath a hydrocarbon reservoir is called a seismic low-frequency shadow and can be used as a hydrocarbon indicator （Tarter et al., 1979） bu the physical mechanism of the observed low-fre...Strong low-frequency energy beneath a hydrocarbon reservoir is called a seismic low-frequency shadow and can be used as a hydrocarbon indicator （Tarter et al., 1979） bu the physical mechanism of the observed low-frequency shadow is still unclear. To stud） the mechanism, we performed seismic numerical simulation of geological models with a hydrocarbon-bearing zone using the 2-D diffusive-viscous wave equation which car effectively model the characteristics of velocity dispersion and transform the seismic dat~ centered in a target layer slice within a time window to the time-frequency domain by usinl time-frequency signal analysis and sort the frequency gathers to common frequency cubes. Then, we observe the characteristics of the seismic low-frequency shadow in the common frequency cubes. The numerical simulations reveal that the main mechanism of seismic lowfrequency shadows is attributed to high attenuation of the medium to high seismic frequency components caused by absorption in the hydrocarbon-filled reservoir. Results from a practical example of seismic low-frequency shadows show that it is possible to identify the reservoir by the low-frequency shadow with high S/N seismic data.展开更多
We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficienc...We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficiency and accuracy between this method and other seismic simulation methods such as finite difference and pseudospec-tral method. Numerical experiments demonstrate that the algorithm constituted by convolutional Forsyte polyno-mial differentiator has high efficiency and accuracy and needs less computational resources, so it is a numerical modeling method with much potential.展开更多
When simulating seismic wave propagation in free space, it is essential to introduce absorbing boundary conditions to eliminate reflections from artificially trtmcated boundaries. In this paper, a damping factor refer...When simulating seismic wave propagation in free space, it is essential to introduce absorbing boundary conditions to eliminate reflections from artificially trtmcated boundaries. In this paper, a damping factor referred to as the Gaussian dmping factor is proposed. The Gaussian damping factor is based on the idea of perfectly matched layers （PMLs）. This work presents a detailed analysis of the theoretical foundations and advantages of the Gaussian damping factor. Additionally, numerical experiments for the simulation of seismic waves are presented based on two numerical models： a homogeneous model and a multi-layer model. The results show that the proposed factor works better. The Gaussian damping factor achieves a higher Signal-to-Noise Ratio （SNR） than previously used factors when using same number of PMLs, and requires less PMLs than other methods to achieve an identical SNR.展开更多
Based on the damage constitutive model for concrete, the Weibull distribution function was used to characterize the random distribution of the mechanical properties of materials by finely subdividing concrete slab ele...Based on the damage constitutive model for concrete, the Weibull distribution function was used to characterize the random distribution of the mechanical properties of materials by finely subdividing concrete slab elements, and a concrete random mesoscopic damage model was established. The seismic response of a 100-m high concrete face rockfill dam(CFRD), subjected to ground motion with different intensities, was simulated with the three-dimensional finite element method(FEM), with emphasis on exploration of damage and the cracking process of concrete slabs during earthquakes as well as analysis of dynamic damage and cracking characteristics during strong earthquakes. The calculated results show that the number of damaged and cracking elements on concrete slabs grows with the duration of earthquakes. With increasing earthquake intensity, the damaged zone and cracking zone on concrete slabs grow wider. During a 7.0-magnitude earthquake, the stress level of concrete slabs is low for the CFRD, and there is almost no damage or slight damage to the slabs. While during a 9.0-magnitude strong earthquake, the percentages of damaged elements and macrocracking elements continuously ascend with the duration of the earthquake, peaking at approximately 26% and 5% at the end of the earthquake, respectively. The concrete random mesoscopic damage model can depict the entire process of sprouting, growing, connecting, and expanding of cracks on a concrete slab during earthquakes.展开更多
To study the bubble plume's seismic response characteristics,the model of a plume water body has been built in this article using the bubble-contained medium acoustic velocity model and the stochastic medium theor...To study the bubble plume's seismic response characteristics,the model of a plume water body has been built in this article using the bubble-contained medium acoustic velocity model and the stochastic medium theory based on an analysis of both the acoustic characteristics of a bubble-contained water body and the actual features of a plume.The finite difference method is used for forward modelling,and the single-shot seismic record exhibits the characteristics of a scattered wave field generated by a plume.A meaningful conclusion is obtained by extracting seismic attributes from the pre-stack shot gather record of a plume.The values of the amplitude-related seismic attributes increase greatly as the bubble content goes up,and changes in bubble radius will not cause seismic attributes to change,which is primarily observed because the bubble content has a strong impact on the plume's acoustic velocity,while the bubble radius has a weak impact on the acoustic velocity.The above conclusion provides a theoretical reference for identifying hydrate plumes using seismic methods and contributes to further study on hydrate decomposition and migration,as well as on distribution of the methane bubble in seawater.展开更多
In the present study seismic wave propagation in heterogeneous media is numerically simulated by using the pseudospectral method with the staggered grid RFFT differentiation in order to clarify the cause for the compl...In the present study seismic wave propagation in heterogeneous media is numerically simulated by using the pseudospectral method with the staggered grid RFFT differentiation in order to clarify the cause for the complicated distribution characteristics of strong ground motion in regions with basin structure. The results show that the maximum amplitudes of simulated ground acceleration waveforms are closely related to the basin structure. Interference of seismic waves in the basin strongly affects the distribution of maximum seismic waveforms, which may result in peak disasters during earthquakes. Peak disasters might be away from basin boundaries or earthquake faults. Seismic energy transmitted into the basin from the bedrock can hardly penetrate the bottom of the basin and then travel back into the bedrock region. The seismic energy is absorbed by basin media, and transferred into the kinematical energy of seismic waves with great amplitude in the basin. Seismic waves between basins may result in serious damage to buildings over the basin. This is significant for aseismatic research. Geological surveys in and around urban areas would benefit aseismatic research and mitigation of seismic disasters of a city. Such geological surveys should involve seismic velocity structure in the media above the bedrock besides such subjects as active faults and geological structure.展开更多
Hybrid simulation can be a cost effective approach for dynamic testing of structural components at full scale while capturing the system level response through interactions with a numerical model.The dynamic response ...Hybrid simulation can be a cost effective approach for dynamic testing of structural components at full scale while capturing the system level response through interactions with a numerical model.The dynamic response of a seismically isolated structure depends on the combined characteristics of the ground motion,bearings,and superstructure.Therefore,dynamic full-scale system level tests of isolated structures under realistic dynamic loading conditions are desirable towards a holistic validation of this earthquake protection strategy.Moreover,bearing properties and their ultimate behavior have been shown to be highly dependent on rate-of-loading and scale size effects,especially under extreme loading conditions.Few laboratory facilities can test full-scale seismic isolation bearings under prescribed displacement and/or loading protocols.The adaptation of a full-scale bearing test machine for the implementation of real-time hybrid simulation is presented here with a focus on the challenges encountered in attaining reliable simulation results for large scale dynamic tests.These advanced real-time hybrid simulations of large and complex hybrid models with several thousands of degrees of freedom are some of the first to use high performance parallel computing to rapidly execute the numerical analyses.Challenges in the experimental setup included measured forces contaminated by delay and other systematic control errors in applying desired displacements.Friction and inertial forces generated by the large-scale loading apparatus can affect the accuracy of measured force feedbacks.Reliable results from real-time hybrid simulation requires implementation of compensation algorithms and correction of these various sources of errors.Overall,this research program confirms that real-time hybrid simulation is a viable testing method to experimentally assess the behavior of full-scale isolators while capturing interactions with the numerical models of the superstructure to evaluate system level and in-structure response.展开更多
Based on a Chinese national high arch dam located in a meizoseismal region, a nonlinear numerical analysis model of the damage and failure process of a dam-foundation system is established by employing a 3-D deformabl...Based on a Chinese national high arch dam located in a meizoseismal region, a nonlinear numerical analysis model of the damage and failure process of a dam-foundation system is established by employing a 3-D deformable distinct element code(3DEC) and its re-development functions. The proposed analysis model considers the dam-foundation-reservoir coupling effect, infl uence of nonlinear contact in the opening and closing of the dam seam surface and abutment rock joints during strong earthquakes, and radiation damping of far fi eld energy dissipation according to the actual workability state of an arch dam. A safety assessment method and safety evaluation criteria is developed to better understand the arch dam system disaster process from local damage to ultimate failure. The dynamic characteristics, disaster mechanism, limit bearing capacity and the entire failure process of a high arch dam under a strong earthquake are then analyzed. Further, the seismic safety of the arch dam is evaluated according to the proposed evaluation criteria and safety assessment method. As a result, some useful conclusions are obtained for some aspects of the disaster mechanism and failure process of an arch dam. The analysis method and conclusions may be useful in engineering practice.展开更多
US Vice President Al Gore's vision of Digital Earth applies us with prospects for brand-new ways of solving problems the earth is facing such as seismic disaster. ms paper first briefly introduces the concept of ...US Vice President Al Gore's vision of Digital Earth applies us with prospects for brand-new ways of solving problems the earth is facing such as seismic disaster. ms paper first briefly introduces the concept of Digital Earth. Then in the context of Digital Earth. the Origin, concept and application of Virtual Reality technology are reviewed. After that we present in detail our preliminary case study--CVR-USD (Computer Virtual Reality for Urban Seismic Disaster Simulation) System which aims to simulate and manage seismic disaster through integrating RS, GIS and VR technologies. For this system, we've built USD subsystem, developed SMVR software to implement CVR. and also developed a Spatial Dare Analysis Package to handle spatial data related to earthquake disaster.展开更多
Variable curvature friction pendulum bearings(VCFPB)effectively reduce the dynamic response of storage tanks induced by earthquakes.Shaking table testing is used to assess the seismic performance of VCFPB isolated sto...Variable curvature friction pendulum bearings(VCFPB)effectively reduce the dynamic response of storage tanks induced by earthquakes.Shaking table testing is used to assess the seismic performance of VCFPB isolated storage tanks.However,the vertical pressure and friction coefficient of the scaled VCFPB in the shaking table tests cannot match the equivalent values of these parameters in the prototype.To avoid this drawback,a real-time hybrid simulation(RTHS)test was developed.Using RTHS testing,a 1/8 scaled tank isolated by VCFPB was tested.The experimental results showed that the displacement dynamic magnification factor of VCFPB,peak reduction factors of the acceleration,shear force,and overturning moment at bottom of the tank,were negative exponential functions of the ratio of peak ground acceleration(PGA)and friction coefficient.The peak reduction factors of displacement,acceleration,force and overturning moment,which were obtained from the experimental results,are compared with those calculated by the Housner model.It can be concluded that the Housner model is applicable in estimation of the acceleration,shear force,and overturning moment of liquid storage tank,but not for the sliding displacement of VCFPBs.展开更多
The energy caused by the dynamic impact in mining engineering forth release and spread by the way of seismic waves, monitoring is an effective way for forecasting mine dynamical disasters, such as rockburst and coal a...The energy caused by the dynamic impact in mining engineering forth release and spread by the way of seismic waves, monitoring is an effective way for forecasting mine dynamical disasters, such as rockburst and coal and gas outburst. Three-dimensional dynamic model was built to simulate the propagating progress of seismic waves in the elastoplastic tunnel rock and analyzed the propagating law of perturbation acceleration around tunnel, based on the finite element dynamic analysis software ANSYS/L S-DYNA. The simulation results indicate that： （1） The propagation attenuation of seismic wave is a negative index relationship; （2） The acceleration amplitude of seismic wave decays rapidly in near-field and decays slowly in far-field; （3） When the perturbation is generated in the dead ahead of tunnel, the acceleration of seismic wave become smaller and smaller away from the roadway-rib;（4） The elastic and plastic stress state of tunnel rock is also an important factor for propagation process of wave, the energy of seismic wave is mainly consumed for geometric spreading and plastic deformation in propagation in the elastoplastic medium model.展开更多
In order to carry out numerical simulation using geologic structural data obtained from Landmark(seismic interpretation system), underground geological structures are abstracted into mechanical models which can reflec...In order to carry out numerical simulation using geologic structural data obtained from Landmark(seismic interpretation system), underground geological structures are abstracted into mechanical models which can reflect actual situations and facilitate their computation and analyses.Given the importance of model building, further processing methods about traditional seismic interpretation results from Landmark should be studied and the processed result can then be directly used in numerical simulation computations.Through this data conversion procedure, Landmark and FLAC(the international general stress software) are seamlessly connected.Thus, the format conversion between the two systems and the pre-and post-processing in simulation computation is realized.A practical application indicates that this method has many advantages such as simple operation, high accuracy of the element subdivision and high speed, which may definitely satisfy the actual needs of floor grid cutting.展开更多
There is little low-and-high frequency information on seismic data in seismic exploration,resulting in narrower bandwidth and lower seismic resolution.It considerably restricts the prediction accuracy of thin reservoi...There is little low-and-high frequency information on seismic data in seismic exploration,resulting in narrower bandwidth and lower seismic resolution.It considerably restricts the prediction accuracy of thin reservoirs and thin interbeds.This study proposes a novel method to constrain improving seismic resolution in the time and frequency domain.The expected wavelet spectrum is used in the frequency domain to broaden the seismic spectrum range and increase the octave.In the time domain,the Frobenius vector regularization of the Hessian matrix is used to constrain the horizontal continuity of the seismic data.It eff ectively protects the signal-to-noise ratio of seismic data while the longitudinal seismic resolution is improved.This method is applied to actual post-stack seismic data and pre-stack gathers dividedly.Without abolishing the phase characteristics of the original seismic data,the time resolution is signifi cantly improved,and the structural features are clearer.Compared with the traditional spectral simulation and deconvolution methods,the frequency distribution is more reasonable,and seismic data has higher resolution.展开更多
基金Scientific Research Fund of Institute of Engineering Mechanics,China Earthquake Administration under Grant No.2023C01National Natural Science Foundation of China under Grant No.52478570Distinguished Young Scholars Program of the Natural Science Foundation of Heilongjiang Province,China under Grant No.JQ2024E002。
文摘The 2025 M_(w)7.7 Myanmar earthquake highlighted the challenge of near-fault seismic intensity field reconstruction due to sparse seismic networks.To address this limitation,a framework was proposed integrating seismic wave simulation with a data-constrained finite-fault rupture model.The constraint is implemented by identifying the optimal ground motion models(GMMs)through a scoring system that selects the best-fit GMMs to mid-and far-field China Earthquake Networks Center(CENC)seismic network data;and applying the optimal GMMs to refine the rupture model parameters for near-fault intensity field simulation.The simulated near-fault seismic intensity field reproduces seismic intensities collected from Myanmar’s sparse seismic network and concentrated in≥Ⅷintensity zones within 50 km of the projected fault plane;and identifies abnormal intensity regions exhibiting≥Ⅹintensity along the Meiktila-Naypyidaw corridor and near Shwebo that are attributed to soft soil amplification effects and near-fault directivity.This framework can also be applied to post-earthquake assessments in other similar regions.
文摘We develop and implement a Stochastic Discrete Event Simulation(SDES)algorithm to model the housing re-covery trajectory after an extreme event.The algorithm models discrete events and their underlying uncertainties in each construction phase.Specifically,the algorithm is developed for the Government Assisted Owner Driven(GAOD)reconstruction system to simulate long-term recovery trajectory.SDES,as a flexible modeling approach,can simulate any housing recovery scenario that follows phased reconstruction.The 2015 M 7.8 Gorkha earthquake sequence in Nepal is considered the extreme event,with 796,245 buildings requiring reconstruction.We present some recovery trajectories from severely hit,crisis hit,and earthquake hit parishes,comparing them with the actual reconstruction progress.We also assess quality and improvement of reconstructed buildings using seismic fragility functions,compared to pre-earthquake constructions.Housing recovery uncertainties are dissected in relation to reconstruction pace.We conclude that the vast majority of the reconstructed buildings followed the Build Back Better(BBB)approach and missed the opportunity to pursue the Build Back Resilient(BBR)approach due to multifaceted challenges ranging from unclear policies to economic constraints.We critically assess the GAOD vs Owner Driven(OD)recovery framework and conclude that insurance-supported and technically assisted OD approach could be the most suitable model for post extreme event housing recovery.
基金financially supported by the National Key Research and Development Program (No. 2016YFC0801402)the National Natural Science Foundation of China (No. 51374236)Chongqing Research Program of Basic Research and Frontier Technology of China (No. cstc2015jcyj BX0076)
文摘Gas extraction practice has been proven for the clear majority of coal mines in China to be unfavorable using drill holes in the coal seam. Rather, mining-induced fractures in the goaf should be utilized for gas extraction. To study gas migration in mining-induced fractures, one mining face of 10 th Mine in Pingdingshan Coalmine Group in Henan, China, has been selected as the case study for this work. By establishing the mathematical model of gas migration under the influence of coal seam mining, discrete element software UDEC and Multiphysics software COMSOL are employed to model gas migration in mining-induced fractures above the goaf. The results show that as the working face advances, the goaf overburden gradually forms a mining-induced fracture network in the shape of a trapezoid, the size of which increases with the distance of coal face advance. Compared with gas migration in the overburden matrix, the gas flow in the fracture network due to mining is far greater. The largest mining-induced fracture is located at the upper end of the trapezoidal zone, which results in the largest gas flux in the network. When drilling for gas extraction in a mining-induced fracture field, the gas concentration is reduced in the whole region during the process of gas drainage, and the rate of gas concentration drops faster in the fractured zone. It is shown that with gas drainage, the gas flow velocity in the mininginduced fracture network is faster.
基金supported by the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(No.2021QNLM020001)the Major Scientific and Technological Projects of Shandong Energy Group(No.SNKJ2022A06-R23)+1 种基金the Funds of Creative Research Groups of China(No.41821002)the Major Scientific and Technological Projects of CNPC(No.ZD2019-183-003).
文摘To reduce the spatial simulation error generated by the finite difference method,previous researchers compute the optimal finite-difference weights always by minimizing the error of spatial dispersion relation.However,we prove that the spatial simulation error of the finite difference method is associated with the dot product of the spatial dispersion relation of the finite-difference weights and the spectrum of the seismic wavefield.Based on the dot product relation,we construct a L_(2) norm cost function to minimize spatial simulation error.For solving this optimization problem,the seismic wavefield infor-mation in wavenumber region is necessary.Nevertheless,the seismic wavefield is generally obtained by costly forward modeling techniques.To reduce the computational cost,we substitute the spectrum of the seismic wavelet for the spectrum of the seismic wavefield,as the seismic wavelet plays a key role in determining the seismic wavefield.In solving the optimization problem,we design an exhaustive search method to obtain the solution of the L_(2) norm optimization problem.After solving the optimization problem,we are able to achieve the finite-difference weights that minimize spatial simulation error.In theoretical error analyses,the finite-difference weights from the proposed method can output more accurate simulation results compared to those from previous optimization algorithms.Furthermore,we validate our method through numerical tests with synthetic models,which encompass homogenous/inhomogeneous media as well as isotropic and anisotropic media.
基金National Key R&D Program of China under Grant No.2022YFC3003601。
文摘The March 28,2025 Myanmar earthquake generated ground shaking that was perceptible throughout Myanmar and adjacent regions.This study simulated three-component ground motions across the affected region using an improved stochastic finite-fault method to systematically assess seismic impacts.Observed near-field recordings at MM.NGU station was used to determine the reliability of the theoretically derived stress drop as input for simulation.Far-field recordings constrained the frequency-dependent S-wave quality factors(Q(f)=283.305f^(0.588))for anelastic attenuation modeling.Comparisons of peak accelerations between simulation and empirical ground-motion models showed good agreement at moderate-to-large distances.However,lower near-fault simulations indicate a weaker-than-average source effect.Analysis of simulated instrumental seismic intensity revealed key patterns.Maximum intensity(Ⅹ)occurred in isolated patches within the ruptured fault projection,correlating with shallow high-slip areas.TheⅨ-intensity zone formed a north-south elongated band centered on fault projection.Significant asymmetry inⅧ-intensity distribution perpendicular to the fault strike was observed,with a wider western extension attributed to lower shear-wave velocities west of the fault.Supershear rupture behavior enhanced ground motions,expanding intensity ranges by~20%compared to sub-shear rupture.This study reveals the integrated effects of fault geometry,slip spatial distribution,rupture velocity,and site condition in governing ground motion patterns.
基金financially supported by the State Key Research Development Program of China(Grant No.2016YFC0600701)the National Natural Science Foundation of China(Grant No.51674170)
文摘In this study, the spatial distributions of stress and fracture fields for three typical underground coal mining layouts, Le, non-pillar mining （NM）, top-coal caving mining （TCM） and protective coal-seam mining （PCM）, are modeled using discrete element software UDEC, The numerical results show that different mining layouts can lead to different mining-induced stress fields, resulting in diverse fracture fields, For the PCM, the mining influenced area in front of the mining faces is the largest, and the stress concentration factor in front of the mining faces is the lowest, The spatial shapes of the mining-induced fracture fields under NM, TCM and PCM differ, and they are characterized by trapezoidal, triangular and tower shapes, respectively, The fractal dimensions of mining-induced fractures of the three mining layouts decrease in the order of PCM, TCM and NM, It is also shown that the PCM can result in a better gas control effect in coal mines with high outburst potential, The numerical results are expected to provide a basis for understanding of mining-induced gas seepage fields and provide a reference for high- efficiency coal mining,
文摘Seismic quantitative reservoir simulations and characterizations have played a vital role in exploring stratigraphic traps,such as lateaggradational lowstands prograding wedge systems(LPWS)within lowstands systems tracts(LST).However,seismic data acquisition operations are always dominated by exceptional seismic coherent noise events,e.g.,multiples,which reduce the signal strengths of the sourcegenerated incident seismic waves within vertically and laterally heterogeneous earth systems.Hence,these noise events create hurdles in predicting paleo-depositional impedance(PDI),paleo-thickness(PTS),paleo-dense fractured networks,erosional and depositional zones,faultcontrolled migrations,and types of seismic reflection configurations(SRFC),which are key elements in developing stratigraphic pinch-out traps.This research utilizes the state-of-the-art technologies of spectral wavelet-based instantaneous time-frequency analysis and seismic waveform frequency-controlled porosity-constrained static reservoir simulation(FDPVS)tools to quantify the LPWS inside the Onshore Basin,Pakistan.The use of conventional amplitude-based seismic attributes,such as the average energy,remained a better tool for deciphering the overall geological architecture of the LPWS.Conventional FDPVS realizations resolved a PDI of−1.391 gm./c.c.^(*)m/s to−0.97 gm./c.c.^(*)m/s for LPWS with PTS of 12 and 20 m,respectively.A 0.9 km lateral extent of paleo-dense fractured networks(PDFN)with a strong linear regression R^(2)=0.93 was also resolved.Average energy attribute-based instantaneous frequency FDPVS realizations enabled the imaging of parallel-toprograding SRFC with resolved magnitudes of−0.259 gm./c.c.^(*)m/s for PDI,20 m for PTS,and 0.73 km for PDFN with linear regression transforms at R^(2)=0.92,which indicates the deposition of onlap fill facies inside the LPWS during extensive sea-level fall.These realizations have also resolved frequency-controlled fault migrations on 27-Hz spectral waveform-based amplitude plots with 2.174 gm./c.c.^(*)m/s PDI for conduit fault systems and 27-Hz with 0.585 gm./c.c.^(*)m/s PDI for sealing fault systems.All these structural configurations are completely sealed up by transgressive seals of transgressive systems tracts and,hence,developed into pure stratigraphic-based oil and gas plays.This research has strong implications for side-tracking drilling locations and provides an analogue for basins with similar geology and stratigraphy worldwide.
基金supported by the National Hi-tech Research and Development Program of China (863 Program) (Grant No. 2006AA0AA 02 - 2).
文摘Strong low-frequency energy beneath a hydrocarbon reservoir is called a seismic low-frequency shadow and can be used as a hydrocarbon indicator （Tarter et al., 1979） bu the physical mechanism of the observed low-frequency shadow is still unclear. To stud） the mechanism, we performed seismic numerical simulation of geological models with a hydrocarbon-bearing zone using the 2-D diffusive-viscous wave equation which car effectively model the characteristics of velocity dispersion and transform the seismic dat~ centered in a target layer slice within a time window to the time-frequency domain by usinl time-frequency signal analysis and sort the frequency gathers to common frequency cubes. Then, we observe the characteristics of the seismic low-frequency shadow in the common frequency cubes. The numerical simulations reveal that the main mechanism of seismic lowfrequency shadows is attributed to high attenuation of the medium to high seismic frequency components caused by absorption in the hydrocarbon-filled reservoir. Results from a practical example of seismic low-frequency shadows show that it is possible to identify the reservoir by the low-frequency shadow with high S/N seismic data.
基金Open Fund of State Key Laboratory of Geological Processes and Mineral Resources, China University of Geo-sciences (GPMR0750)National Natural Science Foundation of China (40437018)
文摘We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficiency and accuracy between this method and other seismic simulation methods such as finite difference and pseudospec-tral method. Numerical experiments demonstrate that the algorithm constituted by convolutional Forsyte polyno-mial differentiator has high efficiency and accuracy and needs less computational resources, so it is a numerical modeling method with much potential.
基金supported by the National Natural Science Foundation of China(No. 61072118)
文摘When simulating seismic wave propagation in free space, it is essential to introduce absorbing boundary conditions to eliminate reflections from artificially trtmcated boundaries. In this paper, a damping factor referred to as the Gaussian dmping factor is proposed. The Gaussian damping factor is based on the idea of perfectly matched layers （PMLs）. This work presents a detailed analysis of the theoretical foundations and advantages of the Gaussian damping factor. Additionally, numerical experiments for the simulation of seismic waves are presented based on two numerical models： a homogeneous model and a multi-layer model. The results show that the proposed factor works better. The Gaussian damping factor achieves a higher Signal-to-Noise Ratio （SNR） than previously used factors when using same number of PMLs, and requires less PMLs than other methods to achieve an identical SNR.
基金supported by the Key Laboratory of Failure Mechanism and Safety Control Techniques of Earth-rock Dams of the Ministry of Water Resources(Grant No.YK914019)the CRSRI Open Research Program(Grant No.CKWV2016376/KY)the National Natural Science Foundation of China(Grant No.51009055)
文摘Based on the damage constitutive model for concrete, the Weibull distribution function was used to characterize the random distribution of the mechanical properties of materials by finely subdividing concrete slab elements, and a concrete random mesoscopic damage model was established. The seismic response of a 100-m high concrete face rockfill dam(CFRD), subjected to ground motion with different intensities, was simulated with the three-dimensional finite element method(FEM), with emphasis on exploration of damage and the cracking process of concrete slabs during earthquakes as well as analysis of dynamic damage and cracking characteristics during strong earthquakes. The calculated results show that the number of damaged and cracking elements on concrete slabs grows with the duration of earthquakes. With increasing earthquake intensity, the damaged zone and cracking zone on concrete slabs grow wider. During a 7.0-magnitude earthquake, the stress level of concrete slabs is low for the CFRD, and there is almost no damage or slight damage to the slabs. While during a 9.0-magnitude strong earthquake, the percentages of damaged elements and macrocracking elements continuously ascend with the duration of the earthquake, peaking at approximately 26% and 5% at the end of the earthquake, respectively. The concrete random mesoscopic damage model can depict the entire process of sprouting, growing, connecting, and expanding of cracks on a concrete slab during earthquakes.
基金financially supported by the National Natural Science Foundation of China(No.41306050)the Science and Technology Project of Guangdong Province of China(No.2014A010103030)+1 种基金the Natural Science Foundation of Guangdong Province under contract(No.2015A030313617)the National Marine Important Charity Special Foundation of China(No.201305019)
文摘To study the bubble plume's seismic response characteristics,the model of a plume water body has been built in this article using the bubble-contained medium acoustic velocity model and the stochastic medium theory based on an analysis of both the acoustic characteristics of a bubble-contained water body and the actual features of a plume.The finite difference method is used for forward modelling,and the single-shot seismic record exhibits the characteristics of a scattered wave field generated by a plume.A meaningful conclusion is obtained by extracting seismic attributes from the pre-stack shot gather record of a plume.The values of the amplitude-related seismic attributes increase greatly as the bubble content goes up,and changes in bubble radius will not cause seismic attributes to change,which is primarily observed because the bubble content has a strong impact on the plume's acoustic velocity,while the bubble radius has a weak impact on the acoustic velocity.The above conclusion provides a theoretical reference for identifying hydrate plumes using seismic methods and contributes to further study on hydrate decomposition and migration,as well as on distribution of the methane bubble in seawater.
基金This study is financiallysupported by the 973 Research Project of the Ministry ofScience and Technology (No．2003CB716505) aresearch proiect of Chinese Continental Science DrillingCenter(No．2002207)．
文摘In the present study seismic wave propagation in heterogeneous media is numerically simulated by using the pseudospectral method with the staggered grid RFFT differentiation in order to clarify the cause for the complicated distribution characteristics of strong ground motion in regions with basin structure. The results show that the maximum amplitudes of simulated ground acceleration waveforms are closely related to the basin structure. Interference of seismic waves in the basin strongly affects the distribution of maximum seismic waveforms, which may result in peak disasters during earthquakes. Peak disasters might be away from basin boundaries or earthquake faults. Seismic energy transmitted into the basin from the bedrock can hardly penetrate the bottom of the basin and then travel back into the bedrock region. The seismic energy is absorbed by basin media, and transferred into the kinematical energy of seismic waves with great amplitude in the basin. Seismic waves between basins may result in serious damage to buildings over the basin. This is significant for aseismatic research. Geological surveys in and around urban areas would benefit aseismatic research and mitigation of seismic disasters of a city. Such geological surveys should involve seismic velocity structure in the media above the bedrock besides such subjects as active faults and geological structure.
文摘Hybrid simulation can be a cost effective approach for dynamic testing of structural components at full scale while capturing the system level response through interactions with a numerical model.The dynamic response of a seismically isolated structure depends on the combined characteristics of the ground motion,bearings,and superstructure.Therefore,dynamic full-scale system level tests of isolated structures under realistic dynamic loading conditions are desirable towards a holistic validation of this earthquake protection strategy.Moreover,bearing properties and their ultimate behavior have been shown to be highly dependent on rate-of-loading and scale size effects,especially under extreme loading conditions.Few laboratory facilities can test full-scale seismic isolation bearings under prescribed displacement and/or loading protocols.The adaptation of a full-scale bearing test machine for the implementation of real-time hybrid simulation is presented here with a focus on the challenges encountered in attaining reliable simulation results for large scale dynamic tests.These advanced real-time hybrid simulations of large and complex hybrid models with several thousands of degrees of freedom are some of the first to use high performance parallel computing to rapidly execute the numerical analyses.Challenges in the experimental setup included measured forces contaminated by delay and other systematic control errors in applying desired displacements.Friction and inertial forces generated by the large-scale loading apparatus can affect the accuracy of measured force feedbacks.Reliable results from real-time hybrid simulation requires implementation of compensation algorithms and correction of these various sources of errors.Overall,this research program confirms that real-time hybrid simulation is a viable testing method to experimentally assess the behavior of full-scale isolators while capturing interactions with the numerical models of the superstructure to evaluate system level and in-structure response.
基金National Natural Science Foundation of China under Grant No.90510017
文摘Based on a Chinese national high arch dam located in a meizoseismal region, a nonlinear numerical analysis model of the damage and failure process of a dam-foundation system is established by employing a 3-D deformable distinct element code(3DEC) and its re-development functions. The proposed analysis model considers the dam-foundation-reservoir coupling effect, infl uence of nonlinear contact in the opening and closing of the dam seam surface and abutment rock joints during strong earthquakes, and radiation damping of far fi eld energy dissipation according to the actual workability state of an arch dam. A safety assessment method and safety evaluation criteria is developed to better understand the arch dam system disaster process from local damage to ultimate failure. The dynamic characteristics, disaster mechanism, limit bearing capacity and the entire failure process of a high arch dam under a strong earthquake are then analyzed. Further, the seismic safety of the arch dam is evaluated according to the proposed evaluation criteria and safety assessment method. As a result, some useful conclusions are obtained for some aspects of the disaster mechanism and failure process of an arch dam. The analysis method and conclusions may be useful in engineering practice.
基金National Natural Science Foundation of China (NSFC)!project No. 69896250 Basic Research Project of CAS !project No. KJ951-B1
文摘US Vice President Al Gore's vision of Digital Earth applies us with prospects for brand-new ways of solving problems the earth is facing such as seismic disaster. ms paper first briefly introduces the concept of Digital Earth. Then in the context of Digital Earth. the Origin, concept and application of Virtual Reality technology are reviewed. After that we present in detail our preliminary case study--CVR-USD (Computer Virtual Reality for Urban Seismic Disaster Simulation) System which aims to simulate and manage seismic disaster through integrating RS, GIS and VR technologies. For this system, we've built USD subsystem, developed SMVR software to implement CVR. and also developed a Spatial Dare Analysis Package to handle spatial data related to earthquake disaster.
基金Scientific Research Fund of Institute of Engineering Mechanics,China Earthquake Administration under Grant No.2018D03the National Natural Science Foundation of China under Grant Nos.51608016 and 51421005。
文摘Variable curvature friction pendulum bearings(VCFPB)effectively reduce the dynamic response of storage tanks induced by earthquakes.Shaking table testing is used to assess the seismic performance of VCFPB isolated storage tanks.However,the vertical pressure and friction coefficient of the scaled VCFPB in the shaking table tests cannot match the equivalent values of these parameters in the prototype.To avoid this drawback,a real-time hybrid simulation(RTHS)test was developed.Using RTHS testing,a 1/8 scaled tank isolated by VCFPB was tested.The experimental results showed that the displacement dynamic magnification factor of VCFPB,peak reduction factors of the acceleration,shear force,and overturning moment at bottom of the tank,were negative exponential functions of the ratio of peak ground acceleration(PGA)and friction coefficient.The peak reduction factors of displacement,acceleration,force and overturning moment,which were obtained from the experimental results,are compared with those calculated by the Housner model.It can be concluded that the Housner model is applicable in estimation of the acceleration,shear force,and overturning moment of liquid storage tank,but not for the sliding displacement of VCFPBs.
文摘The energy caused by the dynamic impact in mining engineering forth release and spread by the way of seismic waves, monitoring is an effective way for forecasting mine dynamical disasters, such as rockburst and coal and gas outburst. Three-dimensional dynamic model was built to simulate the propagating progress of seismic waves in the elastoplastic tunnel rock and analyzed the propagating law of perturbation acceleration around tunnel, based on the finite element dynamic analysis software ANSYS/L S-DYNA. The simulation results indicate that： （1） The propagation attenuation of seismic wave is a negative index relationship; （2） The acceleration amplitude of seismic wave decays rapidly in near-field and decays slowly in far-field; （3） When the perturbation is generated in the dead ahead of tunnel, the acceleration of seismic wave become smaller and smaller away from the roadway-rib;（4） The elastic and plastic stress state of tunnel rock is also an important factor for propagation process of wave, the energy of seismic wave is mainly consumed for geometric spreading and plastic deformation in propagation in the elastoplastic medium model.
基金Projects 50221402, 50490271 and 50025413 supported by the National Natural Science Foundation of Chinathe National Basic Research Program of China (2009CB219603, 2009 CB724601, 2006CB202209 and 2005CB221500)+1 种基金the Key Project of the Ministry of Education (306002)the Program for Changjiang Scholars and Innovative Research Teams in Universities of MOE (IRT0408)
文摘In order to carry out numerical simulation using geologic structural data obtained from Landmark(seismic interpretation system), underground geological structures are abstracted into mechanical models which can reflect actual situations and facilitate their computation and analyses.Given the importance of model building, further processing methods about traditional seismic interpretation results from Landmark should be studied and the processed result can then be directly used in numerical simulation computations.Through this data conversion procedure, Landmark and FLAC(the international general stress software) are seamlessly connected.Thus, the format conversion between the two systems and the pre-and post-processing in simulation computation is realized.A practical application indicates that this method has many advantages such as simple operation, high accuracy of the element subdivision and high speed, which may definitely satisfy the actual needs of floor grid cutting.
基金supported by the PetroChina Prospective,Basic,and Strategic Technology Research Project(No.2021DJ0606).
文摘There is little low-and-high frequency information on seismic data in seismic exploration,resulting in narrower bandwidth and lower seismic resolution.It considerably restricts the prediction accuracy of thin reservoirs and thin interbeds.This study proposes a novel method to constrain improving seismic resolution in the time and frequency domain.The expected wavelet spectrum is used in the frequency domain to broaden the seismic spectrum range and increase the octave.In the time domain,the Frobenius vector regularization of the Hessian matrix is used to constrain the horizontal continuity of the seismic data.It eff ectively protects the signal-to-noise ratio of seismic data while the longitudinal seismic resolution is improved.This method is applied to actual post-stack seismic data and pre-stack gathers dividedly.Without abolishing the phase characteristics of the original seismic data,the time resolution is signifi cantly improved,and the structural features are clearer.Compared with the traditional spectral simulation and deconvolution methods,the frequency distribution is more reasonable,and seismic data has higher resolution.
