潮间带地热田水化学特征及成因机制分析.pdf

2023年10 月引用格式：刘润方，李理.潮间带地热田水化学特征及成因机制分析J.水利水电快报,2 0 2 3,44（10）：10 4-10 9.水利水电快报EWRHI第44卷第10 期潮间带地热田水化学特征及成因机制分析刘润方,李理（长江岩土工程有限公司，湖北武汉430 0 10)摘要：广东省江门台山汶村潮间带地热田地热资源储量丰富，研究其水化学特点、成因机制和热储温度可为后续开发利用提供参考。通过Piper三线图分析地热水水化学特点，根据代表组分离子比例系数和其日变化、年变化规律分析地热水成因机制；通过氢氧同位素分析并定量估算丰水期与枯水期地热水中咸水淡水混合比例；利用Na-K-Mg三角图和地热温标公式法估算地热田热储温度。结果表明：研究区地热水属海相沉积地热水，热储温度约150，为深部热源沿断层上升过程中受浅表淡水和海水混合作用的结果，丰水期和枯水期的淡水比例分别为7 7.6 0%和58.7 6%。关键词：水化学特征；潮间带；地热田；淡水混合比例；热储温度中图法分类号：P641文章编号：10 0 6-0 0 8 1(2 0 2 3)10-0 10 4-0 60 引 言广东省地热资源丰富，近年众多学者集中对地热田开展了水化学特点和成因机制等研究工作1-3。国内外目前的研究方法主要有水化学特征元素法、同位素法和特征元素比例系数法等4-1。张光辉等利用同位素手段对黑河流域地下水的补、径、排的时空演变特征进行了研究8 。冯建伟等利用地下水水化学特征和代表性矿物离子演化特性对西非凯凯水电站地下水水化学特征及演化趋势研究 。李鹭运用水化学特征系数法对江西省部分温泉的起源、热水年龄、水化学特征等方面进行了研究10 。高胜强等通过针对 Na,Cl,Ca,Mg等特征元素不同组合的比例系数分析和 7 Sr,Sr同位素含量及其比例系数的分析,得到了明月山地区地热水的物质来源和后期演化特点。李静荣等通过分析D,3o,C,sC和 HeHe等多种同位素，全面认识了碳酸泉形成过程12 前人研究普遍认为，广东地区地热系统是在正常或略为偏高的区域热背景下形成的中低温地热系统13，即地下水通过发育的断层深循环对流传热14-15，而沿海地热也可能源于海水。收稿日期:2 0 2 3-0 3-0 1作者简芥：刘润方，男，工程师，硕士，主要从事水利水电工程地质勘察、地热资源研究工作。E-mail:.:104:文献标志码：AD0I:10.15974/ki.slsdkb.2023.10.018一般前人针对地热水的研究，仅通过水化学特点定性分析了地热水成因机制17-19,对地热资源开发规划的参考价值有限。广东省江门台山汶村地热田钻孔热水日涌量达30 0 50 0 t/d，储量丰富，地热资源潜力巨大，但仍尚未被开发利用。本文在定性分析地热水水化学特点和成因机制的基础上，利用氢氧同位素2 0 1定量计算了汶村地热水中不同季节淡水、咸水的混合比例,并估算了研究区地热田热储温度，为各行业合理开发利用汶村地热资源提供重要依据。1研究区概况1.1自然概况研究区位于汶村镇南部镇海弯及海晏之间的海漫滩,东邻珠海特区，北靠江门,西连开平、恩平、阳江三市,南临南海。研究区属于亚热带海洋性气候，温湿多雨，水系发育，河流呈树枝状，具有明显的夏雨型和暴涨暴落特征。区域年最高气温36.9，最低气温-0.1；年平均降雨量2 297.3 mm,降雨不均。1.2地质概况研究区位于新华夏构造带第二隆起带的南部，刘润方等潮间带地热田水化学特征及成因机制分析区域上发育北西向、北东向断裂，出露现代海相沉积（Q ）层,分布于滨海平原和海漫滩，岩性为黄褐、浅灰色黏土、淤泥和砂等。研究区北部为山区,岩性为燕山期花岗岩，见图1。镇海Y5湾Qen南海4kmQ海相沉积层0图1区域地质示意Fig.1 Regional geological overview区内的地下水总体自北部山区汇流补给，沿山区至平原区为径流，最终向南部滨海区排泄。研究区发育近南北走向断裂，是地下水的深部循环通道，为地热水的形成、运移和出露提供水文地质条件。2地下水水化学特征2.1水样分布与采集本文针对研究区地热钻孔开展多次水样采集工作，并进行室内全分析试验，钻孔位置分布见图2。2.2化学组分含量特征根据室内试验成果,研究区地下水化学组分含量见表1。Tab.1Statistics of chemical composition content of groundwater电导率/TNa*/指标H(ms cml)(mg:L-1)(mgL-l)(mgL-1)(mgL-1)(mg-L-l)(mg.L-l)(mgL-1)(mg-L-1)SZ-01(ZK1)7.35SZ-02(ZK2)7.71SZ-03(ZK3)7.60SZ-04(ZK4)7.74SZ-05(雨水)7.35SZ-06(海水)7.84XZ(沟水)6.75268.91NQm海相沉积层ZK1ZK1钻孔位置及编号Qem海宴镇O研究区Y5燕山期花岗岩表1地下水水化学组分含量统计K+/Ca2*/Mg?+/Si4+/S2+/Li+/14.2081.421164.2514.7282.351168.7114.8079.831082.5114.1281.231098.9315.3425.0025.9727.006575.0025.0010.37ZK2L二小路ZK3Q.m?ZK4海漫滩Y5图2 钻孔位置分布示意Fig.2Borehole location distribution map研究区地热水、雨水和海水均呈弱碱性，溪沟水呈弱酸性。依据水化学组分含量绘制Piper三线图，见图3。研究区钻孔地热水出露温度在7983,属中低温热水，水样均落在7 区,说明水样碱金属大于碱土金属，强酸根大于弱酸根,非碳酸盐碱 50%。水化学类型均为Na-Cl型水,尤其是C1-，占阴离子的95%以上，受海水影响大。淡水一般为CaH CO，型水或Ca-Cl型水。研究区雨水及海水NO；含量最大为1.39mg/L,地表沟水NO；含量达2 5.7 2 mg/L。地热水NO；含量为3.52 9.96 mg/L，在淡水、降雨及海水之间，表明地热水受到地表沟水等的混合作用。地热水电导率达14ms/cm以上，一方面由于热水由深部地热流体通过断裂带上升而来；另一方面也由于其受海水侵袭的影响较大。地热水中Ca+含量达6 0 0 8 0 0 mg/L,远高于海水、雨水和地下淡水中Ca*含量。与此同时，地热水中Mg?*，Na*含量远低于海水。这是因为地热水先与海水混合，使Mg+，Na*含量增大，后又与岩层经历阳离子交换吸附的作用,热水中的Mg+,Na*被岩层中Ca+交换，导致热水中的Mg+，Na*的含量低于海水,Ca+含量高于海水。CI-/SO2-/(mg:L-1)(mgL-l)116.16643.99134.09750.77142.65826.32117.11659.971.740.42247.404.90Br/HCO;/(mg:L-l)（(mgL-1)(mg:L-1)7.9148.529.5756.3412.7158.928.4554.883.160.14260.70756.0017.532.75F-/NO;/17.986.1818.836.0017.303.4818.756.829.154.936.760.085742.164920.905890.495260.202.871.52914245.910.00817.39226.27207.14276.40199.481.901667.788.4330.0022.0323.0526.54121.270.021877.2277.7985.7379.4321.59121.0165.103.183.204.603.310.18:105.3.744.069.963.521.3925.722023年10 月20%钠十钅镁(Na+80%20%60%40%40%60%20%480%80%60%40%20%钙(Ca)图3地下水Piper三线图Fig.3Piper three-line diagram of groundwater地热水中K*,SiO2含量远高于海水,这是由于研究区地层为花岗岩，其所含主要矿物为钾长石（K,0 Al,0，6Si0,）和石英，含量分别达2 5%58%、8%2 8%。高温流体流经花岗岩岩层时溶滤含钾硅酸盐或石英，使K*,SiO，含量增高。花岗岩Sr2+,F-含量较一般沉积岩高，所以地热水中 Sr+，F-的含量也有同样特点。综上,研究区地热水为Na-Cl型水,受海水、降雨和地表淡水等的混合作用，并在径流过程中经历水一岩相互作用。3地热水成因机制和混合作用分析3.1地热水代表组分离子比例系数分析根据水化学试验成果绘制不同组分间的离子比例系数,见图4。CI/Br3001.0yMg/CaYNa/YCI2500.8200-0.6150-0.41000.2注：水平线表示一般海水或与海水相关的沉积水中的离子比例系数。图4地热水水化学离子比例系数分布Fig.4Distribution diagram of chemical ion proportioncoefficient of geothermal water(1）C I/Br 系数。CI/Br为地下水中Cl,Br质量浓度比值,是判定地下水盐分来源的重要依据,尤其:106水利水电快报EWRHI是在区分陆相蒸发源和海相源方面2 1。大洋水中的Cl/Br一般约为30 0 2 2 。淡水的Br含量低，017.80%80%+60%60%C40%40%20%20%40%60%80%20%40%60%80%氯(CI)60,6555504453402Ca/Sr35130上取样点第44卷第10 期7-0203040506010203-04ZX艺-0 580%60%硫酸40%(SO.)20%0Cl/Br一般大于8 8 0。由图4可知,本区地热水Cl/Br平均值为2 2 2，接近30 0。这说明地热水很可能是由海相沉积水组成,残余海水在长期地质历史过程中由于浓缩而产生NaCl沉淀，而溴化物溶解度更大，故残余海水或海相沉积水Cl/Br300。(2）Na/c i 系数。Na/ci为地下水中Na*,Cl-毫克当量浓度比值。海相沉积水中Na/c值一般小于0.8 52 3。本区地热水Nac 值平均为0.34，小于0.8 5,与海相沉积水特点相符。(3）Ca/Sr 系数。Ca/Sr为地下水中 Ca,Sr?+质量浓度比值。当海水浓缩产生盐类沉淀时,SrSO4的沉淀在CaCO，之后,故海水中Ca/Sr系数较小,盐类浓缩沉淀后更小，约为33。本区地热水Ca/Sr系数为34 36。因此可得知汶村地热水是与海水沉积有关的地下热水。（4）M g/c a 系数。Mg/ca为地下水中Mg,Ca毫克当量浓度比值。海水中Mg含量一般比Ca大，Mg/ca系数约为5.5;淡水中Mg含量一般远比Ca小，Mgc a 系数远小于1。因此，Mg/ca系数常被用做计算判断海水的人侵范围和程度。本区地热水mgc.系数为0.15 0.2 6,平均为0.2,小于0.55。这是由于淡水混合的作用，同时也是由于水岩相互作用导致热水中Mg含量降低、Ca含量增加的结果。3.2地质条件分析区域内断裂发育,对热水出露有显著的控制作用：断层及由断层形成岩体破碎带均是地热水的储存空间；热源可通过深大断裂从地球深部往地表运移。因此，地热温泉多沿断裂带呈条带状分布。研究区发育北西向断层f1，包含各次级断层,其组合构成了汶村地热田的储热构造和运移通道。综上，汶村地热水属海相沉积水，与深部热源的混合加热后在上升过程中经历水岩相互作用,并受到地下淡水和海水混合共同作用，主要补给来源为降雨，主要补给区为北部及东北部的山区，见图5。3.3地热水混合作用分析Cl-是地下水中最稳定的元素，常用来分析地下水的混合作用及程度。本文对2 4h和1个水文年内的钻孔地热水取样分析，CI含量日变化、年变化分别见图6,7。刘润方等潮间带地热田水化学特征及成因机制分析400南降雨海研究区40-80Q海相沉积层燕山期花岗岩口推测断层国裂隙网络地下水流向图5汶村地热水形成机制概化模型Fig.5 Schematic model of formation mechanism ofWencun geothermal water700065006000(,T.3u)/鲁号,105500500045004000350030000图6 地热水CI-含量随潮汐日变化动态Fig.6IDiurnal variation of geothermal water Cl-with tide7000110006500CI(.T.3ul)/吾&,ID60005500500045004000关关系，随降雨量的增加而减少,表明了地热水有淡水的混合作用。Q3.4地热水咸、淡水混合比例计算ZK4地热水2.55.0510时间/h降雨60-1057.510.0距离/kmCI潮汐15D和18 0 为水中较为稳定的同位素,对于确定地12.515.0543220258006004002000下水成因类型等具有重要意义。根据广州大气降水方程线，研究区地热水属大气降水成因2 4-2 5。地热水的淡水来源主要为降雨，咸水来源为海水、古海水。分别将淡水与咸水作为混合的两个端元,地热水为两个端元混合的结果，利用D、10 计算不同端元的混合比例，见表2。表2 汶村地热田热水氢氧同位素含量Tab.2Hydrogen and oxygen isotope content of hot water inWencun geothermal field%o9月2月钻孔8180SZ-02-5.81SZ-03-5.94SZ-04-5.97海水4.03降雨-5.88注：SZ-0204为钻孔热水；降雨a180、8 D 为广州降水多年平均值2 5。由于18 0 值受混合作用以外的其他因素影响程度高于SD,本文取&D值作为计算混合比例的指标，计算公式为D(淡水)X+D(海水)(1X）=D(地热水)式中：X为混人淡水的比例,%；8 D（淡水）取值-40.93%o;8D（海水)为研究区海水中的8 D值,取值-2 3.7 4%o;8D（地热水)为研究区地热水中的8 D值。分别取2 月、9月水样氢氧同位素数据进行咸、淡水混合比例估算,结果见表3。由计算结果可知，8D-36.1136.8637.0823.7440.938180-4.99-5.31-5.42SD-33.8435.2634.872淡水混合比例高于海水。其中,9 月丰水期的淡水豆图7 地热水CI含量年变化动态Fig.7Dynamic chart of Cl-content of geothermalwaterwith year由图6,7 可知,研究区地热水CI-含量与潮汐呈正相关关系。在潮汐减小时,Cl-含量减小明显，但潮汐增大时，Cl-含量增大相对不明显。这是由于海水Cl-含量为142 45.9 1mg/L,降雨与地下冷水分别为2.8 7 mg/L和17.39 mg/L，地热水为5500mg/L左右,相比海水的咸化作用,淡水稀释作用更明显,因此在潮汐增强时,Cl-含量变化相对不明显。Cl-含量与降雨量在年内呈较明显的负相三EI-1O日期/月-日)混合比例达7 7.6 0%,2 月枯水期的淡水混合比例降至58.7 6%。4地热田热储温度估算热储温度一般难以直接测量，常用地热温标法估算，主要分两大类：SiO2地热温标，适用于150以上热水；玉髓温标适用于温度低于150 的低温热水。阳离子地热温标，适用于水岩完全平衡的地热田,常见的有Na-K温标、K-Mg温标等,各温标公式见表4。目前,国内外常用Na-K-Mg三角图判断水-岩平衡样点,见图8。由图8 可知，研究区地热水处于部分平衡区，更适宜采用石英温标或玉髓温标估算。:107.2023年10 月表3汶村潮间带地热水淡水与海水混合比例Tab.3 Mixing ratio of geothermal water,fresh water andsea water in Wencun intertidal zone钻孔淡水SZ-0271.96SZ-0376.32SZ-0477.60淡水平均混合比例77.60表4地热温标公式Tab.4Geothermal temperature scale formula地热温标类型编号石英温标T1T=1309/(5.205-lgSi0,)273.15玉髓温标T2T=1032/(4.69-1gSi02)-273.15钠-钾温标T3T=1390/1g(Na/K)+1.75-273.15钾-镁温标T4T=4410/14-lg(K/Mg)-273.15注：Na,K,Mg,SiO2代表水中对应离子或化合物的质量浓度，单位为mg/L;温度单位为。淡水热水240260280300钾(K)/100图8 地下水Na-K-Mg三角图Fig.8 Na-K-Mg triangle diagram of groundwater根据不同钻孔水样数据计算得到热储温度，见表5。结果表明：玉髓温度本更适合计算热储温度低于150 的低温热储，但计算结果达17 0，明显偏大。石英热储温度分布均匀，均在150 左右。因此,综合Na-K-Mg三角图可以发现,由石英温标计算得到的热储温度更准确，平均热储温度为154.62 。表5热储温度估算Tab.5Thermal storage temperature estimation钻孔石英温标SZ-01153.29SZ 02152.73SZ-03156.09SZ-04156.375 结 论（1）研究区地热水水化学类型为Na-Cl型,受海水影响大;淡水水化学类型一般为Ca-HCO，型、:108.水利水电快报EWRHICa-Cl 型。（2）受潮汐影响，研究区地热水Cl-含量与潮%汐呈一定的正相关关系;受降雨影响,Cl-含量与降9月2月海水淡水28.0458.7623.6867.0222.4064.7558.76公式钠(Na)/1000完全平衡区10Q220未平衡区第44卷第10 期海水雨量呈较明显的负相关关系。41.24（3）研究区地热水为海相沉积水，与深部热源32.98混合加热后在上升过程中经历水岩相互作用,并35.25受到地下淡水和咸水共同混合作用。（4）研究区地热水为淡水、咸水混合结果，其中，丰水期淡水混合比例为7 7.6 0%，枯水期淡水混合比例为58.7 6%。（5）根据石英温标法，估算汶村潮间带地热田热储温度为154.6 2。参考文献：1赵慈平.腾冲火山区现代慢源氮释放特征及深部岩浆活动研究D.北京：中国地震局地质研究所,2 0 10.2GUO Q.Hydrogeochemistry of high-temperature geother-mal systems in China:A reviewJ.Applied Geochemis-try,2012,27(10):1887-1898.3 GUO Q,WANG Y,LIU W.Hydrogeochemistry and envi-ronmental impact of geothermal waters from Yangyi of Ti-bet,China J.Journal of Volcanology and GeothermalResearch,2009,180(1):9-20.4GUO Q,WANG Y.Geochemistry of hot springs in theTengchong hydrothermal areas,Southwestern China J.镁(Mg)001/2玉髓温标169.90170.50170.61170.68Journal of Volcanology and Geothermal Research,2012,215:61 73.5GUO Q,WANG Y,LIU W.Major hydrogeochemicalprocesses in the two reservoirs of the Yangbajing geother-mal field,Tibet,China J.Journal of Volcanology andGeothermal Research,2007,166(3):255-268.6刘虹，张国平，金志升，等云南腾冲地区地热流体的地球化学特征J.矿物学报,2 0 0 9（4）：496-50 1.7朱贵麟，张林友，李旭峰，等.青海共和盆地地下水演化特征及成因机制J.人民长江，2 0 2 2,53（8）：36 42.8 张光辉，费宇红，聂振龙，等.区域地下水演化与评价理论方法M.北京：科学出版社,2 0 14.9冯建伟，李林，李辰舟.西非凯凯水电站地下水水化学特征及演化趋势研究J.水利水电快报，2 0 2 2,43(6):125-131.10李鹭.江西省宜春地区温泉水文地球化学特征研究D.南昌：东华理工大学,2 0 18.11 高胜强，刘凯，孙文洁，等.花岗岩带状热储水化学与锶同位素地球化学特征J.人民长江，2 0 2 2,53（8）：36 42,49.【12 李静荣,王中正,王亚，等.广东省河源断裂带碳酸泉刘润方等潮间带地热田水化学特征及成因机制分析水化学特征及形成机制J.中山大学学报（自然科学版),2 0 18,57(5):19-2 8.13王东升，王经兰.中国地下热水的基本类型和成因特征J.第四纪研究,19 9 6,16(2)：139-146.14张珂,马浩明，蔡剑波.华南沿海温泉成因探讨J.中山大学学报（自然科学版）,2 0 0 2,41（1）：8 2 8 6.15CAPACCIONI B,VASELLI O,TASSI F,et al.Hydrogeo-chemistry of the thermal waters from the Sciacca Geo-thermal Field(Sicily,southern Italy)J.Journal of Hy-drology,2011,396(3-4):292-301.16AWALEH M O,HOCH F B,KADIEH I H,et al.The ge-othermal resources of the Republic of Djibouti-I:Hydrogeochemistry of the Obock coastal hot springs JJournal of Geochemical Exploration,2015,152:54-66.17李皓婷.水化学方法在地热资源勘查中的应用综述J.农业与技术,2 0 2 0,40(10)：145-146.18邓亚东，陈伟海，罗书文，等.重庆巴南区热洞温泉成因探讨J.人民长江,2 0 13,19（15）：56-58.19于晓红，赵宾，邹凤钗，等.贵州省铜仁地区地下热水的水化学特征研究J.吉林化工学院学报,2 0 11,2 8Anlysis of hydrochemical characteristics and genetic mechanism(5):86-90.20 LEYBOURNE M I,GOODFELLOW W D.Br/Cl ratiosand O,H,C,and Bisotopic constraints on the origin ofsa-line waters from eastern Canada J.Geochimica et Cos-mochimica Acta,2007,71(9):2209-2223.21张人权,梁杏，靳孟贵，等.水文地质学基础M.武汉：中国地质大学出版社，2 0 2 2.22高胜强，刘凯，孙文洁.花岗岩带状热储水化学与锶同位素地球化学特征J.人民长江，2 0 2 2,53（8）：3642.23 王现国，张慧，张娟娟.开封凹陷区地热水水化学特征及同位素分析J.安全与环境工程,2 0 12,19（6：88 92.24 林韵，高磊，李绍恒，等.广东江门地热水水文地球化学特征及来源分析J.环境化学，2 0 2 0,39（2）：512523.25 刘进达，赵迎昌，刘恩凯，等.中国大气降水稳定同位素时一空分布规律探讨J.勘察科学技术,1997,3(7):34-39.(编辑:李慧)of intertidal geothermal fieldLIU Runfang,LI Li(Changjiang Geotechnical Engineering Company Limited,Wuhan 430010,China)Abstract:The intertidal geothermal field in Wencun,Taishan,Jiangmen,Guangdong Province has abundant re-serves of geothermal resources.The study of its hydrochemical characteristics,genetic mechanism and thermal stor-age temperature can provide a reference for subsequent development and utilization.The hydrochemical characteris-tics of geothermal water were analyzed by using Pipers three-line diagram,and the genesis mechanism of geother-mal water was studied according to the proportion coefficient of representative component ions and its daily and an-nual changes.The mixture ratio of salt water-fresh water in geothermal water of wet and dry seasons was quantita-tively estimated by innovative hydrogen and oxygen isotopes.The thermal storage temperature of geothermal field wasestimated by using Na-K-Mg trigonometric diagram and geothermal temperature scale formula.The results showedthat the geothermal water in the study area belonged to marine sedimentary geothermal water,and the thermal stor-age temperature was about 150 C.These attributed to mix of shallow fresh water and seawater during the risingprocess of deep heat source of the geothermal field along the fault,and the proportion of fresh water in the wet sea-son and dry season was 77.60%and 58.76%,respectively.Key words:hydrochemical characteristics;inter-tidal zone;geothermal field;mixed ratio of fresh water;thermalstorage temperature109.