，，

(1.山東省魯北地質(zhì)工程勘察院，山東 德州 253072；2.山東省武城縣國土資源局，山東 德州 253072)

當(dāng)前嚴(yán)峻的大氣污染形勢再次將地?zé)豳Y源的開發(fā)利用推到了前臺，作為一種綠色可再生能源，地?zé)崮荛_發(fā)成為了解決北方地區(qū)冬季清潔供暖問題的最佳選擇。武城縣作為山東省最早開展地?zé)峁┡牡貐^(qū)之一，具備得天獨厚的地?zé)豳Y源優(yōu)勢。本文旨在對武城縣地?zé)豳Y源的供暖開發(fā)利用前景進行評價，為武城縣的清潔供暖提供依據(jù)。

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

1.1 地質(zhì)構(gòu)造

武城縣在大地構(gòu)造上位于渤海灣盆地的南部，滄縣隆起與臨清坳陷的相接部位(見圖1)。渤海灣盆地的斷裂構(gòu)造走向以近EW向與NNE-NE向為主，斷裂性質(zhì)主要為正斷層[1]。在此背景下，區(qū)內(nèi)主要發(fā)育有NNE走向的滕莊-老城斷裂，EW走向的恩城斷裂。

滕莊-老城斷裂為滄東斷裂的南段，區(qū)域上為滄縣隆起與臨清坳陷的分界斷裂，斷層面呈上陡下緩的鏟狀特征，走向NE，為西盤上升、東盤下降的正斷裂[2]，傾角10°～80°。斷裂于印支-燕山期開始活動[3]，至喜馬拉雅運動第一幕[4](始新世)時活動強度達到最大，新近紀(jì)時活動消亡[5]，該斷裂具有較強的導(dǎo)水能力[6]，有利于深部高溫流體沿斷裂向淺部越流補給。

恩城斷裂呈近EW走向，為南升北降的正斷層，斷裂活動開始于古新世，始新世Es2+3時期活動最為強烈，新近紀(jì)時期活動消亡[7]。該斷裂僅局部切穿了館陶組的底部，館陶組下伏的東營組與沙河街組泥巖發(fā)育，在斷面附近常形成斷層泥，導(dǎo)致斷裂面密封性較好[8]，不利于下部流體對館陶組熱儲的補給。

1.2 地層

研究區(qū)地層分布特征受構(gòu)造的控制明顯，滕莊-老城斷裂以西的老城潛凸起區(qū)缺乏古近紀(jì)地層，該斷裂以東的德州潛凹陷區(qū)古近紀(jì)地層發(fā)育齊全(見圖2)。

古生界自下而上可劃分為下古生代寒武-奧陶系，巖性主要為碳酸鹽巖；上古生代石炭-二迭系，巖性為砂巖或泥巖。中生界主要發(fā)育有侏羅-白堊系，巖性為砂巖或泥巖。古近系自下而上可分為孔店組、沙河街組、東營組，巖性為砂巖或泥巖；其中沙河街組厚度巨大，根據(jù)其巖性特征進一步劃分為四段。新近系自下而上劃分為館陶組與明化鎮(zhèn)組，巖性為砂巖或泥巖，其中館陶組下段砂巖是武城縣地?zé)豳Y源開發(fā)的目標(biāo)熱儲。第四系為未成巖的粘性土或砂性土。

圖1 研究區(qū)構(gòu)造位置示意圖

圖2 武城縣地層空間分布剖面圖

1.3 巖漿巖

德州潛凹陷內(nèi)石油勘探發(fā)現(xiàn)，在孔店組與沙河街組四段中發(fā)育有多層厚度巨大的安山巖、玄武巖[9]。其中位于研究區(qū)西南邊界處的德參2井巖漿巖尤為發(fā)育，巖漿巖累計厚度近200 m。

2 地?zé)岬刭|(zhì)條件

區(qū)內(nèi)上古生代、中生代與古近紀(jì)地層中砂巖的巖性以粉砂巖為主，結(jié)構(gòu)致密，孔隙度低，當(dāng)前技術(shù)條件下開采成本較高，暫不作為地?zé)衢_發(fā)的目標(biāo)熱儲；老城潛凸起區(qū)下古生代寒武-奧陶系碳酸鹽巖埋深2 000 m左右，可作為下步開發(fā)的重點。新近紀(jì)館陶組熱儲最具開發(fā)利用價值[10-11]，也是區(qū)內(nèi)現(xiàn)有地?zé)峋哪繕?biāo)取水段，本次僅對該熱儲中地?zé)豳Y源的形成條件進行評述。

2.1 蓋層

區(qū)內(nèi)館陶組地層之上覆蓋有厚1 000～1 200 m的明化鎮(zhèn)組與第四系，巖性為泥巖或砂巖，熱導(dǎo)率低，成為了館陶組熱儲的良好蓋層。

2.2 熱源與地溫場特征

區(qū)內(nèi)地?zé)豳Y源類型為沉積盆地傳導(dǎo)型地?zé)豳Y源，熱源主要來自地球內(nèi)部的大地?zé)崃鱗12-13]。以往研究表明區(qū)內(nèi)平均大地?zé)崃髦禐?7.45 mw/m2[14]。地球內(nèi)部產(chǎn)生的熱量通過大地?zé)崃餍问皆丛床粩嗟叵蛲廨椛洌?jīng)歷漫長的地質(zhì)歷史時期在熱儲中淀積形成可供開發(fā)的地?zé)豳Y源。

以地?zé)峋膶崪y井口溫度與取水段平均埋深為參數(shù)，計算館陶組及其上覆地層的平均地溫梯度。結(jié)果表明，全區(qū)為地溫梯度>3.0℃/100m的地?zé)岙惓^(qū)(見圖3)，其中滕莊-老城斷裂具有明顯的導(dǎo)熱特征，沿斷裂分布有地溫梯度>3.4℃/100m的高異常帶。

圖3 武城縣館陶組及上覆地層平均地溫梯度等值線圖

圖4 館陶組熱儲分布與W18井身結(jié)構(gòu)圖表1 地?zé)峋疁囟扰c出水量一覽表

編號地點孔深/m溫度/℃降深/m出水量/m3/hW34魯騰家園1 411.715735.772W16利城金岸1 5916234.585W45張莊村1 574.46625890W20濱湖麗都1 558.4623880W19金馬商城1 5826145.495W1時代花園1 5706237.890

2.3 流體通道與熱儲特征

館陶組地層中分布有厚度巨大的砂巖、砂礫巖層，即是流體的通道，又是地?zé)崃黧w的儲集空間。根據(jù)鉆探資料，區(qū)內(nèi)館陶組熱儲可分為上下兩段，館上段熱儲巖性較細，以細砂巖、中細砂巖為主，單層厚度薄，砂巖與地層的厚度比25%左右，熱儲孔隙度約30%。館下段熱儲巖性較粗，以粗砂巖、含礫砂巖及砂礫巖為主，單層厚度大，砂巖與地層的厚度比50%左右，熱儲孔隙度約28%，該段熱儲是目前區(qū)內(nèi)地?zé)峋闹饕_采層段(見圖4)，熱儲溫度57℃～62℃，單井涌水量72～95 m3/h(見表1)。

圖5 館陶組熱儲底板埋深等值線圖

館陶熱儲底板埋深在德州潛凹陷區(qū)1 600～1 800 m，南深北淺、東深西淺；在老城潛凸起區(qū)1 450～1 500 m，南深北淺。其中滕莊-老城斷裂兩側(cè)館陶組底板落差約150 m左右，東深西淺。恩城斷裂兩側(cè)館陶組底板落差約50 m 左右，

北深南淺(見圖5)。館下段熱儲的累計厚度110～150 m，南厚北薄(見圖6)。

圖6 館陶組熱儲底板埋深等值線圖

3 地球化學(xué)特征

區(qū)內(nèi)館下段熱儲地?zé)崴嘘栯x子以Na+為主，當(dāng)量濃度百分比為82.73%～90.83%。陰離子以Cl-為主，當(dāng)量濃度百分比為70.62%～74.65%；其次為SO42-，當(dāng)量濃度百分比為20.61%～25.09%。水化學(xué)類型為Cl-Na或Cl·SO4-Na 型。地?zé)崴V化度4.67～4.81 g/L(見表2)。

表2 地?zé)崴饕x子含量一覽表

采用美國地調(diào)局免費軟件PHREEQC軟件對地?zé)崴兄饕V物的飽和指數(shù)進行計算[15]，結(jié)果表明，當(dāng)設(shè)定的反應(yīng)溫度接近熱儲溫度時，方解石的飽和指數(shù)達到“0”(見表3)，說明地?zé)崴兄饕x子濃度受方解石的水巖平衡反應(yīng)控制。

4 地?zé)豳Y源開發(fā)利用現(xiàn)狀

武城縣現(xiàn)有地?zé)峋?4眼(見圖7)，其中有25眼集中分布在武城縣城區(qū)，多建于2012年以前，由于當(dāng)時城市供暖基礎(chǔ)設(shè)施不完善，房地產(chǎn)商為解決新建住宅冬季供暖問題建設(shè)的。其余19眼地?zé)峋疄榻鉀Q各鄉(xiāng)鎮(zhèn)新建農(nóng)村社區(qū)供暖問題而建設(shè)的。2016年底，鑒于地?zé)峁┡菜欧诺沫h(huán)保壓力，武城縣對城市集中供暖設(shè)施可以輻射到的地?zé)峋畬嵭辛岁P(guān)停，城區(qū)內(nèi)僅保留了對地?zé)峁┡菜M行了回灌的6眼地?zé)峋?，農(nóng)村社區(qū)的19眼地?zé)峋笙奁陂_展地?zé)嵛菜毓唷?/p>

5 地?zé)豳Y源評價

5.1 地?zé)豳Y源量

本次采用《地?zé)豳Y源地質(zhì)勘查規(guī)范》[16]中推薦的熱儲法對館下段熱儲中蘊藏的熱資源量進行計算，由于區(qū)內(nèi)主要開采館下段熱儲中地?zé)崴M行住宅供暖，計算的溫度下限取地?zé)峁┡菜疁囟?0℃。通過非穩(wěn)定流抽水試驗求取的館下段熱儲彈性釋水系數(shù)約為4.2×10-4，比熱儲的孔隙度小2～3個數(shù)量級，流體彈性儲存量中蘊含的熱資源量在計算中予以忽略，計算公式如下：

Q=A·d·[ρs·Cs·(1-φ)+ρw·Cw·φ]·(tr-t0)

(1)

式中：Q為熱資源量(J)；A為熱儲面積(m2)；d為熱儲厚(m)；ρs為熱儲巖石密度(kg/m3)；Cs為巖石的比熱[J/(kg·℃)]；φ為熱儲孔隙度，無量綱；ρw為地?zé)崴芏?kg/m3)；Cw為地?zé)崴谋葻醄J/(kg·℃)]；tr為熱儲平均溫度(℃)；t0為溫度下限(℃)。

經(jīng)計算，武城縣館下段熱儲中儲存的熱資源量為8.11×1018J，折合標(biāo)準(zhǔn)煤2.77億 t。

5.2 地?zé)崴o儲存量

熱儲中地?zé)崴o儲存量包括熱儲孔隙中容積儲量與高壓下地?zé)崴詮椥詨好苄问酱嬖诘膹椥詢α?，計算公式如下?/p>

Q=A·d·φ+A·μe·H

(2)

式中：Q為地?zé)崴o儲存量(m3)；μe為熱儲彈性貯水系數(shù)，數(shù)值等于熱儲的新熱的彈性釋水系數(shù)，無量綱；H為熱儲中點至地?zé)崴杂伤娴木嚯x(m)。

經(jīng)計算，武城縣館下段熱儲中蘊藏的地?zé)崴o儲量為2.75×1010m3。

5.3 地?zé)崴|(zhì)量評價

區(qū)內(nèi)館下段熱儲地?zé)崴饕糜诮ㄖ锕┡?，針對供暖需求進行地?zé)崴慕Y(jié)垢性與腐蝕性評價。因地?zé)崴新入x子含量高，采用拉申指數(shù)判斷地?zé)崴妓徕}結(jié)垢趨勢與腐蝕性，公式如下：

(3)

式中：LI為拉申指數(shù)，無量綱；Cl為氯離子濃度(mg/L)；SO4為硫酸根離子濃度(mg/L)；ALK為總堿度(mg/L )。經(jīng)計算，LI為12.61～14.86，地?zé)崴唤Y(jié)垢，具有強腐蝕性。

表3 不同溫度下地?zé)崴兄饕V物飽和指數(shù)一覽表

圖7 武城縣地?zé)峋植紙D

6 地?zé)豳Y源開發(fā)利用前景

2016年，武城縣對城區(qū)的6眼地?zé)峋畬嵤┝伺鋵毓嗑責(zé)崴墓嗖杀葹?3.98%～97.94%，證明了地?zé)峁┡菜毓嗍乔袑嵖尚械摹Ｔ诨毓?開采的地?zé)豳Y源利用方案下，將供暖尾水全部回注入熱儲中，只取熱不取水，地?zé)豳Y源的開發(fā)利用前景主要受熱儲中蘊含的熱資源量限制。根據(jù)中低溫地?zé)豳Y源開發(fā)利用年限100年計算，武城縣館下段熱儲中蘊藏的年可利用熱資源量為8.11×1016 J/a。以單位面積住宅供暖定額50 W估算，館下段熱儲可滿足1.56×108m2的建筑面積供暖需求。

國家統(tǒng)計局公布的2016年全國人均居住面積40.8 m2，武城縣2010年統(tǒng)計人口為390 412，人口自然增長率7.98‰，推測2017年武城縣人口約為41.27萬人，則住宅的總供暖面積需求約為1 688萬 m2。100 a后武城縣預(yù)計人口數(shù)為86.44萬人，供暖總需求約為3 527萬 m2。由此可見，館下段熱儲中蘊含的熱資源量完全可滿足未來100 a內(nèi)武城縣居民住宅的供暖需求。

7 結(jié)語

(1)武城縣在地質(zhì)構(gòu)造上可劃分為老城潛凸起與德州潛凹陷兩個構(gòu)造單元，館陶熱儲底板埋深在德州潛凹陷區(qū)1 600～1 800 m，南深北淺、東深西淺；在老城潛凸起區(qū)1 450～1 500 m，南深北淺。

(2)館陶組熱儲可劃分為上、下兩段，其中館下段熱儲是武城縣地?zé)豳Y源開采的主要目標(biāo)熱儲，熱儲巖性以粗砂巖、含礫砂巖及砂礫巖為主，熱儲厚110～150 m；熱儲溫度57℃～62℃，單井涌水量72～95 m3/h。

(3)武城縣館下段熱儲中儲存的熱資源量為8.11×1018J，折合標(biāo)準(zhǔn)煤2.77億 t；蘊藏的地?zé)崴o儲量為2.75×1010m3。

(4)武城縣館下段熱儲可支持建筑供暖面積1.56×108m2，完全可滿足未來100年內(nèi)武城縣居民住宅的供暖需求。

[1]漆家福, 楊橋, 陸克政,等. 渤海灣盆地基巖地質(zhì)圖及其所包含的構(gòu)造運動信息[J]. 地學(xué)前緣.2004.11(3):299-307.

[2]高戰(zhàn)武, 徐杰, 宋長青. 華北滄東斷裂的構(gòu)造特征[J]. 地震地質(zhì).2000.22(4):395-404.

[3]地質(zhì)礦產(chǎn)部地質(zhì)專報.天津市區(qū)域地質(zhì)志[M].北京：地質(zhì)出版社.1982.198-199.

[4]賈承造, 何登發(fā), 陸潔民. 中國喜馬拉雅運動的期次及其動力學(xué)背景[J]. 石油與天然氣地質(zhì).2004.25(2):121-125.

[5]侯旭波. 臨清坳陷東部中、新生代盆地演化分析[D]. 中國石油大學(xué).2007.

[6]趙蘇民, 高寶珠, 黎雪梅,等. 滄東斷裂(天津段)特征及導(dǎo)水導(dǎo)熱性質(zhì)分析[J]. 地質(zhì)調(diào)查與研究.2007.30(2):121-127.

[7]張存霞. 臨清坳陷東部地區(qū)地層格架及展布特征[D]. 中國石油大學(xué)(華東).2008.

[8]陳清華, 王紹蘭, 吳孔友. 臨清坳陷莘縣凹陷油氣勘探遠景評價[J]. 中國石油大學(xué)學(xué)報:自然科學(xué)版.1998(2):8-10.

[9]劉桂儀, 孟慶蜂. 德州市低溫地?zé)豳Y源及開發(fā)利用研究[C]. “九五”全國地質(zhì)科技重要成果學(xué)術(shù)交流會. 2000.

[10]王彥俊, 劉桂儀, 胡松濤. 魯北地區(qū)地?zé)豳Y源區(qū)劃研究[J]. 地質(zhì)調(diào)查與研究.2008.31(3):5-12.

[11]王彥俊, 王貞國, 王巖,等. 德州市地?zé)豳Y源開發(fā)與保護[J]. 山東國土資源.2005.21(5):31-35.

[12]趙書泉, 周紹智, 鄒祖光,等. 山東省魯北地?zé)崽锛捌淇沙掷m(xù)開發(fā)利用對策[C]. 全國地?zé)岙a(chǎn)業(yè)可持續(xù)發(fā)展學(xué)術(shù)研討會.2005.

[13]陳墨香. 華北地?zé)醄M]. 科學(xué)出版社.1988.

[14]吳乾蕃, 謝毅真, 祖金華,等. 華北地?zé)釄鲅芯縖J]. 中國地震.1988(1):43-50.

[15]趙季初, 高宗軍. 水-巖平衡在館陶組熱儲溫度估算中的應(yīng)用[J]. 地下水.2015(5):19-21.

[16]GB/T 11615-2010.地?zé)豳Y源地質(zhì)勘查規(guī)范[S].