CN112131990B - Millimeter wave network rainfall inversion model construction method suitable for complex scene - Google Patents
Millimeter wave network rainfall inversion model construction method suitable for complex scene Download PDFInfo
- Publication number
- CN112131990B CN112131990B CN202010964571.2A CN202010964571A CN112131990B CN 112131990 B CN112131990 B CN 112131990B CN 202010964571 A CN202010964571 A CN 202010964571A CN 112131990 B CN112131990 B CN 112131990B
- Authority
- CN
- China
- Prior art keywords
- link
- millimeter wave
- rainfall
- attenuation
- data
- 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
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/14—Rainfall or precipitation gauges
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/20—Monitoring; Testing of receivers
- H04B17/21—Monitoring; Testing of receivers for calibration; for correcting measurements
Abstract
The invention discloses a millimeter wave network rainfall inversion model construction method suitable for a complex scene, which comprises the following specific steps: firstly, obtaining a signal attenuation value A of each link in a millimeter wave network; secondly, establishing a millimeter wave rainfall inversion model considering a complex scene, and thirdly, for each link, actually measuring attenuation value A according to historyCAnd constructing a historical data set by using the historical rainfall data R of the simultaneous scale rainfall station, and determining model parameters s, d, m, c and tau of each link in the fifth step in a parameter calibration mode. The millimeter wave network rainfall inversion model construction method suitable for the complex scene can be suitable for millimeter wave rainfall inversion in different weather areas and different link conditions, noise processing steps such as wet antenna influence and clutter removal are not needed in the early stage, and the reliability of a training data set is improved.
Description
Technical Field
The invention relates to a millimeter wave network rainfall inversion model construction method suitable for complex scenes, and belongs to the technical field of meteorological element monitoring.
Background
The real-time rainfall monitoring by utilizing the millimeter wave network is a novel rainfall monitoring technology. Currently, inversion models are mainly classified into two types: an ITU rain failure model with a physical foundation and a model established based on a machine learning method. The ITU model does not need previous label data, but has lower inversion accuracy along with different link characteristics and application areas, and cannot be well applied to millimeter wave networks in various areas; the machine learning method is to perform learning training by combining data of the rainfall station with surrounding millimeter wave metadata and attenuation data, and the model obtained by the machine learning method has strong ambiguity because the data set is constructed by using the whole millimeter wave network data, and the inversion accuracy of each link still needs to be improved. On the other hand, the prior processing of the data by the existing millimeter wave rainfall inversion method comprises the steps of dry-wet period discrimination, basic attenuation deduction, wet antenna influence deduction and the like, new errors are introduced by the addition of the steps, and the methods cannot effectively process periodic signals existing in millimeter wave signals.
Disclosure of Invention
The purpose of the invention is as follows: in order to overcome the defects in the prior art, the invention provides the millimeter wave network rainfall inversion model construction method suitable for the complex scene, which can be suitable for millimeter wave rainfall inversion in different climatic regions and different link conditions, and noise processing steps such as wet antenna influence, clutter removal and the like are not required to be performed in the early stage, so that the reliability of a training data set is improved.
The technical scheme is as follows: in order to solve the technical problem, the invention provides a millimeter wave network rainfall inversion model construction method suitable for complex scenes, which comprises the following specific steps:
the method comprises the steps of firstly, obtaining signal strength TSL of each link signal transmitting end and signal strength RSL of a receiving end in a millimeter wave network, preprocessing data, interpolating lost signal strength, removing singular values and interpolating to obtain synchronous data of rainfall stations nearby links;
step two, subtracting the signal strength RSL of the receiving end from the signal strength TSL of the transmitting end of each link after preprocessing to obtain a signal attenuation value A;
thirdly, obtaining the attenuation A of the basic signal by using the minimum value of the sliding windowBLAnd subtracting the base signal attenuation A from the signal attenuation value ABLObtaining a corrected attenuation value ACThe size of the sliding window is set to be tau, and for the t time of the ith link, the specific form is as follows:
by calculating correction values for all T momentsThen, a corrected attenuation value A of the ith link is obtainedC;
Fourthly, the attenuation value A after the millimeter wave link data is correctedCTime resampling is carried out, the time precision is consistent with rainfall data of a rainfall station, and the corrected attenuation value after resampling is recorded as A'CThe attenuation value for the ith link is
Fifthly, respectively establishing a millimeter wave rainfall inversion model considering a complex scene for each link in the millimeter wave network,
for the ith link, i is a positive integer, and the model concrete form is as follows:
wherein R isiThe rainfall intensity of the ith link; si、di、mi、ci、τiFor the ith link model parameter, si、diDepending on the link frequency, polarization, mi、ciRelated to regional environment, tau is a sliding window parameter; a'CCorrecting the ith link obtained in the fourth step by resampling to obtain a basic attenuation amount; l isiIs the ith link length;
sixthly, for each link, actually measuring an attenuation value A 'according to the history'CAnd constructing a historical data set by historical rainfall data R of the rainfall station in the same period, and determining model parameters s, d, m, c and tau of each link in the fifth step in a parameter calibration mode.
Preferably, in the fourth step, the average value of the modified attenuations in the time step is selected as the modified attenuation value at the end of the time period by resampling the modified attenuation data of the millimeter wave link.
Preferably, in the sixth step, the SCE-UA algorithm is used to calibrate the model parameters, and the objective function is the nash efficiency coefficient NSE.
Has the advantages that: according to the millimeter wave network rainfall inversion model construction method suitable for the complex scene, the influence of regional environment elements on millimeter wave attenuation is considered in the model structure, corresponding inversion models can be provided for different regions and different link attributes, the requirement of millimeter wave rainfall inversion under the complex scene can be met, and the method has stronger applicability compared with other technologies; according to the method, the millimeter wave attenuation data are corrected by using the minimum value of the sliding window, and noise is filtered through local attenuation characteristics, so that compared with other technologies, the data processing steps are simplified, the introduction of errors is avoided, the millimeter wave periodic signal change can be effectively removed, and the reliability of a data set is improved; the sliding window parameters adopted by the millimeter wave attenuation data are determined by adopting a calibration mode, and the requirements of objective facts of different link frequencies, lengths and polarization modes of a millimeter wave network are met.
Drawings
FIG. 1 is a schematic flow chart of the present invention.
FIG. 2 is an example of link basis attenuation, attenuation correction, and rainfall comparison determined by the method of the present invention.
Detailed Description
The present invention will be further described with reference to the accompanying drawings.
As shown in fig. 1, a millimeter wave network rainfall inversion model construction method suitable for complex scenes includes the following steps,
the method comprises the steps of firstly, obtaining signal intensity TSL of each link signal transmitting end and signal intensity RSL of a receiving end in a millimeter wave network, collecting microwave data once every 10s, preprocessing the data, interpolating lost signal intensity, removing and interpolating singular values, and performing linear interpolation by using data before and after the lost values; acquiring rainfall intensity data of rainfall stations nearby a link in the same period of 15 min;
step two, subtracting the signal strength RSL of the receiving end from the signal strength TSL of the transmitting end of each link after pretreatment to obtain a signal attenuation value A, and obtaining the attenuation value A of the ith linki=TSLi-RSLi;
Thirdly, obtaining the attenuation of the basic signal by using the minimum value of the sliding windowABLAnd subtracting the base signal attenuation A from the signal attenuation value ABLObtaining a corrected attenuation value ACThe size of the sliding window is set to be τ equal to 91, that is, the minimum of 91 data is used as the base attenuation of the 91 th time, and for the t time of the ith link, the specific form is:
by calculating correction values for all T momentsThen, a corrected attenuation value A of the ith link is obtainedC;
Fourthly, the attenuation value A after the millimeter wave link data is correctedCTime resampling is carried out, the time precision is 15min, the corrected attenuation values within 15min are averaged to be used as the corrected attenuation values at the time of the end of the period and are kept consistent with rainfall data of a rainfall station, and the corrected attenuation values after resampling are recorded as A'C;
And fifthly, respectively establishing a millimeter wave rainfall inversion model considering a complex scene for each link in the millimeter wave network, wherein for the ith link, the model concrete form is as follows:
wherein R isiThe rainfall intensity of the ith link; si、di、mi、ci、τiFor the ith link model parameter, si、diDepending on the link frequency, polarization, mi、ciRelated to regional environment, tau is a sliding window parameter; a'CPassing the ith link obtained in the fourth step through a duplicateThe sampled corrected base attenuation; l isiIs the ith link length;
sixthly, for each link, actually measuring an attenuation value A according to the historyCAnd constructing a data set by using the historical rainfall data R of the time scale rainfall station, carrying out rating on the model parameters s, d, m, c and tau of each link by using an SCE-UA algorithm and combining a rating data set, and selecting a Nash efficiency coefficient NSE by using an objective function.
In the invention, through the fourth step, all parameters in the model in the third step are obtained, and in practice, only the resampled corrected attenuation value A 'needs to be utilized'CAnd obtaining the link rainfall intensity data R. When the millimeter wave link is used for rain measurement, due to the fact that the millimeter wave link is located in different climatic regions and different combination modes of link frequency and length, data of all links cannot be calculated completely by using one formula, no matter an inversion model with a physical basis or a machine learning model, differential modeling is not performed according to the millimeter wave network condition in a complex scene, the inversion model provided by the invention not only considers the attribute of the millimeter wave link, but also introduces the influence of environmental elements, and the millimeter wave link has stronger applicability compared with other models. Before the millimeter wave attenuation data is subjected to rainfall inversion, the steps of dry-wet period discrimination, basic attenuation determination, wet antenna influence removal and the like are often required to be carried out, the processing modes of the steps are different, great uncertainty is brought, unreasonable processing often brings great errors, the existing method uses the same processing mode for each link, the requirement of a complex millimeter wave network cannot be met, in addition, obvious periodicity exists in millimeter wave signals, other methods cannot be effectively solved, in the invention, the millimeter wave attenuation data is corrected by using the minimum value of a sliding window, the minimum value of the sliding window is the minimum value in a period of time, for the millimeter wave attenuation variable caused by rainfall, a random process can be considered in a local range, the minimum value of the random process meets the expectation and the variance is 0, and the attenuation caused by non-rainfall environmental factors such as humidity and wet antennas is not a random variable, but rather a local constant, the sliding window minimum may therefore represent attenuation induced by non-rainfall elements, i.e. the minimum already contains wet antenna effects and other non-rainfallCompared with other technologies, the method simplifies the data processing steps, avoids the introduction of errors, can effectively remove millimeter wave periodic signal changes, and improves the reliability of a data set; for different links, because the link frequency, the length and the polarization mode are different, and the window size is also different, the sliding window parameters adopted by the millimeter wave attenuation data in the invention are determined by adopting a calibration parameter mode, thereby meeting the objective fact requirements of different link frequencies, lengths and polarization modes of the millimeter wave network.
As shown in fig. 2, the invention can effectively filter the periodic attenuation of the millimeter wave signal, and the attenuation obtained by the invention after correction has a very high matching degree with the rainfall data of the actual rainfall station.
The above description is only of the preferred embodiments of the present invention, and it should be noted that: it will be apparent to those skilled in the art that various modifications and adaptations can be made without departing from the principles of the invention and these are intended to be within the scope of the invention.
Claims (3)
1. A millimeter wave network rainfall inversion model construction method suitable for complex scenes comprises the following specific steps:
the method comprises the steps of firstly, obtaining signal strength TSL of each link signal transmitting end and signal strength RSL of a receiving end in a millimeter wave network, preprocessing data, interpolating lost signal strength, removing singular values and interpolating to obtain synchronous data of rainfall stations nearby links;
step two, subtracting the signal strength RSL of the receiving end from the signal strength TSL of the transmitting end of each link after preprocessing to obtain a signal attenuation value A;
thirdly, obtaining the attenuation A of the basic signal by using the minimum value of the sliding windowBLAnd subtracting the base signal attenuation A from the signal attenuation value ABLObtaining a corrected attenuation value ACThe size of the sliding window is set to τ, which is specific to time t of the ith linkThe form is as follows:
by calculating correction values for all T momentsThen, a corrected attenuation value A of the ith link is obtainedC;
Fourthly, the attenuation value A after the millimeter wave link data is correctedCTime resampling is carried out, the time precision is consistent with rainfall data of a rainfall station, and the corrected attenuation value after resampling is recorded as A'C;
Fifthly, respectively establishing a millimeter wave rainfall inversion model considering a complex scene for each link in the millimeter wave network
For the ith link, i is a positive integer, and the model concrete form is as follows:
wherein R isiThe rainfall intensity of the ith link; si、di、mi、ci、τiFor the ith link model parameter, si、diDepending on the link frequency, polarization, mi、ciRelated to regional environment, tau is a sliding window parameter; a'CCorrecting the ith link obtained in the fourth step by resampling to obtain a basic attenuation amount; l isiIs the ith link length;
sixthly, for each link, actually measuring an attenuation value A 'according to the history'CAnd historical rainfall data R of the same-period rainfall station constructs a historical data set which is calibrated through parametersAnd determining the model parameters s, d, m, c and tau of each link in the fifth step.
2. The millimeter wave network rainfall inversion model construction method suitable for complex scenes of claim 1, characterized in that: and in the fourth step, the average value of the attenuation after correction in the time step is selected as the attenuation value after correction at the end of the time period by resampling the attenuation data after the millimeter wave link correction.
3. The millimeter wave network rainfall inversion model construction method suitable for complex scenes of claim 1, characterized in that: and in the sixth step, the SCE-UA algorithm is utilized to calibrate the model parameters, and the objective function selects a Nash efficiency coefficient NSE.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010964571.2A CN112131990B (en) | 2020-09-15 | 2020-09-15 | Millimeter wave network rainfall inversion model construction method suitable for complex scene |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010964571.2A CN112131990B (en) | 2020-09-15 | 2020-09-15 | Millimeter wave network rainfall inversion model construction method suitable for complex scene |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112131990A CN112131990A (en) | 2020-12-25 |
CN112131990B true CN112131990B (en) | 2021-05-14 |
Family
ID=73845076
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010964571.2A Active CN112131990B (en) | 2020-09-15 | 2020-09-15 | Millimeter wave network rainfall inversion model construction method suitable for complex scene |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112131990B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113240148B (en) * | 2021-01-27 | 2022-02-08 | 河海大学 | Precipitation type prediction method based on millimeter wave data layered intelligent sampling |
CN113114388B (en) * | 2021-04-13 | 2022-02-15 | 中国人民解放军国防科技大学 | Antenna wetting loss determination method based on level time sequence dynamic memory |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101543400A (en) * | 2009-04-23 | 2009-09-30 | 深圳先进技术研究院 | Aanimal behavior detection and automatic analysis system and animal behavior analysis methods |
CN110031916A (en) * | 2019-03-07 | 2019-07-19 | 中国人民解放军国防科技大学 | Rainfall intensity measurement method based on satellite-ground link attenuation effect |
CN110133602A (en) * | 2018-02-08 | 2019-08-16 | 英飞凌科技股份有限公司 | Radar sensing with phasing |
CN110163472A (en) * | 2019-04-11 | 2019-08-23 | 中国水利水电科学研究院 | A wide range of extreme drought emergency monitoring and impact evaluation method and system |
CN110490228A (en) * | 2019-07-15 | 2019-11-22 | 中山大学 | A kind of Hydro-Model Parameter Calibration Technology dynamic rating method based on CPP |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109255100B (en) * | 2018-09-10 | 2020-07-28 | 河海大学 | Urban rainfall inversion method based on microwave attenuation characteristic response fingerprint identification |
-
2020
- 2020-09-15 CN CN202010964571.2A patent/CN112131990B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101543400A (en) * | 2009-04-23 | 2009-09-30 | 深圳先进技术研究院 | Aanimal behavior detection and automatic analysis system and animal behavior analysis methods |
CN110133602A (en) * | 2018-02-08 | 2019-08-16 | 英飞凌科技股份有限公司 | Radar sensing with phasing |
CN110031916A (en) * | 2019-03-07 | 2019-07-19 | 中国人民解放军国防科技大学 | Rainfall intensity measurement method based on satellite-ground link attenuation effect |
CN110163472A (en) * | 2019-04-11 | 2019-08-23 | 中国水利水电科学研究院 | A wide range of extreme drought emergency monitoring and impact evaluation method and system |
CN110490228A (en) * | 2019-07-15 | 2019-11-22 | 中山大学 | A kind of Hydro-Model Parameter Calibration Technology dynamic rating method based on CPP |
Non-Patent Citations (2)
Title |
---|
Extraordinary capacitive deionization performance of highly-ordered mesoporous carbon nanopolyhedra for brackish water desalination;杨涛 等;《Environmental Science Nano》;20190228;第1-2页 * |
基于微波链路的降雨场反演方法研究;姜世泰 等;《物理学报》;20130815;第62卷(第15期);第1-8页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112131990A (en) | 2020-12-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112131990B (en) | Millimeter wave network rainfall inversion model construction method suitable for complex scene | |
CN109581546B (en) | Rainfall type identification method based on microwave link attenuation and polarization information | |
CN110650436B (en) | WiFi data-based position fingerprint database establishing and fitting method | |
CN113139700B (en) | River flow prediction method, device, equipment and storage medium | |
CN116546431B (en) | Beidou all-network communication-based multi-network fusion data acquisition communication system and method | |
CN110543893A (en) | Microwave attenuation precipitation particle type identification method based on BP neural network | |
CN113281754B (en) | WRF-Hydro key parameter calibration method for quantitatively estimating rainfall by integrating rainfall station with radar | |
CN112131989B (en) | Millimeter wave rain measurement model parameter obtaining method based on space rainfall data | |
CN117053838A (en) | Ranging result calibration method for ranging sensor | |
CN110111377B (en) | Evaluation method for risk of regional landslide after earthquake by considering earthquake displacement field | |
CN116776094A (en) | Crystal oscillator temperature test data intelligent analysis storage system | |
CN114675277B (en) | Near-ground atmosphere refractive index profile monitoring method based on commercial microwave return link | |
CN114091362B (en) | Water level prediction method, system, storage medium and equipment | |
CN105338462A (en) | Implementation method for reproducing hearing-aid insertion gain | |
CN115656626A (en) | New energy station terminal dynamic frequency space-time correlation analysis method | |
CN112700039B (en) | Steady state detection and extraction method for load operation data of thermal power plant | |
CN114679216A (en) | Algorithm and device for improving analysis precision based on OTDR | |
CN114157023A (en) | Distribution transformer early warning information acquisition method | |
CN113077083B (en) | Method for monitoring pollution condition based on satellite remote sensing technology | |
CN113364541B (en) | Wireless signal main path transmission delay measurement method based on multi-frequency point signal intensity | |
CN112444705A (en) | Regression correction method for wavelet transformation fault location | |
CN113114388B (en) | Antenna wetting loss determination method based on level time sequence dynamic memory | |
CN116761202B (en) | Network stability monitoring maintenance system based on 5G communication technology | |
US9316715B2 (en) | Process and system for calibrating a first loop feature value estimation method using a first locally measurable loop characteristic and a first set of parameters | |
CN117574305B (en) | Real-time monitoring method and system for running state of equipment |
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 |