CN113592365B - Energy optimization scheduling method and system considering carbon emission and green electricity consumption - Google Patents

Energy optimization scheduling method and system considering carbon emission and green electricity consumption Download PDF

Info

Publication number
CN113592365B
CN113592365B CN202111005567.4A CN202111005567A CN113592365B CN 113592365 B CN113592365 B CN 113592365B CN 202111005567 A CN202111005567 A CN 202111005567A CN 113592365 B CN113592365 B CN 113592365B
Authority
CN
China
Prior art keywords
energy
service provider
park
ith
buy
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111005567.4A
Other languages
Chinese (zh)
Other versions
CN113592365A (en
Inventor
李锋
陆仕荣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University
Chongqing Institute of Green and Intelligent Technology of CAS
Original Assignee
Chongqing University
Chongqing Institute of Green and Intelligent Technology of CAS
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 Chongqing University, Chongqing Institute of Green and Intelligent Technology of CAS filed Critical Chongqing University
Priority to CN202111005567.4A priority Critical patent/CN113592365B/en
Publication of CN113592365A publication Critical patent/CN113592365A/en
Application granted granted Critical
Publication of CN113592365B publication Critical patent/CN113592365B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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
    • 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
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S50/00Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
    • Y04S50/16Energy services, e.g. dispersed generation or demand or load or energy savings aggregation

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Game Theory and Decision Science (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Health & Medical Sciences (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Water Supply & Treatment (AREA)
  • Primary Health Care (AREA)
  • Public Health (AREA)
  • General Health & Medical Sciences (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides an energy optimization scheduling method and system considering carbon emission and green electricity consumption, which comprises the steps of constructing an output model of energy supply equipment according to the running characteristics of the energy supply equipment in a park and a carbon emission right; configuring carbon emission rights and green electricity consumption indexes for different user types, and creating energy demand constraints of output models of energy supply equipment according to different user types; according to the energy demand constraint and the energy prices and powers of the current park energy suppliers and other park energy suppliers, an outer energy supplier optimized scheduling model is constructed by taking the optimal income of the energy suppliers as a target, and the outer energy supplier optimized scheduling model is solved through an optimization algorithm to obtain an energy optimal scheduling scheme; according to the invention, the carbon emission right and the green electricity consumption are fully considered, and the accuracy and the applicability of energy optimization scheduling are improved.

Description

Energy optimization scheduling method and system considering carbon emission and green electricity consumption
Technical Field
The invention relates to the field of new energy application, in particular to an energy optimization scheduling method and system considering carbon emission and green electricity consumption.
Background
The traditional single energy system is limited to single energy forms such as electricity, gas and heat, cannot give full play to advantage complementation among energy sources, and cannot effectively coordinate comprehensive utilization of multiple energy sources, renewable energy consumption and the like. With the rapid development of economy, energy and environmental problems are increasingly prominent, and how to realize the clean and efficient utilization of multiple energy sources becomes a key point of research in recent years, and the concepts of energy internet and comprehensive energy systems are generated at the same time.
Many experts and scholars have conducted a great deal of research on electricity-heat-gas comprehensive energy systems in certain parks (such as residential parks, industrial parks, etc.) and parks, and have achieved remarkable results. But the capacity of a single campus is limited and the types of users in a single campus are relatively single. The green electricity consumption index and the carbon emission rate have increasingly more influence on the energy demand of different users in the garden, and the status of the user side on the optimal scheduling of the comprehensive energy system in the garden or a plurality of parks in the area is also increasingly important. In the existing optimal scheduling scheme of the comprehensive energy of the park, the optimal configuration of a regional multi-park comprehensive energy system considering the green electricity requirements of different users does not exist, and the important role played by the users in the system optimal scheduling of the energy Internet and the comprehensive energy system can not be effectively reflected.
Disclosure of Invention
In view of the problems in the prior art, the invention provides an energy optimization scheduling method and system considering carbon emission and green electricity consumption, and mainly solves the problems that the green electricity consumption of users cannot be accurately and effectively reflected by the traditional energy scheduling, and the applicability is poor.
In order to achieve the above and other objects, the present invention adopts the following technical solutions.
An energy optimization scheduling method considering carbon emission and green electricity consumption comprises the following steps:
constructing a power output model of the energy supply equipment according to the operating characteristics of the energy supply equipment in the park and the carbon emission right; configuring carbon emission rights and green electricity consumption indexes for different user types, and creating energy demand constraints of output models of energy supply equipment according to different user types; wherein, the output model of the energy supply equipment comprises: the system comprises a gas turbine output model, a boiler equipment output model, a new energy power generation equipment output model, an electric heating equipment output model, an electric refrigerating equipment output model and an electric hydrogen production equipment output model; the user types include: large industrial users, general industrial and commercial users, and residential users;
and constructing an outer energy service provider optimized dispatching model by taking the optimal income of the energy provider as a target according to the energy demand constraint and the energy prices and powers of the current park energy provider and other park energy providers, and solving the outer energy service provider optimized dispatching model through an optimization algorithm to obtain an energy optimal dispatching scheme.
Optionally, the gas turbine output model is:
V GT ≤δ(C lim_GT +C lim_GT_buy )
C GT =aC lim_GT_buy
in the formula: v GT Natural gas consumption for gas turbines; delta is the carbon emission coefficient, C lim_GT To allocate carbon emission limits; c GT Represents the extra carbon cost for the energy supplier to generate electricity using the gas turbine; c lim_GT_buy Represents the extra carbon emission rights purchased by the energy supplier to ensure normal production; a is a carbon unit price;
the boiler equipment comprises a gas boiler, and the output model of the boiler equipment is as follows:
V GB ≤δ(C lim_GB +C lim_GB_buy )
C GB =aC lim_GB_buy
wherein, V GB Consuming power for gas boiler natural gas; c lim_GB A carbon emission allowance assigned to the gas boiler plant; c GB Represents the extra carbon cost for the supplier to supply heat with the gas boiler; c lim_GB_buy Represents the extra carbon emission rights purchased by the energy supplier to ensure normal production heating;
the new energy power generation equipment comprises solar power generation equipment and wind power generation equipment, and the output model of the new energy power generation equipment is as follows:
C PV_WT =aδ PV_WT (P WT +P PV )Δt
wherein, C PV_WT Income from new energy generation, delta PV_WT A carbon emission factor; (ii) a P WT Outputting power for wind power generation; p PV Outputting power for photovoltaic power generation; Δ t is the time span;
the electric heating equipment comprises electric boiler equipment, and the output model of the electric heating equipment is as follows:
C EB =aδ EB P EB
wherein, P EB Consuming electrical energy for the electrical heating device; c EB Additional carbon cost, δ, of using electrical heating equipment for non-green electricity EB Energy consumption equivalent carbon emission factors for the electric heating equipment;
the output model of the electric refrigeration equipment is as follows:
C HP =aδ HP P HP
wherein, P HP Consuming electrical energy for the electrical refrigeration equipment; c HP The extra carbon cost, δ, of using electric refrigeration equipment for non-green electricity HP Energy-consuming equivalent carbon emission factor for electric refrigeration equipment;
The output model of the electrical hydrogen production equipment is as follows:
C SOEC =δ SOEC P SOEC
wherein, P SOEC Electrical energy consumed for the electrolytic cell; c SOEC The additional carbon cost of using electrical hydrogen production equipment for non-green electricity; delta SOEC The unit of the electric hydrogen production equipment consumes energy and carbon cost.
Optionally, configuring carbon emission rights and green electricity consumption indexes for different user types, and creating energy demand constraints of output models of each energy supply device according to different user types, including:
according to the carbon emission right and the renewable energy power consumption responsibility weight index of the large industrial user, creating an energy demand constraint of the large industrial user, wherein the energy demand constraint is expressed as:
C em,user1 =δ 1 P user1
P user1 =P user1_PV_WT_buy +P user1_PV_WT +P user1_buy
0≤C em,user1 ≤C em,user1_lim +C em,user1_buy
0≤ε 1 P user1 ≤P user1_PV_WT_buy +P user1_PV_WT
C user1 =aC em,user1_buy +bP user1_PV_WT_buy +cP user1_buy +d
wherein the carbon emissions, C, are estimated from the total electricity consumption of large industrial users em,user1 Represents the carbon emissions of the large industrial user; delta. for the preparation of a coating 1 Represents a carbon emission conversion factor; p user1 Representing the total electric energy consumption of the large industrial user; p user1_PV_WT_buy Representing green electric energy purchased by the large industrial user from an integrated energy service provider;
Figure GDA0003549211490000031
green electric energy representing self distributed photovoltaic and wind power generation of the large industrial users; p user1_buy Represents conventional electrical energy purchased by the large industrial user; c em,user1_lim Representing carbon emission rights assigned to the large industrial user; c em,user1_buy Representing the carbon emission rights purchased by the large industrial user from the comprehensive energy service provider, and when the value is negative, representing that the surplus carbon emission rights are traded with the energy service provider; epsilon 1 Representing a renewable energy power consumption liability weight for the large industrial user; c user1 Represents the total cost of production energy purchase for the large industrial user; aC em,user1_buy Represents a cost of purchasing carbon emissions rights; bP user1_PV_WT_buy Represents the cost of purchasing green electricity; cP (personal computer) user1_buy Represents the cost of purchasing traditional electrical energy; d represents other costs for the large industrial user to purchase energy costs;
the energy demand constraints of the general industrial and commercial users are as follows:
C em,user2 =δ 2 P user2
P user2 =P user2_PV_WT_buy +P user2_PV_WT
0≤C em,user2 ≤C em,user2_lim +C em,user2_buy
0≤ε 2 P user2 ≤P user2_PV_WT_buy +P user2_PV_WT
C user2 =aC em,user2_buy +bP user2_PV_WT_buy +c
wherein, C em,user2 Represents the general industrial and commercial user carbon emission; delta 2 Representing the carbon emission conversion coefficient of the general industrial and commercial users; p user2 Representing the total electric energy consumption of the general industrial and commercial users; p user2_PV_WT_buy The green electric energy which is purchased from the comprehensive energy service provider by the general industrial and commercial user is represented, the value is 0 or a negative value, and the green electric energy is not purchased from the comprehensive energy service provider or surplus green electric power is sold to the comprehensive energy service provider; p user2_PV_WT Green electric energy representing the self distributed photovoltaic and wind power generation of the general industrial and commercial users; c em,user2_lim Representing carbon emission rights assigned to the general industrial and commercial user; c em,user2_buy Represents the carbon purchased from the general industrial and commercial users to the comprehensive energy service providerThe emission right indicates that the surplus carbon emission right trades with the energy service provider when the value is negative; epsilon 2 A renewable energy power consumption responsibility weight representing the general industrial and commercial user; c user2 Representing the total cost of the general industrial and commercial user production energy purchase; aC em,user2_buy Represents a cost of purchasing carbon emissions rights; bP user2_PV_WT_buy Represents the cost of purchasing green electricity; c represents other costs for the general industry and commerce users to purchase energy;
the energy demand constraint of the residential user is as follows:
P user3 =P user3_PV_WT +P user3_PV_WT_buy
C user3 =aP user3_PV_WT_buy +b
not considering the carbon emission right and the green electricity consumption index of the residential user, wherein P user3 Representing the total power consumption of the residential user; p user3_PV_WT_buy The method comprises the steps that green electric energy purchased from a comprehensive energy service provider by a resident user is represented, if the value is positive, the distributed energy of the resident user cannot meet the self demand, the energy is required to be purchased from the energy service provider, if the value is zero, the energy is not purchased from the comprehensive energy service provider, and if the value is negative, the surplus green electric power of the resident user is sold to the comprehensive energy service provider; p user3_PV_WT Representing the green electric energy generated by the distributed photovoltaic and wind power of the resident users; c user3 Representing the total cost of energy purchase by the residential user; aP user3_PV_WT_buy Representing the cost of purchasing green electricity, and representing the surplus energy profit when the value is negative; b represents other energy cost of the residential user.
Optionally, according to the energy demand constraint and the energy prices and powers of the current park energy provider and other park energy providers, with the optimal income of the energy provider as a target, an outer energy provider optimized scheduling model is constructed, and the outer energy provider optimized scheduling model is solved through an optimization algorithm to obtain an energy optimal scheduling scheme, including:
setting comprehensive energy service providers of all parks in the area to perform external energy trade through information sharing, and constructing an outer-layer energy service provider optimized scheduling model; each regional park service provider carries out optimization regulation and control by taking the benefit maximization of the regional park service provider as a target, and the optimization target function is as follows:
Figure GDA0003549211490000051
wherein the content of the first and second substances,
Figure GDA0003549211490000052
Figure GDA0003549211490000053
Figure GDA0003549211490000054
Figure GDA0003549211490000055
Figure GDA0003549211490000056
Figure GDA0003549211490000057
wherein N is a regional park set, which indicates that N parks of the regional comprehensive energy network participate in energy trading, and N is the number of the parks; f is the total income of the park comprehensive energy service provider;
Figure GDA0003549211490000058
the electricity selling income of the comprehensive energy service provider of the ith park is obtained for the time period t;
Figure GDA0003549211490000059
for a time period t of the ith campusThe electricity selling income among the comprehensive energy service providers;
Figure GDA00035492114900000510
the comprehensive energy service provider of the ith park sells electricity to the power grid for the time t;
Figure GDA00035492114900000511
purchasing energy costs from the ith campus energy supplier for time period t;
Figure GDA00035492114900000512
the energy interaction cost of the ith park and other parks in the area is t;
Figure GDA00035492114900000513
the electricity purchasing cost from the power grid for the ith park comprehensive energy service provider in the time period t;
Figure GDA00035492114900000514
the unit electric energy income for the ith energy service provider and the nth park service provider,
Figure GDA00035492114900000515
The unit heat energy profit for the ith energy service provider and the nth park service provider,
Figure GDA00035492114900000516
The cost for the interaction of the energy service provider and the service providers of other parks in the area,
Figure GDA00035492114900000517
Heat energy interaction cost is provided for the ith energy service provider and the nth park service provider in the region; p in Electric energy interaction quantity H for ith energy service provider and nth regional park service provider in The heat energy interaction amount is provided for the ith energy service provider and the nth park service provider in the region;
Figure GDA0003549211490000061
the energy sale cost of the ith park energy service provider from the power grid unit,
Figure GDA0003549211490000062
Purchasing energy cost from a power grid unit for the ith park energy service provider in the time period t;
Figure GDA0003549211490000063
the electric energy interaction quantity between the energy service provider of the ith park and the power grid is t;
Figure GDA0003549211490000064
the unit electric energy cost purchased by the ith park energy service provider and the own park energy provider in the time period t,
Figure GDA0003549211490000065
The unit heat energy cost purchased by the ith park energy service provider and the energy supplier of the park at the time t;
Figure GDA0003549211490000066
the electric energy purchased by the energy service provider of the park for the ith time period t through the energy provider of the park,
Figure GDA0003549211490000067
The heat energy purchased by the energy service provider of the ith park in the time period t through the energy provider of the park;
Figure GDA0003549211490000068
the unit price of electricity sold to the user by the ith campus energy service provider for the time period t,
Figure GDA0003549211490000069
The unit heat energy price sold to the user by the ith park energy service provider for the time period t; Δ t is the time span.
Optionally, the setting of the energy external transaction by the integrated energy service provider of each park in the area through information sharing includes:
creating external trading power demand constraints and energy supply constraints for the campus and other camps in the campus, expressed as maximum sales or purchase revenue objective function for the campus and other camps
Figure GDA00035492114900000610
Wherein the content of the first and second substances,
Figure GDA00035492114900000611
Figure GDA00035492114900000612
Figure GDA00035492114900000613
Figure GDA00035492114900000614
Figure GDA00035492114900000615
Figure GDA00035492114900000616
wherein, P out ,H out Respectively carrying out external transaction on electric energy and heat energy, including transaction with other regional energy service providers and energy transaction with a network;
Figure GDA00035492114900000617
representing the thermal energy purchased by the ith campus energy service provider from the grid unit for time period t.
Optionally, the optimization algorithm comprises a particle swarm algorithm or an ant colony algorithm.
An energy optimization scheduling system considering carbon emission and green electricity consumption, comprising:
the model and constraint creation module is used for constructing a power output model of the energy supply equipment according to the operating characteristics of the energy supply equipment in the park and the carbon emission right; configuring carbon emission rights and green electricity consumption indexes for different user types, and creating energy demand constraints of output models of energy supply equipment according to different user types; wherein, the output model of the energy supply equipment comprises: the system comprises a gas turbine output model, a boiler equipment output model, a new energy power generation equipment output model, an electric heating equipment output model, an electric refrigerating equipment output model and an electric hydrogen production equipment output model; the user types include: large industrial users, general industrial and commercial users, and residential users;
and the optimization calculation module is used for constructing an outer energy service provider optimization scheduling model by taking the optimal income of the energy provider as a target according to the energy demand constraint and the energy prices and powers of the current park energy provider and other park energy providers, and solving the outer energy service provider optimization scheduling model through an optimization algorithm to obtain an energy optimal scheduling scheme.
As described above, the energy optimization scheduling method and system considering carbon emission and green power consumption according to the present invention have the following advantages.
The comprehensive energy system established by the invention considers the carbon emission index factor of the energy industry, adds the carbon transaction into the production energy of power generation enterprises or users, and perfects the traditional equivalent model;
users in the traditional model are divided into large industrial users, general industrial and commercial users and residential users, so that the carbon emission rights and green electricity consumption requirements of different users are fully reflected, and the important position of the users in the comprehensive energy system is highlighted;
through cooperation and transaction among energy service providers, energy interaction among multiple regional garden intervals is promoted, and the energy utilization efficiency of the comprehensive energy system of the garden is further improved.
Drawings
FIG. 1 is a topological diagram of an integrated energy system for a single campus in accordance with the present invention.
FIG. 2 is a topological diagram of the regional multi-park integrated energy system of the present invention.
FIG. 3 is a flow chart of the optimal scheduling solution of the regional multi-park integrated energy system of the present invention.
FIG. 4 is a flow chart of solving particle swarm optimization used in the present invention.
Detailed Description
The embodiments of the present invention are described below with reference to specific embodiments, and other advantages and effects of the present invention will be easily understood by those skilled in the art from the disclosure of the present specification. The invention is capable of other and different embodiments and of being practiced or of being carried out in various ways, and its several details are capable of modification in various respects, all without departing from the spirit and scope of the present invention. It is to be noted that the features in the following embodiments and examples may be combined with each other without conflict.
It should be noted that the drawings provided in the following embodiments are only for illustrating the basic idea of the present invention, and the components related to the present invention are only shown in the drawings rather than drawn according to the number, shape and size of the components in actual implementation, and the type, quantity and proportion of the components in actual implementation may be changed freely, and the layout of the components may be more complicated.
The regional multi-park comprehensive energy system optimization scheduling method considering carbon emission and green electricity requirements of different users adds carbon emission rights into an equivalent model, divides users into large industrial users, general industrial and commercial users and residential users, and constructs a regional multi-park comprehensive energy system optimization scheduling model by taking an energy service provider as a button and taking energy price and power provided by the service provider as links. The model comprises an inner-layer single-park comprehensive energy network system and an outer-layer regional multi-park energy service provider cooperation transaction model.
In the model solving process, firstly, different users propose energy demand according to self carbon emission rights, energy demand and current energy price; then, the optimal scheduling model of the inner layer of the park performs optimal scheduling on the energy output of the energy supplier, and simultaneously returns the output result to the energy service provider; and then the energy service provider integrally plans the energy output condition and the green electricity consumption condition of the energy provider in the park and the energy prices and power of the energy providers in other parks, carries out regional multi-park coordination solving through the outer layer optimization scheduling model, obtains a regional multi-park comprehensive energy system optimization scheduling scheme, and finally determines the energy prices of all parks.
The invention discloses an energy optimization scheduling method considering carbon emission and green electricity consumption, which comprises the following steps:
1) constructing a power output model of the energy supply equipment according to the operating characteristics of the energy supply equipment in the park and the carbon emission right; configuring carbon emission rights and green electricity consumption indexes for different user types, and creating energy demand constraints of output models of energy supply equipment according to different user types; wherein, the output model of the energy supply equipment comprises: the system comprises a gas turbine output model, a boiler equipment output model, a new energy power generation equipment output model, an electric heating equipment output model, an electric refrigerating equipment output model and an electric hydrogen production equipment output model; the user types include: large industrial users, general industrial and commercial users, and residential users;
energy supply equipment include gas turbine, gas boiler, electricity heating equipment, electric refrigeration plant, absorption chiller, electricity hydrogen manufacturing equipment, new forms of energy power generation facility etc. each power consumption equipment considers carbon emission right, and new forms of energy power generation facility considers carbon income, the user considers green electricity consumption index and carbon income according to the type difference, specific equivalent model is as follows:
(1.1) wherein the gas turbine generates high-temperature gas to expand and do work through the mixed combustion of compressed air and natural gas so as to generate electric energy, the high-temperature flue gas after combustion meets the cold and hot load requirements of users through an absorption refrigerator or a heat exchanger, and the equivalent mathematical expression of the gas turbine output model is as follows:
G GT,t =V GT,t
P GT,t =G GT,tGT
Figure GDA0003549211490000091
P GT,t,min ≤P GT,t ≤P GT,t,max
V GT ≤δ(C lim_GT +C lim_GT_buy )
C GT =aC lim_GT_buy
in the formula: g GT,t Consumption of input power, V, of natural gas for gas turbines GT,t Natural gas consumption for gas turbines; h gas Is the heat value of natural gas; p GT,t For the power, eta, of the gas turbine GT Efficiency of power supply to the gas turbine; q WHB,t Is the heat recovery power, eta of the waste heat boiler WHB The heat recovery efficiency of the waste heat boiler; p GT,t,max Is the upper limit, P, of the electrical output of the gas turbine GT,t,min A lower limit for the electrical output power of the gas turbine; is the carbon emission coefficient, C lim_GT To allocate carbon emission limits. C GT Represents the extra carbon cost for the supplier to generate electricity using the gas turbine; c lim_GT_buy Represents the extra carbon emission rights purchased by the supplier to ensure normal production;
(1.2) the refrigeration power of the absorption chiller is related to the heat energy discharged from the gas turbine, so that the carbon emission power is not separately restricted, and the mathematical model is as follows
Q ac =Q p COP ac
Wherein Q is ac The refrigeration power of the absorption refrigerator; q p Consuming thermal power for the absorption chiller; COP (coefficient of Performance) ac The conversion efficiency of the absorption refrigerator is obtained.
(1.3) the boiler equipment is a gas boiler, converts chemical energy into heat energy, acts as a standby heat source when the comprehensive energy gas turbine can not meet the heat load demand of a user, adds carbon emission right constraint, shows that the heat output of the gas boiler is not only limited by the rated heat power limit of the gas boiler equipment, but also limited by the carbon emission amount distributed to the gas boiler equipment, and the specific equivalent mathematical model is as follows:
Q GB,t =η GB V GB,t
0≤V GB,t ≤V GB,t,max
V GB ≤δ(C lim_GB +C lim_GB_buy )
C GB =aC lim_GB_buy
wherein Q is GB,t Outputting thermal power for the gas boiler; eta GB The conversion efficiency of the gas boiler is obtained; v GB,t Consuming power for gas boiler natural gas; v GB,t,max The maximum air inlet power of the gas boiler; c lim_GB Is the carbon emission limit allocated to the gas boiler plant.
C GB Represents the extra carbon cost for the supplier to supply heat with the gas boiler; c lim_GB_buy Represents the extra carbon emissions purchased by the supplier to ensure normal production heating; a is a carbon unit price.
(1.4) the mathematical expression of the equivalent model of the new energy power generation (green power) equipment is as follows:
(1.4.1) for photovoltaic power generation equipment, the output power mainly depends on the intensity of solar radiation, and the modeling process is as follows:
the distribution of the solar radiation intensity approximates to a Beta distribution with a distribution function of:
Figure GDA0003549211490000101
wherein G is t Is the illumination intensity at the time t (W/m 2); g max The maximum illumination intensity in the time period; b (alpha, Beta) is a Beta function, and alpha and Beta are shape parameters of a distribution function.
The photovoltaic power generation output mainly depends on the illumination intensity, and is influenced by the area of a battery plate, the temperature and the like, and the output power of the photovoltaic power generation can be estimated by the following formula:
P pv =A{P sc G t [1+k(T c -T r )]/G scPV
a is the area of the photovoltaic array; eta PV The photovoltaic conversion efficiency of the photovoltaic module is obtained; p sc The maximum output power of the photovoltaic power generation system under standard conditions (the illumination intensity is 1000W/m2, the ambient temperature is 25 ℃); k is a temperature coefficient, and is-0.47%/K; g sc The illumination intensity under the standard condition is 1000W/m 2; t is r For reference temperature, 25 ℃ is generally adopted;T c the calculation formula of the working temperature of the photovoltaic cell assembly is as follows:
T c =T a +30G t /1000
T a is ambient temperature
(1.4.2) for the wind power generation equipment, the output power of the wind power generation equipment mainly depends on the wind speed, the calculation of the wind speed is generally carried out by using a Weber distribution, and the mathematical expression of the wind power generation equipment is as follows:
Figure GDA0003549211490000111
wherein v is the wind speed; λ is a proportional parameter; k is a shape parameter
The wind power generation has three typical wind speed values, i.e. cut-in wind speed v in Rated wind speed v r And cut-out wind speed v out When the output of the wind generating set is modeled in a digital mode, the power curve needs to be approximately simulated, and when the rated power of a fan is P r,WT The equivalent numerical model is as follows:
Figure GDA0003549211490000112
for a new energy generation (green electricity) plant, considering its carbon yield, it is shown as follows:
C PV_WT =aδ PV_WT (P WT +P PV )Δt
C PV_WT revenue for new energy generation (green electricity), delta PV_WT A carbon emission factor.
(1.5) the electric heating equipment, namely the electric boiler equipment, wherein the mathematical expression of the output model of the electric heating equipment is as follows:
Q EB =η EB P EB
P EB_min ≤P EB ≤P EB_max
C EB =aδ EB P EB
wherein, P EB 、η EB Are respectively electric boilersThe heating equipment consumes electric energy and has conversion efficiency; c EB Additional carbon costs, delta, generated by using electric boiler plants for non-green electricity EB The carbon emission factor is equivalent to the energy consumption of the electric boiler.
(1.6) the electric refrigeration equipment utilizes electric energy to do work to transfer heat, in recent years, a Heat Pump (HP) attracts wide attention, can be directly installed on a user side of an integrated energy system to achieve the effects of heating in winter and cooling in summer, and can meet the requirements of cold and hot loads together with other energy supply equipment in the system, and the output model of the electric refrigeration equipment is as follows:
Q HP =η HP P HP
P HP_min ≤P HP ≤P HP_max
C HP =aδ HP P HP
wherein, P HP 、η HP Respectively consuming electric energy and converting efficiency for the electric refrigeration equipment; c HP The extra carbon cost, δ, of using electric refrigeration equipment for non-green electricity HP The carbon emission factor is equivalent to the energy consumption of the electric refrigeration equipment.
(1.7) the electrical hydrogen production equipment uses electrical energy to electrolyze water to generate hydrogen and oxygen through an electrolytic cell (SOEC), and the mathematical expression of the output model of the electrical hydrogen production equipment is as follows:
Figure GDA0003549211490000121
C SOEC =δ SOEC P SOEC
wherein, g SOEC The amount of hydrogen produced by the electrolytic cell in the current time interval;
Figure GDA0003549211490000122
is a low calorific value of hydrogen; p SOEC The electric energy consumed by the electrolytic cell in the current time interval; eta SOEC The electrolysis efficiency of the electrolytic cell; c SOEC The additional carbon cost of using electrical hydrogen production equipment for non-green electricity; delta SOEC Energy-consuming carbon for unit electric hydrogen production equipmentCost; according to the current practical production and application conditions, in order to ensure the maximum benefit, the hydrogen production by green electricity is only carried out when surplus green electricity exists.
(1.8) because the energy consumption factors which need to be considered by different users are different, and the influence on the optimization scheduling of the comprehensive energy system is different, the traditional user model is divided into large-scale industrial users, general industrial and commercial users and residential users, and different energy constraints are set.
(1.8.1) for large industrial users, there may be a carbon emission right requirement in addition to the regular energy requirement, i.e. additional carbon emission rights need to be purchased from the comprehensive energy service provider for normal production; in addition, the power consumption of the power generation system has a renewable energy power consumption responsibility weight index, namely the proportion of green power (new energy power generation) in the total power consumed by the power generation system needs to reach a specified responsibility weight index. According to the carbon emission right and the renewable energy power consumption responsibility index of the large industrial user, creating the energy demand constraint of the large industrial user, which can be expressed as:
C em,user1 =δ 1 P user1
P user1 =P user1_PV_WT_buy +P user1_PV_WT +P user1_buy
0≤C em,user1 ≤C em,user1_lim +C em,user1_buy
0≤ε 1 P user1 ≤P user1_PV_WT_buy +P user1_PV_WT
C user1 =aC em,user1_buy +bP user1_PV_WT_buy +cP user1_buy +d
the carbon emission of large industrial users is related to the production process, and the model estimates the carbon emission by using the total power consumption of enterprises, wherein C em,user1 Represents the carbon emissions of the large industrial user; delta 1 Represents a carbon emission conversion factor; p user1 Representing the total power consumption of the large industrial user; p user1_PV_WT_buy Representing green electric energy purchased by the large industrial user from the comprehensive energy service provider; p user1_PV_WT Represents the distributed photovoltaic and wind power generation of the large industrial userGreen electric energy is generated; p user1_buy Represents the conventional power purchased by the large industrial user; c em,user1_lim Indicating carbon emission rights assigned to the large industrial user; c em,user1_buy Representing the carbon emission right purchased by the large industrial user from the comprehensive energy service provider, and representing that the surplus carbon emission right is traded with the energy service provider when the value is negative; epsilon 1 Representing a renewable energy power consumption responsibility weight for the large industrial user; c user1 Represents the total cost of the production energy purchase of the large industrial user; aC em,user1_buy Represents a cost of purchasing carbon emissions rights; bP user1_PV_WT_buy Represents the cost of purchasing green electricity; cP (personal computer) user1_buy Represents the cost of purchasing traditional electrical energy; d represents the additional cost of energy expenditure for a large industrial user to purchase.
(1.8.2) for general industrial and commercial users, due to the relatively low energy consumption level and the large application of distributed energy, the spontaneous green power can basically meet the electric energy requirement, and the full green power is realized, and the energy requirement constraint of the general industrial and commercial users can be expressed as follows:
C em,user2 =δ 2 P user2
P user2 =P user2_PV_WT_buy +P user2_PV_WT
0≤C em,user2 ≤C em,user2_lim +C em,user2_buy
0≤ε 2 P user2 ≤P user2_PV_WT_buy +P user2_PV_WT
C user2 =aC em,user2_buy +bP user2_PV_WT_buy +c
the carbon emission of general industrial and commercial users is related to the production process, and the model estimates the carbon emission by using the total power consumption of enterprises, wherein C em,user2 Represents the carbon emission of the general industrial and commercial users; delta 2 Representing the carbon emission conversion coefficient thereof; p user2 Representing the total electric energy consumption of the general industrial and commercial users; p user2_PV_WT_buy The general industrial and commercial user is represented with the green electric energy purchased from the comprehensive energy service provider, and the value is generally 0 or negative according to the user type, which means that the general industrial and commercial user does not take the comprehensive energy serviceBuying green electric energy or selling surplus green electric energy to a comprehensive energy service provider by a service provider; p user2_PV_WT Representing the green electric energy generated by the distributed photovoltaic and wind power of the large industrial user; c em,user2_lim Indicating the carbon emission rights assigned to the general industrial and commercial user; c em,user2_buy The carbon emission right purchased from the general industrial and commercial user to the comprehensive energy service provider is represented, and when the value is negative, the surplus carbon emission right is represented to be traded with the energy service provider; epsilon 2 Representing a renewable energy power consumption responsibility weight of the general industrial and commercial user; c user2 Represents the total cost of the production energy purchase of the general industrial and commercial users; aC em,user2_buy Represents a cost of purchasing carbon emissions rights; bP user2_PV_WT_buy Represents the cost of purchasing green electricity; and c represents other costs for the general industry and commerce users to purchase energy.
(1.8.3) for the residential users, the carbon emission right and the green electricity consumption index are not considered, and due to the condition difference, the difference of the residential users' own distributed energy is large, the bidirectional flow of electric energy exists between the residential users and the energy service providers, the residential surplus green electric energy can be traded with the energy service providers, and the energy demand constraint of the residential users can be expressed as:
P user3 =P user3_PV_WT +P user3_PV_WT_buy
C user3 =aP user3_PV_WT_buy +b
not considering the carbon emission right and the green electricity consumption index of the residential user, wherein P user3 Representing the total electric energy consumption of the resident user; p user3_PV_WT_buy The method represents green electric energy purchased by the resident user from the comprehensive energy service provider, and the values are different according to different user types: if the value is positive, the resident user self distributed energy cannot meet self requirements and needs to purchase energy from the energy service provider, if the value is zero, the resident user does not purchase energy from the comprehensive energy service provider, and if the value is negative, the resident user self surplus green power is sold to the comprehensive energy service provider; p user3_PV_WT Representing the green electric energy generated by the distributed photovoltaic and wind power of the resident user; c user3 Indicating the total cost of energy purchase of the resident user;aP user3_PV_WT_buy Representing the cost of purchasing green electricity, and representing the surplus energy profit when the value is negative; b represents other energy cost of the residential user.
2) Constructing an energy supply network structure model of the single-park comprehensive energy system;
an energy supply network structure model of a single-park integrated energy system is constructed, as shown in the attached drawing 1, the integrated energy system network comprises four kinds of energy sources of cold, heat, electricity and hydrogen, and relates to equipment and technologies such as a gas turbine, a gas boiler, new energy power generation, absorption refrigeration, electricity hydrogen production, electricity heating, electricity refrigeration and the like, and energy demand constraints of integrated energy suppliers and service providers and different energy demands of large-scale industrial users, general industrial and commercial users and residential users can be comprehensively considered.
3) An energy service provider is used as a button, and energy price and power provided by the service provider are used as links to construct an outer-layer energy service provider optimized scheduling model;
and (3) constructing an outer energy service provider optimized dispatching model, as shown in the attached figure 2, taking each energy service provider in the park as a representative of the energy in the corresponding park, comprehensively summarizing the energy supply and demand conditions in the park, making energy transaction prices, and performing cooperative transaction with other energy service providers to form a regional multi-park comprehensive energy network.
4) And initializing an energy supply network model of the integrated energy system in each single park according to a preset value, carrying out optimized scheduling on the integrated energy system of the inner single park, optimizing the benefit of an energy supplier in the park by using an objective function of the optimized scheduling model of the inner park, and feeding back the energy price and the energy supply quantity to an outer energy supplier.
Figure GDA0003549211490000151
Figure GDA0003549211490000152
Figure GDA0003549211490000153
Figure GDA0003549211490000154
F is the total cost of the park comprehensive energy supplier;
Figure GDA0003549211490000155
cost of purchasing gas for the outside;
Figure GDA0003549211490000156
the unit price of the purchased gas; p gas,t The gas amount is purchased;
Figure GDA0003549211490000157
operating and maintaining costs for energy suppliers;
Figure GDA0003549211490000158
the operation and maintenance costs of each unit of output of the gas turbine and the gas boiler are respectively;
Figure GDA0003549211490000159
electric power and thermal power of a gas turbine and a gas boiler of an energy supplier in a time period t respectively;
Figure GDA00035492114900001510
the carbon cost for the production and operation of suppliers specifically comprises the carbon cost generated by the production of gas turbines, gas boilers, new energy power generation, electric heating and electric refrigeration equipment.
5) The energy service provider comprehensively manages the energy demand of users in the park, the energy price and power provided by the energy provider and other park energy service providers, constructs an outer energy service provider optimized dispatching model by taking the benefit optimization of the energy service provider as a target, and solves the optimized model by utilizing a particle swarm algorithm to obtain an optimal dispatching scheme.
And (5.1) comprehensive energy service providers of all parks in the area form an outer-layer energy service provider optimized scheduling model through information sharing, energy cooperation and transaction. According to the information provided by the comprehensive energy service provider, the energy demand and the energy price of each park in the regional system are known. And (3) aiming at maximizing the benefits of each park in the area, carrying out single-park energy optimization regulation and control, wherein the optimization objective function is as follows:
Figure GDA00035492114900001511
wherein the content of the first and second substances,
Figure GDA00035492114900001512
Figure GDA00035492114900001513
Figure GDA00035492114900001514
Figure GDA00035492114900001515
Figure GDA00035492114900001516
Figure GDA00035492114900001517
wherein N is a regional park set, which indicates that N parks of the regional comprehensive energy network participate in energy trading, and N is the number of the parks; f is the total income of the park comprehensive energy service provider;
Figure GDA00035492114900001518
the electricity selling income of the comprehensive energy service provider of the ith park is obtained for the time period t;
Figure GDA0003549211490000161
the electricity selling income of the integrated energy service business of the ith park is obtained for the time period t;
Figure GDA0003549211490000162
the comprehensive energy service provider of the ith park sells electricity to the power grid for the time t;
Figure GDA0003549211490000163
(ii) purchase energy costs from the ith campus energy provider for time period t;
Figure GDA0003549211490000164
the energy interaction cost of the ith park and other parks in the area is t;
Figure GDA0003549211490000165
the electricity purchasing cost from the power grid for the ith park comprehensive energy service provider in the time period t;
Figure GDA0003549211490000166
the unit electric energy income for the ith energy service provider and the nth regional park service provider,
Figure GDA0003549211490000167
The unit heat energy profit for the ith energy service provider and the nth park service provider,
Figure GDA0003549211490000168
The cost for the interaction of the energy service provider and the service providers of other parks in the area,
Figure GDA0003549211490000169
Heat energy interaction cost is provided for the ith energy service provider and the nth park service provider in the region; p in Electric energy interactive quantity H for ith energy service provider and nth regional park service provider in The heat energy interaction amount is provided for the ith energy service provider and the nth park service provider in the region;
Figure GDA00035492114900001610
time t the ith campusThe energy service provider can sell energy cost from the power grid unit,
Figure GDA00035492114900001611
Purchasing energy cost from a power grid unit for the ith park energy service provider in the time period t;
Figure GDA00035492114900001612
the electric energy interaction quantity between the energy service provider of the ith park and the power grid is t;
Figure GDA00035492114900001613
the unit electric energy cost purchased by the ith park energy service provider and the own park energy provider in the time period t,
Figure GDA00035492114900001614
The unit heat energy cost purchased by the ith park energy service provider and the park energy supplier in the time period t;
Figure GDA00035492114900001615
the electric energy purchased by the energy service provider of the park for the ith time period t through the energy provider of the park,
Figure GDA00035492114900001616
The heat energy purchased by the energy service provider of the ith park in the time period t through the energy provider of the park;
Figure GDA00035492114900001617
the unit price of electricity sold to the user by the ith campus energy service provider for the time period t,
Figure GDA00035492114900001618
The unit heat energy price sold to the user by the ith park energy service provider for the time period t; Δ t is the time span.
And (5.2) determining the external trading power composition of the park according to the external trading power demand constraint and the energy supply constraint of other parks in the region by using the external trading (selling or purchasing) maximum profit of the park as an objective function through the external trading optimization strategy. The objective function is:
Figure GDA00035492114900001619
wherein the content of the first and second substances,
Figure GDA00035492114900001620
Figure GDA0003549211490000171
Figure GDA0003549211490000172
Figure GDA0003549211490000173
Figure GDA0003549211490000174
Figure GDA0003549211490000175
wherein, P out ,H out Respectively carrying out external transaction on electric energy and heat energy, including transaction with other regional energy service providers and energy transaction with a network;
Figure GDA0003549211490000176
representing the thermal energy purchased by the ith campus energy service provider from the grid unit for time period t.
And (5.3) solving the optimized dispatching model of the external energy service provider by utilizing a particle swarm algorithm to obtain the optimal dispatching scheme of the comprehensive energy system.
In the optimization solution of the objective function, common optimization algorithms include a Particle Swarm Optimization (PSO), a Genetic Algorithm (GA), an improved non-dominated sorting genetic algorithm (NSGA-II), an ant colony algorithm, and the like, where the particle swarm optimization algorithm has a simple structure, a fast convergence speed, and is easy to implement and intelligent, and can perform optimization solution through iteration according to a multidimensional constraint condition in a calculation process, and the method is widely applied to solving optimization problems of load economic distribution, power grid planning, and the like of a power system.
Particle swarm optimization designs a particle without mass, and the particle has only two properties: speed, which represents how fast the movement is, and position, which represents the direction of the movement. And each particle independently searches an optimal solution in the search space, records the optimal solution as a current individual extremum, shares the individual extremum with other particles in the whole particle swarm, finds the optimal individual extremum as a current global optimal solution of the whole particle swarm, and adjusts the speed and the position of each particle in the particle swarm according to the found current individual extremum and the shared current global optimal solution of the whole particle swarm.
Initially, a group of random particles (random solution) is present, and in each iteration, the particles update their speed and position by tracking extreme values (pbest, gbest) according to the following formula, and then find the optimal solution by iteration.
v i =ω*v i +c 1 *rand()*(pbest i -x i )+c 2 *rand()*(gbest i -x i )
x i =x i +v i
Wherein i is 1,2,3 … … N, N is the total number of the particle groups; omega is an inertia factor and is more than or equal to 0, when the value of omega is larger, the global optimizing capability of the algorithm is strong, the local optimizing capability of the algorithm is weak, and when the value of omega is smaller, the global optimizing capability of the algorithm is weak, and the local optimizing capability of the algorithm is strong; v. of i Is the velocity of the particle, and v i ≤v max (ii) a rand () is a random number between (0, 1); x is the number of i Is the current position of the particle; c. C 1 、c 2 Are learning factors, usually equal and take the value 2;
in the practical application of the particle swarm algorithm, the inertia factor in the formula mostly adopts a dynamic value to adjust the global and local searching capability of the algorithm, and the determining method is as follows:
ω (t) =(ω intend )(G k -g)/G kend
wherein ω is int The value of the initial inertia factor is 0.9 as a typical value; omega end The value of the inertia factor is the value when the iteration reaches the maximum algebra, and the typical value is 0.4; g k Is the maximum number of iterations.
The solving process of the regional multi-park energy dispatching optimization model is shown in the attached drawing 3, and the specific solving process of the particle swarm algorithm is shown in the attached drawing 4.
The embodiment provides an energy optimization scheduling system considering carbon emission and green power consumption, which is used for executing the energy optimization scheduling method considering carbon emission and green power consumption in the foregoing method embodiments. Since the technical principle of the system embodiment is similar to that of the method embodiment, repeated description of the same technical details is omitted.
In one embodiment, an energy optimization scheduling system considering carbon emissions and green power consumption includes:
the model and constraint creation module is used for constructing a power output model of the energy supply equipment according to the operating characteristics of the energy supply equipment in the park and the carbon emission right; configuring carbon emission rights and green electricity consumption indexes for different user types, and creating energy demand constraints of output models of energy supply equipment according to different user types; wherein, the functional device output model comprises: the method comprises the following steps of (1) generating a new energy power generation system by using a gas turbine output model, a boiler equipment output model, a new energy power generation equipment output model, an electric heating equipment output model, an electric refrigerating equipment output model and an electric hydrogen production equipment output model; the user types include: large industrial users, general industrial and commercial users, and residential users;
and the optimization calculation module is used for constructing an outer energy service provider optimization scheduling model by taking the optimal income of the energy provider as a target according to the energy demand constraint and the energy prices and powers of the current park energy provider and other park energy providers, and solving the outer energy service provider optimization scheduling model through an optimization algorithm to obtain an energy optimal scheduling scheme.
The foregoing embodiments are merely illustrative of the principles and utilities of the present invention and are not intended to limit the invention. Any person skilled in the art can modify or change the above-mentioned embodiments without departing from the spirit and scope of the present invention. Accordingly, it is intended that all equivalent modifications or changes which can be made by those skilled in the art without departing from the spirit and technical spirit of the present invention be covered by the claims of the present invention.

Claims (5)

1. An energy optimization scheduling method considering carbon emission and green electricity consumption is characterized by comprising the following steps:
constructing a power output model of the energy supply equipment according to the operating characteristics of the energy supply equipment in the park and the carbon emission right; configuring carbon emission rights and green electricity consumption indexes for different user types, and creating energy demand constraints of output models of energy supply equipment according to different user types; wherein, the output model of the energy supply equipment comprises: the system comprises a gas turbine output model, a boiler equipment output model, a new energy power generation equipment output model, an electric heating equipment output model, an electric refrigerating equipment output model and an electric hydrogen production equipment output model; the user types include: large industrial users, general industrial and commercial users, and residential users; the gas turbine output model is as follows:
V GT ≤δ(C lim_GT +C lim_GT_buy )
C GT =aC lim_GT_buy
in the formula: v GT Natural gas consumption for gas turbines; is the carbon emission coefficient, C lim_GT To allocate carbon emission limits; c GT Represents the extra carbon cost for the energy supplier to generate electricity using the gas turbine; c lim_GT_buy Represents the extra carbon emission rights purchased by the energy supplier to ensure normal production; a is a carbon unit price;
the boiler equipment comprises a gas boiler, and the output model of the boiler equipment is as follows:
V GB ≤δ(C lim_GB +C lim_GB_buy )
C GB =aC lim_GB_buy
wherein, V GB Consuming power for gas boiler natural gas; c lim_GB A carbon emission allowance assigned to the gas boiler plant; c GB Represents the extra carbon cost for the supplier to supply heat with the gas boiler; c lim_GB_buy Represents the extra carbon emission rights purchased by the energy supplier to ensure normal production heating;
the new energy power generation equipment comprises solar power generation equipment and wind power generation equipment, and the output model of the new energy power generation equipment is as follows:
C PV_WT =aδ PV_WT (P WT +P PV )Δt
wherein, C PV_WT Income from new energy generation, delta PV_WT A carbon emission factor; p WT Outputting power for wind power generation; p PV Outputting power for photovoltaic power generation; Δ t is the time span;
the electric heating equipment comprises electric boiler equipment, and the output model of the electric heating equipment is as follows:
C EB =aδ EB P EB
wherein, P EB Consuming electrical energy for the electrical heating device; c EB Additional carbon cost, δ, of using electrical heating equipment for non-green electricity EB Energy consumption equivalent carbon emission factors for the electric heating equipment;
the output model of the electric refrigeration equipment is as follows:
C HP =aδ HP P HP
wherein, P HP Consuming electrical energy for the electrical refrigeration equipment; c HP Additional carbon cost, δ, of using electric refrigeration equipment for non-green electricity HP Energy consumption equivalent carbon emission factors for the electric refrigeration equipment;
the output model of the electrical hydrogen production equipment is as follows:
C SOEC =δ SOEC P SOEC
wherein, P SOEC Electrical energy consumed by the electrolytic cell; c SOEC Electric hydrogen production facility for non-green electricityExtra carbon cost for production; delta SOEC Energy consumption and carbon cost of unit electric hydrogen production equipment;
according to the energy demand constraint and the energy prices and powers of the current park energy suppliers and other park energy suppliers, an outer energy service provider optimized dispatching model is constructed by taking optimal income of the energy suppliers as a target, and the outer energy service provider optimized dispatching model is solved through an optimization algorithm to obtain an energy optimal dispatching scheme, which comprises the following steps: setting comprehensive energy service providers of all parks in the area to perform external energy trade through information sharing, and constructing an outer-layer energy service provider optimized scheduling model; each regional park service provider carries out optimization regulation and control by taking the benefit maximization of the regional park service provider as a target, and the optimization target function is as follows:
Figure FDA0003676106300000021
wherein the content of the first and second substances,
Figure FDA0003676106300000022
Figure FDA0003676106300000023
Figure FDA0003676106300000024
Figure FDA0003676106300000025
Figure FDA0003676106300000026
Figure FDA0003676106300000027
the method comprises the following steps that N is a regional park set, N shows that N parks of a regional comprehensive energy network participate in energy trading, and N is the number of the parks; f is the total income of the park comprehensive energy service provider;
Figure FDA0003676106300000028
the electricity selling income of the comprehensive energy service provider of the ith park is obtained for the time period t;
Figure FDA0003676106300000029
the electricity selling income of the integrated energy service business of the ith park is obtained for the time period t;
Figure FDA00036761063000000210
the comprehensive energy service provider of the ith park sells electricity to the power grid for the time t;
Figure FDA00036761063000000211
purchasing energy costs from the ith campus energy supplier for time period t;
Figure FDA00036761063000000212
the energy interaction cost of the ith park and other parks in the area is t;
Figure FDA0003676106300000031
the electricity purchasing cost from the power grid for the ith park comprehensive energy service provider in the time period t;
Figure FDA0003676106300000032
the unit electric energy income for the ith energy service provider and the nth park service provider,
Figure FDA0003676106300000033
The unit heat energy profit for the ith energy service provider and the nth park service provider,
Figure FDA0003676106300000034
The cost for the interaction of the energy service provider and the service providers of other parks in the area,
Figure FDA0003676106300000035
Heat energy interaction cost is provided for the ith energy service provider and the nth park service provider in the region; p in Electric energy interaction quantity H for ith energy service provider and nth regional park service provider in The heat energy interaction amount is provided for the ith energy service provider and the nth park service provider in the region;
Figure FDA0003676106300000036
the energy selling cost of the ith park energy service provider from the power grid unit for the time period t,
Figure FDA0003676106300000037
Purchasing energy cost from a power grid unit for the ith park energy service provider in the time period t;
Figure FDA0003676106300000038
the electric energy interaction quantity between the energy service provider of the ith park and the power grid is t;
Figure FDA0003676106300000039
the unit electric energy cost purchased by the ith park energy service provider and the own park energy provider in the time period t,
Figure FDA00036761063000000310
The unit heat energy cost purchased by the ith park energy service provider and the energy supplier of the park at the time t;
Figure FDA00036761063000000311
the electric energy purchased by the energy service provider of the park for the ith time period t through the energy provider of the park,
Figure FDA00036761063000000312
For a time period t the ith energy of the parkThe source service provider purchases heat energy through the energy supplier of the park;
Figure FDA00036761063000000313
the unit price of electricity sold to the user by the ith campus energy service provider for the time period t,
Figure FDA00036761063000000314
The unit price of heat energy sold to the user by the ith campus energy service provider for time period t.
2. The energy optimization scheduling method considering carbon emission and green electricity consumption according to claim 1, wherein carbon emission rights and green electricity consumption indexes are configured for different user types, and energy demand constraints of output models of energy supply devices are created according to different user types, and the method comprises the following steps:
according to the carbon emission weight and the renewable energy power consumption responsibility weight index of the large industrial user, creating an energy demand constraint of the large industrial user, wherein the energy demand constraint is represented as:
C em,user1 =δ 1 P user1
P user1 =P user1_PV_WT_buy +P user1_PV_WT +P user1_buy
0≤C em,user1 ≤C em,user1_lim +C em,user1_buy
0≤ε 1 P user1 ≤P user1_PV_WT_buy +P user1_PV_WT
C user1 =aC em,user1_buy +bP user1_PV_WT_buy +cP user1_buy +d
wherein the carbon emissions, C, are evaluated on the basis of the total power consumption of large industrial users em,user1 Represents the carbon emissions of the large industrial user; delta 1 Represents a carbon emission conversion factor; p user1 Representing the total electric energy consumption of the large industrial user; p user1_PV_WT_buy Representing green electric energy purchased by the large industrial user from an integrated energy service provider;
Figure FDA0003676106300000041
P user1_PV_WT green electric energy representing self distributed photovoltaic and wind power generation of the large industrial users; p user1_buy Represents conventional electrical energy purchased by the large industrial user; c em,user1_lim Representing carbon emission rights assigned to the large industrial user; c em,user1_buy Representing the carbon emission rights purchased by the large industrial user from the comprehensive energy service provider, and when the value is negative, representing that the surplus carbon emission rights are traded with the energy service provider; epsilon 1 Representing a renewable energy power consumption liability weight for the large industrial user; c user1 Represents the total cost of production energy purchase for the large industrial user; aC em,user1_buy Represents a cost of purchasing carbon emissions rights; bP user1_PV_WT_buy Represents the cost of purchasing green electricity; cP (personal computer) user1_buy Represents the cost of purchasing traditional electrical energy; d represents other costs for the purchase of energy by large industrial users;
the energy demand constraints of the general industrial and commercial users are as follows:
C em,user2 =δ 2 P user2
P user2 =P user2_PV_WT_buy +P user2_PV_WT
0≤C em,user2 ≤C em,user2_lim +C em,user2_buy
0≤ε 2 P user2 ≤P user2_PV_WT_buy +P user2_PV_WT
C user2 =aC em,user2_buy +bP user2_PV_WT_buy +c
wherein, C em,user2 Represents the general industrial and commercial user carbon emission; delta 2 Representing the carbon emission conversion coefficient of the general industrial and commercial users; p user2 Representing the total electric energy consumption of the general industrial and commercial users; p user2_PV_WT_buy The green electric energy which is purchased from the comprehensive energy service provider by the general industrial and commercial user is represented, the value is 0 or a negative value, and the green electric energy is not purchased from the comprehensive energy service provider or surplus green electric power is sold to the comprehensive energy service provider; p user2_PV_WT Green electric energy representing the self distributed photovoltaic and wind power generation of the general industrial and commercial users; c em,user2_lim Representing carbon emission rights assigned to the general industrial and commercial user; c em,user2_buy The carbon emission right purchased from the general industrial and commercial user to the comprehensive energy service provider is represented, and when the value is negative, the surplus carbon emission right is represented to be traded with the energy service provider; epsilon 2 A renewable energy power consumption responsibility weight representing the general industrial and commercial user; c user2 Representing the total cost of the general industrial and commercial user production energy purchase; aC em,user2_buy Represents a cost of purchasing carbon emissions rights; bP user2_PV_WT_buy Represents the cost of purchasing green electricity; c represents other costs for the general industry and commerce users to purchase energy;
the energy demand constraint of the residential user is as follows:
P user3 =P user3_PV_WT +P user3_PV_WT_buy
C user3 =aP user3_PV_WT_buy +b
not considering the carbon emission right and the green electricity consumption index of the residential user, wherein P user3 Representing the total power consumption of the residential user; p user3_PV_WT_buy The method comprises the steps that green electric energy purchased from a comprehensive energy service provider by a resident user is represented, if the value is positive, the distributed energy of the resident user cannot meet the self demand, the energy is required to be purchased from the energy service provider, if the value is zero, the energy is not purchased from the comprehensive energy service provider, and if the value is negative, the surplus green electric power of the resident user is sold to the comprehensive energy service provider; p user3_PV_WT Representing the green electric energy generated by the distributed photovoltaic and wind power of the resident users; c user3 Representing the total cost of energy purchase by the residential user; aP user3_PV_WT_buy Representing the cost of purchasing green electricity, and representing the surplus energy profit when the value is negative; b represents other energy cost of the residential user.
3. The energy optimization scheduling method considering carbon emission and green electricity consumption according to claim 1, wherein the energy external trade is performed by the comprehensive energy service providers of the parks in the set area through information sharing, and comprises the following steps:
creating external trading power demand constraints and energy supply constraints for the campus and other camps in the campus, expressed as maximum sales or purchase revenue objective function for the campus and other camps
Figure FDA0003676106300000051
Wherein the content of the first and second substances,
Figure FDA0003676106300000052
Figure FDA0003676106300000053
Figure FDA0003676106300000054
Figure FDA0003676106300000055
Figure FDA0003676106300000056
Figure FDA0003676106300000057
wherein, P out ,H out Respectively carrying out external transaction on electric energy and heat energy, including transaction with other regional energy service providers and energy transaction with a network;
Figure FDA0003676106300000061
to representTime t thermal energy purchased by the ith campus energy service provider from the grid unit.
4. The energy optimization scheduling method considering carbon emission and green power consumption according to claim 1, wherein the optimization algorithm comprises a particle swarm algorithm or an ant colony algorithm.
5. An energy optimization scheduling system considering carbon emission and green electricity consumption, comprising:
the model and constraint creation module is used for constructing a power output model of the energy supply equipment according to the operating characteristics of the energy supply equipment in the park and the carbon emission right; configuring carbon emission rights and green electricity consumption indexes for different user types, and creating energy demand constraints of output models of energy supply equipment according to different user types; wherein, the output model of the energy supply equipment comprises: the system comprises a gas turbine output model, a boiler equipment output model, a new energy power generation equipment output model, an electric heating equipment output model, an electric refrigerating equipment output model and an electric hydrogen production equipment output model; the user types include: large industrial users, general industrial and commercial users, and residential users; the gas turbine output model is as follows:
V GT ≤δ(C lim_GT +C lim_GT_buy )
C GT =aC lim_GT_buy
in the formula: v GT Natural gas consumption for gas turbines; is the carbon emission coefficient, C lim_GT To allocate carbon emission limits; c GT Represents the extra carbon cost for the energy supplier to generate electricity using the gas turbine; c lim_GT_buy Represents the extra carbon emission rights purchased by the energy supplier to ensure normal production; a is a carbon unit price;
the boiler equipment comprises a gas boiler, and the output model of the boiler equipment is as follows:
V GB ≤δ(C lim_GB +C lim_GB_buy )
C GB =aC lim_GB_buy
wherein, V GB Is natural for gas boilerGas consumption power; c lim_GB A carbon emission allowance assigned to the gas boiler plant; c GB Represents the extra carbon cost for the supplier to supply heat with the gas boiler; c lim_GB_buy Represents the extra carbon emission rights purchased by the energy supplier to ensure normal production heating;
the new energy power generation equipment comprises solar power generation equipment and wind power generation equipment, and the output model of the new energy power generation equipment is as follows:
C PV_WT =aδ PV_WT (P WT +P PV )Δt
wherein, C PV_WT Revenue of new energy generation, delta PV_WT A carbon emission factor; p is WT Outputting power for wind power generation; p is PV Outputting power for photovoltaic power generation; Δ t is the time span;
the electric heating equipment comprises electric boiler equipment, and the output model of the electric heating equipment is as follows:
C EB =aδ EB P EB
wherein, P EB Consuming electrical energy for the electrical heating device; c EB Additional carbon cost, δ, of using electrical heating equipment for non-green electricity EB Energy consumption equivalent carbon emission factors for the electric heating equipment;
the output model of the electric refrigeration equipment is as follows:
C HP =aδ HP P HP
wherein, P HP Consuming electrical energy for the electrical refrigeration equipment; c HP The extra carbon cost, δ, of using electric refrigeration equipment for non-green electricity HP Energy consumption equivalent carbon emission factors for the electric refrigeration equipment;
the output model of the electrical hydrogen production equipment is as follows:
C SOEC =δ SOEC P SOEC
wherein, P SOEC Electrical energy consumed by the electrolytic cell; c SOEC The additional carbon cost of using electrical hydrogen production equipment for non-green electricity; delta SOEC Energy consumption and carbon cost of unit electric hydrogen production equipment;
the optimization calculation module is used for constructing an outer energy service provider optimization scheduling model by taking optimal income of the energy provider as a target according to the energy demand constraint and the energy prices and powers of the current park energy provider and other park energy providers, and solving the outer energy service provider optimization scheduling model through an optimization algorithm to obtain an energy optimal scheduling scheme, and the optimization calculation module comprises the following steps: setting comprehensive energy service providers of all parks in the area to perform external energy trade through information sharing, and constructing an outer-layer energy service provider optimized scheduling model; each regional park service provider carries out optimization regulation and control by taking the benefit maximization of the regional park service provider as a target, and the optimization target function is as follows:
Figure FDA0003676106300000071
wherein the content of the first and second substances,
Figure FDA0003676106300000072
Figure FDA0003676106300000073
Figure FDA0003676106300000074
Figure FDA0003676106300000075
Figure FDA0003676106300000076
Figure FDA0003676106300000077
wherein N is a set of regional parks representing the regional integrationThe energy network has N parks participating in energy trading, and N is the number of the parks; f is the total income of the park comprehensive energy service provider;
Figure FDA0003676106300000081
the electricity selling income of the comprehensive energy service provider of the ith park is obtained for the time period t;
Figure FDA0003676106300000082
the electricity selling income of the integrated energy service business of the ith park is obtained for the time period t;
Figure FDA0003676106300000083
the comprehensive energy service provider of the ith park sells electricity to the power grid for the time t;
Figure FDA0003676106300000084
purchasing energy costs from the ith campus energy supplier for time period t;
Figure FDA0003676106300000085
the energy interaction cost of the ith park and other parks in the area is t;
Figure FDA0003676106300000086
the electricity purchasing cost from the power grid for the ith park comprehensive energy service provider in the time period t;
Figure FDA0003676106300000087
the unit electric energy income for the ith energy service provider and the nth park service provider,
Figure FDA0003676106300000088
The unit heat energy profit for the ith energy service provider and the nth park service provider,
Figure FDA0003676106300000089
The cost for the interaction of the electric energy of the energy service provider and the service providers of other parks in the area,
Figure FDA00036761063000000810
Heat energy interaction cost is provided for the ith energy service provider and the nth park service provider in the region; p in Electric energy interaction quantity H for ith energy service provider and nth regional park service provider in The heat energy interaction amount is provided for the ith energy service provider and the nth park service provider in the region;
Figure FDA00036761063000000811
the energy selling cost of the ith park energy service provider from the power grid unit for the time period t,
Figure FDA00036761063000000812
Purchasing energy cost from a power grid unit for the ith park energy service provider in the time period t;
Figure FDA00036761063000000813
the electric energy interaction quantity between the energy service provider of the ith park and the power grid is t;
Figure FDA00036761063000000814
the unit electric energy cost purchased by the ith park energy service provider and the own park energy provider in the time period t,
Figure FDA00036761063000000815
The unit heat energy cost purchased by the ith park energy service provider and the park energy supplier in the time period t;
Figure FDA00036761063000000816
the electric energy bought by the energy supplier of the park for the ith park energy service provider in the time period t,
Figure FDA00036761063000000817
The heat energy purchased by the energy service provider of the ith park in the time period t through the energy provider of the park;
Figure FDA00036761063000000818
the unit price of electricity sold to the user by the ith campus energy service provider for the time period t,
Figure FDA00036761063000000819
The unit price of heat energy sold to the user by the ith campus energy service provider for time period t.
CN202111005567.4A 2021-08-30 2021-08-30 Energy optimization scheduling method and system considering carbon emission and green electricity consumption Active CN113592365B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111005567.4A CN113592365B (en) 2021-08-30 2021-08-30 Energy optimization scheduling method and system considering carbon emission and green electricity consumption

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111005567.4A CN113592365B (en) 2021-08-30 2021-08-30 Energy optimization scheduling method and system considering carbon emission and green electricity consumption

Publications (2)

Publication Number Publication Date
CN113592365A CN113592365A (en) 2021-11-02
CN113592365B true CN113592365B (en) 2022-08-05

Family

ID=78240210

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111005567.4A Active CN113592365B (en) 2021-08-30 2021-08-30 Energy optimization scheduling method and system considering carbon emission and green electricity consumption

Country Status (1)

Country Link
CN (1) CN113592365B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115577234B (en) * 2022-11-21 2023-04-07 中国电力科学研究院有限公司 Node power supply emission factor calculation method and system based on power flow distribution
CN117220346B (en) * 2023-07-27 2024-04-16 河海大学 Comprehensive energy service business electricity-carbon-green certificate double-layer distributed scheduling method
CN116780535B (en) * 2023-08-16 2024-01-02 国网浙江省电力有限公司金华供电公司 Light-storage collaborative optimization scheduling method based on ladder-type carbon transaction mechanism

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100042453A1 (en) * 2008-08-12 2010-02-18 Efficiency 2.0, LLC. Methods and apparatus for greenhouse gas footprint monitoring
EP2387776A4 (en) * 2009-01-14 2013-03-20 Integral Analytics Inc Optimization of microgrid energy use and distribution
CN107194514B (en) * 2017-05-27 2020-08-18 重庆大学 Demand response multi-time scale scheduling method for wind power prediction error
CN110889549B (en) * 2019-11-21 2022-08-05 国网江苏省电力有限公司经济技术研究院 Multi-objective optimization scheduling method of comprehensive energy system considering human comfort
CN112036747A (en) * 2020-08-31 2020-12-04 国网河南省电力公司经济技术研究院 Evaluation method of park comprehensive energy system multi-demand response implementation model
CN112116476B (en) * 2020-09-23 2024-03-01 中国农业大学 Comprehensive energy system simulation method considering wind power and carbon transaction mechanism

Also Published As

Publication number Publication date
CN113592365A (en) 2021-11-02

Similar Documents

Publication Publication Date Title
Wang et al. Operation optimization of regional integrated energy system based on the modeling of electricity-thermal-natural gas network
Luo et al. Multi-objective capacity optimization of a distributed energy system considering economy, environment and energy
Liu et al. Two-phase collaborative optimization and operation strategy for a new distributed energy system that combines multi-energy storage for a nearly zero energy community
CN113592365B (en) Energy optimization scheduling method and system considering carbon emission and green electricity consumption
CN108206543B (en) Energy router based on energy cascade utilization and operation optimization method thereof
Shen et al. Multi-objective capacity configuration optimization of an integrated energy system considering economy and environment with harvest heat
CN111463836B (en) Comprehensive energy system optimal scheduling method
CN108009693A (en) Grid-connected micro-capacitance sensor dual blank-holder based on two-stage demand response
Guo et al. Two-layer co-optimization method for a distributed energy system combining multiple energy storages
CN109861302B (en) Master-slave game-based energy internet day-ahead optimization control method
Fan et al. Two-layer collaborative optimization for a renewable energy system combining electricity storage, hydrogen storage, and heat storage
Zhang et al. Energy scheduling optimization of the integrated energy system with ground source heat pumps
CN113779783A (en) Multi-uncertainty-considered planning and operation joint optimization method for regional comprehensive energy system
Song et al. A fuzzy‐based multi‐objective robust optimization model for a regional hybrid energy system considering uncertainty
CN116432824A (en) Comprehensive energy system optimization method and system based on multi-target particle swarm
Wang et al. Co-optimization of configuration and operation for distributed multi-energy system considering different optimization objectives and operation strategies
Li et al. Intraday multi-objective hierarchical coordinated operation of a multi-energy system
Dong et al. Hierarchical multi-objective planning for integrated energy systems in smart parks considering operational characteristics
CN116468215A (en) Comprehensive energy system scheduling method and device considering uncertainty of source load
Wu et al. Multi-parameter cooperative optimization and solution method for regional integrated energy system
Wu et al. Multi-parameter optimization design method for energy system in low-carbon park with integrated hybrid energy storage
CN114742276A (en) Multi-objective optimization scheduling method of park comprehensive energy system with ORC (organic Rankine cycle) considering exergy efficiency
Mahoor et al. Smart energy management for a micro-grid with consideration of demand response plans
Wang et al. Strategy and capacity optimization of renewable hybrid combined cooling, heating and power system with multiple energy storage
Ma et al. Performance optimization of phase change energy storage combined cooling, heating and power system based on GA+ BP neural network algorithm

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