層狀熱儲地?zé)峋畽?quán)益保護半徑計算探討

劉帥， 馮守濤， 劉志濤， 黃松， 黃星， 朱利民

(1.山東省地勘局第二水文地質(zhì)工程地質(zhì)大隊(山東省魯北地質(zhì)工程勘察院)，德州 253015； 2.山東省地?zé)崆鍧嵞茉刺綔y開發(fā)與回灌工程技術(shù)研究中心，德州 253015)

0 引言

積極推進冬季清潔能源取暖是改善空氣質(zhì)量的重要措施之一，地?zé)豳Y源作為集熱能和水資源為一體的可再生新能源和清潔能源，具有清潔環(huán)保、儲量大、分布廣、易開發(fā)、利用成本低、穩(wěn)定可靠等優(yōu)點，其開發(fā)利用受到越來越多學(xué)者們的關(guān)注[1-4]。20世紀(jì)90年代以來，地?zé)豳Y源開發(fā)利用發(fā)展迅猛，形成了以天津、陜西、河北、山東為代表的地?zé)峁┡瘏^(qū)和以北京、東南沿海為代表的溫泉旅游與療養(yǎng)等中深層地?zé)豳Y源直接利用區(qū)。至2015年底，全國水熱型地?zé)崮芄┡娣e達1.02億 m2[5]，實現(xiàn)年替代煤290萬 t，減排二氧化碳750萬 t，對減輕北方地區(qū)冬季霧霾發(fā)揮了積極作用。

地?zé)豳Y源大規(guī)模集中開采將導(dǎo)致地?zé)崴畨毫Σ粩嘞陆担纬刹煌?guī)模的地?zé)崴德渎┒?，從而影響地?zé)豳Y源的有效開發(fā)與利用，地?zé)嵛菜毓嗍墙鉀Q以上問題的有效途徑[6-7]。但由于采、灌井布局不合理，在回灌過程中可能會發(fā)生熱突破現(xiàn)象，造成開采井水溫降低，達不到供暖效果。因此，研究地?zé)衢_采井和回灌井的合理化布局非常重要。

地?zé)峋畽?quán)益保護半徑是指以地?zé)峋?開采井或回灌井)為中心的圓形區(qū)域半徑，在這個圓形區(qū)域內(nèi)，在整個開采(回灌)期內(nèi)，地?zé)峋纯砷_采量正常開采(回灌)，地?zé)岙a(chǎn)出量與流體溫度都不會受開采(回灌)的影響。地?zé)峋畽?quán)益保護半徑分為開采井權(quán)益保護半徑和回灌井權(quán)益保護半徑。地?zé)峋畽?quán)益保護半徑為確定開采井和回灌井合理化布局提供了量化依據(jù)。

《GB/T 11615—2010地?zé)豳Y源地質(zhì)勘查規(guī)范》[8]明確提出了單井開采權(quán)益保護半徑的概念和計算方法，但在實際工作中發(fā)現(xiàn)，f(水比熱與熱儲巖石比熱的比值)的定義及給定參考值不太合理，從而造成開采井權(quán)益保護半徑計算結(jié)果往往超出了批復(fù)的礦權(quán)范圍，導(dǎo)致地方政府、地?zé)衢_發(fā)主管部門及社會公眾產(chǎn)生疑慮： 是開采井權(quán)益保護半徑的計算結(jié)果偏大，還是批復(fù)的礦權(quán)范圍太小？計算結(jié)果不能對地?zé)豳Y源的合理開發(fā)、利用起到技術(shù)指導(dǎo)作用。此外，隨著地?zé)嵛菜毓喙ぷ鞯拇笠?guī)模開展，一些學(xué)者也提出了不同的回灌井權(quán)益保護半徑的計算公式[9-11]，但并未列出適用條件，計算結(jié)果以誰為準(zhǔn)也有待商榷。為保證地?zé)豳Y源的可持續(xù)開發(fā)及利用，本文參照以往研究成果，給出了明確的開采井和回灌井權(quán)益保護半徑的計算公式。

1 開采井權(quán)益保護半徑

1.1 現(xiàn)有計算方法

《GB/T 11615—2010地?zé)豳Y源地質(zhì)勘查規(guī)范》[8]給出了單井開采權(quán)益保護半徑的定義及計算方法,計算公式為

(1)

式中：R采為地?zé)峋_采100 a排出熱量對熱儲的影響半徑，m；Q為地?zé)崴砷_采量，m3/d；f為水比熱與熱儲巖石比熱的比值，為3～5；H為熱儲層厚度，m。

1.2 原理分析

《GB/T 11615—2010地?zé)豳Y源地質(zhì)勘查規(guī)范》[8]給出的開采井權(quán)益保護半徑公式是根據(jù)熱量均衡原理推導(dǎo)得出的，即地?zé)峋_采熱量=熱儲層損失的熱量。若已知物體比熱，可根據(jù)公式(2)推導(dǎo)出物體的放熱量，計算方法為

Q熱=cmΔt，

(2)

式中：Q熱為物體放出(吸收)的熱量，MJ；c為物體比熱，MJ/(kg·℃)；m為物體質(zhì)量，kg； Δt為溫度差，℃；

因此，地?zé)峋_采熱量可表示為cω(ρωQt)(T1-T0)，一定范圍內(nèi)熱儲層損失的熱量可表示為cα(λραπR開2H)(T1-T0)，則存在以下數(shù)學(xué)關(guān)系式

(3)

式中：Q采為地?zé)崴砷_采量，m3/d；λ為熱儲回收率，無量綱，松散巖類孔隙熱儲回收率取值0.25，巖溶裂隙熱儲回收率取值0.15[12]；H為地?zé)峋脽醿穸?，m；t為地?zé)峋_采總天數(shù)，d；R采為開采井權(quán)益保護半徑，m；ρωcω為流體密度與比熱的乘積，MJ/(m3·℃)；ραcα為熱儲層平均密度與平均比熱的乘積，MJ/(m3·℃);ρrcr為巖石密度與比熱的乘積，MJ/(m3·℃);n為有效孔隙度，無量綱；T1為開采井流體溫度，℃；T0為基準(zhǔn)溫度，一般取恒溫層溫度或多年平均氣溫，℃。

由此得出開采井權(quán)益保護半徑計算公式為

(4)

式中符號含義同公式(3)。

1.3 存在的問題

通過對比公式(1)和公式(4)，發(fā)現(xiàn)公式(1)存在以下3個問題。

(1)公式(1)中f為水比熱與熱儲巖石比熱的比值，即cω/cr，不合適，應(yīng)改為ρωcω與ραcα的比值。

(2)公式(1)中的0.15指的是熱儲回收率，但不同巖性的熱儲回收率取值并不相同，《DZ 40—85 地?zé)豳Y源評價方法》[12]中給出松散巖類孔隙熱儲回收率取值0.25，巖溶裂隙熱儲回收率取值0.15，中生代砂巖和花崗巖等火成巖類熱儲回收率根據(jù)裂隙發(fā)育情況取值0.05～0.1。因此，公式(1)中規(guī)定熱儲回收率取值0.15比較片面。

(3)公式(1)中的36 500指的是100 a的開采天數(shù)，但絕大部分地?zé)峋辉诙竟┡瘯r使用，開采天數(shù)一般為120 d/a，100 a開采天數(shù)應(yīng)為12 000 d。

綜上，公式(4)推導(dǎo)過程合理，表達更為準(zhǔn)確、清晰，可根據(jù)不同熱儲層及具體開采時間，計算合理的開采井權(quán)益保護半徑。因此，建議利用公式(4)進行開采井權(quán)益保護半徑的計算。

2 回灌井權(quán)益保護半徑

進行回灌井權(quán)益保護半徑計算時，在未發(fā)生熱突破的情況下，存在以下兩種假設(shè)情況： 開采井和回灌井相距較近時，可將開采井和回灌井看作一個整體系統(tǒng)； 開采井和回灌井相距較遠時，回灌井可看作一個相對獨立的系統(tǒng)。

2.1 開采和回灌系統(tǒng)

當(dāng)開采井和回灌井相距較近時，可將開采井和回灌井作為一個整體系統(tǒng)，假設(shè)除了開采和回灌的熱量外，系統(tǒng)與外界沒有能量交換。此種情況下，劉志明等[9-10]推導(dǎo)出以下公式

(5)

2.2 回灌單系統(tǒng)

當(dāng)開采井和回灌井井距較大時，可將回灌井看作一個獨立系統(tǒng)。Gringarten等[11]研究發(fā)現(xiàn)，回灌水吸收的熱量等于熱儲層補給的熱量。假設(shè)回灌井向熱儲最高溫度過渡為線性，則存在以下公式

(6)

整理后獲得的公式為

(7)

式中：x為地?zé)崴鞣较蛏夏骋稽c距回灌井底部平面投影的距離，m; 其他符號含義同公式(5)。

2.3 回灌井權(quán)益保護半徑的確定

根據(jù)開采井和回灌井距離的遠近，分別對采、灌系統(tǒng)和回灌單系統(tǒng)進行了回灌井權(quán)益保護半徑計算公式的推導(dǎo)，由于二者界定值的不確定性，按保守性原則，應(yīng)分別按公式(5)、公式(7)進行計算，取較大值作為回灌井權(quán)益保護半徑。二者計算結(jié)果均為正數(shù)，將公式(5)和公式(7)做除法運算，得出以下公式

(8)

式中符號含義同公式(5)。

因此，對于層狀熱儲，當(dāng)(β-1-α)>3λ時，公式(8)大于1，回灌井權(quán)益保護半徑取較大值，應(yīng)按公式(5)計算； 當(dāng)(β-1-α)<3λ時，公式(8)小于1，回灌井權(quán)益保護半徑取較大值，應(yīng)按公式(7)計算。

2.4 開采井與回灌井權(quán)益保護半徑對比

將上述所確定的開采井和回灌井權(quán)益保護半徑計算公式進行對比計算。

跟人學(xué)人樣，跟鬼學(xué)鬼樣。這之后，王祥跟著錢總學(xué)起了古錢幣生意，不過學(xué)的不是識貨認(rèn)價，而是學(xué)起錢總一肚子的壞水。日后時機成熟，這兩人還要在古玩市場里鬧出新花樣。這就是后話了。

(2)公式(7)和公式(4)進行對比。對公式(7)和公式(4)進行除法運算，得出以下公式

。

(9)

式中符號含義同公式(5)。

在地?zé)嵛菜毓噙^程中，回灌量必定≤開采量，則Q回/Q采≤1，λ取值0.15或0.25，公式(9)恒小于1。因此，公式(7)計算得出的回灌井權(quán)益保護半徑R回恒小于公式(4)計算得出的開采井權(quán)益保護半徑R采。

綜上，開采井權(quán)益保護半徑恒大于回灌井權(quán)益保護半徑。

3 計算實例

3.1 研究區(qū)概況

德城區(qū)為山東省德州市的中心城區(qū)，北依京津地區(qū)，南靠濟南，地處天津濱海新區(qū)和環(huán)渤海經(jīng)濟圈，總面積為231 km2。該區(qū)主要開采熱儲為館陶組孔隙型沉積盆地型碎屑巖孔隙-裂隙熱儲，地表無熱流顯示，地?zé)豳Y源類型屬熱傳導(dǎo)型。底板埋深1 200～1 800 m，坳陷區(qū)中心埋深較深，厚度大，四周埋深較淺，厚度薄。熱儲巖性主要為細砂巖、粗砂巖、含礫砂巖和砂礫巖。熱儲厚度為80～120 m，占地層厚度的37%～45%，在取水段1 000～1 500 m深度內(nèi)，單井出水量為70～120 m3/h。熱水礦化度為6～10 g/L，水化學(xué)類型為Cl-Na型，井口水溫為45～65 ℃，屬低溫地?zé)豳Y源中的溫?zé)崴?熱水型地?zé)豳Y源[13-15]。

3.2 計算參數(shù)

根據(jù)巖性，熱儲層巖石的密度和比熱采用《GB/T 11615—2010地?zé)豳Y源地質(zhì)勘查規(guī)范》[8]中給出的經(jīng)驗值。熱儲層厚度、回灌前溫度和抽水量采用山東省地勘局第二水文地質(zhì)工程地質(zhì)大隊在德城區(qū)建設(shè)的砂巖熱儲尾水回灌示范工程(水文家園小區(qū))中獲取的數(shù)據(jù)。具體計算參數(shù)見表1。

表1 地?zé)峋畽?quán)益保護半徑計算參數(shù)

3.3 計算結(jié)果及討論

(1)開采井權(quán)益保護半徑。按照《GB/T 11615—2010地?zé)豳Y源地質(zhì)勘查規(guī)范》[8]提出的公式(1)和本文提出的公式(4)，得出開采井權(quán)益保護半徑分別為1 958.82 m和482.71 m。公式(1)計算結(jié)果遠大于地?zé)峋_采實際影響范圍，且遠大于地?zé)岬V區(qū)批復(fù)范圍； 公式(4)，計算結(jié)果貼近地?zé)峋_采實際影響范圍，比較合理。

(2)回灌井權(quán)益保護半徑。經(jīng)計算，(β-1-α)=0.54(<3λ)。因此，回灌井權(quán)益保護半徑按公式(7)進行計算，結(jié)果為418.04 m。為驗證結(jié)果的準(zhǔn)確性，回灌井權(quán)益保護半徑按公式(5)和公式(7)分別進行計算，結(jié)果為353.59 m和418.04 m，取較大值，為418.04 m。因此，根據(jù)(β-1-α)與3λ的大小關(guān)系確定回灌井權(quán)益保護半徑計算公式的結(jié)論正確。

此外，根據(jù)計算結(jié)果，回灌井權(quán)益保護半徑為418.04 m，小于開采井權(quán)益保護半徑482.71 m。因此，提出的“開采井權(quán)益保護半徑恒大于回灌井權(quán)益保護半徑”的結(jié)論正確。

4 結(jié)論

(1)地?zé)豳Y源開采存在回灌井時，應(yīng)分別計算開采井和回灌井的權(quán)益保護半徑。

(2)參照《GB/T 11615—2010地?zé)豳Y源地質(zhì)勘查規(guī)范》中的規(guī)定，對開采井權(quán)益保護半徑計算公式進行了優(yōu)化，重新給定了參數(shù)f的定義，并將熱儲回收率和開采時間以變量符號替代，可根據(jù)不同熱儲層和具體開采時間加以確定。

(3)對采、灌系統(tǒng)和回灌單系統(tǒng)分別進行了回灌井權(quán)益保護半徑計算公式的推導(dǎo)，按保守性原則，應(yīng)取2個計算結(jié)果的較大值。當(dāng)(β-1-α)>3λ時，回灌井權(quán)益保護半徑按公式(5)進行計算； 當(dāng)(β-1-α)<3λ時，回灌井權(quán)益保護半徑按公式(7)進行計算。

(4)由公式對比與實例計算得出，開采井權(quán)益保護半徑恒大于回灌井權(quán)益保護半徑。