Geological characteristics and potential of geothermal resources in the uplift mountain-type of Hainan Island
-
摘要:研究目的
海南岛隆起山地型地热资源十分丰富，发现了46个地热田、温泉。研究其分布的地质背景、地热地质特征和地质成因及资源潜力，对地热资源赋存规律认识和勘探、开发利用具有重要意义。
研究方法综合区域地质调查、地球物理、水文地质、地热田（温泉）勘查等资料，开展海南岛隆起山地型地热资源地热地质综合研究，系统分析地热地质要素，明确地热地质特征和成因地质模式及资源潜力，指出地热资源勘探开发有利区。
研究结果海南岛中南部隆起山地型地热田（温泉）主要赋存于东西向和北东向深大断裂附近或断裂交会部位，处于大地热流值高值区和恒温层高异常区，受断裂破碎的中生代花岗岩是最主要的热储；地热田（温泉）地下水补给均为大气降水，沿控热断裂与深部大地热流对流，在深循环过程中吸热逐渐形成地热流体，于断裂破碎带热储中富集为热矿水体，在盖层薄弱处排泄出地表则形成热水泉；海南岛隆起山地型地热资源总的可采热水量约8.1×104 m3/d，可采资源量26.17×108 kcal/d，折合126.81 MW，相当于标煤373.83 t/d，按地热资源温度和规模分类，分别划分为中温、低温2类和中、小型2类地热资源。
结论海南岛东南部和西北部深大断裂附近或交会部位、大地热流高值区、中生代花岗岩侵入体出露（或隐伏）区、具有高的地面水温的地区是隆起山地型地热资源勘探开发的有利区。
Abstract:ObjectiveHainan Island is very rich in geothermal resources of the uplifted mountain−type, with 46 geothermal fields and hot springs discovered. Studying the geological background, geothermal geological characteristics, genesis, and resource potential is of great significance for understanding the occurrence regularity, exploration, and development of geothermal resources.
MethodsBased on comprehensive regional geological surveys, geophysics, hydrogeological, regional geological surveys and geothermal fields (hot springs) survey data, to carry out a comprehensive study on geothermal geology of the uplifted mountain−type geothermal resources in Hainan Island . a systematic analysis of geothermal geological elements, basically clear the geothermal geological characteristics, genetic model, and the potential of the resources, pointed out a favorable area for geothermal resource exploration and development.
ResultsGeothermal fields (hot springs) in the south−central uplifted mountainous−type of Hainan Island.are mainly located in the east−west and north−east near the deep−seated faults or the intersection of fractures, in the high value of the heat flow and high anomalies of constant temperature layer, broken by the fracture of the Mesozoic granite is the most important heat storage. Geothermal fields (hot springs) groundwater recharge are atmospheric precipitation, along the heat−controlled fractures and deep earth heat flow convection, heat absorption in the deep cycle process gradually formed geothermal fluid, in the fracture zone thermal storage for the enrichment of thermal mineral water body, in the cover of the weak drainage out of the surface is the formation of hot water springs. Hainan Island uplifted mountain−type geothermal resources, the total amount of recoverable thermal water of about 8.1×104 m3/d, the recoverable resources 26.17×108 kcal/d, equivalent to 126.81 MW, equivalent to the standard coal of 373.83 t/d, according to the temperature of the geothermal resources and the scale of classification, respectively, are divided into medium−temperature, low−temperature and medium−scale, small−scale geothermal resources.
ConclusionsThe areas near or at the intersection of deep−seated faults in the southeast and northwest of Hainan Island, the high value area of geothermal heat flow, the exposed (or hidden) area of Mesozoic granite intrusion, and the areas with high surface water temperature are favorable areas for exploration and development of uplifted mountain−type geothermal resources.
创新点首次系统开展了海南岛王五-文教断裂以南隆起山地型地热资源形成的区域地质构造背景、地热地质特征研究及地热资源综合评价，基本明确了海南岛东南部地热资源分布规律、成因机制及地热资源潜力，指明了地热资源勘探开发有利区，提出了开发利用建议。
-
硒（Se）是人体必需的微量元素（Daniel，2008），被誉为“抗癌之王”（毛香菊等，2021；周墨等，2021）。人体补充Se元素最安全最有效的途径是通过天然富硒农产品等膳食摄入（Banuelos et al., 2015；杨奎等，2018；徐雪生等，2022），而富硒土壤是生产天然富硒产品的先决条件，富硒土地资源的开发利用也越来越受到关注（廖启林等，2020），已成为当前农业地质研究的一大热点（刘健等，2022）。江苏的富硒土壤资源稀缺，土壤Se分布极不均匀，含量最高值与最低值相差数百倍，全省表层土壤的Se含量均值为0.18 mg/kg，约1/2的市县存在不同程度的缺硒现象，成片富硒区主要分布于苏南宜溧低山—丘陵区，多出露煤系地层或产有煤矿，苏北大部分平原区基本未发现大面积富硒土地（廖启林等，2020）。廖启林等（2020）认为，江苏全省富硒来源主要可分为2类，一是与富硒岩石有关的“先天”来源，二是人类活动引起的“后天”来源。对于远离母岩的平原区富硒土壤的成因，大多研究表明（杨泽等，2021；柴冠群等，2022；王仁琪等，2022）与Se元素在地表的迁移密切相关，影响土壤硒富集的因素主要有地形、气候、成土母质、土壤质地、有机质、人为因素等（Umesh et al., 2010；Supriatin et al., 2015; Jones et al., 2017）。前人研究结果对于富硒土壤成因研究具有重要的意义，但对于成陆时间短，成陆期表生环境对于现今地表土壤硒富集的影响鲜见报道。
研究区位于江苏省里下河平原的边缘部位，现为江苏省粮食主产区，被称为江苏省的“粮仓”。江苏省1∶25万多目标区域地球化学调查结果揭示（廖启林等，2007），在里下河边缘潟湖相沉积区土壤Se含量多介于0.3~0.5 mg/kg之间，相比于周围地区，呈明显的相对富集特点，具有发展富硒特色农业的潜力。因此，在该区域寻找和挖掘更多的富硒土壤资源用于富硒农产品的开发具有重要的现实意义。该区第四系覆盖数百米厚，成土母质主要为河流冲积物和海相沉积物，这些物质远离母岩，在长期的风化搬运过程中，理化性质与母岩有很大差异。近期在里下河地区发现了约60 km2的硒富集区，本文在1∶5万土地质量调查的基础上，开展了平原区富硒土壤成因及控制因素研究，为类似地区富硒土壤的寻找和开发提供参考。
1. 材料与方法
1.1 研究区概况
研究区位于江苏省南通市西北部（图1），北纬32°27′~32°44′、东经120°11′~120°31′ 之间，总面积约500 km2（其中农用地面积约430 km2）。研究区属北亚热带季风气候，四季分明，年平均气温14.5℃，年平均降水量1025 mm，年平均相对湿度为60%~80%。区内大体属淮河水系，正常河网水位1.2 m；南部（通扬运河以南）属长江水系，正常河网水位2.5 m；主要河流有新通扬运河（东西向）、串场河、通榆河（均南北向）。地貌上，主体为里下河浅洼圩田平原，南部和东部跨少量高沙平原和滨海平原，其中里下河浅洼区为古潟湖相沉积和黄泛冲积平原区，沉积物主要为粉砂质粘土，多夹泥炭或淤泥质亚粘土；高沙平原为海积冲积相沙嘴沙洲堆积平原，河口相沉积的扬泰古沙嘴东缘由长江冲积而成；东部滨海平原为海积相滨海堆积平原。土壤类型主要为水稻土、滨海盐土、高沙土，其中里下河浅洼区主要土壤类型为水稻土，以渗育型、潴育型、脱潜型为主；南部为高沙土分布区，主要土壤类型有泡沙土、板而沙、夹沙土等；东部滨海平原为潮盐土地区，主要土壤类型为夹沙土、板而沙、粘土等。研究区第四系沉积物覆盖厚度在300 m左右，全新世以来，研究区经历了从海湾到三角洲，又从三角洲到陆地，以及几千年的沉积和地壳升降运动的过程（姜鹏等，2006）。境内成陆最早的是扬泰古沙嘴（扬泰岗地），在南莫镇青墩村、海安镇隆政村吉家墩、曲塘镇刘圩村等多处发现新石器时代遗址和哺乳动物化石，据同位素测定约为6000 a。
1.2 样品采集与测试
根据土地利用方式，结合遥感影像资料，在研究区农用地中采集了1800个表层土壤样品（0~20 cm），平均点密度约4个/ km2；深层土壤柱（0~180 cm）样品200个，平均点密度约1个/2 km2，每柱分5层采样：表层(0~20 cm)、犁底层(20~40 cm)、心土层(50~80 cm)、底土层(80~100 cm)和深层土(150~180 cm)。土壤样品采集时间为2022年1—2月。野外样品采集方法参照《土地质量地球化学评价规范》（DZ/T 0295—2016）（中华人民共和国国土资源部，2016）。表层土壤样以2021年遥感影像上确定的地块图斑中心位置为主样点，同时向四周辐射20~30 m，确定4个子样点（子样点均位于同一土地利用类型范围），混合均匀为1个样品；土壤柱采用深层取土器分土层采集样品。按照采样点位3%的比例采集重复样，由检查人员根据采样航点航迹在原点位处进行采集。样品加工：首先在制样室风干样品，然后用木锤碾压进行粗磨，过孔径10目尼龙筛，过筛后的样品一部分用于测试土壤pH值和粒度，另一部分继续进行细磨（过60目尼龙筛），用于测试有机质和Se含量。样品测试方法见表1。实验室内部质控为每隔12个样品插入1个国家一级标准物质，每10个样品做1个平行双样分析，标准物质测试相对偏差（RE）及平行双样相对偏差（RD）均在规范允许范围内，合格率均为100%。
表 1 土壤样品测试方法Table 1. Test method of soil sample测试指标 样品粒度 测试方法 检出限 RE RD 参照标准 pH ≤2 mm ISE 0.1 0.7~7.9% 1.6~15.8% 《土壤中pH的测定》
（NY/T 1377—2007）粒度分析 ≤2 mm 激光衍射法 0.01% 0.6~15% 2.1~19.5% 《粒度分布 激光衍射法》
（GB/T 19077—2016）有机质（SOM） ≤0.149 mm 重铬酸钾法 0.20‰ 0~4.1% 0.1~10% 《土壤检测 第6部分：土壤有机质的测定》（NY/T 1121.6—2006） 全Se ≤0.149 mm 原子荧光光谱法 0.01
mg/kg0.2~4.8% 3.0~12.8% 《土壤中全硒的测定》
（NY/T1104—2006)1.3 数据分析方法
本文数据分析方法如下所述：①采用Excel进行一般数据统计分析，采用SPSS 17软件进行相关性分析。地球化学分布图采用Surfer 11软件绘制，插值方法选用克里金插值法。②粒径公式：Ф= −log2D，D为颗粒直径(μm)，Ф值越大，表明粒径越细，反之亦然。③母质层沉积相划分方法：以深层土壤(150~180 cm)段平均粒径划分成陆期母质沉积相，其中平均粒径≤5 Ф（粘粒含量低于5%）划为砂坝，5~7 Ф（粘粒含量5%~20%）为潟湖—砂坝过渡相，＞7 Ф（粘粒含量高于20%）为潟湖相。④土壤质地类型划分方法：参照国际制土壤质地分类标准（农业大词典编辑委员会，1998；吴克宁等，2019）进行土壤质地类型分类。⑤多元线性回归分析方法：采用SPSS 17软件回归分析完成，根据土壤Se与其影响因子相关性结果，以土壤Se为因变量，选择与其出现极显著相关的影响因子为自变量，为消除不同因子之间数值的差异，对因变量和自变量进行对数转换；然后通过SPSS 17软件线性回归功能，按照相关系数大小依次引入自变量构建回归方程，检验方法选用德宾-沃森。
2. 结果与讨论
2.1 土壤硒地球化学分布特征
研究区土壤Se及主要理化指标统计结果如表2所示。土壤硒、有机质（SOM）含量由深层至表层逐渐增加，表生富集效应较明显。pH值由深层至表层逐渐减小，至犁底层中值约为7.89，总体上呈碱性，至耕作层pH值下降，总体略呈中性。平均粒径Ф值由深层至表层逐渐减小，总体为上部略粗下部较细。
表 2 土壤Se及主要理化指标特征统计结果Table 2. Characteristics of soil selenium and its main physicochemical indicators土壤层位 耕作层(0~20 cm) 犁底层(20~40 cm) 心土层(50~80 cm) 底土层(80~100 cm) 深土层(150~180 cm) 样本数 n=1800 n=200 n=200 n=200 n=200 Se/(mg·kg−1) 范围 0.07~0.63 0.05~0.37 0.04~0.32 0.04~0.28 0.05~0.23 均值 0.23 0.18 0.13 0.12 0.11 SOM/‰ 范围 6.2~64.5 3.6~28.4 3.2~23.1 1.3~25.2 0.3~25.4 均值 26.9 13.2 8.7 6.7 6.5 pH 范围 4.90~8.48 6.80~8.50 6.89~8.74 7.19~9.07 7.11~9.25 中值 7.24 7.89 8.01 8.23 8.41 平均粒径/Ф 范围 4.97~8.04 5.00~8.71 4.67~8.72 4.91~9.05 4.49~9.23 均值 5.91 6.18 6.27 6.43 6.47 参照《天然富硒土地划定与标识( 试行)》(DZ/T0380—2021)（中华人民共和国自然资源部，2021），对于碱性土壤Se含量高于0.3 mg/kg，即达到了富硒土地Se含量阈值，结合江苏全省和南通市表层土壤Se平均含量分别为0.18 mg/kg和0.16 mg/kg，在《土地质量地球化学评价规范》（DZ/T0295—2016）基础上，对土壤Se元素富集等级进行了一定的微调，分级标准如表3所示。参照表3，研究区各层位土壤Se分级如图2所示。
表 3 土壤Se分级标准Table 3. Soil selenium classification standard指标 缺乏 边缘 足硒 富集区 过剩 Se/（mg·kg−1） ≤0.125 0.125~0.175 0.175~0.3 0.3~3 ＞3 主要特征如下所述：①从深层土—底土层—心土层—犁底层—表层自下而上土壤Se逐渐出现富集，足硒和硒富集区逐渐扩大，尤其是表层和犁底层，足硒级以上区域已占主导地位。②深层土壤中约90%的区域处于边缘－缺硒级，仅在乐百年小镇（海安里下河湿地公园）周边约20 km2处于足硒级；底土层足硒区范围逐渐扩大，至心土层已在墩头－仇湖—吉庆一带形成连片区。③犁底层硒富集区扩大至40 km2（约占全区10%），足硒区约200 km2（约占全区50%），缺硒地区大幅收缩。表层土壤Se富集区（≥ 0.3 mg/kg）面积约66.25 km2，主要分布于北部里下河地区。足硒区域约占全区70%，边缘－缺硒地区主要分布于南部高沙平原、东部滨海平原及北部蒋阳村西南部。④从各层位土壤硒等级分区图看（图2），研究区土壤硒出现强烈的富集效应，尤其是在犁底层—表层土表聚性十分明显。从平面上看，研究区自深层土—表土层，各层位土壤Se含量均表现出不均匀性，总体为北部相对富集，南部和东部较贫化。
2.2 成土母质沉积环境
粒度特征是沉积物物源、地形地貌、水动力条件、搬运距离等综合作用的结果（Mclaren et al., 1985；Gao et al., 2019；梅西等，2020），可反映沉积介质的流体力学性质和能量（Liu，et al.，2006），平均粒径（Mz）指示了沉积物粒径频率分布的中心趋向，大小反映了沉积物的平均动能，故通过沉积物粒度特征可以反演母质层的沉积环境。
2.2.1 母质层粒度特征
据粒度分析结果，研究区深土层和底土层粒度特征如图3所示，主要特征如下。①研究区位于里下河古潟湖的边缘，在成壤前存在明显的沉积环境分异，沉积物粒度特征总体表现为古潟湖—砂坝沉积体系，二者粒度特征存在显著差异。母质层粒度特征均揭示，里下河区存有沉积物较细的潟湖相沉积物，中间为较粗的砂坝分割，其中东侧潟湖发育较完整，西侧则有较粗物质淤塞的趋势。②在潟湖区粒度明显细于砂坝区域（表4），平均粒径6.88 ~9.23 Ф，均值为8.14 Ф，其中在乐百年小镇一带出现8 Ф（约合3.906 μm）以上的细粒物质连片区。粘粒（粒径＜2 μm，Ф≥9）含量一般在20%以上，平均为34.5%；粉粒（2~20 μm，粒径5.75 ~9 Ф）含量30%~51.29%，平均为44%；砂粒（20~2000 μm，粒径−1~5.75 Ф）含量9.28%~39.13%，平均为21.5%。该区域SOM含量相对较高，一般在10‰以上，局部区域达到20‰。Se含量为0.06~0.23 mg/kg，平均为0.14 mg/kg，为区域内含量最高的沉积单元。③潟湖之间为砂坝所分割，该砂坝平均粒径为4.56~7.36 Ф，粘粒含量平均约为2.1%，Se含量平均为0.1 mg/kg。在南部和东部均有较粗的砂坝分隔高沙平原和滨海平原，高沙平原区和滨海平原区平均粒径为5~6 Ф，粘粒含量均低于2%，砂粒含量一般大于50%，该区域SOM含量低，一般在3‰左右，Se平均含量在0.06 mg/kg左右。④对照深土层和底土层粒度特征，东侧潟湖仍较稳定，但沉积物粒度略变粗，平均粒径8 Ф以上较细物质分布区逐渐减小；西侧逐渐分割成大小湖荡，其中西北部的白甸—青墩一带已被较粗物质覆盖，今南莫河一带已与北部彻底被较粗物质分割。上述特征总体上表明，随着海平面下降，海岸线东移，区内受到长江、江淮等河流冲积物的影响，逐渐被较粗物质覆盖。
![]()
图 3 深土层和底土层平均粒径(Mz)平面分布图Figure 3. Mz distribution map of deep soil layer and subsoil layer表 4 研究区各沉积单元粒度参数特征Table 4. Characteristics of particle size parameters of each sedimentary unit in the study area地貌分区 沉积单元 参数 平均粒径/Ф 粘粒/% 粉粒/% 砂粒/% SOM/‰ Se/（mg·kg−1） 里下河
浅洼区潟湖 范围 6.88~9.23 20~51.5 30~51.29 9.28~39.13 1.26~25.36 0.06~0.23 平均值 8.14 34.5 44.04 21.45 10.08 0.14 砂坝 范围 4.56~7.36 0.34~4.22 17.5~79.05 18.29~82.17 1.2~16.76 0.05~0.16 平均值 5.84 2.14 50.89 46.94 5.79 0.1 高沙平原 砂坝 范围 4.49~6.53 0.42~3.48 15.87~71.55 24.97~83.71 1.26~10.36 0.04~0.11 平均值 5.61 1.83 45.66 52.48 3.11 0.06 滨海平原 砂坝 范围 4.93~6.32 0.6~3.16 28.44~64.57 32.27~70.96 1.2~7.26 0.04~0.11 平均值 5.48 1.34 41.59 57.07 3.07 0.06 2.2.2 母质层沉积环境分区
参照国际制土壤分类标准（农业大词典编辑委员会，1998；吴克宁等，2019），研究区深层土土壤类型如图4所示，共识别出4种土壤类型：粘土、粘壤土、粉砂质壤土和砂质壤土，其中粘壤土的分布范围极小，仅零星分布于粘土覆盖区。粘土主要分布于里下河潟湖沉积区，粉砂质壤土全区均有分布，砂质壤土主要分布于南部高沙平原和东部的滨海平原。参照研究区地表高程、成陆历史及地质背景，深土层土壤大体为区内成陆期的地表层位。研究区沉积分区如图5所示，共分布有3种沉积相：潟湖相、砂坝—潟湖过渡相和砂坝相。
表层土壤Se富集区（≥0.3 mg/kg）主要分布于潟湖相及少量潟湖—砂坝过渡相沉积物上，无论是南部里下河地区还是东部滨海平原区的砂坝表层土壤均未形成土壤Se的富集，同时里下河地区南片区表层土壤也未形成Se的富集。因此，表层土壤Se的富集与潟湖相沉积母质密切相关，古潟湖的沉积环境是区内土壤Se富集的基础，一定程度上控制了Se的富集范围。在潟湖分布区的南、北部并未出现土壤Se的富集，表明土壤Se的富集还受其他因素的影响。
2.3 土壤Se富集成因分析
2.3.1 成土母质
土壤Se的含量虽然受多种因素的影响，但在很大程度上取决于成土母质的组成和性质（王莹，2008；杨泽等，2021）。在地表元素迁移的作用下，地球化学元素在局部区域会出现相对富集。有研究认为，成土母质与土壤硒有较大的相关性（章海波，2005；周殷竹等，2020），因此本文分析了表层土Se含量与下部各土层Se含量的相关性（表5）。结果表明，表层土壤Se含量与下部各土层均呈现出极显著的正相关性，相关系数分别为0.79、0.71、0.65和0.6(p＜0.01)，相关程度为强相关，与表层土壤Se含量的相关性总体表现为犁底层＞心土层＞底土层＞深土层，随着土层深度的增加相关系数略有下降，这也反映表层土壤Se含量与母质层有极密切的关系，在母质层富集的区域表层土壤出现富集。因此，母质层的沉积环境对于表层土壤Se的富集有一定的控制作用。
表 5 表层土壤Se含量与下部各土壤Se含量相关系数（n=200）Table 5. Relationship between Se content in surface soil and Se content in lower soil土层/cm 犁底层（20~40） 心土层（50~80） 底土层（80~100） 深土层（150~180） 表层（0～20） 0.79** 0.71** 0.65** 0.60** 注：**表示在0.01水平（双侧）上显著相关 2.3.2 地形条件
地形条件是地表元素迁移重要的影响因素。依据本次200个点土壤剖面获得的地表高程（基准面为1985国家高程基准）数据，编制了研究区地形图（图6）。从图6可以看出，表层土富集区与地表低洼地带高度相关，即表土硒的富集范围为地势低的区域。地表高程与潜水位埋深呈显著正相关（图7），相关系数为0.56 (p＜0.01，n=200)，表明地势越高潜水位埋深越深。表6为地表高程与各层位Se、SOM、粘粒、Mz相关系数，结果表明：①地表高程与各土层Se含量均出现极显著负相关，相关系数均在−0.6(p＜0.01，n=200)左右，地势越高土壤Se含量越低，反之亦然，表明地形条件对Se的富集有重要影响。②地表高程与各土层SOM含量也出现极显著的负相关性，相关系数在−0.5(p＜0.01，n=200)左右，地势越高土壤SOM含量越低，表明地形条件对SOM也有显著影响。③地表高程与各土层粘粒、Mz的关系较复杂，在80 cm的上部土壤，呈微弱的负相关，在80 cm以下土层，相关系数为 −0.3~−0.4(p＜0.01，n=200)，表明地形条件对粘粒物质产生了一定的影响，总体趋势为地势越低的区域，粘粒含量越高。
![]()
图 7 地表高程-潜水位埋深相关关系Figure 7. Correlation between surface elevation and groundwater depth表 6 地表高程与各层位Se、SOM、粘粒、Mz相关系数（n=200）Table 6. Correlation coefficient between surface elevation and Se, SOM, Clay, Mz of each layer项目 土层/cm 0~20 20~40 50~80 80~100 150~180 地表高程-Se −0.600** −0.646** −0.590** −0.621** −0.577** 地表高程-SOM −0.401** −0.519** −0.488** −0.492** −0.539** 地表高程-Clay −0.116 −0.101 −0.247** −0.246** −0.318** 地表高程-Mz −0.178* −0.164* −0.296** −0.281** −0.405** 注：**表示在0.01 水平（双侧）上显著相关；*表示在 0.05 水平（双侧）上显著相关（下同） 综上所述，地形条件对研究区Se的富集也有较明显的影响，同时影响SOM及粘粒物质的含量。
2.3.3 土壤理化指标
土壤理化指标是硒富集重要的影响因素（戴慧敏等，2015；陈锦平等，2018；柴龙飞等，2019），研究区土壤Se含量与主要理化指标土壤有机质（SOM）、pH值、粘粒（Clay）、粉粒（Silt）、砂粒（Sand）及平均粒径（Mz））的相关系数如表7所示。
表 7 各土层Se全量与土壤主要理化性质的相关系数Table 7. Correlation coefficient between total soil Se content and main physicochemical properties土层/cm 有机质 pH 粘粒 粉粒 砂粒 平均粒径 0~20 0.46** −0.04 0.20** 0.23** −0.23** 0.15* 20~40 0.82** −0.64** 0.16* 0.31** −0.39** 0.28** 50~80 0.83** −0.55** 0.26** 0.34** −0.52** 0.39** 80~100 0.86** −0.65** 0.47** 0.15* −0.62** 0.54** 150~180 0.87** −0.64** 0.55** 0.27** −0.67** 0.65** 0~180 0.86** −0.70** 0.36** 0.30** −0.29** 0.11* 分析结果表明：①研究区土壤全硒含量均与有机质（SOM）呈极显著正相关，表土层、犁底层、心土层、底土层、深土层及全土层（0~180 cm）相关系数分别为0.46、0.82、0.83、0.86、0.87、0.86(p＜0.01)。总体表现为在犁底层以下的土层，全硒含量均与SOM呈极显著正相关，说明在成壤过程中有机质对硒具有一定的吸附和固定作用（吴俊，2018；柴冠群等，2022），硒能够以腐殖质缔合的形态存在并在土壤中固定下来，SOM对于土壤Se的富集起决定性的作用。也有研究发现，土壤全硒含量与土壤有机质无明显的相关性（陈俊坚等，2012；周越等，2014）。王松山等（2011）研究认为，在相同母质下，有机质含量越高，Se含量越高。②各土层中，除表土层土壤Se全量与pH值未出现显著相关性外，其他土层均出现显著负相关，犁底−心土层、心土层、底土层、深土层及全土层相关系数分别为−0.64、−0.55、−0.65、−0.64、−0.70 (p＜0.01)。有研究表明，土壤pH 值对硒的影响主要是影响其形态及生物有效性，与其全量相关性不大（周越等，2014）。但也有研究认为， pH 值与Se全量有一定的负相关性，pH可以通过影响土壤中硒的形态和价态，改变土壤硒的迁移转化能力，进而影响土壤Se含量（黄春雷等，2013；刘永贤等，2018）。本文相关性揭示出的下层土壤与pH有一定的负相关，而表层几乎没有关系，中下层土壤环境较稳定，一般不易流失。研究区母质环境为碱性—强碱性土壤环境，SOM与pH有较显著的负相关，而SOM与Se全量有一定的正相关性，中下层土壤出现Se全量与pH值的负相关，很可能是在SOM高含量区，即Se全量高值区，其土壤受SOM的影响pH值相对较低，因此在相关性分析过程中出现了一定的负相关。③土壤全Se与土壤质地出现显著相关性，与粘粒（Clay）的关系在底土层（80~100 cm）和深土层（150~180 cm）出现中等强度正相关，相关系数分别为0.47和0.55 (p＜0.01)，其他层位均为弱相关；与粉粒均为中等强度正相关；与砂粒均为负相关，在心土层（50~80 cm）、底土层（80~100 cm）、深土层（150~180 cm）3个层位出现中等—强相关；与平均粒径均为正相关，其中在底土层（80~100 cm）和深土层（150~180 cm）为强相关，其他层位为弱相关。总体而言，在深层土壤，土壤全硒与质地有较好的相关性，总体表现为细粒物质Se全量高。而在表层土壤，受农业生产及高沙土治理等影响，表层土壤质地已与深层出现明显变化，土壤质地不是影响Se含量的单一因素。
2.3.4 土壤硒富集主要控制因素
研究区潟湖-砂坝沉积体系控制了土壤硒的初始富集范围，地形条件和有机质对硒的富集起决定性作用。研究区地处古潟湖边缘水陆交接处，地势低洼处沉积物粒度较细，有机质含量较高，细粒物质中粘土矿物含量高，有机质和粘土矿物均对土壤硒有一定的固定和富集作用。平面上，在地表元素迁移作用下，砂性土壤Se易流失，Se元素往地势低的区域迁移，而地势低的区域是SOM含量高的地区，SOM对Se有较好的吸附作用和固定作用，硒能够以腐殖质缔合的形态存在并在土壤中固定下来。垂向上由于上部土壤SOM含量总体高于下部，在SOM的吸附作用下，Se元素向表层不断聚集迁移。
为揭示研究区SOM对土壤Se富集的控制程度，本次利用多元线性回归技术分析了SOM对Se富集的控制作用。前文分析表明，研究区与土壤Se富集相关的关键土壤理化因子包括SOM、pH、质地等。土壤质地包括粘粒（Clay）、粉粒（Silt）、砂粒（Sand）及平均粒径（Mz），由于具有一定的同质性，本次选择Mz代表质地指标。以获得的土壤柱5个土层1000个样本数据为基础，通过对数转换，然后以土壤Se为因变量，以SOM、pH、Mz为自变量，应用多元回归分析技术，构建土壤Se与其影响因子的线性模型（表8）。
表 8 全土层（0 ～ 180 cm）Se含量多元回归线性模型Table 8. Multiple regression linear model for Se content in the entire soil layer (0~180 cm)因子 预测方程 (n=1000) R2 p 德宾-沃森值 F 单因子 SOM lg Sesoil＝ 0.480 lg SOM – 1.298 0.747 ＜0.01 1.225 322.52 pH lg Sesoil＝ - 3.96 lg pH + 2.703 0.456 ＜0.01 1.075 84.22 Mz lg Sesoil＝ 0.682 lg Mz – 1.391 0.046 ＜0.01 0.477 17.53 双因子 SOM+pH lg Sesoil＝ 0.449 lg SOM – 0.443 lg pH – 0.873 0.764 ＜0.01 1.203 163.50 SOM+Mz lg Sesoil＝ 0.476 lg SOM + 0.157 lg Mz – 1.419 0.765 ＜0.01 1.238 162.79 三因子 SOM+Mz+pH lg Sesoil＝ 0.436 lg SOM + 0.213 lg Mz
– 0.561 lg pH – 0.9220.779 ＜0.01 1.220 110.87 注：式中，Sesoil为土壤Se含量（单位为mg/kg）；SOM为土壤有机质含量（单位：‰）；Mz为平均粒径（单位：Ф）。R2为回归分析中自变量变异对因变量的解释度，即相对控制程度，范围为0~1；p为显著性水平，p < 0.05为有显著性，p < 0.01为极显著性，p＞0.05为没有显著性；德宾-沃森值为检验变量自相关性的指标，若在0~4之间，符合数据独立性；F值为组间和组内的离差平方和与自由度的比值，在p < 0.01下其值越大表明模型越具统计学意义 （1）分别引入单因子SOM、pH、Mz构建模型，R2分别为0.747、0.456和0.046（p＜0.01），表明SOM、pH、Mz单因子对土壤Se的控制精度分别为74.7%、45.6%和4.6%。
（2）由于研究区土壤剖面pH与SOM呈极显著负相关，且研究区成土第一环境以碱性—强碱性土为主体，pH下降有很大一部分因素是受SOM升高造成的，因此，虽然pH单因子对土壤Se的控制程度为45.6%，但并不代表完全因酸碱度的影响造成Se的富集。为此分别引入SOM + pH、SOM + Mz进入模型，2组模型R2分别为0.764和0.765（p＜0.01），表明双因子下对土壤Se的控制程度分别为76.4%和76.5%，扣除SOM单因子的控制精度，pH和Mz对Se的贡献度分别为1.5%和1.8%。
（3）当同时引入SOM + pH + Mz三因子时，模型R2为0.779（p＜0.01），说明上述3个指标对土壤Se的控制程度为77.9%。对于土壤Se的控制，主体是受SOM的影响，SOM对土壤Se的控制程度为74.7%。
（4）上述模型德宾-沃森值均在0~4之间，表明自变量之间符合独立性，F值均大于10且很小的p值意味着至少有一个自变量对应因变量是显著的。
（5）为验证该模型准确度，本文从1800个表层土壤样品中选取了100个点位实测数据对三因子lg Sesoil＝ 0.436 lg SOM + 0.213 lg Mz – 0.561 lg pH – 0.922 方程进行验证，即以实测的SOM、Mz、pH值经对数转换后代入上述方程，得出该点位土壤Se含量预测值，与实测值进行对比，统计结果如表9所示，模型含量最大偏差为0.13 mg/kg，最大相对偏差约为30%，相对偏差绝对值平均约为17.2%，总体上该模型精度较可靠。
表 9 土壤Se实测值与模型预测值偏差统计Table 9. Deviation between measured soil Se values and model predictions项目 含量偏差/
（mg·kg−1）实测值/
（mg·kg−1）预测值/
（mg·kg−1）相对偏差/% 最大正偏差 0.07 0.22 0.29 31.8 最大负偏差 −0.13 0.44 0.31 −29.5 均值 −0.001 0.24 0.23 −1.1 绝对值均值 0.035 0.24 0.23 17.2 注：含量偏差=预测值−实测值；相对偏差=含量偏差/实测值 3. 结 论
（1）江苏海安里下河地区土壤硒出现了强烈的表生富集效应，富集区主要分布于潟湖相沉积低洼区的表层土壤，富集面积达66 km2，Se最高含量达0.63 mg/kg。
（2）研究区成陆期母质沉积环境主要受潟湖-砂坝体系控制，据粒度特征可分为潟湖相、潟湖-砂坝过渡相、砂坝相，其中潟湖相低洼区出现硒的相对富集，该区域有机质含量高，土壤质地较细，在此基础上发育的土壤进一步富集硒，至表层已达富硒土地标准。
（3）土壤硒与深土母质、有机质、地形条件、土壤质地等均出现极显著的相关性，总体表现为母质Se和有机质含量越高、地势越低、土壤质地越细，土壤中Se含量越高，母质沉积环境和地形控制了土壤硒富集边界，有机质对硒具有重要的吸附作用和固定作用。
（4）研究区富硒土壤是发育于古潟湖相的母质经表生富集作用形成的，成因类型为沉积型（古潟湖相），有机质对于研究区Se的富集起到了决定性作用，对土壤硒的贡献率达70%。区内有机质丰富，土壤硒有稳定的来源，具有发展富硒产业的资源基础。
致谢：工作中得到江苏省地质局海洋院“海安富硒土壤调查项目部”野外一线人员的大力支持和审稿专家的悉心指导，在此一并表示衷心感谢。
-
![]()
图 1 海南岛断裂系统、区域构造单元划分与地热田（温泉）分布图（图中断裂分布据海南省地质调查院，2017修改)
Ⅰ—华南褶皱系；Ⅱ—南海地台；Ⅰ-1—雷琼断陷；Ⅰ-2—五指山褶皱带；Ⅱ-1 —三亚台缘拗陷带；Ⅰ-2（1）—抱板隆起区；Ⅰ-2（2）—白沙坳陷；Ⅰ-2（3）—王五凹陷；Ⅰ-2（4）—五指山隆起区；Ⅰ-2（5）—雷鸣凹陷；Ⅰ-2（6）—阳江凹陷；Ⅰ-2（7）—长昌坳陷；F1—王五-文教断裂；F2—昌江-琼海断裂；F3—尖峰-吊罗断裂；F4—九所-陵水断裂带；F5—定安-乐东断裂；F6—琼海-三亚断裂
Figure 1. Distribution map of fault system, regional tectonic units and geothermal fields (hot springs) in Hainan Island
![]()
图 4 海南岛典型地热田水化学类型Piper图
Figure 4. Piper diagram of hydrochemical types in typical geothermal fields in Hainan Island
![]()
图 5 官塘地热田地热地质图（a）和成因模式图（b）
1—第四系；2—白垩系下统鹿母湾组；3—三叠纪中期花岗岩；4—断层；5—地质界线；6—不整合界线；7—地下水补给；8—大地热对流；9—地下水流向；10—地面温泉及编号；11—模式图剖面位置；12—三叠纪中期花岗岩断裂破碎带热储；13—三叠纪中期花岗岩基岩盖层；14—地下水补给方向；15—地热水排泄方向；16—热水钻孔及编号；17—常温钻孔及编号；18—地热钻孔及编号；19—地热田边界；20—水温等值线；21—村镇
Figure 5. Geothermal geological map (a) and genetic model ideograph of Guantang geothermal field (b)
表 1 海南岛地热资源按温度分类标准
Table 1 Classification standard of geothermal resources by temperature in Hainan Island
温度分级 温度（t）区间/℃ 主要用途 GB/T11615-2010 海南岛 高温地热资源 t≥150 t≥150 发电、烘干、采暖 中温地热资源 90≤t＜150 90≤t＜150 烘干、发电、采暖、制冷 低温地热资源 热水 60≤t＜90 70≤t＜90 采暖、制冷、理疗、洗浴、温室 温热水 40≤t＜60 40≤t＜70 理疗、洗浴、采暖、温室、养殖 温水 25≤t＜40 27≤t＜40 洗浴、温室、养殖、农灌 
下载: 导出CSV
表 2 海南岛地热田、温泉水温统计及分类
Table 2 Statistic of water temperature and classification table of geothermal fields and hot springs in Hainan Island
序号 地热田、温泉名称 水温/℃ 温度分类 序号 地热田、温泉名称 水温/℃ 温度分类 1 七仙岭地热田 95 中温 24 上安温泉 51 低温温热水 2 蓝洋地热田 93 25 二甲温泉 51 3 高坡岭温泉 78 低温热水 26 兴隆玫瑰园地热田 50 4 南平温泉 77 27 千家温泉 48 5 官新温泉 75 28 木棉温泉 48 6 半岭温泉 77 29 陀烈温泉 47 7 高峰温泉 72 30 七叉温泉 47 8 红鞋温泉 68 低温温热水 31 福报地热田 47 9 官塘地热田 68 32 乌坡温泉 46 10 林旺温泉 67 33 海坡地热田 44 11 兴隆地热田 66 34 新村地热田 42 12 九曲江地热田 64 35 桂根温泉 42 13 南田地热田 60 36 五指山市温泉 40.1 14 崖城温泉 59 37 莺歌岭-九所地热田 40 15 油甘温泉 58 38 光雅温泉 40 16 蓝山温泉 58 39 邦溪温泉 40 17 沙田温泉 57 40 八所-感城地热田 40 18 西达地热田 55 41 中沙温泉 38 低温温水 19 石硐地热田 54 42 红岗温泉 37 20 加答温泉 54 43 茄新地热田 36.1 21 凤凰山庄地热田 53 44 石门山温泉 36 22 大田温泉 52 45 石壁温泉 36 23 新街温泉 51 46 南岛温泉 34.8
下载: 导出CSV
表 3 海南岛典型地热田地球化学组成
Table 3 Summary table of hydrochemical analysis of typical geothermal fields in Hainan Island
地热田
名称阳离子/（mg·L−1） 阴离子/（mg·L−1） pH值 矿化度/(g·L−1) 水化学类型 K+ Na+ Ca2+ Mg2+ Fe3+ Cl− SO42− HCO3− SiO42− NO3− F− 七仙岭 2.8 67.2 2.8 0.07 <0.02 12.0 25.9 58.6 128 0.94 12.0 8.98 0.299 HCO3−Na 蓝洋 5.7 84.4 19.2 1.8 <0.02 23.6 66.1 134 151 0.29 14.0 8.32 0.467 HCO3•SO4−Na 官塘 7.9 185 8.8 0.34 <0.02 115 73.8 164 113 0.50 24.0 8.39 0.678 HCO3•Cl−Na 九曲江 58.7 1590 664 4.5 1.0 3410 292 73.2 126 17.3 2.4 7.47 6.21 Cl−Na•Ca 南田 16.0 393 190 2.6 <0.02 741 255 48.8 100 1.7 4.8 7.89 1.73 Cl−Na•Ca
下载: 导出CSV
表 4 海南岛典型地热田放射性同位素分析及热水循环深度估算
Table 4 Radioisotope analysis and the hot water circulation depth estimation of typical geothermal fields in Hainan Island
地热田名称 放射性同位素/（Bq·L−1） 阳离子/（mg·L−1） 热储温度/℃ 热水循环深度估算/km 226Ra/10−3 220Ra 总α 总β K+ Na+ Mg2+ K−Mg温标 K−Na温标 七仙岭 2.8 67.2 0.07 99 171 4.11 蓝洋 14.2 17.78 0.512 0.334 5.7 84.4 1.8 76 203 5.02 官塘 4.9 4.65 0.088 0.414 7.9 185 0.34 106 172 4.16 九曲江 23.42 39.91 0.059 0.073 58.7 1590 4.5 127 164 3.90 南田 78.2 50.57 0.58 0.615 16.0 393 2.6 97 169 4.07
下载: 导出CSV
表 5 海南岛地热田资源量估算及规模分类
Table 5 Resources estimation and classification of geothermal fields in Hainan Island
序号 地热田名称 水温/℃ 可采热水储量/（m3·d−1） 可采资源量/（108kcal·d−1） 可采资源量折合标煤/t 热能/MW 规模分类 1 蓝洋地热田 93 7000 4.4740 63.91 21.680 中型 2 西达地热田（含九乐宫温泉） 55 1200 0.3312 4.73 1.605 小型 3 官塘地热田 68 6700 2.6849 38.36 13.011 中型 4 九曲江地热田 64 7800 2.8297 40.42 13.712 中型 5 兴隆地热田 66 7708 2.9475 42.11 14.283 中型 6 兴隆玫瑰园地热田 50 2367.36 0.5380 7.69 2.607 小型 7 茄新地热田 36.1 5774.98 0.5224 7.46 2.532 小型 8 新村地热田 42 2265.28 0.3371 4.82 1.634 小型 9 七仙岭地热田 95 4040 2.6420 37.74 12.803 中型 10 石硐地热田 54 1032 0.2747 3.92 1.331 小型 11 南田地热田 60 8030 2.6051 37.22 12.624 中型 12 凤凰山庄地热田 53 4640 1.1891 16.99 5.762 小型 13 海坡地热田 44 2972.5 0.5004 7.15 2.425 小型 14 福报地热田 47 1385.47 0.2740 3.91 1.328 小型 15 莺歌岭-九所地热田 40 3934 0.5074 7.25 2.459 小型 16 八所-感城地热田 40 2887.48 0.3724 5.32 1.805 小型 合计 6.97×104 23.03 329 111.6
下载: 导出CSV
表 6 海南岛温泉资源量估算及规模分类
Table 6 Resources estimation and classification of hot springs in Hainan Island
序号 地热田名称 水温/℃ 可采热水储量/（m3·d−1） 可采资源量/（108kcal·d−1） 可采资源量折合标煤/t 热能/MW 规模分类 1 沙田温泉 57 541.44 0.1598 2.28 0.775 小型 2 加答温泉 54 70.85 0.0189 0.27 0.091 小型 3 桂根温泉 42 216 0.0321 0.46 0.156 小型 4 红岗温泉 37 603.36 0.0600 0.86 0.291 小型 5 官新温泉 75 296 0.1385 1.98 0.671 小型 6 蓝山温泉 58 25.92 0.0079 0.11 0.038 小型 7 石壁温泉 36 61.92 0.0055 0.08 0.027 小型 8 乌坡温泉 46 494.4 0.0930 1.33 0.451 小型 9 上安温泉 51 463.68 0.1099 1.57 0.533 小型 10 油甘温泉 58 26 0.0079 0.11 0.038 小型 11 南平温泉 77 1228.8 0.5976 8.54 2.896 中型 12 红鞋温泉 68 174.53 0.0700 1.00 0.339 小型 13 高峰温泉 72 401.76 0.1768 2.53 0.857 小型 14 林旺温泉 67 1054 0.4134 5.91 2.003 中型 15 五指山市地热点 40.1 384 0.0499 0.71 0.242 小型 16 半岭温泉 77 563.24 0.2742 3.92 1.329 中型 17 南岛温泉 34.8 850.18 0.0659 0.94 0.319 小型 18 崖城温泉 59 157.92 0.0497 0.71 0.241 小型 19 千家温泉 48 426 0.0885 1.26 0.429 小型 20 石门山温泉 36 1382 0.1236 1.77 0.599 小型 21 中沙温泉 38 25.92 0.0028 0.04 0.014 小型 22 陀烈温泉 47 75.12 0.0149 0.21 0.072 小型 23 高坡岭温泉 78 667.9 0.3318 4.74 1.608 中型 24 二甲温泉 51 129.6 0.0307 0.44 0.149 小型 25 大田温泉 52 254.02 0.0627 0.90 0.304 小型 26 新街温泉 51 174.53 0.0414 0.59 0.201 小型 27 七叉温泉 47 157.25 0.0311 0.44 0.151 小型 28 邦溪温泉 40 256.56 0.0331 0.47 0.160 小型 29 木棉温泉 48 38.88 0.0081 0.12 0.039 小型 30 光雅温泉 40 298.94 0.0386 0.55 0.187 小型 合计 1.15×104 3.14 44.84 15.21
下载: 导出CSV
表 7 水热型中低温地热田和温泉规模分类标准
Table 7 Classification standards of low temperature geothermal fields and hot springs
地热资源
规模地热田分类 温泉分类 热能/MW 保证开采
年限/年热能/MW 保证开采
年限/年大型 > 50 100 >5 100 中型 10∼50 100 1∼5 100 小型 < 10 100 < 1 100 注：地热田规模分类标准按《地热资源地质勘查规范》(GB/T11615—2010)；温泉规模分类标准按同级地热田的1/10
下载: 导出CSV
-
Chao H X, Han X H, Yang Z H, et al. 2016. New exploration of geotectonic characteristics of Hainan Island[J]. Earth Science Frontier, 23(4): 200−211 (in Chinese with English abstract).
Chen Y M. 2008. Present Situation of Geothermal Resource in Hainan Island and Suggestions for Development and Exploitation[J]. Scientific and Technological Management of Land and Resources, 25(6): 61−65 (in Chinese with English abstract).
Gao F L, Yang X Q, Wu G A, et al. 2009. Characteristics of thermal springs and genesis of thermal underground waters in Hainan Island[J]. Journal of Jilin University (Earth Science Edition), 39(2): 281−287.
Gao W, Guo Z H, Zhou J X, et al. 2020. High precision aeromagnetic characteristics and Curie depth analysis of the Hainan Island[J]. Acta Geologica Sinica, 94(11): 3249−3262 (in Chinese with English abstract).
Hainan Geological Survey. 2017. Regional geology of Hainan Province[M]. Beijing: Geological Publishing House (in Chinese).
Hu S B, He L J, Wang J Y. 2001. Compilation of heat flow data in The China Continental area (3rd edition)[J]. Chinese Journal of Geophysics, 44(5): 611−626 (in Chinese with English abstract).
Jiang G Z, Gao P, Rao S, et al. 2016. Compilation of heat flow data in the continental area of China (4th edition)[J]. Chinese Journal of Geophysics, 59(8): 2892−2910 (in Chinese with English abstract).
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[J]. Acta Geologica Sinica, 90(8): 2043−2058 (in Chinese with English abstract).
Wu X J, Zhang Q, Chen C L, et al. 2015. The application of comprehensive electrical prospecting to exploration of thermal mineral water in Wuzhishan area[J]. Chinese Journal of Engineering Geophysics, 12(3): 348−353 (in Chinese with English abstract).
Xia B D, Shi G Y, Fang Z, et al. 1991. The Late Palaeozoic rifting in Hainan Island, China[J]. Acta Geologica Sinica, 65(2): 103−115 (in Chinese with English abstract).
Xu D. 2017. Hydrogeochemistry of geothermal field of Long Mu Bay, Hainan Province[D]. Master Thesis of East China University of Technology: 1−65(in Chinese with English abstract).
Yang F, Ruan M, Zhang D Q, et al. 2018. Study on the geochemical characteristics of hot mineral water isotope in Haipo district, Sanya City, Hainan Province[J]. Ground water, 40(4): 15−17 (in Chinese with English abstract).
Yang X M, Liu C N, Yang Z Z. 2004. The application of three−dimensional resistivity view to disclosing structural form of the Nantian geothermal field[J]. Geophysical & Geochemical Exploration, 28(1): 62−64 (in Chinese with English abstract).
Zhang Y. 2019. A study of the characteristics and formation of the hot springs in Hainan Island[D]. Master Thesis of China University of Geosciences (Beijing) : 1−56(in Chinese with English abstract).
Zhang Y L. 2019. Analysis of geothermal geological conditions in the area of Liudongying of Huairen sag in Datong Basin[J]. Underground Water, 41(5): 15−17(in Chinese with English abstract).
Zhao T. 2016. Characteristics and evaluation of geothermal resources in Hainan Province [D]. Master Thesis of Guilin University of Technology: 1−79(in Chinese with English abstract).
晁会霞, 韩孝辉, 杨志华, 等. 2016. 对海南岛大地构造特征的新探索[J]. 地学前缘, 23(4): 200−211. 陈颖民. 2008. 海南岛地热资源现状及勘查开发利用建议[J]. 国土资源科技管理, 25(6): 61−65. doi: 10.3969/j.issn.1009-4210.2008.06.013 高芳蕾, 杨小强, 吴国爱, 等. 2009. 海南岛温泉特征与地下热水成因[J]. 吉林大学学报(地球科学版), 39(2): 281−287. 高维, 郭志宏, 周坚鑫, 等. 2020. 海南岛高精度航磁特征与居里等温面深度分析[J]. 地质学报, 94(11): 3249−3262. doi: 10.3969/j.issn.0001-5717.2020.11.006 海南省地质调查院. 2017. 中国区域地质志·海南志[M]. 北京: 地质出版社. 胡圣标, 何丽娟, 汪集旸. 2001. 中国大陆地区大地热流数据汇编(第三版)[J]. 地球物理学报, 44(5): 611−626. doi: 10.3321/j.issn:0001-5733.2001.05.005 姜光政, 高堋, 饶松, 等. 2016. 中国大陆地区大地热流数据汇编(第四版)[J]. 地球物理学报, 59(8): 2892−2910. 蔺文静, 甘浩男, 王贵玲, 等. 2016. 我国东南沿海干热岩赋存前景及与靶区选址研究[J]. 地质学报, 90(8): 2043−2058. 吴小洁, 张前, 陈长亮, 等. 2015. 综合电法勘探在五指山地区找热矿水中的应用[J]. 工程地球物理学报, 12(3): 348−353. 夏邦栋, 施光宇, 方中, 等. 1991. 海南岛晚古生代裂谷作用[J]. 地质学报, 65(2): 103−115. 徐单. 2017. 海南省龙沐湾地热田的水文地球化学研究[D]. 华东理工大学硕士学位论文: 1−65. 杨峰, 阮明, 张东强, 等. 2018. 海南省三亚市海坡地区热矿水同位素地球化学特征研究[J]. 地下水, 40(4): 15−17. 杨兴沐, 刘川宁, 杨忠政. 2004. 三维电阻率视图揭示南田地热田构造形态[J]. 物探与化探, 28(1): 62−64. 张颖. 2019. 海南岛温泉特征及成因研究[D]. 中国地质大学(北京)硕士学位论文: 1−56. 张玉良. 2019. 大同盆地怀仁凹陷柳东营一带地热地质条件分析[J]. 地下水, 41(5): 15−17. 赵童. 2016. 海南省地热资源特征及评价[D]. 桂林理工大学硕士学位论文: 1−79. 中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会. 2011. 地热资源地质勘查规范(GB/T11615—2010)[S]. 北京: 中国标准出版社. -
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 1209
- HTML全文浏览量: 257
- PDF下载量: 185
- 被引次数: 1
