
Inversion of seismic data: accuracy and convergence of an iterative scheme based on acoustic imaging
Author(s) -
Carrion Philip M.,
Kuo John T.,
Patton Waldo A.
Publication year - 1984
Publication title -
geophysical journal of the royal astronomical society
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0016-8009
DOI - 10.1111/j.1365-246x.1984.tb02233.x
Subject(s) - slowness , inversion (geology) , seismogram , geology , synthetic seismogram , synthetic data , wave equation , seismic wave , mathematical analysis , plane wave , acoustic wave equation , continuation , geodesy , mathematics , acoustic wave , algorithm , seismology , acoustics , physics , computer science , optics , tectonics , programming language
Summary This paper describes the accuracy, stability and convergence of a new iterative inversion scheme based on acoustic imaging. The scheme first decomposes the observed data (seismogram) into plane waves (PWD); then it continues the PWD data downwards. The downward continuation can take place in the time or the frequency domains. The data are continued by finding the solution of the acoustic wave equation for a trial velocity selected from a class of admissible functions. Next the downward continued data are projected on to the plane of horizontal slowness p versus depth z . The imaging principle used to recover the velocity is based on finding the maximum of the projected wavefield energy, which corresponds to the critically reflected and refracted arrivals. This maximum must coincide with the true velocity profile in the p‐z domain. If not, the trial velocity function is adjusted to be closer to the maximum of the transformed data and the procedure is restarted. Examples with synthetic and real data are presented.