南阳市碳酸盐岩地热流体同位素特征分析
2019-06-26龚晓凌
龚晓凌,葛 雁
(河南省地质矿产勘查开发局第二地质环境调查院,河南 郑州 450003)
同位素具有化学稳定性高、不易被岩土吸附、不易生成沉淀化合物、检测灵敏度高的特性。尤其氢氧同位素,它们本身就是水分子的组成部分,其组份的变化真正反映水的演化环境,是理想的示踪剂。因此,同位素在近年来越来越多用于地下水的来源研究。
南阳市在控热构造的作用下有丰富的地热资源储备,为系统研究南阳市地热流体的径流特征、热储温度等特征,采用稳定同位素(D、18O)的组成变化判别地热流体的来源及补给源,采用化学温标计算循环深度,判断热储温度,藉此对南阳市地热重点开采层位流体资源进行综合研究。
1 地质背景概述
1.1 区域地质背景
南阳市位于河南省深大断裂构造集中发育地区,断裂主要展布方位为西北、西向,见图1。根据南阳市地热开发现状,主要开采层位为古元古界秦岭岩群(Pt1Q)雁岭沟岩组碳酸盐岩类热储,其呈西北—东南向夹持于朱夏断裂与商丹断裂之间,东界大致为南阳- 方城断裂,平面形态上呈现北部开放东南收缩的簸箕状,埋深受控热断裂及其次生支脉的控制,埋深从山前的数十米向南部缓倾,在断裂影响带附近陡倾,到南阳市区南部达到1400 m左右[1]。
图1 南阳市地热地质图
1.2 古元古界碳酸盐岩热储岩性特征
雁岭沟组大理岩出露于南阳盆地以西及桐柏县北部,以碳酸盐夹钙泥质碎屑岩沉积为主,岩性主要为各种大理岩夹少量斜长角闪岩及长英质片麻岩等。含钙质岩石(大理岩、变粒岩等)的原岩为白云岩、泥灰岩和钙质砂岩。根据西峡县米坪子沟—千脚沟、内乡县蚌峪沟—西阳山沟实测剖面数据③雁岭沟岩组大理岩占比在为86.47%~89.98%。下伏地层为秦岭岩群郭庄岩组,两者之间整合接触。雁岭沟岩组大理岩多发育两到三组张裂隙,部分地段可见溶蚀现象。张裂隙贯通良好,充填物多为泥质,偶见硅质,裂隙宽度1 cm~25 cm 不等,具备较好的地下水储存、运移条件,地下水回灌条件较好,具有作为良好供暖载体的条件[2]。
1.3 古元古界碳酸盐岩热储分布特征
南阳市已有钻孔没有揭穿该层,最大揭露厚度345 m(体育中心地热井),见表1。
表1 南阳市地热井基本情况表
其厚度由西部山区向东部盆地渐薄,热储结构见图2。
图2 热储结构图
2 地热流体同位素特征
2.1 取样分析情况
2016 年4 月在工作区进行了水样采集,其同位素分析由中国科学院地理科学与资源研究所承担。该研究所是中国最早的国家重点实验室之一,致力于地球信息系统科学的基础理论和方法研究。
2.2 同位素测试结果
南阳市及周边不同热储层地热流体的氢氧同位素测试数据见表2。
表2 南阳及周边地区地热水同位素测试数据一览表
2.3 地热流体来源分析
根据表2 数据绘制稳定同位素散点图,见图3。
图3 南阳及周边地区地热流体稳定同位素散点图
由图3 知:①所有的温泉、地热井均位于大气降水线下,说明补给来源均为大气降水;②鲁山五大温泉表现出从上游到下游,逐渐偏离降水线,显示其补给来源距离逐渐增大;③大理岩地热水偏离降水线,分析认为是因为含氧岩石与热水发生同位素交换反应,使地下水中的δ18O 增高所产生的氧漂移;④氘过量参数d 值越小说明地下水滞留时间越长,地下水封存条件好,还原作用、生物化学作用慢,地下水温度高;同一地区,d 值小意味着地下水的补给量大;源于大气降水的同一含水层中,从补给区到排泄区,d 值逐渐降低。即南阳市碳酸盐岩地热流体源自大气降水,形成于封闭的高温高压环境[3]。
2.4 地热流体补给高程分析
据式(1)计算地热水补给高程:
式中:H 为补给高程;δS 为取样点的地下水δ18O 值;δP 为取样点附近的大气降水的δ18O 值;K 为同位素高程梯度,南阳地区取0.225;H0为采样点高程。
计算结果:古近系砂岩地热流体的补给区高程为1013 m~1131 m 左右,雁岭沟岩组大理岩地热流体的补给区高程为1438 m~1547 m 左右,温泉的补给区高程为3220 m~3618 m左右,显示出补给来源及运移距离的差异。
据补给高程并结合水文地质条件和当地地理条件判断知,南阳地区古近系和元古界大理岩热储地热流体补给区为西部及北部的山区。
3 地热流体循环深度确定
根据地下水水化学分析指标,进行地球化学温标计算,据此推算深部热储的温度及地下水循环深度。
3.1 地热温标计算
根据南阳市地热水多属于中低温热水的实际情况,分别采用Giggenbach1988 年建立的K-Mg 温标及SiO2温标进行地热水温标计算,计算结果见表3。
表3 地热水化学温标计算一览表
从表3 知,研究区古元古界大理岩深部热储温度应在63.42 ℃~79.09 ℃。估算温度高于井底温度的地段,继续向深部钻进有可能达到估算温度,且估算温度接近热源温度。
3.2 循环深度计算
据下式进行循环深度计算:
式中:D 为循环深度;tR为热储温度,采用SiO2温标温度;t0为恒温带温度;g 为地温梯度。
据此计算各井地热流体循环深度,见表4。
表4 地热水循环深度计算一览表
由表4 知,化学温标计算结果与地温梯度推算的相应深度处的温度相比偏高,原因如下:①可能与大理岩的地温梯度取值偏低有关,因本区仅有2 孔有大理岩垂向测温数据,取二者均值并借鉴相关岩性经验数值;②利用化学温标推算的是热储温度,地热井水温多为混合开采的水温,一般会低于井底温度。
4 结论
(1)南阳市碳酸盐岩地热流体的补给来源为西部、北部山区大气降水。
(2)补给区高程为1438 m~1547 m。
(3)南阳市碳酸盐岩地热流体的循环深度为1500 m~2500 m。
(4)碳酸盐岩的热储温度为63.42 ℃~79.09 ℃。