油气界面临界运移速度计算方法、装置、电子设备及介质.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN114065655A(43)申请公布日2022.02.18(21)申请号202010791072.8(22)申请日2020.08.07(71)申请人中国石油化工股份有限公司地址100027北京市朝阳区朝阳门北大街22号申请人中国石油化工股份有限公司石油勘探开发研究院(72)发明人廉培庆段太忠赵华伟赵磊(74)专利代理机构北京思创毕升专利事务所11218代理人孙向民廉莉莉(51)Int.Cl.G06F30/28(2020.01)G01V99/00(2009.01)G06F119/14(2020.01)权利要求书2页说明书11页附图3页(54)发明名称油气界面临界运移速度计算方法、装置、电子设备及介质(57)摘要公开了一种油气界面临界运移速度计算方法、装置、电子设备及介质。该方法可以包括：建立不等径并联毛细管模型，获取临界界面运移速度的计算公式；根据孔隙大小划分多个储层类型；根据计算公式计算每一个储层类型的临界界面运移速度。本发明通过微观孔隙结构定量分析结果，结合注气稳定重力驱可动用孔隙界限计算公式，确定注气稳定重力驱油气界面临界运移速度，将临界界面运移速度与孔喉最小流动界限相结合，在保持油气界面稳定运移的同时，动用绝大多数的孔喉剩余油，提高油藏采收率。CN114065655ACN114065655A权利要求书1/2页1.一种油气界面临界运移速度计算方法，其特征在于，包括：建立不等径并联毛细管模型，获取临界界面运移速度的计算公式；根据孔隙大小划分多个储层类型；根据所述计算公式计算每一个储层类型的临界界面运移速度。2.根据权利要求1所述的油气界面临界运移速度计算方法，其中，所述临界界面运移速度的计算公式为：其中，v界面为临界界面运移速度，r1为最小孔喉半径，r2为最大孔喉半径，ρg、ρo分别为气、油密度，μg、μo分别为气、油黏度，g为重力加速度。3.根据权利要求1所述的油气界面临界运移速度计算方法，其中，多个储层类型包括：微孔储层、小孔储层、中孔储层与大孔储层。4.根据权利要求2所述的油气界面临界运移速度计算方法，其中，根据所述计算公式计算每一个储层类型的临界界面运移速度包括：计算每一个储层类型的平均孔喉分布特征曲线，计算不同孔喉对应的渗透率贡献程度，进而计算最小流动孔喉半径；以最小流动孔喉半径为最小孔喉半径；根据所述计算公式计算每一个储层类型的临界界面运移速度。5.根据权利要求4所述的油气界面临界运移速度计算方法，其中，以累计渗透率贡献程度为99.95％时的孔喉半径为最小流动孔喉半径。6.根据权利要求2所述的油气界面临界运移速度计算方法，其中，还包括：计算孔喉半径下限，使所述最小孔喉半径大于所述孔喉半径下限。7.根据权利要求6所述的油气界面临界运移速度计算方法，其中，通过公式(2)计算孔喉半径下限：其中，r’为孔喉半径下限，ρg、ρo分别为气、油密度，μg、μo分别为气、油黏度，g为重力加速度，v为界面运移速度、θ为油藏倾角。8.一种油气界面临界运移速度计算装置，其特征在于，包括：公式建立模块，建立不等径并联毛细管模型，获取临界界面运移速度的计算公式；储层分类模块，根据孔隙大小划分多个储层类型；计算模块，根据所述计算公式计算每一个储层类型的临界界面运移速度。9.一种电子设备，其特征在于，所述电子设备包括：存储器，存储有可执行指令；处理器，所述处理器运行所述存储器中的所述可执行指令，以实现权利要求1-7中任一2CN114065655A权利要求书2/2页项所述的油气界面临界运移速度计算方法。10.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1-7中任一项所述的油气界面临界运移速度计算方法。3CN114065655A说明书1/11页油气界面临界运移速度计算方法、装置、电子设备及介质技术领域[0001]本发明涉及油气田开发工程领域，更具体地，涉及一种油气界面临界运移速度计算方法、装置、电子设备及介质。背景技术[0002]顶部注气重力驱技术是一种较常见的气驱提高原油采收率技术，在石油开采行业已有较长的应用历史。由于油藏储层往往有较强的非均质和各向异性，以及气体粘度低于原油粘度，常规气驱技术中注入气体容易沿单层或渗流能力较强的方向优先推进。油藏储层顶部注气重力驱技术则有效利用了气体的重力分离作用，将常规水平方向的驱替转变为垂向驱替，以增大气体的波及体积，提高石油采收率。[0003]确定油藏储层内油气界面临界运移速度是顶部注气重力驱技术的关键。如果注气速度过高，影响油气界面运移的稳定性，加快形成驱替气指进和舌进现象，气体过早突破从而降低注气效果。目前，现有技术是利用Dietz模型等方法判断界面稳定性