US20060158184A1 - Multiple echo train inversion - Google Patents

Multiple echo train inversion Download PDF

Info

Publication number
US20060158184A1
US20060158184A1 US11/037,834 US3783405A US2006158184A1 US 20060158184 A1 US20060158184 A1 US 20060158184A1 US 3783405 A US3783405 A US 3783405A US 2006158184 A1 US2006158184 A1 US 2006158184A1
Authority
US
United States
Prior art keywords
echo
echo trains
bins
trains
earth formation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Abandoned
Application number
US11/037,834
Inventor
Carl Edwards
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Baker Hughes Holdings LLC
Original Assignee
Baker Hughes Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Baker Hughes Inc filed Critical Baker Hughes Inc
Priority to US11/037,834 priority Critical patent/US20060158184A1/en
Assigned to BAKER HUGHES INCORPORATED reassignment BAKER HUGHES INCORPORATED ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: EDWARDS, CARL M.
Priority to GB0526441A priority patent/GB2422198B/en
Publication of US20060158184A1 publication Critical patent/US20060158184A1/en
Priority to US11/781,522 priority patent/US7812602B2/en
Abandoned legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/18Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
    • G01V3/32Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with electron or nuclear magnetic resonance
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/081Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity

Definitions

  • This invention is related to methods for acquiring nuclear magnetic resonance (NMR) measurements for determination of petrophysical properties of formations and properties of fluids therein. Specifically, the invention deals with simultaneous inversion of multiple echo trains for processing and interpreting Nuclear Magnetic Resonance (NMR) log data that exhibit relaxation and/or polarization.
  • NMR Nuclear Magnetic Resonance
  • Nuclear magnetic resonance is used in the oil industry, among others, and particularly in certain oil well logging tools.
  • NMR instruments may be used for determining, among other things, the fractional volume of pore space and the fractional volume of mobile fluid filling the pore space of earth formations. Methods of using NMR measurements for determining the fractional volume of pore space and the fractional volume of mobile fluids are described, for example, in “Spin Echo Magnetic Resonance Logging: Porosity and Free Fluid Index Determination,” M. N. Miller et al., Society of Petroleum Engineers paper no. 20561, Richardson, Tex., 1990. Further description is provided in U.S. Pat. No. 5,585,720, of Carl M. Edwards, issued Dec. 17, 1996 and having the same assignee as the present application, entitled “Signal Processing Method For Multiexponentially Decaying Signals And Applications To Nuclear Magnetic Resonance Well Logging Tools.” The disclosure of that patent is incorporated herein by reference.
  • the total porosity includes clay-bound water (CBW), capillary bound water (also known as Bulk volume Irreducible or BVI), movable water and hydrocarbons.
  • CBW clay-bound water
  • BVI capillary bound water
  • movable water movable water
  • hydrocarbons hydrocarbons.
  • the effective porosity a quantity of interest to production engineers, is the sum of the last three components and does not include the clay bound water.
  • Accurate spectra are also essential to estimate the irreducible and the movable fluid volumes; distortion of partial porosity distributions that has been commonly observed for a variety of reasons, affects the estimates of these quantities. The reasons for the distortions to occur are mainly due to poor signal-to-noise ratio (SNR) and poor resolution in the time domain of the NMR data.
  • SNR signal-to-noise ratio
  • the most common NMR log acquisition and core measurement method employs T 2 measurements using CPMG (Carr, Purcell, Meiboom and Gill) sequence, as taught by Meiboom and Gill in “Modified Spin-Echo Method for Measuring Nuclear Relaxation Time,” Rev. Sci. Instrum. 1958, 29, pp. 688-691.
  • the echo data in any given echo train are collected at a fixed time interval, the interecho time (TE).
  • TE interecho time
  • TE interecho time
  • a few hundred to a few thousand echoes are acquired to sample relaxation decay. Determining a light oil component, which has long relaxation time, requires taking several hundreds of ms of data while determination of CBW, which decays very fast, can be done with echo sequences of as short as a few tens of milliseconds.
  • NMR logging techniques used for obtaining information about earth formations and fluids.
  • MWD measurement-while-drilling
  • measurements are made while the wellbore is being drilled while in wireline logging, measurements are made after a wellbore has been drilled.
  • the logging tools are lowered into the borehole and NMR signals are obtained using different configurations of magnets, transmitter coils and receiver coils.
  • a static magnetic field is produced in the formation using permanent or electro-magnets. The static field aligns nuclear spins within the formation parallel to the static field.
  • a pulsed RF field is applied using a transmitter on the logging tool and the nuclear magnetization signals produced by the pulsed RF field are analyzed to determine formation properties.
  • the prior art shows different radio frequency (RF) pulsing schemes for generating RF fields in the formation.
  • the most commonly used pulsing schemes are variations of the CPMG sequence denoted by ( TW i ,90 ⁇ /2 ,( ⁇ ,180 ⁇ ,echo) j ) i (1)
  • TW is a wait time
  • 90 is a tipping pulse that tips the nuclear spins by an angle substantially equal to 90°
  • 180 is a refocusing pulse that tips the nuclear spins by an angle substantially equal to 180°
  • echo is a spin echo.
  • the time interval between successive refocusing pulses is 2 ⁇
  • the number of echoes is j
  • i denotes repetitions of the basic pulse sequence.
  • a variation of the CPMG sequence is taught in U.S. Pat. No. 6,163,153 to Reiderman in which the use of a refocusing pulse with a tipping angle less than 180° is disclosed.
  • Rig time is expensive, so that the general objective in wireline logging is to obtain interpretable data within as short a time as possible.
  • MWD logging on the other hand, no additional rig time is involved.
  • the parameters that may be varied are the acquisition frequencies and the number of different frequencies, the tip angles, the wait time, the number of pulses within a CPMG sequence, and the time interval between the pulses.
  • Long wait times are needed for proper evaluation of formation fluids that have long relaxation times, e.g., gas reservoirs while short wait times and/or short pulse spacings are used for evaluating faster relaxing components, e.g., irreducible fluid (BVI) and clay bound water (CBW).
  • BVI irreducible fluid
  • CBW clay bound water
  • U.S. Pat. No. 6,331,775 to Thern et al. having the same assignee as the present application discusses the use of a dual wait time acquisition for determination of gas saturation in a formation.
  • U.S. Pat. No. 5,023,551 to Kleinberg et al discusses the use of CPMG sequences in well logging.
  • U.S. Pat. No. 6,069,477 to Chen et al teaches the use of pulse sequences with different pulse spacings to determine CBW.
  • NMR fluid typing applications often involve acquiring multiple echo trains to exploit the diffusion and/or polarization contrasts between the water and hydrocarbon phases. Although proven successful in wells where these contrasts are large, the NMR-based techniques are challenging when the contrasts are small. Carbonates generally exhibit smaller NMR surface relaxivity than clastics, which reduces the relaxation time contrast between movable water and light oil. This difficulty is exacerbated by the small difference in the diffusivities of water and very light oil. In such cases, the methods for processing data are critical and the introduction of a priori information important.
  • New generation NMR well logging tools can acquire multiple echo trains with different TWs and TEs in a single logging pass. Simultaneously inverting all echo data to obtain the different fluid T 2 spectra is a logical method for the processing and interpretation of multiple echo trains acquired with different acquisition parameters. In fact, echo trains acquired with different TEs and TWs can not simply be averaged together to increase the signal-to-noise and inverted because they may exhibit very different T 2 decay behaviors. Thus, it is desirable to use an analysis technique that accounts for the different acquisition parameters in the individual echo trains.
  • Multiecho sequences are acquired from a first and second region of interest using a first and second radio frequency (RF) pulse sequence.
  • RF radio frequency
  • a correction factor depending at least in part on a diffusivity of a fluid in the earth formation is determined, and the first and second multiecho sequences are combined using the correction factor to obtain a combined multiecho sequence.
  • the method taught therein is difficult to extend to multiple echo trains, and the specific problem of small contrasts is not addressed.
  • a similar approach (combining echo trains without gradient variations) is discussed in U.S. Pat. No. 6,377,042 issued to Menger.
  • EP 0886792 to Bonnie et al. discloses a method in which NMR echo signals are acquired with different combinations of TW, TE and field gradient (produced by a gradient coil), and a curve fitting procedure is used to match the resulting signals to a predetermined model.
  • the model itself is limited in scope.
  • One embodiment of the present invention is a method of characterizing an earth formation having a fluid therein.
  • the method includes obtaining a plurality of nuclear magnetic resonance (NMR) echo trains, each of the echo trains resulting from pulsing of the earth formation with an associated pulse sequence having an associated wait time.
  • a plurality of bins of a T 2 distribution of the earth formation is defined.
  • an associated cutoff time for full polarization is determined, and a value associated with each of the plurality of bins is determined based at least in part on the plurality of echo trains and the associated cutoff time.
  • At least one of NMR echo trains may be a fully polarized echo train.
  • the echo trains may be fully polarized echo trains, a clay bound water echo train, or a bulk volume irreducible echo train.
  • Each of the associated cutoff times may be determined from a wait time of the associated pulse sequence and a maximum longitudinal relaxation time of the earth formation.
  • a least squares inversion may be performed.
  • the inversion may be a weighted inversion. Determination of the value associated with each of the plurality of bins involves a simultaneous fitting over each of the echo trains for a bin that is fully polarized for each of the plurality of echo trains and may further involve a simultaneous fitting of a subset of the echo trains for a bin that is fully polarized for the subset of echo trains.
  • At least one characteristic of the earth formation selected from (i) clay bound water, (ii) bulk volume irreducible, and, (iii) porosity may be determined.
  • the apparatus includes a nuclear magnetic resonance (NMR) tool conveyed in a borehole in the earth formation.
  • the NMR tool pulses the earth formation with a plurality of radio frequency (RF) magnetic pulse sequences and receives a plurality of associated echo trains.
  • RF radio frequency
  • a processor determines from the plurality of echo trains a value associated with each of a plurality of bins of a T 2 distribution, the determination being based at least in part on cutoff times for full polarization associated with each of the plurality of echo trains.
  • the echo trains may be a fully polarized echo trains, a clay bound water echo trains or a BVI echo train.
  • the associated cutoff times may be determined from a wait time of the associated pulse sequence and a maximum longitudinal relaxation time of the earth formation.
  • the processor may determine the value associated with each of the plurality of bins by further performing a least squares inversion.
  • the inversion may be a weighted inversion.
  • the inversion may involve a simultaneous fitting over each of the echo trains for a bin that is fully polarized for each of the plurality of echo trains and may further involve a simultaneous fitting of a subset of the echo trains for a bin that is fully polarized for the subset of echo trains.
  • the processor may be at a surface location, a downhole location, or a remote location.
  • the NMR tool may be conveyed in the borehole on a wireline, a drillstring, or coiled tubing.
  • the processor may determine a characteristic of the earth formation such as clay bound water, bulk volume irreducible, and porosity.
  • the medium includes instructions that enable a nuclear magnetic resonance (NMR) tool conveyed in a borehole to acquire a plurality of echo trains.
  • the echo trains result from pulsing of the earth formation by an associated pulse sequence that has a wait time.
  • the instructions also enable determination from the plurality of echo trains a value associated with each of a plurality of bins of a T 2 distribution. The determination is based at least in part on cutoff times for full polarization associated with each of the plurality of echo trains.
  • the instructions enable the processor to process a fully polarized echo train, a CBW echo train, and a BVI echo train.
  • the instructions further enable performing performing a least squares inversion.
  • the machine readable medium may be a ROM, an EPROM, an EAROMs, a Flash Memory, or an optical disk.
  • the instructions may further enable determination of a characteristic of the formation such as clay bound water, bulk volume irreducible, or porosity.
  • FIG. 1 shows the different constituents of a fluid filled rock
  • FIG. 2 depicts diagrammatically a NMR logging tool conveyed in a borehole in an earth formation
  • FIG. 3 (prior art) show configurations of magnets, antenna and shield of a multifrequency NMR logging tool suitable for use with the present invention
  • FIG. 4 (prior art) is an example of a NMR pulse sequence at three frequencies
  • FIGS. 5 a and 5 b show a flow chart illustrating the method of the present invention
  • FIG. 6 shows an exemplary comparison of a T 2 distribution obtained using the method of the present invention with results obtained using prior art methods on a set of synthetic examples
  • FIG. 7 shows an exemplary comparison of the CBW distribution obtained using the method of the present invention with results obtained using prior art methods on a set of synthetic examples.
  • FIG. 2 depicts a borehole 10 which has been drilled in a typical fashion into a subsurface geological formation 12 to be investigated for potential hydrocarbon producing reservoirs.
  • An NMR logging tool 14 suitable for use with the present invention has been lowered into the hole 10 by means of a cable 16 and appropriate surface equipment represented diagrammatically by a reel 18 and is being raised through the formation 12 comprising a plurality of layers 12 a through 12 g of differing composition, to log one or more of the formation's characteristics.
  • the NMR logging tool is provided with bowsprings 22 to maintain the tool in an eccentric position within the borehole with one side of the tool in proximity to the borehole wall.
  • the permanent magnets used for providing the static magnetic field are indicated by 23 and the magnet configuration is that of a line dipole. Signals generated by the tool 14 are passed to the surface through the cable 16 and from the cable 16 through another line 19 to appropriate surface equipment 20 for processing, recording and/or display or for transmission to another site for processing, recording and/or display.
  • FIG. 3 schematically illustrates a preferred embodiment of the present invention wherein this shaping of the static and RF fields is accomplished.
  • the tool cross-sectional view in FIG. 3 illustrates a main magnet 217 , a second magnet 218 , and a transceiver antenna, comprising wires 219 and core material 210 .
  • the arrows 221 and 223 depict the polarization (e.g., from the South pole to the North pole) of the main magnet 217 and the secondary magnet 218 .
  • a noteworthy feature of the arrangement shown in FIG. 3 is that the polarization of the magnets providing the static field is towards the side of the tool, rather than towards the front of the tool (the right side of FIG. 3 ) as in prior art devices. The importance of this rotated configuration is discussed below.
  • the second magnet 218 is positioned to augment the shape of the static magnetic field by adding a second magnetic dipole in close proximity to the RF dipole defined by the wires 219 and the soft magnetic core 210 . This moves the center of the effective static dipole closer to the RF dipole, thereby increasing the azimuthal extent of the region of examination.
  • the second magnet 218 also reduces the shunting effect of the high permeability magnetic core 210 on the main magnet 217 : in the absence of the second magnet, the DC field would be effectively shorted by the core 210 .
  • the second magnet besides acting as a shaping magnet for shaping the static field to the front of the tool (the side of the main magnet) also acts as a bucking magnet with respect to the static field in the core 210 .
  • the bucking function and a limited shaping could be accomplished simply by having a gap in the core; however, since some kind of field shaping is required on the front side of the tool, in a preferred embodiment of the invention, the second magnet serves both for field shaping and for bucking. If the static field in the core 210 is close to zero, then the magnetostrictive ringing from the core is substantially eliminated.
  • the NMR instrument described above makes measurements of nuclear spin characteristics of a portion of the earth formation.
  • the magnet arrangement produces a static magnetic field that is spatially varying, particularly as a function of distance from the instrument into the earth formation.
  • the instrument is thus suited for making simultaneous measurements from a plurality of regions, each region characterized by a different field strength (and Larmor frequency) and field gradient. It is to be noted that by using a field shifting magnet, it is also possible to make measurements from the same region at different frequencies.
  • the block 401 represents a single CPMG sequence (or modified CPMG sequence with shortened refocusing pulse) at a first frequency f 1 .
  • This sequence has a typically TE 1 range from 0.2-1.0 ms with a total length of 0.5-1.0 s.
  • the block 403 represents a series of trainlets with a TE 2 being the shortest possible value by the instrument, usually in 0.2-0.5 ms range, NE 2 *TE 2 ⁇ 8 ms, TW ⁇ 30 ms, and number of trainlets NS>>1. This is acquired at the same frequency as block 401 and the phase is alternated between successive trainlets.
  • the block 405 at a second frequency f 2 comprises a short CPMG (or modified CPMG) sequence with TE 3 preferably selected to be the same as TE 1 , NE 3 *TE 3 ⁇ 40 ms, that are fully polarized. This is referred to hereafter as a short sequence.
  • the block 407 at the second frequency f 2 comprises a short sequence similar to 405 with TE 3 and identical acquisition length NE 3 *TE 3 ⁇ 40 ms, that are only partially polarized (Note the short wait time between the end of 405 and the beginning of 407 ) and that has a phase that is reversed relative to the phase of 405 .
  • the next component at the second frequency f 2 is the block 409 that is similar to block 403 but does not necessarily contain same number of trainlets, NS.
  • blocks 411 and 413 that are similar to 405 and 407 , respectively (i.e., short CPMG sequences with a long and short wait time respectively, polarity flipped between 411 and 413 ) and a block 415 similar to 409 .
  • this is merely an example of the variety of pulse sequences that may be processed using the method of the present invention and is not intended to be a limitation.
  • N echoes and M relaxation components There are N echoes and M relaxation components.
  • d Ax, (3), where d is an N ⁇ 1 column vector; A is an N ⁇ M matrix; and x is an M ⁇ 1 column matrix.
  • the column vector, T contains the corresponding relaxation times for the amplitudes in x.
  • T (T 21 ,T 22 , . . . ,T 2M ) T
  • This expression can be solved using many different methods from conventional inversion to singular value decomposition.
  • the weighting matrices are chosen to be any that satisfy the problem and a priori knowledge of the problem we are trying to solve.
  • the simultaneous problem is solved with M independent amplitudes.
  • a different subset of these amplitudes is used to fit the each one of the echo trains.
  • the T 2 for each amplitude does not necessarily need to be unique.
  • One echo train might use an amplitude with a T 2 of 64 msec and a different amplitude is used for a different echo train but with the same T 2 of 64 msec.
  • x i M i x, (9), where x i is a M i ⁇ 1 column matrix and M i is a M i ⁇ M with only a single one in each row that corresponds to an amplitude in the i th subspace and zeros elsewhere. In addition, M i ⁇ M.
  • Echo1 is the long T2 echo train. It is fully polarized. Echo2 is echo train from the CBW trainlets. Echo3 is the bound water echo trains. The bin times are 2 (2n ⁇ 3)/2 Parameter Echo1 Echo2 Echo3 Units NE 833 16 83 TE 0.6 0.6 0.6 Ms NA 1 50 4 TW Fully polarized 30 100 Ms Noise 2 0.28 1 Pu Bins 2 ( ⁇ 1,1,3, . . . 11) 2 ( ⁇ 1,1,3, . . . 7) 2 ( ⁇ 1.5, ⁇ 1,0, . . .
  • Echo1 All the relaxation components in Echo1 are fully polarized. Only some of the components in Echo2 and Echo3 are fully polarized.
  • the bin times are 2 (2n ⁇ 3)/2 ms where n is a positive integer.
  • the amplitudes to be fit simultaneously are found by comparing the polarization of each amplitude.
  • the CBW trainlets (Echo2) have the smallest TW and therefore the fewest components are fully polarized.
  • the relaxation times for those components must satisfy T 2 ⁇ TW 1.5 ⁇ 3 .
  • the factor 1.5 is the T 1 /T 2 ratio and it is assumed that TW>3T 1 to be fully polarized.
  • this ratio is not a limitation to the present invention and can be different for the different echo trains.
  • TW is the minimum TW and corresponds to the CBW echo trainlets or Echo2.
  • the cutoff value is 6.67. Therefore the components with T 2 ⁇ 6.67 msec for all echo trains are fully polarized and a common set of partial porosities is used to fit these components.
  • the remainder of the components that fit Echo2 are not fully polarized and require independent amplitudes.
  • Echo1 is fit by 7 amplitudes
  • Echo2 is fit by 4 amplitudes
  • Echo3 is fit by 5 amplitudes.
  • the total number of amplitudes is 16, but only 11 are independent.
  • eqns. (20)-(22) a superscript has been used to avoid any possible confusion with actual longitudinal and transverse relaxation times.
  • the largest time in eqn (20)-(22) is at least as long as the pulse sequence.
  • the elements of eqn. (18) are: (i) all the elements of the times from eqn. (20), (ii) the elements of the times from eqn. (21) that exceed the cutoff time of 6.67 ms, and (iii) the elements of the times from eqn. (22) that exceed the cutoff time of 22.2 ms.
  • the regularization matrix minimizes the curvature.
  • the optimum regularization coefficient is given as ⁇ i ⁇ ( A i ′ ⁇ x - d i ) T ⁇ ( A i ′ ⁇ x - d i ) NE i ⁇ ( M i ⁇ x ) T ⁇ ( M i ⁇ x ) .
  • the T 2 bins are defined 451 and depicted by 471 , 473 , 475 , 477 , 479 and 481 . To simplify the illustration, the bins are shown of equal size, but typically they are defined on a logarithmic scale as given by eqn. (18).
  • a plurality of echo trains 453 , 455 , . . . 457 is acquired.
  • the range of bin sizes for which a T 2 distribution may be obtained using the method of the present invention is determined from the characteristics of the individual echo trains. For the example of Table 1, the largest T 2 is determined from the time at which the first (fully polarized) echo train amplitude has an acceptable signal to noise ratio. The smallest T 2 is controlled by the smallest TE in any of the echo trains.
  • echo train 1 For echo train 1, a determination of a cutoff time is made 453 using eqn. (17). It should be noted that in eqn. (17), the factor of 1.5 is for exemplary purposes only and other factors could be used based on a priori knowledge. In this regard, the method of the present invention is quite different from prior art methods where a constant ratio of T 1 /T 2 is assumed. Based on this determined cutoff, echo train 1 is used only for bins 471 , 473 , 475 , 477 , 479 and 481 . This may be considered as a threshold test for the contribution of echo train to the final inversion. This is repeated for the other echo trains.
  • the cutoff corresponds to bin 475 while the maximum time for echo train 2 corresponds to bin 479 .
  • the bin vector for echo train has independent components for bins 475 , 477 and 479 . This process is repeated for the remaining echo trains.
  • echo train 1 In the case of the echo trains of Table 1, echo train 1 is fully polarized, hence it will contribute to all bins up to its maximum time.
  • the bin time vector given by eqn. (18) is a concatenation of:
  • Bin times for the fully polarized echo data were (0.35, 0.5, 0.71, . . . , 2048).
  • the bin times for Echo2 were (0.35, 0.5, 0.71 . . . 128).
  • the bin times for the bound water echo trains or Echo3 were (0.35, 0.5, 0.71 . . . 512).
  • the regularization matrix used was the curvature smoothing with zero amplitude boundary conditions.
  • the e n,j are noise amplitudes randomly generated from Gaussian distributions with RMS widths of (2, 0.28, 1) for echo train j.
  • T 2 was varied from 0.1 to 1000 ms.
  • One hundred different instances of the echo trains were generated for each T 2 in order to compute statistical measures of the inversion results.
  • the curve 501 is the result of separately inverting the three echo trains and splicing the results.
  • the curve 505 is another prior art method referred to as a joint inversion.
  • the Joint Inversion method uses a common set of partial porosities to invert the echo trains simultaneously. The difference in the polarization of the echo trains is accounted for by multiplying each partial porosity by a polarization factor. Because the polarization factor depends on T 1 instead of T 2 , a new parameter is introduced call the T 1 /T 2 ratio. The T 1 /T 2 ratio is optimized along with the partial porosities. Unlike the Separate Inversion method, it produces a smooth T 2 distribution.
  • the limitation to the joint inversion method is the assumption that the T 1 /T 2 ratio is constant.
  • the T 1 /T 2 ratio varies from about 1 to more than 3. It should be expected that this ratio will vary within the sensitive volume of an NMR logging tool.
  • FIG. 7 shows the CBW estimate for the separate inversion 551 , the multiple echo-train inversion of the present method 553 and the joint inversion 555 .
  • the three echo train inversion CBW estimate has improved characteristics at 0.5 ms and goes to zero above 4 ms faster than the separate or joint inversion. Having the complete T 2 distribution enables determination of CBW, BVI and porosity.
  • the processing of the measurements made by the probe in wireline applications may be done by the surface processor 20 or may be done by a downhole processor (not shown).
  • the processing may be done by a downhole processor that is part of a bottomhole assembly BHA conveyed on a tubular such as a drillstring or coiled tubing. This downhole processing reduces the amount of data that has to be telemetered. Alternatively, some or part of the data may be telemetered to the surface.
  • the measurements may be stored on a suitable memory device downhole and processed when the drillstring is tripped out of the borehole. Part of the processing may also be done at a remote location.
  • the operation of the NMR sensor may be controlled by the downhole processor and/or the surface processor. Implicit in the control and processing of the data is the use of a computer program implemented on a suitable machine readable medium that enables the processor to perform the control and processing.
  • the machine readable medium may include ROMs, EPROMs, EAROMs, Flash Memories and Optical disks.

Abstract

A method for inversion of multiple echo trains with different wait times uses a cutoff times for each of the echo trains for full polarization. Simultaneous inversion is carried out for T2 bins where full polarization exists

Description

    BACKGROUND OF THE INVENTION
  • 1. Field of the Invention
  • This invention is related to methods for acquiring nuclear magnetic resonance (NMR) measurements for determination of petrophysical properties of formations and properties of fluids therein. Specifically, the invention deals with simultaneous inversion of multiple echo trains for processing and interpreting Nuclear Magnetic Resonance (NMR) log data that exhibit relaxation and/or polarization.
  • 2. Description of the Related Art
  • Nuclear magnetic resonance is used in the oil industry, among others, and particularly in certain oil well logging tools. NMR instruments may be used for determining, among other things, the fractional volume of pore space and the fractional volume of mobile fluid filling the pore space of earth formations. Methods of using NMR measurements for determining the fractional volume of pore space and the fractional volume of mobile fluids are described, for example, in “Spin Echo Magnetic Resonance Logging: Porosity and Free Fluid Index Determination,” M. N. Miller et al., Society of Petroleum Engineers paper no. 20561, Richardson, Tex., 1990. Further description is provided in U.S. Pat. No. 5,585,720, of Carl M. Edwards, issued Dec. 17, 1996 and having the same assignee as the present application, entitled “Signal Processing Method For Multiexponentially Decaying Signals And Applications To Nuclear Magnetic Resonance Well Logging Tools.” The disclosure of that patent is incorporated herein by reference.
  • Deriving accurate relaxation spectra from nuclear magnetic resonance (NMR) data from logging subterranean formations is critical to determining total and effective porosities, irreducible water saturations, and permeabilities of the formations. U.S. Pat. No. 6,069,477 to Chen et al having the same assignee as the present application discusses the constituents of a fluid saturated rock and various porosities of interest. Referring to FIG. 1, the solid portion of the rock is made up of two components, the rock matrix and dry clay. The total porosity as measured by a density logging tool is the difference between the total volume and the solid portion. The total porosity includes clay-bound water (CBW), capillary bound water (also known as Bulk volume Irreducible or BVI), movable water and hydrocarbons. The effective porosity, a quantity of interest to production engineers, is the sum of the last three components and does not include the clay bound water. Accurate spectra are also essential to estimate the irreducible and the movable fluid volumes; distortion of partial porosity distributions that has been commonly observed for a variety of reasons, affects the estimates of these quantities. The reasons for the distortions to occur are mainly due to poor signal-to-noise ratio (SNR) and poor resolution in the time domain of the NMR data.
  • The most common NMR log acquisition and core measurement method employs T2 measurements using CPMG (Carr, Purcell, Meiboom and Gill) sequence, as taught by Meiboom and Gill in “Modified Spin-Echo Method for Measuring Nuclear Relaxation Time,” Rev. Sci. Instrum. 1958, 29, pp. 688-691. In this method, the echo data in any given echo train are collected at a fixed time interval, the interecho time (TE). Depending on the relaxation rate of the nuclear species under investigation in the underlying system, usually, a few hundred to a few thousand echoes are acquired to sample relaxation decay. Determining a light oil component, which has long relaxation time, requires taking several hundreds of ms of data while determination of CBW, which decays very fast, can be done with echo sequences of as short as a few tens of milliseconds.
  • There are numerous examples of NMR logging techniques used for obtaining information about earth formations and fluids. In measurement-while-drilling (MWD) operation, measurements are made while the wellbore is being drilled while in wireline logging, measurements are made after a wellbore has been drilled. The logging tools are lowered into the borehole and NMR signals are obtained using different configurations of magnets, transmitter coils and receiver coils. A static magnetic field is produced in the formation using permanent or electro-magnets. The static field aligns nuclear spins within the formation parallel to the static field. A pulsed RF field is applied using a transmitter on the logging tool and the nuclear magnetization signals produced by the pulsed RF field are analyzed to determine formation properties. The prior art shows different radio frequency (RF) pulsing schemes for generating RF fields in the formation. The most commonly used pulsing schemes are variations of the CPMG sequence denoted by
    (TW i,90±π/2,(τ,180τ,echo)j)i  (1)
    where TW is a wait time, 90 is a tipping pulse that tips the nuclear spins by an angle substantially equal to 90°, 180 is a refocusing pulse that tips the nuclear spins by an angle substantially equal to 180°, and echo is a spin echo. The time interval between successive refocusing pulses is 2τ, the number of echoes is j, and i denotes repetitions of the basic pulse sequence. A variation of the CPMG sequence is taught in U.S. Pat. No. 6,163,153 to Reiderman in which the use of a refocusing pulse with a tipping angle less than 180° is disclosed.
  • Rig time is expensive, so that the general objective in wireline logging is to obtain interpretable data within as short a time as possible. In MWD logging, on the other hand, no additional rig time is involved. However, when more measurements can be acquired in a given time, the data quality can be improved. The parameters that may be varied are the acquisition frequencies and the number of different frequencies, the tip angles, the wait time, the number of pulses within a CPMG sequence, and the time interval between the pulses. Long wait times are needed for proper evaluation of formation fluids that have long relaxation times, e.g., gas reservoirs while short wait times and/or short pulse spacings are used for evaluating faster relaxing components, e.g., irreducible fluid (BVI) and clay bound water (CBW). For example, U.S. Pat. No. 6,331,775 to Thern et al., having the same assignee as the present application discusses the use of a dual wait time acquisition for determination of gas saturation in a formation. U.S. Pat. No. 5,023,551 to Kleinberg et al discusses the use of CPMG sequences in well logging. U.S. Pat. No. 6,069,477 to Chen et al teaches the use of pulse sequences with different pulse spacings to determine CBW.
  • NMR fluid typing applications often involve acquiring multiple echo trains to exploit the diffusion and/or polarization contrasts between the water and hydrocarbon phases. Although proven successful in wells where these contrasts are large, the NMR-based techniques are challenging when the contrasts are small. Carbonates generally exhibit smaller NMR surface relaxivity than clastics, which reduces the relaxation time contrast between movable water and light oil. This difficulty is exacerbated by the small difference in the diffusivities of water and very light oil. In such cases, the methods for processing data are critical and the introduction of a priori information important.
  • U.S. patent application Ser. No. 10/288,115 of Chen et al., having the same assignee as the present invention and the contents of which are fully incorporated herein by reference, teaches the use of multifrequency NMR acquisition using various combinations of wait times, interecho times in pulse sequences of different lengths in an objective oriented method for formation and reservoir analysis. The objectives of Chen '115 may include formation evaluation, fluid typing, diffusivity determination The pulse sequences taught therein are suitable for use with the present invention. The mention of the Chen '115 application is not intended to be a limitation on the present invention and other types of pulse sequences could be used.
  • New generation NMR well logging tools can acquire multiple echo trains with different TWs and TEs in a single logging pass. Simultaneously inverting all echo data to obtain the different fluid T2 spectra is a logical method for the processing and interpretation of multiple echo trains acquired with different acquisition parameters. In fact, echo trains acquired with different TEs and TWs can not simply be averaged together to increase the signal-to-noise and inverted because they may exhibit very different T2 decay behaviors. Thus, it is desirable to use an analysis technique that accounts for the different acquisition parameters in the individual echo trains. U.S. patent application Ser. No. 10/435,419 of Chen, having the same assignee as the present invention and the contents of which are incorporated herein by reference, discloses an apparatus and a method of combining echo trains acquired with different parameters. Multiecho sequences are acquired from a first and second region of interest using a first and second radio frequency (RF) pulse sequence. A correction factor depending at least in part on a diffusivity of a fluid in the earth formation is determined, and the first and second multiecho sequences are combined using the correction factor to obtain a combined multiecho sequence. The method taught therein is difficult to extend to multiple echo trains, and the specific problem of small contrasts is not addressed. A similar approach (combining echo trains without gradient variations) is discussed in U.S. Pat. No. 6,377,042 issued to Menger. EP 0886792 to Bonnie et al. discloses a method in which NMR echo signals are acquired with different combinations of TW, TE and field gradient (produced by a gradient coil), and a curve fitting procedure is used to match the resulting signals to a predetermined model. The model itself is limited in scope.
  • Others have processed multiple echo trains of NMR data by separately inverting them and then by splicing the two separate inversion results. See, for example, U.S. Pat. No. 6,005,389 to Prammer. The methods of splicing generally assume that the ratio T1/T2 is constant.
  • There is a need for an apparatus and method of inversion of multiple echo trains of NMR data. Such a method should be robust and efficient. The present invention satisfies this need.
  • SUMMARY OF THE INVENTION
  • One embodiment of the present invention is a method of characterizing an earth formation having a fluid therein. The method includes obtaining a plurality of nuclear magnetic resonance (NMR) echo trains, each of the echo trains resulting from pulsing of the earth formation with an associated pulse sequence having an associated wait time. A plurality of bins of a T2 distribution of the earth formation is defined. For each of the plurality of echo trains, an associated cutoff time for full polarization is determined, and a value associated with each of the plurality of bins is determined based at least in part on the plurality of echo trains and the associated cutoff time. At least one of NMR echo trains may be a fully polarized echo train. The echo trains may be fully polarized echo trains, a clay bound water echo train, or a bulk volume irreducible echo train. Each of the associated cutoff times may be determined from a wait time of the associated pulse sequence and a maximum longitudinal relaxation time of the earth formation. A least squares inversion may be performed. The inversion may be a weighted inversion. Determination of the value associated with each of the plurality of bins involves a simultaneous fitting over each of the echo trains for a bin that is fully polarized for each of the plurality of echo trains and may further involve a simultaneous fitting of a subset of the echo trains for a bin that is fully polarized for the subset of echo trains. At least one characteristic of the earth formation selected from (i) clay bound water, (ii) bulk volume irreducible, and, (iii) porosity may be determined.
  • Another embodiment of the present invention is an apparatus for characterizing an earth formation having a fluid therein. The apparatus includes a nuclear magnetic resonance (NMR) tool conveyed in a borehole in the earth formation. The NMR tool pulses the earth formation with a plurality of radio frequency (RF) magnetic pulse sequences and receives a plurality of associated echo trains. A processor determines from the plurality of echo trains a value associated with each of a plurality of bins of a T2 distribution, the determination being based at least in part on cutoff times for full polarization associated with each of the plurality of echo trains. The echo trains may be a fully polarized echo trains, a clay bound water echo trains or a BVI echo train. The associated cutoff times may be determined from a wait time of the associated pulse sequence and a maximum longitudinal relaxation time of the earth formation. The processor may determine the value associated with each of the plurality of bins by further performing a least squares inversion. The inversion may be a weighted inversion. The inversion may involve a simultaneous fitting over each of the echo trains for a bin that is fully polarized for each of the plurality of echo trains and may further involve a simultaneous fitting of a subset of the echo trains for a bin that is fully polarized for the subset of echo trains. The processor may be at a surface location, a downhole location, or a remote location. The NMR tool may be conveyed in the borehole on a wireline, a drillstring, or coiled tubing. The processor may determine a characteristic of the earth formation such as clay bound water, bulk volume irreducible, and porosity.
  • Another embodiment of the present invention is a machine readable medium for use with an apparatus for characterizing an earth formation having a fluid therein. The medium includes instructions that enable a nuclear magnetic resonance (NMR) tool conveyed in a borehole to acquire a plurality of echo trains. The echo trains result from pulsing of the earth formation by an associated pulse sequence that has a wait time. The instructions also enable determination from the plurality of echo trains a value associated with each of a plurality of bins of a T2 distribution. The determination is based at least in part on cutoff times for full polarization associated with each of the plurality of echo trains. The instructions enable the processor to process a fully polarized echo train, a CBW echo train, and a BVI echo train. The instructions further enable performing performing a least squares inversion. The machine readable medium may be a ROM, an EPROM, an EAROMs, a Flash Memory, or an optical disk. The instructions may further enable determination of a characteristic of the formation such as clay bound water, bulk volume irreducible, or porosity.
  • BRIEF DESCRIPTION OF THE FIGURES
  • The present invention is best understood with reference to the accompanying figures in which like numerals refer to like elements and in which:
  • FIG. 1 (prior art) shows the different constituents of a fluid filled rock;
  • FIG. 2 (prior art) depicts diagrammatically a NMR logging tool conveyed in a borehole in an earth formation
  • FIG. 3 (prior art) show configurations of magnets, antenna and shield of a multifrequency NMR logging tool suitable for use with the present invention;
  • FIG. 4 (prior art) is an example of a NMR pulse sequence at three frequencies;
  • FIGS. 5 a and 5 b show a flow chart illustrating the method of the present invention;
  • FIG. 6 shows an exemplary comparison of a T2 distribution obtained using the method of the present invention with results obtained using prior art methods on a set of synthetic examples; and
  • FIG. 7 shows an exemplary comparison of the CBW distribution obtained using the method of the present invention with results obtained using prior art methods on a set of synthetic examples.
  • DETAILED DESCRIPTION OF THE INVENTION
  • FIG. 2 depicts a borehole 10 which has been drilled in a typical fashion into a subsurface geological formation 12 to be investigated for potential hydrocarbon producing reservoirs. An NMR logging tool 14 suitable for use with the present invention has been lowered into the hole 10 by means of a cable 16 and appropriate surface equipment represented diagrammatically by a reel 18 and is being raised through the formation 12 comprising a plurality of layers 12 a through 12 g of differing composition, to log one or more of the formation's characteristics. The NMR logging tool is provided with bowsprings 22 to maintain the tool in an eccentric position within the borehole with one side of the tool in proximity to the borehole wall. The permanent magnets used for providing the static magnetic field are indicated by 23 and the magnet configuration is that of a line dipole. Signals generated by the tool 14 are passed to the surface through the cable 16 and from the cable 16 through another line 19 to appropriate surface equipment 20 for processing, recording and/or display or for transmission to another site for processing, recording and/or display.
  • FIG. 3 schematically illustrates a preferred embodiment of the present invention wherein this shaping of the static and RF fields is accomplished. The tool cross-sectional view in FIG. 3 illustrates a main magnet 217, a second magnet 218, and a transceiver antenna, comprising wires 219 and core material 210. The arrows 221 and 223 depict the polarization (e.g., from the South pole to the North pole) of the main magnet 217 and the secondary magnet 218. A noteworthy feature of the arrangement shown in FIG. 3 is that the polarization of the magnets providing the static field is towards the side of the tool, rather than towards the front of the tool (the right side of FIG. 3) as in prior art devices. The importance of this rotated configuration is discussed below.
  • The second magnet 218 is positioned to augment the shape of the static magnetic field by adding a second magnetic dipole in close proximity to the RF dipole defined by the wires 219 and the soft magnetic core 210. This moves the center of the effective static dipole closer to the RF dipole, thereby increasing the azimuthal extent of the region of examination. The second magnet 218 also reduces the shunting effect of the high permeability magnetic core 210 on the main magnet 217: in the absence of the second magnet, the DC field would be effectively shorted by the core 210. Thus, the second magnet, besides acting as a shaping magnet for shaping the static field to the front of the tool (the side of the main magnet) also acts as a bucking magnet with respect to the static field in the core 210. Those versed in the art would recognize that the bucking function and a limited shaping could be accomplished simply by having a gap in the core; however, since some kind of field shaping is required on the front side of the tool, in a preferred embodiment of the invention, the second magnet serves both for field shaping and for bucking. If the static field in the core 210 is close to zero, then the magnetostrictive ringing from the core is substantially eliminated.
  • Additional details of the logging tool and the field shaping are disclosed in U.S. Pat. No. 6,348,792 of Beard et al., having the same assignee as the present invention and the contents of which are fully incorporated herein by reference. As discussed in Beard, within the region of investigation, the static field gradient is substantially uniform and the static field strength lies within predetermined limits to give a substantially uniform Larmor frequency. Those versed in the art would recognize that the combination of field shaping and bucking could be accomplished by other magnet configurations than those shown in FIG. 3.
  • The NMR instrument described above makes measurements of nuclear spin characteristics of a portion of the earth formation. The magnet arrangement produces a static magnetic field that is spatially varying, particularly as a function of distance from the instrument into the earth formation. The instrument is thus suited for making simultaneous measurements from a plurality of regions, each region characterized by a different field strength (and Larmor frequency) and field gradient. It is to be noted that by using a field shifting magnet, it is also possible to make measurements from the same region at different frequencies.
  • Turning now to FIG. 4, an example of a three frequency acquisition is shown. The block 401 represents a single CPMG sequence (or modified CPMG sequence with shortened refocusing pulse) at a first frequency f1. This sequence has a typically TE1 range from 0.2-1.0 ms with a total length of 0.5-1.0 s. The block 403 represents a series of trainlets with a TE2 being the shortest possible value by the instrument, usually in 0.2-0.5 ms range, NE2*TE2˜8 ms, TW˜30 ms, and number of trainlets NS>>1. This is acquired at the same frequency as block 401 and the phase is alternated between successive trainlets. This is referred to hereafter as a trainlet sequence. The block 405 at a second frequency f2 comprises a short CPMG (or modified CPMG) sequence with TE3 preferably selected to be the same as TE1, NE3*TE3˜40 ms, that are fully polarized. This is referred to hereafter as a short sequence. The block 407 at the second frequency f2 comprises a short sequence similar to 405 with TE3 and identical acquisition length NE3*TE3˜40 ms, that are only partially polarized (Note the short wait time between the end of 405 and the beginning of 407) and that has a phase that is reversed relative to the phase of 405. The next component at the second frequency f2 is the block 409 that is similar to block 403 but does not necessarily contain same number of trainlets, NS. Moving to the third frequency f3, shown are blocks 411 and 413 that are similar to 405 and 407, respectively (i.e., short CPMG sequences with a long and short wait time respectively, polarity flipped between 411 and 413) and a block 415 similar to 409. It should be noted that this is merely an example of the variety of pulse sequences that may be processed using the method of the present invention and is not intended to be a limitation.
  • In the inversion process, the inversion involves finding amplitudes {xi} such that d ( t j ) = i = 1 M x i - t j / T 2 i , ( 2 )
    where d(tj) is the echo amplitude at time tj, and {T2i} are the relaxation times associated with the {xi}. There are N echoes and M relaxation components. In matrix notation
    d=Ax,  (3),
    where d is an N×1 column vector; A is an N×M matrix; and x is an M×1 column matrix. The column vector, T, contains the corresponding relaxation times for the amplitudes in x.
    T=(T21,T22, . . . ,T2M)T
    The elements of A are computed as
    A i,j =e −t i /T 2j   (4).
    The inversion involves minimization of
    σ2=(Ax−d)T W(Ax−d)+αx T W m x,  (5),
    where W is a data weighting matrix and Wm is a regularization matrix representing are prior knowledge about the result x, and a is the regularization parameter. The resulting expression for x is
    0=A T WAx+αW m x−A T Wd  (6).
    This expression can be solved using many different methods from conventional inversion to singular value decomposition. The weighting matrices are chosen to be any that satisfy the problem and a priori knowledge of the problem we are trying to solve.
  • Separate inversion solves two or more of the above expressions. Joint inversion concatenates two or more of the above expressions. If two echo trains are to be inverted, then the concatenation/minimization may be represented by A = [ A 1 A 2 ] ; W = [ W 1 0 0 W 2 ] ; d = [ d 1 d 2 ] . ( 7 )
    For additional echo trains, the expressions are easy to generalize. For partially polarized amplitudes where the TW are small,
    A (1,2)i,j =e −t I /T 2,j (1−e −TW (1,2) /rT 2,j )  (8).
  • In the present invention, the simultaneous problem is solved with M independent amplitudes. A different subset of these amplitudes is used to fit the each one of the echo trains. The T2 for each amplitude does not necessarily need to be unique. One echo train might use an amplitude with a T2 of 64 msec and a different amplitude is used for a different echo train but with the same T2 of 64 msec. Let the vector of independent amplitudes be x. Then we can map the subset of xi used for each echo train by a mapping matrix.
    x i =M i x,  (9),
    where xi is a Mi×1 column matrix and Mi is a Mi×M with only a single one in each row that corresponds to an amplitude in the ith subspace and zeros elsewhere. In addition, Mi≦M. The column vector T contains the corresponding relaxation times
    T=(T2,1,T2,2, . . . ,T2,M)T.
  • Substituting from eqn. (9) into eqn. (5) gives σ 2 = i ( A i x i - d i ) T W i ( A i x i - d i ) + α i x i T W m , i x i , = i ( A i M i x - d i ) T W i ( A i M i x - d i ) + α i ( M i x ) T W m , i M i x . ( 10 )
    Each Ai is an Ni×Mi matrix where each echo train has Ni echoes. It has elements
    (A k)i,j =e −t i /(T k ) j   (11),
    where the relaxation times are taken from the vector
    T k =M k T.  (12).
    Taking the derivative of eqn. (10) with respect to x results in an equation that can be solved for the amplitudes: i M i T A i T W i d i = [ i ( M i T A i T W i A i M i + α i M i T W m , i M i ) ] x . ( 13 )
    This can be further simplified by setting
    A i′ =A i M i  (14),
    and
    W m,i ′=M i T W m,i M i  (15).
    Eqn. (13) then becomes i A i T W i d i = [ i ( A i T W i A i + α i A m , i ) ] x . ( 16 )
    Eqn. (16) can be inverted using prior art methods.
  • A specific example is given to illustrate the multiple echo train inversion and the structure of the various matrices. The acquisition parameters for a three-frequency acquisition are shown in Table 1.
    TABLE 1
    Echo1 is the long T2 echo train. It is fully polarized.
    Echo2 is echo train from the CBW trainlets. Echo3 is the
    bound water echo trains. The bin times are 2(2n−3)/2
    Parameter Echo1 Echo2 Echo3 Units
    NE 833 16 83
    TE 0.6 0.6 0.6 Ms
    NA
    1 50 4
    TW Fully polarized 30 100 Ms
    Noise
    2 0.28 1 Pu
    Bins
    2(−1,1,3, . . . 11) 2(−1,1,3, . . . 7) 2(−1.5,−1,0, . . . 9)

    All the relaxation components in Echo1 are fully polarized. Only some of the components in Echo2 and Echo3 are fully polarized. We seek to fit the echo trains with amplitudes that have the bin times given in the table. The bin times are 2(2n−3)/2 ms where n is a positive integer.
  • The amplitudes to be fit simultaneously are found by comparing the polarization of each amplitude. The CBW trainlets (Echo2) have the smallest TW and therefore the fewest components are fully polarized. Thus, the relaxation times for those components must satisfy T 2 TW 1.5 · 3 . ( 17 )
    Here, the factor 1.5 is the T1/T2 ratio and it is assumed that TW>3T1 to be fully polarized. As noted above, this ratio is not a limitation to the present invention and can be different for the different echo trains. For this particular case TW is the minimum TW and corresponds to the CBW echo trainlets or Echo2. The cutoff value is 6.67. Therefore the components with T2<6.67 msec for all echo trains are fully polarized and a common set of partial porosities is used to fit these components. The remainder of the components that fit Echo2 are not fully polarized and require independent amplitudes.
  • To find the amplitudes that are common between Echo1 and Echo3 we simply substitute the TW of Echo3 or the BW echo trains. The cutoff is 22.2 ms for Echo3. Thus components with 6.67<T2<22.2 ms use a common set of partial porosities to fit both Echo1 and Echo3. The remainder of components with T2>22.2 ms that fit Echo3 are not fully polarized and the partial porosities are independent. The remainder of the components that fit Echo1 are fully polarized and fit with independent partial porosities. To obtain the fully polarized T2 spectrum we take the fully polarized independent amplitudes. These amplitudes fit Echo1.
  • Echo1 is fit by 7 amplitudes, Echo2 is fit by 4 amplitudes, and Echo3 is fit by 5 amplitudes. The total number of amplitudes is 16, but only 11 are independent. The bin time vector becomes
    T=(0.5, 2, 8, 32, 128, 512, 2048, 8, 32, 32, 128)T  (18).
    The mapping matrices are: M 1 = ( 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 ) M 2 = ( 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 ) M 3 = ( 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 ) . ( 19 )
    The elements of the forward models are ( A 1 ) i , j = - TE 1 / ( T 1 ) j ( 20 ) ( T 1 ) j { .5 , 2 , 8 , 32 , 128 , 512 , 2048 } i { 1 , 2 , 3 , , NE 1 } , ( A 2 ) i , j = - TE 2 / ( T 2 ) j ( 21 ) ( T 2 ) j { .5 , 2 , 8 , 32 } i { 1 , 2 , 3 , , NE 2 } , and ( A 3 ) i , j = - TE 3 / ( T 3 ) j ( 22 ) ( T 3 ) j { .5 , 2 , 8 , 32 , 128 } i { 1 , 2 , 3 , , NE 3 } .
    In eqns. (20)-(22), a superscript has been used to avoid any possible confusion with actual longitudinal and transverse relaxation times. The largest time in eqn (20)-(22) is at least as long as the pulse sequence. Note that the elements of eqn. (18) are:
    (i) all the elements of the times from eqn. (20),
    (ii) the elements of the times from eqn. (21) that exceed the cutoff time of 6.67 ms, and
    (iii) the elements of the times from eqn. (22) that exceed the cutoff time of 22.2 ms.
    The weighting matrices are NEi×NEi diagonal matrices with values W = 1 σ e 2 [ σ e 2 / σ 1 2 1 1 ] , ( 23 )
    where σ1 is the noise on the first echo and σe is the noise on the remainder of the echoes. In our case the ratio is σ1e=1.6. The regularization matrix minimizes the curvature. The curvature matrices with zero amplitude boundary conditions are Mi×Mi tridiagonal matrices with values C = [ - 2 1 1 - 2 1 1 - 2 1 - 1 2 ] , ( 24 ) W m = C T C . ( 25 )
    The optimum regularization coefficient is given as α i ( A i x - d i ) T ( A i x - d i ) NE i ( M i x ) T ( M i x ) . ( 26 )
    This result follows from a result in Butler et al., “Estimating Solutions of First Kind Integral Equations with Non-negative Constraints and Optimal Smoothing,” SIAM J. Numerical Analysis, vol. 18, no. 3, pp 381-397.
  • To summarize, in the multiple echo train inversion, an assumption is made that some upper limit exists for the T1/T2 ratio and is known a priori. There is no assumption that the ratio T1/T2 is constant. In addition, all the available data are used to compute an optimized total porosity T2 spectrum. The criterion for determining which components are fit simultaneously compares TW to the component's relaxation time. If the component is fully polarized for all the echo trains it is fit simultaneously. If the component is fully polarized for a subset of the echo trains, it is fit simultaneously only to those echo trains. The method uses non-negative least-square formalism to obtain the T2 spectrum. The formalism demonstrates that only an M×M matrix need be inverted, where M is the number of independent amplitudes in the inversion. This matrix is the sum of matrices calculated in the separate inversion process.
  • This is depicted schematically in FIG. 5 a. The T2 bins are defined 451 and depicted by 471, 473, 475, 477, 479 and 481. To simplify the illustration, the bins are shown of equal size, but typically they are defined on a logarithmic scale as given by eqn. (18). A plurality of echo trains 453, 455, . . . 457 is acquired. The range of bin sizes for which a T2 distribution may be obtained using the method of the present invention is determined from the characteristics of the individual echo trains. For the example of Table 1, the largest T2 is determined from the time at which the first (fully polarized) echo train amplitude has an acceptable signal to noise ratio. The smallest T2 is controlled by the smallest TE in any of the echo trains.
  • For echo train 1, a determination of a cutoff time is made 453 using eqn. (17). It should be noted that in eqn. (17), the factor of 1.5 is for exemplary purposes only and other factors could be used based on a priori knowledge. In this regard, the method of the present invention is quite different from prior art methods where a constant ratio of T1/T2 is assumed. Based on this determined cutoff, echo train 1 is used only for bins 471, 473, 475, 477, 479 and 481. This may be considered as a threshold test for the contribution of echo train to the final inversion. This is repeated for the other echo trains.
  • Turning to FIG. 5 b, for echo train 2, the cutoff corresponds to bin 475 while the maximum time for echo train 2 corresponds to bin 479. Hence the bin vector for echo train has independent components for bins 475, 477 and 479. This process is repeated for the remaining echo trains.
  • In the case of the echo trains of Table 1, echo train 1 is fully polarized, hence it will contribute to all bins up to its maximum time. The bin time vector given by eqn. (18) is a concatenation of:
  • (i) all bins corresponding to the fully polarized echo trains,
  • (ii) partially polarized bins for echo train 2, i.e, the range between its cutoff time and its maximum time; and
  • (iii) partially polarized bins for echo train 3, i.e., the range between its cutoff time and its maximum time.
  • Once the vector T has been defined, the rest is straightforward.
  • To evaluate the characteristics of the method of the present invention, a synthetic example is shown. Bin times for the fully polarized echo data were (0.35, 0.5, 0.71, . . . , 2048). The bin times for Echo2 were (0.35, 0.5, 0.71 . . . 128). The bin times for the bound water echo trains or Echo3 were (0.35, 0.5, 0.71 . . . 512). The regularization matrix used was the curvature smoothing with zero amplitude boundary conditions. Finally, the synthetic data were constructed using a single exponential such that
    d n,j=10 exp {−nTE/T 2}(1−exp {−TW j/1.5T 2})+e n,j,
    where j represents the echo train number. The en,j are noise amplitudes randomly generated from Gaussian distributions with RMS widths of (2, 0.28, 1) for echo train j. T2 was varied from 0.1 to 1000 ms. One hundred different instances of the echo trains were generated for each T2 in order to compute statistical measures of the inversion results.
  • Shown in FIG. 6 are the total porosity estimates 503 using the method of the present invention. The curve 501 is the result of separately inverting the three echo trains and splicing the results. The curve 505 is another prior art method referred to as a joint inversion. The Joint Inversion method uses a common set of partial porosities to invert the echo trains simultaneously. The difference in the polarization of the echo trains is accounted for by multiplying each partial porosity by a polarization factor. Because the polarization factor depends on T1 instead of T2, a new parameter is introduced call the T1/T2 ratio. The T1/T2 ratio is optimized along with the partial porosities. Unlike the Separate Inversion method, it produces a smooth T2 distribution. However, the limitation to the joint inversion method is the assumption that the T1/T2 ratio is constant. For core plugs, the T1/T2 ratio varies from about 1 to more than 3. It should be expected that this ratio will vary within the sensitive volume of an NMR logging tool.
  • Still with reference to the synthetic model, FIG. 7 shows the CBW estimate for the separate inversion 551, the multiple echo-train inversion of the present method 553 and the joint inversion 555. The three echo train inversion CBW estimate has improved characteristics at 0.5 ms and goes to zero above 4 ms faster than the separate or joint inversion. Having the complete T2 distribution enables determination of CBW, BVI and porosity.
  • An alternate embodiment of the invention starts with the same cost function as that given by eqn. (5):
    σ2=(Ax−d)T W(Ax−d)+αx T W m x,  (27),
    where A i , j = - t i T 2 , j . ( 28 )
    Now suppose that we want to use the same coefficient for multiple inversions. The fully polarized components are included in all the inversions. Thus, any bin time less than 3*TW is included all inversions. Thus, bins of the T2 distribution for which T2,i≦4 ms are included in the CBW, BVI, and FP inversions. Bins of the T2 distribution for which 4<T2,i≦16 ms are included in the BVI and FP inversions. Finally, only those bins larger than 16 ms are included in the FP inversion. Also important is the notion, that to completely fit the partially polarized data we need to include bins for the partially polarized part of the echoes. Thus
    T2cbw,i={0.35,0.5, . . . ,4,5.6, . . . ,256}
    T 2bvi,i={0.35,0.5, . . . ,16,22.6, . . . ,256}
    T 2fp,i={0.35,0.5, . . . ,2048}
    and
    xi={x0.35, . . . ,x2048,xbvi,22.6, . . . ,xbvi,256,xcbw,5.6, . . . xcbw,256}
    T2,i={0.35, . . . ,2048,22.6, . . . ,256,5.6, . . . ,256}
    ti={tfp,1, . . . ,tfp,N fp ,tbvi,1, . . . ,tbvi,N bvi ,tcbw,1, . . . ,tcbw,N cbw}
    Those labeled with CBW and BVI are used only to fit the data from those echo trains and discarded. The final product are the unlabeled bins. Let x1 to xJ be the final product and let xJ+1 to xJ+K be those fit only to the BVI echo train and xJ+K+1 to xJ+K+L be fit only to the CBW echo trains. Thus the A matrix becomes A = ( - t 1 / T 2 , 1 - t 1 / T 2 , J - t N fp / T 2 , 1 - t N fp / T 2 , J 0 0 - t N fp / T 2 , 1 - t N fp + 1 / T 2 , J bwi - t N bwi + N fp / T 2 , 1 - t N bwi + N fp / T 2 , J bwi 0 - t N + t / T 2 , J + 1 - t N + 1 / T 2 , J + K - t N bwi + N fp / T 2 , J + 1 - t N bwi + N fp / T 2 , J + K 0 - t 1 / T 2 , 1 - t 1 / T 2 , J , cbw - t N cbw + N bwi + N fp / T 2 , 1 - t N cbw + N bwi + N fp / T 2 , J , cbw 0 0 - t N bwi + N fp + 1 / T 2 , J + K + 1 - t N bwi + N fp + 1 / T 2 , J , + K + 1 - t N cbw + N bwi + N fp / T 2 , J + K + 1 - t N cbw + N bwi + N fp / T 2 , J , + K + 1 ) A = ( [ N fp × J ] 0 0 [ N bwi × J bwi ] 0 [ N bwi × K ] 0 [ N cbw × J cbw ] 0 0 [ N cbw × L ] )
    Constraints such as minimum norm or curvature may be added. The minimum norm matrix for a single inversion is given by W m = ( 1 0 0 1 0 0 1 0 0 1 0 0 1 0 0 1 )
    This could be used straightforwardly. Different weighting factors could be applied to those coefficients that are used in the fit of more than one echo train. The curvature smoothing constraint is less easy. For zero amplitude boundary conditions, C = ( - 2 1 1 - 2 1 1 - 2 1 1 - 2 1 1 - 2 1 1 - 2 ) W m = C T C
    The regularization conditions need to be kept separate because they are squared before they are summed. Thus we get to solve this equation ( i α i W m , i + A T WA ) x = A T Wd
  • The processing of the measurements made by the probe in wireline applications may be done by the surface processor 20 or may be done by a downhole processor (not shown). For MWD applications, the processing may be done by a downhole processor that is part of a bottomhole assembly BHA conveyed on a tubular such as a drillstring or coiled tubing. This downhole processing reduces the amount of data that has to be telemetered. Alternatively, some or part of the data may be telemetered to the surface. In yet another alternative, the measurements may be stored on a suitable memory device downhole and processed when the drillstring is tripped out of the borehole. Part of the processing may also be done at a remote location.
  • The operation of the NMR sensor may be controlled by the downhole processor and/or the surface processor. Implicit in the control and processing of the data is the use of a computer program implemented on a suitable machine readable medium that enables the processor to perform the control and processing. The machine readable medium may include ROMs, EPROMs, EAROMs, Flash Memories and Optical disks.
  • While the foregoing disclosure is directed to the specific embodiments of the invention, various modifications will be apparent to those skilled in the art. It is intended that all such variations within the scope and spirit of the appended claims be embraced by the foregoing disclosure.

Claims (29)

1. A method of characterizing an earth formation, the method comprising:
(a) obtaining a plurality of nuclear magnetic resonance (NMR) echo trains,
(b) defining a plurality of bins of a T2 distribution of the earth formation; and
(c) determining a value associated with each of the plurality of echo trains by simultaneously inverting the plurality of echo trains.
2. The method of claim 1 wherein at least one of NMR echo trains is a fully polarized echo train.
3. The method of claim 1 wherein the plurality of echo trains are selected from the group consisting of (i) a fully polarized echo train, (ii) a CBW echo train, and, (iii) a BVI echo train.
4. The method of claim 27 wherein each of the associated cutoff times is determined from a wait time of an associated pulse sequence and a maximum longitudinal relaxation time of the earth formation.
5. The method of claim 1 wherein determining the value associated with each of the plurality of bins further comprises performing a least squares inversion.
6. The method of claim 5 wherein the inversion comprises a weighted inversion.
7. The method of claim 1 wherein determining the value associated with each of the plurality of bins further comprises a simultaneous fitting over each of the echo trains for a bin that is fully polarized for each of the plurality of echo trains.
8. The method of claim 1 wherein determining the value associated with each of the plurality of bins further comprises a simultaneous fitting of a subset of the echo trains for a bin that is fully polarized for the subset of echo trains.
9. The method of claim 2 wherein determining the value associated with each of the plurality of bins further comprises defining a vector that includes bins associated with the fully polarized echo train.
10. The method of claim 1 further comprising determining at least one characteristic of the earth formation selected from (i) clay bound water, (ii) bulk volume irreducible, and, (iii) porosity.
11. An apparatus for characterizing an earth formation, the apparatus comprising:
(a) a nuclear magnetic resonance (NMR) tool conveyed in a borehole in the earth formation, the NMR tool pulsing the earth formation with a plurality of radio frequency (RF) magnetic pulse sequences and receiving a plurality of associated echo trains;
(b) a processor which determines from the plurality of echo trains a value associated with each of a plurality of bins of a T2 distribution based at least in Dar be simultaneously inverting the echo trains.
12. The apparatus of claim 11 wherein the plurality of echo trains are selected from the group consisting of (i) a fully polarized echo train, (ii) a CBW echo train, and, (iii) a BVI echo train.
13. The apparatus of claim 29 wherein each of the associated cutoff times is determined from a wait time of the associated pulse sequence and a maximum longitudinal relaxation time of the earth formation.
14. (canceled)
15. The apparatus of claim 11 wherein the inversion comprises a weighted inversion.
16. The apparatus of claim 11 wherein the processor determines the value associated with each of the plurality of bins by performing a simultaneous fitting over each of the echo trains for a bin that is fully polarized for each of the plurality of echo trains.
17. The apparatus of claim 11 wherein the processor determines the value associated with each of the plurality of bins by performing a simultaneous fitting of a subset of the echo trains for a bin that is folly polarized for the subset of echo trains.
18. The apparatus of claim 11 wherein the processor determines the value associated with each of the plurality of bins by defining a vector that includes bins associated with a fully polarized echo train.
19. The apparatus of claim 11 wherein a location of the processor is selected from the group consisting of (i) a surface location, (ii) a downhole location, and, (iii) a remote location.
20. The apparatus of claim 11 further comprising a conveyance device for the NMR tool, the conveyance device selected from the group consisting of: (i) a wireline, (ii) a drillstring, and, (iii) coiled tubing.
21. The apparatus of claim 11 wherein the processor further determines a characteristic of the earth formation selected from (i) clay bound water, (ii) bound volume irreducible, and, (iii) porosity.
22. A machine readable medium for use with an apparatus for characterizing an earth formation, the apparatus comprising
(a) a nuclear magnetic resonance (NMR) tool conveyed in a borehole in the earth formation to acquire a plurality of echo trains, each of the echo trains resulting from pulsing of the earth formation with an associated pulse sequence having an associated wait time;
the medium comprising instructions which enable a processor to:
(b) simultaneously invert the plurality of echo trains and determine a value associated with each of a plurality of bins of a T2 distribution.
23. The machine readable medium of claim 22 wherein instructions further comprise enabling the processor to process echo trains are selected from the group consisting of (i) a fully polarized echo train, (ii) a CBW echo train, and, (iii) a BVI echo train.
24. (canceled)
25. The machine readable medium of claim 22 wherein the medium is selected from the group consisting of (i) ROMs, (ii) EPROMs, (iii) EAROMs, (iv) Flash Memories, and, (v) Optical disks.
26. The machine readable medium of claim 22 wherein the medium further comprises instructions for determining a characteristic of the formation selected from (i) clay bound water, (ii) bulk volume irreducible, and, (iii) porosity.
27. The method of claim 1 wherein inverting the plurality of echo trains further comprises:
(i) for each of the plurality of echo trains, determining an associated cutoff time for full polarization, and
(ii) determining the value associated with each of the plurality of bins based at least in part on the associated cutoff time.
28. The method of claim 1 wherein each of the echo trains results from pulsing of the earth formation with an associated pulse sequence having an associated wait time.
29. The apparatus of claim 11 wherein the processor inverts the plurality of echo trains by further:
(i) determining a cutoff time for full polarization of each of the plurality of echo trains, and
(ii) determining the value associated with each of the plurality of bins using the associated cut off time.
US11/037,834 2005-01-18 2005-01-18 Multiple echo train inversion Abandoned US20060158184A1 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US11/037,834 US20060158184A1 (en) 2005-01-18 2005-01-18 Multiple echo train inversion
GB0526441A GB2422198B (en) 2005-01-18 2005-12-23 Multiple echo train inversion
US11/781,522 US7812602B2 (en) 2005-01-18 2007-07-23 Multiple echo train inversion

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US11/037,834 US20060158184A1 (en) 2005-01-18 2005-01-18 Multiple echo train inversion

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US11/781,522 Continuation US7812602B2 (en) 2005-01-18 2007-07-23 Multiple echo train inversion

Publications (1)

Publication Number Publication Date
US20060158184A1 true US20060158184A1 (en) 2006-07-20

Family

ID=35841233

Family Applications (2)

Application Number Title Priority Date Filing Date
US11/037,834 Abandoned US20060158184A1 (en) 2005-01-18 2005-01-18 Multiple echo train inversion
US11/781,522 Active 2026-05-28 US7812602B2 (en) 2005-01-18 2007-07-23 Multiple echo train inversion

Family Applications After (1)

Application Number Title Priority Date Filing Date
US11/781,522 Active 2026-05-28 US7812602B2 (en) 2005-01-18 2007-07-23 Multiple echo train inversion

Country Status (2)

Country Link
US (2) US20060158184A1 (en)
GB (1) GB2422198B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009038744A1 (en) * 2007-09-18 2009-03-26 Baker Hughes Incorporated Nuclear magnetic resonance evaluation using independent component analysis (ica)-based blind source separation
US20140278113A1 (en) * 2013-03-14 2014-09-18 Hamed Chok Real-Time Determination of Formation Fluid Properties Using Density Analysis
WO2018039054A1 (en) * 2016-08-23 2018-03-01 Baker Hughes, A Ge Company, Llc Simultaneous inversion of nmr multiple echo trains and conventional logs
US20180308502A1 (en) * 2017-04-20 2018-10-25 Thomson Licensing Method for processing an input signal and corresponding electronic device, non-transitory computer readable program product and computer readable storage medium
CN111980663A (en) * 2020-07-21 2020-11-24 中海油田服务股份有限公司 Multi-frequency multi-dimensional nuclear magnetic logging method and device

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7565246B2 (en) * 2007-03-22 2009-07-21 Baker Hughes Incorporated Determination of gas saturation radial profile from multi-frequency NMR data
US9551213B2 (en) * 2009-04-07 2017-01-24 Baker Hughes Incorporated Method for estimation of bulk shale volume in a real-time logging-while-drilling environment
US8614578B2 (en) * 2009-06-18 2013-12-24 Schlumberger Technology Corporation Attenuation of electromagnetic signals passing through conductive material
US9671483B2 (en) * 2014-02-26 2017-06-06 Baker Hughes Incorporated T2 inversions with reduced motion artifacts
US10061053B2 (en) 2015-04-30 2018-08-28 Baker Hughes, A Ge Company, Llc NMR T2 distribution from simultaneous T1 and T2 inversions for geologic applications
US10739489B2 (en) * 2016-01-15 2020-08-11 Baker Hughes, A Ge Company, Llc Low gradient magnetic resonance logging for measurement of light hydrocarbon reservoirs
US10429536B2 (en) 2016-04-04 2019-10-01 Baker Hughes, A Ge Company, Llc T2 inversions with reduced motion artifacts

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5023551A (en) * 1986-08-27 1991-06-11 Schlumberger-Doll Research Nuclear magnetic resonance pulse sequences for use with borehole logging tools
US5585720A (en) * 1995-10-23 1996-12-17 Western Atlas International, Inc. Signal processing method for multiexponentially decaying signals and application to nuclear magnetic resonance well logging tools
US6069477A (en) * 1997-09-05 2000-05-30 Western Atlas International, Inc. Method for improving the accuracy of NMR relaxation distribution analysis with two echo trains
US6163153A (en) * 1998-09-11 2000-12-19 Western Atlas International, Inc. Nuclear magnetic resonance pulse sequence for optimizing instrument electrical power usage
US6331775B1 (en) * 1999-09-15 2001-12-18 Baker Hughes Incorporated Gas zone evaluation by combining dual wait time NMR data with density data
US6366087B1 (en) * 1998-10-30 2002-04-02 George Richard Coates NMR logging apparatus and methods for fluid typing
US6377042B1 (en) * 1998-08-31 2002-04-23 Numar Corporation Method and apparatus for merging of NMR echo trains in the time domain
US20030214286A1 (en) * 2002-05-16 2003-11-20 Ralf Heidler Method for the inversion of CPMG measurements enhanced by often repeated short wait time measurements
US20050264285A1 (en) * 2004-05-27 2005-12-01 Baker Hughes Incorporated Method of detecting, quantifying and correcting borehole contaminations from multi-frequency, multi-sensitive-volume NMR logging data

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
MY122012A (en) 1996-03-14 2006-03-31 Shell Int Research Determining a fluid fraction in an earth formation
US6232778B1 (en) * 1998-06-11 2001-05-15 Schlumberger Technology Corporation Method for obtaining NMR bound fluid volume using partial polarization
US6255819B1 (en) * 1999-10-25 2001-07-03 Halliburton Energy Services, Inc. System and method for geologically-enhanced magnetic resonance imaging logs
US6348792B1 (en) 2000-07-27 2002-02-19 Baker Hughes Incorporated Side-looking NMR probe for oil well logging

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5023551A (en) * 1986-08-27 1991-06-11 Schlumberger-Doll Research Nuclear magnetic resonance pulse sequences for use with borehole logging tools
US5585720A (en) * 1995-10-23 1996-12-17 Western Atlas International, Inc. Signal processing method for multiexponentially decaying signals and application to nuclear magnetic resonance well logging tools
US6069477A (en) * 1997-09-05 2000-05-30 Western Atlas International, Inc. Method for improving the accuracy of NMR relaxation distribution analysis with two echo trains
US6377042B1 (en) * 1998-08-31 2002-04-23 Numar Corporation Method and apparatus for merging of NMR echo trains in the time domain
US6163153A (en) * 1998-09-11 2000-12-19 Western Atlas International, Inc. Nuclear magnetic resonance pulse sequence for optimizing instrument electrical power usage
US6366087B1 (en) * 1998-10-30 2002-04-02 George Richard Coates NMR logging apparatus and methods for fluid typing
US20030016012A1 (en) * 1998-10-30 2003-01-23 Coates George Richard NMR logging apparatus and methods for fluid typing
US6825658B2 (en) * 1998-10-30 2004-11-30 George Richard Coates NMR logging apparatus and methods for fluid typing
US6331775B1 (en) * 1999-09-15 2001-12-18 Baker Hughes Incorporated Gas zone evaluation by combining dual wait time NMR data with density data
US20030214286A1 (en) * 2002-05-16 2003-11-20 Ralf Heidler Method for the inversion of CPMG measurements enhanced by often repeated short wait time measurements
US6714009B2 (en) * 2002-05-16 2004-03-30 Schlumberger Technology Corporation Method for the inversion of CPMG measurements enhanced by often repeated short wait time measurements
US20050264285A1 (en) * 2004-05-27 2005-12-01 Baker Hughes Incorporated Method of detecting, quantifying and correcting borehole contaminations from multi-frequency, multi-sensitive-volume NMR logging data

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009038744A1 (en) * 2007-09-18 2009-03-26 Baker Hughes Incorporated Nuclear magnetic resonance evaluation using independent component analysis (ica)-based blind source separation
US20140278113A1 (en) * 2013-03-14 2014-09-18 Hamed Chok Real-Time Determination of Formation Fluid Properties Using Density Analysis
US10400595B2 (en) * 2013-03-14 2019-09-03 Weatherford Technology Holdings, Llc Real-time determination of formation fluid properties using density analysis
WO2018039054A1 (en) * 2016-08-23 2018-03-01 Baker Hughes, A Ge Company, Llc Simultaneous inversion of nmr multiple echo trains and conventional logs
US10209391B2 (en) 2016-08-23 2019-02-19 Baker Hughes, A Ge Company, Llc Simultaneous inversion of NMR multiple echo trains and conventional logs
US10605952B2 (en) 2016-08-23 2020-03-31 Baker Hughes, A Ge Company, Llc Simultaneous inversion of NMR multiple echo trains and conventional logs
US20180308502A1 (en) * 2017-04-20 2018-10-25 Thomson Licensing Method for processing an input signal and corresponding electronic device, non-transitory computer readable program product and computer readable storage medium
CN111980663A (en) * 2020-07-21 2020-11-24 中海油田服务股份有限公司 Multi-frequency multi-dimensional nuclear magnetic logging method and device
CN111980663B (en) * 2020-07-21 2023-08-15 中海油田服务股份有限公司 Multi-frequency multi-dimensional nuclear magnetic logging method and device

Also Published As

Publication number Publication date
GB2422198A (en) 2006-07-19
US7812602B2 (en) 2010-10-12
GB2422198B (en) 2009-05-27
US20070290684A1 (en) 2007-12-20
GB0526441D0 (en) 2006-02-08

Similar Documents

Publication Publication Date Title
US7812602B2 (en) Multiple echo train inversion
US5680043A (en) Nuclear magnetic resonance technique for determining gas effect with borehole logging tools
US7199580B2 (en) System and methods for T1-based logging
Freedman Advances in NMR logging
US6972564B2 (en) Objective oriented methods for NMR log acquisitions for estimating earth formation and fluid properties
US6344744B2 (en) Multiple frequency method for nuclear magnetic resonance longitudinal relaxation measurement and pulsing sequence for power use optimization
US6331775B1 (en) Gas zone evaluation by combining dual wait time NMR data with density data
US7804297B2 (en) Methodology for interpretation and analysis of NMR distributions
US8330460B2 (en) Method and apparatus for determining multiscale similarity between NMR measurements and a reference well log
US6703832B2 (en) Method for detecting hydrocarbons by comparing NMR response at different depths of investigation
US6023163A (en) Well logging method and apparatus for determining gas and diffusion coefficient using NMR
AU2126699A (en) A method for determining a characteristic of a gas-bearing formation traversed by a borehole
WO2001042817A1 (en) Nuclear magnetic resonance method and logging apparatus
CA2563698A1 (en) Use of measurements made in one echo train to correct ringing in second to avoid use of phase alternated pair in the second
US8781745B2 (en) NMR-DNA fingerprint
AU2009215414B2 (en) Echo-decay-acceleration data acquisition method for gas identification using a low-field gradient
Hürlimann et al. NMR well logging
US6452389B1 (en) NMR pulse sequences for increasing the efficiency of acquisition
CN1474199A (en) Method for detecing hydrocarbon compund from NMR data
US7336071B2 (en) Enhancement of NMR vertical resolution using walsh function based inversion
GB2434875A (en) T1 based logging

Legal Events

Date Code Title Description
AS Assignment

Owner name: BAKER HUGHES INCORPORATED, TEXAS

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:EDWARDS, CARL M.;REEL/FRAME:016197/0763

Effective date: 20050113

STCB Information on status: application discontinuation

Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION