You must log in to edit PetroWiki. Help with editing

Content of PetroWiki is intended for personal use only and to supplement, not replace, engineering judgment. SPE disclaims any and all liability for your use of such content. More information


PEH:Nuclear Magnetic Resonance Applications in Petrophysics and Formation Evaluation

PetroWiki
Revision as of 15:57, 26 April 2017 by Denise Watts (Denisewatts) (talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Publication Information

Vol5REPCover.png

Petroleum Engineering Handbook

Larry W. Lake, Editor-in-Chief

Volume V – Reservoir Engineering and Petrophysics

Edward D. Holstein, Editor

Chapter 3E - Nuclear Magnetic Resonance Applications in Petrophysics and Formation Evaluation

Stephen Prensky, SPE, Consultant and Jack Howard, SPE, Halliburton

Pgs. 289-355

ISBN 978-1-55563-120-8
Get permission for reuse


Nuclear magnetic resonance (NMR) has been, and continues to be, widely used in chemistry, physics, and biomedicine and, more recently, in clinical diagnosis for imaging the internal structure of the human body. The same physical principles involved in clinical imaging also apply to imaging any fluid-saturated porous media, including reservoir rocks. The petroleum industry quickly adapted this technology to petrophysical laboratory research and subsequently developed downhole logging tools for in-situ reservoir evaluation (see the next section of this chapter).

NMR logging, a subcategory of electromagnetic logging, measures the induced magnet moment of hydrogen nuclei (protons) contained within the fluid-filled pore space of porous media (reservoir rocks). Unlike conventional logging measurements (e.g., acoustic, density, neutron, and resistivity), which respond to both the rock matrix and fluid properties and are strongly dependent on mineralogy, NMR-logging measurements respond to the presence of hydrogen protons. Because these protons primarily occur in pore fluids, NMR effectively responds to the volume, composition, viscosity, and distribution of these fluids (i.e., oil, gas, and water). NMR logs provide information about the quantities of fluids present, the properties of these fluids, and the sizes of the pores containing these fluids. From this information, it is possible to infer or estimate the volume (porosity) and distribution (permeability) of the rock pore space, rock composition, type and quantity of fluid hydrocarbons, and hydrocarbon producibility (Fig. 3E.1).

NMR logging provides measurements of a variety of critical rock and fluid properties in varying reservoir conditions (e.g., salinity, lithology, and texture), some of which are unavailable using conventional logging methods (Fig. 3E.1) and without requiring radioactive sources (Table 3E.1). Whether run independently as a standalone service or integrated with conventional log and core data for advanced formation and fluid analyses, NMR logging has significantly contributed to the accuracy of hydrocarbon-reservoir evaluation. During the past decade, a new generation of wireline-logging devices has been introduced into commercial service. In the past few years, logging-while-drilling (LWD) devices and downhole NMR spectrometers have also been introduced.




Historical Development


Within a few years after the first successful observations of NMR in 1946, and the demonstration of free-precession NMR in the earth’s magnetic field in 1948, the petroleum industry recognized the potential of NMR measurements for evaluating reservoir rocks, pore fluids, and fluid displacement (flow).

In the early 1950s, several companies—particularly California Research (Chevron), Magnolia (Mobil), Texaco, Schlumberger, and Shell—began extensive investigations to understand the NMR properties of fluids in porous media for the purpose of characterizing reservoir rocks (porosity, permeability, and fluid content).[1][2][3] In addition to laboratory research, these investigations included proposals for logging devices and the development of well-logging methods to permit formation evaluation in situ.[1][4] Although a number of patents for logging tools were issued in the 1950s, it was not until Chevron completed an experimental Earth’s field nuclear-magnetic-log (NML) logging device in 1958 that a functioning device was actually developed.[1][5] Limited commercial service of these devices was introduced in 1962 by Atlas, using the Chevron centralized design, and followed in 1965 by Schlumberger, using a pad-type tool of its own design. An improved version of the Schlumberger tool was introduced in 1978. Although the potential applications for this measurement were significant, particularly in the shallow, heavy-oil fields of the San Joaquin Valley,[6] in general, they did not live up to expectations and were not commercially successful.[7] Tool reliability and operational limitations proved to be major obstacles: the tool was not combinable, it required high (surface) power; the signal level varied geographically and was generally very low as a result of the low-operational frequency (2 kHz); and the borehole had to be doped with powdered magnetite to suppress the proton signal from the mud.[1][2][8] The final version of the Schlumberger NML tool—a centralized tool introduced in 1984—proved reliable and commercially successful and was in service until the advent of modern pulse-echo tools in 1994.

In 1978, the Los Alamos Natl. Laboratory developed a logging tool that employed permanent magnets and used a pulsed radio frequency (RF) (pulse-echo) NMR method. Although this particular design had serious limitations—such as a low signal-to-noise ratio (S/N) that prohibited continuous, nonstationary logging—the concept set the stage for the development of modern commercial NMR tools. This advance was soon followed by improvements in magnet and coil design that enabled continuous logging. During the 1980s, while developing a commercial logging tool, industry also carried out laboratory experiments to further understand NMR behavior in fluid-filled porous media and to develop petrophysical interpretations from these data. Ultimately, two wireline tools using different magnet and coil configurations emerged from these efforts: Numar’s mandrel device (MRIL) and Schlumberger’s skid (sometimes called "pad") design [Combinable Magnetic Resonance (CMR) tool]. Commercial logging began with these tools in 1991 and 1995, respectively. These wireline-tool designs continue to evolve (see section on Tool Design in this chapter). Recent improvements allow simultaneous acquisition of more measurements, operation in a wider range of borehole conditions, and faster logging speeds. Detailed accounts of the historical development of NMR and NMR logging are available in several published references.[4][9][10][11]

In addition to improvements in wireline tools, new acquisition schemes and processing methods have improved the resolution, quality, and utility of the acquired data and enabled enhanced interpretation methods and data analysis. Concurrent with wireline improvements, LWD NMR logging services were being developed and have been introduced in the past few years. In a related development, a downhole NMR spectrometer is now available for use with a formation-testing tool for in-situ fluid analysis.

NMR Physics


Atomic nuclei spin, and this angular moment produces a magnetic moment (i.e., a weak magnetic field). The NMR technique measures the magnetic signal emitted by spinning protons (hydrogen nuclei are the protons of interest in NMR logging) as they return to their original state following stimulation by an applied magnetic field and pulsed RF energy. These signals, which are observed (measured) as parallel or perpendicular to the direction of the applied magnetic field, are expressed as time constants that are related to the decay of magnetization of the total system.

NMR devices—both laboratory spectrometers and logging tools—use strong magnets to create a static magnetic field, B0, that aligns (polarizes) the protons in the pore fluid from their resting (random) state to the direction of the imposed magnetic field (Fig. 3E.2).

Polarization is not instantaneous—it grows with a time constant, which is called the longitudinal relaxation time, denoted as T1. Once full polarization (magnetic equilibrium) has been achieved, the applied static magnetic field, B0, is turned off.

The protons begin to lose energy as the imposed magnetization, M0, decays and the protons fall out of alignment, back to their original orientation and low-energy state. The protons’ angular momentum causes them to behave like tiny gyroscopes, and the loss of energy occurs during a wobbling or axial rotation (called precession) in the direction of the applied magnetic field. M0, also known as the bulk magnetization, provides the signals measured by NMR devices. The frequency at which the energy is emitted or is initially absorbed, f, called the Larmor or resonance frequency, is proportional to the strength of the external magnetic field, B0, (Fig. 3E.3). The Larmor frequency is used to tune a NMR probe, permitting it to image very thin slices of a sample at different distances from the tool.

An antenna detects and records the decaying magnetic field generated by the precessing nuclei. At any given time, t , the strength of this magnetic field, Mz, is proportional to the number of protons, the magnitude of B0, and the inverse of the absolute temperature (Eq. 3E.1):

RTENOTITLE....................(3E.1)

where Mz(t) = the magnitude of magnetization at t, M0 = the final and maximum magnetization at a given magnetic field, and t = the time that the protons are exposed to the B0 field.

The signal recorded parallel to the direction of the applied magnetic field (z plane) is called T1, or longitudinal (spin-lattice) relaxation. T1 describes how quickly the protons align within the static magnetic field. The T1 curve is an exponential curve that characterizes the rate of change of the proton magnetization (Fig. 3E.4).

T1 is the time at which the magnetization reaches 63% of its final value, and three times T1 is the time at which 95% polarization is achieved. Full polarization of typical reservoir-pore fluids may take several seconds. Large values of T1 (measured in milliseconds) correspond to weak coupling between the fluid and its surrounding environment and a slow approach to magnetic equilibrium, whereas, small T1 values represent strong coupling and a rapid approach to equilibrium.[1] Different fluids, such as water, oil, and gas, have very different T1 values. T1 is directly related to pore size and viscosity.

Pulse NMR devices use precisely timed bursts (pulse sequences) of RF energy that generate an oscillating magnetic field (B1) that tilts or "tips" the aligned protons perpendicular (x-y plane) to the direction of the applied magnetic field. The application of B1 results in a change in energy state that causes the protons to precess in phase to one another. These changes are known as NMR.

When the B1 field is turned off, the precessions of the protons are no longer in phase with one another, and the net magnetization decreases. In this situation, a receiver coil (antenna) that measures magnetization in the transverse direction will detect an exponential decaying signal called free-induction decay (FID); see Fig. 3E.5. NMR-logging tools use the same antenna to transmit the RF pulse (kilowatt scale) and receive the decay signal (nanovolt scale).

The FID signal measured in the x-y plane is called T2 —the transverse or spin-spin relaxation. In contrast to T1, T2 of hydrocarbons is much shorter (see Table 3E.2).

The primary objectives in NMR logging are measuring T1 signal amplitude (as a function of polarization), T2 signal amplitude and decay, and their distributions. The total signal amplitude is proportional to the total hydrogen content and is calibrated to give formation porosity independent of lithology effects. Both relaxation times can be interpreted for pore-size information and pore-fluid properties, especially viscosity.

In the laboratory, T1 is generally measured by either of two pulse sequences: inversion recovery or saturation recovery. Inversion recovery consists of a 180° spin inversion followed by a variable recovery time and then a 90° read pulse. The magnetization vector is entirely in the longitudinal range and, thus, has a higher dynamic range than the other method. Saturation recovery uses a 90° pulse, followed by a 90° read pulse. Saturation recovery is generally considered the more robust and efficient method. Although the actual T1 sampling sequence is very short—involving several short echoes trains, each of which requires only a few milliseconds—the total amount of time required to obtain the number of samples sufficient to define the T1 spectrum is significantly greater.

Depending on the activation used, the computation of a T1 spectrum requires at least 25% more, and sometimes double, the time needed for the computation of a T2 spectrum. In NMR logging, T1 measurement initially required either a stationary mode or very slow logging speeds. With the latest multifrequency tools, a technique used for speeding up T1 measurements is to make simultaneous measurements of the individual steps observed during a T1 recovery experiment in adjacent volumes; at least two such volumes are required. This technique enables T1 acquisition in less time, thereby permitting faster logging speeds.

T2 measurement uses the spin-echo technique,[12] in which the protons are first tipped into the transverse (x-y) plane by a 90° RF pulse and then inverted (flipped) by a subsequent 180° RF pulse at a fixed-time interval to rephase the dephasing protons. Rephasing the protons creates a detectable signal called a spin echo (Fig. 3E.6).

In practice, a sequence of pulses is used to generate a series of spin echoes (echo train) in which echo amplitude decreases exponentially with the time constant, T2. A variety of multiple-echo pulse sequences have been developed for different purposes.[11] In well logging and petrophysical studies, the most widely used is the Carr-Meiboom-Purcell-Gill (CMPG) sequence.[13][14] A polarization period is followed by a 90° tip pulse, which in turn is followed by a series of alternating RF pulses and measurements of echo amplitudes detected by the logging-tool antenna. Successive 180° pulses are applied at a fixed-time interval (echo spacing, TE), and the echoes are recorded between the pulses (Fig. 3E.7). By recording an echo train, T2 can be calculated from the decay in the height (amplitude) of successive echoes[11] using Eq. 3E.2:

RTENOTITLE....................(3E.2)

where Mx(t) = the amplitude of the transverse magnetization (i.e., the amplitude of the spin-echo train) at time t, and M0x = the magnitude of the transverse magnetization when t = 0 (i.e., the time at which the 90° pulse stops).

A single T2-pulse sequence may involve several hundred or thousand echoes. Only the amplitude (peak) of each spin echo is measured and stored. A series of echo trains is recorded and the signals stacked to improve S/N, especially at shorter relaxation times.

When recording multiple CMPG sequences, the time period between spin-echo recovery and the next 90° CMPG excitation—during which the protons are repolarized by the static magnetic field—is called the wait time, TW (Fig. 3E.8). Each CMPG sequence may use a different wait time, echo spacing, and number of echoes. An additional advantage of the CMPG sequence is that a small echo spacing, TE, in the CMPG sequence can minimize the diffusion effect on T2. CMPG measurement sets are always collected in phase-alternate pairs (PAP) to preserve the signal and to eliminate low-frequency electronic offsets. In general, pulse NMR offers better methods to measure relaxation times and quantify liquid displacement in rock.[15]

In-Gradient Diffusion

FID is caused by inhomogeneities in the magnetic field that are primarily caused by the existence of magnetic-field gradients. Gradients in the magnetic field occur, in part, because of the distance from the magnet to the sensitive (measurement) volume. For a given geometry, the gradient is inversely related to magnetic-field strength. Compared to laboratory and clinical NMR devices, NMR-logging tools produce a relatively weak and inhomogeneous static magnetic field. In the case of reservoir rocks, differences between the magnetic properties of the rock matrix and pore fluids may also contribute to a magnetic-field gradient. T2, but not T1, is affected by this phenomenon, which is called diffusion. In the presence of high magnetic-field gradients, diffusion effects make T2 interpretation difficult. However, because the gradients produced by NMR-logging tools are relatively constant, they can be accounted for in T2 interpretation. In fact, the existence of these field gradients has actually proved beneficial in NMR logging. Magnetic resonance imaging (MRI) is the process by which NMR measurements are obtained in a gradient magnetic field.

NMR Petrophysics

Laboratory Studies

Extensive laboratory studies on NMR behavior and on the properties of fluid-saturated porous media have been conducted since the inception of NMR and throughout the development of NMR-logging tools. The results from these investigations have provided the petrophysical foundation for understanding the logging measurements and for developing interpretation models and applications.

Low-field, bench-top pulse-NMR spectrometers were developed concurrently with logging tools so that wellbore measurements could be duplicated on core samples in the laboratory.[16][17] These instruments operate and record data in the same manner as NMR-logging tools.[18] Because NMR analysis is nondestructive, NMR and conventional capillary-pressure measurements can be performed on the same samples, in both the saturated and partially saturated states. Low-field spectrometers provide the ability to make repeatable measurements of rock- and fluid-NMR properties. This ability, in turn, permitted correlation and calibration of laboratory and field measurements and also permitted direct transfer of interpretation models developed in the laboratory to logging data. Where core is unavailable for NMR-log calibration, new technology and methods now allow NMR petrophysical measurements on drill cuttings.[19]

Laboratory NMR studies are routinely conducted for the following purposes:

  • Verifying formation porosity.
  • Evaluating textural effects, such as microporosity, on NMR-log response.
  • Determining formation-specific models that enhance the accuracy of determining bulk-volume-irreducible (BVI) water, free-fluid index (FFI), and, ultimately, permeability.
  • Developing models to identify and quantify hydrocarbons, including residual oil.
  • Developing models to predict changes in pore size (facies).


Much of this work is summarized in Kenyon,[20] Murphy,[21] Woessner,[3] and Dunn et al.[11] The most recent laboratory studies have suggested that some established NMR core-log relationships should be further investigated to better account for data scatter.[22] A related area of study not dealt with here is NMR imaging of fluid flow in core.[23]

As in NMR logging, data quality is critical. To achieve the desired objectives, laboratory NMR studies should include a preplanning phase similar to that used in logging. (See the Job Planning section of this chapter).

Petrophysical Properties

The basic petrophysical parameters of porosity, permeability, and producibility can be determined from either T1 or T2 echo-decay data. Until low-field spectrometers were developed, T1 was the preferred acquisition method in the laboratory, where time is not a concern.[18][20] Because T1 measurement requires more time than T2, T2 became the primary acquisition mode in pulse-NMR logging because it allowed logging at speeds that were commercially viable. Fortunately, there is a correlation between T1 and T2[24]; and T1 can be estimated from T2 data by extrapolating the T2 decay-obtained-polarization pulses of different lengths (i.e., using different values of TW; see Fig. 3E.9).

Modern logging tools are capable of operating in either T1- or T2-acquisition modes. The logging mode is dictated by operational factors and job objectives and may, in fact, switch back and forth, as needed. (See further discussion in the LWD Tool section of this chapter.)

NMR Properties of Fluids

T1 relaxation occurs when the precessing proton system transfers energy to its surroundings. T2 relaxation occurs through a similar transfer in energy and also through dephasing. Consequently, transverse relaxation, T2, is always faster than longitudinal relaxation, T1. The emphasis of the proton NMR techniques used in formation evaluation is on the NMR fluid response from the pore fluids where T2T1. The NMR response in solids (i.e., the non-shale/clay component of the rock matrix) is very short compared to the pore-fluid signal, and is generally not measured by laboratory or logging devices.

NMR relaxation of fluids depends on whether the fluid is measured in bulk form as a wetting-pore fluid within a rock matrix or in a gradient magnetic field. Bulk relaxation is the intrinsic relaxation property of a fluid that is controlled by viscosity, chemical composition, temperature, and pressure: T2bulkT1bulk. Fluids contained within pores have different relaxation characteristics, namely those of surface relaxation.

Surface relaxation occurs at the fluid-solid interface between a wetting-pore fluid and the rock-pore walls (Fig. 3E.10) and is different from the relaxation in either the solid or the fluid, individually. Surface relaxation dramatically decreases both T2 and T1 and is the dominant component contributing to T2. When a nonwetting fluid (e.g., oil) is also present in rock pores, the nonwetting fluid may continue to relax at its bulk-relaxation rate.

Surface relaxation is expressed by the following equations (Eqs. 3E.3 and 3E.4):

RTENOTITLE....................(3E.3)

and

RTENOTITLE....................(3E.4)

where ρ2 = T2 surface relaxivity (i.e., T2 relaxing strength of the grain surfaces); ρ1 = T1 surface relaxivity (i.e., T1 relaxing strength of the grain surfaces); and (S/V) pore=ratio of pore surface to fluid volume.

Surface relaxivity varies with mineralogy; and, for simple pore shapes (S is the surface area of a pore and V is the volume of the same pore), S/V is a measure of pore size. In a brine-wet rock, T2 in smaller pores will be less than T2 in large pores; consequently, identical pore water in different rocks can have a wide range of relaxation times because of variations in surface relaxivity. Laboratory studies have demonstrated that in water-wet rocks, the surface and volume ratio (S/V) is also a measure of permeability.

Fluids controlled by surface relaxation exhibit T2 values that are not dependent on temperature and pressure. For this reason, laboratory NMR measurements made at room conditions are commonly used to calibrate formulas used to estimate petrophysical parameters such as permeability and bound water.[25][26]

Diffusion-induced relaxation occurs when a significant gradient exists in the static magnetic field. Molecular diffusion in this gradient causes additional dephasing that contributes to increased T2 relaxation. In addition to the magnetic-field gradient, diffusion is also controlled by inter-echo spacing and fluid diffusivity, viscosity, molecular composition, temperature, and pressure.

Bulk-fluid processes and surface relaxation affect both T1 and T2, while diffusion only affects T2 relaxation. All three processes are independent and act in parallel according to the following equations (Eqs. 3E.5 and 3E.6):

RTENOTITLE....................(3E.5)

and

RTENOTITLE....................(3E.6)

The relative importance of the three diffusion-relaxation mechanisms depends on the fluid type in the pores (e.g., water, oil, or gas), the sizes of the pores, the strength of the surface relaxation, and the wettability of the rock surface (see Table 3E.3).

T2 Decay

Eq. 3E.2 states that the T2 decay associated with a single pore size in water-saturated rocks is proportional to the pore size.[20] In fact, because reservoir rocks typically comprise a distribution of pore sizes and frequently contain more than one fluid type, a CMPG T2 spin-echo train actually consists of a distribution of T2 decays, rather than a single T2 decay. In these cases, the exponential decay is described by Kenyon et al.[27] as follows:

RTENOTITLE....................(3E.7)

where M(t) = measured magnetization at t; Mi(0) = initial magnetization from the ith component of relaxation; and T2i = decay constant of the ith component of transverse relaxation. The summation is over the entire sample (i.e., all pores and all different types of fluid).

Fig. 3E.11 illustrates the multiexponential decay character of a porous medium containing pores of different sizes and a single wetting phase. Surface relaxation dominates when a short inter-echo spacing is used and the formation is only brine saturated. Under this condition, T2 is directly proportional to pore size. When all pores are assumed to have similar geometric shape, the largest pores (see Fig. 3E.11, left column) have the lowest S/V and, thus, the longest T2. Medium-size pores have smaller S/V, yielding shorter T2 values. The smallest pores have the highest S/V and the shortest T2 values.

Eq. 3E.7 can also be expressed as follows: [28]

RTENOTITLE....................(3E.8)

where (S/V)t is the surface-to-volume ratio of the ith pore. When t = 0, the following is true:

RTENOTITLE....................(3E.9)

If the measured magnetization for 100% bulk water, M100%(0), is known, then M(0) and M0i can be calibrated to porosity by the following equation:

RTENOTITLE....................(3E.10)

where ϕ = calibrated porosity of the formation, ϕi = calibrated porosity associated with all pores of the ith pore size (also known as the incremental porosity). Therefore, the T2 distribution (in the form of the amplitudes, M0i, associated with the time constants, T2i, is calibrated to the porosity distribution (i.e., the individual pores ϕi with the associated time constants T2i).

If the rock is water-wet and the pores are partially saturated (i.e., the pores contain oil and/or gas in addition to water), then the total signal comprises contributions from each component in the following equation:

RTENOTITLE....................(3E.11)

where Moil = magnetization produced by oil protons in the pores, Mgas = magnetization produced by gas protons in the pores, T2oil = T2 of oil measured with a CMPG sequence, and T2gas = T2 of gas measured with a CMPG sequence.

In a brine-saturated rock, the T2 decay spectrum represents a pore-size distribution. However, when nonwetting fluids (e.g., oil or gas) are present, the T2 spectrum includes a bulk response from the nonwetting fluid, in addition to the pore-size response. Pores containing the nonwetting fluid either appear in the spectrum at a decay time that is faster than is normally associated with the pores, or do not appear at all if the surface layer is too thin. This behavior affects the appearance of the T2 spectrum and associated T2 distribution but not the total signal amplitude (i.e., porosity).

Data Fit-Inversion

The raw data recorded by an NMR device are a series of spin-echo amplitudes (echo train) as a function of time, usually at fixed time increments (bins). These NMR measurements are statistical, and stacking is required to improve the precision of the log outputs. The precision depends on formation-fluid properties, activation (acquisition) type, and the number of pulses stacked.[29] The data are mathematically inverted (mapped) by use of a best-fit curve to produce a distribution of T2 values as a function of relaxation time (Figs. 3E.12 and 3E.13). Initially, biexponential-fitting algorithms were used[30]; however, T2 decay in fluid-saturated rocks is multiexponential, because of this, multiexponential models were developed and are commonly used for inverting the data.[11] The inverse solution (T2 distribution) is a function of both the measured echo data and the chosen smoothness (regularization) for the inversion. However, because regularization is controlled in part by the S/N, the fit that is actually used is not unique; that is, there can be a number of differently shaped T2 distributions that fit the original echo-decay curve. In general, the area under the T2-distribution curve (porosity) and the general location in time of the high-porosity bins are robust; however, caution is advised when interpreting the fine details of the distribution.[31] Because the spin echoes are measured over a very short time, an NMR tool travels no more than a few inches along the wellbore while recording the spin-echo train; thus, the recorded spin-echo data can be displayed on a log as a function of depth.


T2 Distribution

The mathematical statement (Eq. 3E.10 and Eq. 3E.11) that T2 distribution observed in water-saturated rock represents the pore-size distribution and porosity of the rock has been confirmed using mercury-injection capillary pressure (MICP) methods (see Figs. 3E.14 and 3E.15).[32][33]

NMR responds to pore-body size, and MICP responds to pore-throat size.[32][34][35][36] In clastic rocks in which there is a good correlation between pore-body and pore-throat size, there is often good qualitative agreement between NMR and MICP data, as Figs. 3E.14 and 3E.15 illustrate. The NMR porosity, reflected in the T2 distribution, is a spectrum comprising rock matrix and fluid components (see Fig. 3E.16).

Although matrix minerals and dry clay may contain hydrogen atoms in the form of hydroxyl groups, the T1 relaxation times of these nuclei are too long to be polarized by a moving NMR-logging tool, and their T2 relaxation times are too short to be recorded.[37] The fluid component is subdivided into bound and free subcomponents. The hydrogen nuclei of clay-bound water are adsorbed on the surfaces of clay grains. These hydrogen protons can be polarized by NMR-logging tools and recorded when a sufficiently short TE is used.[38] Similarly, hydrogen protons in capillary-bound water and movable fluids (e.g., free water, mud filtrates, oil, and gas), are polarized and recorded by NMR-logging tools with appropriate values for TE and TW.

The porosity and pore-size information from NMR measurements can be used to estimate producible porosity (i.e., the movable fluids) and permeability and for hydrocarbon identification (see the Applications section of this chapter).

Logging Tools


Compared to laboratory or clinical NMR, the dimensions of a typical borehole and the nature of continuous logging impose severe constraints on the physics, equipment, and operation of NMR-logging tools. Unlike laboratory devices, logging measures an external sample by use of a weaker magnetic field while in motion relative to the sample.

NML Tool

The NML tool—the first generation of NMR logging (1960–1994)—was an Earth’s field device that measured the free-induction decay in the Earth’s magnetic field. Proton polarization (alignment) was achieved using a magnetic field produced by a coil energized with a strong direct current. Each experiment required several seconds to allow complete polarization. The power was turned off, and the same coil was used to receive the free-induction signal.

There were a number of problems inherent in this technology. NML devices measured the proton signal in the borehole fluid as well as in the surrounding formation. Magnetite powder was circulated into the mud system to cancel the borehole hydrogen signal thereby preventing it from overwhelming the formation signal. This was a time-consuming process. The low intensity of the Earth’s magnetic field (0.5 gauss) resulted in a low S/N. The NML had a large dead time (i.e., time between the cessation of the static magnetic field and the first measurement), approximately 20 ms; pore fluids with relaxation times less than the dead time were not measurable. The initial signal amplitude had to be extrapolated backwards from the subsequent free-induction decay data. The NML could not distinguish between oil and water or measure total porosity; however, the dead time was close to the T2 cutoff established for irreducible water saturation. The NML could detect movable fluids and provided a measurement called the Free Fluid Index (FFI).[11] An expanded discussion of the NML technology can be found in Brown and Neuman[39] and Brown.[1]

Pulse NMR

Unlike conventional magnetic resonance imaging (MRI) devices in which the sample is placed inside the stationary coil, a borehole-logging device investigates a sample (rock volume) that is outside the device itself while moving along the borehole. This setup has been called the "inside-out" NMR problem because the geometry of the magnets and coils is inverse to that used in laboratory-NMR spectrometers.[40][41][42][43][44] To obtain useful measurements, a logging tool must generate a large, radially symmetric static magnetic field and also a high-frequency oscillating magnetic field. Each must be capable of penetrating one or more inches into the formation surrounding the borehole. The diameter of a typical well limits the size of the permanent magnets that can be used and, therefore, the strength of the magnetic field that can be generated by a logging device. In comparison with laboratory-NMR and clinical-NMR devices, which may operate at 10 MHz and the static magnetic fields of which are commonly in the range 1 to 2 Tesla (high field), modern logging tools and laboratory spectrometers designed for petroleum investigations are considered low-field devices. They operate at or below 2 MHz and generate relatively weak (< 200 gauss) and inhomogeneous magnetic fields (i.e., gradients up to 20 gauss/cm). By comparison, the NML operating frequency was only 2 kHz, and the Earth’s magnetic field is only 0.5 gauss. These factors limit borehole investigation to protons (hydrogen) and the use of relaxation data only. Chemical shifts (widely used in the biological field) are not observable.[7] Furthermore, to compensate for the lower S/N that results from low-field intensity, logging tools must acquire more echoes and/or stack data to improve S/N. As mentioned previously, with the introduction of pulse-echo tools, T2 became the primary acquisition mode because it permitted faster logging—a major factor in high-cost wells.

This new design provides a number of operational advantages over the NML:

  • Using permanent magnets rather than electromagnets reduces the surface-power requirement.
  • Focusing the sensitive region of the magnetic field at some distance into the formation eliminates the requirement for suppressing the mud signal.
  • Using an RF pulse from a coil tuned to the Larmor frequency ensures that only nuclei in the sensitive region are in resonance.
  • Controlling the pulse duration means shorter dead times that, in turn, allow a better estimate of initial decay amplitude (porosity) measurement for short T2 components (bound-fluid evaluation).
  • Enabling more sophisticated pulse sequences allows for measurement of additional rock and fluid properties.[44]


One initially unwanted product of this design is the creation of gradients in the static magnetic field that causes molecular diffusion. The strength of the magnetic field gradient, G, is controlled by tool design and configuration (i.e., tool size and tool frequency); by environmental conditions such as formation temperature; and by internal gradients induced by the applied field, B0. Subsequent characterization of these gradients has enabled the in-gradient diffusion to be used for hydrocarbon typing.

Wireline-Tool Designs

From the beginning, a debate in NMR logging has been whether to use a centralized or eccentered tool design. There are two different wireline designs in current commercial service: (1) the Numar (now a part of Halliburton) magnetic resonance imaging tool (MRIL), a centralized mandrel device[30][45][46][47][48]; and (2) the Schlumberger CMR tool[49][50][51][52][53][54][55] and the Baker Atlas MREX[56] eccentered devices, both of which require contact with the borehole wall. The general designs of both the MRIL (mandrel-type) and CMR (pad-type) have evolved over the past decade with the addition of new capabilities and faster logging speeds, made possible by improved electronics and improved data acquisition and processing. Computalog is field testing a prototype mandrel device: the Nuclear Magnetic Resonance Tool (NMRT) developed by NPF Karotazh.[57]

A centralized NMR-logging tool like the MRIL must be longer than a pad device, simply to contain magnets of sufficient strength to project the required magnetic field, through the borehole and into the formation. This factor results in a greater sensitivity to borehole salinity than a pad device, which can exclude the mud from its measurement. The MRIL design generates relatively thin (1 to 2 mm) sensitive volumes, but the reduced S/N that accompanies these volumes is compensated by the vertical size of the sensitive area. It also generates a relatively high magnetic gradient. In contrast, a contact device, such as the CMR, can use smaller magnets and electronics, which provide higher vertical resolution but a shallower depth of investigation (DOI) and greater sensitivity to borehole rugosity. In addition to the standard permanent magnets, some designs now include "prepolarization" magnets, which are added to ensure full polarization at typical logging speeds.

The latest wireline NMR-logging tools operate simultaneously in several RF frequencies to measure (image) multiple sample volumes. In the presence of a gradient magnetic field, pulses with different frequencies will cause protons in different (and parallel) regions of space (i.e., measurement or sensitive volumes) to resonate.[31] Cycling through several frequencies excites protons in different cylindrical volumes, allowing measurements to be made more quickly. If the frequencies of multifrequency measurements are very close together, then the sensitive volumes are very close together; and, for practical purposes, the rocks sampled can be considered to be the same. This principle is the same as that used for slice selection in medical MRI imaging. These tools acquire multiple echo trains using different values of TW, TE, and variable magnetic gradients (G) in a single logging pass. The time between measurements made at multiple frequencies can be as little as the time of an echo train, and the time between measurements made at a single frequency is essentially the time needed to repolarize (TW). The thickness of these sensitive volumes may be as small as 1 mm. Furthermore, recent advances in tool design now permit cost-effective T1 acquisition.[58] Table 3E.4 summarizes the capabilities, advantages, and disadvantages of T1 and T2 acquisition. These differences influenced the designs of the LWD tools discussed in the next section. Multifrequency operation provides measurements at multiple DOIs (typically 1 to 4 in.). This allows invasion effects to be accounted for in the data interpretation, thus enabling determination of near-wellbore fluid saturation and oil properties at high resolution.[54][58][59][60]


LWD Tools

The latest entries into NMR logging are LWD tools. The development of LWD-NMR is ongoing and significant changes in hardware design, as well as significant changes and improvements in data acquisition and processing, can be expected in the next few years. The general benefits of LWD have been discussed elsewhere—in particular, NMR-LWD offers a nonradioactive alternative for porosity measurement, an NMR alternative to wireline in high-risk and high-cost wells, and enables high-resolution fluid analysis in thin beds and laminated reservoirs.[61]

By definition, logging tools operating in the drilling environment are built into drill collars and are, therefore, mandrel devices. In contrast to wireline tools, they must be capable of making measurements when the drillstring is stationary, sliding or rotating, centered or eccentered.

LWD-tool measurements are either omnidirectional or azimuthal, depending on tool design. The incorporation of magnetometers allows binning of the data into azimuthal sectors. Omnidirectional measurements can be generated from azimuthal data, but not the reverse. Although current LWD-NMR services provide only omnidirectional data, patents have been issued for tools with azimuthal capability.[62]

A major concern introduced with the advent of LWD-NMR is the effect of drillstring lateral motion on the basic NMR measurement.[63] NMR measurements are not instantaneous, they involve both polarization and decay, time-dependent components. Lateral movement of a wireline tool or LWD drillstring shifts the polarization volume and the sensitive-measurement area, relative to one another, and these shifts may result in incomplete polarization or incomplete measurement of the decay. Low-velocity motion affects only the decay, but high-speed motion can also affect the initial decay amplitude.

The currently available LWD tools offer different solutions to this concern. While both tool designs can operate in either T1 or T2 acquisition mode and both incorporate accelerometers and magnetometers for detecting lateral motion for quality control of the NMR measurements, they differ in their choice of primary measurement mode. Operational factors, such as the slower logging speed (i.e., rate of penetration in drilling), compared with typical wireline logging, also affect the choice of measurement mode.

The NMR acquisition sequences are programmable and interchangeable with those used in wireline tools. Switching between these acquisition modes is accomplished by a variety of methods, including elapsed time, counting measurements, and differentiating between drilling and nondrilling conditions. T1-mode records echo amplitudes as a function of time. Data output consists of porosity, free-fluid, and bound-water volumes, and can provide a quick-look permeability estimate. T2-mode is a multifrequency mode that records multiple wait-time CMPG spin-echo amplitudes and is capable of using all wireline-pulse sequences. T2-mode output includes porosity, free fluid, clay- and capillary-bound water, and differences in the multiple wait-time data are used for hydrocarbon indication.

The Halliburton tool (MRI-LWD)[64][65][66][67][68] uses T1 as its preferred acquisition mode. Halliburton considers T1 more robust for determining porosity and free-fluid volume. The anticipated maximum rates of penetration for LWD—1 to 3 ft/min—allow T1 acquisition. T1 is motion tolerant compared with T2. A sequence of interleaved measurements made at different recovery times is used to construct the T1 relaxation decay (buildup). As long as the sensitive volume (shell) is contained within the much larger volume reached by the saturation pulse, the measurement is valid. During post-processing, drilling and nondrilling periods are identified, and the invalid T1 data recorded during drilling are discarded.

Schlumberger’s tool (ProVision LWD-NMR) [69][70][71] uses T2 for its primary measurement. While both companies agree that T1 is motion tolerant, Schlumberger considers T1 to be less robust for estimating porosity, bound-fluid, and free-fluid volumes because of the poor S/N resulting from the longer time require for equivalent-quality T1 measurements. A reduced S/N impacts data quality (e.g., statistical repeatability and vertical resolution), logging speed, and, ultimately, the results. T2 measurement also enables rapid calibration and correlation with the large body of wireline-NMR data. The multiple wait-time acquisition includes a fully polarizing (3 to 12 seconds), a partially polarizing (normally 0.6 to 1 seconds), and a very fast wait time (typically 0.08 seconds). Porosity and bound-fluid volume are calculated for both the long- and the medium-wait-time measurements, and significant differences between the long-wait-time and medium-wait-time porosity provide real-time hydrocarbon indication. When necessary, logs are corrected for incomplete polarization of the hydrogen nuclei, and T1 distributions are estimated from the measured T2 relaxation.

Baker Atlas is presently field testing a new LWD device (MagTrak) that operates in T2 acquisition mode and has a vertical resolution of less than 3 in.[72] This tool design achieves a motion-tolerant T2 measurement through a combination of a very-low gradient magnetic field, circuitry that permits a TE as low as 0.6 ms, and the use of special stabilizers.

Downhole NMR Spectrometer

Contamination of hydrocarbon reservoirs by oil-based mud (OBM) and synthetic oil-base mud (SOBM) is a significant problem for accurate reservoir and fluid analysis in wells in which OBMs are used, especially offshore. Direct knowledge or an estimate of the OBM’s NMR characteristics is required to distinguish it from connate oil in fluid samples obtained by formation testers. Acquisition of T1, T2, and D0 (the self-diffusion coefficient) permits evaluation of connate oil, oil viscosity, and gas/oil ratio.[31][73] Furthermore, in costly drilling environments, real-time acquisition of NMR properties permits immediate, rather than delayed, fluid evaluation. Information on NMR fluid properties is also valuable for the interpretation of NMR wireline and LWD logs.

The recently introduced downhole NMR spectrometer, incorporated into a formation-testing tool, can obtain NMR measurements of OBM contamination directly, on live samples at in-situ conditions.[74][75][76][77] As the testing tool pumps fluid from the reservoir into the borehole or sample chamber, the spectrometer—using a measurement time of 30 seconds—measures the hydrogen index (HI), T1, T2, and diffusion. The T1 measurement is made while flowing; T2 and diffusion are static measurements. The T1 distribution is important in differentiating between highly refined OBM filtrates and native oil. The T1 characteristics of these common filtrates are measured and cataloged so that the data from oil-based filtrates can be distinguished from native hydrocarbons.

Log-Presentation


NMR-log data are presented in a variety of formats designed to emphasize specific aspects of the data and thus enable rapid visual interpretation of movable and immovable fluids, porosity, and permeability. Data interpretation is further enhanced when additional log and core information are also included in the log presentation.

The T2 distribution is typically displayed in waveform presentation, image (variable-density log, VDL) format, and bin-distribution plot. Each T2 format represents the distribution of the porosity over T2 values, and hence, over the pore sizes (Figs. 3E.17 and 3E.18).


NMR Applications


NMR-log data can be analyzed independently or in combination with conventional-log and core data. As an independent logging service, NMR can provide porosity, permeability index, and complete information on fluid type and saturation of the flushed zone. Some data-interpretation methods operate in the echo-decay time domain, while others operate in the T2-relaxation domain.

Porosity Determination with NMR

The initial amplitude of the spin-echo train is proportional to the number of hydrogen nuclei associated with the fluids in the pores within the sensitive volume. This amplitude is calibrated in porosity units (see Eq. 3E.10). The accuracy of this amplitude measurement depends on three factors: first, a sufficiently long polarization time, TW, is needed to achieve complete polarization of the hydrogen nuclei in the fluids; second, a sufficiently short inter-echo spacing, TE, is needed to record the decays for fluids associated with clay pores and other pores of similar size—if either TW is too short or TE too long, NMR porosity will be underestimated; third, the number of hydrogen nuclei in the fluid should equal the number in an equivalent volume of water (i.e., HI = 1); if the HI of any of the pore fluids is significantly less than 1, a correction to porosity is necessary.

Porosity was one of the earliest NMR measurements and is still an important one. Assuming that the logging tool is properly calibrated and functioning normally (see discussion of Quality Control in this chapter), interpretation of the porosity measurement depends on several other factors, including vintage of the logging tool, whether hydrocarbons are present, and fluid type.

There are two major contrasts between NMR porosity and conventional density and neutron-nuclear-porosity logs. First, NMR porosity does not depend on the mineralogy of the matrix, except in cases in which the formation contains significant amounts of ferromagnetic or paramagnetic materials, and in most cases, it is considered a lithology-independent measurement. Second, NMR porosity is not sensitive to either borehole or mudcake and does not require corrections because the measurement zone (i.e., sensitive volume) is focused within the formation, beyond the borehole wall and these influences. The accuracy and precision of NMR-derived porosity has been confirmed through comparisons with core porosity obtained using conventional laboratory measurement methods (Fig. 3E.19).

Influence of Tool Version on the Porosity Measurement. The original NML tool had such a large dead time (e.g., minimum echo spacing, TE) that it could not measure clay- and capillary-bound fluids; thus, it was only capable of measuring free-fluid porosity, which is now called NMR effective porosity (MPHI). The minimum TE used in the early versions of pulse-NMR tools was still relatively large, and those tools still could not see all the clay-bound fluids. With improvements in logging technology, most NMR logs recorded after 1997 are considered "total-porosity" logs. Current logging tools use a minimum TE of 0.3 ms in continuous mode and a TE as small as 0.1 ms in stationary mode.

In clean, water-filled formations, MPHI should approximately equal neutron-density crossplot porosity. In shaly sands, MPHI should approximately equal density porosity, calculated with the correct grain density; however, the MPHI may not equal effective porosity because of the effects of HI and long T1 components (Eq. 3E.12):

RTENOTITLE....................(3E.12)

where MPHI is measured by the NMR tool; ϕe is effective porosity of the formation; HI is related to the amount of fluid in the effective porosity system; TW is polarization time used during logging; and T1 is longitudinal relaxation time of the fluid in the effective porosity system. Older NMR logs may indicate very low porosity in shales.

MPHI is almost always less than NMR total porosity (MSIG):

RTENOTITLE....................(3E.13)

where MSIG is measured by NMR total-porosity logging, and clay-bound-water porosity (MCBW) is measured by the NMR tool with partial-polarization acquisition. In very clean formations, however, NMR MCBW is virtually zero, and then MPHI equals MSIG (see Fig. 3E.16).

NMR total porosity may approximate density or neutron-density porosity. The reported clay-bound volume should not be considered absolute because the boundary between clay- and capillary-bound porosity components can vary with clay type and distribution. Whenever possible, the relationship of clay type and volume should be confirmed through core analysis and comparison with other log data.

Porosity-Logging Modes. Several modes of porosity logging are currently in use: standard T2 mode, short TW (e.g., bound-fluid logging), intermediate TW (e.g., polarization correction and the "clean, wet formation" method), and long TW, or multiple-wait-time acquisition (e.g., total-porosity logging).

Total-porosity-logging mode provides data to determine porosity, permeability, and productivity (mobile fluids). Total-porosity acquisition is a multiple-wait-time acquisition that acquires two echo trains to obtain the total porosity. One echo-train acquisition uses a long TE (e.g., 0.9 or 1.2 ms) and a long TW (3 to 8 seconds) to achieve complete polarization. This echo train provides the effective porosity. The second echo-train acquisition uses a short TE (0.2 to 0.6 ms) and a short TW (20 ms) that is only long enough to achieve complete polarization of the fluids in the small pores (i.e., clay-bound porosity).[46] Although the long TW acquisition may achieve full polarization in gas zones, the porosity may still require HI correction. The short TW echo trains are stacked to improve S/N. Total porosity is obtained by combining the two separate T2 distributions, usually at 4 ms. The standard T2 acquisition is used in situations in which there is little to no diffusion or T1 contrast. The benefits of standard T2 acquisition are increased logging speed without degradation of data quality; or superior data quality at normal logging speeds.

In bound-fluid logging mode, TW is kept relatively short, ranging from 0.3 to 1.0 seconds. The benefit of using a short TW is faster logging speeds while determining the clay- and capillary-bound fluid volumes. The drawback is underpolarization of the free-fluid volume.

When this logging mode is used, the free-fluid volume is computed from other log data (e.g., density, neutron-density crossplot, or acoustic). In polarization mode (CMR tool), a correction to the free-fluid volume is made using the observed relationship of T1 to T2, on the basis of the T2 acquired while logging. The magnitude of the polarization correction must be monitored because, in addition to porosity, the T1/T2 ratio is also influenced by changes in the types of formation and pore fluid. This method generally works well in clastics, but in conditions of high T1 and T2, such as in carbonates or light oil, the correction may actually introduce additional error.[78]

Another mode of operation is to find a clean, wet formation and perform a sweep of different wait times (TW) ranging from 2 to 5 seconds to determine the one needed for full polarization. Because oil and gas may have very long T1 values or require an HI correction, the use of intermediate TW methods does not guarantee that the NMR porosity will be fully reported.

Optimization of the acquisition is critical to accurate NMR measurements and interpretation (see the Job Planning section of this chapter). The effect of incorrect acquisition parameters on the T2 distribution is illustrated in Fig. 3E.20. In this figure, the sequence of graphs shows the effects of TE, TW, pulsing time (TP), and S/N on the measured NMR characteristics of a core sample:

  • A: Initial acquisition results show an apparent unimodal distribution.
  • B: Increasing the TP shows that this sample is, in fact, bimodal. There is a significant change in the T2 distribution, but little to no change in the cumulative porosity.
  • C: Decreasing the echo spacing (TE) reveals fast components in this sample that were previously masked by a TE that is too long. There is a significant change in the cumulative porosity with a characteristic shift of the fastest components in the distribution to shorter T2 values.
  • D: Increasing the wait time (TW) has significantly increased cumulative porosity. After correction for a TW that is too short, there is a characteristic shift in the slow T2 components to longer T2 values.
  • E: Increasing the S/N results in sharper resolution of the two major components. Although the definition of each component in the T2 distribution improves, little change is noted in the cumulative porosity.
  • F: The final result, in which all acquisition parameters have been optimized (red), is compared with the initial result (blue).


BVI Determination. Determining the BVI of water in a formation is one of the earliest and most widely used applications of NMR logging. BVI refers to the immovable or bound water in a formation. BVI is a function of both the capillary-pressure curve for the rock and the height above free water (Fig. 3E.21).

In practice, BVI serves as an indicator of, and is used for, estimating producibility and permeability. There are two methods currently used for BVI determination. One is the cutoff-BVI (CBVI) model, which is based on a fixed-T2 value (T2cutoff) that divides BVI and FFI. The second is the spectral BVI (SBVI) model, which assumes the coexistence of free and bound fluids in any pore described by the T2 distribution at water saturation (Sw) = 100%. This coexistence is expressed through a weighting factor that defines the fraction of bound water associated with each pore size.[79] The SBVI method is used primarily for quantifying movable water and, secondarily, for estimating permeability.

CBVI (Cutoff) Model. The T2 signal from the rock matrix is actually so rapid that even modern logging tools cannot detect it. Consequently, the recorded T2 distribution represents only the porosity occupied by movable (FFI) water and immovable BVI and clay-bound-water (CBW) components. Assuming that the mobility of reservoir fluids is primarily controlled by pore size (i.e., the producible fluids reside in large pores and the immobile, or bound, fluids reside in small pores), a fixed-T2 value can relate directly to a pore size at or below which fluids will not move. This value (T2cutoff) is used to divide the T2 distribution into movable (i.e., producible or free fluids and FFI) and immovable (i.e., bound-fluid, BVI, and CBVI) components (Fig. 3E.22).[80] T2cutoff is a variable that differs from one formation to another and is influenced by a variety of factors including capillary pressure, lithology, grain size, compaction, and pore characteristics.

In practical use, T2cutoff is either determined in the laboratory or a default value is assumed. If time and expense permit, accurate T2cutoff values can be defined by comparing the T2 distributions obtained on fully and partially water-saturated core samples taken from the logged interval[81][82] (Fig. 3E.23).

In the absence of laboratory data, lithology-dependent default values are used for T2cutoff. Default T2cutoff values are based on local field experience or common practice; for example, typical values used in the Gulf of Mexico are 33 ms for clastics (sandstones) and 92 or 100 ms for carbonates.[83][84][85] When using a default value, remember that, in fact, the T2cutoff can actually vary within a single lithology because the capillary pressure at which irreducible water saturation (Swirr) is achieved varies between the actual rocks. In addition, T2cutoff is also influenced by pore-wall chemistry, the presence of minor paramagnetic or ferromagnetic components, texture, and pore-throat to pore-body ratios. The use of an incorrect value for T2cutoff may result in underestimation of recoverable reserves or in bypassing a pay zone.

SBVI (Spectral) Model. The SBVI model was developed to address the limitations of the CBVI model in some rocks, particularly in very-well-sorted rocks in which there is a very narrow range of grain and pore sizes. In these rocks, the NMR echo-decay can typically be fit by a single-exponential decay that produces a sharp peak in the T2 distribution. When a fixed cutoff is used for determining BVI, it may result in significant error because of the presence of immovable water contained in microporosity associated with pore irregularities.[84] This microporosity component is apparent at Swirr, but not at Sw = 100% (Fig. 3E.24).

A weighting function defines the fraction of bound water in each pore size that is present in the T2 distribution at Sw = 100% (Fig. 3E.25). Several methods have been proposed for obtaining the weighting function.[31][85][86]

In general, the SBVI method is superior to the CBVI method for determining BVI in rocks in which Sw = 100%. The presence of hydrocarbons (i.e., Sw < 100%), however, introduces complications to the SBVI method. To overcome the limitations of both methods for determining BVI, the recommended practice is to compute two bound-water values—one from the cutoff (CBVI) and one from the spectral (SBVI) method—and take the larger of the two.

Permeability Estimation

The ability to estimate formation permeability is one of the earliest benefits of NMR logging and remains the most important application. Laboratory studies demonstrate that pore-water relaxation time is inversely related to the S/V ratio of the pore space (Fig. 3E.26). The NMR estimate of permeability is based on theoretical and core-based models that show that permeability increases with increasing porosity and pore size (S/V).[80][87][88][89][90]

The measurement of formation permeability, in general, is greatly influenced by the method used, the limitations of each, and the scale at which the measurements are made.[91] As stated previously, MICP curves obtained on core samples correlate to pore-throat size, while NMR measures pore-body size.

NMR logging does not provide direct and continuous measurement of permeability; rather, a formation-permeability estimate, or index, is calculated from the spectral-porosity measurements using permeability models that are based on a combination of empirical and theoretical relationships. Several permeability models have been developed, and two are in common use: the free-fluid (Timur-Coates or Coates) model and the mean-T2 [the Schlumberger-Doll-Research (SDR)] model.[27][92][93][94] The free-fluid model can be applied to water-saturated and hydrocarbon-saturated reservoirs, and the mean-T2 model can be applied to water-saturated reservoirs.[95] These permeability models assume that a good correlation exists between porosity, pore-body and pore-throat size, and pore connectivity. This assumption is generally valid in clastic (e.g., sand/shale) sequences, but in carbonates or other lithologies, model-derived permeabilities may not be reliable.

Typically, a permeability model is calibrated over a particular zone of interest and verified, wherever possible, by core or well/formation test data. Once this is done, the NMR log can provide a robust continuous-permeability estimate within the zone of interest.

Both models treat permeability as an exponential function of porosity, ϕ4, and include a parameter to account for the fact that NMR measures pore-body size, not pore-throat size[20] (Fig. 3E.27). In the Coates model, the pore-size parameter enters implicitly through T2cutoff, which determines the ratio of FFI to BVI. In the SDR model, the size parameter enters through the geometrical mean of the relaxation spectra, T2gm. In water-saturated rocks, both models provide similar and good results; however, when hydrocarbons are present, the SDR model fails because T2gm is no longer controlled exclusively by pore size.[96]

Free-Fluid (Timur-Coates or Coates) Model. In the simplest form of the free-fluid model, permeability, kCoates, is expressed as follows (Eq. 3E.14):

RTENOTITLE....................(3E.14)

where ϕ is MSIG, MBVI is obtained through the CBVI or SBVI method, MFFI is the difference between MSIG and MBVI (assuming that there is no clay-bound water, see Fig. 3E.16), and C is a formation-dependent variable. The free-fluid model is very flexible and has been calibrated using core data for successful use in different formations.

To calibrate the model to core, Eq. 3E.14 is solved in the form of a straight line, y = mx + b:

RTENOTITLE....................(3E.15)

Assuming b = 0 in the equation (3E.15), core permeability is substituted for k. The slope of the line, m (i.e., C value in Eq. 3E.15), is determined using a least-squares regression (Fig. 3E.28).

Despite the flexibility of this model there are formation conditions that limit the effectiveness of the model and may require a correction (Table 3E.5). The presence of hydrocarbons (i.e., oil, oil filtrate, or gas) in the BVI component may result in an overestimate of BVI by either the CBVI or SBVI methods, leading to an underestimate of permeability. An HI correction can be applied when gas is present. The very short T2 values associated with heavy oil may be counted in the BVI component and result in an underestimate of permeability.

Mean-T2[Schlumberger-Doll-Research (SDR)] Model. Using the SDR model, permeability is expressed as

RTENOTITLE....................(3E.16)

where ϕ is NMR effective porosity (MPHI), T2gm is the geometric mean of the T2 distribution, and C is a formation-dependent variable. The SDR model works very well in water-saturated zones. In the presence of oil or oil filtrates, the mean T2 is skewed toward the T2bulk, because of the effects of partial polarization, leading to an overestimate of permeability. In unflushed gas zones, mean-T2 values are too low relative to the flushed-gas zone; and permeability is underestimated. Because hydrocarbon effects on T2gm are not correctable, the SDR model fails in hydrocarbon-bearing formations. The Coates and SDR models represent matrix permeability and, therefore, are not applicable to estimation of permeability in fractured formations.

Other Comments on NMR-Based Permeability Estimation. Table 3E.5 compares the Coates and SDR models under different reservoir conditions, and it may be advisable to use both methods in an effort to constrain values for permeability.

There are a number of benefits in having available NMR-derived permeability and BVI. This information enables more-accurate quantification of reservoir heterogeneity and improves estimation of reserves and ultimate recovery. Other applications include optimizing the locations of perforations, well spacing, tailoring completions to maximize recovery rates and efficiencies, and improving primary and secondary recovery design schemes.

Hydrocarbon (Fluid) Typing

The NMR T2-porosity relationship in which T2 is a function of pore size (i.e., S/V ratio, see Eq. 3E.8) holds for water-saturated rocks. The presence of hydrocarbons in water-wet rocks alters the T2 distribution, thus affecting the porosity interpretation. Despite the variability in the NMR properties of fluids, the locations of signals from different types of fluids in the T2 distribution can often be predicted or, if measured data are available, identified (Fig. 3E.29).

Hydrocarbon typing and prediction of fluid properties by NMR logs is predicated on reliable laboratory correlations between NMR measurements (i.e., relaxation times and diffusion) and fluid properties [e.g., specific gravity, viscosity, and gas/oil ratio (GOR)]. Early studies were limited to investigations at ambient conditions; [97] however, using the standard correlations derived from these studies may result in seriously underestimating viscosity. More-recent studies have expanded these correlations to oils and mud filtrates at reservoir conditions.[97][98][99][100]

NMR logging uses specialized CMPG pulse-acquisition sequences to exploit these differences in pore-fluid NMR properties to achieve specific objectives. The CMPG parameters—wait time (TW), echo spacing (TE), the number of echoes (NE), and the number of sequence repetitions—are selected to take advantage of the wide variation in NMR fluid properties (Fig. 3E.30) for estimating total porosity and for hydrocarbon typing. Table 3E.2 illustrates the range of NMR-related properties of fluids for Gulf of Mexico sandstones, for example. Table 3E.6 lists variations in CMPG pulse sequences for different objectives using time-domain analysis (TDA) and Enhanced Diffusion Method (EDM). The multifrequency tools now in service permit the simultaneous acquisition of multiple measurements on the same rock volume using different acquisition sequences during the same logging run. The general approach is to log in a mode that allows gathering the full spectrum of data. Specialized applications, including some direct fluid-identification methods, involve customized-acquisition sequences that require slower logging and acquisition of more echoes.

Determination of the appropriate interpretation method is largely based on the estimated viscosity of the anticipated hydrocarbons (see Fig. 3E.31). These methods will be explained in the following paragraphs.

Selection of the appropriate acquisition sequence and the choice of acquisition parameters depend on the logging objectives and are part of prejob planning. This process considers several factors regarding the anticipated rock and pore-fluid properties. For example, for characterizing large pores (clastics), clean formations, carbonates, and light oil—all of which are associated with long-T2 values[101]—a large number of echoes should be acquired, and a longTW should be used. However, more echoes and longer TW may require reduced logging speeds. In contrast, characterization of small pores (i.e., low permeability), shaly formations, low porosity, and BVI determination involves the short T2 component of the NMR-porosity spectrum and often it can be accomplished using fewer echoes and shorter TW, which may allow normal logging speeds.

Fluid-typing methods fall into two broad categories depending on the NMR properties that are being exploited: (1) T1-weighting mechanisms take advantage of differences in fluid T1 values, and (2) diffusion-weighted mechanisms make use of the diffusion-constant differences between oil and water. There are two general types, or sets, of CMPG acquisitions that are associated with each mechanism: dual TW and dual TE. These two sets cover the range of major fluid-typing objectives; some interpretation techniques can use one or the other, or both. Each serves specific purposes and is optimized to provide data for specialized analysis programs. In general terms, each consists of at least two echo-train acquisitions in which one or more parameters are varied. The total, or the difference between the echo-train signals, provides an estimate or indicator of porosity, light hydrocarbons, or oil (Fig. 3E.32).

Advanced hydrocarbon-typing objectives can involve customized-acquisition sequences[102] that are run at reduced logging speeds or even in stationary mode.

T1-Weighted Mechanism (Dual-TW Acquisition). The dual-TW acquisition method is used primarily to identify and quantify light hydrocarbons (gas and light oil) by separating them from the water signal through T1 weighting. The dual-TW acquisition can also provide the standard-T2 acquisition dataset for determination of spectral NMR porosity, permeability, and productivity (mobile fluids). The TE is kept the same in both echo-train acquisitions (0.9 or 1.2 ms), and a short TW (1 second) and long TW (8 seconds) are used. The water signal is contained in both acquisitions, but because of the light hydrocarbons (which have long T1 values), the signal is suppressed in the acquisition using the short TW. Thus, the presence of a signal in the difference of the measurements (differential spectrum) is considered a strong indicator of gas or light oils[103] (Fig. 3E.33).

Fluid volumes can be quantified by integrating the difference into T2 spectra and correcting for the polarization difference between long and short wait times. Fluid typing and quantification are performed through the differential-spectrum (DSM) method and the TDA method. Two conditions must be met to ensure successful DSM interpretation: T1 contrast between the hydrocarbon and brine phases (i.e., a water-wet formation containing light hydrocarbons), and T2 contrast between the gas and oil phases. DSM assumes that the hydrocarbon phases relax uniexponentially. DSM is commonly used for: (1) hydrocarbon typing in shaly sands, where the the neutron-density crossover may be suppressed because of the high amount of clay minerals in the rock[104][105] to confirm the presence of light oil in fine-grained rock, and (2) for gas detection in the presence of OBM invasion.[106] Ideally, DSM can be used to compensate for the hydrocarbon effects on NMR measurements and thereby enable correction of NMR total porosity and NMR effective porosity; however, because of S/N requirements, TDA is the preferred technique for correcting NMR logs.

In contrast to the DSM approach, the TDA technique works directly with the echo-decay data (i.e., in the time domain rather than the T2 domain; see Fig. 3E.34). The measured long- and short-wait-time decay trains are subtracted into an echo difference that indicates the presence of a long-T2 component (usually a hydrocarbon). By working in the time domain, when the T2 inversion is performed, only the hydrocarbon T2 component will be present.

TDA—an outgrowth of the DSM technique—is based on the fact that different fluids have different rates of polarization (i.e., different T1 relaxation times).[107] TDA provides the following results:

  • Fluid types in the flushed zone.
  • Corrected NMR porosity in gas reservoirs.
  • Corrected NMR porosity in light oils.
  • Complete fluid-saturation analysis in the flushed zone, using only NMR data.


The TDA technique is a more robust method than DSM, in part because it provides better corrections for underpolarized hydrogen and HI effects[108] (Figs. 3E.34 and 3E.35).

Several factors inherent in the dual-TW acquisition require reduced logging speeds: to achieve full polarization (long TW) in one acquisition channel, to acquire the small signal amplitudes associated with the short TW values in another other channel, and to meet the requirement for acceptable S/N levels. A new triple-wait-time method addresses these issues and also allows T1 acquisition.[109]

Diffusion-Weighted Mechanism (Dual-TE Acquisition). NMR-logging tools generate relatively large field gradients, which are a function of the operating frequency and tool type. If the pore fluids have high diffusivity constants (D0), these large gradients may cause diffusion to become the dominant T2-relaxation mechanism, even at a reasonably short interecho time. The T2 contrast between the liquid and gas phases can then be exploited for hydrocarbon typing.

The identification of viscous oil in water-wet rock using T1-weighting mechanisms (dual-TW acquisition) is very difficult because there is little T1 contrast between water and viscous oil. However, a significant difference in the diffusion characteristics of water and viscous oil is used to create a T2 contrast that separates the two NMR signals, thereby permitting identification.[110] The inter-echo spacing, TE, is chosen long enough that the water and oil signals are fully separated in the T2-domain (i.e., water is at lower T2 than oil; see Fig. 3E.36).

Diffusion-weighting techniques, such as shifted-spectrum method (SSM), enhanced-diffusion method (EDM), and diffusion analysis (DIFAN), use the dual-TE acquisition to create this T2 contrast. TE is chosen such that diffusion, rather than surface relaxation, is the dominant relaxation mechanism, and thus, the upper limit of T2 for pore water is T2 diffusion, designated TDW. The dual-TE acquisition set uses long TW and TE values, typically 0.9 or 1.2 ms and 3.6 or 4.8 ms, respectively. In this acquisition, the fluid with the larger diffusion constant (water) has a spectrum shifted more to earlier times than the fluid with the smaller diffusion constant (viscous oil). The presence in the spectra of a minimally shifted portion identifies high-viscosity oil in the formation (see the lower portion of Fig. 3E.37).[110][111] These techniques are primarily used for identifying the presence of viscous oil. They are also used in carbonates in which a reduced surface relaxation (compared to clastics) may also result in a reduced T1 contrast.[112][113]

SSM. This method was the first developed using the diffusion-weighted mechanism and was originally developed for use in gas reservoirs as a qualitative technique (see the upper portion of Fig. 3E.37). However, because of signal-processing difficulties and slower logging speeds, this method was effectively replaced by the EDM.[103]

EDM. This method is the later variation of SSM. Depending on oil-NMR properties and the job objective, EDM data processing can be done in either the T2 domain or time domain.[110] If the EDM objective is to discriminate pay from nonpay, then a single CMPG measurement with long TW (for full polarization) and long TE (for diffusion enhancement) is sufficient; thus, standard-T2 logging with a long TE can be used. If the EDM objective is to quantify fluids in a pay zone, then dual-TE logging will be required. The short-TE measurements will provide correct MPHI and BVI. If the T2 contrast over the zone of interest is not expected to be large enough to separate the T2 values of water and oil, then dual-TW logging with a single, long TE may be required to obtain data for TDA processing.

One caveat for users of the EDM: the flushed-zone saturation (Sxo) obtained from EDM may not reflect the true volume of residual hydrocarbons in the flushed zone if the faster portion of the hydrocarbon’s T2 spectrum gets fully polarized by both wait-time measurements. A correction to Sxo can be applied if the full-hydrocarbon T1 spectrum can be determined by other means, such as laboratory measurements.

DIFAN. DIFAN is an empirical model used for quantitative-diffusion analysis. This method uses a dual-TE acquisition and is useful for typing and quantifying oils with viscosities ranging between 0.5 to 35 cp at reservoir conditions (i.e., temperature ≥ 200°F and pressure ≥ 2,000 psi). It is not recommended for light oils and condensates because the D 0 contrast between hydrocarbons and water is too small.

Neither EDM nor DIFAN is recommended for use in high-viscosity oil (heavy oil) because the difference between the T2 values of dead oil and irreducible water is too small. The DIFAN model generates two T2 distributions using the two echo trains generated from dual-TE logging.

Echo Ratio Method. An apparent diffusion coefficient, calculated from the ratio of long and short TE echo-train data (time domain), serves as a qualitative gas indicator.[114] This method is accurate and provides a direct measurement of gas volume, but requires high HI values for good results.

The dual-TW and dual-TE techniques described above only provide T2 distributions and are considered one-dimensional methods.[115] Oil-saturation determination using these methods may be difficult or even impossible because differential diffusion-weighted methods assume that NMR water and hydrocarbon signals can be cleanly separated by subtracting the acquisitions made at different TE values. In practice, however, partial—or even total—overlap is common.

To avoid these problems, new multifluid-forward models were developed to take advantage of the full suite of data obtained using the latest multifrequency logging tools. These models perform a simultaneous inversion of the NMR fluid properties obtained in multiple parameter-domain dimensions (e.g., T1, T2, and D0).[59][116][117][118][119][120] These methods make use of the enhanced T2 relaxation in a magnetic-field gradient (i.e., incorporate measurement of the diffusion coefficient, D0) and the large contrast between diffusion coefficients of oil and water. In effect, they are simultaneously combining both T1- and diffusion-weighted processing techniques, and thus, provide more accurate and robust results. In addition to porosity and permeability information, the output data from these inversion models are plotted as two-dimensional plots (e.g., D-T2 plot) in which the oil and water signals are clearly separated.[121] Recent two-dimensional NMR techniques include diffusion editing, diffusion mapping, and relaxation-diffusion 2D.[26][118][122][123][124][125]

Multidimensional (2D and 3D) NMR analysis greatly increases the accuracy of fluid typing and saturation determination.[115] These methods are particularly useful for identifying highly diffusive fluids (e.g., gas, condensate, and light oil)[59][125][126][127][128] and for identifying wettability alteration caused by OBM-filtrate invasion or chemical surfactants used in enhanced recovery operations.[118][127][129] This information is essential for reservoir simulation, development, and production.[130]

Residual Oil (Sxo) Calculation

Estimation of residual-oil saturation is one of the oldest applications of NMR logging. Unlike resistivity-log analysis, NMR analysis does not rely on formation-water salinity to obtain water saturation. This feature provides NMR logging with a distinct advantage over conventional resistivity analysis in mixed or unknown salinity conditions—an advantage that can be extremely useful in waterflood or steamflood projects for evaluating residual-oil saturation after the flood or to look for bypassed oil. Initially, NMR evaluation of residual-oil saturation required use of a dopant (e.g., Mn-EDTA and MnCl2) to flush the invaded zone with paramagnetic ions to shorten the bulk relaxation time of the brine. This process enabled separation of the oil and water signals, leading to direct measurement of Sxo. With the combination of modern multifrequency NMR tools and new methods of analysis, such as EDM and TDA, the use of borehole dopant is no longer necessary.[25][111] Standalone NMR interpretation is possible in OBM, but in water-base-mud (WBM), it is possible only when Swirr is known.

Viscosity Evaluation

In water-wet rocks, the NMR-relaxation spectrum is controlled by viscosity[118][131]; relaxation is directly related to viscosity, η, by Eq. 3E.17:

RTENOTITLE....................(3E.17)

For oil and water, the diffusion constant, D0, can be approximated by Eq. 3E.18:

RTENOTITLE....................(3E.18)

where A = 2.5 for water, A = 1.3 for oil, and TK is temperature in K.

NMR properties of gas can be obtained from published charts that relate viscosity to (1) the center of the relaxation curve, (2) to an American Petroleum Inst. (API) standard value,[103] or (3) from published formulas.[108] These published sources assume that methane is the dominant component of the gas. In the absence of laboratory data at in-situ conditions, reservoir NMR properties can be estimated by use of Eq. 3E.17. One note of caution: Eqs. 3E.17 and 3E.18 are based on "dead oil" measurements. As mentioned earlier, the relaxation-time dependence on viscosity/temperature of live crude oil may differ significantly from correlations based on hydrocarbon liquids at ambient conditions.[97][98][99]

Anisotropy and Geomechanics

NMR T1 and T2 relaxometry and derived permeability can be used to determine the Biot elastic constant, which is used for estimating pressures that are critical to sand control, hydraulic fracturing, wellbore stability, and determination of formation stress.[132][133]

Low-Permeability (Tight) Sandstones

Field experience indicates that invasion or imbibition in tight sands is very shallow. Depending on borehole size, the diameter of invasion, and fluid properties, a mandrel tool may have sufficient DOI to measure beyond the flushed zone. A combined interpretation using both T1 and T2 can provide positive identification of fluids in these reservoirs.[134] In addition, in these tight formations, in which formation testers typically may not obtain a fluid sample within a reasonable time period, NMR fluid characterization can separate hydrocarbon from oil filtrate and other pore fluids.[135]

Heavy Oil, Tar Sands, and Tar Mats

The early acquisition of reliable viscosity information is essential to efficient development of heavy-oil reservoirs. NMR logs offer a viable alternative to downhole fluid sampling for determining viscosity information in heavy-oil reservoirs.[136][137][138] The presence of tar mats in a reservoir, commonly near the bottom of the oil column, may form vertical permeability barriers and, thereby, isolate the oil leg from the water-drive aquifer. NMR logs, in conjunction with conventional logs, can provide accurate identification of tar-mat levels and viscosity estimation, from empirical relationships.[139] (See the Job Planning section of this chapter.)

Carbonates and Comples Lithologies

NMR-log evaluation is relatively routine in what might be termed "conventional reservoirs," namely those of homogenous lithology and uniform pore sizes, typically sandstone and chalk reservoirs.[62][140] In contrast, log evaluation of complex and heterogeneous reservoirs, with complex pore geometries, is not straightforward. These reservoirs include, in particular, the highly important Middle East carbonates, as well as other reservoirs comprised of mixed lithologies and mineralogies, or both, in which wettability may also vary. In these reservoirs, there is likely no simple relationship between petrophysical properties and porosity. Instead of a dependence on the volume of pore space, they are dependent on the typically heterogeneous pore distribution, pore types, pore connectivity, and grain sizes. This fundamental difference between siliciclastic and carbonate rocks (primarily the result of diagenetic processes) limits the applicability of routine NMR methods, especially permeability evaluation.[141] Improving NMR evaluation of carbonates has proved challenging and is the subject of a number of recent studies and proposed techniques.[120][141][142][143][144] (See the Job Planning section of this chapter.)

Pseudocapillary-Pressure Curves

As discussed in the Petrophysics section, establishing a correlation between NMR T2 distribution and MICP is fundamental to NMR interpretation and the computation of Sw. Once such correlations are established for a particular reservoir or field, pseudocapillary-pressure curves can be generated directly from the NMR-log relaxation-time distributions.[145][146][147][148]

Productibility

Accurate determination of BVI enables evaluation of reservoir-fluid (e.g., gas, oil, and water) contacts, production characteristics (producibilty), and the determination of net and recoverable reserves.[149][150] Furthermore, NMR-derived permeability can be used to generate idealized flow profiles across a completion interval. These profiles provide a diagnostic tool for identifying nonflowing portions of the zone and the need for remedial work.[151][152]

Combined NMR Applications


NMR tools have shallow depths of investigation and provide results only for the invaded zone. NMR-log data can be integrated with core data and conventional log data (i.e., neutron, density, acoustic, and resistivity) in post-acquisition processing to provide improved determinations of reservoir rock properties, hydrocarbon storage capacity, and reservoir productivity in a variety of environments including gas-bearing and low-resistivity reservoirs. Interpretation models that include NMR data can provide more reliable results than those using only conventional logs.

The combination of NMR and deep-resistivity data provides a complete analysis of the fluids in the uninvaded zone. Resistivity measurements alone cannot distinguish between capillary-bound and movable water, but they do represent deep investigations of fluid saturation. Furthermore, resistivity-based methods are often inadequate or unreliable in reservoirs in which salinity and lithology vary. The addition of NMR-derived BVI and CBW from the flushed or invaded zone can significantly enhance the estimation of resistivity-based fluid saturation, both in clastics and carbonates. The addition of NMR data allows the identification and evaluation of water-free production in low-resistivity formations (Fig. 3E.38).[105][153][154]

The combination of conventional deep-resistivity data with NMR-derived CBW, BVI, FFI, and MPHI can greatly enhance petrophysical estimations of effective pore volume, water cut, and permeability. It is the preferred technique for identifying low-resistivity pay zones. Fig. 3E.38 presents the results of Halliburton’s MRI analysis (MRIAN) service in a turbidite sequence. The sand below XX200 depth has an average resistivity of approximately 0.5 ohm-m (Track 2) and average neutron-density porosity of approximately 38% (Track 4). A quick look, or preliminary analysis, using only the conventional data presented in Tracks 1, 2, and 4, would label this a wet zone. MRIL-NMR data are also presented in Track 1 (T2-distribution bin data), Track 2 [NMR-derived permeability (Coates)], and Track 3 (VDL presentation of T2 distribution). The MRIAN results, BVI, and FFI are presented in Track 4. BVI gradually increases with depth, suggesting that the sand is fining downward (i.e., as sand grains become finer, the volume of capillary-bound water that they hold increases). Comparison of the BVI with the resistivity profile (Track 2) shows that the resistivity decreases where the bound water increases. The MRIAN analysis clearly shows that the zone does not contain any movable water and will produce only oil. The interval below XX200 was tested and produced oil with no water.[32]

The combination of NMR fluid identification with resistivity-derived saturation provides a better understanding of hydrocarbon movement in gas monitoring, as seen in the following field example (Fig. 3E.39).

NMR and LWD logs were run in a deviated Middle East light-oil carbonate reservoir to monitor gas injection. NMR logs were well suited for gas monitoring by use of the TDA technique and were not affected by formation-water salinity.[31] NMR gas-corrected porosity was in close agreement with core porosity; conventional log-porosity was low because of salinity effects. Combined interpretation of water-saturation analyses from tools with different depths of investigation identified the movable fluid present in the reservoir as WBM filtrate. Track 1 contains T2 bin data and conventional SP and gamma ray curves; Track 2 shows NMR-permeability derived from the standard Coates model and LWD resistivity; Track 3 contains a T2 distribution in a variable-density format; Track 4 contains T2 distributions of data acquired with long TW and short TW; Track 5 is the differential spectrum; Track 6 contains TDA results; and Track 7 contains MRIAN results.

NMR Acoustic/Density Combination

Because NMR tools are affected by the lower HI in gas-bearing reservoirs, a moving tool may not fully polarize the gas, which has a long T1. Conventional acoustic and density logs are not influenced by these factors and, when used jointly with NMR logs, they can provide robust porosity evaluation.[155][156] The basic assumptions of these techniques are that the reference porosity is measured correctly (e.g., appropriate matrix density is used) and that the NMR TW is long enough to recover the water but not all of the gas. If the TW is too short, a false gas signal may be seen. One way to check whether the TW is appropriate is to log with a longer wait time (e.g., TW = 2 seconds rather than 1 second) and look for a significant difference in porosity (i.e., the porosity may be slightly higher because of increased polarization of the gas).

NMR-Log Quality Control


Quality control (QC) procedures are especially important during the creation of integrated computed products[37][78][157] and to ensure optimal NMR data acquisition. NMR tools have calibration standards and real-time QC indicators. Standard techniques include examination and verification of service company calibrations and QC curves, and calibrations and examination of the repeatability of the computed porosities. Processed results should agree with other data from logs, core, and/or formation test results. NMR QC includes a series of prejob and post-job checks and calibrations.

Prejob Calibration and Quality Checks During Logging

Several prejob QC steps are necessary to ensure reliable results. These steps are described in this section.

Frequency Check. The amplitude of the NMR measurement is proportional to the square of the magnetic field strength in the measurement region. To attain maximum signal strength, the tool must operate at the correct frequency. Because temperature affects magnetic field strength, a frequency sweep is recorded before the calibration and logging run and is repeated downhole to account for changes in temperature. The CMR tool includes sensors that monitor temperature and changes in the magnetic field strength and adjusts frequency automatically to changes in sonde temperature. The accumulation of metal drilling debris on the magnets can adversely affect the NMR measurements, and the tool may require occasional retuning. The MRIL-gradient field creates a self-correcting effect, and the tool rarely requires retuning after the initial downhole tuning. Drifting off frequency may present itself as a loss of porosity (CMR) or loss of precision (MRIL). Field calibration should be included on the log.

Statistical Calibration. A statistical calibration to 100% porosity (water) should be made at the shop and/or wellsite before and after every logging job for all combinations of TE frequency and for expected quality factor (Q level); calibration should be done at the shop at least monthly. Primary calibration involves either a flask of fluid (for the CMR tool) or a water tank (for the MRIL tool).

HI. NMR-tool calibration uses water where HI = 1. However, when the HI of the formation pore fluids is expected to be less than one, the NMR-porosity readings will be proportionally lower, and correction is necessary. An NMR tool’s processing software will, therefore, compute a correction based on the salinity of the mud filtrate and on the maximum formation pressure and temperature. In general, a mud-filtrate relaxation time of < 200 ms may suppress the long T2 components in a T2 distribution. A specific HI value may be used when residual hydrocarbons or OBM filtrates are present.

System Gain. CMR and MRIL tools monitor and calibrate in real time for downhole changes in system gain. A gain measurement is made as a part of each pulse sequence. Gain indicates the amount of loading applied to the NMR tool’s transmitter circuit by borehole fluid and formation resistivity. Gain is affected by changes in temperature, mud conductivity, and borehole size. Because gain is frequency dependent, the operating frequency of a tool should be set to achieve maximum gain. Abrupt changes or spikes in gain should not be present. For a particular MRIL tool, the gain value determines the appropriate acquisition power (Q) level.

System Noise (Ringing). Noise contamination of the spin-echo signal, known as ringing, may be a remnant of the RF pulsing process. Noise can interfere with the echo trains and, when present, is usually evident in the NMR porosity readings, either as poor repeatability or lack of agreement with other tools. Both tools (CMR tool and MRIL tool) monitor ringing.

χ. χ is a measure of the quality of fit between the calculated decay curve and the recorded echo amplitudes. Problems with the echo data are likely to be reflected in the χ curve. In general, the value of χ should be less than two but may average slightly higher in certain situations. Spikes in χ may correlate with spikes in porosity[78] and usually indicate tool problems, even if χ remains lower than two. χ serves as is a primary MRIL log-quality indicator and is monitored while logging.

Correction for RF-Tipping Pulse. The strength of the CMPG RF pulse (B1) that produces proton tipping and rephasing is measured as part of every pulse sequence and requires correction for changes in borehole temperature. If the pulse angles are either < 90° or > 90°, the magnetization will be undertipped or overtipped, respectively. The measured amplitude will, thus, be too small, and porosity will be underestimated. The B1 curve should be relatively constant but should show some variation with changes in borehole and formation conductivity. B1 will decrease across conductive washouts and conductive formations. Changes in the B1 values should track changes in total conductivity and vary together in the same direction as gain. The need for excessive correction (i.e., > 5% of the optimum shop-peak value) may indicate a tool problem and can result in undertipping or overtipping of the protons, a reduced S/N ratio, and a loss of precision in determining porosity. Sudden changes in the B1 correction curve may also indicate a tool problem.

γ. Regularization methods are used to select a smooth T2 distribution that is consistent with the spin-echo sequence. These methods require a parameter γ that is automatically computed from the raw echo data. Values of γ are dependent on the S/N ratio and the shape of the underlying T2 distribution. In high-S/N environments (i.e., medium-to-high-porosity formations), typically γ < 5; in low-S/N environments (i.e., tight sands and shales), γ > 10.

Repeatability. Whenever possible, a repeat pass should be recorded with parameters identical to those used in the main pass. These parameters include TW, NE, and TE, as well as computation parameters such as the echoes selected for processing, and for T2cutoff. The generally accepted goal for porosity is a standard deviation of 1 porosity unit (p.u.). Repeatability of BVI is usually > 1 p.u., and repeatability of FFI is usually >> 1 p.u.

If repeatability is a concern, decreasing logging speed or increasing the degree of echo data stacking during logging in post-job processing can improve S/N ratio and repeatability. Typically, fluid-typing applications using dual-TW or dual-TE methods are more sensitive to data repeatability than a standard NMR log acquired for porosity, bound fluid, and permeability.

A given data set should agree with similar data acquired by other logs, formation tests, and/or core analysis.

Environmental Corrections. NMR logs are similar to neutron logs in that they respond to the hydrogen volume present in a sample. Because the hydrogen volume changes with temperature and pressure, NMR logs require environmental corrections for temperature and pressure similar to those applied to neutron logs. Also, the magnitude of NMR resonance (and S/N ratio) varies inversely with temperature. NMR-logging tools include a temperature sensor to acquire the data needed for this correction. At high temperatures, data stacking can be increased or logging speeds reduced to compensate for a lower S/N ratio.

Porosity Check. NMR porosity can serve as an important diagnostic of data quality. High porosity readings may result from washouts (e.g., borehole fluid in the sensed volume), tool problems, improper environmental corrections or calibration issues, loss of pad contact (for the CMR tool), tool eccentering, or borehole ellipticity (for the MRIL tool). Low porosity readings may result from insufficient wait time, light-hydrocarbon effects, the presence of heavy oil, improper frequency tuning, and/or calibration error.

Post-Logging Quality Check

NMR-log responses should be checked against conventional logs when they are available (see the Applications section of this chapter).

NMR-Log Job Planning


In some complex reservoirs, low-resistivity/low-contrast pay, low-porosity/low-permeability, and medium-to-heavy oil, NMR-log data—independently or in combination with other log data—provide the best and/or only means of accurate formation and fluid evaluation. Because NMR-log data acquisition is complex, job preplanning is essential to ensure optimal selection of acquisition parameters that will result in reliable and accurate data and in the maximum information possible in any given reservoir and logging environment. A clear understanding of the logging job objectives is necessary for optimizing the NMR acquisition parameters to best achieve these objectives.[31][122][157][158] This process must take place before the actual logging.

Typical preplanning consists of three steps:

  • Define the need for NMR measurements.
  • Collect all available borehole (e.g., diameter, mud, salinity, and temperature) and reservoir (e.g., formation and fluid properties) information needed to assess the expected NMR responses in the zone of interest, and understand what can and cannot be resolved with NMR.
  • Select the appropriate tool (on the basis of operational considerations, borehole size, and condition) and acquisition type (i.e., determining the appropriate acquisition parameters, data resolution, and logging speed) that will provide maximum answers for a given job.[31]


Although the actual in-situ reservoir characteristics may be unknown, estimates of the anticipated fluid properties, based on available information such as reports for nearby wells or fields, are used to define and optimize an acquisition sequence that will provide the data needed to meet the job objectives.

In addition to job objectives, determination of the appropriate NMR-acquisition parameters is also influenced by operational considerations and the anticipated in-situ reservoir properties (Fig. 3E.40). The most critical factors follow.


Lithology

Although reservoir lithology generally plays a minor role in NMR-data acquisition, it does play a significant role in data analysis and interpretation. Aspects of reservoir lithology that influence reservoir T2 values include the following.

Carbonates. Different T2cutoff values are required because surface relaxivity in carbonates is weaker than in sandstones, resulting in slower relaxation rates (longer T2). Longer T1 in carbonates than in sandstones may require longer TW during acquisition.[101] (See the Application section of this chapter.)

Isolated Pores. The presence of relatively isolated pores (e.g., vugs) will not affect NMR porosity, but it will cause the standard permeability equations (Coates and SDR) to overestimate permeability. (See the Permeability section of this chapter.)

Ferromagnetic and Paramagnetic Minerals. The presence of these minerals may enhance surface relaxation significantly, shifting the T2 spectrum to very short relaxation times. Depending on the amount of paramagnetic material, relaxation may become too fast to be detected, and the NMR measurement will underestimate porosity. In these cases, standard cutoff values do not apply.[159]

Heavy Oil and Tar Sands. Intervals containing these types of hydrocarbons have very fast relaxation components and may not be detected using conventional acquisition methods. Special methods have been developed for detection and accurate evaluation of these reservoirs.[136][137][160][161][162][163][164][165][166][167] (See the Applications section of this chapter.)

Wettability

Wettability can have a significant impact on NMR-log response.[168] The use of NMR for determination of wettability has been extensively studied both in laboratory[129][169] and in the field. In general, petrophysical-NMR studies and NMR-logging applications assume reservoir rocks are water wet; however, because mixed-wettability reservoirs do exist (e.g., some carbonates, black shales, and heavy-oil reservoirs),[112] this assumption may lead to incorrect reserve estimates and to unexpected dynamic behavior during waterflood. When a pore is water wet, oil relaxes at its bulk rate. In mixed-wettability reservoirs, the oil and water each relax through a combination of bulk relaxation and surface interaction and depend on the ratio of water-wet surface area to water volume and oil-wet surface area to oil volume. The oil relaxation spectra will be shifted from bulk relaxation into the irreducible-water part, resulting in complex spectra that are difficult to interpret.[170][171] Nevertheless, these shifts provide qualitative wettability indicators that allow NMR logs to provide an early indication of reservoir wetting behavior.[113][118][172][173][174]

The invasion of OBM or SOBM can alter formation wettability and is a significant concern in NMR logging,[113][131][175] which typically measures fluid and formation properties in the flushed zone. Invasion by these muds can alter strongly water-wet sandstones and carbonates to intermediate-wet or oil-wet rocks. In water-wet reservoirs that have undergone OBM or SOBM invasion, the T2cutoff model may significantly underestimate Swirr because wettability alteration changes the water and oil relaxation-time distributions. The magnitude of underestimation depends on the type of OBM surfactants, their concentration in the flushing fluid, and the flushing volume. Controlling the volume of OBM invasion and the concentration of OBM surfactants should minimize the effects of OBM invasion on estimation of Swirr.[176][177]

Borehole Rugosity

NMR-logging tools have relatively shallow DOIs. Pad-type tools (e.g., the CMR tool) are run eccentered and require good contact with the borehole wall for accurate measurements. Measurements made by pad-type tools can be significantly affected by severe borehole rugosity and washouts, resulting in overestimation of porosity.

Mandrel devices are run centered, and DOI can range from 1 to 4 in., depending on borehole size. The sensitive area is normally beyond minor borehole rugosity. When borehole conditions result in an elliptical borehole (e.g., breakouts or erosion) or otherwise inhibit or prevent centralization (e.g., in highly deviated or horizontal boreholes), contact tools may be a better choice if pad alignment and contact with the borehole wall can be maintained.

Mud Type

The quality of NMR data acquired in OBM is generally superior to that acquired in WBM. The conductivity of OBM is lower; lower conductivity reduces loading effects on the transmitter/receiver system, resulting in higher S/N ratio. This issue is of greater concern for mandrel tools (e.g., MRIL tool) because, in conductive muds, the power of the RF pulse is reduced as the pulse is transmitted across the borehole. The use of fluid-excluding sleeves can minimize this problem.

Signal dissipation in conductive mud is not a serious concern for pad-type tools (CMR) that maintain contact with the borehole wall. Because NMR tools read the flushed zone, OBM filtrate invasion produces an additional hydrocarbon signal that may significantly complicate log interpretation.[178] Careful prejob planning can reduce interference of the OBM-filtrate signal and the response from the native fluids. The relatively long T1 relaxation times and diffusivity of OBM make it difficult to differentiate its signal using the shifted-spectrum or differential-spectrum approaches.

Metal Debris

Metal drilling debris in the borehole fluid may affect NMR-measurement quality by distorting and altering the logging tool’s magnetic field. The pad-type tool is more susceptible to field distortion. Metal debris should be removed from the mud, either through the use of the prepolarizing magnets (included in the latest tool designs) or by using magnets at the shale shaker. A new wireline tool uses an autotuning feature to correct the operating frequency for changes in the static magnetic field caused by metallic debris in the borehole.[55]

Logging Speed and Running Average

The logging speed of an NMR tool is influenced by a number of factors, primarily by tool type (e.g., centered or eccentered, number of operating frequencies, or length of antenna), logging objectives (e.g., acquisition type—TW, TE, or NE—sequence repetitions and vertical resolution), and borehole properties (e.g., diameter and mud resistivity). S/N ratio is primarily controlled by borehole size and mud resistivity. As S/N ratio decreases, the running average (RA) needed to maintain a specified error in porosity increases. The general practice is to require a porosity standard deviation of ≤ 1 p.u. The value of RA and the antenna aperture (i.e., length), combined with the logging speed, determines the vertical resolution. Even so, there is always a complex tradeoff in logging speed, accuracy (e.g., S/N ratio or NE), and job objectives. High accuracy and precision require reduced logging speeds. A method for increasing the overall logging speed is to reduce the vertical resolution in zones of secondary or no interest and also to reduce the vertical resolution in homogeneous intervals.

Vertical Resolution

The vertical resolution of NMR-logging tools is primarily a function of antenna length (i.e., tool design) and logging speed. The maximum vertical resolution, usually obtained when the tool is at rest (e.g., in stationary mode), is the length of the antenna. During continuous logging, vertical resolution decreases at a rate proportional to logging speed. Contact-logging tools, in general, use smaller sensors and antennae and, thus, have better vertical resolution than centered tools. The contact-NMR tool (i.e., the CMR tool) has a resolution advantage when bed thickness is in the range between 0.5 and 5 ft. Outside of this range, both designs deliver similar results. Prejob planning includes selecting a logging speed to obtain the optimum resolution.

Post-job data reprocessing to enhance bed resolution may result in a loss of repeatability. Vertical resolution can also be improved by optimizing the NMR signal through the removal of signal noise during data processing.[179][180]

Summary


This chapter has outlined the fundamental properties that are measured by NMR tools and has reviewed the use of these measurements to discern various characteristics of the reservoir rock including its fluid contents. Whether used as a standalone service or in combination with other logs and core data, NMR logs can provide an improved understanding of reservoir petrophysics and producibility. However, NMR logs are the most complex logging service introduced to date and require extensive prejob planning to ensure optimal acquisition of the appropriate data needed to achieve the desired objectives. Review of the references cited in this chapter indicates the rapid pace of advancement in NMR logging research, development, and applications. Users of NMR logging, including engineers, log analysts, petrophysicists, and geologists, should anticipate new developments in this discipline, as reported by SPE, the Soc. of Petrophysicists and Well Log Analysts (SPWLA), and other publications noted in the references.

NMR-Tool Mnemonics


Table 3E.7 presents a cross reference of logging-tool output-data mnemonics for the different service companies currently offering NMR-logging services.

Nomenclature


A = pore-fluid-specific value used to approximate the diffusion constant
a = constant in the mean-T2 and viscosity relationship
B0 = static magnetic field, gauss
B1 = amplitude of the oscillating magnetic field perpendicular to B0, gauss
C = coefficient in the Coates permeability model
D0 = molecular diffusion coefficient, gauss/cm
f = Larmor (precessional) frequency, Hz
G = field-strength gradient, gauss/cm
k = permeability, darcy
kCoates = permeability derived using the Timur-Coates model, darcy
kSDR = permeability derived using the mean-T2 model, darcy
M = magnetization, gauss/cm3
M0 = macroscopic magnetization, gauss/cm3
M0i = magnitude of the initial magnetization from the ith component, gauss/cm3
M0x = magnitude of the transverse magnetization at t = 0, gauss/cm3
M100% = magnitude of the magnetization for 100% bulk water, gauss/cm3
Mi(0) = magnitude of the initial magnetization from the ith component of relaxation gauss/cm3
M(0) = magnitude of the initial magnetization, gauss/cm3
M(t) = measured magnetization at time t, gauss/cm3
Mx(t) = transverse magnetization at time t, gauss/cm3
Mz = strength of magnetic field, gauss/cm3
Mz(t) = longitudinal magnetization at time t, gauss/cm3
NE = number of echoes
Pc = capillary pressure, dynes/cm2
Q = quality factor of a resonant circuit
Rt = true formation resistivity, ohm-m
(S/V)i = ratio of pore surface (S) to fluid volume (V), of the ith pore, 1/cm
(S/V)pore = ratio of pore surface (S) to fluid volume (V), 1/cm
Sw = water saturation, %
Swirr = irreducible water saturation, %
Swb = clay-bound water saturation, %
Sxo = flushed-zone saturation, %
t = time, seconds
T = temperature, °C
T1 = longitudinal relaxation time, seconds
T1bulk = pore-fluid bulk-T1 relaxation time, seconds
T1surface = pore-surface T1 relaxation time, seconds
T2 = transverse relaxation time, seconds
T2bulk = pore-fluid bulk-T2 relaxation time, seconds
T2cutoff = T2 cutoff value, seconds
T2diffusion = pore-fluid T2 relaxation time in a magnetic field gradient, seconds
T2DW = upper limit of the measured T2 for water, seconds
T2gm = T2 geometric mean value, seconds
T2i = pore-fluid surface T2 relaxation time of the ith component, seconds
T2surface = pore-fluid surface T2 relaxation time, seconds
TDW = dual wait time—the value of TE such that diffusion, rather than surface relaxation, is the dominant relaxation mechanism, seconds
TE = CMPG interecho spacing, seconds
TK = absolute temperature, K
Tp = pulsing time, seconds
TW = polarization (wait) time, seconds
TWL = polarization time, long, seconds
TWS = polarization time, short, seconds
Vshale = shale volume, %
x, y, z = cartesian space coordinates
γ = gyromagnetic ratio—the ratio of the magnetic dipole moment to the mechanical angular momentum, Hz/gauss
ΔTE = incremental change in echo spacing, seconds
ΔTW = incremental change in wait time, seconds
η = fluid viscosity, cp
ξ = apparent T1/T2 ratio of fluid
ρ1 = T1 surface relaxivity, cm/sec
ρ2 = T2 surface relaxivity, cm/sec
τ = time over which an oscillating field is applied
ϕ = porosity, %
ϕe = effective porosity, %
ϕeff = effective porosity, %
ϕH = hydrocarbon porosity, %
ϕi = calibrated porosity associated with all pores of the i th pore size
ϕt = total porosity, %
ϕw = water porosity, %
χ = goodness of fit


Acknowledgements


The authors of this chapter wish to acknowledge the help provided by the many former and current staff of Halliburton Energy Services who have been involved with the development of NMR logging. All figures except Fig. 3E.21 were provided courtesy of Halliburton Energy Services.

References


  1. 1.0 1.1 1.2 1.3 1.4 1.5 Brown, R.J.S. 2001. The Earth's-field NML development at Chevron. Concepts in Magnetic Resonance 13 (6): 344-366. http://dx.doi.org/10.1002/cmr.1020
  2. 2.0 2.1 Chandler, R. 2001. Proton free precession (Earth's‐field) logging at Schlumberger (1956–1988). Concepts in Magnetic Resonance 13 (6): 366-367. http://dx.doi.org/10.1002/cmr.1021
  3. 3.0 3.1 Woessner, D.E. 2001. The early days of NMR in the Southwest. Concepts in Magnetic Resonance 13 (2): 77-102. http://dx.doi.org/10.1002/1099-0534(2001)13:2<77::aid-cmr1000>3.0.co;2-c
  4. 4.0 4.1 Jackson, J.A. ed. 2001. Special Issue: The History of NMR Well Logging. Concepts in Magnetic Resonance 13 (6): 340-411.
  5. Brown, R.J.S. and Gamson, B.W. 1960. Nuclear Magnetism Logging. J Pet Technol (August 1960): 210-209. SPE-1305-G.
  6. Nikias, P.A. and Eyraud, L.E. 1963. Some Examples of Nuclear Magnetism Logging in Three San Joaquin Valley Oil Fields. J Pet Technol 15 (1): 23-27. SPE-434-PA. http://dx.doi.org/10.2118/434-PA
  7. 7.0 7.1 Kleinberg, R.L. and Jackson, J.A. 2001. An introduction to the history of NMR well logging. Concepts in Magnetic Resonance 13 (6): 340-342. http://dx.doi.org/10.1002/cmr.1018
  8. Kleinberg, R.L. 2001. NMR well logging at Schlumberger. Concepts in Magnetic Resonance 13 (6): 396-403. http://dx.doi.org/10.1002/cmr.1026
  9. Becker, E.D., Fisk, C.L., and Khetrapal, C.L. 1996. The Development of NMR. In Encyclopedia of Nuclear Magnetic Resonance, ed. D.M. Grant and R.K. Harris, Vol. 1. New York: John Wiley & Sons.
  10. Kleinberg, R.L. 1999. Nuclear Magnetic Resonance. In Methods in the Physics of Porous Media, ed. P. Wong, No. 35, Chap. 9, 337–385. San Diego, California: Experimental Methods in the Physical Sciences Series, Academic Press.
  11. 11.0 11.1 11.2 11.3 11.4 11.5 Dunn, K.-J., Bergman, D.J., and LaTorraca, G.A. ed. 2002. Nuclear Magnetic Resonance—Petrophysical and Logging Applications, Vol. 32. New York: Handbook of Geophysical Exploration: Seismic Exploration, Pergamon Press.
  12. Hahn, E.L. 1950. Spin Echoes. Physical Review 80 (4): 580-594. http://dx.doi.org/10.1103/PhysRev.80.580
  13. Carr, H.Y. and Purcell, E.M. 1954. Effects of Diffusion on Free Precession in Nuclear Magnetic Resonance Experiments. Physical Review 94 (3): 630-638. http://dx.doi.org/10.1103/PhysRev.94.630
  14. Meiboom, S. and Gill, D. 1958. Modified Spin-Echo Method for Measuring Nuclear Relaxation Times. Rev. Sci. Instrum. 29 (8): 688-691. http://dx.doi.org/10.1063/1.1716296
  15. Woessner, D.E. 1996. In a Maze in NMR. In Encyclopedia of Nuclear Magnetic Resonance, D.M. Grant and R.K. Harris, Vol. 1, 700. New York: John Wiley & Sons.
  16. Griffin, D.D., Kleinberg, R.L., and Fukuhara, M. 1993. Low-frequency NMR spectrometer. Meas. Sci. Technol. 4 (9): 968. http://dx.doi.org/10.1088/0957-0233/4/9/009
  17. Taicher, Z., Coates, G., Gitartz, Y. et al. 1994. A comprehensive approach to studies of porous media (rocks) using a laboratory spectrometer and logging tool with similar operating characteristics. Magn. Reson. Imaging 12 (2): 285-289. http://dx.doi.org/10.1016/0730-725X(94)91537-7
  18. 18.0 18.1 Morriss, C., Vinegar, H., Rossini, D. et al. 1997. Core Analysis By Low-field Nmr. The Log Analyst 38 (2). SPWLA-1997-v38n2a3.
  19. Mirotchnik, K., Kryuchkov, S., and Strack, K. 2004. A Novel Method to Determine NMR Petrophysical parameters from Drill Cuttings. Presented at the SPWLA 45th Annual Logging Symposium, Noordwijk, The Netherlands, 6–9 June. SPWLA-2004-MM.
  20. 20.0 20.1 20.2 20.3 Kenyon, W.E. 1997. Petrophysical Principles of Applications of NMR Logging. The Log Analyst 38 (2): 21–43.
  21. Murphy, D.P. 1995. NMR logging and core analysis -- simplified. World Oil 216 (4): 65-70. OSTI ID 39931.
  22. Lonnes, S., Guzman-Garcia, A., and Holland, R. 2003. NMR Petrophysical Predictions on Cores. Presented at the SPWLA 44th Annual Logging Symposium, Galveston, Texas, USA, 22–25 June. SPWLA-2003-DDD.
  23. Chang, C.T.P., Watson, A.T., and Edwards, C.M. 1999. NMR Imaging of Fluids and Flow in Porous Media. In Methods in the Physics of Porous Media, ed. P. Wong, No. 35, 387–423. San Diego, California: Experimental Methods in the Physical Sciences Series, Academic Press.
  24. Kleinberg, R.L., Farooqui, S.A., and Horsfield, M.A. 1993. T1/T2 Ratio and Frequency Dependence of NMR Relaxation in Porous Sedimentary Rocks. J. Colloid Interface Sci. 158 (1): 195-198. http://dx.doi.org/10.1006/jcis.1993.1247
  25. 25.0 25.1 Flaum, C., Bedford, J., and Kleinberg, R.L. 1998. Bound Water Volume, Permeability and Residual Oil Saturation From Incomplete Magnetic Resonance Logging Data. Presented at the SPWLA 39th Annual Logging Symposium, Keystone, Colorado, USA, 26–29 May. SPWLA-1998-UU.
  26. 26.0 26.1 Freedman, R., Lo, S., Flaum, M. et al. 2001. A New NMR Method of Fluid Characterization in Reservoir Rocks: Experimental Confirmation and Simulation Results. SPE J. 6 (4): 452-464. SPE-75325-PA. http://dx.doi.org/10.2118/75325-PA
  27. 27.0 27.1 Kenyon, W.E., Day, P.I., Straley, C. et al. 1988. A Three-Part Study of NMR Longitudinal Relaxation Properties of Water-Saturated Sandstones. SPE Form Eval 3 (3): 622–636. SPE-15643-PA. http://dx.doi.org/10.2118/15643-PA
  28. Kenyon, W.E., Ehrlich, R., Horkowitz, K. et al. 1989. Pore-Size Distribution and NMR in Microporous Cherty Sandstones. Presented at the SPWLA 30th Annual Logging Symposium, Denver, 11–14 June. SPWLA-1989-LL.
  29. Freedman, R. and Morriss, C.E. 1995. Processing of Data From an NMR Logging Tool. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, 22-25 October. SPE-30560-MS. http://dx.doi.org/10.2118/30560-MS
  30. 30.0 30.1 Miller, M.N., Paltiel, Z., Gillen, M.E. et al. 1990. Spin Echo Magnetic Resonance Logging: Porosity and Free Fluid Index Determination. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, 23-26 September. SPE-20561-MS. http://dx.doi.org/10.2118/20561-MS
  31. 31.0 31.1 31.2 31.3 31.4 31.5 31.6 Coates, G.R., Xiao, L.Z., and Prammer, M.G. 1999. NMR Logging: Principles and Applications, 234. Houston: Halliburton Energy Services.
  32. 32.0 32.1 32.2 Marschall, D., Gardner, J.S., Mardon, D. et al. 1995. Method for Correlating NMR Relaxometry and Mercury Injection Data. Presented at the Society of Core Analysts International Symposium, San Francisco, California, USA, 12–14 September. SCA-9511.
  33. Ausbrooks, R., Hurley, N.F., May, A. et al. 1999. Pore-Size Distributions in Vuggy Carbonates From Core Images, NMR, and Capillary Pressure. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 3-6 October. SPE-56506-MS. http://dx.doi.org/10.2118/56506-MS
  34. Pittman, E.D. 1992. Relationship of porosity and permeability to various parameters derived from mercury injection-capillary pressure curves for sandstone. AAPG Bulletin 76 (2): 191–198.
  35. Nelson, P. 1994. Permeability-porosity relationships in sedimentary rocks. The Log Analyst 35 (3): 38–62.
  36. Mirotchnik, K., Allsopp, K., and Kantzas, A. 1997. Combination of NMR and Ultracentrifuge Techniques for Effective Carbonate Reservoir Characterization. Presented at the Society of Core Analysts International Symposium, Calgary, 7–10 September. SCA-9703.
  37. 37.0 37.1 Prammer, M.G., Drack, E.D., Bouton, J.C. et al. 1996. Measurements of Clay-Bound Water and Total Porosity by Magnetic Resonance Logging. Presented at the SPE Annual Technical Conference and Exhibition, Denver, 6-9 October. SPE-36522-MS. http://dx.doi.org/10.2118/36522-MS
  38. Chitale, D.V., Day, P.I., and Coates, G.R. 1999. Petrophysical Implications of Laboratory NMR and Petrographical Investigation on a Shaly Sand Core. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 3-6 October. SPE-56765-MS. http://dx.doi.org/10.2118/56765-MS
  39. Brown, R.J.S. and Neuman, C.H. 1982. The Nuclear Magnetism Log--A Guide For Field Use. The Log Analyst 23 (5): 4. SPWLA-1982-vXXIIIn5a1.
  40. Cooper, R.K. and Jackson, J.A. 1980. Remote (Inside-Out) NMR. I. Remote Production of a Region of Homogeneous Magnetic Field. J. Magn. Reson. 41: 400–405.
  41. Burnett, L.J. and Jackson, J.A. 1980. Remote (Inside-Out) NMR. II. Sensitivity of NMR Detection for External Samples. J. Magn. Reson. 41: 406–410.
  42. ackson, J.A., Burnett, L.J., and Harmon, F. 1980. Remote (Inside-Out) NMR. III. Detection of Nuclear Magnetic Resonance in a Remotely Produced Region of Homogeneous Magnetic Field. J. Magn. Reson. 41: 411–421.
  43. Jackson, J.A. 1984. Nuclear Magnetic Resonance Well Logging. The Log Analyst 25 (5): 16. SPWLA-1984-vXXVn5a3.
  44. 44.0 44.1 Jackson, J.A. 2001. Los Alamos NMR well logging project. Concepts in Magnetic Resonance 13 (6): 368-378. http://dx.doi.org/10.1002/cmr.1022
  45. Chandler, R.N., Drack, E.O., Miller, M.N. et al. 1994. Improved Log Quality With a Dual-Frequency Pulsed NMR Tool. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, 25-28 September. SPE-28365-MS. http://dx.doi.org/10.2118/28365-MS
  46. 46.0 46.1 Prammer, M.G., Drack, E.D., Bouton, J.C. et al. 1996. Measurements of Clay-Bound Water and Total Porosity by Magnetic Resonance Logging. Presented at the SPE Annual Technical Conference and Exhibition, Denver, 6-9 October. SPE-36522-MS. http://dx.doi.org/10.2118/36522-MS
  47. Drack, E., Chandler, R., Bouton, I. et al. 1999. Theory and Operation of a New, Multi-Volume NMR Logging System. Presented at the SPWLA Annual Logging Symposium, Oslo, Norway, 30 May–3 June. SPWLA-1999-DD.
  48. Prammer, M.G., Bouton, J., Drack, E.D. et al. 2001. A New Multiband Generation of NMR Logging Tools. SPE Res Eval & Eng 4 (1): 59-63. SPE-69670-PA. http://dx.doi.org/10.2118/69670-PA
  49. Kleinberg, R.L., Sezginer, A., Griffin, D.D. et al. 1992. Novel NMR apparatus for investigating an external sample. J. Magn. Reson. 97 (3): 466-485. http://dx.doi.org/10.1016/0022-2364(92)90028-6
  50. Morriss, C.E., Kenyon, W.E., Tutunjian, P.N. et al. 1993. Field Test of an Experimental Pulsed Nuclear Magnetism Tool. Presented at the SPWLA Annual Logging Symposium, Calgary, 13–16 June. SPWLA-1993-GGG.
  51. Freedman, R., Morriss, C.E., Flaum, C. et al. 1997. Measurement of Total NMR Porosity Adds New Value to NMR Logging. Presented at the SPWLA Annual Logging Symposium, Houston, 15–18 June. SPWLA-1997-OO.
  52. Mirth, C.C., Willis, D., Hurlimann, M. et al. 1999. An Improved NMR Tool for Faster Logging. Presented at the SPWLA Annual Logging Symposium, Oslo, Norway, 30 May–3 June. SPWLA-1999-CC.
  53. Heaton, N.J., Freedman, R., Karmonik, C. et al. 2002. Applications of a New-Generation NMR Wireline Logging Tool. Presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, USA, 29 September-2 October. SPE-77400-MS. http://dx.doi.org/10.2118/77400-MS
  54. 54.0 54.1 DePavia, L., Heaton, N., Ayers, D. et al. 2003. A Next-Generation Wireline NMR Logging Tool. Presented at the SPE Annual Technical Conference and Exhibition, Denver, 5-8 October. SPE-84482-MS. http://dx.doi.org/10.2118/84482-MS
  55. 55.0 55.1 Minh, C.C., Bordon, E., Hurliman, M. et al. 2005. Field Test Results of the New Combinable Magnetic Resonance Autotune Logging Tool. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, 9-12 October. SPE-96759-MS. http://dx.doi.org/10.2118/96759-MS
  56. Chen, S., Fang, S., Zhang, G. et al. 2003. MR Explorer Log Acquisition Methods: Petrophysical-Objective-Oriented Approaches. Presented at the SPWLA Annual Logging Symposium, Galveston, Texas, USA, 22–25 June. SPWLA-2003-ZZ.
  57. Khamatdinov, R., Mityushin, E., Murtsovkin, V. et al. 2003. Field Test of a New Nuclear Magnetic Resonance Tool. Presented at the SPWLA Annual Logging Symposium, Galveston, Texas, USA, 22–25 June. SPWLA-2003-AAA.
  58. 58.0 58.1 Bonnie, R.J.M., Akkurt, R., Al-Waheed, H. et al. 2003. Wireline T1 Logging. Presented at the SPE Annual Technical Conference and Exhibition, Denver, 5-8 October. SPE-84483-MS. http://dx.doi.org/10.2118/84483-MS Cite error: Invalid <ref> tag; name "r58" defined multiple times with different content
  59. 59.0 59.1 59.2 Heaton, N.J., Minh, C.C., Kovats, J. et al. 2004. Saturation and Viscosity From Multidimensional Nuclear Magnetic Resonance Logging. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 26-29 September. SPE-90564-MS. http://dx.doi.org/10.2118/90564-MS
  60. Bonnie, R.J.M., Galford, J., Fox, P. et al. 2004. Advanced NMR Acquisition Has Improved the Assessment of Hydrocarbons in Place in the Gulf of Mexico. Presented at the SPWLA 45th Annual Logging Symposium, Noordwijk, The Netherlands, 6–9 June. SPWLA-2004-II.
  61. Howard, J., Reppert, M., Bonnie, R. et al. 2005. Porosity and Water Saturation from LWD NMR in a North Sea Chalk Formation. Presented at the SPWLA 46th Annual Logging Symposium, New Orleans, 26–29 June. Paper 2005-D.
  62. 62.0 62.1 Prammer, M., Menger, S., Knizhnik, S. et al. 2003. Directional Resonance: New Applications for MRIL. Presented at the SPE Annual Technical Conference and Exhibition, Denver, 5-8 October. SPE-84479-MS. http://dx.doi.org/10.2118/84479-MS
  63. Flaum, C., Speier, P., Kleinberg, R.L. et al. 1999. Reducing Motion Effects on Magnetic Resonance Bound Fluid Estimates. Presented at the SPWLA 40th Annual Logging Symposium, Oslo, Norway, 30 May-3 June. Paper II.
  64. Drack, E.D., Prammer, M.G., Zannoni, S. et al. 2001. Advances in LWD Nuclear Magnetic Resonance. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, 30 September-3 October. SPE-71730-MS. http://dx.doi.org/10.2118/71730-MS
  65. Prammer, M.G. 2001. NMR logging-while-drilling (1995–2000). Concepts in Magnetic Resonance 13 (6): 409-411. http://dx.doi.org/10.1002/cmr.1029
  66. Prammer, M.G., Menger, S., Akkurt, R. et al. 2002. A New Direction in Wireline and LWD NMR. Presented at the SPWLA 43rd Annual Logging Symposium, Oiso, Japan, 2–5 June. SPWLA-2002-DDD.
  67. Prammer, M.G., Drack, E., Goodman, G. et al. 2001. The Magnetic-Resonance While-Drilling Tool: Theory and Operation. SPE Res Eval & Eng 4 (4): 270-275. SPE-72495-PA. http://dx.doi.org/10.2118/72495-PA
  68. Appel, M., Radcliffe, N.J., Aadireddy, P. et al. 2002. Nuclear Magnetic Resonance While Drilling in the U.K. Southern North Sea. Presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, USA, 29 September-2 October. SPE-77395-PA. http://dx.doi.org/10.2118/77395-MS
  69. Morley, J., Heidler, R., Horkowitz, J. et al. 2002. Field Testing of a New Nuclear Magnetic Resonance Logging-While-Drilling Tool. Presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, USA, 29 September-2 October. SPE-77477-MS. http://dx.doi.org/10.2118/77477-MS
  70. Poitzsch, M., Scheibal, J.R., Hashem, M. et al. 2002. Applications of a New Magnetic Resonance Logging-While-Drilling Tool in a Gulf of Mexico Deepwater Development Project. Presented at the SPWLA 43rd Annual Logging Symposium, Oiso, Japan, 2–5 June. SPWLA-2002-EEE.
  71. Heidler, R., Morriss, C., and Hoshun, R. 2003. Design and Implementation of a New Magnetic Resonance Tool for the While Drilling Environment. Presented at the SPWLA 44th Annual Logging Symposium, Galveston, Texas, USA, 22–25 June. SPWLA-2003-BBB.
  72. Borghi, M., Porrera, F., Lyne, A. et al. 2005. Magnetic Resonance While Drilling Streamlines Reservoir Evaluation. Presented at the SPWLA 46th Annual Logging Symposium, New Orleans, 26–29 June. Paper 2005-HHH.
  73. Kleinberg, R.L. and Vinegar, H.J. 1996. NMR Properties of Reservoir Fluids. The Log Analyst 37 (6): 20–32.
  74. Prammer, M.G., Bouton, J., and Masak, P. 2001. The Downhole NMR Fluid Analyzer. Presented at the SPWLA Annual Logging Symposium, Houston, 16–20 June. SPWLA-2001-N.
  75. Bouton, J., Prammer, M.G., Masak, P. et al. 2001. Assessment of Sample Contamination by Downhole NMR Fluid Analysis. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, 30 September-3 October. SPE-71714-MS. http://dx.doi.org/10.2118/71714-MS.
  76. Masak, P.C., Prammer, M.G., Drack, E. et al. 2002. Field Test Results and Applications of the Downhole Magnetic Resonance Fluid Analyzer. Presented at the SPWLA 43rd Annual Logging Symposium, Oiso, Japan, 2–5 June. SPWLA-2002-GGG.
  77. Akkurt, R., Fransson, C.-M., Witkowsky, J.M. et al. 2004. Fluid Sampling and Interpretation with the Downhole NMR Fluid Analyzer. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 26-29 September. SPE-90971-MS. http://dx.doi.org/10.2118/90971-MS
  78. 78.0 78.1 78.2 Stambaugh, B., Svor, R., and Globe, M. 2000. Quality Control of NMR Logs. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, 1-4 October. SPE-63212-MS. http://dx.doi.org/10.2118/63212-MS
  79. 79.0 79.1 Marschall, D.M. 2000. HBVI: An NMR Method To Determine BVI as a Function of Reservoir Capillarity. Presented at the SPWLA 41st Annual Logging Symposium, Dallas, 4–7 June. SPWLA-2000-KK.
  80. 80.0 80.1 Timur, A. 1969. Effective Porosity and Permeability of Sandstones Investigated Through Nuclear Magnetic Principles. The Log Analyst 10 (1): 3.
  81. Chen, S., Ostroff, G., and Georgi, D.T. 1998. Improving Estimation of NMR Log T2 Cutoff Value With Core NMR and Capillary Pressure Measurements. Presented at the SCA International Symposium, The Hague, 14–16 September. SCA-9822.
  82. Agut, R., Levallois, B., and Klopf, W. 2000. Integrating Core Measurements and NMR Logs in Complex Lithology. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, 1-4 October. SPE-63211-MS. http://dx.doi.org/10.2118/63211-MS
  83. Zhang, Q., Lo, S.-W., Hirasaki, G.J. et al. 1998. Some Exceptions to Default NMR Rock and Fluid Properties. Presented at the SPWLA Annual Logging Symposium, Keystone, Colorado, USA, 26–29 May. SPWLA-1998-FF.
  84. 84.0 84.1 Coates, G.R., Marschall, D.M., Mardon, M. et al. 1998. A New Characterization of Bulk-Volume Irreducible Using Magnetic Resonance. The Log Analyst 39 (1): 51.
  85. 85.0 85.1 Chen, S., Liu, C., Georgi, D. et al. 1998. Methods for Computing SWI and BVI From NMR Logs. Presented at the SPWLA Annual Logging Symposium, Keystone, Colorado, USA, 26–29 May. SPWLA-1998-HH.
  86. Kleinberg, R.L. and Boyd, A. 1997. Tapered Cutoffs for Magnetic Resonance Bound Water Volume. Presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, USA, 5-8 October. SPE-38737-MS. http://dx.doi.org/10.2118/38737-MS
  87. Seevers, D.O. 1966. A Nuclear Magnetic Method for Determining the Permeability of Sandstones. Presented at the SPWLA 7th Annual Logging Symposium, Tulsa, 8–11 May. Paper L.
  88. Timur, A. 1969. Pulsed Nuclear Magnetic Resonance Studies of Porosity, Movable Fluid, and Permeability of Sandstones. J Pet Technol 21 (6): 775-786. SPE-2045-PA. http://dx.doi.org/10.2118/2045-PA
  89. Timur, A. 1968. An Investigation of Permeability, Porosity, and Residual Water Saturation Relationships for Sandstone Reservoirs. The Log Analyst 9 (4): 8.
  90. Ahmed, U., Crary, S.F., and Coates, G.R. 1991. Permeability Estimation: The Various Sources and Their Interrelationships. J Pet Technol 43 (5): 578-587. SPE-19604-PA. http://dx.doi.org/10.2118/19604-PA
  91. Worthington, P.F. 2003. The Effect of Scale on the Petrophysical Estimation of Intergranular Permeability. Presented at the SPWLA 44th Annual Logging Symposium, Galveston, Texas, USA, 22–25 June. SPWLA-2003-A.
  92. Bryant, S.L., Cade, C.A., and Melor, D.W. 1993. Permeability prediction from geological models. AAPG Bulletin 77 (8): 1338–1350.
  93. Chang, D., Vinegar, H., Morris, C. et al. 1997. Effective Porosity, Producible Fluid and Permeability in Carbonates From NMR Logging. The Log Analyst 38 (2): 60-72.
  94. Babadagli, T. and Al-Salmi, S. 2002. Improvement of Permeability Prediction for Carbonate Reservoirs Using Well Log Data. Presented at the SPE Asia Pacific Oil and Gas Conference and Exhibition, Melbourne, Australia, 8-10 October. SPE-77889-MS. http://dx.doi.org/10.2118/77889-MS
  95. Marschall, D., Gardner, J.S., and Curby, F.M. 1997. MR Laboratory Measurements—Requirements to Assure Successful Measurements that Will Enhance MRI Log Interpretation. Presented at the Society of Core Analysts International Symposium, Calgary, 8–10 September. SCA-9704.
  96. Prammer, M.G. 1994. NMR Pore Size Distribution and Permeability at the Well Site. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, 25–28 September. SPE-28368-MS. http://dx.doi.org/10.2118/28368-MS
  97. 97.0 97.1 97.2 Winkler, M., Freeman, J., and Appel, M. 2004. The Limits Of Fluid Property Correlations Used in NMR Well Logging; an Experimental Study of Reservoir Fluids at Reservoir Conditions. Presented at the SPWLA 45th Annual Logging Symposium, Noordwijk, The Netherlands, 6–9 June. SPWLA-2004-DD.
  98. 98.0 98.1 Hirasaki, G.J., Lo, S.-W., and Zhang, Y. 2003. NMR properties of petroleum reservoir fluids. Magn. Reson. Imaging 21 (3–4): 269-277. http://dx.doi.org/10.1016/S0730-725X(03)00135-8
  99. 99.0 99.1 Chen, S., Zhang, G., Kwak, H. et al. 2004. Laboratory Investigation of NMR Crude Oils and Mud Filtrates Properties in Ambient and Reservoir Conditions. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 26-29 September. SPE-90553-MS. http://dx.doi.org/10.2118/90553-MS
  100. Prammer, M.G. 2001. Hydrocarbon Saturation Measurements by NMR. Concepts in Magnetic Resonance 13 (6): 406.
  101. 101.0 101.1 Akkurt, R., Tutunjian, P.N., Guillory, A.J. et al. 1996. NMR Logging of Natural Gas Reservoirs. The Log Analyst 37 (6): 33. SPWLA-1996-v37n6a1.
  102. Coates, G.R., Miller, D.L., and Gardner, J.S. 1994. Applying Pulse-Echo NMR to Shaly Sand Formation Evaluation. Presented at the SPWLA Annual Logging Symposium, Tulsa, 19–22 June. SPWLA-1994-B.
  103. 103.0 103.1 103.2 Chitale, D.V., Day, P.I., and Coates, G.R. 1999. Petrophysical Implications of Laboratory NMR and Petrographical Investigation on a Shaly Sand Core. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 3-6 October. SPE-56765-MS. http://dx.doi.org/10.2118/56765-MS
  104. Moore, M.A. and Akkurt, R. 1996. Nuclear Magnetic Resonance Applied to Gas Detection in a Highly Laminated Gulf of Mexico Turbidite Invaded with Synthetic Oil Filtrate. Presented at the SPE Annual Technical Conference and Exhibition, Denver, 6-9 October. SPE-36521-MS. http://dx.doi.org/10.2118/36521-MS
  105. 105.0 105.1 Buller, D. 2000. Carbonate Evaluation Using NMR Time Domain Analysis. Presented at the SPWLA 41st Annual Logging Symposium, Dallas, 4–7 June. SPWLA-2000-GGG.
  106. Prammer, M.G., Mardon, D., Coates, G.R. et al. 1995. Lithology-Independent Gas Detection by Gradient-NMR Logging. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, 22-25 October. SPE-30562-MS. http://dx.doi.org/10.2118/30562-MS
  107. Hou, B.L. and Miller, D. 2000. Determining Fluid Volume in Gas and/or Light Oil Reservoirs—Using a New Triple-Wait-Time NMR Logging Method. Presented at the SPWLA 41st Annual Logging Symposium, Dallas, 4–7 June. SPWLA-2000-U.
  108. 108.0 108.1 Timur, A. 1972. Nuclear Magnetic Resonance Study of Carbonate Rocks. The Log Analyst 13 (5): 3.
  109. Akkurt, R., Mardon, D., Solanet, F. et al. 1998. Enhanced Diffusion—Expanding the Range of NMR Direct Hydrocarbon-Typing Applications. Presented at the SPWLA 39th Annual Logging Symposium, Keystone, Colorado, USA, 26–29 May. SPWLA-1998-GG.
  110. 110.0 110.1 110.2 Hürlimann, M.D. 2001. Diffusion and Relaxation Effects in General Stray Field NMR Experiments. J. Magn. Reson. 148 (2): 367-378. http://dx.doi.org/10.1006/jmre.2000.2263
  111. 111.0 111.1 Akkurt, R., Marschall, D., Eyvazzadeh, R.Y. et al. 1999. Determination of Residual Oil Saturation by Use of Enhanced Diffusion. SPE Res Eval & Eng 2 (3): 303-309. SPE-56990-PA. http://dx.doi.org/10.2118/56990-PA
  112. 112.0 112.1 Clerke, E.E., Hartman, D.E., Vinegar, H.J. et al. 1997. Residual Oil Saturation Measurements In Carbonates With Pulsed Nmr Logs. The Log Analyst 38 (2): 73. SPWLA-1997-v38n2a2.
  113. 113.0 113.1 113.2 Howard, J.J. 1994. Wettability and Fluid Saturations Determined From NMR T1 Distributions. Magn. Reson. Imaging 12 (2): 197-200. http://dx.doi.org/10.1016/0730-725X(94)91514-8
  114. Flaum, C., Kleinberg, R.L., and Htirlimann, M.D. 1996. Identification of Gas With the Combinable Magnetic Resonance Tool (CMR). Presented at the SPWLA Annual Logging Symposium, New Orleans, 16–19 June. SPWLA-1996-L.
  115. 115.0 115.1 Sun, B., Van Dalen, S.C., Al-Rushaid, A. et al. 2004. Two-dimensional NMR Logging and Field Test Results. Presented at the SPWLA 45th Annual Logging Symposium, Noordwijk, The Netherlands, 6–9 June. SPWLA-2004-KK.
  116. Slijkerman, W.F.J., Looyestijn, W.J., Hofstra, P. et al. 2000. Processing of Multi-Acquisition NMR Data. SPE Res Eval & Eng 3 (6): 492-497. SPE-68408-PA. http://dx.doi.org/10.2118/68408-PA
  117. Chen, S., Hursan, G., Beard, D. et al. 2003. Correction for Processing Multigradient, Multiple- NMR Log Data. Presented at the SPE Annual Technical Conference and Exhibition, Denver, 5-8 October. SPE-84481-MS. http://dx.doi.org/10.2118/84481-MS
  118. 118.0 118.1 118.2 118.3 118.4 Freedman, R., Heaton, N., Flaum, M. et al. 2002. Wettability, Saturation, and Viscosity Using the Magnetic Resonance Fluid Characterization Method and New Diffusion-Editing Pulse Sequences. Presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, USA, 29 September-2 October. SPE-77397-MS. http://dx.doi.org/10.2118/77397-MS
  119. Fang, S., Chen, S., Tauk, R. et al. 2004. Quantification of Hydrocarbon Saturation in Carbonate Formations Using Simultaneous Inversion of Multiple NMR Echo Trains. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 26-29 September. SPE-90569-MS. http://dx.doi.org/10.2118/90569-MS
  120. 120.0 120.1 Sun, B. and Dunn, K.-J. 2005. A global inversion method for multi-dimensional NMR logging. J. Magn. Reson. 172 (1): 152-160. http://dx.doi.org/10.1016/j.jmr.2004.10.003
  121. Sun, B. and Dunn, K.-J. 2005. Two-dimensional nuclear magnetic resonance petrophysics. Magn. Reson. Imaging 23 (2): 259-262. http://dx.doi.org/10.1016/j.mri.2004.11.056
  122. 122.0 122.1 Bonnie, R.J.M., Marschall, D.M., Fam, M.Y. et al. 2001. Advanced Forward Modeling Helps Planning and Interpreting NMR Logs. Presented at the SPE Annual Technical Conference and Exhibition, New Orleans, 30 September-3 October. SPE-71735-MS. http://dx.doi.org/10.2118/71735-MS
  123. Hürlimann, M.D. and Venkataramanan, L. 2002. Quantitative Measurement of Two-Dimensional Distribution Functions of Diffusion and Relaxation in Grossly Inhomogeneous Fields. J. Magn. Reson. 157 (1): 31-42. http://dx.doi.org/10.1006/jmre.2002.2567
  124. Hürlimann, M.D., Venkataramanan, L., Flaum, C. et al. 2002. Diffusion Editing: New NMR Measurement of Saturation and Pore Geometry. Presented at the SPWLA 43rd Annual Logging Symposium, Oiso, Japan, 2–5 June. SPWLA-2002-FFF.
  125. 125.0 125.1 Sun, B. and Dunn, K. 2002. Core Analysis With Two Dimensional NMR. Presented at the 2002 International Symposium of the Society of Core Analysts, Monterey, California, USA, 22-25 September. SCA 2002-38.
  126. Hürlimann, M.D., Flaum, M., Venkataramanan, L. et al. 2003. Diffusion-relaxation distribution functions of sedimentary rocks in different saturation states. Magn. Reson. Imaging 21 (3–4): 305-310.
  127. 127.0 127.1 Toumelin, E., Torres-Verdin, C., Sun, B. et al. 2004. A Numerical Assessment of Modern Borehole NMR Interpretation Techniques. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 26-29 September. SPE-90539-MS. http://dx.doi.org/10.2118/90539-MS
  128. Hursan, G., Chen, S., and Murphy, E. 2005. New NMR Two-Dimensional Inversion of T1/T2app vs. T2app Method for Gas Well Petrophysical Interpretation. Presented at the SPWLA 46th Annual Logging Symposium, New Orleans, 26–29 June. Paper GGG.
  129. 129.0 129.1 Flaum, M., Chen, C., and Hirasaki, G. 2004. NMR Diffusion Editing for D-T2 Maps—Application to Recognition of Wettability Change. Presented at the SPWLA 45th Annual Logging Symposium, Noordwijk, The Netherlands, 6–9 June. SPWLA-2004-JJ.
  130. Sun, B. et al. 2003. NMR Imaging with Diffusion and Relaxation. Presented at the International Symposium of the Society of Core Analysts, Pau, France, 21–24 September. SCA 2003-24.
  131. 131.0 131.1 Cory, D.G., Sen, P.N., Guzmn-Garcia, A.G. et al. 2002. Nmr Identification Of Fluids And Wettability In Situ In Preserved Cores. Petrophysics 43 (1): 13. SPWLA-2002-v43n1a1.
  132. Klimentos, T. 2003. NMR Applications in Petroleum Related Rock-Mechanics: Sand Control, Hydraulic Fracturing, Wellbore Stability. Presented at the SPWLA 44th Annual Logging Symposium, Galveston, Texas, USA, 22–25 June. SPWLA-2003-HHH.
  133. van der Zwaag, C.H., Veliyulin, E., Skjetne, T. et al. 2003. Deformation and Failure of Rock Samples Probed by T1 and T2 Relaxation. Magn. Reson. Imaging 21 (3–4): 405-407. http://dx.doi.org/10.1016/S0730-725X(03)00151-6
  134. Mullen, M. et al. : Mullen, M., Bray, J., and Bonnie, R. 2005. Fluid Typing With T1 NMR: Incorporating T1 and T2 Measurements for Improved Interpretation in Tight Gas Sands and Unconventional Reservoirs. Presented at the SPWLA Annual Logging Symposium, New Orleans, 26–29 June. SPWLA-2005-III.
  135. Longis, C., Vignau, S., and White, J.H. 2005. NMR and Capture Spectroscopy Help Resolve Producibility and Fluid Distribution in the North Alwyn Triassic. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, 9-12 October. SPE-96609-MS. http://dx.doi.org/10.2118/96609-MS
  136. 136.0 136.1 Latorraca, G.A., Carison, R.M., Stonard, S.W. et al. 1999. Heavy Oil Viscosity Determination Using NMR Logs. Presented at the SPWLA Annual Logging Symposium, Oslo, Norway, 30 May–3 June. SPWLA-1999-PPP.
  137. 137.0 137.1 Galford, J.E. and Marschall, D.M. 2000. Combining NMR and Conventional Logs to Determine Fluid Volumes and Oil Viscosity in Heavy-Oil Reservoirs. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, 1-4 October. SPE-63257-MS. http://dx.doi.org/10.2118/63257-MS
  138. Seccombe, J., Bonnie, R.J.M., Smith, M. et al. 2005. Ranking Oil Viscosity in Heavy-Oil Reservoirs. Presented at the SPE/PS-CIM/CHOA International Thermal Operations and Heavy Oil Symposium, Calgary, 1-3 November. SPE-97935-MS. http://dx.doi.org/10.2118/97935-MS
  139. Nascimento, J.d.D.S. and Gomes, R.M.R. 2004. Tar Mats Characterization from NMR and Conventional Logs, Case Studies in Deepwater Reservoirs, Offshore Brazil. Presented at the SPWLA 45th Annual Logging Symposium, Noordwijk, The Netherlands, 6–9 June. SPWLA-2004-FF.
  140. Bustos, U.D., Ortiz, A.C., Breda, E.W. et al. 2004. Detecting Restricted Diffusion Effect when Identifying Hydrocarbons with NMR Logs. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 26-29 September. SPE-90089-MS. http://dx.doi.org/10.2118/90089-MS
  141. 141.0 141.1 Westphal, H., Surholt, I., Kiesl, C. et al. 2005. NMR Measurements in Carbonate Rocks: Problems and an Approach to a Solution. Pure Appl. Geophys. 162 (3): 549-570. http://dx.doi.org/10.1007/s00024-004-2621-3
  142. Rose, D., Hansen, P.M., Damgaard, A.P. et al. 2003. A Novel Approach to Real Time Detection of Facies Changes in Horizontal Carbonate Wells Using LWD NMR. Presented at the SPWLA 44th Annual Logging Symposium, Galveston, Texas, USA, 22–25 June. SPWLA-2003-CCC.
  143. Chen, S., Mette, M., Hursan, G. et al. 2005. Simple, Robust NMR-Based Indicators for Detection of Hydrocarbon Gas or Oil, Borehole Contamination, and Vugs. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, 9-12 October. SPE-96409-MS. http://dx.doi.org/10.2118/96409-MS
  144. Fleury, M., Al-Nayadi, K., and Boyd, D. 2005. Water Saturation from NMR, Resistivity and Oil Base Core in a Heterogeneous Middle East Carbonate Reservoir. Presented at the SPWLA 46th Annual Logging Symposium, New Orleans, 26–29 June. Paper JJJ.
  145. Volokitin, Y., Hofman, J., Slijkermail, W. et al. 2001. Constructing Capillary Pressure Curves from NMR Log Data in the Presence of Hydrocarbons. Petrophysics 42 (4): 334.
  146. Glorioso, J.C., Aguirre, O., Piotti, G. et al. 2003. Deriving Capillary Pressure and Water Saturation from NMR Transversal Relaxation Times. Presented at the SPE Latin American and Caribbean Petroleum Engineering Conference, Port-of-Spain, Trinidad and Tobago, 27-30 April. SPE-81057-MS. http://dx.doi.org/10.2118/81057-MS
  147. Grattoni, C.A. et al. 2003. An Improved Technique for Deriving Drainage Capillary Pressure from NMR T2 Distributions. Presented at the International Symposium of the Society of Core Analysts, Pau, France, 21–24 September. SCA 2003-25.
  148. Dastidar, R., Rai, C., and Sondergeld, C. 2004. Integrating NMR With Other Petrophysical Information to Characterize a Turbidite Reservoir. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 26-29 September. SPE-89948-MS. http://dx.doi.org/10.2118/89948-MS
  149. Corbelleri, A., Ortea, J., Lenge, D. et al. 1996. Application of the Magnetic Resonance Logging in San Jorge Basin (Argentina). Presented at the SPWLA Annual Logging Symposium, New Orleans, 16–19 June. SPWLA-1996-VV.
  150. Davis, B., McDonald, T., and Aly, M. 2002. Practical Applications of NMR Technology Enhance Formation Evaluation, Testing and Completion Decisions. Presented at the SPE Western Regional/AAPG Pacific Section Joint Meeting, Anchorage, Alaska, 20-22 May. SPE-76717-MS. http://dx.doi.org/10.2118/76717-MS
  151. Soliman, M.Y., Creel, P., Rester, S. et al. 2000. Technology Integration of Conformance and Magnetic Resonance Imaging Enhances Stimulation Treatment. Presented at the SPE Permian Basin Oil and Gas Recovery Conference, Midland, Texas, USA, 21-23 March. SPE-59530-MS. http://dx.doi.org/10.2118/59530-MS.
  152. Weiland, J., Teff, C., and Donovan, G. 2004. Integrating NMR, Rock Mechanics, and Production Evaluation to Optimize Hydrocarbon Production in Deepwater GOM. Presented at the SPWLA 45th Annual Logging Symposium, Noordwijk, The Netherlands, 6–9 June. SPWLA-2004-VVV.
  153. Hodgkins, M.A. and Howard, J.J. 1999. Application of NMR logging to reservoir characterization of low-resistivity sands in the Gulf of Mexico. AAPG Bull. 83 (1): 114-127.
  154. Ostroff, G.M. and Shorey, D.S. 1999. Integration of NMR and Conventional Log Data For Improved Petrophysical Evaluation of Shaly Sands. Presented at the SPWLA 40th Annual Logging Symposium, Oslo, Norway, 30 May-3 June. Paper OOO.
  155. Freeman, J.J., Freedman, R., Cao Minh, C. et al. 1998. Combining NMR and Density Logs for Petrophysical Analysis in Gas-Bearing Formations. Presented at the SPWLA 39th Annual Logging Symposium, Denver, 26–29 May. Paper II.
  156. Minh, C.C., Gubelin, G., Ramamoorthy, R. et al. 1999. Sonic-Magnetic Resonance Method: A Sourceless Porosity Evaluation in Gas-Bearing Reservoirs. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 3-6 October. SPE-56767-MS. http://dx.doi.org/10.2118/56767-MS
  157. 157.0 157.1 Chen, S., Georgi, D., Fang, S. et al. 1999. Optimization of NMR Logging Acquisition and Processing. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 3-6 October. SPE-56766-MS. http://dx.doi.org/10.2118/56766-MS
  158. Akkurt, R., Prammer, M.G., and Moore, M.A. 1996. Selection Of Optimal Acquisition Parameters For MRIL Logs. The Log Analyst 37 (6): 43. SPWLA-1996-v37n6a2.
  159. oumelin, E., Torres-Verdin, C., and Chen, S. 2002. Quantification of Multi-Phase Fluid Saturations in Complex Pore Geometries From Simulations of Nuclear Magnetic Resonance Measurements. Presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, USA, 29 September-2 October. SPE-77399-MS. http://dx.doi.org/10.2118/77399-MS
  160. Appel, M., Freeman, J.J., Perkins, R.B. et al. 1999. Fluid Identification in Heavy Oil Reservoirs by Nuclear Magnetic Resonance at Elevated Temperature. Presented at the International Symposium of the Society of Core Analysts, Golden, Colorado, USA, 1–4 August. SCA-9917.
  161. Appel, M., Freeman, J.J., Perkins, R.B. et al. 1999. Fluid Identification in Heavy Oil Reservoirs by Nuclear Magnetic Resonance at Elevated Temperature. Presented at the International Symposium of the Society of Core Analysts, Golden, Colorado, USA, 1–4 August. SCA-9917.
  162. Wyatt, D.F. 1999. Nuclear Magnetic Resonance Logging Aids in Defining Heavy Oil Types. Presented at the International Thermal Operations and Heavy Oil Symposium, Bakersfield, California, USA, 17–19 March. SPE-54121-MS. http://dx.doi.org/10.2118/54121-MS
  163. Carmona, R. and Decoster, S.A.E. 2001. Assessing Production Potential of Heavy Oil Reservoirs From the Orinoco Belt with NMR Logs. Presented at the SPWLA Annual Logging Symposium, Houston, 17–20 June. SPWLA-2001-ZZ.
  164. Padilha, S. and Baudot, P. 2001. Application of NMR Logging for Heavy Oil Identification and Quantification in OBM Environment. Presented at the SPWLA 42nd Annual Logging Symposium, Houston, 17–20 June. SPWLA-2001-FF.
  165. Manalo, F.P., Bryan, J.L., and Kantzas, A. 2002. Addressing the Clay/Heavy Oil Interaction When Interpreting Low Field Nuclear Magnetic Resonance Logs. Presented at the SPE International Thermal Operations and Heavy Oil Symposium and International Horizontal Well Technology Conference, Calgary, 4-7 November. SPE-78971-MS. http://dx.doi.org/10.2118/78971-MS
  166. Najia, W., Shaban, M., Gossa, H. et al. 2002. Nuclear Magnetic Resonance (NMR), a valuable tool for Tar detection in a Carbonate Formation of Abu Dhabi. Presented at the Abu Dhabi International Petroleum Exhibition and Conference, Abu Dhabi, UAE, 13-16 October. SPE-78485-MS. http://dx.doi.org/10.2118/78485-MS
  167. Manalo, F., Ding, M., Bryan, J. et al. 2003. Separating the Signals from Clay Bound Water and Heavy Oil in NMR Spectra of Unconsolidated Samples. Presented at the SPE Annual Technical Conference and Exhibition, Denver, 5-8 October. SPE-84480-MS. http://dx.doi.org/10.2118/84480-MS
  168. Hirasaki, G.J., Zhang, G.Q., and Huang, C.-C. 2000. Interpretation Of Wettability In Sandstones With Nmr Analysis. Petrophysics 41 (3): 223. SPWLA-2000-v41n3a1.
  169. Guan, H., Brougham, D., Sorbie, K.S. et al. 2002. Wettability effects in a sandstone reservoir and outcrop cores from NMR relaxation time distributions. J. Pet. Sci. Eng. 34 (1–4): 35-54. http://dx.doi.org/10.1016/S0920-4105(02)00151-1
  170. Eidesmo, T., Relling, O.M., and Rueslâtten, H. 1996. NMR Logging of Mixed-Wet North Sea Sandstone Reservoirs. Presented at the SPWLA Annual Logging Symposium, New Orleans, 16–19 June. SPWLA-1996-PP.
  171. Akkurt, R., Merkel, R.H., Coates, G.R. et al. 2000. Remaining Oil Saturation from NMR in a Mixed-Wet, Three-Phase Carbonate Reservoir. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, 1-4 October. SPE-63216-MS. http://dx.doi.org/10.2118/63216-MS
  172. Al-Mahrooqi, S.H., Grattoni, C.A., Moss, A.K. et al. 2003. An investigation of the effect of wettability on NMR characteristics of sandstone rock and fluid systems. J. Pet. Sci. Eng. 39 (3–4): 389-398. http://dx.doi.org/10.1016/S0920-4105(03)00077-9
  173. Fleury, M., Al-Nayadi, K., and Boyd, D. 2005. Water Saturation from NMR, Resistivity and Oil Base Core in a Heterogeneous Middle East Carbonate Reservoir. Presented at the SPWLA 46th Annual Logging Symposium, New Orleans, 26–29 June. Paper JJJ.
  174. Looyestijn, W.J. and Hofman, J. 2005. Wettability Index Determination by Nuclear Magnetic Resonance. Presented at the SPE Middle East Oil and Gas Show and Conference, Kingdom of Bahrain, 12–15 March. SPE-93624-MS. http://dx.doi.org/10.2118/93624-MS
  175. Howard, J.J. 1998. Quantitative estimates of porous media wettability from proton NMR measurements. Magn. Reson. Imaging 16 (5): 529-533.
  176. Chen, J., Hirasaki, G.J., and Flaum, M. 2004. Effects of OBM Invasion on Irreducible Water Saturation: Mechanisms and Modifications of NMR Interpretation. Presented at the SPE Annual Technical Conference and Exhibition, Houston, 26-29 September. SPE-90141-MS. http://dx.doi.org/10.2118/90141-MS
  177. Shafer, J., Flaum, M., Boyd, A. et al. 2004. Core And Log NMR Measurements Indicate Reservoir Rock Is Altered By OBM Filtrate. Presented at the SPWLA Annual Logging Symposium, Noordwijk, The Netherlands, 6–9 June. SPWLA-2004-CC.
  178. Rodrigues, P., Alcocer, Y., and Matheus, J. 2000. NMR Signal Pattern Classification for Estimation of Petrophysical Properties. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, 1-4 October. SPE-63258-MS. http://dx.doi.org/10.2118/63258-MS
  179. Sigal, R.F., Miller, D.L., Galford, J.E. et al. 2000. A Method for Enhancing the Vertical Resolution of NMR Logs. Presented at the SPE Annual Technical Conference and Exhibition, Dallas, 1-4 October. SPE-63215-MS. http://dx.doi.org/10.2118/63215-MS
  180. Akkurt, R. and Cherry, R. 2001. The Key to Improving the Vertical Resolution of Multi-Frequency Logging Tools. Presented at the SPWLA 42nd Annual Logging Symposium, Houston, 17–20 June. SPWLA-2001-YY.

SI Metric Conversion Factors


cm2/sec × 1.0* E + 02 = mm2/sec
cp × 1.0* E − 03 = Pa•s
°C 273.15 = K
dynes/cm2 × 1.0* E − 01 = Pa
ft × 3.048* E − 01 = m
ft/min × 5.080* E − 03 = m/sec
°F (°F + 459.67)/1.8 = K
gauss × 1.0* E − 04 = Tesla
in. × 2.54* E − 02 = m
ohm-m × 1.0* E − 02 = ohm-cm
psi × 6.894 757 E + 00 = kPa


*

Conversion factor is exact.