4. 原始记录噪声分析实战

本节把前面学到的噪声识别知识落到一份真实可操作的诊断流程上。目标不是“看一眼记录说噪声大”，而是对一段原始地震记录给出可复核的结论：噪声由哪几类成分构成、各自占据什么频带、沿什么方向传播、在哪些时段和道集上最严重、有效反射还能不能保住，最终形成一份能指导后续去噪参数选择的噪声诊断报告。下面以一个陆上可控震源单炮记录为例展开，示例中的频率、视速度、信噪比数值均为教学假设值，用于说明判读方法，不代表任何实际工区实测结果。

一、先明确分析对象与参照约定

分析前必须固定三件事，否则后续结论无法对齐。第一是坐标与方向约定：本文统一以“炮点位于记录左侧、检波器沿测线向右排列”为参照，视速度取绝对值并注明正负号含义——正号表示视速度方向与排列方向一致（下行/前向），负号表示反向传播（上行/后向）。第二是时间基准：以激发时刻为零时刻，记录长度取6 s，采样间隔2 ms。第三是道集范围：选取近偏移距（0～800 m）、中偏移距（800～2000 m）、远偏移距（2000～3500 m）三个区段分别统计，避免用全排列平均值掩盖空间差异。

原始记录进入分析前只做必要的格式转换与道头检查，不做任何滤波、增益均衡或道编辑，否则会改变噪声的真实频谱与视速度特征，使诊断失真。可以在显示时使用AGC仅用于肉眼观察，但所有定量计算必须回到未均衡数据上完成。

二、噪声类型判别：从记录形态到成因假设

判别的基本方法是“形态—频带—视速度—空间分布”四要素交叉验证。单看波形容易误判，例如低频线性干扰与面波在单道上都表现为低频强振幅，必须结合视速度和频散特征区分。

• 面波（地滚波）：表现为从炮点向大号方向延伸的扇形强能量，低频（示例主频5～18 Hz）、视速度低（示例200～600 m/s）、有明显频散（低频快、高频慢），在近偏移距最强，随偏移距增大逐渐衰减并与反射分离。
• 随机噪声：无固定视速度，道间不相关，频谱宽而平缓，在整张记录上均匀分布，表现为背景“毛刺”，是信噪比的本底限制。
• 50 Hz工频及谐波：单频窄带，道间同相，视速度趋于无穷大（近似垂直同相轴），在频谱上表现为尖锐单峰。
• 线性干扰（如折射波、侧向能量）：视速度高（示例1500～3000 m/s），同相轴平直，频带与有效波重叠，是最难用单纯滤波去除的一类。
• 异常道：单道或少数道出现削顶、死道、脉冲尖峰，需在道编辑阶段单独标记，不参与统计。

实操中建议先在显示上圈出可疑区域，再回到该区域的原始道做频谱与相关分析，形成“假设—验证”的闭环，而不是凭经验直接下结论。

三、频谱分析：分时、分偏移距、分道集统计

频谱分析不能只对整炮做一次FFT，那样会把面波、反射、随机噪声混在一起。正确做法是分窗口统计：在时间上取浅层强噪声段（示例0.3～1.2 s）与中深层反射段（示例1.5～3.5 s），在空间上按近、中、远偏移距分段，对每段做振幅谱并叠加平均，得到“时间—偏移距—频率”的三维能量分布认识。

判读要点有三：一是看有效反射的主频带（示例15～45 Hz）与面波频带（示例5～18 Hz）的重叠区，重叠越宽，单纯频带滤波的伤害越大；二是看工频单峰的幅度与稳定性，判断是否需要陷波；三是看高频端（示例>60 Hz）能量是否被随机噪声主导，决定反褶积和增益策略。若面波与反射在15～18 Hz存在明显重叠，就应优先考虑f-k或视速度滤波，而不是陡截止的带通。

四、视速度估算：从同相轴斜率到f-k能量团

视速度估算有两种互补手段。其一是时距曲线法：在记录上拾取某一道干扰同相轴的到达时间t1、t2，对应偏移距x1、x2，则视速度 v = (x2−x1)/(t2−t1)。示例中面波在x=500 m处t≈1.0 s、x=1000 m处t≈1.9 s，估算视速度约555 m/s，与面波典型值吻合。其二是f-k谱法：对分窗口数据做二维傅里叶变换，在f-k平面上观察能量团的斜率，斜率即视速度倒数，能量团所在频率范围即干扰频带。f-k谱能同时给出频率与视速度的联合分布，是区分面波与线性干扰最直观的工具。

估算时要注意空间假频：道间距过大时，高视速度干扰会在f-k谱上发生折叠，出现虚假的低视速度能量团。示例道间距25 m、采样2 ms，对应Nyquist波数约0.02 周/m，视速度高于约2500 m/s的干扰需警惕假频，必要时抽道或加密分析。

五、信噪比评估：分频带、分时段定量

信噪比不能只给一个全频带数字。推荐做法是：在反射目的层时窗内，选取信噪比相对高的道段估计信号能量S，在无反射的噪声时窗内估计噪声能量N，按频带分别计算SNR=10·lg(S/N)。示例结果：全频带SNR约6 dB，15～45 Hz反射主频带约9 dB，5～18 Hz面波带仅约−2 dB（噪声能量大于信号），>60 Hz高频带约1 dB。这组数字直接告诉我们去噪重点在低频面波带，而高频端主要受随机噪声限制，应靠叠加而非滤波解决。

此外应统计面波影响的最大偏移距（示例约1800 m以内最强）和工频干扰的稳定性，作为后续f-k滤波锥形参数与陷波器设计的依据。

六、形成噪声诊断报告

诊断报告应做到“结论可追溯、参数可执行”。建议包含以下要素：数据基本信息（炮号、道数、道间距、采样率、记录长度）；分析窗口划分（时间窗与偏移距分段）；各类噪声的频带、视速度、空间分布与强度排序；分频带信噪比表；异常道清单；以及针对每类噪声推荐的去噪手段与关键参数初值。例如面波主导时推荐f-k滤波，锥形参数依据实测视速度555 m/s与频带5～18 Hz设定；工频干扰推荐50 Hz及谐波陷波；随机噪声则说明主要靠叠加与反褶积改善，不宜过度滤波。

完成上述步骤后，一份原始记录噪声诊断报告就具备了指导后续去噪的能力：它明确了“压什么、保什么、在哪个频带和视速度区间动手”，也为去噪前后的效果对比提供了量化基线。