CN116609552B - Wind speed measurement uncertainty assessment method, system, storage medium and device - Google Patents
Wind speed measurement uncertainty assessment method, system, storage medium and device Download PDFInfo
- Publication number
- CN116609552B CN116609552B CN202310879686.5A CN202310879686A CN116609552B CN 116609552 B CN116609552 B CN 116609552B CN 202310879686 A CN202310879686 A CN 202310879686A CN 116609552 B CN116609552 B CN 116609552B
- Authority
- CN
- China
- Prior art keywords
- wind speed
- class
- parameters
- uncertainty
- expressed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P21/00—Testing or calibrating of apparatus or devices covered by the preceding groups
- G01P21/02—Testing or calibrating of apparatus or devices covered by the preceding groups of speedometers
- G01P21/025—Testing or calibrating of apparatus or devices covered by the preceding groups of speedometers for measuring speed of fluids; for measuring speed of bodies relative to fluids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明提供了一种风速测量不确定度评定方法、系统、存储介质及设备，该方法通过确定风速测量不确定度的评定参数，并将评定参数分类，得到A类参数和B类参数，其中，评定参数至少包括风洞的稳定性系数和均匀性系数；在检定点风速下，对风洞的稳定性和均匀性进行测试，并计算它们之间的相关系数值；在检定点风速下，分别对A类参数和B类参数的分量值进行计算；根据A类参数和B类参数的分量值计算结果以及相关系数值，计算风速合成标准不确定度评定结果，后转换为扩展不确定度，并输出，具体的，由于充分考虑风洞流场均匀性和风洞流场稳定性这两个输入量之间的相关性，能够有效提高风速测量结果的可信程度。
The invention provides a wind speed measurement uncertainty assessment method, system, storage medium and equipment. The method determines the assessment parameters of the wind speed measurement uncertainty and classifies the assessment parameters to obtain Class A parameters and Class B parameters, where , the evaluation parameters at least include the stability coefficient and uniformity coefficient of the wind tunnel; at the wind speed at the calibration point, test the stability and uniformity of the wind tunnel, and calculate the correlation coefficient value between them; at the wind speed at the calibration point, The component values of Class A parameters and Class B parameters are calculated respectively; based on the component value calculation results of Class A parameters and Class B parameters and the correlation coefficient values, the wind speed composite standard uncertainty assessment results are calculated, and then converted into expanded uncertainties. , and output. Specifically, due to the full consideration of the correlation between the two input quantities of wind tunnel flow field uniformity and wind tunnel flow field stability, the credibility of the wind speed measurement results can be effectively improved.
Description
Technical Field
The invention belongs to the technical field of wind speed measurement uncertainty evaluation, and particularly relates to a wind speed measurement uncertainty evaluation method, a wind speed measurement uncertainty evaluation system, a storage medium and wind speed measurement uncertainty evaluation equipment.
Background
With the development of society, the accuracy of meteorological observation is receiving more and more attention, wherein the assessment of uncertainty of wind speed measurement is also one of the important points of research.
The existing wind speed measurement uncertainty assessment method generally adopts the method of The guide to the expression of uncertainty in measurement (GUM method) published by the International organization for Association of ISO (1993), namely the measurement uncertainty assessment guide method, and the general flow of the method is as follows: analyzing uncertainty sources and building a measurement model, evaluating standard uncertainty, calculating synthetic standard uncertainty, determining extended uncertainty, and reporting measurement results, wherein uncertainty propagation rate is generally used for calculation in the step of calculating synthetic standard uncertainty, while specific flow steps need to be selected according to whether correlation exists between standard uncertainty components.
It should be noted that, for the first step of the procedure of GUM evaluation measurement uncertainty in the wind speed measurement process, when analyzing uncertainty sources and establishing a measurement model, the uncertainty sources generally comprise factors such as a wind speed standard and matched equipment, wherein the uncertainty sources introduced by the flow field of the wind tunnel in the matched equipment generally comprise introduction of uniformity of the flow field of the wind tunnel and introduction of stability of the flow field of the wind tunnel. Whereas existing wind speed measurement uncertainty assessment processes typically consider the relationship between these two inputs as uncorrelated for ease of analysis and computation. However, according to practical situations, the uniformity of the wind tunnel flow field is often changed, and meanwhile, the uniformity of the wind tunnel flow field is also changed, so that the stability of the wind tunnel flow field is changed, namely, certain correlation exists between the wind tunnel flow field and the wind tunnel flow field. If the relation between the uniformity of the wind tunnel flow field and the stability of the wind tunnel flow field is simply considered according to uncorrelation, the finally obtained result of uncertainty of wind speed measurement is necessarily inaccurate.
Disclosure of Invention
Based on the above, the embodiment of the invention provides a method and a system for evaluating uncertainty of wind speed measurement, which aim to solve the problem that in the prior art, the relation between two input quantities, namely wind tunnel flow field uniformity and wind tunnel flow field stability, is considered according to uncorrelation, so that an obtained result of uncertainty of wind speed measurement is inaccurate.
A first aspect of an embodiment of the present invention provides a wind speed measurement uncertainty assessment method, the method including:
determining evaluation parameters of the uncertainty of wind speed measurement according to actual measurement conditions, and classifying the evaluation parameters to obtain class A parameters for evaluating class A uncertainty components and class B parameters for evaluating class B uncertainty components, wherein the evaluation parameters at least comprise the stability and uniformity of a wind tunnel;
testing the stability and uniformity of the wind tunnel under the verification point wind speed, and calculating a correlation coefficient value according to the data of the stability and uniformity test of the wind tunnel;
under the wind speed of a verification point, calculating component values of the class A parameter and the class B parameter according to a class A uncertainty component evaluation standard and a class B uncertainty component evaluation standard;
Calculating a wind speed synthesis standard uncertainty assessment result according to the component value calculation results of the class A parameters and the class B parameters and the correlation coefficient value;
and converting the wind speed synthesis standard uncertainty evaluation result into an expanded uncertainty, and outputting the expanded uncertainty.
Further, the step of determining the evaluation parameters of the uncertainty of the wind speed measurement according to the actual measurement situation, and classifying the evaluation parameters to obtain the class a parameters for evaluating the class a uncertainty component and the class B parameters for evaluating the class B uncertainty component includes:
acquiring evaluation parameters of uncertainty of wind speed measurement, and sequentially judging whether each evaluation parameter is a parameter introduced by repeated measurement;
if yes, the corresponding assessment parameters are divided into class A parameters;
if not, the corresponding assessment parameters are divided into B-class parameters.
Further, the step of testing the stability and uniformity of the wind tunnel at the verification point wind speed and calculating the correlation coefficient value according to the data of the stability and uniformity test of the wind tunnel comprises the following steps:
the pitot tube is arranged at the center point of the section of the selected wind tunnel test section, and a plurality of wind measuring points are respectively arranged at the upstream and downstream equidistant positions of the pitot tube, wherein the wind measuring points positioned at one side of the pitot tube are uniformly distributed on the same plane;
Controlling the blower motor to rotate at a preset rotating speed, and collecting a first wind speed value through a pitot tube arranged at the center point of the section of the selected wind tunnel test section after the flow speed is stable;
controlling the blower motor to rotate at a preset rotating speed, and after the flow speed is stable, sequentially installing the pitot tubes to all wind measuring points to acquire a second wind speed value;
and calculating the correlation coefficient value according to the first wind speed value and the second wind speed value.
Further, in the step of calculating the correlation coefficient value according to the first wind speed value and the second wind speed value, a calculation formula of the correlation coefficient value is:
wherein ,xd Representing the first wind speed value, y, measured the d-th time d A second wind speed value, n, representing the d-th measurement 1 Expressed as a first measurement totalThe times s (x) and s (y) respectively represent the corresponding test standard deviation,expressed as average of all first wind speeds,Expressed as the average of all second wind speeds,Expressed as correlation coefficient values.
Further, in the step of calculating the component values of the class a parameter and the class B parameter according to the class a uncertainty component evaluation standard and the class B uncertainty component evaluation standard at the verification point wind speed, a component value calculation formula of the class a parameter is:
wherein ,component values expressed as class A parameters, +.>The wind speed measurement error value expressed as the ith test,expressed as the average of wind speed measurement errors, n 2 Expressed as the second total number of measurements, +.>Display values of wind speed sensor expressed as ith test, r,/and so on>，Respectively representing the correction coefficient, the Pitot tube coefficient and the micropressure coefficient of the liquid in the digital micropressure meter,Expressed as the actual measured wind pressure of the digital micro manometer,t、p、u、e w expressed as temperature, atmospheric pressure, humidity and saturated water vapor pressure, respectively.
Further, in the step of calculating the component values of the class a parameter and the class B parameter according to the class a uncertainty component evaluation standard and the class B uncertainty component evaluation standard at the verification point wind speed, the component value calculation formula of the class B parameter is as follows:
wherein ,component values expressed as class B parameters, +.>Half-width of interval expressed as class B parameter, < >>Represented as inclusion factors.
Further, in the step of calculating the wind speed synthesis standard uncertainty evaluation result according to the component value calculation results of the class a parameter and the class B parameter and the correlation coefficient value, a calculation formula of the wind speed synthesis standard uncertainty evaluation result is as follows:
wherein ,expressed as windFast synthesis standard uncertainty assessment results, n 3 Expressed as the total number of B-class parameters, u f Component value expressed as f-th parameter of B-class parameters,>component values expressed as stability, < >>Component values expressed as uniformity.
A second aspect of an embodiment of the present invention provides a wind speed measurement uncertainty assessment system, the system comprising:
the evaluation parameter determining module is used for determining an evaluation parameter of the uncertainty of the wind speed measurement according to the actual measurement condition, classifying the evaluation parameter to obtain a class A parameter for evaluating the component of the uncertainty of the class A and a class B parameter for evaluating the component of the uncertainty of the class B, wherein the evaluation parameter at least comprises the stability and the uniformity of the wind tunnel;
the correlation coefficient value calculation module is used for testing the stability and uniformity of the wind tunnel under the wind speed of the verification point and calculating the correlation coefficient value according to the data of the stability and uniformity test of the wind tunnel;
the component value calculation module is used for calculating component values of the class A parameters and the class B parameters respectively according to the class A uncertainty component evaluation standard and the class B uncertainty component evaluation standard under the wind speed of the verification point;
The assessment result calculation module is used for calculating an uncertainty assessment result of the wind speed synthesis standard according to the component value calculation results of the class A parameters and the class B parameters and the correlation coefficient values;
and the output module is used for converting the wind speed synthesis standard uncertainty evaluation result into an expansion uncertainty and outputting the expansion uncertainty.
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 wind speed measurement uncertainty assessment method provided by the first aspect.
A fourth aspect of an embodiment of the invention provides 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 wind speed measurement uncertainty assessment method provided in the first aspect when executing the program.
The wind speed measurement uncertainty assessment method, the system, the storage medium and the device provided by the embodiment of the invention have the following beneficial effects:
determining an evaluation parameter of the uncertainty of wind speed measurement according to the actual measurement condition, and classifying the evaluation parameter to obtain a class A parameter for evaluating a class A uncertainty component and a class B parameter for evaluating a class B uncertainty component, wherein the evaluation parameter at least comprises a stability coefficient and a uniformity coefficient of a wind tunnel; testing the stability and uniformity of the wind tunnel under the verification point wind speed, and calculating a correlation coefficient value according to the data of the stability and uniformity test of the wind tunnel; under the wind speed of a verification point, calculating component values of the class A parameter and the class B parameter according to the class A uncertainty component evaluation standard and the class B uncertainty component evaluation standard; calculating a wind speed synthesis standard uncertainty assessment result according to the component value calculation results and the related coefficient values of the class A parameters and the class B parameters; the wind speed synthesis standard uncertainty assessment result is converted into the expansion uncertainty and output, and specifically, the reliability of the wind speed measurement result can be effectively improved due to the fact that the correlation between the uniformity of the wind tunnel flow field and the stability of the wind tunnel flow field is fully considered, and the wind speed synthesis standard uncertainty assessment method has important significance in improving the quality of wind speed observation data and revealing weather rules and weather characteristics.
Drawings
FIG. 1 is a flowchart of an implementation of a wind speed measurement uncertainty assessment method provided by a first embodiment of the present invention;
FIG. 2 is a schematic diagram of stability and uniformity control points of a wind tunnel wind measuring section;
FIG. 3 is a block diagram of a wind speed measurement uncertainty evaluation system provided by a second embodiment of the present invention;
fig. 4 is a block diagram of an electronic device according to a third 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 wind speed measurement uncertainty assessment method according to a first embodiment of the present invention, and the method specifically includes steps S01 to S05.
And S01, determining evaluation parameters of the uncertainty of wind speed measurement according to actual measurement conditions, and classifying the evaluation parameters to obtain class A parameters for evaluating class A uncertainty components and class B parameters for evaluating class B uncertainty components, wherein the evaluation parameters at least comprise the stability and the uniformity of the wind tunnel.
Specifically, first, according to the actual measurement situation, an evaluation parameter of uncertainty of wind speed measurement should be determined, and in this embodiment, a wind speed standard device adopted in wind speed measurement includes a standard wind speed measurement device and a wind speed generation and control device. The standard wind speed measuring device consists of a pitot tube, a micro-manometer and an environment temperature-humidity-pressure sensor, and the wind tunnel is a wind speed generating and controlling device. When wind speed is measured, the indicated wind speed value of the measured sensor is compared with the standard wind speed value of the standard wind speed measuring device, and whether the measured sensor is qualified is judged by solving an error mode. The evaluation parameters considering the uncertainty of wind speed measurement at least comprise an indication wind speed value of a measured sensor, a standard wind speed value, pitot tube accuracy, a pitot tube installation angle error, a pitot tube blocking error, a micro differential pressure meter uncertainty, a correction error during computer processing, an ambient temperature measurement error, an ambient humidity measurement error, an ambient air pressure measurement error, an error introduced by the stability of a wind tunnel and an error introduced by uniformity. It should be noted that, according to the actual measurement situation, the target measurement model may be expressed as:
wherein ,expressed as a target measurement model, V M Indicating the wind speed value, V, for the sensor to be measured S Expressed as standard wind speed value, V 1 Expressed as pitot tube accuracy introduced error, V 2 Expressed as error introduced by pitot tube mounting angle deviation, V 3 Expressed as pitot tube obstruction induced error, V 4 Expressed as error, V, introduced by micro-differential pressure gauge uncertainty 5 Expressed as error introduced by ambient temperature measurement, V 6 Expressed as error, V, introduced by ambient humidity measurements 7 Expressed as error, V, introduced by ambient air pressure measurement 8 Expressed as a reduction error, V, in computer processing 9 Expressed as an error introduced by the stability of the wind tunnel, V 10 Represented as errors introduced by the uniformity of the wind tunnel.
Further, judging whether the evaluation parameters are parameters introduced by repeated measurement in sequence, if so, dividing the corresponding evaluation parameters into class A parameters; if not, the corresponding assessment parameters are divided into B-class parameters. Specifically, the standard uncertainty introduced by the repeatability measurement is a class A standard uncertainty, "V" in the target measurement model M -V S "components introduced for repeatability measurements, i.e. class A parameters, while" V "in the target measurement model 1 ~V 10 "is a component introduced by uncertainty of class B criteria, i.e., a class B parameter.
And S02, testing the stability and uniformity of the wind tunnel at the verification point wind speed, and calculating a correlation coefficient value according to the stability and uniformity test data of the wind tunnel.
Specifically, a pitot tube is arranged at the center point of the section of a selected wind tunnel test section, and a plurality of wind measuring points are respectively arranged at the upstream and downstream equidistant positions of the pitot tube, wherein the wind measuring points positioned at one side of the pitot tube are uniformly distributed on the same plane; controlling the blower motor to rotate at a preset rotating speed, and collecting a first wind speed value through a pitot tube arranged at the center point of the section of the selected wind tunnel test section after the flow speed is stable; controlling the blower motor to rotate at a preset rotating speed, and after the flow speed is stable, sequentially installing the pitot tubes to all wind measuring points to acquire a second wind speed value; and calculating a correlation coefficient value according to the first wind speed value and the second wind speed value.
Referring to fig. 2, a schematic diagram of stability and uniformity control points of a wind section of a wind tunnel is shown, in this embodiment, according to the requirements of QX/T84-2007 "meteorological low speed wind tunnel performance test specification", stability and uniformity of an HDF500 loop wind tunnel adopted in an experiment are tested, specifically, a pitot tube is installed according to a selected measuring point position, a motor rotation speed is controlled, after a measured section flow velocity is stable, a wind speed indication value starts to be read, when a first wind speed value is acquired, after the measured section flow velocity is stable, a wind speed is recorded every 5 seconds, and 72 times of data are measured in an accumulated manner. When the second wind speed value is acquired, after the flow speed of the section to be measured is stable, the wind speed indication value is read, the wind speed of one test point is recorded every 5 seconds, then the next point is sequentially moved, and 72 times of data are measured in an accumulated mode.
And finally, respectively bringing the stability data and the uniformity data which are measured for 72 times, namely the first wind speed value and the second wind speed value into calculation formulas of a stability evaluation value, a uniformity evaluation value and a related coefficient value, wherein the calculation formulas of the stability evaluation value are as follows:
wherein ,expressed as stability evaluation value, < >>Expressed as the mth instantaneous flow rate value in a certain test,/->The flow velocity is expressed as an average flow velocity value of one minute in a certain test, and it can be understood that the flow velocity is used for calculating the stability of all the obtained first wind velocity values, specifically, a pitot tube is installed at the center O point of the section of a selected test section, the rotating speed of a motor is controlled, after the flow velocity is stable, wind velocity indication values are displayed every 5 seconds, the number of tests is 1 for 12 times, 6 groups are counted, 72 times of data are accumulated, and in addition, an arithmetic average value is calculated from the stability evaluation values measured for 6 times, and the calculated average value is the stability evaluation value under the measured flow velocity.
The calculation formula of the uniformity evaluation value is as follows:
wherein ,expressed as uniformity evaluation value,/->Expressed as the difference between the first windfinding flow rate and the average cross-sectional flow rate, +.>The average flow velocity value of the measured section is represented as the number of wind measuring points, L is represented as the number of wind measuring points, and it can be understood that the average flow velocity value is used for calculating the uniformity of all the obtained second wind velocity values, specifically, the installation section of the wind tunnel main test section instrument is the section where the O point is located, and the sections at the positions of 200mm above and below the wind tunnel main test section instrument are respectively the section A and the section B. The section A and the section B adopt an equidistant method to arrange wind measuring points (36 points are distributed and controlled on each section and 72 wind speed measuring points are totally) and the calculation formula of the correlation coefficient value is as follows:
wherein ,xd Representing the first wind speed value, y, measured the d-th time d A second wind speed value, n, representing the d-th measurement 1 Expressed as the first total number of measurements, n is known from the above 1 72, s (x) and s (y) respectively represent corresponding test standard deviations,expressed as average of all first wind speeds,Expressed as the average of all second wind speeds,Expressed as correlation coefficient values.
In this embodiment, in order to effectively calculate the correlation coefficient between stability and uniformity, the stability test and the uniformity test are performed synchronously, that is, the stability test data (main pitot tube) is read once every 5 seconds, and the uniformity test data (data of pitot tube a or pitot tube B) is read at the same time, and the test data of the stability and uniformity are shown in table 1, taking the verification point wind speed of 10m/s as an example.
Table 1:10m/s stability and uniformity test data (unit: m/s)
The stability and uniformity data measured in table 1 are respectively brought into the calculation formulas of the stability evaluation value, the uniformity evaluation value and the related coefficient value, and the stability evaluation value, the uniformity evaluation value and the related coefficient value of the wind tunnel are obtained through calculation and are shown in table 2. Taking a 10m/s detection point as an example, the correlation coefficient value r (x, y) = 0.1561 of wind tunnel stability and uniformity.
Table 2: wind tunnel stability average value, uniformity evaluation value and correlation coefficient value
And S03, respectively calculating component values of the class A parameters and the class B parameters according to the class A uncertainty component evaluation standard and the class B uncertainty component evaluation standard at the verification point wind speed.
Specifically, the component value calculation formula of the class a parameter is:
wherein ,component values expressed as class A parameters, +.>Wind speed measurement error value denoted as ith test,/">Expressed as the average of wind speed measurement errors, n 2 Expressed as the second total number of measurements, +.>Display values of wind speed sensor expressed as ith test, r,/and so on>，Respectively representing the correction coefficient, the Pitot tube coefficient and the micropressure coefficient of the liquid in the digital micropressure meter,Expressed as the actual measured wind pressure of the digital micro manometer,t、p、u、e w expressed as temperature, atmospheric pressure, humidity and saturated water vapor pressure, respectively, it being understood that references herein are made tot、p、uI.e. data of ambient temperature, ambient pressure, ambient humidity.
In the embodiment of the invention, each verification point reads 10 times, namely n 2 Taking 10m/s as an example of a calibration point, it should be noted that, in general calibration regulations, standard wind speed values all belong to the calibration point near the calibration point, and specifically, the corresponding original records are shown in table 3. Calculation of u by means of a component value calculation formula with class A parameters A =0.0075m/s。
Table 3:10m/s wind speed original record (unit: m/s)
The component value calculation formula of the B-class parameter is as follows:
wherein ,component values expressed as class B parameters, +.>Represented as B-type ginsengHalf width of interval>Expressed as inclusion factors, it is noted that the interval range values of the measured parameters, e.g., [ w-a, w+a ], are first determined based on the relevant historical information (data) or experience]A is interval half width, the probability distribution of the measured value is assumed, k is determined according to the probability distribution and the required probability, and specifically, the probability distribution and the inclusion factor have a mapping relation, namely, when the probability distribution is normal distribution, k=2; when the probability distribution is uniform, k=1.732; when the probability distribution is an arcsine distribution, k=1.414.
In the embodiment of the invention, taking the verification point as 10m/s as an example, the uncertainty components of each standard corresponding to the wind speed sensor are evaluated and obtained as shown in table 4.
Table 4: uncertainty component value at 10m/s (unit: m/s)
Specifically, u 1 Component values expressed as pitot tube accuracy, the accuracy of the etalon should at least meet + -0.05% (interval half width 0.05%) within its speed measurement range, corresponding to interval half width a of 10m/s for example, according to the requirements of the regulations 1 =10m/s×0.0005=0.005 m/s, belongs to a normal distribution, and comprises a factor k=2; u (u) 2 The component value expressed as the pitot tube installation angle error, the maximum deflection angle of the airflow and the pitot tube pressure hole is +/-3 degrees, 3 degrees/360 degrees are=0.8 percent, the introduced error range is-0.8 percent (the interval half width is 0.8 percent), and the corresponding interval half width a is 10m/s as an example 2 =10m/s×0.0008=0.08 m/s, is uniformly distributed and contains a factor k 2 =1.732；u 3 Component values expressed as pitot tube obstruction errors, strut obstruction errors should be less than + -0.5% (interval half-width 0.5%), be uniformly distributed, contain a factor k 3 =1.732, for example 10m/s, corresponding to a half-width a of the interval 3 =10m/s×0.005=0.05m/s；u 4 Score expressed as micro differential pressure meter uncertaintyThe uncertainty of the digital micropressure meter is +/-0.8 Pa, the actual measured wind pressure value is 61.2Pa, taking 10m/s of the detection point as an example, and the allowable maximum relative error value is +/-0.8 Pa/61.28 Pa= +/-1.307% (the half width of the interval is 1.307%), and the corresponding half width of the interval is 10m/s as an examplea 4 =10m/s× 1.307% = 0.1307m/s, is normally distributed and contains a factor k 4 =2；u 5 Component values expressed as ambient temperature, u 6 Component values expressed as ambient humidity, u 7 In the present embodiment, the corresponding MPEV (Maximum Permissible Error absolute value, maximum allowable absolute value) of the ambient temperature, the ambient humidity and the ambient air pressure are ±maximum allowable, ±4%rh, ±0.2hpa, and the relative errors of the measured wind speed values when the environments are t=20h, p=1000hpa and u=50% rh are:
wherein ,a5 I.e. the half width of the interval of the ambient temperature, a 6 I.e. the half width of the interval of the ambient humidity, a 7 Namely the half width of the interval of the ambient air pressure, taking 10m/s as an example, the half width of the interval introduced by the temperature measurement errora 5 Interval half width introduced by humidity measurement error of 10m/s×0.068% =0.0068 m/sa 6 Section half width introduced by air pressure measurement error of =10m/s×0.018% =0.0018 m/sa 7 =10m/s×0.010% =0.0010 m/s, all of which are uniformly distributed and contain a factor k 5 =k 6 =k 7 =1.732；u 8 Component values expressed as reduction errors in computer processing, the resolution corresponding to the reduction results being 0.01m/s, and determining the possible value regions to be measuredM [ x ] Repair about -0.005，x Repair about +0.005]Correspond to a 8 =0.005 m/s, is uniformly distributed and comprises a factor k 8 =1.732；u 9 The component value expressed as the stability of the wind tunnel, namely the uncertainty component introduced by the stability of the wind tunnel, is exemplified by a detection point of 10m/s, the stability evaluation value of the wind tunnel is 0.3014%, the corresponding wind speed is 10m/s× 0.3014% = 0.03014m/s, the corresponding stability range is-0.01507 m/s to 0.01507m/s, the corresponding interval half widtha 9 0.01507m/s, satisfy the arcsine distribution, k in the interval 9 =1.414；u 10 The component value expressed as uniformity of the wind tunnel, namely, uncertainty component introduced by uniformity of the wind tunnel, is exemplified by a detection point of 10m/s, the uniformity evaluation value of the wind tunnel is 0.2467%, the corresponding wind speed is 10m/s× 0.2467% = 0.02467m/s, the corresponding uniformity range is-0.0123 m/s to 0.0123m/s, and the corresponding interval half width a 10 =0.0123 m/s, satisfying a uniform distribution, k 10 As can be seen from the above, the interval half width a= (stability evaluation value×verification point wind speed)/2 of stability, and the interval half width a= (uniformity evaluation value×verification point wind speed)/2 of uniformity.
And S04, calculating a wind speed synthesis standard uncertainty assessment result according to the component value calculation results of the class A parameters and the class B parameters and the correlation coefficient value.
Specifically, the calculation formula of the wind speed synthesis standard uncertainty assessment result is as follows:
wherein ,expressed as wind speed synthesis standard uncertainty assessment result, n 3 Expressed as the total number of B-class parameters, u f Component value expressed as f-th parameter of B-class parameters,>component values expressed as stability, < >>Component values expressed as uniformity, it will be appreciated that, in this embodiment, < +.>Namely u is as above 9 ，Namely u is as above 10 。
And S05, converting the wind speed synthesis standard uncertainty evaluation result into an expanded uncertainty, and outputting the expanded uncertainty.
Specifically, the final output result is composed ofThe product of the factor k and the product of the factor k is shown in Table 5.
Table 5: output result (Unit: m/s)
From the calculation results in table 5, it was found that the expansion uncertainty u=0.18 m/s when the correlation of the wind tunnel flow field was considered.
In summary, according to the wind speed measurement uncertainty evaluation method in the above embodiment of the present invention, by determining an evaluation parameter of wind speed measurement uncertainty according to an actual measurement situation, and classifying the evaluation parameter, a class a parameter for class a uncertainty component evaluation and a class B parameter for class B uncertainty component evaluation are obtained, where the evaluation parameter at least includes a stability coefficient and a uniformity coefficient of a wind tunnel; testing the stability and uniformity of the wind tunnel under the verification point wind speed, and calculating a correlation coefficient value according to the data of the stability and uniformity test of the wind tunnel; under the wind speed of a verification point, calculating component values of the class A parameter and the class B parameter according to the class A uncertainty component evaluation standard and the class B uncertainty component evaluation standard; calculating a wind speed synthesis standard uncertainty assessment result according to the component value calculation results and the related coefficient values of the class A parameters and the class B parameters; the wind speed synthesis standard uncertainty assessment result is converted into the expansion uncertainty and output, and specifically, the reliability of the wind speed measurement result can be effectively improved due to the fact that the correlation between the uniformity of the wind tunnel flow field and the stability of the wind tunnel flow field is fully considered, and the wind speed synthesis standard uncertainty assessment method has important significance in improving the quality of wind speed observation data and revealing weather rules and weather characteristics.
Example two
Referring to fig. 3, fig. 3 is a block diagram of a wind speed measurement uncertainty evaluation system 200 according to a second embodiment of the present invention, where the wind speed measurement uncertainty evaluation system 200 includes: an evaluation parameter determination module 21, a correlation coefficient value calculation module 22, a component value calculation module 23, an evaluation result calculation module 24, and an output module 25, wherein:
the evaluation parameter determining module 21 is configured to determine an evaluation parameter of uncertainty of wind speed measurement according to an actual measurement situation, and classify the evaluation parameter to obtain a class a parameter for evaluating a class a uncertainty component and a class B parameter for evaluating a class B uncertainty component, where the evaluation parameter at least includes stability and uniformity of a wind tunnel;
the correlation coefficient value calculation module 22 is configured to test stability and uniformity of the wind tunnel at the verification point wind speed, and calculate a correlation coefficient value according to data of the stability and uniformity test of the wind tunnel, where a calculation formula of the correlation coefficient value is:
wherein ,xd Representing the first wind speed value, y, measured the d-th time d A second wind speed value, n, representing the d-th measurement 1 Expressed as a first measurement total number of times, s (x) and s (y) respectively represent corresponding test standard deviations, Expressed as average of all first wind speeds,Expressed as the average of all second wind speeds,Expressed as a correlation coefficient value;
the component value calculating module 23 is configured to calculate, at the verification point wind speed, component values of the class a parameter and the class B parameter according to a class a uncertainty component evaluation criterion and a class B uncertainty component evaluation criterion, where a component value calculation formula of the class a parameter is:
wherein ,component values expressed as class A parameters, +.>The wind speed measurement error value expressed as the ith test,expressed as the average of wind speed measurement errors, n 2 Expressed as the second total number of measurements, +.>Display values of wind speed sensor expressed as ith test, r,/and so on>，Respectively representing the correction coefficient, the Pitot tube coefficient and the micropressure coefficient of the liquid in the digital micropressure meter,Expressed as the actual measured wind pressure of the digital micro manometer,t、p、u、e w the component values of the B-class parameters are respectively expressed as temperature, atmospheric pressure, humidity and saturated water vapor pressure, and the calculation formulas are as follows:
wherein ,component values expressed as class B parameters, +.>Half-width of interval expressed as class B parameter, < >>Expressed as inclusion factors;
the evaluation result calculation module 24 is configured to calculate a wind speed synthesis standard uncertainty evaluation result according to the component value calculation results of the class a parameter and the class B parameter and the correlation coefficient value, where a calculation formula of the wind speed synthesis standard uncertainty evaluation result is:
wherein ,expressed as wind speed synthesis standard uncertainty assessment result, n 3 Expressed as the total number of B-class parameters, u f Component value expressed as f-th parameter of B-class parameters,>component values expressed as stability, < >>Component values expressed as uniformity;
and the output module 25 is used for converting the wind speed synthesis standard uncertainty evaluation result into an expanded uncertainty and outputting the expanded uncertainty.
Further, in some alternative embodiments of the present invention, the evaluation parameter determining module 21 includes:
the judging unit is used for acquiring the evaluation parameters of the uncertainty of the wind speed measurement and sequentially judging whether each evaluation parameter is a parameter introduced by the repeatability measurement;
the first dividing unit is used for dividing the corresponding evaluation parameters into class A parameters when judging that the evaluation parameters are parameters introduced by the repeatability measurement;
and the second dividing unit is used for dividing the corresponding evaluation parameters into B-type parameters when judging that the evaluation parameters are not parameters introduced by the repeatability measurement.
Further, in some alternative embodiments of the present invention, the correlation coefficient value calculation module 22 includes:
the arrangement unit is used for installing the pitot tube at the center point of the section of the selected wind tunnel test section, and respectively arranging a plurality of wind measuring points at the upstream and downstream equidistant positions of the pitot tube, wherein the wind measuring points at one side of the pitot tube are uniformly distributed on the same plane;
The first control unit is used for controlling the blower motor to rotate at a preset rotating speed, and collecting a first wind speed value through a pitot tube arranged at the center point of the section of the selected wind tunnel test section after the flow speed is stable;
the second control unit is used for controlling the blower motor to rotate at a preset rotating speed, and after the flow speed is stable, the pitot tube is sequentially installed at each wind measuring point to acquire a second wind speed value;
and the correlation coefficient value calculation unit is used for calculating the correlation coefficient value according to the first wind speed value and the second wind speed value.
Example III
In another aspect, referring to fig. 4, an electronic device according to a third embodiment of the present invention includes a memory 20, a processor 10, and a computer program 30 stored in the memory and executable on the processor, where the processor 10 implements the wind speed measurement uncertainty assessment method described above 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. 4 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 program, when being executed by a processor, implements the wind speed measurement uncertainty assessment method as described above.
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 (4)
1. A method of assessing uncertainty in wind speed measurements, the method comprising:
according to actual measurement conditions, determining evaluation parameters of wind speed measurement uncertainty, and classifying the evaluation parameters to obtain class A parameters for class A uncertainty component evaluation and class B parameters for class B uncertainty component evaluation, wherein the evaluation parameters at least comprise an indication wind speed value of a measured sensor, a standard wind speed value, pitot tube accuracy, a pitot tube installation angle error, a pitot tube blocking error, a micro differential pressure meter uncertainty, a computer processing time repair error, an environment temperature measurement error, an environment humidity measurement error, an environment air pressure measurement error, an error introduced by the stability of a wind tunnel and an error introduced by uniformity;
Testing the stability and uniformity of the wind tunnel under the verification point wind speed, and calculating a correlation coefficient value according to the data of the stability and uniformity test of the wind tunnel;
and under the condition of verification point wind speed, calculating component values of the class A parameters and the class B parameters according to class A uncertainty component evaluation standards and class B uncertainty component evaluation standards, wherein the component value calculation formula of the class A parameters is as follows:
wherein ,component values expressed as class A parameters, +.>Wind speed measurement error value denoted as ith test,/">Expressed as the average of wind speed measurement errors, n 2 Expressed as the second total number of measurements, +.>Display values of wind speed sensor expressed as ith test, r,/and so on>，Respectively representing the correction coefficient, the Pitot tube coefficient and the micropressure coefficient of the liquid in the digital micropressure meter,Expressed as the actual measured wind pressure of the digital micro manometer,t、p、u、e w the component values of the B-class parameters are respectively expressed as temperature, atmospheric pressure, humidity and saturated water vapor pressure, and the calculation formulas are as follows:
wherein ,component values expressed as class B parameters, +.>Half-width of interval expressed as class B parameter, < >>Expressed as inclusion factors;
according to the component value calculation results of the A-type parameter and the B-type parameter and the correlation coefficient value, calculating a wind speed synthesis standard uncertainty evaluation result, wherein the calculation formula of the wind speed synthesis standard uncertainty evaluation result is as follows:
wherein ,expressed as wind speed synthesis standard uncertainty assessment result, n 3 Expressed as the total number of B-class parameters, u f Component value expressed as f-th parameter of B-class parameters,>component values expressed as stability, < >>Component values expressed as uniformity;
converting the wind speed synthesis standard uncertainty evaluation result into an expanded uncertainty, and outputting the expanded uncertainty;
the method for calculating the correlation coefficient value comprises the following steps of testing the stability and uniformity of the wind tunnel at the verification point wind speed and calculating the correlation coefficient value according to the data of the stability and uniformity test of the wind tunnel:
the pitot tube is arranged at the center point of the section of the selected wind tunnel test section, and a plurality of wind measuring points are respectively arranged at the upstream and downstream equidistant positions of the pitot tube, wherein the wind measuring points positioned at one side of the pitot tube are uniformly distributed on the same plane;
controlling the blower motor to rotate at a preset rotating speed, and collecting a first wind speed value through a pitot tube arranged at the center point of the section of the selected wind tunnel test section after the flow speed is stable;
controlling the blower motor to rotate at a preset rotating speed, and after the flow speed is stable, sequentially installing the pitot tubes to all wind measuring points to acquire a second wind speed value;
calculating the correlation coefficient value according to the first wind speed value and the second wind speed value, wherein the calculation formula of the correlation coefficient value is as follows:
wherein ,xd Representing the first wind speed value, y, measured the d-th time d A second wind speed value, n, representing the d-th measurement 1 Expressed as a first measurement total number of times, s (x) and s (y) respectively represent corresponding test standard deviations,expressed as average of all first wind speeds,Expressed as the average of all second wind speeds,Expressed as a correlation coefficient value;
the step of determining the evaluation parameters of the uncertainty of the wind speed measurement according to the actual measurement condition and classifying the evaluation parameters to obtain the class A parameters for evaluating the class A uncertainty component and the class B parameters for evaluating the class B uncertainty component comprises the following steps:
acquiring evaluation parameters of uncertainty of wind speed measurement, and sequentially judging whether each evaluation parameter is a parameter introduced by repeated measurement;
if yes, the corresponding assessment parameters are divided into class A parameters;
if not, the corresponding assessment parameters are divided into B-class parameters.
2. A wind speed measurement uncertainty assessment system, the system comprising:
the evaluation parameter determining module is used for determining an evaluation parameter of the uncertainty of wind speed measurement according to the actual measurement condition, classifying the evaluation parameter to obtain a class A parameter for evaluating a class A uncertainty component and a class B parameter for evaluating a class B uncertainty component, wherein the evaluation parameter at least comprises an indication wind speed value of a sensor to be measured, a standard wind speed value, pitot tube accuracy, a pitot tube installation angle error, a pitot tube blocking error, a micro differential pressure meter uncertainty, a computer processing time reduction error, an environment temperature measurement error, an environment humidity measurement error, an environment air pressure measurement error, an error introduced by the stability of a wind tunnel and an error introduced by uniformity;
The correlation coefficient value calculation module is used for testing the stability and uniformity of the wind tunnel under the wind speed of the verification point and calculating the correlation coefficient value according to the data of the stability and uniformity test of the wind tunnel;
the component value calculation module is used for respectively calculating the component values of the class A parameter and the class B parameter according to the class A uncertainty component evaluation standard and the class B uncertainty component evaluation standard at the verification point wind speed, and the component value calculation formula of the class A parameter is as follows:
wherein ,component values expressed as class A parameters, +.>Wind speed measurement error value denoted as ith test,/">Expressed as the average of wind speed measurement errors, n 2 Expressed as the second total number of measurements, +.>Display values of wind speed sensor expressed as ith test, r,/and so on>，Respectively representing the correction coefficient, the Pitot tube coefficient and the micropressure coefficient of the liquid in the digital micropressure meter,Expressed as the actual measured wind pressure of the digital micro manometer,t、p、u、e w the component values of the B-class parameters are respectively expressed as temperature, atmospheric pressure, humidity and saturated water vapor pressure, and the calculation formulas are as follows:
wherein ,component values expressed as class B parameters, +.>Half-width of interval expressed as class B parameter, < >>Expressed as inclusion factors;
The evaluation result calculation module is used for calculating a wind speed synthesis standard uncertainty evaluation result according to the component value calculation results of the class A parameters and the class B parameters and the related coefficient value, and the calculation formula of the wind speed synthesis standard uncertainty evaluation result is as follows:
wherein ,expressed as wind speed synthesis standard uncertainty assessment result, n 3 Expressed as the total number of B-class parameters, u f Component value expressed as f-th parameter of B-class parameters,>component values expressed as stability, < >>Component values expressed as uniformity;
the output module is used for converting the wind speed synthesis standard uncertainty evaluation result into an expansion uncertainty and outputting the expansion uncertainty;
the correlation coefficient value calculation module includes:
the arrangement unit is used for installing the pitot tube at the center point of the section of the selected wind tunnel test section, and respectively arranging a plurality of wind measuring points at the upstream and downstream equidistant positions of the pitot tube, wherein the wind measuring points at one side of the pitot tube are uniformly distributed on the same plane;
the first control unit is used for controlling the blower motor to rotate at a preset rotating speed, and collecting a first wind speed value through a pitot tube arranged at the center point of the section of the selected wind tunnel test section after the flow speed is stable;
The second control unit is used for controlling the blower motor to rotate at a preset rotating speed, and after the flow speed is stable, the pitot tube is sequentially installed at each wind measuring point to acquire a second wind speed value;
a correlation coefficient value calculating unit, configured to calculate the correlation coefficient value according to the first wind speed value and the second wind speed value, where a calculation formula of the correlation coefficient value is:
wherein ,xd Representing the first wind speed value, y, measured the d-th time d A second wind speed value, n, representing the d-th measurement 1 Expressed as a first measurement total number of times, s (x) and s (y) respectively represent corresponding test standard deviations,expressed as average of all first wind speeds,Expressed as the average of all second wind speeds,Expressed as a correlation coefficient value;
the evaluation parameter determination module includes:
the judging unit is used for acquiring the evaluation parameters of the uncertainty of the wind speed measurement and sequentially judging whether each evaluation parameter is a parameter introduced by the repeatability measurement;
the first dividing unit is used for dividing the corresponding evaluation parameters into class A parameters when judging that the evaluation parameters are parameters introduced by the repeatability measurement;
and the second dividing unit is used for dividing the corresponding evaluation parameters into B-type parameters when judging that the evaluation parameters are not parameters introduced by the repeatability measurement.
3. A computer readable storage medium, on which a computer program is stored, characterized in that the program, when being executed by a processor, implements the wind speed measurement uncertainty assessment method according to claim 1.
4. 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 wind speed measurement uncertainty assessment method of claim 1 when the program is executed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310879686.5A CN116609552B (en) | 2023-07-18 | 2023-07-18 | Wind speed measurement uncertainty assessment method, system, storage medium and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310879686.5A CN116609552B (en) | 2023-07-18 | 2023-07-18 | Wind speed measurement uncertainty assessment method, system, storage medium and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116609552A CN116609552A (en) | 2023-08-18 |
| CN116609552B true CN116609552B (en) | 2023-10-20 |
Family
ID=87682117
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310879686.5A Active CN116609552B (en) | 2023-07-18 | 2023-07-18 | Wind speed measurement uncertainty assessment method, system, storage medium and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116609552B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117150386B (en) * | 2023-10-31 | 2024-02-06 | 江西省气象探测中心 | Assessment method and device for measurement uncertainty of humidity sensor based on self-adaption |
| CN117949692B (en) * | 2024-03-27 | 2024-05-31 | 江西省气象探测中心 | Evaluation method and device for measurement uncertainty of photoelectric cup type anemometer |
| CN120831498B (en) * | 2025-09-18 | 2025-12-09 | 江西省气象探测中心 | Wind direction calibration result measurement uncertainty evaluation method and system |
Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1260481A (en) * | 2000-01-19 | 2000-07-19 | 山西省气象局 | Method and device for automatic detection and control of wind tunnel speed |
| CN1589370A (en) * | 2001-11-27 | 2005-03-02 | 艾劳埃斯·乌本 | Method for monitoring a sensor |
| EP1582876A1 (en) * | 2004-03-31 | 2005-10-05 | Klaus Stengler | Wind speed sensor |
| CN101034033A (en) * | 2007-01-16 | 2007-09-12 | 中国计量学院 | Wind tunnel calibration method for large flow gas pipeline averaging velocity tube flowmeter |
| KR20130099479A (en) * | 2012-02-29 | 2013-09-06 | 목포대학교산학협력단 | Method of sensorless mppt neural control for wind energy conversion systems |
| JP2014174113A (en) * | 2013-03-12 | 2014-09-22 | Tokyo Electric Power Co Inc:The | Method for predicting probability density distribution of repeated load caused to strung wire implement due to natural wind |
| CN204536363U (en) * | 2015-03-20 | 2015-08-05 | 天津量传计量检测技术有限公司 | A kind of standard wind speed detection system |
| CN108693377A (en) * | 2018-05-02 | 2018-10-23 | 中山大学 | A kind of temperature data modification method measuring three-dimensional velocity and ultrasonic temperature instrument |
| CN111090932A (en) * | 2019-12-10 | 2020-05-01 | 华能威宁风力发电有限公司 | A site calibration method for wind farms in operation adapting to moderately complex terrain |
| CN112683564A (en) * | 2020-12-25 | 2021-04-20 | 江苏中科君达物联网股份有限公司 | Method for calculating rated refrigerating capacity deviation of air conditioner at different wind speeds |
| CN112923967A (en) * | 2021-05-10 | 2021-06-08 | 北京康斯特仪表科技股份有限公司 | Instrument calibration system supporting rapid assessment of instrument calibration uncertainty |
| CN113221311A (en) * | 2021-02-02 | 2021-08-06 | 中国人民解放军国防科技大学 | Uncertainty quantification method for wind speed of atmospheric boundary layer |
| CN113790870A (en) * | 2021-11-18 | 2021-12-14 | 中国飞机强度研究所 | Indoor open type wind tunnel test flow field uniformity stability calibration method |
| CN113850440A (en) * | 2021-09-30 | 2021-12-28 | 浙江运达风电股份有限公司 | Wind speed prediction method using MCP based on average wind speed correction |
| JP2022074901A (en) * | 2020-11-05 | 2022-05-18 | 国立研究開発法人宇宙航空研究開発機構 | Information processing device, information processing method, and program |
| CN114707734A (en) * | 2022-04-12 | 2022-07-05 | 华润电力技术研究院有限公司 | Method, device and medium for predicting ultra-short-term wind power |
| CN114923547A (en) * | 2022-06-29 | 2022-08-19 | 成都秦川物联网科技股份有限公司 | An automatic evaluation device and method for error uncertainty of gas representation value |
| CN114962313A (en) * | 2022-06-07 | 2022-08-30 | 重庆通用工业(集团)有限责任公司 | Method for evaluating uncertainty of performance test of ventilator |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8781651B2 (en) * | 2012-09-21 | 2014-07-15 | Georgia Tech Research Corporation | Systems and methods providing a fuel-efficient RTA implementation with uncertain winds |
| DK3807529T3 (en) * | 2018-09-17 | 2024-03-25 | American Superconductor Corp | Yaw autocalibration for a wind turbine generator |
-
2023
- 2023-07-18 CN CN202310879686.5A patent/CN116609552B/en active Active
Patent Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1260481A (en) * | 2000-01-19 | 2000-07-19 | 山西省气象局 | Method and device for automatic detection and control of wind tunnel speed |
| CN1589370A (en) * | 2001-11-27 | 2005-03-02 | 艾劳埃斯·乌本 | Method for monitoring a sensor |
| EP1582876A1 (en) * | 2004-03-31 | 2005-10-05 | Klaus Stengler | Wind speed sensor |
| CN101034033A (en) * | 2007-01-16 | 2007-09-12 | 中国计量学院 | Wind tunnel calibration method for large flow gas pipeline averaging velocity tube flowmeter |
| KR20130099479A (en) * | 2012-02-29 | 2013-09-06 | 목포대학교산학협력단 | Method of sensorless mppt neural control for wind energy conversion systems |
| JP2014174113A (en) * | 2013-03-12 | 2014-09-22 | Tokyo Electric Power Co Inc:The | Method for predicting probability density distribution of repeated load caused to strung wire implement due to natural wind |
| CN204536363U (en) * | 2015-03-20 | 2015-08-05 | 天津量传计量检测技术有限公司 | A kind of standard wind speed detection system |
| CN108693377A (en) * | 2018-05-02 | 2018-10-23 | 中山大学 | A kind of temperature data modification method measuring three-dimensional velocity and ultrasonic temperature instrument |
| CN111090932A (en) * | 2019-12-10 | 2020-05-01 | 华能威宁风力发电有限公司 | A site calibration method for wind farms in operation adapting to moderately complex terrain |
| JP2022074901A (en) * | 2020-11-05 | 2022-05-18 | 国立研究開発法人宇宙航空研究開発機構 | Information processing device, information processing method, and program |
| CN112683564A (en) * | 2020-12-25 | 2021-04-20 | 江苏中科君达物联网股份有限公司 | Method for calculating rated refrigerating capacity deviation of air conditioner at different wind speeds |
| CN113221311A (en) * | 2021-02-02 | 2021-08-06 | 中国人民解放军国防科技大学 | Uncertainty quantification method for wind speed of atmospheric boundary layer |
| CN112923967A (en) * | 2021-05-10 | 2021-06-08 | 北京康斯特仪表科技股份有限公司 | Instrument calibration system supporting rapid assessment of instrument calibration uncertainty |
| CN113850440A (en) * | 2021-09-30 | 2021-12-28 | 浙江运达风电股份有限公司 | Wind speed prediction method using MCP based on average wind speed correction |
| CN113790870A (en) * | 2021-11-18 | 2021-12-14 | 中国飞机强度研究所 | Indoor open type wind tunnel test flow field uniformity stability calibration method |
| CN114707734A (en) * | 2022-04-12 | 2022-07-05 | 华润电力技术研究院有限公司 | Method, device and medium for predicting ultra-short-term wind power |
| CN114962313A (en) * | 2022-06-07 | 2022-08-30 | 重庆通用工业(集团)有限责任公司 | Method for evaluating uncertainty of performance test of ventilator |
| CN114923547A (en) * | 2022-06-29 | 2022-08-19 | 成都秦川物联网科技股份有限公司 | An automatic evaluation device and method for error uncertainty of gas representation value |
Non-Patent Citations (5)
| Title |
|---|
| An optimized deep nonlinear integrated framework for wind speed forecasting and uncertainty analysis;Jujie Wang;Applied Soft Computing;全文 * |
| Estimation of extreme wind speed based on upcrossing rate of mean wind speeds with Weibull distribution;Lin Da;Journal of Wind Engineering and Industrial Aerodynamics;全文 * |
| Wind damage estimation of roof sheathing panels considering directionality: Influences of both correlations of directional wind speeds and multiple response coefficients in each direction;Fengbo Wu;Journal of Wind Engineering & Industrial Aerodynamics;全文 * |
| 基于输入量相关性分析的PTB210型气压传感器的不确定度评定研究;闻春华;电子测量技术;第79页-第83页 * |
| 自动气象站风速传感器测量不确定度评定;曾涛;电子测量技术;第96页-第99页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116609552A (en) | 2023-08-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116609552B (en) | Wind speed measurement uncertainty assessment method, system, storage medium and device | |
| US8875567B2 (en) | Measurement of diesel engine emissions | |
| US20210164826A1 (en) | Detecting a change in a vibratory meter based on two baseline meter verifications | |
| CN114812758B (en) | Liquid level sensor calibration method, device, electronic device and storage medium | |
| US20190378181A1 (en) | Diagnosis cost output device, diagnosis cost output method, and computer-readable recording medium | |
| Gossweiler et al. | On fast-response probes: part 1—technology, calibration, and application to turbomachinery | |
| CN105387841B (en) | Height detecting device, load driving apparatus and height detection method | |
| CN114199318B (en) | Gas volume measurement methods, devices and storage media | |
| CN114923547A (en) | An automatic evaluation device and method for error uncertainty of gas representation value | |
| CN117686140B (en) | Method and system for evaluating measurement uncertainty of empty box barometer | |
| US20200133250A1 (en) | Method and system for providing quality controlled data from a redundant sensor system | |
| CN117949494B (en) | Evaluation method and device for measurement uncertainty of ultrasonic evaporation sensor | |
| CN106370779A (en) | Measurement deviation bridging method based on near disturbance source and environmental change | |
| CN120831498B (en) | Wind direction calibration result measurement uncertainty evaluation method and system | |
| CN120403713B (en) | Wave measuring instrument comparison error analysis method, device, equipment and medium | |
| CN112834208A (en) | An automatic pressure calibration system for safety valve calibration equipment | |
| CN112948770A (en) | Signal stability testing method and device, terminal equipment and system | |
| Shinder et al. | NIST’s new 3D airspeed calibration rig, addresses turbulent flow measurement challenges | |
| JP7518724B2 (en) | Evaluation display method, evaluation display device, and evaluation display program | |
| CN116429193B (en) | A method for measuring air flow in a limited space | |
| JP2002131439A (en) | Deterioration diagnosing device for semiconductor detector | |
| CN120084360B (en) | Radial Positioning Error Evaluation Method for Spatially Stabilized Inertial Navigation System | |
| JPH07225777A (en) | Method and device for determining inspection cycle | |
| Kumar et al. | Quantification and sensitivity analysis of uncertainties in turbocharger compressor gas stand measurements using monte carlo simulation | |
| CN120891221A (en) | Method for measuring flue gas velocity in a circular flue considering the velocity attenuation on the flue wall |
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 |