CN111666691B - Statistical method for atmospheric optical turbulence parameters - Google Patents
Statistical method for atmospheric optical turbulence parameters Download PDFInfo
- Publication number
- CN111666691B CN111666691B CN202010530547.8A CN202010530547A CN111666691B CN 111666691 B CN111666691 B CN 111666691B CN 202010530547 A CN202010530547 A CN 202010530547A CN 111666691 B CN111666691 B CN 111666691B
- Authority
- CN
- China
- Prior art keywords
- pitch angle
- atmospheric optical
- atmospheric
- optical turbulence
- solar
- 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
Images
Classifications
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W2001/003—Clear air turbulence detection or forecasting, e.g. for aircrafts
-
- 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)
- Environmental & Geological Engineering (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a statistical method of atmospheric optical turbulence parameters, which comprises the following steps: the elevation angle values of the sun before and after passing the top are distinguished by adopting a sun azimuth mark, when the sun azimuth angle is 0-180 degrees, the atmospheric optical turbulence parameter is the upper half-day data, and when the sun azimuth angle is 180-360 degrees, the atmospheric optical turbulence parameter is the lower half-day data; normalizing the value of the pitch angle larger than 0 degree and the maximum value of the pitch angle, and normalizing the value of the pitch angle smaller than 0 degree and the minimum value of the pitch angle; calculating a statistical analysis rule of the solar pitch angle alpha after normalization of the atmospheric optical turbulence parameters; and (3) the normalized sun pitch angle alpha is subjected to periodicity by utilizing a trigonometric function to obtain a statistical analysis rule of the atmospheric optical turbulence parameter along with the periodic sun pitch angle alpha, and the method can be used for quantitatively analyzing the atmospheric optical turbulence parameter measured values in different months and different places.
Description
Technical Field
The invention relates to a statistical method, in particular to a statistical method of atmospheric optical turbulence parameters.
Background
The random variation of the refractive index due to the random fluctuations of the atmospheric temperature, atmospheric density is called atmospheric optical turbulence. When light waves propagate in turbulent atmosphere, light beam phase distortion and light intensity random fluctuation occur, and accordingly atmospheric turbulence effects such as light beam expansion, drifting and flickering are caused. Atmospheric coherence length r 0 And near-surface turbulence intensity (atmospheric refractive index junction)Constant of structure) C n 2 Is a fundamental parameter describing the effect of atmospheric turbulence. In the application fields of laser engineering such as astronomical site selection, laser atmospheric communication and the like, the statistical characteristics of the atmospheric turbulence parameters have very important significance.
In recent years, a plurality of documents report the measurement results of atmospheric optical parameters of different sites in China. The national astronomical table is based on the purpose of astronomical site selection, long-term observation is carried out in places such as Yunnan Kunming, lijiang, hebei Xinglong, tibet Ali and the like, and a large number of atmospheric optical parameters are obtained; the institute of optical precision machinery of Anhui of Chinese academy of sciences carries out long-term measurement in Qinghai, xinjiang and the like based on the purpose of atmospheric optical parameter research, and accumulates a large amount of data for analyzing atmospheric optical parameters in different areas and different periods.
In the visible literature, the statistical analysis method for the atmospheric optical turbulence parameter is to analyze the daily variation rule of the atmospheric turbulence parameter according to the Beijing time statistics according to the place and month (or season). Considering that the daily change of the atmospheric turbulence parameter is actually caused by the air density change caused by the sun irradiating the earth surface, the statistical analysis of the daily change rule of the atmospheric turbulence parameter along with the solar elevation angle and azimuth angle can better embody the change characteristics of the atmospheric turbulence parameter in different regions and in different months (or seasons).
In literature 1 (Chen Xiaowei, li Xuebin, etc., typical region of western china, whole-floor optical turbulence observation and analysis, optical science 2016), literature 2 (Song Zhengfang, yang Gaochao, liu Xiaochun, etc., yunnan astronomical stage atmospheric seeing measurement, quantum electronics, 1997), literature 3 (seal double-connected, zhang Zhigang, jiang Xiwen, etc., statistical analysis of gobi near-ground atmospheric refractive index structural constants, intense laser and particle beam, 2012), literature 4 (Lvyu, park moth, et al, near-ground atmospheric refractive index structural constants characteristic analysis in autumn and winter in western region, atmospheric and environmental optics, 2013), statistical analysis of atmospheric optical turbulence parameters at one or more sites over a period of time is involved. As can be seen from the description of the above documents, the data statistics method adopted by the method is the atmospheric optical turbulence parameter (atmospheric coherence length r) which changes with time 0 Near toIntensity of ground atmospheric turbulence C n 2 Etc.) are statistically averaged according to seasons, and then the daily change rule of atmospheric optical turbulence parameters of one or more places along with time is described and analyzed in a comparison mode.
However, when the daily variation of the atmospheric optical turbulence parameter is counted according to time in the process of analyzing the characteristics of the atmospheric optical turbulence parameter, the following defects exist:
1) When the measurement is carried out at the same place, the measurement data of different months (seasons) have the phenomenon of obviously changing and deviating along with the sunrise and sunset moments.
2) When the data are measured at different places, due to different longitudes, sunrise and sunset moments are different, and the measured data have obvious offset.
Due to the above disadvantages, when comparing and analyzing the atmospheric optical turbulence parameter measurement data, qualitative descriptions of sunrise time, sunset time, switching time, day, night, etc. are often used, and then quantitative analysis cannot be performed on the atmospheric optical turbulence parameter measurement values in different months and different places.
Disclosure of Invention
The invention aims to overcome the defects of the prior art and provides a statistical method for atmospheric optical turbulence parameters, which can quantitatively analyze the measured values of the atmospheric optical turbulence parameters in different months and different places.
In order to achieve the above purpose, the statistical method of the atmospheric optical turbulence parameter according to the present invention comprises the following steps:
acquiring atmospheric optical turbulence parameters of an acquisition place, wherein the atmospheric optical turbulence parameters comprise near-ground atmospheric turbulence intensity and atmospheric coherence length; acquiring longitude and latitude and altitude of a collection place, and calculating a solar pitch angle and a solar azimuth angle corresponding to the collection time of the atmospheric optical turbulence parameters according to an astronomy formula;
the elevation angle values of the sun before and after passing the top are distinguished by adopting a sun azimuth mark, when the sun azimuth angle is 0-180 degrees, the atmospheric optical turbulence parameter is upper-half-day data, and when the sun azimuth angle is 180-360 degrees, the atmospheric optical turbulence parameter is lower-half-day data;
normalizing the value of the pitch angle larger than 0 degree and the maximum value of the pitch angle, and normalizing the value of the pitch angle smaller than 0 degree and the minimum value of the pitch angle;
calculating a statistical analysis rule of the solar pitch angle alpha after normalization of the atmospheric optical turbulence parameters;
and (4) utilizing a trigonometric function to carry out periodicity on the normalized solar pitch angle alpha, obtaining a statistical analysis rule of the atmospheric optical turbulence parameter along with the solar pitch angle alpha after the periodicity, and completing the statistics of the atmospheric optical turbulence parameter.
The normalized sun pitch angle alpha is cycled by a trigonometric function, i.e.
The invention has the following beneficial effects:
in the statistical method of the atmospheric optical turbulence parameters, during specific operation, the elevation numerical values before and after the sun passes through the top are distinguished by adopting the sun azimuth mark, so that the problem of the repetition of the elevation numerical values before and after the sun passes through the top is solved; normalizing the value of the pitch angle larger than 0 degree and the maximum value of the pitch angle, and normalizing the value of the pitch angle smaller than 0 degree and the minimum value of the pitch angle so as to solve the problem that the solar pitch angle is inconvenient for daily change rule statistics; and meanwhile, the normalized solar pitch angle is cycled by utilizing a trigonometric function, so that the problem that average data of the same day are difficult to unify in the same coordinate system is solved, the statistical analysis rule of the atmospheric optical turbulence parameter along with the cycled solar pitch angle is obtained, and the quantitative analysis of the atmospheric optical turbulence parameter measured values in different months and different places is realized.
Drawings
FIG. 1 is a comparative analysis chart of atmospheric optical turbulence parameter statistics for different months;
FIG. 2 is a graph of the elevation of the sun as a function of time;
FIG. 3 is a graph of solar azimuth angle over time;
FIG. 4 is a graph of the maximum and minimum solar elevation angles for each day of the year;
FIG. 5a is a statistical analysis rule graph of atmospheric optical turbulence parameters with normalized solar pitch angle over a half-day period of 2016 years 3;
FIG. 5b is a graph of the statistical analysis of the atmospheric optical turbulence parameters with the normalized solar pitch angle in the second half of 3 months in 2016 years;
FIG. 6 is a diagram of a statistical analysis of atmospheric optical turbulence parameters with periodic solar pitch angles;
FIG. 7 is a graph of statistical comparative analysis of atmospheric optical turbulence parameters for different months.
Detailed Description
The invention is described in further detail below with reference to the accompanying drawings:
in the process of analyzing the atmospheric optical turbulence parameter characteristics, the invention provides a statistical method of the atmospheric optical turbulence parameters combining the atmospheric optical parameter characteristics with the normalized solar elevation angle, which specifically comprises the following steps:
acquiring atmospheric optical turbulence parameters of an acquisition place, wherein the atmospheric optical turbulence parameters comprise near-ground atmospheric turbulence intensity and atmospheric coherence length; acquiring longitude and latitude and altitude of a collection place, and calculating a solar pitch angle and a solar azimuth angle corresponding to the collection time of the atmospheric optical turbulence parameters according to an astronomy formula;
in order to overcome the problem of repeated elevation values before and after the sun passes through the top, as shown in points A and B in fig. 2, the sun azimuth mark is adopted for distinguishing, as shown in fig. 3, when the sun azimuth is 0-180 degrees, the atmospheric optical turbulence data is the data of the upper half day, and when the sun azimuth is 180-360 degrees, the atmospheric optical turbulence data is the data of the lower half day.
In addition, in order to overcome the defect that the maximum value and the minimum value of the solar pitch angle change every day, as shown in fig. 4, daily change rule statistics is inconvenient, the method adopts the value of the pitch angle larger than 0 degree to normalize with the maximum value of the pitch angle, and the value of the pitch angle smaller than 0 degree to normalize with the minimum value of the pitch angle.
Averaging the normalized atmospheric optical turbulence parameter data according to a set angle interval to obtain a statistical analysis rule of the atmospheric optical turbulence parameter along with the normalized solar pitch angle alpha, as shown in fig. 5a and 5 b.
Considering that the values of the normalized solar pitch angle have a process from-1 to 1 in the upper half and from 1 to-1 in the lower half in 24 hours per day, it is difficult to unify the average data of one day in one coordinate system, so for the convenience of statistics and understanding, the normalized solar pitch angle α is periodically calculated by using trigonometric function, i.e. the normalized solar pitch angle α is calculated
The statistical law of the normalized solar pitch angle after the cyclic processing is shown in fig. 6.
Using the invention to analyze the near-ground atmospheric turbulence intensity C of different months n 2 The data are measured, and the factors influencing sunrise and sunset, such as the longitude and latitude of the normalized sun pitch angle, the longitude and latitude of a data acquisition place, the month and the like are not related, so that the data in different areas and different periods can be quantitatively statistically analyzed.
Claims (2)
1. A statistical method for atmospheric optical turbulence parameters is characterized by comprising the following steps:
acquiring atmospheric optical turbulence parameters of an acquisition place, wherein the atmospheric optical turbulence parameters comprise near-ground atmospheric turbulence intensity and atmospheric coherence length; acquiring longitude and latitude and altitude of a collection place, and calculating a solar pitch angle and a solar azimuth angle corresponding to the collection time of the atmospheric optical turbulence parameters according to an astronomy formula;
the elevation angle values of the sun before and after passing the top are distinguished by adopting a sun azimuth mark, when the sun azimuth angle is 0-180 degrees, the atmospheric optical turbulence parameter is the upper half-day data, and when the sun azimuth angle is 180-360 degrees, the atmospheric optical turbulence parameter is the lower half-day data;
normalizing the value of the pitch angle larger than 0 degree and the maximum value of the pitch angle, and normalizing the value of the pitch angle smaller than 0 degree and the minimum value of the pitch angle;
calculating a statistical analysis rule of the solar pitch angle alpha after normalization of the atmospheric optical turbulence parameters;
and (4) utilizing a trigonometric function to carry out periodicity on the normalized solar pitch angle alpha, obtaining a statistical analysis rule of the atmospheric optical turbulence parameter along with the solar pitch angle alpha after the periodicity, and completing the statistics of the atmospheric optical turbulence parameter.
2. Statistical method of atmospheric optical turbulence parameters according to claim 1, characterized in that the normalized solar pitch angle α is cycled by a trigonometric function, i.e. it is
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010530547.8A CN111666691B (en) | 2020-06-11 | 2020-06-11 | Statistical method for atmospheric optical turbulence parameters |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010530547.8A CN111666691B (en) | 2020-06-11 | 2020-06-11 | Statistical method for atmospheric optical turbulence parameters |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111666691A CN111666691A (en) | 2020-09-15 |
| CN111666691B true CN111666691B (en) | 2022-10-21 |
Family
ID=72386686
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010530547.8A Active CN111666691B (en) | 2020-06-11 | 2020-06-11 | Statistical method for atmospheric optical turbulence parameters |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111666691B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112948352B (en) * | 2021-02-04 | 2022-08-19 | 中国科学院合肥物质科学研究院 | Method for constructing atmospheric optical turbulence space-time characteristics and probabilistic database |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107121712A (en) * | 2017-04-25 | 2017-09-01 | 中国科学院合肥物质科学研究院 | Difference light beam picture moves the method for measuring turbulent flow profile in real time with flicker laser radar |
| CN109635419A (en) * | 2018-12-10 | 2019-04-16 | 长春理工大学 | Laser atmospheric turbulence Transmission characteristics method based on machine learning |
| CN110954506A (en) * | 2019-11-08 | 2020-04-03 | 南昌大学 | Three-parameter comprehensive measurement method for whole-layer atmospheric optical turbulence |
-
2020
- 2020-06-11 CN CN202010530547.8A patent/CN111666691B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107121712A (en) * | 2017-04-25 | 2017-09-01 | 中国科学院合肥物质科学研究院 | Difference light beam picture moves the method for measuring turbulent flow profile in real time with flicker laser radar |
| CN109635419A (en) * | 2018-12-10 | 2019-04-16 | 长春理工大学 | Laser atmospheric turbulence Transmission characteristics method based on machine learning |
| CN110954506A (en) * | 2019-11-08 | 2020-04-03 | 南昌大学 | Three-parameter comprehensive measurement method for whole-layer atmospheric optical turbulence |
Non-Patent Citations (1)
| Title |
|---|
| 西北高原地区光学湍流的观测与分析;陈小威等;《红外与激光工程》;20160525;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111666691A (en) | 2020-09-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Filonchyk et al. | Combined use of satellite and surface observations to study aerosol optical depth in different regions of China | |
| WO2018153143A1 (en) | Method for measuring mudflat elevation by remotely sensed water content | |
| Wang et al. | Prediction of diffuse solar radiation based on multiple variables in China | |
| Kang et al. | A new characterization and classification method for daily sky conditions based on ground-based solar irradiance measurement data | |
| Qie et al. | A review of atmospheric electricity research in China | |
| CN111666691B (en) | Statistical method for atmospheric optical turbulence parameters | |
| Ogbulezie et al. | A review of regression models employed for predicting diffuse solar radiation in North-Western Africa | |
| Xia et al. | Seasonal variations of atmospheric CH4 at Jingdezhen station in Central China: Understanding the regional transport and its correlation with CO2 and CO | |
| Raisal et al. | The moon phases influence on the beginning of astronomical dawn determination in Yogyakarta | |
| KR102074707B1 (en) | Irradiation measuring sensor using Solar Cell | |
| Shi et al. | Temporal and spatial variations of high-impact weather events in China during 1959–2014 | |
| CN103499390B (en) | Based on monitoring and the analytical approach of the sunlight irradiation of macromolecule material aging test | |
| CN114219176A (en) | Multi-source algorithm integrated solar radiation forecasting system and method | |
| Mucomole et al. | Temporal Variability of Solar Energy Availability in the Conditions of the Southern Region of Mozambique | |
| He et al. | Spatial and Temporal Differences in Surface Albedo over Different Underlying Surfaces in the Badain Fijaran Desert, China | |
| Yang et al. | In-situ measurement of CO2 at the Xinglong regional background station over North China | |
| Liu et al. | Spatiotemporal distribution characteristics and variation trends of hierarchical precipitation in Guangdong province over the past 50 years | |
| Xia et al. | Monthly calibration and optimization of Ångström-Prescott equation coefficients for comprehensive agricultural divisions in China | |
| Kaczorowska et al. | Long term performance of a PV system with monocrystalline PV cells—a case study | |
| Wang et al. | Bias correction and variability attribution analysis of surface solar radiation from MERRA-2 reanalysis | |
| Guha et al. | Lightning electrical characteristics during tropical summer thunderstorm in north-east India | |
| CN106708781B (en) | Photovoltaic power station flicker level measurement method and device | |
| Jiang et al. | Characteristics of cloud-to-ground lightning (CG) and differences between+ CG and-CG strokes in China regarding the China National Lightning Detection Network | |
| Bischoff et al. | Lidar error in complex terrain-case study based on two measurement campaigns | |
| Li et al. | Npp estimation of high heterogeneous region based on spatiotemporal fusion |
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 |