基于边界元法含局部不渗透区域任意形状气藏渗流问题.pdfVIP

大 庆 石 油 学 院 学 报 第 33卷 第 2期 2O。¨9年 4月 J()URNAIOFDAQING PETR()LEUM 1NSTITUTE Vo1．33 N0．2 Apr． 2OO9 基于边界元法的含局部不渗透区域 任意形状气藏渗流 问题 王海涛 ，张烈辉 ，冀秀香 ，史文选 ，王 宇。 (1．西南石油大学 油气藏地质及开发工程国家重点实验室 ，四川 成都 6lO5(】()； 2．中国石油吉林油 田公司 勘探开发 研究院，吉林 松原 138O0l； 3．西南油气 田公司 采气工程研究院，四川 广汉 6183OO) 摘 要：为求解任意形状气藏的渗流问题，建立含局部不渗透区域的任意形状气藏的边界元渗流模型，求解模型并 绘制井底压力动态的典型曲线 ；根据渗流机理和渗流过程 ，分析曲线 的形态 ；将简单气藏模型的边界元解与解析解进行 比较 ，验证该气藏边界元模型及方法的正确性和可靠性．结论有助于认识该类复杂气藏 的渗流机理和渗流过程 ，也可用 于试井分析 ，识别大范围内不渗透区域的存在，为认识和开发气藏提供理论基础． 关 键 词 ：气藏 ；不渗透区域 ；任意形状 ；边界元 ；积分方程 ；广义 函数 ；拉氏变换 中图分类号 ：TE31 文献标识码 ：A 文章编号：1O0O一1891(2∞9)O2一O062一O6 0 引言 与解析法相 比，边界元方法可以处理更复杂的油气藏渗流问题 ，如任意不规则形状油气藏 、井位于内 区或外区的任意不规则形状复合油气藏 、一部分边界不渗透而另一部分边界恒压油气藏 、既非封闭也非恒 压油气藏 、任意形状油气藏中的多井问题等．1973年 ，GringartenAc和RameyJrHJⅢ将格林 函数理论 引入油藏不稳定渗流问题的研究 ，为边界元方法在油藏渗流中的应用打下基础；1977年 ，IiggettJA和 LiuP 首次利用边界元方法研究多孔介质中的渗流问题 ；1992年，KikaniJ和 HorneRN[。应用边界 元方法分析任意形状油藏的不稳定压力动态问题 ，给出在真实空间和Iaplace空间求解均质油藏不稳定 渗流问题的边界元公式 ，并通过与简单模型解析解的对 比，说明2种方法具有较高的计算精度 ，可以应用 于存在多井及复杂边界条件影响的油藏压力动态计算 ，但在 Laplace空间处理更方便 ；SatoK和 HorneR Nf 提出将扰动方程加入控制方程中求解 自由空间的格林 函数 ，利用边界元法求解孔隙度和渗透率发 生变化的渗流问题 ；1993年，张望月和 曾萍__6利用边界元方法在 Iaplace空间研究不规则边界均质油藏 中 含不渗透区域的渗流问题 ；2O06年 ，何应付 等利用边界元方法研究双重介质油藏的压力动态．该领域的 研究都局限于油藏渗流问题 ，研究 的流体为不可压缩或微可压缩流体 ，渗流时其黏度不随压力变化 ，而对 于气藏内可压缩流体渗流的研究还未见报道．拟提 出含局部不渗透区域的任意形状气藏内可压缩流体 的 边界元渗流模型，并绘制井底压力动态的典型曲线． 1 气藏渗流边界元数学模型 假设气藏水平等厚、气体可压缩、等温渗流．根据质量守恒原理 、达西运动定律 、气体状态方程 ，结 合点源的思想和函数 的特性 。，推导得到压力形式的气藏边界元渗流控制方程 ： 收稿 日期 ：2)【(】8一O7一l8；审稿人 ：刘义坤 ；编辑 ：任志平 基金项 目：教育部博士点基金项 目(2()O4O615()(】4)；四川省学术和技术带头人培养基金项 目(22oo320)；高等学校优秀青年教师教学科 研奖励计划项 目(TRAP()YT) 作者简介：王海涛(】98o一)，男，博士生，主要从事油气田开发、试井分析及油气藏数值模拟方面的研究． · 62 · 第 2期 王海涛等 ：基于边界元法 的含局部不渗透区域任意形状气藏渗流问题 ( 户)==： +丽 ～ 一 ， ㈩ 式中：是为气藏渗透率 ；为气藏孑L隙度 ； 为压力 ；为气体黏度 ；z为气体偏差因子 ；f为开井生产时间； 为气藏中的总井数 ；q 为第 口