地震资料采集与处理

序号成绩

中国地质大学（）本科生

实验报告

《地震资料采集与处理》上机实验报告

姓名：建明

班级： 061154

学号：

指导老师：卞爱飞

小组成员：建明，朴青峰

完成日期： 2018年5月11日

目录

1.一维带通滤

波………………………………………………………………………………………………

(1)

2. 动校正与叠加 (10)

3. 偏移算子点脉冲响应 (14)

4. 叠后数据偏移 (17)

5. 总结 (21)

1.一维带通滤波实验

1.1.实验目的

利用一维频率域滤波方法分析实际地震资料中有效信号与干扰波的时空分布特征，掌握低通、带通、高通滤波器的设计方法和相关SU模块的调用方法，设计频率域滤波器进行有效信号与噪音的分离，对滤波前后地震剖面进行处理效果对比显示，分析一维滤波方法的优缺点。

1.2.基本原理

本实验核心处理模块为sufilter常用的模块调用方法为：

sufilter fileout.su [ f=x1,x2,x3,x4 amps=y1,y2,y3,y4 ] & 吧

其中 sufilter 为调用模块名称，filein.su 为输入的 SU 格式时间域多道地震信号文件名，fileout.su 为处理输出的 SU 格式时间域多道地震信号文件名，f 为频率控制点，amps 为对应频率控制点的振幅值，&表示后台运行，[]表示方括号的参数有默认值，可选填。对于标准的频率域带通滤波器，4个控制频点及相应频点振幅谱即可确定一个带通滤波器的形态。

1.3.实验步骤

（1）在当前控制台输入’cd $CWPROOT/demos/nmo’命令进入SU动校正与叠加演示目录，输入命令‘sh pre.sh学号’学号为个人实际学号。

（2）合成演示数据。

（3）原始数据显示。

（4）原始数据增益补偿。

（5）原始数频谱分析。

（6）对比不同频带信号特征。

（7）设计带通滤波器

1.4.实验结果与分析

1.4.1.原始数据

原始数据的图像如图1-1所示，振幅只有在中间的上半部分比较明显，其他地方振幅基

本为0，通过显示模块suximage后得到的图像如图1-2—图1-5所示，从图1-2可以看出在中间部分存在比较明显的扫帚状的面波噪声， perc的选择并不是越大越好，在perc=90时，数据显示最清晰，继续增大，数据显示效果降低。

图1-1：原始数据图1-2：原始数据（perc=85）

图1-3：原始数据（perc=90）图1-4：原始数据（perc=95）

图1-3：原始数据（perc=90

1.4.2原始数据增益补偿

原始数据的图像如图1-1所示，通过增益显示后得到的图像如图1-6—图1-9所示，从图中可以看出在中间部分存在比较明显的扫帚状的面波噪声，地震波能量明显减弱，不同的Tpow值对不同深度的增益补偿效果不同，Tpow值越小，浅部的能量突出效果明显，Tpow值越大，深部的能量突出效果越明显。

图1-6 增益补偿（tpow=1）图1-7 增益补偿（tpow=2）

图1-8 增益补偿（tpow=3）图1-9 增益补偿（tpow=4）

1.4.3原始数据频谱分析

原始数据频谱分析如图1-10所示，从图中可以看出能量主要集中在频率为5-40Hz的低频围，在60-80Hz的频率分布有零星的能量，可能是在对应位置下面存在波阻抗较大的界面，并且随着perc值得增大，浅部能量也被压制。

图1-10频谱分析（perc分别为95,97,99,99.9）

1.4.4低通滤波

对其进行低通滤波（sufilter模块，f=0,15,40,50，amps=1,1,0,0）得到的图像如图1-11所示，与原始图像相差不大，中间的面波噪声依然没有过滤掉，且浅层反射的振幅减弱，分析其原因，主要是因为浅层反射波频率高，被过滤掉了。

图1-11：低通滤波图1-12：高通滤波

1.4.5高通滤波

设计不同的高通滤波器对其进行高通滤波，得到的图像如图1-12所示，通过比较，高通滤波对面波的压制作用较大。

1.1.6带通滤波

设计带通滤波器，得到的图像如图1-13所示，其中滤波器参数为f=15，20，50,70，amps=0,1,1,0时，带通滤波结果可以基本过滤掉面波，得到的带通滤波的结果中在中间某些部分振幅减弱，图像变得不连续。

图1-13：带通滤波

1.1.7带通滤波设计

设计带通滤波器，其中滤波器参数如下图（图1-14），滤波结果如图1-15，其中参数为f=19，20，30,35，amps=0,1,1,0时，带通滤波结果可以基本过滤掉面波，得到的带通滤波的结果中在中间某些部分振幅减弱，图像变得不连续。

图1-14 带通滤波器参数

图1-15 带通滤波结果

2.动校正与叠加

2.1.实验目的

利用合成地震数据制作速度谱进行交互速度分析，掌握SU动校正与叠加模块的原理与使用方法，分析动校正前后CDP道集变化特征与叠加效果。

2.2.基本原理

动校正核心sunmofileout.su par=parf&

其中sunmo为调用模块名称，filein.su为输入的SU格式CMP道集文件名，fileout.su为处理输出的SU格式CMP道集文件名，parf为速度交互分析拾取得到动校正速度参数文件，&表示后台运行。

叠加核心Sustackfileout.su &

其中sustack为调用模块名称，filein.su为输入SU格式动校正后CMP道集文件名，fileout.su为处理输出的SU格式CMP道集叠加文件名，&表示后台运行。

2.3.实验步骤

（1）在当前控制台输入’cd $CWPROOT/demos/nmo’命令进入SU动校正与叠加演示目录，输入命令‘sh pre.sh学号’学号为个人实际学号

（2）层状模型数据合成。在当前控制台输入命令‘sh 学号_job1.sh学号’生成演示模型数据并显示。

（3）叠前炮集数据合成。在当前控制台输入命令‘sh 学号_job2.sh学号’得到合