NMR SPIN-ECHO AMPLITUDE ESTIMATION
20170219675 · 2017-08-03
Inventors
Cpc classification
G01R33/583
PHYSICS
G01R33/4625
PHYSICS
G01R33/448
PHYSICS
G01R33/441
PHYSICS
International classification
G01R33/58
PHYSICS
Abstract
Technologies applicable to NMR spin-echo amplitude estimation are disclosed. Example methods may calibrate for distortion of a shape and estimated amplitude of measured NMR spin or gradient echoes. NMR spin or gradient echo measurements may be performed on a sample. The measured NMR spin or gradient echoes may be used to calculate an echo-shape calibration factor. The echo-shape calibration factor may estimate an effect of echo shape on estimated spin or gradient echo amplitude(s) of the NMR spin or gradient echoes. The echo-shape calibration factor may be used to correct for underestimation or overestimation of the spin or gradient echo amplitude(s).
Claims
1. A method for calibrating for distortion of a shape and estimated amplitude of one or more recorded Nuclear Magnetic Resonance (NMR) spin or gradient echoes, the distortion caused by changes in a magnetic field gradient across or within a sample of investigation, comprising: performing one or more NMR spin or gradient echo measurements on the sample in order to record one or more NMR spin or gradient echoes, wherein the one or more NMR spin or gradient echo measurements resolve a time domain shape of the one or more recorded NMR spin or gradient echoes; using the one or more recorded NMR spin or gradient echoes to calculate an echo-shape calibration factor that estimates an effect of echo shape on estimated spin or gradient echo amplitude(s) of the one or more recorded NMR spin or gradient echoes; and using the echo-shape calibration factor to correct for underestimation or overestimation of the spin or gradient echo amplitude(s).
2. The method of claim 1, wherein the echo-shape calibration factor is calculated as a ratio of a peak measured or estimated amplitude of a spin or gradient echo in the one or more NMR spin or gradient echo measurements, divided by a time-windowed average of the spin or gradient echo within a time-averaging window.
3. The method of claim 1, wherein using the echo-shape calibration factor to correct for the underestimation or overestimation of the spin or gradient echo amplitude(s) comprises multiplying the spin or gradient echo amplitude(s) by the echo-shape calibration factor.
4. The method of claim 1, further comprising performing one or more background NMR spin or gradient echo measurements on a known sample under a background magnetic field gradient condition to resolve a time domain shape of one or more recorded background NMR spin or gradient echoes.
5. The method of claim 4, further comprising normalizing the one or more recorded background NMR spin or gradient echoes to have a first unit peak amplitude, and normalizing the one or more recorded NMR spin or gradient echoes, recorded pursuant to the one or more NMR spin or gradient echo measurements on the sample, to have a second unit peak amplitude.
6. The method of claim 5, wherein using the one or more recorded NMR spin or gradient echoes to calculate the echo-shape calibration factor comprises calculating the echo-shape calibration factor as a ratio of a time-windowed average of the one or more normalized background NMR spin or gradient echo measurements, divided by a time-windowed average of the one or more normalized NMR spin or gradient echo measurements.
7. The method of claim 1, further comprising applying the echo-shape calibration factor to the one or more recorded NMR spin or gradient echoes to correct for errors in estimated echo amplitudes due to unknown external or internal magnetic field gradients across the sample.
8. The method of claim 1, further comprising estimating spin or gradient echo shape(s) in the one or more recorded NMR spin or gradient echoes at least in part by averaging multiple time-offset spin or gradient echoes in order to reduce noise.
9. The method of claim 1, further comprising correcting for NMR spin or gradient shape variations on fluid content measurements of fluid-saturated samples investigated using an NMR spectrometer.
10. The method of claim 1, wherein the one or more recorded NMR spin or gradient echoes are induced and measured between NMR measurements of bulk fluids and NMR measurements of fluid-saturated porous materials, and wherein using the one or more recorded NMR spin or gradient echoes to calculate the echo-shape calibration factor uses different spin or gradient echo shapes induced and measured between the NMR measurements of bulk fluids and the NMR measurements of fluid-saturated porous materials.
11. The method of claim 1, wherein using the one or more recorded NMR spin or gradient echoes to calculate the echo-shape calibration factor comprises calculating an effect of magnetically susceptible earth formations on NMR estimates of volumetric fluid content within said earth formations.
12. The method of claim 1, wherein the sample comprises an earth core sample or a soil sample, wherein the one or more NMR spin or gradient echo measurements are performed with a NMR spectrometer, and further comprising performing one or more background calibration measurements in the NMR spectrometer, using samples of bulk fluids.
13. The method of claim 1, wherein the one or more NMR spin or gradient echo measurements are performed within an earth borehole with an NMR logging tool, and further comprising performing one or more background calibration measurements using the NMR logging tool or a similar NMR logging tool to measure an NMR response of bulk fluids.
14. The method of claim 1, wherein the one or more NMR spin or gradient echo measurements are performed with an earth's field surface NMR instrument, and wherein the spin or gradient echo amplitude(s) corrected using the echo shape calibration factor comprise spin or gradient echo amplitude(s) used for water content and relaxation time estimates.
15. The method of claim 1, further comprising using the one or more recorded NMR spin or gradient echoes or the echo shape calibration factor to determine magnetic susceptibility effects in the sample.
16. The method of claim 15, wherein the one or more NMR spin or gradient echo measurements are performed with a medical Magnetic Resonance Imaging (MM) device.
17. The method of claim 15, wherein the one or more NMR spin or gradient echo measurements are performed within an earth borehole with an NMR logging tool, or at an earth surface with an earth's field surface NMR instrument.
18. The method of claim 15, wherein the one or more NMR spin or gradient echo measurements are performed with a NMR spectrometer.
19. The method of claim 15, wherein the one or more NMR spin or gradient echo measurements comprise a Nuclear Quadrupole Resonance (NQR) measurement.
20. The method of claim 1, wherein the one or more NMR spin or gradient echo measurements comprise a Nuclear Quadrupole Resonance (NQR) measurement.
21. An apparatus adapted to calibrate for distortions of shapes and estimated amplitudes of recorded NMR spin or gradient echoes, the distortion caused by changes in a magnetic field gradient across or within a sample of investigation, the apparatus comprising: a measurement control unit configured to perform one or more NMR spin or gradient echo measurements on the sample in order to record the one or more NMR spin or gradient echoes, wherein the one or more NMR spin or gradient echo measurements resolve a time domain shape of the one or more recorded NMR spin or gradient echoes; and a data processing unit adapted to: use the one or more recorded NMR spin or gradient echoes to calculate an echo-shape calibration factor that estimates an effect of echo shape on estimated spin or gradient echo amplitude(s) of the one or more recorded NMR spin or gradient echoes; and; and use the echo-shape calibration factor to correct for underestimation or overestimation of the spin or gradient echo amplitude(s).
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0009]
[0010]
[0011]
DETAILED DESCRIPTION
[0012] Prior to explaining embodiments of the invention in detail, it is to be understood that the invention is not limited to the details of construction or arrangements of the components and method steps set forth in the following description or illustrated in the drawings. The invention is capable of other embodiments and of being practiced and carried out in various ways. Also, it is to be understood that the phraseology and terminology employed herein are for the purpose of the description and should not be regarded as limiting.
[0013] Technologies directed to NMR spin-echo amplitude estimation are disclosed. In one embodiment, an echo-shape correction factor may be estimated as a ratio of a peak measured or estimated amplitude of an echo in an NMR sample measurement, divided by a time-windowed average of the echo within a time-averaging window. Spin echo amplitudes estimated by time-domain averaging may then be multiplied by this correction factor to correct for the effect of the echo shape.
[0014] In another embodiment, an optional background measurement of the spin echo shape may be normalized to have unit peak amplitude, and a sample measurement of the spin echo shape may also be normalized to have unit peak amplitude. The echo shape correction factor may then be calculated as a ratio of the time-windowed average of the normalized background measurement, divided by the time-windowed average of the normalized sample measurement. The echo-shape calibration factor may then be applied to the sample NMR measurement to correct for errors in the estimated echo amplitudes due to unknown external or internal magnetic field gradients across the sample. This embodiment enables the background measurement to be used as a measurement calibration for spin density under known conditions (i.e. for calibrating the NMR amplitude response for 100% volumetric bulk water), while the echo-shape calibration procedure corrects for differences in the time-windowed-averages of echo amplitudes that are caused by differences in the shapes of the echoes under the background and sample measurement conditions.
[0015] In another embodiment, the echo shape(s) in the sample and optional background conditions may be estimated by averaging multiple time-offset spin echoes in order to reduce the noise on the estimate of the echo shape function. For example, in a Carr-Purcell-Meiboom-Gill (CPMG) measurement that acquires 10 or more echoes in a single scan, the average echo shape can be estimated by averaging the first 10 echoes together, to provide a composite “mean echo shape” with higher signal to noise than the first echo alone.
[0016] In another embodiment, the disclosed method may be used to correct for NMR echo shape variations on fluid content measurements of fluid-saturated samples investigated using an NMR spectrometer. In this embodiment, the disclosed method may be used to calibrate for the different echo shapes that are induced and measured between NMR measurements of bulk fluids and NMR measurements of fluid-saturated porous materials.
[0017] In another embodiment, the disclosed method may be used to calibrate for the effect of magnetically susceptible earth formations on NMR estimates of volumetric fluid content within said earth formations. In one such embodiment, the NMR sample measurements are performed on earth core or soil samples in a NMR spectrometer, and optional background calibration measurements are performed in the same NMR spectrometer using samples of bulk fluids such as water or oil. In another such embodiment, the NMR sample measurements are performed within earth boreholes using an NMR logging tool, and optional background calibration measurements are performed using the same NMR or similar logging tool to measure the NMR response of bulk fluids such as water or oil. In another such embodiment, the sample NMR measurements are performed using an earth's field surface NMR instrument, and the echo shape calibration factor(s) are used to correct the estimates of spin echo amplitudes that are used for water content and relaxation time estimates.
[0018] In another embodiment, the disclosed method may be used to correct for magnetic susceptibility effects in medical MRI images.
[0019] In some embodiments, any of the various methods disclosed herein may be carried out by a computing device. Computing devices comprising appropriate hardware, and computer readable media comprising appropriate instructions to carry out the methods described herein are expressly contemplated. This includes, e.g., NMR measurement acquisition hardware, signal processing software, and any subsequent data processing software.
[0020]
[0021] In
[0022] In general, with regard to
[0023] NMR measurement apparatus 700 may also be configured to receive and record NMR signal data received via the NMR sensor(s) 750. NMR measurement apparatus 700 may be configured to receive and record NMR signal data after one or more excitation pulses. In some embodiments, the computer 710 may be configured to close the receive switch(es) 760 after a pulse. The preamplifier(s) 770 amplify NMR signals received via induction coil(s) 750. The AD converter(s) 720 convert the received and amplified signals to digital NMR signal data, e.g. by sampling received NMR signals at a desired sampling rate, and the computer 710 or other device equipped with storage media may be configured to store the resulting digital NMR signal data.
[0024] In some embodiments, the data processing module(s) 702 may be configured to process NMR measurement data from sample volume 780, generated by operation of the measurement control module 701 and the various other components of NMR measurement apparatus 700. Data processing module(s) 702 may be configured to use measured NMR properties to determine and/or facilitate determination of any of various information used in connection with methods described herein. It will be appreciated that while the computer 710 may be configured to include data processing module(s) 702, in some embodiments NMR measurements and data processing may be performed separately, e.g., by first performing measurements with system 700, then processing acquired NMR data at a later time and/or with a different computing device or by a human operator.
[0025] It will be appreciated that NMR measurement apparatus may be configured differently than illustrated in
[0026] In some embodiments, NMR sensor(s) 750 may comprise an array of coils comprising one or more transmit coils, one or more receive coils, and/or one or more combination transmit and receive coils. For example, NMR sensor(s) 750 may comprise one transmit coil and multiple receive coils. NMR sensor(s) 750 may comprise one combination transmit and receive coil, and multiple receive coils. NMR sensor(s) 750 may comprise multiple combination transmit and receive coils. These and other multicoil arrangements may be configured in some embodiments as will be appreciated. Multicoil arrangements are useful for localization of fluids in sample volume 780, as described for example in U.S. Pat. No. 7,466,128, entitled “Multicoil Data Acquisition and Processing Methods,” issued Dec. 16, 2008, which is incorporated by reference herein.
[0027] Any combination of hardware and software that enables the acquisition and processing of NMR signals is suitable to implement embodiments of this disclosure. An architecture to implement the disclosed methods could comprise, for example, elements illustrated in
[0028] In general, NMR measurements may be collected by transmitting one or more pulses of alternating current through NMR sensor(s) 750. The alternating current may be tuned to the Larmor frequency of hydrogen nuclei, for example, and may generate a magnetic field in the sample volume 780 alternating at the Larmor frequency. The alternating magnetic field radiates into the sample volume 780 and modifies the nuclear magnetization state of hydrogen atoms present in fluids in sample volume 780. The transmitted alternating magnetic field perturbs the magnetization from equilibrium alignment in a static magnetic field, so that some component of the nuclear magnetization rotates into the transverse “xy” plane. Once rotated from equilibrium, the magnetization relaxes over time back to the equilibrium state over time, decaying from the transverse plane and re-growing along the longitudinal axis. The rotation of the magnetization by the transmitted pulse(s) and subsequent relaxation to equilibrium are described by the phenomenological Bloch equations. The evolution of the magnetization under the Bloch equations depends on several variables including the amplitude of the transmitted field, the duration and timing of the transmitted field, the phase of the transmitted field, the longitudinal relaxation time T1, FID relaxation rate T2*, and/or the spin-spin relaxation time T2 of the hydrogen nuclei under investigation. These aspects of NMR measurement may be used in determining the various NMR properties described herein.
[0029]
[0030] In
[0031] It will be appreciated that the apparatus 700 illustrated in
[0032]
[0033] Methods according to
[0034] In some embodiments, methods according to
[0035] Methods according to
[0036] In some embodiments, spin or gradient echo shape(s) in the one or more recorded NMR spin or gradient echoes may be estimated at least in part by averaging multiple time-offset spin or gradient echoes, in order to reduce noise.
[0037] In some embodiments, e.g., for fluid-saturated samples investigated using an NMR spectrometer, methods may correct for NMR spin or gradient shape variations on fluid content measurements. When the sample comprises an earth core sample or a soil sample, and measurements are performed with a NMR spectrometer, and methods may include performing one or more background calibration measurements in the NMR spectrometer, using samples of bulk fluids.
[0038] In some embodiments, the one or more NMR spin or gradient echo measurements may be performed within an earth borehole with an NMR logging tool, and methods may include performing one or more background calibration measurements using the NMR logging tool (or a similar NMR logging tool) to measure an NMR response of bulk fluids.
[0039] In some embodiments, the one or more NMR spin or gradient echo measurements may be performed with an earth's field surface NMR instrument, and the spin or gradient echo amplitude(s) corrected using the echo shape calibration factor may comprise spin or gradient echo amplitude(s) used for water content and relaxation time estimates. In some embodiments, the one or more NMR spin echo or gradient echo measurements comprise a Nuclear Quadrupole Resonance (NQR) measurement.
[0040] In some embodiments, methods according to
[0041] In some embodiments, the one or more recorded NMR spin or gradient echoes may be induced and measured, e.g., at block 301, between NMR measurements of bulk fluids and NMR measurements of fluid-saturated porous materials. Using the one or more recorded NMR spin or gradient echoes to calculate the echo-shape calibration factor at block 303 may use different spin or gradient echo shapes induced and measured between the NMR measurements of bulk fluids and the NMR measurements of fluid-saturated porous materials.
[0042] In some embodiments, using the one or more recorded NMR spin or gradient echoes may be used to calculate the echo-shape calibration factor at block 303 may comprise calculating an effect of magnetically susceptible earth formations on NMR estimates of volumetric fluid content within said earth formations. In general, the one or more recorded NMR spin or gradient echoes, or the echo shape calibration factor, or a combination thereof, may be used to determine magnetic susceptibility effects in the sample. The measurements may be performed with a medical Magnetic Resonance Imaging (MRI) device, and methods may determine magnetic susceptibility effects in medical samples (e.g., human organs and the like). The measurements may be performed within an earth borehole with an NMR logging tool, or at an earth surface with an earth's field surface NMR instrument, and methods may determine magnetic susceptibility effects in subsurface earth formations. The measurements may be performed within a NMR spectrometer, and methods may determine magnetic susceptibility effects in corresponding spectrometer samples. The measurements may comprise Nuclear Quadrupole Resonance (NQR) measurements, and methods may determine magnetic susceptibility effects in samples measured by NQR.
[0043] There are various approaches by which apparatus and methods processes described herein can be implemented (e.g., hardware, software, and/or firmware), and the preferred approach may vary with the context in which the apparatus and methods are deployed. For example, if an implementer determines that speed and accuracy are paramount for operations of a computer or controller, the implementer may opt for a mainly hardware and/or firmware vehicle; if flexibility is paramount, the implementer may opt for a mainly software implementation; or, yet again alternatively, the implementer may opt for some combination of hardware, software, and/or firmware.
[0044] With respect to the use of substantially any plural and/or singular terms herein, those having skill in the art can translate from the plural to the singular and/or from the singular to the plural as is appropriate to the context and/or application. The various singular/plural permutations may be expressly set forth herein for sake of clarity.
[0045] It will be understood by those within the art that, in general, terms used herein, and especially in the appended claims (e.g., bodies of the appended claims) are generally intended as “open” terms (e.g., the term “including” should be interpreted as “including but not limited to,” the term “having” should be interpreted as “having at least,” the term “includes” should be interpreted as “includes but is not limited to,” etc.). It will be further understood by those within the art that if a specific number of an introduced claim recitation is intended, such an intent will be explicitly recited in the claim, and in the absence of such recitation no such intent is present. For example, as an aid to understanding, the following appended claims may contain usage of the introductory phrases “at least one” and “one or more” to introduce claim recitations. However, the use of such phrases should not be construed to imply that the introduction of a claim recitation by the indefinite articles “a” or “an” limits any particular claim containing such introduced claim recitation to inventions containing only one such recitation, even when the same claim includes the introductory phrases “one or more” or “at least one” and indefinite articles such as “a” or “an” (e.g., “a” and/or “an” should typically be interpreted to mean “at least one” or “one or more”); the same holds true for the use of definite articles used to introduce claim recitations. In addition, even if a specific number of an introduced claim recitation is explicitly recited, those skilled in the art will recognize that such recitation should typically be interpreted to mean at least the recited number (e.g., the bare recitation of “two recitations,” without other modifiers, typically means at least two recitations, or two or more recitations). Furthermore, in those instances where a convention analogous to “at least one of A, B, and C, etc.” is used, in general such a construction is intended in the sense one having skill in the art would understand the convention (e.g., “a system having at least one of A, B, and C” would include but not be limited to systems that have A alone, B alone, C alone, A and B together, A and C together, B and C together, and/or A, B, and C together, etc.). In those instances where a convention analogous to “at least one of A, B, or C, etc.” is used, in general such a construction is intended in the sense one having skill in the art would understand the convention (e.g., “a system having at least one of A, B, or C” would include but not be limited to systems that have A alone, B alone, C alone, A and B together, A and C together, B and C together, and/or A, B, and C together, etc.). It will be further understood by those within the art that virtually any disjunctive word and/or phrase presenting two or more alternative terms, whether in the description, claims, or drawings, should be understood to contemplate the possibilities of including one of the terms, either of the terms, or both terms. For example, the phrase “A or B” will be understood to include the possibilities of “A” or “B” or “A and B.”
[0046] While various embodiments have been disclosed herein, other aspects and embodiments will be apparent to those skilled in art.