简介:面向地震勘探与勘查技术工程专业学习者,这是一份以 Vista5.5 为处理工具的高分辨率地震数据处理技术 PDF 文档,适合正在学习地震资料处理课程、需要掌握常规处理流程的学生及入门工程师使用。资源共 1 个 PDF 文件,压缩包大小约 2.57 MB,内容以实验报告形式呈现,包含新建项目、加载 SEGY 数据、建立观测系统、频谱分析、频率滤波、振幅补偿、FK 滤波、CMP 道集抽取、速度分析及动校正叠加等完整操作步骤,并附有处理流程简图和参数设置说明,方便按图实操。文档特别针对实际数据给出滤波参数(20-25-65-70)与面波切除思路,对理解分辨率提升和去噪处理具有直接参考价值。目前已有 380 人学习下载,适合作为课内实验的补充材料或上机操作速查手册。
1. 高分辨地震数据处理方法:一份 Vista 5.5 全流程实验报告的复现价值
地震资料处理这件事,刚上手时最难的不是某个模块不会点,而是不知道一条能出图的高分辨率处理流程该按什么顺序走、每一步参数凭什么这么定。这份“高分辨地震数据资料处理方法技术Vista5.5.pdf”实验报告,正好把答案写了下来:用 Vista 5.5 把 25 炮 SGY 原始数据从新建项目一路处理到偏移剖面,覆盖频谱分析、带通滤波、球面扩散振幅补偿、FK 滤波、CMP 道集抽取、速度分析、NMO 动校正、叠加、反褶积和叠后偏移。实验采用偏移距 16m、道间距 2m、炮间距 2m、24 道接收、12 次覆盖、8 次叠加的观测系统,是一条单测线 48 道的小排列浅层数据。它适合两类人:一类是做地震勘探课程设计或实验报告、需要照着跑通流程并交付成果的学生;另一类是刚拿到野外 SGY 数据、想快速建立第一条 Vista 5.5 处理流程的初级工程师。
2. 数据准备:新建工程、SGY 加载与观测系统参数
Vista 5.5 里所有分析都建立在数据能被正确定位的前提下。与其后面反复调滤波参数,不如先把观测系统搞对。这一章的每一步都在为后面的模块打底,速度分析、动校正、叠加能跑出什么结果,很大程度上取决于这里参数填得对不对。
2.1 新建二维处理项目与 SGY 数据集加载:文件顺序是第一个坑
File → new project 弹出新建项目对话框,选择二维处理项目类型,指定新项目存储位置后单击 OK。接着在 Data List 窗口的数据集区域,按住不放弹出按钮组合,选择二维数据集,再按住不放选择 SEG-Y 文件格式,在弹出的 Input 对话框中选择数据文件 1.sgy 到 25.sgy,单击打开,数据就加入数据集了。
SEG-Y 是地震数据最通用的交换格式,道头里保存炮号、道号、坐标等信息,但不同采集系统的道头填充习惯不完全一致。Vista 5.5 加载后默认看到一个数据集列表,这里有一个非常容易翻车的细节:25 个文件加载进来后,顺序不一定是 1 到 25,需要手动调整,让文件按预期顺序排列后再点击 OK。原报告单独把这一步写出来,说明它真实发生过、也确实影响后续处理。如果文件顺序不对,后面的 CMP 道集和速度谱会全部错位,而且这种错位在早期不容易发现,等叠加剖面出来才发现道序错乱,返工成本极高。
提示:养成加载后先扫一遍数据集列表的习惯。25 个文件确认顺序只需要十几秒,能省掉后面一整轮返工。
2.2 观测系统参数:AUTO-CALCULATE 之前要理解道距和炮检距
在 Data List 窗口的数据集 1 中,依次点击输入道数、道间距、炮检距等参数,键入观测系统炮点参数和接收点参数,然后点击 AUTO-CALCULATE 自动计算出相关参数,再点 OK。接着点击计算 CMP 及检波器叠加次数和炮检距离,以便写入道头相应位置。
AUTO-CALCULATE 听起来很智能,实际上它只是根据已知参数推导派生参数,并不会校验原始输入是否合理。如果道间距填错,或者炮点、检波点关系填反,算出来的 CMP 号、覆盖次数、炮检距都会跟着错。这份实验数据的观测系统参数如下:
| 参数 | 数值 |
|---|---|
| 偏移距 | 16 m |
| 道间距 | 2 m |
| 炮间距 | 2 m |
| 每次放炮接收道数 | 24 道 |
| 覆盖次数 | 12 次 |
| 叠加次数 | 8 次 |
| 单条测线总道数 | 48 道 |
道间距 2m、炮间距 2m、每炮 24 道接收,配合 16m 偏移距,构成了 12 次覆盖的规则观测系统。这类小排列浅层数据在城市物探、工程勘察里很常见,处理结果的深度尺度在几十米范围,跟油气勘探用的大排列数据完全不是一个量级,速度分析时尤其要注意低速层的存在。
2.3 把观测系统写进道头:验证 CMP 与覆盖次数的正确姿势
执行计算后,在 Data List 中查看 CMP 及检波器叠加次数分布,以此判断观测系统是否正确。确认无误后,按住保存按钮不放,在出现的按钮组合中选择保存观测系统,以便下次使用;再选择将相关信息写入道头。
“写入道头”这一步最容易被当成“保存设置”而忽略,但它决定了后续模块能否正确读到炮检距和 CMP 号。观察覆盖次数曲线是快速验证手段:规则观测系统在测线中间段覆盖次数应该是稳定平台,两端逐渐降为低覆盖;如果中间出现周期性缺口或数值跳变,说明炮点或检波点参数与实际采集不符。覆盖次数图正常后再写道头,写完后还可以再查看一次道头内容,确认 CMP 号、炮检距字段非零再继续。
3. 频谱分析与滤波参数:20-25-65-70 这个梯形滤波器怎么来
处理流程里最常被问的一句话是:滤波参数你凭什么选 20-25-65-70?答案不在软件里,在频谱分析的结果里。这一章就是回答这个问题的。
3.1 振幅谱分析:先看主频范围再定通带
在地震剖面图窗口点击频谱分析按钮,选择分析项目,再点加号执行,弹出振幅谱分析窗口。放大振幅谱观察主频率范围,20-70Hz,尤其 25-45Hz 之间能量最强,为下面滤波做准备。除了振幅谱,Vista 5.5 中还可以分析相位特征、信噪比等。
频谱分析这一步是后面所有滤波参数的地基。原报告的分析章节写的是:有效波频率范围在 20-65Hz 之间,尤其 25-45Hz 之间。频谱窗口里看到的 20-70Hz,是包含干扰在内的整体能量范围;放大振幅谱并仔细对比后,有效信号真正稳定占优的频带是 25-45Hz。这个从“宽范围”到“窄频带”的判断过程,比直接读谱面更重要。频谱分析的结果不只服务于带通滤波,反褶积的预测步长、偏移前的道头检查,都要回看频谱来验证参数有没有改变信号的实际频率成分。
3.2 一维带通滤波:Ormsby 与 TVOrmsby 的适用边界
Vista 5.5 中一维滤波模块很多,常用的是 Ormsby Band-Pass(一维时不变带通滤波)和 TVOrmsby(一维时变带通滤波)。本次使用 Ormsby Band-Pass,滤波参数是梯形滤波器 20-25-65-70。
四个频率点构成一个梯形带通:20Hz 和 70Hz 是切除频率,25Hz 和 65Hz 是通带边界。低于 20Hz 的能量逐步压制,20-25Hz 过渡,25-65Hz 全通,65-70Hz 再过渡,高于 70Hz 切除。用渐变过渡代替硬切,可以避免时间域的振铃效应——矩形滤波器在频率域陡变,对应时间域会产生振荡拖尾,把浅层波形拉出假相位。
| 频段 | 作用 |
|---|---|
| <20 Hz | 缓降切除:压制面波等低频干扰 |
| 20-25 Hz | 过渡带:振幅从 0 缓升至 1 |
| 25-65 Hz | 通带:有效信号完整保留 |
| 65-70 Hz | 过渡带:振幅从 1 缓降至 0 |
| >70 Hz | 缓降切除:压制高频环境噪声 |
TVOrmsby 是时变带通,需要设置随时间变化的滤波窗口点,适合浅层高频丰富、深层高频衰减明显的资料。如果整条剖面的频率特征横向稳定,用 Ormsby 即可;如果浅层能看到明显高频、深层明显低频,再用 TVOrmsby。原报告里专门提了一句:一维滤波可以不做,高低频可以在 FK 滤波中切除。这意味着两种滤波手段存在功能重叠:FK 滤波通过视速度差异分离干扰,能同时兼顾频率和空间方向性,实际效果往往比一维滤波更干净。
3.3 球面扩散振幅补偿:为什么先补能量再做 FK
对资料使用 SphDiv 工作流做球面扩散振幅补偿,可以消除地震波在非完全弹性介质中传播的吸收衰减。这一步放在一维滤波之后、FK 滤波之前,顺序有讲究。
球面扩散补偿解决的是几何扩散和吸收衰减造成的深层能量衰减。地震波传播时波前面面积持续增大,单位面积能量按距离衰减;介质非弹性又使高频成分衰减更快。这两者叠加,深层反射波振幅远小于浅层,不补偿的话,FK 域里深层有效波能量会被面波噪声完全盖住,切面波时等于把深层有效信号一并切掉。这就是 FK 滤波前必须先做振幅补偿的原因——宁可先放大噪声,也要保证深层有效能量进入 FK 域时占据一个可识别的地位。
Vista 5.5 里另一个可选模块是地表一致性振幅补偿(Surface Consistent Application),校正震源能量差异和检波器耦合差异造成的横向振幅不一致。球面扩散补偿解决的是随深度变化的能量衰减,地表一致性解决的是随空间变化的能量不均衡,两者不是替代关系。本次数据为小排列浅层观测,单炮能量相对稳定,优先做球面扩散补偿即可;如果发现相邻炮的能量明显跳动,再考虑补一遍地表一致性。
4. FK 滤波、CMP 道集与常见问题:从切面波到道头陷阱
FK 滤波是这套流程里最能体现“二维”价值的一步。一维带通滤波只能按频率切,面波和有效波的频率范围有重叠时,一维滤波会把两者一起保留或一起切掉。FK 滤波在频率-波数域利用视速度差异,把面波和反射波分开。
4.1 生成 FK 滤波文件:选面波最强的一道做切除基准
打开要做 FK 滤波的地震记录图,选择面波较强的一道。点击选择,适当切掉面波区域,保存滤波文件,作为各道滤波的标准。在结果中查看滤波效果,在 FK 域分析能量,可以继续滤波直至满意。
选面波最强的一道作为切除基准,看起来是经验操作,实际有道理:面波的能量和分布随炮点位置变化,选最强的一道作基准,意味着所有其他道在这个切除窗口下都不会残留明显面波。反过来选一道面波弱的作基准,换到面波强的炮集上就会泄漏大量面波能量。
面波视速度低,在 FK 域表现为靠近低频端的低速扇形区;有效反射波视速度高,能量集中在高速区域。两者频率上有重叠,视速度上分得开,这正是 FK 滤波能处理一维带通解决不了的问题的原因。切除时通常沿着低速扇形区的边界画一道,而不是平切一条线。
4.2 FK Filter 工作流:反复滤波与能量检查
生成滤波文件后,打开 FK Filter 工作流。输入选择球面扩散振幅处理结果(或上一轮 FK 输出),滤波文件选择频谱切除的结果文件,执行后回到 FK 域检查能量分布,若仍有面波残留,继续做下一轮。
这里有个翻车点:有人把 FK 滤波当成一次性操作,执行一轮就往下走。浅层地震资料里,面波可能有多次反射或散射,一轮切除后 FK 域仍能看到低速能量残余。正确的节奏是“滤波→回 FK 域看能量→再滤波”,直到能量分布干净。另一条原则是宁可多切两轮也不要一次切太多:切除窗口过大会同时切掉与面波视速度接近的折射波或反射波,对有效信号造成不可逆损伤。判断标准很简单——每轮切除前后做一次 FK 域对比,看到有效反射能量开始损失就停手。
4.3 抽 CMP 道集:输入输出与 cmp-no 排列
抽取 CMP 道集是为速度分析和动校正做准备。以 FK 滤波结果作为输入,炮点排列方式选择 cmp-no,输入输出各一个按钮即可完成,输出即得到按共中心点排列的道集。
cmp-no 的意思是按 CMP 号排列炮点,把不同炮、不同接收道中属于同一中心点的道组合到一起。只有按这个方式组织的道集才能用于速度谱计算和动校正。抽道集本身不改变数据内容,只重新排列道序,却容易出错:输入文件选错、排列方式选错,抽出来的道集在速度分析时会表现出同相轴不收敛。抽完后可以抽查几个 CMP 道集,看同一同相轴是否大致呈双曲线形态,再进入速度分析。
4.4 踩坑记录:道序、道头与 FK 切除范围的三个问题
把整个流程跑下来,最容易出问题的三个环节集中说:
第一条,加载 SGY 文件后顺序错乱。现象:CMP 道集里炮序跳变,速度谱能量团不连续。原因:SEG-Y 道头顺序与文件名称顺序不一致,加载后默认排序不等于采集顺序。解决:加载后在数据集列表里手动调整文件顺序,确认排列无误后再建观测系统。
第二条,观测系统没写进道头。现象:后续模块正常执行,但输出的 CMP 号或炮检距字段全为 0,速度分析时无法正确排序共中心点道集。原因:AUTO-CALCULATE 只完成了参数计算,没有执行“将相关信息写入道头”。解决:保存观测系统后单独执行写道头操作,用覆盖次数显示功能复核,确认道头字段非零再继续。
第三条,FK 滤波切除窗口过大。现象:滤波输出剖面里浅层反射波变弱,波形明显拉伸、频率变低。原因:低速扇形区被整体平切,误伤视速度接近折射波的能量。解决:缩小切除范围,分两轮迭代滤波,每轮在 FK 域对比切除前后能量差,看到有效反射能量开始损失就停手。
5. 速度分析、NMO 与叠加:Velzone 三件套与参数联动
速度分析是整条流程里最靠经验的一环。它不像滤波那样参数直接可见,而是要在速度谱、常速度叠加剖面上人工判断能量团位置,再用拾取结果去校验叠加效果。Vista 5.5 把这件事拆成 Velzone 工作流、速度拾取、动校正三个环节,一环套一环。
5.1 Velzone 工作流:速度谱、CVS 与偏移距道集三件套
打开工作流 Velzone,输入选择抽道集结果 CMP,并将 Velzone 标记为 for execution。Velzone 是专门用于速度分析的数据输入命令,数据输入控制只允许选择共中心点道集。它内部包含三部分:相似(叠加)速度扫描、常速扫描 CVS 和抽共偏移距道集输出。
在右键菜单 Velzone → Date/Header SELECTIONS 中,选取若干个满覆盖次数的 CMP 道集进行速度分析。速度分析的目标是获得最佳叠加速度,为动校正提供速度场,常用速度谱法和常速度叠加法 CVS。速度谱看能量团集中的位置,CVS 看不同速度叠加剖面上同相轴拉平的程度,两相对照,比只看其中一种可靠得多。
CVS 参数里三个数字分别对应起始扫描速度、终止扫描速度和扫描速度间隔,可以通过 NEW 按钮修改。Sembl 参数设置中,start velocity 为扫描起始速度,end velocity 为扫描终止速度,velocity incr 为速度间隔,与 CVS 保持一致。速度间隔太大,速度谱能量团离散,不好拾取;太小则计算量增大,能量团横向拉长,同样不利判断。浅层低速资料的一般做法是先宽范围粗扫,锁定速度区间后再缩小间隔精扫。
| 参数 | CVS | Sembl |
|---|---|---|
| 起始扫描速度 | NEW 按钮修改 | start velocity |
| 终止扫描速度 | NEW 按钮修改 | end velocity |
| 速度间隔 | NEW 按钮修改 | velocity incr |
| 取值范围 | 结合工区先验速度 | 与 CVS 设置相同 |
Offsrts 中的偏移距参数可以默认,也可以根据排列长度调整。共偏移距道集用于观察不同偏移距下的速度响应,如果速度随偏移距明显变化,说明存在速度各向异性或 AVO 效应,叠加速度分析要更谨慎。
5.2 速度拾取与层速度计算:从交互窗口到深度估算
点击 interactive → velocity interactive → velocity analysis 弹出交互窗口,选择速度分析结果 1、2、3 做速度拾取。拾取前要结合工区先验知识预估速度范围:本次实习区表层为沙土及下部砂砾石,速度较低,大致落在 100-500-900m/s 区间。这个先验信息把速度拾取约束在合理区间内,避免在速度谱上乱抓低质量能量团。
拾取之后定量计算各层平均速度:124.0ms 对应 255.0m/s,208.0ms 对应 315.0m/s,据此确定两层深度分别为 15.8m 和 32.8m,分析可知浅层明显界面可能为潜水面。这两个深度可以用速度乘单程时校验:255.0m/s × 124.0ms / 2 ≈ 15.8m,315.0m/s × 208.0ms / 2 ≈ 32.8m,与报告结果一致。速度拾取的质量直接决定动校正效果:拾取偏快,同相轴校正过量上弯;拾取偏慢,同相轴校正不足仍下弯。这两种误差在叠加剖面上都会以剩余时差的形式糊掉有效同相轴,所以拾取时要多选几个 CMP 道集对比,取能量聚焦最稳的层段。
5.3 动校正 NMO 与叠加:Stretch Mute 设为 100% 的含义
选取动校正与叠加工作流 nmo,输入选择抽取了 CMP 道集的结果,动校正文件是速度拾取产生的 vel 文件。其中 Stretch Mute 为动校正切除量,本次设置为 100%,即不切除。
这里的 100% 很容易让人反向理解。在动校正模块里,Stretch Mute 对应的是动校正拉伸比例阈值:100% 意味着对拉伸畸变不做切除,全部保留。动校正的拉伸效应在浅层、大偏移距处最严重,表现为波形周期被拉长、频率降低。切除一部分拉伸区能改善浅层波形,但会减少有效道数、降低覆盖次数。本次数据单测线道数少,选择保留全部数据换取叠加次数是合理的;如果数据信噪比足够高,适度切除反而能获得更干净的浅层叠加结果。
叠加使用 stack 模块做常规水平叠加,也可以选择自适应加权水平叠加。常规叠加把同一 CMP 道集内经动校正后的各道直接相加,信噪比按道数平方根提升;自适应加权叠加按局部信噪比给各道分配权重,对低信噪比道降权。本次数据覆盖次数只有 8 次,常规叠加已能满足要求,加权叠加的优势不明显。
6. 反褶积与叠后偏移:从 decon 参数到 FK 偏移的验证技巧
叠加完成后,真正把“分辨率”往上推的环节是反褶积和偏移。使用 decon 模块,可以选择预测反褶积、脉冲反褶积和零相位反褶积。本实验在叠加后进行预测反褶积处理:先对叠加剖面做一次频谱分析,确认主频位置,设定预测步长,执行后重新查看频谱,观察主频是否往高频方向移动,以及高频段信噪比是否可接受。
预测步长的设定是反褶积里最需要试验的参数。步长太短,相当于只做子波整形,压缩效果有限;步长过长,会把反射系数序列也当作子波的一部分压掉,噪声被明显放大。常见做法是先用自相关估算子波的主瓣宽度,把步长设在主瓣宽度附近,然后做一组对比试验,取频谱改善和信噪比之间最平衡的数值。零相位反褶积对剖面的保幅特性更好,适合后续要做振幅相对保持解释的资料;预测反褶积在压缩子波的同时还能压制多次波,更适合以构造解释为目标的处理流程。
偏移使用 FK 2-D Migration 对叠后剖面做处理,作用是让倾斜同相轴归位到真实空间位置,绕射波收敛。FK 域偏移计算快,但对速度场的精度有一定容忍度;速度偏低时同相轴校正不足,偏高时过度校正,表现为断裂处的异常波形或同相轴翻转。偏移前务必确认速度场来自合理拾取,不能拿未检查的速度文件硬跑。
还有一个习惯值得养成:每个环节做完后,回看频谱和波形变化。从带通滤波前的 20-70Hz,到反褶积后主频是否真的向 25-45Hz 集中,再到偏移后剖面上同相轴的连续性能否与原始地震记录对应,这些对比比单看一张成果图更能说明问题。从那以后我每次处理叠后数据,都强制把原始频谱、滤波后频谱和反褶积后频谱三张图并排放在一起,确认每一步都有实际改善才往下走,这个习惯在 Vista 5.5 和后来的处理软件里都一样管用。希望帮到你。
本文还有配套的精品资源,点击获取