地震资料数字处理实验部分

地震资料数字处理实验部分二、实验目的1、了解地震记录格式（SEG-Y）2、掌握Matlab中用于地震数据分析的有关函数难点：程序设计重点：3D观测系统CDP面元生成和Inline、crossline道头生成3、掌握3D观测系统CDP面元生成和Inline、crossline道头生成三、地震记录格式地震数据文件格式,由勘探地球物理学家协会SEG(SocietyofExplorationGeophysicists)制定。(全世界统一的地震数据格式)1、野外(时序)记录格式:SEGA；SEGB；SEGC；SEGD2、室内(道序)记录格式:SEGY(解编后的炮集记录、CMP道集、水平叠加、叠前/叠后偏移均为此格式)SEG_Y文件结构：三、SEG-Y格式......EndoffileEOF三、SEG-Y格式2)IBM实数格式：1)IEEE实数格式：Where:sisthesignbit=0(+)/1(-)eistheexponent=(p+bias)fisthefraction.bias=127(IEEE)or64(IBM)（1）数据格式三、SEG-Y格式（2）（text

header）40行说明信息（实例）：石东1井西区准格尔盆地腹部石东1井西处理类型:偏移处理系统:FOCUS4.3客户：新疆石油管理局采样间隔2ms,记录长度6000ms空间位置三维三、SEG-Y格式(3)（ReelHeader)卷头信息

(400bytes):三、SEG-Y格式三、SEG-Y格式（4）（TraceHeader）道头信息

(240bytes)线内道顺序号

卷内道顺序号

野外记录号（inline）原始野外记录中的道号炮点号

CDP号(xline)

CDP道集内的道序号道识别码l=地震数据；2=死道；3=无效道(空道）;4=爆炸信号；5=井口道；6＝扫描道；7＝计时信号；三、SEG-Y格式垂直叠加道数

水平叠加道数

数据类型：1=生产；2=试验炮点到接收点的距离(如果排列与激发前进方向相反取负值)(分米)接收点的地面高程。高于海平面为正，低于海平面为负(cm)炮点的地面高程(cm)炮井深度(正数，cm)接收点基准面高程(cm)炮点基准面高程(cm)炮点的水深(cm)接收点的水深(cm)三、SEG-Y格式为（41一68字节中）高程和深度的比例因子=l，土10，土100，土1000或者土10000。如果为正，乘以因子；如果为负，则除以因子为（73—88字节中）坐标的比例因子=1，土10，土100，土1000如果为正，乘以因子；如果为负，除以因子。三、SEG-Y格式三、SEG-Y格式三、SEG-Y格式三、SEG-Y格式每个样值可用4字节实数表示三、SEG-Y格式（5）小结{`ILINENO',189,4,`n/a',`In-linenumber'},...{`XLINENO',193,4,`n/a',`Cross-linenumber'}},四、S4m中地震数据分析常用函数（1）读SEGY格式地震数据及显示地震剖面

f1=‘f:/seismic_data/js02_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）Interactivespectrumcomputation(交互谱计算)（2）Interactivespectrumcomputation(交互谱计算)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变换前后地震数据的对比分析实验目的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;forindexP=1:nPP=P_min+DP*(indexP-1);forindexTau=1:nTauTau=Tau_min+DTau*(indexTau-1);forindextrace=1:trace_numt=Tau+P*(indextrace-1)*Dx;if(t<=interest_tmax&&t>=interest_tmin)indext=fix(t/Dt)+1;TP(indexTau,indexP)=TP(indexTau,indexP)+Rec(indext,indextrace);

时间域代码实现示意工业实现一般用频率域，效率高，但实现较复杂输入地震数据

TAUP变换后一、实验目的让学生熟悉掌握地震数据处理中的速度分析、动校正、剩余静校正叠加的方法的使用速度扫描实验二、实验原理地震数据处理流程图三、主要仪器设备SUN计算机、Linux操作系统、地震处理软件系统四、实验内容与要求实验内容:1、地震数据的速度分析。2、CMP道集动校正的应用。3、剩余静校正量的计算与应用。4、水平叠加处理实验要求:

运用地震数据处理软件对所给地震数据进行抽道、速度分析、动校正、剩余静校正以及叠加处理。实验完成后，要求学生严格按照实验报告有关编写格式进行实验报告编写，提供纸质和电子版报告。五、处理模块及参数描述

经去噪和反褶积处理后的二维地震数据：测线一条，炮数99炮，每炮120道，中间放炮。

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

4、水平叠加处理总结

通过本次地震数据处理课程设计，对地震数据处理的流程，包括速度分析、动校正、剩余静校正及叠加的方法的使用有了更深一步的