A PRELIMINARY STUDY ON THE MECHANICS AND TECTONIC RELATIONSHIP TO THE GEOTHERMAL FIELD OF THE HEYUAN FAULT ZONE IN GUANGDONG PROVINCE
-
摘要: 广东省河源断裂带位于中国东南沿海地热异常区,地热资源十分丰富,但其形成机制和利用前景尚不确定。为此文章开展了多学科综合分析,获得以下初步认识:温泉是断裂带内深循环地下水被地温加热而成,断裂剪切热和花岗岩浆残余热的贡献基本可以排除;沿断裂展布的厚层硅化带是古水热活动的产物;硅化带形成时期的挤压应力方向为北东-南西,与河源断裂及河源盆地晚白垩世以来的伸展活动对应,现代构造应力场为北西西-南东东方向挤压,与古应力场相比发生了明显变化;现今构造应力场使得北东向河源断裂呈右旋挤压运动,而北西向断裂则发生左旋张剪,导致地下水循环格局也发生相应改变;目前温泉沿河源断裂呈带分布,沿北西向断裂呈线性溢出,断裂交汇部位是热泉上升的主要通道。总体而言,河源地区拥有令人鼓舞的地热资源及应用前景,有可能达到建设地热发电厂的目标。建议继续深入开展地质学、地球物理、水文地质学和地热成因机理等多学科综合研究,从而更好地定量评价地热潜能与开发前景。Abstract: Located in the geothermal anomaly area along the southeast coast of China, the Heyuan fault zone in Guangdong Province owns abundant geothermal resources. To investigate its potential, the mechanics and tectonic relationship to the geothermal field of the Heyuan fault zone in Guangdong Province are analysed. It is preliminarily believed that:1) The origin of hot springs is mainly attributed to deep circulation of groundwater along the fault zone heated by the geothermal gradient; the contribution of shear heat and residual heat of granite magma can be eliminated. 2) The thick quartz reef distributed along the fault is the product of the ancient hydrothermal activity. 3) The direction of compressive stress during the formation of quartz reef is NE-SW, corresponding to the extension of the Heyuan fault and the Heyuan Basin since the Late Cretaceous. This is a marked change to the stress direction of the modern tectonic stress field of NWW-SEE. 4) The current tectonic stress field makes the NE Heyuan faults dextral trans-compressional, while the NW faults have left-lateral strike and tension, resulting in a corresponding change in the groundwater circulation pattern. 5) The hot springs are distributed along the Heyuan fault zone and aligned along the NW faults; the intersection of the faults is the channel through which hot springs rise. Overall, the Heyuan area has promising geothermal resources and potential to establish geothermal power plants. It is recommended that further multidisciplinary studies are carried out, including geology, geophysics, hydrogeology and geothermal.
-
Key words:
- Heyuan fault /
- geothermal /
- quartz reef /
- tectonic stress field /
- groundwater circulation
-
近年来, 地球空间信息学和对地观测学迅猛发展[1], 随着计算机图形技术、三维建模技术和可视化技术的迅猛发展, 地下三维可视化建模技术得到了日益广泛的应用。长期以来在传统的地质现象表达中, 通常是以二维平面图和剖面图来表达野外地质的成果[2], 这种方式存在着空间信息的损失与失真、制图过程繁杂及信息更新困难。地质体三维可视化建模技术正是针对传统的地质信息模拟与表达方法的缺陷, 借助计算机和科学计算可视化技术, 直接从三维空间的角度去理解和表达地质体与地质环境。它具有形象、直观、准确、动态、丰富等特点[3], 因此, 三维可视化研究有了越来越重要的现实意义和实用价值。在涌现出的众多地质三维可视化软件中, 处在世界领先地位的Micromine公司开发的Micromine软件在这些方面已取得了可喜的成果。本文利用国际矿业界通用的Micromine软件, 以东昆仑造山带1:5万水泥厂幅为例, 探索应用三维地质实体建模的技术和方法。
1. 研究区地质概况
格尔木南1: 5万水泥厂幅位于东昆仑造山带东段, 大地构造上位于昆中断裂带与昆南断裂带之间的东昆仑南地体[4]。该地区山势雄伟, 相对高差达1000m。该区出露的岩石地层主要有:中-新元古界万保沟群大理岩和砾状灰岩; 奥陶-志留系纳赤台群水泥厂组灰岩、石灰厂组火山岩和碳酸岩、哈拉巴依沟组千枚岩和砂岩等; 下泥盆统牦牛山组砾岩夹砂岩及流纹岩; 上石炭统浩特洛洼组灰岩、板岩和砂岩; 中下三叠统洪水川组砂岩、板岩和砾岩; 新生代沉积沿雪水河和格尔木河谷分布。由于加里东期和印支期强烈造山运动, 上述地层被近东西向的逆冲断层分割成大小不等的构造岩片, 原始沉积层序已被打乱, 形成了一套总体有序、局部无序的构造地层体(图 1)。
2. Micromine软件简介
Micromine软件是澳大利亚Micromine国际矿业软件有限公司开发的从事地质勘探、资源评估、储量计算及露天和地下矿山开采设计的地理信息系统(GIS)三维立体空间平台软件。
软件已经汉化, 适用于WINDOWS95和WINDOWSNT及以上的平台, 自引进中国以来, 在国内多个省(区)的地勘单位、矿山企业取得了较好的使用效果, 受到矿业系统的好评; 软件的记录格式符合中国的格式要求, 具有和其他软件的接口, 使该系统的数据可被其他数据库管理系统及相关软件查询和编辑, 如可以导入CAD、Mapgis、Excel等格式的数据; 软件能够实现各种工程和矿体的三维立体显示和成图, 并根据地质统计学的原理提供进行矿体品位估值和储量计算的各种方法, 同时还提供了矿山开采设计与境界优化等功能; 软件可以制作演示工作成果的三维电影, 以利于进行工作汇报、成果展示和招投标等工作。
3. 三维建模的基本概念及原理
在三维建模的过程中, 地表模型的构建依托于数字地形模型(DTM)的生成, 地下模型的构建实际上是各地质体三维线框的生成以及相互之间布尔运算的结果。
3.1 DTM
DTM (Digital Terrain Model)-数字地形模型是利用一个任意坐标系中大量选择的已知x、y、z的坐标点对连续地面的一个简单的统计表示, 或者说, DTM就是地形表面形态属性信息的数字表达, 是带有空间位置特征和地形属性特征的数字描述。地形表面形态的属性信息一般包括高程、坡度、坡向等。
DTM系统主要是由计算机程序实现的, 应用于各种类型的计算机系统的DTM已经在许多国家开发成功, 尽管使用的方法不同、用户界面各异, 但主要功能都是从离散数据构建出相互连接的网络结构, 以此作为地形的数字模型基础。目前, DTM已经成为GIS的重要组成部分[5]。
数字地形模型通过野外测量、航测、地图数字化等途径可得到设计所需的地形资料, 这些资料经过计算机处理后可以形成数字地形模型[6~7]。通常数字地形模型建立的方法, 一是根据现在测量仪器(全站仪)得到的第一手地形测量数据(通常是点坐标)按照一定格式输入到相关软件后可以直接通过软件功能生成数字地形模型; 另一种是根据该地区提供的地质地形图进行矢量化后导入到软件中进行高程赋值, 然后直接生成数字地形模型。
3.2 三维线框实体基本原理
三维线框模型的构建主要是采用了TIN技术(不规则三角网模型)中的Voronoi图与Delaunay三角形算法。TIN是一种表示数字高程模型(DEM)的方法, 它既减少规则格网方法带来的数据冗余, 同时在计算效率(如坡度)方面又优于纯粹基于等高线的方法。这种表面模型可以进行体积估算、表面渲染、切制剖面、快速三维显示等操作, 比计算机图形学中的表面模型有所扩展, 能满足地质勘探中地质制图的基本要求, 也是建立三维实体模型的基础。地质体的地质形态复杂多变, 很难用规则的几何体来描述, 它需要一种灵活、简便、快速的方法来建立地质体的不规则几何模型, 这种模型正是基于这一需要而提出来的。
TIN不规则三角网技术 不规则三角网TIN属于场的数据模型范畴, 三角形单元遍布整个区域无重叠、无缝隙, 实体间的空间关系通过三角形连接表达。用于支持TIN建立的几何构造-Delaunay三角网因其特有的“外接圆规则”和“最大最小角规则”, 使得其成为空间邻近表达的有力工具。这些模型可以看作是不规则三角网TIN在平面空间场表达上发挥的具体作用。
Voronoi图 基于场论观念, 空白区域可看作是目标的相互“竞争”生成空间, 在群结构中, 每个面目标的“影响区域”通过竞争以中剖分结果达到平衡, 这是Voronoi图几何构造的基本思想。Voronoi图又叫泰森多边形或Dirichlet图, 它是由1组连接2邻点直线的垂直平分线组成的连续多边形组成。N个在平面上有区别的点, 按照最邻近原则划分平面, 每个点与它的最近邻区域相关联。Delaunay三角形是由与相邻Voronoi多边形共享1条边的相关点连接而成的三角形。Delaunay三角形的外接圆圆心是与三角形相关的Voronoi多边形的1个顶点。Voronoi三角形是Delaunay图的偶图[8]。
Delaunay三角形网的通用算法-逐点插入算法 基于散点建立数字地面模型, 常采用在d维的欧几里得空间中构造Delaunay三角形网的通用算法-逐点插入算法。具体算法过程包括:①遍历所有散点, 求出点集的包容盒, 得到作为点集凸壳的初始三角形并放入三角形链表; ②将点集中的散点依次插入, 在三角形链表中找出其外接圆包含插入点的三角形(称为该点的影响三角形), 删除影响三角形的公共边, 将插入点同影响三角形的全部顶点连接起来, 从而完成一个点在Delaunay三角形链表中的插入; ③根据优化准则对局部新形成的三角形进行优化(如互换对角线等), 将形成的三角形放入Delaunay三角形链表; ④循环执行上述第2步, 直到所有散点插入完毕。
3.3 线框模型布尔运算
线框模型布尔运算在复杂地质体建模中主要是对有相交关系的复杂地质体(实体与实体之间、面与实体之间)进行交、差、并等, 基本的布尔运算功能对实体进行交接拼合以得到建模所需的组合实体模型。
目前大多数数字矿山工程软件都可以实现实体模型的并、交和差运算, 表面模型可以进行差、并运算。在建模实际中, 按照软件提供的布尔运算功能进行组合, 能够满足复杂地质体的建模要求。复杂地质体建模中主要遇到的问题是地质体的分支复合和地质体被断层穿插切割等地质现象很难翔实地反映在可视化模型中。在建模中, 实际上2个实体剖面之间是按照直线的方式来进行三角网连接的, 但是遇到呈曲面的断层或者断层破碎带, 实际上这种地质体与断层实体之间的吻合关系是很难在建模实体中反映的。所以只有通过线框模型布尔运算得到的实体组合来逼近地质体实体间的切割和相交等空间形态和关系。
4. 1: 5万水泥厂幅三维地质模型的建立
对于水泥厂幅三维地质模型的建立主要有3方面内容, 数字地形模型(DTM)的生成以及地下地质体模型的构建, 进而将二者无缝结合一体化显示。
4.1 数字地形模型(DTM)
数字地形模型主要由栅格和TIN技术2种数据格式表示。该地质体地表模型的构建采用了TIN技术, 因为TIN模型可以根据地形的具体特点和复杂程度而确定采样点的密度和位置, 可以避免数据的冗余, 并且可以相对准确的描述出地形的结构和一些局部地区, 同时又便于进行地形计算和分析, 能够较好地反映实际地形信息[9]。具体步骤包括:①将1: 5万水泥厂幅地质地形图在Mapgis中进行矢量化(*.wl), 并压缩存盘。操作:打开Mapgis-输入编辑-矢量化-压缩存盘; ②把矢量化生成的Mapgis线文件(*.wl)转换输出为Mapgis明码格式文件(*.wal)。操作:图形处理-文件转换-装入线-输出Mapgis明码格式; ③将Mapgis明码格式文件(*.wal)导入到Micromine中, 转换成Micromine中的线文件(*. str)。操作:在Micromine菜单栏中, 文件-导入-Mapgis文件-线串; ④根据地形图等高线值为线串的每条等高线进行高程赋值。操作:视图管理器-线(加载)-属性窗口-Z值; ⑤对等高线线文件进行坐标矫正。操作:测量-坐标系统转换-平面; ⑥生成DTM, 其空间数据插值主要是通过TIN技术来实现的。操作:工具栏-创建DTM-选择线(图 2a); ⑦生成的DTM根据研究的需要可直接披覆该区域的地质图、地形图以及遥感图像等。操作:视图管理器-线框-属性-叠加图像。
4.2 地下地质体模型
剖面图的绘制 地下地质体模型的建立主要有2种数据, 钻孔数据和剖面图数据, 这里采用的是剖面图数据。根据水泥厂地质图以及地表DTM在软件中的坐标数据, 我们将由西到东绘制12幅地质剖面图, 具体步骤包括:①在Micromine中加载已经完成的地表DTM, 进行相应设置。操作:视图管理器-线框(地表DTM)-显示选项-二维切块; 工具栏-向西看; ②设置完成后, 软件中显示出地形线, 输入第一个剖面的东坐标。操作:显示设置-范围-东2200;③新建线文件, 对该剖面地形线进行捕捉矢量化, 选择合适的底面线, 生成闭合的剖面线轮廓。操作:工具栏-新线、创建线、捕捉线上点; ④依据相关资料, 对地质图中的地质体走向、产状、断层和地质界线等详细分析和合理的解释, 完成此幅剖面图的绘制, 并且保证剖面图中相同岩性的的岩石连接成一个封闭的曲线, 这样就形成反映地质体现状的剖面图(图 2b); ⑤同理, 重复1~4步骤, 完成其余11幅地质剖面图, 最后生成一个12幅地质剖面图的空间三维线文件(图 2c)。
生成线框模型 ①加载地质剖面图线文件, 新建线框, 按照地质体的趋势, 利用地质体轮廓线连接成面技术在相邻同种地质体剖面线之间用三角网连接三维地质体表面, 直到连接完毕, 最后根据具体情况选择将地质体两端封闭终点或者封闭到点, 这样就形成了地质体的实体。操作:工具栏-创建线框-新三角网-封闭终点(图 2c); ②同理连接完成所有该区域的地质体, 根据地质体颜色对照表分别赋予合适的地质体颜色(图 2d); ③对相交的地质体以及高出地表DTM地质体进行布尔运算, 布尔运算能够使不同实体间的吻合空间关系表现的淋漓尽致。当地质体构造很复杂的情况下, 布尔运算将是一个工程量相当大的工作, 而且要求建模人员格外仔细, 每进行一个运算时必须要对比两块岩体, 判断正确的运算结果。操作:线框-布尔运算(图 2e、2f); ④将地下地质体模型线框(3D)另存为2D的地质体模型线框, 用于截取剖面图时的二维剖面显示。操作:线框属性-显示选项-二维切块。
4.3 数字地形模型(DTM)和地下地质体模型的结合
在Micromine同一个三维视图中加载地表DTM和地下地质体线框(2D和3D均加载), 这样就实现了地表三维数据与地下三维地质体无缝集成三维一体化显示(图 2g)。由此, 按需要可截取任一范围、任一方向、任一比例尺的剖面图, 进一步可以进行矿产勘查分析和地质分析(图 2h)。
5. 基本认识
本文探讨了基于Micromine系统的三维可视化建模技术方法流程, 对于复杂地质体的建模, 地质断层、褶皱等复杂地质结构的表示方面仍然是地质体可视化建模研究的重点之一。本文所建立的1: 5万水泥厂幅三维地质模型, 再现了该地区地质体的三维空间形态和接触关系。但是, 深部地质体的分布与相互关系强烈地依赖于对深部地质的认知程度。尽管如此, 三维可视化建模技术在今后的地质应用中潜力巨大。
致谢: 感谢澳大利亚Micromine国际矿业软件及咨询集团对本项研究工作的大力支持。责任编辑:范二平 -
图 1 河源地区地质构造简图和剖面示意图(据广东省地质局、广东省有色地质环境中心内部资料以及野外实测综合)
1—第四系;2—古近系;3—白垩系;4—侏罗系;5—石炭—二叠系(局部少量三叠系);6—泥盆系;7—寒武奥陶系;8—晚元古界;9—侏罗纪花岗岩(局部白垩纪);10—三叠纪花岗岩;11—白垩纪末—古近纪初基性岩;12—地质界线;13—断层及其编号(粗细示断层规模,虚线代表隐伏断裂或线性构造);14—天然温泉;15—热水井;16—钻孔;17—水文取样点;18—取样编号;19—剖面位置;20—城镇
F1—河源断裂,F2—人字石断裂,F3—大坪—岩前断裂,F4—石角—新港—白田断裂,F5—南山—坳头断裂;A-A′—表示剖面位置
a—河源地区相关断裂、地层、主要花岗岩体、温泉和取样位置图;b—横切河源断裂带和河源盆地剖面示意图Figure 1. Overview geological map of the Heyuan area and a cross-section schematic
图 2 河源断裂带野外岩相变化示意图(位置见图 1b)
a、b、c、d分别表示观测点照片及对应位置;露头照片显示断裂带不同岩相(从花岗岩到石英岩脉)脆性变形和石英成分均增加
Figure 2. Cross-section schematic of the Heyuan Fault, showing the relative position of the fault zone facies, based on field observations
表 1 研究区内沿河源断裂带分布的温泉温度和流量
Table 1. Measurements of hot springs along the Heyuan Fault Zone.
点号 数据来源 样品点位 纬度 经度 描述 温度/℃ 流量/(L/s) L1 本次研究采集 Loc.5 23° 51′09″N 114° 47′31″E 自流井(废弃建筑旁) 55.7 4.8 L2 本次研究采集 Loc.5 (距L1点20 m) 23° 51′07″N 114° 47′25″E 当地村民开采井 57.9 2.9 L3 本次研究采集 Loc.6 24° 04′32″N 115° 08′57″E 水池 62.5 5.6 M4 Mao等[21] DG-13 24°04′29.6″N 115°08′59.4″E 泉点(位置同L3) 63.2 M5 Mao等[21] DG-14 23°51′11.2″N 114°47′29.4″E 泉点(位置同L1/L2) 56.7 M6 Mao等[21] DG-15 23°26′49.6″N 115°06′26.9″E 泉点 78.6 M7 Mao等[21] DG-16 23°12′10.5″N 114°21′31.3″E 泉点 59.8 注:表中以Loc.命名的取样点位见图 2a;以DG命名的取样点见文献[21] 表 2 河源断裂带的浅层钻孔测量的地下水温度和流量
Table 2. Groundwater measurements obtained from shallow boreholes along a section of the Heyuan Fault in the Badengcheng area.
井名 纬度 经度 深度/m 孔内最高温度/℃ 孔口温度/℃ 流量/(L/s) ZK1 23°37′35.43″N 114°37′7.33″E 150.81 30.1 27 6.48 ZK2 23°37′43.85″N 114°37′32.11″E 323.25 39.5 34 1.37 ZK3 23°38′8.52″N 114°37′47.01″E 350 43.3 37 0.75 ZK4 23°37′47.60″N 114°37′39.44″E 358 51.8 45.2 5.00 ZK6 23°37′48.00″N 114°37′38.81″E 375 52.0 45.2 13.31 ZK7 23°37′48.41″N 114°37′38.29″E 381.12 64.6 61.0 16.20 表 3 地质模型中的断层及其参数
Table 3. Faults considered in the geomodel, the parameters used, and comparison with other sources.
编号 断裂 走向 倾角及依据 切割深度 F1 河源断裂 NE5°~80° 45°(据[9])35°~50°(据[11])35°~62°(据[12])37°~49°平均43°(据野外观测点测量) -10 km F2 人字石断裂 NE40°~60° 60°(据[9]) -10 km F3 大坪—岩前断裂 NE30~50° 70°(据[9]) -10 km F4 石角—新港—白田断裂 NE25°~50° 70°(据[9]) -10 km F5 南山—坳头断裂 N85°E 90°(据[9]) -15 km F6 未命名断裂 N52°W 80°(据野外观测点测量) -10 km F7 未命名断裂 N60°W 90° (典型的安德森型高倾角走滑断裂) -10 km F8 推测断裂 N52°W 90° (典型的安德森型高倾角走滑断裂) -10 km F9 推测断裂 N23°E 75° (据地质填图显示的小断裂及地表线性构造) -10 km F10 推测断裂 N14°E 70° (据地质填图显示的小断裂及地表线性构造) -10 km -
[1] LUND J W, BOYD T L. Direct utilization of geothermal energy 2015 worldwide review[C]//Proceedings World Geothermal Congress. Melbourne, Australia, 2015. https://www.sciencedirect.com/science/article/pii/S0375650511000344 [2] BERTANI R. Geothermal power generation in the world 2010-2014 update report[J]. Geothermics, 2016, 60:31-43. doi: 10.1016/j.geothermics.2015.11.003 [3] 陈墨香.中国地热资源的分布及其开发利用[J].自然资源, 1992, 7(3):40-46, 58. http://d.old.wanfangdata.com.cn/Conference/6896043CHEN Moxiang. Geothermal distribution and utilities in China[J]. Advances of Earth Sciences, 1992, 7(3):40-46, 58. (in Chinese) http://d.old.wanfangdata.com.cn/Conference/6896043 [4] 汪集旸, 胡圣标, 庞忠和, 等.中国大陆干热岩地热资源潜力评估[J].科技导报, 2012, 30(32):25-31. doi: 10.3981/j.issn.1000-7857.2012.32.002WANG Jiyang, HU Shengbiao, PANG Zhonghe, et al. Estimate of geothermal resources potential for hot dry rock in the continental area of China[J]. Science & Technology Review, 2012, 30(32):25-31. (in Chinese with English abstract) doi: 10.3981/j.issn.1000-7857.2012.32.002 [5] WAN T F. The tectonics of China:data, maps and evolution[M]. Berlin Heidelberg:Springer, 2012. [6] WANG A D, SUN Z X, HU B Q, et al. Guangdong, a potential province for developing hot dry rock geothermal resource[J]. Applied Mechanics and Materials, 2014, 492:583-585. doi: 10.4028/www.scientific.net/AMM.492 [7] LIU R X, XIE G H, ZHOU X H, et al. Tectonic environments of cenozoic volcanic rocks in china and characteristics of the source regions in the mantle[J]. Chinese Journal of Geochemistry, 1995, 14(4):289-302. doi: 10.1007/BF02872628 [8] 王霄飞, 余珊, 龚跃华, 等.华南北东向断裂在南海北部陆架的延伸[J].大地构造与成矿学, 2014, 38(03):557-570. http://d.old.wanfangdata.com.cn/Periodical/ddgzyckx201403006Wang Xiao Fei, Yu Shan, Gong Yue Hua, et al. Extension of NE-trending faults in south china to northern south china sea continental shelf[J]. Geotectonica et Metallogenia, 2014, 38(3):557-570(in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/ddgzyckx201403006 [9] CHENG H H, ZHANG H, ZHU B J, et al. Finite element investigation of the poroelastic effect on the Xinfengjiang reservoir-triggered earthquake[J]. Science China Earth Sciences, 2012, 55(12):1942-1952. doi: 10.1007/s11430-012-4470-8 [10] CHEN L, TALWANI P. Reservoir-induced seismicity in China[J]. Pure and Applied Geophysics, 1998, 153(1):133-149. doi: 10.1007/s000240050188 [11] QIU X, FENTON C. Factors controlling the occurrence of reservoir-induced seismicity[C]//Lollino G. Engineering Geology for Society and Territory. Cham: Springer, 2015, 6: 567-570. doi: 10.1007%2F978-3-319-09060-3_102 [12] LEE C F, YE H, ZHOU Q. On the potential seismic hazard in Hong Kong[J]. Episodes, 1997, 20(2):89-94. http://cn.bing.com/academic/profile?id=8c2a56bd34fe1c1d813736b273b3158a&encoded=0&v=paper_preview&mkt=zh-cn [13] 刘大任.邵武-河源断裂带活动性及分段评价[J].地质力学学报, 1997, 3(2):54-60. http://journal.geomech.ac.cn/ch/reader/view_abstract.aspx?flag=1&file_no=19970221&journal_id=dzlxxbLIU Daren. Segmentation of the Shaowu Heyuan fault zone and their activity assessment[J]. Journal of Geomechanics, 1997, 3(2):54-60. (in Chinese with English abstract) http://journal.geomech.ac.cn/ch/reader/view_abstract.aspx?flag=1&file_no=19970221&journal_id=dzlxxb [14] 邹和平, 彭樊源, 苏章歆, 等.河源伸展剥离断层(博罗-龙川段)及其第四纪活动特征[J].华南地震, 2010, 30(S1):1-9. http://d.old.wanfangdata.com.cn/Periodical/hndz2010z1002ZOU Heping, PENG Fanyuan, SU Zhangxin, et al. Discussions on the Heyuan extensional detachment fault from Boluo to Longchuan and its quaternary activities[J]. South China Journal of Seismology, 2010, 30(S1):1-9. (in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/hndz2010z1002 [15] REGENAUER-LIEB K, VEVEAKIS M, POULET T, et al. Stimulating granites: from synchrotron microtomography to enhancing reservoirs[C]//Proceedings World Geothermal Congress 2015. Melbourne, Australia, 2015. [16] QIU X L, WANG Y, WANG Z Z, et al. Determining the origin, circulation path and residence time of geothermal groundwater using multiple isotopic techniques in the Heyuan Fault Zone of Southern China[J]. Journal of Hydrology, 2018, 567:339-350. doi: 10.1016/j.jhydrol.2018.10.010 [17] HU S B, HE L J, WANG J Y. Heat flow in the continental area of China:A new data set[J]. Earth and Planetary Science Letters, 2000, 179(2):407-419. doi: 10.1016/S0012-821X(00)00126-6 [18] WANG G, LI K, WEN D, et al. Assessment of geothermal resources in China[C]//Thirty-Eighth Workshop on Geothermal Reservoir Engineering. California: Stanford University, Stanford, 2013, 10. [19] 田春艳.广东省中高温地热资源勘查与开发利用建议[J].地下水, 2012, 34(4):61-63. http://d.old.wanfangdata.com.cn/Periodical/dixs201204023TIAN Chunyan. Suggestions on the exploration and development of high temperature geothermal resources in Guangdong province[J]. Groundwater, 2012, 34(4):61-63. (in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/dixs201204023 [20] XI Y, WANG Y, HU X, et al. Geothermal structure revealed by Curie isotherm surface in Guangdong province[C]//International Workshop and Gravity, Electrical & Magnetic Methods and their Applications. Chengdu, China, 2015: 189-192. [21] MAO X M, WANG Y X, ZHAN H B, et al. Geochemical and isotopic characteristics of geothermal springs hosted by deep-seated faults in Dongguan Basin, Southern China[J]. Journal of Geochemical Exploration, 2015, 158:112-121. doi: 10.1016/j.gexplo.2015.07.008 [22] LACHENBRUCH A H, SASS J H. Heat flow and the thermal regime of the crust[C]//HEACOCK J G. The Earth's Crust, Its Nature and Physical Properties. Washington, D. C: American Geophysical Union, 1977: 626-675. [23] VIGNERESSE J L, CUNEY M. Are granites representative of heat flow provinces[A]//AČG ERMÁK V, RYBACH L. Terrestrial heat flow and the lithosphere structure[M]. Berlin: Springer, 1991. [24] SUN Z X, WANG A D, LIU J H, et al. Radiogenic heat production of granites and potential for hot dry rock geothermal resource in Guangdong province, Southern China[C]//Proceedings World Geothermal Congress 2015. Melbourne, Australia, 2015. [25] RYBACH L. Determination of heat production rate[C]//HÄNEL R, RYBACK L, STEGENA L. Handbook of Terrestrial Heat Flow Density Determination. Dordrecht: Kluwer, 1988, 125-142. [26] ZHANG Y, YANG J H, SUN J F, et al. Petrogenesis of Jurassic fractionated I-type granites in Southeast China:Constraints from whole-rock geochemical and zircon U-Pb and Hf-O isotopes[J]. Journal of Asian Earth Sciences, 111:268-283. doi: 10.1016/j.jseaes.2015.07.009 [27] BIRCH F, ROY R F, DECKER E R. Heat flow and thermal history in New York and New England[C]//ZEN F A, WHITE W S, HADLEY J B, et al. Studies of Appalachian Geology: Northern and Maritime. New York: Interscience Jr, 1968: 437-451. [28] WEBB P C, LEE M K, BROWN G C. Heat flow-heat production relationships in the UK and the vertical distribution of heat production in granite batholiths[J]. Geophysical Research Letter, 1987, 14(3):279-282. doi: 10.1029/GL014i003p00279 [29] 钟建强, 周蒂.华南沿海温泉分布与地震活动关系初探[J].华南地震, 1990, 10(4):22-29. http://www.cnki.com.cn/Article/CJFDTOTAL-HNDI199004004.htmZHONG Jianqiang, ZHOU Di. A preliminary study of the relationship between the distribution of hot springs and the activity of earthquakes along south China coast[J]. South China Journal of Seismology, 1990, 10(4):22-29. (in Chinese with English abstract) http://www.cnki.com.cn/Article/CJFDTOTAL-HNDI199004004.htm [30] SIBSON R H. Fault rocks and fault mechanisms[J]. Journal of the Geological Society, 1977, 133(3):191-213. doi: 10.1144/gsjgs.133.3.0191 [31] SIBSON R H. Continental fault structure and the shallow earthquake source[J]. Journal of the Geological Society, 1983, 140(5):741-767. doi: 10.1144/gsjgs.140.5.0741 [32] ZHAO B, BAI Z M, XU T, et al. Lithological model of the South China crust based on integrated geophysical data[J]. Journal of Geophysics and Engineering, 2013, 10(2):25005. http://cn.bing.com/academic/profile?id=ece1d32a35e290b3acf364900bb119c5&encoded=0&v=paper_preview&mkt=zh-cn [33] 广东省地质矿产局.广东省区域地质志[M].北京:地质出版社, 1988.Geological Bureau of Guangdong Province. The geology of Guangdong Province[M]. Beijing:Geological Press, 1998. (in Chinese) [34] 饶春涛, 李平鲁.珠江口盆地热流研究[J].中国海上油气(地质), 1991, 5(6):7-18. http://d.old.wanfangdata.com.cn/Periodical/dxqy200504029RAO Chuntao, LI Pinglu. Study of heat flux in basins around the Pearl river delta[J]. China Offshore Oil and Gas (Geology), 1991, 5(6):7-18. (in Chinese with English abstract) http://d.old.wanfangdata.com.cn/Periodical/dxqy200504029 [35] 王妙月, 杨懋源, 胡毓良, 等.新丰江水库地震的震源机制及其成因初步探讨[J].地球物理学报, 1976, 17(1):1-17. http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX197601000.htmWANG Miaoyue, YANG Maoyuan, HU Yuliang, et al. Mechanism of the reservoir impounding earthquakes at Xinfengjiang and a preliminary endeavour to discuss their cause[J]. Acta Geophysica Sinica, 1976, 19(1):1-17. (in Chinese with English abstract) http://www.cnki.com.cn/Article/CJFDTOTAL-DQWX197601000.htm [36] 陈伟光.华南沿海沉积盆地的新构造运动及其与地震的关系[J].华南地震, 1995, 15(2):55-61. http://www.cnki.com.cn/Article/CJFDTOTAL-HNDI502.008.htmCHEN Weiguang. On the relation between earthquake and neotectonic movement of depositional basins in coastal area of south China[J]. South China Journal of Seismology, 1995, 15(2):55-61. (in Chinese with English abstract) http://www.cnki.com.cn/Article/CJFDTOTAL-HNDI502.008.htm [37] CHADWICK R A, LEONARD R B. Structural controls of hot-spring systems on southwestern Montana[M/OL]. USGS Open-File Report 79-1333. U.S. Geological Survey, 1979. https://pubs.usgs.gov/of/1979/1343/report.pdf. [38] CUREWITZ D, KARSON J A. Structural settings of hydrothermal outflow:Fracture permeability maintained by fault propagation and interaction[J]. Journal of Volcanology and Geothermal Research, 1997, 79(3-4):149-168. doi: 10.1016/S0377-0273(97)00027-9 [39] PERSON M, HOFSTRA A, SWEETKIND D, et al. Analytical and numerical models of hydrothermal fluid flow at fault intersections[J]. Geofluids, 2012, 12(4):312-326. doi: 10.1111/gfl.2012.12.issue-4 [40] 中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会.地热资源地质探勘规范: GB/T 11615-2010[S].北京: 中国标准出版社, 2011.General administration of quality supervision, inspection and quarantine of the People's Republic of China, Standardization administration. Code for geological exploration of geothermal Resources: GB/T 11615-2010[S]. Beijing: Standards Press of China, 2011. (in Chinese) -