三维航磁正演计算及简单球状体反演.docVIP

三维航磁正演计算及简单球状体反演

摘要

在重磁勘探领域中，正演问题是反演问题的基础，本文在MATLAB软件上实现了地下球状磁性体的两种公式的正演计算——体积分公式法和解析公式法，对它们的结果进行了对比。对于实现体积分公式法，先把地下半空间均匀剖分成一系列小正方体单元，任意不规则地下磁性体都能近似为一些小正方体单元的组合。用每个小正方体单元的中心点坐标代入体积分公式，可以算出所有小正方体单元在观测点引起的磁异常，再累加求和，进而能得到整个不规则磁性体在观测点引起的磁异常值。接着采用控制变量法，探究在不同的地下半空间剖分步径、不同的磁倾角、不同测量高度的情况下两种公式正演计算结果的差别，