CN106527142A - CCHP (combined cooling, heating and power) system coordinated scheduling method under active power distribution network environment - Google Patents

Abstract

The invention discloses a CCHP (combined cooling, heating and power) system coordinated scheduling method under an active power distribution network environment. The coordinated scheduling method is applied when a CCHP system is accessed to the active power distribution network and operates to solve the coordination problem between power distribution network reconstruction optimization and CCHP system economical optimization targets. The coordinated scheduling strategy is used for carrying out system optimization on the whole system comprising CCHP system subsystems. Therefore, through the artificial intelligence technology, the method can provide flexible, safe and reliable techniques for the inner portion of each CCHP system, between the CCHP systems and between the CCHP system and the active power distribution network, thereby realizing flexible economical scheduling and ensuring safety operation of a power grid.

Description

A kind of cooling heating and power generation system coordinated scheduling method under active power distribution network environment

Technical field

The present invention relates to the cooling heating and power generation system coordinated scheduling method under a kind of active power distribution network environment, belongs to active and matches somebody with somebody Optimal dispatch field.

Background technology

With the continuous development of global economy, increasingly serious, the cleaning energy of the gradually increase and atmospheric pollution of energy crisis Source and regenerative resource receive greatly attention.Many countries have formulated carbon emission standard, it is therefore intended that reducing energy consumption Amount, improve generating efficiency, the discharge for reducing greenhouse gases and the development of encouragement clean energy resource.Although conventional fossil energy sources generating skill Art is ripe, single-machine capacity is big, operation stability is high, but efficiency of energy utilization is high, peak load regulation ability, pollutant emission The shortcomings of measuring big constrains its further development.Clean energy resource is divided into renewable and non-renewable two class, regenerative resource bag Wind energy, water energy, solar energy, geothermal energy, tide energy etc. is included, non-renewable clean energy resource includes natural gas, hydrogen etc..Renewable energy Source has the advantages that power generation energy resource low cost, discharge are few, but single-machine capacity is little, and it is main lacking that the energy has randomness and undulatory property Fall into.Natural gas capacity of power unit constantly increases, and peak-shaving capability is good, can not only be used for base load longtime running, and can double shift Peaking operation.

In order to solve the problems, such as the above, cooling heating and power generation system is proposed.Cooling heating and power generation system one kind is with supply of cooling, heating and electrical powers Equipment is core, comprising multiple distributed units (generating, energy storage etc.), there is the microgrid form of the polynary energy balance.With biography System coal, the fossil energy such as oil peter out, and environmental problem becomes increasingly conspicuous, used as effective group of distributed energy supply unit Form is knitted, supply of cooling, heating and electrical powers system is widely regarded as solution with primary energy ratio height, the little, high reliability of pollution The certainly effective scheme of energy problem.

Using traditional approach cooling heating and power generation system is incorporated to active power distribution network operation be optimized scheduling when, cool and thermal power system System mainly from active power distribution network power taking, without carrying out sale of electricity, the traffic order that active power distribution network passively receives, and is seldom joined In dispatching to power system global optimization.

The content of the invention

Present invention aim at there is provided the cooling heating and power generation system coordinated scheduling method under a kind of active power distribution network environment, Universal coordination optimization framework is used for reference, is a kind of coordinated scheduling side for considering active power distribution network and cooling heating and power generation system Method.When cooling heating and power generation system is incorporated to active power distribution network and runs, application coordination dispatching method solves power distribution network reconfiguration to the method Coordination problem between optimization and cooling heating and power generation system economic load dispatching optimization aim.The present invention proposes coordinated scheduling strategy, For carrying out system optimization to the total system of the subsystem comprising cooling heating and power generation system.Therefore, it is artificial herein by using Intellectual technology, can be between supply of cooling, heating and electrical powers subsystem internal, cooling heating and power generation system and cooling heating and power generation system with have Flexible, safe and reliable technology is provided between the power distribution network of source, flexible economic load dispatching is realized, and is ensured the safe operation of electrical network.

Cooling heating and power generation system coordinated scheduling method under a kind of active power distribution network environment of the present invention, by active power distribution network with System of the cooling heating and power generation system as a synthesis, with the operation of active power distribution network loss minimization and cooling heating and power generation system into This minimum is coordinated to the integrated system that active power distribution network environment and cooling heating and power generation system are constituted excellent as optimization aim Change, realize the coordinated scheduling of integrated system.

Cooling heating and power generation system coordinated scheduling method under a kind of active power distribution network environment of the present invention comprises the steps：

Step 1, is input into the data of active power distribution network and cooling heating and power generation system；

Step 2, set up respectively active power distribution network optimization object function with active power distribution network loss minimization as optimization aim, With the cooling heating and power generation system optimization object function of the minimum optimization aim of cooling heating and power generation system operating cost；

Step 3, sets up the optimization object function of the integrated system that active power distribution network environment and cooling heating and power generation system are constituted；

Step 4, according to the energy balance of integrated system, sets up each balancing the load equation of integrated system；

Step 5, according to active power distribution network and the service requirement of cooling heating and power generation system, sets up each operation of integrated system about Shu Fangcheng；

Step 6, the optimization object function and equilibrium equation, constraint equation according to integrated system, using particle cluster algorithm The optimal objective value of integrated system and the output of each equipment different periods in cooling heating and power generation system is solved, it is comprehensive so as to draw The optimal scheduling of syzygy system.

As the further prioritization scheme of the present invention, the input active power distribution network and supply of cooling, heating and electrical powers system described in step 1 The data of system include：The structural information of active power distribution network, line parameter circuit value, load data；It is each that cooling heating and power generation system is included The position of device parameter, gas turbine.

Used as the further prioritization scheme of the present invention, the cooling heating and power generation system described in step 1 includes gas turbine, wind Power generator, photovoltaic unit, donkey boiler, heating coil, energy storage device, waste-heat recovery device, electric refrigerating machine, absorption refrigeration Machine.

Used as the further prioritization scheme of the present invention, the active power distribution network optimization object function described in step 2 is：

In formula, f₁For the network loss of active power distribution network；Hop count when T is dispatching cycle；Branch road sums of the N for active power distribution network；R_i For the resistance of i-th branch road；P_i,t、Q_i,tAnd U_i,tThe respectively active power of i-th branch road of t period active power distribution networks, idle Power and branch road head end voltage amplitude.

Used as the further prioritization scheme of the present invention, the cooling heating and power generation system optimization object function described in step 3 is：

min f₂=C_ELC+C_F

In formula, f₂For the operating cost of cooling heating and power generation system；Hop count when T is dispatching cycle；C_ELCFor total power purchase expense；C_e For electricity price；Δ t is duration dispatching cycle；P_grid+,t、P_grid-,tRespectively t period cooling heating and power generation systems are from active power distribution network Upper power taking, the power of sale of electricity；C_FFor buying the expense of natural gas；C_fFor buying the unit calorific value price of natural gas；P_Gt,tFor t The generated output of gas turbine in period cooling heating and power generation system；Q_Boi,tFor donkey boiler in t period cooling heating and power generation systems Heat production power；η_g、η_bThe respectively efficiency of gas turbine, donkey boiler.

Used as the further prioritization scheme of the present invention, the optimization object function of the integrated system described in step 3 is：

F=min [α f₁+(1-α)f₂]

In formula, α is weight, and 0≤α≤1.

Used as the further prioritization scheme of the present invention, each balancing the load equation of the integrated system described in step 4 includes Electric, cold power, the equilibrium equation of thermal power, specially：

(1) electric load equilibrium equation：

P_grid,t+P_PV,t+P_WT,t+P_Gt,t=P_L,t+P_S,t

In formula, P_grid,t、P_PV,t、P_WT,tPower that respectively t period cooling heating and power generation systems are obtained from active power distribution network, The power of the power of photovoltaic unit, wind-driven generator；P_L,t、P_S,tThe respectively electrical load requirement of t period users, energy storage device Power；

(2) heating power balance equation：

Q_Hrs,t+Q_Boi,t+Q_L,t=Q_H,t+Q_ex,t

In formula, Q_Hrs,tFor the thermal power that t period waste-heat recovery devices are provided；Q_Boi,tFor the system of t period donkey boilers Thermal power；Q_L,tPower is heated for t period heating coils；Q_H,t、Q_ex,tRespectively the thermal load demands of t period users, The used heat that cooling heating and power generation system is excluded；

(3) cold power balance equation：

Q_Ac,t+Q_Ec,t=Q_C,t

In formula, Q_Ac,t、Q_Ec,t、Q_C,tThe respectively refrigeration work(of the refrigeration work consumption of t period Absorption Refrigerators, electric refrigerating machine Rate, the refrigeration duty demand of user.

Used as the further prioritization scheme of the present invention, each operation constraint equation of the integrated system described in step 5 includes：

(1) node voltage of active power distribution network and branch current constraint：

In formula, U_n,t、WithThe respectively voltage and voltage upper and lower limit of t n-th node of period；I_i,t、 The respectively electric current of t i-th branch road of period and its higher limit；

(2) wind-driven generator units limits：

0≤P_{WT, t}≤P_WT,on

In formula, P_WT,onFor the rated power of wind-driven generator, P_{WT, t}For the power of t period wind-driven generators；

(3) photovoltaic unit output constraint：

0≤P_{PV, t}≤P_PV,on

In formula, P_PV,onFor the rated power of photovoltaic unit, P_{PV, t}For the power of t period photovoltaic units；

(4) Gas Turbine Output constraint：

0≤P_Gt,t≤P_Gt,on

In formula, P_Gt,onFor the rated power of gas turbine, P_Gt,tFor the power of t period gas turbines；

(5) donkey boiler units limits：

0≤Q_Boi,t≤Q_Boi,on

In formula, Q_Boi,onFor the rated power of donkey boiler；Q_Boi,tFor the power of t period donkey boilers；

(6) electric refrigerating machine constraint：

0≤Q_Ec,t≤Q_Ec,on

In formula, Q_Ec,onFor the rated power of electric refrigerating machine；Q_Ec,tFor the power of t period electric refrigerating machines；

(7) absorption refrigerating machine constraint：

0≤Q_Ac,t≤Q_Ac,on

In formula, Q_Ac,onFor the rated power of absorption refrigerating machine；Q_Ac,tFor the power of t period absorption refrigerating machines；

(8) heating coil constraint：

0≤Q_L,t≤Q_L,on

In formula, Q_L,onFor the rated power of heating coil；Q_L,tFor the power of t period heating coils；

(9) energy storage device constraint：

|P_S,t|≤P_S,on

In formula, SOC^tFor the state-of-charge of t period energy storage devices；SOC、Respectively energy storage device state-of-charge Upper and lower limit；P_S,tFor the power of t period user energy storage devices；P_S,onFor the rated power of energy storage device.

The present invention adopts above technical scheme compared with prior art, with following technique effect：

1st, the present invention based on cooling heating and power generation system economic optimization, power system cleaning, Effec-tive Function can be promoted, it is real Power-balance inside existing cooling heating and power generation system；

2nd, cooling heating and power generation system can be considered as the present invention the single controllable of an access active power distribution network, so as to The grid-connected problem of internal cooling heating and power generation system is solved, the autonomous system with self management and control ability is formed；

3rd, the cooling heating and power generation system that can be applied to various scales of the invention is incorporated to active power distribution network operation, in force It is stable, it is economical optimal.

Description of the drawings

Fig. 1 is method of the present invention flow chart.

Fig. 2 is the integrated system schematic diagram of the embodiment of the present invention.

Specific embodiment

Embodiments of the present invention are described below in detail, the example of the embodiment is shown in the drawings, wherein ad initio Same or similar element or the element with same or like function are represented to same or similar label eventually.Below by ginseng The embodiment for examining Description of Drawings is exemplary, is only used for explaining the present invention, and is not construed as limiting the claims.

The present invention cooling heating and power generation system be by a typhoon power generator (i.e.：WT), a photovoltaic train unit is (i.e.： PV), a miniature gas turbine is (i.e.：) and the part such as household loads constitutes MT.

Technical scheme is further elaborated with reference to specific embodiment, certain typical distribution of the present embodiment Web frame is as shown in Fig. 2 specific implementation process is as shown in Figure 1.

Step 1：The data of input active power distribution network and cooling heating and power generation system, wherein, the data of input include：Power distribution network Structural information, load data, as shown in table 1；Each device parameter of cooling heating and power generation system, as shown in table 2.

The initial parameter of 1 33 node system of table

Each device parameter of 2 cooling heating and power generation system of table

Step 2：Row write the active power distribution network optimization object function with active power distribution network loss minimization as optimization aim：

In formula, f₁For the network loss of active power distribution network；Hop count when T is dispatching cycle；Branch road sums of the N for active power distribution network；R_i For the resistance of i-th branch road；P_i,t、Q_i,tAnd U_i,tThe respectively active power of i-th branch road of t period active power distribution networks, idle Power and branch road head end voltage amplitude.

Step 3：Row write excellent with the cooling heating and power generation system of the minimum optimization aim of the operating cost of cooling heating and power generation system Change object function：

min f₂=C_ELC+C_F

In formula, f₂For the operating cost of cooling heating and power generation system；Hop count when T is dispatching cycle.C_ELCFor total power purchase expense, Unit is unit, if selling electricity to electrical network, need to deduct from expense.C_eFor electricity price；Δ t is duration dispatching cycle；P_grid+,t、P_grid-,t The respectively power of t period cooling heating and power generation systems power taking, sale of electricity from active power distribution network；C_FFor buying taking for natural gas With unit is unit；C_fFor buying the unit calorific value price of natural gas；P_Gt,tFor gas turbine in t period cooling heating and power generation systems Generated output；Q_Boi,tFor the heat production power of donkey boiler in t period cooling heating and power generation systems；η_g、η_bRespectively combustion gas wheel The efficiency of machine, donkey boiler.

Step 4：Set up the optimization object function of the integrated system that active power distribution network environment and cooling heating and power generation system are constituted：

F=min [α f₁+(1-α)f₂]

In formula, α is weight, and 0≤α≤1, and in the present invention, value is 0.4.

Step 5：According to the energy balance of integrated system, row write each balancing the load equation of integrated system, including comprehensive system The electric load equilibrium equation of system and hot and cold balancing the load equation.

(1) electric load equilibrium equation：

P_grid,t+P_PV,t+P_WT,t+P_Gt,t=P_L,t+P_S,t

In formula, P_grid,t、P_PV,t、P_WT,tPower that respectively t period cooling heating and power generation systems are obtained from active power distribution network, The power of the power of photovoltaic unit, wind-driven generator；P_L,t、P_S,tThe respectively electrical load requirement of t period users, energy storage device Power；

(2) heating power balance equation：

Q_Hrs,t+Q_Boi,t+Q_L,t=Q_H,t+Q_ex,t

In formula, Q_Hrs,tFor the thermal power that t period waste-heat recovery devices are provided；Q_Boi,tFor the system of t period donkey boilers Thermal power；Q_L,tPower is heated for t period heating coils；Q_H,t、Q_ex,tRespectively the thermal load demands of t period users, The used heat that cooling heating and power generation system is excluded；

(3) cold power balance equation：

Q_Ac,t+Q_Ec,t=Q_C,t

In formula, Q_Ac,t、Q_Ec,t、Q_C,tThe respectively refrigeration work(of the refrigeration work consumption of t period Absorption Refrigerators, electric refrigerating machine Rate, the refrigeration duty demand of user.

Step 6：According to the service requirement of active power distribution network, row write node current, voltage constraint equation, and according to cold and hot The service requirement of chp system, row write each device constraints of cooling heating and power generation system.

(1) node voltage and branch current constraint：

In formula, U_n,t、WithThe respectively voltage and voltage upper and lower limit of t n-th node of period；I_i,t、 The respectively electric current of t i-th branch road of period and its higher limit；

(2) wind-driven generator units limits：

0≤P_{WT, t}≤P_WT,on

In formula, P_WT,onFor the rated power of wind-driven generator, P_{WT, t}For the power of t period wind-driven generators；

(3) photovoltaic unit output constraint：

0≤P_{PV, t}≤P_PV,on

In formula, P_PV,onFor the rated power of photovoltaic unit, P_{PV, t}For the power of t period photovoltaic units；

(4) Gas Turbine Output constraint：

0≤P_Gt,t≤P_Gt,on

In formula, P_Gt,onFor the rated power of gas turbine, P_Gt,tFor the power of t period gas turbines；

(5) donkey boiler units limits：

0≤Q_Boi,t≤Q_Boi,on

In formula, Q_Boi,onFor the rated power of donkey boiler；Q_Boi,tFor the power of t period donkey boilers；

(6) electric refrigerating machine constraint：

0≤Q_Ec,t≤Q_Ec,on

In formula, Q_Ec,onFor the rated power of electric refrigerating machine；Q_Ec,tFor the power of t period electric refrigerating machines；

(7) absorption refrigerating machine constraint：

0≤Q_Ac,t≤Q_Ac,on

In formula, Q_Ac,onFor the rated power of absorption refrigerating machine；Q_Ac,tFor the power of t period absorption refrigerating machines；

(8) heating coil constraint：

0≤Q_L,t≤Q_L,on

In formula, Q_L,onFor the rated power of heating coil；Q_L,tFor the power of t period heating coils；

(9) energy storage device constraint：

|P_S,t|≤P_S,o_n

In formula, SOC^tFor the state-of-charge of t period energy storage devices；SOC、Respectively energy storage device state-of-charge Upper and lower limit；P_S,tFor the power of t period user energy storage devices；P_S,onFor the rated power of energy storage device.

Step 7：Optimization object function and equilibrium equation, constraint equation according to integrated system, using particle cluster algorithm Solve the desired value of the optimum of integrated system and the output of each equipment different periods, wherein, one dispatching cycle duration Δ t =1h as shown in Tables 3 and 4, so as to draw the optimal scheduling of integrated system.

The optimal objective value of 3 day part of table

Period	f1/ unit	f2/kw	Period	f1/ unit	f2/kw
						1	6037.295	274	13	2322.978	285
2	5761.394	275	14	2608.162	281
						3	5455.036	277	15	2898.501	282
4	4815.92	278	16	3200.46	282
						5	4193.235	279	17	3419.503	281
6	3971.167	280	18	3606.095	279
						7	4364.633	279	19	3788.631	282
8	3610.152	279	20	4137.477	282
						9	3261.305	279	21	4849.47	282
10	3013.868	281	22	5412.891	278
						11	2445.289	282	23	5843.346	275
12	2232.844	284	24	6305.22	274

The output of 4 each equipment different periods of table

The above, the only specific embodiment in the present invention, but protection scope of the present invention is not limited thereto, and appoints What be familiar with the people of the technology disclosed herein technical scope in, it will be appreciated that the conversion expected or replacement, should all cover The present invention include within the scope of, therefore, protection scope of the present invention should be defined by the protection domain of claims.

Claims (8)

1. a kind of cooling heating and power generation system coordinated scheduling method under active power distribution network environment, it is characterised in that supply of cooling, heating and electrical powers System is connected to active power distribution network, minimum as excellent using the operating cost of active power distribution network loss minimization and cooling heating and power generation system Change target, optimization is coordinated to the integrated system that active power distribution network environment and cooling heating and power generation system are constituted, realize comprehensive system The coordinated scheduling of system；

Specific implementation step is as follows：

Step 1, is input into the data of active power distribution network and cooling heating and power generation system；

Step 2, sets up active power distribution network optimization object function with active power distribution network loss minimization as optimization aim, respectively with cold The cooling heating and power generation system optimization object function of the minimum optimization aim of cogeneration system operating cost；

Step 3, sets up the optimization object function of the integrated system that active power distribution network environment and cooling heating and power generation system are constituted；

Step 4, according to the energy balance of integrated system, sets up each balancing the load equation of integrated system；

Step 5, according to active power distribution network and the service requirement of cooling heating and power generation system, sets up each operation constraint side of integrated system Journey；

Step 6, the optimization object function and equilibrium equation, constraint equation according to integrated system, using PSO Algorithm Go out the optimal objective value of integrated system and the output of each equipment different periods in cooling heating and power generation system, so as to draw comprehensive system The optimal scheduling of system.

2. the cooling heating and power generation system coordinated scheduling method under a kind of active power distribution network environment according to claim 1, its It is characterised by, the data of input active power distribution network and cooling heating and power generation system described in step 1 include：The knot of active power distribution network Structure information, line parameter circuit value, load data；Each device parameter that cooling heating and power generation system is included, the position of gas turbine.

3. the cooling heating and power generation system coordinated scheduling method under a kind of active power distribution network environment according to claim 1, its It is characterised by, the cooling heating and power generation system described in step 1 includes gas turbine, wind-driven generator, photovoltaic unit, auxiliary pot Stove, heating coil, energy storage device, waste-heat recovery device, electric refrigerating machine, Absorption Refrigerator.

4. the cooling heating and power generation system coordinated scheduling method under a kind of active power distribution network environment according to claim 3, its It is characterised by, the active power distribution network optimization object function described in step 2 is：

$\min f_1=\underset{t=1}{\overset{T}{\Σ}}\underset{i=1}{\overset{N}{\Σ}}{R}_i\frac{P_{i,t}^2+Q_{i,t}^2}{U_{i,t}^2}$

In formula, f₁For the network loss of active power distribution network；Hop count when T is dispatching cycle；Branch road sums of the N for active power distribution network；R_iFor The resistance of i branch road；P_i,t、Q_i,tAnd U_i,tThe respectively active power of i-th branch road of t period active power distribution networks, reactive power With branch road head end voltage amplitude.

5. the cooling heating and power generation system coordinated scheduling method under a kind of active power distribution network environment according to claim 4, its It is characterised by, the cooling heating and power generation system optimization object function described in step 3 is：

min f₂=C_ELC+C_F

$C_{ELC}=C_e\underset{t=1}{\overset{T}{\Σ}}(P_{grid+,t}-P_{grid-,t})\mathrm{\Δ}t$

$C_F=C_f\underset{t=1}{\overset{T}{\Σ}}(\frac{P_{Gt,t}}{{\mathrm{\η}}_g}+\frac{Q_{Boi,t}}{{\mathrm{\η}}_b})\mathrm{\Δ}t$

In formula, f₂For the operating cost of cooling heating and power generation system；Hop count when T is dispatching cycle；C_ELCFor total power purchase expense；C_eFor electricity Valency；Δ t is duration dispatching cycle；P_grid+,t、P_grid-,tRespectively t period cooling heating and power generation systems are taken from active power distribution network Electricity, the power of sale of electricity；C_FFor buying the expense of natural gas；C_fFor buying the unit calorific value price of natural gas；P_Gt,tFor the t periods The generated output of gas turbine in cooling heating and power generation system；Q_Boi,tFor the product of donkey boiler in t period cooling heating and power generation systems Thermal power；η_g、η_bThe respectively efficiency of gas turbine, donkey boiler.

6. the cooling heating and power generation system coordinated scheduling method under a kind of active power distribution network environment according to claim 5, its It is characterised by, the optimization object function of the integrated system described in step 3 is：

F=min [α f₁+(1-α)f₂]

In formula, α is weight, and 0≤α≤1.

7. the cooling heating and power generation system coordinated scheduling method under a kind of active power distribution network environment according to claim 6, its It is characterised by, each balancing the load equation of the integrated system described in step 4 includes electric, cold power, the equilibrium equation of thermal power, Specially：

(1) electric load equilibrium equation：

P_grid,t+P_PV,t+P_WT,t+P_Gt,t=P_L,t+P_S,t

In formula, P_grid,t、P_PV,t、P_WT,tPower that respectively t period cooling heating and power generation systems are obtained from active power distribution network, photovoltaic The power of the power of unit, wind-driven generator；P_L,t、P_S,tThe respectively electrical load requirement of t period users, the work(of energy storage device Rate；

(2) heating power balance equation：

Q_Hrs,t+Q_Boi,t+Q_L,t=Q_H,t+Q_ex,t

In formula, Q_Hrs,tFor the thermal power that t period waste-heat recovery devices are provided；Q_Boi,tWork(is heated for t period donkey boilers Rate；Q_L,tPower is heated for t period heating coils；Q_H,t、Q_ex,tRespectively thermal load demands of t period users, cold and hot The used heat that chp system is excluded；

(3) cold power balance equation：

Q_Ac,t+Q_Ec,t=Q_C,t

In formula, Q_Ac,t、Q_Ec,t、Q_C,tRespectively the refrigeration work consumption of t period Absorption Refrigerators, the refrigeration work consumption of electric refrigerating machine, The refrigeration duty demand of user.

8. the cooling heating and power generation system coordinated scheduling method under a kind of active power distribution network environment according to claim 7, its It is characterised by, each operation constraint equation of the integrated system described in step 5 includes：

(1) node voltage of active power distribution network and branch current constraint：

$U_n^{\min }\≤U_{n,t}\≤U_n^{\max }$

$I_{i,t}\≤I_i^{\max }$

In formula, U_n,t、WithThe respectively voltage and voltage upper and lower limit of t n-th node of period；I_i,t、Respectively The electric current of t i-th branch road of period and its higher limit；

(2) wind-driven generator units limits：

0≤P_{WT, t}≤P_WT,on

In formula, P_WT,onFor the rated power of wind-driven generator, P_{WT, t}For the power of t period wind-driven generators；

(3) photovoltaic unit output constraint：

0≤P_{PV, t}≤P_PV,on

In formula, P_PV,onFor the rated power of photovoltaic unit, P_{PV, t}For the power of t period photovoltaic units；

(4) Gas Turbine Output constraint：

0≤P_Gt,t≤P_Gt,on

In formula, P_Gt,onFor the rated power of gas turbine, P_Gt,tFor the power of t period gas turbines；

(5) donkey boiler units limits：

0≤Q_Boi,t≤Q_Boi,on

In formula, Q_Boi,onFor the rated power of donkey boiler；Q_Boi,tFor the power of t period donkey boilers；

(6) electric refrigerating machine constraint：

0≤Q_Ec,t≤Q_Ec,on

In formula, Q_Ec,onFor the rated power of electric refrigerating machine；Q_Ec,tFor the power of t period electric refrigerating machines；

(7) absorption refrigerating machine constraint：

0≤Q_Ac,t≤Q_Ac,on

In formula, Q_Ac,onFor the rated power of absorption refrigerating machine；Q_Ac,tFor the power of t period absorption refrigerating machines；

(8) heating coil constraint：

0≤Q_L,t≤Q_L,on

In formula, Q_L,onFor the rated power of heating coil；Q_L,tFor the power of t period heating coils；

(9) energy storage device constraint：

|P_S,t|≤P_S,on

$\underset{\‾}{SOC}\≤{\mathrm{SOC}}^t\≤\stackrel{\‾}{SOC}$

In formula, SOC^tFor the state-of-charge of t period energy storage devices；SOC、Respectively energy storage device state-of-charge it is upper, Lower limit；P_S,tFor the power of t period user energy storage devices；P_S,onFor the rated power of energy storage device.