一种自适应的岩石物理建模方法技术

本发明专利技术涉公开了一种自适应的岩石物理建模方法，对每一个测井深度点计算各种岩石物理模型的优先级参数，最高优先级的岩石物理模型计算纵横波速度参数、密度参数与测井纵横波速度及密度的整体误差目标函数，将误差保持在阈值内，进行岩石物理模型的筛选，特别适用于工区存在纵向变化快的特征，对不同深度点采用最合适的岩石物理模型，提高了岩石物理模型对岩石的描述和预测的准确性。石的描述和预测的准确性。石的描述和预测的准确性。

【技术实现步骤摘要】
一种自适应的岩石物理建模方法


[0001]本专利技术涉及油气地球物理领域，特别是一种自适应的岩石物理建模方法。

技术介绍

[0002]过去的几十年间，在研究与地震勘探学和天然地震学有关的岩石物理特性方面取得了巨大的进步，在此期间，发展了很多的理论，同时进行了很多实验，许多岩石物理理论与实验结果已经在推进地球科学和勘探技术方面起到了重要准则的作用。在勘探地震学中，地震波以旅行时间、反射波振幅及相位变化的形式带来了地下岩石和流体的信息。在早期的勘探地震学中，地震数据主要用作构造解释，这些构造可能含有油气。随着计算能力的提高和地震处理、解释技术的进步，现在对地震数据的分析一般是为了预测岩性、孔隙度、孔隙流体以及饱和度。因为岩石物理学为地震数据与油藏特性和参数之间架起了桥梁，近年来它已在有关新技术的开发中发挥作用，诸如4D地震油藏监测，地震岩性识别，以及“亮点”和反射系数随人射角变化的分析等油气直接检测技术。
[0003]岩石物理是把油藏特性和参数与地震数据相连的桥梁，在许多试图从地震数据中导出储层岩石和流体特性及油藏参数的技术方法中，它起到一种基本准则的作用。随着越来越多的油气田进入成熟期以及许多方法技术的进一步改进和发展，包括孔隙度和流体流的成图，储层压力、温度、饱和度变化的测定，甚至推断含油气的类型。在这样一种发展过程中，岩石物理将是定量解释地震数据的一项十分重要的技术手段。
[0004]如何准确地构建实用和准确的岩石物理模型是储层物性估计的岩石物理新模型构建部分。以横波估计为例，目前，在没有横波信息的情况下，估计横波速度的方法主要有两种，一是统计拟合法；二是理论公式法。但是这些方法往往具有一定的局限性。针对存在的局限性，很多学者基于岩石物理模型，进行了横波速度预测方法的综合研究，Greenberg等假设纵横波速度间有稳健的关系，基于Biot-Gassmann理论预测横波速度，Xu等使用Kutster-Toksoz理论和微分等效介质理论相结合预测横波速度，并运用孔隙纵横比的概念来表征干岩石颗粒的接触关系，等运用有效介质理论预测横波速度，并且认为基于岩石物理的方法的横波预测的准确度要高于统计拟合法得到的横波速度。国内方面，也有很多学者开展了基于岩石物理模型的横波预测研究，并取得了较好的研究成果。
[0005]目前主流的岩石物理模型往往都存在着各自较为适用的范围，如在常规砂岩储层中常使用Xu-White模型，在碳酸盐岩中通常使用Xu-Payne模型，在页岩中常使用各向异性SCA模型等等，同时各模型还存在着更多针对不同特殊情况的变形或者组合形式，针对某一组分、结构岩石来说，总能找到最适用的岩石物理模型。
[0006]然而，目前岩石物理建模过程中，通常针对某一工区只建立一种岩石物理模型，但部分工区存在着纵向变化快的特征，只使用一种岩石物理模型去描述整个工区的岩石变化特征则不尽适用，使得岩石物理模型对岩石的描述和预测的准确性下降。

技术实现思路

[0007]本专利技术的专利技术目的在于：针对现有技术中岩石物理建模存在的对岩石描述和预测的准确性低的问题，提供一种自适应的岩石物理建模方法。
[0008]为了实现上述目的，本专利技术采用的技术方案为：
[0009]一种自适应的岩石物理建模方法，包括：
[0010]S100对测井数据进行解释，获得测井参数；
[0011]S200对测井参数中的测井矿物数据进行规则化处理，处理为囊括该工区内所有矿物类型的标准测井矿物数据；
[0012]S300基于规则化后的测井矿物数据，针对每一个测井深度点计算各种岩石物理模型的优先级参数，并根据优先级参数进行岩石物理模型的优先级排序；
[0013]S400基于最高优先级的岩石物理模型默认矿物模量，计算出岩石基质模量，再根据默认孔隙类型、孔隙纵横比，并通过Gassmann公式获得测井流体状态下的纵横波速度参数和密度参数；
[0014]S500计算步骤S400得到的纵横波速度参数、密度参数与测井纵横波速度及密度的整体误差目标函数，误差高于阈值则对最高优先级的岩石物理模型默认参数进行优化，重复步骤S400重新计算获得新的纵横波速度参数和密度参数；
[0015]S600若最高优先级的岩石物理模型的默认参数优化到默认极限后，误差仍然高于阈值，则选择更低优先级的岩石物理模型重复步骤S400、S500，直到模型匹配性低于阈值时，获得与该深度点测井数据相匹配的岩石物理模型。
[0016]一种自适应的岩石物理建模方法，对每一个测井深度点计算各种岩石物理模型的优先级参数，最高优先级的岩石物理模型计算纵横波速度参数、密度参数与测井纵横波速度及密度的整体误差目标函数，将误差保持在阈值内，进行岩石物理模型的筛选，特别适用于工区存在纵向变化快的特征，对不同深度点采用最合适的岩石物理模型，提高了岩石物理模型对岩石的描述和预测的准确性。
[0017]优选的，还包括步骤S700：若所有备选的岩石物理模型都不匹配，则回到步骤S100进行测井参数的优化再解释。
[0018]优选的，其特征在于，所述测井参数包括密度、孔隙度、含水饱和度、矿物含量。
[0019]优选的，所述测井矿物数据包括石英含量、粘土含量、方解石含量、白云石含量、膏岩含量。
[0020]优选的，步骤S200的规则化处理为：将每个深度点的矿物数据都扩充为该工区所有矿物，其中不存在于对应深度点的矿物的值设置为0。
[0021]优选的，步骤S300中任一深度点的岩石物理模型优先级参数的计算方法：
[0022]S310在岩石物理模型集合中，若测井数据未提供要求提供的数据，则将对应的岩石物理模型从集合中剔除，剩下的即为备选的岩石物理模型；
[0023]S320在备选的岩石物理模型中，将深度点的矿物进行主要成分分析，得到每个岩石物理模型的优先级参数。
[0024]优选的，所述岩石物理模型集合包括DEM微分等效模型，SCA自洽模型，KT包裹体模型，Xu-White模型，Xu-Payne模型，Hudson裂隙介质模型，各向异性SCA-DEM模型，胶结砂岩模型，未胶结砂岩模型。
[0025]优选的，步骤S400中采用Voigt-Reuss-Hill平均计算岩石基质模量：
[0026][0027][0028][0029]其中，M
m
为碳酸盐岩饱和石英和泥质骨架的岩石等效弹性模量，M
i
为碳酸盐岩组成岩石矿物成分第i个组成成分的模量，f
i
为碳酸盐岩组成岩石矿物成分第i个组成成分的体积分量，M
V
为采用Voigt上限方法计算获得的岩石模量，M
R
为采用Reuss下限方法计算获得的岩石模量。
[0030]优选的，步骤S500中的目标函数：
[0031]有横波数据的情况：L
s
＝(V
pi-V
′
pi
)2(V
si-V
′
本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种自适应的岩石物理建模方法，其特征在于，包括：S100对测井数据进行解释，获得测井参数；S200对测井参数中的测井矿物数据进行规则化处理，处理为囊括该工区内所有矿物类型的标准测井矿物数据；S300基于规则化后的测井矿物数据，针对每一个测井深度点计算各种岩石物理模型的优先级参数，并根据优先级参数进行岩石物理模型的优先级排序；S400基于最高优先级的岩石物理模型默认矿物模量，计算出岩石基质模量，再根据默认孔隙类型、孔隙纵横比，并通过Gassmann公式获得测井流体状态下的纵横波速度参数和密度参数；S500计算步骤S400得到的纵横波速度参数、密度参数与测井纵横波速度及密度的整体误差目标函数，误差高于阈值则对最高优先级的岩石物理模型默认参数进行优化，重复步骤S400重新计算获得新的纵横波速度参数和密度参数；S600若最高优先级的岩石物理模型的默认参数优化到默认极限后，误差仍然高于阈值，则选择更低优先级的岩石物理模型重复步骤S400、S500，直到模型匹配性低于阈值时，获得与该深度点测井数据相匹配的岩石物理模型。2.根据权利要求1所述的方法，其特征在于，还包括步骤S700：若所有备选的岩石物理模型都不匹配，则回到步骤S100进行测井参数的优化再解释。3.根据权利要求1-2任一项所述的方法，其特征在于，所述测井参数包括密度、孔隙度、含水饱和度、矿物含量。4.根据权利要求3所述的方法，其特征在于，所述测井矿物数据包括石英含量、粘土含量、方解石含量、白云石含量、膏岩含量。5.根据权利要求4所述的方法，其特征在于，步骤S200的规则化处理为：将每个深度点的矿物数据都扩充为该工区所有矿物，其中不存在于对应深度点的矿物的值设置为0。6.根据权利要求5所述的方法，其特征在于，步骤S300中任一深度点的岩石物理模型优先级参数的计算方法：S310在岩石物理模型集合中，若测井数据未提供要求提供的数据，则将对应的岩石物理模型从集合中剔除，剩下的即为备选的岩石物理模型；S320在备选的岩石物理模型中，将深度点的矿物进行主要成分分析，得到每个岩石物理模型的优先级参数。7.根据权利要求6所述...

【专利技术属性】
技术研发人员：赵爽，丁蔚楠，叶泰然，刘成川，詹国卫，毕有益，董霞，王荐，张岩，喻勤，
申请(专利权)人：中国石油化工股份有限公司西南油气分公司，
类型：发明
国别省市：