[发明专利]一种全自动地震波初至拾取方法及系统

说明书

本发明提出了一种全自动地震波初至拾取算法，通过常速扫描炮集地震数据获取参考初至，然后利用改进的Coppens算法计算初至，通过分析拾取准确的初至与拾取不准确的初至的能量比值的差异，数量化能量比值的变化趋势，计算初至的可信度因子，通过设置可信度阈值约束，可以获取准确的初至数据供后续流程使用。本发明方法不仅拾取初至，并且自动的判断初至拾取的是否准确，可以有效地提高初至拾取效率。通过对不同工区的地震资料进行拾取，结果表明该算法可以快速，准确的拾取初至，极大的降低了人工交互操作的时间，有效的提高了初至拾取的效率，节省大量人力。

技术领域

本发明属于石油地震勘探数据处理领域，具体涉及一种全自动地震波初至拾取方法及系统。

背景技术

在地震勘探中，将从炮点激发并最先到达检波点的地震波称为初至波，并将这个到达时间称为地震波初至时间，简称地震波初至。

在地震资料处理过程中，地震波初至拾取是广泛应用的基础技术，地震波初至拾取准确与否对后续处理过程有极大的影响，对于浅层折射波勘探和层析成像的反演结果的影响尤其明显。

在地震资料处理技术中，地震波初至拾取算法多达上百种，归纳起来大概分为能量比算法、分形维算法、相关法、神经网络法和边缘检测法等几大类。然而大部分的地震波初至拾取算法都是基于首先自动的拾取地震波初至，然后通过人工交互的方式检查拾取的地震波初至是否准确，由人工检查修改后才能供下一处理环节使用。这就使得地震波初至拾取成为地震资料处理中需要进行大量人工交互操作的环节。例如在常用的商业软件Omega和Opencps中都有地震波初至自动拾取这一模块，然而也都都需要进行大量的人工交互操作来检查地震波初至拾取是否准确。

现阶段，随着勘探区域的复杂度以及勘探技术的提高，野外采集的地震数据的数据量激增，一个工区经常会采集多达十几或者几十TB的地震数据。一个处理员一天最多交互检查20GB(100炮，7000道，8000个样点)左右的地震数据，因此，需要20个处理员要持续工作25天才能检查完10TB的地震数据拾取的初至。

因此，现有技术中的震资料处理方法，存在的处理周期长，导致处理效率低，同时需要极大的人工成本等技术问题。

发明内容

为解决现有技术中的震资料处理方法，存在的处理周期长，导致处理效率低，同时需要极大的人工成本等技术问题。本发明提供一种全自动地震波初至拾取方法及系统。

本发明提供一种全自动地震波初至拾取方法，其特征在于，包括以下步骤：

S100；首先扫描炮集数据，获取各道地震数据的参考初至；

S200：基于拾取范围和各道地震数据的参考初至拾取初至；

S300：根据能量比值计算拾取初至的可信度；

S400：根据可信度约束输出拾取准确的初至。

在一个实施例中，在所述步骤S100包括以下子步骤：

S110：用N个常数速度v_k,k＝1，…，N(递增排列)扫描炮集数据，获取N个能量值Q_k，k＝1，…，N；所述能量值Q_k由下式获得：

其中：v_k为常数速度，递增排列，

O(i)为第i道地震数据的偏移距，

S_i(t)为第i道地震数据t时刻的振幅值，

t为地震记录时刻；

S120：计算能量起跳点Q_T，所述能量起跳点由下式获得：