CN103868198B - Energy Efficiency Ratio modification method in a kind of wet and dry bulb temperature franchise - Google Patents

Energy Efficiency Ratio modification method in a kind of wet and dry bulb temperature franchise Download PDF

Info

Publication number
CN103868198B
CN103868198B CN201410081540.7A CN201410081540A CN103868198B CN 103868198 B CN103868198 B CN 103868198B CN 201410081540 A CN201410081540 A CN 201410081540A CN 103868198 B CN103868198 B CN 103868198B
Authority
CN
China
Prior art keywords
amp
energy efficiency
efficiency ratio
bulb temperature
indoor
Prior art date
Application number
CN201410081540.7A
Other languages
Chinese (zh)
Other versions
CN103868198A (en
Inventor
张忠斌
黄虎
张若楠
张敬坤
Original Assignee
南京师范大学
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 南京师范大学 filed Critical 南京师范大学
Priority to CN201410081540.7A priority Critical patent/CN103868198B/en
Publication of CN103868198A publication Critical patent/CN103868198A/en
Application granted granted Critical
Publication of CN103868198B publication Critical patent/CN103868198B/en

Links

Abstract

The invention discloses Energy Efficiency Ratio modification method in a kind of wet and dry bulb temperature franchise, the room air conditioner of different Energy Efficiency Ratio grades Energy Efficiency Ratio test result under cooling condition can be modified: gather 14 groups of room air conditioner Energy Efficiency Ratio data under operating condition of test, by each constant coefficient value in Bessel Formula certain Energy Efficiency Ratio regression equation, so that it is determined that Energy Efficiency Ratio drift correction formula; Gather model machine Energy Efficiency Ratio under any wet and dry bulb temperature in national standard, determine Energy Efficiency Ratio drift correction value by Energy Efficiency Ratio drift correction formula, thus obtaining the actual measurement operating mode recurrence correction value to nominal condition. The present invention provides a kind of Energy Efficiency Ratio Acquisition Error modification method being applied to performance of room air conditioners test and appraisal, the method can improve test and appraisal accuracy, reflect performance of room air conditioners really, it is to avoid manual control room air conditioner efficiency grade occurs predictable mistake to divide; Can instructing the setting of the air-conditioning equipment operational factors such as convertible frequency air-conditioner, the adjustment of space air state is more accurate simultaneously.

Description

Energy Efficiency Ratio modification method in a kind of wet and dry bulb temperature franchise

Technical field

The present invention relates to a kind of Energy Efficiency Ratio Acquisition Error modification method, a kind of method being applied to room air conditioner Energy Efficiency Ratio Acquisition Error correction specifically.

Background technology

Room air conditioner Energy Efficiency Standard is the Energy Efficiency Standard of the electrical equipment that China works out the earliest, and front and back experienced by three revisions, and the related specifications of the U.S., European Union and Japan has been followed in the formulation of Energy Efficiency Standard and revision. The revision of 2000 and twice Energy Efficiency Standard in 2004 is mainly reflected in the rise of energy efficiency market, and existing room air conditioner Energy Efficiency Standard, its efficiency grade is changed to three grades of efficiencies of existing GB GB12021.3-2010 by the Pyatyi efficiency of former GB GB12021.3-2004: existing Energy Efficiency Standard directly delete three in former Energy Efficiency Standard, four and Pyatyi efficiency grade, and former Energy Efficiency Standard I and II efficiency grade is newly defined as the two of existing Energy Efficiency Standard, three grades of efficiency grades. Based on this when unique step, the one-level first order as working standard that two grades of efficiencies of existing Energy Efficiency Standard are boosted.

In order to be appropriately determined the efficiency grade of room air conditioner, introducing thermodynamics consummating degree room air conditioner is carried out efficiency evaluation in " the thermodynamics consummating degree analysis of room air conditioner " literary composition, this is a kind of room air conditioner energy efficiency analysis method for air merited attention. Meanwhile, room air conditioner Energy Efficiency Standard updates revision, and accordingly, its performance standard is also required to revision, and matches, and namely avoids because measuring the former mistake division thus resulting in room air conditioner efficiency grade on to greatest extent.

Specified refrigerating capacity and Energy Efficiency Ratio need to reflect room air conditioner performance under nominal condition point, and this is also the foundation of room air conditioner efficiency ranking. When Cooling Capacity For The Room Air Conditioner and Energy Efficiency Ratio are surveyed, the impact of tested person condition and anthropic factor, the actual condition point of performance of room air conditioners test often offsets in tolerance. The franchise respectively ± 1 DEG C of correlation standard wet and dry bulb temperature and ± 0.5 DEG C.This test result often leading to refrigerating capacity and Energy Efficiency Ratio is not at nominal condition point, but records under a certain operating point in wet and dry bulb tolerance. Based on this test result, room air conditioner is carried out efficiency ranking, it is possible to the mistake division to room air conditioner efficiency grade can be caused. On the one hand, can reach under nominal condition (not up to) room air conditioner of three grades of energy efficiency market, it is likely to be due to actual measurement condition and is positioned at wet and dry bulb temperature tolerance lower limit (upper limit), and be divided into by mistake and can not reach (reaching) three grades of energy efficiency market, and three grades of energy efficiency market are the threshold values of the room air conditioner market access. On the other hand, along with instrument and meter precision improves constantly, the actual measurement condition point of manual control is positioned at wet and dry bulb temperature high tolerance or lower limit, it is also possible to cause that room air conditioner can occur predictable mistake to divide.

Summary of the invention

Technical problem: for the problems referred to above, the present invention provides one to be applied to room air conditioner, can obtain wet and dry bulb temperature in national standard tolerance under any operating point Energy Efficiency Ratio to the recurrence correction value of nominal condition, improve Energy Efficiency Ratio modification method in the wet and dry bulb temperature franchise of test and appraisal accuracy.

Technical scheme: Energy Efficiency Ratio modification method in the wet and dry bulb temperature franchise of the present invention, comprises the following steps:

1) the room air conditioner Energy Efficiency Ratio data under operating condition of test are gathered, described operating condition of test is 35 DEG C for keeping outside dry-bulb temperature, wet bulb temperature is 24 DEG C, inside difference control room, wet and dry bulb temperature is in 3 big classes, under 14 kinds of operating modes, carrying out the measurement of room air conditioner Energy Efficiency Ratio, recording Energy Efficiency Ratio with the deviation of Energy Efficiency Ratio under nominal condition is Δ Ei(i=1,2,14), i represents different measurement condition:

Wherein, first kind operating mode is that indoor wet bulb temperature is constant, change indoor dry-bulb temperature; Operating mode 1~5 indoor wet bulb temperature is 19 DEG C, indoor dry-bulb temperature corresponding 26 DEG C, 26.5 DEG C, 27 DEG C, 27.5 DEG C, 28 DEG C respectively;

Equations of The Second Kind operating mode is that indoor dry-bulb temperature is constant, change indoor wet bulb temperature; Operating mode 6~10 indoor dry-bulb temperature is 27 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 18.8 DEG C, 19 DEG C, 19.3 DEG C, 19.5 DEG C respectively;

3rd class is indoor wet and dry bulb temperature limit couple variations, and operating mode 11,12 indoor dry-bulb temperature is 26 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 19.5 DEG C respectively; Operating mode 13,14 indoor dry-bulb temperature is 28 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 19.5 DEG C respectively;

2) by described step 1) in gather each floor data substitute in following regression equation respectively, obtain the over-determined systems of 14 hexa-atomic linear function compositions, then solve described equation group with Bessel Formula for criterion, obtain every constant coefficient E of regression equation0, a, b, c, d, e value:

Δ E ( t 1 , t 2 ) = E 0 + at 1 + bt 2 + ct 1 2 + dt 2 2 + et 1 t 2 ,

Wherein, t1Represent indoor dry-bulb temperature and the temperature difference of indoor dry-bulb temperature, t under nominal condition point under tested operating mode2Represent indoor wet bulb temperature and the temperature difference of indoor wet bulb temperature under nominal condition point under tested operating mode;

E0Represent that in regression equation, constant term, a represent temperature difference item t of indoor dry-bulb temperature1Coefficient, b represent temperature difference item t of indoor wet bulb temperature2Coefficient, c be indoor dry-bulb temperature secondary temperature difference item t2 1Coefficient, d represent indoor wet bulb temperature secondary temperature difference item t2 2Coefficient, e represent indoor wet and dry bulb temperature coupling temperature difference item t1t2Coefficient;

3) the Energy Efficiency Ratio drift correction amount in wet and dry bulb temperature tolerance, under any practical running operating point, it is possible to calculate according to following formula:

ΔE(t1,t2)=E0+ζΛQt+(Qt)TΛ Qt;

Wherein, t is temperature difference vector, t=(t1,t2)T, T represents the transposition of vector;

Λ is weight matrix, Λ=diag (θ1, θ2), diag represents diagonal matrix, θ1Represent the indoor dry bulb temperature difference weighing factor to Energy Efficiency Ratio deviation, θ2Represent the indoor wet-bulb depression weighing factor to Energy Efficiency Ratio deviation;

Q is temperature transition matrix, Q = α 11 α 12 α 21 α 22 , α11Represent the dry bulb temperature difference influence coefficient to indoor dry-bulb temperature tolerance in indoor under measurement condition, α21Represent the dry bulb temperature difference influence coefficient to indoor wet bulb temperature tolerance in indoor under measurement condition, α12Represent the wet-bulb depression influence coefficient to indoor dry-bulb temperature tolerance in indoor under measurement condition, α22Represent the wet-bulb depression influence coefficient to indoor wet bulb temperature tolerance in indoor under measurement condition;

ζ is temperature difference item relative to the high order temperature difference item relative effect coefficient vector to Energy Efficiency Ratio, ζ=(ζ1, ζ2); ζ1Represent that a temperature term of the indoor dry bulb temperature difference is relative to the high order temperature term relative effect coefficient to Energy Efficiency Ratio deviation, ζ2Represent that a temperature term of indoor wet-bulb depression is relative to the high order temperature term relative effect coefficient to Energy Efficiency Ratio deviation;

4) according to described step 3) the Energy Efficiency Ratio drift correction amount Δ E (t that obtains1,t2) the Energy Efficiency Ratio Acquisition Error of tested room air conditioner is modified:

Et=Em-ΔE(t1,t2);

Wherein, EmFor the Energy Efficiency Ratio measured value under any operating mode of room air conditioner tested in wet and dry bulb tolerance; Δ E (t1,t2) for described step 3) and in Energy Efficiency Ratio drift correction amount under this operating mode calculated, EtCorrection value is returned for Energy Efficiency Ratio under this operating mode.

Energy Efficiency Ratio measured value under operating mode any in wet and dry bulb tolerance is modified by the present invention, obtain the recurrence correction value to nominal condition of the Energy Efficiency Ratio under this operating mode, can be used for the actual operation parameters adjustment of convertible frequency air-conditioner: when Energy Efficiency Ratio returns correction value more than Energy Efficiency Ratio measured value, suitably turn down compressor operating frequency; When Energy Efficiency Ratio returns correction value less than Energy Efficiency Ratio measured value, suitably heighten compressor operating frequency; When Energy Efficiency Ratio returns correction value equal to Energy Efficiency Ratio measured value, it is not necessary to adjust.

In the inventive method, described step 2) in Bessel Formula be:

ε=k σ, wherein ε represents monitoring residual error, and k represents monitoring coefficient, and σ represents the standard deviation of Energy Efficiency Ratio actual measurement deviation, σ = Σ i = 1 14 [ ΔE i - 1 14 Σ i = 1 14 ΔE i ] 13 , ( i = 1 , 2 , ... , 14 ) ;

Bessel Formula decision condition is: monitoring coefficient k value 0.65,2,3 successively, as k=0.65, has at least the Energy Efficiency Ratio actual measurement deviation of 50% to meet relational expression εi≤ ε; As k=2, the Energy Efficiency Ratio actual measurement deviation of 95% is had at least to meet relational expression εi≤ ε; As k=3, the Energy Efficiency Ratio actual measurement deviation of 100% is had at least to meet relational expression εi≤ ε;

εiRepresent the residual error of Energy Efficiency Ratio measurement error and equation deviation, ε under different operating pointi=| Δ Ei-ΔEi(t1,t2) |; Δ EiRepresent the difference of Energy Efficiency Ratio measured value and nominal value; Δ Ei(t1,t2) represent the difference of Energy Efficiency Ratio functional value and nominal value; I represents different measurement condition, i=1,2, and, 14.

In the inventive method, described step 3) in,

α 11 = e | λ 1 | + | λ 2 | / | | p 1 | | , α 12 = ( c - d - ( c - d ) 2 + e 2 ) | λ 1 | + | λ 2 | / | | p 1 | | ;

α 21 = e | λ 1 | + | λ 2 | / | | p 2 | | , α 22 = ( c - d + ( c - d ) 2 + e 2 ) | λ 1 | + | λ 2 | / | | p 2 | | ;

ζ 1 = ( aα 22 - bα 12 ) θ 1 ( α 11 α 22 - α 12 α 21 ) , ζ 2 = ( bα 11 - aα 21 ) θ 2 ( α 11 α 22 - α 12 α 21 ) ;

Wherein, λ1、λ2For the eigenvalue of secondary temperature difference term coefficient matrix B, p1、p2For the characteristic vector of secondary temperature difference term coefficient matrix B, B = c e / 2 e / 2 d , &lambda; 1 = c + d - ( c - d ) 2 + e 2 2 , &lambda; 2 = c + d + ( c - d ) 2 + e 2 2 , | &lambda; 1 | < | &lambda; 2 | ; | | p 1 | | = ( e 2 + &lsqb; c - d - ( c - d ) 2 + e 2 &rsqb; 2 ) 1 2 , | | p 2 | | = ( e 2 + &lsqb; c - d + ( c - d ) 2 + e 2 &rsqb; 2 ) 1 2 ) .

Beneficial effect: the present invention compared with prior art, has the advantage that

Energy Efficiency Ratio measured value in the wet and dry bulb temperature tolerance of the inventive method application finite discrete, describes the Energy Efficiency Ratio continuous print variation characteristic about wet and dry bulb temperature.

The standard deviation of deviation function deviation is surveyed for object with Energy Efficiency Ratio, with the numerical relation of test value and the residual sum standard deviation of equation value for target, each measurement condition is clicked on Mobile state monitoring, to determine the specific features and accuracy thereof that Energy Efficiency Ratio deviation changes with wet-bulb depression seriality.

The inventive method can obtain wet and dry bulb temperature, and in national standard tolerance, under any operating point, Energy Efficiency Ratio is to the recurrence correction value of nominal condition, thus instructing the setting of the air-conditioning equipment carrying out practically parameters such as convertible frequency air-conditioner, equipment runs more accurate;

The inventive method can obtain wet and dry bulb temperature in national standard tolerance under any operating point Energy Efficiency Ratio to the correction value of nominal condition, improve test and appraisal accuracy, reflect performance of room air conditioners more really, and provide reference to the revision of national standard;

Energy Efficiency Ratio drift correction formula is once determined, the room air conditioner Energy Efficiency Ratio measured value of follow-up identical Energy Efficiency Ratio grade is determined either directly through this correction formula, it is to avoid the huge consumption of the human and material resources that repeatedly repeated trials brings;

Along with instrument and meter precision improves constantly, it is to avoid the actual measurement condition point of manual control is positioned at wet and dry bulb temperature high tolerance or lower limit, thus causing that room air conditioner efficiency grade occurs predictable mistake to divide.

Accompanying drawing explanation

Fig. 1 is the process step figure of the inventive method.

Detailed description of the invention

Below in conjunction with drawings and Examples, the present invention will be further described:

In the wet and dry bulb temperature franchise of the present invention, the room air conditioner of different Energy Efficiency Ratio grades Energy Efficiency Ratio test result under cooling condition can be modified by Energy Efficiency Ratio modification method. It is modified to embodiment with the Energy Efficiency Ratio measured result under three class room air conditioner cooling conditions below to illustrate, specifically comprises the following steps that

Step 1) gather the room air conditioner Energy Efficiency Ratio data under operating condition of test, described operating condition of test is 35 DEG C for keeping outside dry-bulb temperature, wet bulb temperature is 24 DEG C, inside difference control room, wet and dry bulb temperature is in 3 big classes, under 14 kinds of operating modes, carrying out the measurement of room air conditioner Energy Efficiency Ratio, recording Energy Efficiency Ratio with the deviation of Energy Efficiency Ratio under nominal condition is

ΔE1=-0.03267kW`kW-1、ΔE2=-0.01971kW`kW-1、ΔE3=0kW`kW-1

ΔE4=0.0329kW`kW-1、ΔE5=0.05315kW`kW-1、ΔE6=-0.02574kW`kW-1

ΔE7=-0.0094kW`kW-1、ΔE8=0kW`kW-1、ΔE9=0.03773kW`kW-1

ΔE10=0.05315kW`kW-1、ΔE11=-0.06462kW`kW-1、ΔE12=0.01385kW`kW-1

ΔE13=-0.05007kW`kW-1、ΔE14=0.10285kW`kW-1:

Wherein, first kind operating mode is that indoor wet bulb temperature is constant, change indoor dry-bulb temperature; Operating mode 1~5 indoor wet bulb temperature is 19 DEG C, indoor dry-bulb temperature corresponding 26 DEG C, 26.5 DEG C, 27 DEG C, 27.5 DEG C, 28 DEG C respectively;

Equations of The Second Kind operating mode is that indoor dry-bulb temperature is constant, change indoor wet bulb temperature; Operating mode 6~10 indoor dry-bulb temperature is 27 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 18.8 DEG C, 19 DEG C, 19.3 DEG C, 19.5 DEG C respectively;

3rd class is indoor wet and dry bulb temperature limit couple variations, and operating mode 11,12 indoor dry-bulb temperature is 26 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 19.5 DEG C respectively; Operating mode 13,14 indoor dry-bulb temperature is 28 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 19.5 DEG C respectively;

Step 2) by described step 1) in gather each floor data substitute into following regression equation respectively:

&Delta; E ( t 1 , t 2 ) = E 0 + at 1 + bt 2 + ct 1 2 + dt 2 2 + et 1 t 2 ,

Wherein, t1Represent indoor dry-bulb temperature and the temperature difference of indoor dry-bulb temperature, t under nominal condition point under tested operating mode2Represent indoor wet bulb temperature and the temperature difference of indoor wet bulb temperature under nominal condition point under tested operating mode;

E0Represent that in regression equation, constant term, a represent temperature difference item t of indoor dry-bulb temperature1Coefficient, b represent temperature difference item t of indoor wet bulb temperature2Coefficient, c be indoor dry-bulb temperature secondary temperature difference itemCoefficient, d represent indoor wet bulb temperature secondary temperature difference itemCoefficient, e represent indoor wet and dry bulb temperature coupling temperature difference item t1t2Coefficient;

The 3rd, under 6 operating modes, t1、t2t1t2、ΔEi(t1,t2) value be 0, therefore obtain the over-determined systems of 12 hexa-atomic linear functions composition:

E0-a+c=-0.02574;

E0-0.5a+0.25c=-0.01278;

E0+ 0.5a+0.25c=0.0329;

E0+ a+c=0.05315;

E0-0.5b+0.25d=-0.02574;

E0-0.2b+0.04=-0.0094;

E0+ 0.3b+0.09d=0.03773;

E0+ 0.5b+0.25d=0.05315;

E0-a-0.5b+c+0.25d+0.5e=-0.06462;

E0-a+0.5b+c+0.25d-0.5e=0.01385;

E0+ a-0.5b+c+0.25d-0.5e=-0.05007;

E0+ a+0.5b+c+0.25d+0.5e=0.10285

Solve above-mentioned equation group, obtain every constant coefficient E of regression equation0, a, b, c, d, e value respectively 0.002,0.044,0.088,0.0049,0.019,0.0069;

Meanwhile, the convergence of Bessel Formula checking Energy Efficiency Ratio regression equation is adopted:

The residual error of test value and equation value is εi=| Δ Ei-ΔEi(t1,t2) |, Δ EiRepresent the difference of Energy Efficiency Ratio measured value and nominal value, Δ Ei(t1,t2) represent the difference of Energy Efficiency Ratio functional value and nominal value; It is respectively as follows: ε according to the residual error that above-mentioned formula tries to achieve each operating point corresponding1=0.00245, ε2=0.00121, ε3=0.000059, ε4=0.00106, ε5=0.001304, ε6=0.00357, ε7=0.00372, ε8=0.00023, ε9=0.00094, ε10=0.00128, ε11=0.00275, ε12=0.0024, ε13=0.00118 ε14=0.00132

σ represents the standard deviation of Energy Efficiency Ratio actual measurement deviation, &sigma; = &Sigma; i = 1 14 &lsqb; &Delta;E i - 1 14 &Sigma; i = 1 14 &Delta;E i &rsqb; 2 13 = 0.002 , ( i = 1 , 2 , ... , 14 ) ;

Bessel Formula is ε=k σ, respectively when monitoring coefficient k value 0.65,2,3 successively, as k=0.65, has the Energy Efficiency Ratio actual measurement deviation of 57% to meet relational expression εi≤ ε, (i represents different measurement condition, i=1,2, and, 14); As k=2, the Energy Efficiency Ratio actual measurement deviation of 100% is had to meet relational expression εi≤ ε; As k=3, the Energy Efficiency Ratio actual measurement deviation of 100% is had to meet relational expression εi≤ ε. Therefore the convergence of Energy Efficiency Ratio regression equation, W0, a, b, c, d, e value respectively 0.0019,0.19,0.37,0.0027,0.011,0.0098;

&Delta; W ( t 1 , t 2 ) = W 0 + at 1 + bt 2 + ct 1 2 + dt 2 2 + et 1 t 2

= 0.002 + 0.044 t 1 + 0.088 t 2 + 0.0049 t 1 2 + 0.019 t 2 2 + 0.0069 t 1 t 2

Step 3) in wet and dry bulb temperature tolerance, a room air conditioner model machine (specified refrigerating capacity and Energy Efficiency Ratio respectively 12kW, 3.32kW kW of the present invention-1) actual measurement Energy Efficiency Ratio Em=3.114kW`kW-1(operating condition of test point: indoor dry-bulb temperature is 26.7 DEG C, wet bulb temperature is 18.9 DEG C), then the Energy Efficiency Ratio drift correction amount under this operating mode, it is possible to calculate according to following formula:

ΔE(t1,t2)=E0+ζΛQt+(Qt)TΛ Qt;

Wherein, t is temperature difference vector, t=(t1,t2)T, T represents the transposition of vector;

Λ is weight matrix, Λ=diag (θ1, θ2), diag represents diagonal matrix, θ1Represent the indoor dry bulb temperature difference weighing factor to Energy Efficiency Ratio deviation, θ2Represent the indoor wet-bulb depression weighing factor to Energy Efficiency Ratio deviation;

Q is temperature transition matrix, Q = &alpha; 11 &alpha; 12 &alpha; 21 &alpha; 22 , α11Represent the dry bulb temperature difference influence coefficient to indoor dry-bulb temperature tolerance in indoor under measurement condition, α21Represent the dry bulb temperature difference influence coefficient to indoor wet bulb temperature tolerance in indoor under measurement condition, α12Represent the wet-bulb depression influence coefficient to indoor dry-bulb temperature tolerance in indoor under measurement condition, α22Represent the wet-bulb depression influence coefficient to indoor wet bulb temperature tolerance in indoor under measurement condition;

ζ is temperature difference item relative to the high order temperature difference item relative effect coefficient vector to Energy Efficiency Ratio, ζ=(ζ1, ζ2); ζ1Represent that a temperature term of the indoor dry bulb temperature difference is relative to the high order temperature term relative effect coefficient to Energy Efficiency Ratio deviation, ζ2Represent that a temperature term of indoor wet-bulb depression is relative to the high order temperature term relative effect coefficient to Energy Efficiency Ratio deviation;

θ1、θ2、α11、α12、α21、α22、ζ1、ζ2Determine by the following method:

Coefficient c, d, e of secondary temperature term in Energy Efficiency Ratio regression equation are carried out quadratic form and changes to obtain Second Order with Constant Coefficients matrix B = c e / 2 e / 2 d , Determine its eigenvalue λ1、λ2And characteristic of correspondence vector p1、p2:

&lambda; 1 = c + d + ( c - d ) 2 + e 2 2 , &lambda; 2 = c + d - ( c - d ) 2 + e 2 2 ,

p 1 = ( e , c - d - ( c - d ) 2 + e 2 ) T , p 2 = ( e , c - d + ( c - d ) 2 + e 2 ) T ;

Then, &theta; 1 = &lambda; 1 | &lambda; 1 | + | &lambda; 2 | , &theta; 2 = &lambda; 2 | &lambda; 1 | + | &lambda; 2 | , Weight matrix Λ=diag (θ1, θ2); &alpha; 11 = e | &lambda; 1 | + | &lambda; 2 | / | | p 1 | | , &alpha; 12 = ( c - d - ( c - d ) 2 + e 2 ) | &lambda; 1 | + | &lambda; 2 | / | | p 1 | | , &alpha; 21 = e | &lambda; 1 | + | &lambda; 2 | / | | p 2 | | , &alpha; 22 = ( c - d + ( c - d ) 2 + e 2 ) | &lambda; 1 | + | &lambda; 2 | / | | p 2 | | , Temperature transition matrix Q = &alpha; 11 &alpha; 12 &alpha; 21 &alpha; 22 ; &zeta; 1 = ( a&alpha; 22 - b&alpha; 12 ) &theta; 1 ( &alpha; 11 &alpha; 22 - &alpha; 12 &alpha; 21 ) , &zeta; 2 = ( b&alpha; 11 - a&alpha; 21 ) &theta; 2 ( &alpha; 11 &alpha; 22 - &alpha; 12 &alpha; 21 ) , Temperature difference item is relative to high order temperature difference item relative effect coefficient vector ζ=(ζ to Energy Efficiency Ratio1, ζ2);

Namely Energy Efficiency Ratio drift correction formula is:

&Delta; E ( t 1 , t 2 ) = 0.002 + ( 0.84 , 0.77 ) 0.16 0.84 0.15 - 0.036 0.036 0.15 t 1 t 2 ( 0.15 - 0.036 0.036 0.15 t 1 t 2 ) T 0.16 0.84 0.16 - 0.036 0.036 0.15 t 1 t 2

The Energy Efficiency Ratio drift correction amount under this operating mode is calculated according to above-mentioned solving result:

&Delta; E ( 26.7 , 18.9 ) = 0.002 + ( 0.84 , 0.77 ) 0.16 0.84 0.15 - 0.036 0.036 0.15 t 1 t 2 + ( 0.15 - 0.036 0.036 0.15 t 1 t 2 ) T 0.16 0.84 0.15 - 0.036 0.036 0.15 t 1 t 2 = - 0.01916 l W ` kW - 1 ;

Step 4) according to described step 3) the Energy Efficiency Ratio drift correction amount Δ E (t that obtains1,t2) the Energy Efficiency Ratio Acquisition Error of tested room air conditioner is modified:

Et=Em-ΔE(t1,t2)=3.0953kW`kW-1+0.01916kW`kW-1=3.1145kW`kW-1

Namely dry-bulb temperature be 26.7 DEG C, wet bulb temperature be that at 18.9 DEG C, Energy Efficiency Ratio returns correction value be 3.1145kW`kW-1

Wherein, EmFor the Energy Efficiency Ratio measured value under any operating mode of tested room air conditioner, for 3.0953kW`kW-1; Δ E (t1,t2) for described step 3) and in Energy Efficiency Ratio drift correction amount under this operating mode calculated, for-0.01916kW`kW-1; EtCorrection value is returned for Energy Efficiency Ratio under this operating mode.

Claims (3)

1. Energy Efficiency Ratio modification method in a wet and dry bulb temperature franchise, it is characterised in that the method comprises the following steps:
1) the room air conditioner Energy Efficiency Ratio data under operating condition of test are gathered, described operating condition of test is 35 DEG C for keeping outside dry-bulb temperature, wet bulb temperature is 24 DEG C, inside difference control room, wet and dry bulb temperature is in 3 big classes, under 14 kinds of operating modes, carrying out the measurement of room air conditioner Energy Efficiency Ratio, recording Energy Efficiency Ratio with the deviation of Energy Efficiency Ratio under nominal condition point is Δ Ei(i=1,2,14), i represents different measurement condition:
Wherein, first kind operating mode is that indoor wet bulb temperature is constant, change indoor dry-bulb temperature; Operating mode 1~5 indoor wet bulb temperature is 19 DEG C, indoor dry-bulb temperature corresponding 26 DEG C, 26.5 DEG C, 27 DEG C, 27.5 DEG C, 28 DEG C respectively;
Equations of The Second Kind operating mode is that indoor dry-bulb temperature is constant, change indoor wet bulb temperature; Operating mode 6~10 indoor dry-bulb temperature is 27 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 18.8 DEG C, 19 DEG C, 19.3 DEG C, 19.5 DEG C respectively;
3rd class is indoor wet and dry bulb temperature limit couple variations, and operating mode 11,12 indoor dry-bulb temperature is 26 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 19.5 DEG C respectively; Operating mode 13,14 indoor dry-bulb temperature is 28 DEG C, indoor wet bulb temperature corresponding 18.5 DEG C, 19.5 DEG C respectively;
2) by described step 1) in gather each floor data substitute in following regression equation respectively, obtain the over-determined systems of 14 hexa-atomic linear function compositions, then solve described equation group with Bessel Formula for criterion, obtain every constant coefficient E of regression equation0, a, b, c, d, e value:
&Delta; E ( t 1 , t 2 ) = E 0 + at 1 + bt 2 + ct 1 2 + dt 2 2 + et 1 t 2 ,
Wherein, t1Represent indoor dry-bulb temperature and the temperature difference of indoor dry-bulb temperature, t under nominal condition point under tested operating mode2Represent indoor wet bulb temperature and the temperature difference of indoor wet bulb temperature under nominal condition point under tested operating mode;
E0Represent that in regression equation, constant term, a represent temperature difference item t of indoor dry-bulb temperature1Coefficient, b represent temperature difference item t of indoor wet bulb temperature2Coefficient, c be indoor dry-bulb temperature secondary temperature difference item t2 1Coefficient, d represent indoor wet bulb temperature secondary temperature difference item t2 2Coefficient, e represent indoor wet and dry bulb temperature coupling temperature difference item t1t2Coefficient;
3) the Energy Efficiency Ratio drift correction amount in wet and dry bulb temperature tolerance, under any practical running operating point, it is possible to calculate according to following formula:
ΔE(t1,t2)=E0+ζΛQt+(Qt)TΛ Qt;
Wherein, t is temperature difference vector, t=(t1,t2)T, T represents the transposition of vector;
Λ is weight matrix, Λ=diag (θ1, θ2), diag represents diagonal matrix, θ1Represent the indoor dry bulb temperature difference weighing factor to Energy Efficiency Ratio deviation, θ2Represent the indoor wet-bulb depression weighing factor to Energy Efficiency Ratio deviation;
Q is temperature transition matrix, Q = &alpha; 11 &alpha; 12 &alpha; 21 &alpha; 22 , α11Represent the dry bulb temperature difference influence coefficient to indoor dry-bulb temperature tolerance in indoor under measurement condition, α21Represent the dry bulb temperature difference influence coefficient to indoor wet bulb temperature tolerance in indoor under measurement condition, α12Represent the wet-bulb depression influence coefficient to indoor dry-bulb temperature tolerance in indoor under measurement condition, α22Represent the wet-bulb depression influence coefficient to indoor wet bulb temperature tolerance in indoor under measurement condition;
ζ is temperature difference item relative to the high order temperature difference item relative effect coefficient vector to Energy Efficiency Ratio, ζ=(ζ1, ζ2); ζ1Represent that a temperature term of the indoor dry bulb temperature difference is relative to the high order temperature term relative effect coefficient to Energy Efficiency Ratio deviation, ζ2Represent that a temperature term of indoor wet-bulb depression is relative to the high order temperature term relative effect coefficient to Energy Efficiency Ratio deviation;
4) according to described step 3) the Energy Efficiency Ratio drift correction amount Δ E (t that obtains1,t2) the Energy Efficiency Ratio Acquisition Error of tested room air conditioner is modified:
Et=Em-ΔE(t1,t2);
Wherein, EmFor the Energy Efficiency Ratio measured value under any operating mode of room air conditioner tested in wet and dry bulb tolerance; Δ E (t1,t2) for described step 3) and in Energy Efficiency Ratio drift correction amount under this operating mode calculated, EtCorrection value is returned for Energy Efficiency Ratio under this operating mode.
2. Energy Efficiency Ratio modification method in wet and dry bulb temperature franchise according to claim 1, it is characterised in that described step 2) in Bessel Formula be: ε=k σ, wherein ε represents monitoring residual error, k represents monitoring coefficient, and σ represents the standard deviation of Energy Efficiency Ratio actual measurement deviation &sigma; = &Sigma; i = 1 14 &lsqb; &Delta;E i - 1 14 &Sigma; i = 1 14 &Delta;E i &rsqb; 2 13 , ( i = 1 , 2 , ... , 14 ) ;
Described Bessel Formula decision condition is: monitoring coefficient k value 0.65,2,3 successively, as k=0.65, has at least the Energy Efficiency Ratio actual measurement deviation of 50% to meet relational expression εi≤ ε; As k=2, the Energy Efficiency Ratio actual measurement deviation of 95% is had at least to meet relational expression εi≤ ε; As k=3, the Energy Efficiency Ratio actual measurement deviation of 100% is had at least to meet relational expression εi≤ ε;
εiRepresent the measurement error of Energy Efficiency Ratio under different operating mode and the residual error of equation deviation, εi=| Δ Ei-ΔEi(t1,t2) |; Δ EiRepresent the difference of Energy Efficiency Ratio measured value and nominal value; Δ Ei(t1,t2) represent the difference of Energy Efficiency Ratio functional value and nominal value; I represents different measurement condition, i=1,2,14.
3. Energy Efficiency Ratio modification method in wet and dry bulb temperature franchise according to claim 1, it is characterised in that described step 3) in, &theta; 1 = &lambda; 1 | &lambda; 1 | + | &lambda; 2 | , &theta; 2 = &lambda; 2 | &lambda; 1 | + | &lambda; 2 | ;
&alpha; 11 = e | &lambda; 1 | + | &lambda; 2 | / | | p 1 | | , &alpha; 12 = ( c - d - ( c - d ) 2 + e 2 ) | &lambda; 1 | + | &lambda; 2 | / | | p 1 | | ;
&alpha; 21 = e | &lambda; 1 | + | &lambda; 2 | / | | p 2 | | , &alpha; 22 = ( c - d - ( c - d ) 2 + e 2 ) | &lambda; 1 | + | &lambda; 2 | / | | p 2 | | ;
&zeta; 1 = ( a&alpha; 22 - b&alpha; 12 ) &theta; 1 ( &alpha; 11 &alpha; 22 - &alpha; 12 &alpha; 21 ) , &zeta; 2 = ( b&alpha; 11 - a&alpha; 21 ) &theta; 2 ( &alpha; 11 &alpha; 22 - &alpha; 12 &alpha; 21 ) ;
Wherein, λ1、λ2For the eigenvalue of secondary temperature difference term coefficient matrix B, p1、p2For the characteristic vector of secondary temperature difference term coefficient matrix B, B = c e / 2 e / 2 d , &lambda; 1 = c + d - ( c - d ) 2 + e 2 2 , &lambda; 2 = c + d + ( c - d ) 2 + e 2 2 , 1| < | λ2|; | | p 1 | | = ( e 2 + &lsqb; c - d - ( c - d ) 2 + e 2 &rsqb; 2 ) 1 2 , | | p 2 | | = ( e 2 + &lsqb; c - d + ( c - d ) 2 + e 2 &rsqb; 2 ) 1 2 ) .
CN201410081540.7A 2014-03-06 2014-03-06 Energy Efficiency Ratio modification method in a kind of wet and dry bulb temperature franchise CN103868198B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410081540.7A CN103868198B (en) 2014-03-06 2014-03-06 Energy Efficiency Ratio modification method in a kind of wet and dry bulb temperature franchise

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410081540.7A CN103868198B (en) 2014-03-06 2014-03-06 Energy Efficiency Ratio modification method in a kind of wet and dry bulb temperature franchise

Publications (2)

Publication Number Publication Date
CN103868198A CN103868198A (en) 2014-06-18
CN103868198B true CN103868198B (en) 2016-06-15

Family

ID=50907007

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410081540.7A CN103868198B (en) 2014-03-06 2014-03-06 Energy Efficiency Ratio modification method in a kind of wet and dry bulb temperature franchise

Country Status (1)

Country Link
CN (1) CN103868198B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107084480A (en) * 2017-04-13 2017-08-22 青岛海尔空调器有限总公司 Air-conditioning heating efficiency when heating capacity online test method
CN107084477B (en) * 2017-04-13 2019-12-31 青岛海尔空调器有限总公司 Heating operation control method for air conditioner

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4344565A (en) * 1979-08-01 1982-08-17 Nippondenso Co., Ltd. Control method and apparatus for air conditioners
CN102788005A (en) * 2012-08-02 2012-11-21 青岛海信日立空调系统有限公司 Method for controlling the operation of the compressors in optimized VRV air conditioning system and device thereof
CN103542489A (en) * 2012-07-12 2014-01-29 珠海格力节能环保制冷技术研究中心有限公司 Control method, device and system of air conditioner system

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7886984B2 (en) * 2007-12-19 2011-02-15 Institute For Information Industry Method of utilizing air conditioner to control thermal comfort level of environment

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4344565A (en) * 1979-08-01 1982-08-17 Nippondenso Co., Ltd. Control method and apparatus for air conditioners
CN103542489A (en) * 2012-07-12 2014-01-29 珠海格力节能环保制冷技术研究中心有限公司 Control method, device and system of air conditioner system
CN102788005A (en) * 2012-08-02 2012-11-21 青岛海信日立空调系统有限公司 Method for controlling the operation of the compressors in optimized VRV air conditioning system and device thereof

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
单元式空调机季节能耗的计算与分析;沈伟刚;《暖通空调》;20120129;第42卷(第3期);全文 *

Also Published As

Publication number Publication date
CN103868198A (en) 2014-06-18

Similar Documents

Publication Publication Date Title
Nguyen et al. A performance comparison of sensitivity analysis methods for building energy models
CN106352484B (en) One drags more air-conditionings and its refrigeration control method
Kissock et al. Ambient-temperature regression analysis for estimating retrofit savings in commercial buildings
Wang et al. Simplified building model for transient thermal performance estimation using GA-based parameter identification
Li et al. Decoupling features and virtual sensors for diagnosis of faults in vapor compression air conditioners
Hydeman et al. Development and testing of a reformulated regression-based electric chiller model
US8738185B2 (en) Altitude adjustment for heating, ventilating and air conditioning systems
CN201673133U (en) Heat transfer coefficient detector of building retaining structure
US20180313559A1 (en) Economizer/dcv controller with manual sensor calibration
Li et al. A review of virtual sensing technology and application in building systems
Ye et al. A new approach for measuring predicted mean vote (PMV) and standard effective temperature (SET∗)
US7200468B2 (en) System for determining overall heating and cooling system efficienies
Li et al. Development, evaluation, and demonstration of a virtual refrigerant charge sensor
Gu Airflow network modeling in EnergyPlus
Iragorry et al. A critical review of properties and models for frost formation analysis
Nassif et al. A new operating strategy for economizer dampers of VAV system
CN202770606U (en) Ventilation air-conditioning system air leakage automatic testing device
CN102486630B (en) Case reasoning technology-based damping machine intelligent optimized control method
Liu et al. Simplified thermal modeling of indirect evaporative heat exchangers
Wemhoff et al. Predictions of energy savings in HVAC systems by lumped models
CN102934036A (en) Model based system and method for estimating parameters and states in temperature controlled spaces
US20130013117A1 (en) Methods and apparatus for differential energy based airside economizer changeover
US20130345995A1 (en) Air Flow Control And Power Usage Of An Indoor Blower In An HVAC System
CN106123255B (en) Air conditioning control method and air-conditioning
CN106324521B (en) A kind of method of Combined estimator electrokinetic cell system parameter and state-of-charge

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model