超级CGG培训资料.doc

CGG LEVEL 5000培训资料CGGLEVEL 5000培训1CGG简介CGG是目前国际较为通用的一款地震数据处理软件，与其他地震数据处理软件相同，它将所有处理功能模块化，通过多个模块的搭配使用来实现某一种处理功能，一系列模块的搭配称之为一个JOB，一系列JOB的组合来完成对整个地震资料的处理。相对PROMAX，CGG具有更高的交互水平，尤其是新版本的CGG，增加了作业流程图和作业文本的同步互译功能，使得操作更加便捷。CGG具有更友好更强大的QC功能，它拥有单独的QC工具，在此工具下能实现各种形式的查看、相关数据运算、滤波、谱分析等功能，使得QC变得更加灵活实用。另外，CGG将被处理数据按照用途的差异进行了详细规划，具有更清晰的数据管理机制，便于数据的调用和管理。但是，相对PROMAX来说，CGG却又显得不够干净利索，对于PROMAX的使用者来说，初接触CGG会给人一种臃肿的、拖泥带水的感觉，尤其是文本作业的形式，让人有一种在写程序的感觉，不过随着学习的深入，我们会逐渐发现其便捷所在、优势所在。在此特别注意的是，5000版本的CGG增添了作业流程图和作业文本的同步互译功能，这一功能存在互译不稳定性，对于一些模块及模块之间的数据链接不能正确翻译为作业文本，此时我们需要以文本为主，严格按照文本格式进行修正，保证文本准确即可。但这一缺陷对CGG的交互水平影响不大，同步互译功能自有其妙用。以上概述，在下面的讲解中会得到详细的体现。2数据文件分类2.1 辅助数据文件XPS，是CGG内部专门为辅助数据而设置的文件存储格式，如切除文件、拾取水底资料、速度文件等地震属性的非地震数据主体文件库，分类细致，具有独立的QC窗口，具备列表和成图功能，便于查看和分析。很多模块中具有调用xps文件的接口，使用方便。XPS文件由PROJECT（预存储工区名）、ID和VERSION组成，要求每一个名称长度不能超过15个字节。2.2 作业文本文件CGG中，作业文件分为3种，最新版本的.job文件、老版本的.xjj文件和.gsl的作业文本文件，存储于工区目录下的JOBS目录下。2.3 库文本文件CGG对库文件具有很详细的划分，每一个库模块对应一个库文件类型，例如，.lmu(切除库)、.lvi（速度库）、.lfd（水深库）、.lcn（滤波库）、.ldb（输出库）等等。有很多库模块在调用库文件时具有双重选择，既可以选XPS文件也可以选择文本库文件。2.4 地震数据文件CGG中，地震数据根据输入输出的不同分为多种，.DAT文件同时用于输入输出，用前缀P和I来区分，分别为外部文件和内部文件，具有非常严格的格式要求，格式如下：外部：内部：外（内）前缀编号工区号.DAT或外（内）编号工区号.DAT， 前缀和编号又有严格规定，或者位前缀个数字，或者直接用个数字。外部的DAT数据可由模块SEGIN读入，外部的.sgy数据需改名整理成标准的外部DAT数据才能够被CGG识别读入。内部的DAT数据可由OUTBD输出或由INPTR读入。DAT数据一般只出现在外部读入和最后输出。.cst文件是内部的临时地震数据体，一般用于QC，主要是测试处理效果等，暂存在临时数据存储的目录下，由模块WUNET和RUNET输出和读入，称其为网络文件，不作为地震数据的主体。SDS文件，是真正的地震数据文件，类似于XPS文件无格式后缀，CGG为SDS单独在工区中开辟一个存储空间，存储路径由PROJECT（当前工区名）、ID和VERSION组成。是在整个流程中处理和输出的对象，包含全部的地震数据。由模块SETRA和GETRA输出和读入。2.5 导航定位文件.P190文件是野外采集的导航定位文件，属于外部文件，包含野外采集各炮点、检波点的坐标及沉放信息以及声纳采集的海底信息，是地震数据的一个重要部分，有专人专门制作，由模块WILOC将其读入并转化为内部格式.DAT文件，CGG对其文件的命名和路径没有特别要求。2.6 其他文件.list是作业执行列表，从该文件中可以查看作业执行进程、错误原因和执行结果信息等，是一个非常重要的文件，每一个作业执行后都会生成一个.list文件。.log系统内部的操作日志，没什么用。.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数据（SGY数据）和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数据时为由导航数据确定的标准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：通常为空值，习惯将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亦可进入。工区界面中一般主要包含以下：属性显示1.FastQC2.MapQC3.Geodesk2 作业管理程序（批处理）4.Jseg2FI 对野外地震数据进行分析5.Sditr 折射波静校正（海上处理时只在OBC做）6.TeamView 显示工具（质量控制，包含部分处理功能，如滤波、自相关等）7.Xps 辅助地震数据库（如切除、滤波、速度等文件都可以存放在这里）8.Chronovista 速度分析9.xjob 交互作业（老版本）10.Jxjob 交互作业（新版本5000）11.Xenu队列管理（选择发送作业的节点）12.Dvimon 作业监控（监控作业运行情况）13.Isub 发送作业（也就是执行作业）14.Xopreg 资源管理器15.Xdoc 说明书可以右键添加（Application）常用的，例如：GRILLTH 建立网格下面简单介绍一下TeamView和Chronovista。3.2 TeamView显示工具（质量控制）TeamView主要是用于质量控制（QC），类似于Promax中的Trace Display，但比Trace Display的功能要强大许多，不仅可以显示各种分选形式的地震数据，还可以对地震数据做一些简单的处理，如：加、减、乘、除、自相关等运算。1.地震数据的显示File-creat project 创建临时工区，当对这个工区做出改变后，退出时才会提示是否保存。点击工区名左边的圆圈，出现下拉菜单，左键点击seismic，右键选择load（一定要先左键点击seismic才行）加载地震数据，由于TeamView主要是用于质量控制（QC），所以加载的数据一般比较小，一般是.cst的文件（WUNET输出的数据），比较大的数据（例如最终的叠加剖面）一般是.sds的文件。数据加载完后，把数据用左键拖到右边的显示窗口即可显示。右键点击剖面可以进行相关显示参数的修改，以达到最佳的显示效果。在进行AGC显示时，把Samples改为Seconds，还要勾上Apply AGC。对数据进行分选显示，左键选中所加载的地震数据，右键选择分选DATA BOSRT，例如按炮集、CDP道集显示。多剖面显示时，如果想把对某一剖面的显示参数应用到其它剖面，可以直接点击传播图标即可。如果对两个剖面作对比显示，例如对去多次波前后的叠加剖面做对比显示，先左键选中两个地震数据，然后拖到显示窗口中或右键Multiple Display，点击左上角的View，再左键选中两个地震数据，右键选择最大化显示，修改显示参数，得到最佳显示效果后然后点击传播图标，再左键选中两个地震数据，右键选择movie（动画显示）即可。2.地震数据的简单处理在TeamView中可以对地震数据做一些简单的处理，如：加、减、乘、除、自相关等运算。自相关主要用于检测反褶积中Gap（间隙值）是否合适，标准是旁瓣越小越好；减法可以检测处理前后两个数据体的差异。在TeamView中还可以对线性噪声做速度的拾取（去线性噪声的时候需要用到），对地震数据做各种切除，拾取水底库等，也可以对地震数据做频谱分析，这里需要注意的是，对滤波后的数据做频谱分析得到的结果与对没有滤波的数据做频谱分析得到的结果是一样的，因为这里的滤波只是一个显示，并没有改变原始的数据。3.3 Chronovista 速度分析在做速度分析之前，必须要先由VSEPA得到一个.velcom文件，类似Promax里速度预计算得到文件，这里主要介绍一下CGG进行速度分析的步骤：Project-Creat（产生一个临时目录，目录的第一个字符必须是字母，路径必须是proj/当前工区名/GEM，系统只认这一目录）-点Geovel（双击）-Warning-OK（在这里定义一个观测系统，对于二维数据Topography是虚拟的，设置Geometry时必须将所有的线号和CDP号包含进来）-OK-Geovel-Datamanage-velocity-Analysis-ImportFrom（Velcom file）-List（找到.velcom文件）-Load-关闭当前窗口-左键选中加载后的.velcom文件，再右键-view-进入速度普拾取界面（每次拾取速度前将右上角所有工具全选中）-Velocity-new-建立一个速度文件（这一步非常的重要，否则无法拾取速度）速度拾取完后先保存，然后需要导出，右键你拾取的速度文件名Export -弹出一个对话框-即可按文本文件输出也可按.xps文件输出，一般选择按.xps文件输出，得到的速度文件一般还需要人工修改，主要包括扩边（扩展到第一个CDP号和最后一个CDP号）、设置初速度（一般为海水的速度）和末速度（最大速度）、允许速度反转（加上TN，不加则默认为反转50%的速度）。3.4 Geodesk2作业管理程序该工具主要用到的是其批量处理的功能，通过表格控制作业的批量执行，该部分在三维资料处理中具有重要意义。详细内容待补充。4常用模块功能简介SEGIN ：外部地震数据的输入。只能输入SEGY或SEGD格式的地震数据，可以同时输入多个外部地震数据。SEGIN对输入进来的地震数据的命名有严格的要求，必须按照CGG的命名规则对输入进来的地震数据重新命名，命名规则：P前缀编号工区号+.DAT 前缀和编号又有严格规定，或者位前缀个数字，或者不要前缀直接用个数字。与SEGIN相对应的还有两个库：格式库（LIBRI FI）和输入库（LIBRI SI）。由于SEGIN读入的地震数据没有给出具体的路径，只给出了盘号和编号，所以还需要设置一个输入环境变量（即该文件的存放路径），例：*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）或WILOC后得到的导航数据（如IWI10013DSRME8.DAT），与之对应的有一个TR库。也需要设置一个环境变量。GETRA：一般读取.sds的地震数据，需要给出Project（工区名）+Identity(目录)+Vesion（文件名），如果不清楚该文件的Project+Identity+Vesion，可以去SDS库中找到该文件，然后双击会出现该文件的属性，从其属性中可以查找到。 可以对数据进行分选读取，例：WHERE SERIES MOT2=(1001-1756,I2)，按炮号读取，从1001炮到1756炮每隔两炮读一次。RUNET：读取网路文件（.cst文件），一般用于QC。OUTBD：输出.DAT格式的地震数据或导航文件，与之对应的有一个BD库，输出的文件名称按照以下规则自动命名：I前缀编号工区号+.DAT 前缀和编号又有严格规定，或者位前缀个数字，或者不要前缀直接用个数字。也需要设置一个环境变量。一般和SEGIN或WILOC搭配使用。SETRA：输出.sds的地震数据，一般和GETRA搭配使用，输出比较大的数据。由于这个版本不能正确翻译SETRA模块，所以每次还得人工修改文本，修改如下（标准格式）： * SETRA AA REPLACE DATA INTO 3DSRME8/Leo/raw_combine USING COMPRESSION=F2WUNET：输出网络文件（.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数据体置2D的CDP号。EGRID：定义三维观测系统。需要定义一个网格库（LIBRI GR）LISTE：道列表。输出一个或多个地震道的值，如炮号。QCXPS：提取各种地震属性，如水深库。RSAMP：重采样。若数据量较大，2ms采样的数据可先重采样为3ms的数据在进行后续处理，以提高工作效率。因为一般海上数据的有效信号的频率在3ms的Niquist频率范围之内，所以重采样不会影响最终处理结果的精度。注意输出最终的结果时必须再重采样为2ms的数据。SCALE：道缩放，按比例显示。BSORT：分选。MUTES：切除。RAMUR：去线性噪声（NA）或着去中远距离多次波（MA）。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+MOT17MNGTY：Y位控制，在进行多道处理时前必须加上MNGTY，一般还可以在MNGTY前加一个SELTR先把数据分选一下。如后面的处理是在CDP道集上，则先用SELTR把数据按号分选，然后再用按做控制。FXINT：道内插值（在F-X域）。一般在做某些处理时先插值再处理，这样做可以防止空间假频。REFOR：振幅补偿。一般是为了计算用的，例如计算RMS（均方根值）。HISTA：静校正（如水深校正）或者用于提高各种变换精度的运算（如在Tau_P变换前可以先做一个HISTA，提高Tau_P变换的精度）。LINTP：线性Tau-P变换。TRITA：反褶积。一般在Tau-P域做反褶积，效果要好于在时间域的。主要去层间短周期多次波。SMMOD：建立多次波模型，用于SRME（主要是去近偏移距多次波）。LASUB：道头一致性检查或做减法处理（原始数据减去模型数据），用于SRME。DUPLI：数据流的复制。把多个流程的数据输入到另一个流程中。DYNQU：增益。类似AGC。RESTR：TRVEL：产生速度道。MODVI：速度场的平滑。两个方向都得平滑，一次按线方向（MOT19,MOT4），一次按CDP方向(MOT4,MOT19)。而MODVI每次这能按一个方向对速度场进行平滑，所以得分两次对速度场做平滑，一般有两种方法实现两个方向的平滑：1， 利用Tabedit打开一个速度场（3D.lvi）文件，缺省情况下，生成表格的第一列是线号（MOT19），第二列是CDP（）号，第三列是时间（），第四列是速度（）。首先按照缺省的情况（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号排列的速度场文件进行平滑，至此速度场实现了两个方向的平滑。2， 利用TRVEL，步骤如下：LIBRI VI-TRVEL-LT速度道（MOT19,MOT4）-MODVI-BSORT（MOT4,MOT19）-MODVITIKIM：偏移。HFNMO：动校正或反动校正。STAPA：保持振幅叠加（有振幅加权）。STACK：水平叠加（无振幅加权）。IFTHN：条件判断。2D SCAC 地表一致性振幅校正QCAMP：计算均方根值。GAINA：计算增益值。GAINX：地表一致性增益补偿。HABIN：面元均化。某一面元内的某共偏移距范围（例如：0-275）内的道数在理想的情况下只有一道或0道，但是实际由于拖缆摆动的原因，共偏移距范围内的道数可能有多个，通过HABIN可以把多余的道删掉，只保留离面元中心最近的那一道。这就保证了某一面元内的某共偏移距范围内的道数只有一道或0道。REG3D：数据规则化。由于HABIN之后某一面元内的某共偏移距范围内的道数可能是0道，REG3D通过插值，是得每一个面元内的每个共偏移距范围内的道数都是一道，并把每一道都移到面元中心上去，这样就达到了覆盖次数均一化的效果，避免偏移噪声。LAYES：HDPIC：高精度速度场的拾取。FXYNA：F-XY域随机噪声衰减（3D）SPARC：PRF3DHABILMULINRSAMPSDICO：球面扩散校正。 3DSRMEEXPEN：剩余振幅补偿。SYMUL：MOMUL：TRASY：TDNFK：5作业格式介绍（交互格式和文本格式）CGG作业流程中所调用的模块分为两类，即库模块和功能模块，当流程中读入或调用数据时需要通过调用功能模块并连接相应库模块来实现。库模块一般只用于标识所调用数据的路径。CGG作业流程根据数据体的形态分为诸多循环，流程中每次分选、域转换、插值、重采样等改编数据体形态的操作之后都要重启一个循环，针对同一形态数据体的所有处理一般情况下均在同一个循环中。不同循环间的数据接口一般用标识循环代号的数字，循环内的数据接口一般用大写字母或+ += =来标识，当用+ += =时，+ +一般标识输出，= =一般标识输入，但因其单一性只能完成上下一对一的数据传递。在编写流程时，我们往往先将所需要的模块按顺序拖入交互界面形成一个完整的流程图，这样在文本中便会自动正确划分循环（但也会有特殊情况的出现，如： 做LASUB时会发生翻译紊乱），自动正确设置数据流接口，循环和数据流接口建立好之后再设置各个模块的参数，在此需要特别注意的是，交互模式填写的参数有时不能正确翻译为文本，反之亦然，这就需要手动修改文本，因为我们的作业是以文本为主的，交互流程图仅仅是为了增强软件的交互能力，因此交互流程的紊乱我们可以不管，但文本一定要正确，这也是CGG流程文本互译时不稳定性的所在，但只要文本正确就不会影响到我们的处理效果。作业文本是DOS格式，*打头表示注释掉后续行，文本框的下拉菜单中有批量注释掉的功能，在此要特别注意的是，该注释并不同于PROMAX中的注释，PROMAX中做完注释后会自动重新规划数据接口，CGG文本中做完注释不会自动重新规划数据接口，需要手动修改接口。而在交互流程图中又没有注释掉某模块的功能，只能进行删除或切除操作，但这两中操作是具有数据接口自动重设功能的。因此需要注释的时候我们一般还是在文本中操作。一个示例作业文本的格式如下：* LEVEL 5000*环境变量设置，CGG版本号* PATCH S_HAL5000_1322* PATCH P_TIKIM5000_7203*环境变量设置，补丁* GVRPA 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,YOR4242958.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,LLINE905,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#*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)去线性噪音，正反两个方向都做，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+GAINX，路上资料用NMBSO模块。10、REFOR，振幅一致性处理，消除单道内的振幅衰减，SDICO，球面扩散校正，消除球面上的振幅衰减。11、PSTM，偏移处理，计算偏移速度场。12、VESPA+PSTM，第一次偏移后重新拾取速度进行第二次偏移速度场计算。13、STACK。14、叠后处理，FXNET+FILTER。15、AGC，LFK去除叠后浅层噪音，如海底抖动。第3步中，如果资料道数有限，无法辨别线性噪音则省去线噪去除，如果叠后资料中存在线性噪音，可以叠后去线性。第7步，同样当资料道数有限时就无远偏移距可言了，如第八节中所讲到的实例中便没用到。在此特别介绍本流程中去海底多次波的思路如下：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，提取水底库。根据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号可以直接用MODET来计算，也可以用GEOUP来划分，因为三维资料每条缆的最小偏移距是不等的，因此推荐用MODET更便捷。8、DECON， T-P预测反褶积消除层间多次波。9、VESPA，获得速度资料V1。10、RAMUR，出去中远偏移距海底多次波。11、2D SCAC，炮集上做地表一致性处理，该处理要对每一条航行线去做，每条航行线上做一个GIANA，应用到对应的航行线上，一条航行线包含所有的CDP轨迹线，二源六缆则包含12条线的数据。12、HABIN，共偏移距组划分，并筛选CDP覆盖道，处理结果是将数据划分为N组，每组内的CDP点上最多存在一道覆盖。该过程将所有数据按照偏移距划分为N组，每组均包含位于一定偏移距间隔范围内的数据，原则上一个偏移距道集内每一道数据均对应一个CDP，即共偏移距道集所能够覆盖到的CDP原则上仅能被覆盖一次，因此我们通常将道集内包含多次覆盖的CDP点作剔道处理，一般保留距CDP中心最近的那一道，剔道完毕对道集中的每一道求取均方根给已建立的共偏移距CDP网格赋值，存在覆盖的便存在值，覆盖不到的便显示空缺（要注意的是所有共偏道集的CDP网格是统一的，即按照最小坐标值和方位角确定的那个网格）。假设资料的满覆盖次数为N次，则理论上应该分为N个共偏移距面（组），组距原则上应该为两倍的炮检距，但实际上野外采集时缆是漂动不定的，因此要按实际坐标的观测系统求取，公式为：，其中XRM、NRM分别为实际的最大、最小偏移距。理论上第一组偏移距范围应为0，NRM，但因为每一条缆的NRM均不同，因此应该选择其中最大的NRM，基本原则就是第一个共偏移距组要将所有缆的首道包含在内，确定了第一个共偏移距面后，后续的共偏面按照划分即可，参考图7-1理解。共偏移距面INLINECROSSLINEOFFSET图7-1 共偏移距面的划分图示13、REG3D，3D数据规则化，在HABIN的基础上进行规则化处理，该处理结果有二，一是将每个共偏移距面上的CDP覆盖道进行规则移位，移动到每个CDP面元的中央；二是将内部零覆盖的CDP面元上进行插值，使其具有一道的覆盖次数；这样就保证了每个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，当然拾取前可以视情况做线性噪音的去除。17、HFNMO，STAPA。18、叠后处理，AGC、FILTER、PRF3D等。在三维资料的处理中实际上在第12步之前做的都是二维处理，将三维数据体抽取成单源单缆数据，真正的三维处理是从HABIN开始的，而前面的EGRID是一个三维观测网格，将CDP号赋值给道头4，因此为保证后面三维处理的准确性，在进行二维处理时我们不能改动道头4，因此需要利用MODET或者GEOUP对抽取的每一个二维数据体定义一个二维网格，一般将二维CDP号赋给道头52。8海上2D资料处理作业实例详解该部分按照一系列实际海上2D资料处理流程的作业文本来介绍作业的详细内容，主要包括循环控制、输入、处理和输出。8.1 野外采集班报勘探方式：拖缆走航式野外炮号：51-3026；道数：48；道间距：12.5m；采样率：2ms；记录长度：6000ms；炮间距：31.9；最小偏移距：135m；因水听器故障导致前16道废道，所有地震数据分为4个sgy文件。坐标信息：无。8.2 数据准备将四个地震数据文件分别rename为标准外部数据格式文件名：P2000012011cgg.DATP2000022011cgg.DATP2000032011cgg.DATP2000042011cgg.DAT存放于路径/data1/a1a2/之下。8.3 处理流程图SEGIN（数据导入）|GEOUP（建观）|PREPROCESS（MUTE、FILTER、RENOISE）（预处理，切除、滤波、去野值和其他随机噪音）|SRME(SMMOD、LASUB)（去近偏移距多次波）|TPDEC（FXINT、LINTP、TRITA）（t-p域反褶积，道内插、T-P转换、T-P预测反褶积）|VESPA（创建速度谱文件，拾取速度谱）|NMO（动校正）|TIKIM（TIKIMIN、TIKIMMI）（偏移处理）|STACK（叠加剖面） 说明：因为该资料最大偏移距稍大于700m，因此在流程中没有用去远偏移距多次波的作业，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 0004, PREFIX=20,F1, STG2,*LIBRI SI即LSI库，为数据路径的调用库模块，01为库模块代号，后续调用标识；*CREW后跟工区名，P为外部文件标识，同时为盘标识，后跟文件编号， 这些确*定要读取数据的文件名。PREFIX为文件编号前缀，F1为盘文件号。CGG在存储数据*是给每个DAT数据选择*一个磁盘空间，P为10G容量的磁盘空间，每个预分配的盘*内只允许存在一个文件。该作业中需要读入四个地震数据，因此给四个文件名标识，执行*结果是按顺序将四个数据体首尾相接在一起，最后的的STG2标识外部磁盘文件。对*应的还有STG和STG1均为内部磁盘文件。* LIBRI 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