CN113777627A - Laser radar ratio region transmission method based on novel foundation laser radar network - Google Patents

Laser radar ratio region transmission method based on novel foundation laser radar network Download PDF

Info

Publication number
CN113777627A
CN113777627A CN202111344183.5A CN202111344183A CN113777627A CN 113777627 A CN113777627 A CN 113777627A CN 202111344183 A CN202111344183 A CN 202111344183A CN 113777627 A CN113777627 A CN 113777627A
Authority
CN
China
Prior art keywords
laser radar
aerosol
ratio
lidar
content
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202111344183.5A
Other languages
Chinese (zh)
Other versions
CN113777627B (en
Inventor
刘�东
童奕澄
陈斯婕
肖达
张凯
方菁
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202111344183.5A priority Critical patent/CN113777627B/en
Publication of CN113777627A publication Critical patent/CN113777627A/en
Application granted granted Critical
Publication of CN113777627B publication Critical patent/CN113777627B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88—Lidar systems specially adapted for specific applications
    • G01S17/95—Lidar systems specially adapted for specific applications for meteorological use
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Electromagnetism (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

The invention discloses a laser radar ratio regional transmission method based on a novel foundation laser radar network, which comprises the steps of calculating the laser radar ratio and the content of aerosol corresponding to a selected position; screening the laser radar ratio and content data of the corresponding aerosol by combining historical data; establishing a nonlinear regression model between the ratio and the content of the laser radar corresponding to the aerosol; calculating the content of the corresponding aerosol observed by all the sunlight meters around the selected position; establishing a mathematical relation between the transmission distance and the laser radar relative error; and determining the application range of the laser radar ratio regional transmission method. By using the invention, the advanced laser radar can realize the transmission of the laser radar ratio region to the Mi scattering laser radar, and the system cost can be controlled while the accuracy of the observation data of the laser radar network is ensured.

Description

Laser radar ratio region transmission method based on novel foundation laser radar network
Technical Field
The invention belongs to the technical field of foundation atmosphere laser radar networking, and particularly relates to a laser radar ratio area transmission method based on a novel foundation laser radar network.
Background
The laser radar has the advantages of high space-time resolution, high precision, continuous monitoring and the like, and is widely applied to aerosol observation. The ground-based laser radar networking observation has important significance for developing accurate research of multi-region aerosol parameters and comprehensively mastering the nature of atmospheric motion change and the time-space distribution evolution rule among different regions. Currently, a plurality of ground-based laser radar networks are established internationally, such as European Aerosol Research laser radar networks (EARLINET), Asian Dust laser radar Observation networks (AD-Net), Micro-Pulse laser radar networks (MPLNET), and the like. The laser radar network represented by EARLINET is mainly based on a Raman laser radar system, and the laser radar of the type can realize accurate inversion of aerosol optical characteristic parameters without assuming the ratio of extinction coefficient to backscattering coefficient (also called laser radar ratio).
In addition, high spectral resolution lidar can also achieve the same function as raman lidar, which can work all day long, whereas raman lidar can only work at night. However, advanced lidar systems like high-spectral resolution lidar and raman lidar have complex structures and higher cost, and are not easy to realize large-scale deployment. The laser radar networks represented by AD-Net and MPLNET mainly use a Mie scattering laser radar system, and the laser radar of the type needs to assume a laser radar ratio when inverting optical characteristic parameters of aerosol, which results in lower accuracy of data inversion. However, the Mie scattering lidar system has a simple structure and low cost, and is easy to realize large-scale deployment.
In summary, the problems of sparse arrangement density of advanced lidar systems and low accuracy of data observed by the Mie scattering lidar exist in the existing lidar networking, so that a technology for realizing parameter accuracy transmission from the advanced lidar to the Mie scattering lidar needs to be developed urgently. The existing direct parameter transmission technology can be realized by tracing through a Hybrid Single-Particle Integrated track (HYSPLIT) model, which requires that an aerosol group observed by an advanced laser radar just moves to the upper part of a Mie scattering laser radar, however, the situations are few, the limitation is large, and the observation requirement of a large-scale laser radar network is difficult to meet.
Disclosure of Invention
The invention provides a laser radar ratio region transmission method based on a novel foundation laser radar network, which can realize the balance between the precision of observation data of the laser radar network and the system cost on the premise of not increasing the cost of additional instruments.
A laser radar ratio area transmission method based on a novel foundation laser radar net comprises a high-spectral-resolution laser radar module arranged at the center and a plurality of meter-scattering laser radar modules arranged around the high-spectral-resolution laser radar module; the distance between each millimeter scattering laser radar module and the high spectral resolution laser radar module is different;
the high-spectral-resolution laser radar module comprises a high-spectral-resolution laser radar and a corresponding sunshine photometer which are cooperatively observed at the same position;
the meter scattering laser radar module comprises a meter scattering laser radar and a corresponding sunshine photometer which are arranged at the same position and used for cooperative observation;
the laser radar specific area transmission method comprises the following steps:
(1) selecting an O position and an aerosol type of an area to be researched, and arranging a novel foundation laser radar net by taking the O position as a center; respectively calculating the aerosol laser radar ratio observed by the high-spectral-resolution laser radar and the aerosol content observed by the corresponding sunshine photometer at the O position at the same moment;
(2) screening the aerosol laser radar ratio obtained in the step (1) and the data of the corresponding aerosol content by combining the historical data of the aerosol at the O position;
(3) after data screening, establishing a nonlinear regression model between the aerosol laser radar ratio at the O position and corresponding aerosol content data;
(4) calculating the content of aerosol observed by all sunlight meters within 500km around the O position by taking the O position as a center;
(5) substituting the aerosol content calculated in the step (4) into the nonlinear regression model in the step (3), and establishing a mathematical relation between relative distances (transmission distances) between a plurality of groups of different sunlight meter positions and a laser radar comparison relative error;
(6) and (5) determining the upper limit of the relative error of the laser radar ratio, and determining the transmission range applicable to the laser radar ratio regional transmission method by combining the established mathematical relation in the step (5).
In the step (1), the same time refers to data that the aerosol laser radar ratio and the aerosol content need to be selected at the integral time (0, 1 … 22, 23) of each day.
The calculation process of the aerosol laser radar ratio is as follows:
the high spectral resolution lidar obtains lidar equations for three channels:
Figure 701168DEST_PATH_IMAGE001
Figure 165779DEST_PATH_IMAGE002
Figure 617620DEST_PATH_IMAGE003
in the formula (I), the compound is shown in the specification,
Figure 684933DEST_PATH_IMAGE004
indicating that the lidar attenuated backscatter signal after system correction,
Figure 436988DEST_PATH_IMAGE005
=1, 2 or 3;
Figure 931555DEST_PATH_IMAGE006
representing the vertical channel backscatter coefficient of aerosol or atmospheric molecules,
Figure 136271DEST_PATH_IMAGE007
represents vertical;
Figure 7275DEST_PATH_IMAGE008
representing the parallel channel backscatter coefficients of aerosol or atmospheric molecules,
Figure 817099DEST_PATH_IMAGE009
which represents a parallel,
Figure 748146DEST_PATH_IMAGE010
and
Figure 440159DEST_PATH_IMAGE011
respectively representing aerosol and atmospheric molecules;
Figure 849274DEST_PATH_IMAGE012
and
Figure 44764DEST_PATH_IMAGE013
respectively representing the detected initial height and the detected end height;
Figure 146712DEST_PATH_IMAGE014
represents the extinction coefficient of an aerosol or atmospheric molecule;
Figure 794862DEST_PATH_IMAGE015
represents the transmittance of an aerosol or molecular channel;
through the laser radar equation of the three channels
Figure 742089DEST_PATH_IMAGE016
、
Figure 588823DEST_PATH_IMAGE017
、
Figure 861672DEST_PATH_IMAGE018
And
Figure 528277DEST_PATH_IMAGE019
the aerosol laser radar ratio observed by the high spectral resolution laser radar is calculated by the following formula:
Figure 279195DEST_PATH_IMAGE020
。
the calculation of the aerosol content is as follows:
the aerosol content is calculated by adopting an observation parameter inversion formula based on a sunshine photometer, which is proposed by Bahadur:
Figure 714856DEST_PATH_IMAGE021
in the formula (I), the compound is shown in the specification,
Figure 158606DEST_PATH_IMAGE022
represents the optical thickness of the absorbed aerosol, and can be obtained from the observation data of a solar photometer;
Figure 46928DEST_PATH_IMAGE023
which represents the wavelength of the light emitted by the light source,
Figure 867116DEST_PATH_IMAGE005
440, 675 or 880 nm;
Figure 157283DEST_PATH_IMAGE024
which represents the reference wavelength(s) of the light,
Figure 771936DEST_PATH_IMAGE025
440, 675 or 880 nm;
Figure 616395DEST_PATH_IMAGE026
representing the optical thickness of the absorbing aerosol based on different aerosol contents at a reference wavelength,
Figure 709116DEST_PATH_IMAGE027
is black carbon
Figure 119368DEST_PATH_IMAGE028
Brown carbon
Figure 170501DEST_PATH_IMAGE029
Or dust and sand
Figure 33415DEST_PATH_IMAGE030
；
Figure 460985DEST_PATH_IMAGE031
Which represents the absorption of the light in the angstrom index,
Figure 194586DEST_PATH_IMAGE032
=1 or 2;
the optical aerosol thickness of the aerosol components is obtained by the following formula:
Figure 416620DEST_PATH_IMAGE033
in the formula (I), the compound is shown in the specification,
Figure 501251DEST_PATH_IMAGE034
represents the optical thickness of the aerosol of black carbon, brown carbon or sand dust;
Figure 201353DEST_PATH_IMAGE035
representing the single scattering albedo of black carbon, brown carbon or sand dust; by calculation of
Figure 258302DEST_PATH_IMAGE034
Account for
Figure 385658DEST_PATH_IMAGE036
The content of the corresponding aerosol can be calculated and obtained.
In the step (2), the screening conditions for the aerosol laser radar ratio comprise height, backscattering ratio, aerosol depolarization ratio and laser radar ratioA range of (d); wherein the backscattering ratio
Figure 668568DEST_PATH_IMAGE037
Deviation ratio of aerosol
Figure 172362DEST_PATH_IMAGE038
The formula of (1) is as follows:
Figure 146134DEST_PATH_IMAGE039
Figure 444391DEST_PATH_IMAGE040
in the formula (I), the compound is shown in the specification,
Figure 769193DEST_PATH_IMAGE041
and
Figure 76678DEST_PATH_IMAGE042
respectively representing the depolarization ratio of the aerosol and the depolarization ratio of atmospheric molecules; the laser radar ratio value is an average value of data with the height of 0.5-4 km after screening by using the conditions of a backscattering ratio, an aerosol depolarization ratio and a laser radar ratio;
the screening condition for the aerosol content is to screen based on the month in which the O position was explored to have a higher aerosol content as a lower limit.
The specific process of the step (5) is as follows:
(5-1) optionally selecting the two sites of the sun photometer mentioned in the step (4), wherein one site is set as an M site, and the other site is set as an N site; converting the content of the aerosol observed by the two stations into a laser radar ratio through the nonlinear regression model in the step (3);
(5-2) assuming that the M site is taken as the center, and the longitude and latitude of the M point and the N point are respectively
Figure 170536DEST_PATH_IMAGE043
And
Figure 639695DEST_PATH_IMAGE044
the relative distance between the two stations is calculated M, N, as follows:
Figure DEST_PATH_IMAGE046AAAA
in the formula (I), the compound is shown in the specification,
Figure 530421DEST_PATH_IMAGE047
representing the radius of the earth, respectively selecting the data of aerosol contents of two stations in the same hour, wherein the data is the data of the whole point time of each day, and the laser radar ratio of the M stations is used as the true value
Figure 376018DEST_PATH_IMAGE048
The laser radar ratio of N sites is a reference value
Figure 589961DEST_PATH_IMAGE049
Calculating the relative error of the laser radar ratio; the laser radar relative error and the transmission distance of the two stations are called a group of corresponding points;
(5-3) repeating the step (5-2), and calculating the laser radar relative error and relative distance between any two stations to obtain a plurality of groups of corresponding points;
and (5-4) establishing a mathematical relation between the laser radar specific relative error and the relative distance according to the plurality of groups of corresponding points calculated in the step (5-3).
In the step (6), an error transfer relation between a backscattering coefficient, an extinction coefficient and an aerosol laser radar ratio needs to be established, and the formula is as follows:
Figure 964442DEST_PATH_IMAGE050
in the formula (I), the compound is shown in the specification,
Figure 732678DEST_PATH_IMAGE051
indicating the total relative error of the aerosol lidar ratio,
Figure 647544DEST_PATH_IMAGE052
representing the statistical standard deviation of the corresponding quantity; order to
Figure 450415DEST_PATH_IMAGE053
And
Figure 261376DEST_PATH_IMAGE054
respectively is
Figure 985750DEST_PATH_IMAGE055
And
Figure 173149DEST_PATH_IMAGE056
then the total relative error of the aerosol lidar ratio is solved as
Figure 96105DEST_PATH_IMAGE057
(ii) a And (5) determining the applicable transmission range of the laser radar ratio region transmission method by combining the upper limit of the aerosol laser radar ratio relative error with the mathematical relation obtained in the step (5).
In addition, the specific number (more than or equal to 1) of the Mi scattering laser radars can be arranged according to the station positions and the number of the sunshine photometers existing in the range of 500km by taking the high-spectral-resolution laser radars as the center.
Compared with the prior art, the invention has the following beneficial effects:
1. the novel laser radar network provided by the invention realizes the balance between the observation data precision and the system cost of the laser radar network on the premise of not increasing the cost of additional instruments, fully exerts the respective advantages of advanced laser radars and meter scattering laser radars, and can promote the large-scale deployment of the laser radar network in China;
2. the method has simple thought, provides an important thought for realizing the precision transmission of the parameters between the laser radar networks, and has strong popularization value.
Drawings
Fig. 1 is a schematic layout of the novel ground-based lidar network of the present invention.
Fig. 2 is a schematic flow chart of a laser radar ratio region transmission method in the present invention.
Fig. 3 is a non-linear regression model between aerosol lidar ratios and corresponding contents established in an embodiment of the present disclosure.
Fig. 4 is a mathematical relationship between the aerosol lidar relative error and the transmission distance established in the embodiment of the present invention.
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, the novel ground-based lidar network in this embodiment includes a high spectral resolution lidar module 1 and a plurality of millimeter-scattering lidar modules 2. The high spectral resolution laser radar module 1 mainly comprises a combination of a high spectral resolution laser radar and a sun photometer, and the meter scattering laser radar module 2 mainly comprises a combination of a meter scattering laser radar and a sun photometer. High spectral resolution lidar module 1 is arranged at the central position, meter scattering lidar module 2 is arranged at the surrounding position, and the distances between different meter scattering lidar modules 2 and high spectral resolution lidar module 1 are different.
As shown in fig. 2, a lidar specific area transmission method based on a novel ground-based lidar network includes the following steps:
s1: the place selected in this example is seoul university, seoul, korea, whose latitude and longitude is (126.95E, 37.56N). The high spectral resolution laser radar is selected as 'AHSRL' developed by Wisconsin university. The time period of the data is selected to be 2016.1.1-2018.12.31, the time-space resolution of the data is 1 hour (hour of the whole point of the day, 0, 1 … 22, 23 hours) and 7.5m, and the height range is 0-10 km. The AHSRL high spectral resolution lidar can obtain lidar equations for three channels:
Figure 343547DEST_PATH_IMAGE001
Figure 351954DEST_PATH_IMAGE002
Figure 608623DEST_PATH_IMAGE003
in the formula (I), the compound is shown in the specification,
Figure 386086DEST_PATH_IMAGE004
indicating that the lidar attenuated backscatter signal after system correction,i=1, 2 or 3;
Figure 70009DEST_PATH_IMAGE006
representing the vertical channel backscatter coefficient of aerosol or atmospheric molecules,
Figure 300133DEST_PATH_IMAGE007
represents vertical;
Figure 94913DEST_PATH_IMAGE058
representing the parallel channel backscatter coefficients of aerosol or atmospheric molecules,
Figure 195725DEST_PATH_IMAGE009
which represents a parallel,
Figure 784969DEST_PATH_IMAGE010
and
Figure 33548DEST_PATH_IMAGE011
respectively representing aerosol and atmospheric molecules;
Figure 100861DEST_PATH_IMAGE012
and
Figure 852916DEST_PATH_IMAGE013
respectively representing the detected initial height and the detected end height;
Figure 878641DEST_PATH_IMAGE014
to representExtinction coefficient of aerosol or atmospheric molecules;
Figure 348937DEST_PATH_IMAGE015
represents the transmittance of an aerosol or molecular channel;
through the laser radar equation of three channels can be solved
Figure 485520DEST_PATH_IMAGE016
、
Figure 826503DEST_PATH_IMAGE059
、
Figure 288708DEST_PATH_IMAGE018
And
Figure 980720DEST_PATH_IMAGE060
the aerosol laser radar ratio observed by the AHSRL high spectral resolution laser radar can be obtained by the following calculation:
Figure 655415DEST_PATH_IMAGE061
。
in this embodiment, the AERONET site name is selected as "Yonseii University". The solar photometer model is CE318, which is manufactured by CIMEL electronics of france. The selected time period is 2016.1.1-2018.12.31. The data selected were 1.5 levels of AERONET version 3 and were averages at each hour time (0, 1 … 22, 23) which had to be the same time as the selection of the lidar ratio described above. The aerosol content is calculated by adopting an inversion formula based on observation parameters of a sunshine photometer, which is proposed by American scientist Bahadur:
Figure 585325DEST_PATH_IMAGE021
in the formula (I), the compound is shown in the specification,
Figure 952853DEST_PATH_IMAGE022
representing the optical thickness of the absorbed aerosol, which can be measured from a solar photometerObtaining;
Figure 866582DEST_PATH_IMAGE023
which represents the wavelength of the light emitted by the light source,
Figure 610547DEST_PATH_IMAGE005
440, 675 or 880 nm;
Figure 191701DEST_PATH_IMAGE062
which represents the reference wavelength(s) of the light,
Figure 198971DEST_PATH_IMAGE025
440, 675 or 880 nm;
Figure 131155DEST_PATH_IMAGE063
representing the optical thickness of the absorbing aerosol based on different aerosol contents at a reference wavelength,
Figure 616494DEST_PATH_IMAGE027
is Black carbon (Black carbon,
Figure 52155DEST_PATH_IMAGE028
) Brown carbon (Brown carbon,
Figure 495906DEST_PATH_IMAGE029
) Or dust (a)
Figure 649807DEST_PATH_IMAGE030
）；
Figure 735574DEST_PATH_IMAGE064
Which represents the absorption of the light in the angstrom index,
Figure 25741DEST_PATH_IMAGE032
=1 or 2; in the present embodiment, the first and second electrodes are,
Figure 905973DEST_PATH_IMAGE065
、
Figure 281590DEST_PATH_IMAGE066
and
Figure 639890DEST_PATH_IMAGE067
the values are respectively 0.55 +/-0.24, 4.55 +/-2.01 and 2.20 +/-0.50;
Figure 784564DEST_PATH_IMAGE068
、
Figure 304538DEST_PATH_IMAGE069
and
Figure 901873DEST_PATH_IMAGE070
the values are respectively 0.85 plus or minus 0.40, 0 and 1.15 plus or minus 0.50.
The optical aerosol thickness of the aerosol components is obtained by the following formula:
Figure 63864DEST_PATH_IMAGE033
in the formula (I), the compound is shown in the specification,
Figure 352060DEST_PATH_IMAGE034
represents the optical thickness of the aerosol of black carbon, brown carbon or sand dust;
Figure 574094DEST_PATH_IMAGE035
the values of the single scattering albedo of black carbon, brown carbon or sand dust are 0.480, 0.772 and 0.870 respectively in this embodiment.
By calculation of
Figure 189884DEST_PATH_IMAGE034
Account for
Figure 889986DEST_PATH_IMAGE036
The content of the corresponding aerosol can be calculated and obtained.
S2: in this embodiment, the sand aerosol is selected, and the screening conditions for the sand lidar ratio include the height, the backscattering ratio, the aerosol depolarization ratio, and the lidar ratio range. Wherein the backscattering ratio
Figure 274831DEST_PATH_IMAGE037
Deviation ratio of aerosol
Figure 933346DEST_PATH_IMAGE038
The formula of (1) is as follows:
Figure 770852DEST_PATH_IMAGE039
Figure 274645DEST_PATH_IMAGE040
in the formula (I), the compound is shown in the specification,
Figure 248418DEST_PATH_IMAGE041
and
Figure 812254DEST_PATH_IMAGE042
respectively representing the depolarization ratio of the aerosol and the depolarization ratio of the atmospheric molecules. Wherein
Figure 402635DEST_PATH_IMAGE041
Can be solved by radar equations of three channels of the high-spectral-resolution laser radar in the embodiment
Figure 444541DEST_PATH_IMAGE071
The value was 0.0143. The screening conditions were as follows: the height is 0.5-4 km; backscatter ratio: 1.2-10; aerosol depolarization ratio: 0.15 to 0.3; laser radar ratio: 30 to 90 Sr. The laser radar ratio value is an average value of data with the height of 0.5-4 km after screening conditions of a backscattering ratio, an aerosol depolarization ratio and a laser radar ratio.
The average value of the dust content of the Yangshan university in 3 and 4 months per year in the period of 2016.1.1-2018.12.31 is calculated, and the calculation result of the embodiment is 21.6%. For convenience, the screening condition of the present example for the content of sand dust is required to be more than 20%.
Through screening, the present embodiment obtains a value 400 group of the sand-dust lidar ratio and the sand-dust content.
S3: as shown in fig. 3, a non-linear regression model between the sand-dust lidar ratio at the O position and 400 groups of data with the screened corresponding content is established. The embodiment adopts a quadratic polynomial fitting method, the decision coefficient of the fitting is 0.78, and the fitting has strong correlation.
S4: considering the data volume problem, the present example selects six sites of maryland center (76.61W, 39.28N), einz germany (8.30E, 50.00N), carl pomtra france (5.06E, 44.08N), taihu china (120.22E, 31.42N), beijing china (116.38E, 39.98N), and korean university of extension (126.95E, 37.56N) as the center, the data selection time ranges from 2001.1.1 to 2021.3.31, and calculates the content of dust aerosol observed by all the sun photometers at every whole time (0, 1 … 22, 23) within 500km around each center site. In this embodiment, a total of 63 solar photometer sites are selected. In addition, as described in step S2, the present embodiment also requires the screening condition for the content of sand dust in step S4 to be more than 20%.
S5: the embodiment specifically implements the following process in this step:
s5-1: two solar photometric sites mentioned in the step S4 are arbitrarily selected, one of which is set as an M site, and the other is set as an N site. Converting the sand and dust content observed by the two stations into a laser radar ratio through the nonlinear regression model in the step S3;
s5-2: suppose that the M site is taken as the center, and the longitude and latitude of the M point and the N point are respectively
Figure 538399DEST_PATH_IMAGE043
And
Figure 538716DEST_PATH_IMAGE044
the relative distance between the two stations is calculated M, N, as follows:
Figure DEST_PATH_IMAGE046_5A
in the formula (I), the compound is shown in the specification,
Figure 429442DEST_PATH_IMAGE047
representing the radius of the earth. Selecting the data of aerosol contents of the two stations in the same hour (integral time of day, 0, 1 … 22, 23 hours) respectively, and taking the laser radar ratio of the M stations as a real value
Figure 540618DEST_PATH_IMAGE048
The laser radar ratio of N sites is a reference value
Figure 754562DEST_PATH_IMAGE049
The relative error of the lidar ratio is calculated. The relative distance between the laser radar phase comparison relative error and the two stations can be called a group of corresponding points;
s5-3: repeating the step S5-2, and calculating the laser radar relative error and the relative distance between any two stations so as to obtain 1953 groups of corresponding points;
s5-4: the embodiment requires that the total number of the calculated average lidar relative errors in each group of corresponding points is more than 50, and the distance between two stations is less than 500 km. After screening, 78 sets of valid corresponding point data are obtained. As shown in fig. 4, a mathematical relationship between the lidar relative error and the relative distance is established.
S6: the error transfer formulas of the backscattering coefficient, the extinction coefficient and the aerosol laser radar ratio are as follows:
Figure 660201DEST_PATH_IMAGE050
in the formula (I), the compound is shown in the specification,
Figure 959595DEST_PATH_IMAGE051
indicating the total relative error of the aerosol lidar ratio,
Figure 343303DEST_PATH_IMAGE052
representing the statistical standard deviation of the corresponding quantity; order to
Figure 677332DEST_PATH_IMAGE053
And
Figure 753873DEST_PATH_IMAGE054
respectively is
Figure 274984DEST_PATH_IMAGE055
And
Figure 462383DEST_PATH_IMAGE056
then the total relative error of the aerosol lidar ratio is solved as
Figure 385339DEST_PATH_IMAGE057
. The Russian scientist Igor et al thinks that when the method is used for inverting the micro-physical characteristics of the aerosol, the acceptable maximum relative errors of the backscattering coefficient and the extinction coefficient are both 20%. Therefore, this embodiment takes
Figure 898360DEST_PATH_IMAGE055
=
Figure 172347DEST_PATH_IMAGE056
=20%, the calculated maximum relative error of the lidar ratio was 28.3%. With reference to fig. 4, the range of the laser radar applicable to the regional transmission can be determined to be 0-128 km. By using the method provided by the invention in the range, the error of the inversion extinction coefficient of the Mie scattering laser radar can be within 20%. The specific number (more than or equal to 1) of the Mi scattering laser radars can be arranged according to the positions and the number of stations of the sunshine photometer existing in the range of 500km by taking the high-spectral-resolution laser radars as the center.
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.

Claims (6)

1. A laser radar ratio regional transmission method based on a novel foundation laser radar net is characterized in that the novel foundation laser radar net comprises a high-spectral-resolution laser radar module arranged at the center and a plurality of meter-scattering laser radar modules arranged around the high-spectral-resolution laser radar module; the distance between each millimeter scattering laser radar module and the high spectral resolution laser radar module is different;
the high-spectral-resolution laser radar module comprises a high-spectral-resolution laser radar and a corresponding sunshine photometer which are cooperatively observed at the same position;
the meter scattering laser radar module comprises a meter scattering laser radar and a corresponding sunshine photometer which are arranged at the same position and used for cooperative observation;
the laser radar specific area transmission method comprises the following steps:
(1) selecting an O position and an aerosol type to be researched, and arranging a novel foundation laser radar net by taking the O position as a center; respectively calculating aerosol laser radar ratios observed by the high-spectral-resolution laser radar at the positions O and aerosol content corresponding to the observation of the sun photometer at the same time;
(2) screening the aerosol laser radar ratio obtained in the step (1) and the data of the corresponding aerosol content by combining the historical data of the aerosol at the O position;
(3) after data screening, establishing a nonlinear regression model between the aerosol laser radar ratio at the O position and corresponding aerosol content data;
(4) calculating the content of aerosol observed by all sunlight meters within 500km around the O position by taking the O position as a center;
(5) substituting the aerosol content calculated in the step (4) into the nonlinear regression model in the step (3), and establishing a mathematical relation between transmission distances between a plurality of groups of different sunlight meter positions and a laser radar ratio relative error; the specific process is as follows:
(5-1) optionally selecting the two sites of the sun photometer mentioned in the step (4), wherein one site is set as an M site, and the other site is set as an N site; converting the content of the aerosol observed by the two stations into a laser radar ratio through the nonlinear regression model in the step (3);
(5-2) assuming that the longitude and latitude of M point and N point are respectively centered on M station
Figure 322812DEST_PATH_IMAGE001
And
Figure 801198DEST_PATH_IMAGE002
the relative distance between the two stations is calculated M, N, as follows:
Figure 116773DEST_PATH_IMAGE003
in the formula (I), the compound is shown in the specification,
Figure 655202DEST_PATH_IMAGE004
representing the radius of the earth, respectively selecting the data of aerosol contents of two stations in the same hour, wherein the data is the data of the whole point time of each day, and the laser radar ratio of the M stations is used as the true value
Figure 56227DEST_PATH_IMAGE005
The laser radar ratio of N sites is a reference value
Figure 603883DEST_PATH_IMAGE006
Calculating the relative error of the laser radar ratio; the laser radar relative error and the transmission distance of the two stations are called a group of corresponding points;
(5-3) repeating the step (5-2), and calculating the laser radar relative error and relative distance between any two stations to obtain a plurality of groups of corresponding points;
(5-4) establishing a mathematical relation between the laser radar relative error and the relative distance according to the plurality of groups of corresponding points calculated in the step (5-3);
(6) and (5) determining the upper limit of the relative error of the laser radar ratio, and determining the transmission range applicable to the laser radar ratio regional transmission method by combining the established mathematical relation in the step (5).
2. The lidar ratio regional transmission method based on the novel ground-based lidar network according to claim 1, wherein in the step (1), the same time refers to data that the aerosol lidar ratio and the aerosol content need to be selected at an integral time of day.
3. The lidar ratio regional transmission method based on the novel ground-based lidar network according to claim 1, wherein in the step (1), the calculation process of the aerosol lidar ratio is as follows:
the high spectral resolution lidar obtains lidar equations for three channels:
Figure 508385DEST_PATH_IMAGE007
Figure 217715DEST_PATH_IMAGE008
Figure 106037DEST_PATH_IMAGE009
in the formula (I), the compound is shown in the specification,
Figure 926225DEST_PATH_IMAGE010
indicating that the lidar attenuated backscatter signal after system correction,
Figure 950813DEST_PATH_IMAGE011
=1, 2 or 3;
Figure 831044DEST_PATH_IMAGE012
representing the vertical channel backscatter coefficient of aerosol or atmospheric molecules,
Figure 675504DEST_PATH_IMAGE013
represents vertical;
Figure 564962DEST_PATH_IMAGE014
representing the parallel channel backscatter coefficients of aerosol or atmospheric molecules,
Figure 709636DEST_PATH_IMAGE015
which represents a parallel,
Figure 760768DEST_PATH_IMAGE016
and
Figure 623682DEST_PATH_IMAGE017
respectively representing aerosol and atmospheric molecules;
Figure 785673DEST_PATH_IMAGE018
and
Figure 784853DEST_PATH_IMAGE019
respectively representing the detected initial height and the detected end height;
Figure 6887DEST_PATH_IMAGE020
represents the extinction coefficient of an aerosol or atmospheric molecule;
Figure 91518DEST_PATH_IMAGE021
represents the transmittance of an aerosol or molecular channel;
through the laser radar equation of the three channels
Figure 322779DEST_PATH_IMAGE022
、
Figure 442045DEST_PATH_IMAGE023
、
Figure 545963DEST_PATH_IMAGE024
And
Figure 117890DEST_PATH_IMAGE025
the aerosol laser radar ratio observed by the high spectral resolution laser radar is calculated by the following formula:
Figure 621684DEST_PATH_IMAGE026
。
4. the lidar network based on the novel ground-based lidar network according to claim 1, wherein in the step (1), the calculation process of the aerosol content is as follows:
the aerosol content is calculated by adopting an observation parameter inversion formula based on a sunshine photometer, which is proposed by Bahadur:
Figure 861035DEST_PATH_IMAGE027
in the formula (I), the compound is shown in the specification,
Figure 424872DEST_PATH_IMAGE028
represents the optical thickness of the absorbed aerosol, and can be obtained from the observation data of a solar photometer;
Figure 218515DEST_PATH_IMAGE029
which represents the wavelength of the light emitted by the light source,
Figure 260421DEST_PATH_IMAGE011
440, 675 or 880 nm;
Figure 354279DEST_PATH_IMAGE030
which represents the reference wavelength(s) of the light,
Figure 557858DEST_PATH_IMAGE031
440, 675 or 880 nm;
Figure 104377DEST_PATH_IMAGE032
representing the optical thickness of the absorbing aerosol based on different aerosol contents at a reference wavelength,
Figure 949973DEST_PATH_IMAGE033
is black carbon
Figure 632758DEST_PATH_IMAGE034
Brown carbon
Figure 538397DEST_PATH_IMAGE035
Or sand dust;
Figure 572213DEST_PATH_IMAGE036
which represents the absorption of the light in the angstrom index,
Figure 487079DEST_PATH_IMAGE037
=1 or 2;
the optical aerosol thickness of the aerosol components is obtained by the following formula:
Figure 24371DEST_PATH_IMAGE038
in the formula (I), the compound is shown in the specification,
Figure 100911DEST_PATH_IMAGE039
represents the optical thickness of the aerosol of black carbon, brown carbon or sand dust;
Figure 622022DEST_PATH_IMAGE040
representing the single scattering albedo of black carbon, brown carbon or sand dust; by calculation of
Figure 75000DEST_PATH_IMAGE039
Account for
Figure 997957DEST_PATH_IMAGE041
The content of the corresponding aerosol can be calculated and obtained.
5. The lidar ratio regional transmission method based on the novel ground-based lidar network according to claim 1, wherein in the step (2), the screening conditions for the aerosol lidar ratio comprise height, backscattering ratio, aerosol depolarization ratio and lidar ratio range; wherein the backscattering ratio
Figure 979819DEST_PATH_IMAGE042
Deviation ratio of aerosol
Figure 253806DEST_PATH_IMAGE043
The formula of (1) is as follows:
Figure 244896DEST_PATH_IMAGE044
Figure 22359DEST_PATH_IMAGE045
in the formula (I), the compound is shown in the specification,
Figure 440702DEST_PATH_IMAGE046
and
Figure 670826DEST_PATH_IMAGE047
respectively representing the depolarization ratio of the aerosol and the depolarization ratio of atmospheric molecules; the laser radar ratio value is an average value of data with the height of 0.5-4 km after screening by using the conditions of a backscattering ratio, an aerosol depolarization ratio and a laser radar ratio;
the screening condition for the aerosol content is to screen based on the month in which the O position was explored to have a higher aerosol content as a lower limit.
6. The lidar ratio regional transmission method based on the novel ground-based lidar network according to claim 1, wherein in the step (6), an error transfer relationship between a backscattering coefficient, an extinction coefficient and an aerosol lidar ratio is established, and the formula is as follows:
Figure 934448DEST_PATH_IMAGE048
in the formula (I), the compound is shown in the specification,
Figure 566418DEST_PATH_IMAGE049
indicating the total relative error of the aerosol lidar ratio,
Figure 421241DEST_PATH_IMAGE050
representing the statistical standard deviation of the corresponding quantity; order to
Figure 138662DEST_PATH_IMAGE051
And
Figure 471554DEST_PATH_IMAGE052
respectively is
Figure 223609DEST_PATH_IMAGE053
And
Figure 452597DEST_PATH_IMAGE054
then the total relative error of the aerosol lidar ratio is solved as
Figure 657313DEST_PATH_IMAGE055
(ii) a And (5) determining the applicable transmission range of the laser radar ratio region transmission method by combining the upper limit of the aerosol laser radar ratio relative error with the mathematical relation obtained in the step (5).
CN202111344183.5A 2021-11-15 2021-11-15 Laser radar ratio region transmission method based on novel foundation laser radar network Active CN113777627B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111344183.5A CN113777627B (en) 2021-11-15 2021-11-15 Laser radar ratio region transmission method based on novel foundation laser radar network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111344183.5A CN113777627B (en) 2021-11-15 2021-11-15 Laser radar ratio region transmission method based on novel foundation laser radar network

Publications (2)

Publication Number Publication Date
CN113777627A true CN113777627A (en) 2021-12-10
CN113777627B CN113777627B (en) 2022-01-18

Family

ID=78957099

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111344183.5A Active CN113777627B (en) 2021-11-15 2021-11-15 Laser radar ratio region transmission method based on novel foundation laser radar network

Country Status (1)

Country Link
CN (1) CN113777627B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114252871A (en) * 2021-12-13 2022-03-29 中国西安卫星测控中心 Radar measurement accuracy compensation method based on machine learning
CN115201791A (en) * 2022-06-08 2022-10-18 中山大学·深圳 Laser radar attenuation estimation method and system based on empirical model

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106019312A (en) * 2016-08-04 2016-10-12 浙江大学 Multi-longitudinal-mode high-spectral-resolution lidar based on interference spectrum frequency detector
CN110161532A (en) * 2019-05-30 2019-08-23 浙江大学 A method of based on multi-wavelength laser radar inverting microfluidic aerosol physical characteristic
CN113189563A (en) * 2021-04-07 2021-07-30 中国科学院大气物理研究所 Method for inverting aerosol components by utilizing laser radar ratio and depolarization ratio

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106019312A (en) * 2016-08-04 2016-10-12 浙江大学 Multi-longitudinal-mode high-spectral-resolution lidar based on interference spectrum frequency detector
CN110161532A (en) * 2019-05-30 2019-08-23 浙江大学 A method of based on multi-wavelength laser radar inverting microfluidic aerosol physical characteristic
CN113189563A (en) * 2021-04-07 2021-07-30 中国科学院大气物理研究所 Method for inverting aerosol components by utilizing laser radar ratio and depolarization ratio

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
CHRISTOPHER MCPHERSON等: "Progress in the validation of dual-wavelength aerosol retrieval models via airborne high spectral resolution lidar data", 《 2010 IEEE INTERNATIONAL GEOSCIENCE AND REMOTE SENSING SYMPOSIUM》 *
FANGWEN BAO等: "High-Spatial-Resolution Aerosol Optical Properties", 《 IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING 》 *
吴仪等: "秋冬季节珠三角典型区域气溶胶浓度垂直分布的激光雷达反演", 《遥感学报》 *
宋跃辉等: "气溶胶激光雷达比的迭代反演", 《中国激光》 *
范学花等: "中国大气气溶胶辐射特性参数的观测与研究进展", 《大气科学》 *
高飞等: "双扫描激光雷达精细探测低层大气气溶胶方法", 《光子学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114252871A (en) * 2021-12-13 2022-03-29 中国西安卫星测控中心 Radar measurement accuracy compensation method based on machine learning
CN115201791A (en) * 2022-06-08 2022-10-18 中山大学·深圳 Laser radar attenuation estimation method and system based on empirical model

Also Published As

Publication number Publication date
CN113777627B (en) 2022-01-18

Similar Documents

Publication Publication Date Title
CN113777627B (en) Laser radar ratio region transmission method based on novel foundation laser radar network
Turnbull et al. Toward quantification and source sector identification of fossil fuel CO2 emissions from an urban area: Results from the INFLUX experiment
CN100362318C (en) Atmosphere correction method of airosol optical thickness of aeronautical high-spectrum remote-sensing inversion boundary layer
Ortega et al. The state of solar energy resource assessment in Chile
CN111965666A (en) Aerosol three-dimensional distribution mapping method
CN102539336A (en) Method and system for estimating inhalable particles based on HJ-1 satellite
Liu et al. Changes in the relationship between solar radiation and sunshine duration in large cities of China
CN114910430B (en) Method for detecting traffic pollution source by trace gas horizontal distribution based on hyperspectral remote sensing
Jiang et al. Hourly 5-km surface total and diffuse solar radiation in China, 2007–2018
Ou et al. Evaluation of MERRA-2 aerosol optical and component properties over China using SONET and PARASOL/GRASP data
CN115015486B (en) Carbon emission measurement and calculation method based on regression tree model
CN103234942B (en) Utilize method and the device of sky scattering photo measure atmospheric horizontal visibility
Assi et al. Prediction of global solar radiation in Abu Dhabi
Tang et al. Cloud scattering influence on satellite laser altimetry data and its correction
CN111898296A (en) Multi-scale simulation method and system for atmospheric diffusion and sedimentation of nuclear substances
CN111914396B (en) Sub-grid terrain three-dimensional earth surface solar radiation forced effect rapid parameterization method based on high-resolution DEM data
Garrison et al. Diffuse–direct ultraviolet ratios with a compact double monochromator
Han et al. Ground-based synchronous optical instrument for measuring atmospheric visibility and turbulence intensity: Theories, design and experiments
Jia et al. Verification of MCD19A2 data and study of aerosol characteristics in Beijing-Tianjin-Hebei region
Peng et al. Visibility measurements using two-angle forward scattering by liquid droplets
Ding et al. A practical approach to improve the MODIS MCD43A products in snow-covered areas
Li et al. The Ultraviolet to Mid-infrared Extinction Law of the Taurus Molecular Cloud Based on the Gaia DR3, GALEX, APASS, Pan-STARRS1, 2MASS, and WISE Surveys
CN105403201B (en) A kind of remote sensing images atmospheric path radiation acquisition methods based on pixel analysis
Nofuentes et al. Experimental evaluation of a spectral index to characterize temporal variations in the direct normal irradiance spectrum
Lu et al. A novel hyperspectral remote sensing technique with hour-hectometer level horizontal distribution of trace gases: To accurately identify emission sources

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant