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

权利要求书

1.一种电磁式可控震源抛物线Radon变换压制并行震源串扰噪声方法，包括以下步骤：

a、在勘探区域内，采用多个震源构成一组并行震源，采用近同时激发方法，令各震源并行工作，获得共炮点地震记录集，此时各地震记录含有来自其它炮点的串扰地震信号；

b、进行共炮点道集到共检波点道集的转换，此时并行震源的串扰噪声表现为随机噪声形式；

c、针对一个共接收点道集数据进行抛物线Radon正变换，如公式

U＝L^T(LL^T+λ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称为变换算子，包括n_q行，n_x列，L的第i行第k列对应的元素i＝1,2…n_q，k＝1,2…n_x，j为虚数单位，ω为角频率，L^T为L共轭转置算子，λ为阻尼因子，一般取0.1≤λ≤1，E为单位矩阵，n_x为检波点的采集道数；n_q采用如下方式定义，任取满足的q值,其中x_max为最大偏移距，Δx为道间距，f为地震波最小频率，令Δq为曲率采样间隔，任取满足的Δq值，定义,q_i＝-|q|+i·Δq,x_k＝10k；

d、根据U,绘制τ-q域图像，其有效信号呈规律的点状或拖长点状，其余为噪声信号，保留U中有效信号，其余数据置0，得到U'；

e、对U'进行抛物线Radon反变换，具体如公式

d'＝LU'(3)

d'为反变换后的t-x域共检波点地震记录；

f、对所有共检波点道集重复步骤c～e，直至完成全部共检波点地震记录处理；

g、完成所有共检波点道集到共炮点道集转换。