一种基于支持矢量机算法的地震相分析方法技术

本发明专利技术提供了一种支持矢量机算法进行地震波形分类的方法。该方法首先从目的层的叠后地震数据体中选择待分类的N道地震数据；然后获取已知的能够表征目的层一段区域内地震信号横向变化的K个类别的样本数据；其中K≤N；然后利用已知的样本数据训练支持矢量机，以确定用于将样本数据分成K个类别的K个分类函数；之后利用训练好的K个分类函数对待分类的N道地震数据进行分类，根据分类结果形成所选岩石物理参数或者沉积相平面分布规律图。理参数或者沉积相平面分布规律图。理参数或者沉积相平面分布规律图。

【技术实现步骤摘要】
一种基于支持矢量机算法的地震相分析方法


[0001]本专利技术属于油气及煤层气地震勘探与开发，尤其是一种基于支持矢量机算法的地震相分析方法。

技术介绍

[0002]在石油、煤炭等地下沉积矿产的勘探开发中，沉积相研究具有极为重要的意义。然而，由于目的层深埋于地下，因此所采用的研究手段和研究方法与露头区的沉积相研究相比有很大不同。
[0003]在地下相分析中只有通过岩石资料才能够观察到目的的沉积相标志，而钻井取心一般都不是连续进行的，并且一口探井的全井取心率往往只有百分之几到百分之十几，这给沉积相研究造成很大困难。利用电测井资料进行测井相分析虽可对全井做出连续的沉积相解释，但其多解性较强，因此除上述两种资料外，还迫切需要从其它资料中获取更多的信息以提高沉积相解释的准确性。
[0004]更重要的是，即使单井相分析的资料足够充分，但采用传统研究方法所得到的毕竟只是一部分信息，而如地层叠置模式、沉积体外形等重要信息并没有利用。进一步看，即使解释完全正确，但毕竟只是“一孔之见”。要想进一步掌握沉积相的平面展布特征就必须有大量的足够密集的钻孔，而这在勘探阶段恰恰难以满足。因此迫切需要一种仅用少量钻孔就能较好地掌握沉积相平面变化特征的新手段、新方法。
[0005]地震相分析正是为满足上述迫切需要而产生的。地震相就是在地震反射时间剖面上所表现出来的反射波的面貌。地震相分析则是根据地震相特征进行沉积相的解释推断。在石油勘探及某些煤田、盐矿勘探中，地震勘探资料是必不可少的重要基础资料。这些资料一般在勘探初期就可获得，且一般都能覆盖整个盆地，其中具有极为丰富的地层、构造和沉积相信息，因此是地下地质分析中极为宝贵的基础资料。地震相分析作为地震地层学的一个重要组成部分，诞生于1977年左右，并在世界上迅速传播。十几年来它在广泛的实践中不断发展完善，已成为地下相分析不可缺少的锐利武器。
[0006]地震相分析的方法就是识别每个层序内独特的地震反射波组特征及其形态组合，并将其赋予一定的地质含义，进而进行沉积相的解释，这一过程称为地震相分析。
[0007]地震相分析和识别有两种方法，第一种方法是通过人工观测地震反射特征，并与所建立的标准地震相特征进行比较，判别属于何种地震相。这种方法一般应用于局部地震资料解释和分析中，解释和识别精度比较低。另外一种方法则是应用地震数据处理技术和计算机技术以及一定的数学方法对地震数据体或者地震属性数据进行分析和计算，提取出能够反映沉积相变化的地震相，这是一种高效、先进、定量的地震相识别方法。
[0008]地震波形是地震数据的基本性质，它包含了所有的定性和定量信息，如反射模式、相位、频率和振幅等信息，是地震信息的总体特征，其动态变化蕴含了丰富的内在信息，能够真实地反映地下结构的特征。波形分类法是最常采用的地震相分析方法，通过对地震信号波形进行分类，可以实现对地震相的划分。
[0009]传统的波形分类方法有K均值和自组织神经网络方法。它们都采用无监督的学习方法。
[0010]本专利技术以新的视角提出一种基于支持矢量机算法的地震相分析方法,该方法属于一种有监督的学习方法，为之后油气储量的评估提供科学依据。

技术实现思路

[0011]针对上述问题，本专利技术提出了一种基于支持矢量机算法的地震相分析方法，以及相应的计算机存储介质和计算机设备。
[0012]根据本专利技术的第一个方面，本专利技术的一种基于支持矢量机算法的地震相分析方法，主要包括以下步骤：
[0013]S100，从目的层的叠后地震数据体中选择待分类的N道地震数据；
[0014]S200，获取已知的能够表征目的层一段区域内地震信号横向变化的K个类别的样本数据，其中K≤N；
[0015]S300，利用已知的样本数据训练支持矢量机，以确定用于将样本数据分成K个类别的K个分类函数；
[0016]S400，利用训练好的K个分类函数对待分类的N道地震数据进行分类，以将N道地震数据划分成K个分区；
[0017]S500，基于分类结果分析目的层沉积相平面分布规律。
[0018]根据本专利技术的实施例，上述步骤S100包括以下步骤：
[0019]根据给定的时窗中心和时窗长度，为目的层的时间域或者深度域的叠后地震数据体设置数据时窗；
[0020]将数据时窗中的N道数据作为待分类的N道地震数据。
[0021]根据本专利技术的另一实施例，上述步骤S100包括以下步骤：
[0022]解释目的层的构造层位，将解释层位作为时窗中心；
[0023]根据给定的时窗中心和时窗长度，为目的层的时间域或者深度域的叠后地震数据体设置数据时窗；
[0024]将数据时窗中的N道数据作为待分类的N道地震数据。
[0025]根据本专利技术的实施例，上述步骤S200包括以下步骤：
[0026]获取目的层的波形分类时窗中心的地震属性切片；
[0027]在波形分类时窗中心的地震属性切片上找出能够表征地震信号横向变化的K个类别的区域；
[0028]记录下该区域内的数据点信息，并为这些数据点赋予相应的类别标签k，k＝1,2,
…
,K，从而获得分为K个类别的样本数据。
[0029]根据本专利技术的实施例，上述步骤S300中，所述K个分类函数中的第k个分类函数定义为：
[0030][0031]且
[0032]式中，y为支持矢量机算法的二元分类标记，W
kT
和b
k
为第k个分类函数的系数，为N道地震数据中的一个数据点的矢量分量。
[0033]根据本专利技术的实施例，上述步骤S300中，通过以下步骤确定用于将样本数据分成K个类别的K个分类函数中的第k个分类函数：
[0034]将样本数据中标签为k的数据与其他数据分开，并将标签为k的数据赋予标记y＝1以及将所有标签不为k的数据赋予标记y＝-1，以形成包含y＝1和y＝-1两种标记类型的数据的训练集；
[0035]基于支持矢量机算法，利用所述训练集训练第k个分类函数，以确定第k个分类函数的系数和b
k
，从而确定第k个分类函数
[0036]根据本专利技术的实施例，上述步骤S400包括以下步骤：
[0037]对于N道地震数据中的第i个数据点，将其代入训练好的各个分类函数中，以计算各个分类函数的函数值；
[0038]确定各个分类函数的函数值中的最大值，与该最大值相对应的k值所代表的类别为第i个数据点所属的类别。
[0039]根据本专利技术的实施例，上述步骤S500包括以下步骤：
[0040]根据分类结果形成平面上离散的地震相，对地震相图进行岩石物理参数或沉积相标定，根据标定的地震相图形成目的层岩石物理参数或者沉积相空间变化的定量图。
[0041]根据本专利技术的另一方面，本专利技术还提供一种计算机存储介质，其中存储有用于实现上述方法的计算机程序。
[0042]根据本专利技术的另一方面，本专利技术还提供一种本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于支持矢量机算法的地震相分析方法，其特征在于，包括以下步骤：S100，从目的层的叠后地震数据体中选择待分类的N道地震数据；S200，获取已知的能够表征目的层一段区域内地震信号横向变化的K个类别的样本数据，其中K≤N；S300，利用已知的样本数据训练支持矢量机，以确定用于将样本数据分成K个类别的K个分类函数；S400，利用训练好的K个分类函数对待分类的N道地震数据进行分类，以将N道地震数据划分成K个分区；S500，基于分类结果分析目的层沉积相平面分布规律。2.根据权利要求1所述的方法，其特征在于，所述步骤S100包括以下步骤：根据给定的时窗中心和时窗长度，为目的层的时间域或者深度域的叠后地震数据体设置数据时窗；将数据时窗中的N道数据作为待分类的N道地震数据。3.根据权利要求1或2所述的方法，其特征在于，所述步骤S100包括以下步骤：解释目的层的构造层位，将解释层位作为时窗中心；根据给定的时窗中心和时窗长度，为目的层的时间域或者深度域的叠后地震数据体设置数据时窗；将数据时窗中的N道数据作为待分类的N道地震数据。4.根据权利要求1或2所述的方法，其特征在于，所述步骤S200包括以下步骤：获取目的层的波形分类时窗中心的地震属性切片；在波形分类时窗中心的地震属性切片上找出能够表征地震信号横向变化的K个类别的区域；记录下该区域内的数据点信息，并为这些数据点赋予相应的类别标签k，k＝1,2,
…
,K，从而获得分为K个类别的样本数据。5.根据权利要求4所述的方法，其特征在于，所述步骤S300中，所述K个分类函数中的第k个分类函数定义为：且...

【专利技术属性】
技术研发人员：郑四连，刘百红，朱海伟，杨子兴，
申请(专利权)人：中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：