热补偿作用下增强型地热系统热流动三维动态模拟.docxVIP

热补偿作用下增强型地热系统热流动三维动态模拟 0 热储裂隙岩体的温度演化影响 近年来，随着世界岩浆岩资源总量的加快和开发利用环境保护的加剧，人类必须尽快开发清洁能源。利用3km低渗绿色晶体热涂层（egs）的能耗目标为km以下的低渗干热系统，正逐渐成为世界各国新能源发展的重点方向之一。egs通过手动手段（如水流激励）在干燥的岩石中建立高渗人工热储量，并向低温液体介质中注入可更换能耗。将热流泵输送至地板并在环境中产生热流。研究表明，egs具有良好的研发价值和广阔的应用前景。 理解地下岩体中的热量迁移过程对于高效开采HDR中的热能具有重要意义.通过水力激发等方式建立的人工热储是一个有限的空间, 在采热过程中, 热储周围岩体中几乎不存在流体的流动, 循环工质采集的热量主要来源于热储.随着EGS的运行, 热储与周围岩石之间会出现越来越大的温差, 周围岩石对热储将产生一定的热补偿作用 (thermal compensation) .热补偿的长期积累最终会对采热性能及EGS的运行寿命产生一定影响.在当前的EGS研究中, 涉及到热补偿作用的研究工作极少, 探究热补偿对地下采热过程的细节影响对提高EGS的经济性有一定的参考价值. 由于人工热储内裂隙网络的极端复杂性, 精确构建与实际情形相同的物理模型非常困难.在对EGS的数值模拟工作中, 依据热储内裂隙的不同物理描述, 当前的数值模型可分为三类:单孔隙率模型 (single porosity model), 双孔隙率模型 (dual porosity model)和多孔隙率模型 (multi-porosity model).单孔隙率模型宏观上均化了裂隙和岩体骨架的物性参数, 数学描述时将热储视为均匀性质的多孔介质、计算量较小.但由于传统的单孔隙率多孔介质采用的是局部热平衡 (local thermal equilibrium) 假设, 不能描述岩石骨架和裂隙流体之间的对流换热.双孔隙率模型将热储分为两个均质的多孔介质子区域:低孔隙率的岩石区和高孔隙率的裂隙网络区.虽然这种模型仍采用局部热平衡假设, 但由于其分别处理岩石区温度和裂隙流体区温度, 因此也可近似处理流体和岩石之间的对流换热, 当前应用较为广泛 (如TOUGH2, TOUGHREACT) .但双孔隙率模型中考虑的裂隙往往只能是简单理想的直通道 (组) , 而且裂隙区和岩石区各自的物性参数无法通过地质勘测或实验数据直接得到, 较难于准确设定.多孔隙率模型考虑了实际热储中裂隙的空间分布与几何形状, 通过设定地质物理参数在空间的分布函数来实现更精确的模拟.但这需要大量的实地数据支持, 当前还未见广泛应用. 我们在综合考虑上述裂隙岩体内流体数值模拟方法的基础上, 提出了单孔隙率双能量方程模型.本文运用该模型对不同地质条件下EGS的长期运行过程进行模拟, 分析热储周围岩石的热补偿对EGS热产出以及热储内岩石、流体温度演化的影响, 为EGS的相关应用提供启示. 1 模型建立及求解 EGS的地下部分由注入井、生产井、热储和热储周围的岩体等三个不同的区域组成 (图1) .模型将EGS的地下部分分为三个性质不同的子区域:①开放流道性质的注入井和生产井;②多孔介质性质的热储;③渗透性可忽略不计的热储周围岩石.模型假设单相流体流动, 不考虑循环流体与岩石的化学作用以及流体物性的变化, 同时忽略岩石热应变导致的热储孔隙率以及渗透率的变化.采用的控制方程如下: 连续性方程 动量方程 流体的能量方程 岩石的能量方程 式中ε, ρ, u, P, μ, K, cp和g分别表示热储孔隙率, 密度, 流体表观速度, 压强, 流体动力粘性系数, 热储渗透率, 比热容和重力加速度.h和a分别表示岩石-流体对流换热系数和热储内裂隙的比表面积, 下标s和f分别代表固体岩石和流体.在EGS的实际运行中, 热储中基岩温度往往较高 (423.15 K) , 而注入的循环流体一般为浅层地下水或地面电厂发电后排出的温水或进一步利用后的低温热水 (353.15 K), 采热流体和热储岩石之间存在较大的温差.为充分考虑岩石和流体之间的对流换热, 模型基于局部非热平衡 (local thermal non-equilibrium) 思想. 数值求解采用通用计算流体软件FLUENT6, 通过加载用户自定义函数 (user-defined functions) 来实现各特殊项处理.模型的详细描述及具体求解方法请参考文献. 2 热储系统构成 为探究热储周围岩石的热补偿对EGS采热过程的影响, 考虑如图2所示的某假想双井EGS.人工热储的体积为500 m×500 m×500 m, 在xy平面呈菱形, 热储中心位于地表以下4 000 m深处.整个计算区域大小为6 000 m×2 000 m×2 000