A multi-garden integrated energy system operation method considering hydrogen energy optimization scheduling

Abstract

The application provides a multi-garden comprehensive energy system operation method considering hydrogen energy optimization scheduling. The method is based on the detailed model of each garden and the characteristics of the hydrogen energy long pipe trailer, and comprehensively considers the energy balance in the garden and the hydrogen energy scheduling optimization problem between gardens. Specifically, first, hydrogen is prepared by using new energy power generation, and hydrogen is used instead of traditional coke for ironmaking, and an innovative new energy steel garden model is constructed. In view of the insufficient hydrogen preparation in the steel garden and the insufficient power consumption capacity of the surrounding new energy power generation garden, a multi-garden hydrogen energy collaborative scheduling scheme is proposed. In the hydrogen transportation process, by constraining the hydrogen storage capacity, transportation path and time of the hydrogen energy long pipe trailer, the efficient distribution of hydrogen energy between gardens is realized. The method reduces carbon emissions by using hydrogen instead of coke for ironmaking, effectively improves the power consumption capacity of the surrounding garden, optimizes the energy utilization efficiency, and significantly reduces the overall operation cost of the system.

Description

Technical Field

[0001] The present invention belongs to the technical field of hydrogen energy optimization scheduling of integrated energy systems, and in particular relates to a multi-park integrated energy system operation method considering hydrogen energy optimization scheduling. Background Art

[0002] As a foundational industry of modern industry, the steel industry plays a vital role in global economic development. However, traditional steel production processes rely heavily on carbon-reduced iron (CRI), a process that generates significant emissions of carbon dioxide and other pollutants. Hydrogen, a replacement for CRI, is emerging as a key path for the steel industry to achieve a green transformation. Hydrogen, a clean energy carrier, is used to reduce iron ore in steel production. Its primary reaction product is water, which is expected to achieve near-zero carbon emissions compared to traditional processes, offering a promising solution for the steel industry's green transformation.

[0003] However, relying solely on the steel park's own electricity generation and hydrogen production still faces many challenges, such as tight energy supply, limited site space, and high hydrogen production costs. In addition, although direct purchase of hydrogen can solve some problems, its economic efficiency is poor. At present, there have been some studies on the energy optimization scheduling of a single steel park, but there are relatively few studies on the coordinated optimization scheduling between multiple parks with different functions and energy needs. In particular, in the areas of new energy hydrogen production, hydrogen energy storage and transportation, and the coupled utilization of multiple energy forms, there are still many technical difficulties that need to be solved, including how to ensure the balance of energy supply and demand in each park, and how to optimize the hydrogen transportation route and time to improve energy utilization efficiency.

[0004] Based on the above description, the present invention proposes a multi-park integrated energy system operation method considering hydrogen energy optimization scheduling, focusing on the following aspects: First, constructing a detailed energy model of the new energy steel park, covering key links such as production capacity, energy conversion, storage and use, to achieve refined energy management; Second, on the basis of the new energy steel park model, combined with other new energy parks, the excess electricity generated by them is used to provide energy for the steel park through electrolysis to produce hydrogen, thereby achieving energy complementarity and deep coupling between parks; Third, in terms of hydrogen energy transportation, the transportation route and time are comprehensively optimized, and hydrogen long-tube trailers are used to efficiently transport hydrogen to ensure that hydrogen is supplied to the steel park and hydrogen refueling station in a timely manner, thereby improving the efficiency of hydrogen energy utilization. Summary of the Invention

[0005] Purpose of the invention: The present invention first proposes an innovative green hydrogen ironmaking steel park model, which uses photovoltaic power generation to power the electrolytic cell to produce green hydrogen, and applies hydrogen to steel production to replace the traditional carbon reduction iron process, thereby effectively restraining carbon emissions in the steel production and park power supply links. On this basis, the present invention further proposes a multi-park integrated energy system operation method that takes into account hydrogen energy optimization scheduling. By unified planning and coordinated scheduling of energy production, conversion, storage, transmission and consumption of multiple new energy parks (including steel parks), it can effectively absorb the excess energy of new energy parks, improve energy utilization efficiency, reduce energy costs, and reduce carbon emissions. In addition, in the hydrogen energy transmission link of the park, a hydrogen long-tube trailer is introduced and its operation path is optimized. At the same time, the coupling relationship between hydrogen production, hydrogen transmission and hydrogen use processes in the time dimension is deeply analyzed to achieve deep energy coupling and coordinated operation between multiple parks.

[0006] Technical solution: In order to solve the above technical problems, the present invention proposes a multi-park integrated energy system operation method considering hydrogen energy optimization scheduling, which includes the following steps:

[0007] Step 1: Obtain the operating parameters of the equipment in each park. For the new energy photovoltaic steel park, this includes the installed capacity of photovoltaic generators, electrolyzer capacity, fuel cell efficiency, energy storage equipment capacity, and charge and discharge loss efficiency; for the new energy park, this includes the installed capacity of wind turbines, photovoltaic power generation capacity, and equipment parameters of electrolyzers; for hydrogen refueling stations, this includes the hydrogen demand for each time period and the capacity of hydrogen storage tanks;

[0008] Step 2: Obtain load demand data for each park, including power load, heat load, and steel demand, and collect geographical distances between parks and wind power output information;

[0009] Step 3: Based on the above information, construct each park with energy supply and demand balance constraints and equipment operation constraints as constraints;

[0010] Step 4. Based on the park constraints and combined with the hydrogen demand of the hydrogen refueling station, a hydrogen energy optimization scheduling model for the multi-park integrated energy system is constructed. The model is restricted by the overall hydrogen energy supply and demand balance, the capacity constraint of the hydrogen tube trailer for hydrogen transportation, the path planning constraint and the time constraint. The objective function is to minimize the overall energy purchase cost of the multi-park and the operating cost of the hydrogen tube trailer. The model is solved using the GUROBI solver to obtain the operation plan of the equipment in each park, including the power generation power of the power generation equipment, the hydrogen production power of the electrolyzer, the charging and discharging power of the energy storage equipment, and the path and time arrangement of hydrogen transportation, so as to achieve the optimal scheduling of the multi-park integrated energy system.

[0011] Furthermore, in step 3, the relevant operation constraints include energy supply and demand balance constraints, which are as follows:

[0012]

[0013]

[0014] Where, subscript t represents the scheduling period, subscript pa represents all park nodes, subscript ic represents the steel park, subscript R represents the new energy hydrogen production park, subscript sD represents the energy storage element, superscript pv represents the photovoltaic generator set, superscript el represents the electrolyzer, superscript fc represents the fuel cell, and superscript gb represents the boiler; P ch e,ic,t With P ch h,ic,t They represent the electricity and hydrogen reserves of the ic park at time t, P ch g,pa,t P represents the power of methane storage in the park at time t, dis e,ic,t With P dis h,ic,t They represent the capacity of electricity and hydrogen released by the IC park at time t, P dis g,pa,t The power of methane gas released by the pa park at time t, P pv ic,t 、P ce ic,t With P fc ic,t They represent the photovoltaic power generation of IC park, additional power purchase and hydrogen fuel cell power generation at time t, respectively. ch ic,t Indicates the hydrogen production state of the electrolyzer in the ic park at time t. If it is in operation state, the value is 1, otherwise it is 0. el ic,t Indicates the power consumption of the electrolytic cell in the ic park at time t, P l ic,t represents the power load of the ic park at time t, P pvmin With P pvmax Respectively represent the maximum and minimum values ​​of photovoltaic output, M el ic,t and M el R,t Respectively represent the hydrogen production of ic and R park electrolyzer at time t, M fc ic,t 、M car t With M dc t They represent the hydrogen consumption of hydrogen fuel cells, new energy hydrogen vehicles for steel transportation, and hydrogen consumption for hydrogen ironmaking, respectively. ch ic,t represents the amount of hydrogen purchased by IC Park at time t, Q gbic,t , Q fc ic,t With Q l ic,t They represent the heat output and heat load demand of the IC park boiler and fuel cell at time t, M cg pa,t With M gb pa,t They represent the gas purchase volume of the park and the methane demand of the boiler at time t, P ch h,R,t 、P dis h,R,t They represent the amount of hydrogen stored and released in R Park at time t, M give R,t represents the excess hydrogen amount in R Park at time t, P ce pa,t 、M ch pa,t 、M cg pa,t They represent the additional electricity purchase, additional hydrogen purchase and additional gas purchase of the park at time t, P cemax 、M chmax With M cgmax They represent the upper limits of electricity, hydrogen and gas purchases from the power grid respectively.

[0015] Furthermore, in step 3, the relevant operation constraints include electrolytic cell operation constraints, which are as follows:

[0016]

[0017] Where, subscript v represents all parks in the interconnected system, M el pa，t represents the hydrogen production of the electrolyzer in the pa park at time t, η el represents the electrolysis efficiency, X ch pa,t Indicates the hydrogen production status of the electrolyzer in the pa park at time t. The value is 1 if the operation status is running, otherwise it is 0. el pa,t P represents the power consumption of the electrolytic cell in the park at time t, elmax Indicates the maximum allowable input power of the electrolyzer.

[0018] Furthermore, in step 3, the relevant operating constraints include energy storage equipment operating constraints, which are as follows:

[0019]

[0020] Where S sD Indicates the rated capacity of the energy storage device, U ch pa,t and U dispa,t They represent whether the park pa stores or releases energy at time t. A value of 1 indicates that energy storage or release begins, while a value of 0 indicates that energy is not stored or released. ch sD,pa,t With P dis sD,pa,t They represent the energy stored and released in the park at time t, P ss sD,pa,t represents the actual capacity of the energy storage equipment in the park at time t, P ssmin sD With P ssmax sD Respectively represent the actual maximum and minimum capacity of the energy storage device, λ ch sD and λ dis sD They represent the efficiency of energy storage and release respectively.

[0021] Furthermore, in step 4, the minimum objective function is:

[0022]

[0023] Where T represents the time section number, subscript D represents the hydrogen refueling station, subscript K represents the set of hydrogen tube trailers, and subscript k represents the specific hydrogen tube trailer vehicle number; C e t represents the time-of-use electricity price, C h Indicates the price of industrial hydrogen per kilogram, C g Indicates the price per kilogram of methane, L R,D , L D,ic , L R,ic A represents the distance between the new energy park R and the hydrogen station D, the hydrogen station D and the steel park ic, and the new energy park R and the steel park ic. R,D,k,t Indicates the running status of vehicle k at time t. If the vehicle starts from the new energy park R and goes to the hydrogen refueling station D, the value is 1, otherwise it is 0. D,ic,k,t 、A R,ic,k,t Indicates the running status of the vehicle from D and R to the steel park ic. If it is in operation, the value is 1, otherwise it is 0. car Indicates vehicle driving expenses.

[0024] Furthermore, in step 4, the hydrogen transport constraints of hydrogen tube trailers, vehicle routing constraints and time constraints, and multi-park hydrogen supply and demand balance constraints are as follows:

[0025] (1) Hydrogen long tube trailer route constraints

[0026]

[0027] (2) Constraints on hydrogen transport by hydrogen tube trailers

[0028]

[0029] (3) Hydrogen delivery time coupling constraints

[0030]

[0031] (4) Conservation constraints on hydrogen supply and demand

[0032]

[0033] Where A i,j,k,t Indicates the running status of the hydrogen vehicle. If vehicle k moves from park i to park j at time t, the value is 1, otherwise it is 0. max i H represents the maximum number of hydrogen delivery vehicles owned by park i. need D,k,t represents the amount of hydrogen supplied to hydrogen refueling station D by vehicle k at time t, M h i,k,t M represents the amount of hydrogen carried by the kth vehicle when leaving the i-th park. h R,k,t 、M h D,k,t M represents the amount of hydrogen carried by the kth vehicle leaving the new energy park R and the hydrogen refueling station D at time t, respectively. hmax Indicates the maximum amount of hydrogen carried by the hydrogen tube trailer, S carH k,t H represents the amount of hydrogen supplied to the park ic by the kth vehicle at time t. give R,k,t Indicates the amount of hydrogen actually delivered to the steel park at time t, represents t-Δt give R,k,t -Δt R,D,k,t -Δt D,ic,k,t -Δt R,ic,k,t The amount of hydrogen supplied to the steel park by the kth vehicle at time Δt i,j,k,t represents the time it takes for the kth vehicle to reach park j from park i, Δt R,D,k,t , Δt D,ic,k,t , Δt R,ic,k,t They represent the time consumed by the kth vehicle from the new energy park R to the hydrogen station D, the hydrogen station D to the steel park ic, and the new energy park R to the steel park ic at time t, respectively. leave i,k,t and t arrive j,k,t They represent the time when the kth vehicle leaves park i and arrives at park j at time t, Δt give R,k,tIndicates the time consumed from R to the steel park, ceil indicates the use of rounded up counting, t arrive i,k,t represents the time when the kth vehicle arrives at park i at time t, t arrivemax t , t arrivemin t , t leavemax t and t leavemin t Indicates the earliest or latest time of arrival or departure from the park within the time period t, t leave R,k,t represents the specific time when vehicle k leaves R park at time t, t arrive R,k,t-1 represents the specific time when vehicle k arrives at park R at time t-1, v h2 Indicates the hydrogen flow rate, which is assumed to be a constant value here.

[0034] Furthermore, in step 4, a new energy park model is constructed through constraints (A-1) to (A-16), (A-25), (A-27), and (A-33), and the objective function (A-17) is solved. By calculating the power generation power of the power generation equipment, the hydrogen production power of the electrolyzer, and the charging and discharging power of the energy storage equipment, the minimum total energy purchase cost of all parks is finally solved; a hydrogen tube trailer model is constructed through constraints (A-19) to (A-33), and the objective function (A-18) is solved. By calculating the route and specific time of hydrogen transportation, the minimum hydrogen tube trailer transportation cost is obtained.

[0035] Beneficial effects: Compared with the prior art, the technical solution of the present invention has the following beneficial technical effects:

[0036] The present invention plans to promote photovoltaic power generation and hydrogen production for steel production in new energy photovoltaic steel parks, and introduce fuel cells to generate electricity and heat, thereby reducing the use of traditional fossil energy such as coal and reducing carbon emissions. At the same time, fuel vehicles used to transport steel are replaced with new energy trucks that transport steel, further reducing carbon emissions in the transportation process. Through collaborative optimization and scheduling among multiple parks, the complementary utilization of different energy forms is achieved, energy conversion and utilization efficiency is improved, energy waste is reduced, and energy costs are reduced. Finally, the time coupling constraints of hydrogen production, transportation, and use, as well as the operating constraints of each device, are considered to ensure the stability and reliability of energy supply and improve the ability of the entire integrated energy system to respond to emergencies. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] Figure 1 It is a flow chart of the method of the present invention.

[0038] Figure 2 It is a flow chart of the hydrogen energy transmission path between the parks.

[0039] Figure 3 This is a comparison chart of carbon emissions from ironmaking using different reducing agents.

[0040] Figure 4 This is a chart showing changes in energy purchases before the steel park is interconnected with other parks.

[0041] Figure 5 This is a chart showing changes in energy purchases after the steel park is interconnected with other parks. DETAILED DESCRIPTION

[0042] The present invention is further illustrated below with reference to the accompanying drawings and specific embodiments. It should be understood that these embodiments are only used to illustrate the present invention and are not used to limit the scope of the present invention. After reading the present invention, various equivalent modifications of the present invention by those skilled in the art fall within the scope defined by the claims attached to this application.

[0043] like Figure 1 As shown, the present invention proposes a multi-park integrated energy system operation method considering hydrogen energy optimization scheduling, which includes the following steps:

[0044] Step 1: Obtain the operating parameters of the equipment in each park. For the new energy photovoltaic steel park, this includes the installed capacity of photovoltaic generators, electrolyzer capacity, fuel cell efficiency, energy storage equipment capacity, and charge and discharge loss efficiency; for the new energy park, this includes the installed capacity of wind turbines, photovoltaic power generation capacity, and equipment parameters of electrolyzers; for hydrogen refueling stations, this includes the hydrogen demand for each time period and the capacity of hydrogen storage tanks;

[0045] Step 2: Obtain load demand data for each park, including power load, heat load, and steel demand, and collect geographical distances between parks and wind power output information;

[0046] Step 3: Based on the above information, construct each park with energy supply and demand balance constraints and equipment operation constraints as constraints;

[0047] Step 4. Based on the park constraints and combined with the hydrogen demand of the hydrogen refueling station, a hydrogen energy optimization scheduling model for the multi-park integrated energy system is constructed. The model is restricted by the overall hydrogen energy supply and demand balance, the capacity constraint of the hydrogen tube trailer for hydrogen transportation, the path planning constraint and the time constraint. The objective function is to minimize the overall energy purchase cost of the multi-park and the operating cost of the hydrogen tube trailer. The model is solved using the GUROBI solver to obtain the operation plan of the equipment in each park, including the power generation power of the power generation equipment, the hydrogen production power of the electrolyzer, the charging and discharging power of the energy storage equipment, and the path and time arrangement of hydrogen transportation, so as to achieve the optimal scheduling of the multi-park integrated energy system.

[0048] Furthermore, in step 3, the relevant operation constraints include energy supply and demand balance constraints, which are as follows:

[0049]

[0050]

[0051] Where, subscript t represents the scheduling period, subscript pa represents all park nodes, subscript ic represents the steel park, subscript R represents the new energy hydrogen production park, subscript sD represents the energy storage element, superscript pv represents the photovoltaic generator set, superscript el represents the electrolyzer, superscript fc represents the fuel cell, and superscript gb represents the boiler; P ch e,ic,t With P ch h,ic,t They represent the electricity and hydrogen reserves of the ic park at time t, P ch g,pa,t P represents the power of methane storage in the park at time t, dis e,ic,t With P dis h,ic,t They represent the capacity of electricity and hydrogen released by the IC park at time t, P dis g,pa,t The power of methane gas released by the pa park at time t, P pv ic,t 、P ce ic,t With P fc ic,t They represent the photovoltaic power generation of IC park, additional power purchase and hydrogen fuel cell power generation at time t, respectively. ch ic,t Indicates the hydrogen production state of the electrolyzer in the ic park at time t. If it is in operation state, the value is 1, otherwise it is 0. el ic,t Indicates the power consumption of the electrolytic cell in the ic park at time t, P l ic,t represents the power load of the ic park at time t, P pvmin With P pvmax Respectively represent the maximum and minimum values ​​of photovoltaic output, M el ic,t and M el R,t Respectively represent the hydrogen production of ic and R park electrolyzer at time t, M fc ic,t 、M car t With M dc t They represent the hydrogen consumption of hydrogen fuel cells, new energy hydrogen vehicles for steel transportation, and hydrogen consumption for hydrogen ironmaking, respectively. ch ic,trepresents the amount of hydrogen purchased by IC Park at time t, Q gb ic,t , Q fc ic,t With Q l ic,t They represent the heat output and heat load demand of the IC park boiler and fuel cell at time t, M cg pa,t With M gb pa,t They represent the gas purchase volume of the park and the methane demand of the boiler at time t, P ch h,R,t 、P dis h,R,t They represent the amount of hydrogen stored and released in R Park at time t, M give R,t represents the excess hydrogen amount in R Park at time t, P ce pa,t 、M ch pa,t 、M cg pa,t They represent the additional electricity purchase, additional hydrogen purchase and additional gas purchase of the park at time t, P cemax 、M chmax With M cgmax They represent the upper limits of electricity, hydrogen and gas purchases from the power grid respectively.

[0052] Furthermore, in step 3, the relevant operation constraints include electrolytic cell operation constraints, which are as follows:

[0053]

[0054] Where, subscript v represents all parks in the interconnected system, M el pa,t represents the hydrogen production of the electrolyzer in the pa park at time t, η el represents the electrolysis efficiency, X ch pa,t Indicates the hydrogen production status of the electrolyzer in the pa park at time t. The value is 1 if the operation status is running, otherwise it is 0. el pa,t P represents the power consumption of the electrolytic cell in the park at time t, elmax Indicates the maximum allowable input power of the electrolyzer.

[0055] Furthermore, in step 3, the relevant operating constraints include energy storage equipment operating constraints, which are as follows:

[0056]

[0057] Where S sD Indicates the rated capacity of the energy storage device, U chpa,t and U dis pa,t They represent whether the park pa stores or releases energy at time t. A value of 1 indicates that energy storage or release begins, while a value of 0 indicates that energy is not stored or released. ch sD,pa,t With P dis sD,pa,t They represent the energy stored and released in the park at time t, P ss sD,pa,t represents the actual capacity of the energy storage equipment in the park at time t, P ssmin sD With P ssmax sD Respectively represent the actual maximum and minimum capacity of the energy storage device, λ ch sD and λ dis sD They represent the efficiency of energy storage and release respectively.

[0058] Furthermore, in step 4, the minimum objective function is:

[0059]

[0060] Where T represents the time section number, subscript D represents the hydrogen refueling station, subscript K represents the set of hydrogen tube trailers, and subscript k represents the specific hydrogen tube trailer vehicle number; C e t represents the time-of-use electricity price, C h Indicates the price of industrial hydrogen per kilogram, C g Indicates the price per kilogram of methane, L R,D , L D,ic , L R,ic A represents the distance between the new energy park R and the hydrogen station D, the hydrogen station D and the steel park ic, and the new energy park R and the steel park ic. R,D,k,t Indicates the running status of vehicle k at time t. If the vehicle starts from the new energy park R and goes to the hydrogen refueling station D, the value is 1, otherwise it is 0. D,ic,k,t 、A R,ic,k,t Indicates the running status of the vehicle from D and R to the steel park ic. If it is in operation, the value is 1, otherwise it is 0. car Indicates vehicle driving expenses.

[0061] Furthermore, in step 4, the hydrogen transport constraints of hydrogen tube trailers, vehicle routing constraints and time constraints, and multi-park hydrogen supply and demand balance constraints are as follows:

[0062] (1) Hydrogen long tube trailer route constraints

[0063]

[0064] (2) Constraints on hydrogen transport by hydrogen tube trailers

[0065]

[0066] (3) Hydrogen delivery time coupling constraints

[0067]

[0068] (4) Conservation constraints on hydrogen supply and demand

[0069]

[0070] Where A i,j,k,t Indicates the running status of the hydrogen vehicle. If vehicle k moves from park i to park j at time t, the value is 1, otherwise it is 0. max i H represents the maximum number of hydrogen delivery vehicles owned by park i. need D,k,t represents the amount of hydrogen supplied to hydrogen refueling station D by vehicle k at time t, M h i,k,t M represents the amount of hydrogen carried by the kth vehicle when leaving the i-th park. h R,k,t 、M h D,k,t M represents the amount of hydrogen carried by the kth vehicle leaving the new energy park R and the hydrogen refueling station D at time t, respectively. hmax Indicates the maximum amount of hydrogen carried by the hydrogen tube trailer, S carH k,t H represents the amount of hydrogen supplied to the park ic by the kth vehicle at time t. give R,k,t Indicates the amount of hydrogen actually delivered to the steel park at time t, represents t-Δt give R,k,t -Δt R,D,k,t -Δt D,ic,k,t -Δt R,ic,k,t The amount of hydrogen supplied to the steel park by the kth vehicle at time Δt i,j,k,t represents the time it takes for the kth vehicle to reach park j from park i, Δt R,D,k,t , Δt D,ic,k,t , Δt R,ic,k,t They represent the time consumed by the kth vehicle from the new energy park R to the hydrogen station D, the hydrogen station D to the steel park ic, and the new energy park R to the steel park ic at time t, respectively. leave i,k,t and t arrive j,k,tThey represent the time when the kth vehicle leaves park i and arrives at park j at time t, Δt give R,k,t Indicates the time consumed from R to the steel park, ceil indicates the use of rounded up counting, t arrive i,k,t represents the time when the kth vehicle arrives at park i at time t, t arrivemax t , t arrivemin t , t leavemax t and t leavemin t Indicates the earliest or latest time of arrival or departure from the park within the time period t, t leave R,k,t represents the specific time when vehicle k leaves R park at time t, t arrive R,k,t-1 represents the specific time when vehicle k arrives at park R at time t-1, Indicates the hydrogen flow rate, which is assumed to be a constant value here.

[0071] Furthermore, in step 4, a new energy park model is constructed through constraints (A-1) to (A-16), (A-25), (A-27), and (A-33), and the objective function (A-17) is solved. By calculating the power generation power of the power generation equipment, the hydrogen production power of the electrolyzer, and the charging and discharging power of the energy storage equipment, the minimum total energy purchase cost of all parks is finally solved; a hydrogen tube trailer model is constructed through constraints (A-19) to (A-33), and the objective function (A-18) is solved. By calculating the route and specific time of hydrogen transportation, the minimum hydrogen tube trailer transportation cost is obtained.

[0072] Case Analysis

[0073] This paper uses a 5-node park for example analysis. Figure 2 As shown. Node Iron is a photovoltaic steel park, equipped with a 120MW photovoltaic generator set and eight 10MW electrolyzers. Node HRSS1 is also a photovoltaic power generation park, with a built-in 40MW photovoltaic generator set, a 6MW traditional gas turbine and a 10MW electrolyzer. Node HRSS2 is a wind power generation park, equipped with a 40MW wind turbine and a 10MW electrolyzer. In addition, two hydrogen refueling stations HRS are added to the system to bear the hydrogen energy load. The present invention is implemented through the GAMS optimization platform, and the GUROBI solver is used to solve the MIQCP problem.

[0074] Based on this example, the present invention first analyzed the differences in cost and carbon emissions between hydrogen ironmaking and traditional coke ironmaking (see Table 1 for results). The analysis results show that the process of replacing carbon with hydrogen can effectively reduce carbon emissions. However, the hydrogen production capacity of a single steel park is limited, and the need to purchase hydrogen will significantly increase costs. Therefore, in the multi-park integrated energy system, the ironmaking cost of the park is reduced by optimizing the hydrogen energy configuration, and the changes in the total energy purchase cost before and after the park interconnection are further compared (see Table 2 for results). The results show that park interconnection can not only ensure the supply of hydrogen energy, but also enhance the system's ability to absorb wind and solar energy, thereby improving the system's operating economy and energy utilization efficiency.

[0075] Table 1 Cost and carbon emissions of different ironmaking methods

[0076] plan Coke ironmaking Purchasing hydrogen for ironmaking Electric hydrogen production and ironmaking Hydrogen ironmaking after park interconnection Ironmaking cost (yuan) 102375.00 365622.08 174633.40 76589.025 Carbon emissions (t) 245.70 0 20.849 0

[0077] Table 2 Optimization results before and after campus interconnection

[0078] plan Before the park interconnection After the park is interconnected Total energy purchase cost of the park (yuan) 271107.27 141275.05

[0079] This invention is based on a detailed model of a new energy steel park. It couples the electric-hydrogen energy network through electric hydrogen production, while taking into account the insufficient hydrogen supply in the steel park and the power consumption problems caused by wind and solar fluctuations in the new energy park. It improves energy utilization and enhances the economic efficiency of system operation through the hydrogen energy interconnection method of multi-park integrated energy systems.

[0080] The above description is only a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any changes or substitutions that can be easily conceived by a technician familiar with this technical field within the technical scope disclosed by the present invention should be covered by the scope of protection of the present invention.

Claims

1. A multi-park integrated energy system operation method considering hydrogen energy optimization scheduling, characterized in that: The method comprises the following steps: Step 1: Obtain the operating parameters of the equipment in each park. For the new energy photovoltaic steel park, this includes the installed capacity of photovoltaic generators, electrolyzer capacity, fuel cell efficiency, energy storage equipment capacity, and charge and discharge loss efficiency; for the new energy park, this includes the installed capacity of wind turbines, photovoltaic power generation capacity, and equipment parameters of electrolyzers; for hydrogen refueling stations, this includes the hydrogen demand for each time period and the capacity of hydrogen storage tanks; Step 2: Obtain load demand data for each park, including power load, heat load, and steel demand, and collect geographical distances between parks and wind power output information; Step 3: Based on the above information, construct each park with energy supply and demand balance constraints and equipment operation constraints as constraints; Step 4. Based on the park constraints and combined with the hydrogen demand of the hydrogen refueling station, a hydrogen energy optimization scheduling model for the multi-park integrated energy system is constructed. This model takes the overall hydrogen energy supply and demand balance, the capacity constraint of hydrogen tube trailers for hydrogen transportation, the path planning constraint, and the time constraint as the constraints. The model takes the minimization of the overall energy purchase cost of the multi-parks and the operating cost of the hydrogen tube trailers as the objective function. The model is solved using the GUROBI solver to obtain the operation plan of the equipment in each park, including the power generation power of the power generation equipment, the hydrogen production power of the electrolyzer, the charging and discharging power of the energy storage equipment, and the path and time arrangement of hydrogen transportation, to achieve the optimal scheduling of the multi-park integrated energy system; In step 4, the minimum objective function is: Where T represents the time section number, subscript D represents the hydrogen refueling station, subscript K represents the set of hydrogen tube trailers, and subscript k represents the specific hydrogen tube trailer vehicle number; C e t represents the time-of-use electricity price, C h Indicates the price of industrial hydrogen per kilogram, C g Indicates the price per kilogram of methane, L R,D , L D,ic , L R,ic A represents the distance between the new energy park R and the hydrogen station D, the hydrogen station D and the steel park ic, and the new energy park R and the steel park ic. R,D,k,t Indicates the running status of vehicle k at time t. If the vehicle starts from the new energy park R and goes to the hydrogen refueling station D, the value is 1, otherwise it is 0. D,ic,k,t 、A R,ic,k,t Indicates the running status of the vehicle from D and R to the steel park ic. If it is in operation, the value is 1, otherwise it is 0. car represents vehicle travel expenses; P ce pa,t 、M ch pa,t 、M cg pa,t They represent the additional electricity purchase, additional hydrogen purchase and additional gas purchase of the pa park at time t respectively.

2. A multi-park integrated energy system operation method considering hydrogen energy optimization scheduling according to claim 1, characterized in that: In step 3, the relevant operating constraints include energy supply and demand balance constraints, as follows: Where, subscript t represents the scheduling period, subscript pa represents all park nodes, subscript ic represents the steel park, subscript R represents the new energy hydrogen production park, subscript sD represents the energy storage element, superscript pv represents the photovoltaic generator set, superscript el represents the electrolyzer, superscript fc represents the fuel cell, and superscript gb represents the boiler; P ch e,ic,t With P ch h,ic,t They represent the electricity and hydrogen reserves of the ic park at time t, P ch g,pa,t P represents the power of methane storage in the park at time t, dis e,ic,t With P dis h,ic,t They represent the capacity of electricity and hydrogen released by the IC park at time t, P dis g,pa,t The power of methane gas released by the pa park at time t, P pv ic,t 、P ce ic,t With P fc ic,t They represent the photovoltaic power generation of IC park, additional power purchase and hydrogen fuel cell power generation at time t, respectively. ch ic,t Indicates the hydrogen production state of the electrolyzer in the ic park at time t. If it is in operation state, the value is 1, otherwise it is 0. el ic,t Indicates the power consumption of the electrolytic cell in the ic park at time t, P l ic,t represents the power load of the ic park at time t, P pvmin With P pvmax Respectively represent the maximum and minimum values ​​of photovoltaic output, M el ic,t and M el R,t Respectively represent the hydrogen production of ic and R park electrolyzer at time t, M fc ic,t 、M car t With M dc t They represent the hydrogen consumption of hydrogen fuel cells, new energy hydrogen vehicles for steel transportation, and hydrogen consumption for hydrogen ironmaking, respectively. ch ic,t represents the amount of hydrogen purchased by IC Park at time t, Q gb ic,t , Q fc ic,t With Q l ic,t They represent the heat output and heat load demand of the IC park boiler and fuel cell at time t, M cg pa,t With M gb pa,t They represent the gas purchase volume of the park and the methane demand of the boiler at time t, P ch h,R,t 、P dis h,R,t They represent the amount of hydrogen stored and released in R Park at time t, M give R,t represents the excess hydrogen amount in R Park at time t, P ce pa,t 、M ch pa,t 、M cg pa,t They represent the additional electricity purchase, additional hydrogen purchase and additional gas purchase of the park at time t, P cemax 、M chmax With M cgmax They represent the upper limits of electricity, hydrogen and gas purchases from the power grid respectively.

3. A multi-park integrated energy system operation method considering hydrogen energy optimization scheduling according to claim 2, characterized in that: In step 3, the relevant operating constraints include the electrolytic cell operating constraints, which are as follows: Where, subscript v represents all parks in the interconnected system, M el pa,t represents the hydrogen production of the electrolyzer in the pa park at time t, η el represents the electrolysis efficiency, X ch pa,t Indicates the hydrogen production status of the electrolyzer in the pa park at time t. The value is 1 if the operation status is running, otherwise it is 0. el pa,t P represents the power consumption of the electrolytic cell in the park at time t, elmax Indicates the maximum allowable input power of the electrolyzer.

4. A multi-park integrated energy system operation method considering hydrogen energy optimization scheduling according to claim 3, characterized in that: In step 3, the relevant operating constraints include energy storage equipment operating constraints, as follows: Where S sD Indicates the rated capacity of the energy storage device, U ch pa,t and U dis pa,t They represent whether the park pa stores or releases energy at time t. A value of 1 indicates that energy storage or release begins, while a value of 0 indicates that energy is not stored or released. ch sD,pa,t With P dis sD,pa,t They represent the energy stored and released in the park at time t, P ss sD,pa,t represents the actual capacity of the energy storage equipment in the park at time t, P ssmin sD With P ssmax sD Respectively represent the actual maximum and minimum capacity of the energy storage device, λ ch sD and λ dis sD They represent the efficiency of energy storage and release respectively.

5. A multi-park integrated energy system operation method considering hydrogen energy optimization scheduling according to claim 4, characterized in that: In step 4, the hydrogen delivery constraints of hydrogen tube trailers, vehicle routing constraints and time constraints, and multi-park hydrogen supply and demand balance constraints are: (1) Hydrogen long tube trailer route constraints (2) Constraints on hydrogen transport by hydrogen tube trailers (3) Hydrogen delivery time coupling constraints (4) Conservation constraints on hydrogen supply and demand Where A i,j,k,t Indicates the running status of the hydrogen vehicle. If vehicle k moves from park i to park j at time t, the value is 1, otherwise it is 0. max i H represents the maximum number of hydrogen delivery vehicles owned by park i. need D,k,t represents the amount of hydrogen supplied to hydrogen refueling station D by vehicle k at time t, M h i,k,t M represents the amount of hydrogen carried by the kth vehicle when leaving the i-th park. h R,k,t 、M h D,k,t M represents the amount of hydrogen carried by the kth vehicle leaving the new energy park R and the hydrogen refueling station D at time t, respectively. hmax Indicates the maximum amount of hydrogen carried by the hydrogen tube trailer, S carH k,t H represents the amount of hydrogen supplied to the park by the kth vehicle at time t. give R,k,t Indicates the amount of hydrogen actually delivered to the steel park at time t, represents t-Δt give R,k,t -Δt R,D,k,t -Δt D,ic,k,t -Δt R,ic,k,t The amount of hydrogen supplied to the steel park by the kth vehicle at time Δt i,j,k,t represents the time it takes for the kth vehicle to reach park j from park i, Δt R,D,k,t , Δt D,ic,k,t , Δt R,ic,k,t They represent the time consumed by the kth vehicle from the new energy park R to the hydrogen station D, the hydrogen station D to the steel park ic, and the new energy park R to the steel park ic at time t, respectively. leave i,k,t and t arrive j,k,t They represent the time when the kth vehicle leaves park i and arrives at park j at time t, Δt give R,k,t Indicates the time consumed from R to the steel park, ceil indicates the use of rounded up counting, t arrive i,k,t represents the time when the kth vehicle arrives at park i at time t, t arrivemax t , t arrivemin t , t leavemax t and t leavemin t Indicates the earliest or latest time of arrival or departure from the park within the time period t, t leave R,k,t represents the specific time when vehicle k leaves R park at time t, t arrive R,k,t-1 represents the specific time when vehicle k arrives at park R at time t-1, Indicates the hydrogen flow rate, which is assumed to be a constant value here.

6. A multi-park integrated energy system operation method considering hydrogen energy optimization scheduling according to claim 5, characterized in that: In step 4, a new energy park model is constructed through constraints (A-1) to (A-16), (A-25), (A-27), and (A-33), and the objective function (A-17) is solved. By calculating the power generation power of the power generation equipment, the hydrogen production power of the electrolyzer, and the charging and discharging power of the energy storage equipment, the minimum total energy purchase cost of all parks is finally solved; a hydrogen tube trailer model is constructed through constraints (A-19) to (A-33), and the objective function (A-18) is solved. By calculating the route and specific time of hydrogen transportation, the minimum hydrogen tube trailer transportation cost is obtained.

Images