CN114675277B - Near-ground atmosphere refractive index profile monitoring method based on commercial microwave return link - Google Patents
Near-ground atmosphere refractive index profile monitoring method based on commercial microwave return link Download PDFInfo
- Publication number
- CN114675277B CN114675277B CN202210305997.6A CN202210305997A CN114675277B CN 114675277 B CN114675277 B CN 114675277B CN 202210305997 A CN202210305997 A CN 202210305997A CN 114675277 B CN114675277 B CN 114675277B
- Authority
- CN
- China
- Prior art keywords
- refractive index
- index profile
- attenuation
- commercial microwave
- ground
- 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
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
-
- 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 application discloses a near-ground atmosphere refractive index profile monitoring method based on a commercial microwave return link, which comprises the following steps: acquiring clear sky actual measurement attenuation of each return link based on the return links of commercial microwaves in a preset area; setting and initializing near-ground atmosphere refractive index profile characterization parameters; obtaining clear-air theoretical attenuation of each link based on the characterization parameters and an electromagnetic propagation calculation model according to the antenna parameters of the commercial microwaves, and the spatial positions of a receiving end and a transmitting end of the link; calculating a clear sky attenuation target function according to the measured clear sky attenuation and the theoretical clear sky attenuation; and obtaining the near-ground atmospheric refractive index profile characterization parameters by using an optimization algorithm based on the clear sky attenuation target function to obtain an optimal solution, and completing the near-ground atmospheric refractive index profile monitoring. The method and the device realize effective monitoring of the refractive index profile of the near-ground atmosphere within the region range, and have the advantages of low cost, high space-time resolution, wide coverage range and the like.
Description
Technical Field
The application belongs to the technical field of near-ground atmospheric environment parameter acquisition, and particularly relates to a near-ground atmospheric refractive index profile monitoring method based on a commercial microwave return link.
Background
The near-ground atmosphere refractive index profile has an important influence on an electromagnetic signal propagation path, and the abnormal atmosphere refractive index vertical gradient often causes that commercial communication equipment cannot work normally, and meanwhile, the over-the-horizon detection of the ground-based radar can also be realized. At present, the monitoring means of the near-ground atmosphere refractive index profile information is limited. The atmospheric temperature, pressure and humidity parameter profiles are telemetered through a ball-mounted radiosonde, and corresponding atmospheric refractive index profile information can be calculated through a theoretical formula. However, in the conventional business, the radiosonde is only flown twice a day at a fixed point, so that fine atmospheric refractive index profile change information cannot be acquired. Furthermore, the vertical resolution of the radiosonde near the ground is also low. The refractometer can directly and accurately measure the atmospheric refractive index, but the instrument is heavy and expensive, and cannot be popularized in a large range.
In addition to the above-mentioned professional atmospheric refractive index detection methods, non-cooperative electromagnetic wave sources that are widely present throughout the world may also be used for the detection of the atmospheric refractive index profile. The most typical example is to calculate the atmospheric refractive index using the amount of atmospheric refraction generated by the GPS signal passing through the troposphere. The method can realize low-cost and high-vertical-resolution global atmospheric refractive index profile monitoring, but is difficult to acquire effective information in a near-ground range. Besides GPS non-cooperative electromagnetic wave sources, a large number of commercial microwave return links widely exist near the ground, the operating frequency of the links is in the range of 6-42 GHz, and the link length is in the range of thousands of meters to dozens of kilometers. The electromagnetic propagation process of the commercial microwave return link is also influenced by the atmospheric refraction profile, so that the signals at the receiving end fluctuate, and therefore, the signals also contain the atmospheric refraction profile information and can be theoretically used for monitoring the atmospheric refraction profile.
Disclosure of Invention
The application provides a near-ground atmosphere refractive index profile monitoring method based on a commercial microwave return link, and the measured attenuation of a multilink in clear sky is calculated according to the power of a transmitting end and a receiving end of the microwave return link in an area; defining characterization parameters of the near-ground atmospheric refractive index profile and setting initial values to realize description of the atmospheric refractive index profile; based on the spatial positions of a receiving end and a transmitting end of a link, combining an electromagnetic propagation model, and calculating the clear-sky theoretical attenuation of the multiple links; and calculating a target function according to the difference between the actual measurement attenuation and the theoretical attenuation, and solving the optimal solution of the near-ground atmosphere refractive index profile characterization parameters through a certain optimization algorithm to realize the effective inversion of the near-ground atmosphere refractive index profile in the region range.
In order to achieve the above purpose, the present application provides the following solutions:
the method for monitoring the near-ground atmosphere refractive index profile based on the commercial microwave return link comprises the following steps:
acquiring clear sky actual measurement attenuation of each commercial microwave return link based on commercial microwave return links in a preset area;
setting and initializing characterization parameters of the near-ground atmosphere refractive index profile;
obtaining clear-sky theoretical attenuation of each commercial microwave return link based on the characterization parameters and the electromagnetic propagation calculation model according to antenna parameters of commercial microwaves, and spatial positions of a link receiving end and a transmitting end;
calculating a clear sky attenuation target function according to the measured clear sky attenuation and the theoretical clear sky attenuation;
and obtaining the near-ground atmosphere refractive index profile characterization parameters by using an optimization algorithm based on the clear sky attenuation objective function to obtain an optimal solution, and completing the near-ground atmosphere refractive index profile monitoring.
Optionally, the method for obtaining the measured attenuation in clear sky of each commercial microwave backhaul link includes:
selecting N commercial microwave return links which stably run in the preset area, and acquiring the transmitting end spatial position, the receiving end spatial position and the antenna parameters of the commercial microwave return links;
respectively calculating the average power of a transmitting end and the average power of a receiving end of each commercial microwave return link according to a preset time interval;
and obtaining the clear sky actual measurement attenuation of each commercial microwave return link according to the average power of the transmitting terminal and the average power of the receiving terminal.
Optionally, the electromagnetic propagation calculation model includes, but is not limited to, a parabolic equation method.
Optionally, the clear-sky theoretical attenuation of each commercial microwave return link is obtained by using a step-by-step wide-angle parabolic equation model based on the characterization parameters and the electromagnetic propagation calculation model according to the antenna parameters of the commercial microwaves, and the spatial positions of the receiving end and the transmitting end of the link.
Optionally, the clear sky attenuation objective function is:
wherein Regularizer (M) is a regularization term.
Optionally, the regularization term includes, but is not limited to, an L1 regularization term.
Optionally, the optimal solution of the near-ground atmosphere refractive index profile characterization parameter is a solution that minimizes the objective function, and the expression is
Optionally, the optimization algorithm includes, but is not limited to, an ant colony algorithm.
The beneficial effect of this application does:
the application discloses a near-ground atmosphere refractive index profile monitoring method based on a commercial microwave return link, which is characterized in that on the basis of extracting attenuation information of the commercial microwave return link, an electromagnetic propagation calculation model is combined to search for optimal atmosphere refractive index profile characterization parameters, so that effective monitoring of the near-ground atmosphere refractive index profile within an area range is realized; meanwhile, the method can realize the acquisition of the near-ground atmosphere refractive index profile information based on the widely existing commercial microwave return link on the ground without investing additional instruments and equipment, and has the advantages of low cost, high space-time resolution, wide coverage range and the like, and has extremely high application value. The method can be applied to actual services as a new method for monitoring the near-ground atmospheric refractive index profile.
Drawings
In order to more clearly illustrate the technical solutions of the present application, the drawings required to be used in the embodiments are briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and it is obvious for those skilled in the art that other drawings can be obtained according to these drawings without inventive labor.
Fig. 1 is a schematic flowchart of a method for monitoring a refractive index profile of a near-ground atmosphere based on a commercial microwave return link according to an embodiment of the present application;
fig. 2 is a schematic diagram illustrating a distribution of a commercial microwave backhaul link according to an embodiment of the present application;
FIG. 3 is a schematic representation of an atmospheric refractive index profile characterizing parameter of an embodiment of the present application;
fig. 4 is a schematic diagram of an inversion result of the atmospheric refractive index profile according to an embodiment of the present application.
Detailed Description
The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
In order to make the aforementioned objects, features and advantages of the present application more comprehensible, the present application is described in further detail with reference to the accompanying drawings and the detailed description.
As shown in fig. 1, a schematic flow chart of a method for monitoring a near-ground atmospheric refractive index profile based on a commercial microwave return link according to an embodiment of the present invention mainly includes the following steps:
firstly, acquiring clear sky actual measurement attenuation of each return link based on the return links of commercial microwaves in a preset area.
In this embodiment, the specific implementation method of this step includes the following 3 steps:
1. selecting 4 commercial microwave return links MBL (minimum Block Limited) stably running in a preset area 1 ,MBL 2 ,MBL 3 ,MBL 4 The spatial position of the transmitting end is S 1 ＝(xs 1 ,ys 1 ,zs 1 ),S 2 ＝(xs 2 ,ys 2 ,zs 2 ),S 3 ＝(xs 3 ,ys 3 ,zs 3 ),S 4 ＝(xs 4 ,ys 4 ,zs 4 ) The spatial position of the receiving end is respectively E 1 ＝(xe 1 ,ye 1 ,ze 1 ),E 2 ＝(xe 2 ,ye 2 ,ze 2 ),E 3 ＝(xe 3 ,ye 3 ,ze 3 ),E 4 ＝(xe 4 ,ye 4 ,ze 4 ) The antenna parameters are respectively A 1 ,A 2 ,A 3 ,A 4 ；
2. According to a preset time interval, calculating the average power TP of the transmitting end of each link 1 ,TP 2 ,TP 3 ,TP 4 Average power RP at the receiving end 1 ,RP 2 ,RP 3 ,RP 4 (ii) a In the present embodiment, 60min is used as the time interval.
3. According to the power difference between the receiving end and the transmitting end, the measured attenuation PL in clear sky of each link is calculated i ＝TP i -RP i (i＝1,2,…,4)。
In practical applications, the information of the near-surface average atmospheric refractive index profile in a certain area may be obtained according to the receiving end power and the transmitting end power recorded by the plurality of commercial microwave return links, and in this embodiment, the distribution of the commercial microwave return links is schematically shown in fig. 2.
Then, as shown in fig. 3, a near-surface atmospheric refractive index profile characterization parameter M = [ k ] is set 1 ,k 2 ,k 3 ,H 1 ,H 2 ]And is initialized to k 1 ＝-10N-unitkm -1 ,k 2 ＝10N-unitkm -1 ,k 3 ＝-10N-unitkm -1 ,H 1 ＝30m,H 2 ＝60m。
Thirdly, according to the antenna parameters obtained in the initial stage, the spatial positions of the receiving end and the transmitting end of the link and the near-ground atmosphere refractive index profile characterization parameters set in the previous step, combining an electromagnetic propagation calculation model and based on a step-by-step wide-angle parabolic equation model, calculating the clear-sky theoretical attenuation PL of each link 01 ,PL 02 ,PL 03 ,PL 04 . In the present embodiment, the electromagnetic propagation calculation model includes, but is not limited toIn the parabolic equation method.
And then, calculating a clear sky attenuation target function according to the measured attenuation of the clear sky and the theoretical attenuation of the clear sky.
In this embodiment, the objective function J (M) is:
where Regularizer (M) is a regularization term, in the present embodiment, the regularization term includes, but is not limited to, an L1 regularization term.
And finally, obtaining a near-ground atmospheric refractive index profile characterization parameter by using an optimization algorithm based on the clear sky attenuation target function to obtain an optimal solution, and completing near-ground atmospheric refractive index profile monitoring, wherein an atmospheric refractive index profile inversion result is shown in fig. 4.
In the embodiment, the optimal solution of the near-ground atmosphere refractive index profile characterization parameter is the solution with the minimum objective function, and the expression is
The optimization algorithm includes, but is not limited to, an ant colony algorithm.
The above-described embodiments are merely illustrative of the preferred embodiments of the present application, and do not limit the scope of the present application, and various modifications and improvements made to the technical solutions of the present application by those skilled in the art without departing from the spirit of the present application should fall within the protection scope defined by the claims of the present application.
Claims (8)
1. The near-ground atmosphere refractive index profile monitoring method based on the commercial microwave return link is characterized by comprising the following steps of:
acquiring clear sky actual measurement attenuation of each commercial microwave return link based on commercial microwave return links in a preset area;
setting and initializing characterization parameters of the near-ground atmosphere refractive index profile;
obtaining clear-sky theoretical attenuation of each commercial microwave return link based on the characterization parameters and the electromagnetic propagation calculation model according to antenna parameters of commercial microwaves, and spatial positions of a link receiving end and a transmitting end;
calculating a clear sky attenuation target function according to the measured clear sky attenuation and the theoretical clear sky attenuation;
and obtaining the near-ground atmosphere refractive index profile characterization parameters by using an optimization algorithm based on the clear sky attenuation objective function to obtain an optimal solution, and completing the near-ground atmosphere refractive index profile monitoring.
2. The method for near-surface atmospheric refractive index profile monitoring based on a commercial microwave return link according to claim 1,
the method for obtaining the measured attenuation in clear sky of each commercial microwave return link comprises the following steps:
selecting N commercial microwave return links which stably run in the preset area, and acquiring the transmitting end spatial position, the receiving end spatial position and the antenna parameters of the commercial microwave return links;
respectively calculating the average power of a transmitting end and the average power of a receiving end of each commercial microwave return link according to a preset time interval;
and obtaining the clear sky actual measurement attenuation of each commercial microwave return link according to the average power of the transmitting terminal and the average power of the receiving terminal.
3. The method according to claim 1, wherein the electromagnetic propagation computational model comprises a parabolic equation method.
4. The method according to claim 1, wherein the clear-sky theoretical attenuation of each commercial microwave return link is obtained by using a step-by-step wide-angle parabolic equation model based on the characterization parameters and the electromagnetic propagation calculation model according to the antenna parameters of commercial microwaves, the spatial positions of the receiving end and the transmitting end of the link.
5. The method according to claim 1, wherein the clear sky attenuation objective function is:
where M represents a near-surface atmospheric refractive index profile characterizing parameter, regularizer (M) is a regularization term.
6. The near-ground atmospheric refractive index profile monitoring method based on a commercial microwave backhaul link according to claim 5, wherein the regularization term comprises an L1 regularization term.
7. The method according to claim 6, wherein the optimal solution of the near-ground atmospheric refractive index profile characterization parameter is a solution that minimizes the objective function, and the expression is expressed as
Wherein M represents a near-surface atmospheric refractive index profile characterizing parameter.
8. The method according to claim 7, wherein the optimization algorithm comprises an ant colony algorithm.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210305997.6A CN114675277B (en) | 2022-03-25 | 2022-03-25 | Near-ground atmosphere refractive index profile monitoring method based on commercial microwave return link |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210305997.6A CN114675277B (en) | 2022-03-25 | 2022-03-25 | Near-ground atmosphere refractive index profile monitoring method based on commercial microwave return link |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114675277A CN114675277A (en) | 2022-06-28 |
| CN114675277B true CN114675277B (en) | 2022-11-04 |
Family
ID=82076365
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210305997.6A Active CN114675277B (en) | 2022-03-25 | 2022-03-25 | Near-ground atmosphere refractive index profile monitoring method based on commercial microwave return link |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114675277B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116338609A (en) * | 2023-05-30 | 2023-06-27 | 中国船舶集团有限公司第七〇七研究所 | AIS data processing method, AIS data processing device, AIS data processing equipment and storage medium |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5686919A (en) * | 1995-06-06 | 1997-11-11 | Jordan; James R. | Process for generating wind profiler data free of fixed ground clutter contamination |
| WO2013148703A1 (en) * | 2012-03-27 | 2013-10-03 | Colorado State University Research Foundation | Robust attenuation correction system for radar reflectivity and differential reflectivity |
| CN103792538A (en) * | 2014-02-24 | 2014-05-14 | 北京航空航天大学 | Atmosphere profile inversion method based on foundation hyperspectral microwave radiometer |
| CN104133216A (en) * | 2014-07-17 | 2014-11-05 | 北京无线电测量研究所 | Method and device for detecting radar acquiring low-altitude wind profiles |
| CN105974433A (en) * | 2016-05-05 | 2016-09-28 | 中国人民解放军陆军军官学院 | Continuous laser based detection method for ground layer aerosol backscattering coefficient profile |
| CN106772386A (en) * | 2016-12-13 | 2017-05-31 | 中国人民解放军理工大学 | One kind is using LPSO algorithms by radar return inverting atmospheric duct method |
| CN112213727A (en) * | 2020-10-15 | 2021-01-12 | 国家卫星气象中心(国家空间天气监测预警中心) | Precipitation correction method of satellite-borne radar based on active and passive microwave combined detection |
| CN113391315A (en) * | 2021-06-11 | 2021-09-14 | 中国人民解放军国防科技大学 | Method for inverting atmospheric waveguide by radar echo data based on electromagnetic wave parabolic equation adjoint mode |
-
2022
- 2022-03-25 CN CN202210305997.6A patent/CN114675277B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5686919A (en) * | 1995-06-06 | 1997-11-11 | Jordan; James R. | Process for generating wind profiler data free of fixed ground clutter contamination |
| WO2013148703A1 (en) * | 2012-03-27 | 2013-10-03 | Colorado State University Research Foundation | Robust attenuation correction system for radar reflectivity and differential reflectivity |
| CN103792538A (en) * | 2014-02-24 | 2014-05-14 | 北京航空航天大学 | Atmosphere profile inversion method based on foundation hyperspectral microwave radiometer |
| CN104133216A (en) * | 2014-07-17 | 2014-11-05 | 北京无线电测量研究所 | Method and device for detecting radar acquiring low-altitude wind profiles |
| CN105974433A (en) * | 2016-05-05 | 2016-09-28 | 中国人民解放军陆军军官学院 | Continuous laser based detection method for ground layer aerosol backscattering coefficient profile |
| CN106772386A (en) * | 2016-12-13 | 2017-05-31 | 中国人民解放军理工大学 | One kind is using LPSO algorithms by radar return inverting atmospheric duct method |
| CN112213727A (en) * | 2020-10-15 | 2021-01-12 | 国家卫星气象中心(国家空间天气监测预警中心) | Precipitation correction method of satellite-borne radar based on active and passive microwave combined detection |
| CN113391315A (en) * | 2021-06-11 | 2021-09-14 | 中国人民解放军国防科技大学 | Method for inverting atmospheric waveguide by radar echo data based on electromagnetic wave parabolic equation adjoint mode |
Non-Patent Citations (2)
| Title |
|---|
| 多波段微波辐射计反演大气温湿廓线性能分析;谭泉 等;《遥感技术与应用》;20150228;第30卷(第1期);170-177 * |
| 降水对Ku/Ka频段星地链路衰减特性的影响研究;咸明皓 等;《电子学报》;20200331;第48卷(第3期);417-725 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114675277A (en) | 2022-06-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109581546B (en) | Rainfall type identification method based on microwave link attenuation and polarization information | |
| US9465106B2 (en) | Robust attenuation correction system for radar reflectivity and differential reflectivity | |
| CN114675277B (en) | Near-ground atmosphere refractive index profile monitoring method based on commercial microwave return link | |
| Grabner et al. | Atmospheric refraction and propagation in lower troposphere | |
| CN109884605B (en) | Method for extracting absorption attenuation and Mie scattering attenuation of cloud rain to radar signals | |
| KR101063318B1 (en) | Position measurement system and method of mobile terminal | |
| Nagao et al. | Study on radio propagation prediction by machine learning using urban structure maps | |
| CN107340519B (en) | Weather radar rainfall detection analysis system and method | |
| CN112733394A (en) | Atmospheric parameter inversion method and device | |
| Pastore et al. | Comparison of atmospheric refractivity estimation methods and their influence on radar propagation predictions | |
| CN112560342A (en) | DNN-based atmospheric waveguide parameter estimation method | |
| JP2012149920A (en) | Precipitation intensity estimation system and precipitation intensity estimation method | |
| CN105301593B (en) | A kind of airbound target Height Estimation method based on double frequency RCS information fusions | |
| Varner et al. | Enhanced RF modeling accuracy using simple minimum mean-squared error correction factors | |
| CN113852433B (en) | Wireless channel shadow fading model prediction method based on computer vision | |
| CN115575914A (en) | Multi-band dual-polarization weather radar observed quantity error quantification method | |
| KR102171276B1 (en) | Method and system for analyzing the propagation path loss for duct atmosphere of aesa radar | |
| Correia et al. | Propagation analysis for wireless sensor networks applied to viticulture | |
| Janco et al. | Rain estimation from smart city’s E-band links | |
| Louf et al. | Real-time monitoring of weather radar network calibration and antenna pointing | |
| CN105933078A (en) | Multi-antenna single-channel radio direction finding system and method having optimal amplitude comparison | |
| Claverie | Vertical refractivity profiles within the marine surface boundary layer | |
| KR102119121B1 (en) | Mrthod for calibrating loss of radio-propagation using clutter-morphology calibrating factor based on measurement data | |
| Vestakis et al. | Multistatic narrowband localization in backscatter sensor networks | |
| CN116184342B (en) | Cloud testing radar data calibration method and system based on multi-radar networking |
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 |