CN115422859B - 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法 - Google Patents
一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法 Download PDFInfo
- Publication number
- CN115422859B CN115422859B CN202211385246.6A CN202211385246A CN115422859B CN 115422859 B CN115422859 B CN 115422859B CN 202211385246 A CN202211385246 A CN 202211385246A CN 115422859 B CN115422859 B CN 115422859B
- Authority
- CN
- China
- Prior art keywords
- steam
- calculating
- oil
- vertical
- thick
- 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.)
- Active
Links
- 238000010793 Steam injection (oil industry) Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 27
- 239000003921 oil Substances 0.000 claims description 151
- 238000004519 manufacturing process Methods 0.000 claims description 57
- 238000004364 calculation method Methods 0.000 claims description 50
- 230000005484 gravity Effects 0.000 claims description 49
- 238000002347 injection Methods 0.000 claims description 48
- 239000007924 injection Substances 0.000 claims description 48
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 20
- 230000035699 permeability Effects 0.000 claims description 18
- 238000006073 displacement reaction Methods 0.000 claims description 11
- 230000001133 acceleration Effects 0.000 claims description 10
- 230000009471 action Effects 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000013508 migration Methods 0.000 claims description 9
- 230000005012 migration Effects 0.000 claims description 9
- 239000011148 porous material Substances 0.000 claims description 9
- 239000010779 crude oil Substances 0.000 claims description 7
- 239000008346 aqueous phase Substances 0.000 claims description 4
- 239000012071 phase Substances 0.000 claims description 4
- 238000011084 recovery Methods 0.000 abstract description 12
- 238000010795 Steam Flooding Methods 0.000 abstract description 8
- 238000013461 design Methods 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract 2
- 239000000295 fuel oil Substances 0.000 description 20
- 238000009826 distribution Methods 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 4
- 230000001186 cumulative effect Effects 0.000 description 4
- 239000000203 mixture Substances 0.000 description 4
- 238000005457 optimization Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000010796 Steam-assisted gravity drainage Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003129 oil well Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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
- E21B43/24—Enhanced recovery methods for obtaining hydrocarbons using heat, e.g. steam injection
- E21B43/2406—Steam assisted gravity drainage [SAGD]
-
- 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
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法,对于目标区块处,包括:S1收集地质油藏基本参数;S2计算蒸汽微元在油层垂直向的受力值;S3计算蒸汽微元在油层水平向的受力值;S4计算蒸汽微元在油藏垂直方向的合力和瞬时渗流速度;S5计算蒸汽微元在油藏水平向的合力和瞬时渗流速度;S6计算垂向分流系数,绘制垂向分流系数与井间距离的图版;S7计算注入蒸汽的超覆程度,绘制超覆程度与井间距离的图版;S8定量计算倾斜厚层稠油注蒸汽纵向波及系数。本发明达到的有益效果是:实现了在海上大井距热采且无观察井的条件下,定量预测倾斜厚层稠油注蒸汽吞吐或者蒸汽驱纵向波及系数,从而有助于稠油热采方案设计。
Description
技术领域
本发明涉及油井勘探开采技术领域,特别是一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法。
背景技术
室内研究和矿场实践表明,热力采油技术(如蒸汽吞吐、蒸汽驱、蒸汽辅助重力泄油等)是一项可以大幅度提高稠油油田采收率的有效技术,是目前陆上稠油开发的有效手段。渤海目前稠油热采处于规模化推广阶段,稠油热采油藏的粘度不断增加,油藏厚度不断扩大(地层条件原油粘度从50mPa•s到50000mPa•s,油藏单砂体厚度从6m扩展到40m),定量获取厚层稠油油藏注蒸汽吞吐或蒸汽驱纵向波及系数,有助于评价注蒸汽热采(包括蒸汽吞吐、蒸汽驱等)的有效加热范围、热采渗流模式,从而指导热采方案设计和注采参数优化。
目前关于注蒸汽的纵向波及系数,陆地油田通常在井间设置观察井,通过观察井适时监控地层中的温度,从而通过地层中的温度判断纵向上的动用情况。而海上稠油热采,由于开发成本高,通常采用大井距热采,且井间没有观察井,因此定量获取厚层稠油油藏注入蒸汽纵向波及系数成为难题。
发明内容
本发明的目的在于克服现有技术的缺点,提供一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法,解决了定量获取厚层稠油油藏注入蒸汽纵向波及系数的难题。
本发明的目的通过以下技术方案来实现:一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法,包括以下步骤:
S1、收集目标区块地质油藏基本参数;
S2、计算目标区块蒸汽微元在油层垂直方向上的受力值;
S3、计算目标区块蒸汽微元在油层水平方向上的受力值;
S4、计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度;
S5、计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度;
S6、计算垂向分流系数,并绘制垂向分流系数与井间距离的图版;
S7、计算注入蒸汽的超覆程度,并绘制超覆程度与井间距离的图版;蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值,用D表示;
S8、定量计算倾斜厚层稠油注蒸汽纵向波及系数。
进一步地,所述S1中的基本参数:是指收集目标区块的基础数据,包括油层厚度、地层倾角、注采井间距离、油的密度、水的密度、气体粘度、界面张力、接触角、孔隙半径、垂直方向渗透率、水平方向渗透率。为通过受力分析和渗流规律计算注入地层中蒸汽的分布做好数据准备。
进一步地,所述S2中计算目标区块蒸汽微元在油层垂直方向上的受力值,包括重力、浮力;
S21、重力的计算
S22、浮力的计算;
在注入蒸汽突破之前,蒸汽微元受到的浮力F浮力看作是恒定的,方向为竖直向上;
在步骤S2中,通过计算的重力和浮力值是决定注入蒸汽微元在垂向上运动的动力,由于考虑了浮力和地层倾角的影响,可使描述倾斜油藏垂直油层方向的蒸汽运动更加精确。
进一步地,所述S3中计算目标区块蒸汽微元在油层水平方向上的受力值,包括毛管力和驱替压力:
S31、计算驱替压力,即计算注采井间主流线上的注采压力梯度,公式为
S32、计算毛管力
在水湿情况下,毛细管力是注入水运移的动力,取正值;在油湿情况下,毛管力是注入水运移的阻力,取负值;
在步骤S3中,通过注入蒸汽在水平方向上受力主要包括驱替力和毛管力,是决定蒸汽在水平方向上运动的动力,由于考虑了毛管力和地层倾角的影响,可使描述倾斜油藏油层水平方向的蒸汽运动更加精确。
进一步地,所述步骤S4计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度;
S41、计算目标区块蒸汽微元在油藏垂直方向的合力
根据蒸汽微元在Y方向的重力分量FY重力,蒸汽微元在Y方向的重力分量FY浮力,垂直油藏方向的垂直毛细管力FY毛管力,计算得到蒸汽微元在油藏垂直方向的合力FY合力;
S42、计算目标区块蒸汽微元在油藏垂直方向的瞬时渗流速度
注入蒸汽沿油层斜上方渗流,可将速度分解为沿油层X方向的速度和Y方向的速度;
蒸汽在重力、浮力及毛管力的共同作用下,根据达西定律,计算得到Y方向瞬时渗流速度,即在油藏垂直方向的瞬时渗流速度VX。
在步骤S4中,由于是倾斜油藏,并且S2、S3计算的各种力在油藏垂直方向均有分力,且作用方向不一致,因此需要求取S2和S3在垂直方向的合力,并根据合力和渗流理论,计算得到油藏垂直方向的瞬时渗流速度。由于考虑了多种作用力的合力和地层倾角的影响,让垂直反向的瞬时渗流速度计算更加精确。
进一步地,所述S5计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度;
S51、计算目标区块蒸汽微元在油藏水平方向的合力
根据蒸汽微元在X方向的重力分量FX重力,蒸汽微元在X方向的重力分量FX浮力,驱替压力,垂直油藏方向的垂直毛细管力FX毛管力,计算得到蒸汽微元在油藏垂直方向的合力FX合力;
S52、计算目标区块蒸汽微元在油藏水平方向的瞬时渗流速度
根据达西定律,油滴在浮力及毛管力的共同作用下,计算得到X方向瞬时渗流速度为VY。
在步骤S5中,由于是倾斜油藏,第S2和S3计算的各种力在油藏水平方向均有分力,且作用方向不一致,因此需要求取S2和S3在水平方向的合力,并根据合力和渗流理论,计算得到油藏水平方向的瞬时渗流速度。由于考虑了多种作用力的合力和地层倾角的影响,让水平方向瞬时渗流速度计算更加精确。
进一步地,所述S6计算垂向分流系数,并绘制垂向分流系数与井间距离的图版;
S61、计算垂向分流系数
垂向分流系数为垂向瞬时渗流速度与总瞬时渗流速度的比值,计算方程为
在步骤S6中,根据S4、S5计算得到的垂直方向和水平方向瞬时渗流速度的计算结果,可以通过计算垂向分流系数,定量计算出注入蒸汽在垂向上的突进情况,解决了海上油田由于没有观察井,难以掌握注入蒸汽油层内分布情况的难题。
进一步地,所述S7计算注入蒸汽的超覆程度,并绘制超覆程度与井间距离的图版;
S71、计算注入蒸汽的超覆程度
蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值,用D表示;
将注采井间距离划分为N块,当蒸汽前缘运移到注采井间 i 块( 0 < i < N)时,利用此处的径向蒸汽运移量与垂向流量比的乘积和第i-1块的超覆程度迭代,计算出此时的蒸汽超覆程度,计算公式为
S72、绘制超覆程度与井间距离的图版
由于S6计算的垂向分流量值是一个变化的量,在注采井间的不同位置具有不同的值,通过S7,可以描绘出注采井间不同位置处的蒸汽分布情况
进一步地,所述S8定量计算倾斜厚层稠油注蒸汽纵向波及系数;
S7计算的蒸汽超覆程度Di,在注采井间的不同位置具有不同的值,通过S8,可以描绘出注采井间任何剖面的蒸汽分布面积,根据面积加权,可以定量计算出注采井间的蒸汽波及厚度。这样得到的波及厚度和波及系数值,考虑了注采井间不同位置的超覆程度差异,使得计算结果更加贴近矿场实际情况。
本发明具有以下优点:
(1)基于注入蒸汽在地层中的受力分析,提出了通过计算注采井间距离处的纵向分流系数和不同位置的蒸汽超覆程度来预测倾斜厚层稠油油藏的纵向波及系数,能有效解决海上大井距热采且无观察井的条件下,定量预测倾斜厚层稠油注蒸汽吞吐或者蒸汽驱纵向波及系数的难题,为海上稠油油田的热采方案设计提供了支撑;
(2)传统的注蒸汽的纵向波及系数,陆地油田通常在井间设置观察井,通过观察井适时监控地层中的温度,从而通过地层中的温度判断纵向上的动用情况;
通俗地说,传统的注蒸汽的纵向波及系数,是通过经验情况进行实验的,最终获得相应系数;根据经验的方式,使得系数的结果非常不准确;若地质情况复杂,以经验的方式测得的数据不可靠;
本发明提供的测量方式,是在对注入蒸汽在地层中受力分析的基础上,定量表征出倾斜厚层稠油油藏注蒸汽开发不同方向上的瞬时渗流速度,根据距离注汽井点的不同距离计算出垂向分流系数,并计算得到注采井间不同位置的蒸汽超覆程度,从而定量计算出倾斜厚层稠油注入蒸汽的纵向波及系数;
本发明提供的测量方式,通过计算的方式获得系数,不仅中间数据计算非常精确,而且最后经计算获得的结果也非常精确;为确定海上厚层稠油油藏注蒸汽吞吐和蒸汽驱方案设计、注采参数优化、水平井纵向位置优化等提供了非常可靠的依据;
(3)本方案的测量方式,不凭借经验,而是实际计算出来的,操作性强,应用广泛。
附图说明
图1为定量预测倾斜厚层稠油注蒸汽纵向波及系数的主要步骤;
图2为实施例1倾斜厚层稠油油藏(润湿性:油湿)地层中蒸汽微元受力示意图;
图3为实施例1注采井间不同距离处的垂向分流量图版;
图4为实施例1注采井间不同距离处的注入蒸汽超覆程度图版;
图5为实施例2中典型油田注采井间不同距离处的垂向分流量图版;
图6为实施例2中典型油田注采井间不同距离处的注入蒸汽超覆程度图版;
表1 为实施例1中A目标区块的基础数据;
表2 为实施例1中A目标区注采井间不同距离处的垂向分流量和蒸汽超覆程度计算表;
表3 为实施例2中B目标区块的基础数据;
表4 为实施2中B目标区块注采井间不同距离处的垂向分流量和蒸汽超覆程度计算表。
具体实施方式
下面结合附图对本发明做进一步的描述,但本发明的保护范围不局限于以下所述。
(实施例1)
一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法,主要包括以下步骤(图1):
S1,收集目标区块地质油藏基本参数等静态资料。收集目标区块的基础数据,包括油层厚度、地层原油粘度、渗透率等参数(表1),为通过受力分析和渗流规律计算注入地层中蒸汽的分布做好数据准备。表1中收集的为A目标区的基础数据。
表1 A目标区块的基础数据
S2,计算目标区块蒸汽微元在油层垂直方向上的重力和浮力值。计算的重力和浮力值是决定注入蒸汽微元在垂向上运动的动力,由于考虑了浮力和地层倾角的影响,可使描述倾斜油藏垂直油层方向的蒸汽运动更加精确。如图2所示,厚层稠油油藏地层中某蒸汽微元主要受力有浮力、毛管力、重力和驱替压力。
对于蒸汽微元,受到的重力为:
在注入蒸汽突破以前,蒸汽微元受到的浮力可以近似看作是恒定的,方向为垂直向上,表达式为:
根据表1的数据,油的密度取值0.95g /cm3,水的密度取值1.0g /cm3,重力加速度取值9.8 m/s2,地层倾角取值30°,蒸汽微元的体积取值1.0cm3。根据式子(A-2)、(A-3)、(A-4)、(A-5)计算得到单位体积(1.0cm3,)蒸汽所受到的X方向重力和浮力分别为:-4.655mN和-0.245mN;Y方向的重力和浮力分别为8.063mN和-0.424mN。
S3,计算目标区块蒸汽微元在油层水平方向上的毛管力和驱替压力梯度值。注入蒸汽在水平方向上受力主要包括驱替力和毛管力,是决定蒸汽在水平方向上运动的动力,由于考虑了毛管力和地层倾角的影响,可使描述倾斜油藏油层水平方向的蒸汽运动更加精确。
用(A-7)计算注采井间主流线上的注采压力梯度。
式中:为生产井产量,cm3 /s;为蒸汽注入量,cm3 /s;为油层厚度,m;r为任意点到生产井的距离,cm;为 r 处的压力梯度,10-1MPa/cm; L 为注采井间距离,cm;K 为渗透率,;为原油粘度,mPa·s。
对于蒸汽,受到的毛管力为:
在水湿情况下,毛细管力是注入水运移的动力,取正值;在油湿情况下,毛管力是注入水运移的阻力,取负值。毛管力的方向与注入水运移方向同向或反向,为便于计算与研究,分解为平行油藏方向的水平毛管力和垂直油藏方向的垂直毛管力。
根据表1的数据,界面张力取值0.015 N/m,接触角取值150°,孔隙半径取值0.00005m。由于接触角较大,可以判定为油湿,毛管力为注入水运移的阻力,取负值。根据式子(A-9)、(A-10)计算得到蒸汽所受到的X和Y方向的毛管力分别为0.0015Pa和-0.003Pa。
S4,计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度。由于是倾斜油藏,S2和S3计算的各种力在油藏垂直方向均有分力,且作用方向不一致,因此需要求取S2和S3在垂直方向的合力,并根据合力和渗流理论,计算得到油藏垂直方向的瞬时渗流速度。由于考虑了多种作用力的合力和地层倾角的影响,让垂直反向的瞬时渗流速度计算更加精确。
联立式子(A-3)、(A-6)、(A-10),得到Y方向合力计算方程:
注入蒸汽沿油层斜上方渗流,可将速度分解为沿油层X方向的速度和Y方向的速度。根据达西定律,蒸汽在重力、浮力及毛管力的共同作用下,Y方向瞬时渗流速度为:
根据表1的数据,垂直方向渗透率取值500ⅹ10-3µm2,气体粘度取值5 mPa·s,根据式子(A-11)、(A-12)计算得到蒸汽Y方向的瞬时渗流速度为763.534 µm2/s。
S5,计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度。由于是倾斜油藏,S2和S3计算的各种力在油藏水平方向均有分力,且作用方向不一致,因此需要求取S2和S3在水平方向的合力,并根据合力和渗流理论,计算得到油藏水平方向的瞬时渗流速度。由于考虑了多种作用力的合力和地层倾角的影响,让水平方向瞬时渗流速度计算更加精确。
联立式子(A-2)、(A-5)、(A-7)、(A-9),得到X方向合力计算方程:
将(A-2)、(A-5)、(A-5)计算结果代入(A-13),得到:
根据达西定律,油滴在浮力及毛管力的共同作用下,X方向瞬时渗流速度为:
根据表1的数据,水平方向渗透率取值2000ⅹ10-3µm2,气体粘度取值5 mPa·s,根据式子(A-13)、(A-14)计算得到蒸汽X方向的瞬时渗流速度计算方程为:
S6,计算垂向分流系数,并绘制垂向分流系数与井间距离的图版。其特点在于,根据S4和S5计算得到的垂直方向和水平方向瞬时渗流速度的计算结果,可以通过计算垂向分流系数,定量计算出注入蒸汽在垂向上的突进情况,解决了海上油田由于没有观察井,难以掌握注入蒸汽油层内分布情况的难题。
根据表1的数据和A区块的实际情况,为生产井产量,取值34722.22cm3 /s(300m3 /d);为蒸汽注入量,取值34722.22cm3 /s(300m3 /d);为油层厚度,取值40m;L为注采井间距离,取值15000cm;K 为渗透率,取值;为蒸汽粘度,取值5mPa·s。将上述取值代入(A-17),计算方程变为:
S7,计算注入蒸汽的超覆程度,并绘制超覆程度与井间距离的图版。S6计算的垂向分流量值是一个变化的量,在注采井间的不同位置具有不同的值,通过S7,可以描绘出注采井间不同位置处的蒸汽分布情况。其特点在于,蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值,用 D表示(表2)。具体计算方法为:将注采井间的距离划分为N块,当蒸汽前缘运移到注采井间第 i 块( 0 < i < N) 时,利用此处的径向蒸汽运移量与垂向流量比的乘积和第i-1块的超覆程度Di-1迭代,并用(A-19)计算出此时的蒸汽超覆程度Di。并以为纵坐标,井间距离为横坐标,绘制超覆程度与注采井间距离L的图版(图4)。其中,表示第块的垂向分流量。
表2 A目标区注采井间不同距离处的垂向分流量和蒸汽超覆程度计算表
S8,定量计算倾斜厚层稠油注蒸汽纵向波及系数。S7计算的蒸汽超覆程度Di,在注采井间的不同位置具有不同的值,通过S8,可以描绘出注采井间任何剖面的蒸汽分布面积,根据面积加权,可以定量计算出注采井间的蒸汽波及厚度。这样得到的波及厚度和波及系数值,考虑了注采井间不同位置的超覆程度差异,使得计算结果更加贴近矿场实际情况。
根据S7的图版,用(A-19)计算厚层稠油蒸汽吞吐初期不同分块的蒸汽超覆程度(表2第5列);通过(1-Di)计算每个区块纵向波及无因次厚度hi(表2第6列);用每个区块的长度(10m)乘以无因次厚度1得到每个区块的面积Ai(表2第7列);用Ai hi计算得到每个区块的每个区块面积厚度乘积(表2第8列)。
根据表2的数据,根据(A-20)计算得到面积加权后的无因次厚度为65.3(表2第8列求和)/150(表2第7列求和)=0.435;无因次厚度h取值为1,根据(A-20)计算得到纵向厚度波及系数为43.5%。
(实施例2)
一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法,对某厚层倾斜稠油油藏B进行纵向波及系数预测,过程如下:
S1,收集目标区块地质油藏基本参数等静态资料。收集目标区块的基础数据,包括油层厚度、地层原油粘度、渗透率等参数(表3)。表3中收集的为B目标区的基础数据。
表3 B目标区的基础数据
S2,计算目标区块蒸汽微元在油层垂直方向上的受力值。对于蒸汽微元,受到的重力为:
在注入蒸汽突破以前,蒸汽微元受到的浮力可以近似看作是恒定的,方向为垂直向上,表达式为:
根据表3的数据,油的密度取值0.95g /cm3,水的密度取值1.0g /cm3,重力加速度取值9.8 m/s2,地层倾角取值20°,蒸汽微元的体积取值1.0cm3。根据式子(A-2)、(A-3)、(A-4)、(A-5)计算得到单位体积(1.0cm3,)蒸汽所受到的X方向重力和浮力分别为:-3.184mN和-0.168mN;Y方向的重力和浮力分别为8.749mN和-0.460mN。
S3,计算目标区块蒸汽微元在油层水平方向上的受力值。主要包括驱替力和毛管力,用(A-7)计算注采井间主流线上的注采压力梯度。
式中:为生产井产量,cm3 /s;为蒸汽注入量,cm3 /s;为油层厚度,m;r为任意点到生产井的距离,cm;为 r 处的压力梯度,10-1MPa/cm; L 为注采井间距离,cm;K 为渗透率,;为原油粘度,mPa·s。
对于蒸汽,受到的毛管力为:
在水湿情况下,毛细管力是注入水运移的动力,取正值;在油湿情况下,毛管力是注入水运移的阻力,取负值。毛管力的方向与注入水运移方向同向或反向,为便于计算与研究,分解为平行油藏方向的水平毛管力和垂直油藏方向的垂直毛管力。
根据表3的数据,界面张力取值0.02 N/m,接触角取值100°,孔隙半径取值0.00005m。由于接触角较大,可以判定为油湿,毛管力为注入水运移的阻力,取负值。根据式子(A-9)、(A-10)计算得到蒸汽所受到的X和Y方向的毛管力分别为0.001Pa和-0.001Pa。
S4,计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度。联立式子(A-3)、(A-6)、(A-10),得到Y方向合力计算方程:
注入蒸汽沿油层斜上方渗流,可将速度分解为沿油层X方向的速度和Y方向的速度。根据达西定律,蒸汽在重力、浮力及毛管力的共同作用下,Y方向瞬时渗流速度为:
根据表3的数据,垂直方向渗透率取值600ⅹ10-3µm2,气体粘度取值5 mPa·s,根据式子(A-11)、(A-12)计算得到蒸汽Y方向的瞬时渗流速度为994.40 µm2/s。
S5,计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度。联立式子(A-2)、(A-5)、(A-7)、(A-9),得到X方向合力计算方程:
将(A-2)、(A-5)、(A-5)计算结果代入(A-13),得到:
根据达西定律,油滴在浮力及毛管力的共同作用下,X方向瞬时渗流速度为:
根据表3的数据,水平方向渗透率取值3000ⅹ10-3µm2,气体粘度取值5 mPa·s,根据式子(A-13)、(A-14)计算得到蒸汽X方向的瞬时渗流速度计算方程为:
S6,计算垂向分流系数,并绘制垂向分流系数与井间距离的图版。垂向分流系数为垂向瞬时渗流速度与总瞬时渗流速度的比值,计算方程为(A-16);以为纵坐标,井间距离L为横坐标,绘制垂向分流系数与井间距离L的图版(图5)。
根据表3的数据和A区块的实际情况,为生产井产量,取值34722.22cm3 /s(300m3 /d);为蒸汽注入量,取值34722.22cm3 /s(300m3 /d);为油层厚度,取值50m;L 为注采井间距离,取值20000cm;K 为渗透率,取值;为蒸汽粘度,取值5mPa·s。将上述取值代入(A-17),计算方程变为:
S7,计算注入蒸汽的超覆程度,并绘制超覆程度与井间距离的图版。蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值,用 D表示(表4)。将注采井间距离划分为N块,当蒸汽前缘运移到注采井间 i 块( 0 < i < N) 时,利用此处的径向蒸汽运移量与垂向流量比的乘积和第i-1块的超覆程度迭代,用(A-19)计算出此时的蒸汽超覆程度。并以为纵坐标,井间距离为横坐标,绘制超覆程度与注采井间距离L的图版(图6)。其中,表示第块的垂向分流量;。
表4 B目标区块注采井间不同距离处的垂向分流量和蒸汽超覆程度计算表
S8,定量计算倾斜厚层稠油注蒸汽纵向波及系数。根据S7的图版,用(A-19)计算厚层稠油蒸汽吞吐初期不同分块的蒸汽超覆程度(表4第5列);通过计算每个区块纵向波及无因次厚度(表4第6列);用每个区块的长度(10m)乘以无因次厚度1得到每个区块的面积(表4第7列);用计算得到每个区块的每个区块面积厚度乘积(表4第8列)。
根据表4的数据,根据(A-20)计算得到面积加权后的无因次厚度为36.5(表4第8列求和)/200(表4第7列求和)=0.1825;无因次厚度h取值为1,根据(A-20)计算得到纵向厚度波及系数为18.25%。
上述实施例仅表达了较为优选的实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。
Claims (9)
1.一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法,其特征在于:包括以下步骤:
S1、收集目标区块地质油藏基本参数;
S2、计算目标区块蒸汽微元在油层垂直方向上的受力值;
S3、计算目标区块蒸汽微元在油层水平方向上的受力值;
S4、计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度;
S5、计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度;
S6、计算垂向分流系数,并绘制垂向分流系数与井间距离的图版;
S7、计算注入蒸汽的超覆程度,并绘制超覆程度与井间距离的图版;蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值,用D表示;
S8、定量计算倾斜厚层稠油注蒸汽纵向波及系数;
S81、根据S7的图版,计算得到厚层稠油蒸汽吞吐初期不同分块的蒸汽超覆程度;
S82、计算每个区块纵向波及无因次厚度hi;
S83、用每个区块的长度乘以无因次厚度hi得到每个区块的面积Ai;
S84、用Aihi计算得到每个区块的每个区块面积厚度乘积;
计算得到面积加权后的纵向波及系数Vk。
2.根据权利要求1所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法,其特征在于:所述S1中的基本参数:
是指收集目标区块的基础数据,包括油层厚度、地层倾角、注采井间距离、油的密度、水的密度、气体粘度、界面张力、接触角、孔隙半径、垂直方向渗透率、水平方向渗透率。
3.根据权利要求2所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法,其特征在于:所述S2中计算目标区块蒸汽微元在油层垂直方向上的受力值,包括重力、浮力;
S21、重力的计算
根据蒸汽微元的密度为ρg、蒸汽微元的体积V,再根据重力加速度g,计算得蒸汽微元的重力F重力;
由于地层有倾角θ倾角,以油层水平方向为X轴、油层垂直方向为Y轴,建立直角坐标系,计算蒸汽微元的重力在X方向的分量Fx重力、Y方向的分量FY重力;
S22、浮力的计算
在注入蒸汽突破之前,蒸汽微元受到的浮力F浮力看作是恒定的,方向为竖直向上;
浮力F浮力,根据油相的密度ρo、水相的密度ρw、蒸汽微元的体积V,重力加速度g,计算得到;
由于地层有倾角θ倾角,以油层水平方向为X轴、油层垂直方向为Y轴,建立直角坐标系,计算浮力F浮力在X方向的分量FX浮力、Y方向的分量FY浮力。
4.根据权利要求3所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法,其特征在于:所述S3中计算目标区块蒸汽微元在油层水平方向上的受力值,包括毛管力和驱替压力:
S31、计算驱替压力,即计算注采井间主流线上的注采压力梯度,公式为
其中,Q1为生产井产量;Q2为蒸汽注入量;h为油层厚度;r为任意点到生产井的距离;dp/dr为r处的压力梯度;L为注采井间距离;K为渗透率;μ为原油粘度;
S32、计算毛管力
对于蒸汽,受到毛管力为F毛管力=2×10-5σcosβ接触角/R孔隙
式中,σ为界面张力;θ接角为接触角;R孔隙为孔隙半径;
在水湿情况下,毛细管力是注入水运移的动力,取正值;在油湿情况下,毛管力是注入水运移的阻力,取负值;
由于地层有倾角θ倾角,根据该倾角,计算平行油藏方向的水平毛管力Fx毛管力和垂直油藏方向的垂直毛管力FY毛管力。
5.根据权利要求4所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法,其特征在于:所述步骤S4计算目标区块蒸汽微元在油藏垂直方向的合力和瞬时渗流速度;
S41、计算目标区块蒸汽微元在油藏垂直方向的合力
根据蒸汽微元在Y方向的重力分量FY重力,蒸汽微元在Y方向的重力分量FY浮力,垂直油藏方向的垂直毛细管力FY毛管力,计算得到蒸汽微元在油藏垂直方向的合力FY合力;
S42、计算目标区块蒸汽微元在油藏垂直方向的瞬时渗流速度
注入蒸汽沿油层斜上方渗流,将速度分解为沿油层X方向的速度和Y方向的速度;
蒸汽在重力、浮力及毛管力的共同作用下,根据达西定律,计算得到Y方向瞬时渗流速度,即在油藏垂直方向的瞬时渗流速度VX。
6.根据权利要求5所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法,其特征在于:所述S5计算目标区块蒸汽微元在油藏水平方向的合力和瞬时渗流速度;
S51、计算目标区块蒸汽微元在油藏水平方向的合力
根据蒸汽微元在X方向的重力分量FX重力,蒸汽微元在X方向的重力分量FX浮力,驱替压力,垂直油藏方向的垂直毛细管力FX毛管力,计算得到蒸汽微元在油藏垂直方向的合力FX合力;
S52、计算目标区块蒸汽微元在油藏水平方向的瞬时渗流速度
根据达西定律,油滴在浮力及毛管力的共同作用下,计算得到X方向瞬时渗流速度为VY。
8.根据权利要求7所述一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法,其特征在于:所述S7计算注入蒸汽的超覆程度,并绘制超覆程度与井间距离的图版;
S71、计算注入蒸汽的超覆程度
蒸汽超覆程度定义为累计运移到垂向上的蒸汽量与累计注入蒸汽量的比值,用D表示;
将注采井间距离划分为N块,当蒸汽前缘运移到注采井间i块时,0<i<N,利用此处的径向蒸汽运移量与垂向流量比的乘积和第i-1块的超覆程度迭代,计算出此时的蒸汽超覆程度,计算公式为
Di=(1-Di-1)fi+Di-1
其中,fi表示第i块的垂向分流量;
S72、绘制超覆程度与井间距离的图版
以Di为纵坐标,井间距离为横坐标,绘制超覆程度Di与注采井间距离L的图版。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211385246.6A CN115422859B (zh) | 2022-11-07 | 2022-11-07 | 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211385246.6A CN115422859B (zh) | 2022-11-07 | 2022-11-07 | 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115422859A CN115422859A (zh) | 2022-12-02 |
CN115422859B true CN115422859B (zh) | 2023-01-24 |
Family
ID=84207638
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211385246.6A Active CN115422859B (zh) | 2022-11-07 | 2022-11-07 | 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115422859B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115859530B (zh) * | 2023-01-18 | 2023-08-11 | 西南石油大学 | 一种基于无因次泵效的抽油泵健康状态评价方法 |
CN116341406B (zh) * | 2023-03-03 | 2024-04-23 | 常州大学 | 一种确定倾斜油藏注采干扰时间的计算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4637466A (en) * | 1986-04-03 | 1987-01-20 | Texaco Inc. | Method of improving conformance in steam floods with carboxylate steam foaming agents |
US5214384A (en) * | 1991-07-24 | 1993-05-25 | Mobil Oil Corporation | Method including electrical self potential measurements for detecting multiphase flow in a cased hole |
CN106150466A (zh) * | 2016-08-23 | 2016-11-23 | 中国石油化工股份有限公司 | 凝胶泡沫抑制底水锥进的稠油热采方法 |
CN107435535A (zh) * | 2016-05-26 | 2017-12-05 | 中国石油大学(北京) | 一种采用平面重力驱开采高倾角稠油油藏的方法 |
CN112036033A (zh) * | 2020-08-31 | 2020-12-04 | 中国海洋石油集团有限公司 | 一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2877640C (en) * | 2014-01-13 | 2021-12-14 | John A. Stanecki | Oil recovery with fishbone wells and steam |
CN104318052B (zh) * | 2014-09-16 | 2017-10-27 | 中国石油大学(北京) | 一种低渗透油藏水驱波及系数的评价方法 |
CN108708699B (zh) * | 2018-05-14 | 2019-04-16 | 中国石油大学(华东) | 一种强化sagd蒸汽腔突破低物性储层的超稠油开发方法 |
-
2022
- 2022-11-07 CN CN202211385246.6A patent/CN115422859B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4637466A (en) * | 1986-04-03 | 1987-01-20 | Texaco Inc. | Method of improving conformance in steam floods with carboxylate steam foaming agents |
US5214384A (en) * | 1991-07-24 | 1993-05-25 | Mobil Oil Corporation | Method including electrical self potential measurements for detecting multiphase flow in a cased hole |
CN107435535A (zh) * | 2016-05-26 | 2017-12-05 | 中国石油大学(北京) | 一种采用平面重力驱开采高倾角稠油油藏的方法 |
CN106150466A (zh) * | 2016-08-23 | 2016-11-23 | 中国石油化工股份有限公司 | 凝胶泡沫抑制底水锥进的稠油热采方法 |
CN112036033A (zh) * | 2020-08-31 | 2020-12-04 | 中国海洋石油集团有限公司 | 一种定量确定蒸汽吞吐转蒸汽驱地层温度时机的方法 |
Non-Patent Citations (3)
Title |
---|
"A New Method for Calculating the Inflection Point Temperature of Heavy-Oil Rheology Transforming From Non-Newton Fluid into the Newton Fluid";Dong Liu等;《Journal of Energy Resources Technology》;20211112;第144卷(第8期);第1-11页 * |
"Sequential Multiwell Steam Huff n Puff in Heavy-Oil Development";W.Shuhong等;《2005,SPE/PS-CIM/CHOA International Thermal Operations and Heavy Oil Symposium》;20051101;第1-6页 * |
蒸汽-泡沫复合驱提高稠油采收率机理研究;周阴国;《中国石油和化工标准与质量》;20130115(第02期);第253页 * |
Also Published As
Publication number | Publication date |
---|---|
CN115422859A (zh) | 2022-12-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115422859B (zh) | 一种定量评价厚层稠油注蒸汽吞吐纵向波及系数的方法 | |
Alvarez et al. | Unified model for steady-state foam behavior at high and low foam qualities | |
Li et al. | Generalized scaling approach for spontaneous imbibition: an analytical model | |
Lian et al. | The characteristics of relative permeability curves in naturally fractured carbonate reservoirs | |
CN105822302A (zh) | 一种基于井地电位法的油水分布识别方法 | |
CN104989341A (zh) | 一种确定低渗透油藏有效驱替注采井距的方法 | |
CN109884269A (zh) | 岩心自发渗吸采收率预测方法及*** | |
CN108301825A (zh) | 一种高温高压储层孔隙压力分析方法 | |
Uchida et al. | Sand production modeling of the 2013 Nankai offshore gas production test | |
CN104975827A (zh) | 预测二氧化碳驱油藏指标的物质平衡方法 | |
Christman et al. | Comparison of laboratory-and field-observed CO2 tertiary injectivity | |
Huang et al. | Effects of water invasion law on gas wells in high temperature and high pressure gas reservoir with a large accumulation of water-soluble gas | |
Shahverdi | Characterization of three-phase flow and WAG injection in oil reservoirs | |
Taheriotaghsara et al. | The impact of spatially correlated heterogeneity and adsorption on modified salinity water in carbonates | |
CN112733412B (zh) | 一种水动力作用滑坡运动机制研究的速度等效表征方法 | |
CN105298478A (zh) | 一种用于断裂构造的地层孔隙压力确定方法 | |
CN112836438B (zh) | 一种非均质边水气藏气井见水时间预测方法及装置 | |
Teeuw et al. | Design of a pilot polymer flood in the Marmul Field, Oman | |
Pursley et al. | Borregos field surfactant pilot test | |
Dimri et al. | Reservoir geophysics: Some basic concepts | |
Masalmeh et al. | The importance of special core analysis in modelling remaining oil saturation in carbonate fields | |
Wood et al. | Determining effective residual oil saturation for mixed wettability reservoirs: Endicott Field, Alaska | |
Moradi | Cost-effective and safe oil production from existing and near-future oil fields | |
Temizel et al. | An analysis of scale buildup in seawater injection of waterflooding operations | |
Zhou et al. | Quantifying viscous cross-flow and its impact on tertiary polymer flooding in heterogeneous reservoirs |
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 |