基于主成分分析的并行震源地震勘探资料强工频噪声压制方法技术

本发明专利技术公开一种基于主成分分析的并行震源地震勘探资料强工频噪声压制方法，并行震源地震勘探是在勘探区域内不同炮点处布置多个震源，采用近同时激发技术采集资料，可提高地震勘探生产效率和地震资料质量。在人文活动区并行震源激发采集得到的原始地震记录中经常混合了强工频噪声，对于干扰道根据相空间理论构造Hankel矩阵，并对其进行主成分分解与恢复有用信号。经验证，本发明专利技术能够实现对并行震源采集得到的地震勘探资料进行有效的强工频噪声压制，与传统的强工频噪声压制方法相比，处理数据快，不需要精确估计强工频噪声的频率，无相位偏移，信号保护性好，能够有效改善并行震源勘探的地震资料质量，降低了并行震源地震勘探资料处理成本。

【技术实现步骤摘要】
基于主成分分析的并行震源地震勘探资料强工频噪声压制方法
：本专利技术涉及一种基于主成分分析的并行震源地震勘探资料强工频噪声压制方法，该方法用于并行震源地震勘探时强工频噪声的压制。
技术介绍
：并行震源地震勘探在采集过程中往往受到强工频噪声的影响，采集得到并行地震资料中常混有强工频噪声，因此得到的地震资料品质不高，进而影响了后期地震资料解释和偏移成像质量。而关于并行震源地震勘探资料强工频噪声压制方法的文献尚少见，目前其他领域对强工频噪声抑制的方法主要有频域法和时域法两大类，频域法如陷波法，小波滤波法等都在频率域进行工频压制，虽然该类方法简单，方便，但当有用信号频率和工频的频率互相混叠时，这类方法会对有用信号造成一定程度的损害，同时也会产生一定的相移。时域法如正余弦逼近法，工频回归相减法，自适应滤波法等是将工频噪声表示为以振幅，频率，相位为变量的函数，再通过正余弦函数逼近法来估计强工频噪声，但这类方法需要对强工频噪声频率进行精确估计，同时这类方法对多个工频谐波存在时的情况需要多次处理，从而算法效率较低。可见，上述两类方法都不适应于并行震源地震勘探资料强工频噪声的压制。
技术实现思路
：本专利本文档来自技高网...

【技术保护点】
基于主成分分析的并行震源地震勘探资料强工频噪声压制方法，其特征在于，包括以下步骤：a、对并行地震勘探资料进行频谱分析，根据强工频噪声在频谱中幅值异常突出和频带窄，且频带在50Hz附近的特点，识别含有强工频噪声的干扰道；b、根据相空间重构理论，对xi(l)构造Hankel矩阵

【技术特征摘要】
1.基于主成分分析的并行震源地震勘探资料强工频噪声压制方法，其特征在于，包括以下步骤：a、对并行地震勘探资料进行频谱分析，根据强工频噪声在频谱中幅值异常突出和频带窄，且频带在50Hz附近的特点，识别含有强工频噪声的干扰道；b、根据相空间重构理论，对xi(l)构造Hankel矩阵该矩阵的行数记为m，列数记为n，若N为偶数，则令m＝N/2，n＝N/2，若N为奇数，则令m＝(N/2)+1，n＝(N/2)+1；c、计算Hi的协方差矩阵Γi，如公式HiT为Hi的转置矩阵，“·”表示矩阵乘法；d、利用奇异值分解法，计算协方差矩阵Γi的特征值矩阵Λi和特征向量矩阵Ri，则存在公式Γi＝Ri·Λi·RiT(3)其中Λi为由大到小排列的特征值矩阵，Ri为各个特征值对应的特征向量矩阵，RiT为Ri的转置矩阵，且满...

【专利技术属性】
技术研发人员：姜弢，汪彦龙，
申请(专利权)人：吉林大学，
类型：发明
国别省市：吉林,22