Modeling the Heat-Hydrogen Balance Characteristic of Hydrogen Energy Storage and Cooperative Dispatch of Wind-Hydrogen Hybrid System

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 . 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 , 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 , compressed air energy storage (Alirahmi et al., 2021), hot dry rock geothermal energy , 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 .
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.  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.  further proposed a distributionally robust chanceconstrained (DRCC) dispatch method for W-HHS in the day-ahead power market and transformed it into a standard secondorder 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 heathydrogen 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.

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. 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 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 P t ele and Q t ele , respectively, represent the electricity consumed and the heat produced by the electrolyzer at time t. _ n t H2 represents the rate of hydrogen production, and H HV represents the higher heating value of hydrogen. η ele 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 P t fuel and Q t fuel represent the electricity and heat produced by the fuel cell at time t, respectively. _ m t H2 represents the hydrogen consumption rate of PEMFC. η fuel 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 SOC t H2 represents the SOC of hydrogen energy stored in the tank at time t. S max H2 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 SOC t th represents the SOC of thermal energy stored in the thermal tank at time t. S max th represents the maximum thermal energy storage capacity of the thermal tank. Q t ele and Q t fuel represent the heat generated by the electrolyzer and the PEMFC, respectively. Q t sys and Q t load represent the heat consumed by the HESS and the heat provided for the extra heat loads, respectively. η ex represents the efficiency of the heat exchanger.
In Eq. 4, Q t sys 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, Q t sys can be expressed as: where T t ele and T t fuel represent the operating temperature of the electrolyzer and PEMFC. T t atm represent the ambient temperature. P ele and P fuel are the installed capacity of electrolyzer and PEMFC. _ n t w , c p , and T t w represent the molar mass flow rate, specific heat capacity and temperature of the water entering the electrolyzer. λ ele and λ fuel 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 A ele /A fuel to unit area thermal resistance R ele /R fuel . λ ele and λ fuel can be expressed as:

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 λ t W , a simplified model of the wind plant output can be expressed asP where P W represents the installed power of the wind plant.P t W 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λ t W represents the predicted power output impact factor. P t W and P t W represent the predicted output and dispatchable power of the wind plant.λ t W reflects the availability of power grid dispatch with wind power forecasting. P t W 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 Δ W 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 P represents the probability distribution of the actual power output of the wind plant.P represents the empirical distribution of the power output. M(Ξ) represents the set formed by all possible distributions satisfying Wasserstein divergence d w , 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 Frontiers in Energy Research | www.frontiersin.org November 2021 | Volume 9 | Article 791829 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 P(P t W ) represents the probability distribution of the actual wind plant output. c e represents the on-grid power tariff. c p represents the penalty coefficient. P t,i gen represents the power output of the units. a i , b i , and c i represent the operation cost coefficients of the conventional unit. i represents the bus index. N represents the set of all buses.

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 T. Thus, we have the hydrogen balance SOC constraint is: where SOC ini H2 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 SOC ini th is the initial value of the thermal tank. Besides of heat SOC constraint, heat balance also includes heat power balance constraint in operation.
where Q t storage 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 P t l and Q t l represent the active and reactive power flow on line l at time t. V t i and θ t i represent the voltage magnitude and phase angle of bus i. P t Ll and Q t Ll represent the active and reactive power flow errors on line l (Guo et al., 2020).
The constraint of line capacity is: where S t l represents the capacity of line l. Eq. 17 can be linearized by the outer approximation approach as The power balance constraint is where P t i and Q t i represent the active and reactive power injected into bus i, that is, the power generated by the generator connected to the bus minus the active/reactive load P t,i ld /Q t,i ld on the same bus. G ij and B ij 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 K is the number of data samples. α is the confidence level. 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.

System Parameters
In this section, we test our models and method on the IEEE 30bus 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.
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.

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.
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.
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.
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 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

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 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.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.