基于高分辨率Radon变换的VSP波场分离方法

讨论了时间域线性Radon变换、频率域最小平方Radon变换和频率域高分辨率Radon变换等方法的基本原理。Radon变换在频率域通常采用最小平方反演的方式来实现,但由于最小平方法对于不同的p值(视慢度)范围采用相同的阻尼因子,得到的解不够精确。高分辨率Radon变换对这一缺陷进行了改进,即通过对不同的p值范围采用不同的加权值来提高解的稀疏性和精确性。具体方法是:在反演迭代过程中,根据前一次迭代的结果,通过Bayes原理将加权矩阵与前一次迭代的结果联系起来,得到新的加权矩阵;然后求解这个加权矩阵方程,得到频率域的稀疏解。设计了一个3层水平地层的地质模型,利用高分辨率Radon变换进行了VSP波场分离。τ-p域的结果表明,能量得到了很好的收敛。在分离出来的上、下行波波场记录上,波形恢复得很好,上、下行波波场的相互干扰被完全消除。