|Publication number||US5796252 A|
|Application number||US 08/783,778|
|Publication date||Aug 18, 1998|
|Filing date||Jan 15, 1997|
|Priority date||Mar 23, 1995|
|Also published as||CA2172424A1, CA2172424C|
|Publication number||08783778, 783778, US 5796252 A, US 5796252A, US-A-5796252, US5796252 A, US5796252A|
|Inventors||Robert L. Kleinberg, Lawrence L. Latour, Abdurrahman Sezginer|
|Original Assignee||Schlumberger Technology Corporation|
|Export Citation||BiBTeX, EndNote, RefMan|
|Patent Citations (19), Non-Patent Citations (32), Referenced by (133), Classifications (4), Legal Events (3)|
|External Links: USPTO, USPTO Assignment, Espacenet|
90-tg -pfg-tr -180-tg -pfg-tr -echo1 - ti -180-ti -echoi+1 !i,
90-tg -pfg-tr -180-tg -pfg-tr -echo1 - ti -180-ti -echoi !i,
90-t1 -90-ts -90-t1 -echo
90-t1 -90-ts -90-t1 -echo
90-t1 -90-ts -90-t1 -echo1 - ti -180-ti -echoi+1 !i
90-tg -pfg-tr -180-tg -pfg-tr -echo1 - ti -180-ti -echoi+1 !i,
90-te -180-te -echo1 - ti -180-ti -echoi+1 !i,
90-tg -pfg-tr -90-ts -90-tg -pfg-tr -echo
90-te -180-te -echo1,
90-tg -pfg-tr -180-tg -pfg-tr -echo1,
This application is a file wrapper continuation of parent application Ser. No. 08/394,994 filed Mar. 23, 1995 now abandoned.
This invention relates to apparatus and techniques for determining characteristics of earth formations surrounding a borehole and, more particularly, to an apparatus and method for nuclear magnetic resonance (NMR) borehole logging to determine residual oil saturation and other formation characteristics.
A major goal in the evaluation of hydrocarbon bearing earth formations is the accurate determination of the volumes of oil and water in the pore space of sedimentary rocks. Measurements made with signals from logging instruments have been used to obtain estimates of these volumes. However, these methods are approximate and are useful as rough guides to production strategies. The most credible measurement of producible fluid volumes is to actually produce fluids from the formation, such as by using a drill stem test or by using a logging device that extracts fluids from the formations.
The above methods are useful for estimating oil and water saturations in formations from which oil has not yet been produced. When a new reservoir is discovered, it is first produced by natural hydraulic pressure and then by flooding the formation with water forced down a neighboring well. These are the primary and secondary recovery techniques. They are relatively inexpensive to implement, but leave a substantial amount of residual oil in place. The residual oil is generally about thirty percent of the total oil originally contained in the reservoir.
Methods used to extract the vast amount of residual oil left after secondary recovery are termed tertiary or enhanced oil recovery (EOR). A variety of EOR methods are employed, including steam flooding, miscible fluid sweeps, and other techniques. All such techniques are relatively expensive. The economics of EOR depend critically on the quantity of residual oil in the formation; that is, the residual oil saturation (ROS). Existing methods of determining ROS are often not precise enough to determine the economic viability of EOR projects, and one of the objects of the present invention is to improve the accuracy of ROS determinations for this and other purposes.
Existing nuclear magnetic resonance logging devices can determine ROS and other formation characteristics, but with limitations. More than twenty years ago, the assignee hereof commercialized a logging device called the "Nuclear Magnetism Tool". This tool has a large coil to which direct current is applied. The resulting static field polarizes nuclear spins (protons) in a region approximately one foot in radius centered on the mandrel on which the coil is mounted. The current is then switched off, and the nuclear spins in the polarized region precess in the earth's magnetic field, inducing a voltage in the coil. The proton signal is selected using a bandpass filter. The length of time the direct current is passed through the coil before being turned off is systematically varied. By correlating the amplitude of the received voltage with the length of time the direct current is left on immediately prior to the measurement, the nuclear magnetic resonance properties of the earth material of the formation can be determined. This type of tool can make a measurement of ROS. Because nuclear magnetic resonance logging devices have relatively shallow depths of investigation, the probe volume is typically flooded to residual oil saturation by invasion of drilling mud filtrate N. R. Morrow, I. Chatzis, D. Siegel, J. J. Taber, "Corrections to In-Situ Measurements of Residual Oil for Flushing at the Wellbore", Journal of Canadian Petroleum Technology vol. 24, no. 4, p. 54 (July 1985)!. After a well is drilled, a water solution of manganese complexed with ethylene diamine tetraacetate (Mn-EDTA) is introduced into the wellbore and allowed to penetrate the formation. The paramagnetic manganese ion reduces the NMR relaxation time of the water fraction to less than 25 ms, so the water fraction becomes invisible to the NMR logging device. Because the Mn-EDTA is not soluble in oil, the nuclear magnetic resonance signal from the oil is unaffected, so the amplitude of the signal detected by the logging device is approximately proportional to the ROS.
When an NMR logging device is employed in the described fashion, there can be uncertainty as to whether the MN-EDTA solution has sufficiently saturated the formation throughout the volume of investigation of the logging device which extends several inches into the formation. If the Mn-EDTA is not present in sufficient concentration, ROS will be overestimated. Therefore, in employing this type of logging approach, users go to extraordinary lengths to ensure that the pore water is thoroughly saturated with Mn-EDTA G. Setser, M. R. Williams, "Measurement of Remaining Oil Saturation in Northern Michigan Using Nuclear Magnetism Log Data and Pressure Core", Society of Petroleum Engineers, Paper 14276 (1985)!. The procedure typically involves running several borehole logging runs before and after application of Mn-EDTA. Also, because manganese is not used in the original drilling process, it is necessary to ream the hole to introduce it. These steps involve considerable expense.
In recent years, improved nuclear magnetic well logging devices and techniques have been proposed and/or developed. A summary of these devices and techniques is presented in U.S. Pat. No. 5,023,551, assigned to the same assignee as the present application. The referenced '551 Patent also reviews basic principles of NMR and NMR logging, and part of that review will be summarized. Reference can be made to the '551 Patent for further detail.
NMR has been a common laboratory technique for forty years, and a theoretical description is available in Abragam, Principles of Nuclear Magnetism, Clarendon Press (Oxford, 1961), and Farrar and Becker, Pulse and Fourier Transform NMR, Academic Press (New York 1971). NMR is based on the fact that the nuclei of many elements have angular momentum ("spin") and a magnetic moment. The nuclear spins align themselves along an externally applied static magnetic field. This equilibrium situation can be disturbed by a pulse of an oscillating magnetic field, which tips the spins away from the static field direction. The angle through which the spins are tipped is under the control of the experimenter, as explained below.
After tipping, two things occur simultaneously. First, the spins precess around the static field at a particular frequency (i.e. the Larmor frequency), given by ωo =γBo where Bo is the strength of the static field and γ is the gyromagnetic ratio, a nuclear constant. Second, the spins return to the equilibrium direction according to a decay time known as the "spin lattice relaxation time" or T1. For hydrogen nuclei, γ/2π=4258 Hz/Gauss, so for a static field of 235 Gauss, the frequency of precession is 1 MHz. T1 is controlled totally by the molecular environment and is typically ten to one thousand milliseconds for water and oil in rocks.
Also associated with the spin of molecular nuclei is a second relaxation time known as the "spin-spin relaxation time", or T2. At the end of a ninety degree tipping pulse, all the spins are pointed in a common direction perpendicular to the static field, and they all precess at the Larmor frequency. However, because of small inhomogeneities in the static field due to imperfect instrumentation or microscopic material heterogeneities, each nuclear spin precesses at a slightly different rate. Hence, after a time that is long compared to the precession period, but shorter than T1, the spins will no longer be precessing in unison. When this "dephasing" is due to static field inhomogeneity of the apparatus, the dephasing is called T2*. When it is due to properties of the material, the dephasing time is called T2. T2 and T2* can be measured independently. For liquids in rocks, T2 is approximately two-thirds of T1.
As aforementioned, the parameters T1 and T2 are sensitive to molecular environment. For example, T2 can be several seconds in an unconfined low viscosity liquid such as water, while it can be as short as ten microseconds in a solid. Liquids confined in the pores of rocks present an intermediate case with T2 in the range of tens to hundreds of milliseconds, depending on pore size and fluid viscosity.
In the basic NMR measurement, a pulse of oscillating field is applied to the sample to tip the spins of the nuclei in the sample. The angle (in radians) through which the spins are tipped is given by the equation
θ=γB1 t.sub.ρ /2 (1)
where γ is the gyromagnetic ratio, B1 is the linearly polarized oscillating field strength, and t.sub.ρ is the duration of the pulse. Tipping pulses of ninety and one hundred and eighty degrees are the most common.
The precessing spins are detected by voltage induced in a coil. Only that component of the nuclear magnetization that is precessing in the plane perpendicular to the static field can be sensed by the coil. Hence, a signal will be generated after a ninety degree tipping pulse but not after a one hundred eighty degree tipping pulse. In fact, after a one hundred eighty degree tipping pulse, the spins do not precess at all, but just slowly return along the Bo axis to the equilibrium direction.
In measuring the spin-lattice relaxation time T1, many different techniques are known both in the material science arts and in the medical arts. The "inversion recovery" technique suggests that after the nuclei have aligned themselves along the static magnetic field, a one hundred eighty degree pulse is applied to reverse the direction of the spins. Over time, the spins decay toward their equilibrium direction according to T1, but no measurement is yet made as the one hundred eighty degree pulse does not induce a signal in the coil. Before the decay is complete, however, it is interrupted by a ninety degree pulse which rotates the spins into the measurement plane (i.e., induces a signal in the coil). However, the measurable signal lasts only as long as the spins precess in unison. As they dephase, the net magnetization decreases, even if all the spins remain in the transverse plane. Therefore, the signal decays exponentially with time constant T2*, also known as the "free induction decay". Fortunately, the information of interest is the amplitude of the signal immediately after the ninety degree "read out" pulse. This amplitude depends on the "recovery time" (τ) between the original one hundred eighty degree pulse and the ninety degree pulse. Following a determination of amplitude, the spin system is permitted to completely relax back to equilibrium, and the pulse sequence is then repeated preferably numerous times with different recovery times. The detected amplitudes are then plotted against τ with the decay typically being expressed as a single exponential.
The inversion recovery technique for measuring T1 has been used in laboratories for about forty years. It is very time consuming, and therefore undesirable for well logging and other material property investigations. To overcome some of the shortcomings of inversion recovery, other techniques such as preparation recovery, steady state, and magnetization conserving techniques have been developed. Reference can be made to the above noted U.S. Pat. No. 5,023,551.
While many different methods for measuring T1 have been developed, a single standard known as the CPMG sequence (Carr-Purcell-Meiboom-Gill) for measuring T2 has evolved. In solids, where T2 is very short, T2 can be determined from the decay of the free induction signal. In liquids, where T2*<<T2, the free induction decay becomes a measurement of the apparatus-induced inhomogeneities. To measure the true T2 in such situations, it is necessary to cancel the effect of the apparatus-induced inhomogeneities. To accomplish the same, a series of pulses is applied. First a ninety degree pulse causes the spins to start precessing. Then a one hundred eighty degree pulse is applied to keep the spins in the measurement plane, but to cause the spins which are dispersing in the transverse plane to reverse direction and to refocus. By repeatedly reversing the spins using one hundred eighty degree pulses, a series of "spin echoes" appear. This succession of one hundred eighty degree pulses after an initial ninety degree pulse is the Carr-Purcell sequence which measures the irreversible dephasing (i.e., T2) due to material properties.
While the Carr-Purcell sequence would appear to provide a solution to eliminating apparatus induced inhomogeneities, it was found by Meiboom and Gill that if the one hundred eighty degree pulses in the Carr-Purcell sequence were slightly misset, the transverse magnetization would steadily be rotated out of the transverse plane. As a result, substantial errors would enter the T2 determination. Meiboom and Gill devised a modification to the Carr-Purcell pulse sequence such that after the spins are tipped by ninety degrees and start to dephase, the carrier of the one hundred eighty degree pulses is phase shifted relative to the carrier of the ninety degree pulse. As a result, any error that occurs during an even pulse of the CPMG sequence is cancelled out by an opposing error in the odd pulse.
The above-referenced U.S. Pat. No. 5,023,551 disclosed, inter alia, an improved pulse sequence for use in NMR borehole logging, called Fast Inversion Recovery/CPMG (FIR/CPMG). This pulse sequence can be set forth in NMR notation as
Wi -180-τi -90-(tcp -180-tcp -echo)j !i
where j=1,2, . . . J, where J is the number of echoes collected in the CPMG sequence and is typically on the order of one or two hundred, but always greater than ten; i=1,2, . . . I, where I is the number of recovery times; Wi are waiting times; τi are recovery times, and tcp is the Carr-Purcell spacing. Measurements of the signals induced in the formation as a result of the magnetic fields are made of each of a predetermined number of echoes. In an embodiment disclosed in the '551 Patent, each echo measurement of the CPMG sequence is a measurement of the integrated amplitude of the echo, rather than a measurement of the greatest amplitude of the received echo. Determinations of T1, T2, and amplitude (Mo) are then made from the measurements. From one or more of the T1, T2, and amplitude determinations, formation parameters such as porosity and permeability can be derived using relationships known in the art.
A technique of interest herein is pulsed field gradient (PFG) nuclear magnetic resonance, which has long been used in the laboratory to measure the diffusion coefficient of materials see e.g., E. O. Stejskal, J. E. Tanner, "Spin Diffusion Measurements: Spin Echoes in the Presence of a Time-Dependent Field Gradient", Journal of Chemical Physics, 42,288 (1965)!. The PFG pulse sequence is illustrated in FIG. 1. After the nuclear spins are polarized by the static field, a 90° RF pulse tips them into the transverse plane. Then a pulse of magnetic field gradient is applied. A 180° RF pulse followed by a second gradient pulse refocuses the spins into an echo. The echo amplitude can be written
A(2τ)=A(0)exp -γ2 Dδ2 g2 (Δ-δ/3)!(3)
A(2τ) is the echo amplitude where 2τ is the time of echo formation, A(0) is proportional to the number of spins resonated, D is the diffusion coefficient, γ is the gyromagnetic ratio of the resonated nucleus, g is the strength of the applied magnetic field gradient, δ is the gradient pulse width, and Δ is the time between the leading edges of the two gradient pulses, as shown in FIG. 1. FIG. 2, which plots phase versus time, illustrates how the pulsed field gradient sequence is useful in measuring the diffusion coefficient of a substance in the gradient field. Consider the difference in behavior between a spin that has moved and one which has not. After application of the 90 degree pulse and prior to the first gradient pulse, both spins exhibit a relatively slow phase change (dephasing) due to background gradient. The gradient pulse results in rapid phase change, and then, after the gradient pulse, slow phase change continues. The 180 degree pulse reverses the spins, and they now begin to change phase in the opposite sense (rephase); slowly, except during application of the second gradient pulse. For a spin that has not moved significantly, the phase will return to zero (its original phase) and contribute to the spin echo. However, a spin that has moved (diffused) to a significantly different position in the gradient field will not rephase to zero, and will not contribute to the echo.
From the viewpoint of nuclear magnetic resonance, the most important difference between oil and water is viscosity. The viscosity of many crude oils is at least an order of magnitude greater than water. The viscosity is related to the diffusion coefficient through the relation
D=c(N/V)1/3 kT/η (4)
where D is the self-diffusion coefficient of a pure fluid, k is Boltzmann's constant, T is the absolute temperature, (N/V) is the molecular number density, η is the viscosity, and c is a constant whose numerical value is around 1/6 for a wide range of liquids. J. C. M. Li, P. Chang, "Self-Diffusion Coefficient and Viscosity in Liquids", Journal of Chemical Physics 23, 518 (1955)!. Where two fluids present in the NMR investigation region have different diffusion coefficients, the signals from the two fluids add in the spin echo. Thus, if there is a large difference in diffusion coefficients, it is possible to find a combination of pulse sequence parameters for which the contribution of the rapidly diffusing liquid is eliminated from the echo signal while the contribution from the slowly diffusing liquid is essentially unaffected. This method has been used to suppress the signal of a solvent so that the NMR spectrum of a more slowly diffusing solute can be studied P. C. M. van Zijl, C. T. W. Moonen, "Complete Water Suppression for Solutions of Large Molecules Based on Diffusional Differences Between Solute and Solvent", Journal of Magnetic Resonance 87, 18 (1990)!. In a technique hereof the relatively large difference between the diffusion coefficients of oil and water is exploited in distinguishing the components of the NMR signal from oil and from water.
In accordance with an embodiment of the invention, there is provided a nuclear magnetic resonance method for determining the volume of hydrocarbons in a region of investigation of earth formations surrounding a borehole, comprising the following steps: a) providing a logging device that is moveable through the borehole; b) generating, from the logging device, a static magnetic field in the region of investigation of the formations; c) generating, from the logging device, a sequence of magnetic field pulses in the region of investigation of the formations, and detecting the nuclear magnetic resonance spin echo strength from the formations; d) varying a parameter of the sequence of magnetic field pulses that causes said echo strength to vary with the diffusion times of fluid species in the formations, and repeating step (c) a plurality of times to detect the echo strength of further nuclear magnetic spin echoes from the formations; and e) determining the volume of hydrocarbons in the region of investigation of the earth formations from the detected spin echo strengths. In a preferred embodiment of the invention, the step (e) comprises determining the volume of hydrocarbons from echo strengths associated with diffusion times longer than a pre-determined diffusion time. In an embodiment of the invention, the step (c) comprises generating a sequence of RF magnetic field pulses and magnetic field gradient pulses. In this embodiment, the step of generating pulses of magnetic field gradient in the formations comprises applying a direct current to an electrically conductive loop in the logging device.
In one embodiment hereof, the sequence of pulses and echoes of step (c) is
90-tg -pfg-tr -180-tg -pfg-tr -echo1 - τi -180-τi -echoi+l !i,
where pfg is a pulsed field gradient, tg is the time between the 90 degree pulse and the pulsed field gradient, tr is the time between the pulsed field gradient and the 180 degree reversing pulse, echo1 is the first spin echo, τi are Carr-Purcell spacings for the subsequent 180 degree reversing pulses and echoes (i=1,2, . . . ).
In another embodiment hereof, the sequence of pulses and echoes is:
90-t1 -90-ts -90-t1 -echo1 - τi -180-τi -echoi+1 !i
where t1 is the time between the initial 90 degree pulse and the second 90 degree pulse, ts is the time between the second and third 90 degree pulses, and τi are the Carr-Purcell spacings (times) for the subsequent 180 degree reversing pulses and echoes (i=1,2 . . . ); d) repeating step (c) a plurality of times with different times t1 or ts ; and e) computing the volume of oil in said region of investigation from the spin echoes detected in steps (c) and (d).
Further features and advantages of the invention will become more readily apparent from the following detailed description when taken in conjunction with the accompanying drawings.
FIG. 1 is a diagram illustrating a prior art pulsed field gradient pulse sequence for nuclear magnetic resonance.
FIG. 2 is a diagram of phase versus time for a moved (diffused) spin and a motionless spin, subjected to the pulse sequence of FIG. 1.
FIG. 3 is a diagram, partially in block form, of a nuclear magnetic resonance logging apparatus which can be used in practicing an embodiment of the invention.
FIG. 4 is a cross-sectional diagram of a portion of the FIG. 3 nuclear magnetic resonance logging apparatus.
FIG. 5 is a perspective view of the RF antenna of the FIG. 3 nuclear magnetic resonance logging apparatus.
FIG. 6 is a cross-sectional view of the RF antenna of FIG. 5, as taken through a section defined by the arrows 6--6 in FIG. 5.
FIG. 7 is a diagram of an embodiment of a nuclear magnetic resonance pulse sequence used in an embodiment of the invention.
FIG. 8A is a diagram illustrating the type of nuclear magnetic resonance data, as a function of a pulsed gradient magnetic field variable, that is useful in understanding an embodiment of the invention.
FIG. 8B shows curves that are useful in understanding FIG. 8A.
FIG. 9 is a perspective diagram of a portion of a logging device that includes an antenna for producing a magnetic field gradient in accordance with an embodiment of the invention.
FIG. 10 is a perspective diagram of a portion of a logging device that includes an antenna for producing a magnetic field gradient in accordance with another embodiment of the invention.
FIG. 11 is a perspective view of the shell of an RF antenna that has been modified so that its structure can also be used as an antenna for producing a magnetic field gradient.
FIG. 12 is a cross-sectional view of the antenna shell of the FIG. 11 embodiment.
FIG. 13 is a block diagram of an embodiment of the electronics which can be used in practicing an embodiment of the invention.
FIG. 14 is a flow diagram of a routine for controlling a processor subsystem in accordance with an embodiment of the invention.
FIG. 15 is a flow diagram of a routine represented by a block of the FIG. 14 routine, for applying a pulse sequence and storing echo responses.
FIG. 16A and 16B are diagrams of further nuclear magnetic resonance pulse sequences.
FIG. 17 is a flow diagram of a routine that can be used for controlling a processor subsystem in accordance with a further embodiment of the invention.
FIG. 18 is a flow diagram of a routine represented by a block of the FIG. 17 routine, for applying a pulse sequence and storing echo responses.
FIG. 19 is a diagram of a nuclear magnetic resonance pulse sequence that can be used in a further embodiment of the invention.
FIG. 20 is a diagram of a Hahn nuclear magnetic resonance pulse sequence.
FIG. 21 shows a CPMG sequence modified to determine diffusion.
FIG. 22 shows Stimulated Echo (STE)
FIG. 23 illustrates repeated refocussing of the stimulated echo (first echo).
FIG. 24 shows a scope trace of STE and refocussed echoes.
FIG. 25 illustrates the approximately quadratic field strength of regions of investigation of borehole NMR tools.
FIG. 26 shows amplitude of stimulated spin echo in arbitrary units as a function of encoding time.
Referring to FIG. 3, there is shown an apparatus for investigating subsurface formations 31 traversed by a borehole 32, which can be used in practicing embodiments of the invention. The borehole 32 is typically, although not necessarily, filled with a drilling fluid or mud which contains finely divided solids in suspension, and mudcake 16 is shown on the walls of the borehole.
An investigating apparatus or logging device 30 is suspended in the borehole 32 on an armored cable 33, the length of which substantially determines the relative depth of the device 30. The cable length is controlled by suitable means at the surface such as a drum and winch mechanism (not shown). Surface equipment, represented at 7, can be of conventional type, and can include a processor subsystem and communicates with the downhole equipment. As described in the U.S. Pat. No. 5,055,788, the device 30 has a face 14 shaped to intimately contact the borehole wall, with minimal gaps or standoff, and a retractable arm 15 which can be activated to press the body of the tool 13 against the borehole wall during a logging run, with the face 14 pressed against the wall's surface. Although the tool 13 is shown as a single body, the tool may alternatively comprise separate components such as a cartridge, sonde or skid, and the tool may be combinable with other logging tools. Also, while a wireline is illustrated, alternative forms of physical support and communicating link can be used, for example in a measurement while drilling system.
The logging device 30 includes a magnet array 17 and an RF antenna 18 positioned between the array 17 and the wall engaging face 14. Magnet array 17 produces a static magnetic field B0 in regions surrounding the tool 13. The antenna 18 produces, at selected times, an oscillating magnetic field B1 which is focussed into formation 31, and is superposed on the static field Bo within those parts of formation opposite the face 14. The "volume of investigation" of the tool, shown in dotted lines in FIG. 3, is a vertically elongated region directly in front of tool face 14 in which there is a point at which the magnetic field produced by the magnet array 17 has a spatial gradient which is approximately zero. As described in the referenced Patent, the tool 13 can make measurements by magnetically tipping the nuclear spins of particles in formation 31 with a pulse of oscillating field B1, and then detecting the precession of the tipped particles in the static, homogeneous field Bo within the volume of investigation over a period of time.
FIG. 4 shows a magnet array 17 of the type disclosed in an embodiment of the referenced '788 Patent. The magnet array includes two permanent magnets 24 and 26, which are mounted generally parallel to each other within a metal alloy body 27. The body 27 should be of a material having low magnetic permeability, so as to not interfere with the static magnetic field. Magnets 24 and 26 are slabs which are elongated in the longitudinal direction of the borehole. The magnetic poles of each magnet are not on the smallest faces of the slab, commonly viewed as the ends of a bar magnet: instead, the poles appear on the two opposing edges of the slab magnet and point to the left and right, respectively, in the Figure. Therefore, within the formation 31, the magnetic field Bo surrounding the magnets remains fairly constant along the longitudinal direction of the borehole axis. In the illustration of FIG. 4, magnets 24, 26 are symmetrically mounted in the two sides of the body 27 with the north poles facing the same direction, that is, the direction of the face 14 of the tool. One or more further permanent magnets can be used.
As described in the referenced '788 Patent, the metal body 27 has, on the front face 14 thereof, a semi-cylindrically shaped cavity or slot 28 which faces formations engaged by the face 14. The cavity 28 is adapted for receiving an RF antenna 18 that is shown in FIGS. 5-6. The antenna 18 is positioned outside of the metal body 27 (FIG. 4) of the tool, and is thereby shielded from electromagnetic communication with regions of the borehole which lie behind the body 27, or regions of other formations in directions intercepted by the body 27. Antenna 18 is thus responsive only to magnetic fields originating in front of the wall engaging face 14, e.g., fields originating in the formation 31 in the vicinity of the antennas 18. In a disclosed embodiment of the referenced Patent, the body 27 is made of metal alloy sheathing, rigidly attached to interior metal bracing, which envelops most components of the tool other than the antenna 18, including the circuitry, the magnet array 17, and the hydraulics system of the arm 15. The Patent points out that the body 27 can alternatively be constructed of other materials, so long as the overall structure is sufficiently strong and the magnetic field of the magnet array 17 can penetrate the body and enter the adjoining formation 31.
In the referenced '788 Patent, the antenna 18 is used both as an RF transmitter to produce an oscillating magnetic field in the formations, and as a receiving antenna to detect coherent magnetic signals emanating from precessing protons (spins) after the oscillating field is terminated. The antenna, which has a body 29 and an elongated center probe 42, across which signals are applied and detected, serves effectively as a current loop which produces an oscillating magnetic field B1 (see FIG. 6) within the volume of investigation that is perpendicular to the static magnetic field, B0 (which is radial in the volume of investigation). The body 29 is trough-shaped and has end plates 40, 41 with the center conductor or probe 42 extending from one end plate 40 to the other end plate 41, parallel to and centered in the semi-cylindrical trough 29. The U.S. Pat. No. 5,153,514 discloses that the trough antenna, which can be filled with a ferrite, can have an inner conductive shell that is separated from a steel body by a rubber layer, which suppresses magnetoacoustic ringing. Before further describing the structure of antennas used herein, some of the underlying theory will be treated.
When two fluids present in the NMR investigation region have different diffusion coefficients, the signals from the two fluids add in the spin echo. Thus, if there is a large difference in diffusion coefficients, it is possible to find a combination of pulse sequence parameters for which the contribution of the rapidly diffusing liquid is eliminated from the echo signal while the contribution from the slowly diffusing liquid is essentially unaffected. In a technique hereof, the relatively large difference between the diffusion coefficients of hydrocarbon and water (oil and water, or gas and water) is exploited in distinguishing the components of the NMR signal from hydrocarbon and from water.
In an embodiment hereof, a pulse sequence that includes pulsed field gradients is used, and an example thereof is illustrated in FIG. 7. Reference can also be made to U.S. Pat. No. 5,212,447 where it is proposed to use a pulse sequence with pulsed field gradients in determining a diffusion coefficient.! After the 90 degree tipping pulse and a gradient pulse, a 180 degree phase reversing pulse is applied (at time t) followed by another gradient pulse. The spin echo (after the rephasing to zero) occurs at time 2t. The 180 degree phase reversing pulses are then applied, in CPMG fashion, every 2τ with spin echoes of successively lower amplitude occurring a time τ after each phase reversing pulse. The first three phase reversing pulses and the first three spin echoes are illustrated in the Figure, it being understood that there may be many more of each. In NMR notation, the pulse sequence is
90-tg -pfg-tr -180-tg -pfg-tr -echo1 - τi -180-τi -echoi+1 !i (5)
where tg is the time between the 90 degree pulse and the pulsed field gradient (about t-Δ/2, in this example), pfg is the pulsed field gradient, tr is the time between the pulsed field gradient and the 180 degree reversing pulse (about Δ/2, in this example), echo1 is the first spin echo, and τi are Carr-Purcell spacings (times) for the subsequent 180 degree reversing pulses and echoes (i=1,2, . . . ). The sequence will be referred to as a "PFG-CPMG" sequence. The initial echo strength can be represented as
M=Mo exp -γ2 g2 δ2 12 !(6)
where γ is the gyromagnetic ratio, g is the pfg amplitude, δ is the pfg pulse width, and l2 is the mean square diffusion length
l2 ≈DΔ(Δ>>δ/3) (7)
where D is the diffusion coefficient and Δ is the time between the gradient pulses. Therefore, equation (6) can be written as
M=Mo exp -γ2 g2 δ2 DΔ!(8)
As previously noted, where two fluids are present which have different diffusion coefficients, the signals from the fluids add in the spin echo. Therefore, for oil and water, equation (8) would be
M=Voil exp -γ2 δ2 g2 Doil Δ!+Vw exp -γ2 δ2 g2 Dw Δ!(9)
where Voil is the volume fraction of oil, Vw is the volume fraction of water, Doil is the diffusion coefficient of oil, and Dw is the diffusion coefficient of water. However, in "Determining Oil and Water Saturations Without Determining Diffusion Coefficients," below, solutions to Voil and Vw are discussed without using diffusion coefficients, which are uncertain, at best.
The values of initial echo strength M as a function of a gradient pulse parameter (again, g, δ, and/or Δ), can also be used to obtain Voil by using values of M that are affected almost exclusively by the oil component of the investigated formation region (and to the exclusion of the water component of the investigated formation region). As described above, a distinction can be based on the difference in diffusion coefficients between oil and water, so by selecting gradient pulse parameters that reduce or eliminate the signal from the more rapidly diffusing liquid (the water), one can use the signal from the oil to determine Voil. FIG. 8A shows the shape of a typical data curve, with the points representing the initial echo strength M(g) for different values of gradient pulse amplitude, g, for a situation where oil and water are present in a sandstone. FIG. 8B illustrates why the Figure A curve has its characteristic shape, with two portions of different slope. In FIG. 8B, the higher slope curve, in dashed line, represents the M(g) vs. g characteristic for a substance (e.g. water) having a relatively high diffusion coefficient, and the lower slope curve, in broken dashed line, represents the M(g) vs. g characteristic for a substance (e.g. oil) having a relatively low diffusion coefficient. The total signal resulting from both substances will add, and result in a characteristic shape like that of FIG. 8A. Referring again to FIG. 8A, by using pulse sequences with values of g that result in data points on the shallower sloped portion of the curve (where the water component has little or no effect), the curve can be extrapolated, as shown in dotted line, to obtain a total signal value representative of the oil-filled porosity or the volume of oil, Voil. At least two points, and preferably more, on the shallower sloped portion of the curve, can be used to establish the extrapolation line. The intercept of the line provides the volume of oil, and the slopes of the two portions respectively are proportional to the diffusion coefficients of the oil and the total of the liquids. The constant of proportionality may be unknown.
While FIGS. 8A and 8B illustrate curves for oil and water, it will be understood that, for example, water and gas will have curves with quite different slopes (in such case, a higher slope for gas, in which molecules diffuse much faster than in a liquid) and could be distinguished in similar manner. A three component (oil, gas, water) model could also be utilized.
Accordingly, the total signal resulting from two or more substances can be used to determine the volume of hydrocarbons in an earth formation. The total signal can be used in at least two ways: to directly indicate the volume of hydrocarbons as a log of the formation, and used in characterizing another parameter of the formation. Specifically, the signal could be used in interpretation techniques to more accurately indicate porosity, for example.
Referring to FIGS. 9-11, there are shown embodiments of an antenna in a logging device that can be used to produce pulse field gradients. The logging device shown in FIG. 9 is of the general type described in U.S. Pat. Nos. 5,055,788 and 5,153,514, and illustrated in part in FIGS. 4-6. In FIGS. 9 and 10, the housing for the permanent magnets is shown at 27, the permanent magnets are shown at 24, 26, and the front face of RF antenna (18) is represented shaded at 518. In the embodiment of FIG. 9, the antenna for generating the magnetic field gradient is a coil 550 (of one turn or a plurality of turns), and the coil 550 is coupled with a current source power supply, as described in conjunction with FIG. 13. The antenna 550 can be wound in an insulating material (such as epoxy) and sandwiched between the front face of the magnet housing and a wear plate (not shown), which is formed of a non-magnetic non-metallic material, such as a heavy plastic, that is mounted on the front face of the logging device. A wear plate is described, for example, in U.S. Pat. No. 5,153,514. In operation, when a DC current is applied to the antenna 550, the magnetic field lines will be generally as indicated by dashed lines in the Figure. If the longitudinal axis of the borehole is considered as the z direction, and the direction the tool is facing is considered as the y direction (see coordinate axes in FIG. 9), then, in the investigation region, the gradient of the magnetic field for the FIG. 9 configuration will be dBy/dy; that is, the gradient will be in the radial direction with respect to the borehole axis. The static field direction in the investigation region is the y direction, and the RF field direction in the investigation region is the x direction.!
In the embodiment of FIG. 10, two coils, 560 and 570 (each of one or more turns), are provided, and the coils are energized so that currents circulate therein in opposite directions, thereby resulting in magnetic field lines as represented in the Figure in dashed line. In this case, the magnetic field gradient in the investigation region will be dBy/dz. Again, these coils can be wound in an insulating medium that is sandwiched between the face of the logging device and a wear plate. It will be understood that the coil or coils could also be oriented to provide a magnetic field gradient dBy/dx in the investigation region.
FIGS. 11 and 12 illustrate another embodiment of an antenna for producing the gradient magnetic field. In this embodiment, a modification of the RF antenna described in U.S. Pat. Nos. 5,055,788 and 5,153,514 is utilized so that the modified antenna can also produce the magnetic field gradient. FIG. 11 illustrates the modified RF antenna which is generally in the shape of a trough and, as disclosed in the referenced Patents, can be filled with a ferrite (not shown). As also described in the '514 Patent, the body of the antenna can have an inner conductive shell separated from the body (which can be formed of steel) by a rubber layer for suppression of magnetoacoustic ringing. In accordance with the modification hereof, the conductive trough is longitudinally split, as illustrated in FIG. 11, and the current source is applied across the two halves 1131, 1132. The other end of the split antenna is shorted, as represented by conductive end plate 1135. The current, which is illustrated in the cross-section of FIG. 12 as flowing into one of the halves of the body and out of the other, results in the type of magnetic field shown in dashed line. In this case, the magnetic field gradient in the investigation region is like that of the FIG. 9 embodiment, and can be designated dBy/dy.
FIG. 13 is a block diagram of electronics that can be utilized in embodiments of the present invention. The electronics can be located, for example, in the body 27 or in a separate cartridge or sonde. The electronics associated with the RF operation of the RF antenna 18 can be of the type described in U.S. Pat. No. 5,055,788, and reference can be made to the '788 Patent for further details. An oscillator circuit 1305 is provided, and the output thereof is coupled to a phase shifter 1310, which is under control of a timing circuit (represented at block 1390) which, in turn, is controlled by a processor subsystem, represented at 1350, to control the phase of signals output therefrom so that, in known manner, the signals used in forming the pulses such as the 90 degree and 180 degree pulses, can be generated. The processor subsystem 1350 may comprise, for example, a microprocessor with associated memory, clock, and input/output circuits (not separately shown). The output of phase shifter 1310 is coupled to an amplitude modulator 1315, which also operates under control of the timing circuit 1390 (which is, in turn, under control of the processor subsystem 1350). The amplitude modulator forms the pulse shapes as described in further detail in the U.S. Pat. No. 5,055,788. The output of amplitude modulator 1315 is coupled to a power amplifier 1325, the output of which is, in turn, coupled with the RF antenna 18. As described in the above-referenced U.S. Pat. No. 5,055,788, a Q-switch circuit 1340 is provided, and is also coupled with the RF antenna, to quickly damp any ringing. The RF antenna 18 is also coupled with receiving circuitry that includes a duplexer 1360 that operates to protect the receiver circuit from the high power signals during transmission. The output of the duplexer 1360 is coupled with a detector circuit 1365 that receives a reference signal which can be derived from the oscillator circuit, as is again described in the referenced '788 Patent. The detector circuit 1365 can be suitably gated. The output of the detector circuit 1365 is coupled to a digitizer circuit 1375 which is under control of the timing circuit, and the output of the digitizer circuit 1375 is coupled with the processor subsystem 1350.
FIG. 13 also shows a further portion of the electronics which is utilized in producing the gradient pulses in accordance with an embodiment hereof. A direct current source power supply 1381 is provided, and may be of the type sold by Techron Company of Elkhart, Ind. The output level of the current source is controlled by processor subsystem 1350. The output of the DC current source 1381 is coupled to the gradient antenna, e.g. an antenna as shown in FIGS. 9-11, via a switching circuit 1383 that is controlled by the processor subsystem 1350, via timing circuit 1390. In operation, the processor subsystem 1350, programmed in accordance with the flow diagrams described below, controls the transmitter circuitry of the electronics to apply the transmitter signals to the RF antenna, to apply gradient pulses to the PFG antenna(s), and to activate the receiver portion of the electronics to detect and process the spin echoes received by the RF antenna.
Referring to FIG. 14, there is shown a routine for controlling a processor, such as the processor subsystem of FIG. 13 and/or, in part, an uphole processor in equipment 7, in accordance with an embodiment of the invention. In the routine of FIG. 14, one or more parameters of the gradient pulses (for example, the amplitude g and/or the gradient pulse duration δ) is varied in order to affect the rapidity of dephasing and rephasing (e.g. of diffusing spins). The time integral of the gradient pulse is approximately gδ, and it is this product that is varied in the present embodiment (by varying g, δ, or both). The block 1420 represents the initializing of gδ to a starting value. The pulse sequence of FIG. 7 is then applied, and the detected spin echoes stored, as represented by the block 1430. The routine represented by the block 1430 is described in conjunction with FIG. 15, and reference will be made for the moment to the FIG. 15. As shown therein, the block 1505 represents application of a 90 degree pulse, using control of the transmitting circuitry as previously described in conjunction with FIG. 13. After a wait of tg (block 1510), the first gradient pulse is applied, using processor control of the circuitry of FIG. 13 to energize the gradient coil (as represented by the block 1525). After a time tr (block 1530) a 180 degree reversing pulse is applied (again, using the transmitter circuitry of FIG. (13), as represented by the block 1540). After a wait of tg (block 1545), the second gradient pulse is applied (block 1550). The circuitry of FIG. 13 is then switched to a receiving mode (as represented by the block 1560) for detection and storage of the spin echo. As described, for example in the above-referenced U.S. Pat. No. 5,023,551, one can detect and store the peak of the spin echo, the area under the spin echo, or both. A reversing pulse index is initialized, as represented by the block 1565. Inquiry is then made (decision block 1570) as to whether the last 180 degree reversing pulse has been transmitted; i.e., whether the reversing pulse index has reached its maximum value. If not, the reversing pulse index is incremented (block 1575) and, after a wait time τ (block 1580--the total wait time of 2τ, since previous 180 degree pulse), the next 180 degree reversing pulse is applied (block 1585), followed by switching to the receiving mode for detection and storage of the next spin echo (block 1590). The decision block 1570 is then reentered, and the loop 1595 is continued until the predetermined number of 180 degree reversing pulse transmissions and spin echo detections have been implemented, whereupon the block 1440 of FIG. 14 is entered.
Referring again to FIG. 14, the block 1440 represents the determination and storage of the first spin echo magnitude (or area, as the case may be) from the stored sequence of spin echoes. Briefly, the amplitude of the first echo is determined by generating the exponential components (Laplace transform coefficients) of the spin echoes and then summing the coefficients to effectively extrapolate to obtain the first echo amplitude. This type of extrapolation is well known, and reference can be made, for example, to Butler et al., "Estimating solutions Of First Kind Integral Equations With Non-Negative Constraints And Optimal Smoothing", SIAM Journal of Numerical Analysis, Vol. 18, No. 3, June, 1981. The decision block 1450 is then entered, and inquiry is made as to whether the last variation of the gradient pulse area (gδ) has been reached. If not, gδ is incremented (block 1460), such as by incrementing g, and the block 1430 is reentered. The loop 1465 then continues until the last value of gδ has been reached, whereupon the block 1480 is entered. This block represents the determination of Voil from the initial spin echo amplitudes (or areas) previously determined in terms of the diagram of FIG. 8, by using values of g that result in data points on the curve, and then interpolating to obtain Voil. Also, as previously described, values for Voil, Vw, can be obtained by simultaneous solution for several values of M as a function of g. Processing can then be implemented for the next depth level.
A so-called "stimulated echo" sequence can be used in practicing an embodiment of the invention. Stimulated echo has been used in laboratory measurement of materials for many years. See, for example, Tanner, "Use Of The Stimulated Echo In NMR Diffusion Studies", Journal Of Chemical Physics, Vol. 52, No. 5, March 1970, and, more recently, Latour et al., "Improved PFG Stimulated Echo Method For The Measurement Of Diffusion In Inhomogeneous Fields", Journal Of Magnetic Resonance, Ser. B. Vol. 101, February, 1993, and Latour et al., "Time-Dependent Diffusion Coefficient of Fluids in Porous Media as a Probe of Surface-to-Volume Ratio", Journal of Magnetic Resonance, Series A, Vol. 101, February, 1993. See also U.S. Pat. No. 5,212,447, where it is proposed to use a stimulated echo sequence in determining a diffusion coefficient.! An example of a stimulated echo pulse sequence, as in Latour et al., is shown in FIG. 16A. The typical 90 degree pulse is applied to tip the spins into the transverse plane. The spins can now be affected by magnetization decay; e.g., from small gradients in the static magnetic field. After a time t1 (which is typically much smaller than diffusion time to be measured) another 90 degree pulse is then used to orient the spins in a direction parallel to the static field. The spins can be thought of as being "stored", since they are in the B0 direction, so they do not precess, and are not affected by the magnetic field gradient. After a period of time, tS, another 90 degree pulse is provided, and this tips the spins back into the transverse plane. Now, the spins behave as though the two latest 90 degree pulses had been a 180 degree reversing pulse, and rephasing begins, so that an echo is obtained after a time which corresponds to the time (t1) between the first two 90 degree pulses. Since diffusion was not suspended during the time that the spins were "stored", the echo will be responsive to diffusion effects. Accordingly, relatively long diffusion times can be measured in the presence of a moderate magnetic field gradient. This is because spins that have not moved a substantial distance in the magnetic field gradient will contribute to the echo (since they will not readily decay during the time they are in "storage"), whereas spins that have moved substantially to a different location in the static magnetic field will not rephase and contribute to the echo. In FIG. 16A, pulsed field gradients are shown during the t1 intervals, as used in Latour et al., supra. FIG. 16B illustrates the sequence without pulsed field gradients, for use when sufficient magnetic field gradient is present without them. After the first spin echo, subsequent cycles of phase reversal and echo can be implemented, as was illustrated above in FIG. 7, and is represented in FIG. 16B.
FIG. 17 is a flow diagram of a routine for controlling the processor subsystem using a pulse sequence of the type shown in FIG. 16B. In the routine of FIG. 17, the time t1 or the time tS (or both) is varied in order to vary the effective diffusion time. The block 1710 represents the initializing of the times t1 and tS. The pulse sequence of FIG. 16B is then applied, and the detected spin echoes stored, as represented by the block 1720. The routine represented by the block 1720 is described in conjunction with FIG. 18, and reference will be made for the moment to the FIG. 18. As shown therein, the block 1805 represents application of a 90 degree pulse, using control of the transmitting circuitry as previously described in conjunction with FIG. 13. After a wait of t1 (block 1810), another 90 degree pulse is applied, as represented by the block 1815. Then, after a wait tS (block 1820), another 90 degree pulse is applied, as represented by the block 1830. The circuitry of FIG. 13 is then switched to a receiving mode (as represented by the block 1850) for detection and storage of the spin echo. As noted above, one can detect and store the peak of the spin echo, the area under the spin echo, or both. If only a single echo is used, return to block 1730 of FIG. 17 can then be implemented. If further spin echoes are to be produced and detected in the pulse sequence, a reversing pulse index is initialized, as represented by the block 1865. Inquiry is then made (decision block 1870) as to whether the last 180 degree reversing pulse has been transmitted; i.e., whether the reversing pulse index has reached its maximum value. If not, the reversing pulse index is incremented (block 1875) and, after a wait time τ (block 1880), a 180 degree reversing pulse is applied (block 1885), followed by switching to the receiving mode for detection and storage of the next spin echo (block 1890). The decision block 1870 is then re-entered, and the loop 1895 is continued until the predetermined number of 180 degree reversing pulse transmissions and spin echo detections have been implemented, whereupon the block 1730 of FIG. 17 is entered.
Referring again to FIG. 17, the block 1730 represents the determination and storage of the amplitude or area (as the case may be) of the echo response from the one or more spin echoes, as above. The decision block 1740 is then entered, and inquiry is made as to whether the last value of t1 or tS has been reached. If not, t1 or tS is incremented (block 1760), and the block 1720 is reentered. The loop 1770 then continues until the last value of t1 or tS has been reached, whereupon the block 1750 is entered. As described above, this block represents the determination of Voil from the echo amplitude response previously determined. Again, for example, simultaneous equations or an extrapolation technique can be used. In FIG. 8A, the horizontal axis can be diffusion time. As also previously described, values for Vw and can also be obtained. Processing can then be implemented for the next depth level.
In employing the techniques described in conjunction with FIGS. 16-18, it will be understood that other types of logging devices, for example the type described in copending U.S. patent application Ser. No. 08/041,643, assigned to the same assignee as the present invention, can be used, if desired, and can provide relatively larger constant magnetic field gradients.
In the foregoing techniques, a magnetic field gradient is used to advantage in obtaining spin echoes that contain useful information about diffusion properties of the formation fluids. In addition to the gradients that can be applied by the apparatus, there are magnetic field gradients within porous media. The magnitude of these gradients is approximately ##EQU1## where Δ.sub.χv is the difference in volumetric magnetic susceptibility between the grain material and the pore fluid, R is a size characterizing the pore geometry, and Bo is the magnitude of the static field applied by the apparatus L. E. Drain, "The Broadening of Magnetic Resonance Lines Due to Field Inhomogeneities in Powdered Samples", Proceedings of the Physical society, 80,1380 (1962); J. A. Glasel, K. H. Lee, "On the Interpretation of Water Nuclear Magnetic Resonance Relaxation Times in Heterogeneous System", Journal of American Chemical Society 96,970 (1974)!. In rocks these gradients can be roughly a few G/cm when the static magnetic field is around 500 Gauss. Internal gradients contribute to the attenuation of signals from rapidly diffusing species E. O. Stejskal, J. E. Tanner, "Spin Diffusion Measurements: Spin Echoes in the Presence of a Time-Dependent Field Gradient", Journal of Chemical Physics, 42,288 (1965)!. For the Stejskal-Tanner sequence, the echo strength for a species with a diffusion coefficient D is ##EQU2## In this equation Go is the internal gradient, which in a simple model is given by equation (10), t1 is the time between the initial 90 pulse and the start of the first gradient pulse, and t2 is the time from the end of the second gradient pulse to the center of the echo. Other symbols were defined in connection with equation (3).
If the internal gradients are large enough, pulse gradients are not needed for diffusion editing. A Hahn spin echo sequence E. L. Hahn, "Spin Echoes", Physical Review 80,580 (1950)! can be used. This sequence is shown in FIG. 20. The echo is attenuated according to ##EQU3## so incrementing τ will quickly destroy the signal from a rapidly diffusing species, if the internal gradient is sufficiently large.
If internal gradients in formations being logged are substantial, a technique of using a pulse sequence with pulsed field gradients of alternating polarity can be utilized to reduce or substantially eliminate the effect of internal gradients in the formations. In the FIG. 13 embodiment, the switching circuit 1383 can be operative, under control of processor subsystem 1350, to reverse the polarity of the current applied to the coil or coils of FIGS. 9-12 for certain gradient field pulses of the pulse sequence. A pulse sequence can be used that is of the type described in Latour et al., "Time-Dependent Diffusion Coefficient Of Fluids In Porous Media As A Probe Of Surface-To-Volume Ratio", Journal Of Magnetic Resonance, Ser. A 101, 342-346, February, 1993. Reference can also be made to Latour et al., "Improved Pulse Field Gradient Stimulated Echo Method For The Measurement Of Diffusion In Inhomogeneous Fields", Journal Of Magnetic Resonance, Ser. B. Vol. 101, 72-77, February, 1993, and to R. M. Cotts et al., "Pulsed Field Gradient Stimulated Echo Methods For Improved NMR Diffusion Measurements In Heterogeneous Systems", Journal Of Magnetic Resonance, Vol. 83, 252, June, 1989. As described in these publications, when using a pulse sequence that has 180 degree pulses followed by gradient pulses that alternate in polarity, the effects of internal gradients in laboratory samples tend to cancel out. In the present embodiment, the effects of internal gradients in earth formations are reduced or eliminated by employing the gradient field coils to implement pulses of alternating polarities as part of the type of pulse sequence illustrated in FIG. 19. In the example of FIG. 19, the time between the initial 90 degree tipping pulse and the first 180 degree reversing pulse is τ, and the time between subsequent 180 degree pulses is 2τ. The gradient pulses are evenly spaced between the 180 degree pulses, and have alternating polarities. In this example, after the initial sequence (which is illustrated as including two pairs of gradient pulses, but may include several more or many more), a 90 degree pulse is used, in the manner of the FIG. 16 example, to achieve a "stimulated echo" type of operation, where the spins are "stored" in the direction of the static field for a time ts. After a further 90 degree pulse to take the spins out of the "stored" state, the initial sequence is repeated, and the spin echo, at a time τ after the final 180 degree pulse, is received and processed. If desired, a CPMG sequence may be appended to this pulse sequence to repeatedly refocus the first echo. In this manner, a more precise estimate of the amplitude of the first echo may be obtained, as described above.
The invention can be implemented with an NMR wireline tool suspended from a cable into a borehole. See U.S. Pat. No. 5,055,787 to Kleinberg et al., for example. The invention can also be implemented with an NMR logging-while-drilling tool mounted to a drill string in a borehole. See U.S. patent application Ser. No. 07/922,254 to Kleinberg, et al., for example. With either type of tool, the pulse sequence can be altered as describe in this application. Also, each tool could house two sets of magnets, RF coil, etc. for separately performing NMR measurements, each set operating at a pulse sequence which is different from the other. Additionally, two tools of one type could be connected in series such that one tool in the series operates at a pulse sequence which is different from that of the other tool.
One of the most important potential applications of nuclear magnetic resonance (NMR) logging is the determination of the volumes of gas, oil, and water in the flushed zone described above.
U.S. Pat. No. 5,212,447 to Paltiel alludes to determining various petrophysical quantities from NMR data. All these methods require the determination of a diffusion coefficient. The computation of a diffusion coefficient is problematic for a number of reasons. In order to compute a diffusion coefficient from static field gradient NMR data, it is necessary to know the strength of the static magnetic field gradient. In some borehole logging tools, permanent magnets create a field that is very nonuniform in space and has no single well-defined magnetic field gradient. In order to compute a diffusion coefficient from pulsed field gradient data, the time integrals of the magnetic field gradient pulses are required. Knowledge of the time integrals may require one or more of the following: components and circuits which are very stable with respect to variations of temperature, extra calibration procedures, or extra sensors to measure the time integral of the gradient pulses during borehole measurements. These requirements increase the cost of the tool, or increase the time required for downhole measurements.
The inventors, however, have discovered how to avoid these uncertainties. It is not necessary to find the diffusion coefficient in order to determine gas, oil, or water saturations from nuclear magnetic resonance logging tool data. The method described here requires that at least one pulse sequence parameter be increased by known multiplicative factors during the measurement cycle. However, no absolute knowledge of pulse sequence or tool parameters is required.
In the following development, it is shown how to solve for the quantities Voil, Vw, and Vg, the volumes of oil, water and gas, respectively, in the formation. These quantities may also be expressed in terms of porosity and fluid saturations:
Voil =αΦSo HIo (13)
Vw =αΦSw HIw (14)
Vg =αΦSg HIg (15)
α is a tool calibration constant, Φ is the porosity of the formation, So, Sw and Sg are the saturations of oil, water and gas, respectively, and HIo, HIw and HIg are the hydrogen indices of oil, water and gas, respectively.
Discussed above are a number of ways to use contrasting diffusions to quantify oil, water, and gas in the zone of investigation of NMR tools. In the following, the determination of Voil and Vw is demonstrated. The same technique can be used for any mixture of two fluids in the earth formation, for example water/gas or oil/gas.
A representative equation is:
M1 =Voil exp(-γ2 δ2 g2 Doil Δ)+Vw exp(-γ2 δ2 g2 Dw Δ)(16)
where M1 is the signal level resulting from the application of a pulse sequence with parameters δ, g, and Δ; the symbols are explained earlier. For concision this may be rewritten as
M1 =Voil A+Vw B (17)
Consider a second measurement, in which one of the following changes is made: (1) δ is increased by a factor of √2; (2) g is increased by a factor of √2; or (3) Δ is increased by a factor of 2. Then the result of this second measurement is
M2 =Voil A2 +Vw B2 (18)
Combining equations (17) and (18), and adopting the definition M0 =Voil +Vw, we can express A and B in terms of the measurements M1 and M2, and the petrophysical quantities of interest Voil and Vw : ##EQU4## Note that the Equations (19a) and (19b) utilize opposite signs between the two terms. The choice is made such that A,B>0, and depends of the ratio Voil /Vw.
Consider now a third measurement, in which a number m, not equal to 1 or 2, is selected, and one of the following changes is made relative to the values used in the first measurement: (1) δ is increased by a factor of √m; (2) g is increased by a factor of √m; or (3) Δ is increased by a factor of m. Then the result of this third measurement is
M3 =Voil Am +Vw Bm (20)
Similarly a fourth measurement is made, in which a number n, not equal to 1, 2, or m, is selected and one of the following changes is made relative to the values used in the first measurement: (1) δ is increased by a factor of √n; (2) g is increased by a factor of √n; or (3) Δ is increased by a factor of n. Then the result of this fourth measurement is
M4 =Voil An +Vw Bn (21)
Equations (20) and (21) are two equations in two unknowns, Voil and Vw. In general, these equations cannot be solved analytically. However, they can be solved numerically.
The above procedure can be simplified if other borehole measurements are made on the same formation. If the temperature, pressure and fluid compositions are known, the hydrogen indices can be readily computed see e.g. Schlumberger Log Interpretation Principles/Applications (1987)!. If the porosity is known from another measurement, the following relationship is useful:
So +Sw +Sg =1 (22)
Thus in the example where only oil and water are present (Sg =0) the above equations can be considerably simplified, and only three measurements are required.
The presence of three fluids complicates the above analysis. The mth measurement can be described by the equation
Mm =Voil Am +Vw Bm +Vg Cm(23)
To solve for Voil, Vw and Vg in the absence of other log information, it is necessary to make six measurements. Each measurement is characterized by an unique factor m, which need not be an integer, such that one of the following is true: (1) δ is increased by a factor of √m; (2) g is increased by a factor of √m; or (3) Δ is increased by a factor of m. The result is solvable using six equations and six unknowns, A,B,C, Voil, Vw, and Vg. Only the last three are of interest, as A, B, and C are functions of tool and pulse sequence parameters. The tool and pulse sequence parameters may be difficult or expensive to determine.
We present pulsed NMR sequences that collect information faster than previously known measurement protocols. We show that neither the ability to pulse the field gradient nor its spatial uniformity are essential.
Nuclear magnetic resonance with pulsed field gradients, discussed above, is a technique for observing molecular diffusion. The basic principle of this measurement is as follows: The rate of nuclear precession is proportional to the strength of the applied static field. If the static field is position dependent, the rate of precession, hence the phase angle of precession after a duration δ, is position dependent, i.e., the position of a spin is encoded in its phase angle of precession. Suppose we invert (multiply by -1) the phase of precession by an RF pulse sequence (such as in the Hahn spin echo) and let the spins precess in the same inhomogeneous field for the same duration δ, thus reversing the phase encoding. At the end of this period, had the spins remained stationary, the phase of each spin would be the same. This is the principle of the Hahn spin echo. A nuclear spin in motion observes different field strengths (precession rates) during forward and reverse phase encoding. Therefore, its phase of precession attains a shift that is proportional to the inner product of the field gradient and the displacement of the spin between the forward and reverse phase encoding. Bulk motion of the sample leads to a phase shift in the total NMR signal; however, molecular diffusion causes each nuclear spin to have a random phase, which leads to a reduction in the signal amplitude.
NMR measurements are notorious for their low signal to noise ratios. Diffusion measurements might further aggravate the situation since diffusion is observed via a reduction in the NMR signal. It is therefore imperative to develop pulse sequences that gather diffusion information as efficiently as possible.
The pulse sequences we shall describe are relevant to an NMR apparatus that can subject the specimen to pulses of RF magnetic field perpendicular to a strong static magnetic field with a spatial gradient. Again, the frequency of the RF magnetic field is γBo, where Bo is the magnitude of the static magnetic field and γ=2.675×104 radian/second/Gauss is the gyromagnetic ratio for protons. For example, an apparatus relevant to this application is described in U.S. Pat. Nos. 5,055,787 and 5,055,788, 8 Oct. 1991. The static field gradient produced by this apparatus is nonuniform in the specimen. The nonuniformity of the gradient introduces a surmountable complication in the NMR measurement. The ability to pulse or otherwise modulate the gradient of the static field is not necessary.
In our notation, a 90°x RF pulse effects a 90° nutation about the x axis in the rotating frame Farrar and Becker, Fourier Transform and Pulse NMR, Academic Press!. This means, the RF pulse is such that γB1 τ/2=π/2 where t and B1 are the duration and amplitude of the RF magnetic field in the specimen. The subscript of 90°x refers to the phase shift applied to the RF pulse: the subscripts x, y, and -x correspond to 0, 90°, and 180° phase shifts. The phase shift determines the direction of the RF pulse in the rotating frame. In the figures, we schematically denote the RF pulses by vertical lines. They are typically 1 to 100 μs-wide pulses containing several cycles of the RF field.
Phase Alternation: The first pulse in the sequences shown in FIGS. 21 to 23, 90°.sub.±x, indicates that the entire sequence is to be performed first by a 90°x pulse and then with a 90°-x pulse. All measurements performed in the latter sequence are to be subtracted from their counterparts in the former sequence. This operation (phase alternation) cancels instrument base line and most of the noise and spurious spin echoes produced by the RF pulses A. Sezginer, R. L. Kleinberg, M. Fukuhara, and L. L. Latour, `Very Rapid Simultaneous Measurement of Nuclear Magnetic Resonance Spin-Lattice Relaxation Time and Spin-Spin Relaxation Time,` Journal of Magnetic Resonance, Vol. 92, p. 504-527, 1991; R. L. Kleinberg, A. Sezginer, and M. Fukuhara, `Nuclear Magnetic Resonance Pulse Sequences for use with Borehole Logging Tools,` U.S. Pat. No. 5,023,551, Jun. 11, 1991!.
The apparent relaxation rate observed by the Carr-Purcell-Meiboom-Gill (CPMG) sequence in bulk liquids has two components: spin-spin relaxation and diffusion. ##EQU5## D is the bulk diffusion coefficient, g=∇∥B0 ∥ is the gradient of the static field intensity, and 2τ is the echo spacing. (For simplicity, we assume that D, g, T2 have single values; they are not distributed over a range of values).
The stimulated spin echo (see FIG. 22) has become the standard sequence in diffusion studies because it allows the observation of diffusion over relatively long times (on the order of the spin-lattice relaxation time, T1) E. L. Hahn, `Spin Echoes,` Physical Review, Vol. 80, No. 4, p. 580-594, 1950!. The amplitude of the STE is: ##EQU6## The amplitude of a stimulated spin echo is only one half of the amplitude of a Hahn echo. However, during the diffusion time Δ, the position information is stored in the amplitude of the longitudinal magnetization; therefore, the position encoding effected by the first two pulses can be retained for as long as a few times T1.
We refocus the stimulated spin echo several hundred times as in the CPMG sequence, thus collecting more signal (see FIG. 23). FIG. 23 illustrates that the stimulated echo (the first echo) can be repeatedly refocused as in the CPMG sequence. The echo spacing is so short that the relaxation rate is not affected by diffusion. This sequence is well suited for observing diffusion over long times (on the order of T1) and collecting signal efficiently. It can be used with either pulsed (dashed) or static (solid) field gradient. The echo spacing (2tcp) is so short that ##EQU7## diffusion does not affect the apparent spin-spin relaxation rate. The amplitude of the n th echo is: ##EQU8## The parameters A, T1, and T2 can be obtained by standard NMR measurements.
FIG. 24 shows an actual scope trace of the stimulated echo (STE) and the first three refocussed echoes (0.4 ms spacing). Typically, hundred to thousand echoes are formed before the signal decays by spin-spin relaxation. The variations in the echo amplitudes are due to field inhomogeneities and they are predictable. The scope trace is suppressed in between the echoes, where 180°y pulses are applied.
In conventional imaging equipment, spatially uniform field gradients are produced by gradient coils driven by pulse generators and amplifiers. In borehole logging, it is not possible to produce either a constant static field or a constant field gradient in the earth formation. Furthermore, typical static field gradients of borehole NMR instruments are comparable to the gradients that can be produced by downhole coils. Using the static gradients leads to a simpler instrument. Below, we shall concentrate on using the static gradients although the pulse sequences described above are equally applicable to pulsed field gradients.
The borehole NMR tool of Schlumberger has approximately quadratic static field strengths in their regions of investigation (see FIG. 25).
B0 ≈constan t+g1 x2 +g2 y2(27)
The magnitude of the field gradient ∇∥B0 ∥ is approximately proportional to the distance from a point at the center of the sensitive zone (when g1 2 ≈g2 2). We shall show below that nonuniform field gradients complicate but do not disable the interpretation of field gradient experiments. In typical borehole logging, water and oil molecules cannot diffuse far enough to sample the quadratic nature of the static field. For example, at 80° C., in one second, a water molecule diffuses 200 μm in bulk water and even less in a water saturated rock. FIG. 25 illustrates that the field strength of some borehole NMR instruments are approximately quadratic in their region of investigation.(Eg., Bo =410G+(y2 -x2) (20 G/cm2). This provides nonuniform gradients.
For typical parameters (tcp =150 μs, T2 =150 ms), refocusing the stimulated echo repeatedly yields a factor of 250 (24 dB) improvement in the signal to noise power ratio. The gradient of the field shown in FIG. 25 does not change significantly in 200 μm:
g2 =4g1 2 x2 +4g2 2 y2 (28)
(-g1=g2=20 G/cm2 in FIG. 25) To calculate the total signal via (24), each infinitesimal part of the sample can be treated as a sample with a constant gradient: ##EQU9## For a rectangular sensitive region of size -a1, a1 !× -a2, a2 !, the stimulated echo amplitude is ##EQU10## where the amplitude A is proportional to the length of the region of investigation in the third dimension and to the porosity of the specimen, and erf(z)=(1/√π)∫-z z e-t2 dt is the error function.
This calculation is shown in FIG. 26, (solid curves) which plots stimulated spin echoes as a function of the phase encoding time δ. The circles denote measurements performed on a water sample, and plus signs denote measurements performed on an oil sample S6 viscosity standard, Cannon Instrument Company, State College, Pa. 16804!, and the solid lines are calculated according to the theory described here. Two regimes are visible in the figure: in the flat region at the left, the random phase shifts due to diffusion are neglible compared to one revolution (2π) and the stimulated echo amplitude reduces to ##EQU11## On the right-most part of FIG. 26, the random phase shifts due to diffusion are large compared to one revolution. In this regime, the value of the error function in (30) approaches unity; thus, the amplitude of the stimulated echo reduces to ##EQU12## Notice that in a quadratically varying magnetic field, the echo amplitude behaves like δ-2 asymptotically, which is less pronouced than the exp(-δ2 . . . ) dependence in a uniform gradient.
This calculation was fitted to measurements performed on oil and water according to the pulse sequence shown in FIG. 23 using 1500 refocused echoes of the stimulated echo. The amplitudes shown in FIG. 26 were calculated by fitting an exponential decay to the train of 1500 spin echoes.
|Cited Patent||Filing date||Publication date||Applicant||Title|
|US3657730 *||Jan 14, 1970||Apr 18, 1972||Shell Oil Co||Method for determining residual hydrocarbons present in a subterranean earth formation|
|US4717877 *||Sep 25, 1986||Jan 5, 1988||Numar Corporation||Nuclear magnetic resonance sensing apparatus and techniques|
|US4719423 *||Aug 13, 1985||Jan 12, 1988||Shell Oil Company||NMR imaging of materials for transport properties|
|US5023551 *||Dec 19, 1989||Jun 11, 1991||Schlumberger-Doll Research||Nuclear magnetic resonance pulse sequences for use with borehole logging tools|
|US5055787 *||Dec 5, 1989||Oct 8, 1991||Schlumberger Technology Corporation||Borehole measurement of NMR characteristics of earth formations|
|US5055788 *||Dec 5, 1989||Oct 8, 1991||Schlumberger Technology Corporation||Borehole measurement of NMR characteristics of earth formations|
|US5153514 *||Feb 19, 1991||Oct 6, 1992||Schlumberger Technology Corp.||Antenna and wear plates for borehole logging apparatus|
|US5212447 *||Dec 3, 1990||May 18, 1993||Numar Corporation||Apparatus and technique for nmr diffusion measurement|
|US5280243 *||Apr 16, 1993||Jan 18, 1994||Numar Corporation||System for logging a well during the drilling thereof|
|US5291137 *||Nov 2, 1992||Mar 1, 1994||Schlumberger Technology Corporation||Processing method and apparatus for processing spin echo in-phase and quadrature amplitudes from a pulsed nuclear magnetism tool and producing new output data to be recorded on an output record|
|US5363041 *||Dec 31, 1992||Nov 8, 1994||Schlumberger Technology Corporation||Determining bound and unbound fluid volumes using nuclear magnetic resonance pulse sequences|
|US5387865 *||Mar 11, 1993||Feb 7, 1995||Exxon Research And Engineering Company||Permeability determination from NMR relaxation measurements for fluids in porous media|
|US5428291 *||Jul 1, 1993||Jun 27, 1995||Exxon Research And Engineering Company||Determination of fluid transport properties in porous media by nuclear magnetic resonance measurements of fluid flow|
|US5486762 *||Aug 15, 1994||Jan 23, 1996||Schlumberger Technology Corp.||Apparatus including multi-wait time pulsed NMR logging method for determining accurate T2-distributions and accurate T1/T2 ratios and generating a more accurate output record using the updated T2-distributions and T1/T2 ratios|
|US5497087 *||Oct 20, 1994||Mar 5, 1996||Shell Oil Company||NMR logging of natural gas reservoirs|
|US5498960 *||Oct 20, 1994||Mar 12, 1996||Shell Oil Company||NMR logging of natural gas in reservoirs|
|EP0489578A1 *||Dec 3, 1991||Jun 10, 1992||Numar Corporation||Apparatus and technique for NMR diffusion measurement|
|EP0581666A2 *||Jul 27, 1993||Feb 2, 1994||Schlumberger Limited||Pulsed nuclear magnetism tool for formation evaluation while drilling|
|GB2279754A *||Title not available|
|1||A. Sezginer, R. L. Kleinberg, M. Fukuhara, and L. L. Latour, "Very Rapid Simultenous Measurement of Nuclear Magnetic Resonance Spin-Lattice Relaxation Time and Spin-Spin Relaxation Time", Journal of Magnetic Resonance, vol. 92, pp. 504-527, 1991.|
|2||*||A. Sezginer, R. L. Kleinberg, M. Fukuhara, and L. L. Latour, Very Rapid Simultenous Measurement of Nuclear Magnetic Resonance Spin Lattice Relaxation Time and Spin Spin Relaxation Time , Journal of Magnetic Resonance, vol. 92, pp. 504 527, 1991.|
|3||*||Abragam, Principles of Nuclear Magnetism, Clarendon Press (Oxford, 1996).|
|4||Butler et al., "Estimating Solutions of First Kind Integral Equations with Non-Negative Constraints and Optimal Smoothing", SIAM Journal of Numcerical Analysis, vol. 18, No. 3, Jun. 1981.|
|5||*||Butler et al., Estimating Solutions of First Kind Integral Equations with Non Negative Constraints and Optimal Smoothing , SIAM Journal of Numcerical Analysis, vol. 18, No. 3, Jun. 1981.|
|6||E. L. Hahn, "Spin Echoes", Physical Review 80,580 (1950).|
|7||*||E. L. Hahn, Spin Echoes , Physical Review 80,580 (1950).|
|8||E. O. Stejskal, J. E. Tanner, "Spin Diffusion Measurements: Spin Echoes in the Presence of a Time-Dependent Field Gradient", Journal of Chemical Physics, 42,288 (1965).|
|9||*||E. O. Stejskal, J. E. Tanner, Spin Diffusion Measurements: Spin Echoes in the Presence of a Time Dependent Field Gradient , Journal of Chemical Physics, 42,288 (1965).|
|10||*||Farrar and Becker, Pulse and Fourier Transform NMR, Academic Press (New York 1971).|
|11||*||G. Setser, M. R. Williams, Measurement of Remaining Oil Saturation in Northern Michigan Using Nuclear Magnetism Log Data and Pressure Core, Society of Petroleum Engineers, Paper 14276 (1985).|
|12||J. A. Glasel, K. H. Lee, "On the Interpretation of Water Nuclear Magnetic Resonance Relaxation Times in Heterogeneous System", Journal of American Chemical Society 96, 970 (1974).|
|13||*||J. A. Glasel, K. H. Lee, On the Interpretation of Water Nuclear Magnetic Resonance Relaxation Times in Heterogeneous System , Journal of American Chemical Society 96, 970 (1974).|
|14||J.C.M. Li, P. Chang, "Self-Diffusion Coefficient and Viscosity in Liquids", Journal of Chemical Physics 23,518 (1955).|
|15||*||J.C.M. Li, P. Chang, Self Diffusion Coefficient and Viscosity in Liquids , Journal of Chemical Physics 23,518 (1955).|
|16||L.E. Drain, "The Broadening of Magnetic Resonance Lines Due to Field Inhomogeneities in Powered Samples", Proceedings of the Physical Society, 80,1380 (1962).|
|17||*||L.E. Drain, The Broadening of Magnetic Resonance Lines Due to Field Inhomogeneities in Powered Samples , Proceedings of the Physical Society, 80,1380 (1962).|
|18||Latour et al., "Time Dependent Diffusion Coefficient of Fluids in Porous Media as a Probe of Surface-to-Volume Ratio", Journal of Magnetic Resonance, Series A, vol. 101, 342-346, Feb. 1993.|
|19||*||Latour et al., Improved PFG Stimulated Echo Method for the Measurement of Diffusion in Inhomogeneous Fields, Journal of Magnetic Resonance, Ser. B., vol. 101, 72 77, Feb. 1993.|
|20||Latour et al., Improved PFG Stimulated Echo Method for the Measurement of Diffusion in Inhomogeneous Fields, Journal of Magnetic Resonance, Ser. B., vol. 101, 72-77, Feb. 1993.|
|21||*||Latour et al., Time Dependent Diffusion Coefficient of Fluids in Porous Media as a Probe of Surface to Volume Ratio , Journal of Magnetic Resonance, Series A, vol. 101, 342 346, Feb. 1993.|
|22||N. R. Morrow, I. Chatzis, D. Siegel, J. J. Taber, "Corrections to In-Situ Measurements of Residual Oil for Flushing at the Wellbore", Journal of Canadian Petroleum Technology, vol. 24, No. 4, p. 54 (Jul. 1985).|
|23||*||N. R. Morrow, I. Chatzis, D. Siegel, J. J. Taber, Corrections to In Situ Measurements of Residual Oil for Flushing at the Wellbore , Journal of Canadian Petroleum Technology, vol. 24, No. 4, p. 54 (Jul. 1985).|
|24||P.C.M. van Zij, C.T.W. Moonen, "Complete Water Suppression for Solutions of Large Molecules Based on Diffisuional Differences Between Solute and Solvent", Journal of Magnetic Resonance 87, 18 (1990).|
|25||*||P.C.M. van Zij, C.T.W. Moonen, Complete Water Suppression for Solutions of Large Molecules Based on Diffisuional Differences Between Solute and Solvent , Journal of Magnetic Resonance 87, 18 (1990).|
|26||R. Akkurt, H. J. Vinegar, P. N. Tutunjian, and A. J. Guillory, "NMR Logging of Natural Gas Reservoirs", SPWLA 36th Annual Logging Symposium, Jun. 26-29, 1995, pp. 1-12.|
|27||*||R. Akkurt, H. J. Vinegar, P. N. Tutunjian, and A. J. Guillory, NMR Logging of Natural Gas Reservoirs , SPWLA 36th Annual Logging Symposium, Jun. 26 29, 1995, pp. 1 12.|
|28||R. M. Cotts et al., "Pulsed Field Gradient Stimulated Echo Methods for IMproved NMR Diffusion Measurements in Heterogeneous Systems" Journal of Magnetic Resonance, vol. 83, 252, Jun. 1989.|
|29||*||R. M. Cotts et al., Pulsed Field Gradient Stimulated Echo Methods for IMproved NMR Diffusion Measurements in Heterogeneous Systems Journal of Magnetic Resonance, vol. 83, 252, Jun. 1989.|
|30||*||Schlumberger Log Interpretation Principles/Applications (1987).|
|31||Tanner, "Use of the Stimulated Echo in NMR Diffusion Studies", Journal of Chemical Physics, vol. 52, No. 5, Mar. 1970.|
|32||*||Tanner, Use of the Stimulated Echo in NMR Diffusion Studies , Journal of Chemical Physics, vol. 52, No. 5, Mar. 1970.|
|Citing Patent||Filing date||Publication date||Applicant||Title|
|US6023164 *||Feb 20, 1998||Feb 8, 2000||Numar Corporation||Eccentric NMR well logging apparatus and method|
|US6051973 *||Dec 23, 1997||Apr 18, 2000||Numar Corporation||Method for formation evaluation while drilling|
|US6097184 *||Dec 31, 1997||Aug 1, 2000||Schlumberger Technology Corporation||Nuclear magnetic resonance well logging to determine gas-filled porosity and oil-filled porosity of earth formations without a constant static magnetic field gradient|
|US6107796 *||Aug 17, 1998||Aug 22, 2000||Numar Corporation||Method and apparatus for differentiating oil based mud filtrate from connate oil|
|US6111408 *||Dec 23, 1997||Aug 29, 2000||Numar Corporation||Nuclear magnetic resonance sensing apparatus and techniques for downhole measurements|
|US6184681||Nov 5, 1998||Feb 6, 2001||Schlumberger Technology Corporation||Apparatus and method for computing a distribution of spin-spin relaxation times|
|US6242913||Apr 18, 2000||Jun 5, 2001||Numar Corporation||Method for formation evaluation while drilling|
|US6246236||Nov 5, 1998||Jun 12, 2001||Schlumberger Technology Corporation||Apparatus and method for obtaining a nuclear magnetic resonance measurement while drilling|
|US6255818||Aug 4, 1999||Jul 3, 2001||Schlumberger Technology Corporation||Method and apparatus for performing magnetic resonance measurements|
|US6255819||Oct 25, 1999||Jul 3, 2001||Halliburton Energy Services, Inc.||System and method for geologically-enhanced magnetic resonance imaging logs|
|US6268726||Jan 15, 1999||Jul 31, 2001||Numar Corporation||Method and apparatus for nuclear magnetic resonance measuring while drilling|
|US6291995||Nov 5, 1998||Sep 18, 2001||Schlumberger Technology Corporation||Apparatus and method for generating a pulse sequence|
|US6316940||Mar 17, 1999||Nov 13, 2001||Numar Corporation||System and method for identification of hydrocarbons using enhanced diffusion|
|US6326784||Nov 5, 1998||Dec 4, 2001||Schlumberger Technology Corporation||Nuclear magnetic resonance logging with azimuthal resolution using gradient coils|
|US6331775 *||Sep 15, 1999||Dec 18, 2001||Baker Hughes Incorporated||Gas zone evaluation by combining dual wait time NMR data with density data|
|US6362619||Jul 31, 2001||Mar 26, 2002||Numar Corporation||Method and apparatus for nuclear magnetic resonance measuring while drilling|
|US6366087||Oct 28, 1999||Apr 2, 2002||George Richard Coates||NMR logging apparatus and methods for fluid typing|
|US6377042||Aug 10, 1999||Apr 23, 2002||Numar Corporation||Method and apparatus for merging of NMR echo trains in the time domain|
|US6429654 *||Sep 28, 2000||Aug 6, 2002||Baker Hughes Incorporated||Nuclear magnetic resonance pulse sequence for improving signal-to-noise ratio|
|US6466013 *||Apr 18, 2000||Oct 15, 2002||Baker Hughes Incorporated||Nuclear magnetic resonance measurements in well logging using an optimized rephasing pulse sequence|
|US6498484 *||Jun 15, 2000||Dec 24, 2002||Schlumberger Technology Corporation||Method for reducing ringing in nuclear magnetic resonance well logging instruments|
|US6512371||Jun 5, 2001||Jan 28, 2003||Halliburton Energy Services, Inc.||System and method for determining oil, water and gas saturations for low-field gradient NMR logging tools|
|US6518755||Apr 10, 2001||Feb 11, 2003||Baker Hughes Incorporated||Measurement technique and apparatus for high-resolution multi-volume NMR well logging|
|US6522136 *||Mar 20, 2000||Feb 18, 2003||Schlumberger Technology Corporation||Well logging technique and apparatus for determining pore characteristics of earth formations using magnetic resonance|
|US6525534||Jun 15, 2001||Feb 25, 2003||Halliburton Energy Services, Inc.||System and methods for NMR signal processing without phase alternated pair stacking|
|US6531868||Jun 5, 2001||Mar 11, 2003||Halliburton Energy Services, Inc.||System and methods for formation evaluation while drilling|
|US6531869 *||Sep 12, 2002||Mar 11, 2003||Schlumberger Technology Corporation||Methods and apparatus for measuring flow velocity in a wellbore using NMR and applications using same|
|US6541969||Dec 14, 2000||Apr 1, 2003||Halliburton Energy Services, Inc.||Method and apparatus for improving the vertical resolution of NMR logs|
|US6545471 *||Dec 11, 2001||Apr 8, 2003||Wollin Ventures, Inc.||Method for resistivity well logging utilizing nuclear magnetic resonance|
|US6570381 *||Mar 25, 1999||May 27, 2003||Schlumberger Technology Corporation||Nuclear magnetic resonance well logging method and apparatus|
|US6570382||Nov 28, 2000||May 27, 2003||Schlumberger Technology Corporation||Nuclear magnetic resonance method and logging apparatus|
|US6577125||Dec 18, 2000||Jun 10, 2003||Halliburton Energy Services, Inc.||Temperature compensated magnetic field apparatus for NMR measurements|
|US6583621||Mar 27, 2002||Jun 24, 2003||Numar Corporation||Method and apparatus for nuclear magnetic resonance measuring while drilling|
|US6586932||Oct 16, 2002||Jul 1, 2003||Schlumberger Technology Corporation||Enhanced performance antennas for NMR logging|
|US6626027||Jun 12, 2001||Sep 30, 2003||Intertech Development Company||Method and apparatus for detecting a gas leak using nuclear magnetic resonance|
|US6630357||Feb 2, 2001||Oct 7, 2003||University Technologies International Inc.||Quantification of bitumen using NMR|
|US6646438 *||Jul 9, 2002||Nov 11, 2003||Baker Hughes Incorporated||NMR data acquisition with multiple interecho spacing|
|US6650114||Jun 28, 2001||Nov 18, 2003||Baker Hughes Incorporated||NMR data acquisition with multiple interecho spacing|
|US6661226||Aug 10, 2000||Dec 9, 2003||Halliburton Energy Services, Inc.||NMR apparatus and methods for measuring volumes of hydrocarbon gas and oil|
|US6686737 *||Oct 12, 2001||Feb 3, 2004||Baker Hughes Incorporated||Amplitude and/or phase modulated NMR pulse sequences|
|US6727696 *||Aug 13, 2001||Apr 27, 2004||Baker Hughes Incorporated||Downhole NMR processing|
|US6774628||Jan 18, 2002||Aug 10, 2004||Schlumberger Technology Corporation||Nuclear magnetic resonance imaging using phase encoding with non-linear gradient fields|
|US6781371||Sep 6, 2002||Aug 24, 2004||Schlumberger Technology Corporation||High vertical resolution antennas for NMR logging|
|US6825657||Oct 30, 2001||Nov 30, 2004||Schlumberger Technology Corporation||Magnetic resonance method for characterizing fluid samples withdrawn from subsurface earth formations|
|US6825658||Apr 2, 2002||Nov 30, 2004||George Richard Coates||NMR logging apparatus and methods for fluid typing|
|US6825659||Jun 24, 2003||Nov 30, 2004||Numar||Method and apparatus for nuclear magnetic resonance measuring while drilling|
|US6833698||May 15, 2002||Dec 21, 2004||Chevrontexaco U.S.A. Inc.||Methods of decoupling diffusion effects from relaxation times to determine properties of porous media containing fluids|
|US6841996||Jan 22, 2003||Jan 11, 2005||Schlumberger Technology Corporation||Nuclear magnetic resonance apparatus and methods for analyzing fluids extracted from earth formation|
|US6850060||Apr 10, 2003||Feb 1, 2005||Schlumberger Technology Corporation||Method and apparatus for rapid characterization of diffusion|
|US6856132||Nov 8, 2002||Feb 15, 2005||Shell Oil Company||Method and apparatus for subterranean formation flow imaging|
|US6883702||Mar 7, 2003||Apr 26, 2005||Schlumberger Technology Corporation||Method and apparatus for NMR measurement of wettability|
|US6891369||Dec 13, 2002||May 10, 2005||Schlumberger Technology Corporation||Nuclear magnetic resonance method and logging apparatus for fluid analysis|
|US6894493||Dec 9, 2003||May 17, 2005||Schlumberger Technology Corporation||Method and apparatus for NMR measurement of magnetic materials|
|US6937013 *||Jun 19, 2002||Aug 30, 2005||Schlumberger Technology Corporation||NMR tool for making formation evaluation measurements using gradient echoes|
|US6937014||Mar 24, 2003||Aug 30, 2005||Chevron U.S.A. Inc.||Method for obtaining multi-dimensional proton density distributions from a system of nuclear spins|
|US6956371||Jan 28, 2003||Oct 18, 2005||Halliburton Energy Services, Inc.||Method and apparatus for detecting diffusion sensitive phases with estimation of residual error in NMR logs|
|US6972564||Nov 5, 2002||Dec 6, 2005||Baker Hughes Incorporated||Objective oriented methods for NMR log acquisitions for estimating earth formation and fluid properties|
|US6984980 *||Apr 26, 2002||Jan 10, 2006||Baker Hughes Incorporated||Method and apparatus for NMR sensor with loop-gap resonator|
|US7053611 *||Jun 4, 2004||May 30, 2006||Schlumberger Technology Corporation||Method and apparatus for using pulsed field gradient NMR measurements to determine fluid properties in a fluid sampling well logging tool|
|US7194698||Mar 25, 2002||Mar 20, 2007||Elliot A. Gottfurcht||Method to advertise and search on television for web content using a simplified interface|
|US7253618||Dec 19, 2006||Aug 7, 2007||Schlumberger Technology Corporation||Method for determining more accurate diffusion coefficient distributions of reservoir fluids using Bi-polar pulsed field gradients|
|US7358725||May 16, 2006||Apr 15, 2008||Baker Hughes Incorporated||Correction of NMR artifacts due to axial motion and spin-lattice relaxation|
|US7412898 *||Jul 28, 2006||Aug 19, 2008||Disney Enterprises, Inc.||Load sensing system including RFID tagged fasteners|
|US7441196||Mar 13, 2006||Oct 21, 2008||Elliot Gottfurcht||Apparatus and method of manipulating a region on a wireless device screen for viewing, zooming and scrolling internet content|
|US7564240 *||Jun 19, 2007||Jul 21, 2009||Schlumberger Technology Corporation||Method and apparatus for measuring free induction decay signal and its application to composition analysis|
|US7703318 *||Jan 19, 2004||Apr 27, 2010||Schlumberger Technology Corporation||Permanently eccentered formation tester|
|US7733086||Oct 30, 2008||Jun 8, 2010||Halliburton Energy Services, Inc.||Systems and methods for deep-looking NMR logging|
|US7741841||Dec 28, 2007||Jun 22, 2010||Schlumberger Technology Corporation||Time-lapsed diffusivity logging for monitoring enhanced oil recovery|
|US7755354||Nov 4, 2008||Jul 13, 2010||Halliburton Energy Services, Inc.||System and methods for T1-based logging|
|US7999542||Apr 30, 2009||Aug 16, 2011||Schlumberger Technology Corporation||Method for determining formation parameter|
|US8191416||Nov 24, 2008||Jun 5, 2012||Schlumberger Technology Corporation||Instrumented formation tester for injecting and monitoring of fluids|
|US8212567 *||Oct 20, 2008||Jul 3, 2012||Baker Hughes Incorporated||Externally mounted band antennae requiring minimal metal cutting on drillstring for reduction of mechanical stresses|
|US8235106||Jan 18, 2010||Aug 7, 2012||Halliburton Energy Services, Inc.||Formation testing and sampling apparatus and methods|
|US8258790 *||Nov 20, 2008||Sep 4, 2012||Baker Hughes Incorporated||Oscillator sensor for determining a property of an earth formation|
|US8370884||Jun 30, 2010||Feb 5, 2013||United Video Properties, Inc.||Interactive television systems with digital video recording and adjustable reminders|
|US8427145||Mar 24, 2010||Apr 23, 2013||Schlumberger Technology Corporation||System and method for emulating nuclear magnetic resonance well logging tool diffusion editing measurements on a bench-top nuclear magnetic resonance spectrometer for laboratory-scale rock core analysis|
|US8452653||Jun 30, 2003||May 28, 2013||Personal Data Network Corporation||Method for providing information and recommendations based on user activity|
|US8522870||Jul 31, 2012||Sep 3, 2013||Halliburton Energy Services, Inc.||Formation testing and sampling apparatus and methods|
|US8680858||Jun 29, 2009||Mar 25, 2014||Schlumberger Technology Corporation||NMR logging of miscible displacement|
|US8686724||Mar 13, 2013||Apr 1, 2014||Schlumberger Technology Corporation||System for emulating nuclear magnetic resonance well logging tool diffusion editing measurements on a bench-top nuclear magnetic resonance spectrometer for laboratory-scale rock core analysis|
|US8799971||May 20, 2011||Aug 5, 2014||United Video Properties, Inc.||Interactive television systems with digital video recording and adjustable reminders|
|US8806546||Apr 19, 2013||Aug 12, 2014||United Video Properties, Inc.||Interactive television systems with digital video recording and adjustable reminders|
|US8816684 *||Oct 28, 2010||Aug 26, 2014||Vista Clara Inc.||Noise canceling in-situ NMR detection|
|US9071872||Jun 24, 2014||Jun 30, 2015||Rovi Guides, Inc.||Interactive television systems with digital video recording and adjustable reminders|
|US9322948||Mar 7, 2012||Apr 26, 2016||Schlumberger Technology Corporation||Method of determining a formation parameter|
|US9369741||Jun 16, 2015||Jun 14, 2016||Rovi Guides, Inc.||Interactive television systems with digital video recording and adjustable reminders|
|US9494665||Jun 23, 2011||Nov 15, 2016||Ramot At Tel-Aviv University Ltd.||Magnetic resonance analysis using a plurality of pairs of bipolar gradient pulses|
|US20020097261 *||Mar 25, 2002||Jul 25, 2002||Gottfurcht Elliot A.||Apparatus and method for simple wide-area network navigation|
|US20030006768 *||Oct 30, 2001||Jan 9, 2003||Schlumberger Technology Corporation||Magnetic resonance method for characterizing fluid samples withdrawn from subsurface earth formations|
|US20030016012 *||Apr 2, 2002||Jan 23, 2003||Coates George Richard||NMR logging apparatus and methods for fluid typing|
|US20030092311 *||Nov 12, 2002||May 15, 2003||France Telecom||Comb and a method for making a branch connection to preexisting cabling|
|US20030107374 *||Nov 5, 2002||Jun 12, 2003||Baker Hughes Incorporated||Objective oriented methods for NMR log acquisitions for estimating earth formation and fluid properties|
|US20030137297 *||Jan 18, 2002||Jul 24, 2003||Krishnamurthy Ganesan||Nuclear magnetic resonance imaging using phase encoding|
|US20030151408 *||Apr 26, 2002||Aug 14, 2003||Baker Hughes Incorporated||Method and apparatus for NMR sensor with loop-gap resonator|
|US20030178994 *||Dec 13, 2002||Sep 25, 2003||Schlumberger Technology Corporation,||Nuclear magnetic resonance method and logging apparatus|
|US20030214287 *||May 15, 2002||Nov 20, 2003||Boqin Sun||Methods of decoupling diffusion effects from relaxation times to determine properties of porous media containing fluids and multi-dimensional representation of those properties|
|US20040124837 *||Jun 24, 2003||Jul 1, 2004||Numar||Method and apparatus for nuclear magnetic resonance measuring while drilling|
|US20040140800 *||Jan 22, 2003||Jul 22, 2004||Schlumberger Technology Corporation||Nuclear magnetic resonance apparatus and methods for analyzing fluids extracted from earth formation|
|US20040150676 *||Jan 21, 2004||Aug 5, 2004||Gottfurcht Elliot A.||Apparatus and method for simple wide-area network navigation|
|US20040150677 *||Jan 21, 2004||Aug 5, 2004||Gottfurcht Elliot A.||Method for navigating web content with a simplified interface using audible commands|
|US20040164735 *||Dec 9, 2003||Aug 26, 2004||Schlumberger Technology Corporation||Method and apparatus for NMR measurement of magnetic materials|
|US20040189296 *||Mar 24, 2003||Sep 30, 2004||Chevron U.S.A. Inc.||Method for obtaining multi-dimensional proton density distributions from a system of nuclear spins|
|US20040199932 *||Jan 21, 2004||Oct 7, 2004||Gottfurcht Elliot A.||Navigating web content with a simplified interface using directional input|
|US20050257610 *||Jul 12, 2005||Nov 24, 2005||Baker Hughes Incorporated||Automatic adjustment of NMR pulse sequence to optimize SNR based on real time analysis|
|US20050270023 *||Jun 4, 2004||Dec 8, 2005||Robert Freedman||Method and apparatus for using pulsed field gradient NMR measurements to determine fluid properties in a fluid sampling well logging tool|
|US20050272158 *||Jan 4, 2005||Dec 8, 2005||Galford James E||Method and apparatus for detecting hydrocarbons with NMR logs in wells drilled with oil-based muds|
|US20060150726 *||Jan 19, 2004||Jul 13, 2006||Aaron Jacobson||Permanently eccentered formation tester|
|US20060273787 *||May 16, 2006||Dec 7, 2006||Baker Hughes Incorporated||Correction of NMR artifacts due to axial motion and spin-lattice relaxation|
|US20080315873 *||Jun 19, 2007||Dec 25, 2008||Schlumberger Technology Corporation||Method and Apparatus for Measuring Free Induction Decay Signal and Its Application to Composition Analysis|
|US20090167302 *||Dec 28, 2007||Jul 2, 2009||John Edwards||Time-Lapsed Diffusivity Logging for Monitoring enhanced Oil Recovery|
|US20100097279 *||Oct 20, 2008||Apr 22, 2010||Baker Hughes Incorporated||Externally mounted band atennae requiring minimal metal cutting on drillstring for reduction of mechanical stresses|
|US20100116494 *||Jan 18, 2010||May 13, 2010||Halliburton Energy Services, Inc.||Formation Testing and Sampling Apparatus and Methods|
|US20100123461 *||Nov 20, 2008||May 20, 2010||Baker Hughes Incorporated||Oscillator sensor for determining a property of an earth formation|
|US20100126717 *||Nov 24, 2008||May 27, 2010||Fikri Kuchuk||Instrumented formation tester for injecting and monitoring of fluids|
|US20100277166 *||Apr 30, 2009||Nov 4, 2010||Schlumberger Technology Corporation||Method for determining formation parameter|
|US20110109311 *||Oct 28, 2010||May 12, 2011||Vista Clara, Inc.||Noise canceling in-situ nmr detection|
|US20110175607 *||Jun 29, 2009||Jul 21, 2011||Schlumberger Technology Corporation||Nmr logging of miscible displacement|
|US20130176026 *||Jan 9, 2013||Jul 11, 2013||Schlumberger Technology Corporation||Magnetic Resonance Imaging Methods|
|US20140203807 *||Mar 26, 2014||Jul 24, 2014||Green Imaging Technologies, Inc.||Methods for performing nmr measurements on porous media|
|CN100401107C||Jun 3, 2005||Jul 9, 2008||施卢默格海外有限公司||Method for using pulsed field gradient NMR measurements to determine fluid properties in a fluid sampling well logging tool|
|CN100409024C||Jun 19, 2003||Aug 6, 2008||施卢默格海外有限公司||NMR spectrum analysis by means of gradient field NMR|
|CN101328805B||Jun 19, 2008||Mar 26, 2014||普拉德研究及开发股份有限公司||Method and apparatus for measuring free induction decay signal and its application in composition analysis|
|CN103221837A *||Jun 23, 2011||Jul 24, 2013||雷蒙特亚特特拉维夫大学有限公司||Magnetic resonance analysis using a plurality of pairs of bipolar gradient pulses|
|EP1305653A1 *||Jun 7, 2001||May 2, 2003||Baker-Hughes Incorporated||Nuclear magnetic resonance tool with active rf spoiler antenna|
|EP1305653B1 *||Jun 7, 2001||Dec 31, 2008||Baker Hughes Incorporated||Nuclear magnetic resonance tool with active rf spoiler antenna|
|EP1397708A2 *||Jun 5, 2002||Mar 17, 2004||Halliburton Energy Services, Inc.||System and method for determining oil, water and gas saturations for low-field gradient nmr logging tools|
|EP1397708A4 *||Jun 5, 2002||Aug 26, 2009||Halliburton Energy Serv Inc||System and method for determining oil, water and gas saturations for low-field gradient nmr logging tools|
|WO2001042817A1 *||Nov 30, 2000||Jun 14, 2001||Schlumberger Limited||Nuclear magnetic resonance method and logging apparatus|
|WO2002099463A2||Jun 5, 2002||Dec 12, 2002||Halliburton Energy Services, Inc.||System and method for determining oil, water and gas saturations for low-field gradient nmr logging tools|
|WO2003003056A1 *||Jun 27, 2002||Jan 9, 2003||Baker Hughes Incorporated||Nmr data acquisition with multiple interecho spacing|
|WO2010059540A2 *||Nov 16, 2009||May 27, 2010||Baker Hughes Incorporated||Oscillator sensor for determining a property of an earth formation|
|WO2010059540A3 *||Nov 16, 2009||Aug 26, 2010||Baker Hughes Incorporated||Oscillator sensor for determining a property of an earth formation|
|WO2011161683A1||Jun 23, 2011||Dec 29, 2011||Ramot At Tel-Aviv University Ltd.||Magnetic resonance analysis using a plurality of pairs of bipolar gradient pulses|
|Jan 24, 2002||FPAY||Fee payment|
Year of fee payment: 4
|Jan 27, 2006||FPAY||Fee payment|
Year of fee payment: 8
|Jan 29, 2010||FPAY||Fee payment|
Year of fee payment: 12