張 霄，孫玉學(xué)，張慶松，李相輝，尹占超，李 壯，周 政

(山東大學(xué) 巖土與結(jié)構(gòu)工程研究中心，山東 濟(jì)南 250061)

近年來(lái)，化石資源臨近枯竭，環(huán)境問(wèn)題日益突出[1]。干熱巖資源因儲(chǔ)量巨大、清潔環(huán)保，決定了其必然會(huì)成為未來(lái)重要能源之一[2]。增強(qiáng)型地?zé)嵯到y(tǒng)（EGS）作為開(kāi)發(fā)干熱巖資源的有效手段，受到了廣泛關(guān)注[3-5]。

研究EGS取熱過(guò)程時(shí)，若利用建造干熱巖電站的方法獲得實(shí)際工程的地質(zhì)資料、工程經(jīng)驗(yàn)和數(shù)據(jù)，對(duì)研究干熱巖儲(chǔ)層熱提取極為有利，但目前的干熱巖儲(chǔ)層改造技術(shù)尚不完善，投資巨大、周期長(zhǎng)。因此，現(xiàn)階段研究大多基于數(shù)值模擬方法建立數(shù)學(xué)模型，利用多場(chǎng)耦合過(guò)程對(duì)EGS取熱過(guò)程進(jìn)行評(píng)估，指導(dǎo)EGS的設(shè)計(jì)與優(yōu)化[6-8]。在熱儲(chǔ)改造和EGS運(yùn)行過(guò)程中，熱儲(chǔ)巖體的物理化學(xué)性質(zhì)都會(huì)隨時(shí)間和空間不斷變化[9-10]，改造后熱儲(chǔ)尺寸多以km為量級(jí)，內(nèi)部裂隙通常以mm為量級(jí)，在利用精細(xì)化網(wǎng)格直接模擬取熱過(guò)程時(shí)計(jì)算量巨大[11]，需要對(duì)熱儲(chǔ)內(nèi)的裂隙進(jìn)行簡(jiǎn)化。目前，常用的簡(jiǎn)化模型有單孔隙模型[12]、雙孔隙模型[13-14]和多孔隙模型。其中，雙孔隙模型和多孔隙模型因精度要求較高，需要大量試驗(yàn)數(shù)據(jù)支撐，在實(shí)際中應(yīng)用困難[15]。

干熱巖熱量資源預(yù)測(cè)方面，目前主要有兩種預(yù)測(cè)方式：一種是利用目標(biāo)熱儲(chǔ)體積及熱儲(chǔ)溫度和某一確定溫度（如地表溫度）的差值計(jì)算，按百分比的可開(kāi)采量預(yù)測(cè)[16-19]，主要有體積法、蒙特卡羅法及地表流量法等[20-23]，但這些方法在應(yīng)用于EGS產(chǎn)能預(yù)測(cè)時(shí)，熱儲(chǔ)開(kāi)采前后溫差的設(shè)定值偏大，導(dǎo)致預(yù)測(cè)不準(zhǔn)確；另一種是數(shù)值法[24]，根據(jù)詳細(xì)的熱儲(chǔ)參數(shù)計(jì)算和評(píng)價(jià)勘查程度較高、具有一定開(kāi)采歷史的地?zé)崽颷23]，此種方法較為繁瑣，計(jì)算周期長(zhǎng)，需要詳細(xì)的熱儲(chǔ)參數(shù)，不適合未開(kāi)采的熱儲(chǔ)產(chǎn)能的預(yù)測(cè)。在EGS運(yùn)行系統(tǒng)可控因素的研究方面，目前的研究多集中于以不同的熱媒介質(zhì)注入溫度、不同的熱媒介質(zhì)注入速度及不同井間距等人為可控因素控制EGS的運(yùn)行[8,25]。以注水速度控制時(shí)，隨著EGS的運(yùn)行，注入井的壓力不斷產(chǎn)生變化，大大增加了EGS運(yùn)行時(shí)的不可控性[26]。EGS取熱效果方面，大多數(shù)研究將熱儲(chǔ)視為單孔隙模型，針對(duì)熱儲(chǔ)改造后的孔隙比、滲透率、注入井溫度、布井方式等多種因素分析EGS取熱效果[27-28]。如：王昌龍[29]考慮了11種因素對(duì)EGS取熱效果的影響。在這些研究中，大多局限于單因素的定性研究，較少考慮多因素組合的定量研究，且鮮有將大地?zé)崃鞯臒嵫a(bǔ)償作為初始條件加入到模型之中[30-33]。有研究發(fā)現(xiàn)，熱儲(chǔ)周?chē)拇蟮責(zé)崃鲗?duì)熱儲(chǔ)層長(zhǎng)期積累的熱補(bǔ)償熱量會(huì)對(duì)EGS運(yùn)行壽命產(chǎn)生不可忽略的影響[34]。

以吸放熱公式[35]和Dupuit公式為基礎(chǔ)，建立EGS產(chǎn)能和壽命控制方程；依據(jù)產(chǎn)能和熱儲(chǔ)參數(shù)，通過(guò)調(diào)節(jié)熱媒介質(zhì)注入溫度和注采壓差兩個(gè)可控因素取值預(yù)測(cè)EGS壽命，增加EGS運(yùn)行時(shí)的可控性。同時(shí)，根據(jù)地下水流動(dòng)方程和熱傳導(dǎo)方程分析出熱儲(chǔ)初始溫度、熱媒介質(zhì)注入溫度、EGS壽命、熱儲(chǔ)體積和熱儲(chǔ)比表面積5個(gè)因素對(duì)發(fā)生熱突破（生產(chǎn)井產(chǎn)出溫度低于423.16 K[36]）時(shí)對(duì)控制方程中熱儲(chǔ)減少熱量、大地?zé)嵫a(bǔ)償熱量、生產(chǎn)井平均產(chǎn)出溫度和熱儲(chǔ)形狀系數(shù)4個(gè)未知參數(shù)的影響。在考慮大地?zé)嵫a(bǔ)償?shù)那闆r下，利用雙井EGS的3維數(shù)值模型，對(duì)不同工況的EGS在發(fā)生熱突破時(shí)的4個(gè)未知參數(shù)的變化規(guī)律進(jìn)行研究，得到4個(gè)未知參數(shù)隨5個(gè)影響因素變化的具體預(yù)測(cè)公式，實(shí)現(xiàn)EGS產(chǎn)能預(yù)測(cè)和多因素定量化分析。

1 流熱耦合下EGS換熱過(guò)程分析

如圖1所示，假設(shè)熱儲(chǔ)周?chē)鷩鷰r滲透率為0，熱媒介質(zhì)在流動(dòng)過(guò)程中無(wú)熱量損失，人工儲(chǔ)層原有熱量Qm和大地補(bǔ)償?shù)臒崃縌t完全由熱媒介質(zhì)攜帶并提取，其中，熱媒介質(zhì)攜帶熱量由Qw表示。

圖1 EGS換熱過(guò)程的流熱耦合機(jī)制Fig.1 Flow-heat coupling mechanism of EGS heat transfer process

根據(jù)熱量平衡，由吸放熱公式（式（1））結(jié)合Dupuit公式（式（2））可得熱媒介質(zhì)攜帶熱量（EGS產(chǎn)出熱量），如式（3）所示：

式（1）～（3）中：t為EGS壽命；k為熱儲(chǔ)滲透率；mw為熱媒介質(zhì)質(zhì)量； ρw為熱媒介質(zhì)密度；L為熱媒介質(zhì)路徑；A1為熱媒介質(zhì)在熱儲(chǔ)中流過(guò)的截面面積；cw為熱媒介質(zhì)的熱容；vw為熱媒介質(zhì)平均流速；P為注采壓差；ΔT1為注采溫差，ΔT1=Tout-Tin。

若熱儲(chǔ)在tn時(shí)刻發(fā)生熱突破，tm時(shí)刻的產(chǎn)出溫度Tout(tm)可作為T(mén)out(t)在t0至tn時(shí)刻的平均值，即：

整理式（3）和（4），可得EGS產(chǎn)能和壽命控制方程：

利用達(dá)西定律和多孔介質(zhì)傳熱方程描述流熱耦合過(guò)程，如式（6）～（9）所示：

連續(xù)性方程：

流體動(dòng)量方程：

流體能量方程：

周?chē)鷩鷰r熱補(bǔ)償能量方程：

牛頓冷卻公式：

式（5）～（10）中，t為熱媒介質(zhì)流動(dòng)時(shí)間，xi為熱媒介質(zhì)在i方向上的位移，ρw、ρm分別為流體、圍巖的熱儲(chǔ)密度，vi為熱儲(chǔ)中流體的流速，F(xiàn)i為體積力， ?2為拉普拉斯算子，cw、cm為流體和圍巖比熱容，A2為熱儲(chǔ)與圍巖接觸面積，η為動(dòng)力黏度，T為熱儲(chǔ)溫度場(chǎng)，T∞為圍巖初始溫度，λw、λp分別為流體和圍巖的導(dǎo)熱系數(shù)，h為表面?zhèn)鳠嵯禂?shù)，Φ為熱流量，A3為流體與固體接觸面積，Ts、Tw分別為固體和流體溫度。

利用式（6）～（9），結(jié)合牛頓冷卻公式（10），對(duì)式（5）中的4個(gè)參數(shù)進(jìn)行分析，可確定出影響4個(gè)參數(shù)的主要因素，如表1所示。

表1 4個(gè)未知變量的影響因素匯總Tab. 1 Summary of influencing factors of four unknown parameters

Qm由熱儲(chǔ)終溫Tml和熱儲(chǔ)初始溫度Tm控制，所以在預(yù)測(cè)時(shí)僅需考慮Tml即可，整理后如式（11）所示。

若利用式（5）達(dá)到討論Tin和P與t之間關(guān)系的目的，可利用數(shù)值模擬預(yù)測(cè)4個(gè)參數(shù)的具體表達(dá)式，并代入式（5）中得到Tin、P、t三者關(guān)系。

2 EGS取熱模型概述

2.1 地質(zhì)模型構(gòu)建

如圖2所示，EGS地下部分由人工熱儲(chǔ)、熱儲(chǔ)周?chē)膰鷰r、注入井和生產(chǎn)井組成。為了更準(zhǔn)確地模擬EGS地下熱交換過(guò)程，在模型中作出主要假設(shè)：1）熱儲(chǔ)裂隙小而均勻，可看做單孔隙率的多孔介質(zhì)；2）熱媒介質(zhì)流動(dòng)為單相流；3）熱儲(chǔ)中熱媒介質(zhì)流動(dòng)是飽和層流流動(dòng)（0.003 m/s以下），初始時(shí)刻熱儲(chǔ)層中充滿(mǎn)了與周?chē)鷰r石溫度相同的熱媒介質(zhì)；4）周?chē)鷩鷰r滲透率為0，熱媒介質(zhì)在流動(dòng)過(guò)程中無(wú)損失；5）熱儲(chǔ)及周?chē)鷩鷰r不會(huì)產(chǎn)生形變和各種化學(xué)反應(yīng)。

圖2 EGS地下部分概念圖Fig.2 Concept picture of EGS underground part

圖3 雙井EGS模型3維圖Fig.3 Three-dimensional pictures of the double well EGS model

如圖3所示，熱儲(chǔ)周?chē)?00 m厚的花崗巖圍巖，在外表面施加以3 K/100 m隨深度線性增加的溫度邊界，模擬大地?zé)崃鹘o予熱儲(chǔ)的熱量補(bǔ)給；分別建立不同邊長(zhǎng)的正六面體，研究體積對(duì)4個(gè)參數(shù)的影響，如圖3（a）所示；分別建立同體積下不同形狀的熱儲(chǔ)，研究比表面積對(duì)4個(gè)參數(shù)的影響，如圖3（b）所示；其余模型邊界條件如表2所示。模型采用“下注上采”的開(kāi)采方式，最大模擬時(shí)間50 a，采用瞬態(tài)求解方式，計(jì)算結(jié)果每年采集一次，且最終結(jié)果取符合EGS運(yùn)行25 a，生產(chǎn)井出水速率達(dá)到80 kg/s以上，保證EGS商業(yè)化利用[35]。

表2 模型邊界條件取值Tab. 2 Values of model boundary conditions

2.2 模型熱物性參數(shù)取值

模型其他熱物性參數(shù)如表3所示[30-31]。

表3 模型熱物性參數(shù)取值Tab. 3 Values of model thermal property parameters

2.3 網(wǎng)格獨(dú)立性驗(yàn)證

本文主要對(duì)計(jì)算結(jié)果中的生產(chǎn)井平均產(chǎn)出溫度和平均速度、熱儲(chǔ)平均溫度和熱儲(chǔ)中流體平均速度進(jìn)行獨(dú)立性檢驗(yàn)。

進(jìn)行網(wǎng)格無(wú)關(guān)性檢驗(yàn)時(shí)，按時(shí)間步長(zhǎng)0.01 s，對(duì)不同自由四面體的網(wǎng)格數(shù)量進(jìn)行無(wú)關(guān)性檢驗(yàn)。采用熱儲(chǔ)邊長(zhǎng)為500 m正六面體，熱儲(chǔ)初始溫度-熱媒介質(zhì)注入溫度-注采壓差為565.66 K-353.16 K-20 MPa的工況，分別選取第15年和第30年的計(jì)算結(jié)果進(jìn)行比較（圖4）。由圖4可知：網(wǎng)格數(shù)量對(duì)計(jì)算結(jié)果有著較大的影響，當(dāng)網(wǎng)格數(shù)量小于25萬(wàn)個(gè)時(shí)，計(jì)算結(jié)果存在較大差異；當(dāng)網(wǎng)格數(shù)量大于30萬(wàn)個(gè)時(shí)，模擬結(jié)果變化不大。

圖4 網(wǎng)格無(wú)關(guān)性檢驗(yàn)結(jié)果Fig.4 Results of the grid independence verification

進(jìn)行時(shí)間獨(dú)立性檢驗(yàn)時(shí)，以第12種（36.76萬(wàn)個(gè)）網(wǎng)格數(shù)量為網(wǎng)格模型，時(shí)間步長(zhǎng)為0.007 5、0.010 0、0.012 5和0.015 0 s，其余條件與網(wǎng)格獨(dú)立性檢驗(yàn)時(shí)相同，結(jié)果如表4所示。由表4可知，當(dāng)網(wǎng)格數(shù)量達(dá)到30萬(wàn)以后，時(shí)間步長(zhǎng)在0.012 5 s以下時(shí)對(duì)模擬結(jié)果影響不大。

由圖4和表4可知，在計(jì)算時(shí)，最大網(wǎng)格取40 m，最小網(wǎng)格取0.3 m，網(wǎng)格增長(zhǎng)率取2.4，單元格數(shù)量41.21萬(wàn)個(gè)，時(shí)間步長(zhǎng)選擇0.01 s較為合理。

表4 時(shí)間獨(dú)立性檢驗(yàn)結(jié)果Tab. 4 Results of the time independence verification

3 結(jié)果分析

圖5 Tm、Tin、V、SSa、t對(duì)Tml、Qt、、 α等參數(shù)的影響結(jié)果Fig.5 Influences of five factors (Tm、Tin、V、SSa、t) to four parameters (Tml、Qt、 Tout、 α)

由圖5（a）、（c）、（d）、（e）可知：當(dāng)Tm和Tin增加時(shí)，Tml逐漸上升，且Tm對(duì)Tml的影響要低于Tin對(duì)Tml的影響；隨著SSa的增加，Tml也會(huì)逐漸升高；隨著V的增大，Tml表現(xiàn)出逐漸減小的趨勢(shì)；以Tm=565.66 K、Tin=293.16 K為例，EGS不同壽命所對(duì)應(yīng)的Tml相差不大，近似服從高斯分布，與時(shí)間項(xiàng)無(wú)關(guān)。由于熱突破條件并非與熱媒介質(zhì)注入溫度相同，所以在注入井與生產(chǎn)井之間必定存在溫度梯度，在均勻介質(zhì)滲流場(chǎng)中，兩井之間必存在流速最快的最優(yōu)通路。在最優(yōu)通路間，溫度梯度介于注入溫度和熱突破溫度之間，在最優(yōu)通路外其溫度梯度上限與熱儲(chǔ)初始溫度有關(guān)。因此，隨著Tm的升高，最優(yōu)通路外的熱儲(chǔ)溫度升高；隨著Tin的升高，整個(gè)熱儲(chǔ)溫度梯度下限升高；隨著SSa的增加，熱儲(chǔ)與圍巖接觸面積增加，熱補(bǔ)償?shù)臒崃吭龆?，都?huì)導(dǎo)致Tml的升高。此外，其他條件相同時(shí)，V的增加使得熱突破時(shí)低溫區(qū)域增大，Tml逐漸減小。

由圖5（f）、（g）、（h） 可知：受?chē)鷰r熱補(bǔ)償?shù)挠绊懀?5年后 ，Qt隨著EGS運(yùn)行時(shí)間逐漸增加；隨著Tm、V和SSa的增大，Qt也逐漸增加；對(duì)于Tin對(duì)Qt的影響，由圖5（f）分析可以發(fā)現(xiàn)，除熱儲(chǔ)初始溫度處于485.66 K時(shí)，Qt隨熱媒介質(zhì)注入溫度逐漸增加，其余工況均無(wú)明顯規(guī)律。分別取第27年和第29年Qt的值，利用極差分析法將Tm與Tin對(duì)Qt的影響進(jìn)行比較，如表5和6所示，發(fā)現(xiàn)熱儲(chǔ)初始溫度Tm對(duì)Qt的影響要遠(yuǎn)遠(yuǎn)大于Tin對(duì)Qt的影響，可近似將Tin對(duì)Qt的影響忽略。因此，可以確定，Qt主要受到Tm、Tin、V和SSa的影響。

表5 第27年Tm與Tin極差分析Tab. 5 Range analysis of Tm and Tin in the 27th year

表6 第29年Tm與Tin極差分析Tab. 6 Range analysis of Tm and Tin in the 29th year

由圖5（g）、（h）可知：隨著V的增大，α表現(xiàn)出先減小后增大的趨勢(shì)，其最低點(diǎn)大致出現(xiàn)在邊長(zhǎng)為400 m的正六面體熱儲(chǔ)附近；隨著SSa的增大，α表現(xiàn)出先增大后減小的趨勢(shì)，其最高點(diǎn)出現(xiàn)在正六面體熱儲(chǔ)附近。

經(jīng)上述分析，將表1中的影響因素予以修正，如表7所示，同時(shí)將式（11）修正為式（12）：

表7 修正后的Tml、Qt、、α影響因素匯總Tab. 7 Summary of revised influencing factors of Tml, Qt,,α

表7 修正后的Tml、Qt、、α影響因素匯總Tab. 7 Summary of revised influencing factors of Tml, Qt,,α

參數(shù) 影響因素Tml Tm、Tin、V、SSa Qt Tm、V、SSa、t Tout Tm、Tin、V、SSa α V、SSa

為得到式（12）的具體表現(xiàn)形式，需要對(duì)上述計(jì)算結(jié)果繼續(xù)分析。

3.1 熱儲(chǔ)終溫Tml 變化規(guī)律

為探究表7中4個(gè)因素與Tml的定量關(guān)系，以Tin為353.16 K時(shí)為基準(zhǔn)，確定出Tml353.16隨Tm的變化；然后，基于此分別討論Tin、SSa與V對(duì)Tml的增量的變化，如圖6所示；最后，將上述4項(xiàng)進(jìn)行線性疊加，如式（13）所示：

如圖6（a）所示，極易確定出Tml353.16隨Tm變化公式。如圖6（b）所示，不同Tm的ΔTmlTin隨著Tin呈一致的線性負(fù)相關(guān)關(guān)系，因此，ΔTmlTin不受Tm的影響。取6種不同Tm的ΔTmlTin的平均值，得到ΔTmlTin隨Tin的變化公式，如圖6（c）所示。由于比表面積在不同體積時(shí)難以一致表示，因此，預(yù)測(cè)比表面積引起的熱儲(chǔ)終溫增量ΔTmlS時(shí)，應(yīng)以在同體積下與比表面積變化規(guī)律一致的量—各表面積與相應(yīng)體積下正六面體表面積的比值S/SR表示，以消除熱儲(chǔ)體積對(duì)ΔTmlS的影響。如圖6（d）所示，得到了ΔTmlS隨S/SR的變化規(guī)律。考慮ΔTmlV隨不同體積熱儲(chǔ)與邊長(zhǎng)為500 m的正六面體熱儲(chǔ)體積比值V/V500的變化規(guī)律，用以量化體積對(duì)Tml的影響，如圖6（e）所示。所以得到式（13）中各分項(xiàng)的具體表達(dá)式為：

圖6 熱儲(chǔ)終溫Tml計(jì)算圖Fig.6 Calculation of final storage temperature (Tml)

3.2 大地?zé)嵫a(bǔ)償熱量 Qt變化規(guī)律

首先，得到前25年的補(bǔ)償熱量隨熱儲(chǔ)初始溫度Tm的變化量Qt25，基于此分別討論Qt隨t和SSa的增量變化ΔQtt和ΔQtS，將上述3個(gè)量進(jìn)行線性疊加。如圖5（g）所示，當(dāng)熱儲(chǔ)體積為0時(shí)，Qt也一定為0。令不同體積得到的補(bǔ)償熱量與邊長(zhǎng)為500 m的正六面體熱儲(chǔ)體積得到的補(bǔ)償熱量的比值為ηV，將ηV作為體積系數(shù)與Qt25、ΔQtt和ΔQtS疊加得到的補(bǔ)償熱量相乘，以考慮體積因素對(duì)Qt的影響，如式（15）所示：

Qt25與Tm的關(guān)系、ΔQtS與S/SR的關(guān)系及ηV與V/V500的關(guān)系如圖7（a）、（b）和（d）所示。由圖7（e）可知：ΔQtt隨著時(shí)間的延長(zhǎng)不斷增加，近似呈線性變化；但不同Tm對(duì)應(yīng)的ΔQtt的增量不同，且隨著Tm的上升，ΔQtt隨著時(shí)間變化的線性曲線斜率不斷增加，截距不斷減小。故在考慮ΔQtt隨著時(shí)間的變化時(shí)，還要考慮Tm的影響。分別計(jì)算出斜率c和截距d隨Tm的變化，如圖7（c）所示。最終得到式（15）中各分項(xiàng)的具體表達(dá)式為：

圖7 大地?zé)嵫a(bǔ)償熱量Qt計(jì)算Fig.7 Calculation of geothermal compensation heat (Qt)

3.3 平均產(chǎn)出溫度 Tout變化規(guī)律

式中，Tout353.16、ΔToutS和ΔToutV極易通過(guò)數(shù)值計(jì)算結(jié)果通過(guò)擬合曲線得到，如圖8（a）、（c）和（d）。如圖8（e）所示，隨著Tin的減小，ΔToutTin逐漸增加，可近似看做線性曲線，且不同Tm的ΔToutTin在同一注入溫度的增量不同，隨著注入溫度的上升而不斷增加。因此，在計(jì)算ΔToutTin時(shí)，需要考慮Tm的影響。與ΔQtt計(jì)算方法類(lèi)似，如圖8（b）所示，分別計(jì)算出ΔToutTin隨Tin變化線性曲線的斜率a和截距b隨Tm的變化，即可得到ΔToutTin的計(jì)算值，則式（17）中各分項(xiàng)的具體表達(dá)式為：

圖8 平均產(chǎn)出溫度 計(jì)算Fig.8 Calculation of average output temperature

3.4 系數(shù) α的變化規(guī)律

由表7可知，系數(shù)α只與V和SSa有關(guān)。為得到三者間的定量關(guān)系，將α分為αV和ΔαS。首先，得到αV隨V的變化，如圖9（a）所示；然后，基于αV計(jì)算由SSa引起的增量ΔαS，如圖9（b）所示；最后，將αV與ΔαS進(jìn)行線性疊加得到系數(shù)α的計(jì)算式，如式（19）所示：

如圖9所示，最終得到式（19）中各分項(xiàng)的具體表達(dá)式為：

圖9 系數(shù) α的計(jì)算圖Fig.9 Calculation of coefficientsα

式（21）中各分式可由式（14）、（16）、（18）、（20）表示。將式（21）代入式（5）中，結(jié)合勘測(cè)到的熱儲(chǔ)滲透率即可得到Tin、P、t三者關(guān)系，指導(dǎo)EGS建設(shè)。

3.5 計(jì)算方程的適用性驗(yàn)證

目前，大多數(shù)研究主要利用工質(zhì)循環(huán)流量表示熱媒介質(zhì)注入的量及熱媒介質(zhì)產(chǎn)出量，利用注采壓差P與已有研究比較來(lái)驗(yàn)證計(jì)算公式適用性的方法難以直接實(shí)現(xiàn)。根據(jù)Dupuit公式可知，通過(guò)熱儲(chǔ)尺寸和滲透率可將注采壓差P換算為熱媒介質(zhì)在生產(chǎn)井的平均流速，由生產(chǎn)井的尺寸即可得到熱媒介質(zhì)產(chǎn)出量，即工質(zhì)循環(huán)流量。為驗(yàn)證計(jì)算方程的可靠性，將其應(yīng)用至文獻(xiàn)[15]的模型中。將文獻(xiàn)[15]中的熱儲(chǔ)尺寸、熱儲(chǔ)初始溫度、熱媒介質(zhì)注入溫度和不同EGS運(yùn)行時(shí)間代入式（5）中，得到不同EGS運(yùn)行時(shí)間對(duì)應(yīng)的注采壓差P；然后，僅使用 COMSOL Multiphysics中達(dá)西定律模塊計(jì)算得到工質(zhì)循環(huán)流量，所得結(jié)果如圖10所示。研究發(fā)現(xiàn)，利用式（5）計(jì)算結(jié)果與文獻(xiàn)[15]所得曲線結(jié)果一致，說(shuō)明式（5）具有良好的適用性。

圖10 計(jì)算公式適用性驗(yàn)證Fig.10 Applicability verification of calculation formula

4 結(jié) 論

利用理論分析與數(shù)值模擬相結(jié)合，描述EGS地下取熱過(guò)程，主要結(jié)論如下：

1）將吸放熱公式與Dupuit公式相結(jié)合，得到了雙井EGS產(chǎn)能與壽命控制方程，揭示了EGS產(chǎn)能隨熱儲(chǔ)滲透率、EGS壽命、熱媒介質(zhì)注入溫度、注采壓差和熱儲(chǔ)形狀系數(shù)的變化關(guān)系。

3）將控制方程和預(yù)測(cè)公式與現(xiàn)有研究進(jìn)行對(duì)比，控制方程和預(yù)測(cè)公式具有良好的適用性，可有效指導(dǎo)EGS建設(shè)，加強(qiáng)EGS運(yùn)行期間的實(shí)時(shí)調(diào)控。