The present invention concerns logging applications and more specifically relates to a method for data acquisition and processing, which uses time-domain averaging to increase the signal to noise ratio (SNR) or to sharpen the spatial resolution of the obtained logging data. In an exemplary application, the method is used in nuclear magnetic resonance (NMR) logging to increase the SNR of single event echo-train for high-speed NMR logging and/or to sharpen the spatial resolution of the received data.
In oil and gas exploration it is desirable to understand the structure and properties of the geological formation surrounding a borehole, in order to determine if the formation contains hydrocarbon resources (oil and/or gas), to estimate the amount and producibility of hydrocarbon contained in the formation, and to evaluate the best options for completing the well in production. A significant aid in this evaluation is the use of wireline logging and/or logging-while-drilling (LWD) measurements of the formation surrounding the borehole (referred to collectively as “logs” or “log measurements”). Typically, one or more logging tools are lowered into the borehole and the tool readings or measurement logs are recorded as the tools traverse the borehole. These measurement logs are used to infer the desired formation properties.
In recent years nuclear magnetic resonance (NMR) logging has become very important for purposes of formation evaluation and is one of the preferred methods for determining formation parameters. Improvements in the NMR logging tools, as well as advances in data analysis and interpretation allow log analysts to generate detailed reservoir description reports, including clay-bound and capillary-bound related porosity, estimates of the amounts of bound and free fluids, fluid types (i.e., oil, gas and water), permeability and other properties of interest.
NMR tools used in practical applications include, for example, the centralized MRIL® tool made by NUMAR Corporation, a Halliburton company, and the sidewall CMR tool made by Schlumberger. The MRIL® tool is described, for example, in U.S. Pat. No. 4,710,713 to Taicher et al. and in various other publications including: “Spin Echo Magnetic Resonance Logging: Porosity and Free Fluid Index Determination,” by Miller, Paltiel, Gillen, Granot and Bouton, SPE 20561, 65th Annual Technical Conference of the SPE, New Orleans, La., Sep. 23–26, 1990; “Improved Log Quality With a Dual-Frequency Pulsed NMR Tool,” by Chandler, Drack, Miller and Prammer, SPE 28365, 69th Annual Technical Conference of the SPE, New Orleans, La., Sep. 25–28, 1994. Certain details of the structure and the use of the MRIL® tool, as well as the interpretation of various measurement parameters are also discussed in U.S. Pat. Nos. 4,717,876; 4,717,877; 4,717,878; 5,212,447; 5,280,243; 5,309,098; 5,412,320; 5,517,115, 5,557,200; 5,696,448; 5,936,405; 6,005,389; 6,023,164; 6,051,973; 6,107,796 and 6,111,408. The structure and operation of the Schlumberger CMR tool is described, for example, in U.S. Pat. Nos. 4,939,648; 5,055,787 and 5,055,788 and further in “Novel NMR Apparatus for Investigating an External Sample,” by Kleinberg, Sezginer and Griffin, J. Magn. Reson. 97, 466–485, 1992; and “An Improved NMR Tool Design for Faster Logging,” D. McKeon et al., SPWLA 40th Annual Logging Symposium, May–June 1999. The content of the above patents is hereby expressly incorporated by reference for all purposes, and all non-patent references are incorporated by reference for background.
NMR tools of the type discussed above generally measure the time for hydrogen nuclei present in the earth formation to realign their spin axes, and consequently their bulk magnetization, either with an externally applied magnetic field, or perpendicularly to the magnetic field, after momentary reorientation due to the application of specific radio frequency (RF) pulses. The externally applied magnetic field is typically provided by a magnet disposed in the tool. The spin axes of the hydrogen nuclei in the earth formation are, in the aggregate, caused to be aligned with the magnetic field induced in the earth formation by the magnet. The NMR tool includes an antenna positioned near the magnet and shaped so that a pulse of radio frequency (RF) power conducted through the antenna induces a magnetic field in the earth formation orthogonal to the field induced by the magnet. The RF pulse has a duration predetermined so that the spin axes of the hydrogen nuclei generally align themselves perpendicular both to the orthogonal magnetic field induced by the RF pulse and to the externally applied magnetic field. After the pulse ends, the nuclear magnetic moment of the hydrogen nuclei gradually relax, i.e., return to their alignment with the externally applied magnetic field; at the same time an antenna, which is typically the same as the one used by the initial pulse, is electrically connected to a receiver, which detects and measures voltages induced in the antenna by precessional rotation of the spin axes of the hydrogen nuclei.
An actual NMR measurement involves a plurality of pulses grouped into pulse sequences, most frequently of the type known in the art as Carr-Purcell-Meiboom-Gill (CMPG) pulsed spin echo sequences. As known in the art, each CPMG sequence consists of a 90-degree (i.e., π/2) pulse followed by a large number of 180-degree (i.e., π) pulses. The 90-degree pulse rotates the proton spins into the transverse plane and the 180-degree pulses generate a sequence of spin echoes by refocusing the transverse magnetization after each spin echo. A wait time interval (TW) of approximately 0.5–10 sec is used first to allow for polarization of the formation by the tool's static magnetic field. Then, a CPMG pulse-echo train is executed, consisting of an excitation pulse and an alternating sequence of refocusing pulses. Following each pair of excitation pulse and a refocusing pulse, an acquisition window is applied next. Complex data points from such a pair of echo trains are co-added on an echo-by-echo basis to remove certain artifacts and to enhance the NMR signal, as known in the art. More pairs may be added to enhance the signal-to-noise ratio. The echo train, consisting of a superposition of exponentially decaying signals is then submitted to a processor, which calculates the underlying decay modes of the received NMR echo signal.
One of the key information sought from the NMR logging data is the T2 decay time distribution, which is related to the pore size distribution of the rock formation. From the T2 decay time distribution one can estimate, for example, the capillary bound water volume. The T2 distribution is obtained by inversion of the echo train data, using well-known methods as described, for example, in U.S. Pat. No. 5,517,115 to the assignee of the present application.
It is also well-known in the art that inversion of NMR echo trains is greatly affected by noise.
In actual logging, the tool is moving at a certain speed, and echo train from each event is not collected at the same logging location or depth. When the logging speed is slow, the tool moves only a short distance. Hence one may assume that the earth formation is of nearly constant property over the short distance, and stacking of multiple events from different depths may be acceptable. However, when the logging speed is relatively fast, the tool moves over longer distances and the formation properties may change significantly over the distance. As a result, stacking of echo train data from multiple events, as generally used in the prior art, reduces apparent spatial resolution of measuremenets beyond the intrinsic resolution of the tool.
In particular, in high-signal formations (e.g., high porosity, oil-or water-filled rocks), where the logging speeds can be comparatively fast, the effective vertical resolution of the NMR log becomes dominated by the tool movement during a single PAP. For example, with a recovery time of 10 seconds between echo-trains in a PAP, with PAPs acquired at all possible frequencies, the elapsed time between the first echo-train in the first-frequency PAP and the second echo-train for the last-frequency PAP, is close to 20 seconds. At a logging speed of 900 ft/hr (15 ft/min), the MRIL tool will move approximately 5 feet during this measurement: when combined with the inherent vertical resolution of the antenna (which is approximately 2 feet), the effective vertical resolution becomes roughly 7 feet.
Enhancing the resolution of the logs is a significant problem, because subsurface formations are generally heterogeneous, so that porosity, saturation and lithology vary with position. A common example of heterogeneity is the presence in the formation of geological layers, or beds. Because logging tools have a nonzero volume of investigation, more than one layer may lie within the volume of investigation of a tool. In such cases, the petrophysical evaluation of one layer may be distorted by the presence of another layer falling within the larger volume of investigation of the tool. The above phenomenon leads to a specific problem in the analysis of subsurface formations that include one or more underground layers, especially when the layers are thin compared with the vertical resolution of the measuring tool. Such layers have become subject to significant commercial interest because of their production potential. Any knowledge about the composition and properties of such layered formations that helps better estimate their production potential has thus become increasingly valuable. It is clear therefore that any mechanism that for a given SNR supported by the formation can increase the vertical resolution of the tool without decreasing the logging speed is highly desirable.
In a preferred embodiment the present invention uses time-domain average of single event data to increase the SNR, but not stacking of data from multiple events, and hence, achieves higher spatial resolution. In accordance with the present invention, the SNR can be increased also by application of both the time-domain-average and stacking techniques. In such case, high spatial resolution can be achieved using a lesser number of stacking events. Various combinations of techniques known in the art with the teachings of the present invention are possible in this regard and will be appreciated by a person of skill in the art.
In particular, according to this invention a method for measuring an indication of attributes of materials containing a fluid state is disclosed, the method comprising the steps of: providing a time-domain signal indicative of attributes of said materials in a single measurement; constructing a time-domain averaged data train from said signal, the averaging being performed over one or more time intervals Δi; and computing an indication of attributes of said materials from the time-domain averaged data train. In different embodiments, the time intervals Δi can be selected constant or different from each other. In a specific preferred embodiment of the method the time-domain signal is an NMR echo train.
In accordance with another aspect of the invention, a method for increasing the spatial resolution of NMR logging measurements is disclosed, comprising the steps of: providing an NMR echo-train indicative of attributes of materials of interest; and constructing a time-domain averaged data train from said NMR echo train, the averaging being performed over one or more time intervals Δi. It should be understood that the method of the present invention is applicable to processing of data from virtually any type of logging device and thus is not limited to NMR logging applications.
The method in accordance with the present invention will be described next in the context of NMR logging and will use terminology generally accepted in this field. It should be understood, however, that the method is equally applicable to processing of data from any number of other logging applications, and more generally to processing of data that can be modeled as an exponential function of time.
With reference to the NMR logging applications discussed above, in accordance with the present invention, the T2 decay time distribution can be obtained by inverting the time-domain echo train data, assuming the following physical model:
where φ (T2) is the porosity (or population) of the pores corresponding to the exponential decay time T2. One such example of echo-train is plotted in
To reduce the effect of random noise, it is commonly practiced to stack echo-trains from a multiple number of events.
As indicated in
In accordance with the present invention, the signal-to-nose ratio (SNR) can be increased without stacking echo train data from multiple events (or by using reduced stacking) and thus without worsening the spatial resolution at high logging speeds, but by taking time-domain-average(s) of the input signal. The principal advantage of this approach is therefore that the spatial resolution of the tool over the coarse of a logging experiment is not worsened at the price of the resolution in the time domain.
With reference to the notations in
It is important to note that the noise in the averaged echo is still random. To ensure that, the averaged echo separated by time interval of Δ is considered. Namely, the averaged echo train at t=t0, t0+Δ, t0+2Δ, . . . , t0+NΔ is used to obtain T2 distribution.
It should be apparent that in different applications the value of the parameter Δ can be changed to conform with the underlying physics of the signal. Furthermore, in alternative embodiments of the invention, instead of a single time period Δ, one can use two or more periods Δi the values of which can generally be different. As in the preceding case, the values for Δi can be selected dependent on the practical application and in particular on the need to increase the SNR of the received signal.
In accordance with another aspect of the present invention, stacking of multiple-event echo-train data and time-domain averaging can be used in combination. For instance, instead of taking 10-event stacking or time-domain averaging over 10 sampling points, one can use 5-event stacking and then apply time-domain averaging to increase the SNR of the signal, without worsening the vertical resolution of the log as much as prior art methods would, which only use stacking.
It should be noted that while the method of the present invention has been described in the context of NMR logging, it is applicable to a multitude of practical applications which need not be limited to NMR logging. Virtually all types of logging measurements with relatively low SNR signals can be enhanced using the method outlined above, without unduly compromising the vertical resolution of the moving tool. It should be apparent that various combinations of prior art measurement methods and techniques can be used along with the time-domain averaging approach of the present invention. It should also be clear to those skilled in the art that the approach used in accordance with the present invention can also be used in measuring properties of samples of porous materials in a laboratory setting, as well as in situ logging-type including logging/measuring while drilling (LWD/MWD) measurements. The use of the method alone or in combination with other techniques, as discussed above is believed to be a significant contribution to the art of logging with wide ranging applications involving virtually all tools, and a broad range of practical applications, including logging while drilling.
While the invention has been described with reference to the preferred embodiments, it will be appreciated by those of ordinary skill in the art that modifications can be made to the structure and form of the invention without departing from its spirit and scope which is defined in the following claims.
Number | Name | Date | Kind |
---|---|---|---|
2702367 | Ergen | Feb 1955 | A |
2912641 | Ruble | Nov 1959 | A |
3213357 | Brown et al. | Oct 1965 | A |
3231356 | Brown et al. | Jan 1966 | A |
3402344 | Brown et al. | Sep 1968 | A |
3617867 | Herzog | Nov 1971 | A |
3638484 | Tixier | Feb 1972 | A |
3667035 | Slichter | May 1972 | A |
3881651 | Wilhelm, Jr. | May 1975 | A |
4389613 | Brown | Jun 1983 | A |
4412178 | Brown | Oct 1983 | A |
4412179 | Brown | Oct 1983 | A |
4536714 | Clark | Aug 1985 | A |
4710713 | Strikman | Dec 1987 | A |
4717876 | Masi et al. | Jan 1988 | A |
4717877 | Taicher et al. | Jan 1988 | A |
4717878 | Taicher et al. | Jan 1988 | A |
4728892 | Vinegar et al. | Mar 1988 | A |
4933638 | Kenyon et al. | Jun 1990 | A |
4939648 | O'Neill et al. | Jul 1990 | A |
5023551 | Kleinberg et al. | Jun 1991 | A |
5055787 | Kleinberg et al. | Oct 1991 | A |
5055788 | Kleinberg et al. | Oct 1991 | A |
5122746 | King et al. | Jun 1992 | A |
5212447 | Paltiel | May 1993 | A |
5280243 | Miller | Jan 1994 | A |
5291137 | Freedman | Mar 1994 | A |
5309098 | Coates et al. | May 1994 | A |
5350925 | Watson | Sep 1994 | A |
5363041 | Sezginer | Nov 1994 | A |
5381092 | Freedman | Jan 1995 | A |
5412320 | Coates | May 1995 | A |
5486762 | Freedman et al. | Jan 1996 | A |
5497087 | Vinegar et al. | Mar 1996 | A |
5498960 | Vinegar et al. | Mar 1996 | A |
5517115 | Prammer | May 1996 | A |
5557200 | Coates | Sep 1996 | A |
5557201 | Kleinberg et al. | Sep 1996 | A |
5680043 | Hurlimann et al. | Oct 1997 | A |
5696448 | Coates et al. | Dec 1997 | A |
5698979 | Taicher et al. | Dec 1997 | A |
5705927 | Sezginer et al. | Jan 1998 | A |
5936405 | Prammer et al. | Aug 1999 | A |
6005389 | Prammer | Dec 1999 | A |
6023164 | Prammer | Feb 2000 | A |
6051973 | Prammer | Apr 2000 | A |
6107796 | Prammer | Aug 2000 | A |
6111408 | Blades et al. | Aug 2000 | A |
6115671 | Fordham et al. | Sep 2000 | A |
6121774 | Sun et al. | Sep 2000 | A |
6133734 | McKeon | Oct 2000 | A |
6140817 | Flaum et al. | Oct 2000 | A |
6163153 | Reiderman et al. | Dec 2000 | A |
6366088 | Hagiwara | Apr 2002 | B1 |
6385548 | Ananthaiyer et al. | May 2002 | B1 |
6452389 | Edwards | Sep 2002 | B1 |
6459263 | Hawkes et al. | Oct 2002 | B1 |
6541969 | Sigal et al. | Apr 2003 | B1 |
6600315 | Heaton et al. | Jul 2003 | B1 |
20010033163 | Sigal et al. | Oct 2001 | A1 |
Number | Date | Country |
---|---|---|
0 649 035 | Apr 1995 | EP |
2 056 082 | Jul 1980 | GB |
WO 9825164 | Jun 1998 | WO |
Number | Date | Country | |
---|---|---|---|
20020163334 A1 | Nov 2002 | US |