[发明专利]一种平面地磁异常场向下延拓递归余弦变换法

摘要

本发明公开了一种平面地磁异常场向下延拓递归余弦变换法，属于地磁位场数据处理领域。具体步骤为：平面地磁异常场测量数据的预处理；平面地磁异常场测量数据的网格化处理；利用递归余弦变换计算平面地磁异常场的余弦变换波谱；根据向下延拓精度的要求设定一个小量ε，选取阻尼因子α的值；由观测面余弦变换波谱和平面地磁异常场向下延拓的补偿因子直接计算延拓平面地磁异常场场的波谱；对波谱进行行向和列向的一维递归反余弦变，得到延拓值。仅需一步波谱变换和反变换能稳定地进行地磁异常场的大距离向下延拓，减少了向下延拓算法的计算量；递归算法计算量小；提高了平面地磁异常场一步波数域补偿的大距离向下延拓的精度。

主权项

一种平面地磁异常场向下延拓递归余弦变换法，其特征在于，包括如下步骤(1)平面地磁异常场测量数据的预处理，选取网格化参数，对平面地磁异常场测量数据进行预处理，剔除超出合理范围的粗大测量误差；(2)平面地磁异常场测量数据的网格化处理，利用径向基函数插值法，将观测平面上的离散地磁异常场测量数据进行网格化，获得观测平面上的网格化地磁异常测量数据ΔT(nx,ny,z0)，nx＝1,2,L,Mx，ny＝1,2，L,My；(3)利用递归余弦变换由ΔT(nx,ny,z0)计算平面地磁异常场的余弦变换波谱ΔTC(ux,uy,z0)；(4)根据向下延拓精度的要求设定一个小量ε，选取阻尼因子α的值，0＜α＜1，通常取α为0.1；(5)由观测面地磁异常场的余弦变换波谱ΔTC(ux,uy,z0)和平面地磁异常场向下延拓的一步波数域补偿因子直接计算延拓平面地磁异常场场的波谱ΔTC(ux,uy,z0+Δz)；(6)对ΔTC(ux,uy,z0‑Δz)进行行向和列向的一维递归反余弦变，得到延拓平面上的地磁异常场的延拓值ΔT(nx,ny,z0+Δz)。