高 勝 強(qiáng),劉 凱,孫 文 潔,郭 朝 斌
(1.中國(guó)礦業(yè)大學(xué)(北京) 地球科學(xué)與測(cè)繪工程學(xué)院,北京 100083; 2.中國(guó)地質(zhì)科學(xué)院,北京 100037)
江西省明月山地區(qū)位于華南加里東褶皺帶中段北緣,新元古代華南與揚(yáng)子兩大陸板塊的碰撞縫合帶。該區(qū)內(nèi)地?zé)豳Y源豐富,熱水出露點(diǎn)較多,其中以溫湯鎮(zhèn)溫泉較為聞名,在地?zé)崮芾蒙嫌芯薮鬂摿?,地?zé)崴Y源具有良好的經(jīng)濟(jì)價(jià)值,揭示其水化學(xué)特征有利于地?zé)豳Y源保護(hù)和合理開(kāi)發(fā)利用。地?zé)崴Y源作為一種新型清潔能源,近幾年受到密切關(guān)注,得到了廣泛的應(yīng)用[1-3]。近半個(gè)世紀(jì)來(lái),水文地球化學(xué)研究在地?zé)豳Y源的勘查與開(kāi)發(fā)利用中發(fā)揮了重要作用[4],國(guó)內(nèi)外大量學(xué)者對(duì)不同地區(qū)的地?zé)崴瘜W(xué)特征及成因做了分析[5-8]。Giggenbach在1988~1992年間建立了一系列水化學(xué)組分三角圖用于研究地?zé)崃黧w[9]。Frengstad等研究討論了Norway地區(qū)地下水中各個(gè)元素含量與pH之間的相關(guān)關(guān)系[10]。李學(xué)禮、李鷺等運(yùn)用水文地球化學(xué)方法和同位素、水文地質(zhì)等法,對(duì)江西省部分溫泉的起源、熱水年齡、水化學(xué)特征等方面進(jìn)行了研究[11-14]。本文通過(guò)樣品采集分析,根據(jù)水化學(xué)常量組分、微量組分測(cè)試結(jié)果,分析研究區(qū)不同水體的化學(xué)組分和水化學(xué)指標(biāo),查明當(dāng)?shù)厮牡厍蚧瘜W(xué)特征、水質(zhì)類型,有助于研究當(dāng)?shù)氐責(zé)崴乃瘜W(xué)成因,評(píng)價(jià)地?zé)崴?巖平衡狀態(tài),了解不同水體的混合作用、分布規(guī)律,綜合分析地?zé)崴黧w特征。
研究區(qū)地處湘贛邊界的羅霄山脈北段,萬(wàn)龍山-溫湯北東向斷裂構(gòu)造帶兩側(cè),位于114.13°E~114.43°E、27.58°N~27.77°N之間的狹長(zhǎng)地帶,北東走向,屬中亞熱帶季風(fēng)氣候類型,兼具山地氣候特征,四季分明,年平均氣溫為15 ℃,年降雨量為1 624.3~1 832.6 mm,年平均蒸發(fā)量為1 319 mm,相對(duì)濕度為82%。水資源以大氣降水為主,地表水為沖溝水系,流向以北西、南東為主,溪流多流入袁水,地下水主要源出花崗巖體斷裂帶內(nèi)及基巖裂隙含水層,研究區(qū)內(nèi)溫湯富硒溫泉即來(lái)自于花崗巖斷裂帶。萬(wàn)龍山-溫湯斷裂呈北東(東)向展布,傾向100°~120°,傾角50°~70°,自宜春溫湯向西南延至蘆溪張家坊,延伸長(zhǎng)52 km,斷裂帶寬45~150 m不等??煞譃樗榱褞r帶、硅化巖帶和片理化帶。本斷裂形成時(shí)期為印支期-喜山期,且具有多期次活動(dòng)的特點(diǎn)。印支期為韌脆性的逆斷層,形成片理化帶;喜山期為脆性正斷層,形成硅化巖帶和碎裂巖帶。
本次研究的采樣工作時(shí)間主要集中在2020年4月29日至2020年5月17日,主要從研究區(qū)內(nèi)地?zé)峋捌渲車乇硭等?,盡可能包含較大水系,同時(shí)考慮不同水系之間的距離。共采集地?zé)崴?組,地表水22組,地下冷水2組,泉水10組。取樣位置見(jiàn)圖1。將本次樣品送圭瑞測(cè)試科技(北京)進(jìn)行全分析檢測(cè),送中科院進(jìn)行同位素檢測(cè)。采用Origin 2017分析處理數(shù)據(jù),利用AquaChem 2014和MapGIS 67軟件繪制圖件。
研究區(qū)內(nèi)水體的陽(yáng)離子主要由K+、Na+、Ca2+、Mg2+組成,見(jiàn)表1,其中Na+濃度較高,其變化范圍為0.40~46.80 mg/L,平均值為10.73 mg/L。陰離子主要由HCO3-、Cl-、SO42-、NO3-組成,其中HCO3-整體來(lái)說(shuō)濃度較高,其變化范圍為8.90~124.00 mg/L,平均值為38.86 mg/L。研究區(qū)內(nèi)所有水樣均為淡水,弱酸或弱堿性,TDS值在16~156 mg/L之間,平均值為52.05 mg/L;pH的范圍為6.53~8.76,平均值為7.36。
圖1 研究區(qū)水文地質(zhì)圖Fig.1 The hydrogeological map of the study area
表1 研究區(qū)水化學(xué)組分統(tǒng)計(jì)特征值Tab.1 Statistical characteristic values of hydrochemical components in the study area
圖2 研究區(qū)水樣schoeller圖Fig.2 Schoeller of water samples in the study area
表2 研究區(qū)水化學(xué)組分相關(guān)關(guān)系Tab.2 Correlation of hydrochemical components in the study area
圖3 研究區(qū)水化學(xué)Piper圖Fig.3 Piper diagram of hydrochemistry in the study area
根據(jù)Piper三線圖可以很明顯地看出研究區(qū)不同水體樣品的水化學(xué)類型(見(jiàn)圖3),三端元圖可以更加直觀地看出不同類型水樣主要陽(yáng)離子類型(見(jiàn)圖4~5)。地表水主要為Ca(HCO3)2和NaHCO3型水,其中Ca(HCO3)2型水又占很大比重;地?zé)崴訬aHCO3型水為主;泉水以Ca(HCO3)2為主;地下冷水2個(gè)樣品均為NaHCO3型水??梢钥闯觯翰还苣姆N水體陰離子均以HCO3-為主,地表水陽(yáng)離子以Na+、Ca2+為主;地?zé)崴偷叵吕渌?yáng)離子以Na+為主;泉水的主要陽(yáng)離子為Ca2+。
為進(jìn)一步探討研究區(qū)內(nèi)水化學(xué)形成機(jī)制,將研究區(qū)內(nèi)樣品水化學(xué)數(shù)據(jù)導(dǎo)入Gibbs模型[17]繪制Gibbs圖(見(jiàn)圖6),該圖可以直觀地反映“大氣降水”“巖石風(fēng)化”“蒸發(fā)結(jié)晶”對(duì)水的主要成分的控制作用。從圖6中可以看出,研究區(qū)內(nèi)水樣Na+/(Na++Ca2+)、Cl-/(Cl-+HCO3-)均落在Gibbs圖中部,遠(yuǎn)離大氣降水與蒸發(fā)結(jié)晶控制端。由此可知,研究區(qū)內(nèi)水體水化學(xué)特征主要受到水-巖作用的影響,研究區(qū)內(nèi)水體主要受大氣降水補(bǔ)給,在水-巖交互作用下對(duì)水化學(xué)類型的形成有較大的影響。
圖4 研究區(qū)Na-K-Mg三端元圖Fig.4 The ternary of Na-K-Mg in the study area
圖5 研究區(qū)Na-Ca-Mg三端元圖Fig.5 The ternary of Na-Ca-Mg in the study area
1988年Giggenbach提出Na-K-Mg三角圖[18],用來(lái)評(píng)價(jià)地下熱水的水-巖平衡狀態(tài),將地下熱水分為完全成熟水、部分成熟水和未成熟的水。圖7為本研究區(qū)的Na-K-Mg三角圖。由圖7可知:研究區(qū)水樣均落在未成熟水區(qū)域,說(shuō)明地下水為淺循環(huán)水體,熱水在運(yùn)移過(guò)程中受到不同程度的淺層地下水的混合作用影響,其水巖反應(yīng)程度不高。
圖6 研究區(qū)水化學(xué)Gibbs圖Fig.6 Gibbs diagram of hydrochemistry in the study area
圖7 研究區(qū)Na-K-Mg三角圖Fig.7 The Triangle diagram of Na-K-Mg in the study area
在水化學(xué)成分中,許多化學(xué)組分含量之間存在著某種關(guān)系,依據(jù)水化學(xué)特征系數(shù)可以對(duì)地下水的成因、所處環(huán)境和水-巖作用強(qiáng)度等作出分析[19]。研究區(qū)水化學(xué)特征系數(shù)見(jiàn)表3。
鈉氯系數(shù)(γNa/γCl)反映了地下水中Na+的富集程度,可以用來(lái)分析地下水的變質(zhì)程度以及地層的封閉性能。標(biāo)準(zhǔn)海水中γNa/γCl的比例系數(shù)平均為0.85,巖鹽層溶濾形成的地下水γNa/γCl約為1,受大氣降水淋濾作用影響的地下水γNa/γCl一般大于1,經(jīng)過(guò)陽(yáng)離子交替吸附和強(qiáng)烈的水-巖相互作用的沉積水的γNa/γCl小于0.85[20]。從表3可以看出:研究區(qū)內(nèi)地?zé)崴?、地下冷水和泉水γNa/γCl的值均大于1,說(shuō)明研究區(qū)內(nèi)的地下水不是海相沉積水,受到大氣降水淋濾作用影響。
由于海水中的Mg總比Ca大的多,其γMg/γCa約等于5.5,一般地下水不可能達(dá)到如此高值[21]。研究區(qū)內(nèi)地下水γMg/γCa值遠(yuǎn)小于海水比值,結(jié)合鈉氯系數(shù)說(shuō)明研究區(qū)內(nèi)地下水成因主要受大陸溶濾水演化的影響。
表3 水化學(xué)特征系數(shù)Tab.3 Characteristic coefficients of water chemistry
脫硫系數(shù)(100×γSO4/γCl)是指示地下水環(huán)境氧化還原程度的重要指標(biāo),脫硫系數(shù)越小說(shuō)明還原環(huán)境越強(qiáng),同時(shí)地下水封閉性越好,脫硫系數(shù)小于1的地下水一般認(rèn)為還原徹底,埋藏環(huán)境封閉性良好;反之,還原作用不徹底,可能受淺表層氧化作用的影響[22]。研究區(qū)地下水脫硫系數(shù)遠(yuǎn)大于1,說(shuō)明地下水含水層封閉性較差,可能與溫湯鎮(zhèn)地表回灌井回灌有關(guān)。
大氣降水中Sr2+濃度通常很低,相比于地下水可以忽略,所以地下水中Sr濃度很大程度上取決于地下水與含水介質(zhì)的水-巖交互作用[23]。從圖8中可以看出:地表水、泉水中Ca2+濃度與Sr2+濃度有較好的正相關(guān)關(guān)系,一方面是鍶、鈣在元素周期表處于同一主族,具有相同的性質(zhì);另一方面說(shuō)明兩者經(jīng)歷了相似的水-巖作用而地?zé)崴蠸r2+/Ca2+值略高,偏離擬合曲線,結(jié)合離子之間的相關(guān)性分析,主要原因是溫度升高,水-巖作用增強(qiáng),水中Ca2+與巖石中K+、Na+發(fā)生置換。
目前的研究普遍認(rèn)為全球范圍內(nèi)硅酸巖風(fēng)化來(lái)源Sr的87Sr/86Sr比值應(yīng)為0.716~0.720,而碳酸鹽巖風(fēng)化來(lái)源Sr的87Sr/86Sr比值通常為0.708~0.709[24]。本次研究區(qū)內(nèi)樣品87Sr/86Sr比值平均為0.7247,全部樣品87Sr/86Sr比值均大于0.716,同時(shí)Sr2+濃度小于0.05 mg/L,表明研究區(qū)內(nèi)Sr同位素來(lái)源符合硅酸巖風(fēng)化特征。研究區(qū)內(nèi)水樣按照Sr2+濃度可以分為3組(低濃度0.001~0.018 mg/L,中濃度0.023~0.031 mg/L,高濃度0.039~0.049 mg/L)(見(jiàn)圖9),低濃度組出現(xiàn)兩個(gè)87Sr/86Sr比值異常點(diǎn)。低Sr2+濃度主要以地表水和泉水為主,其循環(huán)時(shí)間較短,水-巖交互作用弱,高87Sr/86Sr比值的異常點(diǎn)為地下冷水,可能由于其埋深較大,循環(huán)時(shí)間久導(dǎo)致更多的87Sr同位素溶入水中。低87Sr/86Sr比值的異常點(diǎn)為泉水,處于87Sr/86Sr比值海水均值線上,造成這一原因更多是由于沒(méi)有經(jīng)過(guò)過(guò)多水巖作用,直接來(lái)自大氣降水。中濃度組和高濃度組多為地?zé)崴?,Sr2+含量高可能與非放射性成因的Sr含量高有關(guān),也可能是與溫度升高提高了Sr2+的溶解度有關(guān),其形成機(jī)理需要進(jìn)一步研究[25]。
圖8 Ca2+和Sr2+濃度關(guān)系Fig.8 Relationship between Ca2+ and Sr2+ concentrations
圖9 87Sr/86Sr比值和Sr2+濃度關(guān)系Fig.9 Relationship between87Sr/86Sr and Sr2+ concentrations
(1) 通過(guò)對(duì)研究區(qū)39個(gè)水樣品水化學(xué)組分測(cè)試數(shù)據(jù)的分析得出:地表水主要為Ca(HCO3)2和NaHCO3型水,其中Ca(HCO3)2型水占比較大,地?zé)崴訬aHCO3型水為主,泉水以Ca(HCO3)2為主,地下冷水2個(gè)樣品均為NaHCO3型水。研究區(qū)水體水化學(xué)類型主要受到水巖作用的影響,其中地?zé)崴腡DS、pH高于其他水體,且隨著溫度的升高,Ca2+、Mg2+對(duì)TDS的影響減弱,K+、Na+的影響增強(qiáng)。地?zé)崴乃疁厥躍O42-的影響最為強(qiáng)烈,可作為熱水的指示性離子。對(duì)于地表水和泉水,Ca2+、Mg2+、HCO3-兩兩之間相關(guān)性較強(qiáng),說(shuō)明物質(zhì)來(lái)源比較相似。
(2) 水化學(xué)特征系數(shù)研究表明:研究區(qū)內(nèi)地?zé)崴⒌叵吕渌腿肗a/γCl的值均大于1,結(jié)合γMg/γCa系數(shù),說(shuō)明研究區(qū)內(nèi)地下水成因上主要受大陸溶濾水演化的影響。由脫硫系數(shù)可以得出研究區(qū)內(nèi)地下水含水層封閉性較差,可能與溫湯鎮(zhèn)地表回灌井回灌有關(guān)。
(3) 通過(guò)對(duì)鍶同位素的分析,發(fā)現(xiàn)所采水樣Sr2+、Ca2+表現(xiàn)出很強(qiáng)的正相關(guān)關(guān)系,地?zé)崴暗叵吕渌霈F(xiàn)向Sr2+偏移,說(shuō)明溫度、循環(huán)時(shí)間對(duì)Sr2+產(chǎn)生較大影響。根據(jù)Sr2+的濃度,所采水樣可明顯地分為三組,在87Sr/86Sr比值上沒(méi)有出現(xiàn)很大的差異,表明Sr2+含量升高可能與非放射性成因的Sr含量高有關(guān),也可能是與溫度升高從而提高Sr2+的溶解度有關(guān)。
(4) 通過(guò)對(duì)研究區(qū)內(nèi)不同類型水樣水化學(xué)組分的研究,分析不同類型水樣之間的特征聯(lián)系規(guī)律,可以為明月山地區(qū)地?zé)崴掷m(xù)開(kāi)發(fā)利用提供基礎(chǔ)資料,為整個(gè)環(huán)武功山地區(qū)地?zé)豳Y源的研究提供樣本資料。