1. Introduction
In construction, poor-quality concrete can lead to a reduction in the load-carrying capacity of a structure. One common issue is the occurrence of cavities, which are internal voids that are not viewable from outside the structure. These defects can significantly impact the integrity of the concrete and weaken its properties. One way to detect these flaws is through nondestructive testing (NDT), a group of analysis techniques that can evaluate specific conditions and properties of concrete without damaging the specimen.
The American Concrete Institute report 12.2R3 provides valuable insights into the various NDT methods used for concrete analysis. Among these methods, using ultrasonic pulse-echo waves to detect flaws is a highly effective technique [
1]. With pulse-echo equipment, operators can quickly analyze a particular element on a single face without needing to access the concrete element's opposite side. The pulse-echo technique uses the synthetic aperture focusing technique (SAFT), a signal post-processing method that improves the accuracy of an area scan to detect concrete defects. This technique has proved to be highly reliable and has been employed in many studies to locate voids in concrete [
2,
3,
4].
The UPE testing principle involves a transmitter that converts electrical voltage energy into an ultrasonic wave and then travels through the concrete at a velocity determined by its properties. The wave continues to travel until it hits a concrete void or defect, reflecting the signal. The reflected signal then travels back through the material to the receiver, where it is converted into electrical energy, amplified as an echo, and recorded as an A-scan [
5,
6,
7].
Figure 1 shows the schematic of a typical ultrasound reflection principle in a UPE device, with a transmitting and a receiving transducer. Equation 1 allows an accurate calculation of the depth of the reflecting interface of the defect,
Z.
Cs refers to the material's shear-wave speed,
Δt is the travel time measured, and X is the horizontal spacing between the transmitter and the defect location (Equation 1). To understand the physical properties of material better, it is essential to consider its density (
φ), Poisson's ratio (
µ), modulus of elasticity (
E), and shear modulus (
G). By analyzing these factors, we can gain valuable insights into the behavior and characteristics of the material (Equation 2).
If applied effectively, UPE devices can be employed in the construction industry to inspect and evaluate the quality of reinforced concrete structures. Commercially available UPE devices, such as MIRA, EyeCon, and Pundit, utilize the tomography technique to generate images of the interior concrete, which can help identify potential imperfections. In this study, the dry-point-contact transducer array unit was used to enhance the accuracy of the inspection further. The UPE device used in this study is a low-frequency ultrasonic shear-wave tomography device with a surface longitudinal spacing of 30 mm, making it an excellent tool for detecting even the smallest defects in concrete structures.
Conventional software of UPE devices has limited the application to image concrete structures because of the assumption of a single-layer homogeneous, isotropic medium for concrete. A concrete bridge deck is a complex structure with multiple layers of different materials, each with varying shear-wave velocities. Due to this variation, the software used for inspection may provide inaccurate information about the multi-layered bridge deck, including depth, thickness, and internal defects. Furthermore, the software's signal processing may not accurately present reflected waves at the surface and defects. The accurate inhomogeneous materials may result in the generation of nonlinear features in wave propagation. The irregular surface texture may cause sensor contact variation and create a stochastic signal.
Pulse-echo testing is a non-destructive technique that uses ultrasonic waves to detect material defects. The technique works by analyzing the amplitude of the echo signal and the time it takes for the signal to reach the receiver. This information helps determine the defect's presence, size, and location. One of the primary advantages of the pulse-echo technique is its flexibility in testing large and irregularly shaped objects. However, the technique has a significant drawback: the loss of sensitivity near the test surface due to the coupling of the transducer with the test specimen. The ultrasonic signal passes through several materials before reaching the test specimen, including a coupling element and a transducer body. The reflected signals create near-field noise in the A-scan, which shows a reflector at each material interface.
The amplitude of the received echo depends on several factors, such as the transmitter power, direction of the transmission, size of the reflector, surface irregularities of the reflector, the reflector's position and orientation, the receiver's size and orientation, loss of signal at the receiver due to re-reflection and lack of coupling, attenuation of the sound wave due to absorption and scattering, and shadow effects, to name a few [
8,
9,
10,
11,
12,
13].
Ultrasonic testing is an effective method for detecting internal defects in concrete structures. However, interpreting ultrasonic testing data is a crucial task requiring extensive expertise. Various methods have been developed to overcome this challenge, including the widely used mature industrial ultrasonic imaging method, SAFT. These methods have shown promising results in detecting internal defects of concrete structures. However, it is essential to note that most of these methods can only qualitatively determine the presence of defects in concrete and cannot quantitatively detect the position and size of these defects.
Accurate positioning, precise evaluation, and clear visualization of void defects are essential for assessing the safety and performance of concrete structures. Unfortunately, the conventional SAFT method is often plagued by background noise and image artifacts, primarily due to the low-frequency ultrasonic pulse with a long wavelength. Consequently, the resulting image tends to droop, resulting in suboptimal quality outcomes.
Ultrasonic array devices are now widely employed to visualize the insides of concrete structures non-destructively. However, the data collected by these devices may sometimes need to be clarified, requiring the use of image reconstruction algorithms to achieve clear images. Low-frequency UPE devices can combat the issue of signal attenuation in concrete structures. These devices emit low-frequency ultrasonic pulses that reduce signal attenuation, pulse duration, and image sagging, resulting in fewer image distortions and an overall improvement in quality.
The UPE device has an impressive collection of ultrasonic transducers that can measure multiple pulse-echo signals in a single scan. This advanced process results in a comprehensive cross-sectional image of the tested object. The device runs on the cutting-edge SAFT-C algorithm, which uses a time-domain approach that relies on the delay-and-sum method to focus on delayed reflections effectively. This algorithm was selected for its user-friendliness and low-performance requirements, making it an excellent microprocessor option.
While the interpolated image provides an intuitive and immediate visualization of the inside of concrete structures, the reliability of the information in the interpolated area decreases as the spacing between 2D images increases. Using a polarized shear wave by the UPE device makes it highly unlikely for a reflector to be arranged in a specific direction. Additionally, the UPE device software cannot combine data collected from multiple orientations to generate an image.
These UPEs require greater precision in quantitatively measuring internal defects within concrete structures. Additionally, detecting beyond a reflector area is challenging due to the high reflectivity of the shear wave. The UPE devices may not be sensitive enough to vertical crack-type defects because of the shallow frequency, resulting in a longer wavelength. Varying ultrasonic wave propagation characteristics like velocity, amplitude, attenuation, and frequency are utilized to measure the thickness and detect defects in concrete structures.
Commonly used methods for detecting voids, imperfections, and defects in concrete structures utilize ultrasonic body waves. Concrete is a homogeneous medium with a consistent shear-wave velocity, and the SAFT technique has been implemented in the UPE devices to achieve this. However, certain modifications can be introduced to improve the results obtained through SAFT, such as calibrating the sagging in the resulting images due to the long wavelength of the pulse (as suggested in references [
14,
15,
16,
17,
18,
19]). Such enhancements can significantly improve the efficacy of SAFT in detecting and characterizing concrete anomalies, thereby increasing the safety and durability of concrete structures.
Two-time indices are called
t1 and
t2 in the ultrasonic transmission, reflection, and receiving process. By precisely analyzing the designated timings, we can effectively discern reflections originating from the front surface and any imperfections (
Figure 1). This pulse type is instrumental in generating an A-scan that meets rigorous quality standards. The data derived from this scan is the time differential tf, which is represented as
tf =
t2 -
t1 and commonly known as the ultrasonic time-of-flight (TOF). Knowing the propagation speed of bulk waves in the material, the TOF is invaluable in determining the location of the defect.
TOF-based methods have proven effective in detecting, locating, and sizing faults in ultrasonic nondestructive testing and evaluation. However, echo signals in UPC testing often contain significant noise. The intensity of this noise can vary depending on the medium's properties and the distance of ultrasonic wave transmission, resulting in low signal-to-noise ratios and peak variations in the time domain. These factors pose a challenge in accurately locating and extracting peak amplitudes from defect echo signals.
Various signal processing techniques enhance the detection capabilities of ultrasonic NDT applications, where signals are often mixed with noise. The ultrasonic signals should be denoised, which is feasible using short-time Fourier transform, wavelet transform, improved wavelet transforms, adaptive filtering, and empirical mode decomposition (EMD) [
20,
21,
22,
23,
24,
25,
26]. Each technique has its strengths and weaknesses, which are considered when selecting the appropriate method for the given application.
Originating from concrete defect echoes in UPC testing; the resulting intense noise makes identifying and quantifying defect echoes difficult. Therefore, the conventional averaging and filtering techniques could be more helpful for pulse-echo noise reduction. Denoising assists with filtering the noise components and retaining the original signal's details.
The most used denoising algorithm thus far is the wavelet threshold algorithm. However, the wavelet threshold denoising has a severe drawback. It requires appropriate wavelet base, threshold, and decomposition level values, thus allowing different parameters to affect the denoising performance remarkably. Empirical Mode Decomposition (EMD) works for nonlinear, stochastic, and non-stationary signal processing. The EMD denoising algorithm does not require presetting the base and decomposition level and is more adaptable than the wavelet threshold algorithm. Furthermore, the EMD denoising algorithm can result in high-frequency resolution.
The electrical noise and piezoelectric signal primarily determine the position of the initial pulses in the UPE device. The signal can travel along a straight line when the distance is short and the temperature variation is slight. The TOF of the echo wave determines the position of the scatter.
The Hilbert transform has been used in signal processing to map an actual signal into an analytical signal with a complex envelope to obtain specific signal features [
27,
28]. The Hilbert transform converts a given signal into an analytical signal with a complex envelope, which facilitates the evaluation of the signal envelope and the determination of the TOF of the defect echo signal. This transformation is beneficial when the magnitude of the defect echo signal is relatively small and gets merged with a superimposed signal.
This paper proposes an approach to denoise ultrasonic pulse-echo signals using advanced EMD and Hilbert transform to identify and quantify the defect echo signal. This advancement complements the commercial UPE device results and accurately calibrates defect locations. This study's experimental program shows testing to validate the proposed denoizing approach.
3. Formulation for Accurate Assessment of Defects and Voids
Most ultrasonic measurements exhibit perplexity, as the results are subject to severe disturbances such as signal deviation, nonlinearity, stochastic surface status, and mode conversion due to anisotropic and heterogeneous domains in the concrete. These are reflected in the UPE signal as spurious echoes and noise.
The UPE testing data analysis contains ambient noise or signals corrupted by non-stationary acoustic noises, reducing the analysis efficiency. Denoising of the signals appears in several studies reported in the literature, many of which have certain limitations. Most studies relied on prior signal knowledge to enhance denoising, resulting in the loss of important data during filtering.
Typical signal processing methods may not be effective in detecting fault features. However, combining various signal processing techniques can improve fault detection and analysis. To more accurately analyze the behavior of non-linear and non-stationary signals, many researchers have utilized Empirical Mode Decomposition (EMD) to identify faults. The EMD is capable of removing signals from non-linear and non-stationary signals.
Research conducted [
29,
30] highlighted that the EMD denoising method surpasses both median filtering and wavelet denoising in effectiveness. Nonetheless, when heavy noise is a factor, selecting the intrinsic mode function (IMF) in EMD can be difficult. The Hurst exponent is a valuable resource for detecting multifractality that may be disguised in nonlinear and non-stationary signals. The Hurst analysis technique was employed to choose the appropriate IMF with the EMD method [
31,
32].
This study presents a novel approach to removing noise from a UPE signal using data-driven analysis. The method involves breaking down the noisy signal into IMF components using EMD analysis, and then examining each IMF's Hurst exponent to determine which components require filtering. The filtered components are then reconstructed to eliminate the noise from the data.
By utilizing the Hurst exponent to identify the threshold for maximum signal noise suppression, the method applies thresholding techniques to the IMFs for optimal noise removal. This method has proven more effective than traditional techniques, especially for low SN signals.
A comprehensive understanding of IMF statistics in noise-only situations is essential in identifying the significance of a given mode. The estimation of the noise Hurst exponent is achieved using the aggregated variance method. This exponent is then utilized to evaluate the energies of IMFs within the noise model, a crucial step in determining the threshold of IMFs. Understanding the statistics of IMFs in noise-only scenarios can help determine the importance of a specific mode. Thus, the aggregated variance method is employed to estimate the noise Hurst exponent, which is then applied to evaluate the IMF energies and establish the IMF threshold.
Ultrasonic nondestructive testing can present challenges when it comes to detecting echoes caused by defects in backscattered signals due to backscattering and electronic noise. However, there is a solution in the form of the EMD algorithm. The EMD algorithm breaks down a non-linear and non-stationary signal into a series of zero-mean amplitude-modulation and frequency modulation components. By doing so, it can accurately represent the observation's characteristic time scale. This means that a multi-component, non-linear, and non-stationary signal can be represented with precision using the EMD within the context of ultrasonic theory.
To decompose a signal using the EMD approach, we represent the instantaneous amplitude and instantaneous phase of the jth component as and , respectively, and n as the number of components. This is done through an iterative sifting process, resulting in zero-mean AM-FM components called IMFs. These IMFs must satisfy two requirements: (a) The number of extremes and zero crossings in the IMF must be equal or differ at most by one, and (b) The mean value of the envelopes defined by the local maxima and local minima must be zero at any point. In essence, the signal must be locally symmetrical around the time axis. To find the IMFs for a given signal, we conduct the sifting process to find the IMFs for the signal , which involves the following steps:
- (1)
To accurately process the input signal, we must initially identify all local maxima and minima, along with their respective positions and amplitudes. Following this, we can employ cubic spline interpolation to generate an upper envelope consisting of the local maxima and a lower envelope using the local minima.
- (2)
After calculating these envelopes, the envelope mean signal, known as , can be determined by taking their mean. Lastly, to finalize the processing, we need to subtract the envelope mean signal from the original input signal (Equation 4).
- (3)
Check if meets the IMF requirements. Treat the data as new data and repeat the process if it does not meet the IMF requirements (Equation 5).
- (4)
Repeat the sifting procedure k times until the resulting component
is an IMF, which becomes the first IMF (Equation 6).
- (5)
A standard method for extracting and analyzing a signal's underlying components is residual analysis. This process involves subtracting the
component from the input signal and defining the resulting remainder as the first residual. Given that the first residual
may contain information relating to longer-period components, it is treated as a new data stream. The procedure is repeated for this new signal. This process may be iterated j times, resulting in the generation of
j residuals. By following this approach, it is possible to obtain a refined understanding of the signal and identify the underlying components contributing to its overall structure.
- (6)
The sifting process is interrupted once either of the two criteria mentioned above is fulfilled: firstly, when the component
or the residual
is reduced to such a minuscule size that it can be regarded as insignificant, or secondly, when the residual (R) becomes a monotonic function that precludes the extraction of the IMF—the objective IMF can be obtained by adding Equations (4) and (5). The original signal can be expressed as a combination of IMFs and a residual, producing significant implications in signal processing and analysis. This observation has been accepted in the academic community.
- (7)
EMD-based denoising, similar to other decomposition-based denoising techniques like wavelet transforms, requires a reliable and robust threshold to distinguish between noise and authentic signal components. In cases where irregularity or noise is present in a time series, the Hurst exponent plays a crucial role in determining the irregularity in the signal. This methodology is more efficient than traditional approaches such as autocorrelation, ANOVA, and spectral analysis in many applications. The Hurst exponent value, whether greater or less than 0.5, indicates the pattern of the nonlinearity of the data set. Some white noise signals have a flat spectrum and are determined by the Hurst exponent H. The autocorrelation function for a zero-mean Gaussian stationary process is expressed as:
- (8)
Equation 10 represents the interdependence of process variance (σ), Hurst exponent (H), and correlation lag (k). Notably, when H equals 0.5, the process is classified as uncorrelated white noise, whereas for other H values, it is labeled as colored Gaussian noise. Moreover, when a generalized white noise signal is subjected to EMD, it acts as a dyadic filter bank. It is important to note that the log-variance of the IMFs follows a simple linear model, which the Hurst exponent of the process ultimately governs.
- (9)
The energy of each of the IMFs, for k≥ 2 and ρH ≈ 2, can be parameterized as a function of the first IMF energy (Equation 11). The energy of the first IMF is given by Equation 12.
- (10)
This particular model can execute denoising BMD-based techniques. The process entails breaking down the noisy signal into IMFs and gauging their energy levels about the estimated noise-only IMF energies derived from Equation10. From there, the signal reconstruction is accomplished by adding up the IMFs whose energy levels deviate from the expected noise model.
- (11)
Peak detection techniques are typically employed to estimate TOF, thereby differentiating between the reflection signal from the front surface and the reflection signal from the defect. Despite the denoising process, the defect echo signals may still exhibit dispersion and weakness, requiring specialized methods for identification and estimation. Techniques such as filtering, cross-correlation, envelop moment analysis, and matching pursuit decomposition with dispersion compensation are necessary for accurate defect detection in such scenarios. The envelope of an echo signal constitutes a vital characteristic that can be employed to extract information regarding the location of the echo waveform.
The Hilbert transform facilitates the computation of instantaneous features of a time series, including the envelope amplitude and instantaneous frequency. The instantaneous envelope denotes the amplitude of the complex Hilbert transform, which can be utilized to identify the dynamic characteristics of nonlinear systems by improving the accuracy of envelope detection through the use of local-maxima interpolation.
The Hilbert transform involves casting the actual signal
a(t) into an analytical signal (i.e., complex envelope)
ã(t). The ã(t) is a complex-valued time-domain signal, where the original actual signal a(t) is the genuine part, and the Hilbert transformed signal â(t) is the imaginary part (Equation 13). By leveraging this approach, one can obtain a complex-valued time-domain signal that provides a more comprehensive representation of the underlying physical phenomenon.
- (12)
In Equation 14, j is the imaginary number, while H[·] denotes the Hilbert transform operation. In the time domain, the Hilbert transform is defined as the convolution of a(t) with 1/πt, where â(t)= H[a(t)]. The envelope of a signal is the magnitude of the analytical signal, which is the same as the magnitude of the real signal. The complex signal ã(t) is formed by the Hilbert transform, a(t), and â(t) as shown in Equation 15. Then, the envelope of the real signal can be given by Equation 16.