郭 帥，王維逸，王 治，馬靜晨，姜 輝，王卓卓（北京市地質(zhì)工程勘察院，北京 100048）

淺析北京平原區(qū)熱儲溫度與斷裂關(guān)系

郭 帥，王維逸，王 治，馬靜晨，姜 輝，王卓卓
（北京市地質(zhì)工程勘察院，北京 100048）

本文搜集了38個水樣點的水化學(xué)數(shù)據(jù)，這些水樣點位于北京平原區(qū)的5個地?zé)崽?。水樣點的取水層位均為薊縣系霧迷山組，水樣點揭露的地?zé)崴鶎儆谕粋€含水系統(tǒng)。根據(jù)水化學(xué)類型分析，該含水系統(tǒng)地下熱水徑流方向由西北、西南部向東部流動。利用Na-K-Mg三角圖和飽和指數(shù)選取合適的地?zé)釡貥?biāo)，估算了各水樣點的熱儲溫度范圍。熱儲溫度的空間分布形態(tài)呈現(xiàn)北東向為長軸方向的橢圓形，與黃莊—高麗營斷裂、良鄉(xiāng)—前門斷裂、八寶山斷裂、順義斷裂、南口—孫河斷裂和小湯山斷裂等控?zé)釘嗔延嘘P(guān)，且北東向斷裂對于溫度的控制作用更大，熱儲溫度在斷裂交匯處出現(xiàn)高值區(qū)域，這些高值區(qū)域正是尋找地?zé)岬挠欣貐^(qū)。

地下熱水；水化學(xué)；熱儲溫度；斷裂；北京

0　引言

北京平原區(qū)地?zé)崽锒嘌貥?gòu)造斷裂來進(jìn)行劃分，據(jù)前人資料，這些斷裂多為導(dǎo)水、導(dǎo)熱斷裂。在地下熱水的研究和開發(fā)利用工程中，熱儲溫度是評價地?zé)豳Y源潛力不可缺少的重要參數(shù)，北京平原區(qū)地?zé)峋姸?，尤其是早期地?zé)峋臒醿囟葦?shù)據(jù)難以搜集。地?zé)釡貥?biāo)方法是利用地下熱水中的某些化學(xué)組分的含量與溫度的關(guān)系，估算礦物與水反應(yīng)的平衡溫度，也就是深部熱儲的溫度（汪集旸等，1993）。本文在水樣點水化學(xué)分析的基礎(chǔ)上，根據(jù)Na-K-Mg三角圖和飽和指數(shù)計算，結(jié)合地?zé)釡貥?biāo)應(yīng)用規(guī)律估算水樣點熱儲溫度。并對熱儲溫度與斷裂的關(guān)系進(jìn)行分析，為尋找地?zé)岬挠欣貐^(qū)提供依據(jù)。

1　區(qū)域地質(zhì)背景

1.1構(gòu)造

北京市平原區(qū)斷裂構(gòu)造比較發(fā)育，燕山運動晚期與以升降為主的喜馬拉雅山運動都形成了較大規(guī)模的斷裂，展布方向主要以北東向和北西向為主（增瑞祥等，2001）。

北東向主要斷裂，從西向東有八寶山斷裂、黃莊—高麗營斷裂、順義斷裂、良鄉(xiāng)—前門斷裂、南苑—通縣斷裂等5條較大的斷裂。八寶山斷裂走向北東，斷面傾向南東，傾角較小20°～30°，從西南向北東傾角漸陡。黃莊—高麗營斷裂走向北東，斷裂面傾向南東，傾角約80°，性質(zhì)為高度角正斷層，最大斷距在1000m以上。順義斷裂走向北東，關(guān)于其傾向和傾角直接資料不多，地球物理資料顯示，該斷裂傾向南東，傾角較大。良鄉(xiāng)—前門斷裂，走向北東，一般傾向北西，其傾角也較大。南苑—通縣斷裂總的走向北東，斷裂面傾向北西，傾角50°左右，北西盤下降，南東盤上升，性質(zhì)為正斷層。

北西向斷裂，從東向西依次有南口—孫河斷裂、太陽宮斷裂、永定河斷裂等3條較大的斷裂。南口—孫河斷裂總體走向北西，斷裂北段傾向南西，南段傾向北東，傾角均較大。太陽宮斷裂走向北西，傾向北東，傾角陡。永定河斷裂呈北西向展布，傾向南西，傾角較大。

次級斷裂走向主要分為北東向、北北東向、北西向、南北向和東西向。與本文關(guān)系密切的小湯山斷裂屬于南北向次級斷裂。

據(jù)前人資料（增瑞祥等，2001），黃莊—高麗營斷裂、順義斷裂、良鄉(xiāng)—前門斷裂、南口—孫河斷裂、小湯山斷裂等均為導(dǎo)水、導(dǎo)熱斷裂。

1.2地層

北京市平原區(qū)代表性地層，按由老到新的順序組合而成，中上元古界到第四系分述如下：長城系（Ch）、薊縣系（Jx）、青白口系 （Qn）、寒武系（∈）、奧陶系（O）、石炭系 （C）、二疊系 （P）、三疊系（T）、侏羅系（J）、白堊系（K）、古近系（E）、新近系 （N）和第四系 （Q）。

1.3熱儲層

薊縣系霧迷山組是北京地區(qū)分布面積最廣的熱儲層，在平原區(qū)約占90%以上，巖性主要為硅質(zhì)白云巖、燧石條帶白云巖、紋層狀泥晶白云巖等，厚度一般大于2000m。由于受多次地質(zhì)構(gòu)造運動的影響，霧迷山組裂隙較為發(fā)育，為巖溶作用創(chuàng)造了條件，具有較好的儲存空間和滲透能力，是北京地區(qū)主要的熱儲層。

2　地?zé)崴瘜W(xué)特征

2.1水化學(xué)類型

為了更好的揭示地?zé)崴\移過程中各礦物質(zhì)組分含量的變化，本文只取薊縣系霧迷山組做為北京平原區(qū)地?zé)崴瘜W(xué)主要研究地層，即所選取水樣點的取水層位均為薊縣系霧迷山組。

搜集到水樣點水化學(xué)數(shù)據(jù)共38份，水樣點均為地?zé)峋?，這些地?zé)峋挥跂|南城區(qū)地?zé)崽?、良鄉(xiāng)地?zé)崽铩⒕┪鞅钡責(zé)崽?、雙橋地?zé)崽锖托降責(zé)崽铮▓D1）。其中東南城區(qū)地?zé)崽锏責(zé)峋?5眼（下文以J1-15表示）；良鄉(xiāng)地?zé)崽锏責(zé)峋?眼（下文以L1-7表示）；京西北地?zé)崽锏責(zé)峋?眼（下文以S1-7表示）；雙橋地?zé)崽锏責(zé)峋?眼（下文以O(shè)1-3表示）；小湯山地?zé)崽锏責(zé)峋?眼（下文以T1-6表示）。

根據(jù)各水樣點的水化學(xué)數(shù)據(jù)繪制水樣點的Piper圖。從圖2可知，東南城區(qū)地?zé)崽锼瘜W(xué)類型比較豐富，包括了地下水徑流過程中上游至下游的多種水質(zhì)類型。良鄉(xiāng)地?zé)崽锼瘜W(xué)類型為HCO3·SO4-Ca·Na·Mg、HCO3·SO4-Na，主要為地下水徑流過程中上游的水質(zhì)類型。京西北地?zé)崽锼瘜W(xué)類型為HCO3·SO4-Na·Ca·Mg、HCO3-Na·Ca，主要為地下水徑流過程中上游—中游的水質(zhì)類型。雙橋地?zé)崽锼瘜W(xué)類型為HCO3-Na 、Cl·HCO3-Na·Ca，主要為地下水徑流過程中中游—下游的水質(zhì)類型。小湯山地?zé)崽锼瘜W(xué)類型為HCO3·SO4-Na·Ca·Mg、HCO3-Na·Ca，主要為地下水徑流過程中上游—中游的水質(zhì)類型。由于搜集地?zé)峋∷畬游痪鶠樗E縣系霧迷山組，且均位于北京平原區(qū)地質(zhì)單元，假設(shè)水樣點揭露的地?zé)崴鶎儆谕粋€含水系統(tǒng)。根據(jù)水化學(xué)類型分析，該含水系統(tǒng)地下熱水徑流方向為由西北、西南部向東部流動。

圖1　水樣點位置分布圖（改自（增瑞祥，2001））Fig.1 Location distribution map of the water samples

圖2　水樣點的Piper圖Fig.2 Piper diagram of the water samples

2.2Na-K-Mg三角圖

通過水樣點的Na-K-Mg圖（圖3）可知，除J4、J8、J9地?zé)崴疄椴糠制胶馑?，其?5眼井的地?zé)崴鶎儆谖闯墒焖?，且均集中在Mg1/2角附近。這些水樣點的Mg2+含量高，水-巖相互作用的平衡溫度不高，地下熱水可能發(fā)生了與冷水的混合作用，原則上利用陽離子溫標(biāo)估算未成熟水的平衡溫度不合理，適合用SiO2溫標(biāo)來估算這些水樣點的熱儲溫度（于湲，2006；劉久榮等，2002）。J4、J8、J9地?zé)崴畬儆诓糠制胶馑部赡艽嬖诨旌献饔?，可用SiO2溫標(biāo)和陽離子溫標(biāo)來估算熱儲溫度。

圖3　水樣點的Na-K-Mg三角圖Fig.3 Diagram of Na-K-Mg of the water samples

2.3飽和指數(shù)

飽和指數(shù)可以定量的計算深部熱儲某礦物與水的反應(yīng)程度。礦物飽和指數(shù)等于0，礦物處于飽和狀態(tài)，該礦物決定的溫標(biāo)結(jié)果準(zhǔn)確；飽和指數(shù)大于0，礦物處于過飽和狀態(tài)，該礦物決定的溫標(biāo)估算結(jié)果偏高；飽和指數(shù)小于0，礦物處于非飽和狀態(tài)，該礦物決定的溫標(biāo)估算結(jié)果偏低。

據(jù)表1，L2、L5、O1、O3這4個水樣的玉髓飽和指數(shù)小于0，石英飽和指數(shù)小于0，選取SiO2溫標(biāo)進(jìn)行熱儲溫度估算均偏低，采用SiO2溫標(biāo)和陽離子溫標(biāo)進(jìn)行估算。J3、J5、S2、O2這4個水樣的玉髓飽和指數(shù)小于0，石英飽和指數(shù)大于0，適合采用石英溫標(biāo)進(jìn)行熱儲溫度估算。其余30個水樣的玉髓飽和指數(shù)、石英飽和指數(shù)均大于0，適合采用SiO2溫標(biāo)進(jìn)行熱儲溫度估算。

3　熱儲溫度估算與分析

3.1熱儲溫度估算

熱儲溫度估算常用的地?zé)釡貥?biāo)包括SiO2溫標(biāo)和陽離子溫標(biāo)。

估算地?zé)崴臒醿囟鹊腟iO2溫標(biāo)主要有：石英溫標(biāo)1-無蒸汽分離或混合作用、石英溫標(biāo)2-無蒸汽損失、石英溫標(biāo)3-最大蒸汽損失、玉髓溫標(biāo)1-無蒸汽損失和玉髓溫標(biāo)2-最大蒸汽損失（Fournier et al.，1977）。陽離子溫標(biāo)有Na-K溫標(biāo)（Arnorsson，1983）、Na-K-Ca溫標(biāo)（Fournier R O，1973）和K-Mg溫標(biāo)（Giggenbach W F，1988）。

據(jù)Na-K-Mg三角圖和SiO2礦物飽和指數(shù)分析，對各水樣點的熱儲溫度進(jìn)行估算，分為以下4種：①第1種水樣點：SiO2溫標(biāo)和陽離子溫標(biāo)估算都適用，見表2。②第2種水樣點：SiO2溫標(biāo)估算適用，共31個水樣點，見表3。③第3種水樣點：石英溫標(biāo)估算適用，見表4。

（以下陽離子溫標(biāo)，濃度單位為mg/L）根據(jù)表2、表3和表4，各水樣點直接采用石英溫標(biāo)估算，結(jié)果均高于實測溫度。直接采用玉髓溫標(biāo)估算，結(jié)果偏低，近半水樣點的估算結(jié)果小于實測溫度，與實際不符。K-Mg溫標(biāo)估算結(jié)果均小于實測溫度，與實際不符。在地?zé)釡貥?biāo)選取過程中還應(yīng)考慮其他因素，例如礦物飽和指數(shù)、冷水的混合等（張衛(wèi)民，2001；王瑩等，2007），故本次地?zé)釡貥?biāo)估算給出的是熱儲溫度的可能范圍。

3.2熱儲溫度分布與斷裂的關(guān)系

利用各水樣點的熱儲溫度估算結(jié)果，選取熱儲溫度范圍的平均值進(jìn)行分析，繪制圖4。本文選取水樣點均位于北京市平原區(qū)，各水樣點的地下熱源深度和溫度相差不大，熱儲層及下覆地層相近，可對熱儲溫度與斷裂關(guān)系進(jìn)行分析。

表1　水樣點石英和玉髓礦物的飽和指數(shù)表Tab.1 Saturation indices of quartz and chalcedony of the water samples

水樣點熱儲溫度在東南城區(qū)地?zé)崽锍霈F(xiàn)最大值，熱儲溫度約為140℃。熱儲溫度高于100℃等值線區(qū)域呈現(xiàn)北西向為長軸方向的橢圓形，熱儲溫度高于80℃等值線區(qū)域范圍較廣，呈現(xiàn)北東向為長軸方向的橢圓形。

表2　第1種水樣點熱儲溫度估算結(jié)果（℃）Tab.2 Temperature range of geothermal reservoirs of the first type of water samples （℃）

表3　第2種水樣點熱儲溫度估算結(jié)果（℃）Tab.3 Temperature range of geothermal reservoirs of the second type of water samples （℃）

表4　第3種水樣點熱儲溫度估算結(jié)果（℃）Tab.4 Temperature range of geothermal reservoirs of the third type of water samples （℃）

圖4　熱儲溫度分布圖Fig.4　Distribution map of the temperature of geothermal reservoirs

熱儲溫度與斷裂構(gòu)造的關(guān)系密切，導(dǎo)熱斷裂附近熱儲溫度高，其他斷裂對熱儲溫度影響較小。高于80℃等值線區(qū)域熱儲溫度等值線長軸方向，即北東方向與黃莊—高麗營斷裂、良鄉(xiāng)—前門斷裂、八寶山斷裂、順義斷裂的走向一致，這4條北東向的大斷裂為控?zé)釘嗔?，其斷裂下部與深部熱源接觸，為地下水提供熱源。高于100℃的熱儲溫度等值線區(qū)域長軸方向為北西向，與南口—孫河斷裂和小湯山斷裂的走向基本一致，這2條斷裂也屬于控?zé)釘嗔选?/p>

根據(jù)熱儲溫度分布形態(tài)可確定北東向斷裂對于溫度的控制作用較強，北西向斷裂的控制作用相對較弱。且熱儲溫度高值區(qū)通常位于斷裂交匯處，斷裂交匯處導(dǎo)熱性良好。

南口—孫河斷裂北西側(cè)部分穿過的熱儲溫度等值線上呈現(xiàn)凸出形態(tài)，熱儲溫度較其他位置有明顯升高，與斷裂導(dǎo)熱性一致；小湯山斷裂穿過的熱儲溫度等值線上呈現(xiàn)凸出形態(tài)，熱儲溫度較其他位置有明顯升高，與該斷裂導(dǎo)熱性一致。

4　結(jié)論

水樣點取自位于東南城區(qū)地?zé)崽?、良鄉(xiāng)地?zé)崽?、京西北地?zé)崽铩㈦p橋地?zé)崽锖托降責(zé)崽锏牡責(zé)峋?，取水層位均為薊縣系霧迷山組，且均位于北京平原區(qū)地質(zhì)單元，地?zé)峋衣兜牡責(zé)崴鶎儆谕粋€含水系統(tǒng)。根據(jù)水化學(xué)類型分析，該含水系統(tǒng)地下熱水徑流方向為由西北、西南部向東部流動。

熱儲溫度與斷裂構(gòu)造的關(guān)系密切，導(dǎo)熱斷裂附近熱儲溫度高，其他斷裂對熱儲溫度影響較小。熱儲溫度的空間分布形態(tài)呈現(xiàn)北東向為長軸方向的橢圓形，與黃莊—高麗營斷裂、良鄉(xiāng)—前門斷裂、八寶山斷裂、順義斷裂、南口—孫河斷裂和小湯山斷裂均屬于控?zé)釘嗔延嘘P(guān)，且北東向斷裂對于溫度的控制作用更大，熱儲溫度在斷裂交匯處出現(xiàn)高值區(qū)域，這些高值區(qū)域正是尋找地?zé)岬挠欣貐^(qū)。

利用搜集的水化學(xué)數(shù)據(jù)，本文僅對北京平原區(qū)大部分區(qū)域的徑流場和熱儲溫度進(jìn)行了空間分布的定性分析。下一步應(yīng)搜集更多地?zé)崽锼瘜W(xué)資料，進(jìn)行更深入的分析研究。

［1］ 汪集旸，熊亮萍，龐忠和. 中低溫對流型地?zé)嵯到y(tǒng)［M ］. 北京：科學(xué)出版社，1993∶ 67～82.

［2］增瑞祥，張進(jìn)平，王治，等. 北京市多參數(shù)立體地質(zhì)調(diào)查北京市地?zé)豳Y源潛力調(diào)查評價報告［R］北京：北京市地質(zhì)礦產(chǎn)勘查開發(fā)局，北京市水文地質(zhì)工程地質(zhì)大隊. 2001∶ 12～23，60～73.

［3］于湲. 北京城區(qū)地?zé)崽锏叵聼崴乃瘜W(xué)及同位素研究：［碩士學(xué)位論文］.北京：中國地質(zhì)大學(xué)（北京），2006.

［4］劉久榮，潘小平，楊亞軍，等. 北京城區(qū)地?zé)崽锬车責(zé)峋疅崴厍蚧瘜W(xué)研究［J］. 現(xiàn)代地質(zhì)，2002，16（3） ∶ 318～321.

［5］Fournier.R.O. Chemical geothermometers and mixing models from geothermal systems. Geothermics. 1977， 5∶41～50.

［6］Arnorsson.S. Chemical equilibrium in Icelandic geothermal system∶ implications for chemical geothermometry investigation. Geothermics. 1983，12∶119～128.

［7］Fournier R O， Truesdell A H. An empirical Na-KCa geothermometer for natural water. Geochimica et Cosmochimica Acta. 1973， 37∶ 1255～1257.

［8］Giggenbach W F. Geothermal solute equilibria∶Derivation of Na-K-Mg-Ca geoindicators. Geochimica et Cosmochimica.Acta. 1988， 52∶ 2749～2765.

［9］張衛(wèi)民. 應(yīng)用SiO2地?zé)釡囟扔嫻浪愕責(zé)醿囟取在M南橫涇地區(qū)若干溫泉為例［J］. 地球?qū)W報，2001，22（2） ∶ 185～188.

［10］王瑩，周訓(xùn)，于湲，等. 應(yīng)用地?zé)釡貥?biāo)估算地下熱儲溫度［J］. 現(xiàn)代地質(zhì)，2007，21（4） ∶ 605～612.

The Analysis between Temperature of Geothermal Reservoirs and Fractures in Beijing Plain

GUO Shuai， WANG Weiyi， WANG Zhi， MA Jingchen， JIANG Hui， WANG Zhuozhuo
（Beijing Institute of Geological and Prospecting Engineering， Beijing 100048， China）

This paper collects hydrochemical data from 38 geothermal water samples， and these water samples are located in five geothermal fields in Beijing plain. The water samples are in Wumishan Group of the Jixian System，and the geothermal water all belong to the same aquifer system. Depending on the analysis of the water type， the runoff direction of the geothermal water is from northwest and southwest to the east. By diagram of Na-K-Mg and saturation indices， the suitable geothermometers are selected to calculate the temperature range of geothermal reservoirs of each water samples. The spatial distribution form of the temperatures of geothermal reservoirs is an oval with northeast long axis. As heat controlling faults， Huangzhuang-Gaoliying fault， Liangxiang-Qianmen fault，Babaoshan fault， Shunyi fault， Nankou-Sunhe fault and Xiaotangshan fault control the temperatures of geothermal reservoirs， and the northeast trending faults have mainly control in temperature. The temperature of geothermal reservoirs is higher in fracture intersection zones， and these high value areas are favorable areas for geothermal.

Geothermal water； Hydrochemistry； Temperature of geothermal reservoirs； Fault； Beijing

P314

A

1007-1903（2016）01-0042-06

10.3969/j.issn.1007-1903.2016.01.009

郭帥（1989- ），女，碩士，主要從事地?zé)岬刭|(zhì)、水文地質(zhì)研究；E-mail：guoshuai8974@126.com。