CN101382598A

CN101382598A - 一种真三维地震数据线性噪音的压制方法

Info

Publication number: CN101382598A
Application number: CNA2007101212623A
Authority: CN
Inventors: 王梅生; 王彦娟; 李培明; 李伟波
Original assignee: China National Petroleum Corp; BGP Inc
Current assignee: China National Petroleum Corp; BGP Inc
Priority date: 2007-09-03
Filing date: 2007-09-03
Publication date: 2009-03-11
Anticipated expiration: 2027-09-03
Also published as: CN101382598B

Abstract

本发明是石油地球物理勘探三维地震数据处理三维地震数据线性噪音的压制方法。具体步骤：采集并记录高密度地震数据；将数据分成最小数据集数据组；对每一个数据组的所有道按炮检距大小和方位角进行排序，形成三维去噪的最小数据集；根据所形成的频率波数谱(FK谱)或原始地震数据分析线性噪音的速度，去除线性噪音；对还有线性噪音存在反复压制，直到噪音去除。本发明适合不同观测系统所采集的地震资料，特别是非十字交叉排列和非规则束状正交观测系统，对原始资料的线性噪音可以很好地压制。

Description

一种真三维地震数据线性噪音的压制方法

技术领域

本发明属于石油地球物理勘探技术，具体是一种三维地震数据处理中的真三维地震数据线性噪音的压制方法。

背景技术

在石油地震勘探的地震资料采集和处理过程中会遇到一些规则干扰波，如面波或浅层多次折射波，其能量强，对提高地震资料的信噪比带来很大困难，使处理后的资料不真实准确。目前，通常采用两步法噪音压制方法，即在X方向与Y方向分别进行压制。这样做不仅效率比较低，而且由于线性噪音在三维空间内是线性的，在非纵二维线上都不是线性的，因此利用两个方向进行二维去噪其压制效果也不是十分理想。

目前虽有适合于压制三维线性噪音压制(FKxKy)方法，但适合于十字交叉排列或能形成十字交叉排列正交束状宽方位观测系统。对于不规则观测系统，对于非十字交叉排列和非规则束状正交观测系统所采集资料中的线性噪音，难以进行三维线性噪音压制。

发明内容

本发明目的是提供一种适合不同观测系统所采集的地震资料，特别是非十字交叉排列和非规则束状正交观测系统所采集资料的真三维地震数据线性噪音的压制方法。

本发明通过以下具体步骤实现：

1)采用常规地震勘探的方法采集并记录高密度地震数据；

2)将采集的高密度空间数据分成最小数据集数据组；

所述的最小数据集数据组形成方法是将一条接收线与其对应的炮点抽取成一组形成最小数据集数据组。

所述的最小数据集数据组形成方法是利用三维一个模版中的所有接收线与纵向上以炮线距为单位、横向上以接收线距为单位，此范围内的所有炮资料形成一个最小数据集数据组。

3)对每一个数据组的所有道按炮检距大小和方位角进行排序，形成三维去噪的最小数据集；

4)根据所形成的频率波数谱(FK谱)或原始地震数据分析线性噪音的速度，去除线性噪音；

所述的去除线性噪音采用以下公式计算：

F (w, k_{x}, k_{y}) = \{\begin{matrix} 1 & w / {({k_{x}}^{2} + {k_{y}}^{2})}^{\frac{1}{2}} > V_{L} \\ 0 & w / {({k_{x}}^{2} + {k_{y}}^{2})}^{\frac{1}{2}} < V_{L} \end{matrix}

式中F是波场在频率波数域的对应值；w是频率；K_x和K_y是在x方向和y方向的波数；V_L是线性噪音的视速度。w、K_x和K_y是时间域数据经过傅立叶变换后频率波数域参数；V_L是滤波中所用的视速度的门槛值；

5)地震数据上还有线性噪音存在或包含有效信号，再次确定线性噪音的速度并采用步骤4)进行压制，直到噪音去除。

本发明适合不同观测系统所采集的地震资料，特别是非十字交叉排列和非规则束状正交观测系统。本发明在三维空间上线性干扰波为一锥体，锥体在时间切片上为一个圆，通过对构建的最小数据集分析，见附图1。

在进行三维线性压制后，原始资料的线性噪音得到了很好地压制，可以从图2中看出：无论是Inline线，还是Crossline线，原始资料的线性噪音都得到了很好的压制。

附图说明

图1a为构建与经过炮检距与方位角调整的最小数据集对比图；

图1b直接构建数据集的时间切片图；

图2按炮检距方位角调整后的时间切片图；

图3利用原始资料线性噪音的速度图；

图4线性噪音压制不足的噪音压制前、噪音压制后、去除噪音效果图；

图5线性噪音压制过量的噪音压制前、噪音压制后、去除噪音效果图；

图6本发明线性噪音压制良好的效果图。

具体实施方式

本发明通过以下具体题步骤实现：

1)采用常规地震勘探的方法采集并记录高密度地震数据；

2)将采集的高密度空间采样数据分成能够形成三维线性噪音压制最小数据集数据组；

所述的最小数据集数据组形成方法有两种：一是将一条接收线与其对应的炮点抽取成一组形成最小数据集数据组；或利用三维一个模版中的所有接收线与纵向上以炮线距为单位、横向上以接收线距为单位，此范围内的所有炮资料形成一个最小数据集数据组。

图1a为构建数据集对应的炮检点，图1b为检查所构成的数据集的时间切片，从时间切片上可以看出其不是规则的圆，说明该数据集存在一定的问题。

对上面的数据集进行按炮检距与方位角重排，对重排后的数据集进行时间切片，时间切片变为规则的圆，见图2。说明其变为了适合三维去噪的最小数据集。

4)进行三维线性噪音压制。根据所形成的FK谱或原始地震数据分析线性噪音的速度，见图3。选择最大线性噪音速度作为V_L，进行噪音压制。

去除线性噪音的表达式为：

F (w, k_{x}, k_{y}) = \{\begin{matrix} 1 & w / {({k_{x}}^{2} + {k_{y}}^{2})}^{\frac{1}{2}} > V_{L} \\ 0 & w / {({k_{x}}^{2} + {k_{y}}^{2})}^{\frac{1}{2}} < V_{L} \end{matrix}

式中F是波场在频率波数域的对应值；w是频率；K_x和K_y是在x方向和y方向的波数；V_L是线性噪音的视速度。w、K_x和K_y是时间域数据经过傅立叶变换后频率波数域参数；V_L是滤波中所用的视速度的门槛值。

5)对去噪效果进行分析，可以从两方面进行分析：一是检查线性噪音去除后的地震数据上是否还有线性噪音存在，若仍有线性噪音存在，说明压制不足；二是检查去除的噪音是否包含有效信号，如果包含有效信号说明噪音去除过量如图5，上述两种情况均需要再次分析线性噪音的速度进行压制。如果线性噪音压制后的地震数据没有线性噪音且去除的噪音也不包含有效信号，说明噪音去除效果良好见图6。

Claims

1、一种真三维地震数据线性噪音的压制方法，其特征在于过以下具体步骤实现：

1)采用常规地震勘探的方法采集并记录高密度地震数据；

2)将采集的高密度空间数据分成最小数据集数据组；

2、根据原理要求1所述的真三维地震数据线性噪音的压制方法，其特征在于：步骤2)所述的最小数据集数据组形成方法是将一条接收线与其对应的炮点抽取成一组形成最小数据集数据组。

3、根据原理要求1所述的真三维地震数据线性噪音的压制方法，其特征在于：步骤2)所述的最小数据集数据组形成方法是利用三维一个模版中的所有接收线与纵向上以炮线距为单位、横向上以接收线距为单位，此范围内的所有炮资料形成一个最小数据集数据组。

4、根据原理要求1所述的真三维地震数据线性噪音的压制方法，其特征在于：步骤4)所述的去除线性噪音采用以下公式计算：

F (w, k_{x}, k_{y}) = \{\begin{matrix} 1 & w / {({k_{x}}^{2} + {k_{y}}^{2})}^{\frac{1}{2}} > V_{L} \\ 0 & w / {({k_{x}}^{2} + {k_{y}}^{2})}^{\frac{1}{2}} < V_{L} \end{matrix}