Optimal scheduling method for electricity-heat multi-energy flow system based on heat supply phasor model

Abstract

An optimal scheduling method of an electricity-heat multi-energy flow system based on a heat supply phasor model is provided. The method considers a mutual influence of the electricity-heat system, establishes a constraint equation of a heat supply system in a phasor form, considers dynamic characteristics of the heat supply system, and realizes an optimal scheduling of the electricity-heat multi-energy flow system.

Claims

1. An optimal scheduling method for an electricity-heat multi-energy flow system based on a heat supply phasor model, comprising: step 1: converting a load power of a heat supply system in the electricity-heat multi-energy flow system into a phasor form as: ${\overset{}{Q}}_{i, k}^{HL} = \frac{1}{K} v (k) {.Math.}_{w = - K + 1}^{K} q_{i}^{HL,_{w}} e^{- jk_{w}},$ $k = 0, .Math., K, i I^{HL},$ where $v (k) = {\begin{matrix} 1 & k = 1, .Math., K - 1 \\ 1 / 2 & k = 0, K \end{matrix},$ $= \frac{}{K},$ a superscript HI, represents a heat load identifier, {dot over (Q)}.sub.i,k.sup.HL represents a heat consumption power phasor when a frequency of an i.sup.th heat load in the heat supply system is k, q.sub.i.sup.HI,.sup.w represents a load power of the i.sup.th heat load in the heat supply system at a scheduling moment .sub.w, represents a fundamental wave frequency of a phasor, K represents an order of the fundamental wave frequency of the phasor, a value of K is equal to a number of scheduling moments and K=24 hours, and represents a scheduling time interval; step 2: setting constraint conditions of the heat supply system in the electricity-heat multi-energy flow system in the phasor form, comprising: step 2-1: a constraint equation of a heat loss of a heat supply network pipe of the heat supply system in the phasor form being: ${\overset{}{T}}_{e, b, k} = {\begin{matrix} e^{- \frac{_{b}}{C_{w} m_{b}} L_{b}} {\overset{}{T}}_{s, b, k} + 1 - e^{- \frac{_{b}}{C_{w} m_{b}} L_{b}} t_{a m}, k = 0 \\ e^{- \frac{_{b}}{C_{w} m_{b}} L_{b}} e^{- jk \frac{_{w} A_{b}}{m_{b}} L_{b}} {\overset{}{T}}_{s, b, k}, k = 1, .Math., K \end{matrix},$ $b B^{H},$ where, B.sup.H represents a set of numbers of pipes of the heat supply system, {dot over (T)}.sub.s,b,k represents a temperature phasor at a head end when a frequency of a b.sup.th pipe in the heat supply system is k, {dot over (T)}.sub.e,b,k represents a temperature phasor at a tail end when a frequency of the b.sup.th pipe in the heat supply system is k, t.sub.am represents an ambient temperature of the heat supply system, m.sub.b represents a flow of the b.sup.th pipe, L.sub.b represents a length of the b.sup.th pipe, C.sub.w represents a specific heat capacity of water, a value of the specific heat capacity is 4182 joule/(kilogram-degree centigrade), .sub.b represents a heat transfer coefficient per unit length of the b.sup.th pipe, and .sub.b is obtained from an energy management system of the electricity-heat coupled multi-energy flow system; step 2-2: a constraint equation of a temperature of a multi-pipe junction of the heat supply system in the phasor form being: $\underset{i S_{n}^{HS}}{.Math.} {\overset{}{Q}}_{i, k}^{H S} + C_{w} \underset{b S_{n}^{H, in}}{.Math.} m_{b} {\overset{}{T}}_{e, b, k} = \underset{i S_{n}^{HL}}{.Math.} {\overset{}{Q}}_{i, k}^{H L} + C_{w} \underset{b S_{n}^{H, out}}{.Math.} m_{b} {\overset{}{T}}_{n, k},$ $n N^{H},$ $k = 0, .Math., K,$ where, a superscript HS represents an identifier of a heat source in the heat supply system, {dot over (Q)}.sub.i,k.sup.HS represents a heat supply power phasor when a frequency of an i.sup.th heat source in the heat supply system is k, {dot over (T)}.sub.e,b,k represents a temperature phasor at a tail end when a frequency of a b.sup.th pipe in the heat supply system is k, C.sub.w represents a specific heat capacity of water, the value of the specific heat capacity is 4182 joule/(kilogram-degree centigrade), m.sub.b represents a flow of the b.sup.th pipe, {dot over (Q)}.sub.i,k.sup.HL represents a heat consumption power phasor when a frequency of an i.sup.th heat load in the heat supply system is k, {dot over (T)}.sub.n,k represents a temperature phasor when a frequency of an n.sup.th node in the heat supply system is k, S.sub.n.sup.HS represents a set of numbers of heat source nodes of the heat supply system, S.sub.n.sup.HL represents a set of numbers of heat load nodes of the heat supply system, S.sub.n.sup.H,in represents a set of numbers of pipes when a node at a tail end is the n.sup.th node in the heat supply system, and S.sub.n.sup.H,out represents a set of numbers of pipes when a node at a head end is the n.sup.th node in the heat supply system; step 2-3: a constraint equation of a temperature at a head end of the pipe of the heat supply system in the phasor form being:
{dot over (T)}.sub.s,b,k={dot over (T)}.sub.n,k,bS.sub.n.sup.H,out,nN.sup.H,k=0, . . . ,K, where, {dot over (T)}.sub.s,b,k represents a temperature phasor at a head end when a frequency of a b.sup.th pipe in the heat supply system is k, {dot over (T)}.sub.u,k represents a temperature phasor when a frequency of a n.sup.th node in the heat supply system is k, S.sub.n.sup.H,out represents a set of numbers of pipes when a node at a head end is the n.sup.th node in the heat supply system, and N.sup.H represents a set of numbers of nodes of the heat supply system; step 2-4: a constraint equation of a heat source phasor of the heat supply system with a Fourier inverse transform being:
q.sub.i.sup.HS,.sup.w=Re(.sub.k=0.sup.K{dot over (Q)}.sub.i,k.sup.HSe.sup.jk.sup.w),iI.sup.HS,w=1, . . . ,K, where, a superscript HS represents an identifier of a heat source in the heat supply system q.sub.i.sup.HS,.sup.w represents a heat supply power of an i.sup.th heat source in the heat supply system at a scheduling moment .sub.w, a function Re() represents taking a real part of a complex number, {dot over (Q)}.sub.i,k.sup.HS represents a heat supply power phasor when a frequency of the i.sup.th heat source in the heat supply system is k, and I.sup.HS represents a set of numbers of heat sources of the heat supply system; step 2-5: a constraint equation of temperature historical data of a node of the heat supply system being:
t.sub.n,his.sup..sup.w,his=Re(.sub.k=0.sup.K{dot over (T)}.sub.n,ke.sup.jk.sup.w,his),nN.sup.H,w=K+1, . . . ,0, where, t.sub.n,his.sup..sup.w represents a node temperature of an n.sup.th node in the heat supply system at a historical scheduling moment .sub.w,his, a function Re() represents taking a real part of a complex number, {dot over (T)}.sub.n,k represents a temperature phasor when a frequency of a n.sup.th node in the heat supply system is k, and N.sup.H represents a set of numbers of nodes of the heat supply system; step 2-6: a constraint equation of a limit of the temperature of the node of the heat supply system being:
t.sub.nRe(.sub.k=0.sup.K{dot over (T)}.sub.n,ke.sup.jk.sup.w)t.sub.n,nN.sup.H,w=1, . . . ,K, where, t.sub.n represents a lower temperature limit of an n.sup.th node in the heat supply system, t.sub.n represents a upper temperature limit of the n.sup.th node in the heat supply system, a function Re() represents taking a real part of a complex number, {dot over (T)}.sub.n,k represents a temperature phasor when a frequency of the n.sup.th node in the heat supply system is k, and N.sup.H represents a set of numbers of nodes of the heat supply system; step 2-7: constraint equations of a combined heat and power unit in the heat supply system being:
p.sub.i.sup.ES,.sup.w=.sub.,i.sup.NK.sup.i.sub.,i.sup..sup.wP.sub.,i,iI.sup.CHP,w=1, . . . ,K,
q.sub.i.sup.HS,.sup.w=.sub.,i.sup.NK.sup.i.sub.,i.sup..sup.wQ.sub.,i,iI.sup.CHP,w=1, . . . ,K,
.sub.=1.sup.NK.sup.i.sub.,i.sup..sup.wQ.sub.,i,iI.sup.CHP,w=1, . . . ,K,
.sub.,i.sup..sup.w0,=1, . . . ,NK.sub.i,iI.sup.CHP,w=1, . . . ,K, where, a superscript ES represents a power supply identifier, p.sub.i.sup.ES,.sup.w represents a power generation of an i.sup.th combined heat and power unit in the heat supply system at a scheduling moment .sub.w, q.sub.i.sup.HS,.sup.w represents a heat supply power of the i.sup.th combined heat and power unit in the heat supply system at the scheduling moment .sub.w, P.sub.,i represents an abscissa of a .sup.th vertex of an approximate polygon of an running feasible domain of the i.sup.th combined heat and power unit, Q.sub.,i represents an ordinate of the .sup.th vertex of an approximate polygon of the running feasible domain of the i.sup.th combined heat and power unit, .sub.,i.sup..sup.w represents a .sup.th combined coefficient of the i.sup.th combined heat and power unit at the scheduling moment .sub.w, NK.sub.i represents a number of vertexes of the approximate polygon of the running feasible domain of the i.sup.th combined heat and power unit, the approximate polygon of the running feasible domain of the combined heat and power unit is obtained from a factory specification of the combined heat and power unit, and I.sup.CHP represents a set of numbers of combined heat and power units in the heat supply system; step 3: setting constraint conditions of a power system in the electricity-heat multi-energy flow system, comprising: step 3-1: constraint equations of a direct current power flow of the power system being: $\underset{i I^{E S}}{.Math.} p_{i}^{ES,_{w}} = \underset{n N^{E}}{.Math.} p_{n}^{EL,_{w}},$ $w = 1, .Math., K,$ $- F_{b} \underset{n N^{E}}{.Math.}_{b, n} (\underset{i S_{n}^{E S}}{.Math.} p_{i}^{ES,_{w}} - p_{n}^{EL,_{w}}) F_{b},$ $b B^{E},$ $w = 1, .Math., K,$ where, a superscript ES represents a power supply identifier, p.sub.i.sup.ES,.sup.w represents a power generation of an i.sup.th generator unit in the power system at a scheduling moment .sub.w, p.sub.n.sup.EL.sup.w represents an electrical load power of an n.sup.th node in the heat supply system at a scheduling moment .sub.w, I.sup.ES represents a set of numbers of generator units of the power system, N.sup.E represents a set of numbers of nodes of the power system, F.sub.b represents an upper limit of a power of a b.sup.th line in the power system, .sub.b,n represents a transfer distribution factor between the n.sup.th node and the b.sup.th line in the power system, S.sub.n.sup.ES represents a set of generator units on the n.sup.th node in the power system, and B.sup.E represents a set of lines in the power system; step 3-2: a constraint equation of the generator unit in the power system being:
p.sub.ip.sub.i.sup.ES,.sup.wp.sub.i,iI.sup.TU,w=1, . . . ,K, where, a superscript.sup.TU represents an identifier of anther generator unit except the combined heat and power unit in the power system, p.sub.i represents a lower power limit of an i.sup.th generator unit in the power system, p.sub.i represents an upper power limit of the i.sup.th generator unit in the power system, p.sub.i.sup.ES,.sup.w represents a power of the i.sup.th generator unit in the power system at a scheduling moment .sub.w, and I.sup.TU represents a set of numbers of the generator units of the power system; step 4: establishing an objective function for optimal scheduling the electricity-heat multi-energy flow system, wherein the objective function is: ${.Math.}_{w = 1}^{K} (\underset{i I^{CHP}}{.Math.} c_{i}^{CHP,_{w}} + \underset{i I^{TU}}{.Math.} c_{i}^{TU,_{w}}),$ where, c.sub.i.sup.CHP,.sup.w represents an operation cost of an i.sup.th combined heat and power unit in the heat supply system at a scheduling moment .sub.w, c.sub.i.sup.TU,.sup.w represents an operation cost of an i.sup.th generator unit in the power system at the scheduling moment .sub.w, I.sup.CHP represents a set of numbers of combined heat and power units in the heat supply system, I.sup.TU represents a set of numbers of generator units in the power system, and specific expressions of c.sub.i.sup.CHP,.sup.w and c.sub.i.sup.TU,.sup.w are:
c.sub.i.sup.CHP,.sup.w=a.sub.0,i+a.sub.1,ip.sub.i.sup.ES,.sup.w+a.sub.2,iq.sub.i.sup.HS,.sup.w+a.sub.3,i(p.sub.i.sup.ES,.sup.w).sup.2+a.sub.4,i(q.sub.i.sup.HS,.sup.w).sup.2+a.sub.5,ip.sub.i.sup.ES,.sup.wq.sub.i.sup.HS,.sup.w,iI.sup.CHP,
c.sub.i.sup.TU,.sup.w=a.sub.0,i+a.sub.1,ip.sub.i.sup.ES,.sup.w+a.sub.3,i(p.sub.i.sup.ES,.sup.w).sup.2,iI.sup.TU, where, a.sub.0,i, a.sub.1,i, a.sub.2,i, a.sub.3,i, a.sub.4,i, and a.sub.5,i represents cost factors of an i.sup.th combined heat and power unit/generator unit, a.sub.0,i, a.sub.1,i, a.sub.2,i, a.sub.3,i, a.sub.4,i, and a.sub.5,i are obtained from the energy management system of the electricity-heat coupled multi-energy flow system, p.sub.i.sup.ES,.sup.w represents a power generation of the i.sup.th combined heat and power unit or generator unit at a scheduling moment .sub.w, and q.sub.i.sup.HS,.sup.w represents a heat supply power of the i.sup.th combined heat and power unit at the scheduling moment .sub.w; step 5: solving, by using an interior point method, an optimization model consisting of the objective function in step 4 and the constraint conditions in step 2 and step 3, obtaining the power generation of the generator unit, the power generation and the heat supply power of the combined heat and power unit in the electricity-heat multi-energy flow system as optimal scheduling parameters of the electricity-heat multi-energy flow system, to achieve the optimized scheduling of the electricity-heat multi-energy flow system based on the heat supply phasor model, and generating power by the generator unit in the power system in the electricity-heat multi-energy flow system based on the power generation of the generator unit obtained in step 5, and generating power and supplying heat by the combined heat and power unit in the heat supply system in the electricity-heat multi-energy flow system based on the power generation and the heat supply power of the combined heat and power unit obtained in step 5.

Description

DESCRIPTION OF EMBODIMENTS

(1) The present disclosure provides an optimal scheduling method of an electricity-heat multi-energy flow system based on a heat supply phasor model. The method includes the following steps.

(2) Step 1: a load power of a heat supply system in the electricity-heat multi-energy flow system is converted into a phasor form as:

(3) ${\overset{}{Q}}_{i, k}^{HL} = \frac{1}{K} v (k) {.Math.}_{w = - K + 1}^{K} q_{i}^{HL,_{w}} e^{- jk_{w}},$ $k = 0, .Math., K, i I^{HL},$ where

(4) $v (k) = {\begin{matrix} 1 & k = 1, .Math., K - 1 \\ 1 / 2 & k = 0, K \end{matrix},$ $= \frac{}{K},$
a superscript HL represents a heat load identifier, {dot over (Q)}.sub.i,k.sup.HL represents a heat consumption power phasor when a frequency of an i.sup.th heat load in the heat supply system is k, q.sub.i.sup.HL,.sup.w represents a load power of the i.sup.th heat load in the heat supply system at a scheduling moment .sub.w, represents a fundamental wave frequency of a phasor, K represents an order of the fundamental wave frequency of the phasor, a value of K is equal to a number of scheduling moments and K=24 hours, and represents a scheduling time interval.

(5) Step 2: constraint conditions of the heat supply system in the electricity-heat multi-energy flow system are set in the phasor form. Step 2 includes the followings.

(6) Step 2-1: a constraint equation of a heat loss of a heat supply network pipe of the heat supply system in the phasor form is

(7) ${\overset{}{T}}_{e, b, k} = {\begin{matrix} e^{- \frac{_{b}}{C_{w} m_{b}} L_{b}} {\overset{}{T}}_{s, b, k} + 1 - e^{- \frac{_{b}}{C_{w} m_{b}} L_{b}} t_{a m}, k = 0 \\ e^{- \frac{_{b}}{C_{w} m_{b}} L_{b}} e^{- jk \frac{_{w} A_{b}}{m_{b}} L_{b}} {\overset{}{T}}_{s, b, k}, k = 1, .Math., K \end{matrix},$ $b B^{H},$ where, B.sup.H represents a set of numbers of pipes of the heat supply system, {dot over (T)}.sub.s,b,k represents a temperature phasor at a head end when a frequency of a b.sup.th pipe in the heat supply system is k, {dot over (T)}.sub.e,b,k represents a temperature phasor at a tail end when a frequency of the b.sup.th pipe in the heat supply system is k, t.sub.am represents an ambient temperature of the heat supply system, m.sub.b represents a flow of the b.sup.th pipe, L.sub.b represents a length of the b.sup.th pipe, C.sub.w represents a specific heat capacity of water, a value of the specific heat capacity is 4182 joule/(kilogram-degree centigrade), .sub.b represents a heat transfer coefficient per unit length of the b.sup.th pipe, and .sub.b is obtained from an energy management system of the electricity-heat coupled multi-energy flow system.

(8) Step 2-2: a constraint equation of a temperature of a multi-pipe junction of the heat supply system in the phasor form is:

(9) 0 $\underset{i S_{n}^{HS}}{.Math.} {\overset{}{Q}}_{i, k}^{H S} + C_{w} \underset{b S_{n}^{H, in}}{.Math.} m_{b} {\overset{}{T}}_{e, b, k} = \underset{i S_{n}^{HL}}{.Math.} {\overset{}{Q}}_{i, k}^{H L} + C_{w} \underset{b S_{n}^{H, out}}{.Math.} m_{b} {\overset{}{T}}_{n, k},$ $n N^{H},$ $k = 0, .Math., K,$ where, a superscript HS represents an identifier of a heat source in the heat supply system, {dot over (Q)}.sub.i,k.sup.HS represents a heat supply power phasor when a frequency of an i.sup.th heat source in the heat supply system is k, {dot over (T)}.sub.e,b,k represents a temperature phasor at a tail end when a frequency of a b.sup.th pipe in the heat supply system is k, C.sub.w represents a specific heat capacity of water, the value of the specific heat capacity is 4182 joule/(kilogram-degree centigrade), m.sub.b represents a flow of the b.sup.th pipe, represents a heat consumption power phasor when a frequency of an i.sup.th heat load in the heat supply system is k, {dot over (T)}.sub.n,k represents a temperature phasor when a frequency of an nth node in the heat supply system is k, S.sub.n.sup.HS represents a set of numbers of heat source nodes of the heat supply system, S.sub.n.sup.HL represents a set of numbers of heat load nodes of the heat supply system, S.sub.n.sup.H,in represents a set of numbers of pipes when a node at a tail end is the nth node in the heat supply system, and S.sub.n.sup.H,out represents a set of numbers of pipes when a node at a head end is the nth node in the heat supply system.

(10) Step 2-3: a constraint equation of a temperature at a head end of the pipe of the heat supply system in the phasor form is:
{dot over (T)}.sub.s,b,k={dot over (T)}.sub.n,k,bS.sub.n.sup.H,out,nN.sup.H,k=0, . . . ,K where, {dot over (T)}.sub.s,b,k represents a temperature phasor at a head end when a frequency of a b.sup.th pipe in the heat supply system is k, {dot over (T)}.sub.n,k represents a temperature phasor when a frequency of a nth node in the heat supply system is k, S.sub.n.sup.H,out represents a set of numbers of pipes when a node at a head end is the nth node in the heat supply system, and N.sup.H represents a set of numbers of nodes of the heat supply system.

(11) Step 2-4: a constraint equation of a heat source phasor of the heat supply system with a Fourier inverse transform is:
q.sub.i.sup.HS,.sup.w=Re(.sub.k=0.sup.K{dot over (Q)}.sub.i,k.sup.HSe.sup.jk.sup.w),iI.sup.HS,w=1, . . . ,K. where, a superscript HS represents an identifier of a heat source in the heat supply system, q.sub.i.sup.HS,.sup.w represents a heat supply power of an i.sup.th heat source in the heat supply system at a scheduling moment .sub.w, a function Re(.circle-solid.) represents taking a real part of a complex number, {dot over (Q)}.sub.i,k.sup.HS represents a heat supply power phasor when a frequency of the i.sup.th heat source in the heat supply system is k, and I.sup.HS represents a set of numbers of heat sources of the heat supply system.

(12) Step 2-5: a constraint equation of temperature historical data of a node of the heat supply system is:
t.sub.n,his.sup..sup.w,his=Re(.sub.k=0.sup.K{dot over (T)}.sub.n,ke.sup.jk.sup.w,his),nN.sup.H,w=K+1, . . . ,0. where, t.sub.n,his.sup..sup.w represents a node temperature of an nth node in the heat supply system at a historical scheduling moment .sub.w,his, a function Re(.circle-solid.) represents taking a real part of a complex number, {dot over (T)}.sub.n,k represents a temperature phasor when a frequency of a nth node in the heat supply system is k, and N.sup.H represents a set of numbers of nodes of the heat supply system.

(13) Step 2-6: a constraint equation of a limit of the temperature of the node of the heat supply system is:
t.sub.nRe(.sub.k=0.sup.K{dot over (T)}.sub.n,ke.sup.jk.sup.w)t.sub.n,nN.sup.H,w=1, . . . ,K. where, t.sub.n represents a lower temperature limit of an nth node in the heat supply system, t.sub.n represents a upper temperature limit of the nth node in the heat supply system, a function Re(.circle-solid.) represents taking a real part of a complex number, {dot over (T)}.sub.n,k represents a temperature phasor when a frequency of the nth node in the heat supply system is k, and N.sup.H represents a set of numbers of nodes of the heat supply system.

(14) Step 2-7: constraint equations of a combined heat and power unit in the heat supply system is:
p.sub.i.sup.ES,.sup.w=.sub.,i.sup.NK.sup.i.sub.,i.sup..sup.wP.sub.,i,iI.sup.CHP,w=1, . . . ,K.
q.sub.i.sup.HS,.sup.w=.sub.,i.sup.NK.sup.i.sub.,i.sup..sup.wQ.sub.,i,iI.sup.CHP,w=1, . . . ,K.
.sub.=1.sup.NK.sup.i.sub.,i.sup..sup.wQ.sub.,i,iI.sup.CHP,w=1, . . . ,K.
.sub.,i.sup..sup.w0,=1, . . . ,NK.sub.i,iI.sup.CHP,w=1, . . . ,K. where, a superscript ES represents a power supply identifier, p.sub.i.sup.ES,.sup.w represents a power generation of an i.sup.th combined heat and power unit in the heat supply system at a scheduling moment .sub.w, q.sub.i.sup.HS,.sup.w represents a heat supply power of the i.sup.th combined heat and power unit in the heat supply system at the scheduling moment .sub.w, P.sub.,i represents an abscissa of a .sup.th vertex of an approximate polygon of an running feasible domain of the i.sup.th combined heat and power unit, Q.sub.,i represents an ordinate of the .sup.th vertex of an approximate polygon of the running feasible domain of the i.sup.th combined heat and power unit, .sub.,i.sup..sup.w represents a .sup.th combined coefficient of the i.sup.th combined heat and power unit at the scheduling moment .sub.w, NK.sub.i represents a number of vertexes of the approximate polygon of the running feasible domain of the i.sup.th combined heat and power unit, the approximate polygon of the running feasible domain of the combined heat and power unit is obtained from a factory specification of the combined heat and power unit, and I.sup.CHP represents a set of numbers of combined heat and power units in the heat supply system.

(15) Step 3: constraint conditions of a power system in the electricity-heat multi-energy flow system are set. Step 3 includes the followings.

(16) Step 3-1: constraint equations of a direct current power flow of the power system are:

(17) $\underset{i I^{E S}}{.Math.} p_{i}^{ES,_{w}} = \underset{n N^{E}}{.Math.} p_{n}^{EL,_{w}},$ $w = 1, .Math., K,$ $- F_{b} \underset{n N^{E}}{.Math.}_{b, n} (\underset{i S_{n}^{E S}}{.Math.} p_{i}^{ES,_{w}} - p_{n}^{EL,_{w}}) F_{b},$ $b B^{E},$ $w = 1, .Math., K,$ where, a superscript ES represents a power supply identifier, p.sub.i.sup.ES,.sup.w represents a power generation of an i.sup.th generator unit in the power system at a scheduling moment .sub.w, p.sub.n.sup.EL,.sup.w represents an electrical load power of an nth node in the heat supply system at a scheduling moment .sub.w, I.sup.ES represents a set of numbers of generator units of the power system, N E represents a set of numbers of nodes of the power system, F.sub.b represents an upper limit of a power of a b.sup.th line in the power system, .sub.b,n represents a transfer distribution factor between the nth node and the b.sup.th line in the power system, S.sub.n.sup.ES represents a set of generator units on the nth node in the power system, and B.sup.E represents a set of lines in the power system. Step 3-2: a constraint equation of the generator unit in the power system is:
p.sub.ip.sub.i.sup.ES,.sup.wp.sub.i,iI.sup.TU,w=1, . . . ,K. where, a superscript TU represents an identifier of anther generator unit except the combined heat and power unit in the power system, p.sub.i represents a lower power limit of an i.sup.th generator unit in the power system, p.sub.i represents an upper power limit of the i.sup.th generator unit in the power system, p.sub.i.sup.ES,.sup.w represents a power of the i.sup.th generator unit in the power system at a scheduling moment .sub.w, and I.sup.TU represents a set of numbers of the generator units of the power system. Step 4: an objective function for optimal scheduling the electricity-heat multi-energy flow system is established. The objective function is:

(18) ${.Math.}_{w = 1}^{K} (\underset{i I^{CHP}}{.Math.} c_{i}^{CHP,_{w}} + \underset{i I^{TU}}{.Math.} c_{i}^{TU,_{w}}),$ where, c.sub.i.sup.CHP,.sup.w represents an operation cost of an i.sup.th combined heat and power unit in the heat supply system at a scheduling moment .sub.w, c.sub.i.sup.TU,.sup.w represents an operation cost of an i.sup.th generator unit in the power system at the scheduling moment .sub.w, I.sup.CHP represents a set of numbers of combined heat and power units in the heat supply system, I.sup.TU represents a set of numbers of generator units in the power system, and specific expressions of c.sub.i.sup.CHP,.sup.w and c.sub.i.sup.TU,.sup.w are:
c.sub.i.sup.CHP,.sup.w=a.sub.0,i+a.sub.1,ip.sub.i.sup.ES,.sup.w+a.sub.2,iq.sub.i.sup.HS,.sup.w+a.sub.3,i(p.sub.i.sup.ES,.sup.w).sup.2+a.sub.4,i(q.sub.i.sup.HS,.sup.w).sup.2+a.sub.5,ip.sub.i.sup.ES,.sup.wq.sub.i.sup.HS,.sup.w,iI.sup.CHP.
c.sub.i.sup.TU,.sup.w=a.sub.0,i+a.sub.1,ip.sub.i.sup.ES,.sup.w+a.sub.3,i(p.sub.i.sup.ES,.sup.w).sup.2,iI.sup.TU. where, a.sub.0,i, a.sub.1,i, a.sub.2,i, a.sub.3,i, a.sub.4,i, and a.sub.5,i represents cost factors of an i.sup.th combined heat and power unit/generator unit, a.sub.0,i, a.sub.1,i, a.sub.2,i, a.sub.3,i, a.sub.4,i, and a.sub.5,i are obtained from the energy management system of the electricity-heat coupled multi-energy flow system, p.sub.i.sup.ES,.sup.w represents a power generation of the i.sup.th combined heat and power unit or generator unit at a scheduling moment .sub.w, and q.sub.i.sup.HS,.sup.w represents a heat supply power of the i.sup.th combined heat and power unit at the scheduling moment .sub.w.

(19) Step 5: an optimization model consisting of the objective function in step 4 and the constraint conditions in step 2 and step 3 is solved by using an interior point method, and the power generation of the generator unit, the power generation and the heat supply power of the combined heat and power unit in the electricity-heat multi-energy flow system are obtained as optimal scheduling parameters of the electricity-heat multi-energy flow system, to achieve the optimized scheduling of the electricity-heat multi-energy flow system based on the heat supply phasor model.

(20) In the step (5), an Interior Point Method used for solving the equation is a method for solving a linear programming or nonlinear convex optimization problem, and is also a known public technology in the technical field.

Optimal scheduling method for electricity-heat multi-energy flow system based on heat supply phasor model

Assignee

Inventors

Cpc classification

Classification Explorer

G05B13/0235

PHYSICS

Classification Explorer

G05B2219/2639

PHYSICS

Classification Explorer

G05B19/042

PHYSICS

International classification

Classification Explorer

G05B19/042

PHYSICS

Abstract

Claims

Description