Geochronology and geochemistry of intermediate−acid intrusive rocks in Xinlong area, Sichuan Province and its constraints on the evolution of Ganzi−Litang Ocean
-
摘要:
基于详细的野外调查，对四川新龙地区不同类型侵入岩进行了岩石学、地球化学和锆石U−Pb同位素测年，探讨其岩石成因及构造背景。研究表明，古隆巴花岗闪长岩和日果闪长岩同属准铝质系列，具有低硅碱、高铁镁特征，中等的稀土元素总量及明显的Nb、Ta负异常，属于I型花岗岩，为俯冲期岩浆活动产物，获得日果闪长岩锆石U−Pb年龄为217.1 ± 0.8 Ma，形成时代为晚三叠世。日里隆巴花岗闪长岩和花岗岩属于过铝质系列，具有高硅碱、低铁镁及弱的Nb、Ta负异常，均属于S型花岗岩，但后者经历了更高程度的分异，分别为碰撞期和后碰撞期岩浆活动的产物，获得日里隆巴花岗闪长岩锆石U−Pb年龄为203.6 ± 0.5 Ma，形成时代为晚三叠世末。结合区域资料认为，晚三叠世甘孜−理塘洋盆大规模向西俯冲，少量的弧岩浆岩就位于被动大陆边缘雅江残余盆地；晚三叠世末期，由于洋板块的消亡发生弧−陆碰撞，地壳发生部分熔融，形成强过铝质花岗岩。
Abstract:The Yajiang residual basin is located on the southeastern margin of Tibetan Plateau and is an important part of the Yulong−Bayankala foreland basin in the Qiangtang−Sanjiang orogenic system. Based on detailed field investigations, we carried out petrology, geochemistry and zircon U−Pb isotope dating of different types of intrusive rocks in the Xinlong area of Sichuan Province, to discuss their petrogenesis and tectonic background. The results indicate that both the Gulongba granodiorite and the Riguo diorite belong to the metaluminous granites, with low silica, low alkaline, high Fe−Mg characteristics, intermediate REE total contents, and obvious Nb and Ta negative anomalies. They belong to the I−type granites, which are products of subduction magmatic activity. The Riguo diorite was formed at the Late Triassic with zircon U−Pb age of 217.1 ± 0.8 Ma. The Rililongba granodiorite and granite are the peraluminous rocks, with high silica, high alkaline, low Fe−Mg and weak Nb and Ta negative anomalies. Both of them belong to S−type granite, but the latter experienced higher differentiation degree. They were the products of collisional and post−collisional magmatic activity. The zircon U−Pb age of the Rilongba granodiorite is 203.6 ± 0.5 Ma, and its formation period is the end of the Late Triassic. It is believed that the Ganzi−Litang paleo−ocean basin was subducted westward in the Late Triassic, and a small amount of arc intrusions were located in the Yajiang residual basin on the passive continental margin. At the end of the Late Triassic, the arc−land collision occurred due to the extinction of oceanic plates, which caused partial melting of the crust to form strong peraluminous granites.
-
海南岛地热资源非常丰富，按其成因可分为水热型和干热岩型两大类，水热型地热资源又可分为沉积盆地型地热资源和隆起山地型地热资源，隆起山地型地热资源也称构造裂隙型。海南岛隆起山地型地热资源主要分布在中南部丘陵和山地，受导热、导水断裂控制，能量传递方式以热对流为主，热储为呈带状分布的断裂破碎的构造岩。目前，在海南岛中南部已发现隆起山地型地热田16个、温泉30处（陈颖民，2008；高芳蕾等，2009；赵童，2016；张颖，2019），其中，地热田是经过勘查或勘探后证实具有一定规模的地热资源，温泉指出露地表的单个或多个热矿泉，一般勘查程度低，地热资源规模不明确。海南岛隆起山地型地热资源总体开发利用程度较低，仅有10处地热田和温泉得到开发利用，且利用领域单一，主要用于洗浴、旅游酒店，在工业、农业等领域的开发利用基本处于空白。
海南岛隆起山地型地热资源勘查、研究程度差异较大，但总体程度较低。较详细的地热资源调查是1962年由广东省地质局水文工程地质队开展的1∶20万、1∶50万海南岛水文地质概查，对发现的温泉点进行了调查，并提交海南岛热矿泉调查表。1981年由海南地质大队完成的1∶20万海南岛区域水文地质普查，对全岛的地热资源情况进行了初步总结。海南建省后，先后对万宁兴隆、三亚南田、琼海官塘、儋州蓝洋、澄迈西达、保亭七仙岭等9处地热田和13处温泉点进行了勘查评价工作。前人对海南岛隆起山地型地热资源的研究，主要是依据地热田（温泉）勘查评价报告对单个地热田（温泉）的地热地质、水化学特征和地热成因、资源利用和地热勘查方法开展研究（杨兴沐，2004；陈颖民，2008；高芳蕾等，2009；吴小洁，2015；赵童，2016；徐单，2017；杨峰等，2018；张颖，2019）。但对海南岛隆起山地型地热资源形成的地质背景、地热地质特征、资源赋存规律、成因机制等综合研究和评价相对薄弱，在一定程度上制约了海南岛中南部隆起山地型地热资源的潜力认识和开发利用。
本次研究结合区域地质调查、地球物理、水文地质、区域地质调查及地热田（温泉）勘查、调研等资料，开展海南岛隆起山地型地热资源地质综合研究，通过对典型地热田地热地质要素解剖，基本明确了海南岛隆起山地型地热资源地热地质特征和成因地质模型及赋存规律；利用海南岛水文地质勘查钻孔水温资料和地热田、温泉点水温资料，编制了海南岛恒温层水温等值线图，以分析海南岛隆起山地型地热资源分布与恒温层温度的关系；应用前人编制的海南岛居里等温面深度等值线图计算出地温梯度，利用大地热流值定义得到大地热流值，最后用北部福山凹陷油气勘探资料和收集的海南岛46个地热田数据、温泉水温数据、可采储量数据，应用水文地质法估算了地热资源量，评价了海南岛隆起山地型地热资源潜力；依据温度、可采资源规模，分别对海南岛隆起山地型地热田和温泉进行了分类；综合海南岛隆起山地型地热地质特征、赋存规律及资源潜力评价，指出海南岛东南部和西北部是隆起山地型地热资源有利区。
1. 地热地质背景特征
海南岛具有优越的隆起山地型地热资源地质背景，是中国隆起山地型地热资源较丰富的地区。燕山期，海南岛地壳区域性的伸展减薄，地幔物质上涌，深部热载体（岩浆、热液）沿深大断裂带喷溢地表或停留于地壳浅部，从而形成区域性的高热流值（海南省地质调查院，2017）。新生代，海南岛在南北向拉张应力作用下发生持续扩张，形成一系列深大断裂，这些断裂控制了海南岛构造格局、地层、侵入岩、火山岩、地热、地震等分布（海南省地质调查院，2017）。海南岛地热田和温泉均分布于深大断裂附近或断裂交会部位。
海南岛所处的大地构造位置独特，自20世纪50年代以来，不同学者从不同的大地构造观点、不同的研究目的对海南岛的大地构造单元提出了不同的划分方案（夏邦栋，1991；晁会霞等，2016），尽管认识众说纷纭，但都认为东西向九所-陵水构造带南北两侧古生代在沉积、变质、变形、古生物面貌等方面具有明显的差异性，东西向王五-文教断裂带南北两侧新生代火山作用和沉积作用的差异性，都将这2条东西向断裂带作为一、二级构造单元的边界。参考《中国区域地质表·海南志》（海南省地质调查院，2017）的大地构造划分方案，以东西向九所-陵水断裂带为界，划分出北面的华南褶皱系和南面的南海地台2个一级构造单元，华南褶皱系以王五-文教断裂带为界划分为五指山褶皱带、雷琼断陷2个二级构造单元，五指山褶皱带以北东向的定安-乐东断裂带为界划分为抱板隆起区、五指山隆起区和5个中、新生代断陷成因的凹陷共7个三级构造单元，南海地台只划分为三亚台褶带1个二级构造单元（图1）。海南岛五指山褶皱带的抱板隆起区、五指山隆起区和三亚台缘坳陷3个构造单元，多为中生代侵入岩和古生界变质岩裸露区，岩石具有高的热导率，且深大断裂发育，是大地热流高值异常区。已知的隆起山地型地热田、温泉均分布于这3个构造单元（图1）。
![]()
图 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《中国区域地质志·海南志》（海南省地质调查院，2017）利用重、磁异常资料，结合地质调查和火山岩、温泉分布推断的区域断裂有43条，其中东西向5条、南北向4条、北东向18条、北西向16条；按规模和切割深度划分，切割深度达地幔顶部的壳断裂6条、规模较大的切割基底断裂28条、规模较小的一般断裂9条（图1）。海南岛6条壳断裂控制了隆起山地型地热田或温泉的分布，已知的地热田和温泉均分布于东西向王五-文教断裂以南地区，形成3个东西向和2个北东向地热田和温泉带，即沿东西向尖峰-吊罗断裂带、昌江-琼海断裂带、九所-陵水断裂带分布的地热田和温泉带，沿北东向琼海-三亚断裂带和定安-乐东断裂带分布的地热田和温泉带。
海南岛隆起山地型地热田、温泉与中生代花岗岩有密切关系，全岛目前发现的46处隆起山地型地热田、温泉点中有41处直接出露于中生代花岗岩侵入体中，仅二甲温泉出露于中元古代侵入岩中，而新生代火山岩中未发现地热田和温泉出露。其中分布于二叠纪侵入岩中的地热田（温泉）点有11个，分布于三叠纪侵入岩中的地热田（温泉）点有16个，分布于白垩纪侵入岩中的地热田（温泉）点有14个。沉积岩层中温泉和地热田不太发育，仅有九曲江、沙田、蓝洋、桂根4个温泉或地热田分布于石炭系或志留系浅变质砂岩与板岩中。
2. 地热地质特征
本次利用海南岛地球物理、水文地质、区域地质调查、地热田勘查等资料，对海南岛隆起山地型地热田（温泉）地热地质要素进行了系统分析，初步总结了海南岛隆起山地型地热资源的大地热流背景、地温场、储层、导热和导水断层、地热水温、水量、水化学等地热地质特征，并探讨了海南岛隆起山地型地热资源的成因地质模式。
2.1 恒温层温度分布特征
本次收集了海南岛水文地质勘查323口钻孔和已知的46个地热田、温泉点钻孔水温资料（赵童，2016；张颖，2019），编制完成海南岛近地表（相当于恒温层）水温等值线图（图2）。
海南岛东南部和西北部恒温层温度较高，一般大于27℃，中部及北部恒温层温度较低，一般为25～26.5℃。46个地热田、温泉点均位于恒温层温度高异常区。
2.2 居里等温面与大地热流背景
高维等（2020）利用海南岛高精度航磁数据，通过频率域功率谱法反演出居里等温面的深度，海南岛居里等温面的上隆与坳陷大部分呈带状分布，居里等温面存在两隆两坳，即琼中-琼东南隆起区、琼西北隆起区和琼东北坳陷区、乐东-白沙坳陷区。琼中-琼东南隆起区呈北东向分布于琼海、陵水、万宁、保亭、三亚一带，是海南岛居里面埋深最浅的地区，其居里等温面深度为16～24 km，极值点位于陵水北东，深度仅为16 km；琼西北隆起区呈北东向分布于儋州市西、昌江西，至东方市沿海一带，居里等温面深度变化于17～25 km之间；乐东-白沙坳陷区呈北东向分布在海南岛西南部的乐东、中部白沙、儋州市东南、屯昌北一带，居里等温面深度总体变化于26～32 km之间，显示块状沉降区特征，琼东北坳陷区南部边界大致为王五-文教深大断裂带，北至海南岛东北部海岸，对应琼北的新生代断陷区，呈近块状或带状分布，其居里等温面深度为26～35 km，澄迈—临高一带是海南岛居里等温面埋深最大的地区，最大埋深达35 km，在龙州市一带，坳陷带内分布一个北西向次级隆起，居里等温面深度为26 km左右。
大地热流比其他地热单项参数(温度、地温梯度)更能准确地反映一个地区地热场的基本特点，具有十分重要的理论价值和现实作用。本次研究按照大地热流的定义，利用公式（1）计算出海南岛各点的大地热流值。
q=100Kr⋅K (1) 式中：q为大地热流（mW/m2）；Kr为岩石导热率（cal/cm·s·℃），其值用全岛分布最广的花岗岩平均热导率3.1 W/(m·oK)（胡圣标等，2001）,1 cal/cm·s·℃=418.68 W/m·oK；K为地温梯度（℃/100m）。
地温梯度计算公式如下：
K=(T−t)/H×100 (2) 式中：K为平均地温梯度(℃/100m)；T为居里面温度（550℃）；t为恒温层温度（℃），由海南岛水文地质勘查323口钻孔水温资料和已知的46个地热田（温泉）点钻孔水温资料编制（图2）；H为居里等温面深度（m），按高维等（2020）反演的居里等温面深度取值。
北部福山凹陷的大地热流值根据油气钻井实测6口井的地温梯度与25口井75个岩石样品的实测热导率，计算得到128个热流值点，最后根据前人在海南岛实测的大地热流值（胡圣标等，2001；姜光政等，2016）对计算出的大地热流值进行校正。首次编制完成了海南岛大地热流值等值线图（图3）。
海南岛内大地热流值变化于60～92 mW/m2之间，平均值为76 mW/m2，与全球典型克拉通型大陆区和大洋区的大地热流值相比（高维等，2020；蔺文静等，2016），海南岛总体上具有中等—较高的区域大地热流值。海南岛的大地热流值分布与居里等温面起伏具有较好的相关性，呈北东向或北西向展布（图3）。较高的热流值分布于琼中-琼东南和琼西北2个居里等温面隆起区，琼中-琼东南热流高值区呈北东向展布，位于琼海—万宁—陵水—保亭—三亚一带，大地热流值大于80 mW/m2，最大值在兴隆南，达92 mW/m2，处于东西向尖峰-吊罗断裂与北东向琼海-三亚断裂交会部位；琼西北热流高值区亦呈北东向分布于儋州市西—东方市一带，大地热流值也大于80 mW/m2，最大值在儋州市西，达90 mW/m2，处于东西向王五-文教断裂与北东向博厚-三更断裂交会部位；而居里等温面坳陷区，大地热流值相对较低，如琼北居里等温面坳陷区具有较低的大地热流值，介于60～70 mW/m2之间，最低值区位于福山凹陷中部，与居里等温面埋深等值线一样，在龙州市一带分布1个北西向高值区，大地热流值为80～84 mW/m2，其对应北西向的石山-定安断裂，处于福山凹陷与云龙凸起之间；乐东-白沙居里等温面坳陷区，大地热流值为70～80 mW/m2，变化相对较平缓。
从海南岛大地热流等直线图（图3）可以看出，海南岛已知的46个地热田、温泉点主要分布于东南部和西北部2个大地热流值高值区，说明海南岛隆起山地型地热资源具有高的大地热流背景。
2.3 热储层和热盖层
海南岛隆起山地型地热田、温泉热储层均为断裂破碎的构造岩，控制热储分布的断裂均为张性断层，也是地热田和温泉的导热、导水断层，热储分布于断层破碎带中。热储原岩主要为中生代花岗岩，岩性主要为正长花岗岩、二长花岗岩、花岗闪长岩，只有少数为浅变质沉积岩，岩性为下石炭统南好组石英砂岩、砂岩与板岩和志留系千枚岩、板岩及变质砂岩为原岩的断裂破碎带。
热储层均呈带状构造，为视电阻率低异常带，视电阻率一般小于200 Ω·m，岩石破裂，节理裂隙发育，多具碎裂、碎斑、糜棱、片麻状等结构。热储物性具有强的非均质性，断层上盘物性优于下盘，断层中部优于两侧、浅部优于深部。热储分布范围取决于导热、导水断裂规模，统计官塘、九曲江、蓝洋、七仙岭、南田5个典型地热田热储延伸长度为1.3～3.7 km，宽0.4～1.6 km，面积0.5～2.5 km2，分布于2条断层交会处的热储长宽比较小，单条地层控制的热储长宽比较大，热储埋深不超过500 m。
海南岛隆起山地型地热田热盖层主要为第四系松散沉积层和地表致密花岗岩体。第四系松散沉积层主要为河流冲洪积中粗砂、含砾亚粘土、粘土，其他地区为下伏基岩的坡残积亚粘土或含砾亚粘土，厚0～50 m，除个别砂层外，一般隔水隔热性能较好；地表致密花岗岩体盖层岩性主要为块状花岗岩，岩石裂隙不发育或是闭合性的，或被矿物充填，岩体较完整，具有较好的隔水隔热作用。在热盖层不发育或薄弱处以热水泉排泄，形成温泉。
2.4 导热、导水及控储断层
北东向和近东西向基底断裂是形成海南岛隆起山地型地热田和温泉的导热或导水断裂，如官塘地热田的控热断裂为北东向官塘-蓝山断裂，总体走向北东60°，沿万泉河分布，延伸出地热田区外，区内出露长度约为8 km，该断裂是文昌-琼海深大断裂的组成部分，位于其中部；七仙岭地热田导热和导水断层为近东西向，位于东西向尖峰-吊罗断裂南侧，是其旁侧次级断层，出露长度仅2 km左右，地热田温度等值线呈近东西向沿该断层分布，向周边水温降低，说明该断层为七仙岭地热田的导热和导水断层。
海南岛山地隆起型地热资源的导热、导水断裂旁侧的一般断层或更次级断层往往控制地热田（温泉）的热储层分布，从而控制其分布范围和规模，这些断层多为正断层，热储层是受张扭应力作用的断裂破碎带中的构造岩，并且多发育于正断层的上盘。如官塘、蓝洋和南田地热田热储均为北西向张性正断层的上盘破碎带，九曲江地热田热储为北东向张性正断层的上盘破碎带。
2.5 地热水温
海南岛隆起山地型地热田、温泉水温差异较大，地面泉口或井口温度变化范围为35～95℃（赵童，2016；张颖，2019），参考地热资源地质勘查规范(GB/T11615—2010)分类标准，结合海南岛地热资源利用及气候，提出海南岛地热资源温度分类标准（表1）。按照《地热资源地质勘查规范》(GB/T11615—2010)（中华人民共和国国家质量监督检验检疫总局等，2011）分类标准，热水水温为60℃≤t＜90℃，考虑海南地区热水水温70℃以上才可用于制冷，将热水下限提高至70℃；海南岛恒温层温度平均在27℃左右，故将水温大于27℃作为地热资源。
表 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 洗浴、温室、养殖、农灌 按表1的地热资源温度分类，海南岛中南部构造裂隙型地热水资源中，中温地热田有七仙岭和蓝洋2个，水温分别为95℃和93℃；低温44个，其中低温热水5个、低温温热水33个、低温温水6个（表2）。
表 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 海南岛已知的46个地热田、温泉热水水温与控制地热田、温泉的断裂关系密切。一般温度大于70℃的地热田、温泉均分布于海南岛东南部和西北部2个大地热流值高值区；且多分布于6条壳断裂附近。如七仙岭地热田位于东西向尖峰-吊罗断裂带上，官塘、南田等水温较高的地热田位于北东向琼海-三亚断裂上（图1）；单个地热田热水温度由地热田中部向边界逐渐降低，导热断层与导水断层交会部位是水温最高的部位，导水断层上盘水温往往高于下盘。如官塘地热田在北东向导热断层与北西向导热断层交会处的泉口水温最高为84℃，向周边水温降低，地热田边界过渡为正常水温，水温大于70℃的地热井均位于北西向导水断层上盘。
2.6 水化学特征及类型
统计海南岛5个典型地热田水质分析资料（赵童，2016；张颖，2019），海南岛隆起山地型地热水化学具有以下特征（表3）。
表 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 （1）阳离子以Na+、Ca2+、K+为主，Na+含量25.1～1590 mg/L，Ca2+含量2.4～664 mg/L，K+含量2.8～58.7 mg/L。九曲江和南田地热田Na+、Ca2+含量高，是其他3个地热田的几十倍。
（2）阴离子以Cl−、SO42−、HCO3−为主，Cl−含量12～3410 mg/L，SO42−含量25.9～292 mg/L，HCO3−含量48.8～164 mg/L。九曲江和南田地热田Cl−含量高出其他3个地热田几十倍。
（3）水化学类型较复杂，变化较大，主要有HCO3−Na、HCO3•SO4−Na、HCO3•Cl−Na、Cl−Na•Ca、HCO3−Na•Ca等类型，靠近海边多为Cl−Na•Ca型，七仙岭地热田属HCO3−Na型，蓝洋地热田属HCO3•SO4−Na型，官塘地热田属HCO3•Cl−Na型，九曲江和南田地热田属Cl−Na•Ca型（图4；表3）。
![]()
图 4 海南岛典型地热田水化学类型Piper图Figure 4. Piper diagram of hydrochemical types in typical geothermal fields in Hainan Island（4）热水中F−含量较高，含量为2.4～24 mg/L，达到氟水水质标准，具有医疗价值。
（5）矿化度除九曲江和南田地热田较高外，其他均较低。九曲江地热田矿化度达6.2 g/L，南田地热田为1.73 g/L，是由于热矿水在径流过程中混入海水的原因。
（6）pH值介于7～9之间，为碱性热矿水。
（7）九曲江和南田地热田Na+、Ca2+、Cl−含量高异常及高矿化度的原因是其位于海边，地热水补给除大气降水外，还有海水补给，或热矿水在径流过程中混入海水。
（8）七仙岭、蓝洋、官塘3个远离海边的地热田的Na+、Cl−比值为1.61～5.60，大于海水平均值（0.85），说明地热水为大气成因的渗透水。而九曲江、南田地热田的Na+、Cl−比值分别为0.45和0.53，小于海水平均值，也说明地热水有海水的补给（张玉良等，2019）。
2.7 地热成因模式
高芳蕾等（2009）、赵童（2016）、张颖（2019）利用海南岛地热田和温泉热水同位素δD和δ18O数据进行分析，δD和δ18O均为负值，数据点落在大气降水线上，结果表明，地热田和温泉热水的补给主要为大气降水。热水地球化学分析和水化学类型结果表明，临海的地热田和温泉有海水的补给（表3）；根据海南岛典型地热田放射性同位素分析，热水放射性元素含量低（表4），由其所致的热流增量小，因此区内地热田和温泉热源主要来源于地壳深部或地幔的热能；根据海南岛典型地热田热水地球化学分析资料，应用钾镁温标估算的热储温度在74～118℃之间，钾钠温标估算的热储温度在171～202℃之间（表5）。一般认为，钾镁温标公式代表不太深处热储层中的热动力平衡条件，估算的温度为浅层热储温度，而钾钠温标公式是根据水岩平衡和热动力方程推导的计算公式，适用于估算深部热储温度，在略高于海南岛平均地温梯度（3.5℃/100m）的大地热流背景下、恒温层温度取27℃，得到海南岛典型地热田和温泉热储埋深深度为4～5 km（表4）。
表 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 表 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 综合海南岛已发现的地热田和温泉地热地质特征，可以将区内地热资源的成因模式归纳为：热源主要来源于地壳深部或地幔的热能，大气降雨通过断裂在地壳深部形成热对流循环，水体得到能量，水温升高，水压增大，使热流体通过裂隙上升，于构造岩体破碎段富集为热矿水体。热矿水一部分沿断裂垂直向上流动，在盖层薄弱处以热水泉排泄，形成温泉；另一部分在热矿泉高水位区向四周作水平流动，到热储边界与地表常水混合，水温降低。热储主要为以花岗岩或变质岩为原岩的断裂破碎带，盖层为第四系松散层或较完整的岩体（图5）。而邻海的地热田或温泉由于海水的侵入，与热水混合，则形成高矿化度、高盐度的咸热水，如九曲江地热田和南田地热田。
![]()
图 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)3. 地热资源潜力评价及赋存规律
3.1 地热资源潜力评价
将海南岛已发现或探明的46个隆起山地型地热田和温泉依据水温、可采储量（赵童，2016），利用水文地质法估算可采资源量，计算公式为：
QR1=ρwCwQ(tr−tj) (3) 式中：QR1=QR/t，为单位时间的地热资源量（kcal/d）；Q=Adn/t，为热矿水可采资源量（L/d）；t为时间（d）；ρw为热水密度（kg/m3），查表得到；Cw为水的比热容（1 kcal/kg•℃）；tr为热矿水水温（℃）；tj为年平均气温，海南岛取27℃。
由于地热田和温泉在勘探程度、可采水量存在巨大差异，一般同一处温泉经勘探后，用地热井开采的可采水量是温泉的10倍左右，故将地热田与温泉分开估算。
将热能单位1 kcal/d换算成MW的换算关系是：
1 W=3.6 kj/h=3.6/4.1868 kcal/h=0.8598×24 kcal/d=20.64 kcal/d。所以，1 MW=20.64×106 kcal/d。
通过估算，海南岛已发现或探明的46个隆起山地型地热田和温泉的总可采热水量约为8.1×104 m3/d，可采资源量为26.17×108 kcal/d，折合126.81 MW，相当于标煤373.83 t/d（表5、表6）。其中，16个地热田可采热水量约为6.97×104m3/d，可采资源量为23.03×108 kcal/d，折合111.6 MW，相当于标煤329 t/d（表5）；30个温泉可采热水量约为1.15×104 m3/d，可采资源量为3.14×108 kcal/d，折合15.21 MW，相当于标煤44.84 t/d（表6）。
表 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 根据《地热资源地质勘查规范》(GB/T11615—2010)（中华人民共和国国家质量监督检验检疫总局等，2011）分类标准（表7），按同级地热田是温泉热能的10倍，将海南岛中南部构造裂隙型地热水资源划分为中型、小型2类地热资源。其中，中型地热田6个（可采能量大于10 MW），分别为蓝洋、官塘、九曲江、兴隆、七仙岭、南田地热田，小型地热田10个（可采能量小于10 MW）；中型温泉4个（可采能量大于1 MW），分别为南平、林旺、半坡、高坡岭温泉，小型温泉26个（可采能量小于1 MW）。
表 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 3.2 地热资源赋存规律及勘探开发有利目标区
通过对海南岛已知的46个地热田和温泉地热地质综合研究，尤其是对其中5个典型地热田的地质要素进行系统分析，基本明确了海南岛隆起山地型地热资源的赋存规律。
（1）海南岛隆起山地型地热资源均分布于中南部五指山褶皱带的抱板隆起区、五指山隆起区和三亚台缘坳陷3个构造单元。
（2）断裂是控制区内地热田和温泉最主要的地质要素，东西向、北东向深大断裂是主要导热、导水断裂，旁侧的次级断层往往控制地热田（温泉）的热储层分布，热储均为断层破碎带构造岩。
（3）海南岛地热田、温泉多分布于中生代花岗岩区，受断裂破碎的花岗岩是最主要的热储。
（4）海南岛隆起山地型水温高、规模较大的地热田、温泉主要分布于东南部和西北部2个大地热流值高值区。
（5）海南岛隆起山地型地热田、温泉主要分布在恒温层水温高值区。
因此，海南岛东南部和西北部深大断裂附近及其交会部位、大地热流高值区、中生代花岗岩侵入体出露（或隐伏）区、具有高的地面水温的地区是隆起山地型地热资源勘探开发的有利目标区，可能也是隆起山地型干热岩地热资源的有利区。已知地热田及温泉周边是地热资源勘探开发的有利目标区，如官塘地热田F6、F16两条北西向断裂与F1断裂交会部位具有与官塘地热田相似的地热地质条件，是官塘地热田周边寻找隆起山地型地热田的有利目标区（图5−a）；七仙岭地热田水温高，地热田深部热储温度可达170℃以上，具备深层水热型和干热岩地热资源勘探开发潜力。
4. 结 论
（1）海南岛隆起山地型地热资源的主要控制因素是东西向和北东向深大断裂，深大断裂是导热、导水断裂，旁侧的次级断层控制热储分布，受断裂破碎的花岗岩是最主要的热储。
（2）海南岛6条壳断裂控制了隆起山地型地热田或温泉的分布，已知的地热田和温泉均分布于东西向王五-文教断裂以南地区，形成沿尖峰-吊罗、昌江-琼海、九所-陵水断裂带分布的3个东西向地热田（温泉）带和沿琼海-三亚、定安-乐东断裂带分布的2个北东向地热田（温泉）带。
（3）海南岛隆起山地型地热资源的热源主要来源于地壳深部或地幔的热能，热水补给为大气降水，通过断裂在地壳深部形成热对流循环，于构造岩体破碎段富集为热矿水体，在盖层薄弱处以热水泉排泄，形成温泉。邻海的地热田或温泉由于海水的侵入，与热水混合，形成高矿化度、高盐度的咸热水。
（4）海南岛隆起山地型地热资源丰富，已发现和探明的46个地热田和温泉可采热水量约8.1×104 m3/d，可采资源量26.17×108 kcal/d，折合126.81 MW，相当于标煤373.83 t/d，按地热田和温泉水温、规模分类属于中、低温和中、小型地热田或温泉。海南岛东南部和西北部深大断裂附近及交会部位、大地热流高值区、中生代花岗岩侵入体出露（或隐伏）区、具有高的恒温层水温的地区是隆起山地型地热资源勘探开发有利区。
蓝洋和七仙岭地热田水温高，目前利用程度较低且单一，热能盈余多，建议加强地热资源用于制冷的可行性研究；海南岛东南部具有高的大地热流背景，主要为中生代花岗岩出露区，分布有七仙岭地热田和南平、高峰等温度较高的温泉，具有干热岩地热资源潜力，建议加大干热岩地热资源研究评价和勘探工作。
-
![]()
图 1 研究区大地构造位置（a, b, 据Yang et al., 2015修改）和四川新龙地区侵入岩地质简图（c）
Figure 1. Tectonic location map of the study area (a, b) and geological sketch map of the intrusive rocks in the Xinlong area, Sichuan Province (c)
![]()
图 2 四川新龙地区侵入岩锆石阴极发光（CL）图像（a, c）和U−Pb谐和年龄（b, d）
Figure 2. Zircon CL images (a, c) and U−Pb concordia diagrams (b, d) for the intrusive rocks in the Xinlong area, Sichuan Province
![]()
图 3 四川新龙地区侵入岩SiO2−（Na2O+K2O）图解（底图据Middlemost, 1994）
Figure 3. SiO2−(Na2O+K2O) diagram of intrusive rocks in the Xinlong area, Sichuan Province
![]()
图 4 四川新龙地区侵入岩SiO2−K2O图解（a, 底图据Richter, 1989）和A/CNK−A/NK图解（b, 底图据 Maniar et al., 1989）
Figure 4. SiO2−K2O diagram (a) and A/CNK−A/NK diagram (b) of the intrusive rocks in the Xinlong area, Sichuan Province
![]()
图 5 四川新龙地区侵入岩稀土元素球粒陨石标准化图解（a）和微量元素原始地幔标准化蛛网图（b）（标准化值据Sun et al., 1989）
Figure 5. Chondrite−normalized REE (a) and primitive−normalized trace element (b) patterns of the intrusive rocks in Xinlong area, Sichuan Province
![]()
图 6 四川新龙地区侵入岩源区判别图解（底图据Sylvester, 1998）
Figure 6. Diagram of source area discrimination of the intrusive rocks in the Xinlong area, Sichuan Province
![]()
图 7 四川新龙地区侵入岩构造环境判别图解（底图据Pearce et al., 1984; Pearce, 1996）
a—Y−Nb图解；b—Yb−Ta图解；c—(Y+Nb)−Rb图解；d—(Yb+Ta)−Rb图解；VAG—火山弧花岗岩；ORG—洋中脊花岗岩；WPG—板内花岗岩；Syn−COLG—同碰撞花岗岩；Post−COLG—后碰撞花岗岩
Figure 7. Discrimination diagrams for the tectonic settings of the intrusive rocks in the Xinlong area, Sichuan Province
![]()
图 8 四川新龙地区侵入岩构造环境Hf−Ta−Rb图解（底图据Harris et al., 1986）
Figure 8. The Hf−Ta−Rb diagrams of tectonic environment for the intrusive rocks in the Xinlong area, Sichuan Province
![]()
图 9 四川新龙地区侵入岩形成模式图
a—岛弧型岩浆形成模式；b—碰撞型岩浆形成模式
Figure 9. The forming model of the intrusive rocks in Xinlong area, Sichuan province
表 1 四川新龙地区侵入岩锆石U−Th−Pb同位素测试结果
Table 1 Zircon U−Th−Pb dating results of the intrusive rocks in the Xinlong area, Sichuan Province
分析点 含量/10−6 Th/U 同位素比值 年龄/Ma Pb Th U 207Pb/ 1σ 207Pb/ 1σ 206Pb/ 1σ rho 207Pb/ 1σ 207Pb/ 1σ 206Pb/ 1σ 206Pb 235U 238U 206Pb 235U 238U RG−DN1 RG−DN1−1 13 145 323 0.45 0.0501 0.0027 0.2439 0.0146 0.0347 0.0007 0.3212 211 123 222 12 220 4 RG−DN1−2 10 104 264 0.40 0.0568 0.0051 0.2676 0.0236 0.0344 0.0007 0.2377 483 229 241 19 218 5 RG−DN1−3 13 158 331 0.48 0.0500 0.0025 0.2370 0.0124 0.0344 0.0006 0.3475 195 119 216 10 218 4 RG−DN1−4 12 149 315 0.47 0.0483 0.0032 0.2253 0.0142 0.0344 0.0007 0.3119 122 139 206 12 218 4 RG−DN1−5 10 127 280 0.45 0.0576 0.0052 0.2371 0.0217 0.0302 0.0011 0.3934 522 198 216 18 192 7 RG−DN1−6 13 135 329 0.41 0.0493 0.0032 0.2371 0.0157 0.0346 0.0006 0.2599 161 47 216 13 219 4 RG−DN1−7 7 70 205 0.34 0.0575 0.0066 0.2423 0.0238 0.0318 0.0010 0.3273 509 254 220 19 202 6 RG−DN1−8 9 88 250 0.35 0.0558 0.0033 0.2481 0.0151 0.0323 0.0006 0.3206 443 133 225 12 205 4 RG−DN1−9 13 152 341 0.45 0.0497 0.0028 0.2310 0.0126 0.0341 0.0006 0.3088 189 125 211 10 216 4 RG−DN1−10 9 115 228 0.51 0.0571 0.0045 0.2550 0.0191 0.0330 0.0009 0.3446 498 169 231 15 209 5 RG−DN1−11 13 171 311 0.55 0.0493 0.0031 0.2332 0.0143 0.0347 0.0006 0.2976 161 146 213 12 220 4 RG−DN1−12 14 179 345 0.52 0.0514 0.0029 0.2405 0.0129 0.0345 0.0005 0.2935 257 130 219 11 218 3 RG−DN1−14 12 106 302 0.35 0.0528 0.0030 0.2458 0.0130 0.0344 0.0005 0.2685 324 134 223 11 218 3 RG−DN1−15 14 169 337 0.50 0.0521 0.0027 0.2476 0.0131 0.0344 0.0005 0.2852 300 149 225 11 218 3 RG−DN1−16 10 116 246 0.47 0.0493 0.0033 0.2282 0.0150 0.0343 0.0006 0.2491 161 161 209 12 217 4 RG−DN1−17 8 69 198 0.35 0.0547 0.0034 0.2495 0.0149 0.0341 0.0006 0.3090 467 173 226 12 216 4 RG−DN1−18 9 76 229 0.33 0.0475 0.0034 0.2200 0.0143 0.0344 0.0007 0.3060 72 159 202 12 218 4 RG−DN1−19 9 80 221 0.36 0.0497 0.0046 0.2293 0.0209 0.0344 0.0009 0.2770 189 200 210 17 218 5 RG−DN1−20 9 90 242 0.37 0.0464 0.0031 0.2104 0.0146 0.0335 0.0008 0.3320 20 152 194 12 212 5 RG−DN1−21 8 84 214 0.39 0.0487 0.0034 0.2208 0.0142 0.0341 0.0007 0.3059 200 87 203 12 216 4 RG−DN1−22 9 79 233 0.34 0.0464 0.0050 0.2155 0.0215 0.0338 0.0008 0.2406 20 241 198 18 214 5 RG−DN1−23 26 408 587 0.70 0.0494 0.0024 0.2368 0.0114 0.0348 0.0004 0.2510 169 115 216 9 220 3 RG−DN1−24 5 48 139 0.34 0.0551 0.0051 0.2564 0.0216 0.0345 0.0009 0.2995 417 207 232 17 219 5 RG−DN1−25 10 109 245 0.44 0.0520 0.0034 0.2405 0.0147 0.0342 0.0006 0.2685 287 152 219 12 217 4 RG−DN1−26 17 208 437 0.48 0.0514 0.0023 0.2448 0.0110 0.0347 0.0007 0.4671 261 100 222 9 220 5 DN2539 DN2539−1 95 153 287 0.53 0.0479 0.0017 0.2134 0.0080 0.0323 0.0004 0.3083 100 85 196 7 205 2 DN2539−2 77 124 326 0.38 0.0509 0.0016 0.2223 0.0070 0.0318 0.0004 0.4375 235 72 204 6 202 3 DN2539−3 222 371 517 0.72 0.0521 0.0013 0.2281 0.0062 0.0317 0.0003 0.3753 300 57 209 5 201 2 DN2539−4 100 170 263 0.65 0.0494 0.0015 0.2184 0.0064 0.0321 0.0003 0.2951 165 70 201 5 204 2 DN2539−6 179 296 685 0.43 0.0540 0.0015 0.2385 0.0069 0.0320 0.0004 0.4418 369 59 217 6 203 3 DN2539−7 59 100 196 0.51 0.0482 0.0029 0.2130 0.0129 0.0321 0.0005 0.2478 106 150 196 11 203 3 DN2539−8 348 269 464 0.58 0.0587 0.0013 0.4256 0.0146 0.0523 0.0014 0.7904 567 48 360 10 328 9 DN2539−9 275 455 615 0.74 0.0522 0.0015 0.2317 0.0066 0.0321 0.0003 0.3666 295 69 212 5 204 2 DN2539−10 53 76 161 0.47 0.0522 0.0022 0.2284 0.0100 0.0317 0.0004 0.3133 300 98 209 8 201 3 DN2539−11 137 224 618 0.36 0.0501 0.0015 0.2237 0.0070 0.0323 0.0004 0.4441 198 69 205 6 205 3 DN2539−12 123 196 300 0.66 0.0505 0.0018 0.2266 0.0086 0.0325 0.0004 0.3523 217 83 207 7 206 3 DN2539−13 93 156 274 0.57 0.0533 0.0020 0.2360 0.0086 0.0322 0.0004 0.3575 343 88 215 7 204 3 DN2539−14 112 182 255 0.72 0.0476 0.0018 0.2140 0.0091 0.0323 0.0005 0.3642 76 89 197 8 205 3 DN2539−15 187 285 684 0.42 0.0524 0.0012 0.2336 0.0064 0.0322 0.0004 0.4948 302 49 213 5 204 3 DN2539−16 100 167 280 0.60 0.0478 0.0017 0.2125 0.0071 0.0324 0.0004 0.3614 100 −116 196 6 206 2 DN2539−17 70 112 204 0.55 0.0516 0.0023 0.2266 0.0098 0.0321 0.0004 0.3009 333 100 207 8 203 3 DN2539−18 181 299 658 0.45 0.0493 0.0012 0.2203 0.0058 0.0324 0.0004 0.4965 161 56 202 5 206 3 DN2539−19 234 382 609 0.63 0.0489 0.0012 0.2157 0.0056 0.0319 0.0003 0.3549 146 62 198 5 203 2 DN2539−20 118 194 343 0.57 0.0494 0.0017 0.2170 0.0072 0.0320 0.0003 0.3178 165 78 199 6 203 2 DN2539−21 118 194 475 0.41 0.0504 0.0016 0.2223 0.0070 0.0321 0.0004 0.3867 213 72 204 6 204 2 DN2539−22 76 124 268 0.46 0.0502 0.0017 0.2211 0.0076 0.0321 0.0005 0.4339 206 78 203 6 204 3 DN2539−23 96 163 349 0.47 0.0530 0.0018 0.2340 0.0081 0.0320 0.0004 0.3725 332 78 214 7 203 3 DN2539−24 177 283 588 0.48 0.0506 0.0011 0.2255 0.0051 0.0323 0.0003 0.4725 220 52 206 4 205 2 DN2539−25 105 177 435 0.41 0.0515 0.0016 0.2270 0.0076 0.0319 0.0004 0.4174 265 77 208 6 203 3 
下载: 导出CSV
表 2 四川新龙地区侵入岩主量、微量和稀土元素分析结果
Table 2 Major, trace and rare earth elements contents of the intrusive rocks in the Xinlong area, Sichuan province
元素 日里隆巴岩体花岗闪长岩 日里隆巴岩体花岗岩 古隆巴岩体花岗闪长岩 日果岩体闪长岩 FX2539−1 FX2540−1 FX2544−1 FX2544−2 FX2544−3 FX2545−1 FX2542−1 FX2543−1 FX2543−2 FX2539−2 FX2540−2 FX2540−3 FX2540−4 PM031−5FX1 PM031−7FX1 PM031−6FX1 PM031−6FX2 PM003−9FX1 PM003−9FX2 FX1191−1 FX1192−1 RG−FX1 RG−FX2 SiO2 67.96 67.22 68.75 67.38 67.19 68.35 73.66 66.26 73.39 73.59 74.02 75.27 73.65 61.26 58.92 61.86 59.67 58.02 58.87 56.98 56.16 56.30 56.56 Na2O 2.70 3.00 2.66 2.87 2.81 2.53 3.30 3.60 3.12 3.11 2.94 2.12 2.26 1.52 1.66 1.45 1.51 2.05 2.26 2.35 2.13 2.01 1.88 K2O 3.25 2.38 3.24 2.94 3.14 3.16 4.50 3.85 4.84 5.23 5.93 5.13 4.75 1.39 1.38 1.31 1.37 1.74 1.98 1.45 1.45 1.39 1.34 CaO 3.31 3.76 3.22 3.88 3.87 3.59 1.95 3.05 2.10 1.07 0.99 1.07 1.01 5.53 5.43 5.32 5.45 7.75 7.28 7.36 8.80 8.56 8.64 FeO 2.96 3.24 2.92 3.31 3.21 3.14 1.61 2.32 1.47 0.25 0.25 0.21 0.37 5.04 5.30 5.55 5.83 6.00 5.68 6.30 6.33 6.30 6.39 Fe2O3 0.25 0.34 0.16 0.34 0.42 0.31 0.28 0.52 0.36 0.12 0.09 0.12 0.03 0.40 0.41 0.58 0.36 1.49 1.65 1.81 1.69 1.54 0.96 MgO 1.45 1.70 1.40 1.80 1.75 1.58 0.53 1.19 0.73 0.07 0.06 0.08 0.07 3.01 3.41 3.12 3.26 4.27 4.12 4.34 4.36 4.32 4.88 Al2O3 16.30 16.96 16.02 17.64 17.49 16.07 15.80 18.27 16.33 13.91 14.21 12.03 11.59 13.85 14.44 13.36 13.88 15.51 15.93 16.04 16.16 16.59 15.97 P2O5 0.17 0.17 0.16 0.18 0.17 0.18 0.10 0.16 0.12 0.01 0.02 0.01 0.02 0.18 0.19 0.19 0.17 0.12 0.11 0.16 0.13 0.16 0.16 TiO2 0.71 0.77 0.67 0.81 0.78 0.75 0.28 0.55 0.37 0.06 0.06 0.06 0.07 1.16 1.18 1.22 1.17 0.80 0.79 0.94 0.89 0.89 0.82 MnO 0.11 0.12 0.11 0.12 0.12 0.12 0.06 0.10 0.07 0.01 0.02 0.01 0.01 0.19 0.20 0.19 0.18 0.18 0.16 0.20 0.18 0.21 0.21 烧失量 0.85 1.26 1.13 1.08 1.00 0.99 0.67 0.72 0.59 0.58 0.65 0.71 0.69 1.71 2.33 0.82 0.94 2.51 1.70 2.67 1.96 1.54 1.51 La 32.67 34.72 33.38 31.17 36.34 38.25 30.67 29.25 20.05 6.74 7.96 11.42 7.78 18.92 20.81 21.96 20.34 23.16 23.30 28.14 26.79 25.38 21.56 Ce 60.08 61.84 60.89 56.33 62.30 66.19 56.29 53.23 36.31 15.05 20.07 25.87 19.34 41.43 41.51 42.97 39.75 42.84 43.24 53.45 47.38 48.09 44.43 Pr 6.57 6.38 6.41 6.06 6.44 7.12 5.58 5.56 4.26 2.18 2.83 3.28 2.87 4.94 5.16 4.85 4.80 4.87 5.14 6.11 5.21 4.99 4.51 Nd 25.02 23.27 21.57 22.96 27.06 26.38 22.76 20.81 15.42 9.74 12.83 16.04 13.57 18.36 21.37 20.20 19.03 18.29 17.78 24.81 20.42 21.08 19.16 Sm 4.39 4.32 4.05 4.38 4.46 4.42 4.04 3.97 3.55 2.55 4.12 3.74 3.79 3.83 4.27 4.11 3.76 4.05 4.13 4.98 4.38 3.87 3.84 Eu 1.37 1.35 1.37 1.39 1.46 1.36 1.00 1.22 0.82 0.41 0.34 0.53 0.36 1.23 1.18 1.15 1.06 1.07 0.96 1.11 1.09 1.11 1.07 Gd 4.42 4.16 3.97 4.29 4.53 4.36 4.04 3.97 3.50 2.59 3.80 3.23 3.98 4.20 4.39 4.38 4.23 3.73 4.15 4.84 4.13 3.93 3.65 Tb 0.70 0.68 0.61 0.73 0.73 0.69 0.61 0.66 0.66 0.48 0.85 0.56 0.88 0.77 0.85 0.79 0.77 0.72 0.77 0.94 0.77 0.79 0.73 Dy 4.84 4.83 4.27 5.15 5.09 4.71 3.80 4.73 4.55 3.63 6.61 3.78 7.02 5.97 6.16 6.08 5.60 4.76 4.78 6.41 5.24 5.47 4.84 Ho 0.90 0.84 0.80 0.97 0.94 0.93 0.66 0.89 0.88 0.70 1.32 0.70 1.34 1.19 1.20 1.22 1.15 1.01 1.14 1.33 1.10 1.18 1.09 Er 2.74 2.73 2.35 3.06 2.84 2.61 1.80 2.59 2.80 2.15 4.46 2.32 4.29 3.65 3.87 3.85 3.49 3.17 2.98 4.10 3.27 3.19 2.92 Tm 0.44 0.44 0.35 0.48 0.42 0.41 0.27 0.40 0.42 0.38 0.82 0.37 0.78 0.57 0.55 0.58 0.54 0.46 0.48 0.63 0.53 0.53 0.46 Yb 2.68 2.60 2.17 3.16 2.62 2.68 1.65 2.51 2.80 2.56 5.31 2.48 5.07 3.64 3.97 3.92 3.59 3.10 3.30 4.27 3.26 3.58 3.39 Lu 0.41 0.39 0.32 0.45 0.39 0.38 0.25 0.36 0.42 0.38 0.80 0.41 0.79 0.51 0.56 0.56 0.54 0.45 0.47 0.69 0.52 0.51 0.50 Y 29.32 28.44 25.12 32.39 29.79 30.03 21.23 28.78 29.30 23.49 44.08 22.84 44.08 37.29 37.31 39.95 35.17 27.46 30.12 36.40 31.86 35.53 31.98 Li 63.77 71.92 65.00 68.15 71.49 62.75 57.77 81.56 52.55 13.48 12.61 13.61 13.80 26.80 40.56 28.07 26.64 21.55 25.01 24.72 34.10 30.64 28.78 Be 3.99 4.27 3.29 3.97 3.48 3.54 4.54 5.35 4.84 4.48 5.24 3.98 5.38 2.21 2.21 2.22 1.92 1.53 1.44 1.56 1.48 1.56 1.54 V 44.84 50.44 44.09 55.06 56.77 54.65 13.77 35.75 22.54 1.86 2.00 2.25 1.67 158.89 175.93 171.18 148.14 88.20 84.13 80.23 84.02 85.36 80.52 Cr 37.85 40.20 37.48 32.96 44.99 37.71 12.98 21.96 19.87 9.46 16.98 7.30 13.50 156.02 175.34 145.02 136.94 35.42 38.10 51.21 48.06 50.85 50.31 Ni 8.08 8.71 6.73 9.15 8.97 7.89 1.35 3.51 2.41 1.27 1.24 1.21 1.25 20.73 22.66 17.53 16.31 6.04 7.36 11.87 9.01 9.16 6.66 Cu 8.08 12.41 60.59 5.88 5.08 32.00 1.56 3.64 2.63 16.83 15.19 12.51 21.28 10.59 21.35 16.32 16.20 12.46 12.05 10.92 7.08 7.98 7.80 Zn 69.04 71.23 65.80 73.90 72.58 69.52 33.13 50.17 40.55 9.74 7.49 6.23 8.81 92.98 100.14 93.62 92.57 86.55 86.08 95.91 93.25 90.38 87.98 Ga 21.19 22.09 20.14 21.62 22.77 21.33 19.98 22.21 19.10 15.57 16.32 16.13 17.10 21.57 20.62 20.65 19.96 18.38 17.98 19.35 19.80 19.01 19.06 As 5.77 4.01 3.31 2.67 2.37 0.98 2.06 1.22 2.32 2.49 3.59 2.73 4.81 2.32 3.67 3.73 2.59 0.90 0.96 0.74 0.71 0.84 0.82 Rb 147.72 133.20 139.07 126.30 117.41 135.12 186.23 166.41 200.64 166.92 186.48 190.80 173.58 62.25 59.32 58.21 56.66 65.36 71.59 54.25 54.03 57.55 59.16 Sr 238.00 246.69 243.71 238.20 264.43 218.93 194.13 258.17 176.90 66.88 57.35 89.12 53.74 208.81 218.68 191.28 184.97 254.80 252.10 254.10 265.60 259.90 247.30 Y 29.32 28.44 25.12 32.39 29.79 30.03 21.23 28.78 29.30 23.49 44.08 22.84 44.08 37.29 37.31 39.95 35.17 27.46 30.12 36.40 31.86 35.53 31.98 Zr 29.93 31.51 34.11 19.17 22.58 30.08 78.16 47.39 93.83 55.14 46.50 59.86 38.82 26.69 12.38 33.75 35.68 44.33 43.61 43.60 34.79 42.58 44.53 Nb 14.84 18.68 14.09 15.52 16.48 14.27 17.93 20.05 14.90 10.67 17.62 11.20 18.22 10.15 9.16 12.17 11.51 16.04 12.31 14.56 12.81 14.57 13.67 Cd 0.01 0.01 0.04 0.04 0.05 0.02 0.01 0.03 0.02 0.01 0.01 0.02 0.02 0.09 0.09 0.10 0.08 0.13 0.12 0.17 0.13 0.14 0.14 Cs 7.85 9.81 7.08 7.37 6.22 7.62 4.87 8.40 6.16 5.20 5.74 7.34 5.25 1.90 3.38 2.07 2.33 2.47 2.71 2.86 2.18 3.57 3.46 Ba 709.38 483.52 711.35 611.90 667.33 588.03 646.45 557.33 292.97 125.08 81.94 184.31 66.25 380.33 380.29 364.10 340.92 437.50 464.60 425.50 407.40 389.30 383.60 Hf 1.34 1.10 1.46 1.05 1.18 1.32 2.94 1.83 4.73 2.53 2.23 2.62 1.98 1.83 0.67 2.31 2.02 2.72 2.37 2.60 2.05 2.19 2.14 Ta 1.41 1.55 0.44 0.88 1.16 0.64 1.33 1.81 1.15 1.88 3.92 2.61 3.24 0.27 0.16 0.39 0.33 0.82 0.46 0.44 0.39 0.41 0.43 W 0.55 0.58 0.55 0.88 0.72 2.92 0.38 0.46 0.44 1.27 4.47 4.64 2.44 0.39 0.39 1.64 0.28 0.66 0.37 0.32 0.38 0.29 0.23 Pb 17.94 11.86 17.17 15.98 17.88 14.03 22.51 24.61 34.92 42.33 43.93 36.45 44.86 7.59 8.72 6.10 7.08 13.51 10.93 13.52 9.26 8.49 21.30 Th 13.15 13.43 13.28 12.51 13.37 15.56 15.57 14.71 16.00 15.64 21.39 23.25 21.29 5.97 5.34 7.05 7.34 7.64 7.53 7.53 8.09 4.83 4.36 U 4.45 4.46 3.22 3.32 3.36 3.58 3.83 4.65 8.91 12.87 22.62 19.89 24.24 2.00 1.56 1.26 1.75 1.31 1.38 1.41 0.87 0.65 0.56 注：主量元素含量单位为%，微量和稀土元素含量单位为10−6
下载: 导出CSV
-
Burchfiel B C, Zhen Z L, Liu Y P, et al. 1995. Tectonics of the Longmen Shan and adjacent regions, central China[J]. International Geology Review, 663−735.
Butler R W H, Harris N B W, Whittington A G. 1997. Interactions between deformation, magmatism and hydrothermal activity during active crustal thickening: A field example from Nanga Parbat, Pakistan Himalayas[J]. Mineralogical Magazine, 61(1): 37−52.
Harris N B W, Pearce J A, Tindle A G. 1986. Geochemical characteristics of collision−zone magmatism[J]. Geological Society, London, Special Publications, 19: 67−81.
Hsü K J, Pan G T, Sengör A M C. 1995. Tectonic evolution of the Tibetan Plateau: A working hypothesis based on the Archipelago model of orogenesis[J]. International Geology Review, 37(6): 473−508. doi: 10.1080/00206819509465414
Küster D, Harms U. 1998. Post−collisional potassic granitoids from the southern and northwestern parts of the Late Neoproterozoic East African Orogen: a review[J]. Lithos, 45(1/4): 177−195.
Maniar P D, Piccoli P M. 1989. Tectonic discrimination of granitoids[J]. Geological Society of America Bulletin, 101(5): 635−643. doi: 10.1130/0016-7606(1989)101<0635:TDOG>2.3.CO;2
Middlemost E A K. 1994. Naming materials in the magma/igneous rock system[J]. Earth−Science Reviews, 37(3/4): 215−224.
Pearce J A, Harris N B W, Tindle A G. 1984. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks[J]. Journal of Petrology, 25(4): 956−983. doi: 10.1093/petrology/25.4.956
Pearce J A. 1996. Sources and settings of granitic rocks[J]. Episodes, 19(4): 120−125. doi: 10.18814/epiiugs/1996/v19i4/005
Patino D A E, Mccarthy T C. 1998. Melting of crustal rocks during continental collision and subduction[C]// Hacker B R, Liu J G. When continents collide: geodynamics and geochemistry of ultrahigh pressure rocks. Netherlands: Kluwer, Dordrecht, Academic Publishers, 27−55.
Richter F M. 1989. Simple models for trace element fractionation during melt segregation[J]. Earth and Planetary Science Letters, 77(3/4): 333−344.
Rudnick R L, Gao S. 2003. Composition of the continental crust[C]// Turekian K K, Holland H D. Treatise on geochemistry. Oxford: Pergamon: 1−64.
Sun S S, McDonough W F. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes[C]// Saunders A D, Norry M J. Magmatism in the ocean basins. Geological Society, London, Special Publications, 42(1): 313−345.
Sylvester P J. 1998. Post−collisional strongly peraluminous granites[J]. Lithos, 45(1/4): 29−44.
Taylor S R, McLennan S M. 1986. The continental crust: Its composition and evolution[J]. The Journal of Geology, 94(4): 632−633. doi: 10.1086/629067
Tischendorf G, Paelchen W. 1985. Zur Klassfication von granitoiden/classification of granitoids[J]. Zeit Schrift fuer Geologische Wissenschaften, 13(5): 615−627.
Yin A, Harrison T M. 2000. Geologic Evolution of the Himalayan−Tibetan Orogen[J]. Ann. Rev. Earth Planet., 28(1): 211−280. doi: 10.1146/annurev.earth.28.1.211
Yang L Q, Deng J, Yildirim D, et al. 2015. Structure, geochronology, and petrogenesis of the Late Triassic Puziba granitoid dikes in the Mianlue suture zone, Qinling orogen, China[J]. Geological Society of America Bulletin, 127(11/12): 1831−1854.
Zhang H F, Zhang L, Harris N, et al. 2006. U−Pb zircon ages, geochemical and isotopic compositions of granitoids in Songpan−Garze fold belt, eastern Tibetan Plateau: Constraints on petrogenesis and tectonic evolution of the basement[J]. Contrib. Mineral. Petrol., 152(1): 75−88. doi: 10.1007/s00410-006-0095-2
Zhang H F, Harris N, Zhang L. 2007. A−type adakitic magmatism association in Songpan−Garze Fold Belt, eastern granite and Tibetan Plateau: implication for lithospheric delamination[J]. Lithos, 97(3/4): 323−335.
Zhang C Z, Li B, Cai J X, et al. 2007. A−type granite and adakitic magmatism association in Songpan−Garze fold belt, eastern Tibetan Plateau: Implication for lithospheric delamination[J]. Lithos, 103(3): 562−564.
邓红, 唐渊, 骆志红, 等. 2021. 松潘−甘孜造山带东缘塔公岩体岩石学、同位素年代学特征及其构造意义[J]. 地球科学, 46(2): 527−539. 段志明, 张玉修, 祝向平, 等. 2013. 松潘−甘孜南部玛孜措石英闪长岩的地球化学特征、同位素年龄及其构造意义[J]. 地质学报, 87(12): 1874−1886. 段志明, 李勇, 张毅, 等. 2005. 青藏高原唐古拉山中新生代花岗岩锆石U−Pb 年龄、地球化学特征及其大陆动力学意义[J], 地质学报, 79(1): 88−98. 费光春, 李佑国, 温春齐. 2009. 四川乡城—稻城地区花岗岩地球化学特征及构造背景探讨[J]. 矿物岩石, 29(2): 88−95. 侯增谦, 曲晓明, 周继荣, 等. 2001. 三江地区义敦岛弧碰撞造山过程: 花岗岩记录[J]. 地质学报, 75(4): 484−497. 侯可军, 李延河, 田有荣. 2009. LA−MC−ICP−MS锆石微区原位U−Pb定年技术[J]. 矿床地质, 28(4): 481−492. 李平, 陈隽璐, 张越, 等. 2023. 商丹俯冲增生带南缘土地沟—池沟地区侵入岩形成时代及地质意义[J]. 西北地质, 56(2): 10−27. 梁信之, 谭庆鹄, 师常庆, 等. 1984. 1∶ 20万新龙幅、禾尼乡幅、康定幅区域调查报告[R]. 四川省地矿局区调队: 46−57. 潘桂棠, 王立全, 张万平, 等. 2013. 青藏高原及邻区大地构造图及说明书[M]. 北京: 地质出版社: 7−80. 秦蒙, 严松涛, 文浪, 等. 2019. 甘孜−理塘蛇绿混杂岩带晚三叠世构造演化——来自理塘地区勇杰岩体地球化学、年代学的制约[J]. 地质通报, 38(10): 1615−1625. 时章亮, 张宏飞, 蔡宏明. 2009. 松潘造山带马尔康强过铝质花岗岩的成因及其构造意义[J]. 地球科学−中国地质大学学报, 34(4): 569−584. 王全伟, 王康明. 2004. 四川1∶25万石渠县幅区域地质调查报告[R]. 四川省地质调查院: 45−137. 王全伟. 2008. 川西地区花岗岩及其成矿系列[M]. 北京: 地质出版社: 1−305. 许志琴, 侯立玮, 王宗秀. 1992. 中国松潘−甘孜造山带的造山过程[M]. 北京: 地质出版社: 1−190. 严松涛, 秦蒙, 段阳海, 等. 2019a. 四川理塘地区二叠纪洋岛型岩石组合的识别及其构造意义: 来自岩石学、地球化学和年代学证据[J]. 地质学报, 93(2): 381−393. 严松涛, 段阳海, 谭昌海, 等. 2019b. 甘孜−理塘蛇绿混杂岩带中三叠世洋岛型岩石组合的识别及其构造意义——来自岩石学、地球化学和年代学证据[J]. 地球学报, 40(6): 816−826. 严松涛, 谭昌海, 秦蒙, 等. 2020a. 四川理塘地区二叠—三叠纪硅质岩的地球化学特征及地质意义[J]. 地球学报, 41(4): 504−514. 严松涛, 谭昌海, 段阳海, 等. 2020b. 甘孜−理塘蛇绿混杂岩带中段洋岛型岩石组合的发现及其对构造演化的指示意义[J]. 地质学报, 94(2): 439−449. 严松涛, 吴青松, 李虎, 等. 2021. 甘孜−理塘蛇绿混杂岩带中段理塘地区混杂岩物质组成及其洋盆演化史[J]. 中国地质, 48(6): 1875−1895. 严松涛, 吴青松, 谭昌海, 等. 2022. 四川理塘地区花岗闪长岩特征及其增生楔弧岩浆活动[J]. 中国地质, 49(4): 1295−1308. 严松涛, 吴青松, 朱利东, 等. 2023. 甘孜-理塘蛇绿混杂岩带晚三叠世洋岛型岩石组合识别及其对甘孜-理塘洋盆构造演化的制约[J]. 地质通报, 42(10): 1684−1695. 喻光明, 毛世东, 周振菊, 等. 2024. 西南三江甘孜-理塘洋晚古生代构造演化: 来自理塘蛇绿混杂岩堆晶辉长岩U−Pb年龄的约束[J]. 地质通报, 43(1): 61−75.
计量
- 文章访问数: 963
- HTML全文浏览量: 485
- PDF下载量: 422
