IPTC-10546-MS - A Novel Ultrasonic Cased-Hole Imager For Enhanced Cement Evaluation by Vankuijk 2005
IPTC-10546-MS - A Novel Ultrasonic Cased-Hole Imager For Enhanced Cement Evaluation by Vankuijk 2005
Flexural Mode
The calculated phase and group velocities for the flexural mode, where v f and v flex are the fluid velocity and the plate flexural
as a function of frequency, are shown in Fig. 1 for a water-loaded wave phase velocity, respectively. This angle is slightly
8-mm-thick steel plate. The phase velocity is far from constant, dependent on the fluid velocity and ranges from 28− 38 degs
meaning that the flexural wave shape (before taking the envelope) when the fluid velocity covers a typical range of 1,250 −1,650
changes rapidly during propagation along the casing. On the other m/s. The higher the slowness value, the larger is the phase-
hand, the flexural group velocity exhibits a peculiar feature: It matching angle. In a reciprocal way, a flexural wave propagating
goes through a maximum around 200 kHz and remains along the casing reradiates or leaks a wave into the materials
approximately constant over a fairly wide interval of 100 to 400 present on both sides. In particular, this allows the receiving
kHz. This means that all the frequency components of a transducer to pick up the signals. The radiated (or leaky) wave in
broadband wave packet within this frequency range will the fluid propagates at an angle again given by Snell’s law such
propagate at nearly the same velocity; hence, the wave packet that maximizing its reception by a receiving transducer requires
remains temporally compact. This fact has an important the latter to be aligned at the same angle that maximizes the
implication: Using the envelope of the flexural wave packet, it is excitation of this wave by the transmitting transducer.
possible to make accurate time of arrival and amplitude
measurements.
Flexural Mode Attenuation
The radiation into the inner and outer media, with respect to the
casing, is accompanied by attenuation of the flexural wave
amplitude along its path. This attenuation is sensitive to the
mechanical properties of both materials. To a good
approximation, the overall attenuation is the sum of the
attenuations due to the inner fluid and the attenuation due to the
IPTC 10546 3
outer material. The calculated attenuation due to water on one this effect decreases the flexural attenuation. The different cases
side of the casing (with vacuum on the other side to prevent any are depicted on Fig. 3 for fluid, slow cement, and fast cement.
leakage) is plotted in Fig. 2 as a function of frequency. For an 8-
mm-thick plate, this curve goes through a broad minimum of 0.31
dB/cm over the frequency range of flat group velocity (100−400
kHz). The fact that the attenuation is approximately constant with
frequency makes estimating it possible by measuring the decay
rate of the envelope of the received broadband pulse.
Fig. 3—Radiation of the flexural wave into water, slow cement ( v p <
Third-Interface Reflections
formation wall are parallel. However, in other directions, (30 mm). The vertical resolution of the attenuation measurement
they are not parallel, and the TIE is not reflected in the is controlled by the inter-receiver spacing (102 mm), while the
direction of incidence, leading to an amplitude drop (Fig. resolution of TIE-related measurements is controlled by the
6). transmitter-to-receiver spacing, at about 150 mm. The azimuthal
sampling is fixed at 10 degs and the vertical sampling can be
• Attenuation in the annulus material. Although water or adjusted, typically between 75 and 150 mm.
uncontaminated set cements have low attenuation, heavy
muds or contaminated/unset/foamed cement may have Because fluid properties can be estimated directly from the two
high attenuation. sets of measurements, no dedicated hardware for separate fluid
properties measurement is required.
• Acoustic contrast between annulus material and
formation. If the contrast is low, then the TIE amplitude
is low. This is true for example between lightweight
cement and shales. On the contrary, in double strings, the
contrast is very high, leading to rather strong TIE.
Data Processing
• Solid is defined through the expected type of cement. 2 and 10 m. The total width of these areas is also evaluated and
Using a laboratory-measured database, this material displayed as a single curve, which indicates the channeling
selection is converted into acoustic properties according severity.
to Table 1, and provisions are made for some
contamination or cement not fully set. Because the SLG map pertains to material immediately behind
casing, the second objective of the processing is to extract
relevant information from the TIE and quantify the full annulus
between the casing and the formation. First, the TIEs following
Cement Density P Velocity Z (MRayls) casing arrival are identified on the waveform envelope, and their
(kg/m3) (m/s)
arrival time and amplitude are determined. From the time
Class G 1,800 3,000 5.4 difference between TIE and casing arrival, and provided enough
azimuthal coverage is available, it is possible to extract the casing
Light Cement 1,200 2,800 3.4 centering in a dimensionless unit, where 100% represents perfect
Ultralight Cement 900 2,800 2.5
centering, and 0% represents fully eccentered casing (contacting
the formation). If the borehole diameter is known, then the time
Table 1—Acoustic properties of uncontaminated set cements. difference between casing arrival and TIE arrival can be further
converted into material velocity and displayed in map form.
The next step of the processing is to predict the measurements
from the expected acoustic material properties, which is trivial for Log examples
the acoustic impedance, but requires using a model for the
The first example (Fig. 9) is from a test well with 7-in., 23-lbm/ft
flexural attenuation. Then, multiple realizations of the
casing inside of a 9.625-in. casing and Class G cement. The
measurements noise are added to generate three clouds of points
section shown straddles the cement top. The left three tracks
(solid, liquid, and gas) in the two-dimensional measurements
display the SLG map and the two data sets used to compute the
plane. From these clouds of points, it is possible to define an
map, the flexural attenuation, and the pulse-echo cement
occurrence probability for each state (solid, liquid, gas), and the
impedance maps. These two maps are strikingly different; i.e., the
measurement plane can be mapped into different regions (Fig. 8)
flexural attenuation identifies two distinct, well-defined materials
by color. The white area corresponds to inconsistencies between
with high and low impedance. The impedance map, somewhat
the measurements, and may occur for example at collar locations.
affected by formation reflections from the outer casing, is barely
able to distinguish such regions. The interpretation is that higher
impedance cement partly displaced lower impedance cement
between the 330- and 290-m depths. The higher impedance
material is above the critical impedance and has a low
attenuation, actually lower than water, as seen from above the
cement top. The lower impedance material, either by design or
due to contamination, is below critical impedance and its high
attenuation leads to the high contrast. Both measurements are
combined into the SLG map, and a solid is shown where it
belongs, for both cement types, thus simplifying the display and
the low-level interpretation. Collars appear as white streaks
because the two measurements are disturbed enough to be
inconsistent. It can be noted that the flexural attenuation is much
less affected by formation reflection than the impedance map, as
long as the two casings are not touching. Contact actually occurs
at 255 m, and the attenuation shows some effect.
The four tracks on the right are related to the TIE information.
Track 4 is the map of TIE transit time difference with the casing
arrival. The white zones correspond to difficulty in the automatic
TIE picking due to the effects of the annulus material and or the
casing centering. This track is used to identify the sections with
high eccentering and the location of the near side, which is
consistent with the “galaxy” patterns seen on the impedance map.
Fig. 8—SLG mapping of the measurement plane for a Class G The transit time data are used to compute the casing centering
cement. Z is the impedance estimated from the pulse-echo technique, shown in track 7. This centering varies from ~100% (perfect
whereas the attenuation pertains to the flexural wave.
centering) to 0% (contact) within a short distance. Track 5 is the
TIE amplitude map and is affected by casing centering (which
The channel map is a further simplification of the measured affects the reflection geometry at the annulus-to-formation
SLG map focusing on the potential cementing problems. It shows interface and the propagation length inside the annulus) such as in
only the continuous uncemented areas with significant vertical the free-pipe section but also by the material type as seen in the
extent. This vertical extent can be adjusted, for example, between
IPTC 10546 7
two cements region. Casing collars from the outer casing string (tracks 5 to 8) answer the question and illustrate the potential
can also be located, for example at 315 and 325 m. Track 6 is the benefits in complex situations. The annulus velocity data show
computed annulus velocity map that is obtained from the TIE that the section xx68–xx73 m has a velocity ranging from 1,800
time map and the knowledge of the outer casing inner diameter. m/s to 1,900 m/s, clearly outside the range of plausible liquid, and
From the top, a value of 1,500 m/s, characteristic of water, is corroborating the SLG map. Furthermore, this velocity map and
shown, gradually increasing until 300 m where it reaches 2,000 to the amplitude map indicate that this solid material is not
2,500 m/s, consistent with the compressional velocity of a slow homogeneous but is three layers of slightly different materials. At
cement. Below 300 m, the TIEs are mostly propagated as shear, the top (xx30−xx34 m) of the liquid section, the velocity again
and the measured velocity drops back to ~1,500 m/s, consistent reaches 1,900 m/s, indicating the presence of a solid, although
with the shear velocity of a fast cement. even the flexural attenuation barely shows a slight increase above
its fluid values. This solid explains the nonflat upper interface,
The free-pipe section immediately above the cemented interval and since it floats on top of water, its density is probably less than
of the same well is shown in Fig. 10. The white line in the middle 1.0 g/cm3 and the targeted cement density. A small amount of
of the SLG map track is due to a casing groove caused by cable segregation within the cement slurry is a likely explanation for
motion, while the red patches along the same azimuths are due to such solid. The casing centering (track 8) is perfect near the
contact between the two casings. These data can be verified on casing centralizer, but 10 m below that depth quickly drops to
the right tracks that display the flexural waveforms along two about 50 %.
orthogonal diameters. The first echoes from the center are the 7-
in. inner casing arrivals, while the following ones are TIE Fig. 14 is a section from a well with 9.625-in., 47-lbm/ft casing
reflections from the outer 9.625-in. casing. The faint echoes after inside a 12.25-in. borehole. The mud density is 1.30 g/cm3, and
the first TIE are multiple reflections (up to five) within the the cement is Class G with a slurry density of 1.90 g/cm3. The
annulus. Such a display provides a geometrical sketch of the SLG map (track 1) indicates a uniform cement across the zone.
casing strings where contact points or potential geometrical The acoustic impedance ranges from 5 to 7 MRayls with some
defects can be identified. indication of formation reflections at about x810 m. The flexural
attenuation (track 2) exhibits both high (0.9 dB/cm) and low
Fig. 11 is an example of a 9.625-in., 47-lbm/ft single-string values (0.5 dB/cm) with the low values matching the high-
casing in a vertical well cemented with a high-performance, impedance area. The cement impedance value is close to the
lightweight cement (1.08 g/cm3). The mud density is 1.3 g/cm3. critical impedance value where the attenuation peaks. The TIE
The impedance image (track 5) is affected by the low impedance time map (track 4) confirms that the “galaxy” pattern shown at
of the cement and by casing roughness such that locating a depth X810 m is indeed due to the closer proximity of the
channel is made virtually impossible. This situation can be borehole wall. The higher impedance at x725 m and 180 degs
contrasted with the flexural attenuation (track 4) map that shows azimuth on track 3 can also be related to a formation reflection
a channel between x465 and x480 m with excellent contrast. The effect. The casing centering curve (track 8) is 100 % near the
processed SLG map (track 3) maintains the same quality with a centralizers, but eventually drops to below 50 % in the areas
clear liquid channel embedded within the cement. The SLG map between them. The small oscillations on the time map and the
has been further processed by the hydraulic communication centering curve are artifacts created by a corkscrew shape of the
algorithm to produce the channel map (track 6) and the channel borehole with a period of about 2 m. The annulus velocity map
width curve (track 7). This processing cleans the SLG map by (track 7) displays a rather uniform value around 1,750 m/s,
discarding small liquid patches and only keeping the liquid- characteristic of the shear velocity of class G cement. One
connected channels with a significant vertical extend, 2 m in this possible exception is the blue stripe at x775 m. Rather than a
case. A polar plot of the flexural waveforms in a variable density lower cement velocity, this is more likely due to a slight borehole
fashion displays the geometry of the casing within the borehole. enlargement. Assuming a constant cement velocity of 1,750 m/s
An inspection of the curvature of the TIE detected within the would give the possibility to present a cement sheath thickness
channel (Fig. 12) reveals that despite a casing centralizer at x474 over the interval. The TIE amplitude map (track 6) exhibits a
m the casing is slightly eccentered in the borehole and the channel striking correlation with the gamma ray (GR) (track 5) low
is located on the narrow side, contrary to first intuition. The amplitudes in shaly zones, as can be expected from a lower
absence of TIE across the cement azimuth may be due to a low impedance contrast between the cement and shales. Apart from
acoustic contrast between cement and formation. The CBL curve casing collars and centralizers, this map is also affected by casing
(Fig. 11, track 1) shows a fairly high reading due to the low centering that creates the two dark vertical stripes seen at about
cement impedance and a water-filled microanulus. x810 m and above x730 m.
unprecedented way, the flexural mode enables deep imaging of High-performance Lightweight Cement Slurries for
the cement sheath up to the cement-formation interface. Improved Zonal Isolation in Challenging Situations,”
The potential this tool holds has been demonstrated in various paper 47830 presented at the 1998 SPE/IADC meeting,
field tests. Enhanced contrast between light-weight cement and September.
4. Al-Suwaidi, A. et al.: "Light as a Feather, Hard as a
the displaced mud significantly improves data evaluation. When Rock," Oilfield Review (2001), 13, No. 2, 2-15.
reflections from the formation are detected, completely new 5. S. Zeroug and B. Froelich: “Ultrasonic Leaky-Lamb
measurements are possible, including the casing position within Wave Imaging through a Highly Contrasting Layer,”
the borehole, estimation of the cement velocity or the borehole Proc. 2003 IEEE Ultrason. Symp., November, 794-798.
shape, or imaging the formation behind casing. 6. S. Zeroug: “Forward Modeling for Ultrasonic Leaky-
Lamb Wave based Imaging through a Highly Contrasting
References Steel Cylindrical Layer,” Proc. 2004 IEEE Ultrason.
1. Hayman, A.J., Hutin, R., and Wright, P.V.: ”High Symp., August.
Resolution Cementation and Corrosion Imaging By
Ultrasound,” paper KK presented at the 1991 SPWLA
32nd Annual Logging Symposium, June 16-19.
2. Moulin, E., Revil, P., and Jain, B.: “Using Concrete
Technology to Improve the Performance of Lightweight
Cements,” paper 39276 presented at the 1997 SPE/IADC
meeting, November.
3. Revil, P. and Jain, B.: “A New Approach to Designing
IPTC 10546 9
Fig. 9—Top section of the cement in a test well with 7-in. casing inside of 9.625-in.
10 IPTC 10546
Fig. 10—Free-pipe section in a test well with a 7-in. casing inside a 9.625-in. casing
IPTC 10546 11
Fig. 12—Polar plot across channel from the far receiver at a depth of 477 m.
IPTC 10546 13
Fig. 13—A 9.625-in double casing section above the high-performance, light-weight cement
14 IPTC 10546
Fig.14 9.625-in. casing within a 12.25-in. borehole with Class G cement with indications of centralizer positions.