With the continuous improvement of the accuracy of geodetic deformation data,the inversion of seismic source parameters puts forward a higher demand for nonlinear inversion algorithms.In this research,an improved Spar...With the continuous improvement of the accuracy of geodetic deformation data,the inversion of seismic source parameters puts forward a higher demand for nonlinear inversion algorithms.In this research,an improved Sparrow Search Algorithm(SSA)is proposed for the seismic source parameter inversion problem.By replacing the original population generation in the improved algorithm with Latin hypercubic sampling,the Sparrow Search Algorithm reduces the repetition of samples in the population initialization.Subsequently,the algorithm introduces adaptive weights in the discoverer generation phase of the sparrow algorithm and combines the Levy flight strategy to make the algorithm more comprehensive and improve the search accuracy during the whole iteration process.Therefore,the improved Latin hypercube-based sparrow search algorithm(ILHSSA)has better advantages in terms of iterative convergence speed and stability.In order to verify the performance of ILHSSA,the basic genetic algorithm(GA)and sparrow search algorithm(SSA)are examined and compared with ILHSSA by simulated earthquakes of two different earthquake types.The simulation experiments show that the improved algorithm ILHSSA outperforms SSA in accuracy and stability.Compared with the GA algorithm,ILHSSA can achieve the same inversion accuracy as GA,and it even surpasses GA in inversion speed and the inversion results of some parameters,demonstrating better stability.Finally,the improved algorithm is used for the 2017 Bodrum-Cos earthquake and the 2016 Amatrice earthquake in Italy.The inversion results all reflect the practicality and reliability of the improved algorithm.展开更多
Presented in this study is a novel method for estimating the depth of single underwater source in shallow water,utilizing vector sensors.The approach leverages the depth distribution of the broadband Stokes parameters...Presented in this study is a novel method for estimating the depth of single underwater source in shallow water,utilizing vector sensors.The approach leverages the depth distribution of the broadband Stokes parameters to estimate source depth accurately.Unlike traditional matched field processing(MFP)and matched mode processing(MMP),the proposed approach can estimate source depth directly from the data received by sensors without requiring complete environmental information.Firstly,the broadband Stokes parameters(BSP)are established using the normal mode theory.Then the nonstationary phase approximation is used to simplify the theoretical derivation,which is necessary when dealing with broadband integrals.Additionally,range terms of the BSP are eliminated by normalization.By analyzing the depth distribution of the normalized broadband Stokes parameters(NBSP),it is found that the NBSP exhibit extreme values at the source depth,which can be used for source depth estimation.So the proposed depth estimation method is based on searching the peaks of the NBSP.Simulations show that this method is effective in relatively simple shallow water environments.Finally,the effect of source range,frequency bandwidth,sound speed profile(SSP),water depth,and signal-to-noise ratio(SNR)are studied.The findings indicate that the proposed method can accurately estimate the source depth when the SNR is greater than-5 d B and does not need to consider model mismatch issues.Additionally,variations in environmental parameters have minimal impact on estimation accuracy.Compared to MFP,the proposed method requires a higher SNR,but demonstrates superior robustness against fluctuations in environmental parameters.展开更多
An earthquake of M S=6.9 occurred in Gonghe County, Qinghai Province, China on April 26, 1990.This earthquake was followed by three larger aftershocks of M S=5.5 on May 7, 1990, M S=6.0 on Jan.3, 199...An earthquake of M S=6.9 occurred in Gonghe County, Qinghai Province, China on April 26, 1990.This earthquake was followed by three larger aftershocks of M S=5.5 on May 7, 1990, M S=6.0 on Jan.3, 1994, and M S=5.7 on Feb.16, 1994, consecutively. The moment tensors of these earthquakes as function of time were obtained by the technique of moment tensor inversion in frequency domain . The results inverted indicate that these earthquakes had a very similar focal mechanism of predominantly reverse faulting on a plane striking NWW, dipping to SSW.The scalar seismic moments of these earthquakes are M 0=9.4×10 18 Nm for the M S=6.9 event, 8.0×10 16 Nm for the M S=5.5 event, 4.9×10 17 Nm for the M S =6.0 event and 2.9×10 17 Nm for the M S=5.7 event, respectively. The results inverted also show that the source processes of these events were significantly different. The main shock had a very complex process, consisting of two distinct sub events with comparable sizes. The first sub event occurred in the first 12s, having a seismic moment of 4.7×10 18 Nm, and the second one continued from 31s to 41s, having a seismic moment of 2.5×10 18 Nm. In addition, a much smaller sub event, having a seismic moment of about 2.1×10 18 Nm, may exist in the interval of 12 s and 31 s, In contrast, the source processes of the three aftershocks are quite simple. The source time function of each of aftershocks is a single impulse, suggestting that each of aftershocks consists of a mainly uninterrupted rupture. The rise times and total rupture durations are 4 s and 11 s for the M S=5.5 event, 6 s and 16 s for the M S= 6.0 event and 6 s and 13 s for the M S=5.7 event, respectively.展开更多
The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault par...The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault parameters inversion has nonlinear characteristics, and the gradient-based optimizer(GBO) has the characteristics of fast convergence speed and falling into local optimum hardly. This paper applies GBO algorithm to simulated earthquakes and real LuShan earthquakes in the nonlinear inversion of the Okada model to obtain the source parameters. The simulated earthquake experiment results show that the algorithm is stable, and the seismic source parameters obtained by GBO are slightly closer to the true value than the multi peak particle swarm optimization(MPSO). In the 2013 LuShan earthquake experiment, the root mean square error between the deformation after forwarding of fault parameters obtained by the introduced GBO algorithm and the surface observation deformation was 3.703 mm, slightly better than 3.708 mm calculated by the MPSO. Moreover, the inversion result of GBO algorithm is better than MPSO algorithm in stability. The above results show that the introduced GBO algorithm has a certain practical application value in seismic fault source parameters inversion.展开更多
The seismoacoustic analysis method has broad potential applications to source parameter estimation for near-surface explosion events such as industrial explosions and terrorist attacks.In this study,current models wer...The seismoacoustic analysis method has broad potential applications to source parameter estimation for near-surface explosion events such as industrial explosions and terrorist attacks.In this study,current models were improved by modifying the acoustic model and adopting the Bayesian Markov-chain-Monte-Carlo inversion method.The source parameters of near-surface small-yield chemical explosions were analyzed via the improved seismoacoustic analysis model and by the estimation accuracy of seismoacoustic joint inversion.Estimation and analysis results showed that the improved seismoacoustic analysis model considered ground shock coupling and the impact of explosion products ejecting from the surface so that the improved acoustic impulse relation was more consistent with the measured data than the Ford impulse relation.It is suitable for deep-burial,shallow-burial,and near-surface aerial explosions.Furthermore,trade-off relationships were declined through the application of the improved model to source parameter inversion for near-surface small-yield chemical explosions,and source parameter estimation accuracy was improved.展开更多
According to the geological structural features, Beijing and the adjacent areas can be divided into two regions of plain in the east and mountain in the west. Among the stations covered by the telemetered digital seis...According to the geological structural features, Beijing and the adjacent areas can be divided into two regions of plain in the east and mountain in the west. Among the stations covered by the telemetered digital seismic station network of Earthquake Administration of Beijing Municipality, the stations in the plain area are all borehole ones and the stations in the western mountainous region are all located on the surface bedrock. In the paper, 511 wave- form data recorded by the network from Oct. 2001 to Oct. 2004 are used in the researches for the entire Beijing region, the western mountainous region and the eastern plain area, respectively. The Q values are calculated for each area by Atkinson′s method and compared with the existed data. The reliability of the Q values and the reasons for the difference in the Q values are also discussed. Then, the source parameters and site response are inverted by the Moya′s method, in which two models are used. The first model uses the Q values, earthquakes and stations in the sub-areas and the second model uses the Q values, earthquakes and stations in the entire Beijing region. The results indicate that the source parameters and site responses obtained by two models are basically consistent with each other. It also indicates that the source parameters obtained by these methods are not affected by the size of station network.展开更多
A method for simultaneous determination of mixed model parameters,which have different physical dimensions or different responses to data,is presented.Mixed parameter estimation from observed data within a single mode...A method for simultaneous determination of mixed model parameters,which have different physical dimensions or different responses to data,is presented.Mixed parameter estimation from observed data within a single model space shows instabilities and trade-offs of the solutions. We separate the model space into N-subspaces based on their physical properties or computational convenience and solve the N-subspaces systems by damped least-squares and singular-value decomposition. Since the condition number of each subsystem is smaller than that of the single global system,the approach can greatly increase the stability of the inversion. We also introduce different damping factors into the subsystems to reduce the tradeoffs between the different parameters. The damping factors depend on the conditioning of the subsystems and may be adequately chosen in a range from 0.1 % to 10 % of the largest singular value. We illustrate the method with an example of simultaneous determination of source history,source geometry,and hypocentral location from regional seismograms,although it is applicable to any geophysical inversion.展开更多
Fluid identification and anisotropic parameters characterization are crucial for shale reservoir exploration and development.However,the anisotropic reflection coefficient equation,based on the transverse isotropy wit...Fluid identification and anisotropic parameters characterization are crucial for shale reservoir exploration and development.However,the anisotropic reflection coefficient equation,based on the transverse isotropy with a vertical axis of symmetry(VTI)medium assumption,involves numerous parameters to be inverted.This complexity reduces its stability and impacts the accuracy of seismic amplitude variation with offset(AVO)inversion results.In this study,a novel anisotropic equation that includes the fluid term and Thomsen anisotropic parameters is rewritten,which reduces the equation's dimensionality and increases its stability.Additionally,the traditional Markov Chain Monte Carlo(MCMC)inversion algorithm exhibits a high rejection rate for random samples and relies on known parameter distributions such as the Gaussian distribution,limiting the algorithm's convergence and sample randomness.To address these limitations and evaluate the uncertainty of AVO inversion,the IADR-Gibbs algorithm is proposed,which incorporates the Independent Adaptive Delayed Rejection(IADR)algorithm with the Gibbs sampling algorithm.Grounded in Bayesian theory,the new algorithm introduces support points to construct a proposal distribution of non-parametric distribution and reselects the rejected samples according to the Delayed Rejection(DR)strategy.Rejected samples are then added to the support points to update the proposal distribution function adaptively.The equation rewriting method and the IADR-Gibbs algorithm improve the accuracy and robustness of AVO inversion.The effectiveness and applicability of the proposed method are validated through synthetic gather tests and practical data applications.展开更多
Predictions of fluid distribution,stress field,and natural fracture are essential for exploiting unconventional shale gas reservoirs.Given the high likelihood of tilted fractures in subsurface formations,this study fo...Predictions of fluid distribution,stress field,and natural fracture are essential for exploiting unconventional shale gas reservoirs.Given the high likelihood of tilted fractures in subsurface formations,this study focuses on simultaneous seismic inversion to estimate fluid bulk modulus,effective stress parameter,and fracture density in the tilted transversely isotropic(TTI)medium.In this article,a novel PP-wave reflection coefficient approximation equation is first derived based on the constructed TTI stiffness matrix incorporating fracture density,effective stress parameter,and fluid bulk modulus.The high accuracy of the proposed equation has been demonstrated using an anisotropic two-layer model.Furthermore,a stepwise seismic inversion strategy with the L_(P) quasi-norm sparsity constraint is implemented to obtain the anisotropic and isotropic parameters.Three synthetic model tests with varying signal-to-noise ratios(SNRs)confirm the method's feasibility and noise robustness.Ultimately,the proposed method is applied to a 3D fractured shale gas reservoir in the Sichuan Basin,China.The results have effectively characterized shale gas distribution,stress fields,and tilted natural fractures,with validation from geological structures,well logs,and microseismic events.These findings can provide valuable guidance for hydraulic fracturing development,enabling more reliable predictions of reservoir heterogeneity and completion quality.展开更多
Machine learning(ML)has strong potential for soil settlement prediction,but determining hyperparameters for ML models is often intricate and laborious.Therefore,we apply Bayesian optimization to determine the optimal ...Machine learning(ML)has strong potential for soil settlement prediction,but determining hyperparameters for ML models is often intricate and laborious.Therefore,we apply Bayesian optimization to determine the optimal hyperparameter combinations,enhancing the effectiveness of ML models for soil parameter inversion.The ML models are trained using numerical simulation data generated with the modified Cam-Clay(MCC)model in ABAQUS software,and their performance is evaluated using ground settlement monitoring data from an airport runway.Five optimized ML models—decision tree(DT),random forest(RF),support vector regression(SVR),deep neural network(DNN),and one-dimensional convolutional neural network(1D-CNN)—are compared in terms of their accuracy for soil parameter inversion and settlement prediction.The results indicate that Bayesian optimization efficiently utilizes prior knowledge to identify the optimal hyperparameters,significantly improving model performance.Among the evaluated models,the 1D-CNN achieves the highest accuracy in soil parameter inversion,generating settlement predictions that closely match real monitoring data.These findings demonstrate the effectiveness of the proposed approach for soil parameter inversion and settlement prediction,and reveal how Bayesian optimization can refine the model selection process.展开更多
Variation of reservoir physical properties can cause changes in its elastic parameters. However, this is not a simple linear relation. Furthermore, the lack of observations, data overlap, noise interference, and ideal...Variation of reservoir physical properties can cause changes in its elastic parameters. However, this is not a simple linear relation. Furthermore, the lack of observations, data overlap, noise interference, and idealized models increases the uncertainties of the inversion result. Thus, we propose an inversion method that is different from traditional statistical rock physics modeling. First, we use deterministic and stochastic rock physics models considering the uncertainties of elastic parameters obtained by prestack seismic inversion and introduce weighting coefficients to establish a weighted statistical relation between reservoir and elastic parameters. Second, based on the weighted statistical relation, we use Markov chain Monte Carlo simulations to generate the random joint distribution space of reservoir and elastic parameters that serves as a sample solution space of an objective function. Finally, we propose a fast solution criterion to maximize the posterior probability density and obtain reservoir parameters. The method has high efficiency and application potential.展开更多
The three parameters of P-wave velocity, S-wave velocity, and density have remarkable differences in conventional prestack inversion accuracy, so study of the consistency inversion of the ＂three parameters＂ is very ...The three parameters of P-wave velocity, S-wave velocity, and density have remarkable differences in conventional prestack inversion accuracy, so study of the consistency inversion of the ＂three parameters＂ is very important. In this paper, we present a new inversion algorithm and approach based on the in-depth analysis of the causes in their accuracy differences. With this new method, the inversion accuracy of the three parameters is improved synchronously by reasonable approximations and mutual constraint among the parameters. Theoretical model calculations and actual data applications with this method indicate that the three elastic parameters all have high inversion accuracy and maintain consistency, which also coincides with the theoretical model and actual data. This method has good application prospects.展开更多
Based on 310 horizontal-component digital seismograms recorded at 14 seismic stations in Shanxi Digital Seis-mograph Network, the inelastic attenuation coefficient in Shanxi region is studied. By the methods of Atkins...Based on 310 horizontal-component digital seismograms recorded at 14 seismic stations in Shanxi Digital Seis-mograph Network, the inelastic attenuation coefficient in Shanxi region is studied. By the methods of Atkinson and Moya, the site response of each station and several source parameters are obtained and the inversion results from both methods are compared and analyzed. The frequency-dependent inelastic attenuation coefficient Q is estimated as Q( f )=323.2 f 0.506. The site responses of 14 seismic stations do not show significant amplification, which is consistent with their basement on rock. We also found the dependence of corner frequency on seismic moment, seismic moment on stress drop, source radius on stress drop.展开更多
The complexity of the seismicity pattern for the subduction zone along the oceanic plate triggered the outer rise events and revealed cyclic tectonic deformation conditions along the plate subduction zones.The outer r...The complexity of the seismicity pattern for the subduction zone along the oceanic plate triggered the outer rise events and revealed cyclic tectonic deformation conditions along the plate subduction zones.The outer rise earthquakes have been observed along the Sunda arc,following the estimated rupture area of the 2005 M_(W)8.6 Nias earthquakes.Here,we used kinematic waveform inversion(KIWI)to obtain the source parameters of the 14 May 2021 M_(W)6.6 event off the west coast of northern Sumatra and to define the fault plane that triggered this outer rise event.The KIWI algorithm allows two types of seismic source to be configured:the moment tensor model to describe the type of shear with six moment tensor components and the Eikonal model for the rupture of pure double-couple sources.This method was chosen for its flexibility to be applied for different sources of seismicity and also for the automated full-moment tensor solution with real-time monitoring.We used full waveform traces from 8 broadband seismic stations within 1000 km epicentral distances sourced from the Incorporated Research Institutions for Seismology(IRIS-IDA)and Geofon GFZ seismic record databases.The initial origin time and hypocenter values are obtained from the IRIS-IDA.The synthetic seismograms used in the inversion process are based on the existing regional green function database model and were accessed from the KIWI Tools Green's Function Database.The obtained scalar seismic moment value is 1.18×10^(19)N·m,equivalent to a moment magnitude M_(W)6.6.The source parameters are 140°,44°,and−99°for the strike,dip,and rake values at a centroid depth of 10.2 km,indicating that this event is a normal fault earthquake that occurred in the outer rise area.The outer rise events with normal faults typically occur at the shallow part of the plate,with nodal-plane dips predominantly in the range of 30°-60°on the weak oceanic lithosphere due to hydrothermal alteration.The stress regime around the plate subduction zone varies both temporally and spatially due to the cyclic influences of megathrust earthquakes.Tensional outer rise earthquakes tend to occur after the megathrust events.The relative timing of these events is not known due to the viscous relaxation of the down going slab and poroelastic response in the trench slope region.The occurrence of the 14 May 2021 earthquake shows the seismicity in the outer rise region in the strongly coupled Sunda arc subduction zone due to elastic bending stress within the duration of the seismic cycle.展开更多
With a more complex pore structure system compared with clastic rocks, carbonate rocks have not yet been well described by existing conventional rock physical models concerning the pore structure vagary as well as the...With a more complex pore structure system compared with clastic rocks, carbonate rocks have not yet been well described by existing conventional rock physical models concerning the pore structure vagary as well as the influence on elastic rock properties. We start with a discussion and an analysis about carbonate rock pore structure utilizing rock slices. Then, given appropriate assumptions, we introduce a new approach to modeling carbonate rocks and construct a pore structure algorithm to identify pore structure mutation with a basis on the Gassmann equation and the Eshelby-Walsh ellipsoid inclusion crack theory. Finally, we compute a single well＇s porosity using this new approach with full wave log data and make a comparison with the predicted result of traditional method and simultaneously invert for reservoir parameters. The study results reveal that the rock pore structure can significantly influence the rocks＇ elastic properties and the predicted porosity error of the new modeling approach is merely 0.74%. Therefore, the approach we introduce can effectively decrease the predicted error of reservoir parameters.展开更多
For density inversion of gravity anomaly data, once the inversion method is determined, the main factors affecting the inversion result are the inversion parameters and subdivision scheme. A set of reasonable inversio...For density inversion of gravity anomaly data, once the inversion method is determined, the main factors affecting the inversion result are the inversion parameters and subdivision scheme. A set of reasonable inversion parameters and subdivision scheme can, not only improve the inversion process efficiency, but also ensure inversion result accuracy. The gravity inversion method based on correlation searching and the golden section algorithm is an effective potential field inversion method. It can be used to invert 2D and 3D physical properties with potential data observed on flat or rough surfaces. In this paper, we introduce in detail the density inversion principles based on correlation searching and the golden section algorithm. Considering that the gold section algorithm is not globally optimized. we present a heuristic method to ensure the inversion result is globally optimized. With a series of model tests, we systematically compare and analyze the inversion result efficiency and accuracy with different parameters. Based on the model test results, we conclude the selection principles for each inversion parameter with which the inversion accuracy can be obviously improved.展开更多
As a high quality seismic imaging method, full waveform inversion （FWI） can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing i...As a high quality seismic imaging method, full waveform inversion （FWI） can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing is computationally expensive, especially for the three-dimension complex medium inversion. Introducing blended source technology into the frequency-domain FWI can greatly reduce the computational burden and improve the efficiency of the inversion. However, this method has two issues： first, crosstalk noise is caused by interference between the sources involved in the encoding, resulting in an inversion result with some artifacts; second, it is more sensitive to ambient noise compared to conventional FWI, therefore noisy data results in a poor inversion. This paper introduces a frequency-group encoding method to suppress crosstalk noise, and presents a frequency- domain auto-adapting FWI based on source-encoding technology. The conventional FWI method and source-encoding based FWI method are combined using an auto-adapting mechanism. This improvement can both guarantee the quality of the inversion result and maximize the inversion efficiency.展开更多
Cauchy priori distribution-based Bayesian AVO reflectivity inversion may lead to sparse estimates that are sensitive to large reflectivities. For the inversion, the computation of the covariance matrix and regularized...Cauchy priori distribution-based Bayesian AVO reflectivity inversion may lead to sparse estimates that are sensitive to large reflectivities. For the inversion, the computation of the covariance matrix and regularized terms requires prior estimation of model parameters, which makes the iterative inversion weakly nonlinear. At the same time, the relations among the model parameters are assumed linear. Furthermore, the reflectivities, the results of the inversion, or the elastic parameters with cumulative error recovered by integrating reflectivities are not well suited for detecting hydrocarbons and fuids. In contrast, in Bayesian linear AVO inversion, the elastic parameters can be directly extracted from prestack seismic data without linear assumptions for the model parameters. Considering the advantages of the abovementioned methods, the Bayesian AVO reflectivity inversion process is modified and Cauchy distribution is explored as a prior probability distribution and the time-variant covariance is also considered. Finally, we propose a new method for the weakly nonlinear AVO waveform inversion. Furthermore, the linear assumptions are abandoned and elastic parameters, such as P-wave velocity, S-wave velocity, and density, can be directly recovered from seismic data especially for interfaces with large reflectivities. Numerical analysis demonstrates that all the elastic parameters can be estimated from prestack seismic data even when the signal-to-noise ratio of the seismic data is low.展开更多
Multiturn coils is an effective transmitter for transient electromagnetic method(TEM) used in narrow space and complex terrain at presently. However, its high mutual inductance coupling and long turn-off time affect t...Multiturn coils is an effective transmitter for transient electromagnetic method(TEM) used in narrow space and complex terrain at presently. However, its high mutual inductance coupling and long turn-off time affect the quality of later data processing and interpretation. Compared with multiturn coils, the new conical source has low mutual inductance and short turn-off time. Based on the superposition principle, we use Hankel transform and numerical filtering method for forward modelling of the conical source field in the layered-media and explore TEM characteristics excited by this source. We apply improved damped least square inversion to integrated transient electromagnetic(TEM) data. We first invert the induced voltage into similar resistivity and apparent depth, and then use the inverted results as input parameters in the initial model and transform the apparent resistivity data into the frequency domain. Then, damped least square inversion is performed in the frequency domain using the initial model. Subsequently, we use automated model building to search for the extremes and inflection points in the resistivity–depth data that are treated as critical layer parameters. The inversion of theoretical and observed data suggests that the method modifies the resistivity and depth and yields a model of the underground layers.展开更多
基金funded by the National Natural Science Foundation of China(42174011).
文摘With the continuous improvement of the accuracy of geodetic deformation data,the inversion of seismic source parameters puts forward a higher demand for nonlinear inversion algorithms.In this research,an improved Sparrow Search Algorithm(SSA)is proposed for the seismic source parameter inversion problem.By replacing the original population generation in the improved algorithm with Latin hypercubic sampling,the Sparrow Search Algorithm reduces the repetition of samples in the population initialization.Subsequently,the algorithm introduces adaptive weights in the discoverer generation phase of the sparrow algorithm and combines the Levy flight strategy to make the algorithm more comprehensive and improve the search accuracy during the whole iteration process.Therefore,the improved Latin hypercube-based sparrow search algorithm(ILHSSA)has better advantages in terms of iterative convergence speed and stability.In order to verify the performance of ILHSSA,the basic genetic algorithm(GA)and sparrow search algorithm(SSA)are examined and compared with ILHSSA by simulated earthquakes of two different earthquake types.The simulation experiments show that the improved algorithm ILHSSA outperforms SSA in accuracy and stability.Compared with the GA algorithm,ILHSSA can achieve the same inversion accuracy as GA,and it even surpasses GA in inversion speed and the inversion results of some parameters,demonstrating better stability.Finally,the improved algorithm is used for the 2017 Bodrum-Cos earthquake and the 2016 Amatrice earthquake in Italy.The inversion results all reflect the practicality and reliability of the improved algorithm.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.12274348 and 12004335)the National Key Research and Development Program of China(Grant No.2024YFC2813800)。
文摘Presented in this study is a novel method for estimating the depth of single underwater source in shallow water,utilizing vector sensors.The approach leverages the depth distribution of the broadband Stokes parameters to estimate source depth accurately.Unlike traditional matched field processing(MFP)and matched mode processing(MMP),the proposed approach can estimate source depth directly from the data received by sensors without requiring complete environmental information.Firstly,the broadband Stokes parameters(BSP)are established using the normal mode theory.Then the nonstationary phase approximation is used to simplify the theoretical derivation,which is necessary when dealing with broadband integrals.Additionally,range terms of the BSP are eliminated by normalization.By analyzing the depth distribution of the normalized broadband Stokes parameters(NBSP),it is found that the NBSP exhibit extreme values at the source depth,which can be used for source depth estimation.So the proposed depth estimation method is based on searching the peaks of the NBSP.Simulations show that this method is effective in relatively simple shallow water environments.Finally,the effect of source range,frequency bandwidth,sound speed profile(SSP),water depth,and signal-to-noise ratio(SNR)are studied.The findings indicate that the proposed method can accurately estimate the source depth when the SNR is greater than-5 d B and does not need to consider model mismatch issues.Additionally,variations in environmental parameters have minimal impact on estimation accuracy.Compared to MFP,the proposed method requires a higher SNR,but demonstrates superior robustness against fluctuations in environmental parameters.
基金The Project Advanced Buildings technology in a Dense Urban Environment of Hong Kong Polytechnic University the Application Lab of Digital Seismic Wave Data in Center for Analysis and Prediction+1 种基金 China Seismological Bureau and the Earthquake Prediction
文摘An earthquake of M S=6.9 occurred in Gonghe County, Qinghai Province, China on April 26, 1990.This earthquake was followed by three larger aftershocks of M S=5.5 on May 7, 1990, M S=6.0 on Jan.3, 1994, and M S=5.7 on Feb.16, 1994, consecutively. The moment tensors of these earthquakes as function of time were obtained by the technique of moment tensor inversion in frequency domain . The results inverted indicate that these earthquakes had a very similar focal mechanism of predominantly reverse faulting on a plane striking NWW, dipping to SSW.The scalar seismic moments of these earthquakes are M 0=9.4×10 18 Nm for the M S=6.9 event, 8.0×10 16 Nm for the M S=5.5 event, 4.9×10 17 Nm for the M S =6.0 event and 2.9×10 17 Nm for the M S=5.7 event, respectively. The results inverted also show that the source processes of these events were significantly different. The main shock had a very complex process, consisting of two distinct sub events with comparable sizes. The first sub event occurred in the first 12s, having a seismic moment of 4.7×10 18 Nm, and the second one continued from 31s to 41s, having a seismic moment of 2.5×10 18 Nm. In addition, a much smaller sub event, having a seismic moment of about 2.1×10 18 Nm, may exist in the interval of 12 s and 31 s, In contrast, the source processes of the three aftershocks are quite simple. The source time function of each of aftershocks is a single impulse, suggestting that each of aftershocks consists of a mainly uninterrupted rupture. The rise times and total rupture durations are 4 s and 11 s for the M S=5.5 event, 6 s and 16 s for the M S= 6.0 event and 6 s and 13 s for the M S=5.7 event, respectively.
基金the National Natural Science Foundation of China(Nos.42174011and 41874001).
文摘The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault parameters inversion has nonlinear characteristics, and the gradient-based optimizer(GBO) has the characteristics of fast convergence speed and falling into local optimum hardly. This paper applies GBO algorithm to simulated earthquakes and real LuShan earthquakes in the nonlinear inversion of the Okada model to obtain the source parameters. The simulated earthquake experiment results show that the algorithm is stable, and the seismic source parameters obtained by GBO are slightly closer to the true value than the multi peak particle swarm optimization(MPSO). In the 2013 LuShan earthquake experiment, the root mean square error between the deformation after forwarding of fault parameters obtained by the introduced GBO algorithm and the surface observation deformation was 3.703 mm, slightly better than 3.708 mm calculated by the MPSO. Moreover, the inversion result of GBO algorithm is better than MPSO algorithm in stability. The above results show that the introduced GBO algorithm has a certain practical application value in seismic fault source parameters inversion.
基金the National Natural Science Foundation of China(No.12072290).
文摘The seismoacoustic analysis method has broad potential applications to source parameter estimation for near-surface explosion events such as industrial explosions and terrorist attacks.In this study,current models were improved by modifying the acoustic model and adopting the Bayesian Markov-chain-Monte-Carlo inversion method.The source parameters of near-surface small-yield chemical explosions were analyzed via the improved seismoacoustic analysis model and by the estimation accuracy of seismoacoustic joint inversion.Estimation and analysis results showed that the improved seismoacoustic analysis model considered ground shock coupling and the impact of explosion products ejecting from the surface so that the improved acoustic impulse relation was more consistent with the measured data than the Ford impulse relation.It is suitable for deep-burial,shallow-burial,and near-surface aerial explosions.Furthermore,trade-off relationships were declined through the application of the improved model to source parameter inversion for near-surface small-yield chemical explosions,and source parameter estimation accuracy was improved.
基金Project of Center for Analysis and Prediction of CEA, Key Project of State Science and Technology in the Tenth Five-year Plan (02-03-04), Joint Seismological Science Foundation of China (604022), and Natural Science Founda-tion of Beijing Municipality (8022009).
文摘According to the geological structural features, Beijing and the adjacent areas can be divided into two regions of plain in the east and mountain in the west. Among the stations covered by the telemetered digital seismic station network of Earthquake Administration of Beijing Municipality, the stations in the plain area are all borehole ones and the stations in the western mountainous region are all located on the surface bedrock. In the paper, 511 wave- form data recorded by the network from Oct. 2001 to Oct. 2004 are used in the researches for the entire Beijing region, the western mountainous region and the eastern plain area, respectively. The Q values are calculated for each area by Atkinson′s method and compared with the existed data. The reliability of the Q values and the reasons for the difference in the Q values are also discussed. Then, the source parameters and site response are inverted by the Moya′s method, in which two models are used. The first model uses the Q values, earthquakes and stations in the sub-areas and the second model uses the Q values, earthquakes and stations in the entire Beijing region. The results indicate that the source parameters and site responses obtained by two models are basically consistent with each other. It also indicates that the source parameters obtained by these methods are not affected by the size of station network.
基金supported by Innovation Project of Chinese Academy of Sciences
文摘A method for simultaneous determination of mixed model parameters,which have different physical dimensions or different responses to data,is presented.Mixed parameter estimation from observed data within a single model space shows instabilities and trade-offs of the solutions. We separate the model space into N-subspaces based on their physical properties or computational convenience and solve the N-subspaces systems by damped least-squares and singular-value decomposition. Since the condition number of each subsystem is smaller than that of the single global system,the approach can greatly increase the stability of the inversion. We also introduce different damping factors into the subsystems to reduce the tradeoffs between the different parameters. The damping factors depend on the conditioning of the subsystems and may be adequately chosen in a range from 0.1 % to 10 % of the largest singular value. We illustrate the method with an example of simultaneous determination of source history,source geometry,and hypocentral location from regional seismograms,although it is applicable to any geophysical inversion.
基金the sponsorship of the Key Technology for Geophysical Prediction of Ultra-Deep Carbonate Reservoirs(P24240)the National Natural Science Foundation of China(U24B2020)the National Science and Technology Major Project of China for New Oil and Gas Exploration and Development(Grant No.2024ZD1400102)。
文摘Fluid identification and anisotropic parameters characterization are crucial for shale reservoir exploration and development.However,the anisotropic reflection coefficient equation,based on the transverse isotropy with a vertical axis of symmetry(VTI)medium assumption,involves numerous parameters to be inverted.This complexity reduces its stability and impacts the accuracy of seismic amplitude variation with offset(AVO)inversion results.In this study,a novel anisotropic equation that includes the fluid term and Thomsen anisotropic parameters is rewritten,which reduces the equation's dimensionality and increases its stability.Additionally,the traditional Markov Chain Monte Carlo(MCMC)inversion algorithm exhibits a high rejection rate for random samples and relies on known parameter distributions such as the Gaussian distribution,limiting the algorithm's convergence and sample randomness.To address these limitations and evaluate the uncertainty of AVO inversion,the IADR-Gibbs algorithm is proposed,which incorporates the Independent Adaptive Delayed Rejection(IADR)algorithm with the Gibbs sampling algorithm.Grounded in Bayesian theory,the new algorithm introduces support points to construct a proposal distribution of non-parametric distribution and reselects the rejected samples according to the Delayed Rejection(DR)strategy.Rejected samples are then added to the support points to update the proposal distribution function adaptively.The equation rewriting method and the IADR-Gibbs algorithm improve the accuracy and robustness of AVO inversion.The effectiveness and applicability of the proposed method are validated through synthetic gather tests and practical data applications.
基金financially supported by the Natural Science Foundation of Sichuan Province(Grant Nos.2023NSFSC0767 and2024NSFSC0809)the China Postdoctoral Science Foundation(Grant No.2024MF750281)the Postdoctoral Fellowship Program of CPSF(Grant No.GZC20230326)。
文摘Predictions of fluid distribution,stress field,and natural fracture are essential for exploiting unconventional shale gas reservoirs.Given the high likelihood of tilted fractures in subsurface formations,this study focuses on simultaneous seismic inversion to estimate fluid bulk modulus,effective stress parameter,and fracture density in the tilted transversely isotropic(TTI)medium.In this article,a novel PP-wave reflection coefficient approximation equation is first derived based on the constructed TTI stiffness matrix incorporating fracture density,effective stress parameter,and fluid bulk modulus.The high accuracy of the proposed equation has been demonstrated using an anisotropic two-layer model.Furthermore,a stepwise seismic inversion strategy with the L_(P) quasi-norm sparsity constraint is implemented to obtain the anisotropic and isotropic parameters.Three synthetic model tests with varying signal-to-noise ratios(SNRs)confirm the method's feasibility and noise robustness.Ultimately,the proposed method is applied to a 3D fractured shale gas reservoir in the Sichuan Basin,China.The results have effectively characterized shale gas distribution,stress fields,and tilted natural fractures,with validation from geological structures,well logs,and microseismic events.These findings can provide valuable guidance for hydraulic fracturing development,enabling more reliable predictions of reservoir heterogeneity and completion quality.
基金supported by the National Natural Science Foundation of China(Nos.52378419 and 52478368).
文摘Machine learning(ML)has strong potential for soil settlement prediction,but determining hyperparameters for ML models is often intricate and laborious.Therefore,we apply Bayesian optimization to determine the optimal hyperparameter combinations,enhancing the effectiveness of ML models for soil parameter inversion.The ML models are trained using numerical simulation data generated with the modified Cam-Clay(MCC)model in ABAQUS software,and their performance is evaluated using ground settlement monitoring data from an airport runway.Five optimized ML models—decision tree(DT),random forest(RF),support vector regression(SVR),deep neural network(DNN),and one-dimensional convolutional neural network(1D-CNN)—are compared in terms of their accuracy for soil parameter inversion and settlement prediction.The results indicate that Bayesian optimization efficiently utilizes prior knowledge to identify the optimal hyperparameters,significantly improving model performance.Among the evaluated models,the 1D-CNN achieves the highest accuracy in soil parameter inversion,generating settlement predictions that closely match real monitoring data.These findings demonstrate the effectiveness of the proposed approach for soil parameter inversion and settlement prediction,and reveal how Bayesian optimization can refine the model selection process.
基金supported by the National Science and Technology Major Project(No.2011 ZX05007-006)the 973 Program of China(No.2013CB228604)the Major Project of Petrochina(No.2014B-0610)
文摘Variation of reservoir physical properties can cause changes in its elastic parameters. However, this is not a simple linear relation. Furthermore, the lack of observations, data overlap, noise interference, and idealized models increases the uncertainties of the inversion result. Thus, we propose an inversion method that is different from traditional statistical rock physics modeling. First, we use deterministic and stochastic rock physics models considering the uncertainties of elastic parameters obtained by prestack seismic inversion and introduce weighting coefficients to establish a weighted statistical relation between reservoir and elastic parameters. Second, based on the weighted statistical relation, we use Markov chain Monte Carlo simulations to generate the random joint distribution space of reservoir and elastic parameters that serves as a sample solution space of an objective function. Finally, we propose a fast solution criterion to maximize the posterior probability density and obtain reservoir parameters. The method has high efficiency and application potential.
基金sponsored by the National Major Program (No. 2011ZX05006-006)the 973 Program of China (No. 2011CB201104)Technical Research of Elastic Flooding Boundary and Well Network Optimization at the Development Late Stage of Low Permeable Oil Field (No. 2011ZX05009)
文摘The three parameters of P-wave velocity, S-wave velocity, and density have remarkable differences in conventional prestack inversion accuracy, so study of the consistency inversion of the ＂three parameters＂ is very important. In this paper, we present a new inversion algorithm and approach based on the in-depth analysis of the causes in their accuracy differences. With this new method, the inversion accuracy of the three parameters is improved synchronously by reasonable approximations and mutual constraint among the parameters. Theoretical model calculations and actual data applications with this method indicate that the three elastic parameters all have high inversion accuracy and maintain consistency, which also coincides with the theoretical model and actual data. This method has good application prospects.
基金National Scientific Research Target Project of China (2001-BA601B01-01-05) during the tenth Five-year Plan.
文摘Based on 310 horizontal-component digital seismograms recorded at 14 seismic stations in Shanxi Digital Seis-mograph Network, the inelastic attenuation coefficient in Shanxi region is studied. By the methods of Atkinson and Moya, the site response of each station and several source parameters are obtained and the inversion results from both methods are compared and analyzed. The frequency-dependent inelastic attenuation coefficient Q is estimated as Q( f )=323.2 f 0.506. The site responses of 14 seismic stations do not show significant amplification, which is consistent with their basement on rock. We also found the dependence of corner frequency on seismic moment, seismic moment on stress drop, source radius on stress drop.
基金supported by the National Natural Science Foundation of China(Grant No.42130312)。
文摘The complexity of the seismicity pattern for the subduction zone along the oceanic plate triggered the outer rise events and revealed cyclic tectonic deformation conditions along the plate subduction zones.The outer rise earthquakes have been observed along the Sunda arc,following the estimated rupture area of the 2005 M_(W)8.6 Nias earthquakes.Here,we used kinematic waveform inversion(KIWI)to obtain the source parameters of the 14 May 2021 M_(W)6.6 event off the west coast of northern Sumatra and to define the fault plane that triggered this outer rise event.The KIWI algorithm allows two types of seismic source to be configured:the moment tensor model to describe the type of shear with six moment tensor components and the Eikonal model for the rupture of pure double-couple sources.This method was chosen for its flexibility to be applied for different sources of seismicity and also for the automated full-moment tensor solution with real-time monitoring.We used full waveform traces from 8 broadband seismic stations within 1000 km epicentral distances sourced from the Incorporated Research Institutions for Seismology(IRIS-IDA)and Geofon GFZ seismic record databases.The initial origin time and hypocenter values are obtained from the IRIS-IDA.The synthetic seismograms used in the inversion process are based on the existing regional green function database model and were accessed from the KIWI Tools Green's Function Database.The obtained scalar seismic moment value is 1.18×10^(19)N·m,equivalent to a moment magnitude M_(W)6.6.The source parameters are 140°,44°,and−99°for the strike,dip,and rake values at a centroid depth of 10.2 km,indicating that this event is a normal fault earthquake that occurred in the outer rise area.The outer rise events with normal faults typically occur at the shallow part of the plate,with nodal-plane dips predominantly in the range of 30°-60°on the weak oceanic lithosphere due to hydrothermal alteration.The stress regime around the plate subduction zone varies both temporally and spatially due to the cyclic influences of megathrust earthquakes.Tensional outer rise earthquakes tend to occur after the megathrust events.The relative timing of these events is not known due to the viscous relaxation of the down going slab and poroelastic response in the trench slope region.The occurrence of the 14 May 2021 earthquake shows the seismicity in the outer rise region in the strongly coupled Sunda arc subduction zone due to elastic bending stress within the duration of the seismic cycle.
基金sponsored by the National Nature Science Foundation of China (Grant No.40904034 and 40839905)
文摘With a more complex pore structure system compared with clastic rocks, carbonate rocks have not yet been well described by existing conventional rock physical models concerning the pore structure vagary as well as the influence on elastic rock properties. We start with a discussion and an analysis about carbonate rock pore structure utilizing rock slices. Then, given appropriate assumptions, we introduce a new approach to modeling carbonate rocks and construct a pore structure algorithm to identify pore structure mutation with a basis on the Gassmann equation and the Eshelby-Walsh ellipsoid inclusion crack theory. Finally, we compute a single well＇s porosity using this new approach with full wave log data and make a comparison with the predicted result of traditional method and simultaneously invert for reservoir parameters. The study results reveal that the rock pore structure can significantly influence the rocks＇ elastic properties and the predicted porosity error of the new modeling approach is merely 0.74%. Therefore, the approach we introduce can effectively decrease the predicted error of reservoir parameters.
基金supported by Specialized Research Fund for the Doctoral Program of Higher Education of China(20110022120004)the Fundamental Research Funds for the Central Universities
文摘For density inversion of gravity anomaly data, once the inversion method is determined, the main factors affecting the inversion result are the inversion parameters and subdivision scheme. A set of reasonable inversion parameters and subdivision scheme can, not only improve the inversion process efficiency, but also ensure inversion result accuracy. The gravity inversion method based on correlation searching and the golden section algorithm is an effective potential field inversion method. It can be used to invert 2D and 3D physical properties with potential data observed on flat or rough surfaces. In this paper, we introduce in detail the density inversion principles based on correlation searching and the golden section algorithm. Considering that the gold section algorithm is not globally optimized. we present a heuristic method to ensure the inversion result is globally optimized. With a series of model tests, we systematically compare and analyze the inversion result efficiency and accuracy with different parameters. Based on the model test results, we conclude the selection principles for each inversion parameter with which the inversion accuracy can be obviously improved.
基金financially supported by the National Natural Science Foundation of China(No.41074075/D0409)the National Science and Technology Major Project(No.2011ZX05025-001-04)
文摘As a high quality seismic imaging method, full waveform inversion （FWI） can accurately reconstruct the physical parameter model for the subsurface medium. However, application of the FWI in seismic data processing is computationally expensive, especially for the three-dimension complex medium inversion. Introducing blended source technology into the frequency-domain FWI can greatly reduce the computational burden and improve the efficiency of the inversion. However, this method has two issues： first, crosstalk noise is caused by interference between the sources involved in the encoding, resulting in an inversion result with some artifacts; second, it is more sensitive to ambient noise compared to conventional FWI, therefore noisy data results in a poor inversion. This paper introduces a frequency-group encoding method to suppress crosstalk noise, and presents a frequency- domain auto-adapting FWI based on source-encoding technology. The conventional FWI method and source-encoding based FWI method are combined using an auto-adapting mechanism. This improvement can both guarantee the quality of the inversion result and maximize the inversion efficiency.
基金supported by the National High-Tech Research and Development Program of China(863 Program)(No.2008AA093001)
文摘Cauchy priori distribution-based Bayesian AVO reflectivity inversion may lead to sparse estimates that are sensitive to large reflectivities. For the inversion, the computation of the covariance matrix and regularized terms requires prior estimation of model parameters, which makes the iterative inversion weakly nonlinear. At the same time, the relations among the model parameters are assumed linear. Furthermore, the reflectivities, the results of the inversion, or the elastic parameters with cumulative error recovered by integrating reflectivities are not well suited for detecting hydrocarbons and fuids. In contrast, in Bayesian linear AVO inversion, the elastic parameters can be directly extracted from prestack seismic data without linear assumptions for the model parameters. Considering the advantages of the abovementioned methods, the Bayesian AVO reflectivity inversion process is modified and Cauchy distribution is explored as a prior probability distribution and the time-variant covariance is also considered. Finally, we propose a new method for the weakly nonlinear AVO waveform inversion. Furthermore, the linear assumptions are abandoned and elastic parameters, such as P-wave velocity, S-wave velocity, and density, can be directly recovered from seismic data especially for interfaces with large reflectivities. Numerical analysis demonstrates that all the elastic parameters can be estimated from prestack seismic data even when the signal-to-noise ratio of the seismic data is low.
基金supported by the National Natural Science Foundation of China(Nos.41564001,41674133,41572185,and 41604104)the Distinguished Young Talent Foundation of Jiangxi Province(No.20171BCB23068)
文摘Multiturn coils is an effective transmitter for transient electromagnetic method(TEM) used in narrow space and complex terrain at presently. However, its high mutual inductance coupling and long turn-off time affect the quality of later data processing and interpretation. Compared with multiturn coils, the new conical source has low mutual inductance and short turn-off time. Based on the superposition principle, we use Hankel transform and numerical filtering method for forward modelling of the conical source field in the layered-media and explore TEM characteristics excited by this source. We apply improved damped least square inversion to integrated transient electromagnetic(TEM) data. We first invert the induced voltage into similar resistivity and apparent depth, and then use the inverted results as input parameters in the initial model and transform the apparent resistivity data into the frequency domain. Then, damped least square inversion is performed in the frequency domain using the initial model. Subsequently, we use automated model building to search for the extremes and inflection points in the resistivity–depth data that are treated as critical layer parameters. The inversion of theoretical and observed data suggests that the method modifies the resistivity and depth and yields a model of the underground layers.
