地球物理学报  2020, Vol. 63 Issue (7): 2697-2709   PDF    
共和盆地干热岩地热资源的成因机制:来自岩石放射性生热率的约束
张超1,2, 胡圣标3, 宋荣彩1,2, 左银辉1,2, 姜光政3, 雷玉德4, 张盛生5, 王朱亭3     
1. 成都理工大学能源学院, 成都 610059;
2. 油气藏地质及开发工程国家重点实验室(成都理工大学), 成都 610059;
3. 中国科学院地质与地球物理研究所岩石圈演化国家重点实验室, 北京 100029;
4. 青海省环境地质勘查局, 西宁 810007;
5. 青海省水文地质工程地质环境地质调查院, 西宁 810008
摘要:共和盆地位于青藏高原东北缘,是我国重要的干热岩地热资源赋存区之一.成因机制研究是干热岩地热资源研究中最为基础与核心的工作之一,也是干热岩地热资源潜力精细评价和合理开发利用的重要依据.本文基于98块采自共和—贵德盆地岩石样品的放射性生热率数据,分析了盆地主要岩性岩石的放射性生热率分布特征,讨论了共和盆地干热岩地热资源的成因机制,并初步建立了干热岩地热资源的成因模式.研究表明,共和盆地恰卜恰地热区沉积岩(以泥岩和粉砂质泥岩为主)放射性生热率为1.21~2.02 μW·m-3,平均值为1.67±0.29 μW·m-3;以花岗岩和花岗岩闪长岩等为主的基底岩石的放射性生热率介于1.17~5.81 μW·m-3之间,平均值为3.20±1.07 μW·m-3.贵德盆地扎仓寺地热区沉积岩(以砂岩和泥质砂岩为主)放射性生热率为1.83~2.40 μW·m-3,平均值为2.13±0.23 μW·m-3;基底花岗质岩石的放射性生热率介于0.92~6.49 μW·m-3之间,平均值为2.81±1.40 μW·m-3.测试数据显示共和—贵德盆地基底花岗岩不存在高放射性生热率异常.但是,新生代以来印度—欧亚板块持续性陆-陆碰撞作用造成的壳内放射性生热元素富集层增厚,导致了花岗岩放射性生热率的热贡献量同步增大(30.3~40.5 mW·m-2),因此,花岗岩放射性生热可为共和盆地干热岩地热资源提供稳定的壳内热源基础.基于放射性生热率数据和热流配分研究,结合研究区已有地质-地球物理研究资料,本文认为壳内部分熔融层作为附加热源为共和盆地干热岩地热资源提供了重要的附加热流贡献.在此基础上,本文初步构建了共和盆地干热岩地热资源成因模式:加厚地壳花岗岩放射性生热与壳内部分熔融层供热.
关键词: 花岗岩      放射性生热率      部分熔融      干热岩      共和盆地      贵德盆地     
Genesis of the hot dry rock geothermal resources in the Gonghe basin: constraints from the radiogenic heat production rate of rocks
ZHANG Chao1,2, HU ShengBiao3, SONG RongCai1,2, ZUO YinHui1,2, JIANG GuangZheng3, LEI YuDe4, ZHANG ShengSheng5, WANG ZhuTing3     
1. College of Energy, Chengdu University of Technology, Chengdu 610059, China;
2. State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Chengdu University of Technology, Chengdu 610059, China;
3. State Key Laboratory of Lithospheric Evolution, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China;
4. Environmental Geological Prospecting Bureau of Qinghai Province, Xining 810007, China;
5. The Survey Institute of Hydrogeology, Engineering Geology and Environment Geology of Qinghai Province, Xining 810008, China
Abstract: The Gonghe Basin, located in the northeastern margin of the Tibetan Plateau, is one of the important hot dry rock (HDR) geothermal resources occurrence area in China. Among the researches of the HDR geothermal resources, the study of the genesis mechanism is one of the most basic and core works and it is also the important basis for the precise potential evaluation and rational development and utilization of the HDR geothermal resources. Based on the radiogenic heat production rate (RHPR) data measured on the 98 core samples gathered from the Gonghe-Guide basin, the characteristics of the RHPR of the main lithology rock in the basin are analyzed, the genesis mechanism of the HDR geothermal resources in the Qabqa geothermal area of the Gonghe basin is discussed, and the corresponding genesis model of HDR geothermal resources is constructed in this study. The results show that the RHPR of the sedimentary rocks (mainly mudstone and silty mudstone) of the Qabqa geothermal area in the Gonghe basin ranges from 1.21 to 2.02 μW·m-3, with an arithmetic average of 1.67±0.29 μW·m-3, while the RHPR of the basal granites, which are mainly composed of granite and granodiorite, ranges from 1.17 to 5.81μW·m-3, with an arithmetic average of 3.20±1.07μW·m-3. For the Zhacangsi geothermal area in the Guide basin, the RHPR of the sedimentary rocks (sandstone and muddy sandstone) ranges from 1.83 to 2.40 μW·m-3, 2.13±0.23 μW·m-3 on average, while it varies from 0.92 to 6.49 μW·m-3(the average of which is 2.81±1.40μW·m-3) for the basal granitic rocks. The tests show that there is no high RHPR anomaly of the basal granites in the Gonghe-Guide basin. However, the thickening of the radioactively enriched layer caused by the continuous collision between the Indian and Eurasian plates since Cenozoic leads to the synchronous increase of the heat contribution of the RHPR of granites, which is about 30.3~40.5 mW·m-2, and therefore it can serve as a stable heat source for the HDR geothermal resources in the Gonghe basin. Based on the RHPR data and the heat flow partition analysis, as well as the existing geological-geophysical research data in this area, this work suggests that the partial melting layer in the crust beneath the study area may be an additional heat source, providing important additional heat flow contribution to the HDR geothermal resources in the study area. On the basis above, the genesis model of the HDR geothermal resources in the Gonghe basin is preliminarily constructed in this study: the joint effect of the radiogenic heat contribution of the granite in the thickened crust and the heat contribution of the partial melting layer within the crust.
Keywords: Granite    Radiogenic heat production rate    Partial melting    Hot dry rock    Gonghe basin    Guide basin    
0 引言

地热资源是一种清洁、稳定、空间分布广泛且储量巨大的可再生能源.根据其不同成因和产出条件, 地热资源可以分为水热型地热资源和干热岩型地热资源(汪集旸等, 2012).作为现阶段各个国家主要的开发利用对象, 水热型地热资源仅占已探明地热资源中的约10%, 超过90%的可利用的地热资源仍储存于干热岩型地热资源中.保守估计全球地壳3~10 km深度所蕴含的干热岩资源储量相当于全球所有石油、天然气和煤炭所蕴藏能量的30倍(许天福等, 2012).自1973年美国开展全球首个干热岩试验项目(Fenton Hill试验项目)以来, 英国、法国、德国、瑞士、日本、澳大利亚和冰岛等国均投入到干热岩地热资源的研究中(Breede et al., 2013).迄今, 全球共建立了47个增强型地热系统工程项目(王贵玲等, 2020).成因机制的研究是干热岩地热资源研究中最为基础与核心内容之一, 也是干热岩地热资源潜力精细评价和合理开发利用的重要依据.目前, 国际上主要干热岩试验场地多选址于与板块边界相关的地区和新生代裂谷区, 这些地区具有高大地热流特征, 如美国Fenton Hill试验场地大地热流为160 mW·m-2(张杨, 2016), Geysers试验场地大地热流超过168 mW·m-2(Stanley and Blakely, 1995), 位于新生代莱茵地堑的法国Soultz干热岩试验场地大地热流亦高达176 mW·m-2(Baria et al., 1999).

大地热流, 作为地球内部热状态在地表最为直观的指标, 是由地壳岩石中放射性生热元素(铀、钍、钾)衰变产生的热量(即地壳热流)和地幔热流两部分组成(汪集旸, 2015Furlong and Chapman, 2013Jiang et al., 2019).热流配分关系是大地热流研究的重要延伸和扩展, 影响着深部地温场的分布、地壳和上地幔的活动性.地热异常区的热流配分关系是认识地热系统热源构成及其热贡献的重要窗口, 而岩石放射性生热率研究正是热流配分关系研究的基础, 同时也是地热系统成因机制研究的主要内容之一(Abbady and Al-Ghamdi, 2018Zhang et al., 2018).

共和盆地地处青藏高原东北缘, 是中-新生代以来发展起来的菱形断陷盆地.盆地内部已发现水温超过15 ℃的热水或地热异常点达84处之多, 其中6处温泉出口温度大于60 ℃, 扎仓寺温泉和曲乃海温泉温度均超过当地沸点, 分别为93.5 ℃和96.6 ℃(薛建球等, 2013康维海, 2016).近年来, 由中国地质调查局和青海省国土资源厅组织实施的多口钻井(DR3、DR4、GR1和GR2)均揭示出高温现象.测温资料显示DR3、DR4和GR1等钻井在2700 m左右深度地层温度均超过150 ℃, 2886 m处温度最高可达到180 ℃(张超等, 2018Zhang et al., 2018, 2019).基于钻井温度和岩石热导率数据计算的大地热流值表明, 共和盆地大地热流值介于93.3~111.0 mW·m-2之间, 平均值为102.2 mW·m-2, 远高于我国大陆地区平均大地热流值60.4 mW·m-2(姜光政等, 2016张超等, 2018).高温温泉、高温钻井以及高大地热流值一致表明共和盆地具有丰富的干热岩地热资源.目前, 围绕共和盆地干热岩岩体开展了一系列地球物理探测(如航磁测量、天然地震背景噪声层析成像和大地电磁测深等)、地热测量和数值模拟工作, 在大地热流特征、干热岩岩体空间展布以及干热岩地热资源潜力评价等方面形成了一定认识和成果(赵振等, 2009李永革, 2016张超等, 2018张盛生等, 2019岳高凡等, 2015Zhang et al., 2018Hou et al., 2018Xu et al., 2018).但是, 目前仍未形成关于共和盆地干热岩地热资源成因机制的共同认识(孙知新等, 2011王斌等, 2015严维德, 2015李林果和李百祥, 2017).

本文在系统分析共和盆地大地热流特征(张超等, 2018)的基础上, 通过对共和盆地及邻区贵德盆地内主要钻井的岩芯样品开展放射性生热率测试, 分析了研究区主要岩性岩石的放射性生热率特征, 进一步讨论了共和盆地干热岩地热资源的热源机理与成因模式.该研究成果不仅可为共和盆地干热岩地热资源潜力精细评价和合理开发利用提供重要的理论支撑和科学依据, 也可丰富我国干热岩地热资源成因机制的理论认识.

1 地质背景

共和—贵德盆地, 位于青藏高原东北缘, 地处秦岭—祁连—昆仑褶皱系中东昆仑与西秦岭造山带的交接转换部位, 是一呈北西西向展布的菱形山间盆地(石宝颐和张峻太, 1982).盆地四周被断褶带隆起山地围限, 北侧为青海南山和拉鸡山断褶隆起带, 南侧为河卡南山和巴吉山断褶隆起带, 东西两侧分别为扎马山断褶隆起带和鄂拉山构造岩浆带, 瓦里贡山构造岩浆带将盆地划分为共和、贵德两盆地(王斌等, 2015).受盆地边界断裂活动的影响, 挽近期周边山脉不断隆起上升, 全盆地内部相对沉降, 沉积了一套巨厚的新生代沉积层, 主要由第四系冲洪积、河湖相堆积物和新、古近系湖相堆积物等组成(图 1, 薛建球等, 2013).盆地基底主要由三叠纪地层和以花岗岩、花岗闪长岩和斑状花岗岩为主的侵入岩组成.张宏飞等(2006)通过锆石U-Pb LA-ICP-MS年代学手段对区内黑马河岩体和温泉岩体进行定年研究, 指出共和盆地周边花岗岩年龄为218±2~235±2 Ma(即印支期), 并基于主微量元素和Pb-Sr-Nd同位素地球化学研究指出花岗质岩石的原岩来自下地壳的变玄武岩类.

图 1 共和—贵德盆地地质简图及主要钻井位置分布图 Fig. 1 Map shows the major tectonic structures and borehole location of the Gonghe—Guide basin

研究区内的主要地热区为位于共和盆地东北部的恰卜恰地热区(图 1c)和贵德盆地西南部的扎仓寺地热区(图 1d).张森琦等(2018a)基于1:5万重磁和航磁测量指出恰卜恰地热区具有三套地层结构:浅部为约600 m厚的第四纪地层, 中部为古近纪-新近纪地层, 厚度约为800 m, 下部为花岗岩基底段, 基底花岗岩的埋藏深度大致为1200 ~1350 m, 与基于可控源音频大地电磁测深获得的基底埋深情况(1000~1400 m, 薛建球等, 2013)较为接近.扎仓寺地热区位于NNW向压扭性逆断层-热光断裂(F5, 图 1d)与NE向张扭性正断层-扎仓断裂(F2, 图 1d)的交汇部位(郎旭娟等, 2016).大气降水入渗到地下沿扎仓断裂径流被加热, 遇热光断裂下盘泥岩受阻而上升出露地表是该区地热异常和地热资源形成的主要原因(李小林等, 2016Jiang et al., 2018, 2019).

2 样品测试

以共和盆地东北部恰卜恰地热区和贵德盆地西南部扎仓寺地热区现有钻井为依托(图 1), 共采集了98块岩石样品进行岩石放射性生热率测试.其中52块样品采自恰卜恰地热区GR1、GR2和DR3钻井, 46块样品采集自扎仓寺地热区ZR2钻井.岩石放射性生热率的测试分为岩石密度测试和岩石放射性生热元素(铀、钍、钾)含量测试两部分, 密度测试部分依托于中国科学院地质与地球物理研究所岩石热物性实验室, 岩石放射性生热元素含量测试由核工业北京地质研究院分析测试研究中心完成, 其中铀和钍元素含量采用等离子质谱仪(ICP-MS)测定, 钾元素含量利用石墨炉原子吸收分析仪测定.

3 测试结果与讨论 3.1 共和—贵德盆地岩石放射性生热率

基于上述测定的放射性生热元素含量, 结合岩石对应的密度, 可按照如下公式(Rybach, 1976)计算得出岩石放射性生热率:

(1)

式中, A为放射性生热率, 单位μW·m-3ρ为岩石密度, 单位g·cm-3CUCThCK分别为放射性元素铀、钍、钾的含量.

图 2显示出共和盆地恰卜恰地热区GR1钻井岩芯样品放射性生热率与放射性生热元素Th/U比值随深度变化情况.以泥岩和粉砂质泥岩为代表的沉积岩的放射性生热率介于1.66~2.02 μW·m-3, 平均值为1.79 μW·m-3(表 1), Th/U比值介于2.34~6.56之间, 平均值为4.91;基底花岗质岩石的放射性生热率随深度没有明显的关系, 总体介于1.73~4.48 μW·m-3, 平均值为3.04 μW·m-3, Th/U比值介于1.37~16.15之间.从图 2可以看出, 1341 m处岩石样品的Th/U比值明显偏离整体趋势, 表现出高异常现象(值为16.15), 测试数据表明该样品中的U元素含量较低(仅为1.43×10-6), 可能由测试误差引起.若剔除该测试样品, 则基底花岗质岩石样品的平均Th/U比值为4.43.GR2钻井岩芯样品的放射性生热率与放射性生热元素Th/U比值随深度变化情况如图 3所示, 3个沉积岩样品的放射性生热率介于1.21~1.90 μW·m-3, 平均值为1.50 μW·m-3, 略小于GR1钻井沉积岩的放射性生热率, Th/U比值介于3.72~4.94之间, 平均值为4.41.与GR1钻井相类似, GR2钻井中基底花岗质岩石的放射性生热率随深度亦未呈现明显的关系, 总体介于1.17~4.92 μW·m-3, 平均值为2.98 μW·m-3, 与GR1钻井的基底花岗岩放射性生热率基本持平, Th/U比值介于0.32~7.18之间, 平均值为3.55.

图 2 共和盆地GR1钻井岩芯样品放射性生热率及放射性生热元素Th/U比值随深度分布图 (钻井岩性剖面图修改自张森琦等, 2018a) Fig. 2 Map shows that the distribution of RHPR and Th/U ratio with depth of the core samples of the GR1 borehole in the Gonghe basin (the borehole lithology profile is modified from Zhang et al., 2018a)
表 1 共和—贵德盆地岩石放射性生热元素含量和放射性生热率统计表 Table 1 Statistical table of radioactive element contents and RHPR of the rocks samples gathered from the Gonghe-Guide basin
图 3 共和盆地GR2钻井岩芯样品放射性生热率及放射性生热元素Th/U比值随深度分布图 Fig. 3 Map shows that the distribution of RHPR and Th/U ratio with depth of the core samples of the GR2 borehole in the Gonghe basin

同样, 作者对采自DR3钻井的5个花岗质岩石样品的测试结果做了类似分析, 结果显示其放射性生热率介于2.12~5.81 μW·m-3, 平均值为4.70 μW·m-3, Th/U比值介于1.83~6.25之间, 平均值为3.14(表 1).总体看来, 恰卜恰地热区沉积岩样品的放射性生热率介于1.21~2.02 μW·m-3, 平均值为1.67±0.29 μW·m-3;基底火成岩放射性生热率统计直方图(图 4)显示, 其放射性生热率大致呈现出正态分布, 所有样品生热率值介于1.17~5.81 μW·m-3之间, 算术平均值为3.20±1.07 μW·m-3, 略大于全球中-新生代(0~250 Ma)花岗岩放射性生热率平均值3.09±1.62 μW·m-3(Artemieva et al., 2017), 与西藏南部花岗岩放射性生热率3.2±1.5 μW·m-3(n=12, 沈显杰等, 1989)相类似, 但大于藏北花岗岩放射性生热率2.6±1.6 μW·m-3(n=24, 沈显杰等, 1989).

图 4 共和盆地恰卜恰地热区花岗质岩石样品放射性生热率分布直方图 Fig. 4 Histogram of the RHPR of the granitic rock samples in the Qiabuqia geothermal area, Gonghe basin

贵德盆地扎仓寺地热区ZR2钻井岩芯样品放射性生热率随深度变化情况如图 5a所示, 岩石放射性生热率随埋藏深度的关系不大, 所有岩石样品的放射性生热率介于0.92~6.49 μW·m-3之间, 平均值为2.58±1.21 μW·m-3.其中, 砂岩和泥质砂岩等沉积岩的放射性生热率相对较小, 其值介于1.83~2.40 μW·m-3之间, 平均为2.13±0.23 μW·m-3.花岗质岩芯样品的放射性生热率统计直方图(图 5b)表明:与恰卜恰地热区基底火成岩相比(图 4), 该区花岗质岩石样品的放射性生热率正态分布性较差, 其放射性生热率介于0.92~6.49 μW·m-3之间, 平均值为2.81±1.4 μW·m-3, 较共和盆基底火成岩放射性生热率低, 但大于Liu等(2017)基于贵德盆地其他钻井的岩芯样品所获得的放射性生热率(2.04 μW·m-3, 表 1).

图 5 贵德盆地ZR2钻井岩芯样品的放射性生热率统计图 (a)岩芯样品放射性生热率-深度分布图;(b)花岗质岩石样品放射性生热率分布直方图. Fig. 5 Map shows the RHPR of the core samples from the ZR2 borehole in the Guide basin (a) The RHPR versus depth of the core samples; (b) Histogram of the RHPR of the granitic rock samples.
3.2 共和盆地干热岩地热资源成因机制 3.2.1 花岗岩放射性生热率对共和盆地热异常的启示

共和盆地恰卜恰地热区大地热流值介于93.3~111.0 mW·m-2之间, 平均值高达102.2 mW·m-2(张超等, 2018), 远高于我国大陆地区平均大地热流值60.4 mW·m-2(姜光政等, 2016), 也高于青藏高原南部喜马拉雅和拉萨—冈底斯块体平均大地热流(分别为90.8和82.1 mW·m-2, Jiang et al., 2019), 属于典型的地热异常区.花岗岩高放射性生热率一直被认为是该区地热异常和干热岩地热资源形成的主要原因(李林果和李百祥, 2017), 但本文所开展的岩石放射性生热率测试结果表明:基底花岗岩放射性生热率平均值为3.20±1.07 μW·m-3, 基本持平于全球中-新生代花岗岩的放射性生热率, 未出现高异常现象.与澳大利亚中部库珀盆地的Habanero干热岩试验场地的Innamincka花岗岩(其放射性生热率为7~10 μW·m-3, Chen and Wyborn, 2009)相比, 该区花岗岩的放射性生热率同样偏小.采自贵德盆地扎仓寺地热区的同期花岗岩测试数据再次证实该区基底花岗岩的放射性生热率维持在正常水平.

大地热流由来自地壳内部岩石的放射性生热元素衰变所产生的热量(即地壳热流)和地幔热流两部分组成.管彦武等(2012)基于青藏高原现有大地热流数据和多条地震剖面计算的青藏高原地区地幔热流分布指出, 班公—怒江缝合带、印度—雅鲁藏布江缝合带和拉萨地体存在地幔热流较高区域, 部分地区地幔热流可达260 mW·m-2, 而青藏高原东北缘地区地幔热流普遍较低.Zhang等(2018)基于贵德盆地现有大地热流和岩石热物性数据, 指出贵德盆地地幔热流为27 mW·m-2, 与管彦武等(2012)的研究结果较为接近.地幔热流是地下深部热状态的直接表现, 其高低反映着区域构造活动性的强弱, 因此地幔热流具有区域性影响特征(Francheteau et al., 1984张林友等, 2016).因此, 可取共和盆地地幔热流为贵德盆地地幔热流(即27 mW·m-2), 则共和盆地地壳岩石放射性生热量应约为75 mW·m-2.

基于上述岩石放射性生热率测试数据, 取共和盆地花岗岩放射性生热率为3.2 μW·m-3, 地壳放射性生热元素富集层厚度取15~20 km(汪集旸等, 2012), 放射性生热元素垂向分布遵循指数衰减模式(Furlong and Chapman, 2013), 则由壳内花岗岩放射性生热率贡献的热流为30.3~40.5 mW·m-2, 约占该区大地热流的29.6%~39.6%.由此可以得出, 尽管该区花岗岩的放射性生热率不高, 但是新生代印度—欧亚板块持续性陆-陆碰撞作用造成的壳内放射性生热元素富集层增厚, 导致了花岗岩放射性生热率的热贡献量同步增加, 故盆地加厚地壳花岗岩的放射性生热可以为该区地热异常提供稳定的壳内热源基础.

通过上述分析, 如不考虑其他附加热源供热情况, 由共和盆地下地壳岩石放射性生热贡献的热流理应为34.5~44.7 mW·m-2.若取共和盆地下地壳岩石放射性生热率为全球大陆下地壳放射性生热率平均值, 0.68 μW·m-3(Furlong and Chapman, 2013), 则盆地下地壳厚度理应高达51~66 km, 然而这个理论厚度完全超过了共和盆地整个地壳的厚度(51~56 km, Jia et al., 2019), 这显然是不合理的.因此, 本文倾向于该区还存在附加热源的额外热流贡献.同时, 来自共和盆地温泉和地热井地下水的地球化学同位素数据表明, 3He/4He值介于(4.65~9.52)×10-8之间, 反映出热量主要来自地壳部分(Hou et al., 2018).因此, 研究区高大地热流可能来自加厚地壳花岗岩放射性生热和壳内局部附加热源热贡献的双重作用.

3.2.2 壳内部分熔融层的地球物理证据

为探测共和盆地干热岩地热资源潜力及圈定潜在干热岩靶区, Gao等(2018)首次利用大地电磁测深方法获得了该区地壳50 km以浅的三维电阻率分布结构(图 6).Gao等(2018)指出共和盆地下距地表 15~35 km深度区间存在明显的低电阻率异常层, 并将该低电阻率异常层解译为部分熔融层(图 6).基于塔秀—共和东西向天然地震背景噪声层析成像和大地电磁测深勘查剖面结果, 张森琦等(2018a)同样指出共和盆地下距地表 8~32 km深度区间存在东西向长约50 km的部分熔融层(即低阻-高导层).但是, 由于地球内部的物质组成与状态的复杂性, 地球物理勘探方法所获得的信息往往具有多解性和不确定性(郭进京和周安朝, 2008), 因此, 需要与其他研究结果相互佐证.地壳温度信息是判断地壳中部分熔融存在与否的最为直接和关键的证据.张超等(2018)基于共和盆地大地热流和岩石热物性数据, 利用一维稳态热传导方程估算的深部温度数据显示, 共和盆地在14~15 km深度上地层温度达到含水花岗岩的熔融温度(660 ℃~670 ℃), 在19~20 km深度处达到干花岗岩的熔融温度(890 ℃~900 ℃), 表明该区20 km以浅或已发生部分熔融.实际上, 不仅是共和盆地, 青藏高原北部多条地球物理测深剖面均揭示出在地壳范围内存在显著的低速-高导层, 这类低速-高导层也被普遍认为是地壳岩石发生部分熔融后在岩石物理性质上的表现(Royden et al., 1997Clark and Royden, 2000Zhang et al., 2004Le Pape et al., 2012Li et al., 2014, 2019Zheng et al., 2016Wu et al., 2017Wang et al., 2018Deng et al., 2019), 同时也得到了岩石地球化学研究数据的进一步证实(Wang et al., 2012).

图 6 共和盆地大地电磁测点分布(a)及各剖面电阻率分布结构(b—f) (数据引自Gao et al., 2018) Fig. 6 Maps show that the Magnetotelluric sites within the Gonghe basin (a) and the resistivity structures (b—f) along the profiles indicated on (a) (The data was cited from Gao et al., 2018)

因此, 基于岩石放射性生热率测试数据, 结合共和盆地已有大地热流、深部地温场、大地电磁测深等研究资料, 本文倾向于大地电磁测深所显示的低电阻率异常(低阻-高导层)现象很可能是壳内部分熔融层的反映.深部地温场预测的部分熔融层深度与大地电磁剖面所呈现的低电阻率异常层(低阻-高导层)深度之间的差异可能来源于:(1)温度场计算过程中的基本假设条件-稳态热传导, 从而忽略了部分熔融层的瞬态热传导甚至热对流对地温场在空间三维和时间四维尺度上造成的影响;(2)用简单的一维代替复杂但真实的三维, 从而忽略了地层结构和岩石热物性在空间三维上的各向异性.要精确分析部分熔融层对共和盆地地温场演化过程的影响以及部分熔融层形成的深部动力学机制, 需要更多地质、地球物理、地球化学以及低温热年代学约束条件, 这也是我们今后工作的重点方向.

3.2.3 共和盆地干热岩地热资源成因机制

综合上述资料, 本文将共和盆地干热岩地热资源可能的成因机制概括如图 7所示.

图 7 盆地尺度干热岩地热资源成因机制模式图 Fig. 7 Genesis mechanism model of the HDR geothermal resources in the basin scale

共和盆地干热岩地热资源的热源为加厚地壳花岗岩放射性生热和壳内部分熔融层供热.根据单次热事件的热松弛特征时间计算公式τd2/4κ(τ为时间, d为深度, κ为热扩散率, 这里取1 mm2·s-1, Turcotte and Schubert, 2014张超等, 2018), 此部分熔融层形成年龄应该小于9.7 Ma (d取35 km, 对应于大地电磁测深剖面上低电阻率异常层的最大深度).干热岩地热资源的“热储层”为基底印支期花岗岩, 花岗岩具有较低的孔渗性(岳高凡等, 2015庄亚芹等, 2017), 由钻井所揭示的花岗岩埋藏深度约为900~1400 m, 钻井温度数据显示花岗岩层的温度梯度介于39~45.2 ℃/km之间, 平均值为41.3 ℃/km(张超等, 2018).基底花岗岩层之上的早中更新世的共和组、中新世的咸水河组与上新世的临夏组共同构成干热岩的盖层, 盖层厚度约900~1400 m, 盖层岩石以泥岩、砂质泥岩和砂岩为主, 且呈以段级地层为单位粗细相间的韵律性沉积(张森琦等, 2018a), 岩石具有低热导率特征(1.59±0.19 W·m-1·K-1, 张超等, 2018), 因此隔热保温效果较好, 能有效地防止深部热量快速散失.

3.3 我国干热岩地热资源勘探的启示

我国西南部地处地中海—喜马拉雅地热带, 东部邻近环太平洋地热带, 干热岩地热资源极为丰富.汪集旸等(2012)评价的我国大陆地区3~10 km深度区间干热岩地热资源储量显示, 干热岩地热资源总量达20.9×106 EJ, 若以2%作为可开采资源量系数, 相当于我国2010年能量消耗总量的4400倍.如何有效地勘探这类储量巨大的清洁可再生能源是影响我国干热岩地热资源开发利用的重要因素之一.干热岩地热资源成因机制研究作为干热岩研究中最为基础与核心的工作之一, 为更精确的圈定干热岩地热资源潜在靶区提供了可能.

在干热岩地热资源形成中, 热源是最为重要的控制因素之一.国际上许多干热岩试验场地均与岩浆热源有关, 如美国Fenton Hill和Geysers干热岩试验场地(Rutqvist et al., 2015Kelkar et al., 2016).我国著名的西藏羊八井和云南腾冲热海地热田也均与壳内岩浆有关(多吉, 2003Tian et al., 2019).本文基于共和—贵德盆地岩石放射性生热率测试数据, 结合已有大地热流、深部地温场和地球物理探测等资料, 将共和盆地干热岩地热资源的热源机制归纳为加厚地壳花岗岩放射性生热和壳内部分熔融层供热, 亦与岩浆热源有关.这类与岩浆热源相关的干热岩地热资源多赋存于近代火山区和构造强烈活动带.我国长白山、五大连池、云南腾冲、雷琼火山等近代火山区和青藏高原部分地区(如藏南、川西地区以及共和盆地)壳内存在异常热源—低速-高导层(廖志杰和张振国, 1989白登海等, 1994胡久常等, 2007汤吉等, 2001张健等, 2017张森琦等, 2018bGao et al., 2018), 具备干热岩地热资源产出的基础条件, 为我国潜在的岩浆成因的干热岩地热资源赋存区.

花岗岩高放射性生热亦可作为干热岩地热资源形成过程中的重要热源.一方面花岗岩高放射性生热可来自放射性生热元素(铀、钍和钾)异常富集导致的高放射性生热率, 以位于澳大利亚中部库珀盆地Habanero干热岩试验场地为代表, 其花岗岩放射性生热率高达7~10 μW·m-3(Chen and Wyborn, 2009), 约为全球上地壳岩石放射性生热率的3~5倍.另一方面壳内放射性生热元素富集层的增厚, 亦可造成花岗岩高放射性生热.本文热流配分研究指出, 尽管共和盆地基底花岗岩放射性生热率不高, 但其热贡献量仍可达30.3~40.5 mW·m-2, 主要原因是新生代以来印-欧大陆碰撞所导致的地壳增厚(郭晓玉等, 2017), 使得壳内放射性生热元素富集层厚度同步增加, 花岗岩放射性生热率的热贡献量也随之增加.这类高放射性生热成因的干热岩地热资源多赋予于地表及地壳浅部发育有大型的中新生代酸性花岗岩类岩体的地区和由于板块活动引起的地壳增厚区.我国东南沿海地区, 作为我国最主要的花岗岩分布区, 大面积地出露着具有高放射性生热率的花岗岩体(蔺文静等, 2016), 其燕山期花岗岩的放射性生热率最高可达到16.3 μW·m-3(马峰等, 2015), 为干热岩地热资源的形成提供了重要条件.此外, 基底起伏变化引起的侧向热折射(李春荣等, 2017)和构造活跃区的构造变形热(Tian et al., 2018)亦能引起局部地热异常.毫无疑问, 干热岩地热资源勘探及优质靶区圈定是一个复杂的系统性工程, 涉及大地构造学、构造地质学、岩石学、地球化学和地球物理学(重、磁、电、震、热)等多学科交叉融合, 本节旨在解剖共和盆地干热岩地热资源成因机制和蕴藏模式的基础上, 对我国今后干热岩地热资源勘探方向提出几点粗浅的看法.

4 结论

本文基于对共和—贵德盆地恰卜恰和扎仓寺地热区主要钻井岩芯样品的放射性生热率测试, 分析了共和盆地基底花岗岩的放射性生热率特征及其对干热岩地热资源的热贡献, 在此基础上, 结合该区现有研究成果, 初步构建了干热岩地热资源的成因模式.基于以上研究数据, 得出如下主要结论:

(1) 共和盆地恰卜恰地热区沉积岩(以泥岩和粉砂质泥岩为主)放射性生热率为1.21~2.02 μW·m-3, 平均值为1.67±0.29 μW·m-3;以花岗岩和花岗岩闪长岩等为主的基底花岗质岩石的放射性生热率介于1.17~5.81 μW·m-3, 平均值为3.20±1.07 μW·m-3, 略大于全球中-新生代花岗岩放射性生热率平均值.

(2) 贵德盆地扎仓寺地热区砂岩和泥质砂岩等沉积岩放射性生热率为1.83~2.40 μW·m-3, 平均值为2.13±0.23 μW·m-3;基底花岗质岩石的放射性生热率介于0.92~6.49 μW·m-3之间, 平均值为2.81±1.4 μW·m-3.

(3) 岩石放射性生热率数据表明共和—贵德盆地基底花岗质岩石并未出现高放射性生热率异常, 但新生代印-欧大陆碰撞造成的地壳放射性生热元素富集层增厚, 导致花岗岩放射性生热率的热贡献量同步增大(30.3~40.5 mW·m-2).因此, 花岗岩放射性生热可为共和盆地干热岩地热资源提供稳定的壳内热源基础.

(4) 综合分析热流配分结果和已有地质-地球物理资料, 初步认为共和盆地壳内部分熔融层作为附加热源为干热岩地热资源提供了重要的附加热流贡献.在此基础上, 构建了共和盆地加厚地壳花岗岩放射性生热与壳内部分熔融层供热的干热岩地热资源成因模式.

致谢  感谢审稿专家对本文提出的建设性意见和编辑部的大力支持!
References
Abbady A G E, Al-Ghamdi A H. 2018. Heat production rate from radioactive elements of granite rocks in north and southeastern Arabian shield Kingdom of Saudi Arabia. Journal of Radiation Research and Applied Sciences, 11(4): 281-290. DOI:10.1016/j.jrras.2018.03.002
Artemieva I M, Thybo H, Jakobsen K, et al. 2017. Heat production in granitic rocks:Global analysis based on a new data compilation GRANITE2017. Earth-Science Reviews, 172: 1-26. DOI:10.1016/j.earscirev.2017.07.003
Bai D H, Liao Z J, Zhao G Z, et al. 1994. The inference of magmatic heat source beneath the Rehai (Hot Sea) field of Tengchong from the result of magnetotelluric sounding. Chinese Science Bulletin, 39(7): 572-577.
Baria R, Baumgärtner J, Gérard A, et al. 1999. European HDR research programme at Soultz-sous-Forêts (France) 1987-1996. Geothermics, 28(4-5): 655-669. DOI:10.1016/S0375-6505(99)00036-X
Breede K, Dzebisashvili K, Liu X L, et al. 2013. A systematic review of enhanced (or engineered) geothermal systems:past, present and future. Geothermal Energy, 1(4): 1-27.
Chen D, Wyborn D. 2009. Habanero field tests in the Cooper Basin, Australia:A proof-of-concept for EGS. Transactions-Geothermal Resources Council, 33: 140-145.
Clark M K, Royden L H. 2000. Topographic ooze:Building the eastern margin of Tibet by lower crustal flow. Geology, 28(8): 703-706. DOI:10.1130/0091-7613(2000)28<703:TOBTEM>2.0.CO;2
Deng Y F, Li J T, Song X D, et al. 2019. The lithospheric-scale deformation in NE Tibet from joint inversion of receiver function and surface wave dispersion. Terrestrial, Atmospheric and Oceanic Sciences, 30(1): 127-137. DOI:10.3319/TAO.2019.01.18.03
Duo J. 2003. The basic characteristics of the Yangbajing geothermal field-a typical high temperature geothermal system. Engineering Science (in Chinese), 5(1): 42-47.
Francheteau J, Jaupart C, Jie S X, et al. 1984. High heat flow in southern Tibet. Nature, 307(5946): 32-36. DOI:10.1038/307032a0
Furlong K P, Chapman D S. 2013. Heat flow, heat generation, and the thermal state of the lithosphere. Annual Review of Earth and Planetary Sciences, 41: 385-410. DOI:10.1146/annurev.earth.031208.100051
Gao J, Zhang H J, Zhang S Q, et al. 2018. Three-dimensional magnetotelluric imaging of the geothermal system beneath the Gonghe Basin, Northeast Tibetan plateau. Geothermics, 76: 15-25. DOI:10.1016/j.geothermics.2018.06.009
Guan Y W, Guan Y, Gao R, et al. 2012. Crustal radiogenic heat and mantle geothermal status in tibetan plateau based on p-wave velocity. Journal of Jilin University (Earth Science Edition) (in Chinese), 42(2): 562-568, 574.
Guo J J, Zhou A C. 2008. Methods of delineating high temperature rock mass in geothermal abnormal area. Geothermal Energy (in Chinese), (4): 4-7.
Guo X Y, Gao R, Gao J R, et al. 2017. Integrated analysis on the deformation of the Liupan Shan fold-thrust belt, NE Tibet, and its tectonic attribution. Chinese Journal of Geophysics (in Chinese), 60(6): 2058-2067. DOI:10.6038/cjg20170603
Hou Z, Xu T, Jiang Z, et al. 2018. Geochemical and isotopic characteristics of geothermal water in the Gonghe Basin of the northeast Tibetan Plateau, China. PROCEEDINGS.//43rd Workshop on Geothermal Reservoir Engineering Stanford University. Stanford, California.
Hu J C, Bai D H, Wang W H, et al. 2007. Magnetotelluric surveying and electrical structure of the deep underground part in Leiqiong Volcanic area. South China Journal of Seismology (in Chinese), 27(1): 1-7.
Jia S X, Guo W B, Mooney W D, et al. 2019. Crustal structure of the middle segment of the Qilian fold belt and the coupling mechanism of its associated basin and range system. Tectonophysics, 770: 128154. DOI:10.1016/j.tecto.2019.06.024
Jiang G Z, Gao P, Rao S, et al. 2016. Compilation of heat flow data in the continental area of China (4th edition). Chinese Journal of Geophysics (in Chinese), 59(8): 2892-2910. DOI:10.6038/cjg20160815
Jiang Z J, Xu T F, Owen D D R, et al. 2018. Geothermal fluid circulation in the Guide Basin of the northeastern Tibetan Plateau:Isotopic analysis and numerical modeling. Geothermics, 71: 234-244. DOI:10.1016/j.geothermics.2017.10.007
Jiang Z J, Xu T F, Mallants D, et al. 2019. Numerical modelling of stable isotope (2H and 18O) transport in a hydro-geothermal system:Model development and implementation to the Guide Basin, China. Journal of Hydrology, 569: 93-105. DOI:10.1016/j.jhydrol.2018.11.065
Kang W H. 2016. Hot spring with an outlet temperature of 96.6℃ has been discovered in the Guide basin, Qinghai Province. Geothermal Energy (in Chinese), (1): 17.
Kelkar S, WoldeGabriel G, Rehfeldt K. 2016. Lessons learned from the pioneering hot dry rock project at Fenton Hill, USA. Geothermics, 63: 5-14. DOI:10.1016/j.geothermics.2015.08.008
Lang X J, Liu F, Liu Z M, et al. 2016. Terrestrial heat flow in Guide Basin, Qinghai. Geological Science and Technology Information (in Chinese), 35(3): 227-232.
Le Pape F, Jones A G, Vozar J, et al. 2012. Penetration of crustal melt beyond the Kunlun Fault into northern Tibet. Nature Geoscience, 5(5): 330-335. DOI:10.1038/ngeo1449
Li C R, Rao S, Hu S B, et al. 2017. Present-day geothermal field of the Jiaoshiba shale gas area in southeast of the Sichuan basin, SW China. Chinese Journal of Geophysics (in Chinese), 60(2): 617-627. DOI:10.6038/cjg20170216
Li H Y, Shen Y, Huang Z X, et al. 2014. The distribution of the mid-to-lower crustal low-velocity zone beneath the northeastern Tibetan Plateau revealed from ambient noise tomography. Journal of Geophysical Research:Solid Earth, 119(3): 1954-1970. DOI:10.1002/2013JB010374
Li L G, Li B X. 2017. A discussion on the heat source mechanism and geothermal system of Gonghe-Guide basin and mountain geothermal field in Qinghai Province. Geophysical and Geochemical Exploration (in Chinese), 41(1): 29-34.
Li M K, Wu T F, Lin X, et al. 2019. Distribution of the crustal low velocity zones beneath the central and northeastern Tibetan Plateau:Insights from joint analysis of receiver functions and surface wave dispersion observations. Physics of the Earth and Planetary Interiors, 286: 179-189. DOI:10.1016/j.pepi.2018.12.003
Li X L, Wu G L, Lei Y D, et al. 2016. Suggestions for geothermal genetic mechanism and exploitation of Zhacang Temple geothermal energy in Guide county, Qinghai province. Journal of Jinlin University (Earth Science Edition) (in Chinese), 46(1): 220-229.
Li Y G. 2016. Hydrogeochemical characteristics and its origin analysis of Geothermal water in the Qia bu-qia area, Gonghe basin, Qinghai province[Master's thesis] (in Chinese). Nanchang: East China University of Technology.
Liao Z J, Zhang Z G. 1989. Geothermics of southern Tibet. Geological Review (in Chinese), 35(4): 366-373.
Lin W J, Gan H N, Wang G L, et al. 2016. Occurrence prospect of HDR and Target site selection study in southeastern of China. Acta Geologica Sinica (in Chinese), 90(8): 2043-2058.
Liu F, Lang X J, Lu C, et al. 2017. Thermophysical parameters and lithospheric thermal structure in Guide Basin, Northeast Qinghai-Tibet Plateau. Environmental Earth Sciences, 76(5): 1-12.
Ma F, Lin W J, Lang X J, et al. 2015. Deep geothermal structures of potential hot dry rock resources area in China. Geological Science and Technology Information (in Chinese), 34(6): 176-181.
Molnar P, Tapponnier P. 1975. Cenozoic tectonics of Asia:Effects of a continental collision. Science, 189(4201): 419-426. DOI:10.1126/science.189.4201.419
Royden L H, Burchfiel B C, King R W, et al. 1997. Surface deformation and lower crustal flow in Eastern tibet. Science, 276(5313): 788-790. DOI:10.1126/science.276.5313.788
Rutqvist J, Dobson P F, Garcia J, et al. 2015. The northwest geysers EGS demonstration project, california:pre-stimulation modeling and interpretation of the stimulation. Mathematical Geosciences, 47(1): 3-29. DOI:10.1007/s11004-013-9493-y
Rybach L. 1976. Radioactive heat production in rocks and its relation to other petrophysical parameters. Pure and Applied Geophysics, 114(2): 309-317. DOI:10.1007/BF00878955
Shen X J, Yang S Z, Shen J Y, et al. 1989. Experimental study of radiogenic heat production of granitic rocks in Tibet. Acta Petrologica Sinica (in Chinese), 5(4): 83-92.
Shi B Y, Zhang J T. 1982. Discussion on the Gap of Qin-Kun tectonic belt. Management & Strategy of Qinghai Land & Resources (in Chinese), (3): 21-29.
Stanley W D, Blakely R J. 1995. The Geysers-Clear Lake geothermal area, California-An updated geophysical perspective of heat sources. Geothermics, 24(2): 187-221. DOI:10.1016/0375-6505(94)00048-H
Sun Z X, Li B X, Wang Z L. 2011. Exploration of the possibility of hot dry rock occurring in the Qinghai Gonghe Basin. Hydrogeology and Engineering Geology (in Chinese), 38(2): 119-124, 129.
Tang J, Deng Q H, Zhao G Z, et al. 2001. Electric conductivity and magma chamber at the Tianchi volcano area in Changbaishan Mountain. Seismology and Geology (in Chinese), 23(2): 191-200.
Tian J, Pang Z H, Guo Q, et al. 2018. Geochemistry of geothermal fluids with implications on the sources of water and heat recharge to the Rekeng high-temperature geothermal system in the Eastern Himalayan Syntax. Geothermics, 74: 92-105. DOI:10.1016/j.geothermics.2018.02.006
Tian J, Pang Z H, Wang Y C, et al. 2019. Fluid geochemistry of the Cuopu high temperature geothermal system in the eastern Himalayan syntaxis with implication on its genesis. Applied Geochemistry, 110: 104422. DOI:10.1016/j.apgeochem.2019.104422
Turcotte D, Schubert G. 2014. Geodynamics. 3rd ed. Cambridge: Cambridge University Press.
Wang B, Li B X, Ma X H. 2015. The prediction of the depth and temperature of the reservoir in the evaluation of hot dry rock (HDR) for Gonghe-Guide basin. Ground Water (in Chinese), 37(3): 28-30, 72.
Wang J Y, Hu S B, Pang Z H, et al. 2012. Estimate of geothermal resources potential for hot dry rock in the continental Area of China. Science & Technology Review (in Chinese), 30(32): 25-31.
Wang J Y. 2015. Geothermics and Its Applications (in Chinese). Beijing: Science Press.
Wang Q, Chung S L, Li X H, et al. 2012. Crustal melting and flow beneath Northern Tibet:evidence from mid-miocene to quaternary strongly peraluminous rhyolites in the Southern Kunlun Range. Journal of Petrology, 53(12): 2523-2566. DOI:10.1093/petrology/egs058
Wang X, Chen L, Ai Y S, et al. 2018. Crustal structure and deformation beneath eastern and northeastern Tibet revealed by P-wave receiver functions. Earth and Planetary Science Letters, 497: 69-79. DOI:10.1016/j.epsl.2018.06.007
Wu Z B, Tao X, Badal J, et al. 2017. Crustal shear-wave velocity structure of northeastern Tibet revealed by ambient seismic noise and receiver functions. Gondwana Research, 41: 400-410. DOI:10.1016/j.gr.2015.08.009
Xu T F, Zhang Y J, Zeng Z F, et al. 2012. Technology Progress in an Enhanced Geothermal System (Hot Dry Rock). Science & Technology Review (in Chinese), 30(32): 42-45.
Xu T F, Yuan Y L, Jia X F, et al. 2018. Prospects of power generation from an enhanced geothermal system by water circulation through two horizontal wells:A case study in the Gonghe Basin, Qinghai Province, China. Energy, 148: 196-207. DOI:10.1016/j.energy.2018.01.135
Xue J Q, Gan B, Li B X, et al. 2013. Geological-geophysical characteristics of enhanced geothermal systems (Hot Dry Rocks) in Gonghe-Guide Basin. Geophysical and Geochemical Exploration (in Chinese), 37(1): 35-41.
Yan W D. 2015. Characteristics of Gonghe Basin hot dry rock and its utilization prospects. Science & Technology Review (in Chinese), 33(19): 54-57.
Yue G F, Deng X F, Xing L X, et al. 2015. Numerical simulation of hot dry rock exploitation using enhanced geothermal systems in Gonghe Basin. Science & Technology Review (in Chinese), 33(19): 62-67.
Zhang C, Jiang G Z, Shi Y Z, et al. 2018. Terrestrial heat flow and crustal thermal structure of the Gonghe-Guide area, northeastern Qinghai-Tibetan plateau. Geothermics, 72: 182-192. DOI:10.1016/j.geothermics.2017.11.011
Zhang C, Zhang S S, Li S T, et al. 2018. Geothermal characteristics of the Qiabuqia geothermal area in the Gonghe basin, northeastern Tibetan Plateau. Chinese Journal of Geophysics (in Chinese), 61(11): 4545-4557. DOI:10.6038/cjg2018L0747
Zhang C, Jiang G Z, Jia X F, et al. 2019. Parametric study of the production performance of an enhanced geothermal system:A case study at the Qiabuqia geothermal area, northeast Tibetan plateau. Renewable Energy, 132: 959-978. DOI:10.1016/j.renene.2018.08.061
Zhang G W, Guo A L, Yao A P. 2004. Western Qinling-Songpan continental tectonic node in China's continental tectonics. Earth Science Frontiers (in Chinese), 11(3): 23-32.
Zhang H F, Chen Y L, Xu W C, et al. 2006. Granitoids around Gonghe basin in Qinghai province:Petrogenesis and tectonic implications. Acta Petrologica Sinica (in Chinese), 22(12): 2910-2922.
Zhang J, Li W Y, Tang X C, et al. 2017. Geothermal data analysis at the high-temperature hydrothermal area in Western Sichuan. Science China Earth Sciences, 60(8): 1507-1521. DOI:10.1007/s11430-016-9053-2
Zhang L Y, Liu Q Y, He L J. 2016. The different lithospheric thermal structure of North China Craton and its implications. Chinese Journal of Geophysics (in Chinese), 59(10): 3618-3626. DOI:10.6038/cjg20161009
Zhang P Z, Shen Z K, Wang M, et al. 2004. Continuous deformation of the Tibetan Plateau from global positioning system data. Geology, 32(9): 809-812. DOI:10.1130/G20554.1
Zhang S Q, Yan W D, Li D P, et al. 2018a. Characteristics of geothermal geology of the Qiabuqia HDR in Gonghe Basin, Qinghai Province. Geology in China (in Chinese), 45(6): 1087-1102.
Zhang S Q, Jia X F, Li B X, et al. 2018b. Characteristics of Wudalianchi volcanic edifice and hot dry rock geophysical field in Heilongjiang province. Geophysical & Geochemical Exploration (in Chinese), 42(3): 473-480.
Zhang S S, Zhang L, Tian C C, et al. 2019. Occurrence geological characteristics and development potential of hot dry rocks in Qinghai Gonghe Basin. Journal of Geomechanics (in Chinese), 25(4): 501-508.
Zhang Y. 2016. The formation mechanism and development potential of Hot Dry Rock: A Case study of Songliao Basin[Master's thesis] (in Chinese). Xi'an: Chang'an University.
Zhao Z, Chen H J, Ma J Q, et al. 2009. On terrestrial heat resource assessment and its exploitation and utilization in Qiabuqia area, Gonghe Basin, Qinghai. Journal of Qinghai Environment (in Chinese), 19(2): 81-84.
Zheng D, Li H Y, Shen Y, et al. 2016. Crustal and upper mantle structure beneath the northeastern Tibetan Plateau from joint analysis of receiver functions and Rayleigh wave dispersions. Geophysical Journal International, 204(1): 583-590. DOI:10.1093/gji/ggv469
Zhuang Y Q, Zhang C H, Zhu M C, et al. 2017. Laboratory study of improving granite permeability by using mud acid as the chemical stimulus. Safety and Environmental Engineering (in Chinese), 24(2): 16-21, 45.
白登海, 廖志杰, 赵国泽, 等. 1994. 从MT探测结果推论腾冲热海热田的岩浆热源. 科学通报, 39(4): 344-347.
多吉. 2003. 典型高温地热系统-羊八井热田基本特征. 中国工程科学, 5(1): 42-47. DOI:10.3969/j.issn.1009-1742.2003.01.008
管彦武, 管烨, 高瑞, 等. 2012. 根据地震纵波速度分析青藏高原地壳放射性生热率和地幔热状态. 吉林大学学报(地球科学版), 42(2): 562-568, 574.
郭进京, 周安朝. 2008. 地热异常区高温岩体圈定的思路与方法. 地热能, (4): 4-7.
郭晓玉, 高锐, 高建荣, 等. 2017. 综合数据分析青藏高原东北缘六盘山地区构造形变及其构造成因独特性探讨. 地球物理学报, 60(6): 2058-2067. DOI:10.6038/cjg20170603
胡久常, 白登海, 王薇华, 等. 2007. 雷琼火山区地下深部大地电磁探测与电性结构分析. 华南地震, 27(1): 1-7. DOI:10.3969/j.issn.1001-8662.2007.01.001
姜光政, 高堋, 饶松, 等. 2016. 中国大陆地区大地热流数据汇编(第四版). 地球物理学报, 59(8): 2892-2910. DOI:10.6038/cjg20160815
康维海. 2016. 青海贵德发现地表水温达96.6℃优质温泉. 地热能, (1): 17.
郎旭娟, 刘峰, 刘志明, 等. 2016. 青海省贵德盆地大地热流研究. 地质科技情报, 35(3): 227-232.
李春荣, 饶松, 胡圣标, 等. 2017. 川东南焦石坝页岩气区现今地温场特征. 地球物理学报, 60(2): 617-627. DOI:10.6038/cjg20170216
李林果, 李百祥. 2017. 从青海共和-贵德盆地与山地地温场特征探讨热源机制和地热系统. 物探与化探, 41(1): 29-34.
李小林, 吴国禄, 雷玉德, 等. 2016. 青海省贵德扎仓寺地热成因机理及开发利用建议. 吉林大学学报(地球科学版), 46(1): 220-229.
李永革. 2016.青海省共和盆地恰卜恰地区地下热水水文地球化学特征及成因分析[硕士论文].南昌: 东华理工大学.
廖志杰, 张振国. 1989. 西藏南部地热. 地质论评, 35(4): 366-373. DOI:10.3321/j.issn:0371-5736.1989.04.011
蔺文静, 甘浩男, 王贵玲, 等. 2016. 我国东南沿海干热岩赋存前景及与靶区选址研究. 地质学报, 90(8): 2043-2058. DOI:10.3969/j.issn.0001-5717.2016.08.031
马峰, 蔺文静, 郎旭娟, 等. 2015. 我国干热岩资源潜力区深部热结构. 地质科技情报, 34(6): 176-181.
沈显杰, 杨淑贞, 沈继英, 等. 1989. 西藏岩浆岩放射性生热率的实验研究. 岩石学报, 5(4): 83-92. DOI:10.3321/j.issn:1000-0569.1989.04.009
石宝颐, 张峻太. 1982. 论秦-昆构造带的共和"缺口". 青海国土经略, (3): 21-29.
孙知新, 李百祥, 王志林. 2011. 青海共和盆地存在干热岩可能性探讨. 水文地质工程地质, 38(2): 119-124, 129. DOI:10.3969/j.issn.1000-3665.2011.02.021
汤吉, 邓前辉, 赵国泽, 等. 2001. 长白山天池火山区电性结构和岩浆系统. 地震地质, 23(2): 191-200. DOI:10.3969/j.issn.0253-4967.2001.02.008
王斌, 李百祥, 马新华. 2015. 青海共和-贵德干热岩勘查评价中热储温度与深度预测. 地下水, 37(3): 28-30, 72. DOI:10.3969/j.issn.1004-1184.2015.03.011
汪集旸, 胡圣标, 庞忠和, 等. 2012. 中国大陆干热岩地热资源潜力评估. 科技导报, 30(32): 25-31. DOI:10.3981/j.issn.1000-7857.2012.32.002
汪集旸. 2015. 地热学及其应用. 北京: 科学出版社.
许天福, 张延军, 曾昭发, 等. 2012. 增强型地热系统(干热岩)开发技术进展. 科技导报, 30(32): 42-45. DOI:10.3981/j.issn.1000-7857.2012.32.004
薛建球, 甘斌, 李百祥, 等. 2013. 青海共和-贵德盆地增强型地热系统(干热岩)地质-地球物理特征. 物探与化探, 37(1): 35-41.
严维德. 2015. 共和盆地干热岩特征及利用前景. 科技导报, 33(19): 54-57. DOI:10.3981/j.issn.1000-7857.2015.19.008
岳高凡, 邓晓飞, 邢林啸, 等. 2015. 共和盆地增强型地热系统开采过程数值模拟. 科技导报, 33(19): 62-67. DOI:10.3981/j.issn.1000-7857.2015.19.010
张超, 张盛生, 李胜涛, 等. 2018. 共和盆地恰卜恰地热区现今地热特征. 地球物理学报, 61(11): 4545-4557. DOI:10.6038/cjg2018L0747
张国伟, 郭安林, 姚安平. 2004. 中国大陆构造中的西秦岭-松潘大陆构造结. 地学前缘, 11(3): 23-32. DOI:10.3321/j.issn:1005-2321.2004.03.004
张宏飞, 陈岳龙, 徐旺春, 等. 2006. 青海共和盆地周缘印支期花岗岩类的成因及其构造意义. 岩石学报, 22(12): 2910-2922.
张健, 李午阳, 唐显春, 等. 2017. 川西高温水热活动区的地热学分析. 中国科学:地球科学, 47(8): 899-915.
张林友, 刘琼颖, 何丽娟. 2016. 华北克拉通热结构差异性特征及其意义. 地球物理学报, 59(10): 3618-3626. DOI:10.6038/cjg20161009
张森琦, 严维德, 黎敦朋, 等. 2018a. 青海省共和县恰卜恰干热岩体地热地质特征. 中国地质, 45(06): 1087-1102.
张森琦, 贾小丰, 李百祥, 等. 2018b. 黑龙江省五大连池火山机构与干热岩地球物理场特征. 物探与化探, 42(3): 473-480.
张盛生, 张磊, 田成成, 等. 2019. 青海共和盆地干热岩赋存地质特征及开发潜力. 地质力学学报, 25(4): 501-508.
张杨. 2016.干热岩形成机理及开发潜力研究——以松辽盆地为例[硕士论文].西安: 长安大学.
赵振, 陈惠娟, 马建青, 等. 2009. 青海省共和盆地恰卜恰地区地热资源评价与开发利用. 青海环境, 19(2): 81-84. DOI:10.3969/j.issn.1007-2454.2009.02.010
庄亚芹, 张灿海, 朱明成, 等. 2017. 土酸化学刺激剂提高花岗岩渗透率的试验研究. 安全与环境工程, 24(2): 16-21, 45.