文章快速检索    
 
  中国地质 2019, Vol. 46 Issue (4): 927-936  
0
引用本文
付光明, 黄进调, 刘阳, 李晓兵, 陈国玉, 黄艺. 2019. 高密度电阻率法与CSAMT法在江西会昌县坝背地热勘查中的综合探测[J]. 中国地质, 46(4): 927-936.  
Fu Guangming, Huang Jintiao, Liu Yang, Li Xiaobing, Chen Guoyu, Huang Yi. 2019. Multi-electrode resistivity method and CSAMT method in geothermal exploration of Babei area in Huichang County, Jiangxi Province[J]. Geology in China, 46(4): 927-936. (in Chinese with English abstract).  

高密度电阻率法与CSAMT法在江西会昌县坝背地热勘查中的综合探测
付光明1,2,3, 黄进调4, 刘阳4, 李晓兵4, 陈国玉4, 黄艺4    
1. 东华理工大学地球物理与测控技术学院, 江西 南昌 330013;
2. 中国地质科学院, 北京 100037;
3. 中国地质调查局-中国地质科学院地球深部探测中心, 北京 100037;
4. 江西省煤田地质局普查综合大队, 江西 南昌 330001
摘要: 地热资源是一种绿色低碳、极具竞争力的可再生能源。江西省会昌县坝背地区断裂构造密集,地热异常广泛发育。为了查明勘查区各断裂构造的产状、规模及深部变化特征,为寻找地热水钻孔布置提供依据,项目组首先在已知温泉位置布置3条高密度试验测线,接着在已知温泉位置以南布置4条高密度测线,然后在高密度电阻率法资料初步解释成果的基础上再施工可控源音频大地电磁测深法(CSAMT),得出其中两条深大断裂是温泉主要的导水导热构造、断裂深切至基底界面、基底界面起伏明显的结论。同时圈定了4个低阻异常区,其中低阻异常范围最大的区域长约1000 m,宽约240 m,高约200 m,是含水最为有利区域,亦是寻找地热水的直接依据,建议在其正上方布设钻孔进行验证,孔深约800 m,满足深大断裂和基底对低阻区的深度控制。
关键词: 坝背    地热水    高密度电阻率法    可控源音频大地电磁测深    
中图分类号:P622+.2;P314.1            文献标志码:A             文章编号:1000-3657(2019)04-0927-10
Multi-electrode resistivity method and CSAMT method in geothermal exploration of Babei area in Huichang County, Jiangxi Province
FU Guangming1,2,3, HUANG Jintiao4, LIU Yang4, LI Xiaobing4, CHEN Guoyu4, HUANG Yi4    
1. School of Geophysics and Measurement-control Technology, East China Institute of Technology, Nanchang, 330013, Jiangxi, China;
2. Chinese Academy of Geological Sciences, Beijing, 100037, China;
3. China Deep Exploration Center, China Geological Survey & Chinese Academy of Geological Sciences, Beijing, 100037, China;
4. Jiangxi Coal Geology Bureau Comprehensive Survey Party, Nanchang 330001, Jiangxi, China
Abstract: Geothermal resources are a kind of green, low carbon and highly competitive renewable energy. The Babei area of Huichang County in Jiangxi Province is characterized by dense faults. In order to find out the characteristics of the attitude, scale and deep change of the fault structure in the exploration area so as to provide a basis for finding the layout of geothermal water drilling, the authors arranged three multi-electrode resistivity method test lines at the location of the known hot springs, and determined the validity of the geophysical method and the water guide fault in the survey area according to the principle of the shallow to deep exploration from the known to the unknown. On the basis of the preliminary interpretation results of multielectrode resistivity method, the controllable source audio magnetotelluric sounding method (CSAMT) was constructed according to the preliminary interpretation of multi-electrode resistivity method. It is concluded that two of the deep faults are the main thermal conductive structures in hot springs, that the faults cut deeply into the basement interface, and that the basement interface undulates obviously. At the same time, four low resistivity anomalous zones were delineated, of which the largest one is about 1000m in length, 240m in width and 200 m in height. It is the most advantageous water-bearing area and the direct basis in the search for geothermal water. It is suggested that drilling holes should be laid directly above the anomalous zones for verification. The hole depth is about 800m, which meets the depth control of deep and large faults and the base for low resistivity zones.
Key words: Babei    geothermal    multi-electrode resistivity method    CSAMT    

1 引言

能源与环境问题一直是困扰世界各国可持续发展的重大问题,进入21世纪,节能减排、开发新能源和走可持续发展道路受到越来越多的关注(刘凯等,2017)。地球是个巨大的热库,其内核(地核)的温度高达约6000℃(汪集旸等,2012)。地核与地表巨大的温差使得地球在不断地向外(大气层)散发着热量,同时地壳内部放射性衰变热、势能转换热、摩擦热等也在不断生成与供给(Clauser,2005)。而在当前的技术经济和地质环境条件下,可以从地壳表面下一定深度内科学、合理地开发出来的地下岩石中的地热能量和地热流体中的热能量及其伴生的有用组分被称之为地热资源(孔维臻,2013)。地热资源是一种极具竞争力的清洁可再生能源,很大程度上可以缓解因大量使用化石能源所造成严重的空气污染问题(马伟斌等,2016),按其成因和产出条件可分为水热型和干热岩型,利用方式分为直接利用和发电两种。

断裂构造对地热田分布具有直接的控制作用(苗可等,2012),因为它是地热水主要储存和运移通道,要掌握地热水的赋存特征,必须掌握断裂的产状及展布特征。江西省会昌县坝背地区地处武夷山和南岭余脉地带,受太平洋构造域的制约,晚侏罗世发生大规模岩浆侵入和火山喷发(燕山期),早白垩世晚期以后发生了强烈的伸展作用,强烈的地壳运动形成密集的断裂构造,成为地下热源导热构造带,并在断裂带附近形成地热异常区。区内有分布于会昌盆地东侧大富足岩体内的河草坑铀矿田和桂坑岩体中的铀矿化(蔡煜琦等,1997),说明区内具有放射性衰变热的持续供给来源。

地下的热作用能明显改变岩石的地球物理性质,构造带、地热田的生储盖等不同部位均具有较明显的电性差异(Wright et al., 1985),地热区电阻率的空间分布不仅受围岩性质的影响,还与勘探对象——热水的分布直接有关(切列缅斯基,1982)。总体上,随着温度的升高,电阻率在不断的降低,可见在地热田勘查中,低阻体成为最重要的寻找目标,由于它的直接指示作用和地热田不同部位电性差异明显,电法和电磁法也就成为地热资源勘探中最直接有效的方法技术(曾昭发等,2012)。直流电测深法、高密度电法、激发极化法、瞬变电磁法等都不乏成功的范例,而CSAMT更是因为其抗干扰能力强、探测深度大、横向分辨率高、受高阻层屏蔽作用影响小等优点,越来越多地被应用到地热勘查中(周仕新,2013)。为了查明会昌县坝背勘查区各断裂构造的产状、规模及深部变化特征,为寻找地热水钻孔布置提供依据,遵循从已知到未知,由浅至深的勘探原则,项目组首先在已知温泉位置布置3条高密度试验测线,确定地球物理方法在该测区的有效性和导水断层的形态、主要物性参数等,接着在已知温泉位置以南布置4条高密度测线,然后在高密度电阻率法资料初步解释成果的基础上再施工CSAMT。

2 地质概况

会昌县坝背区域上地处武夷山和南岭余脉地带,属低山丘陵地区,地势由西南向东北倾斜,境内峰峦起伏。勘查区最高点在南部海拔标高约389.9 m,最低点海拔标高约190 m,最大相对高差近200 m。工作区及周边地球物理工作程度低,未开展过系统的地球物理工作。区域地质概况详细论述如下(图 1)。

图 1 坝背地区区域地质、勘探范围及测线布置图 Fig. 1 Regional geology, exploration limits and survey line layout of Babei area
2.1 地层

勘查区内地层岩性特征由新到老分别为第四系(Qh1),上白垩统周田组(K2z)、茅店组(K2m2),下白垩统鸡笼嶂组(K1j)、震旦系桃溪岩组(Pt32 - 3 tx)。其中Qh1主要分布在河水两岸及山间沟谷附近,成因类型为冲积作用,岩性具有二元结构,上部为黄褐色、褐色亚黏土、亚砂土,厚度1.0~2.5m;下部为砂砾石层,砾石成分为石英砂岩,粒径一般2~4 cm,大者达10 cm,结构松散,厚度1.5~3.5 m,总厚度2.5~ 6.0 m。K2z大面积出露勘查区内,为紫红色中厚层状钙质粉砂岩夹薄层细砂岩,下部为紫红色细粒石英砂岩,厚度170~450 m。K2m2大面积出露勘查区西侧,为紫红色巨厚层状砾岩夹钙质粉砂岩、泥岩及粗砂岩,砾石成分以晶屑凝灰岩为主,厚度300~ 800 m。K1j出露于勘查区东北侧,为杂色流纹质熔结凝灰岩,厚度大于200 m。Pt32-3tx小面积出露勘查区北东角,为勘查区褶皱基底,岩性主要为杂色黑云斜长变粒岩、石英片岩,厚度大于400m。这些地层的成分与结构特征无疑为地热水的存储提供了良好条件。

2.2 岩浆岩

勘查区出露岩浆岩来自两个不同时期,出露面积都较大。一是分布于勘查区北东角印支期晚三叠世罗珊序列(T3γt)黑云母、二云母花岗岩;二是分布于勘查区南东角加里东期早志留世汤湖序列(S1gnγ-ξγ)二长、正长花岗岩。

2.3 地质构造

本区的区域大地构造单元在加里东期为中南武夷山褶皱带(江西省区域地质志,1984),之后受太平洋构造域的制约,晚侏罗世发生大规模岩浆侵入和火山喷发(燕山期),在武夷山隆起地带形成一系列的NE、NNE向火山岩盆地,到早白垩世晚期以后,受太平洋板块北北西向左行走滑影响,发生了强烈的伸展作用,盆地由坳陷转为断陷,沉积了巨厚的红色碎屑岩系,不整合地叠覆于火山岩盆地之上。

2.4 地球物理特征

周田组为紫红色中厚层状钙质粉砂岩夹薄层细砂岩,紫红色细粒石英砂岩其视电阻率小于200 Ω·m。茅店组为紫红色巨厚层状砾岩夹钙质粉砂岩、泥岩及粗砂岩,砾石成分以晶屑凝灰岩为主其视电阻率为300~400 Ω·m。鸡笼嶂组为杂色流纹质熔结凝灰岩其视电阻率400~600 Ω·m。震旦系桃溪岩组,为勘查区褶皱基底,岩性主要为杂色黑云斜长变粒岩、石英片岩,其视电阻率大于800 Ω·m。断裂破碎带多系充水,其电阻率相对围岩较低,在视电阻率断面图上多呈现低阻特征,异常等值线密集或扭曲。这就为电法探测断裂破碎带提供了物性条件。

3 数据采集及处理 3.1 测网布置

勘探区地球物理测线的布置安排在地表调查后,根据地表地质调查情况来确定地球物理测线的起止点。因勘探区内及相邻勘探区均没开展过地球物理工作,进一步了解地球物理方法在该测区的有效性和导水断层的形态、主要物性参数等,为勘查区地球物理参数的设置和地球物理资料的解释工作提供依据,在勘探区外已知的车心出露点布设一条高密度电阻率测量线,编号为S1线,测线长度为0.3 km。由于车心热水泉出露点走向长约50 m,方位为北东10°左右,S1线施工完后视电阻率断面图反映不够明显,因此在S1线南约240 m处布置了长750 m的S2线,方位角为69°;在与车心热水泉出露点走向斜交处布置了长590 m的S3线,方位角为344°,同时将S1线延长到750 m,受地形限制其方位角调整为93°,3条测线极距均为10 m。然后根据地表地质调查情况在勘探区内共布置4条高密度测线,线距500 m,极距10 m。由北向南分别为G1线、G2线、G3线和G4线,每条测线长约1.65 km,方位角为102°。G1测线东侧遇有陡坎、悬崖等障碍物影响布线,因此,实际布线时G1测线往南偏移150 m左右(图 1)。

在高密度电阻率法资料的初步解释成果的基础上,确定在勘探范围内重点勘探区开展可控源音频大地电磁测深工作,由于在勘探区中部有一由西南-北东走向的220 kV的高压电线,对可控源音频大地电磁测深的野外测量数据采集数据质量会有一定的影响,故勘探区内根据实际情况布置了5条可控源音频大地电磁测深测线,测线由南向北、由西向东分别为K1线、K2线、K3线、K4线、K5线。线距在300~500 m,长度在300~840 m,点距30 m,方位角为92°。K2线、K4线、K5线基本上与高密度电阻率法测线重合(图 1)。

3.2 数据采集

高密度电阻率法使用的仪器是重庆奔腾数控技术研究所生产的WDJD-3型和WGMD-1型高密度电阻率测量系统。在试验线S1线布置了两个排列同时用高密度电阻率法α1排列装置和α2排列装置进行测量,经对资料的分析认为采用高密度电阻率法α1排列装置进行测量效果更佳,因而决定高密度电阻率测量工作采用α1排列装置(图 2),电极距10 m,60根电极,排列长度590 m。

图 2 α装置固定断面扫描示意图 Fig. 2 Scanning sketch of fixed section of α device

CSAMT投入的仪器设备为加拿大凤凰公司的V8网络化多功能电法仪,为保证工作质量,在测量之前标定了仪器和磁棒,并进行了一致性测试,测试结果均符合规范要求。CSAMT测点点距30 m,频率范围1~9600 Hz,有效探测深度大于2.0 km。采集方式为标量模式,接收端和发射端工作示意图如图 3。本项目布设了2个场源,一号场源极距为1316 m,最小收发距为9810 m,在第一天工作结束后,发现频率在100 Hz左右进入近场并且数据信躁比大,于是对发射场源进行了调整,并通过增大收发距,加大供电极距,增大供电电流来保证数据质量。最终选定二号场源位置,极距为2079 m,最小收发距为13269 m,最大供电电流17A。

图 3 CSAMT法V8(一主一辅)网络化野外作业方式示意图 Fig. 3 Schematic diagram of field operation mode of CSAMT method using V8 (one main and one auxiliary) network

本次可控源音频大地电磁法测量水平方向电场(MN)平行于场源(AB),水平磁场垂直于场源布设。共设计55个频点,工作频带为1~9600 Hz(表 1),频点间隔均匀分布,高、中频段适度加密。

表 1 CSAMT工作频率 Table 1 CSAMT working frequency list
3.3 数据质量评价

本项目共获得高密度电阻率法测量坐标测点261个,质量检查点21个,质量检查点占测点总数7.4%,高密度电阻率法测点数据没有废点,检查点相对均方误差均小于5%,数据质量良好;可控源音频大地电磁测点82个,质量检查点5个,质量检查点占测点总数5.7%,不存在三类点,检查点相对均方误差均小7%,数据质量良好。

3.4 数据处理

高密度电法野外采集的数据传入计算机后,对每个排列所测数据经软件(RES2DINV)进行必要的一些编辑,对个别畸变点进行剔除(图 4)。由于测线地形相差大,因此在进行二维反演时必须进行地形改正,本次计算选用Schwarz-Christoffel变换法作地形改正计算。然后使用有限差分法进行反演得到每个排列的视电阻率断面图,根据反演结果所得的视电阻率断面图,进行破碎带等地质体的初步划分。然后对每个排列所测的数据进行数据拼接得到整条测线的数据,经一系列处理后,得到整条测线的视电阻率断面图。

图 4 G1线1排列数据圆滑处理前(a)后(b)对比图 Fig. 4 The contrast diagram of Permutation1, Line G1 data (a) and that processed by smoothing (b)

CSAMT数据处理以测线为单位进行,数据预处理使用的软件为加拿大凤凰公司配套软件CMT Pro Version,通过加载数据;核对点位、极距等信息;查看数据质量来完成预处理过程,最后输出*.AVG文件用于反演。数据反演软件为CSAMT-SW,过程包括:测线记录型文件-(D File)、点位偏差校正*-J*、剔除跳点*-ED*、曲线平滑*-S2*图 5)、删除近场频点*-WR*、删除坏道数据(重复点数据)*- WC*、静态位移校正、Bostick反演生成*.DMT文件等,最终绘制成图。

图 5 未经平滑处理(a)、(c)和处理后(b)、(d)视电阻率和相位曲线对比图 Fig. 5 The contrast diagram of apparent resistivity and phase curve (a), (c) and that processed by smoothing (b), (d)
4 电阻率模型分析

高密度视电阻率断面图是以实测数据经数据处理后按点(线)排列方式形成的地表以下不同深度的视电阻率值绘制而成。纵轴为标高,横轴为测点对应距离,剖面方向为102°(图 6图 7)。通过已收集的钻孔资料可知地下热水是通过岩层裂隙涌向地面的,高密度电阻率法试验S3线就通过了该裂隙带,它在视电阻率断面图上反映为低阻区,据此对全区各线进行解释。鉴于高密度电阻率法的探测深度低于CSAMT,由高密度电阻率法推断的断裂或破碎带在深部的发育和延伸情况还需要CSAMT做进一步的追踪。

图 6 高密度试验测线视电阻率断面图 Fig. 6 Section of high density test line apparent resistivity
图 7 高密度测线电阻率断面图 Fig. 7 Section of high density line apparent resistivity

CSAMT的主要参数为卡尼亚视电阻率,由实测正交电、磁信号振幅计算而来,与常规电法中视电阻率的意义相同,主要反映勘探体积内岩(矿)石的综合导电性能。图 8为实测数据经处理后按点(线)排列方式形成的地表以下不同深度的卡尼亚视电阻率值绘制而成。纵轴为标高,横轴为测点对应距离,剖面方向为92°。各勘探线CSAMT二维反演电阻率断面图解译如下。

图 8 可控源测线视电阻率断面 Fig. 8 Section of apparent resistivity of CSAMT line

由K1线CSAMT二维反演电阻率断面图(图 8)可看出,剖面内推断的基底起伏大,基底以上地层呈不规则分布。剖面0~480 m,标高-250 m以上;剖面750~840 m,标高240~-50 m为高阻区,推断为白垩纪晚期茅店组砂砾岩。高阻体以下至基底界面区域存在明显的低阻异常带,剖面90~360 m,标高在-396~-550 m、埋深在682~814 m视电阻率值异常低,为低阻异常区,编为1号低阻异常区。剖面660~840 m,标高在-82~-176 m,埋深在396~491 m视电阻率值异常偏小,为低阻异常区,编为2号低阻异常区。剖面120 m附近,标高在230~150 m处视电阻值与两侧差异明显,推测此处有一断裂构造,结合地质资料和G1线高密度视电阻率断面图将该断裂命名为F3,该断裂构造将0~480 m的高阻错断。剖面480 m处两侧电阻率值差异显著,推断该位置有一断裂,结合地质资料和G1线高密度视电阻率断面图(图 7)将该断裂命名为F1,此构造深切至基底。剖面750 m处电阻率差异大且地层错动,推断该处有一断裂构造,此构造深切基底,且造成切割部位发生一定的凹陷,两侧则抬升的现象,结合地质资料和G1线高密度视电阻率断面图将该断裂命名为F4。

K2剖面内推断的基底起伏较小,基底以上地层呈不规则分布。剖面0~480 m处,标高在250~130 m和剖面0~270 m处,标高在-50~-170 m为高阻区,这2个高阻区认为系同一电性层被断裂F3错断所致,推断其高阻区为白垩纪晚期茅店组砂砾岩。剖面180~360 m处,标高在-200~-462 m,埋深在490~758 m视电阻率值异常偏小,存在一明显的低阻异常区,为1号低阻异常区往南的延伸。剖面390 m处电阻率两侧的差异大,延伸深,基底略微凹陷,推断存在一断裂构造,结合G1线高密度视电阻率断面图推断为F4断层往南的延伸。K3剖面内推断的基底起伏大,基底以上地层呈不规则分布;剖面0~240 m处,标高在-292~-584 m,埋深在532~ 826 m视电阻率值异常小,为低阻异常区,同样为1号低阻异常区在该处的表现。剖面120~180m处,标高在160~-220 m为高阻区,推断为白垩系晚期茅店组砂砾岩,该砂砾岩近直立,底部被F3断裂错断。剖面270~300 m处两侧的电阻率差异大,推断该位置存在一断裂构造,构造近乎直立,结合地质资料和G2线高密度视电阻率断面图推测为F4在该处的表现,基底凹陷明显。

K4剖面内推断的基底起伏明显,剖面中部位置基底隆起。剖面0~180 m处,标高120~-800 m和240~300 m处,标高50~-50 m为高阻区,这2个高阻区有错动,推断存在一断裂构造,结合地质资料和G1线高密度视电阻率断面图将该断裂命名为F5。K5剖面内推断的基底起伏平缓,基底以上地层分层明显,各种岩性在横向上分布较均匀,地表至标高-500 m范围内存在三层次高阻,其间夹杂2层厚度50~100 m的低阻区,推断此区域为白垩系茅店组砂砾岩,剖面上部的低阻区编为3号低阻异常区,3号低阻异常区顶、底标高186~131 m,埋深在95~156 m。在剖面中部的低阻区编为4号低阻异常区,4号低阻异常区顶、底标高-82~-316 m,埋深在368~600 m。剖面90~150 m处从表层到基底电性层均有错动,推断该位置存在一断裂构造,结合地质资料和G2线高密度视电阻率断面图将该断裂命名为F5。

5 解释与讨论

区内断裂主要由一系列北北东或北东及近东西或北西向规模不等的逆冲、斜冲断层和推滑-推覆断裂构造组成,这些密集的(断裂)成为地下热源导热构造带,并在断裂带附近形成地热导常区。通过对电阻率模型的解释获得了测区高阻基底埋深在750~1250 m,基底起伏较大;同时确定了大断层F1、F3、F4、F5的性质和产状(图 8)。F1正断层分布在测区北部边缘,走向为弧形,总体为北西向,倾向为西—西南,倾角73°,断层深切至基底,为测区内导热构造。F3逆断层分布在测区西部,走向为北东向,倾向南东,倾角55°~65°,测区内走向长1800 m。F4正断层分布在测区中东部,走向在G4线— K3线为近南北向,在K3线—K1线为北东向,倾向西—北西,倾角86°,测区内走向长1800 m,深切至基底,为测区内导热构造。F5逆断层分布在测区东部,走向在G4线—G3线为近南北向,在G3线~G1线为北东向,倾向由东转南东,倾角82°,测区内走向长1300 m,深切至基底面。

本区以侵蚀构造低山地形为主,地势上东高西低,植被繁茂,构造裂隙水发育。地下水主要赋存于断裂破碎带和构造裂隙中,富水性极不均一。地下水总体运动方向由南东向北西径流,一部分在车心村的西部(温泉地)排泄于地表,另一部分沿深部构造径流出勘查区。

联合K1、K2、K3剖面图中来看,1号异常区基本上分布在以K1(270 m处)—K2(240 m处)—K3(180 m处)连线为轴线的左右两侧各约100 m范围内,东部为北东向的高阻条带,为F4深大断裂构造的反映。1号异常区在4个低阻异常区内范围最大,长约1000 m,宽约240 m,高约200 m,由南往北的走向特征为北东转向北西,是地热水赋存最为可能的区域,建议在K1线360 m处布设钻孔进行验证,孔深约800 m,满足深大断裂和基底对低阻区的深度控制。

6 结论

(1)通过高密度和CSAMT电阻率勘测推断出F4、F1构造深切至基底界面,为深大断裂,是温泉主要的导水导热构造。

(2)依据工作区的电性结构特点,厘定了基底的岩性界面,刻画了基底面的起伏情况。

(3)根据CSAMT二维反演的断面图上的低阻异常,圈定了4个低阻异常区,其中1号异常范围最大,为本区含水有利区,亦是寻找地热水的直接依据,建议在K1线360 m处布设钻孔进行验证,孔深约800 m。

(4)验证了高密度和CSAMT电阻率法在坝北区域地热探测中的有效性,为类似地区地热勘查提供了一定的借鉴。

致谢: 本文在写作过程中得到了严加永教授级高工和张昆副研究员的宝贵建议,以及编辑部的大力支持,在此表示衷心的感谢!

参考文献
Cai Yuqi, Chen Zuyi, Shi Zuhai, Feng Quanhong, Chen Yuehui. 1997. Huichang shallow thermo-upwelling extensional tectonnics and uranium metallogenesis[J]. Uranium Geology, (5): 257-263 (in Chinese with English abstract).
Jiangxi Bureau of Geology and Mineral Resources. 1984. Regional Geological Records of Jiangxi Province[M]. Beijing: Geological Publishing House, 1-921 (in Chinese).
Kong Weizhen. 2013. Study on Economic Evaluation of Geothermal Energy Developing and Utilizing[D]. China University of Geosciences(in Chinese with English abstract).
Lei Xiaodong, Guan Wei, Guo Gaoxuan, Zhu Liqiong. 2014. Application of CSAMT in the survey of karst water in the eastern Yanqing basin of Beijing[J]. Geotechnical Investigation & Surveying, (5): 89-98 (in Chinese with English abstract).
Liu Kai, Wang Shanshan, Sun Ying, Cui Wenjun, Zhu Deli. 2017. Characteristics and regionalization of geothermal resources in Beijing[J]. Geology in China, 44(6): 1128-1139 (in Chinese with English abstract).
Ma Weibin, Gong Yulie, Zhao Daiqing, Xu Qionghui, Qing Hanshi, Chen Yong. 2016. Geothermal energy exploitation utilization and its development trend in China[J]. Bulletin of Chinese Academy of Sciences, 31(2): 199-207 (in Chinese with English abstract).
Miao Ke, Cao Xinling. 2012. Study of geothermal method for comprehensive detection of earthquake and CSAMT method[J]. Petroleum Instrumenis, 26(6): 60-62, 104 (in Chinese with English abstract).
Wang Jiyang, Hu Shengbiao, Pang Zhonghe, He Lijuan, Zhao Ping, Zhu Chuanqing, Rao Song, Tang Xiaoyin, Kong Yanlong, Luo Lu, Li Weiwei. 2012. Estimate of geothermal resources potential for hot dry rock in the continental area of China[J]. Science & Technology Review, 30(32): 25-31 (in Chinese with English abstract).
Wright P M, Ward S H, Ross H P. 1985. State-of-the art Geophysical exploration for geothermal resources[J]. Geophysic, 50(12): 2666-2699. DOI:10.1190/1.1441889
Zeng Zhaofa, Chen Xiong, Li Jing, Li Tonglin, Zhang Lianghuan. 2012. Advancement of geothermal geophysics exploration[J]. Progress in Geophysics, 27(1): 168-178 (in Chinese with English abstract).
Zhou Shixin. 2013. The application of CSAMT in the geothermal surveying[J]. Technological development of enterprise, 32(13): 55-57 (in Chinese with English abstract).
Череменский Г А. 1982. Practical Geothermal Science[M]. Beijing: Geological Publishing House.
蔡煜琦, 陈祖伊, 施祖海, 封全宏, 陈跃辉. 1997. 会昌浅层次热隆伸展构造与铀成矿[J]. 铀矿地质, (5): 257-263.
江西省地质矿产局. 1984. 江西省区域地质志[M]. 北京: 地质出版社, 1-921.
孔维臻. 2013.地热资源开发利用经济评价研究[D].中国地质大学. http://www.cnki.com.cn/Article/CJFDTotal-CSDI201522040.htm
雷晓东, 关伟, 郭高轩, 朱丽琼. 2014. 可控源音频大地电磁测深在北京延庆盆地东部岩溶水勘查中的应用[J]. 工程勘察, (5): 89-98.
刘凯, 王珊珊, 孙颖, 崔文君, 朱德莉. 2017. 北京地区地热资源特征与区划研究[J]. 中国地质, 44(6): 1128-1139.
马伟斌, 龚宇烈, 赵黛青, 徐琼辉, 秦汉时, 陈勇. 2016. 我国地热能开发利用现状与发展[J]. 中国科学院院刊, 31(2): 199-207.
苗可, 曹新领. 2012. 地震与CSAMT法综合探测地热方法研究[J]. 石油仪器, 26(6): 60-62, 104. DOI:10.3969/j.issn.1004-9134.2012.06.024
汪集旸, 胡圣标, 庞忠和, 何丽娟, 赵平, 朱传庆, 饶松, 唐晓音, 孔彦龙, 罗璐, 李卫卫. 2012. 中国大陆干热岩地热资源潜力评估[J]. 科技导报, 30(32): 25-31. DOI:10.3981/j.issn.1000-7857.2012.32.002
曾昭发, 陈雄, 李静, 李桐林, 张良怀. 2012. 地热地球物理勘探新进展[J]. 地球物理学进展, 27(1): 168-178.
周仕新. 2013. 采用CSAMT法在地热勘查中的应用[J]. 企业技术开发, 32(13): 55-57. DOI:10.3969/j.issn.1006-8937.2013.13.029
ГА切列缅斯基. 1982. 实用地热学[M]. 北京: 地质出版社.