【技术实现步骤摘要】
一种基于PDE的磁测数据等效源上延拓与下延拓方法
本专利技术涉及地球物理勘测
，尤其涉及一种基于PDE的磁测数据等效源上延拓与下延拓方法。
技术介绍
利用地球物理磁法(如航磁测量、地面磁法测量)进行探测时，测量分布在起伏观测曲面上，由于测量位置与磁异常体的距离不同，磁测数据所反映的主要异常信息内容也有所差异。在实际应用中，需要对数据的测量位置进行转换，以满足磁测数据处理的解释需求，例如，利用磁测数据对深部异常结构进行研究时，往往需要更高平面的观测数据，上延拓即是将观测面上的数据延拓到更高的观测面上，用于消除离观测面较近的异常体对磁场的影响，从而突出深层的磁异常；而下延拓即是将观测面上的数据延拓到更低的观测面上，用于消除离观测面较远的异常体对磁场的影响，从而突出浅层的磁异常。将分布在起伏观测曲面上不规则观测位置的磁测数据点，直接延拓到另一个目标曲面上，比如航测起伏观测曲面下延拓至地面、地面数据上延拓至航磁起伏观测曲面等，相比延拓至平面或利用中间平面实现曲面间的数据延拓计算，具有更高的计算精度。现有文献1“Damp...
【技术保护点】
1.一种基于PDE的磁测数据等效源上延拓与下延拓方法，其特征在于，包括以下步骤：/nS1、获取起伏观测曲面上的磁场数据d
【技术特征摘要】
1.一种基于PDE的磁测数据等效源上延拓与下延拓方法，其特征在于，包括以下步骤：
S1、获取起伏观测曲面上的磁场数据d0，并根据磁场数据所在区域的地形高度信息，建立地形起伏曲面；
S2、根据起伏观测曲面的高程信息以及设定的反演最大深度，确定网格剖分的空间范围，并根据地形起伏曲面，对所述空间范围进行连续的结构化非均匀网格剖分，进一步确定等效源反演网格空间；
S3、根据地磁场的磁倾角、磁偏角、磁感应强度，基于所述等效源反演网格空间对磁场数据d0进行带深度规整化因子、正值约束项以及规整化项的PDE三维反演计算，得到磁异常体的多层等效源模型；
S4、根据延拓后的起伏观测曲面的位置，利用步骤S3得到的多层等效源模型进行基于PDE的磁场正演计算，得到磁异常体产生的上延拓和/或下延拓后的磁场数据d。
2.根据权利要求1所述的基于PDE的磁测数据等效源上延拓与下延拓方法，其特征在于，所述多层等效源模型的模型深度面的层数大于3层。
3.根据权利要求1所述的基于PDE的磁测数据等效源上延拓与下延拓方法，其特征在于，步骤S3中，所述PDE三维反演计算的目标函数为：
其中，
Bs＝F(B0,m)
m≥0
式中，φ表示优化目标，m表示输出的多层等效源模型的磁化率矩阵；F(·)表示多层等效源模型的正演计算，Bs表示正演操作得到的磁异常分量数据，T(·)表示磁异常分量数据到磁异常总场数据的转换函数；Qx、Qy、Qz分别表示北向、东向及垂向上的插值函数；B0表示地磁场强度矢量，由地磁场的磁倾角、磁偏角以及磁感应强度确定，B0x、B0y、B0z分别表示其北向、东向及垂向分量；β表示预设的规整...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。