Underground Hydro-Pumped Energy Storage Using Coal Mine Goafs: System Performance Analysis and a Case Study for China

In response to the Paris climate agreement, the Chinese government has taken actions to improve the energy structure by reducing the share of coal-fired thermal power and increasing the use of clean energy. However, due to the extreme shortage of large-scale energy storage facilities, the utilization efficiency of wind and solar power remains low. This paper proposes to use abandoned coal mine goafs serving as large-scale pumped hydro storage (PHS) reservoir. In this paper, suitability of coal mine goafs as PHS underground reservoirs was analyzed with respects to the storage capacity, usable capacity, and ventilation between goaf and outside. The storage capacity is 1.97 × 106 m3 for a typical mining area with an extent of 3 × 5 km2 and a coal seam thickness of 6 m. A typical goaf-PHS system with the energy type α w = 0.74 has a performance of 82.8% in the case of annual operation, able to regulate solar-wind energy with an average value of 275 kW. The performance of the proposed goaf-PHS system was analyzed based on the reservoir estimation and meteorological information from a typical region in China. It has been found that using abandoned coal mine goafs to develop PHS plants is technically feasible in wind and solar-rich northwestern and southwestern China.


INTRODUCTION
China's coal-heavy primary energy structure causes environmental pollution and massive carbon dioxide emission . In 2016, China produced 10.5 billion tons of carbon dioxide emission by fossil fuel burning, severely affecting the global climate change. In response to the Paris climate agreement, the Chinese government is taking actions to improve the energy structure, that is, to increase the share of clean energy while reducing the use of fossil fuel (Pan et al., 2018;Li and Chen, 2019). However, due to the extreme shortage of large-scale energy storage (LSES) facilities, the utilization efficiency of wind and solar power remains low (Bai et al., 2015). By the end of 2018, 22 operational PHS plants were distributed in the coastal region or in eastern China (Figure 1), with an installed capacity of 19.23 GW, entirely insufficient for the largest wind-solar power market in China (Davidson et al., 2016).
Among the existing energy storage technologies, only compressed air energy storage (CAES) and pumped hydroelectric storage (PHS) are cost-effective at large temporal scales, from several hours to many days (Zhao et al., 2015;Rogeau et al., 2017;Collado et al., 2018;Pali and Vadhera, 2018;Tung et al., 2018). Large-scale CAES is known to require specific geological conditions to store high pressure air/gas (Parkes et al., 2018), while the PHS appears to be relatively more mature and reliable (Hunt et al., 2014;Kapila et al., 2017;Ruppert et al., 2017;Hunt et al., 2018). In the wind and solar energy rich northwestern and southwestern China, scarce PHS plants are available for LSES for the abundant energy ( Figure 1), because of the difficulty in site selection and water shortage.
As of now, the number of abandoned mines, including coal mines, metal mines, and chemical mines (non-metallic mine), is estimated to exceed one million, worldwide. In China, the number of documented closed coal mines reaches 3,868 with more than a productive capacity of 350 million (Xie and Liu, 2018) and the other abandoned coal mines are roughly estimated at 12,000 (Zhongguang et al., 2019). During extraction of mineral resources from underground deposits, shafts and extensive galleries are excavated. These underground voids are usually left to be flooded, but often perpetual costs related to pumping to keep a safe water level or water treatment have to be maintained, becoming long-term liabilities. Some renewable/ sustainable post-mining solutions including underground reservoirs (Andrés et al., 2017), hydro-pumped energy storage (PHES) (Menéndez et al., 2017;Pujades et al., 2017), compressed air energy storage (CAES) (Kim et al., 2012;Fan et al., 2018a), thermal energy storage (Al-Habaibeh et al., 2018), and underground hydrogen storage (Simon et al., 2015) were proposed and proved viable by some researchers. The first documented CAES test using an abandoned mine was conducted in Japan. In the United States, an abandoned limestone mine was planned to be converted into a CAES plant with a capacity of 2700 MW. In South Africa, the viability of deep level gold mines transferred into underground pumped hydroelectric energy storage in the Far West Rand gold field were by Winde (Winde et al., 2017). It is reported that the Prosper-Haniel coal mine, Germany, planned to be converted into a 200 MW PHS plant Fan et al., 2017). Menéndez et al. (2018) proposed to develop underground PHS plants by using underground tunnel/roadway networks of more than 30 closed coal mines in the Asturian Central Coal Basin of Spain. However, the roadway space in a mine is limited and insufficient for LSES.
Considering the goafs' considerable space capable for LSES and simultaneously reduce water evaporation, a concept of PHS system using underground coal mine goafs (goaf-PHS) was put forward (Fan et al., 2020). In arid regions of northwest China, waterpreserved mining is the only option for coal mining. Wang et al. investigated the basic conditions and mechanisms of goaf water storage, providing an effective approach for large-scale goaf water reservoirs in the arid regions . Gu et al. proposed  Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 760464 a technical approach of storage and utilization of mine water through underground reservoirs in coal mine. Currently, 32 coal mine underground reservoirs have been established in the Shendong Mining Area and can effectively alleviate the scarcity of local production and domestic water (Gu, 2015). Aforementioned researches demonstrate the enormous potential capacity of underground coal mine goafs served as water reservoirs for PHS plants. Nevertheless, the details on the goaf reservoirs, such as the storage capacity, the useful volume, environmental issues, and the performance of the goaf-PHS system, are still unclear. This paper proposes a hybrid PHS system using underground coal mine goafs, as shown in Figure 2. The study performed a suitability analysis on goafs serving as PHS with respect to the three issues of most concern: storable volume, usable volume, and fluid exchange. Based on the estimation results of the goaf reservoir and meteorological data (solar radiation and wind speed) in northwestern China, the performance of the goaf-PHS system is evaluated. At last, the feasibility of goaf-PHS plants is discussed on the national strategic level.

SYSTEM DESCRIPTION AND ESTIMATION
The upper reservoir, which provides water storage capacity at a high level, uses underground goafs or a surface reservoir. The lower reservoir is located at the underground coal mine goafs with a greater depth to ensure a suitable water head ( Figure 2). The other components of the hybrid PHS are similar with a conventional PHS plant, including wind turbines, photovoltaic generators, reversible hydraulic pumps/turbines, and penstocks. Wind turbines and photovoltaic generators are the energy receiving-conversing end of the system converting clean power into electricity. Reversible hydraulic pumps/turbines are the energy conversing-transferring devices. When electricity generation is greater than demand, the surplus electricity will drive pumps to elevate the water from the lower to the upper reservoir, storing the potential energy. When the electricity generated by wind and solar power cannot satisfy the demand, the stored water will be released to produce electricity and fill the demand gap. The two reservoirs are the energy storage facility, reserving the surplus energy in the form of water potential energy. The penstocks, surge tank, and a variety of connecting tunnels/ rooms provide ancillary services for the system's smooth operation. The system efficiency is mainly decided by the performance of receiving-conversing-transferring devices.

Wind Turbines and Photovoltaic Generator
The energy receiving-conversing devices are mainly composed of wind turbines, photovoltaic generators, and transformer station. Wind turbines and photovoltaic generators convert the wind and solar power into electricity, respectively. Their power P w and P s can be expressed by (Boumaaraf et al., 2018;Pali and Vadhera, 2018;Yesilbudak, 2018) where c t is the efficiency of wind turbines capturing the wind. ρ a represents the air density. R is the length of the blades on the wind turbines and u a is the wind velocity. α, η so , and η L are solar fraction, efficiency of photovoltaic generator, and heat loss rate. J s is the solar radiation intensity in W/m 2 . η L is the heat loss rate of photovoltaic generator. Before using the electricity with hydraulic pumps, electricity transmission is necessary to obtain an appropriate voltage, which may cause some electricity losses (Han et al., 2017). Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 760464 P tm η tm (P w + P s ), where η tm is the comprehensive efficiency of the electricity transmission. P tm and η tm are the output and overall conversion efficiency of the electricity transmission station.

Hydro-Pump/Turbine and Penstock
The energy transferring devices include hydro-pump/turbine and penstock. The input power of a pump P hp, in is related to the excess energy transmitted via the transmission station. The water flow pumped from lower reservoir to upper reservoir can be calculated by (Novara and McNabola, 2018) where ρ w and g are water density and gravitational acceleration, respectively. H av is the lifting height, i.e., the elevation difference between the two reservoirs. Q hp and Q ht are the flow rates. η hp is the efficiency of hydro-pumps. The hydraulic turbine converts the potential energy of water in the upper reservoir into electricity, in which output power P ht, out can be evaluated by P hp, out η ht ρ w gH av Q ht . (4) The efficiency of a hydro-pump/turbine normally varies from 0.86 to 0.95, with water flow, head, etc. In this paper, for simplicity, it is considered constant and equal to 0.91.
Transmission over long distances results in loss of hydraulic energy, due to the friction between the fluid and pipe wall. The pressure loss P of water flow in the penstocks can be evaluated by (Fan et al., 2018b;Fan et al., 2019a) where λ is the drag efficiency and can be calculated using absolute friction. u w represents the water velocity inside the pipes. L and d are the length and diameter of the penstock. Ru is the absolute roughness of penstock interior walls.

Goaf Reservoir
As coal is mined, the remained space in goafs was filled with rubbles from proof and floor surrounding rocks ( Figure 3). Above the rubbles, the rock formations generate a large number of fractures and cracks due to the large displacement . In the further rock, small displacement just induces deformation without fractures forming. The three overlaying rock zones of different fractures above the mined coal seam are known as "vertical three zones," namely, the caving zone, the fissure zone, and the displacement zone Shu et al., 2019). After coal mining, massive fissures and pores, capable of storing water, form in the goaf. The caving zone and fissure zone contain massive pores and fissures, capable of storing water. The rock in the displacement zone only displays some microfractures, which has no storage capability, but has a considerable permeability due to the formation of microfractures (Liu et al., 2020). The range of the caving zone and the fissure zone would vary with lithology of the overlying rock and thickness of the mined coal seam and can be estimated with numerical simulations (Dong et al., 2016), physical detection (Deng et al., 2018;Ren and Wang, 2020), or empirical models (Peng, 1984). To facilitate calculation, an empirical model was used. If the overlaying rocks are hard rocks, the height of the caving zone (h 1 ) and the height of the fissure zone (h 2 ) can be evaluated by Eqs. 6a,b (Peng, 1984), where M is the thickness of the mined coal seam. The storage space (V sg ) in the caving zone and the fissure zone can be estimated by FIGURE 3 | Diagram for "three zones" of overlaying rock strata.
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 760464 where n m , L g , and L t are the number, the length, and the width of the goaf. F is a usability coefficient for the reservoir, while f represents the goaf storage coefficient, determined by the capacity of pores and fissures.
Due to the stress release in the process of mining, the rock volume would expand with the formation of new fractures and pores. After mining, the rubbles and fracture rocks will be recompacted under the gravity of overlying strata. But there is no way to completely restore the volume. The expansion coefficient K of rock mass is defined as the ratio of the expanded rock volume to the original volume, to characterize the volumetric expansion behavior of the surrounding rock after fracturing. Coefficient K and the zone range vary with the pressure, lithology, and distance to the mined coal. The storage coefficient therefore can be calculated using K, The expansion coefficient K of the rock mass at different positions can be evaluated by (Meng et al., 2016a) Mining activities also influence its underlying rocks, which will move toward the remained space after the coal is excavated, leading to new fractures. Compared with caving zone and fissure zone, the new space created within underlying rock is relatively small, thus not included in the storage space. But its permeability changes a lot, enhancing the water seepage. The influence depth h 3 in underlying rocks is estimated by (Fan et al., 2019b) h 3 x a cos φ 2 cos π where C and φ are the cohesion and friction angle of the underlying rock. φ m is the friction angle of coal. ζ, c, and h 0 are respectively the coefficient of stress concentration, unit weight of overlying rock, and depth of the mined coal seam. Water and air seepage in the goafs through rocks are assumed to be slow and follow the Darcy's law. The water flow velocity v is determined by v k μ wa dp dr , where r is the distance, k is the permeability, and p is the water pressure. Using the finite element method, the seepage velocity of water through the upper reservoir boundary can be determined.

GOAF RESERVOIR FEASIBILITY ANALYSIS
After coal excavation, overburden rocks collapse and fracture, making various sizes of fractured rock blocks to pile up within the goafs. The goaf reservoir can be considered as a type of porous medium. The storage capacity depends not only on the goaf volume, but especially on the interspace between blocks. During charging and discharging, water and air would exchange, and the permeability within the goaf determines whether the water or air can flow in and out smoothly or not. Using the expansion coefficient K, the ratio of rock expanded volume after collapsing to its original volume, usable storage capacity, and ventilation of a goaf were evaluated.

Storable Volume
Expansion coefficient K characterizes the volumetric expansion of the rock mass overlying the coal seam after mining. It varies with the pressure, lithology, and mining activities. Based on the K behavior obtained from the literature (Meng et al., 2016a), Figure 4 shows the goaf volume of a single goaf 200 m wide and 3,500 m long vs water level. As the coal seam thickness increases, the total volume V1 rises almost linearly. Considering the resistance of waterproof walls surrounding the water reservoirs, an admissible capacity V2 is calculated according to the suggested maximum water head (17.8 m) (Fan et al., 2020). Storage ratio R1, defined as the ratio of total capacity to the volume of excavated coal, shows a deceleration downtrend. Based on the production history of the shutdown coal mines issued by the State Administration of Coal Mine Safety and China Coal Industry Association, the storage space of the abandoned coal mine goafs could be calculated. It reaches a considerable value of 4.70 × 108 m 3 . Figure 1 shows its distribution in the Mainland China. Taking a typical abandoned (shut down) coal mine, Yima Qianqiu coal mine in Henan province, with a dimension 3 × 5 km 2 and a coal thickness of 6 m as an example, the goaf storage capacity with a water level of 17.8 m is calculated at 1.97 × 106 m 3 .

Usable Volume
To efficiently pump/inject water during charging/discharging, water-collecting wells are made at the center of the goaf reservoir bottom. Water/air flow velocity within the reservoir is determined by pressure gradient and permeability. To ensure sufficient outputs of turbine and pumps, the maximum water flow should be guaranteed. The water saturation lines within the goafs were marked when the water flow in penstocks was constant at 6.25 m 3 /s. The maximum water level in the WCWs was set at 17.8 m in similar consideration of the maximum water head restriction [48,49]. The minimum water level is 0. Figure 5A, B show the highest saturation line during water injection and the lowest saturation line during water releasing. As the permeability decreases, the saturation line tends to decline during injection and rise during releasing, implying a less water volume injected Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 760464 during filling and released during draining the reservoir, and hence a smaller usable capacity of the goaf reservoirs. Subtracting the corresponding saturation lines in Figures 5A, B, the usable coefficient relying on the permeability would be obtained for the goaf reservoirs, which are shown in Figure 6. It can be calculated that the permeability should be above 10 −7 m 2 to have a considerable usable capacity ( ≥ 80% of storage capacity).
Literatures (Konicek et al., 2013;Meng et al., 2016b) reported that the goaf permeability within the overburden rocks above 0-17.8 m ranges from 10 −10 to 10 −6 m 2 which suggests that goafs are highly likely to be able to serve as PHS reservoirs and have a considerable usable capacity.

Fluid Exchange
Ventilation shafts are excavated to connect the atmosphere with goafs for air smooth exchange during water pumping and injecting ( Figure 7). The location of the shaft is selected at the center of the goaf and its bottom is 2 m above the maximum water level in the reservoir. The shaft dimension closely effects pressure loss of passing air. The length is determined by the nature of goaf (depth of coal seam). Figure 8 shows the pressure loss (the pressure difference between Point P0 and P1) decreases dramatically by five orders of magnitude, as the diameter of the ventilation shaft increases from 0.5 to 5 m with the concrete shaft wall lining material. When the shaft diameter is 0.8 m, pressure loss decreases to a negligible value (∼1 KPa). As the permeability increases from 10 −10 m 2 to 10 −7 m 2 , the pressure loss (between P1 and P2) of air passing through the goaf diminishes, as shown in Figure 8. It is suggested that goafs should have a permeability larger than 3 × 10 −9 m 2 to bring about an acceptable pressure loss during the air passing. This, however, is hardly feasible according to the estimation in Ref. (Alehossein and Poulsen, 2010;Poulsen et al., 2018) that the permeability within the goaf above 19.8 m varies between 10 −10 and 10 −13 m 2 . Dendritic horizontal ventilation tubes (Figure 2) are, therefore, advised for the safety and smoothness of water-air exchange. Figure 1 shows the distribution of PHS plants in Mainland China. The PHS existing and planned plants both are located mainly in eastern China with a purpose of adjusting the power supply from western China. However, in northwestern and southwestern China (marked by ovals), with abundant solar and wind resource and storable goaf space, there are no PHS plants in most provinces, where are massive goafs of abandoned coal which were shut down in recent years. Goaf-PHS plants could be constructed in those regions to receive and enhance the use of solar and wind power, which could help China to reduce carbon emissions, thus benefitting the global climate. The following part will take a typical region in that area for the case study to evaluate the performance of the goaf-PHS system.

Performance Indicator
Except for system efficiency η sys , the following parameters are defined as in Eqs 12, the power type, α w representing the share of imported wind energy in the total imported energy, and the regulated-energy per volume (REPV) representing the regulation ability of the goaf-PHS system.
FIGURE 4 | Goaf storage capacity with varying thickness of coal seam. V1 and V2 represent the total volume of goaf reservoir and admissible capacity when the water level is 17.8 m, respectively; R1 is the storage ratio, the ratio of total capacity to the volume of excavated coal.
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 760464 where P w and P s are the outputs of wind power generator and solar collector. T is a calculation period. V sg, pk and V sg,tr represent the peak and the trough value of water volume in the reservoir.

Parameters
Based on the above analysis, the pressure loss during water-air exchange in the goaf is negligible and the storage capacity of one goaf reservoir is temporarily set at 1.97 × 106 m 3 and the usable coefficient is 0.8. The altitude difference between the upper and lower reservoir is 100 m. The initial water level inside the upper and lower reservoirs is zero and 17.8 m, respectively. The load profile, wind velocity, and solar radiation intensity for 1 year or 1 day in Inner Mongolia, a typical richest area in wind and solar energy resources in north China, were used in the following performance evaluations. Detailed data can be seen in Refs. (Xiaolin et al., 2009;Ruichun and Bin, 2014;Wang et al., 2015).
PHS plants usually have several time scales of operation modes, among which the yearly case and daily case are the most representative. The performance of the goaf-PHS system in both yearly and daily operation cases is shown in Figure 9. The model of the goaf-PHS system is detailed. When the water level in the collection well drops to 0 m; to efficiently pump/inject water during charging/discharging, water-collecting wells are made at the center of the goaf reservoir bottom. Water/air flow velocity within the reservoir is determined by pressure gradient and permeability. To ensure sufficient outputs of turbine and pumps, the maximum water flow should be guaranteed. The water saturation lines within the goafs were marked when the water flow in penstocks was constant at 6.25 m 3 /s. The maximum water level in the WCWs was set at 17.8 m in similar consideration of the maximum water head restriction (Bai et al., 2017;Pang et al., 2018). The minimum water level is 0. Figures 5A,B show the highest saturation line during water injection and the lowest saturation line during water releasing. As the permeability decreases, the saturation line tends to decline during injection and rise during releasing, implying a less water volume injected during filling and released during draining the reservoir and hence a smaller usable capacity of the goaf reservoirs.

Yearly Operation Case
Based on the wind and solar electric capacity generated in 2017 (wind, 2.70 × 10 10 kW·h, and solar, 9.67 × 10 9 kW·h, respectively), we could see that the type of energy emplaced in the goaf-PHS consists of 74% of wind generated electricity and 26% of solar electricity. Using the parameters of a typical abandoned coal mine as the standard, the evaluation results for a yearly operations case are shown in Figure 9. Regarding the share of component supply and consumption of the goaf-PHS system over 8,760 h in 1 year, the goaf-PHS is charged every day throughout the 1 year ( Figure 9B), while the discharge stage is mainly from the middle of May to October ( Figure 9A), indicating that the system stores the surplus energy every day, while releasing that from May to October, during the period of energy deficit. The average regulated load is 275 kW as the altitude difference between the two reservoirs is 100 m. The maximum output of the PHS system is delivered in August, when energy supplied by wind is the weakest ( Figure 9A). The maximum input of the PHS system appears in January ( Figure 9B). The details of several selected days (January 1 st -4 th for winter, April 1 st -4 th for spring, August 1 st -4 th for summer, and October 1 st -4 th for autumn) are shown in Figure 10. It can be seen that the pump units work every day in the hours around noon to stores surplus electricity, even in August ( Figure 10B), when the wind energy is seriously inadequate. In January, the pump units work full time, while the turbine units stay inoperative. The turbine units provide the maximum work output 19 h in 1 day, except for 11:00 to 15:00, in August ( Figure 10A). In April and October, the turbines and pumps FIGURE 6 | Effective capacity coefficient of goaf reservoirs with varying goaf permeability.
FIGURE 7 | Diagram for water-air exchange of the goaf reservoir during water pumping/injecting. P0, P1, and P2 represent air pressure in the atmosphere, at the bottom of ventilation shaft, and at the far endpoint within the goaf, respectively.
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 760464 8 work with roughly the same period. It is therefore the water variation in the reservoirs that changes slowly in April and October, rapidly in January, and especially in August ( Figure 11). The water volume in the upper reservoir increases from October to May, reaching a peak of 1.22 × 106 m 3 and then decreases, down to a trough value of −3.6 × 105 m 3 (Figure 12). The average system efficiency in 1 year is 82.8%. The calculated REPV is 2.82 kW·h/m 3 . Different power types would bring about a different system performance. As seen from Figure 13, with the wind contribution   increasing (α w rising), the REPV for the yearly regulation case increases and then declines, reaching a peak around α w 0.25.

Daily Operation Case
For a typical abandoned coal mine, the system performance was evaluated on two selected days (one in winter, 15th January, shown in Figure 14A; the other one in summer, 15th July, shown in Figure 14B). The chosen energy type is α w 0.8 for 15th January and α w 0.5 for 15th July, according to the daily generated electricity capacity. The comprehensive system efficiency of the goaf-PHS is 73.6 and 77.5%, respectively, which are notably lower than that in yearly operation case, since the water flow is faster in daily regulation mode and more friction produces more energy loss. The average regulated load in 15th January is 160.4 MW, while it is 70.5 on the 15th July. The REPV representing the regulatable energy of a unit volume is 2.50 kW·h/m 3 and 1.07 kW·h/m 3 , respectively, which result from different energy types.
Using the load profile and wind and solar parameters in summer, 15th July, the REPV with various energy types calculated is shown in Figure 13. When the wind contribution increases (α w changing from 0 to 1), the REPV for the daily regulation case increases monotonically. It can be observed that in daily time scale, the pure wind electricity supply is the most

Environmental Issue
Goaf water quality, especially the pH value, is another major substantial concern on the construction of the goaf-PHS system. With a low pH value, acid goaf water can corrode equipment, release metal ions, even some heavy metals, damage the underground structures (like waterproof wall), and pollute the surrounding water bodies (González et al., 2018). As already known, most coal contains various amounts of (0.5-3%) sulfide that exists mainly (60-70%) as pyrite minerals. The elevated pH level can cause the precipitation or co-precipitation of metal ions such as aluminum (Al), cadmium (Cd), cobalt (Co), chromium (Cr), copper (Cu), iron (Fe), manganese (Mn), nickel (Ni), lead (Pb), and zinc (Zn) contained in the mine water and result in sediments in the seepage channels and reduce the permeability of rock mass, thus blocking the fluid exchange within the goaf (Pujades et al., 2019).
As already known, most coal contains various amounts of (0.5-3%) sulfide that exist mainly (60-70%) as pyrite minerals. After repeated pumping and injection, mine water is exposed to oxidation conditions and acidizes following the chemical reaction R1-R3: With a hypothesis of geological medium (residual coal and gangue) containing 1% pyrite, Pujades et al. conducted numerical simulations and show that the pH value would decrease continuously to 3.1-3.3 in both surface reservoir and underground goaf reservoir (which is considered as a porous medium in the study) during 30 days of repeated pumpings and injections (Pujades et al., 2018). Under the condition of 10% FIGURE 14 | Output of goaf-PHS system components and water volume variation for the daily regulation case. (A) Average component output and water volume variation of the system using the data of the 15th of January (α w 0.8); (B) average component output and volume water variation of the system using the data of 15th July (α w 0.5).
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 760464 carbonate such as calcite and 1% pyrite in the initial goaf environment, the research shows calcite would mitigate and precipitate around the surface reservoir (Pujades et al., 2016). Many underground mine reservoirs using coal goafs from Shennan coal mine district, northwest China, have been constructed. The quality analysis shows that the original well water is weakly alkaline (pH 7.1-7.8) and there is any substantial difference from natural water, except some obvious increases in permanent hardness and sulfate radical . A long-term water monitoring shows that with the circulation of water resources, several harmful elements have shown a significant remobilization, in spite of ultra-low sulfur deposited in Jurassic coal of the Shennan coal mine district. Taking purification measures, the value of the mine water could be controlled above the utilization standard (pH 6.9-7.1).
In consequence, hydrogeology and hydrochemistry investigations are also indispensable in the feasibility demonstration of goaf-PHS system. Necessary purification treatments are very important to the water safety. The coal gangue packed in goaf reservoirs, which contain clay mineral contents such as illite and kaolinite, could act as a useful adsorbent to reduce the contents of organic compounds and nitrogen in mine water . In one word, the environmental influence of goaf reservoirs should be a point of concern but also a solvable problem when developing the hybrid-PHS plants using abandoned coal mine goafs.

CONCLUSION
In this study, we proposed to use abandoned coal mine goafs serving as large-scale pumped hydro storage reservoirs. Firstly, the suitability of goafs as PHS underground reservoirs was analyzed with respect to the storage capacity, usable capacity, and ventilation between goafs and the outside. Then, the performance of the proposed goaf-PHS system was analyzed based on the meteorological information in a typical Chinese region.
1) For a typical mining area with an extent of 3 × 5 km 2 and a coal seam thickness of 6 m, the storage capacity is 1.97 × 106 m 3 . To ensure a significant usability coefficient ( ≥ 80%), the permeability within goafs should be above 10 −7 m 2 . Pressure loss of air passing through the vertical ventilation shaft with a diameter of 1 m is negligible. Dendritic tubes are suggested for horizontal ventilation. 2) With the energy type α w 0.74, a typical goaf-PHS system has an efficiency of 82.8% in yearly operation case, able to regulate a solar-wind power with an average value of 275 kW. Under the same conditions, the goaf-PHS system can regulate 239 MW of solar-wind power in daily operation mode. 3) Goaf water quality, especially the pH value, is another major substantial concern on the construction of the goaf-PHS system. Hydrogeology and hydrochemistry investigations are indispensable in the feasibility demonstration of goaf-PHS system. Necessary purification treatments are very important to the water safety.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
DJ: methodology, investigation, data curation, writing-original draft, and software resources, supervision, and funding acquisition. SC: investigation, data curation, and software. WL: conceptualization, methodology, investigation, writing-review, and editing. YR: investigation, writing-review, and editing. PG: resources, validation, software, and data curation. ZL: resources, validation, software, and data curation.