CN115619145B - Cooperative control method and device for comprehensive energy system and computer equipment - Google Patents
Cooperative control method and device for comprehensive energy system and computer equipment Download PDFInfo
- Publication number
- CN115619145B CN115619145B CN202211261103.4A CN202211261103A CN115619145B CN 115619145 B CN115619145 B CN 115619145B CN 202211261103 A CN202211261103 A CN 202211261103A CN 115619145 B CN115619145 B CN 115619145B
- Authority
- CN
- China
- Prior art keywords
- representing
- energy system
- park
- power
- equipment
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000007789 gas Substances 0.000 claims description 90
- 238000005338 heat storage Methods 0.000 claims description 53
- 238000003860 storage Methods 0.000 claims description 43
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 36
- 238000004146 energy storage Methods 0.000 claims description 32
- 230000005611 electricity Effects 0.000 claims description 30
- 230000006870 function Effects 0.000 claims description 23
- 239000003345 natural gas Substances 0.000 claims description 18
- 238000002844 melting Methods 0.000 claims description 15
- 238000001816 cooling Methods 0.000 claims description 14
- 238000005485 electric heating Methods 0.000 claims description 14
- 239000002918 waste heat Substances 0.000 claims description 13
- 230000008018 melting Effects 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 10
- 238000010438 heat treatment Methods 0.000 claims description 10
- 238000012360 testing method Methods 0.000 claims description 9
- 238000005286 illumination Methods 0.000 claims description 6
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 238000005057 refrigeration Methods 0.000 claims description 5
- 238000009825 accumulation Methods 0.000 claims description 4
- 238000010248 power generation Methods 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 241001368098 Capis Species 0.000 claims description 2
- 101100112125 Staphylococcus aureus capI gene Proteins 0.000 claims description 2
- 238000007599 discharging Methods 0.000 claims 1
- 239000003795 chemical substances by application Substances 0.000 description 22
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 230000003993 interaction Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000002737 fuel gas Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000002618 waking effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06313—Resource planning in a project environment
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Marketing (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- Operations Research (AREA)
- Biodiversity & Conservation Biology (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Quality & Reliability (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Feedback Control In General (AREA)
- Engine Equipment That Uses Special Cycles (AREA)
Abstract
The invention discloses a cooperative control method and device for a comprehensive energy system and computer equipment. According to the method, a steady-state model of each energy device is established, profit maximization of a comprehensive energy system is used as a first objective function, comprehensive energy system constraint is used as a first constraint condition, a comprehensive energy system scheduling model is established, the comprehensive energy system scheduling model is solved, the gate power value of each park is obtained, the operation cost minimization of each park is used as a second objective function based on the gate power value of each park, the constraint of each park is used as a second constraint condition, the scheduling model of each park is established, the scheduling model of each park is solved, and scheduling instructions of each energy device are obtained, so that each energy device executes corresponding scheduling instructions, the overall economy of the comprehensive energy system and the economy of each park in the comprehensive energy system can be guaranteed simultaneously, and the difficulty in optimizing control of the comprehensive energy system is reduced.
Description
Technical Field
The present invention relates to the field of energy scheduling technologies, and in particular, to a cooperative control method and apparatus for a comprehensive energy system, and a computer device.
Background
In recent years, aiming at a series of problems of resource shortage, unreasonable structure, low utilization efficiency and the like faced by energy production and consumption, countries around the world have reached a high consensus in the aspect of building a new generation power system (comprehensive energy system) mainly comprising new energy on the basis of continuous intensive research and communication. The comprehensive energy system is a main bearing form of a novel power system integrating multiple energy sources such as electric power, fuel gas, heat supply/cold supply and the like, and the distributed energy sources can be automatically used by plug and play as long as the distributed energy sources meet an access standard, but impact trend fluctuation with multiple space-time distribution can be generated at the same time, so that the system operation is endangered. Under an open plug and play access environment, the method is very critical in flexibly waking up and cooperatively controlling mass distributed controllable resources and realizing the rapid stabilization of energy fluctuation so as to ensure the high-quality energy supply of the system. The comprehensive energy system is large in scale, system equipment is complicated, and the performances of the equipment are various and mutually influenced, so that the dynamic process of the whole system is complex. In the system operation, only a single control is implemented on each device, the mutual coordination among the devices cannot be realized, and the operation targets of the system, such as safety, reliability, economy and the like, are difficult to meet, so that the coordination control problem among the devices needs to be considered.
At present, a great deal of related researches are carried out on the coordinated control of the comprehensive energy system at home and abroad, but the researches are mostly based on the development of an electric power system, the high fusion characteristics of the comprehensive energy system to various energy sources such as electric power, fuel gas, heat supply/cold supply and the like are not considered, and meanwhile, related researches on the coordination control of a large amount of distributed controllable resources in the comprehensive energy system and heterogeneous hybrid groups are fresh.
Therefore, how to study a cooperative control method of an integrated energy system is a problem to be solved by those skilled in the art aiming at the mixed population problem existing in mass distributed controllable resource control.
Disclosure of Invention
The purpose of the invention is that: the cooperative control method, the cooperative control device, the computer equipment and the computer readable storage medium for the comprehensive energy system can simultaneously ensure the overall economy of the comprehensive energy system and the economy of each park in the comprehensive energy system, and reduce the difficulty of optimizing and controlling the comprehensive energy system.
In order to achieve the above object, a first aspect of the present invention provides a cooperative control method for an integrated energy system including a plurality of parks, each park including a plurality of energy devices, the cooperative control method comprising:
Establishing a steady-state model of each energy device;
according to the steady-state model, taking profit maximization of the comprehensive energy system as a first objective function, taking comprehensive energy system constraint as a first constraint condition, establishing a comprehensive energy system scheduling model, and solving the comprehensive energy system scheduling model to obtain the gate power values of all parks, wherein the comprehensive energy system constraint comprises comprehensive energy system electric power balance constraint, comprehensive energy system gas power balance constraint and comprehensive energy system peak power constraint;
according to the gate power value, minimizing the operation cost of each park as a second objective function, taking each park constraint as a second constraint condition, establishing each park scheduling model, and solving each park scheduling model to obtain scheduling instructions of each energy device, wherein the park constraint comprises park electric power balance constraint, park thermal power balance constraint, park cold power balance constraint and energy device operation constraint;
and sending the scheduling instruction to the corresponding energy equipment so that the energy equipment executes the scheduling instruction.
Preferably, the steady-state model comprises a photovoltaic device steady-state model, a fan device steady-state model, a gas turbine device steady-state model, a gas boiler device steady-state model, an electric energy storage device steady-state model, a cold storage device steady-state model, a heat storage device steady-state model, a cold supply device steady-state model and a heat supply device steady-state model.
Preferably, the expression of the steady-state model of the photovoltaic device is:
P PV ＝P STC G Ac [1+k(T c -T r )]/G STC ；
wherein P is PV Representing the actual output of the photovoltaic device, P STC G represents the test output of the photovoltaic equipment under standard experimental conditions AC Representing the actual illumination intensity, G STC Represents standard illumination intensity, k represents temperature coefficient, T c Represents the actual temperature value, T r Representing a standard temperature value;
the expression of the fan equipment steady-state model is as follows:
wherein P is WT Representing the actual output of the fan equipment, v represents the wind speed, v out Represents cut-out wind speed v in Representing cut-in wind speed, v r Represents rated wind speed, P WT (v i ) Representing the wind speed v of the fan equipment i Test force under P WT (v i+1 ) Representing the wind speed v of the fan equipment i+1 Test output;
the expression of the steady state model of the gas turbine plant is:
wherein F is MT Representing the natural gas consumption of a gas turbine plant, P MT Represents the output electric power, eta of the gas turbine plant MT Representing the power generation efficiency of the gas turbine plant;
the expression of the steady-state model of the gas boiler equipment is as follows:
wherein F is GB Represents the natural gas consumption of the gas boiler equipment, Q GB Represents the heating power, eta of the gas boiler equipment GB Representing the corresponding heating value of the input gas boiler plant;
the expression of the steady-state model of the electric energy storage device is:
0≤P ES，C ≤CapEsγ ES，C ；
0≤P ES，D ≤CapEsγE S，D ；
W ES，min ≤W ES ≤W ES ， max ；
Wherein,representing the residual electric quantity of the electric energy storage device before and after charging and discharging sigma ES Representing the self-loss rate of the electric energy storage device, P ES，C Representing the charge power, eta of an electrical energy storage device ES，C Representing the charging efficiency, P, of an electrical energy storage device ES，D Represents the discharge power, eta of an electrical energy storage device ES，D Representing the discharge efficiency of an electrical energy storage device, gamma ES，C Representing the maximum charge rate of the electrical energy storage device, gamma ES，D Representing maximum discharge rate of an electrical energy storage deviceCapes represents the rated capacity, W, of the electrical energy storage device ES，min 、W ES，max Representing a minimum allowable capacity and a maximum allowable capacity of the electrical energy storage device;
the expression of the steady-state model of the cold accumulation device is as follows:
0≤Q IS，C ≤CapIsγ IS，C ；
0≤Q IS，D ≤CapISγ IS，D ；
W IS，min ≤W IS ≤W IS，max ；
wherein,representing the residual energy, sigma, of the cold storage device before and after ice making and ice melting IS Represents the self-loss rate of the cold accumulation device, Q IS，C Ice-making power, eta, representative of cold storage device IS，C Represents ice making efficiency, Q of cold storage device IS，D Ice melting power, eta, representing cold storage device IS，D Represents ice melting efficiency, gamma of cold storage equipment IS，C Represents the maximum ice making multiplying power of cold storage equipment, gamma IS，D Representing the maximum ice-melting rate of the cold storage device, capls representing the rated capacity of the cold storage device, W IS，min 、W IS，max Representing a minimum allowable capacity and a maximum allowable capacity of the cold storage device;
the expression of the steady-state model of the heat storage device is:
0≤Q TS，C ≤CapTsγ TS，C ；
0≤Q TS，D ≤CapTsγ TS，D ；
W TS,min ≤W TS ≤W TS,max ；
Wherein,representing the residual energy, sigma, of the heat storage device before and after heat storage and heat release TS Represents the self-loss rate of the heat storage device, Q TS,C Represents the heat storage power, eta of the heat storage device TS,C Representing the heat storage efficiency of the heat storage device, Q TS,γ Represents the exothermic power, eta of the heat storage device TS,D Representing the exothermic efficiency of the thermal storage device, gamma TS,C Representing the maximum thermal storage rate of the thermal storage device, caps gamma TS,D Representing the maximum heat release rate of the heat storage device, capts representing the rated capacity of the heat storage device, W TS,min 、W TS,max Representing the minimum allowable capacity and the maximum allowable capacity of the heat storage device;
the expression of the steady-state model of the cold supply equipment is as follows:
wherein,representing the cold power output by the cold supply device, COP EC σ IS Representing the energy efficiency ratio of the cooling equipment, P EC Represents the power consumption of the cooling device;
the expression of the steady-state model of the heating equipment is as follows:
wherein,representing the thermal power output by the electric heating device, COP HP Representing the energy efficiency of an electric heating plantRatio, P HP Represents the power consumption of the electric heating device, +.>Represents the heat power input by the waste heat boiler, eta WH Representing the efficiency of the waste heat boiler,representing the thermal power output by the waste heat boiler.
Preferably, the expression of the first objective function is:
max F 1 ＝E E -c grid -C gas ；
wherein F is 1 Represents profit and income of comprehensive energy system, E E Represents the electricity selling income of the comprehensive energy system, C grid Representing the electricity purchasing expense of the comprehensive energy system, C gas Represents the gas purchase cost of the comprehensive energy system, P grid Representing the power purchased by the comprehensive energy system from the external power grid through the power transmission line, c s Representing electricity price of electricity selling of comprehensive energy system to large power grid, c b Representing electricity price of electricity purchase of comprehensive energy system to large power grid, C gas，i Representing the gas purchase cost of the ith park in the integrated energy system.
Preferably, the expression of the integrated energy system electric power balance constraint is:
wherein P is AE,i Representing the gate power value of the ith park in the comprehensive energy system;
the expression of the gas power balance constraint of the comprehensive energy system is as follows:
wherein F is gas Represents the total amount of natural gas input by the comprehensive energy system, F MT，i Representing natural gas consumption of gas turbine plant at ith park in integrated energy system, F GB，i Representing the natural gas consumption of the gas boiler equipment of the ith park in the comprehensive energy system;
the expression of the peak power constraint of the comprehensive energy system is as follows:
P grid ≤P line，max ；
wherein P is line,max Representing the maximum allowable value of the tie-line flow between the integrated energy system and the upper grid.
Preferably, the expression of the second objective function is:
min f i ＝C E,i +C H,i +C C,i +C gas,i ；
wherein f i Representing the operation cost of the ith park in the comprehensive energy system, C E,i Representing electricity purchase cost of ith park in comprehensive energy system, C H,i Representing the heat purchase cost of the ith park in the comprehensive energy system, C C,i Representing the cold purchase cost of the ith park in the integrated energy system,represents the time-of-use electricity price at time t, +.>Representing electricity purchasing power of ith park at t moment in comprehensive energy system, c H Represents the heat supply price, Q H,i Representing the heat purchasing power of the ith park in the comprehensive energy system, c C Represents the price of cooling, Q C,i Representing the cold purchase power of the ith park in the integrated energy system, < >>Represents the price of natural gas at time t, +.>Representing the power generated by the gas turbine plant at time t, η of the ith park in the integrated energy system MT Representing the efficiency of a gas turbine plant, +.>Heat production power eta of gas boiler equipment representing ith park in comprehensive energy system at t moment GB Representing the efficiency of the gas boiler plant.
Preferably, the park electric power balance constraint is expressed as:
P AE,i +P MT,i +P PV,i +P WT,i +P ES,D,i ＝P Eload,i +P HP,i +P ES,C,i +P EC,i ；
wherein,representing the power generated by the gas turbine plant of the ith park in the integrated energy system, P PV,i Representing the power generated by photovoltaic devices in the ith park in the integrated energy system, P WT,i Representing the power generated by fan equipment of the ith park in the integrated energy system, P Eload,i Representing the electric load power of the ith park in the integrated energy system, P HP，i Representing the power consumption of the electric heating equipment of the ith park in the comprehensive energy system, P ES,C,i Representing the charging power of the electric energy storage device of the ith park in the integrated energy system, P ES,D,i Representing the discharge power, P, of an electric energy storage device of an ith park in an integrated energy system EC,i Representing the power consumption of the electric refrigeration equipment of the ith park in the comprehensive energy system;
the park thermal power balance constraint expression is:
wherein Q is GB,i Representing the heat power output by gas boiler equipment of ith park in comprehensive energy system, Q TS,D,i Representing the thermal power output by the heat storage equipment of the ith park in the integrated energy system, Q TS,C,i Representing the thermal power input by the heat storage equipment of the ith park in the integrated energy system,representing the thermal power output by the electric heating equipment of the ith park in the integrated energy system,/>Representing the heat power input by the waste heat boiler equipment of the ith park in the comprehensive energy system, Q HL,i Representing the heat load power of an ith park in the comprehensive energy system;
the park cold power balance constraint expression is:
wherein, Represents the cold power output by the electric refrigeration equipment of the ith park in the comprehensive energy system, Q IS,D,i Ice melting power, Q, representing heat storage equipment of ith park in integrated energy system IS,C,i An ice making power input by a heat storage device representing an ith park in the integrated energy system, Q CL Representing the cooling load power of an ith park in the comprehensive energy system;
the expression of the energy device operation constraint is:
wherein P is min,k Represents the minimum operation allowable electric power, P, of the kth energy device in the integrated energy system max,k Represents the maximum operation allowable electric power of the kth energy device in the integrated energy system, Q min,k Represents the minimum operation allowable thermal power, Q, of the kth energy device in the integrated energy system max,k Represents the maximum operation allowable thermal power, W, of the kth energy device in the integrated energy system min,k Representing the minimum allowable operation capacity, W, of the kth energy device in the integrated energy system max,k Representing the maximum allowable operating capacity of the kth energy device in the integrated energy system.
A second aspect of the present invention provides a cooperative control apparatus for an integrated energy system including a plurality of parks, each park including a plurality of energy devices, the cooperative control apparatus comprising:
the steady-state model building module is used for building a steady-state model of each energy device;
The gateway power value solving module is used for establishing a comprehensive energy system scheduling model by taking profit maximization of the comprehensive energy system as a first objective function and comprehensive energy system constraint as a first constraint condition according to the steady-state model, and solving the comprehensive energy system scheduling model to obtain the gateway power value of each park, wherein the comprehensive energy system constraint comprises comprehensive energy system electric power balance constraint, comprehensive energy system gas power balance constraint and comprehensive energy system peak power constraint;
the scheduling instruction solving module is used for establishing each park scheduling model by taking the lowest running cost of each park as a second objective function and taking each park constraint as a second constraint condition according to the gateway power value, and solving each park scheduling model to obtain the scheduling instruction of each energy device, wherein the park constraint comprises park electric power balance constraint, park thermal power balance constraint, park cold power balance constraint and energy device running constraint;
and the scheduling instruction sending module is used for sending the scheduling instruction to the corresponding energy equipment so as to enable the energy equipment to execute the scheduling instruction.
A third aspect of the invention provides a computer device comprising a memory storing a computer program and a processor implementing the steps of the above-described cooperative control method for an integrated energy system when the computer program is executed by the processor.
A fourth aspect of the present invention provides a computer readable storage medium having a computer program stored thereon, which when executed by a processor, implements the steps of the above-described cooperative control method for an integrated energy system.
The invention has at least the following beneficial effects:
according to the method, the steady-state model of each energy device is built, profit maximization of the comprehensive energy system is used as a first objective function, comprehensive energy system constraint is used as a first constraint condition, the comprehensive energy system scheduling model is built, the comprehensive energy system scheduling model is solved, the gate power value of each park can be obtained, the operation cost minimization of each park is used as a second objective function based on the gate power value of each park, the constraint of each park is used as a second constraint condition, the park scheduling model is built, the park scheduling model is solved, scheduling instructions of each energy device can be obtained, further each energy device executes corresponding scheduling instructions, the overall economy of the comprehensive energy system and the economy of each park in the comprehensive energy system can be guaranteed simultaneously, and the difficulty of optimizing and controlling the comprehensive energy system is reduced.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required to be used in the embodiments or the description of the prior art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments described in the present invention, and other drawings may be obtained according to the drawings without inventive effort to those skilled in the art.
FIG. 1 is a schematic flow chart of a cooperative control method for an integrated energy system according to an embodiment of the present invention;
fig. 2 is a schematic structural diagram of a cooperative control apparatus for an integrated energy system according to an embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments.
It will be understood that when an element is referred to as being "fixed to" another element, it can be directly on the other element or intervening elements may also be present. When an element is referred to as being "connected" to another element, it can be directly connected to the other element or intervening elements may also be present. The terms "vertical," "horizontal," "left," "right," and the like are used herein for illustrative purposes only and are not meant to be the only embodiment.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. The term "and/or" as used herein includes any and all combinations of one or more of the associated listed items.
The embodiment of the invention mainly provides a cooperative control method for a building integrated energy system, and in order to more clearly illustrate the scheme, the structure of the building integrated energy system is firstly illustrated.
The embodiment of the invention provides a cooperative control method for a comprehensive energy system. Referring to fig. 1, fig. 1 is a flow chart of a cooperative control method for an integrated energy system according to an embodiment of the invention. Wherein the integrated energy system includes a plurality of campuses, each park including a plurality of energy devices, the cooperative control method may include:
s110, establishing a steady-state model of each energy device.
In the embodiment of the invention, the comprehensive energy system is composed of a plurality of parks, each park comprises a plurality of energy devices such as cold energy storage devices, heat energy storage devices, electricity energy storage devices, gas turbine devices, load adjustment devices and the like, and the flexible controllable units comprise photovoltaic devices, electric energy storage devices, wind power devices, gas turbine devices, load adjustment devices and the like, wherein the cold energy storage devices and the heat energy storage devices are respectively. Firstly, establishing steady-state models of different energy devices (mainly comprising flexible controllable units) of each park in the comprehensive energy system.
S120, according to the steady-state model, taking profit maximization of the comprehensive energy system as a first objective function, taking comprehensive energy system constraint as a first constraint condition, establishing a comprehensive energy system scheduling model, and solving the comprehensive energy system scheduling model to obtain the gate power value of each park, wherein the comprehensive energy system constraint comprises comprehensive energy system electric power balance constraint, comprehensive energy system gas power balance constraint and comprehensive energy system peak power constraint.
In the embodiment of the invention, the constraint conditions of the comprehensive energy system can be established by analyzing the steady-state models of different energy devices in each park in the comprehensive energy system and combining the distribution conditions of the cold, hot, electric and gas energy devices and the flow relation of energy among the devices in each park, wherein the constraint conditions comprise the electric power balance constraint of the comprehensive energy system, the gas power balance constraint of the comprehensive energy system and the peak power constraint of the comprehensive energy system, and the profit maximization of the comprehensive energy system is used as a first objective function to establish the dispatching model of the comprehensive energy system. And then solving the comprehensive energy system scheduling model, and determining the corresponding gateway power range of each park in the comprehensive energy system.
And S130, according to the gateway power value, taking the operation cost minimization of each park as a second objective function, taking the constraint of each park as a second constraint condition, establishing each park scheduling model, and solving each park scheduling model to obtain the scheduling instruction of each energy device, wherein the park constraint comprises park electric power balance constraint, park thermal power balance constraint, park cold power balance constraint and energy device operation constraint.
In the embodiment of the invention, the steady-state models of different energy devices of each park in the comprehensive energy system are analyzed, the distribution condition of the cold, hot, electric and gas energy devices and the flow relation of energy among the devices of each park are combined, the constraint condition of each park can be established, the gate power value of each park is considered, the operation cost minimization of each park is used as a second objective function, and the scheduling model of each park can be established. And then solving the scheduling model of each park, and determining the scheduling instruction of each energy device in each park.
And S140, sending the scheduling instruction to the corresponding energy equipment so that the energy equipment executes the scheduling instruction.
In the embodiment of the invention, after the scheduling instruction of each energy device in each park is obtained, each energy device executes the corresponding scheduling instruction, thereby realizing the cooperative control of the comprehensive energy system.
As can be seen from the foregoing, in the cooperative control method for a comprehensive energy system provided by the embodiment of the present invention, by establishing a steady-state model of each energy device, then taking profit maximization of the comprehensive energy system as a first objective function, taking comprehensive energy system constraint as a first constraint condition, establishing a comprehensive energy system scheduling model, and solving the comprehensive energy system scheduling model, a gate power value of each park can be obtained, then based on the gate power value of each park, with operation cost minimization of each park as a second objective function, taking constraint of each park as a second constraint condition, establishing a scheduling model of each park, and solving the scheduling model of each park, scheduling instructions of each energy device can be obtained, so that each energy device executes corresponding scheduling instructions, and thus, the overall economy of the comprehensive energy system and the economy of each park in the comprehensive energy system can be ensured, and difficulty in optimizing and controlling the comprehensive energy system can be reduced.
Specifically, in the above embodiment, the steady-state model includes a photovoltaic device steady-state model, a fan device steady-state model, a gas turbine device steady-state model, a gas boiler device steady-state model, an electric energy storage device steady-state model, a cold storage device steady-state model, a heat storage device steady-state model, a cold supply device steady-state model, and a heat supply device steady-state model.
Further, in the above embodiment, the expression of the steady-state model of the photovoltaic device is:
P PV ＝P STC G AC [1+k(T c -T r )]/G STC ；
wherein P is PV Representing the actual output of the photovoltaic device, P STC G represents the test output of the photovoltaic equipment under standard experimental conditions AC Representing the actual illumination intensity, G STC Represents standard illumination intensity, k represents temperature coefficient, T c Represents the actual temperature value, T r Representing a standard temperature value;
the expression of the fan equipment steady-state model is as follows:
wherein P is WT Representing the actual output of the fan equipment, v represents the wind speed, v out Represents cut-out wind speed v in Representing cut-in wind speed, v r Represents rated wind speed, P WT (v i ) Representing the wind speed v of the fan equipment i Test force under P WT (v i+1 ) Representing the wind speed v of the fan equipment i+1 Test output;
the expression of the steady state model of the gas turbine plant is:
wherein F is MT Representing the natural gas consumption of a gas turbine plant, P MT Represents the output electric power, eta of the gas turbine plant MT Representing the power generation efficiency of the gas turbine plant;
the expression of the steady-state model of the gas boiler equipment is as follows:
wherein F is GB Represents the natural gas consumption of the gas boiler equipment, Q GB Represents the heating power, eta of the gas boiler equipment GB Representing the corresponding heating value of the input gas boiler plant;
the expression of the steady-state model of the electric energy storage device is:
0≤P ES，C ≤CapEsγ ES,C ；
0≤P ES，D ≤CapEsγ ES，D ；
W ES，min ≤W ES ≤W ES，max ；
Wherein,representing the residual electric quantity of the electric energy storage device before and after charging and discharging sigma ES Representing the self-loss rate of the electric energy storage device, P ES，C Representing the charge power, eta of an electrical energy storage device ES，C Representing the charging efficiency, P, of an electrical energy storage device ES，D Represents the discharge power, eta of an electrical energy storage device ES，D Representing the discharge efficiency of an electrical energy storage device, gamma ES，C Representing the maximum charge rate of the electrical energy storage device, gamma ES，D Representing the maximum discharge rate of the electrical energy storage device, capes represents the rated capacity of the electrical energy storage device, W ES，min 、W ES，max Representing a minimum allowable capacity and a maximum allowable capacity of the electrical energy storage device;
the expression of the steady-state model of the cold accumulation device is as follows:
0≤Q IS，C ≤CapIsγ IS，C ；
0≤Q IS，D ≤CapIsγ IS，D ；
W IS，min ≤W IS ≤W IS，max ；
wherein,representing the residual energy, sigma, of the cold storage device before and after ice making and ice melting IS Represents the self-loss rate of the cold accumulation device, Q IS，C Ice-making power, eta, representative of cold storage device IS，C Represents ice making efficiency, Q of cold storage device IS，D Ice melting power, eta, representing cold storage device IS，D Represents ice melting efficiency, gamma of cold storage equipment IS，C Represents the maximum ice making multiplying power of cold storage equipment, gamma IS，D Representing the maximum ice-melting rate of the cold storage device, capIs representing the rated capacity of the cold storage device, W IS，min 、W IS，max Representing a minimum allowable capacity and a maximum allowable capacity of the cold storage device;
the expression of the steady-state model of the heat storage device is:
0≤Q TS，C ≤CapTsγ TS，C ；
0≤Q TS，D ≤CapTsγ TS，D ；
W TS，min ≤W TS ≤W TS，max ；
Wherein,representing the residual energy, sigma, of the heat storage device before and after heat storage and heat release TS Represents the self-loss rate of the heat storage device, Q TS,C Represents the heat storage power, eta of the heat storage device TS,C Representing the heat storage efficiency of the heat storage device, Q TS,D Represents the exothermic power, eta of the heat storage device TS,D Representing the exothermic efficiency of the thermal storage device, gamma TS,C Representing the maximum thermal storage rate of the thermal storage device, caps gamma TS,D Representing the maximum heat release rate of the heat storage device, capts representing the rated capacity of the heat storage device, W TS,min 、W TS,max Representing the minimum allowable capacity and the maximum allowable capacity of the heat storage device;
the expression of the steady-state model of the cold supply equipment is as follows:
wherein,representing the cold power output by the cold supply device, COP EC σ IS Representing the energy efficiency ratio of the cooling equipment, P EC Represents the power consumption of the cooling device;
the expression of the steady-state model of the heating equipment is as follows:
wherein,representing the thermal power output by the electric heating device, COP HP Representing the energy efficiency ratio of the electric heating equipment, P HP Represents the power consumption of the electric heating device, +.>Represents the heat power input by the waste heat boiler, eta WH Representing the efficiency of the waste heat boiler,representing the thermal power output by the waste heat boiler.
In the embodiment of the invention, the energy equipment in the comprehensive energy system is complex and various, and can be divided into energy generating equipment, energy storage equipment and energy conversion equipment according to the input-output form of energy. The energy production equipment comprises photovoltaic equipment, fan equipment, gas turbine equipment and gas boiler equipment; the energy storage device comprises an electric energy storage device, a hot energy storage device and a cold energy storage device; the energy conversion equipment comprises a cooling equipment and a heating equipment. It should be noted that, the typical cold energy storage device is composed of a double-station host, a base-load host, an ice storage tank and a corresponding water circulation system, and can work in two modes of ice making and refrigeration; the heat supply equipment comprises electric heating equipment and waste heat boiler equipment.
Further, in the above embodiment, the expression of the first objective function is:
max F 1 ＝E E -C grid -C gas ；
wherein F is 1 Represents profit and income of comprehensive energy system, E E Represents the electricity selling income of the comprehensive energy system, C grid Representing the electricity purchasing expense of the comprehensive energy system, C gas Represents the gas purchase cost of the comprehensive energy system, P grid Representing the power purchased by the comprehensive energy system from the external power grid through the power transmission line, c s Representing electricity price of electricity selling of comprehensive energy system to large power grid, c b Representing electricity price of electricity purchase of comprehensive energy system to large power grid, C gas，i Representing the gas purchase cost of the ith park in the integrated energy system.
Further, in the above embodiment, the expression of the electric power balance constraint of the integrated energy system is:
wherein P is AE,i Representing the gate power value of the ith park in the comprehensive energy system;
the expression of the gas power balance constraint of the comprehensive energy system is as follows:
wherein F is gas Represents the total amount of natural gas input by the comprehensive energy system, F MT，i Representing natural gas consumption of gas turbine plant at ith park in integrated energy system, F GB，i Representing the natural gas consumption of the gas boiler equipment of the ith park in the comprehensive energy system;
the expression of the peak power constraint of the comprehensive energy system is as follows:
P gr，id ≤P line，max ；
Wherein P is line,max Representing the maximum allowable value of the tie-line flow between the integrated energy system and the upper grid.
In the embodiment of the invention, in order to ensure friendly interaction between the comprehensive energy system and the upper power grid, the peak power constraint of the gateway of the comprehensive energy system and the upper power grid tie line is also required.
Further, in the above embodiment, the expression of the second objective function is:
min f i ＝C E,i +C H,i +C C,i +C gas,i ；
wherein f i Representing the operation cost of the ith park in the comprehensive energy system, C E,i Representing electricity purchase cost of ith park in comprehensive energy system, C H,i Representing the heat purchase cost of the ith park in the comprehensive energy system, C C,i Representing the cold purchase cost of the ith park in the integrated energy system,represents the time-of-use electricity price at time t, +.>Representing electricity purchasing power of ith park at t moment in comprehensive energy system, c H Represents the heat supply price, Q H,i Representing the heat purchasing power of the ith park in the comprehensive energy system, c C Represents the price of cooling, Q C,i Representing the cold purchase power of the ith park in the integrated energy system, < >>Represents the price of natural gas at time t, +.>Representing the power generated by the gas turbine plant at time t, η of the ith park in the integrated energy system MT Representing the efficiency of a gas turbine plant, +. >Heat production power eta of gas boiler equipment representing ith park in comprehensive energy system at t moment GB Representing the efficiency of the gas boiler plant.
Further, in the above embodiment, the expression of the campus electric power balance constraint is:
P AE,i +P MT,i +P PV,i +P WT,i +P ES,D,i ＝P Eload,i +P HP,i +P ES,C,i +P EC,i ；
wherein,representing the power generated by the gas turbine plant of the ith park in the integrated energy system, P PV,i Representing the power generated by photovoltaic devices in the ith park in the integrated energy system, P WT,i Representing the power generated by fan equipment of the ith park in the integrated energy system, P Eload,i Representing the electric load power of the ith park in the integrated energy system, P HP，i Representing the power consumption of the electric heating equipment of the ith park in the comprehensive energy system, P ES,C,i Representing the charging power of the electric energy storage device of the ith park in the integrated energy system, P ES,D,i Representing the discharge power, P, of an electric energy storage device of an ith park in an integrated energy system EC,i Representing comprehensive energyThe power consumption of the electric refrigerating equipment of the ith park in the system;
the park thermal power balance constraint expression is:
wherein Q is GB,i Representing the heat power output by gas boiler equipment of ith park in comprehensive energy system, Q TS,D,i Representing the thermal power output by the heat storage equipment of the ith park in the integrated energy system, Q TS,C,i Representing the thermal power input by the heat storage equipment of the ith park in the integrated energy system,representing the thermal power output by the electric heating equipment of the ith park in the integrated energy system,/>Representing the heat power input by the waste heat boiler equipment of the ith park in the comprehensive energy system, Q HL,i Representing the heat load power of an ith park in the comprehensive energy system;
the park cold power balance constraint expression is:
wherein,represents the cold power output by the electric refrigeration equipment of the ith park in the comprehensive energy system, Q IS,D,i Ice melting power, Q, representing heat storage equipment of ith park in integrated energy system IS,C,i An ice making power input by a heat storage device representing an ith park in the integrated energy system, Q CL Representing the cooling load power of an ith park in the comprehensive energy system;
the expression of the energy device operation constraint is:
wherein P is min,k Represents the minimum operation allowable electric power, P, of the kth energy device in the integrated energy system max,k Represents the maximum operation allowable electric power of the kth energy device in the integrated energy system, Q min,k Represents the minimum operation allowable thermal power, Q, of the kth energy device in the integrated energy system max,k Represents the maximum operation allowable thermal power, W, of the kth energy device in the integrated energy system min,k Representing the minimum allowable operation capacity, W, of the kth energy device in the integrated energy system max,k Representing the maximum allowable operating capacity of the kth energy device in the integrated energy system.
In the embodiment of the invention, each energy device in each park in the integrated energy system should also meet each operation constraint.
In the specific implementation, the multi-agent concept can be introduced, and a multi-level agent cooperative control framework of the comprehensive energy system with a three-layer structure of a main station layer, an agent layer and a terminal layer can be constructed. Specifically, the master station layer is mainly responsible for comprehensively considering the overall operation safety and economy of the comprehensive energy system, performing global optimization on the system, receiving information provided by each park of the proxy layer, making decisions on interactions among the parks according to the whole-network operation target, and monitoring the operation condition of the whole network. The agent layer serves as a middle layer of the whole system, the energy management system of each park serves as an agent in the agent layer, each agent performs information interaction with the master station layer and the agents of the adjacent parks, operation scheduling, control and state monitoring of each distributed device in each park are achieved according to signals such as transaction electricity price and electric quantity demand, optimal operation output of each distributed device in each park and interaction power with a power distribution network and a power grid of the adjacent park are determined according to operation targets of the agent layer, and demand information is fed back to the master station layer and the agents of the adjacent parks. The terminal layer is used as the bottommost layer of system operation and mainly comprises various energy devices such as cold, heat, electricity, gas and the like, wherein the flexible controllable unit comprises photovoltaic, energy storage, wind power, a gas turbine, a participatable load adjustment device and the like, the cold and heat devices are respectively cold energy storage devices and heat energy devices, the cold and heat energy devices have the capabilities of optimal power generation/demand control, information storage, communication with an upper agent layer and the like, and receive operation control instructions issued by the upper agent layer, so that the control and acquisition feedback of the operation state of the device are realized.
It can be understood that the master station layer serves as an upper control, and transmits the calculated gateway power value to the proxy layer model serving as a middle control through the communication equipment, and serves as a constraint condition of the proxy layer model. The agent layer combines the constraint of the controllable units of each park, minimizes the purchase energy cost as an objective function, obtains a calculation model of the layer, and transmits the calculated execution result to the terminal-controllable unit of each park through communication. The terminal layer executes the command according to the control requirement, thereby forming the multi-layer cooperative control model of the comprehensive energy system based on multiple intelligent agents.
In specific implementation, the comprehensive energy system scheduling model and each park scheduling model can be solved based on the consistency algorithm, and the difficulty of optimizing control can be reduced by combining distributed control with the consistency algorithm. For example, the solution may be based on a consistency algorithm for information feedback:
wherein N is i(t) Representing the number of neighbor agents of agent i and being equal to its length, i.e., N i(t) ＝d ii ，x i (t) represents agents in each park of the agent layer, k 1i Feedback gain, k, representing agent i 2i Representing the difference gain, u, between agent i and adjacent agent j i (t) represents the optimization result after the agent i performs the algorithm.
It will be appreciated that when k 1i And k 2i As the value of (c) increases, the closed loop pole generated by the system becomes farther from the imaginary axis and the convergence speed of the variable value becomes faster.
The embodiment of the invention also provides a cooperative control device for the comprehensive energy system. Referring to fig. 2, fig. 2 is a schematic structural diagram of a cooperative control apparatus for an integrated energy system according to an embodiment of the present invention. Wherein the integrated energy system comprises a plurality of parks, each park comprising a plurality of energy devices, the optimized scheduling apparatus may comprise:
a steady-state model building module 100, configured to build a steady-state model of each energy device;
the gate power value solving module 200 is configured to establish a comprehensive energy system scheduling model by taking profit maximization of the comprehensive energy system as a first objective function and comprehensive energy system constraint as a first constraint condition according to the steady-state model, and solve the comprehensive energy system scheduling model to obtain the gate power value of each park, where the comprehensive energy system constraint includes comprehensive energy system electric power balance constraint, comprehensive energy system gas power balance constraint and comprehensive energy system peak power constraint;
The scheduling instruction solving module 300 is configured to establish each park scheduling model with the operation cost minimization of each park as a second objective function and each park constraint as a second constraint condition according to the gateway power value, and solve each park scheduling model to obtain the scheduling instruction of each energy device, where the park constraint includes a park electric power balance constraint, a park thermal power balance constraint, a park cold power balance constraint and an energy device operation constraint;
the scheduling instruction sending module 400 is configured to send the scheduling instruction to the corresponding energy device, so that the energy device executes the scheduling instruction.
It should be noted that, in the apparatus provided in the above embodiment, when performing the related operation, only the division of each program module is used as an example, and in practical application, the processing allocation may be performed by different program modules according to needs, that is, the internal structure of the terminal is divided into different program modules to complete all or part of the processing described above. In addition, the apparatus provided in the foregoing embodiments belongs to the same concept as the method embodiments in the foregoing embodiments, and specific implementation processes of the apparatus are detailed in the method embodiments, which are not repeated herein.
Based on the hardware implementation of the program modules, and in order to implement the method according to the embodiment of the present invention, the embodiment of the present invention further provides a computer device, including a memory and a processor, where the memory stores a computer program, and the processor implements the steps of the above-mentioned cooperative control method for the integrated energy system when executing the computer program.
In one embodiment, the present invention also provides a computer readable storage medium having a computer program stored thereon, which when executed by a processor, implements the steps of the above-described cooperative control method for an integrated energy system.
It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, and that the above embodiments and descriptions are merely illustrative of the principles of the present invention, and various changes and modifications may be made without departing from the spirit and scope of the invention, which is defined in the appended claims. The scope of the invention is defined by the appended claims and equivalents thereof.
Claims (3)
1. A cooperative control method for an integrated energy system, the integrated energy system comprising a plurality of parks, each of the parks comprising a plurality of energy devices, the cooperative control method comprising:
Establishing a steady-state model of each energy device;
according to the steady-state model, taking profit maximization of a comprehensive energy system as a first objective function, taking comprehensive energy system constraint as a first constraint condition, establishing a comprehensive energy system scheduling model, and solving the comprehensive energy system scheduling model to obtain the gate power value of each park, wherein the comprehensive energy system constraint comprises comprehensive energy system electric power balance constraint, comprehensive energy system gas power balance constraint and comprehensive energy system peak power constraint;
according to the gate power value, taking the minimum running cost of each park as a second objective function, taking each park constraint as a second constraint condition, establishing each park scheduling model, and solving each park scheduling model to obtain a scheduling instruction of each energy device, wherein the park constraint comprises park electric power balance constraint, park thermal power balance constraint, park cold power balance constraint and energy device running constraint;
the scheduling instruction is sent to the corresponding energy equipment, so that the energy equipment executes the scheduling instruction;
The steady-state model comprises a photovoltaic equipment steady-state model, a fan equipment steady-state model, a gas turbine equipment steady-state model, a gas boiler equipment steady-state model, an electric energy storage equipment steady-state model, a cold storage equipment steady-state model, a heat storage equipment steady-state model, a cold supply equipment steady-state model and a heat supply equipment steady-state model;
the expression of the steady-state model of the photovoltaic equipment is as follows:
P PV ＝P STC G AC [1+k(T c -T r )]/G STC ；
wherein P is PV Representing the actual output of the photovoltaic device, P STC G represents the test output of the photovoltaic equipment under standard experimental conditions AC Representing the actual illumination intensity, G STC Represents standard illumination intensity, k represents temperature coefficient, T c Represents the actual temperature value, T r Representing a standard temperature value;
the expression of the fan equipment steady-state model is as follows:
wherein P is WT Representing the actual output of the fan equipment, v represents the wind speed, v out Represents cut-out wind speed v in Representing cut-in wind speed, v r Represents rated wind speed, P WT (v i ) Representing the wind speed v of the fan equipment i Test force under P WT (v i+1 ) Representing the fan equipment inWind speed v i+1 Test output;
the expression of the steady-state model of the gas turbine equipment is as follows:
wherein F is MT Representing the natural gas consumption of a gas turbine plant, P MT Represents the output electric power, eta of the gas turbine plant MT Representing the power generation efficiency of the gas turbine plant;
the expression of the steady-state model of the gas boiler equipment is as follows:
wherein F is GB Represents the natural gas consumption of the gas boiler equipment, Q GB Represents the heating power, eta of the gas boiler equipment GB Representing the corresponding heating value of the input gas boiler plant;
the expression of the steady-state model of the electric energy storage device is as follows:
0≤P ES，C ≤CapEsγ ES，C ；
0≤P ES，D ≤CapEsγ ES，D ；
W ES，min ≤W ES ≤W ES，max ；
wherein,representing the residual electric quantity of the electric energy storage device before and after charging and discharging sigma ES Representing the self-loss rate of the electric energy storage device, P ES，C Representing the charge power, eta of an electrical energy storage device ES，C Representing the charging efficiency, P, of an electrical energy storage device ES，D Represents the discharge power, eta of an electrical energy storage device ES，D Representing the discharge efficiency of an electrical energy storage device, gamma ES，C Representing the maximum charge rate of the electrical energy storage device, gamma ES，D Representing the maximum discharge rate of the electrical energy storage device, capes represents the rated capacity of the electrical energy storage device, W ES，min 、W ES，max Representing a minimum allowable capacity and a maximum allowable capacity of the electrical energy storage device;
the expression of the steady-state model of the cold accumulation equipment is as follows:
0≤Q IS，C ≤CapIsγ IS，C ；
0≤Q IS，D ≤CapIsγ IS，D ；
W IS，min ≤W IS ≤W IS，max ；
wherein,representing the residual energy, sigma, of the cold storage device before and after ice making and ice melting IS Represents the self-loss rate of the cold accumulation device, Q IS，C Ice-making power, eta, representative of cold storage device IS，C Represents ice making efficiency, Q of cold storage device IS，D Ice melting power, eta, representing cold storage device IS，D Represents ice melting efficiency, gamma of cold storage equipment IS，C Represents the maximum ice making multiplying power of cold storage equipment, gamma IS，D Representing the maximum ice-melting rate of the cold storage device, capIs representing the rated capacity of the cold storage device, W IS，min 、W IS，max Representing a minimum allowable capacity and a maximum allowable capacity of the cold storage device;
the expression of the steady-state model of the heat storage equipment is as follows:
0≤Q TS，C ≤CapTsγ TS，C ；
0≤Q TS，D ≤CapTsγ TS，D ；
W TS，min ≤W TS ≤W TS，max ；
wherein,representing the residual energy, sigma, of the heat storage device before and after heat storage and heat release TS Represents the self-loss rate of the heat storage device, Q TS，C Represents the heat storage power, eta of the heat storage device TS，C Representing the heat storage efficiency of the heat storage device, Q TS，D Represents the exothermic power, eta of the heat storage device TS，D Representing the exothermic efficiency of the thermal storage device, caps gamma TS，C Represents the maximum heat storage rate of the heat storage device, gamma TS，D Representing the maximum heat release rate of the heat storage device, capts representing the rated capacity of the heat storage device, W TS，min 、W TS，max Representing the minimum allowable capacity and the maximum allowable capacity of the heat storage device;
the expression of the steady-state model of the cooling equipment is as follows:
wherein,representing the cold power output by the cold supply device, COP EC Representing the energy efficiency ratio of the cooling equipment, P EC Represents the power consumption of the cooling device;
the expression of the steady-state model of the heating equipment is as follows:
wherein,representing the thermal power output by the electric heating device, COP HP Representing the energy efficiency ratio of the electric heating equipment, P Hp Representing the power consumption of the electric heating equipment, Q MH Represents the heat power input by the waste heat boiler, eta WH Representing the efficiency of the waste heat boiler>Representing the heat power output by the waste heat boiler;
the expression of the first objective function is:
max F 1 ＝E E -C grid -C gas ；
wherein F is 1 Represents profit and income of comprehensive energy system, E E Represents the electricity selling income of the comprehensive energy system, C grid Representing the electricity purchasing expense of the comprehensive energy system, C gas Represents the gas purchase cost of the comprehensive energy system, P grid Representing the power purchased by the comprehensive energy system from the external power grid through the power transmission line, c s Representing electricity price of electricity selling of comprehensive energy system to large power grid, c b Representing electricity price of electricity purchase of comprehensive energy system to large power grid, C gas，i Representing the gas purchase cost of the ith park in the comprehensive energy system;
the expression of the electric power balance constraint of the comprehensive energy system is as follows:
wherein P is AE，i Representing the gate power value of the ith park in the comprehensive energy system;
the expression of the gas power balance constraint of the comprehensive energy system is as follows:
wherein F is gas Represents the total amount of natural gas input by the comprehensive energy system, F MT，i Representing natural gas consumption of gas turbine plant at ith park in integrated energy system, F GB，i Representing the natural gas consumption of the gas boiler equipment of the ith park in the comprehensive energy system;
The expression of the peak power constraint of the comprehensive energy system is as follows:
P grid ≤P line，max ；
wherein P is line，max Representing the maximum allowable value of the tie line tide between the comprehensive energy system and the upper power grid;
the expression of the second objective function is:
min f i ＝C E，i +C H，i +C C，i +C gas，i ；
wherein f i Representing the operation cost of the ith park in the comprehensive energy system, C E，i Representing electricity purchase cost of ith park in comprehensive energy system, C H，i Representing the heat purchase cost of the ith park in the comprehensive energy system, C C，i Representing the cold purchase cost of the ith park in the integrated energy system,represents the time-of-use electricity price at time t, +.>Representing electricity purchasing power of ith park at t moment in comprehensive energy system, c H Represents the heat supply price, Q H，i Representing the heat purchasing power of the ith park in the comprehensive energy system, c C Represents the price of cooling, Q C，i Representing the cold purchase power of the ith park in the integrated energy system, < >>Represents the price of natural gas at time t, +.>Representing the power generated by the gas turbine plant at time t, η of the ith park in the integrated energy system MT Representing the efficiency of a gas turbine plant, +.>Heat production power eta of gas boiler equipment representing ith park in comprehensive energy system at t moment GFB Representing the efficiency of the gas boiler plant;
the park electric power balance constraint expression is:
P AE，i +P MT，i +P PV，i +P WT，i +P ES，D，i ＝P Eload，i +P HP，i +P ES，C，i +P EC，i ；
Wherein the method comprises the steps of，P MT，i Representing the power generated by the gas turbine plant of the ith park in the integrated energy system, P PV，i Representing the power generated by photovoltaic devices in the ith park in the integrated energy system, P WT，i Representing the power generated by fan equipment of the ith park in the integrated energy system, P Eload，i Representing the electric load power of the ith park in the integrated energy system, P HP，i Representing the power consumption of the electric heating equipment of the ith park in the comprehensive energy system, P ES，C，i Representing the charging power of the electric energy storage device of the ith park in the integrated energy system, P ES，D，i Representing the discharge power, P, of an electric energy storage device of an ith park in an integrated energy system EC，i Representing the power consumption of the electric refrigeration equipment of the ith park in the comprehensive energy system;
the park thermal power balance constraint expression is:
wherein Q is GB，i Representing the heat power output by gas boiler equipment of ith park in comprehensive energy system, Q TS，D，i Representing the thermal power output by the heat storage equipment of the ith park in the integrated energy system, Q TSC，i Representing the thermal power input by the heat storage equipment of the ith park in the integrated energy system,representing the thermal power output by the electric heating equipment of the ith park in the integrated energy system,representing the heat power input by the waste heat boiler equipment of the ith park in the comprehensive energy system, Q HL，i Representing the heat load power of an ith park in the comprehensive energy system;
the park cold power balance constraint expression is:
wherein,represents the cold power output by the electric refrigeration equipment of the ith park in the comprehensive energy system, Q IS，D，i Ice melting power, Q, representing heat storage equipment of ith park in integrated energy system IS，C，i An ice making power input by a heat storage device representing an ith park in the integrated energy system, Q CL，i Representing the cooling load power of an ith park in the comprehensive energy system;
the expression of the energy equipment operation constraint is as follows:
wherein P is min，k Represents the minimum operation allowable electric power, P, of the kth energy device in the integrated energy system max，k Represents the maximum operation allowable electric power of the kth energy device in the integrated energy system, Q min，k Represents the minimum operation allowable thermal power, Q, of the kth energy device in the integrated energy system max，k Represents the maximum operation allowable thermal power, W, of the kth energy device in the integrated energy system min，k Representing the minimum allowable operation capacity, W, of the kth energy device in the integrated energy system max，k Representing the maximum allowable operating capacity of the kth energy device in the integrated energy system.
2. A computer device comprising a memory and a processor, the memory storing a computer program, characterized in that the processor, when executing the computer program, implements the steps of the coordinated control method for an integrated energy system according to claim 1.
3. A computer readable storage medium having stored thereon a computer program, characterized in that the computer program when executed by a processor realizes the steps of the cooperative control method for an integrated energy system according to claim 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211261103.4A CN115619145B (en) | 2022-10-14 | 2022-10-14 | Cooperative control method and device for comprehensive energy system and computer equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211261103.4A CN115619145B (en) | 2022-10-14 | 2022-10-14 | Cooperative control method and device for comprehensive energy system and computer equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115619145A CN115619145A (en) | 2023-01-17 |
| CN115619145B true CN115619145B (en) | 2024-03-19 |
Family
ID=84863071
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211261103.4A Active CN115619145B (en) | 2022-10-14 | 2022-10-14 | Cooperative control method and device for comprehensive energy system and computer equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115619145B (en) |
Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109245093A (en) * | 2018-09-28 | 2019-01-18 | 上海电力学院 | A kind of supply of cooling, heating and electrical powers distributed busbar protection collaboration Optimization Scheduling |
| WO2019134532A1 (en) * | 2018-01-02 | 2019-07-11 | 佛山科学技术学院 | Operating method of combined heat and power system containing wind and light renewable energy |
| CN110333660A (en) * | 2019-07-29 | 2019-10-15 | 西安科技大学 | A kind of cooling heating and power generation system Multipurpose Optimal Method |
| CN110990785A (en) * | 2019-11-27 | 2020-04-10 | 江苏方天电力技术有限公司 | Multi-objective-based optimal scheduling method for intelligent park comprehensive energy system |
| CN111342451A (en) * | 2020-03-12 | 2020-06-26 | 四川大学 | Park comprehensive energy system economic configuration method for promoting renewable energy consumption |
| CN112068436A (en) * | 2020-09-16 | 2020-12-11 | 国网山东省电力公司经济技术研究院 | Layered and distributed control method and system for comprehensive energy system of industrial park |
| CN112383053A (en) * | 2020-11-19 | 2021-02-19 | 东北大学 | Comprehensive energy system optimization scheduling method based on improved particle swarm optimization |
| CN112465240A (en) * | 2020-12-03 | 2021-03-09 | 上海电力大学 | Cooperative game-based multi-park energy scheduling optimization method for comprehensive energy system |
| CN112600253A (en) * | 2020-11-18 | 2021-04-02 | 国网浙江省电力有限公司台州供电公司 | Park comprehensive energy collaborative optimization method and equipment based on optimal energy utilization efficiency |
| CN112907098A (en) * | 2021-03-08 | 2021-06-04 | 国网江苏省电力有限公司经济技术研究院 | Multi-stage capacity configuration method and configuration system of park comprehensive energy system |
| CN113393126A (en) * | 2021-06-16 | 2021-09-14 | 沈阳工程学院 | High-energy-consumption park and power grid alternative parallel cooperative optimization scheduling method |
| CN113410873A (en) * | 2021-05-19 | 2021-09-17 | 南京华盾电力信息安全测评有限公司 | Regional comprehensive energy system configuration method considering source load multiple uncertainty |
| CN113469430A (en) * | 2021-06-29 | 2021-10-01 | 华电电力科学研究院有限公司 | Multi-energy complementary capacity configuration method of integrated energy park |
| CN113723793A (en) * | 2021-08-25 | 2021-11-30 | 国网江苏省电力工程咨询有限公司 | Method, device, equipment and medium for realizing park comprehensive energy system |
| CN113762632A (en) * | 2021-09-10 | 2021-12-07 | 国网四川省电力公司经济技术研究院 | Collaborative optimization operation method and system of electrical comprehensive energy system |
| CN113762708A (en) * | 2021-07-01 | 2021-12-07 | 国网江西省电力有限公司赣州供电分公司 | Park level comprehensive energy system planning method considering multi-target cooperation |
| CN113987734A (en) * | 2021-07-22 | 2022-01-28 | 太原理工大学 | Multi-objective optimization scheduling method for park comprehensive energy system under opportunity constraint condition |
| CN114066204A (en) * | 2021-11-11 | 2022-02-18 | 内蒙古科技大学 | Integrated optimization planning and operation method and device of comprehensive energy system |
| CN114565245A (en) * | 2022-02-16 | 2022-05-31 | 石河子大学 | Comprehensive energy system optimization scheduling method and device considering electricity, cold and heat |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7261507B2 (en) * | 2020-09-04 | 2023-04-20 | ノース チャイナ エレクトリック パワー ユニバーシティー | Electric heat pump - regulation method and system for optimizing cogeneration systems |
-
2022
- 2022-10-14 CN CN202211261103.4A patent/CN115619145B/en active Active
Patent Citations (19)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2019134532A1 (en) * | 2018-01-02 | 2019-07-11 | 佛山科学技术学院 | Operating method of combined heat and power system containing wind and light renewable energy |
| CN109245093A (en) * | 2018-09-28 | 2019-01-18 | 上海电力学院 | A kind of supply of cooling, heating and electrical powers distributed busbar protection collaboration Optimization Scheduling |
| CN110333660A (en) * | 2019-07-29 | 2019-10-15 | 西安科技大学 | A kind of cooling heating and power generation system Multipurpose Optimal Method |
| CN110990785A (en) * | 2019-11-27 | 2020-04-10 | 江苏方天电力技术有限公司 | Multi-objective-based optimal scheduling method for intelligent park comprehensive energy system |
| CN111342451A (en) * | 2020-03-12 | 2020-06-26 | 四川大学 | Park comprehensive energy system economic configuration method for promoting renewable energy consumption |
| CN112068436A (en) * | 2020-09-16 | 2020-12-11 | 国网山东省电力公司经济技术研究院 | Layered and distributed control method and system for comprehensive energy system of industrial park |
| CN112600253A (en) * | 2020-11-18 | 2021-04-02 | 国网浙江省电力有限公司台州供电公司 | Park comprehensive energy collaborative optimization method and equipment based on optimal energy utilization efficiency |
| CN112383053A (en) * | 2020-11-19 | 2021-02-19 | 东北大学 | Comprehensive energy system optimization scheduling method based on improved particle swarm optimization |
| CN112465240A (en) * | 2020-12-03 | 2021-03-09 | 上海电力大学 | Cooperative game-based multi-park energy scheduling optimization method for comprehensive energy system |
| CN112907098A (en) * | 2021-03-08 | 2021-06-04 | 国网江苏省电力有限公司经济技术研究院 | Multi-stage capacity configuration method and configuration system of park comprehensive energy system |
| CN113410873A (en) * | 2021-05-19 | 2021-09-17 | 南京华盾电力信息安全测评有限公司 | Regional comprehensive energy system configuration method considering source load multiple uncertainty |
| CN113393126A (en) * | 2021-06-16 | 2021-09-14 | 沈阳工程学院 | High-energy-consumption park and power grid alternative parallel cooperative optimization scheduling method |
| CN113469430A (en) * | 2021-06-29 | 2021-10-01 | 华电电力科学研究院有限公司 | Multi-energy complementary capacity configuration method of integrated energy park |
| CN113762708A (en) * | 2021-07-01 | 2021-12-07 | 国网江西省电力有限公司赣州供电分公司 | Park level comprehensive energy system planning method considering multi-target cooperation |
| CN113987734A (en) * | 2021-07-22 | 2022-01-28 | 太原理工大学 | Multi-objective optimization scheduling method for park comprehensive energy system under opportunity constraint condition |
| CN113723793A (en) * | 2021-08-25 | 2021-11-30 | 国网江苏省电力工程咨询有限公司 | Method, device, equipment and medium for realizing park comprehensive energy system |
| CN113762632A (en) * | 2021-09-10 | 2021-12-07 | 国网四川省电力公司经济技术研究院 | Collaborative optimization operation method and system of electrical comprehensive energy system |
| CN114066204A (en) * | 2021-11-11 | 2022-02-18 | 内蒙古科技大学 | Integrated optimization planning and operation method and device of comprehensive energy system |
| CN114565245A (en) * | 2022-02-16 | 2022-05-31 | 石河子大学 | Comprehensive energy system optimization scheduling method and device considering electricity, cold and heat |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115619145A (en) | 2023-01-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109659927B (en) | Energy storage capacity configuration method of comprehensive energy microgrid considering energy storage participation degree | |
| CN102710013B (en) | Park energy-network energy optimizing management system based on microgrids and implementing method thereof | |
| CN111881616A (en) | Operation optimization method of comprehensive energy system based on multi-subject game | |
| CN106022503A (en) | Micro-grid capacity programming method meeting coupling type electric cold and heat demand | |
| CN106950840A (en) | Towards the integrated energy system layered distribution type control method for coordinating of power network peak clipping | |
| CN109634119B (en) | Energy internet optimization control method based on rolling optimization in day | |
| CN109523065A (en) | A kind of micro- energy net Optimization Scheduling based on improvement quanta particle swarm optimization | |
| CN109636056A (en) | A kind of multiple-energy-source microgrid decentralization Optimization Scheduling based on multi-agent Technology | |
| CN108400593A (en) | Active distribution network electrical model method for building up based on layering multi-agent technology | |
| CN112163968B (en) | Optimized scheduling method for comprehensive energy storage system | |
| CN113471976B (en) | Optimal scheduling method based on multi-energy complementary micro-grid and active power distribution network | |
| CN111737884A (en) | Multi-target random planning method for micro-energy network containing multiple clean energy sources | |
| CN114707289A (en) | Opportunity constraint-based multi-objective optimization method for electrothermal coupling comprehensive energy system | |
| CN107302231A (en) | A kind of small hydropower system accesses the random economic load dispatching method of two targets of power network | |
| CN105956693B (en) | A kind of virtual plant economic load dispatching method based on distributed gradient algorithm | |
| CN111030101B (en) | Clean energy consumption linkage regulation and control method and system based on diversified big data | |
| CN112560221A (en) | Capacity distribution method and device for facility agriculture energy network containing enhanced geothermal system | |
| CN116822831A (en) | Micro-energy network group optimization planning method containing shared energy storage system | |
| CN115149540A (en) | Multi-region virtual power plant collaborative optimization operation method based on partition autonomy | |
| Han et al. | Multi-View clustering and discrete consensus based tri-level coordinated control of wind farm and adiabatic compressed air energy storage for providing frequency regulation service | |
| CN107730032A (en) | A kind of electric cold and hot trilogy supply Optimal Scheduling of pan-energy network and method | |
| CN108736518B (en) | Comprehensive energy supply system and method for urban complex and large public building group | |
| CN115619145B (en) | Cooperative control method and device for comprehensive energy system and computer equipment | |
| CN117081143A (en) | Method for promoting coordination and optimization operation of park comprehensive energy system for distributed photovoltaic on-site digestion | |
| CN116911437A (en) | Data center energy system optimal scheduling method considering load response characteristics |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |