专利名称：用于储层烃类检测的多尺度岩石物理图版方法及装置的制作方法
技术领域：
本发明涉及地球物理勘探中的岩石物理模型及地震资料反演与解释领域，特别涉及一种用于储层烃类检测的多尺度岩石物理图版方法及装置。
背景技术：
非均质油气藏的储层孔渗关系复杂，岩石孔隙中不相混溶的油、气、水一般互呈“斑块状(补丁状)”分布。单个“补丁”或“斑块”的尺寸大于岩石颗粒与单个孔隙，但远低于地震波长，已接近_量级。这种非均匀性一般被定义为中观尺度的非均匀性。流体的非均匀分布特征导致储层的各类流体分异不明显，这为实际的地震勘探工程带来了难度。为发展非均质油气藏的地震检测技术，以期在复杂的储层环境中精确的定位富气“甜点”，首
要需要查明的是流体非均匀分布的储层岩石中的地震波传播及响应规律。Gassmann与Biot最早针对孔隙含流体的多孔介质中的地震波传播机理进行了研究，即Biot-Gassmann理论。然而传统Biot-Gassmann理论主要考虑的是饱和一种流体的均一孔隙结构，不能描述实际储层的复杂情况。White等(1975)考虑了岩石内部同时含气、水的情况，其中气体呈“补丁状”分布于含水岩石的背景下。Dutta等(1979)改进了 White理论的不足，使其在低频极限下预测的纵波速度吻合于由经典Biot理论出发所导出的分析结果。为解决相同的问题，Johnson提出了一种分支函数法，这种方法在低频与高频极限分别对非饱和介质中波的频散与衰减进行分析，而在中间频段则采用复变函数的方法进行近似逼近，实现了对不同频段内非饱和介质波传播规律的模拟。近十年来，Milller、Gurevich与Toms等对一维与三维的随机非饱和介质进行了模拟,并将模拟结果与碳酸盐岩与砂岩的实测数据进行对比，认为这种随机介质数值模拟的方法也能为实际岩石的波响应提供合理的描述及预测。巴晶、曹宏等(2008-2012)对双重孔隙介质中的地震波进行了进一步的模拟与分析，将双孔理论的预测结果与川中须家河组致密砂岩的宽频带实验结果进行了对比，证实了此类储层中弹性波在地震勘探频段内可能出现强烈的频散现象。为进一步推导一种格式尽可能简洁、物理参数尽可能少、各参量均具备物理可实现性的波传播方程，以期满足实际的应用研究与工业生产的需要，巴晶、Carcione等基于Biot理论框架,导出了描述非饱和岩石地震波传播机理的Biot-Rayleigh方程(简称B-R方程),通过与前人理论的对比,成功将这种方程应用于实际气藏勘探的工程问题中。目前，常规工程所使用的岩石物理图版只能适用于单一尺度的波形资料(Avseth等，2010)，并且缺乏足够的信息与办法对初始的岩石物理图版进行校正，不能保证图版的精度，更不做到多尺度资料的联合应用。综上所述，前人的研究中主要存在以下问题(I)在常规的岩石物理图版技术中，一般采用的流体替换方法是Biot-Gassmann方程，这种方程的前提假设是假定岩石内部的所有孔隙是完全均一的，因此这种假设不能满足非均匀复杂地层的工程应用。
(2)常规的岩石物理图版技术往往针对单一尺度下的波形资料，不能适用于多尺度的波形资料，不能联合利用多尺度弹性波信息来服务于工业生产。(3)由于缺乏足够的先验信息与判据，常规岩石物理图版的精度有限，往往只能辅助定性的解释储层岩性及流体，不能做到对储层岩石参数及流体饱和度的定量反演。
发明内容
本发明的目的在于，提供一种用于储层烃类检测的多尺度岩石物理图版方法及装置，以提高岩石物理图版的精度，并实现对储层岩石参数及流体饱和度的定量反演。为达上述目的，一方面，本发明实施例提供了一种用于储层烃类检测的多尺度岩石物理图版方法，所述方法包括以下步骤获得储层岩石基质模型，获得储层岩石骨架模型，根据所述储层岩石基质模型、所 述储层岩石骨架模型和流体参数获得含流体的岩石模型，并根据所述含流体的岩石模型获得多尺度的初始储层岩石物理图版；根据岩石物理实验数据，在多个频段下校正所述初始储层岩石物理图版；根据测井解释结果，在声波测井频段校正所述初始储层岩石物理图版；根据地震解释结果，在地震频段校正所述初始储层岩石物理图版；在所述初始储层岩石物理图版通过校正后，输出最终的储层岩石物理图版，根据所述最终的储层岩石物理图版和常规叠前地震反演得到的数据体，进行储层岩石物理参数的反演。为达上述目的，另一方面，本发明实施例提供了一种用于储层烃类检测的多尺度岩石物理图版装置，所述装置包括初始图版获取单元，用于获得储层岩石基质模型，获得储层岩石骨架模型，根据所述储层岩石基质模型、所述储层岩石骨架模型和流体参数获得含流体的岩石模型，并根据所述含流体的岩石模型获得多尺度的初始储层岩石物理图版；第一图版校正单元，用于根据岩石物理实验数据，在多个频段下校正所述初始储层岩石物理图版；第二图版校正单元，用于根据测井解释结果，在声波测井频段校正所述初始储层岩石物理图版；第三图版校正单元，用于根据地震解释结果，在地震频段校正所述初始储层岩石物理图版；最终图版获取单元，用于在所述初始储层岩石物理图版通过校正后，输出最终的储层岩石物理图版，根据所述最终的储层岩石物理图版和常规叠前地震反演得到的数据体，进行储层岩石物理参数的反演。本发明实施例提供的技术方案的有益效果是原有常规的岩石物理图版技术中，一般采用的流体替换方法是Biot-Gassmann方程，这种方程的前提假设不能满足非均匀复杂地层的工程应用，且往往针对单一尺度下的波形资料，不能适用于多尺度的波形资料，不能联合利用多尺度弹性波信息来服务于工业生产，只能辅助定性的解释储层岩性及流体。本发明采用Biot-Rayleigh方程进行流体替换，形成了多尺度岩石物理图版，并提出了联合应用多尺度波形资料在各尺度下分别校正对应尺度的岩石物理图版的方案，所得到的多尺度岩石物理图版可更准确的描述储层岩石及波响应，可用于定量估算储层的岩石参数及流体参数。
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做一简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图I为本发明实施例的储层岩石物理图版的建立及校正方法的流程图；图2A-图2B为本发明实施例I :岩石干骨架体积模量与剪切模量理论预测结果；图3A-图3B为本发明实施例2 :含气、水岩石纵波速度与衰减的理论计算结果； 图4A-图4B为本发明实施例3 :岩石物理图版与实验数据的对比；图5为本发明实施例4 :岩石物理图版与测井精细解释结果的对比；图6为本发明实施例5 :岩石物理图版与地震数据的对比；图7A-图7C为本发明实施例6 :基于岩石物理图版对地表地震数据反演得到的孔隙度与饱和度剖面；图8为本发明实施例的一种用于储层烃类检测的多尺度岩石物理图版装置的功能框图；图9为本发明实施例的初始图版获取单元210的具体功能框图。
具体实施例方式为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。对于非均匀、非饱和的实际储层，迄今文献中还未出现一个能普遍的被各研究机构、产业部门所接受的有效的地震波模型。如何将复杂的理论与繁琐的公式应用于实际的工程问题，这个过程本身就面临困难。本发明将基于前期提出的Biot-Rayleigh理论，与前人理论进行对比与分析，将结合岩石基质、骨架建模及流体替换三个步骤，提供针对实际储层的多尺度岩石物理图版，结合叠前地震反演结果进行储层孔隙度与饱和度的反演。图I为本发明实施例的一种用于储层烃类检测的多尺度岩石物理图版方法的流程图。如图I所示，该方法包括如下步骤110、获得储层岩石基质模型，获得储层岩石骨架模型，根据所述储层岩石基质模型、所述储层岩石骨架模型和流体参数获得含流体的岩石模型，并根据所述含流体的岩石模型获得多尺度的初始储层岩石物理图版。具体地，在步骤110中，根据研究区地质情况，应用等效介质理论与经验公式，估算储层岩石干骨架模量，应用Biot-Rayleigh方程进行流体替换，得到含流体岩石的地震响应特征，结合工程需求输出多尺度的岩石物理图版。如图I所不，输入地质报告、岩石实验观测数据和经验资料，以便在后续步骤中利用这些输入信息计算获得岩石基质模型、岩石骨架模型以及得到流体参数。120、根据岩石物理实验数据，在多个频段下校正所述初始储层岩石物理图版。其中，上述多个频段可以包括超声波频段，以及声波频段以下的频段。130、根据测井解释结果，在声波测井频段校正所述初始储层岩石物理图版。具体地，在步骤130中，根据测井解释结果，在声波测井频段(IO3 4Hz)校正岩石物理图版，并验证测井解释结果，若前期测井解释结果存在问题，则重新进行测井资料的精细解释。140、根据地震解释结果，在地震频段校正所述初始储层岩石物理图版。 具体地，在步骤140中，由地震数据中提取井旁数据，结合初步地震解释结果，对储层岩石物理图版进行校正，同时可监测地震数据反演的质量。150、在所述初始储层岩石物理图版通过校正后，输出最终的储层岩石物理图版，根据所述最终的储层岩石物理图版和常规叠前地震反演得到的数据体，进行储层岩石物理参数的反演。需要指出的是，步骤120、130、140可以同时并列执行，也可相应地调整/调换执行顺序。继续参阅图1，以下对步骤110的处理过程进行更为详细的描述111、获取包括地质报告、岩石实验观测数据和经验资料在内的输入信息。根据地质报告可了解岩石矿物的主要组分及配比，可了解储层岩石孔隙度的分布范围，这些信息可限定生成岩石物理图版的输入参数(矿物配比、矿物弹性常数、孔隙率、渗透率)的取值范围，使输入参数有合理性，并且，可辅助设定初始输入参数。岩石实验观测数据可直接测量岩石矿物的弹性参数、密度参数，可直接测量流体在储层温度、压力环境下的弹性参数、密度参数，因此实验数据可用作岩石物理图版的输入参数。经验资料包括来自相似地区相似储层的工程勘测报告、岩芯数据及学术论文，基于这些相似地区的资料，可辅助确定目的储层的岩石参数与孔隙度分布范围。112、获得储层岩石基质模型。在本实施例中，已知岩石中各种矿物成分与其弹性模量，而矿物间的组合方式未知的情况下，采用Voigt-Reuss-Hill平均方法计算岩石的基质弹性参数，其表达式如下
M0 =-(V fM. + 人)g1 fjj_}
TtMiMg——待计算的岩石基质的弹性模量；fi7 Mi——表示第i种矿物的体积率与弹性模量；N——岩石中矿物的种类数。113、获得储层岩石骨架模型。在本实施例中，采用Pride公式计算干燥岩石骨架的体积模量和剪切模量
「 I ^ (I -Kb =—~
I+ C^Mb = \ ；A
Kb, U b-岩石干骨架的体积模量和剪切模量；Ks, U s——岩石基质的体积模量和剪切模量；小——岩石的孔隙度；c’ c,——为经验性参数，与岩石的固结程度有关。对于砂岩，c的取值范围在2到20之间，相同孔隙度的碳酸盐岩c取值较低，一般在0. 5到10之间，岩石骨架结构越紧密，c取值越低。c'的取值一般为c的I. 5到2倍之间。c与c'的确定可根据实测数据调整。114、获得流体参数。在本实施例中，可采用实验方案直接测量流体参数，具体方法是在实验室控制流体的温度与压力与储层环境相同，测量其声速与密度参数，并根据声速与密度计算弹性参数。也可采用流体力学的Van derWaals方程近似计算流体参数。115、根据上述储层岩石基质模型、储层岩石骨架模型和流体参数获得含流体的岩石模型。在本实施例中，基于步骤110获取的参数(例如矿物配比、矿物弹性参数、矿物密度参数、孔隙度、渗透率、流体组分配比、流体弹性参数、流体密度)，采用Voigt-Reuss-Hill方法计算得到基质模型，采用Pride公式计算得到干骨架模型，并将所得的模型参数与步骤114获取的流体参数代入Biot-Rayleigh方程，得到含流体的岩石模型。也即采用Biot-Rayleigh方程对干燥岩石骨架进行流体替换以得到含流体的岩石模型。在对干燥岩石骨架进行流体替换时采用的Biot-Rayleigh方程组具体如下
权利要求
1.一种用于储层烃类检测的多尺度岩石物理图版方法，其特征在于，所述方法包括以下步骤 获得储层岩石基质模型，获得储层岩石骨架模型，根据所述储层岩石基质模型、所述储层岩石骨架模型和流体參数获得含流体的岩石模型，井根据所述含流体的岩石模型获得多尺度的初始储层岩石物理图版； 根据岩石物理实验数据，在多个频段下校正所述初始储层岩石物理图版； 根据测井解释结果，在声波测井频段校正所述初始储层岩石物理图版； 根据地震解释结果，在地震频段校正所述初始储层岩石物理图版； 在所述初始储层岩石物理图版通过校正后，输出最終的储层岩石物理图版，根据所述最終的储层岩石物理图版和常规叠前地震反演得到的数据体，进行储层岩石物理參数的反演。·
2.根据权利要求I所述的方法，其特征在于 所述获得储层岩石基质模型包括采用Voigt-Reuss-Hill平均公式计算骨架基质的等效弹性模量； 所述获得储层岩石骨架模型包括采用Pride公式计算干燥岩石的等效体积模量和剪切模量； 所述根据所述储层岩石基质模型、所述储层岩石骨架模型和流体參数获得含流体的岩石模型包括采用Biot-Rayleigh方程对干燥岩石骨架进行流体替换以得到含流体的岩石模型。
3.根据权利要求I所述的方法，其特征在于，所述根据岩石物理实验数据，在多个频段下校正所述初始储层岩石物理图版包括 在不同频段下，对比所述初始储层岩石物理图版与岩石实验中观测到的波响应数据，调节计算所述初始储层岩石物理图版的各输入參数，使所述初始储层岩石物理图版与实验数据的分布规律一致。
4.根据权利要求3所述的方法，其特征在于，所述实验数据包括 通过共振棒法或者应カ-应变法得到的中、低频段的实验数据；和/或 通过透射波法測量得到的超声波频段的实验数据。
5.根据权利要求I所述方法，其特征在于，所述根据测井解释结果，在声波测井频段校正所述初始储层岩石物理图版包括 采用声波测井数据及相应的声波测井精细解释結果，调节计算所述初始储层岩石物理图版的各输入參数，使测井尺度的初始储层岩石物理图版与测井孔隙度解释结果的分布规律一致。
6.根据权利要求I所述的方法，其特征在于，所述根据地震解释结果，在地震频段校正所述初始储层岩石物理图版包括 根据实际反演得到的地震波响应数据体，在已知井附近提取井旁地震道在目的层段的地震数据，调节计算所述初始储层岩石物理图版的各输入參数，使地震尺度的初始储层岩石物理图版与地震数据的分布规律一致。
7.根据权利要求3所述的方法，其特征在于，所述在不同频段下，对比所述初始储层岩石物理图版与岩石实验中观测到的波响应数据，调节计算所述初始储层岩石物理图版的各输入参数，使所述初始储层岩石物理图版与实验数据的分布规律一致包括 在横轴为纵波阻抗、纵轴为纵横波速度比的坐标系中，将岩石样本的实验数据投影到初始储层岩石物理图版上，同一岩芯在不同含水饱和度下的实验数据表现为一系列散点；通过校正的条件为实验数据散点的连线与初始储层岩石物理图版中对应岩芯的孔隙度的经线接近或者重合，以满足纵波阻抗相对误差不超过5% ；以及，完全含气与完全含水的实验数据点接近饱气纬线与饱水纬线，各含水饱和度下实验数据点与对应含水饱和度下的纬线接近，以满足纵横波速度比相对误差不超过15%。
8.根据权利要求5所述方法，其特征在于，所述采用声波测井数据及相应的声波测井精细解释结果，调节计算所述初始储层岩石物理图版的各输入参数，使测井尺度的初始储层岩石物理图版与测井孔隙度解释结果的分布规律一致包括 在横轴为纵波阻抗、纵轴为纵横波速度比的坐标系中，将目的储层的测井数据的声波时差换算为波阻抗与纵横波速度比，投影到初始岩石物理图版上，测井数据表现为一系列散点； 通过校正的条件为85%以上的测井数据分布于所述初始岩石物理图版描述的范围内，各数据点对应的测井解释孔隙度与相近的经线或其延长线的孔隙度接近，以满足平均纵波阻抗相对误差不超过5%。
9.根据权利要求6所述的方法，其特征在于，所述根据实际反演得到的地震波响应数据体，在已知井附近提取井旁地震道在目的层段的地震数据，调节计算所述初始储层岩石物理图版的各输入参数，使地震尺度的初始储层岩石物理图版与地震数据的分布规律一致包括 在横轴为纵波阻抗、纵轴为纵横波速度比的坐标系中，将井旁道提取的目的储层的地震数据投影到初始储层岩石物理图版上，根据实际的试气结果，将地震数据划分为水层、气层与致密层； 通过校正的条件为 80%以上的目的层地震数据分布于初始储层岩石物理图版的描述范围内，若有20%以上的数据点超出初始储层岩石物理图版的边界范围，则视为不通过； 若有15%以上的水层数据点高于水层纬线，则视为不通过，若有30%以上的水层数据点距同波阻抗气层纬线的距离较距水层纬线的距离更近，则视为不通过； 若有15%以上的气层数据点低于气层纬线，则视为不通过，若50%以上的气层数据点距同波阻抗水层纬线的距离较距气层纬线的距离更近，则视为不通过； 95%以上的致密层数据分布于孔隙度下限所对应的经线右端，其中位于该经线左端的数据点不超过总致密层数据点的5%，否则视为不通过。
10.一种用于储层烃类检测的多尺度岩石物理图版装置，其特征在于，所述装置包括 初始图版获取单元，用于获得储层岩石基质模型，获得储层岩石骨架模型，根据所述储层岩石基质模型、所述储层岩石骨架模型和流体参数获得含流体的岩石模型，并根据所述含流体的岩石模型获得多尺度的初始储层岩石物理图版； 第一图版校正单元，用于根据岩石物理实验数据，在多个频段下校正所述初始储层岩石物理图版；第二图版校正单元，用于根据测井解释结果，在声波测井频段校正所述初始储层岩石物理图版； 第三图版校正単元，用于根据地震解释结果，在地震频段校正所述初始储层岩石物理图版； 最终图版获取单元，用于在所述初始储层岩石物理图版通过校正后，输出最終的储层岩石物理图版，根据所述最終的储层岩石物理图版和常规叠前地震反演得到的数据体，进行储层岩石物理參数的反演。
11.根据权利要求10所述的装置，其特征在于，所述初始图版获取単元包括 储层岩石基质模型计算模块，用于采用Voigt-Reuss-Hill平均模型计算骨架基质的等效弹性模量； 储层岩石骨架模型计算模块，用于采用Pride公式计算干燥岩石的等效体积模量和剪切模量； 储层岩石流体替换计算模块，用于采用Biot-Rayleigh方程对干燥岩石骨架进行流体替换以得到含流体的岩石模型。
12.根据权利要求10所述的装置，其特征在于，所述第一图版校正単元，具体用于在不同频段下，对比所述初始储层岩石物理图版与岩石实验中观测到的波响应数据，调节计算所述初始储层岩石物理图版的各输入參数，使所述初始储层岩石物理图版与实验数据的分布规律一致。
13.根据权利要求10所述的装置，其特征在于，所述第二图版校正単元，具体用于采用声波测井数据及相应的声波测井精细解释結果，调节计算所述初始储层岩石物理图版的各输入參数，使测井尺度的初始储层岩石物理图版与测井孔隙度解释结果的分布规律一致。
14.根据权利要求10所述的装置，其特征在于，所述第三图版校正単元，具体用于根据实际反演得到的地震波响应数据体，在已知井附近提取井旁地震道在目的层段的地震数据，调节计算所述初始储层岩石物理图版的各输入參数，使地震尺度的初始储层岩石物理图版与地震数据的分布规律一致。
15.根据权利要求12所述的装置，其特征在于，所述第一图版校正単元，具体用于在横轴为纵波阻抗、纵轴为纵横波速度比的坐标系中，将岩石样本的实验数据投影到初始储层岩石物理图版上，同一岩芯在不同含水饱和度下的实验数据表现为一系列散点；通过校正的条件为实验数据散点的连线与初始储层岩石物理图版中对应岩芯的孔隙度的经线接近或者重合，以满足纵波阻抗相对误差不超过5% ；以及，完全含气与完全含水的实验数据点接近饱气纬线与饱水纬线，各含水饱和度下实验数据点与对应含水饱和度下的纬线接近，以满足纵横波速度比相对误差不超过15%。
16.根据权利要求13所述的装置，其特征在于，所述第二图版校正単元，具体用于在横轴为纵波阻抗、纵轴为纵横波速度比的坐标系中，将目的储层的测井数据的声波时差換算为波阻抗与纵横波速度比，投影到初始岩石物理图版上，测井数据表现为一系列散点；通过校正的条件为85%以上的测井数据分布于所述初始岩石物理图版描述的范围内，各数据点对应的测井解释孔隙度与相近的经线或其延长线的孔隙度接近，以满足平均纵波阻抗相对误差不超过5%。
17.根据权利要求14所述的装置，其特征在于，所述第三图版校正単元，具体用于在横轴为纵波阻抗、纵轴为纵横波速度比的坐标系中，将井旁道提取的目的储层的地震数据投影到初始储层岩石物理图版上，根据实际的试气结果，将地震数据划分为水层、气层与致密层； 通过校正的条件为 80%以上的目的层地震数据分布于图版描述范围内，若有20%以上的数据点超出图版边界范围，则视为不通过； 若有15%以上的水层数据点高于水层纬线，则视为不通过；若有30%以上的水层数据点距同波阻抗气层纬线的距离较距水层纬线的距离更近，则视为不通过； 若有15%以上的气层数据点低于气层纬线，则视为不通过，若50%以上的气层数据点距同波阻抗水层纬线的距离较距气层纬线的距离更近，则视为不通过； 95%以上的致密层数据分布于孔隙度下限所对应的经线右端，其中位于该经线左端的数据点不超过总致密层数据点的5%，否则视为不通过。
全文摘要
本发明提供一种用于储层烃类检测的多尺度岩石物理图版方法及装置，该方法包括获得储层岩石基质模型，获得储层岩石骨架模型，根据储层岩石基质模型、储层岩石骨架模型和流体参数获得含流体的岩石模型，并根据含流体的岩石模型获得多尺度的初始储层岩石物理图版；根据岩石物理实验数据，在多个频段下校正初始储层岩石物理图版根据测井解释结果，在声波测井频段校正初始储层岩石物理图版；根据地震解释结果，对初始储层岩石物理图版进行校正；在通过校正后，输出最终的储层岩石物理图版，以进行储层岩石物理参数的反演。该方法可提高岩石物理图版的精度并实现对岩石参数及流体饱和度的定量反演。
文档编号G01V1/28GK102854531SQ201210335739
公开日2013年1月2日 申请日期2012年9月11日 优先权日2012年9月11日
发明者巴晶, 晏信飞, 曹宏, 徐光成, 卢明辉, 孙卫涛 申请人:中国石油天然气股份有限公司