Note: Descriptions are shown in the official language in which they were submitted.
CA 02285791 2002-05-13
77483-41
APPARATUS AND METHOD FOR GENERATING A PULSE SEQUENCE
Background of the Invention
The present invention relates generally to an apparatus and method for
measuring
nuclear magnetic resonance properties of an earth formation traversed by a
borehole, and
more particularly, to an apparatus and method for generating a pulse sequence.
It is well recognized that atomic particles of an earth formation having non-
zero
to nuclear spin magnetic moment, for example protons, have a tendency to align
with a
static magnetic field imposed on the formation. Such a magnetic field may be
naturally
generated, as is the case for the earth's magnetic field; BE. An RF pulse
applying a second
magnetic field transverse to BE creates a magnetization component in the
transverse plane
(perpendicular to B~ which precesses about the BE vector with a characteristic
resonance
known as the Larmor frequency, cuL, which depends on the strength of the
static magnetic
field and the gyromagnetic ratio of the particle. Hydrogen nuclei (protons)
precessing
about a magnetic field BE of 0:5 gauss; for example, have a characteristic
frequency of
approximately 2kHz. If a population of hydrogen nuclei were made to precess in
phase,
the combined magnetic fields of the protons can generate a detectable
oscillating voltage
in a receiver coil, conditions known to those skilled in the art as free
induction decay or a
spin echo. Hydrogen nuclei of water and hydrocarbons occurring in rock pores
produce
CA 02285791 1999-10-13
24.787-CIP3
nuclear magnetic resonance (NMR) signals distinct from signals arising from
other
solids.
U.S. Pat. Nos. 4,717,878 issued to Taicher et al. and 5,055,787 issued to
Kleinberg et al., describe NMR tools which employ permanent magnets to
polarize
hydrogen nuclei and generate a static magnetic field, Bo, and RF antennas to
excite and
detect nuclear magnetic resonance to determine porosity, free fluid ratio, and
permeability
of a formation. The atomic nuclei align with the applied field, Bo, with a
time constant of
T,. After a period of polarization, the angle between the nuclear
magnetization and the
applied field can be changed by applying an RF field, B, , perpendicular to
the static
l0 field Bo, at the Larmor frequency f~ = yB~/2~, where y is the gyromagnetic
ratio of the
proton and Bo designates the static magnetic field strength. After termination
of the RF
pulse, the protons precess in the plane perpendicular to Bo. A sequence of
refocusing RF
pulses generates a sequence of spin-echoes which produce a detectable NMR
signal in the
antenna.
U. S. Pat. No. 5,280,243 issued to Melvin Miller describes a nuclear magnetic
resonance tool for formation evaluation while drilling. The tool includes a
probe section
consisting of a permanent magnet disposed in a longitudinally extending
annular recess
outside the drill collar and an antenna disposed on a non-conductive magnetic
sleeve
outside the drill collar. The gradient of the static magnetic field magnitude
is in the radial
2o direction. The antenna produces an RF magnetic field substantially
perpendicular to both
the longitudinal axis of the tool and the static field direction. With the
'243 apparatus, the
2
CA 02285791 1999-10-13
24.787-CIP3
magnet must be long in axial extent compared to its diameter for the magnetic
fields to
approximate its intended 2-D dipole behavior.
U. S. 5,757,186 issued to Taicher et al. describes a measurement-while-
drilling
tool which includes a sensing apparatus for making nuclear magnetic resonance
measurements of the earth formation. The NMR sensing apparatus is mounted in
an
annular recess formed into the exterior surface of the drill collar. In one
embodiment, a
flux closure is inserted into the recess. A magnet is disposed on the outer
radial surface of
the flux closure. The magnet is constructed from a plurality of radial
segments which are
magnetized radially outward from the longitudinal axis of the tool. The flux
closure is
to required to provide suitable directional orientation of the magnetic field.
The tools disclosed in the '243 and '186 patents suffer from common problems:
both tools require using a nonconductive magnet and placing the magnet outside
the drill
collar. For the '243 tool, the outside surface of the drill collar must
contain a recessed
area to accommodate the nonconductive magnet. For the ' 186 tool, the outside
surface of
the drill collar must contain a recessed area to accommodate the flux closure,
nonconductive magnet, and antenna. Because the strength of the drill collar is
a function
of its radii, reducing the external diameter to accommodate the magnet only or
the flux
closure, magnet, and antenna results in an unacceptably weak section of the
drill collar
which may bend or break during the drilling operation
2o U. S. Pat. No. 5,557,201 issued to Kleinberg et al. describes a pulsed
nuclear
magnetism tool for formation evaluation while drilling. The tool includes a
drill bit, drill
string, and a pulsed nuclear magnetic resonance device housed within a drill
collar made
of nonmagnetic alloy. The tool includes a channel, within the drill string and
pulsed
3
CA 02285791 1999-10-13
24.787-CIP3
NMR device, through which drilling mud is pumped into the borehole. The pulsed
NMR
device comprises two tubular magnets, which are mounted with like poles facing
each
other, surrounding the channel, and an antenna coil mounted in an exterior
surface of the
drill string between the magnets. This tool is designed to resonate nuclei at
a
measurement region known to those skilled in the art as the saddle point.
U. S. Pat. No. 5,705,927 issued to Sezginer et al. also describes a pulsed
nuclear
magnetism tool for formation evaluation while drilling. The tool includes
shimming
magnets, located either inside or outside the tool, which suppress the
magnetic resonance
signal of the borehole fluids by raising the magnitude of the static magnetic
field in the
1o borehole so that the Larmor frequency in the borehole is above the
frequency of the
oscillating field produced by an RF antenna located in a recessed area of the
tool. The
shimming magnets also reduce the gradient of the static magnetic field in the
region of
investigation.
4
CA 02285791 1999-10-13
24.787-CIP3
Summary of the Invention
A substantially axisymmetric static magnetic field is generated in a formation
traversed by a borehole. A pulse programmer is provided for adaptively
controlling the
creation and application of a sequence of RF magnetic field pulses to the
formation. In
order to control the RF pulses, a portion of a memory device is partitioned
into a plurality
of tables and control information is stored in the tables. The tables may
include the
following: a buffer table which describes the layout of stacking buffers, an
acquisition
table which defines the acquired signals accumulated in bufFers, a filter
coefficient table
which prescribes the detection filter employed with a signal acquisition, a
spin dynamics
1 o correction table which designates the spin dynamics correction to be used
for each buffer,
and a data processing table which designates the nuclear magnetic resonance
characteristic calculated from the acquired buffers. In addition to the pulse
programmer,
other tasks may share the control information stored in the tables.
The pulse programmer further comprises a pulse sequence template, used for
generating pulse sequences, which includes a sequence of states dependent on
repetition
and timing variables. These variables are calculated from sequence co~guration
parameters using a calculation block. The timing variables include the wait
time, the
echo spacing, and the number of acquired echoes. The configuration parameters
include
the duration of an excitation pulse, the pulse amplitude, and the pulse shape.
The
sequence of states further comprises a plurality of alternative states for
portions of the
pulse sequence. In real time, one of the alternative states may be selected
based on
external conditions of the apparatus.
5
CA 02285791 2002-12-13
77483-41
In accordance with a first aspect, there is
provided a method for determining a nuclear magnetic
resonance property in an investigation region of earth
formations surrounding a borehole, comprising the steps of:
a) generating a substantially axisymmetric static magnetic
field into a formation traversed by the borehole with a
logging tool; b) adaptively controlling the creation and
application of a sequence of RF magnetic field pulses to the
formation, comprising the steps of: i) partitioning a
portion of a memory device into a plurality of tables and
storing control information in the table; ii) determining an
operating condition for the tool; iii) selecting the
sequence of RF magnetic pulses based upon the control
information and the operating condition; and, iv) defining a
template for executing the selected sequence of RF magnetic
field pulses; and, c) detecting nuclear magnetic resonance
signals from the formation.
In accordance with a second aspect, there is
provided an apparatus for determining a nuclear magnetic
resonance property in an investigation region of earth
formations surrounding a borehole, comprising: a) a logging
tool having at least one magnet for generating a
substantially axisymmetric static magnetic field into a
formation traversed by the borehole; b) a pulse programmer
for adaptively controlling the creation and application of a
sequence of RF magnetic field pulses to the formation,
comprising: i) a memory device, means for partitioning a
portion of the memory device into a plurality of tables, and
means for storing control information in the tables; ii)
means for determining an operating condition for the tool;
iii) means for selecting the sequence of RF magnetic pulses
based upon the control information and the operation
condition; iv) means for defining a template for executing
5a
CA 02285791 2002-12-13
77483-41
the sequence of RF magnetic field pulses; and c) means for
detecting nuclear magnetic resonance signals from the
formation.
5b
CA 02285791 1999-10-13
24.787-CIP3
Brief Description of the Drawings
The advantages of the present invention will become apparent from the
following
description of the accompanying drawings. It is to be understood that the
drawings are to
be used for the purpose of illustration only, and not as a definition of the
invention.
In the drawings:
Fig. 1 illustrates a logging-while-drilling apparatus;
Fig. 2 depicts the low gradient sonde;
Figs. 2a-2d illustrate the contour lines Bo corresponding to four low
gradient magnet configurations;
1o Figs. 3a-3d represent the contour lines of the gradient wBo~ corresponding
to four low gradient magnet configurations;
Fig. 4 depicts the high gradient sonde;
Fig. 4a represents the contour lines IBoI corresponding to the high gradient
magnet configuration;
Fig. 4b represents the contour lines of the gradient wBo ~ corresponding to
the high gradient magnet configuration;
Fig. 5 depicts the simple data acquisition mode;
Fig. 6 depicts the interleaved data acquisition mode;
Fig. 7 depicts the burst data acquisition mode; and,
2o Fig. 8 represents a block diagram of the pulse programmer.
6
CA 02285791 1999-10-13
24.787-CIP3
Detailed Description of the Preferred Embodiment
Referring to Fig. 1, a nuclear magnetic resonance (NMR) logging-while-drilling
tool 10 is illustrated. The tool 10 comprises a drill bit 12, drill string 14,
a plurality of RF
antennas 36, 38, and at least one gradient coil 56. The tool 10 further
comprises electronic
circuitry 20 housed within the drill collar 22. The electronic circuitry 20
comprises RF
resonance circuitry for the antennas 36, 38, a microprocessor, a digital
signal processor,
and a low voltage bus. The tool 10 further comprises a plurality of tubular
magnets 30,
32, and 34 which are polarized in a direction parallel to the longitudinal
axis of the tool
but opposite to each other, i.e., with like magnetic poles facing each other.
The
1o magnets 30, 32, and 34 comprise either a conductive or nonconductive
material. The
configuration of magnets 30, 32, and 34 and antennas 36, 38 provides for at
least two
NMR regions of investigation 60, 62 with a substantially axisymmetric static
and RF
magnetic field.
A means for drilling a borehole 24 in the formation comprises drill bit 12 and
drill
collar 22. The drill collar 22 may include a stabilizing means (not shown) for
stabilizing
radial motion of the tool 10 in the borehole during drilling, however, the
stabilizing
means is not required; therefore, the tool 10 may operate unstabilized or
stabilized. Mud
flow sleeve 28 defines a channel 90 for carrying the drilling fluid through
the drill string
14. A drive mechanism 26 rotates the drill bit 12 and drill string 14. This
drive
mechanism is adequately described in U.S. Pat. No. 4,949,045 issued to Clark
et al.
However, a downhole mud motor may be placed in the drill string as the drive
mechanism 26.
7
CA 02285791 1999-10-13
24.787-CIP3
It is within contemplation of the subject invention to combine N+1 magnets to
obtain at least N regions of investigation in the formation. The combinations
contemplated by this invention include, but are not limited to, a low gradient-
low
gradient, high gradient-high gradient, high gradient-low gradient, low
gradient-high
gradient, or a combination of high gradient, low gradient, and saddle point
regions. The
combination of high and low gradient static field regions in the formation
offers several
advantages. For example, the high gradient region may have a higher signal-to-
noise -
ratio but may experience signal loss when the tool 10 undergoes lateral motion
in the
borehole. On the other hand, the low gradient region has lower susceptibility
to signal
to loss problems when the tool 10 is in motion. Also, with moderate tool
motion, longer
echo trains can be acquired in the low gradient region than in the high
gradient region
thereby providing more information about permeability, bound and free fluid,
and
hydrocarbon types. Moreover, the combination of data acquired with both
gradient
regions may provide quantitative information about the amount of lateral
motion the tool
10 experiences and can be used to motion correct the NMR data, or, at least,
quality
control the data. Measurements of devices, such as strain gauges,
accelerometers, or
magnetometers, or any combinations of these devices, may be integrated with
NMR
information to quality control the data or make corrections to the spin-echo
train. With
the combination of high and low gradient static magnetic fields, the high
gradient region
exhibits more diffusion effect and therefore is of greater interest for
hydrocarbon typing
techniques than the low gradient region. Finally, the low gradient region has
a static
magnetic field having a low amplitude and therefore, this region with its
lower Larmor
frequency is less affected by formation and borehole fluid conductivity.
8
CA 02285791 1999-10-13
24.787-CIP3
Low Gradient Sonde
Referring to Fig. 2, in one section of the tool, hereinafter referred to as
the low
gradient sonde, a central magnet 30 is axially separated from a lower magnet
32. These
magnets 30, 32 generate a substantially axisymmetric static magnetic field
that is radial in
its polarization and, over a reasonably long cylindrical shell, the static
magnetic field has
a fairly constant magnitude. It is within contemplation of the subject
invention to excite a
plurality of cylindrical shells of spins in the formation where each shell is
resonant at a
different RF frequency, and to sequentially interrogate each shell with
sequences of RF
pulses.
l0 The area between magnets 30, 32 is suitable for housing elements such as
electronic components, an RF antenna, and other similar items. For example, a
plurality
of electronic pockets 70 may form an integral part of the mud sleeve 28. These
pockets
70 may house the RF circuitry (e.g., Q-switch, duplexer, and pre-amplifier),
preferably in
close proximity of the RF antenna. In a preferred embodiment of the invention,
the
pockets 70 form an integral part of magnetically permeable member 16. In that
case, to
maintain the axial symmetry of the magnetic field, a highly magnetically
permeable cover
72 is located over each pocket 70.
The magnetically permeable member 16 is positioned inside the drill collar 22
between the magnets 30, 32. Member 16 may consist of a single piece or a
plurality of
2o sections combined between the magnets 30, 32. Member 16 is constructed of a
suitable
magnetically permeable material, such as fernte, permeable steel or another
alloy of iron
and nickel, corrosion resistant permeable steel, or permeable steel having a
structural role
in the member design, such as 15-5 Ph stainless steel. The magnetically
permeable
9
CA 02285791 1999-10-13
24.787-CIP3
member 16 focuses the magnetic field and may also either carry drilling fluid
through the
drill string or provide structural support to the drill collar. Further,
member 16 improves
the shape of the static magnetic field generated by magnets 30, 32 and
minimizes
variations of the static magnetic field due to vertical and lateral tool
motion during the
period of acquiring the NMR signal. The segment of sleeve 28 between magnets
30, 32
may comprise magnetically permeable member 16. In that case, the segments of
sleeve
28 under magnets 30, 32 shall consist of a non-magnetic member. Alternatively,
a
magnetically permeable chassis surrounding the segment of sleeve 28 between
magnets
30, 32 defines member 16. In this case, the segment may consist of a magnetic
or non-
to magnetic material. It is within contemplation of this invention to
integrate the chassis
and segment to form member 16.
The magnets 30, 32 are polarized in a direction parallel to the longitudinal
axis of
the tool 10 with like magnetic poles facing each other. For each magnet 30,
32, the
magnetic lines of induction travel outward from an end of the magnet 30, 32
into the
formation, along the axis of the tool 10, and travel inward to the other end
of the magnet
30, 32. In the region between central magnet 30 and lower magnet 32, the
magnetic lines
of induction travel from the center outward into the formation, creating a
static field in a
direction substantially perpendicular to the axis of the tool 10. The magnetic
lines of
induction then travel inward symmetrically above the central magnet 30 and
below the
lower magnet 32 and converge in the longitudinal direction inside sleeve 28.
Because of
the separation, the magnitude of the static magnetic field in the central
region between the
central 30 and lower 32 magnet is spatially homogeneous in comparison to a
saddle-point
field.
CA 02285791 1999-10-13
24.787-CIP3
The amount of separation between the magnets 30, 32 is determined based on
several factors: ( 1 ) selecting the requisite magnetic field strength and
homogeneity
characteristics; (2) generating a field having small radial variations in the
region of
interest so that the echoes received during a pulse sequence (i.e., CPMG, CPI,
or other
sequences) are less sensitive to lateral tool motion; (3) depth of
investigation; and (4)
minimizing interference between the resonance circuitry and the low voltage
telemetry
bus in order to improve isolation of the receiving antenna which detects NMR
signals
from the formation. As the separation between the magnets 30, 32 decreases,
the
magnetic field becomes stronger and less homogeneous. Conversely, as the
separation
1o between the magnets 30, 32 increases, the magnetic field becomes weaker and
more
homogenous.
Figs. 2a-2d illustrate the contour lines of IBolcorresponding to four
laboratory
modeled configurations of central 30 and lower 32 magnets. These modeled
results were
computed using a tool having a preselected diameter (a constant diameter was
used for
modeling all configurations). The configuration corresponding to Fig. 2a
comprises a
non-magnetically permeable member separating a central 30 and lower 32 magnet
by 25
inches. The configuration corresponding to Fig. 2b comprises a non-
magnetically
permeable member separating a central 30 and lower 32 magnet by 18 inches. The
configuration corresponding to Fig. 2c comprises a non-magnetically permeable
member
2o separating a central 30 and lower 32 magnet by eight inches. The low
gradient sonde,
corresponding to Fig. 2d, comprises a magnetically permeable member 16
separating a
central 34 and lower 32 magnet by 25 inches. The aforementioned dimensions
were
CA 02285791 1999-10-13
24.787-CIP3
modeled to merely illustrate the affect of distance and/or a magnetically or
non-
magnetically permeable member on IBoI. Figs. 3a-3d represent the contour lines
of the
gradient wBo ~ corresponding respectively to configurations illustrated in
Figs. 2a-2d.
In the low gradient sonde, the magnetically permeable member 16 shunts a
significant portion of the magnetic flux into the center of the tool 10. To
illustrate, the
magnitude of the Bo field shown in Fig. 2d at a distance of approximately
seven inches
radially from the longitudinal axis of toil 10 is twice as large as the Bo
field shown in
Fig. 2a which was generated by the same magnet configuration separated by a
non-
magnetically permeable member. Furthermore, the low gradient sonde produces a
longer
1 o and more uniform extent of the static magnetic field in the axial
direction. The NMR
signal measured in this section of the tool is substantially less sensitive to
the vertical
motion of the tool. Referring to Fig. 3d, with the low gradient sonde, a
relatively small,
approximately 3 Gauss/cm, gradient is measured at a distance of approximately
seven
inches radially from the longitudinal axis of tool. This low gradient results
in a measured
NMR signal which is substantially less sensitive to the lateral motion of the
tool 10.
When motion is moderate, longer echo trains may be acquired in this region
thereby
providing more information about permeability, bound and free fluid, and
hydrocarbon
types. In the case of the low gradient sonde, as with other gradient designs,
the proton
rich borehole region surrounding the tool 10 will resonate only at frequencies
higher than
2o those being applied to the volume of investigation, i.e., there is no
proton borehole signal.
Other NMR sensitive nuclei found in drilling mud, such as sodium-23, resonate
at
significantly higher static magnetic field strengths than hydrogen when
excited at the
12
CA 02285791 1999-10-13
24.787-CIP3
same RF frequency. For the low gradient sonde, these higher field strengths
are not
produced in the borehole region surrounding the tool or near the antenna where
such
unwanted signals could be detected.
High Gradient Sonde
Refernng to Fig. 4, in another section of the tool, hereinafter referred to as
the
high gradient sonde, a central magnet 30 is axially separated from an upper
magnet 34.
The magnets 30, 34 are polarized in a direction parallel to the longitudinal
axis of the tool
with like magnetic poles facing each other. These magnets 30, 34 generate a
substantially axisymmetric static magnetic field that is radial in its
polarization and, over
1 o a reasonably long cylindrical shell, the static magnetic field has a
fairly constant
magnitude. It is within contemplation of the subject invention to excite a
plurality of
cylindrical shells of spins in the formation where each shell is resonant at a
different RF
frequency.
As illustrated in Fig. 2c, if the magnet separation between 30 and 34 is
approximately eight inches, the contour lines of the static magnetic field
strength are
substantially straight and the strength of IBo I is greater than the static
magnetic field
strength of the low gradient region. However, the gradient wBo ~ becomes
larger, as
illustrated in Fig. 3c, at a distance of approximately seven inches radially
from the
longitudinal axis of the tool. The contour lines of wBo ~ are curved denoting
variation of
the gradient in the axial direction.
13
CA 02285791 1999-10-13
24.787-C1,P3
The high gradient sonde is improved by inserting a magnetically permeable
member 16 between magnets 30, 34. Fig. 4a represents contour lines of IBoI
corresponding to a configuration where magnetically permeable member 16
separates the
upper 34 and central 30 magnets by eight inches. The contour lines of Fig. 4a
show a
slightly stronger field indicating a better signal-to-noise ratio and less
curvature in the
axial direction than the contour lines of Fig. 2c. Also, as illustrated in
Fig. 4b, the
magnetically permeable member 16 produces a more constant gradient wBo ~ in
the axial
direction that may simplify interpretation of the NMR measurements influenced
by
diffusion.
1o In the case of the high gradient sonde, as with other gradient designs, the
proton
rich borehole region surrounding the tool 10 will resonate only at frequencies
higher than
those being applied to the volume of investigation, i.e., there is no proton
borehole signal.
The high gradient sonde is sensitive to a small part of the sodium from the
borehole fluid.
For a 30% NaCI concentration borehole fluid, possibly the worst case, the
error of the
estimated porosity due to the sodium signal is approximately 0.08 pu. In the
low gradient
sonde, the sodium signal is substantially smaller than in the high gradient
sonde.
Consequently, the sodium signal is negligible for both NMR sondes.
Antennas and Gradient Coils
Referring to Figs. 2 and 4, an RF magnetic field is created in the regions of
2o investigation by antennas 36, 38 which are provided in recessed areas 50,
52. The RF
field may be produced by one or more RF antenna segments that transmit and/or
receive
14
CA 02285791 1999-10-13
24.787-CIP3
from different circumferential sectors of the logging device. See U.S. Pat.
App. Nos.
08/880,343 and 09/094,201 (Atty. Docket. Nos. 24.784 and 24.784-CIP) assigned
to
Schlumberger Technology Corporation. Preferably, each antenna 36, 38 comprises
a coil
18 wound circumferentially around the recessed area 50, 52. The RF field
created by
such a coil arrangement is substantially axisymmetric. It is within
contemplation of the
subject invention to utilize the antenna 36, 38 for detecting NMR signals.
However, a
separate antenna or receiver may be used to detect the signals. A non-
conductive material
54 is provided in the recessed area 50, 52 beneath the antenna 36, 38. The
material 54 is
preferably a ferrite to increase the efficiency of the antenna 36, 38.
Alternatively, the
t o material 54 may comprise plastic, rubber, or a reinforced epoxy composite
material. The
antennas 36, 38 are resonated by RF circuitry to create an RF magnetic field
in the
regions of investigation.
The recessed area 52 forms a shallow groove in the drill collar without
reducing
the inner diameter of the drill collar, which is ordinarily done to increase
strength in a
region of drill collar where the outer diameter has been recessed to provide
an antenna.
The recessed area 50 has a greater depth than recessed area 52. Due to
mechanical
constraints, it is only possible to have one deeply recessed area where the
drill collar
inner diameter is substantially reduced. It is within contemplation of the
subject invention
for the recessed areas 50, 52 to have substantially the same depth or for
recessed area 52
2o to have a greater depth than area 50.
The cylindrical shells of spins in the region of investigation can be
segmented
axially or, preferably, azimuthally by using at least one directionally
sensitive gradient
coil 56 arranged in the recessed area 50 and/or 52. In a preferred embodiment
of the
CA 02285791 1999-10-13
24.787-CIP3
invention, three gradient coils are positioned circumferentially around the
recessed area
and separated by an angular distance segment of 120'. Other quantities of
gradient coils
may be defined, either lesser or greater in number than three, and such coils
may be
separated by angular distances other than 120' and/or unequal angular
segments. Each
coil 56 is constructed with loops of wire, which conform to the curvature of
the outer
surface of the material 54. The magnetic field produced by each gradient coil
56 in a
region of the formation facing the coil is substantially parallel to the
static magnetic field
produced by the magnets.
As is known to those skilled in the art, in the basic NMR measurement, a pulse
1o sequence is applied to the formation under investigation. In U.S. Pat. No.
5,596,274
issued to Abdurrahman Sezginer and U.S. Pat. No. 5,023,551 issued to Kleinberg
et al., a
pulse sequence, such as the Carr-Purcell-Meiboom-Gill (CPMG) sequence, first
applies
an excitation pulse, a 90' pulse, to the formation that rotates the spins into
the transverse
plane. After the spins are rotated by 90' and start to dephase, the Garner of
the refocusing
pulses, the 180' pulses, is phase shifted relative to the carrier of the 90'
pulse sequence
according to the following relationship:
r9o*1 -ro -~tl$a~ -rl -eChom~ -r,, , where the bracketed expression is
repeated for
n=1,2,...N, where N is the number of echoes collected in a single CPMG
sequence and
the echo spacing 1S recho =2r~P =r,$o~ +r, +r2 . 90ts denotes an RF pulse that
causes the
2o spins to rotate by a 90° angle about the +x-axis, as commonly
defined in the rotating
frame of magnetic resonance measurements (phase alternated). The time between
application of the 90° pulse and the 180° pulse, to, is less
than t~P, half the echo spacing.
16
CA 02285791 1999-10-13
za.~s7-c~3
The CPMG sequence enables acquisition of a symmetric measurement (i.e., a
measurement without using the gradient coils). The exact timing parameters,
to, r,, and t2,
depend on various factors (e.g., the shape of the applied pulses).
In the subject invention, a current pulse applied to gradient coil 56
generates an
additional magnetic field, substantially parallel to the static magnetic
field. The current
pulse is applied between the first 90° and the 180° phase
reversing pulse. This additional
field causes an additional phase shift for the spins. Since the 180°
phase reversing pulse
does not compensate for the additional phase shift, the spins subjected to the
additional
field do not form a spin-echo. However, for spins not subjected to the
additional field, a
to spin-echo occurs at time 2t~ with spin-echoes of successively lower
amplitude occurring
at time t~p after each phase reversing pulse. The pulse sequence is
t9o*~ -ro -s-ro -[r~$oy -r, -echom~ -rz, , where ro is the time between the
90° pulse and
n
the gradient pulse of duration 8 , ro is the time between the gradient pulse
and the 180°
reversing pulse, and ro + 8+ ro = ro . Due to the succeeding I Boy pulses and
the
inhomogeneous fields, the x-component of the NMR signal will decay within a
few
echoes. Therefore, we focus only on the y-component of the signal. Thus,
neglecting
relaxation, the first NMR echo signal can be represented as:
Signal = ~m~~ERl(M° +iM°)(r)exp(-iyG(r)8)dc(r)J,
where i is the imaginary complex unit; y is the gyromagnetic ratio; M°
and My are
2o respectively x and y components of the magnetization at location r at the
time of the first
echo in the absence of the gradient pulse; G(r) is the component of the
gradient field
17
CA 02285791 1999-10-13
24.787-CIP3
parallel to Bo at the same location; 8 is the duration of the gradient pulse;
and dc(r)
denotes the differential sensitivity of the NMR sonde.
The gradient coils 56 offer a number of advantages for obtaining azimuthal
measurements. First, because the axisymmetric antenna detects the spin-echoes,
long
echo trains can be recorded while the tool rotates in the borehole. Second,
the coil 56
simplifies the design of an NMR-LWD tool because the coil 56 does not have the
tuning
requirements of an RF antenna 36, 38. Third, the same antenna 36, 38 can be
used to
make symmetric and axisymmetric measurements. Fourth, the coils 56 can be used
to
obtain NMR measurements with excellent spatial resolution, particularly
vertical
to resolution.
Different modes for obtaining azimuthal NMR measurements are contemplated by
the present invention. For example, a "simple spoiling" mode uses at least one
coil 56 to
spoil the spins in a selected quadrant where a quadrant is defined as an
angular distance
segment about the periphery of the tool 10, however, more coils 56 may be used
to spoil a
plurality of quadrants. In either case, two measurements are obtained: a
symmetric phase
alternated pulse sequence (PADS) with a fixed wait time followed by a gradient
PADS
having a variable wait time, with the selected quadrant spoiled by firing the
coil 56 in the
quadrant. In a preferred embodiment of the invention, the aforementioned
gradient pulse
sequence is used. Subtracting the gradient measurement from the symmetric
2o measurement creates the azimuthal measurement. In this mode, one symmetric
measurement is obtained for every two PADS and one azimuthal scan is obtained
for
every eight PAPS. The measurement noise for the azimuthal measurement is
higher than
~s
i
CA 02285791 2002-05-29
77483-41
the noise in the symmetric or gradient measurement because the two
measurements are
combined.
It is possible to reduce the noise contribution by combining different single
quadrant spoiling measurements. For example, four gradient PADS measurements
may
be obtained by spoiling each quadrant. The measurements are combined to create
a
synthetic azimuthal and symmetric measurement. By combining measurements made
without the gradient coils 56 being fired with measurements made with one or
more
gradient coils 56 being fired, axially or azimuthally resolved "images" of the
formation
can be generated. The acquired data, particularly in the form of azimuthal
images of
to porosity and bound fluid, are very desirable for improved petrophysical
interpretation in
highly deviated and horizontal wellbores and for decision-making while
drilling for
geologically based wellbore placement.
Ontimizin~ the Pulse Length and Oneratine Freauencv
For a chosen operating RF frequency, there is an optimum duration for the
90°
l5 pulse, t~, as well as for the 180° pulses, t,~, which ensures a
desired signal-to-noise ratio.
The search for an optimal pulse length may be performed during the master
calibration of
the tool, so that all pulse lengths will be correctly initialized, or when the
static magnetic
field changes in an unpredictable manner, such as a change due to accumulation
of
magnetic debris during the drilling process. See U.S.
Patent 6,237,404 assigned to Schlumberger Technology
Corporation. This technique may also be used to choose the
appropriate frequency to meet other criteria, such as
keeping the depth of investigation constant.
19
CA 02285791 1999-10-13
24.787-CIP3
The optimal pulse length may be determined by measuring the NMR response of
a sample using at least two different pulse durations and using a predefined
mode
independent of the NMR properties of the formation. Alternatively, the optimal
pulse
length may be determined using at least two different pulse durations and
additionally
using a mode computed from the NMR properties of the formation. In the first
case,
stacking the data improves the signal-to-noise ratio, however, the stacking
procedure may
require a long period of time to acquire data from the formation. Preferably,
the measured
data are accumulated during a stationary time window when the tool 10 pauses
.from the
drilling operation, such as during the time when a new section of drill pipe
is added to the
l0 drill string. In the second case, if the TZ distribution of the formation
is known, a best
acquisition mode may be constructed which provides the largest signal-to-noise
ratio for a
unit of acquisition time and provides an optimal linear combination of the
acquired
echoes. Laboratory simulations show that optimum timing for the best
acquisition mode
is achieved when the duration of the echo train is approximately equal to Tz,
m~, the
dominate TZ of the formation, and when the wait time, t,~" is approximately
equal to 2.5 x
TZ~,~ (assuming a constant T,/TZ ratio of 1.5). The best acquisition mode
determines the
optimal pulse length to within a few percent over several seconds. A similar
technique
may be used to optimize the NMR signal with respect to the frequency (e.g.,
saddle point
design). The Tz distribution effectively aids the efficient tuning of pulse
lengths for the
tool 10.
CA 02285791 1999-10-13
24.787-CIP3
Data AccLuisition Modes
As described above, tool 10 has a plurality of antennas 36, 38. In a preferred
embodiment of the invention, these antennas 36, 38 do not transmit or acquire
data
simultaneously. Preferably, after one antenna 36 acquires data, the other
antenna 38
experiences a minimum wait time while the power supply recharges in order to
transmit
the next pulse sequence. It is within contemplation of the subject invention
to transmit or
acquire data simultaneously. Further, this invention contemplates data
acquisition without
a requisite wait time.
Based on these design preferences, a plurality of data acquisition modes may
be
1o used. By way of example, three representative timings for NMR data
acquisition are
described below: a fast timing suitable for water-wet sandstone zones, a slow
timing
appropriate for carbonate zones, and a very slow timing designed for
hydrocarbon bearing
zones (or invasion of oil based mud). The timings are set forth in Table I.
Wait Time (sec.)Echo spacing Number of Echoes
(msec)
Fast 2.3 0.5 400
Slow 4.6 0.5 800
very slow9.2 1.0 800
Table I
Several different modes may be used with each data acquisition timing,
including,
but not limited to, the following: simple, interleaved, and bursts. The
simplest way to
acquire T, information with the tool 10 is to perform CPMG measurements with
both
antennas 36, 38 using the same timing. Fig. 5 illustrates the simple data
acquisition mode
21
t
CA 02285791 2002-05-29
77483-41
used with the fast decaying, slowly decaying, and very slowly decaying timing
from
Table I. Each antenna 36, 38 alternately acquires a long pulse sequence which
provides
an effective porosity measurement from each antenna 36, 38.
With the interleaved mode, the high gradient antenna measures at least two
cylindrical shells at tvo .different frequencies while the low gradient
antenna obtains a
measurement using a single frequency. Fig. 6 illustrates an interleaved
measurement for
fast decaying samples, slowly decaying components, and very slowly decaying
components using the timing from Table I.
The burst mode enhances the signal-to-noise ratio,
especially for the fast decaying components. In addition,
to
the burst mode provides a useful T1 based bound fluid
measurement. See U.S. Patent 6,232,77881 assigned to
Schlumberger Technology Corporation which describes a method
for polarizing the bound fluid of a formation. Fig. 7
illustrates burst measurements for fast decaying samples,
slowly decaying components, and very slowly decaying
components using slightly modified times from Table I.
In addition to the simple, interleaved, and burst modes, with the subject
invention,
it is possible to optimize formation evaluation measurements by detecting
downhole
conditions which create a pause during the drilling operation, determining the
drilling
mode, and using the mode to control data acquisition. Standard rotary drilling
operations
contain many natural pauses where the tool remains stationary: connection time
as a new
section of drill pipe is added to the drill string, circulation time when mud
is circulated
and the drill pipe is possibly rotated, and fishing or jarring time while the
drill string is
22
i
CA 02285791 2002-05-29
a. -
77483-41
stuck and has to be freed before drilling can resume. These natural pauses,
which occur
without interrupting normal drilling operations, or deliberately initiated
pauses, are
utilized to make NMR measurements. The drilling modes include, but are not
limited to,
drilling, sliding, tripping, circulating, fishing, a short trip (up or down),
and drill pipe
connections. Determining the drilling mode enhances the ability to obtain NMR
measurements that take a long time or that benefit from a quiet environment,
e.g., T~, Ti,
antenna tuning, and hydrocarbon typing. See U.S. Patent
6,237,404 assigned to Schlumberger Technology Corporation. .
It is also possible to adjust acquisition modes based on
changes in the environment (e. g., washout, salinity, etc.)
and/or changes in the formation NMR properties (e.g., long T1
versus short T1) .
The spin-echo amplitudes are obtained by hardware integration of the receiver
voltages over a time window. The tool 10 uses phase sensitive detection to
measure the
in-phase and quadrature components of the spin-echo signal-plus-noise
amplitudes. The
techniques disclosed in U.S. 5,381,092, issued to Robert Freedman may be used
to
compute window sums downhole and transmit the window sums to the surface for
Ti
inversion processing and presentation. Also, the techniques disclosed in U.S.
5,363,041
issued to Abdurrahman Sezginer may be implemented to utilize a linear operator
to map a
relaxation-time distribution to spin-echoes, produce a singular value
decomposition
(SVD) of the linear operator, determine vectors of the SVD, and compress the
spin-echo
data using the vectors. Preferably, the TZ spectrum is computed downhole and
transmitted
to the surface. This offers the advantage of eliminating a telemetry
bottleneck created by
transmitting the data required to compute the T, spectrum to the surface. A
digital signal
23
CA 02285791 1999-10-13
24.787-CIP3
processor may be used to invert the Tz data. The amplitudes, A~, of the spin-
echoes are
characterized by the following relationship:
M
f1 j = ~ xjiai + TIj ,
i=1
where r~j is the noise in the measurement A j , a; is the amplitude of the Tz
distribution
taken at T,.;, xj; = expC- ~~t~ ~ 1- expC- t-'" ~~ represents the elements of
matrix x , where
T2; cT2;
tW is the wait time and c is a constant (the T,/TZ ratio), 0t is the echo
spacing, and
j=1,2,...N, where N is the number of echoes collected in a single pulse
sequence. In
matrix notation, the equation becomes A = xa' + i~ . Since the noise, r~ , is
unknown, a can
be approximated by finding a least squares solution, i.e., a minimum of the
functional
~ = I~t - xa IZ . The solution of this equation is strongly affected by noise
present in the
data and the solution may have negative components even though the T, spectrum
does
not have negative components. To overcome this problem, a regularization term,
~.Ilallz , is
added to the functional and the functional .1,~ (a) _ ~~;~ - xa~~2 + ~.Ilall2
is minimized using a
suitable iterative minimization algorithm (e.g., Conjugated Gradient
Projection Method)
under the constraint that a; >_ ofor i=1...M. See Ron S. Dembo and Ulrich
Tulowitzski,
On the Minimization of Quadratic Functions Subject to Box Constraints, Yale
Department of Computer Science (September 1984) (describes the Conjugated
Gradient
Projection Method). The necessary time for performing the T, inversion using a
digital
signal processor is very reasonable. For example, assuming 1800 echoes and 30
samples
24
CA 02285791 1999-10-13
24.787-CIP3
in the TZ domain, the inversion on a digital signal processor requires less
than two
seconds.
Pulse Programmer
For the basic NMR measurement with tool 10, electronic circuitry applies a
pulse
sequence to the formation under investigation. Tool 10 includes a pulse
programmer 80,
which adaptively selects and controls the pulse sequences applied to the
formation. The
pulse programmer 80 establishes the pulse sequence using information found in
the
Measurement Control Block 82 (See Fig. 8) and the operating conditions of the
tool 10.
to Preferably, the Measurement Control Block 82 is stored in a downhole memory
device.
The structure of block 82 is fixed to allow pulse programmer 80 easily to
adapt and
change the timing of the pulse sequences autonomously downhole. It is
advantageous to
partition a portion of block 82 into a plurality of tables 84, 86, and 88.
Instead of
controlling all tool operations that depend on the pulse sequence from the
pulse
programmer 80, the tables 84, 86, 88 are used to control these operations.
This allows the
pulse programmer 80 to vary the pulse sequences without introducing
contradictions in
the tool configuration. The plurality of tables 84, 86, and 88 may include,
but are not
limited to, the following: a buffer table which describes the layout of
stacking buffers, an
acquisition table which defines the acquired signals accumulated in buffers, a
filter
2o coe~cient table which prescribes the detection filter employed with a
signal acquisition,
a spin dynamics correction table which designates the spin dynamics correction
to be
used for each buffer, and a data processing table which designates the nuclear
magnetic
resonance characteristic calculated from the acquired buffers.
CA 02285791 1999-10-13
24.787-CIP3
The pulse programmer 80 includes a pulse sequence template 94, useful for
generating pulse sequences, which comprises a sequence of states dependent on
repetition
and timing variables. These variables are calculated from sequence
configuration
parameters using the calculation block 92. The calculation, block 92 may be
implemented
as an executable or interpretive structure. Based on the physical quantity
that will be
measured, e.g. T2, timing variables may be defined such as the wait time, t""
the echo
spacing, t~ho, and the number of acquired echoes. The configuration parameters
include,
but are not limited to, t~, pulse amplitude, and pulse shape. These parameters
may be
calculated periodically during calibration of the tool 10 or during operation
of the tool 10
to since those parameters may vary as the operating conditions of the tool 10
vary. For
example, the pulse amplitude and shape depends on the antenna quality factor
and,
therefore, on the conductivity of the formation surrounding the tool 10.
Normally, after pulse programmer 80 initiates a pulse sequence, the sequence
mns
deterministically until it is finished. To implement certain azimuthal
measurement modes
with tool 10, the pulse programmer 80 has the ability to vary the pulse
sequence during
execution of the sequence. Programmer 80 may stop execution of the pulse
sequence and
enter a HALT state until an external signal ends the state at time t~ or until
a maximum
time period, t",~, has expired. As previously discussed in the Data
Acquisition Modes
section of this specification, since at least one of the different modes
(interleaved) which
may be used with the data acquisition timing contemplates interleaving several
measurements, the programmer 80 compensates for the time that passed during
the
HALT state. Preferably, compensation is accomplished by grouping HALT events.
For
example, a grouping may comprise a pair of HALT events where one HALT event
26
CA 02285791 1999-10-13
24.787-CIP3
operates as previously described and the other HALT event is a normal event of
duration
t",~ t~. Grouping events allows the programmer 80 to combine sequences having
variable
and deterministic timing.
In addition, the sequence of states, as defined in the pulse sequence template
94,
may comprise several alternatives for parts of the sequence. In real time, one
of the
alternatives (branching) is chosen dependent on external conditions of the
tool (e.g., the
azimuth of the tool).
The foregoing description of the preferred and alternate embodiments of the
present invention has been presented for purposes of illustration and
description. It is not
1 o intended to be exhaustive nor to limit the invention to the precise form
disclosed.
Obviously, many modifications and variations will be apparent to those skilled
in the art.
The embodiments were chosen and described in order to best explain the
principles of the
invention and its practical application thereby enabling others skilled in the
art to
understand the invention for various embodiments and with various
modifications as are
suited to the particular use contemplated. It is intended that the scope of
the invention be
defined by the accompanying claims and their equivalents.
27
