可选地，在步骤（A)中，可分别读取目标工区中的各个点在至少三个方位的三维 地震数据体中所对应的地震数据。
[0020] 在根据本发明示例性实施例的利用地层叠前纹理属性来检测缝洞的方法中，能够 有效检测较小尺度的缝洞，并在一定程度上降低缝洞检测的多解性。
【附图说明】
[0021] 通过下面结合附图进行的对实施例的描述，本发明的上述和/或其它目的和优点 将会变得更加清楚，其中：
[0022] 图1是示出根据本发明示例性实施例的获取地层纹理属性的方法的流程图；
[0023] 图2示出根据本发明示例性实施例的方位为p度_q度的方位角区间的三维地震 数据体的空间模型；
[0024] 图3示出根据本发明示例性实施例的根据读取的地震数据来获取目标工区的待 处理梯度体结构张量方阵的步骤的流程图；
[0025] 图4示出根据本发明示例性实施例的根据读取的地震数据获取包括线号梯度体、 道号梯度体和时间梯度体的方法的流程图；
[0026] 图5A示出现有技术的缝洞检测效果图；
[0027] 图5B示出采用本发明示例性实施例的利用地层叠前纹理属性值来检测缝洞的方 法的缝洞检测效果图。
【具体实施方式】
[0028] 现将详细描述本发明的示例性实施例，所述实施例的示例在附图中示出，其中，相 同的标号指示相同的部分。以下将通过参照附图来说明所述实施例，以便解释本发明。
[0029] 图1示出根据本发明示例性实施例的利用地层叠前纹理属性值来检测缝洞的方 法的流程图。这里，作为示例，所述方法可由用于利用地层叠前纹理属性值来检测缝洞的设 备来实现，也可完全通过计算机程序来实现。
[0030] 如图1所示，在步骤S100,分别读取目标工区中的各个点在不同方位的三维地震 数据体中所对应的地震数据，其中，所述各个点在每个方位的三维地震数据体中的坐标位 置相同。作为示例，所述目标工区可以是一个地层、层位或一个由多个地层组成的工程工 区。所述方位可指示方位角或方位角的区间。这里，所述方位角是指激发点和检波点的连 线与正北方向的夹角。
[0031] 这里，作为示例，为了保证缝洞检测的准确性，可分别读取目标工区中的各个点在 多个不同方位的三维地震数据体中所对应的地震数据，例如，可分别读取目标工区中的各 个点在至少三个方位的三维地震数据体中所对应的地震数据。具体说来，可通过对采集的 原始地震数据进行地震资料预处理（例如，球面发散校正、对地震吸收的Q补偿、振幅处理、 子波反褶积、地表一致性静校正、速度分析、动校正以及剩余静校正等）来获取同一目标工 区不同方位的三维地震数据体，并进而获取目标工区中的各个点在不同方位的三维地震数 据中所对应的地震数据。
[0032] 此外，作为示例，图2示出方位为p度_q度的方位角区间的三维地震数据体的空 间模型。在三维地震勘探中，得到的三维地震数据体的空间模型如图2所示，图2中的黑点 代表地震数据，可将地震数据表示为D p q(x，y，z)。这里的X表示线号，y表示道号，z表示 时间。地震数据Dpq可理解为目标工区中由线号（inLine)、道号（xLine)、时间（time)所 确定的点（即，采样点）在方位为P度_q度的方位角区间的三维地震数据体中所对应的地 震数据。在这里，地震数据可采用部分方位叠加的地震数据，具体地可以为振幅等。可以理 解，X、y、z为取整数的离散变量。
[0033] 再次参照图1，在步骤S200,根据在步骤SlOO读取的地震数据来获取目标工区的 待处理梯度体结构张量方阵。下面，将结合图3和图4来详细描述如何获取目标工区的待 处理梯度体结构张量方阵的方法。
[0034] 图3示出根据本发明示例性实施例的根据读取的地震数据来获取目标工区的待 处理梯度体结构张量方阵的步骤的流程图。
[0035] 如图3所示，在步骤S201，分别确定目标工区中的各个点在不同方位的三维地震 数据体中所对应的地震数据在线号、道号、时间三个方向的梯度，以分别得到所述不同方位 的三维地震数据体的线号梯度体、道号梯度体和时间梯度体。下面，将结合图4来详细描述 如何得到所述不同方位的三维地震数据体的线号梯度体、道号梯度体和时间梯度体。
[0036] 如图4所示，以一个方位的三维地震数据体为例，在步骤S301，获取目标工区的每 个点在三维地震数据体中所对应的地震数据在线号方向的梯度，并将获取的梯度按各个点 在三维地震数据体中的位置排列，以构成线号梯度体。这里，可通过各种梯度计算方法来 获取各个点在三维地震数据体中所对应的地震数据在线号方向的梯度，从而得到线号梯度 体。
[0037] 优选地，可利用本发明提出的方式来获取所述各个点的地震数据在线号方向的梯 度。具体说来，将与线号方向对应的卷积核与三维地震数据体中每个点所对应的地震数据 沿线号方向进行卷积，以得到各个点在三维地震数据体中所对应的地震数据在线号方向的 梯度，计算式如下式（1): CN 105158795 A ~P 5/8 页
[0038]
(I)
[0039] 这里，haiXUne(x，y，z)为三维地震数据体α中的坐标位置为（X，y，z)的点的地震 数据在线号方向的梯度，为与线号方向对应的卷积核，*为卷积符号，D a (X，y，z)为三 维地震数据体a中由x、y、z的值确定的点的地震数据。
[0040] 这里，与线号方向对应的卷积核fxU JI过将一维零均值离散高斯核函数的导数在 离散变量为与线号方向对应的取值范围内的各个整数值时的函数值按对应的离散变量的 从小到大的顺序排列构成。具体地，对一维零均值离散高斯核函数求导数，得到等式（2):
(2)
[0042] 这里，G'（t)为一维零均值离散高斯核函数的导数，t为离散变量，t的取值范围为
为与线号方向对应的核半径，σ 为与线号方 向对应的预定尺度因子。
[0043] 接下来，将与线号方向对应的一维零均值离散高斯核函数的导数G'（t)在离散变 量t为与线号方向对应的取值范围[_RxUra，+RxunJ内的各个整数值时的函数值按对应的离 散变量t的从小到大的顺序排列得到卷积核。例如，假设与线号方向对应的预定尺度因子 σ xUne为1，则相应的核半径R 等于4。那么，离散变量t的取值范围为[-4, +4]，t的从 小到大的取值为{-4, -3, -2, _1，0, +1，+2, +3, +4}。将离散变量t的各个取值代入到等式 (2)中则得到相应的一维零均值离散高斯核函数的导数的函数值分别为{0. 0002,0. 0053， 0· 043,0· 0965,0, -0· 0965, -0· 043, -0· 0053, -0· 0002}。最后，按对应的离散变量 t 的从 小到大的顺序将所述函数值排列得到与线号方向对应的卷积核(0. 0002,0. 0053， 0· 043,0· 0965,0, -0· 0965, -0· 043, -0· 0053, -0· 0002) 〇
[0044] 在步骤S302,获取目标工区的每个点在三维地震数据体中所对应的地震数据在道 号方向的梯度，并将获取的梯度按各个点在三维地震数据体中的位