[发明专利]初至时间拾取方法

说明书

公开了一种初至时间拾取方法。该方法可以包括：基于原始地震道，计算初至搜索中心时间，获得搜索窗口半径，进而确定初至搜索范围；基于初至搜索中心时间与原始地震道的数据，计算原始地震道的稳定因子；在原始地震道中的初至搜索范围内，针对每个样点，分别在样点的前后划分一个目标时窗，进而计算样点的地震属性值；基于每个样点的地震属性值与稳定因子，计算每个样点对应的能量比；提取具有最大的能量比的样点所对应的样点时间作为初至时间。本发明通过新的地震属性道与自适应求取稳定因子，可以精确高效地拾取初至时间。

技术领域

本发明涉及油气及煤层气地震勘探与开发领域，更具体地，涉及一种初至时间拾取方法。

背景技术

静校正是提高叠加剖面信噪比和垂向分辨率的一项关键技术，常规折射波静校正处理是建立在水平均匀层状介质假设前提下的，当勘探区域的地表起伏大，近地表厚度、速度变化复杂时，折射静校正就无法适用了。层析静校正能反演任意介质的近地表速度模型，是目前解决复杂地表静校正最理想的静校正方法，也是最常用的方法。而初至拾取又是层析静校正处理的第一步，它的准确与否会直接影响到层析反演结果。最初的地震波初至拾取是以人机交互的方式人工拾取的，因此初至拾取就成为一项费时而繁琐的工作。近年来，出现了许多自动初至拾取方法，例如相关法、能量比法、最大振幅法、分形维法以及神经网络法等，其中，能量比法以其简单实用而得到了广泛应用。但是在复杂近地表条件下，由于初至波能量特征、波形特征、相位特征变化大，又常常受噪声影响，初至波自动拾取方法的精度和效率仍然无法满足要求。如何进一步提高初至时间拾取的准确性、稳定性和抗干扰能力仍然需要进一步探索研究。

公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

发明内容

本发明提出了一种初至时间拾取方法，其能够通过新的地震属性道与自适应求取稳定因子，可以精确高效地拾取初至时间。

根据本发明的一方面，提出了一种初至时间拾取方法。所述方法可以包括：基于原始地震道，计算初至搜索中心时间，获得搜索窗口半径，进而确定初至搜索范围；基于所述初至搜索中心时间与所述原始地震道的数据，计算所述原始地震道的稳定因子；在所述原始地震道中的所述初至搜索范围内，针对每个样点，分别在所述样点的前后划分一个目标时窗，进而计算所述样点的地震属性值；基于每个样点的地震属性值与所述稳定因子，计算每个样点对应的能量比；提取具有最大的能量比的样点所对应的样点时间作为初至时间。

优选地，所述初至搜索中心时间为：

其中，m表示初至搜索中心时间，x表示原始地震道的偏移距，v表示地下近地表速度的估计值。

优选地，所述稳定因子为：

其中，β表示稳定因子，s_i表示第i个样点的样点值，i表示样点序号，N表示样点个数。

优选地，所述能量比为：

其中，ER(i)表示第i个样点对应的能量比，pdist(j)表示第j个样点的地震属性值，i1表示第i个样点的前一个目标时窗的起始样点的样点序号，i2表示第i个样点的后一个目标时窗的终止样点的样点序号。

优选地，所述地震属性值为：

其中，pdist1(i)表示第i个样点的地震属性值，N_c表示每个目标时窗包含的样点数，s_l表示第l个样点的样点值，s_l+k表示第l+k个样点的样点值，k为整数。

优选地，所述地震属性值为：