CN107656278A - Based on dense precipitation station Quantitative Precipitation estimating and measuring method - Google Patents
Based on dense precipitation station Quantitative Precipitation estimating and measuring method Download PDFInfo
- Publication number
- CN107656278A CN107656278A CN201710775455.4A CN201710775455A CN107656278A CN 107656278 A CN107656278 A CN 107656278A CN 201710775455 A CN201710775455 A CN 201710775455A CN 107656278 A CN107656278 A CN 107656278A
- Authority
- CN
- China
- Prior art keywords
- mrow
- radar
- region
- precipitation
- msub
- 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.)
- Granted
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
- G01S13/958—Theoretical aspects
-
- 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)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention proposes one kind and is based on dense precipitation station Quantitative Precipitation estimating and measuring method, it is many small polygonal regions by a selected rain detection with radar region division, each region uses its respective Z R relation, when different zones are different precipitation property, avoid using error caused by same Z R relations；When radar data quality occurs wrong, for example, velocity of wave stop caused by radar reflectivity factor it is less than normal etc. situations such as, this method can be according to real time data, and by the adjustment to Z R relations, the Calculation of precipitation deviation to caused by due to radar quality gives certain correct.
Description
Technical field
The present invention relates to weather radar pinch-reflex ion diode field, more particularly to one kind to be estimated based on dense precipitation station Quantitative Precipitation
Method.
Background technology
Radar quantitative estimation precipitation is the important technology reference for closing on early warning in short-term, with reference to radar quantitative estimation precipitation
(QPE), there is the achievement in research with various different idea and technical method, such as calculus of variations, Kalman filtering method and
Probability matched pair technique etc..
At present, the major technique of Hubei Province's radar quantitative estimation precipitation is that the point for being calculated as representing with rainfall measures and with weather
Radar combines for the planar survey of representative, and associated methods (RASIM) are integrated using radar and rainfall gauge real-time synchronization.This method
Theoretic discussion, precipitation estimation technology, quality control, error analysis with evaluation criteria etc. compared with domestic and international other method
With advanced and practicality, and this method is easy to use, easy to spread, has preferable accuracy, has become China
Critical function module in Nowcasting system (SWAN), its extensive use, works well.But the technology is in a thunder
A Z-R relation is used only in up to measurement range, it is impossible to accurately reflect the variation characteristic of Rainfall distribution.
The content of the invention
In view of this, the present invention propose it is a kind of can accurately reflect Rainfall distribution quantified based on dense precipitation station
Precipitation estimation method.
The technical proposal of the invention is realized in this way：Estimated the invention provides one kind based on dense precipitation station Quantitative Precipitation
Survey method, comprises the following steps,
S1, for a selected rain detection with radar region, according to the distribution of the precipitation station covered in the rain detection with radar region,
Using Thiessen polygon method, to the rain detection with radar region, it is divided, and obtains multiple polygonal regions；
S2, calculate the Z-R relations of each polygonal region；
S3, according to the Z-R relations of each polygonal region, calculate the rainfall intensity R in selected rain detection with radar region.
On the basis of above technical scheme, it is preferred that in the step S2, for selected polygonal region S, it covers
The precipitation station of lid has i, calculates the rainfall summation Q of the i precipitation station within the h periods in polygonal region SSh, according to
Lower formula calculates the Z-R relations of the polygonal region, i.e. AShValue,
Wherein, Z is radar reflectivity factor, and R is rainfall intensity, and A is the coefficient of Z-R relations, A=10-2000, b Z-R
The index of relation, b=1.0-2.0；
AShFor selected polygonal region S common coefficient；
ZShFor the selected polygonal region S radar total reflectivity factor；
ZihRadar for the corresponding 0.5 ° of elevation angle within the h periods of the precipitation station i overhead in selected polygonal region S is anti-
Penetrate the rate factor.
It is further preferred that b=1.5.
It is further preferred that in step S3, the rainfall intensity R in selected polygonal region S is calculated according to below equation,
ZjhRadar for the corresponding 0.5 ° of elevation angle within the h periods of the mesh point j overhead in selected polygonal region S is anti-
Penetrate the rate factor.
It is furthermore preferred that in step S3, according to the common coefficient Ash of rain detection with radar polygonal region and each mesh point overhead
0.5 degree of elevation angle reflectivity factor, calculate the rainfall intensity R in all polygonal regions, it is multiple to obtain rain detection with radar region
The rainfall intensity R of polygonal region.
Further preferably, according to the rainfall intensity R of the multiple polygonal regions in rain detection with radar region, rain detection with radar region is drawn
The rainfall intensity R of 1km resolution ratio presses the distribution map of multiple polygonal regions.
The present invention's is had the advantages that based on dense precipitation station Quantitative Precipitation estimating and measuring method relative to prior art：
(1) it is many small polygonal regions by a selected rain detection with radar region division, each region is each using its
From Z-R relations, when different zones are different precipitation property, avoid using error caused by same Z-R relations；
(2) when radar data quality occurs wrong, for example, velocity of wave stop caused by radar reflectivity factor it is less than normal etc.
Situation, this method can be according to real time datas, by the adjustment to Z-R relations, the Calculation of precipitation to caused by due to radar quality
Deviation gives certain correct.
Brief description of the drawings
In order to illustrate more clearly about the embodiment of the present invention or technical scheme of the prior art, below will be to embodiment or existing
There is the required accompanying drawing used in technology description to be briefly described, it should be apparent that, drawings in the following description are only this
Some embodiments of invention, for those of ordinary skill in the art, on the premise of not paying creative work, can be with
Other accompanying drawings are obtained according to these accompanying drawings.
Fig. 1 is the region point for being divided to obtain to Wuhan rain detection with radar region using Thiessen polygon method in embodiment 1
Butut；
Fig. 2 is the distribution map in the Wuhan rain detection with radar region obtained in the embodiment of the present invention 1.
Embodiment
Below in conjunction with the accompanying drawing in embodiment of the present invention, the technical scheme in embodiment of the present invention is carried out clear
Chu, it is fully described by, it is clear that described embodiment only a part of embodiment of the present invention, rather than whole realities
Apply mode.Based on the embodiment in the present invention, those of ordinary skill in the art institute under the premise of creative work is not made
The every other embodiment obtained, belongs to the scope of protection of the invention.
Embodiment 1
Choose validation region based on dense precipitation station Quantitative Precipitation estimating and measuring method of the Wuhan City as the present invention.
Firstly, for Wuhan rain detection with radar region, according to the distribution of the precipitation station of its covering, Thiessen polygon method pair is utilized
It is divided in the rain detection with radar region, obtains 1104 polygonal regions, as shown in Figure 1.
Then, i precipitation station rainfall sum Q in 1 hour in each small polygonal region S is utilizedSh, and for rain
The base reflectivity Z at 0.5 ° of elevation angle of weather radar corresponding to the i overhead of amount stationih, according to following equation, try to achieve each polygon
Region S Z-R relations, that is, AShValue,
Wherein, Z is radar reflectivity factor, and R is rainfall intensity, and A is the coefficient of Z-R relations, b 1.5；
AShFor selected polygonal region S common coefficient；
ZShFor the selected polygonal region S radar total reflectivity factor；
ZihRadar for the corresponding 0.5 ° of elevation angle within the h periods of the precipitation station i overhead in selected polygonal region S is anti-
Penetrate the rate factor.
Finally, the rainfall intensity R in selected polygonal region S is calculated according to below equation,
ZjhRadar for the corresponding 0.5 ° of elevation angle within the h periods of the mesh point j overhead in selected polygonal region S is anti-
Penetrate the rate factor.
Further according to the common coefficient A of rain detection with radar polygonal regionshWith 0.5 degree of elevation angle in each mesh point overhead
Reflectivity factor, the rainfall intensity R in all polygonal regions is calculated, obtains the drop of the multiple polygonal regions in rain detection with radar region
Raininess degree R, finally according to the rainfall intensity R of the multiple polygonal regions in rain detection with radar region, draw rain detection with radar region 1km and differentiate
The rainfall intensity R of rate presses the distribution map of multiple polygonal regions.As shown in Figure 2.
The better embodiment of the present invention is the foregoing is only, is not intended to limit the invention, it is all the present invention's
Within spirit and principle, any modification, equivalent substitution and improvements made etc., it should be included in the scope of the protection.
Claims (6)
1. one kind is based on dense precipitation station Quantitative Precipitation estimating and measuring method, it is characterised in that：Comprise the following steps,
S1, for a selected rain detection with radar region, according to the distribution of the precipitation station covered in the rain detection with radar region, utilize
To the rain detection with radar region, it is divided Thiessen polygon method, obtains multiple polygonal regions；
S2, calculate the Z-R relations of each polygonal region；
S3, according to the Z-R relations of each polygonal region, calculate the rainfall intensity R in selected rain detection with radar region.
2. dense precipitation station Quantitative Precipitation estimating and measuring method is based on as claimed in claim 1, it is characterised in that：The step S2
In, for selected polygonal region S, its precipitation station covered has i, calculates the interior i within the h periods of polygonal region S
The rainfall summation Q of individual precipitation stationSh, the Z-R relations of the polygonal region, i.e. A are calculated according to below equationShValue
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>S</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<msub>
<mi>A</mi>
<mrow>
<mi>S</mi>
<mi>h</mi>
</mrow>
</msub>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>b</mi>
</mfrac>
</mrow>
</msup>
<msup>
<msub>
<mi>Z</mi>
<mrow>
<mi>S</mi>
<mi>h</mi>
</mrow>
</msub>
<mfrac>
<mn>1</mn>
<mi>b</mi>
</mfrac>
</msup>
</mrow>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>S</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>S</mi>
</munderover>
<msup>
<msub>
<mi>Z</mi>
<mrow>
<mi>i</mi>
<mi>h</mi>
</mrow>
</msub>
<mfrac>
<mn>1</mn>
<mi>b</mi>
</mfrac>
</msup>
<mo>&rsqb;</mo>
</mrow>
<mi>b</mi>
</msup>
</mrow>
Wherein, Z is radar reflectivity factor, and R is rainfall intensity, and A is the coefficient of Z-R relations, and A=10-2000, b are Z-R relations
Index, b=1.0-2.0；
AShFor selected polygonal region S common coefficient；
ZShFor the selected polygonal region S radar total reflectivity factor；
ZihFor the radar reflectivity at the corresponding 0.5 ° of elevation angle within the h periods of the precipitation station i overhead in selected polygonal region S
The factor.
3. dense precipitation station Quantitative Precipitation estimating and measuring method is based on as claimed in claim 2, it is characterised in that：B=1.5.
4. dense precipitation station Quantitative Precipitation estimating and measuring method is based on as claimed in claim 2, it is characterised in that：In step S3, root
The rainfall intensity R in selected polygonal region S is calculated according to below equation,
<mrow>
<mi>R</mi>
<mo>=</mo>
<msup>
<msub>
<mi>A</mi>
<mrow>
<mi>s</mi>
<mi>h</mi>
</mrow>
</msub>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>b</mi>
</mfrac>
</mrow>
</msup>
<msup>
<msub>
<mi>Z</mi>
<mrow>
<mi>j</mi>
<mi>h</mi>
</mrow>
</msub>
<mfrac>
<mn>1</mn>
<mi>b</mi>
</mfrac>
</msup>
<mo>,</mo>
</mrow>
ZjhFor the radar reflectivity at the corresponding 0.5 ° of elevation angle within the h periods of the mesh point j overhead in selected polygonal region S
The factor.
5. dense precipitation station Quantitative Precipitation estimating and measuring method is based on as claimed in claim 4, it is characterised in that：In step S3, root
According to the common coefficient A of rain detection with radar polygonal regionshWith the reflectivity factor at 0.5 degree of elevation angle in each mesh point overhead, calculate
Rainfall intensity R in all polygonal regions, obtain the rainfall intensity R of the multiple polygonal regions in rain detection with radar region.
6. dense precipitation station Quantitative Precipitation estimating and measuring method is based on as claimed in claim 5, it is characterised in that：According to rain detection with radar
The rainfall intensity R of the multiple polygonal regions in region, the rainfall intensity R of rain detection with radar region 1km resolution ratio is drawn by multiple polygon
The distribution map in shape region.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710775455.4A CN107656278B (en) | 2017-08-31 | 2017-08-31 | Quantitative precipitation estimation method based on dense rainfall station |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710775455.4A CN107656278B (en) | 2017-08-31 | 2017-08-31 | Quantitative precipitation estimation method based on dense rainfall station |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107656278A true CN107656278A (en) | 2018-02-02 |
| CN107656278B CN107656278B (en) | 2020-09-11 |
Family
ID=61128133
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710775455.4A Active CN107656278B (en) | 2017-08-31 | 2017-08-31 | Quantitative precipitation estimation method based on dense rainfall station |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107656278B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108761576A (en) * | 2018-05-28 | 2018-11-06 | 国网山西省电力公司电力科学研究院 | A kind of X-band weather radar and precipitation station data fusion method and system |
| CN108931774A (en) * | 2018-06-26 | 2018-12-04 | 重庆市气象台 | Convective precipitation based on lightning data identifies examination and test of products method and system |
| CN109799550A (en) * | 2019-03-20 | 2019-05-24 | 北京百度网讯科技有限公司 | Method and apparatus for predicting rainfall intensity |
| CN110895354A (en) * | 2019-12-04 | 2020-03-20 | 中国水利水电科学研究院 | Surface rainfall calculation method based on dynamic adjustment of Thiessen polygon |
| CN112526641A (en) * | 2020-12-10 | 2021-03-19 | 重庆市气象台 | Method, system and equipment for identifying quality of rainfall observed value in real time |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030025628A1 (en) * | 2001-08-06 | 2003-02-06 | Nakagawa Katsuhiro | Method and system for estimation of rainfall intensity in mountainous area |
| CN105759274A (en) * | 2016-04-26 | 2016-07-13 | 南京信息工程大学 | Typhoon attention area radar rainfall estimation method |
-
2017
- 2017-08-31 CN CN201710775455.4A patent/CN107656278B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030025628A1 (en) * | 2001-08-06 | 2003-02-06 | Nakagawa Katsuhiro | Method and system for estimation of rainfall intensity in mountainous area |
| CN105759274A (en) * | 2016-04-26 | 2016-07-13 | 南京信息工程大学 | Typhoon attention area radar rainfall estimation method |
Non-Patent Citations (2)
| Title |
|---|
| 万玉发 等: "长江中游临近预报业务系统（MYNOS）及其应用", 《应用气象学报》 * |
| 卢江涛: "利用天气雷达估测和预报降雨分布的研究", 《中国优秀硕士学位论文全文数据库（电子期刊） 基础科学辑》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108761576A (en) * | 2018-05-28 | 2018-11-06 | 国网山西省电力公司电力科学研究院 | A kind of X-band weather radar and precipitation station data fusion method and system |
| CN108761576B (en) * | 2018-05-28 | 2020-11-13 | 国网山西省电力公司电力科学研究院 | Data fusion method and system for X-band meteorological radar and rainfall station |
| CN108931774A (en) * | 2018-06-26 | 2018-12-04 | 重庆市气象台 | Convective precipitation based on lightning data identifies examination and test of products method and system |
| CN109799550A (en) * | 2019-03-20 | 2019-05-24 | 北京百度网讯科技有限公司 | Method and apparatus for predicting rainfall intensity |
| CN109799550B (en) * | 2019-03-20 | 2022-02-18 | 北京百度网讯科技有限公司 | Method and device for predicting rainfall intensity |
| CN110895354A (en) * | 2019-12-04 | 2020-03-20 | 中国水利水电科学研究院 | Surface rainfall calculation method based on dynamic adjustment of Thiessen polygon |
| CN112526641A (en) * | 2020-12-10 | 2021-03-19 | 重庆市气象台 | Method, system and equipment for identifying quality of rainfall observed value in real time |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107656278B (en) | 2020-09-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107656278A (en) | Based on dense precipitation station Quantitative Precipitation estimating and measuring method | |
| CN106950614B (en) | A kind of region automatic weather station hour rainfall data method of quality control | |
| CN102708307B (en) | Vegetation index construction method applied to city | |
| CN105891833A (en) | Method of identifying warm cloud precipitation rate based on Doppler radar information | |
| CN103197299B (en) | Extraction and quantitative analysis system of weather radar radial wind information | |
| CN103196838B (en) | Hyperspectral remote sensing monitoring method for coastal estuary eutrophication | |
| Esteves et al. | Quantifying shoreline changes along the Sefton coast (UK) and the implications for research-informed coastal management | |
| CN102841385A (en) | Local geomagnetic chart constructing method based on multi-fractal Krigin method | |
| CN103792538A (en) | Atmosphere profile inversion method based on foundation hyperspectral microwave radiometer | |
| Lorente et al. | Accuracy assessment of high frequency radar current measurements in the Strait of Gibraltar | |
| Xie et al. | Calculating NDVI for Landsat7-ETM data after atmospheric correction using 6S model: A case study in Zhangye city, China | |
| CN110261857B (en) | Spatial interpolation method for weather radar | |
| CN101915914A (en) | Lookup table based pixel-by-pixel atmospheric correction method of remote sensing images | |
| CN103868860A (en) | Method for monitoring nitrogen concentration of vegetation canopies in wetland based on hyperspectral vegetation index | |
| CN105759274B (en) | A kind of typhoon pays close attention to area's radar precipitation estimating and measuring method | |
| CN107560599A (en) | A kind of road grade data processing method of feature based point sampling and curve matching | |
| CN102707336A (en) | Novel method of using A-Train series satellite data for synergetic inversion of cloud phase states and cloud parameters | |
| CN105069295A (en) | Assimilation method for satellite and ground rainfall measured values based on Kalman filtering | |
| CN108229803A (en) | A kind of extracting method in feature section, the judgment method of Dangerous Area and system | |
| CN103235890A (en) | System and method for satellite short-time approaching rainfall forecasting | |
| CN104392113B (en) | A kind of evaluation method of COASTAL SURFACE cold reactive antibodies wind speed | |
| CN110288117B (en) | Regional reconstruction method for critical frequency of ionosphere parameters | |
| CN105068075B (en) | A kind of computational methods of strong wind near the ground | |
| CN104198985A (en) | Single antenna direction-finding method | |
| CN112733445A (en) | Large-area-scale soil moisture inversion method based on evapotranspiration vegetation index spatial characteristics |
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 |