超级CGG培训资料55.pdf-淘文阁

资源描述

《超级CGG培训资料55.pdf》由会员分享，可在线阅读，更多相关《超级CGG培训资料55.pdf（70页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、 CGGLEVEL 5000 培训 1 CGG简介 CGG 是目前国际较为通用的一款地震数据处理软件，与其他地震数据处理软件相同，它将所有处理功能模块化，通过多个模块的搭配使用来实现某一种处理功能，一系列模块的搭配称之为一个 JOB，一系列 JOB 的组合来完成对整个地震资料的处理。相对 PROMAX，CGG 具有更高的交互水平，尤其是新版本的 CGG，增加了作业流程图和作业文本的同步互译功能，使得操作更加便捷。CGG 具有更友好更强大的 QC 功能，它拥有单独的 QC 工具，在此工具下能实现各种形式的查看、相关数据运算、滤波、谱分析等功能，使得 QC 变得更加灵活实用。另外，CGG 将被处理

2、数据按照用途的差异进行了详细规划，具有更清晰的数据管理机制，便于数据的调用和管理。但是，相对 PROMAX 来说，CGG 却又显得不够干净利索，对于 PROMAX 的使用者来说，初接触 CGG 会给人一种臃肿的、拖泥带水的感觉，尤其是文本作业的形式，让人有一种在写程序的感觉，不过随着学习的深入，我们会逐渐发现其便捷所在、优势所在。在此特别注意的是，5000 版本的 CGG 增添了作业流程图和作业文本的同步互译功能，这一功能存在互译不稳定性，对于一些模块及模块之间的数据链接不能正确翻译为作业文本，此时我们需要以文本为主，严格按照文本格式进行修正，保证文本准确即可。但这一缺陷对 CGG 的交互水平

3、影响不大，同步互译功能自有其妙用。以上概述，在下面的讲解中会得到详细的体现。2 数据文件分类 2.1 辅助数据文件 XPS，是 CGG 内部专门为辅助数据而设置的文件存储格式，如切除文件、拾取水底资料、速度文件等地震属性的非地震数据主体文件库，分类细致，具有独立的 QC 窗口，具备列表和成图功能，便于查看和分析。很多模块中具有调用 xps 文件的接口，使用方便。XPS文件由PROJECT（预存储工区名）、ID 和 VERSION组成，要求每一个名称长度不能超过 15 个字节。2.2 作业文本文件 CGG 中，作业文件分为 3 种，最新版本的.job文件、老版本的.xjj文件和.gsl的作业文本

4、文件，存储于工区目录下的 JOBS 目录下。2.3 库文本文件 CGG 对库文件具有很详细的划分，每一个库模块对应一个库文件类型，例如，.lmu(切除库)、.lvi（速度库）、.lfd（水深库）、.lcn（滤波库）、.ldb（输出库）等等。有很多库模块在调用库文件时具有双重选择，既可以选 XPS 文件也可以选择文本库文件。2.4 地震数据文件 CGG 中，地震数据根据输入输出的不同分为多种，.DAT文件同时用于输入输出，用前缀P和I来区分，分别为外部文件和内部文件，具有非常严格的格式要求，格式如下：外部：内部：外（内）前缀编号工区号.DAT或外（内）编号工区号.DAT，前缀和编号又有严格规定，

5、或者位前缀个数字，或者直接用个数字。外部的 DAT 数据可由模块 SEGIN 读入，外部的.sgy数据需改名整理成标准的外部 DAT数据才能够被 CGG 识别读入。内部的 DAT 数据可由 OUTBD 输出或由 INPTR 读入。DAT数据一般只出现在外部读入和最后输出。.cst文件是内部的临时地震数据体，一般用于 QC，主要是测试处理效果等，暂存在临时数据存储的目录下，由模块 WUNET 和 RUNET 输出和读入，称其为网络文件，不作为地震数据的主体。SDS文件，是真正的地震数据文件，类似于 XPS 文件无格式后缀，CGG 为 SDS 单独在工区中开辟一个存储空间，存储路径由PROJECT

6、（当前工区名）、ID 和 VERSION组成。是在整个流程中处理和输出的对象，包含全部的地震数据。由模块 SETRA 和 GETRA输出和读入。2.5 导航定位文件.P190文件是野外采集的导航定位文件，属于外部文件，包含野外采集各炮点、检波点的坐标及沉放信息以及声纳采集的海底信息，是地震数据的一个重要部分，有专人专门制作，由模块 WILOC 将其读入并转化为内部格式.DAT文件，CGG 对其文件的命名和路径没有特别要求。2.6 其他文件.list是作业执行列表，从该文件中可以查看作业执行进程、错误原因和执行结果信息等，是一个非常重要的文件，每一个作业执行后都会生成一个.list文件。.log

7、系统内部的操作日志，没什么用。.velcom速度谱文件，由模块 VESPA 生成，用于速度谱的拾取。.lgr三维 CDP 网格文件，内部含有原点坐标、CDP 间隔、CDP 线间隔、CDP 方向方位角及 CDP 线方向方位角等信息，用于 3D 建观。2.7 数据存储规划 CGG 的根目录下设有常用的 proj、data1、stage 等子目录，proj 里面是我们所建立的所有工区，每个工区下又有常用的 DATA、LIBRARY、LIST 等子目录，data1、stage、DATA和 LIBRARY 是各种数据的存储位置，一般 data1 作为工区外部路径用于存储外部数据，如P.DAT 数据（SG

8、Y 数据）和 P1-90 数据；stage 一般用于存储内部的 DAT 数据，如 P1-90 的转换数据；而 DATA 中一般临时存储 CST 数据文件，占用工区空间，一般用于 QC；LIBRARY中只存储库文件，如 LMU（切除库），LVI（速度库）等，VELCOM（速度谱文件）。LIST中存储作业执行列表。当然，除了 LIST 之外的这些存储路径都是可以更改的，在此只是提供一种比较合理的存储规划。2.8 CGG内部格式道头介绍 CGG 中地震数据常用道头如下所示（MOT 或 WORD 均为 CGG 可识别的道头字前缀）MOT2：野外炮号 MOT4：CDP 号（3D 数据时为由导航数据确定的

9、标准 CDP 号，2DGEOUP 建立的 CDP号放置在 MOT52 中）MOT9：采样率 MOT10：记录长度 MOT11：有效/废道标识，7 为有效道 MOT12：逻辑线号 MOT15：道间距 MOT16：延迟时 MOT17：道号 MOT18：检波点号 MOT19：三维线号 MOT20：偏移距 MOT22：静校正量 MOT26：缆号 MOT30：MOT33：枪号 MOT34：HABIN 之后属于同一个偏移距的道标识 MOT36：HANBIN 之后 EVERY 计算的 RMS 值 MOT37：MOT43：CDP 点的 X 坐标 MOT44：CDP 点的 Y 坐标 MOT52：通常为空值，习惯

10、将 3D 数据的 2DCDP 号赋给它 MOT60：检波点 X 坐标 MOT61：检波点 Y 坐标 MOT62：炮点 X 坐标 MOT63：炮点 Y 坐标 MOT64：3 界面介绍及常用工具不同于 PROMAX，CGG 为用户开发了很多用以实现特殊功能的工具，如速度分析、批量处理编辑、QC 等，使得处理流程更加清晰，操作更加便捷。3.1 CGG的启动在 terminal 中键入 kereon 然后回车即可启动 CGG，这时会出现一个 kereon master 界面，新建或双击某个已存在的工区名即可进入，或者左键选择一个工区名，然后在 File 下拉菜单中点击 open kereon 亦可

11、进入。工区界面中一般主要包含以下：1.FastQC 2.MapQC 3.Geodesk2 作业管理程序（批处理）4.Jseg2FI 对野外地震数据进行分析 5.Sditr 折射波静校正（海上处理时只在 OBC 做）6.TeamView 显示工具（质量控制，包含部分处理功能，如滤波、自相关等）7.Xps 辅助地震数据库（如切除、滤波、速度等文件都可以存放在这里）8.Chronovista 速度分析 9.xjob 交互作业（老版本）10.Jxjob 交互作业（新版本 5000）11.Xenu 队列管理（选择发送作业的节点）12.Dvimon 作业监控（监控作业运行情况）13.Isub 发送作业（也

12、就是执行作业）14.Xopreg 资源管理器 15.Xdoc 说明书可以右键添加（Application）常用的，例如：GRILLTH 建立网格下面简单介绍一下 TeamView 和 Chronovista。3.2 TeamView 显示工具（质量控制）TeamView 主要是用于质量控制（QC），类似于 Promax 中的 Trace Display，但比 Trace Display 的功能要强大许多，不仅可以显示各种分选形式的地震数据，还可以对地震数据做一些简单的处理，如：加、减、乘、除、自相关等运算。1.地震数据的显示 File-creat project 创建临时工区，当对这个工区

13、做出改变后，退出时才会提示是否保存。点击工区名左边的圆圈，出现下拉菜单，左键点击 seismic，右键选择 load（一定要先左键点击 seismic 才行）加载地震数据，由于 TeamView 主要是用于质量控制（QC），所以加载的数据一般比较小，一般是.cst 的文件（WUNET 输出的数据），比较大的数据（例如最终的叠加剖面）一般是.sds 的文件。数据加载完后，把数据用左键拖到右边的显示窗口即可显示。右键点击剖面可以进行相关显示参数的修改，以达到最佳的显示效果。在进行 AGC 显示时，把 Samples 改为 Seconds，属性显示还要勾上 Apply AGC。对数据进行分选显示，

14、左键选中所加载的地震数据，右键选择分选 DATA BOSRT，例如按炮集、CDP 道集显示。多剖面显示时，如果想把对某一剖面的显示参数应用到其它剖面，可以直接点击传播图标即可。如果对两个剖面作对比显示，例如对去多次波前后的叠加剖面做对比显示，先左键选中两个地震数据，然后拖到显示窗口中或右键 Multiple Display，点击左上角的 View，再左键选中两个地震数据，右键选择最大化显示，修改显示参数，得到最佳显示效果后然后点击传播图标，再左键选中两个地震数据，右键选择 movie（动画显示）即可。2.地震数据的简单处理在 TeamView 中可以对地震数据做一些简单的处理，如：加、减、乘

15、、除、自相关等运算。自相关主要用于检测反褶积中 Gap（间隙值）是否合适，标准是旁瓣越小越好；减法可以检测处理前后两个数据体的差异。在 TeamView 中还可以对线性噪声做速度的拾取（去线性噪声的时候需要用到），对地震数据做各种切除，拾取水底库等，也可以对地震数据做频谱分析，这里需要注意的是，对滤波后的数据做频谱分析得到的结果与对没有滤波的数据做频谱分析得到的结果是一样的，因为这里的滤波只是一个显示，并没有改变原始的数据。3.3 Chronovista 速度分析在做速度分析之前，必须要先由 VSEPA 得到一个.velcom 文件，类似 Promax 里速度预计算得到文件，这里主要介绍一下

16、 CGG 进行速度分析的步骤：Project-Creat（产生一个临时目录，目录的第一个字符必须是字母，路径必须是 proj/当前工区名/GEM，系统只认这一目录）-点 Geovel（双击）-Warning-OK（在这里定义一个观测系统，对于二维数据 Topography 是虚拟的，设置 Geometry 时必须将所有的线号和 CDP 号包含进来）-OK-Geovel-Datamanage-velocity-Analysis-Import From（Velcom file）-List（找到.velcom 文件）-Load-关闭当前窗口-左键选中加载后的.velcom 文件，再右键-view-进

17、入速度普拾取界面（每次拾取速度前将右上角所有工具全选中）-Velocity-new-建立一个速度文件（这一步非常的重要，否则无法拾取速度）速度拾取完后先保存，然后需要导出，右键你拾取的速度文件名 Export-弹出一个对话框-即可按文本文件输出也可按.xps 文件输出，一般选择按.xps 文件输出，得到的速度文件一般还需要人工修改，主要包括扩边（扩展到第一个 CDP 号和最后一个 CDP 号）、设置初速度（一般为海水的速度）和末速度（最大速度）、允许速度反转（加上 TN，不加则默认为反转 50%的速度）。3.4 Geodesk2 作业管理程序该工具主要用到的是其批量处理的功能，通过表格控制作

18、业的批量执行，该部分在三维资料处理中具有重要意义。详细内容待补充。4 常用模块功能简介 SEGIN：外部地震数据的输入。只能输入 SEGY 或 SEGD 格式的地震数据，可以同时输入多个外部地震数据。SEGIN 对输入进来的地震数据的命名有严格的要求，必须按照 CGG 的命名规则对输入进来的地震数据重新命名，命名规则：P前缀编号工区号+.DAT 前缀和编号又有严格规定，或者位前缀个数字，或者不要前缀直接用个数字。与 SEGIN 相对应的还有两个库：格式库（LIBRI FI）和输入库（LIBRI SI）。由于 SEGIN 读入的地震数据没有给出具体的路径，只给出了盘号和编号，所以还需要设置一个输

19、入环境变量（即该文件的存放路径），例：*GURPA setenv gvr_stage /data1/data_BH 类似的输入输出若只给出了盘号和编号的，都需要设置一个环境变量。如果一个处理流程中既要设置一个输入环境变量也要设置一个输出环境变量，且输入输出路径不一样，则要用 i（输入）和 o（输出）来做区分,例：*GURPA setenv gvr_stagei /data1/data_BH*GURPA setenv gvr_stageo /data1/stage 如果输入输出路径一样，就不用区分了，且设置其中一个即可。INPTR：一般读取 OUTBD 输出的文件（I前缀编号工区号+.DAT）或

20、 WILOC 后得到的导航数据（如 IWI10013DSRME8.DAT），与之对应的有一个 TR 库。也需要设置一个环境变量。GETRA：一般读取.sds 的地震数据，需要给出 Project（工区名）+Identity(目录)+Vesion（文件名），如果不清楚该文件的 Project+Identity+Vesion，可以去 SDS 库中找到该文件，然后双击会出现该文件的属性，从其属性中可以查找到。可以对数据进行分选读取，例：WHERE SERIES MOT2=(1001-1756,I2)，按炮号读取，从 1001 炮到 1756 炮每隔两炮读一次。RUNET：读取网路文件（.cst 文件

21、），一般用于 QC。OUTBD：输出.DAT 格式的地震数据或导航文件，与之对应的有一个 BD 库，输出的文件名称按照以下规则自动命名：I前缀编号工区号+.DAT 前缀和编号又有严格规定，或者位前缀个数字，或者不要前缀直接用个数字。也需要设置一个环境变量。一般和 SEGIN 或 WILOC 搭配使用。SETRA：输出.sds 的地震数据，一般和 GETRA 搭配使用，输出比较大的数据。由于这个版本不能正确翻译 SETRA 模块，所以每次还得人工修改文本，修改如下（标准格式）：*SETRA AA REPLACE DATA INTO 3DSRME8/Leo/raw_combine USING CO

22、MPRESSION=F2 WUNET：输出网络文件（.cst 文件），一般用于 QC，输出比较小的临时文件。由于这个版本不能正确翻译 WUNET 模块，所以每次还得人工修改文本，在前两个 FILE 的最后加上一个“+”号即可，这一点必须注意。修改如下：*WUNET AA FILE=local:+FILE=/data1/Leo_DATA/+FILE=RAW1020.cst,STEP4800,RUNET 也需要这样修改！SEGOU：地震数据的输出。输出 SEGY 或 SEGD 格式的地震数据。也需要设置一个环境变量。WILOC：P1/90 导航数据的输入。GEOUP：定义二维观测系统，或者对 3D

23、数据体置 2D 的 CDP 号。EGRID：定义三维观测系统。需要定义一个网格库（LIBRI GR）LISTE：道列表。输出一个或多个地震道的值，如炮号。QCXPS：提取各种地震属性，如水深库。RSAMP：重采样。若数据量较大，2ms 采样的数据可先重采样为 3ms 的数据在进行后续处理，以提高工作效率。因为一般海上数据的有效信号的频率在 3ms 的 Niquist 频率范围之内，所以重采样不会影响最终处理结果的精度。注意输出最终的结果时必须再重采样为 2ms 的数据。SCALE：道缩放，按比例显示。BSORT：分选。MUTES：切除。RAMUR：去线性噪声（NA）或着去中远距离多次波（MA

24、）。FDNAT：分频去噪。FXNAT：F-X 域随机噪声衰减。SPARN：分频去噪。SPASM：去野值。FILTR：滤波，包括各种滤波，如低切滤波、高通滤波、带通滤波、去陷波等。SELTR：分选。MODET：修改道头值。例：*MOT18=MOT17（把道头 17 的值赋给道头 18）*MOT26=1(把缆号设置为 1)*MOT4=3127，MOT2（当前炮号），MIN（减），6，MULT（乘），MOT17（通道号），PLUS（加）意思是（3127-MOT2）*6+MOT17 MNGTY：Y 位控制，在进行多道处理时前必须加上 MNGTY，一般还可以在 MNGTY 前加一个 SELTR 先把数据

25、分选一下。如后面的处理是在 CDP 道集上，则先用 SELTR把数据按号分选，然后再用按做控制。FXINT：道内插值（在 F-X 域）。一般在做某些处理时先插值再处理，这样做可以防止空间假频。REFOR：振幅补偿。一般是为了计算用的，例如计算 RMS（均方根值）。HISTA：静校正（如水深校正）或者用于提高各种变换精度的运算（如在 Tau_P 变换前可以先做一个 HISTA，提高 Tau_P 变换的精度）。LINTP：线性 Tau-P 变换。TRITA：反褶积。一般在 Tau-P 域做反褶积，效果要好于在时间域的。主要去层间短周期多次波。SMMOD：建立多次波模型，用于 SRME（主要是去近偏

26、移距多次波）。LASUB：道头一致性检查或做减法处理（原始数据减去模型数据），用于 SRME。DUPLI：数据流的复制。把多个流程的数据输入到另一个流程中。DYNQU：增益。类似 AGC。RESTR：TRVEL：产生速度道。MODVI：速度场的平滑。两个方向都得平滑，一次按线方向（MOT19,MOT4），一次按 CDP方向(MOT4,MOT19)。而 MODVI 每次这能按一个方向对速度场进行平滑，所以得分两次对速度场做平滑，一般有两种方法实现两个方向的平滑：1，利用 Tabedit 打开一个速度场（3D.lvi）文件，缺省情况下，生成表格的第一列是线号（MOT19），第二列是 CDP（）号，

27、第三列是时间（），第四列是速度（）。首先按照缺省的情况（MOT19,MOT4）用 MODVI 对该速度场平滑一次，即按线方向平滑了一次。然后把平滑后的速度场用 Tabedit 打开，利用以下步骤把平滑后的速度场按 CDP 方向(MOT4,MOT19)重排列：A，选中第一列（线号）-Edit-Insert Filed-在第一列前家一列；B，Operation-Formula-$1=$3-把第三列数据（CDP 号）复制到第一列；C，Edit-delete-删除第三列（CDP 号）；D，Export-导出，得到一个按 CDP 号排列的速度场文件。再用 MODVI 对按 CDP 号排列的速度场文件进行

28、平滑，至此速度场实现了两个方向的平滑。2，利用 TRVEL，步骤如下：LIBRI VI-TRVEL-LT 速度道（MOT19,MOT4）-MODVI-BSORT（MOT4,MOT19）-MODVI TIKIM：偏移。HFNMO：动校正或反动校正。STAPA：保持振幅叠加（有振幅加权）。STACK：水平叠加（无振幅加权）。IFTHN：条件判断。QCAMP：计算均方根值。GAINA：计算增益值。GAINX：地表一致性增益补偿。HABIN：面元均化。某一面元内的某共偏移距范围（例如：0-275）内的道数在理想的情况下只有一道或 0 道，但是实际由于拖缆摆动的原因，共偏移距范围内的道数可能有多个，通过

29、 HABIN 可以把多余的道删掉，只保留离面元中心最近的那一道。这就保证了某一面元内的某共偏移距范围内的道数只有一道或 0 道。REG3D：数据规则化。由于 HABIN 之后某一面元内的某共偏移距范围内的道数可能是 0 道，REG3D 通过插值，是得每一个面元内的每个共偏移距范围内的道数都是一道，并把每一道都移到面元中心上去，这样就达到了覆盖次数均一化的效果，避免偏移噪声。LAYES：HDPIC：高精度速度场的拾取。FXYNA：F-XY 域随机噪声衰减（3D）SPARC：2D SCAC 地表一致性振幅校正 PRF3D HABIL MULIN RSAMP SDICO：球面扩散校正。3DSRME

30、EXPEN：剩余振幅补偿。SYMUL：MOMUL：TRASY：TDNFK：5 作业格式介绍（交互格式和文本格式）CGG 作业流程中所调用的模块分为两类，即库模块和功能模块，当流程中读入或调用数据时需要通过调用功能模块并连接相应库模块来实现。库模块一般只用于标识所调用数据的路径。CGG 作业流程根据数据体的形态分为诸多循环，流程中每次分选、域转换、插值、重采样等改编数据体形态的操作之后都要重启一个循环，针对同一形态数据体的所有处理一般情况下均在同一个循环中。不同循环间的数据接口一般用标识循环代号的数字，循环内的数据接口一般用大写字母或+=来标识，当用+=时，+一般标识输出，=一般标识输入，但因其

31、单一性只能完成上下一对一的数据传递。在编写流程时，我们往往先将所需要的模块按顺序拖入交互界面形成一个完整的流程图，这样在文本中便会自动正确划分循环（但也会有特殊情况的出现，如：做 LASUB 时会发生翻译紊乱），自动正确设置数据流接口，循环和数据流接口建立好之后再设置各个模块的参数，在此需要特别注意的是，交互模式填写的参数有时不能正确翻译为文本，反之亦然，这就需要手动修改文本，因为我们的作业是以文本为主的，交互流程图仅仅是为了增强软件的交互能力，因此交互流程的紊乱我们可以不管，但文本一定要正确，这也是 CGG 流程文本互译时不稳定性的所在，但只要文本正确就不会影响到我们的处理效果。作业文本是

32、DOS 格式，*打头表示注释掉后续行，文本框的下拉菜单中有批量注释掉的功能，在此要特别注意的是，该注释并不同于 PROMAX 中的注释，PROMAX 中做完注释后会自动重新规划数据接口，CGG 文本中做完注释不会自动重新规划数据接口，需要手动修改接口。而在交互流程图中又没有注释掉某模块的功能，只能进行删除或切除操作，但这两中操作是具有数据接口自动重设功能的。因此需要注释的时候我们一般还是在文本中操作。一个示例作业文本的格式如下：*LEVEL 5000*环境变量设置，CGG 版本号*PATCH S_HAL5000_1322*PATCH P_TIKIM5000_7203*环境变量设置，补丁*GVR

33、PA setenv HAL_OMPI_DISABLE ALL*环境变量设置，启用或关闭并行计算*GVRPA setenv gvr_local/data1/stage*环境变量设置，调用或存储路径，固定格式*XPSID ID=3dmarine1 PROJECT=2011bhh*格式符后跟模块，XPSID用于启动 XPS，标识 XPS 数据的 PROJECT 和 ID *LIBRI VI 01 TN DBVI=PSTMVEL*库文件，速度文件*LIBRI SH 01 DBSH=SH*库文件，各向异性文件*PROCESSING GRID*LIBRI GR 01 XOR602024.6,YOR4242

34、958.3,*库文件，网格文件 AZX2.0612338,AZYMS,DCDP 12.5,DLINE 25.0*ANTI-ALIAS GRID*LIBRI GR 02 XOR602024.6,YOR4242958.3,*库文件，网格文件 AZX2.0612338,AZYMS,DCDP 18.75,DLINE 37.5*DLOOP 1 *DLOOP一个循环的开始标识，也用BOUCL，1循环编号*TIKIM MI BB TRAFILE=TIKIM_IN,NUMTRA1,IDTRA=OFF#1#,IMAGE=(FCDP7,LCDP3027,ICDP1,NUMRUN1,FLINE847,LLINE90

35、5,ILINE1,DT4),BIMG02,WORDX4,NOALIAS=(FMAX120,LGR2)DIPLIM70,APERCDP4500,APERLIN4500 LVI1,LSH01,*RAYTRACE,*注释掉的参数 NOREG,SAVDISTR,MACHINE=PC,DISTRIB,REWRITE,N0,*TIKM偏移模块名，MI模块功能选件，此处表示处理和输出，BB数据流*标识，后面的一堆是模块的参数设置代码。*ENDLP*该循环的结束标识，也用FINBO*DLOOP 2*SETRA BB REPLACE DATA INTO 2011bhh/3dmarine/PSTM_OFF#1#*

36、SETRA输出 SDS 数据，BB接TIKIM的输出数据*ENDLP *PROCS X(YB1)*PROCS流程的结束标识，自动生成，X(YB1)循环控制符，自动生成 6 海上 2D 资料常规处理流程本节介绍一个海上 2D 资料常规处理流程，需要特别注意的是，整个流程中的每个处理环节都需要进行严格的 QC 进行处理质量监测。流程如下：1、SEGIN，SEGIN数据读取和海平面校正。2、GEOMETRY，GEOUP（无 P1-90 文件时用来建规则二维观测系统），WILOC+EGRID（有 P1-90 文件时建精确的二维观测系统）。3、DENOISE，预处理去噪，RAMUR(NA)去线性噪音，

37、正反两个方向都做，SPARM去野值模块，SPARN和FDNAT去规则噪音，分频去噪。4、SRME，AMMOD+LASUB，两步除近偏移距海底多次波。5、DECON，TRINA为单道反褶积，MCDEC为多道反褶积，本流程中做 T-P 预测反褶积消除层间多次波。6、VESPA，VESPA速度谱计算，第一次拾取速度谱。7、HFNMO+RAMUR(MA)，CDP 道集去远偏移距海底多次波，因为 RAMUR 在动校正后的数据上做。8、查看效果，还有多次波残留时重新做速度分析，动校，去多次，直至满意为止。9、2DSCAC，在炮集上做地表一致性处理，主要消除炮间和道间的差异，QCAMP+GAINA+GAIN

38、X，路上资料用 NMBSO 模块。10、REFOR，振幅一致性处理，消除单道内的振幅衰减，SDICO，球面扩散校正，消除球面上的振幅衰减。11、PSTM，偏移处理，计算偏移速度场。12、VESPA+PSTM，第一次偏移后重新拾取速度进行第二次偏移速度场计算。13、STACK。14、叠后处理，FXNET+FILTER。15、AGC，LFK 去除叠后浅层噪音，如海底抖动。第 3 步中，如果资料道数有限，无法辨别线性噪音则省去线噪去除，如果叠后资料中存在线性噪音，可以叠后去线性。第 7 步，同样当资料道数有限时就无远偏移距可言了，如第八节中所讲到的实例中便没用到。在此特别介绍本流程中去海底多次波的思

39、路如下：1、利用去噪之后的资料 A 建立一个海底多次波模型 B。2、用原始资料 A 减去海底多次波模型 B 得到 C，检查 C 效果，质量不好需修改模型。3、拾取 2 过程获得数据的速度资料 V1。4、利用 V1 对资料 C 进行动校正处理，得到数据 D，使有效波校平，多次波明显被识别。5、去除 D 中的中远道偏移距海底多次波，得到 E。6、利用 V1 对 E 数据做反动校，后进行速度分析拾取速度 V2。7、利用 V2 对 E 进行动校正处理，检查效果。7 海上 3D 资料常规处理流程海上 3D 处理的流程如下：1、SEGIN，SEGIN地震数据导入。2、WILOC，导航文件导入，QCXPS

40、，提取水底库。根据 WILOC 的执行列表建立工区的三维 CDP 坐标网格，存为.lgr 格式的网格文件。3、EGRID，定义三维网格，CGG 习惯将网格原点定义在实际工区的左下角。建立LGR 文件包括参数 xor、yor、azx、azy、dcdp 和 dline。4、DENOISE，预处理去噪，高通滤波，RAMUR(NA)去线性噪音，正反两个方向都做），SPARM去野值模块，SPARN和FDNAT去规则噪音，分频去噪。5、若有远场子波的话可以做零相位化处理。6、初至切除。7、2DSRME，抽取单源单缆做近偏移距海底多次波去除，该过程需要在 CDP 道集上做，所用的 CDP 号可以直接用MOD

41、ET来计算，也可以用GEOUP来划分，因为三维资料每条缆的最小偏移距是不等的，因此推荐用 MODET 更便捷。8、DECON，T-P 预测反褶积消除层间多次波。9、VESPA，获得速度资料 V1。10、RAMUR，出去中远偏移距海底多次波。11、2D SCAC，炮集上做地表一致性处理，该处理要对每一条航行线去做，每条航行线上做一个GIANA，应用到对应的航行线上，一条航行线包含所有的 CDP 轨迹线，二源六缆则包含 12 条线的数据。12、HABIN，共偏移距组划分，并筛选 CDP 覆盖道，处理结果是将数据划分为 N组，每组内的 CDP 点上最多存在一道覆盖。该过程将所有数据按照偏移距划分为

42、N组，每组均包含位于一定偏移距间隔范围内的数据，原则上一个偏移距道集内每一道数据均对应一个 CDP，即共偏移距道集所能够覆盖到的 CDP 原则上仅能被覆盖一次，因此我们通常将道集内包含多次覆盖的 CDP 点作剔道处理，一般保留距 CDP 中心最近的那一道，剔道完毕对道集中的每一道求取均方根给已建立的共偏移距 CDP 网格赋值，存在覆盖的便存在值，覆盖不到的便显示空缺（要注意的是所有共偏道集的 CDP网格是统一的，即按照最小坐标值和方位角确定的那个网格）。假设资料的满覆盖次数为 N 次，则理论上应该分为 N 个共偏移距面（组），组距原则上应该为两倍的炮检距，但实际上野外采集时缆是漂动不定的，因此

43、要按实际坐标的观测系统求取，公式为：NNRMXRMX，其中 XRM、NRM 分别为实际的最大、最小偏移距。理论上第一组偏移距范围应为0，NRM，但因为每一条缆的 NRM 均不同，因此应该选择其中最大的 NRM，基本原则就是第一个共偏移距组要将所有缆的首道包含在内，确定了第一个共偏移距面后，后续的共偏面按照X划分即可，参考图 7-1 理解。图 7-1 共偏移距面的划分图示 13、REG3D，3D 数据规则化，在 HABIN 的基础上进行规则化处理，该处理结果有二，一是将每个共偏移距面上的 CDP 覆盖道进行规则移位，移动到每个 CDP 面元的中央；二是将内部零覆盖的 CDP 面元上进行插值，使其

44、具有一道的覆盖次数；这样就保证了每个 CDP 上具有相同的覆盖次数，并使得数据变得整齐。便于后续偏移速度场的计算。14、DENOISE，SPASM、SPARN、FDNAT、其中FDNAT为二维去噪，在三维数据体上应用时应该分别在 CDP 号方向和逻辑线方向两个方向上做，最终获得速度资料 V2。15、PSTM，将 V2 进行两个方向的平滑得到 V3、SH1，输出 V3 的速度线。16、TIKIM(IN)(MI，OU)，计算偏移速度场获得 V4、SH2，进行第二次偏移速度场计算，即 Full PSTM（V4，SH2），输出 CDP 道集之后做 HDPIC（高密度速度拾取）得到 V5、SH3，当然拾

45、取前可以视情况做线性噪音的去除。17、HFNMO，STAPA。18、叠后处理，AGC、FILTER、PRF3D 等。在三维资料的处理中实际上在第 12 步之前做的都是二维处理，将三维数据体抽取成单源单缆数据，真正的三维处理是从 HABIN 开始的，而前面的 EGRID 是一个三维观测网格，将 CDP 号赋值给道头 4，因此为保证后面三维处理的准确性，在进行二维处理时我们不能改动道头 4，因此需要利用MODET或者GEOUP对抽取的每一个二维数据体定义一个二维网格，一般将二维 CDP 号赋给道头 52。8 海上 2D 资料处理作业实例详解该部分按照一系列实际海上 2D 资料处理流程的作业文本来

46、介绍作业的详细内容，主要包括循环控制、输入、处理和输出。共偏移距面 INLINE CROSSLINE OFFSET 8.1 野外采集班报勘探方式：拖缆走航式野外炮号：51-3026；道数：48；道间距：12.5m；采样率：2ms；记录长度：6000ms；炮间距：31.9；最小偏移距：135m；因水听器故障导致前 16 道废道，所有地震数据分为 4 个 sgy 文件。坐标信息：无。8.2 数据准备将四个地震数据文件分别 rename 为标准外部数据格式文件名：P2000012011cgg.DAT P2000022011cgg.DAT P2000032011cgg.DAT P20000420

48、用去远偏移距多次波的作业，8.4 作业详解 8.4.1 SEGIN 作业文本：注释：*GVRPA unsetenv SERVTAPE*环境变量设置，该行说明读取数据源非带文件格式*GVRPA setenv gvr_stage/data1/a1a2/*环境变量设置，该行给出读取数据的路径*LEVEL 5000*软件版本号*LIBRI SI 01 CREW2011cgg,P 0001,PREFIX=20,F1,STG2,CREW2011cgg,P 0002,PREFIX=20,F1,STG2,CREW2011cgg,P 0003,PREFIX=20,F1,STG2,CREW2011cgg,P 00

49、04,PREFIX=20,F1,STG2,*LIBRI SI即 LSI 库，为数据路径的调用库模块，01为库模块代号，后续调用标识；*CREW后跟工区名，P为外部文件标识，同时为盘标识，后跟文件编号，这些确*定要读取数据的文件名。PREFIX为文件编号前缀，F1 为盘文件号。CGG 在存储数据*是给每个 DAT 数据选择*一个磁盘空间，P为 10G 容量的磁盘空间，每个预分配的盘*内只允许存在一个文件。该作业中需要读入四个地震数据，因此给四个文件名标识，执行*结果是按顺序将四个数据体首尾相接在一起，最后的的 STG2标识外部磁盘文件。对*应的还有STG和STG1均为内部磁盘文件。*LIBRI

50、FI 01 SEGY,FOR3,LAB0,CGG22(1,32)=HTR3,CGG30(1,32)=HTR60,*LIBRI FI即 LFI 库，为读取数据的格式控制库模块，SEGY格式标识，FOR3*LAB0是，CGG22(1,32)=HTR3将 SEGY 的 HTR3 赋值给 CGG 内部的道头字 22。.*DLOOP 1*SEGIN AA RL6000,SI2,MAXERR20,LFI1,LSI1,*RL记录长度（ms），SI采样率（ms），MAXERR20允许读取错误的最大值，超*出该值报错，LFI1 LSI1库连接标识*SELTR AA AB SELEC=MOT17=(17-48),

展开阅读全文