CN117686140A - Method and system for evaluating measurement uncertainty of empty box barometer - Google Patents
Method and system for evaluating measurement uncertainty of empty box barometer Download PDFInfo
- Publication number
- CN117686140A CN117686140A CN202410145927.8A CN202410145927A CN117686140A CN 117686140 A CN117686140 A CN 117686140A CN 202410145927 A CN202410145927 A CN 202410145927A CN 117686140 A CN117686140 A CN 117686140A
- Authority
- CN
- China
- Prior art keywords
- uncertainty
- temperature coefficient
- verification
- box
- indication
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000005259 measurement Methods 0.000 title claims abstract description 185
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000012795 verification Methods 0.000 claims abstract description 226
- 238000011156 evaluation Methods 0.000 claims abstract description 60
- 230000000694 effects Effects 0.000 claims description 18
- 230000035945 sensitivity Effects 0.000 claims description 18
- 230000015572 biosynthetic process Effects 0.000 claims description 16
- 238000003786 synthesis reaction Methods 0.000 claims description 16
- 238000003556 assay Methods 0.000 claims description 11
- 238000003860 storage Methods 0.000 claims description 9
- BJQHLKABXJIVAM-UHFFFAOYSA-N bis(2-ethylhexyl) phthalate Chemical compound CCCCC(CC)COC(=O)C1=CC=CC=C1C(=O)OCC(CC)CCCC BJQHLKABXJIVAM-UHFFFAOYSA-N 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 7
- 238000011056 performance test Methods 0.000 claims description 7
- 238000012360 testing method Methods 0.000 description 22
- 238000004364 calculation method Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 8
- 238000009826 distribution Methods 0.000 description 8
- 238000009529 body temperature measurement Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 6
- 238000009827 uniform distribution Methods 0.000 description 5
- 239000000203 mixture Substances 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 239000013307 optical fiber Substances 0.000 description 3
- 239000000523 sample Substances 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052697 platinum Inorganic materials 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
Landscapes
- Measuring Fluid Pressure (AREA)
Abstract
The invention provides a method and a system for evaluating measurement uncertainty of an empty box barometer, wherein the method obtains a temperature coefficient by carrying out temperature coefficient verification on the empty box barometer, wherein the temperature coefficient comprises an influence index of temperature fluctuation of a temperature coefficient verification box on the temperature coefficient and an influence index of temperature uniformity on the temperature coefficient; the uncertainty evaluation of the temperature coefficient is measured, and an uncertainty evaluation result of the temperature coefficient is obtained; performing indicating value verification on the air pressure gauge of the empty box to obtain an indicating value error, wherein the indicating value error comprises an influence index of the air leakage rate of the indicating value verification box on a measurement result and an influence index of fluctuation on the measurement result; and according to the uncertainty evaluation result of the temperature coefficient, the uncertainty evaluation of measurement is carried out on the indication value error, so that an uncertainty evaluation result of the measurement of the empty box barometer is obtained, and the accuracy of the uncertainty evaluation of the measurement of the empty box barometer is effectively improved.
Description
Technical Field
The invention belongs to the technical field of measurement uncertainty evaluation of an empty box barometer, and particularly relates to an empty box barometer measurement uncertainty evaluation method and an empty box barometer measurement uncertainty evaluation system.
Background
When the measurement uncertainty of the empty box barometer is evaluated, the empty box barometer is firstly verified according to the requirements of JJG 272-2007 'empty box barometer and empty box barometer' verification rules. The verification process comprises two steps: 1) Performing temperature coefficient verification on the empty box barometer to obtain a corresponding temperature coefficient; 2) And performing indication verification to obtain a corresponding indication error. In the prior researches, the uncertainty of measurement is evaluated only by specially aiming at the temperature coefficient of the empty box barometer; or only the final indication error of the air box barometer is evaluated, and the influence of the measurement uncertainty of the temperature coefficient on the final measurement result is not considered in the process. Therefore, the obtained final measurement uncertainty assessment results are difficult to accurately reflect the actual measurement reality.
Meanwhile, in the conventional measurement uncertainty study of the temperature coefficient of the empty box barometer, the influence of the performance (temperature uniformity and fluctuation) of the temperature coefficient verification box on the temperature coefficient verification result is generally not considered; in addition, the influence of the performance (air leakage rate and volatility) of the indication verification box on the indication verification is not considered in the measurement uncertainty evaluation study on the indication error of the empty box barometer. However, these factors have an effect on the measurement results of the air box barometer. The main reason is that the air pressure change is sensed by the empty box air pressure gauge through the empty box film due to the internal mechanical structure of the empty box air pressure gauge, and the pointer is driven by the lever to indicate the air pressure value.
Disclosure of Invention
Based on the above, the embodiment of the invention provides a method and a system for evaluating the measurement uncertainty of an empty box barometer, which aim to solve the problem of inaccurate evaluation of the measurement uncertainty of the empty box barometer in the prior art.
According to a first aspect of the embodiment of the invention, a performance test result of a temperature coefficient verification box and an indication verification box is introduced to complete measurement uncertainty assessment of the empty box barometer, wherein the temperature coefficient verification box and the indication verification box are both provided with barometer, at least three temperature measuring heads which are arranged in a step shape are arranged in the temperature coefficient verification box, and the indication verification box is internally provided with the temperature measuring heads which are positioned at a central point. The temperature measuring head in the temperature coefficient verification box is electrically connected with the sensor outside the temperature coefficient verification box to collect the temperatures of different points in the temperature coefficient verification box, and the temperature measuring head in the indication verification box is electrically connected with the sensor outside the indication verification box to collect the temperature of the center point in the indication verification box, and the method comprises the following steps:
placing the empty box barometer in a temperature coefficient verification box, and carrying out temperature coefficient verification on the empty box barometer to obtain a temperature coefficient, wherein the temperature coefficient comprises an influence index of temperature fluctuation of the temperature coefficient verification box on the temperature coefficient and an influence index of temperature uniformity of the temperature coefficient verification box on the temperature coefficient;
Evaluating the uncertainty of the temperature coefficient to obtain an uncertainty evaluation result of the temperature coefficient;
placing the empty box barometer in an indication verification box, and performing indication verification on the empty box barometer to obtain an indication error, wherein the indication error comprises an influence index of the air leakage rate of the indication verification box on a measurement result and an influence index of the fluctuation of the indication verification box on the measurement result;
and according to the uncertainty evaluation result of the temperature coefficient, evaluating the measurement uncertainty of the indication error to obtain an empty box barometer measurement uncertainty evaluation result.
Further, in the step of performing temperature coefficient verification on the air box barometer to obtain a temperature coefficient, an expression of a measurement model of the temperature coefficient is as follows:
；
wherein,K t represented as an empty box barometer temperature coefficient,p s1 represented as the average indication of the corrected n times barometric pressure etalon at the high temperature point,p s2 represented as the average indication of the corrected n times barometric pressure etalon at the low temperature point,p b1 expressed as the average indication of the n times empty box barometer at the high temperature point,p b2 represented as the average indication of the n empty box barometer at the low temperature point,t 1 represented as the average indication of the n times thermometer at the hot spot, t 2 Represented as the average indication of the low temperature point n times thermometer,q 1 expressed as the effect of temperature fluctuations in the temperature coefficient assay chamber on the temperature coefficient,q 2 the effect of temperature uniformity in the temperature coefficient assay chamber on the temperature coefficient is expressed.
Further, the step of placing the empty box barometer in an indication verification box and performing indication verification on the empty box barometer to obtain an indication error includes the following steps:
；
wherein DeltaP is the air pressure error value of the air pressure gauge of the empty box,P b represented as an indication of the empty box barometer,K t represented as an empty box barometer temperature coefficient,Trepresented by the temperature of the environment in question,P s indicated as an indication of the corrected barometric pressure standard,Q 1 expressed as the effect of the indicator verification box leak rate on the measurement results,Q 2 expressed as the effect of the fluctuation of the indicator verification box on the measurement result.
Further, the step of performing the assessment of the measurement uncertainty includes:
calculating sensitivity coefficients corresponding to the input quantities in the measurement model and standard uncertainty of the input quantities according to the measurement model;
determining each standard uncertainty component according to the sensitivity coefficient and the standard uncertainty corresponding to each input quantity;
calculating the uncertainty of the synthetic standard according to the uncertainty components of each standard;
And calculating the extension uncertainty according to the synthesis standard uncertainty.
Further, when the temperature coefficient is evaluated for measurement uncertainty, the calculated standard uncertainty of each input quantity comprises the standard uncertainty corresponding to an air pressure standard in the temperature coefficient verification box, the standard uncertainty corresponding to an empty box air pressure meter in the temperature coefficient verification box, the standard uncertainty corresponding to the thermometer in the temperature coefficient verification box, the standard uncertainty corresponding to the temperature fluctuation in the temperature coefficient verification box and the standard uncertainty corresponding to the temperature uniformity in the temperature coefficient verification box;
and when the indication error is evaluated, the calculated standard uncertainty of each input quantity comprises the standard uncertainty corresponding to an air pressure standard in the indication verification box, the standard uncertainty corresponding to an empty box air pressure meter in the indication verification box, the standard uncertainty corresponding to a temperature coefficient, the standard uncertainty corresponding to a thermometer at the central point in the indication verification box, the standard uncertainty corresponding to the air leakage rate of the indication verification box and the standard uncertainty corresponding to the fluctuation of the indication verification box.
Further, in the step of calculating the uncertainty of the synthetic standard according to each standard uncertainty component, the expression for calculating the uncertainty of the synthetic standard is:
；
Wherein c i The sensitivity coefficient expressed as the i-th input quantity, u (x i ) Standard uncertainty expressed as the i-th input quantity, u c Expressed as synthesis standard uncertainty.
Further, in the step of calculating the extended uncertainty according to the synthesis standard uncertainty, the expression for calculating the extended uncertainty is:
；
or (b)
；
Wherein U is expressed as a result of expanding uncertaintyK is denoted as an inclusion factor and is assigned a value of 2, U P Another result, k, expressed as an extended uncertainty p Represented as inclusion factors determined from the inclusion probability p.
A second aspect of the embodiment of the present invention provides a system for evaluating uncertainty of measurement of an empty box barometer, for implementing the method for evaluating uncertainty of measurement of an empty box barometer according to the first aspect of the embodiment of the present invention, where the system includes:
the temperature coefficient acquisition module is used for placing the empty box barometer in a temperature coefficient verification box, and carrying out temperature coefficient verification on the empty box barometer to obtain a temperature coefficient, wherein the temperature coefficient comprises an influence index of temperature fluctuation of the temperature coefficient verification box on the temperature coefficient and an influence index of temperature uniformity of the temperature coefficient verification box on the temperature coefficient;
The first evaluation module is used for evaluating the uncertainty of the temperature coefficient to obtain an uncertainty evaluation result of the temperature coefficient;
the indication error acquisition module is used for placing the empty box barometer in an indication verification box and carrying out indication verification on the empty box barometer to obtain an indication error, wherein the indication error comprises an influence index of the air leakage rate of the indication verification box on a measurement result and an influence index of the fluctuation of the indication verification box on the measurement result;
and the second evaluation module is used for evaluating the measurement uncertainty of the indication error according to the uncertainty evaluation result of the temperature coefficient so as to obtain an empty box barometer measurement uncertainty evaluation result.
A third aspect of an embodiment of the present invention provides a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the empty box barometer measurement uncertainty assessment method provided by the first aspect.
A fourth aspect of an embodiment of the present invention provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the method for evaluating uncertainty of measurement of a barometer with empty boxes provided in the first aspect when the program is executed by the processor.
The method, the system, the storage medium and the equipment for evaluating the measurement uncertainty of the empty box barometer provided by the embodiment of the invention have the following beneficial effects:
the empty box barometer is placed in a temperature coefficient verification box, and temperature coefficient verification is carried out on the empty box barometer to obtain a temperature coefficient, wherein the temperature coefficient comprises an influence index of temperature fluctuation of the temperature coefficient verification box on the temperature coefficient and an influence index of temperature uniformity of the temperature coefficient verification box on the temperature coefficient; the uncertainty evaluation of the temperature coefficient is measured, and an uncertainty evaluation result of the temperature coefficient is obtained; placing the empty box barometer in an indication verification box, and performing indication verification on the empty box barometer to obtain an indication error, wherein the indication error comprises an influence index of the air leakage rate of the indication verification box on a measurement result and an influence index of the fluctuation of the indication verification box on the measurement result; and according to the uncertainty evaluation result of the temperature coefficient, the uncertainty evaluation of measurement is carried out on the indication value error, so that an uncertainty evaluation result of the measurement of the empty box barometer is obtained, and the accuracy of the uncertainty evaluation of the measurement of the empty box barometer is effectively improved.
Drawings
FIG. 1 is a flowchart of an implementation of a method for evaluating uncertainty of measurement of an empty barometer according to a first embodiment of the present invention;
FIG. 2 is a schematic diagram of a temperature coefficient calibration box arrangement;
FIG. 3 is a schematic diagram of an arrangement of the calibration box;
FIG. 4 is a block diagram of a system for evaluating uncertainty in measurement of air pressure gauge in empty box according to a third embodiment of the present invention;
fig. 5 is a block diagram of an electronic device according to a fourth embodiment of the present invention.
Detailed Description
In order that the invention may be readily understood, a more complete description of the invention will be rendered by reference to the appended drawings. Several embodiments of the invention are presented in the figures. This invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete.
It will be understood that when an element is referred to as being "mounted" on another element, it can be directly on the other element or intervening elements may also be present. When an element is referred to as being "connected" to another element, it can be directly connected to the other element or intervening elements may also be present. The terms "vertical," "horizontal," "left," "right," and the like are used herein for illustrative purposes only.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used herein in the description of the invention is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. The term "and/or" as used herein includes any and all combinations of one or more of the associated listed items.
Example 1
Referring to fig. 1, fig. 1 shows a method for evaluating uncertainty of measurement of an empty box barometer according to a first embodiment of the present invention, by introducing performance test results of a temperature coefficient calibration box and an indication calibration box, the empty box barometer measurement uncertainty evaluation is completed, wherein the temperature coefficient calibration box and the indication calibration box are both provided with barometer, at least three temperature probes arranged in a step-like manner are installed in the temperature coefficient calibration box, and a temperature probe located at a center point is installed in the indication calibration box. When the empty box barometer temperature coefficient is detected, the temperature measuring head is electrically connected with a sensor outside the temperature coefficient detecting box so as to collect temperatures of different points in the temperature coefficient detecting box, please refer to fig. 2, which is a schematic diagram of the arrangement structure of the temperature coefficient detecting box, specifically, 3 test points are arranged in the temperature coefficient detecting box, wherein the point B is located at the center of the whole box body, and the point A and the point C are respectively arranged at the lower left corner and the upper right corner. Fig. 3 is a schematic diagram of an arrangement structure of the indication verification box, specifically, only the center point in the indication verification box is provided with a test point, and the method specifically includes steps S01 to S04.
And S01, placing the empty box barometer in a temperature coefficient verification box, and carrying out temperature coefficient verification on the empty box barometer to obtain a temperature coefficient, wherein the temperature coefficient comprises an influence index of temperature fluctuation of the temperature coefficient verification box on the temperature coefficient and an influence index of temperature uniformity of the temperature coefficient verification box on the temperature coefficient.
In specific implementation, a high temperature point (25-35 ℃) is selected firstly by adjusting a temperature switch of the temperature coefficient verification box, the temperature of the central point reaches a set value, then the temperature is stabilized for 10 minutes, and then the temperature value of each measuring point is read every 5 minutes, and the total reading is 12 times. And then setting the verification box at a low temperature point (0-5 ℃), and respectively reading the temperature values of all the measurement points for 12 times by the same operation.
According to the data acquired by the operation, calculating the temperature fluctuation of the temperature coefficient verification box and the temperature uniformity of the temperature coefficient verification box, wherein the temperature fluctuation of the temperature coefficient verification box is expressed as follows:
；
wherein,t bmax maximum temperature value expressed as 12 measurements at point B, unit: at the temperature of the mixture,t bmin the lowest temperature value expressed as 12 measurements at point B, in units: DEG C.
The expression of the temperature uniformity of the temperature coefficient verification box is as follows:
；
Wherein,t imax expressed as the highest temperature value of A, B, C at the ith measurement, in units of: at the temperature of the mixture,t imin expressed as the lowest temperature value in the ith measurement for the A, B, C three measurement points, in units: DEG C.
In addition, the expression of the measurement model of the temperature coefficient is:
；
wherein,K t expressed as empty box barometer temperature coefficient, units: hPa/. Degree.C,p s1 the average indication value expressed as n corrected barometric pressure standards at high temperature point, unit: the number of hPa,p s2 expressed as average indication of corrected n barometric pressure standards at low temperature point, in: the number of hPa,p b1 expressed as average indication value of the empty box barometer n times at high temperature point, unit: the number of hPa,p b2 expressed as average indication of n empty box barometers at low temperature point, in units: the number of hPa,t 1 expressed as average indication of the temperature at the hot spot n times in units of: at the temperature of the mixture,t 2 expressed as average indication of the low temperature point n thermometer in units: at the temperature of the mixture,q 1 the effect of temperature fluctuation in the temperature coefficient verification box on the temperature coefficient is expressed as the unit: hPa/. Degree.C,q 2 the effect of temperature uniformity in the temperature coefficient verification box on the temperature coefficient is expressed as the unit: hPa/. Degree.C.
And step S02, carrying out measurement uncertainty assessment on the temperature coefficient to obtain an uncertainty assessment result of the temperature coefficient.
In this embodiment, the step of performing the evaluation of the measurement uncertainty includes:
calculating sensitivity coefficients corresponding to the input quantities in the measurement model and standard uncertainty of the input quantities according to the measurement model;
determining each standard uncertainty component according to the sensitivity coefficient and the standard uncertainty corresponding to each input quantity;
calculating the uncertainty of the synthetic standard according to the uncertainty components of each standard;
and calculating the extension uncertainty according to the synthesis standard uncertainty.
When the temperature coefficient is evaluated for measurement uncertainty, the calculated standard uncertainty of each input quantity comprises the standard uncertainty corresponding to an air pressure standard in the temperature coefficient verification box, the standard uncertainty corresponding to an empty box air pressure meter in the temperature coefficient verification box, the standard uncertainty corresponding to the thermometer in the temperature coefficient verification box, the standard uncertainty corresponding to the temperature fluctuation in the temperature coefficient verification box and the standard uncertainty corresponding to the temperature uniformity in the temperature coefficient verification box.
The measurement of repeatability is generally classified into a class a standard uncertainty and a class B standard uncertainty, wherein the measurement of repeatability is classified into a class a uncertainty and the measurement of non-repeatability is classified into a class B standard uncertainty. When evaluating the standard uncertainty, it is first judged whether it satisfies the class a standard uncertainty or the class B standard uncertainty, and then evaluated according to its specific characteristics.
(1) Class a standard uncertainty
According to the requirement of class A standard uncertainty evaluation, the experimental standard deviation is mainly calculated when independent repeated measurement is carried out, and when the independent repeated measurement is carried out, the Bessel formula is generally adopted for calculation. Wherein the Bessel formula is:
。
(2) Class B standard uncertainty
The uncertainty of B class is determined by judging the measured possible value interval x-a, x+a based on the information or experience]Where a denotes the half width of the interval, k is determined from the probability distribution assuming a probability distribution of the measured values (k=when a uniform distribution is satisfied when probabilities of occurrence of the measured possible values at points in the interval are equal)The method comprises the steps of carrying out a first treatment on the surface of the When the measured possible value appears at a larger interval end point value, then when the arcsine distribution is satisfied, k= = ->When the measured likelihood value appears in the intermediate value with a high likelihood, and the normal distribution is satisfied, k=2), the uncertainty of the class B standard can be obtained by the formula u B =a/k。
The expression for calculating the uncertainty of the synthesis standard is:
；
wherein c i The sensitivity coefficient expressed as the i-th input quantity, u (x i ) Standard uncertainty expressed as the i-th input quantity, u c Expressed as synthesis standard uncertainty.
The expression for calculating the extended uncertainty is:
；
Or (b)
；
Wherein U is expressed as a result of the extended uncertainty, k is expressed as an inclusion factor, and a value of 2, U is assigned P Another result, k, expressed as an extended uncertainty p Represented as inclusion factors determined from the inclusion probability p.
From the above description, it can be understood that, in order to obtain the uncertainty evaluation result of the temperature coefficient, first, according to the measurement model, the sensitivity coefficient corresponding to each input quantity in the measurement model and the standard uncertainty of each input quantity are calculated, and the sensitivity coefficient c corresponding to each input quantity in the measurement model is calculated i The following are provided:
；
；
；
；
；
；
；
。
the standard uncertainty of each input is calculated as follows:
1) Standard uncertainty corresponding to air pressure standard in temperature coefficient verification boxu(p s )
The corresponding measurement inaccuracy of the air pressure standard device is mainly introduced, the measurement inaccuracy belongs to the class B standard uncertainty, the measurement uncertainty corresponding to the measured air pressure value of the high-temperature point is obtained by inquiring the traceability certificate (verification or calibration certificate of the upper unit) of the air pressure standard deviceU(p s1 ) And thus it corresponds tou(p s1 )=U(p s1 )/k。
Similarly, the standard uncertainty corresponding to the air pressure measured by the air pressure standard device corresponding to the low temperature pointu(p s2 )=U(p s2 )/k。
2) Standard uncertainty corresponding to empty box barometer in temperature coefficient verification box u(p b )
Repeated measurements, mainly derived from empty-box barometers, belonging to the class A standard uncertainty, by measuring repeatedlynThe numerical value of the secondary empty box barometer is substituted into a Bessel formula to calculate, and the method can obtain when the temperature is highu(p b1 ) Is a value of (2).
Similarly, the standard uncertainty corresponding to the air pressure measured by the empty box air pressure meter corresponding to the low temperature point can be obtained by the same methodu(p b2 )。
3）Standard uncertainty corresponding to thermometer in temperature coefficient verification boxu(t)
Mainly because of inaccurate measurement of the thermometer, the method belongs to class B standard uncertainty, and the extended uncertainty corresponding to the measured temperature value is obtained by inquiring the tracing certificate (verification or calibration certificate of the upper unit) of the thermometerU(t 1 ) And thus it corresponds tou(t 1 )=U(t 1 )/k。
Similarly, the standard uncertainty corresponding to the temperature value measured by the thermometer corresponding to the low temperature pointu(t 2 )=U(t 2 )/k。
4) Standard uncertainty corresponding to temperature fluctuation in temperature coefficient verification boxu(q 1 ):
First, it derives from the delta of the volatility test results of the temperature coefficient verification box while the temperature coefficient verification process is being performedt b (at high temperature)t b1 At low temperature oft b2 ) At high temperature, the corresponding high temperaturet 1 The fluctuation interval of (1) is [t 1 -|△t b1 | ,t 1 +|△t b1 |]In the same way, low temperaturet 2 The fluctuation interval of (1) is [ t 2 -|△t b2 |,t 2 +|△t b2 |]。
Then, upper and lower limit endpoints of the temperature coefficient measurement caused by the fluctuation of the temperature coefficient verification box are known according to a temperature coefficient measurement model: the upper limit endpoint is:
；
the main reason is that when measuring at high temperaturet 1 Lower endpoint in the fluctuation intervalt 1 -|△t b1 I, when measured at low temperaturet 2 And is at the upper end pointt 2 +|△t b2 When I, its corresponding temperature coefficient valueK t Maximum (i.e. upper end point). Similarly, when the temperature is measuredt 1 Upper endpoint in the fluctuation intervalt 1 +|△t b1 I, when measured at low temperaturet 2 And is at the lower end pointt 2 -|△t b2 When I, the measured temperature coefficient valueK t Minimum (i.e., lower endpoint). From the temperature coefficient measurement model, it is known that: the lower limit endpoint is:
；
finally, according to the above conclusion, delta is due to the volatility test resultt b The highest probability of the upper and lower limit endpoints of the interval corresponding to the corresponding temperature measured value is obtained, and therefore, the temperature coefficient is obtainedK t Also falls within the interval [K t2 ,K t1 ]The probability of the upper and lower end points is highest, thus satisfying the arcsine distribution. Corresponding inclusion factork= And the corresponding interval half widtha 1 =(K t1 -K t2 ) And/2, therefore,u(q 1 )=a 1 //>。
5) Standard uncertainty corresponding to temperature uniformity in temperature coefficient verification boxu(q 2 )
First, it derives from the delta of the uniformity test results of the temperature coefficient assay chamber while the temperature coefficient assay process is being performed t u (at high temperature)t u1 At low temperature oft u2 ) At high temperature, the corresponding high temperaturet 1 The uniformity interval of (1) is [t 1 -△t u1 ,t 1 +△t u1 ]In the same way, low temperaturet 2 The uniformity interval of (1) is [t 2 -△t u2 ,t 2 +△t u2 ]。
Then, upper and lower limit endpoints of temperature coefficient measurement caused by uniformity of the temperature coefficient verification box are similarly known according to a temperature coefficient measurement model: the upper limit endpoint is:
；
the lower limit endpoint is:
；
finally, according to the above conclusion, delta is due to the uniformity test resultt u The characteristics of the corresponding temperature measurement value are known to be equal in occurrence probability of each point in the interval corresponding to the corresponding temperature measurement value, so that the temperature coefficient is obtainedK t Also falls within the interval [K t4 ,K t3 ]The probability of each value in the inner is equal, and the uniform distribution is satisfied. Corresponding inclusion factork= And the corresponding interval half widtha 2 =(K t3 -K t4 ) And/2, therefore,u(q 2 )=a 2 //>。
the calculated synthesis standard uncertainty is as follows:
the uncertainty component of each standard isc i |u(x i ) Substituting the uncertainty of the synthesis standard into a calculation formula to obtain:
；
it is arranged as follows:
。
finally calculate the expansion uncertaintyDegree ofUOr (b)U p In whichU p The result of (2) is typically selected to have a inclusion probability of 95% and a corresponding inclusion factork 95 =1.96, thus by calculation, an expansion uncertainty is derivedUOr (b)U 95 The result of (2) is:
U=2×u c ；
U 95 =1.96×u c 。
and S03, placing the empty box barometer in an indication verification box, and carrying out indication verification on the empty box barometer to obtain an indication error, wherein the indication error comprises an influence index of the air leakage rate of the indication verification box on a measurement result and an influence index of the fluctuation of the indication verification box on the measurement result.
When the method is implemented, firstly, the air pressure in the indication verification box is adjusted to a selected test point through an air pressure controller; then, at the time of stabilizing for 10 th minute, the air pressure value P of the air pressure standard is read A . Finally, after 10 minutes, the air pressure value P of the air pressure standard is read again B 。
And calculating the air leakage rate of the indication verification box to be the corresponding air leakage per minute according to the following formula:
；
wherein P is B Air pressure value expressed as air pressure standard after 10 minutes, unit: hPa, P A The air pressure value of the air pressure standard device in the initial state is expressed as the unit: hPa.
Then, after the air pressure in the indication verification box is regulated to a selected test point by an air pressure controller, the air pressure in the box is dynamically stabilized in real time, and after the air pressure in the box is stabilized, the air pressure standard value is read once every 10 seconds, and the air pressure standard value is read altogethernAnd twice.
The volatility of the presentation value calibration box was calculated according to the following formula:
；
wherein,P imax represented asnMaximum barometric pressure value of barometric pressure standard during measurement, unit: the number of hPa,P imin represented asnMinimum barometric pressure value of barometric pressure gauge during measurement, unit: hPa.
In addition, the expression of the measurement model of the indication error is:
；
wherein DeltaP is the air pressure error value of the air pressure gauge of the empty box, P b Represented as an indication of the empty box barometer,K t expressed as an empty box barometer temperature coefficient, T is expressed as ambient temperature,P s indicated as an indication of the corrected barometric pressure standard,Q 1 the effect of the air leakage rate of the indicating value verification box on the measurement result is expressed as the following units: the number of hPa,Q 2 expressed as the impact of the fluctuation of the indicator verification box on the measurement result, in units: hPa.
And step S04, according to the uncertainty evaluation result of the temperature coefficient, evaluating the measurement uncertainty of the indication error to obtain an empty box barometer measurement uncertainty evaluation result.
And when the indication error is evaluated, the calculated standard uncertainty of each input quantity comprises the standard uncertainty corresponding to an air pressure standard in the indication verification box, the standard uncertainty corresponding to an empty box air pressure meter in the indication verification box, the standard uncertainty corresponding to a temperature coefficient, the standard uncertainty corresponding to a thermometer at the central point in the indication verification box, the standard uncertainty corresponding to the air leakage rate of the indication verification box and the standard uncertainty corresponding to the fluctuation of the indication verification box.
In order to obtain the uncertainty evaluation result of the indication error, first, the sensitivity coefficient corresponding to each input in the measurement model and the standard uncertainty of each input are calculated according to the measurement model, and the sensitivity coefficient c corresponding to each input in the measurement model is calculated i The following are provided:
；
；
；
；
；
。
the standard uncertainty of each input is calculated as follows:
1) Standard uncertainty corresponding to air pressure standard in indication verification boxu(P s )
Mainly because the measurement inaccuracy of the air pressure standard device is introduced, the air pressure standard device belongs to the class B standard uncertainty, and the expansion uncertainty corresponding to the measured air pressure value is obtained by inquiring the traceability certificate (verification or calibration certificate of the upper unit) of the air pressure standard deviceU(P s ) And thus it corresponds tou(P s )=U(P s )/k。
2) Standard uncertainty corresponding to empty box barometer in indication verification boxu(P b )
Repeated measurements, mainly derived from empty-box barometers, belonging to the class A standard uncertainty, by measuring repeatedlynSubstituting the value of the secondary empty box barometer into the Bessel formula to calculate to obtainu(P b ) Is a value of (2).
3) Temperature coefficientK t Corresponding standard uncertaintyu(K t )
According to the evaluation result of the corresponding expansion uncertainty during temperature coefficient detection, the corresponding standard uncertainty is calculatedu(K t ). When the expansion uncertainty isUWhen corresponding tou(K t )=U2; when the expansion uncertainty isU 95 When corresponding tou(K t )=U 95 /1.96。
4) Standard uncertainty corresponding to the thermometer at the center point in the indication verification box, namely the temperature value of the standard thermometerTStandard uncertainty corresponding to measurement inaccuracy u(T)
Mainly because of inaccurate measurement of the thermometer, the method belongs to class B standard uncertainty, and the extended uncertainty corresponding to the measured temperature value is obtained by inquiring the tracing certificate (verification or calibration certificate of the upper unit) of the thermometerU(T) And thus it corresponds tou(T)=U(T)/k。
5) Standard uncertainty corresponding to air leakage rate of indication verification boxu(Q 1 )
First, it derives from the delta of the leak rate test results of the prover verification box at the same time as the prover verification processP h (corresponding minute ventilation), therefore, at the firstnWhen the measurement is performed in minutes, the corresponding air leakage interval is [ [P s -n△P h ,P s +n△P h ]。
Then, due to the leak rate test result deltaP h The probability of each value being equal in the interval corresponding to the corresponding temperature measurement value is known, and therefore the uniform distribution is satisfied.
Finally, the corresponding inclusion factork= And the corresponding half width of the interval isn△P b Therefore, the method can be used for manufacturing the optical fiber,u(Q 1 )=n△P b //>/>
6) Standard uncertainty corresponding to fluctuation of indication verification boxu(Q 2 )
First, it derives from the delta of the fluctuation test results of the prover verification box at the same time as the prover verification processP b The corresponding fluctuation interval is [ [P s -|△P b |,P s +|△P b |]。
Then, due to the fluctuation test result deltaP b The highest probability of the upper and lower limit points of the interval corresponding to the corresponding temperature measurement value is known, and thus the arcsine distribution is satisfied.
Finally, the corresponding inclusion factork= And the half width of the corresponding interval is deltaP b I, therefore,u(Q 2 )=△P b //>。
the calculated synthesis standard uncertainty is as follows:
the uncertainty component of each standard isc i |u(x i ) Substituting the uncertainty of the synthesis standard into a calculation formula to obtain:
；
it is arranged as follows:
。
finally, the expansion uncertainty is calculatedUOr (b)U p In whichU p The result of (2) is typically selected to have a inclusion probability of 95% and a corresponding inclusion factork 95 =1.96, thus by calculation, an expansion uncertainty is derivedUOr (b)U 95 The result of (2) is:
U=2×u c ；
U 95 =1.96×u c 。
in summary, in the method for evaluating measurement uncertainty of the empty box barometer in the embodiment of the invention, the empty box barometer is placed in a temperature coefficient verification box, and temperature coefficient verification is performed on the empty box barometer to obtain a temperature coefficient, wherein the temperature coefficient comprises an influence index of temperature fluctuation of the temperature coefficient verification box on the temperature coefficient and an influence index of temperature uniformity of the temperature coefficient verification box on the temperature coefficient; the uncertainty evaluation of the temperature coefficient is measured, and an uncertainty evaluation result of the temperature coefficient is obtained; placing the empty box barometer in an indication verification box, and performing indication verification on the empty box barometer to obtain an indication error, wherein the indication error comprises an influence index of the air leakage rate of the indication verification box on a measurement result and an influence index of the fluctuation of the indication verification box on the measurement result; and according to the uncertainty evaluation result of the temperature coefficient, the uncertainty evaluation of measurement is carried out on the indication value error, so that an uncertainty evaluation result of the measurement of the empty box barometer is obtained, and the accuracy of the uncertainty evaluation of the measurement of the empty box barometer is effectively improved.
Example two
In order to verify the scheme in the implementation of the invention, a DMY3 type empty box barometer is selected as a measured object, and the measurement uncertainty is assessed. Firstly, by adopting an improved verification method, temperature coefficient verification is carried out on the sample, and then indication verification is carried out. Then, the improved verification method is implemented, and the measurement uncertainty of the temperature coefficient is assessed first, and then the measurement uncertainty of the indication error is assessed. And finally, the indication error is used for evaluating the uncertainty of measurement, namely the uncertainty result of the final empty box barometer.
During detection, 745A type digital barometer is selected as an air pressure standard (both temperature coefficient detection and indication detection air pressure standard are adopted), and the standard is obtained by inquiring the tracing certificate (the detection or calibration certificate of the upper unit) of the air pressure standard, and the corresponding expansion in the measurement range is obtainedDegree of uncertainty of exhibitionU=0.02hPa(k=2). The standard thermometer adopts an RCY-1G platinum resistance temperature sensor, and is known by inquiring a tracing certificate (verification or calibration certificate of an upper unit) of the standard thermometer, and the corresponding expansion uncertainty in the measurement range of the standard thermometer U=0.02℃(k=2)。
1) For efficient acquisition of the fluctuation and uniformity of a temperature coefficient verification box versus the measured temperature coefficientK t Is a function of (a) and (b). And carrying out performance test on volatility and uniformity of the temperature coefficient verification box while carrying out temperature coefficient verification on the air box barometer.
The high temperature point is selected to be 30 ℃, after the temperature of the central point reaches a set value, the temperature is stabilized for 10 minutes, and then the temperature values of the air pressure standard, the empty box air pressure meter and each measuring point are read every 5 minutes for 12 times. And then setting the verification box at a low temperature point of 5 ℃, and respectively reading the temperature values of the 12-time barometer, the empty box barometer and each measuring point by the same operation. The corresponding test data are shown in table 1.
Table 1: temperature coefficient verification and temperature coefficient verification box performance test related data
By averaging the corresponding 12 measurementsp s1 =1010.3hPa、p b1 =1010.2hPa、p s2 =1010.2hPa、p b2 =1010.4hPa、t 1 =30.0℃、t 2 Calculation is carried out by substituting the temperature coefficient into a measurement model at the temperature of 5.0 ℃ to obtain the corresponding empty box barometer temperature coefficientK t =0.012hPa/℃。
Substituting 12 times of temperature test data of the corresponding point B into a fluctuation calculation formula of a temperature coefficient verification box to obtain corresponding delta at the high temperature point of 30 DEG Ct b1 = ±0.2 ℃; substituting corresponding A, B, C three-point measurement data into a uniformity calculation formula of the temperature coefficient verification box to obtain delta t u1 =0.24 ℃. Similarly, the pair was calculated in the same manner at a low temperature of 5 DEG CFluctuation of the corresponding temperature coefficient verification box is deltat b2 = ±0.1 ℃ and uniformity deltat u2 =0.21℃。
2) In order to effectively obtain the influence of the air leakage rate and the fluctuation of the indication verification box on the final indication error. And performing performance test on the air leakage rate and uniformity of the indication value verification box while performing indication value verification on the air pressure gauge of the empty box.
Firstly, adjusting the air pressure in the indication verification box to a selected test point through an air pressure controller; then, at the time of stabilization for 10 minutes, 1 air pressure standard value (1 st time is recorded as air pressure value) was read every 5 secondsP A ) Total 12 sets (at this time, just the end of 1 min) of 1 empty box barometer value and 1 standard thermometer value were read, and the corresponding test data are shown in table 2; finally, after 10 minutes, the air pressure value of the air pressure standard device is read againP B Air pressure valueP B =1010.21hPa。
Table 2: indicating value verification and indicating value verification box performance test related data
By averaging the corresponding 12 measurementsP s =1010.20hPa、P b =1010.1hPa、T=20 ℃ and temperature coefficientK t =0.012 hPa/. Degree.C, substituting into the measurement model of the indication verification to calculate to obtain the corresponding air pressure error delta of the empty box barometerP=0.14hPa。
The 1 st air pressure standard value P A After 10 min, the barometric pressure standard value is = 1010.18hPaP B The = 1010.21hPa is substituted into the air leakage rate calculation formula of the indication verification box to obtain the corresponding air leakage rate per minute as deltaP h =0.003hPa/min。
Substituting the 12-time air pressure standard value into the fluctuation calculation formula of the indication verification box to obtain the corresponding fluctuation of the indication verification boxP b =±0.02hPa。
3) Measurement uncertainty assessment of temperature coefficient
The expression of the measurement model of the temperature coefficient is:
；
wherein,K t represented as an empty box barometer temperature coefficient,p s1 represented as the average indication of the corrected n times barometric pressure etalon at the high temperature point,p s2 represented as the average indication of the corrected n times barometric pressure etalon at the low temperature point,p b1 expressed as the average indication of the n times empty box barometer at the high temperature point,p b2 represented as the average indication of the n empty box barometer at the low temperature point,t 1 represented as the average indication of the n times thermometer at the hot spot,t 2 represented as the average indication of the low temperature point n times thermometer,q 1 expressed as the effect of temperature fluctuations in the temperature coefficient assay chamber on the temperature coefficient,q 2 the effect of temperature uniformity in the temperature coefficient assay chamber on the temperature coefficient is expressed.
3.1 Calculating sensitivity coefficients corresponding to each input quantity in the measurement modelc i The following are provided:
；
；
；
；
；
；
；
。
3.2 Standard uncertainty of calculating each input quantity u(x i )
(1) Standard uncertainty corresponding to air pressure standardu(p s )
The corresponding measurement inaccuracy of the air pressure standard device is mainly introduced, the measurement inaccuracy belongs to the class B standard uncertainty, the measurement uncertainty is obtained by inquiring the traceability certificate (verification or calibration certificate of the upper unit) of the air pressure standard device, and the corresponding expansion uncertainty of the measured air pressure value isU=0.02hPa(k=2), and thus it corresponds tou(p s1 ) =0.02 hPa/2=0.01 hPa. Similarly, the standard uncertainty corresponding to the air pressure measured by the air pressure standard device corresponding to the low temperature pointu(p s2 )=0.01hPa。
(2) Standard uncertainty corresponding to empty box barometeru(p b )
Repeated measurements, mainly derived from empty-box barometers, belonging to the class A standard uncertainty, by measuring repeatedly12Substituting the value of the secondary empty box barometer into the Bessel formula to calculate, and obtaining when the temperature is highu(p b1 ) =0.09 hPa; similarly, the standard uncertainty corresponding to the air pressure measured by the empty box air pressure meter corresponding to the low temperature point can be obtained by the same methodu(p b2 )=0.07hPa。
(3) Standard uncertainty corresponding to thermometer in temperature coefficient verification boxu(t)
Mainly because of inaccurate measurement of the thermometer, the method belongs to class B standard uncertainty, and the extended uncertainty corresponding to the measured temperature value is obtained by inquiring the tracing certificate (verification or calibration certificate of the upper unit) of the thermometer U=0.02℃(k=2), and thus it corresponds tou(t 1 ) =0.02 ℃/2=0.01 ℃. Similarly, the standard uncertainty corresponding to the temperature value measured by the thermometer corresponding to the low temperature pointu(t 2 )=0.01℃。
(4) Standard uncertainty corresponding to temperature fluctuation in temperature coefficient verification boxu(q 1 )
First, it derives from the delta of the volatility test results of the temperature coefficient verification box while the temperature coefficient verification process is being performedt b (at high temperature)t b1 = ±0.2 ℃ and delta at low temperaturet b2 = ±0.1 ℃) at high temperatures, the corresponding high temperaturest 1 The fluctuation interval of =30deg.C is [ (30-0.2) °C, (30+0.2) °C)]In the same way, low temperaturet 2 The fluctuation interval of =5℃is [ (5-0.1) ℃and (5+0.1) ℃]。
Then, according to a temperature coefficient measurement model, the following steps are obtained:
the upper limit endpoint is:
；
the lower limit endpoint is:
。
temperature coefficientK t In interval [K t2 ,K t1 ]Satisfy an arcsine distribution, corresponding to the inclusion factork= And the corresponding interval half widtha 1 =(K t1 -K t2 ) /2= 0.000155hPa/°c, and therefore,u(q 1 )=a 1 //>=0.00011hPa/℃。
(5) standard uncertainty corresponding to temperature uniformity in temperature coefficient verification boxu(q 2 )
First, it comes fromUniformity test results delta from a temperature coefficient assay chamber while simultaneously performing a temperature coefficient assay processt u (at high temperature)t u1 =0.24 ℃, delta at low temperaturet u2 =0.21 ℃), at high temperatures, its corresponding high temperature t 1 The uniformity interval of =30deg.C is [ (30-0.24) °C, (30+0.24) °C)]In the same way, low temperaturet 2 The uniformity interval of =5℃is [ (5-0.21) ℃and (5+0.21) ℃]。
Then, according to a temperature coefficient measurement model, the following steps are obtained: the upper limit endpoint is
；
The lower limit endpoint is
；
Temperature coefficientK t In interval [K t4 ,K t3 ]Satisfies a uniform distribution, corresponding inclusion factorsk= And the corresponding interval half widtha 2 =(K t3 -K t4 ) /2= 0.000215hPa/°c, and therefore,u(q 2 )=a 2 //>=0.00012hPa/℃
3.3 Calculating synthetic standard uncertainty
The uncertainty component of each standard isc i |u(x i ) Substituting the uncertainty calculation formula of the synthesis standard to obtain:
；
3.4 Calculating an extended uncertainty
Deriving expansion uncertaintyUOr (b)U 95 The result of (2) is:
U=2×u c =0.009hPa/℃
U 95 =1.96×u c =0.009hPa/℃。
4) Measurement uncertainty assessment of indication error
The expression of the measurement model of the indication error is:
；
wherein DeltaP is the air pressure error value of the air pressure gauge of the empty box,P b represented as an indication of the empty box barometer,K t expressed as an empty box barometer temperature coefficient, T is expressed as ambient temperature,P s indicated as an indication of the corrected barometric pressure standard,Q 1 expressed as the effect of the indicator verification box leak rate on the measurement results,Q 2 expressed as the effect of the fluctuation of the indicator verification box on the measurement result.
4.1 Calculating sensitivity coefficients corresponding to each input quantity in the measurement modelc i The following are provided:
；
；
；
；
；
。
4.2 Standard uncertainty of calculating each input quantityu(x i )
(1) Standard uncertainty corresponding to air pressure value measured by air pressure standard deviceu(P s )
Mainly because the measurement inaccuracy of the air pressure standard device is introduced, the air pressure standard device belongs to the class B standard uncertainty, and the air pressure standard device is known by inquiring the traceability certificate (verification or calibration certificate of a superior unit) of the air pressure standard device, and the expansion uncertainty corresponding to the measured air pressure value is thatU=0.02hPa(k=2), and thus it corresponds tou(P s )=0.02hPa/2=0.01hPa。
(2) Standard uncertainty corresponding to air pressure measured by empty box air pressure meteru(P b )
The repeatability measurement mainly from the empty box barometer belongs to the A standard uncertainty, and the method can obtain the empty box barometer by substituting the numerical value of the empty box barometer representation value verification of the repeated measurement for 12 times into a Bessel formula for calculationu(P b )=0.07hPa。
(3) Temperature coefficientK t Standard uncertainty corresponding to measurement inaccuracy of (2)u(K t )
According to the evaluation result of the corresponding expansion uncertainty during temperature coefficient detection, the corresponding standard uncertainty is calculatedu(K t ). Corresponding tou(K t )=0.0046hPa/℃。
(4) Standard thermometer temperature valueTStandard uncertainty corresponding to measurement inaccuracyu(T)
Mainly because of inaccurate measurement of the thermometer, the method belongs to class B standard uncertainty, and the extended uncertainty corresponding to the measured temperature value is obtained by inquiring the tracing certificate (verification or calibration certificate of the upper unit) of the thermometer U=0.02℃(k=2), and thus it corresponds tou(T)=0.02℃/2=0.01℃。
(5) Standard uncertainty corresponding to air leakage rate of indication verification boxu(Q 1 )
First, it comes from the calibration process of the indication value and the indication valueAir leakage rate test result delta of verification boxP h =0.003 hPa/min, therefore, at the time of measurement at the 1 st minute, the barometerP s 1010.20hPa, the corresponding leakage interval is
[P s -n△P h ,P s +n△P h ]=[1010.20hPa-0.003hPa, 1010.20hPa+0.003hPa]。
Then, due to the leak rate test result deltaP h The probability of each value being equal in the interval corresponding to the corresponding temperature measurement value is known, and therefore the uniform distribution is satisfied.
Finally, the corresponding inclusion factork= And the corresponding interval half width is deltaP b Therefore, the method can be used for manufacturing the optical fiber,u(Q 1 )=△P b /=0.001732hPa
(6) standard uncertainty corresponding to fluctuation of indication verification boxu(Q 2 )
First, it derives from the delta of the fluctuation test results of the prover verification box at the same time as the prover verification processP b = ±0.02hPa, the corresponding fluctuation interval is [P s -|△P b |,P s +|△P b |]=[1010.20hPa-0.02hPa, 1010.20hPa+0.02hPa]。
Then, due to the fluctuation test result deltaP b The highest probability of the upper and lower limit points of the interval corresponding to the corresponding temperature measurement value is known, and thus the arcsine distribution is satisfied.
Finally, the corresponding inclusion factork= And the half width of the corresponding interval is deltaP b I, therefore,u(Q 2 )=△P b //>=0.01414hPa。
4.3 Calculating synthetic standard uncertainty
The uncertainty component of each standard is c i |u(x i ) Substituting the uncertainty calculation formula of the synthesis standard to obtain:
4.4 Calculating an extended uncertainty
Deriving expansion uncertaintyUOr (b)U 95 The result of (2) is:
U=2×0.191hPa=0.382hPa
U 95 =1.96×u c =0.363hPa。
through the operation, the content of the whole measurement uncertainty evaluation accords with the actual measurement condition, so that the measurement uncertainty evaluation result of the empty box barometer can be obtained more truly and reliably through the scheme, and further reliable technical guarantee is provided for improving the quality of measurement data of the empty box barometer.
Example III
Referring to fig. 4, fig. 4 is a block diagram of a system 200 for evaluating uncertainty of air-box barometer measurement according to a third embodiment of the present invention, where the system 200 is configured to implement the method for evaluating uncertainty of air-box barometer measurement, and includes: a temperature coefficient acquisition module 21, a first evaluation module 22, an indication error acquisition module 23 and a second evaluation module 24, wherein:
the temperature coefficient obtaining module 21 is configured to place the empty box barometer in a temperature coefficient verification box, perform temperature coefficient verification on the empty box barometer to obtain a temperature coefficient, where the temperature coefficient includes an influence index of temperature coefficient verification box temperature fluctuation on the temperature coefficient and an influence index of temperature coefficient verification box temperature uniformity on the temperature coefficient, and an expression of a measurement model of the temperature coefficient is:
；
Wherein,K t represented as an empty box barometer temperature coefficient,p s1 represented as the average indication of the corrected n times barometric pressure etalon at the high temperature point,p s2 represented as the average indication of the corrected n times barometric pressure etalon at the low temperature point,p b1 expressed as the average indication of the n times empty box barometer at the high temperature point,p b2 represented as the average indication of the n empty box barometer at the low temperature point,t 1 represented as the average indication of the n times thermometer at the hot spot,t 2 represented as the average indication of the low temperature point n times thermometer,q 1 expressed as the effect of temperature fluctuations in the temperature coefficient assay chamber on the temperature coefficient,q 2 the effect of temperature uniformity in the temperature coefficient verification box on the temperature coefficient is expressed;
a first evaluation module 22, configured to perform an uncertainty measurement evaluation on the temperature coefficient, so as to obtain an uncertainty evaluation result of the temperature coefficient;
the indication error obtaining module 23 is configured to place the empty box barometer in an indication verification box, and perform indication verification on the empty box barometer to obtain an indication error, where the indication error includes an influence index of an indication verification box leakage rate on a measurement result and an influence index of an indication verification box volatility on the measurement result, and an expression of a measurement model of the indication error is:
；
Wherein DeltaP is the air pressure error value of the air pressure gauge of the empty box,P b represented as an indication of the empty box barometer,K t expressed as an empty box barometer temperature coefficient, T is expressed as ambient temperature,P s indicated as an indication of the corrected barometric pressure standard,Q 1 expressed as the effect of the indicator verification box leak rate on the measurement results,Q 2 expressed as a pair of measurements of the fluctuation of an indication verification boxInfluence of the results;
and the second evaluation module 24 is used for evaluating the measurement uncertainty of the indication error according to the uncertainty evaluation result of the temperature coefficient so as to obtain the measurement uncertainty evaluation result of the empty box barometer.
Further, in some alternative embodiments of the present invention, the step of performing the measurement uncertainty assessment in the first assessment module 22 and the second assessment module 24 includes:
calculating sensitivity coefficients corresponding to the input quantities in the measurement model and standard uncertainty of the input quantities according to the measurement model;
determining each standard uncertainty component according to the sensitivity coefficient and the standard uncertainty corresponding to each input quantity;
according to each standard uncertainty component, calculating the uncertainty of the synthetic standard, wherein the expression for calculating the uncertainty of the synthetic standard is as follows:
；
Wherein c i The sensitivity coefficient expressed as the i-th input quantity, u (x i ) Standard uncertainty expressed as the i-th input quantity, u c Expressed as synthesis standard uncertainty;
according to the synthesis standard uncertainty, calculating an extension uncertainty, wherein the expression for calculating the extension uncertainty is as follows:
；
or (b)
；
Wherein U is expressed as a result of the extended uncertainty, k is expressed as an inclusion factor, and a value of 2, U is assigned P Another result, k, expressed as an extended uncertainty p Represented as inclusion factors determined from the inclusion probability p.
When the temperature coefficient is evaluated, the calculated standard uncertainty of each input quantity comprises the standard uncertainty corresponding to an air pressure standard in the temperature coefficient verification box, the standard uncertainty corresponding to an empty box air pressure meter in the temperature coefficient verification box, the standard uncertainty corresponding to the thermometer in the temperature coefficient verification box, the standard uncertainty corresponding to the temperature fluctuation in the temperature coefficient verification box and the standard uncertainty corresponding to the temperature uniformity in the temperature coefficient verification box;
and when the indication error is evaluated, the calculated standard uncertainty of each input quantity comprises the standard uncertainty corresponding to an air pressure standard in the indication verification box, the standard uncertainty corresponding to an empty box air pressure meter in the indication verification box, the standard uncertainty corresponding to a temperature coefficient, the standard uncertainty corresponding to a thermometer at the central point in the indication verification box, the standard uncertainty corresponding to the air leakage rate of the indication verification box and the standard uncertainty corresponding to the fluctuation of the indication verification box.
Example IV
In another aspect, referring to fig. 5, an electronic device according to a fourth embodiment of the present invention includes a memory 20, a processor 10, and a computer program 30 stored on the memory and capable of running on the processor, where the processor 10 implements the above-mentioned method for evaluating uncertainty of air gauge measurement when executing the computer program 30.
The processor 10 may be, among other things, a central processing unit (Central Processing Unit, CPU), a controller, a microcontroller, a microprocessor or other data processing chip for running program code or processing data stored in the memory 20, e.g. executing an access restriction program or the like, in some embodiments.
The memory 20 includes at least one type of readable storage medium including flash memory, a hard disk, a multimedia card, a card memory (e.g., SD or DX memory, etc.), a magnetic memory, a magnetic disk, an optical disk, etc. The memory 20 may in some embodiments be an internal storage unit of the electronic device, such as a hard disk of the electronic device. The memory 20 may also be an external storage device of the electronic device in other embodiments, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash Card (Flash Card) or the like. Further, the memory 20 may also include both internal storage units and external storage devices of the electronic device. The memory 20 may be used not only for storing application software of an electronic device and various types of data, but also for temporarily storing data that has been output or is to be output.
It should be noted that the structure shown in fig. 5 does not constitute a limitation of the electronic device, and in other embodiments the electronic device may comprise fewer or more components than shown, or may combine certain components, or may have a different arrangement of components.
The embodiment of the invention also provides a computer readable storage medium, on which a computer program is stored, which when being executed by a processor, implements the above-mentioned empty box barometer measurement uncertainty assessment method.
Those of skill in the art will appreciate that the logic and/or steps represented in the flow diagrams or otherwise described herein, e.g., a ordered listing of executable instructions for implementing logical functions, can be embodied in any computer-readable medium for use by or in connection with an instruction execution system, apparatus, or device, such as a computer-based system, processor-containing system, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions. For the purposes of this description, a "computer-readable medium" can be any means that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device.
More specific examples (a non-exhaustive list) of the computer-readable medium would include the following: an electrical connection (electronic device) having one or more wires, a portable computer diskette (magnetic device), a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber device, and a portable compact disc read-only memory (CDROM). In addition, the computer readable medium may even be paper or other suitable medium on which the program is printed, as the program may be electronically captured, via, for instance, optical scanning of the paper or other medium, then compiled, interpreted or otherwise processed in a suitable manner, if necessary, and then stored in a computer memory.
It is to be understood that portions of the present invention may be implemented in hardware, software, firmware, or a combination thereof. In the above-described embodiments, the various steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, may be implemented using any one or combination of the following techniques, as is well known in the art: discrete logic circuits having logic gates for implementing logic functions on data signals, application specific integrated circuits having suitable combinational logic gates, programmable Gate Arrays (PGAs), field Programmable Gate Arrays (FPGAs), and the like.
In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiments or examples. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
The foregoing examples illustrate only a few embodiments of the invention, which are described in detail and are not to be construed as limiting the scope of the invention. It should be noted that it will be apparent to those skilled in the art that several variations and modifications can be made without departing from the spirit of the invention, which are all within the scope of the invention. Accordingly, the scope of protection of the present invention is to be determined by the appended claims.
Claims (10)
1. The utility model provides a measurement uncertainty evaluation method of empty box barometer, its characterized in that through introducing the performance test result of temperature coefficient examination case and registration verification case to accomplish the measurement uncertainty evaluation of empty box barometer, wherein temperature coefficient examination case and registration verification case all are provided with the barometer, install in the temperature coefficient examination case and take the form of three temperature gauge head of echelonment arrangement, and the registration verification incasement is installed and is located the temperature gauge head of central point department, temperature gauge head and the outside sensor electric connection of temperature coefficient examination case in temperature coefficient examination case are in order to gather the temperature of different positions in the temperature coefficient examination case, the inside temperature gauge head of registration verification case is in electric connection with the outside sensor of registration verification case to gather the temperature of registration verification incasement central point department, the method includes:
Placing the empty box barometer in a temperature coefficient verification box, and carrying out temperature coefficient verification on the empty box barometer to obtain a temperature coefficient, wherein the temperature coefficient comprises an influence index of temperature fluctuation of the temperature coefficient verification box on the temperature coefficient and an influence index of temperature uniformity of the temperature coefficient verification box on the temperature coefficient;
evaluating the uncertainty of the temperature coefficient to obtain an uncertainty evaluation result of the temperature coefficient;
placing the empty box barometer in an indication verification box, and performing indication verification on the empty box barometer to obtain an indication error, wherein the indication error comprises an influence index of the air leakage rate of the indication verification box on a measurement result and an influence index of the fluctuation of the indication verification box on the measurement result;
and according to the uncertainty evaluation result of the temperature coefficient, evaluating the measurement uncertainty of the indication error to obtain an empty box barometer measurement uncertainty evaluation result.
2. The method for evaluating measurement uncertainty of an empty box barometer according to claim 1, wherein in the step of performing temperature coefficient verification on the empty box barometer to obtain a temperature coefficient, an expression of a measurement model of the temperature coefficient is:
；
Wherein,K t represented as an empty box barometer temperature coefficient,p s1 represented as the average indication of the corrected n times barometric pressure etalon at the high temperature point,p s2 represented as the average indication of the corrected n times barometric pressure etalon at the low temperature point,p b1 expressed as the average indication of the n times empty box barometer at the high temperature point,p b2 represented as the average indication of the n empty box barometer at the low temperature point,t 1 represented as the average indication of the n times thermometer at the hot spot,t 2 represented as the average indication of the low temperature point n times thermometer,q 1 expressed as the effect of temperature fluctuations in the temperature coefficient assay chamber on the temperature coefficient,q 2 the effect of temperature uniformity in the temperature coefficient assay chamber on the temperature coefficient is expressed.
3. The method for evaluating measurement uncertainty of empty box barometer according to claim 2, wherein in the step of placing the empty box barometer in an indication verification box and performing indication verification on the empty box barometer to obtain an indication error, an expression of a measurement model of the indication error is:
；
wherein DeltaP is the air pressure error value of the air pressure gauge of the empty box,P b represented as an indication of the empty box barometer,K t expressed as an empty box barometer temperature coefficient, T is expressed as ambient temperature,P s indicated as an indication of the corrected barometric pressure standard, Q 1 The air leakage rate of the indication verification box is expressed as the influence of the air leakage rate of the indication verification box on the measurement result，Q 2 Expressed as the effect of the fluctuation of the indicator verification box on the measurement result.
4. A method of assessing measurement uncertainty of an empty barometer according to claim 3, wherein said step of assessing measurement uncertainty comprises:
calculating sensitivity coefficients corresponding to the input quantities in the measurement model and standard uncertainty of the input quantities according to the measurement model;
determining each standard uncertainty component according to the sensitivity coefficient and the standard uncertainty corresponding to each input quantity;
calculating the uncertainty of the synthetic standard according to the uncertainty components of each standard;
and calculating the extension uncertainty according to the synthesis standard uncertainty.
5. The method for evaluating measurement uncertainty of empty box barometer according to claim 4, wherein when the measurement uncertainty is evaluated on the temperature coefficient, the calculated standard uncertainty of each input quantity comprises a standard uncertainty corresponding to an air pressure standard in the temperature coefficient verification box, a standard uncertainty corresponding to an empty box barometer in the temperature coefficient verification box, a standard uncertainty corresponding to a thermometer in the temperature coefficient verification box, a standard uncertainty corresponding to temperature fluctuation in the temperature coefficient verification box, and a standard uncertainty corresponding to temperature uniformity in the temperature coefficient verification box;
And when the indication error is evaluated, the calculated standard uncertainty of each input quantity comprises the standard uncertainty corresponding to an air pressure standard in the indication verification box, the standard uncertainty corresponding to an empty box air pressure meter in the indication verification box, the standard uncertainty corresponding to a temperature coefficient, the standard uncertainty corresponding to a thermometer at the central point in the indication verification box, the standard uncertainty corresponding to the air leakage rate of the indication verification box and the standard uncertainty corresponding to the fluctuation of the indication verification box.
6. The method according to claim 5, wherein in the step of calculating the uncertainty of the synthetic standard based on each standard uncertainty component, the expression for calculating the uncertainty of the synthetic standard is:
；
wherein c i The sensitivity coefficient expressed as the i-th input quantity, u (x i ) Standard uncertainty expressed as the i-th input quantity, u c Expressed as synthesis standard uncertainty.
7. The method of evaluating uncertainty of air gauge measurement according to claim 6, wherein in the step of calculating an extended uncertainty based on the synthetic standard uncertainty, an expression for calculating an extended uncertainty is:
；
Or (b)
；
Wherein U is expressed as a result of the extended uncertainty, k is expressed as an inclusion factor, and a value of 2, U is assigned P Another result, k, expressed as an extended uncertainty p Represented as inclusion factors determined from the inclusion probability p.
8. A system for evaluating measurement uncertainty of an empty box barometer, for implementing the method for evaluating measurement uncertainty of an empty box barometer of any of claims 1-7, said system comprising:
the temperature coefficient acquisition module is used for placing the empty box barometer in a temperature coefficient verification box, and carrying out temperature coefficient verification on the empty box barometer to obtain a temperature coefficient, wherein the temperature coefficient comprises an influence index of temperature fluctuation of the temperature coefficient verification box on the temperature coefficient and an influence index of temperature uniformity of the temperature coefficient verification box on the temperature coefficient;
the first evaluation module is used for evaluating the uncertainty of the temperature coefficient to obtain an uncertainty evaluation result of the temperature coefficient;
the indication error acquisition module is used for placing the empty box barometer in an indication verification box and carrying out indication verification on the empty box barometer to obtain an indication error, wherein the indication error comprises an influence index of the air leakage rate of the indication verification box on a measurement result and an influence index of the fluctuation of the indication verification box on the measurement result;
And the second evaluation module is used for evaluating the measurement uncertainty of the indication error according to the uncertainty evaluation result of the temperature coefficient so as to obtain an empty box barometer measurement uncertainty evaluation result.
9. A computer readable storage medium having stored thereon a computer program, which when executed by a processor implements the empty box barometer measurement uncertainty assessment method according to any one of claims 1-7.
10. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, the processor implementing the empty box barometer measurement uncertainty assessment method of any of claims 1-7 when the program is executed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410145927.8A CN117686140B (en) | 2024-02-02 | 2024-02-02 | Method and system for evaluating measurement uncertainty of empty box barometer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410145927.8A CN117686140B (en) | 2024-02-02 | 2024-02-02 | Method and system for evaluating measurement uncertainty of empty box barometer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117686140A true CN117686140A (en) | 2024-03-12 |
| CN117686140B CN117686140B (en) | 2024-04-05 |
Family
ID=90128493
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410145927.8A Active CN117686140B (en) | 2024-02-02 | 2024-02-02 | Method and system for evaluating measurement uncertainty of empty box barometer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117686140B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117949692A (en) * | 2024-03-27 | 2024-04-30 | 江西省气象探测中心 | Evaluation method and device for measurement uncertainty of photoelectric cup type anemometer |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106525334A (en) * | 2016-11-11 | 2017-03-22 | 山东省计量科学研究院 | Automatic calibrating device for aneroid barometers |
| KR20180094190A (en) * | 2017-02-14 | 2018-08-23 | 연세대학교 원주산학협력단 | Pressure sensor calibration management system |
| CN110231122A (en) * | 2018-03-05 | 2019-09-13 | 霍尼韦尔国际公司 | Device and method for verifying the operation of air data probe |
| CN114624789A (en) * | 2022-03-18 | 2022-06-14 | 陕西省地震局 | Measurement model and uncertainty evaluation method for flux gate theodolite instrument difference |
| CN117008502A (en) * | 2023-02-24 | 2023-11-07 | 浙江威星智能仪表股份有限公司 | MCU measuring circuit and device calibrated by external reference voltage |
| CN117150386A (en) * | 2023-10-31 | 2023-12-01 | 江西省气象探测中心 | Assessment method and device for measurement uncertainty of humidity sensor based on self-adaption |
-
2024
- 2024-02-02 CN CN202410145927.8A patent/CN117686140B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106525334A (en) * | 2016-11-11 | 2017-03-22 | 山东省计量科学研究院 | Automatic calibrating device for aneroid barometers |
| KR20180094190A (en) * | 2017-02-14 | 2018-08-23 | 연세대학교 원주산학협력단 | Pressure sensor calibration management system |
| CN110231122A (en) * | 2018-03-05 | 2019-09-13 | 霍尼韦尔国际公司 | Device and method for verifying the operation of air data probe |
| CN114624789A (en) * | 2022-03-18 | 2022-06-14 | 陕西省地震局 | Measurement model and uncertainty evaluation method for flux gate theodolite instrument difference |
| CN117008502A (en) * | 2023-02-24 | 2023-11-07 | 浙江威星智能仪表股份有限公司 | MCU measuring circuit and device calibrated by external reference voltage |
| CN117150386A (en) * | 2023-10-31 | 2023-12-01 | 江西省气象探测中心 | Assessment method and device for measurement uncertainty of humidity sensor based on self-adaption |
Non-Patent Citations (4)
| Title |
|---|
| 党选发;牛玉质;宋媛;黄蕊;: "空盒气压表测量结果的不确定度评定", 气象水文海洋仪器, no. 01, 15 March 2010 (2010-03-15) * |
| 李昕娣;陈武框;: "空盒气压表示值误差的不确定度评定实例", 气象水文海洋仪器, no. 03, 15 September 2015 (2015-09-15) * |
| 王国芳;陈明谦;林展祺;: "空盒气压表温度系数的不确定度评定", 计量与测试技术, no. 03, 30 March 2020 (2020-03-30) * |
| 罗家伟;: "空盒气压表示值修正值的不确定度评定", 中国计量, no. 10, 10 October 2009 (2009-10-10) * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117949692A (en) * | 2024-03-27 | 2024-04-30 | 江西省气象探测中心 | Evaluation method and device for measurement uncertainty of photoelectric cup type anemometer |
| CN117949692B (en) * | 2024-03-27 | 2024-05-31 | 江西省气象探测中心 | Evaluation method and device for measurement uncertainty of photoelectric cup type anemometer |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117686140B (en) | 2024-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN117686140B (en) | Method and system for evaluating measurement uncertainty of empty box barometer | |
| US7532992B2 (en) | Measuring apparatuses and methods of using them | |
| CN116609552B (en) | Wind speed measurement uncertainty assessment method, system, storage medium and device | |
| US10359308B2 (en) | Flow meter and a method of calibration | |
| CN110220945B (en) | Full-range temperature compensation method of semiconductor gas sensor | |
| CN103091366B (en) | Dewpoint calibration testing method used under complicated environments | |
| CN117150386B (en) | Assessment method and device for measurement uncertainty of humidity sensor based on self-adaption | |
| CN112461406A (en) | Calibration method based on fiber grating type temperature sensor | |
| CN114812758A (en) | Calibration method and device of liquid level sensor, electronic equipment and storage medium | |
| CN1335932A (en) | Method and apparatus for calibrating analog sensor measurement | |
| CN106370779A (en) | Measurement deviation bridging method based on near disturbance source and environmental change | |
| CN115956208A (en) | Method and system for universal calibration of a device under test | |
| CN115574980B (en) | Temperature calibration method and electronic equipment | |
| CN113358243A (en) | Temperature measuring sheet calibration method | |
| Castrup et al. | Uncertainty analysis for alternative calibration scenarios | |
| CN113655093B (en) | Gas concentration detection method, device, equipment and medium | |
| CN112366148B (en) | Substrate concentration determination method, substrate concentration determination device, computer equipment and readable storage medium | |
| JP2022047339A (en) | Sensor drift evaluation apparatus, sensor drift evaluation method, and sensor drift evaluation program | |
| CN116380339B (en) | Film gauge vacuum gauge calibration method and related equipment | |
| CN115326276B (en) | Temperature compensation and environment vacuum degree calculation method and equipment for Pirani vacuum gauge | |
| Wang et al. | Calibration of psychrometers at spring Singapore | |
| CN111307865B (en) | Method for improving relative humidity measurement precision | |
| CN110455216B (en) | Method for evaluating effectiveness of surface shape measurement method based on cross correlation | |
| Ainsworth | An independent quality check of secondary platinum resistance thermometer calibrations above 0 C | |
| White et al. | Designing Calibrations in a Quality Environment |
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 |