US20210285324A1 - Method for Determining the Integrity of a Solid Bonding between a Wellbore and a Casing - Google Patents
Method for Determining the Integrity of a Solid Bonding between a Wellbore and a Casing Download PDFInfo
- Publication number
- US20210285324A1 US20210285324A1 US17/159,341 US202117159341A US2021285324A1 US 20210285324 A1 US20210285324 A1 US 20210285324A1 US 202117159341 A US202117159341 A US 202117159341A US 2021285324 A1 US2021285324 A1 US 2021285324A1
- Authority
- US
- United States
- Prior art keywords
- envelope
- attenuation
- integral
- domain
- receiver
- 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.)
- Granted
Links
- 239000007787 solid Substances 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 title claims description 51
- 238000005192 partition Methods 0.000 claims abstract description 38
- 239000000463 material Substances 0.000 claims abstract description 28
- 230000001902 propagating effect Effects 0.000 claims abstract description 9
- 230000000644 propagated effect Effects 0.000 claims abstract description 6
- 238000000926 separation method Methods 0.000 claims abstract description 3
- 238000012545 processing Methods 0.000 claims description 11
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 239000007788 liquid Substances 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 8
- 238000007689 inspection Methods 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 230000003213 activating effect Effects 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims description 2
- 239000004568 cement Substances 0.000 description 18
- 235000019687 Lamb Nutrition 0.000 description 12
- 239000006185 dispersion Substances 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 239000012530 fluid Substances 0.000 description 7
- 230000008878 coupling Effects 0.000 description 5
- 238000010168 coupling process Methods 0.000 description 5
- 238000005859 coupling reaction Methods 0.000 description 5
- 229910000831 Steel Inorganic materials 0.000 description 4
- 238000006073 displacement reaction Methods 0.000 description 4
- 239000003921 oil Substances 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 239000011343 solid material Substances 0.000 description 4
- 239000010959 steel Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 239000003129 oil well Substances 0.000 description 3
- 239000010755 BS 2869 Class G Substances 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 2
- TZCXTZWJZNENPQ-UHFFFAOYSA-L barium sulfate Chemical compound [Ba+2].[O-]S([O-])(=O)=O TZCXTZWJZNENPQ-UHFFFAOYSA-L 0.000 description 2
- 229910052601 baryte Inorganic materials 0.000 description 2
- 239000010428 baryte Substances 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 239000006260 foam Substances 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000011796 hollow space material Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/005—Testing the nature of borehole walls or the formation by using drilling mud or cutting data
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/005—Monitoring or checking of cementation quality or level
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/12—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
- E21B47/14—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling using acoustic waves
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
- E21B49/003—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells by analysing drilling variables or conditions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
- G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson
Definitions
- the present invention relates to evaluating a material on one side of a partition using ultrasonic measurements.
- Deep surface wells are created by drilling a hole into the earth with a drilling rig that rotates a drill string with a bit attached. After the hole is drilled to a prescribed depth, sections of steel tubing known as casing are set in the hole. To be able to lower it down, the casing is slightly smaller than the borehole. The small space between the formation and the casing is cemented in order to prevent oil and gas migrating up to the surface. This is done by pushing the cement down to the end of casing and squeezing it from the prescribed depth and up the annulus. During the drilling, completion and permanent closure of an oil well, it is a challenge to verify if the cement has been squeezed up the annulus and bonded to outside of the casing surface.
- the primary method to validate that the annulus is impermeable is to use ultrasonic borehole logging.
- ultrasonic borehole logging To improve the identification of whether an impermeable solid is bonding onto the casing, efforts have been made to provide solutions for using ultrasonic logging in the field.
- US2006198243 discloses a method and apparatus for determining the integrity of a cement bond log disposed in the annular space between a casing and a wellbore.
- the method and apparatus induce a Lamb wave in the casing and into the wellbore.
- the Lamb wave attenuates upon passage beside the cement bond.
- the integrity of the cement bond log can be determined by analysis and evaluation of the attenuation results.
- U.S. Pat. No. 9,534,487 uses a transmitter and several receivers to measure the zero-order mode (S 0 ) symmetric extensional Lamb waveforms. These are processed to obtain the amplitude attenuation dispersion plot and phase dispersion plot as functions of frequency as means to characterize the physical state of a casing and annular fill outside the casing.
- the technique obtains information by transmitting a pressure wave (T 1 ) with oblique incidence at the casing wall that excites a guided wave propagating upwards in the casing.
- the guided wave leaks off energy as it propagates, and the leaked off pressure wavefront is measured at receivers positioned further up the well. Since Lamb waves leak energy as they propagate, the waves are constantly attenuated as the distance increases.
- U.S. Pat. No. 9,534,487 also estimates the amplitude attenuation dispersion and phase dispersion of the extensional Lamb wave, where the ultrasonic waveforms are processed to obtain (i) an amplitude attenuation dispersion plot of attenuation as a function of frequency and (ii) a phase dispersion plot of phase velocity as a function of frequency.
- Barrier wavespeeds are determined by identifying discontinuities within the amplitude attenuation dispersion plot, which are then related to the discontinuities to barrier wavespeeds using the phase dispersion plot.
- U.S. Pat. No. 10,138,727 processes the zero-order mode antisymmetric Lamb flexural waveforms to identify barrier parameters as a function of azimuth and depth along the borehole, wherein the waveforms comprise at least two of sonic signals, ultrasonic pulse-echo signals, and ultrasonic pitch-catch signals.
- the patent states that the preprocessing extract attributes particular to each measurement such as decay rate for the pulse echo and attenuation for the flexural and extensional modes. These attributes are then compared to a multi-dimensional library of pre-calculated attributes as a function of cement properties and bond properties. A match up to a specific criterion provides the set of inverted cement parameters.
- the attenuation can be estimated across two or more receivers. Two attributes can be used: the energy envelope peak decay as well as the attenuation dispersion (i.e., as a function of frequency).
- Embodiments involve receiving acoustic cement evaluation data from one or more acoustic downhole tools used over a depth interval in a well having a casing. Embodiments also involve determining an extensional mode component and a flexural mode component of the leaky-Lamb wave measurements and determining a flexural wave attenuation based on the flexural mode component of the leaky-Lamb wave measurement.
- US2015218930 discloses systems and methods for estimating an acoustic property of an annulus in a cement evaluation system.
- the present invention differs from known approaches in integrating the envelope of the S 0 mode within a predefined time window or interval in the time domain to aggregate the amount of energy, calculating the attenuation of S 0 mode from one receiver to the next, computing the attenuation parameter against the acoustic impedance, and converting the results to a probability estimate of the likeliness of a solid behind the partition.
- a further feature of the invention is that processing targets the early arrival of the S 0 mode where the phase and group velocity is considered more constant. In addition, the early arrival is less influenced by other Lamb modes or direct waves due to the high velocity of the S 0 mode for the lower frequencies.
- an embodiment of the present invention estimates the probability of a solid material being present, i.e. a method for differentiating a solid material from a fluid.
- the invention quantifies the integral attenuation over a predefined interval as will be discussed below.
- FIG. 1 shows schematically a measuring method and system for characterizing the physical state of a partition installed in a borehole
- FIGS. 2 a and 2 b shows schematically the orientation of the transducer with respect to the partition without and with shear coupling respectively;
- FIGS. 3 a and 3 b shows the phase and group velocity of Lamb modes excited in a steel plate in vacuum
- FIG. 4 is a graphical representation showing numerical waveforms similar as expected from extensional logging in the field
- FIGS. 5 a and 5 b are graphical representations showing the envelope of the extensional pulse when the annulus is filled with a liquid or solid, respectively;
- FIGS. 6 a and 6 b are graphical representations showing the attenuation of the extensional mode according to the invention.
- FIG. 7 is a graphical representation showing the attenuation vs the acoustic impedance for different materials in the annulus estimated according the invention.
- FIG. 8 is a generic plot showing the principal analysis method according to the invention.
- FIG. 9 is a generic flowchart showing the principal analysis method according to the invention.
- FIG. 1 shows schematically a measuring method and system according to an embodiment of the invention for characterizing the physical state of a partition installed in a borehole.
- the partition is constituted by the wall of a pipe known as a “casing” that is installed in a borehole, thus separating a first domain inside the casing from a second domain outside the casing.
- the casing is surrounded by solid geological formation, which in the case of an oil or gas well is typically the seabed.
- the casing usually comprises multiple sections of progressively narrower the deeper they are inserted into the geological formation.
- each section of casing After installation of each section of casing, it is cemented inside the borehole so as to form a secure rigid enclosure around the casing, which serves as a conduit for accommodating the drill and releasing the oil or gas through perforations made in the portion of the casing which passes through the production zone, to provide a path for the oil to flow from the surrounding rock into the production tubing.
- a secure rigid enclosure around the casing which serves as a conduit for accommodating the drill and releasing the oil or gas through perforations made in the portion of the casing which passes through the production zone, to provide a path for the oil to flow from the surrounding rock into the production tubing.
- At least one transmitter T 1 is positioned inside the casing, where the transmitted signal propagates and hits the wall of the casing at an oblique incidence.
- the transmitted signal excites a guided wave inside the wall of the casing formed of a known material and having a known thickness.
- a transducer positioned at an angle inside a liquid-filled pipe can excite guided waves that are similar to Lamb waves.
- the diameter, d c , and the thickness, t c , of the cylindrical pipe determine the deviation between the guided waves and the Lamb waves. If the ratio is above 10, which is typically the case in the field, the difference is negligible.
- the similarity also depends on the frequency, but if the wavelength of the guided wave is much less than the pipe circumference and is greater than 10, the effect of the curvature becomes insignificant.
- the well geometry can be understood in 2-dimensional spatial coordinates, where the length of the partition is the axial direction and the azimuthal direction is simplified as only one direction. Further, in the frequency regime of interest, the two zero-order modes dominate the propagating wave in the partition.
- the zero-order anti-symmetric mode (A 0 ) referred to as the flexural wave has an elliptical particle-motion mainly perpendicular to the partition, i.e. in the axial direction. This creates a ‘bending’ or ‘flexural’ motion.
- the particle motion in the casing is elliptical, with the vertex pointing in the direction of the surrounding material, resulting in waves being emitted from each side of the casing.
- the zero-order symmetric mode (S 0 ) referred to as the extensional wave has an elliptical particle displacement that is mainly parallel to the partition, i.e. in the axial direction, which results in the vertex pointing in the direction of the propagating wave. Consequently, the largest proportion of energy will be transmitted as the tangential displacement vector. Therefore, if the material behind the partition consists of a material without shear modulus, only a small fraction of energy will leak off from the steel plate constituting Domain 3 into Domain 2. If there is a shear coupling at Interface 2, sufficiently more energy will radiate into Domain 2, i.e. if a solid material is bonding to the casing.
- FIG. 2 The excitation of S 0 at the pipe/plate requires an oblique incidence angle around 15° if Domain 1 is filled with water and Domain 3 is made of steel. This is illustrated in FIG. 2 , where FIG. 2 a relates to a fluid in Domain 2 and FIG. 2 b relates to a cement in Domain 2.
- FIG. 3 show the phase and group velocity of some of the low Lamb modes, including A 0 and S 0 .
- the horizontal axis is the frequency (f) multiplied by the casing thickness (d), which makes it convenient to compare different dimensions if the ratio f ⁇ d remains constant.
- the figure demonstrates that S 0 is highly dispersive, however, at the earlier frequencies the group velocity is considered more constant.
- the system also requires at least two receivers, as seen in FIG. 1 , as R 1 and R 2 , both of which are located in Domain 1 and oriented in mirror symmetry to the transducer relative to a line normal to the outer surface of the casing at an angle of reflection that is equal to the incidence angle of the transducer.
- the extensional wave leaks off energy as it propagates, and the leaked off pressure wavefront is measured at receivers positioned further up the well. Since the extensional wave leaks energy as it propagates, the wave is being constantly attenuated as the distance increases.
- the materials on both sides of the casing i.e. in Domains 1 and 2 affect the attenuation. Assuming the material inside Domain 1 is known, measuring the attenuation from one receiver to the next provides information about the acoustic properties of the material constituting the annular space (Domain 2).
- FIG. 4 shows numerical waveforms similar as expected from extensional logging in the field and illustrates the dispersion smearing out the pulse, where the first waveform is from R 1 and the second waveform is from R 2 , as shown in FIG. 1 .
- the first pulse is the arrival of the emitted wavefront from the extensional wave (S 0 ), but as seen in FIG. 3 , the group velocity is highly dispersive, which results in the pulse being smeared out.
- the brace above the pulse in FIG. 4 indicates the range of interest of the S 0 mode, and within this range it can be seen that the pulse is preserved due to approximately constant group velocity in that frequency range (as seen in FIG. 3 at the low frequencies).
- the rest of the waveform is a combination of the S 0 mode, the A 0 mode, and the fluid borne wave (direct wave), and this region is not qualified for evaluating the attenuation of the envelope.
- the region where the S 0 mode is being evaluated has only traveled in domains 1 and 3, but the attenuation within this region is highly affected by whether the material in Domain 2 consists of a solid.
- the extensional wave propagates along the casing, waves are constantly being leaked off. However, in the absence of a shear coupling at interface 2 , the attenuation is significantly reduced.
- a shear modulus in Domain 2 allows coupling of the S 0 mode and the waves leaking off in both the normal direction and the tangential direction, as illustrated in FIG. 2 b by the arrow in the cement layer.
- FIG. 5 shows the envelope of the waveforms and the time axis is minimized to demonstrate the attenuation of the early S 0 mode.
- the waveforms in FIG. 5 a are the same as in FIG. 4
- the envelopes in FIG. 5 b are also from a numerical simulation using cement acoustic properties.
- the attenuation can be found by measuring the peak amplitude at two different distances and divide the peak measured at the further distance by the peak measured at the nearest distance, as illustrated by the receivers in FIG. 1 .
- the S 0 mode presented in FIG. 3 shows that the group velocity is constant only within a small frequency range, so a better technique is to calculate the integral within a predefined time window, as seen in the Eq. (1) and estimate the attenuation.
- the integral attenuation ( ⁇ dB ) in decibel (dB) is then found as given below:
- ⁇ d ⁇ B 20 ⁇ log 1 ⁇ 0 ⁇ ( E a ⁇ m ⁇ p , 2 E a ⁇ m ⁇ p , 1 ) ( 2 )
- equation (2) can be then used to estimate the loss in energy within a predefined time-window between the two waveforms.
- FIG. 6 shows integral attenuation, estimated by using Eq. (2), plotted as a logarithmic function, where FIG. 6 a is the attenuation from the envelopes seen in FIG. 5 a , and FIG. 6 b is the attenuation from the envelopes seen in FIG. 5 b .
- the horizontal axis represents the distance the S 0 mode has propagated inside Domain 3. From the figure one can see that the attenuation is significantly higher (more than twice) when Domain 2 consists of a solid.
- the attenuation presented in FIG. 6 can be quantified by dividing it by the distance propagated to estimate attenuation per centimeter (cm):
- ⁇ z,dB is the attenuation expressed in dB/cm, which is plotted in FIG. 7 .
- the function in FIG. 7 shows ⁇ z,dB plotted against the acoustic impedance for various different solids and fluids in Domain 2. The data seen in the figure were obtained from a study using numerical simulations, where the acoustical properties of the materials are presented in Table 1 below.
- all attenuation measurements within the inspection run are used to create a model, and the minimum and maximum attenuation are used to define the outer boundaries.
- the threshold is set between the outer boundaries, and the value may vary from each inspection run.
- the invention provides a complementary method of evaluating ⁇ z,dB .
- ⁇ z,dB instead of plotting ⁇ z,dB against the acoustic impedance, one can plot it against the shear wave velocity. This is seen in FIG. 8 , where ⁇ z,dB increases with increasing shear wave velocity. From this, one can make a probability estimate of the likeliness of a solid in Domain 2. A higher attenuation then equals a higher probability of a solid.
- FIG. 9 shows a generic flowchart to describe the processing to estimate ⁇ z,dB between the extensional waves excited in the casing.
- the flowchart starts with the input of the extensional wave and knowledge of the distance between the receivers.
- the final output is parameter to quantify the peak envelope attenuation.
- FIG. 9 shows how to pick the beginning and end time of the time window, t 1 and t 2 , based on the amplitude within the envelope of the S 0 mode.
- the start time, t 1 is defined as where the signal is first above the signal noise level, which in some embodiments of the invention is defined as ⁇ 30 dB of the maximum amplitude within the normalized envelope of the whole waveform, although other noise thresholds may be employed. Normalizing the envelope means that the maximum amplitude is always at 0 dB.
- the time window needs to end before the arrival of A 0 mode and/or the direct wave, so it does not influence the integral.
- t 2 is calculated from the second peak in the envelope as seen in FIG. 5 , where the method estimates where the amplitude drops below the half maximum of the second peak.
- the invention has been described primarily with regard to a method because the hardware is known per se.
- the transmitter and receivers When used to monitor oil wells, the transmitter and receivers must be inserted into a hollow casing sunk in the ground that is possibly hundreds of meters deep. Since the distances between adjacent receivers must be known, the transducers are installed into an elongated unit that is then lowered into the casing. This also allows the transducers to be tilted at known angles of incidence and reflection although more advanced transducers are available where the wave angle can be adjusted electronically. However, when applied to other situations where the transducers are more accessible, they may be spaced apart discretely at known intervals.
- processing may be performed by a suitably programmed computer.
- the invention contemplates a computer program being readable by a computer for executing the method of the invention.
- the invention further contemplates a machine-readable memory tangibly embodying a program of instructions executable by the machine for performing the method of the invention.
- inventive concepts encompassed by the invention include the following:
- Inventive concept 1 A method for evaluating a material on a remote side of a partition separating a first domain from a second domain, the method comprising:
- the second domain contains a solid if the measured integral envelope attenuation per distance is higher than a predetermined threshold.
- inventive concept 2 The method according to inventive concept 1, wherein establishing that the second domain contains a solid comprises:
- Inventive concept 3 The method of inventive concept 2, wherein computing the time interval includes (i) defining a start time t 1 as where the signal is first above a signal noise level; (ii) obtaining a second peak amplitude within the envelope (iii) defining an end time t 2 as where the signal drops below half the second peak amplitude; (iv) and using the interval between t 1 and t 2 as boundaries for the integration of the envelope.
- Inventive concept 4 The method of inventive concept 3, wherein the integral of the envelope is calculated according to:
- Inventive concept 5 The method of inventive concept 4, wherein computing the integral attenuation includes dividing the integral estimated within predefined time window of the envelope of the measured at the receiver furthest by the integral estimated within predefined time window of the envelope of the measured at the nearest receiver.
- Inventive concept 6 The method according to inventive concept 5, wherein the integral envelope attenuation per distance ( ⁇ z,dB ) is calculated according to:
- ⁇ z , dB 20 ⁇ log 1 ⁇ 0 ⁇ ( E amp , n E amp , 1 ) z 2 - z 1
- (E amp,n ) is the integral of the envelope of the extensional pulse at the n th receiver
- (P m,1 ) is the integral of the envelope of the extensional pulse at the first receiver
- Inventive concept 7 The method according to inventive concept 3, wherein the integral of the envelope amplitude is estimated from a portion of the time domain wherein the symmetric zero-order wavefront (S 0 ) arrives before a higher energy antisymmetric wavefront (A 0 ) is able to reach said portion thereby eliminating the effect of interference between the symmetric and the antisymmetric wavefronts.
- Inventive concept 8 The method according to any one of inventive concepts 3 to 7, wherein the predetermined noise threshold is ⁇ 30 dB of the peak amplitude within the normalized envelope of the whole waveform.
- Inventive concept 9 The method according to any one of inventive concepts 2 to 8, further including computing the integral envelope attenuation per distance ( ⁇ z,dB ) against a probability chart to determine the likelihood of the second domain containing a solid.
- Inventive concept 10 The method according to claim 9 , wherein if the attenuation is greater than a predetermined threshold, this is indicative of the material being a solid.
- Inventive concept 11 The method according to inventive concept 1, wherein establishing that the second domain contains a solid comprises:
- Inventive concept 12 The method according to inventive concept 11, wherein the predetermined attenuation threshold is in the order of 7 dB/cm.
- Inventive concept 13 The method according to any one of the preceding inventive concepts, wherein the partition comprises a casing and an annular fill installed in a borehole traversing a formation.
- Inventive concept 14 The method according to any one of the preceding inventive concepts, wherein each transmitter and receiver is directed toward the partition at an angle of incidence close to 15°.
- Inventive concept 15 A computer program product comprising a computer readable medium storing program code instructions, which when executed on at least one processor that receives as input data representative of the recorded ultrasonic waveforms from a pair of spaced apart receivers, carries out the processing and computing operations of any one of the preceding inventive concepts.
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Quality & Reliability (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
Material on a remote side of a partition separating first and second domains is evaluated by disposing at least one ultrasonic transmitter and a plurality of spaced ultrasonic receivers along a first side of the partition in the first domain. The transmitter is activated to form ultrasonic waveforms that comprise propagated quasi leaky-Lamb waves constituting extensional waves having symmetrical zero-order modes within the partition. The received ultrasonic waveforms are processed to determine the envelope of the waves propagating through the partition from a first receiver to a second receiver that is located more remote from the transmitter than the first receiver and whose separation from the first receiver is known. The integral envelope attenuation per distance of the wave propagating through the second domain over a predefined interval establishes that the second domain contains a solid when higher than a predetermined threshold.
Description
- This application claims to priority to GB2003621.6, which was filed on Mar. 12, 2020, and is incorporated herein by reference in its entirety.
- The present invention relates to evaluating a material on one side of a partition using ultrasonic measurements.
- Deep surface wells are created by drilling a hole into the earth with a drilling rig that rotates a drill string with a bit attached. After the hole is drilled to a prescribed depth, sections of steel tubing known as casing are set in the hole. To be able to lower it down, the casing is slightly smaller than the borehole. The small space between the formation and the casing is cemented in order to prevent oil and gas migrating up to the surface. This is done by pushing the cement down to the end of casing and squeezing it from the prescribed depth and up the annulus. During the drilling, completion and permanent closure of an oil well, it is a challenge to verify if the cement has been squeezed up the annulus and bonded to outside of the casing surface. If no pressure test is planned, the primary method to validate that the annulus is impermeable, is to use ultrasonic borehole logging. To improve the identification of whether an impermeable solid is bonding onto the casing, efforts have been made to provide solutions for using ultrasonic logging in the field.
- US2006198243 discloses a method and apparatus for determining the integrity of a cement bond log disposed in the annular space between a casing and a wellbore. The method and apparatus induce a Lamb wave in the casing and into the wellbore. The Lamb wave attenuates upon passage beside the cement bond. The integrity of the cement bond log can be determined by analysis and evaluation of the attenuation results.
- U.S. Pat. No. 9,534,487 uses a transmitter and several receivers to measure the zero-order mode (S0) symmetric extensional Lamb waveforms. These are processed to obtain the amplitude attenuation dispersion plot and phase dispersion plot as functions of frequency as means to characterize the physical state of a casing and annular fill outside the casing. The technique obtains information by transmitting a pressure wave (T1) with oblique incidence at the casing wall that excites a guided wave propagating upwards in the casing. The guided wave leaks off energy as it propagates, and the leaked off pressure wavefront is measured at receivers positioned further up the well. Since Lamb waves leak energy as they propagate, the waves are constantly attenuated as the distance increases.
- U.S. Pat. No. 9,534,487 also estimates the amplitude attenuation dispersion and phase dispersion of the extensional Lamb wave, where the ultrasonic waveforms are processed to obtain (i) an amplitude attenuation dispersion plot of attenuation as a function of frequency and (ii) a phase dispersion plot of phase velocity as a function of frequency. Barrier wavespeeds are determined by identifying discontinuities within the amplitude attenuation dispersion plot, which are then related to the discontinuities to barrier wavespeeds using the phase dispersion plot.
- U.S. Pat. No. 10,138,727 processes the zero-order mode antisymmetric Lamb flexural waveforms to identify barrier parameters as a function of azimuth and depth along the borehole, wherein the waveforms comprise at least two of sonic signals, ultrasonic pulse-echo signals, and ultrasonic pitch-catch signals. The patent states that the preprocessing extract attributes particular to each measurement such as decay rate for the pulse echo and attenuation for the flexural and extensional modes. These attributes are then compared to a multi-dimensional library of pre-calculated attributes as a function of cement properties and bond properties. A match up to a specific criterion provides the set of inverted cement parameters. The attenuation can be estimated across two or more receivers. Two attributes can be used: the energy envelope peak decay as well as the attenuation dispersion (i.e., as a function of frequency).
- US2016209539 discloses methods and systems for separating leaky-Lamb modes to evaluate cylindrically layered structures. Embodiments involve receiving acoustic cement evaluation data from one or more acoustic downhole tools used over a depth interval in a well having a casing. Embodiments also involve determining an extensional mode component and a flexural mode component of the leaky-Lamb wave measurements and determining a flexural wave attenuation based on the flexural mode component of the leaky-Lamb wave measurement.
- US2015218930 discloses systems and methods for estimating an acoustic property of an annulus in a cement evaluation system.
- Reference is also made to the following publications:
- Xiao He et al. “Ultrasonic leaky flexural waves in multilayered media: Cement bond detection for cased wellbores”, March 2014, Geophysics 79(2):A7-A11. This publication investigates ultrasonic leaky flexural modes in radially multilayered media through finite-difference modeling to evaluate cementing for the whole annulus of a casing.
- Erlend Magnus Viggen et al., “Simulation and modeling of ultrasonic pitch-catch through-tubing logging,” July 2016, Geophysics 81: D383-D393. This publication discloses the use of a finite-element model of a double-casing geometry with a two-receiver pitch-catch setup for logging in double-casing geometries.
- In accordance with the invention there are provided methods for evaluating a material on a remote side of a partition having the features of the independent claims.
- The present invention differs from known approaches in integrating the envelope of the S0 mode within a predefined time window or interval in the time domain to aggregate the amount of energy, calculating the attenuation of S0 mode from one receiver to the next, computing the attenuation parameter against the acoustic impedance, and converting the results to a probability estimate of the likeliness of a solid behind the partition. A further feature of the invention is that processing targets the early arrival of the S0 mode where the phase and group velocity is considered more constant. In addition, the early arrival is less influenced by other Lamb modes or direct waves due to the high velocity of the S0 mode for the lower frequencies. Since a more constant group velocity results in the mode being less smeared out, calculation of the integral within the predefined time interval and the computation of the attenuation, results in the method being more robust than just evaluating a peak. Since the S0 mode is highly affected by the presence of a shear modulus in the material behind the partition, computing the integral attenuation makes it possible to obtain information about the shear coupling. Thus, an embodiment of the present invention estimates the probability of a solid material being present, i.e. a method for differentiating a solid material from a fluid. The invention quantifies the integral attenuation over a predefined interval as will be discussed below.
- In order to understand the invention and to see how it may be carried out in practice, embodiments will now be described, by way of non-limiting example only, with reference to the accompanying drawings, in which:
-
FIG. 1 shows schematically a measuring method and system for characterizing the physical state of a partition installed in a borehole; -
FIGS. 2a and 2b shows schematically the orientation of the transducer with respect to the partition without and with shear coupling respectively; -
FIGS. 3a and 3b shows the phase and group velocity of Lamb modes excited in a steel plate in vacuum; -
FIG. 4 is a graphical representation showing numerical waveforms similar as expected from extensional logging in the field; -
FIGS. 5a and 5b are graphical representations showing the envelope of the extensional pulse when the annulus is filled with a liquid or solid, respectively; -
FIGS. 6a and 6b are graphical representations showing the attenuation of the extensional mode according to the invention; -
FIG. 7 is a graphical representation showing the attenuation vs the acoustic impedance for different materials in the annulus estimated according the invention; -
FIG. 8 is a generic plot showing the principal analysis method according to the invention; and -
FIG. 9 is a generic flowchart showing the principal analysis method according to the invention. -
FIG. 1 shows schematically a measuring method and system according to an embodiment of the invention for characterizing the physical state of a partition installed in a borehole. The partition is constituted by the wall of a pipe known as a “casing” that is installed in a borehole, thus separating a first domain inside the casing from a second domain outside the casing. The casing is surrounded by solid geological formation, which in the case of an oil or gas well is typically the seabed. The casing usually comprises multiple sections of progressively narrower the deeper they are inserted into the geological formation. After installation of each section of casing, it is cemented inside the borehole so as to form a secure rigid enclosure around the casing, which serves as a conduit for accommodating the drill and releasing the oil or gas through perforations made in the portion of the casing which passes through the production zone, to provide a path for the oil to flow from the surrounding rock into the production tubing. Thus, with reference toFIG. 1 it is possible to define four domains as follows: -
-
Domain 1 is the hollow space inside the casing; -
Domain 2 is the annular space between the outer wall of the casing and the formation. If there is no annular space, i.e. the formation is bonded onto the casing,Domain 2 will vanish; and -
Domain 3 is the wall, which constitutes a partition between the formation and the first domain. In an oil well the partition refers to a casing which is a larger pipe that is assembled and inserted into the recently drilled section of the borehole; -
Domain 4 is the geological formation.
-
- We will refer to these four domains as first, second, third and fourth domains, respectively. Further, “
Interface 1” is defined betweenDomains Interface 2” is betweenDomains Interface 3” betweenDomains FIG. 1 . - At least one transmitter T1 is positioned inside the casing, where the transmitted signal propagates and hits the wall of the casing at an oblique incidence. The transmitted signal excites a guided wave inside the wall of the casing formed of a known material and having a known thickness.
- From theory, a transducer positioned at an angle inside a liquid-filled pipe can excite guided waves that are similar to Lamb waves. The diameter, dc, and the thickness, tc, of the cylindrical pipe determine the deviation between the guided waves and the Lamb waves. If the ratio is above 10, which is typically the case in the field, the difference is negligible. The similarity also depends on the frequency, but if the wavelength of the guided wave is much less than the pipe circumference and is greater than 10, the effect of the curvature becomes insignificant.
- The well geometry can be understood in 2-dimensional spatial coordinates, where the length of the partition is the axial direction and the azimuthal direction is simplified as only one direction. Further, in the frequency regime of interest, the two zero-order modes dominate the propagating wave in the partition. The zero-order anti-symmetric mode (A0) referred to as the flexural wave, has an elliptical particle-motion mainly perpendicular to the partition, i.e. in the axial direction. This creates a ‘bending’ or ‘flexural’ motion. Thus, the particle motion in the casing is elliptical, with the vertex pointing in the direction of the surrounding material, resulting in waves being emitted from each side of the casing. Since the vortex is pointing out in the normal direction of the partition, the largest proportion of energy will be transmitted as the normal displacement vector. This makes the A0 mode less dependent of the material in
Domain 2, since both fluids and solids can couple waves radiating as the normal displacement vector from the partition. - The zero-order symmetric mode (S0) referred to as the extensional wave, has an elliptical particle displacement that is mainly parallel to the partition, i.e. in the axial direction, which results in the vertex pointing in the direction of the propagating wave. Consequently, the largest proportion of energy will be transmitted as the tangential displacement vector. Therefore, if the material behind the partition consists of a material without shear modulus, only a small fraction of energy will leak off from the steel
plate constituting Domain 3 intoDomain 2. If there is a shear coupling atInterface 2, sufficiently more energy will radiate intoDomain 2, i.e. if a solid material is bonding to the casing. - The excitation of S0 at the pipe/plate requires an oblique incidence angle around 15° if
Domain 1 is filled with water andDomain 3 is made of steel. This is illustrated inFIG. 2 , whereFIG. 2a relates to a fluid inDomain 2 andFIG. 2b relates to a cement inDomain 2. - Countless different Lamb modes can be excited inside a partition depending on the transducer frequency, partition thickness and the incidence angle.
FIG. 3 show the phase and group velocity of some of the low Lamb modes, including A0 and S0. The horizontal axis is the frequency (f) multiplied by the casing thickness (d), which makes it convenient to compare different dimensions if the ratio f·d remains constant. The figure demonstrates that S0 is highly dispersive, however, at the earlier frequencies the group velocity is considered more constant. - The system also requires at least two receivers, as seen in
FIG. 1 , as R1 and R2, both of which are located inDomain 1 and oriented in mirror symmetry to the transducer relative to a line normal to the outer surface of the casing at an angle of reflection that is equal to the incidence angle of the transducer. The extensional wave leaks off energy as it propagates, and the leaked off pressure wavefront is measured at receivers positioned further up the well. Since the extensional wave leaks energy as it propagates, the wave is being constantly attenuated as the distance increases. The materials on both sides of the casing i.e. inDomains Domain 1 is known, measuring the attenuation from one receiver to the next provides information about the acoustic properties of the material constituting the annular space (Domain 2). -
FIG. 4 shows numerical waveforms similar as expected from extensional logging in the field and illustrates the dispersion smearing out the pulse, where the first waveform is from R1 and the second waveform is from R2, as shown inFIG. 1 . The first pulse is the arrival of the emitted wavefront from the extensional wave (S0), but as seen inFIG. 3 , the group velocity is highly dispersive, which results in the pulse being smeared out. The brace above the pulse inFIG. 4 indicates the range of interest of the S0 mode, and within this range it can be seen that the pulse is preserved due to approximately constant group velocity in that frequency range (as seen inFIG. 3 at the low frequencies). The rest of the waveform is a combination of the S0 mode, the A0 mode, and the fluid borne wave (direct wave), and this region is not qualified for evaluating the attenuation of the envelope. The region where the S0 mode is being evaluated has only traveled indomains Domain 2 consists of a solid. As the extensional wave propagates along the casing, waves are constantly being leaked off. However, in the absence of a shear coupling atinterface 2, the attenuation is significantly reduced. A shear modulus inDomain 2 allows coupling of the S0 mode and the waves leaking off in both the normal direction and the tangential direction, as illustrated inFIG. 2b by the arrow in the cement layer. - The outcome is demonstrated in
FIG. 5 , whereFIG. 5a has a fluid inDomain 2 whileFIG. 5b has a solid.FIG. 5 shows the envelope of the waveforms and the time axis is minimized to demonstrate the attenuation of the early S0 mode. The waveforms inFIG. 5a are the same as inFIG. 4 , and the envelopes inFIG. 5b are also from a numerical simulation using cement acoustic properties. The double-ended arrows defining Δtn, where n=1, 2, show the predefined time intervals that are analyzed, and from the figure one can see that the attenuation is significantly stronger whenDomain 2 consists of a solid material. The integral can be estimated by the equation below, where Se,n is the envelope of the waveform and [tn,1tn,1] defines the closed interval where n=1,2 and specifies which waveform: -
E amp,n=∫tn,1 tn,2 S e,n(t)dt (1) - The amplitude of the extensional wave decreases exponentially as it propagates along the casing if the group velocity is constant. Further, since the amplitude of the leaked off pressure wavefront is related to the amplitude of the extensional wave, the leaked off wavefront will also decrease exponentially. Hence, the attenuation can be found by measuring the peak amplitude at two different distances and divide the peak measured at the further distance by the peak measured at the nearest distance, as illustrated by the receivers in
FIG. 1 . However, the S0 mode presented inFIG. 3 shows that the group velocity is constant only within a small frequency range, so a better technique is to calculate the integral within a predefined time window, as seen in the Eq. (1) and estimate the attenuation. The integral attenuation (αdB) in decibel (dB) is then found as given below: -
- Where equation (2) can be then used to estimate the loss in energy within a predefined time-window between the two waveforms.
-
FIG. 6 shows integral attenuation, estimated by using Eq. (2), plotted as a logarithmic function, whereFIG. 6a is the attenuation from the envelopes seen inFIG. 5a , andFIG. 6b is the attenuation from the envelopes seen inFIG. 5b . The horizontal axis represents the distance the S0 mode has propagated insideDomain 3. From the figure one can see that the attenuation is significantly higher (more than twice) whenDomain 2 consists of a solid. - The attenuation presented in
FIG. 6 can be quantified by dividing it by the distance propagated to estimate attenuation per centimeter (cm): -
- where zn is the distance propagated inside
Domain 3, and n=1,2 represents the index for R1 and R2. αz,dB is the attenuation expressed in dB/cm, which is plotted inFIG. 7 . Specifically, the function inFIG. 7 shows αz,dB plotted against the acoustic impedance for various different solids and fluids inDomain 2. The data seen in the figure were obtained from a study using numerical simulations, where the acoustical properties of the materials are presented in Table 1 below. -
TABLE 1 ρ νp Z νs Material (kg/m3) (m/s) (MRayl) (m/s) Light mud 1115 1670 1.68 Heavy mud 1260 1950 1.8 Foam cement 1114 2203 2.41 700 Settled barite mud 2052 1340 2.75 300 Light cement 1200 2600 3.12 1000 Shale 2500 2000 5.0 720 Class G cement 11820 4241 5.9 1700 Class G cement 21980 3350 6.63 1920 - Some materials are real values from the inventor's doctoral thesis, Sirevaag, T. PhD thesis A Study of the Ultrasonic Measurements for Logging Behind a Steel Pipe: Expanding the processing to improve the differentiation between heavy fluids and light solids, Norwegian University of Science and Technology (NTNU), which is incorporated herein by reference. It is noted that the method according to the present invention has no connection to the work in the doctoral thesis: only the material properties are used to create realistic simulations. In
FIG. 7 it is seen that solids have a significantly higher attenuation, despite having the same acoustic impedance (see foam cement and settled barite mud). Therefore, if the attenuation is greater than a predetermined threshold e.g., this is indicative of the material being a solid. When analyzing the material behind a partition, all attenuation measurements within the inspection run are used to create a model, and the minimum and maximum attenuation are used to define the outer boundaries. The threshold is set between the outer boundaries, and the value may vary from each inspection run. - Alternatively, the invention provides a complementary method of evaluating αz,dB. Instead of plotting αz,dB against the acoustic impedance, one can plot it against the shear wave velocity. This is seen in
FIG. 8 , where αz,dB increases with increasing shear wave velocity. From this, one can make a probability estimate of the likeliness of a solid inDomain 2. A higher attenuation then equals a higher probability of a solid. -
FIG. 9 shows a generic flowchart to describe the processing to estimate αz,dB between the extensional waves excited in the casing. The flowchart starts with the input of the extensional wave and knowledge of the distance between the receivers. The final output is parameter to quantify the peak envelope attenuation. -
FIG. 9 shows how to pick the beginning and end time of the time window, t1 and t2, based on the amplitude within the envelope of the S0 mode. The start time, t1 is defined as where the signal is first above the signal noise level, which in some embodiments of the invention is defined as −30 dB of the maximum amplitude within the normalized envelope of the whole waveform, although other noise thresholds may be employed. Normalizing the envelope means that the maximum amplitude is always at 0 dB. The time window needs to end before the arrival of A0 mode and/or the direct wave, so it does not influence the integral. Typically, t2 is calculated from the second peak in the envelope as seen inFIG. 5 , where the method estimates where the amplitude drops below the half maximum of the second peak. - It should be noted that modifications may be made to both the method and system as described without departing from the scope of the invention as claimed in the appended claims. In particular, it is noted that while the invention has been described with particular reference to evaluating a material on the outside of a pipe or casing buried in the ground, the invention may find more general application for discriminating between two different materials on opposite sides of a partition. This may be useful for evaluating whether liquid is present in a pipe used to transport liquid, thus allowing a blockage in the pipe to be detected. More generally, the partition does not need to be constituted by the wall of a cylindrical casing or pipe and can be any sheet material having disparate materials on opposing surfaces thereof.
- The invention has been described primarily with regard to a method because the hardware is known per se. When used to monitor oil wells, the transmitter and receivers must be inserted into a hollow casing sunk in the ground that is possibly hundreds of meters deep. Since the distances between adjacent receivers must be known, the transducers are installed into an elongated unit that is then lowered into the casing. This also allows the transducers to be tilted at known angles of incidence and reflection although more advanced transducers are available where the wave angle can be adjusted electronically. However, when applied to other situations where the transducers are more accessible, they may be spaced apart discretely at known intervals.
- It should also be noted that features that are described with reference to one or more embodiments are described by way of example rather than by way of limitation to those embodiments. Thus, unless stated otherwise or unless particular combinations are clearly inadmissible, optional features that are described with reference to only some embodiments are assumed to be likewise applicable to all other embodiments also.
- It will also be understood that the processing may be performed by a suitably programmed computer. Likewise, the invention contemplates a computer program being readable by a computer for executing the method of the invention. The invention further contemplates a machine-readable memory tangibly embodying a program of instructions executable by the machine for performing the method of the invention.
- Without derogating from the above generalizations, the inventive concepts encompassed by the invention include the following:
- Inventive concept 1: A method for evaluating a material on a remote side of a partition separating a first domain from a second domain, the method comprising:
- disposing at least one ultrasonic transmitter and a plurality of ultrasonic receivers in longitudinally spaced-apart relationship along a first side of the partition in the first domain;
- activating the at least one ultrasonic transmitter to form ultrasonic waveforms that comprise propagated quasi leaky-Lamb waves constituting extensional waves having symmetrical zero-order (S0) modes within the partition, wherein the spaced-apart receivers record the ultrasonic waveforms;
- processing the recorded ultrasonic waveforms to determine the envelope of the waves propagating through the partition from a first receiver to a second receiver that is located more remote from the transmitter than the first receiver and whose separation from the first receiver is known;
- integrating the envelope signal over a predefined interval;
- computing an integral envelope attenuation per distance of the wave propagating through the second domain; and
- establishing that the second domain contains a solid if the measured integral envelope attenuation per distance is higher than a predetermined threshold.
- Inventive concept 2: The method according to
inventive concept 1, wherein establishing that the second domain contains a solid comprises: - processing the ultrasonic waveforms to obtain the integral envelope attenuation;
- computing from the envelope a time interval within which to calculate the integral attenuation as a function of time;
- creating a model with all attenuation measurements within a complete inspection run and defining a minimum and maximum attenuation as the boundaries; and
- comparing the attenuations and setting an attenuation threshold to determine if the space contains a liquid or a solid.
- Inventive concept 3: The method of
inventive concept 2, wherein computing the time interval includes (i) defining a start time t1 as where the signal is first above a signal noise level; (ii) obtaining a second peak amplitude within the envelope (iii) defining an end time t2 as where the signal drops below half the second peak amplitude; (iv) and using the interval between t1 and t2 as boundaries for the integration of the envelope.
Inventive concept 4: The method ofinventive concept 3, wherein the integral of the envelope is calculated according to: -
E amp,n=∫tn,1 tn,2 S e,n(t)dt - where Se,n is the envelope of the waveform and [tn,1tn,1] defines a closed interval where n=1,2 and specifies which waveform.
Inventive concept 5: The method ofinventive concept 4, wherein computing the integral attenuation includes dividing the integral estimated within predefined time window of the envelope of the measured at the receiver furthest by the integral estimated within predefined time window of the envelope of the measured at the nearest receiver.
Inventive concept 6: The method according toinventive concept 5, wherein the integral envelope attenuation per distance (αz,dB) is calculated according to: -
- where:
- (Eamp,n) is the integral of the envelope of the extensional pulse at the nth receiver;
- (Pm,1) is the integral of the envelope of the extensional pulse at the first receiver;
- Zn−Z1 the distance between the receivers.
- Inventive concept 7: The method according to
inventive concept 3, wherein the integral of the envelope amplitude is estimated from a portion of the time domain wherein the symmetric zero-order wavefront (S0) arrives before a higher energy antisymmetric wavefront (A0) is able to reach said portion thereby eliminating the effect of interference between the symmetric and the antisymmetric wavefronts. - Inventive concept 8: The method according to any one of
inventive concepts 3 to 7, wherein the predetermined noise threshold is −30 dB of the peak amplitude within the normalized envelope of the whole waveform.
Inventive concept 9: The method according to any one ofinventive concepts 2 to 8, further including computing the integral envelope attenuation per distance (αz,dB) against a probability chart to determine the likelihood of the second domain containing a solid.
Inventive concept 10: The method according toclaim 9, wherein if the attenuation is greater than a predetermined threshold, this is indicative of the material being a solid.
Inventive concept 11: The method according toinventive concept 1, wherein establishing that the second domain contains a solid comprises: - identifying Eamp,1 and Eamp,2 in the measured envelope waveforms, where (Eamp,n) is the integral of the envelope of the extensional pulse at the nth receiver;
- determining the integral attenuation;
- calculating the integral envelope attenuation per distance; and
- comparing the integral envelope attenuation per distance to a predetermined attenuation threshold to determine if the space contains a liquid or a solid.
- Inventive concept 12: The method according to inventive concept 11, wherein the predetermined attenuation threshold is in the order of 7 dB/cm.
Inventive concept 13: The method according to any one of the preceding inventive concepts, wherein the partition comprises a casing and an annular fill installed in a borehole traversing a formation.
Inventive concept 14: The method according to any one of the preceding inventive concepts, wherein each transmitter and receiver is directed toward the partition at an angle of incidence close to 15°.
Inventive concept 15: A computer program product comprising a computer readable medium storing program code instructions, which when executed on at least one processor that receives as input data representative of the recorded ultrasonic waveforms from a pair of spaced apart receivers, carries out the processing and computing operations of any one of the preceding inventive concepts.
Claims (15)
1. A method for evaluating a material between a wellbore and a casing, said casing having a wall forming a partition that separates a first domain inside the casing from a second domain outside the casing, the method comprising:
disposing at least one ultrasonic transmitter and a plurality of ultrasonic receivers in longitudinally spaced-apart relationship along a first side of the partition in the first domain;
activating the at least one ultrasonic transmitter to form ultrasonic waveforms that comprise propagated quasi leaky-Lamb waves constituting extensional waves having symmetrical zero-order (S0) modes within the partition, wherein the spaced-apart receivers record the ultrasonic waveforms;
processing the recorded ultrasonic waveforms to determine an envelope of the waves propagating through the partition from a first receiver to a second receiver that is located more remote from the transmitter than the first receiver and whose separation from the first receiver is known;
integrating the envelope over a predefined interval;
computing an integral envelope attenuation per distance of the wave propagating through the second domain; and
establishing that the second domain contains a solid if the measured integral envelope attenuation per distance is higher than a predetermined threshold.
2. The method according to claim 1 , wherein establishing that the second domain contains a solid comprises:
processing the ultrasonic waveforms to obtain the integral envelope attenuation;
computing from the envelope a time interval within which to calculate the integral attenuation as a function of time;
from all attenuation measurements within a complete inspection run, defining a minimum and maximum attenuation as the boundaries; and
comparing the attenuations and setting an attenuation threshold to determine if the space contains a liquid or a solid.
3. The method of claim 2 , wherein computing the time interval includes (i) defining a start time t1 as where the signal is first above a signal noise level; (ii) obtaining a second peak amplitude within the envelope (iii) defining an end time t2 as where the signal drops below half the second peak amplitude; (iv) and using the interval between t1 and t2 as boundaries for the integration of the envelope.
4. The method of claim 3 , wherein the integral of the envelope is calculated according to:
E amp,n=∫tn,1 t n,2 S e,n(t)dt
E amp,n=∫t
where Se,n is the envelope of the waveform and [tn,1tn,1] defines a closed interval where n=1,2 and specifies which waveform.
5. The method of claim 4 , wherein computing the integral attenuation includes dividing the integral estimated within predefined time window of the envelope of the measured at the receiver furthest by the integral estimated within predefined time window of the envelope of the measured at the nearest receiver.
6. The method according to claim 5 , wherein the integral envelope attenuation per distance (αz,dB) is calculated according to:
where:
(Eamp,n) is the integral of the envelope of the extensional pulse at the nth receiver;
(Pm,1) is the integral of the envelope of the extensional pulse at the first receiver;
Zn−Z1 the distance between the receivers.
7. The method according to claim 3 , wherein the integral of the envelope amplitude is estimated from a portion of the time domain wherein the symmetric zero-order wavefront (S0) arrives before a higher energy antisymmetric wavefront (A0) is able to reach said portion thereby eliminating the effect of interference between the symmetric and the antisymmetric wavefronts.
8. The method according to claim 3 , wherein the predetermined noise threshold is −30 dB of the peak amplitude within the normalized envelope of the whole waveform.
9. The method according to claim 1 , wherein establishing that the second domain contains a solid comprises:
processing the ultrasonic waveforms to obtain the integral envelope attenuation; and
computing the integral envelope attenuation per distance (αz,dB) against a probability chart to determine the likelihood of the second domain containing a solid.
10. The method according to claim 2 , wherein if the attenuation is greater than said attenuation threshold, this is indicative of the material being a solid.
11. The method according to claim 1 , wherein establishing that the second domain contains a solid comprises:
identifying Eamp,1 and Eamp,2 in the measured envelope waveforms, where (Eamp,n) is the integral of the envelope of the extensional pulse at the nth receiver;
determining the integral attenuation;
calculating the integral envelope attenuation per distance; and
comparing the integral envelope attenuation per distance to a predetermined attenuation threshold to determine if the space contains a liquid or a solid.
12. The method according to claim 11 , wherein the predetermined attenuation threshold is in the order of 7 dB/cm.
13. The method according to claim 1 , wherein the partition comprises a casing and an annular fill installed in a borehole traversing a formation.
14. The method according to claim 1 , wherein each transmitter and receiver is directed toward the partition at an angle of incidence equal to 15°.
15. A non-transitory computer program product comprising a computer readable medium storing program code instructions, which when executed on at least one processor that receives as input data representative of the recorded ultrasonic waveforms from a pair of spaced apart receivers, carries out the processing and computing operations of claim 1 .
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
GB2003621.6 | 2020-03-12 | ||
GB2003621 | 2020-03-12 | ||
GB2003621.6A GB2592974A (en) | 2020-03-12 | 2020-03-12 | Method for differentiating materials on a remote side of a partition based on the attenuation of the ultrasonic extensional zero mode |
Publications (2)
Publication Number | Publication Date |
---|---|
US11098583B1 US11098583B1 (en) | 2021-08-24 |
US20210285324A1 true US20210285324A1 (en) | 2021-09-16 |
Family
ID=70453632
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US17/159,341 Active US11098583B1 (en) | 2020-03-12 | 2021-01-27 | Method for determining the integrity of a solid bonding between a wellbore and a casing |
Country Status (3)
Country | Link |
---|---|
US (1) | US11098583B1 (en) |
EP (1) | EP3879311B1 (en) |
GB (1) | GB2592974A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20230082439A1 (en) * | 2021-09-16 | 2023-03-16 | Kabushiki Kaisha Toshiba | Inspection system and inspection method |
US20230194743A1 (en) * | 2021-12-20 | 2023-06-22 | Halliburton Energy Services, Inc. | Through tubing cement evaluation based on casing extensional waves |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11378707B2 (en) * | 2020-11-12 | 2022-07-05 | Schlumberger Technology Corporation | Third interface echo (TIE) detection from flexural data for gas/liquid annulus discrimination |
US11970931B2 (en) * | 2021-06-01 | 2024-04-30 | Halliburton Energy Services, Inc. | Through tubing cement evaluation using borehole resonance mode |
GB2609627A (en) | 2021-08-09 | 2023-02-15 | Equanostic As | Method for determining if a wellbore consists of micro annulus, free pipe or solid bonding between the wellbore and a casing |
US12032113B2 (en) | 2022-01-03 | 2024-07-09 | Halliburton Energy Services, Inc. | Through tubing cement evaluation based on rotatable transmitter and computational rotated responses |
GB202213172D0 (en) | 2022-09-08 | 2022-10-26 | Equanostic As | Dual-sensor head for pipe inspection robot and method of use |
Family Cites Families (42)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3401773A (en) * | 1967-12-04 | 1968-09-17 | Schlumberger Technology Survey | Method and apparatus for cement logging of cased boreholes |
US4255798A (en) * | 1978-05-30 | 1981-03-10 | Schlumberger Technology Corp. | Method and apparatus for acoustically investigating a casing and cement bond in a borehole |
GB2008756B (en) * | 1977-09-30 | 1982-07-14 | Ti Group Services Ltd | Ultrasonic inspection |
US6018496A (en) * | 1991-12-17 | 2000-01-25 | Schlumberger Technology Corporation | Method and apparatus for hydraulic isolation determination |
US5357481A (en) * | 1992-11-04 | 1994-10-18 | Western Atlas International, Inc. | Borehole logging tool |
US5591945A (en) * | 1995-04-19 | 1997-01-07 | Elo Touchsystems, Inc. | Acoustic touch position sensor using higher order horizontally polarized shear wave propagation |
CA2316265C (en) * | 1998-01-06 | 2008-03-25 | Schlumberger Canada Limited | Method and apparatus for ultrasonic imaging of a cased well |
US6681185B1 (en) * | 2002-07-26 | 2004-01-20 | Eseis | Method of seismic signal processing |
DK1505252T3 (en) * | 2003-08-08 | 2006-01-30 | Schlumberger Technology Bv | Acoustic multimode imaging in lined wells |
US7697375B2 (en) * | 2004-03-17 | 2010-04-13 | Baker Hughes Incorporated | Combined electro-magnetic acoustic transducer |
US7663969B2 (en) | 2005-03-02 | 2010-02-16 | Baker Hughes Incorporated | Use of Lamb waves in cement bond logging |
EP1672168B1 (en) * | 2004-12-20 | 2008-02-06 | Services Petroliers Schlumberger | Determination of the impedance of a material behind a casing combining two sets of ultrasonic measurements |
US7681450B2 (en) * | 2005-12-09 | 2010-03-23 | Baker Hughes Incorporated | Casing resonant radial flexural modes in cement bond evaluation |
US7773454B2 (en) * | 2006-02-22 | 2010-08-10 | Baker Hughes Incorporated | Method and apparatus for cement evaluation using multiple acoustic wave types |
US7787327B2 (en) * | 2006-11-15 | 2010-08-31 | Baker Hughes Incorporated | Cement bond analysis |
US20090231954A1 (en) * | 2008-03-17 | 2009-09-17 | Baker Hughes Incorporated | Micro-Annulus Detection Using Lamb Waves |
US8964504B2 (en) * | 2010-04-07 | 2015-02-24 | Baker Hughes Incorporated | Method and apparatus for evaluating a cemented borehole casing |
AU2011382521B2 (en) * | 2011-11-30 | 2014-10-30 | Halliburton Energy Services, Inc. | Acoustic transducer apparatus, systems, and methods |
US20140056111A1 (en) * | 2012-08-21 | 2014-02-27 | Cung Khac Vu | Acoustic detector |
US9273545B2 (en) * | 2012-12-23 | 2016-03-01 | Baker Hughes Incorporated | Use of Lamb and SH attenuations to estimate cement Vp and Vs in cased borehole |
US9534487B2 (en) | 2014-01-16 | 2017-01-03 | Schlumberger Technology Corporation | Cement acoustic properties from ultrasonic signal amplitude dispersions in cased wells |
US9784875B2 (en) * | 2014-01-31 | 2017-10-10 | Schlumberger Technology Corporation | Method to estimate cement acoustic wave speeds from data acquired by a cased hole ultrasonic cement evaluation tool |
US10138727B2 (en) | 2014-01-31 | 2018-11-27 | Schlumberger Technology Corporation | Acoustic multi-modality inversion for cement integrity analysis |
US9772419B2 (en) * | 2014-03-10 | 2017-09-26 | Read As | Decomposing full-waveform sonic data into propagating waves for characterizing a wellbore and its immediate surroundings |
US9732607B2 (en) * | 2014-08-18 | 2017-08-15 | Schlumberger Technology Corporation | Methods and apparatus for evaluating properties of cement utilizing ultrasonic signal testing |
US20160209539A1 (en) * | 2014-11-14 | 2016-07-21 | Schlumberger Technology Corporation | Method for Separating Multi-Modal Acoustic Measurements for Evaluating Multilayer Structures |
GB2531836B (en) * | 2014-12-24 | 2020-10-14 | Equinor Energy As | Logging system and method for evaluation of downhole installation |
WO2016105207A1 (en) * | 2014-12-24 | 2016-06-30 | Statoil Petroleum As | Evaluation of downhole installation |
WO2016187240A1 (en) * | 2015-05-18 | 2016-11-24 | Schlumberger Technology Corporation | Method for analyzing cement integrity in cased wells using sonic logging |
US10858933B2 (en) * | 2015-05-18 | 2020-12-08 | Schlumberger Technology Corporation | Method for analyzing cement integrity in casing strings using machine learning |
US10809405B2 (en) * | 2015-07-06 | 2020-10-20 | Schlumberger Technology Corporation | Measurement and processing to detect weak interfacial layers in hydrocarbon-bearing laminated formations with acoustic logging devices |
EP3151037A1 (en) * | 2015-09-30 | 2017-04-05 | Services Pétroliers Schlumberger | Systems and methods for evaluating annular material using beamforming from acoustic arrays |
EP3151038A1 (en) * | 2015-09-30 | 2017-04-05 | Services Pétroliers Schlumberger | Systems, methods, and computer-readable media for determining shear-wave to compressional-wave velocity ratios in well casings |
US10761064B2 (en) * | 2015-10-02 | 2020-09-01 | Schlumberger Technology Corporation | Detection and characterization of thin weak bedding planes in hydrocarbon-bearing laminated formations |
EP3168412A1 (en) * | 2015-11-16 | 2017-05-17 | Services Pétroliers Schlumberger | Cement evaluation using the integration of multiple modes of acoustic measurements |
EP3290961A1 (en) * | 2016-09-06 | 2018-03-07 | Services Pétroliers Schlumberger | Separation of flexural and extensional modes in multi modal acoustic signals |
US11531132B2 (en) * | 2017-03-24 | 2022-12-20 | Schlumberger Technology Corporation | Guided mode beamforming for probing open-hole and cased-hole well environments |
US11131182B2 (en) * | 2017-03-27 | 2021-09-28 | Schlumberger Technology Corporation | Methods for characterizing multi-string cased wells using wide frequency bandwidth signals |
EP3404202A1 (en) * | 2017-05-15 | 2018-11-21 | Services Petroliers Schlumberger | Flexural wave measurement for thick casings |
US11220897B2 (en) * | 2018-04-12 | 2022-01-11 | Schlumberger Technology Corporation | Evaluating casing cement using automated detection of clinging compression wave (P) arrivals |
US11719090B2 (en) * | 2019-03-22 | 2023-08-08 | Baker Hughes Oilfield Operations Llc | Enhanced cement bond and micro-annulus detection and analysis |
GB2588798B (en) * | 2019-11-07 | 2022-11-16 | Equanostic As | Method for evaluating a material on a remote side of a partition using ultrasonic measurements |
-
2020
- 2020-03-12 GB GB2003621.6A patent/GB2592974A/en active Pending
-
2021
- 2021-01-27 US US17/159,341 patent/US11098583B1/en active Active
- 2021-01-29 EP EP21154203.0A patent/EP3879311B1/en active Active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20230082439A1 (en) * | 2021-09-16 | 2023-03-16 | Kabushiki Kaisha Toshiba | Inspection system and inspection method |
US20230194743A1 (en) * | 2021-12-20 | 2023-06-22 | Halliburton Energy Services, Inc. | Through tubing cement evaluation based on casing extensional waves |
Also Published As
Publication number | Publication date |
---|---|
GB2592974A (en) | 2021-09-15 |
GB202003621D0 (en) | 2020-04-29 |
EP3879311A1 (en) | 2021-09-15 |
US11098583B1 (en) | 2021-08-24 |
EP3879311B1 (en) | 2022-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11098583B1 (en) | Method for determining the integrity of a solid bonding between a wellbore and a casing | |
US9784875B2 (en) | Method to estimate cement acoustic wave speeds from data acquired by a cased hole ultrasonic cement evaluation tool | |
US9829597B2 (en) | Model based inversion of acoustic impedance of annulus behind casing | |
US10345465B2 (en) | Resonance-based inversion of acoustic impedance of annulus behind casing | |
US10012749B2 (en) | Fast model based inversion of acoustic impedance of annulus behind casing | |
AU2021202822B2 (en) | A method of identifying a material and/or condition of a material in a borehole | |
US10344582B2 (en) | Evaluation of downhole installation | |
US9322807B2 (en) | Ultrasonic signal time-frequency decomposition for borehole evaluation or pipeline inspection | |
US11143016B2 (en) | Method for evaluating a material on a remote side of a partition using ultrasonic measurements | |
US10191173B2 (en) | Systems and methods for evaluating annular material using beamforming from acoustic arrays | |
US20160209539A1 (en) | Method for Separating Multi-Modal Acoustic Measurements for Evaluating Multilayer Structures | |
WO2009055209A2 (en) | Measurement of sound speed of downhole fluid utilizing tube waves | |
US20150285936A1 (en) | System and Methodology for Determining Fracture Attributes in A Formation | |
US11566517B2 (en) | Quantifying cement bonding quality of cased-hole wells using a quality index based on frequency spectra | |
US10393905B2 (en) | Torsional wave logging | |
GB2528326A (en) | Method of determining a condition of a borehole and apparatus | |
US11598895B2 (en) | Ultrasonic waveform processing using deconvolution in downhole environments | |
Sirevaag et al. | Laboratory Setup for Improved Behind Casing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO SMALL (ORIGINAL EVENT CODE: SMAL); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
AS | Assignment |
Owner name: EQUANOSTIC AS, NORWAY Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:SIREVAAG, TORE;REEL/FRAME:055438/0140 Effective date: 20210121 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |