| Abstract: | We propose a wave scattering approach to the problem of deconvolution by the inversion of the reflection seismogram. Rather than using the least-squares approach, we study the full wave solution of the one-dimensional wave equation for deconvolution. Randomness of the reflectivity is not a necessary assumption in this method. Both the reflectivity and the section multiple train can be predicted from the boundary data (the reflection seismogram). This is in contrast to the usual statistical approach in which reflectivity is unpredictable and random, and the section multiple train is the only predictable component of the seismogram. The proposed scattering approach also differs from Claerbout's method based on the Kunetz equation. The coupled first-order hyperbolic wave equations have been obtained from the equation of motion and the law of elasticity. These equations have been transformed in terms of characteristics. A finite-difference numerical scheme for the downward continuation of the free-surface reflection seismogram has been developed. The discrete causal solutions for forward and inverse problems have been obtained. The computer algorithm recursively solves for the pressure and particle velocity response and the impedance log. The method accomplishes deconvolution and impedance log reconstruction. We have tested the method by computer model experiments and obtained satisfactory results using noise-free synthetic data. Further study is recommended for the method's application to real data. |
