基于凸集投影的重力同时填充扩边和去噪方法

实测重力数据经常因工区不规则导致数据存在空缺。在波数域进行重力数据处理和转换前必须对原始数据进行填充和扩边处理。数据中存在的高频噪声往往是造成数据后续处理不稳定的主要原因。常规重力数据处理一般将数据的填充、去噪和扩边这三个步骤独立进行。为此,将这三个问题统一考虑,提出一种基于改进凸集投影方法的规则网格重力数据同时填充、去噪和扩边迭代方法。首先通过计算并拟合重力数据的径向平均功率谱确定迭代的最终截止波数,然后利用数据低于截止波数的谱对数据同时进行填充、扩边和去噪,直至达到预设的迭代次数。对理论重力模型和阿富汗实测均衡重力数据进行了扩充和去噪实验,结果表明,所提方法原理简单、操作方便,填充和扩边结果的数据衔接处光滑无畸变,且具有较好的填充和去噪效果,处理结果优于经典的最小曲率和Kriging插值法及小波去噪和余弦扩边组合法。     

基于凸集投影的重力数据扩充下延一体化方法

重力勘探的工区经常是不规则的，有时重力数据还存在空白区。在进行波数域向下延拓处理前，需对数据空白区进行填充。同时，为了提高处理精度，还需要将数据预扩边至快速Fourier变换所需要的长度。常规重力向下延拓方法一般仅考虑向下延拓本身，或将数据填充、扩边与向下延拓这三个步骤独立进行。本文统一考虑这三个不适定问题，提出一种基于凸集投影原理的重力数据填充、扩边、下延一体化方法，即采用Gerchberg-Saxton迭代法对数据进行插值和扩边，直至达到预定的迭代次数，然后下延。利用该方法和经典的插值、扩边、下延组合方法对理论模型和航空实测重力数据进行对比实验，结果表明该方法原理简单、操作方便，扩边结果光滑无畸变，插值和下延结果精度较高，优于经典组合方法。