预览加载中,请您耐心等待几秒...
1/10
2/10
3/10
4/10
5/10
6/10
7/10
8/10
9/10
10/10

亲,该文档总共12页,到这已经超出免费预览范围,如果喜欢就直接下载吧~

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

(19)国家知识产权局(12)发明专利申请(10)申请公布号CN114580134A(43)申请公布日2022.06.03(21)申请号202011385849.7G06F113/08(2020.01)(22)申请日2020.12.01(71)申请人中国石油化工股份有限公司地址257000山东省东营市东营区济南路125号申请人中国石油化工股份有限公司胜利油田分公司勘探开发研究院(72)发明人郭奇杨勇张世明孙业恒黄迎松李伟忠杨盛波(74)专利代理机构济南日新专利代理事务所(普通合伙)37224专利代理师刘亚宁(51)Int.Cl.G06F30/20(2020.01)G06F111/10(2020.01)权利要求书3页说明书7页附图1页(54)发明名称一种疏松砂岩油藏特高含水期渗流单元定量描述方法(57)摘要本发明属于石油开发技术领域,具体涉及一种疏松砂岩油藏渗流单元定量描述方法。本发明方法以室内实验为基础,确定疏松砂岩油藏特高含水期不同物性级别储层参数变化规律及数学模型,形成渗流单元划分的标准,根据高速非达西流特性对高速渗流单元中的非线性渗流区域进行筛选,最终形成一套针对中高渗疏松砂岩油藏的渗流单元评价方法,用以指导特高含水期均衡驱替调整及剩余油挖潜。本发明方法通过渗流单元的划分,可以量化疏松砂岩油藏高含水期各区域开发潜力,为各分区内剩余油挖潜方式提供参考,为老油田提高采收率、改善开发效果起到了重要的经济效益。CN114580134ACN114580134A权利要求书1/3页1.一种疏松砂岩油藏渗流单元定量描述方法,其特征在于,包括以下步骤:建立疏松砂岩油藏物性时变数值模型;疏松砂岩油藏渗流单元定量划分。2.根据权利要求1所述方法,其特征在于,建立能够反映物性变化特征的参数时变数学模型,所述模型由油、水、气三相的连续方程组成;其中,油、水、气三相的连续性方程分别如下:油组分:水组分:气组分:以上方程中,k为油藏绝对渗透率;vo,vg,vw分别为油、气、水相的渗流速度;kro、krg、krw分别为油、气、水相的相对渗透率;μo、μg、μw分别为油、气、水相的粘度;po、pg、pw分别为油、气、水相的压力;So、Sg、Sw分别为油、气、水相的饱和度;ρo、ρg、ρw分别为油、气、水相的密度;g为重力加速度;D为某一基准面算起的深度;Bo、Bg、Bw分别为油、气、水相的体积系数;qo、qg、qw分别为油、水、气相在单位时间、单位体积岩石中注入的质量流量;为孔隙度;Rso为溶解气油比;qvo、qvg、qvw分别为标准条件下油、气、水的体积流量;Rsw为溶解气油比;优选地,基于隐压显饱IMPES方法进行参数时变油藏数值模型各个方程的求解,在每一次时间步长的迭代中,首先计算出各网格块的注入水冲刷孔隙体积倍数,然后依据渗透率、相对渗透率随注入水冲刷孔隙体积倍数的变化规律,通过插值方法求得本次时间步长内各网格块的渗透率及油、水相对渗透率值。3.根据权利要求2所述方法,其特征在于,过水倍数的表达式为:式中,Rw(i)为过水倍数;i为第i个网格数;tD为累积注水时间,s;Qin为网格注水速度,33m/s;Vi为第i个网格的体积,m;为第i个网格孔隙度。4.根据权利要求1所述方法,其特征在于,渗流单元可划分为储量未动用区、高速低潜力区、低速高潜力区、高速高潜力区和低速低潜力区五类分区。5.根据权利要求4所述方法,其特征在于,渗流单元定量划分的方法包括:通过流体强度参数将渗流划分为储量未动用区、高潜力区和低潜力区;若所述区域为高潜力区和低潜力区,再根据渗流单元的渗流速度将所述区域划分为高速区或低速区,即将高潜力区和低2CN114580134A权利要求书2/3页潜力区域划分为高速高潜力区或高速低潜力区,或者划分为低速高潜力区或低速低潜力区。6.根据权利要求5所述方法,其特征在于,采用流体强度划分标准:流场强度在0.6‑1的区域为低潜力区,流场强度在0.1‑0.6的区域为高潜力区,流场强度在0‑0.1的区域为储量未动用区。7.根据权利要求5或6所述方法,其特征在于,采用流场强度划分的具体方法:将通过油藏各方向横截面的累积过水量作为流体强度评价的动态指标,其表达式为:3式中,M为通过油藏横截面流体的累积流量,m;Qx为Δt时间内通过x方向流体的流量,333m;Qy为Δt时间内通过y方向流体的流量,m;Qz为Δt时间内通过z方向流体的流量,m;Dx为x方向截面的长度,m;Dy为y方向截面的长度,m;Dz为z方向截面的长度,m;Δt为单位时间间隔,s;选用对数函数对式对上式计算结果进行归一化,其表达式为:式中,L为流体强度系数;M(i)1为累积过水倍数最小值;M(i)2为累积过水倍数最大值;根据归一化结果,按划分标准进行分区。8.根据权利要求5所述方法,其特征在