一种深海全水柱多环境要素信息融合方法及处理终端技术

本发明专利技术涉及一种深海全水柱多环境要素信息融合方法及处理终端，所述方法包括以下步骤：步骤S1：获取包括压力、温度、盐度、电导率、溶解氧、PH、叶绿素、声速、船载ADCP的声速、垂向速度、速度剪切、相对流速、相对回波强度和生物量数据的原始数据；步骤S2：根据获得的原始数据进行数据处理；步骤S3：对经过步骤S2处理获得的数据进行数据融合，生成统一公共深度轴下的数据表格。本发明专利技术对环境要素信息进行修正，提高各个参数的正确性和垂向分辨率；组合不同设备采集的数据，处理生成新的环境要素，延展了调查数据的科研价值；并将所有环境要素信息融合到统一深度公共深度轴，交互性强。

【技术实现步骤摘要】
一种深海全水柱多环境要素信息融合方法及处理终端
本专利技术涉及海洋调查
，具体是一种深海全水柱多环境要素信息融合方法及处理终端。
技术介绍
深海环境基线调查需要同时获取多种有效、高分辨率的水柱环境要素信息，提高对海洋环境特征的辨识能力。而目前采用温盐深测量、走航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个水化学分析样本由CTD测出的电导率，βi表示第i个水化学分析样本由水化学分析实验获得的电导率，slope表示修正后的电导率，根据公式②获得电导率真实值：电导率真实值＝电导率读数*slope------②电导率读数是由步骤S1中获得的电导率，获得电导率真实值后，根据海水状态方程计算出盐度，该盐度即为修正后的盐度值；004)溶解氧、PH和叶绿素修正，溶解氧的修正首先根据公式③计算出斜率修正参数：式中，slopeDO表示斜率修正参数，DO水样表示在某一深度的水化学样本通过水化学实验分析获得的溶解氧值，DOCTD表示与DO水样同一深度由CTD获得的溶解氧，DO水样为直接从步骤S1中获得的溶解氧原始数据中读取到同一深度下的溶解氧数值；获得斜率修正参数后，按公式④计算：溶解氧真实值＝溶解氧读数*slopeDO------④溶解氧读数是由CTD获得的从海面到海底的所有溶解氧的原始数据，溶解氧真实值即是斜率修正后的溶解氧；PH和叶绿素的修正同样按公式③和④进行计算，获得修正后的PH和叶绿素；005)声速校正，获得校正后的相对流速值，根据公式⑤进行计算：Vcorrected＝Vrelative*(Cctd/CADCP)------⑤式中，Vrelative表示未校正的相对流速，由步骤S1中获得，CADCP表示船载ADCP的声速，由步骤S1中获得，Cctd是步骤S1中获得的声速，由步骤S1中获得，Vcorrected表示因声速引起误差而经过校正后的相对流速值；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按步骤007)计算绝对流速，并按照深度单元进行平均，获得V″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米深度处的绝对水平流速东向分量和绝对水平流速北向分量，un、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按公式获得绝对垂向速度，并按照深度单元进行平均，获得ω′absolute(i)(i＝1,2,3......,n)；同时，将上看LADCP按公式获得绝对垂向速度，并按照深度单元进行平均，获得ω″absolute(λ)(λ＝1,2,3......,m)；将ω″absolute(λ)(λ＝1,2,3......,m)线性插值到下看LADCP的深度单元轴上，然后再对两者求算术平均，获得平均绝对垂向速度m表示上看LADCP的深度单元个数，n表示下看LADCP的深度单元个数；016)绝对回波强度计算，获得绝对回波强度，首先根据公式计算深度单元距离R：式中，B是船载ADCP盲区，为常数；P是宣称脉冲长度，为常数；W是深度单元厚度，N是深度单元个数，θ是波束角，cave是船载ADCP某一深度单元的观测范围内的平均声速；根据公式计算发射本文档来自技高网...

【技术保护点】
1.一种深海全水柱多环境要素信息融合方法，其特征在于：其包括以下步骤：步骤S1：获取包括压力、温度、盐度、电导率、溶解氧、PH、叶绿素、声速、船载ADCP的声速、垂向速度、速度剪切、相对流速、相对回波强度和生物量数据的原始数据，其中压力、温度、盐度、电导率和声速直接由CTD获得，溶解氧、PH和叶绿素均通过水样观测获得，船载ADCP的声速由船载ADCP根据温度探头获得，垂向速度、速度剪切、相对流速由LADCP测得，相对回波强度由LADCP和船载ADCP获得，生物量数据通过生物拖网获得；步骤S2：根据获得的原始数据进行数据处理，包括：001)压力修正，修正公式为：压力真实值＝原始压力值‑甲板压力，其中原始压力值为步骤S1中的压力，甲板压力为CTD回收后在甲板处所测得的压力；002)温度修正，修正公式为：offset＝b*residual/n，式中，*为乘法运算，b表示“从航次校正时间起至当前对温度修正之间的天数”，residual表示“航次前与航次后的温度差”，n表示“航次前与航次后温度校正的间隔天数”；003)盐度修正，首先对电导率数据进行校正，校正公式如①：

【技术特征摘要】
1.一种深海全水柱多环境要素信息融合方法，其特征在于：其包括以下步骤：步骤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个水化学分析样本由CTD测出的电导率，βi表示第i个水化学分析样本由水化学分析实验获得的电导率，slope表示修正后的电导率，根据公式②获得电导率真实值：电导率真实值＝电导率读数*slope------②电导率读数是由步骤S1中获得的电导率，获得电导率真实值后，根据海水状态方程计算出盐度，该盐度即为修正后的盐度值；004)溶解氧、PH和叶绿素修正，溶解氧的修正首先根据公式③计算出斜率修正参数：式中，slopeDO表示斜率修正参数，DO水样表示在某一深度的水化学样本通过水化学实验分析获得的溶解氧值，DOCTD表示与DO水样同一深度由CTD获得的溶解氧，DO水样为直接从步骤S1中获得的溶解氧原始数据中读取到同一深度下的溶解氧数值；获得斜率修正参数后，按公式④计算：溶解氧真实值＝溶解氧读数*slopeDO------④溶解氧读数是由CTD获得的从海面到海底的所有溶解氧的原始数据，溶解氧真实值即是斜率修正后的溶解氧；PH和叶绿素的修正同样按公式③和④进行计算，获得修正后的PH和叶绿素；005)声速校正，获得校正后的相对流速值，根据公式⑤进行计算：Vcorrected＝Vrelative*(Cctd/CADCP)------⑤式中，Vrelative表示未校正的相对流速，由步骤S1中获得，CADCP表示船载ADCP的声速，由步骤S1中获得，Cctd是步骤S1中获得的声速，由步骤S1中获得，Vcorrected表示因声速引起误差而经过校正后的相对流速值；006)低散射校正，获得校正后的速度剪切值，根据方程组⑥进行计算：式中，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按步骤007)计算绝对流速，并按照深度单元进行平均，获得V″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米深度处的绝对水平流速东向分量和绝对水平流速北向分量，un、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按公式获得绝对垂向速度，并按照深度单元进行平均，获得ω′absolute(i)(i＝1,2,3……,n)；同时，将上看LADCP按公式获得绝对垂向速度，并按照深度单元进行平均，获得ω″absolute(λ)(λ＝1,2,3……,m)；将ω″absolute(λ)(λ＝1,2,3……,m)线性插值到下看LADCP的深度单元轴上，然后再对两者求算术平均，获得平均绝对垂向速度m表示上看LADCP的深度单元个数，n表示下看LADCP的深度单元个数；016)绝对回波强度计算，获得绝对回波强度，首先根据公式计算深度单元距离R：式中，B是船载ADCP盲区，为常数；P是宣称脉冲长度，为常数；W是深度单元厚度，N是深度单元个数，θ是波束角，cave是船载ADCP某一深度单元的观测范围内的平均声速；根据公式计算发射能量K1：式中，K1c、a、b和c2均是船载ADCP的出厂参数，均为常数，Vs是整个投放过程中高电压值的平均值，高电压值直接从船载ADCP获得，K1为船载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频率有关的系统固定项，为常数，Tx是当前船载ADCP探头所处位置的海水温度，由步骤S1获得，Ea为步骤S1中LADCP测量获得的相对回波强度，Er为有效剖面末端的RSSI值，为常数，SV是绝对回波强度；017)生物量反演，获得反演生物量，选取相同数量的若干个生物量数据和绝对回波强度数据进行数据拟合，根据拟合公式计算出回归系数项a1和b1：式中，DW为步骤S1中获得的生物量数据，经过公式计算出a1和b1后，输入所有绝对回波强度数据，根据公式得到反演生物量：步骤S3：获得上述经过处理后的各个数据后，对经过步骤S2处理获得的数据进行数据融合，生成统一公共深度轴下的数据表格。2.根据权利要求1所述的深海全水柱多环境要素信息融合方法，其特征在于：所述步骤008)中的按照深度单元平均，其具体过程为：步骤a)，以Z'深值为间隔将从上看LADCP初始位置开始至海底进行划分，得到上看LADCP的n个深度单元，分别为上看LADCP的第一深度单元、第二深度单元……；以及以Z”深值为间隔将从下看LADCP初始位置开始直至海底进行划分，得到m个深度单元，分别为下看LADCP的第一深度单元、第二深度单元……，其中，m和n均为正整数，Z'深和Z”深分别为上看LADCP和下看LADCP采集数据前设置的深度单元厚度；步骤b),首先，获取上看LADCP的n个深度单元的每个深度单元范围内所有的Vabsolute值，进行算术平均得到V′absolute(i)(i＝1,2,3……,n)值；以及获取下看LADCP的m个深度单元的每个深度单元范围内所有的Vabsolute值，进行算术平均得到V″absolute(λ)(λ＝1,2,3……,m)值；按深度将V″absolute(λ)线性插值到对应上看LADCP的各深度单元内，得到V″absolute(i)(i＝1,2,3……,n)；然后，将同一深度单元内包含线性插值后的V″absolute(i)和上看LADCP的V′absolute(i)进行算术平均，得到该深度单元内的平均绝对流速最终得到n个深度单元的各个深度单元的平均绝对流速其中，如果m＝n,则V″absolute(λ)刚好一一线性插值到对应上看LADCP的各深度单元内；如果m≠n，则对于上看LADCP第i个深度单元，在下看LADCP深度单元中选取深度最接近前者的两个深度单元，将该深度单元内的V″absolute(λ)值线性插入至上看LADCP的第i个深度单元。3.根据权利要求1所述的深海全水柱多环境要素信息融合方法，其特征在于：所述步骤015)中的按照深度单元平均，其具体过程为：步骤a)，以Z'深值为间隔将从上看LADCP初始位置开始至海底进行划分，得到上看LADCP的n个深度单元，分别为上看LADCP的第一深度单元、第二深度单元……；以及以Z”深值为间隔将从下看LADCP初始位置开始直至海底进行划分，得到m个深度单元，分别为下看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的各深度单元内；如果m≠n，则对于上看LADCP第i个深度单元，在下看LADCP深度单元中选取深度最接近前者的深度单元，将该深度单元内的ω″absolute(λ)值线性插入至上看LADCP的第i个深度单元。4.根据权利要求1所述的深海全水柱多环境要素信息融合方法，其特征在于：所述步骤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”来表示该深度下无对应测量数据，最终，生成统一公共深度轴下的数据表格。5.一种...

【专利技术属性】
技术研发人员：郭斌斌，梁前勇，董一飞，钟超，
申请(专利权)人：广州海洋地质调查局，
类型：发明
国别省市：广东,44