CN109297530B

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

Info

Publication number: CN109297530B
Application number: CN201810886541.7A
Authority: CN
Inventors: 郭斌斌; 梁前勇; 董一飞; 钟超
Original assignee: Guangzhou Marine Geological Survey
Current assignee: Guangzhou Marine Geological Survey
Priority date: 2018-08-06
Filing date: 2018-08-06
Publication date: 2019-08-23
Anticipated expiration: 2038-08-06
Also published as: CN109297530A

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. 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, part of 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 purposes of the invention is to provide a deep sea full water column multi-environment element information fusion method, which can solve the problem of fusion of multiple environment element information.

The invention also provides a processing terminal, which can solve the problem of fusing multiple environmental element information and the problem of processing and generating high-quality multi-environmental element data meeting 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: 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 ⑥:

V_correction＝V₀+△V ------⑥

Wherein EI represents relative echo intensity obtained in step S1, k is a proportionality coefficient and is 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) 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 of the look-down LADCP, and then the two are arithmetically averaged to obtain the average absolute currentSpeed measuring deviceWherein 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 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 system fixed relative to the frequency of the ADCP on boardConstant term, constant, T_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.

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 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(λ) interpolating exactly one-to-one linearly into each depth unit corresponding to the upper view laccp; if m is not equal to n, for the ith depth unit of the up-looking LADCP, selecting two depth units with the depth closest to the ith depth unit from the depth units of the down-looking LADCPElement, 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:

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.

Further, the step S3 of generating the data table under the unified common depth axis specifically includes the following steps:

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.

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:

the true conductivity value is obtained according to equation ②:

true conductivity reading slope ②

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

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

V_correction＝V₀+△V ------⑥

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;

according to the formulaCalculating the emission energy K₁：

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_iRepresenting the sound velocity, fre, corresponding to the i-th depth element_qRepresenting the acoustic frequency of the onboard ADCP, is a 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;

Wherein DW is the biomass data obtained in step S1,

step b), firstly, acquiring n depths of the upper-looking LADCPAll V's in each depth cell range of a cell_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(λ) 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, 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.

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

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, two depth units with the depth closest to the former depth unit are 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.

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 unified 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 biomass data, wherein the pressure, the temperature, the salinity, the conductivity and the sound velocity are directly obtained by CTD (Chinese name 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 biomass data is obtained by a biological trawling network, and the original data and the obtained way are shown as 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;

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 temperature, residual denotes" the temperature difference between before and after the current time ", n denotes" the number of days between the temperature calibration before and after the current time ", and offset denotes a corrected temperature value, for example, when the temperature calibration is performed once on 1/4 th day before the current time and the temperature calibration is performed on 1/7 th day after the current time, and the temperature data of 1/5 th day during the current time is calibrated, n is 90, b is 30, and residual is the temperature calibrated on 1/7 th day after the calibration;

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 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,

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 a water chemistry experiment_CTDRepresentation and DO_{Water sample}Dissolved oxygen, DO, obtained at the same depth from CTD_{Water sample}The dissolved oxygen value at the same depth can be directly read from the dissolved oxygen raw data obtained in step S1, and the above-mentioned certain depth can be usually found in deep sea below 1000 m, only DO is required to be ensured_{Water sample}And DO_CTDIs dissolved oxygen obtained at the same depth;

after obtaining the slope correction parameter, carrying out 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 herein;

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 error caused by the relative flow velocity caused by the sound velocity needs to be corrected, 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 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) low scattering correction, the echo of the LADCP is significantly reduced in the deep sea low scattering environment, resulting in deviation of measured speed shear, and therefore correction of speed shear error caused by low scattering is required, and the correction formula is:

V_correction＝V₀+△V

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

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 labdc absolute flow rate average: 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 … …, n) by depth unit average, and the upper view laccp was calculated as V ″) in step 007)_absolute(λ) (λ ═ 1,2,3 … …, m), and converting V ″ "_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}The value is that the interval will be divided from the initial position of LADCP viewed from above to the sea bottomN depth units of the up-looking LADCP are obtained, which are the first depth unit and the second depth unit … … of the up-looking LADCP, respectively, and are represented by 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 positive integers, and the value of m is not necessarily equal to the value of n because the initial positions of the upper view LADCP and the lower view LADCP are not necessarily the same;

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 … …, m) value, V ″ ", and_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 the depth units corresponding to the look-up 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 maximum range of the LADCP that can be measured from the current position downwards, for example, the depth unit is 4m, and represents the data that can be measured from the current position downwards to the range of 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, the moisture is removed, 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 of Egbert and Erofeeva 2002 to obtain a positive pressure Tide profile, and obtains the absolute flow velocity V obtained in step 007)_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 often 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 and north components 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：

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

ω＝ω_ADCP-ω_CTD

015) obtaining average absolute vertical velocity by averaging absolute vertical velocity omega data of the double LADCP, calculating the lower LADCP in two devices including the upper LADCP and the lower 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 … …, n), and at the same time, the absolute vertical velocity was calculated from the previous view of laccp in step 014), and averaged in units of depth to obtain ω ″_absolute(λ) (λ ═ 1,2,3 … …, m), and converting ω ″ "_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 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 an acoustic 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 cell;

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 the absolute echo intensity;

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 for 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, according to the formula (5), obtainingCalculation formula (6) to biomass data DW:

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 the 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, 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, 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 east component and the north component of the absolute horizontal flow velocity obtained in step 007), the richardson 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, and 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 a 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;

therefore, the depth axis of the first set of data and the depth axis of the second set of data are not consistent, the first set of data and the second set of data are uniformly fused on a uniform common depth axis for uniform comparison, and in addition, whether the depth unit measured by the CTD or the depth unit measured by the ship-borne ADCP and the LADCP is determined by the respective equipment, the parameter can be directly obtained from the corresponding equipment;

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

the first set of data, which may measure many values at a depth of 0.5-1.5m, is the arithmetic mean of all measured values, the arithmetic mean representing the ratio of the measured values at a depth of 1mN temperature values x are measured, for example, 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 backwards from the first depth unit, and finally a data set is obtained, wherein the data set is measured at the depth interval of 1m, namely, only one data set is measured at each depth of 1m, 2m, 3m … … and the like, so that the resolution of the data set is 1 m;

finally, the obtained first and second sets of data form a table as 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 at 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 a 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:

the true conductivity value is obtained according to equation ②:

true conductivity reading slope ②

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}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 ------④

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

V_correction＝V₀+△V ------⑥

Wherein EI represents relative echo intensity obtained in step S1, k is a proportionality coefficient and is 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;

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;

wherein g represents the acceleration of gravity, m represents the depth of m meters, n represents the depth of n meters, ρ_mRepresenting sea water density at a depth of m metersDegree, rho_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;

012) performing low-pass filtering on the CTD depth data, namely 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;

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 the step 014), averaging the omega' according to depth units, simultaneously calculating omega 'obtained by the calculation of the upper LADCP according to the 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

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₁：

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²，

then, according to the formulaCalculating 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:

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 … … of the upper view LADCP 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,

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 in each depth cell range that take down the m depth cells of the lacdp_absoluteValue, arithmetic mean to obtain V ″)_absolute(λ) (λ ═ 1,2,3 … …, m) values, V ″, in depth units 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(λ) exactly one-to-one linear interpolation to each of the corresponding view LADCPsA depth unit; if m is not equal to n, for the ith depth unit of the up-looking LADCP, selecting two depth units of the down-looking LADCP with the depth closest to the ith depth unit of the up-looking LADCP from the depth units of the down-looking LADCP, and selecting V' of the two depth units of the down-looking LADCP with the depth closest to the ith depth unit of the up-looking LADCP from the depth units of the down-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 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) linear interpolation to corresponding look-up LADCPIn each depth unit, the same depth unit is then included with omega' after linear interpolation_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) is just linearly interpolated into each depth unit corresponding to the up-looking LADCP one by one, if m is not equal to n, for the ith depth unit of the up-looking LADCP, two down-looking LADCP depth units with the depth closest to the ith depth unit of the up-looking LADCP are selected from the down-looking LADCP depth units, and omega' of the two down-looking LADCP depth units closest to the ith depth unit of the up-looking LADCP is added_absoluteThe (λ) value is linearly interpolated to see the ith depth element of the LADCP.

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 into the corresponding depth positions of the chlorophyll table, carrying out arithmetic mean on the absolute echo intensity values in the same depth range, obtaining a measurement result of the absolute echo intensity and placing the measurement result in a corresponding depth position of a table;

Various other changes and modifications to the above-described embodiments and concepts will become apparent to those skilled in the art from the above description, 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

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:

the true conductivity value is obtained according to equation ②:

true conductivity reading slope ②

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

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

V_correction＝V₀+△V------⑥

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

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(λ), λ is 1,2,3 … …, m is linearly interpolated on the depth unit axis of the lower view lacpp, 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;

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(λ), λ is 1,2,3 … …, m is linearly interpolated on the depth unit axis of the lower view lacpp, and then the two are arithmetically averaged to obtain the average absolute vertical 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;

according to the formulaCalculating the emission energy K₁：

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²，

then, according to the formulaCalculating absolute echo intensity:

Wherein DW is the biomass data obtained in step S1,

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 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 … …, values for n; 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 … …, values of m; according to depth will V_absolute(lambda) lineThe linear interpolation is carried out in each depth unit corresponding to the 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

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 b), first, each depth unit of the n depth units of the look-up LADCP is acquiredAll values of ω 'in the range of units are arithmetically averaged to give ω'_absolute(i) I ═ 1,2,3 … …, values for n; 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 … …, values of m; 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, two depth units with the depth closest to the former depth unit are 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.

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 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 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 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 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;

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:

the true conductivity value is obtained according to equation ②:

true conductivity reading slope ②

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

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

V_correction＝V₀+△V------⑥

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(λ), λ is 1,2,3 … …, m is linearly interpolated on the depth unit axis of the lower view lacpp, 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;

according to the formulaCalculating the emission energy K₁：

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²，

then, according to the formulaCalculating absolute echo intensity:

Wherein DW is the biomass data obtained in step S1,

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

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 … …, values for n; and obtaining m depth units of the look-down LADCPAll of V within each depth cell range of_absoluteValue, arithmetic mean to obtain V ″)_absolute(λ), λ ═ 1,2,3 … …, values of m; 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

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:

wherein m and n are both positive integers, Z'_{Deep to}And Z "_{Deep to}Respectively collecting LADCP for upper and lower viewsDepth cell thickness set before data collection;

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 … …, values for n; 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 … …, values of m; 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

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: