Next Article in Journal
The Impact of the Rule of Law on Energy Policy in European Union Member States
Previous Article in Journal
Precision Leak Detection in Supermarket Refrigeration Systems Integrating Categorical Gradient Boosting with Advanced Thresholding
Previous Article in Special Issue
Operation Optimization of Wind/Battery Storage/Alkaline Electrolyzer System Considering Dynamic Hydrogen Production Efficiency
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dual-Stage Optimization Scheduling Model for a Grid-Connected Renewable Energy System with Hybrid Energy Storage

1
Powerchina Huadong Engineering Corporation, Hangzhou 310030, China
2
College of Electrical Engineering, Zhejiang University, Hangzhou 310058, China
*
Author to whom correspondence should be addressed.
Submission received: 29 November 2023 / Revised: 30 January 2024 / Accepted: 2 February 2024 / Published: 4 February 2024
(This article belongs to the Special Issue Wind/PV/Hydrogen Integrated Energy System for a Clean Future)

Abstract

:
To operate the grid-connected renewable energy system economically, this study presents a dual-stage optimization scheduling model for grid-connected systems with hybrid energy storage, including day-ahead and intra-days stages. In the day-ahead stage, an economically optimal scheduling model is developed, considering the price peak-to-valley difference. This model aims to enhance the economic efficiency of the system by utilizing hybrid energy storage. In the intra-day stage, more accurate renewable energy forecasts with a shorter time scale are considered. The objectives are to minimize the curtailment rate of renewable energy and to track the day-ahead scheduling outcomes. The NSGA-II algorithm is employed for multi-objective optimization, achieving equilibrium solutions considering multiple optimization objectives. Compared to other published works, the proposed model achieves a balance between different optimization objectives, enabling the system to operate economically and stably. It provides a comprehensive approach to optimize the scheduling of grid-connected systems with hybrid energy storage by considering both economic and operational aspects. Overall, this proposed dual-stage optimization model presents a viable approach to improve economic efficiency and mitigate renewable energy curtailment in grid-connected systems. By effectively integrating renewable energy sources and optimizing their utilization, this model contributes to enhancing the sustainability and optimal operation of the power grid.

1. Introduction

With the gradual depletion of fossil fuels and increasing calls for environmental protection, the substitution of renewable energy for fossil fuels as the source of the power grid has emerged as a solution to address energy and environmental issues [1]. However, due to the volatility and low inertia of renewable energy sources, their integration into the grid can lead to the source–load mismatch, necessitating the involvement of energy storage to achieve source–load balance and maintain the stability of the power grid [2].
Energy storage can be classified into power-based storage and energy-based storage, based on their characteristics [3]. Generally speaking, batteries are considered as power-type energy storage due to their high power density and short storage duration. On the other hand, energy-type storage methods, like hydrogen storage, have longer storage durations but cannot accommodate rapid power fluctuations. Hydrogen storage, as a form of energy-based storage, has gained increasing attention due to the rise of the hydrogen industry and hydrogen-powered vehicles as a means to balance the volatility of renewable energy [4].
The utilization of electrolyzers for electrochemical hydrogen generation represents an environmentally sustainable approach [5]. Electrolyzers encompass various types, including alkaline electrolyzers, proton exchange membrane electrolyzers, and solid oxide electrolyzers. Among these, alkaline electrolyzers have gained widespread application due to their well-established technology and cost-effectiveness [6]. In recent years, the integration of alkaline electrolyzers with lithium batteries for co-generation and energy storage has garnered considerable research interest. This stems from the complementary nature of hydrogen storage, offering high energy density, and lithium batteries, providing high power density [7]. Furthermore, the rapid reaction rates exhibited by lithium batteries can compensate for the slower response times of electrolyzers, enhancing overall system performance. Researchers have explored power allocation for hybrid energy storage systems (HESSs) from the perspectives of frequency decomposition, cost-effectiveness, and equipment characteristics to achieve optimal objectives. An optimal day-ahead optimization schedule for a gas-electric integrated energy system considering the bidirectional energy flow is proposed in [8]. Furthermore, Hong et al. [9] presented an optimization day-ahead scheduling model of wind-hydrogen systems considering hydrogen production efficiency.
Indeed, the studies mentioned above focus on day-ahead optimization. In day-ahead optimization, the time scale for source–load forecasting is typically one day, which often leads to lower prediction accuracy. This limitation can impact the effectiveness of implementing day-ahead optimization results. Due to the inherent variability and uncertainties associated with renewable energy sources and load patterns, accurately forecasting the source–load dynamics over a one-day horizon can be challenging. Factors such as weather conditions, demand fluctuations, and unforeseen events can introduce deviations between the predicted and actual source–load profiles. The lower prediction accuracy in day-ahead optimization can affect the realization of the intended outcomes, as the system may encounter unexpected variations in renewable energy availability or demand patterns [10]. Therefore, the dual-stage optimization model should be applied to the operation of renewable energy systems.
To date, several studies have addressed the dual-stage scheduling of day-ahead and intra-day operations. For example, Wang et al. [11] proposed a two-stage energy management model for the sustainable wind–PV–hydrogen–storage microgrid based on receding horizon optimization. Similarly, Yuan et al. [12] presented a two-time-scale microgrid energy management model for scheduling with low operational costs and high reliability against uncertainties. In [13], the authors proposed a novel multi-energy systems optimization model to maximize investment and operating synergy in the electricity, heating, and transport sectors.
However, these studies predominantly employed single-objective optimization approaches. In single-objective optimization, pursuing a specific target may lead to adverse effects on other objectives. In practice, the intra-day operational process involves numerous factors that require optimization, including wind and solar curtailment rates, grid power tracking, operational costs, and other relevant variables. In optimizing renewable energy power systems, reducing the curtailment rate of wind and solar energy may increase the demand for energy storage, thus raising operational costs [6]. Treating these factors as components of a multi-objective optimization problem is a more appropriate approach. By incorporating multiple objectives, the optimization framework can better capture the complex trade-offs and interdependencies among different operational aspects, leading to more comprehensive and effective decision-making strategies [14].
Multiple multi-objective optimization algorithms, such as the non-dominated sorting genetic algorithm (NSGA-II), moth–flame optimization (MFO), and multi-objective particle swarm optimization (MOPSO), have been proposed in recent years, providing effective tools for solving complex multi-objective problems and offering a range of optimal solutions based on different preferences [15,16]. The MOPSO is employed to optimize the operation of microgrids in [17]. The MFO algorithm is utilized for wind farm layout optimization in [18]. Anosri et al. [19] compared different classes of multi-objective optimization algorithms in the context of reliability design optimization. Among the multi-objective optimization algorithms, one representative is the non-dominated sorting genetic algorithm II (NSGA-II). NSGA-II, an advanced algorithm capable of handling multi-objective optimization, has been utilized by researchers for economic optimization in power systems-related studies. Rejeb et al. [20] applied NSGA-II to estimate the optimal results of the proposed multi-generation system. In [21], a suggested system producing electricity and hydrogen was planned in a daily triple-periodic framework and solved by NSGA-II. With NSGA-II, Maheri et al. [22] designed reproduction operators for robust exploration and exploitation at both size and configuration levels. Since NSGA-2 has a wide range of applications, it is used in this paper to solve the multi-objective problem in the intra-day stage.
This paper aims to apply HESS for the power balance and economic efficiency of renewable energy systems. The main contributions of this work include the following:
  • We propose a dual-stage optimization scheduling model for grid-connected systems with HESS to balance renewable energy, with day-ahead and intra-day stages.
  • In the day-ahead stage, we propose a scheduling model that utilizes the price peak-to-valley difference and HESS to enhance the system’s efficiency. In the intra-day stage, we introduce a multi-objective model, and NSGA-II is applied as a solution algorithm to achieve the optimization of multiple objectives simultaneously.
  • Case studies show that for a real renewable energy system with HESS, the proposed scheduling model reduces its daily costs by around 6% compared with the rule-based operation. Model comparisons are conducted to compare the proposed model and other existing models.
The remainder of this paper is organized as follows. In Section 2, the system modeling and equipment characteristics are described. Section 3 presents the source–load forecasting method. The dual-stage optimization scheduling model is proposed in Section 4. Case studies are presented in Section 5. Finally, Section 6 concludes the paper.

2. System Modeling

2.1. System Description

Figure 1 depicts the structure of a renewable energy system augmented with HESS. In this system, renewable energy sources include wind and solar energy. The majority of the load is supplied by the output of renewable energy sources, which consists of both electrical load and hydrogen load. The energy balance of the system is maintained by a HESS composed of batteries and electrolyzers. The HESS plays a crucial role in two aspects:
  • Participating in the day-ahead optimal system scheduling, utilizing the price difference in the electricity grid to achieve economic optimization.
  • Balancing the source–load imbalance that may occur during the day-ahead scheduling.
Unlike the transient power characteristics of the source load, the state and actions of energy storage devices, such as battery charging and discharging, and power adjustment of electrolyzers, are time-dependent. Therefore, a comprehensive understanding and optimal utilization of the HESS is essential [23].
This approach assists in resolving the intra-day supply-demand mismatch caused by the prediction errors in source–load forecasting, as the energy exchange between the system and the grid needs to align with the results of the day-ahead scheduling.

2.2. Wind Turbine

The output power of a wind turbine can be expressed as (1) [24].
p w i n d ( v w ) = 0 , 0 v w < v i n v w v i n v r a t e d v i n × p w i n d r a t e d , v i n v w < v r a t e d p w i n d r a t e d , v r a t e d v w < v o u t 0 , v w v o u t
where v i n and v o u t are the cut-in and cut-out wind velocities, respectively. v r a t e d and p w i n d r a t e d represent the rated wind velocity and power.

2.3. Photovoltaic

The output power of a photovoltaic power generation device is related to the solar irradiation intensity and the ambient temperature, which can be calculated as (2) [25]:
p p v ( G p v , T e ) = p p v r a t e d f p v G p v G S T C ( 1 α p ( T e T r e f ) )
where p p v r a t e d represents the rated power output (kW) of the PV under standard test condition (STC), G S T C and T r e f are the STC irradiation and temperature, G p v represents the real-time solar irradiation intensity, f p v is the discount factor of the power output, α p is the temperature coefficient, and T e denotes the real-time panel temperature.
T e ( G p v ) = T a + G p v · T c , N O C T T a , N O C T G N O C T
The panel temperature T e can be determined by (3), where T a is the actual ambient temperature, T c , N O C T and T a , N O C T are the nominal operating cell temperature (NOCT) of the PV panels and the NOCT of the ambient temperature, and G N O C T is the NOCT of irradiation intensity [25].

2.4. Hybrid Energy Storage System

In a hybrid energy storage system, the storage state of the battery is determined by the previous moment’s storage state and the actions taken at the current moment. This paper takes into account the energy dissipation in the battery and expresses its storage state according to (4).
S o C ( t ) = S o C ( t δ t ) · ( 1 γ b ) + ( η b c h · p b c h ( t ) p b d i s ( t ) η b d i s ) · δ t / Q b R a t e d
where S o C ( t ) is the state of charge at time t, γ b is the dissipation rate of the battery, p b c h ( t ) and p b d i s ( t ) are the charge and discharge power of the battery at time t, and Q R a t e d b is the rated storage capacity of the battery, which is one of the variables to be optimized.
Compared to the direct energy storage process of batteries, hydrogen storage systems involve the conversion of electrical energy into hydrogen energy through electrolyzers before storage. This process is slower and less efficient than battery charging and discharging. However, hydrogen storage systems have lower unit energy storage construction costs than batteries and offer excellent environmental benefits. Additionally, they can meet users’ hydrogen energy demands, like fuel cell vehicles or industrial processes [2]. Equations (5) and (6) illustrate the chemical process of electrolysis, converting electricity to hydrogen.
H 2 evolution reaction in the negative pole:
H 2 + 2 O H 2 H 2 O = 2 e
O 2 evolution reaction in the positive pole:
2 O H 1 2 O 2 H 2 O = 2 e
Equation (7) explains the mathematical model for the electrolyzer [6].
m e l z o u t = η e l z p e l z L H V H 2
where m e l z o u t represents the hydrogen generated from the electrolyzer, p e l z represents the electricity consumed by the electrolyzer, and L H V H 2 denotes the low heating value of hydrogen gas.
The hydrogen gas produced by the electrolyzer is compressed by a compressor to the appropriate pressure and then input into the hydrogen storage tank. The storage state of the hydrogen tank is also dependent on the previous moment’s state and the production/output at the current moment, as indicated by (8) [24].
L o H ( t ) = L o H ( t δ t ) · ( 1 γ h ) + ( η h i n m e l z o u t ( t ) m h o u t ( t ) η h o u t ) · δ t / Q h R a t e d
where L o H ( t ) is the state of hydrogen of a storage tank at time t, γ h is the hydrogen dissipation rate of the tank, m h o u t ( t ) is the hydrogen output rate of the tank, δ t is the time step of optimization, and Q h R a t e d is the rated storage capacity of the tank.

3. Source–Load Forecasting Method

The electrical load and renewable energy generation, such as wind and solar power, often exhibit complex temporal characteristics, and a single network structure may not achieve satisfactory predictive performance. This paper proposes the use of a convolutional neural network-long short-term memory (CNN-LSTM) network for source–load forecasting. The CNN network is commonly used to process raw data, automatically extracting internal features and reducing algorithm complexity. The LSTM network, with its memory units, is effective in handling long-term sequential data [26].
In the CNN-LSTM network, historical source–load data and meteorological data are used as inputs. The CNN layer extracts features, followed by pooling operations to reduce dimensionality. The fully connected layer then transforms the features into a one-dimensional structure, extracting the feature vector. The LSTM layer learns the internal patterns of load variations from the extracted features to achieve the forecasting function. Finally, the output layer provides the prediction results, corresponding to the predicted source–load sequence [27].
The proposed forecasting method starts with a four-layer one-dimensional convolutional network for feature extraction from the input source–load sequence. The convolutional network uses larger kernel sizes to extract coarse-grained temporal features initially and gradually reduces the kernel size to extract fine-grained temporal features. The number of channels is increased as the network deepens. The kernel sizes for the 4 convolutional layers are 24, 12, 8, and 2, and the corresponding channel numbers are 132, 232, 332, and 432. Batch normalization is applied after each convolutional layer to stabilize data distribution and accelerate convergence. Additionally, a dropout layer with a dropout rate of 0.25 is added after each convolutional layer to prevent overfitting.
The data extracted by the convolutional layers are then fed into the LSTM network for model training. The LSTM network consists of 2 LSTM layers, with each layer having 256 neurons. After training, the model’s predictions are obtained by passing the results through a fully connected layer.
In this paper, the input for day-ahead source–load forecasting consists of the source–load data and weather data from the previous week, as well as the weather forecast for the target day. The output is the hourly source–load data for the target day.
For intra-day forecasting, the input includes the source–load data and weather data for the next 48 h, along with the weather forecast for the next hour. The output is the minute-level source–load data for the next hour.

4. Dual-Stage Optimization Scheduling Model

The flowchart of the proposed dual-stage optimization scheduling model is presented in Figure 2. The white blocks denote the preparation of the model, including data acquisition, forecasting, and device modeling. The green blocks illustrate the day-ahead scheduling model, while the yellow blocks are the intraday scheduling model. Blue blocks denote the scheduling results. The architecture diagram of the proposed two-stage optimization scheduling model is shown in Figure 3. The model consists of two stages: day-ahead and intra-day optimization. In the day-ahead stage, the time scale is one day with a resolution of one hour and the optimization goal is to minimize the daily operation costs. In the intra-day stage, the time scale is one hour, and it utilizes the day-ahead optimization results as a reference. The intra-day optimization stage aims to compensate for power deviations caused by inaccurate source–load forecasting using HESS balancing. The resolution of the intra-day optimization is one minute. The optimization variables of the proposed model are the power of the electrolyzer, the charging or discharging power of the battery, and the power exchange with the utility grid.

4.1. Day-Ahead Optimal Scheduling Model

4.1.1. Optimization Goal

The objective of the day-ahead optimization stage is to minimize the daily operating costs. In the proposed renewable energy system, the proposed renewable energy system, these costs primarily include electricity price costs, equipment operation and maintenance costs, battery degradation costs, and penalty costs incurred due to wind and solar power curtailment. The operating revenue mainly consists of income from selling electricity to the grid and hydrogen sales. The daily costs can be expressed as (9).
O b j d h = C e p u r + C o m + C b d e + C r e p a l B e s e l l B h s e l l
The electricity cost benefits can be expressed as (10). The impact of electricity prices on costs is significant. This paper aims to optimize the economic efficiency of the system by utilizing the price difference between peak and off-peak hours through day-ahead power scheduling.
C e p u r = t = 1 T c e p u r ( t ) p e p u r ( t )
The equipment operation and maintenance costs can be expressed as (11). The battery degradation costs are related to the battery’s charging and discharging power, which can be expressed as (12).
C o m = t = 1 T o e l z p e l z ( t ) + o c m p p c m p ( t ) + o p v p p v ( t ) + o w i n d p w i n d ( t )
C b d e = t = 1 T d b p b c h ( t ) + d b p b d i s ( t )
Equation (13) illustrates the penalty costs incurred due to wind and solar power curtailment.
C r e p a l = t = 1 T c r e p a l p r e p a l ( t )
Given the interconnected nature of the system, it is imperative to account for the revenue generated by exporting electricity to the grid, as expressed in (14). Nonetheless, it should be noted that the selling price of electricity to the grid is typically subject to temporal variations and tends to be significantly lower than the purchasing price. This discrepancy arises due to the inherent risk associated with integrating electricity from the renewable energy system into the power grid, as it may potentially impact grid stability.
B e s e l l = t = 1 T c e s e l l ( t ) p e s e l l ( t )
Furthermore, the system can also generate revenue by selling hydrogen energy, as shown in (15).
B h s e l l = t = 1 T c h s e l l m h s e l l ( t )

4.1.2. Constraints

Power and hydrogen balances should be satisfied during the operation, as expressed in (16) and (17).
p p v ( t ) + p w i n d ( t ) = p l o a d ( t ) + p e l z ( t ) + p c m p ( t ) + p b c h ( t ) p b d i s ( t ) + p g e x ( t ) + p r e p a l ( t )
m h s e l l ( t ) = m h o u t ( t )
The operational power of the compressor is indeed related to factors such as the hydrogen production rate of the electrolyzer and the pressure of the hydrogen storage tank. In this paper, it is assumed that the operational power of the compressor is directly proportional to the hydrogen production rate, as expressed in (18).
p c m p ( t ) = β m e l z o u t ( t )
In the day-ahead optimization, it is necessary to ensure the safe and stable operation of the equipment. This necessitates operating the devices within certain power ranges.
p e l z min p e l z ( t ) p e l z max
p c m p min p c m p ( t ) p c m p max
p b c h min p b c h ( t ) p b c h max
p b d i s min p b d i s ( t ) p b d i s max
p g e x min p g e x ( t ) p g e x max
The states for both the battery and hydrogen storage tanks also need to be constrained. Overcharging or over-discharging can lead to a decrease in battery lifespan. Additionally, a portion of the storage capacity needs to be reserved for intra-day optimization during day-ahead scheduling. Therefore, the lower and upper limits for the storage state are set at 0.2 and 0.8, respectively [28].
0.2 S o C ( t ) 0.8
0.2 L o H ( t ) 0.8
In order to maximize the energy utilization efficiency of the system, the battery is not allowed to charge and discharge simultaneously.
p b c h ( t ) · p b d i s ( t ) = 0
Ultimately, in order to ensure optimal energy utilization, it is necessary to limit the curtailment rate of wind and solar power within a certain range, as indicated in (27).
t t = T p r e p a l ( t ) / t t = T p w i n d ( t ) + p p v ( t ) ε r e c u l

4.2. Intra-Day Optimal Scheduling Model

4.2.1. Optimization Goal

The main objective of the intra-day optimization stage is to align the energy exchange between the system, the energy flow of EHSS, and the grid with the day-ahead scheduling results while minimizing the curtailment rate of wind and solar power. As a result, the intra-day optimization model becomes a multi-objective model with the objective functions represented by (28).
O b j i d = E g r d , E c u l , E h e s s E g r d = τ τ = λ | p g e x p ˜ g e x ( τ ) | / ( λ p g e x ) E c u l = τ τ = λ | p ˜ r e p a l ( τ ) | / ( λ p p v + λ p w i n d ) E h e s s = τ τ = λ ( | p b p ˜ b | + | p e l z p ˜ e l z | ) / ( λ p b + λ p e l z )
where τ represents the optimization time interval of intra-day operation, while λ represents the optimization time period.

4.2.2. Constraints

Similar to the day-ahead operation, during the intra-day operation of the system, several constraints also need to be considered, including the power and hydrogen energy balance constraint, as specified in (16) and (17), the compressor’s power constraint, as specified in (18), the equipment operating range constraint, as specified in (19)–(23), storage constraints for lithium batteries and hydrogen tanks, as specified in (24) and (25), as well as the constraint that prevents simultaneous charging and discharging of the batteries, as specified in (26).

4.3. Solution Method

The day-ahead optimization scheduling model can be formulated as a mixed-integer linear programming (MILP) problem and is solved directly by applying Matlab R2022a with Gurobi 10.0. However, the intra-day scheduling is a multi-objective optimization problem. In this study, the non-dominated sorting genetic algorithm II (NSGA-II) method is employed to solve the intra-day optimization problem.
The NSGA-II algorithm, originally proposed by Deb et al. [29], is a well-structured optimization algorithm widely utilized for solving complex engineering problems. It incorporates several key components such as fast-non-dominated sorting, a smart maintenance approach, and effective crowding distance assessment, which collectively contribute to improving the convergence rate and maintaining population diversity. The solving process of intra-day multi-objective optimization problems based on NSGA-II is shown in Figure 4. During the initialization stage, careful consideration is given to the limitations of decision variables for optimization. If an original population is available, the algorithm proceeds to the generation of the second population; otherwise, it performs the non-dominated selection, followed by selection, crossover, and mutation steps to create the original population. A merging technique is employed to combine the initial and offspring populations. If a new parent population is generated, the algorithm continues with selection, crossover, and mutation steps. Otherwise, it reverts back to the fast-non-dominated sorting, crowding distance assignment, and selection of suitable individuals before proceeding to the selection, crossover, and mutation steps. This iterative process continues for a specified number of generations. Finally, when the maximum number of generations is reached, the algorithm plots the Pareto frontier and employs an optimal point selection method to determine the final optimal solution [21].
Non-dominated sorting and the elitism preservation strategy contribute to the improved performance of the NSGA-II algorithm in terms of computational efficiency, preservation of good solutions, and maintaining population diversity, thereby enhancing its effectiveness in solving multi-objective optimization problems. With the NSGA-II, a three-dimensional Pareto front can effectively showcase the results of multi-objective optimization. It provides a clear visualization of the trade-off relationship between the three objectives and allows for a comprehensive understanding of the optimal solutions in the problem domain.

5. Case Studies

5.1. Parameters and Data Curation

This study validates the proposed scheduling method in a laboratory’s practical experimental system, of which, the electrolyzer is presented in Figure 5. Due to the limitations of the laboratory environment, the capacity of the equipment in the system is relatively small, within 100 kW. It should be noted that the optimized operation method proposed in this paper has little relationship with the equipment capacity level, and can also be applied to large-scale renewable energy systems. The rated power of the wind turbine, PV, electrolyzer, and compressor are 12 kW, 60 kW, 20 kW, and 5 kW, respectively. The storage capacity of the battery and hydrogen tank are 15 kWh and 5 kg. The cut-in wind speed of the wind turbine is 2.5 m/s, the cut-out wind speed is 25 m/s, and the rated wind speed is 12 m/s. The equipment parameters for the PV and wind system of this study are presented in Table 1. The other critical operation parameters, including unit operation and maintenance costs, are presented in Table 2, where the parameter values related to the equipment are derived from the real settings in the laboratory, and the parameter values related to the price are derived from [5].
Case studies were conducted in four scenarios, of which, the day-ahead source–load forecasts are shown in Figure 6, Figure 7 and Figure 8. As can be seen, the PV energy in Scenarios 1 and 2 is sufficient, while the wind power in Scenario 1 is much less than the others. In Scenarios 3 and 4, the PV energy is lacking, and there are non-steamed bun-shaped PV curves in Scenario 4, which are caused by the gloomy weather. The energy demand of the system performs basically the same in the four scenarios. Additionally, Figure 9a illustrates the electricity price for energy exchange between the system and the grid, while Figure 9b represents the hydrogen energy demand. In this study, it is assumed that the hydrogen energy output of the system is connected to an industrial park, resulting in a relatively stable hydrogen energy load.

5.2. Day-Ahead Optimization Results

The comparison between the proposed method and the conventional operational results is presented in Table 3. The operational strategy of the conventional method is as follows: When the power output from renewable energy sources exceeds the electricity demand, the excess power is primarily utilized by the electrolyzer to produce hydrogen. Once the hydrogen storage tank is full or the remaining power is insufficient to meet the minimum operating power of the electrolyzer, the remaining electricity is absorbed by the battery for storage and later use. It can be observed that—for the four selected typical days—the economic performance of the day-ahead optimization scheduling method proposed in this paper surpasses that of the non-optimized system, with around a 6% cost reduction rate.
These observations highlight the system’s adaptive strategies in response to electricity prices and renewable energy availability. Such strategies involve optimizing hydrogen production based on price differentials and leveraging energy storage capabilities to maximize system efficiency and cost-effectiveness.
The system operation results for the four scenarios are shown in Figure 10a–d. In the figure, “ E L Z ” refers to the electrolyzer power, “ G R D P ” represents the power purchased from the grid, “ G R D S ” represents the power sold to the grid, and “ S O C ” and “ L O H ” indicate the energy storage levels of the battery and hydrogen tank, respectively.
By examining the power exchange between the system and the grid, it can be observed that the system’s power purchase from the grid is considered as the power input (negative value), while the power sold to the grid is considered as the power output (positive value). Comparing the four scenarios, it is evident that almost all the power sold to the grid is concentrated between 11:00 and 15:00, as indicated by the red box in the figure. This is because the electricity selling price to the grid is the highest during this time period, and the PV power generation is also high, resulting in concentrated power selling during this timeframe.
The power sold to the grid is closely related to the system’s source–load balance. When renewable energy generation is low, as shown in Scenario 3, the system only purchases power from the grid and does not sell power to the grid. It is worth noting that the optimization model proposed in this paper takes into account the varying electricity prices. Therefore, the system attempts to avoid purchasing power from the grid during high-priced periods, such as 11:00–15:00 and 18:00–20:00 in Scenario 3. In these two time periods, the system compensates for the low input from the source and the high output from the load in two different ways. In the former time period, the system reduces the electrolyzer’s power and releases hydrogen from the storage tank to minimize energy consumption. In the latter time period, the system achieves energy balance by discharging the battery. These operational strategies contribute to the economic benefits of the system. Similar time periods can be observed in Scenarios 2 and 4, as indicated by the black box in the figure. Furthermore, it can be observed that the variation trend of the hydrogen storage level is generally opposite to the trend of electricity prices. This confirms that the optimization method proposed in this paper can effectively allocate renewable energy power by producing hydrogen during periods of low electricity prices and prioritizing meeting the electricity demand during periods of high electricity prices, thereby reducing system costs.
Finally, it can be seen that the variation rate of the battery’s SOC is relatively small. This also indicates that adopting the HESS can effectively reduce the required battery capacity, which is beneficial for reducing the system’s construction costs. If only a battery energy storage system is used, a significant increase in battery capacity would be required to achieve source–load balance.

5.3. Intra-Day Optimization Results

Figure 11 illustrates the intra-day optimization Pareto frontier for the first hour of a typical day in Scenario 1. The points on the Pareto frontier represent the optimal solution set of the multi-objective optimization model for day-ahead planning. At each point along this boundary, one objective cannot be further improved without sacrificing another objective. In other words, it represents the best trade-off or balance achieved among multiple objectives. The Pareto frontier showcases the range of feasible solutions that offer different trade-offs between objectives, allowing decision-makers to choose the most suitable solution based on their preferences and priorities.
The intra-day operation results of Scenario 1 are presented in Figure 12a. It is evident that the operation of the electrolyzer during the intra-day period closely follows the day-ahead scheduling plan. Furthermore, Figure 12b depicts the power exchange between the system and the grid, demonstrating good tracking performance.

5.4. Model Comparison

We compare the proposed dual-stage optimization model with the other existing models, of which, the optimization goals and solution methods are shown in Table 4. Model A represents a day-ahead model designed to minimize day-ahead costs as its optimization objective. The particle swarm optimization (PSO) algorithm is employed to solve this model. Meanwhile, Model B incorporates minute-level intra-day optimization and encompasses two objectives: operational profit and curtailment rate of renewable energy. Given its multi-objective nature, the study employs NSGA-II to tackle this optimization problem. In Model B, the original optimization objective was operational profit. However, in this paper, the consideration of revenue generated from hydrogen energy is omitted. As a result, the objective is transformed into the operational cost as a substitute measure to facilitate model comparison. Model C shares similarities with the proposed model in this paper, as both are dual-stage optimization models. In the day-ahead optimization, Model C focuses on minimizing operational costs but diverges in the intra-day optimization, where the objective is to achieve a tracking rate based on the day-ahead dispatch results, without explicitly considering the curtailment rate of renewable energy. Notably, the literature does not emphasize the two-stage solving approach for Model C, and this paper replaces it with the PSO algorithm.
Table 5 illustrates the daily operating costs, the deviation of exchanged electricity with the grid compared to the day-ahead schedule, and the actual renewable energy curtailment rate for four models, using Scenario 2 as an example.
From the optimization results of the four models, it can be observed that there is minimal difference in the daily operating results. This is because the daily operation is consistently considered by all the models. However, in terms of tracking the day-ahead operation results and managing renewable energy curtailment, Model D exhibits a more balanced and outstanding performance. This can be attributed to the fact that Model A focuses solely on day-ahead optimization and cannot handle the uncertainties arising from wind and solar forecasts. On the other hand, Model C’s intra-day optimization stage is a single-objective optimization that only considers tracking the day-ahead dispatch results, overlooking the issue of renewable energy curtailment. In contrast, both Model B and Model D adopt multi-objective optimization, providing a more balanced consideration of renewable energy curtailment. However, Model B is an intra-day optimization model and cannot provide operational guidance for the day-ahead operation. Therefore, Model D can be considered a more balanced and effective system operation optimization model among them.
By comparing the solution methods of the four models, it is evident that NSGA-II achieves good results for multi-objective optimization, while PSO is more favored for single-objective optimization due to its powerful search capabilities. The day-ahead optimization in this paper is formulated as a MILP model, which can be solved directly using solvers, resulting in faster computation. However, this approach sacrifices modeling accuracy, such as not considering the variability of electrolyzer efficiency with power. It is important to note that the proposed intra-day optimization stage in this paper helps alleviate the potential issues arising from the lower modeling accuracy to some extent.

6. Conclusions

This paper proposes a dual-stage optimization scheduling method for grid-connected systems with hybrid energy storage to balance renewable energy. The model consists of an economically optimal day-ahead scheduling model and an intra-day stage that uses accurate renewable energy forecasting to minimize curtailment. The proposed model offers several advantages, including enhanced economic efficiency in the day-ahead stage and multi-objective optimization in the intra-day stage. It provides a comprehensive approach to optimize the scheduling of grid-connected systems with hybrid energy storage, achieving a balance between different optimization objectives. Compared with other works, the proposed model balances economic and operational considerations, enabling the system to operate economically and stably. However, the faster solution speed of the day-ahead stage sacrifices some modeling accuracy, making it suitable for scenarios where high modeling accuracy is not critical.
Overall, the proposed dual-stage optimization model offers a promising solution for achieving economic efficiency and minimizing curtailment in grid-connected systems with renewable energy sources, contributing to the effective integration of renewable energy and the sustainable operation of the power grid.
Refining electrolyzer modeling and exploring faster and more effective optimization algorithms are promising directions for future research. Improving electrolyzer modeling accuracy will better capture system dynamics, while advanced algorithms will enhance the efficiency and effectiveness of the optimization process, enabling more comprehensive and optimal system operations. These advancements will contribute to the development of sophisticated and efficient energy management strategies in the future.

Author Contributions

Conceptualization and methodology, D.L.; software and writing, J.S.; resources, Y.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 51877188.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

Author Di Lu was employed by the company Powerchina Huadong Engineering Corporation. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Nomenclature

HESShybrid energy storage system
PVphotovoltaic
PSOparticle swarm optimization
MILPmixed-integer linear programming
NSGA-IInon-dominated sorting genetic algorithm II
CNN-LSTMconvolutional neural network–long short-term memory
STCstandard test condition
NOCTnominal operating cell temperature
p p v solar power output
p w i n d wind power output
p l o a d electricity load
p e l z electrolyzer power
p c m p compressor power
p b c h charging power of battery
p b d i s discharging power of battery
p g e x exchange power with the utility grid
p r e p a l renewable power curtailment
p m o u t hydrogen output of system
m e l z o u t hydrogen generation amount of electrolyzer
m h o u t hydrogen demand
S o C stage of charge of battery
L o H level of hydrogen of storage tank
p ˜ r e p a l intra-day renewable power curtailment
p ˜ g e x intra-day exchange power with the utility grid
p ˜ e l z intra-day electrolyzer power
p ˜ b intra-day charging/discharging power of battery
O b j d h objective function of day-ahead stage
O b j i d objective function of intra-day stage
C e p u r electricity purchasing costs
C o m operation and maintenance costs
C b d e battery degradation costs
C r e p a l penalty costs of renewable power
B e s e l l revenue from selling electricity
B h s e l l revenue from selling hydrogen
p ˜ g e x power exchange with the grid: day-ahead optimization result
p ˜ r e p a l renewable power curtailment: day-ahead optimization result
p ˜ b charging or discharging power of the battery: day-ahead optimization result
p ˜ e l z power of electrolyzer: day-ahead optimization result
p w i n d r a t e d rated power of wind turbine
p p v r a t e d rated power of PV
γ b dissipation rate of the battery
γ h dissipation rate of the hydrogen tank
η b c h charging efficiency of the battery
η b d i s discharging efficiency of the battery
Q b R a t e d rated capacity of the battery
η e l z electrolytic efficiency of the electrolyzer
L H V H 2 low heat value of the hydrogen
η h , i n input efficiency of the hydrogen tank
η h , o u t output efficiency of the hydrogen tank
Q h R a t e d rated capacity of the hydrogen tank
c e p u r / e s e l l unit electricity purchasing/selling prices
o e l z / c m p / p v / w i n d unit operation costs of electrolyzer/compressor/PV/wind turbine
d b unit degradation cost of the battery
c r e p a l unit cost of renewable power curtailment.
c h s e l l unit hydrogen selling price
β the ratio of compressor power consumption to compressed hydrogen energy
p e m a x / m i n l z the maximum/minimum power of electrolyzer
p c m a x / m i n m p the maximum/minimum power of compressor
p b c h m a x / m i n the maximum/minimum charging power of the battery
p b d i s m a x / m i n the maximum/minimum discharging power of the battery
p g e x m a x / m i n the maximum/minimum exchange power with the grid
ϵ r e c u l the upper limit of the renewable power curtailment rate

References

  1. Ikonnikova, S.A.; Scanlon, B.R.; Berdysheva, S.A. A global energy system perspective on hydrogen Trade: A framework for the market color and the size analysis. Appl. Energy 2023, 330, 120267. [Google Scholar] [CrossRef]
  2. Salama, H.S.; Magdy, G.; Bakeer, A.; Vokony, I. Adaptive coordination control strategy of renewable energy sources, hydrogen production unit, and fuel cell for frequency regulation of a hybrid distributed power system. Prot. Control Mod. Power Syst. 2022, 7, 34. [Google Scholar] [CrossRef]
  3. Khan, M.I.; Asfand, F.; Al-Ghamdi, S.G. Progress in research and technological advancements of thermal energy storage systems for concentrated solar power. J. Energy Storage 2022, 55, 105860. [Google Scholar] [CrossRef]
  4. De Abreu, V.H.S.; Pereira, V.G.F.; Proença, L.F.C.; Toniolo, F.S.; Santos, A.S. A Systematic Study on Techno-Economic Evaluation of Hydrogen Production. Energies 2023, 16, 6542. [Google Scholar] [CrossRef]
  5. Abomazid, A.M.; El-Taweel, N.A.; Farag, H.E.Z. Optimal Energy Management of Hydrogen Energy Facility Using Integrated Battery Energy Storage and Solar Photovoltaic Systems. IEEE Trans. Sustain. Energy 2022, 13, 1457–1468. [Google Scholar] [CrossRef]
  6. Lee, M.; Ham, S.H.; Lee, S.; Kim, J.; Kim, Y. Multi-objective optimization of solar-assisted ground-source heat pumps for minimizing life-cycle cost and climate performance in heating-dominated regions. Energy 2023, 270, 126868. [Google Scholar] [CrossRef]
  7. Zheng, Y.; You, S.; Bindner, H.W.; Münster, M. Optimal day-ahead dispatch of an alkaline electrolyser system concerning thermal–electric properties and state-transitional dynamics. Appl. Energy 2022, 307, 118091. [Google Scholar] [CrossRef]
  8. Sun, Y.; Zhang, B.; Ge, L.; Sidorov, D.; Wang, J.; Xu, Z. Day-ahead optimization schedule for gas-electric integrated energy system based on second-order cone programming. CSEE J. Power Energy Syst. 2020, 6, 142–151. [Google Scholar] [CrossRef]
  9. Hong, Z.; Wei, Z.; Han, X. Optimization scheduling control strategy of wind-hydrogen system considering hydrogen production efficiency. J. Energy Storage 2022, 47, 103609. [Google Scholar] [CrossRef]
  10. Amin, M.; Shah, H.H.; Bashir, B.; Iqbal, M.A.; Shah, U.H.; Ali, M.U. Environmental Assessment of Hydrogen Utilization in Various Applications and Alternative Renewable Sources for Hydrogen Production: A Review. Energies 2023, 16, 4348. [Google Scholar] [CrossRef]
  11. Wang, X.; Han, L.; Wang, C.; Yu, H.; Yu, X. A time-scale adaptive dispatching strategy considering the matching of time characteristics and dispatching periods of the integrated energy system. Energy 2023, 267, 126584. [Google Scholar] [CrossRef]
  12. Yuan, Z.P.; Xia, J.; Li, P. Two-Time-Scale Energy Management for Microgrids With Data-Based Day-Ahead Distributionally Robust Chance-Constrained Scheduling. IEEE Trans. Smart Grid 2021, 12, 4778–4787. [Google Scholar] [CrossRef]
  13. Fu, P.; Pudjianto, D.; Zhang, X.; Strbac, G. Integration of Hydrogen into Multi-Energy Systems Optimisation. Energies 2020, 13, 1606. [Google Scholar] [CrossRef]
  14. Liponi, A.; Frate, G.F.; Baccioli, A.; Ferrari, L.; Desideri, U. Impact of wind speed distribution and management strategy on hydrogen production from wind energy. Energy 2022, 256, 124636. [Google Scholar] [CrossRef]
  15. Sharafi, M. Multi-objective optimal design of hybrid renewable energy systems using PSO-simulation based approach. Renew. Energy 2014, 68, 67–79. [Google Scholar] [CrossRef]
  16. Seyyedabbasi, A.; Kiani, F. Sand Cat swarm optimization: A nature-inspired algorithm to solve global optimization problems. Eng. Comput. 2023, 39, 2627–2651. [Google Scholar] [CrossRef]
  17. Luo, G.; Mei, Y. Multi-Objective Optimization Scheduling of Microgrids Based on Particle Swarm Optimization Algorithm. In Proceedings of the 2023 International Conference on Network, Multimedia and Information Technology (NMITCON), Bengaluru, India, 1–2 September 2023; pp. 1–7. [Google Scholar] [CrossRef]
  18. Kunakote, T.; Sabangban, N.; Kumar, S.; Tejani, G.G.; Panagant, N.; Pholdee, N.; Bureerat, S.; Yildiz, A.R. Comparative Performance of Twelve Metaheuristics for Wind Farm Layout Optimisation. Arch. Comput. Methods Eng. 2022, 29, 717–730. [Google Scholar] [CrossRef]
  19. Anosri, S.; Panagant, N.; Bureerat, S.; Pholdee, N. Success history based adaptive multi-objective differential evolution variants with an interval scheme for solving simultaneous topology, shape and sizing truss reliability optimisation. Knowl.-Based Syst. 2022, 253, 109533. [Google Scholar] [CrossRef]
  20. Rejeb, O.; Alirahmi, S.M.; Assareh, E.; El Haj Assad, M.; Jemni, A.; Bettayeb, M.; Ghenai, C. Innovative integrated solar powered polygeneration system for green Hydrogen, Oxygen, electricity and heat production. Energy Convers. Manag. 2022, 269, 116073. [Google Scholar] [CrossRef]
  21. Cao, Y. Predication of the sensitivity of a novel daily triple-periodic solar-based electricity/hydrogen cogeneration system with storage units: Dual parametric analysis and NSGA-II optimization. Renew. Energy 2022, 192, 340–360. [Google Scholar] [CrossRef]
  22. Maheri, A.; Unsal, I.; Mahian, O. Multiobjective optimisation of hybrid wind-PV-battery-fuel cell-electrolyser-diesel systems: An integrated configuration-size formulation approach. Energy 2022, 241, 122825. [Google Scholar] [CrossRef]
  23. Wu, Q.; Li, C. Modeling and operation optimization of hydrogen-based integrated energy system with refined power-to-gas and carbon-capture-storage technologies under carbon trading. Energy 2023, 270, 126832. [Google Scholar] [CrossRef]
  24. Dincer, I.; Ishaq, H. Renewable Hydrogen Production; Elsevier: Amsterdam, The Netherlands, 2023. [Google Scholar]
  25. Ghenai, C.; Salameh, T.; Merabet, A. Technico-economic analysis of off grid solar PV/Fuel cell energy system for residential community in desert region. Int. J. Hydrogen Energy 2020, 45, 11460–11470. [Google Scholar] [CrossRef]
  26. Wang, D.; Gan, J.; Mao, J.; Chen, F.; Yu, L. Forecasting power demand in China with a CNN-LSTM model including multimodal information. Energy 2023, 263, 126012. [Google Scholar] [CrossRef]
  27. Zhang, Y.; Qin, C.; Srivastava, A.K.; Jin, C.; Sharma, R.K. Data-Driven Day-Ahead PV Estimation Using Autoencoder-LSTM and Persistence Model. IEEE Trans. Ind. Appl. 2020, 56, 7185–7192. [Google Scholar] [CrossRef]
  28. Zhang, H.; Wu, Q.; Chen, J.; Lu, L.; Zhang, J.; Zhang, S. Multiple stage stochastic planning of integrated electricity and gas system based on distributed approximate dynamic programming. Energy 2023, 270, 126892. [Google Scholar] [CrossRef]
  29. Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A fast and elitist multiobjective genetic algorithm: NSGA-II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef]
  30. Niu, M.; Li, X.; Sun, C.; Xiu, X.; Wang, Y.; Hu, M.; Dong, H. Operation Optimization of Wind/Battery Storage/Alkaline Electrolyzer System Considering Dynamic Hydrogen Production Efficiency. Energies 2023, 16, 6132. [Google Scholar] [CrossRef]
  31. Wang, J.; Li, D.; Lv, X.; Meng, X.; Zhang, J.; Ma, T.; Pei, W.; Xiao, H. Two-Stage Energy Management Strategies of Sustainable Wind-PV-Hydrogen-Storage Microgrid Based on Receding Horizon Optimization. Energies 2022, 15, 2861. [Google Scholar] [CrossRef]
Figure 1. Renewable energy plants augmented with hybrid energy storage systems.
Figure 1. Renewable energy plants augmented with hybrid energy storage systems.
Energies 17 00737 g001
Figure 2. The flowchart of the proposed dual-stage optimization scheduling model.
Figure 2. The flowchart of the proposed dual-stage optimization scheduling model.
Energies 17 00737 g002
Figure 3. Diagram of the dual-stage power allocation model. The numbers of day-ahead scheduling denote the hours, while the numbers of intra-day scheduling denote the minutes.
Figure 3. Diagram of the dual-stage power allocation model. The numbers of day-ahead scheduling denote the hours, while the numbers of intra-day scheduling denote the minutes.
Energies 17 00737 g003
Figure 4. The solving process with NSGA-II.
Figure 4. The solving process with NSGA-II.
Energies 17 00737 g004
Figure 5. Alkaline water electrolyzer of the system in the laboratory.
Figure 5. Alkaline water electrolyzer of the system in the laboratory.
Energies 17 00737 g005
Figure 6. Wind power forecasting results.
Figure 6. Wind power forecasting results.
Energies 17 00737 g006
Figure 7. Solar power forecasting results.
Figure 7. Solar power forecasting results.
Energies 17 00737 g007
Figure 8. Electricity load forecasting results.
Figure 8. Electricity load forecasting results.
Energies 17 00737 g008
Figure 9. Data input of the renewable energy system.
Figure 9. Data input of the renewable energy system.
Energies 17 00737 g009
Figure 10. Day-ahead optimization results of different scenarios.
Figure 10. Day-ahead optimization results of different scenarios.
Energies 17 00737 g010
Figure 11. Pareto frontier of the first hour in Scenario 1.
Figure 11. Pareto frontier of the first hour in Scenario 1.
Energies 17 00737 g011
Figure 12. Intra-day optimization results of Scenario 1.
Figure 12. Intra-day optimization results of Scenario 1.
Energies 17 00737 g012
Table 1. Parameters of the PV and wind turbines.
Table 1. Parameters of the PV and wind turbines.
ParametersValuesParametersValues
p w i n d r a t e d 12 kW p p v r a t e d 60 kW
v i n 2.5 m/s v o u t 25 m/s
v r a t e d 12 m/s G S T C 1 kW/ m 2
T c , N O C T 20 °C T r e f 25 °C
G N O C T 0.8 kW/ m 2 α p 0.005 °C−1
T a , N O C T 25 °C f P V 80%
Table 2. Other operation parameters of the system.
Table 2. Other operation parameters of the system.
ParametersValuesParametersValues
p e l z m a x  (kW)20 o p v / o w i n d  ($/kW)0.033
η b i n / η b o u t 0.975 o e l z 0.29
p e l z m i n  (kW)5 c r e p a l  ($/kW)0.05
p c m p m a x  (kW)5 o H s  ($/kg)0.017
p b c h m a x  (kW)30 d b  ($/kW)0.020
L H V H 2  kJ/mol242 o C m p  ($/kW)0.012
β 2 γ b / γ h 0.0006%
η e l z 0.65 Q b R a t e d  (kWh)15
Q h R a t e d  (kg)5 η h i n / η h o u t 0.975
Table 3. Day-ahead optimization results of four typical days and the comparison with the conventional operation.
Table 3. Day-ahead optimization results of four typical days and the comparison with the conventional operation.
Scenario 1Scenario 2Scenario 3Scenario 4
Daily costs ($)Conventional operation71.3446.7374.5260.78
Optimized operation66.8244.1070.1056.46
Costs reduction rate (%)6.345.635.937.11
Table 4. Optimization goals and solution methods of the four models.
Table 4. Optimization goals and solution methods of the four models.
Optimization GoalsSolution Methods
Day-AheadIntra-DayDay-AheadIntra-Day
Model A [7]Daily operation costs/Particle swarm algorithm/
Model B [30]/Costs & Power abandonment rate/NSGA-II
Model C [31]Daily operation costsTrace the day-ahead operation schemesParticle swarm algorithmParticle swarm algorithm
Model D (Proposed model)Daily operation costsTrace the day-ahead operation schemes & Power abandonment rateMILPNSGA-II
Table 5. Optimization results of the four models.
Table 5. Optimization results of the four models.
Model AModel BModel CModel D
Daily operating costs ($)45.246.145.344.1
Deviation of exchanged electricity (%)5.3/1.42.1
Renewable energy curtailment rate (%)6.41.77.61.3
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

Lu, D.; Peng, Y.; Sun, J. Dual-Stage Optimization Scheduling Model for a Grid-Connected Renewable Energy System with Hybrid Energy Storage. Energies 2024, 17, 737. https://0-doi-org.brum.beds.ac.uk/10.3390/en17030737

AMA Style

Lu D, Peng Y, Sun J. Dual-Stage Optimization Scheduling Model for a Grid-Connected Renewable Energy System with Hybrid Energy Storage. Energies. 2024; 17(3):737. https://0-doi-org.brum.beds.ac.uk/10.3390/en17030737

Chicago/Turabian Style

Lu, Di, Yonggang Peng, and Jing Sun. 2024. "Dual-Stage Optimization Scheduling Model for a Grid-Connected Renewable Energy System with Hybrid Energy Storage" Energies 17, no. 3: 737. https://0-doi-org.brum.beds.ac.uk/10.3390/en17030737

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