电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN106154330A(43)申请公布日2016.11.23(21)申请号201610439306.6(22)申请日2016.06.20(71)申请人吉林大学地址130012吉林省长春市前进大街2699号(72)发明人姜弢汪彦龙徐学纯贾海青岳永高邹艳艳晁云峰(74)专利代理机构长春吉大专利代理有限责任公司22201代理人王立文(51)Int.Cl.G01V1/36(2006.01)权利要求书1页说明书3页附图2页(54)发明名称电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法(57)摘要本发明涉及一种电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法，是通过将多源地震记录即共炮点道集转换到共接收点道集后，通过Radon正变换后在Radon域里分离有效信号与串扰噪声，去噪后再进行Randon反变换，实现了多源地震记录中包含串扰噪声的随机噪声压制。经试验，本发明公开的一种电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法能够实现在多源地震勘探中压制随机延迟带来的串扰噪声,该算法处理数据快，信号细节保护好，在不增加地震勘探数据采集成本情况下，能够有效改善并行震源地震勘探数据质量。CN106154330ACN106154330A权利要求书1/1页1.一种电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法，包括以下步骤：a、在勘探区域内，采用多个震源构成一组并行震源，采用近同时激发方法，令各震源并行工作，获得共炮点地震记录集，此时各地震记录含有来自其它炮点的串扰地震信号；b、进行共炮点道集到共检波点道集的转换，此时并行震源的串扰噪声表现为随机噪声形式；c、针对一个共接收点道集数据进行抛物线Radon正变换，如公式U＝LT(LLT+λE)-1d(1)这里d是t-x域共接收点地震记录，U是τ-q域数据，其中L称为变换算子，包括nq行，nx列，L的第i行第k列对应的元素TL为L共轭转置算子，λ为阻尼因子，一般取0.1≤λ≤1，E为单位矩阵，nx为检波点的采集道数；nq采用如下方式定义，任取满足的q值,其中xmax为最大偏移距，Δx为道间距，f为地震波最小频率，令Δq为曲率采样间隔，任取满足的Δq值，定义,qi＝-|q|+i·Δq,xk＝10k；d、根据U,绘制τ-q域图像，其有效信号呈规律的点状或拖长点状，其余为噪声信号，保留U中有效信号，其余数据置0，得到U'；e、对U'进行抛物线Radon反变换，具体如公式d'＝LU'(3)d'为反变换后的t-x域共检波点地震记录；f、对所有共检波点道集重复步骤c～e，直至完成全部共检波点地震记录处理；g、完成所有共检波点道集到共炮点道集转换。2CN106154330A说明书1/3页电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法技术领域：[0001]本发明涉及一种电磁式可控震源并行震源地震勘探时压制串扰噪声的方法，尤其是电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法。背景技术：[0002]多震源激发就是在多个不同位置上激发多个震源并行工作，因此得到的地震记录是一种由多个震源产生的混合地震记录，提高并行震源地震勘探数据质量，但多震源地震数据为并行处理问题。博士论文“电磁式可控震源地震信号检测关键技术研究中”运用多方向中值滤波技术压制并行震源带来的干扰噪声，但该方法往往需要多次迭代才能达到较好信噪比，数据处理速度较慢，信号细节有一定损失，CN104765069A公开的“一种压制同步激发采集领炮干扰方法”提出采用alpha-trimmed矢量中值阈值滤波去除随机噪声，该方法算法容易引入高频干扰，CN103675903A公开的“一种随机去噪多震源地震波场分离方法”没有明确指出压制噪声方法，CN104536034A公开的“多震源并行激发采集与混合地震记录分离方法”适于多震源采用互异控制信号时的噪声压制，而当前大部分并行震源采用相同或相似的控制信号。发明内容：[0003]本发明的目的就在于针对上述现有技术的不足，提供一种电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法。[0004]本发明的主要思想是：[0005]多震源并行激发技术提高了生产效率，但采集得到的地震数据同时混合了来自多炮点的多个并行震源地震信号，这样获得的数据并不能直接使用，本发明是通过将多源地震记录即共炮点道集转换到共接收点道集后，通过Radon正变换后在Radon域里分离有效信号与串扰噪声，去噪后再进行Radon反变换，实现了多源地震记录中包含串扰噪声的随机噪声压制。[0006]本发明是通过以下技术方案实现的：[0007]电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法，包括以下步骤：[0008]a