CN-122017846-A - Raman laser radar ozone three-dimensional real-time imaging method based on edge calculation

CN122017846ACN 122017846 ACN122017846 ACN 122017846ACN-122017846-A

Abstract

The invention relates to the technical field of radar detection, in particular to a Raman laser radar ozone three-dimensional real-time imaging method based on edge calculation, which not only can capture the influence of instantaneous temperature fluctuation on signals, but also can accurately evaluate the change of a continuous wind field on the air mass form by respectively and independently deducting the spectrum fidelity and topological deformation potential energy, thereby solving the contradiction of different physical processes on space-time scale; the compression and turbulence intensity of the fluid is directly quantized by using topological deformation potential energy, sensing data is mapped into a standard three-dimensional GIS grid by using a Lagrange locus backtracking algorithm and a Gaussian kernel function, the blind area between sensors is filled, each voxel data of the standard three-dimensional GIS grid is ensured to have a clear observation source, and the essential three-dimensional characteristics of the atmospheric environment are restored.

Inventors

DENG CHEN
YANG BIN
LI QIKUN
WU JIABIN

Assignees

南京新环光电科技有限公司

Dates

Publication Date: 20260512
Application Date: 20260122

Claims (7)

1. The Raman laser radar ozone three-dimensional real-time imaging method based on edge calculation is characterized by comprising the following specific steps of: S1, constructing a standard three-dimensional GIS grid based on an urban geographic information system, marking an ozone concentration detection value, and collecting the real-time photon number, the real-time temperature, the real-time air pressure, the real-time three-dimensional wind speed vector and the beam divergence angle of a laser emission system of a laser radar; S2, calculating a Raman scattering section correction coefficient by utilizing a second-order polynomial fitting formula according to the real-time photon number and the real-time temperature of the laser radar, and calculating the spectrum fidelity distributed along the distance; S3, calculating the change rate of the three-dimensional wind speed vector module length along with time according to the real-time three-dimensional wind speed vector to represent the air mass expansion strength, and simultaneously calculating the change rate of the three-dimensional wind speed vector vertical component along with time to represent the air mass vertical shear strength, so as to further obtain the real-time topological deformation potential energy of the wind field to the standard three-dimensional GIS grid; s4, acquiring a historical three-dimensional wind speed vector and historical topological deformation potential energy, calculating coordinates of a moving track point of the air mass based on a Lagrange track backtracking algorithm by taking a laser radar as a starting point, and calculating kinetic potential energy and physical fidelity aiming at any voxel in a standard three-dimensional GIS grid to further obtain a three-dimensional light-kinetic energy resolvable potential energy field; S5, outputting a preliminary ozone concentration field and calculating the concentration of the ozone vertical column by a differential absorption algorithm aiming at any voxel in the standard three-dimensional GIS grid, comparing the concentration of the ozone vertical column with the marked ozone concentration detection value, calculating the relative deviation, and outputting a final corrected three-dimensional ozone concentration grid matrix.
2. The Raman lidar ozone three-dimensional real-time imaging method based on edge calculation according to claim 1, wherein the S2 comprises the following specific steps: S21, calculating a Raman scattering section correction coefficient by using a second-order polynomial fitting formula based on real-time temperature, wherein the calculation formula of the Raman scattering section correction coefficient is as follows: ; Wherein, the The real-time temperature of the detection node at time K is indicated, in K, The standard temperature is used, the unit is K, The temperature coefficients of the rotating Raman scattering of the ozone molecules are respectively as follows The correction coefficient of the Raman scattering cross section at the moment k; s22, calculating the spectrum fidelity distributed along the distance, wherein the calculation formula of the spectrum fidelity is as follows: ; Wherein, the Representing the spectral fidelity of the r-th range bin in the line-of-sight direction of the lidar at time k, The sensitivity coefficient is a value of1, Representing the number of photons of the r-th distance element in the line-of-sight direction of the laser radar at time k, Is constant and takes the value of Is a normalization constant.
3. The edge-calculation-based Raman lidar ozone three-dimensional real-time imaging method according to claim 2, wherein the step S3 comprises the following specific steps: s31, calculating the change rate of the three-dimensional wind speed vector module length along with time according to the real-time three-dimensional wind speed vector to represent the air mass expansion strength, wherein the calculation formula is as follows: ; Wherein, the The air mass expansion intensity at time k is shown, For the three-dimensional wind speed vector modulo length at time k, For the three-dimensional wind speed vector modulo length at time k-1, Is the time interval between time k and time k-1; S32, calculating the change rate of the vertical component of the three-dimensional wind speed vector along with time according to the real-time three-dimensional wind speed vector to represent the vertical shear strength of the air mass, wherein the calculation formula is as follows: ; Wherein, the The vertical shear strength of the air mass at time k is indicated, For the three-dimensional wind speed vector vertical component of time k, The vertical component of the three-dimensional wind speed vector at the moment k-1; s33, calculating the real-time topological deformation potential energy of the wind field to the three-dimensional space grid based on the air mass expansion strength and the air mass vertical shear strength, wherein the calculation formula is as follows: ; Wherein, the For the topological deformation potential energy of the wind field to the standard three-dimensional GIS grid at the moment k, Is a weight factor of the air mass expansion intensity, Is a weight factor for the vertical shear strength of the air mass.
4. The method for three-dimensional real-time imaging of Raman lidar ozone based on edge calculation according to claim 3, wherein the step S4 comprises the following specific steps: S41, acquiring a historical three-dimensional wind speed vector and historical topological deformation potential energy, calculating coordinates of a moving track point of the air mass based on a Lagrange track backtracking algorithm by taking a laser radar as a starting point, wherein the calculation formula is as follows: ; Wherein, the For the trace point coordinates of the air mass at time k of the historical time k-m, The three-dimensional coordinates are fixed for the lidar, The three-dimensional wind speed vector of the historical moment k-n is the value of n which is 0-m; s42, aiming at any voxel in the standard three-dimensional GIS grid Obtaining center coordinates The kinetic potential energy is calculated by an inverse distance weighting method, and the calculation formula is as follows: ; Wherein, the Is the voxel index of a standard three-dimensional GIS grid, Is the coordinates of the voxel center Is used for the dynamic potential energy of the (a), Representing euclidean distance norms; S43, aiming at any voxel in standard three-dimensional GIS grid Obtaining center coordinates Projection length on laser axis Simultaneously acquiring central coordinates Perpendicular distance to laser axis Voxel calculation using gaussian kernel function The physical fidelity of (2) is calculated by the following formula: ; Wherein, the Is a voxel Is used for the physical fidelity of the (c) to be measured, For time k voxel Corresponding distance unit Is used for the optical spectrum fidelity of (c) in the optical spectrum, Divergence angle of the laser emission system beam; s44, calculating a three-dimensional light-kinetic energy resolvable potential energy field based on kinetic potential energy and physical fidelity, wherein the formula is as follows: 。
5. The edge-calculation-based Raman lidar ozone three-dimensional real-time imaging method according to claim 4, wherein the step S5 comprises the following specific steps: s51, aiming at any voxel in standard three-dimensional GIS grid Acquiring a distance unit in the laser radar sight line direction, mapping the distance unit to a corresponding height in a standard three-dimensional GIS grid, and outputting a preliminary ozone concentration field through a differential absorption algorithm based on the laser radar real-time photon number, real-time temperature and real-time air pressure of the corresponding height Calculating The concentration of the ozone vertical column is calculated according to the formula: ; Wherein, the Is that The concentration of the ozone in the vertical column, Is a grid layer in a standard three-dimensional GIS grid Is defined by the vertical thickness of the substrate, The maximum layer number in the vertical direction of the standard three-dimensional GIS grid; S52, vertical column concentration of ozone And marked ozone concentration detection value Comparing, and calculating relative deviation, wherein the calculation formula of the relative deviation is as follows: ; s53, based on preliminary ozone concentration field Three-dimensional photo-kinetic energy distinguishable potential energy field Relative deviation Outputting a final corrected three-dimensional ozone concentration grid matrix, wherein the calculation formula is as follows: ; Wherein, the Is the voxels in the standard three-dimensional GIS grid The final corrected ozone concentration is then provided to the process, Is a grid layer in a standard three-dimensional GIS grid Is a weight of (2).
6. A computer readable storage medium having stored thereon a computer program which, when executed by a processor, implements the raman lidar ozone three-dimensional real-time imaging method based on edge calculation according to any one of claims 1 to 5.
7. An electronic device comprising a memory for storing instructions and a processor for executing the instructions to cause the device to perform a method for implementing the edge-calculation based raman lidar ozone three-dimensional real-time imaging of any one of claims 1 to 5.

Description

Raman laser radar ozone three-dimensional real-time imaging method based on edge calculation Technical Field The invention relates to the technical field of radar detection, in particular to a Raman laser radar ozone three-dimensional real-time imaging method based on edge calculation. Background Raman lidar is currently the mainstream technology for detecting the vertical distribution of atmospheric ozone. The basic principle is that the concentration of ozone is inverted by utilizing a Raman scattering echo signal generated by the interaction of emitted laser and atmospheric molecules and combining a radar equation. In order to realize large-scale monitoring, the prior art trend is to deploy single-point laser radars on a mobile platform or build a distributed monitoring network and build a three-dimensional imaging system by combining a geographic information system. Although the prior laser radar three-dimensional imaging technology has advanced to a certain extent, the prior art still has the following fundamental technical defects when facing to urban complex terrains and changeable meteorological environments, and the imaging result is physically distorted, wherein the prior art generally adopts static atmosphere assumption or simple linear interpolation to process meteorological data, whereas the microscopic physical process (temperature fluctuation changes a Raman scattering section) is transient and high-frequency, while the macroscopic dynamic process (wind field driving air mass diffusion) is continuous and has space hysteresis, the prior art cannot decouple the two processes with quite different scales in a unified frame, the concentration inversion error is caused by ignoring the temperature transient fluctuation or the position dislocation of an imaging grid and an actual air mass in the fast-changing environment, the traditional GIS grid division is generally static or is only adjusted based on simple concentration gradient, however, the air is fluid, the ozone air mass can be stretched, compressed and distorted under the wind field effect, the prior art lacks real-time perception of the dynamic characteristics of the fluid, the prior art cannot cause the fault deformation in a strong area, the fixed and the effective coverage of the imaging result cannot occur. Disclosure of Invention The invention mainly aims to provide a Raman laser radar ozone three-dimensional real-time imaging method based on edge calculation, which can not only capture the influence of instantaneous temperature fluctuation on signals, but also accurately evaluate the change of continuous wind fields on air mass forms by respectively and independently deducting the optical fidelity and topological deformation potential energy, solve the contradiction of different physical processes on space-time scale, directly quantize the compression and turbulence intensity of fluid by utilizing the topological deformation potential energy, map sensing data into a standard three-dimensional GIS grid by utilizing a Lagrangian track backtracking algorithm and a Gaussian kernel function, fill the blind area among sensors, ensure that each voxel data of the standard three-dimensional GIS grid has a definite observation source, and restore the essential three-dimensional characteristics of the atmosphere environment. The technical scheme of the invention is as follows: in a first aspect, a raman lidar ozone three-dimensional real-time imaging method based on edge calculation is provided, and the method comprises the following steps: S1, constructing a standard three-dimensional GIS grid based on an urban geographic information system, marking an ozone concentration detection value, and collecting the real-time photon number, the real-time temperature, the real-time air pressure, the real-time three-dimensional wind speed vector and the beam divergence angle of a laser emission system of a laser radar; S2, calculating a Raman scattering section correction coefficient by utilizing a second-order polynomial fitting formula according to the real-time photon number and the real-time temperature of the laser radar, and calculating the spectrum fidelity distributed along the distance; S3, calculating the change rate of the three-dimensional wind speed vector module length along with time according to the real-time three-dimensional wind speed vector to represent the air mass expansion strength, and simultaneously calculating the change rate of the three-dimensional wind speed vector vertical component along with time to represent the air mass vertical shear strength, so as to further obtain the real-time topological deformation potential energy of the wind field to the standard three-dimensional GIS grid; s4, acquiring a historical three-dimensional wind speed vector and historical topological deformation potential energy, calculating coordinates of a moving track point of the air mass based on a Lagrange track backtracking algorithm by taking