注蒸汽井筒沿程热物性参数及热损失新算法.pdfVIP

第3l卷 第4期 计 算 物 理 Vo1．31，No．4 2014年 7月 CHINESE J0URNAL OFCOMPUTATIONAL PHYSICS Ju1．，2014 文章编号：1001．246X(2014)04-0449—06 注蒸汽井筒沿程热物性参数及热损失新算法 顾 浩， 程林松 ， 黄世军 (中国石油大学 (北京)石油工程教育部重点实验室，北京 102249) 摘 要：基于干度的定义，对传统蒸汽干度方程进行改进，计算井筒热损失速率时，既按深度分段又按注汽时间分 段，并采用Hasan等提出的两相流漂移流动模型计算井筒压力降．结果表明：改进后的模型用新算法的计算结果与 现场实测值吻合很好，比传统的Beggs-Brill算法更准确．另外，在相同井深处 ，随着注汽时间的延长，蒸汽压力、温度 和干度几乎不变 ，但是热损失速率逐渐降低．新算法为准确计算井筒热损失提供 了参考． 关键词 ：热物性参数 ；热损失；压力降；算法 中图分类号：TE345 文献标志码 ：A O 引言 在稠油热采注蒸汽阶段 ，准确计算井筒沿程热物性参数 (压力 、温度和干度)以及热损失是评价注汽效 果的关键 ，对提高蒸汽的热利用率具有重要指导意义．目前，国内外开展了大量研究 ¨。2]，但是仍存在不少问 题．Ramey 首先建立了沿程温度分布和井筒传热模型，但是当注汽时间小于7d时，无 因次地层导热时间函 数不适用；Chiu等 对地层导热时间函数进行了改进 ，但是在计算沿程压力降时仍采用传统的Beggs—Brill 算法，而该算法以及流动型态划分的标准是基于空气 一水两相，并且在计算任意角度倾斜管的持液率时，都 是在水平管的基础上乘 以倾斜校正系数，用于蒸汽 一水两相垂直管流，误差较大；陈月明 在建立沿程干度 方程时，忽略了摩擦力做功；师耀利等 在计算注过热蒸汽井筒热损失时，虽然对井筒深度进行了分段，但 是忽略了按注汽时间分段计算 ，事实上 ，向井筒注蒸汽过程中，径 向热损失 以及水泥环外缘温度既随深度变 化 ，又随时间变化 ，因此 ，在计算时应对时间和深度都进行分段处理． 本文首先根据能量守恒、质量守恒以及传热学理论建立井筒沿程蒸汽压力、温度、干度和热损失模型．建 立蒸汽干度方程时，抛弃传统的复杂方法 ’ ，而是从干度的定义出发直接求解 ，这种处理方法有两大优 点：①使沿程干度计算更准确 ；②加快了整个模型的求解速度 ，因为传统算法对压力和干度都进行迭代 ，这里 只需对压力迭代 ，而干度是在压力求解后直接按定义计算．其次，在计算沿程压力分布时，舍弃传统 的Beggs． Brill和 Orkiszewski等算法 ，采用 Hasan等 叫提 出的适合于蒸汽 一水两相且垂直 向下流动的漂移流动模 型．最后 ，除了分深度段外还分时间段计算井筒热损失 ，使热损失的计算结果更准确 ，利用现场测试数据验证 了改进后的模型及新算法的可靠性 ，并简要分析注汽时间对蒸汽压力 、温度 、干度和井筒热损失速率 的影响． 1 数学模型 基本假设：①井口注入参数(蒸汽的质量流速、压力和干度)保持不变；②井底使用封隔器，保证蒸汽不窜入 油套环空，油套环空充以低压空气；③油管中心到水泥环外缘的热量传递为一维稳定传热，而水泥环外缘到地 层的热量传递为一维不稳定传热，不考虑沿井深方向的传热；④井内油管、隔热层和套管结构如图1所示． 1．1 沿程压力方程 在建立沿程压力方程时，文[11—12]忽略了蒸汽 一水流动的方 向性，对于向上的两相管流，克服重力增 收稿 日期 ：2013—08—02；修 回日期 ：2013—12—10 基金项 目：国家科技重大专项 (2011ZX05012--004，2011ZX05024-005-006)和辽河油田科技项 目(特超深层高存水油藏提高蒸汽热利用率 技术研究)联合资助 作者简介 ：顾浩 (1989一)，男，硕士，博士生，主要从事热力采油技术、油气渗流理论和油藏效值模拟研究，E-mail：guhaol10110@163．COIn 450 计 算 物 理 第 31卷 加压降