CN106547973A - 一种稠油油藏注蒸汽过程汽窜通道的定量描述方法 - Google Patents
一种稠油油藏注蒸汽过程汽窜通道的定量描述方法 Download PDFInfo
- Publication number
- CN106547973A CN106547973A CN201610958641.7A CN201610958641A CN106547973A CN 106547973 A CN106547973 A CN 106547973A CN 201610958641 A CN201610958641 A CN 201610958641A CN 106547973 A CN106547973 A CN 106547973A
- Authority
- CN
- China
- Prior art keywords
- channeling
- well
- injection
- passage
- producing
- 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
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 DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP 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
Abstract
本发明涉及一种稠油油藏注蒸汽过程汽窜通道的定量描述方法,包括以下步骤:对汽驱井组内各生产井汽窜程度分别进行识别,确定各生产井的汽窜类型;根据汽驱井组内各生产井的汽窜类型,以及不同的生产井特征参数,计算该汽驱井组内的汽窜体积。本发明的一种稠油油藏注蒸汽过程汽窜通道的定量描述方法,对汽驱井组所有生产井的窜流程度进行界定和分类,并最终求得注采井组内多向窜流的汽窜通道体积,给出了定量化、可操作的技术方法和实施步骤来计算汽窜通道体积,可以广泛适用于陆上及海上稠油油藏注蒸汽开发研究领域。
Description
技术领域
本发明属于油气田开发领域,尤其涉及一种稠油油藏注蒸汽过程汽窜通道的定量描述方法。
背景技术
稠油油藏经过多周期蒸汽吞吐后,容易出现由于生产井工作制度不同所造成的蒸汽窜流现象。转入蒸汽驱后,随着汽驱程度的不断加深,注采井间单向或多向窜流现象进一步突显,汽窜成为制约蒸汽驱开采效果的主要矛盾。造成汽窜的原因除了储层非均质性差异、油藏出砂等地质因素外,蒸汽与稠油流度比差异悬殊、粘性指进造成注采井间低流动阻力通道形成以及蒸汽超覆等多方面原因均可加剧蒸汽在储层中的突进。除此之外,由于注采参数不当,开发设计不合理等开采因素造成的汽窜现象也不在少数。
对于汽窜现象的识别和汽窜通道的描述,目前普遍采用的技术和方法有示踪剂判别法、试井资料分析法和基于测井资料处理与解释成果的专家系统模糊识别法等。而对于汽窜通道的计算,除了基于圆柱体地层体积模型的传统设计方法外,余倩等提出了针对水驱油藏窜流通道并基于近井孔隙体积和远井双曲线舌进孔隙体积的设计方法;郑家朋等引入汽窜时间和热侵体积的概念定量描述汽窜发生的快慢及强度;刘丽娜在考虑汽窜通道为蒸汽凝析特征的基础上,推导了非活塞驱替模型中的汽窜体积模型,并用该方法确定井间窜通孔隙体积。
但是,上述方法均暴露出一些不足:1)缺乏对产生汽窜原因的地质因素和生产动态因素的综合考虑;2)缺乏对汽驱井组内生产井汽窜程度的有效判定;3)缺乏注采井组内不同汽窜程度生产井汽窜通道的差异化计算。
发明内容
针对上述问题,本发明的目的是提供一种稠油油藏注蒸汽过程汽窜通道的定量描述方法,在对汽窜现象分析的基础上,优选汽窜判别参数,引入汽窜综合评价指标,利用层次分析法对汽窜井进行判别,并在此基础上利用各生产井的差异化参数对汽窜通道进行定量描述。
为实现上述目的,本发明采取以下技术方案:一种稠油油藏注蒸汽过程汽窜通道的定量描述方法,包括以下步骤:
1)对汽驱井组内各生产井汽窜程度分别进行识别,确定各生产井的汽窜类型:
①采集汽驱井组内各生产井的基本油藏参数和动态开发指标作为汽窜通道识别过程的基本参数,汽窜通道识别过程的基本参数包括渗透率、孔隙度、有效厚度、注入温度、井口出油温度、蒸汽注入量、采液速度和采油速度;
②对各生产井的汽窜通道识别过程的基本参数进行预处理,得到各生产井的各汽窜通道判别因素指标,将各生产井的各汽窜通道判别因素指标进行极差正规化处理,得到各生产井最终的各汽窜通道判别因素指标;
③利用层次分析法,分别计算各生产井的各汽窜通道判别因素指标的权重,得到各生产井有效的各汽窜通道判别因素指标权重;
④根据各生产井有效的各汽窜通道判别因素指标权重,分别计算各生产井的综合评价指标,其中,某一口生产井的综合评价指标计算公式为:
式中,CIP为该生产井的综合评价指标,CIP的取值为0到1;i为汽窜通道判别因素指标编号;wi为该生产井的第i种汽窜通道判别因素指标权重;R′i为该生产井第i种汽窜通道判别因素指标极差正规化处理后的参数值;
⑤根据各生产井的综合评价指标CIP,参考基尼系数取值范围设定综合评价指标CIP的取值,分别判断各生产井的汽窜程度:若CIP>0.5,则该生产井完全汽窜;若0.2≤CIP≤0.5,则该生产井未完全汽窜;CIP<0.2,则该井未汽窜;将各生产井均分为CIP>0.5的完全汽窜井以及CIP≤0.5的未完全汽窜井;
2)根据步骤1)确定的汽驱井组内各生产井的汽窜类型,以及不同的生产井特征参数,计算该汽驱井组内的汽窜体积:
I、确定汽驱井组的注入井附近汽窜半径;
II、依据不同的生产井特征参数,分别对相邻两完全汽窜井与注入井的井间夹角进行劈分,并根据注入井附近汽窜半径,确定不同的劈分点,从而计算出不同的汽窜通道体积;生产井特征参数包括渗透率、含水率、采液强度、井口出油温度以及累积产油量,依据不同生产井间渗透率、含水率、采液强度、井口出油温度以及累积产油量的差异分别对汽窜通道体积进行计算;
III、将依据不同的生产井特征参数分别计算得到的不同的汽窜通道体积求平均,得到该汽驱井组内的汽窜体积。
所述步骤②中,通过对各生产井的汽窜通道识别过程基本参数的组合和处理,每个生产井得到4个汽窜通道判别因素指标,分别是:
式中,R1为注采井渗透率差异表征参数;KI和KP分别为注入井和生产井处的平均渗透率;
式中,R2为注采井温度差异表征参数;TI和TP分别为注入井和生产井处的井口温度;
式中,R3为生产井累积水油比;Nw和No分别为累积产水量和累积产油量;
式中,R4为产液强度,表征生产井的生产能力;ql为产液速度;h为各层的有效厚度;φ为各层的孔隙度;对于多层合采而言,qL为总产液速度,hφ为各层hφ之和;
其中,R1值越大,注入蒸汽越容易沿该注采方向突进,即蒸汽窜流现象更易发生;R2值越大,表明注入井和生产井温度差异越小,也表明蒸汽前缘在该生产井方向上运移速度较快;R3用来评价生产井在产油过程中的产水强度,反映蒸汽的主要注入方向。
所述步骤②中极差正规化处理的公式为:
式中,R′ij为第j口井的第i个汽窜通道判别因素指标极差正规化处理后的参数值;Rimax和Rimin分别为第i个汽窜通道判别因素指标的最大值和最小值。
所述步骤③中利用层次分析法计算各汽窜通道判别因素指标的权重,具体包括以下步骤:
a、利用1~9标度法,依次将各汽窜通道判别因素指标两两进行比对,初步确定各汽窜通道判别因素指标的对比标度权重,得到判别矩阵;
b、求解判别矩阵的最大特征值和对应的特征向量,得到各汽窜通道判别因素指标的相对重要性权重;
c、采用一致性检验指标检验得到的各汽窜通道判别因素指标相对重要性权重的一致性,得到有效的各汽窜通道判别因素指标权重。
所述步骤I中,利用汽驱井组内注入井蒸汽吞吐结束时近井地带原油亏空体积计算亏空半径,将亏空半径作为汽驱井组的注入井附近汽窜半径;亏空半径的计算公式为:
式中,r为注入井蒸汽吞吐结束时近井地带原油亏空体积半径,即注入井附近汽窜半径;Np为注入井吞吐过程中的累积产油量;h为注入井点所在油层的有效厚度;Soi为初始含油饱和度;Soc为蒸汽吞吐结束时剩余油饱和度;φ为各层的孔隙度。
本发明由于采取以上技术方案,其具有以下优点:1、本发明在筛选及预处理汽窜评价指标参数的基础上,利用层次分析法得到汽窜判别综合评价指标,对汽驱井组所有生产井的窜流程度进行界定和分类。2、本发明基于蒸汽吞吐阶段注汽井近井地带原油亏空半径,考虑相邻汽窜生产井特征参数差异,计算注采井间多向汽窜通道面积,并最终求得注采井组内多向窜流的汽窜通道体积。3、本发明给出了定量化、可操作的技术方法和实施步骤,来计算汽窜通道体积。4、本发明适用于陆上及海上稠油油藏注蒸汽开发研究领域。
附图说明
图1是汽驱井组多向窜流汽窜体积的确定示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明的一种稠油油藏注蒸汽过程汽窜通道的定量描述方法,包括以下步骤:
1)对汽驱井组内各生产井汽窜程度分别进行识别,确定各生产井的汽窜类型。具体包括以下步骤:
①采集汽驱井组内各生产井的基本油藏参数和动态开发指标作为汽窜通道识别过程的基本参数。
依据可能造成窜流现象的地质因素,综合矿场实际获得的生产动态参数,选取部分地质和生产参数作为汽窜通道识别过程的基本参数,包括渗透率、孔隙度、有效厚度、注入温度、井口出油温度、蒸汽注入量(冷水当量)、采液速度和采油速度。
②对各生产井的汽窜通道识别过程基本参数进行预处理,得到各生产井的各汽窜通道判别因素指标。
由于考虑的汽窜通道识别过程基本参数较多,而各基本参数在反应不同油藏地质特征和生产动态的同时又各有侧重,并可能出现重叠评价或差异评价,因此需要对基本参数进行预处理,得到能够全面、简洁反映对汽窜现象影响的各汽窜通道判别因素指标。
通过对汽窜通道识别过程基本参数的组合和处理,每个生产井可以得到的4个汽窜通道判别因素指标:
或者
式中,R1为注采井渗透率差异表征参数;KI和KP分别为注入井和生产井处的平均渗透率。R1值越大,注入蒸汽越容易沿该注采方向突进,即蒸汽窜流现象更易发生。
式中,R2为注采井温度差异表征参数;TI和Tp分别为注入井和生产井处的井口温度。R2值越大,表明注入井和生产井温度差异越小,也表明蒸汽前缘在该生产井方向上运移速度较快。
式中,R3为生产井累积水油比;Nw和No分别为累积产水量和累积产油量。R3用来评价生产井在产油过程中的产水强度。在蒸汽注入过程中,生产井的产水量主要来自注入井的注入蒸汽,因此R3在一定程度上可以反映蒸汽的主要注入方向。
式中,R4为产液强度;ql为产液速度;h为各层的有效厚度;φ为各层的孔隙度。对于多层合采而言,ql为总产液速度,hφ为各层hφ之和。R4表征了生产井的生产能力。
由于4个汽窜通道判别因素指标的量纲各不相同,绝对值之差较大,因此将各汽窜通道判别因素指标进行极差正规化处理,得到最终的汽窜通道判别因素指标。极差正规化处理如下式所示:
式中,R′ij为第j口井的第i个汽窜通道判别因素指标极差正规化处理后的参数值;Rimax和Rimin分别为第i个汽窜通道判别因素指标的最大值和最小值。一个汽驱井组有多口生产井,每一生产口井均可得到上述四个判别因素指标,因此,每一个判别因素指标在汽驱井组内均有一个最大值和最小值,进行极差正规化处理后,可以得到每一口井的各判别因素指标的处理值,该处理过程的主要目的是将实际值转化为0~1内的值,使得所有判别因素指标值均位于0~1之间,便于后续对比。
③利用层次分析法,分别计算各生产井的各汽窜通道判别因素指标的权重。
为了将4个汽窜通道判别因素指标综合考虑在内,引入综合评价指标CIP,并利用层次分析法(AHP,Analytic Hierarchy Process)确定各汽窜通道判别因素指标的权重。层次分析法是一种将定性分析和定量计算相结合的系统分析方法,它把复杂的问题分为若干层次,依据一定的判别标准对每一层次各元素的相对重要性进行判别并形成矩阵形式,通过求解该矩阵的最大特征根及对应的特征向量,确定每一层次各元素相对重要性的权重,进而得到综合判别指标。
利用层次分析法计算各汽窜通道判别因素指标的权重,具体包括以下步骤:
a、利用1~9标度法,依次将渗透率差异R1、温度差异R2、累积水油比R3和采液强度R4两两进行比对,初步确定渗透率差异R1、温度差异R2、累积水油比R3和采液强度R4的对比标度权重,得到判别矩阵。
1~9标度法如下表1所示:
表1 1~9标度法
最终得到的判别矩阵如下:
式中,A为判别矩阵;aij为指标Ri与指标Rj的对比标度,i,j=1,2,3,4。
b、求解判别矩阵的最大特征值和对应的特征向量,得到各汽窜通道判别因素指标的相对重要性权重。如下式所示:
AW=λmaxW (7)
W=(w1,w2,w3,w4) (8)
式中,λmax为判别矩阵的最大特征值;W为判别矩阵对应的特征向量;wi为指标Ri的相对重要性权重。
c、采用一致性检验指标检验得到的各汽窜通道判别因素指标相对重要性权重的一致性,得到有效的各汽窜通道判别因素指标权重。
由于判别矩阵存在重复信息,因此层次分析法提供了一致性检验指标:
式中,CR为一致性检验指标;IC为一致性指标;n为判别矩阵A的维数,判别矩阵的维数即为判别矩阵的秩,矩阵的秩意指矩阵中非零子式的最高阶数;IR为一定样本容量(500)时的平均随机一致性指标。
若CR>0.1,则需重新进行对比以减小不一致性;若CR<0.1,则一致性结果可以接受,各汽窜通道判别因素指标权重有效。
④根据各生产井有效的各汽窜通道判别因素指标权重,分别计算各生产井的综合评价指标。综合评价指标计算公式为:
式中,CIP为综合评价指标,CIP的取值为0到1;i为汽窜通道判别因素指标编号。
⑤根据各生产井的综合评价指标CIP,分别判断各生产井的汽窜程度,确定各生产井的汽窜类型。
为了利用综合评价指标来判别汽窜现象的发生与否,需对综合评价指标CIP的取值范围进行确定。为了将生产井划分为完全汽窜井和不完全汽窜井,参考基尼系数的取值范围来设定综合评价指标CIP的取值,即若CIP>0.5,则该生产井完全汽窜;若0.2≤CIP≤0.5,则该生产井未完全汽窜;CIP<0.2,则该井未汽窜。在实际应用过程中,可将所有生产井均分为CIP>0.5的完全汽窜井以及CIP≤0.5的未完全汽窜井。
2)如图1所示,根据步骤1)确定的汽驱井组内各生产井的汽窜类型,以及不同的生产井特征参数,计算汽驱井组内的汽窜体积。
图1中实心三角形表示注入井,实心圆表示完全汽窜井,空心圆表示未完全汽窜井,其中,设注入井O的坐标为(xo,yo),其中两口相邻完全汽窜井A和B的坐标分别为(xa,ya)和(xb,yb),其中一口未完全汽窜井C的坐标为(xc,yc)。
计算汽驱井组的汽窜通道体积,具体包括以下步骤:
I、确定注入井附近汽窜半径。
蒸汽驱井组内的中心注入井首先经过蒸汽吞吐阶段,而其之所以选定为注入井,除了井网形式的限制外,在蒸汽驱阶段转为注入井的生产井往往与周围生产井的连通性好,产液量大,综合含水率高,近井周围控制储量采出程度高等。因此,利用注入井蒸汽吞吐结束时近井地带原油亏空体积计算亏空半径,将其作为注入井附近汽窜半径:
πr2hφ(Soi-Soc)=Np (13)
式中,r为注入井蒸汽吞吐结束时近井地带原油亏空体积半径,即注入井附近汽窜半径;Np为注入井吞吐过程中的累积产油量;h为注入井点所在油层的有效厚度;Soi为初始含油饱和度;Soc为蒸汽吞吐结束时剩余油饱和度。
II、分别依据不同的生产井特征参数,分别对相邻两完全汽窜井与注入井的井间夹角进行劈分,并根据注入井附近汽窜半径r,确定不同的劈分点A’B’,从而计算出不同的汽窜通道体积。生产井特征参数包括渗透率、含水率、采液强度、井口出油温度以及累积产油量,依据不同生产井间渗透率、含水率、采液强度、井口出油温度以及累积产油量的差异分别对汽窜通道体积进行计算。
利用相邻汽窜井间的特征参数差异对汽窜井与注入井井间夹角进行劈分的公式为:
式中,∠AOB为完全汽窜井A和B与注入井O的井间夹角;LOA和LOB分别表示注入井O与完全汽窜井A和B之间的距离;LAB表示完全汽窜井A和B之间的距离;∠AOA′B′k表示特征参数为ck时完全汽窜井A和B与劈分点A’B’的夹角;ckA和ckB分别表示完全汽窜井A和B的特征参数ck的值;表示特征参数为ck时完全汽窜井A和B之间的窜流面积。
从而可以求得汽驱井组内窜流总面积为:
式中,Sck为按照汽窜井特征参数ck计算得到的汽驱井组总窜流面积;Ln为注入井与第n口完全汽窜井间的距离;θm为注入井与第m口完全汽窜井间夹角的劈分角;f为汽驱井组中完全汽窜井的数目;为汽驱井组平均有效厚度;为汽驱井组平均孔隙度;Vck为按照汽窜井特征参数ck计算得到的总窜流体积。
III、将计算得到的不同特征参数的汽窜通道体积求平均,得到该汽驱井组内的汽窜体积:
式中,V为汽驱井组内的汽窜体积;k为特征参数的个数。
下面选取某个稠油油藏汽驱井组为例,具体说明本发明的一种稠油油藏注蒸汽过程汽窜通道的定量描述方法。具体包括以下步骤:
1)对汽驱井组内各生产井汽窜程度分别进行识别。具体包括以下步骤:
①采集井组内各生产井的基本油藏参数和动态开发指标,数据如下表2所示:
表2汽驱井组基础参数
②对各生产井的汽窜通道识别过程基本参数进行预处理,得到各生产井的各汽窜通道判别因素指标。
③分别计算各生产井的各汽窜通道判别因素指标的权重。其中,得到的判别矩阵为:
求解判别矩阵的最大特征值和对应的特征向量,分别为:
λmax=4.06,W=(0.15,0.37,0.37,0.11),CR=0.02<0.1 (22)
④根据各生产井有效的各汽窜通道判别因素指标权重,分别计算各生产井的综合评价指标。
⑤根据各生产井的综合评价指标CIP,分别判断各生产井的汽窜程度。结果如下表3所示:
表3汽驱井组汽窜综合评价参数表
井号 | X1 | X2 | X3 | X4 | X5 | X6 | X7 |
CIP | 0.226 | 0.481 | 0.771 | 0.609 | 0.723 | 0.197 | 0.250 |
汽窜类型 | 未窜 | 未窜 | 窜 | 窜 | 窜 | 未窜 | 未窜 |
可以得出,X3、X4和X5三口井的综合评价指标CIP值均大于0.5,根据判别标准,这三口井视为完全汽窜井,其他生产井则视为未完全汽窜井。
2)根据步骤1)确定的汽驱井组内各生产井的汽窜类型,以及不同的生产井特征参数,计算汽驱井组的汽窜通道体积。
I、确定注入井附近汽窜半径。
II、分别依据生产井的渗透率、含水率、采液强度、井口出油温度以及累积产油量对汽窜通道体积进行计算。计算结果如下表4所示:
表4汽驱井组汽窜综合评价参数表
特征参数 | 累产油量 | 采液强度 | 出油温度 | 渗透率 | 含水率 | 平均 |
汽窜通道体积 | 2178 | 2254 | 2299 | 2481 | 2520 | 2347 |
III、将计算得到的不同特征参数的汽窜通道体积求平均,得到该汽驱井组内的汽窜体积。
上述各实施例仅用于说明本发明,其中各部件的结构、设置位置及其连接方式等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (7)
1.一种稠油油藏注蒸汽过程汽窜通道的定量描述方法,包括以下步骤:
1)对汽驱井组内各生产井汽窜程度分别进行识别,确定各生产井的汽窜类型:
①采集汽驱井组内各生产井的基本油藏参数和动态开发指标作为汽窜通道识别过程的基本参数,汽窜通道识别过程的基本参数包括渗透率、孔隙度、有效厚度、注入温度、井口出油温度、蒸汽注入量、采液速度和采油速度;
②对各生产井的汽窜通道识别过程的基本参数进行预处理,得到各生产井的各汽窜通道判别因素指标,将各生产井的各汽窜通道判别因素指标进行极差正规化处理,得到各生产井最终的各汽窜通道判别因素指标;
③利用层次分析法,分别计算各生产井的各汽窜通道判别因素指标的权重,得到各生产井有效的各汽窜通道判别因素指标权重;
④根据各生产井有效的各汽窜通道判别因素指标权重,分别计算各生产井的综合评价指标,其中,某一口生产井的综合评价指标计算公式为:
式中,CIP为该生产井的综合评价指标,CIP的取值为0到1;i为汽窜通道判别因素指标编号;wi为该生产井的第i种汽窜通道判别因素指标权重;R′i为该生产井第i种汽窜通道判别因素指标极差正规化处理后的参数值;
⑤根据各生产井的综合评价指标CIP,参考基尼系数取值范围设定综合评价指标CIP的取值,分别判断各生产井的汽窜程度:若CIP>0.5,则该生产井完全汽窜;若0.2≤CIP≤0.5,则该生产井未完全汽窜;CIP<0.2,则该井未汽窜;将各生产井均分为CIP>0.5的完全汽窜井以及CIP≤0.5的未完全汽窜井;
2)根据步骤1)确定的汽驱井组内各生产井的汽窜类型,以及不同的生产井特征参数,计算该汽驱井组内的汽窜体积:
I、确定汽驱井组的注入井附近汽窜半径;
II、依据不同的生产井特征参数,分别对相邻两完全汽窜井与注入井的井间夹角进行劈分,并根据注入井附近汽窜半径,确定不同的劈分点,从而计算出不同的汽窜通道体积;生产井特征参数包括渗透率、含水率、采液强度、井口出油温度以及累积产油量,依据不同生产井间渗透率、含水率、采液强度、井口出油温度以及累积产油量的差异分别对汽窜通道体积进行计算;
III、将依据不同的生产井特征参数分别计算得到的不同的汽窜通道体积求平均,得到该汽驱井组内的汽窜体积。
2.如权利要求1所述的一种稠油油藏注蒸汽过程汽窜通道的定量描述方法,其特征在于,所述步骤②中,通过对各生产井的汽窜通道识别过程基本参数的组合和处理,每个生产井得到4个汽窜通道判别因素指标,分别是:
式中,R1为注采井渗透率差异表征参数;KI和KP分别为注入井和生产井处的平均渗透率;
式中,R2为注采井温度差异表征参数;TI和TP分别为注入井和生产井处的井口温度;
式中,R3为生产井累积水油比;Nw和No分别为累积产水量和累积产油量;
式中,R4为产液强度,表征生产井的生产能力;ql为产液速度;h为各层的有效厚度;φ为各层的孔隙度;对于多层合采而言,ql为总产液速度,hφ为各层hφ之和;
其中,R1值越大,注入蒸汽越容易沿该注采方向突进,即蒸汽窜流现象更易发生;R2值越大,表明注入井和生产井温度差异越小,也表明蒸汽前缘在该生产井方向上运移速度较快;R3用来评价生产井在产油过程中的产水强度,反映蒸汽的主要注入方向。
3.如权利要求1或2所述的一种稠油油藏注蒸汽过程汽窜通道的定量描述方法,其特征在于,所述步骤②中极差正规化处理的公式为:
式中,R′ij为第j口井的第i个汽窜通道判别因素指标极差正规化处理后的参数值;Rimax和Rimin分别为第i个汽窜通道判别因素指标的最大值和最小值。
4.如权利要求1或2所述的一种稠油油藏注蒸汽过程汽窜通道的定量描述方法,其特征在于,所述步骤③中利用层次分析法计算各汽窜通道判别因素指标的权重,具体包括以下步骤:
a、利用1~9标度法,依次将各汽窜通道判别因素指标两两进行比对,初步确定各汽窜通道判别因素指标的对比标度权重,得到判别矩阵;
b、求解判别矩阵的最大特征值和对应的特征向量,得到各汽窜通道判别因素指标的相对重要性权重;
c、采用一致性检验指标检验得到的各汽窜通道判别因素指标相对重要性权重的一致性,得到有效的各汽窜通道判别因素指标权重。
5.如权利要求3所述的一种稠油油藏注蒸汽过程汽窜通道的定量描述方法,其特征在于,所述步骤③中利用层次分析法计算各汽窜通道判别因素指标的权重,具体包括以下步骤:
a、利用1~9标度法,依次将各汽窜通道判别因素指标两两进行比对,初步确定各汽窜通道判别因素指标的对比标度权重,得到判别矩阵;
b、求解判别矩阵的最大特征值和对应的特征向量,得到各汽窜通道判别因素指标的相对重要性权重;
c、采用一致性检验指标检验得到的各汽窜通道判别因素指标相对重要性权重的一致性,得到有效的各汽窜通道判别因素指标权重。
6.如权利要求1或2或5所述的一种稠油油藏注蒸汽过程汽窜通道的定量描述方法,其特征在于,所述步骤I中,利用汽驱井组内注入井蒸汽吞吐结束时近井地带原油亏空体积计算亏空半径,将亏空半径作为汽驱井组的注入井附近汽窜半径;亏空半径的计算公式为:
式中,r为注入井蒸汽吞吐结束时近井地带原油亏空体积半径,即注入井附近汽窜半径;Np为注入井吞吐过程中的累积产油量;h为注入井点所在油层的有效厚度;Soi为初始含油饱和度;Soc为蒸汽吞吐结束时剩余油饱和度;φ为各层的孔隙度。
7.如权利要求3所述的一种稠油油藏注蒸汽过程汽窜通道的定量描述方法,其特征在于,所述步骤I中,利用汽驱井组内注入井蒸汽吞吐结束时近井地带原油亏空体积计算亏空半径,将亏空半径作为汽驱井组的注入井附近汽窜半径;亏空半径的计算公式为:
式中,r为注入井蒸汽吞吐结束时近井地带原油亏空体积半径,即注入井附近汽窜半径;Np为注入井吞吐过程中的累积产油量;h为注入井点所在油层的有效厚度;Soi为初始含油饱和度;Soc为蒸汽吞吐结束时剩余油饱和度;φ为各层的孔隙度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610958641.7A CN106547973B (zh) | 2016-11-02 | 2016-11-02 | 一种稠油油藏注蒸汽过程汽窜通道的定量描述方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610958641.7A CN106547973B (zh) | 2016-11-02 | 2016-11-02 | 一种稠油油藏注蒸汽过程汽窜通道的定量描述方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106547973A true CN106547973A (zh) | 2017-03-29 |
CN106547973B CN106547973B (zh) | 2019-12-10 |
Family
ID=58393527
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610958641.7A Active CN106547973B (zh) | 2016-11-02 | 2016-11-02 | 一种稠油油藏注蒸汽过程汽窜通道的定量描述方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106547973B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106909757A (zh) * | 2017-04-01 | 2017-06-30 | 中海石油(中国)有限公司 | 一种低渗透油藏超前注水地层合理压力水平确定方法 |
CN107832481A (zh) * | 2017-08-03 | 2018-03-23 | 中国石油化工股份有限公司 | 一种用于稠油油藏组合蒸汽吞吐的分区方法 |
CN109281662A (zh) * | 2018-08-16 | 2019-01-29 | 中国石油天然气股份有限公司 | 水平井不同节段位置储层相对可压性指数的确定方法 |
CN112796745A (zh) * | 2021-02-18 | 2021-05-14 | 中海油田服务股份有限公司 | 一种蒸汽驱油藏生产动态数据的确定方法和装置 |
CN114198073A (zh) * | 2020-08-31 | 2022-03-18 | 中国石油化工股份有限公司 | 蒸汽吞吐汽窜水平井提高采收率方法 |
CN115341885A (zh) * | 2022-08-29 | 2022-11-15 | 中国石油大学(北京) | 蒸汽驱窜流通道体积的确定方法和装置 |
CN115370334A (zh) * | 2021-05-20 | 2022-11-22 | 中国石油天然气股份有限公司 | 降低油藏汽窜的注汽调整方法 |
US11933149B2 (en) | 2020-06-18 | 2024-03-19 | Cenovus Energy Inc. | Fluid flow control in a hydrocarbon recovery operation |
CN115370334B (zh) * | 2021-05-20 | 2024-05-14 | 中国石油天然气股份有限公司 | 降低油藏汽窜的注汽调整方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104314540A (zh) * | 2014-08-19 | 2015-01-28 | 中国石油天然气股份有限公司 | 一种注蒸汽油藏防治汽窜方法 |
CN104879106A (zh) * | 2014-02-27 | 2015-09-02 | 中国石油化工股份有限公司 | 化学蒸汽驱提高稠油油藏采收率的方法 |
-
2016
- 2016-11-02 CN CN201610958641.7A patent/CN106547973B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104879106A (zh) * | 2014-02-27 | 2015-09-02 | 中国石油化工股份有限公司 | 化学蒸汽驱提高稠油油藏采收率的方法 |
CN104314540A (zh) * | 2014-08-19 | 2015-01-28 | 中国石油天然气股份有限公司 | 一种注蒸汽油藏防治汽窜方法 |
Non-Patent Citations (4)
Title |
---|
CHUAN LU,ET AL: "A new profile control design based on quantitative identification of steam breakthrough channel in heavy oil reservoirs", 《JOURNAL OF PETROLEUM EXPLORATION AND PRODUCTION TECHNOLOGY》 * |
卢川等: "浅薄层稠油油藏氮气泡沫调驱适应性研究", 《油气地质与采收率》 * |
郑强等: "油藏注水开发后期窜流通道定量识别方法", 《石油钻探技术》 * |
郑强等: "蒸汽驱后汽窜通道定量描述", 《中国科学》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106909757B (zh) * | 2017-04-01 | 2020-01-24 | 中海石油(中国)有限公司 | 一种低渗透油藏超前注水地层合理压力水平确定方法 |
CN106909757A (zh) * | 2017-04-01 | 2017-06-30 | 中海石油(中国)有限公司 | 一种低渗透油藏超前注水地层合理压力水平确定方法 |
CN107832481A (zh) * | 2017-08-03 | 2018-03-23 | 中国石油化工股份有限公司 | 一种用于稠油油藏组合蒸汽吞吐的分区方法 |
CN107832481B (zh) * | 2017-08-03 | 2021-07-06 | 中国石油化工股份有限公司 | 一种用于稠油油藏组合蒸汽吞吐的分区方法 |
CN109281662A (zh) * | 2018-08-16 | 2019-01-29 | 中国石油天然气股份有限公司 | 水平井不同节段位置储层相对可压性指数的确定方法 |
CN109281662B (zh) * | 2018-08-16 | 2022-03-29 | 中国石油天然气股份有限公司 | 水平井不同节段位置储层相对可压性指数的确定方法 |
US11933149B2 (en) | 2020-06-18 | 2024-03-19 | Cenovus Energy Inc. | Fluid flow control in a hydrocarbon recovery operation |
US11939847B2 (en) | 2020-06-18 | 2024-03-26 | Cenovus Energy Inc. | Fluid flow control in a hydrocarbon recovery operation |
CN114198073A (zh) * | 2020-08-31 | 2022-03-18 | 中国石油化工股份有限公司 | 蒸汽吞吐汽窜水平井提高采收率方法 |
CN112796745B (zh) * | 2021-02-18 | 2022-10-11 | 中海油田服务股份有限公司 | 一种蒸汽驱油藏生产动态数据的确定方法和装置 |
CN112796745A (zh) * | 2021-02-18 | 2021-05-14 | 中海油田服务股份有限公司 | 一种蒸汽驱油藏生产动态数据的确定方法和装置 |
CN115370334A (zh) * | 2021-05-20 | 2022-11-22 | 中国石油天然气股份有限公司 | 降低油藏汽窜的注汽调整方法 |
CN115370334B (zh) * | 2021-05-20 | 2024-05-14 | 中国石油天然气股份有限公司 | 降低油藏汽窜的注汽调整方法 |
CN115341885A (zh) * | 2022-08-29 | 2022-11-15 | 中国石油大学(北京) | 蒸汽驱窜流通道体积的确定方法和装置 |
CN115341885B (zh) * | 2022-08-29 | 2023-07-18 | 中国石油大学(北京) | 蒸汽驱窜流通道体积的确定方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106547973B (zh) | 2019-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106547973A (zh) | 一种稠油油藏注蒸汽过程汽窜通道的定量描述方法 | |
CN105740563B (zh) | 一种成熟油田二次开发之优势通道识别方法 | |
CN107829718A (zh) | 基于均衡水驱理念的油藏井网及注采方案优化设计方法 | |
CN104989341B (zh) | 一种确定低渗透油藏有效驱替注采井距的方法 | |
CN110288258A (zh) | 一种高含水油藏剩余油挖潜方法 | |
CN104965979A (zh) | 一种致密砂岩有效储层识别方法 | |
CN106150477A (zh) | 一种确定单井控制储量的方法 | |
CN109614736B (zh) | 一种海上稠油油田蒸汽吞吐开发定向井与直井产能倍数的确定方法 | |
CN105735961B (zh) | 一种低渗致密储层分层压裂级数优选方法 | |
CN105868508A (zh) | 一种基于气测录井信息的产能定量预测方法 | |
CN104343445B (zh) | 识别气层与水层、高气油比储层及含气饱和度的方法 | |
CN105005712B (zh) | 灰岩含水层富水性评价方法 | |
CN104989357A (zh) | 一种选择致密油气藏重复压裂目标井的方法 | |
CN108319738A (zh) | 一种页岩气井产量预测方法 | |
CN110608023B (zh) | 稠油分层注汽的适应性界限分析评价方法 | |
Liu et al. | Well type and pattern optimization method based on fine numerical simulation in coal-bed methane reservoir | |
CN113743023B (zh) | 一种二氧化碳驱油气窜通道分级表征方法 | |
CN104632150B (zh) | 海上油田不同井组合理产液量确定的方法 | |
Hou et al. | Numerical simulation and evaluation of the fracturing and tight gas production with a new dimensionless fracture conductivity (f cd) model | |
Xu et al. | A fractal physics-based data-driven model for water-flooding reservoir (FlowNet-fractal) | |
CN108661628A (zh) | 一种基于参数优选的工程甜点定量评价方法 | |
CN106779229A (zh) | 一种高含水水平井堵水潜力预测方法 | |
CN104727789A (zh) | 中高渗砂岩油藏水驱波及系数及过水倍数动态描述方法 | |
Jayne et al. | Permeability correlation structure of the Columbia River Plateau and implications for fluid system architecture in continental large igneous provinces | |
CN107798211A (zh) | 地下烃源岩层中有机质的有机酸生成量的确定方法和装置 |
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 |