[发明专利]一种二度体磁场数值计算方法

摘要

本发明提出了一种二度体磁场数值计算方法，包括复杂二度体模型表示以及矩形模型磁场计算两大部分，本发明在特定的复杂二度体模型表示方式条件下，建立矩形二度体磁场叠加模型，根据设计的一种加权系数计算公式，采用一维离散卷积快速算法，实现了磁场计算在效率和精度上的统一。本发明解决了目前现有二度体磁场数值计算方法不能同时保证计算效率和计算精度，无法满足大规模磁场数据磁化率精细反演成像、人机交互建模和解释的需求问题。

主权项

1.一种二度体磁场数值计算方法，其特征在于，包括以下步骤：第一步复杂二度体模型表示：根据二度体展布范围，建立包含所有目标区域的矩形模型，确定矩形模型在x，z方向的起始位置，使得包含起伏地形的目标区域完全嵌入在该矩形模型中；然后将该矩形模型均匀划分成若干个规则小矩形，确定小矩形的几何尺寸Δx，Δz；最后根据二度体的磁化率分布，对每个小矩形磁化率进行赋值，每个小矩形磁化率为常值，不同矩形磁化率取值不同，以此刻画复杂截面形状、任意磁化率分布的二度体；将位于空气部分的小矩形的磁化率值设为零，以此刻画起伏地形；步骤二、矩形模型磁场计算：步骤一中给出的矩形模型，其磁场计算公式为式中，(xm,z0)表示观测点坐标，z0为常值；L表示矩形模型其z方向的小矩形剖分个数；M表示矩形模型其x方向的小矩形剖分个数；(ξp,ζr)表示编号为(p,r)的小矩形几何中心坐标；mx(ξp,ζr)和mz(ξp,ζr)分别表示编号为(p,r)的小矩形的磁化强度的x分量和z分量；hx(xm‑ξp,z0‑ζr)和hz(xm‑ξp,z0‑ζr)分别表示对应磁化强度的x分量和z分量的加权系数。