基于纹理分析的缝洞储集体形态描述方法技术

本发明专利技术公开了一种基于纹理分析的缝洞储集体形态描述方法。该方法包括：将三维地震数据体映射为三维灰度数据体；沿设定角度、以设定的距离差扫描所述三维灰度数据体，得到相应的灰度共生矩阵；提取相应的灰度共生矩阵的特征纹理属性，所述特征纹理属性包括均匀度、对比度、熵以及能量；采用聚类分析技术分析所述特征纹理属性，实现针对所述三维地震数据体的纹理属性分类；基于所述纹理属性分类描述缝洞空间形态分布。应用本发明专利技术，能够精细描述缝洞空间形态。

Texture analysis based method for reservoir morphology description

The invention discloses a shape description method based on texture analysis. The method includes: 3D seismic data mapping for three-dimensional gray data; along the set angle, scanning the three-dimensional data set with gray distance difference between the corresponding gray level co-occurrence matrix; texture feature attribute extraction corresponding gray co-occurrence matrix, the texture attributes including uniformity, contrast, entropy by using the cluster analysis technology and energy; the analysis of the characteristics of texture properties, for texture attribute classification the 3D seismic data; the texture attribute classification description of fracture cave distribution based on spatial form. The application of the invention can accurately describe the space shape of the seam hole.

【技术实现步骤摘要】
基于纹理分析的缝洞储集体形态描述方法
本专利技术涉及地震资料解释领域，更具体地，涉及一种基于纹理分析的缝洞储集体形态描述方法。
技术介绍
作为塔河油田主力油田的奥陶系碳酸盐岩缝洞型油藏，其缝洞储集体的空间形态、内部结构复杂多样。对缝洞体空间结构的预测及描述不清必然会影响到缝洞单元划分的准确程度和油水关系认识，也会影响油藏开发技术政策的制定，而这些问题已经成为缝洞型碳酸盐岩油藏开发的瓶颈问题，制约了油田的开发效果。塔河油田缝洞型油藏的预测和描述技术经历了二十余年的发展，已经形成了针对缝洞储集体有效实用的地震储层识别和预测方法技术系列，并在实际应用中见到效果，主要包括：建立了缝洞型储层地震识别模式；形成了缝洞型储层预测的地球物理方法系列；初步形成了一套针对碳酸盐岩缝洞体储层空间描述的技术方法和流程。但由于缝洞储层地球物理响应关系的极其复杂性，加之地震分辨率的限制，目前所取得的方法技术成果还不能满足对缝洞储集体高精度描述的需求。如何利用地震资料准确地描述缝洞地质体的边界及空间形态是缝洞型油藏开发的重中之重，因此，借助于图像处理技术开发了地震数据的纹理处理技术来探索缝洞储集体的空间结构精细描述
技术实现思路
本专利技术提出了一种用于探索缝洞储集体的空间结构精细描述的方法。根据本专利技术的一方面，提出了一种基于纹理分析的缝洞储集体形态描述方法，包括：将三维地震数据体映射为三维灰度数据体；沿设定角度、以设定的距离差扫描所述三维灰度数据体，得到相应的灰度共生矩阵；提取相应的灰度共生矩阵的特征纹理属性，所述特征纹理属性包括均匀度、对比度、熵以及能量；采用聚类分析技术分析所述特征纹理属性，实现针对所述三维地震数据体的纹理属性分类；基于所述纹理属性分类描述缝洞空间形态分布。本专利技术中，借助图像处理技术，将三维地震数据体映射为三维灰度数据体，扫描三维灰度数据体得到灰度共生矩阵，并通过分析灰度共生矩阵的特征纹理属性实现对不同缝洞形态及结构的精细描述，有利于准确反映储层纵横向分布及隔层发育特征。附图说明通过结合附图对本专利技术示例性实施方式进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施方式中，相同的参考标号通常代表相同部件。图1示出了根据本专利技术的一个实施例的缝洞储集体形态描述方法的流程图。图2示出了根据本专利技术的灰度映射示意图。图3(a)示出了某实际工区的地震剖面；图3(b)示出了应用本专利技术描述的缝洞空间形态；图3(c)示出了应用本专利技术描述的缝洞形态剖面图。具体实施方式下面将参照附图更详细地描述本专利技术的优选实施方式。虽然附图中显示了本专利技术的优选实施方式，然而应该理解，可以以各种形式实现本专利技术而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本专利技术更加透彻和完整，并且能够将本专利技术的范围完整地传达给本领域的技术人员。图1示出了根据本专利技术的一个实施例的……方法的流程图，该方法包括：步骤101，将三维地震数据体映射为三维灰度数据体；步骤102，沿设定角度、以设定的距离差扫描所述三维灰度数据体，得到相应的灰度共生矩阵；步骤103，提取相应的灰度共生矩阵的特征纹理属性，所述特征纹理属性包括均匀度、对比度、熵以及能量；步骤104，采用聚类分析技术分析所述特征纹理属性，实现针对所述三维地震数据体的纹理属性分类；步骤105，基于所述纹理属性分类描述缝洞空间形态分布。本实施例中，将三维地震数据体映射为三维灰度数据体，扫描三维灰度数据体得到灰度共生矩阵，并通过分析灰度共生矩阵的特征纹理属性实现对不同缝洞形态及结构的精细描述。在上述实施例中，本领域技术人员认为使用的任意聚类分析技术进行聚类分析，在聚类分析中，可对比三维地震数据体的特征纹理和地震数据模型道的特征纹理进行聚类距离分析，得到相应的距离分析结果。在一个示例中，将三维地震数据体映射为三维灰度数据体可以包括：针对所述三维地震数据体中任意数据点A(x,y,z)，可以基于下式得到映射后的灰度值g(x,y,z)：其中，(x,y,z)表示位置坐标，Amin为所述三维地震数据体中的最小数据值，Amax为所述三维地震数据体中的最大数据值，INT表示取整运算。图2示出了根据本专利技术的灰度映射示意图。在一个示例中，沿设定角度、以设定的距离差扫描所述三维灰度数据体可以包括：可以分别沿x轴、y轴和z轴、并且可以均以距离差为1进行扫描，其中，沿x轴、以距离差为1进行扫描，得到的相应灰度共生矩阵中第i行第j列的元素pi,j为：pi,j＝Count({(m,n,o),(p,q,r)||m-p|＝1,n-q＝0,o-r＝0,g(m,n,o)＝i,g(p,q,r)＝j})，沿y轴、以距离差为1进行扫描，得到的相应灰度共生矩阵中第i行第j列的元素pi,j为：pi,j＝Count({(m,n,o),(p,q,r)|m-p＝0,|n-q|＝1,o-r＝0,g(m,n,o)＝i,g(p,q,r)＝j})，沿z轴、以距离差为1进行扫描，得到的相应灰度共生矩阵中第i行第j列的元素pi,j为：pi,j＝Count({(m,n,o),(p,q,r)|m-p＝0,n-q＝0,|o-r|＝1,g(m,n,o)＝i,g(p,q,r)＝j})，其中，g(m,n,o)表示所述三维灰度数据体中位置坐标(m,n,o)处的数据点的灰度值，g(p,q,r)表示所述三维灰度数据体中位置坐标(p,q,r)处的数据点的灰度值Count表示在所述三维灰度数据体中符合相应条件的数据点对(m,n,o),(p,q,r)的数量。众所周知，扫描角度和扫描间隔非常重要，其直接影响了得到的灰度共生矩阵是否能准确有效地表达期望得到的纹理特征。专利技术人经过大量研究和实验后发现，按上述扫描方式扫描得到的灰度共生矩阵有利于准确反映能准确有效地描述缝洞形态和结构的纹理信息。在一个示例中，可以通过下式提取相应的灰度共生矩阵的均匀度：其中，pi,j表示相应的灰度共生矩阵中第i行第j列的元素，L表示设定的时窗长度。在一个示例中，可以通过下式提取相应的灰度共生矩阵的对比度：其中，pi,j表示相应的灰度共生矩阵中第i行第j列的元素，L表示设定的时窗长度。在一个示例中，可以通过下式提取相应的灰度共生矩阵的熵：其中，pi,j表示相应的灰度共生矩阵中第i行第j列的元素，L表示设定的时窗长度。在一个示例中，可以通过下式提取相应的灰度共生矩阵的能量：其中，pi,j表示相应的灰度共生矩阵中第i行第j列的元素，L表示设定的时窗长度。应用示例为便于理解本专利技术实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本专利技术，其任何具体细节并非意在以任何方式限制本专利技术。图3(a)示出了某实际工区的地震剖面；图3(b)示出了应用本专利技术描述的缝洞空间形态；图3(c)示出了应用本专利技术描述的缝洞形态剖面图。应用本专利技术的描述结果与钻测井解释结果具有较高的吻合度，在塔河油田缝洞体形态及空间结构描述中取得了良好的应用效果。本专利技术可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质，其上载有用于使处理器实现本专利技术的各个方面的计算机可读程序指令。计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储本文档来自技高网...

【技术保护点】
一种基于纹理分析的缝洞储集体形态描述方法，包括：将三维地震数据体映射为三维灰度数据体；沿设定角度、以设定的距离差扫描所述三维灰度数据体，得到相应的灰度共生矩阵；提取相应的灰度共生矩阵的特征纹理属性，所述特征纹理属性包括均匀度、对比度、熵以及能量；采用聚类分析技术分析所述特征纹理属性，实现针对所述三维地震数据体的纹理属性分类；基于所述纹理属性分类描述缝洞空间形态分布。

【技术特征摘要】
1.一种基于纹理分析的缝洞储集体形态描述方法，包括：将三维地震数据体映射为三维灰度数据体；沿设定角度、以设定的距离差扫描所述三维灰度数据体，得到相应的灰度共生矩阵；提取相应的灰度共生矩阵的特征纹理属性，所述特征纹理属性包括均匀度、对比度、熵以及能量；采用聚类分析技术分析所述特征纹理属性，实现针对所述三维地震数据体的纹理属性分类；基于所述纹理属性分类描述缝洞空间形态分布。2.根据权利要求1所述的方法，其中，将三维地震数据体映射为三维灰度数据体包括：针对所述三维地震数据体中任意数据点A(x,y,z)，基于下式得到映射后的灰度值g(x,y,z)：其中，(x,y,z)表示位置坐标，Amin为所述三维地震数据体中的最小数据值，Amax为所述三维地震数据体中的最大数据值，INT表示取整运算。3.根据权利要求1所述的方法，其中，沿设定角度、以设定的距离差扫描所述三维灰度数据体包括：分别沿x轴、y轴和z轴、并且均以距离差为1进行扫描，其中，沿x轴、以距离差为1进行扫描，得到的相应灰度共生矩阵中第i行第j列的元素pi,j为：pi,j＝Count({(m,n,o),(p,q,r)||m-p|＝1,n-q＝0,o-r＝0,g(m,n,o)＝i,g(p,q,r)＝j})，沿y轴、以距离差为1进行扫描，得到的相应灰度共生矩阵中第i行第j列的元素pi,j为：pi,j＝Count({(m,n,o),(p,q,r)|m-p＝0,|n-q|＝1,o-r＝0,g(m,n,o)＝i,g(p,q,r)＝j})，沿z轴、以距离差为1进行扫描，得到的相应灰度共生矩阵中第i行第j...

【专利技术属性】
技术研发人员：唐金良，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11