US20140074405A1 - Modified pulse sequence to estimate properties - Google Patents
Modified pulse sequence to estimate properties Download PDFInfo
- Publication number
- US20140074405A1 US20140074405A1 US14/083,375 US201314083375A US2014074405A1 US 20140074405 A1 US20140074405 A1 US 20140074405A1 US 201314083375 A US201314083375 A US 201314083375A US 2014074405 A1 US2014074405 A1 US 2014074405A1
- Authority
- US
- United States
- Prior art keywords
- pulse sequence
- property
- properties
- time
- relaxation
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
- G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
- G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/448—Relaxometry, i.e. quantification of relaxation times or spin density
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/14—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electron or nuclear magnetic resonance
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/32—Electric 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
- G01N11/10—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties by moving a body within the material
- G01N11/16—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties by moving a body within the material by measuring damping effect upon oscillatory body
Definitions
- This patent specification relates to estimating properties from exponentially decaying data. More particularly, this patent specification relates to methods and systems for estimating material properties from a non-uniform sequence of excitations, such as with NMR measurements.
- NMR Nuclear Magnetic Resonance
- NMR magnetization is characterized by relaxation time constants T 1 , T 2 or diffusion D and corresponding amplitudes f T 1 (T 1 ), f T 2 (T 2 ) or f D (D). Since the relaxation times and/or diffusion constants are expected to be continuous and typically span several decades, these amplitudes are often referred to as relaxation-time or diffusion distributions. These distributions provide a wealth of information about both the rock as well as fluid properties. See e.g. D. Allen, S. Crary, and R. Freedman, How to use borehole Nuclear Magnetic Resonance, Oilfield Review, pages 34-57, 1997.
- An NMR measurement consists of a sequence of transverse magnetic pulses transmitted by an antenna.
- the CPMG pulse sequence is used to measure T 2 relaxation. It consists of a pulse that tips hydrogen protons 90° and is followed by several thousand pulses that refocus the protons by tipping them 180°. The data are acquired by an antenna between the evenly spaced 180° pulses and are thus uniformly spaced in time.
- the T 1 relaxation is studied using the inversion-recovery pulse sequence. There is no specific rule for choosing the recovery times. See, Y. Q. Song, L. Venkataramanan, and L. Burcaw, Determining the resolution of Laplace inversion spectrum, Journal of Chemical Physics, page 104104, 2005. Similarly, there is no clear mathematical framework to optimally choose the acquisition times for diffusion measurements.
- the data acquired from all of the above pulse sequences can be represented by a Fredholm integral (see, M. D. Hurlimann and L. Venkataramanan, Quantitative measurement of two dimensional distribution functions of diffusion and relaxation in grossly homogeneous fields, Journal Magnetic Resonance, 157:31-42, July 2002 and L. Venkataramanan, Y. Q. Song, and M. D. Hurlimann, Solving Fredholm integrals of the first kind with tensor product structurein 2 and 2.5 dimensions, IEEE Transactions on Signal Processing, 50:1017-1026, 2002, hereinafter “Venkataramanan et al. 2002”),
- variable x typically refers to T 1 , T 2 , or 1/D. See, Venkataramanan et al. 2002.
- the measured data in eqn. (1) is related to a Laplace transform of the unknown underlying function f x (x).
- a typical crude oil has a number of components, such as methane, ethane etc.
- N i refer to the number of carbon atoms or chain length of the i-th component.
- N denote the molar average chain length of the crude oil, also defined as the harmonic mean of the chain lengths,
- f N (N i ) denotes the mass fraction of component i in the hydrocarbon.
- N i the bulk relaxation time of component i in a crude oil follows a power law and depends inversely on its chain length or number of carbon atoms, N i .
- the relaxation time also depends inversely on the average chain length N of the hydrocarbon, according to
- parameters B, ⁇ and ⁇ are constants at a given pressure and temperature.
- D i refers to the diffusion coefficient of component i
- a method for estimating a property of a material is provided.
- a measurable change in the material is induced and measurements are made of an aspect of the material that tends to exponentially decay following the induction.
- the induction is repeated at predetermined intervals, and the material property is estimated based on measurements of the aspect.
- the intervals are designed so as to improve the estimation of the property.
- the induction can be a transmitted transverse magnetic pulse as with NMR measurements, according to some embodiments.
- the induction is the induction of a pressure or temperature change in the material.
- the interval spacing e.g. in time, in frequency, etc.
- the interval spacing in time follows a power law and allows for the material property to be directly computed from the measurements.
- the method can be carried out as part of an oilfield service, and the measurements can be made downhole.
- the material can be a fluid, a solid or a mixture.
- the material is one or more rock solids and the material properties being estimated are petrophysical properties of the rock solids.
- a system for estimating a property of a material includes an inducer that is adapted to induce a measurable change in the material, and one or more sensors adapted to measure an aspect of the material that tends to exponentially decay following an induction.
- a controller is coupled to the inducer, and adapted and programmed to repeat the induction at predetermined intervals.
- a processing system is adapted to estimate the material property based on measurements. The induction intervals are designed so as to improve the estimation of the property.
- FIG. 1 is a plot showing the values of n and ⁇ for different values of ⁇ as determined by equation (8) herein;
- FIG. 2 is a plot showing an example of a single exponential acquired at uniform time intervals in the time-domain
- FIG. 3 is a plot showing a non-uniform sampling, according to some embodiments.
- FIG. 4 plots examples of non-uniform sampling in time domain that can directly provide moments of relaxation or diffusion distribution, according to some embodiments
- FIG. 5 is a plot showing the ratio of the error-bar in the moments obtained from uniform sampling to non-uniform sampling
- FIG. 6A shows an NMR tool being deployed in a wellbore, according to embodiments
- FIG. 6B is a schematic of a processing system used to control, record and process data from an NMR tool.
- FIG. 7 illustrates another example of a wellsite system in which an NMR tool can be employed, according to some embodiments.
- a process may correspond to a method, a function, a procedure, a subroutine, a subprogram, etc. When a process corresponds to a function, its termination corresponds to a return of the function to the calling function or the main function.
- embodiments of the invention may be implemented, at least in part, either manually or automatically.
- Manual or automatic implementations may be executed, or at least assisted, through the use of machines, hardware, software, firmware, middleware, microcode, hardware description languages, or any combination thereof.
- the program code or code segments to perform the necessary tasks may be stored in a machine readable medium.
- a processor(s) may perform the necessary tasks.
- modified pulse sequences allow for the direct measurements of moments of relaxation-time or diffusion distributions.
- the data from these pulse sequences provides direct estimates of fluid properties such as average chain length and viscosity of a hydrocarbon. In comparison to the uniformly-spaced pulse sequence, these pulse sequences are faster and have a lower error bar in computing the fluid properties.
- variable x corresponds to relaxation times T 1 , T 2 , or diffusion D.
- the moments of relaxation times or diffusion can be directly computed from the measured magnetization data
- ⁇ ( ) represents the Gamma function
- variable ⁇ is in two parts: a real number ⁇ and an integer n where the mathematical operator t ⁇ -1 operates on the n-th derivative of the data. Given a value of ⁇ , the values of n and ⁇ are unique and given as,
- FIG. 1 is a plot 110 showing the values of n and ⁇ for different values of ⁇ as determined by above these equations.
- the computation of moments of relaxation times T 1 , T 2 , or diffusion D, in conjunction with the scaling laws, can lead to greatly improved alternative pulse sequences to specifically target or estimate only certain moments of relaxation time/diffusion.
- targeted measurements are described as measurements that target a certain specified parameter at a particular location in the reservoir space. These targeted measurements are achieved by both hardware and software modifications.
- the concept of targeted measurements is applied to NMR data and alternate pulse sequences are described to directly estimate the moments of relaxation time or diffusion coefficients.
- alternate pulse sequences are designed by taking advantage of the specific exponential-kernel structure in eqn. (1) as well as the computation of moments using eqn. (7).
- the magnetization data from the new pulse sequences can be directly interpreted to obtain fluid properties.
- the non-uniformly spaced pulse sequences are targeted to specifically measure one of the moments of the relaxation time T 1 , T 2 , or D distribution.
- the random variable in eqn. (1) is T 1 and a specific ⁇ -th moment of T 1 is desired, where ⁇ >0. This corresponds to the right-hand side of the ⁇ axis in FIG. 1 .
- eqns. (7) and (8) From eqns. (7) and (8),
- M ⁇ ( ⁇ ) is the magnetization data in the ⁇ domain.
- Eqn. (11) indicates that if the data are acquired such that they are uniformly spaced in the ⁇ domain, then the sum of the measured data directly provides the ⁇ -th moment. From eqn. (10), uniform spacing in the ⁇ domain implies non-uniform spacing in the time (t) domain.
- FIG. 2 is a plot showing an example of a single exponential acquired at uniform time intervals in the time-domain.
- FIG. 3 is a plot showing a non-uniform sampling, according to some embodiments. In particular, using eqn. (10), the data shown by curve 310 are sampled uniformly in the ⁇ square root over (t) ⁇ domain.
- FIG. 4 plots examples of non-uniform sampling in time domain that can directly provide moments of relaxation or diffusion distribution, according to some embodiments.
- plots 410 , 412 , 414 and 416 represent sampling for ⁇ values of 1, 0.75, 0.5, and 0.25 respectively. When 0 ⁇ 1, this leads to finer sampling at the initial part of the data and coarser sampling at the tail end of the data.
- the area under the integrand [M T ⁇ M(0)] can be used to directly estimate the ⁇ -th moment.
- the concept of the pulse sequence for ⁇ >0 and ⁇ 1 ⁇ 0 is valid for other ranges of ⁇ as well.
- pulse sequences that directly measure the ⁇ -th moment are obtained by Taylor series expansion of the data with second-order terms in eqn. (13). Higher-order terms in the Taylor series lead to further negative values of ⁇ .
- T 1 ⁇ can be obtained by sampling uniformly in the ⁇ domain, where
- this modified pulse sequence also enables fluid properties such as the average chain length and viscosity to be directly estimated in a simple manner. Re-writing eqns. (3) and (4), the average chain length and viscosity can be computed from moments of T 1 or T 2 relaxation time or diffusion, as
- N _ [ B - 1 / ⁇ ⁇ ⁇ T 1 1 / ⁇ ⁇ ] - ⁇ ⁇ + ⁇ ( 20 )
- N _ [ A - 1 / v ⁇ ⁇ D 1 / v ⁇ ] - v v + ⁇ ( 21 ) ⁇ ⁇ ⁇ D - 1 ⁇ ⁇ ⁇ ⁇ [ ⁇ T 1 1 ⁇ ⁇ ] ⁇ ⁇ ⁇ - ⁇ ⁇ ⁇ v ⁇ + ⁇ ⁇ ⁇ T 1 - v ⁇ ⁇ ( 22 )
- the average chain length can be estimated from a modified pulse sequence to directly estimate the (1/ ⁇ )-th moment of relaxation time or (1/ ⁇ )-th moment of diffusion.
- viscosity can be estimated from two modified pulse sequences to directly estimate the (1/ ⁇ )-th and (1/ ⁇ )-th moment of relaxation time.
- the kernel in eqn. (1) is e ⁇ Dt .
- these pulse sequences are optimized to estimate any desired moments of relaxation time and/or diffusion. In comparison to pulse sequences where the data is uniformly acquired in time, these pulse sequences can be used to acquire the data faster. For example, given ⁇ >0, the error-bar in the ⁇ -th moment of T 1 relaxation is,
- N e denotes the total number of echoes obtained from uniform sampling between a given lower and upper limit of relaxation time and t E corresponds to the spacing between the echoes.
- t E corresponds to the spacing between the echoes.
- a modified pulse sequence can directly provide moments of relaxation-time or diffusion distributions.
- This pulse sequence can be adapted to the desired moment of relaxation-time or diffusion coefficient.
- the data from this pulse sequence provides direct estimates of fluid properties such as average chain length and viscosity of a hydrocarbon. In general, in comparison with data acquired uniformly in time, this pulse sequence is faster and has a lower error bar in computing the fluid properties.
- the first term in eqn. (26) is zero.
- the first term in eqn. (27) is valid only for negative integer values of ⁇ and corresponds to the
- the first term in eqn. (27) is zero.
- FIG. 6A shows an NMR tool being deployed in a wellbore, according to embodiments.
- Wireline truck 610 is deploying wireline cable 612 into well 630 via well head 620 .
- Wireline tool 640 is disposed on the end of the cable 612 in a subterranean formation 600 .
- Tool 640 includes an NMR tool, such as Schlumberger's Combinable Magnetic Resonance (CMR) Tool, or Schlumberger's MR Scanner Tool.
- CMR Combinable Magnetic Resonance
- the NMR tool is combined with other downhole tools such as Schlumberger's Modular Formation Dynamics Tester downhole sampling tool.
- the NMR tool uses a pulse sequence as described herein.
- Data from the NMR tool can be recorded and/or processed downhole within tool 640 , and/or can be transmitted to wireline truck 610 for recording and/or processing.
- the NMR tool can be controlled, including the selection of pulse sequences, locally with processing systems within tool 640 and/or from the surface using processing systems for example in wireline truck 610 .
- FIG. 6B is a schematic of a processing system used to control, record and process data from an NMR tool.
- Processing system 650 includes one or more central processing units 674 , storage system 672 , communications and input/output modules 670 , a user display 676 and a user input system 678 .
- Input/output modules 670 include modules to communicate with and control the NMR tool such that modified pulse sequences can be used as described herein.
- FIG. 7 illustrates another example of a wellsite system in which an NMR tool can be employed, according to some embodiments.
- the wellsite can be onshore or offshore.
- a borehole 711 is formed in subsurface formations by rotary drilling in a manner that is well known.
- Embodiments of the invention can also use directional drilling, as will be described hereinafter.
- a drill string 712 is suspended within the borehole 711 and has a bottom hole assembly 700 which includes a drill bit 705 at its lower end.
- the surface system includes platform and derrick assembly 710 positioned over the borehole 711 , the assembly 710 including a rotary table 716 , kelly 717 , hook 718 and rotary swivel 719 .
- the drill string 712 is rotated by the rotary table 716 , energized by means not shown, which engages the kelly 717 at the upper end of the drill string.
- the drill string 712 is suspended from a hook 718 , attached to a traveling block (also not shown), through the kelly 717 and a rotary swivel 719 which permits rotation of the drill string relative to the hook.
- a top drive system could alternatively be used.
- the surface system further includes drilling fluid or mud 726 stored in a pit 727 formed at the well site.
- a pump 729 delivers the drilling fluid 726 to the interior of the drill string 712 via a port in the swivel 719 , causing the drilling fluid to flow downwardly through the drill string 712 as indicated by the directional arrow 708 .
- the drilling fluid exits the drill string 712 via ports in the drill bit 705 , and then circulates upwardly through the annulus region between the outside of the drill string and the wall of the borehole, as indicated by the directional arrows 709 .
- the drilling fluid lubricates the drill bit 705 and carries formation cuttings up to the surface as it is returned to the pit 727 for recirculation.
- the bottom hole assembly 700 of the illustrated embodiment a logging-while-drilling (LWD) module 720 , a measuring-while-drilling (MWD) module 730 , a roto-steerable system and motor, and drill bit 705 .
- LWD logging-while-drilling
- MWD measuring-while-drilling
- the LWD module 720 is housed in a special type of drill collar, as is known in the art, and can contain one or a plurality of known types of logging tools. It will also be understood that more than one LWD and/or MWD module can be employed, e.g. as represented at 720 A. (References, throughout, to a module at the position of 720 can alternatively mean a module at the position of 720 A as well.)
- the LWD module includes capabilities for measuring, processing, and storing information, as well as for communicating with the surface equipment.
- the LWD module includes an NMR tool such as Schlumberger's ProVISION NMR tool.
- the MWD module 730 is also housed in a special type of drill collar, as is known in the art, and can contain one or more devices for measuring characteristics of the drill string and drill bit.
- the MWD tool further includes an apparatus (not shown) for generating electrical power to the downhole system. This may typically include a mud turbine generator powered by the flow of the drilling fluid, it being understood that other power and/or battery systems may be employed.
- the MWD module includes one or more of the following types of measuring devices: a weight-on-bit measuring device, a torque measuring device, a vibration measuring device, a shock measuring device, a stick slip measuring device, a direction measuring device, and an inclination measuring device.
- the techniques described herein are applied to estimate properties from other types of excitations than NMR.
- the techniques can be applied to other types of data where excitations induce an exponential decay.
- the techniques are applied to inducing a non-uniform sequence of pressure changes in a material.
- the techniques are applied to induce a non-uniform sequence of temperature changes in a material.
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- General Physics & Mathematics (AREA)
- Geology (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Remote Sensing (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Analytical Chemistry (AREA)
- Geochemistry & Mineralogy (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Geophysics And Detection Of Objects (AREA)
- Earth Drilling (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
Methods and related systems are described for estimating fluid or rock properties from NMR measurements. A modified pulse sequence is provided that can directly provide moments of relaxation-time or diffusion distributions. This pulse sequence can be adapted to the desired moment of relaxation-time or diffusion coefficient. The data from this pulse sequence provides direct estimates of fluid properties such as average chain length and viscosity of a hydrocarbon. In comparison to the uniformly-spaced pulse sequence, these pulse sequences are faster and have a lower error bar in computing the fluid properties.
Description
- This application is a continuation application of U.S. patent application Ser. No. 12/717,512, filed on Mar. 4, 2010, entitled, “MODIFIED PULSE SEQUENCE TO ESTIMATE PROPERTIES,” which issues as U.S. Pat. No. 8,587,302 on Nov. 19, 2013. Both of these are incorporated by reference herein in their entirety.
- 1. Field
- This patent specification relates to estimating properties from exponentially decaying data. More particularly, this patent specification relates to methods and systems for estimating material properties from a non-uniform sequence of excitations, such as with NMR measurements.
- 2. Background
- Traditionally, fluid properties using Nuclear Magnetic Resonance (NMR) are obtained by study of longitudinal or transverse relaxation or diffusion measurements. Longitudinal (T1) relaxation data are acquired using an inversion-recovery pulse sequence where the recovery times are often chosen arbitrarily. Transverse (T2) relaxation data are acquired using the CPMG pulse sequence with pulses that are typically uniformly spaced in the time-domain. There is no clear mathematical methodology for choosing the acquisition times at which the diffusion measurements are obtained. The relaxation or diffusion data thus obtained from these pulse sequences are used to infer the relaxation time or diffusion distributions. In turn, these distributions are used to infer petrophysical or hydrocarbon properties.
- Traditional well-logging using NMR is employed to infer petro-physical and fluid properties of the formation or in-situ hydrocarbon. The multi-exponential time-decay of NMR magnetization is characterized by relaxation time constants T1, T2 or diffusion D and corresponding amplitudes fT
1 (T1), fT2 (T2) or fD(D). Since the relaxation times and/or diffusion constants are expected to be continuous and typically span several decades, these amplitudes are often referred to as relaxation-time or diffusion distributions. These distributions provide a wealth of information about both the rock as well as fluid properties. See e.g. D. Allen, S. Crary, and R. Freedman, How to use borehole Nuclear Magnetic Resonance, Oilfield Review, pages 34-57, 1997. - An NMR measurement consists of a sequence of transverse magnetic pulses transmitted by an antenna. Typically, the CPMG pulse sequence is used to measure T2 relaxation. It consists of a pulse that tips hydrogen protons 90° and is followed by several thousand pulses that refocus the protons by tipping them 180°. The data are acquired by an antenna between the evenly spaced 180° pulses and are thus uniformly spaced in time. On the other hand, the T1 relaxation is studied using the inversion-recovery pulse sequence. There is no specific rule for choosing the recovery times. See, Y. Q. Song, L. Venkataramanan, and L. Burcaw, Determining the resolution of Laplace inversion spectrum, Journal of Chemical Physics, page 104104, 2005. Similarly, there is no clear mathematical framework to optimally choose the acquisition times for diffusion measurements.
- The data acquired from all of the above pulse sequences can be represented by a Fredholm integral (see, M. D. Hurlimann and L. Venkataramanan, Quantitative measurement of two dimensional distribution functions of diffusion and relaxation in grossly homogeneous fields, Journal Magnetic Resonance, 157:31-42, July 2002 and L. Venkataramanan, Y. Q. Song, and M. D. Hurlimann, Solving Fredholm integrals of the first kind with
tensor product structurein 2 and 2.5 dimensions, IEEE Transactions on Signal Processing, 50:1017-1026, 2002, hereinafter “Venkataramanan et al. 2002”), -
- where the variable x typically refers to T1, T2, or 1/D. See, Venkataramanan et al. 2002. Thus the measured data in eqn. (1) is related to a Laplace transform of the unknown underlying function fx(x).
- Scaling laws established on mixtures of alkanes provide a relation between fluid composition and relaxation time and/or diffusion coefficients of the components of the mixture (see D. E. Freed, Scaling laws for diffusion coefficients in mixtures of alkanes, Physical Review Letters, 94:067602, 2005, and D. E. Freed, Dependence on chain length of NMR relaxation times in mixtures of alkanes, Journal of Chemical Physics, 126:174502, 2007). A typical crude oil has a number of components, such as methane, ethane etc. Let Ni refer to the number of carbon atoms or chain length of the i-th component. Let
N denote the molar average chain length of the crude oil, also defined as the harmonic mean of the chain lengths, -
- where fN(Ni) denotes the mass fraction of component i in the hydrocarbon. According to the scaling law, the bulk relaxation time T2,i of component i in a crude oil follows a power law and depends inversely on its chain length or number of carbon atoms, Ni. Thus, the larger the molecule, the smaller the relaxation time T2,i. The relaxation time also depends inversely on the average chain length
N of the hydrocarbon, according to -
T 2,i =BN i −κN −γ (3) - In eqn. (3), parameters B, γ and κ are constants at a given pressure and temperature. Similarly, if Di refers to the diffusion coefficient of component i,
-
D i =AN i −νN −β (4) - where A, ν and β and are constants at a given pressure and temperature. The scaling theory can also be used to derive an expression for fluid viscosity,
- Traditional analysis of the measured NMR data to obtain fluid properties goes as follows. The inverse Laplace transform of eqn. (1) is performed to estimate the probability density function fx(x) from the measured data. Next, using eqns. (3)-(5), the average chain length
N and the chain length distribution fN(N) are computed. Viscosity η of the hydrocarbon is then estimated using eqn. (5). This traditional analysis has some disadvantages. First, the inverse Laplace transform of the measured NMR data in eqn. (1) is non-linear due to the non-negativity constraint on fx(x). See, Venkataramanan et al. 2002, and E. J. Fordham, A. Sezginer, and L. D. Hall, Imaging multiexponential relaxation in the (y, log e T 1) plane, with application to clay filtration in rock cores, J. Mag. Resonance A, 113:139-150, 1995. Next, it is also mathematically ill-conditioned. Often, regularization or prior information about the expected solution is incorporated into the problem formulation to make it better conditioned. However, the choice of the regularization functional as well as the weight given to the prior information is a well-known drawback of this transform due to its non-uniqueness. See, W. H. Press, S. A. Teukolsky, and W. T. Vetterling, Numerical Recipes in C, Cambridge University Press, 1992. Thus errors in the estimation of T1, T2, or D distributions are propagated into errors in the estimation of fluid properties. - According to some embodiments, a method for estimating a property of a material is provided. A measurable change in the material is induced and measurements are made of an aspect of the material that tends to exponentially decay following the induction. The induction is repeated at predetermined intervals, and the material property is estimated based on measurements of the aspect. The intervals are designed so as to improve the estimation of the property.
- The induction can be a transmitted transverse magnetic pulse as with NMR measurements, according to some embodiments. According to some embodiments the induction is the induction of a pressure or temperature change in the material. According to some embodiments, the interval spacing (e.g. in time, in frequency, etc.) is designed according to a formula determined by the property being estimated, and the spacing is non-uniform in both linear and logarithmic scales. According to some embodiments the interval spacing in time follows a power law and allows for the material property to be directly computed from the measurements.
- The method can be carried out as part of an oilfield service, and the measurements can be made downhole. The material can be a fluid, a solid or a mixture. According to some embodiments, the material is one or more rock solids and the material properties being estimated are petrophysical properties of the rock solids.
- According to some embodiments, a system for estimating a property of a material is provided that includes an inducer that is adapted to induce a measurable change in the material, and one or more sensors adapted to measure an aspect of the material that tends to exponentially decay following an induction. A controller is coupled to the inducer, and adapted and programmed to repeat the induction at predetermined intervals. A processing system is adapted to estimate the material property based on measurements. The induction intervals are designed so as to improve the estimation of the property.
- Further features and advantages will become more readily apparent from the following detailed description when taken in conjunction with the accompanying drawings.
- The present disclosure is further described in the detailed description which follows, in reference to the noted plurality of drawings by way of non-limiting examples of exemplary embodiments, in which like reference numerals represent similar parts throughout the several views of the drawings, and wherein:
-
FIG. 1 is a plot showing the values of n and μ for different values of ω as determined by equation (8) herein; -
FIG. 2 is a plot showing an example of a single exponential acquired at uniform time intervals in the time-domain; -
FIG. 3 is a plot showing a non-uniform sampling, according to some embodiments; -
FIG. 4 plots examples of non-uniform sampling in time domain that can directly provide moments of relaxation or diffusion distribution, according to some embodiments; -
FIG. 5 is a plot showing the ratio of the error-bar in the moments obtained from uniform sampling to non-uniform sampling; -
FIG. 6A shows an NMR tool being deployed in a wellbore, according to embodiments; -
FIG. 6B is a schematic of a processing system used to control, record and process data from an NMR tool; and -
FIG. 7 illustrates another example of a wellsite system in which an NMR tool can be employed, according to some embodiments. - The following description provides exemplary embodiments only, and is not intended to limit the scope, applicability, or configuration of the disclosure. Rather, the following description of the exemplary embodiments will provide those skilled in the art with an enabling description for implementing one or more exemplary embodiments. It being understood that various changes may be made in the function and arrangement of elements without departing from the spirit and scope of the invention as set forth in the appended claims.
- Specific details are given in the following description to provide a thorough understanding of the embodiments. However, it will be understood by one of ordinary skill in the art that the embodiments may be practiced without these specific details. For example, systems, processes, and other elements in the invention may be shown as components in block diagram form in order not to obscure the embodiments in unnecessary detail. In other instances, well-known processes, structures, and techniques may be shown without unnecessary detail in order to avoid obscuring the embodiments. Further, like reference numbers and designations in the various drawings indicated like elements.
- Also, it is noted that individual embodiments may be described as a process which is depicted as a flowchart, a flow diagram, a data flow diagram, a structure diagram, or a block diagram. Although a flowchart may describe the operations as a sequential process, many of the operations can be performed in parallel or concurrently. In addition, the order of the operations may be re-arranged. A process may be terminated when its operations are completed, but could have additional steps not discussed or included in a figure. Furthermore, not all operations in any particularly described process may occur in all embodiments. A process may correspond to a method, a function, a procedure, a subroutine, a subprogram, etc. When a process corresponds to a function, its termination corresponds to a return of the function to the calling function or the main function.
- Furthermore, embodiments of the invention may be implemented, at least in part, either manually or automatically. Manual or automatic implementations may be executed, or at least assisted, through the use of machines, hardware, software, firmware, middleware, microcode, hardware description languages, or any combination thereof. When implemented in software, firmware, middleware or microcode, the program code or code segments to perform the necessary tasks may be stored in a machine readable medium. A processor(s) may perform the necessary tasks.
- According to some embodiments, modified pulse sequences are provided that allow for the direct measurements of moments of relaxation-time or diffusion distributions. The data from these pulse sequences provides direct estimates of fluid properties such as average chain length and viscosity of a hydrocarbon. In comparison to the uniformly-spaced pulse sequence, these pulse sequences are faster and have a lower error bar in computing the fluid properties.
- Co-pending U.S. Patent Application No. 61/242,218, filed Sep. 14, 2009 (incorporated herein by reference) describes a method to analyze NMR data that obviates the estimation of the T1, T2, or D distributions. A new mathematical formulation is described where the moments of a random variable can be directly computed from the measurement. The method has two distinct advantages: the computation of moments is linear and does not require prior knowledge or information of the expected distribution. A moment of a random variable x with density function fx(x) is defined as,
-
- For NMR applications, the variable x corresponds to relaxation times T1, T2, or diffusion D. The moments of relaxation times or diffusion can be directly computed from the measured magnetization data,
-
- where Γ( ) represents the Gamma function and
-
- The contribution of variable ω is in two parts: a real number μ and an integer n where the mathematical operator tμ-1 operates on the n-th derivative of the data. Given a value of ω, the values of n and μ are unique and given as,
-
ω>0: n=0,μ=107 -
ω≦0: ω=μ−n, 0<μ≦1, n=[−ω]+1 (8) - where [ω] refers to the integral part of real number ω.
FIG. 1 is aplot 110 showing the values of n and μ for different values of ω as determined by above these equations. - According to some embodiments, the computation of moments of relaxation times T1, T2, or diffusion D, in conjunction with the scaling laws, can lead to greatly improved alternative pulse sequences to specifically target or estimate only certain moments of relaxation time/diffusion. In co-pending U.S. Patent Application Publication No. 2008-0314582 (incorporated herein by reference), “targeted” measurements are described as measurements that target a certain specified parameter at a particular location in the reservoir space. These targeted measurements are achieved by both hardware and software modifications.
- According to some embodiments, the concept of targeted measurements is applied to NMR data and alternate pulse sequences are described to directly estimate the moments of relaxation time or diffusion coefficients. According to some embodiments, alternate pulse sequences are designed by taking advantage of the specific exponential-kernel structure in eqn. (1) as well as the computation of moments using eqn. (7). The magnetization data from the new pulse sequences can be directly interpreted to obtain fluid properties.
- According to some embodiments, further detail will now be provided for modified pulse sequences for measuring NMR data. The non-uniformly spaced pulse sequences are targeted to specifically measure one of the moments of the relaxation time T1, T2, or D distribution. For simplicity, consider that the random variable in eqn. (1) is T1 and a specific ω-th moment of T1 is desired, where ω>0. This corresponds to the right-hand side of the ω axis in
FIG. 1 . In this case, from eqns. (7) and (8), -
- We introduce a new variable where
-
τ=tω (10) - Therefore, eqn. (9) can be re-written as
-
- where Mτ(τ) is the magnetization data in the τ domain. Eqn. (11) indicates that if the data are acquired such that they are uniformly spaced in the τ domain, then the sum of the measured data directly provides the ω-th moment. From eqn. (10), uniform spacing in the τ domain implies non-uniform spacing in the time (t) domain.
- We illustrate this with an example.
FIG. 2 is a plot showing an example of a single exponential acquired at uniform time intervals in the time-domain. In particular,curve 210 plots a single exponential decay of the form M(t)=e−t/T1 in the time-domain as shown in where the data are uniformly sampled in time. In this example, suppose that the 0.5-th moment of T1 is required.FIG. 3 is a plot showing a non-uniform sampling, according to some embodiments. In particular, using eqn. (10), the data shown bycurve 310 are sampled uniformly in the √{square root over (t)} domain.Curve 310 is an example of a single exponential in the time-domain into a Guassian in the τ domain with M(t)=e−r2 /T1 . Using eqn. (11), the sum of measured echoes of the Gaussian directly provides T1 0.5 . Thus, if the ω-th moment is desired, then the data need to be sampled non-uniformly in time according to eqn. (10).FIG. 4 plots examples of non-uniform sampling in time domain that can directly provide moments of relaxation or diffusion distribution, according to some embodiments. In particular, plots 410, 412, 414 and 416 represent sampling for ω values of 1, 0.75, 0.5, and 0.25 respectively. When 0<ω<1, this leads to finer sampling at the initial part of the data and coarser sampling at the tail end of the data. - The concept of non-uniform sampling is also valid when ω<0. Consider the case when −1<ω≦0. In this case, eqn. (7) reduces to
-
- Eqn. (12) can be re-written as
-
- Using integration by parts,
-
- When 0<μ<1, the first term in eqn. (15) is zero. The second term is integrable at t=0 since [M(t)−M(0)]˜t as t→0. In this case, we introduce a new variable τ where
-
τ=tμ=1 (16) - Here, τ→∞ when t→0, and τ→0 when t→∞. Therefore, eqn. (15) becomes
-
- Similar to eqn. (11), the area under the integrand [MT−M(0)] can be used to directly estimate the ω-th moment. As shown herein below, the concept of the pulse sequence for ω>0 and −1<ω≦0 is valid for other ranges of ω as well. For example, when −2<ω≦−1, pulse sequences that directly measure the ω-th moment are obtained by Taylor series expansion of the data with second-order terms in eqn. (13). Higher-order terms in the Taylor series lead to further negative values of ω. As shown herein below, when ω is not a negative integer, T1 ω can be obtained by sampling uniformly in the τ domain, where
-
τ=tω (18) - When ω is a negative integer, from eqns. (7) and (8), the ω-th moment is obtained by the |ω|-th derivative of the data,
-
- In addition to computing the moments directly, this modified pulse sequence also enables fluid properties such as the average chain length and viscosity to be directly estimated in a simple manner. Re-writing eqns. (3) and (4), the average chain length and viscosity can be computed from moments of T1 or T2 relaxation time or diffusion, as
-
- According to eqns. (20) and (21), the average chain length can be estimated from a modified pulse sequence to directly estimate the (1/κ)-th moment of relaxation time or (1/ν)-th moment of diffusion. Similarly, viscosity can be estimated from two modified pulse sequences to directly estimate the (1/κ)-th and (1/ν)-th moment of relaxation time. For diffusion data, the kernel in eqn. (1) is e−Dt. Using eqn. (10) with ω=1, the negative first moment of diffusion and viscosity are obtained.
- In addition to directly providing the parameters of interest, these pulse sequences are optimized to estimate any desired moments of relaxation time and/or diffusion. In comparison to pulse sequences where the data is uniformly acquired in time, these pulse sequences can be used to acquire the data faster. For example, given ω>0, the error-bar in the ω-th moment of T1 relaxation is,
-
- Here Ne denotes the total number of echoes obtained from uniform sampling between a given lower and upper limit of relaxation time and tE corresponds to the spacing between the echoes. In the case of non-uniform sampling, the error-bar in the ω-th moment is,
-
- Here τE corresponds to the spacing between echoes in the τ domain.
FIG. 5 is a plot showing the ratio of the error-bar in the moments obtained from uniform sampling to non-uniform sampling. There are two points to note from thiscurve 510. First, there is a significant advantage to non-uniform sampling as the error-bar in the computation of moments is lower (or at worst equal when ω=1). Second, from eqns. (23) and (24), the standard deviation of noise for uniform sampling has to be smaller than that for non-uniform sampling to obtain a desired error-bar in the moments. Since the magnitude of the noise standard deviation in the data is inversely proportional to the square-root of the acquisition time, to achieve the same uncertainty in the parameter the data for non-uniform sampling can be acquired faster than that for uniform sampling. - Thus, according to some embodiments, a modified pulse sequence is provided that can directly provide moments of relaxation-time or diffusion distributions. This pulse sequence can be adapted to the desired moment of relaxation-time or diffusion coefficient. The data from this pulse sequence provides direct estimates of fluid properties such as average chain length and viscosity of a hydrocarbon. In general, in comparison with data acquired uniformly in time, this pulse sequence is faster and has a lower error bar in computing the fluid properties.
- A general expression for the sampling rate of non-uniformly spaced pulse sequence will now be described in further detail, according to some embodiments. When ω is not a negative integer, <T1 ω> can be obtained by sampling uniformly in the τ domain, where
-
τ=tω (25) - Proof: For ω>0 and 1<ω<0, the proof has already been demonstrated in eqn. (11) and eqns. (12)-(17) respectively. When −2<ω≦1, from eqn. (8), n=2 and μ=ω+2. In this case, from eqn. (7),
-
- where δμ1=1 only when μ=1 and is zero otherwise and
-
- When μ=1 (corresponding to ω=1), the first term in eqn. (26) corresponds to the first derivative of the data at t=0 and is consistent with the definition of the negative first moment. In this case, the second term is zero. When −2<μ<1, the first term in eqn. (26) is zero. The second term is integrable at t=0. Using eqn. (25), the area under the integrand [M(t)−M(0)−tM(1)(0)] can be used to directly estimate the ω-th moment.
- In general, when n≧1, eqn. (7) can be written as
-
- Note that:
-
- As before, the first term in eqn. (27) is valid only for negative integer values of ω and corresponds to the |ω|-th derivative of the data. For negative values of ω where n<ω<(n-1), the first term in eqn. (27) is zero. Using eqn. (25), the area under the integrand
-
- in the τ domain directly provides the ω-th moment of relaxation time.
-
FIG. 6A shows an NMR tool being deployed in a wellbore, according to embodiments.Wireline truck 610 is deployingwireline cable 612 into well 630 viawell head 620.Wireline tool 640 is disposed on the end of thecable 612 in asubterranean formation 600.Tool 640 includes an NMR tool, such as Schlumberger's Combinable Magnetic Resonance (CMR) Tool, or Schlumberger's MR Scanner Tool. According to some embodiments the NMR tool is combined with other downhole tools such as Schlumberger's Modular Formation Dynamics Tester downhole sampling tool. According to some embodiments, the NMR tool uses a pulse sequence as described herein. Data from the NMR tool can be recorded and/or processed downhole withintool 640, and/or can be transmitted towireline truck 610 for recording and/or processing. According to embodiments, the NMR tool can be controlled, including the selection of pulse sequences, locally with processing systems withintool 640 and/or from the surface using processing systems for example inwireline truck 610. -
FIG. 6B is a schematic of a processing system used to control, record and process data from an NMR tool.Processing system 650 includes one or morecentral processing units 674,storage system 672, communications and input/output modules 670, a user display 676 and auser input system 678. Input/output modules 670 include modules to communicate with and control the NMR tool such that modified pulse sequences can be used as described herein. -
FIG. 7 illustrates another example of a wellsite system in which an NMR tool can be employed, according to some embodiments. The wellsite can be onshore or offshore. In this exemplary system, aborehole 711 is formed in subsurface formations by rotary drilling in a manner that is well known. Embodiments of the invention can also use directional drilling, as will be described hereinafter. - A
drill string 712 is suspended within theborehole 711 and has abottom hole assembly 700 which includes adrill bit 705 at its lower end. The surface system includes platform andderrick assembly 710 positioned over theborehole 711, theassembly 710 including a rotary table 716,kelly 717,hook 718 androtary swivel 719. Thedrill string 712 is rotated by the rotary table 716, energized by means not shown, which engages thekelly 717 at the upper end of the drill string. Thedrill string 712 is suspended from ahook 718, attached to a traveling block (also not shown), through thekelly 717 and arotary swivel 719 which permits rotation of the drill string relative to the hook. As is well known, a top drive system could alternatively be used. - In the example of this embodiment, the surface system further includes drilling fluid or
mud 726 stored in apit 727 formed at the well site. Apump 729 delivers thedrilling fluid 726 to the interior of thedrill string 712 via a port in theswivel 719, causing the drilling fluid to flow downwardly through thedrill string 712 as indicated by thedirectional arrow 708. The drilling fluid exits thedrill string 712 via ports in thedrill bit 705, and then circulates upwardly through the annulus region between the outside of the drill string and the wall of the borehole, as indicated by thedirectional arrows 709. In this well known manner, the drilling fluid lubricates thedrill bit 705 and carries formation cuttings up to the surface as it is returned to thepit 727 for recirculation. - The
bottom hole assembly 700 of the illustrated embodiment a logging-while-drilling (LWD)module 720, a measuring-while-drilling (MWD)module 730, a roto-steerable system and motor, anddrill bit 705. - The
LWD module 720 is housed in a special type of drill collar, as is known in the art, and can contain one or a plurality of known types of logging tools. It will also be understood that more than one LWD and/or MWD module can be employed, e.g. as represented at 720A. (References, throughout, to a module at the position of 720 can alternatively mean a module at the position of 720A as well.) The LWD module includes capabilities for measuring, processing, and storing information, as well as for communicating with the surface equipment. In the present embodiment, the LWD module includes an NMR tool such as Schlumberger's ProVISION NMR tool. - The
MWD module 730 is also housed in a special type of drill collar, as is known in the art, and can contain one or more devices for measuring characteristics of the drill string and drill bit. The MWD tool further includes an apparatus (not shown) for generating electrical power to the downhole system. This may typically include a mud turbine generator powered by the flow of the drilling fluid, it being understood that other power and/or battery systems may be employed. In the present embodiment, the MWD module includes one or more of the following types of measuring devices: a weight-on-bit measuring device, a torque measuring device, a vibration measuring device, a shock measuring device, a stick slip measuring device, a direction measuring device, and an inclination measuring device. - According to some embodiments, the techniques described herein are applied to estimate properties from other types of excitations than NMR. For example the techniques can be applied to other types of data where excitations induce an exponential decay. According to some embodiments, the techniques are applied to inducing a non-uniform sequence of pressure changes in a material. According to some other embodiments, the techniques are applied to induce a non-uniform sequence of temperature changes in a material.
- While many examples have been described herein with respect to designing interval spacings in the time domain, according to some embodiments, the techniques described herein are applied to interval spacings in other domains, such as designing interval spacings in the frequency domain.
- Whereas many alterations and modifications of the present disclosure will no doubt become apparent to a person of ordinary skill in the art after having read the foregoing description, it is to be understood that the particular embodiments shown and described by way of illustration are in no way intended to be considered limiting. Further, the disclosure has been described with reference to particular preferred embodiments, but variations within the spirit and scope of the disclosure will occur to those skilled in the art. It is noted that the foregoing examples have been provided merely for the purpose of explanation and are in no way to be construed as limiting of the present disclosure. While the present disclosure has been described with reference to exemplary embodiments, it is understood that the words, which have been used herein, are words of description and illustration, rather than words of limitation. Changes may be made, within the purview of the appended claims, as presently stated and as amended, without departing from the scope and spirit of the present disclosure in its aspects. Although the present disclosure has been described herein with reference to particular means, materials and embodiments, the present disclosure is not intended to be limited to the particulars disclosed herein; rather, the present disclosure extends to all functionally equivalent structures, methods and uses, such as are within the scope of the appended claims.
Claims (22)
1. A method of estimating a property of a material, comprising:
inducing a change in the material comprising a modified pulse sequence;
obtaining a measurement of a response to the inducing, wherein the response comprises exponential decay following the inducing;
taking weighted linear combinations of the measurements to determine moments of a variable from the measured exponentially decaying data; and
determining a property from the moments.
2. The method of claim 1 , wherein a pulse sequence is designed in time domain according to a formula determined by the property being estimated.
3. The method of claim 1 , wherein a pulse sequence in frequency domain is designed according to a formula determined by the property being estimated.
4. The method of claim 1 , wherein a pulse sequence is non-uniform in time in both linear and logarithmic scales.
5. The method of claim 1 , wherein a pulse sequence in time is predetermined using a power law.
6. The method of claim 1 , wherein a pulse sequence in frequency domain follows is predetermined using a power law.
7. The method of claim 1 , wherein a pulse sequence in time is selected for the property to be directly computed from the measurements.
8. The method of claim 1 , wherein a pulse sequence in a frequency domain is selected for the property to be directly computed from the measurements.
9. The method of claim 1 , wherein the weights in the weighted linear combination depend on the moments being estimated.
10. The method of claim 1 , wherein the properties are used to compute petro-physical parameters.
11. The method of claim 1 , wherein the material comprises a subterranean formation.
12. The method of claim 1 , wherein the material comprises water or hydrocarbon or both.
13. The method of claim 1 , wherein the properties comprise the mean chain length and viscosity.
14. The method of claim 1 , further comprising obtaining indications of property changes over temperature and pressure and
15. The method of claim 14 , further comprising determining phase transitions within the material from the changes in properties.
16. The method of claim 1 , further comprising determining an uncertainty in the properties from an uncertainty of the measurement.
17. The method of claim 1 , wherein the obtaining occurs in a subterranean formation traversed by a wellbore.
18. The method of claim 1 , wherein the material is a fluid.
19. The method of claim 1 , wherein the material is a rock solid and the properties are petrophysical properties of the rock solid.
20. The method of claim 1 , wherein the inducing comprises nuclear magnetic resonance
21. The method of claim 1 , wherein the property is a moment of relaxation-time or a diffusion distribution or a combination thereof.
22. A system for estimating a property of a material, comprising:
an inducer that is adapted to induce a measurable change in the material;
a sensor to measure a material response that exponentially decays following an induction;
a controller coupled to the inducer, adapted and programmed to repeat the induction at predetermined intervals; and
a processing system adapted to estimate the property based at least in part on the sensor measurements, wherein the intervals are predetermined comprising the material response.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/083,375 US20140074405A1 (en) | 2010-03-04 | 2013-11-18 | Modified pulse sequence to estimate properties |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/717,512 US8587302B2 (en) | 2010-03-04 | 2010-03-04 | Modified pulse sequence to estimate properties |
US14/083,375 US20140074405A1 (en) | 2010-03-04 | 2013-11-18 | Modified pulse sequence to estimate properties |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/717,512 Continuation US8587302B2 (en) | 2010-03-04 | 2010-03-04 | Modified pulse sequence to estimate properties |
Publications (1)
Publication Number | Publication Date |
---|---|
US20140074405A1 true US20140074405A1 (en) | 2014-03-13 |
Family
ID=44530789
Family Applications (4)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/717,512 Active 2031-10-22 US8587302B2 (en) | 2010-03-04 | 2010-03-04 | Modified pulse sequence to estimate properties |
US13/543,393 Active US9671516B2 (en) | 2010-03-04 | 2012-07-06 | Modified pulse sequence to estimate properties |
US14/083,375 Abandoned US20140074405A1 (en) | 2010-03-04 | 2013-11-18 | Modified pulse sequence to estimate properties |
US15/590,963 Active 2030-11-11 US10557809B2 (en) | 2010-03-04 | 2017-05-09 | Modified pulse sequence to estimate properties |
Family Applications Before (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/717,512 Active 2031-10-22 US8587302B2 (en) | 2010-03-04 | 2010-03-04 | Modified pulse sequence to estimate properties |
US13/543,393 Active US9671516B2 (en) | 2010-03-04 | 2012-07-06 | Modified pulse sequence to estimate properties |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/590,963 Active 2030-11-11 US10557809B2 (en) | 2010-03-04 | 2017-05-09 | Modified pulse sequence to estimate properties |
Country Status (5)
Country | Link |
---|---|
US (4) | US8587302B2 (en) |
CA (1) | CA2776572A1 (en) |
GB (1) | GB2486380B (en) |
NO (1) | NO20120402A1 (en) |
WO (1) | WO2011109207A2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170241922A1 (en) * | 2010-03-04 | 2017-08-24 | Schlumberger Technology Corporation | Modified pulse sequence to estimate properties |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
BR112014010273B1 (en) * | 2011-10-31 | 2022-06-21 | Schlumberger Holdings Limited | Method for petrophysically regularized time domain nuclear magnetic resonance (nmr) inversion |
US10031255B2 (en) * | 2014-03-24 | 2018-07-24 | Schlumberger Technology Corporation | Multi-dimensional nuclear magnetic resonance methods for characterizing fluids |
WO2016014377A2 (en) * | 2014-07-23 | 2016-01-28 | Schlumberger Canada Limited | Tar mat formation prediction in late-charge reservoirs |
US10392936B2 (en) | 2014-07-23 | 2019-08-27 | Schlumberger Technology Corporation | Tar mat formation prediction in late-charge reservoirs |
US10041893B2 (en) * | 2014-12-23 | 2018-08-07 | Schlumberger Technology Corporation | Methods and systems for identifying hydrocarbon fluid transition characteristics using nuclear magnetic resonance |
CA2963129C (en) * | 2016-04-12 | 2020-05-26 | Syncrude Canada Ltd. In Trust For The Owners Of The Syncrude Project As Such Owners Exist Now And In The Future | Low-field time-domain nmr measurement of oil sands process streams |
WO2018163188A1 (en) * | 2017-03-09 | 2018-09-13 | B.G. Negev Technologies And Applications Ltd., At Ben-Gurion University | Generation of nuclear magnetic resonance multidimensional t1(spin-matrix)-t2(spin-spin) energy relaxation maps and uses thereof |
WO2022031533A1 (en) * | 2020-08-05 | 2022-02-10 | Board Of Regents, The University Of Texas System | Systems and methods for automated, real-time analysis and optimization of formation-tester measurements |
Family Cites Families (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5291137A (en) * | 1992-11-02 | 1994-03-01 | Schlumberger Technology Corporation | Processing method and apparatus for processing spin echo in-phase and quadrature amplitudes from a pulsed nuclear magnetism tool and producing new output data to be recorded on an output record |
US5486762A (en) * | 1992-11-02 | 1996-01-23 | Schlumberger Technology Corp. | Apparatus including multi-wait time pulsed NMR logging method for determining accurate T2-distributions and accurate T1/T2 ratios and generating a more accurate output record using the updated T2-distributions and T1/T2 ratios |
US5343151A (en) * | 1993-03-11 | 1994-08-30 | Bruker Instruments, Inc. | Method for automatically shimming a high resolution NMR magnet |
US6005389A (en) * | 1996-03-15 | 1999-12-21 | Numar Corporation | Pulse sequences and interpretation techniques for NMR measurements |
GB9803487D0 (en) * | 1998-02-20 | 1998-04-15 | Inst Of Food Research | Nuclear magnetic resonance spectroscopy |
US6891369B2 (en) * | 1998-08-13 | 2005-05-10 | Schlumberger Technology Corporation | Nuclear magnetic resonance method and logging apparatus for fluid analysis |
US6600315B1 (en) * | 2000-03-03 | 2003-07-29 | Schlumberger Technology Corporation | Method for improving resolution of nuclear magnetic resonance measurements by combining low resolution high accuracy measurements with high resolution low accuracy measurements |
DK1301808T3 (en) | 2000-07-21 | 2010-04-06 | Schlumberger Technology Bv | Method and apparatus for analyzing nuclear magnetic resonance data |
US20020175682A1 (en) * | 2001-05-23 | 2002-11-28 | Songhua Chen | Rapid nmr multi-frequency t1 and t2 acquisition for earth formations evaluation with mwd or wireline tools |
US7049815B2 (en) | 2004-01-27 | 2006-05-23 | Baker Hughes Incorporated | Method and apparatus for multi-frequency NMR diffusion measurements in the presence of internal magnetic field gradients |
US7053611B2 (en) * | 2004-06-04 | 2006-05-30 | Schlumberger Technology Corporation | Method and apparatus for using pulsed field gradient NMR measurements to determine fluid properties in a fluid sampling well logging tool |
US7398159B2 (en) | 2005-01-11 | 2008-07-08 | Schlumberger Technology Corporation | System and methods of deriving differential fluid properties of downhole fluids |
US7305306B2 (en) | 2005-01-11 | 2007-12-04 | Schlumberger Technology Corporation | System and methods of deriving fluid properties of downhole fluids and uncertainty thereof |
MX2008015209A (en) * | 2006-06-06 | 2008-12-15 | Halliburton Energy Serv Inc | Fluid typing. |
US7622919B2 (en) * | 2006-07-31 | 2009-11-24 | Schlumberger Technology Corporation | Nuclear magnetic resonance measurement techniques in non-uniform fields |
US7486070B2 (en) * | 2006-12-18 | 2009-02-03 | Schlumberger Technology Corporation | Devices, systems and methods for assessing porous media properties |
US7301339B1 (en) * | 2006-12-26 | 2007-11-27 | Schlumberger Technology Corporation | Estimating the concentration of a substance in a sample using NMR |
US20080314582A1 (en) | 2007-06-21 | 2008-12-25 | Schlumberger Technology Corporation | Targeted measurements for formation evaluation and reservoir characterization |
US7804296B2 (en) * | 2007-10-05 | 2010-09-28 | Schlumberger Technology Corporation | Methods and apparatus for monitoring a property of a formation fluid |
CA2638949C (en) * | 2008-08-20 | 2011-11-15 | Schlumberger Canada Limited | Methods of and apparatus for determining the viscosity of heavy oil |
WO2011059543A1 (en) * | 2009-11-10 | 2011-05-19 | The Johns Hopkins University | Magnetic resonance imaging and spectroscopy of low concentration solutes with exchangeable protons using label transfer modules: frequency transfer, inversion transfer, and dephasing transfer |
US8587302B2 (en) * | 2010-03-04 | 2013-11-19 | Schlumberger Technology Corporation | Modified pulse sequence to estimate properties |
US9222902B2 (en) * | 2012-01-09 | 2015-12-29 | Schlumberger Technology Corporation | Estimations of nuclear magnetic resonance measurement distributions |
US8965094B2 (en) * | 2012-04-14 | 2015-02-24 | Nocimed, Llc | Magnetic resonance spectroscopy pulse sequence, acquisition, and processing system and method |
WO2014022611A1 (en) * | 2012-08-01 | 2014-02-06 | Schlumberger Canada Limited | Single well inject-produce pilot for eor |
GB2532380B (en) * | 2013-12-31 | 2020-06-24 | Halliburton Energy Services Inc | Encoded driving pulses for a range finder |
-
2010
- 2010-03-04 US US12/717,512 patent/US8587302B2/en active Active
-
2011
- 2011-02-24 CA CA2776572A patent/CA2776572A1/en not_active Abandoned
- 2011-02-24 WO PCT/US2011/026008 patent/WO2011109207A2/en active Application Filing
- 2011-02-24 GB GB1205890.5A patent/GB2486380B/en not_active Expired - Fee Related
-
2012
- 2012-04-03 NO NO20120402A patent/NO20120402A1/en not_active Application Discontinuation
- 2012-07-06 US US13/543,393 patent/US9671516B2/en active Active
-
2013
- 2013-11-18 US US14/083,375 patent/US20140074405A1/en not_active Abandoned
-
2017
- 2017-05-09 US US15/590,963 patent/US10557809B2/en active Active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170241922A1 (en) * | 2010-03-04 | 2017-08-24 | Schlumberger Technology Corporation | Modified pulse sequence to estimate properties |
US10557809B2 (en) * | 2010-03-04 | 2020-02-11 | Schlumberger Technology Corporation | Modified pulse sequence to estimate properties |
Also Published As
Publication number | Publication date |
---|---|
US8587302B2 (en) | 2013-11-19 |
US10557809B2 (en) | 2020-02-11 |
US20110215802A1 (en) | 2011-09-08 |
US9671516B2 (en) | 2017-06-06 |
WO2011109207A2 (en) | 2011-09-09 |
NO20120402A1 (en) | 2012-05-04 |
WO2011109207A3 (en) | 2012-01-12 |
GB2486380B (en) | 2015-12-30 |
GB201205890D0 (en) | 2012-05-16 |
CA2776572A1 (en) | 2011-09-09 |
GB2486380A (en) | 2012-06-13 |
US20130002246A1 (en) | 2013-01-03 |
US20170241922A1 (en) | 2017-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8587302B2 (en) | Modified pulse sequence to estimate properties | |
US5696448A (en) | NMR system and method for formation evaluation using diffusion and relaxation log measurements | |
Hürlimann et al. | Quantitative measurement of two-dimensional distribution functions of diffusion and relaxation in grossly inhomogeneous fields | |
US7495436B2 (en) | Rock properties prediction, categorization, and recognition from NMR echo-trains using linear and nonlinear regression | |
US6859032B2 (en) | Method for determining molecular properties of hydrocarbon mixtures from NMR data | |
US20190271224A1 (en) | Correction of motion effect in nuclear magnetic resonance (nmr) logging | |
US6859034B2 (en) | Time-domain data integration of multiple gradient, multiple TE echo trains | |
US7221158B1 (en) | Permeability determinations from nuclear magnetic resonance measurements | |
US7821260B2 (en) | NMR echo train compression using only NMR signal matrix multiplication to provide a lower transmission bit parametric representation from which estimate values of earth formation properties are obtained | |
US9753176B2 (en) | Estimating adsorbed gas volume from NMR and dielectric logs | |
US10539705B2 (en) | Formation evaluation utilizing dual wait time nuclear magnetic resonance | |
US20130057277A1 (en) | Method for Determining Rock Formation Fluid Interaction Properties Using Nuclear Magnetic Resonance Well Logging Measurements | |
US7224162B2 (en) | System and methods for upscaling petrophysical data | |
EP2062073B1 (en) | Nmr echo train compression | |
US9915750B2 (en) | Methods and apparatuses to remove a net detected residual magnetization in a nuclear magnetic resonance (NMR) operation | |
US20090066327A1 (en) | Method and computer program product for estimating true intrinsic relaxation time and internal gradient from multigradient nmr logging | |
WO2015085292A1 (en) | Downhole monitoring of fluids using nuclear magnetic resonance | |
US20160124109A1 (en) | Semi-Analytic Inversion Method For Nuclear Magnetic Resonance (NMR) Signal Processing | |
US20090195246A1 (en) | Method for estimating insitu fluid velocity from nmr measurements |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |