CN114924241A - Frequency correction method and system for satellite-borne rainfall measurement radar and ground-based weather radar - Google Patents

Frequency correction method and system for satellite-borne rainfall measurement radar and ground-based weather radar Download PDF

Info

Publication number
CN114924241A
CN114924241A CN202210399295.9A CN202210399295A CN114924241A CN 114924241 A CN114924241 A CN 114924241A CN 202210399295 A CN202210399295 A CN 202210399295A CN 114924241 A CN114924241 A CN 114924241A
Authority
CN
China
Prior art keywords
radar
satellite
borne
precipitation
frequency correction
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.)
Pending
Application number
CN202210399295.9A
Other languages
Chinese (zh)
Inventor
吴琼
商建
尹红刚
陈林
谷松岩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National Satellite Meteorological Center
Original Assignee
National Satellite Meteorological Center
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National Satellite Meteorological Center filed Critical National Satellite Meteorological Center
Priority to CN202210399295.9A priority Critical patent/CN114924241A/en
Publication of CN114924241A publication Critical patent/CN114924241A/en
Pending legal-status Critical Current

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
    • 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
    • G01S7/4052—Means for monitoring or calibrating by simulation of echoes
    • 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/955—Radar or analogous systems specially adapted for specific applications for meteorological use mounted on satellite
    • 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/411—Identification of targets based on measurements of radar reflectivity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/10—Devices for predicting weather conditions
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/18—Testing or calibrating meteorological apparatus
    • 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)
  • Environmental & Geological Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Atmospheric Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses a frequency correction method and a frequency correction system for a satellite-borne precipitation measurement radar and a ground-based weather radar. The frequency correction method comprises the following steps: the first step is as follows: extracting inverted raindrop spectrum parameters from the actual measurement data of the satellite-borne radar: median diameter D of precipitation particle m Concentration factor N w And the corresponding phase and temperature T; the second step is that: combining the lookup table of the median diameter of the precipitation particles in the S waveband with the scattering function, and obtaining the phase sum D through the first step m Calculating f corresponding to the phase state by interpolation z (ii) a The third step: by usingObtained f z Calculating Z e As equivalent radar reflectivity factor ZS-DPR of S band; the fourth step: and correcting the radar reflectivity factor of the Ku waveband of the satellite-borne radar based on the equivalent radar reflectivity factor ZS-DPR. By using the method and the device, the difference of radar reflectivity factors caused by different frequencies can be eliminated, and further the cross validation of the satellite-borne radar and the ground-based radar is carried out, so that the accuracy of weather data is ensured.

Description

Frequency correction method and system for satellite-borne rainfall measurement radar and ground-based weather radar
Technical Field
The invention relates to a frequency correction method for a satellite-borne precipitation measurement radar and a ground-based weather radar, and also relates to a corresponding frequency correction system, belonging to the technical field of satellite remote sensing.
Background
The satellite-borne precipitation measuring radar (referred to as satellite-borne radar for short) is an effective means and an important remote sensor for measuring the precipitation of the global scale. At present, the only internationally on-orbit satellite-borne precipitation measurement radar is a dual-frequency precipitation measurement radar (DPR) carried on a GPM satellite transmitted in 2014, and comprises two working frequency points of Ku and Ka. The satellite-borne radar can cover a plurality of ground-based weather radars (ground-based radars for short) in a short time. The combined application of the satellite-borne radar and the ground-based radar is beneficial to improving the networking application level of the radars and improving the radar detection precision. In order to effectively develop the combined application of the satellite-borne radar and the ground-based radar, the effectiveness and the difference degree of detection data between the satellite-borne radar and the ground-based radar need to be determined, the consistency of the satellite-borne radar and the ground-based radar needs to be analyzed, and technical support is provided for developing the combined application of the satellite-borne radar and the ground-based radar and verifying the detection precision of the satellite-borne radar.
For the detection accuracy verification of the satellite-borne radar and the ground-based radar, the compared physical quantities are radar reflectivity factors. However, even for the same observation target, when radar with different wavelengths is used for detection, the detected radar reflectivity factors are different, so when data of the two radars are compared, the difference caused by the different radar frequency points is also needed to be considered.
Taking the S-band used by the ground-based weather radar and the Ku-band used by the satellite-borne radar as examples, under the rayleigh scattering condition, the echo intensities measured by the radars at the two frequency points are approximately equal, and can be directly compared. However, as the particle size increases, Mie scattering (Mie scattering) effect begins to appear, so that the equivalent radar reflectivity factors of the two frequency points deviate to a certain extent. Therefore, during the comparison process, it is necessary to correct the deviation (referred to in the art as frequency correction) either to the Ku band or to the S band.
The current frequency correction method mainly assumes that the radar reflectivity factor of a Ku-band satellite-borne radar is Z under the same raindrop spectrum distribution N (D) and the same environmental conditions eKu And the radar reflectivity factor of the S-band ground radar is Z s And defines a frequency correction factor M f ：
Figure BDA0003598972830000021
Wherein λ is Ku Is the emission wavelength, lambda, of Ku-band space-borne radar S Is the S wave band ground radar emission wavelength, T is the temperature, K is (m) 2 －1)/(m 2 +2), m is the complex refractive index of the precipitation particle at a specific frequency and a specific temperature in the S wave band and the Ku wave band, K varies with frequency and temperature, σ is the backscattering cross section of the precipitation particle, and D is the diameter of the precipitation particle.
The above equation (1) assumes that the raindrop spectrum satisfies the distribution:
N(D)＝N 0 D μ e -ΛD (2)
wherein the concentration parameter N 0 The scale parameter lambda and the shape factor mu are drop spectrum parameters, and D is the diameter of the precipitation particles.
On the basis, a frequency correction relational expression is established by setting different mu values, different temperatures and rainfall and taking the rainfall as 0.1mm/h as stepping, and a Ku and S waveband reflectivity factor lookup table is generated. After the satellite-borne radar and the ground-based radar are matched, the position closest to the reflectivity factor value of the Ku waveband in the lookup table is searched, and the reflectivity factor frequency of the Ku waveband is corrected to be the radar reflectivity factor of the S waveband corresponding to the position.
Thus, with M in formula (1) f Radar reflectivity factor Z of satellite radar (DPR) in Ku wave band eKu Converting the equivalent radar reflectivity factor ZS-DPR of the S wave band of the equivalent ground radar; or converting the radar reflectivity factor Zs of the S wave band detected by the ground-based radar into an equivalent radar reflectivity factor Z of the Ku wave band of an equivalent satellite-borne radar eKu-CIN 。
However, the conventional frequency correction method has the following problems: (1) only liquid precipitation below the 0-degree layer can be corrected, and precipitation and solid precipitation of a melting layer cannot be corrected; (2) the assumed raindrop spectrum distribution type (exponential distribution) and the raindrop spectrum distribution type of the DPR are not consistent; (3) in building the look-up table, the spectral parameter N is typically assumed 0 And the assumption is that the raindrop spectrum distribution characteristics of the laminar cloud precipitation are met and the raindrop spectrum distribution characteristics of the convection precipitation are not met. The above-mentioned look-up table is therefore only suitable for cloud precipitation and not for convective precipitation.
Disclosure of Invention
The invention aims to provide a frequency correction method for a satellite-borne rainfall measurement radar and a ground-based weather radar.
The invention aims to solve another technical problem of providing a frequency correction system of a satellite-borne precipitation measurement radar and a ground-based weather radar.
In order to achieve the purpose, the invention adopts the following technical scheme:
according to a first aspect of the embodiments of the present invention, there is provided a frequency correction method for a satellite-borne precipitation measurement radar and a ground-based weather radar, including the steps of:
the first step is as follows: extracting inverted raindrop spectrum parameters from the actual measurement data of the satellite-borne radar: median diameter D of precipitation particle m Concentration factor N w And the corresponding phase;
the second step is that: median diameter D of precipitation particles combined with S wave band m And the scattering function f z The phase and median diameter D of the precipitation particles obtained in the first step m And calculating the scattering function f corresponding to the phase state by interpolation z ；
The third step: using the obtained scattering function f z Calculating radar reflectivity factor Z of space-borne radar e As equivalent radar reflectivity factor ZS-DPR of S band;
the fourth step: and correcting the radar reflectivity factor of the Ku waveband of the satellite-borne radar based on the equivalent radar reflectivity factor ZS-DPR.
Preferably, the lookup table of the median diameter of the precipitation particles and the scattering function value is a table of the corresponding relation between the median diameter of the precipitation particles and the scattering function value of the satellite-borne radar with the specific waveband under the condition of different phase values through the Mie scattering model calculation in combination with the bright band model of the satellite-borne radar.
Wherein preferably said look-up table comprises a plurality of different look-up tables at different wavelengths and phases.
Preferably, the look-up table is obtained based on the distribution of the raindrop spectrum as a Gamma distribution.
Wherein preferably, the raindrop spectrum distribution function n (d) is:
N(D)＝N w f(D；Dm)
wherein the content of the first and second substances,
Figure BDA0003598972830000031
wherein r is a Gamma function, μ is generally 3, Nw is a concentration factor, D is a precipitation particle diameter, D m Is the median diameter of the precipitation particle, D m The definition is as follows:
Figure BDA0003598972830000032
wherein preferably D is obtained from said satellite-borne radar m 、N w And corresponding phase state and temperature, and calculating radar reflectivity factor Z of the space-borne radar e Comprises the following steps:
Z e ＝N w f z (D m )
Figure BDA0003598972830000041
wherein λ is the wavelength of Ku band radar, Kw is the complex refractive index of water, T is temperature, σ b Is the backscattering cross-section of the rainfall particle, D is the diameter of the rainfall particle, and Dm is the median diameter of the rainfall particle.
According to a second aspect of the embodiments of the present invention, there is provided a frequency correction system for a satellite-borne precipitation measurement radar and a ground-based weather radar, comprising:
the receiving module is used for receiving data from the satellite-borne radar;
the processing module is connected with the receiving module and is used for processing the data provided by the receiving module according to the frequency correction method of the satellite-borne precipitation measurement radar and the ground-based weather radar;
and the display module is used for connecting the processing module so as to display the calculation result.
Compared with the prior art, the equivalent radar reflectivity factor of the ground radar can be deduced from the raindrop spectrum parameters inverted by the satellite-borne rainfall measurement radar, so that the difference of the radar reflectivity factors caused by different frequencies is eliminated, the cross validation of the satellite-borne radar and the ground radar is further carried out, and the accuracy of weather data is ensured.
Drawings
FIG. 1 is a schematic flow chart of a frequency correction method for a satellite-borne precipitation measurement radar and a ground-based weather radar according to the present invention;
FIG. 2 is a diagram of a bright band model according to an embodiment of the present invention;
FIG. 3 is a schematic representation of the relationship between Dm and fz in the Ku band, in accordance with an embodiment of the present invention;
fig. 4 is a schematic structural diagram of a frequency correction system of the satellite-borne precipitation measurement radar and the ground-based weather radar provided by the invention.
Detailed Description
The technical contents of the invention are described in detail below with reference to the accompanying drawings and specific embodiments.
It has been mentioned previously that the satellite-borne radar (DPR) is currently the only dual-frequency precipitation measuring radar operating on orbit. The dual-frequency joint observation has the advantages that the information of the raindrop spectrum can be more accurately inverted compared with a single frequency, and the rainfall rate is calculated through the raindrop spectrum, so that the rainfall inversion precision is improved. Based on the raindrop spectrum parameters, in addition to calculating the rainfall rate, a radar reflectivity factor may also be calculated. The radar reflectivity factor finally inverted by the satellite-borne radar (DPR) is obtained by calculating the raindrop spectrum parameters.
As shown in fig. 1, the frequency correction method for the satellite-borne precipitation measurement radar and the ground-based weather radar provided by the invention at least comprises the following steps:
the first step is as follows: extracting inverted raindrop spectrum parameter precipitation particle median diameter D from actual measurement data of satellite-borne radar m Concentration factor N w And the corresponding phase.
Raindrop spectrum parameter precipitation particle median diameter D m Concentration factor N w And the corresponding phase inversion process can be performed by conventional techniques, such as J O U N A L O F A T M O S P H E R I C A N D O C E A N I C T E C H N O L O G Y, volume 36, page 883-902, LIANG LIAO, ROBERT MENEGHINI, "Physical Evaluation of GPM DPR Single-and Dual-wavelet Algorithms".
The second step is that: look-up table combining median diameter of precipitation particles in S band with scattering function (D) m -f z Look-up table), the phase sum D obtained by the first step m Calculating f corresponding to the phase state by interpolation z。
Pre-established D m -f z The lookup table is D of the satellite-borne radar with a specific waveband (such as a Ku waveband) under different phase numerical values by combining a bright band model of the satellite-borne radar (DPR) and calculating through a Mie scattering model m And f z A table of correspondence relationships. D of S band, as shown in FIG. 3 m -f z The lookup table is a scattering function f corresponding to the median diameter of the precipitation particles of each step in a specific phase state in the S wave band z (D m ) The value of (c). D m -f z The look-up table comprises a plurality of different look-up tables at different wavelengths and phases.
In the following, D is further described m -f z And (5) establishing a lookup table.
As shown in fig. 2, the radar echo map of a large range of precipitation is divided into ice crystal layer, melting layer and precipitation layer according to the bright band model, and is generally represented by 7 phases (50, 100, 125, 150, 175, 200 and 250). Where 50 represents a solid state, corresponding to a temperature of-50 degrees, 250 represents a liquid state, corresponding to a temperature of 50 degrees, 100, 125, 150, 175, and 200 represent mixed phase states of the melting layer, corresponding to a temperature of 0 degrees.
Median diameter D of precipitation particle m Typically in the range of 0.1-5.0mm, the step size being chosen to be 0.001mm in the present embodiment. It should be noted that the above numerical values are only examples and do not limit the present invention.
Scattering function f z (D m ) The derivation process of (c) is as follows:
the raindrop spectrum observation is one of rainfall physics observations, is used for knowing the micro-physical structures such as number concentration, water content, spectrum distribution, rainfall microstructure and the like of rainfall, and can establish the relation between radar reflectivity factors and rainfall intensity. The raindrop spectrum is usually described by using a function model such as an M2P distribution, a log-normal distribution, a Gamma distribution and the like.
Considering that the raindrop spectrum distribution of the satellite-borne radar (DPR) is a Gamma distribution, the raindrop spectrum distribution of the present invention is assumed to be a Gamma distribution. Therefore, the distribution function n (d) of the raindrop spectrum provided by the present invention is:
N(D)＝N wf (D；Dm) (3)
wherein the content of the first and second substances,
Figure BDA0003598972830000061
r is Gamma function, mu is generally 3, Nw is concentration factor, D is precipitation particle diameter, D m Is the median diameter of the precipitation particle, D m The definition is as follows:
Figure BDA0003598972830000062
d obtained from a satellite-borne radar (DPR) m 、N w And corresponding phase state and temperature, can calculate radar reflectivity factor Z of the satellite-borne radar e Can be expressed as:
Z e ＝N w f z (D m ) (6)
Figure BDA0003598972830000063
wherein λ is the wavelength of Ku band radar, Kw is the complex refractive index of water, T is temperature, σ b Is the backscattering cross section of the rainfall particle, D is the diameter of the rainfall particle, and Dm is the median diameter of the rainfall particle.
Therefore, the scattering function f of the satellite-borne radar in the S band can be calculated based on the formula (7) z (D m ) The value is obtained.
Setting a group of wavelengths lambda, temperatures T and phase states, and pre-calculating according to a formula (7) to obtain the median diameter D of precipitation particles in each step m And corresponding scattering function f z (D m ) Value, stored as a sheet D m -f z A lookup table; different wavelengths and phases have different D m -f z Look-up table, as shown in fig. 3.
The third step: using the obtained f z Calculating Z e As equivalent radar reflectivity factor ZS-DPR for the S band.
Will f is z Substitution into the formula (6) to obtain Z e And the equivalent radar reflectivity factor ZS-DPR of the satellite-borne radar in the S wave band is used.
The fourth step: and correcting the radar reflectivity factor of the Ku waveband of the satellite-borne radar based on the equivalent radar reflectivity factor ZS-DPR.
Through the correction factors obtained in the steps, effective conversion between data acquired by the ground-based radar through the S wave band and data acquired by the Ku wave band used by the satellite-borne radar can be realized, and adverse effects caused by different working frequencies of the satellite-borne precipitation measurement radar and the ground-based weather radar are overcome. On the basis, a lookup table of scattering calculation can be established for different wavelengths and phase states, and the equivalent radar reflectivity factor of the S-band foundation radar including stratums and convection precipitation can be quickly obtained from a raindrop spectrum inverted by the satellite-borne radar.
Based on the raindrop spectrum parameters, the radar reflectivity factor of any frequency point including the S wave band can be calculated, so that frequency correction is completed. However, in real-time processing, if the calculation is performed one by relying on Mie scattering (Mie scattering), the speed is very slow, so that the method adopts a mode of establishing a corresponding lookup table in advance in order to improve the calculation efficiency.
For space-borne radars (DPR), not only D m Differences in temperature and the proportions of water and ice in the melt affect the calculation of the scattering cross-section. Therefore, the invention combines the bright band model to establish different wavelengths, different temperatures and different D m In case of D m -f z And (6) looking up a table. Because the invention takes into account different temperatures and different D m In the case of (3), since it is possible to match the conditions of a plurality of phases such as the melting layer and the solid precipitation, it is possible to correct the precipitation and the solid precipitation in the melting layer, and it is possible to improve the prediction accuracy.
Moreover, the method is based on the assumption that the distribution of the Gamma-distributed raindrop spectrum is consistent with the actual situation of a satellite-borne radar (DPR), so that the correction method provided by the invention is more accurate and the prediction precision is higher. Therefore, the method not only combines the raindrop spectrum distribution type inverted by the satellite-borne radar, but also comprehensively considers different phase states and precipitation types, thereby improving the correction precision.
Due to the fact that
Figure BDA0003598972830000071
c (h) is a correction factor that varies with altitude,
Figure BDA0003598972830000072
wherein, under the different precipitation type conditions: the values of r, p and q are different, and the parameter Nw is related to the type of precipitation. Therefore, the parameter Nw introduced by the method enables the frequency correction method of the satellite-borne precipitation measurement radar and the ground-based weather radar provided by the invention to better accord with the raindrop spectrum distribution characteristics of layer cloud precipitation and convection precipitation. Therefore, the frequency correction method of the satellite-borne precipitation measurement radar and the ground-based weather radar extends from a liquid state to a mixed phase state and a solid state, extends from layer cloud precipitation to the condition of convective precipitation, and is higher in prediction precision and wider in application range.
Furthermore, considering the requirement of the practical process on the operation speed, the invention passes through the phase state and D m To interpolate to determine f z And (4) calculating an equivalent radar reflectivity factor ZS-DPR of the S wave band by using a formula (6) to complete frequency correction. Thus, the calculation speed of real-time correction is improved.
Therefore, the method has good application value in the aspects of combined application of the satellite-borne radar and the ground-based radar and verification of the detection precision of the satellite-borne radar, and not only can serve for networking observation of the ground-based weather radar, but also can directly serve for verification and verification of other satellite-borne rainfall measurement radars.
On the basis of the frequency correction method of the satellite-borne precipitation measurement radar and the ground-based weather radar, the invention further provides a frequency correction system of the satellite-borne precipitation measurement radar and the ground-based weather radar. As shown in fig. 4, the frequency correction system includes a receiving module, a processing module and a display module. The receiving module is used for receiving data from the satellite-borne radar. The processing module can be realized by a single chip microcomputer or a microcontroller. The processing module is connected with the receiving module and is used for processing the data provided by the receiving module according to the frequency correction method of the satellite-borne precipitation measurement radar and the ground-based weather radar. The display module can be realized by an LCD or OLED display screen and is used for being connected with the processing module so as to display the calculation result.
Compared with the prior art, the frequency correction method of the satellite-borne precipitation measurement radar and the ground-based weather radar can deduce corresponding equivalent radar reflectivity factors of the ground-based radar from the raindrop spectrum parameters inverted by the satellite-borne precipitation measurement radar, so that the difference of radar reflectivity factors caused by different frequencies is eliminated, cross validation of the satellite-borne radar and the ground-based radar is further carried out, and the accuracy of weather data can be ensured.
The frequency correction method and system for the satellite-borne rainfall measurement radar and the ground-based weather radar provided by the invention are explained in detail above. It will be apparent to those skilled in the art that any obvious modifications thereof can be made without departing from the spirit of the invention, which infringes the patent right of the invention and bears the corresponding legal responsibility.

Claims (7)

1. A frequency correction method for a satellite-borne rainfall measurement radar and a ground-based weather radar is characterized by comprising the following steps:
the first step is as follows: extracting inverted raindrop spectrum parameters from the actual measurement data of the satellite-borne radar: median diameter D of precipitation particle m Concentration factor N w And a corresponding phase state;
the second step: median diameter D of precipitation particles combined with S wave band m And a scattering function f z The phase and median diameter D of the precipitation particles obtained in the first step m And calculating the scattering function f corresponding to the phase state by interpolation z ；
The third step: using the obtained scattering function f z Calculating radar reflectivity factor Z of space-borne radar e As equivalent radar reflectivity factor ZS-DPR of S band;
the fourth step: and correcting the radar reflectivity factor of the Ku waveband of the satellite-borne radar based on the equivalent radar reflectivity factor ZS-DPR.
2. The frequency correction method of claim 1, wherein:
the lookup table of the median diameter of the precipitation particles and the scattering function value is a table of the corresponding relation of the median diameter of the precipitation particles and the scattering function value of the satellite-borne radar with the specific wave band under the condition of different phase values, which is obtained by combining a bright band model of the satellite-borne radar and calculating through a Mie scattering model.
3. The frequency correction method of claim 2, wherein:
the look-up table comprises a plurality of different look-up tables at different wavelengths and phases.
4. The frequency correction method of claim 1, wherein:
the look-up table is derived for a Gamma distribution based on the raindrop spectral distribution.
5. The frequency correction method of claim 4, wherein:
the raindrop spectral distribution function n (d) is:
N(D)＝N w f(D；Dm)
wherein the content of the first and second substances,
Figure FDA0003598972820000011
wherein r is a Gamma function, μ is generally 3, Nw is a concentration factor, D is a precipitation particle diameter, D m Is the median diameter of the precipitation particle, D m The definition is as follows:
Figure FDA0003598972820000021
6. the frequency correction method of claim 5, wherein:
d obtained from the space-borne radar m 、N w And calculating radar reflectivity factor Z of the satellite-borne radar according to the corresponding phase state and temperature e Comprises the following steps:
Z e ＝N w f z (D m )
Figure FDA0003598972820000022
wherein λ is the wavelength of Ku band radar, Kw is the complex refractive index of water, T is temperature, σ b Is the backscattering cross section of the rainfall particle, D is the diameter of the rainfall particle, and Dm is the median diameter of the rainfall particle.
7. A frequency correction system for a satellite-borne precipitation measurement radar and a ground-based weather radar, comprising:
the receiving module is used for receiving data from the satellite-borne radar;
the processing module is connected with the receiving module and is used for processing the data provided by the receiving module according to the frequency correction method of the satellite-borne precipitation measurement radar and the ground-based weather radar in any one of claims 1 to 6;
and the display module is used for connecting the processing module so as to display the calculation result.
CN202210399295.9A 2022-04-15 2022-04-15 Frequency correction method and system for satellite-borne rainfall measurement radar and ground-based weather radar Pending CN114924241A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210399295.9A CN114924241A (en) 2022-04-15 2022-04-15 Frequency correction method and system for satellite-borne rainfall measurement radar and ground-based weather radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210399295.9A CN114924241A (en) 2022-04-15 2022-04-15 Frequency correction method and system for satellite-borne rainfall measurement radar and ground-based weather radar

Publications (1)

Publication Number Publication Date
CN114924241A true CN114924241A (en) 2022-08-19

Family

ID=82806103

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210399295.9A Pending CN114924241A (en) 2022-04-15 2022-04-15 Frequency correction method and system for satellite-borne rainfall measurement radar and ground-based weather radar

Country Status (1)

Country Link
CN (1) CN114924241A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115575914A (en) * 2022-12-07 2023-01-06 成都远望科技有限责任公司 Multi-band dual-polarization weather radar observed quantity error quantification method
CN115980756A (en) * 2023-03-17 2023-04-18 中国人民解放军国防科技大学 Method for identifying type of condensate in precipitation based on satellite-borne dual-frequency radar
CN116449331A (en) * 2023-06-20 2023-07-18 成都远望科技有限责任公司 Dust particle number concentration estimation method based on W-band radar and meteorological satellite
CN117452360A (en) * 2023-09-26 2024-01-26 北京华云星地通科技有限公司 Matching method and system for satellite-borne radar and ground-based radar data

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115575914A (en) * 2022-12-07 2023-01-06 成都远望科技有限责任公司 Multi-band dual-polarization weather radar observed quantity error quantification method
CN115575914B (en) * 2022-12-07 2023-03-10 成都远望科技有限责任公司 Multi-band dual-polarization weather radar observed quantity error quantification method
CN115980756A (en) * 2023-03-17 2023-04-18 中国人民解放军国防科技大学 Method for identifying type of condensate in precipitation based on satellite-borne dual-frequency radar
CN116449331A (en) * 2023-06-20 2023-07-18 成都远望科技有限责任公司 Dust particle number concentration estimation method based on W-band radar and meteorological satellite
CN116449331B (en) * 2023-06-20 2023-08-15 成都远望科技有限责任公司 Dust particle number concentration estimation method based on W-band radar and meteorological satellite
CN117452360A (en) * 2023-09-26 2024-01-26 北京华云星地通科技有限公司 Matching method and system for satellite-borne radar and ground-based radar data

Similar Documents

Publication Publication Date Title
CN114924241A (en) Frequency correction method and system for satellite-borne rainfall measurement radar and ground-based weather radar
US7733264B1 (en) System and method for generating weather radar information
CN112213727B (en) Precipitation correction method of satellite-borne radar based on active and passive microwave combined detection
Kirstetter et al. Comparison of TRMM 2A25 products, version 6 and version 7, with NOAA/NSSL ground radar–based National Mosaic QPE
Tabary The new French operational radar rainfall product. Part I: Methodology
Gourley et al. A method for evaluating the accuracy of quantitative precipitation estimates from a hydrologic modeling perspective
Lu et al. Characterizing the FY-3A microwave temperature sounder using the ECMWF model
Delrieu et al. Quantification of path-integrated attenuation for X-and C-band weather radar systems operating in Mediterranean heavy rainfall
CN105842692B (en) A kind of atmospheric correction method in INSAR measurements
Arriola et al. Variational bias correction of GNSS ZTD in the HARMONIE modeling system
CN113378443B (en) Ground wave radar data fusion assimilation method and computer equipment
CN112764057B (en) Method and device for monitoring quality of geostationary satellite external measurement data and electronic equipment
Sideris et al. The CombiPrecip experience: development and operation of a real-time radar-raingauge combination scheme in Switzerland
CN115980756A (en) Method for identifying type of condensate in precipitation based on satellite-borne dual-frequency radar
Jiang et al. Calibration and validation of reprocessed hy-2a altimeter wave height measurements using data from buoys, jason-2, cryosat-2, and saral/altika
CN115047499A (en) Inversion method and system for satellite-borne GNSS-R soil temperature and humidity
Samuel et al. Development and evaluation of a hydrologic data-assimilation scheme for short-range flow and inflow forecasts in a data-sparse high-latitude region using a distributed model and ensemble Kalman filtering
JP2000009857A (en) Weather radar device
KR101836135B1 (en) The apparatus and method for integrating rainfall data of multi-satelite
JP3851641B2 (en) Rainfall prediction method and rainfall prediction program
Wen et al. Evaluation of a method to enhance real-time, ground radar–based rainfall estimates using climatological profiles of reflectivity from space
CN109541565B (en) Radar echo intensity uniformity detection method and system
CN109444892B (en) Calibration method for radar quantitative measurement of precipitation
CN112836449B (en) Method for calibrating hydrologic model
Arias et al. Cross validation of the network of ground-based radar with GPM during the Remote sensing of Electrification, Lightning, And Mesoscale/microscale Processes with Adaptive Ground Observations (RELAMPAGO) field campaign

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