[发明专利]一种高精度空地井不同空间重磁数据变换方法

说明书

本发明公开了一种高精度空地井不同空间重磁数据变换方法，所述高精度空地井不同空间重磁数据变换方法通过泰勒级数展开式证明了向下延拓计算可通过向上延拓和水平导数之和来完成，向上延拓和水平导数的计算是稳定的，受噪音影响较小，且为了弥补导数阶次为有效项造成的误差，采用一种异常向下延拓的水平导数迭代法，有效地增强了计算结果的稳定性和准确性，该发明通过泰勒级数展开式证明了向下延拓计算可通过向上延拓和水平导数之和来完成，向上延拓和水平导数的计算是稳定的，受噪音影响较小，且为了弥补导数阶次为有效项造成的误差，采用一种异常向下延拓的水平导数迭代法，有效地增强了计算结果的稳定性和准确性。

技术领域

本发明属于地球科学技术领域，更具体地说，尤其涉及一种高精度空地井不同空间重磁数据变换方法。

背景技术

针对矿产、油气资源潜力区域，为了了解地下物质的不同深度上物质的分布，往往在一个地区开展航空、地面、井中重磁测量，从而利用不同空间的数据联合完成地下地质体的深度及物性反演工作，可有效地提高解释结果的精度和可靠性。但是在实际勘探中，不同空间数据的范围和比例尺是不一样的，从而大多时候需要利用其它空间的数据进行换算填补缺失的数据，所以不同空间重磁数据变换技术是进行空-地-井位场数据联合反演的前提，其变换精度也是制约反演精度的主要因素。

不同空间位场数据的获得是采用向上和向下延拓技术来实现，向上延拓采用频率域公式来进行，具有计算速度快、受噪声干扰下小的优点；向下延拓若采用频率域技术进行是不稳定的，且会明显增大噪声的干扰，使测得曲线发生强烈地跳动，从而造成测量精确度低，测得数据不准确的问题，越来越不能满足人们的需求，不利于广泛的推广和普及。

发明内容

本发明的目的是为了解决现有技术中存在的缺点，而提出的一种高精度空地井不同空间重磁数据变换方法。

为实现上述目的，本发明提供如下技术方案：一种高精度空地井不同空间重磁数据变换方法，包括如下步骤：

S1、首先在需要勘测地点的地片面计算再地表面的原始异常数据T，并同时勘测向上延伸的异常数据T_-h；

S2、然后利用泰勒展开式方程来计算异常的二阶和四阶垂直导数，并且利用测得的导数来和泰勒展开式计算向下延伸的初始异常数据

S3、然后根据公式计算下界面向上延伸地平面后的异常数据T₁，并把T₁与原始异常数据T进行比较得出误差数据ΔT₁；

S4、如果ΔT₁均方差大于ε时，则重复步骤2和步骤3来重复计算过程，直到测得的ΔT_m均方差小于ε时输出下延伸异常数据

优选的，所述步骤S1中观测面上重磁异常之间的关系式为：

其中，T₀(x，y，0)和T(x，y，h)表示高差为h的两个观测面上的异常，当h＞0时表示向上延拓，当h＜0时表示向下延拓，对公式(1)进行Fourier变换可得到

其中，表示异常T₀(x，y，0)的Fourier谱，(k_x，k_y)是分别表示x和y方向上的波数，弥为延拓因子。

优选的，所述步骤S2中泰勒展开式计算深度为h界面上的异常表达式为：

其中，T(x，y，h)是观测面h上的异常，h为延拓高度，T(x，y，0)是观测面上异常，向上延拓运算的泰勒展开式可表示为：

其中，T(x，y，-h)为观测面-h上的异常，将式(3)和(4)相加后可得：