可控震源独立同步激发地震数据邻炮干扰压制方法和装置制造方法及图纸

本发明专利技术提供了一种可控震源独立同步激发地震数据邻炮干扰压制方法和装置，该方法包括：在可控震源ISS地震数据中分选出非共炮点道集地震数据；对分选出的非共炮点道集地震数据进行矢量中值滤波，得到邻炮干扰压制后的非共炮点道集地震数据；从邻炮干扰压制后的非共炮点道集地震数据中分选出共炮点道集地震数据。本发明专利技术解决了现有技术中解决了现有技术中获取的地震数据的准确性不高的技术问题，达到了有效消除ISS地震数据中的邻炮干扰，提高地震数据准确性的技术效果。

【技术实现步骤摘要】
可控震源独立同步激发地震数据邻炮干扰压制方法和装置
本专利技术涉及地球物理勘探
，特别涉及一种可控震源独立同步激发（IndependentSimultaneousSweeping，ISS）地震数据邻炮干扰压制方法和装置。
技术介绍
可控震源ISS地震数据在采集的过程中是通过多组可控震源独立作业、同步激发，并通过连续观测的方式将多组可控震源同步激发产生的地震波场记录到一个母记录中的方式获得的。这种获取可控震源ISS地震数据的方法极大地压缩了地震数据的数据量、提高了获取地震数据的效率。然而，可控震源ISS地震数据在采集中可控震源同步激发产生的地震波场之间相互会产生干扰、混叠，从而导致地震数据的准确性不高。
技术实现思路
本专利技术提供了一种可控震源ISS地震数据邻炮干扰压制方法和装置，以达到提高地震数据准确性的目的。本专利技术实施例提供了一种可控震源ISS地震数据邻炮干扰压制方法，包括：在可控震源ISS地震数据中分选出非共炮点道集地震数据；对分选出的非共炮点道集地震数据进行矢量中值滤波，得到邻炮干扰压制后的非共炮点道集地震数据；从邻炮干扰压制后的非共炮点道集地震数据中分选出共炮点道集地震数据。在一个实施例中，对分选出的非共炮点道集地震数据进行矢量中值滤波，包括：在时域对分选出的非共炮点道集地震数据进行矢量中值滤波；或者，在频域对分选出的非共炮点道集地震数据进行矢量中值滤波；或者，在时域与频域对分选出的非共炮点道集地震数据进行矢量中值滤波。在一个实施例中，在时域对分选出的非共炮点道集地震数据进行矢量中值滤波，包括：沿时间轴按预定长度将非共炮点道集地震数据分成段；将各段非共炮点道集地震数据作为矢量，按以下公式按地震数据道空间分布对各段非共炮点道集地震数据进行矢量中值滤波，得到邻炮干扰压制后的非共炮点道集地震数据：Σi=1n||Xm,k-Xi,k||p≤Σi=1n||Xj,k-Xi,k||p]]>其中，Xm,k∈{Xi,k|i=1,...,n}，j=1,...,n，Xm,k为输出的结果，k为对非共炮点道集地震数据分段后各段的顺序号，Xi,k为分段后各段的非共炮点道集地震数据，m、i和j为地震数据道的顺序号，n为进行矢量中值滤波的道数，为矢量Lp范数。在一个实施例中，所述预定长度包括：一个子波或者一个波形对应时间轴上的长度。在一个实施例中，在频域对分选出的非共炮点道集地震数据进行矢量中值滤波，包括：按照以下公式对分选出的非共炮点道集地震数据进行傅里叶变换，得到频域的非共炮点道集地震数据：X(ω)i=FFT[X(t)i]其中，X(t)i为时域的非共炮点道集地震数据，FFT[]为一维的傅里叶变换，ω为角频率，X(ω)i为傅里叶变换后得到的频域的非共炮点道集地震数据；按照以下公式对经过傅里叶变换后的非共炮点道集地震数据，按地震数据道空间分布进行矢量中值滤波：Σi=1n||X(ω)m-x(ω)i||p≤Σi=1n||X(ω)j-X(ω)i||p]]>其中，X(ω)m为对频域的非共炮点道集地震数据进行矢量中值滤波后输出的结果，ω为角频率，m、i和j为地震数据道的顺序号，n为进行矢量中值滤波运算的道数，为矢量Lp范数；将经过矢量中值滤波后的频域的非共炮点道集数据进行反傅里叶变换，得到邻炮干扰压制后的时域的非共炮点道集地震数据。在一个实施例中，非共炮点道集包括以下至少之一：共接收点道集、共偏移距道集、共中心点道集、十字排列道集。本专利技术实施例还提供了一种可控震源独立同步激发ISS地震数据邻炮干扰压制装置，包括：第一分选模块，用于在可控震源ISS地震数据中分选出非共炮点道集地震数据；滤波模块，用于对分选出的非共炮点道集地震数据进行矢量中值滤波，得到邻炮干扰压制后的非共炮点道集地震数据；第二分选模块，用于从邻炮干扰压制后的非共炮点道集地震数据中分选出共炮点道集地震数据。在一个实施例中，所述滤波模块包括：第一滤波单元，用于在时域对分选出的非共炮点道集地震数据进行矢量中值滤波；或者，第二滤波单元，用于在频域对分选出的非共炮点道集地震数据进行矢量中值滤波；或者，第三滤波单元，用于在时域与频域对分选出的非共炮点道集地震数据进行矢量中值滤波。在一个实施例中，所述第一滤波单元包括：分段子单元，用于沿时间轴按预定长度将非共炮点道集地震数据分成段；第一滤波子单元，用于将各段非共炮点道集地震数据作为矢量，按以下公式按地震数据道空间分布对各段非共炮点道集地震数据进行矢量中值滤波，得到邻炮干扰压制后的非共炮点道集地震数据：Σi=1n||Xm,k-Xi,k||p≤Σi=1n||Xj,k-Xi,k||p]]>其中，Xm,k∈{Xi,k|i=1,...,n}，j=1,...,n，Xm,k为输出的结果，k为对非共炮点道集地震数据分段后各段的顺序号，Xi,k为分段后各段的非共炮点道集地震数据，m、i和j为地震数据道的顺序号，n为进行矢量中值滤波的道数，为矢量Lp范数。在一个实施例中，所述预定长度包括：一个子波或者一个波形对应时间轴上的长度。在一个实施例中，所述第二滤波子单元包括：第一变换子单元，用于按照以下公式对分选出的非共炮点道集地震数据进行傅里叶变换，得到频域的非共炮点道集地震数据：X(ω)i=FFT[X(t)i]其中，X(t)i为时域的非共炮点道集地震数据，FFT[]为一维的傅里叶变换，ω为角频率，X(ω)i为傅里叶变换后得到的频域的非共炮点道集地震数据；第二滤波子单元，用于按照以下公式对经过傅里叶变换后的非共炮点道集地震数据，按地震数据道空间分布进行矢量中值滤波：Σi=1n||X(ω)m-x(ω)i||p≤Σi=1n||X(ω)j-X(ω)i||p]]>其中，X(ω)m为对频域的非共炮点道集地震数据进行矢量中值滤波后输出的结果，ω为角频率，m、i和j为地震数据道的顺序号，n为进行矢量中值滤波运算的道数，为矢量Lp范数；第二变换子单元，用于将经过矢量中值滤波后的频域的非共炮点道集数据进行反傅里叶变换得到邻炮干扰压制后的时域的非共炮点道集地震数据。在一个实施例中，非共炮点道集包括以下至少之一：共接收点道集、共偏移距道集、共中心点道集、十字排列道集。在本专利技术实施例中，在可控震源ISS地震数据中分选出非共炮点道集地震数据，然后对分选出的非共炮点道集地震数据进行矢量中值滤波，以压制邻炮干扰，然后再从进行邻炮干扰压制后的非共炮点道集地震数据中分选出共炮点道集地震数据，以得到消除了邻炮干扰的共炮点道集地震数据。通过上述方式解决了现有技术中获取的地震数据的准确性不高的技术问题，达到了有效消除ISS地震数据中的邻炮干扰，提高地震数据准确性的技术效果。附图说明此处所说明的附图用来提供对本专利技术的进一步理解，构成本申请的一部分，并不构成对本专利技术的限定。在附图中：图1是本专利技术实施例的ISS地震数据邻炮干扰压制方法的方法流程图；图2是本专利技术实施例的时域矢量中值滤波压制邻炮干扰的方法流程图；图3是本专利技术实施例模拟可控震源ISS同步激发地震数据采集得到的正演共炮点道集地震本文档来自技高网...

【技术保护点】
一种可控震源独立同步激发地震数据邻炮干扰压制方法，其特征在于，包括：在可控震源独立同步激发ISS地震数据中分选出非共炮点道集地震数据；对分选出的非共炮点道集地震数据进行矢量中值滤波，得到邻炮干扰压制后的非共炮点道集地震数据；从邻炮干扰压制后的非共炮点道集地震数据中分选出共炮点道集地震数据。

【技术特征摘要】
1.一种可控震源独立同步激发地震数据邻炮干扰压制方法，其特征在于，包括：在可控震源独立同步激发ISS地震数据中分选出非共炮点道集地震数据；对分选出的非共炮点道集地震数据进行矢量中值滤波，得到邻炮干扰压制后的非共炮点道集地震数据；从邻炮干扰压制后的非共炮点道集地震数据中分选出共炮点道集地震数据；其中，对分选出的非共炮点道集地震数据进行矢量中值滤波，包括：在时域对分选出的非共炮点道集地震数据进行矢量中值滤波；或者，在频域对分选出的非共炮点道集地震数据进行矢量中值滤波；或者，在时域与频域对分选出的非共炮点道集地震数据进行矢量中值滤波；其中，在频域对分选出的非共炮点道集地震数据进行矢量中值滤波包括：按照以下公式对分选出的非共炮点道集地震数据进行傅里叶变换，得到频域的非共炮点道集地震数据：X(ω)i＝FFT[X(t)i]其中，X(t)i为时域的非共炮点道集地震数据，FFT[]为一维的傅里叶变换，ω为角频率，X(ω)i为傅里叶变换后得到的频域的非共炮点道集地震数据；按照以下公式对经过傅里叶变换后的非共炮点道集地震数据，按地震数据道空间分布进行矢量中值滤波：其中，X(ω)m为对频域的非共炮点道集地震数据进行矢量中值滤波后输出的结果，ω为角频率，m、i和j为地震数据道的顺序号，n为进行矢量中值滤波运算的道数，||.||p为矢量Lp范数；将经过矢量中值滤波后的频域的非共炮点道集数据进行反傅里叶变换，得到邻炮干扰压制后的时域的非共炮点道集地震数据。2.如权利要求1所述的方法，其特征在于，在时域对分选出的非共炮点道集地震数据进行矢量中值滤波，包括：沿时间轴按预定长度将非共炮点道集地震数据分成段；将各段非共炮点道集地震数据作为矢量，按以下公式按地震数据道空间分布对各段非共炮点道集地震数据进行矢量中值滤波，得到邻炮干扰压制后的非共炮点道集地震数据：其中，Xm,k∈{Xi,k|i＝1,...,n}，j＝1,...,n，Xm,k为输出的结果，k为对非共炮点道集地震数据分段后各段的顺序号，Xi,k为分段后各段的非共炮点道集地震数据，m、i和j为地震数据道的顺序号，n为进行矢量中值滤波的道数，||.||p为矢量Lp范数。3.如权利要求2所述的方法，其特征在于，所述预定长度包括：一个子波或者一个波形对应时间轴上的长度。4.如权利要求1至3中任一项所述的方法，其特征在于，非共炮点道集包括以下至少之一：共接收点道集、共偏移距道集、共中心点道集、十字排列道集。5.一种可控震源独立同步激发地震数据邻炮干扰...

【专利技术属性】
技术研发人员：张翊孟，汪长辉，黄明忠，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京;11