漳州地热系统基岩裂隙热水温度场数值模拟.pdfVIP

·230· 中国地球物理2011 (6)地热资源及其开发利用 漳州地热系统基岩裂隙热水温度场数值模拟 of field Numericalsimulation temperature forfissured at GeothermalField groundwaterZhangzhou 曾玉超苏正 吴能友+王晓星 YuchaoSu Wu Zeng ZhengNengyou+WangXiaoxing 中国科学院可再生能源与天然气水合物重点实验室 广州能源研究所广州510640 目前在东南沿海发现的温度最高的热田，具有极大的研究与开发价值。本文根据漳州热田近年获得的勘测 结果，建立了漳州基岩裂隙热水系统的概念模型，并采用Fluent求解了物理模型，旨在为热田的进一步勘 探、开发和利用提供基础数据和科学依据。 1．研究目的 确定热田深部基岩裂隙热水温度场分布，揭示深部断裂区对热田水流与热流的控制作用，为未来进一 步地热勘探和开发奠定基础。 2．概念模型与算法 漳州地热田内断层分布较为复杂，断层控制着水流与热流的分布。洋西一圆山断裂F。和郭坑一九湖断 裂F3深切入莫霍面，为阻水断裂；而天宝一漳州断裂F2(宽约800m)内岩石破碎，富水性好，为良好的 储水、导水和导热断裂。漳州地热田的水流和温度分布主要受到F2、F3的控制，其它次级断裂对热水出露 有一定的控制作用，但其作用可以忽略。 对漳州热田基岩裂隙系统简化如下：孔隙含水层对基岩裂隙内水温分布无影响，且不补给基岩裂隙水； 地下水从西北、东北和西南方向进入热源层后垂直升流进入模拟区块，在上部沿着断裂带区域补给孔隙热 水；基岩裂隙水为液态；渗透率和孔隙度分布是均匀的；基岩和热水始终处于热力学平衡状态；断裂区形 状是长方体，边界平整。 采用Fluent软件进行数值计算，其多孔介质求解器模块给定了单相液态水渗流与传热的输运方程，状 态方程采用IAPWS—IF97标准计算。模拟区域的温度和压力变化范围较大，计算中考虑温度和压力对水密 度的影响，采用可压缩流动模型计算。整个模拟系统共4个控制方程，4个求解变量：对形成的矩阵方程 Gradient 组采用预处理的双共轭梯度算法(Preconditioned Solver)求解。各个方程计算收敛的 Bi．eonjugate 相对误差限值均为1．0×10～，绝对误差限值为1．0。 3．计算结果 在2000m深度内基岩温度场结构较简单，随着深度增大，基岩温度逐渐升高。对应于地面高温中心的 是西北断裂带中的一个高温热柱，它是地下热水在断裂带中强烈主动对流的结果。此外，在(1000m，500m) 处有另一埋深约400m的高温热柱，其向周围基岩传导热能但是未能刺穿基岩顶部形成高温中心。在2000m 深度下基岩温度场变得极为复杂，因为随着温度升高、压力增大水流的密度和粘度也发生强烈的改变，在 深部是主动对流引起的自然对流及热传导同时存在的区域。 模拟结果表明地热高温中心沿着西北一东南断裂分布，而在东北一西南方向由于断裂区渗透率较小没 着两端热流迅速降低：在东北一西南断裂带上X=1080．25m、Y--1200m处热流最大值为368．72mW／m2，向 断裂两端延伸热流值也迅速降低。模拟结果和热流的测量值分布趋势基本吻合。测量热流的最大值比计算 最大值偏小，同时部分区段模拟值也较测量值偏大，这是因为实际热流值的测量是在地表进行的，地表热 流受到第四系含水层的影响使得基岩顶部裂隙水的热流分布较基岩顶部偏小。 4．结论 漳州热田基岩裂隙热水的温度场受到断裂区的强烈控制，在西北一东南断裂中水流形成了两个强烈对 流的高温热柱，一侧的热柱刺穿了基岩顶部形成热田高温中心。基岩底部热流值在断裂区中部取得最大值， 基本对应地