阅读说明：本技术 一种海洋激光雷达系统响应优化处理方法 (Response optimization processing method for marine laser radar system ) 是由 刘�东 陈扬 周雨迪 刘群 于 2019-10-29 设计创作，主要内容包括：本发明公开了一种海洋激光雷达系统响应优化处理方法,首先,优化了海洋激光雷达系统响应的器件选择和计算方法：提出了一种系统响应近似模型,可快速获取激光雷达系统响应脉宽；提出了激光器、探测器和采集卡三个器件之间的匹配选型方法,可发挥不同器件的最佳性能；其次,对因海洋激光雷达系统响应导致的误差进行了抑制处理：通过选择适合不同水体的系统响应脉宽,以减小因水面反射和系统展宽导致的水面失真区；通过修正方程对系统响应导致的浮游植物层的测量误差进行修正。本发明通过优化处理海洋激光雷达系统响应,能够便利于研制海洋激光雷达时的器件选取,降低系统响应对水体光学特性反演的影响。(The invention discloses a response optimization processing method of a marine laser radar system, which comprises the following steps of firstly, optimizing a device selection and calculation method of the response of the marine laser radar system: a system response approximate model is provided, and the response pulse width of a laser radar system can be rapidly obtained; the matching model selection method among the three devices of the laser, the detector and the acquisition card is provided, and the optimal performance of different devices can be exerted; secondly, the error caused by the response of the marine laser radar system is restrained: the system response pulse width suitable for different water bodies is selected to reduce a water surface distortion area caused by water surface reflection and system widening; and correcting the measurement error of the phytoplankton layer caused by the system response through a correction equation. According to the invention, the response of the marine laser radar system is optimized, so that the device selection in the process of developing the marine laser radar can be facilitated, and the influence of the system response on the inversion of the optical characteristics of the water body is reduced.)

1. A marine laser radar system response optimization processing method is characterized in that the specific optimization processing method comprises the following steps:

(1) establishing an approximate model for rapidly acquiring the system response pulse width according to the contribution of a laser, a detector and an acquisition card to the system response, wherein the approximate model adopts a Gaussian model, and the expression is as follows:

wherein, T'_TFor approximate system response pulse width, T_LFor laser pulse width, T_DIs the detector pulse width;

(2) calculating to obtain an approximate system response pulse width according to the approximate model, and estimating a required sampling rate SR of the acquisition card according to the approximate system response pulse width;

(3) sampling by the acquisition card according to the obtained sampling rate SR, and calculating by the system to obtain the actually measured phytoplankton layer thickness Thk;

(4) and correcting the error of the actually measured phytoplankton layer thickness caused by the system response through a correction equation to finally obtain the corrected phytoplankton layer thickness.

2. The marine lidar system response optimization processing method according to claim 1, wherein in the step (2), after the approximate system response pulse width is calculated according to the approximation model, the method further comprises: judging whether the approximate system response pulse width meets the measurement condition or not based on the turbidity of the water body, and if so, continuing to estimate the sampling rate; if not, resetting the pulse width of the laser and the detector.

3. The marine lidar system response optimization processing method according to claim 2, wherein the measurement conditions are:

wherein, T'_TFor approximate system response pulse width, K_dIs the water diffusion attenuation coefficient.

4. The marine laser radar system response optimization processing method according to claim 2, wherein the system response pulse widths corresponding to different water turbidities have a certain fluctuation range, and the fluctuation range of the system response pulse widths gradually becomes smaller as the water diffusion attenuation coefficient becomes larger; the lower boundary of the fluctuation range is the condition that the pulse width of the laser is fixed to the minimum value, and the upper boundary is the condition that the pulse width of the detector is fixed to the minimum value.

5. The marine lidar system response optimization processing method of claim 1, wherein in step (2), the process of estimating the required sampling rate SR of the acquisition card is as follows:

obtaining system response pulse width T by linear fitting of data_TAnd the sampling rate SR, the expression is:

SR＝0.8234/T_T+3.3×10⁶

then the system response pulse width T_TWith approximate System response pulse Width T'_TInstead, the estimated sampling rate of the acquisition card is obtained by using the expression.

6. The marine lidar system response optimization processing method of claim 1 or claim 4, wherein in step (2), the sampling rate SR is greater than or equal to twice the system bandwidth B.

7. The marine lidar system response optimization processing method according to claim 1, wherein in the step (4), the correction equation is:

wherein Thk' is the corrected phytoplankton layer thickness; thk isThe actually measured thickness of the phytoplankton layer; l is_TIs an approximate system response pulse length, which is defined by an approximate system response pulse width T'_TObtained by light speed matching.

8. The marine lidar system response optimization processing method of claim 1, wherein the time domain waveform of the laser is gaussian and the time domain waveform of the detector is heavy-tailed.

Technical Field

The invention belongs to the technical field of laser radars, and particularly relates to a response optimization processing method of an ocean laser radar system.

Background

The oceans occupy 71% of the area of the earth's surface, while the area of oceans currently explored by humans occupies only 5% of the total area. Marine ecosystems are a complex collection of many biological communities and their environments in the ocean, where phytoplankton form the basis of most marine ecosystems, with their net photosynthetic carbon fixation approximately equal to the sum of the carbon fixation of all terrestrial plants. The most basic productivity of the food chain drives atmospheric-ocean carbon dioxide exchange and deep-sea fuel carbon fixation, and therefore, phytoplankton productivity plays an important role in the earth's ocean-air coupled system. Furthermore, global stocks of marine phytoplankton are consumed and regrown every week, and this rapid conversion supports the net of marine phytoplankton as well as fishery resources and global food supply. However, with the increasing development of science and technology, the development demand of economy is more and more intense, and the marine environment is badly affected by pollutant discharge in industrial areas, local sea eutrophication and the like; the climate problems such as carbon dioxide emission, global warming and the like can also cause the consequences of sea level rising, ocean ecosystem imbalance and the like. Therefore, the development of efficient and diversified ocean exploration methods is urgently needed to help people to better understand the ocean and develop the ocean, and to better protect the ocean environment.

The measurement of the seawater optical parameters has important significance for researching global climate change and material circulation, and the instrument which can be effectively used for remotely sensing and detecting the internal characteristics of the ocean has great advantages. The ocean laser radar serving as an active optical remote sensing device can be used for detecting the optical characteristics of a water body, can perform remote sensing measurement on an upper ocean vertical distribution structure, and has the advantages of high space-time resolution, day-night continuous observation, global scale measurement and the like. At present, the marine laser radar can be successfully applied to relevant researches such as biological detection, phytoplankton layer distribution research, shallow sea topographic mapping, marine internal wave detection and the like.

The Chinese patent document with the publication number of CN105486664A discloses a device and a method for detecting marine phytoplankton biomass and granular organic carbon, which are based on a high spectral resolution laser radar and combined with a marine phytoplankton biomass and POC inversion algorithm, and the phytoplankton biomass and POC are synchronously inverted by utilizing a particle volume scattering equation and an effective attenuation coefficient when a scattering angle is pi in the laser radar equation.

Chinese patent publication No. CN107976686A discloses a multi-field angle marine laser radar and a field angle preference method thereof, which adopts a multi-field angle marine laser radar operating in a blue-green light band, establishes a parameter sensitivity evaluation method, preferentially selects a field angle of a multi-field angle laser radar system, and can still telemeter the vertical distribution characteristic of water parameters with high precision under the multiple scattering effect of strong scatterers such as water.

Marine lidar signals can be described by lidar equations, and therefore, the inversion of useful information from the signals is also typically based on lidar equations, such as perturbation, Fernald, and slope methods. However, the actually received laser radar echo signal is composed of a laser radar equation containing target characteristics and a laser radar system time system response convolution, and a laser, a detector and an acquisition card are main devices influencing system response. The system response can cause the original laser radar equation to change, and the inversion accuracy is reduced.

Disclosure of Invention

The invention provides a response optimization processing method for a marine laser radar system, which can improve the inversion accuracy of the marine laser radar, reduce the influence of system response and accurately obtain information in a near-water surface and a phytoplankton layer.

The technical scheme of the invention is as follows:

a marine laser radar system response optimization processing method is disclosed, wherein the marine laser radar system comprises a laser transmitting system, a receiving system and a data acquisition and processing system, the data transmitting system comprises a laser and a beam expander, the receiving system comprises a telescope, a diaphragm, a narrow-band filter, a converging lens and a detector, the data acquisition and processing system comprises an acquisition card and a computer, and the specific optimization processing method comprises the following steps:

(1) establishing an approximate model for rapidly acquiring the system response pulse width according to the contribution of a laser, a detector and an acquisition card to the system response, wherein the approximate model adopts a Gaussian model, and the expression is as follows:

wherein, T'_TFor approximate system response pulse width, T_LFor laser pulse width, T_DIs the detector pulse width;

(2) calculating to obtain an approximate system response pulse width according to the approximate model, and estimating a required sampling rate SR of the acquisition card according to the approximate system response pulse width;

(3) sampling by the acquisition card according to the obtained sampling rate SR, and calculating by the system to obtain the actually measured phytoplankton layer thickness Thk;

(4) and correcting the error of the actually measured phytoplankton layer thickness caused by the system response through a correction equation to finally obtain the corrected phytoplankton layer thickness.

The invention is based on the marine laser radar working in blue-green light wave band, and can accurately obtain the information in the near-water surface and the phytoplankton layer by adopting a response optimization processing method of a marine laser radar system.

In the step (1), the invention provides an approximate model for rapidly acquiring the response pulse width of a system, and by analyzing the broadening influence of three different laser waveforms on the response of the system, a Gaussian model is adopted because the Gaussian model conforms to the conditions of most laser waveforms. Using this model, the approximate system response pulse width T 'can be quickly calculated'_TGreatly reducing the system response pulse width T_TThe calculation time of (2).

In the step (2), after the approximate system response pulse width is calculated according to the approximation model, the method further includes: judging whether the approximate system response pulse width meets the measurement condition or not based on the turbidity of the water body, and if so, continuing to estimate the sampling rate; if not, resetting the pulse width of the laser and the detector. The measurement conditions are as follows:

wherein, T'_TFor approximate system response pulse width, K_dIs the water diffusion attenuation coefficient.

The range of the relative error of the diffusion attenuation coefficient obtained by inversion within a certain distance below the water surface is defined as a Signal Distortion Zone (SDZ), and an echo signal can be distorted to a certain extent in the signal distortion zone. The signal distortion zone is related to the type of water and the lidar system parameters used. In the invention, in order to ensure that the signal distortion zone is smaller than a certain value, the response pulse width of the laser radar system corresponding to different water body types is known to have a certain range by analyzing the relationship between the system response and the water diffusion attenuation coefficient. As the diffusion attenuation coefficient of the water body is continuously increased, the response pulse width range of the system is gradually reduced. When the system response pulse width and the water diffusion attenuation coefficient meet the formula, the ratio of the signal distortion area range to the maximum detectable depth is smaller than 20%, and a more accurate laser radar signal can be obtained so as to invert more accurate water optical parameters.

In the step (2), the process of estimating the required sampling rate SR of the acquisition card is as follows:

obtaining system response pulse width T by linear fitting of data_TAnd the sampling rate SR, the expression is:

SR＝0.8234/T_T+3.3×10⁶

then the system response pulse width T_TWith approximate System response pulse Width T'_TInstead, the estimated sampling rate of the acquisition card is obtained by using the expression. The resulting sampling rate SR requires more than or equal to twice the system bandwidth B.

According to the Nyquist sampling theorem, the signal sampled at the Nyquist rate or higher can ensure that the real water body signal can be accurately inverted. In a laser radar system, if the sampling rate of the acquisition card is too low, this will result inAnd if the sampling time is longer than the total system pulse width, the laser radar system undersamples the echo signal, which results in that the laser radar system cannot detect. As known from the equation, the required sampling rate of the acquisition card can be roughly estimated by the system response pulse width. For estimating the system sampling rate, the true value T_TAnd approximate model T'_TThe corresponding results are nearly identical.

In the step (4), when the phytoplankton layer thickness is corrected, and when the system response pulse width is large, the thickness error of the phytoplankton layer is large, and in the step (4), the correction equation is as follows:

wherein Thk' is the corrected phytoplankton layer thickness; thk is the measured thickness of the phytoplankton layer; l is_TIs an approximate system response pulse length, which is defined by an approximate system response pulse width T'_TObtained by light speed matching.

The laser radar system and the parameters thereof have certain broadening effect on the detection of the phytoplankton layer. The broadening effect is mainly shown in that the thickness of the phytoplankton layer in the echo signal is broadened, and the thickness of the detected phytoplankton layer is increased along with the increase of the pulse width of the laser and the response of the detector. The smaller the system response pulse width is, the greater the detection accuracy of the phytoplankton layer is. The accuracy achieved by the phytoplankton layer can be estimated from the total time response or matched as required by the required phytoplankton layer accuracy.

In the invention, the time domain waveform of the laser is Gaussian, and the time domain waveform of the detector is heavy-tail.

Compared with the prior art, the invention has the following beneficial effects;

1. the invention optimizes the device selection and calculation method of the marine laser radar system response, which comprises the following steps: a system response approximate model is provided, and the response pulse width of a laser radar system can be rapidly obtained; the matching model selection method among the three devices of the laser, the detector and the acquisition card is provided, and the optimal performance of different devices can be exerted;

2. the invention carries out suppression treatment on the error caused by the response of the marine laser radar system: the system response pulse width suitable for different water bodies is selected to reduce a water surface distortion area caused by water surface reflection and system widening; and correcting the measurement error of the phytoplankton layer caused by the system response through a correction equation.

3. According to the invention, the response of the marine laser radar system is optimized, so that the device selection in the process of developing the marine laser radar can be facilitated, and the influence of the system response on the inversion of the optical characteristics of the water body is reduced.

Drawings

FIG. 1 is a schematic flow chart of a response optimization processing method of a marine laser radar system according to the present invention;

FIG. 2 is a diagram illustrating the relationship between the relative error of the approximate model and the laser and the detector in the embodiment of the invention;

FIG. 3 is a diagram showing the relationship between the sampling rate of the lidar and the system time response in an embodiment of the present invention, in which (a) the matching condition of the acquisition card with the laser and the detector is provided, and (b) T is shown_TA relationship with SR;

FIG. 4 is a diagram illustrating response pulse widths of laser radar systems corresponding to different water types in an embodiment of the present disclosure;

FIG. 5 is a graph showing the relationship between the relative error of the thickness of the phytoplankton layer and the actual thickness and the system response pulse width, wherein (a) is the relative error of the measured phytoplankton layer thickness and the actual thickness, and L_TThe relationship (b) is the relative error distribution between the corrected phytoplankton layer thickness and the true value.

Detailed Description

The invention will be described in further detail below with reference to the drawings and examples, which are intended to facilitate the understanding of the invention without limiting it in any way.

As shown in fig. 1, a response optimization processing method for a marine laser radar system is based on a marine laser radar working in a blue-green light band, and information in a near-water surface and a phytoplankton layer is accurately obtained by adopting the response optimization processing method for the marine laser radar system.

Firstly, a device selection and calculation method for marine laser radar system response is optimized: according to the contribution of three key devices such as a laser, a detector and an acquisition card to system response, a system response approximate model is provided, and the response pulse width of a laser radar system is rapidly acquired; a matching model selection method among a laser device, a detector and a collection card is provided, and the optimal performance of different devices is exerted; secondly, the error caused by the response of the marine laser radar system is restrained: according to a quasi-single scattering approximate laser radar equation, a water surface distortion area caused by water surface reflection and system broadening is reduced by selecting system response pulse widths suitable for different water bodies; and correcting the measurement error of the phytoplankton layer caused by the system response through a correction equation.

The shape model of the outgoing laser pulse related by the invention is generally a Gaussian model, a heavy-tail model or a piecewise parabolic model. The gaussian model used in this embodiment can be expressed by formula (1):

the falling edge time of the photomultiplier is the rising edge time T _r2 to 3 times of (1), pulse width T_DIs about 2.5 times of the rising edge. The time response function of the PMT can be modeled with a heavy-tailed model, assuming T_D＝3T_rThe PMT response function may be expressed as equation (2):

laser pulse width T_LThe selection range is 2-18 ns, PMT rising edge time is set to be within 0.6-2.2 ns, and corresponding detector pulse width T_DThe dynamic range of the laser radar signal is 1.8-6.6 ns, the analysis is carried out by adopting a control variable method, the normalization processing is carried out on the laser radar signal, and the dynamic range of the laser radar signal is set to be 3 common orders of magnitude. In the simulation model, the surface reflectance is set to be 0.02, and the diffusion attenuation is setCoefficient of reduction K_dGreater than or equal to 0.1m^-1Atmospheric loss is 1, and the solar radiation reflected by the water surface is 0.14mw/m²srnm. The lidar system parameters used in the detailed simulation are shown in table 1.

TABLE 1

The laser radar system response approximate model used in the invention analyzes the broadening influence of three different laser waveforms on the system response, comprehensively considers, and adopts a Gaussian model because the Gaussian model conforms to the condition of most laser waveforms. T for laser pulse width_LIndicating the pulse width of the detector by T_DExpressing, according to the characteristics of the convolution of the gaussian function, an approximate model of the pulse width of the proposed system response can be expressed as formula (3):

wherein T is_DFor detector pulse width, T_LIs the laser pulse width. System response pulse width T_TIs the result of the interaction of the laser and the detector. Approximation model T'_TWith true value T_TThe error between can be written as equation (4):

FIG. 2 shows an approximate model T 'of the calculated system response for different laser and detector pulse widths'_TAnd true value T_TThe relative error between the pulse width and the pulse width of the laser is within 3 percent, the relative error distribution of the approximate model depends on the absolute difference value of the pulse width of the laser and the pulse width of the detector, and when the absolute difference value is larger, the error is smaller.

The bandwidth matching method of the system device is characterized in thatThe sampling rate SR of the acquisition card must be greater than or equal to twice the system bandwidth B (full width at half maximum in the frequency domain) to detect the signal normally. Carrying out linear fitting on the data to obtain T_TThe relationship with SR can be expressed by the following formula (5):

SR＝0.8234/T_T+3.3×10⁶(5)

fig. 3 shows the relationship between the lidar system sampling rate and the system response components.

In fig. 3, (a) is the relationship between the system response and the pulse width of the detector and the laser pulse width, the sampling rate SR of the acquisition card must be greater than or equal to two times the system bandwidth B (full width half maximum value in the frequency domain) to detect the signal normally. The figure provides the matching condition of the acquisition card, the laser and the detector, and the dotted line represents the corresponding laser pulse width and the detector pulse width when the sampling rate of the system is 120MHz, so that the acquisition card of 120MHz meets the system configuration under most conditions.

In FIG. 3, (b) shows T_TThe relation between SR, T can be seen_THas an excellent inverse proportional linear relationship with SR. Wherein the solid line corresponds to the result of the approximate model calculation and the dotted line corresponds to the result of the true value calculation, it can be seen that for the estimated system sampling rate, T_T、T′_TThe corresponding results are nearly identical.

The invention is based on the system response pulse width selection of the water turbidity, and the appropriate laser radar system response pulse width needs to be selected for different water types. Water diffusion attenuation coefficient K_dMaximum allowable system response pulse width T_TSatisfies the relation (6):

as shown in FIG. 4, further discussing the relationship between the response of the analysis system and the water detection capability, assuming that the ratio of the signal distortion zone to the maximum detectable depth is 20%, the response pulse width T of the lidar system corresponding to different water types is explored_T. The circle line in the figure shows the situation of selecting any laser pulse width and detector pulse width, and the solid line in the figure shows the laserThe detector is fixed at 2ns, and the broken line in the figure indicates that the detector is fixed at 1.8 ns.

Analysis of fig. 4 shows that the response pulse widths of the laser radar systems corresponding to different water body types have a certain fluctuation range, and the response pulse width range of the system becomes smaller and smaller as the diffusion attenuation coefficient of the water body is continuously increased. The lower boundary of the fluctuation range is the case where the laser pulse width is fixed to the minimum value, and the upper boundary is the case where the detector is fixed to the minimum value. Because the time domain response waveforms of the laser and the detector are different (the time domain waveform of the laser is Gaussian type, and the time domain waveform of the detector is heavy-tailed type), when the response of the detector is taken as the main condition, the allowable system response pulse width for detecting different water bodies is smaller than that when the response of the laser is taken as the main condition.

The thickness of the phytoplankton layer is corrected, when the system response pulse width is large, the thickness error of the phytoplankton layer is large, and the thickness can pass through the actually measured thickness Thk and the system response pulse length L_TCan be represented by formula (7) by modification:

the system response influences the detection of the phytoplankton layer, mainly the thickness of the phytoplankton layer, and the system response pulse width T_TIn relation to the thickness of the phytoplankton layer detected, T_TThe smaller the size, the greater the accuracy of detection of the phytoplankton layer. Estimating the accuracy achieved by the phytoplankton layer from the total time response, or matching the required T by the required phytoplankton layer accuracy_T. Approximate model T 'of system response pulse width'_TThe same applies to the detection of evaluation of phytoplankton layers.

The effect of phytoplankton layer thickness error on the system response pulse width T is further illustrated and analyzed by fig. 5_TThe setting range is 2.7-18.6 ns, corresponding to the pulse length L in the water body_TThe range is 0.3-2.1 m, and the thickness of the floating plant layer is 0.2-4.2 m.

In FIG. 5, (a) shows the phases for the measured phytoplankton layer thicknessFor error and true thickness and L_TThe relationship between the thickness and L can be seen from the figure_TWhen the difference is larger, the relative error of the measured values will be smaller and smaller. The relative errors of 10%, 20% and 40% are fitted linearly, and when the relative errors are constant, the phytoplankton layer thickness and the total system pulse width T are obtained_TOr L_TSubstantially satisfying a linear relationship. When the system response pulse width is larger, the thickness error of the phytoplankton layer is larger.

In fig. 5, (b) shows the distribution of the relative error between the corrected phytoplankton layer thickness and the true value, and the average value of the corrected relative error is increased from 55% to 6%. The linear regression method is adopted to fit the cases with relative errors of 5%, 10% and 20%, when Thk is more than 0.5L_TWhen +0.5, the correction error is not more than 5%, when Thk is 0.4L_TIn the vicinity of +0.3, the relative error is slightly large, but the distribution is still around 10%. When Thk is less than 0.3L_TAt +0.1, the corrected thickness error is still greater than 20%. In summary, the phytoplankton layer correction formula of the present invention can calibrate the actually measured thickness of the phytoplankton layer in most cases.

According to the invention, the response of the marine laser radar system is optimized, so that the device selection in the process of developing the marine laser radar can be facilitated, and the influence of the system response on the inversion of the optical characteristics of the water body is reduced.

The embodiments described above are intended to illustrate the technical solutions and advantages of the present invention, and it should be understood that the above-mentioned embodiments are only specific embodiments of the present invention, and are not intended to limit the present invention, and any modifications, additions and equivalents made within the scope of the principles of the present invention should be included in the scope of the present invention.