利用地层叠前纹理属性值来检测缝洞的方法技术

本发明专利技术提供一种利用地层叠前纹理属性值来检测缝洞的方法，所述方法包括：(A)分别读取目标工区中的各个点在不同方位的三维地震数据体中所对应的地震数据；(B)获取目标工区的待处理梯度体结构张量方阵；(C)从待处理梯度体结构张量方阵的每个元素中抽取与目标工区的一个点的坐标位置相对应的位置处的数据，构成一个点的结构张量方阵；(D)确定一个点的结构张量方阵的特征值；(E)获取一个点的叠前纹理属性值；(F)利用目标工区所有点的叠前纹理属性值来产生目标工区的缝洞检测结果。根据所述方法，能够有效地检测具有较小尺度的缝洞。

【技术实现步骤摘要】

本专利技术总体来说涉及油气田勘探领域，更具体地说，涉及一种利用地层叠前纹理属性值来检测缝洞的方法。
技术介绍
据统计，在全球的沉积岩中，碳酸盐岩尽管只占20％左右，但却拥有已探明油气资源的50％以上。碳酸盐岩储层储集类型通常由孔隙型、溶洞型、裂缝型、溶洞-裂缝型综合构成，这些裂缝和孔、洞系统对油气的赋存和运移起着控制作用，因此，从一定意义上说，对缝洞进行准确预测就等于直接在此类储层中找到油气，从而可确定地下地质储层的真实状态并较为精确地确定勘探目标。缝洞具有多尺度性，但地震勘探只有几米至数十米的分辨率，除了大的缝和洞之外，绝大多数单个的缝、洞无法利用地震勘探方法进行分辨和识别，但是，由众多细小的缝洞系统或缝洞却能被检测到。目前国内外常用的缝洞解释技术包括缝洞正演模拟技术、多波多分量地震技术、纵波裂缝检测方法各向异性技术、地震属性分析技术及三维可视化技术等。实际应用中，主要采用一些常规的技术手段，如横波分裂法、相干、地应力、倾角、倾向、曲率、边缘检测等。纹理属性来源于图像处理技术，所谓的纹理是指二维空间变化的灰度和颜色所组成的图案，它是图像区域所具有的重要特征之一。而地震剖面、切片甚至三维数据体本身也是图像，所以，从理论上来说，可利用纹理属性来分析地震数据纹理属性来描述地震数据的结构特征，并且在断裂、河道、乱岗状结构、平行与亚平行结构、倾斜层理、波状层理等地震相的描述中有明确的物理意义和良好的应用效果。Chopra(乔普拉)等人就成功地应用纹理属性来识别断层、河道以及地震相划分，并且，纹理属性应用于缝洞检测并获得成功的实例也屡见不鲜。当前，纹理属性的计算都是针对常规的叠后地震数据来处理的，这种计算方法以叠后地震数据体为基础输入，一定程度上能够预测地层中尺度较大、数量较少的缝洞，但是缺点在于：叠后纹理属性值是基于空间各点与周围点的反射振幅差异形成的，而事实上该差异有可能是岩性横向变化带来的，且叠后地震资料的信息量比较小，缺乏偏移距信息和方位角信息，这样不利于基于各向异性来检测较小尺度缝洞。因此，现有的利用纹理属性值来检测缝洞的方法在检测不同类型的缝洞方面具有局限性。
技术实现思路
本专利技术的目的在于提供一种利用地层叠前纹理属性值来检测缝洞的方法，以克服现有技术中无法检测小尺度的缝洞的缺陷。根据本专利技术的示例性实施例，提供一种利用地层叠前纹理属性值来检测缝洞的方法，所述方法包括：(A)分别读取目标工区中的各个点在不同方位的三维地震数据体中所对应的地震数据，其中，所述各个点在每个方位的三维地震数据体中的坐标位置相同；(B)根据读取的地震数据来获取目标工区的待处理梯度体结构张量方阵；(C)从待处理梯度体结构张量方阵的每个元素中抽取与目标工区的一个点的坐标位置相对应的位置处的数据，并将抽取的数据按所对应的元素在待处理梯度体结构张量方阵中的位置排列，以构成所述一个点的结构张量方阵；(D)确定所述一个点的结构张量方阵的特征值；(E)基于所述特征值，获取所述一个点的叠前纹理属性值；(F)利用目标工区所有点的叠前纹理属性值来产生目标工区的缝洞检测结果，其中，通过重复步骤(C)-(E)来获取目标工区所有点的叠前纹理属性值。可选地，步骤(B)可包括：(B1)分别确定目标工区中的各个点在不同方位的三维地震数据体中所对应的地震数据在线号、道号、时间三个方向的梯度，以分别得到所述不同方位的三维地震数据体的线号梯度体、道号梯度体和时间梯度体；(B2)将所述各个方位的三维地震数据体的线号梯度体、道号梯度体和时间梯度体分别进行加权叠加处理，以获得待处理线号梯度体、待处理道号梯度体和待处理时间梯度体，并形成包括待处理线号梯度体、待处理道号梯度体和待处理时间梯度体的待处理梯度体向量；(B3)基于待处理梯度体向量，构建目标工区的待处理梯度体结构张量方阵。可选地，在步骤(B1)中，得到三维地震数据的线号梯度体的步骤可包括：获取目标工区的每个点在三维地震数据体中所对应的地震数据在线号方向的梯度，并将获取的梯度按各个点在三维地震数据体中的位置排列，以构成线号梯度体；得到三维地震数据的道号梯度体的步骤可包括：获取目标工区的每个点在三维地震数据体中所对应的地震数据在道号方向的梯度，并将获取的梯度按各个点在三维地震数据体中的位置排列，以构成道号梯度体；得到三维地震数据的时间梯度体的步骤可包括：获取目标工区的每个点在三维地震数据体中所对应的地震数据在时间方向的梯度，并将获取的梯度按各个点在三维地震数据体中的位置排列，以构成时间梯度体。可选地，在步骤(B2)中，可将待处理线号梯度体、待处理道号梯度体和待处理时间梯度体的顺序纵向排列，以构成待处理梯度体向量。可选地，确定目标工区的每个点在三维地震数据体中所对应的地震数据在线号方向的梯度的步骤可包括：将与线号方向对应的卷积核与目标工区的每个点在三维地震数据体中所对应的地震数据沿线号方向进行卷积，以得到目标工区的每个点在三维地震数据体中所对应的地震数据在线号方向的梯度；确定目标工区的每个点在三维地震数据体中所对应的地震数据在道号方向的梯度的步骤可包括：将与道号方向对应的卷积核与目标工区的每个点在三维地震数据体中所对应的地震数据沿道号方向进行卷积，以得到目标工区的每个点在三维地震数据体中所对应的地震数据在道号方向的梯度；确定目标工区的每个点在三维地震数据体中所对应的地震数据在时间方向的梯度的步骤可包括：将与时间方向对应的卷积核与目标工区的每个点在三维地震数据体中所对应的地震数据沿时间方向进行卷积，以得到目标工区的每个点在三维地震数据体中所对应的地震数据在时间方向的梯度。可选地，与线号方向对应的卷积核、与道号方向对应的卷积核、与时间方向对应的卷积核分别可通过将一维零均值离散高斯核函数的导数在离散变量为相应的取值范围内的各个整数值时的函数值按对应的离散变量的从小到大的顺序排列构成，所述各卷积核中的函数值的计算式为：G′(t)=-12πσi2te[-t2/(2σi2)],]]>其中，G′(t)为一维零均值离散高斯核函数的导数，t为离散变量，t的取值范围为[-Ri，+Ri]，Ri2＝42σi，Ri为核半径，σi为预定尺度因子，i为方向。可选地，在步骤(B3)中可通过下面的等式构建待处理梯度体结构张量方阵：T=ggT=gxLine2gxlineginLinegx本文档来自技高网...

【技术保护点】
一种利用地层叠前纹理属性值来检测缝洞的方法，所述方法包括：(A)分别读取目标工区中的各个点在不同方位的三维地震数据体中所对应的地震数据，其中，所述各个点在每个方位的三维地震数据体中的坐标位置相同；(B)根据读取的地震数据来获取目标工区的待处理梯度体结构张量方阵；(C)从待处理梯度体结构张量方阵的每个元素中抽取与目标工区的一个点的坐标位置相对应的位置处的数据，并将抽取的数据按所对应的元素在待处理梯度体结构张量方阵中的位置排列，以构成所述一个点的结构张量方阵；(D)确定所述一个点的结构张量方阵的特征值；(E)基于所述特征值，获取所述一个点的叠前纹理属性值；(F)利用目标工区所有点的叠前纹理属性值来产生目标工区的缝洞检测结果，其中，通过重复步骤(C)‑(E)来获取目标工区所有点的叠前纹理属性值。

【技术特征摘要】
1.一种利用地层叠前纹理属性值来检测缝洞的方法，所述方法包括：
(A)分别读取目标工区中的各个点在不同方位的三维地震数据体中所
对应的地震数据，其中，所述各个点在每个方位的三维地震数据体中的坐标
位置相同；
(B)根据读取的地震数据来获取目标工区的待处理梯度体结构张量方
阵；
(C)从待处理梯度体结构张量方阵的每个元素中抽取与目标工区的一
个点的坐标位置相对应的位置处的数据，并将抽取的数据按所对应的元素在
待处理梯度体结构张量方阵中的位置排列，以构成所述一个点的结构张量方
阵；
(D)确定所述一个点的结构张量方阵的特征值；
(E)基于所述特征值，获取所述一个点的叠前纹理属性值；
(F)利用目标工区所有点的叠前纹理属性值来产生目标工区的缝洞检测
结果，其中，通过重复步骤(C)-(E)来获取目标工区所有点的叠前纹理属
性值。
2.如权利要求1所述的方法，其中，步骤(B)包括：
(B1)分别确定目标工区中的各个点在不同方位的三维地震数据体中所
对应的地震数据在线号、道号、时间三个方向的梯度，以分别得到所述不同
方位的三维地震数据体的线号梯度体、道号梯度体和时间梯度体；
(B2)将所述各个方位的三维地震数据体的线号梯度体、道号梯度体和
时间梯度体分别进行加权叠加处理，以获得待处理线号梯度体、待处理道号
梯度体和待处理时间梯度体，并形成包括待处理线号梯度体、待处理道号梯
度体和待处理时间梯度体的待处理梯度体向量；
(B3)基于待处理梯度体向量，构建目标工区的待处理梯度体结构张量
方阵。
3.如权利要求2所述的方法，其中，在步骤(B1)中，
得到三维地震数据的线号梯度体的步骤包括：获取目标工区的每个点在
三维地震数据体中所对应的地震数据在线号方向的梯度，并将获取的梯度按
各个点在三维地震数据体中的位置排列，以构成线号梯度体；
得到三维地震数据的道号梯度体的步骤包括：获取目标工区的每个点在
三维地震数据体中所对应的地震数据在道号方向的梯度，并将获取的梯度按
各个点在三维地震数据体中的位置排列，以构成道号梯度体；
得到三维地震数据的时间梯度体的步骤包括：获取目标工区的每个点在
三维地震数据体中所对应的地震数据在时间方向的梯度，并将获取的梯度按
各个点在三维地震数据体中的位置排列，以构成时间梯度体。
4.如权利要求2所述的方法，其中，在步骤(B2)中，将待处理线号
梯度体、待处理道号梯度体和待处理...

【专利技术属性】
技术研发人员：张洞君，杨晓，徐敏，章雄，罗晶，郑虹，陆林超，范晓晓，邹琴，顾雯，兰馨，
申请(专利权)人：中国石油集团川庆钻探工程有限公司地球物理勘探公司，
类型：发明
国别省市：四川;51