Abstract
The heat and hydrogen balance of the hydrogen energy storage system’s intermittent operation becomes a key factor affecting the performance of the wind-hydrogen hybrid system (W-HHS). This work designed a hydrogen energy storage system (HESS), including waste heat utilization. Then, a dual state of charge (SOC) model is established, in which hydrogen and heat storage is considered. Furthermore, based on the distributionally robust method, an optimal dispatch method of W-HHS is proposed to reduce the operation cost of conventional units in the grid and increase the revenue of the W-HHS. The previously proposed dual SOC model of heat-hydrogen balance is regarded as a constraint in this cooperative dispatch. The effectiveness and efficiency of the dual SOC model were verified on the IEEE 30-bus system with an actual wind plant data set. The results show that the hydrogen-heat dual SOC model can fully reflect the influence of heat and hydrogen balance on the operation of the W-HHS. The cooperative dispatch method improves the reliability of the W-HHS operation under the premise of ensuring the heat-hydrogen balance. When the constraints of hydrogen balance SOC and heat balance SOC are met simultaneously, the available power of the wind plant is 6–8% lower than the ideal situation. Parameter analysis indicates that reducing the heat dissipation coefficient can reduce the influence of the SOC constraint of heat balance on the dispatch strategy and increase the power output of the wind plant. When the heat dissipation coefficient is less than 1/1,200, the heat balance SOC constraint fails.
Introduction
With the change of energy structure, a new power system with a high proportion of renewable energy will become the mainstream for the development of energy system transformation in the future (Zhang and Chen, 2020). Wind and solar have become the primary power sources (Jiayu et al., 2021). However, the wide application of renewable energy raises many problems. The ensuing stability and security issues have become a critical bottleneck restricting the development of a high-proportion renewable power grid (Impram et al., 2020).
Wind power is one of the primary forms of renewable energy, and its participation in the power system dispatch has received extensive attention. Some researchers use stochastic optimization (Zhu et al., 2020), robust optimization (Yu et al., 2020), distributionally robust optimization (Guo et al., 2020), and other optimization methods to try to describe the existing uncertainties in the cooperative dispatch integrated with wind power. Among them, the distributionally robust method is a data-driven optimization method, which constructs the ambiguity set of uncertain parameter probability by using the information implicit in the actual data. Significantly, the distributionally robust method based on waister distribution has been developed rapidly. Ref. (Guo et al., 2020) compared and analyzed the adaptability of distributionally robust optimization methods based on Wasserstein divergence distribution. (Esfahani and Kuhn, 2018) demonstrated this method’s performance guarantee and ease of treatment in detail. At the same time, it should also be pointed out that according to lemma 8 (HOTA et al., 2019), the premise of transforming the distributionally robust model of uncertainty into a set of CVar constraints is that the optimization problem studied is convex, which limits the application of this method to a certain extent. In (Yang et al., 2020), a distributionally robust chance constraint (DRCC) model is proposed for the optimal power gas flow (OPGF) problem with uncertain wind power. The DRCC-OPGF model is reformulated as a treatable mixed-integer convex programming problem.
Besides, some works focus on the different forms of hybrid systems with a wind plant. Different types of energy storage were considered to suppress the wind fluctuation, such as batteries (Xu et al., 2020), compressed air energy storage (Alirahmi et al., 2021), hot dry rock geothermal energy (Si et al., 2021), and hydrogen energy storage system (HESS) (Xiao et al., 2020), etc. constitute the optimal dispatch of the hybrid system. The cooperative dispatch strategy of these different hybrid systems was also well studied.
At present, the hybrid system formed by wind plants and energy storage has become an important form to realize the stable operation with a high proportion of wind power. On this basis, the research about reducing the operating costs of conventional units, increasing the power output of wind plants, and raising the operating efficiency of energy storage systems has become more and more popular. HESS has the advantages of fast response speed, large energy storage capacity, and cross-season energy storage. Also, it can adapt to the operating environment of cold climates and significant temperature differences between day and night (Petkov and Gabrielli, 2020). The application of HESS on the wind-storage hybrid system has attracted scholars’ attention in recent years (Li et al., 2020).
In the existing works, some aim to suppress the uncertainty of wind power and improve the system economy. Ref. (Zhang and Wan, 2014) established a hydrogen energy storage model to reduce power curtailment and gave a scheme for wind/hydrogen production through water electrolysis. Ref. (Qiu et al., 2020) proposed a two-stage robust optimization method with stability constraints considering the HESS’s safety and small disturbance stability and dynamic response characteristics. The power output range of the HESS was then determined. Ref. (Xu et al., 2020) further proposed a distributionally robust chance-constrained (DRCC) dispatch method for W-HHS in the day-ahead power market and transformed it into a standard second-order cone programming problem for a solution.
Other scholars regard the electrolyzer as an adjustable load and research from the perspective of demand response. Ref. (Mirzaei et al., 2018) proposed a safety constrained unit commitment model with high-proportion wind power in coordination with demand response. The HESS was used in this work to reduce energy consumption and improve system reliability. Ref. (Mansour-Saatloo et al., 2020) used the HESS in the energy hub to construct a robust dispatch method, which considered the demand response of the HESS and the characteristics of cogeneration.
The existing research shows that in the operation of the HESS, cogeneration can happen both in the hydrogen production and fuel cell generation stages. Its thermal characteristics have an important influence on the reliable and efficient operation of the HESS. On this basis, some work studied the future development of hybrid systems with ultra-high penetration of renewable energy. A thermoelectric hydrogen model with startup/shutdown constraints and a new seasonal hydrogen energy storage model was proposed (Wen et al., 2020). Ref. (Pan et al., 2020) studied the influence of heat load increase on the operation cost of the hybrid system and proposed a day-ahead dispatch scheme for the power system with high-proportion wind power and HESS. Ref. (Wei et al., 2021) further established a heat-hydrogen efficiency model for the intermittent operation of the electrolyzer by describing the heat exchange process in detail. Ref. (Kovač et al., 2021) carried out research on the thermal management of distributed photovoltaic-hydrogenation stations, put forward the optimized thermal management strategy, and improved the energy utilization efficiency. In order to ensure the energy efficiency of the renewable hydrogen energy system, the thermal management of the metal hydride tank is studied (Endo et al., 2021). Research shows that heat management can improve energy efficiency by 50%.
It can be seen that in the grid-connected operation of the hybrid system with renewable energy, the HESS can use the electrolyzer to transform the excess wind power into hydrogen to provide a down-reserves. At the same time. The fuel cell equipped in the HESS can consume hydrogen to generate power to provide an up-reserves for the hybrid system. The electrolyzer and the fuel cell work alternately and intermittently to suppress the power fluctuation. Notice that the electrolyzer and fuel cell’s operation efficiency and response speed are closely related to their operating temperature. Therefore, the ability of HESS to respond to wind power fluctuations quickly, long-term, and efficiently is restricted by the following two factors. First, the balance between hydrogen production in the electrolyzer and the hydrogen consumption of the fuel cell should always be satisfied. Second, the balance between heat production under operation conditions and heat consumption under reserve mode should be held as well. Thus, the operation of W-HHS should meet both heat and hydrogen production-consumption balance simultaneously when it is connected to the grid.
Existing works have neglected the heat balance between electrolyzers and fuel cells in the reserve mode and the inherent relationship between hydrogen and heat balance. While maintaining heat balance and hydrogen balance is the key to ensuring the economic operation of the hybrid power equipped with HESS. Therefore, based on the operating characteristics of the HESS, this paper first designs a brand new structure of the HESS, including the waste heat utilization system, and gives the heat-hydrogen state of charge (SOC) models. Furthermore, the distributionally robust optimization approach is used to model the uncertainty of wind plant power output. Then, an optimal dispatch method for the W-HHS considering the heat and hydrogen balance is proposed. The models and method are verified with the IEEE 30-bus system. The actual data set of a wind plant in Qinghai Province is used in this simulation. Finally, the key parameters are analyzed, and their impacts on system performance are well studied.
The rest of this article is organized as follows. Mathematical Model of HESS With Waste Heat Utilization System elaborates the overall design scheme and mathematical model of the HESS with waste heat utilization system. Mathematical Models of the Wind Plant introduces the mathematical models of wind plants. The cooperative dispatch method of W-HHS considering dual SOC constraints is proposed in Cooperative Dispatch Method of W-HHS Considering Dual SOC Constraints. Case Study verifies the effectiveness of the proposed method through a case consisting of actual data, followed by the disscusions in Discussions.
Mathematical Model of HESS With Waste Heat Utilization System
Structure of HESS With Waste Heat Utilization System
The structure of a HESS with the waste heat utilization system is shown in Figure 1. The whole system consists of an electrolyzer, fuel cell, hydrogen storage system, and waste heat utilization system. This HESS adopts the alkaline electrolyzer as the water electrolysis device due to its mature technology and fast response speed. The electrolyzer is used to absorb the excess wind power to produce hydrogen. As for the fuel cell, the proton-exchange membrane fuel cell (PEMFC) is chosen to consume hydrogen and generate electricity considering the requirement of the operating temperature match.
FIGURE 1
Diagram of the wind-hydrogen hybrid system.
The hydrogen storage system includes a hydrogen tank, an air/oxygen system, and a water supply cycle. The waste heat utilization system consists of a thermal tank, a heat storage exchanger, and a heat cycle system. The heat cycle system builds a heating/cooling cycle between the electrolyzer, PEMFC, water supply tank, and heat storage tank to stabilize the system temperature. The mathematical models of the various components of the HESS are given in the following sections.
Model of HESS With Waste Heat Utilization System
Model of Alkaline Electrolyzer and PEMFC
The alkaline electrolyzer has the advantages of short response time and mature technology. It is very suitable to be used in a HESS with wind plants. Its simplified model (Clua et al., 2018; Xiong et al., 2021) can be expressed as:where and , respectively, represent the electricity consumed and the heat produced by the electrolyzer at time . represents the rate of hydrogen production, and represents the higher heating value of hydrogen. represents the efficiency of the electrolyzer.
The operating temperature of the PEMFC should match that of the alkaline electrolyzer to facilitate the recovery and utilization of the waste heat of the HESS. Its model can be expressed as:where and represent the electricity and heat produced by the fuel cell at time , respectively. represents the hydrogen consumption rate of PEMFC. represents the efficiency of PEMFC.
Dual SOC Model for Hydrogen and Heat Balance
The HESS stores the hydrogen produced by the electrolyzer in the hydrogen tank and supplies the hydrogen and oxygen to the PEMFC when electricity is needed. At the same time, water is supplied to the electrolyzer, and the electricity generation products of the PEMFC are collected through the water supply cycle. Taking the total energy stored in the hydrogen tank as the parameter to describe the hydrogen energy storage state, we establish SOC model of the hydrogen balance as:where represents the SOC of hydrogen energy stored in the tank at time . represents the maximum capacity of the hydrogen tank, and represents the time interval.
The waste heat utilization system provides thermal energy for each link and collects waste heat. On the one hand, the waste heat is used to preheat the water entering the electrolyzer, and on the other hand, it is used to maintain the rated temperature of the PEMFC in reserve mode. Also, when the PEMFC is in operation mode, the waste heat needs to provide heat to the electrolyzer and recover the heat generated from the PEMFC. After meeting the heat balance requirements of the electrolyzer and PEMFC, the waste heat utilization system can also use the remaining thermal energy to provide heat for extra heat loads in the hybrid system. The SOC model of heat balance is:where represents the SOC of thermal energy stored in the thermal tank at time . represents the maximum thermal energy storage capacity of the thermal tank. and represent the heat generated by the electrolyzer and the PEMFC, respectively. and represent the heat consumed by the HESS and the heat provided for the extra heat loads, respectively. represents the efficiency of the heat exchanger.
In Eq. 4, is composed of three losses: the heat dissipation of the electrolyzer, the heat dissipation of the fuel cell, and the preheating of the water supply cycle. Thus, can be expressed as:where and represent the operating temperature of the electrolyzer and PEMFC. represent the ambient temperature. and are the installed capacity of electrolyzer and PEMFC. , , and represent the molar mass flow rate, specific heat capacity and temperature of the water entering the electrolyzer. and represent the heat dissipation coefficient of the unit capacity electrolyzer and PEMFC (Hwang, 2005), which are defined as the heat dissipation area per unit capacity to unit area thermal resistance . and can be expressed as:
Mathematical Models of the Wind Plant
Models of Wind Plant Output Power
The power output of a wind plant can be modeled with wind speed and installed power (Guo et al., 2020). By further equating the wind speed change as a power output impact factor , a simplified model of the wind plant output can be expressed aswhere represents the installed power of the wind plant. represents the actual power output of the wind plant.
While, in an actual situation, the system operator usually uses the predicted wind plant output to determine the dispatch decisions and the dispatchable power. The predicted power output and dispatchable power can be expressed as:where represents the predicted power output impact factor. and represent the predicted output and dispatchable power of the wind plant. reflects the availability of power grid dispatch with wind power forecasting. represents the dispatchable power, forming the dispatching curve of the wind plant, which is the grid-connected power guaranteed by the wind plant. means availability factor of wind plant. Based on this, we can further define the dispatch tracking deviation index of the wind plant (Guo et al., 2020), which is:where represents the total power deviation in the whole period of dispatch.
A Distributionally Robust Model of Wind Plant Output Uncertainty
Distributionally robust optimization is a data-driven analysis approach and has been widely studied and applied. This method does not require an accurate probability distribution to characterize the uncertainties. Its conservativeness lies between robust optimization and stochastic optimization. Therefore, the distributionally robust optimization approach is applied in this work. In this method, we adopt the Wasserstein divergence (Guo et al., 2020) as the ambiguous set to measure the uncertainties of wind power output. Hence, we have:where represents the probability distribution of the actual power output of the wind plant. represents the empirical distribution of the power output. represents the set formed by all possible distributions satisfying Wasserstein divergence , where is the radius of the ambiguous set. The calculation method of the parameters mentioned above can be found in (Esfahani and Kuhn, 2018).
Cooperative Dispatch Method of W-HHS Considering Dual SOC Constraints
The optimal dispatch of the W-HHS should reduce the operation costs of conventional units in the grid as much as possible and then increase the grid-connected power generated by wind plants to make full use of clean energy and reduce carbon emissions. Hence, we try to optimize the dispatch strategy of the HESS to ensure that the power output of the W-HHS can track the dispatchable power accurately, thereby reducing the impact of wind power fluctuation on the power grid.
Objective Function
In the cooperative dispatch, historical wind power data is commonly used to predict the wind plant output. The model that we establish combines the dispatchable power of wind plants with the availability coefficients to optimize the output of conventional units and increase the grid-connected output of wind plants. The operation cost of conventional units and carbon emissions can also be reduced. Therefore, the goal of the W-HHS’s optimal dispatch is:where represents the probability distribution of the actual wind plant output. represents the on-grid power tariff. represents the penalty coefficient. represents the power output of the units. , , and represent the operation cost coefficients of the conventional unit. represents the bus index. represents the set of all buses.
Constrains
Dual SOC Constrains
In order to ensure that the HESS can dispatch the electrolyzer and PEMFC to provide reserves for the hybrid system, the optimal dispatch strategy must meet the balance of the hydrogen production (generated by electrolyzer) and the hydrogen consumption (consumed by the PEMFC). Also, the strategy should hold the heat balance between heat production and consumption.
The hydrogen balance requires that the system does not need to supplement hydrogen from outside during the operation, and the hydrogen in the storage tank will never be lower than the lower bound of the rated value after operating for a period . Thus, we have the hydrogen balance SOC constraint is:where is the initial value of the hydrogen tank.
Similarly, the heat balance requires that the waste heat utilization system supports the W-HHS to operate continuously and efficiently within the specified operation period without the input of external thermal energy. Its SOC constraint can be expressed as:where is the initial value of the thermal tank.
Besides of heat SOC constraint, heat balance also includes heat power balance constraint in operation.where represents the heat stored in the thermal tank.
Constraints for Grid
The constraints of power grid dispatch include power flow constraints, line capacity constraints, and power balance constraints. In this section, the linear AC power flow model is adopted. Then, we have the power flow constraints as follows:where and represent the active and reactive power flow on line at time . and represent the voltage magnitude and phase angle of bus . and represent the active and reactive power flow errors on line (Guo et al., 2020).
The constraint of line capacity is:where represents the capacity of line . Eq. 17 can be linearized by the outer approximation approach as
The power balance constraint iswhere and represent the active and reactive power injected into bus , that is, the power generated by the generator connected to the bus minus the active/reactive load / on the same bus. and are the real and imaginary parts of elements in the network admittance matrix.
Model Solution
The problem established in Eqs. 12–19 is a max-min optimization problem. The actual power output of the wind plant is modeled by the uncertain ambiguous set given in Eq. 11, which cannot be solved directly. To this end, we adopt the method proposed in (Guo et al., 2020) to transform the set into a set of linear chance constraints, which are:where is the number of data samples. is the confidence level. and are the actual wind plant output and dispatchable grid-connected power represented by samples. and are corresponding dual variables.
At this time, since we use samples to describe the uncertainties of the wind plants outputs, the power balance constraint of each line in Eq. 19 can be rewritten as:
The lower and upper bounds of each variable restricted by the limitations of technology and environment are all considered in modeling, which is not listed due to the space limitation. Besides, we notice that product terms of two decision variables exist in Eq. 5, making the model nonlinear. This term can be piecewise linearized by the Boolean expansion method (Pereira et al., 2005).
So far, the W-HHS proposed cooperative dispatch method is transformed into a mixed-integer linear program, which can be solved with matlab 2016b and Cplex12.6.
Case Study
System Parameters
In this section, we test our models and method on the IEEE 30-bus testbed. The actual historical data obtained from the Qinghai Province’s 90 MW wind plant is used as samples. The parameters of the system are shown in Table 1. The W-HHS is integrated on bus #12, which is also given in Figure 2. The capacity of HESS is selected according to the local energy storage configuration policy, equal to 10% of the wind plant capacity. The specific historical data is selected as the predicted value of the wind plant output. 20% prediction error is set according to the current wind power prediction level.
TABLE 1
| Parameters | Value |
|---|---|
| Wind power installed capacity MW | 90 |
| PEMFC installed capacity MW | 9 |
| Electrolyzer installed capacity MW | 9 |
| Volume of hydrogen storage tank m3 | 6.41 |
| Pressure of hydrogen storage tank MPa | 10 |
| Efficiency of PEMFC % | 60 |
| Efficiency of electrolyzer % | 60 |
| PEMFC operating temperature °C | 80 |
| Electrolyzer operating temperature °C | 80 |
| Heat exchanger efficiency % | 80 |
| Ambient temperature °C | 20 |
| Higher heating value of hydrogen kJ/mol | 282 |
| Heat dissipation coefficient of PEMFC °C−1 | 1/800 |
| Heat dissipation coefficient of electrolyzer °C−1 | 1/800 |
| Water specific heat capacity kJ/(kg⋅°C) | 4.2 |
| Availability coefficient | 0.8–1 |
| Confidence level | 0.05 |
System parameters.
FIGURE 2
Diagram of the IEEE 30 power system.
The Latin hypercube method generates 100 samples (Shu et al., 2014) to construct the proposed distributionally robust optimization model. The upper and lower bounds of predicted output, available power, and uncertainties of wind plants are shown in Figure 3.
FIGURE 3
Power chart of the wind plant.
Cooperative Dispatch Results of W-HHS
To analyze the influence of heat balance and hydrogen balance on the dispatch strategy of HPS, we study four different scenarios, which are listed as follows. The results are shown in Table 2.
TABLE 2
| Scenarios | Operation cost of units ($×106) | Wind plant profits ($×106) | Residual hydrogen (kg) | Residual thermal energy (kW⋅h) | Wind plant availability coefficient |
|---|---|---|---|---|---|
| Case 1 | 0.14223 | 0.21934 | −25.04 | — | 1 |
| Case 2 | 0.14645 | 0.20951 | −19.86 | 126.44 | 0.9552 |
| Case 3 | 0.14564 | 0.20575 | 22.58 | −489.15 | 0.9381 |
| Case 4 | 0.14696 | 0.20546 | 22.50 | 135.95 | 0.9364 |
Cooperative dispatch results.
Case 1: Heat and hydrogen balance constraints are not considered;
Case 2: Only heat balance constraints are considered;
Case 3: Only hydrogen balance constraints are considered;
Case 4: Both hydrogen and heat balance constraints are considered.
It can be seen from Table 2 that where the hydrogen and heat balance constraints are ignored (Case 1), the dispatchable power of the wind plant is equal to the predicted value, and the operation cost of conventional units is the lowest. However, 25.04 kg extra hydrogen is needed to support the 24-h operation of the hybrid system. As constraints are added into the model, the availability coefficient of the wind plant decreases to meet hydrogen and heat balance requirements, and the cost of conventional units increases. Compared with Case 1, the wind plant availability coefficient in Case 4 has dropped by 6.33%, and the cost of conventional units has increased by 3.4%. The results of conventional unit dispatch in the four scenarios are shown in Figure 4.
FIGURE 4
Dispatch results of conventional units.
According to the results in Case 2 and 3, we can figure out that the dispatch strategy that only satisfies the hydrogen balance or the heat balance requires the extra hydrogen or thermal supplement. It can be seen that when the system has sufficient day-ahead initial hydrogen storage or thermal storage, the corresponding dispatch strategy can also be adapted to increase the grid-connected power output of wind plants. When the initial hydrogen energy storage and the thermal storage are insufficient, the wind plant availability coefficient can only be reduced, and the dispatch strategy in Case 4 is adopted to maintain this efficient and economical operation of the HESS. Figure 5 shows the operation of the HESS in the four Cases.
FIGURE 5
Operation of hydrogen energy storage system.
It can be seen from Figure 5 that the PEMFC in Case 2 consumes extra hydrogen to provide an up-reserve for the hybrid system, improving the availability coefficient of wind power and meeting the heat balance constraints. In Case 3, the availability coefficient of the wind plant decreases. We can see that the probability of wind power curtailment increases while the probability of load curtailment decreases. The W-HHS calls the electrolyzer to increase hydrogen production, and the PEMFC is controlled to consume less hydrogen, thereby satisfying the hydrogen balance constraint. The SOC changes of hydrogen balance and heat balance are shown in Figure 6.
FIGURE 6
SOC change of hydrogen energy storage system.
Figure 6 shows that only the dispatch strategy adopted in Case 4 can simultaneously satisfy the heat and hydrogen balance constraints. In this case, the hydrogen tank will have the remaining hydrogen. The remaining hydrogen can be used in the subsequent optimal dispatch to increase the grid-connected output of the wind plant, reflecting the character of hydrogen energy storage’s cross-cycle energy storage. When the heat balance SOC shown in Figure 6 is less than 0, it indicates that the hydrogen energy storage system must rely on additional supplementary heat energy to maintain the regular operation. This working condition means that the hydrogen energy storage system’s thermal energy generated and stored is not enough to maintain the heat balance of intermittent work. In practice, the system’s hydrogen production and power generation efficiency will decline to generate sufficient thermal energy.
Compared with Case 3, the optimal strategy in Case 4 further reduces the wind plant’s available power and grid-connected power to satisfy the heat balance. It can be seen that the heat balance constraint reduces the profits of the hybrid system. Thus, the factors that affect the heat balance constraint should be further analyzed to reduce the impact of heat balance on the operation of HESS.
Analysis of Impact Factors on Heat Balance
From the results in Figures 5, 6, we can easily find out that when the W-HHS is in operation, the electrolyzer and PEMFC run intermittently. The electrolyzer cannot run at full power, and the PEMFC is always in reserve mode for a long time in an operation cycle. According to the heat balance model given in Eq. 5, 15, whether the heat produced can meet the heat lost by the system mainly depends on the heat dissipation coefficient of the electrolyzer and PEMFC. Therefore, the influence of the heat dissipation coefficient on the heat balance SOC constraint should be analyzed. Thus, we take the heat dissipation coefficient in a range of [1/900, 1/1,500] to study the performance of the waste heat utilization system. The results are shown in Table 3.
TABLE 3
| Heat dissipation coefficient (°C−1) | Wind plant availability coefficient | Residual hydrogen (kg) | Residual thermal energy (kW⋅h) |
|---|---|---|---|
| 1/900 | 0.9241 | 27.96 | 61.5 |
| 1/1,000 | 0.9346 | 22.5 | 135.95 |
| 1/1,100 | 0.9379 | 17.9 | 812.5 |
| 1/1,200 | 0.9431 | 13.98 | 996.9 |
| 1/1,300 | 0.9446 | 12.79 | 1,187.7 |
| 1/1,400 | 0.9446 | 12.79 | 1,365.8 |
Influence of heat dissipation coefficient.
Table 3 shows that as the heat dissipation coefficient decreases, the operation is less influenced by the heat balance constraints. The high availability coefficient of the wind plant results in a decreased operation probability of the electrolyzer and an increased working time of the PEMFC. Hydrogen storage drops accordingly. Also, due to the decrease of heat dissipation, the thermal energy stored in the thermal tank gradually increases.
When the heat dissipation coefficient is less than 1/1,200, Case 4 is equivalent to Case 3. The hybrid system is no longer affected by heat balance constraints. Thus, we can conclude that improving the heat insulation performance and reducing the heat dissipation power will help raise the system’s operating efficiency. Meanwhile, the available grid-connected power of the wind plant increased, and the profit of the W-HHS improved.
Discussion
According to the W-HHS energy storage system’s electricity, hydrogen, and heat characteristics, this work designs a HESS with a waste heat utilization system. Then, the hydrogen and heat dual SOC constraints are established to describe the operation of the electrolyzer and PEMFC in the HESS. On this basis, a cooperative dispatch method for the W-HHS is proposed, aiming to reduce the operation cost of conventional units in the grid and increase wind plants’ profit. The distributionally robust optimization approach models the uncertainties in the system. Then, the formulated model is transformed into a set of chance constraints and can be solved with an offline solver. Finally, we test our model with the IEEE 30-bus system, and actual wind power data from a Qinghai province wind plant is used.
The results show that without considering the hydrogen balance, the system has a significant wind power availability factor, which can improve the grid-connected power of the wind plant, but in practice, it needs to consume additional hydrogen to ensure fuel cell power generation. It can be seen that considering the hydrogen balance and heat balance will reduce the grid connection availability factor of the wind plant to generate enough hydrogen and thermal energy. In the operation of W-HHS, the available power of the wind plant is 6–8% lower than the ideal situation to meet the SOC constraints of hydrogen and heat balance. The heat dissipation coefficient of HESS has a significant influence on the SOC constraint of heat balance. When the heat dissipation coefficient is small enough, the SOC constraint of heat balance becomes invalid, and the SOC constraints of hydrogen balance only restrict the system operation strategy.
In addition, it should also be seen that after 24 h of continuous operation of the case4 in this paper, both hydrogen storage tank and heat storage tank have residual energy. This shows that in multi-day dispatching scene when there is hydrogen and heat energy in the hydrogen storage tank and heat storage tank, the system dispatching can also adopt an optimistic availability factor to use the remaining hydrogen and heat energy to improve the grid-connected power of wind plant and reduce the generation cost of conventional units. This problem needs further research in the future.
Statements
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
YS: conceptualization, writing-original draft preparation, software. LC: funding acquisition, validation. LM: funding acquisition, conceptualization, methodology, and editing. HM: validation. MG: project administration. SM: supervision. All authors have read and agreed to the published version of the article.
Funding
This research was funded in part by the Joint Fund Project of National Natural Science Foundation of China (U1766203), in part by Key R and D and Transformation Plan of Qinghai Province (2021-GX-109), and in part by the Basic Research Project of Qinghai Province (2021-ZJ-948).
Acknowledgments
The authors would like to thank Prof. Feng Liu and Wei Wei for very helpful discussions on the cooperative game, Dr. Xiaodai Xue and Xuelin Zhang for discussions on the HDR geothermal generator model, Dr. Zhongjie Guo and Jiayu Bai for suggestions on practical issues.
Conflict of interest
The 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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Nomenclature
- DRCC
Distributionally robust chance-constrained
- HESS
Hydrogen energy storage system
- PEMFC
Proton-exchange membrane fuel cell
- SOC
State of charge
- W-HHS
Wind-hydrogen hybrid system
Variables
electrolyzer heat dissipation area per unit capacity
PEMFC heat dissipation area per unit capacity
- ,,
operation cost coefficients of the conventional unit
imaginary parts of elements in network admittance matrix
on-grid power tariff
penalty coefficient
Wasserstein divergence
real parts of elements in network admittance matrix
higher heating value of hydrogen
number of data samples
hydrogen consumption rate of PEMFC
set formed by all possible distributions
rate of hydrogen production at time
probability distribution of the actual wind plant output
empirical distribution of the wind plant output
electricity consumed by the electrolyzer at time
electricity consumed by the fuel cell at time
power output of the units at time
active power injected into bus
active power flow on line at time
active power flow errors on line
active load into bus
installed power of the wind plant
dispatchable power of the wind plant at time
predicted output power of the wind plant at time
actual power output of the wind plant at time
probability distribution of the actual wind plant output
heat produced by the electrolyzer at time
heat produced by the fuel cell at time
reactive power injected into bus
reactive power flow on line at time
reactive power flow errors on line
reactive load into bus
heat provided for extra heat loads
heat consumed by the HESS
unit area thermal resistance of electrolyzer
unit area thermal resistance of PEMFC
maximum capacity of the hydrogen tank
capacity of line
maximum capacity of the thermal tank
SOC of hydrogen energy stored in the tank at time
SOC of thermal energy stored in the tank at time
ambient temperature
operating temperature of the electrolyzer
operating temperature of the PEMFC
temperature of the water entering the electrolyzer
voltage magnitude of bus .
confidence level
efficiency of the electrolyzer
efficiency of the heat exchanger
efficiency of PEMFC
Phase angle of bus
heat dissipation coefficient of unit capacity electrolyzer
heat dissipation coefficient of unit capacity PEMFC
power output impact factor of wind plant
predicted power output impact factor of wind plant
- ,
dual variables
availability factor of wind plant
radius of the ambiguous set
total power deviation in the whole period of dispatch
time interval
References
1
AlirahmiS. M.RazmiA. R.ArabkoohsarA. (2021). Comprehensive Assessment and Multi-Objective Optimization of a green Concept Based on a Combination of Hydrogen and Compressed Air Energy Storage (CAES) Systems. Renew. Sust. Energ. Rev.142, 110850. 10.1016/j.rser.2021.110850
2
EndoN.GoshomeK.TetsuhikoM.SegawaY.ShimodaE.NozuT. (2021). Thermal Management and Power Saving Operations for Improved Energy Efficiency within a Renewable Hydrogen Energy System Utilizing Metal Hydride Hydrogen Storage. Int. J. Hydrogen Energ.46 (1), 262–271. 10.1016/j.ijhydene.2020.10.002
3
EsfahaniP. M.KuhnD. (2018). Data-driven Distributionally Robust Optimization Using the Wasserstein Metric: Performance Guarantees and Tractable Reformulations. Math. Programming171 (1), 115–166. 10.1007/s10107-017-1172-1
4
García ClúaJ. G.MantzR. J.De BattistaH. (2018). Optimal Sizing of a Grid-Assisted Wind-Hydrogen System. Energ. Convers. Manag.166, 402–408. 10.1016/j.enconman.2018.04.047
5
GuoZ.WeiW.ChenL.XieR.MeiS. (2020). Sizing Energy Storage to Reduce Renewable Power Curtailment Considering Network Power Flows: a Distributionally Robust Optimisation Approach. IET Renew. Power Generation14 (16), 3273–3280. 10.1049/iet-rpg.2020.0354
6
HotaA. R.CherukuriA.LygerosJ. (2019). Data-driven Chance Constrained Optimization under Wasserstein Ambiguity Sets. Am. Control. Conf. (Acc). 10.23919/ACC.2019.8814677
7
HwangJ. J. (2006). Thermal-electrochemical Modeling of a Proton Exchange Membrane Fuel Cell. J. Electrochem. Soc.153 (2), A216. 10.1149/1.2137652
8
ImpramS.Varbak NeseS.OralB. (2020). Challenges of Renewable Energy Penetration on Power System Flexibility: A Survey. Energ. Strategy Rev.31, 100539. 10.1016/j.esr.2020.100539
9
JiayuB.XingangW.ChaoshanX.ZhiyongY.ShoutaoT.HeC. (2021). “Development Status and Measures to Promote the Development of Renewable Energy in China,” in 2021 3rd Asia Energy and Electrical Engineering Symposium (Chengdu, China: AEEES), 1102–1107. 10.1109/AEEES51875.2021.9403132
10
KovačA.MarciušD.ParanosM. (2020). Thermal Management of Hydrogen Refuelling Station Housing on an Annual Level. Int. J. Hydrogen Energ.46 (57), 29400–29410. 10.1016/j.ijhydene.2020.11.013
11
LiZ.ZhangW.ZhangR.SunH. (2020). Development of Renewable Energy Multi-Energy Complementary Hydrogen Energy System (A Case Study in China): A Review. Energy Exploration & Exploitation38 (6), 2099–2127. 10.1177/0144598720953512
12
Mansour-SaatlooA.Agabalaye-RahvarM.MirzaeiM. A.Mohammadi-IvatlooB.AbapourM.ZareK. (2020). Robust Scheduling of Hydrogen Based Smart Micro Energy Hub with Integrated Demand Response. J. Clean. Prod.267, 122041. 10.1016/j.jclepro.2020.122041
13
MirzaeiM. A.YazdankhahA. S.Mohammadi-IvatlooB. (2018).Integration of Demand Response and Hydrogen Storage System in Security Constrained Unit Commitment with High Penetration of Wind Energy. In Electrical Engineering (ICEE), Iranian Conference on. IEEE, 1203–1208. 10.1109/ICEE.2018.8472631
14
PanG.GuW.LuY.QiuH.LuS.YaoS. (2020). Optimal Planning for Electricity-Hydrogen Integrated Energy System Considering Power to Hydrogen and Heat and Seasonal Storage. IEEE Trans. Sustain. Energ.11 (4), 2662–2676. 10.1109/tste.2020.2970078
15
PereiraM. V.GranvilleS.FampaM. H. C.DixR.BarrosoL. A. (2005). Strategic Bidding under Uncertainty: a Binary Expansion Approach. IEEE Trans. Power Syst.20 (1), 180–188. 10.1109/tpwrs.2004.840397
16
PetkovI.GabrielliP. (2020). Power-to-hydrogen as Seasonal Energy Storage: an Uncertainty Analysis for Optimal Design of Low-Carbon Multi-Energy Systems. Appl. Energ.274, 115197. 10.1016/j.apenergy.2020.115197
17
QiuY.YangH.XuY.ChenW.WangP. (2020). Stability-constrained Two-Stage Robust Optimization for Integrated Hydrogen Hybrid Energy System. Csee Jpes7 (1), 162–171. 10.17775/CSEEJPES.2020.00810
18
ShuZ.JirutitijaroenP.Leite da SilvaA. M.SinghC. (2014). Accelerated State Evaluation and Latin Hypercube Sequential Sampling for Composite System Reliability Assessment. IEEE Trans. Power Syst.29 (4), 1692–1700. 10.1109/TPWRS.2013.2295113
19
SiY.ChenL.ZhangX.ChenX.ZhengT.MeiS. (2021). Game Approach to HDR-TS-PV Hybrid Power System Dispatching. Appl. Sci.11 (3), 914. 10.3390/app11030914
20
WeiF.SuiQ.LiX.LinX.LiZ. (2021). Optimal Dispatching of Power Grid Integrating Wind-Hydrogen Systems. Int. J. Electr. Power Energ. Syst.125, 106489. 10.1016/j.ijepes.2020.106489
21
WenT.ZhangZ.LinX.LiZ.ChenC.WangZ. (2020). Research on Modeling and the Operation Strategy of a Hydrogen-Battery Hybrid Energy Storage System for Flexible Wind Farm Grid-Connection. IEEE Access8, 79347–79356. 10.1109/ACCESS.2020.2990581
22
XiaoP.HuW.XuX.LiuW.HuangQ.ChenZ. (2020). Optimal Operation of a Wind-Electrolytic Hydrogen Storage System in the Electricity/hydrogen Markets. Int. J. Hydrogen Energ.45 (46), 24412–24423. 10.1016/j.ijhydene.2020.06.302
23
XiongY.SiY.ZhengT.ChenL.MeiS. (2021). Optimal Configuration of Hydrogen Storage in Industrial Park Integrated Energy System Based on Stackelberg Game. Trans. China Electrotechnical Soc.36 (3), 507–516. 10.19595/j.cnki.1000-6753.tces.200576
24
XuX.HuW.CaoD.HuangQ.LiuZ.LiuW.et al (2020). Scheduling of Wind-Battery Hybrid System in the Electricity Market Using Distributionally Robust Optimization. Renew. Energ.156, 47–56. 10.1016/j.renene.2020.04.057
25
YangL.XuY.GuW.SunH. (2020). Distributionally Robust Chance-Constrained Optimal Power-Gas Flow under Bidirectional Interactions Considering Uncertain Wind Power. IEEE Trans. Smart Grid12 (2), 1. 10.1109/TSG.2020.3029027
26
YuD.ZhangT.HeG.NojavanS.JermsittiparsertK.GhadimiN. (2020). Energy Management of Wind-PV-Storage-Grid Based Large Electricity Consumer Using Robust Optimization Technique. J. Energ. Storage27, 101054. 10.1016/j.est.2019.101054
27
ZhangG.WanX. (2014). A Wind-Hydrogen Energy Storage System Model for Massive Wind Energy Curtailment. Int. J. Hydrogen Energ.39 (3), 1243–1252. 10.1016/j.ijhydene.2013.11.003
28
ZhangQ.ChenW. (2020). Modeling China's Interprovincial Electricity Transmission under Low Carbon Transition. Appl. Energ.279, 115571. 10.1016/j.apenergy.2020.115571
29
ZhuF.ZhongP.-a.XuB.LiuW.WangW.SunY.et al (2020). Short-term Stochastic Optimization of a Hydro-Wind-Photovoltaic Hybrid System under Multiple Uncertainties. Energ. Convers. Manag.214, 112902. 10.1016/j.enconman.2020.112902
Summary
Keywords
hydrogen energy storage, hybrid system, cooperative dispatch, distributionally robust method, heat balance, hydrogen balance
Citation
Si Y, Chen L, Ma L, Gao M, Ma H and Mei S (2021) Modeling the Heat-Hydrogen Balance Characteristic of Hydrogen Energy Storage and Cooperative Dispatch of Wind-Hydrogen Hybrid System. Front. Energy Res. 9:791829. doi: 10.3389/fenrg.2021.791829
Received
09 October 2021
Accepted
03 November 2021
Published
22 November 2021
Volume
9 - 2021
Edited by
Zhenhao Tang, Northeast Electric Power University, China
Reviewed by
Xiao Wang, Wuhan University, China
Bingke Yan, State Grid Hubei Electric Power Company, China
Updates
Copyright
© 2021 Si, Chen, Ma, Gao, Ma and Mei.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yang Si, siyang@qhu.edu.cn; Laijun Chen, chenlaijun@qhu.edu.cn
This article was submitted to Wind Energy, a section of the journal Frontiers in Energy Research
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.