浅谈地震勘探处理方法

1、. 浅谈地震勘探处理方法论文提要地震勘探技术在油气田勘探开发中起着重要的作用。地震勘探包括三大阶段，野在采集，数据处理和室解释。其中地震数据处理的目的是对地震采集数据做各种处理提高反射波数据的信噪比，分辨率和保真度以便于解释。地震数据处理主要包括地震反褶积，叠加和偏移成像三大技术。地震数据处理出现于20时纪20年代初期，随后的40年间是对光点记录和模拟记录进展处理，处理技术开展较慢，进入20世纪60年代以后，计算机的出现把地震勘探处理技术带入了数字时代，数字技术为数据处理的开展提供了广阔的前景。下面简单介绍地震数据处理的流程以及地震数据处理的方法。正文一、地震数据处理流程一地震数据处理的三个阶

2、段1.预处理，预处理就是把野在数据格式转换成适合计算机处理的格式，并对数据做相应编辑和校正。预处理包括数据解编，格式转换，编辑，几何扩散校正，建立野在观测系统，野在静校正2.常规处理，是对预处理后地震数据做必要的根本处理预算。常规处理包括反褶积，道均衡，抽取共中心点道集，速度分析，剩余静校正，切除，叠加，偏移。3.特殊处理，针对不同目的采用不同的特殊的处理手段，包括t-p变换，小波变换，三维叠前深度偏移，子波处理，属性分析，反演。二、数字滤波一数字滤波的有关概念从广义上讲，任何一种对输入信号的改造作用都可看成滤波，实现这种滤波的系统称为滤波器，滤波器分为模拟滤波器和数字滤波器。1.模拟滤波器，

3、也称电滤波器，它由电阻、电感和点容等元器件组成，它组成的是一个低通滤波器LCF如图1。RC图1由于模拟滤波器运算速度快，因此*些具有单一滤波功能的构件可由它来完成，但模拟滤波器一旦固定，不易修改，适应面较窄，本钱也较高，所以模拟滤波器进一步开展成了数字滤波器。2.数字滤波器，数字滤波器主要目的是压制噪声，信号要进展数字滤波，首先要进展采样。抽样过程要满足抽样定理，不然会使频谱混叠，产生假频，抽样定理可由以下两个公式描述。1频率域 Ws=2WN=2Wma*式中Ws称为采样频率，WN称为折叠频率，也称为Niquist频率Wma*称为最高频率。2时间域t=1/2fma*式中t为采样间隔，fma*为信

4、号的最高频率。二一维滤波利用有效波和干扰波的频率差异，采用一维滤波，主要压制面波。1.理想低通滤波器有一些地区存在较强高频干扰，因此需要通过低通滤波把高频干扰除去，这种滤波器的形状像门一样，故又称门式滤波器。如图2-CC|H()|0图22.带通滤波器有效波在有限频带围出现高、低干扰波，仪器车上滤去高频低频，频率扫描，检查剖面分辨率和信噪比。如图3|H()|-2-1210图33.带陷滤波器-CC|H()|10又叫带阻滤波器，主要消除*些特殊的干扰，如工区经过高压线即50Hz工业电干扰。如图4图44.高通滤波器地震资料中，有些干扰频率很低，如面波，可用高通滤波器。如图51|H()|-CC0图5三二

5、维滤波1.二维滤波要满足二维抽样定理即时间采样间隔和空间采样间隔应同时满足。V1*=fkV2*=fkf2f10k2.二维频波图，视速度一样的信号成分在频波图上位于过原点的直线上，而且斜率越大，视速度越大。如图6I图 63.扇形滤波器，在实际地震资料中，由于地层的视速度一般都很大，因此常用的二维频率域是扇形滤波器。如图7NOK/k=-v/k=v-N图74.视速度滤波的处理流程如图8所示。输入地震数据二维傅氏变换扇形滤波二维傅氏反变换图8三、速度分析、动静校正和叠加一速度分析速度分析是为叠加提供最正确叠加速度。通常有两种地震波速度的测量方法，一是声波测井的直接测量法；二是地震勘探数据的间接测量法。

6、根据这两方面的信息，可以得到多种不同物理意义和用途的地震速度，如层速度，平均速度，均方根速度，瞬时速度，相速度，群速度，动校正NMO速度，叠加速度和偏移速度等。1.层速度层速度为两个反射界面之间的平均速度，一定岩石组分岩层的层速度受下面几种因素影响：孔隙形状，孔隙压力，孔隙流体饱合度，围压，温度。2.速度扫描原理对一个CMP道集，用一系列常速度重复做动校正，每校正一次得一图像，当速度小时，校正的结果为过校正，速度大时，欠校正，只有当速度适宜时，校正后道集构成一条水平直线，此时的速度为所求的叠加速度，此方法精度高，工作量大。3.常速度叠加法对多个相邻CMP道集，用一系列速度进展校正，叠加，叠后同

7、相轴振幅最强且是连续性好的所对应的速度即所求叠加速度，此法适合解决复杂构造。4.速度谱对共反射点时距曲线，固定to求取ti，取各道记得振幅进展相加得A(V)，当V=Va时计算时间与实际记得同相轴的相位时间一致，由于同相轴的相位时间不一致，叠加振幅小，一系列叠加速度中，振幅最大的速度即所求叠加速度。图9a表示一个层状地层模型的CMP道集，根据叠加振幅，从速度谱上选取下面几个作为叠加速度函数如图9b：2700,2800,3000m/s分别对应浅、中、深三种同相轴，速度谱不仅能提供叠加速度，使我们由此区分反射波屡次波ab图95.影响速度估算的因素以下因素会限制地震资料速度估算的精度和分辨率：排列长度

8、，叠加次数，信噪比，切除，时窗宽度，速度采样，相干属性量的选择，双曲线正常时差的偏离度，数据的频谱宽度二动校正动校正起消除正常时差的过程。1.动校正量随炮检距、深度和速度的关系。NMO随炮检距，随深度递减，同时，速度越大，NMO值越小。2.动校正速度的选择对校正效果的影响，如图10，采用低速度，得到过校正道集，采用高速度，得到欠校正道集。图103.动校正拉伸动校正结果出现了频率畸变，尤其是对浅层的大炮检距，就是NMO拉伸。通常要考虑信噪比和切除两方面，采取折衷，特别是当信噪比很高时，宁愿多切除一些以保存信号频带宽度。反之，如果信噪比很低，这时应多保存一些拉伸，以便叠加之后得到一切可能得到的同相

9、轴，如图11。图114.各种地层模型的NMO速度1单一水平层，反射界面以上的速度2水平层状地层，由小排列给出的均方根速度函数3单一倾斜层，界面以上介质速度除以倾角的余弦4多层任意倾斜层，由小排列，小倾角定义的均方根速度函数三静校正由于激发点和接收点条件变化引起的时差称静校正量，消除静校正量的过程称作静校正。静校正量随地表高程变化和浅层低、降速带的存在而发生变化，其中激发点的静校正量为正，反之，接收点的静校正量为负。静校正包括基准面校正，就是把所有炮点和接收点都校正到同一水平面上，这个水平面称基准面。基准面校正就是把所有炮点和接收点都校正到同一水平面上的过程。由近地表不规则所造成的畸变绝大局部通

10、过野在静校正予以消除，下面为剩余静校正的处理流程图，如图12.经野外静校正的CMP道集速度分析利用初始拾取的速度进展NMO校正计算剩余静校正量应用剩余静校正速度分析利用最终拾取的速度进展NMO校正叠加图12四叠加1.一般水平叠加，参与叠加道的质量根本一致2.自适应水平叠加，参与叠加道质量差异较大。四、反褶积反褶积是地震资料最常用和最重要的处理方法之一。它可用于叠前，也可用于叠后。反褶积的主要作用是压缩地震子波，提高地震资料的分辨率，从而提高资料解释的精度，为油田精细勘探和开发效劳。另外，反褶积还可以消除短周期鸣震和其它屡次干扰，突出有效波，提高地震资料的信噪比。一反褶积的根本概念假设地震波以脉

11、冲波形式激发，经过地层时无吸收，透射和屡次反射等因素的影响，而且整个传播过程不存在随机干扰，这样可以得到理想的输出，这时得到的输出实际上就是反射系数序列，同实际地震记得相比，它有很高的分辨率。1为什么要做反褶积，实际地震记得由于受复杂子波的作用和干扰的影响，分辨能力较低，地质界面上各反射波互相叠加，彼此干预，成为一复杂的形式。由一函数与地震记得褶积得到反射系数，这一过程被称为反褶积，有时也称为反滤波，即反射系数与地震子波褶积，地震记得与反子波褶积得到反射系数。二最正确维纳滤波维纳滤波即最小平方滤波，是由维纳1942年提出的，这种方法以一种最正确准则来设计滤波器，使滤波器的实际输出与期望输出的差

12、的平方和最小。因为维纳滤波器是一种最正确滤波器，维纳滤波又经常被称为最正确维纳滤波器。1.下面介绍维纳滤波的根本原理，如图13。期望输出输入互相关自相关维纳滤波器褶积实际输出图132.维纳滤波与各种反褶积间的关系假设期望输出是零延迟尖脉冲，则有脉冲反褶积；假设期望输出是时间提前了的输入序列，则有预测反褶积，特别是当预测步长为1时，这时预测反褶积成了脉冲反褶积；假设期望输出是任一延迟尖脉冲或波形，则有子波整形反褶积。下列图则是维纳滤波与反褶积方法关系图，如图14。最正确维纳滤波期望输出向前移时的输入预测距为零延迟尖脉冲任意形式预测距为1子波整形反褶积预测反褶积脉冲反褶积图14三其他反褶积方法1.

13、最小平方反褶积，输入是地震子波，输出是脉冲，用最小平方反褶积。2.预测反褶积，用预测滤波原理来实现反褶积问题的方法就是预测反褶积。下面介绍预测反褶积的根本流程。输入地震记得，计算n+a个延迟的自相关n为预测滤波器长，a为预测步长。用方程如图15计算滤波因子，滤波因子与输入相褶积，结果按预测步长延迟，由输入减去刚刚得到的结果得到输出。用预测反滤波消除屡次波是典型，还可以消除海上鸣震干扰。图 153.子波整形反褶积1子波振幅谱一样，最小相位子波对期望输出为零延迟的反褶积，误差最小；2在子波为混合相位和最大相位时，期望输出的相位应与子波的相位匹配，只有这样才能得到适宜的反褶积结果。五、偏移成像地震偏

14、移技术是现代地震勘探数据处理的三大根本技术之一，它是在过去的古典技术上开展起来的，其它两大技术都是从其它相关学科引进到地震中来的。所以，偏移技术具有地震勘探本身的特征。但是，地震偏移方法本身由于使用计算机而引起了许多革命性的变化。这就使得它从研究简单的探测目标的几何图形进而开展成研究反射界面空间的波场特征，振幅变化和反射率等。一地震偏移效果叠前见图16、叠后见图17。图16 图17二偏移方法分类1.叠加：适用于水平层状介质。2.法向射线深度转换：严格用于没有构造倾角且速度只随深度变化的情况。3.时间偏移：适用于叠加剖面上有绕射波或构造倾角以及速度有垂直变化的情况；速度的横向变化不大时也能用。4

15、.深度偏移：用于叠加剖面上有构造倾角和强横向变速的情况。5.叠前局部偏移：叠后偏移适用于叠加剖面与零炮检距剖面等价的情况，但不适合具有不同叠加速度的地层倾角不一致或强横向变速的地区，叠前局部偏移能够为叠后偏移提供更好的叠加剖面，但叠前偏移只解决具有不同叠加速度的地层倾角不一致的问题。6.叠前全时间偏移：输出偏移剖面，不产生未经偏移的中间叠加剖面，但无论如何这是解决倾角不一致地层问题的最准确方法。叠前局部偏移是这种处理方法的一种简化。7.叠前深度偏移：用于严重横向变速的情况，这时无法做适宜的叠加处理。8.三维叠后时间偏移：用于叠加剖面上出现来自射线平面以外的倾斜同相轴即垂直测线方向的情况，这是叠后最常用的一种三维偏移方法。9.三维叠后深度偏移：用来解释与三维地下复杂构造有关的强横向变速问题。10.三维叠前时间偏移：