搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于星地链路的垂直降雨场反演方法

咸明皓 刘西川 印敏 宋堃 高太长

引用本文:
Citation:

基于星地链路的垂直降雨场反演方法

咸明皓, 刘西川, 印敏, 宋堃, 高太长

Inversion of vertical rainfall field based on earth-space links

Xian Ming-Hao, Liu Xi-Chuan, Yin Min, Song Kun, Gao Tai-Chang
PDF
HTML
导出引用
  • 在分析星地链路几何结构及传播模型的基础上, 研究基于联合代数重建技术的星地链路反演二维垂直降雨场的方法. 利用实测降雨资料构建3类降雨场, 并搭建3条17 GHz垂直极化星地链路进行数值仿真. 实验结果表明: 单星地链路无法实现二维垂直降雨场的重构, 反演场与真实场的相关系数分别为0.556, 0.504和0.364; 基于双星地链路的反演结果和真实场的相关系数均高于0.98, 平均偏差分别为0.122, 0.159和0.537 mm/h, 欧式距离均低于0.9 mm/h, 熵的相对误差均小于1.6%, 基本实现了二维垂直降雨场的反演; 三星地链路的应用进一步提升了反演精度, 相关系数接近于1, 能够精准重构降雨场. 实验结果验证了基于星地链路的垂直降雨场反演方法的可行性、准确性和有效性. 星地链路的架设和维护简单, 探测降雨范围广, 易于在高原、山区、海岛等传统地面观测资料缺失的地区使用, 可以作为已有降水测量手段的补充.
    The accurate measurement of rain intensity and its distribution in vertical direction can not only help to understand the process of rainfall development, but also play an important role in human life such as agriculture, weather forecasting, water resources management, and natural disaster warning. According to the analysis of the geometric structure of earth-space link and propagation model of electromagnetic wave in atmosphere, in this paper we propose a method to reconstruct two-dimensional(2D) vertical rainfall field by using earth-space links. Firstly, the measured data of micro rain radar (MRR) from Nanjing are used to generate three real vertical rainfall fields which are marked as I, II and III respectively. Secondly, based on the analysis of the earth-space link’s geometry and the effect of signal attenuation from other factors such as scintillation, atmosphere gas and cloud, the vertical rainfall field inversion model is established. According to the power-law relationship between rain intensity and rain attenuation, which is given by International Telecommunication Union (ITU), the simultaneous algebraic reconstruction technique (SART) is used to inverse the vertical rainfall field. Then, one earth station which can receive a 17 GHz signal from satellite is employed to detect the vertical rainfall field. However, the simulation results show that it is difficult for one earth station to achieve the inversion of rainfall field, and that the correlation coefficients between rainfall fields and inversed fields are 0.556, 0.504 and 0.364 respectively. Based on the result, two earth stations are jointly used. In this simulation, the result shows that after 500 iterations the correlation coefficients all increase above 0.98, and the average biases between rainfall field I, II, III and their inversed fields are 0.122, 0.159 and 0.537 mm/h, respectively. Meanwhile, the Euclidean distances decrease to 0.246, 0.235 and 0.812 mm/h, and the relative errors of entropy are both less than 2%. It can be seen from the inversion fields that the vertical distribution of rain rate is close to that of the real field, which suggests that the method proposed in this paper can basically achieve the inversion of vertical rainfall field by using earth-space links. In addition, with the combined detection of three earth stations the accuracy of the inversion results is significantly improved. The correlation coefficients are all close to 1 and the mean deviations are all on the order of 10–12 mm/h, indicating that the 2D vertical rainfall fields are accurately reconstructed. In the near future, the satellite constellation system will be globally deployed, which can promote the applications of our method in areas, such as plateaus, mountains and islands, where there exist no traditional observation data, serving as a supplement to existing precipitation measurements.
      通信作者: 刘西川, liuxc2012@hotmail.com
    • 基金项目: 国家自然科学基金(批准号: 41975030, 41505135)资助的课题
      Corresponding author: Liu Xi-Chuan, liuxc2012@hotmail.com
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 41975030, 41505135)
    [1]

    Giuli D, Toccafondi A, Gentili G B, Freni A 1991 J. Appl. Meteoro 30 1323Google Scholar

    [2]

    Messer H, Zinevich A, Alpert P 2006 Science 312 713Google Scholar

    [3]

    Rahimi A 2003 J Geophys Res 108 8

    [4]

    Rafael F, Roger H 2002 IEEE T Geosci Remote 40 760Google Scholar

    [5]

    Rahimi A, Holt A, Upton G, Krämer S, Redder A, Verworn H-R 2006 B Am Meteorol Soc 23 395

    [6]

    Fencl M, Rieckermann J, Schleiss M, Stransky D, Bares V 2013 Water Sci Technol 68 1810Google Scholar

    [7]

    Fencl M, Dohnal M, Rieckermann J, Bareš V 2017 Hydrol Earth Syst Sc 21 617Google Scholar

    [8]

    高太长 2013 气象水文装备 24 1

    Gao T C 2013 Meteol. Hydrol 24 1

    [9]

    高太长, 宋堃, 刘西川, 印敏, 刘磊, 姜世太 2015 物理学报 64 174301Google Scholar

    Gao T C, Song K, Liu X C, Yin M, Liu L, Jiang S T 2015 Acta Phys. Sin. 64 174301Google Scholar

    [10]

    姜世太, 高太长, 刘西川, 刘磊, 刘志田 2013 物理学报 62 154303Google Scholar

    Jiang S T, Gao T C, Liu X C, Liu L, Liu Z T 2013 Acta Phys. Sin. 62 154303Google Scholar

    [11]

    宋堃, 高太长, 刘西川, 印敏, 薛杨 2015 物理学报 64 244301Google Scholar

    Song K, Gao T C, Liu X C, Yin M, Xue Y 2015 Acta Phys. Sin. 64 244301Google Scholar

    [12]

    宋堃, 高太长, 刘西川, 印敏, 薛杨 2017 物理学报 66 154301Google Scholar

    Song K, Gao T C, Liu X C, Yin M, Xue Y 2017 Acta Phys. Sin. 66 154301Google Scholar

    [13]

    肖永伟, 孙晨华, 赵伟松 2018 国际太空 11 24Google Scholar

    Xiao Y W, Sun C H, Zhao W S 2018 Space Int. 11 24Google Scholar

    [14]

    李黄 2006 遥感学报 10 568Google Scholar

    Li H 2006 Journal of J Remote Sens. 10 568Google Scholar

    [15]

    Arslan C H, Aydin K, Urbina J, Dyrud L P 2014 IEEE Geoscience and Remote Sensing Symposium, Quebec City, July 13–18, 2014 p4111

    [16]

    François M, Laurent B, Cécile M 2015 J Atmos Ocean Tech 32 1709Google Scholar

    [17]

    Adirosi E, Facheris L 2018 2 nd URSI Atlantic Radio Science Meeting, Meloneras, 28 May–1 June, 2018 p1

    [18]

    International Telecommunication Union 2017 Rec. ITU-R p618−13

    [19]

    宗蓉 2013 博士学位论文 (南京: 南京信息工程大学)

    Zong R 2013 Ph.D. Dissertation (Nanjing: Nanjing University of Information Science and Technology) (in Chinese)

    [20]

    International Telecommunication Union 2005 Rec. ITU-R p838-3

    [21]

    杨富强, 张定华, 黄魁东, 王鹍, 徐哲 2014 物理学报 63 058701Google Scholar

    Yang F Q, Zhang D H, Huang K D, Wang K, Xu Z 2014 Acta Phys. Sin. 63 058701Google Scholar

    [22]

    Kaczmarz S 1993 Int J Control 57 1269Google Scholar

    [23]

    Hansen, Christian P, Hansen S, Maria 2012 J Comput Appl Math 236 2167Google Scholar

    [24]

    Andersen A H, Kak A C 1984 Ultrasonic Imaging 6 81Google Scholar

  • 图 1  星地链路几何结构

    Fig. 1.  Geometry of earth-space link.

    图 2  星地链路探测垂直降雨场

    Fig. 2.  The detection of vertical rainfall field with earth-space link.

    图 3  星地链路闪烁衰减、云致衰减、气体衰减、雨强为10 mm/h和20 mm/h的雨致衰减

    Fig. 3.  Scintillation attenuation, cloud attenuation, gas attenuation and rain attenuation with 10 mm/h and 20 mm/h rain rate in earth space link.

    图 4  计算星地链路降雨衰减简易图

    Fig. 4.  The sketch of calculation of earth-space link rain attenuation.

    图 5  幂律系数随仰角的变化

    Fig. 5.  The change of power law coefficients with elevation.

    图 6  反演算法流程图

    Fig. 6.  The flow chart of vertical rainfall field inversion method.

    图 7  基于MRR实测资料构建的垂直降雨场 (a) 降雨场I; (b) 降雨场II; (c) 降雨场III

    Fig. 7.  Vertical rainfall field derived from the data of MRR: (a) rainfall field I; (b) rainfall field II; (c) rainfall field III.

    图 8  降雨场网格化和星地链路探测示意图

    Fig. 8.  Sketch of grid for rainfall field and detection of earth-space links.

    图 9  垂直降雨场反演效果图 (a) 单地球站反演降雨场I; (b) 单地球站反演降雨场II; (c) 单地球站反演降雨场III; (d) 双地球站反演降雨场I; (e) 双地球站反演降雨场II; (f) 双地球站反演降雨场III; (g) 三地球站反演降雨场I; (h) 三地球站反演降雨场II; (i) 三地球站反演降雨场III

    Fig. 9.  Inversed vertical rainfall field: (a) Inversed field I with one earth station; (b) inversed field II with one earth station; (c) inversed field III with one earth station; (d) inversed field I with two earth stations; (e) inversed field II with two earth stations; (f) inversed field III with two earth stations; (g) inversed field I with three earth stations; (h) inversed field II with three earth stations; (i) inversed field III with three earth stations.

    图 10  欧氏距离随迭代次数的变化情况 (a) 单地球站反演降雨场I; (b) 单地球站反演降雨场II; (c) 单地球站反演降雨场III; (d) 双地球站反演降雨场I; (e) 双地球站反演降雨场II; (f) 双地球站反演降雨场III;(g) 三地球站反演降雨场I; (h) 三地球站反演降雨场II; (i) 三地球站反演降雨场III

    Fig. 10.  The change of Euclidean distance with iterations: (a) Inversed field I with one earth station; (b) inversed field II with one earth station; (c) inversed field III with one earth station;(d) inversed field I with two earth stations; (e) inversed field II with two earth stations; (f) inversed field III with two earth stations;(g) inversed field I with three earth stations; (h) inversed field II with three earth stations; (i) inversed field III with three earth stations.

    图 11  欧氏距离随迭代次数的变化情况 (a) 单地球站反演降雨场I; (b) 单地球站反演降雨场II; (c) 单地球站反演降雨场III; (d) 双地球站反演降雨场I; (e) 双地球站反演降雨场II; (f) 双地球站反演降雨场III; (g) 三地球站反演降雨场I; (h) 三地球站反演降雨场II; (i) 三地球站反演降雨场III

    Fig. 11.  The change of entropy with iterations: (a) inversed field I with one earth station; (b) inversed field II with one earth station; (c) inversed field III with one earth station; (d) inversed field I with two earth stations; (e) inversed field II with two earth stations; (f) inversed field III with two earth stations; (g) inversed field I with three earth stations; (h) inversed field II with three earth stations; (i) inversed field III with three earth stations.

    表 1  星地链路参数

    Table 1.  Parameters of earth-space link.

    地球站频率/ GHz极化方式Δθ/(°)θMin/(°)C/dB
    117垂直0.10.091105
    217垂直0.10.065105
    317垂直0.11.00105
    下载: 导出CSV
  • [1]

    Giuli D, Toccafondi A, Gentili G B, Freni A 1991 J. Appl. Meteoro 30 1323Google Scholar

    [2]

    Messer H, Zinevich A, Alpert P 2006 Science 312 713Google Scholar

    [3]

    Rahimi A 2003 J Geophys Res 108 8

    [4]

    Rafael F, Roger H 2002 IEEE T Geosci Remote 40 760Google Scholar

    [5]

    Rahimi A, Holt A, Upton G, Krämer S, Redder A, Verworn H-R 2006 B Am Meteorol Soc 23 395

    [6]

    Fencl M, Rieckermann J, Schleiss M, Stransky D, Bares V 2013 Water Sci Technol 68 1810Google Scholar

    [7]

    Fencl M, Dohnal M, Rieckermann J, Bareš V 2017 Hydrol Earth Syst Sc 21 617Google Scholar

    [8]

    高太长 2013 气象水文装备 24 1

    Gao T C 2013 Meteol. Hydrol 24 1

    [9]

    高太长, 宋堃, 刘西川, 印敏, 刘磊, 姜世太 2015 物理学报 64 174301Google Scholar

    Gao T C, Song K, Liu X C, Yin M, Liu L, Jiang S T 2015 Acta Phys. Sin. 64 174301Google Scholar

    [10]

    姜世太, 高太长, 刘西川, 刘磊, 刘志田 2013 物理学报 62 154303Google Scholar

    Jiang S T, Gao T C, Liu X C, Liu L, Liu Z T 2013 Acta Phys. Sin. 62 154303Google Scholar

    [11]

    宋堃, 高太长, 刘西川, 印敏, 薛杨 2015 物理学报 64 244301Google Scholar

    Song K, Gao T C, Liu X C, Yin M, Xue Y 2015 Acta Phys. Sin. 64 244301Google Scholar

    [12]

    宋堃, 高太长, 刘西川, 印敏, 薛杨 2017 物理学报 66 154301Google Scholar

    Song K, Gao T C, Liu X C, Yin M, Xue Y 2017 Acta Phys. Sin. 66 154301Google Scholar

    [13]

    肖永伟, 孙晨华, 赵伟松 2018 国际太空 11 24Google Scholar

    Xiao Y W, Sun C H, Zhao W S 2018 Space Int. 11 24Google Scholar

    [14]

    李黄 2006 遥感学报 10 568Google Scholar

    Li H 2006 Journal of J Remote Sens. 10 568Google Scholar

    [15]

    Arslan C H, Aydin K, Urbina J, Dyrud L P 2014 IEEE Geoscience and Remote Sensing Symposium, Quebec City, July 13–18, 2014 p4111

    [16]

    François M, Laurent B, Cécile M 2015 J Atmos Ocean Tech 32 1709Google Scholar

    [17]

    Adirosi E, Facheris L 2018 2 nd URSI Atlantic Radio Science Meeting, Meloneras, 28 May–1 June, 2018 p1

    [18]

    International Telecommunication Union 2017 Rec. ITU-R p618−13

    [19]

    宗蓉 2013 博士学位论文 (南京: 南京信息工程大学)

    Zong R 2013 Ph.D. Dissertation (Nanjing: Nanjing University of Information Science and Technology) (in Chinese)

    [20]

    International Telecommunication Union 2005 Rec. ITU-R p838-3

    [21]

    杨富强, 张定华, 黄魁东, 王鹍, 徐哲 2014 物理学报 63 058701Google Scholar

    Yang F Q, Zhang D H, Huang K D, Wang K, Xu Z 2014 Acta Phys. Sin. 63 058701Google Scholar

    [22]

    Kaczmarz S 1993 Int J Control 57 1269Google Scholar

    [23]

    Hansen, Christian P, Hansen S, Maria 2012 J Comput Appl Math 236 2167Google Scholar

    [24]

    Andersen A H, Kak A C 1984 Ultrasonic Imaging 6 81Google Scholar

计量
  • 文章访问数:  6361
  • PDF下载量:  38
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-08-15
  • 修回日期:  2019-11-02
  • 上网日期:  2020-01-01
  • 刊出日期:  2020-01-20

/

返回文章
返回