The fi rst arrival waveform inversion(FAWI)has a strong nonlinearity due to the objective function using L2 parametrization.When the initial velocity is not accurate,the inversion can easily fall into local minima.In ...The fi rst arrival waveform inversion(FAWI)has a strong nonlinearity due to the objective function using L2 parametrization.When the initial velocity is not accurate,the inversion can easily fall into local minima.In the full waveform inversion method,adding a cross-correlation function to the objective function can eff ectively reduce the nonlinearity of the inversion process.In this paper,the nonlinearity of this process is reduced by introducing the correlation objective function into the FAWI and by deriving the corresponding gradient formula.We then combine the first-arrival wave travel-time tomography with the FAWI to form a set of inversion processes.This paper uses the limited memory Broyden-Fletcher-Goldfarb-Shanno(L-BFGS)algorithm to improve the computational effi ciency of inversion and solve the problem of the low effi ciency of the FAWI method.The overthrust model and fi eld data test show that the method used in this paper can eff ectively reduce the nonlinearity of inversion and improve the inversion calculation effi ciency at the same time.展开更多
The objective function of full waveform inversion is a strong nonlinear function,the inversion process is not unique,and it is easy to fall into local minimum.Firstly,in the process of wavefield reconstruction,the wav...The objective function of full waveform inversion is a strong nonlinear function,the inversion process is not unique,and it is easy to fall into local minimum.Firstly,in the process of wavefield reconstruction,the wave equation is introduced into the construction of objective function as a penalty term to broaden the search space of solution and reduce the risk of falling into local minimum.In addition,there is no need to calculate the adjoint wavefield in the inversion process,which can significantly improve the calculation efficiency;Secondly,considering that the total variation constraint can effectively reconstruct the discontinuous interface in the velocity model,this paper introduces the weak total variation constraint to avoid the excessive smooth estimation of the model under the strong total variation constraint.The disadvantage of this strategy is that it is highly dependent on the initial model.In view of this,this paper takes the long wavelength initial model obtained by first arrival traveltime tomography as a prior model constraint,and proposes a weak total variation constrained wavefield reconstruction inversion method based on first arrival traveltime tomography.Numerical experimental results show that the new method reduces the dependence on the initial model,the interface description is more accurate,the error is reduced,and the iterative convergence efficiency is significantly improved.展开更多
Wave front healing is a common natural phenomenon.To further investigate wave front healing,we simulated wave propagation in a spherical anomaly surrounded by homogeneous media using a high-order finite difference sol...Wave front healing is a common natural phenomenon.To further investigate wave front healing,we simulated wave propagation in a spherical anomaly surrounded by homogeneous media using a high-order finite difference solution of the acoustic equation.Furthermore,we analyzed the characteristics of the wave propagation in the anomaly,and found that they are related to the dominant frequency of the seismic wave and the dimensions of the anomaly.Through quantitative comparison of the wave front energy of the diffracted wave and transmitted wave,we summarized the influences of the wave front healing on seismic tomography.We conclude that,under the strong scattering condition,only positive anomalies can be inverted by ray-based tomography,only large anomalies can be inverted by finite-frequency tomography,and small negative anomalies cannot be inverted by any first-arrival traveltime tomographic methods.These conclusions are verified by tomographic experiments based on different theoretical models.Finally,we propose that more information besides the first-arrival traveltime should be used to invert the high wave number components of the media.Besides the above acquisitions of wave front healing on seismic tomography,we explain the banana-doughnut phenomena,and offer a new insight into the wave scattering,which should be important for better understanding the wave propagation and seismic inversion.展开更多
基金supported by the Major Scientific and Technological Project of PetroChina (ZD2019-183-003)Project of National Natural Science Foundation of China (42074133)+1 种基金the Fundamental Research Funds for the Central Universities (19CX02056A)Project of State Key Laboratory of Shale Oil and Gas Enrichment Mechanisms and Effective Development (33550000-21-FW0399-0009)
文摘The fi rst arrival waveform inversion(FAWI)has a strong nonlinearity due to the objective function using L2 parametrization.When the initial velocity is not accurate,the inversion can easily fall into local minima.In the full waveform inversion method,adding a cross-correlation function to the objective function can eff ectively reduce the nonlinearity of the inversion process.In this paper,the nonlinearity of this process is reduced by introducing the correlation objective function into the FAWI and by deriving the corresponding gradient formula.We then combine the first-arrival wave travel-time tomography with the FAWI to form a set of inversion processes.This paper uses the limited memory Broyden-Fletcher-Goldfarb-Shanno(L-BFGS)algorithm to improve the computational effi ciency of inversion and solve the problem of the low effi ciency of the FAWI method.The overthrust model and fi eld data test show that the method used in this paper can eff ectively reduce the nonlinearity of inversion and improve the inversion calculation effi ciency at the same time.
基金supported by National Key R&D Program of China under contract number 2019YFC0605503CThe Major projects of CNPC under contract number(ZD2019-183-003)+2 种基金the Major projects during the 14th Five-year Plan period under contract number 2021QNLM020001the National Outstanding Youth Science Foundation under contract number 41922028the Funds for Creative Research Groups of China under contract number 41821002.
文摘The objective function of full waveform inversion is a strong nonlinear function,the inversion process is not unique,and it is easy to fall into local minimum.Firstly,in the process of wavefield reconstruction,the wave equation is introduced into the construction of objective function as a penalty term to broaden the search space of solution and reduce the risk of falling into local minimum.In addition,there is no need to calculate the adjoint wavefield in the inversion process,which can significantly improve the calculation efficiency;Secondly,considering that the total variation constraint can effectively reconstruct the discontinuous interface in the velocity model,this paper introduces the weak total variation constraint to avoid the excessive smooth estimation of the model under the strong total variation constraint.The disadvantage of this strategy is that it is highly dependent on the initial model.In view of this,this paper takes the long wavelength initial model obtained by first arrival traveltime tomography as a prior model constraint,and proposes a weak total variation constrained wavefield reconstruction inversion method based on first arrival traveltime tomography.Numerical experimental results show that the new method reduces the dependence on the initial model,the interface description is more accurate,the error is reduced,and the iterative convergence efficiency is significantly improved.
基金supported by National Natural Science Foundation of China(Grant No. 40804023)National Basic Research Program of China (Grant No. 2006CB202402)+1 种基金Hi-tech R&D Program of China (Grant No.2008AA093001)Project of State Key Laboratory of Marine Geology of China (Grant No. MG200909)
文摘Wave front healing is a common natural phenomenon.To further investigate wave front healing,we simulated wave propagation in a spherical anomaly surrounded by homogeneous media using a high-order finite difference solution of the acoustic equation.Furthermore,we analyzed the characteristics of the wave propagation in the anomaly,and found that they are related to the dominant frequency of the seismic wave and the dimensions of the anomaly.Through quantitative comparison of the wave front energy of the diffracted wave and transmitted wave,we summarized the influences of the wave front healing on seismic tomography.We conclude that,under the strong scattering condition,only positive anomalies can be inverted by ray-based tomography,only large anomalies can be inverted by finite-frequency tomography,and small negative anomalies cannot be inverted by any first-arrival traveltime tomographic methods.These conclusions are verified by tomographic experiments based on different theoretical models.Finally,we propose that more information besides the first-arrival traveltime should be used to invert the high wave number components of the media.Besides the above acquisitions of wave front healing on seismic tomography,we explain the banana-doughnut phenomena,and offer a new insight into the wave scattering,which should be important for better understanding the wave propagation and seismic inversion.
