数据记录中可以先分选出需要处理的地震数据，如按照关键字（炮号 序号或者炮检距大小等）分出需要处理的地震数据。然后可以读取地震道中设置的随机 数，将所述地震道按照所述道头中随机数的数值大小进行升序或者降序排列，形成随机排 列后的待处理道集数据，如图3所示，图3是上述图2中地震道道序随机重排后的局部地震 数据记录图。这样，可以将采集的原始地震数据记录经过地震道道序随机重新排列后形成 新的地震数据记录，实现相干异常振幅转化成随机异常振幅。
[0039] 当然，所述对地震道的随机重排不限于上述实施例中所述的利用随机函数生成随 机数的方法，还可以包括本领域技术人员不需要创造性劳动而采取其他对地震道数据重排 的方法。
[0040] 将地震数据记录中地震道的道序进行随机排列，然后组成待处理道集数据。
[0041] S2 :计算所述待处理道集数据的采样点在预定时空窗内的参考振幅。
[0042] 本实施例中可以设置所述随机重新排列后的待处理道集数据包括M个地震道，每 个地震道中有N个采样点，地震数据记录时长为T0至TL。
[0043] 在判断地震道中的采样点振幅是否为异常振幅时，可以为所述采样点设置参考振 幅。所述参考振幅可以为所述采样点在一定空间内采用一定方式计算得到的值。本实施例 中可以为所述采样点设置一个预定时空窗，所述预定时空窗具体的可以以所述采样点为中 心的时长T*地震道数N的空间范围。例如T可以取值为300毫秒、N取值为5,可以表示当 前采样点的时空窗的范围为当前采样点的采样时刻前后150毫秒、当采样点所在地震道左 右各两个地震道的范围。具体的可以表示为以所述当前采样点为中心的300毫秒*5个地 震道的地震数据范围。
[0044] 然后可以计算采样点在所述预定时空范围内的参考振幅值。所述的参考振幅值可 以根据数据处理需求定义为不同的参数和计算方法。一般的，在本发明异常振幅的数据处 理中，可以以所述预定时空范围内所有参考点振幅绝对值平滑滤波后的值作为所述参考振 幅。具体的本发明所述计算所述待处理道集数据的采样点在预定时空窗内的参考振幅可以 包括：
[0045] S21 :计算所述待处理道集数据的采样点在预定时空窗内所有采样点的平均值，以 所述平均值作为所述参考点在所述预定时空窗内的参考振幅。
[0046] 例如在300毫秒*3个地震道的预定时空窗内，假设当前采样点P的采样时刻T为 200毫秒，处于所述待处理道集数据中第2个地震道中。那么在本实施例中可以计算所述待 处理道集数据中第1个地震道至第3个地震道从采样时刻TU为50毫秒到采样时刻TD为 350毫秒范围内所有采样点的振幅的绝对值的平均值，以所述平均值作为所述采样点P在 所述预定时空窗内的参考振幅。
[0047] 上述所述的方法可以在所述预定时空窗内地震道数较少的应用场景中使用。对于 所述预定时空窗内地震道道数较多的应用场景，本发明可以提供另一种实施方式实现参考 振幅的计算。具体的，所述计算所述待处理道集数据的采样点在预定时空窗内的参考振幅 可以包括：
[0048] S221 :判断所述待处理道集数据采样点的预定时空窗内的地震道道数是否大于预 置阈值；
[0049] S222:在所述判断结果为是时，在所述预定时空窗内获取与所述采样点相同采样 时刻的采样点的振幅值；
[0050] S223 :对所述相同采样时刻的振幅值的绝对值进行排序，取所述排序后中间的K 个振幅值的平均值作为所述采样点在所述预定时空窗内的参考振幅，K多1。
[0051] 所述的预定阈值可以根据数据处理速度或其他需求进行设置。本实施例中可以设 置为3,可以表示如果所述设置的预定时空窗中所包括的地震道数大于3,则可以采取上述 所述的方法计算所述参考振幅。在本实施例中，可以在所述预定时空窗内分别获取与所述 当前采样点相同采样时刻的其他地震道中采样点的振幅值，例如当前采样点的采样时刻T 为200毫秒时，可以获取预定时空窗内的7道地震道中所有在采样时刻为200毫秒的采样 点的振幅值。然后可以对例如所述获取的同一采样时刻的7个采样点的振幅值的绝对值进 行从小到大或者从大到小进行排序，选取排序后的中间的K个振幅值。进一步的以所述K 个振幅值的平均值作为所述采样点在所述预定时空窗内的参考振幅。当然，所述K的取值 也可以根据数据处理需求进行设置，例如可以设置为3或者5。
[0052] 优选的实施例中，所述的预置阈值可以设置为3。所述K的取值范围可以设置为： 3 < K < 5。本实施例中所述的参考振幅的计算方法，在所述预定时空窗内包括的地震道数 较多时，可以保障计算得到合适的参考振幅的同时还可以加快数据处理速度。
[0053] 本发明所述的计算参考振幅的方法并不限于取采样点平均值的方式计算得出。其 他的一些实施例中，也可以采用以所述预定时空窗内采样点均方根振幅值作为所述参考振 幅。具体的本发明所述方法中所述计算所述待处理道集数据的采样点在预定时空窗内的参 考振幅可以包括：
[0054] S23 :计算所述待处理道集数据采样点在预定时空窗内所有采样点均方根振幅，以 所述均方根振幅作为所述参考点在所述预定时空窗内的参考振幅。
[0055] 例如当前所述预定时空窗内共有J个采样点，可以采用下式计算所述预定时空窗 内J个采样点的均方根振幅A，以所述均方根振幅A作为当前采样点P的参考振幅：
[0057] 采样采样点的均方根振幅可以使计算得出的参考振幅更加准确、有效。当然，参照 本发明的其他实施例，例如所述预定时空窗内的包括的地震道数大于预定阈值时，也可以 参照实施例S221-至S223的方式计算所述预定时空窗内采样点的均方根振幅，然后进一步 计算得到参考振幅。基于常规的处理方法及本领域无需创造性劳动的其他数据处理的变形 方式都可以实现本发明参考振幅的计算，也都应属于本发明的保护范围，在此不做逐一举 例说明。
[0058]为采样点设置预定时空窗，计算所述待处理道集数据的采样点在所述预定时空窗 内的参考振幅。
[0059] S3:比较所述采样点的振幅与所述参考振幅，判断所述采样点的振幅是否大于基 于所述参考振幅设置的门槛值；将所述判断结果为是的采样点的振幅作为异常振幅；移动 所述预定时空窗至下一个采样点，判断所述下一个采样点是否为异常振幅。
[0060]计算得到当前采样点在所述预定时空窗内的参考振幅后，可以将所述当前采样点 的振幅与所述参考振幅进行比较。本发明中可以设置基于所述参考振幅设置一个判断所述 当前采样点的振幅是否为异常振幅的门槛值。这个门槛值取值范围可以为所述参考振幅的 倍数，例如可以根据数据处理的迭代次数初始设置为15,或者在迭代数次后设置为4或者2 等，具体的可以数据处理需求进行设置。本实施例中所述基于所述参考振幅设置的门槛值 具体设置的取值范围可以包括：
[0061]所述门槛值为所述参考振幅的Q倍，所述Q的取值范围为：0< Q < 20。
[0062] 判断所述当前采样点的振幅值是否大于所述基于参考振幅设置的门槛值，例如当 前采样点的振幅大于所述参考振幅的10倍时认为所述当前采样点的振幅为异常振幅。一 般的，在地震数据异常振幅处理过程中可以结合迭代算法多次计算采样点的参考振幅，迭 代的次数越多，计算得到的参考振幅越准确、有效。因此，所述的门槛值在实际的数据处理 中可以为随着不同的迭代次数不断变化的数值，通常可以随着迭代次数的增加而减小。
[0063] 可以根据所述门槛值判断当前采样点的振幅是否为异常振幅。在判断完当前采样 点后，可以移动所述预定时空窗，计算下一个采样点的参考振幅，判断所述下一个采样点的 振幅是否为异常振幅。
[0064] S4 :对所述判断为异常振幅的采样点振幅进行振幅衰减。
[0065] 逐一计算判断所述待处理道集数据中每个地震道中采样点的振幅是否为异常振 幅，对判断为异常振幅的采样点的振幅，可以根据预先设置的衰减方法进行振幅衰减，以压 制异常振幅干扰。具体的所述振幅衰减的方法可以包括：
[0066] S41 :将所述异常振幅与预先设置的衰减系数运算后得到的振幅值作为衰减后的 振幅值。
[0067] 可以预先设置对异常振幅进行衰减的衰减系数Coef，所述衰减系数可以根据当 前处理的数据采取一定方法