Pulsed nuclear magnetic resonance (NMR) measurements use bursts of electromagnetic fields (also called rf pulses) at a specific frequency (also called the operating frequency or rf frequency) to induce response echoes in the material under measurement. To obtain the maximum response to a particular element, the operating frequency should be at the “Larmor frequency” (also called the “sweet spot” of the instrument). The Larmor frequency is dependent upon local static magnetic fields, which may change during logging. For example, temperature changes or the accumulation of magnetic debris on the tool will affect the strength of the static magnetic field and influence logging measurements. Because the tool is calibrated relative to the Larmor frequency, changes in the static magnetic field must be accounted for, such as by changing the operating frequency of the pulses and the resonance frequency of the antenna.
Accordingly, two frequencies should be tuned for proper instrument operation: (1) the operating frequency should be set to the Larmor frequency taking into account local static fields and (2) the resonance frequency of the antenna should be set to the rf frequency.
Improper tuning of an NMR logging tool affects the porosity calibration in two distinct ways. First, the strength of the rf pulses in the formation is changed, which leads to a reduced precessing magnetization. Second, a given precessing magnetization will induce a smaller signal in the spectrometer.
The conventional method of tuning Schlumberger's CMR™ tool is described in commonly owned U.S. Pat. No. 5,451,873 to Freedman, et al. (the '873 Patent) (incorporated by reference herein in its entirety). In the method of the '873 Patent, the tool is positioned in front of a high porosity zone, preferably as a station stop and the rf (operating) frequency is estimated based on the temperature reading of the sonde. The antenna, which operates like a parallel-resonant LC circuit, is then tuned to this frequency by the so-called “tune word search task” (TWST). The TWST uses a continuous test signal at the operating frequency injected with constant current into the antenna. The tuning capacitor is then varied while monitoring the detected amplitude of a signal injected into a test loop located on the antenna. The value of capacitance that provides the largest amplitude is selected for the best tuning. The TWST has to be done before driving the antenna with the electromagnetic bursts.
With this antenna tuning, the rf frequency is then adjusted to the proper Larmor frequency at the sweet spot with the so-called “Larmor frequency search task” (LFST). According to LFST, the Larmor frequency is determined by taking a set of measurements of the relative echo signal amplitudes at different operating frequencies. This set of measurements is fitted to a predetermined response curve to obtain the frequency maximum amplitude. Essentially, the rf-frequency is varied to find the maximum NMR signal. This is time consuming and may require that the tool be stationary in the well.
The temperature of the magnet is then monitored during logging and the rf frequency is adjusted based on the known temperature coefficient of the magnet material. A Hall probe may also be used to measure the magnetic field strength inside the sonde and detect accumulation of magnetic debris on the magnets.
The LFST cannot be performed during a NMR measurement and, therefore, needs to be performed beforehand. Accordingly, all other factors that affect amplitude need to remain constant during the LFST.
The LFST is a slow procedure because, for every setting of rf frequency, the amplitude has to be measured with a good signal to noise ratio. This is not always practical. A formation of at least 10 pu is required to accomplish the tuning procedure in which the signal strength for seven (7) different rf frequencies is measured and the optimal rf frequency and corresponding antenna tuning is inferred. During this time, the porosity in front of the tool must be constant, requiring a station stop in most cases. However, in many logging environments, station stops are not possible. Further, in fluid sampling tools, the flowing fluid will not likely have a constant hydrogen index.
During a continuous log it is not possible to verify directly that the operating frequency is tracking any changes in the Larmor frequency. Predictable changes occur due to changes in the temperature of the magnets that create the static field. The change of the magnet's field with temperature is a known magnitude and so the operating frequency can be corrected approximately by a measure of the sonde temperature. However, a sudden change in the magnetic field caused by magnetic particles attracted and deposited on the magnets can cause unknown shifts in the Larmor frequency that invalidate the tool calibration. To attempt detection of such shifts the magnetic field is measured in the vicinity of the zone under measurement. But this determination is coarse and does not allow for quantitative correction of the shifts.
Accordingly, there is a need for a continuous measurement of any deviation of the operating frequency from the Larmor frequency during the NMR measurement to allow the operating frequency to track this deviation by means of a feedback loop and so maintain the instrument calibrated throughout the measurement.
There is presented herein an improved method to continuously tune NMR logging tools, which accounts for changes in static magnetic fields.
In accordance with the present invention, the echo signal is analyzed to detect any mistuning of the rf frequency from the Larmor frequency. Likewise, either an independent measurement or the echo signal phase is used to monitor the detuning of the rf frequency with respect to the resonance frequency of the antenna. For the purposes of the present invention, the term “antenna” is the system component that transmits the rf magnetic field and detects the magnetic field due to the NMR precession, such as an NMR coil. The tuning of the antenna and the tuning of the rf frequency with respect to the Larmor frequency may be adjusted, if necessary. This “retuning” can be done continuously while logging to account for any local changes in the static magnetic field, allowing proper tuning of the tool to be maintained while logging and reducing the need for conventional off-line time-consuming tuning procedures. The resulting improvements are a reduction of systematic errors in porosity determination and the optimization of the signal to noise ratio.
The tuning procedure of the present invention is based on information obtained from the NMR signal and may be implemented using the presently used pulse sequences, in particular sequences based on the CPMG sequence. Therefore, the tuning is measured as seen by the spins without requiring a change in the pulse sequence or the introduction of additional signals into the antenna.
Before this method is implemented for a given NMR logging tool, its electronics should preferably be calibrated and characterized during a master calibration phase. For example, temperature dependence of the electronics, such as phase shifts versus temperature, should be identified.
It is noted that the method of the present invention is best used to maintain the tuning of the tool. However, the method may be used to initially calibrate the system. If the method is used to initially calibrate the system, the higher order terms in Equation (5) (discussed below) may be significant. In this case, it may be necessary to iterate the tuning procedure. Alternatively, it might be more efficient to use the conventional procedure to determine the approximate tuning and then implement the method of the present invention to refine and maintain the tuning of the tool.
The present method to measure detuning of the rf frequency from the Larmor frequency by analyzing the echo signal is based on changes in the phase of the echo signal with respect to the rf phase along the echo interval. In a preferred embodiment, the signal is amplified, converted to an intermediate frequency, detected by a phase-sensitive demodulator and integrated over the duration of the echo. A measurement of the signal phase relative to a reference (the instrument's rf frequency) is implicit in the signal demodulation. Therefore, in addition to integrating the whole echo, the interval of duration of the echo is partitioned into several zones (or intervals), the simplest case being the two halves of the echo. Each of the zones is integrated independently, obtaining for each zone the amplitude and a phase value relative to the operating frequency source. Although any measurements of absolute phase are not stable and, therefore, are of little use, the relative phase differences between the integrated zones change with the difference between operating frequency and Larmor frequency. So the deviation from the Larmor frequency correlates nearly linearly with the phase differences measured between zones along the echo.
These phase differences measured are affected in large extent by the antenna detuning (deviation of the resonance frequency of the antenna from the operating frequency). Therefore, the present invention includes means to continuously measure antenna detuning and correct the phase differences for this effect.
Accordingly, in one embodiment of the present invention, a method for tuning a nuclear magnetic resonance (NMR) tool having an operating frequency and equipped with an antenna, is disclosed comprising: (a) transmitting a rf magnetic field to a sample under investigation; (b) receiving an NMR signal from the sample within a detection window; (c) determining mistuning of the antenna relative to the operating frequency; (d) analyzing the received echo signal to determine mistuning of the received signal from the operating frequency. Alternatively, the complex echo shape of the received echo signal may be analyzed. For the purposes of this patent application, “complex echo shape” includes components from both the in-phase and out-of-phase channels. In a second embodiment, the mistuning of the received signals from the operating frequency is determined by analyzing any changes in phase of the echo along the echo signal. Accordingly, signal shape may be quantified by analyzing changes in phase along the echo signal.
The antenna tuning process may be automated by measuring the amplitudes of calibration signals generated at more than one frequency. These amplitudes may be analyzed to determine mistuning of the antenna relative to the operating frequency. In one embodiment, a maximum amplitude may be identified to determine antenna mistuning. The system tuning may be maintained by repeating (a) through (d) while operating the tool and implementing a feedback loop. Further, the effects of the antenna mistuning may be minimized by performing (c) and (d) nearly simultaneously.
The detection window may be partitioned into more than one interval and in-phase and out-of-phase signals detected in each interval. These signals may be analyzed to determine mistuning of the received echo signal from the operating frequency. Any combination of in-phase signals and out-of-phase signals may be analyzed. For example, two in-phase signals from different intervals may be analyzed or two out-of-phase signals from different intervals may be analyzed. Alternatively, in-phase and out-of-phase signals from the same or different intervals may be analyzed. The slope of a signal detected in the out-of-phase channel may be determined at the signal (echo) center. The tuning of the operating frequency may be inferred from the out-of-phase slope. The out-of-phase slope may be approximated by determining the ratio of the out-of-phase portion of the detection window versus the in-phase portion of the detection window according to the ratio defined below. The measurement of may be refined by taking multiple measurements to determine ratio for more than one echo and calculating an average
Mistuning of the antenna may be determined by comparing the signal phase of the received signal to the signal phase of the transmitted signal. In particular, the tuning of the antenna may be inferred based on the difference in phase between the received signal and the transmitted signal such as by: (a) transmitting an rf magnetic field at the operating frequency modulated by an amount approximately equivalent to the width of the antenna resonance; (b) receiving an NMR signal; and (c) comparing the received signal to the operating frequency of the modulated magnetic field. Further, the intensity of sidebands corresponding to the modulated rf magnetic field and the received signal may be measured to assist in this analysis.
Further features and applications of the present invention will become more readily apparent from the figures and detailed description that follows.
FIGS. 6(a) and (b) are graphs showing “density of states” for logging and fluid sampling tools.
In an embodiment of the present invention, mistuning of the rf frequency from the Larmor frequency is measured wherein the main frequency component of the echo signal (Larmor Frequency) is compared with the rf frequency (frequency of the rf signals transmitted by the instrument and the reference of the phase-sensitive echo demodulation) by two or more phase measurements of the echo signal at fixed time intervals along the echo. This principle is illustrated in
In practice, to measure the amplitude of the echoes (the primary data acquired by this instrument) the signal from the antenna is amplified, converted to an intermediate frequency, filtered and digitized by a fast analog to digital converter (ADC) during the time of the echo. Once digitized, the data is processed by a fast digital signal processor (DSP) that performs phase-sensitive demodulation and digital integration, resulting in a set of values for each echo that reflect its amplitude and phase. By partitioning the measurement time of the echo and performing the above described processing to a certain number of sub-periods, it is possible to determine the change in instantaneous phase of the signal along the echo. In the simplest case, the echo signal is digitally integrated over its two halves, and demodulated into the orthogonal R and X signals for each. From these the values of sine and cosine of the phase angle for each half are computed and the phase difference is extracted from these as described above and illustrated below.
It was verified that the increase in phase along the echo (phase of the later part minus phase of the earlier part) increases in proportion to the difference between the operating frequency and the Larmor frequency. The calculation shown in
The phase values from multiple zones along the echo may be calculated in various ways. For practical purposes, a flexible implementation can be limited to two zones as shown in
By way of example, the present invention was implemented using Schlumberger's CMRT-A version of the CMR™ tool. Results of phase difference measurements taken with Schlumberger's CMRT-A tool using the standard water bottle as the measurement sample are shown in
The set of curves in
Another method of the present invention takes advantage of the information that is contained in the overall phase of the signal and in its out-of-phase (or “noise”) channel. The analysis presented below shows that the signal in this channel is antisymmetric with time relative to the echo center. It is shown here that the amplitude of the antisymmetric signal is related to the detuning of the rf frequency with respect to the Larmor frequency ad to the detuning of the antenna resonance frequency with respect to the rf frequency.
Echo Amplitude
The signal generated from any position in space depends on the timing of the pulse sequence and the following two frequencies: ω0≡γ|{right arrow over (B0)}−ωrf, and
Here {right arrow over (B1)}⊥ is the component of {right arrow over (B1)} that is orthogonal to {right arrow over (B0)}, ω0 measures the detuning of the static field from resonance, and ω1 is proportional to the rf field strength.
The asymptotic expression for the transverse magnetization in the rotating at the time of the formation of the echo for the standard CPMG sequence is given by:
where M0 is the thermal equilibrium magnetization, Ω={square root}{square root over (ω02+ω12)}, tπ/2, tπ are the duration of the π/2 and π pulses, respectively, and tE is the echo spacing. This is the asymptotic form for large number of echoes, but it has been shown that in inhomogeneous fields it is approached by the third echo. Note in particular that for all values of (ω0,ω1), {right arrow over (M)}⊥ has only a component along the ŷ axis. The asymptotic echoes of a CPMG pulse sequence in grossly inhomogeneous fields therefore have a pure absorption spectrum.
The detected signal, centered at an echo time t0, is proportional to:
S0(t0+t)=∫∫dω0dω1(ωrf+ω0)ω1F(ω0,ω1)M⊥(ω0,ω1)ejω
Here F(ω0,ω1) is the distribution function of ω0 and ω1 for the particular tool, defined such that F(ω0,ω1)Δω0Δω1 is proportional to the number of spins that are exposed to a static field in the range of ω0±Δω0/2 and an rf field in the range of ω1±Δω1/2. The factor (ωrf+ω0)≅ωrf is due to Faraday detection. The factor ω1 accounts for the antenna efficiency.
Field Profiles
Equation (2) is written in the rotating frame that rotates with the rf frequency of the pulses. This corresponds directly to the demodulated signal that is being observed in the spectrometer. The signal is affected when the rf frequency is changed by Δωrf.
The excitation spectrum M⊥ (ω0,ω1) is centered at the rf frequency of the pulses and is therefore independent of the rf frequency (except weakly through M0). In contrast, the distribution function F (ω0ω1) depends on the static field profile. In the rotating frame, F (ω0,ω1) changes with the offset between the rf frequency and the Larmor frequency in a zone wherein the static magnetic field is influenced (also known as “a sensitive zone”). When the rf frequency is changed by Δωrf, the distribution function should be translated to a new rotating frame, resulting in F(ω0,ω1)→F(ω0−Δωrf,ω1).
There are two extreme cases. In a typical laboratory set-up both the static and rf fields are very homogeneous. In this case, the distribution function is essentially a delta function: F(ω0,ω1)≈δ(ω0−ω0,lab)δ(ω1−ω1,lab). All the spins have the same Larmor frequency. If the rf frequency deviates strongly from this Larmor frequency, all the spins will be excited off-resonance and there will be a large decrease in the detected signal. On the other extreme, if the applied fields have a constant gradient, the distribution function is flat and independent of offset frequency ω0: F(ω0,ω1)≅constant. For any value of rf frequency, they are always roughly the same number of spins on resonance. Changing the rf frequency does not change the detected signals, it only changes the physical place where the signal is generated from.
Some tools, such as Schlumberger's CMR™ tool, are designed somewhere between these two extremes. The field is very inhomogeneous, but there is a saddle point where the fields are more homogeneous. By contrast, some tools, such as Schlumberger's MRX™ and MDT™ tools, are closer to the second extreme. This is illustrated in FIGS. 6(a) and (b) which are plots of the integral of the net distribution function ω1F(ω0,ω1) over the rf field strength ω1: ρ(f0)=∫dω1F(ω0,ω1)ω1. [These plots include the effect of the antenna efficiency.] The functions were calculated based on the measured field maps. The function ρ(f0) can be treated as a density of states. These plots show the projections of ρ(f0)=∫dω1F(ω0,ω1)ω1 of the distribution function ω1F(ω0,ω1) as a function of the Larmor frequency offset
normalized to 1 at f0 equal to 0. For tools with a saddle point (
Echo Shape: rf/Larmor Frequency Tuning
To the first order, the signal of the spin echo in the rotating frame can be written as:
Here ωc is the resonance frequency of the coil, ωrf is the rf frequency of the NMR pulses, and Q1 is the loaded quality factor of the coil. f(ω0) is the spectrum of the signal. The frequency ω0 is the difference between the local Larmor frequency and the rf frequency. The time t is measured with respect to the center between two adjacent π pulses. The spectrum f(ωo) depends on the details of the pulse sequence and the field maps of the logging tool. For the standard CPMG sequence, the analysis shows that the spectrum f(ω0) is purely real after the first few echoes.
The first order result of Equation (3) shows several important features. The finite bandwidth of the antenna delays the echo by a time
This delay is in addition to any other delays that might be intrinsic to the electronics of the tool. With respect to the new echo center, the signal in the in-phase channel is symmetric and in the out-of-phase channel is antisymmetric with time. Further, the overall phase of the signal is only related to the detuning of the coil from the rf frequency, ωrf−ωc. The measured phase can therefore be used to infer the tuning of the antenna.
The slope of the out-of-phase signal at the echo center is directly related to the first moment of the spectrum {overscore (ω)}0, which is a sensitive indicator of the tuning of the rf frequency with respect to the Larmor frequency at the sweet spot of the magnet. As will be discussed below, a ratio that approximates the slope of the out-of-phase signal can be measured more robustly from the signal. is independent of porosity and can be used to indicate the tuning of the rf frequency of the NMR pulses with respect to the Larmor frequency at the sensitive zone of the formation.
The detuning of the rf frequency of the NMR pulses from the Larmor frequency at the sweet spot can be inferred from the shape of the spin echoes. It is important to use the information from both the in-phase and the out-of-phase channel (i.e. the signal in the two channels after the phasing procedure has been performed).
As mentioned above, it is useful to introduce the following ratio
where Sin-phase and Sout-of-phase are the in- and out-of-phase echo shapes, respectively, t0 is the center of the echo, and T is the total duration of the two detection windows.
The ratio of the out-of-phase signal (Sout-of-phase) to the in-phase signal (Sin-phase) can be expressed in terms of an instantaneous signal phase
The ratio can therefore be rewritten as a phase difference, =tan ψ1−tan ψ2≈ψ1−ψ2 where ψ1, ψ2 are the signal phases averaged over the in-phase and out-of-phase channels (detection windows).
is independent of porosity is related to the peak frequency of the spectrum of the measured signal and therefore indicates the tuning of the rf frequency with respect to the Larmor frequency in the formation. In the current data collection and analysis of a tool with a saddle point, the echo signal is integrated over a window of duration T, centered at the echo, in both the in-phase (A) and out-of-phase (B) channels as shown in
The expected shapes of the acquired echoes do not change significantly after the second echo of a CPMG sequence. Therefore, many echoes may be averaged to get a better determination of similar to the procedure commonly performed to determine phase. In addition, the tuning is not expected to change rapidly during logging. This allows the stacking of echoes of many adjacent CPMG sequences, this further decreasing the uncertainty in the determination of
The ratio defined in Equation (4) is just one particularly simple way to obtain a tuning parameter. The location and duration of the integration intervals has to be optimized for the best performance for a given logging tool.
One generalization of the tuning indicator defined above in Equation (4) is to modify the flat top integrations in Equation (4) with windowing functions. There will be a trade-off in better noise performance versus complexity in signal acquisition and processing.
For the NMR tool used, the Larmor frequency offset is preferably less than ±7.5 kHz for an error of less than 3% in porosity. This corresponds to a 1 pu error in a 30 pu formation. In this range,
The exact value of is not necessarily zero, but will depend on the ({right arrow over (B0)},{right arrow over (B1)}) field maps in the sensitive zone of the specific logging tool. For the NMR tool used, is determined to be equal to 0.08.
When the antenna resonance frequency is off-resonance, the spectrum of the detected signal will be somewhat distorted. This will affect the measured value of Close to resonance can be expanded linearly around :
≅+α(fL−frf)+β(fc−frf) (5)
This second term does not cause any significant problems. Once β has been determined for a given logging tool, the second term can be subtracted, because the detuning (fc−frf) can be determined independently from the measurement of the overall phase, as discussed below. In this way, the detuning (fL−frf) can be extracted from the measured value of
Signal Phase: Antenna Tuning
In NMR logging tools, the data is acquired in two channels. It is a standard procedure to rotate the data by an angle φ so that the integrated echoes show a coherent signal in one channel and noise in the other. (A description of this procedure applied to Schlumberger's CMR™ tool can be found in commonly owned U.S. Pat. No. 5,291,137 to Freedman, incorporated by reference herein in its entirety.)
The analysis of the spin dynamics in inhomogeneous fields coupled to a tuned circuit as mentioned above shows that the phase φ of the signal depends predominantly on the detuning of the resonance frequency of the antenna (or coil),fc, from the rf frequency frf:
Here Q1 is the loaded quality factor of the antenna. Note that to the first order, there is no dependence on the Larmor frequency at the sweet spot. Close to resonance, the phase change can be approximated by:
In
using the independently measured value of Q1 equal to 28. As depicted in this graph, the dependence of the phase on rf frequency shows excellent agreement with the theoretical prediction given above. This demonstrates that the measurement of the phase can be used to extract the antenna tuning through Equation (6). The nominal value of the phase on the resonance depends on the exact implementation of the phase cycling. In the above case, it is 180°.
The electronics of a specific NMR tool might exhibit additional phase shifts that should preferably be calibrated (including any potential temperature dependence). Such extra phase shifts will lead to an offset of the phase. However, with a proper calibration, the change in phase can still be used to infer the tuning of the antenna.
A standard electronic way to tune the antenna is the present TWST procedure: a signal of constant amplitude is injected into the calibration loop. The tuning of the. antenna is changed to maximize the detected signal at the output of the spectrometer.
One modified electrical method to determine the resonance frequency of the antenna is to modulate a signal at the rf frequency ωrf of the NMR pulses with a frequency ωmod. Best sensitivity is achieved if ωmod is comparable to
the width of the antenna resonance. The resulting signal, comprising of two frequencies ωrf±ωmod of equal amplitude, is injected into the calibration loop and detected with the spectrometer. When the tool is on-resonance, the intensities of the two detected sidebands, V− and V+ are equal. A potential detuning of the antenna, δω=ωc-ωrf, can be inferred from the ratio of the detected intensities of the two sidebands, V− and V+. To the first order, the detuning is given by:
Equation (7) can be used to estimate the change in tuning capacitance, and if necessary iteration can be used to find the final value.
Another method for measuring antenna tuning is based on time-multiplexed calibrations of tool responses at different frequencies. This measurement does not use the echo signal, it does not depend on the presence or strength of the echo signal and, therefore, it can drive a continuous feedback loop that keeps the antenna tuned. The amplitude calibration of the NMR measurement is performed periodically in-between measuring cycles, during the bursts-free “wait time” allowed for repolarization. It uses a fixed-amplitude test signal at the operating frequency that is a secondary standard derived from a laboratory calibration of the instrument. To measure antenna tuning, two additional calibration measurements are done but with the frequency of the test signal shifted by a small ±deviation around the operating frequency. This provides a three-point amplitude measurement of resonance and a quantitative value of the detuning from resonance.
With a second-order polynomial approximation of the resonance curve, if y1, y2, y3 are the calibration amplitudes at the lower, center and upper frequencies and D is the frequency deviation, the antenna detuning in the vicinity of resonance is proportional to:
x=D*(y1−y3)/abs(2*y2−y1−y3) (8)
The term (2*y2−y1−y3) reflects the curvature of the resonance curve as shown in
A positive curvature value (convex portion of the curve) gives a solid identification of the central zone of the resonance curve. The sign of the term (y1−y3) reflects on what side of the resonance curve the antenna is operating.
With a detuning measurement performed for each CPMG, a feedback loop can adjust the antenna capacitance to maintain the antenna at resonance. By changing the capacitance of the antenna, a feedback loop based on this method can find the optimal point of antenna resonance over the whole range of capacitance/frequency values. The frequency deviation value D can be selected for optimum resolution versus signal-to-noise ratio of the measurement.
The implementation of this measurement in Schlumberger's CMR™ tool follows the timing as shown in
Alternatively, antenna mistuning may be determined by analyzing the overall phases of the detected signals. As shown in
Results of antenna detuning measurements made with three different frequency deviations D of 1.83, 3.66 and 7.32 kHz are shown in
The measuring of Larmor frequency deviation requires that there is sufficient signal amplitude to make measurements with adequate signal-to-noise ratio. This is not always the case during a continuous log in an oil well. A practical system should keep the current method of adjusting the operating frequency with tool temperature, while the feedback loop here described introduces an additional offset into the operating frequency. The feedback loop should respond depending on the quality of the signal present. The stronger the signal, the stronger/faster the loop response should be. In absence of signal, the loop should become dormant and the operating frequency should continue tracking with tool temperature and the last offset value from the loop until sufficient signal is present again for the loop to resume. Alternatively, the echoes may be stacked as described above.
It is noted that the expressions provided above are generally applicable to all NMR tools that have a maximum density of state. While some tools have density of states that are nearly vertically, one skilled in the art would recognize that perfectly vertical density of states are only achievable in highly controlled laboratory environments. In reality nearly all NMR tools have some variation of a saddle point (or maximum density of state), whether sharply defined (nearly vertical) or smooth. Accordingly, the equations and methodology presented herein may be applied (with minor modifications) to nearly all NMR tools, whether designed for medical, earth logging, or food characterization purposes.
The logging device may include, for example, a permanent magnet or permanent magnet array 45, which may be made of a samarium-cobalt-magnetic material, and one or more rf antennas 50. The investigation region, or sensitivity zone, represented generally at 55, is a region in the formation in which the measurements are made. It will be understood that the present invention may be applicable to other tool configurations.
While the invention has been described herein with reference to certain examples and embodiments, it will be evident that various modifications and changes may be made to the embodiments described above without departing from the scope and spirit of the invention as set forth in the claims.