地震资料数字处理实验部分课件.ppt

1、地震资料数字处理实验部分二、实验目的二、实验目的 1、了解地震记录格式（SEG-Y）2、掌握Matlab中用于地震数据分析的有关函数难点：程序设计 重点：3D观测系统CDP面元生成和Inline、crossline道头生成3、掌握3D观测系统CDP面元生成和Inline、crossline道头生成地震数据文件格式,由勘探地球物理学家协会SEG(Society of Exploration Geophysicists)制定。(全世界统一的地震数据格式)1、野外(时序)记录格式:SEG A；SEG B；SEG C；SEG D2、室内(道序)记录格式:SEG Y (解编后的炮集记录、CMP道集、水平

2、叠加、叠前/叠后偏移 均为此格式)SEG_Y 文件结构：.End of file EOF2)IBM实数格式：1)IEEE实数格式：fIEEEFes.02)1()127(fIBMFes.016)1()64(Where:s is the sign bit=0(+)/1(-)e is the exponent=(p+bias)f is the fraction.bias=127(IEEE)or 64(IBM)（1）数据格式（2）（text header）40行说明信息（实例）：石东1井西区准格尔盆地腹部石东1井西处理类型:偏移处理系统:FOCUS 4.3客户：新疆石油管理局采样间隔2ms,记录长度6

3、000ms空间位置三维(3)（Reel Header)卷头信息(400 bytes):（4）（Trace Header）道头信息(240 bytes)线内道顺序号 卷内道顺序号 野外记录号（inline）原始野外记录中的道号 炮点号 CDP号 (xline)CDP道集内的道序号 道识别码 l=地震数据；2=死道；3=无效道(空道）;4=爆炸信号；5=井口道；6扫描道；7计时信号；垂直叠加道数 水平叠加道数 数据类型：1=生产；2=试验炮点到接收点的距离(如果排列与激发前进方向相反取负值)(分米)接收点的地面高程。高于海平面为正，低于海平面为负(cm)炮点的地面高程(cm)炮井深度(正数，cm)

4、接收点基准面高程(cm)炮点基准面高程(cm)炮点的水深(cm)接收点的水深(cm)为（41一68字节中）高程和深度的比例因子=l，土10，土100，土1000或者 土10000。如果为正，乘以因子；如果为负，则除以因子为（7388字节中）坐标的比例因子=1，土10，土100，土1000如果为正，乘以因子；如果为负，除以因子。每个样值可用4字节实数表示（5）小结ILINE NO,189,4,n/a,In-line number,.XLINE NO,193,4,n/a,Cross-line number,（1）读SEGY格式地震数据及显示地震剖面 f1=f:/seismic_data/js02_

5、L126_stack.segy;%文件路径及文件名%地震数据以单精度存放在sa.traces(样点，道号)中 sa=single(read_segy_file(f1,times,100,2200,traces,1:150);%读地震数据至结构体变量sa中 s_wplot(sa,annotation,cdp);%显示地震剖面s_ispectrum(sa,plottype,wiggle,annotation,cdp,frequencies,0,80);（2）Interactive spectrum computation(交互谱计算)（2）Interactive spectrum computat

6、ion(交互谱计算)fsa=s_filter(sa,ormsby,5,10,20,40);%滤波 s_compare(sa,fsa,annotation,cdp);%比较显示（3）带通滤波及比较显示（4）显示卷头、道头、拉平（5）AGC、振幅谱（6）地震包函数(S4M)（6）地震包函数（6）地震包函数（6）地震包函数（6）地震包函数（6）地震包函数（6）测井包函数（6）测井包函数（6）地震包函数（7）读测井数据（Las 文件）、显示测井曲线1、了解TAUP变换的原理和方法2、掌握Matlab中用于TAUP变换的有关函数难点：程序设计 TAUP变换变换3、TAUP变换前后地震数据的对比分析实验目

7、的实验目的TAUP变换变换上图中，当波到达A点时，在地面上所观察的早已经通过了B点，和地下C点了。很明显，地面和地下所观察到的速度比波本身的速度要快，如果我们知道从震源O出发的波面到达B点和C点的时间，那我们能知道从O点到A点的时间吗？直线映射为点直线映射为点双曲线映射为椭圆双曲线映射为椭圆nTau=(Tau_max-Tau_min)/DTau+1;%nP=(P_max-P_min)/DP+1;%TP(1:nTau,1:nP)=0;for indexP=1:nP P=P_min+DP*(indexP-1);for indexTau=1:nTau Tau=Tau_min+DTau*(indexT

8、au-1);for indextrace=1:trace_num t=Tau+P*(indextrace-1)*Dx;if(t=interest_tmin)indext=fix(t/Dt)+1;TP(indexTau,indexP)=TP(indexTau,indexP)+Rec(indext,indextrace);时间域代码实现示意时间域代码实现示意工业实现一般用频率域，效率高，但实现较复杂工业实现一般用频率域，效率高，但实现较复杂输入地震数据输入地震数据 TAUP变换后变换后一、实验目的一、实验目的让让学生熟悉掌握地震数据处理中学生熟悉掌握地震数据处理中的的速度分析、速度分析、动校正、动

9、校正、剩余静校正剩余静校正叠加的方法的使用叠加的方法的使用速度扫描实验速度扫描实验二、实验原理二、实验原理地震数据处理流程图三、主要仪器设备三、主要仪器设备 SUN计算机、Linux操作系统、地震处理软件系统四、实验内容与要求四、实验内容与要求实验内容实验内容:1、地震数据的速度分析。2、CMP道集动校正的应用。3、剩余静校正量的计算与应用。4、水平叠加处理实验要求实验要求:运用地震数据处理软件对所给地震数据进行抽道、速度分析、动校正、剩余静校正以及叠加处理。实验完成后，要求学生严格按照实验报告有关编写格式进行实验报告编写，提供纸质和电子版报告。五、处理模块及参数描述五、处理模块及参数描述 经

10、去噪和反褶积处理后的二维地震数据：测线一条，炮数99炮，每炮120道，中间放炮。速度分析：为动校正、偏移、时深转换等处理提供速度参数；为岩性分析、物性分析、油气预测等提供速度信息；为构造解释提供速度参数。动校正：反射波时距曲线上非零炮检距上的波至时间与零炮检距道上的时间之差称为改道的正常时差，为校正正常时差所用的速度称为动校正速度。动校正是用来消除炮检距对反射波旅行时的影响。剩余静校正：野外静校正后，在地震数据中仍然残留有各种剩余静态时移，通常在叠前必需估计出这类剩余静态时移值，并在CMP道集中加以校正，这种校正称为剩余静校正。水平叠加：水平叠加是将CMP道集记录经NMO校正后叠加起来，目的是

11、压制随机噪音，提高地震信噪比。六、处理效果分析六、处理效果分析1、地震数据的速度分析2、CMP道集动校正的应用道集动校正的应用3、剩余静校正、剩余静校正 4、水平叠加处理、水平叠加处理总结总结 通过本次地震数据处理课程设计，对地震数据处理的流程，包括速度分析、动校正、剩余静校正及叠加的方法的使用有了更深一步的理解；将理论和实践相结合，熟悉了linux系统的操作环境和promax软件的基本操作方法，提高了理论联系实际的能力。在处理的过程中要仔细认真的做好每一步，避免疏忽错误，得到较大误差数据，不然实验结果就不准确，影响课程设计的效果。从各个方面提高了自己，不仅对上机过程遇到的问题，还是一些理论知识的疏漏，在这次试验中得到了检验，并学会该如何处理问题，也提高了自己的实际操作动手能了，给以后的实际工作打好基础。