3-D converted shear wave rotation with layer stripping

• US 5,610,875 A
• Filed: 09/29/1995
• Issued: 03/11/1997
• Est. Priority Date: 09/29/1995
• Status: Expired due to Term

• Current Assignee: WesternGeco LLC (Schlumberger NV)
• Original Assignee: Western Atlas International Incorporated (Northrop Grumman Corporation)
• Inventors: Gaiser, James E.
• Primary Examiner(s): Lobo, Ian J.
• Application Number:
  US08/536,429
• Time in Patent Office: 529 Days
• Field of Search: 367/37, 367/38, 367/75, 367/56, 364/421, 175/40, 175/45, 175/50
• US Class Current: 367/75
• CPC Class Codes: G01V 1/286 Mode conversion

Abstract

A method for using a compressional-wave source to produce converted shear waves which are subjected to Alford 4-component rotation to align the observation coordinates with the natural coordinates of the principal anisotropic axes of a birefringent formation. The static time shift between the fast and the slow shear wavefields due to shear-wave splitting is determined so that they can be synchronized thereby to isotropize the birefringent formation. From those data, fracture-plane orientation can be determined. Based on those data, the direction of a deviated bore hole is aligned perpendicular to the fracture plane strike. For a deep-seated target formation, shallower layers are isotropized prior to rotation and synchronization of the converted shear wavefields originating from that deeper formation.

Claims

1. A method for defining the geographic orientation of the principal axes of a birefringent medium for use in selecting the direction of deviation of a horizontal borehole relative to the strike of a lithologic fracture plane, comprising:

radiating a compressional wavefield from a first source location;

from a first receiver station spaced apart from said first source location along a first azimuth, detecting first and second seismic signal components representative of fast and slow converted shear-wave energy reflected from a subsurface conversion point associated with an acoustic interface beneath a birefringent medium;

radiating a compressional wavefield from a second source location;

from a second receiver station spaced apart from said second source location along a second azimuth substantially perpendicular to said first azimuth, detecting third and fourth signal components representative of fast and slow converted shear-wave energy reflected from said subsurface conversion point;

programming a digital computer to format said first, second, third, and fourth signal components into a 2×

2 matrix, said matrix including diagonal and off diagonal terms;

with the aid of said computer, applying a four-term rotation operator to said matrix which minimizes the energy content of the off-diagonal terms of said matrix to define the orientation of the principal polarization axes of said birefringent medium as the fracture-plane strike; and

deviating a borehole perpendicular to the so-defined strike.

2. The method as defined by claim 1, comprising:

measuring the time lag between the fast and slow diagonal shear-wave terms of said matrix;

applying a static time shift to one of the source components to compensate for the time lag between said fast and slow wavefields thereby to isotropize said birefringent medium; and

imaging said acoustic interface.

3. The method as defined by claim 2, comprising:

downward continuing the imaging of successively deeper acoustic interfaces by isotropizing shallower birefringent media.

4. The method as defined by claim 2, comprising:

in a target formation characterized by substantially vertical fracture planes, defining the axial orientation of said fracture planes by isotropizing overlying birefringent formations.

5. A method for seismic exploration, comprising:

radiating a compressional wavefield from a first source location;

detecting first and second signal components at a first receiver station spaced-apart from said first source location along a first azimuth, said first and second signal components being representative of fast and slow converted shear-wave energy reflected from a conversion point associated with an acoustic interface beneath a birefringent medium;

radiating a compressional wavefield from a second source location;

detecting third and fourth signal components at a second receiver station spaced-apart from said second source location along a second azimuth substantially perpendicular to said first azimuth, said third and fourth signal components being representative of fast and slow converted shear-wave energy reflected from said conversion point;

programming a computer to create a 2×

2 matrix from said first, second, third and fourth signal components, said matrix being characterized by diagonal and off-diagonal terms;

programming said computer to minimize the off-diagonal terms by applying an Alford rotation operator to said matrix to define the natural coordinate frame of the principal axes of birefringence;

programming said computer to measure the time difference between the fast and the slow shear wavefields represented by the diagonal terms of said rotated matrix and applying a corresponding static time shift to one of the source components to isotropize said birefringent medium; and

imaging said acoustic interface.

6. The method as defined by claim 5, comprising:

downward-continuing the imaging of successively deeper acoustic interfaces by isotropizing superjacent birefringent media.