估算及抑制地震数据中异常干扰幅值的方法

摘要

本发明揭示了一种识别地震资料中异常于扰幅值的方法。对构成所述地震资料的预定数目的轨迹，在预定的时间期间内导出一个实测幅值因子，采用高斯-赛得尔矩阵转换法也可计算出对应于相同的时间期间及轨迹数目的无干扰辐值因子。计算出实测幅值因子与无干扰幅值因子之比；从而确定出一个较正因子，将其与一个预先定义的临界值相比较，来确定地震记录中异常干扰幅值的存在。用户可确定是否需要抑制或校正受到异常干扰幅值影响的轨迹。

主权项

1、一种用于识别地震资料中异常干扰幅值数据的方法，包括如下步骤：(a)对地震资料中的每一轨迹，计算出一个作为时间函数的实测幅值因子；(b)对地震资料中一条轨迹的一部分，计算出一个作为时间函数的无干扰幅值因子；(c)通过无干扰幅值因子和实测幅值因子计算出一个校正比率；(d)将校正比率与一个预定值相比较，从而确定一条轨迹中是否存在一异常干扰幅值，以及在采用传统的地震数据处理方法之前是否需要校正异常值或消去轨迹；其中对地震资料的每一轨迹计算作为时间函数的一个实测幅值因子时采用下面的表达式：Aijh=Si＊Rj＊Gkh＊Olh＊Cc＊N其中：Aijh是实测幅值；Si是与源位置i相关的幅值变化分量；Rj是与探测器位置j相关的幅值变化分量；Gkh是在时限h内与公共深度点位置k(也叫中点)相关的幅值变化分量；Olh是在时限h内与偏移1相关的幅值变化分量；Cc是与记录通道c相关的幅值变化分量；N是干扰量；在计算无干扰幅值因子的步骤中进一步包括根据下面的表达式叠加计算矩阵转换：logAijh=logSi+logRj+logGkh+logOlh+logCc+logN其中计算矩阵转换的方法通过高斯-赛得尔(Gauss-Seidel)矩阵转换计算法完成，其中计算校正比率的步骤通过下面的表达式完成：其中上标n和临界值是由用户定义的。