CN114235111B - Ultrasonic water meter flow calibration method based on model optimization - Google Patents

Ultrasonic water meter flow calibration method based on model optimization Download PDF

Info

Publication number
CN114235111B
CN114235111B CN202210168812.1A CN202210168812A CN114235111B CN 114235111 B CN114235111 B CN 114235111B CN 202210168812 A CN202210168812 A CN 202210168812A CN 114235111 B CN114235111 B CN 114235111B
Authority
CN
China
Prior art keywords
population
individuals
mathematical model
water meter
test
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
Application number
CN202210168812.1A
Other languages
Chinese (zh)
Other versions
CN114235111A (en
Inventor
杨金合
宋冠锋
毕晨家
沈华刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qingdao Topscomm Communication Co Ltd
Original Assignee
Qingdao Topscomm Communication Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Qingdao Topscomm Communication Co Ltd filed Critical Qingdao Topscomm Communication Co Ltd
Priority to CN202210168812.1A priority Critical patent/CN114235111B/en
Publication of CN114235111A publication Critical patent/CN114235111A/en
Application granted granted Critical
Publication of CN114235111B publication Critical patent/CN114235111B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Measuring Volume Flow (AREA)

Abstract

The invention relates to the technical field of intelligent water affairs and discloses an ultrasonic water meter flow calibration method based on model optimization, which comprises the following steps: s1, controlling the flow rate to obtain a standard group and a test group; s2, performing nonlinear fitting on the flow characteristic curve based on the standard group data, and calculating an initial mathematical model; s3, evaluating the mathematical model based on the test panel data; s4, optimizing model parameters based on the test set data; s5, obtaining the optimal sectional calibration coefficient of a single water meter through a differential evolution algorithm based on a mathematical model; and s6, repeating the steps from s3 to s5, and completing the flow calibration of all the water meters. The method directly uses the original data to solve the calibration coefficient, is simple and efficient, avoids accuracy errors caused by unreasonable manual selection of segmentation points, corrects the calibration coefficient from the source, optimizes an independent model for each water meter, eliminates the influence caused by processing or installation errors, controls the accuracy within a required range, and realizes the high consistency of water meters in the same batch or model.

Description

Ultrasonic water meter flow calibration method based on model optimization
Technical Field
The invention relates to the technical field of intelligent water affairs, in particular to an ultrasonic water meter flow calibration method based on model optimization.
Background
The ultrasonic water meter is widely applied to the fields of civilian use, industry and the like due to the advantages of high metering precision, wide range ratio, small pressure loss and the like, the flow rate is calculated by calculating the propagation time difference of upstream and downstream ultrasonic signals, and the average flow rate in a pipeline is obtained based on a calibration coefficient. At present, in the flow calibration process, the flow segmentation calibration is generally carried out by referring to a theoretical formula of fluid mechanics and according to manual experience, the process is iterated repeatedly, and the optimal solution is difficult to ensure. In addition, time cost and economic cost are considered, when the calibration coefficients of a plurality of water meters of the same type are secondarily corrected, the utilized measurement data are often less, and the accuracy of the water meters in a full-flow interval is difficult to guarantee through a traditional secondary correction method.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides an ultrasonic water meter flow calibration method based on model optimization. The method directly uses the original measurement data to solve the calibration coefficient, and avoids the precision error caused by unreasonable manual selection of the segmentation points. In addition, the calibration coefficient is corrected from the source, independent mathematical model parameters are optimized for each water meter on the basis of the initial mathematical model, the influence caused by machining or installation errors is eliminated, and the precision is controlled within a required range.
The purpose of the invention can be realized by the following technical scheme:
a model optimization-based ultrasonic water meter flow calibration method is characterized by comprising the following steps:
s 1: using flow data measured by one calibrated ultrasonic water meter in a full flow interval as a standard group, and using flow data measured by all ultrasonic water meters to be calibrated in the full flow interval as a test group;
s 2: drawing a flow characteristic curve based on the standard group data, and carrying out nonlinear fitting on the flow characteristic curve to obtain an initial mathematical model;
s 3: calculating relative error based on test group dataRDAccording to relative errorRDEvaluating the mathematical model to determine whether to skip step s 4;
s 4: optimizing mathematical model parameters through a differential evolution algorithm based on the test group data;
s 5: based on a mathematical model, obtaining an optimal sectional calibration coefficient of a single ultrasonic water meter to be calibrated through a differential evolution algorithm;
s 6: and (5) repeating the steps s3 to s5 to finish the flow calibration of all the ultrasonic water meters to be calibrated.
Further, the data set in step s1 includes the average flow velocity measured by each ultrasonic water meter and the average value of the propagation time difference between the upstream and downstream received signals
Figure 875268DEST_PATH_IMAGE001
Wherein the average flow velocity
Figure 147111DEST_PATH_IMAGE002
The calculation formula of (c) is as follows:
Figure 334510DEST_PATH_IMAGE003
in the formula (I), the compound is shown in the specification,T j is as followsjThe duration of the secondary measurement is long,m j is as followsjDuration of time of secondary measurementT j The mass of the water flowing into the tank,
Figure 506735DEST_PATH_IMAGE004
is as followsjThe density of the water at the time of the secondary measurement,Sthe sectional area of the pipeline is measured.
Further, the upstream and downstream received signal propagation time difference is calculated by a time difference method with the upstream and downstream received signals as input.
Further, the flow characteristic curve in the step s2 is a curve of the calibration coefficient varying with the average value of the propagation time difference of the upstream and downstream received signals, wherein the calibration coefficientKThe calculation formula is as follows:
Figure 754176DEST_PATH_IMAGE005
Figure 513316DEST_PATH_IMAGE006
in the formula (I), the compound is shown in the specification,
Figure 301143DEST_PATH_IMAGE007
is at the firstjAt the time of secondary measurementiUpstream and downstream of the group the propagation time difference of the received signals, NjIs as followsjThe total number of differences in propagation time between the upstream and downstream received signals at the time of the secondary measurement,
Figure 327874DEST_PATH_IMAGE008
the average flow rate for the standard set obtained for the calibrated ultrasonic water meter for the jth measurement,
Figure 746217DEST_PATH_IMAGE009
is as followsjThe average of the upstream and downstream received signal propagation time differences of the secondary measurements.
Further, the non-linear fitting in the step s2 is a least square fitting, and an initial mathematical model is obtained
Figure 727074DEST_PATH_IMAGE010
Wherein
Figure 521854DEST_PATH_IMAGE011
Are model parameters.
Further, the relative error in the step s3RDThe calculation formula of (a) is as follows:
Figure 403092DEST_PATH_IMAGE012
in the formula (I), the compound is shown in the specification,
Figure 992336DEST_PATH_IMAGE013
for the test grouplThe average flow rate of the group is,Kfor the calibration coefficients calculated by the mathematical model,
Figure 460489DEST_PATH_IMAGE014
is the first in the test grouplAn average of the upstream and downstream received signal propagation time differences of the group;
the mathematical model evaluation principle is as follows:
if it isRD l If the error values are less than or equal to one half of the maximum allowable error, the mathematical model is considered to be consistent with the test set data, and the step s4 is skipped;
if present, isRD l If the maximum allowable error is greater than one half of the maximum allowable error and less than or equal to the maximum allowable error, the mathematical model is considered not to be consistent with the test group data, and the step s4 is required to optimize the parameters of the mathematical model;
if present, isRD l If the error is larger than the maximum allowable error, the mathematical model is considered not to meet the requirement, and the ultrasonic water meter is analyzed independently by departing from the process.
Further, the step s4 includes the following specific steps:
s 4.1: average flow rate values based on test group data
Figure 714752DEST_PATH_IMAGE013
And (a)
Figure 732387DEST_PATH_IMAGE015
) Define an objective function by the Euclidean distance between the two, and use a mathematical model
Figure 243265DEST_PATH_IMAGE016
Middle parameter
Figure 447981DEST_PATH_IMAGE017
The preset allowable range of the model is a constraint condition, and variables are designed based on a mathematical model parameter theta;
s 4.2: setting a search range of a differential evolution algorithm on a mathematical model, and initializing the number of population individuals, the maximum iteration times, cross factors and variation factors;
s 4.3: randomly generating a population as a parent population in the search range of the mathematical model parameter theta, adding 1 to the iteration times, and calculating the fitness of individuals in the parent population according to the target function and the constraint condition;
s 4.4: randomly selecting individuals from the parent population to perform crossing and mutation based on the crossing factors and the mutation factors to generate a test population;
s 4.5: calculating the fitness of individuals in the test population, comparing the fitness of individuals in the parent population with the fitness of individuals in the test population, and generating a child population according to a comparison result;
s 4.6: performing non-dominant sorting of individuals in the offspring population;
s 4.7: calculating crowding distances of individuals in the offspring population, and eliminating individuals with relatively small crowding distances to keep the number of the individuals in the offspring population consistent with the number of the individuals in the initial population;
s 4.8: judging whether the existing iteration times reach the requirement of the maximum iteration times, if so, finishing the parameter optimization of the mathematical model, and ending the process; otherwise jump to step s 4.3.
Further, the step s5 specifically includes the following steps:
s 5.1: average flow rate values based on test group data
Figure 833832DEST_PATH_IMAGE013
And (a)
Figure 174815DEST_PATH_IMAGE015
) Defining an objective function by the Euclidean distance between the two, randomly generating a segmentation point position beta and a segmentation calibration coefficient alpha, taking the allowable range of the segmentation point position beta and the segmentation calibration coefficient alpha as a constraint condition, and taking parameters
Figure 387753DEST_PATH_IMAGE018
Designing variables;
s 5.2: setting a search range of a differential evolution algorithm for the parameter gamma, and initializing the number of population individuals, the maximum iteration times, cross factors and variation factors;
s 5.3: randomly generating a population as a parent population in the search range of the parameter gamma, adding 1 to the iteration number, and calculating the fitness of individuals in the parent population according to the target function and the constraint condition;
s 5.4: randomly selecting individuals from the parent population to perform crossing and mutation based on the crossing factors and the mutation factors to generate a test population;
s 5.5: calculating the fitness of individuals in the test population, comparing the fitness of individuals in the parent population with the fitness of individuals in the test population, and generating a child population according to a comparison result;
s 5.6: performing non-dominated sorting of individuals in the progeny population;
s 5.7: calculating crowding distances of individuals in the offspring population, and eliminating individuals with relatively small crowding distances to keep the number of the individuals in the offspring population consistent with the number of the individuals in the initial population;
s 5.8: judging whether the existing iteration times reach the requirement of the maximum iteration times, if so, completing parameter optimization, and ending the process; otherwise jump to step s 5.3.
The invention has the beneficial technical effects that: the method has the advantages that the calibration coefficient is solved by directly using the original measurement data, simplicity and high efficiency are realized, and the precision error caused by unreasonable manual selection of the segmentation points is avoided. In addition, the calibration coefficient is corrected from the source, independent mathematical model parameters are optimized for each water meter on the basis of the initial mathematical model, the influence caused by machining or installation errors is eliminated, and the precision is controlled within a required range, so that the high consistency of water meters in the same batch or model is realized.
Drawings
FIG. 1 is a general flow diagram of the present invention.
Fig. 2 is a flow characteristic curve of the ultrasonic water meter according to the embodiment of the invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and do not limit the invention.
Example (b):
as shown in fig. 1, a method for calibrating flow of an ultrasonic water meter based on model optimization includes the following steps:
s 1: and adjusting the flow control device, and after the flow is stable, using the flow data measured by one calibrated ultrasonic water meter in the full-flow interval as a standard group, and using the flow data measured by all the calibrated ultrasonic water meters in the full-flow interval as a test group. The data set comprises the average flow speed measured by each ultrasonic water meter and the average value of the propagation time difference of the upstream and downstream received signals
Figure 79765DEST_PATH_IMAGE019
Wherein the average flow velocity
Figure 738148DEST_PATH_IMAGE020
The calculation formula of (a) is as follows:
Figure 622053DEST_PATH_IMAGE003
in the formula (I), the compound is shown in the specification,T j is a firstjThe duration of the secondary measurement is long,m j is a firstjDuration of time of secondary measurementT j The mass of water flowing into the tank,
Figure 989581DEST_PATH_IMAGE004
is a firstjThe density of the water at the time of the secondary measurement,Sthe sectional area of the pipeline at the section is measured.
The propagation time difference of the upstream and downstream received signals is calculated by a time difference method by taking the upstream and downstream received signals as input.
s 2: and drawing a flow characteristic curve based on the standard group data, and carrying out nonlinear fitting on the flow characteristic curve to obtain an initial mathematical model. The flow characteristic curve is a variation curve of a calibration coefficient along with the average value of the propagation time difference of upstream and downstream received signals, wherein the calibration coefficientKThe calculation formula is as follows:
Figure 886998DEST_PATH_IMAGE005
Figure 365384DEST_PATH_IMAGE006
in the formula (I), the compound is shown in the specification,
Figure 708989DEST_PATH_IMAGE021
is at the firstjAt the time of secondary measurementiUpstream and downstream of the group, the difference in propagation time of the received signals, NjIs a firstjThe total number of differences in propagation time between the upstream and downstream received signals at the time of the secondary measurement,
Figure 981839DEST_PATH_IMAGE008
the average flow rate for the standard set obtained for the calibrated ultrasonic water meter for the jth measurement,
Figure 897711DEST_PATH_IMAGE009
is as followsjThe average of the upstream and downstream received signal propagation time differences of the secondary measurements.
The nonlinear fitting is least square fitting to obtain an initial mathematical model
Figure 914209DEST_PATH_IMAGE022
. The mathematical model expression is related to the ultrasonic water meter model, and the parameters in the mathematical model are related to the water meter pipeline structure and the sound path design parameters.
As shown in fig. 2, which is a flow characteristic curve of a certain model DN15 ultrasonic water meter, the calibration coefficient changes with the change of the average value of the propagation time difference of the upstream and downstream received signals. After least squares fitting based on the data of fig. 2, the expression of the initial mathematical model is obtained as follows:
Figure 100602DEST_PATH_IMAGE023
in the formula, a, b, c and d are all adjustable parameters in the mathematical model, and the values are related to the model of the ultrasonic water meter. The value range of the parameter a is [0.02016, 0.02107], the value range of the parameter b is [ -0.0001136, 0.0003671], the value range of the parameter c is [ -0.009227, -0.008075], and the value range of the parameter d is [ -0.2465, -0.1726 ].
s 3: calculating relative error based on test group dataRDAccording to the relative errorRDThe mathematical model is evaluated to determine whether to skip step s4 as a result of the evaluation. Relative errorRDThe calculation formula of (a) is as follows:
Figure 996882DEST_PATH_IMAGE012
in the formula (I), the compound is shown in the specification,
Figure 150783DEST_PATH_IMAGE013
to test grouplThe average flow rate of the group is,Kfor the calibration coefficients calculated by the mathematical model,
Figure 298868DEST_PATH_IMAGE014
is the first in the test grouplThe average of the differences in propagation time of the received signals upstream and downstream of the group.
The mathematical model evaluation principle is as follows:
if it isRD l If the error values are less than or equal to one half of the maximum allowable error, the mathematical model is considered to be consistent with the test group data, and the step s4 is skipped;
if present, isRD l If the maximum allowable error is more than one half of the maximum allowable error and less than or equal to the maximum allowable error, the mathematical model is considered not to be in accordance with the test group data, and the step s4 needs to be carried out to optimize the parameters of the mathematical model;
if present, isRD l If the error is larger than the maximum allowable error, the mathematical model is considered not to meet the requirement, and the ultrasonic water meter is analyzed independently without the process.
In the embodiment, a certain model DN15 ultrasonic water meter is a 2-level meter, and the set maximum allowable error is 2-level precision required by national standards.
s 4: and optimizing the mathematical model parameters through a differential evolution algorithm based on the test group data. The method comprises the following specific steps:
s 4.1: average flow rate value based on test group data
Figure 11871DEST_PATH_IMAGE013
And (a) and
Figure 626523DEST_PATH_IMAGE015
) Define an objective function by the Euclidean distance between the two, and use a mathematical model
Figure 487294DEST_PATH_IMAGE016
Middle parameter
Figure 845594DEST_PATH_IMAGE017
The preset allowable range is a constraint condition, and variables are designed based on a mathematical model parameter theta;
s 4.2: setting a search range of a differential evolution algorithm on a mathematical model, and initializing the number of population individuals, the maximum iteration times, cross factors and variation factors;
s 4.3: randomly generating a population as a parent population in the search range of the mathematical model parameter theta, adding 1 to the iteration number, and calculating the fitness of individuals in the parent population according to the target function and the constraint condition;
s 4.4: randomly selecting individuals from the parent population for crossing and mutation based on the crossing factors and the mutation factors to generate a test population;
s 4.5: calculating the fitness of individuals in the test population, comparing the fitness of individuals in the parent population with the fitness of individuals in the test population, and generating a child population according to a comparison result;
s 4.6: performing non-dominant sorting of individuals in the offspring population;
s 4.7: calculating crowding distances of individuals in the offspring population, and eliminating individuals with relatively small crowding distances to keep the number of the individuals in the offspring population consistent with the number of the individuals in the initial population;
s 4.8: judging whether the existing iteration times reach the requirement of the maximum iteration times, if so, finishing the parameter optimization of the mathematical model, and ending the process; otherwise jump to step s 4.3.
s 5: and based on the mathematical model, obtaining the optimal sectional calibration coefficient of the single ultrasonic water meter to be calibrated through a differential evolution algorithm. Sectional schoolQuasi-coefficients are logarithmic models
Figure 505114DEST_PATH_IMAGE024
The specific steps of s5 are as follows:
s 5.1: average flow rate value based on test group data
Figure 385608DEST_PATH_IMAGE013
And (a) and
Figure 497789DEST_PATH_IMAGE015
) The Euclidean distance between the two points defines an objective function, a segment point position beta and a segment calibration coefficient alpha are randomly generated, the allowable range of the segment point position beta and the segment calibration coefficient alpha is used as a constraint condition, and parameters are used
Figure 659781DEST_PATH_IMAGE018
Designing variables;
s 5.2: setting a search range of a differential evolution algorithm for the parameter gamma, and initializing the number of population individuals, the maximum iteration times, cross factors and variation factors;
s 5.3: randomly generating a population as a parent population in the search range of the parameter gamma, adding 1 to the iteration number, and calculating the fitness of individuals in the parent population according to the target function and the constraint condition;
s 5.4: randomly selecting individuals from the parent population for crossing and mutation based on the crossing factors and the mutation factors to generate a test population;
s 5.5: calculating the fitness of individuals in the test population, comparing the fitness of individuals in the parent population with the fitness of individuals in the test population, and generating a child population according to a comparison result;
s 5.6: performing non-dominant sorting of individuals in the offspring population;
s 5.7: calculating crowding distances of individuals in the offspring population, and eliminating individuals with relatively small crowding distances to keep the number of the individuals in the offspring population consistent with the number of the individuals in the initial population;
s 5.8: judging whether the existing iteration times reach the requirement of the maximum iteration times, if so, completing parameter optimization, and ending the process; otherwise jump to step s 5.3.
s 6: and (5) repeating the steps s3 to s5 to finish the flow calibration of all the ultrasonic water meters to be calibrated.
The above-mentioned embodiments are illustrative of the specific embodiments of the present invention, and not restrictive, and those skilled in the relevant art can make various changes and modifications to the invention without departing from the spirit and scope of the invention, so that all equivalent technical solutions should fall within the scope of the present invention.

Claims (3)

1. A model optimization-based ultrasonic water meter flow calibration method is characterized by comprising the following steps:
s 1: the method comprises the steps that flow data measured by one calibrated ultrasonic water meter in a full-flow interval are used as a standard group, and flow data measured by all to-be-calibrated ultrasonic water meters in the full-flow interval are used as a test group; the finally obtained data set comprises the average flow speed measured by each ultrasonic water meter and the average value of the propagation time difference of the upstream and downstream received signals
Figure 432597DEST_PATH_IMAGE001
Wherein the average flow velocity
Figure 511411DEST_PATH_IMAGE002
The calculation formula of (a) is as follows:
Figure 947071DEST_PATH_IMAGE003
in the formula (I), the compound is shown in the specification,T j is a firstjThe duration of the secondary measurement is long,m j is a firstjDuration of time of secondary measurementT j The mass of water flowing into the tank,
Figure 892287DEST_PATH_IMAGE004
is a firstjSecond measurementThe density of water at the time of measurement,Smeasuring the sectional area of the pipeline at the section;
s 2: drawing a flow characteristic curve based on the standard group data, and carrying out nonlinear fitting on the flow characteristic curve to obtain an initial mathematical model; wherein the flow characteristic curve is a curve in which a calibration coefficient varies with an average value of propagation time differences of upstream and downstream received signals, wherein the calibration coefficientKThe calculation formula is as follows:
Figure 46188DEST_PATH_IMAGE005
Figure 866376DEST_PATH_IMAGE006
in the formula (I), the compound is shown in the specification,
Figure 156543DEST_PATH_IMAGE007
is at the firstjAt the time of secondary measurementiUpstream and downstream of the group, the difference in propagation time of the received signals, NjIs as followsjThe total number of differences in propagation time between the upstream and downstream received signals at the time of the secondary measurement,
Figure 472993DEST_PATH_IMAGE008
the average flow rate for the standard set obtained for the calibrated ultrasonic water meter for the jth measurement,
Figure 848611DEST_PATH_IMAGE009
is as followsjAn average value of the upstream and downstream received signal propagation time differences of the secondary measurements;
s 3: calculating relative error based on test set dataRDAccording to relative errorRDEvaluating the mathematical model to determine whether to skip step s 4; wherein the relative errorRDThe calculation formula of (c) is as follows:
Figure 472490DEST_PATH_IMAGE010
in the formula (I), the compound is shown in the specification,
Figure 882743DEST_PATH_IMAGE011
to test grouplThe average flow rate of the group is,Kfor the calibration coefficients calculated by the mathematical model,
Figure 169761DEST_PATH_IMAGE012
is the first in the test grouplAn average of the upstream and downstream received signal propagation time differences of the group;
the mathematical model evaluation principle is as follows:
if it isRD l If the error values are less than or equal to one half of the maximum allowable error, the mathematical model is considered to be consistent with the test group data, and the step s4 is skipped;
if present, isRD l If the maximum allowable error is greater than one half of the maximum allowable error and less than or equal to the maximum allowable error, the mathematical model is considered not to be consistent with the test group data, and the step s4 is required to optimize the parameters of the mathematical model;
if present, isRD l If the error is larger than the maximum allowable error, the mathematical model is considered not to meet the requirement, the process needs to be separated, and the ultrasonic water meter is analyzed independently;
s 4: optimizing mathematical model parameters by a differential evolution algorithm based on test group data; the method comprises the following specific steps:
s 4.1: average flow rate value based on test group data
Figure 298254DEST_PATH_IMAGE013
And (a) and
Figure 460245DEST_PATH_IMAGE014
) Define an objective function by the Euclidean distance between the two, and use a mathematical model
Figure 725004DEST_PATH_IMAGE015
Middle parameter
Figure 274934DEST_PATH_IMAGE016
The preset allowable range is a constraint condition, and variables are designed based on a mathematical model parameter theta;
s 4.2: setting a search range of a differential evolution algorithm on a mathematical model, and initializing the number of population individuals, the maximum iteration times, cross factors and variation factors;
s 4.3: randomly generating a population as a parent population in the search range of the mathematical model parameter theta, adding 1 to the iteration number, and calculating the fitness of individuals in the parent population according to the target function and the constraint condition;
s 4.4: randomly selecting individuals from the parent population to perform crossing and mutation based on the crossing factors and the mutation factors to generate a test population;
s 4.5: calculating the fitness of individuals in the test population, comparing the fitness of individuals in the parent population with the fitness of individuals in the test population, and generating a child population according to a comparison result;
s 4.6: performing non-dominated sorting of individuals in the progeny population;
s 4.7: calculating crowding distances of individuals in the offspring population, and eliminating individuals with relatively small crowding distances to keep the number of the individuals in the offspring population consistent with the number of the individuals in the initial population;
s 4.8: judging whether the existing iteration times reach the requirement of the maximum iteration times, if so, finishing the parameter optimization of the mathematical model, and ending the process; otherwise, jumping to the step s 4.3;
s 5: based on a mathematical model, obtaining an optimal sectional calibration coefficient of a single ultrasonic water meter to be calibrated through a differential evolution algorithm; the method comprises the following specific steps:
s 5.1: average flow rate value based on test group data
Figure 123679DEST_PATH_IMAGE013
And (a)
Figure 823782DEST_PATH_IMAGE014
) Defining an objective function by the Euclidean distance between the two points, randomly generating a segmentation point position beta and a segmentation calibration coefficient alpha, and using the segmentation point position beta and the segmentation calibration coefficientThe allowable range of alpha is used as a constraint condition, and the parameter is used as
Figure 208627DEST_PATH_IMAGE017
Designing variables;
s 5.2: setting a search range of a differential evolution algorithm for the parameter gamma, and initializing the number of population individuals, the maximum iteration times, cross factors and variation factors;
s 5.3: randomly generating a population as a parent population in the search range of the parameter gamma, adding 1 to the iteration times, and calculating the fitness of individuals in the parent population according to the target function and the constraint condition;
s 5.4: randomly selecting individuals from the parent population to perform crossing and mutation based on the crossing factors and the mutation factors to generate a test population;
s 5.5: calculating the fitness of individuals in the test population, comparing the fitness of individuals in the parent population with the fitness of individuals in the test population, and generating a child population according to a comparison result;
s 5.6: performing non-dominant sorting of individuals in the offspring population;
s 5.7: calculating crowding distances of individuals in the offspring population, and eliminating individuals with relatively small crowding distances to keep the number of the individuals in the offspring population consistent with the number of the individuals in the initial population;
s 5.8: judging whether the existing iteration times reach the maximum iteration time requirement, if so, completing parameter optimization, and ending the process; otherwise, jumping to the step s 5.3;
s 6: and (5) repeating the steps s3 to s5 to finish the flow calibration of all the ultrasonic water meters to be calibrated.
2. The method of claim 1, wherein the difference in propagation time between the upstream and downstream received signals is calculated by a time difference method using the upstream and downstream received signals as inputs.
3. The method for calibrating flow rate of an ultrasonic water meter based on model optimization of claim 1, wherein the step s is performed in a manner of2, obtaining an initial mathematical model by using least square fitting as the nonlinear fitting
Figure DEST_PATH_IMAGE019
Wherein
Figure DEST_PATH_IMAGE021
Are model parameters.
CN202210168812.1A 2022-02-24 2022-02-24 Ultrasonic water meter flow calibration method based on model optimization Active CN114235111B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210168812.1A CN114235111B (en) 2022-02-24 2022-02-24 Ultrasonic water meter flow calibration method based on model optimization

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210168812.1A CN114235111B (en) 2022-02-24 2022-02-24 Ultrasonic water meter flow calibration method based on model optimization

Publications (2)

Publication Number Publication Date
CN114235111A CN114235111A (en) 2022-03-25
CN114235111B true CN114235111B (en) 2022-07-15

Family

ID=80747867

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210168812.1A Active CN114235111B (en) 2022-02-24 2022-02-24 Ultrasonic water meter flow calibration method based on model optimization

Country Status (1)

Country Link
CN (1) CN114235111B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114485865B (en) * 2022-04-15 2022-09-27 青岛鼎信通讯股份有限公司 Ultrasonic water meter flow calibration method based on Shannon entropy
CN115950495B (en) * 2023-02-09 2023-06-09 青岛鼎信通讯科技有限公司 Ultrasonic water meter measurement frequency adjusting method
CN116793462B (en) * 2023-08-21 2023-11-03 成都千嘉科技股份有限公司 Calibration method and calibration device based on ultrasonic gas meter
CN117268512B (en) * 2023-11-23 2024-02-09 青岛鼎信通讯科技有限公司 Consistency optimization method suitable for ultrasonic water meter
CN117367527B (en) * 2023-12-08 2024-03-19 青岛鼎信通讯科技有限公司 Metering method capable of improving reliability of ultrasonic water meter

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101048645A (en) * 2004-08-13 2007-10-03 恩特格里公司 System and method for calibration of a flow device
CN102435265A (en) * 2011-09-19 2012-05-02 苏州聚元微电子有限公司 Precision correction method for transit time ultrasonic flowmeter
CN105181040A (en) * 2015-09-25 2015-12-23 辽宁聚焦科技有限公司 Method for digitized calibration and optimization of difference pressure flow meter
CN107270983A (en) * 2017-07-05 2017-10-20 青岛积成电子股份有限公司 The method that data based on instantaneous delivery build ultrasonic flow model
CN109724657A (en) * 2018-12-29 2019-05-07 杭州莱宸科技有限公司 Watermeter flowing rate metering method and system based on modified Delphi approach
CN113033093A (en) * 2021-03-24 2021-06-25 南京航空航天大学 Simulation model-based system design parameter multi-objective optimization method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101048645A (en) * 2004-08-13 2007-10-03 恩特格里公司 System and method for calibration of a flow device
CN102435265A (en) * 2011-09-19 2012-05-02 苏州聚元微电子有限公司 Precision correction method for transit time ultrasonic flowmeter
CN105181040A (en) * 2015-09-25 2015-12-23 辽宁聚焦科技有限公司 Method for digitized calibration and optimization of difference pressure flow meter
CN107270983A (en) * 2017-07-05 2017-10-20 青岛积成电子股份有限公司 The method that data based on instantaneous delivery build ultrasonic flow model
CN109724657A (en) * 2018-12-29 2019-05-07 杭州莱宸科技有限公司 Watermeter flowing rate metering method and system based on modified Delphi approach
CN113033093A (en) * 2021-03-24 2021-06-25 南京航空航天大学 Simulation model-based system design parameter multi-objective optimization method

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
A Novel Mathematical Model for Transit-time Ultrasonic Flow Measurement;Lei Kang;《2019 IEEE International Ultrasonics Symposium (IUS)》;20191209;全文 *
Neural-network-based error reduction in calibrating utility ultrasonic flow meters;B. Yazdanshenasshad;《Flow Measurement and Instrumentation》;20181011;全文 *
Ultrasonic flow detection and research on regression analysis calibration method;Ming-liang Liu;《2012 International Conference on Systems and Informatics (ICSAI2012)》;20120625;全文 *
基于回波峰值拟合的气体超声流量计信号处理;田雷;《电子测量与仪器学报》;20170731;第31卷(第7期);全文 *

Also Published As

Publication number Publication date
CN114235111A (en) 2022-03-25

Similar Documents

Publication Publication Date Title
CN114235111B (en) Ultrasonic water meter flow calibration method based on model optimization
CN109001136B (en) COD on-line monitoring method based on ultraviolet-visible light absorption spectrum
CN111291937A (en) Method for predicting quality of treated sewage based on combination of support vector classification and GRU neural network
CN112926265A (en) Atmospheric porous probe measurement calibration method based on genetic algorithm optimization neural network
CN107016622B (en) Inversion method for urban water supply network node water demand containing large user water consumption information
CN101473196A (en) Measuring apparatus and methods of using them
CN110346005B (en) Coriolis mass flowmeter digital signal processing method based on deep learning
CN103177289B (en) Modeling method for noise-uncertainty complicated nonlinear dynamic system
CN108362911A (en) A kind of current meter scaling method based on least square method and interpolation method
CN114964429A (en) Ultrasonic water meter calibration method
CN111914487B (en) Data-free regional hydrological parameter calibration method based on antagonistic neural network
CN112925984A (en) GCN recommendation-based sample density aggregation method
CN112818495A (en) Novel dynamic correction method for pipeline pressure drop measurement and calculation algorithm parameters
CN114485865B (en) Ultrasonic water meter flow calibration method based on Shannon entropy
CN114487976B (en) Method and system for evaluating traceability uncertainty of MCM electronic transformer calibrator
CN111351862A (en) Ultrasonic measurement calibration method and thickness measurement method
CN112685954B (en) Method and device for predicting wind speed and fan rotating speed of automobile environment wind tunnel
CN114924489A (en) Model autonomous learning method suitable for process industry prediction control
CN116147724B (en) Metering method suitable for ultrasonic water meter
CN115018137A (en) Water environment model parameter calibration method based on reinforcement learning
CN110321650B (en) Structural reliability analysis method based on novel test design and weight response surface
CN110705186A (en) Real-time online instrument checking and diagnosing method through RBF particle swarm optimization algorithm
CN114485863A (en) Flow error correction method, system, computer and medium for ultrasonic water meter
CN116304968B (en) Ultrasonic water meter flow data fusion method and device based on BP neural network
CN116592976A (en) Method for transmitting redundant measurement value of single-measuring-point water level gauge of lock chamber of ship lock

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