一种从地震记录中去除面波的方法和装置制造方法及图纸

本申请实施例公开了一种从地震记录中去除面波的方法和装置。该方法包括：对待处理的地震记录进行分频扫描，获得面波频段的分布范围；根据面波频段的分布范围，对地震记录进行小波分频，获取不同频段上的地震记录；根据面波频段的分布范围，从上述不同频段上的地震记录中识别出包含面波的地震记录，对包含面波的地震记录进行径向道变换；根据面波频段的分布范围确定高通滤波的频率，对径向道变换后的地震记录在径向道域进行高通滤波；对高通滤波后的地震记录进行反径向道变换；将反径向道变换后的地震记录和不同频段上的地震记录中不包含面波的地震记录进行小波重构，得到去除面波的地震记录。

【技术实现步骤摘要】
-种从地震记录中去除面波的方法和装置
本申请涉及石油勘探地震资料处理
，尤其涉及一种从地震记录中去除面 波的方法和装置。
技术介绍
由于面波存在广泛，能量强，严重影响了中深层的有效反射，降低了地震资料的信 噪比，所W面波消除成为地震资料处理中的一项重要工作，针对面波的特点，可W采取多种 去除面波的方法，如径向道变换、小波变换等方法，但是该些单个消除面波的方法往往达不 到理想的效果。 小波变换是将地震记录在多个不同的尺度上进行小波分解，各个尺度上经过小波 分解所得到的小波变换系数就代表原地震记录在不同的时间或频率上的信息，该样就可W 利用面波与有效波在频率上的差异，将面波与有效波区分开来，进而去除面波，但是利用小 波变换去除的面波有残留，去噪能力弱。 径向道变换是一种简单的几何地震道映射算法，由于面波与地震记录在速度上存 在巨大的差异，该方法通过坐标系的变换，将面波与有效的地震记录区分开来，进而去除面 波，但会有部分有效的地震记录也被去除。
技术实现思路
本申请的目的是提供一种从地震记录中去除面波的方法和装置，为了在减少地震 记录损失的情况下，有效的去除面波。 为了实现上述目的，本申请提供了一种从地震记录中去除面波的方法，该方法包 括： S101 ;对待处理的地震记录进行分频扫描，获得面波频段的分布范围； S102;根据所述面波频段的分布范围，对所述地震记录进行小波分频，获取不同频 段上的地震记录，所述不同频段上的地震记录包括包含面波的地震记录和不包含面波段的 地震记录； S103;根据所述面波频段的分布范围，从所述不同频段上的地震记录中识别出包 含面波的地震记录，对所述包含面波的地震记录进行径向道变换； S104 ;根据所述面波频段的分布范围确定高通滤波的频率，对S103中径向道变换 后的地震记录在径向道域进行高通滤波； S105 ;对S104中高通滤波后的地震记录进行反径向道变换； S106 ;将S105中反径向道变换后的地震记录和S102中所述的不包含面波的地震 记录进行小波重构，得到去除面波的地震记录。 在一个优选的实施例中，所述地震记录记录了多道采样点，每道有多个采样点，所 述采样点记录了信号的振幅值。 在一个优选的实施例中，所述对所述地震记录进行分频扫描，获得面波频段的分 布范围包括： 对所述地震记录进行频谱分析，获得地震记录的频段范围； 将所述地震记录进行频段划分，获得不同的频段的地震记录； 对所述地震记录在每个频段范围内进行频率滤波，获得不同频段的时间域的地震 记录； 根据所述不同频段的时间域的地震记录确定面波，获得面波频段的分布范围。 在一个优选的实施例中，所述根据所述面波频段的分布范围，对所述地震记录进 行小波分频，获取不同频段上的地震记录包括： 根据面波频段的分布范围设置所述小波分频的分频范围； 采用马拉特（Mallat)算法对所述地震记录进行分频，把地震记录分解成地震记 录小波细节和地震记录小波逼近，并将上一次分解的地震记录小波逼近作为下一次输入， 获取不同频段上的地震记录，公式如下： [002引W，，/ = 5^/. *// 5] |/ = 5'] / *左 上式中，Hi曲，Low为一对正交镜像滤波器；j是分解次数，j= 0, 1，......j-1,j =0表示第1次分解皆是分析尺度，是第j+2次地震记录小波细节，*V/是第j+1 次地震记录小波逼近，|/是第j+2次地震记录小波逼近。 [00巧]在一个优选的实施例中，所述对所述包含面波的地震记录进行径向道变换包括： 将t-x域的所述包含面波的地震记录变换到X-V域，并利用H次样条插值来实现 径向道变换。从t-x域的变换到T-V域的映射关系如下： t一X = t-t〇 X-V= (x-x〇) / (t-t〇) 上式中，t为双程旅行时；X为偏移距；X。为径向道变换的原点位置；t。为径向道变 换的原点位置的双程旅行时；T为延迟时；V为速度。 在一个优选的实施例中，所述对S104中高通滤波后的地震记录进行反径向道变 换包括： 将T-V域的所述高通滤波后的地震记录变换到t-x域，从t-x域的变换到X-V 域的映射关系如下： X一t = X +t〇 V一X=V +x〇 上式中，t为双程旅行时；X为偏移距；X。为径向道变换的原点位置；t。为径向道变 换的原点位置的双程旅行时；T为延迟时；V为速度。 [00巧]本申请另一方面还提供一种从地震记录中去除面波的装置，该装置包括： 面波频段确定单元，用于对待处理的地震记录进行分频扫描，获得面波频段的分 布范围； 分频单元，用于根据所述面波频段的分布范围，对所述地震记录进行小波分频，获 取不同频段上的地震记录，所述不同频段上的地震记录包括包含面波的地震记录和不包含 面波段的地震记录； 径向道变换单元，用于根据所述面波频段的分布范围，从所述不同频段上的地震 记录中识别出包含面波的地震记录，对所述包含面波的地震记录进行径向道变换； 滤波单元，用于根据所述面波频段的分布范围确定高通滤波的频率，对径向道变 换单元得到的径向道变换后的地震记录在径向道域进行高通滤波； 反径向道变换单元，用于对滤波单元得到的高通滤波后的地震记录进行反径向道 变换； 小波重构单元，用于将反径向道变换单元中得到的反径向道变换后的地震记录和 分频单元得到的所述的不包含面波的地震记录进行小波重构，得到去除面波的地震记录。 由W上本申请实施例提供的技术方案可见，本申请通过对待处理的地震记录进行 分频扫描，获得面波频段的分布范围，根据面波频段的分布范围，对待处理地震记录进行小 波分频，获得不同频段的地震记录，从不同频段的地震记录识别出包含面波的地震记录，进 行径向道变换，利用面波与有效地震记录在速度上的差异，在径向道域进行高频滤波，再将 滤波后的地震记录进行径向道反变换，将径向道反变换后的地震记录与上述小波分频后得 到的不同频段的地震记录中不包含面波的地震记录进行小波重构，获得去除面波的地震记 录。与现有技术相比，不仅去除了面波，且大大减少了有效地震记录的损失。 【附图说明】 为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提 下，还可W根据该些附图获得其他的附图。 图1是本申请实施例提供的一种从地震记录中去除面波的方法的流程图； 图2是待处理的地震记录剖面图； 图3是图2中的地震记录通过小波分频去除面波的方法得到的地震记录剖面图； 图4是图2中的地震记录通过径向道变换去除面波的方法得到的地震记录剖面 图； 图5是图2中的地震记录通过本申请实施例提供一种从地震记录中去除面波的方 法得到的地震记录剖面图； 图6是本申请实施例提供的一种从地震记录中去除面波的装置的示意图。 【具体实施方式】 为了使本
的人员更好地理解本申请中的本文档来自技高网...

【技术保护点】
一种从地震记录中去除面波的方法，其特征在于，该方法包括：S1：对待处理的地震记录进行分频扫描，获得面波频段的分布范围；S2：根据所述面波频段的分布范围，对所述地震记录进行小波分频，获取不同频段上的地震记录，所述不同频段上的地震记录包括包含面波的地震记录和不包含面波段的地震记录；S3：根据所述面波频段的分布范围，从所述不同频段上的地震记录中识别出包含面波的地震记录，对所述包含面波的地震记录进行径向道变换；S4：根据所述面波频段的分布范围确定高通滤波的频率，对S3中径向道变换后的地震记录在径向道域进行高通滤波；S5：对S4中高通滤波后的地震记录进行反径向道变换；S6：将S5中反径向道变换后的地震记录和S2中所述的不包含面波的地震记录进行小波重构，得到去除面波的地震记录。

【技术特征摘要】
1. 一种从地震记录中去除面波的方法，其特征在于，该方法包括： 51 :对待处理的地震记录进行分频扫描，获得面波频段的分布范围； 52 :根据所述面波频段的分布范围，对所述地震记录进行小波分频，获取不同频段上的 地震记录，所述不同频段上的地震记录包括包含面波的地震记录和不包含面波段的地震记 录； 53 :根据所述面波频段的分布范围，从所述不同频段上的地震记录中识别出包含面波 的地震记录，对所述包含面波的地震记录进行径向道变换； 54 :根据所述面波频段的分布范围确定高通滤波的频率，对S3中径向道变换后的地震 记录在径向道域进行高通滤波； 55 :对S4中高通滤波后的地震记录进行反径向道变换； 56 :将S5中反径向道变换后的地震记录和S2中所述的不包含面波的地震记录进行小 波重构，得到去除面波的地震记录。2. 根据权利要求1所述的从地震记录中去除面波的方法，其特征在于，所述地震记录 记录了多道采样点，每道有多个采样点，所述采样点记录了信号的振幅值。3. 根据权利要求1所述的从地震记录中去除面波的方法，其特征在于，所述对所述地 震记录进行分频扫描，获得面波频段的分布范围包括： 对所述地震记录进行频谱分析，获得地震记录的频段范围； 将所述地震记录进行频段划分，获得不同的频段的地震记录； 对所述地震记录在每个频段范围内进行频率滤波，获得不同频段的时间域的地震记 录； 根据所述不同频段的时间域的地震记录确定面波，获得面波频段的分布范围。4. 根据权利要求1所述的从地震记录中去除面波的方法，其特征在于，所述根据所述 面波频段的分布范围，对所述地震记录进行小波分频，获取不同频段上的地震记录包括： 根据面波频段的分布范围设置所述小波分频的分频范围； 采用马拉特（Mallat)算法对所述地震记录进行分频，把地震记录分解成地震记录小 波细节和地震记录小波逼近，并将上一次分解的地震记录小波逼近作为下一次输入，获取 不同频段上的地震记录，公式如下： W2 ,/ = 52 J S#f = S:.PL 上式中，High, Low为一对正交镜像滤波器；j是分解次数，j = 0,1，......j-1, j = 〇 表示第1次分解；2...

【专利技术属性】
技术研发人员：王万里，杨午阳，高建虎，魏新建，何欣，禄娟，陈德武，杨庆，鄢高韩，李琳，李三家，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京;11