Modeling Method of Combined Heat and Power Optimal Dispatching Model

Abstract

A CHP optimal dispatching model is a mixed integer programming model and is used for a district heating system (DHS) comprising a heat source, a heating network and a heat load, and the heating network comprises a heat transmission network and a heat distribution network. A plurality of heating areas is divided, and one day is divided into a plurality of time periods; the heat transmission loss of the heat distribution network is omitted, and a heat transmission network model taking transmission time delay of the heating network into consideration is established according to the heat transmission network; a terminal heat consumer model capable of reflecting indoor temperature is established; and a combined optimal dispatching model comprising conventional generators, wind power units, CHP units, electric boilers and heat storage tanks is established.

Claims

1. A method for dispatching an optimal combined heat and power (CHP) module, characterized in that the optimal CHP module is used for a district heating system (DHS) comprising a heat source, a heating network and a heat load, the heating network comprises a heat transmission network and a heat distribution network, and the method comprises the following steps: step 1: dividing a plurality of heating areas based on a distance from the heat load to the heat source, and dividing one day into a plurality of time periods; step 2: omitting the heat transmission loss of the heat distribution network, and establishing a heat transmission network module taking transmission time delay of the heating network into consideration based on the heat transmission network; step 3: establishing a terminal heat consumer module that being capable of reflecting indoor temperature according to the heat load; and step 4: establishing the optimal CHP module comprising conventional generators, wind power units, CHP units, electric boilers and heat storage tanks based on the heat source.

2. The method according to claim 1, characterized in that establishing the heat transmission network in the step 2 comprises: calculating the heat loss of each heating area within each time period and establishing a confluence node module of the heat transmission network, and describing a steady state heat transmission of the whole heat transmission network.

3. The method according to claim 2, characterized in that establishing the heat transmission network module in the step 2 comprises: calculating unit pipeline length according to a velocity of a water flow and a dispatching time interval:
L=v.Math.t wherein L is the unit pipeline length, v is the velocity of the water flow, and t is the dispatching time interval; (1) calculation of the thermal resistance of soil and the temperature loss of hot water in a pipeline with unit length according to a basic theory of heat transfer as shown as $\begin{array}{cc}{R}_e=\frac{1}{2.Math.{}_e}.Math.\ln \left(\frac{2.Math.\left(h+{}_e/a_r\right)}{d_{\mathrm{ex}}}+\sqrt{{\left(\frac{2.Math.\left(h+{}_e/a_r\right)}{d_{\mathrm{ex}}}\right)}^2-1}\right)& t\end{array}$ $\begin{array}{cc}.Math..Math.T_t^{S,k}=\frac{\left(1+\right).Math.L}{c_{\mathrm{water}}.Math.{\mathrm{MS}}_t^k}.Math.\frac{\left(T_t^{S,k}-T_t^{\mathrm{surf}}\right).Math.2.Math.{}_b}{\ln \left(d_{\mathrm{ex}}/d_{i.Math..Math.n}\right)+2.Math.{}_b.Math.R_e}& t,kK\end{array}$ $\begin{array}{cc}.Math..Math.T_t^{R,k}=\frac{\left(1+\right).Math.L}{c_{\mathrm{water}}.Math.{\mathrm{MR}}_t^k}.Math.\frac{\left(T_t^{R,k}-T_t^{\mathrm{surf}}\right).Math.2.Math.{}_b}{\ln \left(d_{\mathrm{ex}}/d_{i.Math..Math.n}\right)+2.Math.{}_b.Math.R_e}& t,kK\end{array}$ wherein R.sub.e is the thermal resistance of the soil, T.sub.t.sup.S,k represents for the temperature loss of a water supply network in the k.sup.th heating area within the t.sup.th time period, T.sub.t.sup.R,k represents for the temperature loss of a water return network in the k.sup.th heating area within the t.sup.th time period, h represents for the distance between the center of a heating pipeline and a soil surface, .sub.e represents for the soil thermal conductivity coefficient, .sub.r represents for the heat transmission coefficient of the soil surface, d.sub.in and d.sub.ex respectively represent for the internal and external diameters of the pipeline, represents for an additional factor of heat loss caused by appendixs, .sub.b represents for the thermal conductivity of an insulation material on the surface of the pipeline, c.sub.water represents for the specific heat capacity of water, T.sub.t.sup.surf represents for the temperature of the soil surface at period t, MS.sub.t.sup.k and MR.sub.t.sup.k respectively represent for mass flow rates of hot water in a water supply pipeline and a water return pipeline in the k.sup.th heating area at period t, T.sub.t.sup.S,k and T.sub.t.sup.R,k respectively represent for water temperatures of the water supply pipeline and a water return pipe network in an area K at period t, represents for a set of dispatching periods, and K represents for a set of heating areas of heat loads; (2) continuity constraints of mass flow rates of nodes considering time delay as shown as $\begin{array}{cc}\{\begin{array}{c}{\mathrm{MS}}_t^k=m_{t+1}^k+{\mathrm{MS}}_{t+1}^{k+1}\\ {\mathrm{MS}}_t^{k.\mathrm{end}}=m_t^{k.\mathrm{end}}\end{array}& k,k+1K,\end{array}.Math.t$ $\begin{array}{cc}\{\begin{array}{c}{\mathrm{MR}}_t^k=m_t^k+{\mathrm{MR}}_{t+1}^{k+1}\\ {\mathrm{MR}}_t^{k.\mathrm{end}}=m_t^{k.\mathrm{end}}\end{array}& k,k+1K,\end{array}.Math.t$ wherein k.end represents for the last heating area, and m.sub.t.sup.k represents for the mass flow rate of hot water flowing through the k.sup.th heat-exchange station; (3) water temperature at a confluence node as shown as
MR.sub.t+1.sup.k=T.sub.t+1.sup.R,k=m.sub.t.sup.k.sub.t+1.sup.R,k+MR.sub.t.sup.k+1(T.sub.t.sup.R,k+1T.sub.t.sup.R,k+1) kK,t wherein .sub.t+1.sup.S,k and .sub.t+1.sup.R,k respectively represent for the water supply temperature and the water return temperature of the k.sup.th heat-exchange station at period t; (4) temperature loss of hot water in pipelines in the adjacent heating areas: $\{\begin{array}{cc}\begin{array}{c}{T}_{t+1}^{S,1}=T_t^S-.Math..Math.T_t^{S,1}\\ T_{t+1}^{S,k+1}=T_t^{S,k}-.Math..Math.T_t^{S,k}\\ T_{t+1}^R=T_t^{R,k}-.Math..Math.T_t^{R,k}\\ T_t^{S,k}={}_t^{S,k}\end{array}& kK,t\end{array}$ (5) pipeline flow rate limitation as shown as
0m.sub.t.sup.km.sup.k kK,t
0MS.sub.t.sup.kMS.sup.k kK,t
0MR.sub.t.sup.kMR.sup.k kK,t wherein m.sup.k represents for the maximum mass flow rate flowing through the k.sup.th heat-exchange station, and MS.sup.k and MR.sup.k respectively represent for the maximum mass flow rates of the k.sup.th water supply pipeline and the k.sup.th water return pipeline; (6) pressure loss as shown as $\begin{array}{cc}{\mathrm{PS}}_t^k-{\mathrm{PS}}_t^{k+1}=\frac{.Math.f_k.Math.{\left({\mathrm{MS}}_t^k\right)}^2}{{\left(d_{i.Math..Math.n}^k\right)}^5.Math.{\left({}_t^k\right)}^2}& kK,t\end{array}$ $\begin{array}{cc}{\mathrm{PR}}_t^{k+1}-{\mathrm{PS}}_t^k=\frac{.Math.f_k.Math.{\left({\mathrm{MR}}_t^k\right)}^2}{{\left(d_{i.Math..Math.n}^k\right)}^5.Math.{\left({}_t^k\right)}^2}& kK,t\end{array}$ PS.sub.t.sup.k and PR.sub.t.sup.k respectively represent for pressures of the water supply pipeline and the water return pipeline in the k.sup.th area at period t, represents for a pressure coefficient of the pipeline, f.sub.k represents for a friction coefficient in the k.sup.th heating area, d.sub.in.sup.k represents for the internal diameter of the pipeline in the k.sup.th heating area, and .sub.t.sup.k represents for the hot water density in the k.sup.th heating area at period t; and (7) heat-exchange power of the heat-exchange station as shown as
c.sub.water.Math.m.sub.t.sup.k.Math.(.sub.t.sup.S,k.sub.t.sup.R,k)=H.sub.t,k t,kK
.sup.R,kt.sup.R,k.sup.R,k t,kK wherein H.sub.t,k is the heat-exchange power of the heat-exchange station, .sup.R,k is a lower limit of water return temperature of the k.sup.th heat-exchange station, and .sup.R,k is an upper limit of water return temperature of the k.sup.th heat-exchange station.

4. The method according to claim 1, characterized in that establishing the terminal heat consumer module in the step 3 comprises: omitting the heat transmission loss of the heat distribution network, regarding heat absorbed from the heat transmission network by the heat-exchange station as heat transmitted to a consumer side, and calculating heat loss of a building envelope, cold air invasion heat loss and cold air infiltration heat loss of a building in the heating area based on a heat engineering to obtain the average indoor temperature in one heating area so as to balance the thermal comfort degree of a consumer.

5. The method according to claim 4, characterized in that the average indoor temperature in the heating area is: (1) a real heat transfer process is very complex and comprises convection, conduction and radiation, and the heat loss of the building envelope of the building is solved by only adopting a heat loss calculation formula under a steady state environment herein in order to simplify the heat transfer process:
Q.sub.t,k.sup.BE=A.sub.kF.sub.e,k(.sub.t,k.sup.in.sub.t,k.sup.out)(1+x.sub.k.sup.g)(1+x.sub.k.sup.ch) t,kK
.sub.min.sup.in.sub.t,k.sup.in.sub.max.sup.in
.sub.ch.sub.t+1,k.sup.in.sub.t,k+1.sup.in.sub.ch wherein Q.sub.t,k.sup.BB is the heat loss of the building envelope of the building, x.sub.k.sup.g and x.sub.k.sup.ch are respectively a correction coefficient of building height and a correction coefficient of building orientation in the k.sup.th heat load area, A.sub.k represents for the area of the building envelope of the building in the k.sup.th heat load area, F.sub.e,k represents for the heat transmission coefficient of the building envelope of the building in the k.sup.th heat load area, .sub.t,k.sup.in represents for the average indoor temperature in the k.sup.th heat load area at period t, .sub.t,k.sup.out represents for the outdoor temperature in the k.sup.th heat load area at period t, .sub.min.sup.in and .sub.max.sup.in respectively represent for minimum average indoor temperature and maximum average indoor temperature, .sub.ch represents for the maximum change rate of the average indoor temperature within the unit dispatching time, represents for the set of the dispatching periods, and K represents for a set of the heating areas of the heat loads; (2) cold air infiltrates into a room through gaps such as doors and windows under the action of an indoor and outdoor pressure difference caused by wind power and heat pressure and escapes after being heated, and the heat consumed by heating the part of cold air from the outdoor temperature to the indoor temperature becomes cold air infiltration heat loss:
Q.sub.t,k.sup.Inf=0.278.Math.Num.Math.V.sub.kc.sub.p.sub.t.sup.out(.sub.t,k.sup.in.sub.t,k.sup.out) t,kK wherein Q.sub.t,k.sup.Inf is the cold air infiltration heat loss, 0.278 represents for a unit conversion factor, Num represents for the air change rate per hour, V.sub.k represents for the total indoor volume in the k.sup.th heat load area, c.sub.p represents for constant-pressure specific heat capacity of cold air, and .sub.t.sup.out represents for outdoor air density at period t; (3) the cold air invades into the room through opening external doors under the actions of air pressure and heat pressure in winter, and the heat consumed by heating the part of cold air to the indoor temperature is called as cold air invasion heat loss:
Q.sub.t,k.sup.Ven=0.278.Math.Ven.Math.c.sub.p.sub.t.sup.out(.sub.t,k.sup.in.sub.t,k.sup.out) t,kK wherein Q.sub.t,k.sup.Ven is the cold air invasion heat loss, and Ven represents for cold air invasion volume per hour; (4) in order to simplify calculation, an assumption is made as follows: the air temperature in the building is consistent and can be described by .sub.t,k.sup.in; the thermal mass temperature distribution in the building is uniform, which means that the velocity of the surface heat diffusion process on the thermal mass is far larger than that of convection heat exchange; all heat obtained in the building within the unit time can be uniformly expressed by H.sub.t,k; based on the assumption, the average indoor temperature is:
Q.sub.t,k.sup.BE+Q.sub.t,k.sup.Inf+Q.sub.t,k.sup.Ven+C.sub.MM.sub.k(.sub.t+1,k.sup.in.sub.t,k.sup.in)/t=H.sub.t,k wherein H.sub.t is the average indoor temperature, C.sub.M represents for the specific heat capacity of the thermal mass, M.sub.k represents for the mass of the thermal mass in the k.sup.th heat load area, and the product of C.sub.M and M.sub.k can be obtained by an engineering experiment.

6. The method according to claim 1, characterized in that establishing the optimal CHP module in the step 4 comprises: (1) CHP units as shown as $\begin{array}{cc}0h_{i,t}{\stackrel{\_}{h}}_i& i{I}^{\mathrm{CHP}},t\end{array}$ $\{\begin{array}{cc}\begin{array}{c}{p}_{i,t}{\stackrel{\_}{p}}_i-C_i^1.Math.h_{i,t}\\ p_{i,t}\max .Math.\left\{C_i^{\max }.Math.h_{i,t}+C_i,{\underset{\_}{p}}_i-C_i^2.Math.h_{i,t}\right\}\end{array}& i{I}^{\mathrm{CHP}},t\end{array}$ p.sub.i,t is generated output of the CHP units, p.sub.i is an upper limit of the generated output of the CHP units, p.sub.i is a lower limit of the generated output of the CHP units, C.sub.i.sup.max, C.sub.i.sup.1 and C.sub.i.sup.2 are heat and power coupling property factors of the CHP units, h.sub.i,t represents for heating power of the CHP units, h.sub.i represents for the maximum heating power, and I.sup.CHP represents for a set of serial numbers of CHP units; (2) electric boiler module as shown as
0BH.sub.u,tBH.sub.u uU,t
BH.sub.u,t=.Math.BP.sub.u,t uU,t BH.sub.u,t represents for the heating power of the electric boilers u at period t, BP.sub.u,t represents for electrical power consumed by the electric boilers, represents for the thermal efficiency of the electric boilers, BH.sub.u represents for the maximum heating power of the electric boilers, U represents for a set of serial numbers of the electric boilers, and represents for a set of dispatching periods; (3) establishment of heat storage tank model according to same total heat charge/discharge of heat storage tank within one day as shown as $\begin{array}{cc}{S}_{w,t+1}=S_{w,t}+{\mathrm{CS}}_{w,t}-{\mathrm{DCS}}_{w,t}& wW,t\\ f_{w,t}\left\{0,1\right\}& wW,t\\ 0S_{w,t}{\stackrel{\_}{S}}_w& wW,t\\ 0{\mathrm{CS}}_{w,t}{f}_{w,t}.Math.{\stackrel{\_}{\mathrm{CS}}}_w& wW,t\\ 0{\mathrm{DCS}}_{w,t}\left(1-f_{w,t}\right).Math.{\stackrel{\_}{\mathrm{DCS}}}_w& wW,t\end{array}$ $\underset{wW}{.Math.}.Math.\underset{t}{.Math.}.Math.{\mathrm{CS}}_{w,t}=\underset{wW}{.Math.}.Math.\underset{t}{.Math.}.Math.{\mathrm{DCS}}_{w,t}$ S.sub.w,t represents for the heat reserve capacity of the heat storage tanks w at period t, CS.sub.w,t represents for the heat charging power of the heat storage tanks w at period t, DCS.sub.w,t represents for the heat discharging power of the heat storage tanks w at period t, f.sub.w,t represents for the state of the heat storage tanks, 0 represents for heat discharging of the heat storage tanks, 1 represents for heat storage of the heat storage tanks, S.sub.w represents for the maximum heat reserve capacity of the heat storage tanks w, CS.sub.w represents for the maximum heat storage power of the heat storage tanks, DCS.sub.w represents for the maximum heat discharging power of the heat storage tanks w, and W represents for a set of serial numbers of heat storage tanks; (4) electrical power balance as shown as $\underset{i{I}^{\mathrm{CG}}}{.Math.}.Math.p_{i,t}+\underset{i{I}^{\mathrm{CHP}}}{.Math.}.Math.p_{i,t}+\underset{i{I}^{\mathrm{WP}}}{.Math.}.Math.p_{i,t}=\underset{uU}{.Math.}.Math.{\mathrm{BP}}_{u,t}+\underset{n{I}^{\mathrm{Load}}}{.Math.}.Math.{\mathrm{LD}}_{n,t}$ $t$ I.sup.CG, I.sup.WP and I.sup.Load respectively represent for a serial number of the conventional generators, a serial number of wind farms and a serial number of electrical load nodes, and LD.sub.n,t represents for load power connecting to a bus n at period t; (5) heat power constraints as shown as $\underset{uU}{.Math.}.Math.{\mathrm{BH}}_{u,t}+\underset{i{I}^{\mathrm{CHP}}}{.Math.}.Math.h_{i,t}\underset{kK}{.Math.}.Math.Q_{t+k,k}^{\mathrm{Total}}+\underset{wW}{.Math.}.Math.f_{w,t}.Math.{\mathrm{CS}}_{w,t}+\left(1-f_{w,t}\right).Math.{\mathrm{DCS}}_{w,t}$ $.Math.\begin{array}{cc}{Q}_{t+k,k}^{\mathrm{Total}}=Q_{t+k,k}^{\mathrm{BE}}+Q_{t+k,k}^{\mathrm{Inf}}+Q_{t+k,k}^{\mathrm{Ven}}& t,kK\end{array}$ Q.sub.t+k,k.sup.Total is the total heat consumption power of the building in the k heating areas at period t+k; (6) ramp rate constraints of unit as shown as
rdown.sub.i.Math.tp.sub.i,1p.sub.i,t1rup.sub.i.Math.t
iI.sup.CHPUI.sup.CG,t rup.sub.i and rdown.sub.i respectively represent for the maximum upward ramp rate and downward ramp rate of the unit, and t is the dispatching time interval; (7) spinning reserve constraints of system as shown as $\begin{array}{cc}\underset{i{I}^{\mathrm{CG}}}{.Math.}.Math.\min .Math.\left\{{\stackrel{\_}{p}}_i-p_{i,t},.Math..Math.t.Math.{\mathrm{rup}}_i\right\}{\mathrm{SR}}^{\mathrm{up}}& t\end{array}$ $\begin{array}{cc}\underset{i{I}^{\mathrm{CG}}}{.Math.}.Math.\min .Math.\left\{p_{i,t}-{\underset{\_}{p}}_i,.Math..Math.t.Math.{\mathrm{rdown}}_i\right\}{\mathrm{SR}}^{\mathrm{down}}& t\end{array}$ SR.sup.up and SR.sup.down respectively represent for the maximum upward spinning reserve constraint and downward spinning reserve constraint required by a system; (8) system load flow constraints as shown as $.Math.\underset{n{I}^{\mathrm{bus}}}{.Math.}.Math.{\mathrm{SF}}_{l,n}.Math.\left(\underset{i{S}_n^{\mathrm{CG}}.Math.{\mathrm{US}}_n^{\mathrm{CHP}}}{.Math.}.Math.p_{i,t}+\underset{i{S}_n^{\mathrm{WP}}}{.Math.}.Math.p_{i,t}-{\mathrm{LD}}_{n,t}\right).Math.F_l$ $l{I}^{\mathrm{Line}},t$ wherein SF.sub.l,n represents for a transfer factor from the bus n to a line l, F.sub.l represents for the maximum transmission capacity of the line l, I.sup.Line represents for a set of serial numbers of lines, S.sub.n.sup.CG represents for a set of indices of conventional generators, S.sub.n.sup.CHP represents for a set of indices of CHP units, and S.sub.n.sup.WP represents for a set of indices of wind generation units.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0008] FIG. 1 is a structure diagram of the heat transmission network;

[0009] FIG. 2 is a detailed hardware implementation diagram of the present invention.

[0010] FIG. 3a is an average indoor temperature variation curve obtained by using conventional modeling method;

[0011] FIG. 3b is an average indoor temperature variation curve obtained by using the proposed method;

[0012] FIG. 4a is a comparison diagram of total heat output of heat sources;

[0013] FIG. 4b is a comparison diagram of wind power accommodation;

[0014] FIG. 4c is a comparison diagram of output of the CHP unit 2;

[0015] FIG. 4d is a comparison diagram of superposed generated output of the conventional generator 1 and the conventional generator 2;

[0016] FIG. 5a shows three intraday environment temperature curves;

[0017] FIG. 5b is a wind power prediction curve when the temperature is respectively shown as FIG. 5 (a);

[0018] FIG. 6a shows influences of initial heat storage capacity of the heat storage tanks to wind power absorption; and

[0019] FIG. 6b shows influences of heat charging/discharging rate to wind power absorption.

DETAILED DESCRIPTION OF THE INVENTION

[0020] The present invention is further described below in combination with accompanying drawings and specific embodiments.

[0021] A modeling method of a CHP optimal dispatching model is characterized in that the CHP optimal dispatching model is a mixed integer programming model and is used for the district heating system (DHS) comprising heat sources, heating networks and heat loads, the heating network comprises heat transmission network (primary heating network) and heat distribution network (secondary heating network), the primary heating network is shown as FIG. 1, the circle represents for a heat-exchange station, the rectangle represents for the secondary heating network. The modeling method comprises the following steps:

[0022] Step 1: a plurality of heating areas is divided according to the distance between the heat load and the heat source, and the day is divided into a plurality of time periods.

[0023] Step 2: the heat transmission loss of the heat distribution network is omitted, and a heat transmission network model taking transmission time delay of the heating network into consideration is established according to the heat transmission network.

[0024] The transmission speed of hot water in the heat transmission network is taken into consideration, the heat loss in each heating area within each time period is calculated respectively, and a confluence node model of the heat transmission network is established, so that a steady state heat transmission model of the whole heat transmission network is described.

[0025] Firstly, unit pipeline length is calculated according to the velocity of water flow and dispatching time interval, see a formula (1):

L=v.Math.t(1)

wherein L is the unit pipeline length, v is the velocity of the water flow, and t is the dispatching time interval;
(1) the thermal resistance R.sub.e of soil can be calculated by using a formula (2) according to a basic theory of heat transfer. After the thermal resistance value of the soil is obtained, the temperature loss of hot water in a pipeline within a certain length can be calculated according to formulae (3) and (4):

[00001]$\begin{array}{cc}{R}_e=\frac{1}{2.Math.{}_e}.Math.\ln \left(\frac{2.Math.\left(h+{}_e/a_r\right)}{d_{\mathrm{ex}}}+\sqrt{{\left(\frac{2.Math.\left(h+{}_e/a_r\right)}{d_{\mathrm{ex}}}\right)}^2-1}\right).Math..Math.t& \left(2\right)\\ .Math..Math.T_t^{S,k}=\frac{\left(1+\right).Math.L}{c_{\mathrm{water}}.Math.{\mathrm{MS}}_t^k}.Math.\frac{\left(T_t^{S,k}-T_t^{\mathrm{surf}}\right).Math.2.Math.{}_b}{\ln \left(d_{\mathrm{ex}}/d_{\mathrm{in}}\right)+2.Math.{}_b.Math.R_e}.Math..Math..Math.t,kK& \left(3\right)\\ .Math..Math.T_t^{R,k}=\frac{\left(1+\right).Math.L}{c_{\mathrm{water}}.Math.{\mathrm{MR}}_t^k}.Math.\frac{\left(T_t^{R,k}-T_t^{\mathrm{surf}}\right).Math.2.Math.{}_b}{\ln \left(d_{\mathrm{ex}}/d_{\mathrm{in}}\right)+2.Math.{}_b.Math.R_e}.Math..Math..Math.t,kK& \left(4\right)\end{array}$

wherein R.sub.e is the thermal resistance of the soil, T.sub.t.sup.S,k represents for the temperature loss of a water supply network in the k.sup.th heating area within the t.sup.th time period, T.sub.t.sup.R,k represents for the temperature loss of a water return network in the k.sup.th heating area within the t.sup.th time period, h represents for the distance between the center of a heating pipeline and a soil surface, .sub.e represents for the soil thermal conductivity coefficient, .sub.r represents for the heat transmission coefficient of the soil surface, d.sub.in and d.sub.ex respectively represent for the internal and external diameters of the pipeline, represents for an additional factor of heat loss caused by appendixs, .sub.b represents for the thermal conductivity of an insulation material on the surface of the pipeline, c.sub.water represents for the specific heat capacity of water, T.sub.t.sup.surf represents for the temperature of the soil surface at period t, MS.sub.t.sup.k and MR.sub.t.sup.k respectively represent for mass flow rates of hot water in a water supply pipeline and a water return pipeline in the k.sup.th heating area at period t, T.sub.t.sup.S,k and T.sub.t.sup.R,k respectively represent for water temperatures of the water supply pipeline and a water return pipe network in an area K at period t, represents for a set of dispatching periods, and K represents for a set of heating areas of heat loads;
(2) continuity constraints of mass flow rates of nodes considering time delay:

[00002]$\begin{array}{cc}\{\begin{array}{c}{\mathrm{MS}}_t^k=m_{t+1}^k+{\mathrm{MS}}_{t+1}^{k+1}\\ {\mathrm{MS}}_t^{k.Math.\mathrm{end}}=m_t^{k.Math.\mathrm{end}}\end{array}.Math..Math..Math.k,k+1K,t& \left(5\right)\\ \{\begin{array}{c}{\mathrm{MR}}_t^k=m_t^k+{\mathrm{MR}}_{t-1}^{k+1}\\ {\mathrm{MR}}_t^{k.Math.\mathrm{end}}=m_t^{k.Math.\mathrm{end}}\end{array}.Math..Math..Math.k,k+1K,t& \left(6\right)\end{array}$

wherein k.end represents for the last heating area, and m.sub.t.sup.k represents for the mass flow rate of hot water flowing through the k.sup.th heat-exchange station;
(3) water temperature at a confluence node:

MR.sub.t+1.sup.k=T.sub.t+1.sup.R,k=m.sub.t.sup.k.sub.t+1.sup.R,k+MR.sub.t.sup.k+1(T.sub.t.sup.R,k+1T.sub.t.sup.R,k+1) kK,t(7)

wherein .sub.t+1.sup.S,k and .sub.t+1.sup.R,k respectively represent for the water supply temperature and the water return temperature of the k.sup.th heat-exchange station at period t;
(4) temperature loss of hot water in pipelines in the adjacent heating areas:

[00003]$\begin{array}{cc}\{\begin{array}{c}{T}_{t+1}^{S,1}=T_t^S-.Math..Math.T_t^{S,1}\\ T_{t+1}^{S,k+1}=T_t^{S,k}-.Math..Math.T_t^{S,k}\\ T_{t+1}^R=T_t^{R,k}-.Math..Math.T_t^{R,k}\\ T_t^{S,k}={}_t^{S,k}\end{array}.Math..Math..Math.kK,t& \left(8\right)\end{array}$

(5) pipeline flow rate limitation:

0m.sub.t.sup.km.sup.k kK,t(9)

0MS.sub.t.sup.kMS.sup.k kK,t(10)

0MR.sub.t.sup.kMR.sup.k kK,t(11)

wherein m.sup.k represents for the maximum mass flow rate flowing through the k.sup.th heat-exchange station, and MS.sup.k and MR.sup.k respectively represent for the maximum mass flow rates of the k.sup.th water supply pipeline and the k.sup.th water return pipeline;
(6) pressure loss:

[00004]$\begin{array}{cc}{\mathrm{PS}}_t^k-{\mathrm{PS}}_t^{k+1}=\frac{.Math.f_k.Math.{\left({\mathrm{MS}}_t^k\right)}^2}{{\left(d_{\mathrm{in}}^k\right)}^5.Math.{\left({}_t^k\right)}^2}.Math..Math..Math.kK,t& \left(12\right)\\ {\mathrm{PR}}_t^{k+1}-{\mathrm{PR}}_t^k=\frac{.Math.f_k.Math.{\left({\mathrm{MR}}_t^k\right)}^2}{{\left(d_{\mathrm{in}}^k\right)}^5.Math.{\left({}_t^k\right)}^2}.Math..Math..Math.kK,t& \left(13\right)\end{array}$

PS.sub.t.sup.k and PR.sub.t.sup.k respectively represent for pressures of the water supply pipeline and the water return pipeline in the k.sup.th area at period t, represents for a pressure coefficient of the pipeline, f.sub.k represents for a friction coefficient in the k.sup.th heating area, d.sub.in.sup.k represents for the internal diameter of the pipeline in the k.sup.th heating area, and .sub.t.sup.k represents for the hot water density in the k.sup.th heating area at period t; and
(7) heat-exchange power of the heat-exchange station as shown as

c.sub.water.Math.m.sub.t.sup.k.Math.(.sub.t.sup.S,k.sub.t.sup.R,k)=H.sub.t,k t,kK(14)

.sup.R,kt.sup.R,k.sup.R,k t,kK(15)

wherein H.sub.t,k is the heat-exchange power of the heat-exchange station, .sup.R,k is a lower limit of water return temperature of the k.sup.th heat-exchange station, and .sup.R,k is an upper limit of water return temperature of the k.sup.th heat-exchange station.

[0026] Step 3: a terminal heat consumer model capable of reflecting indoor temperature is established according to the heat load.

[0027] The heat transmission loss of the heat distribution network is omitted, and due to the neglection of the heat transmission loss of the heat distribution network, heat absorbed from the heat transmission network by the heat-exchange station is regarded as heat transmitted to the consumer side. Then, the heat loss of building envelopes, the cold air invasion heat loss and the cold air infiltration heat loss of buildings in the heating area can be calculated respectively. According to the related design formula of heat engineering, the average indoor temperature in one heating area can be obtained. It is used to measure the thermal comfort of residents. In this model, the average indoor temperature is a unique index for measuring the thermal comfort degree.

[0028] A residential building can also be regarded as a heat storage device with huge thermal inertia, heat can be stored in indoor air, doors and windows as well as furniture. The thermal inertia of the building can reduce peak value of heat demand and the variation rate of indoor temperature, and the peak staggered regulation of heat supply and power supply can be realized by utilizing the thermal inertia to promote wind power consumption.

(1) Heat Loss of the Building Envelope

[0029] The real heat transfer process is very complex and comprises convection, conduction and radiation. In order to simplify the heat transfer process, only the thermal loss calculation formula in steady state is used to solve the heat loss of building envelope.

Q.sub.t,k.sup.BE=A.sub.kF.sub.e,k(.sub.t,k.sup.in.sub.t,k.sup.out)(1+x.sub.k.sup.g)(1+x.sub.k.sup.ch) t,kK(16)

.sub.min.sup.in.sub.t,k.sup.in.sub.max.sup.in(17)

.sub.ch.sub.t+1,k.sup.in.sub.t,k+1.sup.in.sub.ch(18)

wherein Q.sub.t,k.sup.BB is the heat loss of the building envelope of the building, x.sub.k.sup.g and x.sub.k.sup.ch are respectively a correction coefficient of building height and a correction coefficient of building orientation in the k.sup.th heat load area, A.sub.k represents for the area of the building envelope of the building in the k.sup.th heat load area, F.sub.e,k represents for the heat transmission coefficient of the building envelope of the building in the k.sup.th heat load area, .sub.t,k.sup.in represents for the average indoor temperature in the k.sup.th heat load area at period t, .sub.t,k.sup.out represents for the outdoor temperature in the k.sup.th heat load area at period t, .sub.min.sup.in and .sub.max.sup.in respectively represent for minimum average indoor temperature and maximum average indoor temperature, .sub.ch represents for the maximum change rate of the average indoor temperature within the unit dispatching time, represents for the set of the dispatching periods, and K represents for a set of the heating areas of the heat loads;

(2) Cold Air Infiltration Heat Loss

[0030] Cold air infiltrates into a room through gaps such as doors and windows under the action of indoor and outdoor pressure difference and the cold air escapes after being heated, and the heat consumed by heating cold air from the outdoor temperature to the indoor temperature becomes cold air infiltration heat loss:

Q.sub.t,k.sup.Inf=0.278.Math.Num.Math.V.sub.kc.sub.p.sub.t.sup.out(.sub.t,k.sup.in.sub.t,k.sup.out) t,kK(19)

wherein Q.sub.t,k.sup.Inf is the cold air infiltration heat loss, Num represents for the air change rate per hour, V.sub.k represents for the total indoor volume in the k.sup.th heat load area, c.sub.p represents for constant-pressure specific heat capacity of cold air, and .sub.t.sup.out represents for outdoor air density at period t;

(3) Cold Air Invasion Heat Loss

[0031] The cold air invades into the room through opening external doors under the actions of air pressure and heat pressure in winter, and the heat consumed by heating the cold air to the indoor temperature is called as cold air invasion heat loss:

Q.sub.t,k.sup.Ven=0.278.Math.Ven.Math.c.sub.p.sub.t.sup.out(.sub.t,k.sup.in.sub.t,k.sup.out) t,kK(20)

wherein Q.sub.t,k.sup.Ven is the cold air invasion heat loss, and Ven represents for cold air invasion volume per hour;

(4) Calculation of Average Indoor Temperature Capable of Representing Thermal Inertia of Building

[0032] In order to simplify calculation, an assumption is made as follows: the air temperature in the building is consistent and can be described by .sub.t,k.sup.in; the thermal mass temperature distribution in the building is uniform, which means that the velocity of the surface heat diffusion process on the thermal mass is far larger than that of convection heat exchange; all heat obtained in the building within the unit time can be uniformly expressed by H.sub.t.sup.k Based on the assumption, the calculation of the average indoor temperature can be expressed by the formula (21):

Q.sub.t,k.sup.BE+Q.sub.t,k.sup.Inf+Q.sub.t,k.sup.Ven+C.sub.MM.sub.k(.sub.t+1,k.sup.in.sub.t,k.sup.in)/t=H.sub.t,k(21)

wherein H.sub.t,k is the average indoor temperature, C.sub.M represents for the specific heat capacity of the thermal mass, M.sub.k represents for the mass of the thermal mass in the k.sup.th heat load area, and the product of C.sub.M and M.sub.k can be obtained by an engineering experiment.

[0033] Step 4: a combined optimal dispatching model comprising conventional generators, wind power units, CHP units, electric boilers and heat storage tanks is established according to the heat source. The ultimate purpose of establishing the heating network model and the district heat consumer model is to consume wind power and reduce the combustion of the fossil fuel. Therefore, the combined optimal dispatching model comprising various power and heat supply devices is finally required to be established to check the effectiveness and practicability of the established heating network model and heat consumer model.

(1) CHP Unit Model

[0034] [00005]$\begin{array}{cc}0h_{i,t}{\stackrel{\_}{h}}_i.Math..Math..Math.i{I}^{\mathrm{CHP}},t& \left(22\right)\\ \{\begin{array}{c}{p}_{i,t}{\stackrel{\_}{p}}_i-C_i^1.Math.h_{i,t}\\ p_{i,t}\max .Math.\left\{C_i^{\max }.Math.h_{i,t}+C_i,{\underset{\_}{p}}_i-C_i^2.Math.h_{i,t}\right\}\end{array}.Math..Math..Math.i{I}^{\mathrm{CHP}},t& \left(23\right)\end{array}$

p.sub.i,t is generated output of the CHP units, p.sub.i is an upper limit of the generated output of the CHP units, p.sub.i is a lower limit of the generated output of the CHP units, C.sub.i.sup.max, C.sub.i.sup.1 and C.sub.i.sup.2 are heat and power coupling property factors of the CHP units, h.sub.i,t represents for heating power of the CHP units, h.sub.i represents for the maximum heating power, and I.sup.CHP represents for a set of serial numbers of CHP units;

(2) Electric Boiler Model

[0035]
0BH.sub.u,tBH.sub.u uU,t(24)

BH.sub.u,t=.Math.BP.sub.u,t uU,t(25)

BH.sub.u,t represents for the heating power of the electric boilers u at period t, BP.sub.u,t represents for electrical power consumed by the electric boilers, represents for the thermal efficiency of the electric boilers, BH.sub.u represents for the maximum heating power of the electric boilers, U represents for a set of serial numbers of the electric boilers, and represents for a set of dispatching periods;
(3) Establishment of Heat Storage Tank Model According to Same Total Heat Charge/Discharge of Heat Storage Tank within One Day

[00006]$\begin{array}{cc}{S}_{w,t+1}=S_{w,t}+{\mathrm{CS}}_{w,t}-{\mathrm{DCS}}_{w,t}.Math..Math..Math.wW,t& \left(26\right)\\ f_{w,t}\left\{0,1\right\}.Math..Math..Math.wW,t& \left(27\right)\\ 0S_{w,t}{\stackrel{\_}{S}}_w.Math..Math..Math.wW,t& \left(28\right)\\ 0{\mathrm{CS}}_{w,t}{f}_{w,t}.Math.{\stackrel{\_}{\mathrm{CS}}}_w.Math..Math..Math.wW,t& \left(29\right)\\ 0{\mathrm{DCS}}_{w,t}\left(1-f_{w,t}\right).Math.{\stackrel{\_}{\mathrm{DCS}}}_w.Math..Math..Math.wW,t& \left(30\right)\\ \underset{wW}{.Math.}.Math.\underset{t}{.Math.}.Math.{\mathrm{CS}}_{w,t}=\underset{wW}{.Math.}.Math.\underset{t}{.Math.}.Math.{\mathrm{DCS}}_{w,t}& \left(31\right)\end{array}$

S.sub.w,t represents for the heat reserve capacity of the heat storage tank w at period t, CS.sub.w,t represents for the heat charging power of the heat storage tank w at period t, DCS.sub.w,t represents for the heat discharging power of the heat storage tank w at period t, f.sub.w,t represents for the state of the heat storage tank w at period t, 0 represents for heat discharging of the heat storage tanks, 1 represents for heat storage of the heat storage tanks, S.sub.w represents for the maximum heat reserve capacity of the heat storage tank w, CS.sub.w represents for the maximum heat storage power of the heat storage tank w, DCS.sub.w represents for the maximum heat discharging power of the heat storage tank w, and W represents for a set of serial numbers of heat storage tanks; and the formula (31) shows that the total heat charge/discharge of the heat storage tanks within one day is the same;

(4) Electrical Power Balance

[0036] [00007]$\begin{array}{cc}\underset{i{I}^{\mathrm{CG}}}{.Math.}.Math.p_{i,t}+\underset{i{I}^{\mathrm{CHP}}}{.Math.}.Math.p_{i,t}+\underset{i{I}^{\mathrm{WP}}}{.Math.}.Math.p_{i,t}=\underset{uU}{.Math.}.Math.{\mathrm{BP}}_{u,t}+\underset{n{I}^{\mathrm{Load}}}{.Math.}.Math.{\mathrm{LD}}_{n,t}.Math..Math.t& \left(32\right)\end{array}$

I.sup.CG, I.sup.WP and I.sup.Load respectively represent for a serial number of the conventional generators, a serial number of wind farms and a serial number of electrical load nodes, and LD.sub.n,t represents for load power connecting to a bus n at period t;

(5) Heat Power Constraints

[0037] [00008]$\begin{array}{cc}\underset{uU}{.Math.}.Math.{\mathrm{BH}}_{u,t}+\underset{i{I}^{\mathrm{CHP}}}{.Math.}.Math.h_{i,t}\underset{kK}{.Math.}.Math.Q_{t+k,k}^{\mathrm{Total}}+\underset{wW}{.Math.}.Math.f_{w,t}.Math.{\mathrm{CS}}_{w,t}+\left(1-f_{w,t}\right).Math.{\mathrm{DCS}}_{w,t}& \left(33\right)\\ .Math.\begin{array}{cc}{Q}_{t+k,k}^{\mathrm{Total}}=Q_{t+k,k}^{\mathrm{BE}}+Q_{t+k,k}^{\mathrm{Inf}}+Q_{t+k,k}^{\mathrm{Ven}}& t,kK\end{array}& \left(34\right)\end{array}$

Q.sub.t+k,k.sup.Total is the total heat consumption power of the building in the k heating areas at period t+k;

(6) Ramp Rate Constraints of Unit

[0038]
rdown.sub.i.Math.tp.sub.i,1p.sub.i,t1rup.sub.i.Math.t

iI.sup.CHPUI.sup.CG,t(35)

rup.sub.i and rdown.sub.i respectively represent for the maximum upward ramp rate and downward ramp rate of the unit, and t is the dispatching time interval;

(7) Spinning Reserve Constraints of System

[0039] [00009]$\begin{array}{cc}\begin{array}{cc}\underset{i{I}^{\mathrm{CG}}}{.Math.}.Math.\min .Math.\left\{{\stackrel{\_}{p}}_i-p_{i,t},.Math..Math.t.Math.{\mathrm{rup}}_i\right\}{\mathrm{SR}}^{\mathrm{up}}& t\end{array}& \left(36\right)\\ \begin{array}{cc}\underset{i{I}^{\mathrm{CG}}}{.Math.}.Math.\min .Math.\left\{p_{i,t}-{\underset{\_}{p}}_i,.Math..Math.t.Math.{\mathrm{rdown}}_i\right\}{\mathrm{SR}}^{\mathrm{down}}& t\end{array}& \left(37\right)\end{array}$

SK.sub.up and SR.sup.down represent for the upward spinning reserve rate and downward spinning reserve rate required by a system, respectively;

(8) System Load Flow Constraints

[0040] [00010]$\begin{array}{cc}.Math.\underset{n{I}^{\mathrm{bus}}}{.Math.}.Math.{\mathrm{SF}}_{l,n}.Math.\left(\underset{i{S}_n^{\mathrm{CG}}.Math.{\mathrm{US}}_n^{\mathrm{CHP}}}{.Math.}.Math.p_{i,t}+\underset{i{S}_n^{\mathrm{WP}}}{.Math.}.Math.p_{i,t}-{\mathrm{LD}}_{n,t}\right).Math.F_l.Math..Math.l{I}^{\mathrm{Line}},t& \left(38\right)\end{array}$

wherein SF.sub.l,n represents for a transfer factor from the bus n to line l, F.sub.l represents for the maximum transmission capacity of the line l, I.sup.Line represents for the set of serial numbers of lines, S.sub.n.sup.CG represents for the set of indices of conventional generators, S.sub.n.sup.CHP represents for the set of indices of CHP cogeneration units, and S.sub.n.sup.WP represents for the set of indices of wind generation units.

[0041] The diagram of detailed hardware implementation is as shown in FIG. 2. Basic parameters, such as the number of residents, heating area and ambient temperature, are collected by data acquisition system. Regional geographic division system divides sub-heating areas according to the distance between heat loads and heat sources. Based on latest state evaluation of heating supply network and the satisfaction data fed back from users, the CHP optimal dispatch system takes responsibility for optimizing the proposed model in this invention.

[0042] All parameters used in a modeling process can be obtained through the existing technology.

[0043] In conclusion, a CHP optimal dispatching model is established by jointly using a plurality of sub-models. It is a mixed integer programming model essentially.

[0044] Part of heating system in the embodiment adopts the DHS in a certain city in northeast area in China, and the power system partially extracts actually operating power grid data in this province in order to realize the coordination of the capacities of power and heat supply system. Table 1 shows information of all devices including a conventional generator 1, a conventional generator 2, a CHP unit 1 and a CHP unit 2. Table 2 shows all basic parameters used in this calculation example. The program is compiled on an Matlab R2014b platform by using Yalmip language, and the problem of mixed integer programming established by using the program is solved by Cplex solver.

TABLE-US-00001 TABLE 1 TYPES AND CAPACITIES OF DEVICES INVOLVED IN DISPATCHING Unit Type Capacity Unit Type Capacity CG1 2 350 Electric boiler 2 140 MW MW CG2 3 350 Heat storage 1 1280 MW MW tank CHP1 2 350 Wind farm 700 MW MW CHP2 2 250 MW

TABLE-US-00002 TABLE 2 TYPES AND CAPACITIES OF DEVICES INVOLVED IN DISPATCHING Parameters Value Parameters Value .sub.e 2.4 W/m C. x.sub.k.sup.g 0.02 .sub.r 15 W/m.sup.2 C. x.sub.k.sup.ch 0.15 d.sub.in 2 m Num 0.1/h d.sub.ex 2.114 m Ven 15 m.sup.2/h T.sub.t.sup.surf 3 C. C.sub.M~M.sub.k 1 GJ/ C. .sub.b 0.023 W/m C. C.sub.p 1 KJ/kg C. 0.1 0.95 R.sub.b 0.384 m C./W .sub.min.sup.in 20 C. R.sub.e 0.049 m C./W .sub.max.sup.in 24 C. F.sub.e 0.44 W/m.sup.2 C. .sub.ch 1 C.

Embodiment 1

[0045] In the embodiment 1, system load flow constraints, the pressure loss of the heating pipe network and the heat storage tanks are not taken into consideration temporarily. The transmission time of hot water is far longer than that of electric power, and therefore, the farther distance between heat sources and heat load areas, the longer time delay it is necessary to be considered. In conventional optimal dispatching, the time delay of the heating network is not taken into consideration, and the output power of the heat source constantly follows the heat load, which makes the time period that the heat load reaches a peak value superposed with the time period that a power load reaches a peak value, and the device has to operate in the overload state to satisfy the demands of consumers, which is easy to cause system faults. In conventional optimal dispatching, heat load balance constraints are expressed by formulae (39)-(40):

[00011]$\begin{array}{cc}\begin{array}{cc}\underset{uU}{.Math.}.Math.{\mathrm{BH}}_{u,t}+\underset{i{I}^{\mathrm{CHP}}}{.Math.}.Math.h_{i,t}=\underset{kK}{.Math.}.Math.Q_{t,k}^{\mathrm{Total}}& t\end{array}& \left(39\right)\\ \begin{array}{cc}{Q}_{t,k}^{\mathrm{Total}}=Q_{t,k}^{\mathrm{BE}}+Q_{t,k}^{\mathrm{Inf}}+Q_{t,k}^{\mathrm{Ven}}& t,kK\end{array}& \left(40\right)\end{array}$

[0046] In the embodiment, the heat load area is divided into four parts, FIG. 3 respectively shows the variations of the average indoor temperature obtained from the conventional model and the model proposed by the present invention, and the two horizontal dotted lines in the figure respectively represent for the upper limit of the temperature and the lower limit of the temperature. As we can see, the further the distance between the heat source and heat load areas, the larger temperature fluctuation the inhabitants will suffer. In some periods, the average indoor temperature of the fourth load area is larger than upper limit. Heat users may feel too hot to open their windows to decrease the indoor temperature which cause a great waste of resource. In the actual situation, the indoor temperature is also often lower than the minimum limit where the heat load areas are far away from the heat source. This is one of the reasons why the heat-supply companies often receive complaints during heating periods. Fortunately, in most cases, heat users can be satisfied thanks to the strong thermal inertia of buildings.

[0047] FIG. 4 shows optimal results obtained by using the conventional optimal dispatching model and the optimal model proposed by the present invention, a slash filling part in FIG. 4 (a) represents for heat output regulation of the heat source, the slash filling part in FIG. 4 (b) represents for multi-adsorptive wind power, the slash filling part in FIG. 4 (c) represents for power output regulation of the heat supply unit, and the slash filling part in FIG. 4 (d) represents for spinning reserve capacity of the power system. In these two methods, the electric boilers work at the maximum power. The optimal results show that the generation cost of the conventional optimal dispatching model is $149 120, and the cost of the proposed model is $142 470. Wind power absorbed by the optimal model proposed by the present invention can be increased by 270 MWh within one day due to the adoption of peak staggered regulation by which the CHP units operate more flexibly. In the model proposed by the present invention, the heat source does not urgently follows up the fluctuation of the heat load, and the variation range of indoor temperature is sufficiently utilized. Due to the variation, the heat and power outputs of the CHP units can be reduced, so that more space is vacated to absorb wind power.

Embodiment 2

[0048] In the embodiment 2, a large heat storage tank is taken into consideration in the optimal dispatching model proposed by the present invention. The influences of relevant parameters of the heat storage tanks to wind power absorption are analyzed and shown in the embodiment. FIG. 5 shows indoor temperature variations and wind power forecast curves in three environments with different average temperatures.

[0049] FIG. 6 shows the influences of the initial heat storage capacity as well as heat charging/discharging rate of the heat storage tanks to wind power absorption. Known from FIG. 6 (a), the initial heat storage capacity of the heat storage tanks can increase the wind power absorption rate.

[0050] However, when the temperature follows T3 and the initial capacity is less than 300 MWh, there is no feasible solution to this optimal problem. This is because the heat sources cannot meet the basic heat demand at the beginning of the simulation. Besides, when the initial capacity of the heat storage tank increases to 1100 MWh, the abandoned wind power cannot be cut down any more, this is because the heat reserve capacity is sufficient for schedule throughout the day, and it is no help to facilitate wind power utilization by increasing the heat storage capacity only.

[0051] Generally, the larger charging/discharging rate is supportive for cutting down the amount of abandoned wind power. However, in this project, when the charging rate larger than 100 MW and discharging rate larger than 450 MW, the charging/discharging rate will do not improve the wind power consumption any more due to the spinning reserve constraints. Hence, the parameters of the tank should be optimized based on the size of DHS and composition of all units. This example demonstrates that electric power dispatch and heat supply is still interdependence, though the strong linkage between heat and power is decoupled by heat storage tank.