CA2386195C - Obtaining an image of an underground formation - Google Patents

Obtaining an image of an underground formation Download PDF

Info

Publication number
CA2386195C
CA2386195C CA002386195A CA2386195A CA2386195C CA 2386195 C CA2386195 C CA 2386195C CA 002386195 A CA002386195 A CA 002386195A CA 2386195 A CA2386195 A CA 2386195A CA 2386195 C CA2386195 C CA 2386195C
Authority
CA
Canada
Prior art keywords
underground
seismic energy
data
component receiver
borehole
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CA002386195A
Other languages
French (fr)
Other versions
CA2386195A1 (en
Inventor
Alphonsus Paulus Eduard Ten Kroode
William Alexander Mulder
Pleun Marinus Van Der Sman
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shell Canada Ltd
Original Assignee
Shell Canada Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shell Canada Ltd filed Critical Shell Canada Ltd
Publication of CA2386195A1 publication Critical patent/CA2386195A1/en
Application granted granted Critical
Publication of CA2386195C publication Critical patent/CA2386195C/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48Processing data
    • G01V1/50Analysing data

Abstract

An image of an underground formation (2) around a borehole (1) is obtained by activating an omnidirectional source (9) and recording with a three-component receiver (10) the components of the reflected energy (15); determining therefrom the components the directions from which the energy arrives at the three-component receiver (10) as a function of two-way travel time;
selecting a first underground position (30); assuming a reflector to be present at position (30) and calculating the arrival direction of a ray (35) extending from the source (9) via the reflector (30) to the receiver (10) and two-way travel time along the ray (35);
accepting the data if the calculated arrival direction is substantially equal to an arrival direction that has the same two-way travel time, and attributing the data on the position (30); and selecting a next position (31 or 32) and repeating steps (e) and (f) until the last underground position to obtain the image of the underground formation comprising a set of reflectors attributed to positions.

Description

OBTAINING AN IMAGE OF AN UNDERGROUND FORMFTIOIQ
BACKGROUND

The present invention relates -:o a method of obtairiing an image of an underground formation around a borehole extending through the underground formation. The image that is to be obtained comprises a set of reflectors attributed to underground positions in the underground formation around the borehole. Such an image is produced to provide detailed information of the underground forrrtation while drilling the borehole. This information allows planning the direction into which the borehole is drilled. This is particularly useful when it is required that a horizontal borehole that is being drilled is kept within a thin formation layer.
USA patent specification No. 5 300 929 relates to a method of delineating an interface between a salt and a sediment. The known method comprises (a) arranging an omnidirectional source at the surface and fixedly arranging a three-component receiver in a borehole that extends through the salt;
(b) activating the omnidirectional source to generate seismic energy and recording with the three-component receiver data in the form of the components of the seismic energy;
(c) determining from the components of the seismic energy the directions from which the seismic energy arrives at the three-component receiver as a function of travel time;
(d) accepting a point as being located at the interface if the ray extending through the point has the corresponding travel time; and (e) selecting a new surface position for the omni-directional source, and repeating steps (b)-(d).
The }.nown method is applied to determine the boundary between a salt dome and sediment surrounding the salt dome, wherein the borehole is dzilled into the salt dome.
Because the omnidirectional source and the three-component receiver are spaced apart at either side of the boundary, the seismic energy passes through the underground formation.

SUMMARY
It is an object of an embodiment of the present invention to provide a method for obtaining an image of the underground formation using an omnidirectional source and an three-component receiver that are both located in a borehole extending through the underground formation, which allows imaging reflectors that reflect the seismic energy emitted by the omnidirectional source, wherein the position of the reflectors can be anywhere around the borehole.
To this end the method of obtaining an image of an underground formation around a borehole extending through the underground formation according to the present invention comprises the steps of (a) selecting a number of positions for an omni-directional source and a three-component receiver in the borehole, selecting a number of underground positions in the formation, and attributing a value of zero to the underground positions;
(b) arranging the omnidirectional source and the three-component receiver in a first position in the borehole;
(c) activating the omnidirectional source to generate seismic energy and recording with the three-component receiver data in the form of components of reflected seismic energy;
(d) determining from the components of the reflected seismic energy directions from which the reflected seismic energy arrives at the three-component receive= as a J'unction of two-way travel time;
(e) selecting a first underground position;
(f) calculating an arrival direction of a ray extending from the omnidirectional source to the first underground position and back to the three-component receiver and the two-way travel time of seismic energy passing along the ray;
(g) accepting the data if the calculated arrival direction is substantially equal to an axrival direction as obtained in step (d) pertaining to reflected seismic energy having the same two-way travel time, and adding the accepted data to the value attributed to the first underground position;
(h) selecting a next underground position and repeating steps (f) and (g) until the last underground position;
and (i) arranging the omnidirectional source and the three-component receiver in a next position in the borehole, and repeating the steps (b) through (h) until the last position along the borehole to obtain the image of the underground formation comprising a set of data mapped on underground positions.
in the specification and in the claims the term `two-way travel time' is used to refer to the time it takes for seismic energy to go from a source via a reflector to a receiver.
It will be understood that in order to perform the calculations in step (f), it is required to know the seismic velocities in the formation. These seismic velocities can be obtained from previous seismic work done in relation to the formation, or they can be obtained from core samples. In addition, sonic measurements can provide information on the seismic velocities.

3a BRIEF DESCRIPTION OF THE DRAWING

The invention will now be described by way of example in more detail with reference to the accompanying Figure which is a schematic view of a borehole being drilled in an underground formation.
DET.~ILED DESCRIPTION

In the Figure is shown the lower end of a borehole 1 that is being drilled in an underground formation 2. In this example the borehole 1 is drilled by means of a drill bit 3 suspended in the borehole by means of a drill string assembly 5, which drill string assembly 5 is rotated.
Near the drill bit 3, the drill string assembly 5 comprises an omnidirectional source 9 and a downhole three-component receiver 10.
During normal operation, the drill bit 3 is used to drill the borehole 1, and in order to obtain the image, drilling is interrupted, and the omnidirectional source 9 is activated. Seismic energy emitted by the ornni-directional source 9 spreads out into the formation 2, and the wave fronts pertaining to the reflected seismic energy at different moments in time are schematically shown by the dashed lines 15.
Assume that a reflector 18 is present in the underground formation 2, which reflector 18 reflects the.
seismic energy. The wave fronts pertaining to the reflected seismic energy at different moments in time are schematically shown by the dashed lines 20. Line 22 represents a ray extending from the omnidirectional source 9 to the reflector 18 and back to the three-component receiver 10_ The data received by the three-component receiver 10 includes the components of the reflected seismic energy in time. From this data the directions from which the seismic energy arrives can be determined as a function of the two-way travel time.
Next a number of underground positions 30, 31 and 32 are selected in the formation, and a first one is chosen, let us assume that it is underground position 30. Knowing the seismic velocities in the underground formation, the arrival direction of the ray extending from the omnidirectional source 9 to the reflector at position 30 and back to the three-component receiver 10 is calculated as well as the two-way travel time of seismic energy passing along the ray. This ray is shown by dashed 5 line 35.
The data is accepted if the calculated arrival direction is substantially equal to an arrival direction that has the same two-way travel time. In this example, this is clearly not the case for underground position 30, so that the data is not accepted.
Then a next underground position 31 is selected. Next the arrival direction of the ray extending from the omnidirectional source 9 to the reflector at position 31 and back to the three-component receiver 10 is calculated as well as the two-way travel time of seismic energy passing along the ray. This ray coincides with line 22.
In this case the calculated arrival direction is substantially equal to the arrival direction of the reflection from the reflector 18, and the ray that coincides with line 22 has the same two-way travel time.
Thus the data is accepted and the accepted data are added to the value attributed to the underground position 31.
This adding is also referred to as migration.
Then the third underground position 32 is selected.
The calculated arrival direction of the ray extending from the omnidirectional source 9 to the reflector at position 32 and back to the three-component receiver 10 shown as dashed line 37 is not substantially equal to the arrival direction of the reflection from the reflector 18. Thus the data is not accepted.
Having treated the three underground positions 30, 31 and 32, an image of the underground formation 2 is obtained, which image comprises a reflector attributed to underground position 31 and no reflectors attributed to underground positions 30 and 32.
Drilling is resumed, and after some distance has been drilled the above described procedure is repeated, and so on.
In this way an accurate image of the underground formation near the drill bit can be obtained, in particular when more than three underground positions are selected for each position along the borehole.
In case there are more reflectors than the reflector 18 shown in the Figure, rays from these reflectors (not shown) will be received at different moments in time.
Determining from the components of the reflected seismic energy the directions from which the reflected seismic energy arrives at the three-component receiver as a function of two-way travel time is known, it can for example be done with techniques outlined in the article `Comparison of signal processing techniques for estimating the effects of anisotropy' by C Macbeth and S Crampin, Geophysical Prospecting, 39, 1991, pages 357-385.
The data is accepted if the calculated arrival direction is substantially equal to an arrival direction that has the same two-way travel time. To do so, suitably, the difference between the calculated arrival direction and the arrival direction pertaining to reflected seismic energy having the same two-way travel time is determined. Then a weight factor is determined using a predetermined function of this difference. The data are multiplied with the weight factor, and the weighted data are mapped on the underground position. The weight function is for example a rectangular window function. The window function or box function is a function of the difference, such that the window function equals 1 if the absolute value of the difference is less than a predetermined value and that it equals 0 every-where else. Consequently for a large difference the weight factor is 0 and the weighted data is 0 so that no data are mapped, and for a relatively small difference the weight factor is 1 so that the data are mapped. An alternative weight function is a cosine squared.
Suitably, the data that is mapped on the underground position is the magnitude of the reflected seismic energy, which is the sum of the components of the reflected seismic energy or the square root of the sum of the squares of the components of the reflected seismic energy. The magnitude of the reflected seismic energy is then determined with the migration technique.
Alternatively, the reflectivity can be determined from the data by comparing the reflected seismic energy with the emitted seismic energy and making a correction for the geometrical spreading.
Although it is possible to apply the migration technique for any position of the omnidirectional source relative to the three-component receiver, it is preferred that the omnidirectional source and the three-component receiver are coincident. In the specification and in the claims the word `coincident' is used as follows. Two devices are said to be coincident when they are as close to each other as is technically feasible, in which case they can be considered as one for calculation purposes.
In that case the reflectivity can be calculated using a zero-offset migration algorithm.
In order to distinguish the arrival of shear waves (or s waves) from the arrival of faster compression waves (or p waves), a sensor, such as a hydrophone or an accelerometer can be included in the three-component receiver.
The reflected seismic energy can be passed to surface by known means of transferring data, so that the analysis part of the method is done at surface. Alternatively, the directions from which the reflected seismic energy arrives at the three-component receiver as a function of two-way travel time are determined in-situ, and the results are transferred to surface where the analysis takes place.
Instead of using an omnidirectional separate source, the drill bit itself can be used as a source, and in that case the seismic energy is the noise generated while drilling.
The present invention provides a simple method for obtaining an image from an underground formation in the neighbourhood of a borehole that is being drilled.

Claims (5)

CLAIMS:
1. A method of obtaining an image of an underground formation around a borehole extending through the underground formation, which method comprises the steps of (a) selecting a number of positions for an omni-directional source and a three-component receiver in the borehole, selecting a number of underground positions in the formation, and attributing a value of zero to the underground positions;
(b) arranging the omnidirectional source and the three-component receiver in a first position in the borehole;
(c) activating the omnidirectional source to generate seismic energy and recording with the three-component receiver data in the form of components of reflected seismic energy;
(d) determining from the components of the reflected seismic energy directions from which the reflected seismic energy arrives at the three-component receiver as a function of two-way travel time;
(e) selecting a first underground position;
(f) calculating an arrival direction of a ray extending from the omnidirectional source to the first underground position and back to the three-component receiver and the two-way travel time of seismic energy passing along the ray;
(g) accepting the data if the calculated arrival direction is substantially equal to an arrival direction as obtained in step (d) pertaining to reflected seismic energy having the same two-way travel time, and adding the accepted data to the value attributed to the first underground position;

(h) selecting a next underground position and repeating steps (f) and (g) until the last underground position;
and (i) arranging the omnidirectional source and the three-component receiver in a next position in the borehole, and repeating the steps (b) through (h) until the last position along the borehole to obtain the image of the underground formation comprising a set of data mapped on underground positions.
2. The method according to claim 1, wherein step (g) comprises determining a difference between the calculated arrival direction and the arrival direction as obtained in step (c) pertaining to reflected seismic energy having the same two-way travel time, multiplying the data with a weight factor which is a predetermined function of this difference, and mapping the weighted data on the first underground position.
3. The method according to claim 1 or 2, wherein the data used in step (g) is the sum of the components of the reflected seismic energy recorded in step (c).
4. The method according to any one of the claims 1-3, wherein the omnidirectional source and the three-component receiver are coincident.
5. The method according to any one of the claims 1-4, wherein the three-component receiver further includes a pressure sensor.
CA002386195A 1999-10-14 2000-10-12 Obtaining an image of an underground formation Expired - Fee Related CA2386195C (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP99308119 1999-10-14
EP99308119.9 1999-10-14
PCT/EP2000/010080 WO2001027657A1 (en) 1999-10-14 2000-10-12 Obtaining an image of an underground formation

Publications (2)

Publication Number Publication Date
CA2386195A1 CA2386195A1 (en) 2001-04-19
CA2386195C true CA2386195C (en) 2008-12-16

Family

ID=8241673

Family Applications (1)

Application Number Title Priority Date Filing Date
CA002386195A Expired - Fee Related CA2386195C (en) 1999-10-14 2000-10-12 Obtaining an image of an underground formation

Country Status (13)

Country Link
US (1) US6510105B1 (en)
EP (1) EP1221059B1 (en)
JP (1) JP2003511708A (en)
AU (1) AU767001B2 (en)
BR (1) BRPI0014706B1 (en)
CA (1) CA2386195C (en)
DE (1) DE60019778T2 (en)
EA (1) EA003526B1 (en)
GC (1) GC0000206A (en)
MY (1) MY123403A (en)
NO (1) NO20021738L (en)
OA (1) OA12511A (en)
WO (1) WO2001027657A1 (en)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6811744B2 (en) * 1999-07-07 2004-11-02 Optomec Design Company Forming structures from CAD solid models
US6510105B1 (en) * 1999-10-14 2003-01-21 Shell Oil Company Obtaining an image of an underground formation
US6904368B2 (en) * 2002-11-12 2005-06-07 Landmark Graphics Corporation Seismic analysis using post-imaging seismic anisotropy corrections
US7207215B2 (en) * 2003-12-22 2007-04-24 Halliburton Energy Services, Inc. System, method and apparatus for petrophysical and geophysical measurements at the drilling bit
US7027926B2 (en) * 2004-04-19 2006-04-11 Pathfinder Energy Services, Inc. Enhanced measurement of azimuthal dependence of subterranean parameters
GB2420881B (en) * 2004-12-01 2008-01-16 Westerngeco Ltd Processing seismic data
US20070223822A1 (en) * 2006-03-20 2007-09-27 Pathfinder Energy Services, Inc. Data compression method used in downhole applications
EP2165218B1 (en) 2007-06-15 2017-02-22 Baker Hughes Incorporated Imaging of formation structure ahead of the drill-bit
US7558675B2 (en) * 2007-07-25 2009-07-07 Smith International, Inc. Probablistic imaging with azimuthally sensitive MWD/LWD sensors
US8428880B2 (en) * 2007-12-14 2013-04-23 Shell Oil Company Method of processing data obtained from seismic prospecting
US7791981B2 (en) * 2008-05-15 2010-09-07 Shell Oil Company Velocity analysis for VSP data
US8600115B2 (en) 2010-06-10 2013-12-03 Schlumberger Technology Corporation Borehole image reconstruction using inversion and tool spatial sensitivity functions
WO2012064610A2 (en) * 2010-11-08 2012-05-18 Schlumberger Canada Limited System and method for communicating data between wellbore instruments and surface devices
US9658360B2 (en) 2010-12-03 2017-05-23 Schlumberger Technology Corporation High resolution LWD imaging
CN102879805B (en) * 2012-10-24 2015-06-24 北京市市政工程研究院 Borehole-based and ground combined seismic wave space exploration method

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5230387A (en) * 1988-10-28 1993-07-27 Magrange, Inc. Downhole combination tool
US5170377A (en) 1990-03-22 1992-12-08 Western Atlas International, Inc. 3-D mapping of salt domes
US5081611A (en) 1991-03-06 1992-01-14 Schlumberger Technology Corporation Methods for determining formation and borehole parameters via two-dimensional tomographic reconstruction of formation slowness
US5300929A (en) 1991-10-04 1994-04-05 Chevron Research And Technology Company Method for delineating an anomalous geologic structure
US6023657A (en) 1997-10-15 2000-02-08 Input/Output, Inc. Seismic imaging using omni-azimuth seismic energy sources and directional sensing
US6510105B1 (en) * 1999-10-14 2003-01-21 Shell Oil Company Obtaining an image of an underground formation

Also Published As

Publication number Publication date
AU1024701A (en) 2001-04-23
NO20021738D0 (en) 2002-04-12
CA2386195A1 (en) 2001-04-19
AU767001B2 (en) 2003-10-30
WO2001027657A1 (en) 2001-04-19
US6510105B1 (en) 2003-01-21
DE60019778D1 (en) 2005-06-02
BR0014706A (en) 2002-06-18
DE60019778T2 (en) 2005-11-17
BRPI0014706B1 (en) 2016-06-28
JP2003511708A (en) 2003-03-25
EP1221059A1 (en) 2002-07-10
NO20021738L (en) 2002-04-12
MY123403A (en) 2006-05-31
EA003526B1 (en) 2003-06-26
EA200200442A1 (en) 2002-10-31
OA12511A (en) 2006-05-29
GC0000206A (en) 2006-03-29
EP1221059B1 (en) 2005-04-27

Similar Documents

Publication Publication Date Title
CA2386195C (en) Obtaining an image of an underground formation
US6026913A (en) Acoustic method of connecting boreholes for multi-lateral completion
US5899958A (en) Logging while drilling borehole imaging and dipmeter device
US7035165B2 (en) Imaging near-borehole structure using directional acoustic-wave measurement
CA2548131C (en) Shear wave velocity determination using evanescent shear wave arrivals
US7924652B2 (en) Method for determining seismic anisotropy
CA2747275C (en) Method and apparatus for measuring formation anisotropy while drilling
US10782433B2 (en) Method for an automatic detection of acoustic reflectors and their parameters from borehole acoustic array data
WO2008083032A2 (en) Imaging near-borehole reflectors using shear wave reflections from a multi-component acoustic tool
GB2333155A (en) Acoustic measurement while drilling system using compressional waves reflected from within formation
CA2708096C (en) Method and system for delineating a second wellbore from a first wellbore
CA2465028C (en) Borehole sonic data processing method
CA2992437C (en) Passive ranging using acoustic energy originating from a target wellbore
Nakken et al. A new MWD concept for geological positioning of horizontal wells
JP2816899B2 (en) Underground structure evaluation method by three-dimensional particle motion analysis of released elastic waves during well drilling
Greenberg Seismic while drilling keeps bit turning to right while acquiring key real-time data
Deltombe et al. Combined processing of BHTV travel time and amplitude images
Hakim et al. Subsurface velocity measurement
Calvert Well, Where Are We Now? How to Position Wells without Good Velocity or Good Migration
Leggett et al. A new method for remote sensing of critical Stratigraphic bed boundaries and reservoir positioning of horizontal wells
Juhlin et al. Vertical seismic profiling and integration with reflection seismic studies at Laxemar, 2000
JPH06342079A (en) Analysis method for s-wave field affected by anisotropy of ground

Legal Events

Date Code Title Description
EEER Examination request
MKLA Lapsed

Effective date: 20181012