本发明涉及矿产、油气、页岩油气及煤层气地震勘探与开发领域，更具体地，涉及一种地震数据废道剔除方法及装置。
背景技术
在地震数据处理中，我们首先会对每个炮集的数据进行废道剔除工作，以便后期的地震数据处理更加准确可靠。
news地震综合解释系统(以下简称news系统)是中国石化自主研发的地震综合解释软件。news系统中的叠前反演会用到共炮集数据或其他一些叠前数据，这些数据往往也会存在一部分废道，因此我们需要研发一个可靠及有效的废道剔除方法以解决上述问题。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现要素：
本发明的目的在于解决上述现有技术中存在的难题，提出了一种基于极值点类比的地震数据废道剔除方法及装置。
根据本发明的一方面，提出一种地震道数据废道剔除方法。该方法可以包括以下步骤：
1)单道数据波峰波谷点识别；
2)统计地震数据中每一道地震道的波峰波谷数量值k以及地震道波峰波谷数量平均值ave；
3)将每一道的波峰波谷数量值k与设定阈值f进行比较，如果k>＝f则将该地震道确定为废道；
4)剔除所确定的废道。
优选地，波峰识别公式为：
其中，xn为单个地震道序列中的任意一个样点数，xn、xn+1、xn+2为连续三个样点数。
优选地，波谷识别公式为：
其中，xn为单个地震道序列中的任意一个样点数，xn、xn+1、xn+2为连续三个样点数。
优选地，步骤2)中地震道波峰波谷数量平均值ave根据以下公式计算：
其中，k为每一道的波峰波谷数量值，n为地震数据体中的地震道的道数值。
优选地，设定阈值f在地震道波峰波谷数量平均值ave的15％～25％之间。
根据本发明的另一方面，提出一种地震道废道剔除装置。该装置包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现以下步骤：
1)单道数据波峰波谷点识别；
2)统计地震数据中每一道地震道的波峰波谷数量值k以及地震道波峰波谷数量平均值ave；
3)将每一道的波峰波谷数量值k与设定阈值f进行比较，如果k>＝f则将该地震道确定为废道；
4)剔除所确定的废道。
优选地，波峰识别公式为：
其中，xn为单个地震道序列中的任意一个样点数，xn、xn+1、xn+2为连续三个样点数。
优选地，波谷识别公式为：
其中，xn为单个地震道序列中的任意一个样点数，xn、xn+1、xn+2为连续三个样点数。
优选地，步骤2)中地震道波峰波谷数量平均值ave根据以下公式计算：
其中，k为每一道的波峰波谷数量值，n为地震数据体中的地震道的道数值。
优选地，设定阈值f在地震道波峰波谷数量平均值ave的15％～25％之间。
根据本发明的地震道数据废道剔除方法和装置可以有效地定位废道，为矿产、常规油气、页岩油气及煤层气的勘探、开发服务。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1为根据本发明的示例性实施方案的地震道数据废道剔除方法的流程图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
以下参考图1详细描述根据本发明示例性实施例的地震道数据废道剔除方法。该方法主要包括以下步骤：
步骤1：单道数据波峰波谷点识别。
单个地震道为一个浮点数序列，此序列中的任意一点可描述为xn。波峰的特点是连续的三个样点数的中间数为三数中最大值，波谷的特点是连续的三个样点数的中间数为三数中最小值。三个数的值相同，不做判断。
波峰识别公式为：
其中，xn为单个地震道序列中的任意一个样点数，xn、xn+1、xn+2为连续三个样点数。
波谷识别公式为：
其中，xn为单个地震道序列中的任意一个样点数，xn、xn+1、xn+2为连续三个样点数。
利用公式(1)识别波峰，利用公式(2)识别波谷，并记录每一道的波峰波谷数量k和位置。在地震数据中用线号、炮号、道号确定数据的具体位置。
步骤2：统计地震数据中每一道地震道的波峰波谷数量值k以及地震道波峰波谷数量平均值ave。
根据以下公式计算地震道波峰波谷数量平均值ave：
其中，k为每一道的波峰波谷数量值，n为地震数据体中的地震道的道数值。
步骤3：将每一道的波峰波谷数量值k与设定阈值f进行比较，如果k>＝f则将该地震道确定为废道。
地震数据中废道的波峰波谷点要明显多于正常道，因此可以通过设定一个阈值f来对废道进行定位。
这个阈值f可以在数据体中地震道波峰波谷数量平均值ave的15％～25％之间，然后用每一道的的波峰波谷数量值k去和f做比较，如果k大于等于f，则这一道为废道；如果k小于f，则这一道为正常道。
例如，将f设定为地震道波峰波谷数量平均值ave的20％，即f＝ave*20％。
如果k>＝f,则这一道为废道；如果k<f,则这一道为正常道。
由于在第一步中我们已记录了每道的位置，因此根据计算结果可以顺利定位废道。
步骤4：剔除所确定的废道。
根据该示例性实施例的地震道废道剔除方法及装置，可以有效地定位废道，为矿产、常规油气、页岩油气及煤层气的勘探、开发服务。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。