CN107561555B - Method, apparatus, computer equipment and the storage medium of inversion boundary layer height - Google Patents
Method, apparatus, computer equipment and the storage medium of inversion boundary layer height Download PDFInfo
- Publication number
- CN107561555B CN107561555B CN201710774967.9A CN201710774967A CN107561555B CN 107561555 B CN107561555 B CN 107561555B CN 201710774967 A CN201710774967 A CN 201710774967A CN 107561555 B CN107561555 B CN 107561555B
- Authority
- CN
- China
- Prior art keywords
- boundary
- layer
- height
- cloud
- stratus
- 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.)
- Active
Links
Classifications
-
- 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
Abstract
The present invention provides the method, apparatus of inversion boundary layer height, computer equipment and storage mediums, it is related to height inversion technique field, to in the case where not by other methods or data, only rely on laser radar NRB data, it is accurate to determine cloud layer and effectively reject influence of the cloud to inversion boundary layer height, so that accurate inverting has the Boundary Layer Height in the case of cloud.Technical solution of the present invention are as follows: obtain the laser radar data comprising standardized backscatter signal NRB information;Cloud layer is found according to the NRB information of laser radar data;Cloud layer is determined, determines that cloud layer is boundary stratus or non-boundary stratus;And determine the upper limit threshold height of inversion boundary layer height respectively for boundary stratus and non-boundary stratus;For cloud layer, in different threshold levels hereinafter, being based on laser radar Data Inversion convective boundary layer height.Invention also provides device, computer equipment and the storage mediums for implementing above-mentioned technical proposal.
Description
Technical field
The present invention relates to height inversion technique fields, the in particular to method based on height threshold inversion boundary layer height.
Background technique
Planetary boundary layer (PBL) is the convection current bottom directly by the atmosphere of earth surface effects, be between the earth and atmosphere into
The bridge of row substance and energy exchange, PBL determine the thickness of steam, heat and the mixing of pollutant vertical turbulence.Therefore, side
For interlayer height in Air Quality Analysis, weather forecast and climatic prediction etc. are a very important variables.It is especially arid
The sensitizing range of area and semiarid region as climate change, boundary layer processes play irreplaceable role for climate change etc..
It is many research shows that in deep convection boundary layer (CBL) convection current conveying and diffusion with local land face interaction, synoptic process, gas
It waits variation and the selection of the local strategy of sustainable development is closely related.Thus, accurate inversion boundary layer height is especially arid
The convective boundary layer height on daytime (CBLH) of semiarid region is most important.However, at present Boundary Layer Height can not direct detection obtain
, can only by using the offers such as radio sounding, wind profile radar, acoustic radar and laser radar such as temperature profile,
Aerosol profile etc. carries out inverting.
Micro-pulse lidar (MPL) by its transmitting to the wave beam of eyes fanout free region the advantages of and be widely deployed in
In remote sensing.Under normal circumstances, the emission source of aerosol is limited in surface layer, and by the top cover inversion layer that boundary layer is pushed up,
Aerosol concentration in boundary layer is significantly larger than the free atmosphere on upper layer.The vertical distribution of aerosol is considered as inversion boundary layer height
One good tracer of degree, and after laser radar can provide troposphere inferior atmospheric layer aerosol with high spatial and temporal resolution
To scattered signal profile (being NRB profile after standardization).In recent years, it is received using laser radar Data Inversion Boundary Layer Height
The extensive concern of many scholars.Specific algorithm includes gradient method (GD), curve-fitting method (CF) and wavelet covariance converter technique
Including the methods of Haar wavelet transform method of changing (H-WCT) and Mexican hat wavelet transform method (MH-WCT).However, the above method is based on
Laser radar data may not be all suitable for using inversion boundary layer height under all weather conditions.Particularly, the presence of cloud can be tight
The inverting for interfering Boundary Layer Height again, leads to have the erroneous judgement of Boundary Layer Height in the case of cloud, because of back scattering caused by cloud layer
Successively decreasing for signal is very much like with successively decreasing for backscatter signal caused by aerosol concentration decaying at the top of boundary layer.When cloud covers
When at the top of boundary layer, Boundary Layer Height is consistent with the coboundary of cloud, can obtain boundary layer by detecting the coboundary of cloud
Thickness.However, Boundary Layer Height and the cloud level are entirely different if cloud layer is located at the free atmosphere in boundary layer or more, it should
The inverting of Boundary Layer Height is carried out below cloud base.Thus, for there is the anti-of convective boundary layer height on daytime (CBLH) in the case of cloud
It drills, it is necessary to consider the influence of cloud.
Aiming at the problem that sky cloud effect utilizes laser Data Inversion Boundary Layer Height, have some researchs at present.Wherein, have
Research inversion boundary layer height under a certain specific threshold level.Cheng et al. is based in 2.5km following with gradient method
MPL Data Inversion Boundary Layer Height, however this specific threshold value different border structure corresponding for other underlying surfaces
It is not necessarily suitable for.Under normal circumstances, CBLH on daytime is probably in 1~2km.But in the Northwest such as Dunhuang of China, research shows that
This area CBLH on daytime can achieve 4km.Therefore, for different boundary layer structures, specific threshold level is not easy to determine.
Haeffelin et al. proposes STRAT-2D method to estimate the Boundary Layer Height in the case of cloud.This method, which passes through to combine, hangs down
Histogram determines mixing height to the NRB gradient with time-domain.For each time interval, STRAT-2D method provides 3
Strong gradient height caused by a different colloidal sol floor: most strong gradient height, secondary strong gradient height, the minimum strong gradient of height are high
Degree, mixing height is one of 3 gradient height.Pal et al. combines variance analysis, the specific ladder determined from STRAT-2D method
Average CBLH is obtained in degree height and the height of cloud base.Cimini et al. utilizes STRAT-2D method, provides in conjunction with microwave radiometer
Bright temperature Data Inversion obtain mixing height.Although STRAT-2D is widely used in the boundary layer for determining and having in the case of cloud
Highly, it is not at present still very that Boundary Layer Height corresponds to which of 3 gradient height that this method determines height eventually
It is bright and clear.Furthermore, it is necessary to finally determine CBLH in conjunction with some additional householder methods or other data.Li et al. people is by using one
A objective height threshold as limit level carry out inverting and have cloud in the case of Boundary Layer Height.Different from grinding for Cheng et al.
Study carefully, the isentropic condensation clevel CCL that Li et al. people research uses the Temperature Datum provided by microwave radiometer to calculate is as inverting CBLH
Threshold level, i.e., only CCL NRB data below is used to inversion boundary layer height.This method weakens to a certain extent
Positioned at free atmosphere cloud layer for Boundary Layer Height inverting interference.However, this method is maximum the disadvantage is that threshold level
The calculating accuracy that depends not only upon other data, and calculated dependent on CCL.
So how in the case where not by other methods or data, laser radar NRB data is only relied on, it is accurate to determine
Cloud layer simultaneously effectively rejects influence of the cloud to inverting CBLH, so that accurate inverting has the CBLH in the case of cloud to be contemplated that at present
Problem.
Summary of the invention
The present invention provides the method, apparatus of inversion boundary layer height, computer equipment and storage mediums, not borrow
In the case where helping other methods or data, laser radar NRB data is only relied on, it is accurate to determine cloud layer and effectively reject cloud to inverting
The influence of Boundary Layer Height, so that accurate inverting has the Boundary Layer Height in the case of cloud.
A kind of method of inversion boundary layer height provided by the present invention, detailed process include:
Obtain the laser radar data comprising standardized backscatter signal NRB information.
Cloud layer is found according to the NRB information of laser radar data.
Cloud layer is determined, determines that cloud layer is boundary stratus or non-boundary stratus;And it is directed to boundary stratus and non-boundary
Stratus determines the upper limit threshold height of inversion boundary layer height respectively.
For cloud layer, in different threshold levels hereinafter, being based on laser radar Data Inversion convective boundary layer height.
In the embodiment of the present invention, cloud layer is found according to laser radar data NRB information, further includes: according to NRB information, meter
It calculates NRB gradient and NRB relative gradient: finding the strong positive negative gradient of NRB in laser radar data to layer, the strong positive and negative ladder of NRB
The gradient layer for spending the relative gradient threshold value for being greater than or equal to setting to the NRB relative gradient at positive gradient in layer is determined as cloud layer.
In the embodiment of the present invention, the relative gradient threshold value that sets is 55%.
In the embodiment of the present invention, the strong positive negative gradient of NRB is found in laser radar data to layer, further includes: in laser
In Radar Data, the most strong negative gradient of NRB is found from top to bottom, is finding maximum just in the first setting range below most strong negative gradient
Gradient constitutes the strong positive negative gradient of NRB to layer.
In the embodiment of the present invention, cloud layer is determined, determines that cloud layer is boundary stratus or non-boundary stratus, further includes:
In the cloud layer that the NRB profile provided based on laser radar data is found, the cloud layer that height above sea level is minimum is taken, determines the minimum cloud of height above sea level
For boundary stratus or non-boundary stratus;Most strong negative gradient is found under the minimum cloud layer of height above sea level;Under the minimum cloud layer of height above sea level
Most strong negative gradient height be z1;The height of cloud base of the minimum cloud layer of height above sea level is z2;The coboundary of the cloud of the minimum cloud layer of height above sea level
Height is z3;The bottom of NRB profile is to z1Between NRB average value be a1;z1And z2Between signal averaging be a2: setting threshold
Value is β;
If β × a1<a2<a1, then the minimum cloud layer of height above sea level is boundary stratus;
If a2≥a1Or a2≤β×a1, then the minimum cloud layer of height above sea level is non-boundary stratus.
In the embodiment of the present invention, the upper limit threshold of inversion boundary layer height is determined respectively for boundary stratus and non-boundary stratus
Value height, further includes: for boundary stratus, threshold level is determined that the positive gradient of the above NRB in boundary above the clouds goes out for the first time
Existing position.For non-boundary stratus, threshold level is defined on NRB in the second setting range of cloud base or less and starts the position increased
It sets.
In the embodiment of the present invention, for cloud layer, in different threshold levels hereinafter, being based on the convection current of laser radar Data Inversion
Boundary Layer Height, further includes: cloud layer is directed to, in different threshold levels hereinafter, utilizing gradient method, Ha Er haar wavelet transformation
Method, sombrero Mexican Hat Wavelet Transform and curve-fitting method distinguish inverting convective boundary layer height.
Another embodiment of the present invention additionally provides a kind of device of inversion boundary layer height, which includes laser radar money
Material obtains module, cloud layer searching module, cloud layer determination module and inverting module:
Laser radar data obtains module, includes swashing for standardized backscatter signal NRB profile information for obtaining
Optical radar data, and laser radar data is sent to cloud layer searching module.
Cloud layer searching module is used to find cloud layer according to NRB information, and cloud layer lookup result is sent into cloud layer determination module.
Cloud layer determination module determines that cloud layer is boundary stratus or non-boundary stratus for determining cloud layer;And it is directed to
Boundary stratus and non-boundary stratus determine the upper limit threshold height of inversion boundary layer height respectively, and cloud layer is determined that result is sent into
Inverting module.
Inverting module is used to be directed to cloud layer, laser radar data is based on, in different threshold levels hereinafter, inverting convection current side
Interlayer height.
Another embodiment of the present invention additionally provides a kind of computer equipment, including memory, processor and is stored in
On reservoir and the computer program that can run on a processor, processor realize that any of the above embodiment such as provides when executing program
Method.
Another embodiment of the present invention additionally provides a kind of computer readable storage medium, is stored thereon with computer program,
The method provided such as any of the above embodiment is provided when computer program is executed by processor.
Other features and advantages of the present invention will be illustrated in the following description, also, part becomes from specification
It is clear that understand through the implementation of the invention.The objectives and other advantages of the invention can by written specification,
Specifically noted structure is achieved and obtained in claims and attached drawing.
Below by drawings and examples, technical scheme of the present invention will be described in further detail.
Detailed description of the invention
Attached drawing is used to provide further understanding of the present invention, and constitutes part of specification, with reality of the invention
It applies example to be used to explain the present invention together, not be construed as limiting the invention.In the accompanying drawings:
Fig. 1 is the method flow schematic diagram of inversion boundary layer height provided by one embodiment of the present invention;
Fig. 2 is the apparatus structure schematic diagram for the inversion boundary layer height that another embodiment of the present invention provides;
Fig. 3 is Potential temperature Profile and the station SACOL of the Yuzhong sounding station of the 20:00BJT of September in 2010 3 days (Beijing time)
Standardized NRB profile and the convective boundary layer height of institute's inverting.Fig. 3 (a) is Potential temperature Profile and with 1.5 degree of position Wen Zengliang
The CBLH of method inverting;Fig. 3 (b) is NRB profile and by GD method, H-WCT method, and MH-WCT method and CF method are not added any
The CBLH of inverting under limited case;Fig. 3 (c) is NRB profile and by GD method, H-WCT method, MH-WCT method and CF method
The CBLH of inverting under threshold level limitation.
Fig. 4 (a), (b), (c), (d) have cloud example 20:00BJT by being utilized respectively GD based on laser radar data for 21,
The CBLH (abscissa) of H-WCT, the MH-WCT and CF inverting and CBLH that 1.5 degree of position temperature method of addition invertings are utilized based on Potential temperature Profile
The Comparative result of (ordinate).
Specific embodiment
Hereinafter, preferred embodiments of the present invention will be described with reference to the accompanying drawings, it should be understood that preferred reality described herein
Apply example only for the purpose of illustrating and explaining the present invention and is not intended to limit the present invention.
Embodiment 1
The embodiment of the invention provides a kind of method of inversion boundary layer height, mainly for the boundary layer having in the case of cloud
The inverting of height, this method flow chart is as shown in Figure 1, specifically comprise the following steps:
S101, the laser radar data comprising standardized backscatter signal NRB information is obtained.
S102, cloud layer is found according to the NRB information of laser radar data.
The step specifically includes as follows: according to NRB information, calculating NRB gradient and NRB relative gradient;It is provided in laser radar
The strong positive negative gradient of NRB is found in material to layer, i.e., in laser radar data, the most strong negative gradient of NRB is found from top to bottom, most
Maximum positive gradient is found in the following setting range of strong negative gradient, constitutes the strong positive negative gradient of NRB to layer.
NRB relative gradient in the strong positive and negative gradient layer of NRB at positive gradient is greater than or equal to the relative gradient threshold value of setting
Layer be determined as cloud layer.
In laser radar data, the most strong negative gradient of NRB is found from top to bottom, the first setting model below most strong negative gradient
The maximum positive gradient of interior searching is enclosed, constitutes positive and negative NRB gradient to layer.In the embodiment of the present invention, the first setting range is 1 km km.
In the embodiment of the present invention, the relative gradient threshold value that sets is 55%.I.e. for strong pairs of positive negative gradient to layer (just
Gradient is under, and negative gradient is upper), the layer knot that the relative gradient at positive gradient reaches 55% is judged as cloud layer in the present invention.
S103, cloud layer is determined, determines that cloud layer is boundary stratus or non-boundary stratus;And for boundary stratus and non-
Boundary stratus determines the upper limit threshold height of inversion boundary layer height respectively.
In the cloud layer that the NRB profile provided based on laser radar data is found, the cloud layer that height above sea level is minimum is taken, determines height above sea level
Minimum cloud is boundary stratus or non-boundary stratus.
Most strong negative gradient is found under the minimum cloud layer of height above sea level.
The height of most strong negative gradient under the minimum cloud layer of height above sea level is z1.The height of cloud base of the minimum cloud layer of height above sea level is z2。
The coboundary height of the cloud of the minimum cloud layer of height above sea level is z3.The bottom of NRB profile is to z1Between NRB average value be a1。z1And z2
Between signal averaging be a2.Given threshold is β.In the embodiment of the present invention, given threshold 75%.
If β × a1<a2<a1, then the minimum cloud layer of height above sea level is boundary stratus;
If a2≥a1Or a2≤β×a1, then the minimum cloud layer of height above sea level is non-boundary stratus.
For boundary stratus, threshold level is determined into the position that the positive gradient of the above NRB in boundary above the clouds occurs for the first time
It sets.
For non-boundary stratus, threshold level is defined on NRB in the second setting range of cloud base or less and starts the position increased
It sets.In the embodiment of the present invention, the second setting range is set small as far as possible, i.e., so that threshold level is defined on the NRB close to cloud base
Start the position increased.
S104, it is directed to cloud layer, it is in different threshold levels hereinafter, high based on laser radar Data Inversion convective boundary layer
Degree.
In the present embodiment, for cloud layer, in different threshold levels hereinafter, utilizing gradient method, Ha Er haar wavelet transformation
Method, sombrero Mexican Hat Wavelet Transform and curve-fitting method distinguish inverting convective boundary layer height.
Embodiment 2
Another embodiment of the presently claimed invention provides a kind of device of inversion boundary layer height, any of the above-described for implementing
Embodiment provide inversion boundary layer height method, the device include laser radar data obtain module, cloud layer searching module,
Cloud layer determination module and inverting module, specific structure are as shown in Figure 2.
Laser radar data obtains module 201, includes standardized backscatter signal NRB profile information for obtaining
Laser radar data, and laser radar data is sent to cloud layer searching module 202;
Cloud layer lookup result for finding cloud layer according to NRB information, and is sent into cloud layer and determined by cloud layer searching module 202
Module 203;
Cloud layer determination module 203 determines that cloud layer is boundary stratus or non-boundary stratus for determining cloud layer;And
It determines the upper limit threshold height of inversion boundary layer height respectively for boundary stratus and non-boundary stratus, and cloud layer is determined into result
It is sent into inverting module 204;
Inverting module 204 is based on laser radar data, in different threshold levels hereinafter, inverting pair for being directed to cloud layer
Laminar boundary layer height.
Embodiment 3
A kind of computer equipment, including memory, processor and storage can be run on a memory and on a processor
Computer program, processor execute program when perform the steps of
Obtain the laser radar data comprising standardized backscatter signal NRB profile information;
Cloud layer is found according to NRB information;
Cloud layer is determined, determines that cloud layer is boundary stratus or non-boundary stratus;And it is directed to boundary stratus and non-boundary
Stratus determines the upper limit threshold height of inversion boundary layer height respectively;
For cloud layer, in different threshold levels hereinafter, being based on laser radar Data Inversion convective boundary layer height.
Embodiment 4
A kind of computer readable storage medium is stored thereon with computer program, when computer program is executed by processor
It performs the steps of
Obtain the laser radar data comprising standardized backscatter signal NRB profile information;
Cloud layer is found according to NRB information;
Cloud layer is determined, determines that cloud layer is boundary stratus or non-boundary stratus;And it is directed to boundary stratus and non-boundary
Stratus determines the upper limit threshold height of inversion boundary layer height respectively;
For cloud layer, in different threshold levels hereinafter, being based on laser radar Data Inversion convective boundary layer height.
Embodiment 5
The embodiment of the present invention proposes the new more simple and effective height threshold method not influenced by underlying surface of one kind to reject
Interference of the cloud to Boundary Layer Height inverting, and based on Lanzhou University's semiarid climate and environmental monitoring station (SACOL;35.95 ° N,
104.14°E;Height above sea level 1965.8m) laser radar data utilize GD, CF, H-WCT and MH-WCT inverting has in the case of cloud
CBLH.The judgement of cloud layer, the determination of threshold level and the determination of final CBLH only rely upon laser radar back scattering letter
Number.Meanwhile on the basis of determining cloud layer, the continuity by analyzing Beneath Clouds NRB, which distinguishes cloud, to be boundary stratus or is located at
Cloud, that is, non-boundary stratus more than boundary layer.It is proposed the threshold value of inversion boundary layer height respectively for different cloud layer position situations
Highly.In addition, using by Yuzhong station (35.87 ° of N, 104.15 ° of E;Height above sea level 1875m;It is with the SACOL horizontal distance stood
The low 90.8m of 8.3km, height above sea level ratio SACOL) the CBLH diagnosis that is calculated based on 1.5 degree of position temperature method of additions of sounding Potential temperature Profile
The CBLH of above-mentioned inverting assesses the reasonability of this method.
Fig. 3 is the Potential temperature Profile of the Yuzhong sounding station of the 20:00BJT on the 3rd of September in 2010 and the standardized NRB at the station SACOL
Profile and the convective boundary layer height of institute's inverting.Scheming (a) is Potential temperature Profile and the CBLH with 1.5 degree of position temperature method of addition invertings
(solid black lines).As can be seen that being 1.46km by the CBLH of 1.5 degree of position temperature method of addition invertings.NRB profile and by GD method,
The CBLH of H-WCT method, MH-WCT method and CF inversion method such as figure (b) (any restrictions are not added in inverting) and figure (c) are (in threshold value
Inverting under height limitation) shown in.It is learnt by NRB profile combination NRB relative gradient, in 2.8km to there is a stratus between 3.2km
In the presence of.Meanwhile being respectively positioned near the coboundary of cloud by the Boundary Layer Height of various inversion methods, average out to 3.15km is much high
In the Boundary Layer Height by 1.5 degree of position temperature method of addition invertings.After giving maximum height limit, by the boundary layer of various inversion methods
Height average out to 1.55km, it is obviously more reasonable compared to without restriction.
Fig. 4, which is 21, has cloud example 20:00BJT by being utilized respectively GD based on laser radar data, H-WCT, MH-WCT with
The CBLH (abscissa) of CF inverting and based on Potential temperature Profile using 1.5 degree of position temperature method of addition invertings CBLH (ordinate) result
Comparison.The CBLH result that any restrictions inverting is not added in triangle point representative compares, and meter Zi Dian represents the inverting given after limiting
CBLH Comparative result.As can be seen that, except a small number of exceptions, most of triangle points are far from 1:1 line zero when any restrictions are not added
The lower right for being distributed in 1:1 line is dissipated, that is, has the Boundary Layer Height of the distinct methods inverting based on laser radar data in the presence of cloud
Most of Boundary Layer Heights being significantly larger than by 1.5 degree of position temperature method of addition invertings.In contrast, anti-under threshold level limitation
Boundary Layer Height (meter Zi Dian) integrated distribution drilled is based on two kinds of prospecting tools, different detection method invertings near 1:1 line
Boundary Layer Height all relatively.Thus illustrate, when by laser radar Data Inversion CBLH, cloud can generate it serious dry
It disturbs, and the threshold level that the invention proposes can effectively propose the interference of cloud, when there is cloud, inverting obtains more accurately closing
The Boundary Layer Height of reason.
Obviously, various changes and modifications can be made to the invention without departing from essence of the invention by those skilled in the art
Mind and range.In this way, if these modifications and changes of the present invention belongs to the range of the claims in the present invention and its equivalent technologies
Within, then the present invention is also intended to include these modifications and variations.
Claims (8)
1. a kind of method of inversion boundary layer height, which is characterized in that this method comprises:
Obtain the laser radar data comprising standardized backscatter signal NRB information;
Cloud layer is found according to the NRB information of the laser radar data;
The cloud layer is determined, determines that cloud layer is boundary stratus or non-boundary stratus;And it is directed to boundary stratus and non-boundary
Stratus determines the upper limit threshold height of inversion boundary layer height respectively;
For the cloud layer, in different threshold levels hereinafter, being based on the laser radar Data Inversion convective boundary layer height;
It is described that the cloud layer is determined, determine that cloud layer is boundary stratus or non-boundary stratus, comprising:
In the cloud layer that the NRB profile provided based on the laser radar data is found, the cloud layer that height above sea level is minimum is taken, determines height above sea level
Minimum cloud is boundary stratus or non-boundary stratus;
Most strong negative gradient is found under the minimum cloud layer of the height above sea level;
The height of most strong negative gradient under the minimum cloud layer of the height above sea level is z1;
The height of cloud base of the minimum cloud layer of the height above sea level is z2;
The coboundary height of the cloud of the minimum cloud layer of the height above sea level is z3;
The bottom of the NRB profile is to z1Between NRB average value be a1;
z1And z2Between signal averaging be a2:
Given threshold is β;
If β × a1<a2<a1, then the minimum cloud layer of the height above sea level is boundary stratus;
If a2≥a1Or a2≤β×a1, then the minimum cloud layer of the height above sea level is non-boundary stratus;
The upper limit threshold height for determining inversion boundary layer height respectively for boundary stratus and non-boundary stratus, comprising:
For boundary stratus, threshold level is determined into the position that the positive gradient of the above NRB in boundary above the clouds occurs for the first time;
For non-boundary stratus, threshold level is defined on NRB in the second setting range of cloud base or less and starts the position increased.
2. the method as described in claim 1, which is characterized in that described to find cloud according to the laser radar data NRB information
Layer, further includes:
According to the NRB information, NRB gradient and NRB relative gradient are calculated:
The strong positive negative gradient of NRB is found in the laser radar data to layer, the strong positive negative gradient of the NRB is to ladder positive in layer
The gradient layer for the relative gradient threshold value that NRB relative gradient at degree is greater than or equal to setting is determined as cloud layer.
3. method according to claim 2, which is characterized in that the relative gradient threshold value set is 55%.
4. method according to claim 2, which is characterized in that described to find the strong just of NRB in the laser radar data
Negative gradient is to layer, further includes:
In the laser radar data, the most strong negative gradient of NRB is found from top to bottom, and first sets below the most strong negative gradient
Determine to find maximum positive gradient in range, constitutes the strong positive negative gradient of the NRB to layer.
5. the method as described in claim 1, which is characterized in that it is described be directed to the cloud layer, in different threshold levels hereinafter,
Based on the laser radar Data Inversion convective boundary layer height, further includes:
For the cloud layer, in different threshold levels hereinafter, utilizing gradient method, Ha Er haar Wavelet Transform, sombrero
Mexican Hat Wavelet Transform and curve-fitting method distinguish inverting convective boundary layer height.
6. a kind of device of inversion boundary layer height, which is characterized in that the device includes that laser radar data obtains module, cloud layer
Searching module, cloud layer determination module and inverting module:
The laser radar data obtains module, includes swashing for standardized backscatter signal NRB profile information for obtaining
Optical radar data, and the laser radar data is sent to the cloud layer searching module;
The cloud layer searching module is used to find cloud layer according to the NRB information, and cloud layer lookup result is sent into the cloud layer
Determination module;
The cloud layer determination module determines that cloud layer is boundary stratus or non-boundary stratus for determining the cloud layer;And
It determines the upper limit threshold height of inversion boundary layer height respectively for boundary stratus and non-boundary stratus, and cloud layer is determined into result
It is sent into the inverting module;
The inverting module is used to be directed to cloud layer, the laser radar data is based on, in different threshold levels hereinafter, inverting pair
Laminar boundary layer height;
It is described that the cloud layer is determined, determine that cloud layer is boundary stratus or non-boundary stratus, comprising:
In the cloud layer that the NRB profile provided based on the laser radar data is found, the cloud layer that height above sea level is minimum is taken, determines height above sea level
Minimum cloud is boundary stratus or non-boundary stratus;
Most strong negative gradient is found under the minimum cloud layer of the height above sea level;
The height of most strong negative gradient under the minimum cloud layer of the height above sea level is z1;
The height of cloud base of the minimum cloud layer of the height above sea level is z2;
The coboundary height of the cloud of the minimum cloud layer of the height above sea level is z3;
The bottom of the NRB profile is to z1Between NRB average value be a1;
z1And z2Between signal averaging be a2:
Given threshold is β;
If β × a1<a2<a1, then the minimum cloud layer of the height above sea level is boundary stratus;
If a2≥a1Or a2≤β×a1, then the minimum cloud layer of the height above sea level is non-boundary stratus;
The upper limit threshold height for determining inversion boundary layer height respectively for boundary stratus and non-boundary stratus, comprising:
For boundary stratus, threshold level is determined into the position that the positive gradient of the above NRB in boundary above the clouds occurs for the first time;
For non-boundary stratus, threshold level is defined on NRB in the second setting range of cloud base or less and starts the position increased.
7. a kind of computer equipment, can run on a memory and on a processor including memory, processor and storage
Computer program, which is characterized in that the processor realizes such as any the method for claim 1-5 when executing described program.
8. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the computer program quilt
Such as claim 1-5 any the method is realized when processor executes.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710774967.9A CN107561555B (en) | 2017-08-31 | 2017-08-31 | Method, apparatus, computer equipment and the storage medium of inversion boundary layer height |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710774967.9A CN107561555B (en) | 2017-08-31 | 2017-08-31 | Method, apparatus, computer equipment and the storage medium of inversion boundary layer height |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107561555A CN107561555A (en) | 2018-01-09 |
CN107561555B true CN107561555B (en) | 2019-07-09 |
Family
ID=60978527
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710774967.9A Active CN107561555B (en) | 2017-08-31 | 2017-08-31 | Method, apparatus, computer equipment and the storage medium of inversion boundary layer height |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107561555B (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109710683B (en) * | 2018-11-23 | 2020-11-03 | 中国气象科学研究院 | Cloud area discrimination result imaging method |
CN109710657B (en) * | 2018-11-23 | 2020-11-03 | 中国气象科学研究院 | Cloud area judgment method |
CN109541639A (en) * | 2018-12-25 | 2019-03-29 | 天津珞雍空间信息研究院有限公司 | A kind of inversion boundary layer height method based on particle cluster |
CN110031868A (en) * | 2019-05-15 | 2019-07-19 | 国耀量子雷达科技有限公司 | A method of based on coherent wind laser radar carrier-to-noise ratio inversion boundary layer height |
CN112348388A (en) * | 2020-11-16 | 2021-02-09 | 中科三清科技有限公司 | Ozone leading precursor identification method and device based on emission data |
CN112526547B (en) * | 2020-11-18 | 2024-02-23 | 南京泰爱信科技有限公司 | Atmospheric boundary layer classification method and device based on wind-measuring laser radar |
CN112526638B (en) * | 2020-11-27 | 2022-07-29 | 中国气象局气象探测中心 | Cloud boundary identification method based on multi-source observation data and related equipment |
CN113189616B (en) * | 2021-05-10 | 2023-05-02 | 北京城市气象研究院 | Method for determining boundary layer height of atmosphere and related equipment |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103135113A (en) * | 2013-02-05 | 2013-06-05 | 中国科学院大气物理研究所 | Method for measuring height of atmospheric boundary layer |
CN106569226A (en) * | 2016-10-28 | 2017-04-19 | 兰州大学 | Method using laser radar data to invert boundary layer height under cloudy condition |
-
2017
- 2017-08-31 CN CN201710774967.9A patent/CN107561555B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103135113A (en) * | 2013-02-05 | 2013-06-05 | 中国科学院大气物理研究所 | Method for measuring height of atmospheric boundary layer |
CN106569226A (en) * | 2016-10-28 | 2017-04-19 | 兰州大学 | Method using laser radar data to invert boundary layer height under cloudy condition |
Non-Patent Citations (3)
Title |
---|
Atmospheric Boundary Layer Height Monitoring Using a Kalman Filter and Backscatter Lidar Returns;Diego Lange et al.;《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》;20140831;第52卷(第8期);第4717-4728页 |
SYNERGETIC USE OF LIDAR AND MICROWAVE RADIOMETER OBSERVATIONS FOR BOUNDARY-LAYER HEIGHT DETECTION;Umar Saeed et al.;《2015 IEEE 》;20151231;第3945-3948页 |
基于激光雷达资料的小波变换法反演边界层高度的方法;李红等;《干旱气象》;20150228;第33卷(第1期);第78-88页 |
Also Published As
Publication number | Publication date |
---|---|
CN107561555A (en) | 2018-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107561555B (en) | Method, apparatus, computer equipment and the storage medium of inversion boundary layer height | |
Zhang et al. | On the computation of planetary boundary-layer height using the bulk Richardson number method | |
Drew et al. | Observations of wind speed profiles over Greater London, UK, using a Doppler lidar | |
LeMone et al. | Objectively determined fair-weather CBL depths in the ARW-WRF model and their comparison to CASES-97 observations | |
Banks et al. | Performance evaluation of the boundary-layer height from lidar and the Weather Research and Forecasting model at an urban coastal site in the north-east Iberian Peninsula | |
Rudnick et al. | Cyclonic eddies in the Gulf of Mexico: Observations by underwater gliders and simulations by numerical model | |
Tang et al. | Lagrangian coherent structure analysis of terminal winds detected by lidar. Part I: Turbulence structures | |
Smyth et al. | Three dimensional airflow patterns within a coastal trough–bowl blowout during fresh breeze to hurricane force winds | |
CN102628944B (en) | Stratus cloud and convective cloud automatic recognition method based on Doppler radar data | |
CN112711899B (en) | Fusion prediction method for height of evaporation waveguide | |
Haza et al. | Drogue-loss detection for surface drifters during the Lagrangian Submesoscale Experiment (LASER) | |
Fan et al. | Observations of Irminger Sea anticyclonic eddies | |
CN105891833A (en) | Method of identifying warm cloud precipitation rate based on Doppler radar information | |
Wang et al. | Robust solution for boundary layer height detections with coherent doppler wind lidar | |
CN110031868A (en) | A method of based on coherent wind laser radar carrier-to-noise ratio inversion boundary layer height | |
Mao et al. | Determination of the boundary layer top from lidar backscatter profiles using a Haar wavelet method over Wuhan, China | |
Fratantoni et al. | Mean structure and dynamics of the shelfbreak jet in the Middle Atlantic Bight during fall and winter | |
US20150051883A1 (en) | Process to identify and classify oil seep areas at the seabed through inverse modeling | |
Alfaro et al. | Thermodynamic constraints on the morphology of simulated midlatitude squall lines | |
Wang et al. | Measuring boundary-layer height under clear and cloudy conditions using three instruments | |
Falocchi et al. | A method to determine the characteristic time‐scales of quasi‐isotropic surface‐layer turbulence over complex terrain: A case‐study in the Adige Valley (Italian Alps) | |
Berg et al. | Lidar-observed stress vectors and veer in the atmospheric boundary layer | |
CN106680833A (en) | Method for recognizing thickness of mixed layer on basis of Logistic atmosphere extinction profile | |
CN113553766A (en) | Method for inverting arctic snow depth by using machine learning | |
CN110132246B (en) | Water depth remote sensing detection method based on residual partitioning |
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 |