一种基于概率极值加权函数的多磁性目标三维重建方法技术

本发明专利技术涉及一种基于概率极值加权函数的多磁性目标三维重建方法，属于磁探测技术领域。该方法包括：S1：获取观测平面内的磁梯度张量数据，将反演空间划分为立方体网格单元，计算磁梯度张量正演的核矩阵；S2：利用相关成像法计算每个网格与测得磁梯度张量数据之间的归一化互相关系数C

【技术实现步骤摘要】
一种基于概率极值加权函数的多磁性目标三维重建方法


[0001]本专利技术属于磁探测
，涉及一种基于概率极值加权函数的多磁性目标三维重建方法。

技术介绍

[0002]磁异常三维重建技术是通过地下或水下磁性目标所引起的磁异常值，重建磁性目标的物性参数及其空间分布。该技术能有效地探测出磁性目标，广泛应用于反潜、未爆弹药探测、海底结构探测、地下铁磁性管道探测等，具有重要的实用价值。特别是针对水下目标的探测，重建的三维轮廓可以辅助目标识别，进一步提升探索海洋环境的能力。目前，对多个磁性目标的三维重建研究仍处于初级阶段。类似的磁反演技术主要研究在地质勘探中，包括形态反演和物理反演。形态学反演需要预先知道磁化率，这就导致了对未知磁性目标的重建有一定的局限性；而由于磁测数据的误差和先验信息的缺乏，传统的物理反演方法效果也并不理想。
[0003]因此，亟需一种能够准确实现多目标情况下的三维重建的方法。

技术实现思路

[0004]有鉴于此，本专利技术的目的在于提供一种基于概率极值加权函数的多磁性目标三维重建方法，针对现有空间分布不均的多个铁磁性目标三维重建方法的不足，解决存在多个目标时，由于磁测数据的误差和先验信息的缺乏情况下的三维重建方面存在准确度低、分辨率差、轮廓不清晰的问题。
[0005]为达到上述目的，本专利技术提供如下技术方案：
[0006]一种基于概率极值加权函数的多磁性目标三维重建方法，具体包括以下步骤：
[0007]S1：获取观测平面内的磁梯度张量数据，将反演空间划分为立方体网格单元，计算磁梯度张量正演的核矩阵；
[0008]S2：计算磁偶极子在每个网格单元上产生的理论磁梯度张量，并利用相关成像法计算其与测得磁梯度张量数据之间的归一化互相关系数C
ab,q
；
[0009]S3：采用26邻域法计算相关系数的极大值所在的网格位置，由于相关系数表征了磁性体存在的概率，也就得到了磁性体存在概率的极大值位置；
[0010]S4：将概率极大值位置作为先验信息，构造概率极值加权函数。按照每个网格到概率极大值位置的距离进行加权，靠近概率极大值位置的网格被赋予更低的权重，从而将反演结果约束在概率极大值位置附近的空间内；
[0011]S5：如果磁性体在x,y,z三个方向上分布不均匀，计算各个网格到概率极大值位置的距离时要根据先验信息予以三个方向不同的权重，即计算再加权距离；
[0012]S6：将概率极值加权函数加入反演目标函数求解，再采用可视化方法得到三维重建结果；
[0013]进一步，步骤S1中，假设立方体的边分别为a，b，c，中心坐标为(x0,y0,z0)，磁性立
方体在任意观测点上的磁感应强度三分量表达式为：
[0014][0015][0016][0017]其中(x,y,z)为测量平面内观测点的坐标，μ0为真空中的磁导率，J为立方体的总磁化强度，L0，M0，N0分别为总磁化强度的方向余弦，(ξ,η,ζ)为立方体场源点的坐标，r＝[(ξ
‑
x)2+(η
‑
y)2+(ζ
‑
z)2]1/2
，表示观测点与立方体内任意点间的距离；
[0018]对磁感应强度三分量表达式沿不同方向求导可得到对应的磁梯度张量分量表达式。
[0019]进一步，步骤S1中，由于位场的可叠加性，整个模型对任意观测点的磁梯度张量分量，等于组成该模型的网格单元在该点所产生磁梯度张量分量的线性叠加。
[0020]进一步，步骤S2中，归一化互相关系数C
ab,q
的计算公式为：
[0021][0022]其中，B
ab
(x
i
,y
i
,z
i
)为测区内任一点的实测磁梯度张量的某一分量，B
ab,q
(x
i
,y
i
,z
i
)为地下磁偶极子q在测区内任一点产生的磁梯度张量的某一个分量，a，b＝x，y，z，N为测区内观测点的总个数。
[0023]进一步，步骤S3中，互相关系数C
ab,q
表示该q点存在偶极子的可能性。互相关系数C
ab,q
表示该q点存在磁性体的可能性。C
ab,q
介于
‑
1到1之间，越接近1，存在磁性体的概率越大，因此相关系数的极大值位置反映了磁性体存在概率的极大值所在的位置。如果存在n个磁性目标，相应的可以得到n个这样概率极大值位置。
[0024]进一步，步骤S4中，利用磁性体概率极大值的位置坐标构造概率极值加权函数，其计算表达式为：
[0025][0026]其中，z0为观测平面相对于水平基准面的高度；z
j
为第j个立方体单元的中心理深；n表示磁性体存在概率的极大值的个数；r
1j
至r
nj
表示第j个网格单元中心到第1至第n个概率极大值位置的距离；r0是一个小的稳定常数，为了避免模型过于集中在概率极大值位置，r0的取值为一个网格单元的边长；
[0027]进一步，步骤S5中，第j个网格单元中心到概率极大值位置的再加权距离r
ij
(i＝1，
2，
…
，n)的计算表达式为：
[0028][0029]其中(x
i
，y
i
，z
i
)为第i个预估的磁性体概率极大值位置的坐标,(x
j
，y
j
，z
j
)为第j个网格单元中心的坐标，a
x
，a
y
，a
z
分别为三个方向上的先验比例权重，若无先验信息其默认值都为1。
[0030]进一步，步骤S6中，反演目标函数的表达式为：
[0031][0032]第一项为数据拟合函数，代表预测数据与实测数据的拟合程度。d代表观测数据，m为待求解的模型物性参数，G为由物性空间到数据空间的映射矩阵，也称为核矩阵；
[0033]第二项为模型目标函数，用来降低反演中的自由度以获得唯一解。其中，m0为参考模型，若无对应的先验信息时可设为0值向量；W
m
为模型加权矩阵，由不同类型的约束矩阵共同构成，如概率极值加权矩阵及粗糙度矩阵；β为正则化参数，作用是调节数据拟合目标函数和模型目标函数之间的权重。
[0034]理论上φ(m)的最小化，即在条件下的解m就是目标函数的最优解。
[0035][0036]对上式化简可得反演目标方程：
[0037](G
T
G+βW
mT
W
m
)m＝(G
T
d+βW
mT
W
m
m0)
[0038]进一步，步骤S6中，根据测量得到的磁梯度张量分量，利用共轭梯度法本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于概率极值加权函数的多磁性目标三维重建方法，其特征在于：该方法包括以下步骤：S1：获取观测平面内的磁梯度张量数据，将反演空间划分为立方体网格单元，计算磁梯度张量正演的核矩阵；S2：计算磁偶极子在每个网格单元上产生的理论磁梯度张量，并利用相关成像法计算其与测得磁梯度张量数据之间的归一化互相关系数C
ab,q
；S3：采用26邻域法计算相关系数的极大值所在的网格位置，相关系数表征磁性体存在的概率，得到磁性体存在概率的极大值位置；S4：将概率极大值位置作为先验信息，构造概率极值加权函数；按照每个网格到概率极大值位置的距离进行加权，靠近概率极大值位置的网格被赋予更低的权重，将反演结果约束在概率极大值位置附近的空间内；S5：如果磁性体在x,y,z三个方向上分布不均匀，计算各个网格到概率极大值位置的距离时要根据先验信息予以三个方向不同的权重，即计算再加权距离；S6：将概率极值加权函数加入反演目标函数求解，再采用可视化方法得到三维重建结果。2.根据权利要求1所述的一种基于概率极值加权函数的多磁性目标三维重建方法，其特征在于：所述S1中，设立方体的边分别为a，b，c，中心坐标为(x0,y0,z0)，磁性立方体在任意观测点上的磁感应强度三分量表达式为：意观测点上的磁感应强度三分量表达式为：意观测点上的磁感应强度三分量表达式为：其中(x,y,z)为测量平面内观测点的坐标，μ0为真空中的磁导率，J为立方体的总磁化强度，L0，M0，N0分别为总磁化强度的方向余弦，(ξ,η,ζ)为立方体场源点的坐标，r＝[(ξ
‑
x)2+(η
‑
y)2+(ζ
‑
z)2]
1/2
，表示观测点与立方体内任意点间的距离；对磁感应强度三分量表达式沿不同方向求导得到对应的磁梯度张量分量表达式。3.根据权利要求2所述的一种基于概率极值加权函数的多磁性目标三维重建方法，其特征在于：所述S1中，根据位场的可叠加性，整个模型对任意观测点的磁梯度张量分量，等于组成该模型的网格单元在该点所产生磁梯度张量分量的线性叠加。4.根据权利要求3所述的一种基于概率极值加权函数的多磁性目标三维重建方法，其特征在于：所述S2中，归一化互相关系数C
ab,q
的计算公式为：
其中，B
ab
(x
i
,y
i
,z
i
)为测区内任一点的实测磁梯度张量的某一分量，B
ab,q
(x
i
,y
i
,z
i
)为地下磁偶极子q在测区内任一点产生的磁梯度张量的某一个分量，a，b＝x，y，N为测区内观测点的总个数。5.根据权利要求4所述的一种基于概率极值加权函数的多磁性目标三维重建方法，其特征在于：所述S3中，互相关系数C
ab,q
表示该q点存在磁性体的可能性；C
ab,q
介于
‑
1到1之间，越接近1，存在磁性体的概率越大，相关系数的极大值位置反...

【专利技术属性】
技术研发人员：邱景，兰天，黄双龙，
申请(专利权)人：重庆大学，
类型：发明
国别省市：