Research on algorithm and simulation of acoustic signal denoising based on improved wavelet threshold

SHI Xuewei; XU Dalin; LIU Zhicheng; XU Zhiyan

doi:10.3969/j.issn.1673-3819.2024.05.010

Command Control and Simulation ›› 2024, Vol. 46 ›› Issue (5) : 69-76. DOI: 10.3969/j.issn.1673-3819.2024.05.010

Multimodal Information Fusion

Research on algorithm and simulation of acoustic signal denoising based on improved wavelet threshold

Author information +

History +

Abstract

Due to the low signal-to-noise ratio of raw data, the reliability of data and the accuracy of acoustic source localization are severely affected by fiber optic acoustic sensing technology. To address this issue, this study optimizes the wavelet thresholding method. Firstly, a novel thresholding function is proposed, which achieves denoising while preserving key information through shape adjustment factors. It combines the advantages of both hard and soft threshold functions and has high flexibility and controllability. Secondly, an adaptive threshold calculation method is introduced, utilizing an improved simulated annealing algorithm to optimize threshold selection, reducing the algorithm’s dependence on threshold parameter selection. Through simulation experiments, it has been verified that this research method effectively suppresses noise in the signal and improves data availability. Compared to the original methods, this approach significantly improves the signal-to-noise ratio and demonstrates robustness in simulated tests of real signals.

Key words

fiber optic acoustic sensing technology; wavelet thresholding denoising; thresholding function; simulated annealing algorithm

Cite this article

EndNote

Ris (Procite)

Bibtex

Download Citations

SHI Xuewei , XU Dalin , LIU Zhicheng , XU Zhiyan. Research on algorithm and simulation of acoustic signal denoising based on improved wavelet threshold. Command Control and Simulation. 2024, 46(5): 69-76 https://doi.org/10.3969/j.issn.1673-3819.2024.05.010

分布式光纤声波传感技术(Distributed Fiber Optic Sensing)已经成为监测和测量领域的一项重要技术,它广泛应用于地震监测、油气勘探、安全监测、周界安防等^[1⇓-3]等多个领域。然而,与许多监测技术一样,分布式光纤声波传感数据常常面临多种振动源产生的不同信号的复杂混合,特别是在诸如隧道挖掘和采矿工程等复杂地下环境中,这些信号容易受到各种随机噪声的污染^[4],这对准确分析和解释声波信号的特征具有很大影响。分布式光纤声波传感技术目前面临的一个重要问题就是数据的信噪比太低。信噪比的低下直接影响声波数据的可靠性,进而影响声波源位置的准确性^[5]。在过去的几年中,已经研究了多种数据去噪方法来增强光纤声波传感信号的信噪比,例如数字滤波器^[6]、经验模态分解^[7]、奇异值分解^[8]、小波变换^[9]和深度学习^[10]。然而,数字滤波器、经验模态分解和奇异值分解方法通常用于滤除高频噪声成分,当信号和噪声共享近似频带时,无法突出有用信息。深度学习方法可以对高质量的微震数据实现有效去噪,但其效果高度依赖训练数据集。

自小波阈值去噪算法提出以来,就因其优越的分解性能而广泛应用于多个领域^[11]。它可以将信号分解为近似信息和详细信息,滤除信号中的噪声。然而,它很难确定信号和噪声系数之间的明确边界,其中的阈值规则和阈值函数对算法的性能影响很大^[12]。因此,针对各种特定应用有不同的改进策略^[13]。为了克服传统和现有改进的小波阈值去噪方法在去噪时的缺陷,本文提出了两点改进:一是利用自适应阈值计算方法优化阈值选择,二是提出并采用一种新的阈值函数。改进算法可以抑制高频和低频噪声,进一步提高声波信号的信噪比,为了验证改进算法的有效性和鲁棒性,笔者设计了两组仿真实验。

1 小波阈值算法

1.1 小波阈值去噪

假设一维噪声信号模型可以表示如下:

x(t)=h(t)+ε(t)

(1)

其中,x(t)是混合信号,h(t)是原始信号,ε(t)是噪声。在最简单的情况下,ε(t)可以假设为高斯白噪声。实际工程中,h(t)通常被认为是低频信号或相对稳定的信号,而噪声信号通常是高频信号。小波变换的目的就是抑制高频的无用信号并恢复低频的有用部分。小波阈值的去噪过程包括小波变换分解、阈值处理和信号重构。其原理如图1所示。

Fig.1 Principle of wavelet threshold denoising

图1 小波阈值去噪原理

Full size|PPT slide

对于确定的分解水平,可以使用离散小波变换将信号分解为近似系数序列和细节系数序列。近似系数代表低频分解,通常相当于有用信号,细节系数代表中高频分解,通常被视作噪声。图2显示了一维信号的典型3级小波分解示意图。其中,Y(t)为原始信号,CA为低频信息,CD为高频小波系数。

Fig.2 Typical 3-level wavelet decomposition process

图2 典型3级小波分解的过程

Full size|PPT slide

经过小波分解处理后,有用信号的小波系数幅度远大于噪声的幅度。当小波分解水平增加时,有用信号的小波系数的幅度很少改变,但噪声的小波系数幅度迅速衰减。小波分解的这个特性保证了其可以构造阈值和阈值函数来处理分解系数,从而达到消除噪声的目的。最后,通过应用逆小波变换重建去噪信号。

小波阈值去噪的过程包括小波分解、阈值收缩过程和小波重构。小波分解涉及小波基和分解水平的选择,阈值收缩处理需要考虑阈值计算和阈值函数。

1.2 改进阈值函数

小波阈值是将小波系数与给定阈值进行比较。如果小波系数小于阈值,则小波系数置零。如果小波系数大于该阈值,则应保留或修改小波系数。常用的阈值函数有硬阈值函数和软阈值函数。

软阈值函数为

{\hat{W}}_{J, K}

\{\begin{array}{l} S g n (W_{J, K}) (|W_{J, K}| - λ), | W_{J, K} | \geq λ \\ 0, | W_{J, K} | < λ \end{array}

(2)

硬阈值函数为

{\hat{W}}_{J, K}

\{\begin{array}{l} W_{J, K}, | W_{J, K} | \geq λ \\ 0, | W_{J, K} | < λ \end{array}

(3)

在式(2)和(3)中,W_J_,_K是处理前的小波系数;

{\hat{W}}_{J, K}

是处理后的小波系数;λ是设置的阈值;Sgn(·)是符号函数;J是最大分解尺度,j=1,2,…,J;K是小波系数的长度,k=1,2,…,n_j;n_j=N/2J-j+1;N是信号的长度。

软阈值函数通常会对小波系数进行平滑处理来减小噪声,然而,平滑性过强会导致信号的细节信息丢失,从而导致过度去噪。并且,实际的小波系数和处理后的小波系数始终存在恒定的偏差,这会降低重建信号的质量。硬阈值函数会将小于阈值的小波系数置零,从而导致信号中的小幅度信号被完全消除,造成有效信息的丢失。此外,由于硬阈值函数直接将小波系数截断,在信号重建过程中会出现振荡现象,进一步产生伪吉布斯现象。

为了克服软阈值函数和硬阈值函数的上述缺点,本文提出了一种改进的阈值函数,这种改进的阈值函数可以通过改变参数因子,进行灵活调整,使函数值位于硬阈值和软阈值函数值之间。改进的阈值函数如下:

{\hat{W}}_{J, K} = {\begin{cases} W_{J, K} - Sgn (W_{J, K}) \times \frac{β λ}{e^{α (| W_{J, K} | - λ)}}, & | W_{J, K} | ⩾ λ \\ 0, & | W_{J, K} | ＜ λ \end{cases}

(4)

式中,α和β是调节因子,用于调节

{\hat{W}}_{J, K}

的变化速度,其中,α用于调节函数向渐近线的逼近程度,取值范围为[0,∞],β用于调节阈值λ阶段点的y值,取值范围为[0,1]。图3展示了当β为1时,函数曲线随α的变化情况,随着α的增大,改进的阈值函数更加逼近软阈值函数。图4展示了当α为1时,函数曲线随β的变化情况,随着β的增大,改进的阈值函数在截止点的阈值更接近硬阈值函数。

Fig.3 The curve follows α changes generated when β is 1

图3 β为1时,曲线随α的变化

Full size|PPT slide

Fig.4 The curve follows β changes generated when α is 1

图4 α为1时,曲线随β的变化

Full size|PPT slide

本文提出的改进阈值函数通过调节α和β两个调节因子,就可以更便捷地控制阈值函数,在消除阈值函数不连续性的同时,使函数更迅速靠近期望的位置。改进方法保证了阈值函数的连续性,不会导致突变的产生,从而避免了伪吉布斯现象的发生,使得降噪后信号能够保持原始信号的平滑性。同时,改进方法可以减小由于小波分解重构对信号压缩带来的偏差,从而提高重构信号与真实信号的逼近程度。

1.3 基于改进模拟退火算法优化的阈值计算

阈值常用的选择方法主要有固定阈值、极值阈值、无偏似然估计阈值、极大极小阈值和启发式阈值。但上述方法需要根据信号和噪声的性质及应用需求进行调整,对于小波系数而言,全局最优的阈值并不一定局部最优。为了克服这一缺点,本文采用了一种基于改进模拟退火算法的阈值计算方法来优化阈值计算。

模拟退火算法是一种启发式优化算法,用于解决各种组合优化和搜索问题。但是传统的模拟退火算法会存在丢失最优解、容易陷入局部最优解等问题。

为了解决上述传统模拟退火算法存在的缺点,进一步提高算法的搜索能力,笔者结合目标小波阈值的特点,对算法做出以下改进:

1) 记录全局最优解,解决原始算法由于随机扰动造成的最优解丢失的问题。

2) 采用正态分布作为扰动的邻域分布,选择当前解和最优解的差值作为方差。在远离最优解时,大方差会带来大扰动,当靠近最优解时,小方差避免跳出最优解。

3) 采用方差阈值的方式,连续多次方差小于阈值则增大方差跳出当前状态,防止陷入局部最优解。增大方差时回火升温,避免因为温度过低导致搜索能力变差。

改进算法的流程图如图5所示。

Fig.5 Flowchart of improved simulated annealing algorithm

图5 改进模拟退火算法流程图

Full size|PPT slide

该研究通过改进模拟退火算法,提高了算法的寻优能力,减少搜索陷入局部最优解的情况;通过将阈值的计算转化为基于信噪比的优化,解决了阈值选择问题;通过对优化算法的改进,增强了其搜索能力。

2 仿真实验及分析

2.1 噪声分析

本次实验将真实声波数据添加噪声,针对改进算法进行仿真验证,对实验结果进行分析,为真实信号的处理提供了参考。为了更好地近似光纤声波传感的真实信号,在仿真实验中以微软可定制语音数据库MS-SNSD中的无噪声波信号作为模拟原始信号。向模拟信号添加噪声,噪声选择数据集NoiseX-92中的白噪声以及分布式光纤传感声波信号稳态值模拟的真实噪声。

白噪声常常用来模拟自然环境中的背景噪声,它包含了所有频率的噪声成分。真实噪声包含系统使用过程中存在的噪声。添加不同类型的噪声有助于验证去噪算法的鲁棒性等性能。不同类型的噪声反映了实际环境中存在的干扰,使算法能够应对各种复杂的去噪场景。通过在仿真实验中添加这些噪声,可以评估去噪算法在不同噪声条件下的效果。

2.2 评价指标

实验中,使用信噪比(SNR)、均方误差(MSE)、波形互相关系数(CC)定量地比较降噪算法的性能。SNR可以反映降噪算法的去噪效果,MSE和CC可以反映去噪后信号和真实信号之间的接近程度。信噪比、均方误差和波形相似度的计算公式如下:

SNR

(d B)

=10log₁₀

\frac{P_{s i g n a l}}{P_{n o i s e}}

=10log₁₀

\frac{\overset{n}{\sum_{i = 1}} Y_{i}^{2}}{\overset{n}{\sum_{i = 1}} (Y_{i} - {\hat{Y}}_{i})^{2}}

(5)

MSE=

\frac{1}{n} \overset{n}{\sum_{i = 1}}

(Y_i-

{\hat{Y}}_{i}

)²

(6)

CC=R_xy(m)=

\frac{1}{N} \overset{N}{\sum_{k = 1}}

x(k)y(k-m)

(7)

其中,Y_i是原始信号幅值,

{\hat{Y}}_{i}

是去噪信号的幅值,x(k)为原始信号序列,y(k-m)为去噪信号序列,m为两个序列的移位长度,N为信号序列的长度。

2.3 实验参数设置

1)小波基选择

不同的小波基具有不同的时频特性。为选择最适合声波信号去噪的小波基,本文对仿真信号选取了不同的小波基函数,计算了降噪后的信噪比,基于上述降噪评价规则对比了最优小波基函数。通过添加高斯白噪声(SNR为5 dB)对信号进行仿真,以测试小波基haar、Sym、db、coif和bior函数对声波信号降噪的影响。各小波基去噪后的信噪比如表1所示。表1中展示了各小波基在相同条件下对仿真信号的去噪结果,各小波系仅记录去噪效果最好的小波。其中,sym小波系中的sym8小波对仿真信号去噪效果最好,因此后续的仿真实验采用sym8小波作为小波阈值算法的小波基。

Tab.1 Denoising effect of various wavelet bases

表1 各小波基去噪效果

小波基	信噪比/dB
haar	5.036 0
sym8	7.775 9
db9	7.705 4
coif2	7.286 6
bior2.4	6.897 5

2)分解层数选择

为选择最适合声波信号去噪的分解层数,本文对仿真信号选取了不同的分解层数,计算了降噪后的信噪比,通过添加高斯白噪声(SNR为5 dB)对信号进行仿真,以测试不同分解层数对声波信号降噪的影响。表2中展示了各分解层数在sym8小波基下对仿真信号的去噪效果。从表2可以看出,随着分解层数的提高,信噪比先增大后减小,这是因为分解层数过少时,信号去噪不彻底,去噪效果不理想。更高的分解层数能够带来更好的去噪效果,但是分解层数过高时,信号的有效信息丢失严重,反而会导致信噪比的降低,还会使得运算量增大,算法的处理速度变慢。对本文仿真信号而言,选择3层分解层数可以取得最高的信噪比,因此后续的仿真实验采用3层分解层数。

Tab.2 Noise reduction effect of each decomposition level

表2 各分解层数去噪效果

分解层数	信噪比/dB
2	9.124 7
3	10.167 0
4	7.775 9
5	4.317 9

2.4 实验结果分析

为了验证本文所提方法对小波阈值去噪算法的提升效果,笔者利用数据集中的声波信号进行了噪声消除的仿真实验,对实验结果进行了分析。为了更好地接近光纤声波传感信号的实际信号,在仿真实验中以数据集声波信号作为模拟信号,向模拟信号添加噪声,将具有高斯分布的白噪声和真实光纤声波传感系统噪声作为噪声添加到模拟信号中。信噪比设置为 5 dB。原始信号及包含噪声的信号波形图如图6所示。原始信号及包含噪声的信号频谱图如图7所示。

Fig.6 Waveform of raw signal and noisy signal

图6 原始信号及含噪信号的波形图

Full size|PPT slide

Fig.7 Spectrum of the original signal and noisy signal

图7 原始信号及含噪信号的频谱图

Full size|PPT slide

与硬阈值函数小波变换和软阈值函数小波变换去噪方法相比,采用改进的阈值函数及改进模拟退火算法优化阈值的小波变换方法进行去噪的对比效果如图8、图9所示。

Fig.8 Effect of 5 dB white noise denoising

图8 5 dB白噪声去噪效果图

Full size|PPT slide

Fig.9 Effect of 5 dB real noise denoising

图9 5 dB真实噪声去噪效果图

Full size|PPT slide

表3和表4中分别展示了加入5 dB信噪比的不同类型噪声后,本文改进方法和硬、软小波阈值去噪算法的信噪比、均方误差和波形互相关系数的值。从表中可以看出,本文所提方法在白噪声和真实噪声的仿真实验中三项数据的表现均优于原始算法。从信噪比的值可以看出,本文所提方法对小波阈值算法的去噪效果提升显著。

Tab.3 Add 5 dB white noise denoising indicators

表3 加入5 dB白噪声去噪指标

去噪算法	信噪比/dB	均方误差	波形互相关系数
小波硬阈值去噪	8.670 4	0.010 2	0.925
小波软阈值去噪	10.167 0	0.008 9	0.934
本文方法	16.758 3	0.003 1	0.962

Tab.4 Adding 5 dB real noise denoising indicators

表4 加入5 dB真实噪声去噪指标

去噪算法	信噪比/dB	均方误差	波形互相关系数
小波硬阈值去噪	5.674 7	0.157 6	0.831
小波软阈值去噪	7.573 4	0.011 9	0.914
本文方法	10.897 2	0.008 4	0.935

从表4中可以看出,在加入真实噪声后,本文所提算法依然能取得更好的效果,说明算法具有更强的鲁棒性,能够应对更加真实的噪声情况。

图10展示了不同信噪下两种噪声去噪后的趋势变化。

Fig.10 Signal-to-noise ratio of two types of noise denoising under different signal-to-noise ratios

图10 不同信噪比下两种噪声去噪后的信噪比

Full size|PPT slide

从图10可以看出,随着信噪比的提高,算法的去噪能力逐渐变弱,这是由于信噪比较低时,信号中含有的噪声分量较多,去噪带来的信噪比提升较大;而信噪比高时,信号中含有的噪声分量较少,去噪时容易过度去噪,造成有效信息的损失,从而使信噪比提升幅度变小,其中的小波硬阈值去噪方法在20 dB信噪比进行去噪处理时反而会导致信噪比的下降。

3 结束语

本文提出了一种新的阈值函数和阈值计算方法用于改进小波阈值去噪算法,旨在对光纤传感声波数据取得更好的去噪效果。新的阈值函数解决了软、硬阈值函数存在的缺点,保证了阈值函数的连续性。该函数结合经典的硬阈值函数和软阈值函数,通过形状调整因子调整整体速度,即用α和β调整

{\hat{W}}_{J, K}

的增长率来使其逼近软阈值函数和硬阈值函数,并且可以指定阈值处的截断点,具有更好的灵活性和可控性。本文提出的新阈值计算方法自适应地求解最优阈值,减小了小波阈值去噪方法对于阈值选择的依赖性,解决了小波阈值去噪算法过度依赖参数选择的情况。综合数据测试结果表明,所提改进方法能够抑制高频和低频噪声,有利于声波信号的准确识别,提高了声波信号的可用性。此外,与传统的小波阈值去噪方法相比,本文所提方法能够显著提高两种仿真声波信号的信噪比。在真实光纤声波传感信号噪声的仿真中进一步证实了所提出的改进小波阈值去噪方法可以有效对真实信号进行处理,证明了所提算法的鲁棒性。

References

Publishing order | Descend order by publishing year | Descend order by cited within

[1]	GONZALEZ-HERRAEZ M. Distributed acoustic sensing for seismic monitoring[C]// 2021 Optical Fiber Communications Conference and Exhibition (OFC), San Francisco, CA, USA, 2021. 本文引用 [1]

[2]	SUN W H, ZHU S, LI W, et al. Noise suppression of distributed acoustic sensing based on f-x deconvolution and wavelet transform[J]. IEEE Photonics Journal, 2020, 12(1): 1-8. 本文引用 [1]

[3]	HICKE K, HUSSELS M-T, EISERMANN R, et al. Condition monitoring of industrial infrastructures using distributed fiber optic acoustic sensors[C]// 2017 25th Optical Fiber Sensors Conference (OFS), Jeju, Korea (South), 2017. 本文引用 [1]

[4]	李柯, 龚岩栋, 张卓. φ-OTDR降噪处理的研究进展[J]. 激光与光电子学进展, 2022, 59(23): 25-33. LI K, GONG Y D, ZHANG Z. φ- Research progress of φ-OTDR noise reduction processing[J]. Laser & Optoelectronics Progress, 2022, 59(23): 25-33. 本文引用 [1]

[5]	张芳. 扩频光时域反射技术理论与空间定位应用研究[D]. 秦皇岛: 燕山大学, 2021. ZHANG F. Research on the theory and spatial positioning application of spread spectrum optical time domain reflection technology[D]. Qinhuangdao: Yanshan University, 2021. 本文引用 [1]

[6]	汪豫. 基于Ф-OTDR的光纤振动传感关键技术研究[D]. 成都: 电子科技大学, 2018. WANG Y. Study on the optical fiber vibration sensing technology based on Ф-OTDR[D]. Chengdu: University of Electronic Science and Technology of China, 2018. 本文引用 [1]

[7]	蒋浩. 基于Φ-OTDR的管道第三方破坏监测技术研究[D]. 杭州: 浙江大学, 2023. JIANG H. Base on Φ- Research on third-party pipeline damage monitoring technology for OTDR[D]. Hangzhou: Zhejiang University, 2023. 本文引用 [1]

[8]	吕龙宝, 庞拂飞, 刘奂奂, 等. 基于奇异值分解的Φ-OTDR相位提取[J]. 光通信技术, 2018, 42(9): 51-54. LYU L B, PANG F F, LIU H H, et al. Phase-extraction Ф-OTDR based on SVD[J]. Optical Communication Technology, 2018, 42(9): 51-54. 本文引用 [1]

[9]	齐晓辉, 宋宛鸿, 李文召. 一种改进的小波去噪方法在OTDR中的应用[J]. 光通信技术, 2021, 45(4): 54-58. QI X H, SONG W H, LI W Z. Application of an improved wavelet denoising method in OTDR[J]. Optical Communication Technology, 2021, 45(4): 54-58. 本文引用 [1]

[10]	WANG P F, LV Y J, WANG Y, et al. Adaptability and anti-noise capacity enhancement for Φ-OTDR with deep learning[J]. Journal of Lightwave Technology, 2020, 38(23): 6 699-6 706. 本文引用 [1]

[11]	陈霄, 徐慨, 王潋, 等. 基于小波去噪和三阶全息图的时延估计算法[J]. 指挥控制与仿真, 2015, 37(4): 69-73. CHEN X, XU K, WANG L, et al. Time delay estimation algorithm based on wavelet denoising and third order hologram[J]. Command Control & Simulation, 2015, 37(4): 69-73. 本文引用 [1]

[12]	WANG C, DENG C X, HU Z B. An improved wavelet threshold function and its application in image edge detection[C]// 2019 International Conference on Wavelet Analysis and Pattern Recognition (ICWAPR), Kobe, 2020: 1-7. 本文引用 [1]

[13]	LI C, CAO D G, YUAN Y Y. Research on improved wavelet denoising method for sEMG signal[C]// 2019 Chinese Automation Congress (CAC), Hangzhou, 2020: 5 221-5 225. 本文引用 [1]

PDF(1962 KB)

1085

Accesses

Citation

Detail

Sections

Recommended

Abstract
Key words
Cite this article
1 小波阈值算法
1.1 小波阈值去噪
Fig.1 Principle of wavelet threshold denoising
Fig.2 Typical 3-level wavelet decomposition process
1.2 改进阈值函数
Fig.3 The curve follows α changes generated when β is 1
Fig.4 The curve follows β changes generated when α is 1
1.3 基于改进模拟退火算法优化的阈值计算
Fig.5 Flowchart of improved simulated annealing algorithm
2 仿真实验及分析
2.1 噪声分析
2.2 评价指标
2.3 实验参数设置
Tab.1 Denoising effect of various wavelet bases
Tab.2 Noise reduction effect of each decomposition level
2.4 实验结果分析
Fig.6 Waveform of raw signal and noisy signal
Fig.7 Spectrum of the original signal and noisy signal
Fig.8 Effect of 5 dB white noise denoising
Fig.9 Effect of 5 dB real noise denoising
Tab.3 Add 5 dB white noise denoising indicators
Tab.4 Adding 5 dB real noise denoising indicators
Fig.10 Signal-to-noise ratio of two types of noise denoising under different signal-to-noise ratios
3 结束语
References

Received	Revised	Published
2023-10-09	2023-10-27	2024-10-10
Online
2024-10-10

Please choose a citation manager

Content to export