CN118010119B - Accurate electromagnetic flow measurement method with pressure compensation - Google Patents
Accurate electromagnetic flow measurement method with pressure compensation Download PDFInfo
- Publication number
- CN118010119B CN118010119B CN202410194019.8A CN202410194019A CN118010119B CN 118010119 B CN118010119 B CN 118010119B CN 202410194019 A CN202410194019 A CN 202410194019A CN 118010119 B CN118010119 B CN 118010119B
- Authority
- CN
- China
- Prior art keywords
- signal
- flow
- signal point
- point
- pressure
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000691 measurement method Methods 0.000 title claims abstract description 12
- 230000035772 mutation Effects 0.000 claims abstract description 61
- 238000001914 filtration Methods 0.000 claims abstract description 22
- 238000012937 correction Methods 0.000 claims description 30
- 238000000034 method Methods 0.000 claims description 17
- 230000003044 adaptive effect Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 abstract description 26
- 238000012545 processing Methods 0.000 description 6
- 239000007788 liquid Substances 0.000 description 5
- 239000012535 impurity Substances 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000012966 insertion method Methods 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000004886 process control Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 229910001220 stainless steel Inorganic materials 0.000 description 1
- 239000010935 stainless steel Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F15/00—Details of, or accessories for, apparatus of groups G01F1/00 - G01F13/00 insofar as such details or appliances are not adapted to particular types of such apparatus
- G01F15/02—Compensating or correcting for variations in pressure, density or temperature
- G01F15/022—Compensating or correcting for variations in pressure, density or temperature using electrical means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/56—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using electric or magnetic effects
- G01F1/58—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by using electric or magnetic effects by electromagnetic flowmeters
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Measuring Volume Flow (AREA)
Abstract
The invention relates to the technical field of electromagnetic flow measurement, in particular to an electromagnetic flow accurate measurement method with pressure compensation, which comprises the following steps: acquiring a pressure signal and a flow signal in a pipeline; setting a reference signal segment of a signal point; respectively analyzing the numerical distribution differences of the flow signal and the pressure signal in the reference signal section of the signal point, and constructing a first mutation degree of the flow signal and a first mutation degree of the pressure signal of the signal point; correcting the signal point flow signals according to the distribution characteristics of pressure changes and flow changes among all adjacent signal points in the reference signal section of the signal point; constructing a second degree of mutation of the signal point; and obtaining the self-adaptive time window size of the flow signal of each signal point, and obtaining the denoised flow signal by combining a Kalman filtering algorithm. The invention aims to track measurement noise in real time and measure electromagnetic flow more accurately.
Description
Technical Field
The application relates to the technical field of electromagnetic flow measurement, in particular to an electromagnetic flow accurate measurement method with pressure compensation.
Background
Electromagnetic flow measurement is based on Faraday's law of electromagnetic induction, and has certain requirement on fluid conductivity, and the electromagnetic flow measurement is generally divided into two parts, namely a sensor and a converter. The sensor coil is applied with exciting current to generate a magnetic field, fluid in a pipeline cuts magnetic lines of force to generate induced electromotive force, an induced electromotive force signal is processed and amplified, and flow is obtained based on an algorithm. The electromagnetic flow measuring instrument has the advantages of wide measuring range, small starting flow, high measuring precision, long service life and the like, and is widely applied to the fields of instruments and meters and process control.
The current electromagnetic flow measurement signal size is determined based on the current applied to the sensor exciting coil, and the common electromagnetic flow meter exciting current can only enable the signal to noise ratio to be larger at 100-250mA, so that the principle power consumption is larger; the induced electromotive voltage is small and is easy to be interfered by the outside, and a plurality of factors can influence signals, such as electromagnetic interference in the environment or impurities contained in liquid, and the impurities can cause signal jump when touching the sensor electrode, so that the measurement accuracy is influenced; on the other hand, instruments using electromagnetic flow measurement principles typically employ long-time filtering to signal interference, resulting in signal lag and slow response.
Disclosure of Invention
In order to solve the technical problems, the invention provides an electromagnetic flow accurate measurement method with pressure compensation, so as to solve the existing problems.
The invention relates to an electromagnetic flow accurate measurement method with pressure compensation, which adopts the following technical scheme:
An embodiment of the present invention provides a method for accurately measuring electromagnetic flow with pressure compensation, the method comprising the steps of:
acquiring a pressure signal and a flow signal in a pipeline;
Setting a reference signal section of each signal point; acquiring a first mutation degree of the flow signals of each signal point according to the numerical distribution difference of the flow signals of the reference signal section of each signal point; acquiring the first mutation degree of the pressure signal of each signal point by adopting a calculation method which is the same as the first mutation degree of the flow distribution of each signal point; acquiring a pressure variation correction factor of each signal point according to the distribution characteristics of pressure variation and flow variation among all adjacent signal points in the reference signal section of each signal point and the first mutation degree of the flow signal of each signal point; acquiring a variation correction factor of each signal point according to the pressure variation correction factor of each signal point and the pressure signal variation of each signal point; correcting the flow signals of the signal points according to the variable quantity correction factors of the signal points to obtain corrected flow values of the signal points; obtaining a second mutation degree of each signal point according to the difference characteristics of the first mutation degree of the flow signals of each signal point and the pressure signals and the difference characteristics of the flow signal distribution of each signal point and the reference signal section before and after correction;
Acquiring the self-adaptive time window size of the flow signal of each signal point according to the second mutation degree of each signal point; and obtaining the denoised flow signal according to the self-adaptive window size of each signal point and a Kalman filtering algorithm.
Preferably, the setting the reference signal of each signal point specifically includes:
centering on each signal point;
When the number of signal points in front of each signal point is greater than or equal to N and the number of signal points in back is greater than or equal to N-1, forming a reference signal section of each signal point by each signal point and each N signal points in front and back of each signal point, wherein N is a preset value.
Preferably, the first mutation degree of the flow signal of each signal point is obtained according to the numerical distribution difference of the flow signal of the reference signal segment of each signal point, specifically:
For each signal point;
Obtaining a maximum value, a minimum value and a mean value of flow signals in a reference signal section of each signal point; taking the absolute difference value of each signal point and the flow signal mean value as a first absolute difference value; acquiring the flow signal variation of each signal point;
Acquiring the average value of all the flow signal variation amounts in the reference signal section of each signal point; taking the absolute value of the difference between the flow signal variation and the mean value of each signal point as a second absolute value of the difference; calculating the ratio of the maximum value of the flow signal to the minimum value of the flow signal;
and taking the result of multiplying the first difference absolute value, the second difference absolute value and the ratio as a first mutation degree of the flow signal of each signal point.
Preferably, the flow signal variation of each signal point is specifically an absolute value of a flow signal difference between each signal point and a previous signal point.
Preferably, the obtaining the pressure variation correction factor of each signal point according to the distribution characteristics of the pressure variation and the flow variation between all adjacent signal points in the reference signal segment of each signal point and combining the first mutation degree of the flow signal of each signal point specifically includes:
Acquiring a u-th pressure signal variation delta Y x,u in a reference signal section of an x-th signal point by adopting a calculation method which is the same as the flow signal variation; calculating a sum FT xd,u of the first mutation degree of the flow signal of the first signal point corresponding to the variation of the u-th flow signal in the reference signal section of the x-th signal point, wherein the pressure variation correction factor YB x of the x-th signal point is expressed as follows:
Where U represents the number of flow signal variations, Δl x,u represents the U-th flow signal variation in the reference signal segment of the x-th signal point.
Preferably, the variation correction factor of each signal point is specifically a product of a pressure signal variation of each signal point and a pressure variation correction factor of each signal point.
Preferably, the flow value corrected by each signal point is specifically a sum of the flow signal of the previous signal point of each signal point and the variation correction factor of each signal point.
Preferably, the second mutation degree of each signal point is obtained according to the difference characteristics of the first mutation degree of the flow signal and the pressure signal of each signal point and the difference characteristics of the flow signal distribution of the reference signal section of each signal point before and after correction, the specific expression is:
Where ST x denotes a second degree of abrupt change of the xth signal point, FT x denotes a first degree of abrupt change of the xth signal point flow rate signal, YT x denotes a first degree of abrupt change of the xth signal point pressure signal, L x denotes the xth signal point flow rate signal, L x,y denotes a flow rate value after correction of the xth signal point, Representing the average of all flow signals in the reference signal segment of the xth signal point.
Preferably, the obtaining the adaptive time window size of the flow signal of each signal point according to the second mutation degree of each signal point specifically includes:
Presetting the initial time window size of each signal point; acquiring a normalized value of a second mutation degree of each signal point; calculating the product of the normalized value and the initial time window size; calculating a sum of the product and 0.5; the rounded value of the sum is taken as the adaptive window size of each signal point.
Preferably, the obtaining the flow signal after denoising according to the adaptive window size of each signal point and the kalman filtering algorithm specifically includes:
And taking all the signal point flow signals as the input of Kalman filtering, and combining the self-adaptive time window size of each signal point flow signal, wherein the output of the Kalman filtering algorithm is the flow signal after denoising.
The invention has at least the following beneficial effects:
the invention is mainly based on the original electromagnetic flow measurement principle, increases the collection of pressure signals in the pipeline, and can effectively solve the problems of small signal-to-noise ratio and impurity in the measured liquid under special working conditions by the signal processing of the analog circuit and the digital signal algorithm processing of the singlechip, thereby improving the reliability of instrument data.
Meanwhile, due to uncertainty of environmental conditions, measurement noise of different flow data points may be different in practical application, and the distribution characteristic of the measurement noise can be more accurately described by adaptively adjusting the time window size of each signal point flow signal for Kalman filtering, so that estimation errors are reduced; the method can track and update the measurement noise in real time, and timely reflect the change of the measurement noise so as to maintain the measurement quality, and has stronger robustness and applicability.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions and advantages of the prior art, the following description will briefly explain the drawings used in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are only some embodiments of the invention, and other drawings can be obtained according to the drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of an electromagnetic flow accurate measurement method with pressure compensation according to the present invention;
fig. 2 is a simplified diagram of an improved flow measurement device.
Fig. 3 is a diagram of the steps for obtaining the adaptive time window size.
Detailed Description
In order to further describe the technical means and effects adopted by the invention to achieve the preset aim, the following is a specific implementation, structure, characteristics and effects of an electromagnetic flow accurate measurement method with pressure compensation according to the invention, which are described in detail below with reference to the accompanying drawings and preferred embodiments. In the following description, different "one embodiment" or "another embodiment" means that the embodiments are not necessarily the same. Furthermore, the particular features, structures, or characteristics of one or more embodiments may be combined in any suitable manner.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
The following specifically describes a specific scheme of the electromagnetic flow accurate measurement method with pressure compensation provided by the invention with reference to the accompanying drawings.
The embodiment of the invention provides an electromagnetic flow accurate measurement method with pressure compensation.
Specifically, the following electromagnetic flow accurate measurement method with pressure compensation is provided, please refer to fig. 1, the method includes the following steps:
Step S001: pressure and flow signals are collected within the pipeline.
The common electromagnetic flow measuring device is divided into a flow sensor and a converter, the pressure sensor is additionally arranged on the basis of original hardware, the pressure sensor is protected by a stainless steel shell, pressure signals can be acquired in real time, the pressure sensor is electrically connected with the converter, pressure signal transmission is guaranteed, and a diffusion silicon sleeve can be used for the pressure sensor to ensure reliability. In the embodiment, pressure and flow signals are synchronously collected, the collection time is thirty minutes, the collection frequency is one second, and a converter provides exciting current and voltage for the flow and pressure sensors. Wherein a simplified diagram of the improved flow measurement device is shown in fig. 2.
The signals generated by the flow sensor and the pressure sensor are designed into a differential input channel formed by a high common mode rejection ratio amplifier through a signal amplifying circuit of a converter. And the filter circuit filters high-frequency interference and low-frequency differential interference signals generated by polarization voltage from the synthesized signals output by the front stage through a band-pass filter. And converting the digital signals into 24-bit sigma-delta ADC, and finally processing the digital signals through a singlechip software algorithm.
So far, the pressure signal and the flow signal in the pipeline at each moment are obtained.
Step S002: respectively analyzing the flow signal and pressure signal distribution characteristics of each signal point to obtain a first mutation degree of the flow signal and pressure signal of each signal point; correcting the flow signal of the signal point by combining the pressure signal; analyzing the difference of the first mutation degree of the flow signal and the pressure signal in the reference signal section of each signal point and the distribution characteristics of the difference of the flow signal before and after correction to construct a second mutation degree of the signal point; and carrying out self-adaptive adjustment on the time window size of Kalman filtering on the flow signal of each signal point in the signal section by combining the second mutation degree.
According to the embodiment, the change characteristics of the flow signal and the pressure signal of each time signal point are mainly analyzed, and the size of the time window of the flow signal of each signal point in the signal section is adaptively adjusted, so that more accurate measurement and denoising of the flow signal are realized. Firstly, analyzing the flow value of each signal point and the change characteristics of surrounding signal points to obtain a first mutation degree of the flow signal of each signal point.
Specifically, for each signal point; taking the signal point as the center, selecting N signal points forwards and N signal points backwards. The signal point and the front and rear N signal points are analyzed as reference signal segments of the signal point, where N in this embodiment is set to 50, and the practitioner can set the signal point according to the actual situation. If the number of signal points is less than N, the average value insertion method may be adopted, and the average value insertion method is a known technique, which is not described in the embodiment. Calculating a first mutation degree of flow signals of all signal points, wherein the expression is as follows:
Where FT x represents the first degree of abrupt change in the flow signal at the x-th signal point, L x represents the flow signal at the x-th signal point, Represents the average value of all flow signals in the reference signal section of the xth signal point, deltaL x represents the flow signal variation of the xth signal point and the xth-1 signal point, U represents the number of the flow signal variation, deltaL x,u represents the U-th flow signal variation in the reference signal section of the xth signal point,The maximum flow value and the minimum flow value in the reference signal section of the xth signal point are respectively represented, and the I represents an absolute value. Will beSave as the absolute value of the first differenceStored as a second absolute difference value.
Wherein the larger the difference between the value of each signal point and the average value of the flow signals of the surrounding signal points, i.e The greater the likelihood that the signal point will be mutated, the greater the first degree of mutation at that point; the mean value of all flow signal variation in the reference signal section representing the xth signal point is the greater the difference between the variation of the value of each flow signal and the mean value of all flow signal variation in the reference signal section, i.e The greater the likelihood of an abnormal mutation at the signal point, the greater the first degree of mutation at that point; the greater the degree of polarization representing the magnitude of the signal points within the signal segment, the greater the magnitude of the signal points within the signal segment corresponding to each signal point, and the greater the likelihood of abrupt changes in the signal segment in which the signal point is located, the greater the first degree of abrupt changes in the signal point.
Because the instrument of the electromagnetic flow measurement principle usually adopts a long-time filtering mode facing signal interference, the signal hysteresis response speed is slow, and the measurement value error of the flow signal is larger. Therefore, pressure signals are required to be introduced to compensate the flow signals to a certain extent, the flow value of each signal point is corrected, and the expression is as follows:
Lx,y＝Lx-1+YBx×ΔYx
Wherein L x,y denotes a flow rate value corrected for the x-th signal point, L x-1 denotes a flow rate value of the x-1 th signal point, YB x denotes a pressure change amount correction factor of the x-th signal point, U denotes the number of flow rate signal change amounts, Δl x,u、ΔYx,u denotes a U-th flow rate signal change amount and a pressure signal change amount in a reference signal segment of the x-th signal point, FT xd,u denotes a sum value of a first mutation degree of a U-th signal point flow rate signal and a preceding signal point flow rate signal in a reference signal segment of the x-th signal point, Δy x denotes a pressure signal change amount of the x-th signal point, exp () denotes an exponential function based on a natural constant. Will be Stored as a variation correction factor. The flow signal variation of each signal point is the absolute value of the difference value of the flow signals of each signal point and the previous signal point; and obtaining the pressure signal variation of each signal point by adopting the same calculation method as the flow signal variation.
The ratio of the flow change to the pressure change of each signal point in the reference signal segment representing the xth signal point is corrected, and the flow value of each signal point needs to be weighted average. When the weighted average is carried out, the corresponding weight of each ratio is determined by the first mutation degree of the flow signals in the corresponding adjacent signal points, if the sum of the first mutation degrees of the corresponding two signal points is larger, the group of adjacent flow signals is abnormal, the weight is thatThe smaller should be. By combining DeltaY x withThe flow value variation after the correction of the xth signal point is obtained after multiplication compared with the last signal point, and the flow value after the correction of the xth signal point is obtained after addition with L x-1.
And then analyzing the correlation of the changes between the two to obtain the second mutation degree of each signal point. First, the first mutation degree of the pressure signal of each signal point is obtained by adopting a calculation method which is the same as the first mutation degree of the flow signal of each signal point. As is known from bernoulli's principle, when a liquid flows through a pipe, the velocity of the liquid increases with increasing flow because the cross-sectional area of the pipe interior is constant. Then, analyzing the first mutation degree of the signal point flow signal and the first mutation degree of the pressure signal at the same moment to obtain the second mutation degree of each signal point, wherein the expression is as follows:
Where ST x denotes a second degree of abrupt change of the xth signal point, FT x denotes a first degree of abrupt change of the xth signal point flow rate signal, YT x denotes a first degree of abrupt change of the xth signal point pressure signal, L x denotes the xth signal point flow rate signal, L x,y denotes a flow rate value after correction of the xth signal point, Representing the average of all flow signals in the reference signal segment of the xth signal point.
The |ft x-YTx | represents the difference between the first mutation degree of the flow signal and the pressure signal at the same time of the xth signal point, and the greater the difference between the first mutation degree of each signal point and the first mutation degree of the pressure signal at the same time corresponding to the xth signal point, the greater the possibility that the signal point belongs to a noise point generated by environmental interference, the greater the second mutation degree of the signal point; the larger the difference between the flow value of each signal point compensated by the pressure data and the average value of the flow signal of the signal section where the signal point is located and the difference between the flow value of the signal point and the average value of the flow signal of the signal section where the signal point is located, the larger the difference between the corrected flow and the average value of the flow signal of the signal section is, the more abnormal the corrected flow value is, and the second mutation degree corresponding to the signal point is larger.
Thus, the second mutation degree of each signal point is obtained.
Step S003: and obtaining the time window size of each signal point for Kalman filtering according to the second mutation degree of each signal point, and denoising the flow signals of the signal points according to a Kalman filtering algorithm.
The second mutation degree of each signal point is obtained through the steps. Then the adjustment of the time window size for each signal point in the entire flow signal segment is based on the second degree of abrupt change of each signal point. Starting from a first signal point flow signal of the whole signal section, taking a reference signal section of each signal point as an initial time window for carrying out Kalman filtering, and acquiring the self-adaptive time window size of the flow signal of each signal point by combining a second mutation degree, wherein the expression is as follows:
Cx＝Round([Co×norm(STx)]+0.5)
Where C x represents the time window size of the x-th signal point traffic signal for kalman filtering, C o represents the initial time window size, ST x represents the second degree of abrupt change of the x-th signal point, norm () represents the normalization function, and Round () represents the rounding function. The step of obtaining the adaptive time window size is shown in fig. 3.
The greater the second degree of abrupt change corresponding to the signal points, the greater the likelihood that the current data point to be filtered belongs to noise data, so a smaller time window is required to respond more quickly to changes in system state, thereby reducing the cumulative effect of noise, and so a reduced kalman filter time window size is required.
And carrying out data denoising on the flow signals by using an optimized Kalman filtering algorithm, taking all signal point flow signals as Kalman filtering inputs, and combining the self-adaptive time window size of the flow signals of all signal points, wherein the output of the Kalman filtering algorithm is the flow signals after denoising. It should be noted that, the kalman filtering algorithm is a known technology, and will not be described in detail in this embodiment.
To this end, measurement of the electromagnetic flow based on pressure compensation is completed.
In summary, the embodiment of the invention mainly increases the collection of the pressure signal in the pipeline based on the original electromagnetic flow measurement principle, and can effectively solve the problems of small signal-to-noise ratio and impurity influence in the measured liquid under special working conditions and improve the reliability of instrument data through the signal processing of the analog circuit and the digital signal algorithm processing of the singlechip.
Meanwhile, due to uncertainty of environmental conditions, measurement noise of different flow data points may be different in practical application, and the distribution characteristic of the measurement noise can be more accurately described by adaptively adjusting the time window size of each signal point flow signal for Kalman filtering, so that estimation errors are reduced; the method can track and update the measurement noise in real time, and timely reflect the change of the measurement noise so as to maintain the measurement quality, and has stronger robustness and applicability.
It should be noted that: the sequence of the embodiments of the present invention is only for description, and does not represent the advantages and disadvantages of the embodiments. And the foregoing description has been directed to specific embodiments of this specification. In addition, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In some embodiments, multitasking and parallel processing are also possible or may be advantageous.
In this specification, each embodiment is described in a progressive manner, and the same or similar parts of each embodiment are referred to each other, and each embodiment mainly describes differences from other embodiments.
The above embodiments are only for illustrating the technical solution of the present application, and not for limiting the same; the technical solutions described in the foregoing embodiments are modified or some of the technical features are replaced equivalently, so that the essence of the corresponding technical solutions does not deviate from the scope of the technical solutions of the embodiments of the present application, and all the technical solutions are included in the protection scope of the present application.
Claims (5)
1. The electromagnetic flow accurate measurement method with pressure compensation is characterized by comprising the following steps of:
acquiring a pressure signal and a flow signal in a pipeline;
Setting a reference signal section of each signal point; acquiring a first mutation degree of the flow signals of each signal point according to the numerical distribution difference of the flow signals of the reference signal section of each signal point; acquiring the first mutation degree of the pressure signal of each signal point by adopting a calculation method which is the same as the first mutation degree of the flow distribution of each signal point; acquiring a pressure variation correction factor of each signal point according to the distribution characteristics of pressure variation and flow variation among all adjacent signal points in the reference signal section of each signal point and the first mutation degree of the flow signal of each signal point; acquiring a variation correction factor of each signal point according to the pressure variation correction factor of each signal point and the pressure signal variation of each signal point; correcting the flow signals of the signal points according to the variable quantity correction factors of the signal points to obtain corrected flow values of the signal points; obtaining a second mutation degree of each signal point according to the difference characteristics of the first mutation degree of the flow signals of each signal point and the pressure signals and the difference characteristics of the flow signal distribution of each signal point and the reference signal section before and after correction;
Acquiring the self-adaptive time window size of the flow signal of each signal point according to the second mutation degree of each signal point; acquiring a denoised flow signal according to the self-adaptive window size of each signal point and a Kalman filtering algorithm;
the first mutation degree of the flow signal of each signal point is obtained according to the numerical distribution difference of the flow signal of the reference signal section of each signal point, specifically:
For each signal point;
Obtaining a maximum value, a minimum value and a mean value of flow signals in a reference signal section of each signal point; taking the absolute difference value of each signal point and the flow signal mean value as a first absolute difference value; acquiring the flow signal variation of each signal point;
Acquiring the average value of all the flow signal variation amounts in the reference signal section of each signal point; taking the absolute value of the difference between the flow signal variation and the mean value of each signal point as a second absolute value of the difference; calculating the ratio of the maximum value of the flow signal to the minimum value of the flow signal;
taking the result of multiplying the first difference absolute value, the second difference absolute value and the ratio as a first mutation degree of the flow signal of each signal point;
the flow signal variation of each signal point is specifically the absolute value of the flow signal difference between each signal point and the previous signal point;
the method comprises the steps of obtaining a pressure change quantity correction factor of each signal point according to the distribution characteristics of pressure change and flow change between all adjacent signal points in a reference signal section of each signal point and combining the first mutation degree of the flow signal of each signal point, wherein the correction factor comprises the following specific steps:
acquiring the variation of the ith pressure signal in the reference signal segment of the xth signal point by adopting the same calculation method as the variation of the flow signal ; Calculating the sum value of the first mutation degree of the u-th flow signal variation quantity corresponding to the flow signals of two adjacent signal points in the reference signal section of the x-th signal pointPressure change correction factor for the xth signal pointThe expression is:
；
In the method, in the process of the invention, The number of the flow rate signal variation amounts is represented,A variation of the u-th flow signal in the reference signal section of the x-th signal point is represented;
The change amount correction factor of each signal point is specifically the product of the pressure signal change amount of each signal point and the pressure change amount correction factor of each signal point;
The second mutation degree of each signal point is obtained according to the difference characteristics of the first mutation degree of the flow signal of each signal point and the pressure signal and the difference characteristics of the flow signal distribution of each signal point and the reference signal section before and after correction, the specific expression is:
；
In the method, in the process of the invention, A second degree of mutation representing the xth signal point,A first degree of abrupt change in the x-th signal point flow signal is represented,A first degree of abrupt change in the pressure signal representing the x-th signal point,Representing the x-th signal point traffic signal,Indicating the flow value corrected for the xth signal point,Representing the average of all flow signals in the reference signal segment of the xth signal point.
2. The method for precisely measuring the electromagnetic flow with pressure compensation according to claim 1, wherein the setting of the reference signals of each signal point is specifically as follows:
centering on each signal point;
When the number of signal points in front of each signal point is greater than or equal to N and the number of signal points in back is greater than or equal to N-1, forming a reference signal section of each signal point by each signal point and each N signal points in front and back of each signal point, wherein N is a preset value.
3. The method for precisely measuring the electromagnetic flow with pressure compensation according to claim 1, wherein the flow value corrected by each signal point is specifically the sum of the flow signal of the previous signal point of each signal point and the variation correction factor of each signal point.
4. The method for precisely measuring the electromagnetic flow with pressure compensation according to claim 1, wherein the adaptive time window size of the flow signal of each signal point is obtained according to the second mutation degree of each signal point, specifically:
Presetting the initial time window size of each signal point; acquiring a normalized value of a second mutation degree of each signal point; calculating the product of the normalized value and the initial time window size; calculating a sum of the product and 0.5; the rounded value of the sum is taken as the adaptive window size of each signal point.
5. The method for precisely measuring the electromagnetic flow with pressure compensation according to claim 1, wherein the method for precisely measuring the electromagnetic flow with pressure compensation is characterized in that the flow signal after denoising is obtained according to the self-adaptive window size of each signal point and a Kalman filtering algorithm, specifically comprises the following steps:
And taking all the signal point flow signals as the input of Kalman filtering, and combining the self-adaptive time window size of each signal point flow signal, wherein the output of the Kalman filtering algorithm is the flow signal after denoising.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410194019.8A CN118010119B (en) | 2024-02-21 | 2024-02-21 | Accurate electromagnetic flow measurement method with pressure compensation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410194019.8A CN118010119B (en) | 2024-02-21 | 2024-02-21 | Accurate electromagnetic flow measurement method with pressure compensation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118010119A CN118010119A (en) | 2024-05-10 |
| CN118010119B true CN118010119B (en) | 2024-10-11 |
Family
ID=90944071
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410194019.8A Active CN118010119B (en) | 2024-02-21 | 2024-02-21 | Accurate electromagnetic flow measurement method with pressure compensation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118010119B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112923985A (en) * | 2021-01-22 | 2021-06-08 | 成都秦川物联网科技股份有限公司 | Electronic diaphragm gas meter adopting integrating instrument mode |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7165464B2 (en) * | 2002-11-15 | 2007-01-23 | Cidra Corporation | Apparatus and method for providing a flow measurement compensated for entrained gas |
| US9151649B2 (en) * | 2012-10-19 | 2015-10-06 | Daniel Measurement And Control, Inc. | Ultrasonic flow metering system with an upstream pressure transducer |
| JP6556709B2 (en) * | 2013-06-28 | 2019-08-07 | プロフタガレン アクチエボラグＰｒｏｖｔａｇａｒｅｎ Ａｂ | How to check the normal functioning of sampling equipment |
| CN204575150U (en) * | 2015-04-28 | 2015-08-19 | 江苏华海测控技术有限公司 | The one-piece type electromagnetic flowmeter of temperature and pressure |
| US10330511B2 (en) * | 2017-06-22 | 2019-06-25 | Saudi Arabian Oil Company | Alternating magnetic field flow meters |
| EP4097427A1 (en) * | 2020-01-31 | 2022-12-07 | Micro Motion, Inc. | Method of correcting flow meter variable |
| CN113155210A (en) * | 2021-03-30 | 2021-07-23 | 合肥工业大学 | Electromagnetic flow measurement system signal processing method based on improved Kalman filtering |
| CN216925676U (en) * | 2021-12-20 | 2022-07-08 | 深圳市拓安信计控仪表有限公司 | Electromagnetic water meter flow and built-in integrated pressure measuring device |
-
2024
- 2024-02-21 CN CN202410194019.8A patent/CN118010119B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112923985A (en) * | 2021-01-22 | 2021-06-08 | 成都秦川物联网科技股份有限公司 | Electronic diaphragm gas meter adopting integrating instrument mode |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118010119A (en) | 2024-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Elster et al. | Analysis of dynamic measurements and determination of time-dependent measurement uncertainty using a second-order model | |
| US5827963A (en) | System and method for determining a density of a fluid | |
| CN108534852A (en) | A kind of electromagnetic flowmeter and detection method based on dynamic detection | |
| CN107478402A (en) | A kind of low-frequency vibration displacement detection system and method | |
| CN117171516B (en) | Data optimization correction method for X-ray thickness gauge | |
| CN116625444B (en) | A method for characteristic wave self-adaptation and flow correction of ultrasonic water meter | |
| CA2456855C (en) | Method and apparatus for reducing skew in a real-time centroid calculation | |
| CN113865659B (en) | High-precision household ultrasonic water meter anti-interference filtering method and equipment | |
| CN102986250A (en) | System and method for measuring and validating the occlusion effect of a hearing aid user | |
| CN118010119B (en) | Accurate electromagnetic flow measurement method with pressure compensation | |
| CN115824331A (en) | Low-power consumption measuring method suitable for ultrasonic water meter | |
| CN108255108B (en) | A MEMS sensor controller circuit | |
| CN115773793B (en) | Dynamic regulation method for signal amplitude of ultrasonic water meter | |
| CN114459338B (en) | Adaptive regulation and control system and method for depth sensing signals of underwater vehicle | |
| US11929757B2 (en) | Propagation delay compensation and interpolation filter | |
| CN117928528B (en) | Ship heave measurement method based on adaptive time-delay-free complementary bandpass filter | |
| CN118051724A (en) | A method, device and medium for measuring open channel flow based on adaptive LMS filtering | |
| CN110987875A (en) | Water turbidity detection method for spatial fusion based on support degree and confidence degree | |
| CN112073035A (en) | Self-adaptive filtering method for PSD (phase-sensitive Detector) signals | |
| CN108008205B (en) | Self-adaptive signal conditioning circuit and method for electromagnetic field detection | |
| CN114063807A (en) | Proximity sensor with improved environmental compensation performance and environmental compensation method therein | |
| US10948552B2 (en) | Magnetometer | |
| RU99623U1 (en) | ADAPTIVE SIGNAL FILTER | |
| JPS6383614A (en) | Electromagnetic flow meter | |
| CN113192659B (en) | Nuclear reactor reactivity measurement method based on primary neutron signal average sampling |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |