Abstract
This paper seeks to reconstruct the parameters of elastic layered media such as P-wave velocity, S-wave velocity, density and thickness from the multioffset seismic reflection data. Since the data are highly non-linear to the low-wavenumber components, the non-linear waveform inversion method, with the aid of generalized ray theory, is proposed to solve this problem in space–time domain. As opposed to the layer-stripping method, the present method attempts to invert all layer parameters simultaneously, thus reducing the cumulative errors resulting from the upper layers. The parameters are inverted by minimizing the weighted square error between the observed data and the calculated data of the layered model, the optimization of which is based on the quasi-Newton method. In synthetic tests, we find that the inverted results are good when the variation of parameters between layers is not too large. The modified method for large variation of parameters is first to fix those of deeper layers and neglect the signals reflected from them, then recover some other parameters simultaneously until those of upper layers attain a stable value, and finally, invert all parameters simultaneously. The results so obtained show a significant improvement. This method was tested to be stable in the presence of noise in seismograms.
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
More From: International Journal for Numerical and Analytical Methods in Geomechanics
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.