1.一种磁测数据多层等效源上延拓与下延拓方法，其特征在于，包括以下步骤：

S1、输入起伏观测曲面上的磁场数据d0，并根据起伏观测曲面所在区域的地形高度信息，建立地形起伏曲面；

S2、根据起伏观测曲面以及设定的反演最大深度，确定网格剖分的空间范围，并根据地形起伏曲面对所述空间范围进行连续的结构化非均匀网格剖分，进一步确定反演模型求解空间；

S3、根据地磁场的磁倾角、磁偏角、磁感应强度，基于反演模型求解空间对磁场数据d0进行带深度规整化因子、正值约束项以及规整化项的积分方程三维反演计算，得到磁异常体的多层等效源模型；

S4、根据延拓后的起伏观测曲面的位置，利用步骤S3得到的多层等效源模型进行基于积分方程的磁场正演计算，得到磁异常体产生的上延拓和/或下延拓后的磁场数据d。

2.根据权利要求1所述的磁测数据多层等效源上延拓与下延拓方法，其特征在于，所述多层等效源模型的模型深度面的数目大于3层。

3.根据权利要求1所述的磁测数据多层等效源上延拓与下延拓方法，其特征在于，步骤S3中，积分方程三维反演计算的目标函数为：其中，m≥0

式中，φ表示优化目标，m表示输出的多层等效源模型的磁化率矩阵；F(·)表示基于积分方程的磁场正演操作，G表示灵敏度矩阵，与起伏观测曲面的位置以及地磁场的磁倾角、磁偏角、磁感应强度相关；β表示预设的规整化因子，mref表示参考模型的磁化率矩阵，Wr表示深度规整化因子。

4.根据权利要求3所述的磁测数据多层等效源上延拓与下延拓方法，其特征在于，所述深度规整化因子为：其中，z表示等效源到地形起伏曲面的距离，z0表示起伏观测曲面的高度，r表示深度系数。

5.根据权利要求1所述的磁测数据多层等效源上延拓与下延拓方法，其特征在于，步骤S2中，所述网格剖分的空间范围包括上顶面和下底面，其中，所述上顶面为地形起伏曲面的最大高度所确定的平面，所述下底面为设定的反演最大深度所确定的平面。

6.根据权利要求5所述的磁测数据多层等效源上延拓与下延拓方法，其特征在于，步骤S2中，根据地形起伏曲面的最低点对网格剖分的空间范围进行划分，其中，对所述最低点以上的空间范围进行均匀网格剖分得到精细网格，对所述最低点以下的空间范围进行非均匀网格剖分得到扩展网格；进一步，所述地形起伏曲面至下底面之间的空间范围构成反演模型求解空间。

7.根据权利要求6所述的磁测数据多层等效源上延拓与下延拓方法，其特征在于，所述扩展网格的垂直边以精细网格的垂直边的α1倍速度增长，且设定其最大增速为α2，其中，α2＞α1＞1。

8.根据权利要求3所述的磁测数据多层等效源上延拓与下延拓方法，其特征在于，步骤S4中，根据延拓后的起伏观测曲面的位置以及地磁场的磁倾角、磁偏角、磁感应强度，确定灵敏度矩阵G′，进一步得到上延拓和/或下延拓后的磁场数据d＝F(m,G′)，其中，m为步骤S3输出的磁化率矩阵。