Geological Characteristics and Distribution of Granite Geothermal Reservoir in Southeast Coastal Areas in 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  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 km 2 , 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.
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.
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/m 2 , while that in the central region is lower than 60 mW/m 2 . 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/m 2 , locally up to 220 mW/m 2 ). The average heat flow in northern Jiangsu Basin is 71.2 mW/m 2 , while that in southern Jiangsu and southern Anhui areas is 62.7 mW/m 2 . 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/m 2 , while that in southern Yunnan is 76.7 mW/m 2 . 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/m 2 , and that in the south is 61.6 mW/m 2 . The average heat flow value of 64.2 mW/m 2 in southern China is higher than that of 61 mW/m 2 in other mainland parts in China (Wang et al., 1990), which is close with that of the globe for 65 mW/m 2 (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 km 2 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 km 2 . 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).
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.
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.

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.
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 206 Pb/ 238 U 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 206 Pb/ 238 U age of 112.92 ± 0.85 Ma (MSWD 1.6), representing the formation Frontiers in Earth Science | www.frontiersin.org August 2021 | Volume 9 | Article 683696 6 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 206 Pb/ 238 U 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     (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). 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.

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 SiO 2 /(Na 2 O + K 2 O) ( 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 SiO 2 − FeO T /(FeO T + MgO) and SiO 2 − Al 2 O 3 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 SiO 2 /(Na 2 O + K 2 O), 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 SiO 2 -FeO T /(FeO T + MgO) and SiO 2 -Al 2 O 3 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 SiO 2 /(Na 2 O + K 2 O), 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 SiO 2 -FeO T /(FeO T + MgO) and SiO 2 -Al 2 O 3 diagram ( Figures 9A,B), most samples are plotted in the POG (post-Orogenic Grabite) zone (Frost et al., 2001). On the Nb-Y and 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 faultblock-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.

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/m 2 , which is higher than the average terrestrial heat flow value of 64.2 mW/m 2 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 238 U, 232 Th, and 40 K (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/m 3 and the maximum value up to 9.15 μW/m 3 . 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/m 3 (Rybach, 1988), 3.3-3.8 μW/m 3 (Wollenberg and Smith, 1987) and 3.0 μW/m 3 (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 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/m 3 . Radiogenic heat is one of the main sources of heat in the lithosphere. U, Th and natural radio isotope 40 K 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/m 3 . 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/m 3 , 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 largescale 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.