The thermal state and geothermal energy accumulation mechanism in the Xiong an New Area, China

Transcription

1 Research Article The thermal state and geothermal energy accumulation mechanism in the Xiong an New Area, China Energy Exploration & Exploitation 0(0) 1 14! The Author(s) 2019 DOI: / journals.sagepub.com/home/eea Yue Gaofan 1, Wang Guiling 1, Ma Feng 1, Zhang Wei 1 and Yang Zhijie 2 Abstract There are abundant geothermal resources in the Xiong an New Area, China. However, the thermal state and geothermal energy accumulation mechanism are not clear. Based on the geological conditions and the characteristics of the present geothermal field, a 2D model was established to analyze the process of mantle-derived heat conduction and to predict the distribution of the deep geothermal field. We calculated the terrestrial heat flow for the Rongcheng uplift and Niutuozhen uplift to be 64 and 75 mw/m 2, respectively. The geothermal resources in this area are controlled by a four-element model comprising heat conduction, structural uplift, large deep fractures, and convection within the reservoir. Keywords Xiong an New Area, thermal state, geothermal energy accumulation mechanism, numerical modeling, four-element model Introduction Geothermal resources are increasingly of interest to diverse communities due to their lowcost, sustainable and eco-friendly characteristics. The development and utilization of 1 Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences, Shijiazhuang, China 2 Key Laboratory of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun, China Corresponding author: Yue Gaofan, Institute of Hydrogeology and Environmental Geology, Chinese Academy of Geological Sciences, Shijiazhuang , China. gaofan3904@163.com Creative Commons CC BY: This article is distributed under the terms of the Creative Commons Attribution 4.0 License ( which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (

2 2 Energy Exploration & Exploitation 0(0) geothermal resources is of great significance for improving the energy structure and solving increasingly serious global environmental problems (Wang et al., 2017a). The Xiong an New Area has the most abundant geothermal resources in eastern China (Wang et al., 2017b; Zheng, 2015). In particular, in the Niutuozhen uplift and Rongcheng uplift, extremely abundant geothermal resources characterize the Xiongxian geothermal field and Rongcheng geothermal field. These resources have created great social and economic benefits with the development of geothermal energy exploitation (Zhang et al., 2013). Scholars have conducted many studies on the geology, geophysics, geochemistry, hydrology, drilling engineering, and other aspects of the geothermal resources in the Xiong an New Area (Chen, 1988; Li et al., 2014, 2017; Pang et al., 2011). The characteristics of the regional geology, deep fractures, and main reservoir have been clarified. On this basis, Pang et al. (2015) proposed two heat transfer mechanisms for the geothermal resources in the Xiong an New Area: heat conduction and heat convection. However, the existing research mainly focuses on shallow heat reservoirs, and the scope is usually limited to the periphery of the geothermal field. There are few studies on the deep thermal state and the heat accumulation mechanism at wider and deeper ranges. Numerical simulation methods can comprehensively study the influence of formation heat transfer, fluid heat transfer and radiative heat generation on the geothermal field and predict the distribution of the deep temperature field. At present, there are few numerical simulation studies on the Xiong an New Area; only Pang et al. (2011) established the horizontal fracture medium model of the Xiongxian geothermal field and predicted the recovery concentration and peak value of a tracer in well group tracer tests. The above studies focus only on the shallow reservoir system. The mechanisms of thermal conduction from the deep mantle source and the thermal conduction of deep faults remain unknown. Based on the geological conditions and the characteristics of the present geothermal field in the Xiong an New Area, we established a two-dimensional hydrothermal coupling model to analyze a four-element model comprising heat conduction, structural uplift, large deep fractures and convection within the reservoir, and we use this model to predict the distribution of the deep geothermal field. The results provide a reference for the development and utilization of the regional geothermal resources in the Xiong an New Area. Geological conditions Geothermal characteristics The study area is located in the north of the Jizhong depression (III level) in the Bohai North China basin (II level) on the Sino Korean paraplatform (I level), and in term of grade IV structural units, the study area includes the Rongcheng uplift, Langfang-gu an sag, Niutuozhen uplift and Baxian sag (Yu et al., 2017). The main thermal conductivity fractures in the region are Niudong and Rongcheng fractures (Pang et al., 2018; Wang et al., 2016; Figure 1). The depth of the Moho in this area is km and gradually becomes shallower from the northwest to the southeast. The depth of the Curie interface is km, with shallow burial in the middle and deeper burial at the edge. Preliminary data show that the terrestrial heat flow is between 79.1 and 90 mw/m 2, and the geothermal gradient is generally C/100 m, with a maximum value of 12.6 C/100 m (Chen et al., 1995).

3 Gaofan et al. 3 Figure 1. Structural map of the study area (from Zhang et al., 2017). The distribution of Cenozoic units in the uplifts and sags shows drape sedimentation. The quaternary and upper tertiary sandstone units are nearly horizontal, and the lower tertiary units have gentle dips. The lower strata mainly include Ordovician, Cambrian, Jixian, and Changcheng limestone and Archean metamorphic rocks (Li et al., 2014); Figure 2. Geothermal reservoir There are two types of reservoirs in this area: pore reservoirs and bedrock fissure karst reservoirs. The thermal pore reservoir is mainly composed of the Neogene Minghuazhen and Guantao formations. The Minghuazhen formation has a wide distribution range, while the Guantao formation is only distributed along the edge of the Niutuozhen uplift. The bedrock fissure karst reservoirs are the main reservoirs. In the Niutuozhen uplift, the Jixian Wumishan formation is the main reservoir, and in the Rongcheng uplift, the Changcheng Gaoyuzhuang formation is the main reservoir (Guo and Li, 2013). The variations in the thermal conductivity of the formations and large deep faults are the main factors controlling the geothermal resources in the study area. The surrounding thick rock cover with low thermal conductivity hinders heat conduction. The high thermal conductivity in the uplift region provides favorable conditions for heat flow and causes the heat

4 4 Energy Exploration & Exploitation 0(0) Figure 2. Geological profile of the study area (from Pang et al., 2018). to accumulate in the uplift region. The Niudong fault and the Rongcheng fault interact with the deep heat source and play a role in heat conduction. Tertiary and quaternary sediments are caprocks with low porosity, permeability, and thermal conductivity. These factors contribute to the rich geothermal resources in the Xiong an New Area. Model building Hydrothermal migration is a necessary condition for shallow geothermal resources. Based on geological conditions, this study mainly considers the mechanism of deep heat accumulation and the natural density differential convection process and establishes a 2D profile model of the Xiong an New Area (Figure 3). The top of the model is a constant temperature boundary, the bottom is a constant heat flux boundary and the side boundaries are open boundaries. Mathematical model The mathematical model of heat conduction is established as ðqcþ þ q wc w u rtþrq ¼ Q q ¼ k eff rt ðqc Þ eff ¼ ð1 hþq p C p þ hq w C w k eff ¼ ð1 hþk p þ hk w

5 Gaofan et al. 5 Figure 3. Conceptual model diagram. ðqcþ eff is the effective volume heat capacity, which is the heat capacity of rock and water. T is the temperature, K. C is the constant pressure specific heat capacity, J/kg. u is the water velocity, which comes from Darcy s law, m/s. k is the heat conductivity, W/mK. k eff is the effective thermal conductivity from the volume average method. Q is the source sink term. h is porosity. Subscripts p and w represent the rock skeleton and water, respectively. The Brinkman equation is used to describe the water convection l j u þrp rl ru þ ru h ð ÞT ¼ qgbðt T 0 Þ ru ¼ 0 l is the kinetic viscosity, Pas. j is the permeability, m 2. p is the hydrodynamic pressure, Pa. q is the water density, kg/m 3. g is the gravitational acceleration, 9.8 m/m 2. b is the water thermal expansion coefficient, with a value of in this case. T 0 is the reference temperature, K. Coupling software COMSOL is a large and advanced numerical simulation software, which is widely used in scientific research and engineering calculation in various fields, and can simulate various physical processes. It has efficient calculation performance and excellent direct coupling analysis ability of multiple fields, and realizes high-precision numerical simulation. In this paper, we defined a temperature-dependent force that causes fluid flow, which in turn causes convective heat transfer and affects temperature distribution. Two-way coupling of flow and heat transfer is achieved by COMSOL software. Critical parameters The data show that the thermal flow value at the bottom of the study area is 90 mw/m 2, which is stable and evenly distributed. The constant temperature is 14.8 C at the top. The thermal conductivity of the bedrock, the main heat reservoir, and the caprocks are 3.2, , and 1.4 W/mK, respectively (Wang et al., 2017b). The faults have a higher thermal conductivity of 4.5 W/mK. The porosity of the main reservoir is approximately

6 6 Energy Exploration & Exploitation 0(0) Table 1. Physical property parameters of the formations. Formation Thermal conductivity (W/mK) Density (kg/m 3 ) Specific heat capacity (J/kgK) Porosity (%) Permeability (m 2 ) Q þ N Ed Es Ek Qn Jx Ch Ar Bedrock %, and the permeability is approximately m 2. Detailed parameters are shown in Table 1. Initial condition and discretization In this paper, the process of heat transfer from the base to the shallow part is studied, focusing on deep heat conduction, intralayer convection, and the thermal control mechanisms of fractures and uplifts under natural conditions and ignoring the impacts of human activities and geological events such as volcanic eruptions. The initial temperature field model is set at 14.8 C, and the initial pressure is obtained by hydrostatic pressure balance. We used a triangulation grid to accurately describe the thermal conduction process in the fractures. The mesh was constructed at the fracture site, and a total of 5925 domain units and 979 boundary elements were obtained (Figure 4). The simulated total time was determined according to the results of the ground temperature field calibration (2 million years in this study), and the time step was 10,000 years. Results Model calibration Zhou et al. (1989) studied the distribution of shallow temperature in the Xiong an New Area and predicted the temperature of important sections based on borehole temperature measurement data. These results can be used as an important basis for our model identification. In general, the shallow temperatures of the Niutuozhen uplift and Rongcheng uplift are higher than those of the Gu an sag and Baxian sag. The temperatures of wells Rong-1, Xiong-102, Xiong-4, and Jia-3 at 1000 m depth are 58 C, 66 C, 57 C, and 43 C, respectively (Figure 5). In 2011, researchers found a temperature of 201 C at a depth of 6027 m in the Huabei oilfield in the Baxian sag (Zhao et al., 2011). These five main temperature measurement points can control the horizontal and vertical directions and can calibrate the model results well.

7 Gaofan et al. 7 Figure 4. Grid division display diagram. Figure 5. Distribution of temperature at 1000 m depth along section I-I 0. Figure 6 shows the temperature distribution at a depth of 1000 m obtained by simulation. The simulated temperature is consistent with the well logging results. The temperatures in the Niutuozhen and Rongcheng uplifts are relatively high, while those in the Gu an and Baxian sags are relatively low. Based on fitting the simulated temperatures to those observed at the five control points, the simulated and well logging temperatures at Rong-1 and Xiong- 4 are consistent, and the temperature differences between the simulated and well logging results at Xiong-102 and Jia-3 are less than 1 C. The simulated temperature in Niudong-1 at

8 8 Energy Exploration & Exploitation 0(0) Figure 6. Comparison of simulated temperature and log temperature (1000 m depth along section I-I 0 ). Figure 7. The distribution of temperature at 100,000 years. a depth of 6027 m is C, resulting in an error of only 2.1 C. Taking measurement error into account, our model fitting effect is good. Temperature distribution Figures 7, 8, and 9 show the temperature distributions at times of 100,000, 1 million, and 2 million years. In general, all the regional energy comes from basal heat conduction, and the temperature gradually increases over time. The influence of fractures is obvious. Due to the high thermal conductivity of fractures, the thermal conduction rate in the fracture channel is faster, and the heat reaches shallow formations more rapidly. Therefore, a geothermal contour bulge forms near the fracture (Figure 9). Figure 10 shows the temperature distribution at depths of 5 and 7 km and demonstrates that the temperature of the fracture is significantly higher than that of the surrounding strata. The temperatures at the Rongcheng fault at 5 and 7 km depth reach 170 C and 225 C, respectively, which are C higher than the surrounding strata.

9 Gaofan et al. 9 Figure 8. The distribution of temperature at 1,000,000 years. Figure 9. The distribution of temperature at 2,000,000 years (present day). Figure 10. The distribution of temperature at 2, 5, and 7 km depth along transect I-I 0. The temperature of the Niudong fault at a depth of 7 km is close to that of the strata in the Baxian sag. This pattern occurs because the Niudong fault is close to the basement of the Baxian sag at a depth of 7 km, and the thermal conductivity of the basement is equal to that of the fault.

10 10 Energy Exploration & Exploitation 0(0) Figure 11. Density flow field distribution at the present day. The temperature distribution at a depth of 2 km shows that the temperature in the shallow uplift region of the model is significantly higher than that in the sag region. The temperatures of the Rongcheng uplift and the Niutuozhen uplift can be up to 83 C, whereas the Gu an sag and Baxian sag have the lowest temperatures, including the minimum Gu an sag temperature of 58 C (Figure 10). This pattern occurs because the upper part of the sag features a very thick low thermal conductivity cap, with a slow heat conduction rate, which at the same time receives little heat energy. In contrast, the formations in the uplifts have high thermal conductivity and receive more heat energy. However, a comparison of temperatures at a depth of 5 km with those at a depth of 7 km reveals that the temperature in the deep uplift region is lower than that in the deep sag region, which is contrary to the pattern observed in the shallow region. The temperature in the deep part of the Baxian sag is 234 C (7 km), which is higher than those in the Rongcheng and Niutuozhen uplifts. This elevated temperature is because all formation heat comes from the bottom of the model, and the high thermal conductivity characteristic of the uplift regions is more conducive to the transport of heat to the shallow part. Under the input of set amount of heat, the temperature of the shallow part of the uplift region is relatively high, while that of the deep part is relatively low. Therefore, the shallow part of the uplift and the deep part of the sag are rich in geothermal resources and may represent target areas for future development. Natural convection in the reservoir The fluid density decreases with increasing temperature, and the density difference causes convection in the geothermal fluid. Based on the changes in the density of pure water, this study describes the flow driving force caused by the density difference using volume force (which is calculated with the following formula: qg T). The fluid streamline distribution in the heat reservoir is shown in Figure 11. The parts in red with larger widths represent high flow rates, while the blue thin lines represent low flow rates. As shown in the figure, under the influence of density differences, the deep hightemperature fluid moves towards shallow formations, but the flow velocity distribution is not uniform. Areas close to a fault, such as the Rongcheng fault and the Niudong fault, have high flow velocities, which can reach m/s, and there is no obvious convection in the reservoir far from the fault. Moreover, areas with high flow velocity have the same strong convection effect and even form vortices of different sizes. This phenomenon

11 Gaofan et al. 11 Figure 12. The near-surface terrestrial heat flow along transect I-I 0. occurs because the fracture acts as a heat upwelling channel, and the temperature differences and density differences near the fracture are large, resulting in a larger flow driving force and a higher flow velocity, thereby easily forming vortices. The vortices of different sizes contribute to making the temperature distribution in the reservoir consistent. Heat flow The near-surface geothermal heat flow value in the study area is calculated, and the results are shown in Figure 12. The heat flow in the uplifts is higher than that in the sag regions. The heat flow values in the Rongcheng uplift, Gu an sag, Niutuozhen uplift, and Baxian sag are 64, 32, 75, and 37 mw/m 2, respectively. According to the research results of Li et al. (2014), the heat flow value in the Niutuozhen uplift is mw/m 2, and that in the Rongcheng uplift is mw/m 2. The calculated heat flow value is slightly smaller than these results because the data are all from geothermal wells without accounting for the influence of convection factors, and the results from Li are relatively high. This model weakens the influence of convection to some extent. Although the results are slightly low, they still have reference significance. Well temperature prediction The Xiong an New Area is rich in geothermal resources. At present, the Rongcheng geothermal fields and Xiongxian geothermal fields have been identified. The thermal reservoirs are mainly in the Minghuazhen and Guantao formations and represent bedrock fissure karst reservoirs. These reservoirs primarily have development temperatures of less than 100 C and development depths of less than 2000 m. To further enhance the utilization of clean energy by exploiting geothermal energy, it is particularly important to develop deeper and higher temperature geothermal resources and establish efficient cascade development and utilization of geothermal energy. However, deep high-temperature geothermal drilling is difficult and costly, so accurate well location can play a positive role in saving on development and utilization costs in deep geothermal projects. Based on the temperature distribution from the model, four boreholes were prepared for the Rongcheng uplift (Z1), Gu an sag (Z2), Niutuozhen uplift (Z3), and Niudong fault (Z4).

12 12 Energy Exploration & Exploitation 0(0) Figure 13. Distribution of borehole temperature and geothermal gradient. Considering the current drilling efficiency limit, the temperature distribution is calculated only to a depth of 6 km. See Figure 5 for each drilling location. Figure 13 shows the temperatures and geothermal gradients of the four wells. There is a positive linear relationship between the well temperature and depth in the four boreholes. Z1 and Z3 have the highest temperature at depths less than 3.5 km, and Z2 has the lowest temperature because it contains the thickest sedimentary cover. Below a depth of 3.5 km, the temperature in well Z4 becomes the highest because well Z4 approaches the Niudong fault, and the high thermal conductivity of the fault transports considerable heat. From the perspective of the geothermal gradient distribution, Z1 and Z3 have the largest geothermal gradients at depths of less than 1 km, which is approximately 5 C/100 m. Below 1 km, the geothermal gradient stabilizes at C/100 m. Z2 and Z4 are affected by faults, and the variation in the geothermal gradient is complex. In particular, the geothermal gradient of Z2 increases to 5.8 C/100 m at a depth of km within the range of the influence of the Rongcheng fault. After passing through the fault, the geothermal gradient decreases. Below 5 km, the well approaches a fault again, and the geothermal gradient increases again. In summary, the most efficient method to obtain high temperatures is to drill the hanging wall of the fracture and to not cross the fracture. The ideal situation is to expose the fracture upon completion of the drilling to obtain the maximum bottom hole temperature. Conclusions We established a 2D profile model of hydrothermal transport in the Xiong an New Area and analyzed the distributions of the geothermal field, thermal flow value, well temperature, and geothermal gradient. The conclusions are as follows: The geothermal accumulation in the Xiong an New Area is a four-element model. Deep heat conduction is the source of heat. Deep faults provide an advantageous channel for the upwelling of heat, and the highest temperatures can be obtained near the faults.

13 Gaofan et al. 13 Geothermal resources are relatively abundant in the shallow parts of uplifts and the deep parts of sags. Intralayer convection occurs at different scales around the fault, which tends to make the temperature distribution in the heat reservoir consistent. The terrestrial heat flow values we obtained in the Xiong an New Area are 64 mw/m 2 for the Rongcheng uplift and 75 mw/m 2 for the Niutuozhen uplift. The best drilling location is located in the hanging wall of the Niudong fault, and the designed well should terminate at the depth of the fault. In this study, the temperature of the Z4 well at 6 km depth reaches 200 C. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was funded by National Key R&D Program of China (No. 2018YFC ), Special Funds for Basic Scientific Research of the Chinese Academy of Geological Sciences (SK201704, SK201501, YYWF ), National Natural Science Foundation of China (NO ), and a project of the China Geological Survey (DD ). ORCID id Yue Gaofan Zhang Wei References Chen MX (1988) Geothermal Resources in North China. Beijing: Science Press. (in Chinese) Chen MX, Wang JY and Deng X (1995) Advances in geothermal in China. Earth Science-Journal of China University Geoscience 20(4): (in Chinese) Guo SY and Li XJ (2013) Reservoir stratum characteristics and geothermal resources potential of Rongcheng uplift geothermal field in Baoding, Hebei. Chinese Journal of Geology 48(03): (in Chinese) Li H, Yu JB, Lv H, et al. (2017) Gravity and aeromagnetic responses and heat-controlling structures of Xiongxian geothermal area. Geophysical and Geochemical Exploration 41(02): (in Chinese) Li WW, Rao S, Tang XY, et al. (2014) Borehloe temperature logging and temperature field in the Xiongxian geothermal field, Hebei Province. Chinese Journal of Geology 49(3): (in Chinese) Pang JM, Pang ZH, Lv M, et al. (2018) Geochemical and isotopic characteristics of fluids in the Niutuozhen geothermal field, North China. Environmental Earth Sciences 77(1): Pang JM, Wang SF, Sun CX, et al. (2011) Interpretation and analysis of tracer test in Xiongxian geothermal field. Urban Geology 6(2): (in Chinese) Pang ZH, Pang JM, Kong YL, et al. (2015) Large karstic geothermal reservoirs in sedimentary basins in China: Genesis. Energy potential and optimal exploitation, In: Proceedings world geothermal congress 2015, Melbourne, Australia, April Wang GL, Zhang W, Liang JY, et al. (2017) Evaluation of geothermal resources potential in China. Acta Geoscientica Sinica 38(04): (in Chinese)

14 14 Energy Exploration & Exploitation 0(0) Wang GL, Zhang W, Lin WJ, et al. (2017) Research on formation mode and development potential of geothermal resources in Beijing Tianjin Hebei region. Geology in China 44(6): (in Chinese) Wang YB, Ding WP, Tian Y, et al. (2016) Genetic analysis on high-temperature geothermal water in Niutuo Geothermal Field, Heibei Province. Urban Geology 11(3): (in Chinese) Yu CC, Qiao RX and Zhang DS (2017) The basement tectonic characteristics from interpretation of aeromagnetic data in Xiong an region. Geophysical and Geochemical Exploration 41(3): (in Chinese) Zhang DZ, Liu ZG and Lu HL (2013) Geothermal Resources in Hebei Province. Beijing: Geological Science Press. (in Chinese) Zhang YM, Chang J and Liu N (2017) Present-day temperature pressure field and its implications for the geothermal resources development in the Baxian area, Jizhong Depression of the Bohai Bay Basin. Natural Gas Industry 37(10): (in Chinese) Zhao XZ, Jin FM, Wang Q, et al. (2011) Niudong 1 ultra-deep and ultra-high temperature subtle buried hill field in Bohai Bay Basin: Discovery and significance. Acta Petrolei Sinica 32(06): (in Chinese). Zheng LY (2015) The characteristic and evaluation of geothermal resources in Beijing-Tianjin-Hebei Region. China University of Geosciences, Beijing. (in Chinese) Zhou RL, Liu QS, Zhang J, et al. (1989) The geological features and exploitive prospects of the geothermal field of salient type of bed rock of Niutuozhen in the fault basin of north China. Bull 562 Compr Geol Brig Chin Acad Geol Sci 7: (in Chinese)