GeoEast-振幅处理－反褶积-新模板.ppt

振幅处理和反褶积 一 振幅处理二 反褶积 内容 振幅处理 在相对保持振幅的前提下 通过振幅处理 解决由于由激发 接收 地层吸收等因素引起的接收道的能量差异 GeoEast解决振幅问题的流程为炮间能量均衡 解决炮间的能量差异几何扩散补偿 解决纵向的能量差异地表一致性振幅补偿 解决横向的能量差异 GeoEast振幅处理模块 振幅均衡 AmpEqu 几何扩散补偿 AmpCompenst 二维地表一致性振幅补偿 SCAmpCom2D 三维地表一致性振幅分析 SCAmpAna3D 三维地表一致性振幅文件合并 SCAmpSum3D 三维地表一致性振幅补偿 SCAmpCom3D 振幅分析 AmpAna 叠后剩余振幅处理 PostRsAmpPro 反射强度增益 RtStrnthGain 时变振幅加权 TVarScal 地震道振幅自适应加权 AMpClipWeight 峰值振幅增益 PeakAmpGain 内容 振幅均衡 振幅补偿 二维地表一致性振幅补偿 振幅处理 三维地表一致性振幅补偿 振幅分析 内容 振幅均衡 振幅补偿 二维地表一致性振幅补偿 振幅处理 三维地表一致性振幅补偿 振幅分析 在野外原始记录或经叠加处理后的记录上 由于波前扩散和地层吸收作用 往往出现浅 中 深层反射波的能量差异较大或者道与道间的能量不均衡情形 能量不均衡问题 一方面会影响叠加效果 另一方面不便于数据的显示 为此 在资料处理中需要进行能量均衡处理 针对波前扩散引起的能量不均衡问题 地震道振幅动态平衡模块以较大的时空步长将记录振幅平衡到人为设定的某一级别上 以便于地震记录的显示和绘图 所谓动态平衡是指平衡因子在时间和空间上是变化的 因子的大小直接与记录振幅级别相联系 Ampequ振幅均衡 equalizationoptions平衡处理选项 有六个 singletrace单道平衡 缺省值 多时窗 entiretrace整道平衡 单个时窗 multi trace多时窗区域平衡RMS均方根平衡instantaneous瞬时平衡singlesource单炮平衡 Ampequ振幅均衡 starttime第一个时窗或区域的起始时间 缺省值 0 endtime各时窗或区域的结束时间 缺省为0 整道作为一个时窗处理 windowlength各时窗或区域的平衡时窗长度 ms 选择此参数前的 则弹出如图所示的参数表 最多分9个时窗或区域 Ampequ振幅均衡 equalizationlevel平衡级别值 用来调整振幅 缺省为5000 用于单道平衡 整道平衡 多时窗区域平衡和单炮平衡 RMSlevel均方根平衡级别值 均方根平衡和瞬时平衡选件用 缺省为2000 mutetablename切除表文件名 如果做切除就输入切除文件名 可选 applywaterdepth应用水深值yes应用 表示切除表中的时间用水深值定义水深值 存放在道头CMP HIGHV TOP中 用于处理海上资料 no不应用 缺省值 Ampequ振幅均衡 applytaper应用斜坡处理到平衡因子上 只用于均方根平衡和瞬时平衡 yes表示时窗两端平衡系数以斜坡延伸no表示时窗两端平衡系数以常数延伸normalizationcoefficient振幅归一系数 其取值范围是 0 30 applytruemute应用真切除yes真切 将输出数据的切除部分数据充零 no不真切 输出数据的切除inputfactorfilename输入平衡因子文件名 如果填文件名 读取该文件数据直接应用平衡处理 本程序不再计算平衡因子 该文件格式与道记录数据相同 该参数为空时 不输入平衡因子 缺省值 Ampequ振幅均衡 outputfactorfilename输出平衡因子文件名 如果填文件名写入平衡因子到该文件中 该文件格式与道记录数据相同 该参数为空时 不输出平衡因子 缺省值 applyequalization应用道平衡yes在输出道上应用平衡因子 缺省值 no在输出道上不应用平衡因子inverselyapplyequalization反应用道平衡 当输入因子存在时 该参数起作用 yes在输出道上反应用平衡因子no在输出道上不反应用平衡因子 缺省值 list列表打印参数表 Ampequ振幅均衡 使用实例 Ampequ振幅均衡 Ampequ振幅均衡 原始振幅 整炮均衡 Ampequ振幅均衡 整道均衡 原始振幅 Ampequ振幅均衡 原始振幅 单道均衡 Ampequ振幅均衡 原始振幅 多道均衡 绘制单炮振幅对比图可选用多道增益 既可以避免调整增益水平带来的麻烦 又可以保持相对振幅关系 Ampequ振幅均衡 内容 振幅均衡 振幅补偿 二维地表一致性振幅补偿 振幅处理 三维地表一致性振幅补偿 振幅分析 功能概述 在实际资料中 影响反射波振幅的因素主要包括 激发与接收因素 波前扩散 地层吸收 地层层状结构 反射界面的形态 各种干扰波等 这些因素造成反射波在浅 中 深不同位置上波形和能量存在较大差异 记录道之间或炮记录之间的波形和能量也不相同 这种波形和能量上的差异严重影响着反褶积 动校正 静校正和速度分析的精度 因此在叠前进行振幅补偿显得尤为重要 AmpCompenst振幅补偿 消除球面扩散对地震波振幅的影响消除地层吸收对地震波振幅的影响消除地震船汽枪等能量不同所造成的地震波振幅的影响对于海上资料 需要消除球面扩散和地层吸收等造成的影响 AmpCompenst振幅补偿 功能选项 其它功能对地震道的振幅取绝对值 将地震道振幅规格化到某一给定值 对地震道振幅作常数加权处理 输入 输出补偿因子 应用或反应用补偿因子等 AmpCompenst振幅补偿 compensationmethodoptions补偿方法选件 有四个选项 sphericaldivergence球面扩散补偿 缺省值 stratumabsorption地层吸收补偿energy db 能量分贝补偿marine海上补偿 AmpCompenst振幅补偿 球面扩散补偿参数卡 1 补偿类型 球面补偿 吸收补偿 2 补偿速度 3 速度内插方式 建议输出补偿因子 以便将来需要的时可进行反补偿 AmpCompenst振幅补偿 地层吸收补偿参数卡 stratumabsorption指数函数的指数 值缺省为1 0 AmpCompenst振幅补偿 energycoefficient能量分贝补偿系数 能量分贝补偿选项参数 缺省为1 energyconstant能量常数补偿值 能量分贝补偿选项参数 缺省为1divergencevalue几何扩散补偿值 海上补偿选项参数 缺省为5 normalizationcoefficient振幅归一化系数 处理后的最大振幅值 缺省为0 不做归一化处理 weightingcoefficient常数加权的权系数 缺省值为1 0 inputfactorfilename输入补偿因子文件名 如果填文件名 读取该文件数据直接应用补偿 本程序不再计算补偿因子 该文件格式与道记录数据相同 该参数为空时 不输入补偿因子 缺省值 outputfactorfilename输出补偿因子文件名 如果填文件名 补偿因子写到该文件中 该文件格式与道记录数据相同 该参数为空时 不输出补偿因子 缺省值 applycompensationfactor应用补偿因子yes应用 缺省值 no不应用inverselyapplycompensationfactor反应用补偿因子yes应用no不应用 缺省值 AmpCompenst振幅补偿 内容 振幅均衡 振幅补偿 二维地表一致性振幅补偿 振幅处理 三维地表一致性振幅补偿 振幅分析 二维地表一致性振幅补偿 SCAmpCom2D二维地表一致性振幅补偿本模块是一个叠前振幅补偿程序 其主要目的是去除由激发 接收等因素引起的接收道的能量差异 它包含两方面的内容 其一 几何扩散和吸收系数补偿 其二 地表一致性振幅道平衡 常规叠前振幅处理 在精度和保真度方面还存在着一定的局限性 主要表现在 叠前单道道平衡不是一种相对保持振幅的处理 由于难以给出正确的吸收系数和速度函数 几何扩散和吸收项不能精确补偿 为了克服它们的不足之处 本模块采用 利用地表一致性统计的方法 使得叠前振幅保持相对道平衡 求取地表一致性低频振幅能量曲线 进行叠前几何发散和吸收项的精确补偿 options处理旗标 有两个选项all做几何扩散和吸收系数补偿及地表一致性振幅道平衡 缺省为all nobalance只作几何扩散和吸收系数补偿 二维地表一致性振幅补偿 内容 振幅均衡 振幅补偿 二维地表一致性振幅补偿 振幅处理 三维地表一致性振幅补偿 振幅分析 功能概述 该模块主要用于消除由于激发 接收等因素引起的地震记录在横向上的能量差异 假设在近地表内 炮点处近地表的影响只表现为由该炮点所造成的能量不同 而与信号的记录位置无关 同样地 检波点处近地表的影响也只表现为由该检波点所造成的能量不同 而与信号的炮点位置无关 根据这一假设进行地表一致性振幅补偿 实现地震数据的横向能量均衡 实现步骤 第一步 SCAmpAna3D产生振幅矩阵库文件 第二步 SCAmpSum3D在处理三维资料时 首先在线束上运行振幅分析模块 SCAmpAna3D 得到多个振幅分析文件 然后运行该模块将多个振幅分析文件合并分解 第三步 SCAmpCom3D 振幅补偿 调用上一步合并分解后的补偿因子 完成地表一致性振幅补偿 对于二维资料也可以使用本方法进行二维地表一致性振幅补偿 只需使用第一步和第三步 三维地表一致性振幅补偿 地表一致性振幅分析 三维地表一致性振幅补偿 地表一致性振幅合并 分解 三维地表一致性振幅补偿 地表一致性振幅补偿 三维地表一致性振幅补偿 GeoEast振幅补偿流程 叠前去噪 炮间均衡 可选 几何扩散补偿 二维可选 地表一致性振幅补偿 三维地表一致性振幅补偿 原始单炮 球面扩散补偿后单炮 地表一致性补偿后单炮 三维地表一致性振幅补偿 内容 振幅均衡 振幅补偿 二维地表一致性振幅补偿 振幅处理 三维地表一致性振幅补偿 振幅分析 振幅分析AmpAna 1 功能概述 该模块对振幅进行定量分析 它计算随炮检距变化的平均振幅值和随时间变化的平均振幅值 用于分析由于介质吸收和球面扩散影响而引起的振幅能量衰减情况 用户通过不同的参数可获得不同的CGM图件 用户通过显示方式 displaymode 参数可决定用对数比例绘图还是常规方式 通过输出图形 outputgraphs 参数可获得所需CGM图件 本模块可输出四种CGM图件 a 综合振幅图 不做均化处理 b 随时间变化的振幅图 均化处理 c 随偏移距变化的平均振幅图 均化处理 d 综合振幅图 均化处理 振幅分析AmpAna 2 1输入数据炮集数据 建议先进行切除处理 炮的范围和炮内地震道的选择由输入模块控制 2 2输出数据CGM格式的图形文件 3 使用说明outputCGMfilename 输出的CGM文件名 不填 缺省 时 放在绘图队列 系统定义的CGM文件目录下 文件名为 作业名 cgm numberofoffsets 要处理的炮检距个数 缺省为120 maximumtraceseachoffset 每个炮检距的最大叠加道数 缺省时 叠加同一炮检距的全部输入道 振幅分析AmpAna isplaymode 显示方式decibel采用对数比例绘图方式normal常规绘图方式outputgraphs 输出图形summarygraph只输出综合振幅图 不做均化处理 threegraphs输出三种图形a 随时间变化的振幅图 均化处理 b 随偏移距变化的平均振幅图 均化处理 c 综合振幅图 均化处理 振幅分析AmpAna 4 作业例子 输出test4 cgm文件到当前目录 输出综合振幅图 显示方式为分贝 不做均化处理 5 限制及注意事项 键入gecgm 启动GeoCGM可查看输出的CGM图件 键入gecgm 启动GeoCGM 通过File下的RasterPlot菜单项光栅CGM图件 再通过绘图命令输出光栅后的文件到绘图仪 振幅分析AmpAna 振幅曲线显示 一 振幅处理二 反褶积 内容 内容 单道反褶积 反褶积 反褶积 ShapOperDsgn整形滤波器设计ShapFilt整形滤波SRDecon两步法统计子波反褶积SCSpecDecon3D三维地表一致性反褶积LogSpectrum对数谱 复赛谱计算SCSpecDecom3D谱分解SCSpecAna3D三维地表一致性谱分析SCDecon3D三维地表一致性反褶积PredictDecon预测反褶积 单道 MultiPrDecon多道预测反褶积SpikDecon脉冲反褶积ZeroPhasDecon零相位反褶积 内容 单道反褶积 反褶积 ShapOperDsgn整形滤波器设计ShapFilt整形滤波同一工区采集的数据 由于采集条件的变化或处理条件不同 都会存在差异 这就需要提取整形因子对数据进行整形滤波来消除数据之间的差异 从而改善处理效果 该模块通过输入需要校正的数据道和指定参考数据道 在一定的时间和空间范围内 计算出参考道的平均响应作为滤波器的期望输出 同理计算出需要校正道的平均响应作为滤波器的输入 再利用维纳滤波求出这个滤波器 即整形滤波因子 该因子保存在由用户指定名字的滤波器数据表中 供整形滤波模块使用 子波整形 ShapOperDsgn整形滤波器设计输入数据需要进行整形校正的叠前或叠后地震数据 参考地震数据 注 参考数据由用户指定磁盘文件名 在模块中自行打开 输出数据整形因子和输入数据 整形因子和经过整形的输入数据 子波整形 子波整形 ShapFilt整形滤波通过读取滤波器数据表内的整形因子 对输入地震数据进行滤波 使被处理数据与参考数据的差异最小 从而改善处理结果 达到子波整形的目的 输入数据叠前或叠后地震数据 输出数据经过整形校正的叠前或叠后地震数据 子波整形 整形后的处理 子波整形后的数据 应进行地表一致性反褶积或多道统计子波反褶积处理 以此来进一步完成子波的统一 因为在预测距离较大时 地表一致性反褶积的主要功能是 子波调整 压缩子波的处理由后续的统计子波反褶积或单道反褶积等模块来完成 即串联反褶积 子波整形 内容 单道反褶积 反褶积 地表一致性反褶积分为 一步完成 SCDecon3D三维地表一致性反褶积两步完成 SCSpecAna3D三维地表一致性谱分析SCSpecDecon3D三维地表一致性反褶积三步完成 LogSpectrum对数谱 复赛谱计算模块SCSpecDecom3D谱分解SCSpecDecon3D三维地表一致性反褶积 地表一致性反褶积 SCDecon3D三维地表一致性反褶积该模块是一个多道地表一致性反褶积模块 它考虑以下四个域求解 炮点域 检波点域 CMP域 共炮检距域 在地表一致性假设下 每一道与炮点 检波点 CMP点 炮检距算子进行褶积 CMP算子可选 因为它代表地下构造 可以多个计算时窗 应用时窗程序自动计算 一步地表一致性反褶积 一步地表一致性反褶积 SCSpecAna3D三维地表一致性谱分析SCSpecDecon3D三维地表一致性反褶积 两步地表一致性反褶积 SCSpecDecon3D SCSpecAna3D三维地表一致性反褶积功能概述利用SCSpecANA3D模块计算的自相关函数 组成自相关函数矩阵 解托布里兹矩阵方程组就可求得相应的反滤波因子 对于每一个道 均可得到相应的炮点反褶积因子 接收点反褶积因子 中点位置反褶积因子以及与炮检距有关的反褶积因子 对数据道相继应用这些反因子 就完成地表一致性反褶积 实现步骤 第一步 SCSpecAna3D三维地表一致性谱分析第二步 SCSpecDecon3D三维地表一致性反褶积 两步地表一致性反褶积 SCSpecAna3D三维地表一致性谱分析功能概述SCSpecAna3D是地表一致性反褶积处理中的第一个模块 该模块对输入的地震数据按时窗求取其对数功率谱 将对数功率谱按同一分量分别进行累加 然后对其进行地表一致性相关分析 最终由SCSpecDecon3D模块求出共炮点 共检波点分量的反褶积算子并应用在地震道 两步地表一致性反褶积 参数界面 两步地表一致性反褶积 参数说明windowdefinitiontype时窗定义方法 有两个选项 tracenumber按道号offset按炮检距parametersbytracenumber当按道号定义时窗时 使用此参数 分别输入道号与时窗起止时间值 ms 因为只定义一个时窗 其时窗范围尽量大一些 以便包括较多的有效层 parametersbyoffset当按炮检距定义时窗时 使用此参数 分别输入炮检距与时窗起止时间值 ms 因为只定义一个时窗 其时窗范围尽量大一些 以便包括较多的有效层 如图3所示 横向分别是炮检距 开始 结束时间 最多可输入10对 两步地表一致性反褶积 headerwordforfirstcomponent相关分析第一分量所选用的道头shotnumber炮号reveivernumber检波点号CMPnumberCMP号offset炮检距tracenumber道号numberofgathersinfirstheaderword第一分量道集个数headerwordforsecondcomponentnumberofgathersinsecondheaderwordautocorrelationfunctionlength定义自相关计算时窗长度 用毫秒表示 spectrumfilename输出谱文件名 maximumfrequencyforcorrelation相关分析的最大频率 缺省为80 0Hz 两步地表一致性反褶积 使用实例 炮集作为相关分析的第一个分量 共检波点集为相关分析的第二个分量 由道号定义时窗 求出的对数功率谱记入表文件scspecana3d 两步地表一致性反褶积 SCSpecDecon3D三维地表一致性反褶积功能概述本模块利用SCSpecANA3D模块计算的自相关函数 组成自相关函数矩阵 解托布里兹矩阵方程组就可求得相应的反滤波因子 对于每一个道 均可得到相应的炮点反褶积因子 接收点反褶积因子 中点位置反褶积因子以及与炮检距有关的反褶积因子 对数据道相继应用这些反因子 就完成地表一致性反褶积 参数界面 两步地表一致性反褶积 powerspectrumfilename存放对数功率谱分量的磁盘文件名 muteparameterfilename切除数据表文件名 operatorlength反褶积算子长度 以毫秒表示 最大为SCSpecAna3D模块中相关分析的长度 numberofcomponents分量数量的选项 one对第一个分量做反褶积two对第一分量和第二分量做反褶积processalltracesoptions是否对SCSpecAna3D分析中未参加运算的地震道做反褶积的选项 zerophaseoperatoroptions滤波算子零相位化deconvolutiontype反褶积类型spike脉冲反褶积prediction预测反褶积 两步地表一致性反褶积 predictivedistance ms 预测距离 缺省为1个采样间隔 该参数只有反褶积类型选择为预测反褶积 prediction 时才使用 band passfilteroptions带通滤波选项 lowcutfrequency带通低截频 缺省为10Hz highcutfrequency带通高截频 缺省为90Hz whitenoisecoefficient白噪百分比 缺省为10 percentageofadditionalpowerspectrumoutsidebandpass加入通放带外的功率谱的百分比 缺省为50 此时兼顾了信噪比与分辨率 该值变小时分辨率提高 变大时信噪比提高 两步地表一致性反褶积 使用实例 两步地表一致性反褶积 三步地表一致性反褶积 三维地表一致性反褶积LogSpectrum对数谱 复赛谱 计算SCSpecDecom3D三维地表一致性谱分解SCSpecDecon3D三维地表一致性反褶积因子计算及应用 三步地表一致性反褶积 功能概述 LogSpectrum是地表一致性反褶积处理中的第一个模块 该模块对输入的地震数据按时窗求取其对数功率谱 然后由三维地表一致性谱分解模块 SCSpecDecom3D 将对数功率谱按同一分量分别进行累加 然后对其进行地表一致性相关分析 最终由三维地表一致性反褶积模块 SCSpecDecon3D 求出共炮点 共检波点分量的反褶积算子并应用在地震道 LogSpectrum模块 按炮检距定义时窗 输出复赛谱 此输出由geodiskout模块输出 第一步 对数谱 复赛谱 计算 三步地表一致性反褶积 SCSpecDecom3D模块 把炮点作为第一个分量 检波点作为第二个分量 对LogSpectrum输出的复赛谱进行谱分解 求出各分量记入表文件scspecdecom3d中 放在 项目 工区 下的datatable目录下 第二步 三维地表一致性谱分解 三步地表一致性反褶积 SCSpecDecon3D模块 调用复赛谱分解文件scspecdecom3d 做地表一致性预测反褶积 算子长为120ms 预测距离为24ms 第三步 三维地表一致性反褶积因子计算及应用 三步地表一致性反褶积 注意事项 1 如果三维地表一致性谱分解模块 SCSpecDecom3D 所要分解的数据量比较大 则建议选用复赛谱选件 complexcepstrum 中的yes选件 当多个数据输入的时候 所输入数据的谱类型及参数必须相同 即在对数谱作业中的参数应该完全相同 3 SCSpecDecom3D模块虽然说明可以将对数谱分解为四个分量 但对应的反褶积模块SCSpecDecon3D最多仅能求和应用两个分量的反褶积算子 故建议实际处理时本模块仅选择两个分量进行分解较为合适 建议分选的第一项为炮点项 第二项为检波点项 4 为了节省时间 推荐选择炮点和检波点分量 三步地表一致性反褶积 SRDecon两步法统计子波反褶积本模块中两步法是指 首先在炮集记录上提取统计子波 极小相位 根据不同的期望子波做反褶积处理 然后在共接收点道集上提取统计子波 极小相位或零相位 根据不同的期望子波进行反褶积处理 本模块也可仅在炮集记录上提取子波 进行反褶积处理 即完成两步法中的第一步 本模块不仅能校正压缩激发震源子波波形 还可以消除地表响应 从这个意义上讲 它具有一定的地表一致性反褶积效果 统计子波反褶积 SRDecon两步法统计子波反褶积根据不同的地震记录 用户可定义不同的期望输出 可以为雷克子波 俞氏子波 带通因子 高通因子 低通因子 对于 在炮集记录上提取子波做零相位反褶积 选件 期望输出为小相位 其它选件时期望输出为零相位 本模块共有四个功能选项 functionoptions功能选件 统计子波反褶积 SRDecon功能选项一在炮集记录上提取子波做反褶积该选件完成两步法中的第一步 即只在炮集记录上提取子波 然后进行反褶积 子波为极小相位 期望输出也是极小相位 用户在第一个流程中使用该选件后 在其下一个流程中必须用 在接收点集记录上提取子波做反褶积 选项继续处理 结果才是正确的 统计子波反褶积 SRDecon功能选项二在接收点集记录上提取子波做反褶积该选项在完成 在炮集记录上提取子波做反褶积 的基础上 进一步在接收点集记录上提取子波 然后进行反褶积处理 子波为小相位 期望输出零相位 统计子波反褶积 SRDecon功能选项三在炮集记录上提取子波做零相位反褶积本选件只在炮集记录上完成统计子波反褶积 子波为小相位 期望输出零相位 用户可在下一个流程中选用 在接收点集记录上提取子波做零相位反褶积 选项继续处理 也可以直接以本选件结果作为最终输出 统计子波反褶积 SRDecon功能选项四在接收点集记录上提取子波做零相位反褶积该选件在完成 在炮集记录上提取子波做零相位反褶积 选项的基础上 进一步在共接收点道集上提取统计子波 期望输出为零相位 注意事项a 对输出道的最大振幅值进行修改 当记录中的强能量干扰的道数大于总道数的一半时 本方法将无法识别噪声 故不宜使用本模块 统计子波反褶积 统计子波反褶积 wavelettype子波类型 期望输出有五种子波类型可选择 Rickerwavelet雷克子波 缺省 若选此项 期望输出的地震子波主频 Hz 参数必填 Yuwavelet俞氏子波 若选此项 F1 F2必填 表示俞氏子波的起止频率 highpassoperator期望输出为高通滤波 若选此项 因子F1 F2必填 表示高通滤波算子过渡带的频率 lowpassoperator期望输出为低通滤波因子 若选此项 F1 F2必填 表示低通滤波算子过渡带的频率 bandpassoperator期望输出为带通滤波因子 若选此项 F1 F2 F3 F4必填 表示带通滤波算子过渡带和频宽的频率 统计子波反褶积 内容 单道反褶积 反褶积 PredictDecon预测反褶积 单道反褶积 PredictDecon预测反褶积 功能概述假设反射系数序列为白噪序列 地震子波为最小相位 在以上两个假设的前提下在时域进行反褶积 该模块既可以用于消除地震记录中的长 短周期多次波 又可以用来提高地震记录的分辨率 该模块既可输出反褶积后的数据又可输出反褶积因子 本模块的参数界面 参数说明 operatorlength算子长度