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

说明书

技术领域

本发明属于数学三维模型技术领域，涉及基于全空间正则化下延数据的视磁化强 度三维反演方法。

背景技术

传统的视磁化强度反演方法，是利用在观测面采集的实际磁力数据来反推地下异 常体的磁性参数以及空间赋存状态的一种地球物理手段，这种方法的基本思想就是对地下 待反演地质体进行网格剖分，通过不停的反演迭代直至达到精度要求，来计算出所有剖分 网格单元的视磁化强度。虽然很容易找到满足数据拟合空间的解模型，但是由于重磁场的 体积效应、观测数据有限性及反演问题的欠定性等几种因素的存在，反演结果往往产生多 解性。为了减小这种多种性的影响，提高反演的精度，前人做了大量的研究并提出了一系列 方法，如已知井资料定量约束、概率成像定性约束等等。这些方法虽然一定程度上克服了多 解性的影响，提高了反演结果的精度，但并没有解决上文提到的导致反演多解性的根本问 题——数据的有限性。这种有限性主要体现在两个方面：(1)磁化强度是地下不同深度场的 累加效应，而在数据采集时，只能测量到在地表观测到的数据；(2)当场源埋深较大时，地表 观测的异常信号较弱，不能有效突出异常信号。

发明内容

本发明的目的在于提供基于全空间正则化下延数据的视磁化强度三维反演方法， 解决了目前三维反演过程中由于数据有限，导致反演精度低的问题。

本发明所采用的技术方案是按照以下步骤进行：

步骤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,]]>δ0={ΣmMΣnN|Smn*-Smn|2}12]]>为偏差泛函，为包含误差的混合值，S_mn为不含误差的有效值；