Next Article in Journal
A Novel Workflow for Electrophysiology Studies in Patients with Brugada Syndrome
Previous Article in Journal
Joint Optimization of Age of Information and Energy Consumption in NR-V2X System Based on Deep Reinforcement Learning
Previous Article in Special Issue
Monitoring of Curing Process of Epoxy Resin by Long-Period Fiber Gratings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An SNR Enhancement Method for Φ-OTDR Vibration Signals Based on the PCA-VSS-NLMS Algorithm

1
School of Electronic Information Engineering, Changchun University of Science and Technology, Changchun 130000, China
2
Shenyang Institute of Engineering, Shenyang 110136, China
3
State Grid Changchun Electric Power Surrly Company, Changchun 130000, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(13), 4340; https://doi.org/10.3390/s24134340
Submission received: 13 June 2024 / Revised: 27 June 2024 / Accepted: 2 July 2024 / Published: 4 July 2024
(This article belongs to the Special Issue Advances in Applications of Optical Fiber Sensors)

Abstract

:
To improve the signal-to-noise ratio (SNR) of vibration signals in a phase-sensitive optical time-domain reflectometer (Φ-OTDR) system, a principal component analysis variable step-size normalized least mean square (PCA-VSS-NLMS) denoising method was proposed in this study. First, the mathematical principle of the PCA-VSS-NLMS algorithm was constructed. This algorithm can adjust the input signal to achieve the best filter effect. Second, the effectiveness of the algorithm was verified via simulation, and the simulation results show that compared with the wavelet denoising (WD), Wiener filtering, variational mode decomposition (VMD), and variable step-size normalized least mean square (VSS-NLMS) algorithms, the PCA-VSS-NLMS algorithm can improve the SNR to 30.68 dB when the initial SNR is −1.23 dB. Finally, the PCA-VSS-NLMS algorithm was embedded into the built Φ-OTDR system, an 11.22 km fiber was measured, and PZT was added at 10.19–10.24 km to impose multiple sets of fixed-frequency disturbances. The experimental results show that the SNR of the vibration signal is 8.77 dB at 100 Hz and 0.07 s, and the SNR is improved to 26.17 dB after PCA-VSS-NLMS filtering; thus, the SNR is improved by 17.40 dB. This method can improve the SNR of the system’s position information without the need to change the existing hardware conditions, and it provides a new scheme for the detection and recognition of long-distance vibration signals.

1. Introduction

The distributed fiber-optic sensor known as the Φ-OTDR system primarily measures phase changes resulting from fiber vibration or stress [1]. This system boasts a long detection range, exceptional sensitivity, and a rapid response time [2], making it a popular choice for applications such as monitoring oil and gas pipelines [3], assessing the structural health of buildings [4], ensuring perimeter security [5], monitoring power and communication safety [6], sensing rail transit systems [7,8], and monitoring engineering geology [9,10].
Utilizing the Φ-OTDR system involves examining the Rayleigh backscattering (RBS) phase alterations within optical fibers. This RBS, being a feeble signal of minimal intensity [11], becomes susceptible to various environmental disturbances, such as polarization fading from the laser itself, interference fading, and external noise. The resulting accumulation of noise within the phase signal, coupled with a diminished signal-to-noise ratio (SNR), ultimately impedes the precision and range capabilities of the sensing operation [12].
One way to enhance hardware design was suggested by Zhu et al. in 2015. These authors introduced an active compensation approach using laser frequency scanning and cross-correlation calculation to mitigate the impact of light source frequency drifting (LSFD). This method effectively detected disturbances at a frequency of 0.5 Hz up to a distance of 5 km along the sensing fiber [13]. Subsequently, in 2016, Baker C et al., proposed a technique utilizing sinusoidally modulated optical signals (SMOSs) to produce high extinction analog pulses. This innovation minimized noise in the backscattering trajectory of the fiber-under-test (FUT) pulses, thus enhancing the signal-to-noise ratio (SNR) of the Φ-OTDR system [14]. Furthermore, Wang et al., presented a novel approach in 2019 that implemented linearization and Golay pulse coding for heterodyne Φ-OTDR. Through experimental validation, they achieved submeter-level measurement accuracy and nano-strain resolution with a sensing range of 10 km [15]. In 2020, Chen et al., proposed a technique involving forced carrier recombination to enhance the extinction ratio (ER) of the semiconductor optical amplifier (SOA). This resulted in a 5.2 dB increase in the SNR of the Φ-OTDR system [12]. Lastly, Li et al., proposed a strategy in 2023 that integrated multi-mode fiber (MMF) and optimized the “n” elastomer to enhance the Φ-OTDR system. This method led to improvements of 10.51 dB and 13.38 dB in position and frequency SNR enhancement, respectively [16].
One way to improve the Φ-OTDR system is to enhance data processing efficiency. Wavelet denoising methods were initially introduced by Qin et al., in 2012 [17], which was followed by the implementation of continuous wavelet transform denoising techniques in the Φ-OTDR system [18]. Wu et al. introduced a trajectory denoising approach for the Φ-OTDR system in 2015 that leveraged multi-scale wavelet decomposition [19]. Moreover, He et al., made strides in 2017 with an adaptive image restoration algorithm employing two-dimensional bi-lateral filtering to boost the SNR of the Φ-OTDR system. By applying this approach to a 27.6 km sensing fiber, the SNR for position information of a signal with an original SNR of 6.43 dB notably increased by more than 14 dB [20]. Lv et al., made further advancements in 2019 by introducing the empirical mode decomposition (EMD) method to address phase drift issues. Notably, when faced with low-frequency external vibrations at 0.5 Hz or 0.3 Hz, the SNR post-phase signal elimination reached impressive levels of 55.58 dB and 44.44 dB [21]. Building upon these innovations, He et al. presented the complete ensemble empirical mode decomposition with adaptive noise (CEEEMDAN) algorithm in 2020, which achieved a significant 52.11 dB increase in the SNR of disturbed positions [22]. In 2021, Jiang et al., introduced a signal enhancement technique for the Φ-OTDR system utilizing deep learning, resulting in an increase in the mean SNR from 13.4 dB to 42.8 dB [23]. Similarly, in 2021, Ma et al., proposed the VSS-NLMS denoising method which enabled the detection of five sets of vibration events ranging from 100 Hz to 500 Hz at the 10.14 km sensing fiber location. The position information SNR for the VSS-NLMS noise reduction approach was enhanced to 59.31 dB, 46.81 dB, 50.14 dB, 34.00 dB, and 67.09 dB [24]. A method based on kurtosis parameter statistical analysis (k-parameter) was introduced by He et al. in 2022, resulting in an improvement in the SNR at 11.9 km on a 12 km fiber to 5.61 dB [25]. Lastly, Liu et al., proposed a genetic least mean square (GLMS) method in 2023, which demonstrated a notable improvement in the SNR ranging from 14.37 dB to 23.60 dB during monotone-scale audio signal testing at 60~1000 Hz. Nevertheless, this method faces challenges such as a low SNR and limited adaptability, leading to subpar denoising effects [26]. In 2023, Turov et al., proposed to present a nonlinear two-dimensional processing method for distributed acoustic sensor data in both time and frequency domains, with an 8-fold increase in operation speed and an improvement in signal-to-noise ratio of 3.7 dB when used alone, and 10.8 dB when used in combination with moving averaging–moving differential (MA-MD) [27].
In order to further improve the vibration signal SNR of the Φ-OTDR system, we proposed a noise reduction algorithm based on PCA-VSS-NLMS to filter the background noise of the heterodyne coherent Φ-OTDR system. In this study, the mathematical principle of the PCA-VSS-NLMS algorithm was constructed in detail, and simulation experiments were performed to verify the filtering effect of the algorithm. The simulation results show that, compared with the WD, Wiener, VMD, and VSS-NLMS algorithms, the PCA-VSS-NLMS algorithm can increase the SNR to 30.68 dB when the initial SNR is −1.23 dB. The proposed algorithm was embedded into the Φ-OTDR system, an 11.22 km fiber was measured, and multiple sets of fixed-frequency disturbances were imposed by adding PZT at 10.19 to 10.24 km. At 100 Hz and 0.07 s, the SNR of the vibration signal is 8.77 dB. After PCA-VSS-NLMS filtering, the SNR is improved to 26.17 dB; thus, the SNR is improved by 17.40 dB. We verified the feasibility of the proposed method from the perspectives of theoretical research, simulation analysis, and laboratory construction and experimentation. By applying the proposed algorithm, the SNR of vibration signals can be improved without changing the existing hardware conditions, which provides a new idea for long-distance vibration signal recognition.

2. Working Principle of PCA-VSS-NLMS

An automatically adaptive filter can adjust its parameters based on the input signal’s characteristics to accommodate signal fluctuations. This enhances its ability to process signals that are not consistent over time [28]. When comparing their structures, the normalized LMS (NLMS) algorithm mirrors the standard LMS algorithm [29]. As depicted in Figure 1, both are finite impulse response (FIR) filters, with the weight controller mechanism being the distinguishing factor. The output signal M × 1 generated by the input signal x ( n ) of the y ( n ) tap will be subtracted from the expected signal d ( n ) to obtain the error signal e ( n ) . In response to the combined action of the input signal x ( n ) and error signal e ( n ) , weight adjustment to the weight controller is applied to the FIR filter. In a large number of adaptive loops, the weight vector of the filter is adjusted repeatedly until the filter reaches a steady state.
The basic principle of the NLMS algorithm is that w ^ ( n ) represents the old weight vector of the filter during the adaptive cycle n and w ^ ( n + 1 ) represents the new weight vector of the filter during the adaptive cycle n + 1 . Then, the design criterion of the NLMS algorithm can be expressed as a constrained optimization problem [30]:
δ w ^ ( n + 1 ) = w ^ ( n + 1 ) w ^ ( n )
w ^ H ( n + 1 ) x ( n ) = d ( n )
where the superscript H stands for conjugate transpose.
We use the method of Lagrange multipliers. For the general case of complex data, the cost function is estimated as follows:
J ( n ) = δ w ^ ( n + 1 ) 2 + Re [ λ ( d ( n ) w ^ H ( n + 1 ) x ( n ) ) ]
where λ is the complex Lagrange multiplier and is the complex conjugate; Re [ ] represents the operation of taking the natural part, and the constraint’s contribution to the cost function is real-valued; δ w ^ ( n + 1 ) 2 represents the square operation of the Euclidean norm, and the result is also real-valued. Thus, the cost function J ( n ) is a real-valued quadratic function w ^ ( n + 1 ) and can be expressed as follows:
J ( n ) = ( w ^ ( n + 1 ) w ^ ( n ) ) H ( w ^ ( n + 1 ) w ^ ( n ) ) + Re [ λ ( d ( n ) w ^ H ( n + 1 ) x ( n ) ) ]
In order to find the most updated weight vector with the smallest cost function J ( n ) , the following steps are taken:
Take the derivative of the cost function J ( n ) with respect to w ^ H ( n + 1 ) .
J ( n ) w ^ H ( n + 1 ) = 2 ( w ^ ( n + 1 ) w ^ ( n ) ) λ x ( n )
Set it to zero, and the optimal solution is
w ^ ( n + 1 ) = w ^ ( n ) + 1 2 λ x ( n )
Bring (6) into (2) and solve for the unknown multiplier λ .
λ = 2 e ( n ) x ( n ) 2
e ( n ) = d ( n ) w ^ H ( n ) x ( n )
From (6) and (7),
δ w ^ ( n + 1 ) = w ^ ( n + 1 ) w ^ ( n ) = 1 x ( n ) 2 x ( n ) e ( n )
To control the incremental change of the tap weight vector from one adaptive loop to the next without changing the direction of the vector, a positive real scale factor μ ˜ is introduced. The increment is defined as follows:
δ w ^ ( n + 1 ) = w ^ ( n + 1 ) w ^ ( n ) = μ ˜ x ( n ) 2 x ( n ) e ( n )
w ^ ( n + 1 ) = w ^ ( n ) + μ ˜ x ( n ) 2 x ( n ) e ( n )
When the tap vector is small, the smaller square norm x ( n ) 2 has to be divided by μ ˜ , which may cause numerical difficulties. To overcome this problem, (11) is revised to read as follows:
w ^ ( n + 1 ) = w ^ ( n ) + μ ˜ x ( n ) 2 + c x ( n ) e ( n )
where c is a small positive integer.
The traditional NLMS algorithm uses a fixed step size, which cannot be dynamically adjusted according to the characteristics of the current signal and has poor adaptability to the signal’s dynamic characteristics. The VSS-NLMS algorithm can solve this problem very well [24]. We use the step-size update equation in the literature [31], as shown in Equation (13):
μ ˜ ( n ) = β ( 1 exp ( α | e ( n ) | 2 ) )
where μ ˜ ( n ) is the step size of the n iteration; α is a constant to control the step-size range; and β is another constant to adjust the speed of the step length. Usually, α > 0 and β > 0 .
In the VSS-NLMS algorithm, the expected signal d ( n ) plays an important role, and the selection of this signal will directly affect the convergence speed, stability, and filtering effect. However, in practical applications, we usually do not know the real expected signal of a noisy signal. Currently, an input signal or its variant is often used as the expected signal of the filtering system [24].
The vibration signals of the sensing fiber obtained by the Φ-OTDR system are shown in Figure 2. Among them, the detection fiber was 11.22 km, the PZT parameters were set, the output waveform was set to a sine wave, the amplitude value was the peak-to-peak value (VPP) of the signal set to 10 V. We designed a PCA-VSS-NLMS filtering algorithm, which uses PCA to extract the position information features of the sensing fiber instead of the expected signals in VSS-NLMS and then splices the extracted main features with the data processed via PCA-VSS-NLMS filtering. The detailed algorithm flow is shown in Figure 3.
The sample matrix X is represented as follows:
X = [ x 11 x 1 l x t 1 x t l ]
where t represents the sampling time, l represents the spatial position of the sampling point, and x t l represents the amplitude of the phase signal at the l point at time t .
The sample matrix is normalized as follows:
X i j = x i j x ¯ j s i   ( i = 1 , 2 , , t ; j = 1 , 2 , , l )
s j 2 = i = 1 t ( x i j x ¯ j ) t 1   ,   x ¯ j = i = 1 t x i j t
where X i j represents the standardized data matrix, x i j represents the standardized value of the j th index of the i th sample, s i represents the standard deviation, s j 2 represents the variance, and x ¯ j represents the mean value.
The covariance matrix Σ is calculated as follows:
Σ = 1 t ( X i j ) T ( X i j )
The eigenvalue decomposition is performed on the given covariance matrix Σ . This will yield eigenvalues r 1 , r 2 , , r n and their corresponding eigenvectors v 1 , v 2 , , v n , where n is the number of features in the data. The eigenvalues are arranged in descending order. The top k eigenvalues are selected along with their corresponding eigenvectors. These eigenvectors will form the columns of the projection matrix W . Let us denote the selected eigenvectors as v 1 , v 2 , , v k . The selected eigenvectors are then arranged as columns in the projection matrix W . Thus, W = [ v 1 , v 2 , , v k ] , where v 1 , v 2 , , v k are the top k eigenvectors obtained from the eigenvalue decomposition. The original data matrix X is multiplied by the projection matrix W , resulting in the following:
X p c a n e w = X W
where X p c a n e w denotes the newly generated data matrix that can be used as the expected signal of the principal component analysis expectation variable step normalized least mean square (PCA-d-VSS-NLMS) filter.
In the Φ-OTDR system, we pay more attention to the amplitude change of the phase signal. To further improve the SNR of the signal, we perform feature splicing of the feature data points in each row of X p c a n e w with the filtered matrix X p c a d v s s n l m s row by row.
X p c a n e w = [ x p n 1   1 x p n 1   2 x p n 1   l x p n 2   1 x p n 2   2 x p n t   1 x p n t   2 x p n t   l ]
X p c a d v s s n l m s = [ x p d v n 1   1 x p d v n 1   2 x p d v n 1   l x p d v n 2   1 x p d v n 2   2 x p d v n t   1 x p d v n t   2 x p d v n t   l ]
Assuming the data at the position of the 800th column of the 1st row in the X p c a n e w matrix are the feature data points, the data at x p n 1   800 are replaced with the data at x p d v n 1   800 to obtain the final filtered matrix X p c a v s s n l m s .
Filter feature data by setting a threshold to remove noisy data with small absolute values of amplitude. Only the peaks exceeding the set threshold are retained as feature data. In this paper, we utilize the statistical properties of noise to set the adaptive threshold for feature point extraction, assuming that the noise obeys a Gaussian distribution, and calculate the adaptive threshold:
T = μ n + k t σ n
where T is the adaptive threshold, μ n is the estimated noise mean, σ n is the estimated noise standard deviation, and k t is the scale factor. In this paper, k t = 2 is set.

3. Experimental Section

3.1. Simulation Experiments and Results

The SNR calculation is defined as follows:
S N R = 20 log 10 ( A s i g n a l A n o s i e )
where A s i g n a l and A n o s i e are the signal and noise amplitudes, respectively, and we calculate the root mean square of the noise signal as A n o s i e .
According to Figure 4 brought into Equation (19), the SNR for each stage of the proposed algorithm in this paper is calculated as:
S N R r a w = 20 log 10 ( A s i g n a l _ r a w A n o s i e _ r a w ) = 15 . 25   dB S N R p c a = 20 log 10 ( A s i g n a l _ p c a A n o s i e _ p c a ) = 23 . 41   dB S N R p c a d v s s n l m s = 20 log 10 ( A s i g n a l _ p d v n A n o s i e _ p d v n ) = 1 . 36   dB S N R p c a v s s n l m s = 20 log 10 ( A p c a v s s n l m s A p c a v s s n l m s ) = 30 . 33   dB
The SNR is improved by 8.16 dB from the original signal to PCA processing, and 15.08 dB from the original signal to PCA-VSS-NLMS processing.
We performed simulation experiments to verify the effectiveness of the algorithm and compared its performance with other algorithms. First, the unit pulse signal was used as the expected signal, with 2000 sampling points, 0.9 at the 1900th sampling point and 0 at the rest. Gaussian white noise was added to the unit pulse signal as the filter input signal, as shown in Figure 5a. The SNR was 7.74 dB.
When comparing the performance of various filtering algorithms, the wavelet basis function of the WD algorithm was sym8, the wavelet order was 2, and the SNR after filtering was 7.78 dB, as shown in Figure 5b. The Wiener algorithm could not find feature points after filtering, as shown in Figure 5c. In the VMD algorithm, the penalty parameter was 24, the number of modal components was 4, and the SNR after filtering was 6.69 dB, as shown in Figure 5d. In the VSS-NLMS filtering algorithm, the step-size update equation α was 8, β was 0.01, the filter order was 100, and the SNR after filtering was 7.83 dB, as shown in Figure 5e. In the proposed PCA-VSS-NLMS filtering algorithm, the step-size update equation α was 8; β was 0.01; the filter order was 100; the number of retained principal components was 6; k = 1 feature vectors were selected as the principal components, and the selected feature vectors were arranged into a projection matrix by column; and the SNR after filtering was 27.29 dB, as shown in Figure 5f. To further verify the effectiveness of the algorithm, Gaussian white noise with different SNRs was added to the unit pulse signal to compare the SNR enhancement effects of the WD, Wiener, VMD, VSS-NLMS, and PCA-VSS-NLMS algorithms, as shown in Table 1.
The results presented in Table 1 demonstrate the impact of various filtering algorithms on noise SNRs. The data indicate that the PCA-VSS-NLMS algorithm outperforms other algorithms in terms of SNR enhancement across different noise levels. In particular, at an SNR of −1.23 dB, the PCA-VSS-NLMS algorithm successfully preserves a signal with an SNR of 30.68 dB, while other algorithms fail to do so. These findings suggest that the filtering algorithm introduced in this study effectively enhances the SNR and mitigates noise interference.

3.2. Experimental Setup and Results

Figure 6 shows the construction of the Φ-OTDR system. A laser with a narrow linewidth of 1550.12 nm central wavelength and a 3 kHz linewidth is used to split the light into two channels of 90% and 10% using coupler-1. Through an acousto-optic modulator (AOM), 90% of the light is converted into a 200 MHz frequency shifted optical pulse, which is then amplified by an erbium-doped fiber amplifier (EDFA) before entering an optical circulator (OC). At coupler two, the RBS light from the sensing fiber interferes with the local light to produce beat light. The output signal from the beat light is detected by a balanced photodetector (BPD), which transfers it to the computer via a data acquisition card (DAQ) for phase signal extraction.
We use the I/Q quadrature demodulation method proposed in the literature [20] to demodulate Rayleigh’s backscattered light amplitude and phase signals.
The beat frequency signals acquired by the Φ-OTDR system are as follows:
P ( t ) = E R ( t ) E L O ( t ) cos [ 2 π Δ f t + Φ ( t ) + Φ 0 ]
where E R is the backward Rayleigh scattered light, E L O is the reference light, Δ f is the AOM frequency shift, Φ ( t ) is the phase change due to the vibration signal, and Φ 0 is the phase change due to the noise signal.
y 1 = cos ( 2 π Δ f t + Φ r )
y 2 = sin ( 2 π Δ f t + Φ r )
where Φ r is the phase noise, multiplying Equation (24) by Equations (25) and (26), respectively, yields:
I ( t ) = E R ( t ) E L O ( t ) cos [ 2 π Δ f t + Φ ( t ) + Φ 0 ] cos ( 2 π Δ f t + Φ r )                                                                     = 1 2 E R ( t ) E L O ( t ) [ cos [ 4 π Δ f t + Φ ( t ) + Φ 0 + Φ r ] + cos [ Φ ( t ) + Φ 0 Φ r ] ]
Q ( t ) = E R ( t ) E L O ( t ) cos [ 2 π Δ f t + Φ ( t ) + Φ 0 ] sin ( 2 π Δ f t + Φ r )                                                                     = 1 2 E R ( t ) E L O ( t ) [ cos [ 4 π Δ f t + Φ ( t ) + Φ 0 + Φ r ] sin [ Φ ( t ) + Φ 0 Φ r ] ]
The orthogonal signals I ( t ) and Q ( t ) are obtained after low-pass filtering.
I ( t ) = 1 2 E R ( t ) E L O ( t ) cos [ Φ ( t ) + Φ 0 Φ r ]
Q ( t ) = 1 2 E R ( t ) E L O ( t ) sin [ Φ ( t ) + Φ 0 Φ r ]
A s I ( t ) 2 + Q ( t ) 2
Φ ( t ) = arctan Q ( t ) I ( t )
where Φ 0 Φ r is the phase noise, A s is the signal amplitude, and Φ ( t ) is the phase change due to the vibration signal.
We connected the Φ-OTDR system to an 11 km G 652D single-mode fiber. The experimental setup is shown in Figure 7. The optical outlet of the Φ-OTDR system was connected to the first 5 km disc fiber, whose end was fused to the second 5 km disc fiber; the second 5 km disc fiber was connected to the PZT device INPUT, and then the PZT device OUTPUT was connected to a 1 km disc fiber. The actual length of the fiber measured by the OTDR was 11.22 km, the PZT inner winding used a 50 m fiber, and the vibration position was 10.19–10.24 km. The Φ-OTDR system parameters were set as follows: the detection range was set to 11 km, and the spatial resolution was set to 10 m. The PZT parameters are set as follows: the output waveform is set as a sinusoidal waveform, and the amplitude value is the peak-to-peak value of the signal VPP, which is set to 10 V. The frequencies of 100 Hz, 200 Hz, 300 Hz, 400 Hz, 500 Hz, 600 Hz, 700 Hz, 800 Hz, and 900 Hz were set to compare the effects of the filtering algorithms. The experimental results of the PCA-VSS-NLMS algorithm proposed in this paper are shown in Figure 8. System components, PZT, are procured from Nanjing Fiber Photonics Technology Co., Ltd, Nanjing, Jiangsu Province, China. It was accessed on 12 June 2024 at http://www.fib-tech.com/.
The experiment was repeated 10 times for each fixed-frequency signal under the same conditions, and the average SNR at 0.07 s was recorded. The comparison data of filtering effects are shown in Table 2. As can be seen from Table 2, which is much higher than that of the other filtering algorithms.
A three-dimensional (3D) spatio-temporal diagram was drawn for the PZT vibration signal frequency at 100 Hz, as shown in Figure 8. In this figure, Figure 8a is the 3D spatio-temporal diagram of the unfiltered measured data; Figure 8b shows the 3D spatio-temporal diagram of the filtered data after applying the PCA-VSS-NLMS algorithm; Figure 8c shows the top view of the 3D spatio-temporal diagram of the unfiltered measured data; and Figure 8d shows the top view of the 3D spatio-temporal diagram of the filtered data after applying the PCA-VSS-NLMS algorithm. It can be seen from Figure 8 that the vibration position is consistent with the actual position at 10.19–10.24 km. The proposed PCA-VSS-NLMS algorithm can effectively filter out the background noise and highlight the characteristics of vibration signals.
The time-domain signal at the fixed position in Figure 8 is selected, as shown in Figure 9. The red solid line is the curve after PCA-VSS-NLMS processing, and the black dashed line is the curve of the measured data, where Figure 9a is the comparison of the vibration signal before and after the filtering at the vibration position of PZT at 10.23 km. From Figure 9a, it can be seen that the features of the vibration signal are well preserved, and Figure 9b is the comparison before and after filtering of the background noise at 10.00 km. From Figure 9b, it can be seen that the background noise is well suppressed.
The PZT vibration frequency is 100 Hz, and the time-frequency domain diagram at 10.23 km of the fiber is shown in Figure 10, where Figure 10a is the time-frequency domain diagram before filtering; Figure 10b is the time-frequency domain diagram after filtering. A comparison of the spectra of Figure 10a and Figure 10b shows that the low-frequency component is effectively suppressed after filtering.
For the dual-point vibration experiment, the Φ-OTDR system was connected to a 12 km G 652D single-mode fiber. The schematic diagram of the dual-point position vibration Φ-OTDR system is shown in Figure 11. The actual length of the fiber measured by OTDR is 12.25 km, and the Φ-OTDR system and PZT parameter settings are the same as in Figure 6. A 100 Hz frequency vibration is applied at 10.19–10.24 km and a 100 Hz frequency vibration is applied at 12.20–11.25 km, respectively. The experimental results of the PCA-VSS-NLMS algorithm proposed in this paper are shown in Figure 12, where Figure 12a is the 3D spatio-temporal map of the unfiltered measurement data, Figure 12b is the 3D spatial-temporal plot of the PCA-VSS-NLMS filtered data, Figure 12c is the unfiltered vibration signal at 0.07 s, and Figure 12d is the PCA-VSS-NLMS filtered vibration signal at 0.07 s. The 100 Hz vibration point SNR was raised from 9.05 dB to 22.26 dB, a mention of 13.21 dB. The SNR of the second 100 Hz vibration point increased from 8.57 dB to 22.23 dB, an improvement of 13.66 dB.
To further verify the practical application effect of the algorithm proposed in this study, on 14 April 2024, we embedded the PCA-VSS-NLMS algorithm into the Φ-OTDR system and installed the system in the communication room of the 500 kv Station A of Tongliao City, Inner Mongolia Province of China, as part of the state grid. We measured a 90 km optical fiber composite overhead ground wire (OPGW) as the optical cable line. The Φ-OTDR system parameters were set as follows: the detection range was set to 90 km, and the spatial resolution was set to 100 m. The installation and analysis results are shown in Figure 13.
Where Figure 13a is a schematic diagram of the OPGW cable’s location; Figure 13b is the installation diagram of the Φ-OTDR system in the station; Figure 13c shows the top view of the 3D spatial spectrum of unfiltered measured vibration data from the 90 km OPGW fiber-optic cable; and Figure 13d shows the top view of the 3D spatial spectrum of vibration data filtered by PCA-VSS-NLMS from the 90 km OPGW cable. As shown in Figure 13c,d, the proposed filtering algorithm can effectively filter out the background noise and highlight the vibration signal position. As can be seen in Figure 13d, there is a strong vibration signal at a distance of 46–49 km, which can be viewed as a multi-point vibration. It was verified that the 46–49 km line was in a level 2 dance zone and that the meteorological data for that day were a southerly wind at level 5.

4. Conclusions

In this study, we enhanced the position information SNR of the Φ-OTDR system by introducing a filtering algorithm based on PCA-VSS-NLMS. The mathematical foundation of the PCA-VSS-NLMS algorithm was elucidated, and its effectiveness was established through simulation experiments. The results from the simulation experiments show that the PCA-VSS-NLMS algorithm achieves a significant improvement in SNR, reaching up to 30.68 dB when the initial SNR is only −1.23 dB. It outperforms existing algorithms such as WD, Wiener, VMD, and VSS-NLMS, highlighting the potential of the proposed algorithm in enhancing the performance of the Φ-OTDR system. The PCA-VSS-NLMS algorithm was embedded into the built Φ-OTDR system, an 11.22 km fiber was measured, and PZT was added at 10.19–10.24 km to impose multiple sets of fixed-frequency disturbances. The experimental results show that the SNR of the vibration signal is 8.77 dB at 100 Hz and 0.07 s, and the SNR is improved to 26.17 dB after PCA-VSS-NLMS filtering; thus, the SNR is improved by 17.40 dB. In addition, we carried out practical application measurements to monitor the vibration of a 90 km OPGW as an optical cable line of the 500 kv Station A in Tongliao City, Inner Mongolia Province of China, as part of the state grid, effectively reducing the background noise. The measurement results are consistent with the actual situation. The proposed algorithm can improve the SNR of the Φ-OTDR system’s position information without changing the existing hardware conditions and provides a new scheme for the detection and recognition of long-distance vibration signals.

Author Contributions

X.C.: conceptualization, supervision, and writing-review and editing. H.Y.: conceptualization, methodology, software, investigation, and writing the original draft. J.X.: software and investigation. F.G.: software and investigation. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the Jilin Province Science and Technology Development Plan Project (20210203044SF) and the Capital Construction Fund in the Jilin Provincial Budget in 2022 (Innovation Capacity Building Project) (2022C045-8).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available from the corresponding author upon request. The data are not publicly available due to potential commercial values.

Conflicts of Interest

Author Funan Gao was employed by State Grid Changchun Electric Power Surrly Company. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Liu, Z.; Lu, W.; Liu, L.; Shieh, W. Phase Noise-Induced Interference for Coherently Detected OTDR Systems. Opt. Lett. 2024, 49, 766. [Google Scholar] [CrossRef] [PubMed]
  2. Yu, F.-H.; Liu, S.; Shao, L.; Xu, W.; Xiao, D.; Zhao, F.; Hu, J.; Lin, W.; Wang, G.; Wang, W.; et al. Ultra-Low Sampling Resolution Technique for Heterodyne Phase-OTDR Based Distributed Acoustic Sensing. Opt. Lett. 2022, 47, 3379. [Google Scholar] [CrossRef] [PubMed]
  3. Ashry, I.; Mao, Y.; Wang, B.; Hveding, F.; Bukhamsin, A.; Ng, T.K.; Ooi, B.S. A Review of Distributed Fiber–Optic Sensing in the Oil and Gas Industry. J. Light. Technol. 2022, 40, 1407–1431. [Google Scholar] [CrossRef]
  4. Zahoor, R.; Cerri, E.; Vallifuoco, R.; Zeni, L.; De Luca, A.; Caputo, F.; Minardo, A. Lamb Wave Detection for Structural Health Monitoring Using a Φ-OTDR System. Sensors 2022, 22, 5962. [Google Scholar] [CrossRef] [PubMed]
  5. Lyu, C.; Niu, Z.; Tian, J.; Jin, J.; Yang, J.; Ge, C. Identification of Intrusion Events Based on Distributed Optical Fiber Sensing in Complex Environment. IEEE Internet Things J. 2022, 9, 24212–24220. [Google Scholar] [CrossRef]
  6. Ding, Z.-W.; Zhang, X.-P.; Zou, N.-M.; Xiong, F.; Song, J.-Y.; Fang, X.; Wang, F.; Zhang, Y.-X. Phi-OTDR Based On-Line Monitoring of Overhead Power Transmission Line. J. Light. Technol. 2021, 39, 5163–5169. [Google Scholar] [CrossRef]
  7. Yang, N.; Zhao, Y.; Chen, J. Real-Time Φ-OTDR Vibration Event Recognition Based on Image Target Detection. Sensors 2022, 22, 1127. [Google Scholar] [CrossRef]
  8. Wang, Z.; Jiang, J.; Wang, Z.; Xiong, J.; Rao, Y.-J. Bandwidth-Enhanced Quasi-Distributed Acoustic Sensing With Interleaved Chirped Pulses. IEEE Sens. J. 2020, 20, 12739–12743. [Google Scholar] [CrossRef]
  9. Alekseev, A.E.; Gorshkov, B.G.; Potapov, V.T.; Taranov, M.A.; Simikin, D.E. A Fiber Phase-Sensitive Optical Time-Domain Reflectometer for Engineering Geology Application. Instrum. Exp. Tech. 2023, 66, 843–848. [Google Scholar] [CrossRef]
  10. Hartlieb, S.; Ringkowski, M.; Haist, T.; Sawodny, O.; Osten, W. Multi-Positional Image-Based Vibration Measurement by Holographic Image Replication. Light Adv. Manuf. 2021, 2, 1. [Google Scholar] [CrossRef]
  11. Horiguchi, T.; Shimizu, K.; Kurashima, T.; Tateda, M.; Koyamada, Y. Development of a Distributed Sensing Technique Using Brillouin Scattering. J. Light. Technol. 1995, 13, 1296–1302. [Google Scholar] [CrossRef]
  12. Chen, Y.; Mao, B.-M.; Zhou, B.; Guo, C.; Lin, Z. Improving the SNR of the Phase-OTDR by Controlling the Carrier in the SOA. J. Mod. Opt. 2020, 67, 1241–1246. [Google Scholar] [CrossRef]
  13. Zhu, F.; Zhang, X.; Xia, L.; Guo, Z.; Zhang, Y. Active Compensation Method for Light Source Frequency Drifting in Φ-OTDR Sensing System. IEEE Photon. Technol. Lett. 2015, 27, 2523–2526. [Google Scholar] [CrossRef]
  14. Baker, C.; Vanus, B.; Wuilpart, M.; Chen, L.; Bao, X. Enhancement of Optical Pulse Extinction-Ratio Using the Nonlinear Kerr Effect for Phase-OTDR. Opt. Express 2016, 24, 19424. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, Z.; Zhang, B.; Xiong, J.; Fu, Y.; Lin, S.; Jiang, J.; Chen, Y.; Wu, Y.; Meng, Q.; Rao, Y. Distributed Acoustic Sensing Based on Pulse-Coding Phase-Sensitive OTDR. IEEE Internet Things J. 2019, 6, 6117–6124. [Google Scholar] [CrossRef]
  16. Qin, Z.; Chen, L.; Bao, X. Wavelet Denoising Method for Improving Detection Performance of Distributed Vibration Sensor. IEEE Photon. Technol. Lett. 2012, 24, 542–544. [Google Scholar] [CrossRef]
  17. Qin, Z.; Chen, L. Continuous Wavelet Transform for Non- Stationary Vibration Detection with Phase-OTDR. Opt. Express 2012, 20, 20459–20465. [Google Scholar] [CrossRef] [PubMed]
  18. Wu, H.; Xiao, S.; Li, X.; Wang, Z.; Xu, J.; Rao, Y. Separation and Determination of the Disturbing Signals in Phase-Sensitive Optical Time Domain Reflectometry (Φ-OTDR). J. Light. Technol. 2015, 33, 3156–3162. [Google Scholar] [CrossRef]
  19. He, H.; Shao, L.; Li, H.; Pan, W.; Luo, B.; Zou, X.; Yan, L. SNR Enhancement in Phase-Sensitive OTDR with Adaptive 2-D Bilateral Filtering Algorithm. IEEE Photonics J. 2017, 9, 1–10. [Google Scholar] [CrossRef]
  20. Li, J.; Liu, H.; Shen, X.; Xiao, Y.; Wu, Z.; Guo, P.; Hu, J.; Liu, Y.; Dang, H.; Sun, Q.; et al. SNR Enhancement of Quasi-Distributed Weak Acoustic Signal Detection by Elastomers and MMF Integrated Φ-OTDR. Opt. Express 2023, 31, 37019. [Google Scholar] [CrossRef]
  21. Lv, Y.; Wang, P.; Wang, Y.; Liu, X.; Bai, Q.; Li, P.; Zhang, H.; Gao, Y.; Jin, B. Eliminating Phase Drift for Distributed Optical Fiber Acoustic Sensing System with Empirical Mode Decomposition. Sensors 2019, 19, 5392. [Google Scholar] [CrossRef] [PubMed]
  22. He, M.; Feng, L.; Zhao, D. A Method to Enhance SNR Based on CEEMDAN and the Interval Thresholding in Φ_OTDR Systems. Appl. Phys. B 2020, 126, 97. [Google Scholar] [CrossRef]
  23. Jiang, F.; Zhang, Z.; Lu, Z.; Li, H.; Tian, Y.; Zhang, Y.; Zhang, X. High-Fidelity Acoustic Signal Enhancement for Phase-OTDR Using Supervised Learning. Opt. Express 2021, 29, 33467. [Google Scholar] [CrossRef] [PubMed]
  24. Ma, Q.; Gao, X.; Gao, Y.; Zhang, X.; Zhong, Z. A Study on Noise Reduction of Φ-OTDR System Based on VSS-NLMS Algorithm. IEEE Sens. J. 2021, 21, 7648–7656. [Google Scholar] [CrossRef]
  25. He, Q.; Zeng, Z.; Zhao, Q.; Shang, X.; Li, T. SNR Improvement of Vibration Sensing in a Conventional Phase-OTDR by k-Parameter Statistical Analysis. Opt. Commun. 2022, 509, 127789. [Google Scholar] [CrossRef]
  26. Liu, X.; Liu, Z.; Zhou, X.; Wang, Y.; Bai, Q.; Jin, B. Signal-to-Noise Ratio Improvement for Phase-Sensitive Optical Time-Domain Reflectometry Using a Genetic Least Mean Square Method. Photonics 2023, 10, 1362. [Google Scholar] [CrossRef]
  27. Turov, A.T.; Barkov, F.L.; Konstantinov, Y.A.; Korobko, D.A.; Lopez-Mercado, C.A.; Fotiadi, A.A. Activation Function Dynamic Averaging as a Technique for Nonlinear 2D Data Denoising in Distributed Acoustic Sensors. Algorithms 2023, 16, 440. [Google Scholar] [CrossRef]
  28. Gueraini, I.; Benallal, A.; Tedjani, A. New Variable Step-Size Fast NLMS Algorithm for Non-Stationary Systems. SIViP 2023, 17, 3099–3107. [Google Scholar] [CrossRef]
  29. Cai, Y.; Yu, Z.; Mo, D.; Liu, R.; Chen, A.; Dai, B.; Li, Y. Noise Reduction with Adaptive Filtering Scheme on Interferometric Fiber Optic Hydrophone. Optik 2020, 211, 164648. [Google Scholar] [CrossRef]
  30. Haykin, S. Adaptive Filter Theory Fifth Edition; Pearson Education: Upper Saddle River, NJ, USA, 2014. [Google Scholar]
  31. Li, Y.; Hamamura, M. Zero-attracting Variable-step-size Least Mean Square Algorithms for Adaptive Sparse Channel Estimation. Adapt. Control Signal 2015, 29, 1189–1206. [Google Scholar] [CrossRef]
Figure 1. Working principle of adaptive filtering.
Figure 1. Working principle of adaptive filtering.
Sensors 24 04340 g001
Figure 2. Vibration signals acquired by the Φ-OTDR system: (a) a 100 Hz vibration signal was applied from 10.19 to 10.24 km; (b) signals taken in absolute value after applying 100 Hz vibration normalization at 10.19 to 10.24 km.
Figure 2. Vibration signals acquired by the Φ-OTDR system: (a) a 100 Hz vibration signal was applied from 10.19 to 10.24 km; (b) signals taken in absolute value after applying 100 Hz vibration normalization at 10.19 to 10.24 km.
Sensors 24 04340 g002
Figure 3. Flowchart of the PCA-VSS-NLMS algorithm.
Figure 3. Flowchart of the PCA-VSS-NLMS algorithm.
Sensors 24 04340 g003
Figure 4. Plot of the effect of each filtering process of PCA-VSS-NLMS at 0.07 s with 100 Hz vibration applied from 10.19 km to 10.24 km: (a) Signal to be filtered; (b) PCA-filtered signals; (c) PCA-d-VSS-NLMS-filtered signals; (d) PCA-VSS-NLMS-filtered signals.
Figure 4. Plot of the effect of each filtering process of PCA-VSS-NLMS at 0.07 s with 100 Hz vibration applied from 10.19 km to 10.24 km: (a) Signal to be filtered; (b) PCA-filtered signals; (c) PCA-d-VSS-NLMS-filtered signals; (d) PCA-VSS-NLMS-filtered signals.
Sensors 24 04340 g004
Figure 5. Comparison of the filtering effects of different algorithms: (a) signal to be filtered; (b) WD denoising; (c) Wiener denoising; (d) VMD denoising; (e) VSS-NLMS denoising; and (f) PCA-VSS-NLMS denoising.
Figure 5. Comparison of the filtering effects of different algorithms: (a) signal to be filtered; (b) WD denoising; (c) Wiener denoising; (d) VMD denoising; (e) VSS-NLMS denoising; and (f) PCA-VSS-NLMS denoising.
Sensors 24 04340 g005aSensors 24 04340 g005b
Figure 6. Schematic diagram of the Φ-OTDR system.
Figure 6. Schematic diagram of the Φ-OTDR system.
Sensors 24 04340 g006
Figure 7. The experimental environment of the Φ-OTDR system.
Figure 7. The experimental environment of the Φ-OTDR system.
Sensors 24 04340 g007
Figure 8. 3D spatio-temporal maps of measured vibration signals: (a) 3D spatio-temporal map of unfiltered vibration signals at 11.22 km; (b) 3D spatio-temporal map of filtered vibration signals at 11.22 km PCA-VSS-NLMS; (c) 3D spatio-temporal map top view of unfiltered vibration signals at 11.22 km; (d) 3D spatio-temporal map top view of filtered vibration signals at 11.22 km PCA-VSS-NLMS.
Figure 8. 3D spatio-temporal maps of measured vibration signals: (a) 3D spatio-temporal map of unfiltered vibration signals at 11.22 km; (b) 3D spatio-temporal map of filtered vibration signals at 11.22 km PCA-VSS-NLMS; (c) 3D spatio-temporal map top view of unfiltered vibration signals at 11.22 km; (d) 3D spatio-temporal map top view of filtered vibration signals at 11.22 km PCA-VSS-NLMS.
Sensors 24 04340 g008
Figure 9. Fixed-position time-domain signal plots: (a) before and after filtering at 10.23 km; (b) before and after filtering at 10.00 km.
Figure 9. Fixed-position time-domain signal plots: (a) before and after filtering at 10.23 km; (b) before and after filtering at 10.00 km.
Sensors 24 04340 g009
Figure 10. Time-frequency domain plot of PZT vibration at 100 Hz, 10.23 km from the optical fiber: (a) Time-frequency domain plot before filtering; (b) Time-frequency domain plot after filtering.
Figure 10. Time-frequency domain plot of PZT vibration at 100 Hz, 10.23 km from the optical fiber: (a) Time-frequency domain plot before filtering; (b) Time-frequency domain plot after filtering.
Sensors 24 04340 g010
Figure 11. Schematic diagram of the two-point vibration experiment of the Φ-OTDR system.
Figure 11. Schematic diagram of the two-point vibration experiment of the Φ-OTDR system.
Sensors 24 04340 g011
Figure 12. 3D spatio-temporal maps of two-point vibration experiments: (a) 3D spatio-temporal maps of unfiltered measured data; (b) 3D spatiotemporal maps of PCA-VSS-NLMS-filtered data; (c) unfiltered signals at 0.07 s; (d) PCA-VSS-NLMS-filtered signals at 0.07 s.
Figure 12. 3D spatio-temporal maps of two-point vibration experiments: (a) 3D spatio-temporal maps of unfiltered measured data; (b) 3D spatiotemporal maps of PCA-VSS-NLMS-filtered data; (c) unfiltered signals at 0.07 s; (d) PCA-VSS-NLMS-filtered signals at 0.07 s.
Sensors 24 04340 g012
Figure 13. Installation and analysis results: (a) schematic diagram of the OPGW cable’s location; (b) installation diagram of the Φ-OTDR system in the station; (c) top view of the 3D spatio-temporal diagrams of unfiltered measured vibration data from the 90 km OPGW fiber-optic cable; and (d) top view of the 3D spatio-temporal diagrams of vibration data filtered by PCA-VSS-NLMS from the 90 km OPGW fiber-optic cable.
Figure 13. Installation and analysis results: (a) schematic diagram of the OPGW cable’s location; (b) installation diagram of the Φ-OTDR system in the station; (c) top view of the 3D spatio-temporal diagrams of unfiltered measured vibration data from the 90 km OPGW fiber-optic cable; and (d) top view of the 3D spatio-temporal diagrams of vibration data filtered by PCA-VSS-NLMS from the 90 km OPGW fiber-optic cable.
Sensors 24 04340 g013
Table 1. Comparison of the filtering effects of SNR algorithms with different noise signals.
Table 1. Comparison of the filtering effects of SNR algorithms with different noise signals.
Noisy Signal SNR
Filtering
Algorithm
−1.23 dB4.88 dB7.74 dB10.24 dB12.51 dB14.96 dB17.19 dB19.82 dB
WDNo signalNo signal7.8812.7215.3617.8519.8722.73
WienerNo signalNo signalNo signal8.3210.4113.2316.7921.64
VMDNo signalNo signal6.696.469.6712.4114.1817.18
VSS-NLMSNo signal5.047.8610.3912.8515.0917.2419.78
PCA-VSS-NLMS30.6826.6527.2927.8928.4929.6030.6231.55
Table 2. Comparison of filtering effects of different filtering algorithms at different frequencies.
Table 2. Comparison of filtering effects of different filtering algorithms at different frequencies.
Filtering
Algorithm
100 Hz
SNR/dB
200 Hz SNR/dB300 Hz SNR/dB400 Hz SNR/dB500 Hz SNR/dB600 Hz SNR/dB700 Hz SNR/dB800 Hz SNR/dB900 Hz SNR/dB
Unfiltered8.779.369.9310.0311.169.277.7010.447.91
WD11.4710.4012.1811.7911.8212.6710.8813.5910.36
Wiener6.027.697.978.179.317.525.079.445.28
VMD6.616.347.348.768.379.466.089.248.57
VSS-NLMS8.929.039.959.9510.929.197.5810.207.27
PCA-VSS-NLMS26.1726.4025.5224.5226.9225.3928.3526.0126.92
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chen, X.; Yu, H.; Xu, J.; Gao, F. An SNR Enhancement Method for Φ-OTDR Vibration Signals Based on the PCA-VSS-NLMS Algorithm. Sensors 2024, 24, 4340. https://doi.org/10.3390/s24134340

AMA Style

Chen X, Yu H, Xu J, Gao F. An SNR Enhancement Method for Φ-OTDR Vibration Signals Based on the PCA-VSS-NLMS Algorithm. Sensors. 2024; 24(13):4340. https://doi.org/10.3390/s24134340

Chicago/Turabian Style

Chen, Xiaojuan, Haoyu Yu, Jingyao Xu, and Funan Gao. 2024. "An SNR Enhancement Method for Φ-OTDR Vibration Signals Based on the PCA-VSS-NLMS Algorithm" Sensors 24, no. 13: 4340. https://doi.org/10.3390/s24134340

APA Style

Chen, X., Yu, H., Xu, J., & Gao, F. (2024). An SNR Enhancement Method for Φ-OTDR Vibration Signals Based on the PCA-VSS-NLMS Algorithm. Sensors, 24(13), 4340. https://doi.org/10.3390/s24134340

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop