US20070244646A1 - Method of formation characterication using tool responses - Google Patents
Method of formation characterication using tool responses Download PDFInfo
- Publication number
- US20070244646A1 US20070244646A1 US11/279,409 US27940906A US2007244646A1 US 20070244646 A1 US20070244646 A1 US 20070244646A1 US 27940906 A US27940906 A US 27940906A US 2007244646 A1 US2007244646 A1 US 2007244646A1
- Authority
- US
- United States
- Prior art keywords
- formation
- tool responses
- measurement data
- error distribution
- measurement error
- 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
- 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/26—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
- G01V3/28—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device using induction coils
Definitions
- the present invention generally relates to the logging of earth formations, and more particularly to methods and apparatus for formation characterization using logging tool responses corresponding to logging measurements.
- Wireline tools have been employed to obtain formation measurement.
- induction-logging instruments are susceptible to a variety of environmental effects due to considerable range of investigation. Therefore, to make useful logs, the effects of volume above and below the layer of the interest must be carefully removed.
- the Array Induction Imager (AIT) or AIT-family tools are popular in moderate conductive to oil-based-mud logging environment.
- the AIT tool consists of eight three-coil arrays, six of which are operated simultaneously at two frequencies (AIT-B) or one frequency (AIT-H/M family).
- the AIT tool can produce wellsite resistivity logs having high resolution and high rejection of borehole rugosity effects. For example, these tools can produce logs that are corrected for the most common environmental effects, e.g. borehole effect, shoulder effect, skin effect, et al.
- the post-processing typically is based on Gaussian error assumption and least-squares inversion, thus the post-processed result is Gaussian mean value.
- the Gaussian error assumption is not always correct, thus the traditional post-processed result might not be ideal either.
- the Gaussian error assumption is violated.
- AIT 10, 20, 30, 60 and 90-in curves can behave strangely, either formation inhomogeneity or vertical bedding being overlooked, or the separation of the five curves can not be explained from logging environment information. In such situation, the traditional AIT well-site post-processing result cannot provide accurate indication of formation characterization.
- One aspect of the invention is directed to a method of deriving formation characterization.
- the method comprises receiving measurement data of a subsurface area; defining a measurement error distribution; creating multiple parameterized formation models of the subsurface area based on the measurement error distribution and the measurement data; generating multiple tool responses corresponding to the multiple parameterized formation models; considering the generated tool responses corresponding to the multiple parameterized formation modes; and deriving the formation characterization based on the result of considering generated tool responses.
- Embodiments of the method may further include performing borehole correction of the measurement data on the multiple tool responses and performing vertical processing on the borehole corrected measurement data.
- the method may further comprise allowing a user to choose the measurement error distribution.
- the method may further comprise creating multiple parameterized formation models by using Monte Carlo sampling.
- the apparatus comprises means for receiving measurement data of a subsurface area; means for defining a measurement error distribution; means for creating multiple parameterized formation models of the subsurface area based on the measurement error distribution and the measurement data; means for generating multiple tool responses corresponding to the multiple parameterized formation models; means for considering the generated tool responses corresponding to the multiple parameterized formation models; and means for deriving the formation characterization based on the result of considering generated tool responses.
- FIG. 1 is a flowchart showing steps associated with the present method, apparatus, and article of manufacture
- FIG. 2 is a schematic illustration of computer hardware associated with the apparatus and article of manufacture
- FIGS. 3-4 are a related group of diagrams used to describe the experimental result based on the inventive method.
- FIGS. 5-7 are another related group of diagrams used to describe the experimental result based on the inventive method.
- FIG. 1 shows several steps associated with the present method, apparatus and article of manufacture and provides a general overview of the invention.
- the Receive Measurement Data Step 10 the Measurement Data of the subsurface area is received.
- the received Measurement Data is raw data 12 from arrays of Array Induction Imager (AIT) family tools.
- the AIT-family tools typically include AIT, PlatformExpress AIT (AIT-H), SlimAccess AIT (SAIT), Xtreme AIT (XAIT), and SlimXtreme AIT (QAIT), etc.
- the received Measurement Data is raw data 12 from arrays of High Resolution Laterolog Array (HRLA) tool.
- the received Measurement Data can be any measurement data (not shown in the drawing) of the subsurface area from a resistivity tool or other tools.
- the measurement error distribution is defined based on Gaussian error assumption.
- the Gaussian error assumption assumes that the error distribution is characterized by a mean value (most possible value in large amount of repeated loggings) and standard deviation ⁇ (the possibility of the measured value drops to 63.2%).
- the measurement error distribution is defined based on homogeneous error assumption. The homogeneous error assumption assumes that the possibilities of measured values are the same for certainty resistivity range.
- the user 16 can choose any existing mathematical functions to define the measurement error distribution, or even create his/own functions to define the measurement error distribution.
- the user 16 can apply different measurement error distribution definitions and may reach different results of the formation characterization after the following steps, and the user can choose one of the results that appears the most accurate one by comparing the different results of the formation characterization among each other.
- the parameterized formation models include a plurality of formation and wellbore related parameters, e.g. formation resistivity, invasion diameter, borehole size, mud resistivity, tool stand-off, Young's modulus, Poisson coefficient, shear modulus etc.
- a large amount of parameterized formation models are created using Monte Carlo sampling, i.e. after one parameterized formation model is created, the parameterized formation model is accepted with the chance defined by the possibility density function.
- the amount of the parameterized formation models can be 1,000 or even more depending on the user.
- multiple tool responses are generated corresponding to the multiple parameterized formation models created in the Create Formation Models Step 18 .
- multiple resistivity tool responses are generated by simulating the tools response from a forward modeling engine.
- the AIT-family log processing includes borehole correction as stated in the Perform Borehole Correction Step 22 .
- the raw data is corrected from all eight arrays for the borehole effect.
- borehole corrections for the AIT tools are based on inversion through an iterative forward model to find the borehole parameters that best reproduce the logs from the four shortest arrays—the 6-, 9-, 12-, and 15-in arrays.
- the borehole forward model is based on a solution to Maxwell's equations in a cylindrical borehole with resistivity Rm surrounded by a homogeneous formation of resistivity Rf.
- the AIT-family log processing also includes vertical processing as stated in the Vertical Processing Step 24 .
- the borehole corrected data are going through a non-linear procedure to compensate for skin-effect and shoulder effect.
- the generated tool response behavior corresponding to the multiple parameterized formation models is examined.
- the generated tool responses are averaged.
- the percentage of certain characters of the generated tool response behavior is examined. The user can define the threshold of percentage, e.g. 60% or 80%, for concluding the existence of certain character.
- the formation characterization is derived based on the examination results of the generated tool response behavior.
- the user 16 can define the percentage threshold for certain formation character. For example, the user may define 80% of the results showing curve separation as indication of a formation bed existence. From the statistics of the Monte Carlo simulated and processed results, if the curve separation is equal or more than 80%, the user can conclude that there exists a formation bed.
- the inventive method allows a user/engineer to define the data measurement error distribution.
- AIT response distribution behavior under a large number of formation model and borehole situations can be examined.
- the inventive method does not reply on the Gaussian error assumption so that when this assumption is violated, the inventive method can give a more reasonable result than traditional well-site post-processing method.
- the user can choose different measurement error distribution definitions and may reach different results of the formation characterization.
- the user can choose one of the results that appears the most accurate one by comparing the different results of the formation characterization among each other.
- the user can choose one of the results that appears the most accurate one by comparing the different results of the formation characterization with the results obtained by other means or tools.
- the inventive method can be helpful in formation uncertainty study and answering questions like: whether the inversion is sensitive enough for that parameter; whether there are other possibilities of formation models; whether more information is needed in drawing conclusions, etc.
- the inventive method can help us to benefit more from traditional AIT measurement interpretation.
- FIG. 2 schematically illustrates computer hardware that may be used to implement the inventive method.
- Computer 30 has a media reading device, such as a CD-ROM Reader 32 , a floppy disk device, or a ZIP drive.
- the media reading device may also be capable of recording the output of the program the computer 30 is running.
- a user of the computer 30 may enter commands using a user input device, such as a keyboard 34 or a mouse, may view output of the program code on a visual display device, such as monitor 36 , and may make hardcopies of output using an output device, such as printer 38 .
- computer 30 (and its associated peripheral devices) is an apparatus for deriving formation characterization in accordance with the present invention.
- Computer media such as a CD-ROM 40 , a floppy disk, or a ZIP disk, may have computer readable program code that allows the computer 30 to derive formation characterization in accordance with the inventive method.
- FIG. 3 shows the comparison results between well-site AIT post-processing result based on Gaussian error assumption and based on the inventive method. Specifically, formation anisotropy and vertical thin bed detection of a borehole and formation in La Copital field, Texas, US are examined. On the left, well-site AIT post-processing result, based on traditional Gaussian error assumption, indicates five curves of depth of investigation 10, 20, 30, 60 and 90-in overlap each other (left), showing no invasion, no anisotropy in the depth range 8403-8428 ft. However, the five curves processed with the inventive method show obvious separation and vertical bedding in this depth range (middle).
- GR nuclear tool that measures formation natural gamma ray radiation
- HRLA five curves overlap each other, showing a water zone with more or less constant resistivity Rt of 3 ohm ⁇ m.
- AIT was logged one day after HRLA logging, but the traditional AIT logging shows strange curve separation, which makes it difficult to estimate formation water resistivity, as shown in FIG. 6 .
- AIT shows obvious resistive invasion, where the deep curve (AT90) indicates a more conductive and variable formation water with resistivity ranging from 0.5 ohm ⁇ m to 2 ohm ⁇ m. Since the mud type is very conductive, by the time AIT was logged, conductive invasion can be strong and un-favorable to AIT measurement. If HRLA provides a better estimation in formation water, AIT deep resistivity curve measurement will give an estimation of formation water resistivity only 33% of what has been estimated from HRLA measurements.
- the discrepancy can be greatly eliminated.
- the uncertainty method result agrees quite well with HRLA measurement.
- the AIT deep curve indicates a water zone with resistivity around 3 ohm ⁇ m.
- the Monte Carlo large sample results show that quite possibly there was something wrong of the AIT hardware during logging.
- HRLA alternative measurement
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Electromagnetism (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
A method of deriving formation characterization comprises receiving measurement data of a substrate area; defining a measurement error distribution; creating multiple parameterized formation models of the subsurface area based on the measurement error distribution and the measurement data; generating multiple tool responses corresponding to the multiple parameterized formation models; considering the generated tool responses corresponding to the multiple parameterized formation models; and deriving the formation characterization based on the result of considering generated tool responses.
Description
- The present invention generally relates to the logging of earth formations, and more particularly to methods and apparatus for formation characterization using logging tool responses corresponding to logging measurements.
- Wireline tools have been employed to obtain formation measurement. In certain prior art apparatus, induction-logging instruments are susceptible to a variety of environmental effects due to considerable range of investigation. Therefore, to make useful logs, the effects of volume above and below the layer of the interest must be carefully removed. In practice, the Array Induction Imager (AIT) or AIT-family tools are popular in moderate conductive to oil-based-mud logging environment. The AIT tool consists of eight three-coil arrays, six of which are operated simultaneously at two frequencies (AIT-B) or one frequency (AIT-H/M family). The AIT tool can produce wellsite resistivity logs having high resolution and high rejection of borehole rugosity effects. For example, these tools can produce logs that are corrected for the most common environmental effects, e.g. borehole effect, shoulder effect, skin effect, et al.
- However, such tools typically need a wellsite post-processing to remove the environmental effects. Due to the high requirement of processing speed at well-sites, as well as the difficulty in quantifying the logging measurement error, the post-processing typically is based on Gaussian error assumption and least-squares inversion, thus the post-processed result is Gaussian mean value. In reality, however, the Gaussian error assumption is not always correct, thus the traditional post-processed result might not be ideal either. For example, if the Gaussian error assumption is violated. AIT 10, 20, 30, 60 and 90-in curves can behave strangely, either formation inhomogeneity or vertical bedding being overlooked, or the separation of the five curves can not be explained from logging environment information. In such situation, the traditional AIT well-site post-processing result cannot provide accurate indication of formation characterization.
- One aspect of the invention is directed to a method of deriving formation characterization. The method comprises receiving measurement data of a subsurface area; defining a measurement error distribution; creating multiple parameterized formation models of the subsurface area based on the measurement error distribution and the measurement data; generating multiple tool responses corresponding to the multiple parameterized formation models; considering the generated tool responses corresponding to the multiple parameterized formation modes; and deriving the formation characterization based on the result of considering generated tool responses.
- Embodiments of the method may further include performing borehole correction of the measurement data on the multiple tool responses and performing vertical processing on the borehole corrected measurement data. The method may further comprise allowing a user to choose the measurement error distribution. The method may further comprise creating multiple parameterized formation models by using Monte Carlo sampling.
- Another aspect of the invention is directed to an apparatus for deriving formation characterization. The apparatus comprises means for receiving measurement data of a subsurface area; means for defining a measurement error distribution; means for creating multiple parameterized formation models of the subsurface area based on the measurement error distribution and the measurement data; means for generating multiple tool responses corresponding to the multiple parameterized formation models; means for considering the generated tool responses corresponding to the multiple parameterized formation models; and means for deriving the formation characterization based on the result of considering generated tool responses.
- Additional objects and advantages of the invention will become apparent to those skilled in the art upon reference to the detailed description taken in conjunction with the provided figures.
- The present invention is illustrated by way of example and not intended to be limited by the figures of the accompanying drawings in which like references indicate similar elements and in which:
-
FIG. 1 is a flowchart showing steps associated with the present method, apparatus, and article of manufacture; -
FIG. 2 is a schematic illustration of computer hardware associated with the apparatus and article of manufacture; -
FIGS. 3-4 are a related group of diagrams used to describe the experimental result based on the inventive method; and -
FIGS. 5-7 are another related group of diagrams used to describe the experimental result based on the inventive method. - This invention is not limited in its application to the details of construction and the arrangement of components set forth in the following description or illustrated in the drawings. The invention is capable of other embodiments and of being practiced or of being carried out in various ways. Also, the phraseology and terminology used herein is for the purpose of description and should not be regarded as limiting. The use of “including,” “comprising,” or “having,” “containing”, “involving”, and variations thereof herein, is meant to encompass the items listed thereafter and equivalents thereof as well as additional items.
-
FIG. 1 shows several steps associated with the present method, apparatus and article of manufacture and provides a general overview of the invention. In the ReceiveMeasurement Data Step 10, the Measurement Data of the subsurface area is received. In one embodiment, the received Measurement Data israw data 12 from arrays of Array Induction Imager (AIT) family tools. The AIT-family tools typically include AIT, PlatformExpress AIT (AIT-H), SlimAccess AIT (SAIT), Xtreme AIT (XAIT), and SlimXtreme AIT (QAIT), etc. In another embodiment, the received Measurement Data israw data 12 from arrays of High Resolution Laterolog Array (HRLA) tool. Still in another embodiment, the received Measurement Data can be any measurement data (not shown in the drawing) of the subsurface area from a resistivity tool or other tools. - In the Define Measurement Error Distribution Step 14, a user/
engineer 16 is allowed here to define the data measurement error distribution. As the induction-logging instruments are susceptible to a variety of environmental effects due to considerable range of investigation, the data measurement error distribution can be very difficult to predict. An expert, with data processing experience, however, can make a reasonable assumption of the data measurement error distribution. A new user can make a wide guess of the data measurement error distribution. In one embodiment, the measurement error distribution is defined based on Gaussian error assumption. The Gaussian error assumption assumes that the error distribution is characterized by a mean value (most possible value in large amount of repeated loggings) and standard deviation σ (the possibility of the measured value drops to 63.2%). The possibility density function is
In another embodiment, the measurement error distribution is defined based on homogeneous error assumption. The homogeneous error assumption assumes that the possibilities of measured values are the same for certainty resistivity range. The possibility density function is p(x)=cons tan t. Still in another embodiment, the measurement error distribution is defined based on Laplacian error assumption. Specifically, the Laplacian distribution is given by its probability density function
In fact, theuser 16 can choose any existing mathematical functions to define the measurement error distribution, or even create his/own functions to define the measurement error distribution. In addition, theuser 16 can apply different measurement error distribution definitions and may reach different results of the formation characterization after the following steps, and the user can choose one of the results that appears the most accurate one by comparing the different results of the formation characterization among each other. - In the Create
Formation Models Step 18, multiple parameterized formation models of the subsurface area are created based on the measurement error distribution and the measurement data. Typically the parameterized formation models include a plurality of formation and wellbore related parameters, e.g. formation resistivity, invasion diameter, borehole size, mud resistivity, tool stand-off, Young's modulus, Poisson coefficient, shear modulus etc. In one embodiment, a large amount of parameterized formation models are created using Monte Carlo sampling, i.e. after one parameterized formation model is created, the parameterized formation model is accepted with the chance defined by the possibility density function. When large amount of parameterized formation models are created, they map out the user defined possibility density function. In one embodiment, the amount of the parameterized formation models can be 1,000 or even more depending on the user. - In the Generate
Tool Responses Step 20, multiple tool responses are generated corresponding to the multiple parameterized formation models created in the CreateFormation Models Step 18. In one embodiment, multiple resistivity tool responses are generated by simulating the tools response from a forward modeling engine. - In one embodiment of the invention, the AIT-family log processing includes borehole correction as stated in the Perform Borehole Correction
Step 22. The raw data is corrected from all eight arrays for the borehole effect. In one embodiment, borehole corrections for the AIT tools are based on inversion through an iterative forward model to find the borehole parameters that best reproduce the logs from the four shortest arrays—the 6-, 9-, 12-, and 15-in arrays. The borehole forward model is based on a solution to Maxwell's equations in a cylindrical borehole with resistivity Rm surrounded by a homogeneous formation of resistivity Rf. - In one embodiment of the invention, the AIT-family log processing also includes vertical processing as stated in the Vertical Processing Step 24. The borehole corrected data are going through a non-linear procedure to compensate for skin-effect and shoulder effect.
- In the Consider Generated Tool Responses Step 26, the generated tool response behavior corresponding to the multiple parameterized formation models is examined. In one embodiment, the generated tool responses are averaged. In another embodiment, the percentage of certain characters of the generated tool response behavior is examined. The user can define the threshold of percentage, e.g. 60% or 80%, for concluding the existence of certain character.
- Finally, in the Derive
Formation Characterization Step 28, the formation characterization is derived based on the examination results of the generated tool response behavior. In one embodiment, theuser 16 can define the percentage threshold for certain formation character. For example, the user may define 80% of the results showing curve separation as indication of a formation bed existence. From the statistics of the Monte Carlo simulated and processed results, if the curve separation is equal or more than 80%, the user can conclude that there exists a formation bed. - Therefore, an embodiment of the inventive method is disclosed. The inventive method allows a user/engineer to define the data measurement error distribution. With the help of this inventive method, AIT response distribution behavior under a large number of formation model and borehole situations can be examined. The inventive method does not reply on the Gaussian error assumption so that when this assumption is violated, the inventive method can give a more reasonable result than traditional well-site post-processing method. In addition, the user can choose different measurement error distribution definitions and may reach different results of the formation characterization. Thus, the user can choose one of the results that appears the most accurate one by comparing the different results of the formation characterization among each other. Further, the user can choose one of the results that appears the most accurate one by comparing the different results of the formation characterization with the results obtained by other means or tools.
- There are multiple advantages in applying the inventive method for well-site post-processing. For example, the inventive method can be helpful in formation uncertainty study and answering questions like: whether the inversion is sensitive enough for that parameter; whether there are other possibilities of formation models; whether more information is needed in drawing conclusions, etc. Thus the inventive method can help us to benefit more from traditional AIT measurement interpretation.
-
FIG. 2 schematically illustrates computer hardware that may be used to implement the inventive method. Computer 30 has a media reading device, such as a CD-ROM Reader 32, a floppy disk device, or a ZIP drive. The media reading device may also be capable of recording the output of the program the computer 30 is running. A user of the computer 30 may enter commands using a user input device, such as akeyboard 34 or a mouse, may view output of the program code on a visual display device, such asmonitor 36, and may make hardcopies of output using an output device, such asprinter 38. When properly configured, computer 30 (and its associated peripheral devices) is an apparatus for deriving formation characterization in accordance with the present invention. Computer media, such as a CD-ROM 40, a floppy disk, or a ZIP disk, may have computer readable program code that allows the computer 30 to derive formation characterization in accordance with the inventive method. - We now provide some experimental results of applying the inventive method to derive the formation characterization. A few of case studies are provided, which yields further insights into this invention.
-
FIG. 3 shows the comparison results between well-site AIT post-processing result based on Gaussian error assumption and based on the inventive method. Specifically, formation anisotropy and vertical thin bed detection of a borehole and formation in La Copital field, Texas, US are examined. On the left, well-site AIT post-processing result, based on traditional Gaussian error assumption, indicates five curves of depth ofinvestigation FIG. 4 (in which Rh stands for horizontal resistivity and Rv stands for vertical resistivity). InFIG. 4 fourth track, R*72H stands for inversed formation model (Rh and Rv) result from 72-in spacing high frequency and fifth track; DIP and AZIM stands for inversed dip and azimuth results for the specific spacing and frequency. Therefore, well-site AIT post-processing result based on the inventive method provides more accurate results of the formation characterization than that based on Gaussian error assumption. - Further example for formation water estimation is illustrated below. Water based mud is used and mud filtrate resistivity is around 0.064 ohm·m. As shown in
FIG. 5 , in the depth range 3240-3248 m, GR (a nuclear tool that measures formation natural gamma ray radiation) indicates a clean zone. HRLA five curves overlap each other, showing a water zone with more or less constant resistivity Rt of 3 ohm·m. AIT was logged one day after HRLA logging, but the traditional AIT logging shows strange curve separation, which makes it difficult to estimate formation water resistivity, as shown inFIG. 6 . Specifically, AIT shows obvious resistive invasion, where the deep curve (AT90) indicates a more conductive and variable formation water with resistivity ranging from 0.5 ohm·m to 2 ohm·m. Since the mud type is very conductive, by the time AIT was logged, conductive invasion can be strong and un-favorable to AIT measurement. If HRLA provides a better estimation in formation water, AIT deep resistivity curve measurement will give an estimation of formation water resistivity only 33% of what has been estimated from HRLA measurements. - By using the inventive method, the discrepancy can be greatly eliminated. As shown in
FIG. 7 , by applying the inventive method, the uncertainty method result agrees quite well with HRLA measurement. The AIT deep curve indicates a water zone with resistivity around 3 ohm·m. Also the Monte Carlo large sample results show that quite possibly there was something wrong of the AIT hardware during logging. When there is alternative measurement, like HRLA, it is better using the alternative one. - Having thus described several aspects of at least one embodiment of this invention, it is to be appreciated various alterations, modifications, and improvements will readily occur to those skilled in the art. Such alterations, modifications, and improvements are intended to be part of this disclosure, and are intended to be within the spirit and scope of the invention. Accordingly, the foregoing description and drawings are by way of example only.
Claims (20)
1. A method of deriving formation characterization comprising:
receiving measurement data of a subsurface area;
defining a measurement error distribution;
creating multiple parameterized formation models of the subsurface area based on said measurement error distribution and said measurement data;
generating multiple tool responses corresponding to said multiple parameterized formation models;
considering the generated tool responses corresponding to said multiple parameterized formation models; and
deriving the formation characterization based on the result of considering generated tool responses.
2. The method of claim 1 further comprising performing borehole correction of said measurement data on said multiple tool responses.
3. The method of claim 2 further comprising performing vertical processing on said borehole corrected measurement data.
4. The method of claim 1 , wherein defining said measurement error distribution is based on Gaussian error assumption.
5. The method of claim 1 , wherein defining said measurement error distribution is based on homogeneous error assumption.
6. The method of claim 1 , wherein defining said measurement error distribution is based on Laplacian error assumption.
7. The method of claim 1 , wherein defining said measurement error distribution allows a user to choose said measurement error distribution.
8. The method of claim 1 , wherein creating multiple parameterized formation models comprises functions computed using Monte Carlo sampling.
9. The method of claim 1 , wherein said multiple parameterized formation models include a large number of formation models.
10. The method of claim 1 , wherein generating multiple tool responses includes generating resistivity tool responses.
11. The method of claim 1 , wherein generating multiple tool responses includes generating Array Induction Imager (AIT) or AIT-family tool responses.
12. The method of claim 1 , wherein generating multiple tool responses includes generating High Resolution Laterolog Array (HRLA) tool responses.
13. The method of claim 1 , wherein considering said measurement data includes using probability functions.
14. The method of claim 1 , wherein considering said measurement data includes using average functions.
15. The method of claim 1 further comprising comparing the formation characterizations derived from multiple defined measurement error distributions among each other.
16. The method of claim 1 , wherein deriving the formation characterization includes defining a percentage threshold for a formation character by a user.
17. An apparatus for deriving formation characterization comprising:
means for receiving measurement data of a subsurface area;
means for defining a measurement error distribution;
means for creating multiple parameterized formation models of the subsurface area based on said measurement error distribution and said measurement data;
means for generating multiple tool responses corresponding to said multiple parameterized formation models;
means for considering the generated tool responses corresponding to said multiple parameterized formation models; and
means for deriving the formation characterization based on the result of considering generated tool responses.
18. The apparatus of claim 17 further comprising means for performing borehole correction of said measurement data on said multiple tool responses.
19. The apparatus of claim 18 further comprising means for performing vertical processing on said borehole corrected measurement data.
20. An article of manufacture, comprising:
a computer usable medium having a computer readable program code means embodied therein for deriving formation characterization, the computer readable program code means in said article of manufacture comprising:
computer-readable program means for receiving measurement data of a subsurface area;
computer-readable program means for defining a measurement error distribution;
computer-readable program means for creating multiple parameterized formation models of the subsurface area based on said measurement error distribution and said measurement data;
computer-readable program means for generating multiple tool responses corresponding to said multiple parameterized formation models;
computer-readable program means for considering the generated tool responses corresponding to said multiple parameterized formation models; and
computer-readable program means for deriving the formation characterization based on the result of considering generated tool responses.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/279,409 US20070244646A1 (en) | 2006-04-12 | 2006-04-12 | Method of formation characterication using tool responses |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US11/279,409 US20070244646A1 (en) | 2006-04-12 | 2006-04-12 | Method of formation characterication using tool responses |
Publications (1)
Publication Number | Publication Date |
---|---|
US20070244646A1 true US20070244646A1 (en) | 2007-10-18 |
Family
ID=38605881
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/279,409 Abandoned US20070244646A1 (en) | 2006-04-12 | 2006-04-12 | Method of formation characterication using tool responses |
Country Status (1)
Country | Link |
---|---|
US (1) | US20070244646A1 (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009070384A1 (en) * | 2007-11-30 | 2009-06-04 | Exxonmobil Upstream Research Company | Method for reservoir fracture and cross beds detection using tri- axial/multi-component resistivity anisotropy measurements |
US20100213943A1 (en) * | 2009-02-23 | 2010-08-26 | Baker Hughes Incorporated | Method for accentuating signal from ahead of the bit |
US20130038463A1 (en) * | 2011-08-09 | 2013-02-14 | Denis Heliot | Interactive Display of Results Obtained from the Inversion of Logging Data |
US9366135B2 (en) | 2013-10-08 | 2016-06-14 | Exxonmobil Upstream Research Company | Automatic dip picking from wellbore azimuthal image logs |
US20170045641A1 (en) * | 2015-04-16 | 2017-02-16 | Halliburton Energy Services, Inc. | Correcting log data of array induction tools |
WO2017058692A1 (en) * | 2015-09-30 | 2017-04-06 | Schlumberger Technology Corporation | Methods and systems to analyze bed boundary detection |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4957172A (en) * | 1989-03-01 | 1990-09-18 | Patton Consulting, Inc. | Surveying method for locating target subterranean bodies |
US5675147A (en) * | 1996-01-22 | 1997-10-07 | Schlumberger Technology Corporation | System and method of petrophysical formation evaluation in heterogeneous formations |
US6366088B1 (en) * | 1999-12-03 | 2002-04-02 | Halliburton Energy Services, Inc. | Method to enhance vertical resolution of time-decay log using forward modeling deconvolution of time-decay spectra |
US6502036B2 (en) * | 2000-09-29 | 2002-12-31 | Baker Hughes Incorporated | 2-D inversion of multi-component induction logging data to resolve anisotropic resistivity structure |
US6646437B1 (en) * | 2000-04-07 | 2003-11-11 | Halliburton Energy Services, Inc. | System and method for clay typing using NMR-based porosity modeling |
US20050240349A1 (en) * | 2004-04-26 | 2005-10-27 | Schlumberger Technology Corporation | Multimodal decomposition |
-
2006
- 2006-04-12 US US11/279,409 patent/US20070244646A1/en not_active Abandoned
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4957172A (en) * | 1989-03-01 | 1990-09-18 | Patton Consulting, Inc. | Surveying method for locating target subterranean bodies |
US5675147A (en) * | 1996-01-22 | 1997-10-07 | Schlumberger Technology Corporation | System and method of petrophysical formation evaluation in heterogeneous formations |
US6366088B1 (en) * | 1999-12-03 | 2002-04-02 | Halliburton Energy Services, Inc. | Method to enhance vertical resolution of time-decay log using forward modeling deconvolution of time-decay spectra |
US6646437B1 (en) * | 2000-04-07 | 2003-11-11 | Halliburton Energy Services, Inc. | System and method for clay typing using NMR-based porosity modeling |
US6502036B2 (en) * | 2000-09-29 | 2002-12-31 | Baker Hughes Incorporated | 2-D inversion of multi-component induction logging data to resolve anisotropic resistivity structure |
US20050240349A1 (en) * | 2004-04-26 | 2005-10-27 | Schlumberger Technology Corporation | Multimodal decomposition |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009070384A1 (en) * | 2007-11-30 | 2009-06-04 | Exxonmobil Upstream Research Company | Method for reservoir fracture and cross beds detection using tri- axial/multi-component resistivity anisotropy measurements |
US20100230095A1 (en) * | 2007-11-30 | 2010-09-16 | Hezhu Yin | Method For Reservoir Fracture and Cross Beds Detection Using Tri-Axial/Multi-Component Resistivity Anisotropy Measurements |
US8360146B2 (en) | 2007-11-30 | 2013-01-29 | Exxonmobil Upstream Research Company | Method for reservoir fracture and cross beds detection using tri-axial/multi-component resistivity anisotropy measurements |
US20100213943A1 (en) * | 2009-02-23 | 2010-08-26 | Baker Hughes Incorporated | Method for accentuating signal from ahead of the bit |
US20130038463A1 (en) * | 2011-08-09 | 2013-02-14 | Denis Heliot | Interactive Display of Results Obtained from the Inversion of Logging Data |
US9846256B2 (en) * | 2011-08-09 | 2017-12-19 | Schlumberger Technology Corporation | Interactive display of results obtained from the inversion of logging data |
US9366135B2 (en) | 2013-10-08 | 2016-06-14 | Exxonmobil Upstream Research Company | Automatic dip picking from wellbore azimuthal image logs |
US20170045641A1 (en) * | 2015-04-16 | 2017-02-16 | Halliburton Energy Services, Inc. | Correcting log data of array induction tools |
WO2017058692A1 (en) * | 2015-09-30 | 2017-04-06 | Schlumberger Technology Corporation | Methods and systems to analyze bed boundary detection |
GB2556609A (en) * | 2015-09-30 | 2018-05-30 | Schlumberger Technology Bv | Methods and systems to analyze bed boundary detection |
GB2556609B (en) * | 2015-09-30 | 2021-07-14 | Schlumberger Technology Bv | Methods and systems to analyze bed boundary detection |
US11307323B2 (en) * | 2015-09-30 | 2022-04-19 | Schlumberger Technology Corporation | Methods and systems to analyze bed boundary detection |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP1328832B1 (en) | Method for determination of apparent resistivities of anisotropic reservoirs | |
US8775084B2 (en) | Adaptive borehole corrections accounting for eccentricity for array laterologs | |
CN104169524B (en) | The correction based on inverting of downhole electromagnetic instrument | |
US8515721B2 (en) | Method for integrated inversion determination of rock and fluid properties of earth formations | |
EP1953571B1 (en) | Nuclear tool used in a borehole to determine a property of the formation | |
US7705592B2 (en) | Two dimensional T1/T2APP-T2APP processing of multi-gradient NMR data | |
US10228484B2 (en) | Robust multi-dimensional inversion from wellbore NMR measurements | |
EP2912496B1 (en) | Identifying formation, matrix and fluid related characteristics from subsurface data using factor analysis | |
US20070244646A1 (en) | Method of formation characterication using tool responses | |
US20050116709A1 (en) | System and methods for upscaling petrophysical data | |
US10605952B2 (en) | Simultaneous inversion of NMR multiple echo trains and conventional logs | |
US7502691B2 (en) | Method and computer program product for determining a degree of similarity between well log data | |
Behroozmand et al. | A comprehensive study of parameter determination in a joint MRS and TEM data analysis scheme | |
WO2017192530A1 (en) | Evaluation of formation mechanical properties using magnetic resonance | |
Trevizan et al. | Method for predicting permeability of complex carbonate reservoirs using NMR logging measurements | |
Zhao et al. | The importance of fracture geometry and matrix data on transient hydraulic tomography in fractured rocks: Analyses of synthetic and laboratory rock block experiments | |
US10302801B2 (en) | Temperature correction of magnetic resonance data | |
CN113189672B (en) | Tunnel advance geological forecast method based on multi-attribute inversion | |
CN105556345A (en) | System and method for estimating porosity distribution in subterranean reservoirs | |
Tabanou et al. | Thinly laminated reservoir evaluation in oil-base mud: high resolution versus bulk anisotropy measurement-a comprehensive evaluation | |
RU2399070C2 (en) | Method for determination of residual water saturation and permeability of bed | |
Zhang et al. | Real time quasi–2-D inversion of array resistivity logging data using neural network | |
Deng et al. | Multi-well interpretation of wireline logs and core data in the Eagle Ford Shale using Bayesian Inversion | |
Bowers et al. | A probabilistic approach to determine uncertainty in calculated water saturation | |
Aykroyd et al. | 3D modelling and depth estimation in archaeological geophysics |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: SCHLUMBERGER TECHNOLOGY CORPORATION, TEXAS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:ZHANG, TIANHUA;STOWE, IAN DAVID;REEL/FRAME:017470/0218;SIGNING DATES FROM 20060330 TO 20060404 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |