GB2581062A - Improved borehole log data processing methods - Google Patents

Abstract

A method of processing borehole log data creates one or more image logs. Borehole log data is modelled as components (e.g. illuminance and surface reflectance) of an image. Multiple versions of the image are created having respectively differing static scales. Multiple image versions are combined to create a composite image. The log data may be resistivity, acoustic, or nuclear log data. Histogram equalisation may be performed. The method may include energising a geological formation using a logging sonde, generating signals from the returned energy, and transmitting or storing the signals as the borehole log data. A computer programmed to perform the method may be housed within the sonde or connecte as part of a logging toolstring including a sonde.

Description

IMPROVED BOREHOLE LOG DATA PROCESSING METHODS

The invention relates to improved borehole log data processing methods.

The logging of geological formations is, as is well known, economically a highly important activity. The invention is of benefit in logging activities potentially in all kinds of mining and especially in the logging of reserves of oil and gas, water or other valuable commodities.

Virtually all commodities used by mankind are either farmed on the one hand or are mined or otherwise extracted from the ground on the other, with the extraction of materials from the ground providing by far the greater proportion of the goods used by humans.

It is extremely important for an entity wishing to extract materials from beneath the ground to have as good an understanding as possible of the lithology of a region from which extraction is to take place.

This is desirable partly so that an assessment can be made of the quantity and quality, and hence the value, of the materials in question; and also because it is important to know whether the extraction of such materials is likely to be problematic.

In consequence a wide variety of logging methods has been developed over the years. The logging techniques exploit physical and chemical properties of a formation usually through the use of a logging tool or sonde that is lowered into a borehole (that typically is, but need not be, a wellbore) formed in the formation by drilling.

Typically the sonde sends energy into the formation and detects the energy returned to it that has been altered in some way by the formation. The nature of any such alteration is processed into electrical signals that are then used to generate logs (i.e. numerical, graphical or tabular representations containing much data about the formation in question).

The signals may be thought of as log data, or at least as representing log data. These data may be stored, processed and transmitted as numerical values e.g. in a database or other data memory format; or they may be processed to form graphical logs referred to as image logs.

In general the values of the log data when converted to forms that are useable in image logs are represented by features of the image logs that a skilled analyst can interpret by eye in order to assess the conditions in the donwhole (i.e. subterranean) environment.

Taking the specific example of resistivity log data, when these are converted to image format the data values are mapped to colours that may then be viewed e.g. on a computer screen or in a coloured print-out. It has been found that the use of colours in this way is somewhat intuitive for human analysts to study. In a resistivity image log that has not undergone any processing that alters the scale from place to place in the image therefore the analyst can know with confidence that all regions of the image that are the same colour represent regions, of the formation, that are of the same resistivity.

Log data from other types of logging tool may be similarly processed in order to give rise to other kinds of image log.

In short therefore an image log appears as lines and regions of colour that represent the shapes and distributions of subterranean features.

It is important in this regard to recognise that the logging sondes do not themselves either illuminate the subterranean formations or indeed capture light-based images. Rather, the log data (that by reason of resulting from penetration by the emitted energy of the geological features surrounding the borehole may relate to locations spaced some distance into the rock from the borehole) may be processed into images in order to facilitate interpretation of the data.

A topic of growing interest is that of machine assessment of log data. In particular the Applicant has pioneered numerous techniques that assist to automate processes of image log interpretation, thereby saving significantly on costs and time. In particular the Applicant has developed highly effective techniques that render image logs suitable for automated analysis; and also has developed highly effective automatic feature recognition and analysis methods.

Despite these advances however some aspects of log data and image logs derived from them remain sub-optimal, as explained below. The methods of the invention seek to improve the processing of image log data so that they more readily can be interpreted (whether by automated machinery such as computers, or by humans).

In part the difficulties derive from techniques of normalisation (i.e. the assignment of the log data values to the range of colours available for inclusion in an image log) that are currently in use.

Normalisation in the context of image log processing is the transformation of the measured values to a colour scale In respect of micro-resistivity images the convention is to represent low resistivity values by black and high resistivity values by white, with a gradational "earth-tones" colour scale between. Typical scales comprise 64 or 128 colours (the human eye having limited ability to differentiate between adjacent colours on gradational scales of more than 128 colours).

Normalisation techniques fall broadly into two categories: static normalisation and dynamic normalisation.

In static normalisation a fixed relationship exists between resistivity values and colour defined for the entire well. The desired maximum and minimum resistivity values for the scale are defined, and the interval is divided either linearly or logarithmically into the number of bins in the colour scale. Static scales allow the gross similarities between intervals to be readily identified. The disadvantage is that the resistivity of rocks spans several orders of magnitude, so small but significant variations in resistivity within a narrow band are not rendered with this method.

As the name implies, in dynamic normalisation a dynamic relationship exists between resistivity values and colours in which the colour mapping applied to each depth frame is driven by the range of resistivity values in a depth window centred on the depth. The depth window is typically of the order of 1 m. All 64 (or 128) colours are used within each window, so the small but significant variations in resistivity that are lost in static images are rendered in a dynamic image. The disadvantage is that the same colour from different parts of the same well will not in general correspond to the same resistivity value.

In more detail, dynamic normalization may be implemented (for example) by computing the mean and standard deviation of the resistivity values in a window, then mapping the colour range to the resistivity range (mean -n x standard deviation) to (mean + n x standard deviation). A common alternative dynamic normalization (referred to as equalization) redistributes resistivity values within the depth window so as to achieve a linear cumulative distribution of normalized resitivity values.

A further aspect of dynamic normalisation is that none of the prior art dynamic normalization methods takes account of the spatial distribution of the resistivity values within the window. This means that the colour mapping applied to the depth frame at the central depth in the window is influenced by resistivity values elsewhere in the window, and this has the potential to create ghosts when the window contains high-contrast transitions. As used herein a "ghost" is a region of an image log created using one of the artificial log data to colour mapping techniques in which one colour may inaccurately spread into a region of the image that does not exhibit the log data value in question. Ghosting therefore leads to significant inaccuracies, and the loss of some features from the image logs altogether.

It is important to recognise some further features of image logs that while not necessarily disadvantages must be borne in mind when considering them.

Firstly image logs are usually coloured, with (as explained) the colours signifying data values such as micro-resistivity values that either apply over the entirety of an image log (in the case of a statically normalised log) or within respective depth windows (in the case of dynamically normalised logs).

Image logs however can with some success be rendered in greyscale. The images herein are so rendered but it is to be understood that the method of the invention is not limited to the production of greyscale images and instead is intended to be applicable to, and indeed is primarily applicable to, coloured image logs.

Secondly if it is required to produce logs that are suitable for automated (machine) processing it is desirable to carry out some form of pre-processing on the image logs to render them suitable for this use.

An example of the kind of pre-processing that the Applicant has developed is so-called "in-painting".

In-painting is a technique that is applicable to resistivity logs, and to logs created from multiple sub-components that need not necessarily be logged at one and the same time.

An example of a logging tool type is the so-called multi-pad micro-resistivity borehole imaging tool, such as the tool 10 illustrated in transversely sectioned view in Figure 1. In this logging tool a circular array of (in the example illustrated) eight pads 11 each in turn supporting typically two lines of surface-mounted resistivity electrodes referred to as "buttons" 12 is supported on a series of calliper arms 13 emanating from a central cylinder 14. During use of the tool 10 the arms 13 press the buttons 12 into contact with the very approximately cylindrical wall of a borehole. The borehole is normally filled with a fluid (such as a water-based mud) that if conductive provides an electrical conduction path from the formation surrounding the borehole to the buttons.

Many variants on the basic imaging tool design shown are known. In some more or fewer of the pads 11 may be present. The numbers and patterns of the buttons 12 may vary and the support arms also may be of differing designs in order to achieve particular performance effects. Sometimes the designers of the tools aim to create e.g. two parallel rows of buttons located on the pad one above the other. The buttons in the lower row are offset slightly to one side relative to their counterparts in the row above. When as described below the signals generated by the buttons are processed the outputs of the two rows of buttons are in effect lain over one another. As a result the circumferential portion of the borehole over which the buttons 12 of a pad 11 extend is logged as though there exists a single, continuous, elongate electrode extending over the length in question.

In general in operation of a tool such as resistivity tools 10 electrical current generated by electronic components contained within the cylinder 14 spreads into the rock and passes through it before returning to the pads 11. The returning current induces electrical signals in the buttons 12.

Changes in the current after passing through the rock may be used to generate measures of the resistivity or conductivity of the rock. The resistivity data may be processed according to known techniques in order to create (typically coloured) image logs that reflect the composition of the solid and fluid parts of the rock. These image logs convey much data to geologists and others having the task of visually inspecting and computationally analysing them in order to obtain information about the subterranean formations.

In use of a tool such as that shown in Figure 1 the tool is initially conveyed to a chosen depth in the borehole before logging operations commence. The deployed location may be many thousands or tens of thousands of feet typically but not necessarily below, and in any event separated by the rock of the formation from, a surface location at which the borehole terminates.

Various means for deploying the tools are well known in the mining and oil and gas industries. One characteristic of most if not all of them is that they can cause a logging tool that has been deployed as aforesaid to be drawn from the deployed location deep in the borehole back towards the surface location. During such movement of the tool it logs the formation, usually continuously. As a result the image logs may extend continuously for great distances.

Although the logs are continuous in the longitudinal sense, notwithstanding the pad offsetting explained above they are azimuthally interrupted by reason of the pads not extending all the way continuously around the circumference of the borehole. The design of the tool prevents this since the arms 13 must be extensible in order to press the pads 11 into contact with the borehole wall. Following extension of the arms there exists a series of gaps between the ends of the pads.

No data can be logged in these gaps, which manifest themselves as elongate spaces in the image logs. Examples of image logs including several of these gaps or discontinuities are visible in Figures 6 and 7a. The discontinuities extend from one end of the reconstructed image log to the other, a distance in some cases of thousands of feet.

Filling in the missing data is advantageous for obvious reasons of the desirability of completeness of information. Moreover it is likely to be required when it is desired to process the image logs using automatic pattern recognition programs in order to try and identify certain features in the logs.

Patent application no. GB 1210533.4 discloses an in-painting technique, the purpose of which is to fill in, in a realistic manner, the discontinuities that result from use of a resistivity imaging tool of the kind illustrated in Figure 1.

Such in-painting is an example of the kind of pre-processing technique that may be applied to image log data.

The method of the invention is applicable to log data that has undergone pre-processing such as that outlined above, or perhaps another kind of pre-processing. Equally, however, the methods of the invention are applicable to, and useful in respect of, log data that have not undergone pre-processing.

Furthermore the method of the invention is applicable both to image logs that are intended for analysis by humans and those that are intended for machine analysis. For the avoidance of doubt the human and machine analysis approaches often are not mutually exclusive in one and the same log.

Yet further for the avoidance of doubt, the method of the invention extends to the processing of data in real time as it is generated for example using an imaging tool as described above (or another type of logging tool). It also extends to the processing of data the generation of which was not contemporaneous with the log data processing steps. Thus the methods of the invention are applicable to e.g. so-called logging while drilling (LWD) data or other real-time generated data; or data the production of which is not to contemporaneously associated with operation of a logging tool.

As stated below, the steps of operating a logging tool therefore optionally form part of the invention.

As disclosed but not claimed herein a method of processing borehole log data, to create one or more image logs, comprises the steps of: a. Modelling the data as components of an image in the form i (x, y) = 1(x, y) x r(x, y) (1) in which i(x, y) is an image representative of the log data; /(x, y) denotes an illumination value of the image at two-dimensional spatial co-ordinates x, y; and r(x, y) denotes a surface reflectance value at the co-ordinates x, y; b. Transforming Equation (1) to a logarithmic domain; c. Obtaining a Fourier transform of the resulting logarithmic domain expression to obtain a Fourier domain expression; d. High-pass filtering the Fourier domain expression; e. Obtaining an inverse Fourier transform of the resulting filtered Fourier domain expression; Performing an exponential operation on the result of Step e. to obtain a filtered image model expression; g. Mapping values of the filtered image model expression to respective colour values across the range of the filtered image log expression values; and h. Displaying, printing, saving and/or transmitting the mapped colour values as one or more image logs.

Such a method advantageously separates in the logarithm domain the illumination values and reflectance values in the log data thereby permitting selective filtering in order to provide an enhanced log in which the effects of the log data components that are modelled as illumination values are reduced. In such a log the dynamic range of the image log is extended compared with a normally processed image log, so that more information may be rendered using a fixed colour scale.

This disclosed method may be carried out in respect of all the data contributing to an image log, in which case the method is a substitute for static normalisation as described above. In such a case the background colours of the image generally represent the same resistivity values across all the log data (which may correspond to many feet of logged well), with dynamic normalisation characteristics being apparent only in respect of thin, small-scale geological features. Thus the disclosed method permits the generation of image logs that are essentially similar to statically normalised ones, but with enhanced detail and information in areas corresponding to particularly fine or otherwise difficult to resolve features.

This simulated combination of static and dynamic normalisation characteristics in one and the same image log is likely to be of very significant benefit in the interpreting of image logs.

Alternatively the method may be applied on a sliding window basis, which concept is broadly familiar to those of skill in the image logging art. Thus the log data may relate to a depth window corresponding to less than the whole depth of the logged data, and the method includes repeating the steps of any preceding claim in respect to at least one further depth window in the logged data.

When so practised the disclosed method optionally may include the step of analysing the spatial distribution of log values within each said depth window, and excluding from the data processed in accordance with Steps a. to h. data identified as being of high contrast.

The disclosed method may be termed "dynamorphic filtering", which is intended to convey both the notion of filtering of image log data to remove the illumination value component and the concept of dynamic normalisation at places in which the geological features would otherwise be difficult to resolve. The extent to which the dynamic aspect is applied occurs automatically as a result of the algorithm used.

As further disclosed but not claimed herein a method of processing borehole log data, to create one or more image logs, comprises the steps of: Modelling the log data as components of an image; k. Specifying positional parameters of a virtual light source as azimuthal and elevation angles; Selecting the set of pixels corresponding to the azimuthal angle: m. Computing a one-dimensional derivative of a sub-set of pixels corresponding to this orientation, in which the sub-set is selected based on the elevation angle; n. Shading the image as though the light source illuminates the image at the elevation and azimuthal angles; and o. Displaying, printing, saving and/or transmitting the shaded image.

This method may be carried out on its own, or in combination with the steps of the first method disclosed above. In either case an enhanced image log results.

According to the invention there is provided a method of processing borehole log data, to create one or more image logs, comprising the steps of modelling the data as components of an image; creating multiple versions of the image having respectively differing static scales; and combining the multiple image versions to create a composite image.

Thus the invention amounts to the application of techniques akin to those known in the image processing art as high dynamic range (HDR) processing to image log data. Further enhancements to the qualities of image logs result from the application of this aspect of the invention.

Preferably the log data are resistivity log data and/or acoustic log data (especially ultrasound log data) and/or nuclear log data.

Also preferably the method of the invention additionally includes histogram equalisation.

Conveniently the method of the invention includes the steps of operating a logging sonde to energise a geological formation; receiving returned energy at the sonde; from the returned energy generating one or more signals; and transmitting or storing the signals as the borehole log data.

The processing steps of the method may be carried out using a programmed computer or processor.

In addition to the foregoing the invention is considered to reside in one or more of: a computer or processor that is programmed to perform the processing steps of a method in accordance with any of Claims 1 to 5 hereof; such a computer or processor when operatively connected to a logging sonde; such a computer or processor when housed within a logging sonde, or connected as part of a logging toolstring including a sonde; log data processed in accordance with the methods hereof; and/or image logs produced in accordance with the methods hereof.

There now follows a description of preferred embodiments of the invention, by way of non-limiting example, with reference being made to the accompanying representations in which: Figure 1 is a schematic, partly sectioned representation of a prior art micro-resistivity logging tool; Figure 2 is a flow chart summarising the method of the first aspect of the invention; 10 15 20 25 Figure 3 is a schematic representation of the steps of an optional windowing process used when applying the method of the first aspect of the invention to resistivity image log data; Figure 4 shows in greyscale a prior art statically normalised resistivity image log; Figure 5 shows in greyscale a resistivity image log that is the result of processing the log data of Figure 4 in accordance with the method of the first aspect of the invention; Figures 6a to 6e are greyscale representations of comparisons between exemplary outputs of the method of the invention with statically and dynamically normalised (prior art method) image logs; to Figures 7a to 7c are greyscale representations of comparisons between untreated image log data (Figure 7a), statically normalised, in-painted log data (Figure 7b) and log data (Figure 7c) that have been the subject of Virtual Light Source treatment in accordance with the method of the second aspect of the invention; Figures 8a to 8c are greyscale representations of comparisons between statically normalised, in-painted log data (Figure 8a), virtual light source treated data (Figure 8b), and log data (Figure 8c) that have been the subject of combined treatment in accordance with the method of both the first and the second aspects of the invention; Figure 9a illustrates in greyscale multiple images created as part of the method of the third aspect of the invention; Figure 9b shows the output of the method of the third aspect of the invention, resulting from combining the images in Figures 9a; and Figure 10 shows in schematic form an apparatus used for performing the steps of the methods of the invention.

For the avoidance of doubt where in this specification greyscale images are shown this is for ease of illustration. It would be expected in the majority of cases of use of the methods of the invention that coloured image logs would be the normal output.

Referring to the representations the methods of the invention are described in the following sections, in which the sub-headings indicate the various aspects of the invention.

Illumination-Reflectance Model An image can be considered as a 2D function of the form i(x,y), the value of which at spatial coordinates (x,y) is a positive scalar quantity the physical meaning of which is determined by the source of the image. In the case of greyscale images, when an image is generated from a physical process its values are proportional to energy radiated by a physical source. In other words, an image is an array of measured light intensities that are a function of the amount of light reflected off the objects in the scene.

The intensity is a product of illumination (the amount of source illumination incident on the scene being viewed) and reflectance (the amount of illumination reflected by the objects in the scene). Note that this is still valid for resistivity log images, even though as explained they are not optical images. Denoting illumination as 1(x, y) and surface reflectance as r(x, y), then an image i(x, y) can be expressed as: i(x, y) = 1(x, y) x r(x, y) (1) The model of image formation in Equation (1) is known as the illumination-reflectance model. It is a simplification of the general model in which the specular reflection, the directly visible light source, and caustics are ignored. This simplification does not hold in general, but it is useful for explaining the principles of the invention. The model can be used to address the problem of improving the quality of an image that has been acquired under poor illumination conditions, but in the context of the invention it is used, as explained herein, to improve the dynamic range of borehole resistivity images.

If i(x, y) is given it may be impossible to retrieve either of its constituent components, 1(x, y) or r(x, y). However, for certain applications (including tone reproduction) it may be desirable to separate surface reflectance from signal. Although this is generally an under-constrained problem, it is possible in accordance with the invention to transform Equation (1) to the log domain, where the multiplication of r(x, y) and 1(x, y) becomes an addition.

Then under specific conditions the two components could be separated.

An image represented in the logarithmic domain is referred to as a density image. The fact that the two entities are added in the logarithmic domain gives a direct result to operate in the log domain. Filtering operations such as tone reproduction may be carried out in this domain. Such processes are referred to herein as dynamorphic filtering.

Dvnamorphic Filtering Concept The dynamorphic concept combines the illumination-reflectance model with the concept of the ghost-free dynamic normalization described above.

Illumination results from the conditions present when the image is captured, and can change when these conditions change (for example, when the borehole environment changes). Illumination variations can be thought of as multiplicative noise, and can be reduced by filtering in the log domain. However, reflectance results from the way the objects in the image reflect light (in the case of a resistivity log it is current flowing in the rock), and is determined by the intrinsic properties of the object itself, which (in this theoretical analysis) do not change.

One may consider that illumination varies slowly in space (slow spatial changes low spatial frequency) while reflectance changes abruptly (high spatial frequencies). When seeking to eliminate apparent changes in the resistivity image appearance with changes in lighting conditions, it is desirable to enhance the reflectance while reducing the contribution of illumination. Hence, it is desirable to separate the two components of Equation (1) and then high pass filter the resulting image in the frequency domain.

Dynamorphic filtering as defined herein is a frequency domain filtering process that achieves this objective by transforming the expression in Equation (1) from multiplication to addition, the problem of high pass filtering then being made trivial as it becomes possible to use the multiplication or convolution property of the Fourier transform 3.

A solution to this problem, within the scope of the invention, is to take a natural logarithm (base e) of both sides of Equation (1): z(x, y) = In(i(x, y) + 1) = In(i(x, y) x r(x, y) + 1) = 141(x, y)) + In (r(x, y)) (2) where the +1 is added to make sure the situation In(0) does not arise. Applying the Fourier transform to Equation (2) fi[z(x, y)] = 3{In(/(x, y))] + filln(r(x, y))] (3) or: Z(p, v) = L(p, v) + R (p, v) (4) where L(p, v) and R(p, v) are the Fourier transforms of In(/(x, y)) and In (r(x, y)) , respectively. Now it is possible to high pass filter Z(p, v) by means of a filter function H(p, v) in the frequency domain and obtain a filtered version S(p, v): = 0 x Z04, 0 = Km, 0 x 1,(p, + H (p, x R (p, (5) Taking an inverse Fourier transform of Equation (5) provides: s(x, y) = 3-1fH (p, v) x L(je, 32)} + 3-1[H(tt, v) X R(u, (6) and finally, one may achieve the desired filtered (enhanced) image t(x,y) by the exponential operation: i(x, y) = es(x -1 (7) The preferred high pass filter normally used in this procedure is the Butterworth filter defined as: H(R, v) - (8) i+r Do int [D(ii.v).1 where n defines the order of the filter. Do is the cut-off distance from the centre and D(A, v) is given by: D GA = 2)2 v 2)2 1/2 (9) where M and N are the number of rows and columns of the original image, respectively. The whole process is summarized in Figure 2.

In the case of borehole resistivity images, the dynamorphic filtering process preferably is followed by a classical histogram equalisation for further improvement.

Enhancement algorithms such as histogram equalisation and dynamorphic filtering as defined and claimed herein are global in nature and are intended to enhance an image and deal with it as a whole. For resistivity images, it is in the alternative possible within the scope of the invention to apply the two algorithms on a windowing basis by splitting the original image in sub-images and filtering each sub-image individually.

In summary, dynamorphic normalization (or dynamic range improvement) seeks to increase the amount of useful information that can be rendered by a fixed colour scale. It may be applied to the whole well, or on a sliding window basis. When applied to the whole well it is an alternative to the prior art conventional static image normalisation. When applied on a window basis, it is an alternative to conventional dynamic normalization and/or dynamic equalization.

The sliding window approach is illustrated schematically in Figure 3, and may be summarised as follows: The algorithm comprises the following steps: Select a depth window; Compute column means within that window; Select the middle value from the mean values vector; Retain lines in the window the means of which are within a defined range of the middle value of the vector of means.

Apply the dynamorphic filtering method described above on the retained data, keeping only the line corresponding to the centre line from the depth window; Go to the next depth line, and repeat; After completing for the whole file, normalise the result to 64 colours.

This set of steps has the beneficial effect of removing high-contrast data.

Implementation of Dynamorphic Filtering Additionally to Exclude Gradient Artefacts Conventional dynamic normalization produces spurious colour gradients adjacent to boundaries with high resistivity contrast. This is because high contrast boundaries cause a sudden change in the window-average resistivity that controls the normalization, and the change falls away only gradually away from the boundary. The method of the invention addresses this problem in the way it implements the Dynamorphic Normalization algorithm. Specifically the method considers the spatial distribution of resistivity values within the depth window, and excludes high contrast data from the average while applying the dynamorphic filtering within that window.

The artefact suppression is demonstrated in Figures 6a to 6e in which it is compared with the industry-standard dynamic normalization as presently implemented in existing software such as Well Imager provided by Weatherford International.

In each of Figures 6a to 6e the left hand track shows a resistivity image log that has been the subject of (prior art) static normalisation; the centre track the same log data following treatment using a prior art dynamic normalisation method; and the right hand track the log data having been processed in accordance with the method of the first aspect of the invention.

In the left hand and centre tracks the occlusions resulting from operation of a micro-resistivity logging tool as described above are apparent as vertical lines in which no data are present. In the right hand track there are no occlusions because the data have been subject to a process of filling in the missing data lines. As stated however the method of the invention is applicable to log data that includes the occlusions or has been the subject of a process having the aim of eliminating them.

The following remarks, that relate as indicated below to the respective images 6a to 6e, identify the defects of the images that are not the result of data processing in accordance with the method of the invention, and also explain how the method of the invention produces improved results: Figure 6a The dynamic image misrepresents the thin beds. The dynamorphic image matches the optimally scaled static in this interval while also matching differently scaled static images in other intervals.

Figure 6b The dynamorphic (inventive) image has more detail in the thinly bedded intervals (ignoring the residual speed correction artefacts).

Figure 6c The dynamic image has introduced beds that do not exist in reality (as confirmed in the static image). All the features in the dynamorphic image are real, and better defined.

Figure 6d The dynamic image introduces false "shadows" below high-contrast boundaries, and adds other detail that is not real.

Figure 6e The dynamic image corrupts the central bed and introduces other distortions. The dynamorphic image improves on the dynamic and static images.

In general the images in the right hand track of Figures 6a to 6e tend to appear similar to those of the left hand track; but this is illusory. In view of the variable application of a dynamic scaling effect that takes place as a result of the inventive method the right hand track logs in fact exhibit significantly greater detail and accuracy at locations where this is required.

The centre track of Figures 6a to 6e graphically illustrates the phenomenon of "ghosting", in which certain colour scales spread to occupy areas beyond those corresponding to the log (resistivity) values they are intended to represent. This in turn derives from the fact that in a dynamically normalised image log one colour scale may represent more than one log value. Very clearly the phenomenon of ghosting renders the log considerably harder for either a human analyst or a computer to interpret correctly than should be the case.

Virtual Light Source (VLSI Research has shown that the way data are visualized has a bearing on how the eye-brain system perceives the information. In particular it has been found that the brain perceives low and high spatial frequencies differently, and whereas colour is useful for displaying low-frequency variations in data, it is not well-suited to the display of high-frequency data. The latter is best displayed using variations in luminosity.

The method of the second aspect of the invention involves treating high spatial frequencies differently, and include a Virtual Light Source (VLS) algorithm which consists of: Specifying a "light source" position in terms of its azimuthal and elevation angles.

Selecting the set of pixels corresponding to the azimuthal angle.

Computing a one-dimensional derivative of a sub-set of pixels corresponding to this orientation, in which the sub-set is selected based on the elevation angle.

Outputting the shaded image within these azimuthal and elevation angles.

This creates a virtual landscape (relief) whose peaks and troughs reflect the amplitude of the underlying resistivity values in a similar way to that in which shining a light onto a landscape casts shadows and gives the result a 3D appearance as shown in Figures 7a to 7c. Figure 7a shows an image derived from untreated log data; Figure 7b shows a statically normalized in-painted result and Figure 7c shows a VLS image. The azimuth and the elevation angle of the light source control the nature of the shadows. Features of interest in the image may be optimally illuminated using particular azimuth and elevation angles, thereby providing a mechanism for the preferential selection of features as part of an automated feature recognition process.

Combination of Dvnamorphic and VLS (DVLS1 VLS emphasises high spatial frequencies while retaining maximal low spatial frequency information. However, in the implementation described above it does not necessarily preserve the colour polarity of the static image. In order to preserve the benefits of VLS and at the same time preserve the colour polarity it is possible within the invention to combine VLS and Dynamorphic Normalization. The result is shown in Figure 8.

High Dynamic Range (HDR) Processing High Dynamic Range processing in accordance with the third aspect of the invention enhances contrast in scenes containing large areas of both light and dark by creating multiple versions of the same image log, each with a different static scale, and then combining the individual images into a single composite as illustrated in Figure 9a (showing in greyscale the multiple image versions mentioned) and Figure 9b (showing the composite image log created from them). This process is simpler and faster than the method of the first aspect of the invention, and results are totally free of any artefacts associated with to dynamic normalisation. It is believed however that in the method of the third aspect of the invention colour control is less easy to achieve.

Figure 10 shows apparatus that may be used in the performance of the methods described herein.

In Figure 10 nuclear log data 72, resistivity log data 74 and/or other types of log data such as acoustic log data are fed to a computer or other programmable processing device 70 either directly from a logging tool or wireline connected thereto (as signified by numerals 72 and 74) or via a memory device or medium 76 (that may be associated with or part of a logging tool, or may be separate from it). The computer 70 may perform the steps of the methods described herein in order to produce the enhanced image logs described.

Overall the methods of the invention amount to very significant improvements, as explained, in the quality and usability of image logs, by enhancing the usefulness of the dynamically applied colour mapping, by adding in a virtual light source effect and by producing enhanced image logs through the use for the first time on log data of HDR processing techniques. The resulting improved logs may be more effectively assessed and analysed both by human analysts and by machine vision equipment.

The listing or discussion of an apparently prior-published document in this specification should not necessarily be taken as an acknowledgement that the document is part of the state of the art or is common general knowledge.

Claims (9)

1. CLAIMS1. A method of processing borehole log data, to create one or more image logs, comprising the steps of: a. modelling the data as components of an image; b. creating multiple versions of the image having respectively differing static scales; and c. combining the multiple image versions to create a composite image.
2. 2. A method according to Claim 1 wherein the log data are resistivity log data and/or acoustic log data and/or nuclear log data.
3. 3. A method according to Claim 2 wherein when the log data are are ultrasonic log data.
4. 4. A method according to any preceding claim additionally including histogram equalisation.
5. 5. A method according to any preceding claim including the steps of operating a logging sonde to energise a geological formation; receiving returned energy at the sonde; from the returned energy generating one or more signals; and transmitting or storing the signals as the borehole log data.
6. 6. A method according to any preceding claim the processing steps of which carried out using a programmed computer or processor.
7. 7. A computer or processor that is programmed to perform the processing steps of a method in accordance with any of Claims 'I to 5.
8. 8. A computer or processor according to Claim 7 when operatively connected to a logging sonde.
9. 9. A computer or processor according to Claim 8 when housed within a logging sonde or connected as part of a logging toolstring including a sonde.