Efficient and stable absorbing boundary condition in finite-difference calculations

Patent metadata

Field	Value
Publication number	US-10317546-B2
Application number	US-201514968961-A
Country	US
Kind code	B2
Filing date	Dec 15, 2015
Priority date	Feb 13, 2015
Publication date	Jun 11, 2019
Grant date	Jun 11, 2019

Abstract

Method for generating an effective, efficient, and stable absorbing boundary condition in finite-difference calculations, such as model-simulation of predicted seismic data. The top surface and optionally the bottom surface of the computational domain or grid are treated with one or more layers of PML (51), preferably 1D PML, assuming an orthorhombic medium in the PML implementation (52). The side surfaces are handled with one or more ABC layers (53). Further advantages may be realized by tapering earth model symmetry axis on the top and bottom of the model toward the vertical (54). The invention provides a beneficial compromise between reducing artifacts in the image or physical property model and computational efficiency and stability.

First claim

The invention claimed is:

1. A method for exploring for hydrocarbons in a subsurface region, comprising:
simulating a seismic waveform, using a computer and a model of subsurface velocity or other physical property affecting propagation of seismic waves in the subsurface region, wherein computations are performed on a finite computational grid representing a subsurface region, said finite computational grid using absorbing condition boundary conditions to minimize non-physical wave reflections at grid boundaries, said absorbing condition boundary conditions being PML boundary conditions for top surface of the finite computational grid, and ABC boundary conditions for side surfaces of the finite computational grid;
using the simulated waveform in performing full waveform inversion or reverse time migration of seismic data;
displaying (i) an inferred physical property model from the inversion or (ii) a subsurface image from the migration; and
using the (i) inferred physical property model or (ii) the subsurface image to prospect for hydrocarbons in the subsurface region.

2. The method of claim 1 , wherein the PML boundary conditions are 1D PML.

3. The method of claim 2 , wherein said model of subsurface velocity or other physical property affecting propagation of seismic waves is anisotropic with an axis of symmetry tilted from vertical, and further comprising padding the finite computational grid at the top or bottom surface, or both, to taper the axis of symmetry at said top or bottom surfaces toward vertical.

4. The method of claim 1 , wherein a plurality of absorbing layers is used on the top and bottom surfaces or the side surfaces or both.

5. The method of claim 1 , further comprising applying ABC or PML boundary conditions to a bottom surface of the finite computational grid.

6. The method of claim 1 , wherein an orthorhombic subsurface medium is assumed for PML implementation at the top surface and, if PML is used at a bottom surface of the finite computational grid, an orthorhombic medium is assumed there as well.

7. The method of claim 1 , further comprising jointly using z—PML and x,y—ABC boundary conditions at corners of the finite computational grid.

8. The method of claim 1 , wherein said model of subsurface velocity or other physical property affecting propagation of seismic waves is assumed to be isotropic in a zone at the top surface, or at a bottom surface of the finite computational grid, or both.

9. The method of claim 1 , wherein said model of subsurface velocity or other physical property affecting propagation of seismic waves is anisotropic with transversely isotropic or orthorhombic symmetry in a zone at the top surface, or at a bottom surface of the finite computational grid, or both.

Assignees

• Brytik Valeriy
• Shaw Jason
• Jing Charlie
• Zhao Hong
• Anderson John E
• Exxonmobil Upstream Res Co

Inventors

• Brytik Valeriy
• Shaw Jason
• Jing Charlie
• Zhao Hong
• Anderson John E

Classifications

• G06F17/10

  Complex mathematical operations {(function generation by table look-up G06F1/03; evaluation of elementary functions by calculation G06F7/544)} · CPC title

• G01V1/282 Primary

  Application of seismic models, synthetic seismograms · CPC title

• G01V99/005

  Physics · mapped topic

• G01V20/00 Primary

  Geomodelling in general · CPC title