[发明专利]基于三边滤波的混合震源数据分离方法

权利要求书

1.一种基于三边滤波的混合震源数据分离方法，包括以下步骤：

a、对实际混合采集地震数据进行预处理得到实际混合震源地震数据U；

b、首先假设初始分离单炮炮集记录P_i＝0(i＝0)；

c、在频率域通过混合震源编码Γ_bl将分离单炮炮集记录混合成混合震源记录P_bl：

Γ_bl＝e^-iωt

P_bl＝ifft[fft(P_i)×Γ_bl]

其中ω为频率，t为混合震源中单炮的延迟时间，fft与ifft分别表示快速傅立叶变换与快速傅立叶逆变换；

d、计算实际混合震源地震数据U与分离单炮炮集记录合成的混合震源记录P_bl之间的差值，即剩余混合信号U_dif＝U-P_bl；

e、在频率域通过混合震源编码Γ_bl对U_dif做伪分离，得到伪分离后的单炮地震数据P′：

P′=ifft[fft(Udif)(ΓblHΓbl)-1ΓblH]]]>

其中H表示共轭转置，伪分离后的单炮地震数据P′含有大量的混合噪声呈相干分布；

f、将伪分离后的单炮地震数据由共炮域转换到共偏移距域，得到伪分离后的共偏移距域地震记录，此时共偏移距域内有效信号呈连续分布，混合噪声以随机脉冲形式分布；

g、使用三边滤波方法去除共偏移距域地震记录中的混合噪声；

h、将三边滤波后的共偏移距域地震记录转换回共炮点域，得到剩余混合信号的分离结果P_dif，那么实际混合震源地震数据U的分离结果P_i+1＝P_dif+P_i；

i、取i＝i+1，重复过程b～h直到最终得到满足精度要求的结果。

2.按照权利要求1所述的基于三边滤波的混合震源数据分离方法，其特征在于，g步骤去除共偏移距域地震记录中混合噪声具体方法包括以下步骤：

(1)任意选取共偏移距域地震记录中的一个点x＝(x₁,x₂)，该点灰度值为u_x，令Ω＝Ω_x(N)是以x点为中心，N为半滤波器宽的邻域；

(2)计算点x的ROAD值

ROAD(x)＝r₁(x)+r₂(x)+…+r_m(x)

其中r_k(x)是y∈Ω，|u_y-u_x|的值按照从小到大排序的第k个值，半滤波器宽N＝1，m＝4；

(3)计算点x的权系数：

ω(x,y)＝ω_S(x,y)ω_R(x,y)^1-J(x,y)ω_I(x,y)^J(x,y)

其中表示距离权重，表示灰度权重，表示脉冲权重，

J(x,y)=1-e-(ROAD(x)+ROAD(y)2)2/2σJ2]]>

当x点和y点都不为脉冲噪声点时，ROAD(x)和ROAD(y)值都比较小，此时J(x,y)≈0，脉冲权重ω_I对权系数不起作用，当x点和y点至少有一个脉冲噪声点时，ROAD(x)+ROAD(y)值远远大于0，那么J(x,y)≈1，脉冲权重ω_I对图像平滑起主要作用，

其中σ_I、σ_J取值范围区间为(0，1)，随着i增加σ_I、σ_J取相同的值分别为：0.2、0.23、0.25、0.26、0.27、0.27、0.27；

(4)根据步骤(3)中给定的σ_I、σ_J的值，计算三边滤波修复后的图像点

u~x=Σy∈Ωω(x,y)uyΣy∈Ωω(x,y).]]>