公开/公告号CN109297530A
专利类型发明专利
公开/公告日2019-02-01
原文格式PDF
申请/专利权人 广州海洋地质调查局;
申请/专利号CN201810886541.7
申请日2018-08-06
分类号
代理机构广州君咨知识产权代理有限公司;
代理人吴敏珍
地址 510000 广东省广州市环市东路477号大院
入库时间 2024-02-19 06:57:51
法律状态公告日
法律状态信息
法律状态
2019-08-23
授权
授权
2019-03-01
实质审查的生效 IPC(主分类):G01D21/02 申请日:20180806
实质审查的生效
2019-02-01
公开
公开
技术领域
本发明涉及海洋调查技术领域,具体是一种深海全水柱多环境要素信息融 合方法及处理终端。
背景技术
深海环境基线调查需要同时获取多种有效、高分辨率的水柱环境要素信息, 提高对海洋环境特征的辨识能力。而目前采用温盐深测量、走航ADCP观测和生 物拖网等常规调查方法获得的数据并没有进行融合,垂向分辨率低,缺少部分关 键的环境要素,获得的数据处理结果可交互性差。
发明内容
针对现有技术的不足,本发明的目的之一是提供一种深海全水柱多环境要 素信息融合方法,其能够解决将多种环境要素信息融合的问题。
本发明的目的之二提供一种处理终端,其能够解决将多种环境要素信息融 合的问题,解决处理生成满足环境基线科研要求的高质量多环境要素数据的问 题。
实现本发明目的之一的技术方案为:一种深海全水柱多环境要素信息融合方 法,其包括以下步骤:
步骤S1:获取包括压力、温度、盐度、电导率、溶解氧、PH、叶绿素、声 速、船载ADCP的声速、垂向速度、速度剪切、相对流速、相对回波强度和生物 量数据的原始数据,其中压力、温度、盐度、电导率和声速直接由CTD获得,溶 解氧、PH和叶绿素均通过水样观测获得,船载ADCP的声速由船载ADCP根据温 度探头获得,垂向速度、速度剪切、相对流速由LADCP测得,相对回波强度由 LADCP和船载ADCP获得,生物量数据通过生物拖网获得;
步骤S2:根据获得的原始数据进行数据处理,包括:
001)压力修正,修正公式为:压力真实值=原始压力值-甲板压力,其中原 始压力值为步骤S1中的压力,甲板压力为CTD回收后在甲板处所测得的压力;
002)温度修正,修正公式为:offset=b*residual/n,式中,*为乘法运算, b表示“从航次校正时间起至当前对温度修正之间的天数”,residual表示“航 次前与航次后的温度差”,n表示“航次前与航次后温度校正的间隔天数”;
003)盐度修正,首先对电导率数据进行校正,校正公式如①:
式中,n表示校正所用的水化学分析样本数,αi表示第i个水化学分析样>i表示第i个水化学分析样本由水化学分析实验获得的>
根据公式②获得电导率真实值:
电导率真实值=电导率读数*slope------②
电导率读数是由步骤S1中获得的电导率,获得电导率真实值后,根据海水 状态方程计算出盐度,该盐度即为修正后的盐度值;
004)溶解氧、PH和叶绿素修正,溶解氧的修正首先根据公式③计算出斜率 修正参数:
式中,slopeDO表示斜率修正参数,DO水样表示在某一深度的水化学样本通>CTD表示与DO水样同一深度由CTD获得的>水样为直接从步骤S1中获得的溶解氧原始数据中读取到同一深度下>
获得斜率修正参数后,按公式④计算:
溶解氧真实值=溶解氧读数*slopeDO------④
溶解氧读数是由CTD获得的从海面到海底的所有溶解氧的原始数据,溶解 氧真实值即是斜率修正后的溶解氧;
PH和叶绿素的修正同样按公式③和④进行计算,获得修正后的PH和叶绿 素;
005)声速校正,获得校正后的相对流速值,根据公式⑤进行计算:
Vcorrected=Vrelative*(Cctd/CADCP)------⑤
式中,Vrelative表示未校正的相对流速,由步骤S1中获得,CADCP表示船载>ctd是步骤S1中获得的声速,由步骤S1中获>corrected表示因声速引起误差而经过校正后的相对流速值;
006)低散射校正,获得校正后的速度剪切值,根据方程组⑥进行计算:
V修正=V0+ΔV------⑥
式中,EI表示相对回波强度,由步骤S1获得,k为比例系数,为常数,ΔV 表示修正值,V0表示LADCP获得的速度剪切,由步骤S1获得,V修正表示经过低>
007)绝对流速计算,首先通过剪切法计算出参考流速,计算公式如⑦:
式中,T表示LADCP下放的总时长,t表示时间,Vship表示船舶航速,Vrefer表示参考流速,
获得Vrefer后按公式⑧计算,获得绝对流速Vabsolute:
Vabsolute=Vrelative+Vrefer------⑧
式中,Vabsolute表示绝对流速;
008)双LADCP绝对流速平均:将下看LADCP按步骤007)计算获得绝对流 速,并按照深度单元进行平均,获得V′absolute(i)(i=1,2,3……,n);同时,将上看LADCP>absolute(λ)(λ=1,2,3…...,m);将V″absolute(λ)(λ=1,2,3......,m)线性插值到下看LADCP>其中,m表示上看LADCP的深度单元个数,n表示下 看LADCP的深度单元个数;
009)去潮:通过软件求取正压潮剖面,从步骤007)获得的绝对流速Vabsolute减去所述正压潮剖面,即获得去潮后的结果;
010)理查森数计算,获得理查森数Ri,按公式⑨进行计算:
式中,g表示重力加速度,m表示m米深度,n表示n米深度,ρm表示m米>n表示n米深度处的海水密度,um、vm分别表示m米深度>n、vn分别表示n米深度>
011)垂向扩散系数计算,获得垂向扩散系数,计算公式如公式⑩所示:
式中,v0取10-2m2/s,α取5,vb取10-4m2/s,kb取10-5m2/s,KT是垂向>
012)CTD深度数据低通滤波:采用截断频率为6Hz的低通滤波器对CTD深 度数据进行低通滤波,得到6Hz的CTD深度数据;
013)CTD垂向运动速度计算,获得CTD的垂向运动速度ωCTD,按公式
式中,z(t)表示012)中获得的6Hz的CTD深度数据,t表示时间;
014)绝对垂向速度计算,获得绝对垂向速度,按公式
式中,ωADCP表示垂向速度,由步骤S1中通过LADCP获得的垂向速度,ω表>
015)双LADCP绝对垂向速度平均,获得平均绝对垂向速度,即将包括上看 LADCP和下看LADCP两台设备中的下看LADCP按公式
016)绝对回波强度计算,获得绝对回波强度,首先根据公式
式中,B是船载ADCP盲区,为常数;P是宣称脉冲长度,为常数;W是深 度单元厚度,N是深度单元个数,θ是波束角,cave是船载ADCP某一深度单元的>
根据公式
式中,K1c、a、b和c2均是船载ADCP的出厂参数,均为常数,Vs是整个投放>1为船载ADCP发>
根据公式
其中,
P2=1-1.37×10-4z+6.2×10-9z2,P3=1-1.383×10-5z+4.9×10-10z2,
式中,PH、S、T和z为船载ADCP在第i个深度单元相应的海水PH、盐度、温度、 深度值,由步骤S1获得,ci表示第i个深度单元对应的声速,freq表示船载ADCP>
然后,根据公式
式中,K2是系统噪声因子,为常数,Ks是与船载ADCP频率有关的系统固>x是当前船载ADCP探头所处位置的海水温度,由步骤S1获得,>a为步骤S1中LADCP测量获得的相对回波强度,Er为有效剖面末端的RSSI值,>V是绝对回波强度;
017)生物量反演,获得反演生物量,选取相同数量的若干个生物量数据和 绝对回波强度数据进行数据拟合,根据拟合公式
式中,DW为步骤S1中获得的生物量数据,
经过公式
步骤S3:获得上述经过处理后的各个数据后,对经过步骤S2处理获得的数 据进行数据融合,生成统一公共深度轴下的数据表格。
进一步地,所述步骤008)中的按照深度单元平均,其具体过程为:
步骤a),以Z'深值为间隔将从上看LADCP初始位置开始至海底进行划分,得>深值为间隔将从下看LADCP初始位置开始直至海底进行划分,得>
其中,m和n均为正整数,Z'深和Z”深分别为上看LADCP和下看LADCP采集数>
步骤b),首先,获取上看LADCP的n个深度单元的每个深度单元范围内所有 的Vabsolute值,进行算术平均得到V′absolute(i)(i=1,2,3......,n)值;以及获取下看LADCP>absolute值,进行算术平均得到>absolute(λ)(λ=1,2,3......,m)值;按深度将V″absolute(λ)线性插值到对应上看LADCP的>absolute(i)(i=1,2,3...…,n);然后,将同一深度单元内包含线>absolute(i)和上看LADCP的V′absolute(i)进行算术平均,得到该深度单元>最终得到n个深度单元的各个深度 单元的平均绝对流速
其中,如果m=n,则V″absolute(λ)刚好一一线性插值到对应上看LADCP的各深>absolute(λ)>
进一步地,所述步骤015)中的按照深度单元平均,其具体过程为:
步骤a),以Z'深值为间隔将从上看LADCP初始位置开始至海底进行划分,得>深值为间隔将从下看LADCP初始位置开始直至海底进行划分,得>
其中,m和n均为正整数,Z'深和Z”深分别为上看LADCP和下看LADCP采集数>
步骤b),首先,获取上看LADCP的n个深度单元的每个深度单元范围内所有 的ω值,进行算术平均得到ω'absolute(i)(i=1,2,3……,n)值;以及获取下看LADCP的m>absolute(λ)(λ=1,2,3……,m)值;按深度将ωa″bsolute(λ)线性插值到对应上看LADCP的>a'bsolute(i)(i=1,2,3……,n);然后,将同一深度单元内包含线>absolute(i)和上看LADCP的ω′absolute(i)进行算术平均,得到该深度单元>即得到n个深度单元的各个深度单 元的平均绝对流速
其中,如果m=n,则ω″absolute(λ)刚好一一线性插值到对应上看LADCP的各深>absolute(λ)值线性插>
进一步地,所述步骤S3中的生成统一公共深度轴下的数据表格,其具体过 程如下:
步骤1),首先将从海平面至海底深度从1m开始,并以间距为1m依次递增 生成公共深度轴作为第一列,用深度表示;
步骤2),定义第一组数据包括步骤002)获得的温度、步骤003)获得的 盐度、步骤004)获得的溶解氧、PH、叶绿素和步骤016)获得的绝对回波强度, 第二组数据包括步骤016)获得的绝对回波强度、步骤007)获得的绝对水平流 速东向分量和绝对水平流速北向分量、步骤010)获得的理查森数、步骤011) 获得的垂向扩散系数、步骤014)获得的绝对垂向流速和步骤017)获得的反演 生物量;
步骤3),获取每隔1m的各深度范围内的第一组数据的温度、盐度、溶解 氧、PH、叶绿素和绝对回波强度的测量值,将同一深度范围内的温度值进行算术 平均,得到温度的测量结果并放入表格对应深度处;将同一深度范围内的盐度值 进行算术平均,得到盐度的测量结果并放入表格对应深度处;将同一深度范围内 的溶解氧值进行算术平均,得到溶解氧的测量结果并放入表格对应深度处;将同 一深度范围内的溶PH值进行算术平均,得到PH的测量结果并放入表格对应深度 处;将同一深度范围内的叶绿素值进行算术平均,得到叶绿素的测量结果并放入 表格对应深度处;将同一深度范围内的绝对回波强度值进行算术平均,得到绝对 回波强度的测量结果并放入表格对应深度处;
步骤4),获取均以深度单元为4m进行测得的第二组数据,并在4的整数 倍的各深度上填写对应的第二组数据的各测量值,其余深度上用“NaN”来表示 该深度下无对应测量数据,最终,生成统一公共深度轴下的数据表格。
实现本发明目的之二的技术方案为:一种处理终端,其包括,
存储器,用于存储程序指令;
处理器,用于运行所述程序指令,以执行如下步骤:
步骤S1:获取包括压力、温度、盐度、电导率、溶解氧、PH、叶绿素、声 速、船载ADCP的声速、垂向速度、速度剪切、相对流速、相对回波强度和生物 量数据的原始数据,其中压力、温度、盐度、电导率和声速直接由CTD获得,溶 解氧、PH和叶绿素均通过水样观测获得,船载ADCP的声速由船载ADCP根据温 度探头获得,垂向速度、速度剪切、相对流速由LADCP测得,相对回波强度由 LADCP和船载ADCP获得,生物量数据通过生物拖网获得;
步骤S2:根据获得的原始数据进行数据处理,包括:
001)压力修正,修正公式为:压力真实值=原始压力值-甲板压力,其中原 始压力值为步骤S1中的压力,甲板压力为CTD回收后在甲板处所测得的压力;
002)温度修正,修正公式为:offset=b*residual/n,式中,*为乘法运算, b表示“从航次校正时间起至当前对温度修正之间的天数”,residual表示“航 次前与航次后的温度差”,n表示“航次前与航次后温度校正的间隔天数”;
003)盐度修正,首先对电导率数据进行校正,校正公式如①:
式中,n表示校正所用的水化学分析样本数,αi表示第i个水化学分析样>i表示第i个水化学分析样本由水化学分析实验获得的>
根据公式②获得电导率真实值:
电导率真实值=电导率读数*slope------②
电导率读数是由步骤S1中获得的电导率,获得电导率真实值后,根据海水 状态方程计算出盐度,该盐度即为修正后的盐度值;
004)溶解氧、PH和叶绿素修正,溶解氧的修正首先根据公式③计算出斜率 修正参数:
式中,slopeDO表示斜率修正参数,DO水样表示在某一深度的水化学样本通>CTD表示与DO水样同一深度由CTD获得的>水样为直接从步骤S1中获得的溶解氧原始数据中读取到同一深度下>
获得斜率修正参数后,按公式④计算:
溶解氧真实值=溶解氧读数*slopeDO------④
溶解氧读数是由CTD获得的从海面到海底的所有溶解氧的原始数据,溶解 氧真实值即是斜率修正后的溶解氧;
PH和叶绿素的修正同样按公式③和④进行计算,获得修正后的PH和叶绿 素;
005)声速校正,获得校正后的相对流速值,根据公式⑤进行计算:
Vcorrected=Vrelative*(Cctd/CADCP)------⑤
式中,Vrelative表示未校正的相对流速,由步骤S1中获得,CADCP表示船载>ctd是步骤S1中获得的声速,由步骤S1中获>corrected表示因声速引起误差而经过校正后的相对流速值;
006)低散射校正,获得校正后的速度剪切值,根据方程组⑥进行计算:
V修正=V0+ΔV------⑥
式中,EI表示相对回波强度,由步骤S1获得,k为比例系数,为常数,ΔV 表示修正值,V0表示LADCP获得的速度剪切,由步骤S1获得,V修正表示经过低>
007)绝对流速计算,首先通过剪切法计算出参考流速,计算公式如⑦:
式中,T表示LADCP下放的总时长,t表示时间,Vship表示船舶航速,Vrefer表示参考流速,
获得Vrefer后按公式⑧计算,获得绝对流速Vabsolute:
Vabsolute=Vrelative+Vrefer------⑧
式中,Vabsolute表示绝对流速;
008)双LADCP绝对流速平均:将下看LADCP按步骤007)计算获得绝对流 速,并按照深度单元进行平均,获得V′absolute(i)(i=1,2,3......,n);同时,将上看LADCP>absolute(λ)(λ=1,2,3......,m);将V″absolute(λ)(λ=1,2,3…...,m)线性插值到下看LADCP>其中,m表示上看LADCP的深度单元个数,n表示下 看LADCP的深度单元个数;
009)去潮:通过软件求取正压潮剖面,从步骤007)获得的绝对流速Vabsolute减去所述正压潮剖面,即获得去潮后的结果;
010)理查森数计算,获得理查森数Ri,按公式⑨进行计算:
式中,g表示重力加速度,m表示m米深度,n表示n米深度,ρm表示m米>n表示n米深度处的海水密度,um、vm分别表示m米深度>n、vn分别表示n米深度>
011)垂向扩散系数计算,获得垂向扩散系数,计算公式如公式⑩所示:
式中,v0取10-2m2/s,α取5,vb取10-4m2/s,kb取10-5m2/s,KT是垂向>
012)CTD深度数据低通滤波:采用截断频率为6Hz的低通滤波器对CTD深 度数据进行低通滤波,得到6Hz的CTD深度数据;
013)CTD垂向运动速度计算,获得CTD的垂向运动速度ωCTD,按公式
式中,z(t)表示012)中获得的6Hz的CTD深度数据,t表示时间;
014)绝对垂向速度计算,获得绝对垂向速度,按公式
式中,ωADCP表示垂向速度,由步骤S1中通过LADCP获得的垂向速度,ω表>
015)双LADCP绝对垂向速度平均,获得平均绝对垂向速度,即将包括上看 LADCP和下看LADCP两台设备中的下看LADCP按公式
016)绝对回波强度计算,获得绝对回波强度,首先根据公式
式中,B是船载ADCP盲区,为常数;P是宣称脉冲长度,为常数;W是深 度单元厚度,N是深度单元个数,θ是波束角,cave是船载ADCP某一深度单元的>
根据公式
式中,K1c、a、b和c2均是船载ADCP的出厂参数,均为常数,Vs是整个投放>1为船载ADCP发>
根据公式
其中,
P2=1-1.37×10-4z+6.2×10-9z2,P3=1-1.383×10-5z+4.9×10-10z2,
式中,PH、S、T和z为船载ADCP在第i个深度单元相应的海水PH、盐度、温度、 深度值,由步骤S1获得,ci表示第i个深度单元对应的声速,freq表示船载ADCP>
然后,根据公式
式中,K2是系统噪声因子,为常数,Ks是与船载ADCP频率有关的系统固>x是当前船载ADCP探头所处位置的海水温度,由步骤S1获得,>a为步骤S1中LADCP测量获得的相对回波强度,Er为有效剖面末端的RSSI值,>V是绝对回波强度;
017)生物量反演,获得反演生物量,选取相同数量的若干个生物量数据和 绝对回波强度数据进行数据拟合,根据拟合公式
式中,DW为步骤S1中获得的生物量数据,
经过公式
步骤S3:获得上述经过处理后的各个数据后,对经过步骤S2处理获得的数 据进行数据融合,生成统一公共深度轴下的数据表格。
6.根据权利要求5所述的处理终端,其特征在于:所述步骤008)中的按照 深度单元平均,其具体过程为:
步骤a),以Z'深值为间隔将从上看LADCP初始位置开始至海底进行划分,得>深值为间隔将从下看LADCP初始位置开始直至海底进行划分,得>
其中,m和n均为正整数,Z'深和Z”深分别为上看LADCP和下看LADCP采集数>
步骤b),首先,获取上看LADCP的n个深度单元的每个深度单元范围内所有 的Vabsolute值,进行算术平均得到V′absolute(i)(i=1,2,3…...,n)值;以及获取下看LADCP>absolute值,进行算术平均得到>absolute(λ)(λ=1,2,3......,m)值;按深度将V″absolute(λ)线性插值到对应上看LADCP的>absolute(i)(i=1,2,3......,n);然后,将同一深度单元内包含线>absolute(i)和上看LADCP的Va′bsolute(i)进行算术平均,得到该深度单元>最终得到n个深度单元的各个深度 单元的平均绝对流速
其中,如果m=n,则V″absolute(λ)刚好一一线性插值到对应上看LADCP的各深>absolute(λ)>
7.根据权利要求5所述的处理终端,其特征在于:所述步骤015)中的按照 深度单元平均,其具体过程为:
步骤a),以Z'深值为间隔将从上看LADCP初始位置开始至海底进行划分,得>深值为间隔将从下看LADCP初始位置开始直至海底进行划分,得>
其中,m和n均为正整数,Z'深和Z”深分别为上看LADCP和下看LADCP采集数>
步骤b),首先,获取上看LADCP的n个深度单元的每个深度单元范围内所有 的ω值,进行算术平均得到ω'absolute(i)(i=1,2,3......,n)值;以及获取下看LADCP的m>absolute(λ)(λ=1,2,3......,m)值;按深度将ω″absolute(λ)线性插值到对应上看LADCP的>absolute(i)(i=1,2,3......,n);然后,将同一深度单元内包含线>absolute(i)和上看LADCP的ω′absolute(i)进行算术平均,得到该深度单元>即得到n个深度单元的各个深度单 元的平均绝对流速
其中,如果m=n,则ω″absolute(λ)刚好一一线性插值到对应上看LADCP的各深>absolute(λ)值线>
8.根据权利要求5所述的处理终端,其特征在于:所述步骤S3中的生成统 一公共深度轴下的数据表格,其具体过程如下:
步骤1),首先将从海平面至海底深度从1m开始,并以间距为1m依次递增 生成公共深度轴作为第一列,用深度表示;
步骤2),定义第一组数据包括步骤002)获得的温度、步骤003)获得的 盐度、步骤004)获得的溶解氧、PH、叶绿素和步骤016)获得的绝对回波强度, 第二组数据包括步骤016)获得的绝对回波强度、步骤007)获得的绝对水平流 速东向分量和绝对水平流速北向分量、步骤010)获得的理查森数、步骤011) 获得的垂向扩散系数、步骤014)获得的绝对垂向流速和步骤017)获得的反演 生物量;
步骤3),获取每隔1m的各深度范围内的第一组数据的温度、盐度、溶解 氧、PH、叶绿素和绝对回波强度的测量值,将同一深度范围内的温度值进行算术 平均,得到温度的测量结果并放入表格对应深度处;将同一深度范围内的盐度值 进行算术平均,得到盐度的测量结果并放入表格对应深度处;将同一深度范围内 的溶解氧值进行算术平均,得到溶解氧的测量结果并放入表格对应深度处;将同 一深度范围内的溶PH值进行算术平均,得到PH的测量结果并放入表格对应深度 处;将同一深度范围内的叶绿素值进行算术平均,得到叶绿素的测量结果并放入 表格对应深度处;将同一深度范围内的绝对回波强度值进行算术平均,得到绝对 回波强度的测量结果并放入表格对应深度处;
步骤4),获取均以深度单元为4m进行测得的第二组数据,并在4的整数 倍的各深度上填写对应的第二组数据的各测量值,其余深度上用“NaN”来表示 该深度下无对应测量数据,最终,生成统一公共深度轴下的数据表格。
本发明的有益效果为:本发明对环境要素信息参数进行修正,提高了各个参 数的正确性,进而提高了水合物环境基线调查与评价的正确性,并将所有环境要 素信息融合到统一深度公共深度轴,交互性强。
附图说明
图1为本发明较佳实施例的流程图;
图2为本发明一种处理终端的结构示意图。
具体实施方式
下面,结合附图以及具体实施方式,对本发明做进一步描述:
如图1和图2所示,一种深海全水柱多环境要素信息融合方法,其包括以下 步骤:
步骤S1:获取包括压力、温度、盐度、电导率、溶解氧、PH、叶绿素、声 速、船载ADCP的声速、垂向速度、速度剪切、相对流速、相对回波强度和生物 量数据的原始数据,其中压力、温度、盐度、电导率和声速直接由CTD(中文名 为温盐深仪)设备获得,其中,溶解氧、PH和叶绿素均通过水样观测获得,船 载ADCP的声速由船载ADCP根据温度探头获得的声速值,垂向速度、速度剪切、 相对流速由LADCP测得,相对回波强度由LADCP和船载ADCP获得,生物量数据 通过生物拖网获得,原始数据及获得的途径如下表所示:
本实施例中,LADCP采用由RDI(Teledyne RD Instruments,简称为RDI) 公司生产的WHS●ADCP,且本实施例中某些步骤中会采用两台LADCP,定义为下 看LADCP和上看LADCP,船载ADCP采用RDI公司生产的OS●ADCP;
上述所有的原始数据均为从海面到海底测出的数据,是一定的海洋深度进行 测量获得的数据,该数据包括了深度信息,如温度值,表示在某个海洋深度对应 的温度;
步骤S2:根据获得的原始数据进行数据处理,包括:
001)压力修正,修正公式为:压力真实值=原始压力值-甲板压力,其中原 始压力值为步骤S1中的压力,甲板压力为CTD回收后在甲板处所测得的压力;
002)温度修正,修正公式为:offset=b*residual/n,式中,*为乘法运算, b表示“从航次校准时间起至当前对温度修正之间的天数”,residual表示“航 次前与航次后的温度差”,n表示“航次前与航次后温度校准的间隔天数”,offset 表示修正后的温度值,比如航次前在4月1日进行一次温度校准,航次回来后再 7月1日进行了温度校准,并对此次航行中的5月1日的温度数据进行校准,则 n=90,b=30,residual=7月1日的校准后的温度-4月1日校准后的温度;
003)盐度修正,首先对电导率数据进行校准,校准公式如下:
式中,n表示校准所用的水化学分析样本数,αi表示第i个水化学分析样>i表示第i个水化学分析样本由水化学分析实验获得的>
根据修正后的电导率获得电导率真实值,计算公式为:电导率真实值=电导 率读数*slope,电导率读数表示从海面开始至海底整个测量过程通过CTD所测 得的所有电导率数据,也即是由步骤S1中获得的电导率,获得电导率真实值后, 根据海水状态方程计算出盐度,该盐度即为修正后的盐度值;
004)溶解氧、PH和叶绿素修正,以溶解氧为例,首先计算斜率修正参数, 计算公式为:
式中,slopeDO表示斜率修正参数,DO水样表示在某一深度的水化学样本通>CTD表示与DO水样同一深度由CTD获得的>水样可直接从步骤S1中获得的溶解氧原始数据中读取到同一深度下>水样和DOCTD是在同一深度下获得的溶解氧即可;
获得斜率修正参数后,对CTD设备测量的溶解氧剖面进行斜率修正,输出 新的溶解氧剖面,其计算公式为:溶解氧真实值=溶解氧读数*slopeDO,溶解氧>
005)声速校正,船载ADCP使用固定的盐度值和自身低准确度的温度数据 计算声速,进而引起相对流速计算误差,因此需要对由声速引起相对流速产生误 差而进行校正,校正公式为:
Vcorrected=Vuncorrected*(Cctd/CADCP)
式中,Vuncorrected表示未校正的相对流速值,由步骤S1中获得,CADCP表示船>ctd是步骤S1中获得的声速,Vcorrected表示>
006)低散射校正,LADCP在深海低散射环境下,回波显著减少,导致所测 速度剪切出现偏差,因而需要进行因低散射引起的速度剪切误差的校正,校正公 式为:
V修正=V0+ΔV
式中,EI表示相对回波强度,由步骤S1获得,k为比例系数,为人为选择 的一个常数,在本实施例中,取值为-0.12,ΔV表示修正值,V0表示LADCP获>修正表示经因低散>
007)绝对流速计算,常规的LADCP所测的相对速度为LADCP与水体之间的 相对速度,并非绝对速度,因此需要计算出绝对流速,可通过剪切法获得:
剪切法,计算出船舶漂移距离,确定参考流速,对相同深度单元的流速垂 向剪切进行平均,得到全水柱的流速垂向剪切,再从海底向海面进行垂向积分, 获得相对流速剖面,并且计算船舶漂移距离,确定参考流速,计算公式如下:
式中,T表示LADCP下放的总时长,t表示时间,Vship表示船舶航速,Vrelative为步骤S1获得的相对流速,Vrefer表示参考流速,
获得Vrefer后按公式(1)计算,获得绝对流速Vabsolute:
Vabsolute=Vrelative+Vrefer(1)
绝对流速Vabsolute是矢量,通过对Vabsolute进行正交分解,可以分解出绝对水平>
008)双LADCP绝对流速平均:获得平均绝对流速 步骤a),获取上看和下看LADCP有相同的深度单元Z深值,并以Z深值为间>深值为间> 其中,m和n均为正整数,由于上看LADCP和下看LADCP的初审位置并不一 定相同,因此m值不一定等于n值; 步骤b),首先,获取上看LADCP的n个深度单元的每个深度单元范围内所有 的V′absolute(i)(i=1,2,3……,n)值,以及获取下看LADCP的m个深度单元的每个深度>absolute(λ)(λ=1,2,3......,m)值,将V″absolute(λ)线性插值到对应上>absolute(λ)>absolute(i)进行算术平均,得到该深度单元内的平均绝对流速 其中,如果m=n,则V″absolute(λ)刚好一一线性插值到对应上看LADCP的各深>absolute(λ)值均线性插入至上看LADCP的第n个深度单元,如果m<n,则上看>absolute(λ)插入, 其中,深度单元的值由LADCP自身确定,上看LADCP有与下看LADCP相同 的深度单元,每个LADCP都会有固定的深度单元,不同的LADCP确定的深度单元 可能会有所不同,这里的深度单元是指LADCP由当前位置向下所能测量最大范 围,比如深度单元为4m,表示LADCP由当前位置向下4m范围内能够测得到的数 据;正如前文所述的,随着LADCP从海面向海底移动的过程中,LADCP得到包括 第一深度单元、第二深度单元等等多个深度单元,每个深度单元的值是相同的, 这样LADCP在每个深度单元内能够获得多个不同值,将该深度单元内所有不同值 进行算术平均,算术平均后获得的值即代表该深度单元内的平均绝对流速; 009)去潮,由于LADCP得到的绝对流速Vabsolute中包含潮,需要扣除,即去>absolute扣除正压潮剖面,即获得去潮后的结果; 010)理查森数(Richardson number)计算,按公式(2)进行计算,获得 理查森数Ri: 式中,g表示重力加速度,取值为9.8,ρ为海水密度,z表示海洋深度, u和v分别表示Vabsolute的绝对水平流速东向分量和绝对水平流速北向分量。 公式(2)通常意义上的理解,是对针对有明确函数表达式(显式、隐函数、 反函数、参数方程等)而运用求导法则进行计算的公式,但实际问题中,往往只 是得到自变量和因变量的一些离散的数据;而为了便于计算,在本实施例中,对 公式(2)变换成如下公式进行计算: 式中,m表示m米深度,n表示n米深度,ρm表示m米深度处的海水密度,>n表示n米深度处的海水密度,um、vm分别表示m米深度处的水平流速东向分>n、vn分别表示n米深度处的绝对水平流速东向分量和> 011)垂向扩散系数计算,计算公式如方程(3)所示: 其中,v0取10-2m2/s,α取5,Ri为010)获得的理查森数,vb取10-4m2/s,>b取10-5m2/s,KT是垂向扩散系数; 012)CTD深度数据低通滤波,通过低通滤波器去除高频数据,具体过程为: 采用截断频率为6Hz的低通滤波器对CTD深度数据进行低通滤波,得到6Hz的 CTD深度数据; 013)CTD垂向运动速度计算,按公式(4)进行计算,获得CTD的垂向运动 速度ωCTD: 式中,z(t)表示012)中获得的6Hz的CTD深度数据,t表示时间; 014)绝对垂向速度计算,计算公式如下: ω=ωADCP-ωCTD 式中,ωADCP表示垂向速度,由步骤S1中通过LADCP获得的垂向速度,ω表> 015)双LADCP绝对垂向速度ω数据平均,获得平均绝对垂向速度,将包括 上看LADCP和下看LADCP两台设备中的下看LADCP按步骤014)计算获得绝对垂 向速度,并按照深度单元进行平均,获得ω′absolute(i)(i=1,2,3......,n),同时,将上>absolute(λ)(λ=1,2,3......,m),并将ω″absolute(λ)(λ=1,2,3......,m)线性插值到下看>具体计算过程与步骤008)类似,在这里就不赘述了; 016)绝对回波强度计算,将LADCP与船载ADCP获得的相对回波强度数据转 化为绝对回波强度,首先计算声吸收系数e、深度单元距离R和发射能量K1,然> 计算深度单元距离R的公式如下: 式中,B是船载ADCP盲区,人为选择的一个常数,P是宣称脉冲长度,为 常数,W是深度单元厚度,N是深度单元个数,θ是波束角,CADCP是船载ADCP>ave是从船载ADCP探头到某一深度单> 计算发射能量K1,计算公式如下: K1=0.34K1C(Vs·a-b)/c2 式中,K1c、a、b和c2均是船载ADCP的出厂参数,在本实施例中,其取值分>s是整个投放过程中高电压值的平均值,高>1为船载ADCP发射能量; 计算声吸收系数e,其计算公式如下: 其中, P2=1-1.37×10-4z+6.2×10-9z2,P3=1-1.383×10-5z+4.9×10-10z2, 式中,PH、S、T和z为船载ADCP在第i个深度单元相应的海水PH、盐度、温度、 深度值,从步骤S1获得,ci表示第i个深度单元对应的声速;freq表示船载ADCP> 然后计算绝对回波强度,其计算公式如下: 式中,K2是系统噪声因子,为常数,Ks是与船载ADCP频率有关的系统固>2和Ks的取值分别为3.6和4.17*105,Tx是当前>a为步骤S1中LADCP测>r为有效剖面末端的RSSI值,为常数,本实施例中,>V是绝对回波强度; 017)生物量反演,具体过程如下:选取相同数量的若干个生物量数据和绝 对回波强度,本实施例,选取50个生物量数据和绝对回波强度数据进行数据拟 合,根据拟合公式计算出回归系数项a1和b1,其拟合公式如下: 式中,DW为步骤S1中获得的生物量数据,也由于DW和SV均与时间和深>V为相同时间下的相同深度的数据,>1和b1后,根据公式(5)得到生物量数据DW的计算公式(6): 经过上述计算,即可选取若干个数据经过拟合后得到公式(6),再根据绝 对回波强度SV代入公式(6)即可计算出从海面到海底整个任意深度的生物量数> 步骤S3:获得上述经过处理后的各个数据后,对经过步骤S2处理获得的数 据进行数据融合,生成统一公共深度轴下的数据表格,其具体过程如下: 第一列为深度,从1m开始以间距为1m依次递增生成公共深度轴,直至到海 底,第一列的深度表示距离海平面的海洋深度,比如2m处表示距离海平面为2m 的海洋深度,第一列的最后一个值表示海底到海平面的深度; 定义第一组数据包括步骤002)获得的温度、步骤003)获得的盐度、步骤 004)获得的溶解氧、PH、叶绿素和步骤016)获得的绝对回波强度,第二组数 据包括步骤016)获得的绝对回波强度、步骤007)获得的绝对水平流速东向分 量和绝对水平流速北向分量、步骤010)获得的理查森数、步骤011)获得的垂 向扩散系数、步骤014)获得的绝对垂向流速和步骤017)获得的反演生物量; 需要说明的是,在实际调查测量过程中,CTD、LADCP和船载ADCP是从海面 到海底进行测量,所获得的所有原始数据均包含了从海面到海底之间对应的深度 和时间信息,即任何一个数据包括了深度和时间的信息,比如温度值,该温度值 表示对应在某个时刻下的某个深度的温度;其中,第一组数据均通过CTD和水样 观测测得原始数据进行处理后的结果,本实施例中,CTD测量的深度单元为1m, 即是指CTD能测得当前CTD所处位置1m以内的数据,这样在一个深度单元内, CTD测得的第一组数据中的任意一个参数的个数与CTD的测量频率和其下降速度 有关,假如CTD的测量频率为24次每秒,下降速度为1m/s,则第一组数据的各 参数能在1m的深度间隔内各测得24个数据,再分别对单个深度单元附近的24 个数据进行算数平均,从而在1m、2m、3m……等各深度均只有一个数据;本实 施例中,船载ADCP和LADCP测量的深度单元为4m,假如LADCP的测量频率同样 为一次每秒,下降速度为1m/s,则第二组数据均以4m的整数倍的深度处各测得 一个数据,即在4m、8m、12m……等各深度均只有一个数据; 因此第一组数据的深度轴和第二组数据的深度轴并不一致,为了统一进行对 比,将第一组数据和第二组数据统一融合到统一公共深度轴上,此外,无论是 CTD测量的深度单元还是船载ADCP和LADCP测量的深度单元均由其各自的设备 自身确定,可直接从对应的设备获得该参数; 由于第二组数据是以深度单元为4m进行测得的数据,因此只在以4m的深度 间隔各深度处有一个数据,即在4m、8m、12m……等各深度均只有一个数据,其 余深度上用“NaN”来表示该深度下无对应测量数据,因此,第二组数据的分辨 率为4m; 第一组数据,在0.5-1.5m的深度范围内可能测得多个值,将这所有测得的 值进行算术平均,算术平均后的结果即表示1m深度处的测量结果,比如0.5-1.5m 的深度范围内测得n个温度值xi(i=1,2,3……,n),则在1m深度处的 最终获得的第一组数据和第二组数据形成的表格如表一所示: 表一 对应于表一可以生成一张质控标记的表格,用“0”表示数据质量优良,“1” 表示数据质量未知,“2”表示数据质量差,初始默认值均为1,工作人员可以 使用这一表格依据经验和实际需求对各深度的环境要素参数进行标记,当判断某 个深度下对应的某个参数数据质量优良时,该处标记为“0”即可,当判断某个 深度下对应的某个参数数据质量差时,该处标记为“2”即可,这样有一个直观 的观察,表二如下所示: 表二 另外,本发明还涉及一种实现以上方法的实体装置的处理终端100,其包括, 存储器101,用于存储程序指令; 处理器102,用于运行所述程序指令,以执行如下步骤: 步骤S1:获取包括压力、温度、盐度、电导率、溶解氧、PH、叶绿素、声 速、船载ADCP的声速、垂向速度、速度剪切、相对流速、相对回波强度和生物 量数据的原始数据,其中压力、温度、盐度、电导率和声速直接由CTD获得,溶 解氧、PH和叶绿素均通过水样观测获得,船载ADCP的声速由船载ADCP根据温 度探头获得,垂向速度、速度剪切、相对流速由LADCP测得,相对回波强度由 LADCP和船载ADCP获得,生物量数据通过生物拖网获得; 步骤S2:根据获得的原始数据进行数据处理,包括: 001)压力修正,修正公式为:压力真实值=原始压力值-甲板压力,其中原 始压力值为步骤S1中的压力,甲板压力为CTD回收后在甲板处所测得的压力; 002)温度修正,修正公式为:offset=b*residual/n,式中,*为乘法运算, b表示“从航次校正时间起至当前对温度修正之间的天数”,residual表示“航 次前与航次后的温度差”,n表示“航次前与航次后温度校正的间隔天数”; 003)盐度修正,首先对电导率数据进行校正,校正公式如①: 式中,n表示校正所用的水化学分析样本数,αi表示第i个水化学分析样>i表示第i个水化学分析样本由水化学分析实验获得的> 根据公式②获得电导率真实值: 电导率真实值=电导率读数*slope------② 电导率读数是由步骤S1中获得的电导率,获得电导率真实值后,根据海水 状态方程计算出盐度,该盐度即为修正后的盐度值; 004)溶解氧、PH和叶绿素修正,溶解氧的修正首先根据公式③计算出斜率 修正参数: 式中,slopeDO表示斜率修正参数,DO水样表示在某一深度的水化学样本通>CTD表示与DO水样同一深度由CTD获得的>水样为直接从步骤S1中获得的溶解氧原始数据中读取到同一深度下> 溶解氧真实值=溶解氧读数*slopeDO------④ 溶解氧读数是由CTD获得的从海面到海底的所有溶解氧的原始数据,溶解 氧真实值即是斜率修正后的溶解氧; PH和叶绿素的修正同样按公式③和④进行计算,获得修正后的PH和叶绿 素; 005)声速校正,获得校正后的相对流速值,根据公式⑤进行计算: Vcorrected=Vrelative*(Cctd/CADCP)------⑤ 式中,Vrelative表示未校正的相对流速,由步骤S1中获得,CADCP表示船载>ctd是步骤S1中获得的声速,由步骤S1中获>corrected表示因声速引起误差而经过校正后的相对流速值; 006)低散射校正,获得校正后的速度剪切值,根据方程组⑥进行计算: V修正=V0+ΔV------⑥ 式中,EI表示相对回波强度,由步骤S1获得,k为比例系数,为常数,ΔV 表示修正值,V0表示LADCP获得的速度剪切,由步骤S1获得,V修正表示经因低> 007)绝对流速计算,首先通过剪切法计算出参考流速,计算公式如⑦: 式中,T表示LADCP下放的总时长,t表示时间,Vship表示船舶航速,Vrefer表示参考流速, 获得Vrefer后按公式⑧计算,获得绝对流速Vabsolute: Vabsolute=Vrelative+Vrefer------⑧ 式中,Vabsolute表示绝对流速; 008)双LADCP绝对流速平均:获得平均绝对流速 009)去潮,通过软件求取正压潮剖面,从步骤007)获得的绝对流速Vabsolute减去所述正压潮剖面,即获得去潮后的结果; 010)理查森数计算,获得理查森数Ri,按公式⑨进行计算: 式中,g表示重力加速度,m表示m米深度,n表示n米深度,ρm表示m米>n表示n米深度处的海水密度,um表示m米深度处的绝对>n表示n米深度处的绝对水平流速北向分量; 011)垂向扩散系数计算,获得垂向扩散系数,计算公式如公式 式中,v0取10-2m2/s,α取5,vb取10-4m2/s,kb取10-5m2/s,KT是垂向> 012)CTD深度数据低通滤波,采用截断频率为6Hz的低通滤波器对CTD深 度数据进行低通滤波,得到6Hz的CTD深度数据; 013)CTD垂向运动速度计算,获得CTD的垂向运动速度ωCTD,按公式 式中,z(t)表示012)中获得的6Hz的CTD深度数据,t表示时间; 014)绝对垂向速度计算,获得绝对垂向速度,按公式 式中,ωADCP表示垂向速度,由步骤S1中通过LADCP获得的垂向速度,ω表> 015)双LADCP绝对垂向速度平均,获得平均绝对垂向速度,将包括上看LADCP 和下看LADCP两台设备中的下看LADCP按步骤014)计算获得的ω′,并将ω′按 照深度单元进行平均,同时,将上看LADCP按步骤014)计算获得的ω″,并将ω″ 线性插值到下看LADCP的深度单元轴上,然后再对两者求算术平均,获得平均绝 对垂向速度 016)绝对回波强度计算,获得绝对回波强度,首先根据公式 式中,B是船载ADCP盲区,为常数,P是宣称脉冲长度,为常数,W是深 度单元厚度,N是深度单元个数,θ是波束角,cave是船载ADCP某一深度单元的> 根据公式 式中,K1c、a、b和c2均是船载ADCP的出厂参数,均为常数,Vs是整个投放>1为船载ADCP发> 根据公式 其中, P2=1-1.37×10-4z+6.2×10-9z2,P3=1-1.383×10-5z+4.9×10-10z2, 式中,PH、S、T和z为船载ADCP在第i个深度单元相应的海水PH、盐度、温度、 深度值,由步骤S1获得,ci表示第i个深度单元对应的声速,freq表示船载ADCP> 然后,根据公式 式中,K2是系统噪声因子,为常数,Ks是与船载ADCP频率有关的系统固>x是当前船载ADCP探头所处位置的海水温度,由步骤S1获得,>a为步骤S1中LADCP测量获得的相对回波强度,Er为有效剖面末端的RSSI值,>V是绝对回波强度; 017)生物量反演,获得反演生物量,选取相同数量的若干个生物量数据和 绝对回波强度进行数据拟合,根据拟合公式 式中,DW为步骤S1中获得的生物量数据, 经过公式 步骤S3:获得上述经过处理后的各个数据后,对经过步骤S2处理获得的数 据进行数据融合,生成统一公共深度轴下的数据表格。 进一步地,其具体过程为: 步骤a)以Z'深值为间隔将从上看LADCP初始位置开始至海底进行划分,得到>深值为间隔将从下看LADCP初始位置开始直至海底进行划分,> 其中,m和n均为正整数,Z'深和Z”深分别为上看LADCP和下看LADCP采集数> 步骤b),首先,获取上看LADCP的n个深度单元的每个深度单元范围内所有 的Vabsolute值,进行算术平均得到V′absolute(i)(i=1,2,3......,n)值,以及获取下看LADCP>absolute值,进行算术平均得到>absolute(λ)(λ=1,2,3......,m)值,按深度单元顺序将V″absolute(λ)线性插值到对应上看>absolute(λ)和>absolute(i)进行算术平均,得到该深度单元内的平均绝对流速>最终得到n个深度单元的各个深度单元的平均绝对流 速 其中,如果m=n,则V″absolute(λ)刚好一一线性插值到对应上看LADCP的各深>absolute(λ)值线性>absolute(λ)值,并将该两个深度单元的V″absolute(λ)值将线> 进一步地,所述步骤015)中的按照深度单元平均,其具体过程为: 步骤a)以Z'深值为间隔将从上看LADCP初始位置开始至海底进行划分,得到>深值为间隔将从下看LADCP初始位置开始直至海底进行划分,> 其中,m和n均为正整数,Z'深和Z”深分别为上看LADCP和下看LADCP采集数> 步骤b),首先,获取上看LADCP的n个深度单元的每个深度单元范围内所有 的ω值,进行算术平均得到ω′absolute(i)(i=1,2,3……,n)值,以及获取下看LADCP的>absolute(λ)(λ=1,2,3……,m)值,按深度将ω″absolute(λ)线性插值到对应上看LADCP的>absolute(λ)和上看>absolute(i)进行算术平均,得到该深度单元内的平均绝对流速>最终得到n个深度单元的各个深度单元的平均绝对 流速 其中,如果m=n,则ω″absolute(λ)刚好一一线性插值到对应上看LADCP的各深>absolute(λ)值线性> 进一步地,所述步骤S3中的生成统一公共深度轴下的数据表格,其具体过 程如下: 步骤1),首先将从海平面至海底深度从1m开始,并以间距为1m依次递增 生成公共深度轴作为第一列,用深度表示; 步骤2),定义第一组数据包括步骤002)获得的温度、步骤003)获得的 盐度、步骤004)获得的溶解氧、PH、叶绿素和步骤016)获得的绝对回波强度; 第二组数据包括步骤016)获得的绝对回波强度、步骤007)获得的绝对水平流 速东向分量和绝对水平流速北向分量、步骤010)获得的理查森数、步骤011) 获得的垂向扩散系数、步骤014)获得的绝对垂向流速和步骤017)获得的反演 生物量; 步骤3),获取每隔1m的各深度范围内的第一组数据的温度、盐度、溶解 氧、PH、叶绿素和绝对回波强度的测量值,将同一深度范围内的温度值进行算术 平均,得到温度的测量结果并放入表格对应深度处,将同一深度范围内的盐度值 进行算术平均,得到盐度的测量结果并放入表格对应深度处,将同一深度范围内 的溶解氧值进行算术平均,得到溶解氧的测量结果并放入表格对应深度处,将同 一深度范围内的溶PH值进行算术平均,得到PH的测量结果并放入表格对应深度 处,将同一深度范围内的叶绿素值进行算术平均,得到叶绿素的测量结果并放入 表格对应深度处,将同一深度范围内的绝对回波强度值进行算术平均,得到绝对 回波强度的测量结果并放入表格对应深度处; 步骤4),获取均以深度单元为4m进行测得的第二组数据,并在4的整数 倍的各深度上填写对应的第二组数据的各测量值,其余深度上用“NaN”来表示 该深度下无对应测量数据,最终,生成统一公共深度轴下的数据表格。 对于本领域的技术人员来说,可根据以上描述的技术方案以及构思,做出其 它各种相应的改变以及变形,而所有的这些改变以及变形都应该属于本发明权利 要求的保护范围之内。
机译: 一种学习方法和学习设备,可以通过使用每个对象的深度预测信息和由V2X信息融合技术获取的每个对象的广泛预测信息来重建3D空间来更新HD映射的学习设备,以及用于由V2X信息融合技术获取的每个对象的类信息。测试方法和测试设备
机译: 提供一种用于通过使用竞争性计算和信息融合来支持车辆的自动驾驶的自动驾驶服务平台的方法,以及使用该方法的服务器。通过竞争性计算和信息融合来支持车辆的自动驾驶,并使用相同的服务器来管理}
机译: 一种通信系统,在发送方和接收方的多个终端之间传输邮件要素和用于邮件要素的通信方法,并使用和扩展Internet