CA2313282C - Method and apparatus for determining logging tool displacements - Google Patents
Method and apparatus for determining logging tool displacements Download PDFInfo
- Publication number
- CA2313282C CA2313282C CA002313282A CA2313282A CA2313282C CA 2313282 C CA2313282 C CA 2313282C CA 002313282 A CA002313282 A CA 002313282A CA 2313282 A CA2313282 A CA 2313282A CA 2313282 C CA2313282 C CA 2313282C
- Authority
- CA
- Canada
- Prior art keywords
- logging tool
- displacements
- time
- gravitational acceleration
- displacement
- 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
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/02—Determining slope or direction
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/02—Determining slope or direction
- E21B47/022—Determining slope or direction of the borehole, e.g. using geomagnetism
Abstract
A method of determining the displacements of a logging tool during a measurement interval of the logging tool in a borehole includes obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval. The method further includes calculating a lower bound for the displacements of the logging tool during the measurement interval when the initial velocity and the gravitational acceleration are unknown. The lower bound on the displacements of the logging tool is used to flag the validity of the measurements made by the logging tool.
Description
METHOD AND APPARATUS FOR DETERMINING LOGGING TOOL
DISPLACEMENTS
BACKGROUND OF THE INVENTION
Well logging involves recording data related to one or more characteristics of a subterranean formation penetrated by a borehole as a function of depth. The record is called a log. Many types of logs are recorded by appropriate downhole instruments placed in a housing called a sonde. The sonde is lowered into the borehole on the end of a cable, and the parameters being logged are measured as the sonde is moved along the borehole. Data signals from the sonde are transmitted through the cable to the surface, where the log is made. Figure 1 shows an example of a sonde 2 that measures properties of formation 4 surrounding a borehole 6 using the principles of nuclear magnetic resonance (NMR). The NMR sonde 2 includes a magnet assembly 8 and an antenna 10. The magnet assembly 8 produces a static magnetic field Bo in all regions surrounding the sonde 2, and the antenna 10 produces an oscillating magnetic field B1 that is perpendicular and superimposed on the static magnetic field Bo. The NMR signal comes from a small resonance volume 12 which has a radial thickness that is proportional to the magnitude of the oscillating magnetic field B1 and inversely proportional to the gradient of the static magnetic field Bo. The NMR
sonde 2 makes measurements by magnetically tipping the nuclear spins of protons in the formation with a pulse of the oscillating magnetic field, and then detecting the precession of the tipped particles in the resonance volume 12.
As the NMR sonde 2 traverses the borehole 6 to make measurements, it experiences random accelerations due to borehole forces acting on it. These random accelerations result in displacements of the sonde, which may adversely affect the quality of the log. To further explain this point, the resonance volume 12 generally consists of la , CA 02313282 2000-06-30 PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/01 1001 ) thin cylindrical shells that define a sensitive region extending along the length of the sonde 2 and having a radial thickness of about 1 millimeter. If the NMR sonde
DISPLACEMENTS
BACKGROUND OF THE INVENTION
Well logging involves recording data related to one or more characteristics of a subterranean formation penetrated by a borehole as a function of depth. The record is called a log. Many types of logs are recorded by appropriate downhole instruments placed in a housing called a sonde. The sonde is lowered into the borehole on the end of a cable, and the parameters being logged are measured as the sonde is moved along the borehole. Data signals from the sonde are transmitted through the cable to the surface, where the log is made. Figure 1 shows an example of a sonde 2 that measures properties of formation 4 surrounding a borehole 6 using the principles of nuclear magnetic resonance (NMR). The NMR sonde 2 includes a magnet assembly 8 and an antenna 10. The magnet assembly 8 produces a static magnetic field Bo in all regions surrounding the sonde 2, and the antenna 10 produces an oscillating magnetic field B1 that is perpendicular and superimposed on the static magnetic field Bo. The NMR signal comes from a small resonance volume 12 which has a radial thickness that is proportional to the magnitude of the oscillating magnetic field B1 and inversely proportional to the gradient of the static magnetic field Bo. The NMR
sonde 2 makes measurements by magnetically tipping the nuclear spins of protons in the formation with a pulse of the oscillating magnetic field, and then detecting the precession of the tipped particles in the resonance volume 12.
As the NMR sonde 2 traverses the borehole 6 to make measurements, it experiences random accelerations due to borehole forces acting on it. These random accelerations result in displacements of the sonde, which may adversely affect the quality of the log. To further explain this point, the resonance volume 12 generally consists of la , CA 02313282 2000-06-30 PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/01 1001 ) thin cylindrical shells that define a sensitive region extending along the length of the sonde 2 and having a radial thickness of about 1 millimeter. If the NMR sonde
2 moves 1 millimeter or more in the radial direction, the measurements of the T2 spin-spin relaxation times of the protons may be corrupted. Also, the time during which the nuclear spins of the protons in the formation 4 are polarized by the applied magnetic fields depend on the motion of the NMR sonde 2. If the NMR sonde 2 sticks and slips while moving along the direction of the borehole, T1 relaxation-time measurements can be compromised. In another logging mode which estimates the bound fluid volume by first saturating the nuclear spins and then letting them recover during a small time, the measurement mode overestimates the bound fluid volume if the tool moves faster than expected along the longitudinal axis of the borehole 6, or if the tool is radially displaced by more than 1 millimeter during the recovery period.
If the displacements of the sonde during the measurement interval are known, then the portions of the NMR measurements that are distorted by motions of the sonde can be identified and discarded or corrected using appropriate compensation methods.
Prior art methods have used a motion detection device, such as a strain gauge, an ultrasonic range finder, an accelerometer, or a magnetometer, to detect the motions of a sonde during a logging operation. In this manner, the motion detection device is used to establish a threshold for evaluating the quality of the log. For example, U.S.
Patent 6,051,973 issued to Prammer discloses using accelerometers to monitor peak acceleration values of a logging tool during a measurement interval of the logging tool.
The quality of the log is improved by discarding the measurements made during the period that the peak accelerations indicate that the logging tool may have been displaced by more than allowable by the extent of the sensitive region.
SUMMARY OF THE INVENTION
In one aspect, the invention is a method for determining the displacements of a logging tool during a measurement interval of the logging tool in a borehole.
The method comprises obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval and double integrating the set of accelerometer signals to obtain corresponding displacements of the logging tool as a function of the initial velocity of the logging tool and the gravitational acceleration, wherein the initial velocity of the logging tool and the gravitational acceleration are unknown. The method further comprises assuming a set of feasible initial velocities for the logging tool. For each feasible initial velocity, the method includes estimating the gravitational acceleration, calculating the displacements of the logging tool using the feasible initial velocity and the estimated gravitational acceleration, and determining the maximum of the calculated displacements.
The lower bound on the displacements of the logging tool is set to the minimum of the maximum of the calculated displacements.
In another aspect, a method for determining the displacements of a logging tool during a measurement interval of the logging tool in a borehole comprises obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval and calculating a tool displacement as a time-series from the accelerometer signals. The method further includes constructing a unique quadratic polynomial of time from the displacement time-series, subtracting the unique quadratic polynomial from the displacement time-series, and setting the lower bound to the maximum of the remainder of the displacement time-series.
According to one aspect of the present invention, there is provided a method for determining displacements of
If the displacements of the sonde during the measurement interval are known, then the portions of the NMR measurements that are distorted by motions of the sonde can be identified and discarded or corrected using appropriate compensation methods.
Prior art methods have used a motion detection device, such as a strain gauge, an ultrasonic range finder, an accelerometer, or a magnetometer, to detect the motions of a sonde during a logging operation. In this manner, the motion detection device is used to establish a threshold for evaluating the quality of the log. For example, U.S.
Patent 6,051,973 issued to Prammer discloses using accelerometers to monitor peak acceleration values of a logging tool during a measurement interval of the logging tool.
The quality of the log is improved by discarding the measurements made during the period that the peak accelerations indicate that the logging tool may have been displaced by more than allowable by the extent of the sensitive region.
SUMMARY OF THE INVENTION
In one aspect, the invention is a method for determining the displacements of a logging tool during a measurement interval of the logging tool in a borehole.
The method comprises obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval and double integrating the set of accelerometer signals to obtain corresponding displacements of the logging tool as a function of the initial velocity of the logging tool and the gravitational acceleration, wherein the initial velocity of the logging tool and the gravitational acceleration are unknown. The method further comprises assuming a set of feasible initial velocities for the logging tool. For each feasible initial velocity, the method includes estimating the gravitational acceleration, calculating the displacements of the logging tool using the feasible initial velocity and the estimated gravitational acceleration, and determining the maximum of the calculated displacements.
The lower bound on the displacements of the logging tool is set to the minimum of the maximum of the calculated displacements.
In another aspect, a method for determining the displacements of a logging tool during a measurement interval of the logging tool in a borehole comprises obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval and calculating a tool displacement as a time-series from the accelerometer signals. The method further includes constructing a unique quadratic polynomial of time from the displacement time-series, subtracting the unique quadratic polynomial from the displacement time-series, and setting the lower bound to the maximum of the remainder of the displacement time-series.
According to one aspect of the present invention, there is provided a method for determining displacements of
3 a logging tool during a measurement interval of the logging tool in a borehole, the method comprising: obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval; double integrating the set of accelerometer signals to obtain corresponding displacements of the logging tool as a function of an initial velocity of the logging tool and a gravitational acceleration, wherein the initial velocity of the logging tool and the gravitational acceleration are unknown; determining a set of initial velocities for the logging tool; for each initial velocity, estimating the gravitational acceleration, calculating the displacements of the logging tool using the initial velocity and the estimated gravitational acceleration, and determining a maximum of the calculated displacements; and setting a lower bound on the displacements of the logging tool to a value which is a minimum value of a set of the maximum calculated displacements for the set of initial velocities.
According to another aspect of the present invention, there is provided a method for improving the quality of measurements made by a logging tool during a measurement interval in a borehole, the method comprising:
obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval; double integrating the set of accelerometer signals to obtain corresponding displacements of the logging tool as a function of an initial velocity of the logging tool and a gravitational acceleration, wherein the initial velocity of the logging tool and the gravitational acceleration are unknown; determining a set of initial velocities for the logging tool; for each initial velocity, 3a estimating the gravitational acceleration, calculating the displacements of the logging tool using the initial velocity and the estimated gravitational acceleration, and determining a maximum of the calculated displacements;
estimating a lower bound for the displacements of the logging tool by selecting a value which is a minimum value of a set of the maximum calculated displacements for the initial velocities; and raising a flag if the lower bound for the displacements of the logging tool exceeds a selected threshold.
According to still another aspect of the present invention, there is provided a method for logging a well, comprising: moving a logging tool along a borehole to make measurements in a formation surrounding the borehole;
recording the measurements made by the logging tool;
measuring a set of accelerometer signals of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval; double integrating the set of accelerometer signals to obtain corresponding displacements of the logging tool as a function of an initial velocity of the logging tool and a gravitational acceleration, wherein the initial velocity of the logging tool and the gravitational acceleration are unknown;
determining a set of initial velocities for the logging tool; for each initial velocity, estimating the gravitational acceleration, calculating the displacements of the logging tool using the initial velocity and the estimated gravitational acceleration, and determining a maximum of the calculated displacements; estimating a lower bound for the displacements of the logging tool by selecting a value which is a minimum value of a set of the maximum calculated displacements for the initial velocities; and 3b raising a flag if the lower bound for the displacements of the logging tool exceeds a selected threshold.
According to yet another aspect of the present invention, there is provided a method far determining displacements of a logging tool during a measurement interval of the logging tool in a borehole, the method comprising: obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval; calculating a tool displacement as a time-series from the accelerometer signals; constructing a unique quadratic polynomial of time from the displacement time-series; and subtracting the unique quadratic polynomial from the displacement time-series; and setting a lower bound to the maximum of the remainder of the displacement time-series.
According to a further aspect of the present invention, there is provided a method for improving the quality of measurements made by a logging tool during a measurement interval in a borehole, the method comprising:
obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval; calculating a tool displacement as a time-series from the accelerometer signals; constructing a unique quadratic polynomial of time from the displacement time-series; subtracting the unique quadratic polynomial from the displacement time-series; and setting a lower bound to the maximum of the remainder of the displacement time-series;
and raising a flag if the lower bound for the displacements of the logging tool exceeds a selected threshold.
3c Other aspects and advantages of the invention will be apparent from the following description and the appended claims.
BRIEF DESCRIPTION OF THE DRAWINGS
Figure 1 shows a logging tool suspended in a borehole.
Figure 2 is a cross section of a logging tool suspended in a borehole according to one embodiment of the invention.
Figure 3 depicts a horizontal cross section of the logging tool shown in Figure 2.
3d PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) Figure 4 is a flow chart illustrating a method for determining the displacements of a logging tool according to one embodiment of the invention.
Figure 5 is a flow chart illustrating a method for determining the displacements of a logging tool according to another embodiment of the invention.
DETAILED DESCRIPTION OF THE INVENTION
Embodiments of the invention provide a method for determining displacements of a logging tool during a measurement interval along three orthogonal axes of the logging tool. In general, an accelerometer is used to measure the accelerations of the logging tool along the three orthogonal axes of the logging tool during the measurement interval. The accelerations acquired by the accelerometer, as will be further explained below, have a gravitational portion that is due to gravitational forces acting on the test-mass of the accelerometer and a kinetic portion that is due to the net force acting on the logging tool.
The displacements of the logging tool are determined from the estimated kinetic portion of the accelerations.
The displacements of the logging tool may be used to assess the quality of the measurements made by the logging tool. For example, pulse-echo nuclear magnetic resonance (NN>R) measurements are time-lapse measurements. For the measurement to be accurate, the sensitive zone of the NMR logging tool needs to substantially overlap with itself through out the measurement duration. Thus, accuracy of NMR
logging tools are sensitive to the displacement of the tool during the measurement interval.
By determining the displacements of the 'logging tool during a measurement interval, the validity of the measurements made can be verified. Of course, the invention is not limited to NMR logging tools, but is generally applicable to any logging tool that makes measurements that are sensitive to tool motion.
Various embodiments of the invention will now be discussed with reference to the accompanying figures. In order to fully understand the invention, it is helpful to consider a specific configuration of a logging tool. However, it should be clear that the invention is not limited to the specific configuration of the logging tool discussed herein. Figure 2
According to another aspect of the present invention, there is provided a method for improving the quality of measurements made by a logging tool during a measurement interval in a borehole, the method comprising:
obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval; double integrating the set of accelerometer signals to obtain corresponding displacements of the logging tool as a function of an initial velocity of the logging tool and a gravitational acceleration, wherein the initial velocity of the logging tool and the gravitational acceleration are unknown; determining a set of initial velocities for the logging tool; for each initial velocity, 3a estimating the gravitational acceleration, calculating the displacements of the logging tool using the initial velocity and the estimated gravitational acceleration, and determining a maximum of the calculated displacements;
estimating a lower bound for the displacements of the logging tool by selecting a value which is a minimum value of a set of the maximum calculated displacements for the initial velocities; and raising a flag if the lower bound for the displacements of the logging tool exceeds a selected threshold.
According to still another aspect of the present invention, there is provided a method for logging a well, comprising: moving a logging tool along a borehole to make measurements in a formation surrounding the borehole;
recording the measurements made by the logging tool;
measuring a set of accelerometer signals of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval; double integrating the set of accelerometer signals to obtain corresponding displacements of the logging tool as a function of an initial velocity of the logging tool and a gravitational acceleration, wherein the initial velocity of the logging tool and the gravitational acceleration are unknown;
determining a set of initial velocities for the logging tool; for each initial velocity, estimating the gravitational acceleration, calculating the displacements of the logging tool using the initial velocity and the estimated gravitational acceleration, and determining a maximum of the calculated displacements; estimating a lower bound for the displacements of the logging tool by selecting a value which is a minimum value of a set of the maximum calculated displacements for the initial velocities; and 3b raising a flag if the lower bound for the displacements of the logging tool exceeds a selected threshold.
According to yet another aspect of the present invention, there is provided a method far determining displacements of a logging tool during a measurement interval of the logging tool in a borehole, the method comprising: obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval; calculating a tool displacement as a time-series from the accelerometer signals; constructing a unique quadratic polynomial of time from the displacement time-series; and subtracting the unique quadratic polynomial from the displacement time-series; and setting a lower bound to the maximum of the remainder of the displacement time-series.
According to a further aspect of the present invention, there is provided a method for improving the quality of measurements made by a logging tool during a measurement interval in a borehole, the method comprising:
obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval; calculating a tool displacement as a time-series from the accelerometer signals; constructing a unique quadratic polynomial of time from the displacement time-series; subtracting the unique quadratic polynomial from the displacement time-series; and setting a lower bound to the maximum of the remainder of the displacement time-series;
and raising a flag if the lower bound for the displacements of the logging tool exceeds a selected threshold.
3c Other aspects and advantages of the invention will be apparent from the following description and the appended claims.
BRIEF DESCRIPTION OF THE DRAWINGS
Figure 1 shows a logging tool suspended in a borehole.
Figure 2 is a cross section of a logging tool suspended in a borehole according to one embodiment of the invention.
Figure 3 depicts a horizontal cross section of the logging tool shown in Figure 2.
3d PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) Figure 4 is a flow chart illustrating a method for determining the displacements of a logging tool according to one embodiment of the invention.
Figure 5 is a flow chart illustrating a method for determining the displacements of a logging tool according to another embodiment of the invention.
DETAILED DESCRIPTION OF THE INVENTION
Embodiments of the invention provide a method for determining displacements of a logging tool during a measurement interval along three orthogonal axes of the logging tool. In general, an accelerometer is used to measure the accelerations of the logging tool along the three orthogonal axes of the logging tool during the measurement interval. The accelerations acquired by the accelerometer, as will be further explained below, have a gravitational portion that is due to gravitational forces acting on the test-mass of the accelerometer and a kinetic portion that is due to the net force acting on the logging tool.
The displacements of the logging tool are determined from the estimated kinetic portion of the accelerations.
The displacements of the logging tool may be used to assess the quality of the measurements made by the logging tool. For example, pulse-echo nuclear magnetic resonance (NN>R) measurements are time-lapse measurements. For the measurement to be accurate, the sensitive zone of the NMR logging tool needs to substantially overlap with itself through out the measurement duration. Thus, accuracy of NMR
logging tools are sensitive to the displacement of the tool during the measurement interval.
By determining the displacements of the 'logging tool during a measurement interval, the validity of the measurements made can be verified. Of course, the invention is not limited to NMR logging tools, but is generally applicable to any logging tool that makes measurements that are sensitive to tool motion.
Various embodiments of the invention will now be discussed with reference to the accompanying figures. In order to fully understand the invention, it is helpful to consider a specific configuration of a logging tool. However, it should be clear that the invention is not limited to the specific configuration of the logging tool discussed herein. Figure 2
4 PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001 ) shows a borehole 14 which traverses a subterranean formation (or formations) 16. A
logging tool 18 is suspended in the borehole 14 on the end of a cable 20. The logging tool 18 includes a sonde 22 which measures characteristics of the formation 16 using NMR principles. An electronics cartridge 24 is mounted on the sonde 22. The electronics cartridge 24 includes a pulse generator 26 and may also include a memory 28 for storing data. In one embodiment, the sonde 22 includes a permanent magnet 30 which produces a static magnetic field Bo and an antenna 32 which produces an oscillating magnetic field B 1. The permanent magnet 30 circumscribes a protective sleeve 33. The sleeve provides a conduit for receiving electrical conductors 35 (shown in Figure 3) which transmit signals to the electronic cartridge 24. In one embodiment, the antenna 32 includes a fernte core 34 on which radio-frequency (RF) coils 36, 38, 40 are mounted.
The RF coil 38 has a variable resonant frequency, or receives a variable frequency RF
power, which may be adjusted to select the depth of investigation of the logging tool 18.
The RF coils 36, 38, 40 generate the oscillating magnetic field B1 in response to signals from the pulse generator 26. The pulse generators 26 may be controlled, for example, to generate NMR detection sequences such as a Carr-Puzcell-Meiboom-Gill (CPMG) sequence (not shown). The NMR detection sequence may be applied, for example, to determine the T2 spin-spin relaxation times of hydrogen nuclei in the formation 16. The static magnetic field Bo produced by the permanent magnet 30 and the oscillating magnetic field B1 produced by the antenna 32 create a resonance volume 42 in which the characteristics of the formation 16 can be investigated. In operation, the pulse generator 26 is controlled to produce a desired NMR detection sequence. The spin echo signals from the resonance volume 42 are received by the RF coils 36, 38, 40.
In one embodiment, the spin echo signals are stored in the memory 28 and later transmitted uphole. The spin echo signals may be transmitted uphole via telemetry, in which case, one or more receivers (not shown) will be provided to receive the signals. The spin echo signals may be amplified by amplifiers (not shown) and stored for further processing by a computer 43. For example, the spin echo signals may be analyzed to produce a distribution of T2 times, and the properties of the formation 16 may be obtained from this
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001 ) shows a borehole 14 which traverses a subterranean formation (or formations) 16. A
logging tool 18 is suspended in the borehole 14 on the end of a cable 20. The logging tool 18 includes a sonde 22 which measures characteristics of the formation 16 using NMR principles. An electronics cartridge 24 is mounted on the sonde 22. The electronics cartridge 24 includes a pulse generator 26 and may also include a memory 28 for storing data. In one embodiment, the sonde 22 includes a permanent magnet 30 which produces a static magnetic field Bo and an antenna 32 which produces an oscillating magnetic field B 1. The permanent magnet 30 circumscribes a protective sleeve 33. The sleeve provides a conduit for receiving electrical conductors 35 (shown in Figure 3) which transmit signals to the electronic cartridge 24. In one embodiment, the antenna 32 includes a fernte core 34 on which radio-frequency (RF) coils 36, 38, 40 are mounted.
The RF coil 38 has a variable resonant frequency, or receives a variable frequency RF
power, which may be adjusted to select the depth of investigation of the logging tool 18.
The RF coils 36, 38, 40 generate the oscillating magnetic field B1 in response to signals from the pulse generator 26. The pulse generators 26 may be controlled, for example, to generate NMR detection sequences such as a Carr-Puzcell-Meiboom-Gill (CPMG) sequence (not shown). The NMR detection sequence may be applied, for example, to determine the T2 spin-spin relaxation times of hydrogen nuclei in the formation 16. The static magnetic field Bo produced by the permanent magnet 30 and the oscillating magnetic field B1 produced by the antenna 32 create a resonance volume 42 in which the characteristics of the formation 16 can be investigated. In operation, the pulse generator 26 is controlled to produce a desired NMR detection sequence. The spin echo signals from the resonance volume 42 are received by the RF coils 36, 38, 40.
In one embodiment, the spin echo signals are stored in the memory 28 and later transmitted uphole. The spin echo signals may be transmitted uphole via telemetry, in which case, one or more receivers (not shown) will be provided to receive the signals. The spin echo signals may be amplified by amplifiers (not shown) and stored for further processing by a computer 43. For example, the spin echo signals may be analyzed to produce a distribution of T2 times, and the properties of the formation 16 may be obtained from this
5 PATENT
ATTORNEY DOCKET NO.: 20.291 S (R&0 REF. 09428/01 1001 ) distribution.
As shown in Figure 3, the resonance volumes 42 are typically shaped like a thin sheet with a thickness on the order of 1 millimeter. A particular resonance volume 42 is excited depending on the frequency of operation. Thus, if the logging tool 18 moves 1 millimeter or more in the radial direction, the T2 spin-spin relaxation times may be corrupted. Other NMR measurements, such as Tl relaxation time measurements, may also be compromised if the logging tool 18 accelerates in a direction along the longitudinal axis of the borehole 14 during a measurement interval. Thus, as shown in Figure 2, an accelerometer 44 is provided to sense the motion of the logging tool 18 during a logging operation. In one embodiment, the accelerometer 44 is mounted in the electronics cartridge 24, but may be mounted elsewhere as long as it is positioned as close as possible to the sonde 22 or the part of the logging tool 18 that is most sensitive to motion. The measurements made by the accelerometer 44 may be transmitted uphole via telemetry and processed, for example, by the computer 43.
For discussion purposes, a Cartesian coordinate system is fixed on the logging tool 18. The coordinate system has three mutually perpendicular axes, including radial (R), tangential (T), and axial (A) axes. The positive axial direction points up along the axis of the borehole 14, and the positive radial direction points into the formation 16.
The tangential axis is perpendicular to both the radial and axial axis and tangent to the wall of the borehole 14 where the logging tool 18 contacts the wall. The logging tool 18 is moved along the axis of the borehole 14 to make measurements. The accelerometer 44 includes, for example, three uniaxial sensors, each of which has a sensitive axis aligned with one of the axes of the logging tool 18. The accelerometer 44 measures instantaneous acceleration of the logging tool 18 along the radial, tangential, and axial directions as the logging tool 18 makes measurements.
When the logging tool 18 is at rest or moving at a constant velocity in the earth's gravitational field, the accelerometer 44 measures the radial component (gR), the tangential component (gT), and the axial component (gA) of the acceleration due to gravity (g = 981 cm/s2). The components of the acceleration due to gravity (g) are
ATTORNEY DOCKET NO.: 20.291 S (R&0 REF. 09428/01 1001 ) distribution.
As shown in Figure 3, the resonance volumes 42 are typically shaped like a thin sheet with a thickness on the order of 1 millimeter. A particular resonance volume 42 is excited depending on the frequency of operation. Thus, if the logging tool 18 moves 1 millimeter or more in the radial direction, the T2 spin-spin relaxation times may be corrupted. Other NMR measurements, such as Tl relaxation time measurements, may also be compromised if the logging tool 18 accelerates in a direction along the longitudinal axis of the borehole 14 during a measurement interval. Thus, as shown in Figure 2, an accelerometer 44 is provided to sense the motion of the logging tool 18 during a logging operation. In one embodiment, the accelerometer 44 is mounted in the electronics cartridge 24, but may be mounted elsewhere as long as it is positioned as close as possible to the sonde 22 or the part of the logging tool 18 that is most sensitive to motion. The measurements made by the accelerometer 44 may be transmitted uphole via telemetry and processed, for example, by the computer 43.
For discussion purposes, a Cartesian coordinate system is fixed on the logging tool 18. The coordinate system has three mutually perpendicular axes, including radial (R), tangential (T), and axial (A) axes. The positive axial direction points up along the axis of the borehole 14, and the positive radial direction points into the formation 16.
The tangential axis is perpendicular to both the radial and axial axis and tangent to the wall of the borehole 14 where the logging tool 18 contacts the wall. The logging tool 18 is moved along the axis of the borehole 14 to make measurements. The accelerometer 44 includes, for example, three uniaxial sensors, each of which has a sensitive axis aligned with one of the axes of the logging tool 18. The accelerometer 44 measures instantaneous acceleration of the logging tool 18 along the radial, tangential, and axial directions as the logging tool 18 makes measurements.
When the logging tool 18 is at rest or moving at a constant velocity in the earth's gravitational field, the accelerometer 44 measures the radial component (gR), the tangential component (gT), and the axial component (gA) of the acceleration due to gravity (g = 981 cm/s2). The components of the acceleration due to gravity (g) are
6 PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) referred to herein as "gravitational accelerations." These gravitational accelerations do not result in displacements of the logging tool 18 because the gravitational force on the logging tool 18 is balanced by the time average of the tension in the cable 20 and the friction with the formation 16 and the fluid in the borehole 14.
During a logging operation, however, the variable stretch in the cable 20 and the rough surface of the wall of the borehole 14 can exert fluctuating forces on the logging tool 18. The fluctuations in the net force acting on the logging tool 18 causes the logging tool 18 to accelerate and decelerate. This acceleration is different from the acceleration due to gravity and is called "kinetic acceleration" because it results in displacements of the logging tool 18. The kinetic acceleration is equal to the second time-derivative of the position of the logging tool 18 measured with respect to an inertial reference. The kinetic acceleration has a radial component zR , tangential component xT , and an axial component XA . Following standard conventions, dots above variables denote time-derivatives. The accelerometer 44 also measures the kinetic accelerations along the three axes of the logging tool 18. The total acceleration measured along the radial, tangential, and axial axes is then the sum of the gravitational and the kinetic accelerations.
The three-axis gravitational acceleration provides information on the orientation of the logging tool 18 with respect to the set of fixed axes XYZ. This information can be used to determine the deviation of the borehole 14 and the relative bearing of the logging tool 18 in the borehole 14. The kinetic acceleration, on the other hand, can be used to determine the displacements the logging tool 18. If the orientation of the logging tool 18 does not change during the data acquisition period, the gravitational accelerations along each axis of the logging tool 18 will remain constant. The kinetic accelerations of the logging tool 18 can then be determined by subtracting a constant from the acceleration data. In reality, however, the orientation of the logging tool 18 is not constant, but is generally slowly varying. Thus, a method for determining the gravitational accelerations of the logging tool 18 is needed. Embodiments of the invention provide a method for estimating the gravitational accelerations and removing the gravitational accelerations
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) referred to herein as "gravitational accelerations." These gravitational accelerations do not result in displacements of the logging tool 18 because the gravitational force on the logging tool 18 is balanced by the time average of the tension in the cable 20 and the friction with the formation 16 and the fluid in the borehole 14.
During a logging operation, however, the variable stretch in the cable 20 and the rough surface of the wall of the borehole 14 can exert fluctuating forces on the logging tool 18. The fluctuations in the net force acting on the logging tool 18 causes the logging tool 18 to accelerate and decelerate. This acceleration is different from the acceleration due to gravity and is called "kinetic acceleration" because it results in displacements of the logging tool 18. The kinetic acceleration is equal to the second time-derivative of the position of the logging tool 18 measured with respect to an inertial reference. The kinetic acceleration has a radial component zR , tangential component xT , and an axial component XA . Following standard conventions, dots above variables denote time-derivatives. The accelerometer 44 also measures the kinetic accelerations along the three axes of the logging tool 18. The total acceleration measured along the radial, tangential, and axial axes is then the sum of the gravitational and the kinetic accelerations.
The three-axis gravitational acceleration provides information on the orientation of the logging tool 18 with respect to the set of fixed axes XYZ. This information can be used to determine the deviation of the borehole 14 and the relative bearing of the logging tool 18 in the borehole 14. The kinetic acceleration, on the other hand, can be used to determine the displacements the logging tool 18. If the orientation of the logging tool 18 does not change during the data acquisition period, the gravitational accelerations along each axis of the logging tool 18 will remain constant. The kinetic accelerations of the logging tool 18 can then be determined by subtracting a constant from the acceleration data. In reality, however, the orientation of the logging tool 18 is not constant, but is generally slowly varying. Thus, a method for determining the gravitational accelerations of the logging tool 18 is needed. Embodiments of the invention provide a method for estimating the gravitational accelerations and removing the gravitational accelerations
7 PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) from the acceleration data so that the displacements of the logging tool 18 can be estimated.
The problem addressed by the invention is akin to a physicist estimating the distance traveled by the elevator in which she is riding. The physicist is reading the apparent weight of an apple of known mass on a balance inside the elevator. As the elevator accelerates going up or decelerates going down, the balance reading increases.
As the elevator decelerates going up or accelerates going down, the balance reading decreases. The physicist could calculate the distance traveled by the elevator if she were not handicapped by two factors: (1) the building has an unknown tilt and (2) she is distracted at the beginning so she does not know the balance reading at rest or the initial velocity of the elevator when she starts her measurements. The physicist can determine the changes in acceleration which tells her the position of the elevator up to an arbitrary quadratic polynomial of time. Given this incomplete information, the physicist can only put a lower bound on how much the elevator might have traveled since she started her measurements.
For discussion purposes, let a(t) be the acceleration measured along any one of the axes of the logging tool 18 at time t >_ t1 , where t, is the time that the data acquisition begins. The acceleration measured by the accelerometer 44 includes the kinetic accelerations and the gravitational accelerations of the logging tool 18. That is, a(t) - x(t) - gx (1) where x(t) is the kinetic acceleration of the logging tool 18 due to all forces acting on the logging tool and gx is the component of the acceleration due to gravity, i.e., gravitational acceleration, in the x-direction, i.e., along one of the axes of the logging tool 18. The tool position x(t) along one of the axes of the logging tool 18, generally denoted as the x-direction, at time t is given by the following expression:
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) from the acceleration data so that the displacements of the logging tool 18 can be estimated.
The problem addressed by the invention is akin to a physicist estimating the distance traveled by the elevator in which she is riding. The physicist is reading the apparent weight of an apple of known mass on a balance inside the elevator. As the elevator accelerates going up or decelerates going down, the balance reading increases.
As the elevator decelerates going up or accelerates going down, the balance reading decreases. The physicist could calculate the distance traveled by the elevator if she were not handicapped by two factors: (1) the building has an unknown tilt and (2) she is distracted at the beginning so she does not know the balance reading at rest or the initial velocity of the elevator when she starts her measurements. The physicist can determine the changes in acceleration which tells her the position of the elevator up to an arbitrary quadratic polynomial of time. Given this incomplete information, the physicist can only put a lower bound on how much the elevator might have traveled since she started her measurements.
For discussion purposes, let a(t) be the acceleration measured along any one of the axes of the logging tool 18 at time t >_ t1 , where t, is the time that the data acquisition begins. The acceleration measured by the accelerometer 44 includes the kinetic accelerations and the gravitational accelerations of the logging tool 18. That is, a(t) - x(t) - gx (1) where x(t) is the kinetic acceleration of the logging tool 18 due to all forces acting on the logging tool and gx is the component of the acceleration due to gravity, i.e., gravitational acceleration, in the x-direction, i.e., along one of the axes of the logging tool 18. The tool position x(t) along one of the axes of the logging tool 18, generally denoted as the x-direction, at time t is given by the following expression:
8 PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) r l~ , x(t) - xi + (t - t~ )xi + j~~ f ~zx(t2 ) r, r, where x, is the initial position and z, is the initial velocity of the logging tool 18 at time t1. When equation (1) is substituted into equation (2), the following expression for the tool displacement is obtained:
x(t)-x~=(t-t~)x~+ Zx(t-t~)2+ jdtlrfdtZa(t2) (3) r, r, where it is assumed that gx is approximately constant over the data acquisition period.
S Because gx depends on the orientation of the axes of the logging tool 18 relative to the set of fixed reference axes XYZ, this assumption is equivalent to assuming that the orientation of the logging tool 18 slowly varies with time. This assumption is sensible for short data acquisition periods, which are typically on the order of 0.6 seconds or shorter for the CPMG measurement sequence in NMR logging.
Two quantities in equation (3 ), gx , the gravitational acceleration, and z1 , the initial velocity of the logging tool 18, are unknown. Because, the parameter of interest is the magnitude of the displacement of the logging tool 18 from an initial position, and not the actual position of the logging tool 18 in the borehole 14, the knowledge of x, is not necessary. The displacement x(t) - x, is, therefore, renamed as x(t) from here on. In other words, the initial position is arbitrarily chosen as the origin of the coordinate system. The notation used for the tool displacement from here on emphasizes its functional dependence on the initial velocity z1 and the gravitational acceleration gx , as shown in equation (4) below.
x(t~ x1 ~ 8X ) = (r - t~ )xi + Zx (t - t~ ) 2 + f ~, r f dF2 a(12 ) (q.) r, r,
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) r l~ , x(t) - xi + (t - t~ )xi + j~~ f ~zx(t2 ) r, r, where x, is the initial position and z, is the initial velocity of the logging tool 18 at time t1. When equation (1) is substituted into equation (2), the following expression for the tool displacement is obtained:
x(t)-x~=(t-t~)x~+ Zx(t-t~)2+ jdtlrfdtZa(t2) (3) r, r, where it is assumed that gx is approximately constant over the data acquisition period.
S Because gx depends on the orientation of the axes of the logging tool 18 relative to the set of fixed reference axes XYZ, this assumption is equivalent to assuming that the orientation of the logging tool 18 slowly varies with time. This assumption is sensible for short data acquisition periods, which are typically on the order of 0.6 seconds or shorter for the CPMG measurement sequence in NMR logging.
Two quantities in equation (3 ), gx , the gravitational acceleration, and z1 , the initial velocity of the logging tool 18, are unknown. Because, the parameter of interest is the magnitude of the displacement of the logging tool 18 from an initial position, and not the actual position of the logging tool 18 in the borehole 14, the knowledge of x, is not necessary. The displacement x(t) - x, is, therefore, renamed as x(t) from here on. In other words, the initial position is arbitrarily chosen as the origin of the coordinate system. The notation used for the tool displacement from here on emphasizes its functional dependence on the initial velocity z1 and the gravitational acceleration gx , as shown in equation (4) below.
x(t~ x1 ~ 8X ) = (r - t~ )xi + Zx (t - t~ ) 2 + f ~, r f dF2 a(12 ) (q.) r, r,
9 PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&O REF. 09428/011001 ) In practice, the output of the accelerometers 44 are not continuously recorded in time, but a finite number of samples are acquired with a constant time interval ~.
Assuming that the accelerometer acquires ns samples in the x-direction, i.e., along one of the axes of the logging tool 18, then an = xn - gx for n =1,. .., ns (5) where an is the acceleration measured in the x-direction at the time t = n0.
The term xn is the kinetic acceleration of the logging tool 18, and gx is the component of the gravitational acceleration in the x-direction. A single integration of the acceleration data gives the set of velocities of the logging tool 18. The acceleration data can be integrated using a variety of numerical methods. One suitable method is the trapezoid rule for numerical integration. When the trapezoid rule is applied to equation (5), the following expression is obtained:
xn+, = xn + 2 (xn+I + xn ) for n =1, . . . , ns -1 (6) Equation (6) gives the velocity at the (n+1)'~' time step in terms of the velocity at the previous time step plus the change in the velocity due to the acceleration.
Repeated application of the recursion relation (6) and use of equation (5) leads, after n time steps, to:
zn = x, + (n -1)gx0 + ~ (al + 2a2 + 2a3 + ~ ~ ~ + 2an_, + an ) for n = 2,3, .
. . ns (7) Using the trapezoid rule a second time to integrate equation (6), the following expression is obtained:
xn+1 = xn + 2 (xn+1 + xn ) for n = 1, . . . , Yls -1 (g) PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) Equations (7) and (8) lead to:
xn (x1 ~ 8X ) _ (n -1)xl ~ + (n 1 ~z OZ gX +
+02C2 4 3a +(n-2)a +(n-3)a +...+3a +2a +a + 1 a (9) 1 2 3 n-3 n-2 n-I 4 n n = 2,3,...,nS
Equation (9) shows the explicit functional dependence of the displacement on the unknown initial velocity x1 and gravitational acceleration gx .
Figure 4 illustrates a method for estimating a lower bound on displacements of the logging tool 18 given that x, and gx are unknown. The method starts by acquiring ns acceleration samples during a measurement interval of the logging tool 18 (shown at 46).
The next step is to determine the particular values of x, and gx that minimize the estimated tool displacement in the following sense:
gX (x1 ) = arg min ~ xn 2 (z, , 8x ) ( 10 8r n=2 x1 = argmin{max ( xn(xl,gx(xl)) I} (11) n i The notation " arg min f ( p) " denotes the value of the parameter p that minimizes the P
expression f ( p) . The gravitational acceleration gx (z1 ) is estimated by minimizing the sum of squares of the displacement time-series. This value is readily calculated by setting the derivative of the sum of squares with respect to gx to zero:
. . CA 02313282 2000-06-30 PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&O REF. 09428/011001 ) ~n, 2L xn (x1 ~0)(n _ 1) 2 DZ
gs(x') = n' n (12) ~~n -1~) n=' The minimization in equation (11) with respect to the initial velocity x, is done by searching for the minimum through a set of user-supplied initial velocities {z,~'~,...,x,~m~}. An i''h initial velocity from the set of user-supplied initial velocities is first obtained (shown at 50). For each i'h initial velocity, an estimate gX'~
is next calculated using equation (12) above (shown at 52). For each i'h initial velocity, there will be a time-series of ns displacements corresponding to the n$ acceleration samples and an estimated value of the gravitational acceleration. In step 54, the maximum of the ns-long displacement time-series is selected. The steps 46-54 are repeated until all the displacements for the set of user-supplied initial velocities have been computed. In step 56, the minimum of the maximum displacements computed in step 54 is selected as the lower bound for the displacement of the logging tool 18 during data acquisition. The initial velocity corresponding to this lower bound is the solution to equation (11). The lower bound for the displacement of the logging tool 18 can be used to assess the measurements made by the logging tool 18. For example, the condition that the lower bound for the peak displacement of the logging tool 18 exceeds a certain fraction of the thickness of the resonance volume 42 can be used to flag the NMR measurement as invalid (shown at 57).
In an alternate embodiment, gx is assumed to be approximately constant during the data acquisition period. In this case, the mean value of the acceleration samples acquired in step 45 may provide another estimate of gz . This mean value gx,mean may replace the estimate gx (x~ ) calculated in step 52.
Figure 5 illustrates an alternative method for estimating a lower bound for the displacement of the logging tool 18. Because the tool displacement is known up to an . . CA 02313282 2000-06-30 PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) arbitrary quadratic polynomial of time, if any quadratic polynomial of time from the displacement time-series is subtracted, the result will also be a displacement time-series that is consistent with the measured acceleration time-series. There is a unique quadratic polynomial that will minimize the sum of squares of the resulting time-series.
This is the well-defined, unique lower bound for the tool displacement in the least-squares sense. In this method, the motion of the logging tool 18 is represented by the following expression:
ate - x(t) (13) where x(t) is the acceleration of the logging tool 18 along any one of the tool axes, denoted by x, at time t. The derivative ~t2 is then replaced by a central-difference approximation, as shown in equation ( 14) below:
xn+~ - 2x" + x"_, ( 14) = xn where O is the time spacing between x"+1 and x" . For n = 1 to ns, where ns is the number of acceleration samples acquired along any one of the tool axes with sample spacing 0, a system of ns equations can be written using equation ( 14) above. The system of equations can be expressed in matrix form as follows:
x~ xaz - xo xz xxOz T . . (15) z x», x", 0 - xn,+~
PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) where -2 1 0 ~
1 -2 1 ~
T= . 1 -2 ~
When equation (5) is substituted into equation (15), the following expression is obtained:
x 2_ 1 ~al + gr )~ x0 x 2 ~a2 +gx)e T . . (16) Z
x», ~a", +gx)~ -x»,+~
In this notation, time-series are represented by column vectors. The solution to the matrix equation (16) above is a tool displacement vector x = {x" x2,..., x"s }
, where xo and x"~+, are the boundary values of the displacements of the logging tool 18.
In the following discussion, it is convenient to use Dirac's notation for ket and bra (see Merzbacher, E., Quantum Mechanics, John Wiley & Sons, 1961). Let ~x) represent the displacement vector and let ~ a) represent the vector on the right-hand side of equation ( 16). Then equation ( 16) can be rewritten as follows:
TI x~=~a~ (17) The solution to equation (17) is obtained by inverting the matrix T and multiplying the vector ~ a) by the inverted matrix T:
(x~=T_ ya~ (18) As shown in equation (15), the matrix T is in tridiagonal form and can be readily inverted. See, for example, Ralston, A. and Wilf, H.S., Editors, Mathematical Methods PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001 ) for Digital Computers, Vol. 2, John Wiley & Sons, 1967. It should be noted that the acceleration data provides the values for the elements of the vector ~ a) .
The boundary conditions xo = x" +~ = 0 are used in computing the vector ~ a) . The result does not depend on the choice of the boundary conditions as the operation of subtracting a quadratic polynomial of time undoes the effect of the boundary values.
The method illustrated in Figure 5 starts by acquiring ns acceleration samples during a measurement interval of the logging tool 18 (shown at 5.8). The next step (shown at 60) involves solving for the displacement vector ~x) using equation (18). The method estimates the displacements of the logging tool 18 by removing the projections of ~ x) onto orthogonal vectors that represent constant, linear, and quadratic time dependencies from ~ x) . Consider a subspace consisting of three linearly independent vectors ~ uo ) , ~ u, ) , and ~ uz ) in an ns-dimensional vector space, where:
3 ' ~uz)= 9 (19) n.z 1 ns These vectors are the samples of elementary polynomials, e.g., 1, t, t~ which are linearly independent. Their linear combinations span samples of any quadratic polynomial of time. Orthonomal vectors can be constructed from the vectors ~ uo ) , ( u~ ) , and ~ uz ) by the Gramm-Schmitt orthogonalization procedure Iuo)=luo) (20a) (20b) ~uo~uo) PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) ~u~ I u2 ~I u~ ) ~uo ~ u2 uo ) (20c) ~uz)=~uz)- ~u,Iu,) - ~tl Itl >
The linear and quadratic time dependencies are removed from the displacement vector ~x~ computed in step 60 by subtracting the projection of the displacement vector ~x~
along the orthogonal vectors ~ u-, ) and ~ u-Z ~ in equations (20b) and (20c), shown at 62.
That is, ~x~=w_~~uflw u~) ~ul lure .(21) where w = I x~ .
The minimum displacements during the data acquisition period are obtained by subtracting the initial position from each element in the displacement vector (shown at 64). That is, x) min =I x > -x. I uo ) (22) where x, is the first entry in ( x) . The operation in step 62 is equivalent to removing the constant dependencies from the displacement vector. The norms of the vectors ~
u, ) are needed in equations (21 ) and (22) and can be computed by straightforward algebra using well known summation formulae. See, for example, Jolley, L.B.W., Summation of Series, Dover Publications, Inc., 1961. The norms of the vectors ~u,) are:
(23a) ny (nJ + 1)(n, -1) (23b) u, u, = 12 PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001 ) ~u2 I u2 ~ _ na ~n~ 2 180n84 4> (23 c) The norms shown in equations (23a) through (23c) do not change and can be calculated prior to starting the process of acquiring the acceleration samples and estimating a lower bound on the displacement of the logging tool 18 (shown at 66). As in the previous method, if the lower bound determined in step 64 exceeds a predetermined threshold, a flag can be raised (shown at 68). The algorithm described in Figure 5 is mathematically equivalent to minimizing the sum of squares of equation (9) with respect to x, and gx .
The lower bounds computed by the methods described in Figures 4 and 5 are comparable.
In operation, the logging tool 18 is moved along the borehole 14 to make measurements. The sonde 22 makes NMR measurements by magnetically tipping the nuclear spins of protons in the formation with pulses of the oscillating magnetic field B1, and then detecting the precession of the tipped particles in the resonance volume 42. The accelerometer 44 measures the acceleration of the logging tool 18 during the NMR
measurements. The acceleration signals from the accelerometer 44 may be transmitted to the surface in real time or stored in a memory and later transmitted to the surface. At the surface, the acceleration signals may be amplified and then processed.- Using the methods described above, the computer 43 computes the true displacements of the logging tool 18 during data acquisition along the three orthogonal axes of the logging tool 18. These true displacements can then be used to isolate portions of the NMR
log that may be distorted by motions of the logging tool 18. For example, for T2 relaxation-time measurements, the true displacements along the radial axis of the logging tool 18 can be used to identify invalid data in the NMR log. For T1 relaxation-time measurements, the true displacements along the axial axis of the logging tool 18 is used to assess the quality of the log. It should be clear that the methods described above are not limited to the specific configuration of the logging tool 18 shown in Figures 2 and 3, but can be used to determine true displacements of any logging tool in general, regardless of whether the logging tool is used alone or is included in other assemblies, e.g., a drill string.
PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001 ) While the invention has been described with respect to a limited number of embodiments, those skilled in the art will appreciate that other embodiments can be devised which do not depart from the scope of the invention as disclosed herein.
Accordingly, the scope of the invention should be limited only by the attached claims.
ATTORNEY DOCKET NO.: 20.2715 (R&O REF. 09428/011001 ) In practice, the output of the accelerometers 44 are not continuously recorded in time, but a finite number of samples are acquired with a constant time interval ~.
Assuming that the accelerometer acquires ns samples in the x-direction, i.e., along one of the axes of the logging tool 18, then an = xn - gx for n =1,. .., ns (5) where an is the acceleration measured in the x-direction at the time t = n0.
The term xn is the kinetic acceleration of the logging tool 18, and gx is the component of the gravitational acceleration in the x-direction. A single integration of the acceleration data gives the set of velocities of the logging tool 18. The acceleration data can be integrated using a variety of numerical methods. One suitable method is the trapezoid rule for numerical integration. When the trapezoid rule is applied to equation (5), the following expression is obtained:
xn+, = xn + 2 (xn+I + xn ) for n =1, . . . , ns -1 (6) Equation (6) gives the velocity at the (n+1)'~' time step in terms of the velocity at the previous time step plus the change in the velocity due to the acceleration.
Repeated application of the recursion relation (6) and use of equation (5) leads, after n time steps, to:
zn = x, + (n -1)gx0 + ~ (al + 2a2 + 2a3 + ~ ~ ~ + 2an_, + an ) for n = 2,3, .
. . ns (7) Using the trapezoid rule a second time to integrate equation (6), the following expression is obtained:
xn+1 = xn + 2 (xn+1 + xn ) for n = 1, . . . , Yls -1 (g) PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) Equations (7) and (8) lead to:
xn (x1 ~ 8X ) _ (n -1)xl ~ + (n 1 ~z OZ gX +
+02C2 4 3a +(n-2)a +(n-3)a +...+3a +2a +a + 1 a (9) 1 2 3 n-3 n-2 n-I 4 n n = 2,3,...,nS
Equation (9) shows the explicit functional dependence of the displacement on the unknown initial velocity x1 and gravitational acceleration gx .
Figure 4 illustrates a method for estimating a lower bound on displacements of the logging tool 18 given that x, and gx are unknown. The method starts by acquiring ns acceleration samples during a measurement interval of the logging tool 18 (shown at 46).
The next step is to determine the particular values of x, and gx that minimize the estimated tool displacement in the following sense:
gX (x1 ) = arg min ~ xn 2 (z, , 8x ) ( 10 8r n=2 x1 = argmin{max ( xn(xl,gx(xl)) I} (11) n i The notation " arg min f ( p) " denotes the value of the parameter p that minimizes the P
expression f ( p) . The gravitational acceleration gx (z1 ) is estimated by minimizing the sum of squares of the displacement time-series. This value is readily calculated by setting the derivative of the sum of squares with respect to gx to zero:
. . CA 02313282 2000-06-30 PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&O REF. 09428/011001 ) ~n, 2L xn (x1 ~0)(n _ 1) 2 DZ
gs(x') = n' n (12) ~~n -1~) n=' The minimization in equation (11) with respect to the initial velocity x, is done by searching for the minimum through a set of user-supplied initial velocities {z,~'~,...,x,~m~}. An i''h initial velocity from the set of user-supplied initial velocities is first obtained (shown at 50). For each i'h initial velocity, an estimate gX'~
is next calculated using equation (12) above (shown at 52). For each i'h initial velocity, there will be a time-series of ns displacements corresponding to the n$ acceleration samples and an estimated value of the gravitational acceleration. In step 54, the maximum of the ns-long displacement time-series is selected. The steps 46-54 are repeated until all the displacements for the set of user-supplied initial velocities have been computed. In step 56, the minimum of the maximum displacements computed in step 54 is selected as the lower bound for the displacement of the logging tool 18 during data acquisition. The initial velocity corresponding to this lower bound is the solution to equation (11). The lower bound for the displacement of the logging tool 18 can be used to assess the measurements made by the logging tool 18. For example, the condition that the lower bound for the peak displacement of the logging tool 18 exceeds a certain fraction of the thickness of the resonance volume 42 can be used to flag the NMR measurement as invalid (shown at 57).
In an alternate embodiment, gx is assumed to be approximately constant during the data acquisition period. In this case, the mean value of the acceleration samples acquired in step 45 may provide another estimate of gz . This mean value gx,mean may replace the estimate gx (x~ ) calculated in step 52.
Figure 5 illustrates an alternative method for estimating a lower bound for the displacement of the logging tool 18. Because the tool displacement is known up to an . . CA 02313282 2000-06-30 PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) arbitrary quadratic polynomial of time, if any quadratic polynomial of time from the displacement time-series is subtracted, the result will also be a displacement time-series that is consistent with the measured acceleration time-series. There is a unique quadratic polynomial that will minimize the sum of squares of the resulting time-series.
This is the well-defined, unique lower bound for the tool displacement in the least-squares sense. In this method, the motion of the logging tool 18 is represented by the following expression:
ate - x(t) (13) where x(t) is the acceleration of the logging tool 18 along any one of the tool axes, denoted by x, at time t. The derivative ~t2 is then replaced by a central-difference approximation, as shown in equation ( 14) below:
xn+~ - 2x" + x"_, ( 14) = xn where O is the time spacing between x"+1 and x" . For n = 1 to ns, where ns is the number of acceleration samples acquired along any one of the tool axes with sample spacing 0, a system of ns equations can be written using equation ( 14) above. The system of equations can be expressed in matrix form as follows:
x~ xaz - xo xz xxOz T . . (15) z x», x", 0 - xn,+~
PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) where -2 1 0 ~
1 -2 1 ~
T= . 1 -2 ~
When equation (5) is substituted into equation (15), the following expression is obtained:
x 2_ 1 ~al + gr )~ x0 x 2 ~a2 +gx)e T . . (16) Z
x», ~a", +gx)~ -x»,+~
In this notation, time-series are represented by column vectors. The solution to the matrix equation (16) above is a tool displacement vector x = {x" x2,..., x"s }
, where xo and x"~+, are the boundary values of the displacements of the logging tool 18.
In the following discussion, it is convenient to use Dirac's notation for ket and bra (see Merzbacher, E., Quantum Mechanics, John Wiley & Sons, 1961). Let ~x) represent the displacement vector and let ~ a) represent the vector on the right-hand side of equation ( 16). Then equation ( 16) can be rewritten as follows:
TI x~=~a~ (17) The solution to equation (17) is obtained by inverting the matrix T and multiplying the vector ~ a) by the inverted matrix T:
(x~=T_ ya~ (18) As shown in equation (15), the matrix T is in tridiagonal form and can be readily inverted. See, for example, Ralston, A. and Wilf, H.S., Editors, Mathematical Methods PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001 ) for Digital Computers, Vol. 2, John Wiley & Sons, 1967. It should be noted that the acceleration data provides the values for the elements of the vector ~ a) .
The boundary conditions xo = x" +~ = 0 are used in computing the vector ~ a) . The result does not depend on the choice of the boundary conditions as the operation of subtracting a quadratic polynomial of time undoes the effect of the boundary values.
The method illustrated in Figure 5 starts by acquiring ns acceleration samples during a measurement interval of the logging tool 18 (shown at 5.8). The next step (shown at 60) involves solving for the displacement vector ~x) using equation (18). The method estimates the displacements of the logging tool 18 by removing the projections of ~ x) onto orthogonal vectors that represent constant, linear, and quadratic time dependencies from ~ x) . Consider a subspace consisting of three linearly independent vectors ~ uo ) , ~ u, ) , and ~ uz ) in an ns-dimensional vector space, where:
3 ' ~uz)= 9 (19) n.z 1 ns These vectors are the samples of elementary polynomials, e.g., 1, t, t~ which are linearly independent. Their linear combinations span samples of any quadratic polynomial of time. Orthonomal vectors can be constructed from the vectors ~ uo ) , ( u~ ) , and ~ uz ) by the Gramm-Schmitt orthogonalization procedure Iuo)=luo) (20a) (20b) ~uo~uo) PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001) ~u~ I u2 ~I u~ ) ~uo ~ u2 uo ) (20c) ~uz)=~uz)- ~u,Iu,) - ~tl Itl >
The linear and quadratic time dependencies are removed from the displacement vector ~x~ computed in step 60 by subtracting the projection of the displacement vector ~x~
along the orthogonal vectors ~ u-, ) and ~ u-Z ~ in equations (20b) and (20c), shown at 62.
That is, ~x~=w_~~uflw u~) ~ul lure .(21) where w = I x~ .
The minimum displacements during the data acquisition period are obtained by subtracting the initial position from each element in the displacement vector (shown at 64). That is, x) min =I x > -x. I uo ) (22) where x, is the first entry in ( x) . The operation in step 62 is equivalent to removing the constant dependencies from the displacement vector. The norms of the vectors ~
u, ) are needed in equations (21 ) and (22) and can be computed by straightforward algebra using well known summation formulae. See, for example, Jolley, L.B.W., Summation of Series, Dover Publications, Inc., 1961. The norms of the vectors ~u,) are:
(23a) ny (nJ + 1)(n, -1) (23b) u, u, = 12 PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001 ) ~u2 I u2 ~ _ na ~n~ 2 180n84 4> (23 c) The norms shown in equations (23a) through (23c) do not change and can be calculated prior to starting the process of acquiring the acceleration samples and estimating a lower bound on the displacement of the logging tool 18 (shown at 66). As in the previous method, if the lower bound determined in step 64 exceeds a predetermined threshold, a flag can be raised (shown at 68). The algorithm described in Figure 5 is mathematically equivalent to minimizing the sum of squares of equation (9) with respect to x, and gx .
The lower bounds computed by the methods described in Figures 4 and 5 are comparable.
In operation, the logging tool 18 is moved along the borehole 14 to make measurements. The sonde 22 makes NMR measurements by magnetically tipping the nuclear spins of protons in the formation with pulses of the oscillating magnetic field B1, and then detecting the precession of the tipped particles in the resonance volume 42. The accelerometer 44 measures the acceleration of the logging tool 18 during the NMR
measurements. The acceleration signals from the accelerometer 44 may be transmitted to the surface in real time or stored in a memory and later transmitted to the surface. At the surface, the acceleration signals may be amplified and then processed.- Using the methods described above, the computer 43 computes the true displacements of the logging tool 18 during data acquisition along the three orthogonal axes of the logging tool 18. These true displacements can then be used to isolate portions of the NMR
log that may be distorted by motions of the logging tool 18. For example, for T2 relaxation-time measurements, the true displacements along the radial axis of the logging tool 18 can be used to identify invalid data in the NMR log. For T1 relaxation-time measurements, the true displacements along the axial axis of the logging tool 18 is used to assess the quality of the log. It should be clear that the methods described above are not limited to the specific configuration of the logging tool 18 shown in Figures 2 and 3, but can be used to determine true displacements of any logging tool in general, regardless of whether the logging tool is used alone or is included in other assemblies, e.g., a drill string.
PATENT
ATTORNEY DOCKET NO.: 20.2715 (R&0 REF. 09428/011001 ) While the invention has been described with respect to a limited number of embodiments, those skilled in the art will appreciate that other embodiments can be devised which do not depart from the scope of the invention as disclosed herein.
Accordingly, the scope of the invention should be limited only by the attached claims.
Claims (11)
1. A method for determining displacements of a logging tool during a measurement interval of the logging tool in a borehole, the method comprising:
obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval;
double integrating the set of accelerometer signals to obtain corresponding displacements of the logging tool as a function of an initial velocity of the logging tool and a gravitational acceleration, wherein the initial velocity of the logging tool and the gravitational acceleration are unknown;
determining a set of initial velocities for the logging tool;
for each initial velocity, estimating the gravitational acceleration, calculating the displacements of the logging tool using the initial velocity and the estimated gravitational acceleration, and determining a maximum of the calculated displacements; and setting a lower bound on the displacements of the logging tool to a value which is a minimum value of a set of the maximum calculated displacements for the set of initial velocities.
obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval;
double integrating the set of accelerometer signals to obtain corresponding displacements of the logging tool as a function of an initial velocity of the logging tool and a gravitational acceleration, wherein the initial velocity of the logging tool and the gravitational acceleration are unknown;
determining a set of initial velocities for the logging tool;
for each initial velocity, estimating the gravitational acceleration, calculating the displacements of the logging tool using the initial velocity and the estimated gravitational acceleration, and determining a maximum of the calculated displacements; and setting a lower bound on the displacements of the logging tool to a value which is a minimum value of a set of the maximum calculated displacements for the set of initial velocities.
2. The method of claim 1, wherein estimating the gravitational acceleration comprises minimizing a sum of a square of the displacements with respect to the unknown gravitational acceleration.
3. The method of claim 1, wherein estimating the gravitational acceleration includes averaging the accelerometer signals.
4. A method for improving the quality of measurements made by a logging tool during a measurement interval in a borehole, the method comprising:
obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval;
double integrating the set of accelerometer signals to obtain corresponding displacements of the logging tool as a function of an initial velocity of the logging tool and a gravitational acceleration, wherein the initial velocity of the logging tool and the gravitational acceleration are unknown;
determining a set of initial velocities for the logging tool;
for each initial velocity, estimating the gravitational acceleration, calculating the displacements of the logging tool using the initial velocity and the estimated gravitational acceleration, and determining a maximum of the calculated displacements;
estimating a lower bound for the displacements of the logging tool by selecting a value which is a minimum value of a set of the maximum calculated displacements for the initial velocities; and raising a flag if the lower bound for the displacements of the logging tool exceeds a selected threshold.
obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval;
double integrating the set of accelerometer signals to obtain corresponding displacements of the logging tool as a function of an initial velocity of the logging tool and a gravitational acceleration, wherein the initial velocity of the logging tool and the gravitational acceleration are unknown;
determining a set of initial velocities for the logging tool;
for each initial velocity, estimating the gravitational acceleration, calculating the displacements of the logging tool using the initial velocity and the estimated gravitational acceleration, and determining a maximum of the calculated displacements;
estimating a lower bound for the displacements of the logging tool by selecting a value which is a minimum value of a set of the maximum calculated displacements for the initial velocities; and raising a flag if the lower bound for the displacements of the logging tool exceeds a selected threshold.
5. A method for logging a well, comprising:
moving a logging tool along a borehole to make measurements in a formation surrounding the borehole;
recording the measurements made by the logging tool;
measuring a set of accelerometer signals of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval;
double integrating the set of accelerometer signals to obtain corresponding displacements of the logging tool as a function of an initial velocity of the logging tool and a gravitational acceleration, wherein the initial velocity of the logging tool and the gravitational acceleration are unknown;
determining a set of initial velocities for the logging tool;
for each initial velocity, estimating the gravitational acceleration, calculating the displacements of the logging tool using the initial velocity and the estimated gravitational acceleration, and determining a maximum of the calculated displacements;
estimating a lower bound for the displacements of the logging tool by selecting a value which is a minimum value of a set of the maximum calculated displacements for the initial velocities; and raising a flag if the lower bound for the displacements of the logging tool exceeds a selected threshold.
moving a logging tool along a borehole to make measurements in a formation surrounding the borehole;
recording the measurements made by the logging tool;
measuring a set of accelerometer signals of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval;
double integrating the set of accelerometer signals to obtain corresponding displacements of the logging tool as a function of an initial velocity of the logging tool and a gravitational acceleration, wherein the initial velocity of the logging tool and the gravitational acceleration are unknown;
determining a set of initial velocities for the logging tool;
for each initial velocity, estimating the gravitational acceleration, calculating the displacements of the logging tool using the initial velocity and the estimated gravitational acceleration, and determining a maximum of the calculated displacements;
estimating a lower bound for the displacements of the logging tool by selecting a value which is a minimum value of a set of the maximum calculated displacements for the initial velocities; and raising a flag if the lower bound for the displacements of the logging tool exceeds a selected threshold.
6. A method for determining displacements of a logging tool during a measurement interval of the logging tool in a borehole, the method comprising:
obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval;
calculating a tool displacement as a time-series from the accelerometer signals;
constructing a unique quadratic polynomial of time from the displacement time-series; and subtracting the unique quadratic polynomial from the displacement time-series; and setting a lower bound to the maximum of the remainder of the displacement time-series.
obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval;
calculating a tool displacement as a time-series from the accelerometer signals;
constructing a unique quadratic polynomial of time from the displacement time-series; and subtracting the unique quadratic polynomial from the displacement time-series; and setting a lower bound to the maximum of the remainder of the displacement time-series.
7. The method of claim 6, wherein calculating a tool displacement as a time-series from the accelerometer signals includes setting a second time-derivative of a position of the logging tool to the acceleration of the logging tool.
8. The method of claim 7, further comprising replacing the second time-derivative of the position of the logging tool with a central-difference approximation.
9. The method of claim 8, further comprising constructing a system of equations from the central-difference approximation and the acceleration of the logging tool and solving a system of equations to obtain the tool displacement.
10. The method of claim 7, wherein constructing a unique quadratic polynomial of time from the displacement-time series comprises combining elementary polynomials.
11. A method for improving the quality of measurements made by a logging tool during a measurement interval in a borehole, the method comprising:
obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval;
calculating a tool displacement as a time-series from the accelerometer signals;
constructing a unique quadratic polynomial of time from the displacement time-series;
subtracting the unique quadratic polynomial from the displacement time-series; and setting a lower bound to the maximum of the remainder of the displacement time-series; and raising a flag if the lower bound for the displacements of the logging tool exceeds a selected threshold.
obtaining a set of accelerometer signals corresponding to accelerations of the logging tool along each of three orthogonal axes of the logging tool during the measurement interval;
calculating a tool displacement as a time-series from the accelerometer signals;
constructing a unique quadratic polynomial of time from the displacement time-series;
subtracting the unique quadratic polynomial from the displacement time-series; and setting a lower bound to the maximum of the remainder of the displacement time-series; and raising a flag if the lower bound for the displacements of the logging tool exceeds a selected threshold.
Applications Claiming Priority (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US14339399P | 1999-07-12 | 1999-07-12 | |
| US60/143,393 | 1999-07-12 | ||
| US09/598,629 | 2000-06-21 | ||
| US09/598,629 US6459992B1 (en) | 1999-07-12 | 2000-06-21 | Method and apparatus for determining logging tool displacements |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CA2313282A1 CA2313282A1 (en) | 2001-01-12 |
| CA2313282C true CA2313282C (en) | 2006-06-06 |
Family
ID=26840992
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CA002313282A Expired - Fee Related CA2313282C (en) | 1999-07-12 | 2000-06-30 | Method and apparatus for determining logging tool displacements |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US6459992B1 (en) |
| CA (1) | CA2313282C (en) |
| GB (1) | GB2352818B (en) |
| NO (1) | NO329427B1 (en) |
Families Citing this family (22)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020077568A1 (en) * | 2000-11-22 | 2002-06-20 | Haddock Thomas F. | Biological vessel volume measurement method and apparatus utilizing micro accelerometer |
| US6769497B2 (en) * | 2001-06-14 | 2004-08-03 | Baker Hughes Incorporated | Use of axial accelerometer for estimation of instantaneous ROP downhole for LWD and wireline applications |
| US6518756B1 (en) * | 2001-06-14 | 2003-02-11 | Halliburton Energy Services, Inc. | Systems and methods for determining motion tool parameters in borehole logging |
| US6651496B2 (en) | 2001-09-04 | 2003-11-25 | Scientific Drilling International | Inertially-stabilized magnetometer measuring apparatus for use in a borehole rotary environment |
| US6859032B2 (en) * | 2001-12-18 | 2005-02-22 | Schlumberger Technology Corporation | Method for determining molecular properties of hydrocarbon mixtures from NMR data |
| US6838875B2 (en) * | 2002-05-10 | 2005-01-04 | Schlumberger Technology Corporation | Processing NMR data in the presence of coherent ringing |
| US6703832B2 (en) * | 2002-08-12 | 2004-03-09 | Schlumberger Technology Corporation | Method for detecting hydrocarbons by comparing NMR response at different depths of investigation |
| US6859033B2 (en) * | 2002-08-28 | 2005-02-22 | Schlumberger Technology Corporation | Method for magnetic resonance fluid characterization |
| US6781371B2 (en) * | 2002-09-06 | 2004-08-24 | Schlumberger Technology Corporation | High vertical resolution antennas for NMR logging |
| US6960913B2 (en) * | 2003-01-14 | 2005-11-01 | Schlumberger Technology Corporation | Multi-measurement NMR analysis based on maximum entropy |
| US7034528B2 (en) * | 2003-02-27 | 2006-04-25 | Schlumberger Technology Corporation | Methods for formation evaluation based on multi-dimensional representation of nuclear magnetic resonance data |
| US7463027B2 (en) * | 2003-05-02 | 2008-12-09 | Halliburton Energy Services, Inc. | Systems and methods for deep-looking NMR logging |
| US20090126486A1 (en) * | 2007-11-20 | 2009-05-21 | Baker Hughes Incorporated | Orientation independent gravity sensor |
| US7819183B2 (en) * | 2008-06-16 | 2010-10-26 | Halliburton Energy Services, Inc. | Work string controller |
| US8016050B2 (en) * | 2008-11-03 | 2011-09-13 | Baker Hughes Incorporated | Methods and apparatuses for estimating drill bit cutting effectiveness |
| US8028764B2 (en) * | 2009-02-24 | 2011-10-04 | Baker Hughes Incorporated | Methods and apparatuses for estimating drill bit condition |
| CH707573A1 (en) * | 2013-02-07 | 2014-08-15 | Thomas Frizlen | Method and system for determining the displacement of an anchor. |
| US20180347339A1 (en) * | 2015-10-06 | 2018-12-06 | Halliburton Energy Services, Inc. | Systems and Methods for Detecting Downhole Tool Location Inside a Borehole |
| US20180003028A1 (en) * | 2016-06-29 | 2018-01-04 | New Mexico Tech Research Foundation | Downhole measurement system |
| CN106321074B (en) * | 2016-10-18 | 2023-10-24 | 海斯比得（武汉）石油科技有限公司 | Optical engine-based inclinometer while drilling and measuring device thereof |
| RU2668654C1 (en) * | 2017-08-24 | 2018-10-02 | Федеральное государственное бюджетное учреждение науки Институт геофизики им. Ю.П. Булашевича Уральского отделения Российской академии наук (ИГФ УрО РАН) | Device for measuring geoacoustic signals in well |
| US10625824B2 (en) | 2018-01-13 | 2020-04-21 | Thomas Frizlen | Method and system for determining displacement of an anchor |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE4029215A1 (en) | 1990-09-14 | 1992-04-23 | Deutsche Forsch Luft Raumfahrt | Accurate measurement of spatial angles, trajectories, contours etc. - uses sequential angle and/or velocity and/or position measuring w.r.t. unknown, but identical starting states |
| US5541587A (en) | 1995-01-19 | 1996-07-30 | Western Atlas International, Inc. | System for determining the true depth of an electrical logging tool within a wellbore |
| US6051973A (en) | 1996-12-30 | 2000-04-18 | Numar Corporation | Method for formation evaluation while drilling |
| EP1057047B1 (en) * | 1998-01-16 | 2008-08-06 | Halliburton Energy Services, Inc. | Method and apparatus for nuclear magnetic resonance measuring while drilling |
-
2000
- 2000-06-21 US US09/598,629 patent/US6459992B1/en not_active Expired - Lifetime
- 2000-06-30 GB GB0015942A patent/GB2352818B/en not_active Expired - Fee Related
- 2000-06-30 CA CA002313282A patent/CA2313282C/en not_active Expired - Fee Related
- 2000-07-11 NO NO20003562A patent/NO329427B1/en not_active IP Right Cessation
Also Published As
| Publication number | Publication date |
|---|---|
| NO329427B1 (en) | 2010-10-18 |
| GB2352818B (en) | 2001-09-26 |
| GB0015942D0 (en) | 2000-08-23 |
| GB2352818A (en) | 2001-02-07 |
| CA2313282A1 (en) | 2001-01-12 |
| US6459992B1 (en) | 2002-10-01 |
| NO20003562L (en) | 2001-01-15 |
| NO20003562D0 (en) | 2000-07-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2313282C (en) | Method and apparatus for determining logging tool displacements | |
| US20190271224A1 (en) | Correction of motion effect in nuclear magnetic resonance (nmr) logging | |
| US7196516B2 (en) | Correction of NMR artifacts due to constant-velocity axial motion and spin-lattice relaxation | |
| US8004279B2 (en) | Real-time NMR distribution while drilling | |
| CA2449787C (en) | System and methods for formation evaluation while drilling | |
| CA2010398C (en) | Method to improve directional survey accuracy | |
| US8022698B2 (en) | Joint compression of multiple echo trains using principal component analysis and independent component analysis | |
| EP3423675B1 (en) | Motion detection and correction of magnetic resonance data | |
| US20130214779A1 (en) | Method and system to characterize a property of an earth formation | |
| EP1332382B1 (en) | Method for processing nmr data without phase-alternating-pair (pap) averaging | |
| US10302801B2 (en) | Temperature correction of magnetic resonance data | |
| US6522138B2 (en) | Resolution enhancement for sequential phase alternated pair nuclear magnetic resonance measurements | |
| EP1384097B1 (en) | Estimation of transversal motion of the nmr tool during logging | |
| US20110050234A1 (en) | Resonance Method for Measuring Water-Oil Ratio, Conductivity, Porosity, Permeability and Electrokinetic Constant in Porous Formations | |
| GB2361773A (en) | Determining logging tool displacements | |
| CA2711494C (en) | Joint compression of multiple echo trains using principal component analysis and independent component analysis | |
| US20210103071A1 (en) | Reducing Data Bandwidth Requirements in Downhole Nuclear Magnetic Resonance Processing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| EEER | Examination request | ||
| MKLA | Lapsed |
Effective date: 20160630 |