周 鑫 ，馬致遠(yuǎn)，席臨平，豆惠萍

(長安大學(xué) 環(huán)境科學(xué)與工程學(xué)院，陜西 西安710054)

咸陽城區(qū)熱儲(chǔ)流體混合作用研究及模擬計(jì)算

周 鑫 ，馬致遠(yuǎn)，席臨平，豆惠萍

(長安大學(xué) 環(huán)境科學(xué)與工程學(xué)院，陜西 西安710054)

應(yīng)用水文地球化學(xué)模擬方法對咸陽城區(qū)熱儲(chǔ)流體混合作用分析和混合比例的計(jì)算.結(jié)合研究區(qū)的地?zé)岬刭|(zhì)和水文地質(zhì)條件，進(jìn)行地下熱水混合作用和混合程度的研究.結(jié)果表明，咸陽城區(qū)地下熱水在不同程度上都存在與較低溫度地下熱水的混合作用。應(yīng)用PHREEQC模擬軟件、穩(wěn)定同位素D和溫度計(jì)算熱水與冷水的混合比例，估算在平衡條件下咸陽城區(qū)地下熱水與冷水的混合比例在0.09～0.63之間。

地?zé)崃黧w;混合作用;模擬計(jì)算;PHREEQC

混合作用是指兩種以上不同化學(xué)成分、不同礦化度的地下水混合后，形成一種與原來兩者都不相同的地下水［1］。就地下水而言，其產(chǎn)生主要有三種條件:(1)當(dāng)存在溝通不同含水層或地表水體與含水層水力聯(lián)系的地質(zhì)構(gòu)造時(shí);(2)當(dāng)發(fā)生區(qū)域性或局部越流時(shí);(3)混合開采或成井工藝不佳，導(dǎo)致來自不同含水層水在井管中混合時(shí)［2］?；旌献饔玫难芯?，對認(rèn)識地?zé)崴男纬蓷l件極其化學(xué)成分的演化規(guī)律、對評價(jià)地?zé)崽锏拈_發(fā)潛力等方面都具有非常重要的意義。長期以來，諸多學(xué)者對混合水的水文地球化學(xué)作用展開了深入地研究，包括海水入侵與淡水發(fā)生混合作用的時(shí)的水化學(xué)特征、水文地球化學(xué)作用［3］;地?zé)崴c地下水混合的水化學(xué)特征［4］;地下水在流動(dòng)過程中發(fā)生的混合作用［5］等。

本文以咸陽地?zé)崽餅檠芯繀^(qū)，以地下熱水為研究對象，以水化學(xué)和PHREEQC軟件為主要研究手段，并結(jié)合混合比例計(jì)算方法［6］，研究咸陽地區(qū)地下熱水混合作用及估算混合比例，為咸陽地區(qū)地下熱水的合理開發(fā)與利用提供科學(xué)依據(jù)。

1 研究區(qū)概況

咸陽地?zé)崽锸顷P(guān)中斷陷盆地地?zé)崽锏囊徊糠?，東起涇河與渭河的交匯處，西至武功與扶風(fēng)交界，南臨渭河，北接寶雞峽引渭工程北干渠。熱田內(nèi)已成功開鑿探采結(jié)合井31眼，井深1 600～4 100 m，井口溫度60℃ ～120℃，井口壓力0.141 ～ 4 Mpa，單井涌水量 55 ～306 m3/h［7］。

咸陽地區(qū)主要構(gòu)造形跡是渭河北岸斷裂。此斷裂西起寶雞，經(jīng)武功、咸陽、渭南向東延入豫西，在省內(nèi)延長320 km，切割深度達(dá)20 km以上為一多期隱伏性深大斷裂，其北為咸禮凸起，南為西安凹陷［8］。咸陽地?zé)崽锬壳耙芽辈旌烷_采的有5個(gè)熱儲(chǔ)層位，即第四系三門組熱儲(chǔ)層，上第三系張家坡組熱儲(chǔ)層，上第三系藍(lán)田—灞河組熱儲(chǔ)層，上第三系高嶺群熱儲(chǔ)層，奧陶系巖溶熱水層。研究區(qū)大部分井為藍(lán)田—灞河組熱儲(chǔ)層熱水，該儲(chǔ)層巖性上部為淺棕色泥巖、砂質(zhì)泥巖夾灰白色細(xì)砂巖;下部為棕褐色泥巖與灰白色、灰黃色細(xì)砂巖、粗砂巖、含礫中砂巖呈不等厚互層［9］。

圖1 咸陽地?zé)崽锝煌ㄎ恢眉熬环植紙D

2 水化學(xué)特征及混合作用分析

2.1 散點(diǎn)圖法

為了揭示咸陽地區(qū)地下水離子的分布特征，現(xiàn)把水樣投在三線圖上進(jìn)行分析。從圖2中可以看出，斷裂北部的熱水全部是Cl－Na型水，而斷裂南部的熱水比較復(fù)雜，主要有HCO3SO4– Na、HCO3－ Na、Cl HCO3SO4－ Na、Cl HCO3－ Na這幾種水的類型，冷水得主要成分是HCO3、SO4、Na離子為主，礦化度都很低。而斷裂南邊的一些點(diǎn)與冷水點(diǎn)非常接近，說明在不同程度上斷裂南邊的水有一定的冷水混入，而斷裂北邊的點(diǎn)在piper圖中無法顯示有冷水的混入，同時(shí)也說明用piper圖來分析混合作用有一定的局限性。

我們從piper圖中可以看出，咸陽地區(qū)的熱水溫度各不相同，同時(shí)Cl離子的含量有明顯的變化，因此用Cl離子分別和其他惰性離子還有TDS進(jìn)行分析。在圖3中可以發(fā)現(xiàn)，冷水點(diǎn)都位于坐標(biāo)原點(diǎn)附近，而高溫?zé)崴诘攸c(diǎn)遠(yuǎn)遠(yuǎn)偏離坐標(biāo)原點(diǎn)，而中低溫的水點(diǎn)位于二者之間，由此可以分析出咸陽地區(qū)的熱水有混合作用的生［10］。

圖2 咸陽城區(qū)piper三線圖

圖3 咸陽城區(qū)水化學(xué)散點(diǎn)圖

2.2 Na－K－Mg三角圖

Na－K－Mg三角圖用以反映水化學(xué)平衡狀況及平衡溫度，水樣點(diǎn)位于完全平衡曲線上則說明地?zé)崴_(dá)到或近于平衡狀態(tài)。除此之外，Na－K－Mg三角圖也可用于識別地下熱水的混合作用［11］。冷水的礦化度相對于熱水的低，且水樣點(diǎn)一般落在靠近Mg離子的部分，而熱水點(diǎn)礦化度高，一般遠(yuǎn)離Mg離子的部分。

圖4可以看出咸陽斷裂南部、斷裂北部冷水點(diǎn)都是靠近了Mg離子，咸陽南部的熱水點(diǎn)中，有個(gè)別的壓在平衡曲線上，大部分點(diǎn)靠近平衡曲線，而咸陽北部的熱水點(diǎn)全部位于半飽和區(qū)域，并且在咸陽斷裂南部熱水點(diǎn)與咸陽斷裂北部冷水點(diǎn)之間，說明北部的熱水有冷水的混入，而南部的熱水有少量的冷水混入，沒有北部的熱水混入量大。

圖4 咸陽城區(qū)Na－K－Mg三角圖

2.3 Na－K－Ca地?zé)釡貥?biāo)法

Na－K－Ca地?zé)釡貥?biāo)對識別地下熱水與冷水的混合作用有良好的效果。Fournier和 Tusedell A.H.在1973年研究結(jié)果表明了，當(dāng)Na－K－Ca地?zé)釡貥?biāo)的計(jì)算溫度較大程度的高于溫泉的出口溫度時(shí)，熱水與冷水就可能發(fā)生了混合作用［12］。現(xiàn)將所選的井利用Na－K－Ca地?zé)釡貥?biāo)進(jìn)行計(jì)算，計(jì)算結(jié)果見下表1。

從Na－K－Ca地?zé)釡貥?biāo)遠(yuǎn)大于井口溫度的程度來看，咸陽地區(qū)無論是斷裂以南還是斷裂以北都有不同程度的混合作用。

2.4 水—巖平衡法

一般情況下，地下熱水長期埋藏在地下，與圍巖中的礦物達(dá)到了平衡的狀態(tài)，當(dāng)有冷水的混入時(shí)，混合水與圍巖礦物沒有達(dá)到新的平衡的情況下，混合作用對水—巖平衡關(guān)系的影響集中地表現(xiàn)為礦物 T—SI關(guān)系曲線的向下平移［12］。因此運(yùn)用PHREEQC軟件對每一口進(jìn)行飽和指數(shù)計(jì)算，并且對每一口井，在不同溫度下不同礦物的飽和指數(shù)進(jìn)行計(jì)算。通過計(jì)算結(jié)果繪制溫度與礦物的折線圖。

圖5 咸陽斷裂北岸SP1號井SI－T對比圖

表1 咸陽城區(qū)地?zé)峋疁囟扔?jì)算結(jié)果表

筆者在斷裂以南和斷裂以北各取一口井進(jìn)行說明，圖5、6 a中是在不考慮有冷水混入的情況下繪制的溫度與各礦物的飽和指數(shù)折線圖，從中可以發(fā)現(xiàn)，各個(gè)礦物曲線交匯的不好，主要是由于水—巖反應(yīng)沒有達(dá)到平衡所造成，而b中的圖是取掉一部分的冷水后，所繪制出的溫度與各礦物的飽和指數(shù)折線圖，當(dāng)去掉一部分的冷水以后，各個(gè)礦物的曲線在溫度大于100℃處交匯，說明在熱儲(chǔ)層，溫度在100℃以上時(shí)達(dá)到水—巖平衡。因此，從圖5、6中可以得出，當(dāng)去掉一部分冷水時(shí)，各礦物達(dá)到平衡狀態(tài)，說明有冷水的混入。

圖6 咸陽斷裂南岸帝都花園井SI－T對比圖

3 混合比例計(jì)算

地下水的混合作用可以用18O或2H關(guān)系來分析和定量表達(dá)［13］。且其結(jié)果較大多數(shù)的化學(xué)溶液具有穩(wěn)定，準(zhǔn)確等特點(diǎn)。由于咸陽地區(qū)18O的漂移作用，研究中用D來計(jì)算混合比例?；旌媳壤挠?jì)算原理如下式所示。

其中δ熱取帝都花園的 δ18D最小值為參考點(diǎn)，δ熱=－73.487，δ冷根據(jù)經(jīng)驗(yàn)計(jì)算，δ= －53.2，計(jì)算結(jié)果如表1。

由于研究區(qū)同位素?cái)?shù)據(jù)不足，因此再選取兩種計(jì)算混合比例的方法:用溫度計(jì)算和PHREEQC軟件模擬。

用溫度來計(jì)算混合比例是另外一種行之有效的計(jì)算方法，也能一定程度反應(yīng)冷熱水的混合情況［4］。計(jì)算原理和同位素方法相同。本研究區(qū)的冷水溫度取平均值，T冷=20℃.T熱的取值，根據(jù)前面章節(jié)對熱儲(chǔ)溫度的估算，我們以水巖反應(yīng)最好的步長集團(tuán)的溫度作為參考點(diǎn)，即T熱=130℃，計(jì)算結(jié)果如表1。

利用PHREEQC軟件在多礦物平衡模擬時(shí)，在有冷水混入的情況下，各種礦物無法達(dá)到平衡，因此認(rèn)為的去除一部分冷水后再進(jìn)行軟件模擬，直至各礦物達(dá)到平衡，記下去除水分的量，結(jié)果見表2。

通過三種計(jì)算混合比例的方法得到，運(yùn)用溫度計(jì)算混合比例方法得到的結(jié)果與利用同位素D和軟件模擬得到的結(jié)果相差很大，而用同位素D計(jì)算混合比例和軟件模擬得到的結(jié)果十分接近。由此可見，利用溫度計(jì)算混合比例得到的結(jié)果有很大的誤差，原因可能是數(shù)據(jù)中的溫度反映的是井口溫度，而不是熱儲(chǔ)層的溫度，在計(jì)算混合比例時(shí)出現(xiàn)了誤差;另一方面，在熱儲(chǔ)層中，可能還有其他的熱水混入。

表2 混合比例計(jì)算結(jié)果對比表

4 結(jié)論

(1)研究區(qū)斷裂以南和斷裂以北封閉條件都很良好，但通過對水化學(xué)組分分析、利用Na－K－Mg三角圖法、Na－K－Ca地?zé)釡貥?biāo)法及水—巖平衡法分析得出，咸陽地?zé)崽餆崴忻黠@的混合作用。造成混合作用的原因一方面可能是在地下熱水開采過程中，止水不足，有冷水的混入;還有另一方面可能由于斷裂的導(dǎo)通，導(dǎo)致其他層的低溫水混入造成混合作用。

(2)利用同位素D值和PHREEQC軟件模擬計(jì)算混合比例可以得到，研究區(qū)混合比例的范圍在0.09～0.63之間。地下熱水中Cl離子含量的變化規(guī)律指示了其對應(yīng)的開采強(qiáng)度下混合比例的變化規(guī)律，為合理開發(fā)利用地下熱水提供了依據(jù)。咸陽斷裂南岸的混合作用較斷裂北岸強(qiáng)烈，主要原因是斷裂北岸的井都打在斷裂帶附近，斷裂北岸的導(dǎo)水性較南岸好。

［1］朱學(xué)愚，錢孝星.地下水水文學(xué)［M］.北京:中國環(huán)境科學(xué)出版社.2005.

［2］張曉倫，梁杏，孫敬.奇村地?zé)崽锼瘜W(xué)特征及其混合作用模擬［J］.水文地質(zhì)工程地質(zhì).2007(6).

［3］鄒勝章，朱遠(yuǎn)峰，陳鴻漢.大連大魏家濱海巖溶區(qū)海水入侵化學(xué)過程［J］.海洋地質(zhì)與第四紀(jì)地質(zhì).2004:2月第24卷(1).

［4］鄭西來，郭建青.地?zé)嵯到y(tǒng)中的混合作用及其應(yīng)用研究［J］.西安地質(zhì)學(xué)院學(xué)報(bào).1996:12月第18卷(4).

［5］馬致遠(yuǎn)，錢會(huì)，黃建勛.關(guān)中盆地南部含水層間相互關(guān)系的環(huán)境同位素水文地球化學(xué)證據(jù)［J］.地球科學(xué)與環(huán)境學(xué)報(bào).2006:第26卷(2).

［6］馬致遠(yuǎn)，錢會(huì).水文地球化學(xué)［M］.北京:地質(zhì)出版社.2005.

［7］陳健.咸陽地?zé)崽锏刭|(zhì)構(gòu)造特征與形成機(jī)制研究.陜西煤炭［J］.2009:(2).

［8］劉文輝，胡雪生，邢憲生.咸陽市區(qū)地下熱水化學(xué)特征［J］.陜西地質(zhì).第17卷(2).

［9］李曉華，南小康，蘇英.渭河北岸斷裂對咸陽地?zé)崽锏挠绊懀跩］.安徽農(nóng)業(yè)科學(xué)第37卷(32).

［10］Arnorsson，S，1985，The use of mixing model and chemical geothermometers for estimating underground temperature in geothermal system［J］.Volc，Geotherm.Res.，23，299－335.

［11］杜少平.利用混合模型評價(jià)西藏羊易地?zé)崽餆醿?chǔ)熱流體溫度［J］.1994 － 2011.

［12］錢會(huì).水溶組分平衡分布計(jì)算及其水文地質(zhì)應(yīng)用［M］.西安地圖出版社.2002.9.

［13］馬致遠(yuǎn)，錢會(huì).環(huán)境同位素地下水文學(xué)［M］.陜西科學(xué)技術(shù)出版社.2004.7.

Studies of Mixing and Model Calculation of the Geothermal Systems in Xianyang

ZHOU Xin，MA Zhi－ yuan，XI Lin － ping，DOU Hui－ ping，HOU Chen，YAN Hua，LITing，HU Wei－ wei
(School of Environmental Sci.and Eng.，Chang an University，Xi’an 710054，Shaanxi)

Application of hydro－geochemical simulation method for geothermal fluid mixing analysis and calculation of the mixture ratio in Xian Yang.Combined with the geothermal geology and hydrogeology conditions，the geothermal water mixing and mixed degree for research.The research shows the geothermal water and the colder water existed different degree the mixing action in Xian Yang city.Application of PHREEQC software simulation，stable isotopes D and temperature calculation with the colder water mixture ratio of geothermal water，under the balance conditions，the geothermal water and the colder water of mixture ratio was between 0.09 and 0.63 proportions in Xian Yang city.

Geothermal water;mixture;simulation and PHREEQC

P314

A

1004－1184(2012)03－0036－04

2012－02－15

周鑫(1986－)，男，陜西咸陽人，在讀碩士研究生，主攻方向:環(huán)境同位素地下水文學(xué)。