METHOD OF AUTOMATIC CHARACTERIZATION AND REMOVAL OF PAD ARTIFACTS IN ULTRASONIC IMAGES OF WELLS

Abstract

The present invention presents a method of automatic characterization and removal of marks and pad artifacts in ultrasonic images of reservoir wells. The method demonstrates the effectiveness of automatic characterization of this noise and its removal by modeling a two-dimensional square wave signal, periodic in the angular axis of the image, and includes: the obtention of the average curve of the one-dimensional power spectrum of the well image for the automatic detection of the artifact noise frequency response peak; the derivation of the geometric parameters of the signal of the artifacts by means of the frequency peak estimated in the previous step; the automatic modeling of the signal of the artifact as a periodic square wave using the parameters obtained in the previous steps; the processing of the original image using the square wave model filter obtained in the previous step.

Claims

1. A method of automatic characterization and removal of pad artifacts in ultrasonic images of wells characterized in that it comprises the following steps: 1) Automatic obtention of the period value of the pad marks; 2) Modeling of pad artifacts; 3) Application of a filter; and 4) Obtaining the final processed image.

2. The method according to claim 1, characterized in that Step 1 is only applied to data with a complete azimuthal sweep, after sweeping 360° of the well wall.

3. The method according to claim 2, characterized in that it additionally obtains the magnitude of the one-dimensional spectrum of the Fourier transform of each of the lines of the input image.

4. The method according to claim 3, characterized in that it calculates the average by frequency for all lines, thus-thereby obtaining the average of the one-dimensional spectrum for all lines of the image.

5. The method according to claim 1, characterized in that Step 2 is a periodic square wave on the azimuthal axis, a square pulse with unit amplitude in the closed interval [−b, +b].

6. The method according to claim 5, characterized in that it defines a square wave as the infinite summation of spaced square pulses of a certain constant value.

7. The method according to claim 6, characterized in that it applies the discrete Fourier transform to the square wave.

8. The method according to claim 1, characterized in that, in Step 3, a filter is applied by multiplying the wave model and the two-dimensional Fourier transform of the input image.

9. The method according to claim 1, characterized in that, in Step 4, it obtains the final image processed by applying the inverse two-dimensional Fourier transform to the result obtained in Step 3.

10. The method according to claim 1, characterized in that it is applied to a tool that performs amplitude measurements by a transducer of emission and reception of ultrasonic waves.

11. The method according to claim 1, characterized in that it is applied to boreholes for any type of reservoir and without casing.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0020] The present invention will be described in more detail below, with reference to the attached figures which, in a schematic form and not limiting the inventive scope, represent examples of its embodiment. In the drawings, there are:

[0021] FIG. 1 illustrates an example of the method of automatic detection of pad repetition period, a necessary step for the automatic modeling of the filter and the characterization of the artifacts;

[0022] FIG. 2 illustrates a detail of the modeling of a square pulse with unity amplitude, between the interval [−b, +b];

[0023] FIG. 3 illustrates the process by which the pad artifacts that appear in the analyzed ultrasonic images can be modeled using a periodic square wave on the angular (horizontal) axis of the well wall, repeated vertically for a certain depth;

[0024] FIG. 4 illustrates an example of the frequency response of a filter generated using the technique proposed in this document, for a specific case of pad artifact composition. The magnitude is shown above and the phase below;

[0025] FIG. 5 illustrates four examples of the final result of the method, where it is possible to visualize the original image with the artifacts, next to the data processed using the method, with the artifacts removed and the geological information preserved.

DETAILED DESCRIPTION OF THE INVENTION

[0026] There follows below a detailed description of a preferred embodiment of the present invention, by way of example and in no way limiting. Nevertheless, it will be clear to a technician skilled on the subject, from reading this description, possible further embodiments of the present invention still comprised by the essential and optional features below.

[0027] The present invention provides a new method for overcoming the limitations of the state of the art for removing pad mark artifact noise in ultrasonic impedance imaging in reservoir wells. Furthermore, it also introduces a method of automatic characterization and modeling of this type of noise. There follows the detailed description of the method presented herein, by way of exemplification, referring to the figures.

[0028] Due to the fact that the method is based on image data that are usually obtained in well drilling maneuvers, it is not necessary to use additional resources to perform this task.

[0029] The method presented here uses acoustic impedance image data as input to characterize the noise of pad marks artifacts, by automatically modeling these artifacts as a periodic square wave on the horizontal axis. This modeling allows specific noise removal from the marks by filtering the model equation and the original image.

[0030] This method assumes carrying out a load of the data to be processed, by any method. Usually, this data is stored in files in DLIS or LAS format. It is not the purpose of this method to describe the method used to read/modify/write this type of files, but it is important to consider that this previous step of importing data will be necessary to use the method. The method requires the following steps:

Step 1

[0031] Automatic obtention of the period value of the pad marks. This step should only be applied to data with a complete azimuthal sweep, that is, data that has been obtained after sweeping 360° of the well wall and not just a portion of the same. For this purpose, it will be necessary to obtain the magnitude of the one-dimensional spectrum of the Fourier transform of each of the lines of the input image. Once the spectra are calculated, it will be necessary to calculate the average by frequency for all the lines, thus obtaining the average of the one-dimensional spectrum for all the lines of the image. The location of the peak in this distribution characterizes the frequency of pad marks. This process is exemplified in FIG. 1, and the calculation can be performed by applying equation (1):

[00001]$\begin{array}{cc}{\omega }_x=\underset{\forall x}{\mathrm{argmax}}.Math.\frac{1}{N_y}{.Math.}_{n=0}^{N_y-1}{\mathrm{fft}}_{1D}\left(I_n\right).Math.& \left(1\right)\end{array}$

where I represents the acoustic impedance image to be processed, n is the line index of the image (from 0 to the total size of lines minus 1), N.sub.y represents the total number of lines in the image, x represents the angular indices of the image (number of columns of the image) and ω.sub.x represents the repetition frequency of the pads.

Step 2

[0032] Modeling of pad artifacts as a periodic square wave on the azimuthal axis. A square pulse with unity amplitude in the closed interval [−b, +b], as shown in FIG. 2, can be expressed according to equation (2):

[00002]$\begin{array}{cc}{f}_x=\{\begin{array}{cc}1& -b\le x\le +b\\ 0& \mathrm{another}\mathrm{value}\mathrm{of}x\end{array}& \left(2\right)\end{array}$

where f represents the equation that defines the pulse amplitude.

[0033] FIG. 3 shows how it is possible to use the definition shown in (2) to define a square wave as the infinite summation of square pulses spaced a certain constant value, as expressed according to equation (3):

[00003]$\begin{array}{cc}{g}_{x,y}={.Math.}_{k=-\infty }^{\infty }f\left(x-\frac{2\mathrm{bk}}{D},y\right)& \left(3\right)\end{array}$

where g.sub.x,y represents the equation that defines the amplitude of the periodic square wave in x, f is given by equation (2), x represents the angular indices of the image (number of columns of the image), y represents the vertical indices (number of lines of the image), k is the index of each of the pulses that make up the signal (theoretically from −∞ to +∞), and D represents the proportion between the spacing between the pulses and the size of the pulse, also called “duty cycle” (a perfect square wave with equal spacing has a value of D=0.5).

[0034] By applying the discrete Fourier transform, it is possible to obtain the equivalent equation, in the frequency space, of (3), which is given by equation (4):

[00004]$\begin{array}{cc}{H}_{{\omega }_x,{\omega }_y}=\sqrt{\frac{2}{\pi }}\frac{\sin \left({\omega }_{x}D\frac{N}{N_p}\right)}{{\omega }_x}{.Math.}_{k=0}^{N_p-1}{e}^{-j{\omega }_x\frac{N}{N_p}\left(\frac{D}{2}+k\right)}\delta \left({\omega }_y\right)& \left(4\right)\end{array}$

where H.sub.ω.sub.x.sub.,ω.sub.y represents the equation defining the two-dimensional square wave given by (3) in the frequency space defined by ω.sub.x and ω.sub.y (horizontal and vertical frequency components, respectively), and N.sub.p is the number of pads, which can be calculated from the image size (number of columns) and the frequency of the pads obtained according to (1).

[0035] The value of D is a constant for the entire well that depends on the tool used. This information is accessible by means of the tool catalogs offered by their manufacturers and can be derived from the information on said tool obtained during the data import. The value can also be modified by the user when applying the method and, in case of indeterminacy, it can be considered D=0.5.

[0036] An example of the frequency response of the filter defined in equation (4) can be seen in FIG. 4. In this figure, the upper part shows the magnitude of this response, while the lower part shows the phase.

Step 3

[0037] Application of the filter by multiplying the wave model (4) and the two-dimensional Fourier transform of the input image.

Step 4

[0038] Obtaining the final image processed by applying the inverse two-dimensional Fourier transform to the result obtained in Step 3. FIG. 5 shows an actual example of the result obtained by processing by this method in ultrasonic impedance images of reservoir wells. In it, the original data are shown in 5.1, 5.3, 5.5, and 5.7, whereas the processed data, after applying the method shown here, appear in 5.2, 5.4, 5.6, and 5.8, respectively.

[0039] The invention disclosed herein is capable of being applied to a tool that performs amplitude measurements by a transducer of emission and reception of ultrasonic waves. It is also capable of being applied to boreholes for any type of reservoir and without casing.

[0040] Those skilled in the art will immediately appreciate the important benefits arising from the use of the present invention. This method does not require the use of any other well tool, so it does not incur any extra cost to the exploratory project. Furthermore, the ability to recover geological and lithological information masked by pad artifacts, due to the use of other well property mapping tools, increases not only the usability of the data in automatic methods (such as, for example, application of the data to characterization processes, classification, use of neural networks, texture detectors, etc.), but also improves the robustness and reliability of the results and the information that can be concluded and extracted therefrom by the specialists who analyze the same.