校正后初始时间；还用于根据所述校正后起始时间和设置的时窗长度确定所述近炮点地 震道的处理时窗。
[0071] 优选的实施例中，所述根据所述炮集数据记录的排列中以所述目标炮的炮点对称 的地震道的炮检距差和高程差计算所述排列的折射速度包括；
[0072] 利用下式计算得到风化层慢度P0和折射层慢度P1 :
[0073]
[0074] 上式中，P。表示为风化层慢度，Pi表示为折射层慢度，Wlik表示为第i个排列中第 k对的最大相似系数，Tlik表示为所述Wlik对应的时差为，Slik表示为所述比对地震道的炮检 距差，Hlik表示为所述比对地震道的高程差，k表示为第i个排列中所述比对地震道的对数；
[0075] 根据计算得到风化层慢度P0和折射层慢度P1计算所述目标炮的中排列的折射速 度。
[0076] 优选的实施例中，所述炮偏估值计算模块包括：
[0077] 炮偏门槛值模块，用于存储预设设置的炮偏门槛值，所述炮偏门槛值的取值范围 为0. 5至1. 5个地震道的道间距；
[0078] 比较模块，用于将所述炮偏估值的均值与所述炮偏门槛值进行比较，根据比较的 结果判断所述目标炮是否存在炮偏。
[0079] 本发明提供的一种炮偏检测方法及装置，利用浅层折射波，对炮点两侧的近炮道 进行相关计算，拾取最大相似系数所对应的时差，结合时差及相似系数，综合进行炮偏检 测、诊断。实施过程中，可以在炮集数据记录中选取一炮的近炮点地震道，这样根据对炮偏 较为敏感的近炮点地震道数据来检测炮偏，可以提高炮偏的检测精度。然后在可以确定的 时窗范围内根据以炮点为对称的比对地震道的最大相似系数计算排列的平均时差，和根据 所述比对地震的炮检距差和高程差计算排列的折射速度，进一步采用预先设定的算法计算 得到所述排列的炮偏估值。本发明中可以以目标炮的炮偏估值的均值作为判断所述目标是 否存在炮偏的判断标准，这样，可以排除人工手动检查不稳定性，提高较小偏离的识别率， 提高检查精度。同时还可以自动快速的进行炮偏检测，降低人工工作强度，提高炮偏数据的 处理效率。
【附图说明】
[0080] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提 下，还可以根据这些附图获得其他的附图。
[0081] 图1是本发明一种炮偏检测方法一种实施例的方法流程示意图；
[0082] 图2是本发明选取近炮点地震道的选取示意图；
[0083] 图3是本发明所述一种炮偏检测装置一种实施例的模块结构示意图；
[0084] 图4是本发明炮偏检测装置中所述近炮道选取模块一种实施例的模块结构示意 图；
[0085] 图5是本发明炮偏检测装置中所述近炮道选取模块另一种实施例的模块结构示 意图；
[0086] 图6是本发明炮偏检测装置中所述平均时差计算模块一种实施例的模块结构示 思；
[0087] 图7是本发明炮偏检测装置中所述时窗设置模块一种实施例的模块结构示意图；
[0088] 图8是本发明炮偏检测装置中所述炮偏估值计算模块一种实施例的模块结构示 意图；
[0089] 图9是利用本发明对理论测试数据进行炮偏检测的炮偏估算值与实际值的对比 示意图；
[0090] 图10是利用本发明对某一区域实际地震资料进行炮偏检测的炮偏估算值与实际 值大小的对比示意图。
【具体实施方式】
[0091] 为了使本技术领域的人员更好地理解本发明中的技术方案，下面将结合本发明实 施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施 例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通 技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护 的范围。
[0092] 在石油勘探采集地震资料时，通常需要将接收地震数据的检波点（也可以称为接 收点或者地震道）按照一定的测线分布进行排列。例如可以排布成以东西水平方向为接收 线方向的纵偏，和垂直所述接收线方向的横偏，当然，所述接收线方向也可以为南北方向。 通常的，以测线接收线方向上的地震道为排列进行的炮偏检测可以称为纵偏检测，同样，若 以与所述测线接收线方向垂直方向上的地震道为排列进行的炮偏检测可以称为横偏检测。 本发明所述的炮偏检测方法或装置在利用地震数据检测炮偏时可以采用上述所述的纵偏 检测方式，在一些应用场景中也可以适用于横偏检测的方式。上述两种检测方式区别通常 在于炮偏检测时处理数据的选择方式不同。
[0093] 下面结合附图以纵偏检测为例对本发明所述的炮偏检测方法及装置进行详细的 说明。图1是本发明所述炮偏检测方法一种实施例的方法流程图。虽然本发明提供了如下 述实施例或附图所示的方法操作步骤，但基于常规或者无需创造性的劳动在所述方法中可 以包括更多或者更少的操作步骤。在逻辑性上不存在必要因果关系的步骤中，这些步骤的 执行顺序不限于本发明实施例提供的执行顺序。所述的方法的在实际中的装置、系统或终 端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行（例如并行处 理器或者多线程处理的环境）。
[0094] 具体的本发明提供的一种炮偏检测方法的一种实施例如图1所示，所述方法可以 包括：
[0095] SI:选取炮集数据记录中目标炮的近炮点地震道。
[0096] -般的，在实际地震勘探数据采集时，可以设置有一个或者多个炮点。然后可以在 不同的炮点进行一次或者多次震源激发，由分布在接收线上的多个检波点位置的检波器接 收地震波信号，形成地震道数据记录。通常情况下，将一炮的所述检波器接收的该炮的地震 数据记录按照一定的规则排列后（如CMP、CDP等）形成的数据记录可以称为单炮集数据 记录。同样，将多点或者多次放炮形成的数据记录可以称为多炮集数据记录。本发明中所 述的炮集数据记录可以为多炮集数据记录。在地震数据采集时，根据炮点震源与检波点的 相对位置分布通常可以分为双边放炮的情况。本发明所述的炮偏检测方法可以处理双边放 炮的地震数据。以纵偏检查为例，本发明所述的炮偏检测方法可以要求沿接收线方向为双 边放炮，垂直于接收线方向则不作限制。横偏检查则相反，要求垂直于接收线方向为双边放 炮，沿接收线方向不作限制
[0097] 本发明提供的炮偏检测方法可以按照预先设定的选取规则从所述炮集数据记录 中选取需要进行炮偏检测的目标炮的近炮点地震道。一般的，在震源激发时，越靠近炮点的 地震道受到的影响越大。在本发明方案中，所述近炮点地震道中的地震数据通常对所述目 标炮的炮偏较为敏感，因此可以利用靠近所述目标炮的近炮点地震道进行对目标炮的炮偏 检测，可以有效提高所述炮偏检测结果的精度。所述近炮点地震道的选取通常是以需要进 行炮偏检测的目标炮为中心，选取周围靠近所述目标炮的地震道。本发明可以提供一种选 取所述近炮点地震道的实施方式，具体的，所述选取炮集数据记录中目标炮的近炮点地震 道可以包括：
[0098] 以所述炮集数据记录中目标炮的炮点为对称中心选取2n个地震道排列，以所述 炮点为对称中心在每一个所述地震道排列中选取2m个地震道形成所述近炮点地震道，n为 设置的在所述炮点单侧选取地震道排列的第一预设值，m为设置的在所述地震道排列中炮 点的单侧选取地震道的第二预设值，n、m为自然数。
[0099] 当然，本发明实施例在选取所述近炮点地震的排列数n和每个排列上的地震道数 m可以根据数据处理需求进行设置。
[0100] 具体的本实施例以纵偏检测为例对选取近炮点地震道的过程进行说明。图2是 选取近炮点地震道的选取示意图。在选取近炮点地震道时，对于所述炮集数据记录中的每 一炮，可以选择炮检距在Y方向（垂直接收线方向）离所述目标炮点较近的若干排列。对 应Y方向双边放炮的情况，例如图2所示，当目标炮的炮点01在两条接收线之间时（图2 中*所在位置），可以以所述目标炮的炮点01为对称中心选择所述近炮点地震道的排列数 为2n，即所述炮点两侧各n个排列。对于所述炮点02落在某条接收线上时（图2中▲所 在位置），则可以以所述炮点02所在地震道排列为对称中心选择所述近炮点地震的排列数 为2n+l，其中可以包括炮点所在的排列。进一步的，对应选取的每个排列，可以选择X方向 (沿接收线方向）距离所述目标炮的炮点最近的2m个地震道，即所述炮点两侧各m个地震 道。若在分布上，所述炮点与接收点位置相同，即震源激发点位置的地震道，例如图2中鲁 所在位置，此时可以选择剔除离所述目标炮的炮点最近的位置的数据值，再按照上述在炮 点两侧选择2m个地震道。
[0101] 当然，本发明在选取所述近炮点地震道时可以根据预先分析、研究设置的特定计 算方法进行选取。本发明的提供了选取所述近炮点地震道的另一种实施方式，具体的选取 所述近炮点地震道时，可以假设所述目标炮的炮点、接收点坐标处于相对坐标系中。所述的 相对坐标系一般的可以指正交坐标系的横轴与接收线方向平行，纵轴垂直于接收线方向， 在所述相对坐标系下，每个排列中接收点的Y坐标基本相同，便于区分不同排列。在所述的 相对坐标系中，假设接收线之间的线距为deltgy，接收线中接收点之间的点距为deltgx， 所述目标炮的第i个接收点在Y方向上的炮检距为Ayi、在X方向上的炮检距为AXi，则可 以根据下述公式选取第i个接收点的地震道作为近炮点地震道：
[0102]
[0103] 上述公式中，L可以表示为所述目标炮的接收点总数，相应的，所述deltgy可以 为所述炮集数据记录接收线之间的线距，deltgx可以为所述接收线中接收点之间的点距， A71可以为所述目标炮的第i个接收点在建立的相对坐标系中Y方向上的炮检距，Ax河 以为所述目标炮的第i个接收点在建立的相对坐标系中X方向上的炮检距。可以以所述满 足上述公式的接收点对应的地震道作为所述近炮点地震道。
[0104] 本实施例根据公式选出所述近炮点地震道后还可以根据所述近炮点地震道对所 述炮集数据记录进行分析，判断是否存在异常单炮数据，以便及时输出供人工检查处理。具 体的，例如纵偏检查时，根据AXi的符号，分别计算出正、负排列的道数ML、MR，当ML与MR 不等时，该排列可能缺道，或者说明该炮在X方向为单边放炮，此时，可以输出该炮数据供 人工检查。若横偏检查，则可以根据Ayi的符号，分别计算出炮点两侧的排列个数NL、NR， 若NL等于零或者NR等于零，则说明Y方向单边放炮时，此时，可以输出该炮供人工检查。
[0105] 本发明进行炮偏检测时可以选取炮集数据记录中需要进行炮偏检测的目标炮的 近炮点地震道。
[0106] S2 :确定所述近炮点地震道的处理时窗。
[0107] 本发明的一种实施例获取所述近炮点地震道后，还可以进一步设置地震数据的处 理时窗范围，例如可以取值在300毫秒至500毫秒的范围。本实施例中所述处理时窗的大 小可以根据实际地震资料、数据处理需求或经验值等进行设定
[0108] 在一些实施例中，为减小浅层噪声对炮偏检测的影响，在实施例中所述的处理时 窗的起始时间可以优选在初至位置附近。具体的可以利用浅层折射波进行互相关计