摘要：菏澤斷裂北部鄄城地?zé)崽锏責(zé)豳Y源豐富，探究區(qū)域地?zé)岢梢驒C(jī)制對菏澤斷裂北部地區(qū)的地?zé)豳Y源開發(fā)利用和經(jīng)濟(jì)發(fā)展具有重要意義。本文通過對采集的樣品進(jìn)行水文地球化學(xué)、氣體及其同位素測試分析，結(jié)果顯示，鄄城地?zé)崽锏責(zé)崴瘜W(xué)類型主要以SO4·HCO3·ClCa·Na、SO4·HCO3Ca·Na和SO4Ca·Na型為主；鄄城地?zé)崽锍尸F(xiàn)出從補(bǔ)給區(qū)的淺循環(huán)－開放式巖溶地下水流系統(tǒng)（Ⅰ），至徑流區(qū)的中循環(huán)—半開放式巖溶地?zé)崴飨到y(tǒng)（Ⅱ），再至排泄滯留區(qū)的深循環(huán)—弱開放式巖溶地?zé)崴飨到y(tǒng)（Ⅲ）的演化路徑。熱源以殼源為主，屬典型的中低溫傳導(dǎo)型地?zé)嵯到y(tǒng)。

關(guān)鍵詞：鄄城地?zé)崽?水化學(xué)類型;同位素化學(xué);成熱模式；菏澤地區(qū)

中圖分類號(hào)：P314""" 文獻(xiàn)標(biāo)識(shí)碼：A""" doi：10.12128/j.issn.16726979.2024.03.013

0 引言

地?zé)崾且环N清潔低碳、分布廣泛、資源豐富、安全穩(wěn)定的優(yōu)質(zhì)可再生能源［1］。菏澤市地?zé)豳Y源豐富，以往研究程度較高。馬振民［23］等利用水文地球化學(xué)方法、同位素分析方法及古水文地質(zhì)條件分析方法探討了菏澤凸起地?zé)崴梢?。康鳳新［4］等對菏澤潛凸起巖溶熱儲(chǔ)地?zé)崽锏乃牡厍蚧瘜W(xué)特征、同位素、地溫場、水動(dòng)力場進(jìn)行了系統(tǒng)分析，揭示了地?zé)崴a(bǔ)給、徑流和富集機(jī)理。史啟朋［5］等分析了菏澤凸起地?zé)崽飵r溶地?zé)崴乃瘜W(xué)演化路徑，初步揭示了地?zé)崴h(huán)機(jī)理。

筆者基于前人研究成果，在菏澤斷裂北部鄄城地?zé)崽锊杉说責(zé)崴皻怏w樣品，對水化學(xué)組分、氫氧同位素、碳同位素、氣體組分及同位素等進(jìn)行了實(shí)驗(yàn)測試，結(jié)合以往收集資料和區(qū)域地質(zhì)條件，分析研究了區(qū)內(nèi)地?zé)崴某梢驒C(jī)理。

1 研究區(qū)概況

研究區(qū)為菏澤斷裂北部鄄城地?zé)崽?，區(qū)內(nèi)地?zé)豳Y源豐富，主要熱儲(chǔ)類型為巖溶熱儲(chǔ)，熱儲(chǔ)層呈層狀，分布廣泛，為一套淺海相沉積地層，巖性以灰色、棕灰色灰?guī)r為主，還含有少量的灰色豹皮狀灰?guī)r、泥灰?guī)r、白云質(zhì)灰?guī)r等，巖溶裂隙較發(fā)育，為熱儲(chǔ)層提供了良好的儲(chǔ)水空間，地?zé)崴饕糜诠┡拖丛。?］。研究區(qū)地處華北板塊（Ⅰ）魯西隆起區(qū)（Ⅱ）魯西南潛隆起區(qū)（Ⅲ）菏澤兗州潛斷隆（Ⅳ）菏澤潛凸起（Ⅴ）構(gòu)造單元內(nèi)［6］。受區(qū)域構(gòu)造控制，區(qū)內(nèi)及附近地區(qū)斷裂大致呈網(wǎng)格狀分布，規(guī)模較大，多為凸起和凹陷的界限性斷裂，延伸較遠(yuǎn)，大致可分為近EW向和近SN向兩組。近SN向斷裂主要有聊考斷裂、古城斷裂、侯集斷裂，近EW向斷裂主要有菏澤斷裂、鄆城斷裂。研究區(qū)屬華北平原地層區(qū)魯西地層分區(qū)，地層發(fā)育較為齊全，自老到新依次分布有古生代寒武奧陶系、石炭系、二疊系和新生代新近系、第四系。

研究區(qū)為巖溶裂隙型層狀熱儲(chǔ)，熱儲(chǔ)為巨厚的寒武奧陶紀(jì)碳酸鹽巖，厚度一般大于800m，富水性較好。東部熱儲(chǔ)蓋層主要為第四系和新近紀(jì)松散沉積物，厚度800～1100m，由東向西逐漸增大，地?zé)崴疁?2～50℃，單井涌水量100～150m3/h；西南部熱儲(chǔ)蓋層主要為第四系、新近紀(jì)松散沉積物和石炭二疊紀(jì)砂巖、泥巖，厚度1000～1500m，蓋層總體上東部薄，西部厚，地?zé)崴疁?5～57℃，單井涌水量100～240m3/h。

2 樣品采集與分析

本次共采集水質(zhì)分析樣12組，氣體分析樣2組，氫氧同位素分析樣2組，碳同位素分析樣1組，取樣位置見圖1。其中常規(guī)分析和微量元素分析由山東省魯南地質(zhì)工程勘察院（山東省地質(zhì)礦產(chǎn)勘查開發(fā)局第二地質(zhì)大隊(duì)）實(shí)驗(yàn)測試中心完成，地?zé)崴嘘庪x子的測定采用離子色譜測定，陽離子的測試采用電感耦合等離子體發(fā)射光譜儀測定，微量元素的測試采用電感耦合等離子體光譜/質(zhì)譜測定。

氫氧及碳同位素樣品由美國Beta實(shí)驗(yàn)室測試完成，氫氧同位素采用液體激光同位素分析儀測定，水中14C采用加速質(zhì)譜測定。

氣體樣品由中國科學(xué)院西北生態(tài)環(huán)境資源研究院完成。氣體組分分析采用氣體質(zhì)譜儀測定。稀有氣體同位素3He、4He和20Ne的分析采用稀有氣體質(zhì)譜儀測定，用3He/4He和4He/20Ne來表示。

收集氫氧同位素?cái)?shù)據(jù)9組［24］，氣體組分及同位素?cái)?shù)據(jù)1組［24］，碳同位素?cái)?shù)據(jù)4組［24］。水樣數(shù)據(jù)見表1，氣體樣品數(shù)據(jù)見表2，同位素?cái)?shù)據(jù)見表3。

3 地?zé)崽锏刭|(zhì)特征

3.1 地溫場特征

對比研究區(qū)內(nèi)ZK1地?zé)峋疁y溫?cái)?shù)據(jù)和ZK1西側(cè)的鄄城DR1地?zé)峋疁y溫?cái)?shù)據(jù)［45］（表4，圖2）。ZK1鉆孔全孔地溫梯度2.43℃/100m，蓋層平均地溫梯度2.53℃/100m；鄄城DR1地?zé)峋椎販靥荻?.15℃/100m，蓋層平均地溫梯度3.6℃/100m?？梢钥闯鏊椒较蛏希芯繀^(qū)總體呈現(xiàn)地溫梯度由東向西逐漸增高的趨勢。

垂直方向上，第四紀(jì)地層上部由于靠近地表，使得地層增溫梯度較低，下部由于黏性土層厚度較大，阻擋地下熱能散失的能力較強(qiáng)，成為新近系熱儲(chǔ)良好的蓋層。新近系成巖程度較差，地?zé)嵘⑹С潭容^差，因而地溫梯度較高，說明該層泥巖、砂巖等熱導(dǎo)率低的巖石，是深部地?zé)醿?chǔ)的良好蓋層。石炭二疊系成巖程度較高，熱導(dǎo)率、散熱率都相對較高，該層既能為上部新近系熱儲(chǔ)層傳導(dǎo)熱能，又是深部奧灰熱儲(chǔ)層的保溫蓋層。寒武奧陶系主要巖性為灰?guī)r，巖石的熱導(dǎo)率高，因此地溫梯度低，隨著深度增加，熱儲(chǔ)層溫度增高的程度不大。奧陶紀(jì)灰?guī)r巖溶、裂隙發(fā)育，透水性強(qiáng)，增加了地?zé)岬膫鲗?dǎo)和熱量的散失，地溫梯度較小，是良好的深部熱儲(chǔ)層。

3.2 水化學(xué)特征

3.2.1 水化學(xué)類型

地下水水化學(xué)類型是了解天然水的成因條件和圈劃地下水系統(tǒng)、地表水系統(tǒng)，以及研究兩者間水力聯(lián)系的重要證據(jù)［67］。本文將采集的鄄城地?zé)崽锏責(zé)崴畼影雌渲饕M分（表1）投影到Piper三線圖上顯示，地?zé)崴捌溲a(bǔ)給區(qū)的巖溶地下水可分為Ⅰ、Ⅱ、Ⅲ區(qū)。

第Ⅰ區(qū)為裸露型巖溶補(bǔ)給山區(qū)：東部梁山—嘉祥一帶，水化學(xué)類型以SO4·HCO3·ClCa·Na和HCO3·SO4·ClNa·Mg為主。

Ⅱ區(qū)為覆蓋型巖溶徑流區(qū)：巨野—鄆城一帶，水化學(xué)類型以SO4Ca·Na和SO4Ca型水為主。

Ⅲ區(qū)為深埋藏型巖溶排泄滯留區(qū)：鄄城一帶，水化學(xué)類型主要以SO4Ca·Na型水為主［5］。

鄄城地?zé)崽飵r溶熱儲(chǔ)地?zé)崴畯蘑瘛颉髤^(qū)，在熱儲(chǔ)溫度、壓力和圍巖組分長期作用下水化學(xué)類型由SO4·HCO3·ClCa·Na型轉(zhuǎn)變?yōu)镾O4Ca·Na型，礦化度增大，HCO3濃度減少；SO24、Cl濃度逐漸增大；Ca2+、Mg2+濃度逐漸增大（圖4）。

3.2.2 鈉氯系數(shù)和脫氧系數(shù)

鈉氯系數(shù)（γNa+/γCl）反映了地?zé)崴疂饪s變質(zhì)程度和儲(chǔ)層水文地球化學(xué)環(huán)境［810］，受大氣降水淋濾作用影響的地下水γNa+/γCl一般大于1［12］。脫硫作用都是在缺氧的還原環(huán)境中進(jìn)行，脫硫系數(shù)（100×γSO24/γC1）越小，反映地?zé)崴€原環(huán)境越強(qiáng)，封閉性越好［710］。鄄城地?zé)崽飵r溶熱水含水巖組由海相沉積的石灰?guī)r、白云巖組成，同一區(qū)域地?zé)崴冑|(zhì)系數(shù)、脫硫系數(shù)相差不大，表明地?zé)崴哂邢嗨频乃瘜W(xué)賦存環(huán)境［1314］。通過對變質(zhì)系數(shù)和脫硫系數(shù)特征分析（表5），鄄城地?zé)崽锛案浇鼛r溶熱水的變質(zhì)系數(shù)2.304～2.725，脫硫系數(shù)537.246～646.176，遠(yuǎn)大于1。表明鄄城地?zé)崽锏責(zé)崴艽髿饨邓転V作用影響大，且含水層封閉性較差，大氣降水逐漸淋濾沉積水，形成滲入起源的溶濾水。

3.2.3 同位素化學(xué)特征

（1）氫氧同位素特征。

熱水的δ18O與δD的關(guān)系可以指示熱水的來源［1011］。將地?zé)崴畼与跬凰販y試數(shù)據(jù)（表3）投點(diǎn)至δDδ18O關(guān)系圖上（圖5），可以看出，所有的數(shù)據(jù)點(diǎn)都落在全球大氣降水線及中國大氣降水線附近，說明地?zé)崴鹪从诖髿饨邓?。巖溶冷水?dāng)?shù)據(jù)點(diǎn)落在大氣降水線上，地?zé)峋當(dāng)?shù)據(jù)點(diǎn)均落在降水線下方，氧同位素明顯發(fā)生了向右漂移，這主要是由于水－巖之間同位素的交換作用所致，熱水在地下的滯留時(shí)間長，與圍巖作用的時(shí)間充分，從巖石中獲得的18O越多，“氧漂移”也就越明顯［12］。通過氫氧同位素計(jì)算補(bǔ)給高程，得到各地?zé)峋械責(zé)崴难a(bǔ)給高程為+88.48m～+116.38m（表3）。補(bǔ)給區(qū)高程范圍基本相當(dāng)于梁山、嘉祥一帶基巖山區(qū)高程，研究區(qū)碳酸鹽巖熱儲(chǔ)地?zé)崴赡軄碜杂跂|北部梁山、東部嘉祥基巖山區(qū)的大氣降水入滲，經(jīng)過較遠(yuǎn)距離徑流后補(bǔ)給形成的。

（2）碳同位素。

14C測年的基本原理是應(yīng)用地下水中的溶解無機(jī)碳作為示蹤劑，以14C測定地下水中溶解的無機(jī)碳的年齡［7］。本次研究對鄄城西賈莊地?zé)崴M(jìn)行了14C測定，經(jīng)校正后鄄城西賈莊的14C年齡為31.977ka。本文對實(shí)測、收集的校正后14C年齡綜合氫氧同位素進(jìn)行分析，可知鄄城地?zé)崽锏責(zé)崴饕从诠糯髿饨邓規(guī)r溶水年齡自東向西年齡逐漸變大（表3），說明鄄城地?zé)崽锏責(zé)崴嘌a(bǔ)給區(qū)越遠(yuǎn)，地?zé)崴难h(huán)更新能力越差，可更新能力越差。

3.3 氣體及其同位素特征

地?zé)釟怏w的研究，較之水相組分，地?zé)釟怏w的組分及其同位素具有其特殊優(yōu)勢［1516］。這是因?yàn)闅怏w源于深部，較少受到近地表過程的影響，氣體的化學(xué)特征反映的是系統(tǒng)的流體相與固相間達(dá)到最終平衡時(shí)的狀態(tài)［16］。地?zé)釟怏w的地球化學(xué)和同位素特征可以用來判斷地?zé)崃黧w的來源，識(shí)別殼源組分和幔源組分的比例，以及評(píng)價(jià)殼源熱量和幔源熱量對地表熱量的貢獻(xiàn)比例［16］。

研究區(qū)碳酸鹽巖巖溶熱儲(chǔ)地?zé)崃黧w中氣體的化學(xué)組分如表2所示，主要?dú)怏w組分為N2、O2、CO2。地?zé)嵯到y(tǒng)中N2的可能來源有大氣，以及沉積巖和變質(zhì)巖中的有機(jī)物質(zhì)發(fā)生熱解反應(yīng)產(chǎn)生的N2［16］。本文根據(jù)N2ArHe三種組分之間的相互關(guān)系來對N2的來源進(jìn)行識(shí)別［1619］。如圖6所示，鄄城地?zé)崽锏責(zé)釟怏w中He組分是以地殼的起源為主導(dǎo)?；?guī)r熱儲(chǔ)層氣體樣品落在HeAir趨勢線和HeGroundwater趨勢線之間，且N2/Ar比值介于69～88，接近大氣的N2/Ar比值（84），表明該儲(chǔ)層地?zé)釟怏w中的N2組分既有淺部接受溶解大氣的流體來源，也有深循環(huán)地下水的來源，為兩者的混合來源［16］。因此，初步判定研究區(qū)灰?guī)r熱儲(chǔ)層中的N2組分為大氣起源。

He同位素是識(shí)別地幔來源組分最有效的示蹤劑［16，21］。RA表示大氣的3He/4He比值，值為1.4×106［22］幔源、殼源各自有特征3He/4He值，殼源3He/4He比值為0.1～0.01RA［16，22］；上地幔來源He的3He/4He比值范圍為（8±1）RA，下地幔來源He的3He/4He比值gt;30RA［16，22］。

若流體中He的3He/4He比值gt;0.1RA，意味著含有地幔起源He組分。研究區(qū)地?zé)釟怏w中的He濃度和3He/4He比值如表3所示。通過測定氣體中Ne的濃度和4He/20Ne比值來校正3He/4He比值［1624］，校正后的3He/4He的比值即（R/RA）c為0.11～0.13RA。樣品的3He/4He比值均高于殼源3He/4He，表明樣品中的He含有地幔起源He的加入。可以據(jù)20Ne/4He3He/4He關(guān)系圖來反映地?zé)釟怏w地幔起源He所占的比例（圖7），地?zé)釟怏w幔源He占總He組分的1.13%～1.38%，說明He的來源仍以地殼來源為主。

綜上分析，鄄城地?zé)崽锍尸F(xiàn)出從補(bǔ)給區(qū)的淺循環(huán)—開放式巖溶地下水流系統(tǒng)（Ⅰ），至徑流區(qū)的中循環(huán)—半開放式巖溶地?zé)崴飨到y(tǒng)（Ⅱ），再至排泄滯留區(qū)的深循環(huán)—弱開放式巖溶地?zé)崴飨到y(tǒng)（Ⅲ）的演化路徑。從Ⅰ區(qū)、Ⅱ區(qū)至Ⅲ區(qū)巖溶地下水化學(xué)環(huán)境逐步由開放狀態(tài)演化為較開放至弱開放狀態(tài)。

Ⅰ區(qū)巖溶地下水徑流途徑短、循環(huán)深度淺，地下水化學(xué)環(huán)境處于開放狀態(tài)，能夠接受裸露型巖溶補(bǔ)給及山區(qū)大氣降水入滲補(bǔ)給，變質(zhì)程度低。Ⅱ區(qū)巖溶地下水徑流途徑較遠(yuǎn)、循環(huán)深度較深，巖溶地下水化學(xué)環(huán)境處于半開放狀態(tài)，能夠接受Ⅰ區(qū)巖溶地下水側(cè)向徑流補(bǔ)給，變質(zhì)程度中等。Ⅲ區(qū)巖溶地下水循環(huán)深度增加，距離補(bǔ)給區(qū)更遠(yuǎn)、徑流途徑更長；受西部聊考斷裂等斷裂阻水影響，徑流運(yùn)移速度滯緩，巖溶地下水逐步演化為巖溶地?zé)崴?，地?zé)崴幩瘜W(xué)環(huán)境相對較封閉，變質(zhì)程度高。

4 成因模式

研究區(qū)巖溶熱儲(chǔ)地?zé)豳Y源可以概化為由熱儲(chǔ)層、蓋層、通道、熱源和水源組成的地?zé)崮Ｐ停?］（圖8）。

（1）熱儲(chǔ)層：研究區(qū)熱儲(chǔ)層主要為寒武奧陶紀(jì)灰?guī)r，灰?guī)r層呈層狀，分布廣泛，為一套淺海相沉積地層，灰?guī)r熱儲(chǔ)巖溶裂隙較發(fā)育，為熱儲(chǔ)層提供了良好的儲(chǔ)水空間。

（2）蓋層：鄄城地?zé)崽餆醿?chǔ)蓋層均為密度小、導(dǎo)熱性能差、熱阻大的砂巖、泥巖，為良好的熱儲(chǔ)保溫層。東部蓋層主要為第四系、新近系，厚度800～1100m，由東向西逐漸增大，局部超過1500m；西南部蓋層主要為第四系、新近系、石炭二疊系，厚度1000～1500m，蓋層總體上東部薄，西部厚，局部超過1500m。

（3）通道：在大氣降水補(bǔ)給區(qū)與巖溶熱水區(qū)之間，斷裂構(gòu)造近EW和近SN交錯(cuò)分布，這些斷裂形成于中生代，屬張性斷層，這種性質(zhì)的斷裂，對地下水的富集、運(yùn)移都是有利的，是良好的導(dǎo)水?dāng)嗔褬?gòu)造。且斷裂的破碎溝通了深部熱源的上涌和地下水的深循環(huán)，同時(shí)熱儲(chǔ)層灰?guī)r較發(fā)育的裂隙和溶洞也是較好水源、熱源通道。

（4）熱源：鄄城地?zé)崽餆嵩粗饕獊碜杂诘厍騼?nèi)部熱能傳導(dǎo)，即大地?zé)崃魈籂罴訜?，區(qū)域熱背景較高，大地?zé)崃髦得黠@高于華北盆地的大地?zé)崃髌骄?，表現(xiàn)在地溫場上，整體地溫梯度相對較高，表明其他附加熱源不明顯［2527］。大地?zé)嵩从蓛刹糠謽?gòu)成，即地殼放射性元素產(chǎn)生的熱量和來自上地幔的熱［26］。由于西邊界聊考斷裂為超殼斷裂，溝通了一部分深部熱源，鄄城地?zé)崽飳俚湫偷闹械蜏貍鲗?dǎo)型地?zé)嵯到y(tǒng)。地下水流經(jīng)過長距離的水平、垂向運(yùn)移，充分吸取圍巖的熱量，此時(shí)熱交換以傳導(dǎo)為主；在深入大斷裂附近與深部熱源溝通，此時(shí)熱交換以對流為主［27］。通過地下水的運(yùn)移，對來自深部較均一的熱流在地殼淺部進(jìn)行再分配，將分散在較大范圍內(nèi)地殼中的熱量集中到滲透性好、熱導(dǎo)率高的寒武奧陶紀(jì)灰?guī)r中，從而形成有開采價(jià)值的熱儲(chǔ)。

（5）水源：大氣降水是熱儲(chǔ)層內(nèi)熱水的主要補(bǔ)給源，垂直入滲經(jīng)深部循環(huán)補(bǔ)給地?zé)崴?，地?zé)崴a(bǔ)給主要來自于東北部梁山、東部嘉祥山區(qū)的大氣降水入滲，受西部聊考斷裂等斷裂阻水影響，徑流運(yùn)移速度滯緩，巖溶地下水逐步演化為巖溶地?zé)崴?/p>

5 結(jié)論

（1）鄄城地?zé)崽锏責(zé)崴瘜W(xué)類型主要以SO4·HCO3·ClCa·Na、SO4·HCO3Ca·Na和SO4Ca·Na型為主。

（2）鄄城地?zé)崽锍尸F(xiàn)出從補(bǔ)給區(qū)的淺循環(huán)，開放式巖溶地下水流系統(tǒng)（Ⅰ），至徑流區(qū)的中循環(huán)，半開放式巖溶地?zé)崴飨到y(tǒng)（Ⅱ），再至排泄滯留區(qū)的深循環(huán)，弱開放式巖溶地?zé)崴飨到y(tǒng)（Ⅲ）的演化特點(diǎn)。水化學(xué)類型由SO4·HCO3·ClCa·Na型轉(zhuǎn)變?yōu)镾O4Ca·Na型。

（3）大氣降水是熱儲(chǔ)層內(nèi)熱水的主要補(bǔ)給源，地?zé)崴a(bǔ)給主要來自于東北部梁山、東部嘉祥山區(qū)的大氣降水入滲。

（4）研究區(qū)蓋層為新生代第四系和新近系、古生代二疊系和石炭系；熱儲(chǔ)為寒武奧陶紀(jì)灰?guī)r，正常的大地?zé)崃鳛槠渲饕獰嵩矗瑢俚湫偷闹械蜏貍鲗?dǎo)型地?zé)嵯到y(tǒng)，大斷裂形成了地?zé)崴c深部熱源溝通的通道。

參考文獻(xiàn)：

［1］ 汪集暘，龐忠和，孔彥龍，等.我國地?zé)崆鍧嵢∨a(chǎn)業(yè)現(xiàn)狀與展望［J］.科技促進(jìn)發(fā)展，2020，16（Z1）：294298.

［2］ 馬振民，何江濤，張錫明.菏澤凸起地下熱水的水文地球化學(xué)特征及成因分析［J］.山東地質(zhì)，2000，16（2）：2430.

［3］ 馬哲民，仝路，賈琛，等.山東省菏澤市城區(qū)巖溶熱儲(chǔ)回灌試驗(yàn)［J］.山東國土資源，2018，34（11）：5258.

［4］ 康鳳新，史啟朋，馬哲民，等.盆地潛凸起巖溶熱儲(chǔ)地?zé)崽锍梢驒C(jī)理：以菏澤潛凸起為例［J］.地質(zhì)學(xué)報(bào)，2023，97（1）：221237.

［5］ 史啟朋，宋帥良，孟甲，等.山東省菏澤凸起地?zé)崽飵r溶地?zé)崴瘜W(xué)水平演化特征［J］.中國巖溶，2021，40（2）：310318.

［6］ 馬振民，何江濤，魏加華.鄄城地下熱水的化學(xué)特征及形成機(jī)理［J］.勘察科學(xué)技術(shù)，1999（4）：1821.

［7］ 薄本玉，王韶霞，馬鵬飛，等.基于水化學(xué)的沂源盆地地?zé)崽锍蔁崮Ｊ教接懀跩］.山東國土資源，2022，38（7）：3137.

［8］ 高志友，李大鵬，顏梅，等.基于水文地球化學(xué)特征的黃河流域（濟(jì)陽段）地?zé)崴梢驒C(jī)制研究［J］.山東國土資源，2023，39（3）：6571.

［9］ 孟甲，秦鵬，史啟朋，等.斷陷盆地碳酸鹽巖熱儲(chǔ)勘查及研究：以魚臺(tái)凹陷為例［J］.地質(zhì)科技通報(bào)，2022，41（4）：3845.

［10］ 李賢慶，侯讀杰，唐友軍，等.地層流體化學(xué)成分與天然氣藏的關(guān)系初探：以鄂爾多斯盆地中部大氣田為例［J］.斷塊油氣田，2002（5）：14.

［11］ 肖睿，龐守吉，祝有海，等.新疆甜水海地區(qū)紅山湖泉水化學(xué)特征及其意義［J］.物探與化探，2023，47（1）：3946.

［12］ 蘇艷.關(guān)中盆地中部地壓型地?zé)崃黧w環(huán)境同位素水文地球化學(xué)特征研究［D］.西安：長安大學(xué)，2009：120.

［13］ 隋海波，康鳳新，李常鎖，等.水化學(xué)特征揭示的濟(jì)北地?zé)崴c濟(jì)南泉水關(guān)系［J］.中國巖溶，2017，36（1）：4958.

［14］ 江露露，隋海波，康鳳新，等.魯中隆起北緣地?zé)釁^(qū)巖溶熱儲(chǔ)水化學(xué)特征及形成機(jī)理［J/OL］.中國巖溶：123［20231204］.http：//kns.cnki.net/kcms/detail/45.1157.P.20230807.1037.002.html.

［15］ 崔洋，康鳳新，鐘振楠，等.魯西北平原地?zé)釤嵩礄C(jī)制的氣體同位素約束［J］.地球?qū)W報(bào)，2023，44（1）：93106.

［16］ 龐菊梅.牛駝鎮(zhèn)地?zé)崽飵r溶熱儲(chǔ)對規(guī)?；晒嗟捻憫?yīng)機(jī)理研究［D］.北京：中國科學(xué)院大學(xué)，2018：120.

［17］ 鄭淑蕙，侯發(fā)高，倪葆齡.我國大氣降水的氫氧穩(wěn)定同位素研究［J］.科學(xué)通報(bào)，1983（13）：801806.

［18］ FISCHER T P，GIGGENBACH W F，SANO Y，et al.Fluxes and sources of volatiles discharged from Kudryavy， a subduction zone volcano， Kurile Islands［J］.Earth Planet.sci.lett，1998，160（1/2）：8196.

［19］ HEATON T H E，VOGEL J C.\"Excess air\" in groundwater［J］.Journal of Hydrology， 1981， 50：201208.

［20］ MAMYRIN B A，ANUFRIEV G S，KAMENSKII I L，et al.Determination of the isotopic composition of atmospheric helium［J］.Geochemistry International，1970（7）：498505.

［21］ OZIMA M，PODOSEK F A.Noble gas geochemistry［M］.Cambridge Univ.Pr，1983：367.

［22］ YU J，SAN O，HIRO S，et al.Geographical distribution of 3He/4He ratios in Japan： Implications for arc tectonics and incipient magmatism［J］.Journal of Geophysical Research，1985，90（B10）：87298741.

［23］ CINTI D，TASSI F，PROCESI M，et al.Quattrocchi F（2014） Fluid geochemistry and geothermometry in the unexploited geothermal field of the Vicano–Cimino Volcanic District （Central Italy）［J］.Chem Geol 371：96114.

［24］ 龐菊梅.牛駝鎮(zhèn)地?zé)崽飵r溶熱儲(chǔ)對規(guī)?；晒嗟捻憫?yīng)機(jī)理研究［D］.北京：中國科學(xué)院大學(xué)，2018：120.

［25］ 宋嘉佳，王貴玲，邢林嘯，等.巖石熱導(dǎo)率校正對大地?zé)崃饔?jì)算值的影響：以渤海灣盆地冀中坳陷為例［J］.地質(zhì)論評(píng)，2023，69（4）：13491364.

［26］ 羅文行，孫國強(qiáng)，周洋，等.試論地球深部地?zé)崮軅鬏敊C(jī)理［J］.地學(xué)前緣，2020，27（1）：1016.

［27］ 王澤龍.北京市小湯山地區(qū)地溫場特征及地下熱水成因模式分析［D］.北京：中國地質(zhì)大學(xué)（北京），2007：120.

Study on Geothermal Chemical Characteristics and Genesis of Juancheng Geothermal Fields in Heze Area

MA Long1，2，3， JIA Chen1，2，3， ZHANG Ye1，3， WANG Huafeng1，3， ZHANG Xiaofei1，3

（1.Lunan Geo-engineering Exploration Institute （No.2 Geological Brigade of Shandong Provincial Bureau of Geology and Mineral Resources）， Shandong Ji'ning 272100， China; 2.Shandong Provincial Engineering Research Center of Geothermal Energy Exploration and Development， Shandong Ji'ning 272100，China; 3.Heze Key Laboratory of Geothermal Clean Energy Exploration， Development and Utilization， Shandong Heze 274000， China）

Abstract：Juancheng geothermal field in the northern of Heze fault is rich in geothermal resources. It is of great significance to explore the formation mechanism of regional geothermal resources for the development and utilization of geothermal resources and economic development in the northern part of Heze fault. In this paper， the collected samples have been analyzed by hydrogeochemistry， gas and isotope tests. It is showed that hydrochemical types of geothermal water in Juancheng geothermal field are mainly SO4·HCO3·ClCa·Na， SO4·HCO3Ca·Na and SO4Ca·Na. Juancheng geothermal field shows the evolution path from the shallow circulationopen karst groundwater flow system （Ⅰ） in the recharge area， to the middle circulationsemiopen karst geothermal water flow system （Ⅱ） in the runoff area， and then to the deep circulationweak open karst geothermal water flow system （ Ⅲ） in the discharge retention area. The heat source is mainly shell source， which is a typical mediumlow temperature conduction geothermal system.

Key word：Juancheng geothermal field；water chemistry type；isotope chemistry；heating model；Hezhe