物探新方法新技术之五相干╱方差体技术(Coherency╱Variance Cube).doc

物探新方法新技术 PAGE PAGE 84 岩性地球物理 5 相干/方差体技术 5.1 相干/方差体技术的基本原理 相干/方差体技术通过量化处理地震数据体的相干属性，生成新的相干/方差体数据体，突出和强调地震数据的不相关性。 相干/方差体技术利用相邻道地震信号之间的相似性来描述地层、岩性等的横向非均匀性，特别是在识别断层以及了解与储集层特征密切相关的砂体展布等方面非常有效。应用三维相干/方差数据体时间切片进行构造解释和岩性解释，可以帮助解释人员迅速认识整个工区断层等构造及岩性的整体空间展布特征，从而达到加快解释速度及提高解释精度、缩短勘探周期的目的。 作为三维地震解释和岩性研究重要的技术手段，相干/方差体技术应用和发展非常迅速，已经成为一种常规手段用于地震构造和岩性解释当中。 我国东部的主要可采煤层大部分赋存在石炭二叠纪地层中，由于不同年代的构造运动，煤系地层中的断层极为发育。同时，我国华北各煤田岩溶陷落柱也比较发育，主要是指石炭二叠系煤系地层下部的奥陶系灰岩中古溶洞(裂隙)的塌陷形成的柱体。断层与陷落柱均是对煤矿开采十分不利的构造。因此，利用三维相干/方差数据体进行构造、岩性解释，主要是针对我国目前煤矿生产中遇到的实际问题。另一方面，这项技术也能够提高煤田地震资料的解释水平和速度。 5.2 相干数据体的算法 相干数据体技术很早就在二维地震勘探中被提出来了，但其应用于三维地震解释中是从Bahorich和Farmer(1995)提出的基于互相关的相干性C1算法开始的。近几年来，国内外文献中有关三维地震相干数据体生成方法的文献很多，形成了很多算法。但是在实际地震勘探资料解释中具有较好的应用效果的主要是基于水平时间切片，并根据资料的信噪比、算法的稳定性、利用特征值分析方法进行相干分析和处理，发展了互相关算法、相似性算法、本征值算法三种应用比较广泛的算法。 三维地震数据体反映地下一个规则网格的反射情况，通过对主测线和联络测线方向计算某一时间域内波的相似性，可获得三维相关图。当地下存在断层和地层不连续性变化时，在局部一些地震道上会表现出与相邻地震道不同的反射特征，因而导致道与道之间相关性方面的极不连续，即断层出现的地震错动，会在相应道的相关曲线中出现极高的不相关性，见图5—1。利用这一原理，通过对三维数据体的不连续性进行分析，便可识别构造和断层的分布，使解释人员在解释之前就能获得研究区概略的构造几何形态及断层分布情况，充分利用三维地震数据体原已存在的空间分布信息，能够减少复杂情况下人为因素造成的误差及多解性。 图5—1 由断层引起的波形变化 5.2.1 协方差矩阵 对于三维地震数据体(坐标为)任一格点上的数据道，在其周围分别沿着与其相邻的inline线、crossline线方向上，定义一个总共包含道的平面数据子集区域，给定某一切片时间t(其中,为时间采样间隔)，在以t为中心，w为半时窗、2w+1时窗长度上，计算数据子集中任意两道的互相关。 Marfurt引入了协方差矩阵来协助完成相干算法的实现，对于每个延迟对形成了如下一个高阶的协方差矩阵。假设以为中心的一对视倾角()的个采样点，这个采样点对应着一个的协方差矩阵 (5-1) 式中，为矩阵，即在进行和相干计算中所要用到的协方差矩阵，其中表示地震道沿着视倾角在处的内插值。 图5—2是两道互相关系数的计算方法，其中设置搜索时窗的目的是按层的走势来寻找最大相关值。具体方法如下： (1) 中心道数据不动，以样点间隔为单位，从上向下移动相邻道，每移动一次，都要计算时窗内波形的相关系数，直至获得如图5—2中右边所示的相关函数，其长度为搜索时窗长度，值在-1到1之间。 (2) 通过求取其中的最大值，即两道最相关处的值来表明在该处两相邻道具有最大的相似度，见图5—3。 图5—2 两道互相关系数算法 图5—3 地震相关值 5.2.2 相干算法 该方法是基于传统的归一化互相关原理来实现相干体运算的。相干算法计算每道的横测线和相邻纵测线的互相关，混合两个结果并用能量进行标准化，其实质是两道相关算法的扩展。 设三维地震数据体中inline线上和两个相邻的地震数据道之间在时间、延迟上的互相关函数为，则有 (5-2) 同样，对于三维地震数据体中crossline线上和两个相邻的地震数据道之间在时间、延迟上的互相关函数为，则有 (5-3) 因此，对于利用三维数据体上相邻三个数据道，进行相干体计算的相干算法而言，沿着inline、crossline线方向的时间视倾角，所计算的相干度为 (5-4) 上式中的和表示延迟和使得最大化，对于高质量的资料，其值近似等于各方向的时间视倾角(对于有相干噪声的资料，仅用两道数据确定视