申请/专利权人：北京阳光杰科科技股份有限公司

申请日：2018-09-29

公开（公告）日：2024-03-22

公开（公告）号：CN110968928B

主分类号：G06F30/20

分类号：G06F30/20;G06T17/05;G01V1/30

优先权：

专利状态码：有效-授权

法律状态：2024.03.22#授权;2020.05.01#实质审查的生效;2020.04.07#公开

摘要：本发明提供一种低渗致密砂岩等效三孔隙岩石物理建模方法及其应用，所述方法通过三维数字岩芯扫描将岩芯孔隙分离并确定不同孔隙类型初始纵横比，再通过等效三孔隙类型纵横比和自洽模型进行岩石物理图版建立，最后通过图版进行随机模拟和计算，输出各孔隙类型的孔隙度和最优化横波速度，通过所述方法，建立地震岩石物理弹性参数与表征油藏属性参数间的关系，从而应用地震资料对油藏特性进行研究，以三维数字岩芯扫描和岩芯岩石物理测试数据为基础，针对低渗致密砂岩岩石建立等效三孔隙岩石物理模型，具有模型精度高、实用性和创新性的特点。

主权项：1.一种低渗致密砂岩等效三孔隙岩石物理建模方法，其特征在于，所述方法通过三维数字岩芯扫描将岩芯孔隙分离并确定不同孔隙类型初始纵横比，再通过等效三孔隙类型纵横比和自洽模型进行岩石物理图版建立，最后通过图版进行随机模拟和计算，输出各孔隙类型的孔隙度和最优化横波速度；所述方法包括以下步骤：1：采用三维数字岩芯扫描，将岩芯孔隙分离，在此基础上，对每个孔隙进行定量分析，根据等效原理统计孔隙的等效孔径，以此划分孔隙类型，根据等效孔径确定不同孔隙类型初始纵横比；2:根据岩芯岩石物理测试数据，采用Reuss-Voigt-Hill平均建立骨架模型，采用Brie模型建立流体模型，根据第一步确定的等效三孔隙类型纵横比和自洽模型，建立等效三孔隙类型岩石物理图版；3:在第二步岩石物理图版基础上，投影区域测井数据，进一步确定硬孔纵横比αstiff和裂隙纵横比αcrack，同时对粒间孔隙纵横比αinter设定一定的值域范围，采用蒙特卡罗随机模拟，计算最大后验概率，输出各孔隙类型的孔隙度和最优化横波速度；所述步骤1划分孔隙类型中等效三孔隙包括硬孔、粒间孔和微裂隙；所述步骤2具体为采用岩石物理实验测试，获得岩石的矿物组分及纵波、横波、密度弹性参数和物性参数，根据岩芯分析的矿物组分，采用Reuss-Voigt-Hill平均建立骨架模型，采用Brie模型建立流体模型，Brie方程表达式为Kfl＝Kwater-KgasSwatere+Kgas，其中Kgas、Kwater分别为气水的体积模量，Swater为含水饱和度，e是经验系数；所述步骤3中采用自洽模型建立等效三孔隙类型岩石物理图版时，当岩芯取样深度变化范围大时，在不同深度段建立多个图版；所述步骤3中由于低渗-致密砂岩孔隙之间连通性差，采用Brie模型混合的流体模量直接加入自洽模型的迭代中，将SPM的等效双孔隙类型岩石物理模型推广到三孔隙类型，更符合岩石的实际情况；所述方法还包括步骤4：在步骤1-步骤3的基础上，固定微裂缝纵横比αcrack和硬孔隙纵横比αstiff，给定粒间孔孔隙纵横比αinter数值区间在0.15-0.37之间，用于表征不同沉积环境下的岩石孔隙形态，本方法更符合地质情况，给定等效粒间孔隙纵横比区间，根据贝叶斯算法建立粒间孔隙纵横比αinter与粒间孔隙度Φinter优化的目标函数，目标函数为在总孔隙度Φtotal、含水饱和度Sw和骨架弹性模量Kmatrix，Umatrix已知条件下的最大后验概率分布，表达式为[αinter,Φinter]＝argMax{P[αinter,Φinter]j[Φtotal，Sw，Kmatrix，Umatrix]}，其中P[αinter,Φinter]j[Φtotal，Sw，Kmatrix，Umatrix]为粒间孔隙纵横比及粒间孔隙度的最大后验分布，意义是已知Φtotal，Sw，Kmatrix，Umatrix，计算αinter,Φinter的后验概率分布，当后验概率取最大值时，该值就是最终的结果，优化算法采用蒙特卡罗随机算法求解目标函数。

全文数据：

权利要求：

百度查询： 北京阳光杰科科技股份有限公司 一种低渗致密砂岩等效三孔隙岩石物理建模方法及其应用

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD