地震勘探数据处理.ppt

反射地震资料的数字处理 反射资料的数字处理流程 水平叠加处理 反射资料的数字滤波和反滤波 速度分析 偏移归位处理 （一）反射资料的数字处理流程 （二）水平叠加处理  （三） 修饰性处理  （四）地震反射资料的数字滤波和反滤波 （五） 速度分析 （六）偏移归位处理 ② 子波未知情况下的基本方程 地震子波往往事先不知道。为了在子波未知情况下求出反滤波因子 a(t) ，必须对子波及反射系数序列加上一定的限制，即假设条件： A、 反射系数序列R(t)是随机的白噪序列，即 B、 假设地震子波是最小相位的。 由A可证明子波自相关rbb=记录自相关rgg。 由B可知m0=0和令a’(t)= a (t)/b(0)，有 4) 用最小平方反滤波消除大地滤波作用 ① 理想的地震记录 震源：δ(t)；反射系数：R1, R2,…, RN-1, RN。 ② 实际的反射地震记录 ③ 反滤波 从实际所得的反射地震记录中去掉地层滤波器的作用，使它变为理想的地震记录，这个过程称为反滤波。其目的是压缩地震子波的延续时间，提高地震记录的分辨率。 具体做法是，根据事先已求出的地震子波，求出其自相关函数，代入方程求解即可得到反滤波因子a(t)，用它与地震记录褶积可得到 即反滤波结果接近于理想地震记录（反射系数序列），因而大大提高了纵向分辨率。  实际中，因地震子波未知，用记录自相关rgg代替子波自相关rbb求解反滤波因子a(t)。 一般我们把压缩信号时间长度的反滤波，叫做脉冲反褶积。 反褶积 ⑵ 预测反滤波 1) 预测滤波的原理 设计一个预测滤波因子c(t)，对输入信号g(t)已知的过去值g(t-m), g(t-m+1),…, g(t-2), g(t-1)和现在值g(t)进行滤波处理，获取未来某时刻(t+a)时预测值（a为预测距离） 同反滤波方法类似有 rgg为记录自相关。 求解方程即可得到预测滤波因子c(t)，用它与地震记录褶积可得到预测值 2) 预测反滤波的基本思想 将上述预测滤波理论用于解决反滤波问题称为预测反滤波。 只要确定了预测步长a之后，由前面方程组就可得到预测滤波因子c(t)，然后就能直接写出预测反滤波因子： 物理量从可否预测角度而言可以分为二大类：可预测量和不可预测量。一个滤波器的输出也由两部分组成：脉冲响应和输入内容，前者为可预测部分，后者为不可预测部分。预测反滤波输入某一滤波器的输出后，希望得到它的输入，即不可预测部分，或预测误差。故预测反滤波又称为预测误差滤波。 a-1个 3) 用预测反滤波提高分辨率 二个假设： B、 假设地震子波是最小相位的。 A、 反射系数序列R(t)是随机的白噪序列 由A可以知道反射系数序列，即理想地震记录为不可预测部分。 由B可以得到稳定的、物理可实现的预测反滤波因子。 预测步长越小，地震子波的压缩越大。因此，预测步长越小，预测反滤波所提高的分辨率越高。最小的预测步长是一个采样间隔，分辨率最高。 4) 用预测反滤波消除海上呜震干扰 呜震干扰是一种多次波，其间隔为τH。因此在预测反滤波中选择 a = τH = 2h海/v海，求出预测反滤波因子a″( t )，与地震记录褶积即可消除海上呜震。 呜震形成的过程可以认为是一次滤波，预测反滤波就是“反”呜震滤波的反滤波。 a=12ms a=32ms a=64ms a=128ms 速度是地震勘探中一个十分重要的参数：动校正、偏移、时深转换等处理都以它为参数，它还可以直接用来进行地质构造以及地层岩性的解释。 获得速度参数的途径主要是： ① 利用地震测井、声波测井等，但这种方法是以具备深井为条件，只能在井中进行，故受到限制，并且获得的速度资料只是点上的，面上的资料难以求取。 ② 利用多次覆盖的资料通过速度分析的方法，来求取速度参数。它所求取的是迭加速度，通过它可进一步求取层速度等资料。 定义：速度谱的概念是仿照频率谱的概念而来的。频率谱表示波的能量相对频率的变化规律，因此，我们把地震波的能量相对波的传播速度的变化规律称为速度谱。 1、速度分析准则 水平层状均匀介质时，反射波旅行时为 正常时差－由炮检距不同而引起的时差为 ?t中包函速度V?的信息。 (1) 平均振幅能