[发明专利]一种地震波时频域初至自动拾取方法

说明书

本发明专利提供了一种地震波时频域初至自动拾取方法。该方法包括对炮集的初至波做时窗截取；通过S变换将地震道转换到时频域，对参考道及目标道的时频谱进行反褶积；在炮域内，对反褶积结果进行叠加生成虚拟时频谱；对虚拟时频谱和参考道时频谱进行褶积；在检波点域内，对褶积结果进行叠加得到超级虚拟时频谱；沿时间轴计算超级虚拟时频谱的能量；采用差分算子获得能量差分曲线；在差分曲线上找到局部最大值作为初至波到时；逐道完成时频域初至自动拾取。本发明抗噪能力强，拾取精度高，不受地震道波形不一致和近地表复杂条件的影响，大大减少人工修改初至波拾取的工作量，能很好地满足各种复杂条件下地震勘探对初至波自动拾取的高精度要求。

技术领域

本发明涉及地球物理勘探地震数据处理技术领域，特别涉及一种高精度地震波初至自动拾取方法。

背景技术

最初地震波初至拾取常采用人机交互的方式，这种传统的手工拾取方式虽然相对精度较高但效率低、费时费力，难以面对采集数据量越来越大的实际生产的需求，目前生产上主要采用初至自动拾取方法。在地震资料处理中，准确拾取初至波是非常重要的，其结果的准确性和质量直接影响到静校正、井间地震勘探和层析成像等后续处理环节的精度。

随着地震勘探工作的不断深入，地震勘探目标区逐渐转向近地表速度和结构变化较为复杂的地区，这种情况下，对初至拾取的精度和质量要求更高。然而，在复杂近地表条件下，地震数据中的初至波形态可能会产生畸变，且受复杂条件影响可能导致各地震道之间的初至波形不一致，此外，实际地震数据中往往存在大量干扰噪声。上述这些问题都使得初至波在高背景噪声情况下难以识别，从而常规初至自动拾取方法难以获得高精度的可靠结果，进而影响地震数据后续相关处理步骤的精度。

为此，亟待需要一种行之有效的初至自动拾取方法，改善常规初至自动拾取方法在各种复杂条件下初至拾取精度低的问题。

发明内容

本发明的目的是提供一种初至自动拾取方法，以提高各种复杂条件包括复杂近地表、各地震道集波形不一致、低信噪比等情况下初至拾取的精度，为后续利用初至波属性来解决问题的各个数据处理步骤提供可靠参数。

为实现上述目的，本发明包括以下步骤：

(1)、对原始炮集数据做时窗截取，提取出含有初至波同相轴的部分；

其中，时窗长度一般为初至波附近的2到3个波长。

进一步地，时窗预估位置的自动确定可通过炮集中各道的偏移距信息除以近地表的大致平均速度来实现。

(2)、将震源位于x处且分别在检波器A和B记录的参考道及目标道通过S变换转换到时频域中，分别表示为S(A|x，τ，f)和S(B|x，τ，f)；

(3)、在时频域对参考道和目标道的时频谱进行反褶积计算；

(4)、在炮域内，对所有位于稳相震源的炮集进行步骤(3)并将所有反褶积计算结果进行叠加求和，生成检波器A和B之间的虚拟时频谱S(B|A，τ，f)^virt；

进一步地，计算公式为：

其中，*代表复共轭，k是波数，N_s是稳相震源的数目，ε为正则化参数。

(5)、在时频域对虚拟时频谱和参考道的时频谱进行褶积计算；

(6)、在检波点域内，对所有位于稳相检波点的参考道进行步骤(5)并将所有褶积计算结果进行叠加求和，得到震源位于x处且在检波点B记录的超级虚拟时频谱S(B|x，τ，f)^super；

进一步地，计算公式为：

其中，N_g是稳相检波器的数目。