盧兆群 宋會軍 程洪柱

中化地質礦山總局山東地質勘查院，山東 濟南 250013

地下熱水是蘊藏在地層與巖石中的流體熱能資源，合理開發(fā)利用必須以地熱地質條件的深入研究為基礎[1]。在地下熱水研究和開發(fā)利用中，熱儲溫度是劃分地熱系統(tǒng)的成因類型和評價地熱資源潛力所不可缺少的重要參數(shù)，在通常情況下難以直接測量，熱儲溫度的確定對于有效利用地熱資源具有非常重要的意義。地熱溫標方法是確定地下深部熱儲溫度的一種經(jīng)濟有效的手段[2]。濟南市平陰縣大孫莊氡溫泉是 1977年在該地區(qū)進行鐵礦普查時發(fā)現(xiàn)的，其熱儲層為新太古代泰山巖群變質巖系，巖性主要為片麻狀黑云角閃二長花崗巖及斜長角閃巖。大孫莊氡溫泉出水水溫28.3℃，屬低溫地熱資源，水中氡、鍶、氟等元素含量均達到了理療熱礦水國家相關標準中的命名礦水濃度，可以命名為氡鍶氟復合型理療熱礦水，尤其是其氡含量之高，在全國罕見，具有很高的理療價值[3]。本文通過對比各種溫標的原理和應用條件，選取合適的溫標離子對大孫莊氡溫泉的熱儲溫度進行估算，獲得地下熱儲溫度的估算結果和相關認識，有助于研究該地區(qū)地下熱水系統(tǒng)溫度場特征和評價地熱資源潛力。

1 常用地熱溫標方法及應用條件

地熱溫標方法主要有二氧化硅地熱溫標、陽離子地熱溫標、同位素地熱溫標和氣體溫標4大類。目前國內外研究較多的地熱溫標有二氧化硅地熱溫標和陽離子地熱溫標，同位素地熱溫標和氣體溫標國內外研究較少[4]。主要對比二氧化硅地熱溫標和陽離子地熱溫標的原理及應用條件，進而選取合適的溫標離子進行熱儲溫度估算。

1.1 二氧化硅溫標

二氧化硅溫標是應用最早也是最常用的地熱溫標，其理論基礎是地熱流體中二氧化硅的含量取決于不同溫度、壓力下石英在水中的溶解度[5]。二氧化硅溶解度隨溫度升高而增加，天然水中溶解的二氧化硅一般不受其它離子的影響，也不受絡合物的形成和揮發(fā)散失的影響，并且沉淀速率隨溫度降低而減慢，因此在地表水中二氧化硅的濃度能很好地指示地下熱儲的溫度[6]。研究表明，溫度小于 110℃時，通常是玉髓控制著溶液中的二氧化硅含量；溫度大于 180℃時，通常是石英控制著溶液中的二氧化硅含量；在 180～110℃間，石英和玉髓都可以和溶液達到平衡[7]。石英溫標要考慮熱水中蒸汽的分離效應和二氧化硅的聚合或沉淀[6]。

常用的二氧化硅地熱溫標溫度相關性公式如下。石英地熱溫標（無蒸汽損失）計算公式：

玉髓地熱溫標計算公式：

式中，t為熱儲溫度，℃；c1為熱水中溶解的H4SiO4形式的SiO2含量，mg/L。

1.2 陽離子溫標

陽離子地熱溫標是基于熱水與固相物質間的K、Na、Ca、Mg等陽離子的交換與溫度的關系建立起來的。所有陽離子溫標方法都是經(jīng)驗性的近似方法[2]，廣泛用于熱儲溫度的評價。常用的有Na-K溫標、Na-K-Ca溫標和K-Mg溫標等。

（1）Na-K溫標僅應用于150℃以上的熱水，尤其是鉆孔中的熱水[2]。低溫條件下水溶Na+/K+一般不受共生堿性長石之間陽離子交換反應的控制，其優(yōu)點是受稀釋和蒸汽分離的影響很小[6]。

Na-K溫標計算公式：

（2）Na-K-Ca溫標是專門用來處理富鈣熱水的。沸騰會使估算值偏高；在許多富Mg2+的中低溫熱水中，Na-K-Ca溫標估算得到的結果也明顯偏高，因此需要進行Mg2+校正[6]。

Na-K-Ca溫標計算公式：

（3）Giggenbach于1988年建立了K-Mg溫標[8]。該溫標適用于低溫地下熱水，估算溫度一般高于熱水井的出水溫度，被認為是繼續(xù)向深部鉆進有可能達到的溫度[9]。

K-Mg溫標計算公式：

式中，t為熱儲溫度，℃；

當 t＜100℃時，β=4/3；當 t＞100℃時，β=1/3；

c2為水中鈉的濃度，mg/L；

c3為水中鉀的濃度，mg/L；

c4為水中鈣的濃度，mg/L；

c5為水中鎂的濃度，mg/L。

除了以上陽離子地熱溫標，還有Mg-Li溫標、K-Li溫標、Na-Li溫標、Na-Ca溫標和K-Ca溫標等，各有不同的應用條件[6]。

2 地熱溫標的選取和計算

2.1 水化學特征

根據(jù)大孫莊氡溫泉井1992年11月～2016年5月期間的6次水化學組分分析數(shù)據(jù)（表1），大孫莊氡溫泉水pH值7.30～7.86（平均7.52），呈弱堿性；礦化度 6411.00～6828.00mg/L（平均6580.09mg/L），總體表現(xiàn)為3g/L≤礦化度≤10 g/L的高礦化度咸水；鍶11.88～14.45mg/L、氟2.07～3.00mg/L、氡 253.23～644.9Bq/L，；水化學類型均為 Cl·SO4-Na·Ca型（舒卡列夫分類）?？梢钥闯?，大孫莊氡溫泉水中各離子成分及微量元素含量多年來變化不大，水質動態(tài)較穩(wěn)定。

表1 大孫莊氡溫泉水不同時期水樣水化學組分分析結果表Table 1 Analysis results of hydrochemical composition of Dasunzhuang radon geothermal spring water samples in different periods

2.2 礦物-流體平衡判斷

地熱流體中溶解物的濃度是熱儲溫度的函數(shù)，使用地熱溫標方法的基本前提是作為地熱溫標的某種溶質或氣體和熱儲中礦物達到了平衡狀態(tài)，因此，必須研究地熱水和礦物的平衡狀態(tài)以檢驗地熱溫標方法的可靠性[10]。

2.2.1 Na-K-Mg三角圖解法

1988年，Giggenbach首先提出了 Na-K-Mg三角圖的方法來計算地熱溫標，在圖中分為完全平衡水、部分平衡水和未成熟水3個區(qū)域，常被用來評價水-巖平衡狀態(tài)和區(qū)分不同類型的水樣，并利用該圖來判斷熱水是否適合Na-K-Mg溫標。其應用原理是，鈉、鉀的平衡調整較緩慢，但鉀、鎂含量的平衡調整的很快，即使在溫度較低時亦如此，因此對中低溫熱田熱儲溫度的計算較為有利。它取決于以下2個依賴于溫度的反應[10]：

其中，三角圖中的坐標計算：

在Na-K-Mg三角圖中，從每個三角點到其對邊，相應組分含量（mg/L）的百分數(shù)由 100%變化到0%，則按Na、K和Mg計算的百分數(shù)作平行于其對邊的兩直線，交點則為各水樣在此圖中位置[10]。

按照上述方法將大孫莊氡溫泉水不同時期的6次水化學組分分析數(shù)據(jù)的Na、K、Mg含量經(jīng)線性轉換后投至Na-K-Mg平衡圖解上（圖1）。

圖1 大孫莊氡溫泉水Na-K-Mg平衡圖解Fig.1 Na-K-Mg equilibrium diagram of Dasunzhuang radon geothermal spring water

從圖1中可以看出，大孫莊氡溫泉水不同時期的6次水化學組分分析數(shù)據(jù)均屬于“部分平衡水”，說明大孫莊氡溫泉水水-巖反應的平衡溫度偏低，水-巖之間尚未達到離子平衡狀態(tài)，溶解作用仍在進行，或是熱水由深部向地表上升的過程中受到了淺層冷水的混合稀釋作用，從而使熱水中元素的含量變低。

因此，陽離子地熱溫標不會計算出合理的平衡溫度值，大孫莊氡溫泉水不適合用陽離子地熱溫標來估算熱儲溫度。

2.2.2 飽和指數(shù)法

飽和指數(shù)是表示地下水對某特定礦物的飽和程度，常用SI表示，它可以用來判斷地下水與某特定礦物的反應狀態(tài)。SI＞0，表示地下水相對某特定礦物是過飽和的，該礦物將從地下水中沉淀出來。SI=0，表示地下水對某特定礦物是平衡的。SI＜0，表示地下水對指定礦物是欠飽和的，地下水有繼續(xù)溶解該礦物的能力[11]。

利用PHREEQC軟件對大孫莊氡溫泉水6次水化學組分分析數(shù)據(jù)進行了水文地球化學模擬，對地下水中幾種常見礦物組分存在形式及礦物飽和指數(shù)進行了計算（表2）。

表2 各礦物飽和指數(shù)Table 2 Saturation index of minerals

從表2中計算結果看，各水樣的幾種常見礦物組分的飽和指數(shù)均不等于0；飽和指數(shù)大于0，礦物處于過飽和狀態(tài)，該礦物決定的溫標將估算出過高的溫度；飽和指數(shù)小于 0，礦物處于非飽和狀態(tài)，估算結果偏低。方解石、玉髓、白云石、螢石、石英的飽和指數(shù)大于 0，表明這些礦物處于過飽和狀態(tài)，在適當條件下將會產(chǎn)生沉淀；硬石膏、石膏的飽和指數(shù)小于 0，說明這些礦物處于非飽和狀態(tài)，不產(chǎn)生沉淀，還會繼續(xù)溶解。

在這幾種礦物中，玉髓和石英的飽和指數(shù)均大于 0，但石英的飽和指數(shù)明顯大于玉髓的，因此，相比較而言，玉髓更接近平衡狀態(tài)，表明玉髓與水反應在多數(shù)情況下要比石英與水反應更接近平衡。

2.3 地熱溫標的選取

準確地估算熱儲溫度的關鍵在于判斷地熱水中相關化學組分是否與圍巖礦物達到化學平衡，Na-K-Mg三角圖解法是目前國內外用以判斷水巖平衡樣點的有效方法。從圖1中各水樣坐標落點來看，6組不同時期的氡溫泉水均處于部分平衡狀態(tài)，不能用陽離子地熱溫標估算熱儲溫度。

飽和指數(shù)法可以用來判斷地下水與某特定礦物的反應狀態(tài)。從表2中幾種礦物的飽和指數(shù)SI來看，方解石、玉髓、白云石、螢石、石英均處于過飽和狀態(tài)，因此可以使用二氧化硅地熱溫標計算公式，但估算結果可能會較實際溫度偏高。

3 熱儲溫度的確定

根據(jù)大孫莊氡溫泉井6組不同時期水樣的水化學組分分析結果（表1），選用二氧化硅地熱溫標對氡溫泉井熱儲溫度進行估算，估算結果見表3。

表3 熱儲溫度估算結果Table 3 Estimated results of the temperature of geothermal reservoir

從表3中的估算結果看，采用石英溫標估算的溫度明顯高于玉髓溫標估算的溫度，這是因為石英的飽和指數(shù)較玉髓的要大很多。由此可見，石英溫標更適用于高溫地熱田。

通過玉髓溫標估算的溫度為 45.7～52.2℃（平均47.8℃），比較接近井口出水水溫28.3℃，因此，利用該計算方法估算熱儲溫度較為合理，估算大孫莊氡溫泉井的熱儲平均溫度為47.8℃。

氡溫泉井目前深度為660m，若加大鉆孔深度，有望達到玉髓溫標估算的溫度。因此，大孫莊地區(qū)具備獲得更高溫度地熱氡溫泉資源的潛力。

4 結論

（1）地熱溫標方法是確定地下深部熱儲溫度的一種經(jīng)濟有效的手段，但是任何溫標在使用前都要進行水-巖平衡判斷，不能直接應用。

（2）通過Na-K-Mg三角圖解法，可以判斷大孫莊氡溫泉水屬于“部分平衡水”，未達到完全平衡狀態(tài)，說明大孫莊氡溫泉水水-巖反應的平衡溫度偏低，水-巖之間尚未達到離子平衡狀態(tài)，溶解作用仍在進行，或是熱水由深部向地表上升的過程中受到了淺層冷水的混合稀釋作用，從而使熱水中元素的含量變低，因此不適合用陽離子地熱溫標估算熱儲溫度。

（3）利用PHREEQC對大孫莊氡溫泉水不同時期的6次水化學組分分析數(shù)據(jù)進行了水文地球化學模擬，對地下水中幾種常見礦物組分存在形式及礦物飽和指數(shù)進行了計算，玉髓和石英的飽和指數(shù)均大于 0，但相比較而言，玉髓更接近平衡狀態(tài)，表明玉髓與水反應在多數(shù)情況下要比石英與水反應更接近平衡。

（4）通過玉髓溫標估算的溫度為45.7～52.2℃（平均47.8℃），比較接近井口出水水溫28.3℃，因此，利用玉髓溫標估算大孫莊氡溫泉井的熱儲溫度較為合理，估算大孫莊氡溫泉井的熱儲平均溫度為47.8℃。

（5）氡溫泉井目前深度為660m，若加大鉆孔深度，有望達到玉髓溫標估算的溫度，大孫莊地區(qū)具備獲得更高溫度地熱氡溫泉資源的潛力。