Herong Zheng
Jun Luo
Jianyun Feng
Yan Zeng
Mingchuan Wang- Sinopec Petroleum Exploration and Production Research Institute, Beijing, China
The southeast coastal areas in China have distributed lots of granite outcrops of different periods. Previous research has shown that granite geothermal reservoirs are also distributed under sedimentary basins in these areas, such as in Zhangzhou basin. Therefore, granites with fractures buried in deep can be used as a potential deep geothermal reservoir in these areas. In order to study geological conditions of the deep granite reservoir and discuss the genesis of the deep granite geothermal system, rock geochemistry and zircon U-Pb chronology from outcrop and parts of the drilling cores of granitic rocks have been analyzed, combined with the anatomy of the deep seismic data and electromagnetic detection data in selected area. Based on the results of geochemistry and zircon U-Pb chronology, most granites in this area are of Yanshanian periods. Based on the seismic data, the thickness of the overlying strata on granite in Huangshadong area of Huizhou City is up to 1.5 km. According to the regional geological survey, multi-stage joints are developed in the granite, and most of hot springs rise from intersection of fracture with different directions to the surface. The heat source in the study area mainly comes from the mantle carried up by the deep NNE-trending faults. There are a large number of thermal springs at the intersection of the surface and the NW-trending fault, and the NW-trending fault provides the drainage conditions for the upwelling of underground thermal springs. There is a huge amount of deep granite geothermal resources in the southeast coastal area. The analysis of deep granite geological conditions and genetic models can provide guidance for the evaluation of deep granite geothermal resources and the further optimization of favorable zones in these areas.
Introduction
Geothermal resources, as a renewable and clean energy, have been widely used globally in recent years, and the utilization has been increasing year by year (Bertani, 2012; Bertani, 2016; Lund and Boyd, 2016). Recently, most of the geothermal resources under exploitation and utilization are hydrothermal geothermal resource, which is widely distributed in southeastern areas of China (Li Dewei and Wang Yanxin, 2015). Deep geothermal resources are considered as an important part of future energy supply due to their huge thermal energy storage and reserves. Most of the deep geothermal reservoir are granite (Brown et al., 2012), and Metamorphic rocks (Li et al., 2019) beneath the sedimentary basins and numerous studies have shown that the granite radioactive heat production has significant contribution to the heat source, such as Australia Cooper basin geothermal field caused by granites with intrusion age less than 0.5 Ma (Goldstein et al., 2008), Rose-manowes geothermal field in England due to Early Permian granite heat generation (Richards et al., 1992), the Soultz geothermal field in France due to Late Paleozoic granite heat generation (Genter et al., 1995).
In the southeast coastal area of China, ground or near surface geothermal manifestations are widely distributed, such as hot springs and other surface heat. They are mainly located in the fractures within igneous rock, such granites. Granites in South China are widely distributed, with an outcropping area of 20,000 km2, accounting for about 1/5 of the area, and it is considered to be formed in three phases, i.e., Caledonian, Indosinian and Yanshanian period (Zhou et al., 2006) (Figure 1). Radioactive elements from those granites such as U, Th, K, are important radioactive elements to generate heat by atom decaying. Meanwhile, Southeast coastal areas in China is the second largest region with high heat flow value (Hu et al., 2000). Xi et al. (2018) analyzed the gravity anomaly in Guangdong province and concluded that the decay of thermal elements in granite may be an important part of the geothermal heat sources in the southeast coastal areas. The region is rich in geothermal resources, covering medium-low temperature, medium-high temperature and high-temperature geothermal resources.
FIGURE 1. The granite distribution map of the study area (Modified after Zhang et al., 2020).
In order to further study geological conditions of the deep granite reservoir and discuss the genesis model of the deep granite geothermal system, this paper adopted rock samples from outcrops and parts of the drilling core rocks to analyze the formation background of granite and discuss its geothermal significance to the southeast region in China.
Geological Setting
Geologically, the southeast coastal region of China is part of the Cathaysia block (Wang et al., 2013; Deng et al., 2019). These areas have experienced multiple periods of intense magmatic events since the Paleozoic, mainly manifested as a large number of Paleozoic and Mesozoic granitic rocks (Wang et al., 2013) and Cenozoic mafic magmatism (Gong and John, 2014). Large-scale structural deformation caused by emplacements are very strong to the stratum reconstruction, forming a series of faults with different scales and directions.
Granite reservoirs are distributed among the southeastern coastal areas, such as in Zhangzhou, Fuzhou, Fengshun, Yangjiang, and southern Hainan Province. In these areas, the distribution of reservoir is mainly controlled by NE-trending faults. Most of these granites were formed during the Yanshanian period. There are two types of heat source: Crust source originating from radiogenic heat from K-, Th- and U-bearing minerals of granites and mantle source connecting to the nearby ground by the deep faults. In the Quaternary, most of the radiogenic heat of granite was lost along the fault zone and could not constitute a source of geothermal heat in the fault. Therefore, it is easy to form deep fracture zone and become the channel for hot water to migrate from the deep mantle. In southeast coastal area, 74% of hot springs are located in fault zone or lithological contact zone in magmatic plutons.
The distribution of most granite geothermal fields are mainly located in basins, such as Fuzhou basin, Zhangzhou basin, Sanshui basin (Figure 1). They are mainly controlled by the deep and large fault in NE direction, and most of them are along the Neo-Cathaysian faults. The zone has experienced many strong tectonic movements and multiple periods of magma intrusion, resulting in secondary faults, rock mass fragmentation and joint fractures near the fault zone, providing space and channels for the storage and migration of geothermal fluids. At the same time, the tectonic activities of deep and large faults not only promote the formation of heat storage space, but also communicate the spatial connection between deep geothermal fluids and shallow geothermal reservoirs, becoming an important heat transfer channel in geothermal fields. From the regional analysis, most of the geothermal fields are linear distribution along the fault zone, mainly exposed on the deep fault axis. The rest are mostly distributed among the deep and large faults, which are locally influenced by the secondary NW tensioned water-conducting faults or pinnate faults. The distribution of heat flow in South China is shown in Figure 2. It contains three obvious thermal structure divisions: Eastern, central and southwestern regions. The eastern and southwestern regions are characterized by high heat flow, while the heat flow in the central region is relatively lower. The heat flow in the eastern and southwestern regions are generally higher than 70 mW/m2, while that in the central region is lower than 60 mW/m2. In the eastern region, there are 2 NE-trending high value abnormal belts along the Fuzhou-Zhangzhou zone. These two anomaly zones are characterized by high heat flow values (>80 mW/m2, locally up to 220 mW/m2). The average heat flow in northern Jiangsu Basin is 71.2 mW/m2, while that in southern Jiangsu and southern Anhui areas is 62.7 mW/m2. In the central region, from Baise to Dabie Mountains, lower heat flow is distributed along the NE-direction. In the southwestern region, the higher heat flow is concentrated only in southwestern Yunnan. The average of heat flow in the north of Sichuan Basin is 53.6 mW/m2, while that in southern Yunnan is 76.7 mW/m2. In the middle of northern Jianghan Basin, the heat flow is lower in the north and higher in the south. The average heat flow in the north of the fault is 47.8 mW/m2, and that in the south is 61.6 mW/m2. The average heat flow value of 64.2 mW/m2 in southern China is higher than that of 61 mW/m2 in other mainland parts in China (Wang et al., 1990), which is close with that of the globe for 65 mW/m2 (Yuan, 2006). The abnormally high value areas in South China are mainly distributed in the suture zone or regional deep fault zone. Some granite geothermal reservoirs in typical areas in South China have been studied by previous researchers, including Fuzhou, Zhangzhou, Huizhou, Yangjiang areas (Zhu 2015; Yang 2016; Zhang et al., 2021; Zhou, 2015). The Fujian granite geothermal reservoir is located around Fuzhou basin (Figure 3A). Most of them are Yanshanian granites. Early detailed geothermal resource exploration shows that the geothermal reservoir is an area of 9 km2 with a north-south length of about 5 km and an east-west width of about 2 km. The basement consists of Yanshanian granodiorite, intrusions of medium-coarse biotite granite and fine-grained granite, as well as a large number of intermediate-acid and intermediate-basic dike rocks (Zhu 2015). Zhangzhou granite geothermal reservoir is located around Zhangzhou Basin (Figure 3B). The Quaternary sediments on granite are a set of Marine and continental deposits with a thickness of 20 m–30 m. There are several granite rock masses, such as Zhangzhou granite, Chengxi granite and Maping granite, which are around Zhangzhou geothermal field (Yang, 2016). In Huizhou area, Huangshadong geothermal field is located in the northeast of Huizhou city (Figure 3C). This area has experienced several periods of intense magmatic events since the Paleozoic. The granites include Caledonian granodiorite, Indo-Chinese granites, Yanshanian granites and Neogene tholeiitic basaltic dikes (Zhang et al., 2021). Xinzhou complex granite is located in the coastal area of west Guangdong, Yunkai block geologically, on the southwest Cathaysia block, with a total outcropping area of 350 km2. The rock is composed of medium-coarse-grained macro porphyritic biotite monzonite and fine-grained porphyritic biotite quartz syenite, surrounded by the Proterozoic Yunkai Group plagioclase gneiss. The latest outcropped strata around this area are Late Cretaceous sandstone and conglomerate, distributed with NE direction. Xinzhou complex granite mainly includes Xinzhou granite, Dongping granite and Naqin rock granite (Zhou, 2015).
FIGURE 2. The Heat flow map distribution of the study area (Modified after Zhang et al., 2020).
FIGURE 3. Simplified geological map of selected regions in Southeast Coastal Area: (A) Fuzhou area; (B) Zhangzhou area; (C) Huizhou area (Modified after Zhu 2015; Yang 2016; Zhang et al., 2021).
Favorable granite geothermal reservoir is distributed in the exposed magmatic rock mass, near the late intrusive dike or the contact zone between rock mass and surrounding rock mass. The distribution of geothermal field is not only related to the fracture and fracture development of primary rock mass and contact zone rock due to the multi-stage intrusion of magma or the influence of deep fault structure, but also related to the geological age of intrusive rock.
Sample Description and Analitical Methods
In order to interpret magmatism in the study area of granite, a total of 20 samples of magmatic rocks were collected in the selected areas, such as Fuzhou, Zhangzhou and Huizhou area (Figure 3).and they were measured by principal trace element analysis, and zircon U-Pb dating analysis was conducted for six granite samples.
The granite samples collected were mainly distributed along a line of about 1,300 km with an exposed area of about 10–784 square kilometers. These granites are mainly granitic (Figure 4), with a few being orthogonality, quartz dimorphism and granitic diorite, with granular structure of different thickness and fineness (Figure 5). The mineralogy of these samples is relatively simple, consisting mainly of plagioclase (30–42%), potash feldspar (30–8%), quartz (20–30%), biotite (5%), and a small amount of amphibole (2–10%). Accessory minerals include sphere, zircon, and titanium-iron oxides.
FIGURE 4. Photograph showing outcrops of granites in different areas, (A) granite outcrop with fractures in Minhou, Fuzhou area; (B) granite outcrops in Qinqian, Fuzhou area; (C) granites with fractures in Huizhou; (D) granites in Zhangzhou.
FIGURE 5. Microphotographs showing mineral assemblage of representative granites. The abbreviations are: Pl-plagioclase, Qz-quartz, Bt-biotite, Kfs-K-feldspar, Amp-amphibole.
U-Pb dating and trace element analyses of zircon were conducted synchronously by LA-ICP-MS at the State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, Wuhan. Detailed operating conditions for the laser ablation system and the ICP-MS instrument and data reduction are the same as description by Liu et al. (2010). Laser sampling was performed using a GeoLas 2005. An Agilent 7500a ICP-MS instrument was used to acquire ion-signal intensities. A “wire” signal smoothing device is included in this laser ablation system, by which smooth signals are produced even at very low laser repetition rates down to 1 Hz (Hu et al., 2012b). Helium was applied as a carrier gas. Argon was used as the make-up gas and mixed with the carrier gas via a T-connector before entering the ICP. Nitrogen was added into the central gas flow (Ar + He) of the Ar plasma to decrease the detection limit and improve precision (Hu et al., 2008a; Liu et al., 2010). Each analysis incorporated a background acquisition of approximately 20–30 s (gas blank) followed by 50 s of data acquisition from the sample. The Agilent Chemstation was utilized for the acquisition of each individual analysis. Off-line selection and integration of background and analyze signals, and time-drift correction and quantitative calibration for trace element analyses and U-Pb dating were performed by ICPMS DataCal (Liu et al., 2008).
Zircon 91,500 was used as external standard for U-Pb dating, and was analyzed twice every five analyses. Time-dependent drifts of U-Th-Pb isotopic ratios were corrected using a linear interpolation (with time) for every five analyses according to the variations of 91,500 (i.e., two zircon 91,500 + 5 samples +2 zircon 91,500) (Liu et al., 2010). Preferred U-Th-Pb isotopic ratios used for 91,500 are from Wiedenbeck et al. (1995). Uncertainty of preferred values for the external standard 91,500 was propagated to the ultimate results of the samples. Concordia diagrams and weighted mean calculations were made using Isoplot/Exver3 (Ludwig, 2003). Trace element compositions of zircons were calibrated against multiple-reference materials (BCR-2G and BIR-1G) combined with internal standardization (Liu et al., 2010). The preferred values of element concentrations for the USGS reference glasses are from the GeoReM database.
Results
U-Pb Dating
The concordance ages and analysis results of typical zircons of six samples are shown in Figure 5. They are collected from different areas, such as Fuzhou area (GA02, QQ02) (Figure 3A), Zhangzhou area (ZZ02, ZZ07, ZZ11) (Figure 3B) and Huizhou area (HZ05) (Figure 3C). It can be seen from the CL image that the zircons generally have good autogenesis (Figure 6). Most of them are columnar or long-column-cone shaped with a length of 50–360 m. The aspect ratio of most of them is about 1.5:1–3:1, which is similar to the characteristics of magmatic zircons (Wu Yuanbao and Zheng Yongfei, 2004). The dating results of six samples show that the rock mass in the study area has a wide age distribution, ranging from 94.7 to 138.4 Ma. The LA-ICP-MS U-Pb zircon analysis spots are presented in Table 1, and the results are illustrated in Figure 6.
FIGURE 6. The 206Pb/238U-207Pb/235U concordia plot and the weighted mean age of representative zircons of the granites.
TABLE 1. Zircon U-Pb data of granites in the study area.
A total of 26 valid data were obtained from GA02. These zircon grains have moderate and variable concentrations of Th (102–200 ppm) and U (125–318 ppm), corresponding Th/U ratios of 0.8–1.14, showing the characteristic of magmatic zircon (Hoskin and Black, 2000; Sun and McDonough, 1989). These data yield a weighted mean 206Pb/238U age of 117.95 ± 0.47 Ma (MSWD = 0.45), representing the formation age of the granite. A total of 28 valid data were obtained from QQ02. These zircon grains have moderate and variable concentrations of Th (127–291 ppm) and U (287–606 ppm), corresponding Th/U ratios of 0.37–0.85, showing the characteristic of magmatic zircon. These data yield a weighted mean 206Pb/238U age of 112.92 ± 0.85 Ma (MSWD = 1.6), representing the formation age of the granite. A total of 25 valid data were obtained from ZZ02. These zircon grains have moderate and variable concentrations of Th (108–535 ppm) and U (154–1,090 ppm), corresponding Th/U ratios of 0.3–0.99, showing the characteristic of magmatic zircon. These data yield a weighted mean 206Pb/238U age of 121.3 ± 1.2 Ma (MSWD = 3.2), representing the formation age of the granite. A total of 29 valid data were obtained from ZZ07. These zircon grains have moderate and variable concentrations of Th (50.5–318 ppm) and U (149–865 ppm), corresponding Th/U ratios of 0.34–0.72, showing the characteristic of magmatic zircon. These data yield a weighted mean 206Pb/238U age of 138.4 ± 1.1 Ma (MSWD = 1.6), representing the formation age of the granite. A total of 25 valid data were obtained from ZZ11. These zircon grains have moderate and variable concentrations of Th (58.7–837 ppm) and U (61–1,124 ppm), corresponding Th/U ratios of 0.62–1.53, showing the characteristic of magmatic zircon. These data yield a weighted mean 206Pb/238U age of 94.7 ± 0.25 Ma (MSWD = 1.01), representing the formation age of the granite. A total of 30 valid data were obtained from HZ05. These zircon grains have moderate and variable concentrations of Th (56.2–448 ppm) and U (88.3–548 ppm), corresponding Th/U ratios of 0.33–0.99, showing the characteristic of magmatic zircon. These data yield a weighted mean 206Pb/238U age of 135.1 ± 1.3 Ma (MSWD = 2.4), representing the formation age of the granite.
There is no obvious pattern of age distribution of such samples in the region. In general, the age of all samples in the study area shows that they are formed in Cretaceous age, indicating that Yanshanian tectonic thermal events are the dominant thermal events in the area. According to the U-Pb dating results, at least three acidic magmatic events occurred during Cretaceous: early stage of Early Cretaceous (about 138–135.1 Ma), late stage of Early Cretaceous (110–121 Ma) and the early stage of Late Cretaceous (about 93–96 Ma).
Geochemistry
20 fresh granite samples are obtained for major and trace element analyses. The results are listed in Table 2. The samples have very low loss on ignition (∼0.85 wt%), indicating little or no modification by weathering. Rock geochemical test results show that the rock mass geochemical composition analyses are as follows: SiO2 content ranged from 65.86 to 83.83% with an average of 73.64%; TiO2 content ranged from 0.08 to 0.55% with an average of 0.23%; Al2O3 content ranged from 7.38 to 15.8% with an average of 13.53%; MgO content ranged from 0.12 to 1.73% with an average of 0.48%; and (K2O + Na2O) content ranged from 4.76 to 8.87% with an average of 7.66%, belonging to high-K calc-alkaline series. The granites are characterized by weakly peraluminous (A/CNK = 1.03). At the diagram of SiO2 (Na2O + K2O) (Figure 7), sample dots fell within the granite region. The aluminum saturation index (A/CNK) ranges from 0.88 to 1.33 and ranges mainly from 1.0 to 1.1 on the diagram of A/NK–A/CNK (Figure 7).
TABLE 2. Analytical data of major (%) and trace element(ppm) of granites in the study area.
FIGURE 7. TAS (A) and A/CNK (B) diagrams of magmatic rocks.
On the primitive mantle normalized trace element spider diagram (Figure 8), Large ionic Lithophile Elements (LILE) such as K and Rb, Th, and U are enriched, and there is a strong loss of Ba, Sr, and High Field Strong Elements (HFSE) such as Nb, Ti, Ta, and P, which are similar to the geochemical characteristics of rocks in the subduction zone (Kelemen et al., 2003). On the Chondrite Normalized REE distribution diagram (Figure 8), the pattern is smooth and show light rare Earth element (LREE) concentration and heavy rare Earth element (HREE) losses. Eu in the pattern shows negative anomaly, indicating that the magma evolution of plagioclase crystallization separation. The value of ΣREE is between 50.50 × 10−6–401.02 × 10−6, with a mean of 183.31 × 10−6. The value of LREE/HREE is 0.84–12.46, with an average of 6.44, showing that obvious fractionation between light and heavy rare Earth. The value of (La/Yb)N is 0.67–19.24, while the average value is 7.18. The value of Eu is between 0.03 and 0.71, and the average value is 0.32. In primitive mantle normalized trace element spider diagram, most of the sample elements distribution features are similar, indicating that most of the samples have the same or similar formation patterns.
FIGURE 8. Chondrite Normalized REE distribution patterns (A), (C), (E) and primitive mantle normalized trace element spider diagram (B), (D), (F) of granitic rocks.
Discussion
Genetic Mechanism of Yanshannian Granite
Large-scale Yanshanian granites constitute the main part of southeastern China, and were formed related to the subduction of the Pacific plate, large-scale basaltic magma invasion, crust-mantle mixing, and extensive melting of crustal rocks (Liu et al., 2016; Pan et al., 2017). Meanwhile, other researchers have proposed the relationship between the lithospheric thinning in eastern China and the subduction of the Pacific plate (Wang and Shen, 2003; Wu et al., 2007). Li et al. (2017) discussed the starting time of palaeo-Pacific subduction. Through the study of Mesozoic tectonic transition and transformation in South China, they found that the palaeo-Pacific plate subduction had already started in the Late Triassic. It is considered that the crustal thickness of the South China Block gradually thinning from west to east is caused by two crustal mantle interactions, which resulted in two lithospheric delamination and subsidence, and both of them are related to the subduction of the palaeo-Pacific Plate, which were caused by plate tearing, subduction retraction and high Angle subduction respectively.
According to the results of this paper, at least three acidic magmatic events occurred during Cretaceous: early stage of Early Cretaceous (about 138–135.1 Ma), late stage of Early Cretaceous (about 112–121 Ma) and the early stage of Late Cretaceous (about 93–96 Ma).
The age of the granite sample (ZZ07) in Zhangzhou area and the granite sample (HZ05) in the south of Huangshadong geothermal field in Huizhou are 138.4 ± 1.1 Ma and 135.1 ± 1.3 Ma, indicating that the granite body is the product of magmatic activity in the early stage of early Cretaceous. At the diagram of SiO2/(Na2O + K2O) (Figure 7A), sample dots fell within the granite region. The granites are characterized by weakly peraluminous (A/CNK = 1.02–1.08) (Figure 7B). The aluminum saturation index (A/CNK) ranges from 1.14 to 1.33. On the SiO2 − FeOT/(FeOT + MgO) and SiO2 − Al2O3 diagram (Figures 9A,B), most samples are plotted in the POG (post-Orogenic Grabite) zone (Frost et al., 2001). On the Nb-Y and (Y + Nb)–Rb diagram (Figures 9C,D), samples are plotted within the scope of the volcanic island arc granite (VAG), within plate granite (WPG) and sys-collision granite (syn⁃COLG). Combined with the trace element characteristics of the samples, it can be inferred that the samples have dual characteristics of volcanic arc and within plate environment. The ages of the granite rock (GA02, QQ02) and the granite rock (ZZ02) around the Fuzhou geothermal field and Zhangzhou geothermal field are 117.9 ± 0.47 Ma, 112.9 ± 0.85 Ma and 121.3 ± 1.2 Ma, respectively, which are the products of magmatic activity in the late Early Cretaceous. It is the product of the same magmatism as most of the granites along the Fujian coast area in Fujian Province. At the diagram of SiO2/(Na2O + K2O), sample dots fell within the granodiorite region (Figure 7A). They are characterized by Metaluminous - peraluminous (A/CNK = 0.95–1.05) (Figure 7B). The aluminum saturation index (A/CNK) ranges from 1.42 to 1.68. On the SiO2 – FeOT/(FeOT + MgO) and SiO2 – Al2O3 diagram (Figures 9A,B), most samples are plotted out of the POG (post-Orogenic Grabite) zone (Frost et al., 2001). On the Nb-Y and (Y + Nb) – Rb diagram (Figures 9C,D), samples are plotted within the scope of the volcanic island arc granite (VAG), indicating rocks in this periods are mainly influenced by subduction related environment. The age of the sample granite rock mass (Z11) in coastal Zhangzhou is 94.7 ± 0.65 Ma, which is similar to the age of the rock mass (93 Ma, 95 Ma) in the south of Fuzhou geothermal field. It is the product of magmatic activity in the early stage of Late Cretaceous. At the diagram of SiO2/(Na2O + K2O), sample dots fell within the granite region (Figure 7A). They are characterized by also weakly peraluminous (A/CNK = 1.03–1.34) (Figure 7B). The aluminum saturation index (A/CNK) ranges from 1.06 to 1.38. On the SiO2 – FeOT/(FeOT + MgO) and SiO2 – Al2O3 diagram (Figures 9A,B), most samples are plotted in the POG (post-Orogenic Grabite) zone (Frost et al., 2001). On the Nb-Y and (Y + Nb)−Rb diagram (Figures 9C,D), samples are plotted within the scope of the volcanic island arc granite (VAG), within plate granite (WPG) and sys-collision granite (syn⁃COLG). Combined with the trace element characteristics of the samples, it can be inferred that the samples of this period is similar to the granite formed in ∼135 Ma. Li Jianhua (2013) divided the acidic magmatic rocks of the Cretaceous into three periods: Adakite and gneis-like granite at 145–137 Ma, calc-alkaline granite at 136-118 Ma and related volcanic rocks, bimodal volcanic rocks at 107–90 Ma and A-type granite at 145–137 Ma. This is similar to what we found in this study.
FIGURE 9. SiO2-FeOT/(FeOT + MgO), (A) SiO2 − Al2O3, (B) Y − Nb, (C) and (Yb +Nb) − Rb, (D) diagrams of granites.
Discussion on the Burial Depth of Deep Granite
Seismic Profiles in Huangshadong Area of Guangdong Province
The two seismic profiles in Huangshadong area of Guangdong Province (Figures 10A,B) (Kuang et al., 2020) reflect the sedimentary strata overlying granites. The sedimentary strata in the study area have undergone complex tectonic evolution, repeated tectonic uplifting, compression folding and fault processes have caused the sedimentary strata to have dramatic occurrence changes in a small range, and multiple periods of magmatic intrusion and eruption have complicated the contact relationship strata and rock mass. Therefore, the spatial continuity of stratigraphic interface and the stability of fault characteristics in seismic profile are poor in different ages.
FIGURE 10. Profiles of seismic inversion, (A), (B) and its geological interpretation in Huangshadong area (Kuang et al., 2020).
According to the profiles of seismic inversion (Figure 10), the buried depth of Tg (interface between sedimentary strata and granites) is generally located between 1,000–1,900 ms (corresponding depth of 1,400–2,200 m), and the seismic phase characteristics of Tg interface are significantly different from each other. The granite core sample from Well Huire one is similar to Yanshanian granite, indicating that seismic facies unit may correspond to the intrusive body in Yanshanian period, emplaced from the deep crust upward into shallow rock bed. 2050–2,900 ms (corresponding to the depth of 3,600–4,800 m) for the bedrock under the bottom of the interface (Td), reflection of geological units in phase axis continuity better, different from the adjacent granite rock mass. It may be older sedimentary rocks or metamorphic rocks with certain stratification.
Apparent Resistivity Inversion Profiles in North Hainan Province
In the area outside the east of Fushan Depression of Hainan Province, wide field electromagnetic detection has been conducted (Figure 11) (Tan et al., 2020). From the resistivity profile of the final inversion processing, a set of high resistance developed in the west of the L1 line and a set of low resistance developed in the east (Figure 11A). The two parts were clearly bounded and featured with great difference. Based on regional geological analysis, the corresponding geological model is established (Figure 11B), and it is preliminarily considered that the high resistivity in the east is a set of intrusive concealed granite, and the low resistivity in the west is a paleoproximal sedimentary stratum, and the middle fault controls its scale. According to the analysis of geological model, the basement of North Hainan area is composed of pre-Mesozoic sedimentary rocks and granite intrusive bodies, and the rock strains develop. The fault-block-graben-barrier is developed in this area, and the structure is complex. The fault-block-graben-barrier structure is developed in the east, which is similar to Fushan Depression. Moreover, the distribution of Paleogene strata is obviously controlled by the size of the fault depression. The L4 line extends southward to the granite outflowing area. According to the resistivity result diagram (Figure 11C), the electrical characteristics of the outflowing granite are similar to those of the latent granite (Figure 11D). By comparing the two survey lines, the latent high-resistivity stability exists. In the intersection area where survey line L1 and survey line L4 meet, the granite has a large scale, a buried depth of about 3000 m, and an overburden, which is a potential favorable area for deep high-temperature geothermal exploration.
FIGURE 11. Profiles of apparent resistivity inversion based on wide field electromagnetic method and corresponding geological interpretation in North Hainan Province (Tan et al., 2020).
The Significance of Granite to the Geothermal Energy
The heat source of geothermal energy mainly include heat transfer in the mantle, local anomalous heat transfers in the crust (magma chamber), and radioactive heat generation from granitic rocks. Surface heat flow, as a direct reflection of deep heat, usually accounts for 40–60% of the heat flow from the mantle. The heat flow caused by radioactive heat generation in the crust is slightly less than that from the mantle, accounting for about 40% and up to 50% in some areas (Wang Shejiao et al., 1999). The terrestrial heat flow value in the study area is about 70 mW/m2, which is higher than the average terrestrial heat flow value of 64.2 mW/m2 in South China (Yuan et al., 2006). The Yanshannian granite in the Southeast Coastal region has a high thermal conductivity, with an average thermal conductivity of 3.15 W/mK. This means that the granite in the deep part with a thickness of up to 3.5 km can better conduct the heat in the deeper part to the near surface. Consequently, the mantle heat gets better performance on the surface, which promotes the formation of the geothermal energy in the area.
Former research suggested that the heat generated by radioactive element decay in rocks is one of the main sources of Earth heat (Huang Shaopeng, 1992), and it is also a major factor controlling the temperature field distribution in the lithosphere (Sclater et al., 1980). Previous scholars proposed that the lower crust has little influence on the surface heat flow value (Hasterok and Chapman, 2007), and the research on the lower crust inclusion also shows that the radioactive heat generation rate of rocks in the lower crust is very low (Rudnick and Fountain, 1995). Thus, the heat source in the area mainly comes from the decay of radioactive elements in the middle and upper crust, and depends on the abundance of thermogenic elements, such as 238U, 232Th, and 40K (Rybach and Buntebarth, 1984).
Previous studies have analyzed the overall heat generation rate of radioactive elements in the study area shown as low heat generation rate of strata and high heat generation rate of granitic rocks (Kuang et al., 2020). The value of heat generation rate from Yanshanian rocks is far higher than the average global continental upper crust rock, with the average value of 5.50 μW/m3 and the maximum value up to 9.15 μW/m3. It is also much higher than the average measured value of the radioactive heat generation rate of granitic rocks in other regions, which was about 2.5 μW/m3 (Rybach, 1988), 3.3–3.8 μW/m3 (Wollenberg and Smith, 1987) and 3.0 μW/m3 (Rybach and Buntebarth, 1984).
This is mainly attributed to Yanshanian granites formed from remelting of pre-Cambian rocks in this region (Xu et al., 2007), The rocks have been experienced the redistribution process, and enriched large ion lithophile elements (LILE) of rocks, lead to extreme enrichment of LILE (including Th, U, K element), and then forming high value of heat generation rate in granites. Consequently, the high thermal conductivity and high radioactive thermogenic elements of Yanshanian granites are the reasons for the formation of high-temperature geothermal in this area.
The Significance of Granite to Deep Geothermal in Southeast China
Large-scale granite events occurred in the coastal area of southeast China during Mesozoic, especially Yanshannian (Li, et al., 2014; Gao et al., 2014; Deng et al., 2016). The Yanshannian granitoids magmatism in the southeast coastal area has a general trend from west to east, and the emplacement age of the rocks from old to young. Due to the high radioactive heat generation rate of granitoids, granitoids provide a direct heat source for crustal heat generation. Therefore, the widespread occurrence of granitoids magmatic activity in Mesozoic provides a significant condition for the formation of geothermal resources. The terrestrial heat flow value is one of the important geothermal geological indexes for the occurrence and distribution of geothermal resources. The terrestrial heat flow value in the southeast coastal area of China is obviously affected by the tectonic background (Hu and Wang, 1994; Zhang et al., 2020): The overall tectonic activity in the northwestern basin is relatively stable, showing a lower terrestrial heat flow value, while the southeastern coast shows a higher terrestrial heat flow value. As the overall heat flow value gradually increases from west to east, the local deep thermal structure controls the ground heat flow value. As shown in Figure 2, Fujian and Yangjiang-Maoming reached the highest ground heat flow value, reaching more than 95 μW/m3. Radiogenic heat is one of the main sources of heat in the lithosphere. U, Th and natural radio isotope 40K are the main thermal elements. According to the Zhao Ping et al. (1995), most of the southeast coastal areas are in the range of high heat generation rate, basically higher than 2.1 μW/m3. Especially in the whole province of Guangdong, southern Jiangxi and southern Fujian, the large area of granite outcrops, the heat generation rate background is more than 2.8 μW/m3, such a large area of high heat generation rate area is rare. Radioactive heat generation is the main source of crustal heat generation, which is of great significance to the occurrence of dry hot rock resources. The southeast coastal area, especially Huizhou and Sanhui of Guangdong Province, has a good prospect of exploration and development of geothermal resources.
Based on produced heat rate of typical magmatic rocks in Guangdong, Jiangxi and southeast of Fujian, previous scholars indicated that the radioactive heat production rate of most granitic rocks is greater than the average heat production rate of the Earth’s crust, signifying that the Yanshanian granite radioactive heat production have great contribution to regional heat (Zhao et al., 1995). Meanwhile, the southeastern coastal area has basically the same tectonic evolution pattern in Mesozoic (Wang et al., 2013), which is controlled by Paleo-Pacific plate subduction during the Yanshanian period and experienced large-scale magmatic events.
Based on seismic inversion in the area, it is clear that in this area and the surrounding granite areas, based on P wave inversion of deep seismic profiles, shows 5–10 km depth of granite. A reasonable speculation in the area is that Yashaninan granites have great potential for the geothermal reservoirs. Moreover, the Moho surface in the whole southeast region is relatively shallow and the temperature of the Moho surface is relatively high (Zhang et al., 2018). The deep granite can promote the heat transfer of the mantle to the surface. In summary, Yanshaninan granite in the southeast coastal area is of great significance to the geothermal energy system.
Conclusion
1) The age of granite reservoirs are mainly 94.7–138.4 Ma, mainly late Yanshanian period. The geochemical characteristics of rocks show they are formed in the subduction zone signature, indicating they are formed related to the extensional environment formed by retractable subduction of the ancient Pacific plate.
2) Deep Yanshanian granites promote mantle heat transfer to the surface and increase the proportion of mantle heat flow in the heat flow value. Yanshanian granite has high content of Th, U, and K, and its radioactive decay heat generation provides heat source for the high-temperature geothermal in the study area. The combined action of the two forms the high-temperature geothermal energy in the study area.
3) The deep granite in the study area has huge potential for geothermal resources. The geothermal genesis model in the study area has reference significance for the understanding of geothermal in the southeast region or the calculation of geothermal reserves.
Data Availability Statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author Contributions
JL is responsible for experimental design, HZ and YZ are responsible for providing overall ideas, JF is responsible for data analysis, YZ and MW is responsible for instrument operation.
Funding
This study was financially supported by National Key Research and Development Project (No. 2019YFC0604903), National Natural Science Foundation of China (Grant No. U20B6001) and consultation research project of Chinese Academy of Engineering (No. 2019-XZ-35-04 and No. 2021-XZ-16-01).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Bertani, R. (2012). Geothermal Power Generation in the World 2005-2010 Update Report. Geothermics 41, 1–29. doi:10.1016/j.geothermics.2011.10.001
Bertani, R. (2016). Geothermal Power Generation in the World 2010–2014 Update Report. Geothermics 60, 31–43. doi:10.1016/j.geothermics.2015.11.003
Brown, D. W., Duchane, D. V., Heiken, G., and Hriscu, V. (2012). Mining the Earth’s Heat: Hot Dry Rock Geothermal Energy. Heidelberg: Springer Science & Business Media. doi:10.1007/978-3-540-68910-2
Deng, J., Feng, Y., Di, Y., Cui, L., Xiao, Q., and Guo, S. S. (2016). The Instrusive Spatial Temporal Evolutional Framework in the Southeast China. Geol. Rev. 62 (1), 3–16.
Deng, Y. F., Li, J. T., and Peng, T. P. (2019). Lithospheric Structure in the Cathaysia Block (South China) and its Implication for the Late Mesozoic Magmatism. Physics of the Earth and Planetary Interiors, 24–34.
Frost, B. R., Barnes, C. G., Collins, W. J., Arculus, R. J., Ellis, D. J., and Frost, C. D. (2001). A Geochemical Classification for Granitic Rocks. J. Pet. 42 (11), 2033–2048. doi:10.1093/petrology/42.11.2033
Gao, W., Wang, Z., Li, C., and Wang, D. (2014). Zircon U-Pb Geochronological, Geochemistry and Tectonic Implication of Lndosinian Granite from southeastern Zhejiang, South China. Acta Geologica Sinica 88 (6), 1055∼1067.
Genter, A., Traineau, H., Dezayes, C., Elsass, P., and Villemin, T. (1995). Fracture Analysis and Reservoir Characterization of the Granitic Basement in the HDR Soultz Project (France). Geothermal Sci. Tech. 4 (3), 189–214.
Goldstein, B. A., Hill, A. J., and Long, A. (2008). Hot Rocks in Australia National Overview. ASEG Extended Abstr. 1, 1.
Gong, J., and John Chen, Y. (2014). Evidence of Lateral Asthenosphere Flow beneath the South China Craton Driven by Both Pacific Plate Subduction and the India-Eurasia continental Collisin. Terra Nova 26 (1), 55–63. doi:10.1111/ter.12069
Hasterok, D., and Chapman, D. S. (2007). Continental Thermal Isostasy: 2. Application to North America. J. Geophys. Res. Solid Earth 112 (B6). doi:10.1029/2006jb004664
Hu, S., He, L., and Wang, J. (2000). Heat Flow in the Continental Area of China: A New Data Set. Earth Planet. Sci. Lett. 179 (2), 407–419. doi:10.1016/s0012-821x(00)00126-6
Hu, S., and Wang, J. (1994). Heat Flow Characteristics of Orogenic Belts in southeastern China. Geol. Rev. 40 (5), 387-394.
Huang, S. P. (1992). Variations of Heat Flow and Crustal Thickness in the Continental Area of China. Chin. J. Geophys. 35 (4), 441–450.
Kelemen, P. B., Hanghøj, K., and Greene, A. R. (2003). One View of the Geochemistry of Subduction - Related Magmatic Arcs, with an Emphasis on Primitive Andesite and Lower Crust. Treatise Geochem. 3, 659.
Kuang, J., Qi, S., and Wang, S. (2020). Granite Intrusion in Huizhou, Guangdong Province and its Geothermal Implications. Earth Sci. 45 (2020), 1466–1480.
Li, D. W., and Wang, Y. X. (2015). Major Issues of Research and Development of Hot Dry Rock Geothermal Energy. Earth Sci. 40 (11), 1858–1869.
Li, J., Zhang, Y., Dong, S., and Johnston, S. T. (2014). Cretaceous Tectonic Evolution of South China: a Preliminary Synthesis. Earth-Science Rev. 134, 98–136. doi:10.1016/j.earscirev.2014.03.008
Li, L., Luo, J., Chen, Q., Shi, D., Yun, J., and Li, T. (2019). Discovery of Silurian schist from Well J6 in Chagan Sag, Yin’e Basin: Evidence from Zircon U-Pb Geochronology. Geology in China 46 (1), 209–201
Li, X. H., Li, Z. X., Li, W., Ying, L., Yuan, C., Wei, G., et al. (2007). U– Pb Zircon, Geochemical and Sr-Nd-Hf Isotopic Constraints on Age and Origin of Jurassic I- and A-Type Granites from Central Guangdong, SE China: A Major Igneous Event in Response to Foundering of a Subducted Flat-Slab? Lithos 96 (1-2), 186–204. doi:10.1016/j.lithos.2006.09.018
Li, Z.-X., and Li, X.-H. (2007). Formation of the 1300-km-wide Intracontinental Orogen and Postorogenic Magmatic Province in Mesozoic South China: A Flat-Slab Subduction Model. Geol. 35 (2), 179. doi:10.1130/g23193a.1
Liu, K., Li, Z., Xu, W., Ye, H., Zhao, X., Hu, Y., et al. (2016). Geochemical Characteristics and Metallogenic Regularity of Mesozoic Magmatic Rocks in South China [J]. J. Mineralogy, Pet. Geochem. 35 (6), 1141–1155.
Liu, Y. S., Hu, Z. C., Gao, S., Günther, D., Xu, J., Gao, C. G., et al. (2008). In Situ analysis of Major and Trace Elements of Anhydrous Minerals by LA-ICP-MS without Applying an Internal Standard. Chem. Geology 257 (1-2), 34–43. doi:10.1016/j.chemgeo.2008.08.004
Ludwig, K. R. (2003). Isoplot 3.00: A Geochronological Toolkit for Microsoft Excel. Berkeley: Berkeley Geochronology Center Special Publication, 4–70.
Lund, J. W., and Boyd, T. L. (2016). Direct Utilization of Geothermal Energy 2015 Worldwide Review. Geothermics 60, 66–93. doi:10.1016/j.geothermics.2015.11.004
Pan, Z., Zhang, Q., Chen, G., Jiao, S., Du, X., Miao, X., et al. (2017). The Relationship between the Mesozoic Magmatic Activity and Plate Subduction in Eastern China: The Relation between the Zhe-Fujian and the Japanese Arc and the Andean Arc and its Significance [J]. Acta Petrologica Sinica 33 (5), 1507–1523.
Pearce, J. A., Harris, N. B. W., and Tindle, A. G. (1984). Trace Element Discrimination Diagrams for the Tectonic Interpretation of Granitic Rocks. J. Pet. 25 (4), 956–983. doi:10.1093/petrology/25.4.956
Richards, H. G., Savage, D., and Andrews, J. N. (1992). Granite-water Reactions in an Experimental Hot Dry Rock Geothermal Reservoir, Rosemanowes Test Site, Cornwall, U.K. Appl. Geochem. 7 (3), 193–222. doi:10.1016/0883-2927(92)90038-5
Rudnick, R. L., and Fountain, D. M. (1995). Nature and Composition of the Continental Crust: A Lower Crustal Perspective. Rev. Geophys. 33 (3), 267. doi:10.1029/95rg01302
Rybach, L., and Buntebarth, G. (1984). The Variation of Heat Generation, Density and Seismic Velocity with Rock Type in the Continental Lithosphere. Tectonophysics 103 (1-4), 335–344. doi:10.1016/0040-1951(84)90095-7
Rybach, L., 1988. Determination of Heat Production Rate. Handbook of Terrestrial Heat Flow Density Determinations.
Sclater, J. G., Jaupart, C., and Galson, D. (1980). The Heat Flow through Oceanic and Continental Crust and the Heat Loss of the Earth. Rev. Geophys. 18 (1), 269–311. doi:10.1029/rg018i001p00269
Sun, S.-s., and McDonough, W. F. (1989). Chemical and Isotopic Systematics of Oceanic Basalts: Implications for Mantle Composition and Processes. Geol. Soc. Lond. Spec. Publications 42 (1), 313–345. doi:10.1144/gsl.sp.1989.042.01.19
Tan, H., Lin, F., Wen, J., Li, J., and Fan, Y. (2020). Application of Wide Field Electromagnetic Method in the Exploration of Deep Dry Hot Rocks in Hainan Province. China Pet. Chem. Stand. Qual. (4), 150–152.
Wang, D., and Shen, W. (2003). Granite Genesis and Crustal Evolution in southeastern China. Earth Sci. Front. 10 (3), 209–220.
Wang, S. J., Hu, S. B., and Wang, S. J. (1999). The Geothermal Effect of Radioactive Heat Generation and its Signif⁃icance to Hydrocarbon Maturation in Tarim Basin. Pet. Exploration Dev. 26 (5), 36–38.
Wang, Y., Fan, W., Zhang, G., and Zhang, Y. (2013). Phanerozoic Tectonics of the South China Block: Key Observations and Controversies. Gondwana Res. 23 (4), 1273–1305. doi:10.1016/j.gr.2012.02.019
Wollenberg, H. A., and Smith, A. R. (1987). Radiogenic Heat Production of Crustal Rocks: An Assessment Based on Geochemical Data. Geophys. Res. Lett. 14 (3), 295–298. doi:10.1029/gl014i003p00295
Wu, F., Li, X., Yang, J., and Zheng, Y. (2007). Discussions on the Petrogenesis of Granites [J]. Acta Petrologica Sinica 23 (6), 1217–1238.
Xi, Y., Wang, G., Liu, S., Zhao, Y., and Hu, X. (2018). The Formation of a Geothermal Anomaly and Extensional Structures in Guangdong, China: Evidence from Gravity Analyses. Geothermics 72, 225–231. doi:10.1016/j.geothermics.2017.11.009
Xu, X. S., O’Reilly, S. Y., Griffin, W. L., Wang, X., Pearson, N., and He, Z. (2007). The Crust of Cathaysia: Age, Assembly and Reworking of Two Terranes. Precambrian Res. 158 (1-2), 51–78. doi:10.1016/j.precamres.2007.04.010
Yuan, Y. S., Ma, Y. S., Hu, S. B., Guo, T., and Fu, X. (2006). Present - Day Geothermal Characteristics in South China. Chin. J. Geophys. 49 (4), 1118–1126. doi:10.1002/cjg2.922
Zhang, J., Wang, B. Y., Tang, X. C., and Dong, M. (2018). Temperature Structure and Dynamic Background of Crust and Mantle beneath the High Heat Flow Area of the South China Continental Margin. Chin. J. Geophys. 61 (10), 3917–3932.
Zhang, Y., Luo, J., and Feng, J. (2020). Characteristics of Geothermal Reservoirs and Utilization of Geothermal Resources in the Southeastern Coastal Areas of China. Journal of Groundwater Science and Engineering 8 (2), 134–142.
Zhao, P., Wang, J., and Wang, J. A. (1995). Characteristics of Heat Production Distribution in SE China. Acta Petrologica Sinica 11 (3), 292–305.
Keywords: southeast coastal areas, granite geothermal reservoir, geothermal resources, geochemistry, geochrohology
Citation: Zheng H, Luo J, Zhang Y, Feng J, Zeng Y and Wang M (2021) Geological Characteristics and Distribution of Granite Geothermal Reservoir in Southeast Coastal Areas in China. Front. Earth Sci. 9:683696. doi: 10.3389/feart.2021.683696
Received: 22 March 2021; Accepted: 30 July 2021;
Published: 13 August 2021.
Edited by:
Shu Jiang, The University of Utah, United StatesCopyright © 2021 Zheng, Luo, Zhang, Feng, Zeng and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jun Luo, luojun.syky@sinopec.com