張 萌，秦祥熙，黃 雯，石 磊，鄭 鋼

1. 浙江陸特能源科技股份有限公司，浙江 杭州 310051；2. 東華大學(xué)，上海 201620；3. 百合花集團(tuán)股份有限公司，浙江 杭州 311228；4. 浙江省國(guó)土整治中心，浙江 杭州 310012；5. 滄州陸特新能源科技有限公司，河北 滄州 061000

0 引言

河南省地?zé)豳Y源稟賦條件良好，開(kāi)采潛力較大，主要分布在以濟(jì)源-開(kāi)封凹陷為代表的東部平原區(qū)［1-2］.濟(jì)源-開(kāi)封凹陷地?zé)豳Y源的開(kāi)發(fā)利用由來(lái)已久，眾多學(xué)者［3-8］已分別對(duì)其內(nèi)的武陟、原陽(yáng)、鄭州、中牟等地區(qū)，開(kāi)展了地?zé)豳Y源的熱儲(chǔ)特征、賦存情況、開(kāi)發(fā)利用等研究. 在其東段的開(kāi)封-蘭考地區(qū)內(nèi)，相關(guān)研究多集中于開(kāi)封市區(qū)及周邊［9-11］，且受限于以往的開(kāi)采深度，又多針對(duì)明化鎮(zhèn)組熱儲(chǔ). 近年來(lái)，該區(qū)的地?zé)豳Y源開(kāi)發(fā)向更深部的館陶組熱儲(chǔ)延伸，開(kāi)發(fā)利用規(guī)模進(jìn)一步增加. 目前，在開(kāi)封市、蘭考縣境內(nèi)已建設(shè)多個(gè)館陶組熱水區(qū)域集中供熱項(xiàng)目，圍繞蘭考境內(nèi)館陶組熱儲(chǔ)層的賦存規(guī)律、資源評(píng)價(jià)等研究也相繼開(kāi)展［12-14］，但是關(guān)于館陶組地?zé)崴男纬蓷l件與成因機(jī)制的研究尚存不足.

本文以濟(jì)源-開(kāi)封凹陷東段的開(kāi)封斷隆與開(kāi)封-蘭考斷陷（涉及的行政區(qū)大致為開(kāi)封市區(qū)及西部至蘭考縣城區(qū)）為研究區(qū)，利用2017—2020 年在蘭考、開(kāi)封投資建設(shè)的水熱型地?zé)崮軈^(qū)域集中供熱項(xiàng)目所施工的館陶組地?zé)峋畔?，以館陶組熱儲(chǔ)為研究對(duì)象，通過(guò)開(kāi)展地下水的化學(xué)特征、形成作用等水文地球化學(xué)研究，探討其賦存環(huán)境及成因特點(diǎn)，以期為該地區(qū)內(nèi)地?zé)豳Y源的開(kāi)發(fā)利用提供科學(xué)依據(jù).

1 區(qū)域地?zé)岬刭|(zhì)條件

濟(jì)源-開(kāi)封凹陷（Ⅰ44）屬于中朝準(zhǔn)地臺(tái)（Ⅰ）之華北拗陷（Ⅰ4），沿黃河呈東西向展布于濟(jì)源-開(kāi)封-民權(quán)一帶. 區(qū)內(nèi)斷裂發(fā)育，近東西向的鄭州-開(kāi)封斷裂（F2）為南部邊界，古近紀(jì)時(shí)期為強(qiáng)活動(dòng)正斷層，造成區(qū)內(nèi)巨厚的新生界沉積，至第四紀(jì)已停止活動(dòng)［15］. 北西向的濟(jì)源-商丘深大斷裂（F1）為北部邊界，具有多期次活動(dòng)特征，且切割較深，為深部地下熱能的傳導(dǎo)及水源、熱源的流通創(chuàng)造了有利條件［16］. 一系列的北北東向斷裂自西向東將濟(jì)源-開(kāi)封凹陷分割為濟(jì)源凹陷、武陟斷凸、原陽(yáng)-封丘斷陷、開(kāi)封斷隆及開(kāi)封-蘭考斷陷等多個(gè)斷塊（圖1）［17］. 該區(qū)屬華北地層區(qū)華北平原分區(qū)，基底主要為奧陶系、石炭系、二疊系、侏羅系及白堊系，新生代地層底板埋深2 000～5 000 m，自下而上為古近系、新近系、第四系.

圖1 研究區(qū)與地?zé)峋恢檬疽鈭D（修改自文獻(xiàn)［17-18］）Fig. 1 Location sketch map of the study area and geothermal wells（Modified from References［17-18］）

區(qū)內(nèi)地?zé)豳Y源豐富，源、通、蓋、儲(chǔ)結(jié)構(gòu)完整. 熱源供給為大地?zé)崃鱾鲗?dǎo)；區(qū)域性深大斷裂為深部熱源的溝通傳導(dǎo)與地?zé)崴倪\(yùn)移創(chuàng)造條件；淺部地層發(fā)育有多層黏土和砂質(zhì)黏土層，具有保溫隔熱作用，形成良好的蓋層；新生界地層厚度大、孔隙度高、富水性好，是良好的熱儲(chǔ)層，從上自下依次發(fā)育新近系明化鎮(zhèn)組熱儲(chǔ)層、新近系館陶組熱儲(chǔ)層和古近系熱儲(chǔ)層，各熱儲(chǔ)層相互獨(dú)立，其間分布著弱透水層隔絕上下的連通［12］.

2 館陶組熱儲(chǔ)基本特征

研究區(qū)包括濟(jì)源-開(kāi)封凹陷東段的開(kāi)封斷隆、開(kāi)封-蘭考斷陷（圖1），區(qū)內(nèi)新近系地層整體自西向東傾斜，除部分地區(qū)因人工開(kāi)采造成局部降落漏斗外，區(qū)域上地?zé)崴晕飨驏|徑流［19］. 館陶組地層頂板埋深在開(kāi)封斷隆內(nèi)自南向西北逐漸變深，開(kāi)封市南部約1 250 m，北部和西部可達(dá)1 700 m. 館陶組厚度從南到北約為200～800 m，蘭考縣境1 300～1 600 m 不等，平均約1 500 m，厚度總體呈東薄西厚、北薄南厚的特點(diǎn)［12，19］. 本研究共選用地?zé)峋?0 口（表1），其中開(kāi)封斷隆2 口、開(kāi)封-蘭考斷陷8 口，地?zé)峋罹┩葛^陶組地層，取水層僅為館陶組含水層.

表1 館陶組地?zé)峋厩闆rTable 1 Basic information of geothermal wells in Guantao Formation

開(kāi)封斷隆內(nèi)館陶組巖性上部為黏土夾細(xì)砂、粉細(xì)砂，下部為泥巖夾細(xì)粒砂巖、粉砂巖，頂板埋深1 450～1 520 m，地層厚度692.48～881.36 m. 開(kāi)封-蘭考斷陷內(nèi)館陶組地層巖性上部為細(xì)-中砂巖及泥巖與細(xì)砂巖互層夾粉砂巖，下部巖性為礫巖、含礫砂巖、砂礫巖、細(xì)砂巖夾粉砂巖、泥巖，頂板埋深1 260.71～1 600 m，地層厚度330.93～672.00 m. 圖2 顯示，兩者館陶組熱儲(chǔ)含水層單層厚度差異不大，但開(kāi)封-蘭考斷陷內(nèi)易出現(xiàn)異常大值. 而兩者的含水層孔隙度差別明顯，開(kāi)封斷隆內(nèi)數(shù)值相對(duì)集中且偏小，大部分集中在10.13%～12.41%區(qū)間，平均11.43%，變化幅度也不大；開(kāi)封-蘭考斷陷內(nèi)大部分?jǐn)?shù)值集中在25.72%～41.27%，變化幅度相對(duì)較大，但各地?zé)峋钠骄迪嗖畈欢啵瑸?1.98%～35.91%，整體平均值為34.28%.

圖2 不同區(qū)域館陶組熱儲(chǔ)含水層單層厚度與孔隙度對(duì)比Fig.2 Comparison of thickness and porosity of geothermal reservoir single aquifer in Guantao Formation of different areas

地?zé)峋疁y(cè)溫結(jié)果（圖3）顯示，垂向上研究區(qū)地溫隨深度加深而逐漸升高，呈現(xiàn)典型的傳導(dǎo)型地溫變化特點(diǎn). 平面上，開(kāi)封斷隆內(nèi)館陶組熱儲(chǔ)的蓋層地溫梯度為2.79～2.92 ℃/hm；開(kāi)封-蘭考斷陷內(nèi)為3.25～4.06 ℃/hm，平均3.45 ℃/hm，高于前者（表1）. 開(kāi)封斷隆內(nèi)館陶組熱儲(chǔ)地?zé)崴o水位埋深58.1～78.9 m，井口水溫82～84 ℃，出水量80～96 m3/h，單位涌水量2.08～2.21 m3（/h·m）；開(kāi)封-蘭考斷陷內(nèi)靜水位埋深38.5～55 m，井口水溫70～76 ℃，出水量97～131 m3/h，單位涌水量4.13～11.82 m3（/h·m）.

圖3 研究區(qū)地?zé)峋疁y(cè)溫曲線Fig. 3 Temperature logging curves of geothermal wells in the study area

綜合以上分析可知，開(kāi)封斷隆與開(kāi)封-蘭考斷陷相較，開(kāi)封斷隆內(nèi)館陶組熱儲(chǔ)層埋藏較深，厚度較大，含水層孔隙度較低，抽水降深較大，單井出水量和單位涌水量較低，出水溫度較高.

3 水化學(xué)特征及形成條件

3.1 水質(zhì)分析結(jié)果

研究區(qū)館陶組地?zé)崴Ａ拷M分的測(cè)試結(jié)果見(jiàn)表2，同時(shí)選取明化鎮(zhèn)組地?zé)崴?、第四系地下水以及地表（黃河）水，用于對(duì)比分析.

表2 地?zé)崴饕x子測(cè)試結(jié)果Table 2 Test results of main ions in geothermal water

由表2 可知，研究區(qū)館陶組地?zé)崴畃H 值均介于6.5～8 之間，屬于中性水；明化鎮(zhèn)組地?zé)崴畃H 則略大于8；第四系地下水及地表水也屬于中性水. 館陶組地?zé)崴目側(cè)芙夤腆w（TDS）介于10～50 g/L，均為鹽水，總硬度（以CaCO3計(jì)）除LK1 外均大于450 mg/L，為極硬水，KF1 的總硬度甚至高達(dá)10 969 mg/L. 開(kāi)封斷隆內(nèi)明化鎮(zhèn)組地?zé)崴偷谒南档叵滤甌DS＜1 g/L，為淡水，總硬度小于150 mg/L，為軟水. 開(kāi)封-蘭考斷陷內(nèi)明化鎮(zhèn)組地?zé)崴畡t為微咸水、軟水，第四系地下水為微咸水、硬水. 地表水則為淡水、硬水.

Schoeller 圖（圖4）顯示，開(kāi)封斷隆與開(kāi)封-蘭考斷陷內(nèi)館陶組地?zé)崴Ａ拷M分的變化趨勢(shì)相同，而與明化鎮(zhèn)組熱水、第四系地下水及地表水相異. Piper 三線圖（圖5）表明館陶組熱水化學(xué)類型明顯異于其他地下水和地表水. 除KF1 外，館陶組熱水水化學(xué)類型均為Cl-Na 型，Cl-、Na+占有絕對(duì)優(yōu)勢(shì)，分別為首要陰離子和首要陽(yáng)離子. KF1 總硬度（以CaCO3計(jì)）與其他館陶組熱水水樣數(shù)值相差一個(gè)數(shù)量級(jí)，水中含有大量的Ca2+，因此Ca2+成為次要陽(yáng)離子，水化學(xué)類型表現(xiàn)為Cl-Na·Ca 型. 明化鎮(zhèn)組熱水中，Na+仍為首要陽(yáng)離子，陰離子組分中HCO3-占比較大，在開(kāi)封斷隆內(nèi)為首要陰離子，在開(kāi)封-蘭考斷陷內(nèi)為次要陰離子. 第四系地下水中Cl-占比均低于25%，不表現(xiàn)在水化學(xué)類型中，Ca2+、Mg2+在陽(yáng)離子中的占比較高，超過(guò)25%. 地表水的水化學(xué)類型則為HCO3·SO4·Cl-Na·Ca 型.由以上分析可知，研究區(qū)館陶組地?zé)崴乃|(zhì)特點(diǎn)與化學(xué)類型與上部的明化鎮(zhèn)組熱水、第四系地下水及地表水差異明顯，表明它們之間的水力聯(lián)系微弱，所處環(huán)境封閉［23-24］. 而開(kāi)封斷隆與開(kāi)封-蘭考斷陷相較，前者的館陶組地?zé)崴哂懈叩腡DS 值與硬度值，除KF1因Ca2+含量較高表現(xiàn)為Cl-Na·Ca 型水外，其余均為Cl-Na 型水.

圖4 研究區(qū)水樣Schoeller 圖Fig. 4 Schoeller diagram of water samples from the study area

圖5 研究區(qū)水樣Piper 三線圖Fig. 5 Piper diagrams of water samples from the study area

3.2 水化學(xué)控制因素分析

Gibbs 圖［25-26］（圖6）顯示研究區(qū)館陶組地?zé)崴奂谡舭l(fā)濃縮作用區(qū)域，遠(yuǎn)離大氣降水控制區(qū)域，表明其控制機(jī)制以蒸發(fā)濃縮作用為主導(dǎo)，受大氣降水作用影響小，可能還存在巖石風(fēng)化作用的影響. Mg2+/Na+與Ca2+/Na+離子比值端元圖［27-28］（圖7）中同樣展現(xiàn)出這樣的特征，研究區(qū)館陶組地?zé)崴饕杏谡舭l(fā)鹽巖溶解區(qū)域，部分向硅酸鹽風(fēng)化區(qū)域偏移. 因此，研究區(qū)館陶組熱水主要受鹽巖的蒸發(fā)濃縮作用影響，開(kāi)封斷隆區(qū)內(nèi)可能還受巖石風(fēng)化作用影響，在開(kāi)封-蘭考斷陷內(nèi)鹽巖的蒸發(fā)濃縮作用對(duì)館陶組熱水的控制影響大于開(kāi)封斷隆區(qū). 蒸發(fā)濃縮作用會(huì)導(dǎo)致水中溶解度低的鹽類離子（如Ca2+、Mg2+、HCO3-、CO32-）濃度降低，溶解度高的鹽類離子（如Na+、Cl-、SO42-）濃度升高［29］，因此研究區(qū)館陶組熱儲(chǔ)地?zé)崴故境銮拔乃龅乃瘜W(xué)類型. 明化鎮(zhèn)組熱儲(chǔ)、第四系地下水和地表水在圖6、7中明顯位于與館陶組不同的區(qū)域，表明它們與館陶組熱儲(chǔ)有著不同的控制因素.

圖6 研究區(qū)水樣Gibbs 圖Fig. 6 Gibbs diagrams of water samples from the study area

圖7 Mg2+/Na+與Ca2+/Na+比值關(guān)系Fig. 7 Relationship between Mg2+/Na+and Ca2+/Na+

地下水的封閉程度、變質(zhì)程度、濃縮程度可用變質(zhì)系數(shù)（γNa/γCl）、脫硫系數(shù)（100×γSO4/γCl）、鹽化系數(shù)［γCl/（γHCO3+γCO3）］等特征系數(shù)來(lái)反映［30］. 研究區(qū)地下水的特征系數(shù)值如表3 所示. 研究區(qū)館陶組熱水的變質(zhì)系數(shù)均小于1，部分接近1，表現(xiàn)出含鹽巖地層溶濾水的特征；且小于明化鎮(zhèn)組地?zé)崴偷谒南档叵滤砻髌涮幱诟忾]、更停滯的賦存環(huán)境，水體變質(zhì)程度更深，水體為還原環(huán)境. 然而，部分館陶組熱水的變質(zhì)系數(shù)則小于標(biāo)準(zhǔn)海水γNa/γCl 的平均值（0.85），展現(xiàn)出海相沉積水的特點(diǎn)［31］，但其Cl/Br 系數(shù)均大于293（本研究?jī)H部分館陶組地?zé)峋畼臃治隽薆r 的含量，其中KF1 Cl/Br 系數(shù)為2958.77、KF1 為1 511.84、LK1 為566.12、LK3 為709），可排除其海相沉積水的成因，應(yīng)屬于陸相含鹽沉積型地下水［32-33］.

表3 研究區(qū)地下水特征系數(shù)Table 3 Characteristic coefficients of groundwater in the study area

研究區(qū)館陶組熱水的脫硫系數(shù)遠(yuǎn)小于明化鎮(zhèn)組地?zé)崴偷谒南档叵滤?，相?～2 個(gè)數(shù)量級(jí)，說(shuō)明與其他地下水相比，館陶組熱水熱儲(chǔ)封閉性良好，熱水還原較徹底. 而開(kāi)封斷隆內(nèi)館陶組熱水的脫硫系數(shù)小于1［34］，表明其環(huán)境更封閉、還原更徹底；開(kāi)封-蘭考斷陷內(nèi)的館陶組熱水大于1，則表明其可能受到淺表氧化作用的影響.

研究區(qū)水樣的鹽化系數(shù)與其TDS 含量呈現(xiàn)良好的正相關(guān)關(guān)系（圖8），即鹽化系數(shù)越高，水的TDS 越大，水的濃縮程度越高［35］. 館陶組地?zé)崴柠}化系數(shù)遠(yuǎn)高于明化鎮(zhèn)組地?zé)崴⒌谒南档叵滤?，表明其具有較強(qiáng)的濃縮程度，且開(kāi)封斷隆區(qū)內(nèi)館陶組地?zé)崴c開(kāi)封-蘭考斷陷區(qū)相比更濃縮.

圖8 鹽化系數(shù)與TDS 關(guān)系圖Fig. 8 Relationship between salinization coefficient and TDS

4 水化學(xué)特征對(duì)成因的啟示

已有的同位素研究成果表明，研究區(qū)內(nèi)地?zé)崴鹪从诖髿饨邓?，且年代久遠(yuǎn)［36］. 區(qū)內(nèi)地?zé)豳Y源類型為沉積盆地傳導(dǎo)型，地?zé)崴纱蟮責(zé)崃鱾鲗?dǎo)增溫. 發(fā)育的斷裂、裂隙，特別是區(qū)域性深大斷裂，為地?zé)崴膶?dǎo)熱和運(yùn)移提供了有利條件. 區(qū)內(nèi)發(fā)育多層熱儲(chǔ)，且垂向上聯(lián)系甚微，上部熱儲(chǔ)可作為下部熱儲(chǔ)層的良好蓋層.

研究區(qū)館陶組熱水賦存于砂巖孔隙中，在漫長(zhǎng)的地質(zhì)歷史時(shí)期不斷與周圍地層發(fā)生溶濾作用，屬陸相含鹽巖地層沉積型溶濾水. 館陶組地層封閉性良好，地?zé)崴幕瘜W(xué)特征、控制因素、特征系數(shù)等與上部明化鎮(zhèn)組地?zé)崴⒌谒南档叵滤偷乇硭耆煌?，進(jìn)一步佐證了區(qū)內(nèi)各層地下水間基本無(wú)水力聯(lián)系. 館陶組熱水化學(xué)特征的控制機(jī)制以蒸發(fā)濃縮作用為主導(dǎo)，館陶組熱水賦存環(huán)境封閉、停滯，水體還原程度高，濃縮程度大，水中聚集大量的Na+、Cl-，水化學(xué)類型總體表現(xiàn)為Cl-Na 型. 而開(kāi)封斷隆與開(kāi)封-蘭考斷陷對(duì)比，其館陶組熱儲(chǔ)含水層孔隙度較低、單井出水量、單位涌水量較低，其熱水的賦存環(huán)境更封閉、更停滯，熱水還原更徹底、濃縮程度更高，熱水TDS 和硬度相對(duì)更高.

5 結(jié)論

通過(guò)對(duì)濟(jì)源-開(kāi)封凹陷東段的開(kāi)封斷隆、開(kāi)封-蘭考斷陷內(nèi)館陶組熱儲(chǔ)層的基本特點(diǎn)、熱水水化學(xué)特征及形成條件進(jìn)行分析研究，得出主要結(jié)論如下：

（1）研究區(qū)館陶組熱水為中性水、鹽水、極硬水，水化學(xué)類型以Cl-Na 型為主. 其中開(kāi)封斷隆與開(kāi)封-蘭考斷陷相較，前者館陶組熱儲(chǔ)層埋藏較深，厚度較大，含水層孔隙度較低，抽水降深較大，單井出水量和單位涌水量較低，出水溫度較高.

（2）研究區(qū)館陶組熱水為陸相含鹽巖地層沉積型溶濾水，賦存介質(zhì)為砂巖孔隙. 區(qū)域發(fā)育的斷裂、裂隙為其提供導(dǎo)熱和運(yùn)移的通道，在大地?zé)崃鞯膫鲗?dǎo)增溫下，不斷發(fā)生溶濾作用、蒸發(fā)濃縮作用等. 同時(shí)熱儲(chǔ)層封閉性良好，熱水運(yùn)移較停滯，使得熱水還原性徹底，濃縮程度高. 館陶組熱水與上部明化鎮(zhèn)組熱水、第四系地下水和地表水的水化學(xué)特征、控制因素、特征系數(shù)等完全不同，垂向間基本無(wú)水力聯(lián)系. 開(kāi)封斷隆與開(kāi)封-蘭考斷陷對(duì)比，其館陶組熱水的賦存環(huán)境更封閉、更停滯，熱水還原更徹底、濃縮程度更高，熱水的TDS和硬度相對(duì)更高.

（3）本文的研究分析結(jié)果指示研究區(qū)內(nèi)的館陶組熱水再生能力差，屬消耗型資源，在開(kāi)采過(guò)程中易出現(xiàn)結(jié)垢、堵塞、回灌能力衰減等問(wèn)題，近年來(lái)的實(shí)際工程項(xiàng)目中確實(shí)已出現(xiàn)上述問(wèn)題. 為維持地?zé)崴目沙掷m(xù)開(kāi)發(fā)，完成100%同層回灌，目前已采取的應(yīng)對(duì)措施中最有效的是增設(shè)回灌井，但地?zé)衢_(kāi)發(fā)企業(yè)的投資建造成本、運(yùn)行維護(hù)成本也同步增加. 這些都對(duì)研究區(qū)館陶組地?zé)崴母咝?yīng)用、可持續(xù)開(kāi)發(fā)、吸引市場(chǎng)資金投入等方面造成了挑戰(zhàn). 建議今后結(jié)合本文的研究成果，在地?zé)崴_(kāi)發(fā)利用中做到合理規(guī)劃、科學(xué)有序，同時(shí)加大對(duì)地?zé)崴氯麢C(jī)理、防垢阻垢、增加熱儲(chǔ)層回灌能力等方面的工程技術(shù)研究，促進(jìn)地?zé)豳Y源經(jīng)濟(jì)可持續(xù)開(kāi)發(fā).