Data processing: structural design – modeling – simulation – and em – Modeling by mathematical expression
Reexamination Certificate
2000-03-24
2004-02-03
Frejd, Russell (Department: 2123)
Data processing: structural design, modeling, simulation, and em
Modeling by mathematical expression
C703S010000, C367S073000, C702S014000, C711S173000
Reexamination Certificate
active
06687659
ABSTRACT:
The present invention relates to migration of seismic reflections performed in a high-speed digital computer, and more particularly to suppressing reflection artifacts imposed on a numerical model by artificial boundaries which limit computational size.
BACKGROUND OF THE INVENTION
For many years seismic exploration for oil and gas has involved the use of a source of seismic energy and its reception by an array of seismic detectors, generally referred to as geophones when used on land, and as hydrophones when used offshore. On land the source of seismic energy can be a high explosive charge electrically detonated in a borehole located at a selected point on a terrain, or another energy source having capacity for delivering a series of impacts or mechanical vibrations to the earth's surface. The acoustic waves generated in the earth by these sources are transmitted back from strata boundaries and reach the surface of the earth at varying intervals of time, depending on the distance and the characteristics of the subsurface traversed. These returning waves are detected by the geophones, which function to transduce such acoustic waves into representative electrical signals. In use an array of geophones is generally laid out along a line to form a series of observation stations within a desired locality, the source injects acoustic signals into the earth, and the detected signals are recorded for later processing using digital computers where the data are generally quantized as digital sample points such that each sample point may be operated on individually. Accordingly, seismic field records are reduced to vertical and horizontal cross sections which approximate subsurface features. The geophone array is then moved along the line to a new position and the process repeated to provide a seismic survey. More recently seismic surveys involve geophones and sources laid out in generally rectangular grids covering an area of interest so as to expand areal coverage and enable construction of three dimensional (3D) views of reflector positions over wide areas.
The general principle of offshore prospecting consists of using a seismic source to create a disturbance in the marine environment (e.g. by releasing a volume of air or steam into the water, by varying the volume of an immersed body, by an implosion, etc), and in using hydrophone detectors towed by a prospecting ship, or by geophones placed on the seabed, to obtain seismic data for extracting useful information concerning the geology of the subsoil. Typically in offshore exploration, a long streamer cable that electrically connects multiple hydrophones is towed behind the prospecting ship. Acoustic waves are periodically generated in the water, reflected from subterranean earth layers, and detected by the hydrophones. The hydrophones convert the detected seismic waves into representative electrical signals, which may be processed on the ship, and/or recorded on a storage medium such as magnetic tape for later processing.
As oil and gas production from known reservoirs and producing provinces declines, explorationists seek new areas in which to find hydrocarbons. Many of the new areas under exploration contain complex geological structures that are difficult to image with 2D techniques. Accordingly, 3D seismic processing has come into common use for mapping subterranean structures associated with oil and gas accumulations. Geophysicists, however, are well aware that a 2D seismic record section or 3D view is not a true reflectivity from the earth, but is instead a transformation of the earth's reflectivity into a plane where each recorded event is located vertically beneath the source/receiver midpoint. Steep dip, curved surfaces, buried foci, faults and other discontinuities in subterranean structure each contribute their unique characteristics to the seismic record and, in complexly faulted and folded areas, make interpretation of the geological layering from the seismic record extremely difficult. Migration is the inverse transformation that carries the plane of recorded events into a true 3D reflectivity of the earth, thereby placing reflections from dipping beds in their correct location and collapsing diffractions.
Of the various available migration methods, wave-equation migration is considered to be superior because it is based on accurate propagation of seismic waves through complex models of the earth. Wave equation computations using numerical techniques have led to a procedure called reverse time migration (RTM). By this procedure, the wavefield recorded at the surface is imaged in depth using a model of earth velocities in a numerical solution of the wave-equation. The wavefield at the surface is used as a boundary condition for the numerical computations. Proceeding by inserting the data at the surface of a computational grid for each record time step, starting with the last recorded time sample, and ending with the first, the wavefield migrates to the position from which the reflections originated.
Recently, numerical finite-difference methods have been widely used for obtaining seismic images based on pre-stack and post-stack RTM techniques, and for obtaining synthetic seismograms for studying exploration related problems. A practical procedure for doing reverse time migration is disclosed in a publication, Mufti, I. R., et al, “Finite-Difference Depth Migration of Exploration-Scale 3D Seismic Data,” Geophysics, Vol. 61. No. 3 (May-June 1996), which is incorporated herein by reference. RTM, which requires enormous computer resources as compared to simpler or less accurate migration algorithms, has recently been applied to 3D seismic data. An improved image results from accuracy (dynamic as well as kinematic) of the finite-difference method over conventional normal-moveout and raytracing based seismic imaging methods. In this RTM procedure a finite-difference earth model, which is based on the best estimate of subsurface velocities, is required. This involves dividing the model simulation space into a large number of elementary grid blocks, and assigning a velocity value to each grid.
A common problem with the finite-difference migration method is that simulation of wave propagation in an extensive portion of the earth must be modeled with limited computer resources, i.e., mainly limited central RAM memory in the computer such that the computational size of the model is truncated due to limited computer memory size. Accordingly, artificial boundaries that act as perfect reflectors are imposed by computer memory limitations, and if not properly handled cause unwanted reflection artifacts. These unwanted reflection artifacts may disrupt the image in complex ways and are not easily removed after the modeling and imaging processes are carried out.
One method of absorbing advancing waves in a 3D computational grid is an extensively used method known as a one-way wave equation (1WWE), proposed by Clayton and Engquist, 1980 “Absorbing Boundary Conditions for Wave-Equation Migration”, Geophysics, Vol. 45 pp 895-904, and Keys, 1985 “Absorbing Boundary Conditions for Acoustic Media”, Geophysics, Vol. 50 pp 89-902. This 1WWE method, however, is effective for absorbing advancing waves arriving only at or near a selective incident angle. More recently, much attention has been given to a method called the perfectly matched layer (PML) method, which was developed by Berenger (1994) “A Perfectly Matched Layer for the Absorpiton of Electromatic Waves”, Journal of Computational Physics Vol. 114, pp 185-200, for use in a finite-difference time domain computations. The PML method, which was first applied to a 2D transverse electric wave problem, can absorb waves propagating at both near normal and high incident angles, and thus provides a reflectionless interface between the propagation region of interest and the PML layers at all incident angles of arriving waves. A PML layer can be conveniently applied on a planar face. However an edge or a corner of a rectangular finite difference numerical model volume requires a special ten
ConocoPhillips Company
Cross Ryan N.
Frejd Russell
LandOfFree
Method and apparatus for absorbing boundary conditions in... does not yet have a rating. At this time, there are no reviews or comments for this patent.
If you have personal experience with Method and apparatus for absorbing boundary conditions in..., we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Method and apparatus for absorbing boundary conditions in... will most certainly appreciate the feedback.
Profile ID: LFUS-PAI-O-3349429