CN114966866A - 一种基于矩形阵列的水中运动磁性目标探测定位系统 - Google Patents

一种基于矩形阵列的水中运动磁性目标探测定位系统 Download PDF

Info

Publication number
CN114966866A
CN114966866A CN202111654555.4A CN202111654555A CN114966866A CN 114966866 A CN114966866 A CN 114966866A CN 202111654555 A CN202111654555 A CN 202111654555A CN 114966866 A CN114966866 A CN 114966866A
Authority
CN
China
Prior art keywords
magnetic
initial
detection
target
moment
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111654555.4A
Other languages
English (en)
Inventor
王作帅
左超
徐林
杨文铁
陈盟
王建勋
陈志伟
肖涵琛
耿攀
周诗颖
魏华
陈涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
719th Research Institute of CSIC
Original Assignee
719th Research Institute of CSIC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 719th Research Institute of CSIC filed Critical 719th Research Institute of CSIC
Priority to CN202111654555.4A priority Critical patent/CN114966866A/zh
Publication of CN114966866A publication Critical patent/CN114966866A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
    • G01V3/081—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices the magnetic field is produced by the objects or geological structures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/04—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
    • G01C21/08—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means involving use of the magnetic field of the earth
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Abstract

本申请公开了一种基于矩形阵列的水中运动磁性目标探测定位系统，该系统包括：探测矩阵，磁通计算单元以及初始探测单元；探测矩阵中设置有多个感应式磁场传感器，感应式磁场传感器用于获取探测范围内因水中磁性目标移动产生的感应电动势；磁通计算单元用于当判定感应电动势符合预设条件时，计算探测矩阵中阵列单元所在位置处的磁通变化量；初始探测单元用于根据磁通变化量，通过导数运算，计算磁场时变量，并结合磁通变化量，计算水中磁性目标的初始位置和初始磁矩，并生成初始定位信息，其中，初始定位信息包括水中磁性目标的位置和类型，初始磁矩用于确定类型。通过本申请中的技术方案，克服了探测定位过程中易受海洋环境磁场干扰的缺点。

Description

一种基于矩形阵列的水中运动磁性目标探测定位系统
技术领域
本申请涉及水中目标探测的技术领域，具体而言，涉及一种基于矩形 阵列的水中运动磁性目标探测定位系统。
背景技术
为实现对我国海湾或海岸军事要塞、重要港口、重要水域的安全保 卫，亟需建设水下探测和定位系统，实现水下安全的常态化警戒。
在对水中目标进行探测时，可通过水中目标的声场、电磁场和光学等 特征进行探测和定位，实现水下安全警戒。对于大多数水中目标来说，因 其由金属制造，在空间中会产生静态的磁场分布，利用其磁场特征信号可 以实现探测和定位。
目前磁性目标定位算法均为以目标在空间产生的静态磁场分布为对 象，通过设置静态磁场测量传感器，实现对水下目标的探测和定位。但由 于地磁场通常为几万纳特，而水中磁性目标在远距离处产生的静态场仅为 纳特级甚至更低，特别是在考虑到海洋环境的纳特级电磁干扰时，基于静 态磁场的水中磁性目标探测方案存在信号提取困难、信噪比低等弊端，限 制了水中磁性目标定位距离，且此类探测定位方法需要布置大量的磁场传感器。
另外，还可以通过设置感应式磁场测量传感器，实现对指定区域内的 水中磁性目标进行探测，但这种方式仅能通过测量环境中的磁通变化量， 以反映是否有磁性目标经过，但并不具备对水中磁性目标的定位跟踪功 能。
发明内容
本申请的目的在于：为了解决现有水中磁性目标探测和定位过程中存 在的至少一个技术问题，克服探测定位过程中易受海洋环境磁场干扰的缺 点，减少磁场传感器数量，提高对水中运动目标探测定位的精度。
本申请的技术方案是：提供了一种基于矩形阵列的水中运动磁性目标 探测定位系统，该探测定位系统包括：探测矩阵，磁通计算单元以及初始 探测单元；探测矩阵中设置有多个感应式磁场传感器，感应式磁场传感器 用于获取探测范围内因水中磁性目标移动产生的感应电动势；磁通计算单 元用于当判定感应电动势符合预设条件时，计算探测矩阵中阵列单元所在 位置处的磁通变化量；初始探测单元用于根据磁通变化量，通过导数运算，计算磁场时变量，并结合磁通变化量，计算水中磁性目标的初始位置 和初始磁矩，并生成初始定位信息，其中，初始定位信息包括水中磁性目 标的位置和类型，初始磁矩用于确定类型。
上述任一项技术方案中，进一步地，探测定位系统还包括：寻优单 元；寻优单元用于根据初始位置和初始磁矩，确定水中磁性目标的磁场定 位初始解集，并根据磁场定位初始解集和预设优化目标函数，计算水中磁 性目标的最优位置和最优磁矩，并生成最优定位信息，其中，最优定位信 息包括水中磁性目标的最优位置和最优类型，最优磁矩用于确定最优类 型。
上述任一项技术方案中，进一步地，预设优化目标函数的计算公式 为：
Obj＝∑[(Uxc-Uxm)2+(Uyc-Uym)2+(Uzc-Uzm)2]
式中，Uxc为感应式磁场传感器中磁通x分量Bx对应的感应线圈棒的 电压计算值，Uxm为磁通x分量Bx对应的电压测量值，Uyc为感应式磁场传 感器中磁通y分量By对应的感应线圈棒的电压计算值，Uym为磁通y分量 By对应的电压测量值，Uzc为感应式磁场传感器中磁通z分量Bz对应的感 应线圈棒的电压的计算值，Uzm为磁通z分量Bz对应的电压测量值。
上述任一项技术方案中，进一步地，寻优单元还用于：根据上一采样 时刻的最优位置和最优磁矩，确定下一采样时刻的磁场定位初始解集，并 根据预设优化目标函数，计算下一采样时刻的最优位置和最优磁矩；探测 定位系统还包括：电动势计算单元以及判断单元；电动势计算单元用于根 据下一采样时刻的最优位置和最优磁矩，计算下一采样时刻对应的感应式 磁场传感器的感应电动势计算值以及误差值，其中，误差值为下一采样时刻感应电动势计算值与感应电动势测量值之间的误差；判断单元用于当判 定至少两个感应式磁场传感器的误差值小于误差阈值，且至少四个以上感 应式磁场传感器的感应电动势大于探测阈值时，根据计算出的下一采样时 刻的最优位置和最优磁矩，生成水中磁性目标下一采样时刻的最优定位信 息。
上述任一项技术方案中，进一步地，阵列单元由探测矩阵中的异常传 感器组成，异常传感器为获取到的感应电动势大于探测阈值的感应式磁场 传感器，预设条件为异常传感器的数量大于或等于数量阈值。
上述任一项技术方案中，进一步地，数量阈值为四个，预设条件还包 括：多个异常传感器中存在至少四个异常传感器可以组成矩形阵列，将组 成的矩形阵列记作阵列单元。
上述任一项技术方案中，进一步地，初始探测单元中包括：导数计算 单元，张量矩阵计算单元以及磁矩计算单元；导数计算单元用于根据水中 磁性目标的移动速度和磁通变化量，计算阵列单元所处位置几何中心处的 磁通关于时间的导数；张量矩阵计算单元用于根据导数，构建时变磁梯度 张量矩阵；磁矩计算单元用于根据时变磁梯度张量矩阵，计算初始位置和 初始磁矩，其中，初始磁矩的计算公式为：
Figure BDA0003447910560000041
Figure BDA0003447910560000042
Figure BDA0003447910560000043
式中，mx0、my0、mz0为水中磁性目标在三维坐标系中的三分量磁矩， R0为水中磁性目标距离探测原点的距离，探测原点定义为矩形阵列的几何 中心，μ0为真空磁导率，x0、y0、z0为水中磁性目标相对探测原点的坐 标，G为时变磁梯度张量矩阵。
上述任一项技术方案中，进一步地，感应式磁场传感器由三个相互正 交的感应线圈棒组成。
本申请的有益效果是：
本申请中的技术方案，利用感应式磁场传感器组成探测矩阵，再通过 筛选感应电动势符合预设条件的感应式磁场传感器组成阵列单元，以对水 中磁性目标进行探测和定位，算法适应性好，同时能满足精度要求，基于 该探测定位系统搭建的系统结构精简，磁场传感器数量较少，效费比高。
本申请在进行探测和定位计算时，引入导数运算，计算由水中磁性目 标的移动速度和磁通变化量构成的时变磁梯度张量矩阵，进而得到水中磁 性目标的初始位置和初始磁矩，以实现水中磁性运动目标引起的磁通量变 化测量其感应量进而进行定位，避免了海洋环境波动的干扰，提高了对水 中运动目标探测定位的精度。
附图说明
本申请的上述和/或附加方面的优点在结合下面附图对实施例的描述中 将变得明显和容易理解，其中：
图1是根据本申请的一个实施例的基于矩形阵列的水中运动磁性目标 探测定位系统的示意框图；
图2(a)是根据本申请的一个实施例的感应式磁场传感器的示意图；
图2(b)是根据本申请的一个实施例的探测矩阵的示意图；
图3是根据本申请的一个实施例的采用差分进化算法进行寻优的示意 图；
图4是根据本申请的一个实施例的精度对比试验结果的示意图。
具体实施方式
为了能够更清楚地理解本申请的上述目的、特征和优点，下面结合附 图和具体实施方式对本申请进行进一步的详细描述。需要说明的是，在不 冲突的情况下，本申请的实施例及实施例中的特征可以相互结合。
在下面的描述中，阐述了很多具体细节以便于充分理解本申请，但 是，本申请还可以采用其他不同于在此描述的其他方式来实施，因此，本 申请的保护范围并不受下面公开的具体实施例的限制。
如图1所示，本实施例提供了一种基于矩形阵列的水中运动磁性目标 探测定位系统，该探测定位系统包括：探测矩阵，磁通计算单元以及初始 探测单元；其中，探测矩阵中设置有多个感应式磁场传感器，感应式磁场 传感器用于获取探测范围内因水中磁性目标移动产生的感应电动势；磁通 计算单元用于当判定感应电动势符合预设条件时，计算探测矩阵中阵列单 元所在位置处的磁通变化量，感应式磁场传感器由三个相互正交的感应线圈棒组成。
进一步的，阵列单元由探测矩阵中的异常传感器组成，异常传感器为 获取到的感应电动势大于探测阈值的感应式磁场传感器，预设条件为异常 传感器的数量大于或等于数量阈值。
优选的，数量阈值为四个，预设条件还包括：多个异常传感器中存在 至少四个异常传感器可以组成矩形阵列，将组成的矩形阵列记作阵列单 元。
具体的，以探测矩阵所在区域的几何中心为原点，以垂直向下的方向为 z轴正方向，基于右手定则建立三维坐标系。
如图2所示，所有的感应式磁场传感器由三个相互正交的感应线圈棒 组成，三个相互正交的感应线圈棒与三维坐标系的x、y、z轴平行，利用感 应式磁场传感器获取探测范围内的感应电动势，作为实时采集数据，在任意 采样时刻，判断是否符合预设条件，预设条件如下：
1、存在四个以上感应式磁场传感器探测到的感应电动势大于探测阈值， 将此类感应式磁场传感器记作异常传感器，若存在四个以上的异常传感器， 则说明有水中磁性目标出现，进行下一步判断；否则，重复执行上述过程。
2、根据异常传感器的数目，判断是否存在四个异常传感器，能够组成 一个矩形阵列，若存在，则将该矩形阵列记作阵列单元；否则，重复执行上 述过程。
需要说明的是，符合预设条件的感应式磁场传感器可形成多个矩形阵 列单元，此时，选取阵列单元中感应电动势最大的矩形阵列单元作为计算依 据。
当判定感应电动势符合预设条件时，根据矩形阵列单元中各个感应式 磁场传感器测量得到感应电动势，基于电磁感应定律，计算各感应式磁场传 感器所处位置的磁通变化量。
本实施例中采用电磁感应定律获得磁通变化量的过程如下：
根据磁感应线圈探测变化磁场的原理，当探测范围内存在移动的水中 磁性目标时，将引起各感应式磁场传感器中磁通量的变化，根据磁通量变化 引起的感应电动势公式和磁通量计算公式，可得：
Figure BDA0003447910560000071
φ＝BS
Figure BDA0003447910560000072
其中，ε为感应电动势，n为线圈匝数，Sn为感应线圈垂直磁场B的切 面面积，nSn为感应线圈的线圈常数，B是磁感应强度，Bt为磁通变化量。
因此，若获得磁感应线圈的感应电动势，即可获得由水中磁性目标体 移动引起的磁通变化量
Figure BDA0003447910560000073
所以有
Figure BDA0003447910560000074
即：
Figure BDA0003447910560000075
式中，Bt1、Bt2为相邻两个采样时刻t1与t2的磁通量。
记构成矩形阵列单元的4个感应式磁场传感器的编号分别为1、2、3、 4，以上述探测矩阵的第一个三维坐标系为基础，以该4个感应式磁场传感 器几何中心为矩阵原点(0,0,0)，建立第二个三维坐标系，根据4个感应式磁 场传感器在第一个三维坐标系中的位置，可计算出其在第二个三维坐标系 中的坐标，此时，4个感应式磁场传感器的空间坐标分别为：(x1,y1,z1)、 (x2,y2,z2)、(x3,y3,z3)、(x4,y4,z4)，具体计算过程不再赘述。
需要说明的是，以下计算以第二个三维坐标系，即矩形阵列单元的三维 坐标系为基准。
基于上述磁感应线圈的感应电动势原理，可获得4个感应式磁场传感 器的磁通量变化量依次为：
Figure BDA0003447910560000081
假设在矩形阵列区域内磁场梯度为均匀梯度，可采用平均值方式，获得 矩形阵列单元几何中心处的磁通变化量
Figure BDA0003447910560000082
将该 磁通变化量作为探测矩阵中阵列单元所在位置处的磁通变化量。
本实施例中，初始探测单元用于根据磁通变化量，通过导数运算，计 算磁场时变量，并结合磁通变化量，计算水中磁性目标的初始位置和初始 磁矩，并生成初始定位信息，其中，初始定位信息包括水中磁性目标的位 置和类型，初始磁矩用于确定类型。
具体的，针对矩形阵列单元几何中心处的磁通变化量
Figure BDA0003447910560000083
在磁梯度法 中对磁场求导，构建时变磁场梯度法，计算磁场时变量，与感应式磁场传感 器的磁通量变化量进行配合，计算得到水中磁性目标的初始位置和初始磁 矩，以便估算水中磁性目标的大小范围和所处位置区间。
其中，时变磁场梯度法及水中磁性目标的初始位置、初始磁矩计算细化 如下：
假设水中磁性目标的移动速度为V＝(Vx,Vy,Vz)，其移动方向的方位角 为ξ(绕z轴旋转)，倾角为η(绕x轴旋转)，进行速度分解，可得到以下关 系：
Vx＝Vsinξcosη
Vy＝Vcosξcosη
Vx＝Vsinξ
其中，V为移动速度V＝(Vx,Vy,Vz)的模。
由水中磁性目标移动引起的位置坐标变化可表示为：
xs＝xs0+Vxt
ys＝ys0+Vyt
zs＝zs0+Vzt
其中，(xs0,ys0,zs0)为0时刻水中磁性目标的位置坐标，(xs,ys,zs)为t 时刻水中磁性目标的位置坐标。
由以上公式可知，t时刻水中磁性目标位置坐标(xs,ys,zs)为时间t的函 数，所以有
Figure BDA0003447910560000091
Figure BDA0003447910560000092
Figure BDA0003447910560000093
进一步的，该初始探测单元中包括：导数计算单元，张量矩阵计算单 元以及磁矩计算单元；其中，导数计算单元用于根据水中磁性目标的移动 速度和磁通变化量，计算阵列单元所处位置几何中心处的磁通关于时间的 导数；张量矩阵计算单元用于根据导数，构建时变磁梯度张量矩阵；磁矩 计算单元用于根据时变磁梯度张量矩阵，计算初始位置和初始磁矩。
具体的，在磁梯度法中对磁场求导，构建时变磁场梯度法，计算磁场时 变量，因此，在三维坐标系中，对矩形阵列单元所处位置几何中心处的磁通 变化量
Figure BDA0003447910560000101
进行分解，分解为磁通x分量Bx、磁通y分量By以及磁通z分量 Bz，然后计算上述三个分量关于时间的导数：
Figure BDA0003447910560000102
Figure BDA0003447910560000103
Figure BDA0003447910560000104
其中，
Figure BDA0003447910560000105
为磁通x分量关于x方向的导数，在矩形阵列中，可表示为 三维坐标系中x方向两个磁感应线圈的测量值之差(磁通量变化)与距离 之差的比值，其余分量定义类似，不再赘述。
因此，将磁通关于时间的导数简化可得：
Figure BDA0003447910560000106
Figure BDA0003447910560000107
Figure BDA0003447910560000108
其中，
Figure BDA0003447910560000109
Figure BDA0003447910560000111
Figure BDA0003447910560000112
Figure BDA0003447910560000113
Figure BDA0003447910560000114
Bzz＝-Bxx-Byy
式中，μ0为真空磁导率，为常数，R0为t时刻水中磁性目标距离探测原 点的距离，
Figure BDA0003447910560000115
Bxx为磁通x分量Bx在x方向的导数， Bxy为磁通x分量Bx在y方向的导数，Bxz为磁通x分量Bx在z方向的导数， 其余参数定义类似，不在赘述。
本实施例中，构建出的时变磁梯度张量矩阵为：
Figure BDA0003447910560000116
具体的，根据磁梯度张量定位原理，已知目标磁梯度张量矩阵以及测量 点位置处的磁场值，可反演出目标位置，类似的，结合目标时变磁梯度张量 矩阵和测量点处磁场变化率，可计算得到目标位置信息，因此，水中磁性目 标初始位置的定位公式为：
Figure BDA0003447910560000121
R0为水中磁性目标距离探测原点的距离，探测原点定义为矩形阵列单 元的几何中心，获得的目标位置初值记作x0，y0，z0。再通过坐标变换的方 式，计算水中磁性目标距离R0相对于探测矩阵所在区域的几何中心的位置， 即水中磁性目标的初始位置。
根据上式确定了水中磁性目标的初始位置(x0，y0，z0)后，再根据下式 推算水中磁性目标的初始磁矩，通过初始磁矩即可对目标的状态、类型和尺 度等信息作出初步判断。
Figure BDA0003447910560000122
其中，mx0、my0、mz0为水中磁性目标在三维坐标系中的三分量磁矩。
具体的，一般由钢铁材料制成的水中磁性目标，其磁矩与目标主尺度强 相关，小型蛙人装备的磁矩约为100A·m2以内，中型无人潜航器目标磁矩 约为3～6万A·m2，大型有人目标磁矩约为40～80万A·m2。因此，可以 根据初始磁矩判断水中磁性目标的类型和尺寸。
本实施例中，该探测定位系统还包括：寻优单元；寻优单元用于根据初 始位置和初始磁矩，确定水中磁性目标的磁场定位初始解集，并根据磁场定 位初始解集和预设优化目标函数，计算水中磁性目标的最优位置和最优磁 矩，并生成最优定位信息，其中，最优定位信息包括水中磁性目标的最优位 置和最优类型，最优磁矩用于确定最优类型。
具体的，获得水中磁性目标的初始位置(x0，y0，z0)、初始磁矩(mx0， my0，mz0)后，以该值为中心，取左右各50％为参数区间，以x0为例，x0的 参数区间即为x∈[0.5x0,1.5x0]，估算水中磁性目标大小范围和位置区间。
根据估算得到的水中磁性目标的大小范围和所处的位置区间，按照均 匀概率原则，随机生成一组水中磁性目标的位置和磁矩值，与计算得到的一 系列水中磁性目标的初始位置和初始磁矩，记作水中磁性目标的磁场定位 初始解集。
如图3所示，将该磁场定位初始解集作为初代种群，采用差分进化算 法进行寻优，得到当前采样时刻的水中磁性目标的位置和磁矩参数，并正演 计算传感器组所处位置处的感应电动势，与传感器组测量值进行对比，优化 目标函数为感应线圈棒电压计算值与测量值之间误差最小。
本实施例中，设定预设优化目标函数的计算公式为：
Obj＝Σ[(Uxc-Uxm)2+(Uyc-Uym)2+(Uzc-Uzm)2]
其中，Uxc为感应式磁场传感器中磁通x分量Bx对应的感应线圈棒的电 压计算值，该计算值由水中磁性目标的磁场定位初始解集通过本实施例中 初始位置和初始磁矩的反向运算计算而得，Uxm为磁通x分量Bx对应的感应 线圈棒的电压的测量值，Uyc为感应式磁场传感器中磁通y分量By对应的感 应线圈棒的电压计算值，Uym为磁通y分量By对应的感应线圈棒的电压测 量值，Uzc为感应式磁场传感器中磁通z分量Bz对应的感应线圈棒的电压计算值，Uzm为磁通z分量Bz对应的感应线圈棒的电压测量值。
正演计算主要依据优化得到的最优位置(x，y，z)计算最优磁矩(mx、 my、mz)，计算公式如下：
Figure BDA0003447910560000141
式中，R由最优位置(x，y，z)确定。
在上述实施例的基础上，该寻优单元还用于：根据上一采样时刻的最 优位置和最优磁矩，确定下一采样时刻的磁场定位初始解集，并根据预设 优化目标函数，计算下一采样时刻的最优位置和最优磁矩；探测定位系统 还包括：电动势计算单元以及判断单元；电动势计算单元用于根据下一采 样时刻的最优位置和最优磁矩，计算下一采样时刻对应的感应式磁场传感 器的感应电动势计算值以及误差值，其中，误差值为下一采样时刻感应电动势计算值与感应电动势测量值之间的误差；判断单元用于当判定至少两 个感应式磁场传感器的误差值小于误差阈值，且至少四个以上感应式磁场 传感器的感应电动势大于探测阈值时，根据计算出的下一采样时刻的最优 位置和最优磁矩，生成水中磁性目标下一采样时刻的最优定位信息。
具体的，到下一个采样时刻，根据上一采样时刻的最优位置和最优磁矩， 随机生成一组水中磁性目标的位置和磁矩值形成的初始解集，作为初代种 群，采用差分进化算法进行寻优，得到下一采样时刻的最优位置和最优磁矩， 进而计算出下一采样时刻感应电动势测量值与感应电动势计算值(正演计 算值)之间的误差。
判断是否存在至少2个感应式磁场传感器感应电动势误差小于误差阈 值，同时，至少四个以上感应式磁场传感器的感应电动势大于探测阈值，若 是，则依据下一采样时刻的最优位置和最优磁矩，生成新的种群的最优定位 信息，再重复计算下一采样时刻；否则，判定水中磁性目标消失，重新由磁 通计算单元判断感应电动势是否符合预设条件。
采用本发明的上述技术方案和简化的传感器布置方式，利用时变磁梯 度法就可实现对水中磁性目标的初步定位。但由于测得的不是严格意义上 的同一位置点的磁场强度和磁场梯度，其解总是距离真实值存在一定差距 的。
在获得了水中磁性目标的位置和磁矩初值后，本发明将在步骤5中， 以磁场传感器组测得的最大感应电动势异常值为依据，通过建立反映磁场 定位解与真实水中磁性目标位置和磁矩参数逼近程度的目标函数，对水中 磁性目标的位置和磁矩进行进一步寻优，以实现目标精确定位。
如图4所示，本实施例中的定位系统与现有的静态磁梯度定位法相比 较，通过仿真测试可以看出，在定位距离超过130m时，本实施例中的定位 精度要明显由于现有的静态磁梯度定位法，本实施例中的定位系统能够克 服探测定位过程中易受海洋环境磁场干扰的缺点。
以上结合附图详细说明了本申请的技术方案，本申请提出了一种基于 矩形阵列的水中运动磁性目标探测定位系统，该系统包括：探测矩阵，磁 通计算单元以及初始探测单元；探测矩阵中设置有多个感应式磁场传感 器，感应式磁场传感器用于获取探测范围内因水中磁性目标移动产生的感 应电动势；磁通计算单元用于当判定感应电动势符合预设条件时，计算探 测矩阵中阵列单元所在位置处的磁通变化量；初始探测单元用于根据磁通变化量，通过导数运算，计算磁场时变量，并结合磁通变化量，计算水中 磁性目标的初始位置和初始磁矩，并生成初始定位信息，其中，初始定位 信息包括水中磁性目标的位置和类型，初始磁矩用于确定类型。通过本申 请中的技术方案，克服了探测定位过程中易受海洋环境磁场干扰的缺点。
本申请中的步骤可根据实际需求进行顺序调整、合并和删减。
本申请装置中的单元可根据实际需求进行合并、划分和删减。
尽管参考附图详地公开了本申请，但应理解的是，这些描述仅仅是示 例性的，并非用来限制本申请的应用。本申请的保护范围由附加权利要求 限定，并可包括在不脱离本申请保护范围和精神的情况下针对发明所作的 各种变型、改型及等效方案。

Claims (8)

1.一种基于矩形阵列的水中运动磁性目标探测定位系统，其特征在于，所述探测定位系统包括：探测矩阵，磁通计算单元以及初始探测单元；
所述探测矩阵中设置有多个感应式磁场传感器，所述感应式磁场传感器用于获取探测范围内因水中磁性目标移动产生的感应电动势；
所述磁通计算单元用于当判定所述感应电动势符合预设条件时，计算所述探测矩阵中阵列单元所在位置处的磁通变化量；
所述初始探测单元用于根据所述磁通变化量，通过导数运算，计算磁场时变量，并结合所述磁通变化量，计算所述水中磁性目标的初始位置和初始磁矩，并生成初始定位信息，其中，所述初始定位信息包括所述水中磁性目标的位置和类型，所述初始磁矩用于确定所述类型。
2.如权利要求1所述的基于矩形阵列的水中运动磁性目标探测定位系统，其特征在于，所述探测定位系统还包括：寻优单元；
所述寻优单元用于根据所述初始位置和初始磁矩，确定所述水中磁性目标的磁场定位初始解集，并根据所述磁场定位初始解集和预设优化目标函数，计算所述水中磁性目标的最优位置和最优磁矩，并生成最优定位信息，其中，所述最优定位信息包括所述水中磁性目标的最优位置和最优类型，所述最优磁矩用于确定所述最优类型。
3.如权利要求2所述的基于矩形阵列的水中运动磁性目标探测定位系统，其特征在于，所述预设优化目标函数的计算公式为：
Obj＝∑[(Uxc-Uxm)2+(UyC-Uym)2+(Uzc-Uzm)2]
式中，Uxc为感应式磁场传感器中磁通x分量Bx对应的感应线圈棒的电压计算值，Uxm为磁通x分量Bx对应的电压测量值，Uyc为感应式磁场传感器中磁通y分量By对应的感应线圈棒的电压计算值，Uym为磁通y分量By对应的电压测量值，Uzc为感应式磁场传感器中磁通z分量Bz对应的感应线圈棒的电压的计算值，Uzm为磁通z分量Bz对应的电压测量值。
4.如权利要求2所述的基于矩形阵列的水中运动磁性目标探测定位系统，其特征在于，所述寻优单元还用于：
根据上一采样时刻的最优位置和最优磁矩，确定下一采样时刻的磁场定位初始解集，并根据所述预设优化目标函数，计算下一采样时刻的最优位置和最优磁矩；
所述探测定位系统还包括：电动势计算单元以及判断单元；
所述电动势计算单元用于根据下一采样时刻的最优位置和最优磁矩，计算下一采样时刻对应的感应式磁场传感器的感应电动势计算值以及误差值，其中，所述误差值为下一采样时刻感应电动势计算值与感应电动势测量值之间的误差；
所述判断单元用于当判定至少两个感应式磁场传感器的误差值小于误差阈值，且至少四个以上感应式磁场传感器的感应电动势大于探测阈值时，根据计算出的下一采样时刻的最优位置和最优磁矩，生成所述水中磁性目标下一采样时刻的最优定位信息。
5.如权利要求1所述的基于矩形阵列的水中运动磁性目标探测定位系统，其特征在于，所述阵列单元由所述探测矩阵中的异常传感器组成，所述异常传感器为获取到的感应电动势大于探测阈值的感应式磁场传感器，所述预设条件为所述异常传感器的数量大于或等于数量阈值。
6.如权利要求5所述的基于矩形阵列的水中运动磁性目标探测定位系统，其特征在于，所述数量阈值为四个，所述预设条件还包括：
多个所述异常传感器中存在至少四个异常传感器可以组成矩形阵列，将所述组成的所述矩形阵列记作所述阵列单元。
7.如权利要求1所述的基于矩形阵列的水中运动磁性目标探测定位系统，其特征在于，所述初始探测单元中包括：导数计算单元，张量矩阵计算单元以及磁矩计算单元；
所述导数计算单元用于根据水中磁性目标的移动速度和所述磁通变化量，计算所述阵列单元所处位置几何中心处的磁通关于时间的导数；
所述张量矩阵计算单元用于根据所述导数，构建时变磁梯度张量矩阵；
所述磁矩计算单元用于根据所述时变磁梯度张量矩阵，计算所述初始位置和初始磁矩，其中，所述初始磁矩的计算公式为：
Figure FDA0003447910550000031
Figure FDA0003447910550000032
Figure FDA0003447910550000033
式中，mx0、my0、mz0为所述水中磁性目标在三维坐标系中的三分量磁矩，R0为水中磁性目标距离探测原点的距离，所述探测原点定义为矩形阵列的几何中心，μ0为真空磁导率，x0、y0、z0为所述水中磁性目标相对所述探测原点的坐标，G为所述时变磁梯度张量矩阵。
8.如权利要求1所述的基于矩形阵列的水中运动磁性目标探测定位系统，其特征在于，所述感应式磁场传感器由三个相互正交的感应线圈棒组成。
CN202111654555.4A 2021-12-30 2021-12-30 一种基于矩形阵列的水中运动磁性目标探测定位系统 Pending CN114966866A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111654555.4A CN114966866A (zh) 2021-12-30 2021-12-30 一种基于矩形阵列的水中运动磁性目标探测定位系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111654555.4A CN114966866A (zh) 2021-12-30 2021-12-30 一种基于矩形阵列的水中运动磁性目标探测定位系统

Publications (1)

Publication Number Publication Date
CN114966866A true CN114966866A (zh) 2022-08-30

Family

ID=82974524

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111654555.4A Pending CN114966866A (zh) 2021-12-30 2021-12-30 一种基于矩形阵列的水中运动磁性目标探测定位系统

Country Status (1)

Country Link
CN (1) CN114966866A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115963038A (zh) * 2022-12-14 2023-04-14 中国科学院空间应用工程与技术中心 基于空间微重力条件的磁性颗粒运动轨迹测量系统及方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115963038A (zh) * 2022-12-14 2023-04-14 中国科学院空间应用工程与技术中心 基于空间微重力条件的磁性颗粒运动轨迹测量系统及方法
CN115963038B (zh) * 2022-12-14 2023-07-28 中国科学院空间应用工程与技术中心 基于空间微重力条件的磁性颗粒运动轨迹测量系统及方法

Similar Documents

Publication Publication Date Title
CN110095738B (zh) 磁梯度仪参数设计方法
EP1092163B1 (en) Apparatus for detecting metals
JP2954063B2 (ja) 改良されたダイポールモーメント検出器およびローカライザ
CN112213768B (zh) 一种联合震源机制反演的地面微地震定位方法及系统
CN109579827B (zh) 一种基于弧形阵列的磁性目标探测和定位方法
CN109425906B (zh) 一种磁异常探测矢量磁目标识别方法
CN108802722A (zh) 一种基于虚拟谱的弱目标检测前跟踪方法
KR20010033029A (ko) 자성 감지 장치의 직접 관측에 기초한 자성체
CN111399066B (zh) 一种基于正交基函数处理标量磁异常梯度信号的方法
CN110850483B (zh) 一种基于电场电极阵列布置的水中目标探测和定位方法
CN111522835A (zh) 一种基于数据库特征匹配的多磁性目标位置探测方法
CN114966866A (zh) 一种基于矩形阵列的水中运动磁性目标探测定位系统
CN110058222A (zh) 一种基于传感器选择的双层粒子滤波检测前跟踪方法
CN114488326B (zh) 一种提高探测航磁数据地质体边界能力的方法及系统
CN112668146B (zh) 一种基于欧拉反褶积法实用性改进的场源位置估算方法
Liu et al. Enhanced magnetic imaging for industrial metal workpiece detection through the combination of electromagnetic induction and magnetic anomalies
Zhang et al. Support vector machine and neural network classification of metallic objects using coefficients of the spheroidal MQS response modes
Šimić et al. Inversion-based magnetic polarizability tensor measurement from time domain EMI data
CN115236755B (zh) 基于张量特征值的航磁异常边界检测方法、装置
CN110161460A (zh) 基于微震空间组网的震源精确定位方法
CN111695087B (zh) 用于复合探雷目标检测的方法及系统
CN113552635B (zh) 固源瞬变电磁数据融合隧道三维超前预报方法及系统
Chen et al. Improved differential evolution algorithm for multi-target response inversion detected by a portable transient electromagnetic sensor
CN109579828B (zh) 一种基于弧形磁场传感器阵列的定位系统及其使用方法
CN117452342B (zh) 一种基于极化特征的箔条干扰检测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination