一种砂泥岩油气储层测井横波预测方法、装置及终端设备

技术领域

本发明属于岩石物理研究领域，尤其涉及一种砂泥岩油气储层测井横波预测方法、装置及终端设备。

背景技术

地震岩石物理方法在非常规油气藏勘探开发中起着重要的技术支持作用，能有效的对储层复杂性孔隙结构、物性和含流体性进行描述。在地震岩石物理研究过程中，发现采用岩石物理模型预测出的测井数据精度较高，符合地震地质规律，应用比较成熟的是Xu-White模型，该模型是基于砂泥岩建立的岩石物理模型，能够对砂泥岩的孔隙特征、矿物成分、流体等进行模拟，但是Xu-White模型，并没有考虑到实际地层中矿物成分对骨架支撑贡献的作用，以及没有考虑频率的影响。Xu and Payne(2009)基于Xu-White模型建立的碳酸盐岩岩石物理模型，该模型能够描述碳酸盐岩复杂的孔隙空间，将岩石中的孔隙考虑为三种孔隙，一种是圆孔，一种是微裂隙，另一种是介于前两者之间的孔隙。这些岩石物理模型能够在已知矿物成分、孔隙度、孔隙结构的情况下对岩石的速度等参数进行估算，在测井中常常用来估算横波速度值。

但是，目前已有的测井数据校正方法和横波速度预测一种是基于经验关系的，一种是基于复杂岩石物理模型来估算的。由于基于经验关系的测井校正方法经验关系的区域局限性导致经验模型的不实用性，而岩石物理模型方法则没有考虑频率对弹性的影响，而频散现象在测井尺度是不可忽略的。上述两种方法获得的结果都会导致预测测井曲线与实际地质信息不匹配，进一步影响井震标定和储层预测。因此，需要构建合适的基于岩石物理的砂泥岩储层测井横波预测方法，为砂泥岩储层预测和开发提供支撑作用。

发明内容

基于此，有必要针对上述技术问题，提供一种砂泥岩油气储层测井横波预测方法、装置及终端设备，能够有效和准确的获取测井的横波信息，同时获取手段简单、经济，为测井横波预测提供了有效的途径。

本发明第一方面提供了一种砂泥岩油气储层测井横波预测方法，所述方法包括：

结合预设岩石模型及岩石基质矿物参数，计算地层岩石模量；

基于预设岩石模型及所述地层岩石模量，并结合孔隙类型参数得到岩石骨架模量，所述孔隙类型参数表征为各个孔隙类型的数量占比，所述岩孔类型适宜通过孔隙纵横比确定；

利用测井测井参数及wood公式计算孔隙流体的混合流体参数；

结合所述岩石骨架模量及所述混合流体参数得到测井频段的纵波速度值及横波速度值；

判断原始纵波速度值与所述纵波速度值的差值是否满足预设条件，若是，则所述横波数据值即为横波预测值；否则，通过调节孔隙纵横比对所述岩石骨架模量进行调节，直至原始纵波速度值与新纵波速度值之间的差值满足预设条件，并将相应的新横波数据值视为横波预测值。

可选的，所述预设岩石模型包括：DEM模型；

其中，所述DEM模型为：

K

μ

式中，K

可选的，在所述结合预设岩石模型及岩石基质矿物参数，计算地层岩石模量之前，所述方法还包括：

根据测井GR计算得到地层黏土含量和石英含量。

可选的，所述结合预设岩石模型及岩石基质矿物参数，计算地层岩石模量，包括：

根据黏土含量和石英含量的占比情况，选出待输入所述DEM模型的相1及相2；

将选出的相1及相2输入所述DEM，得到所述地层岩石模量。

可选的，所述选出待输入所述DEM模型的相1及相2，包括：

在黏土含量大于60％的情况下，则将所述黏土作为待输入所述DEM模型的相1、以及石英作为待输入所述DEM模型的相2；

在石英含量大于60％的情况下，则将所述石英作为待输入所述DEM模型的相1、以及黏土作为待输入所述DEM模型的相2。

可选的，在基于预设岩石模型及所述地层岩石模量，并结合孔隙类型参数得到岩石骨架模量之前，所述方法包括：

基于孔隙纵横比对孔隙类型进行分类，并统计各个孔隙类型的数量占比。

可选的，基于预设岩石模型及所述地层岩石模量，并结合孔隙类型参数得到岩石骨架模量，包括：

基于所述各个孔隙类型参数，将所述地层岩石模量作为相1输入DEM模型，并先后将各个孔隙类型的孔隙作为相2输入DEM模型，得到岩石骨架模量，所述岩石骨架模量包括：包含孔隙度且不含流体的岩石骨架等效体积模量和岩石骨架剪切模量。

可选的，结合所述岩石骨架模量及所述混合流体参数得到测井频段的纵波速度值及横波速度值，包括：

将利用Dvorkin等1995提出的Mavko-Jizba射流公式全频带推广公式计算岩石模型中包含流体时岩石的纵波和横波速度值，计算公式:

/>

其中，

其中，

其中，

这里的W是角频率，J

其中，

其中，K

K

K

K

＝孔隙度；

μ

ρ＝岩石密度。

本发明第二方面提供了一种砂泥岩油气储层测井横波预测装置，所述装置包括：

地层岩石模量模块，用于结合预设岩石模型及岩石基质矿物参数，计算地层岩石模量；

岩石骨架模量模块，用于基于预设岩石模型及所述地层岩石模量，并结合孔隙类型参数得到岩石骨架模量，所述孔隙类型参数表征为各个孔隙类型的数量占比，所述岩孔类型适宜通过孔隙纵横比确定；

混合流体参数模块，用于利用测井测井参数及wood公式计算孔隙流体的混合流体参数；

纵横波速度值模块，用于结合所述岩石骨架模量及所述混合流体参数得到测井频段的纵波速度值及横波速度值；

校验模块模块，用于结合原始纵波速度值与所述纵波速度值的差值是否满足预设条件，若是，则所述横波数据值即为横波预测值；否则，通过调节孔隙纵横比对所述岩石骨架模量进行调节，直至原始纵波速度值与新纵波速度值之间的差值满足预设条件，并将相应的新横波数据值视为横波预测值。

本发明第三方面提供了一种终端设备，包括处理器和存储器；

所述存储器用于存储计算机指令，所述处理器用于运行所述存储器存储的计算机指令，以实现上述的砂泥岩油气储层测井横波预测方法。

本发明第四方面提供了一种计算机可读存储介质，所述计算机可读存储介质存储有一个或者多个程序，所述一个或者多个程序可被一个或者多个处理器执行，以实现上述的砂泥岩油气储层测井横波预测方法。

本发明有益效果如下：将岩石物理理论应用到砂泥岩储层测井的横波预测中，同时考虑频率的影响，实现了针对砂泥岩储层在测井频段上的基于岩石物理模型的横波预测方法，该方法能够有效和准确的获取测井的横波信息，同时获取方法简单、经济，为测井横波预测提供了有效的途径。

附图说明

图1为另一个实施例中一种砂泥岩油气储层测井横波预测方法的流程示意图；

图2为另一个实施例中一种砂泥岩油气储层测井横波预测方法的流程示意图；

图3为另一个实施例中横波预测实际效果图；

图4为另一个实施例中横波速度的误差值；

图5为另一个实施例中一种砂泥岩油气储层测井横波预测装置的结构框图；

图6为另一个实施例中终端设备的内部结构图。

具体实施方式

为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。

因为在砂泥岩储层的油气勘探中测试数据的质量好坏对整个油气勘探开发至关重要。尤其是横波测井数据信息的加入对地震反演及解释具有很大的推力作用。但是，在实际钻井中横波信息的获取难、成本高，因此，如何高质量的得到横波测井信息至关重要。同时从地震频段到实验室的超声频段，频率对弹性参数的影响一直认为是不可忽略的，因此以下实施例是利用理论与实验的结合来对测井频段上横波进行预测。

实施例1：

根据图1所示，本实施例提供了一种砂泥岩油气储层测井横波预测方法，该方法包括以下步骤：

步骤S201：结合预设岩石模型及岩石基质矿物参数，计算地层岩石模量；

步骤S202：基于所述预设岩石模型及所述地层岩石模量，并结合孔隙类型参数得到岩石骨架模量，所述孔隙类型参数表征为各个孔隙类型的数量占比，所述岩孔类型适宜通过孔隙纵横比确定；

步骤S203：利用测井测井参数及wood公式计算孔隙流体的混合流体参数；

步骤S204：结合所述岩石骨架模量及所述混合流体参数得到测井频段的纵波速度值及横波速度值；

步骤S205：判断原始纵波速度值与所述纵波速度值的差值是否满足预设条件，若是，则所述横波数据值即为横波预测值；否则，通过调节孔隙纵横比对所述岩石骨架模量进行调节，直至原始纵波速度值与新纵波速度值之间的差值满足预设条件，并将相应的新横波数据值视为横波预测值。

就此，将岩石物理理论应用到砂泥岩储层测井的横波预测中，同时考虑频率的影响，实现了针对砂泥岩储层在测井频段上的基于岩石物理模型的横波预测方法，该方法能够有效和准确的获取测井的横波信息，同时获取方法简单、经济，为测井横波预测提供了有效的途径。

实施例2：

具体的，根据图1-4所示，本实施例提供的一种砂泥岩油气储层测井横波预测方法中，步骤S201并不对该预设岩石模型进行限定，如：该预设岩石模型包括但不限于DEM模型；也不对岩石机制矿物参数的来源进行限定，如：根据测井GR(伽马曲线)计算得到地层黏土(泥质)含量和石英(砂岩)含量，当然，也可根据实验数据或已知数据得到地层黏土(泥质)含量和石英(砂岩)含量。

而且，其中，所述DEM模型为：

K

μ

式中，K

而且，上述步骤S201包括：根据黏土含量和石英含量的占比情况，选出待输入所述DEM模型的相1及相2；将选出的相1及相2输入所述DEM模型，得到所述地层岩石模量，其中，所述地层岩石模量包括：包含泥质的等效体积模量和剪切模量。

具体的，在获知地层的黏土(泥质)含量和石英(砂岩)含量后，若黏土含量大于60％，则将所述黏土作为待输入所述DEM模型的相1、以及石英作为待输入所述DEM模型的相2；若石英含量大于60％，则将所述石英作为待输入所述DEM模型的相1、以及黏土作为待输入所述DEM模型的相2。

即：在确定黏土含量大于60％的情况下，则认为黏土作为岩石基质矿物，并将黏土的体积模量及剪切模量作为待输入DEM模型的初始主相相1的体积模量K

在确定石英含量大于60％的情况下，则认为石英作为岩石基质矿物，并将石英的体积模量及剪切模量作为待输入DEM模型的初始主相相1的体积模量K

具体的，在通过测井GR(伽马曲线)计算得到地层黏土(泥质)含量和石英(砂岩)含量后，确定岩石基质矿物后(如：黏土含量大于60％则认为黏土作为基质矿物，否则石英作为基质矿物)，然后利用DEM模型公式计算(Berryman,1992)岩石基质的等效模量。如：在黏土作为岩石基质矿物时，假设泥质为针状包含物形状，因此等效体积模量P和剪切模量Q的计算公式为：

其中，下标m和i分别指背景材料和包含物材料(Berryman,1995)。

故利用上述公式计算得到黏土作为岩石基质矿物时的地层岩石模量，即：黏土作为岩石基质矿物时地层岩石基质的等效体积模量和剪切模量。

在得到上述的步骤S201可以得到地层岩石模量。因为接下来的步骤S202为：基于所述预设岩石模型及所述地层岩石模量，并结合孔隙类型参数得到岩石骨架模量。

其中，所述孔隙类型参数表征为各个孔隙类型的数量占比，所述岩孔类型适宜通过孔隙纵横比确定；

故针对该孔隙类型，在本实施例中，可预先基于孔隙纵横比对孔隙类型进行分类，并统计各个孔隙类型的数量占比。如：在步骤201与步骤202之间执行该步骤。而且，该孔隙包括但不限于：圆孔、裂缝。而且，针对确定地层目的层孔隙类型，可以结合实际岩心铸体薄片等微观结构成像实验结果得到，此外，砂泥岩地层中主要以近圆形孔的砂岩孔，及近针状的泥质孔为主，利用孔隙纵横比(短轴与长轴之比)可以用来表征孔隙类型。如：孔隙纵横比大于等于0.01则为圆孔，否则，则为裂缝。

在本实施例中，上述步骤S202具体包括：基于所述各个孔隙类型参数，将所述地层岩石模量作为相1输入DEM模型，并先后将各个孔隙类型的孔隙作为相2输入DEM模型，得到岩石骨架模量，所述岩石骨架模量包括：包含孔隙度且不含流体的岩石骨架等效体积模量和岩石骨架剪切模量。

具体的，在经步骤S201得到岩石基质的体积模量和剪切模量的情况下，再利用DEM模型公式计算双孔岩石模型的等效弹性模量，即：加入两种孔隙(圆孔、裂缝)。先加入圆孔，后加入裂缝孔隙。

其中，DEM模型的表达式为：

K

μ

式中，K

在这里计算包含圆孔时，将圆孔等效为球体形状加入。具体的，将经步骤S201得到的岩石基质的剪切模量和体积模量作为初始主相1的剪切模量μ

接着，在计算包含裂缝时，将裂缝等效为硬币状孔隙，通过将圆孔加入后计算所得的体积模量和剪切模量作为初始主相1的剪切模量μ

公式中α为孔隙的纵横比值：0.01，下标m和i分别指背景材料和包含物材料(Berryman，1995)。即：得到包含孔隙度，且不含流体的等效体积模量

当然，在本实施例中并不对加入圆孔和裂缝的先后顺序进行限定，在方案实施时可以根据实际需要适应性选择。如：先加入裂缝，后加入圆孔。故在这里计算包含裂缝时，将裂缝等效为硬币状孔隙加入。具体的，将经步骤S201得到的岩石基质的剪切模量和体积模量作为初始主相1的剪切模量μ

此外，还需进行步骤S203，及：利用测井测井参数及wood公式计算孔隙流体的混合流体参数；

即：确定砂泥岩地层的流体类型及流体等效岩石物理参数计算，可利用测井曲线中的含水饱和度或者含油饱和度及wood公式计算孔隙流体的岩石物理参数；

此外，在通过上述步骤S202得到包含孔隙度，且不含流体的等效体积模量

K

具体的，该步骤S204为：结合所述岩石骨架模量及所述混合流体参数得到测井频段的纵波速度值及横波速度值，包括：

将利用Dvorkin等1995提出的Mavko-Jizba射流公式全频带推广公式计算岩石模型中包含流体时岩石的纵波和横波速度值，计算公式:

其中，

其中，

其中，

这里的W是角频率，J

其中，

其中，K

K

K

K

＝孔隙度；

μ

ρ＝岩石密度。

经步骤S204得到横波速度值及纵波速度值的情况下，可通过步骤S205得到横波预测值，具体的，该步骤S205包括：判断原始纵波速度值与所述纵波速度值的差值是否满足预设条件，若是，则所述横波数据值即为横波预测值；否则，通过调节孔隙纵横比对所述岩石骨架模量进行调节，直至原始纵波速度值与新纵波速度值之间的差值满足预设条件，并将相应的新横波数据值视为横波预测值。

具体的，在上述步骤S204的计算过程中，首先要找到Z值，初始Z＝0.001，通过计算测井频率(10khz)下的V

更具体的，在该步骤S205中，通过比较原始纵波速度值与所述纵波速度值的差值是否满足预设条件，若是，则所述横波数据值即为横波预测值；否则，则返回上述步骤S202中，通过调节孔隙纵横比对所述岩石骨架模量进行调节，在此基础上再次执行步骤S203及步骤S204，以获得新的纵波速度值和横波速度值，再次比较该新的纵波速度值与原始纵波速度值之间的差值是否满足预设条件，以此循环执行，直至原始纵波速度值与新纵波速度值之间的差值满足预设条件，并将相应的新横波数据值视为横波预测值。

就此，将岩石物理理论应用到砂泥岩储层测井的横波预测中，同时考虑频率的影响，实现了针对砂泥岩储层在测井频段上的基于岩石物理模型的横波预测方法，该方法能够有效和准确的获取测井的横波信息，同时获取方法简单、经济，为测井横波预测提供了有效的途径。

在地震储层预测中横波信息至关重要，横波速度是叠前反演的输入数据，横波数据的准确性直接关系到储层预测及流体识别的精度。然而实际常规测井中没有横波信息，需要横波测井才能获取到横波信息，而获取横波信息的测井手段较常规测井难、成本高，如何快速经济的得到准确的横波是测井研究中的关键问题之一。现认为依靠岩石物理方法的横波预测是相对更符合实际岩石的，因此本专利将岩石物理理论应用到砂泥岩储层测井的横波预测中，同时考虑频率的影响，实现了针对砂泥岩储层在测井频段上的基于岩石物理模型的横波预测方法，该方法能够有效和准确的获取测井的横波信息，同时获取方法简单、经济，为测井横波预测提供了有效的途径。

本实施例构建了一个适用于砂泥岩油气储层在测井频段上横波及其他弹性参数的预测方法。第一步：根据测井GR(伽马曲线)计算得到地层黏土(泥质)含量和石英(砂岩)含量，如果黏土含量大于60％，则认为地层岩石基质矿物为黏土(泥质)，如果石英(砂岩)含量大于60％则认为岩石基质矿物为石英，确定了岩石基质矿物后则利用数据公式计算地层基质的岩石物理参数，可采用DEM模型计算；

第二步：确定地层目的层孔隙类型，结合实际岩心铸体薄片等微观结构成像实验结果得到，砂泥岩地层中主要以近圆形孔的砂岩孔，及近针状的泥质孔为主，利用孔隙纵横比(短轴与长轴之比)来表征孔隙类型，并采用DEM模型将孔隙加入岩石基质中，计算得到岩石骨架的岩石物理等效参数；

第三步：确定砂泥岩地层的流体类型及流体等效岩石物理参数计算，可利用测井曲线中的含水饱和度或者含油饱和度及wood公式计算孔隙流体的岩石物理参数；

第四步：将第三步的流体加入岩石骨架的孔隙中，这里采用Mavko和Jizba射流公式来加入并计算得到测井频段包含流体的岩石的等效岩石物理参数，从而可计算获得纵波速度和横波速度，并调节孔隙纵横比来进行调节预测纵波和横波速度的计算值，并将测井原始纵波速度值作为模型预测结果多解性的约束，得到可靠的横波速度，实现测井的横波预测。具体步骤及涉及的计算公式如下：

第一步，确定了岩石基质矿物后(黏土含量大于60％则认为黏土作为基质矿物，否则石英作为基质矿物)。然后利用DEM模型公式计算(Berryman,1992)岩石基质的等效模量,DEM的表达式为：

K

μ

式中，K

在这里计算包含泥质时，假设泥质为针状包含物形状，因此P、Q的计算公式为：

下标m和i分别指背景材料和包含物材料(Berryman,1995)。

利用上述公式计算得到包含泥质的等效体积模量和剪切模量。

(2)在第(1)步的基础上再利用DEM模型公式计算双孔岩石模型的等效弹性模量，即：加入两种孔隙(圆孔、裂缝)。先加入圆孔，后加入裂缝孔隙。DEM的表达式为：

K

μ

式中，K

在这里计算包含圆孔时，将圆孔等效为球体形状加入，公式中的P、Q的计算公式为：

计算包含裂缝时，将裂缝等效为硬币状孔隙，公式中的P和Q的计算公式为：

公式中α为孔隙的纵横比值：0.01，下标m和i分别指背景材料和包含物材料(Berryman，1995)。

利用上述公式计算得到包含孔隙度，且不含流体的等效体积模量

(3)再利用Dvorkin等1995提出的Mavko-Jizba射流公式全频带推广公式计算岩石模型中包含流体时岩石的纵波和横波速度值，计算公式：

其中，

其中，

其中，

这里的W是角频率，J

其中，

其中，K

K

K

K

＝孔隙度；

μ

ρ＝岩石密度；

在上述计算过程中，首先要找到Z值，初始Z＝0.001，通过计算测井频率(10khz)下的V

则认为吻合。此时得到的V

采用本发明提出的方法对实际砂泥岩地层的A井开展了储层横波预测。图2给出的是实际实现的流程图，图3给出的是A井预测横波、预测纵波和原始纵波、原始横波曲线对比图，从图中可看出预测纵波和横波曲线数据均吻合较好，说明了改方法的可行性。与此同时，图4中给出了预测横波速度与原始横波速度的误差曲线图，误差公式为：

从图2-4中可看出，误差范围整体较小(小于0.25),说明了在预测过程中方法的可靠性和稳定性。流程方法简单，在岩石物理物理模型计算中，采用了Dvorkin,Mavko和Nur(1995)在Mavko-Jizba射流公式基础上进行了全频带的推广方法，实现对测井频率弹性参数的预测，不仅解决了实际勘探中横波资料的缺乏问题，同时消除了常规横波预测预测值频率与实际测井频率不符合的问题。

在另一个实施例中，如图5所示，提供了一种砂泥岩油气储层测井横波预测的装置，所述装置包括：

地层岩石模量模块101，用于结合预设岩石模型及岩石基质矿物参数，计算地层岩石模量；

岩石骨架模量模块102，用于基于预设岩石模型及所述地层岩石模量，并结合孔隙类型参数得到岩石骨架模量，所述孔隙类型参数表征为各个孔隙类型的数量占比，所述岩孔类型适宜通过孔隙纵横比确定；

混合流体参数模块103，用于利用测井测井参数及wood公式计算孔隙流体的混合流体参数；

纵横波速度值模块104，用于结合所述岩石骨架模量及所述混合流体参数得到测井频段的纵波速度值及横波速度值；

校验模块模块105，用于结合原始纵波速度值与所述纵波速度值的差值是否满足预设条件，若是，则所述横波数据值即为横波预测值；否则，通过调节孔隙纵横比对所述岩石骨架模量进行调节，直至原始纵波速度值与新纵波速度值之间的差值满足预设条件，并将相应的新横波数据值视为横波预测值。

可选的，所述预设岩石模型包括：DEM模型；

其中，所述DEM模型为：

/>

K

μ

式中，K

可选的，所述装置还包括：测量模块，用于在所述结合预设岩石模型及岩石基质矿物参数，计算地层岩石模量之前，根据测井GR计算得到地层黏土含量和石英含量。

可选的，所述地层岩石模量模块101具体用于：根据黏土含量和石英含量的占比情况，选出待输入所述DEM模型的相1及相2；将选出的相1及相2输入所述DEM模型，得到所述地层岩石模量。

可选的，所述地层岩石模量模块101还用于：在黏土含量大于60％的情况下，则将所述黏土作为待输入所述DEM模型的相1、以及石英作为待输入所述DEM模型的相2；在石英含量大于60％的情况下，则将所述石英作为待输入所述DEM模型的相1、以及黏土作为待输入所述DEM模型的相2。

可选的，所述装置还包括：分类模块，用于在基于预设岩石模型及所述地层岩石模量，并结合孔隙类型参数得到岩石骨架模量之前，基于孔隙纵横比对孔隙类型进行分类，并统计各个孔隙类型的数量占比。

可选的，该岩石骨架模量模块102具体用于：基于所述各个孔隙类型参数，将所述地层岩石模量作为相1输入DEM模型，并先后将各个孔隙类型的孔隙作为相2输入DEM模型，得到岩石骨架模量，所述岩石骨架模量包括：包含孔隙度且不含流体的岩石骨架等效体积模量和岩石骨架剪切模量。

可选的，该校验模块模块具体用于：将利用Dvorkin等1995提出的Mavko-Jizba射流公式全频带推广公式计算岩石模型中包含流体时岩石的纵波和横波速度值，计算公式:

其中，

/>

其中，

其中，

这里的W是角频率，J

其中，

其中，K

K

K

K

＝孔隙度；

μ

ρ＝岩石密度。

在上述计算过程中，首先要找到Z值，初始Z＝0.001，通过计算测井频率(10khz)下的V

关于上述装置的具体限定可以参见上文中对于相关方法的限定，在此不再赘述。上述装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于终端设备中的处理器中，也可以以软件形式存储于终端设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。

在另一个实施例中，提供了一种终端设备，该终端设备可以是服务器，其内部结构图可以如图6所示。该终端设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中，该终端设备的处理器用于提供计算和控制能力。该终端设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该终端设备的数据库用于存储相关数据。该终端设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种砂泥岩油气储层测井横波预测方法。

本领域技术人员可以理解，图6中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的终端设备的限定，具体的终端设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。

在另一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述的一种砂泥岩油气储层测井横波预测方法。

本实施例中的一种计算机可读存储介质所涉及的名词及实现原理具体可以参照本发明实施例中的的一种砂泥岩油气储层测井横波预测方法，在此不再赘述。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。

以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。