1. Field of the Invention
The present invention relates generally to processing of geological data and more particularly to normalization of well log curves.
2. Description of the Related Art
The geology of a borehole may be explored using borehole logging tools. Measurements of various properties of the subsurface formations through which the borehole penetrates can be made, and a log of the measurements against depth can be produced. This borehole data is collected via a number of techniques including resistivity/conductivity measurements, ultrasound, NMR, and gamma radiation, for example. Generally, borehole data is analyzed by human interpreters in order to characterize a subsurface geological formation to allow decisions to be made regarding potential of the well or to determine information about the nature of the surrounding geologic area. Borehole data of this type may be used to replace or supplement the collection of cores for direct inspection.
As with any measurement, well logging is subject to various causes of noise. For example, sensor noise (including, e.g., shot noise), sensor drift and other effects can cause logs to appear different from a theoretical sensor response.
As a number of different wells are logged, accounting for noise can allow a more accurate comparison between logs, providing a more useful understanding of the subsurface geology. In particular, where automated analysis methods are to be applied to the well logs, such normalization can help to reduce incorrect conclusions resulting from, for example, a gain change of the sensor between logging operations.
Aspects of embodiments of the present invention provide a method of normalizing well log data relating to geology of a subsurface region, includes identifying a pair of markers defining a zone on each well log, determining a statistic for a property within the zone for each well log, fitting a smooth function to the statistics to create a first-iteration trend surface, determining, for each statistic, a deviation from the trend surface and eliminating those statistics that deviate from the trend surface by greater than a threshold value, fitting a smooth function to the remaining statistics to create a second-iteration trend surface, normalizing the property for each well log based on the second-iteration trend surface, and analyzing the normalized properties to determine information relating to the geology of the subsurface region.
Aspects of embodiments of the invention may include a computer-readable medium encoded with computer-executable instructions for performing the foregoing method or for controlling the foregoing system.
Aspects of embodiments of the invention may include a system configured and arranged to provide control of the system in accordance with the foregoing method. Such a system may incorporate, for example, a computer programmed to allow a user to control the device in accordance with the method, or other methods.
These and other objects, features, and characteristics of the present invention, as well as the methods of operation and functions of the related elements of structure and the combination of parts and economies of manufacture, will become more apparent upon consideration of the following description and the appended claims with reference to the accompanying drawings, all of which form a part of this specification, wherein like reference numerals designate corresponding parts in the various FIGS. It is to be expressly understood, however, that the drawings are for the purpose of illustration and description only and are not intended as a definition of the limits of the invention. As used in the specification and in the claims, the singular form of “a”, “an”, and “the” include plural referents unless the context clearly dictates otherwise.
In a given region of an exploration area, there may be a large number of wells drilled. Some wells are drilled with the hope that they will be production wells, others are for fluid injection, or for testing purposes. During, or after, drilling of each well, one or more well logs may be produced, measuring properties of the subsurface formations through which they are drilled.
By way of example, resistivity and gamma ray logs may be taken to provide analysts with information that they can interpret in order to determine information about conditions in the subsurface region. Gamma ray logs, for example, may be used in assigning facies classifications to layers of subsurface material, allowing experts to identify impermeable layers such as shale that may help to indicate locations where hydrocarbon resources are less likely to be found.
As will be appreciated, a number of additional common markers could be located and correlated among the logs, even without any additional pre-processing. On the other hand, it is also clear that, for example, that there is a scaling factor (or gain change) between log 10 and log 12. While the features measured for log 12 appear to be quite well correlated to the features measured for log 10, the amplitudes are considerably smaller. Likewise, the amplitudes of log 16 are considerably larger than either log 10 or log 12, again despite representing measurements of physically similar subsurface features.
A method of performing the normalization on the well logs of
The threshold can be predetermined, or can be calculated during the execution of the method. One example of a suitable threshold is one standard deviation from the mean. For 10-50-90 statistics, each may be evaluated against a respective one standard deviation threshold. Once the thresholds have been applied, a function is fit to the remaining statistics to create a second-iteration trend surface 110. In an embodiment, the fitted function is a smooth function. Finally, the property is normalized based on the second-iteration trend surface 112.
In an embodiment, the normalizing includes shifting the statistic for the property for each well to the trend surface. As will be appreciated, normalization in accordance with embodiments of the invention can reduce changes due to changes in gain and/or shifts (offsets) in output of the logging tool.
One example of a particular method of normalizing the well logs in accordance with an embodiment of the invention is to begin with well logs that include at least an interval of interest. While the logs may be derived from depths that are not exactly equal, they should generally overlap and include common geological features. In an embodiment, portions of the logs are selected by defining top and base markers for each well, thereby defining the zone of calibration.
Each well in the group should have a property that requires normalization. In an embodiment, the property is continuous and is normalizable by either shifting the mean value or shifting the mean and adjusting the standard deviation.
The method may include pre-defined or user-definable parameters including lateral smoothing scale, and maximum deviation. The lateral smoothing scale can be increased for cases in which the well-to-well variation in a mean value of the property is not considered to be reflective of true variations in the reservoir. For example, this might be true when an area-blanketing marine shale shows unexpected variations in log properties. The maximum deviation can be adjusted upward in order to reduce a number of wells that are considered to be outliers. This might be useful where there are relatively few wells with a large variation, such that removal of the outliers would leave a sample size too small to provide significant normalization.
In an example of an embodiment, first an average value of the property is calculated between the top and base marker for each well. A least square fit to the average values can be calculated at a series of spatial locations that can be connected to form a surface with smoothly varying property variations. One embodiment of this process is to use discrete smooth interpolation which is available in the commercial software package Gocad licensed by the vendor Paradigm Geophysical. The smoothness of the fit surface is a function of the input lateral smoothing scale. The larger this input parameter the smoother the fit surface. As the lateral smoothing scale gets larger and approaches infinity, the trend surface becomes a single value.
The average values of the property are used to create a first iteration trend surface using discrete smooth interpolation. Once the first iteration trend surface has been established, a deviation of each mean value of the property for each well is computed. Each well that has a mean value that deviates further than the maximum deviation is considered to be an outlier, and removed from the process prior to performing a second iteration.
All wells having average property values that deviate further than the input maximum deviation are removed prior to performing the second iteration. In the second iteration, the steps are the same, simply excluding the outlier wells. Once the second iteration trend surface is calculated, the well log mean is shifted to the trend surface by subtracting deviation from the trend and storing the new result as a normalized curve property.
A system for performing the method is schematically illustrated in
Once the logs are normalized, they may be manipulated or used as any typical well log might be. The normalized logs may be, for example, more easily tied using machine assisted or automated processing based on similarity of common structures.
Although the invention has been described in detail for the purpose of illustration based on what is currently considered to be the most practical and preferred embodiments, it is to be understood that such detail is solely for that purpose and that the invention is not limited to the disclosed embodiments, but, on the contrary, is intended to cover modifications and equivalent arrangements that are within the spirit and scope of the appended claims. For example, though reference is made herein to a computer, this may include a general purpose computer, a purpose-built computer, an ASIC programmed to execute the methods, a computer array or network, or other appropriate computing device. As a further example, it is to be understood that the present invention contemplates that, to the extent possible, one or more features of any embodiment can be combined with one or more features of any other embodiment.