盧 麗，陳余道，代俊鴿，王 喆，鄒勝章，樊連杰，林永生，周長(zhǎng)松

(1.桂林理工大學(xué) 環(huán)境科學(xué)與工程學(xué)院，廣西 桂林 541006；2.中國(guó)地質(zhì)科學(xué)院 巖溶地質(zhì)研究所，廣西 桂林 541004；3.桂林理工大學(xué)廣西隱伏金屬礦產(chǎn)勘查重點(diǎn)實(shí)驗(yàn)室，廣西 桂林 541006)

0 引 言

地?zé)豳Y源是寶貴的可再生資源，具有較大的地?zé)崮荛_(kāi)發(fā)利用潛力[1]。川西地區(qū)作為我國(guó)主要的地?zé)豳Y源富集區(qū)域之一，其地?zé)豳Y源量極為豐富[2-3]，僅天然出露溫泉就高達(dá)196處，地?zé)豳Y源總量約為4.26×1016kJ[4]。在金沙江斷裂、甘孜—理塘斷裂、鮮水河斷裂等三大主控?cái)嗔褞У目刂葡耓5]，區(qū)內(nèi)溫泉多沿?cái)嗔褬?gòu)造帶集中分布，并形成了三大地?zé)釒?，分別為德格—巴塘—鄉(xiāng)城地?zé)釒?、甘孜—新龍—理塘地?zé)釒?、爐霍—道孚—康定地?zé)釒6]。前人已對(duì)川西地區(qū)開(kāi)展了相關(guān)研究，主要集中在地?zé)豳Y源分布特征[7-8]、熱儲(chǔ)層及成藏條件[9-10]、溫泉成因模式[11-12]、地球物理探測(cè)應(yīng)用[13-14]等方面，但多數(shù)研究仍屬于“定性”分析階段，少數(shù)定量化研究集中在巴塘、康定等地?zé)釒15-16]，作為三大地?zé)釒е坏母首巍慢垺硖恋責(zé)釒У难芯枯^少。此外，許多學(xué)者也對(duì)該區(qū)域的熱水水化學(xué)方面開(kāi)展過(guò)研究，主要集中在地?zé)崴乃牡厍蚧瘜W(xué)及同位素特征分析[17-19]、水文地球化學(xué)過(guò)程模擬[20-21]以及基于水化學(xué)分析的實(shí)際應(yīng)用研究[22]等方面。

四川昭覺(jué)縣作為甘孜—新龍—理塘地?zé)釒У难由靺^(qū)域，也擁有豐富的地?zé)豳Y源，但其利用率卻不到8%。作為國(guó)務(wù)院扶貧辦掛牌督戰(zhàn)的全國(guó)52個(gè)貧困縣之一，如何利用好豐富的地?zé)豳Y源，也是當(dāng)?shù)鼐珳?zhǔn)扶貧的主要任務(wù)?；诖耍疚倪x擇四川省昭覺(jué)縣的竹核溫泉為研究對(duì)象，并結(jié)合前人研究理論基礎(chǔ)，對(duì)昭覺(jué)竹核溫泉進(jìn)行全面的水文地球化學(xué)特征分析；利用水化學(xué)數(shù)據(jù)，探討竹核溫泉的水來(lái)源，分析補(bǔ)給溫度及補(bǔ)給高程，計(jì)算熱儲(chǔ)溫度、冷水混入比例及熱循環(huán)深度，總結(jié)溫泉成因模式，為揭示川西甘孜—新龍—理塘地?zé)釒У臒嵫h(huán)機(jī)理提供技術(shù)依據(jù)，同時(shí)對(duì)當(dāng)?shù)氐責(zé)豳Y源的開(kāi)發(fā)利用、助力扶貧攻堅(jiān)提供規(guī)劃參考。

1 研究區(qū)概況

研究區(qū)位于上揚(yáng)子古陸塊的四川昭覺(jué)斷陷盆地內(nèi)，盆地內(nèi)為志留—二疊系濱淺海相碎屑巖、碳酸鹽巖、火山巖巖系和三疊—侏羅系陸相碎屑巖系，巖漿活動(dòng)強(qiáng)烈，主體構(gòu)造線近南北向，少許北東向，壓性斷裂發(fā)育，背斜緊密而向斜開(kāi)闊，為古生代末形成的斷陷盆地。區(qū)內(nèi)出露地層主要為志留系黃葛溪組(S1hg)、嘶風(fēng)崖組(S2sf)、大路寨組(S2d)，二疊系陽(yáng)新組(P2y)、峨眉山玄武巖組(P3em)，侏羅系沙溪廟組(J2s)，第四系更新統(tǒng)(Qp)、全新統(tǒng)(Qh)，巖性主要以泥晶灰?guī)r、塊狀玄武巖及砂巖為主。區(qū)內(nèi)河流較發(fā)育，主要發(fā)育有東西向的牛洛河、哈啦黨拉打河，以及南北向的打拉古洛河，其中牛洛河和打拉古洛河是哈啦黨拉打河的兩個(gè)支流(圖1)。

圖1 西川昭覺(jué)竹核溫泉地質(zhì)簡(jiǎn)圖Fig.1 Simplified geologic map of Zhuhe hot springs in Zhaojue, Sichuan Province1.志留系黃葛溪組；2.志留系嘶風(fēng)崖組；3.志留系大路寨組；4.二疊系陽(yáng)新組；5.二疊系峨眉山玄武巖組二段；6.二疊系峨眉山玄武巖組三段；7.三疊系飛仙關(guān)組；8.三疊系銅街子組；9.三疊系嘉陵江組；10.三疊系雷口坡組；11.三疊系須家河組一段；12.三疊系須家河組二段；13.三疊系須家河組三段；14.侏羅系自流井組；15.中更新統(tǒng)冰水堆積；16.晚更新統(tǒng)冰水堆積；17.坡積；18.沖洪積；19.沖積；20.斷層；21.河流；22.產(chǎn)狀；23.溫泉；24.冷泉

竹核溫泉位于四川省大涼山?jīng)錾揭妥遄灾沃菡延X(jué)縣竹核鄉(xiāng)東側(cè)，竹核溫泉包括大溫泉(QB13)和小溫泉(QB13-1)，其出露高程分別為1 892 m、1 968 m，水溫分別為48 ℃、43 ℃。竹核溫泉發(fā)育在二疊系峨眉山組玄武巖二段內(nèi)(P3em2)，熱儲(chǔ)為玄武巖的構(gòu)造破碎帶或裂隙帶，主要受木佛山斷層(F2)和竹核斷層(F6)控制，木佛山斷層為逆斷層，北段為東傾，南段為西傾，傾角32°～85°，破碎帶寬5～50 m；竹核斷層為逆斷層，傾向北西，傾角70°，具有明顯的平錯(cuò)性質(zhì)，斷距為0.5～1.5 km。

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

研究區(qū)共采集溫泉水樣2個(gè)、地表冷泉水樣1個(gè)，采集的樣品后送至自然資源部巖溶地質(zhì)資源環(huán)境監(jiān)督檢測(cè)中心進(jìn)行水質(zhì)全分析以及氫氧同位素分析，采樣點(diǎn)的水化學(xué)和同位素測(cè)試結(jié)果見(jiàn)表1和表2。水質(zhì)全分析根據(jù)《國(guó)家標(biāo)準(zhǔn)飲用天然礦泉水檢驗(yàn)方法》(GB/T8 538—2008)進(jìn)行檢測(cè)，其中陰離子采用離子色譜分析(Dionex-500)，陽(yáng)離子采用電感耦合等離子體發(fā)射光譜儀(ICP-OES)進(jìn)行分析，全分析的測(cè)試精度控制在3%以內(nèi)，檢測(cè)限為0.02 mg/L。水同位素利用Picarro L1102-I液體同位素激光測(cè)試儀進(jìn)行檢測(cè)，測(cè)試結(jié)果以維也納標(biāo)準(zhǔn)平均海洋水(V-SMOW)作為標(biāo)準(zhǔn)，用δD和δ18O(δ=(Rsample/Rstandard-1)×1000‰)的形式表示，檢測(cè)精度：δD為±0.5‰，δ18O為±0.1‰?，F(xiàn)場(chǎng)利用便攜式儀器測(cè)試采樣點(diǎn)水溫和pH。

表1 研究區(qū)水樣水化學(xué)測(cè)試數(shù)據(jù)

表2 研究區(qū)水樣水溫及同位素測(cè)試數(shù)據(jù)

3 結(jié)果分析與討論

3.1 水化學(xué)特征

竹核溫泉水質(zhì)良好，總體上為Ⅲ類水質(zhì)，水無(wú)色透明無(wú)異味，水溫常年保持在40～50 ℃之間。水化學(xué)類型為HCO3-Na型，pH值范圍為7.91～8.09，屬弱堿性水。溫泉水中陽(yáng)離子以Na+為主，含量40.74～50.62 mg/L，其次為Ca2+、

K+、Mg2+等元素；陰離子以HCO3-為主，含量155.23～166.87 mg/L，其次為SO42-、Cl-等元素。除此外，溫泉水中還富含人體所必需的多種微量元素，如鍶、偏硅酸等。其中偏硅酸(53.85～57.12 mg/L)和鍶(0.25～0.26 mg/L)均達(dá)到我國(guó)《飲用天然礦泉水》(GB8537—2018)的界值標(biāo)準(zhǔn)，可將其命名為偏硅酸、鍶型復(fù)合天然飲用礦泉水。因QB13-1(小溫泉)水中氟化物(1.59 mg/L)超標(biāo)，可采用黏土、明礬吸附法來(lái)減少溫泉水中的氟含量[23]，使溫泉水達(dá)到飲用水Ⅱ類水質(zhì)標(biāo)準(zhǔn)。

3.2 氫氧穩(wěn)定同位素特征分析

3.2.1 溫泉水來(lái)源

地下水氫氧同位素特征分析是判斷地下水來(lái)源的重要方法之一，依托此方法可以進(jìn)一步地確定地下水的補(bǔ)給條件，以及大氣降水與地下水之間的水循環(huán)方式及循環(huán)程度。

在20世紀(jì)80年代初，我國(guó)水文地質(zhì)專家鄭淑彗等[24]修正了我國(guó)的大氣降水線為δD=7.9δ18O+8.2。但在不同緯度地區(qū)的水循環(huán)過(guò)程中，由于同位素成分的熱力分餾作用程度相同，導(dǎo)致各地區(qū)的大氣降水線也存在較大差異[25]。根據(jù)張貴玲等[26]和李娜等[27]的綜合分析，川西地區(qū)的大氣降水線為：

δD=7.96δ18O+9.52

(1)

竹核溫泉水樣的δD-δ18O關(guān)系如圖2所示，從中可以看出，竹核溫泉點(diǎn)基本落在大氣降水線附近，表明研究區(qū)內(nèi)溫泉水的主要來(lái)源為大氣降水；同時(shí)也可以看出溫泉點(diǎn)出現(xiàn)了輕微的“負(fù)氧漂移”現(xiàn)象，反映地下熱水與圍巖發(fā)生了一定的水-巖作用。

圖2 研究區(qū)水樣δD-δ18O關(guān)系圖Fig.2 δD vs. δ18O plot for the Zhuhe water samples

3.2.2 補(bǔ)給溫度及補(bǔ)給高程

一般來(lái)說(shuō)，來(lái)自補(bǔ)給區(qū)的水體溫度可以影響溫泉水中同位素的組成以及水中具體元素含量大小[28]，因此本文采用大氣降水中δ18O的溫度效應(yīng)理論來(lái)估算研究區(qū)溫泉水補(bǔ)給區(qū)的平均溫度。因?yàn)檠芯繀^(qū)平均海拔較高，且靠近青藏高原，屬于半干旱的高原季風(fēng)氣候區(qū)，因此本文選擇同緯度的青藏高原東部地區(qū)納木錯(cuò)流域的降水δ18O值和平均溫度的關(guān)系式[29]作為經(jīng)驗(yàn)公式。根據(jù)陳飛等[30]的研究成果，納木錯(cuò)流域冰川融水的補(bǔ)給比例僅為10.02%，降水仍是該流域的主要補(bǔ)給來(lái)源，這與研究區(qū)相似，所以該經(jīng)驗(yàn)公式在研究區(qū)內(nèi)有代表性，具體公式為：

δ18O=0.59t-16.9

(2)

式中：t為平均溫度，℃。

通過(guò)上式計(jì)算出來(lái)的補(bǔ)給區(qū)溫度在13.31～16.10 ℃之間，平均值為14.71 ℃。

大氣降水的氫氧穩(wěn)定同位素δD、δ18O具有高程效應(yīng)：當(dāng)降水高程逐漸增大時(shí)，穩(wěn)定同位素逐漸虧損，即δD、δ18O減小?？梢岳眠@種高程效應(yīng)來(lái)計(jì)算地下水的補(bǔ)給高程，進(jìn)而判斷地下水的可能補(bǔ)給區(qū)[31-33]。由于在高溫作用下，溫泉水中的δ18O常會(huì)出現(xiàn)氧漂移現(xiàn)象，因此在計(jì)算溫泉水補(bǔ)給高程的過(guò)程中，一般選擇δD值計(jì)算更為準(zhǔn)確，而不采用δ18O值。溫泉水補(bǔ)給高程的計(jì)算公式為：

δD=-0.02hALT-27

(3)

式中：hALT為補(bǔ)給高程，m。

通過(guò)上式可以計(jì)算出竹核溫泉的補(bǔ)給高程為3 345～3 560 m。在溫泉北西側(cè)方向約13 km處的木佛山一帶的平均高程3 400～3 600 m，出露峨眉山玄武巖組地層(P3em)，基本符合竹核溫泉的補(bǔ)給區(qū)條件，推測(cè)其補(bǔ)給區(qū)位于木佛山一帶。

3.3 冷熱水混合比例

隨著對(duì)淺層地下水與地?zé)崴拈_(kāi)發(fā)利用力度不斷增強(qiáng)，地?zé)崴谏仙倪^(guò)程中經(jīng)常出現(xiàn)被淺層地下水(冷水)混合的情況，從而改變地?zé)崴幕瘜W(xué)組分，增加TDS含量；因此，研究混合比例是十分必要的。

3.3.1 混合模型

為了估算出混入的冷水比例，F(xiàn)ournier等[34]在1974年提出了硅-焓混合模型，即當(dāng)?shù)叵聼崴械娜芙鈶B(tài)SiO2處于飽和狀態(tài)時(shí)，溫泉水的溫度和SiO2含量是深部熱水初焓的兩個(gè)不同函數(shù)：

Sc·X+Sh(1-X)=Ss

(4)

SiO2c·X+SiO2h(1-X)=SiO2s

(5)

式中：X為混入的冷水比例, %;Sc為冷泉水的焓，J/g；Sh是深部初始熱水的焓，J/g；Ss是溫泉水的焓，J/g；SiO2c是冷泉水的SiO2含量，mg/L；SiO2h是深部初始熱水的SiO2含量，mg/L；SiO2s是溫泉水的SiO2含量，mg/L。

對(duì)以上兩個(gè)公式進(jìn)行適當(dāng)變化，可推出方程X的表達(dá)式：

X1=(Sh-Ss)/(Sh-Sc)

(6)

X2=(SiO2h- SiO2s)/(SiO2h- SiO2c)

(7)

其中：X1和X2分別是基于泉水的溫度和SiO2含量的初焓函數(shù)計(jì)算出的冷水比例。根據(jù)表1和表2中竹核溫泉水的溫度、SiO2含量以及冷泉水的溫度、SiO2含量，并結(jié)合深部初始熱水中溫度、焓與SiO2含量的關(guān)系(表3)，估算出竹核溫泉的X1和X2數(shù)據(jù)(表4)，分別繪制X1、X2曲線，曲線相交點(diǎn)所對(duì)應(yīng)的比例即為溫泉水中混合的冷水比例，如圖3所示。

表3 竹核熱水溫度、焓及SiO2含量之間關(guān)系

表4 硅-焓混合模型計(jì)算竹核溫泉冷水混入比例的結(jié)果

由圖3得出，竹核溫泉中大溫泉的混入冷水比例約為77%，深部熱水溫度約為171 ℃；小溫泉的混入冷水比例約為81%，深部熱水溫度約為187.0 ℃。

圖3 竹核溫泉中大溫泉 (a)和小溫泉(b)假設(shè)熱水溫度與冷水混入比例的關(guān)系Fig.3 Relations between hypothesized water temperature and cold water incursion ratio for large hot spring (a) and small hot spring (b) at Zhuhe

3.3.2 硅-焓圖解法

將研究區(qū)的冷泉、竹核溫泉中小溫泉和大溫泉的SiO2濃度及焓值投影到硅-焓曲線圖[35-37]上，分別記作點(diǎn)A、B、C，它們幾乎在一條直線上，其延長(zhǎng)線與石英溶解度曲線交于點(diǎn)D(圖4)。經(jīng)計(jì)算，大溫泉和小溫泉中冷水的混入比例分別為75.95% (CD/AD)、78.61% (BD/AD)。

圖4 竹核溫泉中大溫泉、小溫泉的硅-焓模型Fig.4 Silicon enthalpy model of the large and small hot springs at Zhuhe

綜上所述，利用混合模型和硅-焓圖解法兩種方法算出來(lái)的大溫泉和小溫泉中冷水混入比例都偏高，可能是由于二者均位于牛洛河河邊，受地表水入滲的影響較大。

3.4 熱儲(chǔ)熱循環(huán)深度

川西地區(qū)屬于隆起山地型地?zé)豳Y源，昭覺(jué)地?zé)釁^(qū)域主要受構(gòu)造和地貌控制[5]。資料顯示，昭覺(jué)地區(qū)年平均氣溫為10.9 ℃，地溫梯度為4.75 ℃/hm，恒溫帶厚度為20 m[38]，其估算循環(huán)深度的公式見(jiàn)下式：

H=(t-t0)/T+h

(8)

式中：H為循環(huán)深度，m；t為熱儲(chǔ)溫度，℃；t0為恒溫帶溫度, ℃，取當(dāng)?shù)仄骄鶜鉁兀籬為恒溫帶深度，m；T為地溫梯度，℃/hm，計(jì)算結(jié)果見(jiàn)表5。

由于地?zé)崴谏仙^(guò)程中冷水混入導(dǎo)致采取的溫泉水樣處于不平衡或部分平衡狀態(tài)，根據(jù)地?zé)釡貥?biāo)法的熱儲(chǔ)溫度計(jì)算的熱循環(huán)深度是中部熱儲(chǔ)的熱循環(huán)深度，經(jīng)過(guò)硅-焓模型校正后，此時(shí)熱儲(chǔ)溫度計(jì)算得到的熱循環(huán)深度才是深部熱儲(chǔ)的最終熱循環(huán)深度。從表5可見(jiàn)，大溫泉和小溫泉最終熱儲(chǔ)溫度為171.0 ℃和187.0 ℃，最終熱儲(chǔ)循環(huán)深度為3 426.38 m和3 766.81 m。

表5 竹核溫泉熱儲(chǔ)熱循環(huán)深度計(jì)算結(jié)果

3.5 溫泉成因模式

該地區(qū)主要熱儲(chǔ)層為峨眉山玄武巖(P3em2、P3em3)，整個(gè)研究區(qū)內(nèi)構(gòu)造裂隙較發(fā)育。區(qū)內(nèi)兩條主控?cái)嗔?木佛山斷層和竹核斷層)規(guī)模大、切割深，與其他斷裂相接形成大型斷裂帶，是連通地下水和深部熱源的通道和地?zé)崴倪\(yùn)移賦存空間，斷裂帶上部分巖石硅質(zhì)蝕變，表明斷裂帶地下水熱活動(dòng)頻繁。

溫泉接收北西方向木佛山地區(qū)的大氣降水補(bǔ)給，大氣降水順著斷裂及各處裂縫向地下入滲，在兩條主控?cái)嗔鸭胺种嗔阎猩钛h(huán)，大地?zé)崃鲗?duì)補(bǔ)給水進(jìn)行加熱。兩條主控?cái)嗔鸭词菬崮軆?chǔ)存的場(chǎng)所，也是熱能傳輸和運(yùn)移的主要通道。在入滲或上升過(guò)程中與淺層地下冷水發(fā)生混合，使得熱水水質(zhì)發(fā)生改變。同時(shí)在補(bǔ)給水下滲的過(guò)程中，會(huì)與周?chē)膰鷰r及特殊礦物發(fā)生水-巖反應(yīng)，萃取出圍巖或礦物中的特殊元素及組分，直至補(bǔ)給水到達(dá)循環(huán)深度；在溫壓差的影響下，經(jīng)地?zé)崃黧w加熱的溫泉水向上逸散，并且沿?cái)嗔褞Ц浇奂?，最終由于地形切割出露地表，形成以“大溫泉、小溫泉”為中心的中-低溫地?zé)豳Y源(圖5)。

圖5 竹核溫泉成因模式Fig.5 Genetic model of the Zhuhe hot springs1.致密狀玄武巖;2.杏仁狀玄武巖;3.灰?guī)r;4.頁(yè)巖;5.細(xì)砂巖;6.泥巖;7.斷層;8.熱水流向;9.水-巖反應(yīng)介質(zhì)交換;10.淺層地下冷水;11.地下冷水混合;12.溫泉

4 結(jié) 論

(1)竹核溫泉的主要熱儲(chǔ)層為峨眉山玄武巖，水化學(xué)類型為HCO3-Na型，屬偏硅酸、鍶型復(fù)合天然飲用礦泉水。溫泉的氘剩余量較大，表明該區(qū)地下熱水埋深較大，在含水層中的滯留時(shí)間長(zhǎng)，徑流速度較慢，并與圍巖發(fā)生了強(qiáng)烈的水-巖作用。補(bǔ)給區(qū)的溫度約為14.71 ℃，補(bǔ)給高程為3 345～3 560 m，推測(cè)其補(bǔ)給區(qū)位于北西方向的木佛山一帶。

(2)利用混合模型和硅-焓圖解法估算出大溫泉的混入冷水比例分別為77%、75.95%，小溫泉的混入冷水比例分別為81%、78.61%，溫泉的冷水混入比例都偏高，可能是由于二者均位于牛洛河河邊，受地表水入滲的影響較大。研究區(qū)的熱儲(chǔ)熱循環(huán)的深度范圍為3 426.38～3 766.81 m。

(3)竹核溫泉受控于木佛山斷層和竹核斷層等主控?cái)嗔?，接收北西方向木佛山地區(qū)的大氣降水補(bǔ)給，然后在深循環(huán)的過(guò)程中與淺層冷水發(fā)生混合，與圍巖發(fā)生水-巖反應(yīng)，最后在溫壓差的影響下地?zé)崴蛏弦萆?，出露地表形成以“大溫泉、小溫泉”為中心的?低溫地?zé)豳Y源。

致謝：本文在撰寫(xiě)過(guò)程中得到評(píng)審專家的指正和自然資源部南方石山地區(qū)礦山地質(zhì)環(huán)境修復(fù)工程技術(shù)創(chuàng)新中心的資助，在此表示衷心感謝。