二维处理流程的学习资料课件.ppt

1、GRISYS GRISYS 处理系统处理系统二维地震资料处理二维地震资料处理处理流程的建立处理流程的建立 1 1、了解、了解地质任务地质任务，工区所在位置的地理位置、地表地震条件、，工区所在位置的地理位置、地表地震条件、构造和地层特征，向用户索取必要的相关资料（近地表信息）。构造和地层特征，向用户索取必要的相关资料（近地表信息）。2 2、对将要处理的数据进行、对将要处理的数据进行调查分析调查分析，通过试验初步掌握数据的，通过试验初步掌握数据的信噪比，分析了解野外施工质量及数据品质在空间上的变化情况。信噪比，分析了解野外施工质量及数据品质在空间上的变化情况。3 3、针对数据的实际情况，选定处理系

2、统中对应的处理功能，并、针对数据的实际情况，选定处理系统中对应的处理功能，并对该功能进行试验，最终选定某一功能选件，即对该功能进行试验，最终选定某一功能选件，即针对性措施针对性措施。4 4、结合以往的处理历史档案（新资料首次处理，应结合该地区、结合以往的处理历史档案（新资料首次处理，应结合该地区相临资料的处理历史档案），征求用户的意见，相临资料的处理历史档案），征求用户的意见，建立处理流程建立处理流程，同，同时由处理责任人审核批准，完成时由处理责任人审核批准，完成“处理设计处理设计”。一个好的处理流程应具有：科学性、合理性、针对性。一个好的处理流程应具有：科学性、合理性、针对性。数据带解编数据

3、带解编STEP 1STEP 1GRISYS系统输入输出模块SEGBSEGD ：DBAIFDBAISDEMULTSEGD ：SEGDINSEGY(tape+disk)：SEGYIN、SGYOUTGRISYS(TAPE)：TAPEIN、TAPOUTGRISYS(DISK)：DISKIN、DSKOUT、SORTIN DSKTSO、DATINP、DATOUTCGG ：CGGINP、CGGOUTWGC ：CODEIN、CODOUTSEGDIN SEGD格式磁带解编 为配合海上数据处理，对SEGDIN模块进行了较大规模的改进，增加海上SEGD格式数据读取的功能，并能准确地从原始数据中读取相关数据并存放在对

4、应的GRISYS数据道头字中。GRISYS/WS-V8.0 GRISYS/WS-V8.0 改进内容改进内容 增加增加MPMP参数：参数：为乘为乘MPMP因子输入标志。因子输入标志。这是一个选择性参数，填此选件后，程序这是一个选择性参数，填此选件后，程序从数据的头块中求出因子的实际值与仪器从数据的头块中求出因子的实际值与仪器记录数值相乘。从而，得出近似仪器增益放大记录数值相乘。从而，得出近似仪器增益放大前的数值。一般数据输入时不用此参数。仅当前的数值。一般数据输入时不用此参数。仅当做数值分析输入时慎用，缺省为不应用因做数值分析输入时慎用，缺省为不应用因子。同时，取消子。同时，取消TRUETRUE

5、参数。参数。SEGDINSEGDIN GRISYS/WS-V8.0 GRISYS/WS-V8.0 改进内容改进内容 GRISYS/WS-V8.0 GRISYS/WS-V8.0 改进内容改进内容 道头增加内容为：道头增加内容为：HD80=水、陆检排列的标识。水、陆检排列的标识。HD81=海洋资料采集拖缆的缆号。海洋资料采集拖缆的缆号。HD82=分放分放1 1，2 2，3 3，4 4排列。排列。HD83=接收点站号。接收点站号。HD84=接收线号。接收线号。HD85=炮线号。（炮线号。（SEASEA参数起作用）参数起作用）HD86=炮点站号。（炮点站号。（SEASEA参数起作用）参数起作用）HD8

6、8=磁带卷号。磁带卷号。HD89=通道号。通道号。HD90=唯一道号。唯一道号。SEGDINSEGDIN磁盘和磁带数据分析磁盘和磁带数据分析TPUITAPUILISEGYADISGYTAPE 8.0新增新增TPU、ITAPUIL 磁带数据分析磁带数据分析 两者功能基本相同，两者功能基本相同，TPU为早期的版本，为早期的版本，TAPUIL为交互版本。功能主要有：为交互版本。功能主要有：磁带数据格式自动分析；磁带数据格式自动分析；磁带数据磁带数据DUMP；磁带数据拷贝。磁带数据拷贝。其中磁带数据拷贝又提供多种拷贝方式：带对带、其中磁带数据拷贝又提供多种拷贝方式：带对带、带对盘、盘对带、合并拷贝以及

7、控制拷贝等等。带对盘、盘对带、合并拷贝以及控制拷贝等等。ISEGYAD 交互交互SEGY数据显示数据显示 本模块的主要功能是在本模块的主要功能是在GRISYS系统中，增加对系统中，增加对SEG-Y数据数据磁盘文件的监控、对未知数据文件的格式分析等。磁盘文件的监控、对未知数据文件的格式分析等。具体功能为：具体功能为：（1）对数据文件打开的同时，对数据文件进行格式分析；）对数据文件打开的同时，对数据文件进行格式分析；（2）对）对SEG-Y数据体显示；数据体显示；（3）对）对SEG-Y数据数据“C卡卡”的显示、修改和转存；的显示、修改和转存；（4）对）对SEG-Y数据数据“400字节字节”的文件头块

8、主要内容的显示；的文件头块主要内容的显示；（5）对）对SEG-Y数据数据“240字节字节”的道头主要内容的显示。的道头主要内容的显示。交互SEGY磁带数据检查 ISGYTAPEISGYTAPE 在GRISYS/WS-V7.0中推出了SEGY磁盘数据检查的交互模块。此次版本升级我们将SEGY磁带数据检查补充进来，使得SEGY数据检查功能更加完善。其具有以下几个特点：数据DUMP；3200和400字节头块检查；显示参数可调；道头数据即时显示；任意道头数据曲线绘制。GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容 GRISYS/WS-V8.0 GRISYS/WS-V8.0

9、 新增内容新增内容主界面及磁带机选择磁带数据显示即时道头信息读取交互SEGY磁带数据检查 GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容交互SEGY磁带数据检查 其他功能3200 头块400 头块数据DUMP观测系统定义（专题）观测系统定义（专题）STEP 2STEP 2道编辑道编辑STEP 3STEP 3 随着地震勘探技术和计算机技术的不断发展，地震数随着地震勘探技术和计算机技术的不断发展，地震数据体也变得越来越大。为了适应不断增长得地震数据体对据体也变得越来越大。为了适应不断增长得地震数据体对处理系统的需求，此次版本升级对于大数据量数据的处理处理系统的需求，此

10、次版本升级对于大数据量数据的处理进行了较大的改进，并新开发了两个针对大数据体显示和进行了较大的改进，并新开发了两个针对大数据体显示和编辑的模块。编辑的模块。新版的新版的IMUTEDIIMUTEDI；交互地震数据显示模块交互地震数据显示模块IDATVIEWIDATVIEW。GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容 GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容交互地震道切除与编辑 IMUTEDIIMUTEDI（NEWNEW）新版的IMUTEDI模块除了沿袭了原有模块的主要功能以外，具有使用方便，支持大道集、大数据体的输入和快速定位

11、等特点。索引文件引入，实现大数据体的快速定位检索；支持大道集，缺省以2000道为限显示页面；双滑动杆设计，一个为数据体滑动，另一个为道集内部滑动；切除库和道编辑库直接读入。GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容启动启动IMUTEDI数据加载数据加载建立索引文件建立索引文件数据操作数据操作保存道编辑参数库保存道编辑参数库具体操作流程退出退出 GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容IMUTEDI（新）（新）GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容IMUTEDI（新）（新）GRISYS/W

12、S-V8.0 GRISYS/WS-V8.0 新增内容新增内容IMUTEDI（新）（新）GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容IMUTEDI（新）（新）GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容IDATVIEW 较好地适应大数据体输入；较好地适应大数据体输入；独立的数据索引方式；独立的数据索引方式；数据的快速准确定位；数据的快速准确定位；无极放大功能；无极放大功能；支持双屏独立操作；支持双屏独立操作；自动记忆自动记忆10屏显示。屏显示。交互地震数据显示交互地震数据显示启动启动IDATVIEW数据加载数据加载建立索引文件建立索

13、引文件数据操作数据操作退出退出具体操作流程 GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容IDATVIEWIDATVIEW GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容 GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容IDATVIEW GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容IDATVIEW GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容交互地震数据显示和对比 IDATCOMP 在GRISYS/WS-V7.0中推出了绘图文件对比的交互模块，

14、该模块在地震资料处理中监控处理质量和处理效果起到了很好的作用。但是实际上，地震数据的直接对比更能反映地震数据真实的变化。为此我们此次推出了GRISYS地震数据的对比显示工具。其主要功能是：双数据输入，独立或同步翻动比较；显示参数统一，并支持多种显示方式；鼠标点击切换数据，快捷方便；支持不同的GRISYS数据格式。GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容主界面数据类型定义数据显示及对比显示参数定义不同类型显示方式野外静校正野外静校正折射波静校正折射波静校正(专题专题)STEP 4&5STEP 4&5子波整形子波整形(可选可选)STEP 6STEP 6 GRIS

15、YS处理系统中，子波整形由两个模块来完成，一处理系统中，子波整形由两个模块来完成，一是整形滤波器设计（是整形滤波器设计（SHPFIL），二是），二是 整形反褶积（整形反褶积（SHDECN）。）。SHPFIL：同一工区采集的资料，由于使用不同震源或者使同一工区采集的资料，由于使用不同震源或者使用不同记录系统施工，致使信号在工区内不能很好地连续对用不同记录系统施工，致使信号在工区内不能很好地连续对比，使用子波整形校正可消除子波之间的差异，从而改善处比，使用子波整形校正可消除子波之间的差异，从而改善处理效果。理效果。首先输入参考数据道以及需要校正的数据道首先输入参考数据道以及需要校正的数据道，指定参

16、与，指定参与计算子波的参考道和需要校正道的范围和时窗大小，近似计计算子波的参考道和需要校正道的范围和时窗大小，近似计算出参考道的子波和需要校正的道的子波，利用维纳滤波器算出参考道的子波和需要校正的道的子波，利用维纳滤波器求出这两种子波的差异，即滤波因子。求出这两种子波的差异，即滤波因子。：JOB P：DQ,L：595,T：332：DISKIN DFN:GZHY0001，(参考道参考道)DFN:GZHY0000，NTY2 (要校正的道要校正的道)：SHPFIL MU02,CMP97,TZ4950,CMP100,TZ4950,LW380,ACMP96,ATZ4950,ACMP98,ATZ4950,

17、ALW380,OUT1，SHFN:SHAB1,PSTK,SHOT,MAX4000,MIN100,AMAX4000,AMIN100,L100,TAP50,NOIS1,LIST,FOLD120,HD20,ALL,OUT0,PLOT END 把整形滤波器设计应用模块（把整形滤波器设计应用模块（SHPFIL）求得的整形因子，应用到数据）求得的整形因子，应用到数据上，达到子波一致的目的上，达到子波一致的目的 。：JOB P：DQ,L：595,T：332 ：DBAMU VN01,HD4,(1-500)=M100X50,M1000X3000,TAP10 ：DISKIN DFN:GZHY0000,LTR300

18、0,HD2=(195),NTY2,STOP#：SHDECN (-0.05,0.015,-0.122,-0.083,-0.123,0.0026,-0.049,-0.033,-0.072,-0.004,），），LIST ：SHDECN SHFN:SHAB1，LIST ：DSKOUT DFN：GZHY0001，FMT4 ：END参考炮参考炮子波整形前单炮子波整形前单炮参考炮频谱参考炮频谱子波整形前频谱子波整形前频谱参考炮参考炮子波整形后单炮子波整形后单炮子波整形后频谱子波整形后频谱参考炮参考炮频谱频谱845Hz滤波后滤波后 存在相位差存在相位差子波整形前叠加剖面（子波整形前叠加剖面（INLINE）子

19、波整形后叠加剖面（子波整形后叠加剖面（INLINE）子波整形前叠加剖面（子波整形前叠加剖面（INLINE）整形后的处理：整形后的处理：子波整形后的数据，应进行地表一致子波整形后的数据，应进行地表一致性反褶积或多道统计子波反褶积处理，以性反褶积或多道统计子波反褶积处理，以此来进一步完成子波的统一。因为在预测此来进一步完成子波的统一。因为在预测步长较大时，地表一致性反褶积的主要功步长较大时，地表一致性反褶积的主要功能是能是“子波调整子波调整”，压缩子波的处理由后，压缩子波的处理由后续的统计子波反褶积或单道反褶积等模块续的统计子波反褶积或单道反褶积等模块来完成，即串联反褶积。来完成，即串联反褶积。叠

20、前去噪叠前去噪(专题专题)STEP 7STEP 7 GRISYS/WS-V8.0 GRISYS/WS-V8.0 改进内容改进内容IPRTEST新增7个模块，分别为DBAFL、PVNOEL、LINOEL、FILTER、SFNPRE、POLFIT、SCALER；新增模块试验方式与原有模块基本相同；增加交互拾取拟合参数功能。振幅补偿振幅补偿STEP 8STEP 8NXn n为要处理的炮检距个数为要处理的炮检距个数 NTm m为对每个炮检距所规定的叠加道数为对每个炮检距所规定的叠加道数 DB 填此标识时，采用对数比例绘图填此标识时，采用对数比例绘图 空单道平衡选件A整道平衡选件B多时窗区域平衡C均方根

21、平衡选件D瞬时平衡选件E单炮平衡选件 E：在用户定义的单炮记录道的规定时窗内，求取平均振幅值，再根据提供的区域平衡级别系数计算出相应记录的平衡系数，将其与实际振幅值相乘，求其平衡后的振幅值。timeAtimetimeAtime一次速度分析剩一次速度分析剩余静校正（专题）余静校正（专题）剩余静校正见静校正专题剩余静校正见静校正专题STEP 9STEP 9 速度分析叠加速度分析 常速扫描速度分析速度分析数据解释速度谱解释常速扫描解释 层位拾取 速度处理速度换算：均方根速度换算层速度，等时速 度内插，沿层速度内插DBAPV插值交互速度编辑与修改常速扫描叠加的速度与现有速度场的结合功能介绍功能介绍速度

22、分析所要求的相应处理功能NMO 叠加其它辅助性功能TV 函数加载TV 函数输出Mute 函数加载Mute 函数输出速度处理和速度场质控图件显示TV函数剖面展开显示叠加速度二维彩色剖面显示速度等值线图显示速度谱解释窗口速度谱解释窗口常速扫描解释窗口常速扫描解释窗口7.0 对新增和改进部分的说明 作csum 大道集的选项,便于质量控制;条件约束自动寻找极值点 速度微调 tv函数展开显示 谱线 +变面积显示（同批量谱显示）常速扫描叠加的速度场可选择修改扫描范围内的时间速度对 Svel 简单操作过程l在%状态下或 gscons下启动svel，弹出主界面n按照二级管理（在主界面tool下project

23、manage 项）加载数据，同一地区的数据在load中可同时加载多条测线数据。n主界面File选项，打开本次分析的数据。l主界面Seisfunc功能项中，选VanalyAnalyse弹出作速度谱的参数框。l此参数框中有三个此参数框中有三个 参数框必参数框必须人工干预须人工干预:l CMPs for vel.analyse Parameters lMuting function definition Parameters lRef velocity curves Parameters l以上三个参数框填写好后按以上三个参数框填写好后按OK鍵，其它参数用缺省值，鍵，其它参数用缺省值，在速度谱解释窗

24、口中显示出在速度谱解释窗口中显示出速度谱、大道集和变速叠加速度谱、大道集和变速叠加段段。TV&MU输出 主界面主界面File窗口选窗口选Svel项项,弹出窗口，在弹出窗口，在Data type栏中选择栏中选择保存类型。如保存类型。如:TV库，选库，选TVPair，在，在Name对应的参数框填输对应的参数框填输出文件名后点出文件名后点 OK，存成内部文件。，存成内部文件。主界面：主界面：ToolTV I/O选项，弹出窗口，点击选项，弹出窗口，点击OUTPUT，选，选右边栏中的文件名，变黑后，在右边栏中的文件名，变黑后，在Grisys TV Filename下边的路径下边的路径指向上，填写输出文件

25、名，点指向上，填写输出文件名，点OK，完成输出功能。此时可在当，完成输出功能。此时可在当前目录下找到填写的输出文件。前目录下找到填写的输出文件。注意事项 TV库中第一个T切除库中第一个T时间。输入数据应经过振幅恢复处理的cmp集数据 譜图太窄时，可改变Maxim Offest参数，使其小于实际偏移距。偶而出现运行不正常时，一要查看数入数据正确性，二建立新工区，重新填参数或进入/grisys/oper/optrun/svel目录中，清除本工区的所有中间文件 GRISYS/WS-V8.0 GRISYS/WS-V8.0 改进内容改进内容SVEL 8.0改进内容支持1ms采样数据；此次修改主要面向大数

26、据体；CMP数/号最大可以支持9999999；改进谱点定义方式，增加起止点及增量定义方式；速度谱一次可做512个；对于没有解释的谱点，生成TV文件时，不再输出。GRISYS/WS-V8.0 GRISYS/WS-V8.0 改进内容改进内容滑动条，滑动条，CMPCMP个数个数20002000时自动时自动分屏；滑动条上的数字为当前分屏；滑动条上的数字为当前页面的第一个页面的第一个CMPCMP号；拖动滑号；拖动滑动杆来定位数据。动杆来定位数据。新的速度谱点定义方式新的速度谱点定义方式:除了除了保留原有的定义方式之外保留原有的定义方式之外,新新增了起始点、终止点和增量的增了起始点、终止点和增量的定义方式

27、。定义方式。交互速度分析 GRISYS/WS-V8.0 GRISYS/WS-V8.0 改进内容改进内容DBATV/DBAVI 速度库定义增加速度库的时变空变比例系数加权功能。具体是：去除原有模块中的“PE”参数，新增“TPmPn”参数，从而实现速度的时空变比例加权。GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容速度平滑编辑 DBAVE DBAVE模块根据用户提供的参数对速度库进行如下处理：l 速度平滑处理l 速度值编辑修改l 速度时变比例系数加权：JOBJOB：DBAFL VN01DBAFL VN01：DBAMU VN01DBAMU VN01：DBATV/DBAV

28、I VN01DBATV/DBAVI VN01：DBAPV/DBAVS VN01DBAPV/DBAVS VN01：DISKINDISKIN：FANMOV TV01/VI01FANMOV TV01/VI01，MU01MU01：PASTAK LW300 PASTAK LW300 或者或者 ：GASTAKGASTAK：DSKOUTDSKOUT：FILTER FL01FILTER FL01：AMPEQU L300AMPEQU L300：PLOTTR SCL100PLOTTR SCL100，STP10STP10：ENDEND二维叠加作业二维叠加作业地表一致性反褶积地表一致性反褶积（可选）（可选）STEP

29、10STEP 10地表一致性反褶积地表一致性反褶积 地表一致性反褶积可以消除炮点、检波点随地表地表一致性反褶积可以消除炮点、检波点随地表条件和采集因素的变化对子波带来的影响，可输出稳条件和采集因素的变化对子波带来的影响，可输出稳定和基本一致的地震子波定和基本一致的地震子波 参数说明参数说明:使用实例使用实例:第一步第一步地表一致性反褶积地表一致性反褶积参数说明参数说明:使用实例使用实例:第二步第二步 （算子的长度不能大于（算子的长度不能大于SURSAN中的中的LW参数）参数）二次速度分析二次速度分析剩余静校正剩余静校正同一次速度分析和剩余静校正同一次速度分析和剩余静校正STEP 11STEP

30、11反褶积反褶积STEP 12STEP 12 地震记录的褶积模型地震记录的褶积模型 为地层反射系数；为地层反射系数；为地震子波。为地震子波。地震子波包含震源响应、检波器响应，这些响应是时不变的；现在一地震子波包含震源响应、检波器响应，这些响应是时不变的；现在一般检波器都是组合接收，还有组合响应（若炮组合也有组合响应），还有般检波器都是组合接收，还有组合响应（若炮组合也有组合响应），还有地层吸收，而后二者是时变的。地层吸收，而后二者是时变的。反褶积的目的是想消除地震子波反褶积的目的是想消除地震子波 影响，从而得到地下反射系数影响，从而得到地下反射系数 ，克莱伯特形象地将反褶积比喻为显微镜，使我们

31、能更细致地观察到地下物克莱伯特形象地将反褶积比喻为显微镜，使我们能更细致地观察到地下物性变化（纵向）。如有剩余子波也希望是零相位的，因为零相位子波分辨性变化（纵向）。如有剩余子波也希望是零相位的，因为零相位子波分辨率最高。率最高。iiiwsiiwiwi 计算并显示地震信号的振幅谱和相位谱计算并显示地震信号的振幅谱和相位谱：JOB P：NKK，L：LNK1，T：SPEANA：DISKIN DFN：G0000321，HD4=(150-152)：SPEANA HD4，(150-152，W500-W1000)，AMP2,LTWO：END AMP2，以增益方式绘振幅谱，以增益方式绘振幅谱 GRISYS/

32、WS-V8.0 GRISYS/WS-V8.0 改进内容改进内容SPEANA 频谱分析 改进频谱分析结果的输出方式，使得其结果和交互频谱分析的结果保持一致。GRISYS/WS-V8.0 GRISYS/WS-V8.0 改进内容改进内容1.1.可实时的对地震记录道形式的各种数据进行可实时的对地震记录道形式的各种数据进行FourierFourier谱分析，包括单道谱分析，包括单道一维谱、多道一维谱及多道的二维谱（即一维谱、多道一维谱及多道的二维谱（即FKFK谱）；谱）；2.2.在一维谱（包括单道和多道）分析的基础上进行滤波器设计和试验，滤在一维谱（包括单道和多道）分析的基础上进行滤波器设计和试验，滤波

33、器包括低通、带通、高通、带陷、雷克子波和俞氏子波；波器包括低通、带通、高通、带陷、雷克子波和俞氏子波；3.3.可进行频带扫描分析；可进行频带扫描分析；4.4.可进行时频谱分析；可进行时频谱分析；5.5.可进行多时窗的频谱对比，其中包括不同数据集上相同时窗的频谱对比，可进行多时窗的频谱对比，其中包括不同数据集上相同时窗的频谱对比，相同数据集上不同时窗的频谱对比；相同数据集上不同时窗的频谱对比；6.6.可定义并应用滤波器的作用范围及强度，从而实现时变和空变滤波；可定义并应用滤波器的作用范围及强度，从而实现时变和空变滤波；7.7.可定义切割滤波；可定义切割滤波；8.8.可对得到的滤波器建立相应的可对

34、得到的滤波器建立相应的GRISYSGRISYS滤波库，形成相应批量作业，且用滤波库，形成相应批量作业，且用户可随时查看与编辑该作业文件。户可随时查看与编辑该作业文件。GRISYS/WS-V8.0 GRISYS/WS-V8.0 改进内容改进内容IFSANFL一维谱增加分贝显示；一维谱加盖网格，并可修改显示参数；改变频率扫描结果显示方式，增大显示面积。GRISYS/WS-V8.0 GRISYS/WS-V8.0 改进内容改进内容交互频谱分析显示结果和方式与SPEANA模块一致 GRISYS/WS-V8.0 GRISYS/WS-V8.0 改进内容改进内容交互频谱分析增大了结果显示面增大了结果显示面积，

35、方便对比观察积，方便对比观察分频扫描分频扫描：JOB P：NKK，L：LNK1，T：SPEANA：DISKIN DFN：G0000321，HD2=(150-152),STOP：AMPEQU A:FBSCAN BP(7,10,20,28),BP(10,15,30,42),BP(14,20,40,56),BP(21,30,60,84),BP(28,40,80,112),BP(35,50,100,140),BP(42,60,120,168),BP(50,70,140,196),BP(56,80,160,224),M120:PLOTTR SCL100,STP10,GAP2=HD22,HDA2,HDB1

36、7：END TFSCANTFSCAN 应用一系列具有不同中心频率和宽度的三角形滤波器，应用一系列具有不同中心频率和宽度的三角形滤波器，对地震道进行滤波，将同一个地震道的滤波结果安滤波器对地震道进行滤波，将同一个地震道的滤波结果安滤波器中心频率从小到大顺序排列，可研究地震记录中不同频率中心频率从小到大顺序排列，可研究地震记录中不同频率成分在不同时间上的响应，应用于地震地层学分析，时频成分在不同时间上的响应，应用于地震地层学分析，时频补偿。补偿。：JOB P：NKK，L：LNK1，T：SPEANA：DISKIN DFN：G0000321，HD2=(150-152),STOP：AMPEQU A*:T

37、FSCAN：TFSCAN ENVELP :PLOTTR SCL100,STP10,GAP2=HD22,HDA2,HDB17：END ENVELP:指明输出地震道的包络，缺省输出地震道。指明输出地震道的包络，缺省输出地震道。HD22 本模块用两种典型子波俞氏子波和雷克子波对记录本模块用两种典型子波俞氏子波和雷克子波对记录进行滤波扫描，处理员可分析扫描结果，用于选择适当的滤进行滤波扫描，处理员可分析扫描结果，用于选择适当的滤波参数和子波反褶积的期望输出。波参数和子波反褶积的期望输出。当输入为脉冲函数时，本模块输出为参数对应的子波；当输入为脉冲函数时，本模块输出为参数对应的子波；当输入为地震记录时，

38、用俞氏子波或雷克子波对记录做滤波当输入为地震记录时，用俞氏子波或雷克子波对记录做滤波处理，输出为滤波扫描的结果。处理，输出为滤波扫描的结果。:JOB P:2D,L:001,T:T:DATGEN SI1,LTR4000,NT4,RD,D1000,Z4100,1,ZG1000:YUSCAN YU（10，60），），YU（25，150），），YU（35，300）:AMPEQU A:PLOTTR SCL0600，STP100，LR，HEAD，G0，NR,W0900-W1100:END:DISKIN DFN:HLABEL1,HD2=20,STOP：YUSCAN RICK（10，30，50，80，100，

39、200）*:YUSCAN YU（10，60），），YU（25，150），），YU（35，300）：AMPEQU L300：PLOTTR SCL0100，STP100，LR，HEAD，GAP3HD22 反褶积前准备工作反褶积前准备工作 脉冲反褶积；脉冲反褶积；预测反褶积。预测反褶积。预测反褶积预测反褶积 预测反褶积是通过预测滤波来实现的，预测滤波的基预测反褶积是通过预测滤波来实现的，预测滤波的基本思想是，设计一个滤波器，使得该滤波器具有某种预本思想是，设计一个滤波器，使得该滤波器具有某种预测能力，通过它对信号的当前值和过去值的滤波，预测测能力，通过它对信号的当前值和过去值的滤波，预测未来某个时刻

40、将要出现的信号成分。预测时刻和当前时未来某个时刻将要出现的信号成分。预测时刻和当前时刻的距离称为预测距离。预测滤波器的设计采用了最小刻的距离称为预测距离。预测滤波器的设计采用了最小平方准则，即要求预测结果和实际信号间误差的最小平平方准则，即要求预测结果和实际信号间误差的最小平方和最小。方和最小。PRDCONPRDCON是一个预测反褶积模块，既可以用于消除地是一个预测反褶积模块，既可以用于消除地震记录中的长、震记录中的长、短周期多次波，又可以用来提高地震记短周期多次波，又可以用来提高地震记录的分辨率。录的分辨率。反褶积参数选取反褶积参数选取（反褶积效果分析反褶积效果分析 反褶积的迭代反褶积的迭代

41、在在随地表条件和采集因随地表条件和采集因素的变化对子波带来的影响素的变化对子波带来的影响 后，再做后，再做，不仅考虑，不仅考虑又兼顾又兼顾 GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容 此次版本升级，在反褶积和子波处理方面也增加了一些更加此次版本升级，在反褶积和子波处理方面也增加了一些更加实用的模块，主要包括：实用的模块，主要包括：逐点反褶积逐点反褶积 多道反褶积多道反褶积 频率域加权频率域加权 振幅相位的振幅相位的Q Q补偿补偿 子波调整子波调整反褶积及子波处理反褶积及子波处理 逐点反褶积逐点反褶积 PPDCON 该方法是为了解决地震子波的时变性而设计开发的。

42、该方法是为了解决地震子波的时变性而设计开发的。由于是在逐点滑动、且在很短的一个小时窗内求算子并由于是在逐点滑动、且在很短的一个小时窗内求算子并加以应用，可以说是一种真正意义上的时变反褶积。可加以应用，可以说是一种真正意义上的时变反褶积。可应用于叠前、叠后信噪比较高的记录上，能较好的提高应用于叠前、叠后信噪比较高的记录上，能较好的提高地震记录的纵向分辨率。地震记录的纵向分辨率。GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容 GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容原始单炮原始单炮逐点反褶积单炮逐点反褶积单炮 MTDCON 采用维纳

43、滤波的方法从多道地震记录中提取子波，从采用维纳滤波的方法从多道地震记录中提取子波，从而使求取的子波较接近于实际的子波，使反褶积效果更好。而使求取的子波较接近于实际的子波，使反褶积效果更好。模块提供模块提供预测预测和和脉冲脉冲两种选件。两种选件。多道反褶积多道反褶积 GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容 GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容单道预测单炮单道预测单炮多道预测单炮多道预测单炮 GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容 FREQUE 本模块先用付氏变换把地震记录从时间域变

44、换到频率域，在频率本模块先用付氏变换把地震记录从时间域变换到频率域，在频率域内对原记录主频之后的频段内的振幅谱的振幅值进行加权处理，使域内对原记录主频之后的频段内的振幅谱的振幅值进行加权处理，使高频段的振幅值相对加强，使其频谱向高频方向拓宽，而不是向高频高频段的振幅值相对加强，使其频谱向高频方向拓宽，而不是向高频方向移动。然后再用付氏变换把地震记录从频率域变换到时间域，使方向移动。然后再用付氏变换把地震记录从频率域变换到时间域，使地震记录的高频成份的振幅值有所加强。由于本模块是在频率域内只地震记录的高频成份的振幅值有所加强。由于本模块是在频率域内只对振幅谱进行加工，不改变相位谱，所以经本模块处

45、理后的地震记录对振幅谱进行加工，不改变相位谱，所以经本模块处理后的地震记录相位特征不改变，信噪比不降低，不产生假频，只提高分辨率。相位特征不改变，信噪比不降低，不产生假频，只提高分辨率。频率域加权频率域加权 GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容模块应用前模块应用后三次速度分析三次速度分析剩余静校正剩余静校正同一次速度分析和剩余静校正同一次速度分析和剩余静校正STEP 13STEP 13叠前随机噪音衰叠前随机噪音衰减（可选）减（可选）见去噪专题见去噪专题STEP 14STEP 14多次波压制及内切除（可选）多次波压制及内切除（可选）BFMLTA BFMLT

46、A RADMUSRADMUSDBAMUDBAMUMUTINGMUTING（tail)tail)见去噪专题见去噪专题STEP 15STEP 15正式叠加和正式叠加和DMODMO此处做浮动基准面校正（见静校正专题）此处做浮动基准面校正（见静校正专题）STEP 16STEP 16叠加模块叠加模块 GASTAK 一般叠加（无参数）一般叠加（无参数）PASTAK 保幅叠加（保幅叠加（PASTAK LW300）TFESTK 时频域零炮检距记录拟合时频域零炮检距记录拟合 PRADMO 保幅时差校正叠加保幅时差校正叠加 DIPDMO 大倾角大倾角DMO WLSTAK 宽线叠加宽线叠加DMO叠加叠加 保持振幅保

47、持振幅DMO 作业例子作业例子:GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容 DIPDMO是基于是基于Kirchhoff积分法的积分法的DMO处理模处理模块。可以一次完成块。可以一次完成NMO校正、校正、DMO校正到叠加输出校正到叠加输出的处理过程。该模块既适用于二维处理，又适用于三的处理过程。该模块既适用于二维处理，又适用于三维处理，既可输出维处理，既可输出DMO后的道集数据，又可输出后的道集数据，又可输出DMO后的叠加数据。后的叠加数据。DMO GRISYS/WS-V8.0 GRISYS/WS-V8.0 新增内容新增内容 TFESTK 在分析地震共中心道集的

48、在分析地震共中心道集的AVO特性，以及不同偏特性，以及不同偏移距信号经过不同的路径的地层吸收的基础上，利用移距信号经过不同的路径的地层吸收的基础上，利用某一时刻的不同偏移距地震信号的不同频率成分，可某一时刻的不同偏移距地震信号的不同频率成分，可以用多项式描述其和偏移距以用多项式描述其和偏移距x 的关系，进而提出了基的关系，进而提出了基于信号时频分析和多项式拟合的时频域零偏移距地震于信号时频分析和多项式拟合的时频域零偏移距地震道拟合技术。同时，利用多道信号估计的信噪比设计道拟合技术。同时，利用多道信号估计的信噪比设计的时频域的时频域Wiener滤波器，克服了常规多项式拟合技滤波器，克服了常规多项

49、式拟合技术压制噪声能力不足的缺陷。术压制噪声能力不足的缺陷。时频域零炮检距记录拟合时频域零炮检距记录拟合 叠后处理叠后处理STEP 17STEP 17提高信噪比（去噪部分）提高信噪比（去噪部分）RNATTE POLFITRAPFIL提高分辨率提高分辨率 RNATTE 根据有效反射波同相轴可予测性根据有效反射波同相轴可予测性,分离出随机噪音，达到提高剖分离出随机噪音，达到提高剖面的信噪比与连续性的目的。面的信噪比与连续性的目的。参数说明参数说明:LWaa为时窗长度，以毫秒为时窗长度，以毫秒(ms)为单位为单位LXbb为窗口宽度，即空间分段道数为窗口宽度，即空间分段道数LOPcc为予测滤波算子长度

50、为予测滤波算子长度(道数道数)，缺省时，缺省时c=7LWOVdd为时窗重叠系数，按百分比提供为时窗重叠系数，按百分比提供LXOVe e为空间分段重叠系数，按百分比提供为空间分段重叠系数，按百分比提供(CMPi-CMPj)=(Tt1PEv1，Tt2PEi2,TtnPEvn)i j为为 CMP控制点号和控制范围控制点号和控制范围,T,PE是响应的时间和是响应的时间和混波百分比混波百分比：JOB P：XINJIANG，L：X001，T：RNATTE：DISKIN DFN：GHSQ2174，HD4=(100-300)，LTR4000，STOP：RNATTE LW512，LX32，LOP7，LWOV10