CN104820222B - The real-time dynamic attenuation of rainfall radar is corrected and precipitation computational methods - Google Patents
The real-time dynamic attenuation of rainfall radar is corrected and precipitation computational methods Download PDFInfo
- Publication number
- CN104820222B CN104820222B CN201410815601.8A CN201410815601A CN104820222B CN 104820222 B CN104820222 B CN 104820222B CN 201410815601 A CN201410815601 A CN 201410815601A CN 104820222 B CN104820222 B CN 104820222B
- Authority
- CN
- China
- Prior art keywords
- mrow
- radar
- precipitation
- rainfall
- time
- 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
-
- 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
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
-
- 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)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
Corrected and precipitation computational methods the invention discloses a kind of real-time dynamic attenuation of rainfall radar.Pass through the precipitation particles information obtained in real time and the improvement to conventional droop correction method, obtain dynamic decay in real time and correct the relevant parameter calculated with precipitation, be effectively improved because the influence of there and then air, precipitation attenuation and other meteorological factors etc. estimates the error brought to radar precipitation.
Description
Technical field
The invention belongs to atmospheric science field, and in particular to laser raindrop spectrograph and rainfall radar (precipitation
Radar synchro measure), carries out in real time, dynamically decaying to rainfall radar and corrects and precipitation calculating.
Background technology
Since this century, radar quantitative measurement of rainfall obtains the great attention of many countries, particularly developed country again.
The CASA plans that the U.S. has carried out since 2003, are carried out with 4 X-band radar networkings and trigger the precipitation of mountain torrents, urban waterlogging to determine
Measurement；Japan in 2013 proposed 26 " radar rain gauge ", its emphasis flood of development covering by core of radar
The quantitative precipitation measurement in area.These radars are mainly used in the measurement of precipitation, with obvious rainfall Radar Objective Characteristics.In recent years, I
Relevant department of state has also carried out similar work with basin in urban.
Rainfall radar can realize the Rainfall estimation of high-spatial and temporal resolution, and conventional rainfall gauge spatial representative can be overcome very poor
The problem of, but, because rainfall Radar Rainfall Estimation is by the decay in itself of the air of there and then, precipitation and meteorological factor etc.
Influence, resulting echo or radar reflectivity factor (commonly referred to as Z values) are relative value, cause the precipitation error of quantitative estimation
It is larger.To reduce error, both at home and abroad many scientific and technical personnel by precipitation drop-size distribution parting and to type of precipitation parting, count compared with
To approach the coefficient of real radar attenuation and precipitation estimation, but these relatively-stationary statistical values be difficult to comply fully with the world it is each
Each actual Precipitation Process on ground.
With widely using for laser raindrop spectrograph so that true Z values are obtained with rainfall radar and are possibly realized.Laser raindrop
Spectrometer can obtain the particle diameter and falling speed of each precipitation particles in real time, can be directly obtained and worked as according to raininess I and Z value calculation formula
When local I and Z values.Under conditions of time-space registration, I the and Z values measured with laser raindrop spectrograph are carried out real to rainfall radar
When dynamic decay correct and calculated with rainfall, be remarkably improved the accuracy of radar quantitative measurement of rainfall.
The content of the invention
With the progress of radar and raindrop spectrometer technology, the present invention proposes one kind using laser raindrop spectrograph and rainfall thunder
The synchro measure reached, the method that real-time, dynamic decay is corrected and precipitation is calculated is carried out to rainfall radar, to solve single radar
The problem of quantitative estimation precipitation error is larger.
It is an object of the invention to propose that a kind of real-time dynamic attenuation of rainfall radar is corrected and precipitation computational methods.
The real-time dynamic attenuation of the present invention, which is corrected to calculate with precipitation, to be comprised the following steps：
1) the radar reflectivity factor Zp of laser raindrop spectrograph and rainfall radar synchronous real-time measurement on space-time is obtained
(R)。
2) the decay correction method (abbreviation HB methods) proposed according to HITSCHFELD and BORDAN, c, d are solved according to (1) formula
Coefficient,
Wherein, Zm (R) is the radar reflectivity factor with Zp (R) time-space registration of rainfall radar surveying, and R is distance.
3) by known c, d coefficient, decay is carried out by storehouse to every scan line of rainfall radar by (2) formula and corrected,
Wherein, Zt (N) is radar reflectivity factor value of the rainfall radar through revised n-th storehouse of decaying, Zm (N) rainfall
Radar attenuation correct before n-th storehouse radar reflectivity factor value, N be storehouse number, dN be storehouse it is long, J be total storehouse number.
4) radar reflectivity that time laser raindrop spectrograph is matched in time with rainfall radar when obtaining multiple or many in real time
Factor Z p (R) and corresponding raininess Ip (R).
5) a, b coefficient are solved according to (3) formula,
Zpi(R)=a × Ipi(R)b (3)
Wherein, in the case of when many time, i is the when number of times more than 2；In the case of multiple laser raindrop spectrographs, i is
Laser raindrop spectrograph number more than 2.
6) by a, b coefficient tried to achieve, raininess calculating is carried out by storehouse to every scan line of rainfall radar by (4) formula,
Zt (N)=a × It (N)b (4)
Wherein, It (N) is the raininess value in n-th storehouse on rainfall radar scanning line.
7) decay is corrected to calculate with precipitation and terminated.
Advantages of the present invention：
One is that the calculating to HB methods has carried out larger improvement, according to raindrop spectrometer actual measured value, radar is carried out real-time
Dynamic decay is corrected, and is improved a Precipitation Process and is fixed the error that one group of c, d coefficient tape comes.
Two be the Z values according to the actual measurement of raindrop spectrometer, directly carries out decay and corrects, and reduces calculating or fitting decay system
The uncertainty that number k may be brought.
Three be the real-time measurement values for the Z-I relations and rainfall radar set up in real time using raindrop spectrometer, carries out dynamic in real time
Precipitation is calculated, and is improved a Precipitation Process and is fixed the error that one group of a, b coefficient tape comes.
Four be to utilize (1), the processing of (2) formula, and not only decay is corrected, and also the change of radar constant can be carried out
Indirect calibration.
Brief description of the drawings
Fig. 1 corrects the flow chart with precipitation computational methods for the real-time dynamic attenuation of rainfall radar of the present invention；
Embodiment
Embodiments of the present invention are described in detail with reference to Figure of description.
As shown in figure 1, the real-time dynamic attenuation of rainfall radar of the present invention is corrected and comprised the following steps with precipitation computational methods：
1. obtain the radar reflectivity factor Zm (R) of certain period rainfall scope scanning.
2. read the radar reflectivity factor Zp (R) and raininess Ip (R) with period multiple laser raindrop spectrograph measurements.
3. spatially being matched to Zm (R) and Zp (R), and calculate c, d coefficient in (1) formula with seasonal change.
4. bringing c, d coefficient into (2) formula, rainfall radar is carried out to correct by the real-time decay in storehouse.
5. being matched on space-time to Zp (R) and Ip (R), a, b system with seasonal change are calculated or are fitted according to (3) formula
Number.
6. bringing a, b coefficient into (4) formula, rainfall radar is carried out to calculate by the real-time raininess in storehouse.
Terminate 7. decay is corrected to calculate with precipitation.
It is finally noted that, the purpose for publicizing and implementing mode is that help further understands the present invention, but ability
The technical staff in domain is appreciated that：Without departing from the spirit and scope of the invention and the appended claims, it is various replacement and
Modification is all possible.Therefore, the present invention should not be limited to embodiment disclosure of that, the scope of protection of present invention with
The scope that claims are defined is defined.
Claims (4)
1. a kind of real-time dynamic attenuation of rainfall radar correct with precipitation computational methods, be to utilize laser raindrop spectrograph and rainfall radar
Synchro measure, carry out that dynamic decay is corrected in real time and precipitation is calculated to rainfall radar, it is characterised in that including following step
Suddenly：
1) the radar reflectivity factor Zp (R) of laser raindrop spectrograph and rainfall radar synchronous real-time measurement on space-time is obtained；
2) the decay correction method (abbreviation HB methods) proposed according to HITSCHFELD and BORDAN, c, d system are solved according to (1) formula
Number,
<mrow>
<mi>Z</mi>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>Z</mi>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mi>c</mi>
<mi>d</mi>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>R</mi>
</msubsup>
<msup>
<mi>Z</mi>
<mi>d</mi>
</msup>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>R</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>R</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mi>d</mi>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, Zm (R) is the radar reflectivity factor with Zp (R) time-space registration of rainfall radar surveying, and R is distance；
3) by known c, d coefficient, decay is carried out by storehouse to every scan line of rainfall radar by (2) formula and corrected,
<mrow>
<mi>Z</mi>
<mi>t</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>Z</mi>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<mi>c</mi>
<mi>d</mi>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>N</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>J</mi>
</munderover>
<msup>
<mi>Z</mi>
<mi>d</mi>
</msup>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mi>d</mi>
<mi>N</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mi>d</mi>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, Zt (N) is radar reflectivity factor value of the rainfall radar through revised n-th storehouse of decaying, Zm (N) rainfall radar
The radar reflectivity factor value in the n-th storehouse decayed before correcting, N is storehouse number, and dN is that storehouse is long, and J is total storehouse number；
4) radar reflectivity factor that time laser raindrop spectrograph is matched in time with rainfall radar when obtaining multiple or many in real time
ZpiAnd corresponding raininess Ip (R)i(R)；
Wherein, in the case of when many time, number of times when subscript i is；In the case of multiple laser raindrop spectrographs, subscript i is sharp
Light raindrop spectrometer number；
5) a, b coefficient are solved according to (3) formula,
2pi(R)=a × Ipi(R)b (3)
Wherein, in the case of when many time, number of times when subscript i is；In the case of multiple laser raindrop spectrographs, subscript i is sharp
Light raindrop spectrometer number；
6) by a, b coefficient tried to achieve, raininess calculating is carried out by storehouse to every scan line of rainfall radar by (4) formula,
Zt (N)=a × lt (N)b (4)
Wherein, It (N) is the raininess value in n-th storehouse on rainfall radar scanning line；
7) decay is corrected to calculate with precipitation and terminated.
2. the real-time dynamic attenuation of rainfall radar as claimed in claim 1 is corrected and precipitation computational methods, it is characterised in that described
Rainfall radar includes the radar that can be used for S, C, the various models of X-band and the various storehouse length of survey rain.
3. the real-time dynamic attenuation of rainfall radar as claimed in claim 1 is corrected and precipitation computational methods, it is characterised in that described
Rainfall radar includes the radar surveying of list portion and the networking of multi-section radar is measured.
4. the real-time dynamic attenuation of rainfall radar as claimed in claim 1 is corrected and precipitation computational methods, it is characterised in that described
Decay is corrected to calculate with precipitation and is applied among the effective radius of investigation of various radars more than 30km.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410815601.8A CN104820222B (en) | 2014-12-25 | 2014-12-25 | The real-time dynamic attenuation of rainfall radar is corrected and precipitation computational methods |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410815601.8A CN104820222B (en) | 2014-12-25 | 2014-12-25 | The real-time dynamic attenuation of rainfall radar is corrected and precipitation computational methods |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104820222A CN104820222A (en) | 2015-08-05 |
| CN104820222B true CN104820222B (en) | 2017-08-29 |
Family
ID=53730563
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410815601.8A Active CN104820222B (en) | 2014-12-25 | 2014-12-25 | The real-time dynamic attenuation of rainfall radar is corrected and precipitation computational methods |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104820222B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106324580B (en) * | 2016-11-02 | 2018-11-06 | 中国人民解放军理工大学 | A kind of radar return decaying correction method based on microwave links road network |
| FR3061775B1 (en) * | 2017-01-09 | 2020-03-27 | Novimet | METHOD OF ESTIMATING A RAINFALL, IN PARTICULAR RAIN OR SNOW |
| CN109614744B (en) * | 2018-12-28 | 2022-11-22 | 浙江理工大学 | Big data-based precipitation detection method and system |
| CN110018479A (en) * | 2019-04-28 | 2019-07-16 | 中国气象局广州热带海洋气象研究所 | C-band dual-polarization weather radar reflectivity terrain shading decaying correction method |
| CN113252995B (en) * | 2021-05-13 | 2022-05-31 | 中国人民解放军国防科技大学 | Fusion layer attenuation determination method based on ground rainfall intensity weighting |
| CN113985378B (en) * | 2021-12-30 | 2022-03-11 | 雷象科技（北京）有限公司 | Array weather radar attenuation correction method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101535834A (en) * | 2006-11-17 | 2009-09-16 | 罗伯特·博世有限公司 | Method and device for detecting precipitation using radar |
| CN102353955A (en) * | 2010-06-04 | 2012-02-15 | 罗伯特·博世有限公司 | Radar sensor and method for detecting precipitation using a radar senor |
| CN102662162A (en) * | 2012-02-16 | 2012-09-12 | 邓勇 | Precipitation estimation method using X-band Doppler radar |
| CN103499813A (en) * | 2013-10-11 | 2014-01-08 | 邓勇 | Method for distinguishing cloud and precipitation by radar precipitation factor |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7353121B2 (en) * | 2006-04-28 | 2008-04-01 | Smart Disaster Response Technologies, Inc. | Methods, apparatus, media and signals for facilitating real-time management of a utility supply |
-
2014
- 2014-12-25 CN CN201410815601.8A patent/CN104820222B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101535834A (en) * | 2006-11-17 | 2009-09-16 | 罗伯特·博世有限公司 | Method and device for detecting precipitation using radar |
| CN102353955A (en) * | 2010-06-04 | 2012-02-15 | 罗伯特·博世有限公司 | Radar sensor and method for detecting precipitation using a radar senor |
| CN102662162A (en) * | 2012-02-16 | 2012-09-12 | 邓勇 | Precipitation estimation method using X-band Doppler radar |
| CN103499813A (en) * | 2013-10-11 | 2014-01-08 | 邓勇 | Method for distinguishing cloud and precipitation by radar precipitation factor |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104820222A (en) | 2015-08-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104820222B (en) | The real-time dynamic attenuation of rainfall radar is corrected and precipitation computational methods | |
| CN105629263B (en) | A kind of troposphere atmosphere delay estimation error correcting method and correction system | |
| CN106547036B (en) | A kind of united Regional Precipitation measurement method of multiband microwave link | |
| CN102539336B (en) | Method and system for estimating inhalable particles based on HJ-1 satellite | |
| CN103994976B (en) | Agriculture Drought remote-sensing monitoring method based on MODIS data | |
| CN102540166B (en) | Cross radiation calibration method based on optimization algorithm of hyper-spectral sensor | |
| CN104111467B (en) | Network real time kinematic (RTK) instant locating method based on big dipper tri-band wide-lane combination | |
| CN108446999B (en) | Method for estimating ET (energy-efficiency) of different crops in irrigation area based on crown air temperature difference and remote sensing information | |
| CN113177744A (en) | Urban green land system carbon sink amount estimation method and system | |
| CN101852623A (en) | On-track calibration method for internal element of satellite optical remote sensing camera | |
| CN104965207A (en) | Method for acquiring area troposphere zenith delay | |
| CN105203023A (en) | One-stop calibration method for arrangement parameters of vehicle-mounted three-dimensional laser scanning system | |
| CN103901420A (en) | Method for dynamic threshold method remote sensing data cloud identification supported by prior surface reflectance | |
| CN110389087B (en) | PM2.5 concentration satellite remote sensing estimation method in polluted weather | |
| CN102707336A (en) | Novel method of using A-Train series satellite data for synergetic inversion of cloud phase states and cloud parameters | |
| CN107153186A (en) | Laser radar scaling method and laser radar | |
| CN114564767A (en) | Under-cloud surface temperature estimation method based on sun-cloud-satellite observation geometry | |
| CN110927120A (en) | Early warning method for coverage degree of planting | |
| KR101423278B1 (en) | System for calculating rainrate using Local Gauge Correction and method thereof | |
| CN107656278A (en) | Based on dense precipitation station Quantitative Precipitation estimating and measuring method | |
| CN116519913A (en) | GNSS-R data soil moisture monitoring method based on fusion of satellite-borne and foundation platform | |
| CN103197303A (en) | Earth surface two-direction reflection characteristic retrieval method and earth surface two-direction reflection characteristic retrieval system based on multiple sensors | |
| CN101846512B (en) | Method for implementing land patrol | |
| CN107007998A (en) | A kind of image gridding direct range mensuration of golf course | |
| CN103674904A (en) | Rapid atmospheric transmission correction method in infrared characteristic measurement |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| DD01 | Delivery of document by public notice |
Addressee: Deng Yong Document name: Notification of Decision on Request for Restoration of Right |
|
| DD01 | Delivery of document by public notice |
Addressee: Deng Yong Document name: Notification of Passing Preliminary Examination of the Application for Invention |
|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| DD01 | Delivery of document by public notice |
Addressee: Deng Yong Document name: Notification of Publication and of Entering the Substantive Examination Stage of the Application for Invention |
|
| DD01 | Delivery of document by public notice |
Addressee: Deng Yong Document name: Notification of Passing Examination on Formalities |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |