Next Article in Journal
A New Multi-Objective Optimization Strategy for Improved C3MR Liquefaction Process
Next Article in Special Issue
Special Issue: Advances in Renewable Energy Systems
Previous Article in Journal
Robust Observer-Based Proportional Derivative Fuzzy Control Approach for Discrete-Time Nonlinear Descriptor Systems with Transient Response Requirements
Previous Article in Special Issue
Simulation Experiment Design and Control Strategy Analysis in Teaching of Hydrogen-Electric Coupling System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Multi-Objective Energy Management of Renewable Energy Power Plant with Electrolytic Hydrogen Production

1
College of Automation, Nanjing University of Posts and Telecommunications, Nanjing 210023, China
2
Institute of Advanced Technology for Carbon Neutrality, Nanjing University of Posts and Telecommunications, Nanjing 210023, China
*
Author to whom correspondence should be addressed.
Processes 2024, 12(3), 541; https://doi.org/10.3390/pr12030541
Submission received: 26 February 2024 / Revised: 7 March 2024 / Accepted: 7 March 2024 / Published: 9 March 2024
(This article belongs to the Special Issue Advances in Renewable Energy Systems)

Abstract

:
This study focuses on a renewable energy power plant equipped with electrolytic hydrogen production system, aiming to optimize energy management to smooth renewable energy generation fluctuations, participate in peak shaving auxiliary services, and increase the absorption space for renewable energy. A multi-objective energy management model and corresponding algorithms were developed, incorporating considerations of cost, pricing, and the operational constraints of a renewable energy generating unit and electrolytic hydrogen production system. By introducing uncertain programming, the uncertainty issues associated with renewable energy output were successfully addressed and an improved particle swarm optimization algorithm was employed for solving. A simulation system established on the Matlab platform verified the effectiveness of the model and algorithms, demonstrating that this approach can effectively meet the demands of the electricity market while enhancing the utilization rate of renewable energies.

1. Introduction

China’s “3060 Dual Carbon” target, announced in 2020, has led to a consistent rise in the installed capacity of renewable energy generating units, specifically wind power and photovoltaic systems. As of the end of 2022, China’s wind power and photovoltaic capacity reached 7576.1 MW, representing 26.7% of the total installed power capacity. However, the variability of renewable energy output, including intermittency, randomness, and fluctuation, has posed challenges to the stable operation and reliability of the power grid. The lack of dynamic regulation ability of the power system and the issue of growing challenges in renewable energy consumption are becoming more prominent [1,2]. Electrolytic hydrogen production technology, as a flexible resource that can be adjusted, has advantages such as fast response speed and zero carbon emissions [3,4]; moreover, compared with traditional energy storage media, hydrogen energy offers benefits such as high energy density and extended storage duration [5]. The grid-connected operation of the electrolytic hydrogen production system (EHPS) not only helps to improve the fluctuation of wind and solar power output and reduce the rate of wind and solar abandonment (WSA) caused by insufficient power grid consumption capacity but also can participate in peak shaving auxiliary services (PSAS), enhance the dynamic regulation capability of the system, and ensure the secure, stable, and cost-effective operation of the power grid with a high penetration of renewable energy. Therefore, how to effectively use electrolytic hydrogen production technology, enhance the flexibility and regulation capabilities of renewable energy power plants, meet the demand of multiple target scenarios such as smoothing power fluctuation, increasing renewable energy consumption space, and participate in PSAS have emerged as a pressing issue in need of urgent solutions.
There has been a certain research foundation at home and abroad in the coordinated operation of electrolytic hydrogen production systems and renewable energy generating units. The literature [6] investigates a hybrid renewable energy system incorporating hydrogen energy, proposing an efficient hydrogen production system energy management. It considers two different hydrogen production methods and establishes a mathematical model for energy conversion between hydrogen and electricity. The literature [7] explores potential opportunities and risks associated with hydrogen storage systems, assesses the feasibility of using hydrogen storage technology to address the issue of wind and solar energy curtailment in Western China, and proposes two solutions based on fuel cell hydrogen storage systems and hydrogen–natural gas blending. The literature [8] investigates an isolated DC microgrid with a combined electricity–hydrogen storage system. It proposes a cost-minimizing hybrid storage charging and discharging strategy that ensures voltage and power balance within the microgrid. The literature [9] aims to maximize system benefits and reduce environmental costs, constructs a wind–solar complementary hydrogen production system scheduling model, and uses a multi-objective golden eagle algorithm to solve the model. The literature [10] explores refined modeling of hydrogen storage systems and power-to-gas equipment, effectively enhancing the economic performance of integrated energy systems incorporating hydrogen storage. The literature [11] addresses the optimization and scheduling of hybrid systems with hydrogen storage by considering detailed models of flexible loads and hydrogen storage, aiming for an efficient solution. The objective of [12] is to optimize the net income of the microgrid within the service area, with the power balance of the microgrid system as a constraint, construct an optimized dispatching model of the wind and solar power generation service area microgrid including battery and hydrogen energy storage units, and the results show the effectiveness of the optimized dispatching strategy. The literature [13] delves into the capacity configuration of photovoltaic systems combined with hybrid battery–hydrogen storage systems, catering to both electrical and hydrogen loads. It introduces an optimization framework for the energy system servicing the hydrogen load and proposes two energy management strategies to ensure the stable supply of both electricity and hydrogen to the loads. The study in [14] presents an optimized configuration scheme for a wind–solar–hydrogen storage integrated energy system, considering both economic factors and system output fluctuations. The literature [15] studies the control mode of the hydrogen production system, proposes an optimization model of the integrated energy system with hydrogen production, and verifies the role of hydrogen production in system economy and renewable energy consumption in simulation software. The literature [16] based on the wind–hydrogen coupling system develops a capacity optimization model to maximize the revenue of the system and considers the possibility of the system participating in the auxiliary service market to improve system economy. The literature [17] studies integration of the hydrogen production system into the auxiliary peak shaving service of the thermal power unit and proposes a strategy for its participation. The literature [18] investigates integrated electric–gas energy systems incorporating green hydrogen, considering their participation in ancillary services and support for high-proportion renewable energy integration.
The aforementioned research mainly focuses on the operation optimization under a single scenario. However, a single objective is insufficient to meet the demands of the power system; hence, multi-objective optimization scheduling is necessary. The literature [19] proposes a day-ahead optimization scheduling strategy for wind–solar–hydrogen systems, targeting the system’s curtailment rate of wind and solar, power deficit rate, and economic efficiency. The feasibility of its day-ahead optimization scheduling strategy was verified through simulation software. The literature [20] crafted an adaptable optimization algorithm focused on reducing expenses and enhancing battery storage capacity. Through the creation and testing of a power flow optimization framework, it facilitated increased adoption of photovoltaic power in standalone microgrids, alleviating the difficulty of managing operational expenses alongside microgrid dependability. The literature [21] proposes an optimized operation focusing on Battery Energy Storage Systems (BESS) and Power-to-Gas (P2G) technologies, aiming to maximize the benefits of Renewable Energy Sources, P2G, and BESS while minimizing facility costs. Optimization is conducted through Mixed-Integer Linear Programming (MILP), and Stochastic Programming (SP) is employed to manage the uncertainty of renewable energy output.
The aforementioned literature discusses multi-objective energy management in microgrid applications and does not focus on renewable energy power plants as the research objective. The multiple objectives in renewable energy power plants mentioned in this paper and the multiple objectives in a microgrid are not the same, but such objectives are necessary for the electricity market and the power system and demonstrate potential in energy management for renewable energy power plants with electrolyte hydrogen units, such as smoothing renewable energy output fluctuations, participating in PSAS. Meanwhile, the random changes in the output of a renewable energy-generating unit make it a random optimization problem under uncertain conditions, requiring special quantification of the randomness of renewable energy output. This paper quantifies the uncertainty of renewable energy output in the form of fuzzy chance constraints and adopts the method of crisp equivalent classes to convert fuzzy chance constraints into deterministic constraints during the solution process.
This paper considers the use of an electrolytic hydrogen production system to smooth the power fluctuations of renewable energy power generation, participate in PSAS, increase the space for renewable energy consumption, reduce the amount of wind and solar discarded, and other application target scenarios, while considering the uncertainty in renewable energy output. It proposes a multi-objective energy management model based on the coordinated optimization of an electrolytic hydrogen production system and renewable energy power generation, and a particle swarm algorithm based on piecewise mapping and improved Levy flight algorithm has been introduced. Finally, the proposed algorithm and model’s effectiveness is validated through case simulations. As a whole, the main contributions of this paper are as follows:
  • Multi-objective energy management: this paper presents a multi-objective energy management strategy that utilizes electrolysis hydrogen production technology to mitigate the adverse effects caused by renewable energy generation. This strategy enhances the stability of the power system, improves energy utilization, and offers a viable solution to address the challenges of renewable energy generation.
  • Economic and operational efficiency: this paper’s simulation findings deliver an in-depth assessment of the economic and operational advantages associated with a renewable energy power plant including an electrolytic hydrogen production system. It highlights the potential for considerable progress in energy management for renewable energy generation through the optimization of energy usage, involvement in peak-load adjustment services, and participation in carbon trading activities.
  • Uncertain programming: this paper quantifies the uncertainty of renewable energy output in the form of fuzzy chance constraints and adopts the method of crisp equivalent classes to convert fuzzy chance constraints into deterministic constraints during the solution process. The issue of unpredictability in renewable energy generation’s output has been addressed, offering a dependable framework for the enhancement of renewable energy systems.
  • Improved particle swarm algorithm: in this paper, the adoption of the particle swarm algorithm based on piecewise mapping and improved Levy flight algorithm is introduced for achieving comprehensive and effective optimization of the proposed model. This technique contributes to the improvement of the precision and dependability of the energy management approach.

2. Basic Principles of Hydrogen Production from Renewable Energy

The renewable energy power-generating system with an electrolytic hydrogen production system mainly uses the EHPS and the surplus power in the grid-connected consumption process of wind power and photovoltaic power generation to store and utilize in the form of hydrogen energy and chemical energy. Its typical structure is shown in Figure 1.
This plant consists of wind power generation, photovoltaic power generation, EHPS, and related power electronic interface equipment. Normally, wind power and photovoltaic power are directly consumed by the grid. However, the intermittent, random, and fluctuating nature of wind power and photovoltaic power generation, influenced by the natural environment, poses challenges, when the flexible regulation ability of the power grid is insufficient, it is necessary to use the EHPS to smooth the output power fluctuations of the renewable energy-generating unit, store the surplus output power of the renewable energy-generating unit, and participate in PSAS under certain price incentives.
The EHPS consists of an electrolyzer, a hydrogen compressor, and a hydrogen storage tank. The EHPS is an important device for electricity conversion and energy storage. It receives direct current from the renewable energy-generating unit and converts electrical energy into chemical energy through the electrolysis process. This conversion process allows the electrolyzer to store electrical energy when the demand of the power system is lower than the supply and to release the stored chemical energy for power supply when the demand is higher than the supply.

3. Technical and Economic Characteristics of Hydrogen Production from Renewable Energy

3.1. Characteristics of Renewable Energy-Generating Unit

3.1.1. Characteristics of Wind Power Generation

The output power P wt of the wind turbine is shown in Equation (1):
P w t = 1 2 C p ( λ , β ) ρ π R 2 v 3
In the equation, C P is the wind energy utilization coefficient; λ is the tip speed ratio; β is the pitch angle; ρ is the air density; R is the radius of the wind turbine rotor; and v is the wind speed.
The output power of the fan can be calculated according to Formula (1). This relationship can be described by the operating parameters of the fan. Commonly used parameters include cut-in wind speed, rated wind speed, cut-out wind speed, and rated power.
{ 0 0 v < v i n n ( A + B v + C v 3 ) P w t N v i n v v N n P w t N v N v v o u t 0 v o u t < v
where:
A = 1 v i n v N [ v i n ( v i n + v N ) 4 ( v i n × v N ) ( v i n + v N 2 v N ) 2 ]
B = 1 ( v i n v N ) 2 [ 4 ( v i n + v N ) ( v i n + v N 2 v N ) 3 ( 3 v i n + v N ) ]
C = 1 ( v i n v N ) 2 [ 2 4 ( v i n + v N 2 v N ) 3 ]
When approximated by linearization curve, the piecewise function expression of fan output power characteristics can be written as:
{ 0 0 v < v i n n P w t N ( v 3 v i n 3 ) / ( v N 3 v i n 3 ) v i n v v N n P w t N v N v v o u t 0 v o u t < v
In the equation, v i n and v o u t are the cut-in and cut-out wind speeds of the fan, respectively; n is the number of wind turbines; and P w t N is the rated power of the wind turbine.

3.1.2. Characteristics of Photovoltaic Power Generation

The power output of photovoltaic can be calculated using an engineering approximation formula, as shown below:
P p v = P S T C G A C [ 1 + k ( T c T r ) ] / G S T C
In the equation, P p v is the output power of photovoltaic power generation; G A C is the irradiance; P S T C is the maximum output power under standard test conditions; G S T C is the solar intensity under standard test conditions; k is the power temperature coefficient; T c is the working temperature of the photovoltaic module; and T r is the reference temperature.

3.2. Characteristics of Electrolytic Hydrogen Production and Hydrogen Usage

3.2.1. Electrolytic

The electrical model of electrolytic water hydrogen production can be represented as follows:
U e l e = U r + r 1 + r 2 T e l S I e l + K e l ln ( k t 1 + k t 2 T e l + k t 3 T e l 2 S I e l + 1 )
In the equation, U e l e is the working voltage of the electrolyzer; U r is the reversible voltage of the electrolyzer that varies with temperature; I e l is the working current of the electrolyzer; r 1 , r 2 is the ohmic resistance parameter of the electrolyte; T e l is the working temperature of the electrolyzer; K e l , k t 1 , k t 2 , and k t 3 are the voltage parameters on the electrode of the electrolyzer; and S is the surface area of the electrolyzer motor.
The electrolyzer’s hydrogen production rate model is described by the following equation:
V e l e c h ( t ) = η e l e P e l e ( t )
In the equation, V e l e c h ( t ) is the amount of hydrogen output; η e l e is the hydrogen production efficiency of the electrolyzer; and P e l e is the output power of the electrolyzer.

3.2.2. Hydrogen Storage Tank

According to the different forms of hydrogen, hydrogen storage technology mainly includes solid, gaseous, and liquid hydrogen storage. Industrial hydrogen storage uses gaseous hydrogen storage. The gas pressure P H T ( t ) in the hydrogen storage tank at time t is related to the amount of hydrogen storage n H T ( t ) , calculated using the ideal gas equation.
P H T ( t ) = n H T ( t ) R T H 2 V H T n H T ( t ) = n H T ( t 1 ) + n E L ( t ) n F C ( t ) Δ t
In the equation, T H 2 is the temperature of the hydrogen gas; V H T is the volume of the hydrogen storage tank; and n H T ( t 1 ) is the amount of hydrogen gas in the tank at time t − 1. The hydrogen storage status of the hydrogen storage tank can be calculated as follows:
S O H T ( t ) = P H T ( t ) P N
In the equation, PN is the maximum pressure of the hydrogen storage tank.

3.2.3. Fuel Cell

A fuel cell is similar to the reverse process of hydrogen production in an electrolytic cell. The output voltage can be represented as follows:
u c e l l = u n e r n s t u a c t u o h m u c o n
In the formula, u c e l l is the output voltage; u n e r n s t is the thermodynamic electromotive force; uact is the activation polarization overpotential; uohm is the ohmic polarization overpotential; and ucon is the concentration difference overpotential.

4. Optimization Model of Multi-Objective Energy Management

4.1. Optimization Objectives of Energy Management

4.1.1. Smoothing of Power Fluctuations

The integration of renewable energy with significant fluctuations into the grid can cause harm to the grid; thus, certain requirements need to be proposed for the grid-connected renewable energy:
| P n e t ( t ) P n e t ( t 1 ) | D e t
P n e t ( t ) = P wt ( t ) + P p v ( t )
In the formula, P n e t represents the actual output of renewable energy; Det represents the fluctuation range that the power grid can withstand, calculated as 10% of the total installed capacity of the plant’s renewable energy generating units; P wt represents the output power of wind power generation; and P p v represents the output power of photovoltaic power generation.

4.1.2. Participation in Peak Load Regulation Auxiliary Services

Energy storage facilities absorb power during off-peak periods or during periods of wind, solar, and nuclear power curtailment and release power at other times, thereby providing peak load regulation auxiliary services.
According to the “Compensation Management Method” in the “Operating Rules of the Northeast China Power Auxiliary Service Market” released by the Northeast China Power Auxiliary Service Market, a “step-by-step” compensation strategy is proposed. Peak load regulation pricing details are shown in Table 1. The energy storage facilities built within the metering export of the thermal power plant, in conjunction with the unit, participate in peak shaving and are managed, cost-calculated, and compensated according to the “Compensation Management Method”.
Due to the current rules and regulations regarding energy storage participation in peak shaving not being fully developed or comprehensive, the energy storage facilities built within the metering export of renewable energy power stations participate in PSAS. The “Compensation Management Method”, in which the energy storage facilities are established within the metering export of thermal power plants, is used as a reference to establish a compensation strategy:
E sup = λ sup Q sup
λ sup = { 1000 0 < f i 40 % 400 40 % < f i 50 % 0 50 % f i
In the formula, E sup represents the compensation for participating in peak load regulation auxiliary services; λ sup represents the subsidized electricity price; Q sup represents the amount of electricity for output peak shaving; and f i represents the average load rate of the energy storage facilities. The average load rate is calculated as the power generated by the operating unit divided by the capacity of the operating unit, multiplied by 100%.

4.1.3. Increasing the Absorption Space for Renewable Energy

For renewable energy power plants, there are three factors that cause the system to abandon wind and solar. First, the installed capacity of the electrolytic cell is limited. Second, the capacity of the hydrogen storage tank is limited. Third, the installed capacity of electrochemical energy storage is limited.
The power quantity of the plant’s abandoned wind and solar is:
Q loss ( t ) = [ P n e t , t h e o r y ( t ) P n e t ( t ) P e l e ( t ) | V H < V H max P b a t c h ( t ) ] Δ t
In the formula, Q loss represents the amount of electricity discarded from wind and solar; P n e t , t h e o r y represents the theoretical output of renewable energy; P e l e n e t ( t ) | V H < V H max represents the operating power of the electrolyzer when absorbing the surplus output of renewable energy and when the storage volume of the hydrogen storage tank is less than the maximum storage volume; P b a t c h represents the charging power of electrochemical energy storage; Δ t represents the time interval; and Δ t equals 15 min.

4.2. Objective Function

Total Benefit Model

Max   f = E sup + E co 2 K loss K buy
In the formula, E sup represents the revenue from participating in peak load shifting auxiliary services; E co 2 represents the revenue from carbon emission reduction; K loss represents the cost of wind and solar power curtailment; and K buy represents the cost of purchasing electricity.
1.
Revenue of peak load shifting auxiliary services:
E sup = λ sup P sup ( t ) Δ t
2.
Revenue of carbon emission reduction:
E c o 2 = A c o 2 ( B q B p )
B p = λ p t = 1 T P buy ( t ) Δ t
B q = λ q t = 1 T ( P w t + P p v P loss ) Δ t
In the formula, A c o 2 represents the carbon trading price; B p represents the plant carbon emissions; B q represents the allocated carbon emission quota; λ p represents the equivalent carbon emission coefficient of purchasing electricity from the grid [22]; T represents the optimization period (one day is divided into 96 time periods, that is, T = 96); P b u y represents the power of purchasing electricity; and λ q represents the carbon emission quota per unit of electricity generation.
3.
Cost of purchasing electricity:
K buy = t = 1 T A buy P buy ( t ) Δ t
In the formula, A buy represents the time-of-use electricity price and P b u y represents the power of purchasing electricity.
4.
Cost of wind and solar power curtailment:
K loss = γ loss t = 1 T Q loss ( t )
Q loss ( t ) = [ P n e t , t h e o r y ( t ) P n e t ( t ) P e l e ( t ) | V H < V H max P bat c h ( t ) ] Δ t
In the formula, γ l o s s represents the penalty for the unit amount of wind and solar power curtailment.

4.3. Constraint Conditions

4.3.1. Overall Operational Constraints of Renewable Energy Power Stations

  • Power and energy balance constraints:
P wt ( t ) + P p v ( t ) = P s e l l ( t ) + P e l e ( t ) P b u y ( t ) P sup ( t ) P f c ( t ) + P b a t c h ( t ) P b a t d i s c h ( t ) + P l o s s ( t )
In the formula, P s e l l represents the power of selling electricity; P e l e represents the operating power of the electrolytic cell; P sup represents the power involved in peak shifting; and P b a t c h and P b a t d i s c h represent the charging power and discharging power, respectively.
2.
Hydrogen energy supply and demand balance constraints:
V h s ( t + 1 ) = V h s ( t ) + V e l e c h ( t ) V f c d i s c h ( t )
V e l e c h ( t ) = η e l e P e l e ( t ) Δ t
V f c d i s c h ( t ) = η f c P fc ( t ) Δ t
In the formula, V h s represents the storage of hydrogen in the hydrogen storage tank; V f c d i s c h represents the amount of hydrogen used by the fuel cell per unit time; η f c represents the hydrogen utilization efficiency of the fuel cell; and P f c represents the operating power of the fuel cell.
3.
Power fluctuation constraints:
| P n e t ( t ) P n e t ( t 1 ) | D e t

4.3.2. Operational Constraints of Each Subsystem

  • Operational constraints of renewable energy-generating unit:
P w t min ( t ) P w t ( t ) P w t max ( t )
P p v min ( t ) P p v ( t ) P p v max ( t )
In the formula, P w t min and P w t max represent the minimum and maximum output power of wind power, respectively; and P p v min and P p v max represent the minimum and maximum output power of photovoltaic power, respectively.
2.
Operational Constraints of Hydrogen Energy Systems
  • Operational constraints of electrolytic cell:
P e l e min P e l e ( t ) P e l e max
In the formula, P e l e min represents the lower limit of the power consumed by the electrolytic cell during normal operation, and there is a risk of explosion when the operating power of the electrolytic cell is below 20~25% of the rated power for a long time [23]; P e l e max represents the upper limit of the power consumed by the electrolytic cell during normal operation.
   •
Operational constraints of fuel cells:
0 P f c ( t ) P f c max
In the formula, P f c max represents the upper limit of the power consumed by the fuel cell during normal operation.
   •
Capacity constraints of hydrogen storage tanks:
0 V h s ( t ) V h s max
In the formula, V h s max represents the maximum hydrogen storage capacity of the hydrogen storage tank.
3.
Operational Constraints of Electrochemical Energy Storage Systems
  • Charging and discharging power constraints:
{ 0 P b a t c h ( t ) P b a t r a t e d P b a t r a t e d P b a t d i s c h ( t ) 0 P b a t c h ( t ) P b a t d i s c h ( t ) = 0
In the formula, P b a t r a t e d represents the rated power of the electrochemical energy storage.
   •
Energy storage state constraints:
Considering that overcharging and over-discharging of electrochemical energy storage can significantly reduce the cycle life of the battery, it is significant to impose constraints on the state of charge (SOC) of the energy storage system:
S O C ( t ) = E ( t ) E b a t r a t e d
S O C ( t ) S O C ( t 1 ) = ( P b a t c h ( t ) η b a t c h E b a t r a t e d P b a t d i s c h ( t ) E b a t r a t e d η b a t d i s c h ) Δ t
S O C min S O C S O C max
In the formula, E ( t ) and E b a t r a t e d represent the current energy state and rated energy state of the electrochemical energy storage, respectively; η b a t c h and η b a t d i s c h represent the charging efficiency and discharging efficiency, respectively; S O C min and S O C max represent the maximum capacity and minimum capacity of the electrochemical energy storage, respectively.

4.4. Fuzzy Chance Constraints

4.4.1. Fundamental Principle

The modeling idea of fuzzy chance-constrained programming is it allows the decision made to not satisfy the constraint conditions to some extent, but the probability of the fuzzy constraint conditions being established is not less than the preset confidence level.
The single-objective programming with fuzzy chance constraints is as follows:
{ min f ( x , ζ ) s . t . Pr { g ( x , ζ ) 0 } α
In the formula, x represents the decision variable; ξ represents the fuzzy vector; min f ( x , ζ ) represents the objective function; g ( x , ζ ) 0 represents the constraint conditions; Pr { } represents the probability of a certain event occurring; and α represents the confidence level of the system.

4.4.2. Handling of Power and Energy Balance Constraints

Different from traditional optimization methods, due to the randomness of renewable energy output, when a significant amount of renewable energy is integrated into the grid, the error between actual output and predicted output cannot be ignored. In the energy management model, the output of renewable energy has uncertainty. A model is established based on the theory of Chance Constrained Programming, and the fuzzy parameters P ˜ w t and P ˜ p v of renewable energy are introduced to form fuzzy chance constraints. Therefore, Equation (26) can be relaxed to the power and energy balance constraints under a certain confidence level condition, ensuring that the probability of meeting the balance constraints is not less than α .
Pr { P ˜ wt ( t ) + P ˜ pv ( t ) = P s e l l ( t ) + P e l e ( t ) P buy ( t ) P sup ( t ) P f c ( t ) + P b a t c h ( t ) P b a t d i s c h ( t ) + P loss ( t ) } α

5. Solution Method of Energy Management Model

5.1. Handling Method of Fuzzy Chance Constraints

The methods for handling fuzzy chance-constrained programming are divided into two categories. For simple problems, the fuzzy parameters and decision variables in the constraint conditions can be separated or, when they have a certain linear relationship, they can be transformed into clear equivalent classes for processing and then calculated by traditional methods; for complex problems, random simulation methods can be used for processing. However, the results of random simulation are not accurate and it is not easy to grasp the size of the sample capacity of the simulation, so the first method is chosen [24].
During the operation of the power system in practice, the requirements for safety are high; therefore, the probability of the power and energy balance constraints being established should be above 0.5, that is, α 0.5 , so:
Assume the form of function g ( x , ζ ) is as follows [25]:
g ( x , ζ ) = ψ 1 ( x ) ζ 1 + ψ 2 ( x ) ζ 2 + + ψ t ( x ) ζ t + ψ 0 ( x )
In the formula, represents the trapezoidal fuzzy variable ( γ k 1 , γ k 2 , γ k 3 , γ k 4 ), ψ = 1,2,…, t. If ψ k + 0 and ψ k 0 , k = 1, 2, …, t. So:
The fuzzy chance constraints in Equation (42) can be clearly equivalent to:
( 2 2 α ) k = 1 t ( γ k 3 ψ k + ( x ) γ k 2 ψ k ( x ) ) + ( 2 α 1 ) k = 1 t ( γ k 4 ψ k + ( x ) γ k 1 ψ k ( x ) ) + ψ 0 ( x ) 0
This paper describes the fuzzy parameters P ˜ wt and P ˜ pv of the renewable energy output using trapezoidal fuzzy parameters. Among them, the trapezoidal fuzzy variable is a quadruple ( P w t 1 , P w t 2 , P w t 3 , P w t 4 ) composed of clear numbers and P w t 1 < P w t 2 P w t 3 < P w t 4 . The quadruple values can be determined based on the predicted renewable energy output:
{ P ˜ w t ( ω 1 , ω 2 , ω 3 , ω 4 ) = P w t f o r e ( r 1 , r 2 , r 3 , r 4 ) P ˜ pv ( ω 1 , ω 2 , ω 3 , ω 4 ) = P p v f o r e ( r 1 , r 2 , r 3 , r 4 )
In the formula, ω 1 , ω 2 , ω 3 , ω 4 represents the membership parameter; r 1 , r 2 , r 3 , r 4 represents the proportionality constant.
The wind and solar prediction output are parameterized in Equation (44):
{ P ˜ w t α 1 4 P w t 1 + α 2 P w t 2 + 3 α 4 P w t 3 + 1 α 2 P w t 4 P ˜ p v α 1 4 P p v 1 + α 2 P p v 2 + 3 α 4 P p v 3 + 1 α 2 P p v 4

5.2. Improved Particle Swarm Algorithm

5.2.1. Standard Particle Swarm Optimization (PSO) Algorithm

The update formula for the PSO algorithm is given by the following equation:
v i = ω × v i + c 1 × r a n d ( ) × ( P b e s t i x i ) + c 2 × r a n d ( ) × ( g b e s t i x i )
x i = x i + v i
In the formula, i = 1, 2, …, N, where N is the total number of particles in this group; ω is the inertia factor; v i is the speed of the particle; r a n d ( ) is a random number between (0, 1); x i is the current position of the particle; and c 1 and c 2 are learning factors.

5.2.2. Particle Swarm Algorithm Based on Piecewise Mapping and Improved Levy Flight (PLPSO)

(1)
Piecewise Mapping
Chaos mapping is a method for generating chaotic sequences. Piecewise mapping, as a typical chaos mapping, has a concise mathematical form and good ergodicity and randomness. The definition of piecewise chaotic mapping is as follows:
x ( t + 1 ) = { x ( t ) p , 0 x ( t ) < p x ( t ) p 0.5 p , p x ( t ) < 0.5 1 p x ( t ) 0.5 p , 0.5 x ( t ) < 1 p 1 x ( t ) p , 1 p x ( t ) < 1
In the formula, p = 1 ; x ( 1 ) = r a n d .
(2)
Improved Levi Flight
Levy flight is a method of simulating random motion using heavy-tailed distributions. It can describe the dynamic characteristics of some complex systems, such as phenomena in the fields of ecology, geography, physics, etc. The update formula for Levy flight is as follows:
X t + 1 = X t + a L e v y ( β )
In the formula, X t is the t-th iteration position of X ; a is the step length scaling factor; and L e v y ( β ) is the Levy random path, satisfying: L e v y ~ u = t β ,   1 β 3 .
Levy flight is essentially a random step length, and the distribution of the step length satisfies the Levy distribution. The Levy distribution is complex and is usually implemented using the Mantegna algorithm. For the specific implementation process, please refer to reference [26].
Due to the step length scaling factor of Levy flight usually being set to a fixed value, it leads to slow search speed and low search accuracy during the search process [27]. To address these issues, this paper introduces an enhanced Levy flight approach:
X t + 1 = X t + a sin 2 ( π 2 1 t i t e r ) L e v y ( β )
In the formula, the fixed step length a = 0.01 ; sin 2 ( π 2 1 t i t e r ) is the adaptive step length factor, where t is the current iteration number and iter is the maximum number of iterations.
(3)
Improved Particle Swarm Algorithm
PSO has some disadvantages, such as easy premature convergence, lack of diversity, and susceptibility to initial values. The improved PSO uses chaos mapping to enhance the diversity of particles and combines with improved Levy flight to enhance the vitality and jumping ability of particles [26], realizing the complementary advantages of PSO, chaos mapping, and Levy flight. Due to the complexity of Levy flight, the algorithm’s computation time will experience an increase. However, in the later stage of optimization, the particle swarm algorithm will fall into local optimum, and particles need to perform Levy flight to escape from the local optimum. Therefore, an adaptive random number ( 1 + t i t e r ) r a n d is introduced. Before each update of particles, the generated adaptive random number is compared with 0.5. If it is greater than 0.5, Levy flight is used to update the particles; otherwise, the particles’ velocity and position are updated iteratively through Formulas (46) and (47).
The process description of the improved PSO is as follows:
(a)
Population initialization. Randomly initialize the speed of the population and use chaos mapping to initialize the position of the population and find the current optimal particle position and solution.
(b)
Calculate fitness. Evaluate the fitness value of each particle and store it in pbest. Compare the fitness values of all particles in the population and store the best individual (gbest) in the population.
(c)
Update particle position. Generate an adaptive random number and compare it with 0.5. If it is greater than 0.5, update the particle through Levy flight; otherwise, update the particle through Formulas (46) and (47).
(d)
Update individual best position. Compare the current fitness value of the particle with the fitness value before the update, and use the better particle as the particle after the iteration.
(e)
Update global optimal solution. Compare pbest with gbest. If pbest is better than gbest, assign the value of pbest to gbest.
(f)
Algorithm termination. If the iteration termination condition is met, the algorithm terminates; otherwise, return to step c to continue the search.
The process of the PLPSO algorithm is shown in Figure 2:

6. Simulation Analysis

6.1. Case Description

Based on the typical daily wind power and photovoltaic output prediction data of a renewable energy power station in a province in Northeast China, with a sampling period of 15 min, the obtained forecast output curve of wind power and photovoltaic on a typical day is shown in Figure 3, and the program output and peak shaving command curve are shown in Figure 4. An EHPS and an electrochemical energy storage system are configured; the installed capacities of wind power and photovoltaics in this region are 25 MW and 200 MW, respectively, the capacity of the electrolytic cell is configured as 25 MW, and the fuel cell capacity is 10 MW; the equipment-related parameters are shown in Table 2 and the economic parameters are shown in Table 3.

6.2. Comparative Analysis of PSO Algorithm before and after Improvement

The mathematical model proposed in this paper is solved by the traditional PSO and PLPSO algorithms, respectively. The two algorithms have the same parameters: ω = 0.8 , c 1 = c 2 = 1.2 , the population size is 100, the maximum number of iterations is 15,000, and both algorithms run 100 times under the scenario with a confidence level of 0.9. Table 4 gives a comparison of the results of the traditional PSO and PLPSO algorithms.
As can be seen from Table 4, the PLPSO algorithm is superior to the PSO algorithm in terms of running time, average value, and optimal value. Through the calculation of the average value, it can be known that the fitness value of the PLPSO algorithm is 114,057.37 smaller than that of the traditional PSO algorithm, and the convergence accuracy is improved by 17.4%; the mean absolute deviation (MAD) measures the average distance between data points and their average value. A smaller average absolute deviation value indicates that the data are more stable. The MAD of the data obtained by the PLPSO algorithm is 71.9% smaller than that of the PSO, indicating that the PLPSO algorithm is more stable. The convergence curves of the fitness values of the two algorithms are shown in Figure 5.
As can be seen from Figure 5, due to the use of piecewise mapping by the PLPSO algorithm, the initial fitness value of the PLPSO algorithm exhibits improvement over the initial fitness value of the traditional PSO algorithm; the PLPSO algorithm terminates convergence around the 5000th generation and the traditional PSO algorithm terminates convergence around the 7500th generation. In comparison to the traditional PSO algorithm, the PLPSO algorithm exhibits an approximately 33.3% faster convergence speed. As evident from the zoomed-in figure, at the end of convergence, the convergence accuracy of the PLPSO algorithm is 14.9% higher than that of the traditional PSO algorithm. Therefore, the PLPSO algorithm is superior to the traditional PSO algorithm.

6.3. Analysis of the Impact of Varying Confidence Levels on Decision-Making Outcomes

The confidence level significantly affects the system’s economy and safety, and different confidence levels will lead to different scheduling results. As the confidence level increases, the system’s rate of WSA and comprehensive benefits are shown in Table 5.
As shown in Table 5, as the confidence level continues to increase, the system’s wind and light abandonment rate generally shows a downward trend; under the premise that the installed capacity of the EHPS and the electrochemical energy storage system is certain, when the confidence level is low, due to the insufficient installed capacity of the EHPS and the electrochemical energy storage system, the rate of WSA is larger. To mitigate the rate of WSA, a portion of its capacity is allocated to participate in the PSAS, absorbing excess electricity, thereby reducing the economic benefits brought by engaging in the PSAS. Since the plant still has a high rate of WSA, it causes a large WSA penalty. Therefore, when the confidence level is low, the comprehensive benefits are low; when the confidence level is high, due to the sufficient installed capacity of the EHPS and the electrochemical energy storage system, the plant can easily decrease the rate of WSA and has more capacity to engage in PSAS to obtain the highest peak-shaving subsidy. Therefore, as the confidence level increases, the comprehensive benefits show an increasing trend.
The operation of the EHPS and the electrochemical energy storage system is intricately linked to the confidence level. During the scheduling decision-making process, the scheduler must carefully weigh the plant demand, safety, and economy factors to determine the optimal confidence level.

6.4. Analysis of Energy Management Decision-Making Scheme

Taking the confidence level as an example, using the PLPSO algorithm for solving, we can obtain the scheduling situation of the renewable energy output, the EHPS, and the electrochemical energy storage system 24 h before the day. Among them, the peak-shaving subsidy and the peak-shaving command curve are depicted in Figure 6.
The renewable energy power plant containing electrolytic hydrogen is derived from the planned output reported to the upper-level power grid, taking into account the time-of-use electricity price, power fluctuation suppression, peak-shaving ladder subsidy, carbon-trading benefits, WSA penalties, and various system operation constraints. It ensures that the actual output fluctuates between 90% and 110% of the planned output. While completing the PSAS, the excess electricity is utilized for hydrogen production (charge) through the EHPS and the electrochemical energy storage system. In the case of insufficient renewable energy output, it discharges. The output curves of the EHPS and the electrochemical energy storage system are displayed in Figure 7 and Figure 8, respectively.
As evident from Figure 7, between 0:00~2:45, 11:00~15:45, and 17:00~19:45, the electrolyzer and fuel cell operate at a certain power to maintain a lower level of hydrogen storage, leaving enough hydrogen storage space for the subsequent periods (3:00~4:45, 16:00~16:45, and 20:00~21:45) to prepare for peak shaving. During the peak-shaving period, the maximum operating power of the electrolyzer is 21.94 MW, the minimum operating power is 14.84 MW, and a total of 68,590 N m 3 hydrogen is produced; between 5:00 and 8:45, the fuel cell maintains a low-power operation to maintain a higher storage amount of hydrogen in the hydrogen storage tank, preparing for the next period (9:00~10:45) for peak shaving. During the peak-shaving period, the electrolyzer maintains operation at 25% of the rated power (i.e., 6.25 MW) to ensure the electrolyzer operates safely, and the fuel cell is in a high-power working state, with a maximum discharge power of 10 MW, a minimum discharge power of 3.52 MW, a total discharge of 52.02 MW, and a hydrogen consumption of 16,866.8 N m 3 ; between 22:00 and 24:00, the electrolyzer operates at high power and the fuel cell operates at low power to ensure that there is a higher amount of hydrogen storage in the hydrogen storage tank at 24:00 to obtain higher hydrogen sales benefits. In one operating cycle, the electrolyzer produces a total of 244,691.03 N m 3 hydrogen, the fuel cell consumes a total of 218,985.38 N m 3 hydrogen, and the final hydrogen storage amount is 29,993.85 N m 3 , which is close to the maximum hydrogen storage amount and can obtain higher hydrogen sales benefits.
As shown in Figure 8, during the peak-shaving periods of 3:00~4:45, 16:00~16:45, and 20:00~21:45, the electrochemical energy storage is charged, with a total charge of 23.77 MW in the three peak-shaving periods; during the peak-shaving period of 9:00~10:45, the electrochemical energy storage discharges, with a total discharge of 6.34 MW during the period; the SOC of the electrochemical energy storage is maintained between 0.2 and 0.8, ensuring that the electrochemical energy storage works in a “shallow charge and shallow discharge” state, which helps prolong the service life of the electrochemical energy storage system; the SOC of the electrochemical energy storage is maintained between 0.45 and 0.8 for a total of 66 periods (i.e., 68.75% of a cycle), leaving sufficient power and capacity for scheduling.
Through the operation control of the electrolytic hydrogen production system and the electrochemical energy storage system, the renewable energy power generation system has achieved higher economic benefits under the premise of completing multiple target application scenarios. Among them, the income from PSAS is CNY 524,500, the income from carbon emission reduction is CNY 297,325.2, the income from hydrogen sales is CNY 89,981.55, the cost of electricity purchase is CNY 87,909.82, the penalty for abandoning wind and solar energy is CNY 13,851, and the comprehensive benefit is CNY 810,045.93.
As shown in Figure 9, after calculation, it is found that there is a positive under-compensation of 48.78 MW and a negative under-compensation of 31.18 MW before suppression. The power fluctuates greatly between 9:00 and 7:00, with a total of 14 time periods exceeding the fluctuation limit of renewable energy grid-connected power, and the over-limit probability is 14.58%, which brings a significant negative impact to the power system; after suppression by the EHPS and the electrochemical energy storage system, all the renewable energy outputs participating in the grid connection meet the grid-connected power fluctuation limit.

6.5. Effect of Electrolytic Hydrogen Production on Economic Benefits

The effectiveness analysis of the three scenarios—no electrolytic hydrogen with electrochemical energy storage, electrolytic hydrogen without electrochemical energy storage, and electrolytic hydrogen with electrochemical energy storage—are shown in Table 6. The installed capacity of equipment under the three scenarios is as follows:
Scenario 1 (no electrolytic hydrogen with electrochemical energy storage): no electrolyzer, fuel cell, or hydrogen storage tank is set up, and the electrochemical energy storage is 25 MW;
Scenario 2 (electrolytic hydrogen without electrochemical energy storage): in this scenario, the operation and maintenance expenses of the 25 MW electrochemical energy storage in Scenario 3 is proportionally allocated to the EHPS and the fuel cell, among which the electrolyzer is 30 MW, the fuel cell is 10 MW, the hydrogen storage tank is 36,483 N m 3 , and no electrochemical energy storage is set up;
Scenario 3 (electrolytic hydrogen with electrochemical energy storage): the electrolyzer is 25 MW, the fuel cell is 10 MW, the hydrogen storage tank is 30,000 N m 3 , and the electrochemical energy storage is 20 MW. This scenario is the scenario of the proposed mathematical model.
As evident from Table 6, the comprehensive benefit of Scenario 3 proposed in this paper is optimal. Compared with Scenario 1, due to the constraint that electrochemical energy storage cannot perform simultaneous charging and discharging operations, it is impossible to flexibly absorb or release electricity in Scenario 1, and the hydrogen sale benefits brought by electrolytic hydrogen is lacking. Therefore, compared with Scenario 1, Scenario 3 reduces the amount of wind and solar energy discarded by 393.8 MW, reduces the penalty for discarding wind and solar energy by CNY 118,140, increases carbon trading benefits by CNY 10,218.5, and increases comprehensive benefits by CNY 123,247.6. Compared with Scenario 2, due to the lack of electrochemical energy storage and the lower conversion rate of electrolytic hydrogen, Scenario 3 reduces the amount of wind and solar energy discarded by 262.3 MW compared to Scenario 2, reduces the penalty for discarding wind and solar energy by CNY 78,690, and increases comprehensive benefits by CNY 49,214.4. Scenario 3 combines the flexibility of the electrolytic hydrogen production system and the higher conversion rate of the electrochemical energy storage system. Under the premise of completing the multi-objective application scenarios of smoothing power fluctuations, participating in PSAS, and reducing the amount of wind and solar energy discarded, it can still obtain higher comprehensive benefits.

7. Conclusions

This paper focuses on a renewable energy power plant that includes hydrogen production through electrolysis. It considers the use of electrolytic hydrogen to smooth out fluctuations in renewable energy power generation, participate in peak-shaving auxiliary services, and increase the absorption space of renewable energy. It proposes a multi-objective energy management model and algorithm based on the co-ordinated optimization of electrolytic hydrogen and renewable energy power generation. Through simulation, based on the analysis of the operating conditions and safety factors of various facilities in the power station, the advantages of this model for power system safety were carefully studied. With participation in the power market, the system’s economics have been greatly improved.
At the same time, the PLPSO algorithm improves the convergence accuracy, convergence speed, and stability of the proposed model compared with the basic algorithm.
The research presented in this paper offers a valuable theoretical reference for energy management in renewable energy-generating plants that include an electrolytic hydrogen production system.
The next step will be to further study the capacity configuration of renewable energy power stations that include an electrolytic hydrogen production system, providing technical support for the installation and construction of renewable energy power stations.

Author Contributions

Conceptualization, T.S.; methodology, T.S. and L.G.; software, L.G.; validation, L.G., J.S. and Z.X.; formal analysis, T.S.; investigation, L.G.; resources, T.S.; data curation, L.G.; writing—original draft preparation, L.G.; writing—review and editing, T.S.; visualization, L.G.; supervision, T.S.; project administration, T.S.; funding acquisition, T.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the major project of National Natural Science Foundation of China: “Theory and method of coordinated operation optimization of multi energy supply and demand system containing hydrogen”, grant number 62192753.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Li, Z.; Zhang, R.; Sun, H.X. Review on Key Technologies of Hydrogen Generation, Storage and Transportation Based on Multi-Energy Complementary Renewable Energy. Trans. China Electrotech. Soc. 2021, 36, 446–462. [Google Scholar]
  2. Zhang, P.C.; Xu, J.; Ke, D.P. Optimal Configuration of Energy System in Iron and Steel Park Driven by Hydrogen Energy. Autom. Electr. Power Syst. 2022, 46, 1–10. [Google Scholar]
  3. Zhang, H.; Yuan, T.J.; Tan, J. Medium and Long-term Forecast of Hydrogen Load in Unified Energy System. Proc. CSEE 2021, 41, 3364–3372+3662. [Google Scholar]
  4. Cai, G.W.; Kong, L.G.; Xue, Y. Overview of Research on Wind Power Coupled with Hydrogen Production Technology. Autom. Electr. Power Syst. 2014, 38, 127–135. [Google Scholar]
  5. Jing, P.; Xu, G.Z.; Zhao, B. Large-scale Energy Storage Technology for Global Energy Internet. Smart Grid 2015, 3, 486–492. [Google Scholar]
  6. Fard, H.H.; Tooryan, F.; Collins, E.R.; Jin, S.; Ramezani, B. Design and optimum energy management of a hybrid renewable energy system based on efficient various hydrogen production. Int. J. Hydrogen Energy 2020, 45, 30113–30128. [Google Scholar]
  7. Jin, X.; Zhuang, Y.X.; Wang, H.; Zhou, S.Y. Feasibility Analysis Research on Abandoning Wind and Solar Energy with Hydrogen Energy Storage Technology. Electrotech. Electr. 2019, 4, 63–68. [Google Scholar]
  8. Pu, Y.C.; Li, Q.; Chen, W.R.; Huang, W.Q.; Hu, B.B.; Han, Y.; Wang, X. Energy Management for Islanded DC Microgrid with Hybrid Electric-hydrogen Energy Storage System Based on Minimum Utilization Cost and Energy Storage State Balance. Power Syst. Technol. 2019, 43, 918–927. [Google Scholar]
  9. Liang, T.; Sun, B.; Mi, D.; Sha, Z. Operation Optimization of Wind-solar Complementary Hydrogen Production System Based on MOGEO. In Proceedings of the 7th International Conference on Power and Renewable Energy (ICPRE), Shanghai, China, 23–26 September 2022. [Google Scholar] [CrossRef]
  10. Cui, Y.; Yan, S.; Zhong, W.Z.; Wang, Z.; Zhang, P.; Zhao, Y.T. Optimal Thermoelectric Dispatching of Regional Integrated Energy System with Power-to-gas. Power Syst. Technol. 2020, 44, 4254–4264. [Google Scholar] [CrossRef]
  11. Deng, J.; Jiang, F.; Wang, W.Y.; He, G.X.; Zhang, X.H.; Liu, K.C. Low-carbon Optimized Operation of Integrated Energy System Considering Electric-heat Flexible Load and Hydrogen Energy Refined Modeling. Power Syst. Technol. 2022, 46, 1692–1704. [Google Scholar] [CrossRef]
  12. Shi, R.F.; Ning, J.; Gao, L.Q. Research on Optimal Dispatch Strategy of Wind and Solar Self-consistent Microgrid in Road Transportation System with Hydrogen Energy Storage. Acta Energiae Solaris Sin. 2023, 44, 513–521. [Google Scholar] [CrossRef]
  13. Mah, A.; Ho, W.S.; Hassim, M.H.; Hashim, H.; Ling, G.H.; Ho, C.S.; Muis, A. Optimization of a standalone photovoltaic-based microgrid with electrical and hydrogen loads. Energy 2021, 235, 121218. [Google Scholar] [CrossRef]
  14. Wu, F.B. Optimized Configuration of Wind-Solar Hydrogen Storage Integrated Energy System. Shandong Chem. Ind. 2020, 49, 135–136+138. [Google Scholar] [CrossRef]
  15. Guo, M.J.; Yan, Z.; Zhou, Y. Optimized Operation Design of Integrated Energy System with Wind Power Hydrogen Production. Electr. Power 2020, 53, 115–123+161. [Google Scholar]
  16. Jiang, Y.W.; Deng, Z.H.; Shi, Y. Size optimization and economic analysis of a coupled wind-hydrogen system with curtailment decisions. Int. J. Hydrogen Energy 2019, 44, 19658–19666. [Google Scholar] [CrossRef]
  17. Cao, W.; Zhong, X.; Wang, H.H. Optimal Capacity Allocation of Hydrogen Production System Participating Peak Regulation for Auxiliary with Thermal Power Plant. Distrib. Energy 2020, 5, 15–20. [Google Scholar] [CrossRef]
  18. Cheng, S.; Zhang, J.C.; Han, H.T.; Wei, Z.N.; Sun, G.Q.; Zhou, Y.Z. Optimal Dispatch of Electricity-Gas-Hydrogen Integrated Energy System Considering Auxiliary Service. Autom. Electr. Power Syst. 2023, 47, 110–120. [Google Scholar]
  19. Jia, C.Z.; Wang, L.M.; Meng, E.L. Optimal Capacity Configuration and Day-Ahead Scheduling of Wind-Solar-Hydrogen Coupled Power Generation System. Electr. Power 2020, 53, 80–87. [Google Scholar]
  20. Harasis, S.; Sozer, Y.; Elbuluk, M. Reliable Islanded Microgrid Operation Using Dynamic Optimal Power Management. IEEE Trans. Ind. Appl. 2021, 57, 1755–1766. [Google Scholar] [CrossRef]
  21. Son, Y.-G.; Kim, S.-Y.; Bae, I.-S. Optimal Coordination of Energy Coupling System Considering Uncertainty of Renewable Energy Sources. Energies 2024, 17, 931. [Google Scholar] [CrossRef]
  22. Zhu, Y.; Lan, Z.B.; Kui, Z. Research on optimal operation of wind-PV-ES complementary system considering carbon emission cost. Power Syst. Prot. Control 2019, 47, 127–133. [Google Scholar]
  23. Yuan, T.J.; Guo, J.H.; Yang, Z.J. Optimal Allocation of Power Electric-Hydrogen Hybrid Energy Storage of Stabilizing Wind Power Fluctuation. Proc. CSEE 2023, 1–10. [Google Scholar] [CrossRef]
  24. Cui, Y.; Zhou, H.J.; Zhong, W.Z.; Li, H.B.; Zhao, Y.T. Low-carbon scheduling of power system with wind power considering uncertainty of both source and load sides. Electr. Power Autom. Equip. 2020, 40, 85–93. [Google Scholar] [CrossRef]
  25. Liu, B.D.; Zhao, R.Q.; Wang, G. Uncertain Programming with Applications; Tsinghua University Press: Beijing, China, 2003; pp. 178–198. [Google Scholar]
  26. Wang, Q.X.; Guo, X.B. Particle swarm optimization algorithm based on Levy flight. Appl. Res. Comput. 2016, 33, 2588–2591. [Google Scholar]
  27. Zhang, C.; He, X.S.; Ye, Y.R. Particle swarm optimization algorithm based on elite strategy and levy flight. J. Xi’an Polytech. Univ. 2018, 32, 731–738. [Google Scholar]
Figure 1. Renewable energy power plant with electrolytic hydrogen production system.
Figure 1. Renewable energy power plant with electrolytic hydrogen production system.
Processes 12 00541 g001
Figure 2. The process of the PLPSO algorithm.
Figure 2. The process of the PLPSO algorithm.
Processes 12 00541 g002
Figure 3. Forecast output curve of wind power and photovoltaic on a typical day.
Figure 3. Forecast output curve of wind power and photovoltaic on a typical day.
Processes 12 00541 g003
Figure 4. Program output and peak shaving command curve.
Figure 4. Program output and peak shaving command curve.
Processes 12 00541 g004
Figure 5. Comparison of convergence curves of two algorithms.
Figure 5. Comparison of convergence curves of two algorithms.
Processes 12 00541 g005
Figure 6. Peak-shaving subsidy and peak-shaving electricity curve.
Figure 6. Peak-shaving subsidy and peak-shaving electricity curve.
Processes 12 00541 g006
Figure 7. Operation of electrolytic cell, fuel cell, and hydrogen storage tank.
Figure 7. Operation of electrolytic cell, fuel cell, and hydrogen storage tank.
Processes 12 00541 g007
Figure 8. Operation of electrochemical energy storage.
Figure 8. Operation of electrochemical energy storage.
Processes 12 00541 g008
Figure 9. The grid-connected power fluctuation before and after stabilization.
Figure 9. The grid-connected power fluctuation before and after stabilization.
Processes 12 00541 g009
Table 1. Peak load regulation pricing details.
Table 1. Peak load regulation pricing details.
PeriodQuotation FileLoad Rate of Thermal Power PlantLower Quotation Limit (CNY/kWh)Upper Quotation Limit (CNY/kWh)
non-heating seasonfirst gear40% < load rate 48%00.4
second gearload rate 40%0.41
heating seasonfirst gear40% < load rate 50%00.4
second gearload rate 40%0.41
Table 2. Operation parameters.
Table 2. Operation parameters.
ParameterValue
P wt min /MW0
P wt m ax /MW25
P pv min /MW0
P pv max /MW200
η bat c h , η bat d i s c h 0.9
S O C min 0.2
S O C max 0.8
The initial state of SOC0.5
E b a t r a t e d /MW20
P ele min /MW6.25
P ele m ax /MW25
η ele / ( N m 3 / M W h ) 190
η fc / ( N m 3 / M W h ) 555.56
V h s max / N m 3 30,000
V h s ( 1 ) / N m 3 2000
Table 3. Economic parameters.
Table 3. Economic parameters.
ParameterValue
A hyd /(CNY/ N m 3 )3
A c o 2 /(CNY/kg)0.075
λ p /(kg/MWh)500
λ q /(kg/MWh)798
γ fc /(CNY/MWh·Day)1136
γ hs /(CNY/ N m 3 ·Day)0.5
γ ele /(CNY/MW·Day)1748
γ bat /(CNY/MW·Day)850
γ loss /(CNY/MWh)300
A b u y (0:00 to 8:00)/(CNY/MWh)135.0
A b u y (9:00 to 12:00 and 18:00 to 23:00)/(CNY/MWh)521.4
A b u y (others)/(CNY/MWh)328.2
Table 4. Comparison between traditional PSO and PLPSO.
Table 4. Comparison between traditional PSO and PLPSO.
PSOPLPSO
Number of runs100100
average run time/s261.3248.4
mean value−654,818.07−768,875.44
optimal value−807,479.2−830,172.1
mean absolute deviation (MAD)62,614.617,792.2
Table 5. Scheduling results under different confidence levels.
Table 5. Scheduling results under different confidence levels.
Confidence LevelProportion of the Wind and Solar Energy AbandonComprehensive Benefit/CNY
0.559.97%720,723
0.607.17%714,133
0.655.36%731,836.8
0.702.06%736,225.2
0.752.32%730,571.3
0.801.98%738,359.2
0.851.25%748,917.6
0.901.42%762,298.3
0.951.16%774,511.2
1.000.67%784,582.3
Table 6. Effectiveness analysis in three scenarios.
Table 6. Effectiveness analysis in three scenarios.
Scenario 1Scenario 2Scenario 3
electricity purchasing cost/CNY28,280.571,492.487,909.8
abandoned wind power and solar/MW500.0308.546.2
Hydrogen sales revenue/CNY0109,449.089,981.6
Peaking revenue/CNY524,500512,500524,500
carbon trading income/CNY297,789.9302,926.6297,325.2
comprehensive benefit/CNY644,002.6760,831.5810,045.9
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shi, T.; Gu, L.; Xu, Z.; Sheng, J. Research on Multi-Objective Energy Management of Renewable Energy Power Plant with Electrolytic Hydrogen Production. Processes 2024, 12, 541. https://doi.org/10.3390/pr12030541

AMA Style

Shi T, Gu L, Xu Z, Sheng J. Research on Multi-Objective Energy Management of Renewable Energy Power Plant with Electrolytic Hydrogen Production. Processes. 2024; 12(3):541. https://doi.org/10.3390/pr12030541

Chicago/Turabian Style

Shi, Tao, Libo Gu, Zeyan Xu, and Jialin Sheng. 2024. "Research on Multi-Objective Energy Management of Renewable Energy Power Plant with Electrolytic Hydrogen Production" Processes 12, no. 3: 541. https://doi.org/10.3390/pr12030541

APA Style

Shi, T., Gu, L., Xu, Z., & Sheng, J. (2024). Research on Multi-Objective Energy Management of Renewable Energy Power Plant with Electrolytic Hydrogen Production. Processes, 12(3), 541. https://doi.org/10.3390/pr12030541

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop