Denoising method and device for wire vibration acceleration signal, electronic equipment and storage medium

Abstract

The invention discloses a de-noising method and device for a wire vibration acceleration signal, electronic equipment and a storage medium, and belongs to the technical field of signal de-noising, and the method comprises the steps: collecting a vibration acceleration signal which is not de-noised, and then cyclically executing decomposition and reconstruction; in each round of iteration, determining a modal decomposition order range according to a current signal, and selecting an optimal order with the maximum weighted kurtosis; selecting an optimal penalty factor with the maximum harmonic coefficient from a preset penalty factor range; and decomposing and generating candidate modal components based on the two parameters, and screening effective components to reconstruct a new signal. And judging a convergence state by calculating an energy difference value between the current reconstructed signal and the last reconstructed signal: if the difference value is smaller than a threshold value, outputting a final de-noised signal, otherwise, taking the current reconstructed signal as the next input to continue iteration. Therefore, by implementing the method and the device, the problem that de-noising of the wire vibration acceleration signal is not accurate enough in the prior art can be solved.

Description

Technical Field

[0001] The present invention relates to the technical field of signal denoising, and in particular to a denoising method, device, electronic equipment and storage medium for a wire vibration acceleration signal. Background Art

[0002] Conductor vibration acceleration signals are important physical quantities used to monitor conductor operating conditions in power systems. During operation, conductors are affected by various factors, such as wind, external vibrations, and equipment failures, generating vibrations of varying frequencies and amplitudes. However, due to factors such as environmental noise and equipment interference, collected conductor vibration acceleration signals often contain significant noise. This noise can severely impact the signal's validity, making it difficult to accurately extract valuable physical features from the signal, thus hindering subsequent fault diagnosis and equipment health monitoring.

[0003] However, existing technologies often lack accuracy in the process of denoising wire vibration acceleration signals. Existing methods usually use variational modal decomposition for signal denoising, but often use fixed modal decomposition orders and penalty factors. The modal decomposition order and penalty factor are two important parameters in variational modal decomposition. The modal decomposition order determines the decomposition accuracy of the signal. A lower order may lead to over-simplification of the signal, while a higher order may lead to over-refined decomposition and ineffective noise removal. The penalty factor is used to control the smoothness and frequency range of the modal decomposition. A penalty factor that is too large may cause the signal to be over-smoothed, while a penalty factor that is too small may not be able to effectively remove noise. Therefore, when processing signals, the denoising effect of the variational modal decomposition method may show over-denoising or under-denoising. Summary of the Invention

[0004] The embodiments of the present invention provide a method, device, electronic device and storage medium for denoising a wire vibration acceleration signal, which can solve the problem of inaccurate denoising of a wire vibration acceleration signal in the prior art.

[0005] An embodiment of the present invention provides a method for denoising a wire vibration acceleration signal, comprising:

[0006] Obtaining the non-denoised vibration acceleration signal of the conductor;

[0007] Repeat the variational mode decomposition and reconstruction operation until the energy difference between the vibration acceleration reconstructed signal of the current round and the vibration acceleration reconstructed signal of the previous round is less than a preset convergence threshold, thereby generating a denoised vibration acceleration signal;

[0008] The variational mode decomposition and reconstruction operation includes:

[0009] Determining the range of the modal decomposition order corresponding to the current round according to the vibration acceleration signal of the current round; wherein the vibration acceleration signal of the initial round is the non-denoised vibration acceleration signal;

[0010] From the range of modal decomposition orders corresponding to the current round, the modal decomposition order with the largest weighted kurtosis index is selected as the optimal modal decomposition order corresponding to the current round;

[0011] Based on the optimal modal decomposition order corresponding to the current round, the penalty factor with the largest harmonic coefficient is selected from the preset penalty factor range as the optimal penalty factor corresponding to the current round;

[0012] Decomposing the vibration acceleration signal of the current round according to the optimal modal decomposition order and the optimal penalty factor to generate a number of candidate modal components;

[0013] Screen each candidate modal component to generate the vibration acceleration reconstruction signal of the current round;

[0014] Determine whether the energy difference between the current round of vibration acceleration reconstruction signal and the previous round of vibration acceleration reconstruction signal is less than a preset convergence threshold. If so, use the current round of vibration acceleration reconstruction signal as the denoised vibration acceleration signal; if not, update the current round of vibration acceleration reconstruction signal to the next round of vibration acceleration signal.

[0015] Furthermore, the step of obtaining the non-denoised vibration acceleration signal of the conductor includes:

[0016] Obtaining the original vibration acceleration signal of the conductor;

[0017] Cleaning the original vibration acceleration signal of the conductor to generate a cleaned vibration acceleration signal;

[0018] Performing coordinate transformation processing on the vibration acceleration signal after cleaning to generate a vibration acceleration signal in the vertical direction, and removing the gravity acceleration signal component in the vibration acceleration signal in the vertical direction to generate an effective vibration acceleration signal in the vertical direction;

[0019] An effective vibration acceleration signal within a preset time period is selected from the effective vibration acceleration signals in the vertical direction as the non-denoised vibration acceleration signal of the wire.

[0020] Furthermore, the determining of the range of the modal decomposition order corresponding to the current round according to the vibration acceleration signal of the current round includes:

[0021] Calculate the number of narrowband interference components in the vibration acceleration signal of the current round according to the classical threshold method, and use the number of narrowband interference components as the minimum value of the modal decomposition order corresponding to the current round;

[0022] According to the smaller threshold method, the number of pulse interferences and the total number of frequency band interference components in the vibration acceleration signal of the current round are calculated;

[0023] Generate a maximum value of the modal decomposition order corresponding to the current round according to the number of pulse interferences, the total number of frequency band interference components, and a preset modal margin;

[0024] The range of the modal decomposition order corresponding to the current round is determined according to the minimum value of the modal decomposition order corresponding to the current round and the maximum value of the modal decomposition order corresponding to the current round.

[0025] Furthermore, the step of selecting the modal decomposition order with the largest weighted kurtosis index from the range of modal decomposition orders corresponding to the current round as the optimal modal decomposition order corresponding to the current round includes:

[0026] Decomposing the vibration acceleration signal of the current round in sequence according to each modal decomposition order within the range of the modal decomposition order corresponding to the current round and a preset first penalty factor to generate a modal component group corresponding to each modal decomposition order;

[0027] For each modal decomposition order, the weighted kurtosis index corresponding to the current modal decomposition order is calculated based on the modal component group corresponding to the current modal decomposition order and the vibration acceleration signal of the current round;

[0028] According to the weighted kurtosis index corresponding to each modal decomposition order, the modal decomposition order with the largest weighted kurtosis index is selected as the optimal modal decomposition order corresponding to the current round.

[0029] Furthermore, the step of selecting the modal decomposition order with the largest weighted kurtosis index from the range of modal decomposition orders corresponding to the current round as the optimal modal decomposition order corresponding to the current round includes:

[0030] Decomposing the vibration acceleration signal of the current round in sequence according to each modal decomposition order within the range of the modal decomposition order corresponding to the current round and a preset first penalty factor to generate a modal component group corresponding to each modal decomposition order;

[0031] For each modal decomposition order, the weighted kurtosis index corresponding to the current modal decomposition order is calculated based on the modal component group corresponding to the current modal decomposition order and the vibration acceleration signal of the current round;

[0032] According to the weighted kurtosis index corresponding to each modal decomposition order, the modal decomposition order with the largest weighted kurtosis index is selected as the optimal modal decomposition order corresponding to the current round.

[0033] Furthermore, based on the optimal modal decomposition order corresponding to the current round, the penalty factor with the largest harmonic coefficient is selected from a preset penalty factor range as the optimal penalty factor corresponding to the current round, including:

[0034] According to the optimal modal decomposition order corresponding to the current round and each penalty factor in the preset penalty factor range, the vibration acceleration signal of the current round is decomposed in sequence to generate a modal component group corresponding to each penalty factor;

[0035] For each penalty factor, Fourier transform is performed on the modal component group corresponding to the current penalty factor to generate the spectrum amplitude group corresponding to the current penalty factor;

[0036] Extracting the spectrum amplitude at a preset excitation frequency point and the spectrum amplitudes of a preset number of frequency points adjacent to the excitation frequency point from the spectrum amplitude group corresponding to the current penalty factor as the target spectrum amplitude group corresponding to the current penalty factor;

[0037] Calculate and generate harmonic coefficients corresponding to the current penalty factor based on the spectrum amplitude group corresponding to the current penalty factor and the target spectrum amplitude group corresponding to the current penalty factor;

[0038] According to the harmonic coefficients corresponding to each penalty factor, the penalty factor with the largest harmonic coefficient is selected as the optimal penalty factor corresponding to the current round.

[0039] Furthermore, the screening of the candidate modal components to generate the vibration acceleration reconstruction signal of the current round includes:

[0040] Calculate and generate the Pearson correlation coefficient between each candidate modal component and the vibration acceleration signal of the current round based on each candidate modal component and the vibration acceleration signal of the current round;

[0041] The candidate modal component with a Pearson correlation coefficient greater than a preset coefficient threshold is used as the target modal component of the current round;

[0042] The target modal components of the current round are combined to generate the vibration acceleration reconstruction signal of the current round.

[0043] Furthermore, the energy difference is calculated by the following formula:

[0044]

[0045] Where ΔE kis the energy difference between the vibration acceleration reconstruction signal of the kth round and the k-1th round; N is the total number of sampling points; x k (i) is the i-th sampling point value of the vibration acceleration reconstruction signal of the k-th round; x k-1 (i) is the i-th sampling point value of the vibration acceleration reconstruction signal of the k-1th round.

[0046] Based on the above method embodiments, the present invention provides corresponding device embodiments.

[0047] An embodiment of the present invention provides a denoising device for a wire vibration acceleration signal, comprising: a data acquisition module and a variational mode decomposition and reconstruction module.

[0048] The variational modal decomposition and reconstruction module is used to repeatedly perform the variational modal decomposition and reconstruction operation until the energy difference between the vibration acceleration reconstructed signal of the current round and the vibration acceleration reconstructed signal of the previous round is less than a preset convergence threshold, thereby generating a denoised vibration acceleration signal;

[0049] The variational mode decomposition and reconstruction operation includes:

[0050] Determining the range of the modal decomposition order corresponding to the current round according to the vibration acceleration signal of the current round; wherein the vibration acceleration signal of the initial round is the non-denoised vibration acceleration signal;

[0051] From the range of modal decomposition orders corresponding to the current round, the modal decomposition order with the largest weighted kurtosis index is selected as the optimal modal decomposition order corresponding to the current round;

[0052] Based on the optimal modal decomposition order corresponding to the current round, the penalty factor with the largest harmonic coefficient is selected from the preset penalty factor range as the optimal penalty factor corresponding to the current round;

[0053] Decomposing the vibration acceleration signal of the current round according to the optimal modal decomposition order and the optimal penalty factor to generate a number of candidate modal components;

[0054] Screen each candidate modal component to generate the vibration acceleration reconstruction signal of the current round;

[0055] Determine whether the energy difference between the current round of vibration acceleration reconstruction signal and the previous round of vibration acceleration reconstruction signal is less than a preset convergence threshold. If so, use the current round of vibration acceleration reconstruction signal as the denoised vibration acceleration signal; if not, update the current round of vibration acceleration reconstruction signal to the next round of vibration acceleration signal.

[0056] Based on the above method embodiment, the present invention provides a corresponding electronic device embodiment.

[0057] An embodiment of the present invention provides an electronic device, comprising a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor. When the processor executes the computer program, the method for denoising a wire vibration acceleration signal as described in any one of the above-mentioned method embodiments is implemented.

[0058] Based on the above method embodiment, the present invention provides a corresponding storage medium embodiment.

[0059] An embodiment of the present invention provides a storage medium having a computer program stored thereon, wherein when the computer program is executed, the device where the storage medium is located is controlled to execute the method for denoising a wire vibration acceleration signal as described in any one of the above method embodiments.

[0060] Compared with the prior art, the present invention has the following beneficial effects:

[0061] The embodiment of the present invention provides a denoising method, device, electronic device and storage medium for a wire vibration acceleration signal. The method obtains the undenoised vibration acceleration signal of the wire, repeatedly performs the variational modal decomposition and reconstruction operation until the energy difference between the vibration acceleration reconstruction signal of the current round and the signal of the previous round is less than a preset convergence threshold, and finally generates a denoised signal; in each round, first determine the candidate range of the modal decomposition order based on the current vibration acceleration signal, and then select the optimal modal order based on the weighted kurtosis index; under this order, further select the penalty factor with the largest harmonic coefficient from the penalty factor candidate set as the optimal value to complete the parameter adaptive selection. Subsequently, the signal is decomposed using the current optimal parameters to obtain candidate modal components; by screening the candidate modal components, the vibration acceleration reconstruction signal of the current round is constructed, and it is determined whether it converges.

[0062] Unlike traditional variational modal decomposition denoising methods that usually adopt fixed parameters and complete decomposition in one go, the technical solution of the present invention fully considers the dynamic nature of signal characteristics as the rounds change. In each round of decomposition, weighted kurtosis and harmonic coefficients are used as evaluation indicators to select the modal decomposition order and penalty factor respectively, effectively improving the degree of matching between the decomposition result and the actual signal. At the same time, the candidate modal components are screened to further improve the effectiveness of modal reconstruction. Through multiple rounds of iterative decomposition and reconstruction operations, and using energy difference as the basis for convergence, not only can the problem of excessive or insufficient denoising caused by improper parameter setting in traditional methods be avoided, but the signal can also be gradually purified, and the ability to extract effective signal components is enhanced, thereby significantly improving the accuracy of wire vibration acceleration signal denoising. BRIEF DESCRIPTION OF THE DRAWINGS

[0063] Figure 1 The figure is a flow chart of a method for denoising a wire vibration acceleration signal provided by one embodiment of the present invention.

[0064] Figure 2 It is a flowchart of a variational modal decomposition and reconstruction operation provided by one embodiment of the present invention.

[0065] Figure 3 The figure is a schematic structural diagram of a device for denoising a wire vibration acceleration signal provided by one embodiment of the present invention. DETAILED DESCRIPTION

[0066] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present invention without making any creative efforts shall fall within the scope of protection of the present invention.

[0067] like Figure 1 As shown, in order to solve the problem of inaccurate denoising of wire vibration acceleration signals in the prior art, an embodiment of the present invention provides a denoising method for wire vibration acceleration signals, which includes at least the following steps:

[0068] Step S1, obtaining a non-denoised vibration acceleration signal of the conductor;

[0069] In a preferred embodiment, obtaining the un-denoised vibration acceleration signal of the conductor includes:

[0070] Obtaining the original vibration acceleration signal of the conductor;

[0071] Cleaning the original vibration acceleration signal of the conductor to generate a cleaned vibration acceleration signal;

[0072] Performing coordinate transformation processing on the vibration acceleration signal after cleaning to generate a vibration acceleration signal in the vertical direction, and removing the gravity acceleration signal component in the vibration acceleration signal in the vertical direction to generate an effective vibration acceleration signal in the vertical direction;

[0073] An effective vibration acceleration signal within a preset time period is selected from the effective vibration acceleration signals in the vertical direction as the non-denoised vibration acceleration signal of the wire.

[0074] Specifically, to obtain the conductor's original vibration acceleration signal, a triaxial accelerometer mounted on its surface collects real-time time-series data containing three-dimensional vibration information. The original signal is then cleaned, with transient impact interference eliminated through sliding window outlier detection and high-frequency electromagnetic noise removed using adaptive threshold filtering. This generates a smooth and continuous post-cleaning vibration acceleration signal. The cleaned three-dimensional vibration signal then undergoes coordinate transformation. Based on the relationship between the conductor's suspension orientation and spatial position, vector projection calculations are used to transform the original coordinate system into a local coordinate system based on the vertical direction, extracting the vertical axial vibration component. Gravity acceleration compensation is performed on the vertical signal. The mean offset of the static gravity component is calculated using a sliding time window. Dynamic baseline correction is then used to remove the constant influence of gravity acceleration, ultimately isolating the effective acceleration signal reflecting the conductor's true vibration characteristics. After obtaining the effective vibration signal, continuous time periods with stable vibration amplitude and uniform frequency distribution are selected based on the duration characteristics of the conductor's vibration stability, serving as the undenoised vibration acceleration signal for the final analysis. This step ensures that the selected time period is free of external sudden interference and meets the vibration pattern continuity requirements.

[0075] Step S2: Repeat the variational mode decomposition and reconstruction operation until the energy difference between the vibration acceleration reconstructed signal of the current round and the vibration acceleration reconstructed signal of the previous round is less than a preset convergence threshold, thereby generating a denoised vibration acceleration signal;

[0076] like Figure 2 As shown, in a preferred embodiment, the variational mode decomposition and reconstruction operation includes:

[0077] Step S2.1, determining the range of the modal decomposition order corresponding to the current round based on the vibration acceleration signal of the current round; wherein the vibration acceleration signal of the initial round is the non-denoised vibration acceleration signal;

[0078] In a preferred embodiment, determining the range of the modal decomposition order corresponding to the current round based on the vibration acceleration signal of the current round includes:

[0079] Calculate the number of narrowband interference components in the vibration acceleration signal of the current round according to the classical threshold method, and use the number of narrowband interference components as the minimum value of the modal decomposition order corresponding to the current round;

[0080] According to the smaller threshold method, the number of pulse interferences and the total number of frequency band interference components in the vibration acceleration signal of the current round are calculated;

[0081] Generate a maximum value of the modal decomposition order corresponding to the current round according to the number of pulse interferences, the total number of frequency band interference components, and a preset modal margin;

[0082] The range of the modal decomposition order corresponding to the current round is determined according to the minimum value of the modal decomposition order corresponding to the current round and the maximum value of the modal decomposition order corresponding to the current round.

[0083] Specifically, the modal decomposition order range is determined by combining spectrum analysis with interference component quantification. First, the vibration signal is subjected to spectrum analysis. The amplitude threshold is set to 8%-12% of the fundamental amplitude using the classical threshold method. The number of discrete narrowband interference peaks exceeding the threshold (e.g., 2 power frequency interferences and 1 mechanical resonance peak are detected) is counted and used as the minimum value K of the modal decomposition order. min On this basis, the time domain pulse detection and frequency band energy scanning are performed by using a smaller threshold method (taking 3%-5% of the fundamental amplitude), and the transient pulse events exceeding the time domain kurtosis threshold (such as detecting 3 lightning interferences) and the number of broadband interference areas identified by wavelet packet decomposition (such as 2 wind noise bands) are counted to obtain the total number of interference components N. Finally, combined with the preset modal margin Δ (usually 2-3), the formula K is used. max =N+Δ determines the upper limit of the decomposition order (for example, when N=5, K max =7), thus forming the reasonable order range of the modal decomposition of this round [K min ,K max ].

[0084] This range determination method can avoid modal aliasing caused by insufficient order (such as K min When the order is too high, it can also prevent the introduction of noise pseudo components (such as K max If it is too large, meaningless high-frequency modes will be generated), ensuring that the decomposition dimension matches the actual vibration characteristics.

[0085] Step S2.2: From the range of modal decomposition orders corresponding to the current round, select the modal decomposition order with the largest weighted kurtosis index as the optimal modal decomposition order corresponding to the current round;

[0086] In a preferred embodiment, selecting the modal decomposition order with the largest weighted kurtosis index from the range of modal decomposition orders corresponding to the current round as the optimal modal decomposition order corresponding to the current round includes:

[0087] Decomposing the vibration acceleration signal of the current round in sequence according to each modal decomposition order within the range of the modal decomposition order corresponding to the current round and a preset first penalty factor to generate a modal component group corresponding to each modal decomposition order;

[0088] For each modal decomposition order, the weighted kurtosis index corresponding to the current modal decomposition order is calculated based on the modal component group corresponding to the current modal decomposition order and the vibration acceleration signal of the current round;

[0089] According to the weighted kurtosis index corresponding to each modal decomposition order, the modal decomposition order with the largest weighted kurtosis index is selected as the optimal modal decomposition order corresponding to the current round.

[0090] Specifically, the weighted kurtosis index is calculated by the following formula:

[0091]

[0092] Among them, KEI k is the weighted kurtosis index of the kth round; ρ is the preset adjustment factor; ΔE k The energy difference between the vibration acceleration reconstruction signal of the kth round and the k-1th round; KI k is the kurtosis index of the kth round;

[0093]

[0094] Among them, KI k is the kurtosis index of the kth round; K k is the modal decomposition order of the kth round; u j，k (i) is the value of the jth modal component at the i-th sampling point in the k-th round; μ j，k is the mean of the jth modal component in the kth round; σ j，k is the variance of the jth modal component in the kth round;

[0095]

[0096] Where ΔE k is the energy difference between the vibration acceleration reconstruction signal of the kth round and the k-1th round; N is the total number of sampling points; x k (i) is the i-th sampling point value of the vibration acceleration reconstruction signal of the k-th round; x k-1 (i) is the i-th sampling point value of the vibration acceleration reconstruction signal of the k-1th round.

[0097] Specifically, according to the modal decomposition order range of the current round, the vibration acceleration signal is sequentially subjected to multi-dimensional decomposition and feature evaluation. For each candidate value within the order range, a preset basic penalty factor is used to perform variational modal decomposition to obtain a decomposition result group containing different numbers of modal components. For each decomposition result group, the effective vibration feature weight of each component is determined by calculating the time domain correlation and energy proportion of each modal component with the original signal. On this basis, the kurtosis statistical characteristics of each component (reflecting the signal impact intensity) and the corresponding weight value are combined to generate a weighted kurtosis index that characterizes the effectiveness of the current decomposition order. By horizontally comparing the index values ​​of all candidate orders, the decomposition order that maximizes the weighted kurtosis index is selected as the optimal solution.

[0098] This step ensures that the decomposition process can not only fully capture the main frequency components of the conductor vibration (such as low-frequency dancing and high-frequency breeze vibration), but also avoid the introduction of noise pseudo-components due to too high an order or modal aliasing due to insufficient order, thereby providing an accurate decomposition dimension basis for subsequent penalty factor optimization.

[0099] Step S2.3: Based on the optimal modal decomposition order corresponding to the current round, select the penalty factor with the largest harmonic coefficient from the preset penalty factor range as the optimal penalty factor corresponding to the current round;

[0100] In a preferred embodiment, the method of selecting the penalty factor with the largest harmonic coefficient from a preset penalty factor range based on the optimal modal decomposition order corresponding to the current round as the optimal penalty factor corresponding to the current round includes:

[0101] According to the optimal modal decomposition order corresponding to the current round and each penalty factor in the preset penalty factor range, the vibration acceleration signal of the current round is decomposed in sequence to generate a modal component group corresponding to each penalty factor;

[0102] For each penalty factor, Fourier transform is performed on the modal component group corresponding to the current penalty factor to generate the spectrum amplitude group corresponding to the current penalty factor;

[0103] Extracting the spectrum amplitude at a preset excitation frequency point and the spectrum amplitudes of a preset number of frequency points adjacent to the excitation frequency point from the spectrum amplitude group corresponding to the current penalty factor as the target spectrum amplitude group corresponding to the current penalty factor;

[0104] Calculate and generate harmonic coefficients corresponding to the current penalty factor based on the spectrum amplitude group corresponding to the current penalty factor and the target spectrum amplitude group corresponding to the current penalty factor;

[0105] According to the harmonic coefficients corresponding to each penalty factor, the penalty factor with the largest harmonic coefficient is selected as the optimal penalty factor corresponding to the current round.

[0106] Exemplarily, when the preset number is 2, the harmonic coefficient corresponding to the current penalty factor is calculated and generated by the following formula:

[0107]

[0108] Among them, β is the harmonic coefficient corresponding to the current penalty factor; I f is the spectrum amplitude at the excitation frequency; I f+1 is the spectrum amplitude at the sampling point after the excitation frequency, that is, the adjacent frequency point to the right of the excitation frequency; I f-1 is the spectrum amplitude at the sampling point before the excitation frequency, that is, the adjacent frequency point on the left side of the excitation frequency; M is the total number of frequency points in the current spectrum amplitude group, that is, the total number of spectrum sampling points; I m is the amplitude at the mth frequency point in the spectrum amplitude group.

[0109] Specifically, in this step, a variational modal decomposition operation is performed based on the optimal modal decomposition order corresponding to the current round, mainly to improve the overall computational efficiency while ensuring the decomposition accuracy. Specifically, after determining the modal decomposition order, the system traverses the preset penalty factor range (e.g., 1000-5000), performs variational modal decomposition on the vibration acceleration signal of the current round, and generates the corresponding modal component group. Subsequently, a Fourier transform is performed on each group of decomposition results to extract the amplitude of the preset excitation frequency point and the five frequency points before and after it to form the target spectrum amplitude group, and calculates the proportion of the squared sum of the amplitude of this group in the squared sum of the amplitude of the entire frequency band as the harmonic coefficient corresponding to the current penalty factor. By comparing the harmonic coefficient values ​​of all penalty factors, the penalty factor that can maximize the harmonic coefficient is selected as the optimal penalty factor for the current round (e.g., when α=3000, the coefficient is 0.72, and when α=3500, it is 0.81, so α=3500 is preferred). This optimization mechanism effectively ensures that the inherent vibration characteristics of the wire are retained during the decomposition process by quantifying the energy concentration near the excitation frequency, while reducing the impact of non-characteristic noise on the subsequent reconstruction accuracy.

[0110] Step S2.4: Decompose the vibration acceleration signal of the current round according to the optimal modal decomposition order and the optimal penalty factor to generate a number of candidate modal components;

[0111] Specifically, based on the determined optimal modal decomposition order and penalty factor, the specific process of performing variational modal decomposition to generate candidate modal components is as follows: the vibration signal of the current round is input into the variational modal decomposition algorithm. The algorithm decomposes the signal into multiple modal components with different center frequencies through iterative solution according to the set decomposition order and penalty factor. During the decomposition process, preset convergence conditions (such as a residual energy ratio threshold of 0.1% or a maximum number of iterations of 200) are adopted to ensure the stability of the frequency bandwidth of each modal component. The generated candidate modal components are arranged from low to high frequency, where low-frequency components mainly contain large-scale vibration characteristics such as wire dancing, while high-frequency components contain breeze vibration and noise components. The time domain waveform and frequency domain characteristics of each component are compared and verified with the original signal, and components that obviously deviate from physical laws (such as abnormal modes with drastic instantaneous frequency jumps or energy accounts for less than 1%) are eliminated. Finally, a group of candidate modal components with clear vibration characteristics is retained.

[0112] Step S2.5: Screen each candidate modal component to generate a vibration acceleration reconstruction signal for the current round;

[0113] In a preferred embodiment, screening the candidate modal components to generate the vibration acceleration reconstruction signal of the current round includes:

[0114] Calculate and generate the Pearson correlation coefficient between each candidate modal component and the vibration acceleration signal of the current round based on each candidate modal component and the vibration acceleration signal of the current round;

[0115] The candidate modal component with a Pearson correlation coefficient greater than a preset coefficient threshold is used as the target modal component of the current round;

[0116] The target modal components of the current round are combined to generate the vibration acceleration reconstruction signal of the current round.

[0117] Specifically, effective modal component screening is achieved through quantitative correlation analysis. First, the Pearson correlation coefficient is calculated for each candidate modal component and the original vibration signal. This coefficient reflects the degree of linear correlation between the component and the original signal. This calculation utilizes a point-to-point matching method for time-domain data, with a correlation coefficient threshold typically set between 0.3 and 0.5 (the specific value varies depending on the conductor type, such as 0.35 for ACSR). Modal components with correlation coefficients exceeding the threshold are identified as valid vibration components (e.g., a component with a correlation coefficient of 0.42 is retained), while low-correlation noise components (e.g., a component with a correlation coefficient of 0.25) are eliminated. The filtered target modal components are then superimposed on their time-domain waveforms to reconstruct a denoised signal containing the primary vibration characteristics. The reconstruction process maintains phase consistency between the components to avoid waveform distortion. By quantifying signal correlation, this screening mechanism effectively removes approximately 65% ​​of random noise energy while retaining the characteristic components of conductor galloping (low-frequency 0.1-3Hz) and breeze vibration (high-frequency 10-40Hz), improving the signal-to-noise ratio of the reconstructed signal by 8-12dB.

[0118] Step S2.6, determine whether the energy difference between the vibration acceleration reconstruction signal of the current round and the vibration acceleration reconstruction signal of the previous round is less than the preset convergence threshold. If so, the vibration acceleration reconstruction signal of the current round is used as the denoised vibration acceleration signal; if not, the vibration acceleration reconstruction signal of the current round is updated to the vibration acceleration signal of the next round.

[0119] In a preferred embodiment, the energy difference is calculated by the following formula:

[0120]

[0121] Where ΔE k is the energy difference between the vibration acceleration reconstruction signal of the kth round and the k-1th round; N is the total number of sampling points; x k (i) is the i-th sampling point value of the vibration acceleration reconstruction signal of the k-th round; x k-1 (i) is the i-th sampling point value of the vibration acceleration reconstruction signal of the k-1th round.

[0122] Based on the above method embodiments, the present invention provides corresponding device embodiments.

[0123] like Figure 3 As shown, an embodiment of the present invention provides a denoising device for a wire vibration acceleration signal, comprising: a data acquisition module and a variational mode decomposition and reconstruction module;

[0124] The data acquisition module is used to acquire the non-denoised vibration acceleration signal of the wire;

[0125] The variational modal decomposition and reconstruction module is used to repeatedly perform the variational modal decomposition and reconstruction operation until the energy difference between the vibration acceleration reconstructed signal of the current round and the vibration acceleration reconstructed signal of the previous round is less than a preset convergence threshold, thereby generating a denoised vibration acceleration signal;

[0126] The variational mode decomposition and reconstruction operation includes:

[0127] Determining the range of the modal decomposition order corresponding to the current round according to the vibration acceleration signal of the current round; wherein the vibration acceleration signal of the initial round is the non-denoised vibration acceleration signal;

[0128] From the range of modal decomposition orders corresponding to the current round, the modal decomposition order with the largest weighted kurtosis index is selected as the optimal modal decomposition order corresponding to the current round;

[0129] Based on the optimal modal decomposition order corresponding to the current round, the penalty factor with the largest harmonic coefficient is selected from the preset penalty factor range as the optimal penalty factor corresponding to the current round;

[0130] Decomposing the vibration acceleration signal of the current round according to the optimal modal decomposition order and the optimal penalty factor to generate a number of candidate modal components;

[0131] Screen each candidate modal component to generate the vibration acceleration reconstruction signal of the current round;

[0132] Determine whether the energy difference between the current round of vibration acceleration reconstruction signal and the previous round of vibration acceleration reconstruction signal is less than a preset convergence threshold. If so, use the current round of vibration acceleration reconstruction signal as the denoised vibration acceleration signal; if not, update the current round of vibration acceleration reconstruction signal to the next round of vibration acceleration signal.

[0133] It should be noted that the embodiment of the device described above corresponds to the above-mentioned embodiment of the present invention, and it can implement the denoising method of the wire vibration acceleration signal described in any one of the above-mentioned embodiments of the present invention. In addition, the embodiment of the above-mentioned device is merely schematic, wherein the modules described as separate components may or may not be physically separated, and the components displayed as modules may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the scheme of this embodiment. In addition, in the drawings of the embodiment of the device provided by the present invention, the connection relationship between the modules indicates that there is a communication connection between them, which can be specifically implemented as one or more communication buses or signal lines. A person of ordinary skill in the art can understand and implement it without paying any creative work.

[0134] Based on the above method embodiment of the present invention, a corresponding electronic device embodiment is provided.

[0135] An embodiment of the present invention provides an electronic device, comprising a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor. When the processor executes the computer program, the method for denoising a wire vibration acceleration signal described in any one of the present invention is implemented, or when the processor executes the computer program, the functions of the modules in the above-mentioned device embodiments are implemented.

[0136] Exemplarily, the computer program may be divided into one or more modules, which are stored in the memory and executed by the processor to implement the present invention. The one or more modules may be a series of computer program instruction segments capable of performing specific functions, and the instruction segments are used to describe the execution process of the computer program in the terminal device.

[0137] The terminal device may be a computing device such as a desktop computer, a notebook computer, a PDA, a cloud server, etc. The terminal device may include, but is not limited to, a processor and a memory.

[0138] The processor may be a central processing unit (CPU), other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or any conventional processor, etc. The processor is the control center of the terminal device, connecting various parts of the entire terminal device using various interfaces and lines.

[0139] The memory can be used to store the computer programs and / or modules. The processor implements various functions of the terminal device by running or executing the computer programs and / or modules stored in the memory and calling the data stored in the memory. The memory can mainly include a program storage area and a data storage area, wherein the program storage area can store an operating system, at least one application required for a function, etc.; the data storage area can store data created based on the use of the mobile phone, etc. In addition, the memory can include a high-speed random access memory and can also include a non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a smart memory card (Smart Media Card, SMC), a secure digital (Secure Digital, SD) card, a flash card (Flash Card), at least one disk storage device, a flash memory device, or other volatile solid-state storage device.

[0140] Based on the above method embodiment, the present invention provides a corresponding storage medium embodiment;

[0141] Another embodiment of the present invention provides a storage medium, which includes a stored computer program. When the computer program is executed, the device where the storage medium is located is controlled to execute any one of the above-mentioned wire vibration acceleration signal denoising methods of the present invention.

[0142] The storage medium is a computer-readable storage medium, and the computer program includes computer program code, which may be in source code form, object code form, executable file, or some intermediate form. The computer-readable medium may include any entity or device capable of carrying the computer program code, a recording medium, a USB flash drive, a mobile hard drive, a magnetic disk, an optical disk, a computer memory, a read-only memory (ROM), a random access memory (RAM), an electrical carrier signal, a telecommunications signal, and a software distribution medium.

[0143] In the description of this specification, the reference terms "one embodiment," "some embodiments," "example," "specific example," or "some examples" mean that the specific features, structures, materials, or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. Moreover, the specific features, structures, materials, or characteristics described may be combined in any appropriate manner in any one or more embodiments or examples. In addition, those skilled in the art may combine and integrate different embodiments or examples described in this specification, as well as features of different embodiments or examples, unless they are mutually inconsistent.

[0144] The above is a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principles of the present invention. These improvements and modifications are also considered to be within the scope of protection of the present invention.

Claims

1. A method for denoising a wire vibration acceleration signal, characterized in that: include: Obtaining the non-denoised vibration acceleration signal of the conductor; Repeat the variational mode decomposition and reconstruction operation until the energy difference between the vibration acceleration reconstructed signal of the current round and the vibration acceleration reconstructed signal of the previous round is less than a preset convergence threshold, thereby generating a denoised vibration acceleration signal; The variational mode decomposition and reconstruction operation includes: Determining the range of the modal decomposition order corresponding to the current round according to the vibration acceleration signal of the current round; wherein the vibration acceleration signal of the initial round is the non-denoised vibration acceleration signal; From the range of modal decomposition orders corresponding to the current round, the modal decomposition order with the largest weighted kurtosis index is selected as the optimal modal decomposition order corresponding to the current round; Based on the optimal modal decomposition order corresponding to the current round, the penalty factor with the largest harmonic coefficient is selected from the preset penalty factor range as the optimal penalty factor corresponding to the current round; Decomposing the vibration acceleration signal of the current round according to the optimal modal decomposition order and the optimal penalty factor to generate a number of candidate modal components; Screen each candidate modal component to generate the vibration acceleration reconstruction signal of the current round; Determine whether the energy difference between the current round of vibration acceleration reconstruction signal and the previous round of vibration acceleration reconstruction signal is less than a preset convergence threshold. If so, use the current round of vibration acceleration reconstruction signal as the denoised vibration acceleration signal; if not, update the current round of vibration acceleration reconstruction signal to the next round of vibration acceleration signal.

2. The method for denoising a wire vibration acceleration signal according to claim 1, wherein: The obtaining of the non-denoised vibration acceleration signal of the conductor includes: Obtaining the original vibration acceleration signal of the conductor; Cleaning the original vibration acceleration signal of the conductor to generate a cleaned vibration acceleration signal; Performing coordinate transformation processing on the vibration acceleration signal after cleaning to generate a vibration acceleration signal in the vertical direction, and removing the gravity acceleration signal component in the vibration acceleration signal in the vertical direction to generate an effective vibration acceleration signal in the vertical direction; An effective vibration acceleration signal within a preset time period is selected from the effective vibration acceleration signals in the vertical direction as the non-denoised vibration acceleration signal of the wire.

3. The method for denoising a wire vibration acceleration signal according to claim 2, wherein: Determining the range of the modal decomposition order corresponding to the current round according to the vibration acceleration signal of the current round includes: Calculate the number of narrowband interference components in the vibration acceleration signal of the current round according to the classical threshold method, and use the number of narrowband interference components as the minimum value of the modal decomposition order corresponding to the current round; According to the smaller threshold method, the number of pulse interferences and the total number of frequency band interference components in the vibration acceleration signal of the current round are calculated; Generate a maximum value of the modal decomposition order corresponding to the current round according to the number of pulse interferences, the total number of frequency band interference components, and a preset modal margin; The range of the modal decomposition order corresponding to the current round is determined according to the minimum value of the modal decomposition order corresponding to the current round and the maximum value of the modal decomposition order corresponding to the current round.

4. The method for denoising a wire vibration acceleration signal according to claim 3, wherein: The method of selecting the modal decomposition order with the largest weighted kurtosis index from the range of modal decomposition orders corresponding to the current round as the optimal modal decomposition order corresponding to the current round includes: Decomposing the vibration acceleration signal of the current round in sequence according to each modal decomposition order within the range of the modal decomposition order corresponding to the current round and a preset first penalty factor to generate a modal component group corresponding to each modal decomposition order; For each modal decomposition order, the weighted kurtosis index corresponding to the current modal decomposition order is calculated based on the modal component group corresponding to the current modal decomposition order and the vibration acceleration signal of the current round; According to the weighted kurtosis index corresponding to each modal decomposition order, the modal decomposition order with the largest weighted kurtosis index is selected as the optimal modal decomposition order corresponding to the current round.

5. The method for denoising a wire vibration acceleration signal according to claim 4, wherein: The method of selecting the penalty factor with the largest harmonic coefficient from a preset penalty factor range based on the optimal modal decomposition order corresponding to the current round as the optimal penalty factor corresponding to the current round includes: According to the optimal modal decomposition order corresponding to the current round and each penalty factor in the preset penalty factor range, the vibration acceleration signal of the current round is decomposed in sequence to generate a modal component group corresponding to each penalty factor; For each penalty factor, Fourier transform is performed on the modal component group corresponding to the current penalty factor to generate the spectrum amplitude group corresponding to the current penalty factor; Extracting the spectrum amplitude at a preset excitation frequency point and the spectrum amplitudes of a preset number of frequency points adjacent to the excitation frequency point from the spectrum amplitude group corresponding to the current penalty factor as the target spectrum amplitude group corresponding to the current penalty factor; Calculate and generate harmonic coefficients corresponding to the current penalty factor based on the spectrum amplitude group corresponding to the current penalty factor and the target spectrum amplitude group corresponding to the current penalty factor; According to the harmonic coefficients corresponding to each penalty factor, the penalty factor with the largest harmonic coefficient is selected as the optimal penalty factor corresponding to the current round.

6. The method for denoising a wire vibration acceleration signal according to claim 5, wherein: The step of screening the candidate modal components to generate a vibration acceleration reconstruction signal for the current round includes: Calculate and generate the Pearson correlation coefficient between each candidate modal component and the vibration acceleration signal of the current round based on each candidate modal component and the vibration acceleration signal of the current round; The candidate modal component with a Pearson correlation coefficient greater than a preset coefficient threshold is used as the target modal component of the current round; The target modal components of the current round are combined to generate the vibration acceleration reconstruction signal of the current round.

7. The method for denoising a wire vibration acceleration signal according to claim 6, wherein: The energy difference is calculated using the following formula: Where ΔE k is the energy difference between the vibration acceleration reconstruction signal of the kth round and the k-1th round; N is the total number of sampling points; x k (i) is the i-th sampling point value of the vibration acceleration reconstruction signal of the k-th round; x k-1 (i) is the i-th sampling point value of the vibration acceleration reconstructed signal in the k-1th round.

8. A denoising device for a wire vibration acceleration signal, characterized in that: include: Data acquisition module and variational mode decomposition and reconstruction module; The data acquisition module is used to acquire the non-denoised vibration acceleration signal of the wire; The variational modal decomposition and reconstruction module is used to repeatedly perform the variational modal decomposition and reconstruction operation until the energy difference between the vibration acceleration reconstructed signal of the current round and the vibration acceleration reconstructed signal of the previous round is less than a preset convergence threshold, thereby generating a denoised vibration acceleration signal; The variational mode decomposition and reconstruction operation includes: Determining the range of the modal decomposition order corresponding to the current round according to the vibration acceleration signal of the current round; wherein the vibration acceleration signal of the initial round is the non-denoised vibration acceleration signal; From the range of modal decomposition orders corresponding to the current round, the modal decomposition order with the largest weighted kurtosis index is selected as the optimal modal decomposition order corresponding to the current round; Based on the optimal modal decomposition order corresponding to the current round, the penalty factor with the largest harmonic coefficient is selected from the preset penalty factor range as the optimal penalty factor corresponding to the current round; Decomposing the vibration acceleration signal of the current round according to the optimal modal decomposition order and the optimal penalty factor to generate a number of candidate modal components; Screen each candidate modal component to generate the vibration acceleration reconstruction signal of the current round; Determine whether the energy difference between the current round of vibration acceleration reconstruction signal and the previous round of vibration acceleration reconstruction signal is less than a preset convergence threshold. If so, use the current round of vibration acceleration reconstruction signal as the denoised vibration acceleration signal; if not, update the current round of vibration acceleration reconstruction signal to the next round of vibration acceleration signal.

9. An electronic device, characterized in that: The method comprises a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, wherein when the processor executes the computer program, the method for denoising a wire vibration acceleration signal according to any one of claims 1 to 7 is implemented.

10. A storage medium, characterized in that: The storage medium includes a stored computer program, wherein when the computer program is executed, the device where the storage medium is located is controlled to execute the method for denoising a wire vibration acceleration signal according to any one of claims 1 to 7.

Images