三维非稳态热传导问题的边界元法.pdfVIP

第23 卷 第 18 期 岩石力学与工程学报 23(18) ：3168～3173 2004 年9 月 Chinese Journal of Rock Mechanics and Engineering Sept.，2004 三维非稳态热传导问题的边界元法 刘承论 秦忠诚 ( 山东科技大学资源与环境工程学院 泰安 271019) 摘要 为预测高温岩体热液储层抽取热能的生产能力和寿命，需建立一套数值模拟系统，而人工裂隙与高温岩体 间的三维非稳态热传导边界元数值模拟系统就是其中的一个子系统。为此，研究开发了该子系统，将之与解析解 验证的结果说明，该数值系统可靠、精度高；另外，针对 GAUSS 数值积分法的局限性，开发了两种高精度的特 殊的积分方法。 关键词 数值分析，地热能，热液储层，非稳定热传导，边界元法，高精度积分法 分类号 O 414.14，O 241.83 文献标识码 A 文章编号 1000-6915(2004)18-3168-06 BOUNDARY ELEMENT METHOD FOR 3D NON-STEADY HEAT CONDUCTION Liu Chenglun ，Qin Zhongcheng (College of Resources and Environmental Engineering，Shandong University of Science and Technology ， Tai′an 271019 China) Abstract A set of numerical simulation system should be established in order to predict the output and life of underground geothermal energy generation. One of the numerical subsystems is developed for 3D non-steady heat conduction in this chapter，and is verified by analytical solution. Aimed at the limitation of Gauss integral in this system，two special integral methods with high precision are developed. Key words numerical analysis ，geothermal energy，thermal storage beds ，non-steady heat conduction ，BEM ， integral method with high precision 具有解析解的球模型和由作者推导获得解析解的无 1 绪 言 限域内的空球模型进行了解析和数值模拟，验证了 该数值系统的可靠性和模拟精度；最后，对接近实 本研究课题的最终目的是研究开发地下干燥高 际人工裂隙的圆板状裂隙用TIME INTEGRAL 法进 温岩体人工热液储层抽取热能的生产能力及其寿命 行了数值模拟与分析。 预测数值模拟系统。本文主要基于直接边界元法就 其中人工裂隙与高温岩体间的三维非稳态热传导问 2 三维非稳态热传导边界元系统 题[1，2]开发了TIME INTEGRAL 数值系统，并有针 对性地开发了两种特殊的积分方法；然后对选取的 2.1 积分方程式 2003 年8 月4 日收