This application is a U.S. National Stage Filing under 35 U.S.C. 371 from International Application No. PCT/US2014/034314, filed on 16 Apr. 2014; which application is incorporated herein by reference in its entirety.
After drilling operations, borehole casings can be used to help stabilize the well. The casing is cemented to assure zonal isolation and sustain its stability. Over time deterioration and/or imperfections in the casing and/or cement can cause the borehole to collapse and/or fluids to migrate therefore compromising the zonal isolation.
Fluid (e.g., oil, gas) pipelines can be used to transport the fluids for long distances from one location to another. Environmental factors and fluid chemicals can result in corrosion of the pipeline inside and outside, thus weakening the pipeline walls. Improper welding of the pipeline sections can also result in leaking of the fluid into the surrounding environment. It can thus be desirable to be able to determine the integrity of a borehole casing or a fluid pipeline.
Ultrasound wave propagation can be used as part of an interrogation method for borehole evaluation or fluid pipeline inspection. The interrogation method can be used in either a pulse-echo configuration or a pitch-catch configuration.
In order to obtain borehole information from measured ultrasonic signals, such as casing inside diameter (ID), casing thickness, cement bonding impedance, and a borehole image, the identity or estimate of discrete reflected/reverberation echoes from the measured ultrasonic signals can be determined. Inspection of pipeline integrity, including pipeline ID, pipeline thickness, corrosion, and weld condition, can be accomplished by separating ultrasonic echoes and/or wave packages from the measured signals. Embodiments disclosed herein use ultrasonic signal time-frequency decomposition (e.g., a matching pursuit signal decomposition, wavelet transform method), in order to decompose a measured ultrasonic signal into discrete wavelets. Each discrete wavelet can represent values of arrival time, frequency, and amplitude for received physical echoes (e.g., wave packages). The arrival time and/or frequency and/or amplitude represented by the wavelets can then be used to calculate information regarding the physical condition of the casing or pipeline and the cement casing bonding impedance. The decomposed wavelets can also be used to compress the measured waveforms for the purposes of data storage and telemetry. The compressed waveforms can then be stored.
As used herein, a wavelet can be defined as a received echo signal that results from a transmitted ultrasonic signal. The echo signal (e.g., wavelet) may have bounced back from the inside wall of a casing or pipeline, the outside wall of the casing or pipeline, and/or imperfections in the casing or pipeline. For example,
The received ultrasonic signal can then be decomposed 204. The decomposition can be a time-frequency decomposition, such as matching pursuit or traditional wavelet decomposition. The decomposition can have a set of basis functions that are substantially similar to the physical echoes (e.g., wavelets) in waveform shape.
Matching pursuit is a signal decomposition method whereby a signal is decomposed into a linear combination of basis functions that are selected from a redundant dictionary. The purpose of using redundant basis functions is to maximize the waveform matching between the basis functions and the signal, which may not be achieved if the basis functions are restricted by orthogonality, as occurs with discrete wavelet decomposition.
The following description of the decomposition of the received signal into wavelets uses Gabor functions for purposes of illustration only. The present embodiments are not limited to any one basis function type for decomposing the received signal. Other embodiments can use more than one basis function type.
A family of basis functions (e.g., Gabor functions) are shown subsequently, and each received signal can be matched with the best function from the family of functions. The family of modified Gabor functions can be used as the basis functions for methods described herein since Gabor functions can be qualitatively and quantitatively very similar to ultrasonic echoes and wave packages. Modified Gabor functions retain the waveform shape of their respective Gabor functions, but render the proposed decomposition translation invariant (i.e., a shift of the original signal in time does not change the decomposed wavelets except for the same amount of shift in time in each decomposed wavelet). This is a desirable property for stable decomposition of a measured signal in the presence of measurement errors and noise.
Each modified Gabor function is a product of a Gaussian function with an in-phase, complex sinusoid and is specified by arrival time, scaling factor, and frequency, denoted by u, s, and ω, respectively. For example, one such Gabor function can be represented by:
The M decomposed wavelets can be denoted by (Am, φm, um, sm, ωm) where mεM 205 if the measured signal is decomposed into M wavelets where (Am, φm) are the weights (amplitude and phase) for each decomposed wavelet.
By performing a cross-correlation between each wavelet and the Gabor basis function, the arrival time (u) of the wavelet at the receiver, the phase (φ) of the wavelet, and the amplitude (A) of the wavelet can be determined. These parameters (A, φ, u) can be used for borehole evaluation and pipeline inspection.
Referring again to
The (Ai, φi, ui) parameters of reverberation wavelets (e.g., see reflected echoes 102-107 of
R=R12+T12R23T21e−i2k
where Rij and Tij are the reflection and transmission coefficients, respectively; d refers to the casing thickness and kc is the wavenumber. The i and j indices (i.e. 1, 2 and 3) refer to the layers of fluid, steel and cement respectively. Each R and T is defined by the impedance contrast between the layers. As the impedance contrast between cement and steel changes, the rate of decay of the reverberations amplitudes and energy change. Thus, to determine the impedance, the wavelet is constrained and the energy and estimate of the energy decay is determined.
The (Ai, φi, ui) parameters of guided wave arrivals have been determined 230 as discussed previously. The guided wave arrivals can result from transmission of the ultrasonic signal within the wall of a pipeline. Because the pipe wall is a wave guide, the waves propagating within the pipeline wall can be referred to as guided waves and the properties of the guided waves can be affected by the wave guide (i.e., the pipe wall). This is the principle that can be used for pipeline inspection 231. The matching pursuit method above can be used to extract useful information (e.g., Ai, φi, ui) from the received guided waves to infer the integrity of the pipe. For example, the pipeline inspection can include pipeline welding integrity (i.e., determined according to reflected arrivals of guided wavelets) and pipeline corrosion (i.e., determined according to speed change in the guided wavelets).
The M decomposed wavelet parameters (Am, φm, um, sm, ωm), where mεM, have been determined 240 as discussed previously. The wavelet parameters can be stored and/or propagated via telemetry 241 to be used in analyzing the wavelets by one or more processors located in a wireline tool or on the surface. As discussed previously, u is the arrival time of a wavelet at the receiver, φ is the phase of the wavelet at the receiver, A is the amplitude of the wavelet at the receiver, ω is the frequency of the wavelet, and s is a scaling factor.
The method calculates a synthetic reverberation signal 300 using a theoretical model based on certain prior information. For example, the prior information can include a physical model of the pipeline, including a layer of fluid in the pipe followed by a layer of steel followed by a layer of fluid outside the pipe. Each layer can be characterized by specific acoustic properties including speed and impedance. The physical model may also incorporate a specific radiation pattern of the transducer and curvature of the pipe. The synthetic signal can then be decomposed into M wavelets 303. In an embodiment, the decomposition can be performed using matching pursuit decomposition as previously described. The resulting M decomposed wavelet parameters include (Am, φm, um, sm, ωm), where mεM 305. The parameters (φm, sm, ωm), where mεM, of the decomposed wavelets, define their respective waveform shapes and (Am, um), where mεM, define the magnitudes and arrival times.
The phase, scaling factor, and frequency (φm, sm, ωm), where mεM, from the synthetic signal are determined to be constraint parameters 307. When a measured signal is received 310, the constraint parameters can be used during the decomposition of the measured signal 312 so that the shapes of the basis function are kept the same for the decomposition of the measured signal. Such constraints can be indicated as: (φ′m=φm, s′m=sm, ω′m=ωm), where mεM. With these constraints, the decomposition leads to a set of decomposed wavelets 314 that can be quantified by (A′m, φm, u′m, sm, φm), where mεM, where only the amplitude (A′m) and arrival time (u′m) are fitted from the measured signal.
In another embodiment, an additional constraint, |u′m−um|=δm, can be placed on the time shift u′m of each modified Gabor function. The value of δm is set to prevent a matching Gabor function from jumping to another portion of the measure signal, so that its value stays in the vicinity of values produced by the original Gabor function. The determined values may then be used for pipeline and borehole inspection.
The inspection tool 402 may comprise an ultrasonic tool with a housing 404 coupled to a rotatable mount 424 for borehole evaluation (e.g., casing 408 and cement 412 condition) or pipeline inspection. The rotatable mount 424 is attached in turn to a pulse-echo device 440 (e.g., a pulse-echo circumferential acoustic scanning tool (CAST) transducer), and a pitch-catch array 416 that, for purposes of illustration, includes two transmitters 448 and five receivers 452. Other embodiments may use either the pulse-echo device 400 or the pitch-catch device 416.
The rotatable mount 424 is also attached to the apparatus 400, which comprises another pulse-echo device 444 (e.g., a second pulse-echo transducer) for velocity measurement in a fluid 436 (e.g., mud), and a dual-faced target 456 that provides at least two different propagation distances through the fluid 436, which can be averaged to provide a better measurement of the fluid velocity.
The pitch-catch array 416 and the pulse-echo device 440 can receive reflected signals (e.g., measured signals) resulting from transmitted ultrasonic signals. The measured signal can be decomposed by the above-described method in order to perform the borehole evaluation or pipeline inspection.
The array 416 can include transmitters 448 and receivers 452 that extend along the length of the axis of rotation L, which is substantially aligned with the longitudinal axis of the housing 404. The rotatable mount 424 (e.g., an in-line rotatable head) is shown with in-line centralizers 428, which allows the rotatable mount 424 to be located at any point in a logging tool string. One or more inline centralizers 428 may be disposed between the housing 404 and the rotatable mount 424. The rotatable mount 424 can comprise a framework 432 having two or more arms (two arms are shown—a first arm 432′ and a second arm 432″) approximately horizontally opposed across the longitudinal axis of the housing 404. The pulse-echo devices 440, 444 may be attached to the first arm 432′, and the array 416 may be attached to the second arm 432″.
The weight of the individual arms 432′, 432″ may be adjusted to balance the weight of various components (e.g., transducers, such as the array 416) attached to the arms 432′, 432″, to impart rotational stability to the rotatable mount 424. The arms 432′, 432″ may be provided with a fluid dynamic profile design to reduce drag in the borehole fluid 436 while rotating. The framework 432 may permit radial expansion and contraction, to operate more easily in different sizes of the borehole.
By using the array 416, cement evaluation with imaging behind the casing 408 may be accomplished more easily. Both pitch-catch (e.g., in the array 416) and pulse-echo devices 440 (e.g., a pulse-echo transducer) can be used to assess the integrity of the cement 412 behind the casing 408. For example, the pulse-echo device 440 can be used to estimate the thickness of the casing 408, the borehole shape, and the impedance of the cement 412 behind the casing 408 (from the normal incidence of the ultrasonic wave). The pitch-catch array 416 can be used to help identify reflectors behind the casing 408 as part of imaging behind the casing 408, including reflectors that are not normal to the incidence wave provided by the pulse-echo device 440.
Array focusing may be used to further enhance the received reflected signals. In some embodiments, the transmitters 448 operate with a center frequency of about 50 kHz to about 150 kHz, or a broad band (impulse or chirp) pulse from about 50 kHz to about 200 kHz. The response of the receivers 452 may be selected to peak over a wider range, of about 20 kHz to about 250 kHz.
While a rotatable mount 424 is shown, the pitch-catch concept can also be applied to a non-rotating segmented transducer array assembly to image behind casing. In this case, a non-rotating array of receivers 452 with five azimuthal segments can map the interface between the cement 412 and the formation 420, although perhaps with less resolution than the rotatable mount 424 provides.
Additional transducers, including an array, may be mounted to the same framework 432 as the rotating pulse-echo device 444 that is used to determine the fluid velocity. Thus, the housing 404 may be attached to an array 416 of pitch-catch transducers mounted to a second arm 432″ of a framework 432 forming a part of the rotatable mount 424, the pitch-catch transducers facing away from the axis of rotation L.
The array 416 may include multiple transmitters and receivers, wherein the receivers are shared by the transmitters. As shown in
The pulse-echo device 444 and the array 416 may be mounted on opposite sides of the framework. Thus, the housing 404 may be fashioned so that the pulse-echo device 444 (e.g., an acoustic transducer) is mounted to a first arm 432′ of the framework 432, the first arm 432′ substantially opposing the second arm 432″ with respect to the axis of rotation L.
Yet another device 440, perhaps operating in pulse-echo mode, may be mounted to the rotating framework 432, facing away from the axis of rotation L. Thus, the tool 402 may comprise a pulse-echo device 440 (e.g., an acoustic transducer) mounted to the first arm 432′ of the framework 432, the pulse-echo device 440 facing away from the axis of rotation L. The pulse-echo device 440 may comprise a wide band transmitter so that a single transducer can be applied for a wide range of casing sizes. Further embodiments may be realized.
Other components of the system 502 may include memory 550 and data acquisition logic 540. The memory 550 is communicatively coupled to the processor(s) 530. The processor(s) 530, memory 550, and logic 540 may form part of a data acquisition system, where a database 534 is used to store waveform telemetry that results from the above-described method for ultrasonic signal decomposition for borehole evaluation or pipeline inspection.
The processor(s) 530 can execute instructions stored in the memory to receive a reflected signal from an acoustic receiver, decompose the reflected signal into at least one wavelet using the time-frequency decomposition and set of basis functions as described previously, determine parameters associated with each wavelet, and perform the borehole evaluation or pipeline inspection based on the parameters. The decomposition can be performed at a surface processor once data from the tool is received.
A data transmitter 524 may be used to transmit acquired data values and/or processing results to the surface 566. Thus, the system 502 may comprise a data transmitter 524 (e.g., a telemetry transmitter) to transmit acquired echo data to a surface workstation 556. The data transmitter 524 can be either wireless or wireline (e.g., wire, fiber optics). Additional embodiments may be realized.
Drilling oil and gas wells is commonly carried out using a string of drill pipes connected together so as to form a drilling string that is lowered through a rotary table 610 into a wellbore or borehole 612. Here it is assumed that the drilling string has been temporarily removed from the borehole 612 to allow a wireline logging tool body 670, such as a probe or sonde, to be lowered by wireline or logging cable 674 into the borehole 612. Typically, the wireline logging tool body 670 is lowered to the bottom of the region of interest and subsequently pulled upward at a substantially constant speed.
During an upward or downward trip, at a series of depths, the instruments (e.g., the apparatus 400 of
In some embodiments, the tool body 670 comprises the apparatus 400 for obtaining and analyzing acoustic noise measurements from a subterranean formation through a borehole. The tool body 670 is suspended in the wellbore by a wireline cable 674 that connects the tool body 670 to a surface control unit (e.g., comprising a workstation 654). The tool body 670 with the apparatus 400 may be deployed in the borehole 612 on coiled tubing, jointed drill pipe, hard wired drill pipe, or any other suitable deployment technique.
The bottom hole assembly 720 may include drill collars 722, a down hole tool 724, and a drill bit 726. The drill bit 726 may operate to create the borehole 712 by penetrating the surface 704 and the subsurface formations 714. The downhole tool 724 may comprise any of a number of different types of tools including MWD tools, LWD tools, and others.
During drilling operations, the drill string 708 (perhaps including the Kelly 716, the drill pipe 718, and the bottom hole assembly 720) may be rotated by the rotary table 710. Although not shown, in addition to, or alternatively, the bottom hole assembly 720 may also be rotated by a motor (e.g., a mud motor) that is located down hole. The drill collars 722 may be used to add weight to the drill bit 726. The drill collars 722 may also operate to stiffen the bottom hole assembly 720, allowing the bottom hole assembly 720 to transfer the added weight to the drill bit 726, and in turn, to assist the drill bit 726 in penetrating the surface 704 and subsurface formations 714.
During drilling operations, a mud pump 732 may pump drilling fluid (sometimes known by those of ordinary skill in the art as “drilling mud”) from a mud pit 734 through a hose 736 into the drill pipe 718 and down to the drill bit 726. The drilling fluid can flow out from the drill bit 726 and be returned to the surface 704 through an annular area 740 between the drill pipe 718 and the sides of the borehole 712. The drilling fluid may then be returned to the mud pit 734, where such fluid is filtered. In some embodiments, the drilling fluid can be used to cool the drill bit 726, as well as to provide lubrication for the drill bit 726 during drilling operations. Additionally, the drilling fluid may be used to remove subsurface formation cuttings created by operating the drill bit 726.
The downhole tool 724 can comprise the apparatus 400 as discussed previously with reference to
For example, in some embodiments, a system 502, 664, 764 may comprise a housing and one or more apparatus 400 attached to the housing with a rotatable mount associated with an axis of rotation (e.g., the axis L shown in
The housing can be used in wireline or MWD/LWD applications. Thus, the housing may comprise one of a wireline tool body, a measurement while drilling downhole tool, or a logging while drilling downhole tool.
Such modules may include hardware circuitry, and/or a processor and/or memory circuits, software program modules and objects, and/or firmware, and combinations thereof, as desired by the architect of the apparatus 400 and systems 502, 664, 764 and as appropriate for particular implementations of various embodiments. For example, in some embodiments, such modules may be included in an apparatus and/or system operation simulation package, such as a signal propagation simulation package, a power usage and distribution simulation package, a power/heat dissipation simulation package, and/or a combination of software and hardware used to simulate the operation of various potential embodiments.
It should also be understood that the apparatus and systems of various embodiments can be used in applications other than for logging operations, and thus, various embodiments are not to be so limited. The illustrations of apparatus 400 and systems 502, 664, 764 are intended to provide a general understanding of the structure of various embodiments, and they are not intended to serve as a complete description of all the elements and features of apparatus and systems that might make use of the structures described herein.
The systems 502, 664, 764 can include a non-transitory machine-readable storage device can comprise instructions stored thereon, which, when performed by a machine, cause the machine to perform operations, the operations comprising one or more features similar to or identical to features of methods and techniques related to ultrasonic signal decomposition for borehole evaluation or pipeline inspection. These operations include any one or all of the operations forming the methods shown in
A machine-readable storage device, herein, is a physical device that stores data represented by physical structure within the device. Examples of non-transitory machine-readable storage devices can include, but are not limited to, read only memory (ROM), random access memory (RAM), a magnetic disk storage device, an optical storage device, a flash memory, and other electronic, magnetic, and/or optical memory devices.
Traditional methods to determine the thickness of a casing is to use a peak frequency of a recorded wave reverberation. This approach is typically limited to thin casing and relative large diameter pipelines. For evaluating the cement-casing bond, the traditional cement bond log uses the normal-incident plane-wave theory to calculate the acoustic impedance behind the casing. This theory does not consider either the curvature of the casing wall or the transducer's beam spreading. In addition, the cement bond impedance is calculated based on a resonance response. This can limit the area of applicability of traditional methods to relatively thin casings.
The apparatus, systems, and methods disclosed herein may provide a new signal analysis approach for ultrasonic borehole evaluation and pipeline inspection based on signal decomposition. The adopted translation-invariant matching pursuit decomposition decomposes a measured ultrasonic waveform into distinct reflection and reverberation echoes and/or wave packages (e.g., wavelets). This can be accomplished robustly in the presence of measurement errors and noise. The decomposed wavelets can be directly used for a variety of purposes in borehole evaluation and pipeline inspection, including casing/pipeline ID, casing/pipeline thickness, cement bonding impedance, borehole image, pipeline corrosion, and pipeline welding condition.
The accompanying drawings that form a part hereof, show by way of illustration, and not of limitation, specific embodiments in which the subject matter may be practiced. The embodiments illustrated are described in sufficient detail to enable those skilled in the art to practice the teachings disclosed herein. Other embodiments may be utilized and derived therefrom, such that structural and logical substitutions and changes may be made without departing from the scope of this disclosure. This Detailed Description, therefore, is not to be taken in a limiting sense, and the scope of various embodiments is defined only by the appended claims, along with the full range of equivalents to which such claims are entitled.
In the foregoing Detailed Description, it can be seen that various features are grouped together in a single embodiment for the purpose of streamlining the disclosure. This method of disclosure is not to be interpreted as reflecting an intention that the claimed embodiments require more features than are expressly recited in each claim. Rather, as the following claims reflect, inventive subject matter lies in less than all features of a single disclosed embodiment. Thus the following claims are hereby incorporated into the Detailed Description, with each claim standing on its own as a separate embodiment.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US2014/034314 | 4/16/2014 | WO | 00 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2015/160340 | 10/22/2015 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
4709357 | Maki, Jr. | Nov 1987 | A |
5216638 | Wright | Jun 1993 | A |
5970798 | Gleman | Oct 1999 | A |
5995447 | Mandal | Nov 1999 | A |
6470749 | Han | Oct 2002 | B1 |
6745129 | Li et al. | Jun 2004 | B1 |
6856936 | Chen et al. | Feb 2005 | B1 |
7603238 | Valero et al. | Oct 2009 | B2 |
20030114989 | Blanch et al. | Jun 2003 | A1 |
20040100282 | Christensen | May 2004 | A1 |
20080210010 | Orth | Sep 2008 | A1 |
20130247673 | Davis | Sep 2013 | A1 |
20140107954 | Lih | Apr 2014 | A1 |
Number | Date | Country |
---|---|---|
1040251 | Jan 2012 | EP |
WO-0169286 | Sep 2001 | NO |
WO-2007050680 | May 2007 | WO |
WO-2013081608 | Jun 2013 | WO |
WO-2014011180 | Jan 2014 | WO |
WO-2015160340 | Oct 2015 | WO |
Entry |
---|
Lu, Yinghui, et al., “Numerical Implementation of Matching Pursuit for the Analysis of Complex Ultrasonic Signals”, IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 55(1), (2008), 173-182. |
Lu, Yinghui, et al., “Ultrasonic Signal Decompensation Via Matching Pursuit With an Adaptive and Interpolated Dictionary”, Review of Progress in Quantitative Nondestructive Evaluationvol. 26A, Thompson, D.O., et al., Editors, (2007), 579-586. |
“International Application Serial No. PCT/US2014/034314, International Search Report mailed Jan. 14, 2015”, 3 pgs. |
“International Application Serial No. PCT/US2014/034314, Written Opinion mailed Jan. 14, 2015”, 6 pgs. |
Number | Date | Country | |
---|---|---|---|
20150338378 A1 | Nov 2015 | US |