CN115239538A - Comprehensive energy park low-carbon operation optimization method based on dynamic carbon emission factor - Google Patents

Comprehensive energy park low-carbon operation optimization method based on dynamic carbon emission factor Download PDF

Info

Publication number
CN115239538A
CN115239538A CN202210677812.4A CN202210677812A CN115239538A CN 115239538 A CN115239538 A CN 115239538A CN 202210677812 A CN202210677812 A CN 202210677812A CN 115239538 A CN115239538 A CN 115239538A
Authority
CN
China
Prior art keywords
carbon
carbon emission
energy
park
dynamic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210677812.4A
Other languages
Chinese (zh)
Inventor
姚力
陆春光
刘炜
肖涛
李亦龙
王朝亮
王佳颖
宋磊
黄荣国
蒋群
俞佳莉
韩鑫泽
迟长云
蒋雪冬
李晓波
汪超群
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Marketing Service Center of State Grid Zhejiang Electric Power Co Ltd
Original Assignee
Zhejiang University ZJU
Marketing Service Center of State Grid Zhejiang Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU, Marketing Service Center of State Grid Zhejiang Electric Power Co Ltd filed Critical Zhejiang University ZJU
Priority to CN202210677812.4A priority Critical patent/CN115239538A/en
Publication of CN115239538A publication Critical patent/CN115239538A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10—Services
    • G06Q50/26—Government or public services
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/80—Management or planning
    • Y02P90/84—Greenhouse gas [GHG] management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Theoretical Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Mathematical Analysis (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Development Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Public Health (AREA)
  • Computing Systems (AREA)
  • Game Theory and Decision Science (AREA)
  • Educational Administration (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Water Supply & Treatment (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention belongs to the technical field of carbon emission, and particularly relates to a comprehensive energy park low-carbon operation optimization method based on a dynamic carbon emission factor. Aiming at the defect that the carbon emission factor is considered to be rough and shallow by the existing comprehensive energy park operation optimization method, the invention adopts the following technical scheme: a low-carbon operation optimization method for a comprehensive energy park based on dynamic carbon emission factors comprises the following steps: s1, acquiring dynamic carbon elimination factors of province areas or local cities of parks, which change along with time; s2, establishing a target function related to the dynamic carbon emission factor by taking the minimum carbon emission of the park as a target; and S3, optimizing according to the relevant constraint conditions. The beneficial effects of the invention are: according to the equipment operation and the energy storage equipment charge-discharge of more accurate instruction comprehensive energy garden of developments carbon row factor, use clean electric energy more when electric power carbon row factor is lower in order to reduce the carbon total emission volume in garden.

Description

Comprehensive energy park low-carbon operation optimization method based on dynamic carbon emission factor
Technical Field
The invention belongs to the technical field of carbon emission, and particularly relates to a comprehensive energy park low-carbon operation optimization method based on a dynamic carbon emission factor.
Background
China clearly proposes ambitious targets of 'carbon peak reaching' and 'carbon neutralization', and the 'double-carbon' target promotes industrial structure and energy structure adjustment, and changes the mode of green low-carbon energy supply and consumption while considering economic development. The aim of 'double carbon' also puts higher requirements on accelerating diversification, cleanness, low carbonization, high efficiency of energy consumption and the like of the propulsion energy supply of a power grid company.
More than 70% of industrial energy in China is concentrated in a comprehensive energy park, the power consumption is high, the energy utilization modes of users are diversified, and a large amount of cold and heat demands are met. How to strengthen the multi-energy cooperative supply, the energy cascade utilization and the low-carbon high-efficiency operation of the comprehensive energy park becomes a new focus and research focus.
In the existing comprehensive energy park operation optimization method, the aim of optimizing the economical efficiency is always taken, a model considering the carbon emission only directly substitutes the regional power carbon emission factor or converts the carbon value into the economic aim, and the dynamic change characteristic of the carbon emission factor under the high-speed development of new energy is not considered, so that the operation cannot be optimized according to the dynamic change characteristic of the carbon emission factor.
Disclosure of Invention
The invention provides a comprehensive energy park low-carbon operation optimization method based on dynamic carbon emission factors, aiming at the defect that the existing comprehensive energy park operation optimization method considers the carbon emission factors roughly and shallowly, and aiming at the lowest carbon emission, considering the equipment operation constraint and the energy balance constraint, optimizing the equipment operation mode of the comprehensive energy park by using the dynamic carbon emission factors to obtain the unit operation mode combination with the lowest carbon emission.
In order to realize the purpose, the invention adopts the following technical scheme: a comprehensive energy park low-carbon operation optimization method based on dynamic carbon emission factors comprises the following steps:
s1, acquiring dynamic carbon emission factors of provinces or local cities of parks along with time change;
s2, establishing a target function related to the dynamic carbon emission factor by taking the minimum carbon emission of the garden as a target;
and S3, optimizing according to the relevant constraint conditions.
According to the comprehensive energy park low-carbon operation optimization method based on the dynamic carbon emission factor, the dynamic carbon emission factor is obtained, an objective function about the dynamic carbon emission factor is established, the dynamic carbon emission factor can more accurately depict the total carbon emission amount of the comprehensive energy park, and the phenomenon of the low-power carbon emission factor when a new energy generating set such as photovoltaic and wind power generation sets generates high power is described; according to the equipment operation and the energy storage equipment charge-discharge of more accurate instruction comprehensive energy garden of developments carbon row factor, use clean electric energy more when electric power carbon row factor is lower in order to reduce the carbon total emission volume in garden.
Further, in step S1, the calculation formula of the provincial dynamic carbon row factor is as follows:
Figure BDA0003697108800000021
in the formula: f is dynamic carbon emission factor of provincial power supply side and has a unit of tCO 2 /kWh; i is the energy type of the power-saving power supply side regulating unit, and comprises coal-fired power generation, gas power generation, hydroelectric power generation, photovoltaic power generation, wind power generation and nuclear power generation; w is an external electricity region; c is the type of the coal burner group; g is a gas engine set type; g is generated energy, the unit is kWh, and the value is the accumulated value of the gateway metering for 15 min; EF c Is the carbon emission factor of the coal-fired unit and has the unit of tCO 2 /kWh；EF g Is the carbon emission factor of the gas turbine set and has the unit of tCO 2 /kWh；EF w Average carbon exclusion factor for foreign electric regions; t represents time.
Further, after the provincial region dynamic carbon emission factor is obtained, the carbon flow conditions of various cities are continuously refined according to the provincial power topological graph, and the dynamic carbon emission factor of the city to which the comprehensive energy park belongs is obtained.
Further, in step S2, the objective function is:
Figure BDA0003697108800000022
in the formula: f is the carbon emission of the park;
Figure BDA0003697108800000023
carbon emission coefficient per unit of electrical energy purchased from the grid for IES, which dynamically changes over time, tco 2 /kWh；P FromGrid To purchase electricity from the grid, kWh;
Figure BDA0003697108800000024
carbon emission coefficient, tco, for gas-fired plants 2 /kWh；
Figure BDA0003697108800000025
Carbon emission coefficient, tCO, for coal burning plants 2 /kWh；
Figure BDA0003697108800000026
kW is the output power of the gas equipment;
Figure BDA0003697108800000027
kW is the output power of the coal burning equipment; t represents time; i is the device type; j is the equipment serial number; the carbon emission coefficient of each device is the attribute of the device, and has specificity.
Figure BDA0003697108800000028
Namely the provincial dynamic carbon elimination factor or the local city dynamic carbon elimination factor.
Further, in step S3, the relevant constraints include energy balance constraints and equipment operation constraints.
Further, in step S3, the energy balance constraints include an electrical balance constraint, a thermal balance constraint, and a cold balance constraint.
Further, in step S3, the expression of the electrical balance constraint is:
Figure BDA0003697108800000029
in the formula: p is load Representing the total electrical load of the park; p toIES Indicating electricity purchased in the park, P fromIES Indicating the electric quantity of the online in the park;
Figure BDA00036971088000000210
representing the output electrical power of the electrical energy-related device;
Figure BDA00036971088000000211
representing the input electrical power of the electrical energy-related device.
Further, in step S3, the expression of the thermal balance constraint is:
Figure BDA00036971088000000212
in the formula: h load Representing the total heat load of the park;
Figure BDA00036971088000000213
representing the output power of the thermal energy related device;
Figure BDA00036971088000000214
representing the input power of the thermal energy related device;
in step S3, the expression of the cold balance constraint is:
Figure BDA0003697108800000031
in the formula: q Load Indicating the total cooling load of the park, kW;
Figure BDA0003697108800000032
representing the output power of the thermal energy related device;
Figure BDA0003697108800000033
representing the input power of the thermal energy related device.
Further, in step S3, the expression of the operation constraint of the chiller in the equipment operation constraint is:
Figure BDA0003697108800000034
Figure BDA0003697108800000035
EnforceRun≤OperatState(t)≤1-EnforceOutage
in the formula:
Figure BDA0003697108800000036
representing the input electrical power of the chiller plant;
Figure BDA0003697108800000037
the output cold power of the water chilling unit equipment is represented; COP j Representing the refrigeration energy efficiency ratio of the water chilling unit equipment;
Figure BDA0003697108800000038
the installation capacity of the water chilling unit equipment is represented; enforceRun represents whether the equipment is forced to run or not, and is a variable of 0/1; the EnforceOutage indicates whether the equipment is overhauled, shut down and standby, and is a variable of 0/1; operatState represents the operational state of the device, where 1 represents run and 0 represents shut down.
Further, in the step S3, after optimization, more power is used when the clean energy is generated greatly, and less power is used when the traditional energy power generation occupation ratio is high.
Further, in step S3, the integrated energy park includes an energy storage device, the energy storage device includes at least one of an energy storage battery, a cold storage device and a heat storage device, when the dynamic carbon emission factor is low, the energy storage device stores energy, and when the dynamic carbon emission factor is high, the energy storage device releases energy. The high-low division of the dynamic carbon rejection factor can be adjusted according to actual data.
Further, based on the consideration of the park economy, the low-carbon operation optimization is performed when the importance of the low-carbon emission is higher than the park economy, or when the influence of the low-carbon operation on the economy is within an acceptable range. It can also be said that optimizing operation also takes into account park economics constraints.
The comprehensive energy park low-carbon operation optimization method based on the dynamic carbon emission factor has the beneficial effects that: acquiring a dynamic carbon emission factor, establishing a target function about the dynamic carbon emission factor, wherein the dynamic carbon emission factor can more accurately depict the total carbon emission amount of a comprehensive energy park and describe the phenomenon of low power carbon emission factor when a new energy generating set such as photovoltaic power generation, wind power generation and the like generates high power; the equipment operation and the energy storage equipment charge and discharge of comprehensive energy garden are more accurately guided according to the dynamic carbon emission factor, and when the electric power carbon emission factor is lower, the clean electric energy is used more to reduce the total carbon emission amount of the garden.
Drawings
FIG. 1 is a flow chart of a method for optimizing low-carbon operation of an integrated energy park based on dynamic carbon emission factors, according to an embodiment of the invention.
FIG. 2 is a schematic diagram of a method for calculating a dynamic carbon rejection factor according to an embodiment of the present invention.
Detailed Description
The technical solutions of the embodiments of the present invention will be explained and explained below with reference to the drawings of the embodiments of the present invention, but the embodiments described below are only preferred embodiments of the present invention, and are not all embodiments. Other embodiments obtained by persons skilled in the art without any inventive work based on the embodiments in the embodiment belong to the protection scope of the invention.
Referring to fig. 1 and fig. 2, the method for optimizing low-carbon operation of the integrated energy park based on the dynamic carbon emission factor according to the present invention comprises:
s1, acquiring dynamic carbon elimination factors of province areas or local cities of parks, which change along with time;
s2, establishing a target function related to the dynamic carbon emission factor by taking the minimum carbon emission of the garden as a target;
and S3, optimizing according to the relevant constraint conditions.
The comprehensive energy park low-carbon operation optimization method based on the dynamic carbon emission factors obtains the dynamic carbon emission factors, establishes an objective function about the dynamic carbon emission factors, and can more accurately depict the total carbon emission amount of the comprehensive energy park and describe the low-power carbon emission factor phenomenon when a new energy generating set such as a photovoltaic generator set and a wind power generator set generates high power; the equipment operation and the energy storage equipment charge and discharge of comprehensive energy garden are more accurately guided according to the dynamic carbon emission factor, and when the electric power carbon emission factor is lower, the clean electric energy is used more to reduce the total carbon emission amount of the garden.
Examples
Referring to fig. 1 and fig. 2, in an embodiment of the present invention, a method for optimizing low-carbon operation of an integrated energy park based on a dynamic carbon emission factor includes:
s1, acquiring dynamic carbon elimination factors of province areas or local cities of parks, which change along with time;
s2, establishing a target function related to the dynamic carbon emission factor by taking the minimum carbon emission of the garden as a target;
and S3, optimizing according to the relevant constraint conditions.
Wherein the relevant constraints include expected operating curves of the relevant devices.
In this embodiment, in step S1, the calculation formula of the provincial dynamic carbon elimination factor is as follows:
Figure BDA0003697108800000041
in the formula: f is dynamic carbon emission factor of provincial power supply side and has the unit of tCO 2 kWh; i is the energy type of the power-saving power supply side regulating unit, and comprises coal-fired power generation, gas power generation, hydroelectric power generation, photovoltaic power generation, wind power generation and nuclear power generation; w is an external electricity region; c is the type of the coal burner group; g is a gas engine set type; g is the generated energy with the unit of kWh,taking the value as the accumulated value of the gateway metering for 15 min; EF c Is the carbon emission factor of the coal-fired unit and has the unit of tCO 2 /kWh；EF g Is the carbon emission factor of the gas turbine set and has the unit of tCO 2 /kWh；EF w Average carbon emission factor of foreign electric regions; t represents time.
In this embodiment, after the provincial dynamic carbon emission factor is obtained, the carbon flow conditions of each local city are continuously refined according to the provincial power topological graph, and the dynamic carbon emission factor of the local city to which the integrated energy park belongs is obtained.
In this embodiment, in step S2, the objective function is:
Figure BDA0003697108800000051
in the formula: f is the carbon emission of the park;
Figure BDA0003697108800000052
carbon emission coefficient per unit of electrical energy purchased from the grid for IES, which dynamically changes over time, tco 2 /kWh；P FromGrid To purchase electricity from the grid, kWh;
Figure BDA0003697108800000053
carbon emission coefficient, tCO, of gas-fired plants 2 /kWh；
Figure BDA0003697108800000054
Carbon emission coefficient, tCO, for coal burning plants 2 /kWh；
Figure BDA0003697108800000055
The output power of the gas equipment is kW;
Figure BDA0003697108800000056
the output power of the coal-fired equipment is kW; t represents time; i is the device type; j is the equipment serial number; the carbon emission coefficient of each device is the property of the device, and has specificity.
Figure BDA0003697108800000057
Namely the provincial dynamic carbon elimination factor or the local city dynamic carbon elimination factor.
In this embodiment, in step S3, the relevant constraint conditions include an energy balance constraint and an equipment operation constraint.
In this embodiment, in step S3, the energy balance constraint includes an electrical balance constraint, a thermal balance constraint, and a cold balance constraint.
In this embodiment, in step S3, the expression of the electrical balance constraint is:
Figure BDA0003697108800000058
in the formula: p load Representing the total electrical load of the park; p toIES Indicating electricity purchased in the park, P fromIES Indicating the electric quantity of the online in the park;
Figure BDA0003697108800000059
representing the output electrical power of the electrical energy-related device;
Figure BDA00036971088000000510
representing the input electrical power of the electrical energy-related device.
In this embodiment, in step S3, the expression of the thermal balance constraint is:
Figure BDA00036971088000000511
in the formula: h load Representing the total heat load of the park;
Figure BDA00036971088000000512
representing the output power of the thermal energy related device;
Figure BDA00036971088000000513
representing the input power of the thermal energy related device;
in step S3, the expression of the cold balance constraint is:
Figure BDA00036971088000000514
in the formula: q Load Indicating the total cooling load of the park, kW;
Figure BDA00036971088000000515
representing the output power of the thermal energy related device;
Figure BDA00036971088000000516
representing the input power of the thermal energy related device.
In this embodiment, in step S3, an expression of the operation constraint of the chiller in the equipment operation constraint is as follows:
Figure BDA0003697108800000061
Figure BDA0003697108800000062
EnforceRun≤OperatState(t)≤1-EnforceOutage
in the formula:
Figure BDA0003697108800000063
representing the input electrical power of the chiller plant;
Figure BDA0003697108800000064
the output cold power of the water chilling unit equipment is represented; COP j Representing the refrigeration energy efficiency ratio of the water chilling unit equipment;
Figure BDA0003697108800000065
the installation capacity of the water chilling unit equipment is represented; the EnforceRun indicates whether the equipment is forced to run, and is a variable of 0/1; enforceOutage indicates whether the equipment is overhauled, shut down and standby, and is changed into 0/1An amount; operatState represents the operational state of the device, where 1 represents run and 0 represents shut down.
In this embodiment, in step S3, after optimization, the power is used more when the clean energy is generated greatly, and the power is used less when the power generation occupancy ratio of the conventional energy is high.
The basic equipment of the comprehensive energy park comprises a transformer, a photovoltaic device, a fan, an energy storage battery, a cold accumulation device, a heat accumulation device, a coal-fired unit, a gas unit, an electric boiler, a heat pump, an absorption refrigerator, a water chilling unit and the like. The electric energy related equipment comprises a transformer, a photovoltaic, a fan, an energy storage battery, a coal-fired unit, a gas unit, an electric boiler, a heat pump, a water chilling unit and the like; the heat energy related equipment comprises heat storage equipment, a coal-fired unit, a gas unit, an electric boiler, a heat pump, an absorption refrigerator and the like; the cold energy related equipment comprises cold storage equipment, a heat pump, an absorption refrigerator, a water chilling unit and the like. The electric energy related equipment and the heat energy related equipment are partially crossed, and the electric energy related equipment and the cold energy related equipment are partially crossed.
For a more specific description, a certain park is taken as an example, and the park is assumed to have only one type of electric equipment, one set of energy storage equipment and a transformer. The optimization result is as follows: when the factor is lower to electric power carbon row, energy storage battery begins to charge, when the factor is higher to electric power carbon row, energy storage equipment begins to discharge, and the electric quantity is purchased from the electric wire netting to the reduction of garden to satisfy consumer's electric load demand. Similarly, the park comprising various equipment types is optimized, and after optimization, the carbon emission factor of the electric power is lower when the clean energy is generated greatly, at the moment, the energy storage battery and the cold and heat storage equipment start to store energy, and the power supply of the park is mainly provided by external power purchase; when traditional energy power generation occupation ratio, electric power carbon emission factor is higher, and energy storage equipment output power increases this moment, reduces outsourcing electricity occupation ratio, satisfies the power consumption, the heat with, with cold demand in garden.
The comprehensive energy park low-carbon operation optimization method based on the dynamic carbon emission factor has the beneficial effects that: acquiring a dynamic carbon emission factor, and establishing an objective function about the dynamic carbon emission factor, wherein the dynamic carbon emission factor can more accurately depict the total carbon emission amount of a comprehensive energy park and describe the phenomenon of a low-power carbon emission factor when a new energy generating set such as photovoltaic power, wind power and the like generates high power; the equipment operation and the energy storage equipment charging and discharging of the comprehensive energy park are more accurately guided according to the dynamic carbon emission factor, outsourcing power is added to store energy for the energy storage equipment when the electric carbon emission factor is low, clean electric energy is used more to reduce the total carbon emission amount of the park, and the energy storage equipment discharges energy when the electric carbon emission factor is high to reduce the outsourcing power; park economics are taken into constraints to minimize negative effects on park economics.
While the invention has been described with reference to specific embodiments, it will be understood by those skilled in the art that the invention is not limited thereto but is intended to be limited only by the foregoing description. Any modification which does not depart from the functional and structural principles of the invention is intended to be included within the scope of the following claims.

Claims (10)

1. A comprehensive energy park low-carbon operation optimization method based on dynamic carbon emission factors is characterized by comprising the following steps: the comprehensive energy park low-carbon operation optimization method based on the dynamic carbon emission factor comprises the following steps:
s1, acquiring dynamic carbon emission factors of provinces or local cities of parks along with time change;
s2, establishing a target function related to the dynamic carbon emission factor by taking the minimum carbon emission of the park as a target;
and S3, optimizing according to the relevant constraint conditions.
2. The method for optimizing the low-carbon operation of the integrated energy park based on the dynamic carbon emission factor according to claim 1, wherein the method comprises the following steps: in step S1, the calculation formula of the provincial domain dynamic carbon emission factor is as follows:
Figure FDA0003697108790000011
in the formula: f is dynamic carbon elimination factor of power supply side in provincial regionIn the position tCO 2 /kWh; i is the energy type of the power-saving power supply side regulating unit, and comprises coal-fired power generation, gas power generation, hydroelectric power generation, photovoltaic power generation, wind power generation and nuclear power generation; w is an external electricity region; c is the type of the coal burner group; g is a gas engine set type; g is generated energy, the unit is kWh, and the value is the accumulated value of the gateway metering for 15 min; EF c Is the carbon emission factor of the coal-fired unit and has the unit of tCO 2 /kWh；EF g Is the carbon emission factor of the gas turbine set and has the unit of tCO 2 /kWh；EF w Average carbon emission factor of foreign electric regions; t represents time.
3. The method for optimizing the low-carbon operation of the integrated energy park based on the dynamic carbon emission factor as claimed in claim 2, wherein the method comprises the following steps: and after the provincial region dynamic carbon emission factor is obtained, according to the provincial power topological graph, the carbon flow conditions of various cities are continuously refined, and the dynamic carbon emission factor of the city to which the comprehensive energy park belongs is obtained.
4. The comprehensive energy park low-carbon operation optimization method based on the dynamic carbon emission factor according to claim 2 or 3, characterized in that: in step S2, the objective function is:
Figure FDA0003697108790000012
in the formula: f is the carbon emission of the park;
Figure FDA0003697108790000013
carbon emission coefficient per unit of electrical energy purchased from the grid for IES, which dynamically changes over time, tco 2 /kWh；P FromGrid To purchase electricity from the grid, kWh;
Figure FDA0003697108790000014
carbon emission coefficient, tCO, of gas-fired plants 2 /kWh；
Figure FDA0003697108790000015
Carbon emission coefficient, tCO, for coal burning plants 2 /kWh；
Figure FDA0003697108790000016
kW is the output power of the gas equipment;
Figure FDA0003697108790000017
the output power of the coal-fired equipment is kW; t represents time; i is the device type; j is the equipment serial number; the carbon emission coefficient of each device is the property of the device, and has specificity.
5. The method for optimizing the low-carbon operation of the integrated energy park based on the dynamic carbon emission factor according to claim 4, wherein the method comprises the following steps: in step S3, the relevant constraint conditions include an energy balance constraint and a device operation constraint, and the energy balance constraint includes an electrical balance constraint, a thermal balance constraint and a cold balance constraint.
6. The comprehensive energy park low-carbon operation optimization method based on the dynamic carbon emission factor according to claim 5, characterized in that: in step S3, the electrical balance constraint expression is:
Figure FDA0003697108790000021
in the formula: p load Representing the total electrical load of the park; p is toIES Indicating electricity purchased in the park, P fromIES Indicating the electric quantity of the online in the park;
Figure FDA0003697108790000022
represents the output electric power of the electric energy related device;
Figure FDA0003697108790000023
representing the input electrical power of the electrical energy-related device.
7. The comprehensive energy park low-carbon operation optimization method based on the dynamic carbon emission factor according to claim 5, characterized in that: in step S3, the expression of the thermal balance constraint is:
Figure FDA0003697108790000024
in the formula: h load Representing the total heat load of the park;
Figure FDA0003697108790000025
representing the output power of the thermal energy related device;
Figure FDA0003697108790000026
representing the input power of the thermal energy related device;
in step S3, the expression of the cold balance constraint is:
Figure FDA0003697108790000027
in the formula: q Load Indicating the total cooling load of the park, kW;
Figure FDA0003697108790000028
representing the output power of the thermal energy related device;
Figure FDA0003697108790000029
representing the input power of the thermal energy related device.
8. The comprehensive energy park low-carbon operation optimization method based on the dynamic carbon emission factor according to claim 5, characterized in that: in step S3, the expression of the operation constraint of the chiller in the equipment operation constraint is:
Figure FDA00036971087900000210
Figure FDA00036971087900000211
EnforceRun≤OperatState(t)≤1-EnforceOutage
in the formula:
Figure FDA00036971087900000212
representing the input electrical power of the chiller plant;
Figure FDA00036971087900000213
the output cold power of the water chilling unit equipment is represented; COP j Representing the refrigeration energy efficiency ratio of the water chilling unit equipment;
Figure FDA00036971087900000214
the installation capacity of the water chilling unit equipment is represented; the EnforceRun indicates whether the equipment is forced to run, and is a variable of 0/1; the EnforceOutage indicates whether the equipment is overhauled, shut down and standby, and is a variable of 0/1; operatState represents the operational state of the device, where 1 represents run and 0 represents shut down.
9. The comprehensive energy park low-carbon operation optimization method based on the dynamic carbon emission factor as claimed in claim 1, characterized in that: in the step S3, after optimization, the power is used more when the clean energy is generated greatly, and the power is used less when the traditional energy power generation occupation ratio is high; in the step S3, the comprehensive energy park comprises energy storage equipment, the energy storage equipment comprises at least one of an energy storage battery, cold accumulation equipment and heat accumulation equipment, when the dynamic carbon emission factor is low, the energy storage equipment stores energy, and when the dynamic carbon emission factor is high, the energy storage equipment discharges energy.
10. The method for optimizing the low-carbon operation of the integrated energy park based on the dynamic carbon emission factor according to claim 1, wherein the method comprises the following steps: and step S3, combining the economy of the park, and performing low-carbon operation optimization when the importance of low-carbon emission is higher than the economy of the park, or performing low-carbon operation when the influence of the low-carbon operation on the economy is within an acceptable range.
CN202210677812.4A 2022-06-16 2022-06-16 Comprehensive energy park low-carbon operation optimization method based on dynamic carbon emission factor Pending CN115239538A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210677812.4A CN115239538A (en) 2022-06-16 2022-06-16 Comprehensive energy park low-carbon operation optimization method based on dynamic carbon emission factor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210677812.4A CN115239538A (en) 2022-06-16 2022-06-16 Comprehensive energy park low-carbon operation optimization method based on dynamic carbon emission factor

Publications (1)

Publication Number Publication Date
CN115239538A true CN115239538A (en) 2022-10-25

Family

ID=83669507

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210677812.4A Pending CN115239538A (en) 2022-06-16 2022-06-16 Comprehensive energy park low-carbon operation optimization method based on dynamic carbon emission factor

Country Status (1)

Country Link
CN (1) CN115239538A (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115545551A (en) * 2022-11-04 2022-12-30 北京如实智慧电力科技有限公司 Photovoltaic online carbon asset checking system and calculation method
CN117171949A (en) * 2023-07-18 2023-12-05 南京电力设计研究院有限公司 Method for deducting carbon emission situation of digital park
CN117745109A (en) * 2024-02-21 2024-03-22 新奥数能科技有限公司 Low-carbon optimized energy supply mode determining method and system based on multi-energy complementation

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115545551A (en) * 2022-11-04 2022-12-30 北京如实智慧电力科技有限公司 Photovoltaic online carbon asset checking system and calculation method
CN117171949A (en) * 2023-07-18 2023-12-05 南京电力设计研究院有限公司 Method for deducting carbon emission situation of digital park
CN117171949B (en) * 2023-07-18 2024-04-05 南京电力设计研究院有限公司 Method for deducting carbon emission situation of digital park
CN117745109A (en) * 2024-02-21 2024-03-22 新奥数能科技有限公司 Low-carbon optimized energy supply mode determining method and system based on multi-energy complementation

Similar Documents

Publication Publication Date Title
CN115239538A (en) Comprehensive energy park low-carbon operation optimization method based on dynamic carbon emission factor
CN103580063B (en) A kind of method of large-scale grid connection wind-powered electricity generation of dissolving based on demanding party's response
CN104616208B (en) A kind of supply of cooling, heating and electrical powers type micro-capacitance sensor operation method based on Model Predictive Control
CN109474025B (en) Optimized dispatching model of park level comprehensive energy system
CN108808663B (en) Industrial user heat demand response method based on multi-energy complementation
CN113762708A (en) Park level comprehensive energy system planning method considering multi-target cooperation
Chen et al. Optimal planning and operation of CCHP system considering renewable energy integration and seawater desalination
CN114626721A (en) Agricultural industrial park near-zero carbon implementation method based on time-shifting load scheduling
CN114580746A (en) Comprehensive energy station composite energy storage configuration optimization method based on low-carbon economic benefit quantification
CN108197412B (en) Multi-energy coupling energy management system and optimization method
CN114255137A (en) Low-carbon comprehensive energy system optimization planning method and system considering clean energy
CN114396822B (en) Energy comprehensive utilization configuration and operation method
Huan et al. Retracted: Research on Optimized Operation of Electrothermal Combined System for Enhancing Wind Power Consumption
CN114330835A (en) Optimal configuration method of electricity/heat hybrid energy storage system in comprehensive energy microgrid
CN112465228A (en) User-side comprehensive energy system optimal configuration method
Shi et al. Research on Market Mechanism of Thermal Power Flexible Transformation to Promote Clean Energy Consumption
CN113807746B (en) Comprehensive operation optimization method of combined cooling heating power system
Chen et al. Flexible transformation cost analysis of Thermal power in Jiangxi Province for depth peak-shaving
CN112907147B (en) Distributed energy system operation method and system based on optimal working condition point
Xie et al. Planning Method for Power-Gas-Heat-Steam Integrated Energy System in Industrial Park Considering Annual Comprehensive Cost
Zhao et al. Optimal Capacity Allocation of Integrated Energy System Considering Multiple Uncertainties
Gong et al. Combined Cooling Heating and Power System Design and Capacity Configuration taking into account Solar Photovoltaic
Zhao et al. CCHP Capacity Optimization with User Demand Characteristics
Liu et al. Research on capacity allocation of community integrated energy system planning based on multi-objective optimization
Chen et al. Combined Heat and Power Optimal Scheduling of Virtual Power Plant Considering Carbon Capture and Demand Response

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