【技术实现步骤摘要】
一种基于PDE的磁测数据等效源化极与类型转换方法
本专利技术涉及地球物理勘测
，尤其涉及一种基于PDE的磁测数据等效源化极与类型转换方法。
技术介绍
在地磁探测中，人们通常对磁场的总场异常进行探测，但在实际的数据解释中，又往往需要将这些数据转换成所需要的参数和类型，比如不同磁化方向的换算(化极数据)、磁场分量、张量数据的换算等，磁场数据类型转换的主要挑战在于航磁(地面磁测)的起伏测量曲面以及数据的不规则测量位置(非网格化测量)，传统的数据类型转换以及化极计算无法处理这些数据，而采用传统单层、双层、三层等效源的数据转换方法，处理结果精度较低，无法满足应用需求。现有文献1“Dampney,C.N.G.THEEQUIVALENTSOURCETECHNIQUE[J].geophysics,1969,34(1):39.”首次提出了等效源方法，利用虚拟场源模拟实测异常，可用于位场数据的空间延拓(包括曲面延拓)、梯度计算以及分量转换等；选择单层等效源且将其布置于近地表是等效源方法的主要特点，比如,文献2“LiY.Reduct...
【技术保护点】
1.一种基于PDE的磁测数据等效源化极与类型转换方法，其特征在于，包括以下步骤：/nS1、获取起伏观测曲面上的磁场数据d
【技术特征摘要】
1.一种基于PDE的磁测数据等效源化极与类型转换方法，其特征在于，包括以下步骤：
S1、获取起伏观测曲面上的磁场数据d0，并根据磁场数据所在区域的地形高度信息，建立地形起伏曲面；
S2、根据起伏观测曲面的高程信息以及设定的反演最大深度，确定网格剖分的空间范围，并根据地形起伏曲面，对所述空间范围进行连续的结构化非均匀网格剖分，进一步确定等效源反演网格空间；
S3、根据地磁场的磁倾角、磁偏角和磁感应强度，基于所述等效源反演网格空间对磁场数据d0进行带深度规整化因子、正值约束项以及规整化项的PDE三维反演计算，得到磁异常体的多层等效源模型；
S4、给定不同的地磁极参数，利用多层等效源模型进行非线性PDE的磁场三维正演计算，得到磁场数据d0的化极数据Bs′、磁异常体产生的磁异常总场数据、磁异常分量数据Bs以及磁梯度张量数据中的一种或多种。
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分别表示其北向、东向及垂向分量；β表示预设的规整化因子；mref表示参考等效源模型的磁化率矩阵，Wr表示深度规整化因子。
4.根据权利要求3所述的基于PDE的磁测数据等效源化极与类型转换方法，其特征在于，所述深度规整化因子为：
其中，z表示等...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。