西藏日喀则市查孜地热田水化学及同位素特征研究

罗绍强，徐琳，唐华，肖进，胡林

(四川省地质矿产勘查开发局川西北地质队，绵阳 621000)

0 引言

地热资源作为一种绿色清洁能源，其开发利用越来越受到人们的重视[1]。西藏自治区地热资源丰富，地热储量位居我国首位[1-3]。笔者及团队在西藏日喀则地区进行野外地质调查，共发现地热显示区106处[4]。查孜地热田位于日喀则地区西北部，以温泉为主，地表出露2处规模较大的温泉点，最高水温87 ℃，汇总流量达20 000 m3/d[4]，具有较大的开发潜力和较好的利用价值。本文结合前人研究成果并利用现代分析测试手段，在总结温泉形成条件的基础上，研究西藏日喀则市查孜地热田地下热水的水化学及同位素地球化学特征，进一步探讨温泉起源、年龄、热储温度及天然放热量等问题，为该地区地热资源的进一步开发及利用提供基础资料。

1 区域地质背景

研究区位于青藏高原西南部的阿尔卑斯—喜马拉雅巨型构造带东段，大地构造位置属冈底斯—腾冲陆块冈底斯岩浆弧带(图1)。区域上，该带呈EW向展布，南侧以雅鲁藏布江缝合带与弧前盆地接触，主要由白垩纪、古近纪基性、中酸性侵入岩和古近纪中酸性火山岩以及大竹卡组磨拉石建造组成。侵入岩分布于冈底斯岩浆弧带南部，中酸性火山岩呈面状分布，以中心式喷发为主，发育环带状、放射状火山机构，构成一个完整的岩浆岩序列[5-6]。

图1 研究区大地构造位置(a)及地质简图(b)Fig.1 Geotectonic location map(a) and geological sketch (b) of the study area

根据西藏自治区新构造分区[4]，研究区属冈底斯—念青唐古拉断块隆起区当惹雍措—许如措SN向地堑区，为新生代断陷盆地。盆地以第四纪河-湖相堆积物为主，有少量新近纪中新世火山岩沿盆地两侧的活动断裂带分布。该区主要地貌标志是以现代湖泊为中心的向心水系格局。断块山地分布现代冰川和晚更新世冰川堆积物，表明自中更新世末以来，该区一直处于强烈的断块上升阶段，形成多则布(6 436 m)等极高山，耸立于高原腹地。

在许如措东岸，中更新世湖积层存在大量火山灰沉积物，表明新近纪—中更新世，该区存在多次火山活动。地热资源的存在进一步反映该区自中更新世以来新构造运动不但未停止，且有加剧的趋势[7]。

2 样品采集及分析测试

共采集水化学样品6件，其中地热水样品2件，河水样品3件，井水样品1件。具体采样位置见图1。采样严格按照《GB/T 11615—2010地热资源地质勘查规范》[8]执行，样品在西藏自治区地质矿产勘查开发局中心实验室完成测试。主要仪器为Optima 5300电感耦合等离子发射光谱仪、ELAN DRC-e电感耦合等离子体质谱仪、AFS-820原子荧光光谱仪、L2130i同位素分析仪等。分析结果经阴、阳离子平衡检验，误差符合规定，可以满足研究要求。样品主要离子及同位素测试结果见表1。

表1 查孜地热田水化学样品主要离子及同位素含量Tab.1 Main ion and isotope contents of the hydrochemical samples from Chazi geothermal field

3 水化学特征

图2 研究区地下热水Piper图Fig.2 Piper map of the underground hot water in the study area

该地热田冷水和热水离子浓度存在明显差异，但同时具有一定的水力联系。地下热水主要来源于大气降水，降水入渗到地下，在深循环过程中，与沿断裂破碎带运移的地热流体混合，形成地下热水，再沿断裂破碎带排泄至地表形成温泉。因此，断裂破碎带既是温泉的主要通道，也是地热流体的升流储集场所。

4 同位素特征

地下水在演化中，除了形成一般的物理、化学踪迹外，还形成了微观同位素踪迹，记录了地下水的起源及其演化过程[11]。因稳定同位素和放射性同位素对水具有标记和计时作用[12-13]，故利用同位素方法对研究区温泉的起源及年龄进行分析和计算。

4.1 温泉起源

4.1.1 温泉补给来源

1961年，Craig[14]提出了全球大气降水线的概念，通过总结全球各地区大气降水中的δD和δ18O变化范围及二者之间的关系，得出了用于判断地下水现代补给来源的Craig降水线。研究表明，绝大部分地热系统中的热水来自大气降水[13-15]。由表1可知，查孜地热田热水中δD、δ18O及标准平均海水的偏差值分别为-163‰和-20.5‰～-20.6‰，均低于标准平均海水的δD及δ18O值，出现了负偏差，表明该地热水具有现代大气降水及地表水渗入起源的特征。通过与全球大气降水线(Craig降水线)对比，发现研究区2个样品数据投影点均位于当地及Craig降水线下方(图3)，说明查孜地热田地下水的主要补给来源是大气降水和冰雪融水。

图3 氘氧同位素分析值与克雷格雨水线关系[16]Fig.3 Relationship between D,O isotope analysis and Craig′s rain water line[16]

4.1.2 温泉补给高度

根据氧同位素含量与海拔的关系可计算温泉补给高度。其公式为

H=(δs-δp)/K+h，

(1)

式中：H为同位素入渗高度或补给区高度，m；h为取样点(泉)标高，m；δs为地下水的δ18OVSMOW值，‰；δp为取样点附近大气降水的δ18OVSMOW值，‰；K为同位素高度梯度，即海拔每变化100 m时，δ18OVSMOW值的变化量，取值0.3‰/100 m。

经计算，查孜地热田温泉补给高度为海拔5 652 m以上(表2)，属于山地对流型。该区地下热水是山地大气降水和冰雪融水沿断裂带渗入地下后经深循环，被深部岩浆热源加热对流的结果。

表2 温泉补给高度计算结果Tab.2 Calculation of the hot spring replenishment height

4.2 温泉年龄

根据氚的放射性衰变原理[13]可计算温泉年龄。其公式为

N=N0×e-λt，

(2)

式中：λ为氚的衰变常数，由于氚的半衰期为12.43 a，故λ=0.056 a-1；N为氚的当前值，TU；N0为氚的背景值，TU。

由于自然界中的氚含量一般>10 TU，故N0取值10 TU。研究区样品中氚的当前含量<1.0 TU(表1)，N取值1.0 TU，由此计算温泉年龄约为41 a。通常氚的背景值均大于当前背景值[17-18]，故推断该地热田热水在地下运移滞留至少41 a。

5 热储温度估算

地下热储温度是评价地热资源形成机制和开发利用潜力的重要参数，通常选用地球化学温标技术估算热储温度[19]。本次采样为地表水，由于浅层冷水的混合作用，热水体系中化学组分并未达到水-岩平衡状态，因此，在估算地热显示区热储温度之前，通常利用Na-K-Mg三角图解(图4)判别地下水的水-岩平衡状态。由图4可知，研究区仅热水为部分平衡水，河水及井水均为未成熟水。因此，该区不适宜利用阳离子浓度估算热储温度，而应选择SiO2温标进行热储温度计算。石英温标分为有最大蒸汽损失温标和无最大蒸汽损失温标，由于本次热水采样温度均达到当地沸点温度，因此，使用有最大蒸汽损失温标计算热储温度较合理[20]。其计算温度为148.18～153.49 ℃(表3)。

表3 热储温度计算结果Tab.3 Calculation results of the thermal storage temperature

图4 Na-K-Mg三角图解Fig.4 Na-K-Mg triangular diagram

6 天然放热量估算

温泉的天然放热量可根据以下公式[21]估算。

Q=365×q×ρ×c×(t-t0)，

(3)

式中：Q为天然放热量，J/a；q为温泉流量，m3/d，该地热田汇总流量为20 000 m3/d；ρ为热水的密度，988 kg/m3；c为热水的比热容，4 186 J/(kg·℃)；t为热水的温度，取平均值80 ℃；t0为非热异常区恒温层温度，根据该地热田浅层测温结果，取值5 ℃。

根据上述数据，计算查孜地热田温泉的天然放热量为2 264.33×1012J/a，折合标准煤约为77 360 t/a。

7 结论

(2)查孜地热田地下热水的主要补给来源是大气降水和冰雪融水。降水入渗到地下，在深循环过程中，地热流体沿断裂破碎带运移，升流过程中有冷水混入，断裂破碎带为温泉的主要通道，是地热流体的升流储集场所。地下热水的补给高度为海拔5 652 m以上，属于山地对流型。查孜地热田热水在地下运移滞留至少41 a。

(3)查孜地热田热储温度为148.18～153.49 ℃，天然放热量为2 264.33×1012J/a，具有较好的开发利用潜力。