【论文】单丹丹等:增强型地热系统采热性能评价——以共和盆地恰卜恰地区干热岩储层为例

0 　引言

干热岩（ HDR ）地热资源清洁环保可再生，储量丰富 ^[1] ，增强型地热系统（ EGS ）技术是从 HDR 中提取热能的有效手段 ^[2] ，从 EGS 中提取的热能主要用于发电和供暖 ^[3] 。 EGS 是从低渗透性岩体中经济地采出深层地热的人工地热系统 ^[4] 。地热流体的流动过程包含温度场、压力场、力学场和化学场等多物理场的耦合作用 ^[5-6] ，并在时空上发生复杂的变化，数值模拟能够全面精确地反映这种复杂耦合过程 ^[7-14] ，在 EGS 研究中发挥着巨大的作用 ^[15] 。对于完整的 EGS 采热评价，不仅要考虑井筒与地层的传热过程 ^[16] ，还要对裂隙储层进行精细模拟 ^[17] 。对于热储内裂隙岩体热流耦合的数学模型可分为等效连续介质模型与离散裂隙网络模型两大类，后者不仅与现实更接近， 且能更好地模拟水热迁移过程，因而受到广泛关注， 本研究在模拟裂隙热储层的过程中采用了蒙特卡洛方法随机生成的三维离散裂隙网络。与以往不同的是，模拟采用间接耦合的方式，将 EGS 分成注入井、含裸眼段裂隙网络热储层与采出井 3 个部分，分别对其进行模拟计算，以减少计算量，提高计算效率。分 3 段后在注采井井底处，重要模拟参数（质量流量、流体温度）保持一致，以期得到与真实情形更接近的结果，为日后理论研究与工程应用提供参考。

笔者以青海共和盆地恰卜恰干热岩储层为研究对象，建立表征 EGS 水热耦合的数值模型，研究 EGS 开发过程中的温度场与压力场的时空分布特征， 分析 EGS 工程产热特征，同时开展了裂隙属性及操作参数对采热性能的影响研究，得出影响系统产能与寿命的各项因素，以期为商业化开采提供理论依据。

1 　共和盆地地质背景

1.1 　地质概况及勘查现状

共和盆地 ^[18] 所处地区位于青藏高原东北部，青海省东中部，是我国目标地热开发示范区 ^[19] 。大地构造上，经历了早古生代—晚古生代的裂谷坳陷与中生代以来的陆内构造演化两个阶段 ^[20] 。印支期是共和盆地构造转变的重要时期，同时，盆地岩浆侵入活动强烈，在环盆地区出露大量侵入岩，后经燕山期盆地西部隆起和喜马拉雅期青藏高原东北缘差异性隆升作用的叠加，最终形成现今的构造格局 ^[21] 。盆地内的热流值较高，基底花岗岩平均地温梯度大于 5 ℃ /100 m ，热异常明显。早期地球物理资料显示， 共和盆地共划分为塘格木坳陷、贵德坳陷、贵南坳陷、祁家隆起以及黄河隆起 ^[22] 5 个次级构造单元。恰卜恰地热区处于塘格木坳陷与黄河隆起的过渡斜坡带处， 该区地热能属传导型地热，热源来自地下深处，井测温曲线具有线性特征 ^[23] 。目前，恰卜恰地区钻遇干热岩体的勘探井共有 7 口，其中有 4 口井均钻获温度高于 180 ℃的干热岩，再次印证了共和盆地底部存在干热岩且分布广泛 ^[24] 。

选取深度为 2 700 ～ 3 200 m 的裂隙网络储层作为地热能提取的目标储层，假设已经进行较为理想的人工压裂，形成了 500 m × 500 m × 500 m 的裂隙网络储层，裂隙特征依统计参数来定，先给定裂隙数目、裂隙平均迹长、裂隙组方向，再利用蒙特卡洛方法随机生成裂隙网络，为提高计算效率，裂隙假设为无厚度面单元，应用 MATLAB 软件编程实现图形展现，然后采用坐标转换方式将其生成 AutoCAD 脚本文件，进而对生成的裂隙网络进行加工编辑，最终导入 COMSOL 中进行模拟计算，计算时间为 25 年。所选取的裂隙系统模型参数如表 1 所示。系统为双垂直井，注采井井眼半径为 150 mm ，由油管、环空、套管及水泥环构成，其中，油管、套管及水泥参数如表 2 所示，两井间距为 279.11 m 。研究区的 EGS 物理模型如图 1 所示。

层 1 对应于新近系岩层，该层主要由砂岩和砾石持石组成；层 2 代表古近系，该层主要由细砂和粗砂岩组成。基底分为 3 层和 4 层，主要由三叠纪砂质板岩和印支期酸性火山岩组成 ^[25] 。为简化，每一层都看作是单一孔隙率和各向同性的均匀多孔介质，其岩性特征如表 3 所示。一般情况下，岩石的密度随深度的增加而增大，而渗透率和孔隙率则随深度的增加而减小。初始地层条件如下：地表温度， 15 ℃；地表压力， 101 325 Pa ；底部热通量， 0.176 W/m ² ；初始温度， 6.07 z ＋ 15 ℃（ z 表示垂直深度， m ）；初始压力， 0.009 8 z ＋ 0.1 MPa 。

1.3 　数学模型

井筒内流体的能量守恒方程为 ^[26] ：

水流到水泥环的热量损失为 ^[26] ：

井筒总传热系数的表达式为 ^[26] ：

水泥环外缘向地层无限远处的热量损失为 ^[26] ：

无量纲温度 T _D 表达式为 ^[27] ：

岩石基质中的流体流动遵循质量守恒方程 ^[28] ：

裂隙中的流体流动方程描述如下 ^[28] ：

对于传热过程，固体相 s 与流体相 w 的热容和热导率不同，这里假设局部热平衡 ^[29] ，其中引入有效比热容和热导率，并由加权体积和来计算如下 ^[28] ：

同理，裂隙中的能量守恒方程可描述如下 ^[28] ：

方程（ 12 ）、（ 13 ）中的 Φ _T 项代表基岩与裂隙之间的传热量。

在此选用水作为工质流体，是由于其广泛存在性以及高热容性。在热储层高温条件下，水的物理性质是温度的函数。当温度处于 273 ～ 533 K 时，水的物性参数描述如下 ^[30] ：

应用式（ 14 ）～（ 17 ）来计算流体水的物理性质随温度的变化。

2 　边界条件

流量边界条件为由注入井井筒注入的流量值， 设定为 40 kg/s ，压力边界条件为采出井裸眼段压力值，设定为 26 MPa ，裂隙网络储层注入井裸眼段以注入井井底、采出井井底以裂隙网络储层采出井裸眼段的操作参数为准，从而保持整体一致性。边界条件共包含注入井、采出井及裂隙网络储层这三个部分。

2.1 　温度场边界条件

注入井注水温度： T _inj ( t )= T _inj =50 ℃ ；注入井裸眼段水温： T _bor1 ( t )= T _bot1 ( t ) ；裂隙网络储层初始温度： T _lay ( t )= T _ei ( z ) ；采出井井底水温： T _bot2 ( t )= T _bor2 ( t ) 。

2.2 　渗流场边界条件

采出井裸眼段： p _pro ( t )= p _wf =26 MPa ；注入井裸眼段： Q _inj ( t )= Q _inj =40 kg/s ；裂隙网络储层顶边界： p _up =0.009 8 z _on +0.1=0.009 8 × 2 700+0.1=26.6 MPa ；裂隙网络储层底边界： p _un =0.009 8 z _un +0.1=0.009 8 × 3 200+0.1=31.5 MPa 。

若使得地热能发电设备得以发挥最大效率，对于一完整的 EGS ，应设法满足采出井井口温度下降尽量小，同时采出流量尽量大，还要能够维持系统寿命尽可能长，即采出程度尽可能高。因此确定采出井井口温度、采出流量、热提取速率、注入井井底压力与流动阻抗为系统采热性能评价指标。

采出井井口温度计算公式如下 ^[17] ：

热提取速率计算公式如下 ^[12] ：

在干热岩工程中，常用流动阻抗值来评价储层的渗流能力，其定义是指驱动单位质量流体流经热储所需消耗的能量，计算公式如下 ^[12] ：

3 　模拟结果与分析

图 2 给出的是储层裂隙面温度随时间变化云图，可知储层裂隙面温度降低区由注入井逐渐向采出井蔓延，深度方向上，储层底面压力最大，渗流流量最大，与基岩热交换速度也最快，故使得底面裂隙面温度降低最快，从而引起底面处裸眼段采出井采出温度降低很快，影响整体采热温度。图 3 是储层整体压力随时间变化云图，由图可见，储层整体压力场不再随深度均匀增大，注入井端压力大于采出井，在注采井间形成压力梯度，导致水流由注入井流向采出井，以实现热能开采。考虑顶底面流量损失后，在注采井压差及顶底面压差作用下，会引起水流经顶底面渗流流出热储层，进而经采出井采出流量减少，造成储层流量与热量损失。

由图 4-a 采出井井底与井口温度对比曲线可见，开采 5 年后，井底温度 178.56 ℃，井口则 170.52 ℃，二者相差 8.04 ℃，若不考虑井筒热损失，采出温度会被高估 5% ；开采 10 年后，井底温度 171.33 ℃，井口则 164.22 ℃，相差了 7.11 ℃，与不考虑井筒热损失相比，也将被高估 4% ；可见，考虑初前期采出井井筒热损失对于深层地热能开采更具有现实指导意义。图 4-b 为采出井井口水温与热焓变化曲线，可知，二者变化规律相同，都随时间的增长而降低，且降低的趋势是先快再有所减慢最后再加快，这是由于热焓是温度的函数。图 4-c 是采出流量及热提取速率变化曲线，二者都随时间增长而减小，且 5 年内降低最快，时间越长降低越缓慢，由图可见，开采 5 年后流量由 23.80 kg/s 降至 20.12 kg/s ， 10 年后降至 19.55 kg/s ， 25 年后降至 18.51 kg/s ，与注入流量 40 kg/s 相比，降幅超一半，说明不考虑储层流量与热量损失必然会导致采出流量过高评价，热提取速率也过高评价。这里考虑流量、热量损失后，热提取速率由最初的 14.50 MW 降低至 6.90 MW 。图 4-d 为注入井井底压力与流动阻抗随时间的变化规律，由图可知，二者都随时间的增长而增大，且开始增大较快，而后逐渐趋于稳定，呈现指数分布规律。

4 　系统采热性能影响因素分析

4.1 　注入流量

图 5 为不同注入流量下各数据结果对比分析曲线。从图可知，对于开采 15 年以前，随注入流量的增大，井口采出温度并未降低，但采出流量与热提取速率均增大，这既保证了采出高温热流体又提高了热提取速率，从而有利于热能的提取，但注入流量增大会引起注入压力增大，从而使流动阻抗增大，进而影响系统效率。对于开采 15 年以后，井口采出温度随注入流量的增大而降低较快，这时会缩短储层寿命。因此，提高注入流量，会在整个开采期限内持续一段时间维持高温高提取率开采，但之后将会出现系统效率降低，储层寿命降低的情况。

4.2 　生产压力

为探究生产压力的变化对高温热储层开采效果的影响，这里选取 3 种生产压力工况进行模拟研究，分别为 p _pro =22 MPa 、 p _pro =24 MPa 以及 p _pro =26 MPa ，其余模拟参数保持不变，经计算，得出了如下图 6 所示的不同采出井压力 EGS 采热性能评价指标曲线图。从图可知，生产压力的变化对采出井井底与井口温度的影响较大，随着生产压力的增大，井底与井口温度也越高，井口热焓也随之增大，生产压力每提高 2 MPa ，将导致井口温度在热开采 25 年后提高 10 ℃左右，井口热焓提高 35 kJ/kg 以上。生产压力变化对采出流量的影响也较大，随着生产压力的增大，采出流量不断减小，说明流量损失越多。井口温度与采出流量共同影响着热提取速率，由于流量引起的热提取速率降低速度，要大于井口温度引起的热提取速率升高速度，从而使得随着生产压力增大，热提取速率降低，对于开采初期，降低的更快。生产压力的变化对注入井井底压力的影响非常小，而对流动阻抗的影响却较大，原因在于采出流量受生产压力的影响较大，生产压力越大，采出流量越小，流动阻抗就会越大。

以裂隙渗透率为对比分析工况进行研究，这里选取3种不同裂隙渗透率工况来模拟计算，分别为 k _f =0.5×10 ^{－ 10} m ² 、 k _f =1.0×10 ^{－ 10} m ² 以及 k _f =2.0×10 ^{－ 10} m ² ， 对于其他模拟参数均保持不变，经计算，得出如下图7 所示的不同裂隙渗透率EGS 采热性能评价指标曲线图。可以得知，裂隙渗透率变化，则井底与井口温度、井口热焓都有较大变化，随着裂隙渗透率的增大，井底与井口温度、井口热焓则降低，且随裂隙渗透率的增大，温度与热焓的降低程度越大。裂隙渗透率增大，则采出流量与热提取速率也增大， 但热提取速率是由采出流量与井口热焓共同决定的， 且二者对热提取速率的影响效果又是相反的，因此会导致热提取速率的变化开始较大，随着时间的增长， 逐渐减小，进而使系统寿命缩短。裂隙渗透率增大， 注入井井底压力会减小，由注入井向采出井流动的水流速度会加快，从而使得水流阻抗也相应减小。

4.4 　裂隙开度

为了分析裂隙开度的变化对储层开采效果的影响，这里选取 3 种裂隙开度的工况进行对比研究， 分别为 d _f =1.5 mm 、 d _f =2.0 mm 以及 d _f =2.5 mm ，其余模拟参数均保持不变，对其进行模拟计算，最终得出了如下图 8 所示的不同裂隙开度 EGS 采热性能评价指标曲线图。由图可知，对于开采初期，裂隙开度变化几乎对采出温度没有影响，在开采 5 年后，温度曲线开始出现变化，当裂隙开度增大时，储层采出温度随之降低，对于井口热焓而言，其变化规律与温度相一致，温度总体上是采出井井底裸眼高于井口，采出井处于热量损失状态，对于采出流量，当裂隙开度增大时，采出井采出流量也增大，并且裂隙开度越大，这种增大效应越明显，当裂隙开度由 d _f =1.5 mm 增大至 d _f =2.0 mm 时，储层的采出流量增大不到 0.5 kg/s ，而当裂隙开度由 d _f =2.0 mm 增大至 d _f =2.5 mm 时，储层的采出流量却增大 1 kg/s 以上。井口热焓与采出流量对热提取速率的影响相反，在二者的共同作用下，热提取速率变化不大。裂隙开度的微小变化对注入井井底压力与流动阻抗的影响也不大。

5 　结论

1）不考虑注采井井筒部分热损失，将会使开采初前期采出温度被高估4% ～ 5% ；考虑了顶底面的流量与热量损失后，会使储层较绝热状态下的采热温度相比有所降低；同时，在注采井及顶底面压力差的作用下，地层流场会发生变化，一部分水流经顶底面渗流会流出热储层，从而经采出井采出的流量会减少，造成储层流量与热量的损失。

2 ）对比分析注入流量、生产压力的变化对 EGS 各采热性能评价指标的影响，得出注入流量增大会使系统维持一段时间高温高热提取率开采，但之后储层寿命将降低；生产压力增大会降低热提取速率， 且流动阻抗增大，不利于热开采。

3 ）对比分析裂隙渗透率、裂隙开度的变化对 EGS 各采热性能评价指标的影响，得出裂隙渗透率对各产能评价指标都有较大影响，且渗透率越大，影响越明显，裂隙渗透率增大虽能提高热提取速率，却会引起储层寿命降低；当裂隙开度变化较小时，对系统采热效果影响较小，主要在于裂隙开度的细微变化对裂隙渗流流量的影响不大。