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

权利要求书

1.一种高精度空地井不同空间重磁数据变换方法，其特征在于,包括如下步骤：

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

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

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

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

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

其中，T₀(x,y,0)和T(x,y,h)表示高差为h的两个观测面上的异常，当h0时表示向上延拓，当h0时表示向下延拓，对公式(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)相加后可得：

所述步骤S2中对初始异常数据T_i^h进行技术的方程式为：

且采用Laplace方程来计算异常的二阶和四阶垂直导数：

所述步骤S3中计算T₁与原始异常数据T的公式为：

ΔT₁(x，y，0)＝T(x，y，0)-T₁(x，y，0) (8)；

利用公式(6)将△T₁向下延拓高度h，可以得到：

其中，ΔT₁(x,y,h)为异常ΔT₁(x,y,0)下延h后异常，ΔT₁(x,y,-h)为异常ΔT₁(x,y,0)上延h后异常，ΔT₁(x,y,h)作下界面异常的一个修正量，则下界面异常可改写为：

T₁(x，y，h)＝T(x，y，h)+ΔT₁(x，y，h) (10)；

所述步骤S4中计算过程，直至ΔT_m(x，y，0)的均方差小于给定值ε，因此下界面h上的最终异常为：

T_m(x，y，h)＝T(x，y，h)+ΔT₁(x，y，h)+…+ΔT_m(x，y，h) (11)。