Total organic carbon (TOC) prediction with elastic parameter inversions has been widely used in the identification and evaluation of source rocks. However, the elastic parameters used to predict TOC are not only deter...Total organic carbon (TOC) prediction with elastic parameter inversions has been widely used in the identification and evaluation of source rocks. However, the elastic parameters used to predict TOC are not only determined by TOC but also depend on the other physical properties of source rocks. Besides, the TOC prediction with the elastic parameters inversion is an indirect method based on the statistical relationship obtained from well logs and experiment data. Therefore, we propose a rock physics model and define a TOC indicator mainly affected by TOC to predict TOC directly. The proposed rock physics model makes the equivalent elastic moduli of source rocks parameterized by the TOC indicator. Combining the equivalent elastic moduli of source rocks and Gray’s approximation leads to a novel linearized approximation of the P-wave reflection coefficient incorporating the TOC indicator. Model examples illustrate that the novel reflectivity approximation well agrees with the exact Zoeppritz equation until incident angles reach 40°. Convoluting the novel P-wave reflection approximation with seismic wavelets as the forward solver, an AVO inversion method based on the Bayesian theory is proposed to invert the TOC indicator with seismic data. The synthetic examples and field tests validate the feasibility and stability of the proposed AVO inversion approach. Using the inversion results of the TOC indicator, TOC is directly and accurately estimated in the target area.展开更多
Although the ambiguity of seismic inversion is widely recognized in both theory and practice, so far as a concrete inversion example is concerned, there is not any objective, controllable method or any standard for ho...Although the ambiguity of seismic inversion is widely recognized in both theory and practice, so far as a concrete inversion example is concerned, there is not any objective, controllable method or any standard for how to evaluate and determine its ambiguity and reliability, especially for the high frequency components beyond the effective seismic frequency band. Taking log-constrained impedance inversion as an example, a new appraisal method is proposed on the basis of analyzing a simple geological model. Firstly, the inverted impedance model is transformed to a reflection coefficient series. Secondly, the maximum effective frequency of the real seismic data is chosen as a cutoff point and the reflection coefficient series is decomposed into two components by low-pass and high-pass filters. Thirdly, the geometrical reflection characteristics of the high-frequency components and that of the real seismic data are compared and analyzed. Then, the reliability of the inverted impedance model is appraised according to the similarity of geometrical characteristics between the high-frequency components and the real seismic data. The new method avoids some subjectivity in appraising the inverted result, and helps to enhance the reliability of reservoir prediction by impedance inversion technology.展开更多
Comprehensive inversion of logging and seismic data presented in this paper is a method to improve seismic data resolution. It involves using ample high-frequency information and complete low-frequency information of ...Comprehensive inversion of logging and seismic data presented in this paper is a method to improve seismic data resolution. It involves using ample high-frequency information and complete low-frequency information of known logging to make up for the lack of limited bandwidth of practical seismic recording, obtaining an approximate reflection coefficient sequence (or wave impedance) of high resolution by iterative inversion and providing more reliable seismic evidence for further lithologic interpretation and lateral tracking, correlation and prediction of thin reservoir. The comprehensive inversion can be realized in the following steps: (1) to establish an initial model of higher resolution; (2) to obtain wavelets, and (3) to constrain iterative inversion. The key to this inversion lies in building an initial model. It is assumed from our experience that when the initial model is properly given, iterative inversion can be quickly converged to the ideal result.展开更多
Deep learning has achieved great success in a variety of research fields and industrial applications.However,when applied to seismic inversion,the shortage of labeled data severely influences the performance of deep l...Deep learning has achieved great success in a variety of research fields and industrial applications.However,when applied to seismic inversion,the shortage of labeled data severely influences the performance of deep learning-based methods.In order to tackle this problem,we propose a novel seismic impedance inversion method based on a cycle-consistent generative adversarial network(Cycle-GAN).The proposed Cycle-GAN model includes two generative subnets and two discriminative subnets.Three kinds of loss,including cycle-consistent loss,adversarial loss,and estimation loss,are adopted to guide the training process.Benefit from the proposed structure,the information contained in unlabeled data can be extracted,and adversarial learning further guarantees that the prediction results share similar distributions with the real data.Moreover,a neural network visualization method is adopted to show that the proposed CNN model can learn more distinguishable features than the conventional CNN model.The robustness experiments on synthetic data sets show that the proposed method can achieve better performances than other methods in most cases.And the blind-well experiments on real seismic profiles show that the predicted impedance curve of the proposed method maintains a better correlation with the true impedance curve.展开更多
Seismic amplitude variation with offset(AVO) inversion is an important approach for quantitative prediction of rock elasticity,lithology and fluid properties.With Biot-Gassmann's poroelasticity,an improved statist...Seismic amplitude variation with offset(AVO) inversion is an important approach for quantitative prediction of rock elasticity,lithology and fluid properties.With Biot-Gassmann's poroelasticity,an improved statistical AVO inversion approach is proposed.To distinguish the influence of rock porosity and pore fluid modulus on AVO reflection coefficients,the AVO equation of reflection coefficients parameterized by porosity,rock-matrix moduli,density and fluid modulus is initially derived from Gassmann equation and critical porosity model.From the analysis of the influences of model parameters on the proposed AVO equation,rock porosity has the greatest influences,followed by rock-matrix moduli and density,and fluid modulus has the least influences among these model parameters.Furthermore,a statistical AVO stepwise inversion method is implemented to the simultaneous estimation of rock porosity,rock-matrix modulus,density and fluid modulus.Besides,the Laplace probability model and differential evolution,Markov chain Monte Carlo algorithm is utilized for the stochastic simulation within Bayesian framework.Models and field data examples demonstrate that the simultaneous optimizations of multiple Markov chains can achieve the efficient simulation of the posterior probability density distribution of model parameters,which is helpful for the uncertainty analysis of the inversion and sets a theoretical fundament for reservoir characterization and fluid discrimination.展开更多
Yushu Ms7.1 earthquake occurred on the Ganzi-Yushu fault zone, across which we carried out a joint relative-gravity and seismic-reflection survey, and then performed a gravity inversion constrained by the seismic-refl...Yushu Ms7.1 earthquake occurred on the Ganzi-Yushu fault zone, across which we carried out a joint relative-gravity and seismic-reflection survey, and then performed a gravity inversion constrained by the seismic-reflection result. Based on the data of complete Bouguer gravity anomaly and seismic reflection, we obtained a layered interface structure in deep crust down to Moho. Our study showed that the inversion could reveal the interfaces of strata along the survey profile and the directions of regional faults in two-dimension. From the characteristics of the observed topography of the Moho basement, we tentatively confirmed that the uplift of eastern edge of Qinghai-Tibet plateau was caused by the subduetion of the Indian plate.展开更多
To improve the accuracy of inversion results,geological facies distributions are considered as additional constraints in the inversion process.However,the geological facies itself also has its own uncertainty.In this ...To improve the accuracy of inversion results,geological facies distributions are considered as additional constraints in the inversion process.However,the geological facies itself also has its own uncertainty.In this paper,the initial sedimentary facies maps are obtained by integrated geological analysis from well data,seismic attributes,and deterministic inversion results.Then the fi rst iteration of facies-constrained seismic inversion is performed.According to that result and other data such as geological information,the facies distribution can be updated using cluster analysis.The next round of facies-constrained inversion can then be performed.This process will be repeated until the facies inconsistency or error before and after the inversion is minimized.It forms a new iterative facies-constrained seismic inversion technique.Compared with conventional facies-constrained seismic inversion,the proposed method not only can reduces the non-uniqueness of seismic inversion results but also can improves its resolution.As a consequence,the sedimentary facies will be more consistent with the geology.A practical application demonstrated that the superposition relationship of sand bodies could be better delineated based on this new seismic inversion technique.The result highly increases the understanding of reservoir connectivity and its accuracy,which can be used to guide further development.展开更多
We investigate the accuracy and robustness of moment tensor(MT)and stress inversion solutions derived from acoustic emissions(AEs)during the laboratory fracturing of prismatic Barre granite specimens.Pre-cut flaws in ...We investigate the accuracy and robustness of moment tensor(MT)and stress inversion solutions derived from acoustic emissions(AEs)during the laboratory fracturing of prismatic Barre granite specimens.Pre-cut flaws in the specimens introduce a complex stress field,resulting in a spatial and temporal variation of focal mechanisms.Specifically,we consider two experimental setups:(1)where the rock is loaded in compression to generate primarily shear-type fractures and(2)where the material is loaded in indirect tension to generate predominantly tensile-type fractures.In each test,we first decompose AE moment tensors into double-couple(DC)and non-DC terms and then derive unambiguous normal and slip vectors using k-means clustering and an unstructured damped stress inversion algorithm.We explore temporal and spatial distributions of DC and non-DC events at different loading levels.The majority of the DC and the tensile non-DC events cluster around the pre-cut flaws,where macro-cracks later develop.Results of stress inversion are verified against the stress field from finite element(FE)modeling.A good agreement is found between the experimentally derived and numerically simulated stress orientations.To the best of the authors’knowledge,this work presents the first case where stress inversion methodologies are validated by numerical simulations at laboratory scale and under highly heterogeneous stress distributions.展开更多
Using the technique of seismic moment tensor inversion, the source mechanisms of 10 earthquakes with Ms5.2that occurred in China from November 1996 to January 1998 were determined rapidly. The determined resultswere s...Using the technique of seismic moment tensor inversion, the source mechanisms of 10 earthquakes with Ms5.2that occurred in China from November 1996 to January 1998 were determined rapidly. The determined resultswere sent as 'Bulletins of Source Mechanism Parameters of Earthquakes' to the Seismic Regime Guards' Office,China Seismological Bureau, and the relevant provincial seismological bureaus. These bulletins have played rolein the fast response to large earthquakes.展开更多
Deep learning is widely used for seismic impedance inversion,but few work provides in-depth research and analysis on designing the architectures of deep neural networks and choosing the network hyperparameters.This pa...Deep learning is widely used for seismic impedance inversion,but few work provides in-depth research and analysis on designing the architectures of deep neural networks and choosing the network hyperparameters.This paper is dedicated to comprehensively studying on the significant aspects of deep neural networks that affect the inversion results.We experimentally reveal how network hyperparameters and architectures affect the inversion performance,and develop a series of methods which are proven to be effective in reconstructing high-frequency information in the estimated impedance model.Experiments demonstrate that the proposed multi-scale architecture is helpful to reconstruct more high-frequency details than a conventional network.Besides,the reconstruction of high-frequency information can be further promoted by introducing a perceptual loss and a generative adversarial network from the computer vision perspective.More importantly,the experimental results provide valuable references for designing proper network architectures in the seismic inversion problem.展开更多
On September 16,2021,a MS6.0 earthquake struck Luxian County,one of the shale gas blocks in the Southeastern Sichuan Basin,China.To understand the seismogenic environment and its mechanism,we inverted a fine three-dim...On September 16,2021,a MS6.0 earthquake struck Luxian County,one of the shale gas blocks in the Southeastern Sichuan Basin,China.To understand the seismogenic environment and its mechanism,we inverted a fine three-dimensional S-wave velocity model from ambient noise tomography using data from a newly deployed dense seismic array around the epicenter,by extracting and jointly inverting the Rayleigh phase and group velocities in the period of 1.6–7.2 s.The results showed that the velocity model varied significantly beneath different geological units.The Yujiasi syncline is characterized by low velocity at depths of~3.0–4.0 km,corresponding to the stable sedimentary layer in the Sichuan Basin.The eastern and western branches of the Huayingshan fault belt generally exhibit high velocities in the NE-SW direction,with a few local low-velocity zones.The Luxian MS6.0 earthquake epicenter is located at the boundary between the high-and low-velocity zones,and the earthquake sequences expand eastward from the epicenter at depths of 3.0–5.0 km.Integrated with the velocity variations around the epicenter,distribution of aftershock sequences,and focal mechanism solution,it is speculated that the seismogenic mechanism of the main shock might be interpreted as the reactivation of pre-existing faults by hydraulic fracturing.展开更多
A genetic algorithm of body waveform inversion is presented for better understanding of crustal and upper mantle structures with deep seismic sounding (DSS) waveform data. General reflection and transmission synthet...A genetic algorithm of body waveform inversion is presented for better understanding of crustal and upper mantle structures with deep seismic sounding (DSS) waveform data. General reflection and transmission synthetic seismogram algorithm, which is capable of calculating the response of thin alternating high and low velocity layers, is applied as a solution for forward modeling, and the genetic algorithm is used to find the optimal solution of the inverse problem. Numerical tests suggest that the method has the capability of resolving low-velocity layers, thin alternating high and low velocity layers, and noise suppression. Waveform inversion using P-wave records from Zeku, Xiahe and Lintao shots in the seismic wide-angle reflection/refraction survey along northeastern Qinghai-Xizang (Tibeteau) Plateau has revealed fine structures of the bottom of the upper crust and alternating layers in the middle/lower crust and topmost upper mantle.展开更多
Under the condition of thin interbeds with great lateral changes in terrestrial basins,a seismic meme inversion method is established based on the analysis of seismic sedimentology technology.The relationship between ...Under the condition of thin interbeds with great lateral changes in terrestrial basins,a seismic meme inversion method is established based on the analysis of seismic sedimentology technology.The relationship between seismic waveform and high-frequency well logs is established through dynamic clustering of seismic waveform to improve the vertical and horizontal resolution of inversion results;meanwhile,by constructing the Bayesian inversion framework of different seismic facies,the real facies controlled inversion is realized.The forward model verification results show that the seismic meme inversion can realize precise prediction of 3 m thick thin interbeds,proving the rationality and high precision of the method.The application in the Daqing placanticline shows that the seismic meme inversion could identify 2 m thin interbeds,and the coincidence rates of inversion results and drilling data were more than 80%.The seismic meme inversion method can improve the accuracy of reservoir prediction and provides a useful mean for thin interbeds prediction in terrestrial basins.展开更多
2-D velocity structure and tectonics of the crust and upper mantle is revealed by inversion of seismic refraction and wide-angle reflection traveltimes acquired along the profile L1 in the Changbaishan-Tianchi volcani...2-D velocity structure and tectonics of the crust and upper mantle is revealed by inversion of seismic refraction and wide-angle reflection traveltimes acquired along the profile L1 in the Changbaishan-Tianchi volcanic region. It is used in this study that seismic traveltime inversion for simultaneous determination of 2-D velocity and interface structure of the crust and upper mantle. The result shows that, under Changbaishan-Tianchi crater, there exists a low-velocity body in the shape of an inverted triangle, and the crustal reflecting boundaries and Moho all become lower by a varying margin of 2-6 km, forming a crustal root which is assumed to be the Changbaishan-Tianchi volcanic system. Finally, we make a comparison between our 2-D velocity model and the result from the studies by using trial-and-error forward modeling with SEIS83.展开更多
The firework algorithm(FWA) is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model...The firework algorithm(FWA) is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model show that the FWA performs well in both solution quality and efficiency. We apply the FWA in this study to crustal velocity structure inversion using regional seismic waveform data of central Gansu on the northeastern margin of the Qinghai-Tibet plateau. Seismograms recorded from the moment magnitude(MW) 5.4 Minxian earthquake enable obtaining an average crustal velocity model for this region. We initially carried out a series of FWA robustness tests in regional waveform inversion at the same earthquake and station positions across the study region,inverting two velocity structure models, with and without a low-velocity crustal layer; the accuracy of our average inversion results and their standard deviations reveal the advantages of the FWA for the inversion of regional seismic waveforms. We applied the FWA across our study area using three component waveform data recorded by nine broadband permanent seismic stations with epicentral distances ranging between 146 and 437 km. These inversion results show that the average thickness of the crust in this region is 46.75 km, while thicknesses of the sedimentary layer, and the upper, middle, and lower crust are 3.15,15.69, 13.08, and 14.83 km, respectively. Results also show that the P-wave velocities of these layers and the upper mantle are 4.47, 6.07, 6.12, 6.87, and 8.18 km/s,respectively.展开更多
With the objective of establishing the necessary conditions for 3D seismic data from mountainous areas in western China, we compared the application results of wave impedance technology in the lithology and exploratio...With the objective of establishing the necessary conditions for 3D seismic data from mountainous areas in western China, we compared the application results of wave impedance technology in the lithology and exploration of coal fields. First, we introduce principles and features of three kinds of inversion methods. i.e., Model-Based Inversion, Constrained Sparse Spike Inversion (CSSI) and Geology-Seismic Feature Inversion. Secondly, these inversion methods are contrasted in their application to 3D seismic data from some coalfields in western China. The main information provided by the research includes: improving the vertical resolution of coal deposit strata, inferring lateral variation of the lithology and predicting coal seams and their roof lithology. Finally, the comparison between the three methods shows that the model-based inversion has the higher resolution, while CSSI inversion has better waveform continuity. The geology-seismic feature inversion requires information from a large number of wells and many types of logging curves of good quality. All three methods can meet the requirements of seismic exploration for lithological exploration in coal fields.展开更多
In Recent years,seismic data have been widely used in seismic oceanography for the inversion of oceanic parameters represented by conductivity temperature depth(CTD).Using this technique,researchers can identify the w...In Recent years,seismic data have been widely used in seismic oceanography for the inversion of oceanic parameters represented by conductivity temperature depth(CTD).Using this technique,researchers can identify the water structure with high horizontal resolution,which compensates for the deficiencies of CTD data.However,conventional inversion methods are modeldriven,such as constrained sparse spike inversion(CSSI)and full waveform inversion(FWI),and typically require prior deterministic mapping operators.In this paper,we propose a novel inversion method based on a convolutional neural network(CNN),which is purely data-driven.To solve the problem of multiple solutions,we use stepwise regression to select the optimal attributes and their combination and take two-dimensional images of the selected attributes as input data.To prevent vanishing gradients,we use the rectified linear unit(ReLU)function as the activation function of the hidden layer.Moreover,the Adam and mini-batch algorithms are combined to improve stability and efficiency.The inversion results of field data indicate that the proposed method is a robust tool for accurately predicting oceanic parameters.展开更多
In order to improve the fine structure inversion ability of igneous rocks for the exploration of underlying strata, based on particle swarm optimization(PSO), we have developed a method for seismic wave impedance inve...In order to improve the fine structure inversion ability of igneous rocks for the exploration of underlying strata, based on particle swarm optimization(PSO), we have developed a method for seismic wave impedance inversion. Through numerical simulation, we tested the effects of different algorithm parameters and different model parameterization methods on PSO wave impedance inversion, and analyzed the characteristics of PSO method. Under the conclusions drawn from numerical simulation, we propose the scheme of combining a cross-moving strategy based on a divided block model and high-frequency filtering technology for PSO inversion. By analyzing the inversion results of a wedge model of a pitchout coal seam and a coal coking model with igneous rock intrusion, we discuss the vertical and horizontal resolution, stability and reliability of PSO inversion. Based on the actual seismic and logging data from an igneous area, by taking a seismic profile through wells as an example, we discuss the characteristics of three inversion methods, including model-based wave impedance inversion, multi-attribute seismic inversion based on probabilistic neural network(PNN) and wave impedance inversion based on PSO.And we draw the conclusion that the inversion based on PSO method has a better result for this igneous area.展开更多
Non-liner wave equation inversion,wavelet analysis and artificial neural networks were used to obtain stratum parameters and the distribution of thin coal seams.The lithology of the water-bearing/resisting layer in th...Non-liner wave equation inversion,wavelet analysis and artificial neural networks were used to obtain stratum parameters and the distribution of thin coal seams.The lithology of the water-bearing/resisting layer in the Quaternary system was also predicted.The implementation process included calculating the well log parameters,stratum contrasting the seismic data and the well logs,and extracting,studying and predicting seismic attributes.Seismic inversion parameters,including the layer velocity and wave impedance,were calculated and effectively used for prediction and analysis.Prior knowledge and seismic interpretation were used to remedy a dearth of seismic data during the inversion procedure.This enhanced the stability of the inversion method.Non-linear seismic inversion and artificial neural networks were used to interpret coal seismic lithology and to study the water-bearing/resisting layer in the Quaternary system.Interpretation of the 1~2 m thin coal seams,and also of the water-bearing/resisting layer in the Quaternary system,is provided.The upper mining limit can be lifted from 60 m to 45 m.The predictions show that this method can provide reliable data useful for thin coal seam exploitation and for lifting the upper mining limit,which is one of the principles of green mining.展开更多
基金The authors acknowledge the sponsorship of National Natural Science Foundation of China(42174139,41974119,42030103)Laoshan Laboratory Science and Technology Innovation Program(LSKj202203406)Science Foundation from Innovation and Technology Support Program for Young Scientists in Colleges of Shandong Province and Ministry of Science and Technology of China(2019RA2136).
文摘Total organic carbon (TOC) prediction with elastic parameter inversions has been widely used in the identification and evaluation of source rocks. However, the elastic parameters used to predict TOC are not only determined by TOC but also depend on the other physical properties of source rocks. Besides, the TOC prediction with the elastic parameters inversion is an indirect method based on the statistical relationship obtained from well logs and experiment data. Therefore, we propose a rock physics model and define a TOC indicator mainly affected by TOC to predict TOC directly. The proposed rock physics model makes the equivalent elastic moduli of source rocks parameterized by the TOC indicator. Combining the equivalent elastic moduli of source rocks and Gray’s approximation leads to a novel linearized approximation of the P-wave reflection coefficient incorporating the TOC indicator. Model examples illustrate that the novel reflectivity approximation well agrees with the exact Zoeppritz equation until incident angles reach 40°. Convoluting the novel P-wave reflection approximation with seismic wavelets as the forward solver, an AVO inversion method based on the Bayesian theory is proposed to invert the TOC indicator with seismic data. The synthetic examples and field tests validate the feasibility and stability of the proposed AVO inversion approach. Using the inversion results of the TOC indicator, TOC is directly and accurately estimated in the target area.
基金supported by the Major Basic Research Development Program of China’s 973 Project(grant No.2007CB209608)the Science and Technology Innovation Foundation of CNPC(grant No.2010D-5006-0301)
文摘Although the ambiguity of seismic inversion is widely recognized in both theory and practice, so far as a concrete inversion example is concerned, there is not any objective, controllable method or any standard for how to evaluate and determine its ambiguity and reliability, especially for the high frequency components beyond the effective seismic frequency band. Taking log-constrained impedance inversion as an example, a new appraisal method is proposed on the basis of analyzing a simple geological model. Firstly, the inverted impedance model is transformed to a reflection coefficient series. Secondly, the maximum effective frequency of the real seismic data is chosen as a cutoff point and the reflection coefficient series is decomposed into two components by low-pass and high-pass filters. Thirdly, the geometrical reflection characteristics of the high-frequency components and that of the real seismic data are compared and analyzed. Then, the reliability of the inverted impedance model is appraised according to the similarity of geometrical characteristics between the high-frequency components and the real seismic data. The new method avoids some subjectivity in appraising the inverted result, and helps to enhance the reliability of reservoir prediction by impedance inversion technology.
文摘Comprehensive inversion of logging and seismic data presented in this paper is a method to improve seismic data resolution. It involves using ample high-frequency information and complete low-frequency information of known logging to make up for the lack of limited bandwidth of practical seismic recording, obtaining an approximate reflection coefficient sequence (or wave impedance) of high resolution by iterative inversion and providing more reliable seismic evidence for further lithologic interpretation and lateral tracking, correlation and prediction of thin reservoir. The comprehensive inversion can be realized in the following steps: (1) to establish an initial model of higher resolution; (2) to obtain wavelets, and (3) to constrain iterative inversion. The key to this inversion lies in building an initial model. It is assumed from our experience that when the initial model is properly given, iterative inversion can be quickly converged to the ideal result.
基金financially supported by the NSFC(Grant No.41974126 and 41674116)the National Key Research and Development Program of China(Grant No.2018YFA0702501)the 13th 5-Year Basic Research Program of China National Petroleum Corporation(CNPC)(2018A-3306)。
文摘Deep learning has achieved great success in a variety of research fields and industrial applications.However,when applied to seismic inversion,the shortage of labeled data severely influences the performance of deep learning-based methods.In order to tackle this problem,we propose a novel seismic impedance inversion method based on a cycle-consistent generative adversarial network(Cycle-GAN).The proposed Cycle-GAN model includes two generative subnets and two discriminative subnets.Three kinds of loss,including cycle-consistent loss,adversarial loss,and estimation loss,are adopted to guide the training process.Benefit from the proposed structure,the information contained in unlabeled data can be extracted,and adversarial learning further guarantees that the prediction results share similar distributions with the real data.Moreover,a neural network visualization method is adopted to show that the proposed CNN model can learn more distinguishable features than the conventional CNN model.The robustness experiments on synthetic data sets show that the proposed method can achieve better performances than other methods in most cases.And the blind-well experiments on real seismic profiles show that the predicted impedance curve of the proposed method maintains a better correlation with the true impedance curve.
基金the sponsorship of National Grand Project for Science and Technology(2016ZX05024004,2017ZX05009001,2017ZX05032003)the Fundamental Research Funds for the Central Universities(20CX06036A)+1 种基金the Postdoctoral Applied Research Project of Qingdao(QDYY20190040)the Science Foundation from SINOPEC Key Laboratory of Geophysics(wtyjy-wx2019-01-04)。
文摘Seismic amplitude variation with offset(AVO) inversion is an important approach for quantitative prediction of rock elasticity,lithology and fluid properties.With Biot-Gassmann's poroelasticity,an improved statistical AVO inversion approach is proposed.To distinguish the influence of rock porosity and pore fluid modulus on AVO reflection coefficients,the AVO equation of reflection coefficients parameterized by porosity,rock-matrix moduli,density and fluid modulus is initially derived from Gassmann equation and critical porosity model.From the analysis of the influences of model parameters on the proposed AVO equation,rock porosity has the greatest influences,followed by rock-matrix moduli and density,and fluid modulus has the least influences among these model parameters.Furthermore,a statistical AVO stepwise inversion method is implemented to the simultaneous estimation of rock porosity,rock-matrix modulus,density and fluid modulus.Besides,the Laplace probability model and differential evolution,Markov chain Monte Carlo algorithm is utilized for the stochastic simulation within Bayesian framework.Models and field data examples demonstrate that the simultaneous optimizations of multiple Markov chains can achieve the efficient simulation of the posterior probability density distribution of model parameters,which is helpful for the uncertainty analysis of the inversion and sets a theoretical fundament for reservoir characterization and fluid discrimination.
基金supported by the Key Foundation of Institute of Seismology,China Earthquake Administration( IS200916004)
文摘Yushu Ms7.1 earthquake occurred on the Ganzi-Yushu fault zone, across which we carried out a joint relative-gravity and seismic-reflection survey, and then performed a gravity inversion constrained by the seismic-reflection result. Based on the data of complete Bouguer gravity anomaly and seismic reflection, we obtained a layered interface structure in deep crust down to Moho. Our study showed that the inversion could reveal the interfaces of strata along the survey profile and the directions of regional faults in two-dimension. From the characteristics of the observed topography of the Moho basement, we tentatively confirmed that the uplift of eastern edge of Qinghai-Tibet plateau was caused by the subduetion of the Indian plate.
基金This research is supported by the Joint Funds of the National Natural Science Foundation of China(No.U20B2016)the National Natural Science Foundation of China(No.41874167)the National Natural Science Foundation of China(No.41904130).
文摘To improve the accuracy of inversion results,geological facies distributions are considered as additional constraints in the inversion process.However,the geological facies itself also has its own uncertainty.In this paper,the initial sedimentary facies maps are obtained by integrated geological analysis from well data,seismic attributes,and deterministic inversion results.Then the fi rst iteration of facies-constrained seismic inversion is performed.According to that result and other data such as geological information,the facies distribution can be updated using cluster analysis.The next round of facies-constrained inversion can then be performed.This process will be repeated until the facies inconsistency or error before and after the inversion is minimized.It forms a new iterative facies-constrained seismic inversion technique.Compared with conventional facies-constrained seismic inversion,the proposed method not only can reduces the non-uniqueness of seismic inversion results but also can improves its resolution.As a consequence,the sedimentary facies will be more consistent with the geology.A practical application demonstrated that the superposition relationship of sand bodies could be better delineated based on this new seismic inversion technique.The result highly increases the understanding of reservoir connectivity and its accuracy,which can be used to guide further development.
基金supported by National Key Basic Research Development Program (Grant No. 2007CB209600)National Major Science and Technology Program (Grant No. 2008ZX05010-002)
文摘We investigate the accuracy and robustness of moment tensor(MT)and stress inversion solutions derived from acoustic emissions(AEs)during the laboratory fracturing of prismatic Barre granite specimens.Pre-cut flaws in the specimens introduce a complex stress field,resulting in a spatial and temporal variation of focal mechanisms.Specifically,we consider two experimental setups:(1)where the rock is loaded in compression to generate primarily shear-type fractures and(2)where the material is loaded in indirect tension to generate predominantly tensile-type fractures.In each test,we first decompose AE moment tensors into double-couple(DC)and non-DC terms and then derive unambiguous normal and slip vectors using k-means clustering and an unstructured damped stress inversion algorithm.We explore temporal and spatial distributions of DC and non-DC events at different loading levels.The majority of the DC and the tensile non-DC events cluster around the pre-cut flaws,where macro-cracks later develop.Results of stress inversion are verified against the stress field from finite element(FE)modeling.A good agreement is found between the experimentally derived and numerically simulated stress orientations.To the best of the authors’knowledge,this work presents the first case where stress inversion methodologies are validated by numerical simulations at laboratory scale and under highly heterogeneous stress distributions.
文摘Using the technique of seismic moment tensor inversion, the source mechanisms of 10 earthquakes with Ms5.2that occurred in China from November 1996 to January 1998 were determined rapidly. The determined resultswere sent as 'Bulletins of Source Mechanism Parameters of Earthquakes' to the Seismic Regime Guards' Office,China Seismological Bureau, and the relevant provincial seismological bureaus. These bulletins have played rolein the fast response to large earthquakes.
基金supported by the National Natural Science Foundation of China under Grant No.42050104
文摘Deep learning is widely used for seismic impedance inversion,but few work provides in-depth research and analysis on designing the architectures of deep neural networks and choosing the network hyperparameters.This paper is dedicated to comprehensively studying on the significant aspects of deep neural networks that affect the inversion results.We experimentally reveal how network hyperparameters and architectures affect the inversion performance,and develop a series of methods which are proven to be effective in reconstructing high-frequency information in the estimated impedance model.Experiments demonstrate that the proposed multi-scale architecture is helpful to reconstruct more high-frequency details than a conventional network.Besides,the reconstruction of high-frequency information can be further promoted by introducing a perceptual loss and a generative adversarial network from the computer vision perspective.More importantly,the experimental results provide valuable references for designing proper network architectures in the seismic inversion problem.
基金This work was supported by the Special Fund of the Institute of Geophysics,China Earthquake Administration(Nos.DQJB22B19,DQJB22R29 and DQJB22B26)the National Natural Science Foundation of China(Nos.41974066,U1839209 and 42074053)。
文摘On September 16,2021,a MS6.0 earthquake struck Luxian County,one of the shale gas blocks in the Southeastern Sichuan Basin,China.To understand the seismogenic environment and its mechanism,we inverted a fine three-dimensional S-wave velocity model from ambient noise tomography using data from a newly deployed dense seismic array around the epicenter,by extracting and jointly inverting the Rayleigh phase and group velocities in the period of 1.6–7.2 s.The results showed that the velocity model varied significantly beneath different geological units.The Yujiasi syncline is characterized by low velocity at depths of~3.0–4.0 km,corresponding to the stable sedimentary layer in the Sichuan Basin.The eastern and western branches of the Huayingshan fault belt generally exhibit high velocities in the NE-SW direction,with a few local low-velocity zones.The Luxian MS6.0 earthquake epicenter is located at the boundary between the high-and low-velocity zones,and the earthquake sequences expand eastward from the epicenter at depths of 3.0–5.0 km.Integrated with the velocity variations around the epicenter,distribution of aftershock sequences,and focal mechanism solution,it is speculated that the seismogenic mechanism of the main shock might be interpreted as the reactivation of pre-existing faults by hydraulic fracturing.
基金National Nature Science Foundation of China (40334040) & Joint Seismological foundation of CEA (101026)
文摘A genetic algorithm of body waveform inversion is presented for better understanding of crustal and upper mantle structures with deep seismic sounding (DSS) waveform data. General reflection and transmission synthetic seismogram algorithm, which is capable of calculating the response of thin alternating high and low velocity layers, is applied as a solution for forward modeling, and the genetic algorithm is used to find the optimal solution of the inverse problem. Numerical tests suggest that the method has the capability of resolving low-velocity layers, thin alternating high and low velocity layers, and noise suppression. Waveform inversion using P-wave records from Zeku, Xiahe and Lintao shots in the seismic wide-angle reflection/refraction survey along northeastern Qinghai-Xizang (Tibeteau) Plateau has revealed fine structures of the bottom of the upper crust and alternating layers in the middle/lower crust and topmost upper mantle.
文摘Under the condition of thin interbeds with great lateral changes in terrestrial basins,a seismic meme inversion method is established based on the analysis of seismic sedimentology technology.The relationship between seismic waveform and high-frequency well logs is established through dynamic clustering of seismic waveform to improve the vertical and horizontal resolution of inversion results;meanwhile,by constructing the Bayesian inversion framework of different seismic facies,the real facies controlled inversion is realized.The forward model verification results show that the seismic meme inversion can realize precise prediction of 3 m thick thin interbeds,proving the rationality and high precision of the method.The application in the Daqing placanticline shows that the seismic meme inversion could identify 2 m thin interbeds,and the coincidence rates of inversion results and drilling data were more than 80%.The seismic meme inversion method can improve the accuracy of reservoir prediction and provides a useful mean for thin interbeds prediction in terrestrial basins.
基金Key Project (95-11-02-01) from China Seismological Bureau.Contribution No. RCEG200129, Research Center of Exploration Geophysi
文摘2-D velocity structure and tectonics of the crust and upper mantle is revealed by inversion of seismic refraction and wide-angle reflection traveltimes acquired along the profile L1 in the Changbaishan-Tianchi volcanic region. It is used in this study that seismic traveltime inversion for simultaneous determination of 2-D velocity and interface structure of the crust and upper mantle. The result shows that, under Changbaishan-Tianchi crater, there exists a low-velocity body in the shape of an inverted triangle, and the crustal reflecting boundaries and Moho all become lower by a varying margin of 2-6 km, forming a crustal root which is assumed to be the Changbaishan-Tianchi volcanic system. Finally, we make a comparison between our 2-D velocity model and the result from the studies by using trial-and-error forward modeling with SEIS83.
基金supported by the National Natural Science Foundation of China (No. 41174034)
文摘The firework algorithm(FWA) is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model show that the FWA performs well in both solution quality and efficiency. We apply the FWA in this study to crustal velocity structure inversion using regional seismic waveform data of central Gansu on the northeastern margin of the Qinghai-Tibet plateau. Seismograms recorded from the moment magnitude(MW) 5.4 Minxian earthquake enable obtaining an average crustal velocity model for this region. We initially carried out a series of FWA robustness tests in regional waveform inversion at the same earthquake and station positions across the study region,inverting two velocity structure models, with and without a low-velocity crustal layer; the accuracy of our average inversion results and their standard deviations reveal the advantages of the FWA for the inversion of regional seismic waveforms. We applied the FWA across our study area using three component waveform data recorded by nine broadband permanent seismic stations with epicentral distances ranging between 146 and 437 km. These inversion results show that the average thickness of the crust in this region is 46.75 km, while thicknesses of the sedimentary layer, and the upper, middle, and lower crust are 3.15,15.69, 13.08, and 14.83 km, respectively. Results also show that the P-wave velocities of these layers and the upper mantle are 4.47, 6.07, 6.12, 6.87, and 8.18 km/s,respectively.
基金part of an ongoing project of the National Important Industry Technological Development Project (High Precision 3D Seismic Technology of Coal Resources of Western China)the financial support from the National Basic Research Program of China (No.2009CB 219603)the National Key Scientific and Technological Project of China (No.2008ZX05035-005-003HZ)
文摘With the objective of establishing the necessary conditions for 3D seismic data from mountainous areas in western China, we compared the application results of wave impedance technology in the lithology and exploration of coal fields. First, we introduce principles and features of three kinds of inversion methods. i.e., Model-Based Inversion, Constrained Sparse Spike Inversion (CSSI) and Geology-Seismic Feature Inversion. Secondly, these inversion methods are contrasted in their application to 3D seismic data from some coalfields in western China. The main information provided by the research includes: improving the vertical resolution of coal deposit strata, inferring lateral variation of the lithology and predicting coal seams and their roof lithology. Finally, the comparison between the three methods shows that the model-based inversion has the higher resolution, while CSSI inversion has better waveform continuity. The geology-seismic feature inversion requires information from a large number of wells and many types of logging curves of good quality. All three methods can meet the requirements of seismic exploration for lithological exploration in coal fields.
基金This research is jointly funded by the National Key Research and Development Program of China(No.2017 YFC0307401)the National Natural Science Foundation of China(No.41230318)+1 种基金the Fundamental Research Funds for the Central Universities(No.201964017)and the National Science and Technology Major Project of China(No.2016ZX05024-001-002).
文摘In Recent years,seismic data have been widely used in seismic oceanography for the inversion of oceanic parameters represented by conductivity temperature depth(CTD).Using this technique,researchers can identify the water structure with high horizontal resolution,which compensates for the deficiencies of CTD data.However,conventional inversion methods are modeldriven,such as constrained sparse spike inversion(CSSI)and full waveform inversion(FWI),and typically require prior deterministic mapping operators.In this paper,we propose a novel inversion method based on a convolutional neural network(CNN),which is purely data-driven.To solve the problem of multiple solutions,we use stepwise regression to select the optimal attributes and their combination and take two-dimensional images of the selected attributes as input data.To prevent vanishing gradients,we use the rectified linear unit(ReLU)function as the activation function of the hidden layer.Moreover,the Adam and mini-batch algorithms are combined to improve stability and efficiency.The inversion results of field data indicate that the proposed method is a robust tool for accurately predicting oceanic parameters.
基金provided by the National Science and Technology Major Project(No.2011ZX05004-004)China National Petroleum Corporation Key Projects(No.2014E2105)the Priority Academic Program Development of Jiangsu Higher Education Institutions(PAPD)
文摘In order to improve the fine structure inversion ability of igneous rocks for the exploration of underlying strata, based on particle swarm optimization(PSO), we have developed a method for seismic wave impedance inversion. Through numerical simulation, we tested the effects of different algorithm parameters and different model parameterization methods on PSO wave impedance inversion, and analyzed the characteristics of PSO method. Under the conclusions drawn from numerical simulation, we propose the scheme of combining a cross-moving strategy based on a divided block model and high-frequency filtering technology for PSO inversion. By analyzing the inversion results of a wedge model of a pitchout coal seam and a coal coking model with igneous rock intrusion, we discuss the vertical and horizontal resolution, stability and reliability of PSO inversion. Based on the actual seismic and logging data from an igneous area, by taking a seismic profile through wells as an example, we discuss the characteristics of three inversion methods, including model-based wave impedance inversion, multi-attribute seismic inversion based on probabilistic neural network(PNN) and wave impedance inversion based on PSO.And we draw the conclusion that the inversion based on PSO method has a better result for this igneous area.
基金Projects 40574057 and 40874054 supported by the National Natural Science Foundation of ChinaProjects 2007CB209400 by the National Basic Research Program of ChinaFoundation of China University of Mining and Technology (OF4471)
文摘Non-liner wave equation inversion,wavelet analysis and artificial neural networks were used to obtain stratum parameters and the distribution of thin coal seams.The lithology of the water-bearing/resisting layer in the Quaternary system was also predicted.The implementation process included calculating the well log parameters,stratum contrasting the seismic data and the well logs,and extracting,studying and predicting seismic attributes.Seismic inversion parameters,including the layer velocity and wave impedance,were calculated and effectively used for prediction and analysis.Prior knowledge and seismic interpretation were used to remedy a dearth of seismic data during the inversion procedure.This enhanced the stability of the inversion method.Non-linear seismic inversion and artificial neural networks were used to interpret coal seismic lithology and to study the water-bearing/resisting layer in the Quaternary system.Interpretation of the 1~2 m thin coal seams,and also of the water-bearing/resisting layer in the Quaternary system,is provided.The upper mining limit can be lifted from 60 m to 45 m.The predictions show that this method can provide reliable data useful for thin coal seam exploitation and for lifting the upper mining limit,which is one of the principles of green mining.