The parameter reconstruction of strong-scattering media is a challenge for conventional full waveform inversion(FWI).Direct envelope inversion(DEI)is an effective method for large-scale and strongscattering structures...The parameter reconstruction of strong-scattering media is a challenge for conventional full waveform inversion(FWI).Direct envelope inversion(DEI)is an effective method for large-scale and strongscattering structures imaging without the need of low-frequency seismic data.However,the current DEI methods are all based on the acoustic approximation.Whereas,in real cases,seismic records are the combined effects of the subsurface multi-parameters.Therefore,the study of DEI in elastic media is necessary for the accurate inversion of strong-scattering structures,such as salt domes.In this paper,we propose an elastic direct envelope inversion(EDEI)method based on wave mode decomposition.We define the objective function of EDEI using multi-component seismic data and derive its gradient formulation.To reduce the coupling effects of multi-parameters,we introduce the wave mode decomposition method into the gradient calculation of EDEI.The update of Vp is primarily the contributions of decomposed P-waves.Two approaches on Vs gradient calculation are proposed,i.e.using the petrophysical relation and wave mode decomposition method.Finally,we test the proposed method on a layered salt model and the SEG/EAGE salt model.The results show that the proposed EDEI method can reconstruct reliable large-scale Vp and Vs models of strong-scattering salt structures.The successive elastic FWI can obtain high-precision inversion results of the strong-scattering salt model.The proposed method also has a good anti-noise performance in the moderate noise level.展开更多
The multi-parameter inverse scattering problem of elastic waveequation with single fre- quency is investigated within Bornapproximation. By use of a wideband measuring scheme in which bothtransmitters and receivers sc...The multi-parameter inverse scattering problem of elastic waveequation with single fre- quency is investigated within Bornapproximation. By use of a wideband measuring scheme in which bothtransmitters and receivers scan over the half-space surface, theformula of the scattering field of elastic wave is derived. Fourtypes of mode conversion of elastic wave(P→P, P→S, S→P, S→S)areseparated from the scattering field. These components containsufficient information for usto recon- struct the configuration ofthe density and Lame parameters of the medium.展开更多
In this paper we calculate a synthetic medium surface displacement response that is consistent with real measurement data by applying the least-square principle and a niche genetic algorithm to the parameters inversio...In this paper we calculate a synthetic medium surface displacement response that is consistent with real measurement data by applying the least-square principle and a niche genetic algorithm to the parameters inversion problem of the wave equation in a two-phase medium. We propose a niche genetic multi-parameter （including porosity, solid phase density and fluid phase density） joint inversion algorithm based on a two-phase fractured medium in the BISQ model. We take the two-phase fractured medium of the BISQ model in a two- dimensional half space as an example, and carry out the numerical reservoir parameters inversion. Results show that this method is very convenient for solving the parameters inversion problem for the wave equation in a two-phase medium, and has the advantage of strong noise rejection. Relative to conventional genetic algorithms, the niche genetic algorithm based on a sharing function can not only significantly speed up the convergence, but also improve the inversion precision.展开更多
The waveform inversion method is applied-- using synthetic ocean-bottom seismometer （OBS） data--to study oceanic crust structure. A niching genetic algorithm （NGA） is used to implement the inversion for the thickn...The waveform inversion method is applied-- using synthetic ocean-bottom seismometer （OBS） data--to study oceanic crust structure. A niching genetic algorithm （NGA） is used to implement the inversion for the thickness and P-wave velocity of each layer, and to update the model by minimizing the objective function, which consists of the misfit and cross-correlation of observed and synthetic waveforms. The influence of specific NGA method parameters is discussed, and suitable values are presented. The NGA method works well for various observation systems, such as those with irregular and sparse distribu- tion of receivers as well as single receiver systems. A strategy is proposed to accelerate the convergence rate by a factor of five with no increase in computational complex- ity; this is achieved using a first inversion with several generations to impose a restriction on the preset range of each parameter and then conducting a second inversion with the new range. Despite the successes of this method, its usage is limited. A shallow water layer is not favored because the direct wave in water will suppress the useful reflection signals from the crust. A more precise calculation of the air-gun source signal should be considered in order to better simulate waveforms generated in realistic situa- tions; further studies are required to investigate this issue.展开更多
Using a bottom simulating reflector(BSR)on a seismic profile to identify marine gas hydrate is a traditional seismic exploration method.However,owing to the abundance differences between the gas hydrate and free gas i...Using a bottom simulating reflector(BSR)on a seismic profile to identify marine gas hydrate is a traditional seismic exploration method.However,owing to the abundance differences between the gas hydrate and free gas in different regions,the BSR may be unremarkable on the seismic profile and invisible in certain cases.With the improvement of exploration precision,difficulty arises in meeting the requirements of distinguishing the abundance differences in the gas hydrate based on BSR.Hence,we studied other sensitive attributes to ascertain the existence of gas hydrate and its abundance variations,eventually improving the success rate of drilling and productivity.In this paper,we analyzed the contradiction between the seismic profile data and drilling sampling data from the Blake Ridge.We extracted different attributes and performed multi-parameter constraint analysis based on the prestack elastic wave impedance inversion.Then,we compared the analysis results with the drilling sampling data.Eventually,we determined five sensitive attributes that can better indicate the existence of gas hydrate and its abundance variations.This method overcomes the limitations of recognizing the gas hydrate methods based on BSR or single inversion attribute.Moreover,the conclusions can notably improve the identification accuracy of marine gas hydrate and provide excellent reference significance for the recognition of marine gas hydrate.Notably,the different geological features of reservoirs feature different sensitivities to the prestacking attributes when using the prestack elastic inversion in different areas.展开更多
基金financial support jointly provided by the National Key R&D Program of China under contract number 2019YFC0605503Cthe Major Projects during the 14th Five-year Plan period under contract number 2021QNLM020001+2 种基金the National Outstanding Youth Science Foundation under contract number 41922028the Funds for Creative Research Groups of China under contract number 41821002the Major Projects of CNPC under contract number ZD2019-183-003。
文摘The parameter reconstruction of strong-scattering media is a challenge for conventional full waveform inversion(FWI).Direct envelope inversion(DEI)is an effective method for large-scale and strongscattering structures imaging without the need of low-frequency seismic data.However,the current DEI methods are all based on the acoustic approximation.Whereas,in real cases,seismic records are the combined effects of the subsurface multi-parameters.Therefore,the study of DEI in elastic media is necessary for the accurate inversion of strong-scattering structures,such as salt domes.In this paper,we propose an elastic direct envelope inversion(EDEI)method based on wave mode decomposition.We define the objective function of EDEI using multi-component seismic data and derive its gradient formulation.To reduce the coupling effects of multi-parameters,we introduce the wave mode decomposition method into the gradient calculation of EDEI.The update of Vp is primarily the contributions of decomposed P-waves.Two approaches on Vs gradient calculation are proposed,i.e.using the petrophysical relation and wave mode decomposition method.Finally,we test the proposed method on a layered salt model and the SEG/EAGE salt model.The results show that the proposed EDEI method can reconstruct reliable large-scale Vp and Vs models of strong-scattering salt structures.The successive elastic FWI can obtain high-precision inversion results of the strong-scattering salt model.The proposed method also has a good anti-noise performance in the moderate noise level.
基金Foundation of Ph.D Program of the State Education Commission of China
文摘The multi-parameter inverse scattering problem of elastic waveequation with single fre- quency is investigated within Bornapproximation. By use of a wideband measuring scheme in which bothtransmitters and receivers scan over the half-space surface, theformula of the scattering field of elastic wave is derived. Fourtypes of mode conversion of elastic wave(P→P, P→S, S→P, S→S)areseparated from the scattering field. These components containsufficient information for usto recon- struct the configuration ofthe density and Lame parameters of the medium.
基金sponsored by the National Science and Technology Major Project(Grant No.2011ZX05025-001-07)
文摘In this paper we calculate a synthetic medium surface displacement response that is consistent with real measurement data by applying the least-square principle and a niche genetic algorithm to the parameters inversion problem of the wave equation in a two-phase medium. We propose a niche genetic multi-parameter （including porosity, solid phase density and fluid phase density） joint inversion algorithm based on a two-phase fractured medium in the BISQ model. We take the two-phase fractured medium of the BISQ model in a two- dimensional half space as an example, and carry out the numerical reservoir parameters inversion. Results show that this method is very convenient for solving the parameters inversion problem for the wave equation in a two-phase medium, and has the advantage of strong noise rejection. Relative to conventional genetic algorithms, the niche genetic algorithm based on a sharing function can not only significantly speed up the convergence, but also improve the inversion precision.
基金supported by the National Natural Science Foundation grant No.41174034the Major State Basic Research Development Program of China(973 Program)
文摘The waveform inversion method is applied-- using synthetic ocean-bottom seismometer （OBS） data--to study oceanic crust structure. A niching genetic algorithm （NGA） is used to implement the inversion for the thickness and P-wave velocity of each layer, and to update the model by minimizing the objective function, which consists of the misfit and cross-correlation of observed and synthetic waveforms. The influence of specific NGA method parameters is discussed, and suitable values are presented. The NGA method works well for various observation systems, such as those with irregular and sparse distribu- tion of receivers as well as single receiver systems. A strategy is proposed to accelerate the convergence rate by a factor of five with no increase in computational complex- ity; this is achieved using a first inversion with several generations to impose a restriction on the preset range of each parameter and then conducting a second inversion with the new range. Despite the successes of this method, its usage is limited. A shallow water layer is not favored because the direct wave in water will suppress the useful reflection signals from the crust. A more precise calculation of the air-gun source signal should be considered in order to better simulate waveforms generated in realistic situa- tions; further studies are required to investigate this issue.
基金supported by the National Natural Science Foundation of China (No. 41230318)
文摘Using a bottom simulating reflector(BSR)on a seismic profile to identify marine gas hydrate is a traditional seismic exploration method.However,owing to the abundance differences between the gas hydrate and free gas in different regions,the BSR may be unremarkable on the seismic profile and invisible in certain cases.With the improvement of exploration precision,difficulty arises in meeting the requirements of distinguishing the abundance differences in the gas hydrate based on BSR.Hence,we studied other sensitive attributes to ascertain the existence of gas hydrate and its abundance variations,eventually improving the success rate of drilling and productivity.In this paper,we analyzed the contradiction between the seismic profile data and drilling sampling data from the Blake Ridge.We extracted different attributes and performed multi-parameter constraint analysis based on the prestack elastic wave impedance inversion.Then,we compared the analysis results with the drilling sampling data.Eventually,we determined five sensitive attributes that can better indicate the existence of gas hydrate and its abundance variations.This method overcomes the limitations of recognizing the gas hydrate methods based on BSR or single inversion attribute.Moreover,the conclusions can notably improve the identification accuracy of marine gas hydrate and provide excellent reference significance for the recognition of marine gas hydrate.Notably,the different geological features of reservoirs feature different sensitivities to the prestacking attributes when using the prestack elastic inversion in different areas.
