CN115034551A - Energy interaction optimization method and device for multi-park comprehensive energy system - Google Patents

Energy interaction optimization method and device for multi-park comprehensive energy system Download PDF

Info

Publication number
CN115034551A
CN115034551A CN202210440078.XA CN202210440078A CN115034551A CN 115034551 A CN115034551 A CN 115034551A CN 202210440078 A CN202210440078 A CN 202210440078A CN 115034551 A CN115034551 A CN 115034551A
Authority
CN
China
Prior art keywords
energy
park
energy system
ith
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.)
Pending
Application number
CN202210440078.XA
Other languages
Chinese (zh)
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.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Economic and Technological Research Institute of State Grid Fujian Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Economic and Technological Research Institute of State Grid Fujian Electric Power 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 State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, Economic and Technological Research Institute of State Grid Fujian Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202210440078.XA priority Critical patent/CN115034551A/en
Publication of CN115034551A publication Critical patent/CN115034551A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06312Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315Needs-based resource requirements planning or analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P80/00Climate change mitigation technologies for sector-wide applications
    • Y02P80/10Efficient use of energy, e.g. using compressed air or pressurized fluid as energy carrier

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Tourism & Hospitality (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • Primary Health Care (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention relates to the technical field of comprehensive energy systems, and particularly provides a method and a device for optimizing energy interaction of a multi-park comprehensive energy system, wherein the method comprises the following steps: acquiring the load capacity of equipment of each park comprehensive energy system in the multi-park comprehensive energy system energy sharing system in the energy interaction process; substituting the load amount into a pre-constructed multi-park comprehensive energy system energy interaction optimization model and solving to obtain an optimization result; and obtaining an energy interaction optimization scheme of the multi-park comprehensive energy system based on the optimization result. The technical scheme provided by the invention realizes the cooperative optimization of the energy supply side and the energy demand side of the comprehensive energy systems of a plurality of parks, and improves the local consumption capacity of new energy output and the utilization efficiency of various energy sources, thereby realizing the overall optimal operation scheduling of the comprehensive energy service provider and the comprehensive energy systems of the plurality of parks.

Description

Multi-park comprehensive energy system energy interaction optimization method and device
Technical Field
The invention relates to the technical field of comprehensive energy systems, in particular to a method and a device for optimizing energy interaction of a multi-park comprehensive energy system.
Background
The comprehensive energy system is the organic coordination of energy network, energy production, conversion, storage, consumption and end use participants, and is a complex multi-energy coupling system. Park-level integrated energy system (PIES) systems are usually located at the end of an energy network to meet the energy demand of end users in a park, and are widely applied to the fields of industrial/commercial parks, schools, buildings, residential communities and the like. Due to the diversified energy requirements of a plurality of users in the comprehensive energy system of the park, the energy sources of different types are closely coupled and flexibly converted, and great challenges are brought to the cooperative operation and energy management of the system.
The comprehensive energy systems in different parks have different new energy output levels, load consumption capacities and energy response characteristics, and the problems of low energy utilization level, small adjustable potential, high operation cost and the like can occur during independent operation. With the rise of sharing economy, energy sharing in adjacent parks is becoming an effective means for improving energy utilization efficiency. Therefore, a method for optimizing energy interaction of a multi-park integrated energy system is needed.
Disclosure of Invention
In order to overcome the defects, the invention provides a method and a device for optimizing energy interaction of a multi-park comprehensive energy system.
In a first aspect, a method for optimizing energy interaction of a multi-park integrated energy system is provided, where the method includes:
acquiring the load capacity of equipment of each park integrated energy system in the multi-park integrated energy system energy sharing system in the energy interaction process;
substituting the load amount into a pre-constructed multi-park comprehensive energy system energy interaction optimization model and solving to obtain an optimization result;
and obtaining an energy interaction optimization scheme of the multi-park comprehensive energy system based on the optimization result.
Preferably, the load amount includes at least one of: electrical load, gas load, thermal load, cold load.
Preferably, the energy interaction optimization scheme of the multi-park integrated energy system comprises at least one of the following: EB power consumption power of the park comprehensive energy system in a time period t, EC equipment power consumption power of the park comprehensive energy system in the time period t, electric load power participating in comprehensive demand response, electric energy storage discharge power of the park comprehensive energy system in the time period t, electric energy sharing quantity of the park comprehensive energy system in the time period t, AC equipment heat consumption power of the park comprehensive energy system in the time period t, heat load power participating in comprehensive demand response, heat energy sharing quantity of the park comprehensive energy system in the time period t, CCHP gas consumption power of the park comprehensive energy system in the time period t, gas load power participating in comprehensive demand response and cold energy sharing quantity of the park comprehensive energy system in the time period t.
Further, the pre-constructed multi-park integrated energy system energy interaction optimization model comprises: and an objective function and a constraint condition are established for the energy interaction optimization of the multi-park comprehensive energy system.
Further, the objective function is calculated as follows:
Figure BDA0003613528410000021
in the above formula, C is the total cost of energy interaction optimization of the park comprehensive energy system,
Figure BDA0003613528410000022
in order to integrate the operating costs of the energy service provider,
Figure BDA0003613528410000023
the operation cost of the ith park comprehensive energy system is shown, and I is the total number of the park comprehensive energy systems.
Further, the mathematical model of the constraint condition is as follows:
Figure BDA0003613528410000024
Figure BDA0003613528410000025
Figure BDA0003613528410000026
Figure BDA0003613528410000027
Figure BDA0003613528410000028
Figure BDA0003613528410000029
in the above formula, the first and second carbon atoms are,
Figure BDA00036135284100000210
for the net demand of electrical energy for the ith campus complex energy system during time period t,
Figure BDA00036135284100000211
the electric energy balance of the ith park comprehensive energy system in the time period t, P e,i,t Electric load of the i-th park integrated energy system in the period of t, P PV,i,t The photovoltaic output level P of the ith park integrated energy system in the period of t WT,i,t For the wind power output level of the ith campus comprehensive energy system in the time period t,
Figure BDA00036135284100000212
the CCHP electricity generation power of the ith park integrated energy system in the time period t,
Figure BDA00036135284100000213
the power consumption of the EB of the ith park comprehensive energy system in the time period t,
Figure BDA00036135284100000214
the power consumption of EC equipment in the ith park comprehensive energy system in the period of t is P IDR,t For the time period t the electrical load power participating in the integrated demand response,
Figure BDA00036135284100000215
the electric energy storage and discharge power of the ith park comprehensive energy system in the time period t,
Figure BDA00036135284100000216
electric energy sharing amount, P, of the ith park integrated energy system in the period of t h,i,t For the heat load of the ith campus integrated energy system during the period t,
Figure BDA00036135284100000217
for the CCHP heat-generating power of the ith park integrated energy system in the period t,
Figure BDA0003613528410000031
is the ith parkThe EB of the comprehensive energy system generates heat power in a period t,
Figure BDA0003613528410000032
for the heat consumption power of the AC equipment of the ith park comprehensive energy system in the time period t, H IDR,t The thermal load power participating in the integrated demand response for the period t,
Figure BDA0003613528410000033
the heat energy sharing amount of the ith park integrated energy system in the time period t,
Figure BDA0003613528410000034
for the net demand of gas energy of the ith park energy system at time t,
Figure BDA0003613528410000035
for the CCHP gas consumption power, P, of the ith park integrated energy system in the period of t g,i,t Air load of the i-th park integrated energy system at time t, G IDR,t For the time period t the gas load power participating in the integrated demand response,
Figure BDA0003613528410000036
cooling power of EC equipment of the ith campus integrated energy system in the period of t,
Figure BDA0003613528410000037
cooling power of AC equipment in time t period for ith park comprehensive energy system c,i,t For the cooling load of the ith campus energy system during the time period t,
Figure BDA0003613528410000038
for the cooling energy sharing amount of the ith campus comprehensive energy system in the time period t,
Figure BDA0003613528410000039
is the energy sharing upper limit of m types of energy sources in the unit time interval of the comprehensive energy system of the ith park,
Figure BDA00036135284100000310
the shared energy of the m types of energy sources in the unit time interval of the comprehensive energy system of the ith park is obtained.
Further, after obtaining the energy interaction optimization scheme of the multi-park integrated energy system based on the optimization result, the method includes:
solving a pre-constructed maximum utility model to obtain a subsidy result;
and obtaining the energy sharing excitation subsidy of each park integrated energy system in the multi-park integrated energy system energy sharing system based on the subsidy result.
Further, the pre-constructed maximum utility model comprises: and the maximum utility objective function and the constraint condition are used for exciting the comprehensive energy systems of all the parks to participate in energy sharing.
Further, the maximum utility objective function is calculated as follows:
Figure BDA00036135284100000311
in the above formula, f is the maximum utility objective function value, Δ IESP For the economic benefit of the energy sharing of the comprehensive energy service provider,
Figure BDA00036135284100000312
for the economic benefit of the i-th park comprehensive energy system participating in energy sharing, omega i And sigma is the return coefficient of the ith park integrated energy system.
Further, the calculation formula of the economic benefit of the multi-park integrated energy system energy sharing system participating in energy sharing is as follows:
Figure BDA00036135284100000313
the calculation formula of the economic benefit of the ith park integrated energy system participating in energy sharing is as follows:
Figure BDA0003613528410000041
the calculation formula of the market contribution degree of the ith park integrated energy system is as follows:
Figure BDA0003613528410000042
in the above formula, the first and second carbon atoms are,
Figure BDA0003613528410000043
net operating cost, ζ, for the integrated energy facilitator after participation in energy sharing i To be subsidized for energy sharing and excitation of the ith park comprehensive energy system,
Figure BDA0003613528410000044
in order to participate in the operation cost of the ith campus comprehensive energy system after energy sharing,
Figure BDA0003613528410000045
and the shared value of the ith campus comprehensive energy system in the time period T, wherein T is the total time period.
Further, the calculation formula of the shared value of the ith park integrated energy system in the time period t is as follows:
Figure BDA0003613528410000046
in the above formula, the first and second carbon atoms are,
Figure BDA0003613528410000047
and
Figure BDA0003613528410000048
the local marginal prices for electrical, thermal, and cold loads, respectively.
In a second aspect, an energy interaction optimization device for a multi-park integrated energy system is provided, which includes:
the acquisition module is used for acquiring the load of equipment of each park integrated energy system in the multi-park integrated energy system energy sharing system in the energy interaction process;
the solving module is used for substituting the load into a pre-constructed multi-park comprehensive energy system energy interaction optimization model and solving to obtain an optimization result;
and the analysis module is used for obtaining an energy interaction optimization scheme of the multi-park comprehensive energy system based on the optimization result.
Preferably, the load amount includes at least one of: electrical load, gas load, thermal load, cold load.
Preferably, the energy interaction optimization scheme of the multi-park integrated energy system comprises at least one of the following: EB power consumption power of the park comprehensive energy system in a time period t, EC equipment power consumption power of the park comprehensive energy system in the time period t, electric load power participating in comprehensive demand response, electric energy storage discharge power of the park comprehensive energy system in the time period t, electric energy sharing quantity of the park comprehensive energy system in the time period t, AC equipment heat consumption power of the park comprehensive energy system in the time period t, heat load power participating in comprehensive demand response, heat energy sharing quantity of the park comprehensive energy system in the time period t, CCHP gas consumption power of the park comprehensive energy system in the time period t, gas load power participating in comprehensive demand response and cold energy sharing quantity of the park comprehensive energy system in the time period t.
Further, the pre-constructed multi-park integrated energy system energy interaction optimization model comprises: and an objective function and a constraint condition are established for the energy interaction optimization of the multi-park comprehensive energy system.
Further, the objective function is calculated as follows:
Figure BDA0003613528410000051
in the above formula, C is the total cost of energy interaction optimization of the park comprehensive energy system,
Figure BDA0003613528410000052
in order to integrate the operating costs of the energy service provider,
Figure BDA0003613528410000053
the operation cost of the ith park comprehensive energy system is shown, and I is the total number of the park comprehensive energy systems.
Further, the mathematical model of the constraint condition is as follows:
Figure BDA0003613528410000054
Figure BDA0003613528410000055
Figure BDA0003613528410000056
Figure BDA0003613528410000057
Figure BDA0003613528410000058
Figure BDA0003613528410000059
in the above formula, the first and second carbon atoms are,
Figure BDA00036135284100000510
for the net demand of electrical energy for the ith campus complex energy system during time period t,
Figure BDA00036135284100000511
the electric energy balance of the ith park comprehensive energy system in the time period t is ensured,P e,i,t electric load of the i-th park integrated energy system in the period of t, P PV,i,t The photovoltaic output level P of the ith park integrated energy system in the period of t WT,i,t For the wind power output level of the ith park comprehensive energy system in the time period t,
Figure BDA00036135284100000512
the CCHP electricity generation power of the ith park integrated energy system in the time period t,
Figure BDA00036135284100000513
the EB of the ith park comprehensive energy system consumes power in the time period t,
Figure BDA00036135284100000514
the power consumption of EC equipment in the ith park comprehensive energy system in the period of t is P IDR,t For the time period t the electrical load power participating in the integrated demand response,
Figure BDA00036135284100000515
the electric energy storage and discharge power of the ith park comprehensive energy system in the time period t,
Figure BDA00036135284100000516
the electric energy sharing quantity P of the i-th park comprehensive energy system in the time period t h,i,t For the thermal load of the ith campus energy system during time t,
Figure BDA00036135284100000517
for the CCHP heat-generating power of the ith park integrated energy system in the period t,
Figure BDA00036135284100000518
for the EB heat generation power of the ith park comprehensive energy system in the period of t,
Figure BDA00036135284100000519
for the heat consumption power of the AC equipment of the ith park comprehensive energy system in the time period t, H IDR,t Participating in the Hot Duck of Integrated demand response for time tThe power of the load is controlled by the power of the load,
Figure BDA00036135284100000520
the heat energy sharing amount of the ith park integrated energy system in the time period t,
Figure BDA00036135284100000521
for the net demand of gas energy of the ith park energy system at time t,
Figure BDA00036135284100000522
for the CCHP gas consumption power, P, of the ith district comprehensive energy system in the period of t g,i,t Gas load of the i-th park integrated energy system at time t, G IDR,t For the time period t the gas load power participating in the integrated demand response,
Figure BDA0003613528410000061
for the cooling power of the EC equipment of the ith park integrated energy system in the t period,
Figure BDA0003613528410000062
AC equipment refrigeration power, P, for the ith park energy system during the period t c,i,t For the cooling load of the ith campus energy system during the time period t,
Figure BDA0003613528410000063
for the cooling energy sharing amount of the ith campus comprehensive energy system in the time period t,
Figure BDA0003613528410000064
is the energy sharing upper limit of m types of energy sources in the unit time interval of the comprehensive energy system of the ith park,
Figure BDA0003613528410000065
the shared energy of the m types of energy sources in the unit time interval of the comprehensive energy system of the ith park is obtained.
In a third aspect, a computer device is provided, comprising: one or more processors;
the processor to store one or more programs;
when the one or more programs are executed by the one or more processors, the method for energy interactive optimization of the multi-campus integrated energy system is implemented.
In a fourth aspect, a computer readable storage medium is provided, on which a computer program is stored, which, when executed, implements the method for energy interactive optimization of a multi-campus integrated energy system.
One or more technical schemes of the invention at least have one or more of the following beneficial effects:
the invention provides a method and a device for optimizing energy interaction of a multi-park comprehensive energy system, wherein the method comprises the following steps: acquiring the load capacity of equipment of each park comprehensive energy system in the multi-park comprehensive energy system energy sharing system in the energy interaction process; substituting the load amount into a pre-constructed multi-park comprehensive energy system energy interaction optimization model and solving to obtain an optimization result; and obtaining an energy interaction optimization scheme of the multi-park comprehensive energy system based on the optimization result. The technical scheme provided by the invention has the advantages that by means of the difference of net loads among different parks, each independent park is connected based on the comprehensive energy service provider, the interaction and sharing of various energy sources among the parks are promoted, the energy operation scheme of the comprehensive energy system of each park is favorably optimized in a combined manner, the local consumption capacity of new energy sources and the utilization efficiency of various energy equipment are improved, and the operation cost of the comprehensive energy service provider and the operation cost of the comprehensive energy system of the park are reduced.
Drawings
FIG. 1 is a schematic flow chart of the main steps of the energy interaction optimization method of the multi-park integrated energy system according to the embodiment of the invention;
FIG. 2 is a block diagram of an embodiment of the invention for an energy sharing system for a multi-campus integrated energy system;
figure 3 is a diagram of an exemplary architecture of a campus complex energy system according to an embodiment of the present invention;
fig. 4 is a main structural block diagram of the energy interaction optimization device of the multi-park integrated energy system according to the embodiment of the invention.
Detailed Description
The following describes embodiments of the present invention in further detail with reference to the accompanying drawings.
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As introduced in the background, an integrated energy system is an organic coordination of energy network, energy production, conversion, storage, consumption and end use participants, and is a complex multi-energy coupled system. Park-level integrated energy system (PIES) systems are usually located at the end of an energy network to meet the energy demand of end users in a park, and are widely applied to the fields of industrial/commercial parks, schools, buildings, residential communities and the like. Due to the diversified energy requirements of a plurality of users in the comprehensive energy system of the park, the energy sources of different types are closely coupled and flexibly converted, and great challenges are brought to the cooperative operation and energy management of the system.
The comprehensive energy systems in different parks have different new energy output levels, load consumption capacities and energy response characteristics, and the problems of low energy utilization level, small adjustable potential, high operation cost and the like can occur during independent operation. With the rise of sharing economy, energy sharing in adjacent parks is becoming an effective means for improving energy utilization efficiency. Therefore, a method for optimizing the energy interaction of the multi-park integrated energy system is needed.
In order to solve the problems, the invention provides an energy interaction optimization method for a multi-park integrated energy system, which is used for acquiring the load capacity of equipment of each park integrated energy system in an energy sharing system of the multi-park integrated energy system in the energy interaction process; substituting the load of the equipment of each park integrated energy system in the multi-park integrated energy system energy sharing system in the energy interaction process into a pre-constructed multi-park integrated energy system energy interaction optimization model and solving to obtain an optimization result; and obtaining an energy interaction optimization scheme of the multi-park integrated energy system based on the optimization result, so that the cooperative optimization of the energy supply side and the energy demand side of the multi-park integrated energy system is realized, the local consumption capability of new energy output and the utilization efficiency of various energy sources are improved, and the overall optimal operation scheduling of the integrated energy service provider and the multi-park integrated energy system is realized.
The energy interactive optimization method of the multi-park comprehensive energy system is described in detail below.
Example 1
Referring to fig. 1, fig. 1 is a schematic flow chart of main steps of a method for energy interactive optimization of a multi-park integrated energy system according to an embodiment of the present invention. As shown in fig. 1, the method for optimizing energy interaction of a multi-park integrated energy system in the embodiment of the present invention mainly includes the following steps:
step S101: acquiring the load capacity of equipment of each park comprehensive energy system in the multi-park comprehensive energy system energy sharing system in the energy interaction process;
step S102: substituting the load amount into a pre-constructed multi-park comprehensive energy system energy interaction optimization model and solving to obtain an optimization result;
step S103: and obtaining an energy interaction optimization scheme of the multi-park comprehensive energy system based on the optimization result.
Wherein the load amount comprises at least one of: electrical load, gas load, thermal load, cold load.
In this embodiment, as shown in fig. 2, the energy sharing system of the multi-park integrated energy system mainly includes an external energy network, an Integrated Energy Service Provider (IESP), and a plurality of PIES. Wherein there is an energy interaction between the IESP and the external energy network and the PIES, and also between the PIES.
Fig. 3 is a typical architecture of the campus integrated energy system used in the present invention, which includes three components of an energy network, a campus side device and a user side load according to an energy flow direction. The park side equipment comprises distributed new energy equipment, energy coupling equipment and energy storage equipment, and relates to links such as energy production, conversion and storage. The customer side loads include electrical, thermal, gas and cold loads for participating in the integrated demand response.
It should be noted that the schematic diagram shown in fig. 3 is only an exemplary illustration of the embodiment of the present invention, and is not limited to the integrated energy system architecture.
Specifically, the energy interaction optimization scheme of the multi-park integrated energy system comprises at least one of the following: EB power consumption power of the park comprehensive energy system in a time period t, EC equipment power consumption power of the park comprehensive energy system in the time period t, electric load power participating in comprehensive demand response, electric energy storage discharge power of the park comprehensive energy system in the time period t, electric energy sharing quantity of the park comprehensive energy system in the time period t, AC equipment heat consumption power of the park comprehensive energy system in the time period t, heat load power participating in comprehensive demand response, heat energy sharing quantity of the park comprehensive energy system in the time period t, CCHP gas consumption power of the park comprehensive energy system in the time period t, gas load power participating in comprehensive demand response and cold energy sharing quantity of the park comprehensive energy system in the time period t.
In this embodiment, the IESP will associate the PIES together to perform the shared optimal scheduling of the surplus energy. Energy sharing can deviate the PIES from the individual optimal scheduling to meet the energy demand of the adjacent PIES, so the IESP needs to allocate sharing benefits to the PIES to enable the PIES to share energy, therefore, the pre-constructed multi-park integrated energy system energy interaction optimization model comprises the following steps: and an objective function and a constraint condition are established for the energy interaction optimization of the multi-park comprehensive energy system.
Further, the objective function is calculated as follows:
Figure BDA0003613528410000081
in the above formula, C is a circleThe total cost of energy interaction optimization of the regional integrated energy system,
Figure BDA0003613528410000082
in order to integrate the operating costs of the energy service provider,
Figure BDA0003613528410000083
the operation cost of the ith park comprehensive energy system is shown, and I is the total number of the park comprehensive energy systems.
Further, the mathematical model of the constraint is as follows:
Figure BDA0003613528410000084
Figure BDA0003613528410000085
Figure BDA0003613528410000091
Figure BDA0003613528410000092
Figure BDA0003613528410000093
Figure BDA0003613528410000094
in the above formula, the first and second carbon atoms are,
Figure BDA0003613528410000095
for the net demand of electrical energy for the ith campus complex energy system during time period t,
Figure BDA0003613528410000096
the electric energy balance of the ith park comprehensive energy system in the time period t is P e,i,t Electric load, P, for the ith park energy system during the period t PV,i,t The photovoltaic output level P of the ith park integrated energy system in the period of t WT,i,t For the wind power output level of the ith park comprehensive energy system in the time period t,
Figure BDA0003613528410000097
the CCHP electricity generation power of the ith park integrated energy system in the time period t,
Figure BDA0003613528410000098
the EB of the ith park comprehensive energy system consumes power in the time period t,
Figure BDA0003613528410000099
the power consumption of EC equipment in the ith park comprehensive energy system in the period of t is P IDR,t For the time period t the electrical load power participating in the integrated demand response,
Figure BDA00036135284100000910
the electric energy storage and discharge power of the ith park comprehensive energy system in the time period t,
Figure BDA00036135284100000911
the electric energy sharing quantity P of the i-th park comprehensive energy system in the time period t h,i,t For the heat load of the ith campus integrated energy system during the period t,
Figure BDA00036135284100000912
for the CCHP heat-generating power of the ith park integrated energy system in the period t,
Figure BDA00036135284100000913
for the EB heat production power of the ith campus comprehensive energy system in the period of t,
Figure BDA00036135284100000914
consuming heat for AC equipment of ith park integrated energy system in t periodPower, H IDR,t The thermal load power participating in the integrated demand response for the period t,
Figure BDA00036135284100000915
the heat energy sharing amount of the ith park integrated energy system in the time period t,
Figure BDA00036135284100000916
for the net demand of gas energy of the ith park energy system at time t,
Figure BDA00036135284100000917
for the CCHP gas consumption power, P, of the ith park integrated energy system in the period of t g,i,t Gas load of the i-th park integrated energy system at time t, G IDR,t For the time period t the gas load power participating in the integrated demand response,
Figure BDA00036135284100000918
for the cooling power of the EC equipment of the ith park integrated energy system in the t period,
Figure BDA00036135284100000919
cooling power of AC equipment in time t period for ith park comprehensive energy system c,i,t For the cooling load of the ith campus energy system during the time period t,
Figure BDA00036135284100000920
for the cold energy sharing amount of the ith park integrated energy system in the time period t,
Figure BDA00036135284100000921
is the energy sharing upper limit of m types of energy sources in the unit time interval of the comprehensive energy system of the ith park,
Figure BDA00036135284100000922
the shared energy of the m types of energy sources in the unit time interval of the comprehensive energy system of the ith park is obtained.
By the step, an energy scheduling scheme and the overall operation cost of the comprehensive energy system of each park can be obtained after the energy sharing is considered, at the moment, the scheduling plan of each energy device is changed to meet the overall optimization, and at the moment, the economic benefit obtained by the energy sharing needs to be redistributed to each park so as to encourage the parks to participate in the energy sharing.
Further, since energy sharing requires that the campus deviate from the individual optimal schedule, thereby increasing individual costs, a fair and efficient benefit allocation mechanism is needed to encourage the campus to participate in energy sharing. The main idea is to distribute sharing benefits according to market contribution degrees of parks participating in energy sharing, and therefore, after obtaining an energy interaction optimization scheme of the multi-park integrated energy system based on the optimization result, the method includes:
solving a pre-constructed maximum utility model to obtain a subsidy result;
and obtaining the energy sharing excitation subsidy of each park integrated energy system in the multi-park integrated energy system based on the subsidy result.
In one embodiment, the pre-constructed maximum utility model comprises: and the maximum utility objective function and the constraint condition are used for exciting the comprehensive energy systems of all the parks to participate in energy sharing.
In one embodiment, the maximum utility objective function is calculated as follows:
Figure BDA0003613528410000101
in the above formula, f is the maximum utility objective function value, Δ IESP For the economic benefit of the energy sharing of the comprehensive energy service provider,
Figure BDA0003613528410000102
for the economic benefit of the i-th park comprehensive energy system participating in energy sharing, omega i And sigma is the return coefficient of the ith park integrated energy system.
In one embodiment, the economic benefit of the multi-campus energy system energy sharing system participating in energy sharing is calculated as follows:
Figure BDA0003613528410000103
the calculation formula of the economic benefit of the ith park integrated energy system participating in energy sharing is as follows:
Figure BDA0003613528410000104
the calculation formula of the market contribution degree of the ith park integrated energy system is as follows:
Figure BDA0003613528410000105
in the above formula, the first and second carbon atoms are,
Figure BDA0003613528410000106
net operating cost, ζ, for the integrated energy facilitator after participation in energy sharing i Is subsidized for energy sharing and excitation of the integrated energy system of the ith park,
Figure BDA0003613528410000107
in order to participate in the operation cost of the ith campus comprehensive energy system after energy sharing,
Figure BDA0003613528410000108
and the shared value of the ith campus comprehensive energy system in the time period T, wherein T is the total time period.
The calculation formula of the shared value of the ith park integrated energy system in the t period is as follows:
Figure BDA0003613528410000111
in the above formula, the first and second carbon atoms are,
Figure BDA0003613528410000112
and
Figure BDA0003613528410000113
the local marginal prices for electrical, thermal, and cold loads, respectively.
Further, the process of solving the pre-constructed maximum utility model may be:
and (3) carrying out logarithm solving on two sides of a maximized utility objective function equation to obtain:
Figure BDA0003613528410000114
it is a concave function for ζ i Solving a first-order KKT condition:
Figure BDA0003613528410000115
summing the above N equations can be solved:
Figure BDA0003613528410000116
wherein the content of the first and second substances,
Figure BDA0003613528410000117
in the technical scheme provided by the invention, the constraint condition established for the energy interaction optimization of the multi-park comprehensive energy system can also be expressed as the change of the quantity and form of energy consumed by a user in unit time by the fundamental principle of comprehensive demand response, and mainly comprises the following three loads: load shedding, transferable loads and alternative loads are possible. The multipotent relationship of the integrated demand response is represented as:
Figure BDA0003613528410000118
wherein, P IDR,t 、H IDR,t 、G IDR,t And C IDR,t Respectively representing the power of electric, heat, gas and cold loads participating in comprehensive demand response in the time period t; p re,t 、P tr,t And P rp,t The reducible load, the transferable load and the alternative load of the electric energy in the t period are respectively expressed, and the heat load, the air load and the cold load are treated in the same way.
According to the actual needs of the system, the load adjustment amount allowed by three demand response modes is restricted, namely an adjustable load proportion range is given:
Figure BDA0003613528410000119
wherein delta re 、δ tr And delta rp The proportionality coefficients respectively representing reducible load, transferable load and substitutive load on the demand side are determined by specific load characteristics; x t Load level at time t before implementing the integrated demand response.
Example 2
In a second aspect, an energy interaction optimization apparatus for a multi-park integrated energy system is provided, as shown in fig. 4, the energy interaction optimization apparatus for a multi-park integrated energy system includes:
the acquisition module is used for acquiring the load of equipment of each park integrated energy system in the multi-park integrated energy system energy sharing system in the energy interaction process;
the solving module is used for substituting the load into a pre-constructed multi-park comprehensive energy system energy interaction optimization model and solving to obtain an optimization result;
and the analysis module is used for obtaining an energy interaction optimization scheme of the multi-park comprehensive energy system based on the optimization result.
Preferably, the load amount includes at least one of: electrical load, gas load, thermal load, cold load.
Preferably, the energy interaction optimization scheme of the multi-park integrated energy system comprises at least one of the following: EB power consumption power of the park comprehensive energy system in a time period t, EC equipment power consumption power of the park comprehensive energy system in the time period t, electric load power participating in comprehensive demand response, electric energy storage discharge power of the park comprehensive energy system in the time period t, electric energy sharing quantity of the park comprehensive energy system in the time period t, AC equipment heat consumption power of the park comprehensive energy system in the time period t, heat load power participating in comprehensive demand response, heat energy sharing quantity of the park comprehensive energy system in the time period t, CCHP gas consumption power of the park comprehensive energy system in the time period t, gas load power participating in comprehensive demand response and cold energy sharing quantity of the park comprehensive energy system in the time period t.
Further, the pre-constructed multi-park integrated energy system energy interaction optimization model comprises: and (4) constructing an objective function and constraint conditions for the energy interaction optimization of the multi-park comprehensive energy system.
Further, the objective function is calculated as follows:
Figure BDA0003613528410000121
in the above formula, C is the total cost of energy interaction optimization of the park comprehensive energy system,
Figure BDA0003613528410000122
in order to integrate the operating costs of the energy service provider,
Figure BDA0003613528410000123
the operation cost of the ith campus comprehensive energy system is shown, and I is the total number of the campus comprehensive energy systems.
Further, the mathematical model of the constraint is as follows:
Figure BDA0003613528410000124
Figure BDA0003613528410000125
Figure BDA0003613528410000126
Figure BDA0003613528410000127
Figure BDA0003613528410000128
Figure BDA0003613528410000131
in the above formula, the first and second carbon atoms are,
Figure BDA0003613528410000132
for the net demand of electrical energy for the ith campus complex energy system during time period t,
Figure BDA0003613528410000133
the electric energy balance of the ith park comprehensive energy system in the time period t is P e,i,t Electric load, P, for the ith park energy system during the period t PV,i,t The photovoltaic output level P of the ith park integrated energy system in the period of t WT,i,t For the wind power output level of the ith park comprehensive energy system in the time period t,
Figure BDA0003613528410000134
the CCHP electricity generation power of the ith park integrated energy system in the time period t,
Figure BDA0003613528410000135
the power consumption of the EB of the ith park comprehensive energy system in the time period t,
Figure BDA0003613528410000136
the power consumption of EC equipment in the ith park comprehensive energy system in the period of t is P IDR,t For participating in a combined demand response for a period tThe power of the electric load is controlled,
Figure BDA0003613528410000137
the electric energy storage and discharge power of the ith park comprehensive energy system in the time period t,
Figure BDA0003613528410000138
the electric energy sharing quantity P of the i-th park comprehensive energy system in the time period t h,i,t For the thermal load of the ith campus energy system during time t,
Figure BDA0003613528410000139
the CCHP heat production power of the ith campus integrated energy system in the period t,
Figure BDA00036135284100001310
for the EB heat generation power of the ith park comprehensive energy system in the period of t,
Figure BDA00036135284100001311
for the heat consumption power of the AC equipment of the ith park comprehensive energy system in the time period t, H IDR,t The thermal load power participating in the integrated demand response for the period t,
Figure BDA00036135284100001312
the heat energy sharing amount of the ith park integrated energy system in the time period t,
Figure BDA00036135284100001313
for the net demand of gas energy of the ith park energy system at time t,
Figure BDA00036135284100001314
for the CCHP gas consumption power, P, of the ith park integrated energy system in the period of t g,i,t Gas load of the i-th park integrated energy system at time t, G IDR,t For the time period t the gas load power participating in the integrated demand response,
Figure BDA00036135284100001315
is the ith gardenThe district integrated energy system has the refrigeration power of the EC equipment in the t period,
Figure BDA00036135284100001316
cooling power of AC equipment in time t period for ith park comprehensive energy system c,i,t For the cooling load of the ith campus energy system during the time period t,
Figure BDA00036135284100001317
for the cold energy sharing amount of the ith park integrated energy system in the time period t,
Figure BDA00036135284100001318
is the energy sharing upper limit of m types of energy sources in the unit time interval of the comprehensive energy system of the ith park,
Figure BDA00036135284100001319
the shared energy of the m types of energy sources in the unit time interval of the integrated energy system of the ith park is obtained.
Further, after obtaining the energy interaction optimization scheme of the multi-park integrated energy system based on the optimization result, the method includes:
solving a pre-constructed maximum utility model to obtain a subsidy result;
and obtaining the energy sharing excitation subsidy of each park integrated energy system in the multi-park integrated energy system based on the subsidy result.
Further, the pre-constructed maximum utility model comprises: and the maximum utility objective function and the constraint condition are used for exciting the comprehensive energy systems of all the parks to participate in energy sharing.
Further, the maximum utility objective function is calculated as follows:
Figure BDA0003613528410000141
in the above formula, f is the maximum utility objective function value, Δ IESP Economic benefits for energy sharing for integrated energy facilitators,
Figure BDA0003613528410000142
For the economic benefit of the i-th park comprehensive energy system participating in energy sharing, omega i And sigma is the return coefficient of the ith park integrated energy system.
Further, the calculation formula of the economic benefit of the multi-park integrated energy system energy sharing system participating in energy sharing is as follows:
Figure BDA0003613528410000143
the calculation formula of the economic benefit of the ith park integrated energy system participating in energy sharing is as follows:
Figure BDA0003613528410000144
the calculation formula of the market contribution degree of the ith park integrated energy system is as follows:
Figure BDA0003613528410000145
in the above-mentioned formula, the compound has the following structure,
Figure BDA0003613528410000146
net operating cost, ζ, for the integrated energy facilitator after participation in energy sharing i To be subsidized for energy sharing and excitation of the ith park comprehensive energy system,
Figure BDA0003613528410000147
in order to participate in the operation cost of the ith campus comprehensive energy system after energy sharing,
Figure BDA0003613528410000148
and the shared value of the ith park integrated energy system in the time period T is shown, and T is the total time period.
Further, the calculation formula of the shared value of the ith park integrated energy system in the time period t is as follows:
Figure BDA0003613528410000149
in the above formula, the first and second carbon atoms are,
Figure BDA00036135284100001410
and
Figure BDA00036135284100001411
the local marginal prices for electrical, thermal, and cold loads, respectively.
Example 3
Based on the same inventive concept, the present invention also provides a computer apparatus comprising a processor and a memory, the memory being configured to store a computer program comprising program instructions, the processor being configured to execute the program instructions stored by the computer storage medium. The Processor may be a Central Processing Unit (CPU), or may be other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf Programmable Gate Array (FPGA) or other Programmable logic device, a discrete Gate or transistor logic device, a discrete hardware component, etc., which is a computing core and a control core of the terminal, and is specifically adapted to implement one or more instructions, and specifically adapted to load and execute one or more instructions in a computer storage medium so as to implement a corresponding method flow or a corresponding function, so as to implement the steps of the energy interaction optimization method of the multiple park Integrated energy system in the foregoing embodiments.
Example 4
Based on the same inventive concept, the present invention further provides a storage medium, in particular, a computer-readable storage medium (Memory), which is a Memory device in a computer device and is used for storing programs and data. It is understood that the computer readable storage medium herein can include both built-in storage media in the computer device and, of course, extended storage media supported by the computer device. The computer-readable storage medium provides a storage space storing an operating system of the terminal. Also, one or more instructions, which may be one or more computer programs (including program code), are stored in the memory space and are adapted to be loaded and executed by the processor. It should be noted that the computer-readable storage medium may be a high-speed RAM memory, or may be a non-volatile memory (non-volatile memory), such as at least one disk memory. One or more instructions stored in the computer-readable storage medium may be loaded and executed by the processor to implement the steps of the energy interaction optimization method for the multi-campus integrated energy system according to the above embodiments.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting the same, and although the present invention is described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.

Claims (19)

1. A multi-park integrated energy system energy interaction optimization method is characterized by comprising the following steps:
acquiring the load capacity of equipment of each park comprehensive energy system in the multi-park comprehensive energy system energy sharing system in the energy interaction process;
substituting the load amount into a pre-constructed multi-park comprehensive energy system energy interaction optimization model and solving to obtain an optimization result;
and obtaining an energy interaction optimization scheme of the multi-park comprehensive energy system based on the optimization result.
2. The method of claim 1, wherein the load amount comprises at least one of: electrical load, gas load, thermal load, cold load.
3. The method of claim 1, wherein the multi-campus energy integrated system energy interaction optimization scenario comprises at least one of: EB power consumption power of the park comprehensive energy system in a time period t, EC equipment power consumption power of the park comprehensive energy system in the time period t, electric load power participating in comprehensive demand response, electric energy storage discharge power of the park comprehensive energy system in the time period t, electric energy sharing quantity of the park comprehensive energy system in the time period t, AC equipment heat consumption power of the park comprehensive energy system in the time period t, heat load power participating in comprehensive demand response, heat energy sharing quantity of the park comprehensive energy system in the time period t, CCHP gas consumption power of the park comprehensive energy system in the time period t, gas load power participating in comprehensive demand response and cold energy sharing quantity of the park comprehensive energy system in the time period t.
4. The method of claim 3, wherein the pre-constructed multi-campus energy integrated system energy interaction optimization model comprises: and an objective function and a constraint condition are established for the energy interaction optimization of the multi-park comprehensive energy system.
5. The method of claim 4, wherein the objective function is calculated as follows:
Figure FDA0003613528400000011
in the above formula, C is the total cost of energy interaction optimization of the park comprehensive energy system,
Figure FDA0003613528400000012
in order to integrate the operating costs of the energy service provider,
Figure FDA0003613528400000013
the operation cost of the ith campus comprehensive energy system is shown, and I is the total number of the campus comprehensive energy systems.
6. The method of claim 5, wherein the mathematical model of the constraints is as follows:
Figure FDA0003613528400000014
Figure FDA0003613528400000015
Figure FDA0003613528400000016
Figure FDA0003613528400000017
Figure FDA0003613528400000018
Figure FDA0003613528400000021
in the above formula, the first and second carbon atoms are,
Figure FDA0003613528400000022
for the net demand of electric energy of the integrated energy system of the ith park at the time period t,
Figure FDA0003613528400000023
electric energy net surplus of the ith park integrated energy system in time period tAmount, P e,i,t Electric load of the i-th park integrated energy system in the period of t, P PV,i,t The photovoltaic output level P of the ith park integrated energy system in the period of t WT,i,t For the wind power output level of the ith park comprehensive energy system in the time period t,
Figure FDA0003613528400000024
the CCHP electricity generation power of the ith park integrated energy system in the time period t,
Figure FDA0003613528400000025
the power consumption of the EB of the ith park comprehensive energy system in the time period t,
Figure FDA0003613528400000026
the power consumption of EC equipment in the ith park comprehensive energy system in the period of t is P IDR,t For the time period t the electrical load power participating in the integrated demand response,
Figure FDA0003613528400000027
the electric energy storage and discharge power of the ith park comprehensive energy system in the time period t,
Figure FDA0003613528400000028
the electric energy sharing quantity P of the i-th park comprehensive energy system in the time period t h,i,t For the heat load of the ith campus integrated energy system during the period t,
Figure FDA0003613528400000029
for the CCHP heat-generating power of the ith park integrated energy system in the period t,
Figure FDA00036135284000000210
for the EB heat generation power of the ith park comprehensive energy system in the period of t,
Figure FDA00036135284000000211
for the i-th park comprehensive energyAC appliance heat power consumption of source system during time period t, H IDR,t The thermal load power participating in the integrated demand response for the period t,
Figure FDA00036135284000000212
for the heat energy sharing amount of the ith park integrated energy system in the t period,
Figure FDA00036135284000000213
for the net demand of gas energy of the ith park energy system at time t,
Figure FDA00036135284000000214
for the CCHP gas consumption power, P, of the ith park integrated energy system in the period of t g,i,t Gas load of the i-th park integrated energy system at time t, G IDR,t For the time period t the gas load power participating in the integrated demand response,
Figure FDA00036135284000000215
for the cooling power of the EC equipment of the ith park integrated energy system in the t period,
Figure FDA00036135284000000216
cooling power of AC equipment in time t period for ith park comprehensive energy system c,i,t For the cooling load of the ith campus energy system during the time period t,
Figure FDA00036135284000000217
for the cooling energy sharing amount of the ith campus comprehensive energy system in the time period t,
Figure FDA00036135284000000218
is the energy sharing upper limit of m types of energy sources in the unit time interval of the comprehensive energy system of the ith park,
Figure FDA00036135284000000219
for the m types of energy sources in the unit time interval of the ith park integrated energy systemTo share energy.
7. The method of claim 6, wherein the obtaining the multi-campus energy system energy interaction optimization scheme based on the optimization result comprises:
solving a pre-constructed maximum utility model to obtain a subsidy result;
and obtaining the energy sharing excitation subsidy of each park integrated energy system in the multi-park integrated energy system energy sharing system based on the subsidy result.
8. The method of claim 7, wherein the pre-constructed maximized utility model comprises: and (4) maximizing a utility objective function and constraint conditions for exciting the comprehensive energy system of each park to participate in energy sharing.
9. The method of claim 8, wherein the maximizing utility objective function is calculated as follows:
Figure FDA0003613528400000031
in the above formula, f is the maximum utility objective function value, Δ IESP Economic benefits, Δ, for integrated energy providers to participate in energy sharing i PIES For the economic benefit of the i-th park comprehensive energy system participating in energy sharing, omega i And sigma is the return coefficient of the ith park integrated energy system.
10. The method of claim 9, wherein the economic benefit of the multi-campus energy system energy sharing system participating in energy sharing is calculated as follows:
Figure FDA0003613528400000032
the calculation formula of the economic benefit of the ith park integrated energy system participating in energy sharing is as follows:
Figure FDA0003613528400000033
the calculation formula of the market contribution degree of the ith park integrated energy system is as follows:
Figure FDA0003613528400000034
in the above-mentioned formula, the compound has the following structure,
Figure FDA0003613528400000035
net operating cost, ζ, for the integrated energy facilitator after participation in energy sharing i To be subsidized for energy sharing and excitation of the ith park comprehensive energy system,
Figure FDA0003613528400000036
in order to participate in the operation cost of the ith campus comprehensive energy system after energy sharing,
Figure FDA0003613528400000037
and the shared value of the ith park integrated energy system in the time period T is shown, and T is the total time period.
11. The method of claim 10, wherein the shared value of the ith campus energy system during time t is calculated as follows:
Figure FDA0003613528400000038
in the above formula, the first and second carbon atoms are,
Figure FDA0003613528400000039
and
Figure FDA00036135284000000310
the local marginal price of the electrical load, the thermal load, the cold load, respectively.
12. An energy interaction optimization device for a multi-park integrated energy system, the device comprising:
the acquisition module is used for acquiring the load capacity of equipment of each park integrated energy system in the multi-park integrated energy system energy sharing system in the energy interaction process;
the solving module is used for substituting the load quantity into a pre-constructed multi-park comprehensive energy system energy interaction optimization model and solving to obtain an optimization result;
and the analysis module is used for obtaining an energy interaction optimization scheme of the multi-park comprehensive energy system based on the optimization result.
13. The apparatus of claim 12, wherein the amount of load comprises at least one of: electrical load, gas load, thermal load, cold load.
14. The apparatus of claim 12, wherein the multi-campus energy integrated system energy interaction optimization scheme comprises at least one of: EB power consumption power of the park comprehensive energy system in a time period t, EC equipment power consumption power of the park comprehensive energy system in the time period t, electric load power participating in comprehensive demand response, electric energy storage discharge power of the park comprehensive energy system in the time period t, electric energy sharing quantity of the park comprehensive energy system in the time period t, AC equipment heat consumption power of the park comprehensive energy system in the time period t, heat load power participating in comprehensive demand response, heat energy sharing quantity of the park comprehensive energy system in the time period t, CCHP gas consumption power of the park comprehensive energy system in the time period t, gas load power participating in comprehensive demand response and cold energy sharing quantity of the park comprehensive energy system in the time period t.
15. The apparatus of claim 14, wherein the pre-constructed multi-campus energy integrated system energy interaction optimization model comprises: and an objective function and a constraint condition are established for the energy interaction optimization of the multi-park comprehensive energy system.
16. The apparatus of claim 15, wherein the objective function is calculated as follows:
Figure FDA0003613528400000041
in the above formula, C is the total cost of energy interaction optimization of the park comprehensive energy system,
Figure FDA0003613528400000042
in order to integrate the operating costs of the energy service provider,
Figure FDA0003613528400000043
the operation cost of the ith park comprehensive energy system is shown, and I is the total number of the park comprehensive energy systems.
17. The apparatus of claim 16, wherein the mathematical model of the constraints is as follows:
Figure FDA0003613528400000044
Figure FDA0003613528400000045
Figure FDA0003613528400000046
Figure FDA0003613528400000047
Figure FDA0003613528400000048
Figure FDA0003613528400000049
in the above formula, the first and second carbon atoms are,
Figure FDA00036135284000000410
for the net demand of electrical energy for the ith campus complex energy system during time period t,
Figure FDA00036135284000000411
the electric energy balance of the ith park comprehensive energy system in the time period t, P e,i,t Electric load of the i-th park integrated energy system in the period of t, P PV,i,t The photovoltaic output level P of the ith park integrated energy system in the period of t WT,i,t For the wind power output level of the ith campus comprehensive energy system in the time period t,
Figure FDA0003613528400000051
the CCHP power generation power of the ith garden integrated energy system in the t period,
Figure FDA0003613528400000052
the power consumption of the EB of the ith park comprehensive energy system in the time period t,
Figure FDA0003613528400000053
the power consumption of EC equipment in the ith park comprehensive energy system in the period of t is P IDR,t For the time period t the electrical load power participating in the integrated demand response,
Figure FDA0003613528400000054
the electric energy storage and discharge power of the ith park comprehensive energy system in the time period t,
Figure FDA0003613528400000055
the electric energy sharing quantity P of the i-th park comprehensive energy system in the time period t h,i,t For the thermal load of the ith campus energy system during time t,
Figure FDA0003613528400000056
for the CCHP heat-generating power of the ith park integrated energy system in the period t,
Figure FDA0003613528400000057
for the EB heat generation power of the ith park comprehensive energy system in the period of t,
Figure FDA0003613528400000058
AC equipment heat power consumption H for ith district integrated energy system in t time period IDR,t The thermal load power participating in the integrated demand response for the period t,
Figure FDA0003613528400000059
the heat energy sharing amount of the ith park integrated energy system in the time period t,
Figure FDA00036135284000000510
for the net demand of gas energy of the ith park energy system at time t,
Figure FDA00036135284000000511
for the CCHP gas consumption power, P, of the ith park integrated energy system in the period of t g,i,t Air load of the i-th park integrated energy system at time t, G IDR,t For the time period t the gas load power participating in the integrated demand response,
Figure FDA00036135284000000512
is a firstThe cooling power of the EC equipment of the i park integrated energy systems in the time period t,
Figure FDA00036135284000000513
cooling power of AC equipment in time t period for ith park comprehensive energy system c,i,t For the cooling load of the ith campus energy system during the period t,
Figure FDA00036135284000000514
for the cold energy sharing amount of the ith park integrated energy system in the time period t,
Figure FDA00036135284000000515
is the energy sharing upper limit of m types of energy sources in the unit time interval of the comprehensive energy system of the ith park,
Figure FDA00036135284000000516
the shared energy of the m types of energy sources in the unit time interval of the comprehensive energy system of the ith park is obtained.
18. A computer device, comprising: one or more processors;
the processor to store one or more programs;
the one or more programs, when executed by the one or more processors, implement the method for energy interactive optimization of a multi-campus energy complex according to any one of claims 1 to 11.
19. A computer-readable storage medium having stored thereon a computer program which, when executed, implements a method for energy interactive optimization of a multi-campus energy complex according to any one of claims 1 to 11.
CN202210440078.XA 2022-04-25 2022-04-25 Energy interaction optimization method and device for multi-park comprehensive energy system Pending CN115034551A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210440078.XA CN115034551A (en) 2022-04-25 2022-04-25 Energy interaction optimization method and device for multi-park comprehensive energy system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210440078.XA CN115034551A (en) 2022-04-25 2022-04-25 Energy interaction optimization method and device for multi-park comprehensive energy system

Publications (1)

Publication Number Publication Date
CN115034551A true CN115034551A (en) 2022-09-09

Family

ID=83119317

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210440078.XA Pending CN115034551A (en) 2022-04-25 2022-04-25 Energy interaction optimization method and device for multi-park comprehensive energy system

Country Status (1)

Country Link
CN (1) CN115034551A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117273810A (en) * 2023-11-03 2023-12-22 连云港智源电力设计有限公司 Comprehensive energy sharing scheduling method and system with excitation compatibility

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117273810A (en) * 2023-11-03 2023-12-22 连云港智源电力设计有限公司 Comprehensive energy sharing scheduling method and system with excitation compatibility
CN117273810B (en) * 2023-11-03 2024-04-05 连云港智源电力设计有限公司 Comprehensive energy sharing scheduling method and system with excitation compatibility

Similar Documents

Publication Publication Date Title
Wei et al. A bi-level scheduling model for virtual power plants with aggregated thermostatically controlled loads and renewable energy
Moghaddam et al. A comprehensive model for self-scheduling an energy hub to supply cooling, heating and electrical demands of a building
Liu et al. Hybrid energy sharing for smart building cluster with CHP system and PV prosumers: A coalitional game approach
Iwafune et al. Cooperative home energy management using batteries for a photovoltaic system considering the diversity of households
CN112180741B (en) Comprehensive demand response method and system with participation of load aggregators
CN111008739B (en) Optimal regulation and control and income distribution method and system for cogeneration virtual power plant
Lim et al. Optimal allocation of energy storage and solar photovoltaic systems with residential demand scheduling
Harb et al. Decentralized scheduling strategy of heating systems for balancing the residual load
Belli et al. A unified model for the optimal management of electrical and thermal equipment of a prosumer in a DR environment
Li et al. Electricity scheduling optimisation based on energy cloud for residential microgrids
CN105528658B (en) Resident user cooperation game power utilization arrangement method considering household distributed power supply
CN113746089A (en) Multi-user-oriented multi-time-scale power package and family energy optimization method
Diekerhof et al. Distributed optimization for scheduling electrical demand in complex city districts
Fang et al. Energy management strategy based on energy storage equalization technology and transferable load
Wang et al. Event-triggered online energy flow control strategy for regional integrated energy system using Lyapunov optimization
Cai et al. Optimal scheduling of a hybrid AC/DC multi-energy microgrid considering uncertainties and Stackelberg game-based integrated demand response
Arnold et al. Investigating renewable infeed in residential areas applying model predictive control
CN115034551A (en) Energy interaction optimization method and device for multi-park comprehensive energy system
Chen et al. Multi-objective robust optimal bidding strategy for a data center operator based on bi-level optimization
CN113487188A (en) Comprehensive energy system optimal scheduling method considering electric and gas joint price guide mechanism
Siepermann et al. The potential of shifting residential energy consumption for the energy transition
CN115730747A (en) Multi-subject benefit distribution method of comprehensive energy system and application thereof
CN115659585A (en) Micro-energy network low-carbon cooperative scheduling method and device considering demand response, memory and equipment
CN115481858A (en) Electric-gas system optimal power flow dispersion calculation method and system considering demand response
Lei et al. Two‐stage game framework for energy management in islanded multi‐microgrid system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication