基于最小二乘两步法的地磁矢量测量多源误差校正方法与流程

1.本发明属于地磁矢量测量技术领域，具体涉及一种基于最小二乘两步法的地磁矢量测量多源误差校正方法。


背景技术：

2.地磁场是地球固有的一种矢量场，进行地磁矢量场测量获得的信息更加丰富，其具有地磁标量测量不可比拟的优势。在进行移动式地磁矢量测量过程中，需要将地磁矢量测量系统搭载于测量载体之上(载体可以为测量船/艇、有人/无人飞机、拖体、auv等)，通常地磁矢量测量系统包含三分量磁力仪、高精度惯导(姿态仪)、数据采集模块、gnss模块、数据存储模块等。在载体行进过程中，三分量磁力仪和高精度惯导分别测量不同点位的三分量磁场值和姿态角信息，然后根据姿态角信息对三分量磁场进行坐标旋转，可得到地理坐标系下的地磁矢量。
3.当前，用于三分量磁力测量的仪器主要为三轴磁通门传感器，受限于材料和加工工艺的限制，其存在零偏误差、标度系数误差、三轴非正交误差，这些误差会影响三轴磁通门传感器的测量精度。理想情况下，三分量磁力仪和姿态仪坐标系处于对准状态，可将三分量磁场准确地变换至地理坐标系下；在实际工程中，三分量磁力仪坐标系和姿态仪坐标系均不可视，单纯依靠机械安装的方式难以保证二者坐标系精确对准，即存在非对准误差；非对准误差会导致三分量磁场姿态坐标变换出现偏差，影响地理坐标系下矢量磁场的测量精度。通常的测量载体，在其动力装置、电气设备、载体结构中不可避免含有铁磁性材料，其会产生包含硬磁、感应磁场和涡流磁场的干扰磁场。
4.上述三轴磁通门传感器转向差、非对准误差、载体干扰磁场，会严重影响地磁矢量测量的精度，对这些误差和干扰磁场进行校正和补偿是移动式地磁矢量测量中的关键技术问题，众多学者开展了相关的研究。
5.在三轴磁通门传感器转向差校正方面，庞学亮等提出了基于遗传算法三轴磁传感器误差校正方法，首先构造目标函数，利用遗传算法求取目标函数的最优解，然后对测量数据进行校正，通过数值仿真表明，该方法可以稳定地收敛并可以达到很高的误差参数估计精度(庞学亮,林春生.基于遗传算法三轴磁传感器校正系数求解[j].探测与控制学报,2017,39(1):42-45,51.)。在非对准误差校正方面，庞鸿锋等人提出了一种基于绕轴法的地磁矢量系统非对准校正方法，该该方法要求将三轴磁传感器和加速度计封装在无磁正六面体中，将无磁正六面体放置在无磁转台上，分别绕无磁正六面体x轴和z轴旋转，分别计算出无磁正六面体到三轴磁传感器、加速度计的非对准角，然后确定三轴磁传感器和加速度计之间的坐标转换关系(庞鸿锋，张琦，万成彪等.基于绕轴法的地磁矢量系统非对准校正方法，专利申请号：201410740758.9)。在载体磁干扰补偿方面，刘继昊等人提出了一种载体磁干扰矢量补偿方法，该方法提出了一种改进的磁干扰参数判定准则，但在半实物仿真实验中，补偿后的分量误差仍高达千nt量级，其需要将三轴磁传感器、三轴陀螺仪安装在铝板上，且其假定二者坐标系对准，然后在无干扰源和有干扰源条件下完成测量，在实际载体测
量条件下并不易实现(刘继昊,李夕海,张莹,等.一种通用的地磁测量中载体干扰矢量补尝方法[j].仪器仪表学报,2020,41(11):112-118.)。
[0006]
上述方法均是对单项测量误差进行校正或补偿，在实际测量中如逐项误差进行校正和补偿，不仅程序复杂、可操作性较低，而且各项误差可能存在耦合关系，分别校正可能效果较差。在误差综合校正方面，张琦等提出了一种基于拉格朗日乘数法的地磁矢量测量系统误差校正方法，该方法借助无磁正六面箱体和无磁平台，建立包含磁传感器误差、非对准误差、软磁和硬磁的线性方程，利用磁总场约束的方法求解误差参数，实现误差校正(张琦，万成彪，潘孟春，基于拉格朗日乘数法的地磁矢量测量系统误差校正方法，专利申请号：cn201710067382.3)。李婷等提出了基于阻尼粒子群优化的地磁场测量误差补偿方法，首先建立了包含磁传感器误差和载体磁干扰得误差补偿模型，利用阻尼粒子群算法和总场约束的方法对模型参数进行估计，半实物试验表明可对磁总场误差取的较好的补偿效果(李婷,张金生,王仕成,等.基于阻尼粒子群优化的地磁场测量误差补偿[j].仪器仪表学报,2017,38(10):2446-2452)。
[0007]
综上可知，现有地磁矢量测量误差校正和补偿还存在一些问题：1)针对单项误差开展校正补偿，不能实现多项误差的综合校正补偿，在实际中的应用效率和可操作性较低；2)利用磁总场作为误差参数计算的约束条件，其针对磁总场误差校正效果较好，但对地磁矢量误差校正效果较差；3)方法局限性较高，需要复杂的设备工装，在实际测量条件下难以开展应用。
[0008]
本背景技术所公开的上述信息仅仅用于增加对本技术背景技术的理解，因此，其可能包括不构成本领域普通技术人员已知的现有技术。


技术实现要素：

[0009]
本发明针对现有技术中存在的上述问题，提供一种基于最小二乘两步法的地磁矢量测量多源误差校正方法。
[0010]
为达到上述技术目的，本发明采用以下技术方案实现：一种基于最小二乘两步法的地磁矢量测量多源误差校正方法，所述方法为：
[0011]
s1、建立测量模型：
[0012]bm
＝hbe+b0(公式2)；
[0013]
其中，h＝kc(i+m)，b0＝kcbh+b0，i为单位矩阵，bm＝[b
mx
,b
my
,b
mz
]
t
为三轴磁通门传感器的输出值，be＝[b
ex
,b
ey
,b
ez
]
t
为地磁场在正交磁力仪坐标系下的投影矢量，为载体感应磁场系数矩阵，为三轴磁通门传感器非正交误差矩阵，α、β、γ为非正交角，为三轴磁通门传感器标度系数矩阵，bh＝[b
hx
,b
hy
,b
hz
]
t
为载体硬磁在正交磁力仪坐标系下的投影，b0＝[b
0x
,b
0y
,b
0z
]
t
为三轴磁通门传感器零偏误差向量；
[0014]
s2、根据测量模型得到误差校正模型：
[0015][0016]
其中，是惯导坐标系下的地磁投影矢量，是正交磁力仪坐标系至惯导坐标系的旋转矩阵，g为系数矩阵，b0为系数向量；
[0017]
s3、采集标定数据；
[0018]
s4、利用igrf模型计算得到背景磁场，并将其旋转至惯导坐标系；
[0019]
利用igrf模型计算的地理坐标系下的地磁矢量值，来近似表示此处真实背景地磁场值；设利用igrf模型计算出的标定区域地理坐标系下的磁场矢量为根据标定过程中姿态仪输出的姿态角，对进行坐标旋转，可得到惯导坐标系下的磁场矢量：
[0020]
其中，表示惯导坐标系下的地磁矢量，表示地理坐标系到惯导坐标系的转换矩阵，
[0021]
r、p、h分别表示姿态仪输出的横滚角、俯仰角、航向角；
[0022]
s5、采用地磁分量约束法，利用最小二乘获得系数矩阵g
′
，进而获得矩阵g：
[0023]
令有
[0024]
其中，系数矩阵g为3
×
3矩阵，向量b0为3
×
1矩阵，设1矩阵，设得到：
[0025]
其中，
[0026]
由于采集多组标定数据：bm(i)、(i＝1,2,
…
,n,n≥12)，将(公式8)写为：
[0027]
其中，
根据最小二乘法原理，利用(公式9)获得g
′
的最佳估计；
[0028]
获得系数矩阵g；
[0029]
s6、根据最小二乘原理，获得向量b0：
[0030]
令g4＝[g
14
,g
24
,g
34
]
t
，其为已知量，根据(公式8)，得到gb0＝g4(公式11)；
[0031]
根据最小二乘法原理获得b0的估计值；
[0032]
s7、根据(公式5)计算惯导坐标系下的地磁矢量
[0033]
s8、将惯导坐标系下的地磁矢量旋转至地理坐标系：
[0034]
根据姿态仪输出的姿态角对进行坐标旋转，得到地理坐标系下地磁矢量进行坐标旋转，得到地理坐标系下地磁矢量
[0035]
其中，是从惯导坐标系到地理坐标系的旋转矩阵：
[0036][0037]
如上所述的基于最小二乘两步法的地磁矢量测量多源误差校正方法，地磁矢量测量系统安装至移动载体后，载体坐标系与姿态仪坐标系一致。
[0038]
如上所述的基于最小二乘两步法的地磁矢量测量多源误差校正方法，所述方法建立了包含三轴磁通门传感器误差、非对准误差、软磁和硬磁干扰的误差校正模型。
[0039]
如上所述的基于最小二乘两步法的地磁矢量测量多源误差校正方法，采集标定数据的方法为：测区附近选择一个地磁梯度低于设定值、周围磁干扰小于设定值的开阔区域，选择靠近该区域中心的某一点，以该点为中心，令载体进行机动航行采集标定数据。
[0040]
如上所述的基于最小二乘两步法的地磁矢量测量多源误差校正方法，所述机动方式为：以所选点为中心，分别令载体以测量航速和最小转弯半径绕该点做顺时针和逆时针航行，各旋转一圈。
[0041]
如上所述的基于最小二乘两步法的地磁矢量测量多源误差校正方法，所述机动方式为：以所选点为中心，令载体做“十”字往返航行。
[0042]
如上所述的基于最小二乘两步法的地磁矢量测量多源误差校正方法，所述“十”字每条边长度为20倍载体长度。
[0043]
与现有技术相比，本发明的优点和积极效果是：本发明基于最小二乘两步法的地
磁矢量测量多源误差校正方法，首先建立包含三轴磁通门传感器误差、非对准误差、软磁和硬磁的误差综合补偿模型，然后令测量载体按照一定的动作航行采集标定数据，采用地磁分量约束法，利用两步最小二乘算计算出各项补偿参数，最后对实测数据进行补偿。本发明可对上述各项误差进行综合补偿，不需要其它辅助设备、操作简单、计算速度快，可在实际移动式地磁矢量测量条件下开展应用。本发明可对地磁矢量测量中的三分量磁力仪误差、非对准误差、载体磁场干扰实现快速校正补偿。
附图说明
[0044]
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0045]
图1为地磁矢量测量系统坐标系示意图。
[0046]
图2为磁力仪坐标系相互关系图。
[0047]
图3为载体机动示意图。
[0048]
图4为基于最小二乘两步法的地磁矢量测量多源误差校正方法的流程图。
具体实施方式
[0049]
下面参照附图来描述本发明的优选实施方式。本领域技术人员应当理解的是，这些实施方式仅仅用于解释本发明的技术原理，并非旨在限制本发明的保护范围。
[0050]
需要说明的是，在本发明的描述中，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方向或位置关系的术语是基于附图所示的方向或位置关系，这仅仅是为了便于描述，而不是指示或暗示所述装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
[0051]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。
[0052]
地磁矢量测量系统中有如下坐标系，惯导坐标系ox
pypzp
、理想正交磁力仪坐标系oxmymzm，此二者均为正交坐标系，二者示意图如图1；另有，由三轴磁通门传感器三个敏感轴构成的实际磁力仪坐标系ox
′mym′z′m，由于三轴磁通门传感器存在非正交误差，实际磁力仪坐标系为非正交坐标系，其与理想正交磁力仪坐标系的对应关系如图2。地磁矢量测量系统安装至移动载体后，设载体坐标系与惯导坐标系一致。
[0053]
如图4所示，本实施例的方法包括如下步骤：
[0054]
s1、建立测量模型：
[0055]
在上述坐标系关系下，当只考虑三轴磁通门传感器误差、软磁和硬磁干扰时，可建立移动式地磁矢量测量条件下的测量模型：
[0056]bm
＝kc(be+mbe+bh)+b0ꢀꢀꢀ
(公式1)；
[0057]
其中，bm＝[b
mx
,b
my
,b
mz
]
t
为三轴磁通门传感器的输出值，be＝[b
ex
,b
ey
,b
ez
]
t
为地磁
场在正交磁力仪坐标系下的投影矢量，为载体感应磁场系数矩阵，为三轴磁通门传感器非正交误差矩阵(α、β、γ为非正交角，如图2中所示)，为三轴磁通门传感器标度系数矩阵，bh＝[b
hx
,b
hy
,b
hz
]
t
为载体硬磁在正交磁力仪坐标系下的投影，b0＝[b
0x
,b
0y
,b
0z
]
t
为三轴磁通门传感器零偏误差向量。
[0058]
对(公式1)简化得到：bm＝hbe+b0ꢀꢀꢀ
(公式2)；
[0059]
其中，h＝kc(i+m)，b0＝kcbh+b0。
[0060]
s2、根据测量模型得到误差校正模型：
[0061]
对(公式2)进行变换，得到误差校正模型：be＝h-1
(b
m-b0)
ꢀꢀꢀ
(公式3)；
[0062]
其中，be为地磁场在正交磁力仪坐标系下的投影矢量，实际上，正交磁力仪坐标系和惯导坐标系之间存在非对准误差，要想通过高精度惯导(姿态仪)输出的姿态角和航向角把be为旋转至地理坐标系，首先应该使正交磁力仪坐标系和惯导坐标系对准，即首先把be旋转至惯导坐标系：
[0063]
其中，是地磁场旋转至惯导坐标系下的地磁矢量，是正交磁力仪坐标系至惯导坐标系的旋转矩阵。
[0064]
令则有：
[0065]
根据(公式5)，如果求得系数矩阵g和系数向量b0就可以通过三轴磁通门传感器的输出值，获得惯导坐标系下的地磁矢量真实值。
[0066]
在实际的移动式地磁矢量测量中，在测量任务开始前或结束后，需要令载体按照一定动作进行航行，测量系统开机采集标定数据，用于误差校正参数的计算。
[0067]
s3、选择标定区域，采集标定数据bm(i)及姿态角a
(r,p,h)
(i)，(i＝1,2,...,n,n≥12)；。
[0068]
在实际操作中，在测区附近选择一个地磁梯度较小(地磁梯度小于设定值)、周围无明显磁干扰源(磁干扰小于设定值)的开阔区域，选择靠近该区域中心的某一点，以该点为中心，令载体进行机动航行采集标定数据，机动方式为：1)以所选点为中心，分别令载体以测量航速和最小转弯半径绕该点做顺时针和逆时针航行，各旋转一圈；2)以所选点为中心，令载体做“十”字往返航行，“十”字每条边长度按照载体长度定义，一般选择20倍载体长度，图1中“l”表示载体长度。载体机动轨迹示意图如图3。
[0069]
s4、利用igrf模型计算得到背景磁场，并将其旋转至惯导坐标系。
[0070]
由于标定区域地磁梯度较小、周围无明显磁干扰源，可利用igrf模型计算的地理
坐标系下的地磁矢量值，来近似表示此处真实背景地磁场值；设利用igrf模型计算出的标定区域地理坐标系下的磁场矢量为根据标定过程中高精度惯导输出的姿态角，对进行坐标旋转，可得到惯导坐标系下的磁场矢量，过程如下：
[0071][0072]
其中，表示背景磁场转转至惯导坐标系下的地磁矢量，表示地理坐标系到惯导坐标系的转换矩阵，有：
[0073]
r、p、h分别表示高精度惯导输出的横滚角、俯仰角、航向角。
[0074]
s5、采用地磁分量约束法，利用最小二乘，获得系数矩阵g
′
，进而获得矩阵g：
[0075]
联合(公式5)、(公式6)，若有则说明系数矩阵g和向量b0是正确的；但实际中，磁传感器和高精度惯导均存在测量噪声，和不会完全相等，可令有：
[0076][0077]
系数矩阵g为3
×
3矩阵，向量b0为3
×
1矩阵，设1矩阵，设将其展开，得到：
[0078][0079]
其中，
[0080]
在实际标定过程中，会采集多组标定数据，有：bm(i)、(i＝1,2,...,n,n≥12)。
[0081]
(公式8)可以改写为：
[0082]
其中，
根据最小二乘法原理，利用(公式9)可获得g
′
的最佳估计，即
[0083][0084]
获得g
′
后，即获得了系数矩阵g。
[0085]
s6、根据最小二乘原理，获得向量b0：
[0086]
令g4＝[g
14
,g
24
,g
34
]
t
，其为已知量，根据(公式8)，有：gb0＝g4(公式11)
[0087]
同样，根据最小二乘法原理，可获得b0的估计值，即
[0088]
b0＝(g
t
g)-1gt
g4ꢀꢀꢀ
(公式12)。
[0089]
s7、根据(公式5)计算惯导坐标系下的地磁矢量根据(公式5)对实测数据进行校正。
[0090]
通过上面推导和计算，已经得到系数矩阵g和向量b0，可根据(公式5)对实际测量作业获得的地磁矢量数据bm进行实时在线补偿或离线补偿，获得惯导坐标系下的地磁矢量
[0091]
s8、将惯导坐标系下的地磁矢量旋转至地理坐标系：
[0092]
根据高精度惯导输出的姿态角对进行坐标旋转，便可得到地理坐标系下地磁矢量bg＝[b
gx
,b
gy
,b
gz
]
t
，过程如下：
[0093][0094]
其中，是从惯导坐标系到地理坐标系的旋转矩阵：
[0095]
r、p、h分别表示高精度惯导输出的横滚角、俯仰角、航向角。
[0096]
本实施例梳理了各坐标系的对应关系，以惯导坐标系为基准，建立了包含三轴磁通门传感器误差、非对准误差、软磁和硬磁干扰的误差校正模型；本实施例可一次性对多项误差和干扰进行校正补偿，相比单项误差分别校正，可操作性更强、效率更高。
[0097]
本实施例在标定数据采集过程中，令载体分别以顺、逆时针绕“o”航行和往返“十”航行的方式进行；其相比传统“八”方位法更节省时间，可减小补偿过程中地磁日变的影响；相比单独绕“o”法更能覆盖测量中可能出现的姿态，采集的标定数据更有效。
[0098]
本实施例在误差参数求解中采用地磁分量约束法，求解出的误差参数对分量误差补偿精度更高；相比地磁总场约束，计算的补偿参数更准确，对地磁矢量数据的补偿效果更好。
[0099]
本实施例采用两步最小二乘算法分别求解校正矩阵g和向量b0，本实施例采用最小二乘法计算补偿参数，相比一般迭代计算方法、卡尔曼滤波等方法，该方法计算速度更快，效率更高。
[0100]
本实施例不需要辅助设备，实施过程简单，实用性较强。
[0101]
本实施例通用性较强，可应用于测量船/艇、有人/无人飞机、拖体、auv等载体条件下的移动式地磁矢量测量的多源误差综合补偿。
[0102]
以上实施例仅用以说明本发明的技术方案，而非对其进行限制；尽管参照前述实施例对本发明进行了详细的说明，对于本领域的普通技术人员来说，依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或替换，并不使相应技术方案的本质脱离本发明所要求保护的技术方案的精神和范围。

技术特征：
1.一种基于最小二乘两步法的地磁矢量测量多源误差校正方法，其特征在于，所述方法为：s1、建立测量模型：b
m
＝hb
e
+b0ꢀꢀꢀꢀ
(公式2)；其中，h＝kc(i+m)，b0＝kcb
h
+b0，i为单位矩阵，b
m
＝[b
mx
,b
my
,b
mz
]
t
为三轴磁通门传感器的输出值，b
e
＝[b
ex
,b
ey
,b
ez
]
t
为地磁场在正交磁力仪坐标系下的投影矢量，为载体感应磁场系数矩阵，为三轴磁通门传感器非正交误差矩阵，α、β、γ为非正交角，为三轴磁通门传感器标度系数矩阵，b
h
＝[b
hx
,b
hy
,b
hz
]
t
为载体硬磁在正交磁力仪坐标系下的投影，b0＝[b
0x
,b
0y
,b
0z
]
t
为三轴磁通门传感器零偏误差向量；s2、根据测量模型得到误差校正模型：其中，是惯导坐标系下的地磁投影矢量，是正交磁力仪坐标系至惯导坐标系的旋转矩阵，g为系数矩阵，b0为系数向量；s3、采集标定数据；s4、利用igrf模型计算得到背景磁场，并将其旋转至惯导坐标系；利用igrf模型计算的地理坐标系下的地磁矢量值，来近似表示此处真实背景地磁场值；设利用igrf模型计算出的标定区域地理坐标系下的磁场矢量为根据标定过程中姿态仪输出的姿态角，对进行坐标旋转，可得到惯导坐标系下的磁场矢量：其中，表示惯导坐标系下的地磁矢量，表示地理坐标系到惯导坐标系的转换矩阵，r、p、h分别表示姿态仪输出的横滚角、俯仰角、航向角；s5、采用地磁分量约束法，利用最小二乘获得系数矩阵g
′
，进而获得矩阵g：令有
其中，系数矩阵g为3
×
3矩阵，向量b0为3
×
1矩阵，设1矩阵，设得到：其中，由于采集多组标定数据：b
m
(i)、将(公式8)写为：其中，其中，根据最小二乘法原理，利用(公式9)获得g
′
的最佳估计；获得系数矩阵g；s6、根据最小二乘原理，获得向量b0：令g4＝[g
14
,g
24
,g
34
]
t
，其为已知量，根据(公式8)，得到gb0＝g4(公式11)；根据最小二乘法原理获得b0的估计值；s7、根据(公式5)计算惯导坐标系下的地磁矢量s8、将惯导坐标系下的地磁矢量旋转至地理坐标系：根据姿态仪输出的姿态角对进行坐标旋转，得到地理坐标系下地磁矢量b
g
＝[b
gx
,b
gy
,b
gz
]
t
，其中，是从惯导坐标系到地理坐标系的旋转矩阵：
2.根据权利要求1所述的基于最小二乘两步法的地磁矢量测量多源误差校正方法，其特征在于，地磁矢量测量系统安装至移动载体后，载体坐标系与姿态仪坐标系一致。3.根据权利要求1所述的基于最小二乘两步法的地磁矢量测量多源误差校正方法，其特征在于，所述方法建立了包含三轴磁通门传感器误差、非对准误差、软磁和硬磁干扰的多源误差校正模型。4.根据权利要求1所述的基于最小二乘两步法的地磁矢量测量多源误差校正方法，其特征在于，采集标定数据的方法为：测区附近选择一个地磁梯度低于设定值、周围磁干扰小于设定值的开阔区域，选择靠近该区域中心的某一点，以该点为中心，令载体进行机动航行采集标定数据。5.根据权利要求3所述的基于最小二乘两步法的地磁矢量测量多源误差校正方法，其特征在于，所述机动方式为：以所选点为中心，分别令载体以测量航速和最小转弯半径绕该点做顺时针和逆时针航行，各旋转一圈。6.根据权利要求3或4所述的基于最小二乘两步法的地磁矢量测量多源误差校正方法，其特征在于，所述机动方式为：以所选点为中心，令载体做“十”字往返航行。7.根据权利要求6所述的基于最小二乘两步法的地磁矢量测量多源误差校正方法，其特征在于，所述“十”字每条边长度为20倍载体长度。

技术总结
本发明提供了一种基于最小二乘两步法的地磁矢量测量多源误差校正方法，首先建立包含三轴磁通门传感器误差、非对准误差、软磁和硬磁的误差综合补偿模型，然后令测量载体按照一定的动作航行采集标定数据，利用地磁分量约束法,采用两步最小二乘算计算出各项补偿参数，最后对实测数据进行补偿。本发明可对上述各项误差进行综合补偿，不需要其它辅助设备、操作简单、计算速度快，可在实际移动式地磁矢量测量条件下开展应用。量条件下开展应用。量条件下开展应用。


技术研发人员：罗建刚 李海兵 刘静晓 张峰 李海虎 郭子伟 付碧波 王友东
受保护的技术使用者：青岛海洋科技中心
技术研发日：2023.02.23
技术公布日：2023/7/11