CN113435095A - Method and system for optimizing scheduling of comprehensive energy system - Google Patents

Method and system for optimizing scheduling of comprehensive energy system Download PDF

Info

Publication number
CN113435095A
CN113435095A CN202110655638.9A CN202110655638A CN113435095A CN 113435095 A CN113435095 A CN 113435095A CN 202110655638 A CN202110655638 A CN 202110655638A CN 113435095 A CN113435095 A CN 113435095A
Authority
CN
China
Prior art keywords
power
heat
energy
energy system
objective function
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
CN202110655638.9A
Other languages
Chinese (zh)
Inventor
牛耕
刘晓娟
吴鸣
季宇
李蕊
寇凌峰
张智慧
王文博
孙丽敬
张海
屈小云
丁保迪
张颖
郑楠
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Xinjiang Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
State Grid Xinjiang Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI, State Grid Xinjiang Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202110655638.9A priority Critical patent/CN113435095A/en
Publication of CN113435095A publication Critical patent/CN113435095A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/25—Design optimisation, verification or simulation using particle-based methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/12—Computing arrangements based on biological models using genetic models
    • G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
    • 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
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06315—Needs-based resource requirements planning or analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • 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
    • G06Q30/00—Commerce
    • G06Q30/02—Marketing; Price estimation or determination; Fundraising
    • G06Q30/0201—Market modelling; Market analysis; Collecting market data
    • 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
    • G06Q30/00—Commerce
    • G06Q30/02—Marketing; Price estimation or determination; Fundraising
    • G06Q30/0201—Market modelling; Market analysis; Collecting market data
    • G06Q30/0206—Price or cost determination based on market factors
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • 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
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS 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
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • 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
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS 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/00—Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
    • Y04S50/14—Marketing, i.e. market research and analysis, surveying, promotions, advertising, buyer profiling, customer management or rewards

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Human Resources & Organizations (AREA)
  • General Physics & Mathematics (AREA)
  • Development Economics (AREA)
  • Economics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Accounting & Taxation (AREA)
  • Finance (AREA)
  • Marketing (AREA)
  • General Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Evolutionary Computation (AREA)
  • Game Theory and Decision Science (AREA)
  • Biophysics (AREA)
  • Data Mining & Analysis (AREA)
  • Tourism & Hospitality (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Educational Administration (AREA)
  • Mathematical Physics (AREA)
  • Biomedical Technology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computational Linguistics (AREA)
  • Geometry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Evolutionary Biology (AREA)
  • Computer Hardware Design (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Primary Health Care (AREA)
  • Genetics & Genomics (AREA)

Abstract

The invention discloses a method and a system for optimizing and scheduling of an integrated energy system, wherein the method comprises the following steps: establishing an objective function for optimizing scheduling of an integrated energy system, the objective function comprising: an objective function based on operating cost, an objective function based on environmental cost, and an objective function based on reliability; determining constraints for optimal scheduling of the integrated energy system, the constraints comprising: the system comprises a cold-heat-electricity load balance constraint, a power generation device constraint, a refrigeration and heating device constraint, a device climbing rate constraint, an electric energy transmission power constraint and an energy storage constraint between a comprehensive system and a power grid; and acquiring the operation parameters of the comprehensive energy system, solving the objective function and the constraint condition based on the operation parameters of the comprehensive energy system, and generating an operation scheme for optimizing and scheduling the comprehensive energy system.

Description

Method and system for optimizing scheduling of comprehensive energy system
Technical Field
The invention relates to the technical field of comprehensive energy system optimization, in particular to a method and a system for optimizing and scheduling a comprehensive energy system.
Background
Energy is the basis for human survival and development, and how to fully utilize renewable energy and improve the utilization efficiency of comprehensive energy is a hot content of research in the related field of comprehensive energy. How to optimally control the operation of the electricity/gas/heat and other comprehensive energy systems is an important problem through coupling complementation and cascade utilization of various energy forms, the impact of distributed energy fluctuation on a power grid is reduced, the development and application of renewable energy sources are promoted, the shortage of fossil energy is relieved, and the environmental pollution is reduced. From the perspective of energy utilization, the comprehensive energy system coupled by multiple energy sources has correlation and complementarity on different time scales, and can store and supply energy on multiple time scales. Therefore, the optimal scheduling of the integrated energy system becomes a fundamental problem in the generation and utilization of energy by the integrated energy system, and intensive research is required.
In the prior art, scholars experts at home and abroad carry out a great deal of research on the aspect of optimizing and scheduling of the comprehensive energy system, but most of the existing related research is directed at the comprehensive energy system with a cogeneration system or a combined cooling heating and power system. However, there are still few studies on the optimized scheduling aspect of the integrated energy systems of electricity, gas and heat that comprehensively consider the renewable energy power generation equipment such as wind and light, the micro-combustion unit, the electric heat energy storage, the cooling, heating and power load requirements and the energy interaction with the power grid side, and at the same time, such integrated energy systems more conform to the actual situation of most integrated energy systems in China.
Therefore, intensive research needs to be carried out on an optimized scheduling model suitable for the electric, gas and thermal comprehensive energy system.
Disclosure of Invention
The technical scheme of the invention provides a method and a system for optimizing and scheduling a comprehensive energy system, which aim to solve the problem of how to optimize and schedule the comprehensive energy system.
In order to solve the above problems, the present invention provides a method for optimal scheduling of an integrated energy system, the method comprising:
determining an objective function of optimal scheduling of the comprehensive energy system and a constraint condition; and acquiring the operation parameters of the comprehensive energy system, solving the objective function and the constraint condition based on the operation parameters, and generating an operation scheme for optimizing and scheduling the comprehensive energy system.
Preferably, the integrated energy system comprises at least one of: electrical energy, gas energy, and thermal energy.
Preferably, the solving the objective function and the constraint condition based on the operation parameters includes:
and solving the objective function and the constraint condition through one algorithm of a genetic algorithm, a particle swarm algorithm or a differential evolution algorithm.
Preferably, the establishing an objective function of the optimized scheduling of the integrated energy system includes: an objective function based on operating cost, an objective function based on environmental cost, and an objective function based on reliability, wherein:
the objective function based on the running cost is:
F1＝minCop＝Cop-E+Cop-H+Cop-NG
Figure BDA0003113569280000021
Figure BDA0003113569280000022
Figure BDA0003113569280000023
in the formula:
F1an integrated energy system operating cost objective function;
Copthe operating cost of the comprehensive energy system is reduced;
Cop-Ethe operating cost of the power supply module;
Cop-Hthe operating cost of the heat supply module;
Cop-NGthe acquisition cost of the natural gas consumed in the operation process of the comprehensive energy system;
Cithe operation cost of the distributed power generation equipment i under unit power;
Pi(t) is the output power of the distributed generation equipment i at the moment t;
n is the total number of distributed power generation equipment;
t is the total running time of the equipment;
Cgrid(t) is the real-time electricity price at time t;
Pgridsell(t) the power of the comprehensive energy system selling electricity to the power grid at the moment t;
Pgridbuy(t) the power purchased by the comprehensive energy system to the power grid at the moment t;
Celec-gridthe fixed operation maintenance cost of the electric energy supply module network in a unit time interval can be generally calculated by taking one hour as the unit time interval;
CHB-girdfor the fixed operation and maintenance cost of the heat supply module network in a unit time interval, one can be usedCalculating by taking the hour as a unit time interval;
αh，jthe operating cost of the heat supply equipment j under unit power;
Ph,j(t) is the thermal output power of the heat supply device j at time t;
r is the total number of heat supply devices;
aNG，sthe cost of gas usage per unit power for a natural gas consuming device s;
PNG，s(t) the output power of the plant s consuming natural gas at time t;
g is the total number of devices consuming natural gas;
the objective function based on the environmental cost is:
Figure BDA0003113569280000031
in the formula:
F2an environmental cost objective function of the integrated energy system;
CEenvironmental cost for the integrated energy system;
Pk(t) power of discharge source k at time t;
ωk,jis the emission coefficient of pollutant j from emission source k;
δE,jis a pollutant emission level;
ρEC-ppenalizing costs for pollutant emissions;
m is the total number of types of emissions;
w is the total number of emission sources;
the reliability-based objective function is:
Figure BDA0003113569280000041
in the formula:
F3a power supply reliability objective function for the comprehensive energy system;
the LESP is the energy supply loss rate of the comprehensive energy system;
t is the total operation time of the comprehensive energy system;
Edemand,tthe total energy demand in the comprehensive energy system at the time t;
Esupply,tand supplying energy demand in the comprehensive energy system at the moment t.
Preferably, the constraint condition for determining the optimal scheduling of the integrated energy system comprises: the energy storage system comprises a cold-heat-electricity load balance constraint, a power generation device constraint, a refrigeration and heating device constraint, a device climbing rate constraint, an electric energy transmission power constraint and an energy storage constraint between the comprehensive energy system and a power grid, wherein:
the cooling, heating and power load balance constraint is calculated according to the following formula:
Ppv(t)+Pwt(t)+Pmt(t)+Pgridbuy(t)+Pbat,d(t)＝Pec(t)+Pgridsell(t)+Pbat,c(t)+Pload(t)
in the formula: ppv(t)、Pwt(t)、Pmt(t)、Pgridbuy(t)、Pbat,d(t) photovoltaic power generation, wind power generation, micro gas turbine power generation, power purchase to a power grid and electric energy storage and discharge capacity are respectively carried out; pec(t)、Pgridsell(t)、Pbat,c(t)、Pload(t) the power consumption of the electric refrigerator, the power selling to the power grid, the charging amount of the electric energy storage and the demand of the electric load are respectively;
the output of the heating equipment in the comprehensive energy system is the same as the heat load demand, and the calculation is carried out according to the following formula:
Qrec+Qgb+Qhs,d＝Qh,lond+Qc+Qhs,c
in the formula: qrec、Qgb、Qhs,dThe heat output by the waste heat boiler of the micro-gas turbine, the heat emitted by the gas boiler and the heat emitted by the heat energy storage device are respectively; qh,lond、Qc、Qhs,cThe heat load demand, the heat consumed by the absorption refrigerator and the heat storage quantity of the heat storage energy are respectively;
the output of the refrigeration equipment in the comprehensive energy system is the same as the cold load demand, and the calculation is carried out according to the following formula:
Qac+Qec＝Qc,lond
in the formula: qac、QecRespectively the cold energy output by the absorption refrigerator and the cold energy generated by the electric refrigerator; qc,loadIs the cooling load demand;
the power plant constraints are calculated according to the following formula:
Figure BDA0003113569280000051
Figure BDA0003113569280000052
Figure BDA0003113569280000053
in the formula:
Figure BDA0003113569280000054
the maximum values of the output power of the wind power, the photovoltaic and the micro-combustion engine are respectively;
the constraints of the refrigerating and heating equipment are calculated according to the following formula:
Figure BDA0003113569280000055
Figure BDA0003113569280000056
Figure BDA0003113569280000057
in the formula: qgb、Qac、QecAre respectively gas-fired boilersOutput power of the furnace, the absorption refrigerator, the electric refrigerator;
Figure BDA0003113569280000058
the maximum output power values of the gas boiler, the absorption refrigerator and the electric refrigerator are respectively;
the equipment ramp rate constraint is calculated according to the following formula:
-Ri,down＜Pi(t)-Pi(t-1)＜Ri,up
in the formula: pi(t)、Pi(t-1) power of the equipment at time t and (t-1) respectively; ri,down、Ri,upMaximum downward and maximum upward climbing rates;
the electric energy transmission power constraint between the comprehensive energy system and the power grid is calculated according to the following formula:
Figure BDA0003113569280000059
Figure BDA00031135692800000510
in the formula:
Figure BDA00031135692800000511
the upper limit values of electricity selling and electricity purchasing between the comprehensive energy system and the power grid are respectively set;
the energy storage constraints include electrical energy storage constraints and thermal energy storage constraints, wherein the electrical energy storage constraints are calculated according to the following formula:
Figure BDA00031135692800000512
0≤Pbat,c≤Srated·γbat,c
0≤Pbat,d≤Srated·γbat,d
in the formula:
Figure BDA00031135692800000513
respectively storing minimum and maximum energy of electric energy; sratedRated capacity for electrical energy storage; gamma raybat,c、γbat,dThe maximum charging rate and the maximum discharging rate of the electric energy storage are respectively; pbat，d、Pbat，cDischarge power and charging power which are electric energy storage respectively; the charging and discharging of the electric energy storage belong to a coupling relation, and the electric energy storage and the discharging can not be carried out at the same time;
Pbat，c·Pbat，d＝0
the thermal energy storage constraint is calculated according to the following formula:
Figure BDA0003113569280000061
Figure BDA0003113569280000062
Figure BDA0003113569280000063
in the formula:
Figure BDA0003113569280000064
the minimum and maximum heat storage quantity of the heat storage device are respectively;
Figure BDA0003113569280000065
the maximum heat storage power and the maximum heat release power of the heat storage device are respectively; qhs，c、Qhs，dThe heat storage power and the heat release power of the heat storage device are respectively; the heat storage device can not store heat and release heat at the same time, and the following constraint is added:
Qhs，c·Qhs，d＝0。
in accordance with another aspect of the present invention, there is provided a system for optimized scheduling of an integrated energy system, the system comprising:
the system comprises an establishing unit, a scheduling unit and a scheduling unit, wherein the establishing unit is used for determining an objective function of optimal scheduling of the comprehensive energy system and a constraint condition; and the result unit is used for acquiring the operation parameters of the comprehensive energy system, solving the objective function and the constraint condition based on the operation parameters and generating an operation scheme for optimizing and scheduling the comprehensive energy system.
Preferably, the integrated energy system comprises at least one of: electrical energy, gas energy, and thermal energy.
Preferably, the result unit is specifically configured to solve the objective function and the constraint condition through one of a genetic algorithm, a particle swarm algorithm, or a differential evolution algorithm.
Preferably, the establishing unit is configured to establish an objective function of the optimized scheduling of the integrated energy system, where the objective function includes: an objective function based on operating cost, an objective function based on environmental cost, and an objective function based on reliability, wherein:
the objective function based on the running cost is:
F1＝minCop＝Cop-E+Cop-H+Cop-NG
Figure BDA0003113569280000071
Figure BDA0003113569280000072
Figure BDA0003113569280000073
in the formula:
F1an integrated energy system operating cost objective function;
Copthe operating cost of the comprehensive energy system is reduced;
Cop-Ethe operating cost of the power supply module;
Cop-Hthe operating cost of the heat supply module;
Cop-NGthe acquisition cost of the natural gas consumed in the operation process of the comprehensive energy system;
Cithe operation cost of the distributed power generation equipment i under unit power;
Pi(t) is the output power of the distributed generation equipment i at the moment t;
n is the total number of distributed power generation equipment;
t is the total running time of the equipment;
Cgrid(t) is the real-time electricity price at time t;
Pgridsell(t) the power of the comprehensive energy system selling electricity to the power grid at the moment t;
Pgridbuy(t) the power purchased by the comprehensive energy system to the power grid at the moment t;
Celec-gridthe fixed operation maintenance cost of the electric energy supply module network in a unit time interval can be generally calculated by taking one hour as the unit time interval;
CHB-girdthe fixed operation maintenance cost of the heat supply module network in a unit time interval can be generally calculated by taking one hour as the unit time interval;
ah，jthe operating cost of the heat supply equipment j under unit power;
Ph,j(t) is the thermal output power of the heat supply device j at time t;
r is the total number of heat supply devices;
aNG，sthe cost of gas usage per unit power for a natural gas consuming device s;
PNG，s(t) the output power of the plant s consuming natural gas at time t;
g is the total number of devices consuming natural gas;
the objective function based on the environmental cost is:
Figure BDA0003113569280000081
in the formula:
F2a cost objective function for the integrated system environment;
CEfor integrated system environmental costs;
Pk(t) power of discharge source k at time t;
ωk,jis the emission coefficient of pollutant j from emission source k;
δE,jis a pollutant emission level;
ρEC-ppenalizing costs for pollutant emissions;
m is the total number of types of emissions;
w is the total number of emission sources;
the reliability-based objective function is:
Figure BDA0003113569280000082
in the formula:
F3a power supply reliability objective function for the comprehensive energy system;
the LESP is the energy supply loss rate of the comprehensive energy system;
t is the total operation time of the comprehensive energy system;
Edemand,tthe total energy demand in the comprehensive energy system at the time t;
Esupply,tand supplying energy demand in the comprehensive energy system at the moment t.
Preferably, the determining unit is configured to determine a constraint condition for optimizing scheduling of the integrated energy system, where the constraint condition includes: the energy storage system comprises a cold-heat-electricity load balance constraint, a power generation device constraint, a refrigeration and heating device constraint, a device climbing rate constraint, an electric energy transmission power constraint and an energy storage constraint between the comprehensive energy system and a power grid, wherein:
the cooling, heating and power load balance constraint is calculated according to the following formula:
Ppv(t)+Pwt(t)+Pmt(t)+Pgridbuy(t)+Pbat,d(t)＝Pec(t)+Pgridsell(t)+Pbat,c(t)+Pload(t)
in the formula: ppv(t)、Pwt(t)、Pmt(t)、Pgridbuy(t)、Pbat,d(t) photovoltaic power generation, wind power generation, micro gas turbine power generation, power purchase to a power grid and electric energy storage and discharge capacity are respectively carried out; pec(t)、Pgridsell(t)、Pbat,c(t)、Pload(t) the power consumption of the electric refrigerator, the power selling to the power grid, the charging amount of the electric energy storage and the demand of the electric load are respectively;
the output of the heating equipment in the comprehensive energy system is the same as the heat load demand, and the calculation is carried out according to the following formula:
Qrec+Qgb+Qhs,d＝Qh,lond+Qc+Qhs,c
in the formula: qrec、Qgb、Qhs,dThe heat output by the waste heat boiler of the micro-gas turbine, the heat emitted by the gas boiler and the heat emitted by the heat energy storage device are respectively; qh,lond、Qc、Qhs,cThe heat load demand, the heat consumed by the absorption refrigerator and the heat storage quantity of the heat storage energy are respectively;
the output of the refrigeration equipment in the comprehensive energy system is the same as the cold load demand, and the calculation is carried out according to the following formula:
Qac+Qec＝Qc,lond
in the formula: qac、QecRespectively the cold energy output by the absorption refrigerator and the cold energy generated by the electric refrigerator; qc,loadIs the cooling load demand;
the power plant constraints are calculated according to the following formula:
Figure BDA0003113569280000091
Figure BDA0003113569280000092
Figure BDA0003113569280000093
in the formula:
Figure BDA0003113569280000094
the maximum values of the output power of the wind power, the photovoltaic and the micro-combustion engine are respectively;
the constraints of the refrigerating and heating equipment are calculated according to the following formula:
Figure BDA0003113569280000095
Figure BDA0003113569280000096
Figure BDA0003113569280000097
in the formula: qgb、Qac、QecThe output powers of the gas boiler, the absorption refrigerator and the electric refrigerator are respectively;
Figure BDA0003113569280000098
the maximum output power values of the gas boiler, the absorption refrigerator and the electric refrigerator are respectively;
the equipment ramp rate constraint is calculated according to the following formula:
-Ri,down＜Pi(t)-Pi(t-1)＜Ri,up
in the formula: pi(t)、Pi(t-1) power of the equipment at time t and (t-1) respectively; ri,down、Ri,upMaximum downward and maximum upward climbing rates;
the electric energy transmission power constraint between the comprehensive energy system and the power grid is calculated according to the following formula:
Figure BDA0003113569280000101
Figure BDA0003113569280000102
in the formula:
Figure BDA0003113569280000103
the upper limit values of electricity selling and electricity purchasing between the comprehensive energy system and the power grid are respectively set;
the energy storage constraints include electrical energy storage constraints and thermal energy storage constraints, wherein the electrical energy storage constraints are calculated according to the following formula:
Figure BDA0003113569280000104
0≤Pbat,c≤Srated·γbat,c
0≤Pbat,d≤Srated·γbat,d
in the formula:
Figure BDA0003113569280000105
respectively storing minimum and maximum energy of electric energy; sratedRated capacity for electrical energy storage; gamma raybat,c、γbat,dThe maximum charging rate and the maximum discharging rate of the electric energy storage are respectively; pbat，d、Pbat，cDischarge power and charging power which are electric energy storage respectively; the charging and discharging of the electric energy storage belong to a coupling relation, and the electric energy storage and the discharging can not be carried out at the same time;
Pbat，c·Pbat，d＝0
the thermal energy storage constraint is calculated according to the following formula:
Figure BDA0003113569280000106
Figure BDA0003113569280000107
Figure BDA0003113569280000108
in the formula:
Figure BDA0003113569280000109
the minimum and maximum heat storage quantity of the heat storage device are respectively;
Figure BDA00031135692800001010
the maximum heat storage power and the maximum heat release power of the heat storage device are respectively; qhs.c、Qhs.dThe heat storage power and the heat release power of the heat storage device are respectively; the heat storage device can not store heat and release heat at the same time, and the following constraint is added:
Qhs，c·Qhs，d＝0。
the technical scheme of the invention provides a method and a system for optimizing and scheduling a comprehensive energy system, wherein the method comprises the following steps: establishing an objective function of optimized scheduling of the integrated energy system, wherein the objective function comprises the following steps: an objective function based on operating cost, an objective function based on environmental cost, and an objective function based on reliability; determining constraint conditions of optimal scheduling of the integrated energy system, wherein the constraint conditions comprise: the system comprises a cold-heat-electricity load balance constraint, a power generation device constraint, a refrigeration and heating device constraint, a device climbing rate constraint, an electric energy transmission power constraint and an energy storage constraint between a comprehensive system and a power grid; and acquiring the operation parameters of the comprehensive energy system, solving the objective function and the constraint condition based on the operation parameters of the comprehensive energy system, and generating an operation scheme for optimizing and scheduling the comprehensive energy system. The technical scheme of the invention provides an optimal scheduling method and system suitable for an electric, gas and thermal comprehensive energy system, which can coordinate, optimize and regulate multiple energy flows in the electric, gas and thermal comprehensive distributed energy system, and realize economic, environment-friendly and reliable operation of the whole system.
Drawings
A more complete understanding of exemplary embodiments of the present invention may be had by reference to the following drawings in which:
FIG. 1 is a flow diagram of a method for optimal scheduling of an integrated energy system in accordance with a preferred embodiment of the present invention;
FIG. 2 is a flow chart of a method for optimal regulation of an integrated energy system according to a preferred embodiment of the present invention; and
fig. 3 is a system configuration diagram for optimal scheduling of an integrated energy system according to a preferred embodiment of the present invention.
Detailed Description
The exemplary embodiments of the present invention will now be described with reference to the accompanying drawings, however, the present invention may be embodied in many different forms and is not limited to the embodiments described herein, which are provided for complete and complete disclosure of the present invention and to fully convey the scope of the present invention to those skilled in the art. The terminology used in the exemplary embodiments illustrated in the accompanying drawings is not intended to be limiting of the invention. In the drawings, the same units/elements are denoted by the same reference numerals.
Unless otherwise defined, terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Further, it will be understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the relevant art and will not be interpreted in an idealized or overly formal sense.
Fig. 1 is a flowchart of a method for optimal scheduling of an integrated energy system according to a preferred embodiment of the present invention. The invention provides an optimized dispatching model suitable for an electric, gas and thermal comprehensive energy system, in particular to a complex electric, gas and thermal comprehensive energy system containing equipment elements such as photovoltaic equipment, a fan, a micro-gas turbine, a gas boiler, electric energy storage, thermal energy storage, refrigerating/heating equipment, cooling, heating and thermal power loads and the like, and can realize coordinated optimized dispatching control of various energy flows in the comprehensive energy system. The invention provides an optimized dispatching model suitable for an electric, gas and thermal comprehensive energy system, which mainly comprises an objective function and constraint conditions and is a complex nonlinear multi-objective optimization model.
As shown in fig. 1, the present invention provides a method for optimized scheduling of an integrated energy system, the method comprising:
step 101: establishing an objective function of optimized scheduling of the integrated energy system, wherein the objective function comprises the following steps: an objective function based on operating cost, an objective function based on environmental cost, and an objective function based on reliability; preferably, the integrated energy system comprises: electrical energy, gas energy, and thermal energy.
Preferably, an objective function of the optimized scheduling of the integrated energy system is established, and the objective function includes: an objective function based on operating cost, an objective function based on environmental cost, and an objective function based on reliability, wherein:
the objective function based on the running cost is:
F1＝minCop＝Cop-E+Cop-H+Cop-NG
Figure BDA0003113569280000121
Figure BDA0003113569280000122
Figure BDA0003113569280000123
in the formula:
F1an integrated energy system operating cost objective function;
Copthe operating cost of the comprehensive energy system is reduced;
Cop-Ethe operating cost of the power supply module;
Cop-Hthe operating cost of the heat supply module;
Cop-NGthe acquisition cost of the natural gas consumed in the operation process of the comprehensive energy system;
Cithe operation cost of the distributed power generation equipment i under unit power;
Pi(t) is the output power of the distributed generation equipment i at the moment t;
n is the total number of distributed power generation equipment;
t is the total running time of the equipment;
Cgrid(t) is the real-time electricity price at time t;
Pgridsell(t) the power of the comprehensive energy system selling electricity to the power grid at the moment t;
Pgridbuy(t) the power purchased by the comprehensive energy system to the power grid at the moment t;
Celec-gridthe fixed operation maintenance cost of the electric energy supply module network in a unit time interval can be generally calculated by taking one hour as the unit time interval;
CHB-girdthe fixed operation maintenance cost of the heat supply module network in a unit time interval can be generally calculated by taking one hour as the unit time interval;
ah，jthe operating cost of the heat supply equipment j under unit power;
Ph,j(t) is the thermal output power of the heat supply device j at time t;
r is the total number of heat supply devices;
αNG，sthe cost of gas usage per unit power for a natural gas consuming device s;
PNG，s(t) the output power of the plant s consuming natural gas at time t;
g is the total number of devices consuming natural gas;
the objective function based on the environmental cost is:
Figure BDA0003113569280000131
in the formula:
F2a cost objective function for the integrated system environment;
CEfor integrated system environmental costs;
Pk(t) power of discharge source k at time t;
ωk,jis the emission coefficient of pollutant j from emission source k;
δE,jis a pollutant emission level;
ρEC-ppenalizing costs for pollutant emissions;
m is the total number of types of emissions;
w is the total number of emission sources;
the reliability-based objective function is:
Figure BDA0003113569280000141
in the formula:
F3a power supply reliability objective function for the comprehensive energy system;
the LESP is the energy supply loss rate of the comprehensive energy system;
t is the total operation time of the comprehensive energy system;
Edemand,tthe total energy demand in the comprehensive energy system at the time t;
Esupply,tand supplying energy demand in the comprehensive energy system at the moment t.
Step 102: determining constraint conditions of optimal scheduling of the integrated energy system, wherein the constraint conditions comprise: the system comprises a cold-heat-electricity load balance constraint, a power generation device constraint, a refrigeration and heating device constraint, a device climbing rate constraint, an electric energy transmission power constraint and an energy storage constraint between a comprehensive system and a power grid;
preferably, constraints for the optimal scheduling of the integrated energy system are determined, the constraints comprising: the energy storage system comprises a cold-heat-electricity load balance constraint, a power generation device constraint, a refrigeration and heating device constraint, a device climbing rate constraint, an electric energy transmission power constraint and an energy storage constraint between the comprehensive energy system and a power grid, wherein:
the cold, heat and electricity load balance constraint is as follows:
supposing that the heat load of a user is zero in summer, the waste heat of the micro-combustion engine is completely used for refrigerating by the absorption refrigerator; the user cold load is assumed to be zero in winter. The output electric power of each micro-source of the combined cooling heating and power system is the same as the demand of the electric load, namely:
Ppv(t)+Pwt(t)+Pmt(t)+Pgridbuy(t)+Pbat,d(t)＝Pec(t)+Pgridsell(t)+Pbat,c(t)+Pload(t)
in the formula: ppv(t)、Pwt(t)、Pmt(t)、Pgridbuy(t)、Pbat,d(t) photovoltaic power generation, wind power generation, micro gas turbine power generation, power purchase to a power grid and electric energy storage and discharge capacity are respectively carried out; pec(t)、Pgridsell(t)、Pbat,c(t)、Pload(t) the power consumption of the electric refrigerator, the power selling to the power grid, the charging amount of the electric energy storage and the demand of the electric load are respectively;
the output of the heating equipment in the comprehensive energy system is the same as the heat load demand, namely:
Qrec+Qgb+Qhs,d＝Qh,lond+Qc+Qhs,c
in the formula: qrec、Qgb、Qhs,dThe heat output by the waste heat boiler of the micro-gas turbine, the heat emitted by the gas boiler and the heat emitted by the heat energy storage device are respectively; qh,lond、Qc、Qhs,cThe heat load demand, the heat consumed by the absorption refrigerator and the heat storage quantity of the heat storage energy are respectively;
the output of the refrigeration equipment in the comprehensive energy system is the same as the cold load demand, namely:
Qac+Qec＝Qc,lond
in the formula: qac、QecRespectively the cold energy output by the absorption refrigerator and the cold energy generated by the electric refrigerator; qc,loadIs the cooling load demand;
the power generation equipment constraints are:
Figure BDA0003113569280000151
Figure BDA0003113569280000152
Figure BDA0003113569280000153
in the formula:
Figure BDA0003113569280000154
the maximum values of the output power of the wind power, the photovoltaic and the micro-combustion engine are respectively;
the constraint of the refrigerating and heating equipment is as follows:
Figure BDA0003113569280000155
Figure BDA0003113569280000156
Figure BDA0003113569280000157
in the formula: qgb、Qac、QecThe output powers of the gas boiler, the absorption refrigerator and the electric refrigerator are respectively;
Figure BDA0003113569280000158
the maximum output power values of the gas boiler, the absorption refrigerator and the electric refrigerator are respectively;
the equipment climbing rate constraint is as follows:
the output increasing rate of each device in the comprehensive energy system is smaller than the maximum upward climbing rate of the device; the rate of decrease of the output of the device is less than the maximum rate of downward hill climbing of the device;
-Ri,down＜Pi(t)-Pi(t-1)＜Ri,up
in the formula: pi(t)、Pi(t-1) power of the equipment at time t and (t-1) respectively; ri,down、Ri,upMaximum downward and maximum upward climbing rates;
the electric energy transmission power constraint between the comprehensive energy system and the power grid is as follows:
Figure BDA0003113569280000159
Figure BDA00031135692800001510
in the formula:
Figure BDA00031135692800001511
the upper limit values of electricity selling and electricity purchasing between the comprehensive energy system and the power grid are respectively set;
the energy storage constraint includes an electrical energy storage constraint and a thermal energy storage constraint, wherein the electrical energy storage constraint:
Figure BDA0003113569280000161
0≤Pbat,c≤Srated·γbat,c
0≤Pbat,d≤Srated·γbat,d
in the formula:
Figure BDA0003113569280000162
respectively storing minimum and maximum energy of electric energy; sratedRated capacity for electrical energy storage; gamma raybat,c、γbat,dThe maximum charging rate and the maximum discharging rate of the electric energy storage are respectively; pbat，d、Pbat，cDischarge power and charging power which are electric energy storage respectively; the charging and discharging of the electric energy storage belong to a coupling relation, and the electric energy storage and the discharging can not be carried out at the same time;
Pbat，c·Pbat，d＝0
the thermal energy storage constraint is:
Figure BDA0003113569280000163
Figure BDA0003113569280000164
Figure BDA0003113569280000165
in the formula:
Figure BDA0003113569280000166
the minimum and maximum heat storage quantity of the heat storage device are respectively;
Figure BDA0003113569280000167
the maximum heat storage power and the maximum heat release power of the heat storage device are respectively; qhs，c、Qhs，dThe heat storage power and the heat release power of the heat storage device are respectively; the heat storage device can not store heat and release heat at the same time, and the following constraint is added:
Qhs，c·Qhs，d＝0。
step 103: and acquiring the operation parameters of the comprehensive energy system, solving the objective function and the constraint condition based on the operation parameters of the comprehensive energy system, and generating an operation scheme for optimizing and scheduling the comprehensive energy system.
Preferably, solving the objective function and the constraint condition based on the operation parameters of the integrated energy system includes:
and solving the objective function and the constraint condition through one algorithm of a genetic algorithm, a particle swarm algorithm or a differential evolution algorithm.
The invention relates to an optimized dispatching model suitable for an electric, gas and thermal comprehensive energy system, which mainly comprises an objective function and constraint conditions and is a complex nonlinear multi-objective optimization model.
System optimization scheduling model objective function
The invention constructs the objective function of the system optimization scheduling model respectively from the three aspects of the running cost, the environmental cost and the reliability of the system. The concrete description is as follows:
(1) based on an operating cost objective function
The goal of integrated energy system economic dispatch is to minimize energy costs. The economic operating cost is mainly composed of the trading cost of natural gas, electric energy and heat energy. In addition to this, the operating costs of the system itself are included. The objective function is as follows:
F1＝minCop＝Cop-E+Cop-H+Cop-NG (1)
Figure BDA0003113569280000171
Figure BDA0003113569280000172
Figure BDA0003113569280000173
in the formula:
F1an integrated energy system operating cost objective function;
Copthe operating cost of the comprehensive energy system is reduced;
Cop-Ethe operating cost of the power supply module;
Cop-Hthe operating cost of the heat supply module;
Cop-NGthe acquisition cost of the natural gas consumed in the operation process of the comprehensive energy system;
Cithe operation cost of the distributed power generation equipment i under unit power;
Pi(t) is the output power of the distributed generation equipment i at the moment t;
n is the total number of distributed power generation equipment;
t is the total running time of the equipment;
Cgrid(t) is the real-time electricity price at time t;
Pgridsell(t) the power of the comprehensive energy system selling electricity to the power grid at the moment t;
Pgridbuy(t) the power purchased by the comprehensive energy system to the power grid at the moment t;
Celec-gridthe fixed operation maintenance cost of the electric energy supply module network in a unit time interval can be generally calculated by taking one hour as the unit time interval;
CHB-girdthe fixed operation maintenance cost of the heat supply module network in a unit time interval can be generally calculated by taking one hour as the unit time interval;
ah，jthe operating cost of the heat supply equipment j under unit power;
Ph,j(t) is the thermal output power of the heat supply device j at time t;
r is the total number of heat supply devices;
aNG，sthe cost of gas usage per unit power for a natural gas consuming device s;
PNG，s(t) the output power of the plant s consuming natural gas at time t;
g is the total number of equipment consuming natural gas.
(2) Objective function based on environmental cost
Micro-combustion engines and gas boilers using natural gas as fuel are important power and heat supply units of the system and are also important sources of pollutant emission of the system. The environmental cost of operating an integrated energy system mainly includes the following two aspects: environmental losses and non-environmental losses due to energy production pollutants; the pollution discharge fee charged by the relevant departments. The environmental cost minimization model is shown below.
Figure BDA0003113569280000181
In the formula:
F2a cost objective function for the integrated system environment;
CEfor integrated system environmental costs;
Pk(t) power of discharge source k at time t;
ωk,jis the emission coefficient of pollutant j from emission source k;
δE,jis a pollutant emission level;
ρEC-ppenalizing costs for pollutant emissions;
m is the total number of types of emissions;
w is the total number of emission sources.
(3) Reliability-based objective function
The system power shortage is a power supply reliability index commonly used by the power system. According to the correlation between the regional comprehensive energy system and the electric power system, the energy supply reliability index of the comprehensive energy system is constructed: loss of Energy Supply (LESP), which represents the ratio of the Energy Supply gap of the system to the total Energy demand over a certain period of time.
Figure BDA0003113569280000191
In the formula:
F3a power supply reliability objective function for the comprehensive energy system;
the LESP is the energy supply loss rate of the comprehensive energy system;
t is the total operation time of the comprehensive energy system;
Edemand,tthe total energy demand in the comprehensive energy system at the time t;
Esupply,tand supplying energy demand in the comprehensive energy system at the moment t.
Constraint condition of (II) system optimization scheduling model
The constraint conditions of the system optimization scheduling model mainly comprise cold, heat and power load balance constraint, output constraint of various types of equipment, power grid connection constraint and energy storage constraint. The concrete description is as follows:
(1) cold, heat, electric load balance constraint
Supposing that the heat load of users is zero in summer, the waste heat of the micro-combustion engine is completely used for refrigeration of the absorption refrigerator. The user cold load is assumed to be zero in winter. The output electric power of each micro-source of the combined cooling heating and power system is the same as the demand of the electric load, namely:
Ppv(t)+Pwt(t)+Pmt(t)+Pgridbuy(t)+Pbat,d(t)＝Pec(t)+Pgridsell(t)+Pbat,c(t)+Pload(t) (7)
in the formula: ppv(t)、Pwt(t)、Pmt(t)、Pgridbuy(t)、Pbat,d(t) photovoltaic power generation, wind power generation, micro gas turbine power generation, power purchase to a power grid and electric energy storage and discharge capacity are respectively carried out; pec(t)、Pgridsell(t)、Pbat,c(t)、PloadAnd (t) the power consumption of the electric refrigerator, the power selling amount to the power grid, the charging amount of the electric energy storage and the demand of the electric load are respectively.
The output of the heating equipment in the system is the same as the heat load demand, namely:
Qrec+Qgb+Qhs,d＝Qh,lond+Qc+Qhs,c (8)
in the formula: qrec、Qgb、Qhs,dThe heat output by the waste heat boiler of the micro-gas turbine, the heat emitted by the gas boiler and the heat emitted by the heat energy storage device are respectively; qh,lond、Qc、Qhs,cAbsorption refrigeration machines with heat load demand respectivelyHeat consumed by refrigeration, and heat stored by heat storage.
The output of the refrigeration equipment in the system is the same as the cold load demand, namely:
Qac+Qec＝Qc,lond (9)
in the formula: qac、QecRespectively the cold energy output by the absorption refrigerator and the cold energy generated by the electric refrigerator; qc,loadIs the cooling load demand.
(2) Power plant restraint
Figure BDA0003113569280000201
In the formula:
Figure BDA0003113569280000202
the maximum values of the output power of the wind power, the photovoltaic and the micro-combustion engine are respectively.
(3) Refrigeration and heating equipment restraint
Figure BDA0003113569280000203
In the formula: qgb、Qac、QecThe output powers of the gas boiler, the absorption refrigerator and the electric refrigerator are respectively;
Figure BDA0003113569280000204
the maximum output power values of the gas boiler, the absorption refrigerator and the electric refrigerator are respectively.
(4) Equipment ramp rate constraint
The output increasing rate of each device in the system is smaller than the maximum upward climbing rate of the device; the output decrease rate is less than the maximum ramp down rate.
-Ri,down＜Pi(t)-Pi(t-1)＜Ri,up (12)
In the formula: pi(t)、Pi(t-1) power of the equipment at time t and (t-1) respectively; ri,down、Ri,upThe maximum downward and upward climbing rates are respectively. Because of uncertainty of fan and photovoltaic output, the output is output according to the predicted power, so the climbing rate constraint of the micro-combustion engine and the gas boiler is mainly considered.
(5) Electric energy transmission power constraint between system and power grid
Figure BDA0003113569280000205
In the formula:
Figure BDA0003113569280000206
the upper limit values of electricity selling and electricity purchasing between the comprehensive energy system and the power grid are respectively.
(6) Restraint of stored energy
The electrical energy storage constraints are:
Figure BDA0003113569280000207
in the formula:
Figure BDA0003113569280000208
respectively storing minimum and maximum energy of electric energy; sratedRated capacity for electrical energy storage; gamma raybat,c、γbat,dThe maximum charging rate and the maximum discharging rate of the electric energy storage are respectively; pbat，d、Pbat，cRespectively the discharge power and the charging power of the electric energy storage. Furthermore, the charging and discharging of the electrical energy storage belong to a coupled relationship, and cannot both store and discharge electricity at the same time, so the following constraint is added:
Pbat，c·Pbat，d＝0 (15)
the thermal energy storage constraint is:
Figure BDA0003113569280000211
in the formula:
Figure BDA0003113569280000212
the minimum and maximum heat storage quantity of the heat storage device are respectively;
Figure BDA0003113569280000213
the maximum heat storage power and the maximum heat release power of the heat storage device are respectively; qhs.c、Qhs.dThe heat storage power and the heat release power of the heat storage device are respectively. Similar to electrical energy storage, the thermal storage device is at the same time, heat storage and heat release cannot be performed simultaneously, adding the following constraints:
Qhs，c·Qhs，d＝0 (17)
the system optimization scheduling model provided by the invention is a complex nonlinear multi-objective optimization model, and when the model is optimized and solved, an intelligent optimization algorithm can be adopted for solving, such as a genetic algorithm, a particle swarm algorithm, a differential evolution algorithm and the like.
The invention provides an optimized dispatching model suitable for an electric, gas and thermal comprehensive energy system, and compared with the prior art, the optimized dispatching model has the following beneficial effects: the model provided by the invention can solve the problem of coordinated optimization scheduling of various energy flows in a complex electrical, gas and thermal integrated energy system containing equipment elements such as photovoltaic, fans, micro-gas turbines, gas boilers, electrical energy storage, thermal energy storage, refrigeration/heating equipment, cooling and heating loads and the like, comprehensively considers the optimization of the running cost, the emission cost and the energy supply reliability of the system, and combines the consideration of the energy interaction with a main network, thereby realizing the comprehensive optimization of the whole system on the running economy, the environmental protection and the reliability.
Fig. 2 is a flowchart of a specific implementation method for performing optimal regulation and control on an integrated energy system by applying an optimal scheduling model applicable to an electrical, gas and thermal integrated energy system, which is provided by the present invention, and as shown in fig. 2, an embodiment of the method of the present invention totally comprises the following 3 steps:
step 1, inputting the operation related parameters of the comprehensive energy system into the comprehensive energy system optimization scheduling model provided by the invention. Wherein, mainly include: the method comprises the following steps of (1) defining parameters of a model, predicting data of photovoltaic, wind power and load in the day ahead, and integrating data information such as operation limit values of various devices in an energy system;
and 2, carrying out optimization solution on the system optimization scheduling model based on an intelligent optimization algorithm. Wherein, an artificial intelligent optimization algorithm such as a genetic algorithm, a particle swarm algorithm and the like can be selected to optimize and solve the model;
and 3, forming and outputting an optimized dispatching operation scheme of the comprehensive energy system.
Fig. 3 is a system configuration diagram for optimal scheduling of an integrated energy system according to a preferred embodiment of the present invention. As shown in fig. 3, the present invention provides a system for optimized scheduling of an integrated energy system, the system comprising:
the establishing unit 301 is configured to establish an objective function of optimized scheduling of the integrated energy system, where the objective function includes: an objective function based on operating cost, an objective function based on environmental cost, and an objective function based on reliability; preferably, the integrated energy system comprises: electrical energy, gas energy, and thermal energy.
Preferably, the establishing unit 301 is configured to establish an objective function of the optimized scheduling of the integrated energy system, where the objective function includes: an objective function based on operating cost, an objective function based on environmental cost, and an objective function based on reliability, wherein:
the objective function based on the running cost is:
F1＝minCop＝Cop-E+Cop-H+Cop-NG
Figure BDA0003113569280000221
Figure BDA0003113569280000222
Figure BDA0003113569280000223
in the formula:
F1an integrated energy system operating cost objective function;
Copthe operating cost of the comprehensive energy system is reduced;
Cop-Ethe operating cost of the power supply module;
Cop-Hthe operating cost of the heat supply module;
Cop-NGthe acquisition cost of the natural gas consumed in the operation process of the comprehensive energy system;
Cithe operation cost of the distributed power generation equipment i under unit power;
Pi(t) is the output power of the distributed generation equipment i at the moment t;
n is the total number of distributed power generation equipment;
t is the total running time of the equipment;
Cgrid(t) is the real-time electricity price at time t;
Pgridsell(t) the power of the comprehensive energy system selling electricity to the power grid at the moment t;
Pgridbuy(t) the power purchased by the comprehensive energy system to the power grid at the moment t;
Celec-gridthe fixed operation maintenance cost of the electric energy supply module network in a unit time interval can be generally calculated by taking one hour as the unit time interval;
CHB-girdthe fixed operation maintenance cost of the heat supply module network in a unit time interval can be generally calculated by taking one hour as the unit time interval;
ah，jthe operating cost of the heat supply equipment j under unit power;
Ph,j(t) is the thermal output power of the heat supply device j at time t;
r is the total number of heat supply devices;
aNG，scost of gas consumption per unit power for natural gas consuming units s；
PMG，s(t) the output power of the plant s consuming natural gas at time t;
g is the total number of devices consuming natural gas;
the objective function based on the environmental cost is:
Figure BDA0003113569280000231
in the formula:
F2a cost objective function for the integrated system environment;
CEfor integrated system environmental costs;
Pk(t) power of discharge source k at time t;
ωk,jis the emission coefficient of pollutant j from emission source k;
δE,jis a pollutant emission level;
ρEC-ppenalizing costs for pollutant emissions;
m is the total number of types of emissions;
w is the total number of emission sources;
the reliability-based objective function is:
Figure BDA0003113569280000232
in the formula:
F3a power supply reliability objective function for the comprehensive energy system;
the LESP is the energy supply loss rate of the comprehensive energy system;
t is the total operation time of the comprehensive energy system;
Edemand,tthe total energy demand in the comprehensive energy system at the time t;
Esupply,tand supplying energy demand in the comprehensive energy system at the moment t.
A determining unit 302, configured to determine constraints for optimal scheduling of the integrated energy system, where the constraints include: the system comprises a cold-heat-electricity load balance constraint, a power generation device constraint, a refrigeration and heating device constraint, a device climbing rate constraint, an electric energy transmission power constraint and an energy storage constraint between a comprehensive system and a power grid;
and the result unit 303 is configured to obtain the operation parameters of the integrated energy system, solve the objective function and the constraint condition based on the operation parameters of the integrated energy system, and generate an operation scheme for optimizing scheduling of the integrated energy system.
Preferably, the determining unit 302 is configured to determine constraints for optimizing scheduling of the integrated energy system, where the constraints include: the energy storage system comprises a cold-heat-electricity load balance constraint, a power generation device constraint, a refrigeration and heating device constraint, a device climbing rate constraint, an electric energy transmission power constraint and an energy storage constraint between the comprehensive energy system and a power grid, wherein:
the cold, heat and electricity load balance constraint is as follows:
supposing that the heat load of a user is zero in summer, the waste heat of the micro-combustion engine is completely used for refrigerating by the absorption refrigerator; the user cold load is assumed to be zero in winter. The output electric power of each micro-source of the combined cooling heating and power system is the same as the demand of the electric load, namely:
Ppv(t)+Pwt(t)+Pmt(t)+Pgridbuy(t)+Pbat,d(t)＝Pec(t)+Pgridsell(t)+Pbat,c(t)+Pload(t)
in the formula: ppv(t)、Pwt(t)、Pmt(t)、Pgridbuy(t)、Pbat,d(t) photovoltaic power generation, wind power generation, micro gas turbine power generation, power purchase to a power grid and electric energy storage and discharge capacity are respectively carried out; pec(t)、Pgridsell(t)、Pbat,c(t)、Pload(t) the power consumption of the electric refrigerator, the power selling to the power grid, the charging amount of the electric energy storage and the demand of the electric load are respectively;
the output of the heating equipment in the comprehensive energy system is the same as the heat load demand, namely:
Qrec+Qgb+Qhs,d＝Qh,lond+Qc+Qhs,c
in the formula: qrec、Qgb、Qhs,dThe heat output by the waste heat boiler of the micro-gas turbine, the heat emitted by the gas boiler and the heat emitted by the heat energy storage device are respectively; qh,lond、Qc、Qhs,cThe heat load demand, the heat consumed by the absorption refrigerator and the heat storage quantity of the heat storage energy are respectively;
the output of the refrigeration equipment in the comprehensive energy system is the same as the cold load demand, namely:
Qac+Qec＝Qc,lond
in the formula: qac、QecRespectively the cold energy output by the absorption refrigerator and the cold energy generated by the electric refrigerator; qc,loadIs the cooling load demand;
the power generation equipment constraints are:
Figure BDA0003113569280000251
Figure BDA0003113569280000252
Figure BDA0003113569280000253
in the formula:
Figure BDA0003113569280000254
the maximum values of the output power of the wind power, the photovoltaic and the micro-combustion engine are respectively;
the constraint of the refrigerating and heating equipment is as follows:
Figure BDA0003113569280000255
Figure BDA0003113569280000256
Figure BDA0003113569280000257
in the formula: qgb、Qac、QecThe output powers of the gas boiler, the absorption refrigerator and the electric refrigerator are respectively;
Figure BDA0003113569280000258
the maximum output power values of the gas boiler, the absorption refrigerator and the electric refrigerator are respectively;
the equipment climbing rate constraint is as follows:
the output increasing rate of each device in the comprehensive energy system is smaller than the maximum upward climbing rate of the device; the rate of decrease of the output of the device is less than the maximum rate of downward hill climbing of the device;
-Ri,down＜Pi(t)-Pi(t-1)＜Ri,up
in the formula: pi(t)、Pi(t-1) power of the equipment at time t and (t-1) respectively; ri,down、Ri,upMaximum downward and maximum upward climbing rates;
the electric energy transmission power constraint between the comprehensive energy system and the power grid is as follows:
Figure BDA0003113569280000259
Figure BDA00031135692800002510
in the formula:
Figure BDA00031135692800002511
the upper limit values of electricity selling and electricity purchasing between the comprehensive energy system and the power grid are respectively set;
the energy storage constraint includes an electrical energy storage constraint and a thermal energy storage constraint, wherein the electrical energy storage constraint:
Figure BDA00031135692800002512
0≤Pbat,c≤Srated·γbat,c
0≤Pbat,d≤Srated·γbat,d
in the formula:
Figure BDA0003113569280000261
respectively storing minimum and maximum energy of electric energy; sratedRated capacity for electrical energy storage; gamma raybat,c、γbat,dThe maximum charging rate and the maximum discharging rate of the electric energy storage are respectively; pbat，d、Pbat，cDischarge power and charging power which are electric energy storage respectively; the charging and discharging of the electric energy storage belong to a coupling relation, and the electric energy storage and the discharging can not be carried out at the same time;
Pbat，c·Pbat，d＝0
the thermal energy storage constraint is:
Figure BDA0003113569280000262
Figure BDA0003113569280000263
Figure BDA0003113569280000264
in the formula:
Figure BDA0003113569280000265
the minimum and maximum heat storage quantity of the heat storage device are respectively;
Figure BDA0003113569280000266
the maximum heat storage power and the maximum heat release power of the heat storage device are respectively; qhs.c、Qhs.dRespectively the heat storage power of the heat storage deviceHeat release power; the heat storage device can not store heat and release heat at the same time, and the following constraint is added:
Qhs，c·Qhs，d＝0。
preferably, the result unit 303 is configured to solve the objective function and the constraint condition based on the operation parameters of the integrated energy system, and includes:
and solving the objective function and the constraint condition through one algorithm of a genetic algorithm, a particle swarm algorithm or a differential evolution algorithm.
The system 300 for optimal scheduling of an integrated energy system according to a preferred embodiment of the present invention corresponds to the method 100 for optimal scheduling of an integrated energy system according to another preferred embodiment of the present invention, and will not be described herein again.
The invention has been described with reference to a few embodiments. However, other embodiments of the invention than the one disclosed above are equally possible within the scope of the invention, as would be apparent to a person skilled in the art from the appended patent claims. As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting the same, and although the present invention is described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.
Generally, all terms used in the claims are to be interpreted according to their ordinary meaning in the technical field, unless explicitly defined otherwise herein. All references to "a/an/the [ device, component, etc ]" are to be interpreted openly as referring to at least one instance of said device, component, etc., unless explicitly stated otherwise. The steps of any method disclosed herein do not have to be performed in the exact order disclosed, unless explicitly stated.

Claims (10)

1. A method for optimized scheduling of an integrated energy system, the method comprising:
determining an objective function of optimal scheduling of the comprehensive energy system and a constraint condition;
and acquiring the operation parameters of the comprehensive energy system, solving the objective function and the constraint condition based on the operation parameters, and generating an operation scheme for optimizing and scheduling the comprehensive energy system.
2. The method of claim 1, the integrated energy system comprising at least one of: electrical energy, gas energy, and thermal energy.
3. The method of claim 1, the solving the objective function and constraints based on the operating parameters, comprising:
and solving the objective function and the constraint condition through one algorithm of a genetic algorithm, a particle swarm algorithm or a differential evolution algorithm.
4. The method of claim 1, the establishing an objective function for optimized scheduling of an integrated energy system, the objective function comprising: an objective function based on operating cost, an objective function based on environmental cost, and an objective function based on reliability, wherein:
the objective function based on the running cost is:
F1＝minCop＝Cop-E+Cop-H+Cop-NG
Figure FDA0003113569270000011
Figure FDA0003113569270000012
Figure FDA0003113569270000013
in the formula:
F1an integrated energy system operating cost objective function;
Copthe operating cost of the comprehensive energy system is reduced;
Cop-Ethe operating cost of the power supply module;
Cop-Hthe operating cost of the heat supply module;
Cop-NGthe acquisition cost of the natural gas consumed in the operation process of the comprehensive energy system;
Cithe operation cost of the distributed power generation equipment i under unit power;
Pi(t) is the output power of the distributed generation equipment i at the moment t;
n is the total number of distributed power generation equipment;
t is the total running time of the equipment;
Cgrid(t) is the real-time electricity price at time t;
Pgridsell(t) the power of the comprehensive energy system selling electricity to the power grid at the moment t;
Pgridbuy(t) the power purchased by the comprehensive energy system to the power grid at the moment t;
Celec-gridthe fixed operation maintenance cost of the electric energy supply module network in a unit time interval can be generally calculated by taking one hour as the unit time interval;
CHB-girdthe fixed operation maintenance cost of the heat supply module network in a unit time interval;
αh，jthe operating cost of the heat supply equipment j under unit power;
Ph,j(t) is the thermal output power of the heat supply device j at time t;
r is the total number of heat supply devices;
αNG，sthe cost of gas usage per unit power for a natural gas consuming device s;
PNG，s(t) the output power of the plant s consuming natural gas at time t;
g is the total number of devices consuming natural gas;
the objective function based on the environmental cost is:
Figure FDA0003113569270000021
in the formula:
F2a cost objective function for the integrated system environment;
CEfor integrated system environmental costs;
Pk(t) power of discharge source k at time t;
ωk,jis the emission coefficient of pollutant j from emission source k;
δE,jis a pollutant emission level;
ρEC-ppenalizing costs for pollutant emissions;
m is the total number of types of emissions;
w is the total number of emission sources;
the reliability-based objective function is:
Figure FDA0003113569270000031
in the formula:
F3a power supply reliability objective function for the comprehensive energy system;
the LESP is the energy supply loss rate of the comprehensive energy system;
t is the total operation time of the comprehensive energy system;
Edemand,tthe total energy demand in the comprehensive energy system at the time t;
Esupply,tand supplying energy demand in the comprehensive energy system at the moment t.
5. The method of claim 1, the determining constraints for optimal scheduling of the integrated energy system, the constraints comprising: the energy storage system comprises a cold-heat-electricity load balance constraint, a power generation device constraint, a refrigeration and heating device constraint, a device climbing rate constraint, an electric energy transmission power constraint and an energy storage constraint between the comprehensive energy system and a power grid, wherein:
the cooling, heating and power load balance constraint is calculated according to the following formula:
Ppv(t)+Pwt(t)+Pmt(t)+Pgridbuy(t)+Pbat,d(t)＝Pec(t)+Pgridsell(t)+Pbat,c(t)+Pload(t)
in the formula: ppv(t)、Pwt(t)、Pmt(t)、Pgridbuy(t)、Pbat,d(t) photovoltaic power generation, wind power generation, micro gas turbine power generation, power purchase to a power grid and electric energy storage and discharge capacity are respectively carried out; pec(t)、Pgridsell(t)、Pbat,c(t)、Pload(t) the power consumption of the electric refrigerator, the power selling to the power grid, the charging amount of the electric energy storage and the demand of the electric load are respectively;
the output of the heating equipment in the comprehensive energy system is the same as the heat load demand, and the calculation is carried out according to the following formula:
Qrec+Qgb+Qhs,d＝Qh,lond+Qc+Qhs,c
in the formula: qrec、Qgb、Qhs,dThe heat output by the waste heat boiler of the micro-gas turbine, the heat emitted by the gas boiler and the heat emitted by the heat energy storage device are respectively; qh,lond、Qc、Qhs,cThe heat load demand, the heat consumed by the absorption refrigerator and the heat storage quantity of the heat storage energy are respectively;
the output of the refrigeration equipment in the comprehensive energy system is the same as the cold load demand, and the calculation is carried out according to the following formula:
Qac+Qec＝Qc,lond
in the formula: qac、QecRespectively the cold energy output by the absorption refrigerator and the cold energy generated by the electric refrigerator; qc,loadIs the cooling load demand;
the power plant constraints are calculated according to the following formula:
Figure FDA0003113569270000041
Figure FDA0003113569270000042
Figure FDA0003113569270000043
in the formula:
Figure FDA0003113569270000044
the maximum values of the output power of the wind power, the photovoltaic and the micro-combustion engine are respectively;
the constraints of the refrigerating and heating equipment are calculated according to the following formula:
Figure FDA0003113569270000045
Figure FDA0003113569270000046
Figure FDA0003113569270000047
in the formula: qgb、Qac、QecThe output powers of the gas boiler, the absorption refrigerator and the electric refrigerator are respectively;
Figure FDA0003113569270000048
the maximum output power values of the gas boiler, the absorption refrigerator and the electric refrigerator are respectively;
the equipment ramp rate constraint is calculated according to the following formula:
-Ri,down＜Pi(t)-Pi(t-1)＜Ri,up
in the formula: pi(t)、Pi(t-1) power of the equipment at time t and (t-1) respectively; ri,down、Ri,upMaximum downward and maximum upward climbing rates;
the electric energy transmission power constraint between the comprehensive energy system and the power grid is calculated according to the following formula:
Figure FDA0003113569270000049
Figure FDA00031135692700000410
in the formula:
Figure FDA00031135692700000411
the upper limit values of electricity selling and electricity purchasing between the comprehensive energy system and the power grid are respectively set;
the energy storage constraints include electrical energy storage constraints and thermal energy storage constraints, wherein the electrical energy storage constraints are calculated according to the following formula:
Figure FDA00031135692700000412
0≤Pbat,c≤Srated·γbat,c
0≤Pbat,d≤Srated·γbat,d
in the formula:
Figure FDA00031135692700000413
respectively storing minimum and maximum energy of electric energy; sratedRated capacity for electrical energy storage; gamma raybat,c、γbat,dRespectively for storing energyMaximum charge rate, maximum discharge rate; pbat，d、Pbat，cDischarge power and charging power which are electric energy storage respectively; the charging and discharging of the electric energy storage belong to a coupling relation, and the electric energy storage and the discharging can not be carried out at the same time;
Pbat，c·Pbat，d＝0
the thermal energy storage constraint is calculated according to the following formula:
Figure FDA0003113569270000051
Figure FDA0003113569270000052
Figure FDA0003113569270000053
in the formula:
Figure FDA0003113569270000054
the minimum and maximum heat storage quantity of the heat storage device are respectively;
Figure FDA0003113569270000055
the maximum heat storage power and the maximum heat release power of the heat storage device are respectively; qhs，c、Qhs，dThe heat storage power and the heat release power of the heat storage device are respectively; the heat storage device can not store heat and release heat at the same time, and the following constraint is added:
Qhs，c·Qhs，d＝0。
6. a system for optimized scheduling of an integrated energy system, the system comprising:
the system comprises an establishing unit, a scheduling unit and a scheduling unit, wherein the establishing unit is used for determining an objective function of optimal scheduling of the comprehensive energy system and a constraint condition;
and the result unit is used for acquiring the operation parameters of the comprehensive energy system, solving the objective function and the constraint condition based on the operation parameters and generating an operation scheme for optimizing and scheduling the comprehensive energy system.
7. The system of claim 6, the integrated energy system comprising at least one of: electrical energy, gas energy, and thermal energy.
8. The system of claim 6, the result unit, in particular, configured to solve the objective function and the constraint condition by one of a genetic algorithm, a particle swarm algorithm, or a differential evolution algorithm.
9. The system of claim 8, the establishing unit configured to establish an objective function for optimal scheduling of an integrated energy system, the objective function comprising: an objective function based on operating cost, an objective function based on environmental cost, and an objective function based on reliability, wherein:
the objective function based on the running cost is:
F1＝minCop＝Cop-E+Cop-H+Cop-NG
Figure FDA0003113569270000061
Figure FDA0003113569270000062
Figure FDA0003113569270000063
in the formula:
F1an integrated energy system operating cost objective function;
Copthe operating cost of the comprehensive energy system is reduced;
Cop-Ethe operating cost of the power supply module;
Cop-Hthe operating cost of the heat supply module;
Cop-NGthe acquisition cost of the natural gas consumed in the operation process of the comprehensive energy system;
Cithe operation cost of the distributed power generation equipment i under unit power;
Pi(t) is the output power of the distributed generation equipment i at the moment t;
n is the total number of distributed power generation equipment;
t is the total running time of the equipment;
Cgrid(t) is the real-time electricity price at time t;
Pgridsell(t) the power of the comprehensive energy system selling electricity to the power grid at the moment t;
Pgridbuy(t) the power purchased by the comprehensive energy system to the power grid at the moment t;
Celec-gridthe fixed operation maintenance cost of the electric energy supply module network in a unit time interval;
CHB-girdthe fixed operation maintenance cost of the heat supply module network in a unit time interval;
ah，jthe operating cost of the heat supply equipment j under unit power;
Ph,j(t) is the thermal output power of the heat supply device j at time t;
r is the total number of heat supply devices;
aNG，sthe cost of gas usage per unit power for a natural gas consuming device s;
PNG，s(t) the output power of the plant s consuming natural gas at time t;
g is the total number of devices consuming natural gas;
the objective function based on the environmental cost is:
Figure FDA0003113569270000071
in the formula:
F2an environmental cost objective function of the integrated energy system;
CEenvironmental cost for the integrated energy system;
Pk(t) power of discharge source k at time t;
ωk,jis the emission coefficient of pollutant j from emission source k;
δE,jis a pollutant emission level;
ρEC-ppenalizing costs for pollutant emissions;
m is the total number of types of emissions;
w is the total number of emission sources;
the reliability-based objective function is:
Figure FDA0003113569270000072
in the formula:
F3a power supply reliability objective function for the comprehensive energy system;
the LESP is the energy supply loss rate of the comprehensive energy system;
t is the total operation time of the comprehensive energy system;
Edemand,tthe total energy demand in the comprehensive energy system at the time t;
Esupply,tand supplying energy demand in the comprehensive energy system at the moment t.
10. The system of claim 6, the validation unit to determine constraints for optimal scheduling of the integrated energy system, the constraints comprising: the energy storage system comprises a cold-heat-electricity load balance constraint, a power generation device constraint, a refrigeration and heating device constraint, a device climbing rate constraint, an electric energy transmission power constraint and an energy storage constraint between the comprehensive energy system and a power grid, wherein:
the cooling, heating and power load balance constraint is calculated according to the following formula:
Ppv(t)+Pwt(t)+Pmt(t)+Pgridbuy(t)+Pbat,d(t)＝Pec(t)+Pgridsell(t)+Pbat,c(t)+Pload(t)
in the formula: ppv(t)、Pwt(t)、Pmt(t)、Pgridbuy(t)、Pbat,d(t) photovoltaic power generation, wind power generation, micro gas turbine power generation, power purchase to a power grid and electric energy storage and discharge capacity are respectively carried out; pec(t)、Pgridsell(t)、Pbat,c(t)、Pload(t) the power consumption of the electric refrigerator, the power selling to the power grid, the charging amount of the electric energy storage and the demand of the electric load are respectively;
the output of the heating equipment in the comprehensive energy system is the same as the heat load demand, and the calculation is carried out according to the following formula:
Qrec+Qgb+Qhs,d＝Qh,lond+Qc+Qhs,c
in the formula: qrec、Qgb、Qhs,dThe heat output by the waste heat boiler of the micro-gas turbine, the heat emitted by the gas boiler and the heat emitted by the heat energy storage device are respectively; qh,lond、Qc、Qhs,cThe heat load demand, the heat consumed by the absorption refrigerator and the heat storage quantity of the heat storage energy are respectively;
the output of the refrigeration equipment in the comprehensive energy system is the same as the cold load demand, and the calculation is carried out according to the following formula:
Qac+Qec＝Qc,lond
in the formula: qac、QecRespectively the cold energy output by the absorption refrigerator and the cold energy generated by the electric refrigerator; qc,loadIs the cooling load demand;
the power plant constraints are calculated according to the following formula:
Figure FDA0003113569270000081
Figure FDA0003113569270000082
Figure FDA0003113569270000083
in the formula:
Figure FDA0003113569270000084
the maximum values of the output power of the wind power, the photovoltaic and the micro-combustion engine are respectively;
the constraints of the refrigerating and heating equipment are calculated according to the following formula:
Figure FDA0003113569270000085
Figure FDA0003113569270000086
Figure FDA0003113569270000087
in the formula: qgb、Qac、QecThe output powers of the gas boiler, the absorption refrigerator and the electric refrigerator are respectively;
Figure FDA0003113569270000088
the maximum output power values of the gas boiler, the absorption refrigerator and the electric refrigerator are respectively;
the equipment ramp rate constraint is calculated according to the following formula:
-Ri,down＜Pi(t)-Pi(t-1)＜Ri,up
in the formula: pi(t)、Pi(t-1) power of the equipment at time t and (t-1) respectively; ri,down、Ri,upRespectively, maximum downwardThe maximum upward climbing speed;
the electric energy transmission power constraint between the comprehensive energy system and the power grid is calculated according to the following formula:
Figure FDA0003113569270000091
Figure FDA0003113569270000092
in the formula:
Figure FDA0003113569270000093
the upper limit values of electricity selling and electricity purchasing between the comprehensive energy system and the power grid are respectively set;
the energy storage constraints include electrical energy storage constraints and thermal energy storage constraints, wherein the electrical energy storage constraints are calculated according to the following formula:
Figure FDA0003113569270000094
0≤Pbat,c≤Srated·γbat,c
0≤Pbat,d≤Srated·γbat,d
in the formula:
Figure FDA0003113569270000095
respectively storing minimum and maximum energy of electric energy; sratedRated capacity for electrical energy storage; gamma raybat,c、γbat,dThe maximum charging rate and the maximum discharging rate of the electric energy storage are respectively; pbat，d、Pbat，cDischarge power and charging power which are electric energy storage respectively; the charging and discharging of the electric energy storage belong to a coupling relation, and the electric energy storage and the discharging can not be carried out at the same time;
Pbat，c·Pbat，d＝0
the thermal energy storage constraint is calculated according to the following formula:
Figure FDA0003113569270000096
Figure FDA0003113569270000097
Figure FDA0003113569270000098
in the formula:
Figure FDA0003113569270000099
the minimum and maximum heat storage quantity of the heat storage device are respectively;
Figure FDA00031135692700000910
the maximum heat storage power and the maximum heat release power of the heat storage device are respectively; qhs，c、Qhs，dThe heat storage power and the heat release power of the heat storage device are respectively; the heat storage device can not store heat and release heat at the same time, and the following constraint is added:
Qhs，c·Qhs，d＝0。
CN202110655638.9A 2021-06-11 2021-06-11 Method and system for optimizing scheduling of comprehensive energy system Pending CN113435095A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110655638.9A CN113435095A (en) 2021-06-11 2021-06-11 Method and system for optimizing scheduling of comprehensive energy system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110655638.9A CN113435095A (en) 2021-06-11 2021-06-11 Method and system for optimizing scheduling of comprehensive energy system

Publications (1)

Publication Number Publication Date
CN113435095A true CN113435095A (en) 2021-09-24

Family

ID=77755824

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110655638.9A Pending CN113435095A (en) 2021-06-11 2021-06-11 Method and system for optimizing scheduling of comprehensive energy system

Country Status (1)

Country Link
CN (1) CN113435095A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113872192A (en) * 2021-09-26 2021-12-31 国网电力科学研究院武汉能效测评有限公司 Hospital power grid load optimization control system and control method
CN114545878A (en) * 2022-02-22 2022-05-27 山东大学 Optimized scheduling method and system for comprehensive energy system
CN115630434A (en) * 2022-12-25 2023-01-20 国网上海能源互联网研究院有限公司 Building carbon emission prediction method and device based on multi-agent simulation
CN117611234A (en) * 2023-11-17 2024-02-27 山东正晨科技股份有限公司 Prediction optimization method, system, equipment and medium for comprehensive energy system of expressway

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113872192A (en) * 2021-09-26 2021-12-31 国网电力科学研究院武汉能效测评有限公司 Hospital power grid load optimization control system and control method
CN113872192B (en) * 2021-09-26 2024-03-12 国网电力科学研究院武汉能效测评有限公司 Hospital power grid load optimization control system and control method
CN114545878A (en) * 2022-02-22 2022-05-27 山东大学 Optimized scheduling method and system for comprehensive energy system
CN114545878B (en) * 2022-02-22 2023-09-01 山东大学 Comprehensive energy system optimal scheduling method and system
CN115630434A (en) * 2022-12-25 2023-01-20 国网上海能源互联网研究院有限公司 Building carbon emission prediction method and device based on multi-agent simulation
CN117611234A (en) * 2023-11-17 2024-02-27 山东正晨科技股份有限公司 Prediction optimization method, system, equipment and medium for comprehensive energy system of expressway

Similar Documents

Publication Publication Date Title
CN107832979B (en) Factory comprehensive energy system economic optimization scheduling method considering energy cascade utilization
CN113435095A (en) Method and system for optimizing scheduling of comprehensive energy system
CN109146182A (en) The economic load dispatching method of meter and the distributed triple-generation system of a variety of energy storage
CN113112087A (en) Comprehensive energy system operation cost optimization method considering electric heating load demand response
CN109523065B (en) Micro energy network optimization scheduling method based on improved quantum particle swarm algorithm
CN112464477A (en) Multi-energy coupling comprehensive energy operation simulation method considering demand response
CN107807523A (en) Consider the Regional Energy internet multi-source coordination optimization operation reserve of tou power price
CN106022503A (en) Micro-grid capacity programming method meeting coupling type electric cold and heat demand
CN106786509B (en) Based on the thermo-electrically combined scheduling method of more scenario simulations under large-scale wind power integration
CN110288152A (en) Consider electricity/thermal flexibility load regional complex energy resource system energy storage configuration method
CN110138004A (en) One kind is provided multiple forms of energy to complement each other system optimized operation method
CN111737884B (en) Multi-target random planning method for micro-energy network containing multiple clean energy sources
Ren et al. A novel planning method for design and dispatch of hybrid energy systems
CN110361969B (en) Optimized operation method of cooling, heating and power comprehensive energy system
CN111008739A (en) Optimal regulation and control and income distribution method and system for cogeneration virtual power plant
CN110391655B (en) Multi-energy-coupling micro-energy-network economic optimization scheduling method and device
CN109543889A (en) A kind of regional complex energy resource system cooperates with optimizing operation method a few days ago
CN108596442A (en) The integrated energy system economic load dispatching method of meter and Conditional Lyapunov ExponentP
CN113256045A (en) Park comprehensive energy system day-ahead economic dispatching method considering wind and light uncertainty
CN113128868B (en) Regional comprehensive energy system scheduling optimization method and device
Guo et al. Operation optimization of integrated energy system from the perspective of sustainable development
CN109255487A (en) A kind of integrated energy system optimization method based on normalized matrix model
CN117611382A (en) Rural comprehensive energy system optimization method considering multilayer cooperation and demand response
Yang et al. Integrated energy scheduling using wind-storage cascade control strategy in dynamic peak regulation right
CN113128799A (en) Energy management and control method and device, electronic equipment and computer storage medium

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