吳清華，馮穎，劉帥，張靜，王玲麗

(山東省地勘局第二水文地質(zhì)工程地質(zhì)大隊(duì)(山東省魯北地質(zhì)工程勘察院)，山東 德州 253015)

0 引言

莘縣地區(qū)地?zé)豳Y源豐富，開(kāi)采條件好，開(kāi)發(fā)潛力大，利用前景廣闊。正在進(jìn)行的高鐵新城項(xiàng)目是莘縣正在重點(diǎn)推進(jìn)的兩大戰(zhàn)略項(xiàng)目之一，項(xiàng)目涉及約400hm2基本農(nóng)田需調(diào)整為一般農(nóng)田后轉(zhuǎn)為建設(shè)用地，會(huì)建設(shè)大量的公共、商業(yè)、辦公等建筑，因此查清莘縣地區(qū)地?zé)豳Y源的賦存條件并對(duì)地?zé)豳Y源進(jìn)行計(jì)算和評(píng)價(jià)，對(duì)解決新城區(qū)供暖問(wèn)題、優(yōu)化能源結(jié)構(gòu)、減少污染物排放具有十分重要意義。

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

1.1 地層

研究區(qū)被第四系覆蓋，自中新世以來(lái)，受差異性升降運(yùn)動(dòng)的影響，一直緩慢下沉，沉積了厚大新生代地層。地層由老至新有古近系、新近系和第四系(圖1)。

1—古近紀(jì)東營(yíng)組；2—古近紀(jì)孔店-沙河街組；3—侏羅紀(jì)三臺(tái)組；4—二疊紀(jì)石盒子群； 5—二疊紀(jì)山西組；6—石炭-二疊紀(jì)太原組；7—奧陶紀(jì)馬家溝群；8—物探解譯斷層；9—隱伏平行不整合線；10—研究區(qū)范圍圖1 莘縣地區(qū)區(qū)域地質(zhì)圖

(1)古近系

古近系包括孔店組、沙河街組和東營(yíng)組?？椎杲M主要為砂巖與泥巖互層，底部為礫巖，厚度大于400m；沙河街組巖性主要為灰?guī)r、泥巖和泥砂巖，底部為泥巖及砂礫巖，厚度145～789m；東營(yíng)組巖性主要為砂巖、含礫砂巖、泥巖和細(xì)礫巖，與下伏沙河街組呈整合接觸，厚度100～500m左右。

(2)新近系

新近系包括明化鎮(zhèn)組和館陶組。館陶組巖性主要為細(xì)砂巖、中砂巖、砂礫巖和泥巖，砂巖及砂礫巖分選性較差，磨圓度中等，膠結(jié)性差，與下伏東營(yíng)組呈不整合接觸，厚度300～500m；明化鎮(zhèn)組巖性主要為細(xì)砂巖和泥巖為主，與下伏館陶組呈整合接觸，厚度600～1000m。

(3)第四系

以平原組為主，局部地表可能有黃河組。巖性主要為粉砂、粉細(xì)砂、細(xì)砂、粉質(zhì)黏土、砂質(zhì)黏土為主，下部層底結(jié)構(gòu)致密，含大量鈣質(zhì)、鐵、錳質(zhì)結(jié)核以及鈣質(zhì)沉淀層。厚度120～290m，與下伏新近紀(jì)明化鎮(zhèn)組呈不整合接觸[1]。

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

研究區(qū)大地構(gòu)造位于華北板塊(Ⅰ級(jí))華北凹陷區(qū)(Ⅱ級(jí))臨清凹陷(Ⅲ級(jí))東明-莘縣潛斷陷(Ⅳ級(jí))莘縣潛凹陷(Ⅴ級(jí))范圍內(nèi)。

區(qū)內(nèi)重要斷裂為聊考斷裂，聊考斷裂南起河南蘭考，北至山東禹城，全長(zhǎng)大于300km，走向20°～40°，傾向NW，傾角60°左右，其北端落差變小漸而消失。斷裂帶寬4～8km，落差大于2500m，它控制了中生代以來(lái)的地質(zhì)發(fā)展史，在重力異常圖上表現(xiàn)為明顯的NNE向重力密集帶，電測(cè)深剖面上也有明顯的反映，范縣至聊城，兩側(cè)電法標(biāo)志層(奧陶系頂面)埋深相差2000～3000m。聊考斷裂是一個(gè)新構(gòu)造活動(dòng)帶，1502年—1948年間僅5級(jí)以上地震就發(fā)生過(guò)5次。

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

莘縣地?zé)衢_(kāi)采區(qū)處于聊城西部地?zé)崽?，熱?chǔ)類(lèi)型為層狀孔隙型，主要開(kāi)采新近紀(jì)館陶組砂巖熱儲(chǔ)和古近紀(jì)東營(yíng)組砂巖熱儲(chǔ)。

2.1 熱儲(chǔ)特征

(1)館陶組

該熱儲(chǔ)受區(qū)域構(gòu)造和基底起伏的控制，頂板埋藏深、厚度大，自東向西由淺變深，厚度由小變大(圖2)。頂板埋深約1000m，底板埋深1300～1500m，地層厚度300～500m，含水層累計(jì)厚度155～175m。砂層巖性主要為河流相、沖積扇相的細(xì)砂巖、粗砂巖、含礫砂巖、砂礫巖，礫石呈半圓狀，磨圓度中等。在1000～1500m取水段，單井涌水量1500～2000m3/d。地?zé)崴V化度4～6g/L，水化學(xué)類(lèi)型為Cl-Na和SO4·Cl-Na型，井口出水溫度45～60℃，屬溫?zé)嵝汀獰崴偷蜏氐責(zé)豳Y源。

1—第四系；2—新近紀(jì)明化鎮(zhèn)組；3—新近紀(jì)館陶組；4—古近系；5—侏羅系-白堊系；6—石炭系-二疊系；7—奧陶系；8—斷裂圖2 地?zé)岬刭|(zhì)剖面圖

(2)東營(yíng)組

受區(qū)域構(gòu)造和基底起伏的控制，東營(yíng)組底板埋深在2000m左右，厚度小于500m，在凹陷盆地的中心厚度最大，在盆地邊緣最薄，分布不穩(wěn)定，呈現(xiàn)自西向東、自南向北變厚趨勢(shì)?；讛嗔褬?gòu)造較發(fā)育，構(gòu)成良好的熱源上涌通道。熱儲(chǔ)巖性為細(xì)砂巖、砂礫巖，砂巖熱儲(chǔ)累計(jì)厚度100m左右，單井涌水量720～1440m3/d，地?zé)崴V化度6～10g/L，水化學(xué)類(lèi)型為Cl-Na，Cl·SO4-Na型，井口出水溫度50～60℃，屬溫?zé)崴?-熱水型低溫地?zé)豳Y源。

2.2 地溫梯度特征

為了解研究區(qū)地?zé)崴乃瘜W(xué)特征，本次對(duì)施工地?zé)峋乃|(zhì)分別進(jìn)行了常規(guī)離子、微量元素、污染組分、放射性元素及同位素等方面的檢測(cè)。

(1)水平變化特征

研究區(qū)內(nèi)地溫梯度總體大于3.4℃/100m(圖3)。在水平方向上總的規(guī)律是中間低兩邊高，莘縣東南部分地溫梯度最高。地溫場(chǎng)的展布方向呈帶狀NNE向，地溫在平面上的變化規(guī)律與區(qū)內(nèi)基巖埋藏深淺、凹凸基底的構(gòu)造輪廓、斷裂活動(dòng)的分布具一致性。

1—地溫梯度小于3.6℃/100m；2—地溫梯度3.6～3.8℃/100m；3—地溫梯度大于3.8℃/100m；4—地溫梯度分區(qū)界線；5—施工地?zé)峋畧D3 蓋層地溫梯度分區(qū)圖

(2)垂向變化特征

地溫場(chǎng)的垂向變化規(guī)律主要由井內(nèi)測(cè)溫方法獲得。本次施工地?zé)峋删疃?900m，井底溫度58℃，取水段為1340～1830m，取水段以上熱儲(chǔ)蓋層地溫梯度較高，約3.60℃/100m，熱儲(chǔ)層地溫梯度相對(duì)降低，約3.42℃/100m，綜合計(jì)算地溫梯度為3.51℃/100m。測(cè)溫曲線總體上是一條斜直線(圖4)，反應(yīng)在垂直方向上，地溫隨井深度的增加而遞增，但曲線上也有波折，分析其曲線形態(tài)，主要與地層巖性變化和砂泥巖比例有關(guān)[2-8]，綜合計(jì)算地溫梯度為3.57℃/100m。

圖4 蓋層地溫梯度分區(qū)圖

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

(1)常規(guī)離子成分

根據(jù)測(cè)試結(jié)果，莘縣地區(qū)地?zé)崴V化度介于7163.03～7776.14mg/L之間(表1)，水化學(xué)類(lèi)型屬于Cl-Na或Cl·SO4-Na型。陰離子以氯離子為主，含量介于2957.88～3350.03mg/L；陽(yáng)離子以鈉離子為主，含量介于2105.25～2343mg/L?？傆捕冉橛?220.98～1246.0mg/L之間，pH介于7.46～7.47之間。

表1 地?zé)崴谐Ａ吭睾恳挥[表 mg/L

(2)微量元素

地?zé)崃黧w由于深層循環(huán)過(guò)程中的水—巖相互作用，同時(shí)受到高溫對(duì)溶濾作用的影響，地?zé)崴泻胸S富的微量元素組分。根據(jù)測(cè)試結(jié)果，研究區(qū)地?zé)崴衅鹚岷?2.04mg/L，偏硅酸含量35.75mg/L，鍶含量8.9mg/L，碘含量1.70mg/L，鋰含量0.32mg/L，溴含量5.50mg/L(表2)。

表2 地?zé)崴形⒘吭亟M分含量一覽表 mg/L

(3)污染組分含量及其特征

根據(jù)測(cè)試結(jié)果，地?zé)崴杏泻Τ煞趾臀廴疚锖繕O微，大部分低于檢出限(表3)。

表3 污染成分含量一覽表 mg/L

(4)放射性元素及其特征

地?zé)崴蟹派湫晕镔|(zhì)皆高于一般地下水，放射性元素含量差別較大。根據(jù)其測(cè)試結(jié)果，區(qū)內(nèi)熱水中總α含量<4.3×10-2βq/L，總β含量0.027βq/L，222Rn含量6558mβq/L。

(5)同位素及其特征

根據(jù)測(cè)試結(jié)果，區(qū)內(nèi)地?zé)崴笑腄×10-3含量為-72，δ18O×10-3含量為-9.1，其D，18O關(guān)系點(diǎn)均稍位于全球雨水線附近(δD=8δ18O+10)(圖5)，說(shuō)明區(qū)內(nèi)地?zé)崴鹪从诖髿饨邓?/p>

圖5 地?zé)崴腄和δ18O關(guān)系圖

3 地?zé)崃黧w質(zhì)量評(píng)價(jià)

莘縣地區(qū)地?zé)崴泻卸喾N對(duì)人體有益的微量元素，偏硼酸、偏硅酸、氟含量均達(dá)到了《地?zé)豳Y源地質(zhì)勘查規(guī)范》(GB/T11615-2010)中《理療熱礦水水質(zhì)標(biāo)準(zhǔn)》礦水濃度值，可命名為含硼、硅、氟型熱礦水(表4)。由于地?zé)崴械V化度、總硬度、氯化物、硫酸鹽等離子含量超標(biāo)，不適宜作為農(nóng)田灌溉用水和飲用水水源(表5、表6)。

表4 理療熱礦泉水水質(zhì)標(biāo)準(zhǔn)對(duì)比表 mg/L

表5 農(nóng)業(yè)灌溉用水水質(zhì)標(biāo)準(zhǔn)對(duì)比表 mg/L

表6 生活飲用水水質(zhì)標(biāo)準(zhǔn)對(duì)比表 mg/L

根據(jù)《地?zé)豳Y源地質(zhì)勘查規(guī)范》(GB/T 11615-2010)對(duì)研究區(qū)地?zé)崴M(jìn)行腐蝕性和結(jié)垢評(píng)價(jià)。研究區(qū)地?zé)崴疄楦g性鍋垢很多的地?zé)崴虼嗽诘責(zé)崴_(kāi)發(fā)利用過(guò)程中要考慮地?zé)崴畬?duì)金屬的腐蝕性[9-15]。

4 地?zé)豳Y源計(jì)算

4.1 地?zé)豳Y源量計(jì)算

本次采用“熱儲(chǔ)法”計(jì)算研究區(qū)館陶組和東營(yíng)組熱儲(chǔ)地?zé)豳Y源量。計(jì)算公式如下：

Q=Qr+Qw

Qr=AdρrCr(1-φ)(tr-t0)

QL=Q1+Q2

Q1=Aφd

Q2=ASH

Qw=QLCwρw(tr-t0)

式中：Q—熱儲(chǔ)中儲(chǔ)存的熱量(J)；Qr—巖石中儲(chǔ)存的熱量(J)；QL—熱儲(chǔ)中儲(chǔ)存的水量(m3)；Q1—截止到計(jì)算時(shí)刻，熱儲(chǔ)孔隙中熱水的靜儲(chǔ)量(m3)；Q2—水位降低到目前取水能力極限深度時(shí)熱儲(chǔ)所釋放的水量(m3)；Qw—水中儲(chǔ)存的熱量(J)；A—計(jì)算區(qū)面積(m2)；d—熱儲(chǔ)層厚度(m)；ρr,ρw—熱儲(chǔ)巖石密度、地?zé)崴芏?kg/m3)；Cr,Cw—熱儲(chǔ)巖石比熱、地?zé)崴谋葻?J/(kg·℃))；φ—熱儲(chǔ)巖石的空隙度，無(wú)量綱；tr—熱儲(chǔ)溫度(℃)；t0—當(dāng)?shù)啬昶骄鶜鉁?℃)；S—彈性釋水系數(shù)，無(wú)量綱；H—計(jì)算起始點(diǎn)以上高度(m)。

經(jīng)計(jì)算，本區(qū)館陶組地?zé)豳Y源量為224.46×1016J，折合標(biāo)準(zhǔn)煤7.66×107t；東營(yíng)組地?zé)豳Y源量為191.66×106J，折合標(biāo)準(zhǔn)煤6.55×107t。

4.2 地?zé)崴刹闪坑?jì)算

本次采用“開(kāi)采強(qiáng)度法”計(jì)算工作區(qū)館陶組和東營(yíng)組熱水可采量。計(jì)算公式如下：

Q可=t·ε·4bL

式中：ε—開(kāi)采強(qiáng)度(m3/d·m2)；t—開(kāi)采時(shí)間(d)；2L—開(kāi)采區(qū)長(zhǎng)度(m)；2b—開(kāi)采區(qū)寬度(m)。

經(jīng)計(jì)算，館陶組地?zé)崴刹闪繛?1.38×108m3，放熱量折合標(biāo)準(zhǔn)煤6.71×106t；東營(yíng)組地?zé)崴刹闪繛?.66×108m3，放熱量折合標(biāo)準(zhǔn)煤6.09×106t。

4.3 回灌條件下地?zé)峥刹闪坑?jì)算

根據(jù)《地?zé)豳Y源評(píng)價(jià)方法及估算規(guī)程》(DZ/T0331--2020)，考慮回灌條件，采用熱突破公式計(jì)算地?zé)峄毓鄺l件下流體可開(kāi)采量。計(jì)算公式如下：

ρeCe=φρwCw+(1-φ)ρrCr

式中：Qa—回灌條件下允許開(kāi)采量(m3/d)；t—熱突破時(shí)間(d)；R1—井間距(m)；A—評(píng)價(jià)面積(m2)；M—熱儲(chǔ)層厚度(m)；φ—熱儲(chǔ)巖石孔隙度；ρw，ρr—熱儲(chǔ)水的密度，巖石的密度(kg/m3)；Cw，Cr—熱儲(chǔ)水的比熱，巖石的比熱(kJ/(kg·℃))。熱突破時(shí)間t取100年，即36500d。

經(jīng)計(jì)算，館陶組可采資源量為18.40×106m3/a，東營(yíng)組可采資源量為13.89×106m3/a。

5 結(jié)論

(1)莘縣地區(qū)地?zé)豳Y源主要賦存于新生代新近紀(jì)和古近紀(jì)碎屑沉積巖中。

(2)莘縣地區(qū)館陶組地?zé)豳Y源量為224.46×1016J，地?zé)崴刹闪繛?1.38×108m3，回灌條件下地?zé)崴刹闪繛?8.40×106m3/a；東營(yíng)組地?zé)豳Y源量為191.66×1016J，地?zé)崴刹闪繛?.66×108m3，回灌條件下地?zé)崴刹闪繛?3.89×106m3/a。

(3)莘縣地區(qū)地?zé)崴贿m宜作為農(nóng)田灌溉和飲用水源，且具有腐蝕性。