理工类毕业设计论文范例

1、.wd辽宁石油化工大学毕业设计 论文 GraduationProject (Thesis) for Undergraduate of LSHU题题 目目T TI IT TL LE E学 院School专业班级Major&Class姓 名Name指导教师Supervisor20 年 月 日论文独创性声明论文独创性声明本人所呈交的论文，是在指导教师指导下，独立进展研究和开发工作所取得的成果。除文中已特别加以注明引用的内容外，论文中不包含任何其他个人或集体已经发表或撰写过的研究成果。对本文的工作做出重要奉献的个人和集体，均已在文中以明确方式标明并致谢。本声明的法律结果由本人承担。.wd特此声

2、明。论文作者签名： 年 月 日.wd摘 要针对稠油热采过程中注入大量高温蒸汽而造成井筒局部热损失难以控制的问题，从热工角度初步分析了高温蒸汽在井筒局部的传热过程，建设了井筒及周围地层局部的二维、瞬态传热数学模型，根据所求得的温度场以及定量热损失经历公式，探讨了隔热油管视导热系数和蒸汽温度对单位长度井筒热损失的影响。研究说明：井筒单位长度热损失随着隔热油管视导热系数和蒸汽温度的增加而增加，且隔热油管视导热系数对井筒单位长度热损失的影响更大。高温高压蒸汽流经井筒段后进入油层段，蒸汽通过射孔进入油层，并与油层进展传质传热，从而形成油层热影响区域。本文利用有限元分析软件 ANSYS，对辽河某热采井建设

3、了井筒、地层和油层局部的二维、非稳态数学模型，并模拟分析了油层和周围地层温度场，在此根基之上探讨了注汽时间对油层热影响区域的影响。研究说明：油层热影响区域先随着注汽时间的增大而增大，而后 基本不变，且注汽时间在 35 天最正确。关键词：视导热系数；注汽温度；热损失；热影响区域；ANSYS.wdAbstractAbstractFor the purpose of reducing the heat loss in the mineshaft during the steam injection process, a two dimension, unsteady state heat trans

4、fer model of the mineshaft and surrounding soil was established and then the temperature field was obtained. Base the Heat loss experience formula, the apparent heat conductivity and steam temperature to heat loss per length of the mineshaft influence were discussed. The results shown that the heat

5、loss per length of the mineshaft increased with the apparent heat conductivity and steam temperature increase, furthermore, the apparent heat conductivity affected larger than that of steam temperature.High temperature and high pressure steam flows through the wellbore segment backward, such as the

6、reservoir section, the steam through the perforation into the reservoir, and reservoir mass and heat transfer, and thus the formation of a reservoir heat-affected zone is build. Using the finite element analysis software ANSYS, a thermal recovery wellbore, formation and reservoir part of the two non

7、-steady-state mathematical model of Liaohe was build, and simulation Analysis of the temperature field of the reservoir and surrounding strata was made, and the relationship between the time of steam injection and the reservoir heat-affected zone was discussed. The studies have shown that: at first,

8、reservoir heat-affected regionincreased with the increase of the steam injection time, and then kept install, the time of steam injection in 3 to 5 days was the best.KeywordsKeywords: Thermalconductivity; Steamtemperature; Heat loss; Heat affected zone;ANSYS.wd目录1 1 绪论绪论 1 11.1 目的意义 11.2 国内外的研究现状 21

9、.2.1 井筒热损失的研究 21.2.2 油层热影响区域的研究 31.3 本文的研究内容及创新点 41.3.1 研究内容 41.3.2 创新点 52 2 研究方法研究方法 6 62.1 有限元分析法简介 62.2 ANSYS 简介 62.2.1 ANSYS 开展过程 72.2.2 ANSYS 技术特点 72.2.3 ANSYS 使用环境 82.2.4 ANSYS 程序功能 82.2.5 ANSYS 热分析概述 92.2.6 ANSYS 求解步骤 93 3 井筒热损失井筒热损失 11113.1 问题描述 113.2 数学模型的建设 113.2.1 能量平衡方程 113.2.2 连续介质的守恒方程

10、 133.2.3 传热过程分析 14.wd3.2.4 基本假设 163.2.5 控制方程 163.2.6 边界条件 163.3 ANSYS 模拟 163.3.1 建设有限元模型 163.3.2 施加载荷计算 173.3.3 求解 183.3.4 后处理 183.4 模拟结果与分析 183.4.1 隔热油管视导热系数对井筒热损失的影响 183.4.2 注汽温度对单位长度井筒热损失的影响 204 4 油层热影响区域油层热影响区域 22224.1 问题描述 224.2 数学模型的建设 224.2.1 连续介质的守恒方程 234.2.2 多孔介质的守恒方程 254.2.3 传热过程分析 264.2.4

11、 基本假设 274.2.5 边界条件 274.3 ANSYS 模拟 274.3.1 建设有限元模型 274.3.2 施加载荷计算 284.3.3 求解 294.3.4 后处理 29.wd4.4 模拟结果与分析 295 5 结论结论 3232参考文献参考文献 3333谢辞谢辞 3535.wd1 绪论1.1 目的意义稠油在世界石油资源中储量丰富，据统计，其总量约为 1000108吨，中国是世界上稠油资源丰富的国家之一，预测资源量约为 198 亿吨。稠油由于粘度大，流动性差等特点，常规开采方法并不适用，基于稠油对温度的敏感性特点，热力采油已成为国内外普遍认同的开采方法，即向地下油层内注入高温高压蒸汽

12、，一方面可以提高油层温度、降低稠油粘度；另一方面可使油层压力升高、增加驱油能力。常见的热力采油方法有“蒸汽吞吐采油和“蒸汽驱采油两种方式。 “蒸汽吞吐采油是先向油井注入一定量的蒸汽，焖井一段时间，待蒸汽的热能向油层扩散后，再开井生产的一种开采稠油的增产方法。蒸汽吞吐作业的过程分为三个阶段，即注汽、焖井及回采。“蒸汽驱采油是稠油油藏经过蒸汽吞吐采油之后，为进一步提高采收率而采取的一项热采方法，因为蒸汽吞吐采油只能采油各个油井附近油层中的原油，在油井与油井之间还留有大量的死油区。 “蒸汽驱采油是由注入井连续不断地往油层中注入高干度的蒸汽，蒸汽不断地加热油层，从而大大降低了地层原油的粘度。注入的蒸汽

13、在地层中变为热的流体，将原油驱赶到生产井的周围，并被采到地面上来。在注蒸汽的过程中，由于蒸汽与地层之间存在温差，蒸汽流动过程中又有摩擦阻力，因此注汽过程中存在着能量损失，注汽过程中的能量损失，特别是井筒中的能量损失，直接影响热采效果。为了有效地减少这局部能量损失，提高热采效果，需要对注汽井筒进展传热分析，进而探讨隔热油管视导热系数和蒸汽温度对单位长度井筒热损失的影响，这为热采过程中热量的高效、合理利用提供了一定的理论依据。蒸汽流经井筒后进入目的层即油层，高温高压的蒸汽会通过射孔与油层进展传质传热，将热量释放给原油，使其温度升高，粘度降低，同时增加油层压力，到达驱油效果。但随着蒸汽与油层热量交换

14、的进展，蒸汽自身温度降低，逐渐冷凝成液体，阻止热量继续扩散，会形成热影响区域。研究油层热影响区域对合理确定注汽时间和钻井位置具有重要指导意.wd义。1.2 国内外的研究现状1.2.1 井筒热损失的研究自本世纪初以来，就有许多研究者致力于井筒热损失的研究。井筒传热研究工作的先驱 Ramey 以能量平衡为根基，忽略了流体流动过程中动能和摩擦的影响，最先提出了考虑时间和井深影响时的井筒液体温度计算的理论模型。该模型作为井筒热损研究的雏型，给出了井筒内注入气体和流体两种工质关于井深和时间的温度解析表达式，为后来的井筒热损研究奠定了根基。但该模型中的时间函数要求注入流体的时间不得小于 7 天，所以此模型

15、的精度在时间小于 7 天的时候受到限制。另外该模型仅适用于单相流，并在建设模型过程中假设井筒流体的物性参数不随温度变化，这与实际物理过程不相符1。国内外许多研究者在 Ramey 模型的根基上继续做了多方面的研究，如 Hasan 和 Kabir 建设起井筒内两相流模型，改善了 Ramey 模型仅适用于单相流的缺乏，此模型可用于计算井筒液体温度的分布26。Pacheco 和 Farouq Ali 以流体流动的能量平衡方程为根基，建设了能同时计算蒸汽压力和蒸汽干度的数学模型7。该模型弥补了 Ramey 模型中没有考虑水泥壳的影响的缺点。另外此模型构造形式较复杂。但它可以用来较准确地计算蒸汽物性和热损

16、失随深度和时间的变化。不过此模型还用到了 Ramey 模型中不完善的理论，如地层时间函数f ( )的计算，所以适用精度也受到一定限制。Satter 在使用 Ramey 模型时，考虑了注汽过程中的蒸汽相变影响，进而提出了一种预测冷凝流体干度随深度及时间变化的计算方法，模型中提出的关系式8可用于计算饱和蒸汽注入时的井筒热损失。Squier 等人对热水通过井筒注入时提出了一个完整的分析方法9；Huygen 和 Huitt10给出了在蒸汽注入期间的井筒热损失理论和实验结果，同时指出辐射热损失的重要性；Willhite 提出了完整的计算总传热系数的方法11。Holst和 Flock 以及 Earloug

17、her 总结前面的一些成果，并考虑了蒸汽压力的计算12,13；Eickmeier等人采用有限差分模型计算注入初期过渡过程特性14。Hasan 等人又对井筒压力降和热损失计算进展了完善。但这些模型对于目前井筒注汽热损失的计算尚有一些缺乏，首先它们在.wd计算过程中，把某些参数看作是不随温度变化的常量，这与实际并不完全相符；另外由于隔热管的开展，上述模型中的总传热系数不能较真实地反映井筒现在构造，给计算精度的提高带来困难。我国对稠油热采的相关研究比较晚。胡智勉、王弥康对注汽井筒总传热系数的计算方法提出了一些见解15,16。王弥康在 Ramey 模型的根基上，从传热机理出发，给出了蒸汽井井筒热传递定

18、量计算的数学模型，该模型防止了 Ramey 公式仅适用于长时间注汽的缺乏，扩大了计算时间的适用范围，从而改善了井筒传热计算的精度17。沈惠坊将地面管线与井筒进展了综合考虑，探讨了计算蒸汽压力降、干度降及热量损失的计算模型18。李景勤、陈艳华、王志国等人建设了计算井筒热损数学模型，并提出了热损合理分布的概念1921。倪学锋提出了井筒内参数计算模型22。刘文章提出了用物理模拟方法确定井筒总传热系数的方法23。1.2.2 油层热影响区域的研究1968 年 SpilletteNielsen 通过建设一个二维模型，对油层注冷水和注热水后温度分布进展了模拟计算，把流场分布和温度分布藕合在了一起进展迭代24

19、。1970 年埃克米尔等人在研究井筒换热和温度分布时，建设了一套关于注液和产液期间，液体与井筒周围地层间热交换的有限差分模型。计算中将油管、环空液、套管、水泥环及地层情况全部考虑在内。该模型可以有效地模拟实际井筒条件，计算短时间和长时间注液、产液时井筒内温度的变化，模型只是基于热传导的模式进展分析的，没有考虑注入层的对流传热和井筒中的对流传热。1970 年，Smith，R.C.，Steffensen，R.J.等人提出了利用温度测井的计算机分析估计流动剖面25。他们通过改变影啊温度特性的参数，对注水井的温度测井进展了详细的研究，研究结果说明温度测井对地层的热学特性以及井眼的几何形状具有敏感性:另

20、外，指出了注入水的温度和小流量的漏失对温度测井的影响，并对温度测井解释提出了指导性意见。1972 年维脱霍特与梯克萨通过改变注入速率、注入深度和注入时间，来研究注入期间井下的温度分布，并对注入层中的温度分布进展了研究。1973 年 Steffenson.RJ，SmlthR.C.等人提出在油层.wd中，由于流动的存在，势必发生热交换，要计算井眼的温度必须确定井筒周围油藏的大片区域内的温度分布，并提出了著名的 Joul 一 Thomson 效应。1982 年 John Fagley，H.SeottFogler，C.BrentDavenport，R.S.Millhone 等人对关井前后井筒的热传递进

21、展了数值模拟26，提出关井前一段时间的热注是一种解释老井注入剖面的潜在的重要工具。1988 年袁益壮针对在油藏中水驱油的流动和传热问题提出了多孔介质中完全可压缩、可混溶驱动问题的差分方法。1994 年 FredericMaubeuge，MichelP.Didek，MartynB.Beardsen 等人提出了一个多层注入模型27。1995 年胡耀江等人建设了水饱和的一维多孔介质在一端注入蒸汽后，热液区温度分布的数学模型，提出了运用拉氏变换及数值反演求解此问题的计算方法28。1996 年冯恩民等人在建设的温度模型中，除了考虑热传导外，特别强调了水井中与油层内传质传热，构造了关于导热系数、比热等参变

22、量的优化识别模型29。还有局部学者利用解析法求解地层中的温度分布。1.3 本文的研究内容及创新点1.3.1 研究内容本文研究内容主要分为两个局部，即井筒热损失和油层热影响区域。A.井筒热损失利用传热学、热力学、流体力学等知识分析隔热油管内管至地层的传热过程，建设连续介质的守恒方程，并建设数学模型；利用有限元软件 ANSYS 建设隔热油管内管至地层的有限元模型，并利用其热分析模块对模型进展热分析，求得注汽过程中井筒的温度场分布云图；探讨注汽温度和油管视导热系数对井筒热损失的影响。B.油层热影响区域利用传热学、热力学、流体力学等知识分析套管至油层的传热过程，建设连续介质的守恒方程和多孔介质的守恒方

23、恒，并建设数学模型；利用有限元软件 ANSYS 建设套管至油层的有限元模型，并利用其热分析模块对模型.wd进展热分析，求得不同时刻油层的温度场分布云图；探讨注汽时间与油层热影响区域之间的关系。1.3.2 创新点1在以往的井筒温度场及热损失研究中，以数值解法为主，本文利用有限元软件ANSYS 的热分析模块求得了井筒温度场分布云图，并结合热损失计算经历公式，探讨了注汽温度和油管视导热系数对热损失的影响。与数值法相比，本文研究方法可以清楚、直观地观察到温度场的分布情况。2利用有限元软件 ANSYS 的热分析模块模拟分析了油层热影响区域，并探讨了注汽时间与油层热影响区域之间的关系。.wd2 研究方法2

24、.1 有限元分析法简介有限元法finite element method是一种高效能、常用的计算方法。有限元法在早期是以变分原理为根基开展起来的，所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中这类场与泛函的极值问题有着严密的联系 。自从 1969 年以来，某些学者在流体力学中应用加权余数法中的迦辽金法Galerkin或最小二乘法等同样获得了有限元方程，因而有限元法可应用于以任何微分方程所描述的各类物理场中，而不再要求这类物理场和泛函的极值问题有所联系。其原理是将连续的求解域离散为一组单元的组合体，用在每个单元内假设的近似函数来分片的表示求解域上待求的未知场函数，近似函数通常由未

25、知场函数及其导数在单元各节点的数值插值函数来表达。从而使一个连续的无限自由度问题变成离散的有限自由度问题。有限元法作为目前工程中应用较为广泛的一种数值计算方法，以其独有的计算优势得到了广泛的开展和应用，ANSYS 软件以其多物理场耦合分析功能而成为 CAE工程设计中的计算机辅助工程软件的应用主流，在热分析工程应用中得到了较为广泛的应用。2.2 ANSYS 简介ANSYS 程序是融构造、热、流体、电磁、声学于一体的大型通用有限元商用分析软件，可广泛应用于核工业、铁道、石油化工、航空航天、机械制造、能源、电子、造船、汽车交通、国防军工、土木工程、生物医学、轻工、地矿、水利、日用家电等各个领域。该软

26、件可在大多数计算机及操作系统中运行，从 PC 到工作站直至巨型计算机，ANSYS 文件在其所有的产品系列和工作平台上均兼容；该软件基于 Motif 的菜单系统使用户能够通过对话框、下拉式菜单和子菜单进展数据输入和功能选择，此举大大方便了用户操作。它由世界上著名的有限元分析软件公司美国 ANSYS 公司开发的，能与多数 CAD 软件配合使用，实现数据的共享和交换，是现代产品设计中的高级 CAD 工具之一。在 ANSYS 公司相继收购.wdICEM、CENTURY DYNAMICS、AAVID THERMAL、FLUENT 等世界著名有限元分析程序制造公司并将产品整合之后，ANSYS 实际上已成为

27、世界上最通用和有效的商用有限元软件。2.2.1 ANSYS 开展过程ANSYS 公司是由美国匹兹堡大学的 John Swanson 博士在 1970 年创立的，其总部位于美国宾夕法尼亚州的匹兹堡，目前是世界 CAE 行业最大的公司。自创立以来，ANSYS 公司一直致力于设计分析软件的开发，不断吸取新的计算方法和计算技术，领导着世界有限元技术的开展。ANSYS 程序的最初版本与本文用的版本 ANSYS10.0 相比有很大区别。最初版本仅仅提供了热分析和线性分析功能，是一个批处理程序，而且只能在大型计算机上使用。20 世纪70 年代初，随着非线性、子构造以及更多的单元类型的参加。ANSYS 程序发

28、生了很大的变化，新技术的融入进一步满足了用户的需求。20 世纪 70 年代末。交互方式的参加是该软件最为显著的变化，此举使得模型生成和结果评价大为简化。2.2.2 ANSYS 技术特点与其他有限元计算软件相比，ANSYS 具有以下技术特征： 能实现多场及多场耦合功能。 集前后处理、分析求解及多场分析于一体。 独一无二的优化功能，唯一具有流场优化功能的 CFD 软件。 具有强大的非线性分析功能。 具备快速求解器。 最早采用并行计算技术的 FEA 软件。 支持从微机、工作站到巨型机的所有硬件平台。 从微机、工作站、大型机直至巨型机，与所有硬件平台的全部数据文件兼容。 从微机、工作站、大型机直至巨型

29、机，所有硬件平台具有统一用户界面。.wd 可与大多数的 CAD 软件集成并具备接口。 可进展智能网格划分。 具有多层次多框架的产品系列。 具备良好的用户开发环境。2.2.3 ANSYS 使用环境ANSYS 程序是一个功能强大的有限元计算分析软件包。它可运行于 PC、NT 工作站、UNIX 工作站以及巨型计算机等各类计算机及操作系统中，其数据文件在其所有的产品系列和工作平台上均兼容。该软件有多种不同版本，本文应用版本为 ANSYS10.0 版，其微机版本要求的操作系统为 Windows95/98/2000/XP，也可运行与 UNIX 系统中；微机版的 基本硬件要求为：显示分辨率为 1024768

30、 像素，显示内存为 2MB 以上，硬盘大于 350MB，推荐使用 17 英寸显示器。ANSYS 多物理场耦合的功能，允许在同一模型上进展各种耦合计算，例如：热-构造耦合、热-电耦合、磁-构造耦合以及热-电-磁-流体耦合，同时在 PC 上生成的模型可运行于工作站及巨型计算机上，所有这一切就保证了 ANSYS 用户对多领域工程问题的求解。ANSYS 可与多种先进的 CAD如 AutoCAD、Pro/Engineer、NASTRAN、Alogor、I-DEAS 等软件共享数据，利用 ANSYS 的数据接口，可以准确地将在 CAD 系统下生成的几何模型数据传输到 ANSYS，并通过必要的修补可准确地在

31、模型上划分网格并进展求解，这样就可以节省用户在创立模型的过程中所花费的大量时间，使用户的工作效率大幅度提高。2.2.4 ANSYS 程序功能ANSYS 程序主要包括 3 个局部30：前处理模块、求解模块和后处理模块。前处理模块提供了一个强大的实体建模及网格划分工具，用户可以方便地构造有限元模型；求解模块包括构造分析构造线性分析、构造非线性分析和构造高度非线性分析 、热分析、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理.wd介质的相互作用，具有灵敏度分析及优化分析能力；后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流显示、立体切片显示、透

32、明及半透明显示等图形方式显示出来，也可将计算结果以图表、曲线形式显示或输出。ANSYS 程序提供了 100 种以上的单元类型，用来模拟实际工程中的各种构造和材料。2.2.5 ANSYS 热分析概述在实际生产过程中，常常会遇到多种多样的热量传递问题：如计算某个系统或部件的温度分布、热量的获取或损失、热梯度、热流密度、热应力、相变等。所涉及的部门包括：能源、化工、冶金、建筑、电子、航空航天、农业、制冷、船舶等。以机械加工为例，往往需要估算和控制工件温度场；分析不同工作条件下，不同材料及几何形状对温度场变化的影响以及防止加工过程工件中缺陷的产生。ANSYS 作为新颖的有限元分析软件在热分析问题方面具

33、有强大的功能，而且界面友好，易于掌握。用户可以随心所欲地选择图形用户界面方式GUI或命令流方式进展计算。ANSYS 的热分析模块一般包括ANSYS/Muliphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED 五种，其中 ANSYS/FLOTRAN 不含相变热分析。ANSYS 热分析基于能量守恒原理的热平衡方程，用有限元法计算物体内部各节点的温度，并导出其他热物理参数。运用 ANSYS软件可进展热传导、热对流、热辐射、相变、热应力及接触热阻等问题的分析求解。2.2.6 ANSYS 求解步骤ANSYS 求解模型的过程中的主要步

34、骤：定义单元类型及选项，根据所分析的问题在单元库中选择分析问题的单元类型；定义单元是常数，单元是常数是与特定单元有关的量；定义材料属性，在材料属栏中选择性的输入热导率、泊松比、密度等所需参数；建设几何模型；定义网格并划分参数；.wd网格划分以创立的模型。ANSYS 在后处理器中可以查看结果：显示等值线；显示分析结果数据；对结果数据进展计算和路径操作；绘制云图。结合本文研究内容及 ANSYS10.0 软件的 ANSYS/Thermal 热分析模块功能，ANSYS10.0适用于本文研究内容。.wd3 井筒热损失3.1 问题描述以辽河油田某钻井为研究对象，注汽时间T=5d，油套环空导热系数yt为0.

35、023W/mK，套管导热系数tg为 43.2 W/mK，水泥环导热系数cem为 0.81 W/mK，地层导热系数cm为 1.745 W/mK，地层热扩散系数c为 0.027 m2/h，地层初始温度Tcm=40oC，蒸汽与隔热油管内壁对流换热系数为 10 kw/m2oC。3.2 数学模型的建设隔热油管内管内半径r1为 0.031m，外半径r2为 0.0365m；隔热油管外管内半径r3为0.05015m，外半径r4为 0.05715m；套管内半径r5为 0.0809m，外半径r6为 0.0889m；水泥环外半径r7为 0.1239m；地层半径r8为 0.2478m，详见图 3-1 井筒及周围地层截

36、面示意图。图图 3-13-1 井筒及其周围地层截面示意图井筒及其周围地层截面示意图3.2.1 能量平衡方程依据热平衡原理，选取井筒内的一微元流体单元体为研究对象，如图 3-2 所示。r8r7r6r5r2r1r4r3蒸 汽隔热层环 空 水泥环地层 .wd图图 3-23-2 流体单元能量平衡模型流体单元能量平衡模型当蒸汽纵向流入或流出该单元体的控制面 1、2 时，携带的能量可表示为：3-12 1i ).2(2，iiioigzvhmE式中：蒸汽携带的能量，W；iE蒸汽流量，kg/s；om蒸汽焓值，J/kg ；ih蒸汽流速，m/s ；iv 高度，m ；iz重力加速度，9.81m/s2。g该流体单元的能

37、量平衡表达式为：3-2)2()2(1211022220zgvhmzgvhmQ式中：单位时间的径向散热量，W。Q蒸汽的物性参数在井筒截面上是不变的，但不同截面处的物性参数是随温度变化的。E1E2123Q.wd由式3-2看出蒸汽沿井筒纵向的焓值在减少，蒸汽的能量是按递减规律传至地层的。3.2.2 连续介质的守恒方程A.井筒内流体的能量平衡方程井筒油管内高温蒸汽在注入过程中遵循热力学第一规律，其能量平衡方程为：3-3dzdzzdzzzzzQvgHvgzH2dz)-(z22-dz-z-2式中：、z及z-dz处的焓，kJ；zHdzzH-vz、 vz-dz z及z-dz处的速度，m/h；、z及z-dz处的

38、流体密度，kg/m3 ；zdz-zQdz蒸汽散热量，kJ 。公式3-3表示包括蒸汽焓、势能和动能在z处的蒸汽能量与z + dz处的蒸汽能量及蒸汽散热量的平衡方程。B.井筒内流体的连续性方程高温蒸汽在油管内流动时，质量流量保持不变，因此满足流体的连续性方程，即 3-40mvAvAdzzdzzdzzzzz式中：、z及z+dz处的横截面积，m2；zAdzzAvz、 vz+dz z及z+dz处的速度，m/h；、z及z+dz处的流体密度，kg/m3 。zdzzC.井筒内流体的动量平衡方程把蒸汽在油管内的流动看作是两相流的流动，为描述其流动特性，通常从两相流压力降的研究出发。管内稳态均质两相流满足动量平衡

39、方程，即3-5vdvdzdvfgdzdPmmm-2-12式中：段内的压力降，MPa；dPdz蒸汽密度 kg/m3；m.wd摩阻系数；f井筒内径，m。1d表示段内的势能变化， 表示段内的摩擦损失，表示段内gdzmdzdzvdvmdz的动能变化。也就是说管内两相流的压力降是摩擦损失、势能变化和动能变化的综合结果。D.油管内湿蒸汽的焓由于湿蒸汽为汽水两相混合物，所以它的焓可用一样温度或压力下的饱和水、饱和蒸汽及湿蒸汽的干度来表示，即3-6)-1 (xhxhhlg式中：蒸汽干度；x湿蒸汽的焓，kJ/kg；h饱和蒸汽的焓，kJ/kg；gh饱和水的焓，kJ/kg。lh3.2.3 传热过程分析根据热力学31

40、和传热学32知识对井筒传热进展分析，由于水泥环外侧的地层体积和热容量比水泥环内侧的井筒大得多，水泥环与地层交界处的温度 随注入时间的推移缓慢上kt升。表征热量在地层中扩散速率的参数热扩散系数 的数量级仅为 10-7 2m /s ，因而与未受热影响的远处地层温度的差值不大；相反，隔热管内蒸汽温度与的差异较ktcitftkt大，有一个数量级的差异。因而，可以假设水泥环内侧简称井筒内热量传递是稳态的，而水泥层外侧简称地层内热量传递那么按非稳态处理。A.井筒内稳态传热井筒内传热过程主要包括：蒸汽与隔热油管内壁之间的对流换热、隔热油管多孔隔热材料中伴有热辐射的导热、油套管环空中对流、辐射同时作用的换热以

41、及套管、水泥环和地层的导热等环节。dzdvfm122.wd对于注汽井一微元段dz 的径向热损失，可表示为： 3-7)-(22kfttdzkrdQ式中：段内井筒的径向热损失，W；dQdz隔热油管外半径，m；2r以隔热油管内管外外表为基准面积的总传热系数，W/m2k 。k当井筒内隔热管抽真空时， 可表示为：k3-817654321111RRRRRRRk式中：蒸汽凝结换热热阻；1R隔热油管热阻；2R隔热管内隔热层换热热阻；3R环空辐射换热热阻；4R环空对流热热阻；5R套管壁导热热阻；6R水泥环导热热阻。7RB.地层内非稳态传热在井筒稳态传热的同一微元段dz 内，地层非稳态导热方程如下：3-9eeet

42、artrrt1122初始条件时：03-10ZatttTbcie边界条件时：03-110rte3-12772-rrrtdzrdQee式中：地层内某处温度， oC；et， )(rfte.wd地层导热系数，W/mK ；e地层扩散系数，。/smc/2p，ee3.2.4 基本假设1从油管内的蒸汽到地层的热量传递为一维稳态传热，地层为非稳态导热；2视油管内的蒸汽流为一维均质两相流，且无相变及泄露现象；3忽略地层温度梯度的影响。3.2.5 控制方程根据传热 基本理论，二维瞬态传热问题的控制方程为：3-13TyTxT12222式中：x，y分别指的是x和y方向；是时间，s；是热扩散率，m2/s。3.2.6 边界

43、条件蒸汽与油管内壁之间属于对流换热，所以油管内壁为第三类边界，应满足：3-14）（fTTanTk在不考虑地温梯度的情况下，地层温度为定值，因此地层边界处为第一类边界条件，其他边界均为绝热边界。3.3 ANSYS 模拟利用 ANSYS10.0 热分析软件，做出在不同隔热油管视导热系数和注汽温度时的井筒温度场分布情况。3.3.1 建设有限元模型1翻开 ANSYS 10.0 有限元分析软件，过滤图形用户界面；2确定工作文件名“heat loss in well，单位制默认为国际单位制；3进入 PREP7 前处理器，选择 plane55 热分析单元，并设置成轴对称单元；4将物性参数定义为材料热物性参数

44、；.wd5根据图 3-1 建设有限元模型，如图 3-3；6采用自由网格划分方法对模型进展网格划分，结果如图 3-4。3.3.2 施加载荷计算A.定义分析类型定义新的分析类型为瞬态分析“Transient。B.设置瞬态热分析的初始条件由于初始温度场是不均匀的且又是未知的，所以首先要进展稳态热分析，为稳态热分析建设初始条件：在井筒壁处施加对流载荷，在地层处施加原始地层温度场载荷，其他边界为绝热边界。图图 3-33-3 井筒有限元模型井筒有限元模型图图 3-43-4 网格划分结果网格划分结果.wd3.3.3 求解采用 GUI 命令：MainSolutionLoad Step OptsTime Fre

45、quencTime-Time Step，设置计算终止时间为 0.01s，时间步长为 0.01s，并将时间积分效果设定为 Off，对模型进展稳态热分析。稳态热分析完毕后，将时间积分效果设定为 On，并删除所有节点的温度载荷，计算终止时间设定为 5 天，时间步长设定为 1 小时，对模型进展瞬态热分析。3.3.4 后处理对于瞬态热分析，ANSYS 提供了两种后处理方式，即 POST 1 和 POST 26。本文采用POST 1 进展后处理，利用 GUI 命令：Main MenuGeneral PostprocPlot Results 显示温度场彩色云图，并拾取套管内径边缘处温度。 3.4 模拟结果与

46、分析取水泥环与地层之间截面的平均温度作为温度参数，并结合文献33中给出的热损失经历公式，分别算出在不同的隔热油管视导热系数及注汽温度下井筒热损失情况。3.4.1 隔热油管视导热系数对井筒热损失的影响现在使用的隔热油管多为抽真空隔热油管，结合预应力隔热油管标准SY/T 5324-94中提出的视导热系数概念，把隔热管的隔热层内以导热、对流和辐射三种方式传递的热量视为以纯导热方式传递的热量。其中，视导热系数数值划分为五个范围等级，如表 3-1 所示，因此，分别选取视导热系数为 0.07、0.05、0.03、0.01 和 0.005 W/mK。 表表 3-13-1 隔热油管等级分类隔热油管等级分类 单

47、位：W/mK隔热性能等级ABCDE视导热系数0.060.080.040.060.020.040.0060.020.0020.006在实际的稠油热采过程中，注汽温度一般为 300350之间，基于此，本文以蒸汽温度为 312情况下进展模拟，图 3-5 为在视导热系数为 0.07 时的温度分布云图，结合热损失经历公式，得到单位长度井筒热损失与视导热系数之间的关系，如表 3-2 和图 3-6 所示。.wd图图 3-53-5 视导热系数为视导热系数为 0.070.07 W/mW/m K K 时的温度分布云图时的温度分布云图表表 3-23-2 不同视导热系数下的计算结果不同视导热系数下的计算结果序号视导热

48、系数W/mK套管温度oC热损失W/m123450.0050.010.030.050.0741.58842.84846.04647.79748.9022.7154.87010.33813.33215.222.wd图图 3-63-6 隔热油管视导热系数对井筒热损失的影响隔热油管视导热系数对井筒热损失的影响根据图 3-6 分析可知，井筒单位长度热损失整体上随着隔热油管视导热系数的增大而增大。其中，对于 A 级的隔热油管，单位长度井筒热损失约达 15W/m 左右，而对于 E 级的隔热油管，单位长度井筒热损失仅约为 3.5W/m，两者单位长度井筒热损失相差约11.5W/m，因此，选用视导热系数高级别的隔

49、热油管对于降低井筒单位长度的热损失具有明显的影响。此外，当视导热系数从 0.07W/mK 提高至 0.03 W/mK 后，单位长度热损失的变化速率随着视导热系数的减小而出现明显的降低。综上所述，在实际的应用过程中，建议选用隔热油管视导热系数在 0.030.005W/mK。3.4.2 注汽温度对单位长度井筒热损失的影响上述视导热系数的建议范围内，在视导热系数为 0.03 W/mK 情况下，分别研究以注汽温度分别为 312、324、333、346 和 350oC 五种情况下，注汽温度对单位长度井筒热损失的影响，详细情况如表 3-3 和图 3-7 所示，图 3-8 为当注汽温度为 324oC 时的温

50、度分布云图。表表 3-33-3 不同注汽温度下的计算结果不同注汽温度下的计算结果序号注汽温度oC套管温度oC热损失W/m1231232448.90249.29715.22215.897.wd34533334635049.59249.98750.08216.40117.07717.239根据上图分析所知，单位长度井筒的热损失随着注汽温度的升高几乎呈线性增大的趋势。在注汽温度为 312oC 情况下，单位长度井筒热损失约为 15.2W/m，在注汽温度在350oC 情况下，单位长度井筒热损失约为 17.2W/m，由此可见，注汽温度的提高对于井筒单位长度热损失的影响较小，每提高 50oC，井筒单位长度热

51、损失相差约 2W/m。图图 3-73-7 注汽温度对井筒热损失的影响注汽温度对井筒热损失的影响.wd图图 3-83-8 注汽温度为注汽温度为 324324o oC C 时的温度分布云图时的温度分布云图.wd4 油层热影响区域4.1 问题描述高温高压的蒸汽流经井筒后进入封隔器以下的油层区，通过射孔与油层进展传质传热，传热过程如图 4-1 所示。图图 4-14-1 油层传热过程示意图油层传热过程示意图以辽河油田某钻井为研究对象，注汽时间T=10d，套管导热系数tg为 43.2 W/mK，水泥环导热系数cem为 1.2 W/mK，油层综合导热系数cm为 2.53 W/mK，地层导热系数为1.968

52、W/mK，地层热扩散系数c为 0.027 m2/h，地层初始温度Tcm=40oC，蒸汽与隔热油管内壁对流换热系数为 10 kw/m2oC。4.2 数学模型的建设套管内半径r1为 0.0809m，外半径r2为 0.0889m；水泥环外半径r3为 0.1239m；地层和.wd油层半径r4为 20m，详见图 4-2 油层截面示意图。蒸 汽水 泥 环地 层油 层套 管r1r2r30.5mm1m图图 4-24-2 油层截面示意图油层截面示意图4.2.1 连续介质的守恒方程求解与流体中动量、热量和质量传递有关的问题时，首先应根据 基本的守恒定律对该物理系统建设一组方程。这些方程包括:能量平衡方程；连续性方

53、程；焓方程。A.连续介质的能量平衡方程井筒油管内高温蒸汽在注入过程中遵循热力学第一规律，其能量平衡方程为：4-1dzdzzdzzzzzQvgHvgzH2dz)-(z22-dz-z-2式中：Hz、Hz-dzz及z-dz处的焓，kJ；vz、 vz-dz z及z-dz处的速度，m/h；、z及z-dz处的流体密度，kg/m3；zdz-zQdz蒸汽散热量，kJ 。上式表示包括蒸汽焓、势能和动能在z处的蒸汽能量与z + dz处的蒸汽能量及蒸汽散热量的平衡方程。B.连续介质的连续性方程高温蒸汽在油管内流动时，质量流量保持不变，因此满足流体的连续性方程，即 4-20mvAvAdzzdzzdzzzzz式中：、z

54、及z+dz处的横截面积，m2；zAdzzAvz、 vz+dz z及z+dz处的速度，m/h；、z及z+dz处的流体密度，kg/m3。zdzz.wdC.连续介质的动量平衡方程把蒸汽在油管内的流动看作是两相流的流动，为描述其流动特性，通常从两相流压力降的研究出发。管内稳态均质两相流满足动量平衡方程，即4-3vdvdzdvfgdzdPmmm-2-12式中：段内的压力降，MPa；dPdz蒸汽密度 kg/m3；m摩阻系数；f井筒内径，m。1d表示段内的势能变化， 表示段内的摩擦损失，表示段内gdzmdzdzvdvmdz的动能变化。D.油管内湿蒸汽的焓由于湿蒸汽为汽水两相混合物，所以它的焓可用一样温度或压

55、力下的饱和水、饱和蒸汽及湿蒸汽的干度来表示，即 4-4)-1 (xhxhhlg式中：蒸汽干度；x湿蒸汽的焓，kJ/kg；h饱和蒸汽的焓，kJ/kg；gh饱和水的焓，kJ/kg。lh4.2.2 多孔介质的守恒方程A.多孔介质的质量守恒方程多孔介质的质量守恒方程为：4-501vt式中：密度，kg/m3；dzdvfm122.wd介质的孔隙度；渗流速度，m/d。vB.多孔介质的动力学方程a.单相流在只有饱和牛顿流体的空隙介质中，达西定律如下：4-6)-(-Dgpkv式中：流体的渗流速度，m/d；v渗透率，；k2m粘度，Pas；流体的密度，kg/m3；压力，MPa；pD垂直高度，m；G重力加速度，m/s

56、2。b.多相流多相流的渗流微分方程为： (4-7)()-(11llllrlstqDgpkk式中：sll相流体的饱和度；ql流体流入或流出的流体质量流量，kg/s。4.2.3 传热过程分析A.井筒内稳态传热井筒内传热过程主要包括：蒸汽与套管内壁之间的对流换热、套管与水泥环和水泥环与油层的导热等环节。对于注汽井一微元段dz 的径向热流量，可表示为： .wd 4-8)-(22kfttdzkrdQ式中：段内井筒的径向热流量，W；dQdz套管外半径，m；2r以套管外外表为基准面积的总传热系数，W/m2k 。k的表达式为：k4-91321RRRk式中：蒸汽与套管内壁对流换热热阻；1R套管壁导热热阻；2R水

57、泥环导热热阻。3RB.油层传热分析 井筒蒸汽与套管通过对流换热向井筒周围散热，同时产生从射孔到油层的传质热传导。由于在有限元中无法解决传质问题，因此一般简化为具有一定半径的加热区模型。非稳态无内热源轴对称温度场的微分方程为：4-100c-),(p2222tTrTrTrxTrktrxTD应用 Galerkin 法可得到轴对称温度场整体区域“变分计算的 基本方程为: (4-11)TlDlplllPdSnTrkWdxdrtTrWcrTrWxTxWkrTJ0-利用温度场变分方程对模型进展剖别离散、单元分析及最终总体合成，即可完成模型温度场的求解。4.2.4 基本假设沿井深方向上的热能不发生传递；视蒸汽

58、流为一维均质两相流，且无相变；忽略油层温度梯度的影响；以井眼中心为中心的同心圆上的温度是一样的。.wd4.2.5 边界条件蒸汽与套管内壁之间属于对流换热，所以套管内壁为第三类边界，应满足：4-12）（fTTanTk加热区为第一类边界条件，由于不受油层温度梯度的影响，所以油层为第一类边界条件，其他为绝热边界。4.3 ANSYS 模拟4.3.1 建设有限元模型翻开 ANSYS 10.0 有限元分析软件，过滤图形用户界面；确定工作文件名“you ceng re ying xiang qu yu，单位制默认为国际单位制；进入 PREP7 前处理器，选择 plane55 热分析单元，并设置成轴对称单元；

59、将物性参数定义为材料热物性参数；根据问题描述建设有限元模型，结果如图 4-3 所示；图图 4-34-3 油层有限元模型油层有限元模型局部放大图局部放大图.wd采用自由网格划分方法对模型进展网格划分，结果如图 4-4 所示。图图 4-44-4 网格划分结果网格划分结果局部放大图局部放大图4.3.2 施加载荷计算A.定义分析类型定义新的分析类型为瞬态分析“transient。B.设置瞬态热分析的初始条件由于初始温度场是不均匀的且又是未知的，所以首先要进展稳态热分析，为稳态热分析建设初始条件：在套管内壁处施加对流载荷，在油层和地层处施加原始温度场载荷，其他边界为绝热边界。4.3.3 求解采用 GUI

60、 命令：MainSolutionLoad Step OptsTime FrequencTime-Time Step，设置计算终止时间为 0.01s，时间步长为 0.01s，并将时间积分效果设定为 Off，对模型进展稳态热分析。稳态热分析完毕后，将时间积分效果设定为 On，并删除所有节点的温度载荷，计算终止时间设定为 10 天，时间步长设定为 1 小时，对模型进展瞬态热分析。.wd4.3.4 后处理采用 POST 1 进展后处理，利用 GUI 命令：Main MenuGeneral PostprocPlot Results显示温度场彩色云图。4.4 模拟结果与分析模拟油层温度场分布云图如图 4-5 至图