一种压制同步激发采集邻炮干扰的方法技术

技术编号:11701726 阅读:143 留言:0更新日期:2015-07-09 01:15
一种压制同步激发采集邻炮干扰的方法。本发明专利技术是地球物探勘探数据处理中的噪声压制技术,将地震数据按照共炮检距道集记录方式排列,计算地震数据时间空间窗中的中值矢量,将时空窗中每一个地震数据道到其他地震数据道的距离和从小到大进行排序,然后舍弃掉L个距离和最大的地震道,计算时空窗中剩下的地震数据道的中心样点振幅值的平均值,用平均值的整数倍作为阈值对时空窗的中心样点进行滤波,所有样点都处理后完成压制同步激发采集邻炮干扰。本发明专利技术相对于矢量中值滤波,在压制邻炮干扰的同时,能更好地保留主炮能量,具有更高的计算效率,满足海量数据的需求,有很好的鲁棒性,适用于各类数据。

【技术实现步骤摘要】
一种压制同步激发采集邻炮干扰的方法
本专利技术属于地球物探勘探的地震数据处理技术,进一步细分可以作为地震数据处理中的噪声压制技术。
技术介绍
高密度地震勘探是物探采集技术的发展方向之一,使用传统的采集方式进行高密度勘探需要很高的勘探成本,而高效采集技术可以大幅提高生产效率,降低勘探成本。同步激发采集是一种非常有效的高效采集技术,但是其数据存在严重的邻炮能量干扰,有效压制这些邻炮能量可以得到更好的处理效果。目前邻炮压制技术分为两大类,第一类是基于迭代反演算法,第二类是基于压噪类算法;基于迭代反演算法需要同步激发采集时记录每炮的延时数据,然后用此延时数据构造混叠参数矩阵,使用该矩阵可以把邻炮数据分离抽象成一个矩阵求解问题。由于混叠矩阵的逆往往不存在,所以不能通过矩阵求逆直接得到分离后的数据,而是通过混叠参数矩阵的共轭转置与混叠数据矩阵的相乘结果作为反演初值进行迭代反演,这种迭代反演算法效率很低,并且这种算法对数据信噪比要求较高,所以很难适用于陆上低信噪比的海量数据。根据同步激发采集方式的特点,相同的排列接收了不同震源在不同的位置以不同的时间起震所激发的反射波、直达波等能量。在共激发点道集本文档来自技高网...
一种压制同步激发采集邻炮干扰的方法

【技术保护点】
一种压制同步激发采集邻炮干扰的方法,特点是采用以下具体实施步骤:1)使用同步激发采集方式采集地震数据,按照共炮检距道集记录方式排列;2)计算地震数据时间空间窗中的中值矢量;所述的中值矢量计算公式为:Σi=1N||Amv-Ai||1≤Σi=1N||Aj-Ai||1]]>式中:N是时间空间窗中空间方向的道数,Aj和Ai是时间空间窗中的任意两个不同的地震数据道,时空窗中每个数据道在时间方向有M个样点,i和j的取值范围是1到N,表示时空窗中任意一个地震数据道到其他地震数据道的距离的和,Amv是时空窗中的中值矢量;所述的任意两个不同地震数据道的距离是通过两个地震数据道中所有...

【技术特征摘要】
1.一种压制同步激发采集邻炮干扰的方法,特点是采用以下具体实施步骤:1)使用同步激发采集方式采集地震数据,按照共炮检距道集记录方式排列;2)计算地震数据时间空间窗中的中值矢量;所述的中值矢量计算公式为:式中:N是时间空间窗中空间方向的道数,Aj和Ai是时间空间窗中的任意两个不同的地震数据道,时空窗中每个数据道在时间方向有M个样点,i和j的取值范围是1到N,表示时空窗中任意一个地震数据道到其他地震数据道的距离的和,Amv是时空窗中的中值矢量;所述的任意两个不同地震数据道的距离是通过两个地震数据道中所有对应序号样点的差的绝对值的和得到的,计算公式为:式中:ajk表示数据道Aj的第k个样点,aik表示数据道Ai的第k个样点,k的取值范围是1到M;3)把步骤2)时间空间窗中每一个地震数据道到其他地震数据道的距离和从小到大进行排序,然后舍弃掉L个距离和最大的地震道;所述L个距离的值通过下式计算:L=ceil(alpha×N)式中:alpha的取值范围为0到0.5,N时间空间窗中空间方向的道数,ceil()为对括号中的数取大于等于它的最小整数;4)计算时空窗中剩下的地震数据道的中心样点振幅值的平均值,计算公式为:式中:N是时间空间窗中空间方向的道数,L为步骤3)中得到的整数,x的取值范围是1...

【专利技术属性】
技术研发人员:王文闯李合群赵波罗国安钱忠平张旭东
申请(专利权)人:中国石油集团东方地球物理勘探有限责任公司中油油气勘探软件国家工程研究中心有限公司
类型:发明
国别省市:河北;13

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1