方茜娟，劉運(yùn)濤，左 偉，李 川，周稱稱

1.河南省國(guó)土資源開發(fā)投資管理中心，河南 鄭州 450046；2.河南省地質(zhì)礦產(chǎn)勘查開發(fā)局 第五地質(zhì)勘查院，河南 鄭州 450001；3.河南省水文地質(zhì)工程地質(zhì)應(yīng)用技術(shù)中心，河南 鄭州 450001

0 引言

中國(guó)地?zé)豳x存類型分為沉積盆地型和隆起山地型，后者可分為火山型和構(gòu)造型[1]。沉積盆地型地?zé)崃黧w化學(xué)、成因研究程度較高[2-6]，地?zé)崃黧w的水化學(xué)類型主要為Cl-Na、Cl·SO4-Na型等，主要補(bǔ)給來(lái)源為大氣降水，通過(guò)C14測(cè)年顯示形成年齡20～50 ka[4-7]?；鹕叫偷?zé)崽镏饕栽颇向v沖和臺(tái)灣大屯地區(qū)為代表，其具有強(qiáng)烈的地?zé)犸@示，如沸泉、酸性泉、水熱爆炸和硅華等[8]；構(gòu)造型地?zé)嵩谖覈?guó)中西部山區(qū)以及東南沿海分布較廣，水化學(xué)類型多為HCO3·SO4-Na、HCO3-Na型，如河南魯山[9]、陜西眉縣[10]、河北遵化[11]、四川巴塘[12]、廣東惠州[13]等，3H、C14測(cè)年結(jié)果為0.01～0.5 ka[14-16]。這些研究工作解釋了地?zé)崽锪黧w化學(xué)特征和來(lái)源(大氣降水等)，但對(duì)帶狀熱儲(chǔ)補(bǔ)給源與徑流分析不足。此次研究利用流體化學(xué)、區(qū)域地?zé)岬刭|(zhì)構(gòu)造構(gòu)建補(bǔ)給、徑流、排泄全過(guò)程，進(jìn)而系統(tǒng)分析帶狀熱儲(chǔ)地?zé)崽锏某梢驒C(jī)制。

1 地?zé)岬刭|(zhì)條件

1.1 區(qū)域溫泉分布

新安—平頂山斷裂帶為區(qū)域控?zé)針?gòu)造，走向130°～140°，傾向北東，形成于燕山期，喜山期復(fù)活，至今仍在活動(dòng)，沿線分布溫泉泉群3處(圖1，引自河南省基巖地質(zhì)圖)，分別是洛陽(yáng)新安縣暖泉溝溫泉群(2處，29℃)、洛南龍門山溫泉群(19處，24℃～43℃，通過(guò)草店、魏灣斷層與西南平頂山斷裂交匯)，汝州溫泉鎮(zhèn)溫泉群。

1.2 研究區(qū)地?zé)岣艣r與斷裂構(gòu)造

溫泉鎮(zhèn)地?zé)崽锸堑湫偷臉?gòu)造型地?zé)崽?。溫泉呈泉群排泄，大體北西南東向出露，共10個(gè)泉點(diǎn)，泉流量8.294～14.82 m3/h，溫度48℃～71℃?，F(xiàn)有地?zé)峋?眼，開采利用7眼。研究區(qū)主要斷裂構(gòu)造有：

(1)北西向斷裂：新安—平頂山斷裂帶，上盤為中元古界熊耳群許山組、汝陽(yáng)群云夢(mèng)山組；下盤為寒武系張夏組、石炭系、二疊系。斷裂帶寬30～50 m，構(gòu)造碎裂巖、斷層泥發(fā)育；斷面產(chǎn)狀40°∠50°～70°，斷距>1 000 m，正斷層。主要有F2、F3、F4、F5四個(gè)次級(jí)斷層。

(2)近東西向斷裂：九皋山—溫泉街?jǐn)嗔褞В蠪0、F1兩個(gè)次級(jí)斷層，F(xiàn)1斷層走向南東，為正斷裂，傾向北東，傾角約80°，上盤為中元古界安山玢巖，下盤為寒武系灰?guī)r。

(3)北東向斷裂：F6斷層、F7斷層，均為北東走向，傾向南東，正斷層，寬度約10 m，前者傾角約80°、后者傾角約70°，形成晚于F1斷層、F2斷層、F3斷層。

2 樣品采集與測(cè)試

2018年8月采集同位素分析測(cè)試樣品7組(圖2)，其中地表水樣品3組、地?zé)崃黧w樣品1組(ZK1井)、淺層地下水樣品1組(<20 m冷水)、深層地下水樣品1組(300 m冷水)，大氣降水樣品1組；同期采集測(cè)試地?zé)崃黧w全分析樣品2組，對(duì)應(yīng)WR04、WR08井。

圖2 研究區(qū)構(gòu)造與采樣點(diǎn)分布圖Fig.2 Geological structural and sample point distributionmap of the study area

硫氧同位素前處理方法：現(xiàn)場(chǎng)分裝1.5 L水樣，24 h內(nèi)用孔徑0.45 μm醋酸纖維濾膜過(guò)濾，加入1 mL優(yōu)級(jí)純濃鹽酸，混合均勻后，靜置0.5 h后，加入10 mL飽和氯化鋇溶液，靜置過(guò)夜；實(shí)驗(yàn)室內(nèi)用0.22 μm孔徑PES濾膜過(guò)濾生成的硫酸鋇沉淀，置于850℃馬弗爐內(nèi)烘干2 h，待測(cè)。測(cè)試方法：Flash2000HT元素分析儀和Delta V穩(wěn)定同位素質(zhì)譜儀測(cè)定，采用國(guó)際標(biāo)準(zhǔn)NBS127(硫同位素值為20.3‰，氧同位素值為8.6‰)和實(shí)驗(yàn)室標(biāo)準(zhǔn)校準(zhǔn)樣品，其中硫同位素分析精度小于0.2‰，氧同位素分析精度小于0.5‰，測(cè)試在中國(guó)科學(xué)院南京地質(zhì)古生物研究所和河南理工大學(xué)生物遺跡與成礦過(guò)程省級(jí)重點(diǎn)實(shí)驗(yàn)室完成。水樣分析結(jié)果見表1。

表1 水化學(xué)與穩(wěn)定同位素分析結(jié)果

3 水化學(xué)與同位素特征

3.1 水化學(xué)組成特征

地?zé)崃黧w水化學(xué)類型HCO3·SO4·Cl-Na·Ca型，三組樣品SO4、HCO3、Cl摩爾當(dāng)量較為接近，陽(yáng)離子以Na+為主(圖3)。在指示理療熱礦水相關(guān)因子方面，F(xiàn)、H2SiO3含量較其它水體明顯偏高。同處新安—平頂山斷裂帶的暖泉溝、龍門地?zé)崽?，地?zé)崃黧w水化學(xué)類型SO4·Cl-Na·Ca，與研究區(qū)水化學(xué)類型和組分基本相似。

圖3 Piper三線圖Fig 3 Piper map of water samples in studied area

相對(duì)地?zé)崃黧w，非封閉地表水體礦化度較低，汝河水(SY-1)水化學(xué)類型為HCO3-Ca型，陽(yáng)離子以Ca2+為主，陰離子以HCO3-為主，體現(xiàn)了碳酸鹽巖地區(qū)地表水的特征；澗山口水庫(kù)(SY-2)水化學(xué)類型為HCO3·SO4·Cl-Na·Ca型，指示了水域的半封閉性；坑塘(SY-3)礦化度及主要陰陽(yáng)離子均遠(yuǎn)大于其它水體，具封閉、人為影響強(qiáng)烈的水體特征。地下水(冷水)TDS值介于地表水和地?zé)崃黧w之間，主要陰陽(yáng)離子為HCO3-、Ca2+，水化學(xué)類型HCO3·Cl-Na·Ca型，體現(xiàn)了淋濾水的特征?？偟恼f(shuō)來(lái)，較其它水體，地?zé)崃黧w是高礦化度、富鈉、貧鎂、硫酸鹽豐富的氟·偏硅酸型理療熱礦水。

3.2 同位素組成特征

研究區(qū)不同水體氫氧、硫酸鹽硫氧同位素組成差異明顯(見表1)，但地?zé)崃黧w與同處新安—平頂山斷裂帶的安暖泉溝、洛陽(yáng)煤田二隊(duì)(洛陽(yáng)龍門)地?zé)崃黧w氫氧同位素較為接近。2013年12月，河南省地質(zhì)調(diào)查院測(cè)試值分別為-67‰和-9.2‰、-66‰和-9.2‰、-67‰和-8.9‰，可暫時(shí)推斷三者補(bǔ)給為同一來(lái)源。

4 討論

4.1 地?zé)崃黧w來(lái)源

4.1.1 同位素差異

研究區(qū)水體氫氧同位素組成(圖4)顯示其來(lái)源及所經(jīng)歷蒸發(fā)過(guò)程的差異性。各水體氫氧同位素位于鄭州西南大氣降水線附近[17]，表明其主要受大氣降水補(bǔ)給。地?zé)崃黧wSY-4、深層地下水SY-6存在δ18O值偏負(fù)的漂移，表明與圍巖中氧同位素發(fā)生同位素交換。地表水SY1、SY2、SY3隨封閉程度增加(河流—水庫(kù)—坑塘)蒸發(fā)過(guò)程愈強(qiáng)，淺層地下水SY-5亦顯示經(jīng)歷蒸發(fā)過(guò)程。

圖4 研究區(qū)不同水體氫氧同位素組成圖Fig.4 Water isotope compositions in studied area

圖5 溶解性硫酸鹽硫和氧同位素組成Fig.5 Composition of dissolved sulfate sulfur and oxygen isotopes

不同類型水體硫酸鹽硫和氧同位素的差異同樣顯著(表1，圖5)，大氣降水SY7硫酸鹽含量低，對(duì)地表水SY1、SY2、SY3、淺層地下水SY5硫酸鹽貢獻(xiàn)有限，SY1、SY2、SY5硫和氧同位素值接近，主要來(lái)自于土壤有機(jī)硫氧化以及人類活動(dòng)影響，坑塘地表水SY3處于封閉環(huán)境，存在硫酸鹽細(xì)菌還原過(guò)程，導(dǎo)致其硫同位素異常偏高。深層地下水SY-6硫和氧同位素值較低，說(shuō)明硫酸鹽來(lái)源受硫化物氧化控制比例較大。地?zé)崃黧w與其它水體硫氧同位素的差異顯示了補(bǔ)給不是來(lái)源于區(qū)內(nèi)其它水體，是遠(yuǎn)程補(bǔ)給的。

4.1.2 硫同位素來(lái)源

地?zé)崃黧w中硫同位素值較高，而氧同位素值較低，推測(cè)原因可能包括硫酸鹽細(xì)菌還原作用黃鐵礦及石膏溶解作用。硫酸鹽細(xì)菌還原作用的基本反應(yīng)方程式如下：

(1)

溫泉地?zé)崽餆醿?chǔ)主要是石炭系、寒武系石灰?guī)r，常伴生石膏、黃鐵礦等礦物，其溶解會(huì)改變地?zé)崃黧w中硫酸鹽硫同位素構(gòu)成。奧陶系、寒武系石膏中硫同位素通常20‰～32‰，煤系地層(石炭系)黃鐵礦硫同位素多低于4‰[18]，故分析推斷硫更多來(lái)源于石膏溶解。地?zé)崃黧w中硫酸鹽氧同位素值降低可能與周圍水體的氧同位素平衡交換作用有關(guān)，在高溫環(huán)境中地?zé)崃黧w硫酸鹽氧同位素與水氧同位素發(fā)生同位素交換，造成流體硫酸鹽氧同位素值降低[19]。

4.2 地?zé)崽锍梢蚍治?/h3>

4.2.1 補(bǔ)給高程

利用大氣降水δD和δ18O隨高程增加呈線性降低規(guī)律[20]，可以計(jì)算補(bǔ)給高程，通過(guò)對(duì)水文地質(zhì)單元地質(zhì)構(gòu)造、補(bǔ)徑排條件分析，進(jìn)而判斷補(bǔ)給區(qū)位置。計(jì)算公式：

H=H0+(δD-δD1)/g

(2)

式(2)中，H地?zé)崃黧w補(bǔ)給高程(m)；H0-樣品采集高程(m)；δD-補(bǔ)給水D同位素含量(‰，VSMOW)；δD1送檢樣品D同位素體積比(‰，VSMOW)；g為δD隨高程遞減梯度(‰/100 m)。

以洛陽(yáng)市大氣降水δD(‰，VSMOW)均值-52.87[21-22]作為補(bǔ)給水的計(jì)算值。洛陽(yáng)市δD高程遞減梯度介于1.35‰～2.35‰/100 m，計(jì)算補(bǔ)給高程791～566 m。

表2 補(bǔ)給高程計(jì)算結(jié)果表

4.2.2 溫標(biāo)計(jì)算

常用的地?zé)釡貥?biāo)有陽(yáng)離子溫標(biāo)和二氧化硅溫標(biāo)，使用的前提是作為溫標(biāo)計(jì)算的溶質(zhì)與圍巖達(dá)到了溶解-平衡的狀態(tài)。汝州市地?zé)崃黧w屬未成熟水，不宜使用與Na、K、Mg相關(guān)的Na-K、K-Mg等陽(yáng)離子溫標(biāo)。常用的二氧化硅溫標(biāo)為石英溫標(biāo)、玉髓溫標(biāo)，這類溫標(biāo)計(jì)算公式是建立在熱水中SiO2溶解度與溫度的關(guān)系之上的。石英地?zé)釡貥?biāo)適用于熱交換溫度120～180℃的條件，玉髓應(yīng)用于<110℃的條件[23]。

1.石英溫標(biāo)(傳導(dǎo)冷卻)[24]：

(3)

2.石英溫標(biāo)(最大蒸汽損失)[25]：

(4)

3.玉髓溫標(biāo)[24]：

(5)

4.計(jì)算結(jié)果分析

溫標(biāo)計(jì)算結(jié)果見表3。前文已分析洛陽(yáng)龍門與研究區(qū)地?zé)崃黧w水化學(xué)、氫氧同位素是相似的，推斷可能是同源的，而洛陽(yáng)龍門實(shí)測(cè)最高溫度97.2℃，遠(yuǎn)高于玉髓計(jì)算值，故選擇二氧化硅溫標(biāo)計(jì)算值。

表3 溫標(biāo)計(jì)算結(jié)果Table 3 Results of temperature scale calculation

4.2.3 水熱循環(huán)深度推斷

新安縣—平頂山斷裂帶洛陽(yáng)段，石臺(tái)階組平均地溫梯度3.40℃/100 m，蟒川組3.57℃/100 m，三疊系3.88℃/100 m，石炭、奧陶、寒武約4℃/100 m，自上而下地溫梯度逐呈漸增高趨勢(shì)。常溫層溫度14℃，推斷水熱循環(huán)深度2 812～3 109 m。

4.2.4 成因推斷

新安縣—平頂山斷裂帶體系中，僅新安縣曹村鄉(xiāng)西北低山區(qū)高程(500～800 m)與計(jì)算補(bǔ)給高程值相符，可初步推斷為補(bǔ)給水來(lái)源，盡管其地層為汝陽(yáng)群、洛峪群石英砂巖、礫巖，不利于大氣降水下滲進(jìn)入地下水系統(tǒng)，但通過(guò)地表徑流形式快速進(jìn)入下游的曹村鄉(xiāng)、石井鎮(zhèn)、北冶鎮(zhèn)、石寺鎮(zhèn)(350～500 m)一帶(圖6)，這些地區(qū)以O(shè)+∈碳酸鹽巖為主，溶蝕較強(qiáng)，且處于兩大構(gòu)造交匯處，裂隙發(fā)育，上游徑流易進(jìn)入地下水系統(tǒng)。入滲后向南東運(yùn)移增溫，在新安縣城區(qū)一帶局部受阻上涌，出露暖泉溝溫泉(29℃)，其余大部運(yùn)移至洛陽(yáng)盆地西緣洛河褶斷束，此處三疊系(砂巖)、古近系(蟒川組、石臺(tái)階組致密砂巖)沉積厚度約2 000 m，地?zé)崃黧w受阻，向下進(jìn)入下部巖溶發(fā)育的碳酸鹽巖層(埋藏深度約3 000 m)，這與計(jì)算推斷的熱循環(huán)深度是吻合的。穿越洛河褶斷束后，地層抬升，碳酸鹽巖熱儲(chǔ)層埋深1 000～2 000 m，地?zé)崃黧w向上運(yùn)移，此時(shí)溫度與圍巖交換逐漸下降，至鳳翔山莊、龍門山一帶，揭露最高95.2℃。在臨汝鎮(zhèn)一帶，寒武系—石炭系碳酸鹽巖地層再次抬升，流體溫度持續(xù)冷卻，至溫泉鎮(zhèn)遇阻出露，溫度約60℃。

圖6 地?zé)崃黧w成因示意簡(jiǎn)圖Fig.6 diagram of the causes of geothermal fluid

5 結(jié)論

(1)汝州溫泉地?zé)崃黧w水化學(xué)類型SO4·HCO3·Cl-Na·Ca型，是高礦化度、富鈉、貧鎂、硫酸鹽豐富的氟·偏硅酸型復(fù)合理療熱礦水。

(2)氫氧同位素研究表明地?zé)崃黧w主要來(lái)源于大氣降水。溶解性硫酸鹽同位素研究顯示地?zé)崃黧w同其它水體的顯著差異，表明地?zé)崃黧w的補(bǔ)給與研究區(qū)無(wú)關(guān)，來(lái)源于遠(yuǎn)程補(bǔ)給。計(jì)算補(bǔ)給高程為566～791 m，結(jié)合新安—平頂山斷裂帶地形地貌和地質(zhì)條件，確定補(bǔ)給來(lái)源于新安縣石井、石寺盆地，該盆地出露碳酸鹽巖，溶蝕較強(qiáng)，且斷裂發(fā)育，高程雖低于計(jì)算值，但受納北麓低山區(qū)降水徑流。

(3)研究區(qū)地?zé)崃黧w屬未成熟的水，石英溫標(biāo)計(jì)算值最可靠，形成溫度112.49℃～138.36℃，推斷水熱循環(huán)深度為2 812～3 109 m，與新安—平頂山斷裂帶穿越洛陽(yáng)盆地西緣洛河褶斷束的碳酸鹽巖地層埋深吻合。隨斷層上涌，與低溫圍巖物質(zhì)和溫度交換，分別在龍門山、溫泉鎮(zhèn)遇阻上涌，以溫泉形式排泄，形成地?zé)崽铩?/p>