[发明专利]一种噪声压制方法、系统、存储介质及电子设备在审
申请号: | 202210651638.6 | 申请日: | 2022-06-09 |
公开(公告)号: | CN115016004A | 公开(公告)日: | 2022-09-06 |
发明(设计)人: | 王明星;张庆淮;王国来;路交通;张强;赵配叶;吴安楚;李华喜;唐旷怡 | 申请(专利权)人: | 中石化石油工程技术服务有限公司;中石化石油工程地球物理有限公司 |
主分类号: | G01V1/36 | 分类号: | G01V1/36;G01V1/30;G01V1/38 |
代理公司: | 北京聿宏知识产权代理有限公司 11372 | 代理人: | 陈超德;吴昊 |
地址: | 100029 北京*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 噪声 压制 方法 系统 存储 介质 电子设备 | ||
1.一种噪声压制方法,其特征在于,所述方法包括:
对地震数据进行预处理;
基于所述预处理后的地震数据,确定所述地震数据中的舒尔特波的加权系数;
基于所述舒尔特波的加权系数,确定所述地震数据的第一偏振特征信息;
基于所述地震数据的第一偏振特征信息,压制所述地震数据中的舒尔特波;
基于压制舒尔特波后的地震数据,确定所述地震数据中的水层波的加权系数;
基于所述水层波的加权系数,确定所述地震数据的第二偏振特征信息;
基于所述地震数据的第二偏振特征信息,压制所述地震数据中的水层波,以得到压制噪声后的地震信号。
2.根据权利要求1所述的方法,其特征在于,所述地震数据包括水听器数据和地震检波器数据,所述基于所述预处理后的地震数据,确定舒尔特波加权系数,包括:
分别获取所述水听器数据和所述地震检波器中的舒尔特波振幅;
基于所述地震检波器中的舒尔特波振幅和所述水听器数据中的舒尔特波振幅,通过加权系数计算式确定所述舒尔特波加权系数;
所述加权系数计算式为:
ws=Gscholte/Hscholte;
其中,Gscholte为地震检波器中舒尔特波振幅,Hscholte为水听器中舒尔特波振幅,ws为舒尔特波加权系数。
3.根据权利要求2所述的方法,其特征在于,所述偏振特征信息包括:椭圆率,所述基于所述舒尔特波的加权系数,确定所述地震数据的第一偏振特征信息,包括:
基于所述舒尔特波的加权系数,分别确定所述水听器数据在预设时窗内的加权平均值和所述地震检波器数据在预设时窗内的的平均值;
基于所述水听器数据的加权平均值和所述地震检波器数据的平均值,确定所述地震数据的特征值;
基于所述特征值,通过椭圆率计算式确定所述地震数据的第一椭圆率;
所述椭圆率计算式为:
其中,εs为所述地震数据的椭圆率,λ1和λ2分别为所述地震数据的特征值。
4.根据权利要求3所述的方法,其特征在于,所述基于所述舒尔特波的加权系数,分别确定所述水听器数据在预设时窗内的加权平均值和所述地震检波器数据在预设时窗内的的平均值,包括:
根据加权平均值计算式,确定所述水听器数据中地震数据在预设时窗内的加权平均值
所述加权平均值计算式为:
根据平均值计算式,确定所述地震检波器数据中地震数据在预设时窗内的平均值
所述平均值计算式为:
其中,Hi表示第i时刻所述水听器获得的地震数据的振幅值,Gi表示第i时刻所述地震检波器获得的地震数据的振幅值,t1~t2为预设时窗。
5.根据权利要求4所述的方法,其特征在于,所述基于所述水听器数据的加权平均值和所述地震检波器数据的平均值,确定所述地震数据的特征值,包括:
基于所述水听器数据、所述地震检波器数据、所述水听器数据的加权平均值和所述地震检波器数据的平均值,构建所述地震数据中所述水听器数据和所述地震检波器数据的协方差函数;
所述协方差函数为:
基于所述协方差函数,确定协方差矩阵为:
基于所述协方差矩阵,确定特征值条件关系式:Cu=λu,其中λ为特征值,u为特征向量;
在λ满足所述特征值条件关系式的情况下,确定所述地震数据的特征值λ1和λ2。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中石化石油工程技术服务有限公司;中石化石油工程地球物理有限公司,未经中石化石油工程技术服务有限公司;中石化石油工程地球物理有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202210651638.6/1.html,转载请声明来源钻瓜专利网。