测井速度预测、基于井震融合的孔隙度确定方法及装置制造方法及图纸

本说明书公开了测井速度预测、基于井震融合的孔隙度确定方法及装置，其中测井速度预测方法包括：获取目标储层的声波测井数据；根据目标储层中各类矿物的含量计算得到第一体积模量和第一剪切模量；根据储层中目标深度位置处的孔隙度、第一体积模量和第一剪切模量，生成以孔隙纵横比为变量的第一目标函数；将第一目标函数取值最小时的孔隙纵横比作为目标孔隙纵横比；根据目标孔隙纵横比计算得到第二体积模量和第二剪切模量；根据第二体积模量、第二剪切模量计算得到频率的声波测井信号在目标储层中传播时的横波速度和/或纵波速度。本方案能够对稠油储层中测井信号的传播速度进行预测，为基于传播速度尤其是横波速度的地震岩石属性反演建立了基础。岩石属性反演建立了基础。岩石属性反演建立了基础。

【技术实现步骤摘要】
测井速度预测、基于井震融合的孔隙度确定方法及装置


[0001]本申请涉及石油、天然气地震勘探反演和定量解释
，特别涉及测井速度预测、基于井震融合的孔隙度确定方法及装置。

技术介绍

[0002]地震岩石物理学通过分析测井资料、地震资料和岩石岩性、物性资料，研究岩石样品的弹性参数、物性参数等基本交汇关系，分析岩石的岩石物理规律。当地震波在地下岩石中传播时，岩石中的流体与颗粒骨架之间的相互作用将改变地震波的旅行时间、振幅以及相位等信息。通过地震波的信息来分析地下多孔岩石以及流体的方法也越来越受到重视。地震岩石物理学的发展与进步也使得通过地震数据预测储层物性参数(孔隙度、渗透率)得以实现。
[0003]横波是地震波中前进方向与质点振动方向垂直的波，横波在判断岩性、裂缝以及含油气性方面有其固有的优点。横波测井曲线对于地震资料反演及准确求取地层的物性与弹性参数有重要作用。在石油天然气的地震勘探中，往往需要输入横波测井曲线作为地质模型的参数用于地震波场的数值模拟或地震属性的参数反演。
[0004]通常情况下，开发井中只有极少数井有横波测井资料，而开发时间较长的老油田则横波测井资料更少。在缺少横波测井资料的情况下，从纵波测井资料或其它测井数据中预测横波测井资料就显得较为重要。在缺少横波测井资料的老油田，横波测井预测能够为储层岩性识别、岩石弹性参数计算和地震流体识别等技术的研究和应用提供横波测井曲线，有利于提高地震储层预测及油气检测的精度。此外，在地震岩石属性反演研究中，尤其是地震数据叠前反演中，也需要使用横波速度测井资料。

技术实现思路

[0005]本申请实施方式的目的是提供一种测井速度预测、基于井震融合的孔隙度确定方法及装置，以基于测井数据预测测井信号的传播速度。
[0006]为解决上述技术问题，本说明书第一方面提供一种测井速度预测方法，包括：获取目标储层的声波测井数据，所述声波测井数据包括：声波测井信号的频率、储层中目标深度位置处的孔隙度、各类矿物的含量；根据所述目标储层中各类矿物的含量计算得到第一体积模量和第一剪切模量；根据所述储层中目标深度位置处的孔隙度、所述第一体积模量和所述第一剪切模量，生成以孔隙纵横比为变量的第一目标函数；将所述第一目标函数取值最小时的孔隙纵横比作为目标孔隙纵横比；根据所述目标孔隙纵横比计算得到第二体积模量和第二剪切模量；根据所述第二体积模量、所述第二剪切模量计算得到所述频率的声波测井信号在所述目标储层中传播时的横波速度和/或纵波速度。
[0007]本说明书第二方面提供一种基于井震融合的孔隙度确定方法，包括：获取目标储层的声波测井数据，所述声波测井数据包括：声波测井信号的频率、储层中目标深度位置处的孔隙度、各类矿物的含量；根据所述目标储层中各类矿物的含量计算得到第一体积模量
和第一剪切模量；根据所述储层中目标深度位置处的孔隙度、所述第一体积模量和所述第一剪切模量，生成以孔隙纵横比为变量的第一目标函数；将所述第一目标函数取值最小时的孔隙纵横比作为目标孔隙纵横比；根据所述目标孔隙纵横比计算得到第二体积模量和第二剪切模量；根据所述第二体积模量、所述第二剪切模量计算得到所述频率的声波测井信号在所述目标储层中传播时的纵波速度；根据目标储层中各个深度位置的声波测井数据计算得到各个深度位置对应的孔隙纵横比；根据目标储层中多个深度位置对应的孔隙纵横比确定目标储层的孔隙纵横比的替代值；根据所述目标储层的孔隙纵横比的替代值，生成以孔隙度为变量的第二目标函数；将所述第二目标函数取值最小时的孔隙度作为目标孔隙度。
附图说明
[0008]为了更清楚地说明本申请实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0009]图1示出了对声波测井数据和地震资料进行岩石物理的应用的方法示意图；
[0010]图2示出了不同温度条件下CCM模型计算的稠油复剪切模量随频率的变化关系；
[0011]图3A示出了利用CPA模型模拟不同孔隙纵横比条件下，体积模量随孔隙度的变化关系；
[0012]图3B示出了利用CPA模型模拟不同孔隙纵横比条件下，剪切模量随孔隙度的变化关系；
[0013]图4A示出了岩石的体积模量随频率的变化关系示意图；
[0014]图4B示出了岩石的剪切模量随频率的变化关系示意图；
[0015]图5A示出了改进的CPA模型模拟的不同含油饱和度条件下岩石的体积模量随频率的变化关系；
[0016]图5B示出了改进的CPA模型模拟的不同含油饱和度条件下岩石的剪切模量随频率的变化关系；
[0017]图6示出了根据本说明书实施例的一种测井速度预测方法的流程图；
[0018]图7A示出了X井的纵波速度的预测值和实测值；
[0019]图7B示出了X井的横波速度的预测值；
[0020]图7C示出了X井的孔隙纵横比的预测值；
[0021]图7D示出了Y井的纵波速度的预测值和实测值；
[0022]图7E示出了Y井的横波速度的预测值和实测值；
[0023]图7F示出了Y井的孔隙纵横比的预测值；
[0024]图8示出了根据本说明书实施例的一种基于井震融合的孔隙度确定方法的流程图；
[0025]图9示出了预测得到的孔隙度曲线图；
[0026]图10示出了将声波纵波速度的实测值校正到地震频段的方法流程图；
[0027]图11A示出了根据本说明书实施例的一种测井速度预测装置的原理框图；
[0028]图11B示出了根据本说明书实施例的一种基于井震融合的孔隙度确定装置的原理框图；
[0029]图12示出了根据本说明书实施例的电子设备的原理框图。
具体实施方式
[0030]为了使本
的人员更好地理解本申请中的技术方案，下面将结合本申请实施方式中的附图，对本申请实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式仅仅是本申请一部分实施方式，而不是全部的实施方式。基于本申请中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都应当属于本申请保护的范围。
[0031]本专利技术采用CCM模型表征饱和稠油岩石的粘弹性性质，引入饱和度的概念改进了CPA模型，使CPA模型能够表征稠油、水两种流体分布。改进的CPA模型具有频变属性，利用CPA模型同时对声波测井数据和地震资料进行岩石物理的应用。
[0032]具体地，如图1所示，该应用方法对工区样品进行了高频(例如0
‑
1kHz)的声波测井试验和低频(例如0
‑
30Hz)的地震试验，利用实验数据获得岩石样品的孔隙纵横比的范围。基于改进的CPA模型对稠油储层进行横波速度预测，并将声波测井数据校正到地震频段，利用地震频本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种测井速度预测方法，其特征在于，包括：获取目标储层的声波测井数据，所述声波测井数据包括：声波测井信号的频率、储层中目标深度位置处的孔隙度、各类矿物的含量；根据所述目标储层中各类矿物的含量计算得到第一体积模量和第一剪切模量；根据所述储层中目标深度位置处的孔隙度、所述第一体积模量和所述第一剪切模量，生成以孔隙纵横比为变量的第一目标函数；将所述第一目标函数取值最小时的孔隙纵横比作为目标孔隙纵横比；根据所述目标孔隙纵横比计算得到第二体积模量和第二剪切模量；根据所述第二体积模量、所述第二剪切模量计算得到所述频率的声波测井信号在所述目标储层中传播时的横波速度和/或纵波速度。2.根据权利要求1所述的方法，其特征在于，所述根据所述储层中目标深度位置处的孔隙度、所述第一体积模量和所述第一剪切模量，生成以孔隙纵横比为变量的第一目标函数，包括：根据孔隙纵横比、所述孔隙度、所述声波测井信号的频率，建立以孔隙纵横比为变量的预测第二体积模量的表达式和预测第二剪切模量的表达式；将所述预测第二体积模量的表达式和预测第二剪切模量的表达式组成的数据点作为预测数据点，将所述第一体积模量、所述第一剪切模量组成的数据点作为实测数据点，根据实测数据点与预测数据点之间的距离，生成以孔隙纵横比为变量的第一目标函数。3.一种基于井震融合的孔隙度确定方法，其特征在于，包括：获取目标储层的声波测井数据，所述声波测井数据包括：声波测井信号的频率、储层中目标深度位置处的孔隙度、各类矿物的含量；根据所述目标储层中各类矿物的含量计算得到第一体积模量和第一剪切模量；根据所述储层中目标深度位置处的孔隙度、所述第一体积模量和所述第一剪切模量，生成以孔隙纵横比为变量的第一目标函数；将所述第一目标函数取值最小时的孔隙纵横比作为目标孔隙纵横比；根据所述目标孔隙纵横比计算得到第二体积模量和第二剪切模量；根据所述第二体积模量、所述第二剪切模量计算得到所述频率的声波测井信号在所述目标储层中传播时的纵波速度；根据目标储层中各个深度位置的声波测井数据计算得到各个深度位置对应的孔隙纵横比；根据目标储层中多个深度位置对应的孔隙纵横比确定目标储层的孔隙纵横比的替代值；根据所述目标储层的孔隙纵横比的替代值，生成以孔隙度为变量的第二目标函数；将所述第二目标函数取值最小时的孔隙度作为目标孔隙度。4.根据权利要求3所述的方法，其特征在于，根据所述目标储层的孔隙纵横比的替代值，生成以孔隙度为变量的第二目标函数，包括：根据孔隙度、所述孔隙纵横比的替代值、所述声波测井信号的频率，建立以孔隙度为变量的预测第三体积模量的表达式和预测第三剪切模量的表达式；根据所述预测第三体积模量的表达式和所述预测第三剪切模量的表达式，计算纵波速
度的预测值；获取地震资料中纵波速度的实测值；根据所述纵波速度的预测值与实测值，生成以孔隙度为变量的第二目标函数。5.根据权利要求4所述的方法，其特征在于，在所述纵波速度的实测值与地震频段的测井信号对应，且声波测井信号的频段大于地震信号的频段的情况下，在所述根据所述目标储层的孔隙纵横比的替代值，生成以孔隙度为变量的第二目标函数；将所述第二目标函数取值最小时的孔隙度作为目标孔隙度之前，还包括：将声...

【专利技术属性】
技术研发人员：唐跟阳，韩旭，王尚旭，张郑玉成，
申请(专利权)人：中国石油大学北京，
类型：发明
国别省市：