碳酸盐岩的参数预测方法以及装置制造方法及图纸

本申请公开了一种碳酸盐岩的参数预测方法，涉及地质勘测领域。该方法包括：通过测井数据中每个深度点的至少一种矿物基质的元素含量以及随机生成的每个深度点的溶蚀结构的预测含量，计算每个深度点的岩石中的第二纵波速度；判断测井数据中第一纵波速度和计算得到的第二纵波速度的误差是否大于阈值；当上述误差大于阈值时通过模拟退火算法生成新的溶蚀结构的预测含量，并再次计算上述每个深度点的岩石中的第二纵波速度；否则，确定每个深度点的溶蚀结构的预测含量为每个深度点的溶蚀结构的含量。该技术方案解决了碳酸盐岩的溶蚀孔、缝、洞参数含量的预测问题。

【技术实现步骤摘要】
碳酸盐岩的参数预测方法以及装置
本专利技术涉及地质勘测领域，特别涉及一种碳酸盐岩的参数预测方法以及装置。
技术介绍
碳酸盐岩中溶蚀孔、缝、洞参数是碳酸盐岩储层中一项很重要的参数，比如，该参数可以用于石油的勘探。在各向异性岩石物理模型的理论建立阶段，本领域的技术人员相继提出了椭球体的有效弹性模量的计算方法，裂缝介质的杨氏模量、剪切模量和泊松比的计算方法，以及针对单套裂缝的理论。而到了各向异性岩石物理模型的理论发展阶段，技术研究的方向从各项同性背景中的单套裂缝转入了多套裂缝的研究，主要研究了基于不同介质的裂缝模型，比如基于各项同性(TransverseIsotropy，TI)介质、垂直的各项同性(VerticalTransverseIsotropy，VTI)介质、以及水平的各项同性(HorizontalTransverseIsotropy，HTI)介质的裂缝模型的研究。针对碳酸盐岩储层，提出了碳酸盐岩中印模孔的预测方法；反演得出基质孔隙度、次生孔隙度和孔隙类型等物性参数；还提出了裂缝方位各项异性的预测方法。但是上述研究的结果都是关于碳酸盐岩的孔隙结构的，并没有一种方法能够得出碳酸盐岩的溶蚀孔、缝、洞在碳酸盐岩中的含量。
技术实现思路
本申请实施例提供了一种碳酸盐岩的参数预测方法以及装置，可以解决碳酸盐岩的溶蚀孔、缝、洞参数含量的预测问题。所述技术方案如下：根据本申请的一个方面，提供了一种碳酸盐岩的参数预测方法，应用于服务器中，该方法包括：获取测井数据中每个深度点的至少一种矿物基质的元素含量和岩石中第一纵波速度；根据上述至少一种矿物基质的元素含量，计算每个深度点的岩石中矿物基质的弹性模量；随机生成每个深度点的溶蚀结构的预测含量；根据上述岩石中矿物基质的弹性模量和上述溶蚀结构的预测含量，计算每个深度点的岩石中的第二纵波速度；当上述第一纵波速度与上述第二纵波速度之间的误差大于阈值时，通过模拟退火算法生成上述溶蚀结构的预测含量，并再次根据上述岩石中矿物基质的弹性模量和上述溶蚀结构的预测含量计算每个深度点的岩石中的第二纵波速度；当误差小于或等于阈值时，确定上述溶蚀结构的预测含量为每个深度点的岩石中溶蚀结构的含量。在一些实施例中，上述溶蚀结构的预测含量包括溶蚀孔的预测含量、溶蚀洞的预测含量和溶蚀缝的预测含量；根据上述岩石中矿物基质的弹性模量和上述溶蚀结构的预测含量，计算每个深度点的岩石中的第二纵波速度，包括：根据上述岩石中矿物基质的弹性模量、上述溶蚀孔的预测含量和上述溶蚀洞的预测含量计算每个深度点的干燥岩石的弹性模量；根据上述干燥岩石的弹性模量和上述溶蚀缝的预测含量计算每个深度点的干燥岩石的第一弹性张量；根据上述干燥岩石的第一弹性张量计算每个深度点的饱和岩石的第二弹性张量；根据上述饱和岩石的第二弹性张量计算每个深度点的岩石中的第二纵波速度。在一些实施例中，上述岩石中矿物基质的弹性模量包括上述岩石中矿物基质的体积模量和剪切模量；根据上述岩石中矿物基质的弹性模量、上述溶蚀孔的预测含量和上述溶蚀洞的预测含量计算每个深度点的干燥岩石的弹性模量，包括：根据上述岩石中矿物基质的体积模量和剪切模量、上述溶蚀孔的预测含量和上述溶蚀洞的预测含量，通过双相介质包含物(Kuster-Toksoz，KT)模型计算每个深度点的干燥岩石的弹性模量。在一些实施例中，根据上述岩石中矿物基质的体积模量和剪切模量、上述溶蚀孔的预测含量和上述溶蚀洞的预测含量，通过双相介质包含物KT模型计算每个深度点的干燥岩石的弹性模量，包括：将上述岩石中矿物基质的体积模量和剪切模量代入包含体Berryman模型的三维孔洞形状公式中，得到KT模型的模型参数；将上述岩石中矿物基质的体积模量和剪切模量、上述溶蚀孔的预测含量、上述溶蚀洞的预测含量以及上述模型参数代入KT模型中，计算每个深度点的干燥岩石的弹性模量。在一些实施例中，上述干燥岩石的弹性模量包括上述干燥岩石的体积模量和剪切模量；根据上述干燥岩石的弹性模量、上述溶蚀缝的预测含量计算每个深度点的干燥岩石的第一弹性张量，包括：将上述干燥岩石的体积模量和剪切模量、以及上述溶蚀缝的预测含量代入裂缝介质Hudson模型中，计算每个深度点的干燥岩石的第一弹性张量。在一些实施例中，根据上述至少一种矿物基质的元素含量，计算每个深度点的岩石中矿物基质的弹性模量，包括：将上述至少一种矿物基质的元素含量代入混合矿物的弹性模量VRH模型中，计算每个深度点的岩石中矿物基质的弹性模量。在一些实施例中，该方法还包括：采用相对误差计算公式计算上述第一纵波速度与上述第二纵波速度之间的误差。根据本申请的另一个方面，提供了一种碳酸盐岩的参数预测装置，该装置包括：获取模块，用于获取测井数据中每个深度点的至少一种矿物基质的元素含量和岩石中第一纵波速度；计算模块，用于根据上述至少一种矿物基质的元素含量，计算每个深度点的岩石中矿物基质的弹性模量；生成模块，用于随机生成每个深度点的溶蚀结构的预测含量；计算模块，用于根据上述岩石中矿物基质的弹性模量和上述溶蚀结构的预测含量，计算每个深度点的岩石中的第二纵波速度；确定模块，用于当上述第一纵波速度与上述第二纵波速度之间的误差大于阈值时，通过模拟退火算法生成上述溶蚀结构的预测含量，并再次根据上述岩石中矿物基质的弹性模量和上述溶蚀结构的预测含量计算每个深度点的岩石中的第二纵波速度；当误差小于或等于阈值时，确定上述溶蚀结构的预测含量为每个深度点的岩石中溶蚀结构的含量。在一些实施例中，上述溶蚀结构的预测含量包括溶蚀孔的预测含量、溶蚀洞的预测含量和溶蚀缝的预测含量；计算模块，用于根据上述岩石中矿物基质的弹性模量、上述溶蚀孔的预测含量和上述溶蚀洞的预测含量计算每个深度点的干燥岩石的弹性模量；根据上述干燥岩石的弹性模量和上述溶蚀缝的预测含量计算每个深度点的干燥岩石的第一弹性张量；根据上述干燥岩石的第一弹性张量计算每个深度点的饱和岩石的第二弹性张量；根据上述饱和岩石的第二弹性张量计算每个深度点的岩石中的第二纵波速度。在一些实施例中，上述岩石中矿物基质的弹性模量包括上述岩石中矿物基质的体积模量和剪切模量；计算模块，用于根据上述岩石中矿物基质的体积模量和剪切模量、上述溶蚀孔的预测含量和上述溶蚀洞的预测含量，通过双相介质包含物KT模型计算每个深度点的干燥岩石的弹性模量。在一些实施例中，计算模块，用于将上述岩石中矿物基质的体积模量和剪切模量代入包含体Berryman模型的三维孔洞形状公式中，得到KT模型的模型参数；将上述岩石中矿物基质的体积模量和剪切模量、上述溶蚀孔的预测含量、上述溶蚀洞的预测含量以及上述模型参数代入KT模型中，计算每个深度点的干燥岩石的弹性模量。在一些实施例中，上述干燥岩石的弹性模量包括上述干燥岩石的体积模量和剪切模量；计算模块，用于将上述干燥岩石本文档来自技高网...

【技术保护点】
1.一种碳酸盐岩的参数预测方法，应用于服务器中，其特征在于，所述方法包括：/n获取测井数据中每个深度点的至少一种矿物基质的元素含量和岩石中第一纵波速度；/n根据所述至少一种矿物基质的元素含量，计算所述每个深度点的岩石中矿物基质的弹性模量；/n随机生成所述每个深度点的溶蚀结构的预测含量；/n根据所述岩石中矿物基质的弹性模量和所述溶蚀结构的预测含量，计算所述每个深度点的岩石中的第二纵波速度；/n当所述第一纵波速度与所述第二纵波速度之间的误差大于阈值时，通过模拟退火算法生成所述溶蚀结构的预测含量，并再次根据所述岩石中矿物基质的弹性模量和所述溶蚀结构的预测含量计算所述每个深度点的岩石中的第二纵波速度；当所述误差小于或等于阈值时，确定所述溶蚀结构的预测含量为所述每个深度点的岩石中溶蚀结构的含量。/n

【技术特征摘要】
1.一种碳酸盐岩的参数预测方法，应用于服务器中，其特征在于，所述方法包括：
获取测井数据中每个深度点的至少一种矿物基质的元素含量和岩石中第一纵波速度；
根据所述至少一种矿物基质的元素含量，计算所述每个深度点的岩石中矿物基质的弹性模量；
随机生成所述每个深度点的溶蚀结构的预测含量；
根据所述岩石中矿物基质的弹性模量和所述溶蚀结构的预测含量，计算所述每个深度点的岩石中的第二纵波速度；
当所述第一纵波速度与所述第二纵波速度之间的误差大于阈值时，通过模拟退火算法生成所述溶蚀结构的预测含量，并再次根据所述岩石中矿物基质的弹性模量和所述溶蚀结构的预测含量计算所述每个深度点的岩石中的第二纵波速度；当所述误差小于或等于阈值时，确定所述溶蚀结构的预测含量为所述每个深度点的岩石中溶蚀结构的含量。


2.根据权利要求1所述的方法，其特征在于，所述溶蚀结构的预测含量包括溶蚀孔的预测含量、溶蚀洞的预测含量和溶蚀缝的预测含量；
所述根据所述岩石中矿物基质的弹性模量和所述溶蚀结构的预测含量，计算所述每个深度点的岩石中的第二纵波速度，包括：
根据所述岩石中矿物基质的弹性模量、所述溶蚀孔的预测含量和所述溶蚀洞的预测含量计算所述每个深度点的干燥岩石的弹性模量；
根据所述干燥岩石的弹性模量和所述溶蚀缝的预测含量计算所述每个深度点的干燥岩石的第一弹性张量；
根据所述干燥岩石的第一弹性张量计算所述每个深度点的饱和岩石的第二弹性张量；
根据所述饱和岩石的第二弹性张量计算所述每个深度点的岩石中的第二纵波速度。


3.根据权利要求2所述的方法，其特征在于，所述岩石中矿物基质的弹性模量包括所述岩石中矿物基质的体积模量和剪切模量；
所述根据所述岩石中矿物基质的弹性模量、所述溶蚀孔的预测含量和所述溶蚀洞的预测含量计算所述每个深度点的干燥岩石的弹性模量，包括：
根据所述岩石中矿物基质的体积模量和剪切模量、所述溶蚀孔的预测含量和所述溶蚀洞的预测含量，通过双相介质包含物KT模型计算所述每个深度点的干燥岩石的弹性模量。


4.根据权利要求3所述的方法，其特征在于，所述根据所述岩石中矿物基质的体积模量和剪切模量、所述溶蚀孔的预测含量和所述溶蚀洞的预测含量，通过双相介质包含物KT模型计算所述每个深度点的干燥岩石的弹性模量，包括：
将所述岩石中矿物基质的体积模量和剪切模量代入包含体Berryman模型的三维孔洞形状公式中，得到所述KT模型的模型参数；
将所述岩石中矿物基质的体积模量和剪切模量、所述溶蚀孔的预测含量、所述溶蚀...

【专利技术属性】
技术研发人员：彭达，肖富森，冉崎，谢冰，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京;11