1.强磁性二度体的勘探方法,其特征在于,包括以下步骤:设置地下目标区域、强磁性二度体的展布范围、强磁性二度体的磁化率分布数据、地面观测高度和水平观测点坐标;
根据地下目标区域、强磁性二度体的展布范围、强磁性二度体的磁化率分布数据,计算地下目标区域强磁性二度体的磁场;
根据地面观测高度、水平观测点坐标、地下目标区域强磁性二度体的磁场,计算水平观测点位置处强磁性二度体的磁场;
判断计算得到的水平观测点位置处强磁性二度体的磁场与仪器测量得到的水平观测点位置处强磁性二度体的磁场相同,如果相同,则将强磁性二度体的磁化率分布数据作为所述强磁性二度体的实际磁化率分布数据以用于勘探所述强磁性二度体。
2.根据权利要求1所述的强磁性二度体的勘探方法,其特征在于,计算地下目标区域强磁性二度体的磁场,包括:
(a)根据所述地下目标区域和所述强磁性二度体的展布范围建立初始二维矩形模型;
(b)将所述初始强磁性二度体的展布范围均匀划分成多个单元矩形;
(c)根据所述强磁性二度体的磁化率分布数据对每个所述单元矩形的磁化率进行赋值,得到对应所述初始强磁性二度体的目标二维矩形模型;
(d)根据所述目标二维矩形模型计算得到其模型加权系数;
(e)根据地球主磁场模型,计算所述目标二维矩形模型中每个所述单元矩形中心处的地球主磁场;
(f)将每个所述单元矩形中心处的地球主磁场作为其磁场初始值;
(g)根据所述目标二维矩形模型、所述磁场初始值以及模型加权系数计算得到空间域异常磁场;
(h)根据所述磁场初始值和所述空间域异常磁场计算得到总磁场;
(i)设置迭代收敛条件,判断所述总磁场是否满足迭代收敛条件,若所述总磁场满足迭代收敛条件,将所述总磁场作为地下目标区域强磁性二度体的磁场,若所述总磁场不满足所述给定迭代收敛条件,将所述总磁场作为所述磁场初始值,并重复执行步骤(g)至步骤(i),直至满足所述给定迭代收敛条件。
3.根据权利要求2所述的强磁性二度体的勘探方法,其特征在于,步骤(d)中,模型加权系数ω(x1‑ξm,z1‑ζn)包括ωx(xi‑ξm,zj‑ζn)和ωz(xi‑ξm,zj‑ζn)两个分量,分别为:其中:目标二维矩形模型所在区域即地下目标区域强磁性二度体的磁场所在区域,地下目标区域强磁性二度体的磁场所在区域与地下目标区域强磁性二度体的磁场观测点区域重合,各单元矩形的中心同时作为地下目标区域强磁性二度体的磁场观测点区域的观测点,(xi,zj)和(ξm,ζn)分别表示目标二维矩形模型中的单元矩形的中心坐标以及地下目标区域强磁性二度体的磁场观测点区域的观测点坐标,i=1,2,…,Nx,j=1,2,…,Nz,m=1,
2,…,Nx,n=1,2,…,Nz,Nx和Nz分别为二维矩形模型x和z方向上划分的单元矩形的数量,Δx和Δz分别为单元矩形在二维矩形模型x和z方向的尺寸,arctan表示反正切运算,ln表示对数运算。
4.根据权利要求3所述的强磁性二度体的勘探方法,其特征在于,步骤(g)中,空间域异常磁场Ha(xi,zj)由x方向上和z方向上的两个分量构成,如下:其中:mx(ξm,ζn)和mz(ξm,ζn)分别表示(ξm,ζn)处的空间域磁化强度M(xi,zj)的x分量和z(0) (0)
分量,M(xi,zj)=χ(xi,zj)H (xi,zj),H (xi,zj)=Hb(xi,zj),Hb(xi,zj)为单元矩形中心处的地球主磁场。
5.根据权利要求4所述的强磁性二度体的勘探方法,其特征在于,步骤(g)中,采用一种二维离散卷积快速算法,计算空间域异常磁场Ha(xi,zj)的两个分量Hax(xi,zj)和Haz(xi,zj)。
6.根据权利要求4所述的强磁性二度体的勘探方法,其特征在于,步骤(h)中,总磁场H(1)
(xi,zj)=Ha(xi,zj)+Hb(xi,zj)。
7.根据权利要求6所述的强磁性二度体的勘探方法,其特征在于,步骤(i)中设置的迭代收敛条件为:
其中,ε0为期望的数值精度。
8.根据权利要求3、4、5、6或7所述的强磁性二度体的勘探方法,其特征在于,计算水平观测点位置处强磁性二度体的磁场,包括:计算观测高度磁场的观测点加权系数ω(Xp‑ξm,Z0‑ζn):其中:ωx(Xp‑ξm,Z0‑ζn)和ωz(Xp‑ξm,Z0‑ζn)分别表示观测高度磁场的观测点加权系数的两个分量,Z0表示观测高度,Xp表示水平观测点的坐标,P=1,2,…,N0,N0表示观测点个数,水平观测点坐标间隔为Δx;
计算水平观测点位置处强磁性二度体的磁场Ha(Xp,Z0)其x方向上和z方向上两个分量,如下:
9.根据权利要求8所述的强磁性二度体的勘探方法,其特征在于,采用一维离散卷积快速算法,反复调用Nz次,计算水平观测点位置处强磁性二度体的磁场Ha(Xp,Z0)其x方向上和z方向上两个分量。
10.强磁性二度体的勘探系统,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1所述强磁性二度体的勘探方法的步骤。