地震解释的现状及发展趋势

地震波地质信息综合解释摘要：地震解释质量决定了一个区块勘探开发的方向和进程，地震解释的发展对解释人员提出了更高的要求，即要求解释人员通晓地质知识，同时具有物探知识。本文主要从现今已经在应用的解释技术和方法以及近年来涌现出来的一些新思路、新方法展开论述。分别包括三维可视化技术、构造解释、构造解释和利用振幅属性预测含烃概率、利用波峰瞬时频率计算薄层厚度、多子波地震道分解和重构等。关键字：地震解释、构造解释、振幅属性、波峰瞬时频率引言：地震资料解释是勘探和开发地震的最后环节，其功能是将地震信息翻译成地质语言或符号;其目的是直接服务于勘探和开发。因此解释质量决定了一个区块勘探开发的方向和进程。地震勘探开发技术发展的目标都是为了提供更好的易于解释的具更高可信度的地震资料。地震解释现在更多地强调综合性和在地质规律控制下的地震解释。这对解释人员提出了更高的要求，即要求解释人员通晓地质知识，同时具有物探知识。地震解释从来就不是从事物探方法研究人员单纯可以从事的工作。地震解释已经开始从注重地震解释方法向注重多学科综合性的转变，现在更为明显!地震解释的另一个明显的趋势是强调在地质规律认识下的地震解释，即地震和地质的紧密结合。一、地震综合解释的现今技术及方法在地震综合解释方面，主要是以地震反演技术、多种属性分析技术及三维解释为主体的地震综合储层预测技术，通过与层序地层学、测井和地质等其他测量解释成果的结合给出地震资料综合解释的应用实例。例如AmoutColpaert应用神经网络将地震解释数据和井中岩石物理特性分析联合实现多属性分析，从而进行岩相预测。靶区的目标地层是岩溶发育的斜坡形向陆架坡过渡的碳酸盐岩地层，探区内井资料很少或几乎没有，作者综合应用了基于井资料的层序地层分析、岩石物理分析和多属性地震分析，对无井控制区的岩相进行了预测。其基本流程见图1。而涉及储层预测的主要是综合应用地震、地质、钻并和测井等各项资料对地下储层的分布、厚度及岩性和物理性质变化进行追踪和预测的实例。GuoYabin给出了利用对苏里格气田的3C一2D的P波及转换波波形特征分析、AVO分析、地震反演、振幅分析、vp/vs分析结果进行了气藏预侧的实例。Kontorovichdarya给出了综合地震、测井和钻井数据来进行储层预侧和复杂圈闭描述的实例。Caolaisheng通过联片3D地震数据和关键井来进行地震相分析和盆地填充分析，具体方法包括:合成记录与水平切片进行相关、同一位置上2D剖面和3D重叠剖面的比较、两步f-k滤波、地震相干体扫描以及用于确定“采集足迹”的时间切片，最后用地震属性分析和波形分类来确定徽相。BaiYuhua给出的是基于高分辨率地震数据，联合测井、地质及钻井数据对河流相砂岩储层进行了预测的实例，具体为:用3D地震相干技术进行沉积亚相的预侧，用层切片技术预侧砂体展布，用水平切片技术确定砂体边界，用3D可视化技术描述砂体轮廓，用谱分解技术解析河道结构。三维可视化技术三维可视化技术是全三维解释技术的核心技术，它借助于无限实时图形工作站的先进显示技术和先进的可视化软件，对地震数据体进行全方位的透视，用图形直观地显示地质、侧井、地震等数据，因而可以直观、高效和便捷地震示三维地质信息，并在三维立体空间显示。它不仅使地球科学家们能更深刻地理解各种地质现象的发生、发展及影响，而且使他们的想象力更加丰富多彩，以直观的方式使人们加深对数据的理解，从中得到新的见解。比较新颖的LaurentCastanie提出的先进而稳健的体可视化算法。该算法基于可编程的图形处理器(GPU)和转换函数预积分，即假定视线)(X是由到视点的距离参数化的，)()(XXc和分别定义了点X处的颜色和消光系数，则可定义体绘制积分为ddXXcID})]([exp{)]([00。其中，颜色和消光系数是由转换函数确定，该转换函数将标量属性空间映射到了RGBA(红、绿、蓝、透明度)空间。然后对所有可能的切片的前景和背景，在转换函数空间对上述积分进行计算，计算结果分别存贮为2D的纹理bfss和sf和几，在体绘制时再将这些纹理提取出来，提取是由运行在GPU上的“片断程序”完成的(图2)。构造解释PierreJaequemin提出了基于两次霍夫变换的自动断层提取技术:断层面近似地由(x，z)面上的一系列直线构成，在第1个参数空间用霍夫变换把每条直线变为点，则在参数空间，断层变成了由这些点构成的直线，然后再用霍夫变换把该直线变为第2个参数空间中的点。用反霍夫变换就可以把每个断层变为一系列点。IsraelCohen提出了在3D地震数据中检测和提取断层面的方法。具体步骤如下:(1)解释人员在一个分析点上倾斜并旋转一个3D数据体来选出一个3D分析数据体:tNkxyttxytxyttxytxytdNddEdd11}{ˆ式中，xytd和xytdˆ分别是某一道在点（x，y）处的第t个原始的样点和修改后的样点;Nt是每一道的总采样数。这样就选择了一个相对小的3D分析数据体。(2)把该3D分析数据体移遍整个3D数据体，并在每个点上产生一个归一化差分谱(NDE)。NDE值代表断层面的相似性即在分析点上与分析数据体具有相似的倾角和方位角。这个分析数据体是由主轴长度1L、次轴长度t2N12、L个采样点的持续时间、方位角和倾角所决定(图3)。这样，与关于分析点),,(tyx处倾斜并旋转的数据体相比，在2个子空间上的采样点就变成了2个列矢量),(),(2,22,1VV和。在分析点),,(tyx处的NDE由公式ppkVVVVN),(),(),(),(),(,2,1,2,1决定。如果没有断层，则两个列矢量),(),(,2,1VV和相等，从而该点的NDE为零。(3)用离散的墨西哥帽函数对NDE进行对比增强滤波，再对经过增强滤波的NDE进行方向滤波来提取断层面片段。将经过方向滤波后的NDE系数设为阀值并再次进行滤波得到方向局部断层提取(LFT，)数据体，LFT为解释人员提供了一个非常清晰断层面的可视化标识，最后由最大方向LFT给出LFT分布。油气检测油气检测就是综合利用烃类存在的多种地震特性参数(速度、频率、振幅和相位等)来确定油气富集带。这类技术有许多种，比较常用的有亮点技术和AVO技术等。TracyJ.Stark提出了借助可视化软件，用交汇图、异常体及颜色叠加来确定AVO异常的实例。对于一个3D数据体，在1000ms的窗口内，分别实现近、中、远道的叠加，并用不同颇色标示(图4，图5)。然后用测线上不同距离益加后的数据的时间切片而不是时间窗来产生AVO交汇图，从而提高了纵向分辨率(图6，图7)，最后借助可视化软件将异常体显示出来(图8).SatinderChoPra将纹理属性分析的统计方法用于3D地震数据体来刻画储层的油气分布。从3D地震数据体中选一个小数据体(类似于纹理分析中的纹理元)，由这个小数据体产生灰度矩阵，用统计方法研究其空间上的分布特征—灰度级同现矩阵，及其各种其他纹理特性:能量、摘、对比度、均匀度，从而提取出局部数据的纹理特性;将上述做法遍历整数据体，则可以提取出整个3D数据体的纹理特性。最后锐化整个数据体纹理特性来凸现储层的油气分布。FredAminzadeh提出了基于神经网络的地震目标探测技术。神经网络的输人是不同偏移距叠加道集的AVO属性等地震属性，用已知地质信息和井的信息来训练网络，这样在训练网络时就可以加人专家知识从而提高预测成功率。TracyJ.Stark等由高分辨率年代数据体在3D地震数据体中直接产生地层年代，这样一与产生它的地震和年代数姻体具有相同的侧向尺度.纵轴代表相对地质年代，绿色区是局部的不确定区，而年代体中的“无致据区”则体现在品红色区方面可以对整个3D数据体而不是只对单个剖面进行标定，另一方面也避免了在地震剖面上手工标出地层年代的麻烦.他首先是产生一个地质年代的数据体，然后把某一旅行时的振幅值加到与具有相同旅行时的地质年代上，从而产生地震轮幅体(见图9)。JianleiLiu等用Morlet小波对地震数据进行匹配追踪分解从而分析其瞬时频率特性，其流程见图10.假定地震信号可以分解为Morlet小波的形式Ntfmatddmi),,()(式中，d(t)是地震道;im为第i个Morlet小波;a为振幅;mf为平均频率;dt为膨胀时间;为相角;N为噪声。可以用复数属性分析获得振幅和膨胀时间。膨胀时间是包络的峰值时间，相角为瞬时相位角，平均频率与瞬时频率相等。接下来计算复数地震道的Morlet小波的复数属性，并提取与Morlet小波相匹配的地震道的实部。上述过程反复进行，直到能量小于门槛值，此时就能得到与Morlet小波相匹配的模型道。最后用模型道就可以计算时频分布。此外，HongliuZeng用90度相位子波进行了薄层的地震解释;RockyRoden在分析了不同含油气砂岩振幅响应特征的基础上讨论了地震波振幅在勘探风险评估中的作用;SteveCarlson用模式识别工具来刻画地震信号及含油砂岩的AVA响应;KathleenBaker分析了露头在测井、地震及直接油气检测方面的相似性，从而研究了岩性对地震信号的影响。构造建模构造建模是地震解释工作的一个重要阶段和成果。目前越来越多的构造建模都是综合应用地表地质、遥感、钻井、重、磁、电和地震等资料，将地表构造、浅层构造与深层构造有机结合起来从而实现对构造形态、特征和要素的精细描述。E，Mon*n提出了一个新的一体化的体建模流程，该模型可以在多个尺度上同时满足几何、地质和地球物理需要，用该流程可以建立起蓦本模型的分级结构，并能同时保持和处理的分辨率，但是其中心思想仍然是基于体的解释。利用重、磁、地震和侧井数据来反演基地构造。首先是把这些测量信息融合到GIS中，然后用地球物理技术建立一个对各种侧量信息都适合的基本构造，然后用该基本构造来对区域内的局部地质体进行取舍而产生较为精细的构造，再将该精细的构造应用到所侧量的信息，从而实现对该构造的细化。二、地震综合解释的进展及发展趋势利用振幅属性预测含烃概率利用振幅和井资料预测含烃概率最简单的方法是选取门槛值，这在实际工作中经常使提是工区的振幅变化规律强。Denham和Johnson等提出利用高质量的侧并资料获得砂岩振幅随深度变化的信息，应用Greenberg一Castagi扭横波估算技术(Greenberg和Castagna，1992)和Gassmann方程(Gassmann，1951)将含烃砂岩的振幅同含水砂岩的振幅区分开，并对一定范围内的井取平均值，进而预测未钻井区块砂岩的含烃性。该技术应用到的数据库中包括墨西哥清5500多口井的岩石物性信息。按200ft的间隔对砂岩和泥岩物性取平均值。泥值含量和温度用侧井获得数值进行插值。含砂量大于55%的岩石认为是砂岩，含砂量小于55%的岩石认为是泥岩。运用统计的方法三维空间(地理位置和海平面深度)指定点的含烃概率可以被计算出来。一般而言，目标点是没有含烃量信息的，但可以通过在数据库中根据指定的搜索半径搜索出井资料及目标深度处的含烃量信息，当目标点满足以下条件并根据以下步骤，每个点的含烃概率都可以被刻画出来。(l)如果该点的资料信息满足不了法向人射时页岩一水层的反射系数计算、密度求取和P波速度求取，则该点不能参与计算;(2)计算水层的法向反射系数;(3)应用Greenberg和Castagna(1992)提出的方法计算砂岩和泥岩的横波速度;(4)应用Gassmann方程计算干层体积模t，基质体积模童，干层硬度，岩石基质密度和孔隙度;如果某点Gassmann反演结果不符合物理规律，则该点在进一步计算时应抛弃;(5)Gasslnann被用来计算当岩石分别为水渗层、油层和气层时，砂岩速度和密度的变化规律，(6)计算法向入射时，油层和气层的反射系数。利用波峰瞬时频率计算薄层厚度JianleiLiu等发表了文章“利用波峰瞬时频率计算薄层厚度”，认为波峰瞬时频率反比于地层厚度