Method and system for eliminating pulsating pressure interference peaks in a pump source pipeline

Abstract

The present invention provides a method and system for eliminating pulsating pressure interference peaks in a pump source pipeline, including: using the zero-phase filtering principle and the total variation filtering method to establish a sparse interference peak filtering model for the logarithmic spectrum of the pulsating pressure amplitude, and then forming an iterative calculation method for the effective signal of the pulsating pressure amplitude spectrum through the principal minimization algorithm. Through iterative calculation, the effective signal of the pulsating pressure can be extracted from the amplitude spectrum containing sparse interference peaks, improving the accuracy of the pulsating pressure signal measurement result. Aiming at the problem of signal distortion in the result obtained by the total variation filtering model, the present invention proposes to perform secondary processing on the pulsating pressure obtained by the filtering model using the moving least squares method, solves the divergence problem of the extraction result of the filtering model, and enhances the reliability of the pulsating pressure filtering algorithm.

Description

Technical Field

[0001] The present invention relates to the technical field of fluid pressure, and in particular to a method and system for eliminating pulsating pressure interference peaks in a pump source pipeline. Background Art

[0002] In daily environments, the sound source pressure wave induced by a pump source and the turbulent pulsating pressure in the fluid boundary layer, as important excitation sources of ship pipelines, often induce pipeline vibrations and reduce their acoustic performance. Therefore, how to accurately obtain the pulsating pressure characteristics of the inner wall of the pipeline driven by the pump source is crucial for understanding the flow-acoustic load excitation characteristics of the pump source and solving pipeline vibration problems.

[0003] Currently, hydrophones or pulsating pressure sensors are often used in engineering to measure the pulsating pressure on the inner wall of the pipeline, and then the excitation characteristics of the flow-acoustic load induced by the pump source are obtained through cross-spectrum analysis and coherence analysis methods. However, since pulsating pressure sensors or hydrophones often need to be installed on the pipeline in an invasive manner, during the measurement process, the pulsating pressure sensor itself will be affected by the vibration of the pipeline structure, resulting in the signal received by the sensor containing not only the pulsating pressure signal of the pipe wall but also the resonance response signal of the local structure. In addition, the dynamic characteristics of the fluid boundary layer are affected not only by the evolution characteristics of the pipeline flow field but also by the structure vibration. There is a certain difference between the pressure results measured by the pulsating pressure sensor and the state of an ideal rigid wall. From the perspective of signal processing, the error in the pulsating pressure test results caused by structural vibration is mainly manifested in the fact that the pulsating pressure spectrum shows line spectrum or bulge characteristics at the natural frequencies of the structure. Since these line spectra and bulges are caused by structural effects and are not induced by the evolution characteristics of the flow field boundary layer, they belong to a kind of sparse interference peaks. There is a large error in the results when calculating pipeline vibrations using these pulsating pressures containing pseudo-peaks. Therefore, how to eliminate the interference peaks in the pulsating pressure spectrum has become the key to accurately understanding the flow-acoustic load excitation characteristics induced by the pump source and improving the prediction of pipeline vibration response.

[0004] Currently, no one has proposed a filtering method for the measured pulsating pressure interference peaks. In engineering, smoothing and smoothing methods such as Savitzky-Golay, adjacent averaging, and percentile filtering are often used to smooth the measured pulsating pressure test data to solve the test noise problem of the measured pulsating pressure data. The above methods can greatly suppress the noise problem of the pulsating pressure test results, but they are ineffective for the interference peak problem in the pulsating pressure spectrum. Summary of the Invention

[0005] The present invention provides a method and system for eliminating pulsating pressure interference peaks in a pump source pipeline, aiming to solve the defect in the prior art that the acquisition of pulsating pressure in a ship pipeline is easily interfered by the vibration of the inner wall structure of the pipeline and the measuring equipment, resulting in a large error in the measurement result. By using the zero-phase filtering principle and the total variation filtering method, a sparse interference peak filtering model for the logarithmic spectrum of pulsating pressure amplitude is established, and the effective signal of pulsating pressure can be extracted from the amplitude spectrum containing sparse interference peaks through simple iterative calculation, improving the accuracy of the measurement result of the pulsating pressure signal.

[0006] In a first aspect, the present invention provides a method for eliminating pulsating pressure interference peaks in a pump source pipeline, including:

[0007] Collect the measured logarithmic spectrum of pulsating pressure amplitude in the pump source pipeline;

[0008] Construct a filtering vector for the measured logarithmic spectrum of pulsating pressure amplitude by using the zero-phase filtering principle;

[0009] Based on the constructed difference matrix, establish a total variation filtering model for the logarithmic spectrum of pulsating pressure amplitude for the filtering vector, and solve the total variation filtering model for the logarithmic spectrum of pulsating pressure amplitude through the principal minimization algorithm to obtain the total variation filtered logarithmic spectrum of pulsating pressure amplitude;

[0010] Obtain the smoothed logarithmic spectrum of pulsating pressure amplitude from the total variation filtered logarithmic spectrum of pulsating pressure amplitude;

[0011] Based on the measured logarithmic spectrum of pulsating pressure amplitude and the total variation filtered logarithmic spectrum of pulsating pressure amplitude, separate the test noise spectrum and the sparse interference peak spectrum, extract the effective amplitude spectrum of pulsating pressure of the total variation filtered logarithmic spectrum of pulsating pressure amplitude, and obtain the measured effective signal of pulsating pressure.

[0012] According to the method for eliminating pulsating pressure interference peaks in a pump source pipeline provided by the present invention, collecting the measured logarithmic spectrum of pulsating pressure amplitude in the pump source pipeline includes:

[0013] Install pulsating pressure sensors in the pump source pipeline to measure the measured pulsating pressure signals at each specified position in the pipeline;

[0014] Obtain the time series vector of the pulsating pressure correlation function corresponding to the measured pulsating pressure signal, and perform Fourier transform on the time series vector of the pulsating pressure correlation function to obtain the measured logarithmic spectrum of pulsating pressure amplitude.

[0015] According to the method for eliminating pulsating pressure interference peaks in a pump source pipeline provided by the present invention, constructing a filtering vector for the measured logarithmic spectrum of pulsating pressure amplitude by using the zero-phase filtering principle includes:

[0016] Obtain the vector dimension of the measured pulsating pressure amplitude logarithmic spectrum N , and the filter cut-off frequency;

[0017] Obtain the transfer coefficient according to the cosine function of the filter cut-off frequency, determine the first filter coefficient matrix from the vector dimension, and determine the second filter coefficient matrix from the vector dimension and the transfer coefficient. Both the first filter coefficient matrix and the second filter coefficient matrix are N × N dimensional banded matrices;

[0018] Obtain the filtered vector from the inverse matrix of the first filter coefficient matrix, the second filter coefficient matrix, the transpose matrix of the second filter coefficient matrix, and the measured pulsating pressure amplitude logarithmic spectrum.

[0019] A method for eliminating the pulsating pressure interference peak in the pump source pipeline provided by the present invention, based on the constructed difference matrix, establishes a total variation filtering model of the pulsating pressure amplitude logarithmic spectrum for the filtered vector, and solves the total variation filtering model of the pulsating pressure amplitude logarithmic spectrum through the principal minimization algorithm to obtain the total variation filtered pulsating pressure amplitude logarithmic spectrum, including:

[0020] Based on the vector dimension N Construct a difference matrix, and the difference matrix includes ( N -1)× N dimensional first-order difference matrix and ( N -2)× N dimensional second-order difference matrix;

[0021] Iterate the measured pulsating pressure amplitude logarithmic spectrum to obtain the iterated measured pulsating pressure amplitude logarithmic spectrum;

[0022] Initialize the iterated measured pulsating pressure amplitude logarithmic spectrum to obtain the initialized pulsating pressure amplitude logarithmic spectrum, and the initialized pulsating pressure amplitude logarithmic spectrum corresponds to the first step of iteration;

[0023] Calculate the penalty matrix of the product of the iterated measured pulsating pressure amplitude logarithmic spectrum and the difference matrix. The penalty matrix is N dimensional diagonal matrix, and each diagonal element is obtained from the iterated measured pulsating pressure amplitude logarithmic spectrum and a preset fixed value;

[0024] Determine the first penalty coefficient corresponding to the first-order difference matrix and the second penalty coefficient corresponding to the second-order difference matrix. Based on the first penalty coefficient, the second penalty coefficient, the difference matrix, the transpose matrix of the difference matrix, and the penalty matrix, obtain the iteration matrix;

[0025] Obtain the filtering vector of the current iteration step. From the filtering vector of the current iteration step, the iteration matrix, the first filtering coefficient matrix, and the second filtering coefficient matrix, obtain the logarithmic spectrum of the measured pulsating pressure amplitude of the next iteration step;

[0026] Construct a residual function from the logarithmic spectrum of the measured pulsating pressure amplitude of the iteration, the logarithmic spectrum of the measured pulsating pressure amplitude, the first filtering coefficient matrix, the second filtering coefficient matrix, the first penalty coefficient, the first-order difference matrix, the second penalty coefficient, and the second-order difference matrix. Perform iterative calculations on the residual function until the preset convergence condition is met to obtain the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering.

[0027] According to a method for eliminating pulsating pressure interference peaks in a pump source pipeline provided by the present invention, from the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering, obtain a smoothed logarithmic spectrum of the pulsating pressure amplitude, including:

[0028] Use the moving least squares method to perform fitting and smoothing on the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering to obtain the smoothed logarithmic spectrum of the pulsating pressure amplitude.

[0029] According to a method for eliminating pulsating pressure interference peaks in a pump source pipeline provided by the present invention, based on the logarithmic spectrum of the measured pulsating pressure amplitude and the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering, separate to obtain a test noise spectrum and a sparse interference peak spectrum, and extract the effective amplitude spectrum of the pulsating pressure of the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering to obtain a measured effective pulsating pressure signal, including:

[0030] Separate the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering in the logarithmic spectrum of the measured pulsating pressure amplitude, and combine the first filtering coefficient matrix and the second filtering coefficient matrix to obtain the test noise spectrum;

[0031] Separate the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering and the test noise spectrum in the logarithmic spectrum of the measured pulsating pressure amplitude in sequence to obtain the sparse interference peak spectrum;

[0032] Determine the pulsating pressure reference value, and based on the pulsating pressure reference value and the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering, obtain the effective amplitude spectrum of the pulsating pressure;

[0033] Use the smoothed logarithmic spectrum of the pulsating pressure amplitude as the measured effective pulsating pressure signal.

[0034] In a second aspect, the present invention also provides a system for eliminating pulsating pressure interference peaks in a pump source pipeline, including:

[0035] An acquisition module for acquiring the logarithmic spectrum of the measured pulsating pressure amplitude in the pump source pipeline;

[0036] A construction module for constructing a filtering vector of the logarithmic spectrum of the measured pulsating pressure amplitude by using the zero-phase filtering principle;

[0037] A calculation module for establishing a total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude for the filtering vector based on the constructed difference matrix, and solving the total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude through the principal minimization algorithm to obtain the total variation filtered logarithmic spectrum of the pulsating pressure amplitude;

[0038] A smoothing module for obtaining a smoothed logarithmic spectrum of the pulsating pressure amplitude from the total variation filtered logarithmic spectrum of the pulsating pressure amplitude;

[0039] An elimination module for separating and obtaining a test noise spectrum and a sparse interference peak spectrum based on the logarithmic spectrum of the measured pulsating pressure amplitude and the total variation filtered logarithmic spectrum of the pulsating pressure amplitude, extracting the effective amplitude spectrum of the pulsating pressure of the total variation filtered logarithmic spectrum of the pulsating pressure amplitude, and obtaining the measured effective pulsating pressure signal.

[0040] In a third aspect, the present invention further provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor. When the processor executes the program, the method for eliminating the pulsating pressure interference peak in the pump source pipeline as described in any one of the above is implemented.

[0041] In a fourth aspect, the present invention further provides a non-transitory computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the method for eliminating the pulsating pressure interference peak in the pump source pipeline as described in any one of the above is implemented.

[0042] A method and system for eliminating the pulsating pressure interference peak in a pump source pipeline provided by the present invention utilize the zero-phase filtering principle and the total variation filtering method to establish a sparse interference peak filtering model for the logarithmic spectrum of the pulsating pressure amplitude, and then form an iterative calculation method for the effective signal of the pulsating pressure amplitude spectrum through the principal minimization algorithm. Through iterative calculation, the effective signal of the pulsating pressure can be extracted from the amplitude spectrum containing sparse interference peaks, improving the accuracy of the measurement result of the pulsating pressure signal. Aiming at the problem of signal distortion in the result obtained by the total variation filtering model, the present invention proposes to perform secondary processing on the pulsating pressure obtained by the filtering model by using the moving least squares method, solves the divergence problem of the extraction result of the filtering model, and enhances the reliability of the pulsating pressure filtering algorithm. Description of the Drawings

[0043] To more clearly illustrate the technical solutions in the present invention or the prior art, the following will briefly introduce the drawings required for the description of the embodiments or the prior art. Obviously, the drawings in the following description are some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0044] Figure 1 is a schematic flow chart of a method for eliminating pulsating pressure interference peaks in a pump source pipeline provided by the present invention;

[0045] Figure 2 is a diagram of the straight pipe structure and measuring point positions of a pump source pipeline provided by the present invention;

[0046] Figure 3 is a diagram of the filtering decomposition result provided by the present invention;

[0047] Figure 4 is a comparison diagram of the filtering result and the theoretical value of pulsating pressure provided by the present invention;

[0048] Figure 5 is a schematic structural diagram of a system for eliminating pulsating pressure interference peaks in a pump source pipeline provided by the present invention;

[0049] Figure 6 is a schematic structural diagram of an electronic device provided by the present invention. Detailed Embodiments

[0050] To make the objectives, technical solutions, and advantages of the present invention clearer, the following will clearly and completely describe the technical solutions in the present invention in conjunction with the drawings in the present invention. Obviously, the described embodiments are some, but not all, of the embodiments of the present invention. Based on the embodiments in the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts fall within the scope of protection of the present invention.

[0051] Figure 1 is a schematic flow chart of eliminating pulsating pressure interference peaks in a pump source pipeline provided by an embodiment of the present invention. As Figure 1 shown, it includes:

[0052] Step 100: Collect the measured pulsating pressure amplitude logarithmic spectrum in the pump source pipeline;

[0053] Step 200: Construct a filtering vector for the measured pulsating pressure amplitude logarithmic spectrum using the zero-phase filtering principle;

[0054] Step 300: Based on the constructed difference matrix, establish a total variation filtering model for the logarithmic spectrum of the pulsating pressure amplitude of the filtering vector, and solve the total variation filtering model for the logarithmic spectrum of the pulsating pressure amplitude through the principal minimization algorithm to obtain the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering;

[0055] Step 400: Obtain the smoothed logarithmic spectrum of the pulsating pressure amplitude from the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering;

[0056] Step 500: Based on the measured logarithmic spectrum of the pulsating pressure amplitude and the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering, separate the test noise spectrum and the sparse interference peak spectrum, extract the effective amplitude spectrum of the pulsating pressure of the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering, and obtain the measured effective pulsating pressure signal.

[0057] Aiming at the problems that the pulsating pressure test results in current engineering are affected by structural vibration, there are interference peaks and the results are inaccurate, an embodiment of the present invention proposes a method for eliminating wave of sparse interference peaks in the frequency domain of measured pulsating pressure for the pump source pipeline in a ship. This method uses the zero-phase filtering principle and the total variation filtering method to establish a sparse interference peak filtering model for the logarithmic spectrum of the pulsating pressure amplitude, and then through the principal minimization algorithm, forms an iterative calculation method for the effective signal of the pulsating pressure amplitude spectrum. Through iterative calculation, the effective signal of the pulsating pressure can be extracted from the amplitude spectrum containing sparse interference peaks, improving the accuracy of the measurement result of the pulsating pressure signal.

[0058] Based on the above embodiment, step 100 includes:

[0059] Arrange pulsating pressure sensors in the pump source pipeline to test the measured pulsating pressure signals at each specified position in the pipeline;

[0060] Obtain the time series vector of the pulsating pressure correlation function corresponding to the measured pulsating pressure signal, and perform Fourier transform on the time series vector of the pulsating pressure correlation function to obtain the measured logarithmic spectrum of the pulsating pressure amplitude.

[0061] Specifically, for the pump source pipeline in actual engineering, by arranging pulsating pressure sensors to test the pulsating pressure signals at typical positions on the pipeline, solve the time series vector of the pulsating pressure correlation function at each position (denoted as p ), and perform Fourier transform on it to obtain the measured logarithmic spectrum of the pulsating pressure amplitude .

[0062] Taking Figure 2Taking the straight pipe as an example, the specific implementation process of the present invention is described. The outer diameter of the straight pipe is 52 mm, the inner diameter is 49 mm, and the length is 1580 mm. The two ends of the bent pipe are fixed, filled with water, and the flow velocity of the water medium is about 10 m / s. An intrusive pulsating pressure sensor is installed at point P1 of the pipeline, and the logarithmic spectrum data of the pulsating pressure amplitude at point P1 is measured as shown in Figure 3 as shown in (a) of

[0063] Based on the above embodiments, step 200 includes:

[0064] Obtain the vector dimension of the measured logarithmic spectrum of the pulsating pressure amplitude N , and the filter cut-off frequency;

[0065] Obtain the transfer coefficient according to the cosine function of the filter cut-off frequency, determine the first filter coefficient matrix from the vector dimension, and determine the second filter coefficient matrix from the vector dimension and the transfer coefficient. Both the first filter coefficient matrix and the second filter coefficient matrix are N × N dimensional banded matrices;

[0066] Obtain the filtered vector from the inverse matrix of the first filter coefficient matrix, the second filter coefficient matrix, the transpose matrix of the second filter coefficient matrix, and the measured logarithmic spectrum of the pulsating pressure amplitude.

[0067] Specifically, in the embodiments of the present invention, the filtered vector of the logarithmic spectrum of the pulsating pressure amplitude is constructed using the zero-phase filtering principle :

[0068]

[0069] In the formula, is the measured logarithmic spectrum of the pulsating pressure amplitude, N is the vector dimension (sequence length), is the angular frequency; A and B are both N × N dimensional banded matrices, representing the first filter coefficient matrix and the second filter coefficient matrix respectively.

[0070]

[0071]

[0072] In the formula, α is the transfer coefficient, which can be expressed as

[0073]

[0074] In the formula, represents the filter cut-off frequency, and it is recommended to take 0.01π in the embodiments of the present invention.

[0075] Based on the above embodiments, step 300 includes:

[0076] Based on the vector dimension N Construct a difference matrix, the difference matrix includes( N -1)× N dimensional first-order difference matrix and( N -2)× N dimensional second-order difference matrix;

[0077] Iterate the measured pulsating pressure amplitude logarithmic spectrum to obtain the iterated measured pulsating pressure amplitude logarithmic spectrum;

[0078] Initialize the iterated measured pulsating pressure amplitude logarithmic spectrum to obtain the initialized pulsating pressure amplitude logarithmic spectrum, and the initialized pulsating pressure amplitude logarithmic spectrum corresponds to the first step of iteration;

[0079] Calculate the penalty matrix of the product of the iterated measured pulsating pressure amplitude logarithmic spectrum and the difference matrix, and the penalty matrix is N dimensional diagonal matrix, and each diagonal element is obtained from the iterated measured pulsating pressure amplitude logarithmic spectrum and a preset fixed value;

[0080] Determine the first penalty coefficient corresponding to the first-order difference matrix and the second penalty coefficient corresponding to the second-order difference matrix. Based on the first penalty coefficient, the second penalty coefficient, the difference matrix, the transpose matrix of the difference matrix, and the penalty matrix, obtain the iteration matrix;

[0081] Obtain the filtering vector of the current iteration step. From the filtering vector of the current iteration step, the iteration matrix, the first filtering coefficient matrix, and the second filtering coefficient matrix, obtain the measured pulsating pressure amplitude logarithmic spectrum of the next iteration step;

[0082] Construct a residual function from the iterated measured pulsating pressure amplitude logarithmic spectrum, the measured pulsating pressure amplitude logarithmic spectrum, the first filtering coefficient matrix, the second filtering coefficient matrix, the first penalty coefficient, the first-order difference matrix, the second penalty coefficient, and the second-order difference matrix, and perform iterative calculation on the residual function until the preset convergence condition is satisfied to obtain the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering.

[0083] Specifically, first construct the first-order difference matrix D1 and the second-order difference matrix D2:

[0084]

[0085]

[0086] Wherein, D1 represents( N -1)×N The first-order difference matrix of dimension; D2 represents ( N -2)× N The second-order difference matrix of dimension.

[0087] Then, establish the total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude, and solve it through the following principal minimization algorithm:

[0088] (1) Initialization, input , when , , represents the number of iteration steps, is a positive integer greater than or equal to 1. When , represents the measured logarithmic spectrum of the pulsating pressure amplitude in the iteration.

[0089] (2) Calculate the penalty matrix ( i = 1, 2). Here, is the simplified

[0090]

[0091] In the formula, is

[0092]

[0093] In the embodiments of the present invention It is recommended to take .

[0094] (3) Calculate the iteration matrix M:

[0095]

[0096] In the formula, is the penalty coefficient. In the embodiments of the present invention, it is recommended , .

[0097] (4) Calculate :

[0098]

[0099] (5) Calculate the residual function :

[0100]

[0101] Here, is the simplified ;

[0102] Continue the above iteration until the convergence condition is met, i.e.,

[0103]

[0104] In the formula, is 0, ε is the convergence determination criterion, and in the embodiments of the present invention, it is recommended to take .

[0105] Through the calculation of the main minimization algorithm, the measured logarithmic spectrum of the pulsating pressure amplitude after total variation filtering can be obtained .

[0106] Based on the above embodiments, step 400 includes:

[0107] Use the moving least squares method to perform fitting and smoothing processing on the filtered logarithmic spectrum of the pulsating pressure amplitude to obtain the smoothed logarithmic spectrum of the pulsating pressure amplitude .

[0108] It should be noted that the fitting and smoothing processing adopted in the embodiments of the present invention is an existing algorithm, which is a geometric processing method combining curve / surface fitting and smoothing optimization, mainly used to improve the smoothness of the model and maintain geometric features. Its core is to eliminate noise or local irregularities while approximating the original data through mathematical optimization means, and will not be elaborated here.

[0109] Based on the above embodiments, step 500 includes:

[0110] Separate the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering in the measured logarithmic spectrum of the pulsating pressure amplitude, and combine the first filter coefficient matrix and the second filter coefficient matrix to obtain the test noise spectrum;

[0111] Separate the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering and the test noise spectrum in the measured logarithmic spectrum of the pulsating pressure amplitude in sequence to obtain the sparse interference peak spectrum;

[0112] Determine the pulsating pressure reference value, and based on the pulsating pressure reference value and the logarithmic spectrum of the pulsating pressure amplitude after total variation filtering, obtain the effective amplitude spectrum of the pulsating pressure;

[0113] Use the smoothed logarithmic spectrum of the pulsating pressure amplitude as the measured effective signal of the pulsating pressure.

[0114] Specifically, the test noise spectrum separated in the embodiments of the present invention is , the obtained frequency-domain sparse interference peak spectrum is , and the logarithmic spectrum of the effective amplitude of the pulsating pressure is .

[0115] By the above process, the influence of sparse interference peaks and test noise can be completed, and the effective amplitude spectrum y( ω ) of the pulsating pressure can be obtained as

[0116]

[0117] In the formula, is the reference value of the pulsating pressure.

[0118] Using the method proposed by the present invention to separate the logarithm of the measured pulsating pressure, the effective amplitude logarithm spectrum, the sparse interference peak spectrum, and the test noise spectrum of the pulsating pressure are respectively as Figure 3 shown in (b), (c), and (d) in. To prove the effectiveness of the pulsating pressure spectrum obtained by the proposed method, the effective amplitude logarithm spectrum of the pulsating pressure obtained by filtering is compared with the theoretical value of the logarithm spectrum of the pulsating pressure amplitude in a straight pipe under the same flow rate. The results are as Figure 4 shown. It can be seen from Figure 4 that the error between the filtering result of the proposed method and the theoretical value of the pulsating pressure is very small, less than 1 dB (10 - 1000 Hz), which proves the effectiveness of the method proposed by the present invention.

[0119] Next, the system for eliminating the pulsating pressure interference peak in the pump source pipeline provided by the present invention will be described. The system for eliminating the pulsating pressure interference peak in the pump source pipeline described below can be mutually referred to the method for eliminating the pulsating pressure interference peak in the pump source pipeline described above.

[0120] Figure 5 is a schematic structural diagram of the system for eliminating the pulsating pressure interference peak in the pump source pipeline provided by an embodiment of the present invention. As Figure 5 shown, it includes: an acquisition module 51, a construction module 52, a calculation module 53, a smoothing module 54, and an elimination module 55, where:

[0121] The acquisition module 51 is used to acquire the logarithmic spectrum of the measured pulsating pressure amplitude in the pump source pipeline; the construction module 52 is used to construct the filtering vector of the logarithmic spectrum of the measured pulsating pressure amplitude by using the zero-phase filtering principle; the calculation module 53 is used to establish a total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude for the filtering vector based on the constructed difference matrix, solve the total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude through the principal minimization algorithm, and obtain the total variation filtered logarithmic spectrum of the pulsating pressure amplitude; the smoothing module 54 is used to obtain the smoothed logarithmic spectrum of the pulsating pressure amplitude from the total variation filtered logarithmic spectrum of the pulsating pressure amplitude; the elimination module 55 is used to separate and obtain the test noise spectrum and the sparse interference peak spectrum based on the measured logarithmic spectrum of the pulsating pressure amplitude and the total variation filtered logarithmic spectrum of the pulsating pressure amplitude, extract the effective amplitude spectrum of the pulsating pressure of the total variation filtered logarithmic spectrum of the pulsating pressure amplitude, and obtain the measured effective pulsating pressure signal.

[0122] Figure 6 An example of the physical structure diagram of an electronic device is shown in Figure 6 As shown, the electronic device may include: a processor 610, a communication interface 620, a memory 630, and a communication bus 640. Among them, the processor 610, the communication interface 620, and the memory 630 complete communication with each other through the communication bus 640. The processor 610 can call the logical instructions in the memory 630 to execute the method for eliminating the pulsating pressure interference peak in the pump source pipeline. The method includes: acquiring the logarithmic spectrum of the measured pulsating pressure amplitude in the pump source pipeline; constructing the filtering vector of the logarithmic spectrum of the measured pulsating pressure amplitude by using the zero-phase filtering principle; establishing a total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude for the filtering vector based on the constructed difference matrix, solving the total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude through the principal minimization algorithm, and obtaining the total variation filtered logarithmic spectrum of the pulsating pressure amplitude; obtaining the smoothed logarithmic spectrum of the pulsating pressure amplitude from the total variation filtered logarithmic spectrum of the pulsating pressure amplitude; separating and obtaining the test noise spectrum and the sparse interference peak spectrum based on the measured logarithmic spectrum of the pulsating pressure amplitude and the total variation filtered logarithmic spectrum of the pulsating pressure amplitude, extracting the effective amplitude spectrum of the pulsating pressure of the total variation filtered logarithmic spectrum of the pulsating pressure amplitude, and obtaining the measured effective pulsating pressure signal.

[0123] In addition, when the logical instructions in the above-mentioned memory 630 are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present invention. The aforementioned storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), magnetic disks, or optical discs that can store program codes.

[0124] On the other hand, the present invention also provides a computer program product. The computer program product includes a computer program that can be stored on a non-transitory computer-readable storage medium. When the computer program is executed by a processor, the computer can execute the method for eliminating the pulsating pressure interference peak in the pump source pipeline provided by the above-mentioned various methods. The method includes: collecting the measured pulsating pressure amplitude logarithmic spectrum in the pump source pipeline; constructing a filtering vector of the measured pulsating pressure amplitude logarithmic spectrum using the zero-phase filtering principle; based on the constructed difference matrix, establishing a total variation filtering model of the pulsating pressure amplitude logarithmic spectrum for the filtering vector, and solving the total variation filtering model of the pulsating pressure amplitude logarithmic spectrum through the principal minimization algorithm to obtain the total variation filtered pulsating pressure amplitude logarithmic spectrum; obtaining the smoothed pulsating pressure amplitude logarithmic spectrum from the total variation filtered pulsating pressure amplitude logarithmic spectrum; based on the measured pulsating pressure amplitude logarithmic spectrum and the total variation filtered pulsating pressure amplitude logarithmic spectrum, separating to obtain the test noise spectrum and the sparse interference peak spectrum, and extracting the effective pulsating pressure amplitude spectrum of the total variation filtered pulsating pressure amplitude logarithmic spectrum to obtain the measured effective pulsating pressure signal.

[0125] In another aspect, the present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, it realizes the method for eliminating the pulsating pressure interference peaks in the pump source pipeline provided by the above-mentioned various methods. The method includes: collecting the measured pulsating pressure amplitude logarithmic spectrum in the pump source pipeline; constructing a filtering vector of the measured pulsating pressure amplitude logarithmic spectrum by using the zero-phase filtering principle; based on the constructed difference matrix, establishing a total variation filtering model of the pulsating pressure amplitude logarithmic spectrum for the filtering vector, and solving the total variation filtering model of the pulsating pressure amplitude logarithmic spectrum through the principal minimization algorithm to obtain the total variation filtered pulsating pressure amplitude logarithmic spectrum; obtaining the smoothed pulsating pressure amplitude logarithmic spectrum from the total variation filtered pulsating pressure amplitude logarithmic spectrum; based on the measured pulsating pressure amplitude logarithmic spectrum and the total variation filtered pulsating pressure amplitude logarithmic spectrum, separating to obtain the test noise spectrum and the sparse interference peak spectrum, extracting the pulsating pressure effective amplitude spectrum of the total variation filtered pulsating pressure amplitude logarithmic spectrum, and obtaining the measured pulsating pressure effective signal. The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. Those of ordinary skill in the art can understand and implement it without creative labor.

[0126] Through the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on such an understanding, the essence of the above technical solution, or the part that contributes to the prior art, can be embodied in the form of a software product. The computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0127] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit them. Although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments or perform equivalent replacements for some of the technical features. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for eliminating the pulsating pressure interference peak in a pump source pipeline, characterized in that, Including: Collecting the logarithmic spectrum of the measured pulsating pressure amplitude in the pump source pipeline; Constructing a filtering vector for the logarithmic spectrum of the measured pulsating pressure amplitude by using the zero-phase filtering principle; Based on the constructed difference matrix, establishing a total variation filtering model for the logarithmic spectrum of the pulsating pressure amplitude with respect to the filtering vector, and solving the total variation filtering model for the logarithmic spectrum of the pulsating pressure amplitude through the principal minimization algorithm to obtain the total variation filtered logarithmic spectrum of the pulsating pressure amplitude; Obtaining the smoothed logarithmic spectrum of the pulsating pressure amplitude from the total variation filtered logarithmic spectrum of the pulsating pressure amplitude; Based on the logarithmic spectrum of the measured pulsating pressure amplitude and the total variation filtered logarithmic spectrum of the pulsating pressure amplitude, separating to obtain the test noise spectrum and the sparse interference peak spectrum, extracting the effective pulsating pressure amplitude spectrum of the total variation filtered logarithmic spectrum of the pulsating pressure amplitude to obtain the measured effective pulsating pressure signal; Constructing a filtering vector for the logarithmic spectrum of the measured pulsating pressure amplitude by using the zero-phase filtering principle, including: Obtain the vector dimension of the measured pulsating pressure amplitude logarithmic spectrum N , and the filter cut-off frequency; The transfer coefficient is obtained according to the cosine function of the filtering cut-off frequency, the first filtering coefficient matrix is determined by the vector dimension, and the second filtering coefficient matrix is determined by the vector dimension and the transfer coefficient. Both the first filtering coefficient matrix and the second filtering coefficient matrix are N × N dimensional banded matrices; Obtaining the filtering vector from the inverse matrix of the first filtering coefficient matrix, the second filtering coefficient matrix, the transpose matrix of the second filtering coefficient matrix, and the logarithmic spectrum of the measured pulsating pressure amplitude.

2. The method for eliminating the pulsating pressure interference peak in the pump source pipeline according to claim 1, characterized in that Collecting the logarithmic spectrum of the measured pulsating pressure amplitude in the pump source pipeline, including: Deploying pulsating pressure sensors in the pump source pipeline to measure the measured pulsating pressure signals at each specified position in the pipeline; Obtaining the time series vector of the pulsating pressure correlation function corresponding to the measured pulsating pressure signal, and performing Fourier transform on the time series vector of the pulsating pressure correlation function to obtain the logarithmic spectrum of the measured pulsating pressure amplitude.

3. The method for eliminating the pulsating pressure interference peak in the pump source pipeline according to claim 1, wherein, Based on the constructed difference matrix, establishing a total variation filtering model for the logarithmic spectrum of the pulsating pressure amplitude with respect to the filtering vector, and solving the total variation filtering model for the logarithmic spectrum of the pulsating pressure amplitude through the principal minimization algorithm to obtain the total variation filtered logarithmic spectrum of the pulsating pressure amplitude, including: Based on the vector dimension N Construct a difference matrix, the difference matrix including( N -1)× N a first-order difference matrix of dimension and( N -2)× N a second-order difference matrix of dimension; Iterating the logarithmic spectrum of the measured pulsating pressure amplitude to obtain the iterated logarithmic spectrum of the measured pulsating pressure amplitude; Initializing the iterated logarithmic spectrum of the measured pulsating pressure amplitude to obtain the initialized logarithmic spectrum of the pulsating pressure amplitude, and the initialized logarithmic spectrum of the pulsating pressure amplitude corresponds to the first step of the iteration; Calculate the penalty matrix of the product of the measured pulsating pressure amplitude logarithmic spectrum of the iteration and the difference matrix, where the penalty matrix is N a diagonal matrix of dimension, and each diagonal element is obtained from the measured pulsating pressure amplitude logarithmic spectrum of the iteration and a preset fixed value; Determining the first penalty coefficient corresponding to the first-order difference matrix and the second penalty coefficient corresponding to the second-order difference matrix, and obtaining the iteration matrix based on the first penalty coefficient, the second penalty coefficient, the difference matrix, the transpose matrix of the difference matrix, and the penalty matrix; Obtaining the filtering vector of the current iteration step, and obtaining the logarithmic spectrum of the measured pulsating pressure amplitude of the next iteration step from the filtering vector of the current iteration step, the iteration matrix, the first filtering coefficient matrix, and the second filtering coefficient matrix; Constructing a residual function from the iterated logarithmic spectrum of the measured pulsating pressure amplitude, the logarithmic spectrum of the measured pulsating pressure amplitude, the first filtering coefficient matrix, the second filtering coefficient matrix, the first penalty coefficient, the first-order difference matrix, the second penalty coefficient, and the second-order difference matrix, and performing iterative calculation on the residual function until the preset convergence condition is satisfied to obtain the total variation filtered logarithmic spectrum of the pulsating pressure amplitude.

4. The method for eliminating the pulsating pressure interference peak in the pump source pipeline according to claim 1, wherein From the logarithm spectrum of the pulsating pressure amplitude after the total variation filtering, obtain the smoothed logarithm spectrum of the pulsating pressure amplitude, including: Perform fitting and smoothing processing on the logarithm spectrum of the pulsating pressure amplitude after the total variation filtering by using the moving least squares method to obtain the smoothed logarithm spectrum of the pulsating pressure amplitude.

5. The method for eliminating the pulsating pressure interference peak in the pump source pipeline according to claim 3, characterized in that Based on the measured logarithm spectrum of the pulsating pressure amplitude and the logarithm spectrum of the pulsating pressure amplitude after the total variation filtering, separately obtain the test noise spectrum and the sparse interference peak spectrum, and extract the effective amplitude spectrum of the pulsating pressure of the logarithm spectrum of the pulsating pressure amplitude after the total variation filtering to obtain the measured effective pulsating pressure signal, including: Separate the logarithm spectrum of the pulsating pressure amplitude after the total variation filtering from the measured logarithm spectrum of the pulsating pressure amplitude, and combine the first filter coefficient matrix and the second filter coefficient matrix to obtain the test noise spectrum; Successively separate the logarithm spectrum of the pulsating pressure amplitude after the total variation filtering and the test noise spectrum from the measured logarithm spectrum of the pulsating pressure amplitude to obtain the sparse interference peak spectrum; Determine the pulsating pressure reference value, and based on the pulsating pressure reference value and the logarithm spectrum of the pulsating pressure amplitude after the total variation filtering, obtain the effective amplitude spectrum of the pulsating pressure; Use the smoothed logarithm spectrum of the pulsating pressure amplitude as the measured effective pulsating pressure signal.

6. A system for eliminating pulsating pressure interference peaks in a pump source pipeline, based on the method for eliminating pulsating pressure interference peaks in a pump source pipeline according to any one of claims 1 to 5, characterized in that, Include: An acquisition module for acquiring the measured logarithm spectrum of the pulsating pressure amplitude in the pump source pipeline; A construction module for constructing a filter vector of the measured logarithm spectrum of the pulsating pressure amplitude by using the zero-phase filtering principle; A calculation module for establishing a total variation filtering model of the logarithm spectrum of the pulsating pressure amplitude for the filter vector based on the constructed difference matrix, and solving the total variation filtering model of the logarithm spectrum of the pulsating pressure amplitude by using the principal minimization algorithm to obtain the logarithm spectrum of the pulsating pressure amplitude after the total variation filtering; A smoothing module for obtaining the smoothed logarithm spectrum of the pulsating pressure amplitude from the logarithm spectrum of the pulsating pressure amplitude after the total variation filtering; An elimination module for separately obtaining the test noise spectrum and the sparse interference peak spectrum based on the measured logarithm spectrum of the pulsating pressure amplitude and the logarithm spectrum of the pulsating pressure amplitude after the total variation filtering, and extracting the effective amplitude spectrum of the pulsating pressure of the logarithm spectrum of the pulsating pressure amplitude after the total variation filtering to obtain the measured effective pulsating pressure signal.

7. An electronic device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein When the processor executes the program, it implements the method for eliminating the pulsating pressure interference peak in the pump source pipeline according to any one of claims 1 to 5.

8. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the method for eliminating the pulsating pressure interference peak in the pump source pipeline according to any one of claims 1 to 5.

9. A computer program product, comprising a computer program, characterized in that, When the computer program is executed by the processor, it implements the method for eliminating the pulsating pressure interference peak in the pump source pipeline according to any one of claims 1 to 5.

Images