The Concept, Technical System and Heat Transfer Analysis on Phase-Change Heat Storage Backﬁll for Exploitation of Geothermal Energy

: In view of high ground stress, high geothermal temperature, and thermal hazard during deep mineral resource exploitation, the concept of phase-change heat storage backﬁll was put forward in this study. Further, the corresponding technical system was constructed and the main content involved in technical system, which is the optimized proportion of the backﬁll slurry added with phase-change materials (PCMs), was examined. Moreover, we elaborated upon the collaborative optimization of a backﬁll body’s mechanical and thermal properties and the mutual cooperation on backﬁll mining, geothermal energy exploitation, and simultaneous stope cooling. The heat transfer behavior of a backﬁll body plays a key role in technology system. We numerically simulated the heat transfer among a backﬁll body, surrounding rock, and airﬂow in the heat storage process, as well as the heat transfer between backﬁll body and cold ﬂuid during the heat release process. The temperature distribution of a backﬁll body at di ﬀ erent heat storage / heat release times—i.e., the temperature distribution and its evolution—with heat transfer were revealed and analyzed. This study can provide theoretical guidance for a phase-change heat storage backﬁll, as it has an important signiﬁcance for the collaborative exploitation of mineral resources and geothermal energy. the heat storage/heat release process of the backfill body), the heat transfer model of the phase-change backfill body was established and the temperature distribution and its evolution in the backfill body was numerically simulated and analyzed during the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy. The backﬁll body is a homogeneous and isotropic solid, and the thermophysical parameters remain constant. the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy. The backﬁll body is a porous medium with constant porosity, and all phases are in local thermodynamic equilibrium.


Introduction
The exploitation and utilization of deep mineral resources has become an inevitable trend with the gradual reduction and exhaustion of shallow mineral resources. The exploitation and utilization of deep mineral resources is also in line with the strategic requirements of "deep space, deep sea, deep blue, and deep land" put forward in "The National Medium-and Long-Term Program for Science and Technology Development (2006-2020)". The mining technology of deep mineral resources is an important field of "deep land" exploration and it is the high-end development of modern mining technology. Deep mining is in the special environment of "three high", which is high ground stress, high geothermal temperature, and high mine depth. Thus, it is bound to bring many scientific and technological problems to the mining of deep mineral resources and also seriously affect the safety, efficiency, and cost of mine production, even related to the feasibility of resource exploitation [1,2].
Backfill mining technology is an effective way to solve the problem of high ground stress, control the ground subsidence, and improve the mining rate of resources when mining deep mineral It is unnecessary for drilling engineering to occur on the earth's surface, as it can save the cost of geothermal energy exploitation.
Energies 2020, 13, x FOR PEER REVIEW 2 of 22 Backfill mining technology is an effective way to solve the problem of high ground stress, control the ground subsidence, and improve the mining rate of resources when mining deep mineral resources. Moreover, it does not cause serious damage to the ecological environment. With the rapid development of the mining industry, the situation of mineral resources is becoming more and more complex, and the geographical location and environmental conditions are also different for different mines. Therefore, it is necessary not only to select the corresponding backfill mining technology according to the specific mining conditions, but also to reasonably select the backfill materials in combination with the actual situation, so as to ensure its quality and performance can meet the mining operation need. The innovation and development of backfill materials explores a new way for resource reuse of solid wastes in mines, such as tailings, fly ash, and waste rocks. Under the premise of controlling ground subsidence, the solid waste of mines should be selected as backfill materials so backfill mining is also an effective means to make solid wastes harmless and reusable, and realize "green mining" in mines [3,4].
For deep mining, the ground temperature rises and heat releases to the stope from surrounding rocks; this is the main source of thermal hazard in the mine. In addition, other heat sources, such as air compression and mechanical equipment, lead to an increase of stope temperature. The deeper the mining depth, the higher the temperature of surrounding rock, and the more serious the thermal hazard. In many countries, thermal hazard in mines appears one after another [5][6][7][8].
In South Africa, gold mines are about 2800 m deep and the temperature of rock reaches 75 °C, whereas the Mponeng gold mine is about 4100 m deep and the temperature of rock reaches 66 °C. In Japan, the Chang-pan coal mines rock temperature reaches 42-48 °C and the water temperature reaches 72 °C. In the Feng-yu lead-zinc mine in Hokkaido, the temperature of rock reaches 69-130 °C. The average geothermal temperature per kilometer is 30-40 °C in Russia, with individual areas reaching 52 °C. In India, where gold mine are about 3000 m deep, the geothermal temperature reaches 70 °C. In China, some mines are one km deep, the rock temperature has exceeded 50 °C. In addition, many countries (e.g., Zambia, Mexico, Nicaragua, Poland, Germany) have thermal hazards in mines. However, the high geothermal temperature, which leads to the occurrence of thermal hazard, provides a good condition for geothermal energy exploitation and geothermal energ, as a clean and renewable energy, which belongs to the category of mining [9][10][11]. For mineral resources and geothermal energy to coexist in the deep layer, it is necessary to build respective mining systems. Based on this common characteristics, if geothermal energy can be simultaneously exploited with the help of the deep mineral resource systems and the collaborative exploitation of geothermal energy with mineral resources can be realized, then many advantages will be had [12].  It is unnecessary for drilling engineering to occur on the earth's surface, as it can save the cost of geothermal energy exploitation.  Deep space has a considerable scale. High temperature rocks provide a continuous heat source for geothermal energy exploitation.  Deep mining has complete elevation and transportation, as well as complete electric power, water supply, and drainage system, which provides pipeline layout safeguard and power supply for geothermal energy exploitation.  After mining is completed, the abandoned mine can be used as a geothermal energy plant.
By focusing on the collaborative exploitation of deep mineral resources and geothermal energy, the concept of phase-change heat storage backfill was put forward. The corresponding technical system was constructed and several parts involved in the technical system were analyzed. Aiming at the key of the technical system (I.E., the heat storage/heat release process of the backfill body), the heat transfer model of the phase-change backfill body was established and the temperature distribution and its evolution in the backfill body was numerically simulated and analyzed during the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy.
Deep space has a considerable scale. High temperature rocks provide a continuous heat source for geothermal energy exploitation.
Energies 2020, 13, x FOR PEER REVIEW 2 of 22 Backfill mining technology is an effective way to solve the problem of high ground stress, control the ground subsidence, and improve the mining rate of resources when mining deep mineral resources. Moreover, it does not cause serious damage to the ecological environment. With the rapid development of the mining industry, the situation of mineral resources is becoming more and more complex, and the geographical location and environmental conditions are also different for different mines. Therefore, it is necessary not only to select the corresponding backfill mining technology according to the specific mining conditions, but also to reasonably select the backfill materials in combination with the actual situation, so as to ensure its quality and performance can meet the mining operation need. The innovation and development of backfill materials explores a new way for resource reuse of solid wastes in mines, such as tailings, fly ash, and waste rocks. Under the premise of controlling ground subsidence, the solid waste of mines should be selected as backfill materials so backfill mining is also an effective means to make solid wastes harmless and reusable, and realize "green mining" in mines [3,4].
For deep mining, the ground temperature rises and heat releases to the stope from surrounding rocks; this is the main source of thermal hazard in the mine. In addition, other heat sources, such as air compression and mechanical equipment, lead to an increase of stope temperature. The deeper the mining depth, the higher the temperature of surrounding rock, and the more serious the thermal hazard. In many countries, thermal hazard in mines appears one after another [5][6][7][8].
In South Africa, gold mines are about 2800 m deep and the temperature of rock reaches 75 °C, whereas the Mponeng gold mine is about 4100 m deep and the temperature of rock reaches 66 °C. In Japan, the Chang-pan coal mines rock temperature reaches 42-48 °C and the water temperature reaches 72 °C. In the Feng-yu lead-zinc mine in Hokkaido, the temperature of rock reaches 69-130 °C. The average geothermal temperature per kilometer is 30-40 °C in Russia, with individual areas reaching 52 °C. In India, where gold mine are about 3000 m deep, the geothermal temperature reaches 70 °C. In China, some mines are one km deep, the rock temperature has exceeded 50 °C. In addition, many countries (e.g., Zambia, Mexico, Nicaragua, Poland, Germany) have thermal hazards in mines. However, the high geothermal temperature, which leads to the occurrence of thermal hazard, provides a good condition for geothermal energy exploitation and geothermal energ, as a clean and renewable energy, which belongs to the category of mining [9][10][11]. For mineral resources and geothermal energy to coexist in the deep layer, it is necessary to build respective mining systems. Based on this common characteristics, if geothermal energy can be simultaneously exploited with the help of the deep mineral resource systems and the collaborative exploitation of geothermal energy with mineral resources can be realized, then many advantages will be had [12].  It is unnecessary for drilling engineering to occur on the earth's surface, as it can save the cost of geothermal energy exploitation.  Deep space has a considerable scale. High temperature rocks provide a continuous heat source for geothermal energy exploitation.  Deep mining has complete elevation and transportation, as well as complete electric power, water supply, and drainage system, which provides pipeline layout safeguard and power supply for geothermal energy exploitation.  After mining is completed, the abandoned mine can be used as a geothermal energy plant.
By focusing on the collaborative exploitation of deep mineral resources and geothermal energy, the concept of phase-change heat storage backfill was put forward. The corresponding technical system was constructed and several parts involved in the technical system were analyzed. Aiming at the key of the technical system (I.E., the heat storage/heat release process of the backfill body), the heat transfer model of the phase-change backfill body was established and the temperature distribution and its evolution in the backfill body was numerically simulated and analyzed during the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy.
Deep mining has complete elevation and transportation, as well as complete electric power, water supply, and drainage system, which provides pipeline layout safeguard and power supply for geothermal energy exploitation. Backfill mining technology is an effective way to solve the problem of high ground stress, control the ground subsidence, and improve the mining rate of resources when mining deep mineral resources. Moreover, it does not cause serious damage to the ecological environment. With the rapid development of the mining industry, the situation of mineral resources is becoming more and more complex, and the geographical location and environmental conditions are also different for different mines. Therefore, it is necessary not only to select the corresponding backfill mining technology according to the specific mining conditions, but also to reasonably select the backfill materials in combination with the actual situation, so as to ensure its quality and performance can meet the mining operation need. The innovation and development of backfill materials explores a new way for resource reuse of solid wastes in mines, such as tailings, fly ash, and waste rocks. Under the premise of controlling ground subsidence, the solid waste of mines should be selected as backfill materials so backfill mining is also an effective means to make solid wastes harmless and reusable, and realize "green mining" in mines [3,4].
For deep mining, the ground temperature rises and heat releases to the stope from surrounding rocks; this is the main source of thermal hazard in the mine. In addition, other heat sources, such as air compression and mechanical equipment, lead to an increase of stope temperature. The deeper the mining depth, the higher the temperature of surrounding rock, and the more serious the thermal hazard. In many countries, thermal hazard in mines appears one after another [5][6][7][8].
In South Africa, gold mines are about 2800 m deep and the temperature of rock reaches 75 °C, whereas the Mponeng gold mine is about 4100 m deep and the temperature of rock reaches 66 °C. In Japan, the Chang-pan coal mines rock temperature reaches 42-48 °C and the water temperature reaches 72 °C. In the Feng-yu lead-zinc mine in Hokkaido, the temperature of rock reaches 69-130 °C. The average geothermal temperature per kilometer is 30-40 °C in Russia, with individual areas reaching 52 °C. In India, where gold mine are about 3000 m deep, the geothermal temperature reaches 70 °C. In China, some mines are one km deep, the rock temperature has exceeded 50 °C. In addition, many countries (e.g., Zambia, Mexico, Nicaragua, Poland, Germany) have thermal hazards in mines. However, the high geothermal temperature, which leads to the occurrence of thermal hazard, provides a good condition for geothermal energy exploitation and geothermal energ, as a clean and renewable energy, which belongs to the category of mining [9][10][11]. For mineral resources and geothermal energy to coexist in the deep layer, it is necessary to build respective mining systems. Based on this common characteristics, if geothermal energy can be simultaneously exploited with the help of the deep mineral resource systems and the collaborative exploitation of geothermal energy with mineral resources can be realized, then many advantages will be had [12].  It is unnecessary for drilling engineering to occur on the earth's surface, as it can save the cost of geothermal energy exploitation.  Deep space has a considerable scale. High temperature rocks provide a continuous heat source for geothermal energy exploitation.  Deep mining has complete elevation and transportation, as well as complete electric power, water supply, and drainage system, which provides pipeline layout safeguard and power supply for geothermal energy exploitation.  After mining is completed, the abandoned mine can be used as a geothermal energy plant.
By focusing on the collaborative exploitation of deep mineral resources and geothermal energy, the concept of phase-change heat storage backfill was put forward. The corresponding technical system was constructed and several parts involved in the technical system were analyzed. Aiming at the key of the technical system (I.E., the heat storage/heat release process of the backfill body), the heat transfer model of the phase-change backfill body was established and the temperature distribution and its evolution in the backfill body was numerically simulated and analyzed during the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy.
After mining is completed, the abandoned mine can be used as a geothermal energy plant.
By focusing on the collaborative exploitation of deep mineral resources and geothermal energy, the concept of phase-change heat storage backfill was put forward. The corresponding technical system was constructed and several parts involved in the technical system were analyzed. Aiming at the key of the technical system (I.E., the heat storage/heat release process of the backfill body), the heat transfer model of the phase-change backfill body was established and the temperature distribution and its evolution in the backfill body was numerically simulated and analyzed during the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy.

The Control of Thermal Hazard Development from Single Treatment to Combined with Geothermal Energy Exploitation
With the continuous extension of mineral resource exploitation into the deep layer, thermal hazards have become a major problem that restricts mining safety [5,13]. For thermal hazards, the control methods have gradually developed from initial ventilation, heat source control, personal protection, and other non-mechanical refrigeration to mechanical refrigeration technology, which can produce cold air, cold water, and cold ice to cooling the airflow in stope [5,14]. However, the mechanical refrigeration leads to an increase in mining cost and has high power consumption [15,16]. For the exploitation of mineral resources in the deep layer, the geothermal heat deteriorates working conditions and is a main thermal hazard heat source [14]. However, due to the problems of environmental pollution and energy shortage caused by excessive fossil fuel exploitation and consumption, geothermal energy has attracted attention as a clean and renewable green energy source [17,18]. The exploitation and utilization of geothermal energy can weaken the main thermal hazard heat source in deep mines, make passive prevention and control of thermal hazard gradually develop towards proactive geothermal energy exploitation combined with mine cooling, and provide a feasible technical way to prevent and control thermal hazards in deep mines.

The Extension and Integration from Deep Mineral Mining to Deep Geothermal Energy Exploitation
High geothermal temperature, high ground stress, rock physical property deterioration, mine lifting, and other problems in deep mines increases the cost of mines. However, the geothermal energy stored in the deep layer is a clean and renewable energy. Moreover, the reserve is very abundant, the advantage of environmental protection is obvious, the application field is wide, and it is a reliable, stable, and impervious to climatic conditions during exploitation and utilization [10,19,20]. The exploitation of geothermal energy requires a heat transfer fluid (cold fluid) to bring thermal energy to the ground. Therefore, it is necessary to drill well, reinforce well, and complete well when extracting thermal energy from the subsurface target layer. The depth of drilling well for extracting geothermal energy generally is between 1000 m and 3000 m. Most, however, are around 2000 m and a few are more than 3000 m. At present, the cost of drilling and completing well nearly accounts for the half of the total cost in geothermal energy exploitation [21]. With the development of exploitation technology, drilling well can extend to the deeper layer, which increases the cost. If the exploitation of geothermal energy can be combined with mineral resources in the deep layer, underground tunnels, complete elevation, transportation, power, water supply, and drainage systems can provide good conditions for piping layout, power supply, and the cycle of the heat transfer fluid in geothermal energy exploitation. In such a scenario, the cost of geothermal energy exploitation is be greatly reduced [22][23][24].

The Connotation of Phase-Change Heat Storage Backfill
Backfill mining technology has a long history in domestic and foreign metal mines, with the development of backfill mining theory and increasing requirements of deep mining. Backfill mining technology has developed from simple accumulative backfill (e.g., dry back-filling of waste rock, water-sand backfill) to structural backfill (e.g., cemented paste backfill), and then gradually transitioned to functional backfill [12]. So-called functional backfill means that the backfill body is formed by the solidification of backfill slurry and should not only have the basic function of realizing ground stress management and preventing ground subsidence, but also have some special function beyond transition, thus realizing the organic unification of each function. Based on the combination of controlling thermal hazards, exploiting geothermal energy, and extending from deep mineral mining to deep geothermal energy exploitation, that latent heat storage with phase-change materials is a promising technique to handle the intermittent issue of renewable energy [25][26][27][28]. This paper puts forth the concept of phase-change heat storage backfill ( Figure 1). When backfill mining is conducted, a certain dosage of PCM (phase-change materials) is added into backfill slurry under the condition of satisfying the requirements of mechanical properties for backfill mining. The backfill body is formed via solidified backfill slurry with good thermal properties. Geothermal energy can be exploited by heat storage and release from the backfill body during heat storage. The backfill body absorbs heat from surrounding rocks and stope airflow, reducing heat from surrounding rock to airflow, so it is realized of simultaneous stope cooling.
is conducted, a certain dosage of PCM (phase-change materials) is added into backfill slurry under the condition of satisfying the requirements of mechanical properties for backfill mining. The backfill body is formed via solidified backfill slurry with good thermal properties. Geothermal energy can be exploited by heat storage and release from the backfill body during heat storage. The backfill body absorbs heat from surrounding rocks and stope airflow, reducing heat from surrounding rock to airflow, so it is realized of simultaneous stope cooling.
As the core of backfill mining technology, backfill materials directly determine the quality of the backfill body, the cost of backfill mining, and the effect of controlling surface subsidence. Moreover, a phase-change heat storage backfill puts forth a new requirement for backfill materials, i.e., when a certain dosage of PCM is added into backfill material, the flow characteristics of backfill slurry must meet the requirements of pipeline transportation for backfill mining. The backfill body formed by solidified backfill slurry also should have good mechanical strength required by backfill mining and good thermal properties required by geothermal energy exploitation. Shape-stabilized PCMs can overcome negative influence caused by seepage or large volume changes during the solid-liquid and solid-gas phase-change. It does not need encapsulation as it increases the safety of material, reduces the heat transfer resistance of container, and is more conducive to heat transfer [29]. Shape-stabilized PCM is composed of a carrier substrate (supporting material) and phase-change medium. The carrier substrate can maintain the stability of the shape and physicochemical properties in the phase-change process of PCM, and has good thermal conductivity and mechanical properties, which can carry and limit the leakage of liquid PCMs [29,30]. The relevant research shows that tailings, fly ash, and waste rocks contained in backfill material can provide a carrier substrate for shape-stabilized PCMs [31,32], as they provide good conditions for implementing phase-change heat storage backfill technology.
According to the predetermined proportion, during the phase-change heat storage backfill process tailings, fly ash, waste rocks (or other aggregates), cementing materials, additive agents, and PCM are mixed with water and stirred to form a homogeneous backfill slurry. Then, the backfill slurry is transported via pipeline and poured into goaf, in which the cold fluid in the heat transfer tube is prearranged. Finally, the phase-change heat storage backfill body with inner buried tube is formed by solidified backfill slurry. The working process of the backfill body is shown in Figure 2. During the of heat storage process, the backfill body absorbs heat from surrounding rock and airflow in stopes. The temperature increases when the phase-change temperature is reached and As the core of backfill mining technology, backfill materials directly determine the quality of the backfill body, the cost of backfill mining, and the effect of controlling surface subsidence. Moreover, a phase-change heat storage backfill puts forth a new requirement for backfill materials, i.e., when a certain dosage of PCM is added into backfill material, the flow characteristics of backfill slurry must meet the requirements of pipeline transportation for backfill mining. The backfill body formed by solidified backfill slurry also should have good mechanical strength required by backfill mining and good thermal properties required by geothermal energy exploitation.
Shape-stabilized PCMs can overcome negative influence caused by seepage or large volume changes during the solid-liquid and solid-gas phase-change. It does not need encapsulation as it increases the safety of material, reduces the heat transfer resistance of container, and is more conducive to heat transfer [29]. Shape-stabilized PCM is composed of a carrier substrate (supporting material) and phase-change medium. The carrier substrate can maintain the stability of the shape and physicochemical properties in the phase-change process of PCM, and has good thermal conductivity and mechanical properties, which can carry and limit the leakage of liquid PCMs [29,30]. The relevant research shows that tailings, fly ash, and waste rocks contained in backfill material can provide a carrier substrate for shape-stabilized PCMs [31,32], as they provide good conditions for implementing phase-change heat storage backfill technology.
According to the predetermined proportion, during the phase-change heat storage backfill process tailings, fly ash, waste rocks (or other aggregates), cementing materials, additive agents, and PCM are mixed with water and stirred to form a homogeneous backfill slurry. Then, the backfill slurry is transported via pipeline and poured into goaf, in which the cold fluid in the heat transfer tube is prearranged. Finally, the phase-change heat storage backfill body with inner buried tube is formed by solidified backfill slurry. The working process of the backfill body is shown in Figure 2.
During the of heat storage process, the backfill body absorbs heat from surrounding rock and airflow in stopes. The temperature increases when the phase-change temperature is reached and PCM melts. Then, the temperature of the backfill body continually increases near surrounding rock. Heat is collected and stored in the backfill body by means of sensible and latent heat. In this case, due to the heat absorption of the backfill body from the surrounding rock and airflow, this reduced heat release and so stope cooling was simultaneously realized. During the heat release process, the heat transfer tube is filled with cold fluid, the backfill body releases the heat stored in itself to cold fluid, and the temperature decreases. PCM solidifies after losing heat and then the backfill body temperature continually decreases near the cold fluid. This provides the conditions of sustainable collection for geothermal energy.
Energies 2020, 13, x FOR PEER REVIEW 5 of 22 PCM melts. Then, the temperature of the backfill body continually increases near surrounding rock.
Heat is collected and stored in the backfill body by means of sensible and latent heat. In this case, due to the heat absorption of the backfill body from the surrounding rock and airflow, this reduced heat release and so stope cooling was simultaneously realized. During the heat release process, the heat transfer tube is filled with cold fluid, the backfill body releases the heat stored in itself to cold fluid, and the temperature decreases. PCM solidifies after losing heat and then the backfill body temperature continually decreases near the cold fluid. This provides the conditions of sustainable collection for geothermal energy.
(a) Heat storage process (b) Heat release process

The Construction of Technical System
Backfill mining involves pouring backfill slurry into a goaf. The backfill slurry then solidifies to form a cemented backfill body, which can control ground stress, as well as reduce, delay, and prevent the destruction and movement of surrounding rock after mining. Further, it can prevent surface subsidence. In preparation of backfill slurry, it is important to make use of solid wastes in order to realize green mining. Moreover, it is necessary to add a certain dosage of PCM into backfill material and to prearrange the heat transfer tube in the backfilled region. The backfill body has good thermal properties while meeting the requirements of mechanical strength for backfill mining. The

The Construction of Technical System
Backfill mining involves pouring backfill slurry into a goaf. The backfill slurry then solidifies to form a cemented backfill body, which can control ground stress, as well as reduce, delay, and prevent the destruction and movement of surrounding rock after mining. Further, it can prevent surface subsidence. In preparation of backfill slurry, it is important to make use of solid wastes in order to realize green mining. Moreover, it is necessary to add a certain dosage of PCM into backfill material Energies 2020, 13, 4755 6 of 20 and to prearrange the heat transfer tube in the backfilled region. The backfill body has good thermal properties while meeting the requirements of mechanical strength for backfill mining. The extraction and utilization of geothermal energy can be realized by heat storage/heat release of the backfill body and cold fluid circulation. Therefore, the technical system of phase-change heat storage backfill involves optimized backfill slurry proportions added with PCM. The collaborative optimization for mechanical and thermal properties of the backfill body has mutual cooperation for backfill mining, geothermal energy exploitation, and simultaneous stope cooling.
In different mines, the geological conditions and ore body characteristics are different. The characteristics of geothermal temperature and ground stress are different with mining depth. Therefore, the types of solid wastes, PCM selections, and requirements of mechanical properties for backfill body are different. The backfill body is placed in high ground stress and high ground temperature environment in the deep mine. The complexity of the geological and ore body occurring in the mine has a diversity of solid wastes and PCMs and ground temperature and pressure variability. Moreover, there is an expansion and contraction of volume, and change of pore structure and seepage flow that occurs in the backfill body. All of these directly affect each segment of the technical system for phase-change heat storage backfill. These influencing factors have a complex relationship, mutual influence, and mutual restraint with the proportion of backfill materials, the mechanical/thermal properties, and the heat storage/heat release characteristics of the backfill body. The technical system for the phase-change heat storage backfill is shown in Figure 3.
In different mines, the geological conditions and ore body characteristics are different. The characteristics of geothermal temperature and ground stress are different with mining depth. Therefore, the types of solid wastes, PCM selections, and requirements of mechanical properties for backfill body are different. The backfill body is placed in high ground stress and high ground temperature environment in the deep mine. The complexity of the geological and ore body occurring in the mine has a diversity of solid wastes and PCMs and ground temperature and pressure variability. Moreover, there is an expansion and contraction of volume, and change of pore structure and seepage flow that occurs in the backfill body. All of these directly affect each segment of the technical system for phase-change heat storage backfill. These influencing factors have a complex relationship, mutual influence, and mutual restraint with the proportion of backfill materials, the mechanical/thermal properties, and the heat storage/heat release characteristics of the backfill body. The technical system for the phase-change heat storage backfill is shown in Figure 3.

The optimized proportion of Backfill Slurry
Mining gradually develops to the deep layer. It is the goal of backfill mining technology to reduce backfill cost and ensure mining safety. The proportion of backfill slurry directly affects the mechanical and thermal properties of the backfill body, as well as the cost of backfill mining. Considering the specific conditions of various mines, making full use of solid wastes, selecting the heat storage materials with a suitable phase-change temperature, selecting for good comprehensive performance, determining the optimal proportion, preparing the backfill slurry with flow characteristics that meet the requirements of pipeline transportation, and forming a backfill body with mechanical strength that meets the requirements of backfill mining technology with good heat storage/heat release characteristics that meet the requirements of geothermal energy exploitation. All of the aforementioned qualities are important parts of phase-change heat storage backfill technology.

The Optimized Proportion of Backfill Slurry
Mining gradually develops to the deep layer. It is the goal of backfill mining technology to reduce backfill cost and ensure mining safety. The proportion of backfill slurry directly affects the mechanical and thermal properties of the backfill body, as well as the cost of backfill mining. Considering the specific conditions of various mines, making full use of solid wastes, selecting the heat storage materials with a suitable phase-change temperature, selecting for good comprehensive performance, determining the optimal proportion, preparing the backfill slurry with flow characteristics that meet the requirements of pipeline transportation, and forming a backfill body with mechanical strength that meets the requirements of backfill mining technology with good heat storage/heat release Energies 2020, 13, 4755 7 of 20 characteristics that meet the requirements of geothermal energy exploitation. All of the aforementioned qualities are important parts of phase-change heat storage backfill technology.
When preparing backfill slurry, solid wastes are selected-e.g., tailings, fly ash, waste rocks, coal gangue, slag, etc. [33,34]-and there are many types of heat storage materials available. Shape-stabilized PCMs are the main choice for heat storage materials, as their supporting materials mainly include cross linked high-polymer resin materials, porous inorganic materials, and inorganic nanomaterials. Phase-change mediums typically use solid-liquid PCMs with a high latent heat and nearly constant phase-change temperature, such as paraffin, fatty acids, high aliphatic hydrocarbons, esters, and polybasic alcohol, all of which have been widely utilized in thermal energy storage [35,36]. When heat storage materials are added into backfill materials, PCMs are required to have good thermal properties-such as suitable phase-change temperature, high latent heat, high thermal conductivity, large specific heat capacity, phase-change reversible, and repeatable cycle-as well as good environmental and economic performance [37] (i.e., no toxicity, no combustion, no pollution, convenient source, and a low cost). Various types of shape-stabilized PCMs, such as microcapsules, inorganic porous matrix, polymer matrix, nanometer matrix, and a functional additive type, can be obtained by synthesizing PCM supporting materials via different preparation methods [38,39].
To determine the optimal proportion of backfill slurry, an orthogonal or uniform test can be selected to conduct the design and optimization on proportion. The rheological parameters of backfill slurry (i.e., yield stress, shear stress, viscosity, slump, diffusivity, and precipitation) and the mechanical and thermal parameters of the backfill body (i.e., such as compressive strength, density, specific heat capacity, and thermal conductivity) were obtained by experimental testing. Based on the experimental results, the composition and proportion of backfill materials, the types and dosage of PCMs, cement-sand ratio, and slurry concentration on the rheological properties of backfill slurry can be analyzed and theoretically discussed. The optimization proportion model can be established for backfill slurry. The theoretical estimation method for thermal properties can be developed for backfill materials according to experimental data. Backfill slurry proportion optimization is shown in Figure 4. When preparing backfill slurry, solid wastes are selected-e.g., tailings, fly ash, waste rocks, coal gangue, slag, etc. [33,34]-and there are many types of heat storage materials available. Shape-stabilized PCMs are the main choice for heat storage materials, as their supporting materials mainly include cross linked high-polymer resin materials, porous inorganic materials, and inorganic nanomaterials. Phase-change mediums typically use solid-liquid PCMs with a high latent heat and nearly constant phase-change temperature, such as paraffin, fatty acids, high aliphatic hydrocarbons, esters, and polybasic alcohol, all of which have been widely utilized in thermal energy storage [35,36]. When heat storage materials are added into backfill materials, PCMs are required to have good thermal properties-such as suitable phase-change temperature, high latent heat, high thermal conductivity, large specific heat capacity, phase-change reversible, and repeatable cycle-as well as good environmental and economic performance [37] (i.e., no toxicity, no combustion, no pollution, convenient source, and a low cost). Various types of shape-stabilized PCMs, such as microcapsules, inorganic porous matrix, polymer matrix, nanometer matrix, and a functional additive type, can be obtained by synthesizing PCM supporting materials via different preparation methods [38,39].
To determine the optimal proportion of backfill slurry, an orthogonal or uniform test can be selected to conduct the design and optimization on proportion. The rheological parameters of backfill slurry (i.e., yield stress, shear stress, viscosity, slump, diffusivity, and precipitation) and the mechanical and thermal parameters of the backfill body (i.e., such as compressive strength, density, specific heat capacity, and thermal conductivity) were obtained by experimental testing. Based on the experimental results, the composition and proportion of backfill materials, the types and dosage of PCMs, cement-sand ratio, and slurry concentration on the rheological properties of backfill slurry can be analyzed and theoretically discussed. The optimization proportion model can be established for backfill slurry. The theoretical estimation method for thermal properties can be developed for backfill materials according to experimental data. Backfill slurry proportion optimization is shown in Figure 4.

The Collaborative Optimization on Mechanical and Thermal Properties of the Backfill Body
The technical system for a phase-change heat storage backfill should meet the requirements of mechanical properties for the backfill body (e.g., compressive strength) in order to control ground pressure, support surrounding rocks, and prevent surface subsidence. Yet it should also meet the

The Collaborative Optimization on Mechanical and Thermal Properties of the Backfill Body
The technical system for a phase-change heat storage backfill should meet the requirements of mechanical properties for the backfill body (e.g., compressive strength) in order to control ground pressure, support surrounding rocks, and prevent surface subsidence. Yet it should also meet the requirements of thermal properties for the backfill body (i.e., specific heat capacity, thermal conductivity, and enthalpy) in order to gather, store, and extract geothermal energy. Thermal properties of the backfill body can be improved by adding PCMs into backfill materials, while mechanical strength can be weakened when adding PCMs [40]. Therefore, it is important to test and analyze the mechanical and thermal parameters of the backfill body with different materials, proportions, and PCMs dosages to obtain the appropriate PCM dosage.
When preparing backfill slurry, solid wastes and PCMs can be selected. Moreover, due to the influence of backfill materials, PCM types and dosages, cement-sand ratio, and slurry concentration, phase-change heat storage backfill bodies can be formed via a solidified backfill slurry. As such, their mechanical and thermal properties are quite different as they meet the technical requirements of backfill mining and geothermal energy exploitation under various conditions. When the backfill body is placed in the high pressure and temperature deep mining environment, it not only needs to finish heat storage and heat release through temperature variation and a PCM phase-change cycle to realize the continuous exploitation and utilization of geothermal energy, but it also needs to maintain its mechanical strength to control ground pressure and support surrounding rocks during temperature variation and PCM phase-change. During the continuous of heat storage process, the backfill body temperature gradually increases and PCM absorbs heat and melts, leading to a decrease in the backfill body's mechanical strength, which affects the backfill body's supporting role on stope. Therefore, it is necessary to test and analyze the mechanical and thermal properties of the backfill body during its temperature variation and phase-change.
During the heat storage and heat release cycle under different temperature environments, the mechanical properties, thermal properties, and thermal stability of the backfill body with different materials and proportions need to be experimentally tested and theoretically analyzed. Moreover, a thermal-mechanical response model for phase-change heat storage also needs to be established. The characteristics and evolutionary rules of the temperature field and stress-strain field need to be revealed using numerical simulation, and the variation in mechanical properties, thermal properties and the coupling mechanism on the backfill body's thermal-mechanical properties needs to be obtained. Through our tests, simulations, and analyses, the backfill body can meet the requirements of mechanical properties for controlling ground pressure and maintaining good thermal properties for geothermal energy exploitation under thermal behavior. It can not only provide basic data for collaborative optimization on mechanical and thermal properties but also provide theoretical reference for the optimized proportion of backfill slurry. This is because the backfill body with good mechanical and thermal properties is an important objective of the optimized proportion of backfill slurry. The backfill body's collaborative optimization on mechanical and thermal properties is shown in Figure 4.

The Mutual Cooperation on Backfill Mining, Geothermal Energy Exploitation, and Simultaneous Stope Cooling
In the backfill mining, the heat transfer tube for geothermal energy exploitation is prearranged in the backfill body, as shown in Figure 2. During the heat storage process, the backfill body absorbs heat from surrounding rocks and stope airflow, the temperature increases, and the heat is stored in the backfill body via sensible heat and latent heat. In this case, stope cooling is simultaneously realized. During the heat release process, the cold fluid absorbs heat from a backfill body during temperature difference. The backfill body releases the stored heat to cold fluid and the temperature decreases. The backfill body accumulates and stores heat via the heat storage process and transfers the stored heat to cold fluid through the heat release process [41]. Cold fluid obtains heat and the temperature increases. Fluid heat can be extracted and utilized using a heat pump unit to realize the collaboration between backfill mining and geothermal energy exploitation.
The thermal behavior of the backfill body mainly includes the convective heat transfer with airflow and cold fluid, as well as heat conduction between the backfill body and surrounding rock, and heat Energies 2020, 13, 4755 9 of 20 conduction and phase-change heat transfer. Therefore, the stope airflow, geometric parameters, heat transfer tube arrangement, cold fluid temperature/flow characteristics, surrounding rock temperature, and the thermal properties of backfill material all have an influence on heat transfer. The heat transfer model for the backfill body should be constructed under the influence of multiple factors. Three-dimensional unsteady heat conduction, convective heat transfer, and phase-change heat transfer involved in them should be studied via numerical simulation and an experimental test. According to our simulation and test results, the heat transfer mechanism in the backfill body, surrounding rock, airflow, and cold fluid can be revealed. The influence of surrounding rock, airflow, cold fluid, and thermal properties on heat storage/heat release behavior can be analyzed and discussed. This can provide a theoretical basis when studying the thermal control system of "backfill mining-geothermal energy exploitation-simultaneous stope cooling". The mutual cooperation of the three aspects is shown in Figure 5. surrounding rock temperature, and the thermal properties of backfill material all have an influence on heat transfer. The heat transfer model for the backfill body should be constructed under the influence of multiple factors. Three-dimensional unsteady heat conduction, convective heat transfer, and phase-change heat transfer involved in them should be studied via numerical simulation and an experimental test. According to our simulation and test results, the heat transfer mechanism in the backfill body, surrounding rock, airflow, and cold fluid can be revealed. The influence of surrounding rock, airflow, cold fluid, and thermal properties on heat storage/heat release behavior can be analyzed and discussed. This can provide a theoretical basis when studying the thermal control system of "backfill mining-geothermal energy exploitation-simultaneous stope cooling". The mutual cooperation of the three aspects is shown in Figure 5.

The Heat Transfer Analysis of a Phase-Change Heat Storage Backfill
The heat transfer behavior of the backfill body plays a key role in geothermal energy exploitation and stope cooling. It is the core of the phase-change heat storage backfill, which includes heat transfer among the backfill body, surrounding rock, and airflow during the heat storage process and heat transfer between the backfill body and cold fluid. During the heat release process, the phase-change heat transfer is in the backfill body.

The Heat Storage/Heat Release Model of the backfill body
During backfill mining, backfill slurry mixed with PCM is transported to a deep stope through a pipeline, solidifies to form the backfill body, and the cold fluid in the heat transfer tube is prearranged inside the backfill body. In the case of stratified backfill mining, mining is done layer by layer, and the backfill also is done layer by layer. The layers of the backfill body are divided into three cases: the top layer adjoins with stope, its upper surface absorbs heat from airflow by convection heat transfer, the bottom layer adjoins with surrounding rock, several intermediate layers are located between the top and the bottom layers, the sides of all layers adjoin with the surrounding rock, and the bottom and the sides both can absorb heat from surrounding rock. In describing the heat transfer problem, the physical model is shown in Figure 6. During the heat storage process, the backfill body absorbs heat from surrounding rocks through heat conduction, which can reduce heat dissipation from surrounding rocks to airflow, and absorbs the heat from

The Heat Transfer Analysis of a Phase-Change Heat Storage Backfill
The heat transfer behavior of the backfill body plays a key role in geothermal energy exploitation and stope cooling. It is the core of the phase-change heat storage backfill, which includes heat transfer among the backfill body, surrounding rock, and airflow during the heat storage process and heat transfer between the backfill body and cold fluid. During the heat release process, the phase-change heat transfer is in the backfill body.

The Heat Storage/Heat Release Model of the Backfill Body
During backfill mining, backfill slurry mixed with PCM is transported to a deep stope through a pipeline, solidifies to form the backfill body, and the cold fluid in the heat transfer tube is prearranged inside the backfill body. In the case of stratified backfill mining, mining is done layer by layer, and the backfill also is done layer by layer. The layers of the backfill body are divided into three cases: the top layer adjoins with stope, its upper surface absorbs heat from airflow by convection heat transfer, the bottom layer adjoins with surrounding rock, several intermediate layers are located between the top and the bottom layers, the sides of all layers adjoin with the surrounding rock, and the bottom and the sides both can absorb heat from surrounding rock. In describing the heat transfer problem, Energies 2020, 13, 4755 10 of 20 the physical model is shown in Figure 6. During the heat storage process, the backfill body absorbs heat from surrounding rocks through heat conduction, which can reduce heat dissipation from surrounding rocks to airflow, and absorbs the heat from airflow in stope through convection heat transfer, so as to decrease airflow temperature of airflow, i.e., realize simultaneous stope cooling. During the heat release process, the backfill body transfers the stored heat to cold fluid. There are complex three-dimensional unsteady heat transfer processes among the backfill body, surrounding rocks, and airflow. Moreover, phase-change heat transfer exists in the backfill body. In order to simplify the calculation, the following assumptions are made: geothermal temperature per kilometer is 30-40 °C in Russia, with individual areas reaching 52 °C. In India, where gold mine are about 3000 m deep, the geothermal temperature reaches 70 °C. In China, some mines are one km deep, the rock temperature has exceeded 50 °C. In addition, many countries (e.g., Zambia, Mexico, Nicaragua, Poland, Germany) have thermal hazards in mines. However, the high geothermal temperature, which leads to the occurrence of thermal hazard, provides a good condition for geothermal energy exploitation and geothermal energ, as a clean and renewable energy, which belongs to the category of mining [9][10][11]. For mineral resources and geothermal energy to coexist in the deep layer, it is necessary to build respective mining systems. Based on this common characteristics, if geothermal energy can be simultaneously exploited with the help of the deep mineral resource systems and the collaborative exploitation of geothermal energy with mineral resources can be realized, then many advantages will be had [12].
 It is unnecessary for drilling engineering to occur on the earth's surface, as it can save the cost of geothermal energy exploitation.  Deep space has a considerable scale. High temperature rocks provide a continuous heat source for geothermal energy exploitation.  Deep mining has complete elevation and transportation, as well as complete electric power, water supply, and drainage system, which provides pipeline layout safeguard and power supply for geothermal energy exploitation.  After mining is completed, the abandoned mine can be used as a geothermal energy plant.
By focusing on the collaborative exploitation of deep mineral resources and geothermal energy, the concept of phase-change heat storage backfill was put forward. The corresponding technical system was constructed and several parts involved in the technical system were analyzed. Aiming at the key of the technical system (I.E., the heat storage/heat release process of the backfill body), the heat transfer model of the phase-change backfill body was established and the temperature distribution and its evolution in the backfill body was numerically simulated and analyzed during the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy.
The backfill body is a homogeneous and isotropic solid, and the thermophysical parameters remain constant. coal mines rock temperature reaches 42-48 °C and the water temperature reaches 72 °C. In the Feng-yu lead-zinc mine in Hokkaido, the temperature of rock reaches 69-130 °C. The average geothermal temperature per kilometer is 30-40 °C in Russia, with individual areas reaching 52 °C. In India, where gold mine are about 3000 m deep, the geothermal temperature reaches 70 °C. In China, some mines are one km deep, the rock temperature has exceeded 50 °C. In addition, many countries (e.g., Zambia, Mexico, Nicaragua, Poland, Germany) have thermal hazards in mines. However, the high geothermal temperature, which leads to the occurrence of thermal hazard, provides a good condition for geothermal energy exploitation and geothermal energ, as a clean and renewable energy, which belongs to the category of mining [9][10][11]. For mineral resources and geothermal energy to coexist in the deep layer, it is necessary to build respective mining systems. Based on this common characteristics, if geothermal energy can be simultaneously exploited with the help of the deep mineral resource systems and the collaborative exploitation of geothermal energy with mineral resources can be realized, then many advantages will be had [12].  It is unnecessary for drilling engineering to occur on the earth's surface, as it can save the cost of geothermal energy exploitation.  Deep space has a considerable scale. High temperature rocks provide a continuous heat source for geothermal energy exploitation.  Deep mining has complete elevation and transportation, as well as complete electric power, water supply, and drainage system, which provides pipeline layout safeguard and power supply for geothermal energy exploitation.  After mining is completed, the abandoned mine can be used as a geothermal energy plant.
By focusing on the collaborative exploitation of deep mineral resources and geothermal energy, the concept of phase-change heat storage backfill was put forward. The corresponding technical system was constructed and several parts involved in the technical system were analyzed. Aiming at the key of the technical system (I.E., the heat storage/heat release process of the backfill body), the heat transfer model of the phase-change backfill body was established and the temperature distribution and its evolution in the backfill body was numerically simulated and analyzed during the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy.
The backfill body is a porous medium with constant porosity, and all phases are in local thermodynamic equilibrium.
gold mines are about 2800 m deep and the temperature of rock reaches 75 °C, whereas the Mponeng gold mine is about 4100 m deep and the temperature of rock reaches 66 °C. In Japan, the Chang-pan coal mines rock temperature reaches 42-48 °C and the water temperature reaches 72 °C. In the Feng-yu lead-zinc mine in Hokkaido, the temperature of rock reaches 69-130 °C. The average geothermal temperature per kilometer is 30-40 °C in Russia, with individual areas reaching 52 °C. In India, where gold mine are about 3000 m deep, the geothermal temperature reaches 70 °C. In China, some mines are one km deep, the rock temperature has exceeded 50 °C. In addition, many countries (e.g., Zambia, Mexico, Nicaragua, Poland, Germany) have thermal hazards in mines. However, the high geothermal temperature, which leads to the occurrence of thermal hazard, provides a good condition for geothermal energy exploitation and geothermal energ, as a clean and renewable energy, which belongs to the category of mining [9][10][11]. For mineral resources and geothermal energy to coexist in the deep layer, it is necessary to build respective mining systems. Based on this common characteristics, if geothermal energy can be simultaneously exploited with the help of the deep mineral resource systems and the collaborative exploitation of geothermal energy with mineral resources can be realized, then many advantages will be had [12].  It is unnecessary for drilling engineering to occur on the earth's surface, as it can save the cost of geothermal energy exploitation.  Deep space has a considerable scale. High temperature rocks provide a continuous heat source for geothermal energy exploitation.  Deep mining has complete elevation and transportation, as well as complete electric power, water supply, and drainage system, which provides pipeline layout safeguard and power supply for geothermal energy exploitation.  After mining is completed, the abandoned mine can be used as a geothermal energy plant.
By focusing on the collaborative exploitation of deep mineral resources and geothermal energy, the concept of phase-change heat storage backfill was put forward. The corresponding technical system was constructed and several parts involved in the technical system were analyzed. Aiming at the key of the technical system (I.E., the heat storage/heat release process of the backfill body), the heat transfer model of the phase-change backfill body was established and the temperature distribution and its evolution in the backfill body was numerically simulated and analyzed during the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy.
During the heat storage process, the influences of the heat transfer tube and cold fluid are ignored.
In South Africa, gold mines are about 2800 m deep and the temperature of rock reaches 75 °C, whereas the Mponeng gold mine is about 4100 m deep and the temperature of rock reaches 66 °C. In Japan, the Chang-pan coal mines rock temperature reaches 42-48 °C and the water temperature reaches 72 °C. In the Feng-yu lead-zinc mine in Hokkaido, the temperature of rock reaches 69-130 °C. The average geothermal temperature per kilometer is 30-40 °C in Russia, with individual areas reaching 52 °C. In India, where gold mine are about 3000 m deep, the geothermal temperature reaches 70 °C. In China, some mines are one km deep, the rock temperature has exceeded 50 °C. In addition, many countries (e.g., Zambia, Mexico, Nicaragua, Poland, Germany) have thermal hazards in mines. However, the high geothermal temperature, which leads to the occurrence of thermal hazard, provides a good condition for geothermal energy exploitation and geothermal energ, as a clean and renewable energy, which belongs to the category of mining [9][10][11]. For mineral resources and geothermal energy to coexist in the deep layer, it is necessary to build respective mining systems. Based on this common characteristics, if geothermal energy can be simultaneously exploited with the help of the deep mineral resource systems and the collaborative exploitation of geothermal energy with mineral resources can be realized, then many advantages will be had [12].  It is unnecessary for drilling engineering to occur on the earth's surface, as it can save the cost of geothermal energy exploitation.  Deep space has a considerable scale. High temperature rocks provide a continuous heat source for geothermal energy exploitation.  Deep mining has complete elevation and transportation, as well as complete electric power, water supply, and drainage system, which provides pipeline layout safeguard and power supply for geothermal energy exploitation.  After mining is completed, the abandoned mine can be used as a geothermal energy plant.
By focusing on the collaborative exploitation of deep mineral resources and geothermal energy, the concept of phase-change heat storage backfill was put forward. The corresponding technical system was constructed and several parts involved in the technical system were analyzed. Aiming at the key of the technical system (I.E., the heat storage/heat release process of the backfill body), the heat transfer model of the phase-change backfill body was established and the temperature distribution and its evolution in the backfill body was numerically simulated and analyzed during the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy.
During the heat release process, there is no heat transfer in the backfill body to surrounding rocks and airflow.
air compression and mechanical equipment, lead to an increase of stope temperature. The deeper the mining depth, the higher the temperature of surrounding rock, and the more serious the thermal hazard. In many countries, thermal hazard in mines appears one after another [5][6][7][8].
In South Africa, gold mines are about 2800 m deep and the temperature of rock reaches 75 °C, whereas the Mponeng gold mine is about 4100 m deep and the temperature of rock reaches 66 °C. In Japan, the Chang-pan coal mines rock temperature reaches 42-48 °C and the water temperature reaches 72 °C. In the Feng-yu lead-zinc mine in Hokkaido, the temperature of rock reaches 69-130 °C. The average geothermal temperature per kilometer is 30-40 °C in Russia, with individual areas reaching 52 °C. In India, where gold mine are about 3000 m deep, the geothermal temperature reaches 70 °C. In China, some mines are one km deep, the rock temperature has exceeded 50 °C. In addition, many countries (e.g., Zambia, Mexico, Nicaragua, Poland, Germany) have thermal hazards in mines. However, the high geothermal temperature, which leads to the occurrence of thermal hazard, provides a good condition for geothermal energy exploitation and geothermal energ, as a clean and renewable energy, which belongs to the category of mining [9][10][11]. For mineral resources and geothermal energy to coexist in the deep layer, it is necessary to build respective mining systems. Based on this common characteristics, if geothermal energy can be simultaneously exploited with the help of the deep mineral resource systems and the collaborative exploitation of geothermal energy with mineral resources can be realized, then many advantages will be had [12].  It is unnecessary for drilling engineering to occur on the earth's surface, as it can save the cost of geothermal energy exploitation.  Deep space has a considerable scale. High temperature rocks provide a continuous heat source for geothermal energy exploitation.  Deep mining has complete elevation and transportation, as well as complete electric power, water supply, and drainage system, which provides pipeline layout safeguard and power supply for geothermal energy exploitation.  After mining is completed, the abandoned mine can be used as a geothermal energy plant.
By focusing on the collaborative exploitation of deep mineral resources and geothermal energy, the concept of phase-change heat storage backfill was put forward. The corresponding technical system was constructed and several parts involved in the technical system were analyzed. Aiming at the key of the technical system (I.E., the heat storage/heat release process of the backfill body), the heat transfer model of the phase-change backfill body was established and the temperature distribution and its evolution in the backfill body was numerically simulated and analyzed during the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy.
During the heat release process, the heat transfer along the axial direction of the heat transfer tube is ignored, and the heat transfer of the backfill body is a two-dimensional unsteady heat conduction.
For deep mining, the ground temperature rises and heat releases to the stope from surrounding rocks; this is the main source of thermal hazard in the mine. In addition, other heat sources, such as air compression and mechanical equipment, lead to an increase of stope temperature. The deeper the mining depth, the higher the temperature of surrounding rock, and the more serious the thermal hazard. In many countries, thermal hazard in mines appears one after another [5][6][7][8].
In South Africa, gold mines are about 2800 m deep and the temperature of rock reaches 75 °C, whereas the Mponeng gold mine is about 4100 m deep and the temperature of rock reaches 66 °C. In Japan, the Chang-pan coal mines rock temperature reaches 42-48 °C and the water temperature reaches 72 °C. In the Feng-yu lead-zinc mine in Hokkaido, the temperature of rock reaches 69-130 °C. The average geothermal temperature per kilometer is 30-40 °C in Russia, with individual areas reaching 52 °C. In India, where gold mine are about 3000 m deep, the geothermal temperature reaches 70 °C. In China, some mines are one km deep, the rock temperature has exceeded 50 °C. In addition, many countries (e.g., Zambia, Mexico, Nicaragua, Poland, Germany) have thermal hazards in mines. However, the high geothermal temperature, which leads to the occurrence of thermal hazard, provides a good condition for geothermal energy exploitation and geothermal energ, as a clean and renewable energy, which belongs to the category of mining [9][10][11]. For mineral resources and geothermal energy to coexist in the deep layer, it is necessary to build respective mining systems. Based on this common characteristics, if geothermal energy can be simultaneously exploited with the help of the deep mineral resource systems and the collaborative exploitation of geothermal energy with mineral resources can be realized, then many advantages will be had [12].  It is unnecessary for drilling engineering to occur on the earth's surface, as it can save the cost of geothermal energy exploitation.  Deep space has a considerable scale. High temperature rocks provide a continuous heat source for geothermal energy exploitation.  Deep mining has complete elevation and transportation, as well as complete electric power, water supply, and drainage system, which provides pipeline layout safeguard and power supply for geothermal energy exploitation.  After mining is completed, the abandoned mine can be used as a geothermal energy plant.
By focusing on the collaborative exploitation of deep mineral resources and geothermal energy, the concept of phase-change heat storage backfill was put forward. The corresponding technical system was constructed and several parts involved in the technical system were analyzed. Aiming at the key of the technical system (I.E., the heat storage/heat release process of the backfill body), the heat transfer model of the phase-change backfill body was established and the temperature distribution and its evolution in the backfill body was numerically simulated and analyzed during the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy.
The contact thermal resistance between the backfill body and heat transfer tube is ignored. and realize "green mining" in mines [3,4].
For deep mining, the ground temperature rises and heat releases to the stope from surrounding rocks; this is the main source of thermal hazard in the mine. In addition, other heat sources, such as air compression and mechanical equipment, lead to an increase of stope temperature. The deeper the mining depth, the higher the temperature of surrounding rock, and the more serious the thermal hazard. In many countries, thermal hazard in mines appears one after another [5][6][7][8].
In South Africa, gold mines are about 2800 m deep and the temperature of rock reaches 75 °C, whereas the Mponeng gold mine is about 4100 m deep and the temperature of rock reaches 66 °C. In Japan, the Chang-pan coal mines rock temperature reaches 42-48 °C and the water temperature reaches 72 °C. In the Feng-yu lead-zinc mine in Hokkaido, the temperature of rock reaches 69-130 °C. The average geothermal temperature per kilometer is 30-40 °C in Russia, with individual areas reaching 52 °C. In India, where gold mine are about 3000 m deep, the geothermal temperature reaches 70 °C. In China, some mines are one km deep, the rock temperature has exceeded 50 °C. In addition, many countries (e.g., Zambia, Mexico, Nicaragua, Poland, Germany) have thermal hazards in mines. However, the high geothermal temperature, which leads to the occurrence of thermal hazard, provides a good condition for geothermal energy exploitation and geothermal energ, as a clean and renewable energy, which belongs to the category of mining [9][10][11]. For mineral resources and geothermal energy to coexist in the deep layer, it is necessary to build respective mining systems. Based on this common characteristics, if geothermal energy can be simultaneously exploited with the help of the deep mineral resource systems and the collaborative exploitation of geothermal energy with mineral resources can be realized, then many advantages will be had [12].  It is unnecessary for drilling engineering to occur on the earth's surface, as it can save the cost of geothermal energy exploitation.  Deep space has a considerable scale. High temperature rocks provide a continuous heat source for geothermal energy exploitation.  Deep mining has complete elevation and transportation, as well as complete electric power, water supply, and drainage system, which provides pipeline layout safeguard and power supply for geothermal energy exploitation.  After mining is completed, the abandoned mine can be used as a geothermal energy plant.
By focusing on the collaborative exploitation of deep mineral resources and geothermal energy, the concept of phase-change heat storage backfill was put forward. The corresponding technical system was constructed and several parts involved in the technical system were analyzed. Aiming at the key of the technical system (I.E., the heat storage/heat release process of the backfill body), the heat transfer model of the phase-change backfill body was established and the temperature distribution and its evolution in the backfill body was numerically simulated and analyzed during the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy.
In a given temperature range, PCM only undergoes a solid-liquid phase-change, without supercooling and performance degradation. the premise of controlling ground subsidence, the solid waste of mines should be selected as backfill materials so backfill mining is also an effective means to make solid wastes harmless and reusable, and realize "green mining" in mines [3,4].
For deep mining, the ground temperature rises and heat releases to the stope from surrounding rocks; this is the main source of thermal hazard in the mine. In addition, other heat sources, such as air compression and mechanical equipment, lead to an increase of stope temperature. The deeper the mining depth, the higher the temperature of surrounding rock, and the more serious the thermal hazard. In many countries, thermal hazard in mines appears one after another [5][6][7][8].
In South Africa, gold mines are about 2800 m deep and the temperature of rock reaches 75 °C, whereas the Mponeng gold mine is about 4100 m deep and the temperature of rock reaches 66 °C. In Japan, the Chang-pan coal mines rock temperature reaches 42-48 °C and the water temperature reaches 72 °C. In the Feng-yu lead-zinc mine in Hokkaido, the temperature of rock reaches 69-130 °C. The average geothermal temperature per kilometer is 30-40 °C in Russia, with individual areas reaching 52 °C. In India, where gold mine are about 3000 m deep, the geothermal temperature reaches 70 °C. In China, some mines are one km deep, the rock temperature has exceeded 50 °C. In addition, many countries (e.g., Zambia, Mexico, Nicaragua, Poland, Germany) have thermal hazards in mines. However, the high geothermal temperature, which leads to the occurrence of thermal hazard, provides a good condition for geothermal energy exploitation and geothermal energ, as a clean and renewable energy, which belongs to the category of mining [9][10][11]. For mineral resources and geothermal energy to coexist in the deep layer, it is necessary to build respective mining systems. Based on this common characteristics, if geothermal energy can be simultaneously exploited with the help of the deep mineral resource systems and the collaborative exploitation of geothermal energy with mineral resources can be realized, then many advantages will be had [12].  It is unnecessary for drilling engineering to occur on the earth's surface, as it can save the cost of geothermal energy exploitation.  Deep space has a considerable scale. High temperature rocks provide a continuous heat source for geothermal energy exploitation.  Deep mining has complete elevation and transportation, as well as complete electric power, water supply, and drainage system, which provides pipeline layout safeguard and power supply for geothermal energy exploitation.  After mining is completed, the abandoned mine can be used as a geothermal energy plant.
By focusing on the collaborative exploitation of deep mineral resources and geothermal energy, the concept of phase-change heat storage backfill was put forward. The corresponding technical system was constructed and several parts involved in the technical system were analyzed. Aiming at the key of the technical system (I.E., the heat storage/heat release process of the backfill body), the heat transfer model of the phase-change backfill body was established and the temperature distribution and its evolution in the backfill body was numerically simulated and analyzed during the heat storage process/heat release. This study can provide the theoretical framework for phase-change heat storage backfill technology and technical guidance for the collaborative exploitation of deep mineral resources and geothermal energy.
We did not consider the natural convection heat transfer of liquid PCM.
Energies 2020, 13, x FOR PEER REVIEW 10 of 22 airflow in stope through convection heat transfer, so as to decrease airflow temperature of airflow, i.e., realize simultaneous stope cooling. During the heat release process, the backfill body transfers the stored heat to cold fluid. There are complex three-dimensional unsteady heat transfer processes among the backfill body, surrounding rocks, and airflow. Moreover, phase-change heat transfer exists in the backfill body. In order to simplify the calculation, the following assumptions are made:  The backfill body is a homogeneous and isotropic solid, and the thermophysical parameters remain constant.  The backfill body is a porous medium with constant porosity, and all phases are in local thermodynamic equilibrium.  During the heat storage process, the influences of the heat transfer tube and cold fluid are ignored.  During the heat release process, there is no heat transfer in the backfill body to surrounding rocks and airflow.  During the heat release process, the heat transfer along the axial direction of the heat transfer tube is ignored, and the heat transfer of the backfill body is a two-dimensional unsteady heat conduction.  The contact thermal resistance between the backfill body and heat transfer tube is ignored.  In a given temperature range, PCM only undergoes a solid-liquid phase-change, without super-cooling and performance degradation.  We did not consider the natural convection heat transfer of liquid PCM. The control equation for the backfill body during the heat storage process is shown as follows: The control equation for the backfill body during the heat storage process is shown as follows: Energies 2020, 13, 4755 11 of 20 the initial condition: T(x, y, z, τ)| τ=0 = T ini the boundary condition of surrounding rocks: T s (x, y, z, τ) = T b the boundary condition of stope: −λ ∂T where T is the temperature ( • C); H is the total enthalpy equal to the sum of sensible and latent heat in the backfill body (kJ/kg); λ is the thermal conductivity (W/(m. • C)); ρ is the density (kg/m 3 ); T ini is the initial temperature of the backfill body at the beginning of heat storage ( • C); h is the convection heat transfer coefficient (W/(m 2 . • C)). The subscript s expresses the boundary between the backfill body and surrounding rocks, w expresses the boundary between the backfill body and stope, and af expresses airflow. n is the normal direction of w; T b is the temperature of surrounding rock at the boundary s ( • C); T w is the temperature of the backfill body at the boundary w, ( • C).

The Mathematical Model of the Backfill Body during the Heat Release Process
For the complex three-dimensional unsteady heat transfer between the backfill body and cold fluid, the heat transfer models at different directions can solve the problem. The axial one-dimensional heat transfer model is established for cold fluid and the two-dimensional heat transfer model perpendicular to the heat transfer tube is established for the backfill body for the acting region of cold fluid. The heat transfer between cold fluid at the axial direction of the heat transfer tube and its surrounding backfill body can be calculated. The heat transfer coupling between the two regions can be realized.
(1) The axial one-dimensional heat transfer model for cold fluid.
We take the micro-segment of the heat transfer tube, ∆x, as the object of study, according to the energy balance: where q x is the heat transfer capacity of unit length heat transfer tube (W/m); m is the mass flow of cold fluid (kg/s); and c is the specific heat capacity at constant pressure (J/(kg. • C)). The subscript cf expresses cold fluid; in expresses inlet; out expresses outlet.
(2) The two-dimensional heat transfer model on the plane perpendicular to the heat transfer tube for the backfill body.
the initial condition: T c f = T p = T b f = T ini the boundary condition: (1) the boundary condition of the cold fluid's acting region: q l = 0 (2) the boundary condition of the contact surface for the outer wall of the heat transfer tube and backfill body: the boundary condition of the cold fluid inlet: the boundary condition of the contact surface for the heat transfer tube's inner wall and cold fluid: where T ini is the initial temperature of the backfill body at the beginning of heat release ( • C); q is the heat flux density (W/m 2 ); n is the normal direction of the contact surface for the outer wall of the heat transfer tube and backfill body. The subscript bf expresses the backfill body, p expresses the wall of the heat transfer tube, and l expresses the acting region boundary of cold fluid on the plane perpendicular to the heat transfer tube.

The Phase-Change Heat Transfer Model for the Backfill Body
In this study, the backfill body formed by a solidified backfill slurry, which was added with a certain dosage of PCM and uniformly mixed. It was regarded as a porous medium. In the numerical simulation, the phase-change mechanism of the backfill body during the heat storage process/heat release was analyzed based on the enthalpy-porosity method. The heat balance between PCM and a solid aggregate of the backfill body was divided into local thermal equilibrium and non-local thermal equilibrium. Based on the assumptions made in Section 4.1, the natural convection heat transfer of liquid PCM was not considered, so local thermal equilibrium model was adopted. The expression is shown as follows [42,43]: The parameters in Equation (4) are expressed as follows: where ε is the porosity of the backfill body and is the volume fraction of PCM, %; L f is the latent heat of PCM, kJ/kg. The subscript eff expresses equivalent; bfs expresses the solid aggregate of the backfill body; pcm expresses phase-change material; l expresses the liquid phase of PCM; and s expresses the solid phase of PCM. γ is the liquid fraction of PCM, wherein the value depends on the temperature of the backfill body and can be expressed as follows: where T s is the phase-change temperature of solid PCM, • C; T l is the phase-change temperature of liquid PCM, • C.
The melting and solidification model of PCM is shown as follows: The continuity equation: The energy equation: where H is the enthalpy of PCM, i.e., the sum of sensible heat enthalpy and latent heat enthalpy ∆H, kJ/kg; v i is the flow velocity of liquid phase when phase-change occurs, m/s; → u is velocity vector.

Grid Generation and Independence Verification
According to the heat transfer characteristics of the backfill body, a three dimensional model was established using Solidworks Software and a grid was generated using ICEM Software (Figure 7). In the whole grid generation, the structural hexahedron was adopted for the backfill body and cold Energies 2020, 13, 4755 13 of 20 fluid. In order to reduce calculation error and enhance the heat transfer simulation between cold fluid and the backfill body, the cold fluid was divided by the inner O block, and the wall of the heat transfer tube and U-tube bending were encrypted. The simulation result was feasible when the U-tube bending was treated according to adiabatic condition, i.e., the energy change was ignored [44].

Grid Generation and Independence Verification
According to the heat transfer characteristics of the backfill body, a three dimensional model was established using Solidworks Software and a grid was generated using ICEM Software ( Figure  7). In the whole grid generation, the structural hexahedron was adopted for the backfill body and cold fluid. In order to reduce calculation error and enhance the heat transfer simulation between cold fluid and the backfill body, the cold fluid was divided by the inner O block, and the wall of the heat transfer tube and U-tube bending were encrypted. The simulation result was feasible when the U-tube bending was treated according to adiabatic condition, i.e., the energy change was ignored [44].
. In order to determine the optimal grid number and time step for numerical simulation, the average temperature of the backfill body and cold fluid outlet temperature were taken as the reference for heat storage and release, respectively. Numerical simulation was conducted over 5 h. During the heat storage process, grid numbers were 0.1296, 0.2176, and 0.3563 million. During the heat release process, grid numbers were 0.2915, 0.4656, and 0.6886 million (Figure 8a and Figure 9a). Time step independence was determined by taking 1 s, 5 s, and 10 s, respectively (Figure 8b and Figure 9b). As can be seen, the grid number and time step had little effect on the simulation results, as the average temperature of the backfill body and cold fluid outlet temperature had no obvious difference in grid numbers or time steps. To ensure simulation result accuracy, we reduced calculation workload and improved calculation efficiency. Thus, grid numbers of 0.2176 million and 0.4656 million were adopted for heat storage and heat release. Further, a time step of 5 s were adopted for simulation. In order to determine the optimal grid number and time step for numerical simulation, the average temperature of the backfill body and cold fluid outlet temperature were taken as the reference for heat storage and release, respectively. Numerical simulation was conducted over 5 h. During the heat storage process, grid numbers were 0.1296, 0.2176, and 0.3563 million. During the heat release process, grid numbers were 0.2915, 0.4656, and 0.6886 million (Figures 8a and 9a). Time step independence was determined by taking 1 s, 5 s, and 10 s, respectively (Figures 8b and 9b). As can be seen, the grid number and time step had little effect on the simulation results, as the average temperature of the backfill body and cold fluid outlet temperature had no obvious difference in grid numbers or time steps. To ensure simulation result accuracy, we reduced calculation workload and improved calculation efficiency. Thus, grid numbers of 0.2176 million and 0.4656 million were adopted for heat storage and heat release. Further, a time step of 5 s were adopted for simulation.

Setting Simulation Conditions
As shown in Figure 7, the size of the backfill body was 500 mm × 400 mm × 130 mm. Tailings of a metal mine were used as backfill aggregate. Paraffin microcapsules were selected as heat storage materials and uniformly mixed with backfill slurry (slurry concentration = 70%; cement-sand ratio of 1:4; paraffin percentage = 5% (the volume fraction)). The thermophysical parameters of the backfill body were measured after curing for 28 days. The inner diameter of the heat transfer tube was 4 mm and the distance between the adjacent tube centers was 80 mm. The initial temperature of the backfill body at the beginning of heat storage was 18 °C. The temperature of surrounding rocks at the boundary of the backfill body was 35 °C. The temperature was 26 °C and velocity was 2 m/s for airflow in stope. The calculated convective heat transfer coefficient on the surface of the backfill body near stope was 11.104 W/(m 2 ⋅°C). The initial temperature of the backfill body at the beginning of heat release was 35 °C. The inlet temperature was 18 °C and the inlet velocity was 0.2 m/s for cold fluid. The thermophysical parameters of the backfill body, PCM, cold fluid, and heat transfer tube are shown in Table 1.

Setting Simulation Conditions
As shown in Figure 7, the size of the backfill body was 500 mm × 400 mm × 130 mm. Tailings of a metal mine were used as backfill aggregate. Paraffin microcapsules were selected as heat storage materials and uniformly mixed with backfill slurry (slurry concentration = 70%; cement-sand ratio of 1:4; paraffin percentage = 5% (the volume fraction)). The thermophysical parameters of the backfill body were measured after curing for 28 days. The inner diameter of the heat transfer tube was 4 mm and the distance between the adjacent tube centers was 80 mm. The initial temperature of the backfill body at the beginning of heat storage was 18 °C. The temperature of surrounding rocks at the boundary of the backfill body was 35 °C. The temperature was 26 °C and velocity was 2 m/s for airflow in stope. The calculated convective heat transfer coefficient on the surface of the backfill body near stope was 11.104 W/(m 2 ⋅°C). The initial temperature of the backfill body at the beginning of heat release was 35 °C. The inlet temperature was 18 °C and the inlet velocity was 0.2 m/s for cold fluid. The thermophysical parameters of the backfill body, PCM, cold fluid, and heat transfer tube are shown in Table 1.

Setting Simulation Conditions
As shown in Figure 7, the size of the backfill body was 500 mm × 400 mm × 130 mm. Tailings of a metal mine were used as backfill aggregate. Paraffin microcapsules were selected as heat storage materials and uniformly mixed with backfill slurry (slurry concentration = 70%; cement-sand ratio of 1:4; paraffin percentage = 5% (the volume fraction)). The thermophysical parameters of the backfill body were measured after curing for 28 days. The inner diameter of the heat transfer tube was 4 mm and the distance between the adjacent tube centers was 80 mm. The initial temperature of the backfill body at the beginning of heat storage was 18 • C. The temperature of surrounding rocks at the boundary of the backfill body was 35 • C. The temperature was 26 • C and velocity was 2 m/s for airflow in stope. The calculated convective heat transfer coefficient on the surface of the backfill body near stope was 11.104 W/(m 2 · • C). The initial temperature of the backfill body at the beginning of heat release was 35 • C. The inlet temperature was 18 • C and the inlet velocity was 0.2 m/s for cold fluid. The thermophysical parameters of the backfill body, PCM, cold fluid, and heat transfer tube are shown in Table 1.  Figure 10 shows the temperature distribution on section X = 0 of the backfill body at different heat storage times. It shows that the higher temperature region gradually expanded from the surrounding rocks to the interior of the backfill body. The expansion was faster in the early stage and gradually slowed in the later stage, leading to a continuous reduction in the lower temperature region. The temperature region of T bf = 18 • C at the direction of Z reduced from 300 mm to 220 mm within 20 min-1 h (the decrement was 80 mm within 40 min). It reduced to 140 mm within 1 to 2 h (the decrement was 80 mm within 1 h) and then almost stayed at 100 mm within the following 3-5 h. The temperature region of T bf = 18 • C at the direction of Y also reduced and gradually moved toward stope. It reduced from 70 mm to 40 mm within 20 min-1 h (the decrement was 30 mm within 40 min) and reduced to 30 mm within 1-2 h (the decrement is 10 mm within 1 h). In addition, there was a temperature region rise on the side near the stope. The backfill body absorbed heat from the airflow. Within the following 3-5 h, the temperature region T bf = 18 • C at the direction of Y reduced from 40 mm to 25 mm (the decrement is 15 mm within 2 h), and the temperature region rise on the side near the stope disappeared, i.e., the backfill body no longer absorbed heat from airflow. Moreover, the temperature of the backfill body near the surrounding rocks was always higher than the adjacent stope. The influence of airflow on heat storage was less than the surrounding rock. Therefore, the heat transfer of surrounding rock was the main heat storage mode for the backfill body. The closer to the surrounding rock, the more obvious the temperature rise of the backfill body.    In the early stage of heat storage, the temperature difference of the backfill body with surrounding rocks and airflow was larger, and the heat transfer rate was higher. However, during the heat storage process, backfill body temperature gradually increased, temperature difference became smaller, and heat transfer rate gradually decreased. Therefore, the expansion of a higher temperature region was faster in the early stage and became slower in the late stage. Moreover, when the temperature of the backfill body was higher than the airflow, the backfill body no longer absorbed heat from airflow, resulting in the disappearance of the temperature rise region on the side near the stope. Figure 11 shows the temperature distribution for section Y = 0 from the backfill body at different heat release times. The temperature of the backfill body in the heat transfer tube gradually decreased and the lower temperature region gradually expanded from the heat transfer tube to the periphery during heat release. Expansion was faster in the early stage and gradually slowed down in the late stage due to a temperature difference. The heat transfer rate was larger in the early stage, then gradually decreased. Figure 11 (Point O) shows that the temperature was about 34.02 • C, 31.93 • C, 30.17 • C, and 29.09 • C at the heat release of 20 min, 1 h, 2 h, and 3 h, respectively. At the inlet of cold fluid, the temperature region of T bf = 27.55 • C around the heat transfer tube expanded from 1.5 mm to 18 mm within 20 min-2 h (the increment was 16.5 mm within 100 min), and expanded from 18 mm to 27 mm within 2 h-5 h (the increment is 9 mm within 180 min). the heat release of 2 h (Figure 11c), the width of the temperature region of Tbf =30 °C between tubes (the yellow region) was 10 mm, 18 mm, and 25 mm along the cold fluid's flow direction. The smaller the width of this region, the more obvious the heat release interaction between tubes was. Until the heat release of 5 h, the temperature of most regions within the influence range was close to each other (about 27 °C). There was still a temperature difference between the backfill body and cold fluid, so the heat release continued.

Conclusions
(1) We investigated the technical system of a phase-change heat storage backfill. The main content involved an optimized proportion of backfill slurry added with phase-change materials (PCMs). There was collaborative optimization on the mechanical and thermal properties of the backfill body, and mutual cooperation on backfill mining, geothermal energy exploitation, and simultaneous stope cooling.
(2) The optimized proportion made backfill slurry with flow characteristics that met pipeline transportation requirements. The backfill body was formed by solidified backfill slurry and its mechanical strength met the requirements of backfill mining technology, with good heat storage/heat release characteristics that met the requirements of geothermal energy exploitation. The heat transfer behavior of the backfill body plays a key role in geothermal energy exploitation and Figure 11. The temperature distribution of the backfill body during the heat release process. Figure 11 shows that the lower temperature region was larger around the middle tubes than around the inlet and outlet tube at the same heat release time. This phenomenon gradually became obvious during the heat release process. At the heat release of 2 h, the temperature region of T bf = 28.5 • C around the heat transfer tube was 17 mm, 21 mm, 18 mm, and 13 mm along the cold fluid's flow direction. This was because the backfill body temperature in the middle tubes was affected by the cold fluid on both sides and there was an interactive heat release region. Thus, the lower temperature Energies 2020, 13, 4755 18 of 20 region around the middle tubes was larger. Moreover, the closer to the inlet the lower the cold fluid temperature, and the more obvious the heat release interaction was between tubes. At the heat release of 2 h (Figure 11c), the width of the temperature region of T bf = 30 • C between tubes (the yellow region) was 10 mm, 18 mm, and 25 mm along the cold fluid's flow direction. The smaller the width of this region, the more obvious the heat release interaction between tubes was. Until the heat release of 5 h, the temperature of most regions within the influence range was close to each other (about 27 • C). There was still a temperature difference between the backfill body and cold fluid, so the heat release continued.

Conclusions
(1) We investigated the technical system of a phase-change heat storage backfill. The main content involved an optimized proportion of backfill slurry added with phase-change materials (PCMs). There was collaborative optimization on the mechanical and thermal properties of the backfill body, and mutual cooperation on backfill mining, geothermal energy exploitation, and simultaneous stope cooling.
(2) The optimized proportion made backfill slurry with flow characteristics that met pipeline transportation requirements. The backfill body was formed by solidified backfill slurry and its mechanical strength met the requirements of backfill mining technology, with good heat storage/heat release characteristics that met the requirements of geothermal energy exploitation. The heat transfer behavior of the backfill body plays a key role in geothermal energy exploitation and stope cooling, and it is the core of the phase-change heat storage backfill system.
(3) We simulated the backfill body's heat storage process. The results showed that the higher temperature region gradually expanded from surrounding rocks to the interior of the backfill body, and the expansion was faster in the early stage and slower in the later stage. Moreover, the influence of the airflow on heat storage was obviously less than the surrounding rock. Therefore, the heat transfer of the surrounding rock was the main heat storage mode for the backfill body. However, the closer to the surrounding rock, the more obvious the temperature rise of the backfill body.
(4) We simulated the backfill body's heat release process. The results showed that the lower temperature region gradually expanded from the heat transfer tube to the periphery during the heat release process. The expansion was faster in the early stage and gradually slowed down in the late stage. The temperature of the backfill body around the middle tubes was affected by the cold fluid on both sides and there was an interaction region of heat release, so the lower temperature region around the middle tubes was larger. Moreover, the closer to the inlet, the more obvious the heat release interaction was between the tubes.