Method to correct eccentricity in ultrasonic image profiles

Abstract

This invention relates to a method of correcting eccentricity of ultrasonic image profiles, measured along a section of rock, including the steps of measuring emitted amplitude values (A.sub.0) and measuring transit time values (t.sub.) of ultrasonic acoustic pulses emitted at a range of default angles (), where the amplitude measured at each angle (A.sub.) is determined by an amplitude decay model in relation to the transit time (t.sub.) defined by A.sub.=A.sub.0e.sup.t.sup..sup./I.sub. where is the decay correction factor; and I.sub. is the reflection coefficient of the rock wall.

Claims

1. A method for correcting eccentricity of ultrasonic image profiles measured along a rock section wherein the rock section is included in a borehole of an uncoated reservoir, the method comprising: dropping a tool including an ultrasonic wave transducer in a center of an oil well; emitting from the tool an ultrasonic acoustic pulse at an angle , the ultrasonic acoustic pulse having an emitted amplitude value A.sub.0, wherein the ultrasonic pulse propagates through drilling fluid until reaching a wall of the oil well; receiving a reflected ultrasonic acoustic pulse reflected from a rock wall of the rock section; measuring a transit time value t.sub. corresponding to the reflected ultrasonic acoustic pulse; determining a measured amplitude A.sub. of the reflected ultrasonic acoustic pulse using an amplitude decay model in relation to the transit time t.sub.; obtaining a corrected amplitude A.sub.C, where the measured amplitude A.sub. is amplified by an inverse of the amplitude decay model, and re-scaled to an average value by an average (c) of corrections applied to a plurality of measured amplitudes corresponding to a plurality of angles included in the ultrasonic image profiles, wherein the amplitude decay model is defined by
A.sub.=A.sub.0e.sup.t.sup..sup./I.sub., wherein is a decay correction factor calculated by an optimization process, wherein I.sub. is a reflection coefficient of the rock wall, wherein the corrected amplitude A.sub.C is defined by: $A_c=A_{}.Math.e^{\frac{t_{}}{}}/\stackrel{\_}{c},$ wherein the average of the corrections applied (c) is defined by: $\stackrel{\_}{c}=\frac{1}{N}{.Math.}_{=1}^N\left(\frac{1}{q}{.Math.}_{\mathrm{prof}=1}^q{e}^{-t_{}\left(\mathrm{prof}\right)/}\right),$ and wherein N is a number of the plurality of angles, q is a number of measurement lines in a layer sampled by the ultrasonic acoustic pulse, and prof refers to a depth to which the measurement lines are measured.

2. The method of claim 1, further comprising calculating the decay correction factor () by the optimization process, wherein objective functions related to statistical standard deviation, skewness and kurtosis parameters are used.

3. The method of claim 2, wherein the objective functions comprise a multivariate objective function.

Description

BRIEF DESCRIPTION OF THE FIGURES

(1) The detailed description provided below is for the attached figures and their respective reference numbers, showing the modalities of this invention. The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.

(2) FIG. 1 illustrates a cross section of a well in which an ultrasonic image profile is to be generated, in which an ultrasonic pulse generating tool is positioned.

(3) FIG. 2a illustrates the first graphical representation of a well suffering a spiraling effect, a first cause of eccentricity of ultrasonic measurement tools.

(4) FIG. 2b illustrates a second graphical representation of a well suffering a spiraling effect, a first cause of eccentricity of ultrasonic measurement tools.

(5) FIG. 3 illustrates a second cause of eccentricity of ultrasonic measurement tools, which is related to the oval shape of a well.

(6) FIG. 4 illustrates a third cause of eccentricity of ultrasonic measurement tools, which is related to the eccentric positioning of the tool itself.

(7) FIG. 5 illustrates amplitude data of a homogeneous layer with an eccentricity effect (average amplitudes) and the same data after applying the method of this invention (average amplitudes corrected).

(8) FIGS. 6a-6c illustrate a comparison of a first uncorrected ultrasonic amplitude image profile after being corrected by a known method of the prior art, and after being corrected by the correction method of an ultrasonic images profile of this invention.

(9) FIGS. 7a-c illustrate a comparison of a second uncorrected ultrasonic amplitude image profile after being corrected after undergoing correction by a known method of the prior art, after being corrected by the method for correcting a profile ultrasonic images by this invention.

DETAILED DESCRIPTION OF THE INVENTION

(10) This invention involves a method to correct ultrasonic image profiles, preferably ultrasonic borehole image profiles for any type of uncoated reservoir, more preferably of oil wells. The correction in question is necessary due to eccentricity issues, due to several factors, among them well spiraling, the oval shape of the well and the eccentricity of the tool itself.

(11) Ultrasonic image profiles are acquired using a tool capable of emitting acoustic wave pulses that propagate through the drilling fluid until they reach the wall of the well, where they are reflected and bounce back to the tool that registers transit time and amplitude properties through the transducer.

(12) Sound waves are reflected due to the contrasting impedance between the fluid-rock medium. The higher the contrast in impedance, the greater the amplitude of the reflected waves. Stacking the amplitude and transit time at depth measurements generates a 360 image of the profiled section of the well.

(13) As discussed, the amplitude measured is greatly affected by the eccentricity since the acoustic waves suffer attenuation along their trajectory in the drilling fluid. The greater the distance traveled, the greater the attenuation, and vice versa. Thus varying distances between transducer and rock along the circumference of the well result in variable attenuation and, consequently, in variable amplitude along this circumference.

(14) These variations in amplitude represent a sum of the effect of geometric attenuation (eccentricity) and of the impedance contrast between the fluid and the reservoir. The transit time measurements provide a model of the variation of the distance between the transducer and the rock, and consequently the eccentricity. So, knowing that the control of transit time in the amplitude background measured in the profile is a function of an attenuation process, according to this invention, this phenomenon is modeled by a decay equation of amplitude to transit time.

(15) Therefore, in summary, the invention provides a method to correct eccentricity of ultrasonic image profiles, measured along a rock section, including the steps of measuring amplitude values emitted A.sub.0 and measuring transit time values t.sub. of the ultrasonic acoustic pulses emitted in a range of default angles, in which the amplitude measured at each angle A.sub. is determined by a decay model for amplitude to transit time t.sub., defined by the equation:
A.sub.=A.sub..Math.e.sup.t.sup..sup./.Math.I.sub.,(1)

(16) where is the decay correction factor; and I.sub. is the reflection coefficient of the rock wall.

(17) In other words, the amplitude A.sub. measured in each direction is a function of the amplitude emitted A.sub.0 multiplied by the exponential of the ratio of transit time t.sub. to decay correction factor and multiplied by the reflection coefficient of the rock wall I.sub..

(18) Preferably, ultrasonic acoustic pulses are emitted in a range of default angles, as described above. The default of such angles may be random so as to be dependent on software, for example. In other configurations, the default may be set at intervals of angles, such as at every 1 or 2, amongst other angle values. Optionally, the range of angles is 360, covering the whole circumference of a reservoir.

(19) After obtaining the decay factor of , the eccentricity of the amplitude image is corrected, which is amplified by the inverse of the decay equation and re-scaled to the average value by an average of corrections applied c in accordance with the equations (2) and (3) below.

(20) $\begin{array}{cc}{A}_c=A_{}.Math.e^{\left(\frac{t_{}}{}\right)}/\stackrel{\_}{c}& \left(2\right)\\ \stackrel{\_}{c}=\frac{1}{N}{.Math.}_{=1}^N\left(\frac{1}{q}{.Math.}_{\mathrm{prof}=1}^q{e}^{t_{}\left(\mathrm{prof}\right)/}\right)& \left(3\right)\end{array}$ Where: A.sub.c=Corrected Amplitude; c=Average of applied corrections; N=number of radial measurements (180); q=number of measurement lines in the layer sampled; and prof=each of the data lines referenced to a certain depth.

(21) Calculation of the decay correction factor in a complete image file, where the layers present their intrinsic heterogeneities is optionally achieved through an optimization process. To do so, the method of the invention uses objective functions with statistical parameters to determine the decay correction factor . It was found that the smaller the standard deviations, the higher the kurtosis deviations, and more negative values of asymmetry are consistent with the best correction results. Optionally, the objective function is a multivariate objective function.

(22) Note that the eccentricity correction method described applies to correcting the ultrasonic amplitude image profile of any type of tool.

(23) Although the method is described as being applicable to oil wells, it is important to point out that such an eccentricity correction method for ultrasonic imaging profiles can be applied to boreholes for any type of uncoated reservoir including, but not limited to, oil wells.

(24) Preferably, the eccentricity correction method can be applied to a tool that measures amplitude and transit time by a transducer emitting and receiving ultrasonic waves.

(25) FIG. 5 illustrates amplitude values extracted from data from a layer of homogeneous rock with acoustic properties. The Average Amplitudes curve is the oscillatory behavior of the amplitude controlled by varying the distance between the transducer and the well wall. The Corrected Average Amplitudes curve was generated after applying the method of the invention, where one can observe removal of this oscillation and conservation of small variations intrinsic to the rock.

(26) For additional tests, the correction method for ultrasonic image profiles of this invention was applied to the acoustic image profiles of two wells with different degrees of eccentricity, and compared with the known correction method of the prior art. The results are illustrated in FIGS. 6 and 7.

(27) As will be clear from the explanations below, note that the corrections obtained were very satisfactory and even allowed for the extraction of permo-porous properties of the reservoirs.

(28) FIG. 6 illustrates the comparison of a first uncorrected ultrasonic image profile, with a known correction method of the prior art and with the proposed correction method, where the ultrasonic amplitude image profile of a first reservoir without any correction method applied is initially illustrated 6a. Then the same profile is illustrated after undergoing a correction process known from art 6b. Finally, the same profile 6a is shown after undergoing the 6c correction process shown here. By comparing the images shown, it is observed that after application of the eccentricity correction method of this invention 6c the eccentricity effect is substantially reduced, while the amplitude properties of the reservoir layers are preserved. Which is not obtained by the known correction method of the prior art.

(29) FIG. 7 illustrates the ultrasonic amplitude image profile of a reservoir that has both cave and large vug geological structures, in which the first image 7a illustrates the profile without any correction. Then the same profile 7a is shown after undergoing a known correction process of the prior art 7b. Finally, the same profile 7a is shown after undergoing the correction process that is revealed here 7c. Note that after applying the eccentricity correction method of this invention 7c, said geological structures (caves and large vugs) keep their amplitude characteristics significantly lower than the matrix, which allows for distinguishing and segmenting them for quantitative analysis. The image corrected by the known method of the prior art 7b, in turn, offers general homogenization of the image, which distorts the results of the correction, hindering and even making it impossible to distinguish the illustrated vugs and caves.

(30) Thus, in summary, this invention provides an eccentricity correction method of ultrasonic amplitude image profiles after acquisition. The method uses the transit time measurements at each angle as a model of the eccentricity of the tool. Variations in amplitudes generated by the geometric effect of the eccentricity are removed by amplifying the measured amplitude using the reverse decay equation and the amplitude image is rescaled by the average of corrections applied.

(31) $A_c=A_{}.Math.e^{\left(\frac{t_{}}{}\right)}/\stackrel{\_}{c}.$

(32) It is therefore clear that the eccentricity correction method here revealed corrects ultrasonic image profiles, greatly reducing the eccentricity effect, which makes qualitative and quantitative interpretation of geological data from the ultrasound images reliable.

(33) Thus, for the oil and gas industry, the correction method of ultrasonic image profiles shown here is a huge breakthrough for the study and characterization of oil reservoirs, in which ultrasound image profile data is quantitatively analyzed to characterize the reservoir flow capacity at high resolution at the scale of the well, especially in reservoirs with geological structures such as caves, extended fractures and large vugs, since it is widely recognized that the presence of these structures in oil wells represents a significant risk of circulation losses, loss of wells and may even cause problems in production and injection tests.

(34) The quantitative use of these profiles also allows characterization of these structures so that they are incorporated into formation test planning and into the geological model of the reservoirs.