Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based ...Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based on the optimal steep descent methods, we present an algorithm which combines the preconditioned bi-conjugated gradient stable method and the multi-grid method to compute the wave propagation and the gradient space. The multiple scale prosperity of the waveform inversion and the multi-grid method can overcome the inverse problems local minima defect and accelerate convergence. The local inhomogeneous three-hole model simulated results and the Marmousi model certify the algorithm effectiveness.展开更多
The complexity of an elastic wavefield increases the nonlinearity of inversion, To some extent, multiscale inversion decreases the nonlinearity of inversion and prevents it from falling into local extremes. A multisca...The complexity of an elastic wavefield increases the nonlinearity of inversion, To some extent, multiscale inversion decreases the nonlinearity of inversion and prevents it from falling into local extremes. A multiscale strategy based on the simultaneous use of frequency groups and layer stripping method based on damped wave field improves the stability of inversion. A dual-level parallel algorithm is then used to decrease the computational cost and improve practicability. The seismic wave modeling of a single frequency and inversion in a frequency group are computed in parallel by multiple nodes based on multifrontal massively parallel sparse direct solver and MPI. Numerical tests using an overthrust model show that the proposed inversion algorithm can effectively improve the stability and accuracy of inversion by selecting the appropriate inversion frequency and damping factor in low- frequency seismic data.展开更多
Amplitude variations with offset or incident angle （AVO/AVA） inversion are typically combined with statistical methods, such as Bayesian inference or deterministic inversion. We propose a joint elastic inversion met...Amplitude variations with offset or incident angle （AVO/AVA） inversion are typically combined with statistical methods, such as Bayesian inference or deterministic inversion. We propose a joint elastic inversion method in the time and frequency domain based on Bayesian inversion theory to improve the resolution of the estimated P- and S-wave velocities and density. We initially construct the objective function using Bayesian inference by combining seismic data in the time and frequency domain. We use Cauchy and Gaussian probability distribution density functions to obtain the prior information for the model parameters and the likelihood function, respectively. We estimate the elastic parameters by solving the initial objective function with added model constraints to improve the inversion robustness. The results of the synthetic data suggest that the frequency spectra of the estimated parameters are wider than those obtained with conventional elastic inversion in the time domain. In addition, the proposed inversion approach offers stronger antinoising compared to the inversion approach in the frequency domain. Furthermore, results from synthetic examples with added Gaussian noise demonstrate the robustness of the proposed approach. From the real data, we infer that more model parameter details can be reproduced with the proposed joint elastic inversion.展开更多
The regularization contributes to the resolution and stability in geophysical inversion. The authors apply dual-parameter shaping regularization to full waveform inversion, aiming at two points ： （ 1 ） improving th...The regularization contributes to the resolution and stability in geophysical inversion. The authors apply dual-parameter shaping regularization to full waveform inversion, aiming at two points ： （ 1 ） improving the boundary resolution, and （2） increasing convergence. Firstly, the forward modeling is done, and the inversion is processed with the optimal solution. Compared with classical Tikhonov regularization scheme, the method re fleets better resolution and stronger convergence. Then, Marmousi model is experimented and inversed, and the deep structure has a sharper outline. The phase residual comparison illustrates weaker cycle-slipping. And a choice scheme of parameter is applied in FWI.展开更多
Full waveform inversion( FWI) is a high resolution inversion method,which can reveal detailed information of the structure and lithology under complex geological background. It is limited by many kinds of noises when ...Full waveform inversion( FWI) is a high resolution inversion method,which can reveal detailed information of the structure and lithology under complex geological background. It is limited by many kinds of noises when the method applied to the real seismic data. Based on Huber function criterion,the objective function combinates the anti-noise of L1 norm and the stability of L2 norm in theory,the authors derive the gradient formula of the Huber function by using L-BFGS algorithm for FWI. The new method is proved by synthetic seismic data with the Gaussian noise and the impulse noise. Numerical test results show that L-BFGS algorithm is applied to the frequency domain FWI with the convergence speed and high calculation accuracy,and can effectively reduce computer memory usage; and the Huber function is more robust and stable than L2 norm even with the noises.展开更多
基金supported by the China State Key Science and Technology Project on Marine Carbonate Reservoir Characterization (No. 2011ZX05004-003)the Basic Research Programs of CNPC during the 12th Five-Year Plan Period (NO.2011A-3603)+1 种基金the Natural Science Foundation of China (No.41104066)the RIPED Young Professional Innovation Fund (NO.2010-13-16-02, 2010-A-26-02)
文摘Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based on the optimal steep descent methods, we present an algorithm which combines the preconditioned bi-conjugated gradient stable method and the multi-grid method to compute the wave propagation and the gradient space. The multiple scale prosperity of the waveform inversion and the multi-grid method can overcome the inverse problems local minima defect and accelerate convergence. The local inhomogeneous three-hole model simulated results and the Marmousi model certify the algorithm effectiveness.
基金supported by the Natural Science Foundation of China(No.41374122)
文摘The complexity of an elastic wavefield increases the nonlinearity of inversion, To some extent, multiscale inversion decreases the nonlinearity of inversion and prevents it from falling into local extremes. A multiscale strategy based on the simultaneous use of frequency groups and layer stripping method based on damped wave field improves the stability of inversion. A dual-level parallel algorithm is then used to decrease the computational cost and improve practicability. The seismic wave modeling of a single frequency and inversion in a frequency group are computed in parallel by multiple nodes based on multifrontal massively parallel sparse direct solver and MPI. Numerical tests using an overthrust model show that the proposed inversion algorithm can effectively improve the stability and accuracy of inversion by selecting the appropriate inversion frequency and damping factor in low- frequency seismic data.
基金supported by the National Nature Science Foundation Project(Nos.41604101 and U1562215)the National Grand Project for Science and Technology(No.2016ZX05024-004)+2 种基金the Natural Science Foundation of Shandong(No.BS2014NJ005)Science Foundation from SINOPEC Key Laboratory of Geophysics(No.33550006-15-FW2099-0027)the Fundamental Research Funds for the Central Universities
文摘Amplitude variations with offset or incident angle （AVO/AVA） inversion are typically combined with statistical methods, such as Bayesian inference or deterministic inversion. We propose a joint elastic inversion method in the time and frequency domain based on Bayesian inversion theory to improve the resolution of the estimated P- and S-wave velocities and density. We initially construct the objective function using Bayesian inference by combining seismic data in the time and frequency domain. We use Cauchy and Gaussian probability distribution density functions to obtain the prior information for the model parameters and the likelihood function, respectively. We estimate the elastic parameters by solving the initial objective function with added model constraints to improve the inversion robustness. The results of the synthetic data suggest that the frequency spectra of the estimated parameters are wider than those obtained with conventional elastic inversion in the time domain. In addition, the proposed inversion approach offers stronger antinoising compared to the inversion approach in the frequency domain. Furthermore, results from synthetic examples with added Gaussian noise demonstrate the robustness of the proposed approach. From the real data, we infer that more model parameter details can be reproduced with the proposed joint elastic inversion.
文摘The regularization contributes to the resolution and stability in geophysical inversion. The authors apply dual-parameter shaping regularization to full waveform inversion, aiming at two points ： （ 1 ） improving the boundary resolution, and （2） increasing convergence. Firstly, the forward modeling is done, and the inversion is processed with the optimal solution. Compared with classical Tikhonov regularization scheme, the method re fleets better resolution and stronger convergence. Then, Marmousi model is experimented and inversed, and the deep structure has a sharper outline. The phase residual comparison illustrates weaker cycle-slipping. And a choice scheme of parameter is applied in FWI.
基金Supported by the National "863" Project(No.2014AA06A605)
文摘Full waveform inversion( FWI) is a high resolution inversion method,which can reveal detailed information of the structure and lithology under complex geological background. It is limited by many kinds of noises when the method applied to the real seismic data. Based on Huber function criterion,the objective function combinates the anti-noise of L1 norm and the stability of L2 norm in theory,the authors derive the gradient formula of the Huber function by using L-BFGS algorithm for FWI. The new method is proved by synthetic seismic data with the Gaussian noise and the impulse noise. Numerical test results show that L-BFGS algorithm is applied to the frequency domain FWI with the convergence speed and high calculation accuracy,and can effectively reduce computer memory usage; and the Huber function is more robust and stable than L2 norm even with the noises.
