唐海1)，张代磊2, 3)，周文纳1)
1) 兰州大学地质科学与矿产资源学院，甘肃省西部矿产资源重点实验室，兰州, 730000；2) 中国地质科学院，北京，100037；3) 自然资源部深地科学与探测技术实验室，北京，100094
内容提要: 地热能是一种安全、清洁、稳定的可再生能源，广东阳江地区的地热资源开发潜力巨大，但研究程度较低。热红外遥感技术能够高效地圈定地热异常区并对隐伏断裂进行预测，笔者等选取研究区landsat8遥感数据，通过分析地表温度、地热异常点及人类活动的关系，利用大气校正方法对该区潜在的地热异常进行了预测。在排除人为热影响下得到高温地热异常区9处，排除湖泊水域以及海拔的影响下得到低温地热异常区4处。在此基础上，利用地表温度与归一化植被指数计算研究区的温度植被干旱指数值(TVDI)，预测可能存在的隐伏断裂23 处。最后结合岩性、断裂构造有效地圈定了5处地热靶区，主要分布在研究区西南区域和新洲镇附近。研究结果有助于预测阳江地区具有开发前景的地热异常区、识别圈定具有勘探意义的靶区以及预测隐伏断裂，对该区后续地热资源的进一步研究具有重要的意义。
关键词：地热能源；热红外遥感；地表温度反演；断裂预测；靶区圈定
减少煤炭、石油等化石能源的使用是实现“碳达峰、碳中和”国家重大战略目标最直接有效的途径之一。为了降低化石能源的使用，新的清洁能源的大力开发和利用将成为重要的突破口。地热能是新清洁能源的重要组成部分，其储量大、分布广，是我国能源供应体系的重要分支，也是极具潜力的战略性可再生能源。
广东省地处欧亚大陆东南边缘，长期受板块运动综合作用的影响，境内断裂构造发育伴随多期次强烈岩浆活动(闫晓雪等，2019)，是我国地热储量最丰富的省份之一。其内温泉密布，温泉点数量仅次于西藏和云南，全省业已调査出露的温泉达320处，可开采量约23.8万m3/d(袁建飞，2013)。位于广东省的阳江盆地在晚近时期有着强烈的新构造运动，地热资源量丰富，拥有巨大的地热能开采潜力。仅阳江的新洲地热田热功率就达到15.47 MW(王双，2013)，一年热产量约等于19000 t 标准煤(周佐民，2015)。但是目前阳江地区地热研究程度依然较低，仅对新洲地热田进行了勘探，计算了其岩体岩石密度、U、Th、K2O 含量和放射性生热率值(周毅等，2016)，以及对其附属的地下热水的地球化学特征与深部循环进行了研究(汪啸，2018)。地热潜能总体开发程度不高，新洲地热田开采面积不足1 km2(袁建飞，2013)。因此，进一步详细研究阳江盆地的地热分布，进行地热靶区的圈定和预测，对整个阳江地区乃至整个广东省地热的勘查研究都具有重要意义。
作为便捷高效的技术手段，在区域性研究中，遥感手段一直占据着重要的地位。在大面积区域上应用热红外遥感技术可以高效、低成本地圈定地热异常区并划分地热能源靶区，因其覆盖范围广、检测精度高 、数据源丰富以及受地面条件限制小的技术优势，在地热资源勘查上具有较高的实用价值(边宇等，2021)。多名研究者在多地应用热红外遥感技术勘探地热资源取得了较为理想的成果。如任正情等(2021)用热红外遥感进行了江西寻乌地区的地温反演，圈定了11个潜在地热异常区。王康(2020)用热红外遥感对辽宁丹东地区进行了地热反演发现了3处与地热存在有关的地热异常区。杨俊颖等(2021)通过热红外遥感技术对西藏尼木地区进行了地热反演，划分出了21个潜在地热异常区。辛磊等(2021)用热红外遥感进行了石家庄地区的地热反演，圈定1处山地对流型地热远景区和2处沉积盆地型地热远景区。可以看出，热红外遥感技术在地热资源勘探上有着较为普遍的应用，在地热靶区圈定中扮演着重要角色。
为了便捷高效地提取阳江盆地的地热异常信息，判断隐伏断裂位置，帮助识别圈定具有勘探意义的地热能开发靶区笔者等通过收集研究区遥感卫星影像、区域地质图以及大气剖面参数，对研究区遥感影像进行了归一化植被指数(NDVI)计算以及地表比辐射率计算，加入大气剖面参数利用大气校正法对地表温度进行反演，利用反演结果进行了阳江盆地地热靶区的圈定以及隐伏断裂预测。
1 研究区概况及地热条件 1.1 研究区概况研究区位于广东省阳江市东部，包括那吉镇、那龙镇、合山镇、横陂镇 、北惯镇、大沟镇和新洲镇。地理坐标: 东经112°02′～112°22′，北纬 21°48′～22°09′，面积1334.27 km2。总体地势西北东北东南高中间及两边低，地形起伏明显，地表水系发育，有多个湖泊分布。在大地构造上地处华南褶皱系西南段，云开隆起区东南缘，吴川—四会断裂的东侧，地质构造复杂，新构造运动主要以周期性的垂直升降运动为主(广东省地质矿产局，1988)。
如图1所示，研究区在区域地质背景上属于海南岛五指山地层分区(中南地区区域地层编写组，1972)，地层发育，从中元古界到新生界均有地层出露，其地层、和侵入岩由老到新为云开群(Pt2-3Y)、蓟县纪晚期片麻状花岗闪长岩(Jx2γδ)、坝里组(Z1b)、牛角河组(1n)、高滩组(2g)、桂头群(D1-2G)、帽子峰组(D3C1m)、春湾组(D3c)、大湖组(C1d)、孤峰组(P2g)、童子岩组(P2t)、金鸡组(J1j)、中侏罗世第一阶段花岗闪长岩(J21γδ)、晚侏罗世第二阶段二长花岗岩(J23ηγ)、早白垩世第三阶段二长花岗岩(K13ηγ)、百足山组(K1b)、白垩系晚白垩世第一阶段花岗斑岩(K21γπ)、黄岗组(Qphg)、桂洲组(Qhg)、大湾镇组(Qhd)。研究区出露的主要岩性为二长花岗岩和花岗闪长岩，占研究区陆地总面积70%以上(广东省地质矿产局，1988)。
1.2 研究区地热条件地热能源有利开采区具有地质学、地球物理、地球化学以及地热学等多方面指标，具体主要包括火山岩浆活动、大地热流值、热储岩性以及低导热率盖层等(罗璐等，2019；钟振楠等，2021；吴姗姗等，2022)。广东阳江地区是我国主要的近代火山活动区，构造活动发育具有多期大规模的幕式岩浆活动(蔺文静等，2016)，其位于东南沿海又是我国最主要的花岗岩分布区(王德滋和沈渭洲，2003)。广泛存在的花岗岩基底属于较好热储岩性，且一般来说酸性岩石富集放射性元素U、Th、K(汪集暘，2014)，其放射衰变释放出来的能量也是大地热流的主要组成部分，具有较高的大地热流值。此外研究区由于位于东南沿海地带，气候湿润，岩石风化剥蚀剧烈，植被发育，在地层表面拥有一层较厚的风化堆积物，拥有较高的热惯性以及较低的热导率。因此，研究区属于地热能源有利开采区，地热异常区的识别对于进一步勘察地热资源靶区有重要意义。
图1 广东省阳江市东部区域地质图Fig. 1 Regional geological map of eastern Yangjiang City, Guangdong Province
图2 数据处理流程图Fig. 2 Data processing flow chart
2 数据源及研究方法利用遥感进行地表温度反演计算其常用的有大气校正法、单通道算法和分裂窗算法。大气校正法无任何遥感数据的限制，适用于反演绿洲地表温度，但计算过程较为复杂(朱贞榕等，2016)，单通道算法适用于多种传感器，较为简便，不需要大气模拟，操作方法上简单易懂，数据易获得(武姿廷，2021)，但大气水汽含量较高时，单通道算法的反演误差较大(张舒婷，2020)。分裂窗算法是目前热红外遥感反演地表温度中发展最成熟、精度较高的地表温度反演算法(宋挺等，2014)，但由于Landsat8 遥感影像的 TIRS 11 波段存在定标参数误差，对于Landsat8 遥感影像的反演结果不稳定，精度差(杨亮彦，2019)。由于研究区位于东南沿海，大气水含量较高，且使用的数据为Landsat8 遥感数据，因此笔者等选取大气校正法对地表温度进行计算。通过收集研究区的遥感数据影像、大气剖面参数以及区域地质图，掌握研究区的基础地质信息，再对收集到的遥感数据进行辐射定标、大气校正、图像裁剪等预处理，校正因大气吸收、散射、传感器定标、地形等因素造成的误差。利用大气剖面参数以及处理后的遥感数据进行地表温度反演和TDVI指数计算，圈定地温异常区和断裂，最后对数据结果进行分析和讨论(图2)。
2.1 大气校正法基本原理大气校正法的基本原理是通过估计大气对地表热辐射的影响，把这部分大气影响从卫星传感器所观测到的热辐射总量中减去，从而得到地表热辐射强度，再把这一热辐射强度通过普朗克公式换算为地表温度(丁凤和徐涵秋，2006)。卫星传感器接收到的热红外辐射亮度值Lλ的表达式可写为：
Lλ=[εB(TS)+(1-ε)L↓]τ+L↑
(1)
式中，ε为地表比辐射率，TS为地表真实温度(K)，B(TS)为黑体热辐射亮度，L↓为大气向下辐射亮度，τ为大气在热红外波段的透过率，L↑为大气向上辐射亮度。通过简单的公式变换即可得到黑体热辐射亮度B(TS)的表达式，再将其带入普朗克公式：
(2)
即可得：
(3)
K1、K2为常数，对于Landsat 8 TIRS Band10，K1= 774.89 W/(m2·μm·sr)，K2 = 1321.08 K(李壮等，2022)。
τ、L↑、L↓ 3个参数可直接通过NASA公布的网站查询(http://atmcorr.gsfc.nasa.gov)。地表比辐射率ε的估算较为复杂，由于其值大小主要取决于地表的物质结构和遥感器的波段区间(覃志豪，2004)，因此需要将地表结构进行分类。从卫星像元的尺度来看，可以大体视作由水体、裸土、城镇、自然表面和植被5类型构成(何好等，2020)。
其中Landsat8 热红外波段典型地物ε值已知，纯像元可直接赋值计算。水体比辐射率εw=0.99383，裸土比辐射率εs=0.96767，植被比辐射率εv=0.98672，建筑物比辐射率εm=0.964885(何好等，2020)。混合像元需要分为相应地物比例进行计算，将自然地表像元视为纯植被与纯裸土的混合，城镇像元视为建筑物与纯植被像元的混合。如此只需将地物通过NDVI指数和植被覆盖度进行分类，计算出对应区域的地表比辐射率即可根据公式(3)计算出地表温度。
2.2 地表温度反演通过大气校正法对地表温度进行反演，并对地表温度反演结果按不同温度进行分级设色，得到地表温度分级图(图3所示)。
图3 广东省阳江市东部区域地温反演结果Fig. 3 Inversion results of ground temperature in eastern Yangjiang area
研究区内最高温度39℃，最低温度14℃，大部分区域温度在20℃到30℃之间，大致呈正态分布，其中高于30℃和低于18℃区域皆成零星散点分布。22℃到24℃为研究区内最普遍的温度，其区域占研究区总面积的43.72%。将反演图像中温度范围在22℃～24℃的区域视为正常地热区域，温度
研究区内存在大量村镇以及河流湖泊等区域，对于地热异常区的划分有着较为严重的误导。24℃以上区域主要集中在西南部和东北部，都在人类居住区，地表温度30℃以上区域多出现在城镇中央。为了正确划分地温正异常需要排除由人类活动引起伪热异常。由于人类活动热影响随着距离增加而减少，其经验距离2000 m为排除人为热的最小必要距离(杨俊颖等，2021)。以城镇中心为圆心，2000 m为半径，划分研究区人为热缓冲区，再进行地热异常识别。西南部北津港—三丫港一带港口码头密布，频繁活动的船只向外散发的高热量是该区域高温异常的主要原因(汪啸，2018)，因此其地热异常也需要排除。
如图4所示，经过排除人为伪地热影响后得到9处地热正异常区。低温区在研究区西北角、东南角以及中心区域较为集中，中心区域为湖泊水域温度20℃到22℃，由于水的比热容较大，其温度比正常地温低1℃～2℃属于正常范围，因此不认为其区域为地热负异常区域。西北、东北与东南角为海拔较高的山区，温度16℃到18℃。虽然随着海拔升高，气温会下降，但是幅度平均为海拔每上升100 m，气温下降0.6℃，三边山区海拔最高不到200 m，研究区平均海拔为18 m，因此根据计算，其地表温度低于正常地热温度属于负异常，经过划分得到4处地热负异常区(图4中红色方框)。
图4 阳江市东部区域地热异常区圈定图Fig. 4 Delineation map of geothermal anomaly area in eastern Yangjiang area
2.3 温度植被干旱指数计算(TVDI)活动断裂可使地表温度升高造成遥感热红外异常但如果断层规模较小并已处于不活动期，则可能并不会使地面温度升高。并且断层破裂带由于上覆第四纪松散沉积物高孔隙度的影响以及自身破碎的关系将存在高含水率并导致高于基底的热惯性。高热惯量意味着断层破碎带与基岩相比会形成低温区，在热红外影像上反映为低值线状异常区(史超等，2020)，因此地温负异常是指示隐伏断层的一种可能标志。
隐伏断裂同时还具有湿度和植被密度相对较高的特点(史超等，2020)。温度干旱植被指数是通过陆地表面温度反演结果获得的地表温度和归一化植被指数间的相关关系。这一指数同时也能反映地表干旱情况，即地表温度和土壤水分的归一化结果。TVDI的值域为[0,1]，TVDI值越大，土壤湿度越低，TVDI值越小，土壤湿度越高。TVDI负异常区即为低地温和高土壤含水率与高植被区。因此TVDI线状低值区可作为指示隐伏断层的标志(Sandholt et al.，2002)，在缺乏目视解译标志且无法直接观察的覆盖区可用于确定断层位置和走向。表达式为:
(4)
LSTmax和 LSTmin分别是最高和最低地面温度。通过预先处理得到的研究区归一化植被指数影像与地表温度影像，可得到研究区 LST 与NDVI 的相关曲线(图5所示)。
对应点的 NDVI 值以及 LSTmax，LST min投影在LST—NDVI 坐标平面上，形成干边和湿边两条边线。最终点位分布结果近似于三角形。D为干边，W为湿边。其中三角形区域的上下两边分别对应 LST max和 LST min，可分别用直线拟合:
LSTmax= a1+ b1× NDVI
(5)
LSTmin= a2+b2×NDVI
(6)
式中，未知系数a1，a2，b1，b2可用最小二乘法确定。
通过计算出的影像各像元的 TVDI 值，生成 TVDI灰度影像。根据TVDI线状负异常识别隐伏断裂如图6所示结果。
图5 LST—NDVI 特征空间散点图及拟合曲线Fig. 5 LST—NDVI characteristic space scatter diagram and fitting curve
图6 温度植被干旱指数反演结果Fig. 6 Inversion results of temperature vegetation drought index
除去研究区湖泊河流等水域引起的负异常，结合区域地质资料(图1)，对其进行解译，如图6所示，得到23条线性负异常，即推测研究区可能存在23条隐伏断裂，断裂主要位于研究区四角，中心东湖水库附近相对较少，其中深大断裂沿北东向展布绝大多数小型断裂和隐伏断裂沿东西向展布。隐伏断裂的识别对地热靶区的圈定具有重要的参考意义。
3 讨论综合地温反演结果和温度植被干旱指数(TVDI) 反演结果，将地温异常圈定的区域和断裂构造投影到遥感影像中(图7所示)，对靶区进行预测。靶区预测的标准采用研究区内仅开采的新洲地热田的岩性与构造特征。新洲地热田区基底岩性为花岗岩，岩体放射性生热率平均值为 4.0 μW/m3，平均大地热流80 mW/m2(周毅等，2016)，其附近断裂发育且多呈北东向和东西向，这些破碎带是深部热能上涌的良好通道，可在破碎带发育处形成一定宽度的地热水储集带(汪啸，2018)。因此将基底为花岗岩、有断裂发育的地热异常区判定为1级靶区；将基底不是花岗岩、断裂发育和基底是花岗岩但断裂不发育的地热异常区判定为2级靶区；将基地既不是花岗岩、断裂也不发育的地区判定为3级靶区。1级靶区存在地热资源的可能性最大，2级次之，3级可能性最小。
图7 地热异常圈定与断裂构造在遥感中投影Fig. 7 Delineation of geothermal anomaly and projection of fault structure in remote sensing
① 号地热区域表层为二长花岗岩和花岗斑岩，断裂不发育；②号区域表层为桂洲组和大湾镇组和二长花岗岩，断裂发育，桂洲组和大湾镇组的岩性为淤泥质粘土、细砂；③号区域表层为坝里组，岩性为浅海相类复理石泥沙质碎屑岩，断裂发育；④号区域表层为二长花岗岩和坝里组浅海相类复理石泥沙质碎屑岩，断裂发育；⑤号区域表层为桂洲组和云开群组和二长花岗岩，断裂发育，桂洲组和云开群组的岩性为淤泥质粘土、浅海相类复理石碎屑岩；⑥号区域表层为春晚组和大湾镇组，岩性为钙质细砂岩、长石石英砂岩和淤泥质粘土，断裂发育；⑦号区域表层为百足山组、帽子峰组和大湾镇组，断裂发育，岩性为凝灰质砂砾岩、灰绿色砂泥质岩、细砂；⑧号区域表层为桂头群粗碎屑岩，断裂不发育；⑨号区域表层为牛角河组、高滩组和二长花岗岩，断裂发育，牛角河组、高滩组岩性为砂、板岩互层夹硅质岩和变余杂砂岩。
由于花岗岩为侵入岩，其出露地表即代表其基底也为花岗岩，因此将②④⑤⑨号区域判定为一级地热靶区；①③⑥⑦号区域断裂发育但基底岩性未知判定为二级靶区；⑧号区域基底岩性未知且断裂不发育判定为三级靶区。其中新洲地热田与③号区域重叠，代表③号区域基底也是花岗岩，因此一级地热靶区有②③④⑤⑨ 5个区域。
4 结论广东阳江盆地位于东南沿海地带，拥有较高的地热潜能，并且具有较好的热储岩性和低热导率盖层，在商业上属于地热有利开采区。热红外遥感有着覆盖范围广，不受地表观测条件限制，成本低，数据源丰富等优点，不仅能够高效精确的圈定地热异常区，还能对隐伏断裂进行识别。通过热红外遥感进行地温反演提取地热异常区是一种经济高效的方法，因此