基于最优化原理的热流体循环数学模型求解方法.pdfVIP

2013年 9月 西安石油大学学报(自然科学版) Sep．2013 第28卷第5期 JournalofXianShiyouUniversity(NaturalScienceEdition) Vol_28No．5 文章编号：1673-064X(2013)05-0075-05 基于最优化原理的热流体循环数学模型求解方法 杨元亮 ，王 辉 ，张 建 ，李清方 ，庞会中 ，王增林。，汪俊义4 (1．中石化胜利油田新春采油厂，山东 东营257000；2．胜利油田胜利勘察设计研究院，山东 东营257026； 3．中石化 胜利油田分公司，山东 东营 257000；4．辽河油田欢喜岭采油厂，辽宁 盘锦 124114) 摘要：井筒热流体循环采油工艺的数学模型是一种具有半隐式边界条件的非线性常微分方程组边 值问题，采用打靶法或共轭函数法求解时会出现数值不稳定和结果失真，基于二分原理的隐式条件 显示化方法计算时又较为费时，应用有限差分方程迭代求解时虽然节省时间但精度较低．为提高运 算速度和改善结果精度，依据方程组的边界条件提出一种以待求初值为设计变量，由结果误差构造 目标函数并求其最小值的直接最优化算法．通过几个典型算例对本文方法和文献已有解法进行 了 考核，分析比较了不同方法的计算精度和CPU运算时间，结果表明本文方法是高效准确的． 关键词：热流体循环；非线性常微分方程组；半隐式边界条件；最优化算法 中图分类号：TE345 文献标识码 ：A 井筒热流体循环数学模型的建立和准确求解对 固有不稳定性导致结果失真．于是在此基础上提出 于稠油开发的方案设计、动态预测和效果评价显得 了基于二分原理的隐式条件显示化求初值的求解方 尤为重要．任瑛 ¨J、赵刚 分别针对开、闭式热流体 案，取得了较好的结果，但其运算耗时较长，不利于 循环建立了相应的数学模型并推导出了相应的解析 参数优化设计计算．张建 采用上游加权格式离散 公式，盖平原 针对空心杆中管掺热液体闭式正、 微分方程并采用高斯消去法进行了求解，计算表明 反2种循环工艺分别给出了温度计算解析式，该类 该法具有较好的收敛性和稳定性．吴晗 则采用下 方法简单准确且运算快捷 ，但只对传热系数和水当 游加权离散方程并采用迭代法进行了求解，该法 由 量为常数的情况有效，仅适合给出井筒温度分布的 于原理简单且收敛速度快，在超深井热流循环工艺 定性规律．热流体循环工艺的数学模型是一种具有 计算中得到较好应用．由于文献 [5-6]的差分方程均 半隐式边界条件的非线性常微分方程组边值问题， 为一阶格式，因此结果精度较低．为保证计算过程的 工程实用上需要对该模型建立更为准确可靠的数值 稳定性同时提高运算速度，改善数值结果的精度，本 解法，以得到井筒温度场的定量分布规律．文世鹏等 文依据微分方程组的边界条件并结合最优化原理， 人_4分析指出，采用多点打靶法求解该模型时需要 提出了一种新的热流体循环数学模型求解方法．通 从2个方 向对方程组积分，有可能出现指数增长太 过几个典型算例对本文方法进行了考核，并在计算 快而造成求解过程不稳定，且当积分区间过长时，还 精度和运算时间上与现有解法进行了比较，验证了 可能超出计算机的容量 ，造成数值解的严重失真，而 本文方法的可行性和有效性． 利用共轭梯度法求解时也会因为原方程组有严重的 收稿 日期：2012—12-08 基金项 目：中石化股份有限公司科技攻关项 目“胜利西部新区污水资源化配套技术研究”；中石化胜利油田重点科技攻 关项 目“胜利油田西部浅层稠油热采地面工程配套技术研究” 作者简介：杨元亮 (1967一)，男 ，高级工程师，博士生，主要从事油气 田开发研究．E—mail：bncyc@163．com 一 76一 西安石油大学学报 (自然科学版)