A sensor calibration method for refrigeration room based on logical self-consistency

Abstract

The present invention provides a sensor calibration method for a refrigeration room based on logical self-consistency, comprising: step 10, according to the sensor layout structure in the refrigeration room, establishing a sensor truth value logic association constraint system, and constructing a system logic inconsistency calculation function; step 20, constructing a sensor correction function; step 30, collecting sensor measurement data in the refrigeration room within a preset time period, and constructing a steady-state measurement data set; step 40, based on the steady-state measurement data set, taking the system logic inconsistency as the optimization target, optimizing the sensor correction function, and stopping the optimization after reaching the optimization cutoff condition, and obtaining the optimized sensor correction function; step 50, taking the correction value output by the optimized sensor correction function as the calibrated data. The present invention provides a sensor calibration method for a refrigeration room based on logical self-consistency, which corrects the sensor measurement data in the refrigeration room and improves the measurement accuracy of the sensor in the refrigeration room.

Description

Technical Field

[0001] The invention belongs to the technical field of sensor calibration, and in particular relates to a sensor calibration method for a refrigeration room based on logical self-consistency. Background Art

[0002] Refrigeration rooms need to continuously monitor various physical measurement data (such as temperature, pressure, flow, etc.), electrical signal data (such as current, power, etc.) and operating status data (such as control signals, valve opening, etc.). After a long period of operation, various sensors working under complex working conditions (such as high temperature, humidity, corrosion, etc.) are easily affected by factors such as equipment aging or noise interference, and are prone to typical failures such as complete failure, fixed deviation, drift deviation and accuracy reduction, which seriously affect the data quality of the monitoring system. For example, the temperature and flow sensors on the refrigeration system side are particularly important parameters in the operation of the refrigeration room. The accuracy of the sensor measurement directly affects the equipment fault diagnosis and operation control of the refrigeration room. The key to sensor calibration is to establish an accurate correspondence between the measured value and the true value (or reference value). At present, there are two mainstream methods: one is to calibrate through a standard sensor; the other is to estimate using a mathematical model. The former uses a "zero-dimensional" single-point constraint, which is usually carried out in a static test environment. Its accuracy guarantee is limited to specific conditions and is not enough to fully reflect the accuracy of the sensor under all operating conditions. The latter relies on "one-dimensional" equations or model constraints, and its accuracy is highly dependent on the correctness of the input parameters. The errors in the model itself may also affect the reliability of the benchmark value. Therefore, although the method of "static calibration and dynamic use" of benchmark values ​​or models seems feasible in some scenarios, this method may lead to significant deviations in complex systems operating under constantly changing internal and external environments, such as refrigeration rooms. Summary of the invention

[0003] The technical problem to be solved by the present invention is to provide a refrigeration room sensor calibration method based on logical self-consistency, which can correct the sensor measurement data in the refrigeration room and improve the measurement accuracy of the sensors in the refrigeration room.

[0004] To solve the above technical problems, the embodiments of the present invention adopt the following technical solutions:

[0005] The embodiment of the present invention provides a refrigeration room sensor calibration method based on logical self-consistency, comprising the following steps:

[0006] Step 10, according to the sensor layout structure in the refrigeration room, establish a sensor truth value logic association constraint system and construct a system logic inconsistency calculation function;

[0007] Step 20, constructing a sensor correction function;

[0008] Step 30, collecting measurement data of sensors in the refrigeration room within a preset time period to construct a steady-state measurement data set;

[0009] Step 40, based on the steady-state measurement data set, taking the system logic inconsistency as the optimization target, optimizing the sensor correction function, and stopping the optimization after reaching the optimization cutoff condition, to obtain the optimized sensor correction function;

[0010] Step 50: Using the correction value output by the optimized sensor correction function as calibrated data.

[0011] As a further improvement of the embodiment of the present invention, the sensor truth value logic association constraint system includes a chilled water return temperature consistency association constraint relationship, a chilled water quality conservation association constraint relationship and a chilled water energy conservation association constraint relationship.

[0012] As a further improvement of the embodiment of the present invention, the sensor truth value logic association constraint system also includes a cooling water return temperature consistency association constraint relationship, a cooling water quality conservation association constraint relationship and a cooling water energy conservation association constraint relationship.

[0013] As a further improvement of an embodiment of the present invention, the sensor truth value logic association constraint system includes a cooling water return temperature consistency association constraint relationship, a cooling water quality conservation association constraint relationship and a cooling water energy conservation association constraint relationship.

[0014] As a further improvement of the embodiment of the present invention, the expression of the chilled water return temperature consistency association constraint relationship is formula (1):

[0015] T chw,r1 =…=T chw,ri =…=T chw ,rn Formula (1)

[0016] Where, T chw,ri represents the chilled water return temperature of the i-th chiller, and n represents the number of chillers;

[0017] The expression of the conservation constraint relationship of chilled water quality is formula (2):

[0018]

[0019] Where M chw,0 Indicates the chilled water flow rate of the chilled water main, M chw,i represents the chilled water flow rate of the chilled water branch pipe where the i-th chiller is located, M chw,p Indicates the chilled water flow rate of the chilled water bypass pipe;

[0020] The expression of the associated constraint relationship of chilled water energy conservation is formula (3):

[0021]

[0022] Where, T chw,s0 Indicates the chilled water temperature of the chilled water supply main, T chw,r0 Indicates the chilled water temperature of the chilled water return main, T chw,si represents the chilled water supply temperature of the i-th chiller, Represents the average value of the chilled water supply temperature of all chillers, It represents the average value of chilled water return temperature of all chillers.

[0023] As a further improvement of the embodiment of the present invention, the expression of the cooling water return temperature consistency association constraint relationship is formula (4):

[0024] T cw,r0 =T cw,r1 =…=T cw,ri =…=T cw ,rn Formula (4)

[0025] Where, T cw,r0 Indicates the cooling water temperature of the cooling water return main, T cw,ri represents the cooling water return temperature of the i-th chiller;

[0026] The expression of the cooling water quality conservation constraint relationship is formula (5):

[0027]

[0028] Where M cw,0 Indicates the cooling water flow rate of the cooling water main, M cw,i represents the cooling water flow rate of the cooling water branch pipe where the i-th chiller is located, M cw,p Indicates the cooling water flow rate of the cooling water bypass pipe;

[0029] The expression of the cooling water energy conservation constraint relationship is formula (6):

[0030]

[0031] Where, T cw,s0 Indicates the cooling water temperature of the cooling water supply main, T cw,r0 Indicates the cooling water temperature of the cooling water return main, T cw,si represents the cooling water supply temperature of the i-th chiller, Represents the average cooling water supply temperature of all chillers, It represents the average value of cooling water return temperature of all chillers.

[0032] As a further improvement of the embodiment of the present invention, the expression of the system logic inconsistency calculation function is formula (7):

[0033]

[0034] Where D(LN) represents the system logic mismatch value, ω e represents the weight coefficient of the e-th association constraint relationship, E represents the total number of association constraint relationships in the sensor truth value logic association constraint system, and Y L,e The left side of the expression representing the calculated value of the e-th association constraint relationship based on the sensor correction value, Y R,e The right side of the expression representing the calculated value of the e-th association constraint relationship based on the sensor correction value, Y L,e,max The maximum value of the left side of the expression of the e-th association constraint relationship calculated based on the sensor correction value at all times, Y R,e,min Represents the minimum value among the calculated values ​​on the right side of the expression for the e-th association constraint relationship calculated based on the sensor correction value at all times.

[0035] As a further improvement of the embodiment of the present invention, the expression of the sensor correction function is formula (8):

[0036] Y=aX+b Formula (8)

[0037] Where Y represents the sensor correction value, a represents the linear deviation calibration factor, X represents the sensor measurement value, and b represents the fixed deviation calibration factor.

[0038] As a further improvement of the embodiment of the present invention, the step 30 specifically includes:

[0039] Step 301, collecting measurement data of sensors in the refrigeration room within a preset time period, and using a rolling time window to pre-process the measurement data of each sensor in different time periods;

[0040] Step 302, using formula (9) to calculate the coefficient of variation of the data of each sensor in each time period:

[0041]

[0042] In the formula, CV represents the coefficient of variation of the data, σ represents the standard deviation of the data, and μ represents the mean value;

[0043] Step 303 : Select the measurement data of the time period in which the data variation coefficients of all sensors are less than the coefficient threshold to form the steady-state measurement data.

[0044] As a further improvement of the embodiment of the present invention, in step 50, a gradient descent method or a genetic algorithm is used to iteratively optimize the sensor correction function to obtain an optimized sensor correction function.

[0045] Compared with the prior art, the technical solution of the present invention has the following beneficial effects:

[0046] An embodiment of the present invention provides a sensor calibration method for a refrigeration room based on logical self-consistency. First, according to the sensor layout structure in the refrigeration room, a sensor truth value logic association constraint system is established to construct a system logic incongruity calculation function; then a sensor correction function is constructed; then, measurement data of sensors in the refrigeration room within a preset time period are collected to construct a steady-state measurement data set; finally, based on the steady-state measurement data set, the system logic incongruity is used as an optimization target to optimize the sensor correction function, and the optimization is stopped after the optimization cutoff condition is reached to obtain an optimized sensor correction function; the correction value output by the optimized sensor correction function is used as the calibrated data. The method of the embodiment of the present invention constructs a sensor correction function to achieve data association between all sensor measurement values ​​and correction values, constructs a system logic incongruence calculation function, optimizes the correction parameters of the sensor correction function with the goal of minimizing the system logic incongruence value, obtains intermediate correction values ​​of all sensors, until these intermediate correction values ​​satisfy the sensor true value logical association constraint system to the greatest extent, determines the most suitable correction parameters, that is, obtains the most suitable sensor correction function, thereby correcting the sensor measurement data in the refrigeration room, obtaining calibrated data that is as close to the true value as possible, and improving the measurement accuracy of the sensors in the refrigeration room. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for use in the embodiments of the present invention will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative work.

[0048] Figure 1 is a flow chart of a method according to an embodiment of the present invention;

[0049] Figure 2 The diagram is a schematic diagram of a sensor layout structure of a refrigeration room in an embodiment of the present invention. DETAILED DESCRIPTION

[0050] The technical solution of the present invention is described in detail below in conjunction with the accompanying drawings.

[0051] Sensor calibration involves comparing the sensor measurement value with the "true value" (reference value) to determine and correct the error so that the output data meets the measurement requirements. In the normal operation of the computer room system, the true value and measurement value of the sensor are constantly changing, and the deviation relationship between the two is uncertain, including drift error, fixed error and random error. At present, the research on the calibration of a small number of sensors with fixed deviations is relatively sufficient, but in actual measurement scenarios, there are many types and quantities of sensors and the type and degree of deviation of each sensor are uncertain, which constitutes a major practical challenge for automatic sensor calibration.

[0052] To this end, an embodiment of the present invention provides a refrigeration room sensor calibration method based on logical self-consistency, such as Figure 1 As shown, the following steps are included:

[0053] Step 10: According to the sensor layout structure in the refrigeration room, a sensor truth value logic association constraint system is established, and a system logic inconsistency calculation function is constructed.

[0054] Step 20: construct a sensor correction function.

[0055] Step 30, collecting measurement data of sensors in the refrigeration room within a preset time period to construct a steady-state measurement data set.

[0056] Step 40, based on the steady-state measurement data set, taking the system logic inconsistency as the optimization target, optimizing the sensor correction function, and stopping the optimization after reaching the optimization cutoff condition, thereby obtaining the optimized sensor correction function.

[0057] Step 50: Using the correction value output by the optimized sensor correction function as the calibrated data.

[0058] The working conditions in the refrigeration room change over time, while the sensor error characteristics are fixed. The embodiment of the present invention constructs a sensor correction function to achieve data association between all sensor measurement values ​​and correction values, constructs a system logic inconsistency calculation function, takes the minimum system logic inconsistency value as the goal, optimizes the correction parameters of the sensor correction function, obtains the intermediate correction values ​​of all sensors, until these intermediate correction values ​​meet the sensor true value logical association constraint system to the greatest extent, determines the most suitable correction parameters, that is, obtains the most suitable sensor correction function, thereby calibrating the sensor measurement data in the refrigeration room, and obtaining calibrated data that is as close to the true value as possible, thereby improving the measurement accuracy of the sensors in the refrigeration room.

[0059] like Figure 1As shown, the sensor layout structure of the refrigeration room is as follows: in the chilled water system, the chilled water main includes a chilled water return main and a chilled water supply main. A chilled water branch pipe that is the same as the number of chillers is arranged between the chilled water return main and the chilled water supply main, and the evaporator of the chiller is arranged on the chilled water branch pipe. The section of the chilled water branch pipe connected to the chilled water return main is the chilled water return branch pipe, and the section connected to the chilled water supply main is the chilled water supply branch pipe. Temperature sensors and flow sensors are arranged on the chilled water return main and the chilled water return branch pipe, and temperature sensors and flow sensors are arranged on the chilled water supply branch pipe and the chilled water supply main pipe. A chilled water bypass pipe is arranged between the chilled water supply main and the chilled water return main, and the flow direction is from the chilled water supply main to the chilled water return main. A flow sensor is arranged on the chilled water bypass pipe.

[0060] In the cooling water system, the cooling water main includes a cooling water return main and a cooling water supply main. A number of cooling water branches that are the same as the number of chillers are arranged between the cooling water return main and the cooling water supply main. The condenser of the chiller is arranged on the cooling water branch. The section of the cooling water branch connected to the cooling water return main is the cooling water return branch, and the section connected to the cooling water supply main is the cooling water supply branch. Temperature sensors and flow sensors are provided on the cooling water return main and the cooling water return branch, and temperature sensors and flow sensors are provided on the cooling water supply branch and the cooling water supply main. A cooling water bypass is provided between the cooling water supply main and the cooling water return main. The flow direction is from the cooling water return main to the cooling water supply main, and a flow sensor is provided on the cooling water bypass.

[0061] In step 10, in the chilled water system of the refrigeration room, since the chilled water of the chilled water branch pipes of each chiller is from the same fully mixed chilled water return main pipe, theoretically, when the sensor has no error, the chilled water return temperature values ​​of all chillers turned on at the same time should be equal. Establish the chilled water return temperature consistency association constraint relationship shown in formula (1):

[0062] T chw,r1 =…=T chw,ri =…=T chw ,rn Formula (1)

[0063] Where, T chw,ri represents the chilled water return temperature of the i-th chiller, and n represents the number of chillers.

[0064] In the absence of leakage, the sum of the chilled water flow through each chilled water branch pipe plus the chilled water flow through the chilled water bypass pipe must match the flow in the chilled water main pipe, that is, the sum of the flow of all chilled water branches and the flow of the chilled water bypass pipe must be equal to the flow of the chilled water main pipe. Establish the chilled water mass conservation association constraint relationship shown in formula (2):

[0065]

[0066] Where M chw,0 Indicates the chilled water flow rate of the chilled water main, M chw,i represents the chilled water flow rate of the chilled water branch pipe where the i-th chiller is located, M chw,p Indicates the chilled water flow rate of the chilled water bypass pipe.

[0067] The chilled water from each chiller is mixed into the chilled water main and should satisfy the energy balance law:

[0068]

[0069] In the formula, Q chw,0 represents the cooling load of the chilled water main, Q chw,i represents the cooling load of the i-th chilled water branch, Q chw,p Indicates the cooling load of the chilled water bypass pipe.

[0070] Since the chilled water return temperature of the chiller is consistent and the flow rate satisfies the mass conservation, the associated constraint relationship of chilled water energy conservation shown in formula (3) is established:

[0071]

[0072] Where, T chw,s0 Indicates the chilled water temperature of the chilled water supply main, T chw,r0 Indicates the chilled water temperature of the chilled water return main, T chw,si represents the chilled water supply temperature of the i-th chiller, Represents the average value of the chilled water supply temperature of all chillers, It represents the average value of chilled water return temperature of all chillers.

[0073] In the cooling water system of the refrigeration room, since the cooling water of the cooling water branch pipes where each chiller is located comes from the same fully mixed cooling water return main pipe, theoretically, when the sensor has no error, the cooling water return temperature values ​​of all chillers turned on at the same time should be equal, and also equal to the cooling water return temperature of the cooling water return main pipe. Establish the cooling water return temperature consistency association constraint relationship shown in formula (4):

[0074] T cw,r0 =T cw,r1 =…=T cw,ri =…=T cw ,rn Formula (4)

[0075] Where, T cw,r0 Indicates the cooling water temperature of the cooling water return main, Tcw,ri Represents the cooling water return temperature of the i-th chiller.

[0076] In the absence of leakage, the sum of the cooling water flow through each cooling water branch pipe plus the cooling water flow of the cooling water bypass pipe must match the flow in the cooling water main pipe, that is, the sum of the flow of all cooling water branches and the flow of the cooling water bypass pipe must be equal to the flow of the cooling water main pipe. Establish the cooling water mass conservation constraint relationship shown in formula (5):

[0077]

[0078] Where M cw,0 Indicates the cooling water flow rate of the cooling water main, M cw,i represents the cooling water flow rate of the cooling water branch pipe where the i-th chiller is located, M cw,p Indicates the cooling water flow rate of the cooling water bypass pipe.

[0079] The cooling water from each chiller is mixed into the cooling water main and should satisfy the energy balance law:

[0080]

[0081] In the formula, Q cw,0 represents the cooling load of the cooling water main, Q cw,i represents the cooling load of the cooling water branch where the i-th chiller is located, Q cw,p Indicates the cooling load of the cooling water bypass pipe.

[0082] Since the chilled water return temperature of the chiller is consistent and the flow rate satisfies the mass conservation, the cooling water energy conservation associated constraint relationship shown in formula (6) is established:

[0083]

[0084] Where, T cw,s0 Indicates the cooling water temperature of the cooling water supply main, T cw,r0 Indicates the cooling water temperature of the cooling water return main, T cw,si represents the cooling water supply temperature of the i-th chiller, Represents the average cooling water supply temperature of all chillers, It represents the average value of cooling water return temperature of all chillers.

[0085] If only the temperature sensor and flow sensor in the chilled water system are calibrated, the chilled water return temperature consistency association constraint, chilled water quality conservation association constraint and chilled water energy conservation association constraint constitute the sensor truth value logic association constraint system.

[0086] If only the temperature sensor and flow sensor in the cooling water system are calibrated, the cooling water return temperature consistency association constraint, cooling water quality conservation association constraint and cooling water energy conservation association constraint constitute a multi-dimensional sensor truth value logic association constraint.

[0087] If both the temperature sensor and flow sensor on the chilled water system side and the temperature sensor and flow sensor on the cooling water system side are calibrated, the chilled water return temperature consistency association constraint, the chilled water quality conservation association constraint, the chilled water energy conservation association constraint, the cooling water return temperature consistency association constraint, the cooling water quality conservation association constraint and the cooling water energy conservation association constraint constitute a multi-dimensional sensor truth value logic association constraint.

[0088] Construct a system logic inconsistency calculation function. Specifically, in the sensor truth value logic association constraint system, the inconsistency of an association constraint relationship refers to the relative deviation of the unequal values ​​on the left and right sides caused by substituting the correction value into the association constraint relationship, as shown in the following formula:

[0089]

[0090] In the formula, x l,1 ,x l,2 ,…,x l,a The sensor correction value on the left side of the expression represents the association constraint relationship, x l,1 ,x l,2 ,…,x l,b The sensor correction value on the right side of the expression represents the association constraint relationship, Y L The calculated value on the left side of the expression representing the association constraint relationship, Y R The right side of the expression representing the association constraint relationship is calculated, D(LN e ) represents the logical inconsistency of the e-th association constraint relationship.

[0091] The total inconsistency of the entire refrigeration room is obtained by weighting the inconsistency of all associated constraints, that is, the expression of the system logic inconsistency calculation function is shown in formula (7):

[0092]

[0093] Where D(LN) represents the system logic mismatch value, ω e represents the weight coefficient of the e-th association constraint relationship, E represents the total number of association constraint relationships in the sensor truth value logic association constraint system, and Y L,e The left side of the expression representing the calculated value of the e-th association constraint relationship based on the sensor correction value, Y R,eThe right side of the expression representing the calculated value of the e-th association constraint relationship based on the sensor correction value, Y L,e,max The maximum value of the left side of the expression of the e-th association constraint relationship calculated based on the sensor correction value at all times, Y R,e,min Represents the minimum value among the calculated values ​​on the right side of the expression for the e-th association constraint relationship calculated based on the sensor correction value at all times.

[0094] In step 20, a sensor correction function is established to compensate for the linear drift deviation and fixed deviation in the sensor measurement, as shown in formula (8):

[0095] Y=aX+b Formula (8)

[0096] Where Y represents the sensor correction value, a represents the linear deviation calibration factor, X represents the sensor measurement value, and b represents the fixed deviation calibration factor.

[0097] The method of this embodiment integrates the measurement values ​​of all sensors to be calibrated into a measurement value matrix, the linear deviation calibration factors of all sensors to be calibrated into a linear deviation calibration matrix, and the fixed deviation calibration factors of all sensors to be calibrated into a fixed deviation calibration matrix to obtain a correction value matrix of all sensors to be calibrated.

[0098] Preferably, step 30 specifically includes:

[0099] Step 301, collect the measurement data of the sensors in the refrigeration room within a preset time period, and use a rolling time window to pre-process the measurement data of each sensor in different time periods. Preferably, the time period of the rolling time window is 30 minutes.

[0100] Step 302, using formula (9) to calculate the coefficient of variation of the data of each sensor in each time period:

[0101]

[0102] In the formula, CV represents the coefficient of variation of the data, σ represents the standard deviation of the data, and μ represents the mean value.

[0103] Step 303, select the measurement data of the time period in which the data variation coefficient of all sensors is less than the coefficient threshold to form a steady-state measurement data set. That is, only when the data variation coefficient of all sensors' measurement values ​​in the same time period is less than the coefficient threshold, the measurement values ​​of all sensors in the time period are added to the steady-state measurement data set. The coefficient threshold is set based on system experience and expected stability. Preferably, the coefficient threshold is 0.1.

[0104] The above embodiment uses the coefficient of variation as an indicator to measure the volatility of the measurement data, and the measurement data with a coefficient of variation less than the coefficient threshold constitute a steady-state measurement data set. The measurement data in the steady-state measurement data set has a small fluctuation, indicating that the system has reached a state of equilibrium and remains unchanged without external disturbances such as load changes, ambient temperature changes, etc., that is, a steady-state condition. The measurement parameters of the system in a steady-state condition satisfy the laws of physics. The steady-state measurement data set reflects the true performance of the system in an undisturbed state, and its data is more representative, which improves the accuracy of true value judgment.

[0105] In step 40, based on the steady-state measurement data set, with the system logic mismatch as the optimization target, a suitable machine learning optimization algorithm (such as gradient descent method or genetic algorithm) is used to continuously iterate and optimize the correction coefficient matrix until the system logic mismatch is minimized, at which time the sensor correction value is closest to its true value. The correction value output by the optimized sensor correction function is used as the calibrated data.

[0106] A specific example is provided below.

[0107] Figure 2 The refrigeration room shown in the figure has three chillers, the chilled water return main pipe and the three chilled water return branches are equipped with temperature sensors and flow sensors, and the chilled water supply branches and the three chilled water supply main pipes are equipped with temperature sensors. Figure 2 The fault-free simulation data of the temperature sensor and the flow sensor in the chilled water system shown in the figure are used. Errors (fixed deviation and linear deviation) are set on the basis of the fault-free simulation data. The simulation data after setting the errors are calibrated using the method of the present invention. After calibration, the data are compared with the fault-free simulation data. The results of error correction are shown in Table 1.

[0108] Table 1 Error settings and calibration results

[0109]

[0110] It can be seen from Table 1 that the method of the present invention can better identify the fixed deviation and linear deviation of the sensor, calibrate the sensor detection data, and the average calibration accuracy is greater than 90%, and the calibration effect is good.

[0111] The basic principles, main features and advantages of the present invention are shown and described above. Those skilled in the art should understand that the present invention is not limited to the above specific embodiments, and the descriptions in the above specific embodiments and the specification are only for further illustrating the principles of the present invention. Without departing from the spirit and scope of the present invention, the present invention may have various changes and improvements, and these changes and improvements all fall within the scope of the present invention claimed for protection.

Claims

1. A sensor calibration method for a refrigeration room based on logical self-consistency, characterized in that: The following steps are involved: Step 10, according to the sensor layout structure in the refrigeration room, establish a sensor truth value logic association constraint system and construct a system logic inconsistency calculation function; Step 20, constructing a sensor correction function; Step 30, collecting measurement data of sensors in the refrigeration room within a preset time period to construct a steady-state measurement data set; Step 40, based on the steady-state measurement data set, taking the system logic disharmony as the optimization target, optimizing the sensor correction function, and stopping the optimization after reaching the optimization cutoff condition, to obtain the optimized sensor correction function; Step 50, using the correction value output by the optimized sensor correction function as calibrated data; The sensor truth value logic association constraint system includes a chilled water return temperature consistency association constraint relationship, a chilled water quality conservation association constraint relationship, and a chilled water energy conservation association constraint relationship; Alternatively, the sensor truth value logic association constraint system includes a cooling water return temperature consistency association constraint relationship, a cooling water quality conservation association constraint relationship, and a cooling water energy conservation association constraint relationship; Alternatively, the sensor truth value logic association constraint system includes a chilled water return temperature consistency association constraint relationship, a chilled water quality conservation association constraint relationship, and a chilled water energy conservation association constraint relationship, and also includes a cooling water return temperature consistency association constraint relationship, a cooling water quality conservation association constraint relationship, and a cooling water energy conservation association constraint relationship; The expression of the consistency constraint relationship of the chilled water return temperature is formula (1): Formula (1) In the formula, Indicates The chilled water return temperature of each chiller, Indicates the number of chillers; The expression of the associated constraint relationship of chilled water quality conservation is formula (2): Formula (2) In the formula, Indicates the chilled water flow rate of the chilled water main. Indicates The chilled water flow rate of the chilled water branch pipe where the chiller is located, Indicates the chilled water flow rate of the chilled water bypass pipe; The expression of the associated constraint relationship of chilled water energy conservation is formula (3): Formula (3) In the formula, Indicates the chilled water temperature of the chilled water supply main. Indicates the chilled water temperature of the chilled water return main. Indicates The chilled water supply temperature of each chiller, Represents the average value of the chilled water supply temperature of all chillers, Indicates the average value of the chilled water return temperature of all chillers; The expression of the cooling water return temperature consistency constraint relationship is formula (4): Formula (4) In the formula, Indicates the cooling water temperature of the cooling water return main. Indicates Cooling water return temperature of each chiller; The expression of the cooling water quality conservation constraint relationship is formula (5): Formula (5) In the formula, Indicates the cooling water flow rate of the cooling water main. Indicates The cooling water flow rate of the cooling water branch pipe where the chiller is located, Indicates the cooling water flow rate of the cooling water bypass pipe; The expression of the cooling water energy conservation constraint is formula (6): Formula (6) In the formula, Indicates the cooling water temperature of the cooling water supply main. Indicates the cooling water temperature of the cooling water return main. Indicates The cooling water supply temperature of each chiller, Represents the average cooling water supply temperature of all chillers, It represents the average value of cooling water return temperature of all chillers.

2. The method for calibrating a sensor in a refrigeration room based on logical self-consistency according to claim 1, characterized in that: The expression of the system logic disharmony calculation function is formula (7): Formula (7) In the formula, Indicates the system logic inconsistency value, Indicates The weight coefficient of the associated constraint relationship, Represents the total number of associated constraint relationships in the sensor truth value logical associated constraint system, Indicates the first The left side of the expression of the associated constraint relationship is calculated. Indicates the first The right side of the expression of the associated constraint relationship is calculated. Indicates the first value calculated based on the sensor correction value at all times The maximum value among the values ​​calculated on the left side of the expression of the associated constraint relationship. Indicates the first value calculated based on the sensor correction value at all times The minimum value among the values ​​calculated on the right side of the expressions of the associated constraints.

3. The method for calibrating a sensor in a refrigeration room based on logical self-consistency according to claim 1, characterized in that: The expression of the sensor correction function is formula (8): Formula (8) In the formula, Indicates the sensor correction value, represents the linear deviation calibration factor, Represents the sensor measurement value, Represents the fixed bias calibration factor.

4. The method for calibrating a sensor in a refrigeration room based on logical self-consistency according to claim 1, characterized in that: The step 30 specifically includes: Step 301, collecting measurement data of sensors in the refrigeration room within a preset time period, and using a rolling time window to pre-process the measurement data of each sensor in different time periods; Step 302, using formula (9) to calculate the coefficient of variation of the data of each sensor in each time period: Formula (9) In the formula, represents the coefficient of variation of the data, represents the standard deviation of the data, represents the average value; Step 303 : Select the measurement data of the time period in which the data variation coefficients of all sensors are less than the coefficient threshold to form the steady-state measurement data.

5. The method for calibrating a sensor in a refrigeration room based on logical self-consistency according to claim 1, characterized in that: In step 50, the sensor correction function is iteratively optimized using a gradient descent method or a genetic algorithm to obtain an optimized sensor correction function.

Images