CN107609684B - Comprehensive energy system economic optimization scheduling method based on micro-grid - Google Patents

Comprehensive energy system economic optimization scheduling method based on micro-grid Download PDF

Info

Publication number
CN107609684B
CN107609684B CN201710738641.0A CN201710738641A CN107609684B CN 107609684 B CN107609684 B CN 107609684B CN 201710738641 A CN201710738641 A CN 201710738641A CN 107609684 B CN107609684 B CN 107609684B
Authority
CN
China
Prior art keywords
energy
power
ice storage
constraint
time period
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710738641.0A
Other languages
Chinese (zh)
Other versions
CN107609684A (en
Inventor
董树锋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wanke Energy Technology Co ltd
Original Assignee
Wanke Energy Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wanke Energy Technology Co ltd filed Critical Wanke Energy Technology Co ltd
Priority to CN201710738641.0A priority Critical patent/CN107609684B/en
Publication of CN107609684A publication Critical patent/CN107609684A/en
Application granted granted Critical
Publication of CN107609684B publication Critical patent/CN107609684B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses a comprehensive energy system economic optimization scheduling method based on a micro-grid. It comprises the following steps: independently modeling energy production equipment, energy conversion equipment and energy storage equipment in a factory, and constructing an energy supply structure of a comprehensive energy system based on an energy exchange network; the ice storage air conditioning system serial and parallel working modes are considered, so that an economic optimization scheduling model of the comprehensive energy system is further perfected, and the optimization model is more in line with the actual requirements of the project; the method is characterized in that the minimum annual operation cost formed by operation maintenance cost, electricity purchasing cost, fuel cost and energy storage depreciation cost is taken as an optimization target, cold-hot electrical balance constraint, equipment operation constraint and energy storage equipment constraint are considered, the micro-grid is optimally scheduled, and self-optimization of a factory is realized. The invention has the beneficial effects that: the energy utilization efficiency of the user side is improved, the energy utilization cost of the user is reduced, the economic benefit and the energy utilization rate are improved, and the method is suitable for the industrial park comprehensive energy systems of different types.

Description

Comprehensive energy system economic optimization scheduling method based on micro-grid
Technical Field
The invention relates to the technical field of comprehensive energy and power demand response correlation, in particular to a comprehensive energy system economic optimization scheduling method based on a micro-grid.
Background
Integrated Energy Systems (IES) are the next generation of intelligent energy systems, so that the energy production, transmission, storage and use of the energy systems have systematic, integrated and refined operation and management. The comprehensive energy system is an important physical carrier of an energy internet and is a key for realizing technologies such as multi-energy complementation, energy gradient utilization and the like. The industrial park is a complex energy system mainly based on industrial load, comprises various energy production/utilization devices, has high requirement on power supply reliability, but has the problems of low energy utilization rate, unreasonable energy structure, large peak-valley power difference, environmental pollution and the like. From the energy consumption situation of each industry in China, industrial energy consumption dominates the energy consumption in China and accounts for about 70% of the total energy consumption of the whole society, so that energy consumption optimization management needs to be carried out on a factory, and the economic benefit and the energy utilization rate of the factory are improved.
Disclosure of Invention
The invention provides a comprehensive energy system economic optimization scheduling method based on a micro-grid, which is used for improving economic benefits and energy utilization rate and aims to overcome the defects in the prior art.
In order to achieve the purpose, the invention adopts the following technical scheme:
a comprehensive energy system economic optimization scheduling method based on a micro-grid comprises the following steps:
(1) independently modeling energy production equipment, energy conversion equipment and energy storage equipment in a factory, and constructing an energy supply structure of a comprehensive energy system based on an energy exchange network;
(2) the ice storage air conditioning system serial and parallel working modes are considered, so that an economic optimization scheduling model of the comprehensive energy system is further perfected, and the optimization model is more in line with the actual requirements of the project;
(3) the minimum annual operation cost formed by operation maintenance cost, electricity purchasing cost, fuel cost and energy storage depreciation cost is taken as an optimization target, cold-hot electrical balance constraint, equipment operation constraint and energy storage equipment constraint are considered, the micro-grid is optimally scheduled, and self-optimization of a factory is realized:
Min CATC＝COM+CES+CBW+CF
wherein, COMReferred to as operating maintenance costs, CESRefers to the cost of electricity purchase, CBWRefers to the fuel cost, CFWhich refers to the energy storage depreciation cost.
The method comprises the steps of firstly modeling energy production equipment, energy conversion equipment and energy storage equipment in a factory, and building an energy supply structure of the comprehensive energy system based on an energy exchange network. Based on the method, two working modes of series connection and parallel connection of the ice cold storage air conditioning system are considered, under the conditions of cold-heat-electricity balance constraint and multiple equipment constraint of equipment, the minimum annual operation cost of a user is taken as a target, a microgrid economic optimization scheduling model is constructed and considered, and self-optimization scheduling of a factory is realized. The method can be applied to different types of industrial park comprehensive energy systems. On one hand, the multi-energy coupling of cold, heat and electricity is considered, the multiple energy sources are cooperatively complemented, a user is guided to make a reasonable energy utilization scheme, the energy utilization efficiency of the user side is improved, the energy utilization cost of the user is reduced, and therefore the economic benefit and the energy utilization rate are improved; on the other hand, the working modes of different devices in a factory are considered, the economic optimization scheduling model of the comprehensive energy system is further perfected, and the control precision of the model is optimized.
Preferably, in step (1), the energy supply structure of the integrated energy system comprises the following:
(a) gas turbine
The gas turbine is a core device in a combined cooling heating and power system, and the electric power and the recovered thermal power of the gas turbine are as follows:
Figure GDA0003252979050000021
Figure GDA0003252979050000031
in the formula:
Figure GDA0003252979050000032
and
Figure GDA0003252979050000033
electric power and gas consumption power respectively output by the ith gas turbine in the time period t; lambda [ alpha ]gasIs the heat value of natural gas;
Figure GDA0003252979050000034
representing the thermal power output by the waste heat boiler;
Figure GDA0003252979050000035
and
Figure GDA0003252979050000036
respectively the power generation efficiency of the gas turbine and the heat recovery efficiency of the waste heat boiler;
(b) gas boiler
Figure GDA0003252979050000037
In the formula:
Figure GDA0003252979050000038
and
Figure GDA0003252979050000039
respectively being the ith gas boilerThermal power and gas consumption rate output by the section t;
Figure GDA00032529790500000310
the heat supply efficiency of the gas boiler is improved;
(c) photovoltaic unit
Figure GDA00032529790500000311
In the formula:
Figure GDA00032529790500000312
electric power output for the ith photovoltaic unit in a time period t;
Figure GDA00032529790500000313
solar panel efficiency; s is the area of the cell panel;
Figure GDA00032529790500000314
the illumination intensity of the ith photovoltaic unit in unit area;
(d) absorption refrigerator
Figure GDA00032529790500000315
In the formula:
Figure GDA00032529790500000316
the cooling power of the ith absorption refrigerator in the time period t;
Figure GDA00032529790500000317
the refrigeration efficiency of the absorption refrigerator;
Figure GDA00032529790500000318
thermal power output for the ith gas turbine during time period t;
(e) heat pump
Figure GDA00032529790500000319
In the formula:
Figure GDA00032529790500000320
and
Figure GDA00032529790500000321
the thermal power and the consumed electric power of the ith heat pump in the time period t respectively;
Figure GDA00032529790500000322
the heat supply efficiency of the heat pump;
(f) household air conditioner
The electric refrigeration/heating user uses the air conditioner to generate cold or heat by using the refrigerating machine under the condition of consuming electric energy:
Figure GDA0003252979050000041
Figure GDA0003252979050000042
in the formula:
Figure GDA0003252979050000043
and
Figure GDA0003252979050000044
electric power respectively representing cooling power, heating power, cooling and heating consumption of the ith household air conditioner in the time period t;
Figure GDA0003252979050000045
and
Figure GDA0003252979050000046
respectively representing the refrigeration energy efficiency ratio and the heating energy efficiency ratio of a household air conditioner;
(g) heat storage device
Figure GDA0003252979050000047
In the formula:
Figure GDA0003252979050000048
and
Figure GDA0003252979050000049
respectively representing the heat storage quantity, the heat storage power and the heat supply power of the ith heat storage device in a time period t;
Figure GDA00032529790500000410
is the self-loss factor of the thermal storage device;
Figure GDA00032529790500000411
and
Figure GDA00032529790500000412
respectively showing the heat storage efficiency and the heating efficiency of the cold storage device; t is the number of time periods, T is the unit period length,
Figure GDA00032529790500000413
(h) battery energy storage
Figure GDA00032529790500000414
In the formula:
Figure GDA00032529790500000415
and
Figure GDA00032529790500000416
respectively representing the stored energy, the charging power and the discharging power of the ith battery during the time t;
Figure GDA00032529790500000417
self-loss coefficient for stored energy;
Figure GDA00032529790500000418
and
Figure GDA00032529790500000419
respectively representing the charging efficiency and the discharging efficiency of the stored energy.
Preferably, in the step (2), the ice storage air conditioner refrigerates at night during the power utilization valley, stores cold energy by using a cold storage medium, and releases the cold energy at daytime during the power utilization peak to meet the cold supply demand of a factory, and according to the connection condition and the working mode of the refrigerator and the ice storage equipment, the ice storage air conditioning system can be divided into a parallel type and a serial type, and according to two working modes of serial connection and parallel connection of the ice storage air conditioning system, an economic optimization scheduling model of the comprehensive energy system is further perfected, so that the optimized scheduling model is more in accordance with the actual demand of the engineering, and specifically comprises the following two working modes:
(i) a parallel ice storage air conditioner based on a dual-working condition refrigerator comprises: the refrigerator and the ice storage tank of the parallel ice storage air conditioning system are in parallel connection in the system, wherein the refrigerator and the ice storage tank can jointly supply cold and can also independently supply cold load, and the refrigerator can simultaneously make and supply cold;
Figure GDA0003252979050000051
Figure GDA0003252979050000052
Figure GDA0003252979050000053
Figure GDA0003252979050000054
Figure GDA0003252979050000055
Figure GDA0003252979050000056
Figure GDA0003252979050000057
in the formula:
Figure GDA0003252979050000058
and
Figure GDA0003252979050000059
respectively representing the refrigerating power of the ith refrigerating machine and the ice storage tank in the time t;
Figure GDA00032529790500000510
and
Figure GDA00032529790500000511
respectively representing the maximum refrigerating power of the ith refrigerating machine and the ice storage tank;
Figure GDA00032529790500000512
and
Figure GDA00032529790500000513
respectively representing the electric power of the ith refrigerating machine and the ice storage tank in the time t;
Figure GDA00032529790500000514
and
Figure GDA00032529790500000515
respectively representing the maximum electric power of the ith refrigerating machine and the ice storage tank;
Figure GDA00032529790500000516
and
Figure GDA00032529790500000517
respectively representing the ith ice in the time period tThe total electric power, the maximum electric power and the refrigeration power of the cold accumulation air-conditioning system; t ismeltIndicating being in the ice-melting period, TrefIndicating the ice accumulation period, TmeltAnd TrefThe formula shows that the ice storage and the ice melting operation of the ice storage tank can not be carried out at the same time;
Figure GDA0003252979050000061
representing the refrigeration energy efficiency ratio of the refrigerator;
Figure GDA0003252979050000062
and
Figure GDA0003252979050000063
respectively representing the ice making energy efficiency ratio and the ice melting efficiency of the ice storage tank;
Figure GDA0003252979050000064
and
Figure GDA0003252979050000065
respectively representing the ice storage capacity of the ith ice storage tank time period t +1 and the time period t;
Figure GDA0003252979050000066
is the self-loss coefficient of the ice storage tank;
(ii) a series type ice storage air conditioner based on a dual-working condition refrigerator comprises: the refrigerating machine and the ice storage tank of the tandem type ice storage air conditioning system are in tandem position in the system, the cold quantity distribution of the refrigerating machine and the ice storage tank meets a certain proportional relation, and the proportional relation met by the cold quantity distribution of the refrigerating machine and the ice storage tank is mainly embodied in the following two stages:
(I) in the stage of ice storage, the refrigerating machine produces cold and stores the cold in the ice storage tank, at the moment, the ice storage tank does not participate in the refrigeration operation, the refrigerating machine participates in the refrigeration operation, and the ice making power of the ice storage tank
Figure GDA0003252979050000067
And cooling power of refrigerator
Figure GDA0003252979050000068
The relationship is as follows:
Figure GDA0003252979050000069
in the formula:
Figure GDA00032529790500000610
and
Figure GDA00032529790500000611
the temperature difference of the ice storage tank and the glycol at the inlet and the outlet of the refrigerator is t;
(II) in the cold supply stage, the ice storage tank and the refrigerator must be supplied with cold simultaneously, and the cold quantity distribution of the ice storage tank and the refrigerator meets a certain proportional relation:
Figure GDA00032529790500000612
Figure GDA00032529790500000613
in the formula: epsilons.iThe distribution coefficient of the cooling capacity of the ith ice storage air conditioning system.
Preferably, in the operation mode (i), the specific control variable of the ice storage air conditioning system is the flow rate of the circulating glycol passing through the ice storage tank and the refrigerating machine, and the cooling capacity of the ice storage tank and the refrigerating machine and the flow rate of the circulating glycol have the following relationship:
Figure GDA0003252979050000071
Figure GDA0003252979050000072
in the formula:
Figure GDA0003252979050000073
and
Figure GDA0003252979050000074
respectively representing the flow rates of the circulating glycol passing through the ice storage tank and the refrigerating machine in a time period t; cgly、ρglyAnd Δ TglyThe specific heat capacity, the liquid density and the return water supply temperature difference of the ethylene glycol solution are respectively;
Figure GDA0003252979050000075
to improve the cooling efficiency of the refrigerator.
Preferably, in step (3),
(A) the operation and maintenance cost is as follows:
Figure GDA0003252979050000076
in the formula: xiOM.iThe operating maintenance cost per unit output power of the equipment i;
Figure GDA0003252979050000077
represents the output power of the ith device during time period t;
(B) the electricity purchasing cost is as follows:
Figure GDA0003252979050000078
in the formula:
Figure GDA0003252979050000079
and
Figure GDA00032529790500000710
the electricity purchase price and the electricity purchase power of the time period t are respectively;
Figure GDA00032529790500000711
and
Figure GDA00032529790500000712
electricity selling prices andselling electricity power;
(C) fuel cost:
Figure GDA00032529790500000713
in the formula:
Figure GDA00032529790500000714
and
Figure GDA00032529790500000715
gas consumption rates of the ith gas turbine and the ith gas boiler, respectively, for a time period tth;
Figure GDA00032529790500000716
is the gas price;
(D) energy storage depreciation cost:
along with the deepening of the discharge depth, the number of times of charge and discharge of the stored energy of the battery is reduced, but the total amount of charge and discharge of the battery is basically unchanged, if the total amount of charge and discharge of the stored energy of the battery in the whole life cycle is constant, the depreciation cost of obtaining the accumulated discharge of the stored energy of the battery by 1kWh is as follows:
Figure GDA0003252979050000081
in the formula: cbat.repReplacement cost for stored energy, qlifetimeOutputting total quantity for the whole service life of the energy storage monomer;
the depreciation cost of stored energy is:
Figure GDA0003252979050000082
wherein:
Figure GDA0003252979050000083
the discharge power of the ith battery in the time period t is stored.
Preferably, in the step (3), the cooling, heating and power balance constraints include an electric power balance constraint, a heating power balance constraint and a cooling power balance constraint; the electric power balance constraint comprises an alternating current bus total load constraint, an alternating current-direct current converter efficiency constraint, a direct current bus total load constraint, a tie line constraint and a power purchasing and selling state constraint, and the specific constraint conditions are as follows:
the total load constraint of the alternating current bus:
Figure GDA0003252979050000084
in the formula:
Figure GDA0003252979050000085
an AC load for a time period t;
Figure GDA0003252979050000086
electric power for the ac-dc converter;
Figure GDA0003252979050000087
total electric power for the user air conditioner;
(II) efficiency constraint of the AC-DC converter:
Figure GDA0003252979050000088
in the formula:
Figure GDA0003252979050000089
the total load of the direct current bus is time t; etaA/DThe conversion efficiency from alternating current to direct current; etaD/AThe conversion efficiency from direct current to alternating current;
and (III) total load constraint of the direct current bus:
Figure GDA0003252979050000091
in the formula:
Figure GDA0003252979050000092
a direct current load for a time period t;
and (IV) connecting line constraint and electricity purchasing and selling state constraint:
Figure GDA0003252979050000093
Figure GDA0003252979050000094
Figure GDA0003252979050000095
in the formula:
Figure GDA0003252979050000096
and
Figure GDA0003252979050000097
respectively the upper power limit of electricity purchasing and electricity selling to the power grid;
Figure GDA0003252979050000098
and
Figure GDA0003252979050000099
respectively 0-1 state variables of purchasing and selling electricity in the time period t,
Figure GDA00032529790500000910
taking 1 indicates that electricity is purchased,
Figure GDA00032529790500000911
and 1 is taken to represent electricity sales, and the condition that electricity cannot be purchased at the same time is limited.
Preferably, the thermal power balance constraint is as follows:
Figure GDA00032529790500000912
Figure GDA00032529790500000913
in the formula:
Figure GDA00032529790500000914
and
Figure GDA00032529790500000915
space thermal load and hot water load of plant equipment, respectively.
Preferably, the constraint conditions of the cold power balance constraint are as follows:
Figure GDA00032529790500000916
in the formula:
Figure GDA00032529790500000917
is the cooling load.
Preferably, in step (3), the constraint conditions of the plant operation constraints are as follows:
Figure GDA00032529790500000918
Figure GDA00032529790500000919
in the formula:
Figure GDA00032529790500000920
and
Figure GDA00032529790500000921
respectively representing the input and output power of the device i in a time period t;
Figure GDA00032529790500000922
and
Figure GDA00032529790500000923
respectively representing the upper and lower limits of the output power of the device i in the time period t;
Figure GDA0003252979050000101
and
Figure GDA0003252979050000102
respectively representing the upper and lower limits of the input power of the device i during the time period t.
Preferably, in the step (3), the energy storage device constraint needs to satisfy an energy storage state constraint and an energy charging and discharging power constraint, and in order to ensure the continuity of scheduling, the energy storage of the energy storage device should be kept consistent before and after the scheduling period; the constraint conditions of the energy storage device constraint are as follows:
Figure GDA0003252979050000103
SL.i＝ST.i
Figure GDA0003252979050000104
Figure GDA0003252979050000105
Figure GDA0003252979050000106
wherein:
Figure GDA0003252979050000107
and
Figure GDA0003252979050000108
representing the maximum and minimum storage capacities of the energy storage device, respectively; sL.iAnd ST.iFor the beginning of energy storageStarting capacity and capacity at the end of the scheduling period;
Figure GDA0003252979050000109
and
Figure GDA00032529790500001010
representing the maximum charge and discharge power of the energy storage device, respectively;
Figure GDA00032529790500001011
and
Figure GDA00032529790500001012
respectively representing the energy storage device in a 0-1 state variable for charging and discharging energy during time period t,
Figure GDA00032529790500001013
taking 1 as the energy to be charged,
Figure GDA00032529790500001014
and 1 is taken to represent energy release, so that the equipment cannot be charged and released simultaneously.
The invention has the beneficial effects that: on one hand, the multi-energy coupling of cold, heat and electricity is considered, the cooperative complementation of various energy sources is realized, a user is guided to formulate a reasonable energy utilization scheme, the energy utilization efficiency of the user side is improved, the energy utilization cost of the user is reduced, and therefore the economic benefit and the energy utilization rate are improved; on the other hand, the working modes of different devices in a factory are considered, the economic optimization scheduling model of the comprehensive energy system is further perfected, and the control precision of the model is optimized; the method is suitable for different types of industrial park comprehensive energy systems.
Drawings
Fig. 1 is a schematic view of the structure of the microgrid according to the present invention.
Detailed Description
The invention is further described with reference to the following figures and detailed description.
In the embodiment shown in fig. 1, the integrated energy system comprises 4 energy forms of cold, heat, electricity and gas, the load in the system is various and the functional equipment is rich, and the main equipment is a micro gas turbine, a photovoltaic cell, a waste heat boiler, an absorption refrigerator, a household air conditioner, a gas boiler, a cell energy storage device, a heat energy storage device and a cold storage device. The system exchanges electric power with a public power grid through a centralized electric power bus, adopts an operation mechanism of 'spontaneous self-use and surplus network access', preferentially meets various local load requirements, and simultaneously allows surplus electric quantity to be transmitted to a power distribution system. Meanwhile, no gas is produced in the comprehensive energy system, and only one-way purchasing behavior exists between the comprehensive energy system and a gas company.
A comprehensive energy system economic optimization scheduling method based on a micro-grid comprises the following steps:
(1) independently modeling energy production equipment, energy conversion equipment and energy storage equipment in a factory, and constructing an energy supply structure of a comprehensive energy system based on an energy exchange network;
(a) gas turbine
The gas turbine is a core device in a combined cooling heating and power system, and the electric power and the recovered thermal power of the gas turbine are as follows:
Figure GDA0003252979050000111
Figure GDA0003252979050000112
in the formula:
Figure GDA0003252979050000113
and
Figure GDA0003252979050000114
electric power and gas consumption power respectively output by the ith gas turbine in the time period t; lambda [ alpha ]gasIs the heat value of natural gas;
Figure GDA0003252979050000115
representing the thermal power output by the waste heat boiler;
Figure GDA0003252979050000116
and
Figure GDA0003252979050000117
respectively the power generation efficiency of the gas turbine and the heat recovery efficiency of the waste heat boiler;
(b) gas boiler
Figure GDA0003252979050000121
In the formula:
Figure GDA0003252979050000122
and
Figure GDA0003252979050000123
respectively outputting thermal power and gas consumption rate of the ith gas boiler in the time period t;
Figure GDA0003252979050000124
the heat supply efficiency of the gas boiler is improved;
(c) photovoltaic unit
Figure GDA0003252979050000125
In the formula:
Figure GDA0003252979050000126
electric power output for the ith photovoltaic unit in a time period t;
Figure GDA0003252979050000127
solar panel efficiency; s is the area of the cell panel;
Figure GDA0003252979050000128
the illumination intensity of the ith photovoltaic unit in unit area;
(d) absorption refrigerator
Figure GDA0003252979050000129
In the formula:
Figure GDA00032529790500001210
the cooling power of the ith absorption refrigerator in the time period t;
Figure GDA00032529790500001211
the refrigeration efficiency of the absorption refrigerator;
Figure GDA00032529790500001212
thermal power output for the ith gas turbine during time period t;
(e) heat pump
Figure GDA00032529790500001213
In the formula:
Figure GDA00032529790500001214
and
Figure GDA00032529790500001215
the thermal power and the consumed electric power of the ith heat pump in the time period t respectively;
Figure GDA00032529790500001216
the heat supply efficiency of the heat pump;
(f) household air conditioner
The electric refrigeration/heating user uses the air conditioner to generate cold or heat by using the refrigerating machine under the condition of consuming electric energy:
Figure GDA00032529790500001217
Figure GDA00032529790500001218
in the formula:
Figure GDA00032529790500001219
and
Figure GDA00032529790500001220
electric power respectively representing cooling power, heating power, cooling and heating consumption of the ith household air conditioner in the time period t;
Figure GDA00032529790500001221
and
Figure GDA00032529790500001222
respectively representing the refrigeration energy efficiency ratio and the heating energy efficiency ratio of a household air conditioner;
(g) heat storage device
Figure GDA0003252979050000131
In the formula:
Figure GDA0003252979050000132
and
Figure GDA0003252979050000133
respectively representing the heat storage quantity, the heat storage power and the heat supply power of the ith heat storage device in a time period t;
Figure GDA0003252979050000134
is the self-loss factor of the thermal storage device;
Figure GDA0003252979050000135
and
Figure GDA0003252979050000136
respectively showing the heat storage efficiency and the heating efficiency of the cold storage device; t is the number of time periods, T is the unit period length,
Figure GDA0003252979050000137
(h) battery energy storage
Figure GDA0003252979050000138
In the formula:
Figure GDA0003252979050000139
and
Figure GDA00032529790500001310
respectively representing the stored energy, the charging power and the discharging power of the ith battery during the time t;
Figure GDA00032529790500001311
self-loss coefficient for stored energy;
Figure GDA00032529790500001312
and
Figure GDA00032529790500001313
respectively representing the charging efficiency and the discharging efficiency of the stored energy.
(2) The ice storage air conditioning system serial and parallel working modes are considered, so that an economic optimization scheduling model of the comprehensive energy system is further perfected, and the optimization model is more in line with the actual requirements of the project;
the ice cold-storage air conditioner refrigerates when the power consumption valley at night, utilize cold-storage medium to store cold volume, and release cold volume when the power consumption peak daytime, in order to satisfy the cooling demand of mill, according to connection condition and the mode of refrigerator and ice storage equipment, ice cold-storage air conditioning system can divide into parallel and tandem type two kinds, establish ties and parallelly connected two kinds of modes according to ice cold-storage air conditioning system, further perfect comprehensive energy system economic optimization scheduling model, make the scheduling model after the optimization more agree with the actual demand of engineering, specifically include following two modes:
(i) a parallel ice storage air conditioner based on a dual-working condition refrigerator comprises: the refrigerator and the ice storage tank of the parallel ice storage air conditioning system are in parallel connection in the system, wherein the refrigerator and the ice storage tank can jointly supply cold and can also independently supply cold load, and the refrigerator can simultaneously make and supply cold;
Figure GDA0003252979050000141
Figure GDA0003252979050000142
Figure GDA0003252979050000143
Figure GDA0003252979050000144
Figure GDA0003252979050000145
Figure GDA0003252979050000146
Figure GDA0003252979050000147
in the formula:
Figure GDA0003252979050000148
and
Figure GDA0003252979050000149
respectively representing the refrigerating power of the ith refrigerating machine and the ice storage tank in the time t;
Figure GDA00032529790500001410
and
Figure GDA00032529790500001411
respectively representing the maximum refrigerating power of the ith refrigerating machine and the ice storage tank;
Figure GDA00032529790500001412
and
Figure GDA00032529790500001413
respectively representing the electric power of the ith refrigerating machine and the ice storage tank in the time t;
Figure GDA00032529790500001414
and
Figure GDA00032529790500001415
respectively representing the maximum electric power of the ith refrigerating machine and the ice storage tank;
Figure GDA00032529790500001416
and
Figure GDA00032529790500001417
respectively representing the total electric power, the maximum electric power and the refrigerating power of the ith ice storage air-conditioning system in a time period tth; t ismeltIndicating being in the ice-melting period, TrefIndicating the ice accumulation period, TmeltAnd TrefThe formula shows that the ice storage and the ice melting operation of the ice storage tank can not be carried out at the same time;
Figure GDA00032529790500001418
representing the refrigeration energy efficiency ratio of the refrigerator;
Figure GDA00032529790500001419
and
Figure GDA00032529790500001420
respectively representing the ice making energy efficiency ratio and the ice melting efficiency of the ice storage tank;
Figure GDA00032529790500001421
and
Figure GDA00032529790500001422
respectively representing the ice storage capacity of the ith ice storage tank time period t +1 and the time period t;
Figure GDA00032529790500001423
is the self-loss coefficient of the ice storage tank;
in engineering practice, the specific control variable of the ice storage air conditioning system is the flow of circulating glycol passing through an ice storage tank and a refrigerating machine, and the cooling capacity of the ice storage tank and the refrigerating machine and the flow of circulating glycol have the following relationship:
Figure GDA0003252979050000151
Figure GDA0003252979050000152
in the formula:
Figure GDA0003252979050000153
and
Figure GDA0003252979050000154
respectively representing the flow rates of the circulating glycol passing through the ice storage tank and the refrigerating machine in a time period t; cgly、ρglyAnd Δ TglyThe specific heat capacity, the liquid density and the return water supply temperature difference of the ethylene glycol solution are respectively;
Figure GDA0003252979050000155
to improve the cooling efficiency of the refrigerator.
(ii) A series type ice storage air conditioner based on a dual-working condition refrigerator comprises: the refrigerating machine and the ice storage tank of the tandem type ice storage air conditioning system are in tandem position in the system, the cold quantity distribution of the refrigerating machine and the ice storage tank meets a certain proportional relation, and the proportional relation met by the cold quantity distribution of the refrigerating machine and the ice storage tank is mainly embodied in the following two stages:
(I) in the ice storage stage, the refrigerating machine produces cold and stores the cold in the ice storage tank, at the moment, the ice storage tank does not participate in the refrigeration operation, the refrigerating machine participates in the refrigeration operation, and the ice storage tank stores the coldIce making power of ice tank
Figure GDA0003252979050000156
And cooling power of refrigerator
Figure GDA0003252979050000157
The relationship is as follows:
Figure GDA0003252979050000158
in the formula:
Figure GDA0003252979050000159
and
Figure GDA00032529790500001510
the temperature difference of the ice storage tank and the glycol at the inlet and the outlet of the refrigerator is t;
(II) in the cold supply stage, the ice storage tank and the refrigerator must be supplied with cold simultaneously, and the cold quantity distribution of the ice storage tank and the refrigerator meets a certain proportional relation:
Figure GDA00032529790500001511
Figure GDA00032529790500001512
in the formula: epsilons.iThe distribution coefficient of the cooling capacity of the ith ice storage air conditioning system.
(3) The minimum annual operation cost formed by operation maintenance cost, electricity purchasing cost, fuel cost and energy storage depreciation cost is taken as an optimization target, cold-hot electrical balance constraint, equipment operation constraint and energy storage equipment constraint are considered, the micro-grid is optimally scheduled, and self-optimization of a factory is realized:
Min CATC＝COM+CES+CBW+CF
wherein, COMReferred to as operating maintenance costs, CESRefers to the cost of electricity purchase, CBWRefers to the fuel cost, CFWhich refers to the energy storage depreciation cost.
(A) The operation and maintenance cost is as follows:
Figure GDA0003252979050000161
in the formula: xiOM.iThe operating maintenance cost per unit output power of the equipment i;
Figure GDA0003252979050000162
represents the output power of the ith device during time period t;
(B) the electricity purchasing cost is as follows:
Figure GDA0003252979050000163
in the formula:
Figure GDA0003252979050000164
and
Figure GDA0003252979050000165
the electricity purchase price and the electricity purchase power of the time period t are respectively;
Figure GDA0003252979050000166
and
Figure GDA0003252979050000167
the price and power of selling electricity in time t;
(C) fuel cost:
Figure GDA0003252979050000168
in the formula:
Figure GDA0003252979050000169
and
Figure GDA00032529790500001610
gas consumption rates of the ith gas turbine and the ith gas boiler, respectively, for a time period tth;
Figure GDA00032529790500001611
is the gas price;
(D) energy storage depreciation cost:
along with the deepening of the discharge depth, the number of times of charge and discharge of the stored energy of the battery is reduced, but the total amount of charge and discharge of the battery is basically unchanged, if the total amount of charge and discharge of the stored energy of the battery in the whole life cycle is constant, the depreciation cost of obtaining the accumulated discharge of the stored energy of the battery by 1kWh is as follows:
Figure GDA0003252979050000171
in the formula: cbat.repReplacement cost for stored energy, qlifetimeOutputting total quantity for the whole service life of the energy storage monomer; the depreciation cost of stored energy is:
Figure GDA0003252979050000172
wherein:
Figure GDA0003252979050000173
the discharge power of the ith battery in the time period t is stored.
1) The cold and hot electrical balance constraints include an electrical power balance constraint, a thermal power balance constraint, and a cold power balance constraint.
i) Electric power balance constraint:
the method comprises the following steps of AC bus total load constraint, AC-DC converter efficiency constraint, DC bus total load constraint, tie line constraint and electricity purchasing and selling state constraint, wherein the specific constraint conditions are as follows:
the total load constraint of the alternating current bus:
Figure GDA0003252979050000174
in the formula:
Figure GDA0003252979050000175
an AC load for a time period t;
Figure GDA0003252979050000176
electric power for the ac-dc converter;
Figure GDA0003252979050000177
total electric power for the user air conditioner;
(II) efficiency constraint of the AC-DC converter:
Figure GDA0003252979050000178
in the formula:
Figure GDA0003252979050000181
the total load of the direct current bus is time t; etaA/DThe conversion efficiency from alternating current to direct current; etaD/AThe conversion efficiency from direct current to alternating current;
and (III) total load constraint of the direct current bus:
Figure GDA0003252979050000182
in the formula:
Figure GDA0003252979050000183
a direct current load for a time period t;
and (IV) connecting line constraint and electricity purchasing and selling state constraint:
Figure GDA0003252979050000184
Figure GDA0003252979050000185
Figure GDA0003252979050000186
in the formula:
Figure GDA0003252979050000187
and
Figure GDA0003252979050000188
respectively the upper power limit of electricity purchasing and electricity selling to the power grid;
Figure GDA0003252979050000189
and
Figure GDA00032529790500001810
respectively 0-1 state variables of purchasing and selling electricity in the time period t,
Figure GDA00032529790500001811
taking 1 indicates that electricity is purchased,
Figure GDA00032529790500001812
and 1 is taken to represent electricity sales, and the condition that electricity cannot be purchased at the same time is limited.
ii) the thermal power balance constraint is as follows:
Figure GDA00032529790500001813
Figure GDA00032529790500001814
in the formula:
Figure GDA00032529790500001815
and
Figure GDA00032529790500001816
space thermal load and hot water load of plant equipment, respectively.
iii) the constraints of the cold power balance constraint are as follows:
Figure GDA00032529790500001817
in the formula:
Figure GDA00032529790500001818
is the cooling load.
2) And (4) equipment operation constraint:
Figure GDA00032529790500001819
Figure GDA0003252979050000191
in the formula:
Figure GDA0003252979050000192
and
Figure GDA0003252979050000193
respectively representing the input and output power of the device i in a time period t;
Figure GDA0003252979050000194
and
Figure GDA0003252979050000195
respectively representing the upper and lower limits of the output power of the device i in the time period t;
Figure GDA0003252979050000196
and
Figure GDA0003252979050000197
respectively representing the upper and lower limits of the input power of the device i during the time period t.
3) Energy storage equipment restraint:
the energy storage state constraint and the energy charging and discharging power constraint are required to be met, and in order to ensure the scheduling continuity, the energy storage of the energy storage equipment is kept consistent before and after the scheduling period; the constraints of the energy storage device constraints are as follows:
Figure GDA0003252979050000198
SL.i＝ST.i
Figure GDA0003252979050000199
Figure GDA00032529790500001910
Figure GDA00032529790500001911
wherein:
Figure GDA00032529790500001912
and
Figure GDA00032529790500001913
representing the maximum and minimum storage capacities of the energy storage device, respectively; sL.iAnd ST.iThe initial capacity of the stored energy and the capacity at the end of the scheduling period;
Figure GDA00032529790500001914
and
Figure GDA00032529790500001915
representing the maximum charge and discharge power of the energy storage device, respectively;
Figure GDA00032529790500001916
and
Figure GDA00032529790500001917
respectively representing the energy storage device in a 0-1 state variable for charging and discharging energy during time period t,
Figure GDA00032529790500001918
taking 1 as the energy to be charged,
Figure GDA00032529790500001919
and 1 is taken to represent energy release, so that the equipment cannot be charged and released simultaneously.
And according to the optimization result, outputting a self-optimization-trending energy utilization scheme of the factory, and reducing the operation cost of industrial users by adjusting the operation mode and the working state of each device in the system. The method comprises the steps of firstly modeling energy production equipment, energy conversion equipment and energy storage equipment in a factory, and building an energy supply structure of the comprehensive energy system based on an energy exchange network. Based on the method, two working modes of series connection and parallel connection of the ice cold storage air conditioning system are considered, under the conditions of cold-heat-electricity balance constraint and multiple equipment constraint of equipment, the minimum annual operation cost of a user is taken as a target, a microgrid economic optimization scheduling model is constructed and considered, and self-optimization scheduling of a factory is realized. The method can be applied to different types of industrial park comprehensive energy systems. On one hand, the multi-energy coupling of cold, heat and electricity is considered, the multiple energy sources are cooperatively complemented, a user is guided to make a reasonable energy utilization scheme, the energy utilization efficiency of the user side is improved, the energy utilization cost of the user is reduced, and therefore the economic benefit and the energy utilization rate are improved; on the other hand, the working modes of different devices in a factory are considered, the economic optimization scheduling model of the comprehensive energy system is further perfected, and the control precision of the model is optimized.

Claims (7)

1. A comprehensive energy system economic optimization scheduling method based on a micro-grid is characterized by comprising the following steps:
(1) independently modeling energy production equipment, energy conversion equipment and energy storage equipment in a factory, and constructing an energy supply structure of a comprehensive energy system based on an energy exchange network;
(2) the ice storage air conditioning system serial and parallel working modes are considered, so that an economic optimization scheduling model of the comprehensive energy system is further perfected, and the optimization model is more in line with the actual requirements of the project;
(3) the minimum annual operation cost formed by operation maintenance cost, electricity purchasing cost, fuel cost and energy storage depreciation cost is taken as an optimization target, cold-hot electrical balance constraint, equipment operation constraint and energy storage equipment constraint are considered, the micro-grid is optimally scheduled, and self-optimization of a factory is realized:
Min CATC＝COM+CES+CBW+CF
wherein, COMReferred to as operating maintenance costs, CESRefers to the cost of electricity purchase, CBWRefers to the fuel cost, CFRefers to energy storage depreciation cost;
wherein, in the step (1), the energy supply structure of the integrated energy system comprises the following steps:
(a) gas turbine
The gas turbine is a core device in a combined cooling heating and power system, and the electric power and the recovered thermal power of the gas turbine are as follows:
Figure FDA0003252979040000011
Figure FDA0003252979040000012
in the formula:
Figure FDA0003252979040000013
and
Figure FDA0003252979040000014
electric power and gas consumption power respectively output by the ith gas turbine in the time period t; lambda [ alpha ]gasIs the heat value of natural gas;
Figure FDA0003252979040000015
representing the thermal power output by the waste heat boiler;
Figure FDA0003252979040000016
and
Figure FDA0003252979040000017
respectively the power generation efficiency of the gas turbine and the heat recovery efficiency of the waste heat boiler;
(b) gas boiler
Figure FDA0003252979040000018
In the formula:
Figure FDA0003252979040000021
and
Figure FDA0003252979040000022
respectively outputting thermal power and gas consumption rate of the ith gas boiler in the time period t;
Figure FDA0003252979040000023
the heat supply efficiency of the gas boiler is improved;
(c) photovoltaic unit
Figure FDA0003252979040000024
In the formula:
Figure FDA0003252979040000025
electric power output for the ith photovoltaic unit in a time period t;
Figure FDA0003252979040000026
solar panel efficiency; s is the area of the cell panel;
Figure FDA0003252979040000027
the illumination intensity of the ith photovoltaic unit in unit area;
(d) absorption refrigerator
Figure FDA0003252979040000028
In the formula:
Figure FDA0003252979040000029
the cooling power of the ith absorption refrigerator in the time period t;
Figure FDA00032529790400000210
the refrigeration efficiency of the absorption refrigerator;
Figure FDA00032529790400000211
thermal power output for the ith gas turbine during time period t;
(e) heat pump
Figure FDA00032529790400000212
In the formula:
Figure FDA00032529790400000213
and
Figure FDA00032529790400000214
the thermal power and the consumed electric power of the ith heat pump in the time period t respectively;
Figure FDA00032529790400000215
the heat supply efficiency of the heat pump;
(f) household air conditioner
The electric refrigeration/heating user uses the air conditioner to generate cold or heat by using the refrigerating machine under the condition of consuming electric energy:
Figure FDA00032529790400000216
Figure FDA00032529790400000217
in the formula:
Figure FDA00032529790400000218
and
Figure FDA00032529790400000219
electric power respectively representing cooling power, heating power, cooling and heating consumption of the ith household air conditioner in the time period t;
Figure FDA00032529790400000220
and
Figure FDA00032529790400000221
respectively representing the refrigeration energy efficiency ratio and the heating energy efficiency ratio of a household air conditioner;
(g) heat storage device
Figure FDA0003252979040000031
In the formula:
Figure FDA0003252979040000032
and
Figure FDA0003252979040000033
respectively representing the heat storage quantity, the heat storage power and the heat supply power of the ith heat storage device in a time period t;
Figure FDA0003252979040000034
is the self-loss factor of the thermal storage device;
Figure FDA0003252979040000035
and
Figure FDA0003252979040000036
respectively showing the heat storage efficiency and the heating efficiency of the cold storage device; t is the number of time periods, T is the unit period length,
Figure FDA0003252979040000037
(h) battery energy storage
Figure FDA0003252979040000038
In the formula:
Figure FDA0003252979040000039
and
Figure FDA00032529790400000310
respectively representing the stored energy, the charging power and the discharging power of the ith battery during the time t;
Figure FDA00032529790400000311
self-loss coefficient for stored energy;
Figure FDA00032529790400000312
and
Figure FDA00032529790400000313
respectively representing the charging efficiency and the discharging efficiency of the stored energy;
in step (2), the ice storage air conditioner refrigerates when the power consumption is low at night, utilize the cold-storage medium to store cold volume, and release cold volume when the power consumption is high daytime, in order to satisfy the cooling demand of mill, according to connection condition and the mode of refrigerator and ice storage equipment, ice storage air conditioning system can be divided into parallel and series connection two kinds, establish ties and parallelly connected two kinds of modes according to ice storage air conditioning system, further perfect comprehensive energy system economic optimization dispatch model, make the dispatch model after optimizing more agree with the actual demand of engineering, specifically include following two modes:
(i) a parallel ice storage air conditioner based on a dual-working condition refrigerator comprises: the refrigerator and the ice storage tank of the parallel ice storage air conditioning system are in parallel connection in the system, wherein the refrigerator and the ice storage tank can jointly supply cold and can also independently supply cold load, and the refrigerator can simultaneously make and supply cold;
Figure FDA00032529790400000314
Figure FDA0003252979040000041
Figure FDA0003252979040000042
Figure FDA0003252979040000043
Figure FDA0003252979040000044
Figure FDA0003252979040000045
in the formula:
Figure FDA0003252979040000046
and
Figure FDA0003252979040000047
respectively representing the refrigerating power of the ith refrigerating machine and the ice storage tank in the time t;
Figure FDA0003252979040000048
and
Figure FDA0003252979040000049
respectively representing the maximum refrigerating power of the ith refrigerating machine and the ice storage tank;
Figure FDA00032529790400000410
and
Figure FDA00032529790400000411
respectively representing the electric power of the ith refrigerating machine and the ice storage tank in the time t;
Figure FDA00032529790400000412
and
Figure FDA00032529790400000413
respectively representing the maximum electric power of the ith refrigerating machine and the ice storage tank;
Figure FDA00032529790400000414
and
Figure FDA00032529790400000415
respectively representing the total electric power, the maximum electric power and the refrigerating power of the ith ice storage air-conditioning system in a time period tth; t ismeltIndicating being in the ice-melting period, TrefIndicating the ice accumulation period, TmeltAnd TrefThe formula shows that the ice storage and the ice melting operation of the ice storage tank can not be carried out at the same time;
Figure FDA00032529790400000416
representing the refrigeration energy efficiency ratio of the refrigerator;
Figure FDA00032529790400000417
and
Figure FDA00032529790400000418
respectively representing the ice making energy efficiency ratio and the ice melting efficiency of the ice storage tank;
Figure FDA00032529790400000419
and
Figure FDA00032529790400000420
respectively representing the ice storage capacity of the ith ice storage tank time period t +1 and the time period t;
Figure FDA00032529790400000421
is the self-loss coefficient of the ice storage tank; in the working mode of (i), the specific control variable of the ice storage air conditioning system is the flow of the circulating glycol passing through the ice storage tank and the refrigerating machine, and the cooling capacity of the ice storage tank and the refrigerating machine and the flow of the circulating glycol have the following relationship:
Figure FDA0003252979040000051
Figure FDA0003252979040000052
in the formula:
Figure FDA0003252979040000053
and
Figure FDA0003252979040000054
respectively representing the flow rates of the circulating glycol passing through the ice storage tank and the refrigerating machine in a time period t; cgly、ρglyAnd Δ TglyThe specific heat capacity, the liquid density and the return water supply temperature difference of the ethylene glycol solution are respectively;
Figure FDA0003252979040000055
the refrigeration efficiency of the refrigerator is improved;
(ii) a series type ice storage air conditioner based on a dual-working condition refrigerator comprises: the refrigerating machine and the ice storage tank of the tandem type ice storage air conditioning system are in tandem position in the system, the cold quantity distribution of the refrigerating machine and the ice storage tank meets the preset proportional relation, and the proportional relation met by the cold quantity distribution of the refrigerating machine and the ice storage tank is mainly embodied in the following two stages:
(I) in the stage of ice storage, the refrigerating machine produces cold and stores the cold in the ice storage tank, at the moment, the ice storage tank does not participate in the refrigeration operation, the refrigerating machine participates in the refrigeration operation, and the ice making power of the ice storage tank
Figure FDA0003252979040000056
And cooling power of refrigerator
Figure FDA0003252979040000057
The relationship is as follows:
Figure FDA0003252979040000058
in the formula:
Figure FDA0003252979040000059
and
Figure FDA00032529790400000510
the temperature difference of the ice storage tank and the glycol at the inlet and the outlet of the refrigerator is t;
(II) in the cold supply stage, the ice storage tank and the refrigerator must be supplied with cold simultaneously, and the cold quantity distribution of the ice storage tank and the refrigerator meets the preset proportional relation:
Figure FDA00032529790400000511
Figure FDA00032529790400000512
in the formula: epsilons.iThe distribution coefficient of the cooling capacity of the ith ice storage air conditioning system.
2. The micro-grid based economic optimization scheduling method for the integrated energy system, according to claim 1, wherein in the step (3),
(A) the operation and maintenance cost is as follows:
Figure FDA0003252979040000061
in the formula: xiOM.iThe operating maintenance cost per unit output power of the equipment i;
Figure FDA0003252979040000062
represents the output power of the ith device during time period t;
(B) the electricity purchasing cost is as follows:
Figure FDA0003252979040000063
in the formula:
Figure FDA0003252979040000064
and
Figure FDA0003252979040000065
the electricity purchase price and the electricity purchase power of the time period t are respectively;
Figure FDA0003252979040000066
and
Figure FDA00032529790400000614
the price and power of selling electricity in time t;
(C) fuel cost:
Figure FDA0003252979040000067
in the formula:
Figure FDA0003252979040000068
and
Figure FDA0003252979040000069
gas consumption rates of the ith gas turbine and the ith gas boiler, respectively, for a time period tth;
Figure FDA00032529790400000610
is the gas price;
(D) energy storage depreciation cost:
along with the deepening of the discharge depth, the number of times of charge and discharge of the stored energy of the battery is reduced, but the total amount of charge and discharge of the battery is basically unchanged, if the total amount of charge and discharge of the stored energy of the battery in the whole life cycle is constant, the depreciation cost of obtaining the accumulated discharge of the stored energy of the battery by 1kWh is as follows:
Figure FDA00032529790400000611
in the formula: cbat.repReplacement cost for stored energy, qlifetimeOutputting total quantity for the whole service life of the energy storage monomer;
the depreciation cost of stored energy is:
Figure FDA00032529790400000612
wherein:
Figure FDA00032529790400000613
the discharge power of the ith battery in the time period t is stored.
3. The microgrid-based economic optimization scheduling method for the integrated energy system, as recited in claim 2, wherein in the step (3), the cooling, heating and power balance constraints comprise an electric power balance constraint, a thermal power balance constraint and a cold power balance constraint; the electric power balance constraint comprises an alternating current bus total load constraint, an alternating current-direct current converter efficiency constraint, a direct current bus total load constraint, a tie line constraint and a power purchasing and selling state constraint, and the specific constraint conditions are as follows:
the total load constraint of the alternating current bus:
Figure FDA0003252979040000071
in the formula:
Figure FDA0003252979040000072
an AC load for a time period t;
Figure FDA0003252979040000073
electric power for the ac-dc converter;
Figure FDA0003252979040000074
total electric power for the user air conditioner;
(II) efficiency constraint of the AC-DC converter:
Figure FDA0003252979040000075
in the formula:
Figure FDA0003252979040000076
the total load of the direct current bus is time t; etaA/DThe conversion efficiency from alternating current to direct current; etaD/AThe conversion efficiency from direct current to alternating current;
and (III) total load constraint of the direct current bus:
Figure FDA0003252979040000077
in the formula:
Figure FDA0003252979040000078
a direct current load for a time period t;
and (IV) connecting line constraint and electricity purchasing and selling state constraint:
Figure FDA0003252979040000079
Figure FDA00032529790400000710
Figure FDA00032529790400000711
in the formula:
Figure FDA0003252979040000081
and
Figure FDA0003252979040000082
respectively the upper power limit of electricity purchasing and electricity selling to the power grid;
Figure FDA0003252979040000083
and
Figure FDA0003252979040000084
respectively 0-1 state variables of purchasing and selling electricity in the time period t,
Figure FDA0003252979040000085
taking 1 indicates that electricity is purchased,
Figure FDA0003252979040000086
and 1 is taken to represent electricity sales, and the condition that electricity cannot be purchased at the same time is limited.
4. The micro-grid based economic optimization scheduling method for the integrated energy system, as claimed in claim 3, wherein the thermal power balance constraint conditions are as follows:
Figure FDA0003252979040000087
Figure FDA0003252979040000088
in the formula:
Figure FDA0003252979040000089
and
Figure FDA00032529790400000810
space thermal load and hot water load of plant equipment, respectively.
5. The micro-grid based economic optimization scheduling method for the integrated energy system as claimed in claim 4, wherein the constraint conditions of the cold power balance constraint are as follows:
Figure FDA00032529790400000811
in the formula:
Figure FDA00032529790400000812
is the cooling load.
6. The microgrid-based economic optimization scheduling method for the integrated energy system, as recited in claim 1, wherein in the step (3), the constraint conditions of the equipment operation constraints are as follows:
Figure FDA00032529790400000813
Figure FDA00032529790400000814
in the formula:
Figure FDA00032529790400000815
and
Figure FDA00032529790400000816
respectively representing the input and output power of the device i in a time period t;
Figure FDA00032529790400000817
and
Figure FDA00032529790400000818
respectively representing the upper and lower limits of the output power of the device i in the time period t;
Figure FDA00032529790400000819
and
Figure FDA00032529790400000820
respectively representing the upper and lower limits of the input power of the device i during the time period t.
7. The microgrid-based economic optimization scheduling method for the comprehensive energy system is characterized in that in the step (3), the energy storage device constraint needs to meet the energy storage state constraint and the charge-discharge energy power constraint, and in order to ensure the scheduling continuity, the energy storage of the energy storage device is kept consistent before and after the scheduling period; the constraint conditions of the energy storage device constraint are as follows:
Figure FDA0003252979040000091
SL.i＝ST.i
Figure FDA0003252979040000092
Figure FDA0003252979040000093
Figure FDA0003252979040000094
wherein:
Figure FDA0003252979040000095
representing the storage capacity of the energy storage device over time period t;
Figure FDA0003252979040000096
and
Figure FDA0003252979040000097
representing the maximum and minimum storage capacities of the energy storage device, respectively; sL.iAnd ST.iThe initial capacity of the stored energy and the capacity at the end of the scheduling period;
Figure FDA0003252979040000098
and
Figure FDA0003252979040000099
representing the charging and discharging power of the energy storage device, respectively, over time period t;
Figure FDA00032529790400000910
and
Figure FDA00032529790400000911
representing the maximum charge and discharge power of the energy storage device, respectively;
Figure FDA00032529790400000912
and
Figure FDA00032529790400000913
respectively representing the energy storage device in a 0-1 state variable for charging and discharging energy during time period t,
Figure FDA00032529790400000914
taking 1 as the energy to be charged,
Figure FDA00032529790400000915
and 1 is taken to represent energy release, so that the equipment cannot be charged and released simultaneously.
CN201710738641.0A 2017-08-24 2017-08-24 Comprehensive energy system economic optimization scheduling method based on micro-grid Active CN107609684B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710738641.0A CN107609684B (en) 2017-08-24 2017-08-24 Comprehensive energy system economic optimization scheduling method based on micro-grid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710738641.0A CN107609684B (en) 2017-08-24 2017-08-24 Comprehensive energy system economic optimization scheduling method based on micro-grid

Publications (2)

Publication Number Publication Date
CN107609684A CN107609684A (en) 2018-01-19
CN107609684B true CN107609684B (en) 2021-12-03

Family

ID=61065897

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710738641.0A Active CN107609684B (en) 2017-08-24 2017-08-24 Comprehensive energy system economic optimization scheduling method based on micro-grid

Country Status (1)

Country Link
CN (1) CN107609684B (en)

Families Citing this family (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108206543B (en) * 2018-02-05 2021-06-04 东北大学 Energy router based on energy cascade utilization and operation optimization method thereof
CN109375588B (en) * 2018-02-12 2021-01-01 浙江大学 Factory multi-energy complementary optimization scheduling method considering generalized demand response
CN108487994B (en) * 2018-02-28 2019-11-05 中国科学院广州能源研究所 A kind of micro- energy net composite energy storage system
CN108830743B (en) * 2018-05-25 2021-10-15 天津大学 Optimal scheduling method of park comprehensive energy system considering various cold accumulation devices
CN110739710A (en) * 2018-07-20 2020-01-31 中国农业大学 Method and device for coordinated scheduling of multiple energy types based on optimization algorithm
CN109146182A (en) * 2018-08-24 2019-01-04 南京理工大学 The economic load dispatching method of meter and the distributed triple-generation system of a variety of energy storage
CN109345012B (en) * 2018-08-29 2021-09-21 华南理工大学 Park energy Internet operation optimization method based on comprehensive evaluation indexes
CN109376912B (en) * 2018-09-29 2021-07-02 东南大学 Operation optimization method of combined cooling heating and power system based on thermal inertia of civil building
WO2020081003A1 (en) * 2018-10-17 2020-04-23 Agency For Science, Technology And Research Cooling plant system and method of operating said system
CN111091404A (en) * 2018-10-24 2020-05-01 中国电力科学研究院有限公司 Micro-energy network energy selling price determining method and system
CN109345030B (en) * 2018-10-26 2022-02-15 南方电网科学研究院有限责任公司 Multi-microgrid comprehensive energy system thermoelectric energy flow distribution type optimization method and device
CN109829643A (en) * 2019-01-28 2019-05-31 山东大学 A kind of the new energy cooling heating and power generation system integrated evaluating method and system of multi-level simulation tool
CN109993345B (en) * 2019-01-29 2022-08-09 国网江苏省电力有限公司经济技术研究院 Garden-oriented dynamic economic dispatching method for multi-energy complementary system for island operation
CN109962476A (en) * 2019-02-01 2019-07-02 中国电力科学研究院有限公司 Source net lotus storage interaction energy management method and device in a kind of micro-capacitance sensor
CN110070216B (en) * 2019-04-11 2021-02-26 河海大学 Economic operation optimization method for industrial park comprehensive energy system
CN112000146B (en) * 2019-05-27 2022-04-19 南京南瑞继保电气有限公司 Scheduling method and system of air temperature adjusting system
CN110163443B (en) * 2019-05-27 2022-09-09 西南石油大学 Natural gas pressure regulating station micro-energy network optimization scheduling method considering electricity-gas comprehensive demand response
CN110361969B (en) * 2019-06-17 2021-01-05 清华大学 Optimized operation method of cooling, heating and power comprehensive energy system
CN110426590B (en) * 2019-07-15 2022-01-25 国电南瑞科技股份有限公司 Multi-energy information interaction device suitable for comprehensive energy system
CN110416992B (en) * 2019-07-24 2022-03-18 东北电力大学 Comprehensive energy optimization energy utilization method suitable for direct current power utilization users
CN110570010B (en) * 2019-07-31 2023-01-17 中国科学院广州能源研究所 Energy management method of distributed system containing heat storage device
CN110516868B (en) * 2019-08-21 2022-05-10 广东电网有限责任公司 Comprehensive energy system optimization operation model considering network constraints
CN110598313B (en) * 2019-09-10 2023-05-16 国网河北省电力有限公司 Comprehensive energy system optimal configuration method considering energy storage full life cycle operation and maintenance
CN111091227B (en) * 2019-11-14 2023-04-18 中国电建集团西北勘测设计研究院有限公司 Comprehensive energy system dispatching management platform
CN110766241B (en) * 2019-11-27 2022-05-03 广西电网有限责任公司 Demand response control method, apparatus, device and storage medium
CN110912124A (en) * 2019-12-05 2020-03-24 深圳供电局有限公司 Multi-energy complementary microgrid system
CN111191820B (en) * 2019-12-17 2023-05-09 国网浙江省电力有限公司 Site selection and volume fixation optimization planning method for energy storage device in comprehensive energy system
CN112257899A (en) * 2020-09-22 2021-01-22 国网河北省电力有限公司营销服务中心 CCHP system optimal scheduling method and terminal equipment
CN112398164B (en) * 2020-10-30 2022-06-28 东南大学 Micro-energy-source network group optimization operation and cost distribution method containing shared energy storage system
CN112803401B (en) * 2021-01-31 2022-12-27 国网黑龙江省电力有限公司 Regulation and control method and device of virtual distributed energy cluster and terminal equipment
CN112861335B (en) * 2021-02-01 2022-11-15 昆明理工大学 P2G and energy storage-containing low-carbon economic dispatching method for comprehensive energy system
CN113378409A (en) * 2021-07-06 2021-09-10 国网江苏省电力有限公司营销服务中心 Comprehensive energy system multi-energy complementary optimization scheduling method and system

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103246263B (en) * 2013-04-22 2015-04-15 天津大学 General optimized dispatching strategy for combined supply of cooling, heating and power microgrid system
CN103617460A (en) * 2013-12-06 2014-03-05 天津大学 Double-layer optimization planning and designing method for combined cooling, heating and power micro-grid system
CN105303247A (en) * 2015-09-16 2016-02-03 北京国电通网络技术有限公司 Garden type hot and cold energy mixed application energy network regulation method and system
CN105869075A (en) * 2016-04-19 2016-08-17 东南大学 Economic optimization scheduling method for cold, heat and electricity combined supply type miniature energy grid
CN106786793B (en) * 2016-12-14 2019-04-09 东南大学 A kind of supply of cooling, heating and electrical powers type microgrid operation method based on robust optimization
CN106786753B (en) * 2016-12-29 2019-08-06 上海博翎能源科技有限公司 The system and its adjusting method of the Regional Energy internet of multi-user
CN106709610B (en) * 2017-01-12 2020-04-21 浙江大学 Micro-grid electricity energy storage and ice storage combined optimization scheduling method

Also Published As

Publication number Publication date
CN107609684A (en) 2018-01-19

Similar Documents

Publication Publication Date Title
CN107609684B (en) Comprehensive energy system economic optimization scheduling method based on micro-grid
CN107832979B (en) Factory comprehensive energy system economic optimization scheduling method considering energy cascade utilization
CN108717594B (en) Economic optimization scheduling method for combined cooling heating and power type multi-microgrid system
He et al. A review on the influence of intelligent power consumption technologies on the utilization rate of distribution network equipment
CN106651026B (en) Multi-time scale microgrid energy management optimization scheduling method
CN109919478B (en) Comprehensive energy microgrid planning method considering comprehensive energy supply reliability
CN109523052B (en) Virtual power plant optimal scheduling method considering demand response and carbon transaction
CN104716644B (en) Renewable energy source cooling, heating and power microgrid system and control method
CN103617460A (en) Double-layer optimization planning and designing method for combined cooling, heating and power micro-grid system
CN106022503A (en) Micro-grid capacity programming method meeting coupling type electric cold and heat demand
CN107560034B (en) Optimal scheduling method of ice storage air conditioning system
CN109474025B (en) Optimized dispatching model of park level comprehensive energy system
CN110619110A (en) Coordinated operation optimization method for comprehensive energy system with heat pump
Tasdighi et al. Energy management in a smart residential building
CN102593855B (en) Method for stabilizing fluctuation of output power of renewable energy power supply in power system
CN114430175A (en) Peak clipping, valley stopping and power distribution method based on fused salt energy storage system
CN107565605A (en) A kind of shop equipment based on micro-capacitance sensor tends to the method for optimization automatically
Zhou et al. Demand response control strategy of groups of central air-conditionings for power grid energy saving
Teng et al. A novel economic analyzing method for CCHP systems based on energy cascade utilization
CN108736518B (en) Comprehensive energy supply system and method for urban complex and large public building group
Xu et al. Optimal intraday rolling operation strategy of integrated energy system with multi-storage
CN113128868B (en) Regional comprehensive energy system scheduling optimization method and device
CN115081700A (en) Comprehensive energy storage technology-based data center multi-energy collaborative optimization method and system
Pasban-Gajan et al. Optimal scheduling of renewable-based energy hubs considering time-of-use pricing scheme
Shi et al. Economic operation of industrial microgrids with multiple kinds of flexible loads

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 310052 Room 1708, Hanshi Building, 1786 Binsheng Road, Changhe Street, Binjiang District, Hangzhou City, Zhejiang Province

Applicant after: Wanke Energy Technology Co., Ltd.

Address before: 310000 Room 1708, Hanshi Building, 1786 Binsheng Road, Changhe Street, Binjiang District, Hangzhou City, Zhejiang Province

Applicant before: Zhejiang Wanke Amperex Technology Limited

GR01 Patent grant
GR01 Patent grant