技术领域
[0001] 本
发明涉及油气地球物理勘探领域，更具体地，涉及一种地震采集激发点空间位置检测方法及系统。
背景技术
[0002] 当前
数据采集都是使用高
密度、小道距、高
覆盖的海量数据，生产中对施工要求越来越高，在实际施工区域地表障碍物增多，导致发生偏差的激发点数量增加，发生偏差的方向性也复杂起来，这会影响剖面
叠加的
精度，常规的线性动校正已不能满足生产需求。因此，有必要开发一种地震采集激发点空间位置检测方法及系统。
[0003] 公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的
现有技术。
发明内容
[0004] 本发明提出了一种地震采集激发点空间位置检测方法及系统，其能够通过获取偏离点误差，确定校正量进行位置校正，可以快速检测判断激发点的实际位置。
[0005] 根据本发明的一方面，提出了一种地震采集激发点空间位置检测方法。所述方法可以包括：在
地震道上划分多个时窗，计算相邻时窗的
能量比；获得能量比最大的后一个时窗，进而确定初至时间；根据所述初至时间进行反演，计算近地表速度，进而计算检波点误差；设定误差
阈值，若所述检波点误差大于所述误差阈值，则对应的炮点为偏移点；重新计算所述检波点误差，选取检波点误差最小的点作为实际炮点位置。
[0006] 优选地，计算相邻时窗的能量比为：
[0007]
[0008] 其中，A为时窗内的初至前后能量比，x(t)为
信号振幅值，T0为前一个时窗的起点，T1为前一个时窗的终点与后一个时窗的起点，T2为后一个时窗的终点。
[0009] 优选地，根据所述初至时间进行反演，计算近地表速度，进而计算检波点误差包括：根据所述初至时间进行反演，计算所述近地表速度；根据所述近地表速度，计算检波点平均迭加
波速；根据所述检波点平均迭加波速，计算检波点滑动平均速度；根据所述检波点滑动平均速度，计算所述检波点误差。
[0010] 优选地，所述近地表速度为：
[0011]
[0012] 其中，vi为近地表速度，Di为偏移距，ti为初至时间，i表示地震道序号，i＝1,2,…,n。
[0013] 优选地，所述检波点平均迭加波速为：
[0014]
[0015] 其中，vs[j]为第j号检波点的平均迭加波速，vj[n]表示第j号检波点的迭加波速，M为第j号检波点处迭加的常规波速的次数。
[0016] 优选地，所述检波点滑动平均速度为：
[0017]
[0018] 其中，va[j]表示第j号检波点的滑动平均速度。
[0019] 优选地，所述检波点误差为：
[0020]
[0021] 其中，e为检波点误差。
[0022] 优选地，还包括：根据最小检波点误差，计算校正量；根据所述校正量校正其他偏移点。
[0023] 根据本发明的另一方面，提出了一种地震采集激发点空间位置检测系统，其特征在于，该系统包括：
存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：在地震道上划分多个时窗，计算相邻时窗的能量比；获得能量比最大的后一个时窗，进而确定初至时间；根据所述初至时间进行反演，计算近地表速度，进而计算检波点误差；设定误差阈值，若所述检波点误差大于所述误差阈值，则对应的炮点为偏移点；重新计算所述检波点误差，选取检波点误差最小的点作为实际炮点位置。
[0024] 优选地，计算相邻时窗的能量比为：
[0025]
[0026] 其中，A为时窗内的初至前后能量比，x(t)为信号振幅值，T0为前一个时窗的起点，T1为前一个时窗的终点与后一个时窗的起点，T2为后一个时窗的终点。
[0027] 优选地，根据所述初至时间进行反演，计算近地表速度，进而计算检波点误差包括：根据所述初至时间进行反演，计算所述近地表速度；根据所述近地表速度，计算检波点平均迭加波速；根据所述检波点平均迭加波速，计算检波点滑动平均速度；根据所述检波点滑动平均速度，计算所述检波点误差。
[0028] 优选地，所述近地表速度为：
[0029]
[0030] 其中，vi为近地表速度，Di为偏移距，ti为初至时间，i表示地震道序号，i＝1,2,…,n。
[0031] 优选地，所述检波点平均迭加波速为：
[0032]
[0033] 其中，vs[j]为第j号检波点的平均迭加波速，vj[n]表示第j号检波点的迭加波速，M为第j号检波点处迭加的常规波速的次数。
[0034] 优选地，所述检波点滑动平均速度为：
[0035]
[0036] 其中，va[j]表示第j号检波点的滑动平均速度。
[0037] 优选地，所述检波点误差为：
[0038]
[0039] 其中，e为检波点误差。
[0040] 优选地，还包括：根据最小检波点误差，计算校正量；根据所述校正量校正其他偏移点。
[0041] 本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的
附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
[0042] 通过结合附图对本发明示例性
实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
[0043] 图1示出了根据本发明的地震采集激发点空间位置检测方法的步骤的
流程图。
具体实施方式
[0044] 下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。
[0045] 图1示出了根据本发明的地震采集激发点空间位置检测方法的步骤的流程图。
[0046] 在该实施例中，根据本发明的地震采集激发点空间位置检测方法可以包括：在地震道上划分多个时窗，计算相邻时窗的能量比；获得能量比最大的后一个时窗，进而确定初至时间；根据初至时间进行反演，计算近地表速度，进而计算检波点误差；设定误差阈值，若检波点误差大于误差阈值，则对应的炮点为偏移点；重新计算检波点误差，选取检波点误差最小的点作为实际炮点位置。
[0047] 在一个示例中，计算相邻时窗的能量比为：
[0048]
[0049] 其中，A为时窗内的初至前后能量比，x(t)为信号振幅值，T0为前一个时窗的起点，T1为前一个时窗的终点与后一个时窗的起点，T2为后一个时窗的终点。
[0050] 在一个示例中，根据初至时间进行反演，计算近地表速度，进而计算检波点误差包括：根据初至时间进行反演，计算近地表速度；根据近地表速度，计算检波点平均迭加波速；根据检波点平均迭加波速，计算检波点滑动平均速度；根据检波点滑动平均速度，计算检波点误差。
[0051] 在一个示例中，近地表速度为：
[0052]
[0053] 其中，vi为近地表速度，Di为偏移距，ti为初至时间，i表示地震道序号，i＝1,2,…,n。
[0054] 在一个示例中，检波点平均迭加波速为：
[0055]
[0056] 其中，vs[j]为第j号检波点的平均迭加波速，vj[n]表示第j号检波点的迭加波速，M为第j号检波点处迭加的常规波速的次数。
[0057] 在一个示例中，检波点滑动平均速度为：
[0058]
[0059] 其中，va[j]表示第j号检波点的滑动平均速度。
[0060] 在一个示例中，检波点误差为：
[0061]
[0062] 其中，e为检波点误差。
[0063] 在一个示例中，还包括：根据最小检波点误差，计算校正量；根据校正量校正其他偏移点。
[0064] 具体地，根据本发明的地震采集激发点空间位置检测方法可以包括：
[0065] 地震记录前后的信号有非常大的差异，利用选取时窗进行初至时间拾取，初至前后时窗能量比的极大值是检查初至时间的方法。在地震道上划分多个时窗，通过公式(1)计算相邻时窗的能量比；获得最大能量比对应的两个时窗的后一个时窗，在该时窗内确定能量最大的时间，即为初至时间。
[0066] 根据初至时间进行反演，对于面元范围内不同偏移距的地震道，通过公式(2)计算近地表速度；炮检点位置的偏移会给反演精度产生误差，目前可以通过对各个检波点处迭
加速度进行平均，减少炮点位置偏差而产成的误差，根据近地表速度，通过公式(3)计算检波点平均迭加波速；用滑动平均波速修正炮点和各检波点间的近地表速度，减少近地表因素对波速造成的影响，根据某个检波点和周围4个检波点的平均迭加波速，通过公式(4)计算检波点滑动平均速度；用地表平均速度计算直达波时间，再计算直达波时间与初至时间之间的误差，即可以得到该检波点处的时间误差，即根据检波点滑动平均速度，通过公式(5)计算检波点误差，误差值越大，检波点偏移的位置就越远。
[0067] 设定误差阈值，若检波点误差小于等于误差阈值，则对应的炮点为正常点，若检波点误差大于误差阈值，则对应的炮点为偏移点；重新计算检波点误差，选取检波点误差最小的点作为实际炮点位置；根据最小检波点误差，计算校正量；根据校正量校正其他偏移点。
[0068] 本方法通过获取偏离点误差，确定校正量进行位置校正，可以快速检测判断激发点的实际位置。
[0069] 应用示例
[0070] 为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
[0071] 根据本发明的地震采集激发点空间位置检测方法包括：
[0072] 在地震道上划分多个时窗，通过公式(1)计算相邻时窗的能量比；获得最大能量比对应的两个时窗的后一个时窗，在该时窗内确定能量最大的时间，即为初至时间。
[0073] 根据初至时间进行反演，对于面元范围内不同偏移距的地震道，通过公式(2)计算近地表速度；炮检点位置的偏移会给反演精度产生误差，目前可以通过对各个检波点处迭加速度进行平均，减少炮点位置偏差而产成的误差，根据近地表速度，通过公式(3)计算检波点平均迭加波速；用滑动平均波速修正炮点和各检波点间的近地表速度，减少近地表因素对波速造成的影响，根据某个检波点和周围4个检波点的平均迭加波速，通过公式(4)计算检波点滑动平均速度；用地表平均速度计算直达波时间，再计算直达波时间与初至时间之间的误差，即可以得到该检波点处的时间误差，即根据检波点滑动平均速度，通过公式(5)计算检波点误差，误差值越大，检波点偏移的位置就越远。
[0074] 设定误差阈值，若检波点误差大于误差阈值，则对应的炮点为偏移点；重新计算检波点误差，选取检波点误差最小的点作为实际炮点位置；根据最小检波点误差，计算校正量；根据校正量校正其他偏移点。
[0075] 野外采集现场的激发点空间位置检测是实时的，要求校正量计算必须迅速准确。以某工区某日施工记录为例，当天采集放炮752炮，用此方法检查出位置偏移点一共8炮，偏移率1.06％，并且迅速计算出偏移的方向和距离。这样，可以对其进行实时校正，不用
修改观测系统，减少了重新加载观测系统的麻烦，确保了野外采集激发点偏离点位置校正工作的精准。本方法已经在野外采集现场多次使用，均能够自动计算偏移点距离，判断偏移方向，减少了人为因素的影响，激发点检查结果更加准确可靠。
[0076] 本发明通过获取偏离点误差，确定校正量进行位置校正，可以快速检测判断激发点的实际位置。
[0077] 本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
[0078] 根据本发明的实施例，提供了一种地震采集激发点空间位置检测系统，其特征在于，该系统包括：存储器，存储有计算机可执行指令；处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：在地震道上划分多个时窗，计算相邻时窗的能量比；获得能量比最大的后一个时窗，进而确定初至时间；根据初至时间进行反演，计算近地表速度，进而计算检波点误差；设定误差阈值，若检波点误差大于误差阈值，则对应的炮点为偏移点；重新计算检波点误差，选取检波点误差最小的点作为实际炮点位置。
[0079] 在一个示例中，计算相邻时窗的能量比为：
[0080]
[0081] 其中，A为时窗内的初至前后能量比，x(t)为信号振幅值，T0为前一个时窗的起点，T1为前一个时窗的终点与后一个时窗的起点，T2为后一个时窗的终点。
[0082] 在一个示例中，根据初至时间进行反演，计算近地表速度，进而计算检波点误差包括：根据初至时间进行反演，计算近地表速度；根据近地表速度，计算检波点平均迭加波速；根据检波点平均迭加波速，计算检波点滑动平均速度；根据检波点滑动平均速度，计算检波点误差。
[0083] 在一个示例中，近地表速度为：
[0084]
[0085] 其中，vi为近地表速度，Di为偏移距，ti为初至时间，i表示地震道序号，i＝1,2,…,n。
[0086] 在一个示例中，检波点平均迭加波速为：
[0087]
[0088] 其中，vs[j]为第j号检波点的平均迭加波速，vj[n]表示第j号检波点的迭加波速，M为第j号检波点处迭加的常规波速的次数。
[0089] 在一个示例中，检波点滑动平均速度为：
[0090]
[0091] 其中，va[j]表示第j号检波点的滑动平均速度。
[0092] 在一个示例中，检波点误差为：
[0093]
[0094] 其中，e为检波点误差。
[0095] 在一个示例中，还包括：根据最小检波点误差，计算校正量；根据校正量校正其他偏移点。
[0096] 本系统通过获取偏离点误差，确定校正量进行位置校正，可以快速检测判断激发点的实际位置。
[0097] 本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
[0098] 以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
