技术领域
[0001] 本
发明涉及
地震波初至检测技术领域，特别是涉及一种地震波初至检测方法及系统。
背景技术
[0002] 地震波初至在矿物资源开发、地震
定位以及地震
层析成像中都起着十分关键的作用。在石油勘探技术中，地震波
信号可以用来了解地质构造，例如石油勘探就需要利用人工地震信号，其中准确检测到地震波到达检波点的初至时刻尤为重要。
[0003] 从根本上说，地震信号的初至波既包含地震信号又包含噪声信号，是两者的
叠加，初至拾取的精确性是地震资料处理中
基础而又重要的环节。传统的初至拾取主要依靠手工，但人工拾取耗时且精确度低。为了提高地震波初至的提取效率与提取速度，现在已经出现了自动拾取的方法，例如
能量比法、最大振幅法、分形维法等。
[0004] 初至自动拾取方法中最常用的是基于能量比的方法。根据能量比达到最大值时的时间点就可以估算出地震波初至到来的时刻。但值得注意的是，当
信噪比相对较低时这种基于能量比的方法就十分不理想。与能量比法原理相类似，振幅特性也常被用来检测初至，利用这种方法，根据地震波初至前后的振幅比的变化来拾取初至。一般来说，当信噪比足够高时这些方法是有效的。但我国石油勘探环境并不都是理想的，地表复杂的地区噪声干扰非常严重，这导致了现有方法对低信噪比情况下的初至拾取并不理想，
精度并不高。
发明内容
[0005] 本发明的目的是提供一种地震波初至检测方法及系统，有效避免噪声对检测精度的影响，提高检测准确度。
[0006] 为实现上述目的，本发明提供了如下方案：
[0007] 一种地震波初至检测方法，包括：
[0008] 获取
检波器采集的待检测的地震记录序列；
[0009] 建立第一
时间窗和第二时间窗，所述第一时间窗和所述第二时间窗的容量均为L个数据，其中L为奇数；所述第二时间窗紧随所述第一时间窗；
[0010] 对所述地震记录序列的左侧和右侧进行数据填充，使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动至所述第一时间窗的中心对准所述地震记录序列的最后一个数据时，所述第一时间窗和所述第二时间窗内均填充满L个数据；
[0011] 从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动所述第一时间窗和所述第二时间窗，直至所述第一时间窗的中心对准所述地震记录序列的最后一个数据，每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述地震记录序列中的各个数据处是否有地震
波动到达；
[0012] 记录所有有地震波动到达的数据在所述地震记录序列中的序号，并计算每个有地震波动到达的数据对应的初至。
[0013] 可选的，所述对所述地震记录序列的左侧和右侧进行数据填充，具体包括：
[0014] 将所述地震记录序列的左侧补充至少 个数据；左侧补充的至少 个数据均与所述地震记录序列最左侧的数据相同；
[0015] 将所述地震记录序列的右侧补充至少 个数据；右侧补充的至少 个数据均与所述地震记录序列最右侧的数据相同。
[0016] 可选的，所述每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述地震记录序列中的各个数据处是否有地震波动到达，具体包括：
[0017] 求取所述第一时间窗内所有数据的最大值和最小值，得到取值区间；
[0018] 将所述取值区间划分为K段，得到K段小区间；
[0019] 计算所述第一时间窗内的数据介于每个所述小区间内的数量F(n)P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F(n)Q(k)；其中n表示位于所述第一时间窗中心的数据在所述地震记录序列中的序号；k为小区间序号，k＝1,2,…,K；
[0020] 计算数量差值比率：
[0021]
[0022] 根据参数K和
置信度α，通过查询卡方表获取卡方值χ2(α,K-1)；
[0023] 判断所述数量差值比率是否大于或等于所述卡方值，得到第一判断结果；
[0024] 若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值，则确定在所述地震记录序列的第n个数据处有地震波动到达；若所述第一判断结果表示所述数量差值比率小于所述卡方值，则确定在所述地震记录序列的第n个数据处没有地震波动到达。
[0025] 本发明还公开一种地震波初至检测系统，包括：
[0026] 序列获取模
块，用于获取检波器采集的地震记录序列；
[0027] 时间窗建立模块，用于建立第一时间窗和第二时间窗，所述第一时间窗和所述第二时间窗的容量均为L个数据，其中L为奇数；所述第二时间窗紧随所述第一时间窗；
[0028] 数据填充模块，用于对所述地震记录序列的左侧和右侧进行数据填充，使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动至所述第一时间窗的中心对准所述地震记录序列的最后一个数据时，所述第一时间窗和所述第二时间窗内均填充满L个数据；
[0029] 地震波初至检测模块，用于从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动所述第一时间窗和所述第二时间窗，直至所述第一时间窗的中心对准所述地震记录序列的最后一个数据，每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述地震记录序列中的各个数据处是否有地震波动到达；
[0030] 记录与初至计算模块，用于记录所有有地震波动到达的数据在所述地震记录序列中的序号，并计算每个有地震波动到达的数据对应的初至。
[0031] 可选的，所述数据填充模块，具体包括：
[0032] 序列左填充单元，用于将所述地震记录序列的左侧补充至少 个数据；左侧补充的至少 个数据均与所述地震记录序列最左侧的数据相同；
[0033] 序列右填充单元，用于将所述地震记录序列的右侧补充至少 个数据；右侧补充的至少 个数据均与所述地震记录序列最右侧的数据相同。
[0034] 可选的，所述地震波初至检测模块，具体包括：
[0035] 区间计算单元，用于求取所述第一时间窗内所有数据的最大值和最小值，得到取值区间；
[0036] 区间划分单元，用于将所述取值区间划分为K段，得到K段小区间；
[0037] 区间数据量计算单元，用于计算所述第一时间窗内的数据介于每个所述小区间内的数量F(n)P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F(n)Q(k)；其中n表示位于所述第一时间窗中心的数据在所述地震记录序列中的序号；k为小区间序号，k＝1,2,…,K；
[0038] 比率计算单元，用于计算数量差值比率：
[0039]
[0040] 卡方值查询单元，用于根据参数K和置信度α，通过查询卡方表获取卡方值χ2(α,K-1)；
[0041] 判断单元，用于判断所述数量差值比率是否大于或等于所述卡方值，得到第一判断结果；
[0042] 地震初至确定单元，用于若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值，则确定在所述地震记录序列的第n个数据处有地震波动到达；若所述第一判断结果表示所述数量差值比率小于所述卡方值，则确定在所述地震记录序列的第n个数据处没有地震波动到达。
[0043] 根据本发明提供的具体
实施例，本发明公开了以下技术效果：本发明的地震波初至检测方法及系统，利用卡方检验的方法识别地震波初至，关注地震记录数据落入某个区间的数量，能够避免依赖于地震数据本身的大小，从而有效避免噪声的影响，提高检测准确度。
附图说明
[0044] 为了更清楚地说明本发明实施例或
现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0045] 图1为本发明地震波初至检测方法实施例的方法
流程图；
[0046] 图2为本发明地震波初至检测系统实施例的系统结构图。
具体实施方式
[0047] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0048] 本发明的目的是提供一种地震波初至检测方法及系统，在低信噪比的情况下有效避免噪声对检测精度的影响，提高检测准确度。
[0049] 为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
[0050] 图1为本发明地震波初至检测方法实施例的方法流程图。
[0051] 本发明将时间窗口内
采样值落入某一区间的次数作为指标，统计次数差值之和，并将此差值之和与卡方值相比较，将比较结果作为地震初至是否到达的标志。
[0052] 参见图1，该地震波初至检测方法，包括：
[0053] 步骤101：获取检波器采集的待检测的地震记录序列。
[0054] 待检测的地震记录序列为p0,p1,…,pN-1，其中N为待检测的地震记录序列的长度，即待检测的地震记录序列所包含的地震数据的数量。
[0055] 步骤102：建立第一时间窗和第二时间窗，所述第一时间窗和所述第二时间窗的容量均为L个数据，其中L为奇数，且L满足31≤L≤41；所述第二时间窗紧随所述第一时间窗。
[0056] 步骤103：对所述地震记录序列的左侧和右侧进行数据填充，使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动至所述第一时间窗的中心对准所述地震记录序列的最后一个数据时，所述第一时间窗和所述第二时间窗内均填充满L个数据。
[0057] 其中，对所述地震记录序列的左侧和右侧进行数据填充，具体包括：将所述地震记录序列的左侧补充至少 个数据；左侧补充的至少 个数据均与所述地震记录序列最左侧的数据相同；将所述地震记录序列的右侧补充至少 个数据；右侧补充的至少个数据均与所述地震记录序列最右侧的数据相同。
[0058] 步骤104：从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动所述第一时间窗和所述第二时间窗，直至所述第一时间窗的中心对准所述地震记录序列的最后一个数据，每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述地震记录序列中的各个数据处是否有地震波动到达。
[0059] 所述第一时间窗内的数据序列为 其中n表示位于所述第一时间窗中心的数据在所述地震记录序列中的序号；pn为位于所述第一时间窗中心的数据；P表示第一时间窗的数据序列，p表示地震记录序列内的各个数据；由于第二时间窗紧随第一时间窗，则所述第二时间窗内的数据序列为
[0060] 其中，每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述地震记录序列中的各个数据处是否有地震波动到达，具体包括：
[0061] 求取所述第一时间窗Pn内所有数据的最大值和最小值，得到取值区间；其中最小值记为p(n)min＝min[Pn]，最大值记为p(n)max＝max[Pn]，则取值区间为[p(n)min,p(n)max]。
[0062] 将所述取值区间划分为K段，得到K段小区间。其中K的取值一般为10；各个小区间的序号为k，划分之后的每个小区间可以表示为：
[0063]
[0064] k＝1,2,…,K
[0065] n＝1,2,…,N
[0066] 计算所述第一时间窗内的数据介于每个所述小区间内的数量F(n)P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F(n)Q(k)。
[0067] 计算数量差值比率：
[0068]
[0069] 根据参数K和置信度α(α一般取为95％)，通过查询卡方表获取卡方值χ2(α,K-1)；
[0070] 判断所述数量差值比率是否大于或等于所述卡方值，得到第一判断结果；
[0071] 若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值，即γ(n)≥χ2(α,K-1)，则确定在所述地震记录序列的第n个数据处有地震波动到达；若所述第一判断结果表示所述数量差值比率小于所述卡方值，即γ(n)＜χ2(α,K-1)，则确定在所述地震记录序列的第n个数据处没有地震波动到达。
[0072] 步骤105：记录所有有地震波动到达的数据在所述地震记录序列中的序号nFB，并计算每个有地震波动到达的数据对应的初至。各个有地震波动到达的数据所对应的地震波初至为t＝(nFB-1)*T。其中T为检波器采样时间间隔，单位为秒。
[0073] 地震波初至检测所面临的最大问题是噪声对检测精度的不良影响。噪声超过一定程度，会造成检测精度极速下降。而本发明则考察地震记录数据落入区间的次数，这种统计方法可以在一定程度上避免噪声的影响。
[0074] 因此，本发明是对相邻的两个数据序列(P和Q)的分布进行统计刻画，如果两者差值之和小于卡方值，则认为这两个数据序列(P和Q)具有相同的分布，没有地震波动到达；如果大于卡方值，则认为这两个数据序列(P和Q)的分布不同，意味着有地震波动到达。利用分布是否相同，可以在一定程度上避免噪声的影响，提高检测精度，并且计算方法十分简单，原理非常明确。
[0075] 图2为本发明地震波初至检测系统实施例的系统结构图。
[0076] 参见图2，该地震波初至检测系统，包括：
[0077] 序列获取模块201，用于获取检波器采集的待检测的地震记录序列。
[0078] 时间窗建立模块202，用于建立第一时间窗和第二时间窗，所述第一时间窗和所述第二时间窗的容量均为L个数据，其中L为奇数；所述第二时间窗紧随所述第一时间窗。
[0079] 数据填充模块203，用于对所述地震记录序列的左侧和右侧进行数据填充，使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动至所述第一时间窗的中心对准所述地震记录序列的最后一个数据时，所述第一时间窗和所述第二时间窗内均填充满L个数据。
[0080] 所述数据填充模块203，具体包括：序列左填充单元和序列右填充单元。其中序列左填充单元，用于将所述地震记录序列的左侧补充至少 个数据；左侧补充的至少个数据均与所述地震记录序列最左侧的数据相同；序列右填充单元，用于将所述地震记录序列的右侧补充至少 个数据；右侧补充的至少 个数据均与所述地震记录序列最右侧的数据相同。
[0081] 地震波初至检测模块204，用于从所述第一时间窗的中心对准所述地震记录序列的第一个数据滑动所述第一时间窗和所述第二时间窗，直至所述第一时间窗的中心对准所述地震记录序列的最后一个数据，每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述地震记录序列中的各个数据处是否有地震波动到达。
[0082] 所述地震波初至检测模块204，具体包括：
[0083] 区间计算单元，用于求取所述第一时间窗内所有数据的最大值和最小值，得到取值区间；区间划分单元，用于将所述取值区间划分为K段，得到K段小区间；区间数据量计算单元，用于计算所述第一时间窗内的数据介于每个所述小区间内的数量F(n)P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F(n)Q(k)；其中n表示位于所述第一时间窗中心的数据在所述地震记录序列中的序号；k为小区间序号，k＝1,2,…,K；比率计算单元，用于计算数量差值比率： 卡方值查询单元，用于根据参数K和2
置信度α，通过查询卡方表获取卡方值χ(α,K-1)；判断单元，用于判断所述数量差值比率是否大于或等于所述卡方值，得到第一判断结果；地震初至确定单元，用于若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值，则确定在所述地震记录序列的第n个数据处有地震波动到达；若所述第一判断结果表示所述数量差值比率小于所述卡方值，则确定在所述地震记录序列的第n个数据处没有地震波动到达。
[0084] 记录与初至计算模块205，用于记录所有有地震波动到达的数据在所述地震记录序列中的序号，并计算每个有地震波动到达的数据对应的初至。
[0085] 根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明的地震波初至检测方法及系统，利用卡方检验的方法识别地震波初至，关注地震记录数据落入某个区间的数量，能够避免依赖于地震数据本身的大小，从而有效避免噪声的影响，提高提高检测准确度。
[0086] 对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
[0087] 本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本
说明书内容不应理解为对本发明的限制。
