裂缝性油藏逆向渗吸采收率预测方法

文档序号:9909761阅读:629来源:国知局
裂缝性油藏逆向渗吸采收率预测方法
【技术领域】
[0001 ]本发明属于石油开采技术领域,具体地说,涉及一种预测裂缝性油藏逆向渗吸采 收率的方法,特别是针对存在孤立裂缝的基质岩块。
【背景技术】
[0002] 裂缝性油藏中原油的采收率低、见水早,基质内原油的驱替难度大。对于裂缝性油 藏,特别是低渗透油藏,通常是见高含水后,采取闷井停产一段时间后再生产,依靠毛管力 作用,即渗吸作用来驱替基质内部的原油。闷井时间的长短、渗吸作用下的采收率变化等方 面的研究对裂缝性油藏提高开采率具有重要的意义。
[0003] 现有的基质内渗吸采收率预测方法主要有三种:室内实验、解析解分析和数值解 分析。大多集中在室内小尺度岩样的实验分析和确定采收率与无因次时间t D之间的变化曲 线上(Mattax、Kyte(1962) ;Bourbiaux、Kalaydjian( 1990) ;Zhang( 1996)) eZhang等人 (1996)的实验表明:不同的边界条件下得到的采收率与tD之间的变化曲线,均逼近统一表 达式。Cilet等人(1998)实验分析了逆向渗吸过程中的影响因素,提出了预测采收率的关系
中R为采收率,R〇〇为最终采收率,常数ω D和η取0.0033和0.6,tD为 无量纲时间。还有研究者通过数学计算的方法来研究渗吸过程中的采收率,给出了逆向渗 吸过程中预测采收率的解析表达式、〇 = < ?^;:αωΖ) = adf2, 2 其中 拉为无量纲时间(Zimmerman和Bodvasson(1989) ; Zimmerman(1990) ;Chen(1995); Tavassoil(2005) ;Li(2011)) Ja等人(1997)通过不同亲水类型的岩石实验结果分析,认为 不同亲水类型的岩心,其采收率与无因次时间tD的变化曲线相似,最后均可通过采收率与 无因次时间t D来表示,并给出关系式i?=圪(1 -),其中,R为采收率,1^为最终采收率,α 为常数,当取α = 0.05时适用于大多数实验结果。Ma等人给出的关系式适用于强亲水介质, 且忽略了润湿性的影响。Gupta和Civan(1994)和Cil(1998)引入了油水接触角Θ的余弦项, 代表介质的润湿性影响。当介质亲水时,则cos0 = 1,当忽略润湿性的影响,则α是接触角分 布的函数。Zhou(2002)通过高孔隙度低渗透率的岩心的物理实验,研究了较大流度比范围 的渗吸过程,并给出了新的表达式
其中tD为无量纲时间,t为时间,k 为渗透率,φ为孔隙度,σ为表面张力,Θ为接触角,ynw为流体粘度,L。特征长度。采用数值模 拟的方法来研究渗吸过程(Pooladi-Darvish、Firoozabadi (2000); Behbahani (2006); Bagherinezhad、Pishvaie(2014)),大多是建立在连续介质模型基础上,Behbahani (2006) 应用Edips/MOO模拟软件,采用细网格模拟渗吸过程。上述三种方法各有优势和不足:
[0004] (1)室内物理实验是地下基质岩心在室内条件下物理过程的真实再现。但室内岩 心大多是圆柱形或者长方体等规则形状;鉴于物理过程的缓慢性和实验条件的有限性,大 多实验岩样尺寸为小岩心岩样;在室内实验过程中,借助于测量仪器仪表实现实验数据测 量,一方面测量数据具有局限性,另一方面测量数据容易受到人为干扰和存在测量误差;物 理实验岩心的加工制作时间和物理实验过程耗时比较长,需要实验周期和人力时间较长。
[0005] (2)渗吸作用的解析解是对渗吸过程采用数学方程进行描述,求解控制方程的解 析解。但要求岩心模型为均匀各向同性,难以适用于非均质和各向异性情况,对于岩样内部 发育裂缝的情况更无法获得其解析解;方程是建立在圆柱坐标或者笛卡尔坐标下,解析解 可用于分析圆柱或者长方体等规则形状的渗吸规律,难以适用于不规则形状的岩样模型。
[0006] (3)数值模拟是采用数值法求解控制方程来预测渗吸过程。数值模拟可以适用于 不规则形状的岩样模型;数值模拟可以方便的模拟不同边界条件、不同岩样模型大小、不同 时刻的渗吸状况,可以获得模型中任意网格点处的数值,具有省时边界、获得数据量丰富的 优点。当前数值模拟主要基于连续介质模型。
[0007] 目前,在逆向渗吸研究中,研究者认为,无论哪一种方法获得的采收率R和无因次 时间tD之间的变化关系均可表示成:
[0008]
(1)
[0009] 其中,R为采收率,R〇〇为最终采收率,α为常数,当取α = 〇.05时适用于大多数实验结 果。
[0010] 对于裂缝性油藏而言,不仅存在连通性较好的裂缝网络,还存在裂缝网络连通性 差的孤立裂缝或裂缝盲端,由于具有一定的开度,其毛管力要小于基质孔隙的毛管力,使得 基质岩块表相处强裂缝具有非均质性,进而影响了基岩逆向渗吸采收率变化。因此,基质岩 块内有无裂缝、裂缝的多少和形态对整个岩块的采收率变化影响不同,目前还未见如何预 测发育孤立裂缝的基质内的渗吸采收率的相关报道。

【发明内容】

[0011] 本发明针对现有裂缝性油藏预测逆向渗吸采收率时存在的结果不准确等上述不 足,提供一种裂缝性油藏逆向渗吸采收率预测方法,该方法考虑孤立裂缝、基质岩块切割的 不规则性以及裂缝与基质之间接触的多样性,全面模拟裂缝性油藏的渗吸过程和渗吸规 律,研究预测渗吸作用下裂缝性介质内剩余油变化,实现了对裂缝的精确刻画,能够模拟出 孤立裂缝的影响,预测结果准确。
[0012] 本发明的技术方案是:一种裂缝性油藏逆向渗吸采收率预测方法,含有以下步骤:
[0013] ( - )考虑裂缝对基质岩块切割的不规则性、裂缝与基质之间接触的多样性以及基 质中存在的孤立裂缝,建立裂缝性介质模型,该裂缝性介质模型为基质岩块切割后任一岩 块边缘裂缝首尾相接组成;
[0014](二)采用非结构化网格对裂缝性介质模型进行离散化处理,获得离散裂缝模型;
[0015] (三)根据离散裂缝模型确定基质的渗吸方程和裂缝的渗吸方程,并根据含水饱和 度和流量关系进行基质有限单元之间、裂缝有限单元之间、基质有限单元与裂缝有限单元 之间的耦合,耦合后获得整体的有限元方程;
[0016] (四)数值模拟裂缝性介质模型的渗吸过程与规律,预测渗吸作用下饱和度变化和 采收率变化。
[0017] 在根据本发明上述预测方法中,步骤(一)中,建立裂缝性介质模型的步骤为:
[0018] ( - )搜集反应裂缝发育状况的岩心、地震、成像测井的资料;
[0019] (二)根据岩心统计结果、地震预测结果、成像测井检测结果,统计裂缝发育状况、 裂缝发育密度、裂缝长度、裂缝开度、裂缝倾角的信息,获得大尺度裂缝网络的连通性信息、 裂缝切割基质岩块形状信息以及裂缝与基质的接触关系信息;
[0020] (三)通过CT扫描、扫描电镜、高压压汞和核磁共振分析小尺度裂缝的连通性,确定 基质内部连通性差、或孤立裂缝发育状况、或连通性差和孤立裂缝发育状况;
[0021] (四)根据裂缝网络对基质岩块的切割尺寸和接触关系,以及连通性差、或孤立裂 缝发育状况、或连通性差和孤立裂缝发育状况,确定裂缝性介质模型。
[0022]在根据本发明上述预测方法中,步骤(三)中,基质的渗吸方程表达式为:
[0023]
(2)
[0024] 式中,vm, w为基质内水相流体流速,Pm,。为基质中毛管力,Km为基质的绝对渗透率, km, 为基质中水相流体相对渗透率,km, r。为基质中油相流体相对渗透率,μ?为水相流体年 度,μ。为油相流体粘度;
[0025] 进一步表示为:
[0026]
(3)
[0027] 其中,Sm,w为基岩含水饱和度,基岩内分流量fm,w为:
[0028]
⑷[0029]定义基岩内毛管力扩散系数
[0030] (5):
[0031]
[0032] (6)
[0033] 同时满足质量守恒,则有:
[0034]
(7)
[0035]采用混合有限元法求解基质的渗吸方程式(6)和(7)获得基质有限单元上的渗吸 方程表达式为:
[0036]
当前第1页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1