邱 輝，朱育坤，李 朋，羅 強(qiáng)，張軒國(guó)

(廣東省地球物理探礦大隊(duì)，廣東廣州 510800)

0 引言

廣東河源地區(qū)地?zé)豳Y源豐富(袁建飛，2013)，以水熱型(溫泉)為主。綠色環(huán)保的溫泉地?zé)豳Y源是國(guó)際能源發(fā)展的新趨勢(shì)，溫泉結(jié)合旅游產(chǎn)業(yè)亦被稱作無(wú)煙的朝陽(yáng)產(chǎn)業(yè)(鄭志忠，2017)。2020年廣東省地質(zhì)局開展了廣東省地?zé)豳Y源勘查工作，將省內(nèi)的地?zé)豳Y源分為隆起山地形和沉積盆地型兩大類(甘浩男等，2020)，但未深入研究其水文地球化學(xué)特征，本文選取黃村地?zé)崽?，開展了系統(tǒng)的地?zé)岬刭|(zhì)研究(譚夢(mèng)如等，2019；王貝貝等，2019)，采用地?zé)岬刭|(zhì)調(diào)查、物探勘探、鉆探、地?zé)崃黧w地球化學(xué)分析等多種勘查研究手段，以查明黃村地?zé)崽锏牡責(zé)豳Y源“家底”和開發(fā)利用現(xiàn)狀。

黃村地?zé)崽餅楹釉磾嗔褞媳姸嗟責(zé)崽镏休^為典型的一處，位于河源斷裂帶中段，具有熱儲(chǔ)溫度高、天然熱泉出露多的特點(diǎn)。黃村地?zé)崽镆酝鶅H開展過地?zé)峥辈楣ぷ鳎催M(jìn)行過系統(tǒng)的研究。本次研究在廣東省地?zé)豳Y源勘查黃村重點(diǎn)區(qū)勘查的基礎(chǔ)上，通過分析總結(jié)黃村地?zé)崽锏牡責(zé)岬刭|(zhì)和地?zé)崃黧w化學(xué)特征，為今后本區(qū)的找熱和科研工作提供技術(shù)支撐。

1 區(qū)域地?zé)岬刭|(zhì)背景

廣東河源區(qū)域上位于華南皺褶系南武夷山隆起帶的南部，在地質(zhì)構(gòu)造史上經(jīng)歷過多次復(fù)雜的構(gòu)造階段，處在北東向河源深斷裂帶、紫金大斷裂的夾持部位(邱輝等，2019)，東西向佛岡-豐良深斷裂帶貫穿全區(qū)。區(qū)域巖漿活動(dòng)強(qiáng)烈，出露大量的燕山中晚期花崗巖，屬于南嶺東段構(gòu)造巖漿巖亞帶的一部分。前人研究表明(趙振華等，1998；邱檢生等，2005)，自三疊紀(jì)開始，區(qū)域內(nèi)的巖石圈處于一個(gè)伸展-減薄的拉張構(gòu)造環(huán)境，以白石岡巖體的形成為代表，有利于地球內(nèi)部能量的釋放。多期次、不同時(shí)代且性質(zhì)各異的構(gòu)造形跡為區(qū)域內(nèi)中低溫構(gòu)造裂隙型地?zé)豳Y源的形成提供了前提(史猛等，2019)。

區(qū)域地?zé)釄?chǎng)研究表明(陳墨香，1991)，廣東大地?zé)崃魈幱?1.5～98.2 mW/m2，平均74.0 mW/m2，高于全國(guó)大陸地區(qū)實(shí)測(cè)熱流值平均值(60.9±15.5 mW/m2)。研究區(qū)所在的河源地區(qū)大地?zé)崃髦翟?0～70 mW/m2，區(qū)域大地?zé)崃髦蹈哂谌珖?guó)均值，但低于廣東省均值。區(qū)域熱儲(chǔ)以花崗巖為主，火山巖、變質(zhì)巖和沉積巖次之，大多缺乏蓋層，屬于開啟式型地?zé)嵯到y(tǒng)(陳墨香和鄧孝，1996)。

2 地?zé)崽锏責(zé)岬刭|(zhì)-地球物理特征

2.1 地質(zhì)特征

2.1.1 蓋層及熱儲(chǔ)層特征

結(jié)合前人研究和本次調(diào)查成果(圖1)，黃村地?zé)崽锏纳w層為第四系大灣鎮(zhèn)組(Qdw)河流相碎屑沉積層，主要分布于黃村河流域的沖積平原和河流階地，巖性自下往上為卵石、中粗砂、細(xì)砂、粘土層，厚度只有1～3 m。黃村地?zé)崽飪?nèi)蓋層極薄，甚至不發(fā)育蓋層，導(dǎo)致局部熱儲(chǔ)層裸露，地?zé)崽飪?nèi)已知的熱泉(井)均分布于第四系內(nèi)。

圖1 黃村地?zé)崽锏責(zé)岬刭|(zhì)圖(大地構(gòu)造位置據(jù)趙正等，2012)

黃村地?zé)崽锏臒醿?chǔ)層主要由寒武紀(jì)片麻花崗巖和三疊紀(jì)二長(zhǎng)巖組成，熱儲(chǔ)呈北西向的帶狀分布，受節(jié)理裂隙帶控制(田春艷，2012)。現(xiàn)有地?zé)徙@孔揭露熱儲(chǔ)層埋深2.7～41.9 m，平均埋深21.2 m，控制熱儲(chǔ)層視厚度0.5～20.64 m，平均視厚度13.9 m，熱儲(chǔ)總體屬于淺層熱儲(chǔ)，且厚度較小。片麻狀花崗巖熱儲(chǔ)內(nèi)的地?zé)峋谒疁?5～91.5℃，平均67.8℃，二長(zhǎng)巖熱儲(chǔ)內(nèi)的地?zé)峋谒疁?0～92℃，平均86.5℃，可見二長(zhǎng)巖熱儲(chǔ)層的隔熱性能比片麻花崗巖更好，二長(zhǎng)巖作為圍巖能更好地將溫度保存于構(gòu)造破碎帶中。這是由這兩種巖石的礦物成分差異決定的，二長(zhǎng)巖的主要成分為長(zhǎng)石、普通輝石、普通角閃石、黑云母等低熱導(dǎo)率礦物，石英等高熱導(dǎo)率礦物極少(含量少于5%)，而片麻狀花崗巖則相反，所以其保溫隔熱性能差于二長(zhǎng)巖。如果二長(zhǎng)巖之下還存在深部熱儲(chǔ)，二長(zhǎng)巖就是一個(gè)穩(wěn)定、熱阻率高，具有較強(qiáng)隔熱性能的良好熱儲(chǔ)蓋層(王少輝等，2020)。

2.1.2 構(gòu)造特征

黃村地?zé)崽飻嗔褬?gòu)造極其發(fā)育，主要為NE向、NW和EW向，形成了“棋盤式”的構(gòu)造骨架。從斷裂切錯(cuò)關(guān)系來(lái)看，NE向斷裂為最早期的構(gòu)造，后被NW向的斷裂錯(cuò)斷，最后EW向的構(gòu)造又切斷了NW向構(gòu)造。

NE向斷裂構(gòu)造表現(xiàn)為規(guī)模大、延伸長(zhǎng)的區(qū)域性壓扭性逆斷層，以F1為代表，具有韌性剪切帶的特征，局部可見明顯的熱水蝕變(王海平等，2002；孔令添等，2019)，發(fā)育皮殼狀、鮞狀蛋白石，是地?zé)崽飪?nèi)主要的控?zé)針?gòu)造。

NW向斷裂構(gòu)造主要是F3～F5，表現(xiàn)為近平行排列的密集節(jié)理裂隙帶。圖2顯示，節(jié)理傾向南西為主，性質(zhì)為張扭性，裂隙內(nèi)充填了石英脈，而后石英脈又發(fā)生熱水蝕變，形成晶洞、晶芽。圖1顯示，現(xiàn)存地?zé)峋黠@沿F3～F4集中分布，在斷裂的傾向上盤均能揭露到熱礦水。NW向斷裂構(gòu)造是地?zé)崽飪?nèi)最主要的導(dǎo)熱導(dǎo)水構(gòu)造。

EW向斷裂構(gòu)造應(yīng)為區(qū)域佛岡-豐良斷裂的一部分，表現(xiàn)為平移走滑斷裂，形成了許多密閉的剪節(jié)理(圖2)，且圍巖發(fā)生了硅化，形成了致密的石英巖，EW向斷裂是一組隔熱構(gòu)造。

圖2 地?zé)崽锕?jié)理傾向玫瑰花圖

NE向的F2、NW向的F3、F5和EW向的F7交織成了一個(gè)梯形的地?zé)崽镞吔?，?jīng)多個(gè)鉆孔驗(yàn)證，邊界外均未揭露到中高溫?zé)岬V水。

2.2 地球物理特征

為了進(jìn)一步了解黃村地?zé)崽锉蔽飨蚩氐V斷裂(F3～F5)在深部的規(guī)模和產(chǎn)狀，在地?zé)崽锖髰徱粠ч_展AMT測(cè)量，共2條北東向的測(cè)線(圖1)，采用Aether大地電磁測(cè)深儀器系統(tǒng)，以四分量(Ex、Ey、Hx、Hy)矢量方式進(jìn)行數(shù)據(jù)采集。

2.2.1 物性特征分析

地?zé)崽飪?nèi)巖性主要第四系大灣鎮(zhèn)組砂層、上侏羅統(tǒng)上龍水組粉砂巖、三疊世二長(zhǎng)巖和早寒武世片麻狀花崗巖。松散巖類的砂層中富含孔隙水，電阻率一般小于100 Ω·m；層狀巖類的粉砂巖中層狀構(gòu)造，發(fā)育密閉不含水的節(jié)理，電阻率一般在1000～2000 Ω·m；塊狀巖類的二長(zhǎng)巖和片麻狀花崗巖淺部強(qiáng)-中風(fēng)化層富含風(fēng)化裂隙水，電阻率一般在300～1000 Ω·m，深部微風(fēng)化巖體為致密塊狀構(gòu)造，電阻率一般3000～10000 Ω·m。構(gòu)造破碎帶中由于富含裂隙水，電阻率一般小于100 Ω·m，具有明顯的低阻特征。

2.2.2 AMT測(cè)量異常特征

圖3顯示，AMT1和AMT2兩條平行測(cè)線的電性結(jié)構(gòu)較為一致，自SW至NE表現(xiàn)為高低阻相間排列，深部高阻、中淺部低阻的異常特征。

圖3 音頻大地電磁法反演電阻率斷面及綜合解釋成果圖

兩條測(cè)線中不連續(xù)分布的橢圓狀、株?duì)罡咦璁惓?ρ>3000 Ω·m)應(yīng)為二長(zhǎng)巖體和片麻狀花崗巖體的反映，與地表出露的巖體位置亦一致。

2.3 地溫場(chǎng)特征

黃村地?zé)崽飪?nèi)現(xiàn)有27口地?zé)峋?，均為自流井，由于水頭較高出現(xiàn)噴涌現(xiàn)象，所有地?zé)峋诰捎昧朔忾]處理，故不能對(duì)各井進(jìn)行詳細(xì)的地溫?cái)?shù)據(jù)采集，僅對(duì)各井井口出水溫度采用溫度計(jì)進(jìn)行測(cè)量。

27個(gè)地?zé)峋谒疁?5～91.5℃，平均69.7℃。圖4顯示，水溫各溫度段均有分布，以中高溫的地?zé)崃黧w為主，即水溫80～90℃段的地?zé)峋疃?，?shù)量達(dá)11個(gè)；水溫大于90°的地?zé)峋?個(gè)，分別分布于地?zé)崽锏谋蔽髋抛雍湍蠔|后崗兩端。

圖4 現(xiàn)存地?zé)峋谒疁胤植贾狈綀D

利用各地?zé)峋谒疁乩L制出等溫線圖(圖5)，由于地?zé)崽餆醿?chǔ)埋深較淺，平均埋深21.2 m，且地表調(diào)查發(fā)現(xiàn)，地?zé)峋罹∮?00 m，地?zé)崃黧w來(lái)源基本相同，故井口水溫等溫線圖能大致反映地?zé)崽锏牡販貓?chǎng)分布特征(張偉等，2020)。

圖5 現(xiàn)存地?zé)峋谒疁氐葴鼐€圖

圖5顯示，等溫線長(zhǎng)軸呈北西-南東走向，軸線與F4斷裂基本重合，整體呈一個(gè)多峰橢圓狀，受斷裂控制明顯(任戰(zhàn)利等，2020)。兩個(gè)峰值分別位于排子和后崗，溫度均為91.5℃，70℃和80℃等溫線沿F4斷裂呈串珠狀排列，地?zé)峋xF4斷裂越遠(yuǎn)其溫度越低。

3 地?zé)崃黧w化學(xué)特征

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

由表1可知，黃村地?zé)崽餆岬V水pH值9.3～10.1，平均9.6，屬于堿性水。TDS值在289.2～348.0 mg/L之間，平均289.3 mg/L，屬于淡水(張人權(quán)等，2011)。SiO2含量108.6～117.1mg/L，換算為偏硅酸133.6～144 mg/L，可作為理療礦泉水(楊曉飛和范二川，2020)。F-含量17.66～19.48mg/L，嚴(yán)重超出標(biāo)準(zhǔn)(1mg/L)。野外調(diào)查發(fā)現(xiàn)黃村地?zé)崽镏苓吘幼〉拇迕穸喟l(fā)地氟病中的斑釉癥(潘柏炳和楊天學(xué)，1992)，對(duì)人體健康不利。

表1 熱礦水樣的化學(xué)特征參數(shù)表

3.2 水化學(xué)類型

從圖6以看出，各樣品的投影點(diǎn)在piper圖中非常集中，反映熱礦水均來(lái)自同一循環(huán)通道，具有同源性(白博文等，2020)。地?zé)崃黧w水化學(xué)類型單一，為HCO3-Na型強(qiáng)循環(huán)水(周海燕等，2007)，此類型的熱流其成因與地?zé)崽飪?nèi)的侵入巖有密切關(guān)系。在溫壓較高的條件下，流體分解溶濾巖體中的鈉長(zhǎng)石、云母(楊曉飛和范二川，2020)，使其釋出Na+離子，使地?zé)崃黧w富含Na+離子，形成HCO3-Na型。這種類型的地?zé)崃黧w性質(zhì)活潑，有利于通過深循環(huán)獲取高溫。

圖6 地?zé)崃黧w水化學(xué)Piper圖

3.3 水-巖平衡狀態(tài)分析

Na-K-Mg三角圖解能直觀地反映地下熱礦水的水-巖平衡狀態(tài)(郭寧等，2020；袁星芳等，2020)。利用Aquachem軟件將表1中的有關(guān)數(shù)據(jù)經(jīng)線型處理(史猛等，2019)后投影至Na-K-Mg三角圖，繪制成圖7。

圖7直觀地顯示了5個(gè)樣品均落在未成熟水區(qū)域中，非常接近部分成熟水，表明了地?zé)崽飪?nèi)的地?zé)崃黧w尚未達(dá)到水-巖平衡狀態(tài)。各樣點(diǎn)均十分靠近鎂離子的極大值，表明地?zé)崃黧w處于水巖作用的初級(jí)階段(袁建飛，2013)，其原因可能是由于地?zé)崽飪?nèi)破碎帶、節(jié)理裂隙發(fā)育，各深度間的連通性好，地?zé)崽餅橐粋€(gè)開放型地?zé)嵯到y(tǒng)(史猛等，2019)，深部地?zé)崃黧w在循環(huán)過程中混入了非平衡的淺部冷水。因此，陽(yáng)離子溫標(biāo)不適合用作估算熱儲(chǔ)溫度(譚夢(mèng)如等，2019)。

圖7 Na-K-Mg圖解(據(jù)Giggenbach，1988)

3.4 熱儲(chǔ)溫度估算

3.4.1 SiO2溫標(biāo)估算熱儲(chǔ)溫度

由于溫度、壓力和TDS幾乎不對(duì)SiO2的溶解度產(chǎn)生影響(周訓(xùn)等，2010)，SiO2溫標(biāo)被廣泛用于估算熱儲(chǔ)溫度。

表2中列舉了石英溫標(biāo)-無(wú)蒸汽損失(0～250℃)、玉髓溫標(biāo)-無(wú)蒸汽損失(0～250℃)、α-方英石溫標(biāo)計(jì)算公式(王貝貝等，2019；李狀等，2020)，以此來(lái)進(jìn)行熱儲(chǔ)溫度估算。

表2 熱儲(chǔ)溫度估算結(jié)果表

從表2可知，估算得到的熱儲(chǔ)溫度：T石英：141.9～146.2℃，T玉髓：217.4～223.5℃，T方晶石：91.3～95.7℃，α-方英石地?zé)釡貥?biāo)估算結(jié)果太過于接近井口水溫，整體明顯偏低，此溫標(biāo)不適用于黃村地?zé)崽铮挥袼铚貥?biāo)估算結(jié)果偏離井口水溫過大，亦不合理；當(dāng)出現(xiàn)上述這種情況的時(shí)候，前人研究(王瑩等，2007；霍冬雪等，2019；郭寧等，2020)認(rèn)為采用石英溫標(biāo)作為估算溫度最為合適。

3.4.2 冷熱水混合作用對(duì)地?zé)崃黧w溫度的影響

圖8 地?zé)崃黧wCl-SO4-HCO3圖解(據(jù)黃珣，2019)

如圖9所示，地?zé)崃黧w經(jīng)過深循環(huán)在地表排泄，這個(gè)過程中，由于地?zé)崽飪?nèi)導(dǎo)水構(gòu)造的聯(lián)通，淺部冷水混入了深部的熱流，這種混入作用可能不止發(fā)生一次，因此我們需要對(duì)冷熱水混合的過程加以甄別。

圖9 冷熱水混合模型(據(jù)宋凱，2011)

3.4.3 硅-焓混合模型估算熱儲(chǔ)溫度及冷熱水比例

由于地?zé)崃黧w尚未達(dá)到水-巖平衡狀(圖7)，且混入了不同程度的冷水(圖8～9)，采用硅-焓混合模型既能定量計(jì)算出冷熱水混合比例，又能獲得地?zé)崃黧w的原始溫度。這種方法可用于與SiO2溫標(biāo)做系統(tǒng)比對(duì)，避免采用地?zé)釡貥?biāo)計(jì)算獲得的熱儲(chǔ)溫度偏低這一問題(黃珣，2019；郭寧等，2020)。

硅-焓混合模型原理是：混合前深部熱流的原始溫度(初焓)和冷熱水的比例為未知數(shù)，通過硅-焓方程(1)、(2)求解未知數(shù)。

X1=(S初焓-S終焓)/(S初焓-S冷水焓)

(1)

X2=(SiO2初-SiO2終)/(SiO2初-SiO2冷)

(2)

式(1)和(2)中，X1為焓變化比率；X2為SiO2變化比率；S初焓為深部熱流初始焓；S終焓為混合后熱流的最終焓；S冷水焓為地表冷水焓值；SiO2初深部熱流SiO2的初始含量；SiO2終為混合后熱流SiO2的最終含量；SiO2冷為地表冷水SiO2含量；焓的單位為×4.1868 J/g，SiO2含量單位為mg/L。

小于100℃的飽和水，其焓值等于水的攝氏溫度數(shù)；大于100℃時(shí)，水的溫度、焓值與SiO2的含量關(guān)系可從表3中可以讀出(宋凱，2011)。黃村地?zé)崽飪?nèi)冷水溫度取調(diào)查數(shù)據(jù)平均值20.5 ℃，即S冷水焓值為20.5；因地表冷水TDS多小于20，故SiO2冷取前人經(jīng)驗(yàn)值5 mg/L。R1～R5井口水溫均大于75℃，設(shè)熱流的初始溫度在100～300 ℃。將表3中各溫度熱流的焓值與SiO2含量依次代入方程式(1)和(2)中，求出X1和X2的系列值，然后把系列值在深部熱水溫度與混入冷水比例圖解中進(jìn)行投影，X1和X2曲線的交點(diǎn)在溫度軸上的投影即為深部地?zé)崃黧w的熱儲(chǔ)溫度，在混合比例軸上的投影即為冷熱水的混合比例。

表3 地?zé)崃黧w溫度、焓和SiO2含量關(guān)系

圖10展示了R1～R5樣點(diǎn)經(jīng)過計(jì)算后得出的值投影到圖解中得出的結(jié)果。圖10顯示，由硅焓混合模型得到的熱儲(chǔ)溫度T硅焓=235～242℃，冷水混入比例為72%～75%。硅-焓混合模型估算的熱儲(chǔ)溫度明顯高于石英溫標(biāo)的估算結(jié)果，接近玉髓溫標(biāo)的估算結(jié)果，這樣的估算結(jié)果可能比實(shí)際值高，究其原因可能是由于硅-焓混合模型是一個(gè)理想狀態(tài)下的模式，僅考慮單因子(如單股冷水混入)，實(shí)際情況可能是多種因子相互影響下的結(jié)果，這與前人同類型的研究成果是一致的(宋凱，2011；閆曉雪等，2019)。

圖10 深部熱流溫度與混入冷水比例圖解

3.5 循環(huán)深度估算

深部地?zé)崃黧w通過深循環(huán)和地?zé)嵩鰷孬@得溫度(周訓(xùn)等，2010)，根據(jù)以下公式可以推算出循環(huán)深度H：

H=(T熱儲(chǔ)-T0)/G+h

(3)

式(3)中，T熱儲(chǔ)選取石英溫標(biāo)估算的結(jié)果最為合適(王貝貝等，2019)：T熱儲(chǔ)=T石英=141.9～146.2 ℃；T0為地?zé)崽锬昃鶜鉁兀?0 ℃，G為地溫梯度。由于黃村地?zé)崽飪?nèi)地?zé)峋鶠槊芊饩?，未能系統(tǒng)測(cè)定地溫梯度，故本研究選取前人研究成果4.45℃/100 m，h為恒溫帶深度：20 m。

通過式(3)估算出黃村地?zé)崽锏牡責(zé)崃黧w的循環(huán)深度2759～2856 m。

4 結(jié)論

(1)黃村地?zé)崽锾幱诤釉瓷顢嗔褞Ш妥辖饠嗔褞У膴A持部位，北東向、北西向和東西向斷裂構(gòu)相互交切，具有明顯的控?zé)針?gòu)造(F1)和導(dǎo)水導(dǎo)熱構(gòu)造(F3～F5)，頻繁的巖漿活動(dòng)提供了充足的熱源。

(2)黃村地?zé)崽锏責(zé)崃黧w為堿性淡水，水化學(xué)類型為HCO3-Na型，具有強(qiáng)循環(huán)特點(diǎn)，其SiO2含量達(dá)到理療礦泉水標(biāo)準(zhǔn)，但F元素含量嚴(yán)重超標(biāo)，導(dǎo)致地?zé)崽锔浇用穸鄶?shù)患有斑釉癥。

(3)通過系統(tǒng)的地?zé)崃黧w化學(xué)分析可知，地?zé)崽锏牡責(zé)崃黧w處于水-巖作用的處級(jí)階段，且混入了大比例冷水，冷水混入比例為72%～75%；熱儲(chǔ)溫度T=141.9～146.2℃，進(jìn)一步計(jì)算出熱流的循環(huán)深度為2759～2856 m。

綜上，黃村地?zé)崽锞哂械锰飒?dú)厚的地?zé)岬刭|(zhì)條件，地?zé)崃黧w溫度高，地?zé)崮荛_發(fā)前景廣闊，宜作為地?zé)崮芫C合開發(fā)利用試點(diǎn)進(jìn)行推廣。

[附中文參考文獻(xiàn)]

白博文,平建華,趙繼昌,楊振威,寧藝武,熊超凡.2020.河南淮陽(yáng)縣地?zé)崃黧w化學(xué)特征及其成因分析[J].中國(guó)地質(zhì),47(2)：406-425.

陳墨香.1991.中國(guó)地?zé)豳Y源的分布及其開發(fā)利用[J].自然資源,7(3)：40-46,58．

陳墨香,鄧孝.1996.中國(guó)地?zé)嵯到y(tǒng)類型圖及其簡(jiǎn)要說(shuō)明[J].地質(zhì)科學(xué),31(2)：114-121.

甘浩男,王貴玲,藺文靜,張薇,秦向輝,張重遠(yuǎn).2020.增強(qiáng)型地?zé)嵯到y(tǒng)環(huán)境地質(zhì)影響現(xiàn)狀研究與對(duì)策建議[J].地質(zhì)力學(xué)學(xué)報(bào),26(2)：211-220.

郭寧,劉昭,男達(dá)瓦,孫會(huì)肖,李皓婷,趙海華.2020.西藏昌都覺擁溫泉水化學(xué)特征及熱儲(chǔ)溫度估算[J].地質(zhì)論評(píng),66(2)：499-509.

霍冬雪,周訓(xùn),劉海生,余鳴瀟,張彧齊.2019.云南祥云縣王家莊堿性溫泉水化學(xué)特征與成因分析[J].現(xiàn)代地質(zhì),33(3)：680-690.

黃珣.2019.鮮水河斷裂帶中段溫泉水化學(xué)特征及溫泉成因研究[D].成都：成都理工大學(xué)：1-64.

李狀,周訓(xùn),方斌,沈曄,徐艷秋,陳柄樺,王蒙蒙,隋麗嬡.2020.安徽大別山區(qū)溫泉的水化學(xué)與同位素特征及成因[J].地質(zhì)通報(bào),38(5)：621-630.

孔令添,黎敦朋,張森琦,羅坤.2019.青海貴德曲乃亥溫泉熱水沉積及水熱蝕變研究[J].青海大學(xué)學(xué)報(bào),37(4)：83-89.

潘柏炳,楊天學(xué).1992.地方性氟中毒斑釉癥和齲病的調(diào)查[J].臨床口腔醫(yī)學(xué)雜志,(2)：118-119.

邱輝,謝浴根,李朋,李劍,張軒國(guó).2019.綜合找礦方法在粵東單竹坑鉛鋅銀礦勘查中的應(yīng)用[J].地質(zhì)與勘探,55(6)：1394-1403.

邱檢生,胡建,王孝磊,蔣少涌,王汝成,徐夕生.2005.廣東河源白石岡巖體：一個(gè)高分異的I型花崗巖[J].地質(zhì)學(xué)報(bào),79(4)：503-514.

任戰(zhàn)利,劉潤(rùn)川,任文波,祁凱,楊桂林,崔軍平,楊鵬,張園園.2020.渭河盆地地溫場(chǎng)分布規(guī)律及其控制因素[J].地質(zhì)學(xué)報(bào),94(7)：1938-1949.

史猛,張杰,殷燾,楊寧,江海洋.2019.膠東半島中低溫對(duì)流型地?zé)豳Y源水化學(xué)特征分析[J].地質(zhì)學(xué)報(bào),93(S1)：138-148.

宋凱.2011.西藏沃卡溫泉形成條件及對(duì)隧道工程影響研究[D].成都：成都理工大學(xué)：1-51.

譚夢(mèng)如,周訓(xùn),張彧齊,劉海生,余鳴瀟,海闊.2019.云南勐海縣勐阿街溫泉水化學(xué)和同位素特征及成因[J].水文地質(zhì)工程地質(zhì),46(3)：70-80.

田春艷.2012.廣東省中高溫地?zé)豳Y源勘查與開發(fā)利用建議[J].地下水,34(4)：61-63．

王貝貝,盧國(guó)平,胡曉農(nóng),歐浩.2019.粵西深大斷裂溫?zé)崛瘜W(xué)分析[J].環(huán)境化學(xué),38(5)：1150-1160．

王海平,張?zhí)鞓?王宗良.2002.騰沖地?zé)嵝篃針?gòu)造系統(tǒng)與水熱蝕變作用[J].礦床地質(zhì),21(S1)：1031-1033.

王少輝,王琦,李哲,閆晉龍,楊家林.2020.河南通許凸起尉氏段地?zé)崽餆醿?chǔ)特征及資源評(píng)價(jià)[J].礦產(chǎn)勘查,11(9)：2022-2027.

王瑩,周訓(xùn),于湲,柳春暉,周海燕.2007.應(yīng)用地?zé)釡貥?biāo)估算地下熱儲(chǔ)溫度[J].現(xiàn)代地質(zhì),21(4)：605-612.

楊曉飛,范二川.2020.貴州苗匡熱礦水地?zé)釁^(qū)水化學(xué)特征及熱儲(chǔ)條件分析[J].地質(zhì)與勘探,56(1)：94-101.

閆曉雪,甘浩男,岳高凡.2019.廣東惠州-從化典型地?zé)崽锼牡厍蚧瘜W(xué)特征及成因分析[J].地質(zhì)論評(píng),65(3)：743-754.

袁建飛.2013.廣東沿海地?zé)嵯到y(tǒng)水文地球化學(xué)研究[D].武漢：中國(guó)地質(zhì)大學(xué)：1-131．

袁星芳,王敬,霍光,宋明忠,王京樓,韓忠,鐘振楠,柳祿湧,隋來(lái)倫.2020.膠東半島洪水嵐湯溫泉水化學(xué)特征與成因分析[J].地質(zhì)與勘探,56(2)：427-437.

張偉,蔡鐵剛,王瑞君.2020.新安煤田中深部地溫場(chǎng)特征及其影響因素[J].地質(zhì)與勘探,56(4)：802-808.

張人權(quán),梁杏,靳孟貴,萬(wàn)力,于青春.2011.水文地質(zhì)學(xué)基礎(chǔ)[M].北京：地質(zhì)出版社：1-100.

趙振華,包志偉,張伯友.1998.湘南中生代玄武巖類地球化學(xué)特征[J].中國(guó)科學(xué)(D輯：地球科學(xué)),28(S2)：7-14.

趙正,陳毓川,陳鄭輝,王登紅,曾載淋,趙斌,張家菁.2012.贛南銀坑礦田高山角花崗閃長(zhǎng)巖SHRIMP U-Pb定年及其與成礦的關(guān)系[J].巖礦測(cè)試,31(3)：536-542.

鄭志忠.2017.水熱型地?zé)豳Y源特征研究[J].能源與環(huán)境,(1)：54-55.

周海燕,周訓(xùn),姚錦梅.2007.廣東從化溫泉的水文地球化學(xué)模擬[J].現(xiàn)代地質(zhì),21(4)：619-623.

周訓(xùn),金曉媚,梁四海,沈曄,張紅梅.2010.地下水科學(xué)專論[M].北京：地質(zhì)出版社：32-74.