[发明专利]基于空间波数域的强磁介质的磁感应强度计算方法

权利要求书

1.基于空间波数域的强磁介质的磁感应强度计算方法，其特征在于，包括：

获取包含异常体在内的三维目标区域，建立包含三维目标区域的初始棱柱体模型；

对初始棱柱体模型沿x、y、z方向分别进行等间隔剖分，在x、y、z方向上得到一系列的节点；根据磁化率分布数据对每个节点进行磁化率赋值，磁化率为标量，用表示，其中在异常体部分的节点，根据异常体的磁化率值赋给每个节点；无异常部分的节点的磁化率为0；

根据棱柱体模型以及给定的x方向和y方向上的高斯参数，计算x和y方向离散偏移波数；

根据地球主磁场模型，计算每个节点处的地球主磁场强度，将其作为空间域背景场磁场强度；

根据空间域背景场磁场强度、空间域异常场磁场强度，得到磁化强度的计算模型；

利用二维傅里叶变换将空间域异常场磁位和磁化强度满足的三维拉普拉斯方程转为空间波数域一维常微分方程；

基于空间波数域一维常微分方程，并结合设定的空间波数域异常场磁位需满足的边界条件，将空间波数域异常场磁位满足的边值问题模型转化为等价的变分问题模型；

通过求解变分问题模型，得到空间波数域异常场磁位；

基于空间波数域异常场磁位，求得空间波数域异常场磁场强度，并通过反傅里叶变换将空间波数域异常场磁位、空间波数域异常场磁场强度转换为空间域异常场磁位、空间域异常场磁场强度；

空间域总磁场强度为空间域异常场磁场强度和空间域背景场磁场强度之和，对空间域总磁场强度运用紧算子进行收紧，得到收紧后的空间域总磁场强度；

判断当前是否满足迭代收敛条件，如满足则输出当前计算得到的收紧后的空间域总磁场强度所对应的空间域异常场磁位、空间域异常场磁场强度；

根据输出的空间域异常场磁场强度求解得到空间域磁感应强度。

2.根据权利要求1所述的基于空间波数域的强磁介质的磁感应强度计算方法，其特征在于，高斯参数包括：

x方向的高斯点个数N_x，区间[-1,1]上高斯点t_a、高斯系数A_a，其中，a=1,2,...,N_x；

y方向的高斯点个数N_y，区间[-1,1]上高斯点t_b、高斯系数A_b，其中，b=1,2,...,N_y。