comsol 裂隙 耦合_单丹丹等：随机裂隙网络储层与井筒热流耦合数值模拟

本文版权归天然气工业杂志社所有

未经允许，不得转载

单丹丹^1,2　闫    铁^1,2　李   玮^1,2

孙士慧^1,2　逯广东^1,2　赵   欢^1,2

1. 东北石油大学石油工程学院

2. 提高油气采收率教育部重点实验室·东北石油大学

摘　要    增强型地热系统(Enhanced Geothermal System，EGS)的渗流传热通道主要由注采井井筒与随机裂隙网络储层两大部分组成，过去对裂隙网络储层进行热流耦合模拟研究大多数都忽视了井筒壁的换热，导致模拟结果的准确性欠佳。为了更准确地评价EGS 出力、寿命等性能指标，开展了井筒—随机裂隙网络储层热流耦合数值模拟研究，基于商业有限元软件COMSOLMultiphysics 对井筒与储层的渗流场、温度场进行耦合求解，分析了影响EGS 采出温度与热开采速率的各项要素。研究结果表明：①注、采井的开孔长度(L₀)对EGS 的产能与寿命具有重要影响，400 m 为最佳开孔长度， 其EGS 具有最佳出力与寿命；②在井筒壁上加保温材料可以有效提高开采初、前期的采出温度，减少热损失，提高开采速率；③随着开采的进行，注入井周围出现明显的低温区，并沿裂隙通道向采出井方向推移，这将导致系统达到开采寿命而衰竭，应停止开采一段时间后再进行热能开采；④裂隙渗透率、裂隙宽度等参数对开采速率的影响都呈现正相关性，参数值增大会提高开采速率、缩短开采年限。结论认为，井筒壁的换热对于EGS 出力与寿命的综合评价具有重要的意义，考虑井筒壁热损失的井筒—热储耦合模拟能够实现对EGS 的完整性评价。

关键词   增强型地热系统　井筒　随机裂隙网络　出力与寿命　热流耦合　COMSOL Multiphysics　采出温度　开采速率

0　引言

干热岩是一种深层热岩体，其资源储量极其丰富^[1-2]。增强型地热系统(Enhanced GeothermalSystem，EGS)是目前开采干热岩资源最有效的技术手段^[3-4]。EGS 的设计与优化主要依据两方面重要因素——出力与寿命^[5]。对系统出力与寿命进行评价，要反复模拟EGS 热流耦合过程。针对这一耦合过程，建立相应的数学模型，并进行井筒—随机裂隙网络储层耦合数值模拟，可以预测EGS 运行状态、利用效率与使用寿命，对系统进行采热评价^[6]。对于热储内裂隙岩体热流耦合的数学模型可分为等效连续介质模型与离散裂隙网络模型两大类^[7]，后者不仅与现实更接近，且能更好地模拟水热迁移过程，因而受到广泛的关注^[8-14]。显式模拟裂隙的离散裂隙网络模型虽接近于真实情况但计算量较大^[15]。由于裂隙宽度方向相对于其他两个方向，其尺度小到可以忽略，Juanes 等^[16] 提出在描述裂隙的过程中可以采用低一维的单元来进行。在裂隙网络生成上，大部分研究都依蒙特卡罗方法随机生成^[17]。关于随机裂隙网络储层的热流耦合研究很少考虑到井筒部分的热损失，而在对井筒—热储耦合模拟的研究中，大多是将储层看作是单孔隙等效多孔介质模型。Jiang 等^[18] 建立了井筒与热储循环通道热流耦合瞬态模型，其中储层被视为单一孔隙度的等效多孔介质。Zeng 等^[19] 提取了DP23-1 井的地质资料，研究了干热岩储层的热潜能，其储层也被看作是等效多孔介质。曹文炅等^[20]在局部非热平衡假设的基础上，用理想的垂直井模型模拟了EGS 的采热能力，热储部分同样以多孔介质来代替，忽略了裂隙的存在。长深井筒的热损失及随机裂隙网络储层的热流耦合过程在循环系统的整体采热评价中都占据着主导作用，为此有必要开展二者的热流耦合模拟。由于COMSOL 软件在多物理场耦合与裂隙细微结构、井筒细长结构的多尺度耦合上具有强大优势，故选用其进行二维有限元数值模拟，可形象展现系统渗流传热过程，在此基础上分析热物性参数与几何参数变化对裂隙面温度、采出温度、热开采速率等的影响，得出影响系统产能与寿命的各项因素，以期为相关领域设计与施工提供理论依据。