CN108035708A - 一种去除地层界面反射波的方法及装置 - Google Patents

Abstract

本申请实施例公开了一种去除地层界面反射波的方法及装置。所述方法提供有目的工区中指定钻井位置处的偶极横波数据；所述方法包括：确定所述偶极横波数据对应的偶极横波时差；对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据；基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据。本申请提供的技术方案，可以实现对地层界面反射波的有效去除。

Description

一种去除地层界面反射波的方法及装置

技术领域

本申请涉及复杂岩性储层油气勘探技术领域，特别涉及一种去除地层界面反射波的方法及装置。

背景技术

复杂岩性储层具有强烈的非均质性和各向异性，目前常规测井系列和微电阻率成像测井系列仅能准确评价井周两米范围内储层发育情况，井周储层发育较差并不能完全表示井旁储层发育较差。但是，偶极横波远探测技术能够反映井旁3m-30m范围内地层或缝洞的声阻抗异常，从而可以有效发现井壁外围存在的缝洞体，可以为钻井井位设计和开采方案的缺点提供科学依据。

唐晓明于2004年首次把偶极子用于反射波测井中，并提出了偶极横波法的远探测声波成像方法，该方法采用偶极子声源在井中向外发射横波，并在井中接收由井外地质体反射回来的横波。2009年唐晓明利用四分量偶极数据对井旁裂缝和盐丘内部的精细构造进行了成像，并取得了显著的应用效果。

目前偶极横波远探测技术的实施过程主要分为三个步骤：第一步是横波反射波提取，第二步是SH、SV波分离，第三步是偏移成像。其中，横波反射波提取是一个非常关键的处理环节，该步骤主要是做两方面的工作：一是尽可能压制井孔直达波信号，二是滤除掉无效反射波信号，比如地层界面反射波等。

当地层弹性性质或井孔尺寸在轴向上变化较大时，便会产生幅度很强的地层界面反射波，地层界面反射波在波形中持续时间长，覆盖面积大，可能淹没有效横波反射波信号。因此，在获取有效横波反射波之前，需要剔除地层界面反射波。迄今，国内外学者针对横波反射波提取，主要是在压制井孔直达波方面做工作，比如利用F-K滤波来滤除井孔直达波信号，而对地层界面反射波消除方面的研究较少。为此，有必要研究一种去除地层界面反射波的方法，以实现对地层界面反射波的有效去除。

发明内容

本申请实施例的目的是提供一种去除地层界面反射波的方法及装置，以实现对地层界面反射波的有效去除。

为解决上述技术问题，本申请实施例提供一种去除地层界面反射波的方法及装置是这样实现的：

一种去除地层界面反射波的方法，提供有目的工区中指定钻井位置处的偶极横波数据；所述方法包括：

确定所述偶极横波数据对应的偶极横波时差；

对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据；

基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据。

优选方案中，在确定所述偶极横波数据对应的偶极横波时差之前，所述方法还包括：

对所述偶极横波数据进行预处理，得到预处理后的偶极横波数据；

相应的，所述确定所述偶极横波数据对应的偶极横波时差，包括：确定所述预处理后的偶极横波数据对应的偶极横波时差；

所述对所述偶极横波数据进行垂直中值滤波，包括：对所述预处理后的偶极横波数据进行垂直中值滤波。

优选方案中，所述对所述偶极横波数据进行预处理，包括：

对所述偶极横波收进行增益恢复处理，得到增益恢复后的偶极横波数据；

对所述增益恢复后的偶极横波数据进行延迟恢复处理，得到所述预处理后的偶极横波数据。

优选方案中，所述对所述预处理后的偶极横波数据进行垂直中值滤波，包括：

沿与水平面垂直的同相轴方向对所述预处理后的偶极横波数据进行中值滤波。

优选方案中，所述根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据，包括：

将所述预处理后的偶极横波数据减去所述初始偶极横波直达波数据，得到所述指定钻井位置处的初始偶极横波反射波数据。

优选方案中，所述确定所述预处理后的偶极横波数据对应的偶极横波时差，包括：

采用时间-慢度相关方法确定所述预处理后的偶极横波数据对应的偶极横波时差；具体地，在时间和慢度两个维度上计算所述预处理后的偶极横波数据对应的相关函数；将所述相关函数的最大函数值对应的慢度作为所述偶极横波时差。

优选方案中，采用下述公式表征所述相关函数：

其中，Corr(DT,Time)表示所述相关函数的函数值；X_m[t+DT(m-1)·RRSP]表示所述预处理后的偶极横波数据中第m道声波波形数据中时间为t+DT(m-1)·RRSP处的振幅值，t表示时间，DT表示慢度，RRSP表示偶极横波测井仪器中两个相邻接收器之间的距离，N表示所述预处理后的偶极横波数据中声波波形数据的道数；Time表示预设对比时窗的起始时间，CWLENGTH表示所述预设对比时窗的窗长；符号|·|为取绝对值运算符。

优选方案中，所述地层界面反射波数据包括：地层界面上行反射波数据和地层界面下行反射波数据；所述基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据，包括：

基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜上行中值滤波，得到倾斜上行中值滤波后的初始偶极横波反射波数据；

将所述初始偶极横波反射波数据减去所述倾斜上行中值滤波后的初始偶极横波反射波数据，得到第一偶极横波反射波数据；所述第一偶极横波反射波数据表示去除所述地层界面上行反射波数据后的初始偶极横波反射波数据；

基于所述偶极横波时差，对所述第一偶极横波反射波数据进行倾斜下行中值滤波，得到倾斜下行中值滤波后的第一偶极横波反射波数据；

将所述第一偶极横波反射波数据减去所述倾斜下行中值滤波后的第一偶极横波反射波数据，得到第二偶极横波反射波数据；所述第二偶极横波反射波数据表示去除所述地层界面上行反射波数据和所述地层界面下行反射波数据后的初始偶极横波反射波数据。

优选方案中，所述基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜上行中值滤波，得到倾斜上行中值滤波后的偶极横波反射波数据，包括：

将斜率为所述偶极横波时差的倒数的两倍的同相轴作为倾斜上行同相轴；

沿所述倾斜上行同相轴对所述初始偶极横波反射波数据进行倾斜上行中值滤波，得到所述倾斜上行中值滤波后的初始偶极横波反射波数据。

优选方案中，所述基于所述偶极横波时差，对所述第一偶极横波反射波数据进行倾斜下行中值滤波，得到倾斜下行中值滤波后的第一偶极横波反射波数据，包括：

将斜率为所述偶极横波时差的倒数的两倍的负数的同相轴作为倾斜下行同相轴；

沿所述倾斜下行同相轴对所述第一偶极横波反射波数据进行倾斜下行中值滤波，得到所述倾斜下行中值滤波后的第一偶极横波反射波数据。

一种去除地层界面反射波的装置，所述装置提供目的工区中指定钻井位置处的偶极横波数据；所述装置包括：时差确定模块、垂直中值滤波模块和倾斜中值滤波模块；其中，

所述时差确定模块，用于确定所述偶极横波数据对应的偶极横波时差；

所述垂直中值滤波模块，用于对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据；

所述倾斜中值滤波模块，用于基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据。

本申请实施例提供了一种去除地层界面反射波的方法及装置，可以确定所述偶极横波数据对应的偶极横波时差；可以对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据；可以基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据。如此，可以实现对地层界面反射波的有效去除。

附图说明

为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1是本申请一种去除地层界面反射波的方法实施例的流程图；

图2是本申请实施例中采用本申请方法去除地层界面反射波的剖面示意图；

图3是本申请实施例中垂直中值滤波过程的示意图；

图4是本申请实施例中倾斜上行中值滤波过程的示意图；

图5是本申请实施例中倾斜下行中值滤波过程的示意图；

图6是本申请去除地层界面反射波的装置实施例的流程图。

具体实施方式

本申请实施例提供一种去除地层界面反射波的方法及装置。

为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。

本申请提供了一种去除地层界面反射波的方法。所述去除地层界面反射波的方法可以有目的工区中指定钻井位置处的偶极横波数据。

在本实施方式中，所述目的工区可以是存在复杂岩性储层的工区。所述目的工区中指定钻井位置可以指通过在所述目的工区中复杂岩性储层中已经开设的钻井，能够采集测井数据的位置。

在本实施方式中，可以通过对在所述钻井位置处获取的常规测井数据、电成像测井数据和气测录井数据进行综合分析，从偶极横波测井数据中提取出在偶极横波测井过程中易发生地层界面反射的目的层段处的偶极横波数据。

图1是本申请一种去除地层界面反射波的方法实施例的流程图。如图1所示，所述去除地层界面反射波的方法，包括以下步骤。

步骤S101：确定所述偶极横波数据对应的偶极横波时差。

在本实施方式中，在确定所述偶极横波数据对应的偶极横波时差之前，所述方法还可以包括：对所述偶极横波数据进行预处理，得到预处理后的偶极横波数据。

相应的，所述确定所述偶极横波数据对应的偶极横波时差，可以包括：可以确定所述预处理后的偶极横波数据对应的偶极横波时差。具体地，在井下采集偶极横波数据过程中，为了保证偶极横波数据中的振幅值处理检测的最高精度，可以对所述偶极横波收进行增益恢复处理，得到增益恢复后的偶极横波数据。同时为了减少记录数据的存储量，通常不采集在首波到达之前的指定时间段的偶极横波数据，需要对所述增益恢复后的偶极横波数据进行延迟恢复处理，得到所述预处理后的偶极横波数据。

在本实施方式中，所述增益恢复处理的具体过程是将所述偶极横波数据除以增益参数，以将所述偶极横波数据中的振幅值从小数放大至整数，从而便于存储记录。其中，可以采用下述公式表征所述增益参数：

AGN＝10^GN*0.05

其中，AGN表示所述增益参数，GN表示增益指数。通常所述增益指数的取值范围为-69～-21分贝。

在本实施方式中，所述延迟恢复处理的具体过程是在所述偶极横波数据的零时刻之前补零，即将所述指定时间段的偶极横波数据的振幅值作为零，并补充至所述偶极横波数据的零时刻之前，以获取准确时间的偶极横波数据。

在本实施方式中，确定所述预处理后的偶极横波数据对应的偶极横波时差，具体可以包括，可以采用时间-慢度相关方法确定所述预处理后的偶极横波数据对应的偶极横波时差。具体地，可以在时间和慢度两个维度上计算所述预处理后的偶极横波数据对应的相关函数。可以将所述相关函数的最大函数值对应的慢度作为所述偶极横波时差。

在本实施方式中，可以采用下述公式表征所述相关函数：

其中，Corr(DT,Time)表示所述相关函数的函数值；X_m[t+DT(m-1)·RRSP]表示所述预处理后的偶极横波数据中第m道声波波形数据中时间为t+DT(m-1)·RRSP处的振幅值，t表示时间，DT表示慢度，RRSP表示偶极横波测井仪器中两个相邻接收器之间的距离，N表示所述预处理后的偶极横波数据中声波波形数据的道数；Time表示预设对比时窗的起始时间，CWLENGTH表示所述预设对比时窗的窗长；符号|·|为取绝对值运算符。

步骤S102：对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据。

在本实施方式中，在确定所述偶极横波数据对应的偶极横波时差之前，所述方法还可以包括：对所述偶极横波数据进行预处理，得到预处理后的偶极横波数据。相应的，对所述偶极横波数据进行垂直中值滤波，具体可以包括：可以对所述预处理后的偶极横波数据进行垂直中值滤波。

在本实施方式中，对所述预处理后的偶极横波数据进行垂直中值滤波，具体可以包括：沿与水平面垂直的同相轴方向对所述预处理后的偶极横波数据进行中值滤波。

在本实施方式中，根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据，具体可以包括：可以将所述预处理后的偶极横波数据减去所述偶极横波直达波数据，得到所述指定钻井位置处的初始偶极横波反射波数据。

在本实施方式中，可以选取深度上相邻的三道偶极横波数据进行所述垂直中值滤波。

步骤S103：基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据。

在本实施方式中，所述地层界面反射波数据可以包括：地层界面上行反射波数据和地层界面下行反射波数据。

在本实施方式中，基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据，具体可以包括，可以基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜上行中值滤波，可以得到倾斜上行中值滤波后的初始偶极横波反射波数据。可以将所述初始偶极横波反射波数据减去所述倾斜上行中值滤波后的初始偶极横波反射波数据，得到第一偶极横波反射波数据。所述第一偶极横波反射波数据表示去除所述地层界面上行反射波数据后的初始偶极横波反射波数据。可以基于所述偶极横波时差，对所述第一偶极横波反射波数据进行倾斜下行中值滤波，得到倾斜下行中值滤波后的第一偶极横波反射波数据。可以将所述第一偶极横波反射波数据减去所述倾斜下行中值滤波后的第一偶极横波反射波数据，得到第二偶极横波反射波数据。所述第二偶极横波反射波数据表示去除所述地层界面上行反射波数据和所述地层界面下行反射波数据后的初始偶极横波反射波数据。

在本实施方式中，基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜上行中值滤波，得到倾斜上行中值滤波后的偶极横波反射波数据，具体可以包括，可以将斜率为所述偶极横波时差的倒数的两倍的同相轴作为倾斜上行同相轴。可以沿所述倾斜上行同相轴对所述初始偶极横波反射波数据进行倾斜上行中值滤波，得到所述倾斜上行中值滤波后的偶极横波反射波数据。

在本实施方式中，基于所述偶极横波时差，对所述第一偶极横波反射波数据进行倾斜下行中值滤波，得到倾斜下行中值滤波后的第一偶极横波反射波数据，具体可以包括，可以将斜率为所述偶极横波时差的倒数的两倍的负数的同相轴作为倾斜下行同相轴。可以沿所述倾斜下行同相轴对所述第一偶极横波反射波数据进行倾斜下行中值滤波，得到所述倾斜下行中值滤波后的第一偶极横波反射波数据。

在本实施方式中，可以选取深度上相邻的三道偶极横波数据进行所述倾斜上行中值滤波和所述倾斜下行中值滤波。

例如，图2是本申请实施例中采用本申请方法去除地层界面反射波的剖面示意图。图2中(a)、(b)、(c)和(d)分别为预处理后的偶极横波数据、初始偶极横波反射波数据、第一偶极横波反射波数据和第二偶极横波反射波数据的剖面示意图。从图2中可以看出，通过去除图2中(a)中的偶极横波直达波数据，可以得到如图2中(b)所示的初始偶极横波反射波数据，接着通过去除图2中(b)中的地层界面上行反射波数据，可以得到如图2中(c)所示的第一偶极横波反射波数据，最后通过去除图2中(c)中的地层界面下行反射波数据，可以得到如图2中(d)所示的第二偶极横波反射波数据，从而可以有效去除地层界面反射波。

图3是本申请实施例中垂直中值滤波过程的示意图。图4是本申请实施例中倾斜上行中值滤波过程的示意图。图5是本申请实施例中倾斜下行中值滤波过程的示意图。图3至图5中的横坐标和纵坐标分别是时间和深度。图3至图5中实线表示井孔直达波，点划线表示地层界面上行反射波，较密集的虚线表示地层界面下行反射波，以及较稀疏的虚线表示有效反射波。如图3至图5所示，通过图3中的垂直中值滤波过程，可以去除如图3中(a)的预处理后的偶极横波数据中的偶极横波直达波数据，即如图3中(b)所示的井孔直达波数据，得到如图3中(c)所示的初始偶极横波反射波数据。通过图4中的倾斜上行中值滤波过程，可以去除如图4中(a)的初始偶极横波反射波数据中的地层界面上行反射波数据，即如图4中(b)所示的地层界面上行反射波数据，得到如图4中(c)所示的第一偶极横波反射波数据。通过图5中的倾斜下行中值滤波过程，可以去除如图5中(a)所示的第一偶极横波反射波数据中的地层界面下行反射波数据，即如图5中(b)所示的地层界面下行反射波，从而得到如图5中(c)所示的有效反射波数据。

所述去除地层界面反射波的方法实施例，可以确定所述偶极横波数据对应的偶极横波时差；可以对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据；可以基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据。如此，可以实现对地层界面反射波的有效去除。

图6是本申请去除地层界面反射波的装置实施例的流程图。所述去除地层界面反射波的装置提供目的工区中指定钻井位置处的偶极横波数据。如图6所示，所述去除地层界面反射波的装置可以包括：时差确定模块100、垂直中值滤波模块200和倾斜中值滤波模块300。

所述时差确定模块100，可以用于确定所述偶极横波数据对应的偶极横波时差。

所述垂直中值滤波模块200，可以用于对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据。

所述倾斜中值滤波模块300，可以用于基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据。

所述去除地层界面反射波的装置实施例与所述去除地层界面反射波的方法实施例相对应，可以实现去除地层界面反射波的方法实施例的技术方案，并取得方法实施例的技术效果。

在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable GateArray，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware DescriptionLanguage)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(RubyHardware Description Language)等，目前最普遍使用的是VHDL(Very-High-SpeedIntegrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。

本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。

上述实施例阐明的装置、模块，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。

为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。

通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中，内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括短暂电脑可读媒体(transitory media)，如调制的数据信号和载波。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于装置实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。

本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。

虽然通过实施例描绘了本申请，本领域普通技术人员知道，本申请有许多变形和变化而不脱离本申请的精神，希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (11)

1.一种去除地层界面反射波的方法，其特征在于，提供有目的工区中指定钻井位置处的偶极横波数据；所述方法包括：

确定所述偶极横波数据对应的偶极横波时差；

对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据；

基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据。

2.根据权利要求1所述的一种去除地层界面反射波的方法，其特征在于，在确定所述偶极横波数据对应的偶极横波时差之前，所述方法还包括：

对所述偶极横波数据进行预处理，得到预处理后的偶极横波数据；

相应的，所述确定所述偶极横波数据对应的偶极横波时差，包括：确定所述预处理后的偶极横波数据对应的偶极横波时差；

所述对所述偶极横波数据进行垂直中值滤波，包括：对所述预处理后的偶极横波数据进行垂直中值滤波。

3.根据权利要求2所述的一种去除地层界面反射波的方法，其特征在于，所述对所述偶极横波数据进行预处理，包括：

对所述偶极横波收进行增益恢复处理，得到增益恢复后的偶极横波数据；

对所述增益恢复后的偶极横波数据进行延迟恢复处理，得到所述预处理后的偶极横波数据。

4.根据权利要求2所述的一种去除地层界面反射波的方法，其特征在于，所述对所述预处理后的偶极横波数据进行垂直中值滤波，包括：

沿与水平面垂直的同相轴方向对所述预处理后的偶极横波数据进行中值滤波。

5.根据权利要求4所述的一种去除地层界面反射波的方法，其特征在于，所述根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据，包括：

将所述预处理后的偶极横波数据减去所述初始偶极横波直达波数据，得到所述指定钻井位置处的初始偶极横波反射波数据。

6.根据权利要求2所述的一种去除地层界面反射波的方法，其特征在于，所述确定所述预处理后的偶极横波数据对应的偶极横波时差，包括：

采用时间-慢度相关方法确定所述预处理后的偶极横波数据对应的偶极横波时差；具体地，在时间和慢度两个维度上计算所述预处理后的偶极横波数据对应的相关函数；将所述相关函数的最大函数值对应的慢度作为所述偶极横波时差。

7.根据权利要求6所述的一种去除地层界面反射波的方法，其特征在于，采用下述公式表征所述相关函数：

<mrow> <mi>C</mi> <mi>o</mi> <mi>r</mi> <mi>r</mi> <mrow> <mo>(</mo> <mi>D</mi> <mi>T</mi> <mo>,</mo> <mi>T</mi> <mi>i</mi> <mi>m</mi> <mi>e</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mi>T</mi> <mi>i</mi> <mi>m</mi> <mi>e</mi> </mrow> <mrow> <mo>(</mo> <mi>T</mi> <mi>i</mi> <mi>m</mi> <mi>e</mi> <mo>+</mo> <mi>C</mi> <mi>W</mi> <mi>L</mi> <mi>E</mi> <mi>N</mi> <mi>G</mi> <mi>T</mi> <mi>H</mi> <mo>)</mo> </mrow> </msubsup> <mo>|</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>X</mi> <mi>m</mi> </msub> <mo>&amp;lsqb;</mo> <mi>t</mi> <mo>+</mo> <mi>D</mi> <mi>T</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mi>R</mi> <mi>R</mi> <mi>S</mi> <mi>P</mi> <mo>&amp;rsqb;</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> <mi>d</mi> <mi>t</mi> </mrow> <mrow> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mi>T</mi> <mi>i</mi> <mi>m</mi> <mi>e</mi> </mrow> <mrow> <mo>(</mo> <mi>T</mi> <mi>i</mi> <mi>m</mi> <mi>e</mi> <mo>+</mo> <mi>C</mi> <mi>W</mi> <mi>L</mi> <mi>E</mi> <mi>N</mi> <mi>G</mi> <mi>T</mi> <mi>H</mi> <mo>)</mo> </mrow> </msubsup> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mo>|</mo> <msub> <mi>X</mi> <mi>m</mi> </msub> <mo>&amp;lsqb;</mo> <mi>t</mi> <mo>+</mo> <mi>D</mi> <mi>T</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mi>R</mi> <mi>R</mi> <mi>S</mi> <mi>P</mi> <mo>&amp;rsqb;</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> </mrow>

其中，Corr(DT,Time)表示所述相关函数的函数值；X_m[t+DT(m-1)·RRSP]表示所述预处理后的偶极横波数据中第m道声波波形数据中时间为t+DT(m-1)·RRSP处的振幅值，t表示时间，DT表示慢度，RRSP表示偶极横波测井仪器中两个相邻接收器之间的距离，N表示所述预处理后的偶极横波数据中声波波形数据的道数；Time表示预设对比时窗的起始时间，CWLENGTH表示所述预设对比时窗的窗长；符号|·|为取绝对值运算符。

8.根据权利要求1所述的一种去除地层界面反射波的方法，其特征在于，所述地层界面反射波数据包括：地层界面上行反射波数据和地层界面下行反射波数据；所述基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据，包括：

基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜上行中值滤波，得到倾斜上行中值滤波后的初始偶极横波反射波数据；

将所述初始偶极横波反射波数据减去所述倾斜上行中值滤波后的初始偶极横波反射波数据，得到第一偶极横波反射波数据；所述第一偶极横波反射波数据表示去除所述地层界面上行反射波数据后的初始偶极横波反射波数据；

基于所述偶极横波时差，对所述第一偶极横波反射波数据进行倾斜下行中值滤波，得到倾斜下行中值滤波后的第一偶极横波反射波数据；

将所述第一偶极横波反射波数据减去所述倾斜下行中值滤波后的第一偶极横波反射波数据，得到第二偶极横波反射波数据；所述第二偶极横波反射波数据表示去除所述地层界面上行反射波数据和所述地层界面下行反射波数据后的初始偶极横波反射波数据。

9.根据权利要求8所述的一种去除地层界面反射波的方法，其特征在于，所述基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜上行中值滤波，得到倾斜上行中值滤波后的偶极横波反射波数据，包括：

将斜率为所述偶极横波时差的倒数的两倍的同相轴作为倾斜上行同相轴；

沿所述倾斜上行同相轴对所述初始偶极横波反射波数据进行倾斜上行中值滤波，得到所述倾斜上行中值滤波后的初始偶极横波反射波数据。

10.根据权利要求8所述的一种去除地层界面反射波的方法，其特征在于，所述基于所述偶极横波时差，对所述第一偶极横波反射波数据进行倾斜下行中值滤波，得到倾斜下行中值滤波后的第一偶极横波反射波数据，包括：

将斜率为所述偶极横波时差的倒数的两倍的负数的同相轴作为倾斜下行同相轴；

沿所述倾斜下行同相轴对所述第一偶极横波反射波数据进行倾斜下行中值滤波，得到所述倾斜下行中值滤波后的第一偶极横波反射波数据。

11.一种去除地层界面反射波的装置，其特征在于，所述装置提供目的工区中指定钻井位置处的偶极横波数据；所述装置包括：时差确定模块、垂直中值滤波模块和倾斜中值滤波模块；其中，

所述时差确定模块，用于确定所述偶极横波数据对应的偶极横波时差；

所述垂直中值滤波模块，用于对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据；

所述倾斜中值滤波模块，用于基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据。