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

1、GRISYS 处理系统二维地震资料处理处理流程的建立 1、了解地质任务，工区所在位置的地理位置、地表地震条件、构造和地层特征，向用户索取必要的相关资料（近地表信息）。 2、对将要处理的数据进行调查分析，通过试验初步掌握数据的信噪比，分析了解野外施工质量及数据品质在空间上的变化情况。 3、针对数据的实际情况，选定处理系统中对应的处理功能，并对该功能进行试验，最终选定某一功能选件，即针对性措施。 4、结合以往的处理历史档案（新资料首次处理，应结合该地区相临资料的处理历史档案），征求用户的意见，建立处理流程，同时由处理责任人审核批准，完成“处理设计”。 一个好的处理流程应具有：科学性、合理性、针对性

2、。二维处理常用模块 观测系统定义 DBAGE PLOTGE HLABEL GRISPS IDFGEDB 野外静校正应用 DBAST HIFIST 道分选 GASORT SORTIN DATINP DATOUT 二维叠前去噪 LINOEL DEGROR HNOIBP COHNMV（2D）RELNOI RNAB2D 振幅补偿 AMPANA AMPCOM SCAMPC 地表一致性反褶积 SURSAN SURDEC SCDCON 速度分析 SVEL 二维剩余静校正 AUREST IRESTAT(交互) 反褶积 TSDCON PRDCON TVDCON MTDCON 叠加 PASTAK GASTAK 叠

3、前时间偏移 HCPSTM DMO PRADMO DIPDMO 叠后修饰 RNATTE POLFIT BLUEFI RAPFIL SMDCOM ZPDCON 偏移 WEMIGR 最终显示及合成图 DBAFL AMPEQU DBASI RESECT数据带解编STEP 1GRISYS系统输入输出模块SEGBSEGD ： DBAIFDBAISDEMULTSEGD ： SEGDINSEGY(tape+disk) ： SEGYIN、SGYOUTGRISYS(TAPE) ： TAPEIN、TAPOUTGRISYS(DISK) ： DISKIN、DSKOUT、SORTIN DSKTSO、DATINP、DATO

4、UTCGG ： CGGINP、CGGOUTWGC ： CODEIN、CODOUTSEGDIN SEGD格式磁带解编 为配合海上数据处理，对SEGDIN模块进行了较大规模的改进，增加海上SEGD格式数据读取的功能，并能准确地从原始数据中读取相关数据并存放在对应的GRISYS数据道头字中。 GRISYS/WS-V8.0 改进内容 增加MP参数：为乘MP因子输入标志。 这是一个选择性参数，填此选件后，程序从数据的头块中求出因子的实际值与仪器记录数值相乘。从而，得出近似仪器增益放大前的数值。一般数据输入时不用此参数。仅当做数值分析输入时慎用，缺省为不应用因子。同时，取消TRUE参数。 SEGDIN G

5、RISYS/WS-V8.0 改进内容 GRISYS/WS-V8.0 改进内容 道头增加内容为： HD80=水、陆检排列的标识。 HD81=海洋资料采集拖缆的缆号。 HD82=分放1，2，3，4排列。 HD83=接收点站号。 HD84=接收线号。 HD85=炮线号。（SEA参数起作用） HD86=炮点站号。（SEA参数起作用） HD88=磁带卷号。 HD89=通道号。 HD90=唯一道号。 SEGDIN磁盘和磁带数据分析TPUITAPUILISEGYADISGYTAPE 8.0新增TPU、ITAPUIL 磁带数据分析 两者功能基本相同，TPU为早期的版本，TAPUIL为交互版本。功能主要有： 磁

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

7、 在GRISYS/WS-V7.0中推出了SEGY磁盘数据检查的交互模块。此次版本升级我们将SEGY磁带数据检查补充进来，使得SEGY数据检查功能更加完善。其具有以下几个特点： 数据DUMP； 3200和400字节头块检查； 显示参数可调； 道头数据即时显示； 任意道头数据曲线绘制。 GRISYS/WS-V8.0 新增内容 GRISYS/WS-V8.0 新增内容主界面及磁带机选择磁带数据显示即时道头信息读取交互SEGY磁带数据检查 GRISYS/WS-V8.0 新增内容交互SEGY磁带数据检查 其他功能3200 头块400 头块数据DUMP观测系统定义（专题）STEP 2道编辑STEP 3 IM

8、UTEDI 交互地震道切除与编辑 IMUTEDI 交互地震道切除与编辑(新） IDATVIEW 交互地震数据显示 IDATCOMP 交互数据对比显示 IPRTEST 交互参数试验IMUTEDI主 要 功 能显示地震道道集数据及道头字交互定义切除交互定义整道充零和开时窗充零交互定义地震道极性反转按道头字18交互定义整道编辑库分析显示指定波形的线性速度IMUTEDI道头字显示IMUTEDI交互定义切除imutedi.old交互道编辑线性速度分析imutedi.old退出程序后,程序将用户定义的所有库以标准GRISYS编码格式存入filename.tool_code中将切除库的二进制数据存入muli

9、b.bin中，以备下次使用imutedi.old 随着地震勘探技术和计算机技术的不断发展，地震数据体也变得越来越大。为了适应不断增长得地震数据体对处理系统的需求，此次版本升级对于大数据量数据的处理进行了较大的改进，并新开发了两个针对大数据体显示和编辑的模块。 新版的IMUTEDI； 交互地震数据显示模块IDATVIEW。 GRISYS/WS-V8.0 新增内容 GRISYS/WS-V8.0 新增内容交互地震道切除与编辑 IMUTEDI（NEW） 新版的IMUTEDI模块除了沿袭了原有模块的主要功能以外，具有使用方便，支持大道集、大数据体的输入和快速定位等特点。 索引文件引入，实现大数据体的快速

10、定位检索； 支持大道集，缺省以2000道为限显示页面； 双滑动杆设计，一个为数据体滑动，另一个为道集内部滑动； 切除库和道编辑库直接读入。 GRISYS/WS-V8.0 新增内容启动IMUTEDI数据加载建立索引文件数据操作保存道编辑参数库具体操作流程退出 GRISYS/WS-V8.0 新增内容IMUTEDI（新）主菜单及子菜单项 GRISYS/WS-V8.0 新增内容IMUTEDI（新）支持多编辑库同时操作，并不互相干扰。 GRISYS/WS-V8.0 新增内容IMUTEDI（新）多种道编辑库的文本文件输入和输出，无需二进制文件。 GRISYS/WS-V8.0 新增内容IMUTEDI（新）其

11、他各项功能及对话框 GRISYS/WS-V8.0 新增内容IDATVIEW 较好地适应大数据体输入； 独立的数据索引方式； 数据的快速准确定位； 无极放大功能； 支持双屏独立操作； 自动记忆10屏显示。交互地震数据显示启动IDATVIEW数据加载建立索引文件数据操作退出具体操作流程 GRISYS/WS-V8.0 新增内容主界面数据加载建立索引数据操作IDATVIEW无极放大功能IDATVIEW GRISYS/WS-V8.0 新增内容 GRISYS/WS-V8.0 新增内容IDATVIEW支持双屏操作；彼此独立。 GRISYS/WS-V8.0 新增内容IDATVIEW数据直接定位多种显示方式；自

12、定义显示道头。显示方式调整 GRISYS/WS-V8.0 新增内容交互地震数据显示和对比 IDATCOMP 在GRISYS/WS-V7.0中推出了绘图文件对比的交互模块，该模块在地震资料处理中监控处理质量和处理效果起到了很好的作用。但是实际上，地震数据的直接对比更能反映地震数据真实的变化。为此我们此次推出了GRISYS地震数据的对比显示工具。其主要功能是： 双数据输入，独立或同步翻动比较； 显示参数统一，并支持多种显示方式； 鼠标点击切换数据，快捷方便； 支持不同的GRISYS数据格式。 GRISYS/WS-V8.0 新增内容主界面数据类型定义数据显示及对比显示参数定义不同类型显示方式野外静校

13、正折射波静校正(专题)STEP 4&5子波整形(可选)STEP 6子波整形 GRISYS处理系统中，子波整形由两个模块来完成，一是整形滤波器设计（SHPFIL），二是 整形反褶积（ SHDECN）。SHPFIL：同一工区采集的资料，由于使用不同震源或者使用不同记录系统施工，致使信号在工区内不能很好地连续对比，使用子波整形校正可消除子波之间的差异，从而改善处理效果。 首先输入参考数据道以及需要校正的数据道，指定参与计算子波的参考道和需要校正道的范围和时窗大小，近似计算出参考道的子波和需要校正的道的子波，利用维纳滤波器求出这两种子波的差异，即滤波因子。子波整形：JOB P：DQ, L：595, T

14、：332：DISKIN DFN: GZHY0001， (参考道) DFN: GZHY0000，NTY2 (要校正的道)：SHPFIL MU02, CMP97, TZ4950, CMP100, TZ4950, LW380, ACMP96, ATZ4950, ACMP98, ATZ4950, ALW380, OUT1，SHFN:SHAB1, PSTK, SHOT, MAX4000, MIN100, AMAX4000, AMIN100, L100, TAP50, NOIS1, LIST, FOLD120, HD20, ALL, OUT0, PLOT END子波整形 把整形滤波器设计应用模块（SHPF

15、IL）求得的整形因子，应用到数据上，达到子波一致的目的 。 ：JOB P：DQ, L：595, T：332 ：DBAMU VN01, HD4, (1-500)=M100X50, M1000X3000, TAP10 ：DISKIN DFN: GZHY0000, LTR3000, 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

16、，FMT4 ： END参考炮子波整形前单炮GRISYS常用的反褶积模块参考炮频谱子波整形前频谱GRISYS常用的反褶积模块参考炮子波整形后单炮GRISYS常用的反褶积模块子波整形后频谱参考炮频谱GRISYS常用的反褶积模块原始CMP道集845Hz滤波后 存在相位差子波整形校正两种检波器GRISYS常用的反褶积模块匹配因子对炸药震源激发单炮进行匹配滤波对比GRISYS常用的反褶积模块可控震源激发的原始单炮匹配滤波后的单炮GRISYS常用的反褶积模块子波整形前叠加剖面（INLINE）子波整形后叠加剖面（INLINE）子波整形前叠加剖面（INLINE）子波整形整形后的处理： 子波整形后的数据，应进行

17、地表一致性反褶积或多道统计子波反褶积处理，以此来进一步完成子波的统一。因为在预测步长较大时，地表一致性反褶积的主要功能是“子波调整”，压缩子波的处理由后续的统计子波反褶积或单道反褶积等模块来完成，即串联反褶积。叠前去噪(专题)STEP 7IPRTEST主 要 功 能可进行频谱分析、俞氏子波参数选取、视速度拾取等分析DEGROR 内切滤波法压制面波ATTGRO 自适应衰减面波RELNOI 叠前线性干扰滤除HNOIBP 强能量干扰的分频压制TSDCON 两步法统计子波反褶积PRDCON 预测反褶积参数试验DEGRORIPRTEST参数试验IPRTESTDEGROR去面波前后单炮参数试验HNOIBP

18、IPRTEST参数试验HNOIBP-压制高能干扰前后单炮IPRTEST参数试验IPRTESTRELNOI参数试验IPRTESTRELNOI-去线性干扰前后单炮参数试验IPRTESTPRDCON参数试验IPRTESTPRDCON-预测反褶积前后单炮参数试验IPRTEST频谱分析参数试验IPRTEST俞氏子波参数选取参数试验 GRISYS/WS-V8.0 改进内容IPRTEST新增7个模块，分别为DBAFL、PVNOEL、LINOEL、FILTER、SFNPRE、POLFIT、SCALER；新增模块试验方式与原有模块基本相同；增加交互拾取拟合参数功能。振幅补偿STEP 8用于分析由于 地表因素 球

19、面扩散 介质吸收 影响而引起的振幅能量衰减情况AMPANA振幅分析AMPANA: JOB P：PPP，L：L302，T：AMPANA: DBAFL VN01，I(10): TAPEIN UNIT01，G302，F1，HD2=(2-11)，STOP: AMPANANX60，NT10，DB，W500-W5000， FL01，(TR1-TR60，I3) : END NXn n为要处理的炮检距个数 NTm m为对每个炮检距所规定的叠加道数 DB 填此标识时，采用对数比例绘图 振幅分析AMPANA不同炮检距平均振幅曲线振幅分析振幅处理 AMPCOM 振幅补偿 SCAMPC 地表一致性振幅补偿 AMPEQ

20、U 振幅动态平衡（E选件）AMPEQU空单道平衡选件A整道平衡选件B多时窗区域平衡C均方根平衡选件D瞬时平衡选件E单炮平衡选件 振幅处理AMPEQU 炮间均衡： DISKIN DFN：GHIFIST1，NTY68（NTY2）： AMPEQU E，NTA60： DSKOUT DFN：GAMPEQU1 E：在用户定义的单炮记录道的规定时窗内，求取平均振幅值，再根据提供的区域平衡级别系数计算出相应记录的平衡系数，将其与实际振幅值相乘，求其平衡后的振幅值。 振幅处理AMPEQU炮间能量均衡前 炮间能量均衡后振幅处理AMPCOM主 要 功 能消除地震船汽枪能量不同所造成的对地震波振幅的影响 取地震道振幅

21、绝对值 消除球面扩散对振幅的影响消除地层吸收对振幅的影响将地震道振幅规格化到一给定值 地震道振幅常数加权 振幅处理AMPCOM： AMPCOM T： AMPCOM T1.5： AMPCOM T2，ABS： AMPCOM ABC，C2： AMPCOM ABC，A3，B4，C2，NR3000： AMPCOM ABS： AMPCOM SDC，TV02，WGT1000 振幅处理AMPCOM振幅统计分析timeAtime振幅处理AMPCOM振幅统计分析timeAtime振幅处理AMPCOM振幅补偿前振幅处理振幅补偿后AMPCOM振幅处理AMPCOM振幅补偿前后单炮振幅处理SCAMPC去除由激发、接收等因

22、素引起 的接收道的能量差异,包括 几何扩散和吸收系数补偿 地表一致性振幅道平衡 振幅处理SCAMPC: JOB P：PPP，L：LLL，T：SCAMPC: SCAMPCFP1，LIST1，TI300，TA1200， LW400，NT120，T1500， V2500，MIX300，MAX3250: DSKOUT .: PLOTTRSCL75，STP15，TW50: END 只做几何扩散和吸收系数补偿 振幅处理SCAMPC: JOB P：PPP，L：LLL，T：SCAMPC: SCAMPCTI300，TA1200，LW400，TCI400，TCA1300， NT120，T1500，V2500，MI

23、X300，MAX3250， STI61，STA152，I1，SHOT87 : DSKOUT .: PLOTTRSCL75，STP15，TW50: END 做几何扩散和吸收系数补偿及地表一致性道振幅平衡 振幅处理SCAMPC地表一致性振幅补偿前单炮振幅处理地表一致性振幅补偿后单炮SCAMPC振幅处理一次速度分析剩余静校正（专题）剩余静校正见静校正专题STEP 9 速度分析叠加速度分析 常速扫描速度分析速度分析数据解释速度谱解释常速扫描解释 层位拾取 速度处理速度换算： 均方根速度换算层速度，等时速 度内插，沿层速度内插DBAPV插值交互速度编辑与修改常速扫描叠加的速度与现有速度场的结合功能介绍速

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

25、，同一地区的数据在load中可同时加载多条测线数据。主界面File选项，打开本次分析的数据。主界面Seisfunc功能项中，选VanalyAnalyse弹出作速度谱的参数框 。此参数框中有三个 参数框必须人工干预 : CMPs for vel. analyse Parameters Muting function definition Parameters Ref velocity curves Parameters 以上三个参数框填写好后按OK鍵，其它参数用缺省值，在速度谱解释窗口中显示出速度谱、大道集和变速叠加段 。TV&MU输出 主界面File窗口选Svel项, 弹出窗口，在Data t

26、ype栏中选择保存类型。如: TV库，选TVPair，在Name对应的参数框填输出文件名后点 OK，存成内部文件。 主界面：ToolTV I/O选项，弹出窗口，点击OUTPUT，选右边栏中的文件名，变黑后，在Grisys TV Filename下边的路径指向上，填写输出文件名，点OK，完成输出功能。此时可在当前目录下找到填写的输出文件。第一次迭代第二次迭代第四次迭代优势频带的速度分析与剩余静校正的迭代注意事项 TV库中第一个T切除库中第一个T时间。 输入数据应经过振幅恢复处理的cmp集数据 譜图太窄时，可改变Maxim Offest参数，使其小于实际偏移距。 偶而出现运行不正常时，一要查看数入

27、数据正确性，二建立新工区，重新填参数或进入/grisys/oper/optrun/svel目录中，清除本工区的所有中间文件 GRISYS/WS-V8.0 改进内容SVEL 8.0改进内容支持1ms采样数据；此次修改主要面向大数据体；CMP数/号最大可以支持9999999；改进谱点定义方式，增加起止点及增量定义方式；速度谱一次可做512个；对于没有解释的谱点，生成TV文件时，不再输出。 GRISYS/WS-V8.0 改进内容滑动条，CMP个数2000时自动分屏；滑动条上的数字为当前页面的第一个CMP号；拖动滑动杆来定位数据。新的速度谱点定义方式:除了保留原有的定义方式之外,新增了起始点、终止点和

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

29、W300 或者 ：GASTAK： DSKOUT： FILTER FL01： AMPEQU L300：PLOTTR SCL100，STP10：END二维叠加作业地表一致性反褶积（可选）STEP 10地表一致性反褶积（第一步） 地表一致性反褶积可以消除炮点、检波点随地表条件和采集因素的变化对子波带来的影响，可输出稳定和基本一致的地震子波 SURSAN对输入的地震数据按时窗求取其对数功率谱，将对数功率谱按同一分量分别进行累加，对共炮点、共检波点、共偏移距、共中心点等五个分量进行地表一致性相关分析。SURSAN参数说明 :NINTXa, Wb-Wc a为炮检距，b、c为时窗起止时间值(ms)，因为只定

30、义一个时窗，其时窗范围尽量大一些，以便包括较多的有效层。 由道号定义时窗 TRd, We-Wf d为道号，e、f为时窗起止时间值(ms)。 PK(g, h) g相关分析时第一分量所选用的道头字，h为g所指定的道集在输入数据中的道集数。 SK(i,j) i为相关分析时第二个分量所选用的道头字，j为i所指定的道 集在输入数据中的道集数。 通常情况下取g=2，i=18进行相关分析，所得效果较为理想。 LWl l为相关长度(ms)。 SPFN: Gxxxxxx 为存放PK分量的对数功率谱的磁盘文件 GPFN: Gxxxxxx 为存放SK分量的对数功率谱的磁盘文件MAXFRf f为相关分析的最大频率，缺

31、省为80.0Hz。 NITERSn n为迭代次数，缺省为10。 NTm m为输入道数，缺省为全部输入道SURSAN使用实例 :第一步：JOB P: LNC, L: SGS, T: DEGROR ：DISKIN DFN: GDBA02, HD2=(1-300),NTY2： FILTER FL01： AMPEQU L300：MUTING BGIN,MU01 ：SURSAN X100，W400W2500， X3050，W370024500，NT30000， MAXFR120，LW200，NITERS20， PK(2, 300)，SK(18, 130)， SPFN：GLF00001，GPFN：GLF0

32、0002 ：END SURSAN地表一致性反褶积（第二步） SURDEC对输入利用SURSAN模块算出的地表一致性反褶积的谱分析结果对输入的地震 ,道先进行地表一致性反褶积算子的求取，然后利用该算子进行地表一致性反褶积。 SURDEC参数说明 :SPFN: Gxxxxxx 为存放PK分量的对数功率谱的磁盘文件。 GPFN: Gxxxxxx 为存放SK分量的对数功率谱的磁盘文件。 Ll l为反褶积算子的长度，MUn n为切除库号。 SPIKE 脉冲反褶积，为缺省方式。 GAP 预测反褶积。 BPSPK 带通脉冲反褶积。 BPGAP 带通预测反褶积。 TWO 填此选件时，既对pk分量方向做反褶积，

33、同时也对sk分量方向做反褶积,缺省为仅做pk分量的反褶积。GAPLNg g为预测步长，缺省为1个采样间隔。 FLa a为带通低截频，缺省为5.0(Hz)。 FHb b为带通高截频，缺省为60.0(Hz)。 a, b的范围应大于有效波的频率范围。 CFc c为白噪百分比，缺省为0.1。 CFOUTd d为加入的通放带外的功率谱的百分比，d缺省为50。此时兼顾了信 噪比与分辨率，d变小时分辨率提高，d变大的信噪比提高SURDEC 使用实例 :第二步 预测反褶积，算子长为120，预测步长为30 （算子的长度不能大于SURSAN中的LW参数）：JOB P: LNC, L: SGS, T: SURDEC

34、：DISKIN DFN: GDBA02, HD2=(1-300),NTY2: SURDEC TWO, GAP, L120，GAPLN30，FL6, MU01, FH60，CF20，CFOUT40.0， SPFN：GLF00001， GPFN：GLF00002 : ENDSURDEC原始单炮地表一致性反褶积后的单炮二次速度分析剩余静校正同一次速度分析和剩余静校正STEP 11反褶积STEP 12反褶积：改造地震记录子波的一种运算。理想情况为把子波改造成单位脉冲。基本假设： 子波是小相位的 反射系数序列为白噪反褶积的基本概念 地震记录的褶积模型 为地层反射系数； 为地震子波。 地震子波包含震源响应

35、、检波器响应，这些响应是时不变的；现在一般检波器都是组合接收，还有组合响应（若炮组合也有组合响应），还有地层吸收，而后二者是时变的。 反褶积的目的是想消除地震子波 影响，从而得到地下反射系数 ，克莱伯特形象地将反褶积比喻为显微镜，使我们能更细致地观察到地下物性变化（纵向）。如有剩余子波也希望是零相位的，因为零相位子波分辨率最高。反褶积的基本概念单道反褶积： 脉冲、预测 (TVDCON、PRDCON) 前提均是子波最小相位,并假设反射系数为白噪,但是实际情况地震子波是混合相位,且频带有限,通过实验得出脉冲反褶积压缩子波功能最强,预测其次.多道反褶积： TSDCON、SURDEC、SCDCON 前

36、提是假设子波具有地表一致性，子波最小相位和反射系数为白噪。 SURDEC、SCDCON以子波调整为主，TSDCON以压缩子波，提高分辨率为主。 反褶积的基本概念频率分析 SPEANA 频谱分析 IFANSFL 交互频谱分析 FBSCAN 频率域滤波与频带扫描 计算并显示地震信号的振幅谱和相位谱：JOB P：NKK，L：LNK1，T：SPEANA ：DISKIN DFN：G0000321， HD4=(150-152) ：SPEANA HD4， (150-152， W500-W1000)， AMP2,LTWO ：END AMP2，以增益方式绘振幅谱SPEANA频率分析SPEANA频率分析 GRIS

37、YS/WS-V8.0 改进内容SPEANA 频谱分析 改进频谱分析结果的输出方式，使得其结果和交互频谱分析的结果保持一致。 GRISYS/WS-V8.0 改进内容交互频谱分析和滤波器设计频率分析IFSANFL 可实时的对地震记录道形式的各种数据进行Fourier谱分析，包括单道一维谱、多道一维谱及多道的二维谱（即FK谱）； 在一维谱（包括单道和多道）分析的基础上进行滤波器设计和试验，滤波器包括低通、带通、高通、带陷、雷克子波和俞氏子波； 可进行频带扫描分析；可进行时频谱分析；可进行多时窗的频谱对比，其中包括不同数据集上相同时窗的频谱对比，相同数据集上不同时窗的频谱对比；可定义并应用滤波器的作用

38、范围及强度，从而实现时变和空变滤波；可定义切割滤波；可对得到的滤波器建立相应的GRISYS滤波库，形成相应批量作业，且用户可随时查看与编辑该作业文件。 频率分析主界面数据属性及显示参数数据显示IFSANFL频率分析IFSANFL频率分析选定分析范围频谱单道频谱滤波器设计器IFSANFL频率分析选定分析范围IFSANFL设计滤波器应用效果效果分析数据应用频率分析IFSANFL多时窗对比分析 GRISYS/WS-V8.0 改进内容IFSANFL一维谱增加分贝显示；一维谱加盖网格，并可修改显示参数；改变频率扫描结果显示方式，增大显示面积。 GRISYS/WS-V8.0 改进内容交互频谱分析显示结果和

39、方式与SPEANA模块一致 GRISYS/WS-V8.0 改进内容交互频谱分析增大了结果显示面积，方便对比观察分频扫描频率分析FBSCAN 最常用也是最直观的原始资料频率特征分析手段。常用于现场处理的野外施工参数对比分析以及室内处理中原始资料的频率分析。 实际上就是对输入数据进行多个不同的滤波处理，并将滤波后的数据显示在同一张剖面上进行对比分析。通过剖面上随着滤波门的变化而产生的数据信噪比的变化来了解资料中有效波和干扰波的频率分布范围和基本特征。 需要配合PLOTTR模块使用。FBSCAN：JOB P：NKK，L：LNK1，T：SPEANA ：DISKIN DFN：G0000321， HD2=

40、(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,HDB17：END 频率分析FBSCAN10-20HZ15-30HZ20-40HZ频率分析FBSCAN30-60HZ40-80HZ50-10

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

42、道的包络，缺省输出地震道。HD22频率分析地震道 时频分析结果 时频分析结果ENVELP频率分析地震道 时频分析结果频谱分析频率分析 本模块用两种典型子波俞氏子波和雷克子波对记录进行滤波扫描，处理员可分析扫描结果，用于选择适当的滤波参数和子波反褶积的期望输出。 当输入为脉冲函数时，本模块输出为参数对应的子波；当输入为地震记录时，用俞氏子波或雷克子波对记录做滤波处理，输出为滤波扫描的结果。频率分析YUSACN 俞氏及雷克子波扫描YUSCAN: JOB P:2D,L:001,T:T: DATGEN SI1,LTR4000,NT4,RD,D1000,Z4100,1,ZG1000: YUSCAN YU

43、（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，200）*: YUSCAN YU（10，60），YU（25，150），YU（35，300）：AMPEQU L300：PLOTTR SCL0100，STP100，LR，HEAD，GAP3HD22输入脉冲函数输入地震道频率分析YUSCAN原始单炮RICK（20）RICK（30）RICK（40）频

44、率分析YUSCANRICK（40）RICK（50）RICK（60）原始单炮频率分析 反褶积前准备工作 尽可能压制各种干扰 振幅补偿 反Q滤波补偿 做好野外静校、折射静校正反褶积的基本概念(子波整形 SHPFIL+SHDECN)脉冲反褶积 TVDCON预测反褶积 PRDCON两步法统计子波反褶积 TSDCON地表一致性反褶积 SURSAN+SUEDEC叠后提高分辨率：SMDCON、ZPDCON GRISYS常用的反褶积模块 预测反褶积 脉冲反褶积 前提均是子波最小相位，并假设反射系数为白噪，但是实际情况地震子波是混合相位，且频带有限,通过实验得出脉冲反褶积压缩子波功能最强，预测其次。 地下地震子

45、波是混合相位用脉冲反褶积； 目标为层间多次和海底全程多次，且信噪比低时，用预测反褶积。 GRISYS常用的反褶积模块预测反褶积 预测反褶积是通过预测滤波来实现的，预测滤波的基本思想是，设计一个滤波器，使得该滤波器具有某种预测能力，通过它对信号的当前值和过去值的滤波，预测未来某个时刻将要出现的信号成分。预测时刻和当前时刻的距离称为预测距离。预测滤波器的设计采用了最小平方准则，即要求预测结果和实际信号间误差的最小平方和最小。 PRDCON是一个预测反褶积模块，既可以用于消除地震记录中的长、 短周期多次波，又可以用来提高地震记录的分辨率。 GRISYS常用的反褶积模块GRISYS常用的反褶积模块原始

46、单炮脉冲反褶积预测反褶积 先在炮集记录上提取统计子波，利用该子波做反褶积处理,消除炮点对子波的影响,第一步期望输出波形为最小相位； 然后在共接收点道集上提取统计子波，利用该子波进行反褶积处理,消除共接收点对子波的影响,第二步期望输出为零相位。 TSDCON不仅能校正、压缩激发震源子波波形，还可以消除地表接收响应。从这个意义上讲，它具有一定的地表一致性反褶积效果。两步法子波反褶积 TSDCONTSDCONGRISYS常用的反褶积模块：JOB P：MAQIN，L：301，T：OPTIONSG ：DBAMU VN01，HD2，(1-100)=M100X30， M500X570，M1000X1455，

47、: DISKI N DFN:GG05854，NTY2（NTY68），HD2=(1-350)， : TSDCON SG，NT96，XMIN25，XMAX1455 , MU01， (WN400-WN1100，WF600-WF1300， XN25-XF1455， T100-T1500，L120，WC1010，BR10-BR65， EP0)， (WN1200-WN2300，WF1400-WF2500， XN25-XF1455， T1700-T5000，L120，WC1015，BR10-BR55， EP0) : DSKOUT DF:GTSDCON1，FMT4: ENDTSDCON第一步GRISYS常用的

48、反褶积模块：JOB P：MAQIN，L：301，T：OPTIONSG ：DBAMU VN01，HD2，(1-100)=M100X30， M500X570，M1000X1455，: DISKIN DFN:GTSDCON2，NTY18: TSDCON ZR，NT48，XMIN25，XMAX1455, MU01， (WN400-WN1100，WF600-WF1300， XN25-XF1455， T100-T1500，L120，WC1010，BR10-BR65， EP0)， (WN1200-WN2300，WF1400-WF2500， XN25-XF1455， T1700-T5000，L120，WC10

49、15，BR10-BR55， EP0) : DSKOUT DF:GTSDCON4: ENDTSDCON第二步GRISYS常用的反褶积模块原始单炮脉冲反褶积预测反褶积GRISYS常用的反褶积模块统计子波反褶积 反褶积参数选取（以两步法为例） 1）计算时窗,应用时窗 2）算子长度 1.52倍子波长度 3）白噪系数 10051020 4）期望输出的子波相位 零相位、小相位 5）期望输出的子波类型 带通子波 脉冲 FC85Hz 雷克子波 FCpickup启动拾取状态并弹出层控制界面 GRISYS/WS-V8.0 改进内容添加层拾取控制点填写拟合参数保存参数交互拟合参数拾取 GRISYS/WS-V8.0 改进内容交互拟合参数拾取参数拾取完成后，传递给参数试验对话框 GRISYS/WS-V8.0 改进内容交互拟合参数拾取拟合处理后对比RAPFIL 径向预测滤波 是一个多道径向预测滤波模块，通过时移求出最大相关值所对