CN102486390B - Method for steam pipe network of iron and steel enterprise to correct metering data - Google Patents

Method for steam pipe network of iron and steel enterprise to correct metering data Download PDF

Info

Publication number
CN102486390B
CN102486390B CN201010572857.2A CN201010572857A CN102486390B CN 102486390 B CN102486390 B CN 102486390B CN 201010572857 A CN201010572857 A CN 201010572857A CN 102486390 B CN102486390 B CN 102486390B
Authority
CN
China
Prior art keywords
data
value
flow
flowmeter
steam pipe
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.)
Expired - Fee Related
Application number
CN201010572857.2A
Other languages
Chinese (zh)
Other versions
CN102486390A (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.)
Shenyang Institute of Automation of CAS
Original Assignee
Shenyang Institute of Automation of CAS
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 Shenyang Institute of Automation of CAS filed Critical Shenyang Institute of Automation of CAS
Priority to CN201010572857.2A priority Critical patent/CN102486390B/en
Publication of CN102486390A publication Critical patent/CN102486390A/en
Application granted granted Critical
Publication of CN102486390B publication Critical patent/CN102486390B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention relates to a method for a steam pipe network of an iron and steel enterprise to correct metering data, which comprises the following steps of: judging whether the working state of a flow meter is normal or not; if YES, judging whether flow data received by the flow meter has abnormal data; if YES, removing the abnormal data; accumulating the flow data after removing the abnormal data into a certification cycle, and carrying out certification to obtain certification data; and comparing the accumulated data and the certification data, if the accumulated data and the certification data have obvious differences, correcting the errors of the flow meter; accurately correcting the corrected flow data according to production consumption balance by a method of data coordination to obtain reasonable metering data of the steam pipe network. According to the method disclosed by the invention, when the flow meter has faults, the data of a point to be metered can be estimated, abnormal values in the measurement values can be timely found and removed, and the accuracy, the reliability and the rationality of the values of the steam flow meter can be improved. The method has an important significance for the iron and steel enterprise to improve the steam management level and reduce the energy consumption.

Description

A kind of method that is applicable to iron and steel enterprise's steam pipe system continuous data correction
Technical field
the present invention relates to a kind ofcontinuous data bearing calibration, specifically a kind of method that is applicable to iron and steel enterprise's steam pipe system continuous data correction.
Background technology
In iron and steel enterprise, in order to realize the comprehensive utilization of residual heat and energy, generally there is complicated steam pipe system that multiple steaming a little with vapour point connected together.Monitor in real time and the generation of accurate-metering steam and the basis of the steam scheduling of use amount ，Shi iron and steel enterprise and energy management.
Traditional pipe network measuring system, relies on the variable of flowmeter completely to the flow measurement of pipeline each point.But due to flowmeter calibration difficulties, with lengthening precision service time, reduce, and when there is the situations such as Causing Flowmeter Damages, converting means fault, transmission channel break down, the continuous data of steam pipe system differs and correctly reflects surely the real traffic of steam in pipe network, thereby has a strong impact on production scheduling and managing power consumption and examination.
The scheduling of iron and steel enterprise production of steam and management, always wish the actual value based on steam pipe system, what can access is variable (measured value), and between continuous data and actual value, always have error.Data correction is to adopt certain algorithm to eliminate the appreciable error in measured value to measurement data, and reduces stochastic error, and the data that make to obtain more approach a kind of method of actual value.The method is proposed in 1961 by Kuehn and Davidson the earliest, and the correlative study of Chinese scholars afterwards, makes progressively systematization of the method, but apply in actual production, also has all difficulties.
The difficult point of iron and steel enterprise's steam measurement Data correction shows the following aspects.
1. the characteristic of steam is complicated, is difficult to accurate-metering
Steam is different in different pressure and temperature lower densities, and the thermal loss of steam in transmission, causes producing two-phase fluid, and the amount that produces condensate water is also difficult to estimate.With traditional flowmeter, be difficult to steam to implement temperature and pressure compensation, and estimate the loss of pipe network.Therefore be there is to relatively large deviation in steam flow measurement.
2. the precision of flowmeter itself is not high, and with lengthening service time, error can strengthen
The main orifice plate of flowmeter, V cone are main, and it is obvious that the precision of this type flow meter is subject to the impact of environment, tenure of use, and measuring error can progressively expand.Once and very difficult calibration is installed at the flowmeter of iron and steel enterprise's steam pipe system.
3. steam pipe system complex structure, the situation of have stoichiometric point not, diffusing and leaking,
The steam pipe system complex structure of iron and steel enterprise, outer contact is many.And steam production and the steam consumption of major equipment are often only concerned about in the production run of enterprise, for little, with vapour point, generally all there is no metering.Therefore the metering redundance of iron and steel enterprise's steam pipe system is very little, to Data correction, brings difficulty.And diffusing and leakage problem of steam pipe system also aggravates the gap between the metering of Liao Qi source and user's variable.Can not simply according to material balance principle, carry out Data correction, but according to practical experience, on-the-spot data rationally be processed.
4. stop the generation of the situation of measuring, aggravation cumulative flow value departs from its actual value.
When flowmeter fault or communication disruption, in the long-range measurement centre that does cummulative metering, just there will be the situation that stops metering, make the generating capacity of steam and use amount depart from actual value farther.
Summary of the invention
For above shortcomings part in prior art, the technical problem to be solved in the present invention is to provide a kind of precision of steam-flow meter value, method that is applicable to iron and steel enterprise's steam pipe system continuous data correction of reliability of improving.
For solving the problems of the technologies described above, the technical solution used in the present invention is:
The method that the present invention is applicable to iron and steel enterprise's steam pipe system continuous data correction comprises the following steps:
Judge that whether flowmeter duty is normal;
If flowmeter duty is normal, judge in the data on flows that flowmeter receives whether have abnormal data;
If there is abnormal data, remove;
The data on flows of removing after abnormal data is accumulate to an authentication period, authenticates, obtain verify data;
The data of accumulation and verify data are contrasted, if the two exists significant difference, flowmeter error is proofreaied and correct;
The method that data on flows after proofreading and correct is coordinated by data according to product consumption balance is carried out fine correction, obtains rational steam pipe system continuous data.
The data of accumulation and verify data are contrasted, if the two does not exist significant difference, the method for directly data on flows being coordinated by data according to product consumption balance is carried out fine correction, obtains rational steam pipe system continuous data.
If there is no abnormal data, continues data on flows is accumulate to an authentication period, authenticates, and obtains verify data step.
If flowmeter working state abnormal, estimates flowmeter variable, obtain the estimated value of data on flows, continue and judge in the data on flows that flowmeter receives whether have abnormal data step.
It is described that flowmeter survey value is carried out to estimation procedure is as follows:
According to pipe network structure, determine and flowmeter in normal metering state nearest with this flowmeter measuring point, be designated as REF;
If:
Figure 154668DEST_PATH_IMAGE001
In formula,
Figure 552151DEST_PATH_IMAGE002
for the flow actual value of this flowmeter position,
Figure 281073DEST_PATH_IMAGE003
the flow actual value of WeiREF position, kfor correlation coefficient.Choosing available N history samples value of up-to-date preservation calculates k:
Figure 422204DEST_PATH_IMAGE004
Wherein i=1 ~ nfor integer, represent sampled value sequence number; nlength for sample sequence.
Figure 982499DEST_PATH_IMAGE005
for this flowmeter i flow measurements in N history samples value,
Figure 488566DEST_PATH_IMAGE006
for REF exists ni measured value of individual history samples value;
The estimated value of this measurement point present flow rate is:
Figure 767101DEST_PATH_IMAGE007
In formula
Figure 714853DEST_PATH_IMAGE008
measured value for current sampling instant REF.
Describedly judge flowmeter duty whether normal processes is as follows;
The absolute undulate quantity of definition data on flows is:
Figure 67337DEST_PATH_IMAGE009
Wherein
Figure 72202DEST_PATH_IMAGE010
(r=1 ~ m) for from current time, the data on flows sequence that is M to length by some cycles continuous acquisition; R is integer, represents the sequence number of sampled value, mlength for sample sequence;
During flowmeter normal operation vthe minimum of variation range and maximum statistical value are with ;
When time, think this flowmeter communication disruption;
When
Figure 375827DEST_PATH_IMAGE014
time, think this flowmeter fault;
Judge whether that the process that has abnormal data and remove is as follows:
If the measured value of t sampling instant is
Figure 363375DEST_PATH_IMAGE015
, normal
Figure 650000DEST_PATH_IMAGE015
should meet the following conditions:
i)
Figure 977076DEST_PATH_IMAGE016
;
ii)
Figure 323743DEST_PATH_IMAGE017
；
In formula,
Figure 736270DEST_PATH_IMAGE018
, represent respectively the minimum and maximum value that this measurement point steam flow can reach in technological design;
Figure 801977DEST_PATH_IMAGE020
represent the sampling period;
Figure 257229DEST_PATH_IMAGE021
be illustrated in the maximum possible rate of change of steam flow in the sampling period;
As measured value is
Figure 219369DEST_PATH_IMAGE015
i satisfies condition when different), ii), think exceptional value, this measuring point data is removed, and by the alternative current measurement data of measurement data in the previous moment.
By removing data on flows after abnormal data, be accumulate to an authentication period to carry out verification process as follows:
The data on flows of judgement accumulation is from vapour source or user, if from vapour source, authentication method is as follows:
Adopt formula
Figure 785479DEST_PATH_IMAGE022
the integrated flux that Dui Qi source produces authenticates, and obtains flow authentication value
Figure 883885DEST_PATH_IMAGE023
, in formula
Figure 510038DEST_PATH_IMAGE024
=1 ~ m (in pipe network, steam source adds up to m), the numbering in expression vapour source.
In formula, ,
Figure 63697DEST_PATH_IMAGE026
,
Figure 282188DEST_PATH_IMAGE027
zhi Qi source, in an authentication period, adds integrated flux, equipment self use or the consumed flow of water and the integrated flux diffusing respectively; , ,
Figure 986205DEST_PATH_IMAGE030
be respectively the correction factor of above three kinds of flows to the actual flow that steams.
There for some steel enterprise, do not have to be right
Figure 996886DEST_PATH_IMAGE026
,
Figure 27159DEST_PATH_IMAGE027
the situation of metering, adopts experimental formula estimation:
Figure 123291DEST_PATH_IMAGE031
Figure 634563DEST_PATH_IMAGE032
In formula
Figure 765330DEST_PATH_IMAGE033
,
Figure 966504DEST_PATH_IMAGE034
be respectively the instant flow of two kinds of situations,
Figure 612249DEST_PATH_IMAGE035
,
Figure 861965DEST_PATH_IMAGE036
be respectively in the lasting time of this metering section.
Figure 909555DEST_PATH_IMAGE037
be respectively flow formula calculation constant, the gauge pressure in P Wei Qi source.
If the data on flows of accumulation is from user, authentication method is as follows:
If b the user that this measuring point is steam pipe system (b=1 ~ n, in pipe network, steam total number of users is n), total p of this user uses vapour point, and each flow measurements with vapour point is (j=1,2 ..., p)
Figure 405761DEST_PATH_IMAGE040
In formula,
Figure 307858DEST_PATH_IMAGE041
(j=1,2 ..., p) represent respectively with vapour point, to consider the flow modificatory coefficient after loss.
If sub-metering is inaccurate or while there is no sub-metering, adopt following method validation at all:
In formula
Figure 143276DEST_PATH_IMAGE043
represent this with vapour o'clock for the production of e(e=1,2 ..., the q) quantity of kind product,
Figure 62691DEST_PATH_IMAGE044
represent the quantity of steam of the e kind product consumption of every production unit quantity.
The data of accumulation and verify data are contrasted, if the two exists significant difference, non intelligent flowmeter error are carried out to trimming process and be:
There is the situation identification of directional bias:
When
Figure 756977DEST_PATH_IMAGE045
, and the difference of measured value and authentication value is always greater than zero or while being less than zero, thinks and have directional bias that u=1 ~ n+m in formula represents the outer contact numbering of steam pipe system.
Figure 530243DEST_PATH_IMAGE046
be the measured value of u outer contact at the accumulated value of an authentication period, authentication value for this node flow accumulated value.
The method of estimation of directional bias Wu:
Suppose u the continuous H of a node China Metrology Accreditation cycle T to pipe network, measured value is
Figure 33085DEST_PATH_IMAGE048
, authentication value is
Figure 644195DEST_PATH_IMAGE049
, h=1,2 ..., H
Figure 794554DEST_PATH_IMAGE050
The accumulated value of this flow after secondary correction is:
Figure 61587DEST_PATH_IMAGE051
To the data on flows after proofreading and correct, according to producing method that consumption balance coordinates by data, to carry out fine correction process as follows:
1) by outer contact numbering, and form incidence matrix a;
　 A＝（1,...,-1,…）；
In A, 1 number equals vapour source sum m;-1 number equals the total n with vapour point; 2) rule of thumb determine steam pipe system loss per hour
Figure 57225DEST_PATH_IMAGE052
and total losses
Figure 726103DEST_PATH_IMAGE053
:
Figure 47363DEST_PATH_IMAGE054
, tfor authentication period;
3) data are coordinated:
31) establishing each outer contact flow measurements of steam pipe system is X(u, y), unit: ton/hour;
U=1 ~ m+n, the outer contact numbering of steam pipe system.Y=1 ~ Z, represents y measured value in a metering cycle, and each authentication period is measured Z time;
If instrument fault or occur exceptional value situation, presses after preceding method is processed and substitutes measured value by estimated value.
32) calculate the instant flow sample average of the outer contact of u
Figure 601022DEST_PATH_IMAGE056
33) calculate in a metering cycle integrated flux of steam pipe system u outer contact:
Figure 124407DEST_PATH_IMAGE057
In formula tsfor traffic sampling cycle, unit hour; If authentication is not passed through, with the secondary correction integrated flux after authentication, substitute this value.
Note
Figure 353918DEST_PATH_IMAGE058
for the cumulative flow value vector after secondary correction;
Note
Figure 657861DEST_PATH_IMAGE059
for the flow vector after data are coordinated:
The objective function that fetches data and coordinate:
Figure 198564DEST_PATH_IMAGE060
Constraint condition is:
Figure 638772DEST_PATH_IMAGE061
Get weighting matrix:
Figure 239518DEST_PATH_IMAGE062
In formula:
Figure 765177DEST_PATH_IMAGE063
, the precision of the measurement that R (u) is u flowmeter;
34) the poor d of flow totalizer volume production consumption:
Figure 109571DEST_PATH_IMAGE064
35) calculate matrix:
Figure 238250DEST_PATH_IMAGE066
36) calculate Lagrange multiplier battle array K:
Figure 454467DEST_PATH_IMAGE067
37) computational data is coordinated residual matrix
Figure 399290DEST_PATH_IMAGE068
:
Figure 751774DEST_PATH_IMAGE069
38) computational data is coordinated rear flow :
Figure 928994DEST_PATH_IMAGE071
The data coordination value of trying to achieve is carried out to boundary value verification and equation of constraint verification: verification is not passed through as boundary value, reset Weighting Matrices; As reset, can't pass through, again authenticate flow and pipe network waste.
The present invention has following beneficial effect and advantage:
1. the inventive method can, when flowmeter breaks down, be estimated the data of this measuring point;
2. the inventive method can be found in time in measured value indivedual exceptional values and be removed.
3. the inventive method can be according to a China Metrology Accreditation cycle, and the data of cummulative metering are implemented to authentication according to the production process of steam and the method that is associated with steel manufacture process.Analyze the deviation of the result of metering result and authentication, the appreciable error that may exist while determining flowmeter survey (directional bias), and in metering, remove appreciable error.
4. for authentication period, occur total generation of steam and the difference between total flow, the method that adopts data to coordinate is rationally removed.
5. the inventive method has improved precision, reliability and the rationality ，Dui iron and steel enterprise raising vapor management level of steam-flow meter value, reduces energy consumption significant.
Accompanying drawing explanation
Fig. 1 is typical iron and steel enterprise steam pipe system observing and controlling structural drawing.
Fig. 2 is principle of work sketch of the present invention.
Fig. 3 is that flowmeter state judgement working properly, estimation and exceptional value are removed journey figure.
Fig. 4 is the elimination process flow diagram of authentication and the appreciable error (directional bias) of steam integrated flux.
Fig. 5 is the data tuning algorithm process flow diagram that produces consumption difference.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further detail.
As shown in Figure 1, the TT&C system of iron and steel enterprise's steam pipe system is mainly to take Distributed Control System (DCS) as basis at present, the state variable relevant with steam pipe system (temperature, pressure, flow) and control variable (valve opening) are passed through to data acquisition unit and industrial network from the flowmeter side of pipe network, be transferred to energy management center or measurement centre.And administrative center or measurement centre can need to be sent to scene by the operational order of the setting of sensing transducer, diagnostic instruction and electric control valve by network and controller according to what produce.
The informationization of steam pipe system is to implement basic condition of the present invention.
As shown in Figure 2, the present invention is applicable to the method that iron and steel enterprise's steam pipe system continuous data is proofreaied and correct, and comprises the following steps:
Judge that whether flowmeter duty is normal;
If flowmeter duty is normal, judge in the data on flows that flowmeter receives whether have abnormal data;
If there is abnormal data, remove;
The data on flows of removing after abnormal data is accumulate to an authentication period, authenticates, obtain verify data;
The data of accumulation and verify data are contrasted, if the two exists significant difference, flowmeter error is proofreaied and correct;
The method that data on flows after proofreading and correct is coordinated by data according to product consumption balance is carried out fine correction, obtains rational steam pipe system continuous data.
As shown in Figure 3, the present invention is according to on-site experience, and the data characteristics that flowmeter is transferred to measurement centre (or command centre) when improper metering shows the curvilinear characteristic that comes, design discharge meter state discriminant function, under function meets certain condition, automatically judge that this flowmeter metering is abnormal.Now the reading of flowmeter is unreliable, need to estimate data on flows.
According to on-site experience, when flowmeter is abnormal, main manifestations goes out two kinds of features, and a kind of is exactly that data do not change (being generally to be caused by communication disruption) for a long time, and another kind is that measured value is in short time fluctuation (sensor or flowmeter fault) up and down by a relatively large margin.In addition, for pipe network a certain vapour source or user, the statistical law during according to normal operation, is removed some obvious abnormal data.
1. the automatic judging method of flow flowmeter operation irregularity state.
If the actual value of the instant flow variable that a certain flowmeter is measured is X(t), to this measuring point continuous acquisition M time (sample frequency 1Hz, or select according to actual conditions) data, complete a whether normal duty of flowmeter and identify.
Definition: the absolute undulate quantity of data on flows is:
Figure 698015DEST_PATH_IMAGE009
Wherein
Figure 967322DEST_PATH_IMAGE010
(r=1 ~ M) is from current time, the data on flows sequence that is M to length by some cycles continuous acquisition; R is integer, represents the sequence number of sampled value, the length that M is sample sequence;
When flowmeter works, for different measurement points, V value is specific interval in one.
That is:
Figure 80772DEST_PATH_IMAGE072
In formula,
Figure 68319DEST_PATH_IMAGE011
with
Figure 292627DEST_PATH_IMAGE012
while working for flowmeter, the minimum of V variation range and maximum statistical value.
When time, think this flowmeter communication disruption.
When
Figure 966371DEST_PATH_IMAGE014
time, think this flowmeter fault.
In order to reduce the probability of erroneous judgement,
Figure 237953DEST_PATH_IMAGE011
with need to determine according to on-the-spot actual operating data.
2. monolithic flowmeter is extremely or in short-term during communication disruption, the method for estimation to the instant data on flows of this measuring point.
According to pipe network structure, determine and flowmeter in normal metering state nearest with this flowmeter measuring point, be designated as REF;
If:
Figure 509851DEST_PATH_IMAGE001
In formula,
Figure 965103DEST_PATH_IMAGE002
for the flow actual value of this flowmeter position,
Figure 927243DEST_PATH_IMAGE003
the flow actual value of WeiREF position, k is correlation coefficient.Choose available N history samples value of up-to-date preservation, adopt the correlation coefficient k in least square method estimation formula.
Figure 493353DEST_PATH_IMAGE004
Wherein i=1 ~ N is integer, represents sampled value sequence number; N is the length of sample sequence.
Figure 594689DEST_PATH_IMAGE005
for this flowmeter i flow measurements in N history samples value,
Figure 283159DEST_PATH_IMAGE006
for i the measured value of REF N history samples value;
The estimated value of this measurement point present flow rate is:
In formula
Figure 836818DEST_PATH_IMAGE008
measured value for current sampling instant REF.
For long-time flowmeter, abnormal or communication disruption, just adopts the flow authentication method in the present invention, uses authentication value as measured value.
3. the identification of exceptional value and removal in measured value
When flowmeter works, there will be and be subject to the point that outer signals disturbing effect obviously departs from other measured values, be i.e. so-called wild value.
If the measured value of t sampling instant is
Figure 992992DEST_PATH_IMAGE015
, normal
Figure 852364DEST_PATH_IMAGE015
should meet the following conditions:
i)
Figure 461200DEST_PATH_IMAGE016
;
ii)
Figure 431430DEST_PATH_IMAGE017
；
In formula,
Figure 504428DEST_PATH_IMAGE018
,
Figure 472384DEST_PATH_IMAGE019
represent respectively the minimum and maximum value that this measurement point steam flow can reach in technological design; represent the sampling period;
Figure 342437DEST_PATH_IMAGE021
be illustrated in the maximum possible rate of change of steam flow in the sampling period;
As measured value is
Figure 535521DEST_PATH_IMAGE015
i satisfies condition when different), ii), think exceptional value, this measuring point data is removed, and by the alternative current measurement data of measurement data in the previous moment.
Above condition shows, to the numerical value of steam pipe system flow measurement, between a set scope, and rate of change is also subject to the restriction of pipe network and technique must be less than a definite value.
One of above condition does not meet, and this measuring point data can be removed, and substitute current measurement data by the measurement data in the previous moment.
For reducing the probability of erroneous judgement,
Figure 674378DEST_PATH_IMAGE018
,
Figure 51614DEST_PATH_IMAGE019
,
Figure 566909DEST_PATH_IMAGE021
three parameters need to be according to historical data and experience adjustments.
As shown in Figure 4, the elimination process flow diagram of the authentication of steam integrated flux and appreciable error (directional bias).
By differentiation and data estimation to flowmeter abnormal work state, obtain the good data of availability, by these data, accumulate, just can obtain the integrated flux of China Metrology Accreditation cycle (as an order of classes or grades at school).
But in actual applications, often occur that some flowmeter, when measuring steam flow, exists appreciable error unidirectional and that numerical value is substantially constant always between measured value and actual value, this appreciable error is a kind of directional bias.
The measurement point of u steam flow on steam pipe system (the vapour source of steam or the user of steam) at an authentication period integrated flux actual value is
Figure 614500DEST_PATH_IMAGE073
, measured value is
Figure 986575DEST_PATH_IMAGE074
.Have
Figure 791720DEST_PATH_IMAGE075
In formula,
Figure 907444DEST_PATH_IMAGE076
for directional bias,
Figure 747224DEST_PATH_IMAGE077
for stochastic error, and:
；
Figure 848221DEST_PATH_IMAGE079
；
Figure 502056DEST_PATH_IMAGE080
Wu is the directional bias between measured value and actual value.
When iron and steel enterprise normally produces, according to certain measurement period, record its measured value series
Figure 461922DEST_PATH_IMAGE074
(u=1,2 ..., m+n), and right by the following method (i=1,2 ..., m+n) authenticate.
1. the feed flow authentication of steam source to pipe network
Steam source feed flow authentication formula:
Figure 17854DEST_PATH_IMAGE022
The integrated flux that Dui Qi source produces authenticates, and obtains flow authentication value
Figure 413063DEST_PATH_IMAGE023
, in formula
Figure 292682DEST_PATH_IMAGE024
=1 ~ m (in pipe network, steam source adds up to m), the numbering in expression vapour source.
In formula,
Figure 115145DEST_PATH_IMAGE025
,
Figure 444495DEST_PATH_IMAGE026
,
Figure 643395DEST_PATH_IMAGE027
zhi Qi source, in an authentication period, adds integrated flux, equipment self use or the consumed flow of water and the integrated flux diffusing respectively;
Figure 374591DEST_PATH_IMAGE028
,
Figure 367954DEST_PATH_IMAGE029
,
Figure 184601DEST_PATH_IMAGE030
be respectively the correction factor of above three kinds of flows to the actual flow that steams.
There for some steel enterprise, do not have to be right
Figure 249509DEST_PATH_IMAGE026
,
Figure 772894DEST_PATH_IMAGE027
the situation of metering, can adopt experimental formula estimation.
Figure 999476DEST_PATH_IMAGE031
Figure 241101DEST_PATH_IMAGE032
In formula ,
Figure 222013DEST_PATH_IMAGE034
be respectively the instant flow of two kinds of situations,
Figure 885075DEST_PATH_IMAGE035
,
Figure 407805DEST_PATH_IMAGE036
be respectively in the lasting time of this metering section.
Figure 752198DEST_PATH_IMAGE037
Figure 46914DEST_PATH_IMAGE038
be respectively flow formula calculation constant.The gauge pressure in P Wei Qi source.
In above formula, ,
Figure 831516DEST_PATH_IMAGE038
under the condition that need to can correctly measure in pressure, flow and time, determine.It is pointed out that
Figure 776338DEST_PATH_IMAGE037
,
Figure 128822DEST_PATH_IMAGE038
can there is little variation with temperature, the pressure of steam.When integrated flux authenticates, its impact can be ignored.
2. the authentication of steam pipe system to user's operational throughput
If b the user that this measuring point is steam pipe system (b=1 ~ n, in pipe network, steam total number of users is n), total p of this user uses vapour point, and each flow measurements with vapour point is
Figure 133687DEST_PATH_IMAGE039
(j=1,2 ..., p)
Figure 571622DEST_PATH_IMAGE040
In formula,
Figure 320135DEST_PATH_IMAGE041
(j=1,2 ..., p) represent respectively with vapour point, to consider the flow modificatory coefficient after loss.
If sub-metering is inaccurate or while there is no sub-metering, adopt following method validation at all:
Figure 527125DEST_PATH_IMAGE042
In formula
Figure 702892DEST_PATH_IMAGE043
represent this with vapour o'clock for the production of e(e=1,2 ..., the q) quantity of kind product,
Figure 628122DEST_PATH_IMAGE044
represent the quantity of steam of the e kind product consumption of every production unit quantity.
In order to improve reliability and the rationality of authentication, , all need to, at flowmeter under the condition through calibration, in conjunction with historical data and experience, confirm.
3. the identification of directional bias and removal
According to the practical experience of metering, when occurring larger difference at certain section of time measured value and authentication value, and this species diversity is unidirectional all the time, and measured value is always greater than authentication value, or is always less than authentication value, it is generally acknowledged and has directional bias, be removed.
When
Figure 388158DEST_PATH_IMAGE081
, and the difference of measured value and authentication value is always greater than zero or while being less than zero, thinks and have directional bias that u=1 ~ n+m in formula represents the outer contact numbering of steam pipe system.
Figure 863002DEST_PATH_IMAGE046
be the measured value of u outer contact at the accumulated value of an authentication period,
Figure 891001DEST_PATH_IMAGE047
authentication value for this node flow accumulated value.
The method of estimation of directional bias Wu:
Suppose u the continuous H of a node China Metrology Accreditation cycle T to pipe network, measured value is
Figure 134901DEST_PATH_IMAGE048
, authentication value is , h=1,2 ..., H
Figure 286713DEST_PATH_IMAGE050
The accumulated value of this flow after secondary correction is:
Figure 852824DEST_PATH_IMAGE051
to assert the important parameter that whether has appreciable error between measured value and authentication value, generally need to be according to on-site experience, the factors such as the precision of instrument are determined between 10%-20%.
For intelligent flow meter, can remove measurement deviation by biasing setting, to non intelligent flowmeter, on production commander or measurement centre's computing machine, process.
According to order of classes or grades at school or some cycles (China Metrology Accreditation cycle T hour) as shown in Figure 5, the principle according to total generation with wastage in bulk or weight balance, implementation data is coordinated (adjustment) algorithm flow chart.
The method design that this method solves according to least-squares estimation and Lagrange multiplier, its method step is expressed as follows:
1. user's input parameter
1) externally contact is numbered, and forms incidence matrix A by certain requirement.
To each vapour source with by vapour point numbering, and require vapour source front during numbering, with vapour point rear.
Obtain matrix A=(1 ... ,-1 ...)
In A, 1 number equals vapour source sum m.
In A ,-1 number equals the total n with vapour point.
2) steam pipe system loss per hour
Figure 905279DEST_PATH_IMAGE052
and total losses
Figure 26819DEST_PATH_IMAGE053
:
Be the condensate water of steam pipe system generation per hour, the leakage of pipe network, the normal amount of diffusing.Unit: ton/hour.The foundation of this value: the calculating of pipe network model, historical data statistics, experience are estimated.
Figure 458937DEST_PATH_IMAGE054
tfor authentication period.
2. data coordination process.
1) establishing each outer contact flow measurements of steam pipe system is X(u, y), unit: ton/hour;
U=1 ~ m+n, the outer contact numbering of steam pipe system.Y=1 ~ Z, represents y measured value in a metering cycle, and each authentication period is measured Z time;
If instrument fault or occur exceptional value situation, presses after preceding method is processed and substitutes measured value by estimated value.
2) calculate the instant flow sample average of the outer contact of u
Figure 736816DEST_PATH_IMAGE055
3) in a metering cycle, the integrated flux of steam pipe system u outer contact:
Figure 267341DEST_PATH_IMAGE057
In formula, Ts is the traffic sampling cycle, unit hour.
In formula tsfor traffic sampling cycle, unit hour; If authentication is not passed through, with the secondary correction integrated flux after authentication, substitute this value.
Note
Figure 175254DEST_PATH_IMAGE058
for the cumulative flow value vector after secondary correction;
Note
Figure 248252DEST_PATH_IMAGE059
for the flow vector after data are coordinated:
The objective function that fetches data and coordinate:
Figure 216208DEST_PATH_IMAGE083
Constraint condition is:
Get weighting matrix:
Figure 86261DEST_PATH_IMAGE062
In formula:
Figure 13766DEST_PATH_IMAGE063
, R (u), for for to u flowmeter dimension inspection experience, calibration and authentication, confirms the precision (relative error) of the measurement of this flowmeter.
Solving data coordination value (essence calibration) is to solve under constraint condition, and making objective function is minimum least square Solve problems.The process solving is according to the 4th) to 11) step implements according to this.
4) the poor d of flow totalizer volume production consumption:
Figure 214940DEST_PATH_IMAGE064
　
5) calculate
Figure 798368DEST_PATH_IMAGE065
matrix:
Figure 378910DEST_PATH_IMAGE066
6) Lagrange multiplier battle array K:
Figure 98604DEST_PATH_IMAGE067
7) data are coordinated residual matrix
Figure 470679DEST_PATH_IMAGE068
:
Figure 541404DEST_PATH_IMAGE069
8) data are coordinated rear flow
Figure 657127DEST_PATH_IMAGE070
:
Figure 559224DEST_PATH_IMAGE071
9) the data coordination value of trying to achieve is carried out to boundary value verification and equation of constraint verification:
As verification is not passed through, need Weighting Matrices to reset.
As reset, can't pass through, again authenticate flow and pipe network waste.

Claims (11)

1. be applicable to the method that iron and steel enterprise's steam pipe system continuous data is proofreaied and correct, it is characterized in that comprising the following steps:
Judge that whether flowmeter duty is normal;
If flowmeter duty is normal, judge in the data on flows that flowmeter receives whether have abnormal data;
If there is abnormal data, remove;
The data on flows of removing after abnormal data is accumulate to an authentication period, authenticates, obtain verify data;
The data of accumulation and verify data are contrasted, if the two exists significant difference, flowmeter error is proofreaied and correct;
The method that data on flows after proofreading and correct is coordinated by data according to product consumption balance is carried out fine correction, obtains rational steam pipe system continuous data;
If flowmeter working state abnormal, estimates flowmeter survey value, obtain the estimated value of data on flows, continue and judge in the data on flows that flowmeter receives whether have abnormal data step;
It is described that flowmeter survey value is carried out to estimation procedure is as follows:
According to pipe network structure, determine and flowmeter in normal metering state nearest with this flowmeter measuring point, be designated as REF;
If:
Figure FDA0000385312280000011
In formula,
Figure FDA0000385312280000012
for the flow actual value of this flowmeter position, the flow actual value of WeiREF position, k is correlation coefficient; Choose available N history samples value of up-to-date preservation and calculate k
Figure FDA0000385312280000014
Wherein i=1～N is integer, represents sampled value sequence number; N is the length of sample sequence, X ifor this flowmeter i flow measurements in N history samples value, X rEFifor i the measured value of REF N history samples value;
The estimated value of this measurement point present flow rate is:
Figure FDA0000385312280000015
X in formula rEFmeasured value for current sampling instant REF.
2. by the method that is applicable to iron and steel enterprise's steam pipe system continuous data correction claimed in claim 1, it is characterized in that: the data of accumulation and verify data are contrasted, if there is not significant difference in the two, the method of directly data on flows being coordinated by data according to product consumption balance is carried out fine correction, obtains rational steam pipe system continuous data.
3. by claimed in claim 1, be applicable to the method that iron and steel enterprise steam pipe system continuous data is proofreaied and correct, it is characterized in that: if there is no abnormal data, continue data on flows is accumulate to an authentication period, authenticate, obtain verify data step.
4. by claimed in claim 1, be applicable to the method that iron and steel enterprise steam pipe system continuous data is proofreaied and correct, it is characterized in that: describedly judge flowmeter duty whether normal processes is as follows;
The absolute undulate quantity of definition data on flows is:
Figure FDA0000385312280000021
X wherein r(r=1～M) is from current time, the data on flows sequence that is M to length by some cycles continuous acquisition; R is integer, represents the sequence number of sampled value, the length that M is sample sequence;
During flowmeter normal operation, minimum and the maximum statistical value of V variation range are V land V h;
As V<V<sub TranNum="525">l</sub>time, think this flowmeter communication disruption;
Work as V>V htime, think this flowmeter fault.
5. by the method that is applicable to iron and steel enterprise's steam pipe system continuous data correction claimed in claim 1, it is characterized in that: judge whether that the process that has abnormal data and remove is as follows:
If the measured value of t sampling instant is X t, normal X tshould meet the following conditions:
i)X L≤X t≤X H;
ii)
Figure FDA0000385312280000022
In formula, X l, X hrepresent respectively the minimum and maximum value that this measurement point steam flow can reach in technological design; T srepresent the sampling period; Δ is illustrated in the maximum possible rate of change of steam flow in the sampling period;
If measured value is X ti satisfies condition when different), ii), think exceptional value, this measuring point data is removed, and by the alternative current measurement data of measurement data in the previous moment.
6. by claimed in claim 1, be applicable to the method that iron and steel enterprise steam pipe system continuous data is proofreaied and correct, it is characterized in that being accumulate to an authentication period to carry out verification process as follows removing data on flows after abnormal data:
The data on flows of judgement accumulation is from vapour source or user, if from vapour source, authentication method is as follows:
Adopt formula
Figure FDA0000385312280000023
the integrated flux that Dui Qi source produces authenticates, and obtains flow authentication value
Figure FDA0000385312280000024
, a=1～m in formula, the numbering in expression vapour source, in pipe network, steam source adds up to m;
In formula, F w, F c, F rzhi Qi source, in an authentication period, adds integrated flux, equipment self use or the consumed flow of water and the integrated flux diffusing respectively; k 1, k 2, k 3be respectively the correction factor of above three kinds of flows to the actual flow that steams.
7. by the method that is applicable to iron and steel enterprise's steam pipe system continuous data correction claimed in claim 6, it is characterized in that:
For some steel, look forward to not to F c, F rthe situation of metering, adopts experimental formula estimation:
Figure FDA0000385312280000031
Figure FDA0000385312280000032
X in formula c, X rbe respectively the instant flow of two kinds of situations, t c, t rbe respectively at lasting time of this metering section, k ck rbe respectively flow formula calculation constant, the gauge pressure in P Wei Qi source.
8. by the method that is applicable to iron and steel enterprise's steam pipe system continuous data correction claimed in claim 6, it is characterized in that: if the data on flows of accumulation is from user, authentication method is as follows:
If b the user that measuring point is steam pipe system, b=1～n, in pipe network, steam total number of users is n, it is F that this user total p uses vapour point, each flow measurements with vapour point bj(j=1,2 ..., p)
Figure FDA0000385312280000033
In formula,
Figure FDA0000385312280000034
for p the average by the flow measurements of vapour point, k j(j=1,2 ..., p) represent respectively with vapour point, to consider the flow modificatory coefficient after loss.
9. by the method that is applicable to iron and steel enterprise's steam pipe system continuous data correction claimed in claim 8, it is characterized in that:
If sub-metering is inaccurate or while there is no sub-metering, adopt following method validation at all:
Figure FDA0000385312280000035
M in formula berepresent this with vapour o'clock for the production of e, e=1,2 ..., q, plants the quantity of product, λ erepresent the quantity of steam of the e kind product consumption of every production unit quantity.
10. by the method that is applicable to iron and steel enterprise's steam pipe system continuous data correction claimed in claim 1, it is characterized in that: the data of accumulation and verify data are contrasted, if the two exists significant difference, non intelligent flowmeter error is carried out to trimming process and is:
There is the situation identification of directional bias:
When
Figure FDA0000385312280000036
and the difference of measured value and authentication value is always greater than zero or while being less than zero, thinks and have directional bias that u=1～n+m in formula represents the outer contact numbering of steam pipe system, F ube the measured value of u outer contact at the accumulated value of an authentication period,
Figure FDA00003853122800000310
authentication value for this node flow accumulated value;
Directional bias W umethod of estimation:
Suppose u the continuous H of a node China Metrology Accreditation cycle T to pipe network, measured value is F uh, authentication value is
Figure FDA0000385312280000037
h=1,2 ..., H
Figure FDA0000385312280000038
The accumulated value of this flow after secondary correction is:
Figure FDA0000385312280000039
11. are applicable to by claimed in claim 1 the method that iron and steel enterprise steam pipe system continuous data is proofreaied and correct, and it is characterized in that: to the data on flows after proofreading and correct, according to producing method that consumption balance coordinates by data, to carry out fine correction process as follows:
1) by outer contact numbering, and form incidence matrix A;
A＝（1,...,-1,...）；
In A, 1 number equals vapour source sum m;-1 number equals the total n with vapour point;
2) rule of thumb determine steam pipe system loss σ per hour and total losses Σ:
Σ=σ * T, T is authentication period;
3) data are coordinated:
31) establishing each outer contact flow measurements of steam pipe system is X(u, y), unit: ton/hour;
U=1～m+n, the outer contact numbering of steam pipe system, y=1～Z, represents y measured value in a metering cycle, each authentication period is measured Z time;
If instrument fault or occur exceptional value situation, presses after preceding method is processed and substitutes measured value by estimated value;
32) calculate the outer instant flow sample average of contact of u
Figure FDA0000385312280000041
Figure FDA0000385312280000042
33) calculate in a metering cycle integrated flux of steam pipe system u outer contact:
Figure FDA0000385312280000043
In formula, Ts is the traffic sampling cycle, and unit is hour; If authentication is not passed through, with the secondary correction integrated flux after authentication, substitute this value;
Note for the cumulative flow value vector after secondary correction;
Note
Figure FDA0000385312280000045
for the flow vector after data are coordinated:
The objective function that fetches data and coordinate:
Figure FDA00003853122800000411
Constraint condition is:
Get weighting matrix: Q=diag (β (1), β (2) ..., β (m+n))
In formula:
Figure FDA0000385312280000047
the precision of the measurement that R (u) is u flowmeter;
34) the poor d:d=AF-Σ of flow totalizer volume production consumption
35) calculate Q ematrix: Q e=AQA t
36) calculate Lagrange multiplier battle array K:
Figure FDA0000385312280000048
37) computational data is coordinated residual matrix Γ: Γ=QA tk
38) computational data is coordinated rear flow
Figure FDA0000385312280000049
:
Figure FDA00003853122800000410
The data coordination value of trying to achieve is carried out to boundary value verification and equation of constraint verification: verification is not passed through as boundary value, reset weighting matrix; As reset, can't pass through, again authenticate flow and pipe network waste.
CN201010572857.2A 2010-12-05 2010-12-05 Method for steam pipe network of iron and steel enterprise to correct metering data Expired - Fee Related CN102486390B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201010572857.2A CN102486390B (en) 2010-12-05 2010-12-05 Method for steam pipe network of iron and steel enterprise to correct metering data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010572857.2A CN102486390B (en) 2010-12-05 2010-12-05 Method for steam pipe network of iron and steel enterprise to correct metering data

Publications (2)

Publication Number Publication Date
CN102486390A CN102486390A (en) 2012-06-06
CN102486390B true CN102486390B (en) 2014-01-15

Family

ID=46151912

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010572857.2A Expired - Fee Related CN102486390B (en) 2010-12-05 2010-12-05 Method for steam pipe network of iron and steel enterprise to correct metering data

Country Status (1)

Country Link
CN (1) CN102486390B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104931722B (en) * 2015-03-23 2018-02-16 杭州电子科技大学 A kind of measuring circuit of drainage pipeline networks flow velocity
CN104950852A (en) * 2015-06-16 2015-09-30 河南中烟工业有限责任公司 Data verification system and data verification method of energy management and control system of tobacco enterprise
CN105426429B (en) * 2015-11-04 2019-03-26 中国联合网络通信集团有限公司 Data processing method, induction element data processing equipment, data processing system
CN106168771A (en) * 2016-07-08 2016-11-30 江苏永钢集团有限公司 Analog signals identifing source system
CN109932024B (en) * 2017-12-15 2020-11-17 南京南瑞继保电气有限公司 Remote fault judgment method and system for industrial flow instrument
CN108763774B (en) * 2018-05-30 2022-06-07 新奥泛能网络科技有限公司 Heat loss calculation method and device for overhead steam pipe network
CN109684605B (en) * 2018-12-25 2023-05-23 阳谷祥光铜业有限公司 Data error processing method and device for metal balance
CN109686412B (en) * 2018-12-25 2023-05-26 阳谷祥光铜业有限公司 Data coordination processing method and device for metal balance
CN110375787B (en) * 2019-07-25 2021-10-15 重庆市计量质量检测研究院 Online evaluation method for running state of metering device
CN111324072A (en) * 2020-03-05 2020-06-23 烨辉（中国）科技材料有限公司 Energy consumption monitoring system of steel strip production line
CN112380491A (en) * 2020-11-16 2021-02-19 中冶赛迪技术研究中心有限公司 Calculation method for fluid pipe network instrument inspection

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000162012A (en) * 1998-09-22 2000-06-16 Akitoshi Kitano Method for correcting instrumental error of volume type flowmeter
CN1694109B (en) * 2005-06-20 2011-05-18 浙江大学 Material data correction method in chemical and oil refinement process
CN101286059B (en) * 2008-05-29 2010-06-02 上海交通大学 On-line measurement data obtaining, correction and gross error serial compensation process
CN101424550B (en) * 2008-09-03 2010-06-23 中国海洋石油总公司 Instrument meter freezing fault rapid detecting method
CN101477369B (en) * 2009-01-19 2011-01-05 浙江大学 Data correction method used for varying data of chemical engineering flow system

Also Published As

Publication number Publication date
CN102486390A (en) 2012-06-06

Similar Documents

Publication Publication Date Title
CN102486390B (en) Method for steam pipe network of iron and steel enterprise to correct metering data
CN109813400B (en) Boiler main steam flowmeter online checking and testing system based on internet of things technology
US7937247B2 (en) Virtual metering
CN110375787B (en) Online evaluation method for running state of metering device
CN108197725B (en) Water demand prior information-based water supply network node water demand checking method
CN108593055B (en) Online automatic calibration method and system for pulverized coal mass flow meter
CN114047472B (en) Metering error monitoring system of intelligent ammeter and monitoring method and device thereof
CN104764504B (en) A kind of flow correction method of saturation and superheated steam
CN104748810A (en) Fluid flow accumulating system
CN105316613A (en) Zinc coating thickness control method and system based on time-varying delay deviation correction technology
CN104515566A (en) Steam drum liquid level measuring device, system and method
CN103175582B (en) Fluid metering method and apparatus
CN108106679B (en) Method and system for measuring inlet air volume of coal mill of power station
CN111444478B (en) Comprehensive plant power consumption statistics method
CN210487022U (en) Steam flow measurement check-up equipment
CN102999028B (en) Gas medium continuous data early warning system digital-to-analogue constructing method
CN114937512A (en) Method and system for flow compensation of coolant of nuclear power unit primary loop
CN115219106A (en) Compressed air pipe network leakage dynamic measurement method based on cloud computing
CN106893793A (en) The monitoring method of real-time coal powder injection rate
CN104165711A (en) Verification table temperature measuring method based on Kalman filtering
CN102645315A (en) Automatic, fast and accurate detection method for air resistance characteristics of large heat exchanger
WO2023135876A1 (en) Measuring instrument inaccuracy assessment system and measuring instrument inaccuracy assessment method
CN110161996A (en) Method and system for the analysis of power plant units consumption
CN113340360B (en) Secondary network system balance flow measurement method
CN204556043U (en) A kind of fluid flow accumulation system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140115

Termination date: 20171205