CN109297530A - A kind of multi-environment element information fusion method of the full water column in deep-sea and processing terminal - Google Patents

Abstract

The present invention relates to a kind of multi-environment element information fusion method of the full water column in deep-sea and processing terminals, the described method comprises the following steps: step S1: obtain include pressure, temperature, salinity, conductivity, dissolved oxygen, PH, chlorophyll, the velocity of sound, the velocity of sound of Ship-Mounted ADCP, vertical velocity, speed shearing, relative velocity, relative return intensity and biomass data initial data；Step S2: data processing is carried out according to the initial data of acquisition；Step S3: data fusion is carried out to the data obtained by step S2 processing, generates the data form under unified common depth axis.The present invention is modified environmental element information, improves the correctness and vertical resolution of parameters；The data of distinct device acquisition are combined, processing generates new environmental element, the scientific research value for the survey data that extended；And all environmental element information are fused to uniform depth common depth axis, interactivity is strong.

Description

Deep-sea full-water-column multi-environment-factor information fusion method and processing terminal

Technical Field

The invention relates to the technical field of marine investigation, in particular to a deep-sea full-water-column multi-environment-factor information fusion method and a processing terminal.

Background

Deep sea environment baseline survey needs to acquire multiple effective and high-resolution water column environment element information at the same time, and improves the identification capability of marine environment characteristics. However, data obtained by conventional investigation methods such as thermohaline depth measurement, navigation ADCP observation, biological trawling and the like are not fused, the vertical resolution is low, partial key environmental elements are lacked, and the obtained data processing results are poor in interchangeability.

Disclosure of Invention

In view of the defects of the prior art, one of the objectives of the present invention is to provide a deep-sea full-water-column multi-environment element information fusion method, which can solve the problem of fusing multiple environment element information.

The second purpose of the present invention is to provide a processing terminal that can solve the problem of merging multiple types of environmental element information and the problem of processing and generating high-quality multi-environmental element data that meets environmental baseline scientific research requirements.

The technical scheme for realizing one purpose of the invention is as follows: a deep sea full water column multi-environment element information fusion method comprises the following steps:

step S1: obtaining original data including pressure, temperature, salinity, conductivity, dissolved oxygen, PH, chlorophyll, sound velocity of shipborne ADCP, vertical velocity, velocity shear, relative flow velocity, relative echo intensity and biological data, wherein the pressure, the temperature, the salinity, the conductivity and the sound velocity are directly obtained by CTD, the dissolved oxygen, the PH and the chlorophyll are obtained by water sample observation, the sound velocity of the shipborne ADCP is obtained by the shipborne ADCP according to a temperature probe, the vertical velocity, the velocity shear and the relative flow velocity are measured by LADCP, the relative echo intensity is obtained by the LADCP and the shipborne ADCP, and the biological data is obtained by a biological trawl;

step S2: and performing data processing according to the obtained original data, wherein the data processing comprises the following steps:

001) and pressure correction, wherein the correction formula is as follows: the actual pressure value is the original pressure value — the deck pressure, where the original pressure value is the pressure in step S1, and the deck pressure is the pressure measured at the deck after the CTD is recovered;

002) temperature correction, wherein the correction formula is as follows: offset is multiplied, b represents "the number of days from the time of correction of the current temperature correction", residual represents "the temperature difference between before and after the current trip", and n represents "the number of days between the temperature correction before and after the current trip";

003) salinity correction, firstly, the conductivity data is corrected, and the correction formula is as follows ①:

wherein n represents the number of samples of the water chemistry analysis used for the calibration, α_iIndicating the conductivity of the ith water chemistry sample as measured by CTD, β_iIndicates the conductivity of the ith water chemistry analysis sample obtained from the water chemistry analysis experiment, slope indicates the corrected conductivity,

the true conductivity value is obtained according to equation ②:

true conductivity reading slope ②

The conductivity reading is the conductivity obtained in the step S1, and after the actual conductivity value is obtained, the salinity is calculated according to the seawater state equation, and the salinity is the corrected salinity value;

004) and (3) correcting the dissolved oxygen, the PH and the chlorophyll, wherein the correction of the dissolved oxygen firstly calculates a slope correction parameter according to a formula ③:

in the formula, slope_DORepresenting a slope correction parameter, DO_{Water sample}Representing dissolved oxygen, DO, obtained by analysis of a water chemistry sample at a certain depth by means of a water chemistry experiment_CTDRepresentation and DO_{Water sample}Dissolved oxygen, DO, obtained at the same depth from CTD_{Water sample}Reading the dissolved oxygen value at the same depth directly from the dissolved oxygen raw data obtained in the step S1;

after obtaining the slope correction parameter, it is calculated according to formula ④:

true dissolved oxygen value ═ dissolved oxygen reading ═ slope_DO------④

The dissolved oxygen reading is the original data of all dissolved oxygen from the sea surface to the seabed obtained by CTD, and the true value of the dissolved oxygen is the dissolved oxygen after slope correction;

correcting the PH and the chlorophyll according to formulas ③ and ④ to obtain corrected PH and chlorophyll;

005) and correcting the sound velocity to obtain a corrected relative flow velocity value, and calculating according to a formula ⑤:

V_corrected＝V_relative*(C_ctd/C_ADCP)------⑤

in the formula, V_relativeRepresenting the uncorrected relative flow rate, obtained in step S1, C_ADCPSound velocity representing the onboard ADCP obtained in step S1, C_ctdIs the sound velocity obtained in step S1, obtained in step S1, V_correctedRepresenting the corrected relative flow velocity value due to the error caused by the sound velocity;

006) low scatter correction, obtaining corrected velocity shear values, calculated according to equation ⑥:

V_correction＝V₀+ΔV------⑥

Wherein EI represents relative echo intensity and is obtained in step S1, k is a proportionality coefficient and is a constant, Δ V represents a correction value, and V is₀Representing the speed shear obtained by LADCP, obtained by step S1, V_CorrectionRepresents the velocity shear value corrected for relative flow rate deviations due to low scattering;

007) absolute flow rate calculation, first calculate the reference flow rate by shear method, formula ⑦:

in the formula, T representsThe total time of the lowering of LADCP, t represents time, V_shipIndicating the speed, V, of the vessel_referA reference flow rate is indicated and,

obtaining V_referThen calculating according to formula ⑧ to obtain the absolute flow velocity V_absolute：

V_absolute＝V_relative+V_refer------⑧

In the formula, V_absoluteRepresents the absolute flow rate;

008) double labdc absolute flow rate average: calculating under-view LADCP according to step 007) to obtain absolute flow rate, and averaging according to depth units to obtain V'_absolute(i) (i ═ 1,2,3 … …, n); meanwhile, LADCP from Shangqi was calculated as absolute flow rate in step 007) and averaged in units of depth to obtain V ″_absolute(λ) (λ ═ 1,2,3 ….., m); will V ″)_absolute(λ) (λ 1,2,3.... times.m) is linearly interpolated on the depth unit axis of the lower view lacpb, and then the two are arithmetically averaged to obtain the average absolute flow velocityWherein m represents the number of depth units of the upper-view LADCP, and n represents the number of depth units of the lower-view LADCP;

009) moisture removal: obtaining the absolute flow velocity V obtained from the step 007) by software_absoluteSubtracting the positive pressure tide profile to obtain a result after the tide is removed;

010) calculating the Richch numbers to obtain the Richch numbers R_iCalculated according to equation ⑨:

wherein g represents the acceleration of gravity, m represents the depth of m meters, n represents the depth of n meters, ρ_mRepresenting the sea water density at a depth of m meters, p_nRepresenting sea water density at a depth of n metersDegree u_m、v_mRespectively representing an east component and a north component of the absolute horizontal flow velocity at a depth of m meters, u_n、v_nRespectively representing an east component and a north component of the absolute horizontal flow velocity at the depth of n meters;

011) and calculating the vertical diffusion coefficient to obtain the vertical diffusion coefficient, wherein the calculation formula is shown as formula ⑩:

in the formula, v₀Get 10^-2m²(ii) s, α takes 5, v_bGet 10^-4m²/s，k_bGet 10^-5m²/s，K_TIs the vertical diffusivity;

012) low-pass filtering the CTD depth data: carrying out low-pass filtering on the CTD depth data by adopting a low-pass filter with the cutoff frequency of 6Hz to obtain the CTD depth data of 6 Hz;

013) calculating the vertical motion speed of the CTD to obtain the vertical motion speed omega of the CTD_CTDAccording to the formulaAnd (4) calculating:

where z (t) represents CTD depth data of 6Hz obtained in 012), and t represents time;

014) calculating absolute vertical velocity to obtain absolute vertical velocity according to formulaAnd (3) calculating:

in the formula, ω_ADCPRepresents the vertical velocity, ω represents the absolute vertical velocity, obtained by the laccp in step S1;

015) obtaining average absolute vertical velocity by averaging absolute vertical velocities of the two LADCPs, namely, obtaining the lower LADCP in two devices including the upper LADCP and the lower LADCP according to a formulaAbsolute vertical velocities were obtained and averaged in units of depth to obtain ω'_absolute(i) (i 1,2,3.. n.); meanwhile, the upper view LADCP is expressed according to the formulaThe absolute vertical velocity is obtained and averaged in units of depth to obtain ω ″_absolute(λ) (λ 1,2,3.. said., m); will omega ″)_absolute(λ) (λ 1,2,3.. said., m) is linearly interpolated on the depth unit axis of the lower view lacpb, and then the two are arithmetically averaged to obtain the average absolute vertical velocitym represents the number of depth units of the upper-view LADCP, and n represents the number of depth units of the lower-view LADCP;

016) calculating absolute echo intensity to obtain absolute echo intensity according to formulaCalculating the depth unit distance R:

wherein B is a shipborne ADCP blind area which is a constant; p is the asserted pulse length, which is a constant; w is depth cell thickness, N is number of depth cellsTheta is the beam angle, c_aveIs the average sound velocity within the observation range of a certain depth unit of the shipborne ADCP;

according to the formulaCalculating the emission energy K₁：

In the formula, K_1cA, b and c₂Are all the delivery parameters of the ship-borne ADCP, are all constants, V_sIs the average value of the high voltage values in the whole putting process, the high voltage values are directly obtained from the ship-borne ADCP, K₁Transmitting energy for the onboard ADCP;

according to the formulaCalculating the sound absorption coefficient e:

wherein,

P₂＝1-1.37×10^-4z+6.2×10^-9z²，P₃＝1-1.383×10^-5z+4.9×10^-10z²，

wherein PH, S, T and z are corresponding seawater PH, salinity, temperature and depth values of the ith depth unit of the ship-borne ADCP, and are obtained in step S1, c_iRepresents the sound velocity corresponding to the ith depth cell, freq represents the sound frequency of the ship-borne ADCP, and is constant,

then, according to the formulaCalculating absolute echo intensity:

in the formula, K₂Is the system noise factor, constant, K_sIs a system fixed term, constant, T, related to the frequency of the onboard ADCP_xIs the seawater temperature at the current position of the onboard ADCP probe, obtained in step S1, E_aRelative echo intensities obtained for the LADCP measurement in step S1, E_rIs the RSSI value at the end of the effective profile, which is constant, S_VIs the absolute echo intensity;

017) biomass inversion is carried out to obtain inversion biomass, a plurality of biomass data and absolute echo intensity data with the same quantity are selected for data fitting, and the biomass inversion is carried out according to a fitting formulaCalculating a regression coefficient term a₁And b₁：

Wherein DW is the biomass data obtained in step S1,

by the formulaCalculate a₁And b₁Then, all absolute echo intensity data are input according to the formulaObtaining an inversion biomass:

step S3: after the processed data are obtained, the data obtained through the processing in step S2 are subjected to data fusion, and a data table under the unified common depth axis is generated.

Further, the averaging according to depth unit in the step 008) specifically includes:

step a), with Z'_{Deep to}Dividing the initial position of the up-looking LADCP to the seabed at intervals to obtain n depth units of the up-looking LADCP, wherein the n depth units are … … of a first depth unit and a second depth unit of the up-looking LADCP respectively; and with Z "_{Deep to}The values are such that the interval will be divided starting from the initial position of the look-down LADCP until the seafloor to obtain m depth units, respectively the first depth unit, the second depth unit … … of the look-down LADCP,

wherein m and n are both positive integers, Z'_{Deep to}And Z "_{Deep to}The depth unit thicknesses are respectively set before data acquisition of the upper-view LADCP and the lower-view LADCP;

step b), first, all V within each depth unit range of n depth units of the up-looking LADCP are acquired_absoluteValue, obtaining V 'by arithmetic average'_absolute(i) (i 1,2,3.. n.) values; and obtaining all V within each depth unit range of m depth units of the lower-view LADCP_absoluteValue, arithmetic mean to obtain V ″)_absolute(λ) (λ ═ 1,2,3.. said, m) values; according to depth will V_absolute(lambda) linearly interpolating to each depth unit corresponding to the upper display LADCP to obtain V ″)_absolute(i)(i＝1,2,3...…N); then, the same depth unit contains V' after linear interpolation_absolute(i) And see V 'of LADCP'_absolute(i) Performing arithmetic mean to obtain average absolute flow velocity in the depth unitFinally obtaining the average absolute flow velocity of each depth unit of the n depth units

Wherein if m ═ n, then V ″)_absolute(λ) interpolating exactly one-to-one linearly into each depth unit of the corresponding upper view LADCP; if m is not equal to n, for the ith depth unit of the upper-view LADCP, selecting two depth units with the depth closest to the former from the depth units of the lower-view LADCP, and combining V ″' in the two depth units_absoluteThe (λ) value is linearly interpolated to see the ith depth element of the LADCP.

Further, the averaging according to the depth unit in the step 015) includes:

step a), with Z'_{Deep to}Dividing the initial position of the up-looking LADCP to the seabed at intervals to obtain n depth units of the up-looking LADCP, wherein the n depth units are … … of a first depth unit and a second depth unit of the up-looking LADCP respectively; and with Z "_{Deep to}The values are such that the interval will be divided starting from the initial position of the look-down LADCP until the seafloor to obtain m depth units, respectively the first depth unit, the second depth unit … … of the look-down LADCP,

wherein m and n are both positive integers, Z'_{Deep to}And Z "_{Deep to}The depth cell thicknesses set before data acquisition for the upper view LADCP and the lower view LADCP, respectively.

Step b), first, all ω values within each depth unit range of the n depth units of the top view LADCP are acquired, and arithmetic mean is performed to obtain ω'_absolute(i) (i ═ 1,2,3 … …, n) values; and obtaining m depth sheets for reading LADCPAll omega values in each depth unit range of the element are arithmetically averaged to obtain omega ″_absolute(λ) (λ ═ 1,2,3 … …, m) values; according to depth will be omega_a″_bsolute(lambda) is linearly interpolated into each depth cell corresponding to the upper viewing LADCP to obtain omega'_a'_bsolute(i) (i ═ 1,2,3 … …, n); then, the same depth unit contains omega' after linear interpolation_absolute(i) And see omega of LADCP'_absolute(i) Performing arithmetic mean to obtain average absolute flow velocity in the depth unitI.e. to obtain the average absolute flow rate of each of the n depth cells

Wherein if m ═ n, then ω ″ ", and_absolute(λ) interpolating exactly one-to-one linearly into each depth unit of the corresponding upper view LADCP; if m is not equal to n, for the ith depth unit of the upper-view LADCP, selecting the depth unit with the depth closest to the former from the depth units of the lower-view LADCP, and dividing omega' in the depth unit_absoluteThe (λ) value is linearly inserted up to the ith depth element of the LADCP.

Further, the step S3 generates a data table under the unified common depth axis, which includes the following specific processes:

step 1), firstly, starting from the sea level to the seabed depth by 1m, and sequentially increasing the distance by 1m to generate a common depth axis as a first column which is expressed by the depth;

step 2), defining a first group of data comprising the temperature obtained in the step 002), the salinity obtained in the step 003), the dissolved oxygen, the PH and the chlorophyll obtained in the step 004), and the absolute echo intensity obtained in the step 016), and a second group of data comprising the absolute echo intensity obtained in the step 016), the eastern component and the northern component of the absolute horizontal flow velocity obtained in the step 007), the Richardson number obtained in the step 010), the vertical diffusion coefficient obtained in the step 011), the absolute vertical flow velocity obtained in the step 014) and the inversion biomass obtained in the step 017);

step 3), obtaining the measured values of the temperature, salinity, dissolved oxygen, PH, chlorophyll and absolute echo intensity of the first group of data in each depth range every 1m, carrying out arithmetic average on the temperature values in the same depth range to obtain the temperature measuring result, and placing the temperature measuring result in the corresponding depth position of the table; carrying out arithmetic mean on salinity values in the same depth range to obtain salinity measurement results, and placing the salinity measurement results at the depth corresponding to the table; carrying out arithmetic mean on the dissolved oxygen values in the same depth range to obtain a measurement result of the dissolved oxygen and placing the measurement result into a corresponding depth position of a table; carrying out arithmetic mean on the solution PH values within the same depth range to obtain a PH measurement result and placing the PH measurement result in a corresponding depth position of a table; carrying out arithmetic mean on chlorophyll values in the same depth range to obtain a chlorophyll measurement result and putting the chlorophyll measurement result into a table at a corresponding depth; carrying out arithmetic mean on the absolute echo intensity values in the same depth range to obtain the measurement result of the absolute echo intensity and placing the measurement result into the depth position corresponding to the table;

and 4) acquiring a second group of data which are measured by taking the depth unit as 4m, filling each measured value of the corresponding second group of data in each depth which is an integer multiple of 4, using NaN to represent that no corresponding measured data exists in the depth in the rest depths, and finally generating a data table under the unified common depth axis.

The second technical scheme for realizing the aim of the invention is as follows: a processing terminal, comprising,

a memory for storing program instructions;

a processor for executing the program instructions to perform the steps of:

step S1: obtaining original data including pressure, temperature, salinity, conductivity, dissolved oxygen, PH, chlorophyll, sound velocity of shipborne ADCP, vertical velocity, velocity shear, relative flow velocity, relative echo intensity and biological data, wherein the pressure, the temperature, the salinity, the conductivity and the sound velocity are directly obtained by CTD, the dissolved oxygen, the PH and the chlorophyll are obtained by water sample observation, the sound velocity of the shipborne ADCP is obtained by the shipborne ADCP according to a temperature probe, the vertical velocity, the velocity shear and the relative flow velocity are measured by LADCP, the relative echo intensity is obtained by the LADCP and the shipborne ADCP, and the biological data is obtained by a biological trawl;

step S2: and performing data processing according to the obtained original data, wherein the data processing comprises the following steps:

001) and pressure correction, wherein the correction formula is as follows: the actual pressure value is the original pressure value — the deck pressure, where the original pressure value is the pressure in step S1, and the deck pressure is the pressure measured at the deck after the CTD is recovered;

002) temperature correction, wherein the correction formula is as follows: offset is multiplied, b represents "the number of days from the time of correction of the current temperature correction", residual represents "the temperature difference between before and after the current trip", and n represents "the number of days between the temperature correction before and after the current trip";

003) salinity correction, firstly, the conductivity data is corrected, and the correction formula is as follows ①:

wherein n represents the number of samples of the water chemistry analysis used for the calibration, α_iIndicating the conductivity of the ith water chemistry sample as measured by CTD, β_iIndicates the conductivity of the ith water chemistry analysis sample obtained from the water chemistry analysis experiment, slope indicates the corrected conductivity,

the true conductivity value is obtained according to equation ②:

true conductivity reading slope ②

The conductivity reading is the conductivity obtained in the step S1, and after the actual conductivity value is obtained, the salinity is calculated according to the seawater state equation, and the salinity is the corrected salinity value;

004) and (3) correcting the dissolved oxygen, the PH and the chlorophyll, wherein the correction of the dissolved oxygen firstly calculates a slope correction parameter according to a formula ③:

in the formula, slope_DORepresenting a slope correction parameter, DO_{Water sample}Representing dissolved oxygen, DO, obtained by analysis of a water chemistry sample at a certain depth by means of a water chemistry experiment_CTDRepresentation and DO_{Water sample}Dissolved oxygen, DO, obtained at the same depth from CTD_{Water sample}Reading the dissolved oxygen value at the same depth directly from the dissolved oxygen raw data obtained in the step S1;

after obtaining the slope correction parameter, it is calculated according to formula ④:

true dissolved oxygen value ═ dissolved oxygen reading ═ slope_DO------④

The dissolved oxygen reading is the original data of all dissolved oxygen from the sea surface to the seabed obtained by CTD, and the true value of the dissolved oxygen is the dissolved oxygen after slope correction;

correcting the PH and the chlorophyll according to formulas ③ and ④ to obtain corrected PH and chlorophyll;

005) and correcting the sound velocity to obtain a corrected relative flow velocity value, and calculating according to a formula ⑤:

V_corrected＝V_relative*(C_ctd/C_ADCP)------⑤

in the formula, V_relativeRepresenting the uncorrected relative flow rate, obtained in step S1, C_ADCPSound velocity representing the onboard ADCP obtained in step S1, C_ctdIs the sound velocity obtained in step S1, obtained in step S1, V_correctedIndicating errors due to speed of soundThe difference and the corrected relative flow rate value;

006) low scatter correction, obtaining corrected velocity shear values, calculated according to equation ⑥:

V_correction＝V₀+ΔV------⑥

Wherein EI represents relative echo intensity and is obtained in step S1, k is a proportionality coefficient and is a constant, Δ V represents a correction value, and V is₀Representing the speed shear obtained by LADCP, obtained by step S1, V_CorrectionRepresents the velocity shear value corrected for relative flow rate deviations due to low scattering;

007) absolute flow rate calculation, first calculate the reference flow rate by shear method, formula ⑦:

wherein T represents the total time length of the LADCP downward placing, T represents the time, and V_shipIndicating the speed, V, of the vessel_referA reference flow rate is indicated and,

obtaining V_referThen calculating according to formula ⑧ to obtain the absolute flow velocity V_absolute：

V_absolute＝V_relative+V_refer------⑧

In the formula, V_absoluteRepresents the absolute flow rate;

008) double labdc absolute flow rate average: calculating under-view LADCP according to step 007) to obtain absolute flow rate, and averaging according to depth units to obtain V'_absolute(i) (i 1,2,3.. n.); meanwhile, the LADCP from the previous view was calculated in step 007) as an absolute flow rate and averaged in units of depth,obtaining V_absolute(λ) (λ 1,2,3.. said., m); will V ″)_absolute(λ) (λ 1,2,3 ….., m) is linearly interpolated on the depth unit axis of the downward-looking LADCP, and then the two are arithmetically averaged to obtain the average absolute flow velocityWherein m represents the number of depth units of the upper-view LADCP, and n represents the number of depth units of the lower-view LADCP;

009) moisture removal: obtaining the absolute flow velocity V obtained from the step 007) by software_absoluteSubtracting the positive pressure tide profile to obtain a result after the tide is removed;

010) calculating the Richch numbers to obtain the Richch numbers R_iCalculated according to equation ⑨:

wherein g represents the acceleration of gravity, m represents the depth of m meters, n represents the depth of n meters, ρ_mRepresenting the sea water density at a depth of m meters, p_nDenotes the sea water density at a depth of n meters, u_m、v_mRespectively representing an east component and a north component of the absolute horizontal flow velocity at a depth of m meters, u_n、v_nRespectively representing an east component and a north component of the absolute horizontal flow velocity at the depth of n meters;

011) and calculating the vertical diffusion coefficient to obtain the vertical diffusion coefficient, wherein the calculation formula is shown as formula ⑩:

in the formula, v₀Get 10^-2m²(ii) s, α takes 5, v_bGet 10^-4m²/s，k_bGet 10^-5m²/s，K_TIs the vertical diffusivity;

012) low-pass filtering the CTD depth data: carrying out low-pass filtering on the CTD depth data by adopting a low-pass filter with the cutoff frequency of 6Hz to obtain the CTD depth data of 6 Hz;

013) calculating the vertical motion speed of the CTD to obtain the vertical motion speed omega of the CTD_CTDAccording to the formulaAnd (4) calculating:

where z (t) represents CTD depth data of 6Hz obtained in 012), and t represents time;

014) calculating absolute vertical velocity to obtain absolute vertical velocity according to formulaAnd (3) calculating:

in the formula, ω_ADCPRepresents the vertical velocity, ω represents the absolute vertical velocity, obtained by the laccp in step S1;

015) obtaining average absolute vertical velocity by averaging absolute vertical velocities of the two LADCPs, namely, obtaining the lower LADCP in two devices including the upper LADCP and the lower LADCP according to a formulaAbsolute vertical velocities were obtained and averaged in units of depth to obtain ω'_absolute(i) (i 1,2,3.. n.); meanwhile, the upper view LADCP is expressed according to the formulaThe absolute vertical velocity is obtained and averaged in units of depth to obtain ω ″_absolute(λ) (λ 1,2,3.. said., m); will omega ″)_absolute(λ) (λ 1,2,3.. said., m) is linearly interpolated on the depth unit axis of the lower view lacpb, and then the two are arithmetically averaged to obtain the average absolute vertical velocitym represents the number of depth units of the upper-view LADCP, and n represents the number of depth units of the lower-view LADCP;

016) calculating absolute echo intensity to obtain absolute echo intensity according to formulaCalculating the depth unit distance R:

wherein B is a shipborne ADCP blind area which is a constant; p is the asserted pulse length, which is a constant; w is depth cell thickness, N is number of depth cells, θ is beam angle, c_aveIs the average sound velocity within the observation range of a certain depth unit of the shipborne ADCP;

according to the formulaCalculating the emission energy K₁：

In the formula, K_1cA, b and c₂Are all the delivery parameters of the ship-borne ADCP, are all constants, V_sIs the average value of the high voltage values in the whole putting process, the high voltage values are directly obtained from the ship-borne ADCP, K₁Transmitting energy for the onboard ADCP;

according to the formulaCalculating the sound absorption coefficient e:

wherein,

P₂＝1-1.37×10^-4z+6.2×10^-9z²，P₃＝1-1.383×10^-5z+4.9×10^-10z²，

wherein PH, S, T and z are corresponding seawater PH, salinity, temperature and depth values of the ith depth unit of the ship-borne ADCP, and are obtained in step S1, c_iRepresents the sound velocity corresponding to the ith depth cell, freq represents the sound frequency of the ship-borne ADCP, and is constant,

then, according to the formulaCalculating absolute echo intensity:

in the formula, K₂Is the system noise factor, constant, K_sIs andthe frequency-dependent systematic fixed term, T, of the onboard ADCP is constant_xIs the seawater temperature at the current position of the onboard ADCP probe, obtained in step S1, E_aRelative echo intensities obtained for the LADCP measurement in step S1, E_rIs the RSSI value at the end of the effective profile, which is constant, S_VIs the absolute echo intensity;

017) biomass inversion is carried out to obtain inversion biomass, a plurality of biomass data and absolute echo intensity data with the same quantity are selected for data fitting, and the biomass inversion is carried out according to a fitting formulaCalculating a regression coefficient term a₁And b₁：

Wherein DW is the biomass data obtained in step S1,

by the formulaCalculate a₁And b₁Then, all absolute echo intensity data are input according to the formulaObtaining an inversion biomass:

step S3: after the processed data are obtained, the data obtained through the processing in step S2 are subjected to data fusion, and a data table under the unified common depth axis is generated.

6. The processing terminal of claim 5, wherein: the averaging according to depth units in the step 008) specifically comprises the following steps:

step a), with Z'_{Deep to}Dividing the initial position of the up-looking LADCP to the seabed at intervals to obtain n depth units of the up-looking LADCP, wherein the n depth units are … … of a first depth unit and a second depth unit of the up-looking LADCP respectively; and with Z "_{Deep to}The values are such that the interval will be divided starting from the initial position of the look-down LADCP until the seafloor to obtain m depth units, respectively the first depth unit, the second depth unit … … of the look-down LADCP,

wherein m and n are both positive integers, Z'_{Deep to}And Z "_{Deep to}The depth unit thicknesses are respectively set before data acquisition of the upper-view LADCP and the lower-view LADCP;

step b), first, all V within each depth unit range of n depth units of the up-looking LADCP are acquired_absoluteValue, obtaining V 'by arithmetic average'_absolute(i) (i ═ 1,2,3 ….., n) values; and obtaining all V within each depth unit range of m depth units of the lower-view LADCP_absoluteValue, arithmetic mean to obtain V ″)_absolute(λ) (λ ═ 1,2,3.. said, m) values; according to depth will V_absolute(lambda) linearly interpolating to each depth unit corresponding to the upper display LADCP to obtain V ″)_absolute(i) (i 1,2,3.. n.); then, the same depth unit contains V' after linear interpolation_absolute(i) And see V of LADCP_a′_bsolute(i) Performing arithmetic mean to obtain average absolute flow velocity in the depth unitFinally obtaining the average absolute flow velocity of each depth unit of the n depth units

Wherein if m ═ n, then V ″)_absolute(λ) interpolating exactly one-to-one linearly into each depth unit of the corresponding upper view LADCP; if m ≠ n, then for the ith depth unit of the look-up LADCP, atTwo depth units with the depth closest to the LADCP depth units are selected from the lower view LADCP depth units, and V' in the two depth units_absoluteThe (λ) value is linearly interpolated to see the ith depth element of the LADCP.

7. The processing terminal of claim 5, wherein: the specific process of the unit averaging according to depth in the step 015) is as follows:

step a), with Z'_{Deep to}Dividing the initial position of the up-looking LADCP to the seabed at intervals to obtain n depth units of the up-looking LADCP, wherein the n depth units are … … of a first depth unit and a second depth unit of the up-looking LADCP respectively; and with Z "_{Deep to}The values are such that the interval will be divided starting from the initial position of the look-down LADCP until the seafloor to obtain m depth units, respectively the first depth unit, the second depth unit … … of the look-down LADCP,

wherein m and n are both positive integers, Z'_{Deep to}And Z "_{Deep to}The depth cell thicknesses set before data acquisition for the upper view LADCP and the lower view LADCP, respectively.

Step b), first, all ω values within each depth unit range of the n depth units of the top view LADCP are acquired, and arithmetic mean is performed to obtain ω'_absolute(i) (i 1,2,3.. n.) values; and obtaining all omega values in each depth unit range of m depth units of the LADCP, and carrying out arithmetic mean to obtain omega ″)_absolute(λ) (λ ═ 1,2,3.. said, m) values; according to depth will omega ″)_absolute(lambda) linearly interpolating to each depth unit corresponding to the upper view LADCP to obtain omega ″)_absolute(i) (i 1,2,3.. n.); then, the same depth unit contains omega' after linear interpolation_absolute(i) And see omega of LADCP'_absolute(i) Performing arithmetic mean to obtain average absolute flow velocity in the depth unitI.e. to obtain the average absolute flow rate of each of the n depth cells

Wherein if m ═ n, then ω ″ ", and_absolute(λ) interpolating exactly one-to-one linearly into each depth unit of the corresponding upper view LADCP; if m is not equal to n, for the ith depth unit of the upper-view LADCP, selecting two depth units with the depth closest to the former from the depth units of the lower-view LADCP, and dividing omega' in the depth unit_absoluteThe (λ) value is linearly inserted up to the ith depth cell of the LADCP.

8. The processing terminal of claim 5, wherein: the specific process of generating the data table under the unified common depth axis in step S3 is as follows:

step 1), firstly, starting from the sea level to the seabed depth by 1m, and sequentially increasing the distance by 1m to generate a common depth axis as a first column which is expressed by the depth;

step 2), defining a first group of data comprising the temperature obtained in the step 002), the salinity obtained in the step 003), the dissolved oxygen, the PH and the chlorophyll obtained in the step 004), and the absolute echo intensity obtained in the step 016), and a second group of data comprising the absolute echo intensity obtained in the step 016), the eastern component and the northern component of the absolute horizontal flow velocity obtained in the step 007), the Richardson number obtained in the step 010), the vertical diffusion coefficient obtained in the step 011), the absolute vertical flow velocity obtained in the step 014) and the inversion biomass obtained in the step 017);

step 3), obtaining the measured values of the temperature, salinity, dissolved oxygen, PH, chlorophyll and absolute echo intensity of the first group of data in each depth range every 1m, carrying out arithmetic average on the temperature values in the same depth range to obtain the temperature measuring result, and placing the temperature measuring result in the corresponding depth position of the table; carrying out arithmetic mean on salinity values in the same depth range to obtain salinity measurement results, and placing the salinity measurement results at the depth corresponding to the table; carrying out arithmetic mean on the dissolved oxygen values in the same depth range to obtain a measurement result of the dissolved oxygen and placing the measurement result into a corresponding depth position of a table; carrying out arithmetic mean on the solution PH values within the same depth range to obtain a PH measurement result and placing the PH measurement result in a corresponding depth position of a table; carrying out arithmetic mean on chlorophyll values in the same depth range to obtain a chlorophyll measurement result and putting the chlorophyll measurement result into a table at a corresponding depth; carrying out arithmetic mean on the absolute echo intensity values in the same depth range to obtain the measurement result of the absolute echo intensity and placing the measurement result into the depth position corresponding to the table;

and 4) acquiring a second group of data which are measured by taking the depth unit as 4m, filling each measured value of the corresponding second group of data in each depth which is an integer multiple of 4, using NaN to represent that no corresponding measured data exists in the depth in the rest depths, and finally generating a data table under the unified common depth axis.

The invention has the beneficial effects that: the method corrects the environmental element information parameters, improves the correctness of each parameter, further improves the correctness of the baseline survey and evaluation of the hydrate environment, fuses all the environmental element information to the uniform depth public depth axis, and has strong interactivity.

Drawings

FIG. 1 is a flow chart of a preferred embodiment of the present invention;

fig. 2 is a schematic structural diagram of a processing terminal according to the present invention.

Detailed Description

The invention will be further described with reference to the accompanying drawings and the detailed description below:

as shown in fig. 1 and fig. 2, a deep sea full water column multi-environment element information fusion method comprises the following steps:

step S1: obtaining original data including pressure, temperature, salinity, conductivity, dissolved oxygen, PH, chlorophyll, sound velocity of shipborne ADCP, vertical velocity, velocity shear, relative flow velocity, relative echo intensity and biological quantity data, wherein the pressure, the temperature, the salinity, the conductivity and the sound velocity are directly obtained by CTD (Chinese name is thermohaloscope) equipment, the dissolved oxygen, the PH and the chlorophyll are obtained by observing a water sample, the sound velocity of the shipborne ADCP is obtained by the shipborne ADCP according to the sound velocity value obtained by a temperature probe, the vertical velocity, the velocity shear and the relative flow velocity are measured by LADCP, the relative echo intensity is obtained by the LADCP and the shipborne ADCP, the biological quantity data is obtained by a biological trawl, and the original data and the obtained way are shown in the following table:

in this embodiment, the LADCP is WHS ● ADCP manufactured by RDI (Teledyne RD Instruments, abbreviated as RDI), and in this embodiment, two LADCPs are used in some steps, which are defined as a lower view LADCP and an upper view LADCP, and the shipborne ADCP is OS ● ADCP manufactured by RDI;

all the original data are data measured from the sea surface to the seabed and are data obtained by measuring a certain ocean depth, and the data comprise depth information, such as temperature values, which represent the temperature corresponding to the certain ocean depth;

step S2: and performing data processing according to the obtained original data, wherein the data processing comprises the following steps:

001) and pressure correction, wherein the correction formula is as follows: the actual pressure value is the original pressure value — the deck pressure, where the original pressure value is the pressure in step S1, and the deck pressure is the pressure measured at the deck after the CTD is recovered;

002) temperature correction, wherein the correction formula is as follows: b, residual/n, where b is a multiplication, b denotes "the number of days from the time of calibration of the current trip to the temperature correction, residual denotes" the temperature difference between before and after the trip ", n denotes" the number of days between the temperature calibration before and after the trip ", and offset denotes a corrected temperature value, for example, when the temperature calibration is performed once on 1/4 th day before the trip and the temperature calibration is performed on 1/7 th day after the trip, and the temperature data of 1/5 th day during the trip is calibrated, n is 90, b is 30, and residual is the temperature after the calibration on 1/7 th day;

003) correcting salinity, firstly calibrating the conductivity data, wherein the calibration formula is as follows:

wherein n represents the number of samples of water chemistry analysis used for calibration, α_iIndicating the conductivity of the ith water chemistry sample as measured by CTD, β_iIndicates the conductivity of the ith water chemistry analysis sample obtained from the water chemistry analysis experiment, slope indicates the corrected conductivity,

obtaining a conductivity true value according to the corrected conductivity, wherein a calculation formula is as follows: the actual conductivity value is the conductivity reading slope, the conductivity reading indicates all the conductivity data measured by the CTD in the whole measuring process from the sea surface to the seabed, that is, the conductivity obtained in step S1, after the actual conductivity value is obtained, the salinity is calculated according to the seawater state equation, and the salinity is the corrected salinity value;

004) dissolved oxygen, PH and chlorophyll are corrected, taking the dissolved oxygen as an example, firstly, a slope correction parameter is calculated, and the calculation formula is as follows:

in the formula, slope_DORepresenting a slope correction parameter, DO_{Water sample}Representing dissolved oxygen, DO, obtained by analysis of a water chemistry sample at a certain depth by means of a water chemistry experiment_CTDRepresentation and DO_{Water sample}Dissolved oxygen, DO, obtained at the same depth from CTD_{Water sample}Can be directly obtained from step S1The dissolved oxygen value under the same depth can be read from the original dissolved oxygen data, the certain depth can be generally in the deep sea below 1000 meters, and only DO (dissolved oxygen) is ensured_{Water sample}And DO_CTDIs dissolved oxygen obtained at the same depth;

after obtaining the slope correction parameter, performing slope correction on the dissolved oxygen profile measured by the CTD equipment, and outputting a new dissolved oxygen profile, wherein the calculation formula is as follows: true dissolved oxygen value ═ dissolved oxygen reading ═ slope_DOThe dissolved oxygen reading is the original data of all dissolved oxygen from the sea surface to the seabed obtained by CTD, and the true value of the dissolved oxygen is the dissolved oxygen after slope correction; the correction process of PH and chlorophyll is similar, and the specific process is not described in detail;

005) and (3) sound velocity correction, wherein the sound velocity is calculated by using a fixed salinity value and self-low-accuracy temperature data by the ship-borne ADCP, and further a relative flow velocity calculation error is caused, so that the relative flow velocity caused by the sound velocity needs to be corrected by generating an error, and the correction formula is as follows:

V_corrected＝V_uncorrected*(C_ctd/C_ADCP)

in the formula, V_uncorrectedA value representing the uncorrected relative flow rate, obtained in step S1, C_ADCPRepresenting the value of the sound velocity, C, obtained by the onboard ADCP from the temperature probe_ctdIs the sound velocity, V, obtained in step S1_correctedRepresenting the calibrated relative flow velocity value caused by error caused by sound velocity;

006) in the low scattering correction, under the deep sea low scattering environment, the echo of the LADCP is significantly reduced, which causes deviation of the measured speed shear, so that the correction of the speed shear error caused by low scattering needs to be performed, and the correction formula is as follows:

V_correction＝V₀+ΔV

Where EI represents the relative echo intensity, obtained in step S1, k is a proportionality coefficient, a constant is artificially selected, and in this embodiment, the value is-0.12, Δ V represents a correction value, and V represents a constant₀Representing the velocity shear obtained by the LADCP, i.e. representing the current gradient, V, between adjacent depth horizons_CorrectionRepresents the velocity shear value corrected for relative flow rate deviations due to low scatter;

007) calculating absolute flow velocity, wherein the relative velocity measured by the conventional LADCP is the relative velocity between the LADCP and a water body, and is not the absolute velocity, so the absolute flow velocity needs to be calculated, and can be obtained by a shearing method:

the shearing method comprises the steps of calculating a ship drift distance, determining a reference flow velocity, averaging the flow velocity vertical shearing of units with the same depth to obtain the flow velocity vertical shearing of a full water column, performing vertical integration from the sea bottom to the sea surface to obtain a relative flow velocity profile, calculating the ship drift distance, and determining the reference flow velocity, wherein the calculation formula is as follows:

wherein T represents the total time length of the LADCP downward placing, T represents the time, and V_shipIndicating the speed, V, of the vessel_relativeThe relative flow rate, V, obtained in step S1_referA reference flow rate is indicated and,

obtaining V_referThen calculating according to the formula (1) to obtain the absolute flow velocity V_absolute：

V_absolute＝V_relative+V_refer(1)

Absolute flow velocity V_absoluteIs a vector, passing through pair V_absoluteOrthogonal decomposition is carried out to decompose an east component u and a north component v of the absolute horizontal flow velocity,

008) double LADCP insulatorAnd (4) averaging the convection velocity: obtaining an average absolute flow velocityIn the actual survey of the marine environmental element information, two devices including the up-looking LADCP and the down-looking LADCP are generally adopted, and the down-looking LADCP is calculated according to the step 007) to obtain V'_absolute(i) (i ═ 1,2,3.... times.n) averaged by depth unit, and the resulting V ″ "was calculated from the previous view of the lacdp in step 007)_absolute(λ) (λ 1,2,3.. said., m), and converting V ″ ", to a corresponding one_absolute(λ) (λ ═ 1,2,3 … …, m) is linearly interpolated on the depth cell axis of the view-down LADCP, and then the two are averaged as follows:

step a), obtaining the same depth cell Z for the up-looking and down-looking LADCPs_{Deep to}A value of and with Z_{Deep to}Dividing the LADCP from the initial position of the upper view to the seabed by the value of interval to obtain n depth units of the upper view LADCP, namely a first depth unit and a second depth unit … … of the upper view LADCP, and Z depth units_{Deep to}The values are such that the interval will be divided starting from the initial position of the look-down ladp until the seafloor to obtain m depth units, respectively the first depth unit, the second depth unit … … of the look-down ladp,

the m and the n are positive integers, and the initial positions of the upper view LADCP and the lower view LADCP are not necessarily the same, so that the value of m is not necessarily equal to the value of n;

step b), first, all V 'within each depth unit range of the n depth units of the top view LADCP are acquired'_absolute(i) (i ═ 1,2,3 … …, n) values, and all V ″'s in each depth cell range obtained for m depth cells of the view LADCP_absolute(λ) (λ 1,2,3.. said., m) value, V ″ ", is calculated using a method of calculating a value of (λ ═ 1,2,3.. said., m)_absolute(lambda) linearly interpolating each depth unit corresponding to the upper display level LADCP, and including the linearly interpolated V ″' in the same depth unit_absolute(λ) and see V 'of LADCP above'_absolute(i) Performing arithmetic mean to obtain average absolute flow velocity in the depth unitFurther obtaining the average absolute flow velocity of each depth unit of the n depth units

Wherein if m ═ n, then V ″)_absolute(λ) is exactly one-to-one linearly interpolated into each depth cell corresponding to the upper view LADCP, if m>n, all V' from the nth depth unit to the mth depth unit of the LADCP are seen_absoluteThe (λ) values are all linearly inserted up to the nth depth cell of the LADCP, if m<n, no V ″' of the LADCP is seen from the m +1 th depth unit to the nth depth unit of the LADCP_absolute(lambda) inserting the (alpha) into the (alpha) substrate,

the value of the depth unit is determined by the LADCP, the upper-view LADCP has the same depth unit as the lower-view LADCP, each LADCP has a fixed depth unit, the depth units determined by different LADCPs may be different, the depth unit here refers to the largest range that the LADCP can measure from the current position downwards, for example, the depth unit is 4m, and represents data that the LADCP can measure from the current position downwards within 4 m; as described above, as the labdc moves from the sea surface to the sea bottom, the labdc obtains a plurality of depth units including a first depth unit, a second depth unit, and so on, and the value of each depth unit is the same, so that the labdc can obtain a plurality of different values in each depth unit, and arithmetically average all the different values in the depth unit, and the value obtained after arithmetically averaging represents the average absolute flow velocity in the depth unit;

009) dehumidification, due to the absolute flow velocity V obtained by LADCP_absoluteContains moisture, and needs to be deducted, namely, removed moisture, and the process is as follows: acquiring longitude, latitude and time data of LADCP data by Earth&Tide Model Driver (TMD) software, v2.05 version, produced by Space Research corporation, calls the regional Tide mode (Egbert and Erofeeva 2002) to obtain a positive pressure Tide profile, and obtains the absolute flow rate obtained in step 007)V_absoluteDeducting the positive pressure tide profile to obtain a result after moisture removal;

010) calculating the Richardson number according to a formula (2) to obtain the Richardson number R_i：

Wherein g represents gravity acceleration, the value is 9.8, rho is sea water density, z represents sea depth, u and V represent V respectively_absoluteAn east component of absolute horizontal flow velocity and a north component of absolute horizontal flow velocity.

The formula (2) is generally understood as a formula that calculates a definite function expression (explicit, implicit, inverse, parametric, etc.) by applying a derivation rule, but in practical problems, only some discrete data of independent variables and dependent variables are obtained; for the sake of convenience of calculation, in the present embodiment, equation (2) is converted into the following equation for calculation:

wherein m represents a depth of m meters, n represents a depth of n meters, ρ_mRepresenting the sea water density at a depth of m meters, p_nDenotes the sea water density at a depth of n meters, u_m、v_mRespectively representing the east component and the north component of the horizontal flow velocity at a depth of m meters, u_n、v_nRespectively representing an east component and a north component of the absolute horizontal flow velocity at the depth of n meters;

011) calculating the vertical diffusion coefficient, wherein the calculation formula is shown as an equation (3):

wherein v is₀Get 10^-2m²(ii)/s, α is 5, R_i010) obtained Richcson number, v_bGet 10^-4m²/s， k_bGet 10^-5m²/s，K_TIs the vertical diffusivity;

012) CTD depth data low-pass filtering, and removing high-frequency data through the low-pass filter, wherein the specific process comprises the following steps: carrying out low-pass filtering on the CTD depth data by adopting a low-pass filter with the cutoff frequency of 6Hz to obtain the CTD depth data of 6 Hz;

013) calculating the vertical movement speed of the CTD according to a formula (4) to obtain the vertical movement speed omega of the CTD_CTD：

Where z (t) represents CTD depth data of 6Hz obtained in 012), and t represents time;

014) calculating the absolute vertical speed according to the following calculation formula:

ω＝ω_ADCP-ω_CTD

in the formula, ω_ADCPRepresents the vertical velocity, ω represents the absolute vertical velocity, obtained by the laccp in step S1;

015) obtaining average absolute vertical velocity by averaging absolute vertical velocity omega data of the double-LADCP, calculating the lower-view LADCP in two devices including the upper-view LADCP and the lower-view LADCP according to step 014) to obtain the absolute vertical velocity, and averaging according to depth units to obtain omega'_absolute(i) (i ═ 1,2,3.... times.n), while the upper view laccp was calculated in step 014) to obtain the absolute vertical velocity, and averaged in units of depth to obtain ω ″_absolute(λ) (λ ═ 1,2,3.. said., m), and converting ω ″ ", to ω ″_absolute(λ) (λ 1,2,3.... times.m) is linearly interpolated onto the depth unit axis of the look-down LADCP, and then two pairs are aligned againThe arithmetic mean is obtained to obtain the mean absolute vertical velocityThe specific calculation process is similar to that of the step 008), and is not described herein;

016) absolute echo intensity calculation, namely converting relative echo intensity data obtained by the LADCP and the ship-borne ADCP into absolute echo intensity, and firstly calculating a sound absorption coefficient e, a depth unit distance R and emission energy K₁Then, the absolute echo intensity is obtained through a sound wave calibration equation, and the detailed calculation process is as follows:

the formula for calculating the depth cell distance R is as follows:

where B is the shipborne ADCP dead zone, an artificially selected constant, P is the declared pulse length, a constant, W is the depth cell thickness, N is the number of depth cells, θ is the beam angle, C is the number of depth cells_ADCPIs the default sound velocity for the onboard ADCP, obtained directly from the onboard ADCP, c_aveIs the average sound velocity over the observation range from the shipborne ADCP probe to a certain depth unit;

calculating the emission energy K₁The calculation formula is as follows:

K₁＝0.34K_1C(V_s·a-b)/c₂

in the formula, K_1cA, b and c₂All are shipping parameters of the shipborne ADCP, and in the embodiment, the values are respectively 3.9, 0.17, 5.767 and 34.247, V_sIs the average value of the high voltage values in the whole putting process, the high voltage values are directly obtained from the ship-borne ADCP, K₁Transmitting energy for the onboard ADCP;

calculating the sound absorption coefficient e, wherein the calculation formula is as follows:

wherein,

P₂＝1-1.37×10^-4z+6.2×10^-9z²，P₃＝1-1.383×10^-5z+4.9×10^-10z²，

wherein PH, S, T and z are corresponding seawater PH, salinity, temperature and depth values of the ith depth unit of the ship-borne ADCP, and are obtained from the step S1, c_iRepresenting the sound velocity corresponding to the ith depth unit; freq represents the sound frequency of the ADCP on board, and is constant,

then, the absolute echo intensity is calculated according to the following formula:

in the formula, K₂Is the system noise factor, constant, K_sIs a system fixed term related to the frequency of the ADCP on board the ship, and is a constant, in the embodiment, K₂And K_sAre taken to be 3.6 and 4.17 x 105, respectively, T_xIs the seawater temperature at the current position of the shipborne ADCP probe, which can be obtained by CTD, E_aRelative echo intensities obtained for the LADCP measurement in step S1, E_rThe RSSI value at the end of the effective profile is a constant, which is 40S in this embodiment_VIs an absolute echoStrength;

017) biomass inversion, which comprises the following steps: selecting a plurality of biomass data and absolute echo intensities with the same quantity, in the embodiment, selecting 50 biomass data and absolute echo intensity data to perform data fitting, and calculating a regression coefficient term a according to a fitting formula₁And b₁The fitting formula is as follows:

in the formula, DW is the biomass data obtained in step S1, and DW and S are also included_VBoth time and depth are related, so when calculating the above equation, DW and S_VA is calculated for data with the same depth at the same time through the formula₁And b₁Then, a calculation formula (6) of the biomass data DW is obtained according to formula (5):

through the calculation, a plurality of data can be selected, a formula (6) is obtained after fitting, and then the absolute echo intensity S is obtained_VThe biomass data of the whole arbitrary depth from the sea surface to the seabed can be calculated by substituting the formula (6), namely the biomass data calculated by the formula (6) is inversion biomass, so that the workload of the whole actual measurement from the sea surface to the seabed through the biological trawl is reduced;

step S3: after the processed data are obtained, data fusion is performed on the data obtained through the processing in step S2, and a data table under the unified common depth axis is generated, which includes the following specific steps:

the first column is depth, and the common depth axis is generated by sequentially increasing from 1m to 1m until the sea floor, the depth of the first column represents the sea depth from the sea level, for example, 2m represents the sea depth from the sea level by 2m, and the last value of the first column represents the depth from the sea floor to the sea level;

defining a first set of data comprising the temperature obtained in step 002), the salinity obtained in step 003), the dissolved oxygen, the PH, the chlorophyll obtained in step 004), and the absolute echo intensity obtained in step 016), and a second set of data comprising the absolute echo intensity obtained in step 016), the eastern component and the northern component of the absolute horizontal flow velocity obtained in step 007), the Richardson number obtained in step 010), the vertical diffusion coefficient obtained in step 011), the absolute vertical flow velocity obtained in step 014), and the inverted biomass obtained in step 017);

it should be noted that, in the actual survey and measurement process, the CTD, the LADCP and the shipborne ADCP are measured from the sea surface to the seabed, and all the obtained raw data include depth and time information corresponding to the distance from the sea surface to the seabed, that is, any one of the data includes depth and time information, such as a temperature value, which represents a temperature corresponding to a certain depth at a certain time; in this embodiment, the depth unit measured by the CTD is 1m, that is, the CTD can measure data within 1m of the current position of the CTD, so that in one depth unit, the number of any parameter in the first set of data measured by the CTD is related to the measurement frequency of the CTD and the descending speed thereof, if the measurement frequency of the CTD is 24 times per second and the descending speed is 1m/s, each parameter of the first set of data can measure 24 data within a depth interval of 1m, and then the 24 data near the single depth unit are respectively subjected to arithmetic averaging, so that only one data is present at each depth of 1m, 2m, 3m … …, and the like; in this embodiment, the depth unit of the shipborne ADCP and LADCP measurement is 4m, and if the measurement frequency of the LADCP measurement is once per second and the descent speed is 1m/s, the second group of data respectively measures one data at the depth of integral multiple of 4m, that is, only one data is respectively measured at each depth of 4m, 8m, 12m … …, and the like;

the depth axis of the first set of data and the depth axis of the second set of data are therefore not identical, and for a uniform comparison, the first set of data and the second set of data are fused uniformly onto a uniform common depth axis, and furthermore, the parameters can be obtained directly from the corresponding devices, whether the depth unit measured by CTD or the depth units measured by shipborne ADCP and LADCP is determined by their respective devices themselves;

since the second set of data is data measured in depth units of 4m, there is only one data at each depth at intervals of 4m, i.e. there is only one data at each depth of 4m, 8m, 12m … …, etc., and the remaining depth is represented by "NaN" to indicate that there is no corresponding measured data at that depth, therefore, the resolution of the second set of data is 4 m;

the first set of data, which may have a plurality of values measured over a depth range of 0.5-1.5m, is the arithmetic mean of all the measured values, the arithmetic mean representing the measurement at a depth of 1m, for example, n temperature values x measured over a depth range of 0.5-1.5m_i(i 1,2,3 … …, n), at a depth of 1m Similarly, the arithmetic mean of all values measured in the depth range of 1.5-2.5m at the depth of 2m represents the result at the depth of 2m, and the calculation is carried out progressively from the first depth unit backwards, and finally a first set of data is obtained, wherein the data are measured at the depth interval of 1m, namely, only one data is measured at each depth of 1m, 2m, 3m … … and the like, so that the resolution of the first set of data is 1 m;

the table formed by the first set of data and the second set of data obtained finally is shown in table one:

watch 1

A table with quality control marks can be generated corresponding to the first table, data quality is represented by '0', data quality is represented by '1', data quality is unknown by '2', data quality is poor, initial default values are all 1, a worker can use the table to mark environment element parameters of various depths according to experience and actual requirements, when the data quality of a certain parameter corresponding to a certain depth is judged to be good, the position is marked as '0', and when the data quality of a certain parameter corresponding to a certain depth is judged to be poor, the position is marked as '2', so that visual observation is provided, and the second table is as follows:

watch two

In addition, the present invention also relates to a processing terminal 100 of a physical device implementing the above method, which comprises,

a memory 101 for storing program instructions;

a processor 102 for executing the program instructions to perform the steps of:

step S1: obtaining original data including pressure, temperature, salinity, conductivity, dissolved oxygen, PH, chlorophyll, sound velocity of shipborne ADCP, vertical velocity, velocity shear, relative flow velocity, relative echo intensity and biological data, wherein the pressure, the temperature, the salinity, the conductivity and the sound velocity are directly obtained by CTD, the dissolved oxygen, the PH and the chlorophyll are obtained by water sample observation, the sound velocity of the shipborne ADCP is obtained by the shipborne ADCP according to a temperature probe, the vertical velocity, the velocity shear and the relative flow velocity are measured by LADCP, the relative echo intensity is obtained by the LADCP and the shipborne ADCP, and the biological data is obtained by a biological trawl;

step S2: and performing data processing according to the obtained original data, wherein the data processing comprises the following steps:

001) and pressure correction, wherein the correction formula is as follows: the actual pressure value is the original pressure value — the deck pressure, where the original pressure value is the pressure in step S1, and the deck pressure is the pressure measured at the deck after the CTD is recovered;

002) temperature correction, wherein the correction formula is as follows: offset is multiplied, b represents "the number of days from the time of correction of the current temperature correction", residual represents "the temperature difference between before and after the current trip", and n represents "the number of days between the temperature correction before and after the current trip";

003) salinity correction, firstly, the conductivity data is corrected, and the correction formula is as follows ①:

wherein n represents the number of samples of the water chemistry analysis used for the calibration, α_iIndicating the conductivity of the ith water chemistry sample as measured by CTD, β_iIndicates the conductivity of the ith water chemistry analysis sample obtained from the water chemistry analysis experiment, slope indicates the corrected conductivity,

the true conductivity value is obtained according to equation ②:

true conductivity reading slope ②

The conductivity reading is the conductivity obtained in the step S1, and after the actual conductivity value is obtained, the salinity is calculated according to the seawater state equation, and the salinity is the corrected salinity value;

004) and (3) correcting the dissolved oxygen, the PH and the chlorophyll, wherein the correction of the dissolved oxygen firstly calculates a slope correction parameter according to a formula ③:

in the formula, slope_DORepresenting a slope correction parameter, DO_{Water sample}Representing dissolved oxygen, DO, obtained by analysis of a water chemistry sample at a certain depth by means of a water chemistry experiment_CTDRepresentation and DO_{Water sample}Dissolved oxygen, DO, obtained at the same depth from CTD_{Water sample}In order to directly read the dissolved oxygen value at the same depth from the original dissolved oxygen data obtained in step S1, after obtaining the slope correction parameter, calculating according to the formula ④:

true dissolved oxygen value ═ dissolved oxygen reading ═ slope_DO------④

The dissolved oxygen reading is the original data of all dissolved oxygen from the sea surface to the seabed obtained by CTD, and the true value of the dissolved oxygen is the dissolved oxygen after slope correction;

correcting the PH and the chlorophyll according to formulas ③ and ④ to obtain corrected PH and chlorophyll;

005) and correcting the sound velocity to obtain a corrected relative flow velocity value, and calculating according to a formula ⑤:

V_corrected＝V_relative*(C_ctd/C_ADCP)------⑤

in the formula, V_relativeRepresenting the uncorrected relative flow rate, obtained in step S1, C_ADCPSound velocity representing the onboard ADCP obtained in step S1, C_ctdIs the sound velocity obtained in step S1, obtained in step S1, V_correctedRepresenting the corrected relative flow velocity value due to the error caused by the sound velocity;

006) low scatter correction, obtaining corrected velocity shear values, calculated according to equation ⑥:

V_correction＝V₀+ΔV------⑥

Wherein EI represents relative echo intensity and is obtained in step S1, k is a proportionality coefficient and is a constant, Δ V represents a correction value, and V is₀Representing the speed shear obtained by LADCP, obtained by step S1, V_CorrectionRepresents the velocity shear value corrected for relative flow rate deviations due to low scattering;

007) absolute flow rate calculation, first calculate the reference flow rate by shear method, formula ⑦:

wherein T represents the total time length of the LADCP downward placing, T represents the time, and V_shipIndicating the speed, V, of the vessel_referA reference flow rate is indicated and,

obtaining V_referThen calculating according to formula ⑧ to obtain the absolute flow velocity V_absolute：

V_absolute＝V_relative+V_refer------⑧

In the formula, V_absoluteRepresents the absolute flow rate;

008) double labdc absolute flow rate average: obtaining an average absolute flow velocityCalculating V 'obtained by using under-view LADCP in two devices including under-view LADCP and up-view LADCP according to step 007)'_absoluteAnd is prepared from V'_absoluteAveraging by depth unit, and calculating the calculated V ″' in step 007) from LADCP_absoluteAnd will V ″)_absoluteLinearly interpolating to the depth unit axis of the lower viewing LADCP, and then arithmetically averaging the two to obtain the average absolute flow velocity

009) Removing tide, calculating the positive pressure tide profile through software, obtaining the absolute flow velocity V from the step 007)_absoluteSubtracting the positive pressure tide profile to obtain a result after the tide is removed;

010) calculating the Richch numbers to obtain the Richch numbers R_iCalculated according to equation ⑨:

wherein g represents the acceleration of gravity, m represents the depth of m meters, n represents the depth of n meters, ρ_mRepresenting the sea water density at a depth of m meters, p_nDenotes the sea water density at a depth of n meters, u_mRepresenting the east component, u, of the absolute horizontal flow velocity at a depth of m meters_nRepresenting the north component of absolute horizontal flow velocity at a depth of n meters;

011) calculating vertical diffusion coefficient to obtain vertical diffusion coefficientShown in the figure:

in the formula, v₀Get 10^-2m²(ii) s, α takes 5, v_bGet 10^-4m²/s，k_bGet 10^-5m²/s，K_TIs the vertical diffusivity;

012) performing low-pass filtering on the CTD depth data by adopting a low-pass filter with the cutoff frequency of 6Hz to obtain the CTD depth data of 6 Hz;

013) calculating the vertical motion speed of the CTD to obtain the vertical motion speed omega of the CTD_CTDAccording to the formulaAnd (4) calculating:

where z (t) represents CTD depth data of 6Hz obtained in 012), and t represents time;

014) calculating absolute vertical velocity to obtain absolute vertical velocity according to formulaAnd (3) calculating:

in the formula, ω_ADCPRepresents the vertical velocity, ω represents the absolute vertical velocity, obtained by the laccp in step S1;

015) averaging absolute vertical velocities of the two LADCPs to obtain an average absolute vertical velocity, calculating omega 'obtained by the calculation of the lower LADCP in the two devices including the upper LADCP and the lower LADCP according to step 014), averaging the omega' according to depth units, calculating omega 'obtained by the calculation of the upper LADCP according to step 014), linearly interpolating omega' on a depth unit axis of the lower LADCP, and arithmetically averaging the two to obtain the average absolute vertical velocity

016) Calculating absolute echo intensity to obtain absolute echo intensity according to formulaCalculating the depth unit distance R:

where B is the shipborne ADCP dead zone and is a constant, P is the asserted pulse length and is a constant, W is the depth cell thickness, N is the number of depth cells, θ is the beam angle, c_aveIs the average sound velocity within the observation range of a certain depth unit of the shipborne ADCP;

according to the formulaCalculating the emission energy K₁：

In the formula, K_1cA, b and c₂Are all the delivery parameters of the ship-borne ADCP, are all constants, V_sIs the average value of the high voltage values in the whole putting process, the high voltage values are directly obtained from the ship-borne ADCP, K₁Transmitting energy for the onboard ADCP;

according to the formulaCalculating the sound absorption coefficient e:

wherein,

P₂＝1-1.37×10^-4z+6.2×10^-9z²，P₃＝1-1.383×10^-5z+4.9×10^-10z²，

wherein PH, S, T and z are corresponding seawater PH, salinity, temperature and depth values of the ith depth unit of the ship-borne ADCP, and are obtained in step S1, c_iRepresents the sound velocity corresponding to the ith depth cell, freq represents the sound frequency of the ship-borne ADCP, and is constant,

then, according to the formulaCalculating absolute echo intensity:

in the formula, K₂Is the system noise factor, constant, K_sIs a system fixed term, constant, T, related to the frequency of the onboard ADCP_xIs the seawater temperature at the current position of the onboard ADCP probe, obtained in step S1, E_aRelative echo intensities obtained for the LADCP measurement in step S1, E_rIs the RSSI value at the end of the effective profile, which is constant, S_VIs the absolute echo intensity;

017) biomass inversion is carried out to obtain inversion biomass, a plurality of biomass data and absolute echo intensity with the same quantity are selected for data fitting, and the biomass inversion is carried out according to a fitting formulaCalculating a regression coefficient term a₁And b₁：

Wherein DW is the biomass data obtained in step S1,

by the formulaCalculate a₁And b₁Then according to the formulaObtaining inversion biomass:

step S3: after the processed data are obtained, the data obtained through the processing in step S2 are subjected to data fusion, and a data table under the unified common depth axis is generated.

Further, the specific process is as follows:

step a) is prepared from Z'_{Deep to}Dividing the LADCP from the initial position of the upper view to the seabed by the interval to obtain n depth units of the upper view LADCP, namely a first depth unit … …, a second depth unit … … and Z'_{Deep to}The values are such that the interval will be divided from the initial position of the look-down LADCP to the seafloor to obtain m depth units, respectively the first depth unit and the second depth unit … … of the look-down LADCP,

wherein m and n are both positive integers, Z'_{Deep to}And Z "_{Deep to}The depth unit thicknesses are respectively set before data acquisition of the upper-view LADCP and the lower-view LADCP;

step b), first, all V within each depth unit range of n depth units of the up-looking LADCP are acquired_absoluteValue, obtaining V 'by arithmetic average'_absolute(i) (i 1,2,3.... n) values, and all V's within each of the m depth units of the look-down LADCP are obtained_absoluteValue of, goArithmetic mean gives V ″)_absolute(λ) (λ 1,2,3.. said., m) values, V ″, in depth unit order_absolute(lambda) linearly interpolating each depth unit corresponding to the upper display level LADCP, and including the linearly interpolated V ″' in the same depth unit_absolute(λ) and see V 'of LADCP above'_absolute(i) Performing arithmetic mean to obtain average absolute flow velocity in the depth unitFinally obtaining the average absolute flow velocity of each depth unit of the n depth units

Wherein if m ═ n, then V ″)_absolute(λ) interpolating exactly one-to-one linearly into each depth unit of the corresponding upper view LADCP; if m is not equal to n, for the ith depth unit of the up-looking LADCP, selecting two down-looking LADCP depth units with the depths closest to the ith depth unit of the up-looking LADCP from the up-looking LADCP depth units, and selecting V' of the two down-looking LADCP depth units closest to the ith depth unit of the up-looking LADCP_absoluteThe (λ) value is linearly interpolated to the ith depth element of the up-looking LADCP, where "closest" means: sequentially carrying out subtraction on the ith depth unit of the upper view LADCP and all the depth units of the lower view LADCP, and taking the absolute value of the subtraction result to obtain two minimum absolute values so as to obtain V' of the two depth units of the lower view LADCP corresponding to the two minimum absolute values_absolute(λ) value and the V ″' of the two depth cells_absoluteThe (λ) value will be linearly inserted into the ith depth cell of the look-up LADCP.

Further, the averaging according to the depth unit in the step 015) includes:

step a) is prepared from Z'_{Deep to}Dividing the LADCP from the initial position of the upper view LADCP to the seabed at intervals to obtain n depth units of the upper view LADCP, namely a first depth unit and a second depth unit … … of the upper view LADCP respectively, so as to obtain a first depth unit and a second depth unit of the upper view LADCPAnd with Z "_{Deep to}The values are such that the interval will be divided from the initial position of the look-down LADCP to the seafloor to obtain m depth units, respectively the first depth unit and the second depth unit … … of the look-down LADCP,

wherein m and n are both positive integers, Z'_{Deep to}And Z "_{Deep to}The depth unit thicknesses are respectively set before data acquisition of the upper-view LADCP and the lower-view LADCP;

step b), first, all ω values within each depth unit range of the n depth units of the top view LADCP are acquired, and arithmetic mean is performed to obtain ω'_absolute(i) (i ═ 1,2,3 … …, n) values, and all ω values within each depth unit range of the m depth units of the view LADCP were obtained, and arithmetic mean was performed to obtain ω ″_absolute(λ) (λ ═ 1,2,3 … …, m) values, dividing ω ″, by depth_absolute(lambda) linearly interpolating to each depth unit corresponding to the upper view LADCP, and then including omega' after linear interpolation in the same depth unit_absolute(λ) and ω 'of top view LADCP'_absolute(i) Performing arithmetic mean to obtain average absolute flow velocity in the depth unitFinally obtaining the average absolute flow velocity of each depth unit of the n depth units

Wherein if m ═ n, then ω ″ ", and_absolute(lambda) just linearly interpolating to each depth unit corresponding to the upper view LADCP one by one, if m is not equal to n, selecting two lower view LADCP depth units with the depths closest to the i-th depth unit of the upper view LADCP from the lower view LADCP depth units, and enabling the two lower view LADCP depth units closest to the i-th depth unit of the upper view LADCP to have omega ″' of the two lower view LADCP depth units_absoluteThe (λ) value is linearly interpolated to see the ith depth element of the LADCP.

Further, the step S3 generates a data table under the unified common depth axis, which includes the following specific processes:

step 1), firstly, starting from the sea level to the seabed depth by 1m, and sequentially increasing the distance by 1m to generate a common depth axis as a first column which is expressed by the depth;

step 2), defining a first group of data comprising the temperature obtained in step 002), the salinity obtained in step 003), the dissolved oxygen, the PH, the chlorophyll obtained in step 004), and the absolute echo intensity obtained in step 016); the second group of data comprises the absolute echo intensity obtained in the step 016), the east component and the north component of the absolute horizontal flow velocity obtained in the step 007), the Richardson number obtained in the step 010), the vertical diffusion coefficient obtained in the step 011), the absolute vertical flow velocity obtained in the step 014) and the inversion biomass obtained in the step 017);

step 3), obtaining the measured values of temperature, salinity, dissolved oxygen, PH, chlorophyll and absolute echo intensity of the first group of data in each depth range every 1m, carrying out arithmetic mean on the temperature values in the same depth range to obtain the measured results of temperature and putting the measured results into the corresponding depth positions of the table, carrying out arithmetic mean on the salinity values in the same depth range to obtain the measured results of salinity and putting the measured results into the corresponding depth positions of the table, carrying out arithmetic mean on the dissolved oxygen values in the same depth range to obtain the measured results of dissolved oxygen and putting the measured results into the corresponding depth positions of the table, carrying out arithmetic mean on the dissolved PH values in the same depth range to obtain the measured results of PH and putting the measured results into the corresponding depth positions of the table, carrying out arithmetic mean on the chlorophyll values in the same depth range to obtain the measured results of chlorophyll and putting the measured results of chlorophyll into the corresponding depth positions of the table, carrying out arithmetic mean on the absolute echo intensity values in the same depth range to obtain the measurement result of the absolute echo intensity and placing the measurement result into the depth position corresponding to the table;

and 4) acquiring a second group of data which are measured by taking the depth unit as 4m, filling each measured value of the corresponding second group of data in each depth which is an integer multiple of 4, using NaN to represent that no corresponding measured data exists in the depth in the rest depths, and finally generating a data table under the unified common depth axis.

Various other changes and modifications to the above-described embodiments and concepts may occur to those skilled in the art, and all such changes and modifications are intended to be included within the scope of the present invention as defined in the appended claims.

Claims (8)

1. A deep sea full water column multi-environment element information fusion method is characterized by comprising the following steps: which comprises the following steps:

step S1: acquiring original data including pressure, temperature, salinity, conductivity, dissolved oxygen, PH, chlorophyll, sound velocity of shipborne ADCP, vertical velocity, velocity shear, relative flow velocity, relative echo intensity and biomass data, wherein the pressure, the temperature, the salinity, the conductivity and the sound velocity are directly obtained by CTD, the dissolved oxygen, the PH and the chlorophyll are obtained by water sample observation, the sound velocity of the shipborne ADCP is obtained by the shipborne ADCP according to a temperature probe, the vertical velocity, the velocity shear and the relative flow velocity are measured by LADCP, the relative echo intensity is obtained by the LADCP and the shipborne ADCP, and the biomass data is obtained by a biological trawl;

step S2: and performing data processing according to the obtained original data, wherein the data processing comprises the following steps:

001) and pressure correction, wherein the correction formula is as follows: the actual pressure value is the original pressure value — the deck pressure, where the original pressure value is the pressure in step S1, and the deck pressure is the pressure measured at the deck after the CTD is recovered;

002) temperature correction, wherein the correction formula is as follows: offset is multiplied, b represents "the number of days from the time of correction of the current temperature correction", residual represents "the temperature difference between before and after the current time", and n represents "the number of days between the temperature correction before and after the current time";

003) salinity correction, firstly, the conductivity data is corrected, and the correction formula is as follows ①:

wherein n represents the number of samples of the water chemistry analysis used for the calibration, α_iIndicating the conductivity of the ith water chemistry analysis sample as measured by the CTD, β_iIndicates the conductivity of the ith water chemistry analysis sample obtained from the water chemistry analysis experiment, slope indicates the corrected conductivity,

the true conductivity value is obtained according to equation ②:

true conductivity reading slope ②

The conductivity reading is the conductivity obtained in the step S1, and after the actual conductivity value is obtained, the salinity is calculated according to the seawater state equation, and the salinity is the corrected salinity value;

004) and (3) correcting the dissolved oxygen, the PH and the chlorophyll, wherein the correction of the dissolved oxygen firstly calculates a slope correction parameter according to a formula ③:

in the formula, slope_DORepresenting a slope correction parameter, DO_{Water sample}Representing dissolved oxygen, DO, obtained by analysis of a water chemistry sample at a certain depth by a water chemistry experiment_CTDRepresentation and DO_{Water sample}Dissolved oxygen, DO, obtained at the same depth from CTD_{Water sample}Reading the dissolved oxygen value at the same depth directly from the dissolved oxygen raw data obtained in step S1;

after obtaining the slope correction parameter, it is calculated according to formula ④:

true dissolved oxygen value ═ dissolved oxygen reading ═ slope_DO------④

The dissolved oxygen reading is the original data of all dissolved oxygen from the sea surface to the seabed obtained by CTD, and the true value of the dissolved oxygen is the dissolved oxygen after slope correction;

correcting the pH and the chlorophyll according to formulas ③ and ④ to obtain corrected pH and chlorophyll;

005) and correcting the sound velocity to obtain a corrected relative flow velocity value, and calculating according to a formula ⑤:

V_corrected＝V_relative*(C_ctd/C_ADCP)------⑤

in the formula, V_relativeRepresenting the uncorrected relative flow rate, obtained in step S1, C_ADCPSound velocity representing the onboard ADCP obtained in step S1, C_ctdIs the sound velocity obtained in step S1, obtained in step S1, V_correctedRepresenting the corrected relative flow velocity value due to the error caused by the sound velocity;

006) low scatter correction, obtaining corrected velocity shear values, calculated according to equation ⑥:

wherein EI represents relative echo intensity and is obtained in step S1, k is a proportionality coefficient and is a constant, Δ V represents a correction value, and V is₀Representing the speed shear obtained by LADCP, obtained by step S1, V_CorrectionCorrected for low scatter induced relative flow rate deviationA speed shear value;

007) calculating absolute flow rate by calculating a reference flow rate by a shearing method according to the formula ⑦:

wherein T represents the total time length of the LADCP downward placing, T represents the time, and V_shipIndicating the speed, V, of the vessel_referA reference flow rate is indicated and,

obtaining V_referThen calculating according to formula ⑧ to obtain the absolute flow velocity V_absolute：

V_absolute＝V_relative+V_refer------⑧

In the formula, V_absoluteRepresents the absolute flow rate;

008) double labdc absolute flow rate average: calculating under-looking LADCP as step 007) to obtain absolute flow rates and averaging in units of depth to obtain V'_absolute(i) (i ═ 1,2,3 … …, n); meanwhile, LADCP from Shangqi was calculated as absolute flow rate in step 007) and averaged in units of depth to obtain V ″_absolute(λ) (λ ═ 1,2,3 … …, m); will V ″)_absolute(λ) (λ ═ 1,2,3 … …, m) is linearly interpolated on the depth cell axis looking down at the LADCP, and then the two are arithmetically averaged to obtain the mean absolute flow velocityWherein m represents the number of depth units of the upper-view LADCP, and n represents the number of depth units of the lower-view LADCP;

009) moisture removal: obtaining the absolute flow velocity V obtained from the step 007) by software_absoluteSubtracting the positive pressure tide profile to obtain a result after the tide is removed;

010) calculating the Richch numbers to obtain the Richch numbers R_iCalculated according to equation ⑨:

wherein g represents the acceleration of gravity, m represents the depth of m meters, n represents the depth of n meters, ρ_mRepresenting the sea water density at a depth of m meters, p_nDenotes the sea water density at a depth of n meters, u_m、v_mRespectively representing an east component and a north component of the absolute horizontal flow velocity at a depth of m meters, u_n、v_nRespectively representing an east component and a north component of the absolute horizontal flow velocity at the depth of n meters;

011) and calculating the vertical diffusion coefficient to obtain the vertical diffusion coefficient, wherein the calculation formula is shown as formula ⑩:

in the formula, v₀Get 10^-2m²(ii) s, α takes 5, v_bGet 10^-4m²/s，k_bGet 10^-5m²/s，K_TIs the vertical diffusivity;

012) low-pass filtering the CTD depth data: carrying out low-pass filtering on the CTD depth data by adopting a low-pass filter with the cutoff frequency of 6Hz to obtain the CTD depth data of 6 Hz;

013) calculating the vertical motion speed of the CTD to obtain the vertical motion speed omega of the CTD_CTDAccording to the formulaAnd (3) calculating:

where z (t) represents CTD depth data of 6Hz obtained in 012), and t represents time;

014) calculating absolute vertical velocity to obtain absolute vertical velocity according to formulaAnd (3) calculating:

in the formula, ω_ADCPRepresents the vertical velocity, ω represents the absolute vertical velocity, obtained by the laccp in step S1;

015) obtaining average absolute vertical velocity by averaging absolute vertical velocities of the two LADCPs, namely, obtaining the lower LADCP in two devices including the upper LADCP and the lower LADCP according to a formulaAbsolute vertical velocities were obtained and averaged in units of depth to obtain ω'_absolute(i) (i ═ 1,2,3 … …, n); meanwhile, the upper view LADCP is expressed according to the formulaThe absolute vertical velocity is obtained and averaged in units of depth to obtain ω ″_absolute(λ) (λ ═ 1,2,3 … …, m); will omega ″)_absolute(λ) (λ ═ 1,2,3 … …, m) is linearly interpolated on the depth cell axis of the bottom view laccp, and then the two are arithmetically averaged to obtain the mean absolute vertical velocitym represents the number of depth units of the upper-view LADCP, and n represents the number of depth units of the lower-view LADCP;

016) calculating absolute echo intensity to obtain absolute echo intensity according to formulaCalculating the depth cell distance R:

wherein B is a shipborne ADCP blind area which is a constant; p is the asserted pulse length, which is a constant; w is depth cell thickness, N is number of depth cells, θ is beam angle, c_aveIs a ship-borne ADCPAn average speed of sound within an observation range of a depth cell;

according to the formulaCalculating the emission energy K₁：

In the formula, K_1cA, b and c₂Are all the delivery parameters of the ship-borne ADCP, are all constants, V_sIs the average value of the high voltage values in the whole putting process, the high voltage values are directly obtained from the ship-borne ADCP, K₁Transmitting energy for the onboard ADCP;

according to the formulaCalculating the sound absorption coefficient e:

wherein, P₂＝1-1.37×10^{- 4}z+6.2×10^-9z²，P₃＝1-1.383×10^-5z+4.9×10^-10z²，

wherein PH, S, T and z are corresponding seawater PH, salinity, temperature and depth values of the ith depth unit of the ship-borne ADCP, and are obtained in step S1, c_iRepresents the sound velocity corresponding to the ith depth cell, freq represents the sound frequency of the ship-borne ADCP, and is constant,

then, according to the formulaCalculating absolute echo intensity:

in the formula, K₂Is the system noise factor, constant, K_sIs a system fixed term, constant, T, related to the frequency of the onboard ADCP_xIs the seawater temperature at the current position of the onboard ADCP probe, obtained in step S1, E_aRelative echo intensities obtained for the LADCP measurement in step S1, E_rIs the RSSI value at the end of the effective profile, which is constant, S_VIs the absolute echo intensity;

017) biomass inversion is carried out to obtain inversion biomass, a plurality of biomass data and absolute echo intensity data with the same quantity are selected for data fitting, and the biomass inversion is carried out according to a fitting formulaCalculating a regression coefficient term a₁And b₁：

Wherein DW is the biomass data obtained in step S1,

by the formulaCalculate a₁And b₁Then, all absolute echo intensity data are input according to the formulaObtaining inversion biomass:

step S3: after the processed data are obtained, the data obtained through the processing in step S2 are subjected to data fusion, and a data table under the unified common depth axis is generated.

2. The deep-sea full-water-column multi-environment-element information fusion method according to claim 1, characterized in that: the averaging according to depth units in the step 008) specifically comprises the following steps:

step a), with Z'_{Deep to}Dividing the initial position of the up-looking LADCP to the seabed at intervals to obtain n depth units of the up-looking LADCP, wherein the n depth units are a first depth unit and a second depth unit … … of the up-looking LADCP respectively; and with Z "_{Deep to}The values are such that the interval will be divided from the initial position of the look-down LADCP to the seafloor to obtain m depth units, respectively the first depth unit and the second depth unit … … of the look-down LADCP,

wherein m and n are both positive integers, Z'_{Deep to}And Z "_{Deep to}The depth unit thicknesses are respectively set before data acquisition of the upper-view LADCP and the lower-view LADCP;

step b), first, all V within each depth unit range of n depth units of the up-looking LADCP are acquired_absoluteValue, obtaining V 'by arithmetic average'_absolute(i) (i ═ 1,2,3 … …, n) values; and obtaining all V within each depth unit range of m depth units of the lower-view LADCP_absoluteValue, arithmetic mean to obtain V ″)_absolute(λ) (λ ═ 1,2,3 … …, m) values; according to depth will V_absolute(lambda) linearly interpolating to each depth unit corresponding to the upper display LADCP to obtain V ″)_absolute(i) (i ═ 1,2,3 … …, n); then, the same depth unit contains V' after linear interpolation_absolute(i) And see V 'of LADCP'_absolute(i) Performing arithmetic mean to obtain average absolute flow velocity in the depth unitTo obtain each of n depth cellsMean absolute flow velocity

Wherein if m ═ n, then V ″)_absolute(λ) interpolating exactly one-to-one linearly into each depth unit corresponding to the upper view laccp; if m ≠ n, for the ith depth unit of the upper-view LADCP, two depth units with the depth closest to the former depth unit are selected from the depth units of the lower-view LADCP, and V ″' in the depth unit is used_absoluteThe (λ) value is linearly interpolated to see the ith depth element of the LADCP.

3. The deep-sea full-water-column multi-environment-element information fusion method according to claim 1, characterized in that: the specific process of the unit averaging according to depth in the step 015) is as follows:

step a), with Z'_{Deep to}Dividing the initial position of the up-looking LADCP to the seabed at intervals to obtain n depth units of the up-looking LADCP, wherein the n depth units are a first depth unit and a second depth unit … … of the up-looking LADCP respectively; and with Z "_{Deep to}The values are such that the interval will be divided from the initial position of the look-down LADCP to the seafloor to obtain m depth units, respectively the first depth unit and the second depth unit … … of the look-down LADCP,

wherein m and n are both positive integers, Z'_{Deep to}And Z "_{Deep to}The depth unit thicknesses are respectively set before data acquisition of the upper-view LADCP and the lower-view LADCP.

Step b), first, all ω values within each depth unit range of the n depth units of the top view LADCP are acquired, and arithmetic mean is performed to obtain ω'_absolute(i) (i ═ 1,2,3 … …, n) values; and obtaining all omega values in each depth unit range of m depth units of the LADCP, and carrying out arithmetic mean to obtain omega ″)_absolute(λ) (λ ═ 1,2,3 … …, m) values; according to depth will omega ″)_absolute(lambda) linearly interpolating to each depth unit corresponding to the upper view LADCP to obtain omega ″)_absolute(i) (i ═ 1,2,3 … …, n); then, the same depth unit contains omega' after linear interpolation_absolute(i) And see omega of LADCP'_absolute(i) Performing arithmetic mean to obtain average absolute flow velocity in the depth unitI.e. to obtain the average absolute flow rate of each of the n depth cells

Wherein if m ═ n, then ω ″ ", and_absolute(λ) interpolating exactly one-to-one linearly into each depth unit corresponding to the upper view laccp; if m ≠ n, for the ith depth unit of the up-view LADCP, the depth unit with the depth closest to the former depth is selected from the depth units of the down-view LADCP, and omega ″' in the depth unit is used_absoluteThe (λ) value is linearly interpolated to see the ith depth element of the LADCP.

4. The deep-sea full-water-column multi-environment-element information fusion method according to claim 1, characterized in that: the specific process of generating the data table under the unified common depth axis in step S3 is as follows:

step 1), firstly, starting from the sea level to the seabed depth by 1m, and sequentially increasing the distance by 1m to generate a common depth axis as a first column which is expressed by the depth;

step 2), defining a first group of data comprising the temperature obtained in the step 002), the salinity obtained in the step 003), the dissolved oxygen, the PH and the chlorophyll obtained in the step 004), and the absolute echo intensity obtained in the step 016), and a second group of data comprising the absolute echo intensity obtained in the step 016), the northward component of the absolute horizontal flow velocity and the northward component of the absolute horizontal flow velocity obtained in the step 007), the Richardson number obtained in the step 010), the vertical diffusion coefficient obtained in the step 011), the absolute vertical flow velocity obtained in the step 014) and the inversion biomass obtained in the step 017);

step 3), obtaining the measured values of the temperature, salinity, dissolved oxygen, PH, chlorophyll and absolute echo intensity of the first group of data in each depth range every 1m, carrying out arithmetic average on the temperature values in the same depth range to obtain the temperature measuring result, and placing the temperature measuring result in the corresponding depth position of the table; carrying out arithmetic mean on salinity values in the same depth range to obtain salinity measurement results, and placing the salinity measurement results at the depth corresponding to the table; carrying out arithmetic mean on the dissolved oxygen values in the same depth range to obtain a measurement result of the dissolved oxygen and placing the measurement result into a corresponding depth position of a table; carrying out arithmetic mean on the solution PH values within the same depth range to obtain a PH measurement result and placing the PH measurement result in a corresponding depth position of a table; carrying out arithmetic mean on chlorophyll values in the same depth range to obtain a chlorophyll measurement result and placing the chlorophyll measurement result in a table at a corresponding depth; carrying out arithmetic mean on the absolute echo intensity values in the same depth range to obtain the measurement result of the absolute echo intensity and placing the measurement result into the depth position corresponding to the table;

and 4) acquiring a second group of data which are measured by taking the depth unit as 4m, filling each measured value of the corresponding second group of data in each depth of integral multiple of 4, using NaN to represent that no corresponding measured data exists in the depth in the rest depths, and finally generating a data table under the unified common depth axis.

5. A processing terminal, characterized by: which comprises the steps of preparing a mixture of a plurality of raw materials,

a memory for storing program instructions;

a processor for executing the program instructions to perform the steps of:

step S1: acquiring original data including pressure, temperature, salinity, conductivity, dissolved oxygen, PH, chlorophyll, sound velocity of shipborne ADCP, vertical velocity, velocity shear, relative flow velocity, relative echo intensity and biomass data, wherein the pressure, the temperature, the salinity, the conductivity and the sound velocity are directly obtained by CTD, the dissolved oxygen, the PH and the chlorophyll are obtained by water sample observation, the sound velocity of the shipborne ADCP is obtained by the shipborne ADCP according to a temperature probe, the vertical velocity, the velocity shear and the relative flow velocity are measured by LADCP, the relative echo intensity is obtained by the LADCP and the shipborne ADCP, and the biomass data is obtained by a biological trawl;

step S2: and performing data processing according to the obtained original data, wherein the data processing comprises the following steps:

001) and pressure correction, wherein the correction formula is as follows: the actual pressure value is the original pressure value — the deck pressure, where the original pressure value is the pressure in step S1, and the deck pressure is the pressure measured at the deck after the CTD is recovered;

002) temperature correction, wherein the correction formula is as follows: offset is multiplied, b represents "the number of days from the time of correction of the current temperature correction", residual represents "the temperature difference between before and after the current time", and n represents "the number of days between the temperature correction before and after the current time";

003) salinity correction, firstly, the conductivity data is corrected, and the correction formula is as follows ①:

wherein n represents the number of samples of the water chemistry analysis used for the calibration, α_iIndicating the conductivity of the ith water chemistry analysis sample as measured by the CTD, β_iIndicates the conductivity of the ith water chemistry analysis sample obtained from the water chemistry analysis experiment, slope indicates the corrected conductivity,

the true conductivity value is obtained according to equation ②:

true conductivity reading slope ②

The conductivity reading is the conductivity obtained in the step S1, and after the actual conductivity value is obtained, the salinity is calculated according to the seawater state equation, and the salinity is the corrected salinity value;

004) and (3) correcting the dissolved oxygen, the PH and the chlorophyll, wherein the correction of the dissolved oxygen firstly calculates a slope correction parameter according to a formula ③:

in the formula, slope_DORepresenting a slope correction parameter, DO_{Water sample}Representing dissolved oxygen, DO, obtained by analysis of a water chemistry sample at a certain depth by a water chemistry experiment_CTDRepresentation and DO_{Water sample}Dissolved oxygen obtained from CTD at the same depth，DO_{Water sample}Reading the dissolved oxygen value at the same depth directly from the dissolved oxygen raw data obtained in step S1;

after obtaining the slope correction parameter, it is calculated according to formula ④:

true dissolved oxygen value ═ dissolved oxygen reading ═ slope_DO------④

The dissolved oxygen reading is the original data of all dissolved oxygen from the sea surface to the seabed obtained by CTD, and the true value of the dissolved oxygen is the dissolved oxygen after slope correction;

correcting the pH and the chlorophyll according to formulas ③ and ④ to obtain corrected pH and chlorophyll;

005) and correcting the sound velocity to obtain a corrected relative flow velocity value, and calculating according to a formula ⑤:

V_corrected＝V_relative*(C_ctd/C_ADCP)------⑤

in the formula, V_relativeRepresenting the uncorrected relative flow rate, obtained in step S1, C_ADCPSound velocity representing the onboard ADCP obtained in step S1, C_ctdIs the sound velocity obtained in step S1, obtained in step S1, V_correctedRepresenting the corrected relative flow velocity value due to the error caused by the sound velocity;

006) low scatter correction, obtaining corrected velocity shear values, calculated according to equation ⑥:

wherein EI represents relative echo intensity and is obtained in step S1, k is a proportionality coefficient and is a constant, Δ V represents a correction value, and V is₀Representing the speed shear obtained by LADCP, obtained by step S1, V_CorrectionRepresenting the velocity shear value corrected for relative flow rate deviations due to low scattering;

007) absolute flow rate calculation, first calculate the reference flow rate by shear method, formula ⑦:

wherein T represents the total time length of the LADCP downward placing, T represents the time, and V_shipIndicating the speed, V, of the vessel_referA reference flow rate is indicated and,

obtaining V_referThen calculating according to formula ⑧ to obtain the absolute flow velocity V_absolute：

V_absolute＝V_relative+V_refer------⑧

In the formula, V_absoluteRepresents the absolute flow rate;

008) absolute flow rate averaging of Dual LADCP, calculating the Absolute flow rate of the next LADCP as step 007), and averaging in units of depth to obtain V'_absolute(i) (i ═ 1,2,3 … …, n); meanwhile, LADCP from Shangqi was calculated as absolute flow rate in step 007) and averaged in units of depth to obtain V ″_absolute(λ) (λ ═ 1,2,3 … …, m); will V ″)_absolute(λ) (λ ═ 1,2,3 … …, m) is linearly interpolated on the depth cell axis looking down at the LADCP, and then the two are arithmetically averaged to obtain the mean absolute flow velocityWherein m represents the number of depth units of the upper-view LADCP, and n represents the number of depth units of the lower-view LADCP;

009) moisture removal: obtaining the absolute flow velocity V obtained from the step 007) by software_absoluteSubtracting the positive pressure tide profile to obtain a result after the tide is removed;

010) calculating the Richch numbers to obtain the Richch numbers R_iCalculated according to equation ⑨:

wherein g represents the acceleration of gravity, m represents the depth of m meters, n represents the depth of n meters, ρ_mRepresenting the sea water density at a depth of m meters, p_nDenotes the sea water density at a depth of n meters, u_m、v_mRespectively representing the east component of the absolute horizontal flow velocity at a depth of m metersAnd the north component of the absolute horizontal flow velocity, u_n、v_nRespectively representing an east component and a north component of the absolute horizontal flow velocity at the depth of n meters;

011) and calculating the vertical diffusion coefficient to obtain the vertical diffusion coefficient, wherein the calculation formula is shown as formula ⑩:

in the formula, v₀Get 10^-2m²(ii) s, α takes 5, v_bGet 10^-4m²/s，k_bGet 10^-5m²/s，K_TIs the vertical diffusivity;

012) low-pass filtering the CTD depth data: carrying out low-pass filtering on the CTD depth data by adopting a low-pass filter with the cutoff frequency of 6Hz to obtain the CTD depth data of 6 Hz;

013) calculating the vertical motion speed of the CTD to obtain the vertical motion speed omega of the CTD_CTDAccording to the formulaAnd (3) calculating:

where z (t) represents CTD depth data of 6Hz obtained in 012), and t represents time;

014) calculating absolute vertical velocity to obtain absolute vertical velocity according to formulaAnd (3) calculating:

ω＝ω_ADCP-ω_CTD------

in the formula, ω_ADCPIndicates the vertical velocity, which is obtained by LADCP in step S1, and ω indicates absolute velocityFor vertical velocity;

015) obtaining average absolute vertical velocity by averaging absolute vertical velocities of the two LADCPs, namely, obtaining the lower LADCP in two devices including the upper LADCP and the lower LADCP according to a formulaAbsolute vertical velocities were obtained and averaged in units of depth to obtain ω'_absolute(i) (i ═ 1,2,3 … …, n); meanwhile, the upper view LADCP is expressed according to the formulaThe absolute vertical velocity is obtained and averaged in units of depth to obtain ω ″_absolute(λ) (λ ═ 1,2,3 … …, m); will omega ″)_absolute(λ) (λ ═ 1,2,3 … …, m) is linearly interpolated on the depth cell axis of the bottom view laccp, and then the two are arithmetically averaged to obtain the mean absolute vertical velocitym represents the number of depth units of the upper-view LADCP, and n represents the number of depth units of the lower-view LADCP;

016) calculating absolute echo intensity to obtain absolute echo intensity according to formulaCalculating the depth cell distance R:

wherein B is a shipborne ADCP blind area which is a constant; p is the asserted pulse length, which is a constant; w is depth cell thickness, N is number of depth cells, θ is beam angle, c_aveIs the average sound velocity within the observation range of a certain depth unit of the shipborne ADCP;

according to the formulaCalculating the emission energy K₁：

In the formula, K_1cA, b and c₂Are all the delivery parameters of the ship-borne ADCP, are all constants, V_sIs the average value of the high voltage values in the whole putting process, the high voltage values are directly obtained from the ship-borne ADCP, K₁Transmitting energy for the onboard ADCP;

according to the formulaCalculating the sound absorption coefficient e:

wherein, P₂＝1-1.37×10^{- 4}z+6.2×10^-9z²，P₃＝1-1.383×10^-5z+4.9×10^-10z²，

wherein PH, S, T and z are corresponding seawater PH, salinity, temperature and depth values of the ith depth unit of the ship-borne ADCP, and are obtained in step S1, c_iRepresents the sound velocity corresponding to the ith depth cell, freq represents the sound frequency of the ship-borne ADCP, and is constant,

then, according to the formulaCalculating absolute echo intensity:

in the formula, K₂Is the system noise factor, constant, K_sIs a system fixed term, constant, T, related to the frequency of the onboard ADCP_xIs the seawater temperature at the current position of the onboard ADCP probe, obtained in step S1, E_aRelative echo intensities obtained for the LADCP measurement in step S1, E_rIs the RSSI value at the end of the effective profile, which is constant, S_VIs the absolute echo intensity;

017) biomass inversion is carried out to obtain inversion biomass, a plurality of biomass data and absolute echo intensity data with the same quantity are selected for data fitting, and the biomass inversion is carried out according to a fitting formulaCalculating a regression coefficient term a₁And b₁：

Wherein DW is the biomass data obtained in step S1,

by the formulaCalculate a₁And b₁Then, all absolute echo intensity data are input according to the formulaObtaining inversion biomass:

step S3: after the processed data are obtained, the data obtained through the processing in step S2 are subjected to data fusion, and a data table under the unified common depth axis is generated.

6. The processing terminal of claim 5, wherein: the averaging according to depth units in the step 008) specifically comprises the following steps:

step a), with Z'_{Deep to}Dividing the initial position of the up-looking LADCP to the seabed at intervals to obtain n depth units of the up-looking LADCP, wherein the n depth units are a first depth unit and a second depth unit … … of the up-looking LADCP respectively; and with Z "_{Deep to}The values are such that the interval will be divided from the initial position of the look-down LADCP to the seafloor to obtain m depth units, respectively the first depth unit and the second depth unit … … of the look-down LADCP,

wherein m and n are both positive integers, Z'_{Deep to}And Z "_{Deep to}The depth unit thicknesses are respectively set before data acquisition of the upper-view LADCP and the lower-view LADCP;

step b), first, all V within each depth unit range of n depth units of the up-looking LADCP are acquired_absoluteValue, obtaining V 'by arithmetic average'_absolute(i) (i ═ 1,2,3 … …, n) values; and obtaining all V within each depth unit range of m depth units of the lower-view LADCP_absoluteValue, arithmetic mean to obtain V ″)_absolute(λ) (λ ═ 1,2,3 … …, m) values; according to depth will V_absolute(lambda) linearly interpolating to each depth unit corresponding to the upper display LADCP to obtain V ″)_absolute(i) (i ═ 1,2,3 … …, n); then, the same depth unit contains V' after linear interpolation_absolute(i) And see V 'of LADCP'_absolute(i) Performing arithmetic mean to obtain average absolute flow velocity in the depth unitFinally obtaining the average absolute flow velocity of each depth unit of the n depth units

Wherein if m ═ n, then V ″)_absolute(λ) exactly one-to-one linear interpolation to pairsLook up in each depth unit of the LADCP; if m ≠ n, for the ith depth unit of the up-view LADCP, two depth units with the depth closest to the former depth unit are selected from the depth units of the down-view LADCP, and V ″' in the two depth units_absoluteThe (λ) value is linearly interpolated to see the ith depth element of the LADCP.

7. The processing terminal of claim 5, wherein: the specific process of the unit averaging according to depth in the step 015) is as follows:

step a), with Z'_{Deep to}Dividing the initial position of the up-looking LADCP to the seabed at intervals to obtain n depth units of the up-looking LADCP, wherein the n depth units are a first depth unit and a second depth unit … … of the up-looking LADCP respectively; and with Z "_{Deep to}The values are such that the interval will be divided from the initial position of the look-down LADCP to the seafloor to obtain m depth units, respectively the first depth unit and the second depth unit … … of the look-down LADCP,

wherein m and n are both positive integers, Z'_{Deep to}And Z "_{Deep to}The depth unit thicknesses are respectively set before data acquisition of the upper-view LADCP and the lower-view LADCP.

Step b), first, all ω values within each depth unit range of the n depth units of the top view LADCP are acquired, and arithmetic mean is performed to obtain ω'_absolute(i) (i ═ 1,2,3 … …, n) values; and obtaining all omega values in each depth unit range of m depth units of the LADCP, and carrying out arithmetic mean to obtain omega ″)_absolute(λ) (λ ═ 1,2,3 … …, m) values; according to depth will omega ″)_absolute(lambda) linearly interpolating to each depth unit corresponding to the upper view LADCP to obtain omega ″)_absolute(i) (i ═ 1,2,3 … …, n); then, the same depth unit contains omega' after linear interpolation_absolute(i) And see omega of LADCP'_absolute(i) Performing arithmetic mean to obtain average absolute flow velocity in the depth unitI.e. obtaining respective ones of the n depth cellsAverage absolute flow velocity of

Wherein if m ═ n, then ω ″ ", and_absolute(λ) interpolating exactly one-to-one linearly into each depth unit corresponding to the upper view laccp; if m ≠ n, for the ith depth unit of the up-view LADCP, the depth unit with the depth closest to the two former depth units is selected from the depth units of the down-view LADCP, and omega ″' in the two depth units is used_absoluteThe (λ) value is linearly interpolated to see the ith depth element of the LADCP.

8. The processing terminal of claim 5, wherein: the specific process of generating the data table under the unified common depth axis in step S3 is as follows:

step 1), firstly, starting from the sea level to the seabed depth by 1m, and sequentially increasing the distance by 1m to generate a common depth axis as a first column which is expressed by the depth;

step 2), defining a first group of data comprising the temperature obtained in the step 002), the salinity obtained in the step 003), the dissolved oxygen, the PH and the chlorophyll obtained in the step 004), and the absolute echo intensity obtained in the step 016), and a second group of data comprising the absolute echo intensity obtained in the step 016), the northward component of the absolute horizontal flow velocity and the northward component of the absolute horizontal flow velocity obtained in the step 007), the Richardson number obtained in the step 010), the vertical diffusion coefficient obtained in the step 011), the absolute vertical flow velocity obtained in the step 014) and the inversion biomass obtained in the step 017);

step 3), obtaining the measured values of the temperature, salinity, dissolved oxygen, PH, chlorophyll and absolute echo intensity of the first group of data in each depth range every 1m, carrying out arithmetic average on the temperature values in the same depth range to obtain the temperature measuring result, and placing the temperature measuring result in the corresponding depth position of the table; carrying out arithmetic mean on salinity values in the same depth range to obtain salinity measurement results, and placing the salinity measurement results at the depth corresponding to the table; carrying out arithmetic mean on the dissolved oxygen values in the same depth range to obtain a measurement result of the dissolved oxygen and placing the measurement result into a corresponding depth position of a table; carrying out arithmetic mean on the solution PH values within the same depth range to obtain a PH measurement result and placing the PH measurement result in a corresponding depth position of a table; carrying out arithmetic mean on chlorophyll values in the same depth range to obtain a chlorophyll measurement result and placing the chlorophyll measurement result in a table at a corresponding depth; carrying out arithmetic mean on the absolute echo intensity values in the same depth range to obtain the measurement result of the absolute echo intensity and placing the measurement result into the depth position corresponding to the table;

and 4) acquiring a second group of data which are measured by taking the depth unit as 4m, filling each measured value of the corresponding second group of data in each depth of integral multiple of 4, using NaN to represent that no corresponding measured data exists in the depth in the rest depths, and finally generating a data table under the unified common depth axis.