1.复杂强磁性体的航空勘探方法，其特征在于，包括以下步骤：根据地下目标区域、强磁性体的展布范围、强磁性体的磁化率分布数据，计算地下目标区域强磁性体的磁场，包括：

(a)根据所述地下目标区域和所述强磁性体的展布范围建立初始三维棱柱体模型；

(b)将所述初始三维棱柱体模型均匀划分成多个规则的小棱柱体；

(c)根据所述强磁性体的磁化率分布数据对每个所述小棱柱体的磁化率进行赋值，得到对应所述强磁性体的目标三维棱柱体模型；

(d)根据所述目标三维棱柱体模型计算得到模型磁场加权系数；

对于三维棱柱体模型中的每个小棱柱体，模型磁场加权系数包括6个，分别是：其中：三维棱柱体模型所在区域即地下目标区域强磁性体的磁场所在区域，地下目标区域强磁性体的磁场所在区域与地下目标区域强磁场观测点区域重合，各小棱柱体的几何中心同时作为地下目标区域强磁场观测点区域中的观测点，(xi,yj,zk)和(ξm,ηn,ζl)分别表示三维棱柱体模型中的小棱柱体的几何中心坐标和地下目标区域强磁场观测点区域中的观测点坐标，i＝1,2,…,Nx，j＝1,2,…,Ny，k＝1,2,…,Nz，m＝1,2,…,Nx，n＝1,2,…,Ny，l＝

1,2,…,Nz，此处Nx，Ny和Nz分别为三维棱柱体模型x，y和z方向上小棱柱体的数量，Δx，Δy和Δz分别为小棱柱体在x，y和z方向的尺寸，arctan表示反正切运算，ln表示对数运算，ξ，η，ζ分别表示加权系数计算公式中的积分上下限，上限取ξ＝xi‑ξm+

0.5Δx，η＝yj‑ηn+0.5Δy，ζ＝zk‑ζl+0.5Δz，下限取ξ＝xi‑ξm‑0.5Δx，η＝yj‑ηn‑0.5Δy，ζ＝zk‑ζl‑0.5Δz；

(e)根据所述目标三维棱柱体模型计算得到紧算子；

紧算子包括：

其中：α(xi,yj,zk)和β(xi,yj,zk)表示以(xi,yj,zk)为几何中心坐标的小棱柱体的紧算子，χ(xi,yj,zk)表示以(xi,yj,zk)为几何中心坐标的小棱柱体的磁化率；

(f)根据地球主磁场模型，计算每个所述小棱柱体几何中心处的地球主磁场；

(g)将每个所述小棱柱体几何中心处的所述地球主磁场作为其对应的磁场初始值；

(h)根据所述目标三维棱柱体模型、所述磁场初始值以及模型磁场加权系数计算得到空间域异常磁场；

空间域异常磁场Ha(xi,yj,zk)的三个分量如下：其中：mx(ξm,ηn,ζl)，my(ξm,ηn,ζl)和mz(ξm,ηn,ζl)分别表示小棱柱体的几何中心坐标(ξm,ηn,ζl)处的空间域磁化强度M(xi,yj,zk)的x分量，y分量和z分量；M(xi,yj,zk)＝χ(xi,yj,zk)H(0) (0) (0)

(xi,yj,zk)，H (xi,yj,zk)为以(xi,yj,zk)为几何中心坐标的小棱柱体的磁场初值，H(xi,yj,zk)＝Hb(xi,yj,zk)，Hb(xi,yj,zk)为小棱柱体的几何中心坐标(xi,yj,zk)处的地球主磁场，由地球主磁场模型计算得到；m＝1,2,…,Nx，n＝1,2,…,Ny，l＝1,2,…,Nz；

(1)

根据所述磁场初始值、所述紧算子和所述空间域异常磁场计算得到总磁场H (xi,yj,zk)，为：

(1) (0)

H (xi,yj,zk)＝α(xi,yj,zk)(Ha(xi,yj,zk)+Hb(xi,yj,zk))+β(xi,yj,zk)H (xi,yj,zk)(j)若所述总磁场满足给定的迭代收敛条件，将所述总磁场作为所述地下目标区域强磁性体的所述磁场；

(k)若所述总磁场不满足所述给定迭代收敛条件，将所述总磁场作为所述磁场初始值，并重复执行步骤(h)至步骤(k)；

根据地上观测高度、水平观测点坐标、地下目标区域强磁性体的磁场，计算地上观测高度强磁性体的磁场，包括：

根据所述地上观测高度、水平观测点坐标计算六个地上观测高度磁场加权系数：其中：ω7(Xp‑ξm,Yq‑ηn,Z0‑ζl)，ω8(Xp‑ξm,Yq‑ηn,Z0‑ζl)，ω9(Xp‑ξm,Yq‑ηn,Z0‑ζl)，ω10(Xp‑ξm,Yq‑ηn,Z0‑ζl)，ω11(Xp‑ξm,Yq‑ηn,Z0‑ζl)和ω12(Xp‑ξm,Yq‑ηn,Z0‑ζl)表示六个地上观测高度磁场加权系数，Z0表示观测高度，Xp，Yq分别表示观测高度水平观测点坐标，p＝1,2,…,Np，q＝1,2,…,Nq，Np和Nq分别表示观测高度水平观测点在x和y方向上的数量，水平观测点坐标间隔分别为Δx和Δy；

根据所述地上观测高度磁场加权系数，计算观测高度水平观测点强磁性体的异常磁场三个分量：

根据地上观测高度、水平观测点坐标、地下目标区域强磁性体的磁场，计算地上观测高度强磁性体的磁场梯度张量，包括：

根据所述地上观测高度、水平观测点坐标计算十个地上观测高度磁场梯度张量加权系数；

其中：ω13(Xp‑ξm,Yq‑ηn,Z0‑ζl)，ω14(Xp‑ξm,Yq‑ηn,Z0‑ζl)，ω15(Xp‑ξm,Yq‑ηn,Z0‑ζl)，ω16(Xp‑ξm,Yq‑ηn,Z0‑ζl)，ω17(Xp‑ξm,Yq‑ηn,Z0‑ζl)，ω18(Xp‑ξm,Yq‑ηn,Z0‑ζl)，ω19(Xp‑ξm,Yq‑ηn,Z0‑ζl)，ω20(Xp‑ξm,Yq‑ηn,Z0‑ζl)，ω21(Xp‑ξm,Yq‑ηn,Z0‑ζl)和ω22(Xp‑ξm,Yq‑ηn,Z0‑ζl)表示十个地上观测高度磁场梯度张量加权系数，Z0表示观测高度，Xp，Yq分别表示观测高度水平观测点坐标，p＝1,2,…,Np，q＝1,2,…,Nq，Np和Nq分别表示观测高度水平观测点在x和y方向上的数量，水平观测点坐标间隔分别为Δx和Δy；

根据所述观测高度磁场梯度张量加权系数，计算观测高度水平观测点强磁性体的磁场梯度张量六个分量：

若地上观测高度强磁性体的磁场、磁场梯度张量分别与仪器测量得到的地上观测高度强磁性体的实际磁场、实际磁场梯度张量相同，将强磁性体的磁化率分布数据作为所述强磁性体的实际磁化率分布数据以用于航空勘探所述强磁性体。

2.根据权利要求1所述的复杂强磁性体的航空勘探方法,其特征在于，步骤(j)中，给定的迭代收敛条件是：

ε0为期望的数值精度。

3.一种复杂强磁性体的航空磁法勘探系统，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1所述复杂强磁性体的航空勘探方法的步骤。