Leakage sound wave noise reduction method and system based on wavelet packet threshold and VMD

Abstract

The invention provides a leakage sound wave noise reduction method and system based on a wavelet packet threshold and VMD, and belongs to the field of filtering noise reduction. The method comprises: acquiring a leakage sound wave acquisition signal; screening key factors of a wavelet packet threshold method based on a high signal-to-noise ratio, and performing multi-scale decomposition on the leakage sound wave acquisition signal based on each factor to obtain a preprocessed signal; an adaptive weight factor and a local disturbance and global jump variation strategy are introduced into the original whale algorithm to form an enhanced whale algorithm; taking the minimum envelope entropy as a fitness function, and adopting an enhanced whale algorithm to optimize a VMD method to obtain an optimal parameter combination; optimizing variational mode decomposition by adopting an optimal parameter combination, and decomposing the preprocessed signal to obtain an intrinsic mode function component; and screening and reconstructing the intrinsic mode function component to realize noise reduction of leakage information. By combining a wavelet packet threshold method with variational mode decomposition, noise can be efficiently reduced, leakage sound wave features can be accurately extracted, and the signal analysis quality is improved.

Description

Technical Field

[0001] The present invention relates to the technical field of pipeline leakage signal filtering and noise reduction, and in particular to a leakage sound wave noise reduction method and system based on wavelet packet threshold and VMD. Background Art

[0002] There are many leakage monitoring methods that can be applied to oil and gas pipelines. Among them, the acoustic wave method has many advantages over the traditional mass balance method, negative pressure wave method, transient model method, etc.: high sensitivity, high positioning accuracy, low false alarm rate, short detection time, strong adaptability; it measures the slight dynamic pressure changes in the pipeline fluid, which is independent of the absolute value of the pipeline operating pressure; the response frequency is wider, the detection range is wider, etc. However, in the actual operation process, the leakage is inevitably affected by the pipeline operation noise, background noise, instrument and equipment noise or other working noise, which makes it difficult to extract effective leakage attenuation characteristics from the pipeline signal. In response to this problem, domestic and foreign scholars have conducted a lot of research. According to the survey results, the current domestic and foreign patents involving oil and gas pipeline leakage acoustic wave noise reduction methods are mainly: US application US6389881 discloses a pipeline leak detection technology based on acoustic wave technology. This technology uses sensors to collect dynamic pressure data in the pipeline and uses pattern matching filtering technology to process the signal to effectively suppress noise and interference, thereby improving the accuracy of leak location; Chinese application ZL2022104774909 discloses a pipeline leakage acoustic signal detection method that optimizes ICEEMDAN. This technology uses the ICEEMDAN algorithm to decompose the collected acoustic signal, calculate the multi-scale permutation entropy MPE of each modal function IMF, select the IMF that meets the requirements according to the obtained MPE value for reconstruction, and obtain the denoised acoustic signal; Chinese application ZL2015100201556 discloses a method for locating oil and gas pipeline leaks based on acoustic wave amplitude. This method uses low-frequency acoustic wave amplitudes obtained after wavelet analysis to detect and locate leaks, and proposes a method for locating leaks that does not consider sound speed and time difference. Chinese application ZL2017114441573 discloses a method for identifying multiphase pipeline leakage acoustic wave signals based on wavelet signals. The technology uses wavelet packet decomposition to decompose the collected multiphase pipeline leakage acoustic wave signals to obtain wavelet signals. The energy extreme value curve of the wavelet signal under each frequency band is calculated; the energy extreme value on the energy extreme value curve under the same conditions is judged to determine whether there is a sudden change, and the leakage signal is identified.

[0003] However, the existing technology only applies the single processing method of wavelet packet threshold method or variational mode decomposition, and lacks the fusion technology of the two methods. At the same time, the traditional wavelet packet threshold method and variational mode decomposition method cannot eliminate the influence that may be caused by the empirical setting parameters. Specifically, it is manifested as: (1) The wavelet packet threshold method can obtain high-frequency signals while extracting low-frequency signals, solving the problem of some features of leaked information being lost in high-frequency components during wavelet decomposition. At the same time, the signal can be set to reduce noise by setting a threshold. However, if the key factors of the wavelet packet threshold method - wavelet basis function, decomposition level, threshold rule and threshold function - are not properly selected, it will directly lead to poor decomposition results.

[0004] (2) The variational mode decomposition method can adaptively realize the frequency domain division of the signal and the effective separation of each component, so it can be applied to the noise reduction of oil and gas pipeline leakage signals. In practical applications, the noise intensity is relatively large, and the single use of the variational mode decomposition method will result in a large computational cost.

[0005] (3) The penalty factor α and the modal decomposition number K in the variational mode decomposition method have the most significant impact on the accuracy of signal decomposition. If the parameters are improperly preset manually, the target feature information will be masked or lost. At the same time, the most commonly used solution is to use a meta-heuristic algorithm based on the component eigenvalue as the criterion for judging the decomposition effect and continuously iterate and optimize to obtain the optimal parameter combination. However, this type of algorithm has many problems of its own. For example, the original whale algorithm has problems with low population richness and premature convergence. Summary of the invention

[0006] In order to solve the above problems, the present invention proposes a leakage sound wave denoising method and system based on wavelet packet threshold and VMD. The key factors of the wavelet packet threshold method are screened based on the high signal-to-noise ratio rule, and the original signal is preprocessed using the wavelet packet threshold method, which can effectively remove the noise components in the signal while retaining the key signal characteristics; the original whale algorithm is introduced with adaptive weight factors, local perturbations and global jump mutation strategies to form an enhanced whale algorithm, and it is introduced into the variational mode decomposition process, with the minimum envelope entropy as the fitness function, and the key parameters such as the decomposition layer number and the penalty factor are adaptively adjusted to ensure that the variational mode decomposition achieves the best performance.

[0007] In order to achieve the above object, the present invention adopts the following technical solution: In a first aspect, the present invention provides a leakage sound wave denoising method based on wavelet packet threshold and VMD, comprising: Acquire leakage acoustic wave acquisition signals; Based on high signal-to-noise ratio, the key factors of the wavelet packet threshold method are selected, and the key factors include the optimal wavelet packet basis function, the number of decomposition layers, the threshold function and the threshold rule; based on the various factors in the wavelet packet threshold method, the leakage acoustic wave acquisition signal is multi-scale decomposed, and the correlation coefficient method is used to select and reconstruct the node coefficients of the last layer to obtain the preprocessed signal; Taking the minimum envelope entropy as the fitness function, the enhanced whale algorithm is used to optimize the variational mode decomposition method to obtain the optimal parameter combination; The optimal parameter combination is used to construct an optimized variational mode decomposition, and the preprocessed signal is decomposed to obtain the intrinsic mode function components; the intrinsic mode function components are screened and reconstructed using the correlation coefficient method to achieve denoising of leakage information.

[0008] Preferably, the composition parameters of the wavelet packet threshold method based on high signal-to-noise ratio screening specifically include: Obtain the original signal of the sample, and perform three rounds of wavelet decomposition on the original signal of the sample based on different wavelet packet basis functions, decomposition levels, threshold functions and threshold rules according to the control variable idea, and select the optimal parameters based on the maximum signal-to-noise ratio of the reconstructed signal and the original signal of the sample; Among them, in the first round, the number of decomposition layers, threshold functions and threshold rules are fixed, and the initial wavelet packet basis functions are screened out based on a variety of wavelet packet basis functions; the second round fixes the threshold functions and threshold rules, and the optimal wavelet packet basis functions and the optimal number of decomposition layers are screened out based on the initial wavelet packet basis functions and different numbers of decomposition layers; the third round fixes the optimal wavelet packet basis functions and the optimal number of decomposition layers, and the optimal threshold functions and the optimal threshold rules are screened out based on a variety of threshold functions and a variety of threshold rules.

[0009] Preferably, the multi-scale decomposition of the leakage acoustic wave acquisition signal based on various parameters in the wavelet packet threshold method, and the selection and reconstruction of the last layer of node coefficients by the correlation coefficient method to obtain the preprocessed signal specifically includes: The optimal wavelet packet basis function, decomposition layer number, threshold function and threshold rule are used to perform multi-scale decomposition on the leakage acoustic wave acquisition signal, and the Pearson correlation coefficient between the node coefficients of the last layer and the leakage acoustic wave acquisition signal is calculated; The nodes whose correlation is greater than or equal to a preset value are reconstructed to obtain a preprocessed signal.

[0010] Preferably, the method uses the minimum envelope entropy as the fitness function and adopts the enhanced whale algorithm to optimize the variational mode decomposition method to obtain the optimal parameter combination, which specifically includes: Initialize the enhanced whale algorithm and determine the initial whale individual based on the minimum envelope entropy; each initial whale individual represents a set of parameter combinations ; The enhanced whale algorithm is used to update the initial individual whale positions; Variational modal decomposition is performed on each whale individual at each updated position, the sum of the envelope entropy of the modal components after decomposition of each whale individual is calculated, and the whale individual with the smallest envelope entropy is recorded; when the maximum number of iterations is reached, the whale individual with the smallest sum of envelope entropy is output as the optimal whale individual; the optimal whale individual represents the optimal parameter combination, including the optimal number of decomposition layers and the penalty factor.

[0011] Preferably, the step of initializing the enhanced whale algorithm and determining the initial whale individual based on the minimum envelope entropy comprises: Set the number of whale populations, generate individual whales, and optimize the target parameters Range restriction; the individual whales represent a group Parameter combination, where is the number of decomposition layers, is the penalty factor; Initialize the parameters of the variational mode decomposition method, perform variational mode decomposition on each whale individual, calculate the sum of the envelope entropies of the modal components of the decomposed signal, and record and save the whale individual with the minimum envelope entropy as the initial whale individual.

[0012] Preferably, the use of the enhanced whale algorithm to update the initial individual whale positions specifically includes: Generate a random number p, and determine whether the random number p is greater than or equal to 0.5; If yes, a weight factor is introduced and a spiral update mechanism is selected to update the individual positions of the whales; If not, continue to judge Is it greater than 1? If so, a random search mechanism is selected to update the individual positions of the whales; if not, a differential mutation perturbation factor is introduced, and a shrinkage and encirclement mechanism is selected to update the individual positions of the whales.

[0013] Preferably, the optimal parameter combination is used to construct an optimized variational mode decomposition, the preprocessed signal is decomposed to obtain the intrinsic mode function components; the intrinsic mode function components are screened and reconstructed by the correlation coefficient method to achieve denoising of leakage information, specifically including: Inputting the optimal parameter combination into a variational mode decomposition algorithm, decomposing the preprocessed signal, and calculating the Pearson correlation coefficient between each intrinsic mode function component and the leakage acoustic wave acquisition signal; The modal function components whose correlation is greater than or equal to a preset value are reconstructed to obtain the noise-reduced leakage sound wave signal.

[0014] In a second aspect, the present invention provides a leakage sound wave noise reduction system based on wavelet packet threshold and VMD, comprising: The original signal acquisition module is used to obtain the leakage sound wave acquisition signal; The signal preprocessing module is used to screen the key factors of the wavelet packet threshold method based on a high signal-to-noise ratio, wherein the key factors include the optimal wavelet packet basis function, the number of decomposition layers, the threshold function and the threshold rule; the leakage acoustic wave acquisition signal is decomposed at multiple scales based on the various factors in the wavelet packet threshold method, and the last layer node coefficients are screened and reconstructed using the correlation coefficient method to obtain the preprocessed signal; The algorithm optimization module is used to use the minimum envelope entropy as the fitness function and the enhanced whale algorithm to optimize the variational mode decomposition method to obtain the optimal parameter combination; The denoising module is used to construct an optimized variational mode decomposition using the optimal parameter combination, decompose the preprocessed signal to obtain the intrinsic mode function components; and screen and reconstruct the intrinsic mode function components using the correlation coefficient method to achieve denoising of leakage information.

[0015] In a third aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of a leakage acoustic wave denoising method based on wavelet packet threshold and VMD described in the first aspect.

[0016] In a fourth aspect, the present invention provides a computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the steps of a leakage acoustic wave denoising method based on wavelet packet threshold and VMD described in the first aspect are implemented.

[0017] Compared with the prior art, the present invention has the following beneficial effects: 1. The present invention is based on the wavelet packet threshold method and the optimized variational mode decomposition fusion method. It has high accuracy in processing target signals, and the correlation with the original component is 99.9%. It can realize the estimation and identification of oil and gas pipeline leakage signals, and provides a new idea for noise reduction of leakage acoustic wave signals. 2. The improved enhanced whale algorithm of the present invention better balances the global exploration search and local development search capabilities, effectively solving the problems of low group richness and premature convergence. At the same time, only one set of parameters needs to be applied to the enhanced whale algorithm to obtain the optimal parameters, which can greatly reduce the time to find the optimal parameters; 3. Compared with the traditional variational mode decomposition method, the variational mode decomposition method after enhanced whale optimization in the present invention can adaptively select the optimal parameter combination, better divide the time domain and frequency domain of the signal, and obtain the effective decomposition components of the leakage acquisition signal; 4. The signal preprocessing method based on the wavelet packet threshold method provided by the present invention solves the difficult problems of determining the wavelet packet basis function, the number of decomposition levels, the threshold function and the threshold rule, and at the same time frees up a large amount of data cost for the optimized variational mode decomposition method.

[0018] Advantages of additional aspects of the present invention will be given in part in the following description, and in part will become obvious from the following description, or will be learned through practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0019] The accompanying drawings constituting a part of the present invention are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their description are used to explain the present invention but do not constitute a limitation of the present invention.

[0020] Figure 1 A main flow chart of a leakage sound wave denoising method based on wavelet packet threshold and VMD provided by an embodiment of the present invention; Figure 2 A detailed flow chart of a leakage sound wave denoising method based on wavelet packet threshold and VMD provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0021] The present invention will be further described below in conjunction with the accompanying drawings and embodiments.

[0022] Embodiment 1 like Figure 1 As shown, this embodiment discloses a leakage sound wave denoising method based on wavelet packet threshold and VMD, comprising the following steps: S1: Acquire leakage acoustic wave acquisition signal; S2: Based on high signal-to-noise ratio, the key factors of the wavelet packet threshold method are selected, and the key factors include the optimal wavelet packet basis function, the number of decomposition layers, the threshold function and the threshold rule; based on the various factors in the wavelet packet threshold method, the leakage acoustic wave acquisition signal is decomposed at multiple scales, and the correlation coefficient method is used to select and reconstruct the node coefficients of the last layer to obtain the preprocessed signal; S3: Taking the minimum envelope entropy as the fitness function, the enhanced whale algorithm is used to optimize the variational mode decomposition method to obtain the optimal parameter combination; S4: Using the optimal parameter combination to construct an optimized variational mode decomposition, decompose the preprocessed signal to obtain the intrinsic mode function components; using the correlation coefficient method to screen and reconstruct the intrinsic mode function components to achieve noise reduction of leakage information.

[0023] Next, combine Figure 2 , a leakage sound wave denoising method based on wavelet packet threshold and VMD disclosed in this embodiment is described in detail.

[0024] In S1, a leakage acoustic wave acquisition signal is obtained, which specifically includes: setting an acoustic wave sensor on the measured pipeline, collecting signals from the leakage point through the sensor, and obtaining a leakage acoustic wave acquisition signal.

[0025] In S2, preferably, the key factors of the wavelet packet threshold method are screened based on a high signal-to-noise ratio, the leakage acoustic wave acquisition signal is decomposed at multiple scales based on the various factors in the wavelet packet threshold method, and the correlation coefficient method is used to screen and reconstruct the last layer of node coefficients to obtain a preprocessed signal; the specific steps are: S201: First, obtain the sample original signal; S202: Optimizing wavelet packet basis functions: fixing the number of decomposition layers, threshold function and threshold rule, respectively using N wavelet packet basis functions to perform wavelet packet decomposition on the original signal of the sample, and reconstructing the last layer of low-frequency signal obtained by decomposition to obtain a first reconstructed signal; calculating the signal-to-noise ratio of the first reconstructed signal and the original signal of the sample, and selecting the first M wavelet packet basis functions as the initial wavelet packet basis functions based on the maximum signal-to-noise ratio; wherein 2≤M<N.

[0026] As an implementation method, the number of decomposition layers is fixed to 4, the threshold function is a hard threshold, the threshold rule is heursure, and multiple groups of wavelet packet bases of sym, db, coif, fk, bior, and rbio clusters are used for wavelet packet decomposition, and the low-frequency signal of the last layer of decomposition is reconstructed, and the signal-to-noise ratio of the reconstructed signal and the original signal is calculated. The first five wavelet packet basis functions are selected according to the maximum signal-to-noise ratio rule, and the signal-to-noise ratio calculation formula is as follows: (1) Where: To decompose the power of the low-frequency signal in the last layer, is the power of the original sample signal.

[0027] S203: Optimizing the number of decomposition layers: fixing the threshold function and the threshold rule, and setting the number of decomposition layers to 1 to a layers; using the M initial wavelet packet basis functions, performing wavelet packet decomposition on the original signal of the sample based on the number of decomposition layers of 1 to a layers, and reconstructing the last layer of low-frequency signal obtained by decomposition to obtain a second reconstructed signal; calculating the signal-to-noise ratio of the second reconstructed signal and the original signal of the sample, and selecting the optimal wavelet packet basis function and the optimal number of decomposition layers based on the maximum signal-to-noise ratio; wherein a>1.

[0028] As an implementation method, the fixed threshold function is a hard threshold, the threshold rule is heursure, the wavelet basis function selects the 5 wavelet packet basis functions selected in S202, the decomposition layer number is set to 1 to 10 layers, and the 5 wavelet packet basis functions are processed by the wavelet packet threshold method with the decomposition layer number of 1 to 10 layers respectively, the low-frequency signal of the last decomposition layer is reconstructed, the signal-to-noise ratio of the reconstructed signal to the original signal is calculated, and the optimal wavelet packet basis function and decomposition layer number are selected according to the maximum signal-to-noise ratio rule.

[0029] S204: Optimizing threshold functions and threshold rules: setting b threshold functions and c threshold rules; using the optimal wavelet packet basis function and the optimal number of decomposition layers, performing wavelet decomposition on the original signal of the sample based on different threshold functions and threshold rules, and reconstructing the last layer of low-frequency signal obtained by decomposition to obtain a third reconstructed signal; calculating the signal-to-noise ratio of the third reconstructed signal and the original signal of the sample, and screening out the optimal threshold function and the optimal threshold rule based on the maximum signal-to-noise ratio; wherein b>1, c>1.

[0030] As an implementation method, the optimal wavelet packet basis function and the number of decomposition layers in S203 are fixed, the threshold function selects a soft threshold and a hard threshold respectively, and the threshold rules select four threshold rules of rigrsure, heursure, minimaxi, and sqtwolog respectively. The low-frequency signals of the last layer of decomposition under the two threshold functions and the four threshold rules are reconstructed respectively, and the signal-to-noise ratio of the reconstructed signal and the original signal is calculated, and the threshold function and the threshold rule are selected according to the rule of maximum signal-to-noise ratio.

[0031] Through the above three rounds of optimization process, this embodiment can construct a wavelet packet threshold method based on high signal-to-noise ratio screening. By accurately selecting the wavelet packet basis function, the number of decomposition levels, the threshold function and the threshold rule, the noise component in the signal can be more effectively removed while retaining the useful information in the signal, thereby improving the quality and availability of the signal.

[0032] S205: The optimal wavelet packet basis function and decomposition level in S203 and the optimal threshold function and threshold rule in step S204 form the parameters of the wavelet packet threshold method, and perform multi-scale decomposition on the collected signal.

[0033] S206: Calculate the Pearson correlation coefficient between the coefficients of each node in the last layer and the leakage acoustic wave acquisition signal. The calculation formula is: (2) In the formula, is the node coefficient, Leakage sound wave acquisition signal, , is the corresponding mathematical expectation, and the value range of the Pearson correlation coefficient is between 0 and 1; S207: Reconstruct the nodes whose correlation is greater than or equal to a preset value to obtain a preprocessed signal. The preset value is preferably 0.7.

[0034] In S3, the minimum envelope entropy is used as the fitness function, and the enhanced whale algorithm is used to optimize the variational mode decomposition method to obtain the optimal parameter combination. The specific steps are: S301: Set the number of whale populations and optimize the target parameters The range is limited and the initial whale individuals are generated, and the variational mode decomposition parameters are initialized. After calculating the variational mode decomposition of each whale individual, the sum of the envelope entropies of the signal modal components is calculated, and the whale individual with the minimum envelope entropy is recorded. The envelope entropy calculation formula is: (3) In the formula, Represented as preprocessed signal The envelope signal sequence obtained after Hilbert demodulation; yes The normalized form of is the corresponding envelope entropy value.

[0035] S302: Generate a random number P in the range of 0 to 1 and update the individual position of the whale; S3021: If P < 0.5, introduce a weight factor to the original whale algorithm, enter global exploration, and calculate parameters , when | When |>1, the whale enters the random search phase. The random search formula is: (4) Where: t is the current iteration number; is the currently randomly selected position vector; is the position vector of the individual whale; is the distance between the individual whale and its prey; | | is the absolute value; · is the element-by-element multiplication; and The size is controlled by the parameter To decide, w is the weight factor that changes with the number of iterations. The specific calculation is as follows: (5) (6) (7) (8) Where: and is a random real number in (0, 1); is a factor that decreases gradually from 2 to 0; is the maximum number of group iterations; S3022: If P < 0.5, | When |≤1, the original whale algorithm enters the shrinkage and encirclement stage of introducing differential mutation perturbation factors. The formula is as follows: (9) (10) Where: is the currently randomly selected position vector, is the differential variation perturbation factor, F is the variation scale factor; S3023: If P ≥ 0.5, the optimization algorithm enters the local optimization process, and changes the position update method of the original whale algorithm from spiral update to Archimedean spiral update. The individual whale moves towards the current best individual whale position. The formula is as follows: (11) Where: b is a constant coefficient, l is a random number between [-1, 1].

[0036] In this embodiment, the adaptive selection of the three modes of global exploration, shrinking and encircling phase (introducing differential mutation) and local optimization is determined based on the current state and parameters of the algorithm (such as the distance between the whale and the prey, the value of the random number P, etc.). This adaptive mechanism enables the algorithm to flexibly adjust its search strategy in different search phases and scenarios, so as to solve the optimization problem more effectively. Therefore, it can be said that the three modes constitute an adaptive search mechanism, forming an enhanced whale algorithm (EWOA) that can maintain global search capabilities while also having strong local search capabilities.

[0037] Dynamic selection is made based on the algorithm's operating status and the relationship between whales and prey. Global exploration and local optimization are responsible for the algorithm's global search and local search capabilities respectively, while the shrinking and encircling phase enhances the diversity of the population by introducing differential mutation. This design enables the EWOA algorithm to maintain global search capabilities while also having strong local search capabilities, thereby effectively solving optimization problems.

[0038] S303: Update the latest whale individual, perform variational modal decomposition on each whale individual, calculate the sum of the envelope entropy of the modal components after decomposition of each whale individual, and record the whale individual with the minimum envelope entropy. When the maximum number of iterations is reached, stop the iteration and output the whale individual with the minimum sum of envelope entropy. The individual contains the optimal number of decomposition layers. K and penalty factor α The variational mode decomposition steps are: S3031: Convert the decomposition of the original signal into a variational problem, perform Hilbert transform on the mode, and calculate its analytical signal to obtain a single-sided spectrum. The signal is then shifted to baseband to align with its spectrum by multiplying by a set of exponential factors to adjust to the estimated center frequency: (12) The square of the L2 norm of the demodulated signal gradient is calculated as an estimate of the bandwidth of each mode, and a Gaussian smoothing exponent is applied. The constrained variational model of VMD is obtained as: (13) (14) In the formula, represents the original input signal, , Respectively represent the decomposition K（ The number of decomposition levels in the current whale individual) modal functions and the central frequency corresponding to each mode.

[0039] S3032: By introducing a quadratic penalty factor α and the augmented Lagrangian function λ , transforming the constrained variational problem into an unconstrained variational problem, thereby solving the above-mentioned constrained optimization problem.

[0040] (15) S3033: Use the alternating direction multiplier method to find the "saddle point" minimum of the extended Lagrangian expression. The optimal solution is the minimum of each mode. , The iteration ends when the following constraints are met: (16) S3034: Divide the frequency bands based on the frequency domain characteristics of the signal to achieve adaptive decomposition of the signal; (17) In S4, the optimal parameter combination is used to construct an optimized variational mode decomposition, and the preprocessed signal is decomposed to obtain the intrinsic mode function components; the intrinsic mode function components are screened and reconstructed by using the correlation coefficient method to achieve denoising of leakage information. The specific steps are: S401: inputting the optimal parameter combination into a variational mode decomposition algorithm, decomposing the multi-scale acoustic wave signal, and calculating the Pearson correlation coefficient between each intrinsic mode function component and the leakage acoustic wave acquisition signal; S402: Reconstruct the modal function components whose correlation is greater than or equal to a preset value to obtain a noise-reduced leakage sound wave signal. The preset value is preferably 0.7.

[0041] This embodiment creatively integrates the wavelet packet threshold method and the variational mode decomposition method. First, through the high signal-to-noise ratio screening mechanism of the wavelet packet threshold method, the key parameters of the wavelet packet decomposition are accurately determined, and the multi-scale fine decomposition of the collected signal is realized. The last layer of highly correlated nodes are screened and reconstructed through the correlation coefficient method to obtain the preprocessed signal, effectively removing the noise component and retaining the key characteristics of the signal. Subsequently, the enhanced whale algorithm is used to optimize the parameters of the variational mode decomposition, further improving the accuracy and efficiency of the decomposition, making the inherent mode function components obtained by the decomposition more accurate. Finally, the correlation coefficient method is used to screen and reconstruct the modal components with high correlation, which achieves deep noise reduction of the leakage information and significantly improves the quality and analyzability of the signal.

[0042] Embodiment 2 This embodiment provides a leakage sound wave noise reduction system based on wavelet packet threshold and VMD, including: The original signal acquisition module is used to obtain the leakage sound wave acquisition signal; The signal preprocessing module is used to screen the key factors of the wavelet packet threshold method based on a high signal-to-noise ratio, wherein the key factors include the optimal wavelet packet basis function, the number of decomposition layers, the threshold function and the threshold rule; the leakage acoustic wave acquisition signal is decomposed at multiple scales based on the various factors in the wavelet packet threshold method, and the last layer node coefficients are screened and reconstructed using the correlation coefficient method to obtain the preprocessed signal; The algorithm optimization module is used to use the minimum envelope entropy as the fitness function and the enhanced whale algorithm to optimize the variational mode decomposition method to obtain the optimal parameter combination; The denoising module is used to construct an optimized variational mode decomposition using the optimal parameter combination, decompose the preprocessed signal to obtain the intrinsic mode function components; and screen and reconstruct the intrinsic mode function components using the correlation coefficient method to achieve denoising of leakage information.

[0043] Embodiment 3 This embodiment provides a computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the steps in the leakage sound wave denoising method based on wavelet packet threshold and VMD as described in the first embodiment above are implemented.

[0044] Embodiment 4 This embodiment provides a computer device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the steps in the leakage sound wave denoising method based on wavelet packet threshold and VMD as described in the first embodiment are implemented.

[0045] The steps or modules involved in the above embodiments 2 to 4 correspond to those in embodiment 1. For the specific implementation, please refer to the relevant description of embodiment 1. The term "computer-readable storage medium" should be understood as a single medium or multiple media including one or more instruction sets; it should also be understood to include any medium that can store, encode or carry an instruction set for execution by a processor and enable the processor to execute any method in the present invention.

[0046] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A leakage sound wave denoising method based on wavelet packet threshold and VMD, characterized in that: include: Acquire leakage acoustic wave acquisition signals; Based on high signal-to-noise ratio, the key factors of the wavelet packet threshold method are selected, and the key factors include the optimal wavelet packet basis function, the number of decomposition layers, the threshold function and the threshold rule; based on the various factors in the wavelet packet threshold method, the leakage acoustic wave acquisition signal is multi-scale decomposed, and the correlation coefficient method is used to select and reconstruct the node coefficients of the last layer to obtain the preprocessed signal; Taking the minimum envelope entropy as the fitness function, the enhanced whale algorithm is used to optimize the variational mode decomposition method to obtain the optimal parameter combination; The optimal parameter combination is used to construct an optimized variational mode decomposition, and the preprocessed signal is decomposed to obtain the intrinsic mode function components; The correlation coefficient method is used to screen and reconstruct the intrinsic mode function components to achieve denoising of leakage information.

2. The method for reducing leakage sound wave noise based on wavelet packet threshold and VMD as claimed in claim 1, characterized in that: The composition parameters of the wavelet packet threshold method based on high signal-to-noise ratio screening specifically include: Obtain the original signal of the sample, and perform three rounds of wavelet decomposition on the original signal of the sample based on different wavelet packet basis functions, decomposition levels, threshold functions and threshold rules according to the control variable idea, and select the optimal parameters based on the maximum signal-to-noise ratio of the reconstructed signal and the original signal of the sample; Among them, in the first round, the number of decomposition layers, threshold functions and threshold rules are fixed, and the initial wavelet packet basis functions are screened out based on a variety of wavelet packet basis functions; the second round fixes the threshold functions and threshold rules, and the optimal wavelet packet basis functions and the optimal number of decomposition layers are screened out based on the initial wavelet packet basis functions and different numbers of decomposition layers; the third round fixes the optimal wavelet packet basis functions and the optimal number of decomposition layers, and the optimal threshold functions and the optimal threshold rules are screened out based on a variety of threshold functions and a variety of threshold rules.

3. The method for reducing leakage sound wave noise based on wavelet packet threshold and VMD as claimed in claim 1, characterized in that: The method of performing multi-scale decomposition on the leakage acoustic wave acquisition signal based on various parameters in the wavelet packet threshold method, and using the correlation coefficient method to screen and reconstruct the last layer of node coefficients to obtain the preprocessed signal specifically includes: The optimal wavelet packet basis function, decomposition layer number, threshold function and threshold rule are used to perform multi-scale decomposition on the leakage acoustic wave acquisition signal, and the Pearson correlation coefficient between the node coefficients of the last layer and the leakage acoustic wave acquisition signal is calculated; The nodes whose correlation is greater than or equal to a preset value are reconstructed to obtain a preprocessed signal.

4. The method for reducing leakage sound wave noise based on wavelet packet threshold and VMD as claimed in claim 1, characterized in that: The method uses the minimum envelope entropy as the fitness function and adopts the enhanced whale algorithm to optimize the variational mode decomposition method to obtain the optimal parameter combination, which specifically includes: Initialize the enhanced whale algorithm and determine the initial whale individual based on the minimum envelope entropy; each initial whale individual represents a set of parameter combinations ; The enhanced whale algorithm is used to update the initial individual whale positions; Variational modal decomposition is performed on each whale individual at each updated position, the sum of the envelope entropy of the modal components after decomposition of each whale individual is calculated, and the whale individual with the smallest envelope entropy is recorded; when the maximum number of iterations is reached, the whale individual with the smallest sum of envelope entropy is output as the optimal whale individual; the optimal whale individual represents the optimal parameter combination, including the optimal number of decomposition layers and the penalty factor.

5. The method for reducing leakage sound wave noise based on wavelet packet threshold and VMD as claimed in claim 4, characterized in that: Initializing the enhanced whale algorithm and determining the initial whale individual based on the minimum envelope entropy includes: Set the number of whale populations, generate individual whales, and optimize the target parameters Range restriction; the individual whales represent a group Parameter combination, where is the number of decomposition layers, is the penalty factor; Initialize the parameters of the variational mode decomposition method, perform variational mode decomposition on each whale individual, calculate the sum of the envelope entropies of the modal components of the decomposed signal, and record and save the whale individual with the minimum envelope entropy as the initial whale individual.

6. The method for reducing leakage sound wave noise based on wavelet packet threshold and VMD as claimed in claim 4, characterized in that: The enhanced whale algorithm is used to update the initial individual whale positions, specifically including: Generate a random number p, and determine whether the random number p is greater than or equal to 0.5; If yes, a weight factor is introduced and a spiral update mechanism is selected to update the individual positions of the whales; If not, continue to judge Is it greater than 1? If so, a random search mechanism is selected to update the individual positions of the whales; if not, a differential mutation perturbation factor is introduced, and a shrinkage and encirclement mechanism is selected to update the individual positions of the whales.

7. The method for reducing leakage sound wave noise based on wavelet packet threshold and VMD as claimed in claim 1, characterized in that: The optimal parameter combination is used to construct an optimized variational mode decomposition, and the preprocessed signal is decomposed to obtain the intrinsic mode function component; The correlation coefficient method is used to screen and reconstruct the intrinsic mode function components to achieve noise reduction of leakage information, including: The optimal parameter combination is input into the variational mode decomposition algorithm to decompose the multi-scale acoustic wave signal, and the Pearson correlation coefficient between each intrinsic mode function component and the leakage acoustic wave acquisition signal is calculated; The modal function components whose correlation is greater than or equal to a preset value are reconstructed to obtain the noise-reduced leakage sound wave signal.

8. A leakage sound wave noise reduction system based on wavelet packet threshold and VMD, characterized in that: include: The original signal acquisition module is used to obtain the leakage sound wave acquisition signal; The signal preprocessing module is used to screen the key factors of the wavelet packet threshold method based on a high signal-to-noise ratio, wherein the key factors include the optimal wavelet packet basis function, the number of decomposition layers, the threshold function and the threshold rule; the leakage acoustic wave acquisition signal is decomposed at multiple scales based on the various factors in the wavelet packet threshold method, and the last layer node coefficients are screened and reconstructed using the correlation coefficient method to obtain the preprocessed signal; The algorithm optimization module is used to use the minimum envelope entropy as the fitness function and the enhanced whale algorithm to optimize the variational mode decomposition method to obtain the optimal parameter combination; A noise reduction module, used for constructing an optimized variational mode decomposition using the optimal parameter combination, decomposing the preprocessed signal, and obtaining an intrinsic mode function component; The correlation coefficient method is used to screen and reconstruct the intrinsic mode function components to achieve denoising of leakage information.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the steps in a leakage sound wave denoising method based on wavelet packet threshold and VMD as described in any one of claims 1 to 7 are implemented.

10. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the steps in the leakage sound wave denoising method based on wavelet packet threshold and VMD as described in any one of claims 1 to 7 are implemented.

Images