Method and apparatus for migration by finite differences
AssigneePhillips Petroleum Company
InventorBrackin A. Smith
Life of this patent
Estimating the term…
Abstract
An improved method for eliminating reflections from artificial model boundaries encountered in finite-difference acoustic wave propagation computations is based on applying two different boundary conditions at the artificial boundary. In this method one boundary condition generates reflections having the same polarity as the incident acoustic wave, and the other boundary condition generates reflections with opposite polarity. To apply the improved method simply compute both boundary conditions on a small narrow region adjacent the artificial boundary on a separate set of grids and average the two solutions.
Claims14
A computer implemented method for canceling reflections from an artificial boundary when computing propagation of an advancing acoustic wave over a large scale earth velocity model limited by said artificial boundary, wherein space in said velocity model is divided into a large number of elementary grid blocks, said method comprising:
- (a) defining a staggered grid region within said large scale velocity model for computing propagation of a plurality of wavefields, wherein said plurality of wavefields includes a reflected scalar pressure wavefield, and a reflected vector particle velocity wavefield;
- (b) wherein said staggered grid region is located adjacent to said artificial boundary, and is much smaller compared to the interior of said velocity model;
- (c) applying a first boundary condition of pressure equal zero at said artificial boundary for solution of an acoustic wave equation, thereby generating said reflected scalar pressure wavefield in-phase with said advancing acoustic wave and said reflected vector particle velocity wavefield out-of-phase with said advancing acoustic wave;
- (d) applying a second boundary condition of particle velocity equal zero at said artificial boundary for solution of said acoustic wave equation, thereby generating said reflected scalar pressure wavefield out-of-phase with said advancing acoustic wave and said reflected vector particle velocity wavefield in-phase with said advancing acoustic wave;
- and
- (e) averaging said reflected pressure wavefield, and averaging said vector particle velocity wavefield generated in steps (c) and (d) over said staggered grid region to essentially cancel reflections from said artificial boundary.
Apparatus for canceling reflections from an artificial boundary when computing propagation of an advancing acoustic wave over a large scale earth velocity model, wherein space in said velocity model is divided into a large number of elementary grid blocks, said apparatus comprising:
- a computer programmed for:
- i) storing said large scale earth velocity model in a computer memory, wherein said earth velocity model includes a staggered grid region, and wherein said staggered grid region is effective for computing propagation of a plurality of reflected wavefields including a reflected scalar pressure wavefield, and a reflected vector particle velocity wavefield;
- (ii) applying a first boundary condition of pressure equal zero at said artificial boundary for solution of an acoustic wave equation, thereby generating said reflected scalar pressure wavefield in-phase with said advancing acoustic wave, and said reflected vector particle velocity wavefield out-of-phase with said advancing acoustic wave;
- (iii) applying a boundary condition of particle velocity equal zero at said artificial boundary for solution of said acoustic wave equation, thereby generating said reflected scalar pressure wavefield out-of-phase with said advancing acoustic wave, and said reflected vector particle velocity wavefield in-phase with said advancing acoustic wave and
- (iv) averaging said reflected pressure wavefields and averaging said particle velocity wavefields generated in steps (ii) and (iii) over said staggered grid region to essentially cancel reflections from said artificial boundary.
A program storage device, readable by a computer, tangibly embodying a program of instructions executable by said computer to perform method steps for canceling reflections from an artificial boundary when computing propagation of an advancing acoustic wave over a large scale earth velocity model in said computer, wherein said velocity model is limited by said artificial boundary, and wherein space in said velocity model is divided into a large number of elementary grid blocks, said method comprising:
- (a) defining a staggered grid region within said large scale velocity model for computing propagation of a plurality of wavefields, wherein said plurality of wavefields includes a reflected scalar pressure wavefield, and a reflected vector particle velocity wavefield;
- (b) wherein said staggered grid region is located adjacent to said artificial boundary, and is much smaller compared to the interior of said velocity model;
- (c) applying a first boundary condition of pressure equal zero at said artificial boundary for solution of an acoustic wave equation, thereby generating said reflected scalar pressure wavefield in-phase with said advancing acoustic wave and said reflected vector particle velocity wavefield out-of-phase with said advancing acoustic wave;
- (d) applying a second boundary condition of particle velocity equal zero at said artificial boundary for solution of said acoustic wave equation, thereby generating said reflected scalar pressure wavefield out-of-phase with said advancing acoustic wave and said reflected vector particle velocity wavefield in-phase with said advancing acoustic wave;
- and
- (e) averaging said reflected pressure wavefield, and averaging said vector particle velocity wavefield generated in steps (c) and (d) over said staggered grid region to essentially cancel reflections from said artificial boundary.
Dependencies are inferred from the claim text.