東營市南展區(qū)砂巖熱儲地?zé)峄毓嗔颗c溫度的關(guān)系探討

賀淼，張樂，袁一鳴，張建英，高延華

(1.東營市國土資源局，山東 東營 257091；2.東營市國土資源局東營分局,山東 東營 257091；3.山東省地礦工程集團(tuán)有限公司，山東 濟(jì)南 250014)

水文地質(zhì)環(huán)境地質(zhì)

東營市南展區(qū)砂巖熱儲地?zé)峄毓嗔颗c溫度的關(guān)系探討

賀淼1，張樂1，袁一鳴1，張建英2，高延華3

(1.東營市國土資源局，山東 東營 257091；2.東營市國土資源局東營分局,山東 東營 257091；3.山東省地礦工程集團(tuán)有限公司，山東 濟(jì)南 250014)

東營地區(qū)地?zé)豳Y源豐富，主要開采層位為新近紀(jì)館陶組、古近紀(jì)東營組，熱儲類型為低溫溫?zé)崴?、熱水型層狀孔?裂隙型砂巖熱儲，地下熱水補(bǔ)給條件較差。近年來，隨著地?zé)豳Y源的開發(fā)利用規(guī)模的加大，開采密集地帶的熱水頭明顯下降，存在資源枯竭，引發(fā)地質(zhì)環(huán)境問題的威脅。采用回灌方式補(bǔ)充地下水源，是實現(xiàn)地?zé)豳Y源可持續(xù)利用和地質(zhì)環(huán)境保護(hù)的有效措施。為研究該區(qū)地?zé)峄毓嗫尚行耘c回灌制約因素，在南展區(qū)開展了多次回灌試驗，試驗表明，該層熱儲回灌空間大，回灌條件較好，壓力對回灌效果影響明顯，回灌水體溫度對回灌效果影響較大。利用回灌量與回灌水頭之間的影響關(guān)系，模擬計算不同溫度下回灌水頭的變化情況，探討了回灌量與回灌水體溫度的影響關(guān)系，對回灌研究具有一定的指導(dǎo)意義。

地?zé)豳Y源；砂巖熱儲；溫度；回灌量；南展區(qū)；東營市

地?zé)峄毓嗍侨斯⒗眠^的地?zé)嵛菜虺氐叵滤?、地表水灌回到熱儲中[1-3]，是維持地?zé)豳Y源可持續(xù)開采和預(yù)防地面沉降等地質(zhì)環(huán)境問題[4-8]的有效措施，但是地?zé)峄毓嗉夹g(shù)復(fù)雜,回灌過程中可能引起回灌井井管及其周圍熱儲的物理、化學(xué)阻塞,使得回灌效率降低,并有可能導(dǎo)致回灌井損壞等。為此，在東營市開展了多次回灌試驗研究①，對回灌影響因素[9-10]進(jìn)行了分析與研究，提出了回灌水頭(壓力)、回灌井工藝、回灌水源等對回灌效果的影響。根據(jù)以往研究成果，回灌壓力影響回灌效果，溫度直接影響回灌水密度，進(jìn)而影響回灌壓力，因此，回灌水溫度對回灌效果的影響同樣不容忽視。南展區(qū)回灌試驗采用同層對井回灌模式，試驗持續(xù)時間較長，取得了豐富的動態(tài)數(shù)據(jù)，通過研究該區(qū)回灌試驗[11-14]中回灌水頭高度、回灌溫度、回灌量三者之間的關(guān)系，模擬計算不同溫度下的水頭高度，進(jìn)一步分析回灌溫度與回灌量的關(guān)系。

1 回灌試驗概況

1.1 熱儲條件

南展區(qū)位于東營市六戶鎮(zhèn)東南，行政區(qū)劃隸屬東營市東營區(qū)。構(gòu)造單元劃分上，該區(qū)位于華北坳陷濟(jì)陽坳陷區(qū)東營凹陷東南部的牛莊洼陷內(nèi)，區(qū)域上發(fā)育有王崗斷裂構(gòu)造背斜帶。區(qū)內(nèi)熱儲類型為層狀熱儲[15]，屬低溫地?zé)豳Y源中的溫?zé)崴獰崴偷責(zé)豳Y源[1],主要熱儲層位為古近紀(jì)東營組及新近紀(jì)館陶組。

(1)東營組熱儲底板埋深約1600～1900m，與下伏沙河街組呈整合或假整合接觸，砂層厚度130～200m，該層巖性為砂泥巖互層，半膠結(jié)狀，正韻律沉積，下部顆粒粗，以含礫砂為主，孔隙度較大，平均孔隙度約29%，富水性較強(qiáng)，單井涌水量大，可達(dá)80～100m3/h，水溫53～76℃，礦化度18.0～20.0g/L，水化學(xué)類型主要為Cl-Na·Ca型和Cl-Na型。

(2)館陶組熱儲頂板埋深1000m左右，底板埋深1350～1500m，與下伏古近紀(jì)東營組呈不整合接觸關(guān)系，砂層厚度在120～150m之間，蓋層分上下兩段，下段巖性為灰色、灰白色塊狀礫巖、含礫砂巖、砂巖。上段巖性為紫紅色、灰綠色泥巖與粉砂巖互層，孔隙度大，平均孔隙度為31%,熱儲溫度一般為53～55℃，單井涌水量大，一般在80m3/h左右，礦化度17.8～20.0g/L，水化學(xué)類型為Cl-Na型。

1.2 試驗場地

圖1 回灌試驗平面位置示意圖

地?zé)峄毓嘣囼灢捎昧嗽搮^(qū)域內(nèi)一對地?zé)峥碧骄脑囼灁?shù)據(jù)，其中一口為水源井(ZK01)，另一口為回灌井(ZK02)，水源井成井深度1716m，井壁管規(guī)格為177.8mm×8.05mm，濾水管長128.91m，取水段區(qū)間為1463.33～1704.91m，降深21.55m時，涌水量為90.8m3/h；回灌井成井深度1750m，井壁管規(guī)格為177.8mm×8.05mm，濾水管長128.64m，取水段區(qū)間為1553.6～1745.9m，降深22.20m時，涌水量為100m3/h。兩井取水層均為東營組，兩井之間間距532.3m?；毓嘣囼灢捎米冾l、回灌加壓泵控制回灌量大小，除砂器、粗(精)效過濾器對回灌水進(jìn)行二級過濾，以及排氣裝置、反沖儲水箱、反沖水泵等回灌設(shè)備，回灌井井口安裝可密封井口裝置，兩井之間通過鍍鋅管、法蘭進(jìn)行連接，回灌水由水源井經(jīng)輸水管道注入回灌井。試驗前期，采取小流量進(jìn)行回灌，回灌水溫度逐漸上升然后趨于穩(wěn)定，回灌井水位逐漸上漲，當(dāng)水頭突破初始水位與井口的高度差產(chǎn)生噴涌后，進(jìn)行密封回灌，回灌壓力、回灌溫度及回灌量同步進(jìn)行觀測，根據(jù)試驗情況，逐步提高回灌量，最終達(dá)到設(shè)計的穩(wěn)定回灌壓力，并維持一定的回灌時間,試驗共持續(xù)50.5d，總回灌量約5.98萬m3。

圖2 試驗工藝流程示意圖

2 回灌量與回灌溫度關(guān)系分析

2.1 回灌量與水頭高度關(guān)系分析

為研究回灌量同回灌壓力的關(guān)系[16],該試驗進(jìn)行了4個流量段回灌，試驗初期水頭高度通過人工測量得出，密封回灌階段水頭高度通過壓力表讀數(shù)進(jìn)行換算，回灌過程中的水頭高度以回灌前靜水位為起算點(diǎn)。根據(jù)回灌試驗數(shù)據(jù)，當(dāng)回灌量為42.6m3/h時，水頭高度46.45m，隨著回灌量增加，水頭高度不斷升高，當(dāng)回灌量增加至57.7m3/h時，水頭高度85.21m。結(jié)合回灌量與水頭高度關(guān)系圖(圖3)，可以看出回灌量與回灌壓力之間相互影響，二者呈正比例關(guān)系，即回灌量與水頭高度之間呈正比例關(guān)系，兩者線性相關(guān)關(guān)系式近似為：

H≈2.56Q-66.73

式中：H—水頭高度(m)；Q—回灌量(m3/h)。

圖3 回灌量-水頭高度關(guān)系圖

2.2 回灌溫度與水頭高度關(guān)系分析

回灌壓力與回灌量相互關(guān)聯(lián)，而使回灌壓力發(fā)生變化的因素除了回灌量以外，流體重力也是重要影響因素，根據(jù)物理學(xué)原理，流體重力與其體積、密度有關(guān)，根據(jù)地?zé)峋删Y(jié)構(gòu)，回灌井井筒[17-18]內(nèi)總體積不發(fā)生變化，密度是影響回灌水重力的唯一因素。流體溫度影響密度，進(jìn)而影響重力，因此，可借助回灌量與水頭高度的關(guān)系，研究溫度與水頭高度的關(guān)系，進(jìn)一步分析溫度與回灌量的關(guān)系。該次地?zé)崴芏葏⒄铡兜責(zé)豳Y源地質(zhì)勘查規(guī)范》(GB/T11615-2010)不同溫度下水的密度，進(jìn)行適應(yīng)性插分取值。

圖4 回灌井剖面圖

討論二者之間的關(guān)系，暫不考慮回灌目的層壓力的變化，只研究井管內(nèi)由回灌水溫度改變引起的壓力變化。假設(shè)回灌井井管內(nèi)原有地?zé)崴换毓嗨赐耆脫Q，置換后井管內(nèi)水溫一致，以南展區(qū)地?zé)峄毓嘣囼炛禐閰⒖?，模擬計算在保持回灌量不變，即保持井底壓力[19-20](濾水管底壓力)不變的條件下，以不同溫度(50℃，40℃，30℃)水源進(jìn)行回灌的回灌井水頭高度情況。以井口處為分界面，井底承受壓力為外界壓力與井筒內(nèi)熱水產(chǎn)生的壓力之和，井底壓力計算公式為：

F=F外+F內(nèi)=(P加+ρ0gh0)×S=ρ1g×(H+h0-h靜)

得出水頭高度計算公式：

式中：F外—密封回灌回灌壓力；F內(nèi)—井筒內(nèi)熱水產(chǎn)生壓力；S—井底面積；P加—密封回灌產(chǎn)生壓強(qiáng)(MPa)(數(shù)據(jù)來源于東營市南展區(qū)回灌試驗資料)；ρ0—實際回灌試驗中熱水密度(kg/m3)；ρ1—模擬回灌熱水密度(kg/m3)；g—重力值，取9.8N/kg；h0—回灌井濾水管埋深(m)；h靜—靜水位埋深(m)；H—回灌后水頭高度(m)。

計算結(jié)果如表1所示。

表1 回灌溫度與水頭高度對應(yīng)關(guān)系

計算結(jié)果顯示，在保持一定回灌量的情況下，隨著回灌水溫度降低，水頭高度也在下降，為確定二者關(guān)系，繪制溫度-水頭高度曲線(圖5)，該曲線近拋物線形式，建立如下公式:

H=at2+bt+c

式中：H—水頭高度(m)；t—回灌溫度(℃)。

根據(jù)表1計算結(jié)果，代入法求得a≈0.01，b≈-0.25,c為常數(shù)N。得出水頭高度與回灌溫度之間的關(guān)系式：

H=0.01t2-0.25t+N

圖5 溫度-水頭高度關(guān)系曲線圖

2.3 回灌量與回灌溫度關(guān)系分析

2.3.1 回灌量與回灌溫度關(guān)系分析

在分析回灌溫度同水頭高度的公式關(guān)系時，發(fā)現(xiàn)常數(shù)N只受回灌量影響，且與回灌量呈正相關(guān)，當(dāng)回灌量一定時，N為定值，不妨將其當(dāng)作回灌量，結(jié)合水頭高度與回灌量呈正相關(guān)的關(guān)系，分析回灌溫度與回灌量的關(guān)系。在保持回灌量不變的情況下，隨著溫度升高，水頭高度增大。反之，回灌溫度變小，回灌井水頭高度變小，井內(nèi)剩余回灌體積空間增加，要達(dá)到原來的水頭高度，回灌量會相應(yīng)增加。即：兩者之間呈反比例關(guān)系。

2.3.2 回灌量與回灌溫度干擾因素分析

根據(jù)以往地?zé)峄毓嘟?jīng)驗，影響回灌效果的因素較多，包括回灌井成井結(jié)構(gòu)、熱儲層特征、回灌壓力大小、水化學(xué)特征等。該次研究回灌量與回灌溫度之間關(guān)系過程中忽略了熱儲層特征變化及兩種地?zé)崴旌虾笏瘜W(xué)特征的影響，因此，確定二者為該次研究的干擾因素。

(1)熱儲特征因素分析

熱儲特征通過熱儲特征參數(shù)體現(xiàn)，主要包括回灌滲透系數(shù)、導(dǎo)水系數(shù)、導(dǎo)壓系數(shù)(表2)。隨著回灌試驗進(jìn)行，熱儲特征參數(shù)發(fā)生改變，熱儲層滲透性能發(fā)生改變。熱儲層結(jié)構(gòu)復(fù)雜，不同溫度水體進(jìn)入后，對熱儲層的影響及其反作用目前難以明確，只研究井管內(nèi)水頭高度變化，可降低熱儲特征因素的干擾。

表2 回灌熱儲特征參數(shù)

(2)水化學(xué)因素分析

圖6 CaCO3溶解度曲線

3 結(jié)論

以水頭高度為結(jié)合點(diǎn)，利用回灌量、回灌溫度同水頭高度關(guān)系，創(chuàng)建回灌量與回灌溫度的關(guān)系式，經(jīng)分析得出：

(1)回灌溫度通過改變回灌壓力對回灌量產(chǎn)生影響，二者呈反比例關(guān)系，即低溫水源回灌效果理論上要優(yōu)于高溫水源。

(2)回灌溫度與回灌量的關(guān)系式是在南展區(qū)回灌試驗基礎(chǔ)上模擬推導(dǎo)得出，對于其他干擾因素的影響關(guān)系尚不明確，有待進(jìn)一步優(yōu)化。根據(jù)單井試驗建立溫度、回灌量之間的關(guān)系，對今后的回灌研究具有一定的指導(dǎo)意義。

[1] 周世海,楊詢昌,梁偉,等.德州市城區(qū)地?zé)崴斯せ毓嘣囼炑芯縖J].山東國土資源,2007,23(9):11-14.

[2] 張新文,胡彩萍,胡松濤,等.東營市中心城區(qū)地?zé)崽锏責(zé)峄毓嗫尚行苑治鯷J].山東國土資源,2009,25(8):39-42.

[3] 趙季初.魯北砂巖熱儲地?zé)嵛菜毓嘣囼炑芯縖J].山東國土資源，2013,29(9):24-31.

[4] 李蓮花,張建斌.地?zé)崴Y源開發(fā)引起的環(huán)境問題分析[J].地下水，2004,26(1):194-195．

[5] 李曉華,馬馳.地?zé)崴Y源開發(fā)引起的環(huán)境問題分析[J].安徽農(nóng)業(yè)科學(xué)，2009,37(18):8675-8677.

[6] 胡彩萍,張景燕,郝夢圓,等.黃河三角洲生態(tài)區(qū)孔隙熱儲地?zé)衢_發(fā)對地面沉降的影響分析[J].山東國土資源,2017,33(2):39-42.

[7] 胡彩萍,潘擁軍.東營市地?zé)豳Y源開發(fā)保護(hù)與行政保障措施探討[J].山東國土資源,2012,28(9):81-83.

[8] 王彥俊,王貞國,王巖,等.德州市地?zé)豳Y源開發(fā)與保護(hù)[J].山東國土資源,2005,21(5):31-33．

[9] 高寶珠,曾梅香.地?zé)釋\(yùn)行系統(tǒng)中回灌井堵塞原因淺析及預(yù)防措施[J].水文地質(zhì)工程地質(zhì),2007,34(2):75-80.

[10] 林黎,王連成,趙蘇民,等.天津地區(qū)孔隙型熱儲層地?zé)崃黧w回灌影響因素探討[J].水文地質(zhì)工程地質(zhì),2008,35(6):125-128.

[11] 王坤,朱家玲.中低溫孔隙型地?zé)崽锘毓嘣囼炑芯縖J].太陽能學(xué)報,2001,22(2):236-238．

[12] 高宗軍,郭加朋,李哲,等.東營市城區(qū)地?zé)醿θ斯せ毓鄺l件及分區(qū)研究[J].地下水,2009,31(5):4-8.

[13] 譚志容.東營市城區(qū)館陶組熱儲回灌性能分析[J].山東國土資源，2010,26(8):13-17.

[14] 張平平,王秀芹.回灌水溫對砂巖熱儲回灌效果的影響分析——以德州市地?zé)崽镳^陶組熱儲為例[J].山東國土資源,2015,31(6),68-71.

[15] 徐希強(qiáng),劉善軍,王偉德,等.山東省地?zé)豳Y源及勘查開發(fā)對策[J].山東國土資源,2015,31(5):31-35.

[16] 楊詢昌，王成明，馮守濤．德州市城區(qū)地?zé)崴畡討B(tài)與開采量關(guān)系[J].山東國土資源，2011,27(7):25-28．

[17] 陳國富,陳虹雁,王學(xué)工.井筒內(nèi)水頭高度與井底壓力對比試驗[J].工程勘察,1997,(4):32-35．

[18] 劉志濤,張志勝,胡松濤,等.地?zé)峋矡嵝?yīng)研究[J].山東國土資源,2014,30(8):54-57.

[19] 周萌,任曉霞,楊淑云.地?zé)峄毓嗑讐簭?qiáng)新解[J].山東國土資源,2016,32(1):38-41.

[20] 張春志,李朋,張秋鋒.魯東構(gòu)造型地?zé)崽锘毓噙^程中壓力對回灌量的影響研究——以威海市寶泉湯為例[J].山東國土資源,2015,31(2):17-20.

StudyontheRelationshipbetweenRe-injectionVolumeandTemperatureofSandstoneGeothermalReservoirinNanzhanAreaofDongyingCity

HE Miao1, ZHANG Le1, YUAN Yiming1, ZHANG Jianying2, GAO Yanhua3

(1.Dongying Bureau of land and resources, Shandong Dongying 257091, China;2.Dongying Branch Bureau of Dongying Bureau of Land and Resources , Shandong Dongying 257091, China;3. Shandong Geo-engineering Exploration Limited Corporation, Shandong Jinan 250014, China)

There is rich geothermal resources in Dongying area, and the main mining layers are Neogene Guantao formation and Paleogene Dongying formation. Thermal reservoir are warm -hot water type and layered pore fissure type sandstone reservoir. The geothermal water supply condition is poor. In recent years, accompanying with the increase of the scale of exploitation and utilization of geothermal resources, geothermal water head in the mining area has been decreased obviously, and there is a threat to resource depletion and geological environment. Re-injection is an effective way to realize the sustainable utilization of geothermal resources and the protection of geological environment. In order to study the feasibility of geothermal reinjection and controlling factors for re-injection in Dongying city, recharge tests have been carried out in Nanzhan area. It is showed that the reservoir has a large space and the recharge condition is good, pressure on the recharge effect is obvious, but the influence of temperature on the water recharge recharge effect can not be ignored. Based on the test process, by using the relation between recharge and recharge of water head, the changes under different temperatures have been simulated and calculated, the influence of the injection volume and injection water temperature has been studide. It will provide some references for re-charge study in future.

Geothermal resources; sandstone geothermal reservoir; temperature; re-injection volume; Nanzhan area; Dongying city

2017-04-24；

2017-09-04；

曹麗麗

賀淼(1984一)，男，山東東營人，工程師，主要從事礦產(chǎn)勘查工作；E-mail:dydk05@126.com

山東省地礦工程集團(tuán)有限公司，山東省東營市南展區(qū)地?zé)豳Y源普查報告，2014年12月。

P314.1

B

賀淼，張樂，袁一鳴,等.東營市南展區(qū)砂巖熱儲地?zé)峄毓嗔颗c溫度的關(guān)系探討[J].山東國土資源，2018,34(1)：44-48.HE Miao, ZHANG Le, YUAN Yiming, etc. Study on the Relationship between Re-injection Volume and Temperature of Sandstone Geothermal Reservoir in Nanzhan Area of Dongying City[J].Shandong Land and Resources, 2018,34(1)：44-48.