CN108229051A - 一种预测油藏空气泡沫驱采收率的方法 - Google Patents

一种预测油藏空气泡沫驱采收率的方法 Download PDF

Info

Publication number
CN108229051A
CN108229051A CN201810049707.XA CN201810049707A CN108229051A CN 108229051 A CN108229051 A CN 108229051A CN 201810049707 A CN201810049707 A CN 201810049707A CN 108229051 A CN108229051 A CN 108229051A
Authority
CN
China
Prior art keywords
foam
gas
oil
recovery ratio
water
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201810049707.XA
Other languages
English (en)
Other versions
CN108229051B (zh
Inventor
尚庆华
赵习森
王玉霞
江绍静
杨永超
黄春霞
张冠华
陈龙龙
张建成
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shaanxi Yanchang Petroleum Group Co Ltd
Original Assignee
Shaanxi Yanchang Petroleum Group Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shaanxi Yanchang Petroleum Group Co Ltd filed Critical Shaanxi Yanchang Petroleum Group Co Ltd
Priority to CN201810049707.XA priority Critical patent/CN108229051B/zh
Publication of CN108229051A publication Critical patent/CN108229051A/zh
Application granted granted Critical
Publication of CN108229051B publication Critical patent/CN108229051B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00—Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16—Enhanced recovery methods for obtaining hydrocarbons

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mining & Mineral Resources (AREA)
  • Theoretical Computer Science (AREA)
  • Geology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Fluid Mechanics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Treating Waste Gases (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明属于油田开发技术领域，特别涉及一种预测油藏空气泡沫驱采收率的计算方法。一种预测油藏空气泡沫驱采收率的方法，包括如下步骤：（1）假设注入流体在地层内会逐渐形成三个区域，即泡沫区、液相区和气相区，并将驱替阶段分为气体突破前以及气体突破后，判断当前阶段属于驱替阶段的某一阶段；（2）待确定当前阶段为气体突破前或气体突破后，分别计算泡沫区、液相区和气相区各区域在当前阶段下的驱替效率和波及系数。（3）采收率R=驱替效率×波及系数。本发明将油藏中的驱替流体划分为三个区域，然后对每个区域的波及系数和驱油效率分别计算，最终得到空气泡沫驱的采收率，可行性和可操作性强。

Description

一种预测油藏空气泡沫驱采收率的方法
技术领域
本发明属于油田开发技术领域，特别涉及一种预测油藏空气泡沫驱采收率的计算方法。
背景技术
油藏采收率是油田开发中的一项重要评价指标。通过采收率预测可以评价油藏在不同驱动方式下的开采效果并计算其可采储量，对于油田开发规划和方案部署具有重要意义。根据国内外大量油藏开发的经验表明，影响油藏采收率的因素是多方面的，除油藏特征、天然能量、地层流体性质、井网类型等因素外，还与开发方式和驱动类型密切相关。
空气泡沫驱作为一项提高油藏采收率的有效技术方法，目前已在国内、外许多油田开发中取得了良好的应用效果。调研可知，目前预测空气泡沫驱采收率的方法主要有物理实验法和数值模拟法，缺乏完善有效的理论计算方法。物理实验法通过模拟油藏储层和流体条件，进行室内岩心驱替实验，根据实验数据计算预测采收率，耗费时间长，需要的成本高，且通常室内物理实验得到的只是一维驱替采收率，也即驱油效率，与实际油藏驱替存在波及系数情况下的采收率差别较大。而数值模拟法首先要建立精细准确的油藏地质模型，并在与生产历史数据进行良好拟合的基础上进行采收率的预测，需要大量详实有效的数据资料，前期准备工作量和后期计算处理工作量都很大，同样耗时也比较长。
发明内容
本发明旨在针对上述问题，提出一种将驱替阶段分为气体突破前以及气体突破后分阶段定量计算采收率的方法。
本发明的技术方案在于：
一种预测油藏空气泡沫驱采收率的方法，包括如下步骤：
（1）假设注入流体在地层内会逐渐形成三个区域，即泡沫区、液相区和气相区，并将驱替阶段分为气体突破前以及气体突破后，判断当前阶段属于驱替阶段的某一阶段；
（2）待确定当前阶段为气体突破前或气体突破后，分别计算泡沫区、液相区和气相区各区域在当前阶段下的驱替效率和波及系数；
（3）采收率R=驱替效率×波及系数。
当A btfoam +A btwater +A btgas< A total 时，即三个区域波及面积之和小于整个油藏面积时，认为气体尚未突破；当A btfoam +A btwater +A btgas= A total 时，即三个区域波及面积之和占满整个油藏时，则认为此时气体开始突破；当A btfoam +A btwater= A total 时，即水基液相区域与泡沫区域波及面积之和占满整个油藏时，则认为此时气相区完全突破；根据定义及渗流力学，上述各区域的面积可通过以下式计算：
（1）
（2）
（3）
其中，①C为和井网有关的系数；当采用五点井网时，C=0.718；当采用反七点井网时，C=0.743；当采用反九点井网时，C=0.525；
②、及分别为泡沫和原油的拟流度比、泡沫液和原油的拟流度比及气体和原油的拟流度比，其值分别取决于对应两种流体在均质条件下的流度比M和渗透率变异系数V，根据非均质条件下不同流度比流体的波及系数拟合关系，由以下式分别计算得到：
当V≤0.7时，
（4）
当V＞0.7时,
(5)
③计算A foam ：首先采用泡沫半衰期的数据计算泡沫剩余量：
（6）
其中,V left 为某个完整段塞经过时间t后的泡沫剩余体积；V i 为一个泡沫段塞刚形成时的体积，计算中假设其为地层压力下的注入空气和起泡体系溶液体积之和，因此也即每个段塞的流体总注入体积；T为泡沫半衰期；
泡沫区的波及体积V foam 可计算如下：
（7）
其中，S wc 为束缚水饱和度；S orfoam 为泡沫驱的残余油饱和度，其值等于泡沫与原油相对渗透率曲线中原油相对渗透率为0时的端点值；
泡沫区的波及面积如下：
（8）
水基液相区的波及面积因泡沫溃破而产生，根据其物理意义可由下式计算：
（9）
其中，φ为油藏孔隙度，h为油藏厚度，λ water 表示每个注入段塞中起泡体系溶液所占的体积比，S orwater 为泡沫液驱的残余油饱和度，其值等于泡沫液与原油相对渗透率曲线中原油相对渗透率为0时的端点值；
气相区的波及面积来自流体注入后未能形成泡沫的空气和泡沫溃破后产生的气体在前缘部分形成的波及区域，根据其物理意义可表示为：
（10）
其中，λ gas 表示每个注入段塞中气体所占的体积比，S orgas 为气驱的残余油饱和度，其值等于气体与原油相对渗透率曲线中原油相对渗透率为0时的端点值，R表示气体的反应系数，由低温氧化反应知R=0.996。
突破前采收率计算包括三部分，即泡沫区采收率、液相区采收率及气相区采收率；
其中，泡沫区采收率计算过程如下：
①计算波及系数
Ev foam = Ez foam ×Es foam （11）
其中，Es foam 为泡沫区平面波及系数，Ez foamm 为泡沫区纵向波及系数
（12）
（13）
其中，M为泡沫与原油流度比，为粘性力和重力的比值
（14）
式中，μ t 表示泡沫渗流速度，μ 0 为地层原油的粘度，χ为油藏沿注水方向的距离，K x 为横向渗透率，Δρ为驱替流体与被驱替流体的密度差；
②计算驱替效率
（15）
式中，S oi 为原始含油饱和度；
③计算采收率：
RF foam =Ed foam ×Ev foam （16）；
同理，液相区采收率过程如下：
Ev water =Ez water ×Es water （17）
其中，Es water 为泡沫液平面波及系数，Ez water 为泡沫液纵向波及系数；Ez water 由起泡体系溶液与原油的流度比代入式（13）计算得到；
（18）
（19）
液相区采收率RF water 为：
RF water = Ed water ×Ev water （20）
同理，气相区采收率过程如下：
Ev gas =Ez gas ×Es gas （21）
其中，Es gas 为气体平面波及系数，Ez gas 为气体纵向波及系数；Ez gas 由气体与原油的流度比代入式（13）计算得到；
（22）
（23）
气相区的采收率RF gas 为：
RF gas = Ed gas ×Ev gas （24）
气体突破前的总体采收率为：
RF= RF foam ＋RF water ＋RF gas 。
突破后采收率包括三部分，即泡沫区采收率、液相区采收率及气相区采收率；
（1）计算体积波及系数：
针对均质油藏，突破时的平面波及系数Es bt 如下式：
（25）
突破后的平面波及系数Es afterbt 如下式：
（26）
其中，D为关于井网类型的系数，当采用五点井网时，D=0.2749；当采用反七点井网时，D =0.2351；当采用反九点井网时，D=0.201；V ibt 为突破时注入流体的体积；
针对非均质油藏，突破时的平面波及系数Es bt 如下式：
（27）
突破后的平面波及系数Es afterbt 如下式：
（28）
其中，系数用来限制突破后的平面波及系数，使其趋向于当前流度比和渗透率变异系数条件下的最大平面波及系数Es max ，其计算方法如下：
（29）
突破后的纵向波及系数Ez afterbt 计算方法同上述式（13）；
体积波及系数
Ev afterbt =Ez afterbt ×Es afterbt （33）；
（2）计算动态驱油效率，动态驱油效率E d 为：
（34）
其中，S d 为驱替流体动态饱和度，为动态驱替下的残余油饱和度，；
根据渗流力学前缘含水饱和度理论，累计注入孔隙体积倍数PV为含水率对含水饱和度变化率的倒数，即
（35）
同时，根据渗流力学理论，含水率f w 可以表示为：
（36）
根据油和泡沫液两相相对渗透率K ro 和K rw 是含水饱和度的函数，通常又可以表示为：
（37）
其中a、b为系数，可以通过和的半对数坐标图线性回归求得；
将式（37）代入式（36）并求导，得到
（38）
（39）
上述累计注入孔隙体积倍数可根据累计注入流体的量确定，即：
（40）
其中B w 为地层水体积系数，A为含油面积，h为油藏厚度，φ为油藏孔隙度，ρ w 为泡沫液密度；
利用以上式子即可求得给定PV数的流度比，再由流度比可求得此时对应的含水率；然后根据渗流力学水相分流公式，可求得此时的动态平均含水饱和度，也即驱替流体动态饱和度S d ，即：
（41）
其中，μ w 为泡沫液粘度，μ 0 为地层原油粘度；
将其带入驱油效率的公式（34），即可求得此时油藏的动态平均驱油效率：
（42）
（3）计算采收率
根据计算得到的体积波及系数和驱油效率，得到采收率RF计算如下：
RF=Ed×Ev （43）。
所述的Es max 的计算方法如下：
引入渗透率变异系数由低到高变化时分界点处的值V c ，根据其与流度比的关系可由下式计算得到：
（30）
当渗透率变异系数V＜V c 时，最大平面波及系数为，
（31）
当渗透率变异系数V≥V c 时，最大平面波及系数为
（32）。
本发明的技术效果在于：
本发明根据空气泡沫驱油的机理将油藏中的驱替流体划分为三个区域，然后对每个区域的波及系数和驱油效率分别计算，最终得到空气泡沫驱的采收率。计算公式中的相关参数易于获得，相比通过室内驱替实验和油藏数值模拟获得空气泡沫驱采收率的方法，需要的工作量小，可行性和可操作性强。可以为油田开发效果评价及规划提供一种有效的技术方法。
附图说明
图1为空气泡沫驱采收率计算分区示意图。
具体实施方式
在空气泡沫驱时，由于泡沫自身的不稳定性，组成泡沫的气体与起泡体系会逐渐分离。随着不断地交替注入气体和水基起泡体系溶液，由于流体在地层内流度比的差别，假设注入流体会在地层内逐渐形成三个区域，即泡沫区、液相区和气相区。由于气体流度大于水基溶液流度，水基溶液流度大于泡沫流度，因此距注入井由近到远的顺序依次为泡沫区、液相区和气相区，具体如图1所示。图1中，A foam 、A water 和A gas 分别表示泡沫区、液相区和气相区形成时所占据的面积，A btfoam 、A btwater 和A btgas 则分别表示三个区域在当前的注入量下各自可以扩展占据的面积，A total 表示一个完整井网的面积。
根据上述分区和定义，将驱替阶段分为气体突破前和突破后，然后分别考虑各区域的驱替效率和波及系数的计算，最终得到空气泡沫驱不同阶段的采收率。在计算阶段采收率前，首先判断驱替处于哪一阶段。
如图1所示，当A btfoam +A btwater +A btgas< A total 时，即三个区域波及面积之和小于整个油藏面积时，认为气体尚未突破；当A btfoam +A btwater +A btgas= A total 时，即三个区域波及面积之和占满整个油藏时，则认为此时气体开始突破；当A btfoam +A btwater= A total 时，即水基液相区域与泡沫区域波及面积之和占满整个油藏时，则认为此时气相区完全突破；根据定义及渗流力学，上述各区域的面积可通过以下式计算：
（1）
（2）
（3）
其中，①C为和井网有关的系数；当采用五点井网时，C=0.718；当采用反七点井网时，C=0.743；当采用反九点井网时，C=0.525；
②、及分别为泡沫和原油的拟流度比、泡沫液和原油的拟流度比及气体和原油的拟流度比，其值分别取决于对应两种流体在均质条件下的流度比M和渗透率变异系数V，根据非均质条件下不同流度比流体的波及系数拟合关系，由以下式分别计算得到：
当V≤0.7时，
（4）
当V＞0.7时,
(5)
③计算A foam ：首先采用泡沫半衰期的数据计算泡沫剩余量：
（6）
其中,V left 为某个完整段塞经过时间t后的泡沫剩余体积；V i 为一个泡沫段塞刚形成时的体积，计算中假设其为地层压力下的注入空气和起泡体系溶液体积之和，因此也即每个段塞的流体总注入体积；T为泡沫半衰期；公式（6）表示单一段塞泡沫剩余量的计算方法，而现场采用分段注入的工作制度，因此采用将多个经过了不同时间的泡沫体积累加的方式计算泡沫剩余量。每个泡沫段塞的形成时刻，按形成该段塞的两种流体都完成注入的时刻来计算。因此由物理意义，泡沫区的波及体积V foam 可计算如下：
（7）
其中，S wc 为束缚水饱和度；S orfoam 为泡沫驱的残余油饱和度，其值等于泡沫与原油相对渗透率曲线中原油相对渗透率为0时的端点值；
得到泡沫区的波及体积之后，泡沫区的波及面积如下：
（8）
水基液相区的波及面积因泡沫溃破而产生，根据其物理意义可由下式计算：
（9）
其中，φ为油藏孔隙度，h为油藏厚度，λ water 表示每个注入段塞中起泡体系溶液所占的体积比，S orwater 为泡沫液驱的残余油饱和度，其值等于泡沫液与原油相对渗透率曲线中原油相对渗透率为0时的端点值；
气相区的波及面积来自流体注入后未能形成泡沫的空气和泡沫溃破后产生的气体在前缘部分形成的波及区域，根据其物理意义可表示为：
（10）
其中，λ gas 表示每个注入段塞中气体所占的体积比，S orgas 为气驱的残余油饱和度，其值等于气体与原油相对渗透率曲线中原油相对渗透率为0时的端点值，R表示气体的反应系数，由低温氧化反应知R=0.996。
在确定了驱替阶段后，分别考虑和计算突破前和突破后的采收率。
（一）突破前采收率
突破前采收率包括三部分，即泡沫区采收率、液相区采收率及气相区采收率；
（1）泡沫区采收率计算过程如下：由于泡沫稳定区域远远小于油藏面积，且其调剖封堵作用明显，故而泡沫区为以井底为圆心半径很小的圆环区域，始终不会突破。且由于突破前，注入体积倍数等于波及系数，所以其波及体积等于注入泡沫体积；
①计算波及系数
Ev foam = Ez foam ×Es foam （11）
其中，Es foam 为泡沫区平面波及系数，Ez foamm 为泡沫区纵向波及系数
（12）
（13）
其中，M为泡沫与原油流度比，为粘性力和重力的比值
（14）
式中，μ t 表示泡沫渗流速度，μ 0 为地层原油的粘度，χ为油藏沿注水方向的距离，K x 为横向渗透率，Δρ为驱替流体与被驱替流体的密度差；
②计算驱替效率
由于泡沫在空气泡沫驱中可以起到明显的封堵调剖作用，因此假设在泡沫区域，泡沫驱替为活塞式驱替，其驱替效率Ed foam 为定值，由如下公式确定
（15）
式中，S oi 为原始含油饱和度；
③计算采收率：
RF foam =Ed foam ×Ev foam （16）；
（2）液相区采收率计算过程如下：
液相区的体积波及系数如下所示：
Ev water =Ez water ×Es water （17）
其中，Es water 为泡沫液平面波及系数，Ez water 为泡沫液纵向波及系数；Ez water 由起泡体系溶液与原油的流度比代入式（13）计算得到；
（18）
根据前缘理论，在前缘到达生产井之前，前缘之前的范围，平均驱替流体饱和度总是相等，因而在突破前泡沫液动态驱替效率Ed water 应为常数，由下式确定
（19）
液相区采收率RF water 为：
RF water = Ed water ×Ev water （20）
（3）气相区采收率计算过程如下：
气相区的体积波及系数如下所示：
Ev gas =Ez gas ×Es gas （21）
其中，Es gas 为泡沫液平面波及系数，Ez gas 为泡沫液纵向波及系数；Ez gas 由气体与原油的流度比代入式（13）计算得到；
（22）
（23）
气相区的采收率RF gas 为：
RF gas = Ed gas ×Ev gas （24）
气体突破前的总体采收率为：
RF= RF foam ＋RF water ＋RF gas 。
（二）突破后采收率
（1）计算体积波及系数：
A：针对均质油藏，由油藏工程方法采用下式计算突破时的平面波及系数Es bt 如下式：
（25）
突破后的平面波及系数Es afterbt 如下式：
（26）
其中，D为关于井网类型的系数，当采用五点井网时，D=0.2749；当采用反七点井网时，D =0.2351；当采用反九点井网时，D=0.201；V ibt 为突破时注入流体的体积；
B：针对非均质油藏，采用对均质油藏公式中的流度比进行修正的方法计算其平面波及系数，引入拟流度比M*， 分别包括泡沫和原油的拟流度比、泡沫液和原油的拟流度比及气体和原油的拟流度比 ，分别计算突破时的泡沫区、液相区以及气相区的平面波及系数
其计算同上述式（4）和（5）；突破时的平面波及系数Es bt 如下式：
（27）
突破后的平面波及系数Es afterbt 如下式：
（28）
其中，系数用来限制突破后的平面波及系数，使其趋向于当前流度比和渗透率变异系数条件下的最大平面波及系数Es max ，其计算方法如下：
（29）
计算最大平面波及系数Es max 之前，首先引入渗透率变异系数由低到高变化时分界点处的值V c ，根据其与流度比的关系可由下式计算得到：
（30）
当渗透率变异系数V＜V c 时，最大平面波及系数为，
（31）
当渗透率变异系数V≥V c 时，最大平面波及系数为
（32）。
其中Es c 为渗透率变异系数等于V c 时的平面波及系数。
纵向波及系数Ez afterbt 计算方法同上述式（13）；
体积波及系数
Ev afterbt =Ez afterbt ×Es afterbt （33）；
（2）计算动态驱油效率，动态驱油效率E d 为：
（34）
其中，S d 为驱替流体动态饱和度，为动态驱替下的残余油饱和度，；
根据渗流力学前缘含水饱和度理论，累计注入孔隙体积倍数PV为含水率对含水饱和度变化率的倒数，即
（35）
同时，根据渗流力学理论，含水率f w 可以表示为：
（36）
根据油和泡沫液两相相对渗透率K ro 和K rw 是含水饱和度的函数，通常又可以表示为：
（37）
其中a、b为系数，可以通过和的半对数坐标图线性回归求得；
将式（37）代入式（36）并求导，得到
（38）
（39）
上述累计注入孔隙体积倍数可根据累计注入流体的量确定，即：
（40）
其中B w 为地层水体积系数，A为含油面积，h为油藏厚度，φ为油藏孔隙度，ρ w 为泡沫液密度；
利用以上式子即可求得给定PV数的流度比，再由流度比可求得此时对应的含水率；然后根据渗流力学水相分流公式，可求得此时的动态平均含水饱和度，也即驱替流体动态饱和度S d ，即：
（41）
其中，μ w 为泡沫液粘度，μ 0 为地层原油粘度；
将其带入驱油效率的公式（34），即可求得此时油藏的动态平均驱油效率：
（42）
（3）计算采收率
根据计算得到的体积波及系数和驱油效率，得到采收率RF计算如下：
RF=Ed×Ev （43）。
实施例1
应用油田实际生产数据对本一种预测油藏空气泡沫驱采收率的方法进行验证：
唐80井区位于陕北延安市甘谷驿镇境内，主要含油层位为三叠系延长组长6油层组，油藏平均埋深为441米，平均孔隙度7.9％，平均渗透率0.82×10-3μm2，原始地层压力系数0.95，油层温度26~30℃，属于超低渗、低压、低温岩性油藏。经过油藏筛选和分析，决定在该区开展空气泡沫驱油矿场试验。在矿场试验前，曾针对该区开展了大量的室内测试和评价研究，取得了比较全面的参数资料，详见表1。
表1 唐80井区室内测试参数值（地层条件下）
在室内综合研究的基础上，优选该区丛54井组和丛55井组开展了注空气泡沫驱油矿场试验，井组含油面积0.58km2，油层组有效厚度10m。在累计注入泡沫液494.27m3，累计注空气528.34m3（地层条件下体积）后受益油井开始见气，此时根据现场累计产油量数据计算空气泡沫驱阶段采收率为1.79%。截至统计数据时，累计注入泡沫液5909.32m3，累计注空气20209.24m3（地层条件下体积），井组一直处于见气阶段，根据现场累计产油量数据计算空气泡沫驱阶段采收率达到7.61%。
根据表1及上述相关参数，利用本一种预测油藏空气泡沫驱采收率的方法对上述井组注空气泡沫驱油不同阶段的采收率进行预测计算，在见气前即累计注入泡沫液494.27m3，累计注空气528.34m3（地层条件下体积）时的采收率计算结果为1.92%；在累计注入泡沫液5909.32m3，累计注空气20209.24m3（地层条件下体积）时的采收率计算结果为8.26%。两次预测计算结果与现场实际数据计算结果的绝对误差分别为0.13%和0.65%，相对误差分别为7.26%和8.54%，满足工程计算中相对误差控制在10%以内的精度要求。说明了本一种预测油藏空气泡沫驱采收率的方法的准确性和适用性，可以为油藏空气泡沫驱效果评价及相关规划提供依据和参考。
上述中已对涉及到的各个物理量所代表的含义进行了说明，其具体单位分别如下： A foam 、A water 、A gas 、A btfoam 、A btwater 、A btgas 、A total 、A等面积的单位均为㎡；V i 、V left 、V ibt 、V foam 等体积的单位均为㎥；μ w 和μ 0 等粘度的单位均为；Δρ和ρ w 等密度的单位均为㎏/㎥；h和χ等空间尺度的单位均为m；μ t 单位为m/d；K x 的单位为md；T的单位为min。流度比参数、、、M，饱和度参数S wc 、S oi 、S orfoam 、S owater 、S orgas 、、S d 、，波及系数Ev foam 、Es foam 、Ez foam 、Ev water 、Es water 、Ez water 、Es gas 、Ez gas 、Es bt 、Es afterbt 、Es max 、Es c 、Ev、Es、Ez，驱油效率Ed foam 、Ed water 、Ed gas 、Ed，采收率RF foam 、RF water 、RF gas 、RF，及B w 、φ、PV、K ro 、K rw 、fw等参数均为无因次。

Claims (5)

1.一种预测油藏空气泡沫驱采收率的方法，其特征在于：包括如下步骤：
（1）假设注入流体在地层内会逐渐形成三个区域，即泡沫区、液相区和气相区，并将驱替阶段分为气体突破前以及气体突破后，判断当前阶段属于驱替阶段的某一阶段；
（2）待确定当前阶段为气体突破前或气体突破后，分别计算泡沫区、液相区和气相区各区域在当前阶段下的驱替效率和波及系数；
（3）采收率R=驱替效率×波及系数。
2.根据权利要求1所述的一种预测油藏空气泡沫驱采收率的方法，其特征在于：A btfoam +A btwater +A btgas< A total 时，即三个区域波及面积之和小于整个油藏面积时，认为气体尚未突破；A btfoam +A btwater +A btgas= A total 时，即三个区域波及面积之和占满整个油藏时，则认为此时气体开始突破；当A btfoam +A btwater= A total 时，即水基液相区域与泡沫区域波及面积之和占满整个油藏时，则认为此时气相区完全突破；根据定义及渗流力学，上述各区域的面积可通过以下式计算：
（1）
（2）
（3）
其中，①C为和井网有关的系数；当采用五点井网时，C=0.718；当采用反七点井网时，C=0.743；当采用反九点井网时，C=0.525；
②、及分别为泡沫和原油的拟流度比、泡沫液和原油的拟流度比及气体和原油的拟流度比，其值分别取决于对应两种流体在均质条件下的流度比M和渗透率变异系数V，根据非均质条件下不同流度比流体的波及系数拟合关系，由以下式分别计算得到：
当V≤0.7时，
（4）
当V＞0.7时,
（5）
③计算A foam ：首先采用泡沫半衰期的数据计算泡沫剩余量：
（6）
其中,V left 为某个完整段塞经过时间t后的泡沫剩余体积；V i 为一个泡沫段塞刚形成时的体积，计算中假设其为地层压力下的注入空气和起泡体系溶液体积之和，因此也即每个段塞的流体总注入体积；T为泡沫半衰期；
泡沫区的波及体积V foam 可计算如下：
（7）
其中，S wc 为束缚水饱和度；S orfoam 为泡沫驱的残余油饱和度，其值等于泡沫与原油相对渗透率曲线中原油相对渗透率为0时的端点值；
泡沫区的波及面积如下：
（8）
水基液相区的波及面积因泡沫溃破而产生，根据其物理意义可由下式计算：
（9）
其中，φ为油藏孔隙度，h为油藏厚度，λ water 表示每个注入段塞中起泡体系溶液所占的体积比，S orwater 为泡沫液驱的残余油饱和度，其值等于泡沫液与原油相对渗透率曲线中原油相对渗透率为0时的端点值；
气相区的波及面积来自流体注入后未能形成泡沫的空气和泡沫溃破后产生的气体在前缘部分形成的波及区域，根据其物理意义可表示为：
（10）
其中，λ gas 表示每个注入段塞中气体所占的体积比，S orgas 为气驱的残余油饱和度，其值等于气体与原油相对渗透率曲线中原油相对渗透率为0时的端点值，R表示气体的反应系数，由低温氧化反应知R=0.996。
3.根据权利要求2所述的一种预测油藏空气泡沫驱采收率的方法，其特征在于：突破前采收率计算包括三部分，即泡沫区采收率、液相区采收率及气相区采收率；
其中，泡沫区采收率计算过程如下：
①计算波及系数
Ev foam = Ez foam ×Es foam （11）
其中，Es foam 为泡沫区平面波及系数，Ez foam 为泡沫区纵向波及系数
（12）
（13）
其中，M为泡沫与原油流度比，为粘性力和重力的比值
（14）
式中，μ t 表示泡沫渗流速度，μ 0 为地层原油的粘度，χ为油藏沿注水方向的距离，K x 为横向渗透率，Δρ为驱替流体与被驱替流体的密度差；
②计算驱替效率
（15）
式中，S oi 为原始含油饱和度；
③计算采收率：
RF foam =Ed foam ×Ev foam （16）；
同理，液相区采收率过程如下：
Ev water =Ez water ×Es water （17）
其中，Es water 为泡沫液平面波及系数，Ez water 为泡沫液纵向波及系数；Ez water 由起泡体系溶液与原油的流度比代入式（13）计算得到；
（18）
（19）
液相区采收率RF water 为：
RF water = Ed water ×Ev water （20）
同理，气相区采收率过程如下：
Ev gas =Ez gas ×Es gas （21）
其中，Es gas 为气体平面波及系数，Ez gas 为气体纵向波及系数；Ez gas 由气体与原油的流度比代入式（13）计算得到；
（22）
（23）
气相区的采收率RF gas 为：
RF gas = Ed gas ×Ev gas （24）
气体突破前的总体采收率为：
RF= RF foam ＋RF water ＋RF gas 。
4.根据权利要求3所述的一种预测油藏空气泡沫驱采收率的方法，其特征在于：突破后采收率包括三部分，即泡沫区采收率、液相区采收率及气相区采收率；
（1）计算体积波及系数：
针对均质油藏，突破时的平面波及系数Es bt 如下式：
（25）
突破后的平面波及系数Es afterbt 如下式：
（26）
其中，D为关于井网类型的系数，当采用五点井网时，D=0.2749；当采用反七点井网时，D =0.2351；当采用反九点井网时，D=0.201；V ibt 为突破时注入流体的体积；
针对非均质油藏，突破时的平面波及系数Es bt 如下式：
（27）
突破后的平面波及系数Es afterbt 如下式：
（28）
其中，系数用来限制突破后的平面波及系数，使其趋向于当前流度比和渗透率变异系数条件下的最大平面波及系数Es max ，其计算方法如下：
（29）
突破后的纵向波及系数Ez afterbt 计算方法同上述式（13）；
体积波及系数
Ev afterbt =Ez afterbt ×Es afterbt （33）；
（2）计算动态驱油效率，动态驱油效率E d 为：
（34）
其中，S d 为驱替流体动态饱和度，为动态驱替下的残余油饱和度，；
根据渗流力学前缘含水饱和度理论，累计注入孔隙体积倍数PV为含水率对含水饱和度变化率的倒数，即
（35）
同时，根据渗流力学理论，含水率f w 可以表示为：
（36）
根据油和泡沫液两相相对渗透率K ro 和K rw 是含水饱和度的函数，通常又可以表示为：
（37）
其中a、b为系数，可以通过和的半对数坐标图线性回归求得；
将式（37）代入式（36）并求导，得到
（38）
（39）
上述累计注入孔隙体积倍数可根据累计注入流体的量确定，即：
（40）
其中B w 为地层水体积系数，A为含油面积，h为油藏厚度，φ为油藏孔隙度，ρ w 为泡沫液密度；
利用以上式子即可求得给定PV数的流度比，再由流度比可求得此时对应的含水率；然后根据渗流力学水相分流公式，可求得此时的动态平均含水饱和度，也即驱替流体动态饱和度S d ，即：
（41）
其中，μ w 为泡沫液粘度，μ 0 为地层原油粘度；
将其带入驱油效率的公式（34），即可求得此时油藏的动态平均驱油效率：
（42）
（3）计算采收率
根据计算得到的体积波及系数和驱油效率，得到采收率RF计算如下：
RF=Ed×Ev （43）。
5.根据权利要求4所述的一种预测油藏空气泡沫驱采收率的方法，其特征在于：所述的Es max 的计算方法如下：
引入渗透率变异系数由低到高变化时分界点处的值V c ，根据其与流度比的关系可由下式计算得到：
（30）
当渗透率变异系数V＜V c 时，最大平面波及系数为，
（31）
当渗透率变异系数V≥V c 时，最大平面波及系数为
（32）。
CN201810049707.XA 2018-01-18 2018-01-18 一种预测油藏空气泡沫驱采收率的方法 Expired - Fee Related CN108229051B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810049707.XA CN108229051B (zh) 2018-01-18 2018-01-18 一种预测油藏空气泡沫驱采收率的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810049707.XA CN108229051B (zh) 2018-01-18 2018-01-18 一种预测油藏空气泡沫驱采收率的方法

Publications (2)

Publication Number Publication Date
CN108229051A true CN108229051A (zh) 2018-06-29
CN108229051B CN108229051B (zh) 2021-05-11

Family

ID=62667969

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810049707.XA Expired - Fee Related CN108229051B (zh) 2018-01-18 2018-01-18 一种预测油藏空气泡沫驱采收率的方法

Country Status (1)

Country Link
CN (1) CN108229051B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110656924A (zh) * 2019-08-29 2020-01-07 长江大学 一种超低渗透油藏分类方法
CN111350498A (zh) * 2020-04-13 2020-06-30 西南石油大学 一种中高渗油藏特高含水开发期弱驱分布特征描述的方法
CN111561301A (zh) * 2020-07-01 2020-08-21 西南石油大学 一种确定稠油油藏注水开发的原油黏度界限方法
CN110242263B (zh) * 2019-05-23 2021-03-16 中国石油天然气股份有限公司 二三结合开发模式下的采收率计算方法
CN112796717A (zh) * 2021-04-12 2021-05-14 西安石油大油气科技有限公司 井下稠油开采用减少气窜的空气泡沫驱方法及其起泡剂
CN113356806A (zh) * 2020-03-04 2021-09-07 中国石油化工股份有限公司 基于气驱波及系数调整注气量的方法、装置及电子设备
CN114113550A (zh) * 2021-08-17 2022-03-01 大连理工大学 一种利用二氧化碳微气泡提高石油采收率的实验方法
CN114136838A (zh) * 2021-11-19 2022-03-04 中国海洋石油集团有限公司 确定海上稠油不同含水阶段注热水驱原油粘度界限的方法
CN114135258A (zh) * 2021-11-30 2022-03-04 常州大学 一种顶部注气过程中气顶形成及扩大规律的室内判定方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010048715A1 (en) * 2008-10-30 2010-05-06 St. Francis Xavier University Gemini surfactants
CN102020981A (zh) * 2010-10-27 2011-04-20 山东大学 一种耐温抗盐低张力泡沫驱油剂及其制备方法
CN102434137A (zh) * 2011-12-16 2012-05-02 中国石油天然气股份有限公司 超低界面张力耦合式空气泡沫驱油方法
CN102606117A (zh) * 2012-04-05 2012-07-25 中国石油天然气股份有限公司 中深层低渗透稀油油藏的空气氮气交替驱开采方法
EP2456951B1 (en) * 2009-07-22 2014-05-21 Bergen Teknologioverføring AS Method for integrated enhanced oil recovery from heterogeneous reservoirs
CN104091069A (zh) * 2014-07-07 2014-10-08 中国海洋石油总公司 确定非均质储层各层位和位置驱油效率和波及系数的方法
CN104437236A (zh) * 2014-11-07 2015-03-25 陕西延长石油(集团)有限责任公司研究院 季铵盐型双子表面活性剂的制备及提纯方法
CN104832144A (zh) * 2015-05-25 2015-08-12 南开大学 一种微生物辅助的空气泡沫驱提高石油采收率的方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010048715A1 (en) * 2008-10-30 2010-05-06 St. Francis Xavier University Gemini surfactants
EP2456951B1 (en) * 2009-07-22 2014-05-21 Bergen Teknologioverføring AS Method for integrated enhanced oil recovery from heterogeneous reservoirs
CN102020981A (zh) * 2010-10-27 2011-04-20 山东大学 一种耐温抗盐低张力泡沫驱油剂及其制备方法
CN102434137A (zh) * 2011-12-16 2012-05-02 中国石油天然气股份有限公司 超低界面张力耦合式空气泡沫驱油方法
CN102606117A (zh) * 2012-04-05 2012-07-25 中国石油天然气股份有限公司 中深层低渗透稀油油藏的空气氮气交替驱开采方法
CN104091069A (zh) * 2014-07-07 2014-10-08 中国海洋石油总公司 确定非均质储层各层位和位置驱油效率和波及系数的方法
CN104437236A (zh) * 2014-11-07 2015-03-25 陕西延长石油(集团)有限责任公司研究院 季铵盐型双子表面活性剂的制备及提纯方法
CN104832144A (zh) * 2015-05-25 2015-08-12 南开大学 一种微生物辅助的空气泡沫驱提高石油采收率的方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
PENGCHENGLIU ET AL: "Enhanced oil recovery by air-foam flooding system in tight oil reservoirs: Study on the profile-controlling mechanisms", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 *
何金钢: "泡沫物理性能表征和泡沫驱油效果研究", 《中国博士学位论文全文数据库电子期刊 工程科技I辑》 *
李娜: "孤东油田空气泡沫驱提高采收率技术研究", 《中国优秀硕士学位论文全文数据库电子期刊 工程科技I辑》 *
王其伟: "泡沫驱提高原油采收率及对环境的影响研究", 《中国博士论文全文数据库电子期刊 工程科技I辑》 *

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110242263B (zh) * 2019-05-23 2021-03-16 中国石油天然气股份有限公司 二三结合开发模式下的采收率计算方法
CN110656924B (zh) * 2019-08-29 2023-08-22 长江大学 一种超低渗透油藏分类方法
CN110656924A (zh) * 2019-08-29 2020-01-07 长江大学 一种超低渗透油藏分类方法
CN113356806A (zh) * 2020-03-04 2021-09-07 中国石油化工股份有限公司 基于气驱波及系数调整注气量的方法、装置及电子设备
CN111350498B (zh) * 2020-04-13 2022-05-03 西南石油大学 一种中高渗油藏特高含水开发期弱驱分布特征描述的方法
CN111350498A (zh) * 2020-04-13 2020-06-30 西南石油大学 一种中高渗油藏特高含水开发期弱驱分布特征描述的方法
CN111561301A (zh) * 2020-07-01 2020-08-21 西南石油大学 一种确定稠油油藏注水开发的原油黏度界限方法
CN112796717A (zh) * 2021-04-12 2021-05-14 西安石油大油气科技有限公司 井下稠油开采用减少气窜的空气泡沫驱方法及其起泡剂
CN112796717B (zh) * 2021-04-12 2021-07-09 西安石油大油气科技有限公司 井下稠油开采用减少气窜的空气泡沫驱方法及其起泡剂
CN114113550A (zh) * 2021-08-17 2022-03-01 大连理工大学 一种利用二氧化碳微气泡提高石油采收率的实验方法
CN114113550B (zh) * 2021-08-17 2023-04-25 大连理工大学 一种利用二氧化碳微气泡提高石油采收率的实验方法
CN114136838A (zh) * 2021-11-19 2022-03-04 中国海洋石油集团有限公司 确定海上稠油不同含水阶段注热水驱原油粘度界限的方法
CN114136838B (zh) * 2021-11-19 2023-11-17 中国海洋石油集团有限公司 确定海上稠油不同含水阶段注热水驱原油粘度界限的方法
CN114135258A (zh) * 2021-11-30 2022-03-04 常州大学 一种顶部注气过程中气顶形成及扩大规律的室内判定方法
CN114135258B (zh) * 2021-11-30 2023-10-20 常州大学 一种顶部注气过程中气顶形成及扩大规律的室内判定方法

Also Published As

Publication number Publication date
CN108229051B (zh) 2021-05-11

Similar Documents

Publication Publication Date Title
CN108229051A (zh) 一种预测油藏空气泡沫驱采收率的方法
US7289942B2 (en) Performance prediction method for hydrocarbon recovery processes
Stone et al. Thermal simulation with multisegment wells
CN104989341B (zh) 一种确定低渗透油藏有效驱替注采井距的方法
CN105089595B (zh) 水平压裂裂缝导流作用下的油藏数值模拟方法及装置
CN105089582B (zh) 基于井下流量控制设备的油藏数值模拟方法及装置
WO2013003269A2 (en) Method for generating a general enhanced oil recovery and waterflood forecasting model
CN105243182A (zh) 致密油压裂水平井的动态储量计算方法
CN112392472B (zh) 确定页岩与邻近油层一体化开发方式的方法及装置
WO2011149779A1 (en) System and method for enhancing oil recovery from a subterranean reservoir
CN104389568B (zh) 蒸汽辅助重力泄油过程中气体辅助用量的获取方法及装置
CN105317407B (zh) 一种特高含水期表外储层的开发方法
CN106951649A (zh) 一种测定水平井sagd蒸汽腔扩展速度的方法
CN109858177A (zh) 一种基于拟流线法的底水油藏水平井水驱数值模拟方法
You et al. Production data analysis of shale gas using fractal model and fuzzy theory: Evaluating fracturing heterogeneity
CN104989385A (zh) 基于表皮系数计算的高温高压油气直井射孔参数优化方法
Shahverdi Characterization of three-phase flow and WAG injection in oil reservoirs
Wu et al. Inflow performance of a cyclic-steam-stimulated horizontal well under the influence of gravity drainage
CN104727789B (zh) 中高渗砂岩油藏水驱波及系数及过水倍数动态描述方法
Hsu et al. Field-Scale CO2-FIood Simulations and Their Impact on the Performance of the Wasson Denver Unit
Mathews et al. Fractal methods improve Mitsue miscible predictions
Rodriguez Inferences of two dynamic processes on recovery factor and well spacing for a shale oil reservoir
RU2558549C1 (ru) Способ исследования и интерпретации результатов исследования скважины
Chen et al. Effects of capillary and viscous forces on two-phase fluid displacement in the microfluidic model
Rassenfoss Finding pathways to produce heavy oil from Canadian carbonates

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210511

CF01 Termination of patent right due to non-payment of annual fee