地热井下换热器一种简化模型研究.docVIP

地热井下换热器一种简化模型研究 摘要：针对较高地层温度梯度情况建立了地热井下换热 器的一维非稳态导热数值模拟演化方程，并采用格子 -Boltzmann方法进行了数值求解，模型中考虑了地层温度分 布、井管材料热物性、进水流量以及进水温度对出水温度和 取热量的影响，并通过地热井下换热器的实验数据进行了理 论与实验对比分析。 关键词：地热井下换热器；格子-波尔兹曼方法；导热 模型 分类号：G353文献标志码：A 引百 我国以地源热泵系统实现建筑空调的建筑面积迗到了 一亿平米以上，合理选择地下换热器的形式与设计对减少用 户投资至关重要。地热井下换热器通过U形管或同轴套管实 现管内水与地下土壤或含水层的热量交换，成井后无需回 填，一般应用于有浅层地热条件的供热工况，而地埋管换热 器较适宜于既冬季供热又夏季制冷，有关这两种浅层地热能 利用模式的热力学区别可参见文献［1］。对于井下换热装置 的机理研究基本上采用的都是有恒定热源的非稳态柱坐标 导热问题，除数值方法求解之外［2］，一般解析方法求解主 要分为线热源法［3］和柱热源法［4］。 本文提出的地热井下换热器一般安装在有较高温度梯 度的含水层，换热器周围地层温度将不是恒定温度，采用恒 定热源法将不再适合地热井下换热器的传热性能分析。本文 利用恒定温度圆柱在无限大空间内的非稳态导热公式，考虑 管内对流传热、内热源、管壁热阻等因素建立了一个相对简 单的一维地埋管非稳态导热模型，通过格子-Boltzmann方法 进行了数值求解，并与地热井下换热器实验数据进行了比 较。 地热井下换热器LBM模型的建立［收稿日期： 通讯作者：焦文静(1988-)，女，硕士研究生，主要从 事暖通、热工设计研究HYPERLINK mailto:apple335840@163.comapple335840@163.com; 地热井下换热器的传热性能受地下井成井结构、井周 含水层的水文参数、热物性参数、以及地上设备运行参数的 影响。本文通过建立一维导热LBM模型，可以实现对地热井 下换热器的一维导热数值计算，由于是一维模拟研究，对井 下换热器的几何结构、物性参数等进行了简化处理，主要假 设如下： 井下换热器为U型管，且两管之间传热忽略不计； 井下换热器周围地层为均匀介质； 忽略井下换热器周围地层可能的自然对流，只考虑 导热; 井壁径向瞬态热流密度只与时间项有关，而与温度 的变化过程无关。 假设井管内水为不可压缩、比热为常数，那么带有体积 内热源一维管内流动换热的能量守恒微分方程式如下：（1） 给定无量纲温度、速度、几何尺寸、时间的无量纲参数 及Re数和Pr数定义如下： ，，，， 式中，Ti为入口水温度，Tmax为受干扰热水层最高温 度，U为管内水平均流速，D为U型管内径，a为水的热扩散 系数。对于地热井下换热器，给定的是线热源的热流密度而 不是体积热流密度。由，方程式（1）可以变为：（2） 式中，，并给定ti* = 1，上式即为有线热源无量纲一维 能量守恒微分方程。 格子-Boltzmann方法求解 格子-Boltzmann法（LBM）是一种应用非连续介质（粒 子碰撞）思想研究宏观物理现象，以Boltzmann方程为基本 动力学方程，求解数值流体力学问题的新方法。近年来在模 拟复杂流体、相变传热、多孔介质中都有广泛的应用，具有 适宜于快速并行计算、简便易行等优点［5］。为了利用LBM 方法模拟，将U型的进出水管看成一根直管，相同深度下进 出水管所处的无限远处热储层温度相同。 为了适应在LBM方法计算，取，演化方程式（2）的格式 可写成: 式中，速度U是LBM方法尺度下的给定流体速度，一般 要求U0. 1。，是D1Q3模型下三方向平衡态粒子分布加权系 数。上述模拟中是个随时间、管内水与远端热储层温度相关 的量。为了近似地确定，采用等温圆环初始条件下无限大恒 温空间内的非稳态导热问题的解，其解的积分格式如下：（4） 式中，p为管外流体或多孔介质地层的导热系数，q为 单位面积下的热流密度，（x）为无穷远端热储温度，（x）为管 外壁表面温度，A令的大小仅与有关。 另外通过传热原理有：（5） 式中，为管内对流热阻与管壁导热热阻的和，为管内流 体平均温度。联立式（4）和式（5），可消去，同时考虑进出水 管的相互影响，最后可得：（6） 式中，是远处随地层深度x的温度分布，是井管内水的 平均温度，地热井下换热器的入口边界条件为已知测试温 度，出口为绝热边界。 计算结果及与实验结果的对比分析 为了验证计算程序的正确性，对有内热源的稳态微分方 程式（2）进行了 LBM求解。给定入口与出口的温度，流体的 Re、Pr数及可以通过实验数据计算得到，方程式（2）的定常 稳态解是二阶常微分方程，假设，，则的一般形式解为：（7） 式中，A、B是积分常数