[发明专利]电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法

摘要

本发明涉及一种电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法，是通过将多源地震记录即共炮点道集转换到共接收点道集后，通过Radon正变换后在Radon域里分离有效信号与串扰噪声，去噪后再进行Randon反变换，实现了多源地震记录中包含串扰噪声的随机噪声压制。经试验，本发明公开的一种电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法能够实现在多源地震勘探中压制随机延迟带来的串扰噪声,该算法处理数据快，信号细节保护好，在不增加地震勘探数据采集成本情况下，能够有效改善并行震源地震勘探数据质量。

主权项

一种电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法，包括以下步骤：a、在勘探区域内，采用多个震源构成一组并行震源，采用近同时激发方法，令各震源并行工作，获得共炮点地震记录集，此时各地震记录含有来自其它炮点的串扰地震信号；b、进行共炮点道集到共检波点道集的转换，此时并行震源的串扰噪声表现为随机噪声形式；c、针对一个共接收点道集数据进行抛物线Radon正变换，如公式U＝LT(LLT+λE)‑1d  (1)这里d是t‑x域共接收点地震记录，U是τ‑q域数据，其中L=ejωq1x12ejωq1x22...ejωq1xnx2ejωq2x12ejωq2x22...ejωq2xnx2............ejωqnqx12ejωqnqx22...ejωqnqxnx2---(2)]]>L称为变换算子，包括nq行，nx列，L的第i行第k列对应的元素i＝1,2…nq，k＝1,2…nx，j为虚数单位，ω为角频率，LT为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、完成所有共检波点道集到共炮点道集转换。