CN117214461A - Suspension wing type flow velocity meter calibration method based on fusion algorithm - Google Patents
Suspension wing type flow velocity meter calibration method based on fusion algorithm Download PDFInfo
- Publication number
- CN117214461A CN117214461A CN202311111714.5A CN202311111714A CN117214461A CN 117214461 A CN117214461 A CN 117214461A CN 202311111714 A CN202311111714 A CN 202311111714A CN 117214461 A CN117214461 A CN 117214461A
- Authority
- CN
- China
- Prior art keywords
- flow velocity
- meter
- wing type
- angle
- data
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 67
- 230000004927 fusion Effects 0.000 title claims abstract description 17
- 239000000725 suspension Substances 0.000 title claims description 9
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 238000001914 filtration Methods 0.000 claims abstract description 4
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Measuring Volume Flow (AREA)
Abstract
The invention discloses a suspension-wing type flow rate meter calibration method based on a fusion algorithm, which comprises the steps of firstly, performing a median method anti-shake filtering on angle S data actually measured by a suspension-wing type flow rate meter; then, the relation between the flow velocity V data measured by the standard flow velocity meter is analyzed and error calculation is carried out; finally, performing data fitting by adopting a nonlinear least square method to obtain a calibration model corresponding to the S-V curve relation expression between the data; and adjusting the flow velocity of the test point and correcting errors through a calibration model corresponding to the S-V curve relation expression. According to the invention, a calibration model is designed by taking a median method and a nonlinear least square method as fusion algorithms, so that the calibration with higher precision and smaller error on the errors of test point data is realized.
Description
Technical Field
The invention belongs to the technical field of water flow measurement, and particularly relates to a suspension wing type flow meter calibration method based on a fusion algorithm.
Background
With the continuous updating of river flowmeters, higher requirements are also put on the calibration and correction of the measurement data of the river flowmeters. The most common handheld flow meter generally calculates the flow rate based on the number of revolutions of the rotor driven by the water flow. However, in the actual use process, the handheld flow rate meter is not suitable for long-time work due to the complex environment in river water, so that the suspended wing type flow rate meter appears, and in the actual measurement, certain parameter calibration is usually carried out on the measured parameters, namely, a plurality of datum reference points are introduced in the measurement use process.
In order to calibrate the measurement parameters, the traditional calibration method is to perform segmentation data processing on the measured data points (including angle values) and the actual flow rate, solve the slope and offset constants of each adjacent measurement parameter point, analyze and screen according to the slope values of the adjacent measurement parameters, and select a more appropriate reference data point. The traditional calibration method is based on a measuring parameter change rate model, so that the method is only suitable for the working condition that the fluctuation of measured data points is not very large, the angle value and the flow velocity value basically meet the linear relation, the working condition that the fluctuation of the measured data points is very large can not be approximated to the linear relation, the error obtained by calculating the error of the reference point selected after the processing of the method is very large, and the method can not realize automatic correction of the larger error, so that the error correction process becomes complex and complicated, and the method is difficult to directly popularize and apply.
Disclosure of Invention
The invention aims to: in order to overcome the defects in the prior art, the invention provides a suspension wing type flow velocity meter calibration method based on a fusion algorithm, which uses a median method and a nonlinear least square method as a fusion algorithm to design a calibration model, so as to realize the calibration of higher precision and smaller error of test point data.
The technical scheme is as follows: in order to achieve the above purpose, the invention provides a suspension wing type flowmeter calibration method based on a fusion algorithm, wherein the fusion algorithm is composed of a median method and a nonlinear least square method, and the specific steps are as follows:
step S1: firstly, performing a median method jitter elimination filtering on the angle S data actually measured by the suspended wing type flow velocity meter;
step S2: then, the relation between the flow velocity V data measured by the standard flow velocity meter is analyzed and error calculation is carried out;
step S3: finally, performing data fitting by adopting a nonlinear least square method to obtain a calibration model corresponding to the S-V curve relation expression between the data;
step S4: and adjusting the flow velocity of the test point and correcting errors through a calibration model corresponding to the S-V curve relation expression.
Further, the data source based on the flow velocity meter is composed of a plurality of groups of standard data point sets actually measured by the standard flow velocity meter, therefore, all standard data point sets measured by the suspended wing type flow velocity meter are respectively composed of a flow velocity V and an angle S, and a linear relation is satisfied between each adjacent standard data point set:
V real world =K i ×S Angle of +C
Wherein: k represents a linear relationship coefficient, K is according to S Angle of Determining the numerical range of (2); k (K) i Representing the slopes of two adjacent standard data point sets; c represents the offset of the overall linear relation formed by a plurality of groups of standard data point sets and is an offset constant; v (V) Real world Representing a flow velocity value actually measured according to the flow velocity meter; s is S Angle of Indicating the rotor angle.
Further, K i The expression of (2) is as follows:
wherein: (S) i ,V i ) Representing sets of standard data points measured by the flow meter, i=1, 2,3, … …, n.
Further, for any group of test point sets, the angle values S in a period of time are sequenced and then median values are removed, then a nonlinear least square method is adopted for fitting, and a final curve function M-degree polynomial is set:
wherein: x is a univariate input, ω 0 ,ω 1 ,…,ω M Is M+1 parameters;
the least squares method using square as the loss function includes:
for omega in the formula j Solving the bias guide and making the bias guide be 0:
the method can obtain:
fitting polynomial coefficientsSolving a system of linear equations:
and (3) calculating:
substituting the calculation result into a linear equation set for solving, and obtaining the S-V fitting relation.
The beneficial effects are that: according to the invention, a calibration model is designed by taking a median method and a nonlinear least square method as fusion algorithms, so that the calibration with higher precision and smaller error on the errors of test point data is realized; in addition, the method of the invention not only can automatically correct errors, but also can automatically output the datum point fitting curve, meets the error requirement of the test point data set, and has good application value.
Drawings
FIG. 1 is a schematic diagram of a standard data point set curve fitted by a linear fitting method and a least square fitting method;
FIG. 2 is a schematic diagram comparing respective curves of a linear fit and a least squares fit to a flow line.
Detailed Description
The invention will be further described with reference to the accompanying drawings.
A suspension wing type flowmeter calibration method based on a fusion algorithm comprises a median method and a nonlinear least square method, and the method comprises the following specific steps:
step S1: firstly, performing a median method jitter elimination filtering on the angle S data actually measured by the suspended wing type flow velocity meter;
step S2: then, the relation between the flow velocity V data measured by the standard flow velocity meter is analyzed and error calculation is carried out;
step S3: finally, performing data fitting by adopting a nonlinear least square method to obtain a calibration model corresponding to the S-V curve relation expression between the data;
step S4: and adjusting the flow velocity of the test point and correcting errors through a calibration model corresponding to the S-V curve relation expression.
The data source based on the flow velocity meter is composed of a plurality of groups of standard data point sets actually measured by the standard flow velocity meter, therefore, all standard data point sets measured by the suspended wing type flow velocity meter are respectively composed of a flow velocity V and an angle S, and the linear relation is satisfied between every two adjacent standard data point sets:
V real world =K i ×S Angle of +C
Wherein: k represents linearityRelation coefficient, K is according to S Angle of Determining the numerical range of (2); k (K) i Representing the slopes of two adjacent standard data point sets; c represents the offset of the overall linear relation formed by a plurality of groups of standard data point sets and is an offset constant; v (V) Real world Representing a flow velocity value actually measured according to the flow velocity meter; s is S Angle of Indicating the rotor angle.
K i The expression of (2) is as follows:
wherein: (S) i ,V i ) Representing sets of standard data points measured by the flow meter, i=1, 2,3, … …, n.
Because the relation between the flow velocity V and the angle S is nonlinear, and the values of V and S are discrete point sets, for any group of test point sets, the angle values S in a period of time are sequenced and then the median value is removed, then a nonlinear least square method is adopted for fitting, and a final curve function M-degree polynomial is set:
wherein: x is a univariate input, ω 0 ,ω 1 ,…,ω M Is M+1 parameters;
the least squares method using square as the loss function includes:
for omega in the formula j Solving the bias guide and making the bias guide be 0:
the method can obtain:
fitting polynomial coefficientsSolving a system of linear equations:
and (3) calculating:
substituting the calculation result into a linear equation set for solving, and obtaining the S-V fitting relation.
According to the invention, a calibration model is designed by taking a median method and a nonlinear least square method as fusion algorithms, so that the calibration with higher precision and smaller error on the errors of test point data is realized; in addition, the method of the invention not only can automatically correct errors, but also can automatically output the datum point fitting curve, meets the error requirement of the test point data set, and has good application value.
The equation for the least squares linear fit is: y=a×x+b. The principle of the nonlinear fitting of the least square method is as follows: in the nonlinear least squares fitting, the actual meaning is to find a function y=f (x) that is closest to all data points under a certain criterion, and to apply a certain basis function form, such as: y=a0+a1/x, where y is the coefficient a= (a) 0 ,a 1 ,a 2 ...a n-1 ,a n ) T For the above-mentioned basis function model, the nonlinear function of (a)Substituting the model, the model can become +>Calculate +.f. from the transform of the raw data (xi, yi) (i=1, 2,3, …, n)>(i=1, 2,3, …, n), a0 and a1 can be obtained by linear fitting, thereby obtaining a nonlinear fitting relation of y=a0+a1/x.
In order to verify the feasibility based on the median method and the least square method, the invention utilizes standard flow velocity equipment and collects the angle data of the suspended wing type flow velocity meter, and the method comprises the following steps:
and carrying out linear processing on the data point set to obtain a standard data point set S-V fitting relation formula: y=0.01485363x+0.06462, as shown in curve B of fig. 1, it can be found that at an actual flow rate of 0.698m/s, the calculated flow rate is 0.8m/s, and the error is 14.6% at the maximum.
And performing least square fitting on the data point set to obtain a standard data point set S-V fitting relation formula: y=6.9×10 -8 x 4 +4.452×10 -7 x 3 -5.02646×10 -4 x 2 +2.868×10 -2 x+0.0272, substituting the angle into the relation to obtain a C curve as shown in FIG. 1, the C curve is more consistent with the actual flow rate than the linear processing, and when the maximum error point is the flow rate of 0.6m/s, the actual flow rate is 0.572m/s, and the error is 4.7%.
In the practical application process, 681 angle data collected by the suspended wing type flow rate meter are processed by a median method, and then flow rate calculation is performed by using a relation formula fitted by a least square method, as shown in fig. 2: by comparison, the flow velocity mutation and the amplitude are obviously reduced after the median method is combined with the least square nonlinear fitting flow velocity.
The foregoing is only a preferred embodiment of the invention, it being noted that: it will be apparent to those skilled in the art that various modifications and adaptations can be made without departing from the principles of the present invention, and such modifications and adaptations are intended to be comprehended within the scope of the invention.
Claims (4)
1. A suspension wing type flow velocity meter calibration method based on a fusion algorithm is characterized in that: the fusion algorithm consists of a median method and a nonlinear least square method, and comprises the following specific steps:
step S1: firstly, performing a median method jitter elimination filtering on the angle S data actually measured by the suspended wing type flow velocity meter;
step S2: then, the relation between the flow velocity V data measured by the standard flow velocity meter is analyzed and error calculation is carried out;
step S3: finally, performing data fitting by adopting a nonlinear least square method to obtain a calibration model corresponding to the S-V curve relation expression between the data;
step S4: and adjusting the flow velocity of the test point and correcting errors through a calibration model corresponding to the S-V curve relation expression.
2. The fusion algorithm-based suspension wing type flow meter calibration method as set forth in claim 1, wherein the method comprises the following steps: the data source based on the flow velocity meter is composed of a plurality of groups of standard data point sets actually measured by the standard flow velocity meter, therefore, all standard data point sets measured by the suspended wing type flow velocity meter are respectively composed of a flow velocity V and an angle S, and the linear relation is satisfied between every two adjacent standard data point sets:
V real world =K i ×S Angle of +C
Wherein: k represents a linear relationship coefficient, K is according to S Angle of Determining the numerical range of (2); k (K) i Representing the slopes of two adjacent standard data point sets; c represents the offset of the overall linear relation formed by a plurality of groups of standard data point sets and is an offset constant; v (V) Real world Representing a flow velocity value actually measured according to the flow velocity meter; s is S Angle of Indicating the rotor angle.
3. The fusion algorithm-based suspension wing type flow meter calibration method as set forth in claim 2, wherein the method comprises the following steps: k (K) i The expression of (2) is as follows:
wherein: (S) i ,V i ) Represents a plurality of sets of standard data points measured by the flow meter, i=1, 2, 3.
4. A method for calibrating a suspended wing type flow rate meter based on a fusion algorithm according to claim 3, wherein the method comprises the following steps: for any group of test point sets, angle values S in a period of time are sequenced and then median values are removed, then a nonlinear least square method is adopted for fitting, and a final curve function M-degree polynomial is set:
wherein: x is a univariate input, ω 0 ,ω 1 ,...,ω M Is M+1 parameters;
the least squares method using square as the loss function includes:
for omega in the formula j Solving the bias guide and making the bias guide be 0:
the method can obtain:
fitting polynomial coefficientsSolving a system of linear equations:
and (3) calculating:
substituting the calculation result into a linear equation set for solving, and obtaining the S-V fitting relation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311111714.5A CN117214461B (en) | 2023-08-31 | 2023-08-31 | Suspension wing type flow velocity meter calibration method based on fusion algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311111714.5A CN117214461B (en) | 2023-08-31 | 2023-08-31 | Suspension wing type flow velocity meter calibration method based on fusion algorithm |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117214461A true CN117214461A (en) | 2023-12-12 |
CN117214461B CN117214461B (en) | 2024-04-09 |
Family
ID=89047221
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311111714.5A Active CN117214461B (en) | 2023-08-31 | 2023-08-31 | Suspension wing type flow velocity meter calibration method based on fusion algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117214461B (en) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1100202A (en) * | 1993-09-06 | 1995-03-15 | 水利部南京水利水文自动化研究所 | Method and apparatus for inspection of propeller-type flowmeter |
JP2000002713A (en) * | 1998-06-18 | 2000-01-07 | Toyota Central Res & Dev Lab Inc | Three-dimensional flow analytical device |
CN1474187A (en) * | 2003-04-01 | 2004-02-11 | 水利部南京水利水文自动化研究所 | Detecting method of current meter and its detection mounting equipment |
CN102207513A (en) * | 2011-03-15 | 2011-10-05 | 水利部南京水利水文自动化研究所 | Method for calibrating rotor current meter |
CN103995147A (en) * | 2014-05-05 | 2014-08-20 | 河海大学 | Data postprocessing system applicable to acoustic Doppler current meter and application thereof |
CN105043462A (en) * | 2015-07-06 | 2015-11-11 | 中国科学院力学研究所 | Method for measuring gas flow of large-diameter pipe |
CN105954543A (en) * | 2016-06-23 | 2016-09-21 | 江苏南水水务科技有限公司 | Current meter high-speed linear formula integration detection method and detection device |
CN108362911A (en) * | 2018-02-09 | 2018-08-03 | 西京学院 | A kind of current meter scaling method based on least square method and interpolation method |
CN208156031U (en) * | 2018-05-07 | 2018-11-27 | 冉乔文 | A kind of measurement device suitable for liquid surface flow velocity |
CN109990842A (en) * | 2019-04-24 | 2019-07-09 | 水利部南京水利水文自动化研究所 | Hydrology current measuring instrument inclination angle performance detector |
CN110658355A (en) * | 2019-11-01 | 2020-01-07 | 烟台龙源电力技术股份有限公司 | Pipeline air flow velocity measuring device |
CN112862898A (en) * | 2021-02-05 | 2021-05-28 | 慧目(重庆)科技有限公司 | Flow velocity measuring method based on computer vision |
US20210223090A1 (en) * | 2018-06-15 | 2021-07-22 | Horiba Stec, Co., Ltd. | Method for calculating piping capacity and calibrator for flow rate control instrument or flow rate measuring instrument |
CN114778888A (en) * | 2022-04-06 | 2022-07-22 | 中国科学院南海海洋研究所 | Inclined flow meter, fluid flow velocity and fluid flow direction calculation method |
-
2023
- 2023-08-31 CN CN202311111714.5A patent/CN117214461B/en active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1100202A (en) * | 1993-09-06 | 1995-03-15 | 水利部南京水利水文自动化研究所 | Method and apparatus for inspection of propeller-type flowmeter |
JP2000002713A (en) * | 1998-06-18 | 2000-01-07 | Toyota Central Res & Dev Lab Inc | Three-dimensional flow analytical device |
CN1474187A (en) * | 2003-04-01 | 2004-02-11 | 水利部南京水利水文自动化研究所 | Detecting method of current meter and its detection mounting equipment |
CN102207513A (en) * | 2011-03-15 | 2011-10-05 | 水利部南京水利水文自动化研究所 | Method for calibrating rotor current meter |
CN103995147A (en) * | 2014-05-05 | 2014-08-20 | 河海大学 | Data postprocessing system applicable to acoustic Doppler current meter and application thereof |
CN105043462A (en) * | 2015-07-06 | 2015-11-11 | 中国科学院力学研究所 | Method for measuring gas flow of large-diameter pipe |
CN105954543A (en) * | 2016-06-23 | 2016-09-21 | 江苏南水水务科技有限公司 | Current meter high-speed linear formula integration detection method and detection device |
CN108362911A (en) * | 2018-02-09 | 2018-08-03 | 西京学院 | A kind of current meter scaling method based on least square method and interpolation method |
CN208156031U (en) * | 2018-05-07 | 2018-11-27 | 冉乔文 | A kind of measurement device suitable for liquid surface flow velocity |
US20210223090A1 (en) * | 2018-06-15 | 2021-07-22 | Horiba Stec, Co., Ltd. | Method for calculating piping capacity and calibrator for flow rate control instrument or flow rate measuring instrument |
CN109990842A (en) * | 2019-04-24 | 2019-07-09 | 水利部南京水利水文自动化研究所 | Hydrology current measuring instrument inclination angle performance detector |
CN110658355A (en) * | 2019-11-01 | 2020-01-07 | 烟台龙源电力技术股份有限公司 | Pipeline air flow velocity measuring device |
CN112862898A (en) * | 2021-02-05 | 2021-05-28 | 慧目(重庆)科技有限公司 | Flow velocity measuring method based on computer vision |
CN114778888A (en) * | 2022-04-06 | 2022-07-22 | 中国科学院南海海洋研究所 | Inclined flow meter, fluid flow velocity and fluid flow direction calculation method |
Non-Patent Citations (1)
Title |
---|
方文强, 叶庆泰: "峰流速仪标定系统数学模型的研究", 机械, no. 03, 30 March 2005 (2005-03-30) * |
Also Published As
Publication number | Publication date |
---|---|
CN117214461B (en) | 2024-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108362911B (en) | Flow velocity meter calibration method based on least square method and interpolation method | |
CN111998919B (en) | Gas meter calibration method and device | |
CN115687983B (en) | Bridge health state monitoring method and system and electronic equipment | |
CN108197339A (en) | The dynamic displacement monitoring method of real time acceleration integral algorithm based on recurrent least square method | |
CN114235111A (en) | Ultrasonic water meter flow calibration method based on model optimization | |
CN113029512A (en) | Air flow field model correction control method for temporary-impulse transonic wind tunnel | |
CN110553631B (en) | Water level measurement series error analysis method about water level flow relation | |
CN109992748B (en) | Estuary tide level harmonic analysis method | |
CN117214461B (en) | Suspension wing type flow velocity meter calibration method based on fusion algorithm | |
CN117332205B (en) | High-precision automatic optimization method and device for temperature compensation resistance of piezoresistor | |
CN112229454A (en) | System and method for measuring flow differential pressure curve of smoke resistance standard part | |
CN112985692A (en) | Atmospheric pressure sensor error calibration method integrating polynomial and learning model | |
CN114877975A (en) | Ultrasonic water meter flow error correction method based on reynolds number of pipeline inlet | |
CN113701711B (en) | High-precision positioning method and system based on Beidou positioning and barometer | |
CN112747773B (en) | Method for improving precision of gyroscope based on Allan variance and random polynomial | |
CN112967825B (en) | Reactivity measurement method based on correction signal uncertainty analysis | |
CN112284366B (en) | Method for correcting course angle error of polarized light compass based on TG-LSTM neural network | |
CN113468693A (en) | Pressure sensor acquisition system based on big data | |
CN114459573A (en) | Novel calibration model and calibration method of ultrasonic gas meter and ultrasonic gas meter | |
CN112528566A (en) | Real-time air quality data calibration method and system based on AdaBoost training model | |
CN106092178A (en) | Improve the data correcting method of certainty of measurement | |
CN111210877A (en) | Method and device for deducing physical property parameters | |
CN116147724B (en) | Metering method suitable for ultrasonic water meter | |
CN104819703B (en) | A kind of angle of inclination parameter fitness method that fulcrum is sensed for inclination angle | |
CN116147741B (en) | NB-IoT based ultrasonic water meter flow calibration method and system |
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 |