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

Abstract

The invention provides a method and a system for eliminating a pulsating pressure interference peak in a pump source pipeline. A sparse interference peak filtering model for a pulsating pressure amplitude logarithmic spectrum is established by using a zero-phase filtering principle and a total variation filtering method, and then a pulsating pressure amplitude spectrum effective signal iterative calculation method is formed through a main minimization algorithm. The effective signal of the pulsating pressure can be extracted from the amplitude spectrum containing the sparse interference peak through iterative calculation, and the accuracy of a pulsating pressure signal measurement result is improved. In order to solve the problem of signal distortion of a result obtained by a total variation filtering model, the invention proposes to carry out secondary processing on the pulsating pressure obtained by the filtering model by using a moving least square method, so that the divergence problem of the extraction result of the filtering model is solved, and the reliability of a pulsating pressure filtering algorithm is enhanced.

Description

Technical Field

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

[0002] In daily life, the acoustic pressure waves induced by the pump source and the turbulent pulsating pressure in the fluid boundary layer are important excitation sources for ship pipelines, which often induce vibrations of the pipelines and reduce their acoustic performance. Therefore, how to accurately obtain the pulsating pressure characteristics of the inner wall of the pipe driven by the pump source is crucial to mastering the flow-acoustic load excitation characteristics of the pump source and solving the pipeline vibration problem.

[0003] At present, hydrophones or pulsating pressure sensors are often used in engineering to test 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 correlation spectrum analysis and coherence analysis methods. However, since the pulsating pressure sensor or hydrophone often needs to be installed on the pipeline in an invasive manner, the pulsating pressure sensor itself will be affected by the vibration of the pipeline structure during the measurement process, resulting in the signal received by the sensor not only containing the pulsating pressure signal of the pipe wall, but also receiving the resonance response signal of the local structure. In addition, the dynamic characteristics of the fluid boundary layer are not only affected by the evolution characteristics of the pipeline flow field but also by the structural vibration. There is a certain difference between the pressure results tested by the pulsating pressure sensor and the ideal rigid wall state. From the perspective of signal processing, the error of the pulsating pressure test results caused by structural vibration is mainly manifested in the pulsating pressure spectrum showing line spectrum or bulge characteristics at the natural frequency of the structure. Since these line spectra and bulges are caused by structural influences, rather than induced by the evolution characteristics of the flow field boundary layer, they are a kind of sparse interference peaks. When using these pulsating pressures containing pseudo-peaks for pipeline vibration calculation, there is a large error in the results. Therefore, how to eliminate the interference peaks in the pulsating pressure spectrum becomes the key to accurately grasp the flow-acoustic load excitation characteristics induced by the pump source and improve the prediction of pipeline vibration response.

[0004] At present, no one has proposed a filtering method for the interference peak of the measured pulsating pressure. In engineering, smoothing methods such as Savitzky-Golay, adjacent average 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 it is difficult to work 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 the interference peak of pulsating pressure in a pump source pipeline, so as to solve the defect in the prior art that the acquisition of pulsating pressure in a ship pipeline is easily affected by the vibration of the inner wall structure of the pipeline and the interference of the measuring equipment, thereby causing 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 the pulsating pressure amplitude is established, so that the effective signal of the pulsating pressure can be extracted from the amplitude spectrum containing the sparse interference peak through simple iterative calculation, thereby improving the accuracy of the pulsating pressure signal measurement result.

[0006] In a first aspect, the present invention provides a method for eliminating a pulsating pressure interference peak in a pump source pipeline, comprising: Collect the logarithmic spectrum of the measured pulsating pressure amplitude in the pump source pipeline; Constructing a filter vector of the logarithmic spectrum of the measured pulsating pressure amplitude by using the zero-phase filtering principle; Based on the constructed difference matrix, a total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude is established for the filtering vector, and the total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude is solved by the main minimization algorithm to obtain the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering; Obtaining a smoothed logarithmic spectrum of the pulsating pressure amplitude from the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering; Based on the measured pulsating pressure amplitude logarithmic spectrum and the pulsating pressure amplitude logarithmic spectrum after the total variation filtering, the test noise spectrum and the sparse interference peak spectrum are separated, and the pulsating pressure effective amplitude spectrum of the pulsating pressure amplitude logarithmic spectrum after the total variation filtering is extracted to obtain the measured pulsating pressure effective signal.

[0007] According to a method for eliminating a pulsating pressure interference peak in a pump source pipeline provided by the present invention, a logarithmic spectrum of the measured pulsating pressure amplitude in the pump source pipeline is collected, comprising: A pulsating pressure sensor is arranged in the pump source pipeline to test the measured pulsating pressure signal at each designated position in the pipeline; A pulsating pressure correlation function time series vector corresponding to the measured pulsating pressure signal is obtained, and the pulsating pressure correlation function time series vector is subjected to Fourier transformation to obtain the measured pulsating pressure amplitude logarithmic spectrum.

[0008] According to a method for eliminating a pulsating pressure interference peak in a pump source pipeline provided by the present invention, a filtering vector of the logarithmic spectrum of the measured pulsating pressure amplitude is constructed by using the zero-phase filtering principle, comprising: Obtain the vector dimension of the logarithmic spectrum of the measured pulsating pressure amplitude N , and the filter cutoff frequency; The transfer coefficient is obtained according to the cosine function of the filter cutoff frequency, the first filter coefficient matrix is ​​determined by the vector dimension, and the second filter coefficient matrix is ​​determined by the vector dimension and the transfer coefficient, and the first filter coefficient matrix and the second filter coefficient matrix are both N × N dimensional band matrix; A filter vector is obtained from the inverse matrix of the first filter coefficient matrix, the second filter coefficient matrix, the transposed matrix of the second filter coefficient matrix, and the logarithmic spectrum of the measured pulsating pressure amplitude.

[0009] The present invention provides a method for eliminating a pulsating pressure interference peak in a pump source pipeline. Based on a constructed differential matrix, a total variation filtering model of a pulsating pressure amplitude logarithmic spectrum for the filtering vector is established. The total variation filtering model of the pulsating pressure amplitude logarithmic spectrum is solved by a principal minimization algorithm to obtain a pulsating pressure amplitude logarithmic spectrum after total variation filtering, including: Based on vector dimension N Construct a difference matrix, the difference matrix includes ( N -1)× N The first-order difference matrix of dimension and ( N -2)× N dimensional second-order difference matrix; Iterating the measured logarithmic spectrum of the pulsating pressure amplitude to obtain an iterative logarithmic spectrum of the measured pulsating pressure amplitude; Initializing the iterative measured pulsating pressure amplitude logarithmic spectrum to obtain an initialized pulsating pressure amplitude logarithmic spectrum, wherein the initialized pulsating pressure amplitude logarithmic spectrum corresponds to the first step of the iteration; Calculate the penalty matrix of the product of the iterative measured pulsating pressure amplitude log spectrum and the difference matrix, the penalty matrix is N dimensional diagonal matrix, each diagonal element is obtained by the iterative measured pulsating pressure amplitude logarithmic spectrum and a preset fixed value; Determine a first penalty coefficient corresponding to the first-order difference matrix and a second penalty coefficient corresponding to the second-order difference matrix, and obtain an iterative matrix based on the first penalty coefficient, the second penalty coefficient, the difference matrix, a transposed matrix of the difference matrix, and the penalty matrix; Obtaining a filter vector of a current iteration step, and obtaining a logarithmic spectrum of a measured pulsating pressure amplitude of a next iteration step from the filter vector of the current iteration step, the iteration matrix, the first filter coefficient matrix, and the second filter coefficient matrix; A residual function is constructed by the iterative measured pulsating pressure amplitude logarithmic spectrum, the measured pulsating pressure amplitude logarithmic spectrum, the first filter coefficient matrix, the second filter coefficient matrix, the first penalty coefficient, the first-order difference matrix, the second penalty coefficient and the second-order difference matrix, and the residual function is calculated and iterated until a preset convergence condition is met to obtain the pulsating pressure amplitude logarithmic spectrum after total variation filtering.

[0010] According to a method for eliminating a pulsating pressure interference peak in a pump source pipeline provided by the present invention, a smoothed pulsating pressure amplitude logarithmic spectrum is obtained from the pulsating pressure amplitude logarithmic spectrum after total variation filtering, comprising: The moving least square method is used to fit and smooth the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering to obtain the smoothed logarithmic spectrum of the pulsating pressure amplitude.

[0011] According to a method for eliminating a pulsating pressure interference peak in a pump source pipeline provided by the present invention, based on the measured pulsating pressure amplitude logarithmic spectrum and the pulsating pressure amplitude logarithmic spectrum after total variation filtering, a test noise spectrum and a sparse interference peak spectrum are separated, and a pulsating pressure effective amplitude spectrum of the pulsating pressure amplitude logarithmic spectrum after total variation filtering is extracted to obtain a measured pulsating pressure effective signal, including: Separating the pulsating pressure amplitude logarithmic spectrum after the total variation filtering from the measured pulsating pressure amplitude logarithmic spectrum, and combining the first filter coefficient matrix and the second filter coefficient matrix to obtain the test noise spectrum; Separating the pulsating pressure amplitude logarithmic spectrum after the total variation filtering and the test noise spectrum in the measured pulsating pressure amplitude logarithmic spectrum in sequence to obtain the sparse interference peak spectrum; Determine a pulsating pressure reference value, and obtain the pulsating pressure effective amplitude spectrum based on the pulsating pressure reference value and the pulsating pressure amplitude logarithmic spectrum after the total variation filtering; The smoothed pulsating pressure amplitude logarithmic spectrum is used as the measured pulsating pressure effective signal.

[0012] In a second aspect, the present invention further provides a system for eliminating a pulsating pressure interference peak in a pump source pipeline, comprising: An acquisition module, used for acquiring the logarithmic spectrum of the measured pulsating pressure amplitude in the pump source pipeline; A construction module, used to construct a filter vector of the logarithmic spectrum of the measured pulsating pressure amplitude by using a zero-phase filtering principle; A calculation module 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 by a main minimization algorithm, and obtain the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering; A smoothing module, used for obtaining a smoothed logarithmic spectrum of the pulsating pressure amplitude from the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering; The elimination module is used to separate the test noise spectrum and the sparse interference peak spectrum based on the measured pulsating pressure amplitude logarithmic spectrum and the pulsating pressure amplitude logarithmic spectrum after the total variation filtering, extract the pulsating pressure effective amplitude spectrum of the pulsating pressure amplitude logarithmic spectrum after the total variation filtering, and obtain the measured pulsating pressure effective signal.

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

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

[0015] The present invention provides a method and system for eliminating the interference peak of pulsating pressure in a pump source pipeline. A sparse interference peak filtering model for the logarithmic spectrum of the pulsating pressure amplitude is established by using the zero-phase filtering principle and the total variation filtering method. Then, an iterative calculation method for the effective signal of the pulsating pressure amplitude spectrum is formed through the main minimization algorithm. The effective signal of the pulsating pressure can be extracted from the amplitude spectrum containing the sparse interference peak through iterative calculation, thereby improving the accuracy of the pulsating pressure signal measurement results. In view of the problem of signal distortion in the results obtained by the total variation filtering model, the present invention proposes to use the moving least squares method to perform secondary processing on the pulsating pressure obtained by the filtering model, thereby solving the divergence problem of the extraction results of the filtering model and enhancing the reliability of the pulsating pressure filtering algorithm. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0017] Figure 1 It is a flow chart of a method for eliminating a pulsating pressure interference peak in a pump source pipeline provided by the present invention; Figure 2 It is a straight pipe structure of the pump source pipeline and a measurement point location diagram provided by the present invention; Figure 3 is a filtering decomposition result diagram provided by the present invention; Figure 4It is a comparison diagram of the filtering result provided by the present invention and the theoretical value of the pulsating pressure; Figure 5 It is a structural schematic diagram of a system for eliminating pulsating pressure interference peaks in a pump source pipeline provided by the present invention; Figure 6 It is a structural schematic diagram of the electronic device provided by the present invention. DETAILED DESCRIPTION

[0018] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention will be clearly and completely described below in conjunction with the drawings of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0019] Figure 1 is a flow chart of eliminating the pulsating pressure interference peak in the pump source pipeline provided by an embodiment of the present invention, such as Figure 1 As shown, including: Step 100: collecting the logarithmic spectrum of the measured pulsating pressure amplitude in the pump source pipeline; Step 200: constructing a filter vector of the logarithmic spectrum of the measured pulsating pressure amplitude using the zero-phase filtering principle; Step 300: Based on the constructed difference matrix, a total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude is established for the filtering vector, and the total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude is solved by a principal minimization algorithm to obtain the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering; Step 400: Obtain a smoothed pulsating pressure amplitude logarithmic spectrum from the pulsating pressure amplitude logarithmic spectrum after the total variation filtering; Step 500: Based on the measured pulsating pressure amplitude logarithmic spectrum and the pulsating pressure amplitude logarithmic spectrum after the total variation filtering, separate the test noise spectrum and the sparse interference peak spectrum, extract the pulsating pressure effective amplitude spectrum of the pulsating pressure amplitude logarithmic spectrum after the total variation filtering, and obtain the measured pulsating pressure effective signal.

[0020] In view of the problem that the pulsating pressure test results in current projects are affected by structural vibration and have interference peaks and inaccurate results, the embodiment of the present invention proposes a method for eliminating sparse interference peaks in the frequency domain of measured pulsating pressure for the pump source pipeline in ships. The 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 forms an iterative calculation method for the effective signal of the pulsating pressure amplitude spectrum through the main minimization algorithm. Through iterative calculation, the effective signal of the pulsating pressure can be extracted from the amplitude spectrum containing sparse interference peaks, thereby improving the accuracy of the pulsating pressure signal measurement results.

[0021] Based on the above embodiment, step 100 includes: A pulsating pressure sensor is arranged in the pump source pipeline to test the measured pulsating pressure signal at each designated position in the pipeline; A pulsating pressure correlation function time series vector corresponding to the measured pulsating pressure signal is obtained, and the pulsating pressure correlation function time series vector is subjected to Fourier transformation to obtain the measured pulsating pressure amplitude logarithmic spectrum.

[0022] Specifically, for the pump source pipeline in the actual project, the pulsating pressure signal at the typical position of the test pipeline is arranged by arranging the pulsating pressure sensor to solve the pulsating pressure signal at each position (denoted as p ) of the pulsating pressure correlation function time series vector , and perform Fourier transformation to obtain the logarithmic spectrum of the measured pulsating pressure amplitude .

[0023] by Figure 2 The straight pipe in the figure is used as an example to illustrate the specific implementation process of the present invention. The straight pipe has an outer diameter of 52 mm and an inner diameter of 49 mm, and a length of 1580 mm. The two ends of the curved pipe are fixed, and the interior is filled with water, and the flow rate of the water medium is about 10 m / s. An invasive 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 follows Figure 3 As shown in (a) in .

[0024] Based on the above embodiment, step 200 includes: Obtain the vector dimension of the logarithmic spectrum of the measured pulsating pressure amplitude N , and the filter cutoff frequency; The transfer coefficient is obtained according to the cosine function of the filter cutoff frequency, the first filter coefficient matrix is ​​determined by the vector dimension, and the second filter coefficient matrix is ​​determined by the vector dimension and the transfer coefficient, and the first filter coefficient matrix and the second filter coefficient matrix are both N × N dimensional band matrix; A filter vector is obtained from the inverse matrix of the first filter coefficient matrix, the second filter coefficient matrix, the transposed matrix of the second filter coefficient matrix, and the logarithmic spectrum of the measured pulsating pressure amplitude.

[0025] Specifically, in the embodiment of the present invention, the filter vector of the logarithmic spectrum of the pulsating pressure amplitude is constructed by using the zero-phase filtering principle. :

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

[0027]

[0028]

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

[0030] In the formula, represents the filter cutoff frequency, and in the embodiment of the present invention, it is recommended to take 0.01π.

[0031] Based on the above embodiment, step 300 includes: Based on vector dimension N Construct a difference matrix, the difference matrix includes ( N -1)× N The first-order difference matrix of dimension and ( N -2)× N dimensional second-order difference matrix; Iterating the measured logarithmic spectrum of the pulsating pressure amplitude to obtain an iterative logarithmic spectrum of the measured pulsating pressure amplitude; Initializing the iterative measured pulsating pressure amplitude logarithmic spectrum to obtain an initialized pulsating pressure amplitude logarithmic spectrum, wherein the initialized pulsating pressure amplitude logarithmic spectrum corresponds to the first step of the iteration; Calculate the penalty matrix of the product of the iterative measured pulsating pressure amplitude log spectrum and the difference matrix, the penalty matrix is N dimensional diagonal matrix, each diagonal element is obtained by the iterative measured pulsating pressure amplitude logarithmic spectrum and a preset fixed value; Determine a first penalty coefficient corresponding to the first-order difference matrix and a second penalty coefficient corresponding to the second-order difference matrix, and obtain an iterative matrix based on the first penalty coefficient, the second penalty coefficient, the difference matrix, a transposed matrix of the difference matrix, and the penalty matrix; Obtaining a filter vector of a current iteration step, and obtaining a logarithmic spectrum of a measured pulsating pressure amplitude of a next iteration step from the filter vector of the current iteration step, the iteration matrix, the first filter coefficient matrix, and the second filter coefficient matrix; A residual function is constructed by the iterative measured pulsating pressure amplitude logarithmic spectrum, the measured pulsating pressure amplitude logarithmic spectrum, the first filter coefficient matrix, the second filter coefficient matrix, the first penalty coefficient, the first-order difference matrix, the second penalty coefficient and the second-order difference matrix, and the residual function is calculated and iterated until a preset convergence condition is met to obtain the pulsating pressure amplitude logarithmic spectrum after total variation filtering.

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

[0033]

[0034] In the formula, D1 represents ( N -1)× N dimensional first-order difference matrix; D2 represents ( N -2)× N dimensional second-order difference matrix.

[0035] Then, the total variation filter model of the logarithmic spectrum of the pulsating pressure amplitude is established and solved by the following main minimization algorithm: (1) Initialization, input ,when hour, , represents the number of iteration steps, is a positive integer greater than or equal to 1. hour, Represents the iterative logarithmic spectrum of the measured pulsating pressure amplitude.

[0036] (2) Calculate the penalty matrix ( i =1, 2), where Simplified

[0037]

[0038] In the formula, for

[0039] In the formula, in the embodiment of the present invention Recommended .

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

[0041] In the formula, is the penalty coefficient. In the embodiment of the present invention, it is recommended , .

[0042] (4) Calculation :

[0043] (5) Calculate the residual function :

[0044] Here Simplified ; Continue the above iteration until the convergence condition is met, that is,

[0045] In the formula, is 0, ε As the convergence criterion, it is recommended to take .

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

[0047] Based on the above embodiment, step 400 includes: The logarithmic spectrum of the pulsating pressure amplitude after filtering is calculated using the moving least squares method. Perform fitting and smoothing to obtain the smoothed logarithmic spectrum of the pulsating pressure amplitude .

[0048] It should be noted that the fitting and smoothing processing used in the embodiment of the present invention is an existing algorithm, which is a geometric processing method that combines curve / surface fitting and smoothing optimization. It is 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, which will not be repeated here.

[0049] Based on the above embodiment, step 500 includes: Separating the pulsating pressure amplitude logarithmic spectrum after the total variation filtering from the measured pulsating pressure amplitude logarithmic spectrum, and combining the first filter coefficient matrix and the second filter coefficient matrix to obtain the test noise spectrum; Separating the pulsating pressure amplitude logarithmic spectrum after the total variation filtering and the test noise spectrum in the measured pulsating pressure amplitude logarithmic spectrum in sequence to obtain the sparse interference peak spectrum; Determine a pulsating pressure reference value, and obtain the pulsating pressure effective amplitude spectrum based on the pulsating pressure reference value and the pulsating pressure amplitude logarithmic spectrum after the total variation filtering; The smoothed pulsating pressure amplitude logarithmic spectrum is used as the measured pulsating pressure effective signal.

[0050] Specifically, the separated test noise spectrum in the embodiment of the present invention is: The obtained frequency domain sparse interference peak spectrum is , the logarithmic spectrum of the effective amplitude of the pulsating pressure is .

[0051] Through the above process, the influence of sparse interference peaks and test noise can be eliminated, and the required pulsating pressure effective amplitude spectrum y( ω )for

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

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

[0054] The system for eliminating the pulsating pressure interference peak in the pump source pipeline provided by the present invention is described below. The system for eliminating the pulsating pressure interference peak in the pump source pipeline described below and the method for eliminating the pulsating pressure interference peak in the pump source pipeline described above can be referenced to each other.

[0055] Figure 5 is a schematic diagram of the structure of a system for eliminating pulsating pressure interference peaks in a pump source pipeline provided by an embodiment of the present invention, such as Figure 5 As shown, it includes: a collection module 51, a construction module 52, a calculation module 53, a smoothing module 54 and an elimination module 55, wherein: The acquisition module 51 is used to acquire the measured logarithmic spectrum of the pulsating pressure amplitude in the pump source pipeline; the construction module 52 is used to construct the filter vector of the measured logarithmic spectrum of the 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 filter vector based on the constructed difference matrix, and solve the total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude by the principal minimization algorithm to obtain the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering; the smoothing module 54 is used to obtain the smoothed logarithmic spectrum of the pulsating pressure amplitude from the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering; the elimination module 55 is used to separate the test noise spectrum and the sparse interference peak spectrum based on the measured logarithmic spectrum of the pulsating pressure amplitude and the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering, extract the effective amplitude spectrum of the pulsating pressure of the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering, and obtain the measured effective signal of the pulsating pressure.

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

[0057] In addition, the logic instructions in the above-mentioned memory 630 can be implemented in the form of a software functional unit and can be stored in a computer-readable storage medium when it is sold or used as an independent product. Based on such an understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including a number of instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, etc. Various media that can store program codes.

[0058] On the other hand, the present invention also provides a computer program product, the computer program product includes a computer program, the computer program can be stored in 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 methods, the method comprising: collecting the measured pulsating pressure amplitude logarithmic spectrum in the pump source pipeline; constructing a filter vector of the measured pulsating pressure amplitude logarithmic spectrum using the zero-phase filtering principle; based on the constructed difference matrix, establishing a pulsating pressure amplitude for the filter vector A total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude is obtained by solving the total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude through the principal minimization algorithm to obtain the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering; a smoothed logarithmic spectrum of the pulsating pressure amplitude is obtained from the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering; a test noise spectrum and a sparse interference peak spectrum are separated based on the measured logarithmic spectrum of the pulsating pressure amplitude and the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering, the effective amplitude spectrum of the pulsating pressure of the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering is extracted, and the measured effective signal of the pulsating pressure is obtained.

[0059] On the other hand, the present invention also provides a non-transitory computer-readable storage medium having a computer program stored thereon, which is implemented when the computer program is executed by a processor to execute the method for eliminating the pulsating pressure interference peak in the pump source pipeline provided by the above methods, the method comprising: collecting the measured pulsating pressure amplitude logarithmic spectrum in the pump source pipeline; constructing a filter vector of the measured pulsating pressure amplitude logarithmic spectrum using the zero-phase filtering principle; based on the constructed difference matrix, establishing a pulsating pressure amplitude logarithmic spectrum total variation filtering model for the filter vector, The minimization algorithm solves the total variation filtering model of the pulsating pressure amplitude logarithmic spectrum to obtain the pulsating pressure amplitude logarithmic spectrum after the total variation filtering; the pulsating pressure amplitude logarithmic spectrum after the total variation filtering is used to obtain the smoothed pulsating pressure amplitude logarithmic spectrum; based on the measured pulsating pressure amplitude logarithmic spectrum and the pulsating pressure amplitude logarithmic spectrum after the total variation filtering, the test noise spectrum and the sparse interference peak spectrum are separated, and the pulsating pressure effective amplitude spectrum of the pulsating pressure amplitude logarithmic spectrum after the total variation filtering is extracted to obtain the measured pulsating pressure effective signal. The device embodiments described above are only schematic, wherein the units described as separated components may or may not be physically separated, and the components displayed as units 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 may be selected according to actual needs to achieve the purpose of the present embodiment. Ordinary technicians in this field can understand and implement it without creative labor.

[0060] Through the description of the above implementation methods, those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware. Based on this understanding, the above technical solution is essentially or the part that contributes to the prior art can be embodied in the form of a software product, and the computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a disk, an optical disk, etc., including a number of instructions for 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.

[0061] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for eliminating pulsating pressure interference peaks in a pump source pipeline, characterized in that: include: Collect the logarithmic spectrum of the measured pulsating pressure amplitude in the pump source pipeline; Constructing a filter vector of the logarithmic spectrum of the measured pulsating pressure amplitude by using the zero-phase filtering principle; Based on the constructed difference matrix, a total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude is established for the filtering vector, and the total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude is solved by the main minimization algorithm to obtain the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering; Obtaining a smoothed logarithmic spectrum of the pulsating pressure amplitude from the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering; Based on the measured pulsating pressure amplitude logarithmic spectrum and the pulsating pressure amplitude logarithmic spectrum after the total variation filtering, the test noise spectrum and the sparse interference peak spectrum are separated, and the pulsating pressure effective amplitude spectrum of the pulsating pressure amplitude logarithmic spectrum after the total variation filtering is extracted to obtain the measured pulsating pressure effective signal.

2. The method for eliminating the pulsating pressure interference peak in the pump source pipeline according to claim 1, characterized in that: Collect the measured pulsating pressure amplitude logarithmic spectrum in the pump source pipeline, including: A pulsating pressure sensor is arranged in the pump source pipeline to test the measured pulsating pressure signal at each designated position in the pipeline; A pulsating pressure correlation function time series vector corresponding to the measured pulsating pressure signal is obtained, and the pulsating pressure correlation function time series vector is subjected to Fourier transformation to obtain the measured pulsating pressure amplitude logarithmic spectrum.

3. The method for eliminating the pulsating pressure interference peak in the pump source pipeline according to claim 1, characterized in that: The filter vector of the logarithmic spectrum of the measured pulsating pressure amplitude is constructed by using the zero-phase filtering principle, including: Obtain the vector dimension of the logarithmic spectrum of the measured pulsating pressure amplitude N , and the filter cutoff frequency; The transfer coefficient is obtained according to the cosine function of the filter cutoff frequency, the first filter coefficient matrix is ​​determined by the vector dimension, and the second filter coefficient matrix is ​​determined by the vector dimension and the transfer coefficient, and the first filter coefficient matrix and the second filter coefficient matrix are both N × N dimensional band matrix; A filter vector is obtained from the inverse matrix of the first filter coefficient matrix, the second filter coefficient matrix, the transposed matrix of the second filter coefficient matrix, and the logarithmic spectrum of the measured pulsating pressure amplitude.

4. The method for eliminating the pulsating pressure interference peak in the pump source pipeline according to claim 3, characterized in that: Based on the constructed difference matrix, a total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude is established for the filtering vector, and the total variation filtering model of the logarithmic spectrum of the pulsating pressure amplitude is solved by the main minimization algorithm to obtain the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering, including: Based on vector dimension N Construct a difference matrix, the difference matrix includes ( N -1)× N The first-order difference matrix of dimension and ( N -2)× N dimensional second-order difference matrix; Iterating the measured logarithmic spectrum of the pulsating pressure amplitude to obtain an iterative logarithmic spectrum of the measured pulsating pressure amplitude; Initializing the iterative measured pulsating pressure amplitude logarithmic spectrum to obtain an initialized pulsating pressure amplitude logarithmic spectrum, wherein the initialized pulsating pressure amplitude logarithmic spectrum corresponds to the first step of the iteration; Calculate the penalty matrix of the product of the iterative measured pulsating pressure amplitude log spectrum and the difference matrix, the penalty matrix is N dimensional diagonal matrix, each diagonal element is obtained by the iterative measured pulsating pressure amplitude logarithmic spectrum and a preset fixed value; Determine a first penalty coefficient corresponding to the first-order difference matrix and a second penalty coefficient corresponding to the second-order difference matrix, and obtain an iterative matrix based on the first penalty coefficient, the second penalty coefficient, the difference matrix, a transposed matrix of the difference matrix, and the penalty matrix; Obtaining a filter vector of a current iteration step, and obtaining a logarithmic spectrum of a measured pulsating pressure amplitude of a next iteration step from the filter vector of the current iteration step, the iteration matrix, the first filter coefficient matrix, and the second filter coefficient matrix; A residual function is constructed by the iterative measured pulsating pressure amplitude logarithmic spectrum, the measured pulsating pressure amplitude logarithmic spectrum, the first filter coefficient matrix, the second filter coefficient matrix, the first penalty coefficient, the first-order difference matrix, the second penalty coefficient and the second-order difference matrix, and the residual function is calculated and iterated until a preset convergence condition is met to obtain the pulsating pressure amplitude logarithmic spectrum after total variation filtering.

5. The method for eliminating the pulsating pressure interference peak in the pump source pipeline according to claim 1, characterized in that: The smoothed pulsating pressure amplitude logarithmic spectrum is obtained from the pulsating pressure amplitude logarithmic spectrum after the total variation filtering, including: The moving least square method is used to fit and smooth the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering to obtain the smoothed logarithmic spectrum of the pulsating pressure amplitude.

6. The method for eliminating the pulsating pressure interference peak in the pump source pipeline according to claim 4, characterized in that: Based on the measured pulsating pressure amplitude logarithmic spectrum and the pulsating pressure amplitude logarithmic spectrum after the total variation filtering, the test noise spectrum and the sparse interference peak spectrum are separated, and the pulsating pressure effective amplitude spectrum of the pulsating pressure amplitude logarithmic spectrum after the total variation filtering is extracted to obtain the measured pulsating pressure effective signal, including: Separating the pulsating pressure amplitude logarithmic spectrum after the total variation filtering from the measured pulsating pressure amplitude logarithmic spectrum, and combining the first filter coefficient matrix and the second filter coefficient matrix to obtain the test noise spectrum; Separating the pulsating pressure amplitude logarithmic spectrum after the total variation filtering and the test noise spectrum in the measured pulsating pressure amplitude logarithmic spectrum in sequence to obtain the sparse interference peak spectrum; Determine a pulsating pressure reference value, and obtain the pulsating pressure effective amplitude spectrum based on the pulsating pressure reference value and the pulsating pressure amplitude logarithmic spectrum after the total variation filtering; The smoothed pulsating pressure amplitude logarithmic spectrum is used as the measured pulsating pressure effective signal.

7. A system for eliminating pulsating pressure interference peaks in a pump source pipeline, characterized in that: include: An acquisition module, used for acquiring the logarithmic spectrum of the measured pulsating pressure amplitude in the pump source pipeline; A construction module, used to construct a filter vector of the logarithmic spectrum of the measured pulsating pressure amplitude by using a zero-phase filtering principle; A calculation module 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 by a main minimization algorithm, and obtain the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering; A smoothing module, used for obtaining a smoothed logarithmic spectrum of the pulsating pressure amplitude from the logarithmic spectrum of the pulsating pressure amplitude after the total variation filtering; The elimination module is used to separate the test noise spectrum and the sparse interference peak spectrum based on the measured pulsating pressure amplitude logarithmic spectrum and the pulsating pressure amplitude logarithmic spectrum after the total variation filtering, extract the pulsating pressure effective amplitude spectrum of the pulsating pressure amplitude logarithmic spectrum after the total variation filtering, and obtain the measured pulsating pressure effective signal.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the program, the method for eliminating the pulsating pressure interference peak in the pump source pipeline as claimed in any one of claims 1 to 6 is implemented.

9. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method for eliminating the pulsating pressure interference peak in the pump source pipeline as claimed in any one of claims 1 to 6 is implemented.

10. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the method for eliminating the pulsating pressure interference peak in the pump source pipeline as claimed in any one of claims 1 to 6 is implemented.

Images