CN112989594B - Comprehensive energy system operation optimization method considering hydrogen energy - Google Patents

Comprehensive energy system operation optimization method considering hydrogen energy Download PDF

Info

Publication number
CN112989594B
CN112989594B CN202110255688.8A CN202110255688A CN112989594B CN 112989594 B CN112989594 B CN 112989594B CN 202110255688 A CN202110255688 A CN 202110255688A CN 112989594 B CN112989594 B CN 112989594B
Authority
CN
China
Prior art keywords
hydrogen
energy
power
gas
cost
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110255688.8A
Other languages
Chinese (zh)
Other versions
CN112989594A (en
Inventor
李晨
郭祚刚
雷金勇
徐敏
袁智勇
叶琳浩
王�琦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China South Power Grid International Co ltd
China Southern Power Grid Co Ltd
Original Assignee
China South Power Grid International Co ltd
China Southern Power Grid 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 China South Power Grid International Co ltd, China Southern Power Grid Co Ltd filed Critical China South Power Grid International Co ltd
Priority to CN202110255688.8A priority Critical patent/CN112989594B/en
Publication of CN112989594A publication Critical patent/CN112989594A/en
Application granted granted Critical
Publication of CN112989594B publication Critical patent/CN112989594B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0631Resource planning, allocation, distributing or scheduling for enterprises or organisations
    • G06Q10/06313Resource planning in a project environment
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

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

Abstract

The invention provides a comprehensive energy system operation optimization method considering hydrogen energy, which comprises the following steps: aiming at a comprehensive energy system considering hydrogen energy, the comprehensive energy system comprises an electric power system, a natural gas system and a hydrogen energy storage system, wherein equipment operation constraint, electric power network constraint and coupling constraint among subsystems are taken as operation constraint conditions, and minimized operation cost is taken as an objective function, and an operation optimization model of the comprehensive energy system is established; inputting operation parameters of the comprehensive energy system into the operation optimization model, and solving the operation optimization model to obtain an optimized operation scheme of the comprehensive energy system; the operation parameters of the comprehensive energy system comprise system structure parameters, load parameters, equipment operation parameters, energy transmission loss parameters, environment cost parameters and benefit parameters of selling hydrogen. The energy coupling constraint among the power system, the natural gas system and the hydrogen energy storage system is fully considered, the optimization model is more comprehensive and reasonable, and the obtained optimized operation scheme is more practical.

Description

Comprehensive energy system operation optimization method considering hydrogen energy
Technical Field
The invention relates to the technical field of comprehensive energy systems, in particular to an operation optimization method of a comprehensive energy system considering hydrogen energy.
Background
The traditional energy system has the defects of single energy situation and low energy utilization rate. The comprehensive energy system can comprehensively utilize various energy sources such as electricity, heat, cold, gas and the like, and becomes an important point of energy field development in recent years due to higher energy utilization efficiency. However, the environmental protection degree of the comprehensive energy system still has room for improvement, and the comprehensive energy system taking the participation of hydrogen energy into consideration becomes one of key technologies for further improving the environmental protection efficiency.
Hydrogen energy is one of the important directions for the development of new energy in recent years and even in the future. The hydrogen energy has the advantages of high energy utilization rate, low carbon and environmental protection, is highly focused in the energy field, and is one of technical means for realizing the aim of carbon neutralization. Hydrogen energy has the disadvantage of being difficult to store on a large scale and transport over long distances. However, natural gas systems can be transported and stored on a large scale, but their energy conversion efficiency is only 45% -65%. Therefore, the comprehensive energy system effectively coordinates various energy sources and the characteristic of multi-energy complementation, can effectively improve the utilization efficiency of various energy sources and realizes the economic and efficient operation of the system.
Most of the existing comprehensive energy system operation optimization methods only consider participation operation scenes of an electric power system, a natural gas system and the like, do not consider participation operation scenes of a hydrogen energy system, do not fully consider operation cost of the comprehensive energy system, do not consider environmental benefits brought by preparing methane through catalytic reaction of hydrogen energy and carbon dioxide so as to reduce greenhouse gases, do not consider economic benefits brought by supplying hydrogen to an oxyhydrogen fuel cell and a hydrogen energy automobile at a higher price, do not fully consider energy conversion relations among various energy subsystems, and can lead to that the obtained comprehensive energy system operation scheme cannot achieve expected effects in the practical application process, and even cannot be viable.
Disclosure of Invention
The invention aims to provide an integrated energy system operation optimization method considering hydrogen energy, so as to solve the technical problems that the hydrogen energy system is not considered to participate in operation and the operation scheme of the integrated energy system is inconsistent with actual application in the related technology.
The aim of the invention can be achieved by the following technical scheme:
an integrated energy system operation optimization method considering hydrogen energy, comprising:
aiming at a comprehensive energy system considering hydrogen energy, the comprehensive energy system comprises an electric power system, a natural gas system and a hydrogen energy storage system, wherein equipment operation constraint, electric power network constraint and coupling constraint among subsystems are taken as operation constraint conditions, and minimized operation cost is taken as an objective function, and an operation optimization model of the comprehensive energy system is established;
inputting operation parameters of the comprehensive energy system into the operation optimization model, and solving the operation optimization model to obtain an optimized operation scheme of the comprehensive energy system; the operation parameters of the comprehensive energy system comprise system structure parameters, load parameters, equipment operation parameters, energy transmission loss parameters, environment cost parameters and benefit parameters of selling hydrogen.
Further, the objective function is:
z=min(f 1 +f 2 +f 3 +f 4 -f 5 );
wherein z is the running cost of the integrated energy system; f (f) 1 、f 2 、f 3 、f 4 、f 5 Separate tableShowing equipment operating cost, energy conversion cost, environmental cost, transmission loss cost and benefits brought by selling hydrogen;
Figure BDA0002968280140000021
Figure BDA0002968280140000022
the operation cost of the thermal power generating unit, the gas turbine and the fuel cell is respectively;
Figure BDA0002968280140000023
Figure BDA0002968280140000024
The operation cost of the electricity storage equipment, the hydrogen storage equipment and the gas storage equipment is as follows;
Figure BDA0002968280140000025
Costs for electro-hydrogen production and methane synthesis, respectively; f (f) 3 The environmental cost for the operation of the comprehensive energy system;
Figure BDA0002968280140000026
The electric transmission loss of the electricity storage equipment;
Figure BDA0002968280140000027
Gas transmission loss for the hydrogen storage device;
Figure BDA0002968280140000028
The gas transmission loss of the gas storage equipment; f (f) 5 To sell the benefits generated by hydrogen.
Further, the method comprises the steps of,
Figure BDA0002968280140000029
Figure BDA00029682801400000210
wherein T is totalNumber of time periods, N TH Is the number of the thermal power generating units,
Figure BDA00029682801400000211
the active output of the ith thermal power generating unit at the moment t, a i ,b i ,c i Is the cost coefficient of the ith thermal power unit;
c GT for the unit operating cost of the gas turbine, N GT Is the number of the gas turbines to be used,
Figure BDA00029682801400000212
for the i-th gas turbine active power output at time t, mu GT The power generation efficiency of the gas turbine is that deltat is the interval of time;
c FC for the unit operation cost of the fuel cell, N FC Is the number of the fuel cells to be used,
Figure BDA0002968280140000031
for the i-th fuel cell active force at time t, mu H2E Is the power generation efficiency of the fuel cell.
Further, the method comprises the steps of,
Figure BDA0002968280140000032
Figure BDA0002968280140000033
Figure BDA0002968280140000034
wherein ,cB For the unit operation cost of the electricity storage equipment, S (t-delta t) is the energy storage charge state, eta ch 、η dch The charge and discharge state of the electricity storage equipment at the moment t is P ch,t 、P dch,t For the charge and discharge power of the electricity storage equipment at the moment of t, U B For terminal voltage of electric storage equipment, C B Is an electric storage deviceIs a capacity of (2);
c G for the unit operation cost of the gas storage equipment E G,t-1 The residual electric quantity of the gas energy storage system at the time t-1,
Figure BDA0002968280140000035
inputting natural gas power of gas storage equipment at time t, P t G Power of natural gas consumed by power generation of gas turbine at t moment mu H2G Is the synthesis efficiency of natural gas, mu GT N is the power generation efficiency of the gas turbine H2G For the number of methane synthesis units, +.>
Figure BDA0002968280140000036
For the hydrogen energy power consumed by the ith methane synthesizing device at the time t, N GT For the number of gas turbines>
Figure BDA0002968280140000037
The active output of the ith gas turbine at the time t is obtained; />
c H For the unit operation cost of the hydrogen storage device E H,t-1 The residual electric quantity of the hydrogen energy storage system at the time t-1,
Figure BDA0002968280140000038
inputting hydrogen energy power of a hydrogen energy storage system at the moment t, P t H The power mu of hydrogen energy consumed by the fuel cell at the moment t E2H For the efficiency of electrical hydrogen production, mu H2E N is the power generation efficiency of the fuel cell EL For the number of electrolytic cells>
Figure BDA0002968280140000039
For the active power consumed by the ith electrolytic cell at time t, N FC For the number of fuel cells>
Figure BDA00029682801400000310
Is the active force of the ith fuel cell at time t.
Further, the method comprises the steps of,
Figure BDA00029682801400000311
wherein ,cEL For the running cost of the electrolytic cell per unit time, N EL For the number of cells to be used,
Figure BDA00029682801400000312
active power consumed by the ith electrolytic cell at the time t; lambda (lambda) H2G For the unit time operation cost of the methane synthesis equipment, N H2G For the number of methane synthesis plants, +.>
Figure BDA0002968280140000041
The active power consumed by the ith methane synthesis apparatus at time t.
Further, the method comprises the steps of,
Figure BDA0002968280140000042
wherein ,
Figure BDA0002968280140000043
is the carbon emission coefficient of the thermal generator set, +.>
Figure BDA0002968280140000044
For the carbon emission coefficient of the gas turbine, < > for>
Figure BDA0002968280140000045
For the carbon dioxide utilization of the integrated energy system in the methane synthesis process, < > for>
Figure BDA0002968280140000046
Is the carbon emission coefficient of the traditional fuel automobile, N EBUS N is the total number of electric vehicles in the comprehensive energy system HBUS The total number of hydrogen energy automobiles in the comprehensive energy system is calculated;
Figure BDA0002968280140000047
For the power of the ith electric car at time t, < >>
Figure BDA0002968280140000048
The power at t moment of the ith hydrogen energy automobile.
Further, the method comprises the steps of,
Figure BDA0002968280140000049
Figure BDA00029682801400000410
Figure BDA00029682801400000411
wherein ,
Figure BDA00029682801400000412
for electric transmission loss, k of electric storage equipment es Loss coefficient, mu, of electric transmission ch 、μ dch For the charge and discharge efficiency of the electricity storage device, < >>
Figure BDA00029682801400000413
The charge and discharge power of the electricity storage equipment at the moment t;
Figure BDA00029682801400000414
for gas transmission loss, k of gas storage apparatus gs Loss factor for gas transmission, +.>
Figure BDA00029682801400000415
Air inlet and air discharge efficiency of the hydrogen storage device, < >>
Figure BDA00029682801400000416
The air inlet and air discharge rates of the air storage equipment at the moment t are respectively;
Figure BDA00029682801400000417
for gas transmission loss, k of hydrogen storage device hs Loss factor for hydrogen transport, +.>
Figure BDA00029682801400000418
For the air intake and air discharge efficiency of the hydrogen storage device, < >>
Figure BDA00029682801400000419
The air inlet and air discharge rates of the hydrogen storage equipment at the moment t are respectively. />
Further, the method comprises the steps of,
Figure BDA00029682801400000420
wherein ,f5 To sell the benefits of hydrogen, c SH The price is sold for the unit power of hydrogen,
Figure BDA00029682801400000421
is the active power at time t of the ith hydrogen-oxygen fuel cell, +.>
Figure BDA00029682801400000422
The power at t moment of the ith hydrogen energy automobile.
Further, the device operational constraints include: thermal power generating unit, gas turbine, fuel cell, electrolysis trough, storage of electricity equipment, storage of hydrogen equipment and the operation constraint of storage of gas equipment.
Further, the coupling constraint between the subsystems comprises a coupling constraint between a natural gas system and a hydrogen energy storage system, a coupling constraint between the natural gas system and an electric power system and a coupling constraint between the hydrogen energy storage system and the electric power system;
the coupling constraint between the natural gas system and the hydrogen energy storage system is as follows:
Figure BDA0002968280140000051
the coupling constraint between the natural gas system and the power system is as follows:
Figure BDA0002968280140000052
The coupling constraint between the hydrogen energy storage system and the power system is as follows:
Figure BDA0002968280140000053
in the formula ,
Figure BDA0002968280140000054
the gas injection amount θ at time t of the ith methane synthesizer gas Is natural gas with high heat value and mu H2G Efficiency of methane synthesis for hydrogen, +.>
Figure BDA0002968280140000055
The methane synthesis device is injected with the methane power value of the natural gas system at time t, < >>
Figure BDA0002968280140000056
The hydrogen energy power flowing into the methane synthesizing device at the moment t for the hydrogen energy storage system;
Figure BDA0002968280140000057
for the active output of the ith micro gas turbine at the time t, mu GT Is the coefficient of the power generation efficiency of the gas turbine,
Figure BDA0002968280140000058
for the natural gas consumption quantity of the ith gas turbine at the time t, theta gas Is natural gas with high heat value;
Figure BDA0002968280140000059
is the active force of the ith fuel cell at the moment t and mu H2E For the fuel cell power generation efficiency coefficient, < >>
Figure BDA00029682801400000510
And inputting the hydrogen energy power of the fuel cell into the hydrogen energy storage system t at the moment.
The invention provides a comprehensive energy system operation optimization method considering hydrogen energy, which comprises the following steps: aiming at a comprehensive energy system considering hydrogen energy, the comprehensive energy system comprises an electric power system, a natural gas system and a hydrogen energy storage system, wherein equipment operation constraint, electric power network constraint and coupling constraint among subsystems are taken as operation constraint conditions, and minimized operation cost is taken as an objective function, and an operation optimization model of the comprehensive energy system is established;
inputting operation parameters of the comprehensive energy system into the operation optimization model, and solving the operation optimization model to obtain an optimized operation scheme of the comprehensive energy system; the operation parameters of the comprehensive energy system comprise system structure parameters, load parameters, equipment operation parameters, energy transmission loss parameters, environment cost parameters and benefit parameters of selling hydrogen.
The technical scheme of the invention can achieve the following beneficial effects:
(1) The comprehensive energy system structure taking the hydrogen energy into consideration comprises a plurality of energy conversion processes among the hydrogen energy system, an electric power system and a natural gas system, and considers the environmental benefit and the economic benefit of new energy loads of hydrogen energy automobiles, electric automobiles and the like in the comprehensive energy system;
(2) In the comprehensive energy system operation optimization model considering the hydrogen energy, the energy coupling constraint between the hydrogen energy system and the power system, between the hydrogen energy system and the natural gas system and between the power system and the natural gas system is fully considered, and the optimization model is more comprehensive and reasonable;
(3) In the integrated energy system operation optimization model considering hydrogen energy, which is provided by the invention, the operation cost of various devices in the system, the energy conversion (electric hydrogen production process and methane synthesis process) cost among subsystems, the environmental cost, the energy transmission loss cost and the relational expression of benefits brought by selling hydrogen energy are described in detail, and the obtained integrated energy optimization operation scheme is more practical.
Drawings
FIG. 1 is a flow chart of a method for optimizing operation of a comprehensive energy system taking into account hydrogen energy in accordance with the present invention;
FIG. 2 is a diagram of the integrated energy system of the present invention, which is a method for optimizing the operation of the integrated energy system in consideration of hydrogen energy;
fig. 3 is a schematic structural diagram of an integrated energy system operation optimization method considering hydrogen energy according to an embodiment of the present invention.
Detailed Description
The embodiment of the invention provides a comprehensive energy system operation optimization method considering hydrogen energy, which aims to solve the technical problems that the hydrogen energy system is not considered to participate in operation and the operation scheme of the comprehensive energy system is inconsistent with actual application in the related technology.
The present invention will be described more fully hereinafter with reference to the accompanying drawings, in which the invention is shown. Preferred embodiments of the present invention are shown in the drawings. This invention may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete.
Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. The terminology used herein in the description of the invention is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. The term "and/or" as used herein includes any and all combinations of one or more of the associated listed items.
Referring to fig. 1, the following is an embodiment of an operation optimization method of an integrated energy system considering hydrogen energy according to the present invention, including:
aiming at a comprehensive energy system considering hydrogen energy, the comprehensive energy system comprises an electric power system, a natural gas system and a hydrogen energy storage system, wherein equipment operation constraint, electric power network constraint and coupling constraint among subsystems are taken as operation constraint conditions, and minimized operation cost is taken as an objective function, and an operation optimization model of the comprehensive energy system is established;
inputting operation parameters of the comprehensive energy system into the operation optimization model, and solving the operation optimization model to obtain an optimized operation scheme of the comprehensive energy system; the operation parameters of the comprehensive energy system comprise system structure parameters, load parameters, equipment operation parameters, energy transmission loss parameters, environment cost parameters and benefit parameters of selling hydrogen.
Referring to fig. 2, the integrated energy system according to the embodiment of the present invention includes three subsystems, i.e., an electric power system, a natural gas system, and a hydrogen energy storage system (i.e., a hydrogen system). The power generation equipment in the power system mainly comprises a traditional thermal generator set, a gas turbine and the like, and is used for meeting the electric load requirements in the comprehensive energy system; the main loads of the power system comprise power loads, energy storage equipment, electric automobiles and the like. The natural gas system mainly comprises a natural gas network, can directly provide gas for a gas load, and can also be used as gas input of a gas turbine; part of carbon dioxide generated by operation of a traditional thermal power generating unit, a gas generating unit and the like in the electric power system can be reacted with hydrogen to synthesize methane gas through catalysis, and the methane gas is input into a natural gas system. The hydrogen system is mainly hydrogen storage equipment, can directly provide hydrogen for hydrogen energy loads such as hydrogen energy automobiles and the like, and can also be directly used as the input of an oxyhydrogen fuel cell; an electrolytic tank in the electric power system generates hydrogen by electrolyzing water and inputs the hydrogen into the hydrogen system.
According to the comprehensive energy system operation optimization method considering the hydrogen energy, an operation optimization model of the comprehensive energy system considering the hydrogen energy is established, and an objective function of the model is to minimize the operation cost of the comprehensive energy system, wherein the model mainly comprises equipment operation cost (a traditional thermal power unit, a gas turbine, a fuel cell, electric energy storage, hydrogen storage equipment and gas storage equipment), energy conversion cost (an electric hydrogen production process and a methane synthesis process), environmental cost and transmission loss cost (electric transmission loss of the electric storage equipment, gas transmission loss of the hydrogen storage equipment and the gas storage equipment) and benefits brought by selling hydrogen to oxyhydrogen fuel cells and hydrogen energy automobiles.
Specifically, the objective function of running the optimization model is as shown in formulas (1) to (4):
z=min(f 1 +f 2 +f 3 +f 4 -f 5 ) (1)
Figure BDA0002968280140000071
Figure BDA0002968280140000072
Figure BDA0002968280140000073
wherein z is the running cost of the integrated energy system; f (f) 1 、f 2 、f 3 、f 4 、f 5 Respectively representing the running cost of equipment, the energy conversion cost, the environmental cost, the transmission loss cost and the income brought by selling hydrogen;
in the formula ,
Figure BDA0002968280140000081
the running cost of the thermal power generating unit is;
Figure BDA0002968280140000082
The operating cost of the gas turbine;
Figure BDA0002968280140000083
The fuel cell operating cost;
Figure BDA0002968280140000084
The running cost for electric energy storage;
Figure BDA0002968280140000085
The operation cost of the hydrogen storage equipment is;
Figure BDA0002968280140000086
The operation cost of the gas storage equipment is;
Figure BDA0002968280140000087
The cost of hydrogen production by electricity is;
Figure BDA0002968280140000088
Is methane synthesisCost is achieved; f (f) 3 The running environment cost of the comprehensive energy system is realized;
Figure BDA0002968280140000089
The electric transmission loss for the electricity storage equipment;
Figure BDA00029682801400000810
The gas transmission loss of the hydrogen storage equipment;
Figure BDA00029682801400000811
The gas transmission loss is the gas storage equipment; f (f) 5 To sell the benefits of hydrogen.
(1) Cost of operation of equipment of integrated energy system
For the traditional thermal generator set, considering the running cost, the corresponding expression is shown in formula (5):
Figure BDA00029682801400000812
wherein T is the total time period number; n (N) TH The number of thermal power generating units;
Figure BDA00029682801400000813
the active output of the ith thermal power generating unit at the time t is obtained; a, a i ,b i ,c i Is the cost coefficient of the ith thermal power unit; Δt is the interval of the period.
It should be noted that, strictly speaking, t in the present embodiment should be a period, and the time t in the present embodiment refers to the start time of the corresponding period, and considering that the interval of the general period is relatively short, the output of the generator is considered to be unchanged in the period, and the time t refers to the corresponding period, so as to conform to the general description habit: and outputting the XX thermal power unit at the XX moment.
For the gas turbine, the operation cost is related to the active power and the power generation efficiency, and the corresponding expression is shown in the formula (6):
Figure BDA00029682801400000814
in the formula ,cGT The unit operation cost of the gas turbine; n (N) GT The number of the gas turbines;
Figure BDA00029682801400000815
the active output of the ith gas turbine at the time t is obtained; mu (mu) GT Is the power generation efficiency of the gas turbine.
For the fuel cell, the operation cost is related to the active power and the power generation efficiency, and the corresponding expression is shown in the formula (7):
Figure BDA00029682801400000816
in the formula ,cFC The unit operation cost of the fuel cell; n (N) FC The number of the fuel cells;
Figure BDA0002968280140000091
an active force of the ith fuel cell at the time t; mu (mu) H2E Is the power generation efficiency of the fuel cell. />
For the electricity storage equipment, the operation cost can be obtained according to the charge state, the charge and discharge power, the terminal voltage and the capacity, and the corresponding expression is shown in the formula (8):
Figure BDA0002968280140000092
in the formula ,cB The running cost of the electric energy storage unit; s (t-delta t) is the energy storage charge state; η (eta) ch 、η dch The charging and discharging state of the electric energy storage at the time t; p (P) ch,t 、P dch,t Charging and discharging power for the electric energy storage at the time t; u (U) B A terminal voltage for electrical energy storage; c (C) B Is an electrical energy storage capacity.
Regarding the gas storage equipment, considering the gas storage capacity at the previous moment, the synthesis amount of natural gas (methane) at the current moment and the natural gas consumption amount of the gas turbine, the operation cost expression is shown in a formula (9):
Figure BDA0002968280140000093
in the formula ,cG The unit operation cost of the gas storage equipment; e (E) G,t-1 The energy storage electric quantity is the gas at the time t-1;
Figure BDA0002968280140000094
inputting natural gas power of the gas storage equipment at the time t; p (P) t G Natural gas power consumed by the power generation of the gas turbine at the time t; mu (mu) H2G Is the synthesis efficiency of natural gas (methane); mu (mu) GT Generating efficiency for the gas turbine; n (N) H2G The number of the methane synthesis devices;
Figure BDA0002968280140000095
The hydrogen energy power consumed by the ith methane synthesizing device at the time t is obtained; n (N) GT The number of the gas turbines;
Figure BDA0002968280140000096
Is the active output of the ith gas turbine at the time t.
Regarding the hydrogen storage equipment, considering the hydrogen storage amount at the last moment, the electric hydrogen production amount in the electrolytic tank at the current moment and the hydrogen consumption amount of the fuel cell, the operation cost expression is obtained as shown in a formula (10):
Figure BDA0002968280140000097
in the formula ,cH The unit operation cost for hydrogen energy storage; e (E) H,t-1 The hydrogen energy storage electric quantity at the time t-1;
Figure BDA0002968280140000098
inputting hydrogen energy power of a hydrogen energy storage system at the time t; p (P) t H The hydrogen energy power consumed by the fuel cell power generation at the time t; mu (mu) E2H The hydrogen production efficiency is achieved by electricity; mu (mu) H2E Generating efficiency for the fuel cell; n (N) EL The number of the electrolytic cells;
Figure BDA0002968280140000099
Active power consumed by the ith electrolytic cell at the time t; n (N) FC Is the number of fuel cells;
Figure BDA00029682801400000910
Is the active force of the ith fuel cell at time t.
(2) Energy conversion cost of integrated energy system
Considering that the hydrogen energy in the future is gradually increased in the comprehensive energy system, the cost in the electric hydrogen production process cannot be ignored. The cost of the electric hydrogen production process is mainly generated by an electrolytic cell, and the expression of the electric hydrogen production cost is shown in a formula (11):
Figure BDA0002968280140000101
in the formula ,cEL The running cost of the electrolytic cell per unit time; n (N) EL The number of the electrolytic cells;
Figure BDA0002968280140000102
the active power consumed by the ith electrolytic cell at time t.
In order to reduce the carbon emission of the comprehensive energy system, hydrogen and carbon dioxide are utilized to catalyze and react to synthesize methane. The methane synthesis costs are mainly related to the operation of the methane synthesis plant, and the methane synthesis costs are shown in formula (12):
Figure BDA0002968280140000103
in the formula ,λH2G The unit time operation cost of the methane synthesis equipment is set; n (N) H2G The number of methane synthesis equipment;
Figure BDA0002968280140000104
for the ith methane synthesis apparatusActive power consumed at time t.
(3) Environmental cost of integrated energy system
One of the main objectives of the integrated energy system is to reduce carbon emissions and increase the utilization of green energy, so that the environmental costs of the integrated energy system are necessarily considered. In the operation optimization model provided by the embodiment, the main carbon emission of the comprehensive energy system is derived from the traditional thermal generator set and the gas turbine, and the synthesis of methane, the development of hydrogen energy automobiles and electric automobiles are main ways for reducing the carbon emission of the comprehensive energy. The environmental cost expression is shown in formula (13):
Figure BDA0002968280140000105
in the formula ,
Figure BDA0002968280140000106
the carbon emission coefficient of the thermal generator set;
Figure BDA0002968280140000107
Carbon emission coefficient of the gas turbine;
Figure BDA0002968280140000108
The utilization rate of carbon dioxide in the comprehensive energy system in the methane synthesis process is utilized;
Figure BDA0002968280140000109
The carbon emission coefficient of the traditional fuel automobile; n (N) EBUS The total number of the electric automobiles in the comprehensive energy system; n (N) HBUS The total number of hydrogen energy automobiles in the comprehensive energy system is calculated;
Figure BDA0002968280140000111
The power of the ith electric automobile at the moment t;
Figure BDA0002968280140000112
The power at t moment of the ith hydrogen energy automobile.
(4) Transmission loss cost of integrated energy system
The electrical loss expression of the electricity storage device is shown in formula (14):
Figure BDA0002968280140000113
in the formula ,
Figure BDA0002968280140000114
the electric transmission loss for the electricity storage equipment; k (k) es Is the electric transmission loss coefficient; mu (mu) ch 、μ dch The charge and discharge efficiency of the electricity storage equipment is improved;
Figure BDA0002968280140000115
And storing energy to charge and discharge power at the moment t.
The gas transmission loss expression of the gas storage device is shown in formula (15):
Figure BDA0002968280140000116
in the formula ,
Figure BDA0002968280140000117
the gas transmission loss is the gas storage equipment; k (k) gs Is the gas transmission loss coefficient;
Figure BDA0002968280140000118
The air inlet and air discharge efficiency of the hydrogen storage equipment;
Figure BDA0002968280140000119
The air inlet rate and the air outlet rate of the air storage equipment at the moment t are respectively.
The gas transmission loss expression of the hydrogen storage device is shown in formula (16):
Figure BDA00029682801400001110
in the formula ,
Figure BDA00029682801400001111
the gas transmission loss of the hydrogen storage equipment; k (k) hs Is the hydrogen transmission loss coefficient;
Figure BDA00029682801400001112
The air inlet and air discharge efficiency of the hydrogen storage equipment;
Figure BDA00029682801400001113
The air inlet rate and the air discharge rate of the hydrogen storage device at the moment t are respectively.
(5) Benefit brought by selling hydrogen by comprehensive energy system
The profitability expression for hydrogen gas to be sold to hydrogen-oxygen fuel cells and hydrogen energy automobiles is shown in formula (17):
Figure BDA00029682801400001114
in the formula ,f5 To sell hydrogen revenue; c SH Selling prices for hydrogen per unit of power;
Figure BDA00029682801400001115
active power at the time of the ith hydrogen-oxygen fuel cell t;
Figure BDA0002968280140000121
The power at t moment of the ith hydrogen energy automobile.
The operation constraint of the comprehensive energy system operation optimization method considering the hydrogen energy provided by the embodiment comprises equipment operation constraint, power network constraint and coupling constraint among subsystems.
Wherein, (1) the equipment operation constraint comprises the operation constraint of a traditional thermal power generating unit, a gas turbine, a fuel cell and an electrolytic tank and the operation constraint of energy storage equipment (including electricity storage equipment, gas storage equipment and hydrogen storage equipment);
and the operation constraint of the traditional thermal power generating unit, the gas turbine, the fuel cell and the electrolytic tank comprises the upper limit constraint and the lower limit constraint of active power and the upper limit and the lower limit constraint of climbing rate.
The operation constraint of the thermal power generating unit is shown in formulas (18) and (19):
Figure BDA0002968280140000122
Figure BDA0002968280140000123
in the formula ,
Figure BDA0002968280140000124
is the upper limit and the lower limit of the active power of the ith thermal power unit t at the moment, ++>
Figure BDA0002968280140000125
The upper limit and the lower limit of the climbing rate of the ith thermal power generating unit.
The gas turbine operating constraints are shown in formulas (20), (21):
Figure BDA0002968280140000126
Figure BDA0002968280140000127
in the formula ,
Figure BDA0002968280140000128
for the upper and lower limit of the active power of the ith gas turbine t at time, +.>
Figure BDA0002968280140000129
Is the upper and lower limits of the ramp rate of the ith gas turbine.
The fuel cell operation constraints are shown in formulas (22), (23):
Figure BDA00029682801400001210
Figure BDA00029682801400001211
in the formula ,
Figure BDA00029682801400001212
upper and lower limits of active power for the ith fuel cell t>
Figure BDA00029682801400001213
Is the upper and lower limits of the ramp rate of the ith fuel cell.
The cell operating constraints are shown in equation (24):
Figure BDA00029682801400001214
in the formula ,
Figure BDA00029682801400001215
is the upper and lower limits of the active power of the ith electrolytic tank t at the moment.
The energy storage equipment of the comprehensive energy system is divided into electricity storage equipment, gas storage equipment and hydrogen storage equipment, and the operation constraint of the energy storage equipment comprises energy balance constraint, energy upper and lower limit constraint, power upper and lower limit constraint and transmission loss constraint.
The energy balance constraint of the electricity storage device is shown in formula (25):
Figure BDA0002968280140000131
in the formula ,
Figure BDA0002968280140000132
the electricity storage quantity is t time; mu (mu) ch 、μ dch The charge and discharge efficiency of the electricity storage equipment is improved;
Figure BDA0002968280140000133
And storing energy to charge and discharge power for the time t.
The upper and lower limit constraints of the electric storage amount are shown in the formula (26):
Figure BDA0002968280140000134
in the formula ,EB,min 、E B,max Is the upper and lower limit of the electricity storage capacity.
The upper and lower limit constraints of the charge and discharge power of the electricity storage device are shown in formulas (27) and (28):
Figure BDA0002968280140000135
Figure BDA0002968280140000136
in the formula ,Pch,min 、P ch,max Charging the energy storage device with an upper and lower power limit; p (P) dch,min 、P dch,max And the upper limit and the lower limit of the power are amplified for the energy storage device.
The hydrogen amount balance constraint of the gas storage device is shown in a formula (29):
Figure BDA0002968280140000137
in the formula ,Vt G The hydrogen storage quantity is the time t,
Figure BDA0002968280140000138
the air inlet and air discharge efficiency of the hydrogen storage equipment;
Figure BDA0002968280140000139
Figure BDA00029682801400001310
The air inlet rate and the air outlet rate of the air storage equipment at the moment t are respectively.
The upper and lower limit constraints of the gas storage volume are shown in a formula (30):
Figure BDA00029682801400001311
in the formula ,VG,min 、V G,max Is the upper and lower limit of the gas storage amount.
The hydrogen amount balance constraint of the hydrogen storage device is as shown in formula (31):
Figure BDA00029682801400001312
in the formula ,Vt H The hydrogen storage quantity is the time t,
Figure BDA00029682801400001313
the air inlet and air discharge efficiency of the hydrogen storage equipment;
Figure BDA00029682801400001314
Figure BDA00029682801400001315
The air inlet rate and the air discharge rate of the hydrogen storage device at the moment t are respectively.
The upper and lower limit constraints of the hydrogen storage amount are shown in the formula (32):
Figure BDA00029682801400001316
in the formula ,Vt H For the hydrogen storage quantity at time t, V H,min 、V H,max Is the upper and lower limit of the hydrogen storage amount.
(2) Power network constraints for integrated energy systems
The power network flow constraint adopts a linear model, and the line loss is ignored. The linear power flow constraint expressions are shown in formulas (33) to (37):
Figure BDA0002968280140000141
Figure BDA0002968280140000142
Figure BDA0002968280140000143
Figure BDA0002968280140000144
v i,min ≤v i,t ≤v i,max (37)
wherein epsilon is a line set; s is S ij,t ,S ki,t The power flow of the line between the nodes i and j and the nodes k and i at the moment t; s is(s) i,t The injection power of the node i at the moment t;
Figure BDA0002968280140000145
the power of the generator at the node i at the moment t;
Figure BDA0002968280140000146
Rated power of the node i generator;
Figure BDA0002968280140000147
the load power at the node i comprises the load charging power of the electric automobile; re (·) is the real part of the acquired complex; v i,t 、v j,t The squares of the voltage amplitude values at the node i and the node j at the moment t are respectively obtained; v i,min 、v i,max The maximum value and the minimum value of the square of the voltage amplitude at the node i are obtained;
Figure BDA0002968280140000148
Transpose the conjugate of the impedance of the line ij.
(3) Inter-subsystem coupling constraints for integrated energy systems
The natural gas network and the hydrogen network are coupled through the methane synthesis device, and the coupling constraint is the relation between the output power of the hydrogen energy system and the natural gas injection power, as shown in formula (38):
Figure BDA0002968280140000149
in the formula ,
Figure BDA00029682801400001410
the gas injection amount at the time t of the ith methane synthesis device; θ gas Is natural gas with high heat value; mu (mu) H2G The efficiency of synthesizing methane from hydrogen;
Figure BDA00029682801400001411
Injecting a methane power value of a natural gas network into the methane synthesis device at the moment t;
Figure BDA00029682801400001412
And (5) feeding hydrogen energy power into the methane synthesis device at the moment of the hydrogen energy storage system t.
The natural gas network and the electric power network are coupled through the gas turbine, and the coupling constraint is the relation between the output power of the micro gas turbine and the natural gas consumption, as shown in a formula (39):
Figure BDA00029682801400001413
in the formula ,
Figure BDA00029682801400001414
the active power output of the ith micro gas turbine at the t moment is obtained; mu (mu) GT The power generation efficiency coefficient of the gas turbine is;
Figure BDA00029682801400001415
The natural gas consumption of the ith gas turbine at the time t is calculated; θ gas Is natural gas with high heat value.
The hydrogen network and the electric network are coupled through an oxyhydrogen fuel cell and an electrolytic tank, wherein the relation between the output power of the oxyhydrogen fuel cell and the hydrogen injection amount is as shown in a formula (40):
Figure BDA0002968280140000151
in the formula ,
Figure BDA0002968280140000152
the active force of the ith fuel cell at the t moment is obtained; mu (mu) H2E The power generation efficiency coefficient of the fuel cell;
Figure BDA0002968280140000153
and inputting the hydrogen energy power of the fuel cell into the hydrogen energy storage system t at the moment.
The relationship between cell power consumption and the amount of hydrogen produced is shown in equation (41):
Figure BDA0002968280140000154
in the formula ,
Figure BDA0002968280140000155
the active force of the ith electrolytic tank at the time t is obtained; mu (mu) E2H The hydrogen production efficiency is achieved by electricity;
Figure BDA0002968280140000156
Power is consumed at time t for the ith electrolytic cell.
The solving process of the operation optimization model of the embodiment is as follows: inputting operation parameters of the comprehensive energy system into the operation optimization model, and solving the operation optimization model to obtain an optimized operation scheme of the comprehensive energy system; the method specifically comprises three parts of data input, optimization software solving and operation scheme output:
(1) Data input: the operation parameters of the comprehensive energy system are input, including system structure parameters, load (electric load, gas load and hydrogen load) parameters, equipment operation parameters, energy transmission loss parameters, environment cost parameters, hydrogen sales gain parameters and the like.
(2) And (3) solving by using optimization software: the operation optimization model of the comprehensive energy system provided by the embodiment of the invention is a linear optimization model, and the solution of the optimization model can be realized by utilizing mature commercial optimization software programming.
(3) And (3) outputting an operation scheme: and obtaining the active output strategies of all the power supplies (thermal power generating units, gas turbines and fuel cells) at all the moments according to the solving result of the operation optimization model, charging and discharging power of the electric energy storage at all the moments, and air inlet and air outlet rates of the air/hydrogen energy storage equipment at all the moments, so as to obtain the economic operation scheme of the comprehensive energy system.
According to the comprehensive energy system operation optimization method considering the hydrogen energy, the energy conversion relation and the operation scene of the comprehensive energy system with the participation of the hydrogen energy are fully considered, the operation coupling relation among all subsystems of the hydrogen energy system, the natural gas system and the electric power system is fully considered, the energy conversion relation among the electric power system, the natural gas system and the hydrogen energy system is given, the operation cost of various devices of the comprehensive energy system, the electricity-natural gas-hydrogen energy conversion cost, the environmental benefit and the economic benefit brought by the participation of the hydrogen energy are taken as optimization targets, a comprehensive energy economic operation optimization model with the participation of the hydrogen energy is established, and the optimal operation mode of the comprehensive energy system with the participation of the hydrogen energy is obtained by solving the operation optimization model.
Referring to fig. 3, the following is another embodiment of an operation optimization method of a comprehensive energy system considering hydrogen energy, which is a comprehensive energy system comprising an electric power system, a natural gas system, a hydrogen storage device, etc., and a structural schematic diagram of the comprehensive energy system containing hydrogen for energy storage is shown in fig. 3.
Two thermal generators at nodes 1 and 2, two gas turbines at nodes 4 and 7 and a fuel cell at node 3 in the power system are main power generation equipment. The electric load comprises a common electric load at each node, an electrolytic tank and an electric automobile. The electricity storage device may be used as both a power generation device and a load. The node N1 in the natural gas system is connected to a natural gas source, the node N3 is a gas storage tank, and the rest nodes have natural gas loads. The hydrogen system mainly comprises hydrogen storage equipment, and the main load is an oxyhydrogen fuel cell and a hydrogen energy automobile.
The power system is coupled to the hydrogen system via an oxyhydrogen fuel cell at node 3 and an electrolyzer. The hydrogen system injects hydrogen into the fuel cell to produce electric energy for injection into the power system. The electrolyzer is used for injecting hydrogen into the hydrogen system through electrolysis of water.
The power system is coupled to the natural gas system by a gas turbine at node 7. The injection of natural gas from the natural gas system into the combustion turbine reflects the generation of electrical energy into the electrical power system.
The natural gas system and the hydrogen system are coupled by a methane synthesis unit at node N3. And the hydrogen input by the hydrogen system and the captured carbon dioxide are catalyzed and reflected to obtain methane gas.
The power system, the natural gas system and the hydrogen energy storage system in the comprehensive energy system realize energy conversion and operation among all subsystems of the comprehensive energy system through the coupling elements.
In the method for optimizing the operation of the integrated energy system taking the hydrogen energy into consideration provided by the embodiment of the invention, in the structure of the integrated energy system taking the hydrogen energy into consideration,
the method comprises various energy conversion processes among the hydrogen energy system, the electric power system and the natural gas system, and considers the environmental benefit and the economic benefit of new energy loads such as hydrogen energy automobiles, electric automobiles and the like in the comprehensive energy system;
the energy coupling constraint between the hydrogen energy system and the power system, between the hydrogen energy system and the natural gas system and between the power system and the natural gas system are fully considered, and the optimization model is more comprehensive and reasonable;
the method has the advantages that the operation cost of various devices in the system, the energy conversion (electric hydrogen production process and methane synthesis process) cost among subsystems, the environmental cost, the energy transmission loss cost and the relational expression of benefits brought by hydrogen selling energy are described in detail, and the obtained comprehensive energy optimization operation scheme is more practical.
It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the above-described systems, apparatuses and units may refer to corresponding procedures in the foregoing method embodiments, which are not repeated herein.
In the embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other manners. For example, the apparatus embodiments described above are merely illustrative, e.g., the division of the units is merely a logical function division, and there may be additional divisions when actually implemented, e.g., multiple units or components may be combined or integrated into another system, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, which may be in electrical, mechanical or other form.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional units.
The integrated units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied essentially or in part or all of the technical solution or in part in the form of a software product stored in a storage medium, including instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
The above embodiments are only for illustrating the technical solution of the present invention, and are not limiting; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims (8)

1. An integrated energy system operation optimization method considering hydrogen energy is characterized by comprising the following steps:
aiming at a comprehensive energy system considering hydrogen energy, the comprehensive energy system comprises an electric power system, a natural gas system and a hydrogen energy storage system, wherein equipment operation constraint, electric power network constraint and coupling constraint among subsystems are taken as operation constraint conditions, and minimized operation cost is taken as an objective function, and an operation optimization model of the comprehensive energy system is established;
inputting operation parameters of the comprehensive energy system into the operation optimization model, and solving the operation optimization model to obtain an optimized operation scheme of the comprehensive energy system; the operation parameters of the comprehensive energy system comprise system structure parameters, load parameters, equipment operation parameters, energy transmission loss parameters, environment cost parameters and benefit parameters of selling hydrogen;
the objective function is:
z=min(f 1 +f 2 +f 3 +f 4 -f 5 );
wherein z is the running cost of the integrated energy system; f (f) 1 、f 2 、f 3 、f 4 、f 5 Respectively represent the running cost of the equipment,Energy conversion cost, environmental cost, transmission loss cost, and revenue from selling hydrogen;
f 1 =f 1 TH +f 1 GT +f 1 FC +f 1 B +f 1 G +f 1 H
Figure FDA0004121969510000011
f 1 TH 、f 1 GT 、f 1 FC the operation cost of the thermal power generating unit, the gas turbine and the fuel cell is respectively; f (f) 1 B 、f 1 H 、f 1 G The operation cost of the electricity storage equipment, the hydrogen storage equipment and the gas storage equipment is as follows;
Figure FDA0004121969510000012
costs for electro-hydrogen production and methane synthesis, respectively;
Figure FDA0004121969510000013
The electric transmission loss of the electricity storage equipment;
Figure FDA0004121969510000014
Gas transmission loss for the hydrogen storage device;
Figure FDA0004121969510000015
The gas transmission loss of the gas storage equipment;
Figure FDA0004121969510000016
wherein T is the total time period number, N TH Is the number of the thermal power generating units,
Figure FDA0004121969510000017
the active output of the ith thermal power generating unit at the moment t, a i ,b i ,c i Is the cost coefficient of the ith thermal power unit;
c GT for the unit operating cost of the gas turbine, N GT Is the number of the gas turbines to be used,
Figure FDA0004121969510000018
for the i-th gas turbine active power output at time t, mu GT The power generation efficiency of the gas turbine is that deltat is the interval of time;
c FC for the unit operation cost of the fuel cell, N FC Is the number of the fuel cells to be used,
Figure FDA0004121969510000019
is the active force of the ith hydrogen-oxygen fuel cell at the moment t, mu H2E Is the power generation efficiency of the fuel cell.
2. The method for optimizing the operation of a comprehensive energy system taking into account hydrogen energy according to claim 1, wherein,
Figure FDA0004121969510000021
Figure FDA0004121969510000022
Figure FDA0004121969510000023
wherein ,cB For the unit operation cost of the electricity storage equipment, S (t-delta t) is the energy storage charge state, eta ch 、η dch The charge and discharge state of the electricity storage equipment at the moment t is P ch,t 、P dch,t For the charge and discharge power of the electricity storage equipment at the moment of t, U B For terminal voltage of electric storage equipment, C B Is the capacity of the electricity storage device;
c G for the unit operation cost of the gas storage equipment E G,t-1 The residual electric quantity of the gas energy storage system at the time t-1 is P t E2G Inputting natural gas power of gas storage equipment at time t, P t G Power of natural gas consumed by power generation of gas turbine at t moment mu H2G For the efficiency of synthesizing methane from hydrogen, N H2G For the number of methane synthesis units,
Figure FDA0004121969510000024
for the active power consumed by the ith methane synthesizer at time t, N GT The number of the gas turbines;
c H for the unit operation cost of the hydrogen storage device E H,t-1 The residual electric quantity of the hydrogen energy storage system at the time t-1 is P t E2H Inputting hydrogen energy power of a hydrogen energy storage system at the moment t, P t H The power mu of hydrogen energy consumed by the fuel cell at the moment t E2H For the efficiency of electric hydrogen production, N EL For the number of cells to be used,
Figure FDA0004121969510000025
for the active power consumed by the ith electrolytic cell at time t, N FC Is the number of fuel cells.
3. The method for optimizing the operation of a comprehensive energy system taking into account hydrogen energy according to claim 2, wherein,
Figure FDA0004121969510000026
wherein ,cEL For the running cost of the electrolytic cell per unit time, N EL For the number of cells to be used,
Figure FDA0004121969510000027
active power consumed by the ith electrolytic cell at the time t; lambda (lambda) H2G Is the unit time operation cost of the methane synthesis equipment.
4. The method for optimizing the operation of a comprehensive energy system taking into consideration hydrogen energy according to claim 3,
Figure FDA0004121969510000031
wherein ,
Figure FDA0004121969510000032
is the carbon emission coefficient of the thermal generator set, +.>
Figure FDA0004121969510000033
For the carbon emission coefficient of the gas turbine, < > for>
Figure FDA0004121969510000034
For the carbon dioxide utilization of the integrated energy system in the methane synthesis process, < > for>
Figure FDA0004121969510000035
Is the carbon emission coefficient of the traditional fuel automobile, N EBUS N is the total number of electric vehicles in the comprehensive energy system HBUS The total number of hydrogen energy automobiles in the comprehensive energy system is calculated;
Figure FDA0004121969510000036
For the power of the ith electric car at time t, < >>
Figure FDA0004121969510000037
The power at t moment of the ith hydrogen energy automobile.
5. The method for optimizing the operation of a comprehensive energy system taking into account hydrogen energy according to claim 4,
Figure FDA0004121969510000038
Figure FDA0004121969510000039
Figure FDA00041219695100000310
wherein ,
Figure FDA00041219695100000311
for electric transmission loss, k of electric storage equipment es Loss coefficient, mu, of electric transmission ch 、μ dch P is the charge and discharge efficiency of the electricity storage equipment t ch 、P t dch The charge and discharge power of the electricity storage equipment at the moment t;
Figure FDA00041219695100000312
for gas transmission loss, k of gas storage apparatus gs Loss factor for gas transmission, +.>
Figure FDA00041219695100000313
Air inlet and air discharge efficiency of the hydrogen storage device, < >>
Figure FDA00041219695100000314
The air inlet and air discharge rates of the air storage equipment at the moment t are respectively;
Figure FDA00041219695100000315
for gas transmission loss, k of hydrogen storage device hs Loss factor for hydrogen transport, +.>
Figure FDA00041219695100000316
For the air intake and air discharge efficiency of the hydrogen storage device, < >>
Figure FDA00041219695100000317
The air inlet and air discharge rates of the hydrogen storage equipment at the moment t are respectively.
6. The method for optimizing the operation of a comprehensive energy system taking into account hydrogen energy according to claim 5, wherein,
Figure FDA00041219695100000318
wherein ,f5 To sell the benefits of hydrogen, c SH The price is sold for the unit power of hydrogen,
Figure FDA00041219695100000319
the power at t moment of the ith hydrogen energy automobile.
7. The integrated energy system operation optimization method considering hydrogen energy of claim 6, wherein the plant operation constraints include: thermal power generating unit, gas turbine, fuel cell, electrolysis trough, storage of electricity equipment, storage of hydrogen equipment and the operation constraint of storage of gas equipment.
8. The method for optimizing operation of a comprehensive energy system taking into account hydrogen energy according to claim 7, wherein the coupling constraints between the subsystems include coupling constraints between a natural gas system and a hydrogen energy storage system, coupling constraints between a natural gas system and an electric power system, and coupling constraints between a hydrogen energy storage system and an electric power system;
the coupling constraint between the natural gas system and the hydrogen energy storage system is as follows:
Figure FDA0004121969510000041
the coupling constraint between the natural gas system and the power system is as follows:
Figure FDA0004121969510000042
The coupling constraint between the hydrogen energy storage system and the power system is as follows:
Figure FDA0004121969510000043
in the formula ,
Figure FDA0004121969510000044
the gas injection amount θ at time t of the ith methane synthesizer gas Is natural gas with high heat value and mu H2G Efficiency for synthesizing methane from hydrogen, P t H2G Methane power value P of natural gas system is injected into methane synthesis device at time t t 1H The hydrogen energy power flowing into the methane synthesizing device at the moment t for the hydrogen energy storage system;
Figure FDA0004121969510000045
for the i-th gas turbine active power output at time t, mu GT For the power generation efficiency of a gas turbine, +.>
Figure FDA0004121969510000046
The natural gas consumption of the ith gas turbine at the time t is calculated;
Figure FDA0004121969510000047
is the active force of the ith hydrogen-oxygen fuel cell at the moment t, mu H2E P is the power generation efficiency of the fuel cell t 2H And inputting the hydrogen energy power of the fuel cell into the hydrogen energy storage system t at the moment. />
CN202110255688.8A 2021-03-09 2021-03-09 Comprehensive energy system operation optimization method considering hydrogen energy Active CN112989594B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110255688.8A CN112989594B (en) 2021-03-09 2021-03-09 Comprehensive energy system operation optimization method considering hydrogen energy

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110255688.8A CN112989594B (en) 2021-03-09 2021-03-09 Comprehensive energy system operation optimization method considering hydrogen energy

Publications (2)

Publication Number Publication Date
CN112989594A CN112989594A (en) 2021-06-18
CN112989594B true CN112989594B (en) 2023-05-19

Family

ID=76336271

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110255688.8A Active CN112989594B (en) 2021-03-09 2021-03-09 Comprehensive energy system operation optimization method considering hydrogen energy

Country Status (1)

Country Link
CN (1) CN112989594B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113410855A (en) * 2021-08-20 2021-09-17 杭州鸿晟电力设计咨询有限公司 Hydrogen energy storage and hydrogen supply chain equipment configuration method based on comprehensive energy system
CN114142460B (en) * 2021-11-17 2024-03-15 浙江华云电力工程设计咨询有限公司 Energy storage double-layer target optimal configuration method and terminal in comprehensive energy system
CN115079564B (en) * 2022-07-21 2022-11-04 清华四川能源互联网研究院 Decarburization path planning optimization method for regional hydrogen generation system
CN115713438A (en) * 2022-11-25 2023-02-24 国网江苏省电力有限公司常州供电分公司 Optimized scheduling method of comprehensive energy system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110009152A (en) * 2019-04-03 2019-07-12 东南大学 A kind of consideration electricity turns gas and probabilistic regional complex energy system operation robust Optimal methods
CN111738502A (en) * 2020-06-15 2020-10-02 上海交通大学 Multi-energy complementary system demand response operation optimization method for promoting surplus wind power consumption
CN112290533A (en) * 2020-09-17 2021-01-29 国家电网公司西南分部 Method for scheduling comprehensive energy micro-grid for hydrogen energy-natural gas mixed energy storage

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110009152A (en) * 2019-04-03 2019-07-12 东南大学 A kind of consideration electricity turns gas and probabilistic regional complex energy system operation robust Optimal methods
CN111738502A (en) * 2020-06-15 2020-10-02 上海交通大学 Multi-energy complementary system demand response operation optimization method for promoting surplus wind power consumption
CN112290533A (en) * 2020-09-17 2021-01-29 国家电网公司西南分部 Method for scheduling comprehensive energy micro-grid for hydrogen energy-natural gas mixed energy storage

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"大规模综合能源系统电-气-热多能潮流建模与计算方法";郭祚刚 等;《电力系统及其自动化学报》;20191031;第31卷(第10期);96-102 *
"考虑风电出力随机特性的氢储能综合能源系统优化调度";黄涛 等;《电器与能效管理技术》;20200531(第5期);72-77 *
"面向多能互补的分布式光伏与气电混合容量规划方法";马溪原 等;《电力系统自动化》;20180225;第42卷(第4期);55-63 *
考虑氢能-天然气混合储能的电-气综合能源微网日前经济调度优化;刘继春等;《电网技术》;20180131;第42卷(第1期);170-178 *

Also Published As

Publication number Publication date
CN112989594A (en) 2021-06-18

Similar Documents

Publication Publication Date Title
CN112989594B (en) Comprehensive energy system operation optimization method considering hydrogen energy
CN107482688B (en) Scheduling optimization method for carbon capture virtual power plant
CN110009152A (en) A kind of consideration electricity turns gas and probabilistic regional complex energy system operation robust Optimal methods
CN112290533A (en) Method for scheduling comprehensive energy micro-grid for hydrogen energy-natural gas mixed energy storage
CN113078687B (en) Energy optimization scheduling method for island multi-energy complementary electricity-gas coupling system
CN116542370A (en) Park low-carbon economic operation method considering carbon capture and carbon transaction
CN116260166B (en) Electric-gas coupled rural optical biogas elastic comprehensive energy system and operation model thereof
CN117236587A (en) Electric hydrogen-carbon cooperative scheduling optimization method for park comprehensive energy system
CN111126675A (en) Multi-energy complementary microgrid system optimization method
CN116247719A (en) Micro-grid two-stage robust optimal configuration method based on ladder carbon transaction
CN117254502A (en) Multi-objective optimization scheduling method based on electric-hydrogen hybrid energy storage comprehensive energy system
CN114742276A (en) Multi-objective optimization scheduling method of park comprehensive energy system with ORC (organic Rankine cycle) considering exergy efficiency
CN117973886B (en) Comprehensive energy system collaborative planning operation method and system for hydrogen-containing energy full link
CN118504880A (en) Comprehensive energy system low-carbon economic dispatching method based on CCS-P2G green certificate-carbon transaction combined interaction mechanism
CN114836776A (en) New energy coupling coal chemical industry multi-energy system, evaluation method and computer readable storage medium
CN116960970A (en) Scene-driven robust optimization regulation and control method for comprehensive energy system containing electricity and gas
CN116976590A (en) Demand response virtual power plant optimization scheduling method based on photo-thermal power station in market
CN115860178A (en) Regional comprehensive energy operation strategy considering waste treatment and methanol synthesis
CN114463130B (en) Energy system scheduling method based on ladder-type carbon transaction mechanism and demand response
CN115907363A (en) Source-load multi-time scale optimization scheduling method based on comprehensive energy system
Yang et al. Optimization of the Joint Operation of an Electricity–Heat–Hydrogen–Gas Multi-Energy System Containing Hybrid Energy Storage and Power-to-Gas–Combined Heat and Power
CN115640894B (en) Comprehensive energy system optimal scheduling method based on double time scales
CN118761570A (en) Low-carbon economic dispatching method and system for efficient utilization of hydrogen and oxygen enrichment mixture
CN117574684B (en) Sequential production simulation method and system for electric-hydrogen-carbon comprehensive energy system
Han et al. Economic Performance of a Combined SOFC System with Green Hydrogen Synthesis of Methanol and Methane

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant