Abstract

Guided wave technique is an efficient method for monitoring structural integrity by detecting and forecasting possible damages in distributed pipe networks. Efficient detection depends on appropriate selection of guided wave modes as well as signal processing techniques. Fourier analysis and wavelet analysis are two popular signal processing techniques that provide a flexible set of tools for solving various fundamental problems in science and engineering. In this paper, effective ways of using Fourier and Wavelet analyses on guided wave signals for detecting defects in steel pipes are discussed for different boundary conditions. This research investigates the effectiveness of Fourier transforms and Wavelet analysis in detecting defects in steel pipes. Cylindrical Guided waves are generated by piezo-electric transducers and propagated through the pipe wall boundaries in a pitch-catch system. Fourier transforms of received signals give information regarding the propagating guided wave modes which helps in detecting defects by selecting appropriate modes that are affected by the presence of defects. Continuous wavelet coefficients are found to be sensitive to defects. Several types of mother wavelet functions such as Daubechies, Symlet, and Meyer have been used for the continuous wavelet transform to investigate the most suitable wavelet function for defect detection. This research also investigates the effect of different boundary conditions on wavelet transforms for different mother wavelet functions.

1. Introduction

Early forecasting of the degradation process caused by adverse environmental effects or mechanical damages in pipe network systems can avoid many catastrophic accidents. Now a days, the detection of the existing defects in pipes is one of the major challenges for the structural health monitoring of pipes. Propagation of cylindrical guided waves through pipes for damage detection is becoming an increasingly popular technique for pipe inspection.

Gazis [1, 2] first analytically solved the propagation of harmonic waves in an infinitely long elastic hollow cylinder. Many investigators [3, 4] used guided wave modes for detecting wall thinning defects in cylindrical pipes. Rose and coworkers [5, 6] designed a special probe which can be used both as a transmitter and a receiver during the pipe inspection. Guo and Kundu [7, 8] designed a new transducer holder mechanism for pipe inspection using cylindrical guided waves. Using those transducer holders Na et al. [9, 10] generated cylindrical guided waves for detecting delaminations between steel bars and concrete interface. Ahmad and Kundu [11] investigated the influence of water flow on pipelines using guided wave techniques.

Successful damage detection requires not only the appropriate guided wave modes but also proper use of signal processing techniques or tools. In recent years, wavelet analysis has become a popular technique for processing received signals with time-varying spectra. Many investigators have used the wavelet analysis to characterize damages in materials. Cho et al. [12] discussed the detection of subsurface lateral defects using wavelet transform on propagating Lamb waves. Rioul and Vetterli [13] and Abbate et al. [14] used wavelet transform for processing signals with nonstationary spectral contents. Kaya et al. [15] used wavelet decomposition to detect flaws in stainless steel samples. Ahmad et al. [16, 17] used daubechies wavelet transforms for detecting defects in free and embedded pipes. Ahmad and Kundu [18] also used continuous wavelet transforms with different wave functions on cylindrical guided wave signals to identify defects in pipes under different boundary conditions.

Gabor transform can be used as another form of wavelet analysis. Gabor [19] adapted the Fourier transform to analyze only a small section of the signal at a time—a technique called windowing of the signal. Bastiaans [2022] related the Gabor expansion and the short-time Fourier transform. Bastiaans introduced the sampled short time Fourier transform to compute the Gabor coefficients and successfully derived a closed form Gaussian function.

Murase and Kawashima [23] tried out different wavelet transforms and showed that when Gabor functions are used as mother wavelets then one can plot group velocity curves for a thin aluminum plate. Ahmad and Kundu [24] also used the Gabor wavelet transform to plot group velocities for defective and defect-free cylindrical pipes from experimental data. Ahmad et al. [25] used daubechies mother wavelet functions to identify defects in transmission pipes.

2. Theory

A wavelet is described by the function , which is obtained by dilation and translation of a function as defined by: where is the set of real numbers. The function is called mother wavelet. The parameter represents the scale index, and indicates the time shifting or translation. The approximation of a function by wavelets is carried out using the coefficients given by: where defines the scalar product.

The continuous wavelet transform or CWT is given by: where * denotes the complex conjugate of the function.

In simple words, the continuous wavelet transform is defined as the sum over all time of the signal multiplied by scaled, shifted version of the wavelet function, : The results of the CWT are many wavelet coefficients , which are a function of scale and position. Therefore, wavelet analysis is nothing but the breaking up of a signal into shifted and scaled version of the original mother wavelet.

3. Experiments

The primary objective of this research is to investigate how cylindrical guided waves can be used effectively to detect defects in steel pipes under different boundary conditions. In the process, it will be explored how different signal processing techniques like Fourier transform and continuous wavelet analysis (CWA) influence the detection process. Fourier transforms will endow with possible propagating modes, and CWA provides scaled coefficients affected by the presence of defects. CWA was performed using Daubechies, Symlet, and Meyer mother wavelet functions on the signals received for both defect free and defective pipes. The pipes were kept under three different boundary conditions: traction-free boundary condition when pipes were kept in open air; when water flows through the pipes, keeping the inner boundary of the pipes in contact with water and outer boundary in open air; the pipes are embedded in soil keeping the outer surface in contact with soil and the inner surface in traction-free boundary condition. Figures 1(a) and 1(b) show the schematic diagrams for different boundary conditions and instrumental arrangements, respectively. Relatively high-frequency (1 MHz) piezoelectric transducers were used as transmitters as well as receivers. The transmitter was activated by a tone-burst excitation using WAVETEK function generator (Model 395) and a MATEC-gated 310 Amplifier for high-frequency excitation. The received signal was amplified by a MATEC 605 amplifier. The signal amplitude was then displayed on an oscilloscope screen as a function of time and frequency. To hold the transducers onto the test specimens, two transducer holders developed by Guo and Kundu [7, 8] were used. Harmonic 1 MHz piezoelectric transducers with a diameter of 13 mm (0.5 inch) were used for the experiments. Vaseline was used as the couplant.

3.1. Specimens

A specimen set, consisting of three steel pipes, were fabricated. All three pipes were 1200 mm (~4 ft) long and had 21.4 mm (13.5/16 inch) outer and 15.6 mm (10/16 inch) inner diameters. One pipe was defect-free, and the other two had mechanical defects: a gouge and a dent, that were artificially fabricated on the pipe. The gouge anomaly was fabricated by pressing the outer wall of the pipe while keeping the inner diameter unchanged by placing a rigid rod inside. The dent type anomaly was formed by pressing the outer wall, keeping the inner wall free to deform. In both types of defect, the outer walls were cold-pressed. The defect covered a complete (360°) revolution. Table 1 gives dimensions of the pipes and the defects. Table 2 shows the acoustic properties of steel—the pipe material. Properties of the surrounding soil used in this experiment are given in Table 3.

4. Results

Experiments were carried out in three phases for pipes with three boundary conditions. During the experiment, cylindrical guided waves were generated by piezoelectric transducers. Guided waves were generated at one end of the pipe and received at the other end by another receiving transducer in a pitch-catch configuration. The received signals were in the form of curves (i.e., amplitude versus frequency) and time series (amplitude versus time).

The effect of different signal processing techniques for detecting defects in pipes was then investigated. Two separate signal processing techniques were used:(a)the first segment used experimental curves to identify the propagating modes through the pipes;(b)the next segment applied continuous wavelet analysis (CWA) on experimental time series signals, using different types of mother wavelet functions to investigate which mother functions are more effective in identifying defects in pipes for all three boundary conditions.

Effects of these signal processing techniques are discussed in the following sections.

4.1. Experimental Curves and Identification of Propagating Modes

In this section, the sensitivity of pipe defects on the propagating guided wave modes is investigated under different boundary conditions. Generation of guided waves in pipes is very sensitive to the incident angle of the transducers. The first challenge is to find the appropriate incident angles for which strong guided waves can be generated. The incident angles of the transmitter were adjusted experimentally to obtain strong signals. The incident angles of 28° and 57° were found to produce strong guided wave signals. The 57° inclination angle produced the strongest signal and used in this paper to construct the curves.

4.1.1. Traction-Free Boundary Condition

Figure 2(a) shows the experimental curves for an incident angle of 57° when the pipes were kept in open air with traction-free boundary condition. As expected, defect-free pipe produced stronger signals than the defective pipes (Figure 2(a)). Theoretical phase velocity dispersion curves for a cylindrical pipe were plotted using the geometric and material properties of pipes given in Tables 2 and 3 using DISPERSE software. In order to identify the propagating guided wave modes through the defective and defect-free pipes, phase velocities are calculated from the peaks in the experimental curves using Snell’s law [ = /].

Figure 2(b) shows theoretical phase velocity dispersion curves for a cylindrical pipe (defect-free), where both longitudinal and flexural modes are plotted. For the defect free signal, phase velocities were calculated for peaks at 380, 1020, and 1140 kHz plotted in Figure 2(b). The phase velocity corresponding to the peak at 380 kHz matches closely with the theoretical mode. Peaks at 1020 and 1140 kHz match closely with , , and modes for the defect-free pipe. This means that the peaks at 1020 and 1140 kHz represent propagating modes , , or , or all of them together. For the defect-free pipe, the strongest peak occurs at 1140 kHz which means that the strongest propagating modes could be , or . For the gouged pipe, the peaks could be selected at 580, 700, and 1180 kHz (Figure 2(a)), and the corresponding guided wave modes can be identified as , , and (Figure 2(b)). The strong peak at 700 kHz for the gouged pipe corresponds to the propagating mode. For the dent pipe, the peaks are observed at 660 and 1300 kHz that again correspond to , , and modes. The stronger peak at 660 kHz represents mode for the dent pipe. From this investigation, it can be observed that for the traction-free boundary condition, stronger , and modes propagate for the defect-free pipe and for the gouged and dented pipes, strong mode and weak and modes may propagate.

4.1.2. Pipes with Water Flow

Figure 3(a) shows experimental curves obtained for both defect-free and defective pipes when water flows through the pipes when the incident angle is 57°. In this condition, the strength of the signal for the defect-free pipe is stronger than the defective pipes. It can be observed that the strength of the signals (Figure 3(a)) is weaker than that obtained for the traction-free condition (Figure 2(a)). Presence of flowing water reduces the strength of the signals implying that a part of the ultrasonic energy is leaked into the water and does not reach the receiving transducer. Theoretical phase velocity dispersion curves are plotted in Figure 3(b) for the condition when water flows through the pipes. A conceptual sink is placed at the center of the pipe for plotting phase velocity dispersion curves. Due to the presence of water, and are absent from the theoretical phase velocities. From the experimental curves shown in Figure 3(a), peaks are selected for calculating the phase velocities to identify the propagating modes. For the defect-free pipe, the selected peaks are at 580, 1100, and 1140 kHz. The peak at 580 kHz corresponds to mode. The stronger peaks at 1100 and 1140 kHz most likely correspond to and modes, respectively. For the gouged pipe, peaks which are weaker than the defect-free pipe are selected at 620 and 1220 kHz. The possible propagating modes can be or . For the dented pipe, the major peaks occur at 500 and 1180 kHz which corresponds to and/or mode. In summary, it can be stated that when water flows through the pipes strong and modes propagate through the wall of the defect-free pipe at slightly higher frequency (>1000 kHz). For the gouged pipe, the possible propagating modes are weak and . For the dented pipe, the possible propagating modes are and .

4.1.3. Pipes Embedded in Soil

Figure 4(a) shows the experimental curves, when the pipes are embedded in soil. The soil properties are given in Table 3. It can be observed that multiple strong and distinct peaks which referred to as propagating guided wave modes are generated for both the defective and defect free pipes. It is interesting to note that the strengths of the signals are similar for both defect-free and defective pipes. Due to the presence of surrounding soil, the overall strength of the signal is reduced compared to the signal strengths in cases of traction-free boundaries. The dented pipe produces stronger and distinct peaks at 620, 740, 1020, and 1340 kHz. The 620 and 740 kHz peaks match with the mode in Figure 4(b). The peaks at 1020 and 1340 kHz represent either , , or mode. These peaks are strong. In previous cases (Sections 4.1 and 4.2), the experimental curve for the dented pipe produced weaker peaks which means generation of weaker propagating modes. Clearly, the propagating modes decay before traveling long distances when the dented pipe is kept in the open air with water flow or no-flow condition. When the pipe is buried under soil, distinct propagating modes are generated that can travel a longer distance. For buried dented pipes, more modes are generated that are very sensitive to defects.

For the defect-free pipe, peaks are observed at 660, 860, 980, and 1140 kHz. Peaks at 660 and 860 kHz represent mode. Peaks at 980 and 1140 kHz represent possible , , or mode. For the gouged pipe, peaks are generated at 700 and 1340 kHz. These peaks correspond to , , or mode. It should be noted that for the defective pipes like dented and gouged, stronger peaks are obtained at a lower frequency (<1000 kHz) which refers to the mode. It implies that for pipes embedded in soil, strong guided wave modes can be generated at low frequencies which would be efficient in detecting defects.

4.2. Continuous Wavelet Transform (CWT) with Different Mother Functions

This section investigates the effectiveness of continuous wavelet transforms (CWTs) as a signal processing tool in assessing the integrity of pipes for the three stated boundary conditions. CWT has a broad spectrum and is associated with proper choice of mother wave function and appropriate scaling. Wavelet analysis produces a time-scale view of a signal that has a scale aspect and a time aspect. Scaling a wavelet means stretching it. The smaller the scale factor, the more compressed the wavelet; while higher scales correspond to the more stretched wavelets. In this section, CWTs are performed on experimental time-series signals obtained for the 3 sets of boundary conditions. To efficiently use CWT as a signal processing tool for health monitoring, it is important to select appropriate mother function with proper scaling which will facilitate the integrity assessment. In Section 4.2.1, the range of scaling on different mother functions is investigated in order to identify appropriate scale for different types of mother wavelet functions. The time series signals were associated with significant amount of white noise. Using wavelets, noise from the signals are removed by identifying which component or components contain the noise. The signals are then reconstructed omitting those components. The one dimensional model for a noisy signal is basically of the form: where time is equally spaced, is a Gaussian white noise , and noise level is equal to 1. The objective of de-noising is to suppress the noise part of the signal and to recover .

4.2.1. Scaling for Different Wave Functions

Fixing the proper scale for a mother wave function is particularly important for processing a signal. The effectiveness of a CWT depends on the right scaling, which stretches or contracts a wave function. For this investigation, signals received for the defect-free pipe in traction free boundary conditions were used. The wave functions that were used for this investigation are Daubechies (db4), Symlet (sym4), Coiflet (coif2), Gaussian (gauss1), Meyer, and Mexican Hat with scaling of 16, 64 and 256.

Figures 5(a)5(f) show CWT of the signal for the defect-free pipe in traction-free boundary condition for a scaling of 16. The experimental signals which are associated with significant white noise were denoised before applying CWT. Figures 5(a), 5(b), and 5(c) that used db4, sym4, and coif2 mother functions, respectively produced similar plots which are helpful in identifying the difference between defective and defect-free pipes. For gauss1 (Figure 5(d)) and Mexican Hat (Figure 5(f)), the difference is not clearly visible for a scaling of 16. For Meyer wavelet function, a scale of 16 gives clear distinction between the defect-free and defective pipes. For a scaling of 64 (Figures 6(a)6(f)), the wave functions, db4, sym4, coif2, gauss1, Meyer, and Mexican Hat do not produce distinct signals which can clearly show the difference. Similarly, for a higher scale of 256 (Figures 7(d), 7(e), and 6(f)), gauss1, Meyer and Mexican Hat do not exhibit distinct differences. The plots produced by db4, sym4 and coif2 (Figures 7(a), 7(b), and 7(c)) wave functions for a scale of 256 are stretched too long and the difference between the defect-free and defective signals are not so differentiable. From this investigation, it can be concluded that for a lower scale of 16, the CWTs are more effective in assessing the integrity of pipes. It can also be observed that wave functions gauss1 and Mexican Hat do not provide any significant information about the differences between defective and defect free pipes. It can be said from this investigation that db4, sym4, coif2, and Meyer wave functions should produce more meaningful results. It may be noted that sym4 and coif2 functions produce almost similar results as well as plots. In the following, only db4, sym4, and Meyer wavelet functions with a scaling factor of 16 are used to investigate the integrity of pipes.

4.2.2. Continuous Wavelet Analysis for Defect Detection of Pipes under Different Boundary Conditions

The effectiveness of continuous wavelet analysis (CWA) in detecting defects in pipes under different boundary conditions mentioned in the previous section (Section 4.1) is then investigated. Daubechies (db4), Symlet (sym4), and Meyer wave functions were employed with a scaling of 16 for this investigation.

Continuous Wavelet Transform for Traction-Free Boundary Conditions
Figures 8(b), 8(c), and 8(d) show continuous wavelet analyses for defect-free and defective pipes under traction-free boundary conditions using db4, sym4, and Meyer functions, respectively. The signals are associated with significant amount of noise. Using db4 wavelets, noise from the signals is removed by identifying which component or components contain the noise. The signals are then reconstructed omitting those noisy components. Figure 8(a) shows the denoised time series signal for the defect-free, gouged, and dented pipes under traction-free boundary conditions. Figures 8(b), 8(c), and 8(d) show the calculated wavelet coefficients with respect to time for db4, sym4, and Meyer functions, respectively. The time scale is shown in 1200 equal divisions, where each division accounts for 7.2 ns (nano seconds). The time series signal (Figure 8(a)) shows that the first patch of modes appear between 400 and 600 time divisions. When db4 and sym4 wave functions are used (Figures 8(b) and 8(c)) some difference in the wavelet coefficients can be observed between the defect-free and the gouged pipe, but it is not so pronounced. But there is a significant difference in wavelet coefficients for the dented pipe. In Figure 8(d), where Meyer wave function is used for the CWT, significant difference in wavelet coefficients can be seen even between the defect-free and gouged pipes. For the dented pipe where the Meyer function with a scaling of 16 produces prominent difference between the defect-free pipe and the dented pipe also. The calculated coefficients for the defect-free pipe are much higher than those for gouged and dented pipe. Therefore, Meyer function gives a clear difference between defect-free and defective pipes. In other words, Meyer wavelet functions are very sensitive to the defects and produce significant difference between the calculated CWT coefficients for defective and defect-free pipes. It makes Meyer wavelet better suited for structural health monitoring.

Continuous Wavelet Transform When Water Flows through the Pipes
Figure 9(a) shows the denoised signals when the water flows through the pipes. Figures 9(b), 9(c), and 9(d) show continuous wavelet analyses for defect-free and defective pipes using db4, sym4, and Meyer functions, respectively. As expected, the signal strength for the defect-free pipe is stronger than that of the defective pipes (Figure 9(a)). The difference in signal strength is not apparently visible in Figure 9(a). When CWTs were performed on these signals, the difference can be clearly observed. In Figures 9(b) and 9(c), where db4 and sym4 functions are used respectively for the CWT, the difference in wavelet coefficients for the defect-free and defective pipes can be observed. Among the wave functions that are used for CWT, Meyer function provides best option for finding the difference between the defect-free and the defective pipes (Figure 9(d)). The difference in wavelet coefficients for the defect-free and defective pipes, obtained using Meyer function, is significant and can be easily recognized.

Continuous Wavelet Transform When Pipes Are Embedded in Soil
The continuous wavelet transforms for the defect-free and defective pipes embedded in soil are shown in Figures 10(b), 10(c), and 10(d) using db4, sym4, and Meyer functions, respectively. In Figures 10(b) and 10(c), where Daubechies (db4) and Symlet (sym4) functions are used respectively, the wavelet coefficients for the defect-free pipe differs considerably from the defective pipes (gouged and dented) for the first patch of propagating modes. It can be also noted that the patterns for the defective pipes (gouged and dented) are similar, and they differ from the pattern shown by the defect-free pipe. The pattern also differs significantly from the traction-free boundary conditions. The signal strength for the defect-free and the gouged pipe is similar, but the strength for the dented pipe is considerably reduced due to the presence of soil (Figure 10(a)). Wavelet coefficients for the defect-free and the gouged pipe are also similar when db4 and sym4 functions are used (Figures 10(b) and 10(c)). Significant difference in wavelet coefficients can be observed when Meyer function is used (Figure 10(d)). This phenomenon is also observed for the traction-free boundary conditions. Even for the soil embedded condition, where the acoustic signal is expected to loose significant amount of energy, CWT using Meyer’s function provides strong evidence of identifying defects in embedded pipes and can easily be distinguished from the defect-free pipe.

5. Conclusion

This paper investigates the importance of signal processing techniques in guided wave application for detecting defects in pipes. Two different approaches were investigated and compared. Fourier transforms of experimental signals provide information regarding the propagating modes and help in identifying modes that are sensitive to defects. Continuous wavelet analysis is helpful in identifying defective pipes by comparing wavelet coefficients for pipes under different boundary conditions. This research also investigates which mother wave functions are more effective in distinguishing defective pipes from a defect-free pipe. Depending on the application, any or both approaches could be utilized for successful monitoring of structural integrity of pipes.

Acknowledgment

This paper was financially supported from the National Science Foundation, Grant no. CMS-9901221 and CMMI-0530991, Dr. S. C. Liu (program manager).