[发明专利]一种CRP道集去噪、校正与叠加的方法

摘要

本发明公开了一种CRP道集去噪、校正与叠加的方法，特征是包括以下步骤：a、将初始CRP道集进行叠加求取参考道，并可对初始CRP道集做前期去噪处理，前期去噪处理包括频域滤波、小波变换、中值滤波或均值滤波；b、开设滑动时窗，并在时窗口内依据上述参考道求取校正量，然后根据校正量校正初始CRP道集得到一级CRP道集；c、将一级CRP道集作为初始CRP道集，依次重复步骤a、b，获取二级CRP道集；d、将二级CRP道集应用SVD去噪，获取三级CRP道集；e、将三级CRP道集数据进行叠加，得到成果剖面。本发明不但能够使CRP道集同相轴校正到水平，使其信噪比大大提高，而且能够使叠加后剖面同相轴连续性增强，信噪比提高、弱信号得到加强。

主权项

1.一种CRP道集去噪、校正与叠加的方法，其特征在于包括以下步骤：a、将初始CRP道集进行叠加求取参考道，并可对初始CRP道集做前期去噪处理，前期去噪处理包括频域滤波、小波变换、中值滤波或均值滤波；将初始CRP道集设定为G_CRP(t,i),i＝1,2,...,N，其中N为道集中的道数、t为地震道时间、i表示道集中第i道，将参考道设定为S_r(t)，Sr(t)=[Σi=1NGCRP(t,i)]/N]]>b、开设滑动时窗，并在时窗口内依据上述参考道求取校正量，然后根据校正量校正初始CRP道集；在校正量的求取过程中，建立目标函数f(t,i)，f(t,i)=Στ=t-T/2τ+T/2[W(τ)Sr(τ)]·[W(τ+Δt(t,i))GCRP(τ+Δt(t,i),i)]→max]]>其中W(t)为窗函数，通过对该目标函数的最优化求解得到校正量Δt(t,i)，利用校正量Δt(t,i)校正初始CRP道集得到一级CRP道集，将一级CRP道集设定为S~CRP1(t,i)=GCRP(t+Δt(t,i),i)]]>c、将一级CRP道集作为初始CRP道集，依次重复步骤a、b，重复次数为1～2次，获取二级CRP道集，将二级CRP道集设定为d、将二级CRP道集应用SVD去噪，获取三级CRP道集；在SVD去噪过程中，首先将二级CRP道集进行奇异值分解，奇异值分解表达式为，G~CRP2=UEVT]]>其中E为对角矩阵，奇异值分布于对角线上，选取优势能量对应的奇异值组成新的对角矩阵E_p，并重组二级CRP道集获得三级CRP道集，将三级CRP道集设定为重构奇异值表达式为，G~CRP3=UEpVT]]>其中重构奇异值个数选取范围从[1～1]至[1～N/2]；e、将三级CRP道集数据进行叠加，得到成果剖面，将成果剖面设定为S_stack,j(t)，Sstack,j(t)=[Σi=1NG~CRP,j3(t,i)]/N,j=1,2,...,M]]>其中j表示第j个道集，M为道集的个数。