Omega陆上三维地震处理

塔里木油田公司勘探事业部OMEGA陆上三维地震资料处理手册二OO七年四月目录O、引言………………………………………………………1一、OMEGA地震资料处理系统界面………………………1二、OMEGA三维地震资料处理流程举例…………………32.1解编磁带……………………………………………………………32.2观测系统定义……………………………………………………142.3预处理（道头修改、网格定义、噪声压制、几何扩散补偿）…252.4地表一致性振幅补偿和地表一致性反褶积……………………402.5预测反褶积………………………………………………………692.6抽CMP道集……………………………………………………762.7预叠加获取面元属性……………………………………………812.8速度分析…………………………………………………………882.9速度内插…………………………………………………………1262.10第一次剩余静校正………………………………………………1332.11二次速度分析……………………………………………………1922.12第二次剩余静校正（+叠加）……………………………………2052.13叠后时间偏移……………………………………………………223三、处理剖面对比显示……………………………………260OMEGA陆上三维地震资料处理O、引言第1页共261页O、引言Omega地震资料处理系统由西方地球物理公司研究开发，其功能强大，各种功能模块较为齐全，但是参数设置极其繁多，模块搭配，参数设置的原则及技巧更是难以灵活掌握，对初学者来说是非常头痛的事情。鉴于目前除了Omega自带的说明书之外（该说明书还不能导出到PC机里），不管是在网上还是在其他地方有关Omega处理系统的学习资料几乎不能找到，因此编写本手册。通过三维处理实例简单介绍了Omega三维地震资料处理流程，常用处理模块及其参数设置。主要目的是在以后遗忘之时作为唤醒记忆之用，亦可供初学者参考，特别是其中好些参数设置可直接照搬，以期达到事半功倍之效。一、OMEGA地震资料处理系统运行界面其中常用的有：项目管理器作业流程编辑多步作业编辑OMEGA陆上三维地震资料处理一、OMEGA地震资料处理系统运行界面第2页共261页地震数据加载与显示交互速度分析观测系统质量控制网格显示作业运行分布情况显示查看数据信息（ID头、历史和道头）地震数据属性显示质量控制图形显示OMEGA陆上三维地震资料处理二、OMEGA三维地震资料处理流程举例第3页共261页二、OMEGA三维地震资料处理流程举例下面以处理TZ54井南第30束线地震资料为例说明Omega三维地震资料处理流程。2.1解编磁带2.1.1磁带数据输入图2-1磁带解编流程2.1解编磁带2.1.1磁带数据输入第4页共261页（1）INPUT_SEG_D图2-2磁带输入。1-9为磁带号。二、OMEGA三维地震资料处理流程举例第5页共261页图2-3磁带输入2.1解编磁带2.1.1磁带数据输入第6页共261页图2-4磁带输入。610为需要处理的炮数，可大于实际数。二、OMEGA三维地震资料处理流程举例第7页共261页图2-5磁带输入。22和44等分别为每盒磁带最后一炮炮号。2.1解编磁带2.1.1磁带数据输入第8页共261页（2）OUTPUT图2-6数据输出二、OMEGA三维地震资料处理流程举例第9页共261页（3）TRACE_SUMMARY图2-7TRACESUMMARY二、OMEGA三维地震资料处理流程举例第10页共261页2.1.2将解编后数据输出为SEGY格式图2-8输出SEGY格式流程2.1解编磁带2.1.2将解编后数据输出为SEGY格式第11页共261页（1）INPUT图2-9输入解编后的数据二、OMEGA三维地震资料处理流程举例第12页共261页图2-10输入解编后的数据2.1解编磁带2.1.2将解编后数据输出为SEGY格式第13页共261页（2）OUTPUT_SEG_Y图2-11SEGY格式输出二、OMEGA三维地震资料处理流程举例第14页共261页2.2观测系统定义图2-12观测系统定义流程（1）（2）（3）2.2观测系统定义第15页共261页（1）SURVEY_IN（1）图2-13输入关系文件二、OMEGA三维地震资料处理流程举例第16页共261页图2-14输入关系文件。方框为关系文件中各数据所占列数。2.2观测系统定义第17页共261页图2-15输入关系文件二、OMEGA三维地震资料处理流程举例第18页共261页图2-16输入关系文件。MaximumTraceNumber为需要处理的总的道数，这里实际应该等于216*12=2592，多填也不会出错。2.2观测系统定义第19页共261页（2）SURVEY_IN（2）图2-17输入炮点文件。注意方框内设置与SURVEY_IN的不同。二、OMEGA三维地震资料处理流程举例第20页共261页图2-18输入炮点文件2.2观测系统定义第21页共261页图2-19输入炮点文件二、OMEGA三维地震资料处理流程举例第22页共261页（3）SURVEY_IN（3）图2-10输入检波点文件。注意方框内设置与SURVEY（1）（2）的不同。2.2观测系统定义第23页共261页图2-21输入检波点文件二、OMEGA三维地震资料处理流程举例第24页共261页图2-22输入检波点文件二、OMEGA三维地震资料处理流程举例第25页共261页2.3预处理（道头修改、网格定义、噪声压制、几何扩散补偿）图2-23预处理流程（1）（2）2.3预处理第26页共261页（1）INPUT图2-24输入数据二、OMEGA三维地震资料处理流程举例第27页共261页（2）GEOMETRY_UPDATE图2-25置道头2.3预处理第28页共261页（3）TRACE_HEADER_MODIFY（1）图2-26修改道头字。将炮点和检波点野外静校正量分别赋给STACOR_SOURCE和STACOR_DETECT进行静校正。修改所用公式为：CLFLSFLP221110二、OMEGA三维地震资料处理流程举例第29页共261页（4）TRACE_HEADER_MODIFY（2）图2-27将野外文件号赋给IDENT_NUM2.3预处理第30页共261页（5）STRTIMES图2-28定义开始时间。这里为自动搜索，也可由START_TIMES定义。二、OMEGA三维地震资料处理流程举例第31页共261页（6）GRID_DEFINE图2-29定义地下网格2.3预处理第32页共261页图2-30求解地下网格边框求解地下网格边框坐标程序doubleprecisionxr1,yr1,xr2,yr2,xr3,yr3,xr4,yr4doubleprecisionxs1,ys1,xs2,ys2,xs3,ys3,xs4,ys4doubleprecisionxm1,ym1,xm2,ym2,xm3,ym3,xm4,ym4doubleprecisionk1,k3,b1,b2,b3,b4,y2a,y2bxr1=4355405.5yr1=679320.9xr2=4374411.5yr2=693563.2xr3=4370640.0yr3=698606.9xr4=4351613.5yr4=684387.6xs1=4358995.5ys1=683743.6xs3=4367247.0ys3=694166.1k1=(yr2-yr1)/(xr2-xr1)xm1=(xr1+xs1)/2.0（xr1,yr1）（xr2,yr2）（xr3,yr3）（xr4,yr4）（xs1,ys1）（xs3,ys3）二、OMEGA三维地震资料处理流程举例第33页共261页ym1=(yr1+ys1)/2.0b1=ym1-k1*xm1k3=(yr2-yr3)/(xr2-xr3)xm3=(xr3+xs3)/2.0ym3=(yr3+ys3)/2.0b3=ym3-k3*xm3xm2=(b3-b1)/(k1-k3)y2a=k1*xm2+b1y2b=k3*xm2+b3ym2=(y2a+y2b)/2b2=ym3-k1*xm3b4=ym1-k3*xm1xm4=(b2-b4)/(k3-k1)y2a=k1*xm4+b2y2b=k3*xm4+b4ym4=(y2a+y2b)/2write(*,*)'xm1=',xm1,'ym1=',ym1write(*,*)'xm2=',xm2,'ym2=',ym2write(*,*)'xm3=',xm3,'ym3=',ym3write(*,*)'xm4=',xm4,'ym4=',ym4open(2,file='gridm')write(2,'(2f10.1)')xm1,ym1,xm2,ym2,xm3,ym3,xm4,ym4close(2)stopend2.3预处理第34页共261页图2-31地面和地下网格边框二、OMEGA三维地震资料处理流程举例第35页共261页（7）WINDOW_DESIGN图2-32定义面波区域2.3预处理第36页共261页（8）ZONE_FILTER图2-33区域滤波衰减面波二、OMEGA三维地震资料处理流程举例第37页共261页图2-34区域滤波衰减面波2.3预处理第38页共261页（9）GEOSPREAD_COMP图2-35几何扩散补偿二、OMEGA三维地震资料处理流程举例第39页共261页（10）OUTPUT图2-36输出数据供下一步反褶积使用二、OMEGA三维地震资料处理流程举例第40页共261页2.4地表一致性振幅补偿和地表一致性反褶积图2-37地表一致性振幅补偿和地表一致性反褶积流程2.4地表一致性振幅补偿和地表一致性反褶积第41页共261页2.4.1地表一致性振幅补偿图2-38地表一致性振幅补偿流程地表一致性振幅补偿模块组合为的是补偿数据体的由炮点、检波点、偏移距的不同引起的振幅变化。SCAC处理后，任意给定的炮点、检波点、偏移距的振幅应该与其他的炮点、检波点、偏移距的振幅水平相应一致。注意SCAC处理不改变资料的信噪比。假如炮A基本上是噪音而炮B基本上是信号，则SCAC后，炮A噪音的振幅水平与炮B的振幅水平一致。二、OMEGA三维地震资料处理流程举例第42页共261页（1）INPUT图2-39输入数据2.4.1地表一致性振幅补偿第43页共261页图2-40输入数据二、OMEGA三维地震资料处理流程举例第44页共261页（2）WINDOW_DESIGN图2-41定义地表一致性振幅补偿时窗。该时窗由两条双曲线组成，而非矩形。每条双曲线可由两个参数确定：Zero-OffsetBeginTime和BeginTimeMoveoutVelocity。2.4.1地表一致性振幅补偿第45页共261页（3）SCAC_PICK图2-42对给定的时窗计算均方根振幅或平均振幅二、OMEGA三维地震资料处理流程举例第46页共261页（4）SCAC_DECOMPOSITION图2-43将SCAC_PICK模块中计算的振幅地表一致性地分解到炮域、检波点域和偏移距域2.4.1地表一致性振幅补偿第47页共261页（5）SCAC_APPLY图2-44使用SCAC_DECOMPOSITION的输出且求出并应用调整每道振幅的补偿比例因子，使所有道集域（炮点、检波点、偏移距）有相同的均方根振幅水平。二、OMEGA三维地震资料处理流程举例第48页共261页（6）OUTPUT图2-45输出地表一致性补偿后的数据二、OMEGA三维地震资料处理流程举例第49页共261页2.4.2地表一致性反褶积图2-46地表一致性反褶积流程2.4.2地表一致性反褶积2.4.2.1算子设计第50页共261页2.4.2.1地表一致性反褶积算子设计图2-47地表一致性反褶积算子设计流程二、OMEGA三维地震资料处理流程举例第51页共261页（1）INPUT图2-48输入地表一致性振幅补偿后的数据2.4.2地表一致性反褶积2.4.2.1算子设计第52页共261页图2-49输入地表一致性振幅补偿后的数据