CN103868198A - Energy efficiency ratio correcting method in wet and dry bulb temperature tolerance - Google Patents

Abstract

The invention discloses an energy efficiency ratio correcting method in wet and dry bulb temperature tolerance. The method can be used for correcting energy efficiency ratio test results of room air conditioners with different energy efficiency ratios under the refrigeration work conditions and comprises the following steps that the energy efficiency ratio data of fourteen groups of room air conditioners under the test work conditions is collected, and each constant coefficient value in the energy efficiency ratio regression equation is determined through a Bessel formula, so the energy efficiency ratio deviation correction value is determined; the energy efficiency ratios of sample machines at any wet and dry bulb temperature in the national standard are collected, and the energy efficiency ratio deviation correcting value is determined through the energy efficiency ratio deviation correcting equation, so the regression correcting value from the actually measured work condition to the nominal work condition is obtained. The invention provides the energy efficiency ratio collection error correcting method applied to room air conditioner performance evaluation. The method has the advantages that the evaluation accuracy can be improved, the performance of room air conditioners can be really reflected, and the predictable error division of the air conditioner energy efficiency level during artificial control is avoided; meanwhile, the operation parameter setting of air conditioning equipment such as variable frequency air conditioners can be guided, and the regulation of the air state in the room is more precise.

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, is a kind of method that is applied to the correction of room air conditioner Energy Efficiency Ratio Acquisition Error specifically.

Background technology

Room air conditioner Energy Efficiency Standard is the Energy Efficiency Standard of the consumer worked out the earliest of China, and front and back have experienced 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 efficiency limit value, 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 is directly deleted three, four and Pyatyi efficiency grade in former Energy Efficiency Standard, and former Energy Efficiency Standard I and II efficiency grade is redefined to two, three grades of efficiency grades into existing Energy Efficiency Standard.Waiting under the condition of step-length based on this, using the existing Energy Efficiency Standard secondary efficiency first order of one-level as current standard of boosting.

For the correct efficiency grade of judging room air conditioner, in " the thermodynamics sophistication of room air conditioner is analyzed " literary composition, to introduce thermodynamics sophistication room air conditioner is carried out to efficiency evaluation, this is a kind of room air conditioner energy efficiency analysis method for air meriting attention.Meanwhile, room air conditioner Energy Efficiency Standard upgrades revision, corresponding, and its performance standard also needs revision, and matches, and on to greatest extent, avoids because measuring former thereby causing the mistake of room air conditioner efficiency grade to divide.

Specified refrigerating capacity and Energy Efficiency Ratio need reflect the performance of room air conditioner 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 human factor, the actual condition point of performance of room air conditioners test tends to be offset in tolerance.Correlation standard the franchise of wet and dry bulb temperature be respectively ± 1 ℃ and ± 0.5 ℃.This test result that often causes refrigerating capacity and Energy Efficiency Ratio is not at nominal condition point, but under a certain operating point, records in wet and dry bulb tolerance.Based on this test result, room air conditioner is carried out to efficiency ranking, may cause the mistake of room air conditioner efficiency grade to be divided.On the one hand, under nominal condition, can reach the room air conditioner of (not reaching) three grades of efficiency limit values, may test operating mode and be positioned at wet and dry bulb temperature tolerance lower limit (upper limit) due to reality, can not reach (reaching) three grades of efficiency limit values and be divided into by mistake, and three grades of threshold values that efficiency limit value is the room air conditioner market access.On the other hand, along with instrument and meter precision improves constantly, the actual test of manual control operating point is positioned at wet and dry bulb temperature high tolerance or lower limit, also may cause room air conditioner that predictable wrong division can occur.

Summary of the invention

Technical problem: for the problems referred to above, the invention provides one and be applied to room air conditioner, can obtain wet and dry bulb temperature in national standard tolerance arbitrarily under 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 wet and dry bulb temperature franchise of the present invention, comprises the following steps:

1) be captured in the room air conditioner Energy Efficiency Ratio data under operating condition of test, described operating condition of test is that dry-bulb temperature is 35 ℃ outside holding chamber, wet-bulb temperature is 24 ℃, inside difference control room, wet and dry bulb temperature is in 3 large classes, under 14 kinds of operating modes, carry out the measurement of room air conditioner Energy Efficiency Ratio, the deviation that records Energy Efficiency Ratio under Energy Efficiency Ratio and nominal condition is Δ E _i(i=1,2,, 14), i represents different test operating modes:

Wherein, first kind operating mode is that indoor wet-bulb temperature is constant, changes indoor dry-bulb temperature; Operating mode 1～5 indoor wet-bulb temperature is 19 ℃, corresponding 26 ℃, 26.5 ℃, 27 ℃, 27.5 ℃, 28 ℃ respectively of indoor dry-bulb temperatures;

Equations of The Second Kind operating mode is that indoor dry-bulb temperature is constant, changes indoor wet-bulb temperature; Operating mode 6～10 indoor dry-bulb temperatures are 27 ℃, corresponding 18.5 ℃, 18.8 ℃, 19 ℃, 19.3 ℃, 19.5 ℃ respectively of indoor wet-bulb temperature;

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

2) by the following regression equation of each floor data difference substitution gathering in described step 1), obtain the overdetermined equation group of 14 hexa-atomic linear function compositions, then take Bessel Formula as criterion solves described equation group, obtain every constant coefficient E of regression equation ₀, a, b, c, d, e value:

ΔE(t ₁,t ₂)=E ₀+at ₁+bt ₂+ct ₁ ²+dt ₂ ²+et ₁t ₂，

Wherein, t ₁represent indoor dry-bulb temperature and the nominal condition point temperature difference of indoor dry-bulb temperature down under tested operating mode, t ₂represent indoor wet-bulb temperature and the nominal condition point temperature difference of indoor wet-bulb temperature down under tested operating mode;

E ₀represent that in regression equation, constant term, a represent temperature difference item t of indoor dry-bulb temperature ₁coefficient, b represent temperature difference item t of indoor wet-bulb temperature ₂coefficient, c be indoor dry-bulb temperature secondary temperature difference item t ² ₁coefficient, d represent indoor wet-bulb temperature secondary temperature difference item t ² ₂coefficient, e represent indoor wet and dry bulb temperature coupling temperature difference item t ₁t ₂

Coefficient;

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

ΔE(t ₁,t ₂)=E ₀+ζΛQt+(Qt) ^TΛQt；

Wherein, t is temperature difference vector, t=(t ₁, t ₂) ^t, T represents vectorial transposition;

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

Q is temperature transition matrix, $Q=\left(\begin{array}{cc}{\mathrm{\α}}_{11}& {\mathrm{\α}}_{12}\\ {\mathrm{\α}}_{21}& {\mathrm{\α}}_{22}\end{array}\right),$ α ₁₁represent the influence coefficient of the indoor dry bulb temperature difference to indoor dry-bulb temperature tolerance under test operating mode, α ₂₁represent the influence coefficient of the indoor dry bulb temperature difference to indoor wet bulb bulb temperature tolerance under test operating mode, α ₁₂represent the influence coefficient of indoor wet-bulb depression to indoor dry-bulb temperature tolerance under test operating mode, α ₂₂represent the influence coefficient of indoor wet-bulb depression to indoor wet-bulb temperature tolerance under test operating mode;

ζ be temperature difference item with respect to the high order temperature difference item relative effect coefficient vector to Energy Efficiency Ratio, ζ=(ζ ₁, ζ ₂); ζ ₁a temperature term of the expression indoor dry bulb temperature difference is the relative effect coefficient to Energy Efficiency Ratio deviation with respect to high order temperature term, ζ ₂a temperature term of expression indoor wet-bulb depression is the relative effect coefficient to Energy Efficiency Ratio deviation with respect to high order temperature term;

4) the Energy Efficiency Ratio drift correction value Δ E (t obtaining according to described step 3) ₁, t ₂) the Energy Efficiency Ratio Acquisition Error of tested room air conditioner is revised:

E _t=E _m-ΔE(t ₁,t ₂)；

Wherein, E _mfor the Energy Efficiency Ratio measured value under any operating mode of tested room air conditioner in wet and dry bulb tolerance; Δ E (t ₁, t ₂) be the drift correction value under this operating mode calculating in described step 3), E _tfor Energy Efficiency Ratio under this operating mode returns correction value.

The present invention revises the Energy Efficiency Ratio measured value under any operating mode in wet and dry bulb tolerance, obtain Energy Efficiency Ratio under this operating mode recurrence correction value to nominal condition, can be used for the actual operation parameters adjustment of convertible frequency air-conditioner: in the time that Energy Efficiency Ratio recurrence correction value is greater than Energy Efficiency Ratio measured value, suitably turn down compressor operating frequency; In the time that Energy Efficiency Ratio recurrence correction value is less than Energy Efficiency Ratio measured value, suitably heighten compressor operating frequency; In the time that Energy Efficiency Ratio recurrence correction value equals Energy Efficiency Ratio measured value, do not need to adjust.

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

ε=k σ, wherein ε represents to monitor residual error, and k represents to monitor coefficient, and σ represents the standard of Energy Efficiency Ratio actual measurement deviation

Deviation, $\mathrm{\σ}=\sqrt{\frac{\underset{i=1}{\overset{14}{\mathrm{\Σ}}}[\mathrm{\Δ}{E}_i-\frac{1}{14}\underset{i=1}{\overset{14}{\mathrm{\Σ}}}\mathrm{\Δ}{E}_i{]}^2}{13}},(i=\mathrm{1,2},...,14);$

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

ε _irepresent the residual error of Energy Efficiency Ratio measurement error and equation deviation under different operating points, ε _i=Δ E _i-Δ E _i(t ₁, t ₂); Δ E _irepresent the difference of Energy Efficiency Ratio measured value and nominal value; Δ E _i(t ₁, t ₂) represent the difference of Energy Efficiency Ratio functional value and nominal value; I represents different test operating modes, i=1, and 2,, 14.

In the inventive method, in described step 3), ${\mathrm{\θ}}_1=\frac{{\mathrm{\λ}}_1}{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|},{\mathrm{\θ}}_2=\frac{{\mathrm{\λ}}_2}{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|};$

${\mathrm{\α}}_{11}=e\sqrt{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|}/\left|\right|p_1\left|\right|,{\mathrm{\α}}_{12}=(c-d-\sqrt{(c-d{)}^2+e^2})\sqrt{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|}/\left|\right|p_1\left|\right|;$

${\mathrm{\α}}_{21}=e\sqrt{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|}/\left|\right|p_2\left|\right|,{\mathrm{\α}}_{22}=(c-d-\sqrt{(c-d{)}^2+e^2})\sqrt{\left|{\mathrm{\λ}}_1\right|+\left|{\mathrm{\λ}}_2\right|}/\left|\right|p_2\left|\right|;$

${\mathrm{\ζ}}_1=\frac{(a{\mathrm{\α}}_{22}-b{\mathrm{\α}}_{12})}{{\mathrm{\θ}}_1({\mathrm{\α}}_{11}{\mathrm{\α}}_{22}-{\mathrm{\α}}_{12}{\mathrm{\α}}_{21})},{\mathrm{\ζ}}_2=\frac{(b{\mathrm{\α}}_{11}-a{\mathrm{\α}}_{21})}{{\mathrm{\θ}}_2({\mathrm{\α}}_{11}{\mathrm{\α}}_{22}-{\mathrm{\α}}_{12}{\mathrm{\α}}_{21})};$

Wherein, λ ₁, λ ₂for the characteristic value of secondary temperature difference item coefficient matrix B, p ₁, p ₂for secondary temperature difference item coefficient matrix

The characteristic vector of B, $B=\left(\begin{array}{cc}c& e/2\\ e/2& d\end{array}\right),$ ${\mathrm{\λ}}_1=\frac{c+d-\sqrt{(c-d{)}^2+e^2}}{2},$ ${\mathrm{\λ}}_2=\frac{c+d+\sqrt{(c-d{)}^2+e^2}}{2},$

$\left|{\mathrm{\&lambda;}}_1\right|<\left|{\mathrm{\&lambda;}}_2\right|;\left|\right|p_1\left|\right|=(e^2+[c-d-\sqrt{(c-d{)}^2+e^2}{]}^2{)}^{\frac{1}{2}},\left|\right|p_2\left|\right|=(e^2+[c-d+\sqrt{(c-d{)}^2+e^2}{]}^2{)}^{\frac{1}{2}}).$

Beneficial effect: the present invention compared with prior art, has the following advantages:

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

The standard deviation of surveying deviation and function deviation take Energy Efficiency Ratio is as object, take the numerical relation of the residual sum standard deviation of test value and equation value as target, each test operating point is carried out to dynamic monitoring, with specific features and the accuracy thereof of determining that Energy Efficiency Ratio deviation changes with wet-bulb depression continuity.

The inventive method can obtain wet and dry bulb temperature in national standard tolerance arbitrarily under operating point Energy Efficiency Ratio to the recurrence correction value of nominal condition, thereby instruct the setting of the concrete operational factors of air-conditioning equipment such as convertible frequency air-conditioner, equipment operation is more accurate;

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

Energy Efficiency Ratio drift correction formula is once definite, and the room air conditioner Energy Efficiency Ratio measured value of follow-up identical Energy Efficiency Ratio grade is directly determined by this correction formula, avoids repeatedly repeating to test the huge consumption of the human and material resources that bring;

Along with instrument and meter precision improves constantly, avoid the actual test of manual control operating point to be positioned at wet and dry bulb temperature high tolerance or lower limit, thereby cause room air conditioner efficiency grade that predictable wrong division occurs.

Accompanying drawing explanation

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

The specific embodiment

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

In wet and dry bulb temperature franchise of the present invention, Energy Efficiency Ratio modification method can the Energy Efficiency Ratio test result under cooling condition be revised the room air conditioner of different Energy Efficiency Ratio grades.Be modified to embodiment with the Energy Efficiency Ratio measured result under three class room air conditioner cooling conditions below and describe, concrete steps are as follows:

Step 1) is captured in the room air conditioner Energy Efficiency Ratio data under operating condition of test, described operating condition of test is that dry-bulb temperature is 35 ℃ outside holding chamber, wet-bulb temperature is 24 ℃, inside difference control room, wet and dry bulb temperature is in 3 large classes, under 14 kinds of operating modes, carry out the measurement of room air conditioner Energy Efficiency Ratio, the deviation that records Energy Efficiency Ratio under Energy Efficiency Ratio and nominal condition is

ΔE ₁=-0.03267kW`kW <sup>-1</sup>、ΔE <sub>2</sub>=-0.01971kW`kW ^-1、ΔE ₃=0kW`kW ^-1、

ΔE ₄=0.0329kW`kW <sup>-1</sup>、ΔE <sub>5</sub>=0.05315kW`kW ^-1、ΔE ₆=-0.02574kW`kW ^-1、

ΔE ₇=-0.0094kW`kW <sup>-1</sup>、ΔE <sub>8</sub>=0kW`kW ^-1、ΔE ₉=0.03773kW`kW ^-1、

ΔE ₁₀=0.05315kW`kW <sup>-1</sup>、ΔE <sub>11</sub>=-0.06462kW`kW ^-1、ΔE ₁₂=0.01385kW`kW ^-1、

ΔE ₁₃=-0.05007kW`kW <sup>-1</sup>、ΔE <sub>14</sub>=0.10285kW`kW ^-1：

Wherein, first kind operating mode is that indoor wet-bulb temperature is constant, changes indoor dry-bulb temperature; Operating mode 1～5 indoor wet-bulb temperature is 19 ℃, corresponding 26 ℃, 26.5 ℃, 27 ℃, 27.5 ℃, 28 ℃ respectively of indoor dry-bulb temperatures;

Equations of The Second Kind operating mode is that indoor dry-bulb temperature is constant, changes indoor wet-bulb temperature; Operating mode 6～10 indoor dry-bulb temperatures are 27 ℃, corresponding 18.5 ℃, 18.8 ℃, 19 ℃, 19.3 ℃, 19.5 ℃ respectively of indoor wet-bulb temperature;

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

Step 2) the each floor data gathering in described step 1) is distinguished to the following regression equation of substitution:

ΔE(t ₁,t ₂)=E ₀+at ₁+bt ₂+ct ₁ ²+dt ₂ ²+et ₁t ₂，

Wherein, t ₁represent indoor dry-bulb temperature and the nominal condition point temperature difference of indoor dry-bulb temperature down under tested operating mode, t ₂represent indoor wet-bulb temperature and the nominal condition point temperature difference of indoor wet-bulb temperature down under tested operating mode;

E ₀represent that in regression equation, constant term, a represent temperature difference item t of indoor dry-bulb temperature ₁coefficient, b represent temperature difference item t of indoor wet-bulb temperature ₂coefficient, c be indoor dry-bulb temperature secondary temperature difference item t ₁ ²coefficient, d represent indoor wet-bulb temperature secondary temperature difference item

coefficient, e represent indoor wet and dry bulb temperature coupling temperature difference item t ₁t ₂coefficient;

Under the 3rd, 6 operating modes, t ₁, t ₂, t ₁ ², , t ₁t ₂, Δ E _i(t ₁, t ₂) value be 0, therefore obtain the overdetermined equation group of 12 hexa-atomic linear functions compositions:

E ₀-a+c=-0.02574；

E ₀-0.5a+0.25c=-0.01278；

E ₀+0.5a+0.25c=0.0329；

E ₀+a+c=0.05315；

E ₀-0.5b+0.25d=-0.02574；

E ₀-0.2b+0.04=-0.0094；

E ₀+0.3b+0.09d=0.03773；

E ₀+0.5b+0.25d=0.05315；

E ₀-a-0.5b+c+0.25d+0.5e=-0.06462；

E ₀-a+0.5b+c+0.25d-0.5e=0.01385；

E ₀+a-0.5b+c+0.25d-0.5e=-0.05007；

E ₀+a+0.5b+c+0.25d+0.5e=0.10285

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

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

The residual error of test value and equation value is ε _i=Δ E _i-Δ E _i(t ₁, t ₂), Δ E _irepresent the difference of Energy Efficiency Ratio measured value and nominal value, Δ E _i(t ₁, t ₂) represent the difference of Energy Efficiency Ratio functional value and nominal value; Trying to achieve according to above-mentioned formula the residual error that each operating point is corresponding is respectively: ε ₁=0.00245, ε ₂=0.00121, ε ₃=0.000059, ε ₄=0.00106, ε ₅=0.001304, ε ₆=0.00357, ε ₇=0.00372, ε ₈=0.00023, ε ₉=0.00094, ε ₁₀=0.00128, ε ₁₁=0.00275, ε ₁₂=0.0024,, ε ₁₃=0.00118 ε ₁₄=0.00132

σ represents the standard deviation of Energy Efficiency Ratio actual measurement deviation,

$\mathrm{\&sigma;}=\sqrt{\frac{\underset{i=1}{\overset{14}{\mathrm{\&Sigma;}}}[\mathrm{\&Delta;}{E}_i-\frac{1}{14}\underset{i=1}{\overset{14}{\mathrm{\&Sigma;}}}\mathrm{\&Delta;}{E}_i{]}^2}{13}}=0.002,(i=\mathrm{1,2},...,14);$

Bessel Formula is ε=k σ, respectively when monitoring coefficient k value 0.65,2,3 successively, in the time of k=0.65, has 57% Energy Efficiency Ratio actual measurement deviation to meet relational expression ε _i≤ ε, (i represents different test operating modes, i=1, and 2,, 14); In the time of k=2, there is 100% Energy Efficiency Ratio actual measurement deviation to meet relational expression ε _i≤ ε; In the time of k=3, there is 100% Energy Efficiency Ratio actual measurement deviation to meet relational expression ε _i≤ ε.Therefore the convergence of Energy Efficiency Ratio regression equation, W ₀, a, b, c, d, e value be respectively 0.0019,0.19,0.37,0.0027,0.011,0.0098;

ΔW(t ₁,t ₂)=W ₀+at ₁+bt ₂+ct ₁ ²+

+et ₁t ₂

=0.00 ²+0.044t ₁+0.088t ₂+0.0049t ₁ ²+0.019t ₂ ²+0.0069t ₁t ₂

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

ΔE(t ₁,t ₂)=E ₀+ζΛQt+(Qt) ^TΛQt；

Wherein, t is temperature difference vector, t=(t ₁, t ₂) ^t, T represents vectorial transposition;

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

Q is temperature transition matrix, $Q=\left(\begin{array}{cc}{\mathrm{\&alpha;}}_{11}& {\mathrm{\&alpha;}}_{12}\\ {\mathrm{\&alpha;}}_{21}& {\mathrm{\&alpha;}}_{22}\end{array}\right),$ α ₁₁represent the influence coefficient of the indoor dry bulb temperature difference to indoor dry-bulb temperature tolerance under test operating mode, α ₂₁represent the influence coefficient of the indoor dry bulb temperature difference to indoor wet bulb bulb temperature tolerance under test operating mode, α ₁₂represent the influence coefficient of indoor wet-bulb depression to indoor dry-bulb temperature tolerance under test operating mode, α ₂₂represent the influence coefficient of indoor wet-bulb depression to indoor wet-bulb temperature tolerance under test operating mode;

ζ be temperature difference item with respect to the high order temperature difference item relative effect coefficient vector to Energy Efficiency Ratio, ζ=(ζ ₁, ζ ₂); ζ ₁a temperature term of the expression indoor dry bulb temperature difference is the relative effect coefficient to Energy Efficiency Ratio deviation with respect to high order temperature term, ζ ₂a temperature term of expression indoor wet-bulb depression is the relative effect coefficient to Energy Efficiency Ratio deviation with respect to high order temperature term;

θ ₁, θ ₂, α ₁₁, α ₁₂, α ₂₁, α ₂₂, ζ ₁, ζ ₂determine by the following method:

Coefficient c, the d of secondary temperature term in Energy Efficiency Ratio regression equation, e are carried out to quadratic form and change to obtain Second Order with Constant Coefficients matrix $B=\left(\begin{array}{cc}c& e/2\\ e/2& d\end{array}\right),$ Determine its eigenvalue λ ₁, λ ₂and characteristic of correspondence vector p ₁, p ₂:

${\mathrm{\&lambda;}}_1=\frac{c+d+\sqrt{(c-d{)}^2+e^2}}{2},$ ${\mathrm{\&lambda;}}_2=\frac{c+d-\sqrt{(c-d{)}^2+e^2}}{2},$

$p_1=(e,c-d-\sqrt{(c-d{)}^2+e^2}{)}^T,p_2=(e,c-d+\sqrt{(c-d{)}^2+e^2}{)}^T;$

, ${\mathrm{\&theta;}}_1=\frac{{\mathrm{\&lambda;}}_1}{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|},{\mathrm{\&theta;}}_2=\frac{{\mathrm{\&lambda;}}_2}{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|};$ Weight matrix Λ=diag (θ ₁, θ ₂);

${\mathrm{\&alpha;}}_{11}=e\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_1\left|\right|,{\mathrm{\&alpha;}}_{12}=(c-d-\sqrt{(c-d{)}^2+e^2})\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_1\left|\right|,$

${\mathrm{\&alpha;}}_{21}=e\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_2\left|\right|,{\mathrm{\&alpha;}}_{22}=(c-d-\sqrt{(c-d{)}^2+e^2})\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_2\left|\right|,$ Temperature transition matrix $Q=\left(\begin{array}{cc}{\mathrm{\&alpha;}}_{11}& {\mathrm{\&alpha;}}_{12}\\ {\mathrm{\&alpha;}}_{21}& {\mathrm{\&alpha;}}_{22}\end{array}\right);{\mathrm{\&zeta;}}_1=\frac{(a{\mathrm{\&alpha;}}_{22}-b{\mathrm{\&alpha;}}_{12})}{{\mathrm{\&theta;}}_1({\mathrm{\&alpha;}}_{11}{\mathrm{\&alpha;}}_{22}-{\mathrm{\&alpha;}}_{12}{\mathrm{\&alpha;}}_{21})},{\mathrm{\&zeta;}}_2=\frac{(b{\mathrm{\&alpha;}}_{11}-a{\mathrm{\&alpha;}}_{21})}{{\mathrm{\&theta;}}_2({\mathrm{\&alpha;}}_{11}{\mathrm{\&alpha;}}_{22}-{\mathrm{\&alpha;}}_{12}{\mathrm{\&alpha;}}_{21})},$ Temperature difference item is the relative effect coefficient vector ζ=(ζ to Energy Efficiency Ratio with respect to high order temperature difference item ₁, ζ ₂);

Be that Energy Efficiency Ratio drift correction formula is:

$\mathrm{\&Delta;E}(t_1,t_2)=0.002+\left(\mathrm{0.84,0.77}\right)\left(\begin{array}{cc}0.16& \\ & 0.84\end{array}\right)\left(\begin{array}{cc}0.15& -0.036\\ 0.036& 0.15\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right)$

$+(\left(\begin{array}{cc}0.15& -0.036\\ 0.036& 0.15\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right){)}^T\left(\begin{array}{cc}0.16& \\ & 0.84\end{array}\right)\left(\begin{array}{cc}0.15& -0.036\\ 0.036& 0.15\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right)$

Calculate the Energy Efficiency Ratio drift correction amount under this operating mode according to above-mentioned solving result:

$\mathrm{\&Delta;E}(t_1,t_2)0.002+\left(\mathrm{0.84,0.77}\right)\left(\begin{array}{cc}0.16& \\ & 0.84\end{array}\right)\left(\begin{array}{cc}0.15& -0.036\\ 0.036& 0.15\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right)$

$+(\left(\begin{array}{cc}0.15& -0.036\\ 0.036& 0.15\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right){)}^T\left(\begin{array}{cc}0.16& \\ & 0.84\end{array}\right)\left(\begin{array}{cc}0.15& -0.036\\ 0.036& 0.15\end{array}\right)\left(\begin{array}{c}{t}_1\\ t_2\end{array}\right)=-0.01916k{W}^{\&prime;}k{W}^{-1};$

The Energy Efficiency Ratio drift correction value Δ W (t that step 4) obtains according to described step 3) ₁, t ₂) the Energy Efficiency Ratio Acquisition Error of tested room air conditioner is revised:

E _t=E _m-ΔE(t ₁,t ₂)=3.0953kW`kW <sup>-1</sup>+0.01916kW`kW ^-1=3.1145kW`kW ^-1

Be dry-bulb temperature be 26.7 ℃, wet-bulb temperature be at 18.9 ℃ Energy Efficiency Ratio to return correction value be 3.1145kW`kW ^-1.

Wherein, E _mfor the Energy Efficiency Ratio measured value under any operating mode of tested room air conditioner, be 3.0953kW`kW <sup>-1</sup>; Δ E (t <sub>1</sub>, t <sub>2</sub>) be the drift correction value under this operating mode calculating in described step 3), for-0.01916kW`kW ^-1; E _tfor Energy Efficiency Ratio under this operating mode returns correction value.

Claims (3)

1. an Energy Efficiency Ratio modification method in wet and dry bulb temperature franchise, is characterized in that, the method comprises the following steps:

1) be captured in the room air conditioner Energy Efficiency Ratio data under operating condition of test, described operating condition of test is that dry-bulb temperature is 35 ℃ outside holding chamber, wet-bulb temperature is 24 ℃, inside difference control room, wet and dry bulb temperature is in 3 large classes, under 14 kinds of operating modes, carry out the measurement of room air conditioner Energy Efficiency Ratio, the deviation that records Energy Efficiency Ratio under Energy Efficiency Ratio and nominal condition point is Δ E _i(i=1,2,, 14), i represents different test operating modes:

Wherein, first kind operating mode is that indoor wet-bulb temperature is constant, changes indoor dry-bulb temperature; Operating mode 1～5 indoor wet-bulb temperature is 19 ℃, corresponding 26 ℃, 26.5 ℃, 27 ℃, 27.5 ℃, 28 ℃ respectively of indoor dry-bulb temperatures;

Equations of The Second Kind operating mode is that indoor dry-bulb temperature is constant, changes indoor wet-bulb temperature; Operating mode 6～10 indoor dry-bulb temperatures are 27 ℃, corresponding 18.5 ℃, 18.8 ℃, 19 ℃, 19.3 ℃, 19.5 ℃ respectively of indoor wet-bulb temperature;

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

2) by the following regression equation of each floor data difference substitution gathering in described step 1), obtain the overdetermined equation group of 14 hexa-atomic linear function compositions, then take Bessel Formula as criterion solves described equation group, obtain every constant coefficient E of regression equation ₀, a, b, c, d, e value:

ΔE(t ₁,t ₂)=E ₀+at ₁+bt ₂+ct ₁ ²+dt ₂ ²+et ₁t ₂，

Wherein, t ₁represent indoor dry-bulb temperature and the nominal condition point temperature difference of indoor wet and dry bulb temperature down under tested operating mode, t ₂represent indoor wet-bulb temperature and the nominal condition point temperature difference of indoor wet-bulb temperature down under tested operating mode;

E ₀represent that in regression equation, constant term, a represent temperature difference item t of indoor dry-bulb temperature ₁coefficient, b represent temperature difference item t of indoor wet-bulb temperature ₂coefficient, c be indoor dry-bulb temperature secondary temperature difference item t ² ₁coefficient, d represent indoor wet-bulb temperature secondary temperature difference item t ² ₂coefficient, e represent indoor wet and dry bulb temperature coupling temperature difference item t ₁t ₂coefficient;

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

ΔE(t ₁,t ₂)=E ₀+ζΛQt+(Qt) ^TΛQt；

Wherein, t is temperature difference vector, t=(t ₁, t ₂) ^t, T represents vectorial transposition;

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

Q is temperature transition matrix, $Q=\left(\begin{array}{cc}{\mathrm{\&alpha;}}_{11}& {\mathrm{\&alpha;}}_{12}\\ {\mathrm{\&alpha;}}_{21}& {\mathrm{\&alpha;}}_{22}\end{array}\right),$ α ₁₁represent the influence coefficient of the indoor dry bulb temperature difference to indoor dry-bulb temperature tolerance under test operating mode, α ₂₁represent the influence coefficient of the indoor dry bulb temperature difference to indoor wet bulb bulb temperature tolerance under test operating mode, α ₁₂represent the influence coefficient of indoor wet-bulb depression to indoor dry-bulb temperature tolerance under test operating mode, α ₂₂represent the influence coefficient of indoor wet-bulb depression to indoor wet-bulb temperature tolerance under test operating mode;

ζ be temperature difference item with respect to the high order temperature difference item relative effect coefficient vector to Energy Efficiency Ratio, ζ=(ζ ₁, ζ ₂); ζ ₁a temperature term of the expression indoor dry bulb temperature difference is the relative effect coefficient to Energy Efficiency Ratio deviation with respect to high order temperature term, ζ ₂a temperature term of expression indoor wet-bulb depression is the relative effect coefficient to Energy Efficiency Ratio deviation with respect to high order temperature term;

4) the Energy Efficiency Ratio drift correction value Δ E (t obtaining according to described step 3) ₁, t ₂) the Energy Efficiency Ratio Acquisition Error of tested room air conditioner is revised:

E _t=E _m-ΔE(t ₁,t ₂)；

Wherein, E _mfor the Energy Efficiency Ratio measured value under any operating mode of tested room air conditioner in wet and dry bulb tolerance; Δ E (t ₁, t ₂) be the drift correction value under this operating mode calculating in described step 3), E _tfor Energy Efficiency Ratio under this operating mode returns correction value.

2. Energy Efficiency Ratio modification method in wet and dry bulb temperature franchise according to claim 1, is characterized in that described step 2) in Bessel Formula be: ε=k σ, wherein ε represents to monitor residual error, k represents to monitor coefficient, and σ represents the standard deviation of Energy Efficiency Ratio actual measurement deviation $\mathrm{\&sigma;}=\sqrt{\frac{\underset{i=1}{\overset{14}{\mathrm{\&Sigma;}}}[\mathrm{\&Delta;}{E}_i-\frac{1}{14}\underset{i=1}{\overset{14}{\mathrm{\&Sigma;}}}\mathrm{\&Delta;}{E}_i{]}^2}{13}},(i=\mathrm{1,2},...,14);$

Described Bessel Formula decision condition is: monitoring coefficient k value 0.65,2,3 successively, in the time of k=0.65, has at least 50% Energy Efficiency Ratio actual measurement deviation to meet relational expression ε _i≤ ε; In the time of k=2, have at least 95% Energy Efficiency Ratio actual measurement deviation to meet relational expression ε _i≤ ε; In the time of k=3, have at least 100% Energy Efficiency Ratio actual measurement deviation to meet relational expression ε _i≤ ε;

ε _irepresent the measurement error of Energy Efficiency Ratio under different operating modes and the residual error of equation deviation, ε _i=Δ E _i-Δ E _i(t ₁, t ₂); Δ E _irepresent the difference of Energy Efficiency Ratio measured value and nominal value; Δ E _i(t ₁, t ₂) represent the difference of Energy Efficiency Ratio functional value and nominal value; I represents different test operating modes, i=1,2,14.

3. Energy Efficiency Ratio modification method in wet and dry bulb temperature franchise according to claim 1, is characterized in that, in described step 3), ${\mathrm{\&theta;}}_1=\frac{{\mathrm{\&lambda;}}_1}{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|},{\mathrm{\&theta;}}_2=\frac{{\mathrm{\&lambda;}}_2}{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|};$

${\mathrm{\&alpha;}}_{11}=e\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_1\left|\right|,{\mathrm{\&alpha;}}_{12}=(c-d-\sqrt{(c-d{)}^2+e^2})\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_1\left|\right|;$

${\mathrm{\&alpha;}}_{21}=e\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_2\left|\right|,{\mathrm{\&alpha;}}_{22}=(c-d-\sqrt{(c-d{)}^2+e^2})\sqrt{\left|{\mathrm{\&lambda;}}_1\right|+\left|{\mathrm{\&lambda;}}_2\right|}/\left|\right|p_2\left|\right|;$

${\mathrm{\&zeta;}}_1=\frac{(a{\mathrm{\&alpha;}}_{22}-b{\mathrm{\&alpha;}}_{12})}{{\mathrm{\&theta;}}_1({\mathrm{\&alpha;}}_{11}{\mathrm{\&alpha;}}_{22}-{\mathrm{\&alpha;}}_{12}{\mathrm{\&alpha;}}_{21})},{\mathrm{\&zeta;}}_2=\frac{(b{\mathrm{\&alpha;}}_{11}-a{\mathrm{\&alpha;}}_{21})}{{\mathrm{\&theta;}}_2({\mathrm{\&alpha;}}_{11}{\mathrm{\&alpha;}}_{22}-{\mathrm{\&alpha;}}_{12}{\mathrm{\&alpha;}}_{21})};$

Wherein, λ ₁, λ ₂for the characteristic value of secondary temperature difference item coefficient matrix B, p ₁, p ₂for the characteristic vector of secondary temperature difference item coefficient matrix B, $B=\left(\begin{array}{cc}c& e/2\\ e/2& d\end{array}\right),$ ${\mathrm{\&lambda;}}_1=\frac{c+d-\sqrt{(c-d{)}^2+e^2}}{2},$ ${\mathrm{\&lambda;}}_2=\frac{c+d+\sqrt{(c-d{)}^2+e^2}}{2},$

$\left|{\mathrm{\&lambda;}}_1\right|<\left|{\mathrm{\&lambda;}}_2\right|;\left|\right|p_1\left|\right|=(e^2+[c-d-\sqrt{(c-d{)}^2+e^2}{]}^2{)}^{\frac{1}{2}},\left|\right|p_2\left|\right|=(e^2+[c-d+\sqrt{(c-d{)}^2+e^2}{]}^2{)}^{\frac{1}{2}}).$

Images