【技术实现步骤摘要】
本专利技术属于磁性目标探测领域,尤其涉及一种在地磁测量及地磁反演应用中利用小子域识别算法确定磁性目标位置的方法。
技术介绍
地磁场一般可分为地磁正常场和地磁异常场。在实际测量过程中,正常场与异常场的划分往往是相对的。实际上测得的地球磁场强度和理论磁场强度是有区别的,这种区别称为地磁异常。它主要是由磁性物体受地磁场磁化而产生的附加磁场。所以,通过对磁异常的研究可以帮助人们有效的寻找磁性物体。根据磁异常信息,通过一定的数学解释手段,判断磁性目标的空间位置、大小、形状等特征参数,可实现目标监控和识别。Euler反褶积方法是以欧拉齐次方程为基础,推断磁性目标中心的空间位置。然而反演场源水平位置和深度时,Euler齐次方程是一个灵敏度极高的方程式,式中异常导数的计算精度等稍有偏差,便会在反演中产生较强的干扰。小子域识别算法是基于滑动平均法原理进行设计的算法,不仅保留了压缩异常区域的功能,能以较高的分辨率保留异常中的区域特征,对具有梯级带特征的断裂构造位置的确定具有较高的分辨率,同时可以有效的识别窗口是否处于异常区域内,这样就可以在压缩异常区域过程中,将异常区域的中心位置显示出来。举例说明算法的实现过程:当磁场梯度带从左向右增大时,小子域的滤波窗口从左向右滑动,各小子域中以平缓系数最小域内的平均值作为窗口中心位置的值。当窗口右侧的子域处于磁异常的梯度带上,这些子域的均方差较大;而左测的子域处于磁异常的平缓区时,这些子域的均方 ...
【技术保护点】
一种利用小子域识别算法确定磁性目标位置的方法,其特征在于:包括以下步骤,步骤一:利用磁力仪,获取区域的局部地磁图数据;步骤二:构建5X5数据点的窗口,将窗口划分成9个小子域,设定初始阈值K=0.1;步骤三:分别计算9个子域里的磁异常方差σi,将9个子域中的方差最小对应的小子域的磁异常均值Δgmin作为窗口中心的值;其中,Ti(p,q)为第i个子域内第(p,q)个点上的异常值,σi为第i个子域的方差值;步骤四:计算整个窗口的方差Sσ:Sσ=1p×qΣp=1,q=15(T(p,q)-Δgmin)2;]]>步骤五:比较整个窗口的方差Sσ和阈值K的大小,当Sσ<K时,窗口中心位置的识别输出值设置为0;当Sσ>K时,窗口中心位置的识别输出值设置为1;步骤六:窗口滑动到下一点,重复步骤二~步骤五,直至完成全部的区域计算;将整个区域中最大的整个窗口的方差Sσ赋予阈值K;步骤七:重复步骤二~步骤六,判断整个窗口的方差Sσ是否和阈值K相等,如果相等,则进行步骤八;如果不相等,重复步骤二~步骤五,直至整个窗口的方差Sσ和阈值K相等;步骤八:得到地磁图中识别标 ...
【技术特征摘要】
1.一种利用小子域识别算法确定磁性目标位置的方法,其特征在于:包括以下步骤,
步骤一:利用磁力仪,获取区域的局部地磁图数据;
步骤二:构建5X5数据点的窗口,将窗口划分成9个小子域,设定初始阈值K=0.1;
步骤三:分别计算9个子域里的磁异常方差σi,将9个子域中的方差最小对应的小子域的
磁异常均值Δgmin作为窗口中心的值;
其中,Ti(p,q)为第i个子域内第(p,q)个点上的异常值,σi为
第i个子域的方差值;
步骤四:计算整个窗口的方差Sσ:
Sσ=1p×qΣp=1,q=15(T(...
【专利技术属性】
技术研发人员:康崇,樊黎明,张晓峻,王明,郑权,周健,
申请(专利权)人:哈尔滨工程大学,
类型:发明
国别省市:黑龙江;23
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。