复杂磁性体井中磁场三分量计算方法、装置、设备及介质

1.本发明属于磁性体磁场数值模拟技术领域，特别涉及一种任意磁化率分布复杂形状磁性体井中磁场三分量计算方法、装置、设备及介质。


背景技术：

2.相比地面和航空磁法勘探，井中磁法勘探因距离场源体更近，能更好反映场源体位置等信息。伴随井中磁力观测设备的发展，井中磁场测量日渐成为一种重要的勘探方法。为适应地面-井中磁场三分量数据联合反演需求，井中磁场三分量高效计算成为亟待解决的问题。相比地面和航空磁场三分量计算，因井中观测点位于地下场源内部，场源点与观测点位于同一位置时格林函数是奇异的，严重影响计算精度。
3.对于任意磁化率分布、复杂形状的磁性体井中磁场三分量计算，一般采用的计算思路是将地下空间剖分成许多小棱柱体，将小棱柱体产生的磁场进行累加逼近复杂磁性体井中磁场，这种方法的最大缺点是计算速度慢。剖分方式和计算方法共同决定了井中磁场三分量计算的效率和精度。计算效率和计算精度是一对矛盾体，累加方法存在的最大问题是不能同时保证计算效率和计算精度，无法满足地面-井中磁场数据联合反演成像、人机交互建模和解释对计算效率的需求。
4.因此，寻找一种计算效率高、同时能保证计算精度的磁性体井中磁场三分量计算方法具有重要的现实意义。


技术实现要素：

5.针对目前现有磁性体井中磁场三分量计算方法计算精度低的问题，同时考虑满足地面-井中磁场数据联合反演成像、人机交互建模和解释对计算效率的需求，本发明提出了复杂磁性体井中磁场三分量计算方法、装置、设备及介质，其能够适应于任意磁化率分布、复杂形状的磁性体井中磁场三分量快速、高精度计算。
6.为实现上述技术目的，本发明提出的技术方案为：
7.一方面，本发明提供一种复杂磁性体井中磁场三分量计算方法，包括：
8.确定目标区域，构建内部包含所述目标区域的棱柱体模型；对棱柱体模型进行网格剖分，剖分成若干个小棱柱体单元，给每个小棱柱体单元的磁化率赋值；
9.采用三维离散卷积法，计算各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量；
10.基于各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量，得到整个棱柱体模型在井中观测点处产生的井中磁场三分量。
11.进一步地，本发明中目标区域中每个小棱柱体单元的磁化率为常值，不同小棱柱体单元的磁化率取值不同，以此刻画任意磁化率分布、复杂形状的磁性体，位于空气部分的小棱柱体单元的磁化率值设为零，以此刻画起伏地形。
12.进一步地，所有小棱柱体单元在井中观测点(x0,y0,zn)处产生的井中磁场三分量
包括9个分量，分别是b
x1
(x0,y0,zn)、b
x2
(x0,y0,zn)、b
x3
(x0,y0,zn)、b
y1
(x0,y0,zn)、b
y2
(x0,y0,zn)、b
y3
(x0,y0,zn)、b
z1
(x0,y0,zn)、b
z2
(x0,y0,zn)和b
z3
(x0,y0,zn)，计算公式分别为：
[0013][0014][0015][0016][0017][0018][0019][0020][0021][0022]
其中各井中观测点的水平坐标为x0，y0，各井中观测点的垂直坐标为zn，n＝1,2,...,n，n为井中观测点个数，p＝1,2,...,p，p为棱柱体模型x方向小棱柱体单元的剖分个数，q＝1,2,...,q，q为棱柱体模型y方向小棱柱体单元的剖分个数，r＝1,2,...,r，r为棱柱体模型z方向小棱柱体单元的剖分个数，m
x
(ξ
p
,ηq,ζr)、my(ξ
p
,ηq,ζr)和mz(ξ
p
,ηq,ζr)分别表示中心点坐标位置为(ξ
p
,ηq,ζr)的小棱柱体单元的磁化强度的x分量、y分量和z分量；h1(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h2(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h3(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h4(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h5(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h6(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h7(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h8(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h9(x
0-ξ
p
,y
0-ηq,z
n-ζr)为9个加权系数。
[0023]
进一步地，将棱柱体模型中所有小棱柱体单元在井中观测点(x0,y0,zn)处产生的井中磁场三分量的各个分量分别进行累加，得到整个棱柱体模型在井中观测点(x0,y0,zn)处产生的井中磁场三分量，即
[0024]bx
(x0,y0,zn)＝b
x1
(x0,y0,zn) b
x2
(x0,y0,zn) b
x3
(x0,y0,zn)
[0025]by
(x0,y0,zn)＝b
y1
(x0,y0,zn) b
y2
(x0,y0,zn) b
y3
(x0,y0,zn)
[0026]bz
(x0,y0,zn)＝b
z1
(x0,y0,zn) b
z2
(x0,y0,zn) b
z3
(x0,y0,zn)
[0027]
另一方面，本发明提供一种复杂磁性体井中磁场三分量计算装置，包括：
[0028]
第一模块，用于确定目标区域，构建内部包含所述目标区域的棱柱体模型；对棱柱体模型进行网格剖分，剖分成若干个小棱柱体单元，给每个小棱柱体单元的磁化率赋值；
[0029]
第二模块，用于采用三维离散卷积法，计算各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量；
[0030]
第三模块，用于基于各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量，得到整个棱柱体模型在井中观测点处产生的井中磁场三分量。
[0031]
另一方面，本发明提供一种计算机设备，包括存储器和处理器，存储器存储有计算机程序，处理器执行计算机程序时实现以下步骤：
[0032]
确定目标区域，构建内部包含所述目标区域的棱柱体模型；对棱柱体模型进行网格剖分，剖分成若干个小棱柱体单元，给每个小棱柱体单元的磁化率赋值；
[0033]
采用三维离散卷积法，计算各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量；
[0034]
基于各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量，得到整个棱柱体模型在井中观测点处产生的井中磁场三分量。
[0035]
再一方面，本发明还提供一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
[0036]
确定目标区域，构建内部包含所述目标区域的棱柱体模型；对棱柱体模型进行网格剖分，剖分成若干个小棱柱体单元，给每个小棱柱体单元的磁化率赋值；
[0037]
采用三维离散卷积法，计算各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量；
[0038]
基于各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量，得到整个棱柱体模型在井中观测点处产生的井中磁场三分量。
[0039]
与现有技术相比，本发明的优点在于：
[0040]
(1)本发明提出的棱柱体模型表示方法简单、灵活，很容易刻画任意磁化率分布、复杂形状磁性体以及起伏地形；
[0041]
(2)本发明能够实现任意磁化率分布、任意测井测点分布情况下复杂形状磁性体井中磁场三分量的快速、高精度计算，可以满足大规模地面-井中磁场数据联合反演、人机交互建模和解释的需求；
[0042]
(3)在进行大规模计算时，本发明不但计算效率和计算精度高，并且所需计算机内存小。
附图说明
[0043]
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图示出的结构获得其他的附图。
[0044]
图1是本发明一实施例的流程图；
[0045]
图2是本发明一实施例中复杂磁性体模型示意图；
[0046]
图3为本发明一实施例中球状磁性体模型图；
[0047]
图4为本发明一实施例中磁场三分量计算值，其中(a)为磁场三分量其b
x
分量计算值，(b)为磁场三分量其by分量计算值，(c)为磁场三分量其bz分量计算值；
[0048]
图5为本发明一实施例中磁场三分量理论值，其中(a)为磁场三分量其b
x
分量理论值，(b)为磁场三分量其by分量理论值，(c)为磁场三分量其bz分量理论值；
[0049]
图6为本发明一实施例中磁场三分量计算值与理论值的相对误差，其中(a)为磁场三分量其b
x
分量计算值与理论值的相对误差，(b)为磁场三分量其by分量计算值与理论值的相对误差，(c)为磁场三分量其bz分量计算值与理论值的相对误差；
[0050]
图7为本发明一实施例的结构示意图。
具体实施方式
[0051]
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面将以附图及详细叙述来清楚说明本发明所揭示内容的精神，任何所属技术领域技术人员在了解本发明内容的实施例后，当可由本发明内容所教示的技术，加以改变及修饰，其并不脱离本发明内容的精神与范围。本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
[0052]
参照图1，本发明一实施例中，提供一种复杂磁性体井中磁场三分量计算方法，包括：
[0053]
(s1)确定目标区域，构建内部包含所述目标区域的棱柱体模型；对棱柱体模型进行网格剖分，剖分成若干个小棱柱体单元，给每个小棱柱体单元的磁化率赋值。
[0054]
(s2)采用三维离散卷积法，计算各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量。
[0055]
(s3)基于各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量，得到整个棱柱体模型在井中观测点处产生的井中磁场三分量。
[0056]
本发明的(s1)中，实质就是复杂磁性体模型表示过程，具体如下：
[0057]
确定目标区域，构建内部包含所述目标区域的棱柱体模型。确定棱柱体模型在x，y，z方向的起始位置，使得目标区域(包含起伏地形)完全嵌入在该棱柱体模型中。
[0058]
根据实际问题需求，将棱柱体模型进行网格剖分，将其均匀划分成许多规则小棱柱体单元，如图2所示，确定小棱柱体单元的几何尺寸δx、δy、δz，其中δx、δy、δz分别表示小棱柱体单元其x方向、y方向、z方向的边长。
[0059]
根据目标区域的磁化率分布，对每个小棱柱体单元的磁化率进行赋值，目标区域中每个小棱柱体单元的磁化率为常值，不同小棱柱体单元的磁化率取值不同，以此刻画任意磁化率分布、复杂形状的磁性体，位于空气部分的小棱柱体单元的磁化率值设为零，以此刻画起伏地形。
[0060]
至此，完成复杂磁性体模型表示。
[0061]
本发明的(s2)中，所述井中观测点有多个，给定各井中观测点的位置坐标(x0,y0,zn)，井中观测点(x0,y0,zn)的水平坐标为x0，y0，井中观测点的垂直坐标为zn，n＝1,2,...,n，n为井中观测点个数。
[0062]
本发明的(s2)中，计算各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量，过程包括：
[0063]
(s2.1)计算中心点坐标位置为(ξ
p
,ηq,ζr)的小棱柱体单元的磁化强度的x分量m
x
(ξ
p
,ηq,ζr)、y分量my(ξ
p
,ηq,ζr)和z分量mz(ξ
p
,ηq,ζr)。
[0064]
(s2.2)计算加权系数h1(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h2(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h3(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h4(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h5(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h6(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h7(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h8(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h9(x
0-ξ
p
,y
0-ηq,z
n-ζr)。
[0065][0066][0067]
[0068][0069][0070][0071]
[0072][0073][0074]
其中μ0表示真空磁导率；arctan()表示反余切函数运算符，ln()表示自然对数运算符；x1＝ξ
p-0.5δx-x0，x2＝ξ
p
0.5δx-x0，y1＝η
q-0.5δy-y0，y2＝ηq 0.5δy-y0，z1＝ζ
r-0.5δz-zn，z2＝ζr 0.5δz-zn；δx、δy、δz分别表示小棱柱体单元其x方向、y方向、z方向的边长。
[0075]
(s2.3)采用三维离散卷积法，计算所有小棱柱体单元在井中观测点处产生的井中磁场三分量。
[0076]
所有小棱柱体单元在井中观测点(x0,y0,zn)处产生的井中磁场三分量包括9个分量，分别是b
x1
(x0,y0,zn)、b
x2
(x0,y0,zn)、b
x3
(x0,y0,zn)、b
y1
(x0,y0,zn)、b
y2
(x0,y0,zn)、b
y3
(x0,y0,zn)、b
z1
(x0,y0,zn)、b
z2
(x0,y0,zn)和b
z3
(x0,y0,zn)。基于以下计算公式，采用采用三维离散卷积法，分别计算所有小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量：
[0077][0078][0079][0080]
[0081][0082][0083][0084][0085][0086]
其中p＝1,2,...,p，p为棱柱体模型x方向小棱柱体单元的剖分个数，q＝1,2,...,q，q为棱柱体模型y方向小棱柱体单元的剖分个数，r＝1,2,...,r，r为棱柱体模型z方向小棱柱体单元的剖分个数。
[0087]
本发明的(s3)中，将棱柱体模型中所有小棱柱体单元在井中观测点(x0,y0,zn)处产生的井中磁场三分量各个分量分别进行累加，得到整个棱柱体模型在井中观测点(x0,y0,zn)处产生的井中磁场三分量，即
[0088]bx
(x0,y0,zn)＝b
x1
(x0,y0,zn) b
x2
(x0,y0,zn) b
x3
(x0,y0,zn)
[0089]by
(x0,y0,zn)＝b
y1
(x0,y0,zn) b
y2
(x0,y0,zn) b
y3
(x0,y0,zn)
[0090]bz
(x0,y0,zn)＝b
z1
(x0,y0,zn) b
z2
(x0,y0,zn) b
z3
(x0,y0,zn)
[0091]
在本发明一实施例中，对于步骤(s2.1)中，计算中心点坐标位置为(ξ
p
,ηq,ζr)的小棱柱体单元的磁化强度的x分量m
x
(ξ
p
,ηq,ζr)、y分量my(ξ
p
,ηq,ζr)和z分量mz(ξ
p
,ηq,ζr)，包括：
[0092]
(s2.1.1)根据地球主磁场模型igrf，计算中心点坐标位置为(ξ
p
,ηq,ζr)的小棱柱体单元的地球主磁场的x分量t
x
(ξ
p
,ηq,ζr)、y分量ty(ξ
p
,ηq,ζr)和z分量tz(ξ
p
,ηq,ζr)。
[0093]
(s2.1.2)根据中心点坐标位置为(ξ
p
,ηq,ζr)的小棱柱体单元的磁化率κ(ξ
p
,ηq,ζr)，计算中心点坐标位置为(ξ
p
,ηq,ζr)的小棱柱体单元的磁化强度的x分量m
x
(ξ
p
,ηq,ζr)、y分量my(ξ
p
,ηq,ζr)和z分量mz(ξ
p
,ηq,ζr)：
[0094]mx
(ξ
p
,ηq,ζr)＝κ(ξ
p
,ηq,ζr)
·
t
x
(ξ
p
,ηq,ζr)
[0095]my
(ξ
p
,ηq,ζr)＝κ(ξ
p
,ηq,ζr)
·
ty(ξ
p
,ηq,ζr)
[0096]mz
(ξ
p
,ηq,ζr)＝κ(ξ
p
,ηq,ζr)
·
tz(ξ
p
,ηq,ζr)
[0097]
在本发明另一实施例中，步骤(s2)中的三维离散卷积法，其步骤包括：
[0098]
(1)将加权系数h(x
0-ξ
p
,y
0-ηq,z
n-ζr)(表示h1(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h2(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h3(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h4(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h5(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h6(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h7(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h8(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h9(x
0-ξ
p
,y
0-ηq,z
n-ζr)中任意一个)排列成维数为p
×q×
(r n-1)的三维数组t，数组元素t
i,j,k
与加权系数h(x
0-ξ
p
,y
0-ηq,z
n-ζr)存在关系
[0099]
t
i,j,k
＝h(x
0-ξi,y
0-ηj,z
1-zk)
[0100]
式中，zk＝z1 (k-n)δz，k＝1,2,
…
,r n-1，i＝1,2,
…
,p，j＝1,2,
…
,q。p为棱柱体
模型x方向小棱柱体单元的个数，q为棱柱体模型y方向小棱柱体单元的个数，r为棱柱体模型z方向小棱柱体单元的个数，n为井中观测点个数，δz分别表示小棱柱体单元其z方向的边长。
[0101]
(2)将磁化强度m(ξ
p
,ηq,ζr)(p＝1,2,...,p，q＝1,2,...,q，r＝1,2,...,r，表示m
x
(ξ
p
,ηq,ζr)、my(ξ
p
,ηq,ζr)、mz(ξ
p
,ηq,ζr)中任一个)排列成维数为p
×q×
r的三维数组m，三维数组m中的数组元素m
p,q,r
与密度值存在关系
[0102]mp,q,r
＝m(ξ
p
,ηq,ζr)
[0103]
将三维数组m补零扩展成为p
×q×
(r n-1)的三维数组m
ext
，三维数组m位于三维数组m
ext
的前p行、前q列、前r层；
[0104]
(3)计算
[0105]
式中，fft3()表示三维快速傅里叶变换；
[0106]
(4)计算
[0107]
式中，“.*”表示对应元素相乘运算；
[0108]
(5)计算
[0109]
式中，ifft3()表示三维快速傅里叶反变换；
[0110]
(6)提取三维数组f
ext
的前p行、前q列、前r层元素，构成三维数组f，即为三维离散卷积计算结果。
[0111]
本发明是一个有机整体，构建内部包含所述目标区域的棱柱体模型，采用三维离散卷积法，计算所有小棱柱体单元在井中观测点处产生的井中磁场三分量，得到整个棱柱体模型在井中观测点处产生的井中磁场三分量。本发明创新性的提出了加权系数h1(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h2(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h3(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h4(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h5(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h6(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h7(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h8(x
0-ξ
p
,y
0-ηq,z
n-ζr)、h9(x
0-ξ
p
,y
0-ηq,z
n-ζr)的计算方法，利用上述方法计算得到9个加权系数并结合采用三维离散卷积快速计算方法，实现了意磁化率分布复杂磁性体井中磁场三分量在计算精度和计算效率上的统一。本发明解决了现有井中磁性体磁场三分量计算方法计算精度低、计算速度慢，无法满足地面-井中磁场数据联合反演成像、人机交互建模和解释需求的问题。
[0112]
下面对本发明提供的复杂磁性体井中磁场三分量计算方法的精度和效率进行检验。
[0113]
为了说明本发明所提出的方法用于计算任意磁化率分布情况下复杂形状磁性体井中磁场三分量量时的效率和精度，设计了如图3所示的复杂磁性体模型。
[0114]
磁化率为常值的棱柱体区域内嵌一个磁化率均匀的球体，棱柱体和球体的中心重合。棱柱体范围为：x方向从-1000m到1000m，y方向从-1000m到1000m，z方向从0m到1000m(z轴向下为正)；球体半径为400m。棱柱体磁化率为0，球体的磁化率为0.1si。将棱柱体剖分成200
×
200
×
400个大小相同的单元，井中观测点水平坐标为(4,10)，观测点垂直坐标从20m到980m等间隔采样，观测点个数为385。
[0115]
井中磁场三分量算法利用matlab语言编程实现，运行程序所用的个人台式机配置为：cpu为i7-2620，主频为2.7ghz，内存为32gb，四核八线程。计算一个磁场三分量运行所需时间约为10秒，由此可见本发明方法效率很高。磁场三分量计算值如图4所示，磁场三分量
理论值如图5所示，从形态上看，两者是一致的，其中图4中(a)为磁场三分量其b
x
分量计算值，(b)为磁场三分量其by分量计算值，(c)为磁场三分量其bz分量计算值，图5中(a)为磁场三分量其b
x
分量理论值，(b)为磁场三分量其by分量理论值，(c)为磁场三分量其bz分量理论值。相对误差由理论值减去计算值得到差值的绝对值除以理论值得到(图6所示，其中(a)为磁场三分量其b
x
分量计算值与理论值的相对误差，(b)为磁场三分量其by分量计算值与理论值的相对误差，(c)为磁场三分量其bz分量计算值与理论值的相对误差)，对相对误差进行统计，统计结果由表1给出，可知算法精度很高。
[0116]
表1磁场三分量理论值和计算值相对误差统计
[0117] 最大值最小值均值b
x
0.0343.95
×
10-6
0.0012by0.0431.10
×
10-6
0.0019bz1.479.89
×
10-7
0.0057
[0118]
本发明另一实施例中，提供一种复杂磁性体井中磁场三分量计算装置，包括：
[0119]
第一模块，用于确定目标区域，构建内部包含所述目标区域的棱柱体模型；对棱柱体模型进行网格剖分，剖分成若干个小棱柱体单元，给每个小棱柱体单元的磁化率赋值；
[0120]
第二模块，用于采用三维离散卷积法，计算各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量；
[0121]
第三模块，用于基于各小棱柱体单元在井中观测点处产生的井中磁场三分量的各个分量，得到整个棱柱体模型在井中观测点处产生的井中磁场三分量。
[0122]
上述各模块功能的实现方法，可以采用前述各实施例中相同的方法实现，在此不再赘述。
[0123]
在本实施例中，提供了一种计算机设备，该计算机设备可以是服务器，其内部结构图可以如图7所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的数据库用于存储样本数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现上述实施例中复杂磁性体井中磁场三分量计算方法的步骤。
[0124]
本领域技术人员可以理解，图7中示出的结构，仅仅是与本技术方案相关的部分结构的框图，并不构成对本技术方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
[0125]
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，该存储器存储有计算机程序，该处理器执行计算机程序时实现上述实施例中复杂磁性体井中磁场三分量计算方法的步骤。
[0126]
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述实施例中复杂磁性体井中磁场三分量计算方法的步骤。
[0127]
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机
可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本技术所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(rom)、可编程rom(prom)、电可编程rom(eprom)、电可擦除可编程rom(eeprom)或闪存。易失性存储器可包括随机存取存储器(ram)或者外部高速缓冲存储器。作为说明而非局限，ram以多种形式可得，诸如静态ram(sram)、动态ram(dram)、同步dram(sdram)、双数据率sdram(ddrsdram)、增强型sdram(esdram)、同步链路(synchlink)dram(sldram)、存储器总线(rambus)直接ram(rdram)、直接存储器总线动态ram(drdram)、以及存储器总线动态ram(rdram)等。
[0128]
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0129]
以上所述实施例仅表达了本技术的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本技术构思的前提下，还可以做出若干变形和改进，这些都属于本技术的保护范围。因此，本技术专利的保护范围应以所附权利要求为准。