处理系统学习手册陆上三维地震

1、塔里木油田公司勘探事业部OMEGA资料处理手册陆上三维二 OO 七年四月目录O、引 言1一、OMEGA资料处理系统界面1二、OMEGA资料处理流程举例3三维2.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、引言O、引

2、 言Omega资料处理系统由西方地球物理公司研究开发，其功能强大，各种功能模块较为齐全，但是参数设置极其繁多，模块搭配，参数设置的原则及技巧更是难以灵活掌握，对初学者来说是非常头痛的事情。鉴于目前除了 Omega 自带的说明书之外（该说明书还不能导出到 PC 机里），不管是在网上还是在其他地方有关 Omega 处理系统的学习资料几乎不能找到，因此编写本手册。通过三维处理实例简单介绍了 Omega 三维资料处理流程，常用处理模块及其参数设置。主要目的是在以后遗忘之时作为唤醒记忆之用，亦可供初学者参考，特别是其中好些参数设置可直接照搬，以期达到事半功倍之效。一、OMEGA资料处理系统运行界面其中常

3、用的有：项目管理器作业流程编辑多步作业编辑第 1 页 共 261 页OMEGA 陆上三维资料处理一、OMEGA资料处理系统运行界面数据加载与显示交互速度分析观测系统质量网格显示作业运行分布情况显示查看数据信息（ID 头、历史和道头）数据属性显示质量图形显示第 2 页 共 261 页OMEGA 陆上三维资料处理二、OMEGA 三维资料处理流程举例二、OMEGA 三维资料处理流程举例下面以处理 TZ54 井南第 30 束线理流程。资料为例说明 Omega 三维资料处2.1解编磁带2.1.1 磁带数据输入图 2-1 磁带解编流程第 3 页 共 261 页2.1 解编磁带2.1.1 磁带数据输入（1）

4、INPUT_SEG_D图 2-2 磁带输入。1-9 为磁带号。第 4 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-3磁带输入第 5 页 共 261 页2.1 解编磁带2.1.1 磁带数据输入图 2-4 磁带输入。610 为需要处理的数，可大于实际数。第 6 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-5 磁带输入。22 和 44 等分别为每盒磁带最后一号。第 7 页 共 261 页2.1 解编磁带2.1.1 磁带数据输入（2）OUTPUT图 2-6 数据输出第 8 页 共 261 页二、OMEGA 三维资料处理流程举例（3）TRACE_SUMMARY图 2-

5、7TRACE SUMMARY第 9 页 共 261 页二、OMEGA 三维资料处理流程举例2.1.2将解编后数据输出为 SEGY 格式图 2-8输出SEGY 格式流程第 10 页 共 261 页2.1 解编磁带2.1.2 将解编后数据输出为 SEGY 格式（1）INPUT图 2-9 输入解编后的数据第 11 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-10输入解编后的数据第 12 页 共 261 页2.1 解编磁带2.1.2 将解编后数据输出为 SEGY 格式（2）OUTPUT_SEG_Y图 2-11 SEGY 格式输出第 13 页 共 261 页二、OMEGA 三维资料处理

6、流程举例2.2观测系统定义（1）（2）（3）图 2-12观测系统定义流程第 14 页 共 261 页2.2 观测系统定义（1）SURVEY_IN（1）图 2-13 输入关系文件第 15 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-14 输入关系文件。方框为关系文件中各数据所占列数。第 16 页 共 261 页2.2 观测系统定义图 2-15 输入关系文件第 17 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-16 输入关系文件。um Trace Number 为需要处理的总的道数，这里实际应该等于 216*12=2592，多填也出错。第 18 页 共 261

7、页2.2 观测系统定义（2）SURVEY_IN（2）图 2-17 输入 点文件。注意方框内设置与 SURVEY_IN 的不同。第 19 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-18输入 点文件第 20页 共 261 页2.2 观测系统定义图 2-19 输入 点文件第 21 页 共 261 页二、OMEGA 三维资料处理流程举例（3）SURVEY_IN（3）图 2-10 输入检波点文件。注意方框内设置与 SURVEY（1）（2）的不同。第 22 页 共 261 页2.2 观测系统定义图 2-21 输入检波点文件第 23 页 共 261 页二、OMEGA 三维资料处理流程举例

8、图 2-22输入检波点文件第 24 页 共 261 页二、OMEGA 三维资料处理流程举例2.3预处理（道头修改、网格定义、噪声压制、几何扩散补偿）（1）（2）图 2-23 预处理流程第 25 页 共 261 页2.3 预处理（1）INPUT图 2-24 输入数据第 26 页 共 261 页二、OMEGA 三维资料处理流程举例（2）GEOMETRY_UPDATE图 2-25 置道头第 27 页 共 261 页2.3 预处理（3）TRACE_HEADER_MODIFY（1）图 2-26 修改道头字。将点和检波点野外静校正量分别赋给 STACOR_SOURCE 和 STACOR_DETECT 进行

9、静校正。修改所用公式为：L = F × S × LP + F × L + C011122第 28 页 共 261 页二、OMEGA 三维资料处理流程举例（4）TRACE_HEADER_MODIFY（2）图 2-27 将野外文件号赋给 IDENT_NUM第 29 页 共 261 页2.3 预处理（5）STRTIMES图 2-28 定义开始时间。这里为自动搜索，也可由 START_TIMES 定义。第 30 页 共 261 页二、OMEGA 三维资料处理流程举例（6）GRID_DEFINE图 2-29定义网格第 31页 共 261 页2.3 预处理（xr2,yr2）（

10、xr3,yr3）（xs3,ys3）（xs1,ys1）（xr1,yr1）（xr4,yr4）图 2-30 求解网格边框求解网格边框坐标程序double precision xr1,yr1,xr2,yr2,xr3,yr3,xr4,yr4 double precision xs1,ys1,xs2,ys2,xs3,ys3,xs4,ys4double precision xm1,ym1,xm2,ym2,xm3,ym3,xm4,ym4 double precision k1,k3,b1,b2,b3,b4,y2a,y2bxr1 = 4355405.5yr1 = 679320.9xr2 = 4374411.5y

11、r2 = 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第 32 页 共 261 页二、OMEGA 三维资料处理流程举例ym1 = (yr1 + ys1) /2.0 b1= ym1-k1*xm1k3= (yr2-yr3)/(xr2-xr3) xm3 = (xr3 + xs3) /2.0 ym3 = (yr

12、3 + ys3) /2.0b3= ym3-k3*xm3 xm2 = (b3-b1)/(k1-k3) y2a = k1*xm2+b1y2b = k3*xm2+b3 ym2= (y2a+y2b)/2b2 =ym3-k1*xm3 b4 =ym1-k3*xm1xm4 = (b2-b4)/(k3-k1) y2a = k1*xm4+b2y2b = k3*xm4+b4 ym4 = (y2a+y2b)/2write(*,*) 'write(*,*) 'write(*,*) 'write(*,*) 'xm1 xm2 xm3 xm4=', xm1 ,'=',

13、 xm2 ,'=', xm3 ,'=', xm4 ,'ym1 ym2 ym3 ym4=', ym1=', ym2=', ym3=', ym4open (2,file='gridm')write(2,'(2f10.1)') xm1,ym1,xm2,ym2,xm3,ym3,xm4,ym4 close(2)stop end第 33 页 共 261 页2.3 预处理图 2-31 地面和网格边框第 34 页 共 261 页二、OMEGA 三维资料处理流程举例（7）WINDOW_DESIGN图 2-32

14、定义面波区域第 35页 共 261 页2.3 预处理（8）ZONE_FILTER图 2-33 区域滤波衰减面波第 36 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-34区域滤波衰减面波第 37 页 共 261 页2.3 预处理（9）GEOSPREAD_COMP图 2-35 几何扩散补偿第 38 页 共 261 页二、OMEGA 三维资料处理流程举例（10）OUTPUT图 2-36 输出数据供下一步反褶积使用第 39 页 共 261 页二、OMEGA 三维资料处理流程举例2.4地表一致性振幅补偿和地表一致性反褶积图 2-37 地表一致性振幅补偿和地表一致性反褶积流程第 40 页

15、 共 261 页2.4 地表一致性振幅补偿和地表一致性反褶积2.4.1 地表一致性振幅补偿图 2-38地表一致性振幅补偿流程地表一致性振幅补偿模块组合为的是补偿数据体的由点、检波点、偏移距的不同引起的振幅变化。SCAC 处理后，任意给定的点、检波点、偏移距的振幅应该与其他的点、检波点、偏移距的振幅水平相应一致。注意 SCAC 处理不改变资料的信噪比。假如A 基本上是噪音而B 基本上是信号，则 SCAC 后，A 噪音的振幅水平与B 的振幅水平一致。第 41 页 共 261 页二、OMEGA 三维资料处理流程举例（1）INPUT图 2-39 输入数据第 42 页 共 261 页2.4.1 地表一致

16、性振幅补偿图 2-40 输入数据第 43 页 共 261 页二、OMEGA 三维资料处理流程举例（2）WINDOW_DESIGN图 2-41 定义地表一致性振幅补偿时窗。该时窗由两条双曲线组成，而非矩形。每条双曲线可由两个参数确定：Zero-Offset Begin Time 和Begin Time Moveout Velocity。第 44 页 共 261 页2.4.1 地表一致性振幅补偿（3）SCAC_PICK图 2-42对给定的时窗计算均方根振幅或平均振幅第 45 页 共 261 页二、OMEGA 三维资料处理流程举例（4）SCAC_DECOMPOSITION图 2-43 将 SCAC_

17、PICK 模块中计算的振幅地表一致性地分解到域、检波点域和偏移距域第 46 页 共 261 页2.4.1 地表一致性振幅补偿（5）SCAC_APPLY图2-44 使用SCAC_DECOMPOSITION 的输出且求出并应用调整每道振幅的补偿比例因子，使所有道集域（点、检波点、偏移距）有相同的均方根振幅水平。第 47 页 共 261 页二、OMEGA 三维资料处理流程举例（6）OUTPUT图 2-45 输出地表一致性补偿后的数据第 48 页 共 261 页二、OMEGA 三维资料处理流程举例2.4.2 地表一致性反褶积图 2-46地表一致性反褶积流程第 49 页 共 261 页2.4.2 地表一

18、致性反褶积2.4.2.1 算子设计2.4.2.1 地表一致性反褶积算子设计图 2-47 地表一致性反褶积算子设计流程第 50 页 共 261 页二、OMEGA 三维资料处理流程举例（1）INPUT图 2-48 输入地表一致性振幅补偿后的数据第 51 页 共 261 页2.4.2 地表一致性反褶积2.4.2.1 算子设计图 2-49 输入地表一致性振幅补偿后的数据第 52 页 共 261 页二、OMEGA 三维资料处理流程举例（2）TRACE_BALANCE图 2-50 道均衡第 53 页 共 261 页2.4.2 地表一致性反褶积2.4.2.1 算子设计（3）STRTIMES图 2-51 定义

19、开始时间。定义 START TIMES 参数，SEARCH TIME WINDOW 参数无效。第 54 页 共 261 页二、OMEGA 三维资料处理流程举例（4）SC_DCN_SPCTRL_ANL图 2-52 地表一致性反褶析。由该模块得到用于地表一致性反褶积处理的对数功率谱。第 55 页 共 261 页2.4.2 地表一致性反褶积2.4.2.1 算子设计（5）SC_DCN_SPCTRL_DECOMP图 2-53 地表一致性谱分解。该模块应用高斯-塞得尔方法将对数功率谱分解为中心点和偏移距分量。点、检波点、第 56 页 共 261 页二、OMEGA 三维资料处理流程举例（6）SC_DCN_O

20、PR_DESIGN图 2-54 地表一致性反褶积算子设计第 57 页 共 261 页2.4.2 地表一致性反褶积2.4.2.1 算子设计（7）OUTPUT图 2-55 输出反褶积算子第 58 页 共 261 页二、OMEGA 三维资料处理流程举例2.4.2.2 反褶积算子的应用（1）（2）图 2-56 反褶积算子应用流程第 59 页 共 261 页2.4.2 地表一致性反褶积2.4.2.2 算子应用（1）INPUT（1）图 2-57 输入地表一致性振幅补偿后的数据第 60 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-58 输入地表一致性振幅补偿后的数据第 61 页 共 261

21、 页2.4.2 地表一致性反褶积2.4.2.2 算子应用（2）TRACE_BALANCE图 2-59 道均衡第 62 页 共 261 页二、OMEGA 三维资料处理流程举例（3）BPFILTER图 2-60最小相位滤波反褶积的基本假设前提是资料的子波为最小相位，反褶积前滤波做最小相位处理的方法一般是在反褶积前进行一个最小相位滤波，以使资料满足褶积及脉冲反褶积最小相位假设。反目前除了 Omega 处理系统在反褶积前做最小相位滤波（MBPF）之外，其他处理系统通常用零相位滤波。主要考虑如下：如果是老资料处理，不知先前处理是何相位，新旧处理资料层位可能对不上。实践表明，地表一致性反褶积前用最小相位滤

22、波比较好。反褶积前的最小相位滤波选件中，白噪系数对反褶积效果影响很大，建议选0.01%较好。实际数据要根据具体问题进行试验确定（此参数应小于 1%）。第 63 页 共 261 页2.4.2 地表一致性反褶积2.4.2.2 算子应用（4）INPUT（2）图 2-61 输入反褶积算子第 64 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-62输入反褶积算子第 65 页 共 261 页2.4.2 地表一致性反褶积2.4.2.2 算子应用（5）DCN_OPR_APPLY图 2-63 应用反褶积算子第 66 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-64应用反褶积算子

23、第 67 页 共 261 页2.4.2 地表一致性反褶积2.4.2.2 算子应用（6）OUTPUT图 2-65 输出地表一致性反褶积后的数据第 68 页 共 261 页二、OMEGA 三维资料处理流程举例2.5反褶积图 2-66反褶积流程第 69 页 共 261 页2.5反褶积（1）INPUT图 2-67 输入地表一致性反褶积后的数据第 70 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-68 输入地表一致性反褶积后的数据第 71 页 共 261 页2.5反褶积（2）PRD_DCN_SPCTRL_ANL图 2-69反褶析。方框内为所定义的时窗，时窗由零偏移距时间和相应动较正速度

24、（确定一条双曲线）确定。注意，目的层不能划入过渡带之间（比如 2144 和 2678 之间）。第 72 页 共 261 页二、OMEGA 三维资料处理流程举例（3）PRD_DCN_OPR_DESIGN图 2-70反褶积算子设计。方框内分别对应上面 PRD_DCN_SPCTRL_ANL 模块中三个时窗的步长。第 73 页 共 261 页2.5反褶积（4）DCN_OPR_APPLY图 2-71 应用反褶积算子第 74 页 共 261 页二、OMEGA 三维资料处理流程举例（5）OUTPUT图 2-72输出反褶积后的数据第 75 页 共 261 页二、OMEGA 三维资料处理流程举例2.6抽 CMP

25、 道集图 2-73抽CMP 道集流程第 76 页 共 261 页2.6 抽CMP 道集（1）INPUT图 2-74 输入反褶积后的数据第 77 页 共 261 页二、OMEGA 三维资料处理流程举例（2）SORT图 2-75抽CMP 道集第 78 页 共 261 页2.6 抽CMP 道集（3）DATUM_RECOMPUTE图 2-76基准面重新计算第 79 页 共 261 页二、OMEGA 三维资料处理流程举例（4）OUTPUT图 2-77输出CMP 道集第 80 页 共 261 页二、OMEGA 三维资料处理流程举例2.7预叠加获取面元属性图 2-78预叠加获取面元属性流程第 81 页 共

26、261 页2.7 预叠加获取面元属性（1）INPUT图 2-79 输入预处理后的数据第 82 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-80输入预处理后的数据第 83 页 共 261 页2.7 预叠加获取面元属性（2）TRACE_HEADER_MODIFY图 2-81修改道头第 84 页 共 261 页二、OMEGA 三维资料处理流程举例（3）SORT图 2-82抽CMP 道集第 85 页 共 261 页2.7 预叠加获取面元属性（4）STACK图 2-83 叠加第 86 页 共 261 页二、OMEGA 三维资料处理流程举例（5）OUTPUT图 2-84输出叠加后的数据第

27、 87 页 共 261 页二、OMEGA 三维资料处理流程举例2.8 速度分析2.8.1生成交互速度分析数据（2）（1）图 2-85 生成交互速度分析数据流程第 88 页 共 261 页2.8 速度分析2.8.1 生成交互速度分析数据（1）INPUT图 2-86 输入CMP 道集数据第 89 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-87输入CMP 道集数据第90 页 共 261 页2.8 速度分析2.8.1 生成交互速度分析数据（2）VEL_GENERATOR图 2-88 设置VEL_GENERATOR参数。方框内为定义的速度范围。第 91 页共 261 页二、OMEGA

28、 三维资料处理流程举例图 2-89 设置VEL_GENERATOR参数第 92 页 共 261 页2.8 速度分析2.8.1 生成交互速度分析数据图 2-90 设置VEL_GENERATOR 参数第 93 页 共 261 页二、OMEGA 三维资料处理流程举例（3）VELAN_DEFINITION图 2-91 定义要进行速度分析的 CMP 线。40、70 和 100 为选取的 CMP 线；260 为每条CMP线上的首个 CMP 号。该参数的选取由 ATTRIB DISP 加载 2.7 步骤输出数据TZ54S_sw30_bstk0 可得。第 94 页 共 261 页2.8 速度分析2.8.1 生

29、成交互速度分析数据（4）BPFILTER图 2-92 零相位滤波第 95 页 共 261 页二、OMEGA 三维资料处理流程举例（5）INSTANT_GAIN图 2-93 瞬时增益第 96 页 共 261 页2.8 速度分析2.8.1 生成交互速度分析数据（6）GATHER_XCOR图 2-94 道集相关第 97 页 共 261 页二、OMEGA 三维资料处理流程举例（7）ADJACENT_CMP_SUM（1）图 2-95第 98 页 共 261 页2.8 速度分析2.8.1 生成交互速度分析数据（8）NONLINEAR_TO_LINEAR图 2-96第 99 页 共 261 页二、OMEGA

30、 三维资料处理流程举例（9）OUTPUT_IVP图 2-97 输出交互速度分析所用的 SEMBLANCES 数据。注意方框内设置要与交互速度分析时相同。第 100 页 共 261 页2.8 速度分析2.8.1 生成交互速度分析数据（10）TRACE_SELECT图 2-98选取需要进行叠加的 CMP 线第 101 页 共 261 页二、OMEGA 三维资料处理流程举例（11）OUTPUT图 2-99 输出需要叠加的 CMP 线数据第 102 页 共 261 页2.8 速度分析2.8.1 生成交互速度分析数据（12）ADJACENT_CMP_SUM（2）图 2-100ADJACENT_CMP_S

31、UM第 103 页 共 261 页二、OMEGA 三维资料处理流程举例（13）OUT_IVP图 2-101 输出交互速度分析所用的 GATHERS 数据。注意方框内设置要与交互速度分析时相同。该道集数据除了交互分析时要用外，还可供定义切除用。第 104 页 共 261 页二、OMEGA 三维资料处理流程举例2.8.2 交互速度分析（IVP）（略）2.8.3生成叠加剖面对比进行局部速度调整图 2-102 生成叠加剖面对比进行局部速度调整流程第 105 页 共 261 页2.8 速度分析2.8.3 生成叠加剖面对比进行局部速度调整（1）INPUT图 2-103 输入 2.8.1 步骤生成的需要进行

32、叠加的 CMP 线数据第 106 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-104 输入 2.8.1 步骤生成的需要进行叠加的 CMP 线数据第 107 页 共 261 页2.8 速度分析2.8.3 生成叠加剖面对比进行局部速度调整（2）TRACE_SELECT图 2-105 选取某条 CMP 线第 108 页 共 261 页二、OMEGA 三维资料处理流程举例（3）VEL_GENERATOR图 2-106将 txt 速度文件转换成二进制速度文件第 109 页 共 261 页2.8 速度分析2.8.3 生成叠加剖面对比进行局部速度调整图 2-107 将 txt 速度文件转换

33、成二进制速度文件第 110 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-108 将 txt 速度文件转换成二进制速度文件第 111 页 共 261 页2.8 速度分析2.8.3 生成叠加剖面对比进行局部速度调整（4）VEL_SCALE图 2-109 速度调整比例第 112 页 共 261 页二、OMEGA 三维资料处理流程举例（5）NMO图 110 动较正第 113 页 共 261 页2.8 速度分析2.8.3 生成叠加剖面对比进行局部速度调整（6）MUTE图 111 定义切除库。该切除库由在 DATA DRIVEN 中加载 2.8.1 步骤 GATHERS 数据经动校正后切

34、除可得。第 114 页 共 261 页二、OMEGA 三维资料处理流程举例（7）RMS_GAIN图 2-112均方根增益第 115 页 共 261 页2.8 速度分析2.8.3 生成叠加剖面对比进行局部速度调整（8）STACK图 2-113 叠加第 116 页 共 261 页二、OMEGA 三维资料处理流程举例（9）DATUM_CORRECTION图 2-114基准面校正第 117 页 共 261 页2.8 速度分析2.8.3 生成叠加剖面对比进行局部速度调整（10）BPFILTER图 2-115 零相位滤波第 118 页 共 261 页二、OMEGA 三维资料处理流程举例（11）RMS_GA

35、IN图 2-116均方根增益第 119 页 共 261 页2.8 速度分析2.8.3 生成叠加剖面对比进行局部速度调整（12）SEISPLOT图 2-117 生成质量图片第 120 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-118生成质量图片第 121 页 共 261 页2.8 速度分析2.8.3 生成叠加剖面对比进行局部速度调整图 2-119 生成质量图片第 122 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-120生成质量图片第 123 页 共 261 页2.8 速度分析2.8.3 生成叠加剖面对比进行局部速度调整图 2-121 生成质量图片第 124

36、 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-122生成质量图片第 125 页 共 261 页二、OMEGA 三维资料处理流程举例2.9速度内插图 2-123速度内插流程第 126页 共 261 页2.9 速度内插（1）INPUT图 2-124 输入 CMP 道集数据第 127 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-125输入 CMP 道集数据第 128 页 共 261 页2.9 速度内插（2）VEL_GENERATOR图 126 将txt 速度文件转换成二进制速度文件第 129 页 共 261 页二、OMEGA 三维资料处理流程举例图 127 将tx

37、t 速度文件转换成二进制速度文件第 130 页 共 261 页2.9 速度内插（3）SPATIAL_PARM_INTERP图 128速度内插第 131 页 共 261 页二、OMEGA 三维资料处理流程举例图 129 速度内插第 132 页 共 261 页二、OMEGA 三维资料处理流程举例2.10 第一次剩余静校正2.10.1叠加生成剩余静校正分析所需数据（1）（1）（2）（2）（1）（2）（1）（2）图 2-130 叠加生成剩余静校正分析数据流程第 133 页 共 261 页2.10 第一次剩余静校正2.10.1 生成剩余静校正分析所需数据（1）INPUT图 2-131 输入 CMP 道集

38、数据第 134 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-132输入 CMP 道集数据第 135 页 共 261 页2.10 第一次剩余静校正2.10.1 生成剩余静校正分析所需数据（2）TRACE_SELECT（1）图 2-133 道选择第 136 页 共 261 页二、OMEGA 三维资料处理流程举例（3）NMO图 2-134 动校正第 137 页 共 261 页2.10 第一次剩余静校正2.10.1 生成剩余静校正分析所需数据（4）MUTE图 2-135 定义切除库第 138 页 共 261 页二、OMEGA 三维资料处理流程举例（5）RMS_GAIN（1）图 2-1

39、36均方根增益第 139 页 共 261 页2.10 第一次剩余静校正2.10.1 生成剩余静校正分析所需数据（6）STACK图 2-137 叠加第 140 页 共 261 页二、OMEGA 三维资料处理流程举例（7）DATUM_CORRECTION（1）图 2-138基准面校正第 141 页 共 261 页2.10 第一次剩余静校正2.10.1 生成剩余静校正分析所需数据（8）BPFILTER图 2-139 零相位滤波第 142 页 共 261 页二、OMEGA 三维资料处理流程举例（9）RMS_GAIN（2）图 2-140均方根增益第 143 页 共 261 页2.10 第一次剩余静校正2

40、.10.1 生成剩余静校正分析所需数据（10）OUTPUT（1）图 2-141生成供剩余静校正分析用的叠加数据第 144 页 共 261 页二、OMEGA 三维资料处理流程举例（11）TRACE_SELECT（2）图 2-142 选取需要打印的 CMP 线剖面第 145 页 共 261 页2.10 第一次剩余静校正2.10.1 生成剩余静校正分析所需数据（12）SEISPLOT图 2-143 生成质量图片第 146 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-144生成质量图片第 147 页 共 261 页2.10 第一次剩余静校正2.10.1 生成剩余静校正分析所需数据图

41、2-145 生成质量图片第 148 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-146生成质量图片第 149 页 共 261 页2.10 第一次剩余静校正2.10.1 生成剩余静校正分析所需数据图 2-147 生成质量图片第 150 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-148生成质量图片第 151 页 共 261 页2.10 第一次剩余静校正2.10.1 生成剩余静校正分析所需数据图 2-149 生成质量图片第 152 页 共 261 页二、OMEGA 三维资料处理流程举例（13）DATUM_CORRECTION（2）图 2-150 基准面校正第 1

42、53 页 共 261 页2.10 第一次剩余静校正2.10.1 生成剩余静校正分析所需数据（14）OUTPUT（2）图 2-151输出供剩余静校正分析使用的 NMO 数据第 154 页 共 261 页二、OMEGA 三维资料处理流程举例2.10.2 地表一致性剩余静校正分析方法原理一个道的反射时间可由下式表示：T = S + D + G + MX 2ijijkkij其中，Si 为点静校正量；Dj 为检波点静校正量；Gk 为构造项；Xij 为检距；Mk 为剩余动校正因子；MkXij2 为剩余动校正项。实现步骤（1） 用动校正的 CMP 道集和模型道做互相关，求出 Tij；（2） 根据 Tij，按

43、上面公式用高斯-塞得尔迭代分解法求得 Si、Dj。（1）（2）（1）（2）图 2-152 用内部模型计算剩余静校正量的流程第 155 页 共 261 页2.10 第一次剩余静校正2.10.2 地表一致性剩余静校正分析（1）INPUT（1）图 2-153 输入 2.10.1 步骤所得的 NMO 数据第 156 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-154 输入 2.10.1 步骤所得的 NMO 数据第 157 页 共 261 页2.10 第一次剩余静校正2.10.2 地表一致性剩余静校正分析（2）BPFILTER（1）图 2-155 零相位滤波第 158 页 共 261

44、页二、OMEGA 三维资料处理流程举例（3）INPUT（2）图 2-156输入 2.10.1 步骤所得的叠加数据作为内部模型道第 159 页 共 261 页2.10 第一次剩余静校正2.10.2 地表一致性剩余静校正分析图 2-157输入 2.10.1 步骤所得的叠加数据作为内部模型道第 160 页 共 261 页二、OMEGA 三维资料处理流程举例（4）BPFILTER（2）图 2-158零相位滤波第 161 页 共 261 页2.10 第一次剩余静校正2.10.2 地表一致性剩余静校正分析（5）XPERT图 2-159 用动校正的 CMP 道集和模型道做互相关第 162 页 共 261 页

45、二、OMEGA 三维资料处理流程举例图 2-160 用动校正的 CMP 道集和模型道做互相关第 163 页 共 261 页2.10 第一次剩余静校正2.10.2 地表一致性剩余静校正分析图 2-161 用动校正的 CMP 道集和模型道做互相关第 164 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-162 用动校正的 CMP 道集和模型道做互相关第 165 页 共 261 页2.10 第一次剩余静校正2.10.2 地表一致性剩余静校正分析关键参数选择（1）相关时窗因为相关运算要在该时窗内进行，所以该参数的选取好坏直接关系到剩余静校正的效果。建议时窗包含浅-中层反射。如果时窗选取

46、得太浅，浅层反射受速度影响较大，当速度时会有残留的剩余动校正量；如果时窗选取得太深，深层反射对静校正又不太敏感，因此，太浅或太深的相关时窗都会影响本模块的使用效果。（2）模型道集扩展该参数指定参与建立模型道的道集个数。一般情况下，低信噪比资料可适当大些，高信噪比资料可适当小些。（3）起始 CMP 号该参数指示互相关运算是从哪个道集开始的。一般情况下，不一定要求从第一个道集开始，应选取资料品质较好的剖面段开始互相关运算。程序以该 CMP为起点，自前、向后运算。（4）最大时移量为了避免发生窜相位现象，该参数对互相关得到的最大静校正量给予限制。一般情况下，该值不能超过主视周期的一半。第 166 页

47、共 261 页二、OMEGA 三维资料处理流程举例（6）REFLECTION_MISER图 2-163求解反射波剩余静校正量第 167 页 共 261 页2.10 第一次剩余静校正2.10.2 地表一致性剩余静校正分析图 2-164 求解反射波剩余静校正量第 168 页 共 261 页二、OMEGA 三维资料处理流程举例（7）OUTPUT图 2-165输出反射波剩余静校正量第 169 页 共 261 页二、OMEGA 三维资料处理流程举例2.10.3应用剩余静校正量（1）（2）（1）（2）图 2-166应用剩余静校正量流程第 170 页 共 261 页2.10 第一次剩余静校正2.10.3 应

48、用剩余静校正量（1）INPUT（1）图 2-167 输入 CMP 道集第 171 页 共 261 页二、OMEGA 三维资料处理流程举例图 2-168输入 CMP 道集第 172 页 共 261 页2.10 第一次剩余静校正2.10.3 应用剩余静校正量（2）TRACE_SELECT（1）图 2-169 道分选。在上面所选三条 CMP 线前后再分别选三条 CMP 线用于二次速度分析。第 173 页 共 261 页二、OMEGA 三维资料处理流程举例（3）INPUT（2）图 2-170 输入第一次剩余静校正量。该剩余静校正量可由 ATTRIB DISP 加载显示，以实现剩余静校正质量。第 174 页 共 261 页2.10 第一次剩余静校正2.10.3 应用剩余静校正量图 2-171 输入第一次剩余静校正量第 175 页 共 261 页二、OMEGA 三维资料处理流程举例（4）STATICS_UPDATE图 2-172应用剩余静校正量第 176 页 共 261 页2.10 第一次剩余静校正2.10.3 应用剩余静校正量（5）TRACE_SELECT（2）图 2-173分选用于质量的 CMP 线第 177 页 共 261 页二、OMEGA 三维资料处理流