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

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

Method and device for removing reflection wave from formation interface

The embodiment of the application discloses a method and a device for removing reflection waves at the formation interface. The method provides the dipole transverse wave data at the designated drilling position at the target area; the method includes: determining the dipole transverse wave time difference corresponding to the dipole transverse wave data; vertical median filtering of the dipole transverse wave data, and obtaining the direct wave data of the even pole transverse wave at the specified drilling position; The data of the initial dipole shear wave reflection wave at the specified drilling position is determined by the direct wave data of the extreme transverse wave; based on the dipole transverse wave time difference, the initial dipole shear wave reflected wave data is tilted median filtering to remove the reflection wave data at the stratum interface at the specified drilling position. The technical proposal provided by this application can effectively remove the reflected wave from the formation interface.

【技术实现步骤摘要】
一种去除地层界面反射波的方法及装置
本申请涉及复杂岩性储层油气勘探
，特别涉及一种去除地层界面反射波的方法及装置。
技术介绍
复杂岩性储层具有强烈的非均质性和各向异性，目前常规测井系列和微电阻率成像测井系列仅能准确评价井周两米范围内储层发育情况，井周储层发育较差并不能完全表示井旁储层发育较差。但是，偶极横波远探测技术能够反映井旁3m-30m范围内地层或缝洞的声阻抗异常，从而可以有效发现井壁外围存在的缝洞体，可以为钻井井位设计和开采方案的缺点提供科学依据。唐晓明于2004年首次把偶极子用于反射波测井中，并提出了偶极横波法的远探测声波成像方法，该方法采用偶极子声源在井中向外发射横波，并在井中接收由井外地质体反射回来的横波。2009年唐晓明利用四分量偶极数据对井旁裂缝和盐丘内部的精细构造进行了成像，并取得了显著的应用效果。目前偶极横波远探测技术的实施过程主要分为三个步骤：第一步是横波反射波提取，第二步是SH、SV波分离，第三步是偏移成像。其中，横波反射波提取是一个非常关键的处理环节，该步骤主要是做两方面的工作：一是尽可能压制井孔直达波信号，二是滤除掉无效反射波信号，比如地层界面反射波等。当地层弹性性质或井孔尺寸在轴向上变化较大时，便会产生幅度很强的地层界面反射波，地层界面反射波在波形中持续时间长，覆盖面积大，可能淹没有效横波反射波信号。因此，在获取有效横波反射波之前，需要剔除地层界面反射波。迄今，国内外学者针对横波反射波提取，主要是在压制井孔直达波方面做工作，比如利用F-K滤波来滤除井孔直达波信号，而对地层界面反射波消除方面的研究较少。为此，有必要研究一种去除地层界面反射波的方法，以实现对地层界面反射波的有效去除。
技术实现思路
本申请实施例的目的是提供一种去除地层界面反射波的方法及装置，以实现对地层界面反射波的有效去除。为解决上述技术问题，本申请实施例提供一种去除地层界面反射波的方法及装置是这样实现的：一种去除地层界面反射波的方法，提供有目的工区中指定钻井位置处的偶极横波数据；所述方法包括：确定所述偶极横波数据对应的偶极横波时差；对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据；基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据。优选方案中，在确定所述偶极横波数据对应的偶极横波时差之前，所述方法还包括：对所述偶极横波数据进行预处理，得到预处理后的偶极横波数据；相应的，所述确定所述偶极横波数据对应的偶极横波时差，包括：确定所述预处理后的偶极横波数据对应的偶极横波时差；所述对所述偶极横波数据进行垂直中值滤波，包括：对所述预处理后的偶极横波数据进行垂直中值滤波。优选方案中，所述对所述偶极横波数据进行预处理，包括：对所述偶极横波收进行增益恢复处理，得到增益恢复后的偶极横波数据；对所述增益恢复后的偶极横波数据进行延迟恢复处理，得到所述预处理后的偶极横波数据。优选方案中，所述对所述预处理后的偶极横波数据进行垂直中值滤波，包括：沿与水平面垂直的同相轴方向对所述预处理后的偶极横波数据进行中值滤波。优选方案中，所述根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据，包括：将所述预处理后的偶极横波数据减去所述初始偶极横波直达波数据，得到所述指定钻井位置处的初始偶极横波反射波数据。优选方案中，所述确定所述预处理后的偶极横波数据对应的偶极横波时差，包括：采用时间-慢度相关方法确定所述预处理后的偶极横波数据对应的偶极横波时差；具体地，在时间和慢度两个维度上计算所述预处理后的偶极横波数据对应的相关函数；将所述相关函数的最大函数值对应的慢度作为所述偶极横波时差。优选方案中，采用下述公式表征所述相关函数：其中，Corr(DT,Time)表示所述相关函数的函数值；Xm[t+DT(m-1)·RRSP]表示所述预处理后的偶极横波数据中第m道声波波形数据中时间为t+DT(m-1)·RRSP处的振幅值，t表示时间，DT表示慢度，RRSP表示偶极横波测井仪器中两个相邻接收器之间的距离，N表示所述预处理后的偶极横波数据中声波波形数据的道数；Time表示预设对比时窗的起始时间，CWLENGTH表示所述预设对比时窗的窗长；符号|·|为取绝对值运算符。优选方案中，所述地层界面反射波数据包括：地层界面上行反射波数据和地层界面下行反射波数据；所述基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据，包括：基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜上行中值滤波，得到倾斜上行中值滤波后的初始偶极横波反射波数据；将所述初始偶极横波反射波数据减去所述倾斜上行中值滤波后的初始偶极横波反射波数据，得到第一偶极横波反射波数据；所述第一偶极横波反射波数据表示去除所述地层界面上行反射波数据后的初始偶极横波反射波数据；基于所述偶极横波时差，对所述第一偶极横波反射波数据进行倾斜下行中值滤波，得到倾斜下行中值滤波后的第一偶极横波反射波数据；将所述第一偶极横波反射波数据减去所述倾斜下行中值滤波后的第一偶极横波反射波数据，得到第二偶极横波反射波数据；所述第二偶极横波反射波数据表示去除所述地层界面上行反射波数据和所述地层界面下行反射波数据后的初始偶极横波反射波数据。优选方案中，所述基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜上行中值滤波，得到倾斜上行中值滤波后的偶极横波反射波数据，包括：将斜率为所述偶极横波时差的倒数的两倍的同相轴作为倾斜上行同相轴；沿所述倾斜上行同相轴对所述初始偶极横波反射波数据进行倾斜上行中值滤波，得到所述倾斜上行中值滤波后的初始偶极横波反射波数据。优选方案中，所述基于所述偶极横波时差，对所述第一偶极横波反射波数据进行倾斜下行中值滤波，得到倾斜下行中值滤波后的第一偶极横波反射波数据，包括：将斜率为所述偶极横波时差的倒数的两倍的负数的同相轴作为倾斜下行同相轴；沿所述倾斜下行同相轴对所述第一偶极横波反射波数据进行倾斜下行中值滤波，得到所述倾斜下行中值滤波后的第一偶极横波反射波数据。一种去除地层界面反射波的装置，所述装置提供目的工区中指定钻井位置处的偶极横波数据；所述装置包括：时差确定模块、垂直中值滤波模块和倾斜中值滤波模块；其中，所述时差确定模块，用于确定所述偶极横波数据对应的偶极横波时差；所述垂直中值滤波模块，用于对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据；所述倾斜中值滤波模块，用于基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据。本申请实施例提供了一种去除地层界面反射波的方法及装置，可以确定所述偶极横波数据对应的偶极横波时差；可以对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据；可以基于所述偶极横波时差，对所本文档来自技高网...

【技术保护点】
一种去除地层界面反射波的方法，其特征在于，提供有目的工区中指定钻井位置处的偶极横波数据；所述方法包括：确定所述偶极横波数据对应的偶极横波时差；对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据；基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据。

【技术特征摘要】
1.一种去除地层界面反射波的方法，其特征在于，提供有目的工区中指定钻井位置处的偶极横波数据；所述方法包括：确定所述偶极横波数据对应的偶极横波时差；对所述偶极横波数据进行垂直中值滤波，得到所述指定钻井位置处的偶极横波直达波数据；根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据；基于所述偶极横波时差，对所述初始偶极横波反射波数据进行倾斜中值滤波，以去除所述指定钻井位置处的地层界面反射波数据。2.根据权利要求1所述的一种去除地层界面反射波的方法，其特征在于，在确定所述偶极横波数据对应的偶极横波时差之前，所述方法还包括：对所述偶极横波数据进行预处理，得到预处理后的偶极横波数据；相应的，所述确定所述偶极横波数据对应的偶极横波时差，包括：确定所述预处理后的偶极横波数据对应的偶极横波时差；所述对所述偶极横波数据进行垂直中值滤波，包括：对所述预处理后的偶极横波数据进行垂直中值滤波。3.根据权利要求2所述的一种去除地层界面反射波的方法，其特征在于，所述对所述偶极横波数据进行预处理，包括：对所述偶极横波收进行增益恢复处理，得到增益恢复后的偶极横波数据；对所述增益恢复后的偶极横波数据进行延迟恢复处理，得到所述预处理后的偶极横波数据。4.根据权利要求2所述的一种去除地层界面反射波的方法，其特征在于，所述对所述预处理后的偶极横波数据进行垂直中值滤波，包括：沿与水平面垂直的同相轴方向对所述预处理后的偶极横波数据进行中值滤波。5.根据权利要求4所述的一种去除地层界面反射波的方法，其特征在于，所述根据所述偶极横波直达波数据确定所述指定钻井位置处的初始偶极横波反射波数据，包括：将所述预处理后的偶极横波数据减去所述初始偶极横波直达波数据，得到所述指定钻井位置处的初始偶极横波反射波数据。6.根据权利要求2所述的一种去除地层界面反射波的方法，其特征在于，所述确定所述预处理后的偶极横波数据对应的偶极横波时差，包括：采用时间-慢度相关方法确定所述预处理后的偶极横波数据对应的偶极横波时差；具体地，在时间和慢度两个维度上计算所述预处理后的偶极横波数据对应的相关函数；将所述相关函数的最大函数值对应的慢度作为所述偶极横波时差。7.根据权利要求6所述的一种去除地层界面反射波的方法，其特征在于，采用下述公式表征所述相关函数：其中，Corr(DT,Time)表示所述相关函数的函数值；Xm[t+DT(m-1)·RRSP]表示所述预处理后的偶极横波数据中第m道声波波形数据中...

【专利技术属性】
技术研发人员：刘鹏，李宁，武宏亮，冯庆付，王克文，冯周，李雨生，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京,11