[发明专利]基于全空间正则化下延数据的视磁化强度三维反演方法

权利要求书

1.基于全空间正则化下延数据的视磁化强度三维反演方法，其特征在于按照以下步骤 进行反演：

步骤1：采用三维正则化下延算子实现不同深度的延拓得到稳定的磁异常全空间场，采 用的三维正则化算子为：

Hα=e-2πfz/{(1+αeβ(f1-f)λx)·(1+αeβ(f-f2)λy)},]]>其中，α是正则化参数，β是滤波系数， f₁、f₂为通带范围，f为波数，z为深度，λ_x、λ_y为基波波长，u，v为二维频域坐标 系，构建目标方程求取最佳的α、β、f₁、f₂值，目标方程为

φ(α,f1,f2)=Σm=1MΣn=1N|Smn*(1+αeβ(f1-f)λx)(1+αeβ(f-f2)λy)-Smn*|2-δ02=0,]]>

为偏差泛函，为包含误差的混合值，S_mn为不含误差的有 效值；

步骤2：求解目标方程的极小值，得到最小值对应的β、α、f₁、f₂，代入正则化下延算子从 而计算出不同深度的正延场值E_h(x,y,z)＝E₀(x,y,z)*h_α，E₀为地表观测值，h_α为H_α的反傅氏 变换；

步骤3：构建全空间场反演目标函数；

选取反演目标函数为σ2=Σk=1N(fk-f‾)2,σ2y=Σk=1N(Ek-E‾)2]]>式中，为全空间 场值E_k(E_k＝E₀+E_h)的平均值，该观测场包含地表观测值E₀和步骤1下延得到的不同深度的场 值E_h，f_k表示第k个观测点的观测值和正演计算值之间的差；为f_k的平均值，其数学表 达式如下：E‾=1NΣk=1NEk,fk=Ek-E^k,]]>

f‾=1NΣk=1Nfk;]]>

步骤4：采用多维搜索黄金分割算法求解反演目标函数的极小值，提高反演速度。

2.按照权利要求1所述基于全空间正则化下延数据的视磁化强度三维反演方法，其特 征在于：所述步骤1中β取值在0.5-2.5之间，以0.1为单位。