Optimal operation control method of air-source heat pump and gas-fired heater combined heating system

Abstract

The present disclosure discloses an optimal operation control method of the air-source heat pump and gas-fired heater combined heating system, including constructing a mathematical model of the air-source heat pump and gas-fired heater combined heating system to simulate real time energy consumption; determining a comprehensive evaluation index system of the combined heating system, including primary evaluation indexes of energy conservation, environmental protection and economical efficiency, and secondary evaluation indexes of EER, clean energy utilization rate, carbon dioxide emission and operation cost; calculating indexes through an analytic hierarchy process and then constructing a comprehensive objective function, to determine operation mode when the comprehensive objective function is the maximum by taking the priority of meeting the heating requirement as a principle, so as to achieve the purpose of optimizing the combined heating system, and achieve the combined heating system efficient, energy-saving and environmental protection.

Claims

1. An optimal operation control method of the air-source heat pump and gas-fired heater combined heating system, including the following steps: step 1. constructing a mathematical model of the air-source heat pump and gas-fired heater combined heating system, and calculating hourly operation performance of the combined heating system; wherein the mathematical model includes a heating and energy consumption model for an air-source heat pump unit and a heating and energy consumption model for a gas-fired heater unit; step 2. determining a comprehensive evaluation index system of the combined heating system, including primary evaluation indexes and secondary evaluation indexes, and constructing a secondary evaluation index calculation model; wherein 3 primary evaluation indexes are provided, including energy conservation, environmental protection and economical efficiency, and 4 secondary evaluation indexes are provided, including a comprehensive energy efficiency ratio (EER) and a clean energy utilization rate, , that are subordinate to the energy conservation index, a CO2 emission TCO2 that is subordinate to the environment protection index, and operation cost Cr that is subordinate to the economical efficiency index; a calculation formula of the EER of the combined heating system is as follows: $\mathrm{EER}=\frac{Q_p+Q_b}{W_p/{}_e+V_g.Math.H_i}$ a calculation formula of the clean energy utilization rate, , of the combined heating system is as follows: $=\frac{Q_p}{Q_p+Q_b}$ a calculation formula of the CO.sub.2 emission, T.sub.CO2, of the combined heating system is as follows: T.sub.CO2=T.sub.b+T.sub.p, T.sub.b=V.sub.ga.sub.b, T.sub.p=W.sub.pa.sub.p; a calculation formula of the operation cost, C.sub.r, of the combined heating system is as follows:
C.sub.r=W.sub.pP.sub.e/(3.610.sup.6)+V.sub.gP.sub.g where W.sub.p is the power consumption of the air-source heat pump unit; .sub.e is the power generation efficiency of a gas-fired device; T.sub.CO2 is the CO.sub.2 emission of the combined heating system; T.sub.b is the CO.sub.2 emission of the gas-fired heater unit; T.sub.p is the CO.sub.2 emission of the air-source heat pump unit; a.sub.b is a conversion coefficient of the gas consumption to the carbon dioxide emission of the gas-fired heater for heating; .sub.p is a conversion coefficient of the power generation capacity to the carbon dioxide emission; C.sub.r is the operation cost; P.sub.e is an electricity price; and P.sub.g is a fired gas price; Q.sub.p is heating capacity of the air-source heat pump unit; Q.sub.b is heating capacity of the gas-fired heater unit; V.sub.g is the gas consumption per second of the gas-fired heater unit; H.sub.i is a calorific value of fired gas fed into the gas-fired heater for heating; step 3. constructing an index determination matrix through an analytic hierarchy process based on the evaluation index system of the combined heating system, carrying out normalization and consistency check, and determining respective weights of the primary evaluation indexes and the secondary evaluation indexes; step 4. normalizing the secondary evaluation indexes of the combined heating system, and respectively calculating membership functions of the EER, the clean energy utilization rate, the carbon dioxide emission and the operation cost; step 5. a comprehensive objective function calculation model is constructed based on the membership functions of the secondary evaluation indexes, with a calculation formula as follows: =w.sub.1(w.sub.11.sub.1(EER)+w.sub.12.sub.2())+w.sub.2.sub.3 (T.sub.CO2)+w.sub.3.sub.4 (C.sub.r), where w.sub.1, w.sub.2, and w.sub.3 denote the weights of the primary evaluation indexes: the energy conservation, the environmental protection and the economical efficiency, respectively; and w.sub.11 and w.sub.12 denote the weights of the secondary evaluation indexes: the EER and the clean energy utilization rate, respectively; step 6. implementing floor radiant heating at a heating terminal form of the combined heating system, constructing a heating parameter prediction model based on characteristics of the heating terminal form, and performing calculation in a constant flow control method to obtain an actual operation water supply temperature t.sub.g, an actual operation return water temperature t.sub.h and a heating load Q at different outdoor temperatures; and step 7. determining whether the operation of an air-source heat pump alone can meet the requirement of the heating load at the current outdoor temperature in real time by taking the priority of meeting the heating requirement as a principle; if the air-source heat pump cannot meet the requirement, respectively calculating a comprehensive objective function in a combined operation mode of the air-source heat pump and a gas-fired heater for heating and a comprehensive objective function in an operation mode of the gas-fired heater for heating alone, whichever the comprehensive objective function is greater; and if the air-source heat pump can meet the requirement of the heating load, calculating a comprehensive objective function in an operation mode of the air-source heat pump alone and the comprehensive objective function in the operation mode of the gas-fired heater for heating alone, respectively, whichever the comprehensive objective function is greater, and based on the operation mode with the greater comprehensive objective function, obtaining an operation control strategy of the air-source heat pump and gas-fired heater combined heating system wherein the operation control strategy is employed to control and operate the air-source heat pump and gas-fired heater.

2. The optimal operation control method of the air-source heat pump and gas-fired heater combined heating system according to claim 1, wherein a calculation formula of a heating capacity of the air-source heat pump unit in step 1 is as follows: $Q_p=\{\begin{array}{c}\begin{array}{c}{Q}_{\mathrm{rate}}\mathrm{when}\mathrm{the}\mathrm{air}-\mathrm{source}\mathrm{heat}\mathrm{pump}\mathrm{is}\mathrm{under}a\\ \mathrm{non}-\mathrm{frosting}\mathrm{condition}\end{array}\\ \begin{array}{c}{Q}_{\mathrm{real}}\mathrm{when}\mathrm{the}\mathrm{air}-\mathrm{source}\mathrm{heat}\mathrm{pump}\mathrm{is}\mathrm{under}a\\ \mathrm{frosting}\mathrm{condition}\mathrm{and}\mathrm{the}\mathrm{defrosting}\mathrm{condition}\end{array}\end{array}$ a calculation formula of power consumption of the air-source heat pump unit is as follows: $W_p=\{\begin{array}{c}\begin{array}{c}{W}_{\mathrm{rate}}\mathrm{when}\mathrm{the}\mathrm{air}-\mathrm{source}\mathrm{heat}\mathrm{pump}\mathrm{is}\mathrm{under}a\\ \mathrm{non}-\mathrm{frosting}\mathrm{condition}\end{array}\\ \begin{array}{c}{W}_{\mathrm{real}}\mathrm{when}\mathrm{the}\mathrm{air}-\mathrm{source}\mathrm{heat}\mathrm{pump}\mathrm{is}\mathrm{under}a\\ \mathrm{frosting}\mathrm{condition}\mathrm{and}\mathrm{the}\mathrm{defrosting}\mathrm{condition}\end{array}\end{array}$ where Q.sub.p is the heating quantity of the air-source heat pump unit; Q.sub.rate is the heating quantity of the air-source heat pump unit under the non-frosting condition; Q.sub.real is the heating capacity of the air-source heat pump unit under the frosting condition and the defrosting condition; W.sub.p is the power consumption of the air-source heat pump unit; W.sub.rate is the power consumption of the air-source heat pump unit under the non-frosting condition; and W.sub.real is the power consumption of the air-source heat pump unit under the frosting condition; data fitting is carried out according to an air-source heat pump heating performance curve provided by a manufacturer to obtain the heating capacity Q.sub.rate, coefficient of performance COP.sub.rate and power consumption W.sub.rate of the air-source heat pump unit under the non-frosting condition, wherein calculation formulas are as follows: $Q_{rate}=a_0+a_1{t}_a+a_2{t}_a^2+a_3{t}_a^3,$ ${\mathrm{COP}}_{rate}=a_4+a_5{t}_a+a_6{t}_a^2+a_7{t}_a^2,$ $W_{rate}=\frac{Q_{\mathrm{rte}}}{CO{P}_{\mathrm{rte}}};$ where Q.sub.rate is the heating capacity of the air-source heat pump unit under the non-frosting condition; COP.sub.rate is the coefficient of performance of the air-source heat pump unit under the non-frosting condition; W.sub.rate is the power consumption of the air-source heat pump unit under the non-frosting condition; t.sub.a is an ambient temperature; and a.sub.0, a.sub.1, a.sub.2, a.sub.3, a.sub.4, a.sub.5, a.sub.6, and a.sub.7 are coefficients of a fitting equation; according to a heat pump heating capacity correction model under the frosting and defrosting conditions, heating capacity Q.sub.real of the air-source heat pump unit under the frosting condition and the defrosting condition is obtained, wherein a calculation formula is as follows:
Q.sub.real=Q.sub.rate[0.311t.sub.a+0.043t.sub.a.sup.2+0.005t.sub.a.sup.3(0.7831.07210.sup.4t.sub.a.sup.3)RH.sup.0.8461.647] where RH is a relative ambient humidity; according to a heat pump coefficient-of-performance correction model under the frosting condition, when the outdoor temperature is higher and not higher than 7 C., performance coefficient COP.sub.real of the air-source heat pump unit under the frosting condition and the defrosting condition is as follows: ${\mathrm{COP}}_{real}=\{\begin{array}{c}{\mathrm{COP}}_{\mathrm{rate}}\left(1-0.1801e^{-\frac{t_a^2}{5}}\right),t_a>7\\ {\mathrm{COP}}_{\mathrm{rate}}\left[1+0.0027\left(t_a-7\right)-0.1801e^{-\frac{t_a^2}{5}}\right],t_{a}7\end{array};$ and power consumption W.sub.real of the air-source heat pump unit under the frosting condition and the defrosting condition is as follows: $W_{real}=\frac{Q_{real}}{CO{P}_{real}};$ a calculation formula of the gas consumption per second V.sub.g of the heating and energy consumption model for the gas-fired heater unit, namely the gas-fired heater unit, is as follows: V.sub.g=Q.sub.b/(H.sub.i.Math..sub.b); wherein Q.sub.b is the heating capacity of the gas-fired heater unit; H.sub.i is a calorific value of fired gas fed into the gas-fired heater; and .sub.b is the heat efficiency of the gas-fired heater for heating; wherein a calculation formula of Q.sub.b is as follows: $Q_b=\{\begin{array}{c}\begin{array}{c}Q-Q_p\mathrm{in}\mathrm{the}\mathrm{event}\mathrm{of}\mathrm{the}\mathrm{combined}\mathrm{operation}\mathrm{of}\mathrm{the}\\ \mathrm{gas}-\mathrm{fired}\mathrm{heater}\mathrm{and}\mathrm{the}\mathrm{air}-\mathrm{source}\mathrm{heat}\mathrm{pump}\end{array}\\ Q\mathrm{in}\mathrm{the}\mathrm{event}\mathrm{of}\mathrm{the}\mathrm{operation}\mathrm{of}\mathrm{the}\mathrm{gas}-\mathrm{fired}\mathrm{heater}\mathrm{alone}\end{array}$ where Q is an actual heating load; and Q.sub.p is the heating capacity of the air-source heat pump unit; data fitting is carried out according to a load-heat efficiency operation curve provided by a manufacturer of the gas-fired heater for heating to obtain the operation heat efficiency .sub.b of the gas-fired heater unit at different load rates, wherein calculation formulas are as follows: ${}_n=a_8+a_9{}_b+a_{10}{}_b^2+a_{11}{}_b^3,{}_b=\frac{Q_b}{Q_{b0}};$ where a.sub.8, a.sub.9, a.sub.10, and a.sub.11 are coefficients of a fitting equation; .sub.b is a load rate of the gas-fired heater unit; Q.sub.b is the heating capacity of the gas-fired heater unit; and Q.sub.b0 is the rated heating capacity of the gas-fired heater unit, wherein the calculated values employed to control and operate the heat pump and gas-fired heater.

3. The optimal operation control method of the air-source heat pump and gas-fired heater combined heating system according to claim 1, wherein calculation formulas of membership functions of EER, clean energy utilization rate, carbon dioxide emission and operation cost are as follows: $f_1\left(\mathrm{EER}\right)=\{\begin{array}{cc}0& \mathrm{EER}<{\mathrm{EER}}_{\min }\\ \frac{\mathrm{EER}-{\mathrm{EER}}_{\min }}{{\mathrm{EER}}_{\max }-{\mathrm{EER}}_{\min }}& {\mathrm{EER}}_{\min }\mathrm{EER}<{\mathrm{EER}}_{\max }\\ 1& \mathrm{EER}{\mathrm{EER}}_{\max }\end{array};$ $f_2\left(\right)=\{\begin{array}{cc}0& <{}_{\min }\\ \frac{-{}_{\min }}{{}_{\max }-{}_{\min }}& {}_{\min }<{}_{\max }\\ 1& {}_{\max }\end{array};$ $f_3\left(T_{\mathrm{CO}2}\right)=\{\begin{array}{cc}1& T_{\mathrm{CO}2}{T}_{\mathrm{CO}{2}_{\max }}\\ \frac{T_{\mathrm{CO}{2}_{\max }}-T_{\mathrm{CO}2}}{T_{\mathrm{CO}{2}_{\max }}-T_{\mathrm{CO}{2}_{\min }}}& T_{\mathrm{CO}{2}_{\min }}{T}_{\mathrm{CO}2}<T_{\mathrm{CO}{2}_{\max }}\\ 0& T_{\mathrm{CO}2}<T_{\mathrm{CO}{2}_{\min }}\end{array};$ $f_4\left(C_r\right)=\{\begin{array}{cc}1& C_r{C}_{r_{\max }}\\ \frac{C_{r_{\max }}-C_r}{C_{r_{\max }}-C_{r_{\min }}}& C_{r_{\min }}{C}_r<C_{r_{\max }}\\ 0& C_r<C_{r_{\min }}\end{array};$ where EER.sub.max and EER.sub.min are EERs of the combined heating system in the optimal operation state and the worst operation state, respectively; .sub.max and .sub.min represent clean energy utilization rates of the combined heating system in the optimal operation state and the worst operation state, respectively; T.sub.CO2 max and T.sub.CO2 min are carbon dioxide emissions of the combined heating system in the worst operation state and the optimal operation state, respectively; and C.sub.r max and C.sub.r min are operation costs of the combined heating system in the worst operation state and the optimal operation state, respectively; when .sub.1 (EER) and .sub.2 () are 1, the combined heating system is in the optimal operation state, and calculated values of the EER and the clean energy utilization rate are the maximum; and when .sub.1 (EER) and .sub.2 () are 0, the combined heating system is in the worst operation state, and calculated values of the EER and the clean energy utilization rate are the minimum; when .sub.3 (T.sub.CO2) and .sub.4 (C.sub.r) are 1, the combined heating system is in the worst operation state, and calculated values of the carbon dioxide emission T.sub.CO2 and the operation cost C.sub.r are the maximum; and when .sub.3 (T.sub.CO2) and .sub.4 (C.sub.r) are 0, the combined heating system is in the optimal operation state, and calculated values of the carbon dioxide emission T.sub.CO2 and the operation cost C.sub.r are the minimum, wherein the calculated value employed to control and operate the heat pump and gas-fired heater.

4. The optimal operation control method of the air-source heat pump and gas-fired heater combined heating system according to claim 1, wherein the formulas in step 6 are as follows: $\{\begin{array}{c}{t}_g=t_n+0.5\left(t_g^{}+t_h^{}-2t_n\right){\stackrel{\_}{Q}}^{\frac{1}{1+C_1}}+0.5\left(t_g^{}+t_h^{}\right)\stackrel{\_}{Q}\\ t_h=t_n+0.5\left(t_g^{}+t_h^{}-2t_n\right){\stackrel{\_}{Q}}^{\frac{1}{1+C_1}}-0.5\left(t_g^{}+t_h^{}\right)\stackrel{\_}{Q}\end{array},$ $Q=Q_{\mathrm{ratio}}{Q}^{},Q_{\mathrm{ratio}}=\frac{t_n-t_a}{t_n-t_a^{}};$ where t.sub.g is a designed water supply temperature; t.sub.h is a designed return water temperature; t.sub.g is an actual operation water supply temperature; t.sub.h is an actual operation return water temperature; c.sub.1 is a characteristic parameter of a radiant panel in floor radiant heating; Q.sub.ratio is a relative heating load ratio; Q is a designed heating load; t.sub.n is a calculated indoor heating temperature; and t.sub.a is a calculated outdoor heating temperature, wherein the calculated value employed to control and operate the air-source heat pump and gas-fired heater.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is a flow block diagram illustrating an optimal operation control method of an air-source heat pump and gas-fired heater combined heating system according to the present disclosure;

(2) FIG. 2 is a flow block diagram illustrating determination of an operation mode of a combined heating system;

(3) FIG. 3 is a diagram illustrating hourly meteorological data and water supply temperature in Embodiment 1;

(4) FIG. 4 is a diagram illustrating an hourly output of a heat source device in Embodiment 1; and

(5) FIG. 5 is a diagram illustrating hourly evaluation indexes and comprehensive objective function values of a combined heating system in Embodiment 1.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(6) The present disclosure will be further described below with reference to drawings and embodiments. The specific embodiments described herein are only for explaining the present disclosure, but are not for limiting the scope of protection of the present disclosure.

(7) An optimal operation control method of the present disclosure mainly solves the optimal control of a combined heating system, provides comprehensive evaluation indexes of the combined heating system, determines a weight of each evaluation index through an analytic hierarchy process, defines comprehensive objective functions, and determines the operation control method of the combined heating system with the maximum comprehensive objective function as an optimizing purpose.

(8) FIG. 1 is a flow block diagram illustrating an optimal operation control method of an air-source heat pump and gas-fired heater combined heating system according to the present disclosure. As shown in FIG. 1, the optimal operation control method of the air-source heat pump and gas-fired heater combined heating system according to the present disclosure includes the following steps:

(9) Step 1. a mathematical model of the air-source heat pump and gas-fired heater combined heating system is constructed, and hourly operation performance of the combined heating system is calculated, wherein the mathematical model includes a heating capacity and energy consumption model for an air-source heat pump unit, and a heating capacity and energy consumption model for a gas-fired heater unit; a calculation formula of a heating capacity of the air-source heat pump unit is as follows:

(10) $\begin{array}{cc}{Q}_p=\{\begin{array}{c}\begin{array}{c}{Q}_{\mathrm{rate}}\mathrm{when}\mathrm{the}\mathrm{air}\mathrm{source}\mathrm{heat}\mathrm{pump}\mathrm{is}\mathrm{under}a\\ \mathrm{non}\mathrm{frosting}\mathrm{condition}\end{array}\\ \begin{array}{c}{Q}_{\mathrm{real}}\mathrm{when}\mathrm{the}\mathrm{air}\mathrm{source}\mathrm{heat}\mathrm{pump}\mathrm{is}\mathrm{under}a\\ \mathrm{frosting}\mathrm{condition}\mathrm{and}\mathrm{the}\mathrm{defrosting}\mathrm{condittion}\end{array}\end{array}& \left(1\right)\end{array}$ a calculation formula of power consumption of the air-source heat pump unit is as follows:

(11) $\begin{array}{cc}{W}_p=\{\begin{array}{c}\begin{array}{c}{W}_{\mathrm{rate}}\mathrm{when}\mathrm{the}\mathrm{air}\mathrm{source}\mathrm{heat}\mathrm{pump}\mathrm{is}\mathrm{under}\mathrm{the}\\ \mathrm{non}\mathrm{frosting}\mathrm{condition}\end{array}\\ \begin{array}{c}{W}_{\mathrm{real}}\mathrm{when}\mathrm{the}\mathrm{air}\mathrm{source}\mathrm{heat}\mathrm{pump}\mathrm{is}\mathrm{under}a\\ \mathrm{frosting}\mathrm{condition}\mathrm{and}\mathrm{the}\mathrm{defrosting}\mathrm{condittion}\end{array}\end{array}& \left(2\right)\end{array}$

(12) Where Q.sub.p is the heating capacity of the air-source heat pump unit, W; Q.sub.rate is the heating capacity of the air-source heat pump unit under the non-frosting condition, W; Q.sub.real is the heating capacity of the air-source heat pump unit under the frosting condition and the defrosting condition, W; W.sub.p is the power consumption of the air-source heat pump unit, W; W.sub.rate is the power consumption of the air-source heat pump unit under the non-frosting condition, W; and W.sub.real is the power consumption of the air-source heat pump unit under the frosting condition, W.

(13) Data fitting is carried out according to the air-source heat pump heating performance curve provided by a manufacturer to obtain the heating capacity Q.sub.rate, coefficient of performance COP.sub.rate and power consumption W.sub.rate of the air-source heat pump unit under the non-frosting condition, wherein calculation formulas are as follows:

(14) $\begin{array}{cc}{Q}_{\mathrm{rate}}=a_0+a_1{t}_a+a_2{t}_a^2+a_3{t}_a^3& \left(3\right)\end{array}$$\begin{array}{cc}{\mathrm{COP}}_{\mathrm{rate}}=a_4+a_5{t}_a+a_6{t}_a^2+a_7{t}_a^2& \left(4\right)\end{array}$$\begin{array}{cc}{W}_{\mathrm{rate}}=\frac{Q_{rate}}{CO{P}_{rate}}& \left(5\right)\end{array}$

(15) Where Q.sub.rate is the heating capacity of the air-source heat pump unit under the non-frosting condition, W; COP.sub.rate is the coefficient of performance of the air-source heat pump unit under the non-frosting condition; W.sub.rate is the power consumption of the air-source heat pump unit under the non-frosting condition, W; t.sub.a is an ambient temperature, C.; and a.sub.0, a.sub.1, a.sub.2, a.sub.3, a.sub.4, a.sub.5, a.sub.6, and a.sub.7 are coefficients of a fitting equation.

(16) According to a heat pump heating capacity correction model under the frosting and defrosting conditions proposed by the related reference (Energy & Buildings, An equivalent temperature drop method for evaluating the operating performances of ASHP units jointly affected by ambient air temperature and relative humidity, Volume 224, Oct. 1, 2020, 110211), the heating capacity, Q.sub.real, of the air-source heat pump unit under the frosting condition and the defrosting condition is obtained, wherein a calculation formula is as follows:
Q.sub.real=Q.sub.rate[0.311t.sub.a+0.043t.sub.a.sup.2+0.005t.sub.a.sup.3(0.7831.07210.sup.4t.sub.a.sup.3)RH.sup.0.8461.647](6)

(17) Where Q.sub.real is the heating capacity of the air-source heat pump unit under the frosting condition and the defrosting condition, W; Q.sub.rate is the heating capacity of the air-source heat pump unit under the non-frosting condition, W; t.sub.a is the ambient temperature, C.; and RH is a relative ambient humidity, %.

(18) Similarly, according to a heat pump coefficient-of-performance correction model under the frosting condition proposed by the related reference (Energy & Buildings, Towards low carbon homes-A simulation analysis of building-integrated air-source heat pump systems, Volume 48, May 1, 2012, pages 127-136), the coefficient of performance, COP.sub.real, of the air-source heat pump unit under the frosting condition and the defrosting condition is obtained, respectively, when the outdoor temperature is higher and not higher than 7 C., wherein a calculation formula is as follows:

(19) $\begin{array}{cc}{\mathrm{COP}}_{real}=\{\begin{array}{c}{\mathrm{COP}}_{\mathrm{rate}}\left(1-0.1801e^{-\frac{t_a^2}{5}}\right),t_a>7\\ {\mathrm{COP}}_{\mathrm{rate}}\left[1-0.0027\left(t_a-7\right)-0.1801e^{-\frac{t_a^2}{5}}\right],t_{a}7\end{array}& \left(7\right)\end{array}$

(20) Where COP.sub.real is the coefficient of performance of the air-source heat pump unit under the frosting condition and the defrosting condition; COP.sub.rate is the coefficient of performance of the air-source heat pump unit under the non-frosting condition; and t.sub.a is the ambient temperature, C.

(21) Based on this, the power consumption, W.sub.real, of the air-source heat pump unit under the frosting condition and the defrosting condition is obtained, wherein a calculation formula is as follows:

(22) $\begin{array}{cc}{W}_{real}=\frac{Q_{real}}{CO{P}_{real}}& \left(8\right)\end{array}$

(23) Where W.sub.real is the power consumption of the air-source heat pump unit under the frosting condition and the defrosting condition, W; Q.sub.real is the heating capacity of the air-source heat pump unit under the frosting condition and the defrosting condition, W; and COP.sub.real is the coefficient of performance of the air-source heat pump unit under the frosting condition and the defrosting condition.

(24) A calculation formula of the gas consumption per second, V.sub.g, of the heating capacity and energy consumption model for the gas-fired heater unit, is as follows:
V.sub.g=Q.sub.b/(H.sub.i.Math..sub.b)(9)

(25) Where V.sub.g is the gas consumption per second, m.sup.3/s, of the gas-fired heater unit; Q.sub.b is the heating capacity of the gas-fired heater unit, W; H.sub.i is a calorific value, J/m.sup.3, of fired gas fed into the gas-fired heater; and .sub.b is the heat efficiency of the gas-fired heater, %.

(26) Wherein a calculation formula of Q.sub.b is as follows:

(27) $\begin{array}{cc}{Q}_b=\{\begin{array}{c}\begin{array}{c}Q-Q_p\mathrm{in}\mathrm{the}\mathrm{event}\mathrm{of}\mathrm{combined}\mathrm{operation}\mathrm{of}\mathrm{the}\\ \mathrm{gas}\mathrm{fired}\mathrm{heater}\mathrm{and}\mathrm{the}\mathrm{air}\mathrm{source}\mathrm{heat}\mathrm{pump}\end{array}\\ Q\mathrm{in}\mathrm{the}\mathrm{event}\mathrm{of}\mathrm{the}\mathrm{operation}\mathrm{of}\mathrm{the}\mathrm{gas}\mathrm{fired}\mathrm{heater}\mathrm{alone}\end{array}& \left(10\right)\end{array}$

(28) Where Q is an actual heating load, W; and Q.sub.p is the heating capacity of the air-source heat pump unit, W.

(29) Data fitting is carried out according to a load-heat efficiency operation curve provided by a manufacturer of the gas-fired heater to obtain the operation heat efficiency (.sub.b) of the gas-fired heater unit at different load rates, wherein calculation formulas are as follows:

(30) $\begin{array}{cc}{}_b=a_8+a_9{}_b+a_{10}{}_b^2+a_{11}{}_b^3& \left(11\right)\end{array}$$\begin{array}{cc}{}_b=\frac{Q_b}{Q_{b0}}& \left(12\right)\end{array}$

(31) Where .sub.b is the operation heat efficiency of the gas-fired heater unit, %,; a.sub.8, a.sub.9, a.sub.10, and a.sub.11 are coefficients of a fitting equation; .sub.b is a load rate of the gas-fired heater unit; Q.sub.b is the heating capacity of the gas-fired heater unit, W; and Q.sub.b0 is the rated heating capacity of the gas-fired heater unit, W.

(32) Step 2. a comprehensive evaluation index system of the combined heating system is determined, wherein the comprehensive evaluation index system includes primary evaluation indexes and secondary evaluation indexes, and a secondary evaluation index calculation model is constructed; and

(33) 3 primary evaluation indexes are provided, including energy conservation, environmental protection and economical efficiency, and 4 secondary evaluation indexes are provided, including a comprehensive energy efficiency ratio (EER) and a clean energy utilization rate that are subordinate to the energy conservation index, a carbon dioxide emission T.sub.CO2 that is subordinate to the environment protection index, and operation cost C.sub.r that is subordinate to the economical efficiency index.

(34) A calculation formula of the EER of the combined heating system is as follows:

(35) $\begin{array}{cc}\mathrm{EER}=\frac{Q_p+Q_b}{W_P/{}_e+V_g.Math.H_i}& \left(13\right)\end{array}$

(36) A calculation formula of the clean energy utilization rate of the combined heating system is as follows:

(37) $\begin{array}{cc}=\frac{Q_p}{Q_p+Q_b}& \left(14\right)\end{array}$

(38) A calculation formula of the carbon dioxide emission (T.sub.CO2) of the combined heating system is as follows:
T.sub.CO2=T.sub.b+T.sub.p(15)
T.sub.b=V.sub.ga.sub.b(16)
T.sub.p=W.sub.pa.sub.p(17)

(39) A calculation formula of the operation cost C.sub.r of the combined heating system is as follows:
C.sub.r=W.sub.pP.sub.e/(3.610.sup.6)+V.sub.gP.sub.g(18)

(40) Where W.sub.p is the power consumption of the air-source heat pump unit, W; .sub.e is the power generation efficiency of a gas-fired device, which is 40%; T.sub.CO2 is the carbon dioxide emission of the combined heating system, kg/s; T.sub.b is the carbon dioxide emission of the gas-fired heater unit, kg/s; T.sub.p is the carbon dioxide emission of the air-source heat pump unit, kg/s; a.sub.b is a conversion coefficient of the gas consumption to the carbon dioxide emission of the gas-fired heater, which is 1.964 kg/m.sup.3; a.sub.p is a conversion coefficient of the power generation capacity to the carbon dioxide emission, which is 0.728 kg/kWh; C.sub.r is the operation cost, Yuan/s; P.sub.e is an electricity price, Yuan/kWh; and P.sub.g is a fired gas price, Yuan/m.sup.3.

(41) Step 3. an index determination matrix is constructed through an analytic hierarchy process based on the evaluation index system of the combined heating system, normalization and consistency check are carried out, and respective weights, W.sub.i, of the primary evaluation indexes and the secondary evaluation indexes are determined.

(42) The following examples will illustrate the specific steps of the present disclosure:

(43) Firstly, an expert questionnaire survey is conducted to determine and analyze the importance of the primary evaluation indexes, which is divided into nine levels, pairwise comparisons among the primary evaluation indexes are then performed, and based on the importance from each pairwise comparison among the primary indexes, a determination matrix for the primary evaluation indexes for evaluation of the combined heating system is constructed, where levels 1, 3, 5, 7, and 9 represent that the two primary evaluation indexes are equally important, slightly important, quite important, obviously important and absolutely important, respectively, and levels 2, 4, 6, and 8 represent importance between two adjacent importance levels of the indexes. The indexes that are considered relatively less important are assigned reciprocal values, such as , , 1/7, and 1/9. A determination matrix for the secondary indexes is constructed in the same method.

(44) The determination matrix is as follows:
A=(A.sub.ij).sub.nn(19)

(45) Where A.sub.ij is an element of the determination matrix A, a scale indicating the importance degree of a factor i and a factor j relative to each other, and n is an order of the determination matrix A, namely the number of the primary evaluation indexes or the secondary evaluation indexes.

(46) Secondly, normalization is carried out, i.e., calculating a product of elements of the determination matrix at every row, then calculating the nth root of the product, and finally carrying out normalization, wherein calculation formulas are as follows:

(47) 0$\begin{array}{cc}{M}_i={.Math.}_{j=1}^n{A}_{lj}& \left(20\right)\end{array}$$\begin{array}{cc}{m}_i=\sqrt[n]{M_i}& \left(21\right)\end{array}$$\begin{array}{cc}{W}_i=\frac{m_i}{{.Math.}_{i=1}^n{m}_i}& \left(22\right)\end{array}$

(48) Where M.sub.i is a product of elements of the determination matrix A at the i.sup.th row; n is an order of the determination matrix A, namely the number of the primary evaluation indexes or the secondary evaluation indexes; A.sub.ij is an element of the determination matrix A; m.sub.i is the nth root of the product of the determination matrix A at the i.sup.th row (the nth root equals to the order of the determination matrix; n is the number of the primary evaluation indexes in the event of the calculation of the weights of the primary evaluation indexes, and n is the number of the secondary evaluation indexes in the event of the calculation of the weights of the secondary evaluation indexes); and W.sub.i is the weight of the i.sup.th evaluation index.

(49) Thirdly, a maximum eigenvalue .sub.max of the determination matrix based on the weight W.sub.i obtained in the above step, and the consistency check is carried out on the determination matrix, wherein calculation formulas are as follows:

(50) $\begin{array}{cc}R=AW& \left(23\right)\end{array}$$\begin{array}{cc}{}_{m\mathrm{ax}}={.Math.}_{i=1}^n\frac{R_i}{n{W}_i}& \left(24\right)\end{array}$$\begin{array}{cc}{C}_1=\frac{{}_{\mathrm{ma}x}-n}{n-1}& \left(25\right)\end{array}$

(51) $\begin{array}{cc}{C}_R=\frac{c_I}{R_I}& \left(26\right)\end{array}$

(52) Where R is a newly constructed matrix (n1); A is a determination matrix; W is a weight matrix (n1) of the evaluation indexes; .sub.max is the maximum eigenvalue; R.sub.i is the i.sup.th element of the newly constructed matrix R; C.sub.I is a consistency index; C.sub.R is a consistency ratio; R.sub.I is a random consistency index, and when n is 1, 2, 3, 4, 5, 6, 7, 8, and 9, R.sub.I is 0, 0, 0.58, 0.90, 1.12, 1.24, 1.32, 1.41, and 1.45, respectively.

(53) The consistency check is satisfied when the calculated consistency ratio C.sub.R is less than 0.1.

(54) Step 4. The secondary evaluation indexes of the combined heating system are normalized, and calculation formulas of membership functions of EER, clean energy utilization rate, carbon dioxide emission and operation cost are as follows:

(55) $\begin{array}{cc}{f}_1\left(\mathrm{EER}\right)=\{\begin{array}{cc}0& \mathrm{EER}<{\mathrm{EER}}_{min}\\ \frac{\mathrm{EER}-{\mathrm{EER}}_{min}}{{\mathrm{EER}}_{\mathrm{ma}x}-{\mathrm{EER}}_{min}}& {\mathrm{EER}}_{min}\mathrm{EER}<{\mathrm{EER}}_{m\mathrm{ax}}\\ 1& \mathrm{EER}{\mathrm{EER}}_{m\mathrm{ax}}\end{array}& \left(27\right)\end{array}$$\begin{array}{cc}{f}_2\left(\right)=\{\begin{array}{cc}0& <{}_{min}\\ \frac{-{}_{min}}{{}_{\mathrm{ma}x}-{}_{min}}& {}_{min}<{}_{m\mathrm{ax}}\\ 1& {}_{m\mathrm{ax}}\end{array}& \left(28\right)\end{array}$$\begin{array}{cc}{f}_3\left(T_{\mathrm{CO}2}\right)=\{\begin{array}{cc}1& T_{\mathrm{CO}2}{T}_{\mathrm{CO}{2}_{\mathrm{ma}x}}\\ \frac{T_{\mathrm{CO}{2}_{m\mathrm{ax}}}-T_{{\mathrm{CO}}_2}}{T_{\mathrm{CO}{2}_{m\mathrm{ax}}}-T_{\mathrm{CO}{2}_{min}}}& T_{\mathrm{CO}{2}_{min}}{T}_{\mathrm{CO}2}<T_{\mathrm{CO}{2}_{m\mathrm{ax}}}\\ 0& T_{\mathrm{CO}2}<T_{\mathrm{CO}{2}_{min}}\end{array}& \left(29\right)\end{array}$$\begin{array}{cc}{f}_4\left(C_r\right)=\{\begin{array}{cc}1& C_r{C}_{r_{m\mathrm{ax}}}\\ \frac{c_{r_{m\mathrm{ax}}}-c_r}{c_{r_{m\mathrm{ax}}}-c_{r_{min}}}& C_{r_{min}}{C}_r<C_{r_{\mathrm{ma}x}}\\ 0& C_r<C_{r_{min}}\end{array}& \left(30\right)\end{array}$

(56) Where EER.sub.max and EER.sub.min are EERs of the combined heating system in the optimal operation state and the worst operation state, respectively; .sub.max and .sub.min represent clean energy utilization rates of the combined heating system in the optimal operation state and the worst operation state, respectively; T.sub.CO2max and T.sub.CO2min are carbon dioxide emissions of the combined heating system in the worst operation state and the optimal operation state, respectively; and C.sub.max and C.sub.rmin are operation costs of the combined heating system in the worst operation state and the optimal operation state, respectively.

(57) When .sub.1(EER) and .sub.2() are 1, the combined heating system is in the optimal operation state, and calculated values of the EER and the clean energy utilization rate are the maximum; and when .sub.1(EER) and .sub.2() are 0, the combined heating system is in the worst operation state, and calculated values of the EER and the clean energy utilization rate are the minimum.

(58) When .sub.3(T.sub.CO2) and .sub.4(C.sub.r) are 1, the combined heating system is in the worst operation state, and calculated values of the carbon dioxide emission T.sub.CO2 and the operation cost C.sub.r are the maximum; and when .sub.3(T.sub.CO2) and .sub.4(C.sub.r) are 0, the combined heating system is in the optimal operation state, and calculated values of the carbon dioxide emission T.sub.CO2 and the operation cost C.sub.r are the minimum.

(59) Step 5. a comprehensive objective function calculation model is constructed based on the membership functions of the secondary evaluation indexes, wherein a calculation formula is as follows:
=w.sub.1(w.sub.11.sub.1(EER)+w.sub.12.sub.2())+w.sub.2.sub.3(T.sub.CO2)+W.sub.3.sub.4(C.sub.r)(31)

(60) Where w.sub.1, w.sub.2 and w.sub.3 are the weights of the primary evaluation indexes: energy conservation, environmental protection and economical efficiency, respectively; and w.sub.11 and w.sub.12 are the weights of the secondary evaluation indexes: the EER and the clean energy utilization rate, respectively.

(61) Step 6. floor radiant heating is implemented at a heating terminal form of the combined heating system, a heating parameter prediction model is constructed based on characteristics of the heating terminal form, and an actual operation water supply temperature t.sub.g, an actual operation return water temperature t.sub.h and a heating load Q at different outdoor temperatures are obtained through calculation in a constant flow control method that only the water supply temperature of a network is changed at a heat source, while the circulation flow rate of the network remains constant at a design level, wherein formulas are as follows:

(62) $\begin{array}{cc}\{\begin{array}{c}{t}_g=t_n+0.5\left(t_g^{}+t_h^{}-2t_n\right){\stackrel{\_}{Q}}^{\frac{1}{1+C_1}}+0.5\left(t_g^{}-t_h^{}\right)\stackrel{\_}{Q}\\ t_g=t_n+0.5\left(t_g^{}+t_h^{}-2t_n\right){\stackrel{\_}{Q}}^{\frac{1}{1+C_1}}+0.5\left(t_g^{}-t_h^{}\right)\stackrel{\_}{Q}\end{array}& \left(32\right)\end{array}$$\begin{array}{cc}Q=Q_{\mathrm{ratio}}{Q}^{}& \left(33\right)\end{array}$$\begin{array}{cc}{Q}_{\mathrm{ratio}}=\frac{t_n-t_a}{t_n-t_a^{}}& \left(34\right)\end{array}$

(63) The derivation process of the above formula is as follows:

(64) For the convenience of analysis and calculation, a basic formula for achieving a heating adjustment in the heating load is as follows, assuming that the heating load of the air-source heat pump and gas-fired heater combined heating system is in direct proportion to a change in an indoor and outdoor temperature difference:

(65) $\begin{array}{cc}{Q}_{\mathrm{ratio}}=\frac{t_n-t_a}{t_n-t_a^{}}=\frac{{\left(t_g+t_h-2t_n\right)}^{1+C_1}}{{\left(t_g^{}+t_h^{}-2t_n\right)}^{1+C_1}}=G_{\mathrm{ratio}}\frac{t_g-t_h}{t_g^{}-t_n^{}}& \left(35\right)\end{array}$

(66) Where t.sub.n is a calculated indoor heating temperature, C., which is set to 18 C.; t.sub.a is a calculated outdoor heating temperature, C., which can be found in design manuals (Lu Yaoqing, Practical Heating and Air Conditioning Design Manual (in Chinese), China Architecture & Building Press, 2nd Edition, p. 187), with 7 C. for Tianjin; t.sub.a is an outdoor ambient temperature, C., which derives from input outdoor meteorological data; t.sub.g is a designed water supply temperature, C.; t.sub.h is a designed return water temperature, C.; t.sub.g is an actual operation water supply temperature, C.; t.sub.h is an actual operation return water temperature, C.; c.sub.1 is a characteristic parameter of a radiant panel in floor radiant heating, which is set to 0.3; Q is a designed heating load; Q.sub.ratio is a relative heating load ratio; and G.sub.ratio is a relative flow rate ratio.

(67) With the constant flow control method, namely, a ratio of the circulation flow rate to the designed circulation flow rate is 1, G.sub.ratio=1 is substituted into the formula (35), and the actual operation water supply temperature, t.sub.g and the actual operation return water temperature, t.sub.h, at different outdoor temperatures can be obtained, as shown in the calculation formulas (32).

(68) Step 7. In an operation mode of the combined heating system, as shown in FIG. 2, it is determined whether the operation of the air-source heat pump alone can meet the requirement of the heating load at the current outdoor temperature in real time, by taking the priority of meeting the heating requirement as a principle, if not, the comprehensive objective function in a combined operation mode of the air-source heat pump and the gas-fired heater, and the comprehensive objective function in an operation mode of the gas-fired heater alone are respectively calculated, whichever the comprehensive objective function is greater; and if the air-source heat pump can meet the requirement of the heating load, the comprehensive objective function in an operation mode of the air-source heat pump alone and the comprehensive objective function in the operation mode of the gas-fired heater alone are calculated, respectively, whichever the comprehensive objective function is greater, and based on the operation mode with the greater comprehensive objective function, an operation control strategy of the air-source heat pump and gas-fired heater combined heating system is obtained.

Embodiment 1

(69) Taking a residential building in Tianjin City as an example, an operation control method of an air-source heat pump and gas-fired heater combined heating system to be employed is optimally designed, and a designed heating load is 8 kW. Taking meteorological data on a typical day as an example for calculation and analysis, wherein an air-source heat pump used for a heat source device has a rated heating capacity of 7.5 kW, and rated power of 1.85 kW; and a gas-fired heater used for the heat source device has rated heat input power of 20 kW, rated heat output power of 18 kW, and rated heat efficiency of 91.41%.

(70) (1) Mathematical model of air-source heat pump and gas-fired heater combined heating system

(71) Mathematical Model of Air-Source Heat Pump

(72) A calculation formula of the heating capacity of an air-source heat pump unit is as follows:

(73) $Q_p=\{\begin{array}{c}\begin{array}{c}{Q}_{\mathrm{rate}}\mathrm{when}\mathrm{an}\mathrm{air}\mathrm{source}\mathrm{heat}\mathrm{pump}\mathrm{is}\mathrm{under}a\\ \mathrm{non}\mathrm{frosting}\mathrm{condition}\end{array}\\ \begin{array}{c}{Q}_{\mathrm{real}}\mathrm{when}\mathrm{the}\mathrm{air}\mathrm{source}\mathrm{heat}\mathrm{pump}\mathrm{is}\mathrm{under}a\\ \mathrm{frosting}\mathrm{condition}\mathrm{and}\mathrm{the}\mathrm{defrosting}\mathrm{condition}\end{array}\end{array}$

(74) A calculation formula of the power consumption of the air-source heat pump unit is as follows:

(75) $W_p=\{\begin{array}{c}\begin{array}{c}{W}_{\mathrm{rate}}\mathrm{when}\mathrm{the}\mathrm{air}\mathrm{source}\mathrm{heat}\mathrm{pump}\mathrm{is}\mathrm{under}\mathrm{the}\\ \mathrm{non}\mathrm{frosting}\mathrm{condition}\end{array}\\ \begin{array}{c}{W}_{\mathrm{real}}\mathrm{when}\mathrm{the}\mathrm{air}\mathrm{source}\mathrm{heat}\mathrm{pump}\mathrm{is}\mathrm{under}\mathrm{the}\\ \mathrm{frosting}\mathrm{condition}\mathrm{and}\mathrm{the}\mathrm{defrosting}\mathrm{condition}\end{array}\end{array}$

(76) Data fitting is carried out according to operation performance curves provided by an air-source heat pump manufacturer, 3 water outlet temperatures can be set for the air-source heat pump, which are 45 C., 50 C. and 55 C., respectively, the operation performance curves are different at different temperatures, and specific fitting results are as follows:

(77) When the water supply temperature is 45 C., fitting formulas of the heating capacity and COP are as follows:
Q.sub.rate=0.1281t.sub.a.sup.3+4.1075t.sub.a.sup.2+160.9t.sub.a+4017.4 R.sup.2=0.9636
COP.sub.rate=0.00004t.sub.a.sup.3+0.0011t.sub.a.sup.2+0.0821t.sub.a+2.7904 R.sup.2=0.9704

(78) When the water supply temperature is 50 C., fitting formulas of the heating capacity and the COP are as follows:
Q.sub.rate=0.1182t.sub.a.sup.3+3.4609t.sub.a.sup.2+167.45t.sub.a+3992.4 R.sup.2=0.9707
COP.sub.rate=0.00003t.sub.a.sup.3+0.0007t.sub.a.sup.2+0.0811t.sub.a+2.6557 R.sup.2=0.9777

(79) When the water supply temperature is 55 C., fitting formulas of the heating capacity and the COP are as follows:

(80) $Q_{\mathrm{rate}}=-0.1738t_a^3+6.2669t_a^2+149.21t_a+3671.1R^2=0\text{.9408}$${\mathrm{COP}}_{\mathrm{rate}}=-0.00003t_a^3+0.001t_a^2+0.0722t_a+2.4966R^2=0\text{.9766}$$W_{\mathrm{rate}}=\frac{Q_{rate}}{CO{P}_{rate}}$

(81) Where Q.sub.rate is the heating capacity of the air-source heat pump unit under the non-frosting condition, W; COP.sub.rate is a coefficient of performance of the air-source heat pump unit under the non-frosting condition; W.sub.rate is power consumption of the air-source heat pump unit under the non-frosting condition, W; and t.sub.a is the outdoor ambient temperature, C.

(82) According to a correction model of the heating capacity of the air-source heat pump under the frosting condition and the defrosting condition provided by a related reference, the heating capacity, Q.sub.real, of the air-source heat pump unit under the frosting condition and the defrosting condition is obtained, wherein a calculation formula is as follows:
Q.sub.real=Q.sub.rate[0.311t.sub.a+0.043t.sub.a.sup.2+0.005t.sub.a.sup.3(0.7831.07210.sup.4t.sub.a.sup.3)RH.sup.0.8461.647]

(83) Where Q.sub.reaL is the heating capacity of the air-source heat pump unit under the frosting condition, W; Q.sub.rate is the heating capacity of the air-source heat pump unit under the non-frosting condition, W; t.sub.a is the outdoor ambient temperature, C.; and RH is the relative ambient humidity, %.

(84) Similarly, according to a correction model of the coefficient of performance of the air-source heat pump under the frosting condition provided by a related reference, coefficient of performance, COP.sub.real, of the air-source heat pump unit under the frosting condition and the defrosting condition is obtained when the outdoor temperature is higher and not higher than 7 C., wherein a calculation formula is as follows:

(85) ${\mathrm{COP}}_{real}=\{\begin{array}{c}{\mathrm{COP}}_{\mathrm{rate}}\left(1-0.1801e^{-\frac{t_a^2}{5}}\right),t_a>7\\ {\mathrm{COP}}_{\mathrm{rate}}\left[1+0.0027\left(t_a-7\right)-0.1801e^{-\frac{t_a^2}{5}}\right],t_{a}7\end{array}$

(86) Where COP.sub.real is the coefficient of performance of the air-source heat pump unit under the frosting condition; COP.sub.rate is the coefficient of performance of the air-source heat pump unit under the non-frosting condition; and t.sub.a is the outdoor ambient temperature, C.

(87) Based on this, the power consumption, W.sub.real, of the air-source heat pump unit under the frosting condition and the defrosting condition is obtained, wherein a calculation formula is as follows:

(88) 0$W_{real}=\frac{Q_{real}}{CO{P}_{real}}$

(89) Where W.sub.real is the power consumption of the air-source heat pump unit under the frosting condition, W; Q.sub.rate is the heating capacity of the air-source heat pump unit under the frosting condition, W; and COP.sub.real is the coefficient of performance of the air-source heat pump unit under the frosting condition.

(90) Mathematical Model of Gas-Fired Heater

(91) A calculation formula of a heating capacity and energy consumption model for the gas-fired heater unit, namely gas consumption per second, V.sub.g, of the gas-fired heater unit, is as follows:
V.sub.9=Q.sub.b/(H.sub.i.Math..sub.b)

(92) Where V.sub.g is the gas consumption per second of the gas-fired heater unit, m.sup.3/s; Q.sub.b is the heating capacity of the gas-fired heater unit, W; H.sub.i is the calorific value of fired gas fed into the gas-fired heater, J/m.sup.3; and .sub.b is the heat efficiency of the gas-fired heater, %.

(93) Wherein a calculation formula of Q.sub.b is as follows:

(94) $Q_b=\{\begin{array}{c}\begin{array}{c}Q-Q_p\mathrm{in}\mathrm{the}\mathrm{event}\mathrm{of}\mathrm{combined}\mathrm{operation}\mathrm{of}\mathrm{the}\\ \mathrm{gas}\mathrm{fired}\mathrm{heater}\mathrm{and}\mathrm{the}\mathrm{air}\mathrm{source}\mathrm{heat}\mathrm{pump}\end{array}\\ Q\mathrm{in}\mathrm{the}\mathrm{event}\mathrm{of}\mathrm{the}\mathrm{operation}\mathrm{of}\mathrm{the}\mathrm{gas}\mathrm{fired}\mathrm{heater}\mathrm{alone}\end{array}$

(95) Where Q is an actual heating load, kW; and Q.sub.p is the heating capacity of the air-source heat pump unit, W.

(96) Data fitting is carried out according to an operation performance curve provided by a manufacturer of the gas-fired heater to obtain the operation heat efficiency .sub.b of the gas-fired heater unit at different load rates, wherein a calculation formula is as follows:

(97) ${}_b=66.819+80.252.Math.{}_b-93.267.Math.{}_b^2+37.677.Math.{}_b^3{R}^2=0\text{.979}$${}_b=\frac{Q_b}{Q_{b0}}$

(98) Where .sub.b is the operation heat efficiency of the gas-fired heater unit, %; .sub.b is a load rate of the gas-fired heater unit; Q.sub.b is the heating capacity of the gas-fired heater unit, W; and Q.sub.b0 is the rated heating capacity of the gas-fired heater unit, W.

(99) (2) Calculation Model of Evaluation Index

(100) A comprehensive evaluation index system of the combined heating system includes primary evaluation indexes and secondary evaluation indexes, wherein 3 primary evaluation indexes are provided, including energy conservation, environmental protection and economical efficiency, and 4 secondary evaluation indexes are provided, including EER and a clean energy utilization rate () that are subordinate to the energy conservation index, a carbon dioxide T.sub.CO2 that is subordinate to the environment protection index, and operation cost C.sub.r that is subordinate to the economical efficiency index.

(101) A calculation formula of the EER of the combined heating system is as follows:

(102) $\mathrm{EER}=\frac{Q_p+Q_b}{W_P/{}_e+V_g.Math.H_i}$

(103) Where .sub.b is power generation efficiency of the gas-fired device, which is set to 40%;

(104) A calculation formula of the clean energy utilization rate () of the combined heating system is as follows:

(105) $=\frac{Q_p}{Q_p+Q_b}$

(106) Calculation formulas of the carbon dioxide emission, T.sub.CO2, of the combined heating system are as follows:
T.sub.CO2=T.sub.b+T.sub.p
T.sub.b=V.sub.9a.sub.b
T.sub.p=W.sub.pa.sub.p

(107) Where T.sub.CO2 is the carbon dioxide emission of the combined heating system, kg/s; T.sub.b is the carbon dioxide emission of the gas-fired heater unit, kg/s; T.sub.P is the carbon dioxide emission of the air-source heat pump unit, kg/s; a.sub.b is the conversion coefficient of the gas consumption to the carbon dioxide emission of the gas-fired heater, which is 1.964 kg/m.sup.3; and a.sub.p is a conversion coefficient of the power generation capacity to the carbon dioxide emission, which is 0.728 kg/kWh, so that carbon dioxide emission is obtained, which is 0.210.sup.6 kg/(W.Math.s).

(108) A calculation formula of the operation cost, C.sub.r, of the combined heating system is as follows:
C.sub.r=W.sub.pP.sub.e/(3.610.sup.6)+V.sub.gP.sub.g

(109) Where C.sub.r is the operation cost, yuan/s; p.sub.e is an electricity price, Yuan/kWh, which is set to 0.5 Yuan/kWh; and p.sub.g is a fired gas price, Yuan/m.sup.3, which is set to 2.2 Yuan/m.sup.3.

(110) (3) Determination of Weights of Evaluation Indexes

(111) Based on the evaluation indexes of the combined heating system, a determination matrix is constructed via an expert scoring method, and the weights of the evaluation indexes are calculated, wherein calculation results are as follows:

(112) TABLE-US-00001 TABLE 1 Determination matrix of primary index Energy Environmental Economical Primary index conservation protection efficiency Weight Energy conservation 1 3 5 0.6370 Environmental 1/3 1 3 0.2583 protection Economical 1/5 1/3 1 0.1047 efficiency

(113) Maximum eigenvalue .sub.max=3.0385, CI=0.0193, CR=CI/RI=0.0322<0.1, which met the consistency check.

(114) TABLE-US-00002 TABLE 2 Determination matrix of secondary index Clean energy Secondary index EER utilization rate Weight EER 1 3 0.75 Clean energy utilization rate 1/3 1 0.25

(115) Maximum eigenvalue: .sub.max=2.00, CI=0, CR=CI/RI=0<0.1, which met the consistency check.

(116) (4) Determination of Membership Function of Evaluation Index

(117) The EER in the system evaluation index is a ratio of total heat provided by the combined heating system to the gas consumption of the combined heating system, and the maximum value corresponds to the condition where the combined heating system is totally heated by the air-source heat pump that operates in the optimal operation state, as a result, COP is 5.70,and the power generation efficiency of the gas-fired boiler is 40%, namely EER=2.28; and the minimum value corresponds to the condition where the combined heating system is totally heated by the gas-fired heater for heating that operates in the worst operation state, namely EER=0.70, so that a membership function of the EER is determined, wherein a calculation formula is as follows:

(118) $f_1\left(EER\right)=\{\begin{array}{cc}0& \mathrm{EER}<0.7\\ \frac{\mathrm{EER}-0.7}{1.58}& 0.7\mathrm{EER}<2.28\\ 1& \mathrm{EER}2.28\end{array}$

(119) The clean energy utilization rate q in the system evaluation index is a ratio of total heat provided by the air-source heat pump to the total heat provided by the combined heating system, and the maximum value corresponds to the condition where the combined heating system is totally heated by the air-source heat pump, namely =1; and the minimum value corresponds to the condition where the combined heating system is totally heated by the gas-fired heater for heating, namely =0, so that a membership function of the clean energy utilization rate () is determined, wherein a calculation formula is as follows:
.sub.2()=

(120) The carbon dioxide emission T.sub.CO2 in the system evaluation index is the sum of the carbon dioxide emission of the gas-fired heater unit and the carbon dioxide emission of the air-source heat pump unit, and the maximum value corresponds to the condition where the combined heating system is totally heated by the gas-fired heater that operates in the worst operation state, namely, the heating capacity of the heater is 8 kW, the heat efficiency is 70%, the fired gas calorific value is 3.710.sup.7 J/m.sup.3, the gas consumption is obtained through calculation from the formula (12), which is 3.0910.sup.4 m.sup.3/s, and consequently, the carbon dioxide emission can be obtained through calculation based on the above conditions, which is 6.0710.sup.4 kg/s; and the minimum value corresponds to the condition where the combined heating system is totally heated by the air-source heat pump that operates in the optimal operation state, namely, the heating capacity of the heat pump is 8 kW, COP is 5.70, power consumption is obtained through calculation from the formula (3), which is 1.410.sup.3 W, and consequently, the carbon dioxide emission is obtained through calculation from the above conditions, which is 2.810.sup.4 kg/s, so that a membership function of the carbon dioxide emission T.sub.CO2 is determined, wherein a calculation formula is as follows:

(121) $f_3\left(T_{\mathrm{CO}2}\right)=\{\begin{array}{cc}1& T_{\mathrm{CO}2}6.07{10}^{-4}\\ \frac{6.07-T_{\mathrm{CO}2}{10}^4}{3.27}& 2.8{10}^{-4}{T}_{\mathrm{CO}2}<6.07{10}^{-4}\\ 0& T_{\mathrm{CO}2}<2.8{10}^{-4}\end{array}$

(122) The operation cost C.sub.r in the system evaluation index is the sum of the operation cost of the air-source heat pump and the operation cost of the gas-fired heater, and the maximum value corresponds to the condition where the combined heating system is totally heated by the gas-fired heater that operates in the worst operation state, namely, the heating capacity of the heater is 8 kW, the heat efficiency is 83%, the fired gas calorific value is 3.710.sup.7 J/m.sup.3, the gas consumption is obtained through calculation from the formula (12), which is 2.6110.sup.4 m.sup.3/s, the fired gas price is 2.2 Yuan/m.sup.3, and consequently, the operation cost is obtained through calculation from the above conditions, which is 5.2310.sup.4 Yuan/s; and the minimum value corresponds to the condition where the combined heating system is totally heated by the air-source heat pump that operates in the optimal operation state, namely, the heating capacity of the heat pump is 8 kW, COP is 5.70, power consumption is obtained through calculation from the formula (3), which is 1.410.sup.3 W, the electricity price is 0.5 Yuan/kWh, and consequently, the operation cost is obtained through calculation from the above conditions, which is 1.9510.sup.4 Yuan/s, so that a membership function of the operation cost C.sub.r is determined, wherein a calculation formula is as follows:

(123) $f_4\left(C_r\right)=\{\begin{array}{cc}1& C_{r}5.23{10}^{-4}\\ \frac{5.23-C_r{10}^4}{1.95}& 3.28{10}^{-4}{C}_r<\text{.23}{10}^{-4}\\ 0& C_r<3.28{10}^{-4}\end{array}$
(5) Determination of Comprehensive Objective Function Calculation Model

(124) A Comprehensive Objective Function Calculation Model is Constructed Based on the Membership Function and Weight of the Evaluation Index at Each Level, Wherein a Calculation Formula is as Follows:
=0.0637(0.75.sub.1(EER)+0.25.sub.2())+0.2583.sub.3(T.sub.CO2)+0.1047.sub.4(C.sub.r)
(6) Determination of Heating Load

(125) An indoor heating calculating temperature t.sub.n is 18 C., and an outdoor heating calculating temperature t.sub.a in Tianjin City is 7 C., designed water supply and return water temperatures, t.sub.g and t.sub.h, are 50 C. and 40 C., respectively, a designed heating load Q is 8 kW, and consequently, a relative heating load Q.sub.ratio and an actual heating load Q are determined, wherein calculation formulas are as follows:

(126) $Q_{\mathrm{ratio}}=\frac{18-t_a}{18-\left(-7\right)}=0.72-0.04t_a$$Q=8\frac{18-t_a}{18-\left(-7\right)}=5.76-0.32t_a$

(127) With the constant flow control method, namely, a ratio of the circulation flow rate to the designed circulation flow rate of 1, a characteristic parameter c.sub.1 of a radiant panel is 0.3, so that a heating and water supply temperature, a heating and water return temperature and a heating load at different outdoor temperatures are obtained, wherein calculation formulas are as follows:

(128) $\{\begin{array}{c}{t}_g=18+27{\left(0.72-0.04t_a\right)}^{0.77}+0.5\left(0.72-0.04t_a\right)\\ t_h=18+27{\left(0.72-0.04t_a\right)}^{0.77}-0.5\left(0.72-0.04t_a\right)\end{array}$
(7) Determination of Operation Control Strategy of Combined Heating System

(129) It is determined whether the operation of the air-source heat pump alone can meet the requirement of the heating load at the current outdoor temperature in real time by taking the priority of meeting the heating requirement as a principle, if not, a comprehensive objective function in a combined operation mode of the air-source heat pump and the gas-fired heating heater for heating and a comprehensive objective function in an operation mode of the gas-fired heater for heating alone are respectively calculated, whichever the comprehensive objective function is greater; and if the air-source heat pump can meet the requirement of the heating load, a comprehensive objective function in an operation mode of the air-source heat pump alone and the comprehensive objective function in the operation mode of the gas-fired heater for heating alone are calculated, respectively, whichever the comprehensive objective function is greater, and based on the operation mode with the greater comprehensive objective function, an operation control strategy of the air-source heat pump and gas-fired heater combined heating system is obtained.

(130) Taking the meteorological data of Tianjin City on a typical day during the winter season as an example, the operation control strategy for the air-source heat pump and gas-fired heater combined heating system is optimally calculated, and results are shown in the figure. Hourly meteorological data, including outdoor temperature, solar irradiance, and calculated actual water supply temperature are shown in FIG. 3, an hourly output of a heat-source device is shown in FIG. 4, and hourly evaluation indexes and comprehensive objective function values of the combined heating system are shown in FIG. 5.

(131) As shown in FIG. 3, the hourly water supply temperatures of the air-source heat pump and gas-fired heater combined heating system correspond to set water supply temperatures for the air-source heat pump and the gas-fired heater for heating at different moments; as shown in FIG. 4, during the time period from 6:00 to 7:00 a.m., the outdoor temperature was relatively lower, the air-source heat pump was in the poor operation state, and therefore, the gas-fired heater operated exclusively during this time period; during the time period from 13:00 to 17:00 p.m., the outdoor temperature was relatively higher, the air-source heat pump was in the suboptimal operation state and could meet the requirements of heating load, and therefore, the air-source heat pump operated alone during this time period; and at other moments, the air-source heat pump operated together with the gas-fired heate. As shown in FIG. 5, the system evaluation index, including the EER, the clean energy utilization rate (), and the comprehensive objective function (), had a change trend similar to that of the outdoor temperature, i.e., initially increasing with the rise in the outdoor temperature and then decreasing as the outdoor temperature dropped. Conversely, the change trends for the carbon dioxide emissions (T.sub.CO2) and the operation cost (C.sub.r) were opposite to the change trend of the outdoor temperature, i.e., decreasing with the rise in the outdoor temperature and increasing as the outdoor temperature dropped.

(132) It should be noted that although the preferred embodiments of the present disclosure are described above with reference to the drawings, the present disclosure is not merely limited to the above specific implementations, and the above implementations are only schematic instead of restrictive. Those skilled in the art may also make many forms without departing from the purpose of the present disclosure and the scope protected by the claims under the inspiration of the present disclosure, and these forms all belong to the scope of protection of the present disclosure.