CN115079564B - Decarburization path planning optimization method for regional hydrogen generation system - Google Patents
Decarburization path planning optimization method for regional hydrogen generation system Download PDFInfo
- Publication number
- CN115079564B CN115079564B CN202210856118.9A CN202210856118A CN115079564B CN 115079564 B CN115079564 B CN 115079564B CN 202210856118 A CN202210856118 A CN 202210856118A CN 115079564 B CN115079564 B CN 115079564B
- Authority
- CN
- China
- Prior art keywords
- hydrogen
- module
- power
- time section
- constraint
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/0205—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system
- G05B13/024—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system in which a parameter or coefficient is automatically adjusted to optimise the performance
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a method for planning and optimizing a decarburization path of a regional hydrogen generation system, and relates to the technical field of planning and optimizing an energy system. The method comprises the following steps: s1: establishing a framework of a regional electric hydrogen system; s2: selecting a planning period and a time section; s3: establishing variables to be optimized of a regional hydrogen generation system; s4: establishing an optimization objective function; s5: establishing a constraint condition equation; s6: configuring an optimization algorithm; s7: updating the input parameters of the regional hydrogen generating system of the current time section; s8: solving an optimization problem; s9: outputting an optimization result of the regional hydrogen generation system of the current time section, and feeding back each equipment installation in the optimization result to the next time section optimization as an existing equipment installation; s10: and judging whether the current time section is the selected last time section or not, if not, returning to S7, and if so, ending. The method can obtain the cost-optimal capacity configuration and operation parameter development path in the decarburization process of the regional hydrogen generation system.
Description
Technical Field
The invention relates to the technical field of energy system planning optimization, in particular to a method for optimizing decarburization path planning of a regional hydrogen generation system.
Background
Decarburization in the power industry is the most important link in achieving the dual carbon goal. Because new energy such as wind, light have the characteristics that volatility is strong, difficult prediction and regulation, consequently novel electric power system will need to dispose large-scale energy storage and the adjustable load of magnanimity and provide services such as peak shaver, frequency modulation to guarantee electric power system safety, satisfy the electric power load demand.
Hydrogen energy is a clean and efficient secondary energy, can be produced in a large scale by electrolyzing water, can be stored, and can be used for generating electricity by equipment such as a fuel cell, a hydrogen gas turbine and the like. In addition to the electric industry, steel, cement, petrochemical and chemical industries are the fine industries with the highest carbon emission ratio, and are also important for decarburization in these industries. The three industries are difficult to realize large-scale decarburization directly through electrification, while hydrogen green is applied to the steel industry to replace coal/natural gas to provide reducing agent and high-temperature industrial heat, the cement industry to replace coal/natural gas to provide high-temperature industrial heat, and the petrochemical and chemical industries to replace coal/natural gas to prepare hydrogen, and is an important even necessary path for decarburization in the industries. Therefore, the hydrogen energy can support the large-scale development of renewable energy sources at the source end, can realize the deep decarburization of industries difficult to decarbonize, such as steel, cement, petrochemical industry and the like, at the load end, and can realize the effect of 1+1 >.
Both power systems and hydrogen energy systems are very complex systems, and the interaction between the two systems further increases the complexity of the overall system. For the electric power system, a large amount of planning, research and design work can be developed by relevant units every year, and clean, low-carbon, safe and efficient operation and development of the electric power system are powerfully guaranteed. Due to the improvement of the complexity of the system, the electric hydrogen system is planned as a whole, and the electric hydrogen system is more important for clean, low-carbon, safe and efficient operation of the future electric power-hydrogen energy system.
Therefore, an optimization method for planning a decarburization path of a regional hydrogen generation system is urgently needed to obtain a cost optimization path for the collaborative decarburization of the regional hydrogen generation system, and meet the requirements of planning a carbon peak reaching path and a carbon neutralization path.
Disclosure of Invention
The invention aims to provide a regional electric hydrogen system decarburization path planning optimization method which can obtain the most cost-effective capacity configuration and operation parameter development path in the regional electric hydrogen system decarburization process and meet the requirements of carbon peak reaching, carbon neutralization and path planning.
Embodiments of the invention may be implemented as follows:
the invention provides a regional hydrogen generation system decarburization path planning and optimizing method, which comprises the following steps:
s1: establishing a framework of a regional electric hydrogen system;
s2: selecting a planning period and a time section;
s3: establishing variables to be optimized of a regional hydrogen generation system;
s4: establishing an optimization objective function;
s5: establishing a constraint condition equation;
s6: configuring an optimization algorithm;
s7: updating the input parameters of the regional hydrogen generating system of the current time section;
s8: solving an optimization problem;
s9: outputting an optimization result of the regional hydrogen generation system of the current time section, and feeding back each equipment installation in the optimization result to the next time section optimization as an existing equipment installation;
s10: and judging whether the current time section is the selected last time section or not, if not, returning to S7, and if so, ending the process.
In an alternative embodiment, in S1, the regional electric hydrogen system includes a power supply unit, a power input unit, a power balancing unit, a power load unit, an electricity storage unit, an electric hydrogen production unit, a hydrogen power generation unit, a hydrogen input unit, a hydrogen storage unit, a hydrogen balancing unit, a hydrogen load unit, a power output unit, and a hydrogen output unit;
the power supply unit, the power input unit, the power output unit, the power load unit, the electricity storage unit, the hydrogen production unit and the hydrogen generation unit are all connected with the power balance unit, and the hydrogen production unit, the hydrogen input unit, the hydrogen output unit, the hydrogen generation unit, the hydrogen storage unit and the hydrogen load unit are all connected with the hydrogen balance unit.
In an optional embodiment, in S1, there are one or more power balancing units, the power balancing units are connected to each other, and there is bidirectional or unidirectional power transmission between the power balancing units connected to each other; the hydrogen balance units are connected with each other, and bidirectional or unidirectional hydrogen transmission is carried out between the hydrogen balance units connected with each other.
In an alternative embodiment, in S2, the planning period spans a plurality of consecutive natural years, the single time section is one natural year in the planning period, and the time sections are one or more.
In an alternative embodiment, in S3, the variables to be optimized of the regional electro-hydrogen system include real variables and one-dimensional array variables.
In an alternative embodiment, the real-type variables include installed capacity of each power module, installed capacity of each power input module, installed capacity of each power output module, installed capacity of each power storage module, installed capacity of each hydrogen production module, installed capacity of each hydrogen generation module, installed capacity of each hydrogen input module, installed hydrogen amount of each hydrogen load module, installed hydrogen amount of each hydrogen output module, and installed capacity of each hydrogen storage module for a selected time cross-section in the regional electric hydrogen system;
each one-dimensional array variable comprises 8760 real-type variables, and the one-dimensional array variables comprise the hourly power generation power of each power module, the hourly power rejection power of each power module, the hourly charge and discharge power of each power storage module, the hourly charge/energy storage state of each power storage module, the hourly running power of each hydrogen production module, the hourly hydrogen production rate of each hydrogen production module, the hourly power generation power of each hydrogen generation module, the hourly hydrogen consumption rate of each hydrogen generation module, the hourly charge and discharge hydrogen rate of each hydrogen storage module and the hourly hydrogen storage state of each hydrogen storage module of a selected time section in the regional electric hydrogen system.
In an alternative embodiment, in S4, the objective function is optimized to minimize the total system cost in the time section, including the allocation of fixed investment, fixed operation and maintenance costs, variable operation costs, and carbon emission costs of each module in the system, which are expressed as:
in the formula:the modules indicating the running cost without change in the system comprise a power module, an electricity storage module, an electricity hydrogen production module, a hydrogen storage module, a hydrogen generation module and the like which do not need to consume fuel;the modules indicating the variable running cost in the system comprise a coal power module, a gas power module, a nuclear power module and other power modules which need to consume fuel;the system refers to input and output modules of electric power and hydrogen energy;is a moduleNewly adding an installation machine from the last time section to the current time section;the module is put into operation from the last time section to the current time sectionAverage unit capital investment of (1);is put into operation from the last time section to the current time sectionAverage design life of (c);is a moduleInstalling the existing machine on the section at the last time;is a moduleAn installation to be decommissioned from a previous time section to the current time section;is a moduleThe ratio of year fixed operation and maintenance cost to fixed investment;is a moduleIn the first placeGenerated power within an hour;is a moduleAverage electrical fuel cost from last time section to present time section;is the firstInput or output of electrical/hydrogen energy in hours;is the average unit price of input or output electric power/hydrogen energy from the last time section to the present time section;the annual carbon emission of the system;is the average carbon number from the last time slice to this time slice.
In an alternative embodiment, in S5, the constraint equations include a system balance constraint equation, a system module constraint equation, and a system design constraint equation.
In an alternative embodiment, in S5, the system balance constraint equations include a power balance node equation and a hydrogen balance node equation;
the power balance node equation is as follows:
in the formula:is the other power balance node connected with the power balance nodeIn the first placePower exchange between an hour and the power balancing node;is a power supply module connected with the power balance nodeIn the first placePower in hours;is a power input module connected with the power balance nodeIn the first placeHourly delivered power;is a power load module connected with the power balance nodeIn the first placePower in hours;is a power output module connected with the power balance nodeIn the first placeHourly delivered power;is an electricity storage module connected with the power balance nodeIn the first placePower in hours;is an electric hydrogen production module connected with the electric power balance nodeIn the first placePower in hours;is a hydrogen power generation module connected with the power balance nodeIn the first placeHourly power. The power flowing into the power balance node is positive, and the power flowing out of the power balance node is negative;
the hydrogen equilibrium node equation is:
in the formula:is the other hydrogen balance node connected to the hydrogen balance nodeIn the first placeHydrogen mass exchange between hours and the hydrogen balance node;is an electric hydrogen production module connected with the hydrogen balance nodeIn the first placeHydrogen production in hours;is a hydrogen power generation module connected with the hydrogen balance nodeIn the first placeHydrogen consumption in hours;is a hydrogen input module connected with the hydrogen balance nodeIn the first placeHydrogen input per hour;is a hydrogen load module connected with the hydrogen balance nodeIn the first placeHydrogen usage in hours;is a hydrogen output module connected with the hydrogen balance nodeIn the first placeHydrogen output per hour;is andthe hydrogen storage module connected with the hydrogen balance nodeIn the first placeHydrogen charge/discharge amount per hour. The amount of hydrogen flowing into the hydrogen balance node is positive and the amount of hydrogen flowing out of the hydrogen balance node is negative.
In an alternative embodiment, in S5, the system module constraint equations include a power module constraint equation, an electricity storage module constraint equation, an electricity hydrogen generation module constraint equation, a hydrogen load module constraint equation, and a hydrogen storage module constraint equation; the power module constraint equation comprises power generation power constraint and renewable energy power abandon rate constraint;
the generated power constraint is as follows:
in the formula:non-renewable power moduleIn the first placeHourly generated power;non-renewable power moduleInstalled capacity of (d);
in the formula:is a renewable power moduleIn the first placeHourly generated power;is a renewable power moduleIn the first placeElectric power discard of hours;is a renewable power moduleInstalled capacity of (a);is a renewable power moduleIn the first placeThe ratio of maximum hourly output to installed capacity;
the renewable energy power abandon rate constraint is as follows:
the electric storage module constraint equation comprises charge/energy storage state constraint, charge and discharge power constraint and capacity balance constraint;
the charge/storage state constraints are:
in the formula:andare respectively an electricity storage moduleLower and upper limits of the state of charge/energy storage;is an electricity storage moduleIn the first placeHourly charge/energy storage state;
the charge and discharge power constraint is as follows:
in the formula:is an electricity storage moduleThe capacity of (a);is an electricity storage moduleMaximum charge rate of;is an electricity storage moduleMaximum discharge rate of (d);
the capacity balance constraint is:
In the formula:is an electricity storage moduleEnergy conversion efficiency of charging;is an electricity storage moduleEnergy conversion efficiency of the discharge;is an electricity storage moduleA ratio of leakage to capacity per hour;is an electricity storage moduleInitial state of charge/energy storage;
the constraint equation of the electrohydrogen production module comprises operation power constraint and energy conversion constraint;
the operating power constraints are:
the energy conversion constraint is:
in the formula:Is an electric hydrogen production moduleElectricity consumption per unit hydrogen production;
the hydrogen power generation module constraint equation comprises an operating power constraint and an energy conversion constraint;
the operating power constraints are:
the energy conversion constraint is:
in the formula:is a hydrogen power generation moduleHydrogen consumption per unit of generated energy;
the hydrogen load module constraint equation comprises a hydrogen usage constraint;
the hydrogen dosage constraints are:
in the formula:is a hydrogen load moduleThe maximum possible amount of hydrogen used per year in the time section;is a hydrogen load moduleAnnual production/traffic in time section;is a hydrogen load moduleAdopting the carbon emission factor of the prior art;is a hydrogen load moduleThe proportion of emission reduction can be realized by replacing hydrogen in the carbon emission in the prior art;is a hydrogen load moduleCarbon reduction per unit hydrogen application;
the hydrogen storage module constraint equation comprises hydrogen storage state constraint, hydrogen charging and discharging rate constraint and capacity balance constraint;
the hydrogen storage state constraints are as follows:
in the formula:andare each a hydrogen storage moduleLower and upper limits of the hydrogen storage state of (a);is a hydrogen storage moduleIn the first placeHydrogen storage status of hours;
the hydrogen charging and discharging rate constraint is as follows:
in the formula:is a hydrogen storage moduleThe capacity of (a);is a hydrogen storage moduleMaximum hydrogen charge rate of;is a hydrogen storage moduleThe maximum hydrogen release rate;
the capacity balance constraints are:
In the formula:is a hydrogen storage moduleThe efficiency of the charging process;is a hydrogen storage moduleThe efficiency of the hydrogen discharge process;is a hydrogen storage moduleThe ratio of hydrogen leakage to capacity per hour;is a hydrogen storage moduleInitial hydrogen storage state.
In an optional embodiment, in S5, the system design constraint equation includes system carbon emission constraints, upper and lower limit constraints of installed capacity of the power supply module, upper and lower limit constraints of annual input power of the power input module, upper and lower limit constraints of annual output power of the power output module, upper and lower limit constraints of capacity of the power storage module, upper and lower limit constraints of installed capacity of the hydrogen production module, upper and lower limit constraints of installed capacity of the hydrogen generation module, upper and lower limit constraints of annual input hydrogen amount of the hydrogen input module, upper and lower limit constraints of annual hydrogen amount of the hydrogen load module, upper and lower limit constraints of annual output hydrogen amount of the hydrogen output module, and upper and lower limit constraints of capacity of the hydrogen storage module;
wherein the system carbon emission constraint is:
in the formula:is a power supply moduleCarbon emission factor of (a);is a power input moduleCarbon emission factor of (a);is a power output moduleThe carbon emission factor of (c);is a hydrogen input moduleThe carbon emission factor of (c);is a hydrogen output moduleCarbon emission factor of (a);is the upper carbon emission limit of the regional electro-hydrogen system over the time section.
In an optional implementation manner, in S6, the configuration optimization algorithm selects a linear solver in the MATLAB according to the optimization objective function and the constraint condition equation.
In an alternative embodiment, in S7, the input parameters of the regional electric hydrogen system of the current time section include an initial existing installed capacity of each power module, upper and lower limits of the installed capacity, a carbon emission factor, a unit fixed investment, an average design life, a ratio of an annual fixed operation and maintenance cost to a fixed investment, a fuel cost of the electric fuel, and a ratio of a time-by-time maximum output of renewable energy to the installed capacity;
the input parameters of the regional electro-hydrogen system of the current time section further comprise the ratio of the time-by-time transmission power to the peak transmission power of the power input module;
the input parameters of the regional hydrogen power system of the current time section further comprise the time-by-time power of the power load module;
the input parameters of the regional electric hydrogen system of the current time section further comprise the ratio of the time-by-time transmission power to the peak transmission power of the power output module;
the input parameters of the regional hydrogen-electricity system of the current time section further comprise the initial existing installed capacity of each electricity storage module, the upper and lower limits of the installed capacity, the maximum charging multiplying factor, the maximum discharging multiplying factor, the energy conversion efficiency in the charging process, the energy conversion efficiency in the discharging process, the ratio of the leakage quantity per hour to the capacity, the unit fixed investment, the average design life and the ratio of the annual fixed operation and maintenance cost to the fixed investment;
the input parameters of the regional electric hydrogen system of the current time section further comprise the initial existing installed capacity of the electric hydrogen production module, the upper and lower limits of the installed capacity, the power consumption of unit hydrogen production, unit fixed investment, the average design life and the ratio of annual fixed operation and maintenance cost to fixed investment;
the input parameters of the regional power-hydrogen system of the current time section further comprise the initial existing installed capacity of each hydrogen power generation module, the upper limit and the lower limit of the installed capacity, the hydrogen consumption of unit power generation, unit fixed investment, the average design life and the ratio of annual fixed operation and maintenance cost to fixed investment;
the input parameters of the regional electric hydrogen system of the current time section further comprise the ratio of the hourly hydrogen output to the peak hydrogen output of the hydrogen input module;
the input parameters of the regional electric hydrogen system of the current time section further comprise the ratio of the hourly hydrogen consumption to the peak hydrogen consumption of each hydrogen load module, the annual output/transport capacity, a carbon emission factor adopting the prior art, the ratio of realizing emission reduction by hydrogen substitution in the carbon emission of the prior art, the carbon reduction amount brought by each unit of hydrogen application and the ratio of the hourly hydrogen consumption to the peak hydrogen consumption of the hydrogen output module;
the input parameters of the regional hydrogen power system of the current time section further comprise the ratio of the hourly hydrogen output quantity to the peak hydrogen output quantity of each hydrogen output module;
the input parameters of the regional hydrogen power system of the current time section further comprise initial existing installed capacity of each hydrogen storage module, upper and lower limits of the installed capacity, maximum hydrogen charging rate, maximum hydrogen discharging rate, efficiency of a hydrogen charging process, efficiency of a hydrogen discharging process, a ratio of hydrogen leakage amount per hour to capacity, unit fixed investment, average design life and a ratio of annual fixed operation and maintenance cost to fixed investment.
The method for optimizing the decarburization path planning of the regional hydrogen generation system provided by the embodiment of the invention has the beneficial effects that:
1. the method and the tool are provided for planning the decarburization path of the regional hydrogen generation system;
2. the method can realize the optimization of the capacity allocation of each component element with the optimal total system cost in the process of advancing carbon neutralization of the electro-hydrogen system in the region;
3. the method can obtain the hourly operation parameters of each component element of the hydrogen power system in the region in the selected planning period and the time section, and can guide the selection of the future operation mode and the operation scheduling arrangement of the hydrogen power system;
4. the method has higher flexibility, and can select and determine the types and the number of elements and modules considered in the system according to the characteristics of the planning area and the requirement on the spatial resolution.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and for those skilled in the art, other related drawings can be obtained according to the drawings without inventive efforts.
Fig. 1 is a flowchart of a method for optimizing a decarburization path planning of a regional hydrogen generation system according to an embodiment of the present invention;
fig. 2 is a schematic diagram of a regional electro-hydrogen system according to an embodiment of the present invention.
Icon: 100-regional hydrogen generation system; 1-a power supply unit; 11-a coal electric module; 12-a gas electric module; 13-a nuclear power module; 14-a hydro-electric module; 15-offshore wind power module; 16-an onshore wind power module; 17-a photovoltaic module; 2-a power input unit; 3-a power balancing unit; 4-an electrical load unit; 5-an electricity storage unit; 51-a battery module; 52-a water pumping energy storage module; 6-an electrohydrogen production unit; 7-a hydrogen power generation unit; 71-a fuel cell module; 72-a hydrogen turbine module; an 8-hydrogen input unit; 9-a hydrogen storage unit; 91-salt cavern hydrogen storage module; 92-a high pressure tank module; a 10-hydrogen equilibrium unit; an 18-hydrogen load cell; 181-steel module; 182-a cement module; 183-chemical module; 184-road freight module; 185-shipping module; 186-an aviation module; 19-a power output unit; 20-hydrogen output unit.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. The components of embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.
Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be obtained by a person skilled in the art without inventive step based on the embodiments of the present invention, are within the scope of protection of the present invention.
It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.
In the description of the present invention, it should be noted that if the terms "upper", "lower", "inside", "outside", etc. indicate an orientation or a positional relationship based on that shown in the drawings or that the product of the present invention is used as it is, this is only for convenience of description and simplification of the description, and it does not indicate or imply that the device or the element referred to must have a specific orientation, be constructed in a specific orientation, and be operated, and thus should not be construed as limiting the present invention.
Furthermore, the appearances of the terms "first," "second," and the like, if any, are used solely to distinguish one from another and are not to be construed as indicating or implying relative importance.
It should be noted that the features of the embodiments of the present invention may be combined with each other without conflict.
Referring to fig. 1, the present embodiment provides a method for optimizing a decarburization path planning of a local hydrogen generation system 100, which includes the following steps:
s1: the architecture of the regional power hydrogen system 100 is established.
Referring to fig. 2, the regional hydrogen generation system 100 includes a power supply unit 1, a power input unit 2, a power balance unit 3, a power load unit 4, a power storage unit 5, an electrical hydrogen production unit 6, a hydrogen power generation unit 7, a hydrogen input unit 8, a hydrogen storage unit 9, a hydrogen balance unit 10, a hydrogen load unit 18, a power output unit 19, and a hydrogen output unit 20.
The number of the power balance units 3 is one or more, the power balance units 3 are connected with each other, and bidirectional or unidirectional power transmission is performed between the power balance units 3 connected with each other; the number of the hydrogen balance units 10 is one or more, the hydrogen balance units 10 are connected with each other, and bidirectional or unidirectional hydrogen transmission exists between the hydrogen balance units 10 connected with each other.
The power supply unit 1 may include a coal electric module 11, a gas electric module 12, a nuclear electric module 13, a hydroelectric module 14, an offshore wind electric module 15, an onshore wind electric module 16, a photovoltaic module 17, and other power supply module types considered within the planning area. Each power module type may comprise one or more modules, each power module being connected to one power balancing unit 3.
The power input unit 2 may comprise one or more power input modules, each connected to one power balancing unit 3;
the power load unit 4 may comprise one or more power load modules, each connected to one of the power balancing units 3.
The power output unit 19 may comprise one or more power output modules, each connected to one of the power balancing units 3.
The electrical storage units 5 may include battery modules 51, pumped-water energy storage modules 52, flywheel modules, and other electrical storage unit types considered within the planning area. Each power storage unit may contain one or more modules, each connected to the power balancing unit 3.
The electrical hydrogen production units 6 include alkaline water electrolysis hydrogen production, proton exchange membrane water electrolysis hydrogen production and other electrical hydrogen production module types considered in the planned region, each electrical hydrogen production unit 6 may include one or more modules, and each electrical hydrogen production module is respectively connected with one power balance unit 3 and one hydrogen balance unit 10.
The hydrogen power generation unit 7 may include a fuel cell module 71, a hydrogen turbine module 72, and other hydrogen power generation module types considered within the planning area. Each hydrogen generating unit 7 may comprise one or more modules, each of which is connected to the power balancing unit 3 and the hydrogen balancing unit 10, respectively.
The hydrogen input unit 8 may comprise one or more hydrogen input modules, each of which is connected to a hydrogen balancing unit 10.
The hydrogen output unit 20 may comprise one or more hydrogen output modules, each of which is connected to one of the hydrogen balancing units 10.
The hydrogen load cells 18 may include steel modules 181, cement modules 182, chemical modules 183, road freight modules 184, shipping modules 185, aviation modules 186, and other hydrogen load module types contemplated within the planning area. Each hydrogen load unit 18 may comprise one or more modules, each hydrogen load module being connected to the hydrogen balancing unit 10.
The hydrogen storage unit 9 may include a salt cavern hydrogen storage module 91, a high pressure storage tank module 92, and other hydrogen storage module types contemplated within the planned region. Each hydrogen storage unit 9 may comprise one or more modules, each connected to a hydrogen balancing unit 10.
S2: and selecting a planning period and a time section.
The planning period may span 2020-2060 years, with a single time slice being a natural year within the planning period, with time slices selected to be 2025, 2030, 2035, 2040, 2050, 2060 years.
S3: variables to be optimized for the regional electro-hydrogen system 100 are established.
The variables to be optimized of the regional electric hydrogen system 100 include two types of real variables and one-dimensional array variables.
The real-type variables include the installed capacity of each power module, the annual input power of each power input module, the annual output power of each power output module, the capacity of each power storage module, the installed capacity of each hydrogen production module, the installed capacity of each hydrogen generation module, the annual input hydrogen amount of each hydrogen input module, the annual hydrogen usage amount of each hydrogen load module, the annual output hydrogen amount of each hydrogen output module, and the capacity of each hydrogen storage module for a selected time section in the regional electric hydrogen system 100.
Each one-dimensional array variable consists of 8760 real type variables. The one-dimensional array variables include the hourly power generation power of each power module of the selected time section in the regional electrical hydrogen system 100, the hourly discarded electric power of each power module (this variable exists only for renewable energy), the hourly charge-discharge power of each power storage module, the hourly charge/energy storage state (ratio of stored energy to capacity) of each power storage module, the hourly operating power of each electrical hydrogen production module, the hourly hydrogen production rate of each electrical hydrogen production module, the hourly power generation power of each hydrogen generation module, the hourly hydrogen consumption rate of each hydrogen storage module, and the hourly hydrogen storage state (ratio of stored hydrogen to charged-discharged capacity) of each hydrogen storage module.
S4: and establishing an optimization objective function.
The objective function is optimized to minimize the total system cost in the time section, which includes the allocation of fixed investment, fixed operation and maintenance cost, variable operation cost, and carbon emission cost of each module in the system, and can be expressed as:
in the formula:the modules indicating the running cost without change in the system comprise a power module, an electricity storage module, an electricity hydrogen production module, a hydrogen storage module, a hydrogen generation module and the like which do not need to consume fuel;the modules refer to the modules with variable running cost in the system, and comprise a coal power module 11, a gas power module 12, a nuclear power module 13 and other power modules which need to consume fuel;an input module for power and hydrogen energy;is a moduleA newly-added installation machine is needed from the last time section to the current time section;is a moduleA newly-added installation machine is needed from the last time section to the current time section;the module is put into operation from the last time section to the current time sectionAverage unit capital investment of (1);is the module commissioned from the last time section to the current time sectionAverage design life of;is the module commissioned from the last time section to the current time sectionAverage design life of (c);is a moduleInstalling the existing machine on the section at the last time;is a moduleAn installation to be decommissioned from a previous time section to the current time section;is a moduleThe ratio of year fixed operation and maintenance cost to fixed investment;is a moduleIn the first placeGenerated power within an hour;is a moduleAverage electrical fuel cost during the last time section to the present time section;is the firstInput or output of electrical/hydrogen energy in hours;is the average unit price of input or output electric power/hydrogen energy from the last time section to the present time section;the annual carbon emission of the system;is the average carbon number from the last time slice to this time slice.
S5: and establishing a constraint condition equation.
The constraint condition equations comprise a system balance constraint equation, a system module constraint equation and a system design constraint equation. The system balance constraint equation comprises a power balance node equation and a hydrogen balance node equation.
The power balance node equation is as follows:
in the formula:is the other power balance node connected with the power balance nodeIn the first placePower exchange between an hour and the power balancing node;is a power supply module connected with the power balance nodeIn the first placePower in hours;is a power input module connected with the power balance nodeIn the first placeHourly delivered power;is a power load module connected with the power balance nodeIn the first placePower in hours;is a power output module connected with the power balance nodeIn the first placeHourly delivered power;is an electricity storage module connected with the power balance nodeIn the first placePower in hours;is an electric hydrogen production module connected with the electric power balance nodeIn the first placePower in hours;is a hydrogen power generation module connected with the power balance nodeIn the first placeHourly power. The power flowing into the power balancing node is positive and the power flowing out of the power balancing node is negative.
The hydrogen balance node equation is:
in the formula:is the other hydrogen balance node connected with the hydrogen balance nodeIn the first placeHydrogen mass exchange between hours and the hydrogen balance node;is an electric hydrogen production module connected with the hydrogen balance nodeIn the first placeHydrogen production in hours;is a hydrogen power generation module connected with the hydrogen balance nodeIn the first placeHydrogen consumption in hours;is a hydrogen input module connected to the hydrogen balance nodeIn the first placeHydrogen output per hour;is a hydrogen load module connected with the hydrogen balance nodeIn the first placeHydrogen usage in hours;is a hydrogen output module connected with the hydrogen balance nodeIn the first placeHydrogen input per hour;is a hydrogen storage module connected with the hydrogen balance nodeIn the first placeHydrogen charge/discharge amount per hour. The amount of hydrogen flowing into the hydrogen balance node is positive and the amount of hydrogen flowing out of the hydrogen balance node is negative.
The system module constraint equations comprise a power module constraint equation, an electricity storage module constraint equation, an electricity hydrogen production module constraint equation, a hydrogen generation module constraint equation, a hydrogen load module constraint equation and a hydrogen storage module constraint equation.
The power module constraint equation comprises a power generation power constraint and a renewable energy curtailment rate constraint.
The generated power constraint is as follows:
in the formula:non-renewable power moduleIn the first placeHourly generated power;non-renewable power moduleInstalled capacity of (c).
In the formula:is a renewable power moduleIn the first placeHourly generated power;is a renewable power moduleIn the first placeElectric power discard per hour;is a renewable power moduleInstalled capacity of (d);is a renewable power moduleIn the first placeThe ratio of maximum hourly output to installed capacity.
The renewable energy power abandon rate constraint is as follows:
The energy storage module constraint equation comprises charge/energy storage state constraint, charge and discharge power constraint and capacity balance constraint.
The charge/storage state constraints are:
in the formula:andare respectively an electricity storage moduleLower and upper limits of the state of charge/energy storage;is an electricity storage moduleIn the first placeCharge/storage state in hours.
The charge and discharge power constraint is as follows:
in the formula:is an electricity storage moduleThe capacity of (a);is an electricity storage moduleMaximum charge rate of (d);is an electricity storage moduleMaximum discharge rate of (d);is an electricity storage moduleThe charge and discharge power of (1).
The capacity balance constraint is:
In the formula:is an electricity storage moduleEnergy conversion efficiency of charging;is an electricity storage moduleEnergy conversion efficiency of the discharge;is an electricity storage moduleA ratio of leakage to capacity per hour;is an electricity storage moduleInitial state of charge/energy storage.
The electrohydrogen production module constraint equations include an operating power constraint and an energy conversion constraint.
The operating power constraints are:
in the formula:is an electric hydrogen production moduleInstalled capacity of (d);is an electric hydrogen production moduleThe operating power of (c).
The energy conversion constraint is:
in the formula:is an electric hydrogen production moduleElectricity consumption per unit hydrogen production;is an electric hydrogen production moduleThe amount of energy conversion.
The hydrogen power generation module constraint equations include an operating power constraint and an energy conversion constraint.
The operating power constraints are:
The energy conversion constraint is:
in the formula:is a hydrogen power generation moduleHydrogen consumption per unit of generated energy;is a hydrogen power generation module。
The hydrogen load module constraint equation includes a constraint with an amount of hydrogen.
The hydrogen dosage constraints are:
in the formula:is a hydrogen load moduleThe maximum possible amount of hydrogen used per year in the time section;is a hydrogen load moduleAnnual production/traffic in time section;is a hydrogen load moduleAdopting the carbon emission factor of the prior art;is a hydrogen load moduleThe proportion of emission reduction can be realized by replacing hydrogen in the carbon emission in the prior art;is a hydrogen load moduleCarbon reduction per unit hydrogen application;。
the hydrogen storage module constraint equation comprises hydrogen storage state constraint, hydrogen charging and discharging rate constraint and capacity balance constraint.
The hydrogen storage state constraints are as follows:
in the formula:andare each a hydrogen storage moduleThe lower limit and the upper limit of the hydrogen storage state of (a);is a hydrogen storage moduleIn the first placeHydrogen storage state in hours.
The hydrogen charging and discharging rate constraint is as follows:
in the formula:is a hydrogen storage moduleThe capacity of (a);is a hydrogen storage moduleMaximum hydrogen charge rate of;is a hydrogen storage moduleThe maximum hydrogen release rate;。
the capacity balance constraint is:
In the formula:is a hydrogen storage moduleThe efficiency of the charging process;is a hydrogen storage moduleThe efficiency of the hydrogen discharge process;is a hydrogen storage moduleThe ratio of hydrogen leakage to capacity per hour;is a hydrogen storage moduleInitial hydrogen storage state.
The system design constraint equation comprises system carbon emission constraint, upper and lower limit constraint of installed capacity of a power supply module, upper and lower limit constraint of annual input electric quantity of a power input module, upper and lower limit constraint of annual output electric quantity of a power output module, upper and lower limit constraint of capacity of an electricity storage module, upper and lower limit constraint of installed capacity of an electricity hydrogen production module, upper and lower limit constraint of installed capacity of a hydrogen power generation module, upper and lower limit constraint of annual input hydrogen quantity of a hydrogen input module, upper and lower limit constraint of annual hydrogen quantity of a hydrogen load module, upper and lower limit constraint of annual output hydrogen quantity of a hydrogen output module and upper and lower limit constraint of capacity of a hydrogen storage module.
The system carbon emission constraints are:
in the formula:is a power supply moduleThe carbon emission factor of (c);is a power input moduleThe carbon emission factor of (c);is a power output moduleThe carbon emission factor of (c);is a hydrogen input moduleCarbon emission factor of (a);is a hydrogen output moduleThe carbon emission factor of (c);is the upper carbon emission limit of the regional electro-hydrogen system over the time section.
S6: and configuring an optimization algorithm.
The configuration optimization algorithm comprises the steps of selecting a solver in commercial software or self-writing a solving algorithm according to an optimization objective function and a constraint condition equation (linear programming problem). In this embodiment, the configuration optimization algorithm selects a linear solver in the MATLAB according to the optimization objective function and the constraint equation.
S7: the input parameters of the regional electro-hydrogen system 100 for the current time section are updated.
The input parameters of the regional electric hydrogen system 100 of the current time section include the initial existing installed capacity of each power module, the upper and lower limits of the installed capacity, carbon emission factors, unit fixed investment, average design life, the ratio of annual fixed operation and maintenance cost to fixed investment, the cost of electric fuel, and the ratio of the hourly maximum output of renewable energy to the installed capacity;
the input parameters of the regional electro-hydrogen system 100 for the current time section also include the ratio of the time-to-time delivered power to the peak delivered power of the power input module;
the input parameters of the regional electro-hydrogen system 100 for the current time section also include the time-by-time power of the power load modules;
the input parameters of the regional electro-hydrogen system 100 for the current time section also include the ratio of the time-to-time delivered power to the peak delivered power of the power output module;
the input parameters of the regional hydrogen generation system 100 of the current time section further include the initial existing installed capacity of each electricity storage module, the upper and lower limits of the installed capacity, the maximum charging rate, the maximum discharging rate, the energy conversion efficiency in the charging process, the energy conversion efficiency in the discharging process, the ratio of the leakage amount per hour to the capacity, the unit fixed investment, the average design life, and the ratio of the annual fixed operation and maintenance cost to the fixed investment;
the input parameters of the regional hydrogen generation system 100 of the current time section further include the initial existing installed capacity of the hydrogen generation module, the upper and lower limits of the installed capacity, the power consumption of the unit hydrogen generation amount, the unit fixed investment, the average design life and the ratio of the annual fixed operation and maintenance cost to the fixed investment;
the input parameters of the regional hydrogen generation system 100 of the current time section further include the initial existing installed capacity of each hydrogen generation module, the upper and lower limits of the installed capacity, the hydrogen consumption of unit generated energy, unit fixed investment, average design life, and the ratio of annual fixed operation and maintenance cost to fixed investment;
the input parameters of the regional electrical hydrogen system 100 of the current time section further include the ratio of the hourly hydrogen output to the peak hydrogen output of the hydrogen input module;
the input parameters of the regional electric hydrogen system 100 of the current time section further include the ratio of the hourly hydrogen consumption to the peak hydrogen consumption of each hydrogen load module, the annual output/transport volume, the carbon emission factor of the existing process, the proportion of the existing process for realizing emission reduction by hydrogen substitution in carbon emission, the carbon reduction amount brought by each unit of hydrogen application, and the ratio of the hourly hydrogen consumption to the peak hydrogen consumption of the hydrogen output module;
the input parameters of the regional electric hydrogen system 100 of the current time section further include the ratio of the hourly hydrogen output to the peak hydrogen output of each hydrogen output module;
the input parameters of the regional hydrogen generation system 100 of the current time section further include the initial installed capacity of each hydrogen storage module, the upper and lower limits of the installed capacity, the maximum hydrogen charging rate, the maximum hydrogen discharging rate, the efficiency of the hydrogen charging process, the efficiency of the hydrogen discharging process, the ratio of the hydrogen leakage amount per hour to the capacity, the unit fixed investment, the average design life, and the ratio of the annual fixed operation and maintenance cost to the fixed investment. S8: and (5) solving an optimization problem.
S9: and outputting an optimization result of the regional hydrogen generation system 100 of the current time section, and feeding back each equipment installation in the optimization result to the next time section optimization as an existing equipment installation.
S10: and judging whether the current time section is the selected last time section or not, if not, returning to S7, and if so, ending the process.
The method for optimizing the decarburization path planning of the regional hydrogen generation system 100 provided by the embodiment has the beneficial effects that:
1. the method has the advantages that various elements such as a power supply, a load, energy storage, hydrogen production and hydrogen storage in the system are fully considered, the fixed investment cost, the operation and maintenance cost and the carbon emission cost in the whole system are considered, the constraints of the aspects such as electric balance, hydrogen balance, carbon emission and equipment operation are considered, the optimal cost capacity configuration and the optimal cost operation parameter development path in the decarburization process of the regional hydrogen-electricity system 100 are obtained, and the method has certain practical significance;
2. a set of methods and tools are provided for planning a decarbonization path for the regional electro-hydrogen system 100;
3. the method can realize the optimization of the capacity allocation of each component element with the optimal total system cost in the process of advancing carbon neutralization of the electro-hydrogen system in the region;
4. the method can obtain the hourly operation parameters of each component element of the hydrogen power system in the region in the selected planning period and the time section, and can guide the selection of the future operation mode and the operation scheduling arrangement of the hydrogen power system;
5. the method has higher flexibility, and the types and the number of elements and modules considered in the system can be selected and determined according to the characteristics of the planning area and the requirement on the spatial resolution.
The above description is only for the specific embodiments of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (8)
1. A regional electric hydrogen system decarburization path planning optimization method is characterized by comprising the following steps:
s1: establishing a framework of a regional electric hydrogen system (100);
s2: selecting a planning period and a time section;
s3: establishing variables to be optimized for a regional electro-hydrogen system (100);
s4: establishing an optimization objective function;
s5: establishing a constraint condition equation;
s6: configuring an optimization algorithm;
s7: updating input parameters of the regional electric hydrogen system (100) of the current time section;
s8: solving an optimization problem;
s9: outputting an optimization result of the regional hydrogen generation system (100) of the current time section, wherein each equipment installation in the optimization result is used as an existing equipment installation and fed back to the next time section for optimization;
s10: judging whether the current time section is the selected last time section or not, if not, returning to S7, and if so, ending the process;
in S4, the optimization objective function is to minimize the total system cost in the time section, including the allocation of fixed investment, fixed operation and maintenance cost, variable operation cost, and carbon emission cost of each module in the system, and is expressed as:
in the formula:modules indicating operating costs within the system that are not varied include a power module, an electricity storage module, an electrical hydrogen production module, a hydrogen storage module, and a hydrogen generation module that do not require fuel consumption;the modules refer to the modules with variable running cost in the system, and comprise a coal-electricity module (11), a gas-electricity module (12) and a nuclear power module (13);the input and output module refers to electric power and hydrogen energy;is a moduleA newly-added installation machine is needed from the last time section to the current time section;is a moduleNewly adding an installation machine from the last time section to the current time section;the module is put into operation from the last time section to the current time sectionAverage unit fixed investment of (2);the module is put into operation from the last time section to the current time sectionAverage unit capital investment of (1);is put into operation from the last time section to the current time sectionAverage design life of (c);is put into operation from the last time section to the current time sectionAverage design life of (c);is a moduleInstalling the existing machine on the section at the last time;is a moduleAn installation to be decommissioned from a previous time section to the current time section;is a moduleYear fixed operation and maintenance cost and fixed investment ratioIs a moduleThe ratio of year fixed operation and maintenance cost to fixed investment;is a moduleIn the first placeGenerated power in hours;is a moduleAverage electrical fuel cost during the last time section to the present time section;is the firstInput or output of electrical/hydrogen energy in hours;is the average unit price of input or output electric power/hydrogen energy from the last time section to the present time section;the annual carbon emission of the system;is the average carbon number from the last time slice to the present time slice;
in S5, the constraint condition equations comprise a system balance constraint equation, a system module constraint equation and a system design constraint equation;
in S5, the system balance constraint equation comprises a power balance node equation and a hydrogen balance node equation;
the power balance node equation is as follows:
in the formula:is the other power balance node connected with the power balance nodeIn the first placePower exchange between an hour and the power balancing node;is a power supply module connected with the power balance nodeIn the first placePower in hours;is a power input module connected with the power balance nodeIn the first placeHourly delivered power;is a power load module connected with the power balance nodeIn the first placePower in hours;is a power output module connected with the power balance nodeIn the first placeHourly delivered power;is an electricity storage module connected with the power balance nodeIn the first placePower in hours;is an electric hydrogen production module connected with the electric power balance nodeIn the first placePower in hours;is a hydrogen generating module connected with the power balance nodeIn the first placePower in hours; the power flowing into the power balance node is positive, and the power flowing out of the power balance node is negative;
the hydrogen equilibrium node equation is:
in the formula:is the other hydrogen balance node connected with the hydrogen balance nodeIn the first placeHydrogen mass exchange between hours and the hydrogen balance node;is an electric hydrogen production module connected with the hydrogen balance nodeIn the first placeHydrogen production in hours;is a hydrogen generating module connected with the hydrogen balance nodeIn the first placeHydrogen consumption in hours;is a hydrogen input module connected with the hydrogen balance nodeIn the first placeHydrogen output per hour;is the hydrogen load connected to the hydrogen balance nodeModuleIn the first placeHydrogen usage in hours;is a hydrogen output module connected with the hydrogen balance nodeIn the first placeHydrogen input per hour;is a hydrogen storage module connected with the hydrogen balance nodeIn the first placeHydrogen charge/discharge amount per hour; the amount of hydrogen flowing into the hydrogen balance node is positive and the amount of hydrogen flowing out of the hydrogen balance node is negative;
in S5, the system module constraint equations comprise a power module constraint equation, an electricity storage module constraint equation, an electricity hydrogen production module constraint equation, a hydrogen generation module constraint equation, a hydrogen load module constraint equation and a hydrogen storage module constraint equation;
the power module constraint equation comprises power generation power constraint and renewable energy power curtailment constraint;
the generated power constraint is as follows:
in the formula:non-renewable power moduleIn the first placeHourly generated power;non-renewable power moduleInstalled capacity of (d);
in the formula:is a renewable power moduleIn the first placeHourly generated power;is a renewable power moduleIn the first placeElectric power discard per hour;is a renewable power moduleInstalled capacity of (d);is a renewable power moduleIn the first placeThe ratio of maximum hourly output to installed capacity;
the renewable energy power abandon rate constraint is as follows:
the electric storage module constraint equation comprises charge/energy storage state constraint, charge and discharge power constraint and capacity balance constraint;
the state of charge/energy storage constraints are:
in the formula:andare respectively an electricity storage moduleLower and upper limits of the state of charge/energy storage;is an electricity storage moduleIn the first placeHourly charge/energy storage state;
the charge and discharge power constraint is as follows:
in the formula:is an electricity storage moduleThe capacity of (a);is an electricity storage moduleMaximum charge rate of;is an electricity storage moduleMaximum discharge rate of (d);is an electricity storage moduleThe charging and discharging power of (1);
the capacity balance constraint is:
In the formula:is an electricity storage moduleEnergy conversion efficiency of charging;is an electricity storage moduleEnergy conversion efficiency of the discharge;is an electricity storage moduleA ratio of leakage to capacity per hour;is an electricity storage moduleInitial state of charge/energy storage;
the constraint equation of the electrical hydrogen production module comprises an operation power constraint and an energy conversion constraint;
the operating power constraint is:
in the formula:is an electric hydrogen production moduleInstalled capacity of (d);is an electric hydrogen production moduleThe operating power of (c);
the energy conversion constraint is:
in the formula:is an electric hydrogen production moduleElectricity consumption per unit hydrogen production;is an electric hydrogen production moduleThe amount of energy conversion of (a);
the hydrogen power generation module constraint equation comprises an operating power constraint and an energy conversion constraint;
the operating power constraint is:
the energy conversion constraint is:
in the formula:is a hydrogen power generation moduleHydrogen consumption per unit of generated energy;is a hydrogen power generation moduleThe amount of energy conversion of (a);
the hydrogen load module constraint equation comprises a hydrogen usage constraint;
the hydrogen usage constraint is:
in the formula:is a hydrogen load moduleThe maximum possible amount of hydrogen used per year within said time section;is a hydrogen load moduleAnnual production/traffic within the time section;is a hydrogen load moduleAdopting the carbon emission factor of the prior art;is a hydrogen load moduleThe proportion of emission reduction is realized by hydrogen substitution in the existing process carbon emission;is a hydrogen load moduleCarbon reduction per unit hydrogen application;;
the hydrogen storage module constraint equation comprises hydrogen storage state constraint, hydrogen charging and discharging rate constraint and capacity balance constraint;
the hydrogen storage state constraints are:
in the formula:andare each a hydrogen storage moduleThe lower limit and the upper limit of the hydrogen storage state of (a);is a hydrogen storage moduleIn the first placeHydrogen storage status of hours;
the hydrogen charging and discharging rate constraint is as follows:
in the formula:is a hydrogen storage moduleThe capacity of (c);is a hydrogen storage moduleMaximum hydrogen charge rate of;is a hydrogen storage moduleThe maximum hydrogen release rate;;
the capacity balance constraint is:
In the formula:is a hydrogen storage moduleThe efficiency of the charging process;is a hydrogen storage moduleThe efficiency of the hydrogen discharge process;is a hydrogen storage moduleThe ratio of hydrogen leakage to capacity per hour;is a hydrogen storage moduleThe initial hydrogen storage state of (a);
in S5, the system design constraint equation comprises system carbon emission constraints, upper and lower limit constraints of installed capacity of a power supply module, upper and lower limit constraints of annual input electric quantity of a power input module, upper and lower limit constraints of annual output electric quantity of a power output module, upper and lower limit constraints of capacity of an electricity storage module, upper and lower limit constraints of installed capacity of an electrical hydrogen production module, upper and lower limit constraints of installed capacity of a hydrogen generation module, upper and lower limit constraints of annual input hydrogen quantity of a hydrogen input module, upper and lower limit constraints of annual hydrogen quantity of a hydrogen load module, upper and lower limit constraints of annual output hydrogen quantity of a hydrogen output module, and upper and lower limit constraints of capacity of a hydrogen storage module;
wherein the system carbon emission constraints are:
in the formula:is a power supply moduleCarbon emission factor of (a);is a power input moduleCarbon emission factor of (a);is a power output moduleThe carbon emission factor of (c);is a hydrogen input moduleCarbon emission factor of (a);is a hydrogen output moduleCarbon emission factor of (a);is the upper carbon emission limit of the regional electro-hydrogen system at the time section;efficiency;is a power input moduleCarbon emission efficiency of;is a power output moduleCarbon emission efficiency of;is a hydrogen input moduleCarbon emission efficiency of;is a hydrogen output moduleCarbon emission efficiency of;is a hydrogen load moduleAnnual production/traffic in time section;is a hydrogen load moduleAdopting the carbon emission factor of the prior art;is a hydrogen load moduleThe proportion of emission reduction can be realized by replacing hydrogen in the carbon emission in the prior art;is a hydrogen load moduleThe amount of carbon reduction per unit of hydrogen application.
2. The regional electric hydrogen system decarburization path planning optimization method according to claim 1, wherein in S1, the regional electric hydrogen system (100) comprises a power supply unit (1), a power input unit (2), a power balancing unit (3), a power load unit (4), an electricity storage unit (5), an electricity hydrogen production unit (6), a hydrogen generation unit (7), a hydrogen input unit (8), a hydrogen storage unit (9), a hydrogen balancing unit (10), a hydrogen load unit (18), a power output unit (19) and a hydrogen output unit (20);
the power supply unit (1), the power input unit (2), the power output unit (19), the power load unit (4), the electricity storage unit (5), the electricity hydrogen production unit (6), the hydrogen power generation unit (7) all link to each other with the power balance unit (3), the electricity hydrogen production unit (6), the hydrogen power generation unit (7), the hydrogen input unit (8), the hydrogen storage unit (9), the hydrogen output unit (20), the hydrogen load unit (18) all link to each other with the hydrogen balance unit (10).
3. The regional electric hydrogen system decarburization path planning optimization method according to claim 2, wherein in S1, the number of the power balance units (3) is one or more, a plurality of the power balance units (3) are connected with each other, and the power balance units (3) are connected with each other in a bidirectional or unidirectional manner; the number of the hydrogen balance units (10) is one or more, the hydrogen balance units (10) are connected with each other, and bidirectional or unidirectional hydrogen transmission exists between the hydrogen balance units (10) which are connected with each other.
4. The regional power and hydrogen system decarburization path planning optimization method of claim 1, wherein in S2, the planning period spans a plurality of consecutive natural years, a single time section is a natural year in the planning period, and the time sections are one or more.
5. The method for optimizing a regional electro-hydrogen system decarburization path planning of claim 1, wherein in S3, the variables to be optimized of the regional electro-hydrogen system (100) include real variables and one-dimensional array variables.
6. The regional electric hydrogen system decarburization path planning optimization method according to claim 5, wherein the real type variables include installed capacity of each power module, annual input power of each power input module, annual output power of each power output module, capacity of each power storage module, installed capacity of each electric hydrogen production module, installed capacity of each hydrogen generation module, annual input hydrogen of each hydrogen input module, annual hydrogen usage of each hydrogen load module, annual output hydrogen of each hydrogen output module, and capacity of each hydrogen storage module for a selected time section in the regional electric hydrogen system (100);
each of the one-dimensional array variables includes 8760 real-type variables, and the one-dimensional array variables include a hourly power generation power of each power module, a hourly discarded electric power of each power module, a hourly charge and discharge electric power of each power module, a hourly charge/energy storage state of each power storage module, a hourly operation power of each hydrogen production module, a hourly hydrogen production rate of each hydrogen production module, a hourly power generation power of each hydrogen generation module, a hourly hydrogen consumption rate of each hydrogen generation module, a hourly charge and discharge hydrogen rate of each hydrogen storage module, and a hourly hydrogen storage state of each hydrogen storage module in a selected time section of the regional hydrogen power generation system (100).
7. The regional electro-hydrogen system decarburization path planning and optimizing method according to claim 1, wherein in S6, the configuration optimization algorithm selects a linprog solver in MATLAB according to the optimization objective function and the constraint condition equation.
8. The regional electric hydrogen system decarburization path planning optimization method according to claim 1, wherein in S7, the input parameters of the regional electric hydrogen system (100) of the current time section include initial existing installed capacity of each power module, upper and lower limits of installed capacity, carbon emission factor, unit fixed investment, average design life, ratio of annual fixed operation and maintenance cost to fixed investment, electric fuel cost, ratio of hourly maximum output of renewable energy to installed capacity;
the input parameters of the regional electro-hydrogen system (100) of the current time section further comprise a ratio of a time-wise delivered power to a peak delivered power of the power input module;
the input parameters of the regional electro-hydrogen system (100) of the current time section further comprise the time-by-time power of the power load module;
the input parameters of the regional electro-hydrogen system (100) of the current time section further comprise a ratio of a time-wise delivered power to a peak delivered power of the power output module;
the input parameters of the regional hydrogen generation system (100) of the current time section further comprise the initial existing installed capacity of each electricity storage module, the upper and lower limits of the installed capacity, the maximum charging multiplying factor, the maximum discharging multiplying factor, the energy conversion efficiency in the charging process, the energy conversion efficiency in the discharging process, the ratio of the hourly leakage amount to the capacity, the unit fixed investment, the average design life and the ratio of the annual fixed operation and maintenance cost to the fixed investment;
the input parameters of the regional hydrogen generation system (100) of the current time section further comprise the initial existing installed capacity of the hydrogen generation module, the upper and lower limits of the installed capacity, the power consumption of unit hydrogen production, unit fixed investment, average design life and the ratio of annual fixed operation and maintenance cost to fixed investment;
the input parameters of the regional hydrogen power system (100) of the current time section further comprise the initial existing installed capacity of each hydrogen power generation module, the upper limit and the lower limit of the installed capacity, the hydrogen consumption of unit generated energy, unit fixed investment, average design life and the ratio of annual fixed operation and maintenance cost to fixed investment;
the input parameters of the regional electric hydrogen system (100) of the current time section further comprise the ratio of the hourly hydrogen output to the peak hydrogen output of the hydrogen input module;
the input parameters of the regional electro-hydrogen system (100) of the current time section further comprise the ratio of the hourly hydrogen consumption to the peak hydrogen consumption of each hydrogen load module, the annual output/transport capacity, a carbon emission factor adopting the prior art, the ratio of realizing emission reduction by hydrogen substitution in the carbon emission of the prior art, the carbon reduction amount brought by each unit of hydrogen application and the ratio of the hourly hydrogen consumption to the peak hydrogen consumption of the hydrogen output module;
the input parameters of the regional electric hydrogen system (100) of the current time section further comprise the ratio of the hourly hydrogen output quantity to the peak hydrogen output quantity of each hydrogen output module;
the input parameters of the regional hydrogen power system (100) of the current time section further comprise initial existing installed capacity of each hydrogen storage module, upper and lower limits of the installed capacity, maximum hydrogen charging rate, maximum hydrogen discharging rate, efficiency of a hydrogen charging process, efficiency of a hydrogen discharging process, a ratio of hydrogen leakage amount per hour to capacity, unit fixed investment, average design life and a ratio of annual fixed operation and maintenance cost to fixed investment.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210856118.9A CN115079564B (en) | 2022-07-21 | 2022-07-21 | Decarburization path planning optimization method for regional hydrogen generation system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210856118.9A CN115079564B (en) | 2022-07-21 | 2022-07-21 | Decarburization path planning optimization method for regional hydrogen generation system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115079564A CN115079564A (en) | 2022-09-20 |
CN115079564B true CN115079564B (en) | 2022-11-04 |
Family
ID=83258868
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210856118.9A Active CN115079564B (en) | 2022-07-21 | 2022-07-21 | Decarburization path planning optimization method for regional hydrogen generation system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115079564B (en) |
Family Cites Families (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2640910T3 (en) * | 2000-10-27 | 2017-11-07 | Air Products And Chemicals, Inc. | Systems and processes to provide hydrogen to fuel cells |
US7803489B2 (en) * | 2007-03-26 | 2010-09-28 | Advanced Hydrogen Power Systems, Inc. | Hydrogen mobile power plant that extracts hydrogen fuel from water |
US7888540B2 (en) * | 2008-04-11 | 2011-02-15 | General Electric Company | Integrated system and method for producing fuel composition from biomass |
CN101832467B (en) * | 2010-06-03 | 2011-09-14 | 浙江大学 | Hydrogen refueling station hydrogen optimizing filling method based on temperature rise and energy consumption control and system thereof |
CN101958423A (en) * | 2010-09-25 | 2011-01-26 | 华为技术有限公司 | Hydrogen supply system, system and method for supplying hydrogen and hydrogen fuel cell system |
EP2695855B1 (en) * | 2012-08-08 | 2017-11-01 | Airbus Defence and Space GmbH | Hydrogen generator |
CN103346562B (en) * | 2013-07-11 | 2015-06-17 | 江苏省电力设计院 | Multi-time scale microgrid energy control method considering demand response |
CN109179320B (en) * | 2018-11-08 | 2020-02-07 | 新地能源工程技术有限公司 | Natural gas on-site hydrogen production device and method |
CN109494809B (en) * | 2019-01-07 | 2020-12-29 | 中国矿业大学(北京) | Electric-gas interconnection system comprising electric gas conversion device and operation method |
CN109741110B (en) * | 2019-01-07 | 2022-07-01 | 福州大学 | Opportunity constraint planning-based wind-hydrogen system joint optimization modeling method |
CN110020506B (en) * | 2019-03-08 | 2023-04-07 | 东南大学 | Differential format selection method based on operation optimization of electric heating type comprehensive energy system |
CN111027846B (en) * | 2019-12-06 | 2023-06-16 | 东南大学溧阳研究院 | Electricity-hydrogen comprehensive energy system considering heat-hydrogen co-production and capacity configuration method thereof |
CN111383719B (en) * | 2019-12-23 | 2023-10-20 | 浙江中控技术股份有限公司 | Synergistic optimization method for hydrogen and gas in oil refining process |
CN112366728B (en) * | 2020-10-27 | 2022-02-01 | 清华四川能源互联网研究院 | Combined frequency modulation testing method, device and control system for battery energy storage and hydrogen production equipment |
CN112952805A (en) * | 2021-02-08 | 2021-06-11 | 国网江苏省电力有限公司电力科学研究院 | Electric hydrogen energy system scheduling method considering flexible hydrogen demand |
CN112989594B (en) * | 2021-03-09 | 2023-05-19 | 南方电网科学研究院有限责任公司 | Comprehensive energy system operation optimization method considering hydrogen energy |
CN114240010A (en) * | 2021-07-26 | 2022-03-25 | 国网山西省电力公司经济技术研究院 | Scheduling method for cogeneration gas turbine unit coping with uncertainty of demand side |
CN114065483A (en) * | 2021-10-20 | 2022-02-18 | 清华大学 | Hydrogen production system efficiency optimization method and device, computer equipment and storage medium |
CN114565204A (en) * | 2021-11-23 | 2022-05-31 | 重庆大学 | Electricity-hydrogen coordination planning method for distributed double-layer hydrogen supply network and power distribution network |
CN114386263A (en) * | 2021-12-31 | 2022-04-22 | 上海电力大学 | Medium-term and long-term planning method of comprehensive energy system considering carbon emission reduction constraint and hydrogen energy |
CN114298440A (en) * | 2022-01-04 | 2022-04-08 | 西安交通大学 | Supply and demand collaborative operation optimization method and control device for hydrogen-containing multi-energy system of building |
CN114386839A (en) * | 2022-01-12 | 2022-04-22 | 国网浙江省电力有限公司 | Method, system, medium and equipment for evaluating economy of electric-hydrogen complementary energy system |
CN114757052A (en) * | 2022-05-05 | 2022-07-15 | 大连理工大学 | Method for evaluating whole life cycle economy of PEM (proton exchange membrane) electric hydrogen production system |
-
2022
- 2022-07-21 CN CN202210856118.9A patent/CN115079564B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN115079564A (en) | 2022-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Parra et al. | A review on the role, cost and value of hydrogen energy systems for deep decarbonisation | |
Jaramillo et al. | Optimal microgrid scheduling with peak load reduction involving an electrolyzer and flexible loads | |
Chua et al. | Cost-benefit assessment of energy storage for utility and customers: A case study in Malaysia | |
US20110106322A1 (en) | Network connection manner of microgrid energy storage backup power source, and method for dispatching the same | |
García et al. | Energy management system based on techno-economic optimization for microgrids | |
Rad et al. | Excess electricity problem in off-grid hybrid renewable energy systems: A comprehensive review from challenges to prevalent solutions | |
Mladenov et al. | Characterisation and evaluation of flexibility of electrical power system | |
Nozari et al. | Development of dynamic energy storage hub concept: A comprehensive literature review of multi storage systems | |
Paniyil et al. | Batteries and free fuel based photovoltaics and complimentary wind energy based DC power networks as 100% source of electric power around the globe | |
CN117973886A (en) | Comprehensive energy system collaborative planning operation method and system for hydrogen-containing energy full link | |
Jin et al. | Coordinated usage of distributed sources for energy cost saving in micro-grid | |
Siti Khodijah et al. | Design of compensation battery for tidal power‐photovoltaics‐SOFC microgrids in Ternate and Pulau‐Tidore Islands | |
CN115079564B (en) | Decarburization path planning optimization method for regional hydrogen generation system | |
CN110648031A (en) | Energy router optimal operation scheduling method based on reverse order dynamic programming | |
Kudria et al. | Hydrogen technologies as a method of compensation for inequality of power generation by renewable energy sources | |
Liu et al. | A power-to-gas integrated microgrid optimal operation strategy based on rolling horizon | |
Langella et al. | Analysis of the cost of electric energy discharged to the grid by some energy storage systems | |
CN112054553A (en) | Coordinated optimization operation method, system, medium and equipment for electric-heat-gas interconnection system | |
Zigurs et al. | National Energy and Climate Plans: Importance of Synergy | |
Tebibel et al. | Levelized Cost of Hydrogen Analysis for Alkaline Water Electrolysis-Photovoltaic Energy Technology in The Near Future (2025–2050) of Algeria | |
Mu et al. | An Operation Scheduling Model for Carbon Neutrality in Industrial Integrated Energy System | |
Seedahmed et al. | Control-oriented model of an optimally designed hybrid storage system for a standalone microgrid | |
CN118214091B (en) | Sea island micro-grid layered coordination method based on electric-hydrogen hybrid energy storage | |
Kiviluoma et al. | Analysing wind and solar power integration with a multi-carrier energy system model of ireland | |
Sevgen | Optimal sizing and location on energy storage systems |
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 |