本申请实施例提供一种拾取初至波的方法和装置。该方法包括对采集的地震数据进行卷积变换，得到卷积后的地震数据；利用所述采集的地震数据和所述卷积后的地震数据计算得到地震道中采样点的包络；根据所述地震道中采样点预设位置的参考采样点的包络确定所述地震道中采样点的瞬时振幅；利用地震道中采样点的瞬时振幅和预设能量比值公式计算所述地震道中采样点的能量比值；计算地震道中采样点的能量比值与所述采样点后一个采样点的能量比值的差值；比较所述差值的绝对值的大小；将数值最大的绝对值对应的采样点作为地震道的初至波的样点。利用本申请实施例提供的技术方案可以大大提高拾取初至波的准确率。
【技术实现步骤摘要】
本专利技术涉及地球物理勘探数据处理
，尤其涉及一种拾取初至波的方法和 装置。
技术介绍
地震勘探中的静校正处理技术可以用来解决地震波在近地表介质旅行时间变化 影响反射波叠加成像的问题。在静校正处理技术中，初至是地震资料的一个重要属性。在 地震资料处理过程中，通常可以利用初至的属性解决地震数据的静校正问题。当地震发生 后，由于各种地震波的传播速度不同，传播到观测点的时间也就有先后。地震观测点最先接 收到的波就可以称为初至波。 现有技术中拾取初至波的方法主要包括边界能量检测的方法拾取初至波。具体 的，可以通过对采集的地震数据进行能量检测，将道集中能量高的采样点对应的波作为所 述道集的初至波。但近年来随着可控震源的广泛应用，采集的地震数据一般具有大数据量、 低信噪比的特点。这就导致现有技术中通过边界能量检测的方法从大数据量、低信噪比的 地震数据中拾取初至波时，地震数据本身具有大量能量干扰信息，比如续至波等。从而降低 了初至波拾取的准确率。 因此，现有技术中亟需一种有效的拾取初至波的方法，从而为后续利用初至波的 属性解决地震数据的静校正问题提供有效的依据。
技术实现思路
本申请的目的是提供一种拾取初至波的方法和装置，以提高拾取初至波的准确 率，为后续利用初至波的属性解决地震数据的静校正问题提供有效的依据。 为了实现上述目的，本申请提供了一种拾取初至波的方法，所述方法包括： 对采集的地震数据进行卷积变换，得到卷积后的地震数据； 利用所述采集的地震数据和所述卷积后的地震数据计算得到地震道中采样点的 包络； 根据所述地震道中采样点预设位置的参考采样点的包络确定所述地震道中采样 点的瞬时振幅； 利用地震道中采样点的瞬时振幅和预设能量比值公式计算所述地震道中采样点 的能量比值； 计算地震道中采样点的能量比值与所述采样点后一个采样点的能量比值的差 值； 比较所述差值的绝对值的大小； 将数值最大的绝对值对应的采样点作为所述地震道的初至波的样点。 在一个优选的实施例中，所述利用所述采集的地震数据和所述卷积后的地震数据 计算得到地震道中采样点的包络的计算公式如下： ， 上式中，民代表所述地震道中第r个采样点的包络；KrSn，n代表所述地震 道的采样点数量；&代表所述地震道中第r个采样点的地震数据；笔代表所述地震道中第 r个采样点的卷积后的地震数据。 在一个优选的实施例中，所述根据所述地震道中采样点预设位置的参考采样点的 包络确定所述地震道中采样点的瞬时振幅包括： 获取所述地震道的时窗采样点数量，将所述时窗采样点数量作为预设数量； 利用所述地震道中第一采样点前的预设数量的采样点的包络计算得到所述第一 采样点的能量均值和所述第一采样点的能量标准方差； 利用所述第一采样点的能量均值和所述第一采样点的能量标准方差计算得到所 述地震道中采样点的瞬时振幅。 在一个优选的实施例中，所述利用所述地震道中第一采样点前的预设数量的采样 点的包络计算得到所述第一采样点的能量均值和所述第一采样点的能量标准方差的计算 公式如下： 上式中，Ajf表所述第一米样点的能量均值；1 <g<n，g彡W，g代表所述第一米 样点的点号，n代表所述地震道的采样点数量；民代表所述地震道中第r个采样点的包络； W代表时窗米样点数量，为正整数； 上式中，Dg代表所述第一采样点的能量标准方差；E^代表所述地震道中第r个采 样点的包络；W代表时窗采样点数量，为正整数；Ag代表所述第一采样点的能量均值。 在一个优选的实施例中，所述利用所述第一采样点的能量均值和所述第一采样点 的能量标准方差计算得到所述地震道中采样点的瞬时振幅的计算公式如下： 上式中，Hg代表所述第一采样点的瞬时振幅;Ag代表所述第一采样点的能量均值； Dg代表所述第一采样点的能量标准方差。 在一个优选的实施例中，所述预设能量比值公式如下： 上式中，民代表所述地震道中第i个采样点的能量比值；1彡i彡n，i彡W，n代 表所述地震道的采样点数量；氏代表所述地震道中第r采样点的瞬时振幅；Hi代表所述第 一米样点的瞬时振幅；W代表时窗米样点数量，为正整数。 -种拾取初至波的装置，所述装置包括： 卷积变换模块，用于对采集的地震数据进行卷积变换，得到卷积后的地震数据； 第一计算模块，用于利用所述采集的地震数据和所述卷积后的地震数据计算得到 地震道中采样点的包络； 确定模块，用于根据所述地震道中采样点预设位置的参考采样点的包络确定所述 地震道中采样点的瞬时振幅； 第二计算模块，用于利用地震道中采样点的瞬时振幅和预设能量比值公式计算所 述地震道中采样点的能量比值； 第三计算模块，用于计算地震道中采样点的能量比值与所述采样点后一个采样点 的能量比值的差值； 比较模块，用于比较所述差值的绝对值的大小； 数据处理模块，用于将数值最大的绝对值对应的采样点作为所述地震道的初至波 的样点。 在一个优选的实施例中，所述第一计算模块的计算公式如下： 上式中，民代表所述地震道中第r个采样点的包络；1 <r$n，n代表所述地震道 的采样点数量；\代表所述地震道中第r个采样点的地震数据；毛代表所述地震道中第r个 采样点的卷积后的地震数据。 在一个优选的实施例中，所述确定模块包括： 数据获取单元，用于获取所述地震道的时窗采样点数量，将所述时窗采样点数量 作为预设数量； 第一计算单元，用于利用所述地震道中第一采样点前的预设数量的采样点的包络 计算得到所述第一采样点的能量均值和所述第一采样点的能量标准方差； 第二计算单元，用于利用所述第一采样点的能量均值和所述第一采样点的能量标 准方差计算得到所述地震道中采样点的瞬时振幅。 在一个优选的实施例中，所述第一计算单元的计算公式可以如下： 上式中，Ajf表所述第一米样点的能量均值；1 <g<n，g彡W，g代表所述第一米 样点的点号，n代表所述地震道的采样点数量；民代表所述地震道中第r个采样点的包络； W代表时窗米样点数量，为正整数； 上式中，Dg代表所述第一采当前第1页1 2 3 4 本文档来自技高网...
【技术保护点】
一种拾取初至波的方法，其特征在于，所述方法包括：对采集的地震数据进行卷积变换，得到卷积后的地震数据；利用所述采集的地震数据和所述卷积后的地震数据计算得到地震道中采样点的包络；根据所述地震道中采样点预设位置的参考采样点的包络确定所述地震道中采样点的瞬时振幅；利用地震道中采样点的瞬时振幅和预设能量比值公式计算所述地震道中采样点的能量比值；计算地震道中采样点的能量比值与所述采样点后一个采样点的能量比值的差值；比较所述差值的绝对值的大小；将数值最大的绝对值对应的采样点作为所述地震道的初至波的样点。
【技术特征摘要】
【专利技术属性】
技术研发人员：许银坡，杨海申，邹雪峰，唐东磊，何永清，徐丽军，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京;11
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。