井控稀疏表征大规模谱聚类地震相划分方法技术

本发明专利技术公开了一种井控稀疏表征大规模谱聚类地震相划分方法，包括以下步骤：S1、输入测井数据和测井数据所对应的岩相标签；S2、输入地震数据所得属性体，对该数据按照地震数据采样方式进行网格化，将网格化坐标进行随机采样，得到k个采样点及其对应的属性曲线；S3、以测井数据作为标签数据，采用支持向量机或神经网络方法训练分类器；S4、以随机选择点和测井数据点作为地标点，进行基于地标点的谱聚类操作，可得到井控谱聚类地震相划分结果。本发明专利技术在基于地标点的大规模数据谱聚类方法中，通过加入测井(钻井)数据作为地标点，可以明确谱聚类方法应用于地震相划分的地质含义，提高地震相划分精度。

【技术实现步骤摘要】
井控稀疏表征大规模谱聚类地震相划分方法
本专利技术属于石油勘探
，特别涉及一种井控稀疏表征大规模谱聚类地震相划分方法。
技术介绍
油藏描述的核心技术是随机建模。随机建模就是利用计算机，以地质统计学为基础，结合三维可视化技术，将储层的空间结构真实的展示出来。由于测井数据、地震资料等原始数据的局限性，使得对地下的复杂结构的认知具有一定程度的不确定性。传统的地质建模技术以测井数据为核心，辅以地震资料、层位解释数据等，在建立合适的变差函数之后使用序贯模拟、指示模拟等方法建立随机模型。若在此过程中运用地震相分析技术提前对区域进行分类，在不同地震相区域上单独建模，并且充分利用测井数据和地震资料，就可以更好的反映真实的地质情况。谱聚类方法是近年来一种流行的机器学习聚类方法，对数据可以进行有效的聚类。谱聚类方法应用于对地震相的分类，可以取得与实际地震资料吻合的结果。因此有足够理由相信，谱聚类方法将成为地震相划分技术中一种重要、有效的方法。但是地震相划分中如果直接采用谱聚类方法不具有直观的地质意义，难以解释。
技术实现思路
本专利技术的目的在于克服现有技术的不足，提供一种在基于地标点的大规模数据谱聚类方法中，通过加入测井(钻井)数据作为地标点，可以明确谱聚类方法应用于地震相划分的地质含义，提高地震相划分精度的井控稀疏表征大规模谱聚类地震相划分方法。本专利技术的目的是通过以下技术方案来实现的：井控稀疏表征大规模谱聚类地震相划分方法，包括以下步骤：S1、输入测井数据和测井数据所对应的岩相标签；S2、输入地震数据所得属性体，对该数据按照地震数据采样方式进行网格化，将网格化坐标进行随机采样，得到k个采样点及其对应的属性曲线，随机采样点不包含井点所在位置；S3、以测井数据作为标签数据，采用支持向量机或神经网络方法训练分类器，并将该分类器用于对前述步骤随机选择点位置提取出的曲线进行分类操作；S4、以随机选择点和测井数据点作为地标点，进行基于地标点的谱聚类操作，可得到井控谱聚类地震相划分结果。进一步地，所述步骤S4具体实现方法为：选择一定数量代表性数据点作为地标，并将原始数据点表示为这些地标的稀疏线性组合；然后利用基于地标的稀疏表示来高效地计算谱聚类中的特征向量从而加速谱聚类；具体包括以下子步骤：S41、设计具有如下属性的邻接矩阵：W＝ZTZ，其中Z∈Rp×n是稀疏的且p远小于n；S42、导出Z，在O(p3+p2n)中计算图拉普拉斯矩阵的特征向量；Z的每一列都是原始输入相对于新基向量组U的p维表示；因此，W＝ZTZ是用新的表示法表示的数据X＝UZ的相似矩阵；S43、对于任一数据点,近似表示成其中uj是U矩阵的第j个列向量，zji是矩阵Z的第j行i列的元素；用下式计算zji,上式中K是一个计算相似度的核函数；S44、通过SVD方法分解矩阵：Z＝A∑BT，中间矩阵为p乘p对角矩阵，其组成元素为被分解矩阵的p个奇异值，即特征值的算术平方根；A是p*p的左奇异矩阵，其列向量称为左奇异向量，B是n*p的右奇异矩阵，其列向量为右奇异向量；根据数学定理可知B的列向量是W的特征向量；A的列向量是ZZT的特征向量；ZZT的维度是p乘p，计算A只需O(p3)的时间；有了A之后，B通过下式被计算出来：BT＝∑-1ATZ通过B求解W的特征向量。本专利技术的有益效果是：本专利技术在基于地标点的大规模数据谱聚类方法中，通过加入测井(钻井)数据作为地标点，可以明确谱聚类方法应用于地震相划分的地质含义，提高地震相划分精度。所提方案充分利用了已知实际测井数据，选取的地标点更具真实性和代表性。分类结果不但与实际地震资料吻合，分辨率也明显高于随机选取地标点和k-means聚类中心作为地标点的方法。附图说明图1为本专利技术的井控稀疏表征大规模谱聚类地震相划分方法的流程图。具体实施方式下面结合附图进一步说明本专利技术的技术方案。谱聚类算法是一种重要的无监督聚类算法。相比传统的k-means算法，谱聚类对数据分布的适应性更强，适合处理非凸线性不可分问题以及稀疏和高维数据。谱聚类算法是近年来聚类学习领域的一大热点，其以谱图理论为基础，将聚类问题转化为图的最优划分问题。与传统聚类方法相比，谱聚类能对任意形状的数据样本进行最优聚类的特点使其近年来应用广泛。主要思想是：将数据看作空间中的点，点之间可以连边，边的权重有高有低。对所有数据点组成的图进行切图，要求切图原则为尽可能让不同子图(由原图点的集合的子集和边集的子集构成的图)间边权重和尽可能低，子图内边权重和尽可能高，每一个子图代表一类。下面为常规谱聚类算法步骤：输入：数据集合D，构造邻接矩阵，降维后维度k，聚类个数k输出：聚类结果C(1)构建相似矩阵S、邻接矩阵W和度矩阵D；(2)计算拉普拉斯矩阵L；(3)计算L的最小的k个特征值及对应的特征向量矩阵A；(4)对矩阵A进行标准化得到B；(5)对标准化后的矩阵B进行k均值聚类，得到聚类结果C。如图1所示，本专利技术的一种井控稀疏表征大规模谱聚类地震相划分方法，包括以下步骤：S1、输入测井数据和测井数据所对应的岩相标签；S2、输入地震数据所得属性体，对该数据按照地震数据采样方式进行网格化，将网格化坐标进行随机采样，得到k个采样点及其对应的属性曲线，随机采样点不包含井点所在位置；S3、以测井数据作为标签数据，采用支持向量机或神经网络方法训练分类器，并将该分类器用于对前述步骤随机选择点位置提取出的曲线进行分类操作；S4、以随机选择点和测井数据点作为地标点，进行基于地标点的谱聚类操作，可得到井控谱聚类地震相划分结果。所述步骤S4具体实现方法为：选择一定数量代表性数据点作为地标，并将原始数据点表示为这些地标的稀疏线性组合；然后利用基于地标的稀疏表示来高效地计算谱聚类中的特征向量从而加速谱聚类；具体包括以下子步骤：S41、设计具有如下属性的邻接矩阵：W＝ZTZ，其中Z∈Rp×n是稀疏的且p远小于n；S42、导出Z，在O(p3+p2n)中计算图拉普拉斯矩阵的特征向量，而无需将整个相似性矩阵明确地存储在主存储器中，这节省了大量的计算量和存储空间。观察可知Z的每一列都是原始输入相对于新基向量组U的p维表示。因此，W＝ZTZ是用新的表示法表示的数据(X＝UZ)的相似矩阵，并且自然可以被用作邻接矩阵来加速谱聚类。然而计算Z和U十分耗时，特别X矩阵十分稠密时更是如此。为了解决这个问题(如何避免迭代求U和Z)，我们注意到基向量(U的列向量)与原始数据点具有相同的维数m。将基向量作为数据集的地标点(能表征原始数据集特征的点)，即地标点作为基向量，避免了迭代。从数据集中选择地标点的最有效方法是随机抽样。除了随机选择之外，还可以使用k-means算法首先对所有数据点进行聚类，然后使用聚类中心作为地标点。假设已经有了地标矩阵U，可以来计算矩阵Z.通过修正U，优化问题就变成了一个稀疏约束线性回归问题。有很多算法可以解决这个问题，此处可以简单地使用Nadaraya-Watson核回归来计算表示矩阵Z。S43、对于任一数据点,近似表示成其中uj是U矩阵的第j个列向量，zji是矩阵Z的第j行i列的元素；一个自然地假设是如果数据点Xi和地标点uj离得近的话zji应该大一点。我们可以强化这一假设，如果地标点uj不是离数据点xi的最近的r(远小于本文档来自技高网...

【技术保护点】
1.井控稀疏表征大规模谱聚类地震相划分方法，其特征在于，包括以下步骤：S1、输入测井数据和测井数据所对应的岩相标签；S2、输入地震数据所得属性体，对该数据按照地震数据采样方式进行网格化，将网格化坐标进行随机采样，得到k个采样点及其对应的属性曲线，随机采样点不包含井点所在位置；S3、以测井数据作为标签数据，采用支持向量机或神经网络方法训练分类器，并将该分类器用于对前述步骤随机选择点位置提取出的曲线进行分类操作；S4、以随机选择点和测井数据点作为地标点，进行基于地标点的谱聚类操作，可得到井控谱聚类地震相划分结果。

【技术特征摘要】
1.井控稀疏表征大规模谱聚类地震相划分方法，其特征在于，包括以下步骤：S1、输入测井数据和测井数据所对应的岩相标签；S2、输入地震数据所得属性体，对该数据按照地震数据采样方式进行网格化，将网格化坐标进行随机采样，得到k个采样点及其对应的属性曲线，随机采样点不包含井点所在位置；S3、以测井数据作为标签数据，采用支持向量机或神经网络方法训练分类器，并将该分类器用于对前述步骤随机选择点位置提取出的曲线进行分类操作；S4、以随机选择点和测井数据点作为地标点，进行基于地标点的谱聚类操作，可得到井控谱聚类地震相划分结果。2.根据权利要求1所述的井控稀疏表征大规模谱聚类地震相划分方法，其特征在于，所述步骤S4具体实现方法为：选择一定数量代表性数据点作为地标，并将原始数据点表示为这些地标的稀疏线性组合；然后利用基于地标的稀疏表示来高效地计算谱聚类中的特征向量从而加速谱聚类；具体包括以下子步骤：S41、设计具有如下属性...

【专利技术属性】
技术研发人员：王峣钧，李昊阳，胡光岷，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川,51