王先超

山東泰山資源勘查有限公司 山東 濟南 250104

地?zé)豳Y源作為一種可再生能源，與化石能源相比，具有強大的社會經(jīng)濟與環(huán)境優(yōu)勢，與風(fēng)能和太陽能等相比，其可利用系數(shù)高、連續(xù)穩(wěn)定、開發(fā)利用方便，是最具競爭力的資源之一[1-3]；地?zé)豳Y源還具有儲量大、分布廣和開發(fā)利用安全、高效的特點[4]。合理可持續(xù)開發(fā)地?zé)豳Y源有助于我國雙碳目標(biāo)的實現(xiàn)。雄安新區(qū)建設(shè)綠色智慧生態(tài)城市歸功于豐富的地?zé)豳Y源以及良好的開采條件[5]。雄安新區(qū)尤其是牛駝鎮(zhèn)凸起具有多年的地?zé)衢_發(fā)歷史[6-8]，雄縣地?zé)岬睦靡殉梢?guī)?；?，是我國第一個用地?zé)崽娲济汗┡摹盁o煙城”[9]，形成了政企合作、統(tǒng)一管理、集中回灌的“雄縣模式”[7]。

雄安新區(qū)范圍內(nèi)高陽凸起區(qū)目前地?zé)豳Y源開發(fā)利用的主要是館陶組砂巖熱儲水熱型地?zé)豳Y源[10]，回灌對于水熱型地?zé)豳Y源的可持續(xù)開發(fā)有著非常重要的作用。首先它可以保持維持熱儲壓力不下降，其次可以防止熱水對地表及淺層地下水體產(chǎn)生污染[11]。但是，地?zé)峄毓嗉夹g(shù)復(fù)雜難操作，需要考慮多方面的因素影響。地?zé)崽锏膲勖才c回灌過程息息相關(guān)。為此，需要采用數(shù)值模擬技術(shù)，模擬熱儲的溫度與壓力在不同采灌情景下的響應(yīng)。利用數(shù)值模擬研究地下熱儲中熱水系統(tǒng)的運動規(guī)律，已成為地?zé)崃黧w開發(fā)和評價的重要技術(shù)方法[12-16]。雄安地區(qū)地?zé)豳Y源的規(guī)模開發(fā)需要制定科學(xué)合理的采灌方案，是保障地?zé)豳Y源可持續(xù)利用的有效途徑。本文建立雄安新區(qū)高陽低凸起區(qū)館陶組熱儲水-熱數(shù)值模型，根據(jù)現(xiàn)狀開釆現(xiàn)狀，設(shè)計兩個對井方案，一個群井方案，分析不同采灌方案對熱儲層溫度場和壓力場的影響程度，為雄安地區(qū)地?zé)豳Y源合理開發(fā)、可持續(xù)開發(fā)提供依據(jù)。

1 地?zé)岬刭|(zhì)特征

1.1 地?zé)岬刭|(zhì)概況

雄安新區(qū)安新縣中南部地區(qū)地質(zhì)構(gòu)造上為高陽臺凸的一部分，也屬于高陽地?zé)崽铩０残驴h有華北最大的淡水湖泊白洋淀，地表水、地下水資源均豐富[10]（劉偉坡等，2013），為開發(fā)地?zé)崽峁┝肆己玫臈l件。安新主體位于高陽低凸起和-保定凹陷，北部一部分位于牛駝鎮(zhèn)凸起、容城凸起和徐水凹陷[8,17]（圖1）。

1.2 地?zé)醿?/h3>

本區(qū)熱儲類型分為兩種類型：孔隙型熱儲（第三系熱儲）和深循環(huán)構(gòu)造裂隙型熱儲（基巖熱儲），前者是以熱傳導(dǎo)為主的，在大地?zé)崃骺刂谱饔孟滦纬傻模笳呤窃跓醾鲗?dǎo)與熱對流共同作用下形成的。目前4000m以淺新近系明化鎮(zhèn)組、館陶組孔隙熱儲與基巖熱儲為可開采熱儲層[18]。

1.3 成藏模式

該地區(qū)沉積盆地傳導(dǎo)型地?zé)豳Y源成藏模式如圖2所示，淺部地層的熱來源于上地幔生熱與基巖的放射性衰變，主要是通過傳導(dǎo)的方式，上部的沉積蓋層可以起到保溫的作用，平原區(qū)的熱對流導(dǎo)水導(dǎo)熱通道傳輸，起到額外補充熱量的作用，熱儲層形成于高孔隙度和高滲透性含水層處[5]。

2 雄安新區(qū)高陽低凸區(qū)熱儲數(shù)值模型描述

對于大區(qū)域熱儲層數(shù)值模擬，為提高計算速度，對研究區(qū)內(nèi)熱儲層及地?zé)崃黧w做如下概化：（1）非穩(wěn)定流、不可壓、三維流動；（2）均質(zhì)各向同性；（3）流體自然對流滿足達(dá)西定律；（4）水與飽和多孔介質(zhì)層處于熱力學(xué)平衡狀態(tài)[19-20]。

2.1 3D模型建立

整個高陽地?zé)崽镌u價面積為727.15 km2，建模時，邊界條件根據(jù)該區(qū)構(gòu)造特征和水文地質(zhì)條件進(jìn)行確定，補給、排泄邊界均作已知流量邊界處理，流量隨時間變化，變化情況參考現(xiàn)有動態(tài)資料。通過不斷調(diào)整邊界流量，直到區(qū)域水位規(guī)律基本符號實際情況，也就確定了相應(yīng)的邊界流量。隔水邊界為模型的頂部和底部水流邊界，定熱流量邊界為模型頂部熱量邊界。水平方向上，研究區(qū)模型范圍以高陽低凸起區(qū)為邊界；西部為保定斷裂；北部為牛南斷裂；東部和南部為雄安新區(qū)邊界。根據(jù)多年氣象及有關(guān)地?zé)岬刭|(zhì)資料確定初始條件，確定恒溫帶深度為25m，恒溫帶溫度為14.5℃，安新縣霧迷山組地溫梯度為3.7℃/100m。根據(jù)研究區(qū)地溫梯度及恒溫層的溫度以及水位統(tǒng)測結(jié)果，計算得到初始溫度和壓力并且一并賦值給模型相應(yīng)的各個單元格。模型初始壓力場見圖3。

2.2 模型結(jié)果驗證

2.2.1 溫度擬合

本研究采用試估-校正法對模型參數(shù)進(jìn)行識別和驗證。研究區(qū)溫度垂直分布基本上是隨深度增加而增加,館陶組溫度分布50℃～70℃左右（圖4），對比圖5研究區(qū)主要地?zé)峋疁胤植紙D，模型分析結(jié)果與地?zé)峋壳伴_采溫度基本一致。

2.2.2 水位擬合

研究區(qū)內(nèi)有一口監(jiān)測井，位于研究區(qū)中北部的D32井，D32井在2019年1月15日至1月21日期間進(jìn)行過一次抽水試驗，歷時158 h。結(jié)果顯示，計算值與觀測值很接近，模擬期、驗證期內(nèi)模擬數(shù)位與觀測水位曲線基本擬合（圖6）。所建數(shù)值模型能夠代表真實的地?zé)嵯到y(tǒng)，可以對未來熱儲層壓力變化進(jìn)行預(yù)報。

3 地下熱水回灌模擬研究

3.1 回灌方案

采灌模式根據(jù)回灌井分布方式主要分為兩種，一種為分散采灌，第二種為外圍集中采灌，本文分別對兩種模式進(jìn)行模擬，分析對比評估兩種模式的效果以及熱儲的響應(yīng)。

3.1.1 對井方案

對井采灌是當(dāng)前水熱型開采最常用的布局模式。根據(jù)文獻(xiàn)對井采灌最優(yōu)布局的研究[17]，采灌井距有一個邊界距離，當(dāng)采灌井距大于400m的時候，溫度變化對井距變化不敏感，當(dāng)大于這個距離的時候，注入的冷水在給定的時間內(nèi)，對開采井溫度影響不大。因此將對井井間距設(shè)置為400m，回灌井位于開采井下游(圖7)。

3.1.2 群井方案

對丼方案由于布置的井?dāng)?shù)最偏少，導(dǎo)致可開采熱量受到限制，現(xiàn)提出圖10所示的群井布置方案。根據(jù)雄安新區(qū)地下水穩(wěn)定同位素分析，地下水的運動方向是由西部和北部向東南運移，為了減少集中回灌給開采井帶來的熱突破風(fēng)險，擬計劃在安置區(qū)東北方向集中布置回灌井，在西南方向集中布置開采井（圖8）。開釆井維持現(xiàn)狀開采，將各井的回灌流量與開采井相同，回灌水的溫度為45°C，回灌井采暖期回灌。

3.2 回灌模擬結(jié)果分析

3.2.1 對井方案Ⅰ

研究區(qū)開采模式采用“一采一灌”進(jìn)行，利用尾水進(jìn)行回灌，保持現(xiàn)有地?zé)峋_采格局和強度，開采井全年開采，在非采暖期用于生活熱水的采量很少，僅為采暖期總采量的5-10% ；回灌井在每年11月中旬至第二年3月中旬的供暖期進(jìn)行回灌，在其余時間即停暖期回灌井停灌。

在模型邊界條件和參數(shù)等配置好的基礎(chǔ)上，根據(jù)方案Ⅰ建立預(yù)測模型，運行得到研究區(qū)100年后預(yù)測的壓力場和溫度場如圖9，10所示，熱儲壓力漏斗中心在大規(guī)模集中開采的影響下不斷擴大，同時，由于在西部和北部開采規(guī)模較大，漏斗中心向西北發(fā)展。而在研究區(qū)范圍內(nèi)及其附近地區(qū)，由于回灌井對熱儲層壓力的支撐作用，水位相對增加，下降速率也相對緩慢。

地?zé)峋患皽囟入S時間變化曲線如圖11、12所示。在方案Ⅰ采灌條件下，根據(jù)預(yù)測模型得到水流場與驗證模型中現(xiàn)狀開采水流場差別不大，井底壓力下降0.18MPa, 相當(dāng)于水位下降18m，全區(qū)水位年降幅在2-3m左右。根據(jù)預(yù)測模型得到溫度場與驗證模型中現(xiàn)狀開采溫流場對比可知，在方案Ⅰ采灌條件下開采前60年，溫度場變化不大，60年之后溫度降低明顯。監(jiān)測井井底溫度下降12℃，全區(qū)溫度年降幅在0.15℃左右。

3.2.2 對井方案Ⅱ

方案Ⅱ中，為了進(jìn)一步分析熱儲層壓力變化趨勢受回灌的影響，開釆井維持現(xiàn)狀開采，調(diào)整回灌量大小，回灌量與開采量比例相同，回灌水的溫度為45°C，回灌井采暖期回灌。

根據(jù)方案Ⅱ建立預(yù)測模型，運行得到研究區(qū)地?zé)峋蛔兓€和溫度變化曲線如圖11、12所示。依據(jù)預(yù)測模型得到的水流場與現(xiàn)狀開采水流場比較可知，在方案Ⅱ采灌條件下壓力下降趨勢有了進(jìn)一步的改善，井底壓力下降0.14MPa，相當(dāng)于水位下降13.66m，全區(qū)水位年降幅在1.5m左右。依據(jù)預(yù)測模型得到溫度場與現(xiàn)狀開采溫流場對比可知，由于回灌溫度較開采方案Ⅰ低，溫度下降明顯，井底溫度下降14℃，全區(qū)水位年降幅在0.15℃左右。

3.2.3 群井方案

根據(jù)方案Ⅲ建立預(yù)測模型，運行得到研究區(qū)100年后地?zé)峋蛔兓€和溫度變化曲線如圖11、12所示。在集中采灌模式下，熱儲壓力在城區(qū)大部基本不變。在水流場作用下，回灌區(qū)注入的回灌水向城區(qū)集中開采區(qū)運動，有效地維持了熱儲壓力。從回灌區(qū)到開采區(qū)，壓力由上升逐漸過度到下降，在其中間地區(qū)形成了壓力穩(wěn)定帶。壓力下降的區(qū)域位于開采區(qū)；而東北部熱儲壓力受到回灌的影響有所上升。從圖11可以看到，與對井布并方案相比，開采溫度下降相對緩慢，但是，由于采灌井間距太遠(yuǎn)，集中回灌對穩(wěn)定和恢復(fù)熱儲壓力的效果明顯減弱了。

4 討論：對井方案與群井方案局限性分析

方案預(yù)測結(jié)果對比表1為三種方案熱儲層預(yù)測壓力和預(yù)測溫度，明顯地看出對井方案Ⅰ與群井方案的灌采比相當(dāng)，50%左右，對井方案Ⅱ的回灌率將近94%，三種方案的水位變化曲線和溫度變化曲線基本相同，熱流體壓力略有下降，但降幅不大。除了群井方案的流體壓力降大于前兩個開采方案，對井方案Ⅰ和對井方案Ⅱ壓力下降幅度基本相同；群井方案與現(xiàn)狀開采壓力下降幅度一致（圖11），這跟與集中回灌區(qū)與開采井的距離有直接的關(guān)系。

表1 三種方案熱儲層預(yù)測壓力和預(yù)測溫度

溫度場的動態(tài)變化：由于該層三個開采方案設(shè)計回灌量都較大，動態(tài)預(yù)測期末三個開采方案溫度都有變化，下降幅度基本以群井方案→對井方案Ⅰ→對井方案Ⅱ遞增。相對于群井方案，對井方案可以有效減少熱儲壓力的降低，但是熱儲溫度降低較大，很容易過早出現(xiàn)熱突破現(xiàn)象，群井方案有利于保持熱儲溫度，但是對于維持熱儲壓力效果不明顯。

5 結(jié)論與建議

（1）建立了高陽地?zé)崽镄郯残聟^(qū)范圍內(nèi)熱儲地?zé)岬刭|(zhì)模型，該區(qū)地?zé)豳Y源是中低溫地?zé)豳Y源，大氣降水深循環(huán)傳熱形成，聚熱機理以傳導(dǎo)為主，局部有熱對流發(fā)生，由蓋層、熱儲層、熱水通道和熱源構(gòu)成，形成了一個完整的地?zé)嵯到y(tǒng)。

（2）建立的水-熱數(shù)值模擬模型，對于預(yù)測地?zé)崽飳﹂_采量的反應(yīng)的比較準(zhǔn)確的。根據(jù)現(xiàn)狀開釆現(xiàn)狀，設(shè)計兩個對井方案，一個群井方案。預(yù)測結(jié)果顯示，三種開采方案的熱流體壓力動態(tài)變化基本相同，略有下降，但降幅不大。動態(tài)預(yù)測期末三個開采方案溫度由于較大的回灌量都有變化，下降幅度基本以群井方案→對井方案Ⅰ→對井方案Ⅱ遞增。相對于群井方案，對井方案可以有效減少熱儲壓力的降低，但是熱儲溫度降低較大，很容易過早出現(xiàn)熱突破現(xiàn)象，群井方案有利于保持熱儲溫度，但是對于維持熱儲壓力效果不明顯。

（3）高陽地?zé)崽镄郯残聟^(qū)范圍內(nèi)地?zé)豳Y源開釆以水熱型地?zé)衢_發(fā)為主。在采灌均衡的條件下，結(jié)合區(qū)域規(guī)劃建設(shè)部署合理布設(shè)釆灌模式，在分散式釆灌的開采布局下，應(yīng)重點考慮相鄰井的熱突破問題；在集中釆灌布局下，應(yīng)重點考慮集中開采區(qū)水位下降問題。

（4）建議加強地?zé)衢_發(fā)利用過程中對水量、水質(zhì)、水溫的長期動態(tài)監(jiān)測，形成年度評價機制和預(yù)測預(yù)警機制，每年對區(qū)域地?zé)豳Y源條件以及可供暖面積進(jìn)行評價，結(jié)合能源站布設(shè)和需求情況，形成長期動態(tài)調(diào)整機制，保障熱儲平衡。