针对强非均质地层的油藏数值模拟方法技术

本发明专利技术提供一种针对强非均质地层的油藏数值模拟方法，该针对强非均质地层的油藏数值模拟方法包括:步骤1，确定强间断存在区域并划分网络；步骤2，计算梯度强度因子并得到局部解析解；步骤3，应用局部解析解计算网格间等效渗透率；步骤4，构建跨界面相流量计算公式；步骤5，求解待求变量的代数方程组；步骤6，进行算法效率及网格无关性验证。该针对强非均质地层的油藏数值模拟方法有助于深入了解强非均质油藏地层中的驱油机理，提高油藏数值模拟的效率和精度，指导油藏开发方式的调整与优化。

Reservoir numerical simulation method for strongly heterogeneous formation

【技术实现步骤摘要】
针对强非均质地层的油藏数值模拟方法
本专利技术涉及油藏数值模拟理论研究、技术应用及油气田开发提高油藏采收率领域，特别是涉及到一种针对强非均质地层的油藏数值模拟方法。
技术介绍
油气藏所在地层，在形成过程中由于受到成岩作用、沉积环境、水动力条件等因素的影响，其地质构造以及内部物性都存在着不均匀分布和显著差异；同时，地质运动所产生的岩层褶皱和断裂也会加剧油藏地层的非均质性。除自然因素外，随着油气田开发过程的深入，长期的水驱开采使得国内许多油田已经进入高含水时期，储层的非均质性十分严重。在这些强非均质储层中渗透率分布极度不均匀，其大小可从几个毫达西变化到几千个毫达西，渗透率级差极大。从室内试验观测发现，强非均质条件下油、气、水的运动将出现明显的高渗区之间的斜方向窜流；同时，压力梯度场在渗透率非连续点附近也呈现趋于无穷大的趋势。如何准确地描述渗流场在强非均质条件下的这些新特性，给油藏数值模拟提出了极大的挑战。然而，传统算法基于变量连续假设，在模拟强非均质条件下的渗透问题时，其网格间等效渗透率计算方法失效，致使计算效率和计算精度极低。迄今为止，针对强非均本文档来自技高网...

【技术保护点】
1.针对强非均质地层的油藏数值模拟方法，其特征在于，该针对强非均质地层的油藏数值模拟方法包括:/n步骤1，确定强间断存在区域并划分网络；/n步骤2，计算梯度强度因子并得到局部解析解；/n步骤3，应用局部解析解计算网格间等效渗透率；/n步骤4，构建跨界面相流量计算公式；/n步骤5，求解待求变量的代数方程组；/n步骤6，进行算法效率及网格无关性验证。/n

【技术特征摘要】
1.针对强非均质地层的油藏数值模拟方法，其特征在于，该针对强非均质地层的油藏数值模拟方法包括:
步骤1，确定强间断存在区域并划分网络；
步骤2，计算梯度强度因子并得到局部解析解；
步骤3，应用局部解析解计算网格间等效渗透率；
步骤4，构建跨界面相流量计算公式；
步骤5，求解待求变量的代数方程组；
步骤6，进行算法效率及网格无关性验证。


2.根据权利要求1所述的针对强非均质地层的油藏数值模拟方法，其特征在于，在步骤1，对强非均质油气藏储层区域划分计算网格，判断渗透率发生强间断的非连续点，间断点邻域被分为若干子区域，各子区域内渗透率不变，但子区域间渗透率差异很大；在该强间断点附近，利用边界元方法进行精细数值模拟，发现压力场在该间断点附近的幂律分布特性，且其梯度场在该间断点处发散，即压力梯度在该点趋于无穷大。


3.根据权利要求1所述的针对强非均质地层的油藏数值模拟方法，其特征在于，在步骤2，压力梯度场在渗透率非连续点附近趋于无穷大的特性通过相应的本征指数表征，称之为梯度强度因子；从渗流方程出发通过理论分析，推导得到梯度强度因子的计算方法，进而得到渗透率非连续点附近渗流方程解析形式的通解。


4.根据权利要求3所述的针对强非均质地层的油藏数值模拟方法，其特征在于，在步骤2，二维渗流场在非连续点附近的典型幂律解析解表示为：



下标n表示对应的各子区域，α是梯度强度因子，Mn是各子区域所对应的复系数；Pn表示第n个子区...

【专利技术属性】
技术研发人员：杨勇，杜玉山，邴绍献，孙业恒，于金彪，曹伟东，赵莹莹，施安峰，刘志峰，王晓宏，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司胜利油田分公司勘探开发研究院，
类型：发明
国别省市：山东;37