The invention discloses a time migration method for a three-dimensional pre-stack
wave equation in an imaging
ray domain, which relates to the technical field of
exploration geophysics and comprises the following steps of: deducing a bidirectional
acoustic wave equation in an imaging
ray coordinate
system based on a Riemann
wave field extrapolation theory in a three-dimensional semi-orthogonal coordinate
system; an equation is directly solved by adopting a
finite difference method, extrapolation of a seismic source and a receiving end
wave field is achieved, and errors caused by Green
function approximation in a traditional migration method are avoided; and generating a pre-stack time migration result by applying a zero-
lag cross-correlation
imaging condition. Through first-order linear disturbance approximation, an accurate three-dimensional time-
depth conversion partial differential equation set is simplified into an optimized layer stripping
algorithm, and a corresponding equation suitable for a weak
transverse velocity change medium is deduced. According to the invention, the working process of three-dimensional
wave equation time migration is realized in an imaging
ray coordinate
system, and an accurate seismic imaging result is generated; two-dimensional plane
hypothesis limitation is avoided, and high-precision imaging technical guarantee is provided for seismic
data processing of a complex structure area.