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

GRISYS 处理系统 二维地震资料处理,处理流程的建立,1、了解地质任务，工区所在位置的地理位置、地表地震条件、构造和地层特征，向用户索取必要的相关资料（近地表信息）。 2、对将要处理的数据进行调查分析，通过试验初步掌握数据的信噪比，分析了解野外施工质量及数据品质在空间上的变化情况。 3、针对数据的实际情况，选定处理系统中对应的处理功能，并对该功能进行试验，最终选定某一功能选件，即针对性措施。 4、结合以往的处理历史档案（新资料首次处理，应结合该地区相临资料的处理历史档案），征求用户的意见，建立处理流程，同时由处理责任人审核批准，完成“处理设计”。 一个好的处理流程应具有：科学性、合理性、针对性。,二维处理常用模块,观测系统定义 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 叠前时间偏移 HCPSTM DMO PRADMO DIPDMO 叠后修饰 RNATTE POLFIT BLUEFI RAPFIL SMDCOM ZPDCON 偏移 WEMIGR 最终显示及合成图 DBAFL AMPEQU DBASI RESECT,数据带解编,STEP 1,GRISYS系统输入输出模块,SEGBSEGD ： DBAIFDBAISDEMULT SEGD ： SEGDIN SEGY(tape+disk) ： SEGYIN、SGYOUT GRISYS(TAPE) ： TAPEIN、TAPOUT GRISYS(DISK) ： DISKIN、DSKOUT、SORTIN DSKTSO、DATINP、DATOUT CGG ： CGGINP、CGGOUT WGC ： CODEIN、CODOUT,SEGDIN SEGD格式磁带解编 为配合海上数据处理，对SEGDIN模块进行了较大规模的改进，增加海上SEGD格式数据读取的功能，并能准确地从原始数据中读取相关数据并存放在对应的GRISYS数据道头字中。,GRISYS/WS-V8.0 改进内容,增加MP参数：为乘MP因子输入标志。 这是一个选择性参数，填此选件后，程序从数据的头块中求出因子的实际值与仪器记录数值相乘。从而，得出近似仪器增益放大前的数值。一般数据输入时不用此参数。仅当做数值分析输入时慎用，缺省为不应用因子。同时，取消TRUE参数。,SEGDIN,GRISYS/WS-V8.0 改进内容,GRISYS/WS-V8.0 改进内容,道头增加内容为： HD80=水、陆检排列的标识。 HD81=海洋资料采集拖缆的缆号。 HD82=分放1，2，3，4排列。 HD83=接收点站号。 HD84=接收线号。 HD85=炮线号。（SEA参数起作用） HD86=炮点站号。（SEA参数起作用） HD88=磁带卷号。 HD89=通道号。 HD90=唯一道号。,SEGDIN,磁盘和磁带数据分析 TPU ITAPUIL ISEGYAD ISGYTAPE 8.0新增,TPU、ITAPUIL 磁带数据分析 两者功能基本相同，TPU为早期的版本，TAPUIL为交互版本。功能主要有： 磁带数据格式自动分析； 磁带数据DUMP； 磁带数据拷贝。 其中磁带数据拷贝又提供多种拷贝方式：带对带、带对盘、盘对带、合并拷贝以及控制拷贝等等。,ISEGYAD 交互SEGY数据显示 本模块的主要功能是在GRISYS系统中，增加对SEG-Y数据磁盘文件的监控、对未知数据文件的格式分析等。 具体功能为： （1）对数据文件打开的同时，对数据文件进行格式分析； （2）对SEG-Y数据体显示； （3）对SEG-Y数据 “C卡”的显示、修改和转存； （4）对SEG-Y数据 “400字节”的文件头块主要内容的显示； （5）对SEG-Y数据 “240字节”的道头主要内容的显示。,交互SEGY磁带数据检查,ISGYTAPE 在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,IMUTEDI 交互地震道切除与编辑 IMUTEDI 交互地震道切除与编辑(新） IDATVIEW 交互地震数据显示 IDATCOMP 交互数据对比显示 IPRTEST 交互参数试验,IMUTEDI,主 要 功 能,显示地震道道集数据及道头字,交互定义切除,交互定义整道充零和开时窗充零,交互定义地震道极性反转,按道头字18交互定义整道编辑库,分析显示指定波形的线性速度,IMUTEDI,道头字显示,IMUTEDI,交互定义切除,imutedi.old,交互道编辑,线性速度分析,imutedi.old,退出程序后,程序将用户定义的所有库以标准GRISYS编码格式存入filename.tool_code中 将切除库的二进制数据存入mulib.bin中，以备下次使用,imutedi.old,随着地震勘探技术和计算机技术的不断发展，地震数据体也变得越来越大。为了适应不断增长得地震数据体对处理系统的需求，此次版本升级对于大数据量数据的处理进行了较大的改进，并新开发了两个针对大数据体显示和编辑的模块。 新版的IMUTEDI； 交互地震数据显示模块IDATVIEW。,GRISYS/WS-V8.0 新增内容,GRISYS/WS-V8.0 新增内容,交互地震道切除与编辑,IMUTEDI（NEW） 新版的IMUTEDI模块除了沿袭了原有模块的主要功能以外，具有使用方便，支持大道集、大数据体的输入和快速定位等特点。 索引文件引入，实现大数据体的快速定位检索； 支持大道集，缺省以2000道为限显示页面； 双滑动杆设计，一个为数据体滑动，另一个为道集内部滑动； 切除库和道编辑库直接读入。,GRISYS/WS-V8.0 新增内容,GRISYS/WS-V8.0 新增内容,IMUTEDI（新）,主菜单及子菜单项,GRISYS/WS-V8.0 新增内容,IMUTEDI（新）,支持多编辑库同时操作，并不互相干扰。,GRISYS/WS-V8.0 新增内容,IMUTEDI（新）,多种道编辑库的文本文件输入和输出，无需二进制文件。,GRISYS/WS-V8.0 新增内容,IMUTEDI（新）,其他各项功能及对话框,GRISYS/WS-V8.0 新增内容,IDATVIEW 较好地适应大数据体输入； 独立的数据索引方式； 数据的快速准确定位； 无极放大功能； 支持双屏独立操作； 自动记忆10屏显示。,交互地震数据显示,GRISYS/WS-V8.0 新增内容,主界面,数据加载,建立索引,数据操作,IDATVIEW,无极放大功能,IDATVIEW,GRISYS/WS-V8.0 新增内容,GRISYS/WS-V8.0 新增内容,IDATVIEW,支持双屏操作； 彼此独立。,GRISYS/WS-V8.0 新增内容,IDATVIEW,数据直接定位,多种显示方式； 自定义显示道头。,显示方式调整,GRISYS/WS-V8.0 新增内容,交互地震数据显示和对比,IDATCOMP 在GRISYS/WS-V7.0中推出了绘图文件对比的交互模块，该模块在地震资料处理中监控处理质量和处理效果起到了很好的作用。但是实际上，地震数据的直接对比更能反映地震数据真实的变化。为此我们此次推出了GRISYS地震数据的对比显示工具。其主要功能是： 双数据输入，独立或同步翻动比较； 显示参数统一，并支持多种显示方式； 鼠标点击切换数据，快捷方便； 支持不同的GRISYS数据格式。,GRISYS/WS-V8.0 新增内容,主界面,数据类型定义,数据显示及对比,显示参数定义,不同类型显示方式,野外静校正 折射波静校正 (专题),STEP 4&5,子波整形 (可选),STEP 6,子波整形,GRISYS处理系统中，子波整形由两个模块来完成，一是整形滤波器设计（SHPFIL），二是 整形反褶积（ SHDECN）。 SHPFIL：同一工区采集的资料，由于使用不同震源或者使用不同记录系统施工，致使信号在工区内不能很好地连续对比，使用子波整形校正可消除子波之间的差异，从而改善处理效果。 首先输入参考数据道以及需要校正的数据道，指定参与计算子波的参考道和需要校正道的范围和时窗大小，近似计算出参考道的子波和需要校正的道的子波，利用维纳滤波器求出这两种子波的差异，即滤波因子。,子波整形,：JOB P：DQ, L：595, T：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,子波整形,把整形滤波器设计应用模块（SHPFIL）求得的整形因子，应用到数据上，达到子波一致的目的 。 ：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，FMT4 ： END,GRISYS常用的反褶积模块,GRISYS常用的反褶积模块,GRISYS常用的反褶积模块,GRISYS常用的反褶积模块,GRISYS常用的反褶积模块,匹配因子,对炸药震源激发单炮进行匹配滤波对比,GRISYS常用的反褶积模块,可控震源激发的原始单炮,匹配滤波后的单炮,GRISYS常用的反褶积模块,子波整形,整形后的处理： 子波整形后的数据，应进行地表一致性反褶积或多道统计子波反褶积处理，以此来进一步完成子波的统一。因为在预测步长较大时，地表一致性反褶积的主要功能是“子波调整”，压缩子波的处理由后续的统计子波反褶积或单道反褶积等模块来完成，即串联反褶积。,叠前去噪(专题),STEP 7,IPRTEST,主 要 功 能,可进行频谱分析、俞氏子波参数选取、视速度拾取等分析,DEGROR 内切滤波法压制面波,ATTGRO 自适应衰减面波,RELNOI 叠前线性干扰滤除,HNOIBP 强能量干扰的分频压制,TSDCON 两步法统计子波反褶积,PRDCON 预测反褶积,参数试验,DEGROR,IPRTEST,参数试验,IPRTEST,DEGROR去面波前后单炮,参数试验,HNOIBP,IPRTEST,参数试验,HNOIBP-压制高能干扰前后单炮,IPRTEST,参数试验,IPRTEST,RELNOI,参数试验,IPRTEST,RELNOI-去线性干扰前后单炮,参数试验,IPRTEST,PRDCON,参数试验,IPRTEST,PRDCON-预测反褶积前后单炮,参数试验,IPRTEST,频谱分析,参数试验,IPRTEST,俞氏子波参数选取,参数试验,GRISYS/WS-V8.0 改进内容,IPRTEST 新增7个模块，分别为DBAFL、PVNOEL、LINOEL、FILTER、SFNPRE、POLFIT、SCALER； 新增模块试验方式与原有模块基本相同； 增加交互拾取拟合参数功能。,振幅补偿,STEP 8,用于分析由于 地表因素 球面扩散 介质吸收 影响而引起的振幅能量衰减情况,AMPANA,振幅分析,AMPANA,: JOB P：PPP，L：L302，T：AMPANA : DBAFL VN01，I(10) : TAPEIN UNIT01，G302，F1，HD2=(2-11)，STOP : AMPANA NX60，NT10，DB，W500-W5000， FL01，(TR1-TR60，I3) : END NXn n为要处理的炮检距个数 NTm m为对每个炮检距所规定的叠加道数 DB 填此标识时，采用对数比例绘图,振幅分析,AMPANA,不同炮检距平均振幅曲线,振幅分析,振幅处理,AMPCOM 振幅补偿,SCAMPC 地表一致性振幅补偿,AMPEQU 振幅动态平衡（E选件）,AMPEQU,空 单道平衡选件 A 整道平衡选件 B 多时窗区域平衡 C 均方根平衡选件 D 瞬时平衡选件 E 单炮平衡选件,振幅处理,AMPEQU 炮间均衡,： DISKIN DFN：GHIFIST1，NTY68（NTY2） ： AMPEQU E，NTA60 ： DSKOUT DFN：GAMPEQU1,E：在用户定义的单炮记录道的规定时窗内，求取平均振幅值，再根据提供的区域平衡级别系数计算出相应记录的平衡系数，将其与实际振幅值相乘，求其平衡后的振幅值。,振幅处理,AMPEQU,炮间能量均衡前 炮间能量均衡后,振幅处理,AMPCOM,主 要 功 能,消除地震船汽枪能量不同所造成的对地震波振幅的影响,取地震道振幅绝对值,消除球面扩散对振幅的影响,消除地层吸收对振幅的影响,将地震道振幅规格化到一给定值,地震道振幅常数加权,振幅处理,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,振幅统计分析,time,A,time,振幅处理,AMPCOM,振幅统计分析,time,A,time,振幅处理,AMPCOM,振幅补偿前,振幅处理,振幅补偿后,AMPCOM,振幅处理,AMPCOM,振幅补偿前后单炮,振幅处理,SCAMPC,去除由激发、接收等因素引起 的接收道的 能量差异,包括 几何扩散和吸收系数补偿 地表一致性振幅道平衡,振幅处理,SCAMPC,: JOB P：PPP，L：LLL，T：SCAMPC : SCAMPC FP1，LIST1，TI300，TA1200， LW400，NT120，T1500， V2500，MIX300，MAX3250 : DSKOUT : PLOTTR SCL75，STP15，TW50 : END 只做几何扩散和吸收系数补偿,振幅处理,SCAMPC,: JOB P：PPP，L：LLL，T：SCAMPC : SCAMPC TI300，TA1200，LW400，TCI400，TCA1300， NT120，T1500，V2500，MIX300，MAX3250， STI61，STA152，I1，SHOT87 : DSKOUT : PLOTTR SCL75，STP15，TW50 : END 做几何扩散和吸收系数补偿及地表一致性道振幅平衡,振幅处理,SCAMPC,地表一致性振幅补偿前单炮,振幅处理,地表一致性振幅补偿后单炮,SCAMPC,振幅处理,一次速度分析剩余静校正（专题）,剩余静校正见静校正专题,STEP 9,速度分析 叠加速度分析 常速扫描速度分析 速度分析数据解释 速度谱解释 常速扫描解释 层位拾取 速度处理 速度换算： 均方根速度换算层速度，等时速 度内插，沿层速度内插 DBAPV插值 交互速度编辑与修改 常速扫描叠加的速度与现有速度场的结合,功能介绍,速度分析所要求的相应处理功能 NMO 叠加 其它辅助性功能 TV 函数加载 TV 函数输出 Mute 函数加载 Mute 函数输出 速度处理和速度场质控图件显示 TV函数剖面展开显示 叠加速度二维彩色剖面显示 速度等值线图显示,速度谱解释窗口,常速扫描解释窗口,7.0 对新增和改进部分的说明 作csum 大道集的选项,便于质量控制; 条件约束自动寻找极值点 速度微调 tv函数展开显示 谱线 + 变面积显示（同批量谱显示） 常速扫描叠加的速度场可选择修改扫描范围内的时间速度对,Svel 简单操作过程,在%状态下或 gscons下启动svel，弹出主界面 按照二级管理（在主界面tool下project manage 项）加载数据，同一地区的数据在load中可同时加载多条测线数据。 主界面File选项，打开本次分析的数据。 主界面Seisfunc功能项中，选VanalyAnalyse弹出作速度谱的参数框 。,此参数框中有三个 参数框必须人工干预 : CMPs for vel. analyse Parameters Muting function definition Parameters Ref velocity curves Parameters 以上三个参数框填写好后按OK鍵，其它参数用缺省值，在速度谱解释窗口中显示出速度谱、大道集和变速叠加段 。,TV&MU输出,主界面File窗口选Svel项, 弹出窗口，在Data type栏中选择保存类型。如: TV库，选TV Pair，在Name对应的参数框填输出文件名后点 OK，存成内部文件。 主界面：ToolTV I/O选项，弹出窗口，点击OUTPUT，选右边栏中的文件名，变黑后，在Grisys TV Filename下边的路径指向上，填写输出文件名，点OK，完成输出功能。此时可在当前目录下找到填写的输出文件。,第一次迭代,第二次迭代,第四次迭代,优势频带的速度分析与剩余静校正的迭代,注意事项,TV库中第一个T切除库中第一个T时间。 输入数据应经过振幅恢复处理的cmp集数据 譜图太窄时，可改变Maxim Offest参数，使其小于实际偏移距。 偶而出现运行不正常时，一要查看数入数据正确性，二建立新工区，重新填参数或进入/grisys/oper/optrun/svel目录中，清除本工区的所有中间文件,GRISYS/WS-V8.0 改进内容,SVEL 8.0改进内容 支持1ms采样数据； 此次修改主要面向大数据体； CMP数/号最大可以支持9999999； 改进谱点定义方式，增加起止点及增量定义方式； 速度谱一次可做512个； 对于没有解释的谱点，生成TV文件时，不再输出。,GRISYS/WS-V8.0 改进内容,滑动条，CMP个数2000时自动分屏；滑动条上的数字为当前页面的第一个CMP号；拖动滑动杆来定位数据。,新的速度谱点定义方式:除了保留原有的定义方式之外,新增了起始点、终止点和增量的定义方式。,交互速度分析,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 LW300 或者 ：GASTAK ： DSKOUT ： FILTER FL01 ： AMPEQU L300 ：PLOTTR SCL100，STP10 ：END,二维叠加作业,地表一致性反褶积 （可选）,STEP 10,地表一致性反褶积（第一步）,地表一致性反褶积可以消除炮点、检波点随地表条件和采集因素的变化对子波带来的影响，可输出稳定和基本一致的地震子波 SURSAN对输入的地震数据按时窗求取其对数功率谱，将对数功率谱按同一分量分别进行累加，对共炮点、共检波点、共偏移距、共中心点等五个分量进行地表一致性相关分析。,SURSAN,参数说明 : NINTXa, Wb-Wc a为炮检距，b、c为时窗起止时间值(ms)，因为只定义一个时窗，其时窗范围尽量大一些，以便包括较多的有效层。 由道号定义时窗 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为相关分析的最大频率，缺省为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：GLF00002 ：END,SURSAN,地表一致性反褶积（第二步）,SURDEC对输入利用SURSAN模块算出的地表一致性反褶积的谱分析结果对输入的地震 ,道先进行地表一致性反褶积算子的求取，然后利用该算子进行地表一致性反褶积。,SURDEC,参数说明 : SPFN: Gxxxxxx 为存放PK分量的对数功率谱的磁盘文件。 GPFN: Gxxxxxx 为存放SK分量的对数功率谱的磁盘文件。 Ll l为反褶积算子的长度， MUn n为切除库号。 SPIKE 脉冲反褶积，为缺省方式。 GAP 预测反褶积。 BPSPK 带通脉冲反褶积。 BPGAP 带通预测反褶积。 TWO 填此选件时，既对pk分量方向做反褶积，同时也对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 ：DISKIN DFN: GDBA02, HD2=(1-300),NTY2 : SURDEC TWO, GAP, L120，GAPLN30，FL6, MU01, FH60，CF20，CFOUT40.0， SPFN：GLF00001， GPFN：GLF00002 : END,SURDEC,原始单炮,地表一致性反褶积后的单炮,二次速度分析剩余静校正,同一次速度分析和剩余静校正,STEP 11,反褶积,STEP 12,反褶积： 改造地震记录子波的一种运算。理想情况为把子波改造成单位脉冲。 基本假设： 子波是小相位的 反射系数序列为白噪,反褶积的基本概念,地震记录的褶积模型 为地层反射系数； 为地震子波。 地震子波包含震源响应、检波器响应，这些响应是时不变的；现在一般检波器都是组合接收，还有组合响应（若炮组合也有组合响应），还有地层吸收，而后二者是时变的。 反褶积的目的是想消除地震子波 影响，从而得到地下反射系数 ，克莱伯特形象地将反褶积比喻为显微镜，使我们能更细致地观察到地下物性变化（纵向）。如有剩余子波也希望是零相位的，因为零相位子波分辨率最高。,反褶积的基本概念,单道反褶积： 脉冲、预测 (TVDCON、PRDCON) 前提均是子波最小相位,并假设反射系数为白噪,但是实际情况地震子波是混合相位,且频带有限,通过实验得出脉冲反褶积压缩子波功能最强,预测其次. 多道反褶积： TSDCON、SURDEC、SCDCON 前提是假设子波具有地表一致性，子波最小相位和反射系数为白噪。 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,频率分析,GRISYS/WS-V8.0 改进内容,SPEANA 频谱分析 改进频谱分析结果的输出方式，使得其结果和交互频谱分析的结果保持一致。,GRISYS/WS-V8.0 改进内容,交互频谱分析和滤波器设计,频率分析,IFSANFL 可实时的对地震记录道形式的各种数据进行Fourier谱分析，包括单道一维谱、多道一维谱及多道的二维谱（即FK谱）； 在一维谱（包括单道和多道）分析的基础上进行滤波器设计和试验，滤波器包括低通、带通、高通、带陷、雷克子波和俞氏子波； 可进行频带扫描分析； 可进行时频谱分析； 可进行多时窗的频谱对比，其中包括不同数据集上相同时窗的频谱对比，相同数据集上不同时窗的频谱对比； 可定义并应用滤波器的作用范围及强度，从而实现时变和空变滤波； 可定义切割滤波； 可对得到的滤波器建立相应的GRISYS滤波库，形成相应批量作业，且用户可随时查看与编辑该作业文件。,频率分析,主界面,数据属性及显示参数,数据显示,IFSANFL,频率分析,IFSANFL,频率分析,选定分析范围,频谱,单道频谱,滤波器设计器,IFSANFL,频率分析,选定分析范围,IFSANFL,设计滤波器,应用效果,效果分析,数据应用,频率分析,IFSANFL,多时窗对比分析,GRISYS/WS-V8.0 改进内容,IFSANFL 一维谱增加分贝显示； 一维谱加盖网格，并可修改显示参数； 改变频率扫描结果显示方式，增大显示面积。,GRISYS/WS-V8.0 改进内容,交互频谱分析,显示结果和方式与SPEANA模块一致,GRISYS/WS-V8.0 改进内容,交互频谱分析,增大了结果显示面积，方便对比观察,分频扫描,频率分析,FBSCAN,最常用也是最直观的原始资料频率特征分析手段。常用于现场处理的野外施工参数对比分析以及室内处理中原始资料的频率分析。 实际上就是对输入数据进行多个不同的滤波处理，并将滤波后的数据显示在同一张剖面上进行对比分析。通过剖面上随着滤波门的变化而产生的数据信噪比的变化来了解资料中有效波和干扰波的频率分布范围和基本特征。 需要配合PLOTTR模块使用。,FBSCAN,：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,HDB17 ：END,频率分析,FBSCAN,10-20HZ,15-30HZ,20-40HZ,频率分析,FBSCAN,30-60HZ,40-80HZ,50-100HZ,60-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:指明输出地震道的包络，缺省输出地震道。 HD22,频率分析,频率分析,频率分析,本模块用两种典型子波俞氏子波和雷克子波对记录进行滤波扫描，处理员可分析扫描结果，用于选择适当的滤波参数和子波反褶积的期望输出。 当输入为脉冲函数时，本模块输出为参数对应的子波；当输入为地震记录时，用俞氏子波或雷克子波对记录做滤波处理，输出为滤波扫描的结果。,频率分析,YUSACN 俞氏及雷克子波扫描,YUSCAN,: 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，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）,频率分析,YUSCAN,RICK（40）,RICK（50）,RICK（60）,原始单炮,频率分析,反褶积前准备工作 尽可能压制各种干扰 振幅补偿 反Q滤波补偿 做好野外静校、折射静校正,反褶积的基本概念,(子波整形 SHPFIL+SHDECN) 脉冲反褶积 TVDCON 预测反褶积 PRDCON 两步法统计子波反褶积 TSDCON 地表一致性反褶积 SURSAN+SUEDEC 叠后提高分辨率：SMDCON、ZPDCON,GRISYS常用的反褶积模块,预测反褶积 脉冲反褶积 前提均是子波最小相位，并假设反射系数为白噪，但是实际情况地震子波是混合相位，且频带有限,通过实验得出脉冲反褶积压缩子波功能最强，预测其次。 地下地震子波是混合相位用脉冲反褶积； 目标为层间多次和海底全程多次，且信噪比低时，用预测反褶积。,GRISYS常用的反褶积模块,预测反褶积,预测反褶积是通过预测滤波来实现的，预测滤波的基本思想是，设计一个滤波器，使得该滤波器具有某种预测能力，通过它对信号的当前值和过去值的滤波，预测未来某个时刻将要出现的信号成分。预测时刻和当前时刻的距离称为预测距离。预测滤波器的设计采用了最小平方准则，即要求预测结果和实际信号间误差的最小平方和最小。 PRDCON是一个预测反褶积模块，既可以用于消除地震记录中的长、 短周期多次波，又可以用来提高地震记录的分辨率。,GRISYS常用的反褶积模块,GRISYS常用的反褶积模块,先在炮集记录上提取统计子波，利用该子波做反褶积处理,消除炮点对子波的影响,第一步期望输出波形为最小相位； 然后在共接收点道集上提取统计子波，利用该子波进行反褶积处理,消除共接收点对子波的影响,第二步期望输出为零相位。 TSDCON不仅能校正、压缩激发震源子波波形，还可以消除地表接收响应。从这个意义上讲，它具有一定的地表一致性反褶积效果。,两步法子波反褶积 TSDCON,TSDCON,GRISYS常用的反褶积模块,：JOB P：MAQIN，L：301，T：OPTIONSG ：DBAMU VN01，HD2，(1-100)=M100X30， M500X570，M1000X1455， : 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 : END,TSDCON,第一步,GRISYS常用的反褶积模块,：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-WN23