[发明专利]一种多次波组合自适应衰减方法

摘要

一种多次波组合自适应衰减方法，属于地震多次波衰减领域，具体步骤如下：首先以原始记录和由自由界面多次波预测方法(SRMP)得到的多次波记录作为初始记录，利用同相轴优化追踪技术和短时窗FK视速度滤波法提取和压制多次波同相轴，从而得到准一次波记录和准多次波记录；然后利用扩展的维纳滤波法进一步消除准一次波记录中的残余多次波，并从准多次波记录中还原出原损伤的一次波信息。该方法综合利用了基于同相轴优化追踪的多次波剔除方法和与扩展维纳滤波法各自的优势，其无论是对常见的双曲线型多次波还是其它类型多次波都可进行有效压制，并可最大限度地减少一次波的损伤。

主权项

一种多次波组合自适应衰减方法，其特征在于它具体包括以下步骤：1)基于同相轴优化追踪的多次波匹配衰减过程，包括以下几步：(1)利用自由界面多次波预测方法得到多次波记录m(x,t)，将其进行动校正和叠加处理，得到多次波叠加速度谱，相应的计算公式为$E_m(v,\mathrm{\τ})=\left|\underset{n=1}{\overset{N}{\Σ}}m(x_n,t=\sqrt{{\mathrm{\τ}}^2+x_n^2/v^2})\right|$式中，n为道号，1≤n≤N；x_n为第n道的偏移距，对E_m(v,τ)进行适度平滑后，m(x,t)中叠加速度值v₀、零偏移距时τ₀的双曲线同相轴，将在谱E_m(v,τ)上形成以(v₀,τ₀)为中心极值的团状结构能量；(2)应用等值线追踪方法求出谱E_m(v,τ)上以(v₀,τ₀)为中心极值的团状结构能量的分布范围，并搜索出其极值点位置，然后根据该点的坐标(v₀,τ₀)拟合出时空域中的相应同相轴，即其所经各道的旅行时t_n，表示为$t_n=\sqrt{{{\mathrm{\τ}}_0}^2+{x_n}^2/{v_0}^2}$然后以原追踪的极值位置(v₀,τ₀)为中心创建原始记录d(x,t)的局部速度谱，在一定的速度与时间范围内扫描则可获得该同相轴的准确参量v′₀与τ′₀，从而实现多次波同相轴的优化追踪过程；(3)在得到各多次波同相轴所经地震道的旅行时t_m后，利用短时窗FK视速度滤波法从原始记录d(x,t)和多次波记录m(x,t)中将其消除，得到去除多次波同相轴后的准一次波记录d_p(x,t)与剩余多次波记录m_l(x,t)，而去除掉的多次波同相轴合成为准多次波记录d_m(x,t)与双曲多次波记录m_h(x,t)；2)利用扩展的维纳滤波法针对准多次波记录d_m(x,t)与准一次波记录d_p(x,t)进行再次自适应衰减，以滤出准多次波记录d_m(x,t)中的一次波信息，并进一步消除准一次波记录d_p(x,t)中的残余多次波，相应的处理过程表示为p₁(x,t)＝ELSF[d_m(x,t),m_h(x,t)]p₂(x,t)＝ELSF[d_p(x,t),m_l(x,t)]其中，ELSF表示扩展最小平方滤波算子；p₁(x,t)中包含的是由记录d_m(x,t)中滤出的一次波信息；p₂(x,t)为消除记录d_p(x,t)中残余多次波的结果，最终衰减多次波后的记录如下式所示：p(x,t)＝p₁(x,t)+p₂(x,t)式中，p(x,t)为组合自适应衰减的最终结果。