混叠噪声压制方法及装置制造方法及图纸

本发明专利技术提供一种混叠噪声压制方法及装置，其中，该方法包括以下步骤：将混叠的共炮点域数据转换为共偏移距域混叠数据；对共偏移距域混叠数据进行混叠噪声压制；若压制混叠噪声后的共偏移距域混叠数据的信噪比小于预设值，则对所述压制混叠噪声后的共偏移距域混叠数据进行混叠计算，得到混叠数据，将所述混叠的共炮点域数据更新为所述混叠数据，重复执行以上步骤，直到压制混叠噪声后的共偏移距域混叠数据的信噪比大于等于预设值；将信噪比大于等于预设值的压制混叠噪声后的共偏移距域混叠数据，转换为压制混叠噪声后的共炮点域数据。本发明专利技术实现了对混叠的共炮点域数据中混叠噪声的压制，提高了地震数据的信噪比，提升了混叠噪声压制效果。噪声压制效果。噪声压制效果。

【技术实现步骤摘要】
混叠噪声压制方法及装置


[0001]本专利技术涉及地球物理勘探
，特别涉及一种混叠噪声压制的方法及装置。

技术介绍

[0002]近年来，随着高效采集技术的发展和在实际中的广泛应用，使得陆上宽方位和高密度采集成为可能。高效采集技术的应用不仅大幅度提高了生产效率，能够获得高覆盖次数的地震资料，而且极大地改善了地震资料品质。混叠采集技术采用多组可控震源独立激发，相互间没有等待时间，对可控震源数量没有明显限制，震源越多，效率越高。根据混叠采集的特点，相同的排列接收了不同震源在不同的位置以不同的时间激发的能量，在共炮点道集上具有相干性；而在共偏移距道集、共接收点道集或共中心点道集上，只有来自主激发点的能量具有相干性，来自邻炮干扰的能量则都表现为随机性强振幅或尖脉冲干扰，所以多组震源同时激发会在原始单炮记录上产生很强的邻炮干扰，会产生严重的混叠噪声，严重降低了混叠采集的地震数据的信噪比。目前混叠噪声的压制方法主要包括基于信号域的滤波方法。而基于信号域的滤波方法原理简单，运算效率快，但是压制效果不佳。

技术实现思路

[0003]本专利技术实施例提供了一种混叠噪声压制方法，用以提升混叠噪声压制效果，采集的地震数据的信噪比同时保护有效信号，该方法包括：
[0004]将混叠的共炮点域数据转换为共偏移距域混叠数据；
[0005]对共偏移距域混叠数据进行混叠噪声压制；
[0006]若压制混叠噪声后的共偏移距域混叠数据的信噪比小于预设值，则对所述压制混叠噪声后的共偏移距域混叠数据进行混叠计算，得到混叠数据，将所述混叠的共炮点域数据更新为所述混叠数据；
[0007]重复执行以上步骤，直到压制混叠噪声后的共偏移距域混叠数据的信噪比大于等于预设值；
[0008]将信噪比大于等于预设值的压制混叠噪声后的共偏移距域混叠数据，转换为压制混叠噪声后的共炮点域数据。
[0009]本专利技术实施例还提供了一种混叠噪声压制装置，用以提升混叠噪声压制效果，该装置包括：
[0010]共偏移距域混叠数据转换模块，用于将混叠的共炮点域数据转换为共偏移距域混叠数据；
[0011]混叠噪声压制模块，用于对共偏移距域混叠数据进行混叠噪声压制；
[0012]混叠数据更新模块，用于若压制混叠噪声后的共偏移距域混叠数据的信噪比小于预设值，则对所述压制混叠噪声后的共偏移距域混叠数据进行混叠计算，得到混叠数据，将所述混叠的共炮点域数据更新为所述混叠数据；
[0013]迭代触发模块，用于重复触发所述共偏移距域混叠数据转换模块、混叠噪声压制
模块和混叠数据更新模块执行，直到压制混叠噪声后的共偏移距域混叠数据的信噪比大于等于预设值；
[0014]压制混叠噪声后的共炮点域数据模块，用于将信噪比大于等于预设值的压制混叠噪声后的共偏移距域混叠数据，转换为压制混叠噪声后的共炮点域数据。
[0015]本专利技术实施例还提供了一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述所述方法。
[0016]本专利技术实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述所述方法的计算机程序。
[0017]本专利技术实施例中，将混叠的共炮点域数据转换为共偏移距域混叠数据，对共偏移距域混叠数据进行混叠噪声压制，若压制混叠噪声后的共偏移距域混叠数据的信噪比小于预设值，则对所述压制混叠噪声后的共偏移距域混叠数据进行混叠计算，得到混叠数据，将所述混叠的共炮点域数据更新为所述混叠数据，重复执行以上步骤，最终可以得到一信噪比大于预设值的压制混叠噪声后的共偏移距域混叠数据，进而确定了压制混叠噪声后的共炮点域数据。本专利技术实施例实现了对混叠的共炮点域数据中混叠噪声的压制，提高了地震数据的信噪比，提升了混叠噪声压制效果。
附图说明
[0018]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0019]图1是本专利技术实施例提供的一种混叠噪声压制方法的流程示意图；
[0020]图2是本专利技术实施例提供的一种原始不含混叠噪声的模拟数据的示意图；
[0021]图3是本专利技术实施例提供的一种含混叠噪声的模拟数据的示意图；
[0022]图4是本专利技术实施例提供的一种压制混叠噪声后的数据的示意图；
[0023]图5是本专利技术实施例提供的一种去除的混叠干扰数据的示意图；
[0024]图6是本专利技术实施例提供的一种不含混叠噪声的实际数据的示意图；
[0025]图7是本专利技术实施例提供的一种模拟的含混叠噪声的实际数据的示意图；
[0026]图8是本专利技术实施例提供的一种压制混叠噪声后的实际数据的示意图；
[0027]图9是本专利技术实施例提供的一种压制的混叠干扰数据的示意图；
[0028]图10是本专利技术实施例提供的一种混叠噪声压制装置的结构示意图。
具体实施方式
[0029]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0030]目前混叠噪声的压制方法主要包括基于信号域的滤波方法。而基于信号域的滤波
方法原理简单，运算效率快，但是压制效果不佳。因此，本专利技术设计了一种迭代的奇异值约束的混叠噪声压制方法。奇异值分解(SVD，Singular Value Decomposition)方法本质上是一种正交分解方法，它利用二维空间中信号相关性压制不相关的噪声，提高地震资料的信噪比。
[0031]奇异值分解(SVD，Singular Value Decomposition)是利用特征值或奇异值作为正交基在信号空间正交分解的特征以增强相干能量，压制干扰信号。设M
×
N的地震数据D，其中M为道数，N为采样点数，则该地震数据的SVD分解表示为：
[0032][0033]其中，u
i
为矩阵DD
T
第i个特征向量，v
i
为矩阵D
T
D第i个特征向量，U＝[u1,...,u
M
]M
×
M
和V＝[v1,...,v
N
]N
×
N
均为正交矩阵，矩阵∑由D的奇异值σ
i
构成，r为地震数据D的秩，奇异值σ1≥σ2...≥σ
r
≥0由大到小排列在其主对角线上。
[0034]地震数据经SVD分解后，根据能量大小可分为若干个奇异值，大奇异值主要反映原始记录本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种混叠噪声压制方法，其特征在于，包括：将混叠的共炮点域数据转换为共偏移距域混叠数据；对共偏移距域混叠数据进行混叠噪声压制；若压制混叠噪声后的共偏移距域混叠数据的信噪比小于预设值，则对所述压制混叠噪声后的共偏移距域混叠数据进行混叠计算，得到混叠数据，将所述混叠的共炮点域数据更新为所述混叠数据；重复执行以上步骤，直到压制混叠噪声后的共偏移距域混叠数据的信噪比大于等于预设值；将信噪比大于等于预设值的压制混叠噪声后的共偏移距域混叠数据，转换为压制混叠噪声后的共炮点域数据。2.如权利要求1所述的方法，其特征在于，对共偏移距域混叠数据进行混叠噪声压制，包括：对共偏移距域混叠数据进行奇异值分解，得到共偏移距域混叠数据的奇异值向量；从共偏移距域混叠数据中选取一仅包括混叠噪声的矩形窗，对矩形窗中的数据进行奇异值分解，得到混叠噪声的奇异值向量；根据共偏移距域混叠数据的奇异值向量和混叠噪声的奇异值向量，确定压制混叠噪声后的共偏移距域混叠数据。3.如权利要求2所述的方法，其特征在于，按如下公式，对共偏移距域混叠数据进行奇异值分解，得到共偏移距域混叠数据的奇异值向量：d
offset
＝U
bl V
bl
W
Tbl
其中，V
bl
为共偏移距域混叠数据的奇异值向量；d
offset
为共偏移距域混叠数据；U
bl
为共偏移距域混叠数据的左特征向量；W
Tbl
为共偏移距域混叠数据的右特征向量。4.如权利要求2所述的方法，其特征在于，从共偏移距域混叠数据中选取一仅包括混叠噪声的矩形窗，包括：从共偏移距域混叠数据中，选取距离激发点最近的共偏移距域混叠数据，确定有效波最先到达的时间；根据有效波最先到达的时间，从距离激发点最近的共偏移距域混叠数据中，选取一仅包括混叠噪声的矩形窗。5.如权利要求2所述的方法，其特征在于，按如下公式，对选取的矩形窗中的混叠数据进行奇异值分解，得到混叠噪声的奇异值向量：d
noise
＝U
noise
V
noise
W
Tnoise
其中，V
noise
为混叠噪声的奇异值向量；d
noise
为选取的仅含有混叠噪声的数据；U
noise
为混叠噪声的左特征向量；W
Tnoise
为混叠噪声的右特征向量。6.如权利要求2所述的方法，其特征在于，根据共偏移距域混叠数据的奇异值向量和混叠噪声的奇异值向量，确定压制混叠噪声后的共偏移距域混叠数据，包括：取混叠噪声的奇异值向量的最大值为混叠噪声的约束值；采用硬阈值函数，根据混叠噪声的约束值，对共偏移距域混叠数据的奇异值向量进行更新计算，得到压制混叠噪声后的共偏移距域混叠数据的奇异值向量；根据压制混叠噪声后的共偏移距域混叠数据的奇异值向量，确定压制混叠噪声后的共
偏移距域混叠数据。7.如权利要求6所述的方法，其特征在于，按如下公式，取混叠噪声的奇异值向量的最大值为混叠噪声的约束值：λ＝max(V
noise
)其中，λ为混叠噪声的约束值；V
noise
为混叠噪声的奇异值向量。8.如权利要求6所述的方法，其特征在于，按如下公式，采用硬阈值函数，根据混叠噪声的约束值对共偏移距域混叠数据的奇异值向量进行更新计算，得到压制混叠噪声后的共偏移距域混叠数据的奇异值向量：V1＝max(0,V
bl-λ)
×
V
bl
其中，V1为压制混叠噪声后的共偏移距域混叠数据的奇异值向量；λ为混叠噪声的约束值；V
bl
为共偏移距...

【专利技术属性】
技术研发人员：董烈乾，张慕刚，汪长辉，骆飞，周大同，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：