留言板

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

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

中国遥感卫星辐射校正场敦煌戈壁场区光环境变化研究

李元,张勇,胡丽琴,陆其峰,卢乃锰

downloadPDF
李元, 张勇, 胡丽琴, 陆其峰, 卢乃锰. 中国遥感卫星辐射校正场敦煌戈壁场区光环境变化研究[J]. , 2021, 14(5): 1231-1242. doi: 10.37188/CO.2020-0129
引用本文: 李元, 张勇, 胡丽琴, 陆其峰, 卢乃锰. 中国遥感卫星辐射校正场敦煌戈壁场区光环境变化研究[J]. , 2021, 14(5): 1231-1242.doi:10.37188/CO.2020-0129
LI Yuan, ZHANG Yong, HU Li-qin, LU Qi-feng, LU Nai-meng. Investigation of optical environment changes in the Dunhuang gobi site of the Chinese radiometric calibration sites[J]. Chinese Optics, 2021, 14(5): 1231-1242. doi: 10.37188/CO.2020-0129
Citation: LI Yuan, ZHANG Yong, HU Li-qin, LU Qi-feng, LU Nai-meng. Investigation of optical environment changes in the Dunhuang gobi site of the Chinese radiometric calibration sites[J].Chinese Optics, 2021, 14(5): 1231-1242.doi:10.37188/CO.2020-0129

中国遥感卫星辐射校正场敦煌戈壁场区光环境变化研究

doi:10.37188/CO.2020-0129
基金项目:国家重点研发计划课题(No. 2018YFB0504601);国家自然基金委员会面上项目(No. 41271373);民用航天“十三五”技术预先研究项目(No. D040401);上海航天科技创新基金资助项目(No. SAST2019-044)
详细信息
    作者简介:

    李 元(1978—),女,北京人,工学博士,国家卫星气象中心副研究员,2001年,2006年于北京理工大学分别获学士,博士学位,主要从事在轨遥感器定标与产品反演等方面的研究。E-mail:liyuan@cma.gov.cn

    张 勇(1977—),男,陕西镇巴人,理学博士,国家卫星气象中心研究员,2000年,2003年于西北大学分别获学士,硕士学位,2006年于中国科学院遥感应用研究所获博士学位,主要从事卫星遥感器辐射定标与真实性检验方面的研究。E-mail:zhangyong@cma.gov.cn

  • 中图分类号:O432.1; P422.1

Investigation of optical environment changes in the Dunhuang gobi site of the Chinese radiometric calibration sites

Funds:Supported by National Key Research and Development Project (No. 2018YFB0504601); National Natural Science Foundation of China General Project (No. 41271373); Civil Aerospace “Thirteenth Five-Year” Technology Pre-research Project (No. D040401); Shanghai Aerospace Science and Technology Innovation Fund Project (No. SAST2019-044)
More Information
  • 摘要:为有效评估集热塔散射辐射对敦煌场区光环境的影响程度,本文采用Monte Carlo三维辐射传输模型模拟与CE318多通道光度计等高线实测分析相结合的定量分析方法,以解决散射辐射交融于背景辐射中难以定量评估的问题。通过使用新型的ASC200云量自动观测仪,提高晴空辩识精度。通过开发CE318四象限定位修正算法,有效提高观测数据质量。2020年1~3月收集到的有效数据显示除了550 nm通道,集热塔未对天空漫射辐射产生明显影响。对于500 nm通道,在有效数据对应的观测几何下(距离0.87~3.07 km,观测天顶角为77.30°~51.32°),集热塔吸热器对天空漫射辐射的影响不超过0.93%。与模型模拟结果相结合进行分析,得出如下结论:当距离电站2 km时大电站散射辐射带来的天空漫射辐射相对变化<1.62%;当与电站距离≥3 km时相对变化<0.93%。本项研究成果对利用敦煌场开展遥感定量化应用、准确评估发电站引进的不确定度因素具有积极意义。

  • 图 1在自动观测基地附近观察到的工作状态的集热塔

    Figure 1.Heat collection tower in working condition observed from the automatic observation base

    图 2等高线测量角度示意图

    Figure 2.Schematic diagram of the almucantar measurement angle

    图 3卫星观测到的集热塔反射光的辐射能量水平分布示例

    Figure 3.An example of the radiant energy level distribution of the light reflected by the heat collection tower observed by the satellite

    图 42020年3月2日UTC 7:00:05时刻ASC200云量自动观测仪红外云图

    Figure 4.Infrared cloud image obtained from ASC200 cloud cover automatic observer at 7:00:05 UTC on March 2, 2020

    图 5归一化的全天空漫射辐射光谱与CE318不同通道光谱响应函数

    Figure 5.Normalized full sky diffuse irradiance and CE318 spectral response function

    图 6440 nm通道RdevBr为最小值(−2.07%)时对应时刻(03−09 09:59)的辐亮度(a)与相对偏差(b)

    Figure 6.Radiance (a) and relative deviation (b) at the corresponding time (03−09 09:59) when the 440nm channel RdevBris the minimum (−2.07%)

    图 7500 nm通道RdevBr为最小值(−3.65%)时对应时刻(03−03 08:03)的辐亮度(a)与相对偏差(b)

    Figure 7.Radiance (a) and relative deviation (b) at the corresponding time (03−03 08:03) when the 500 nm channel RdevBris the minimum (−3.65%)

    图 8500 nm通道RdevBr为最大值(0.93%)时对应时刻(01−25 09:04)的辐亮度(a)与相对偏差(b)

    Figure 8.Radiance (a) and relative deviation (b) at the corresponding time (01−25 09:04) when the 500 nm channel RdevBr is the maximum (0.93%)

    表 1模拟的相关参数定义

    Table 1.Definitions of parameters in the simulation

    相关参数 取值范围
    模拟区域范围 11 km×11 km
    (设定集热塔位于模拟区域中心像元的中央)
    水平空间分辨率 1 km
    反射光源性质 离地面高度260 m;向上各角度随机发射
    计算波段 0.55 μm
    大气背景 中纬夏季大气考虑H2O、O2、O3气体吸收和
    大气分子瑞利散射
    气溶胶特性 光学厚度0.1;单次散射反照率0.938;
    不对称因子0.764
    地表特性 地表反照率0.2
    下载: 导出CSV

    表 2CE318通道设置

    Table 2.CE318 channel specifications

    编号 中心波长/nm 带宽/nm 观测模式
    1 1020 10 S、A、K
    2 1640 25 S、A、K
    3 870 10 S、A、K
    4 675 10 S、A、K
    5 440 10 S、A、K
    6 500 10 S、A、K
    7 1020i 10 S
    8 936 10 S
    9 380 2 S
    10 340 2 S
    下载: 导出CSV

    表 3CE318性能参数

    Table 3.CE318 performance parameters

    项目 内容
    视场 1.2°
    探测器 铟镓砷(InGaSn)探测器:1020i nm,
    1640 nm;硅探测器:其他
    工作温度 −30~+ 60 °C
    太阳追踪方法 四象限探测器主动跟踪
    跟踪精度 优于0.1°
    天空漫射辐射
    观测频率
    固定时刻与大气质量数
    动态范围 增益可调
    量化等级 15
    下载: 导出CSV

    表 4ASC200性能参数

    Table 4.ASC200 performance parameters

    项目 内容
    观测指标 可见光云量、红外云量、综合云量
    红外波长 8~14 μm
    视场角 可见光180°,红外160°
    采样周期 10 min
    工作温度 −45~55 °C
    下载: 导出CSV

    表 5HIM性能参数

    Table 5.HIM performance parameters

    项目 内容
    观测内容 总天空辐照度,天空漫射辐照度,漫射/总漫射比
    波长 400~2400 nm
    半高宽 4 nm @ 400~950 nm,15 nm @ 950~1700 nm,
    20 nm @ 1700~2400 nm
    采样期 6 min
    工作温度 −30~60 °C
    信噪比 ≥600@ 400~1000 nm,≥300@ 1000~1700 nm,
    ≥200@ 1700~2400 nm(M= 2,非吸收通道)
    下载: 导出CSV

    表 6观测数据的质量控制阈值

    Table 6.Quality control threshold of observation data

    项目 数值
    云量 0
    光学厚度 <0.2
    ALL与ALR相对偏差的标准差 <1%
    相对偏差计算点与大小电站的投影距离 >5 km
    太阳与大小电站的投影距离 >2.5 km
    ALL与ALR观测时间差异 <1 min
    下载: 导出CSV

    表 7卫星观测到的集热塔反射光的辐射能量相对分布情况

    Table 7.The relative distribution of the radiant energy of the light reflected by the heat collection tower observed by the satellite

    距离/km 10° 20° 30° 40° 50°
    0 1.000 0.780 0.549 0.381 0.277 0.202
    1 0.434 0.369 0.285 0.194 0.179 0.162
    2 0.035 0.033 0.033 0.032 0.037 0.047
    3 0.020 0.019 0.020 0.021 0.023 0.027
    4 0.016 0.016 0.016 0.017 0.020 0.023
    5 0.015 0.015 0.015 0.016 0.018 0.021
    下载: 导出CSV

    表 8各级筛选后有效观测数据量值统计(组)

    Table 8.Statistics of effective observation data after being selected at different levels (group)

    筛选项目 有效观测数据值 合计
    1020 nm 1640 nm 870 nm 675 nm 440 nm 500 nm
    (a) 804 778 611 466 787 565 4011
    (b) 135 132 81 54 135 70 607
    (c) 45 77 47 44 43 44 300
    (d) 0 3 1 6 27 22 59
    (e) 0 2 1 2 22 16 43
    (f) 0 2 0 1 14 9 26
    (g) 0 2 0 1 14 9 26
    (h) 0 2 0 1 14 9 26
    下载: 导出CSV

    表 9吸热器散射辐射带来的相对偏差(440 nm)

    Table 9.Relative deviation due to scattered radiation from the heat absorber (440 nm)

    UTC Time AOD Sz/(°) Sa/(°) Beta/(°) HB/km on DegB/(°) DdegB/(°) RdevB/% RdevBr/% HS/km on DegS/(°) DdegS/(°) RdevS/% RdevSr/%
    03−10 08:59 0.15 61.83 235.34 0.07 2.06 1 85.37 4.63 −1.18 −1.08 1.54 0 124.87 4.87 −0.29 −0.31
    03−10 07:59 0.14 53.28 220.75 0.11 2.86 1 70.77 0.77 −1.14 −0.92 2.14 0 110.27 9.73 −1.57 −1.49
    03−10 06:59 0.15 47.20 202.90 0.13 3.55 1 52.92 2.92 −0.04 0.31 2.65 0 92.42 2.42 −1.37 −1.18
    03−09 10:27 0.17 77.14 252.07 0.09 0.88 1 102.09 2.09 −0.70 −0.72 0.66 0 141.59 1.59 0.24 0.14
    03−09 09:59 0.15 72.12 247.04 0.12 1.24 1 97.06 2.94 −2.11 −2.07 0.93 0 136.56 3.44 −0.63 −0.74
    03−02 09:02 0.15 64.41 233.46 0.06 1.84 1 83.48 3.48 −0.51 −0.42 1.38 1 122.98 2.98 −1.23 −1.24
    03−02 07:02 0.18 50.22 202.09 0.13 3.19 1 52.12 2.12 −0.20 0.19 2.39 0 91.62 1.62 −1.47 −1.28
    02−27 10:14 0.14 77.10 246.64 0.09 0.88 1 96.66 3.34 −0.89 −0.87 0.66 0 136.16 3.84 −0.75 −0.84
    02−27 09:02 0.16 65.23 232.61 0.16 1.77 1 82.63 2.63 −1.47 −1.23 1.33 0 122.13 2.13 0.47 0.39
    02−27 08:02 0.19 57.00 218.41 0.09 2.49 1 68.43 1.57 0.09 0.30 1.86 0 107.93 7.93 1.33 1.38
    02−27 07:02 0.18 51.24 201.49 0.08 3.08 1 51.51 1.51 −0.51 −0.29 2.30 0 91.02 1.02 −0.61 −0.50
    01−25 09:24 0.16 76.61 230.44 0.11 0.92 1 80.46 0.46 −0.18 0.00 0.69 0 119.97 0.03 −0.81 −0.91
    01−25 09:03 0.17 73.56 226.29 0.06 1.13 1 76.32 3.68 −0.14 −0.03 0.85 0 115.82 4.18 −0.60 −0.64
    01−10 07:56 0.19 68.63 211.65 0.04 1.50 0 61.67 1.67 −1.02 −0.89 1.12 1 101.17 1.17 −0.44 −0.43
    下载: 导出CSV

    表 10吸热器散射辐射带来的相对偏差(500 nm)

    Table 10.Relative deviation due to scattered radiation from the heat absorber (500 nm)

    UTC Time AOD Sz/(°) Sa/(°) Beta/(°) HB/km on DegB/(°) DdegB/(°) RdevB/% RdevBr/% HS/km on DegS/(°) DdegS/(°) RdevS/% RdevSr/%
    03−10 09:00 0.14 62.01 235.59 0.12 2.04 1 85.61 4.39 −0.85 −0.63 1.53 0 125.11 5.11 0.88 0.87
    03−10 08:00 0.14 53.42 221.04 0.25 2.85 1 71.06 1.06 −0.74 −0.14 2.13 0 110.56 9.44 −1.19 −0.97
    03−03 08:03 0.15 55.92 219.77 0.14 2.60 1 69.80 0.20 −4.01 −3.65 1.94 0 109.30 9.30 −1.21 −1.10
    02−27 10:16 0.14 77.30 246.84 0.08 0.87 1 96.86 3.14 0.29 0.32 0.65 0 136.36 3.64 −0.10 −0.18
    02−27 07:03 0.16 51.32 201.84 0.14 3.07 1 51.86 1.86 −2.24 −1.81 2.29 0 91.36 1.36 −0.46 −0.22
    01−25 09:25 0.15 76.80 230.69 0.17 0.90 1 80.71 0.71 0.57 0.88 0.68 0 120.21 0.21 −0.09 −0.21
    01−25 09:04 0.16 73.72 226.51 0.11 1.12 1 76.54 3.46 0.68 0.93 0.84 0 116.04 3.96 −0.20 −0.25
    01−10 08:58 0.17 75.87 224.52 0.08 0.97 0 74.54 4.54 0.49 0.71 0.72 1 114.04 5.96 0.31 0.28
    01−10 07:57 0.17 68.75 211.92 0.07 1.49 0 61.94 1.94 −0.50 −0.25 1.12 1 101.44 1.44 −0.02 0.01
    下载: 导出CSV

    表 11大电站散射辐射带来的天空漫射辐射相对变化随距离与观测角度变化(%)

    Table 11.The relative change of the sky diffuse radiation caused by the scattered radiation from the large power station changes with distance and the observation angle (%)

    距离/km 天空漫射辐射相对变化(%)
    10° 20° 30° 40° 50°
    0 34.44 26.87 18.91 13.12 9.54 6.96
    1 14.95 12.71 9.82 6.68 6.17 5.58
    2 1.21 1.14 1.14 1.10 1.27 1.62
    3 0.69 0.65 0.69 0.72 0.79 0.93
    4 0.55 0.55 0.55 0.59 0.69 0.79
    5 0.52 0.52 0.52 0.55 0.62 0.72
    下载: 导出CSV
  • [1] 张勇, 祁广利, 戎志国. 卫星红外遥感器辐射定标模型与方法[M]. 北京: 科学出版社, 2015.

    ZHANG Y, QI G L, RONG ZH G.Models and Methodologies of Radiometric Calibrations for Satellite Infrared Remote Sensors[M]. Beijing: Science Press, 2015. (in Chinese)
    [2] 卢乃锰, 谷松岩. 气象卫星发展回顾与展望[J]. 遥感学报,2016,20(5):832-841.

    LU N M, GU S Y. Review and prospect on the development of meteorological satellites[J].Journal of Remote Sensing, 2016, 20(5): 832-841. (in Chinese)
    [3] ZAHNG Y, LI ZH L, LI J. Comparisons of emissivity observations from satellites and the ground at the CRCS Dunhuang Gobi site[J].Journal of Geophysical Research:Atmospheres, 2014, 119(22): 13026-13041.doi:10.1002/2014JD022216
    [4] HU X Q, LIU J J, SUN L,et al. Characterization of CRCS Dunhuang test site and vicarious calibration utilization for Fengyun (FY) series sensors[J].Canadian Journal of Remote Sensing, 2010, 36(5): 566-582.doi:10.5589/m10-087
    [5] 张勇, 李元, 戎志国, 等. 中国遥感卫星辐射校正场陆表热红外发射率光谱野外测量[J]. 光谱学与光谱分析,2009,29(5):1213-1217.doi:10.3964/j.issn.1000-0593(2009)05-1213-05

    ZHANG Y, LI Y, RONG ZH G,et al. Field measurement of gobi surface emissivity spectrum at Dunhuang calibration site of China[J].Spectroscopy and Spectral Analysis, 2009, 29(5): 1213-1217. (in Chinese)doi:10.3964/j.issn.1000-0593(2009)05-1213-05
    [6] 张勇, 戎志国, 闵敏. 中国遥感卫星辐射校正场热红外通道在轨场地辐射定标方法精度评估[J]. 地球科学进展,2016,31(2):171-179.doi:10.11867/j.issn.1001-8166.2016.02.0171

    ZHANG Y, RONG ZH G, MIN M. Accuracy evaluations of the CRCS In-orbit field radiometric calibration methods for thermal infrared channels[J].Advances in Earth Science, 2016, 31(2): 171-179. (in Chinese)doi:10.11867/j.issn.1001-8166.2016.02.0171
    [7] ZHANG Y X, ZHANG G SH, LIU ZH Q,et al. Spectral reflectance measurements at the China radiometric calibration test site for the remote sensing satellite sensor[J].Acta Meteorologica Sinica, 2001, 15(3): 377-382.
    [8] 李元, 戎志国, 郑照军, 等. FY-3A扫描辐射计的可见近红外通道在轨场地定标[J]. 光学 精密工程,2009,17(12):2966-2974.

    LI Y, RONG ZH G, ZHENG ZH J,et al. Post launch site calibration of visible and near-infrared channels of FY-3A visible and infrared radiometers[J].Optics and Precision Engineering, 2009, 17(12): 2966-2974. (in Chinese)
    [9] 李元, 张勇, 刘京晶, 等. 风云二号静止气象卫星可见光通道辐射校正场定标方法研究[J]. 光学学报,2009,29(1):41-46.doi:10.3788/AOS20092901.0041

    LI Y, ZHANG Y, LIU J J,et al. Calibration of the visible and near-infrared channels of the FY-2C/FY-2D GEO meteorological satellite at radiometric site[J].Acta Optica Sinica, 2009, 29(1): 41-46. (in Chinese)doi:10.3788/AOS20092901.0041
    [10] 孙凌, 郭茂华, 徐娜, 等. 基于敦煌场地定标的FY-3 MERSI反射太阳波段在轨响应变化分析[J]. 光谱学与光谱分析,2012,32(7):1869-1877.doi:10.3964/j.issn.1000-0593(2012)07-1869-09

    SUN L, GUO M H, XU N,et al. On-orbit response variation analysis of FY-3 MERSI reflective solar bands based on Dunhuang site calibration[J].Spectroscopy and Spectral Analysis, 2012, 32(7): 1869-1877. (in Chinese)doi:10.3964/j.issn.1000-0593(2012)07-1869-09
    [11] HAN Q J, ZHANG X W, LIU L,et al. Lifetime radiometric calibration of HJ-1A/B CCD sensor using Dunhuang Gobi site[J].Proceedings of SPIE, 2014, 9299: 92990B.
    [12] ZHANG X W, HAN Q J, LIU L. Using the Dunhuang test site to monitor the radiometric stability of the ZY-3 multispectral sensor[J].Proceedings of SPIE, 2014, 9298: 92980H.
    [13] CHEN ZH CH, ZHANG B, ZHANG H,et al. Vicarious calibration of Beijing-1 multispectral imagers[J].Remote Sensing, 2014, 6(2): 1432-1450.doi:10.3390/rs6021432
    [14] LI Y, RONG ZH G, LI Y Q,et al. Solar reflection band site automatic calibration by the Dunhuang site automatic observation radiometric calibration operational system[J].Proceedings of SPIE, 2019, 11127: 111271M.
    [15] LI Y, RONG Z G, LI Y Q,et al.. Introduction of Dunhuang site automatic observation radiometric calibration operational system (DARCOS)[C].Proceedings of CEOS-WGCV-IVOS Meeting 31, 2019.
    [16] 吕佳彦, 何明元, 陈林, 等. 基于敦煌辐射校正场的自动化辐射定标方法[J]. 光学学报,2017,37(8):0801003.doi:10.3788/AOS201737.0801003

    LÜ J Y, HE M Y, CHEN L,et al. Automated radiation calibration method based on Dunhuang radiometric calibration site[J].Acta Optica Sinica, 2017, 37(8): 0801003. (in Chinese)doi:10.3788/AOS201737.0801003
    [17] 邱刚刚, 李新, 韦玮, 等. 基于场地自动化观测技术的遥感器在轨辐射定标试验与分析[J]. 光学学报,2016,36(7):0701001.doi:10.3788/AOS201636.0701001

    QIU G G, LI X, WEI W,et al. Experiment and analysis of on-orbit radiometric calibration for remote sensors based on in-site automated observation technology[J].Acta Optica Sinica, 2016, 36(7): 0701001. (in Chinese)doi:10.3788/AOS201636.0701001
    [18] ZHANG Y, LI X, RONG ZH G,et al. China radiometric calibration sites ground-based automatic observing systems for CAL/VA[J].Proceedings of SPIE, 2015, 9639: 96391E.doi:10.1117/12.2194633
    [19] 敦煌市委宣传部. 敦煌首航节能10兆瓦塔式熔盐光热发电项目有序推进[EB/OL]. (2015-05-04).http://www.cspplaza.com/article-5110-1.html.

    The Publicity Department of Dunhuang municipal Party committee. The energy-saving 10 MW tower molten salt solar thermal power generation project of Dunhuang Capital Airlines was promoted in an orderly manner[EB/OL]. (2015-05-04).http://www. cspplaza.com/article-5110-1.html.
    [20] 环球网. 敦煌100 MW塔式熔盐光热电站土建完成进入安装环节[EB/OL]. (2018-05-18).https://baijiahao.baidu.com/s?id=1600755797116589317&wfr=spider&for=pc.

    Global network. The civil engineering of Dunhuang 100MW tower molten salt optical thermal power station has been completed and entered the installation link[EB/OL]. (2018-05-18).https://baijiahao.baidu. com/s?id=1600755797116589317&wfr=spider&for=pc.
    [21] 胡丽琴. 风云四号卫星闪电辐射观测的辐射传输模拟研究[C]. 2011年第二十八届中国气象学会年会论文集, 中国气象学会, 2011: 1-10.

    HU L Q. Radiative transfer simulation of the Fengyun-4 lightning radiation observation[C].Proceedings of the 28th Annual Meeting of the Chinese Meteorological Society, Chinese Meteorological Society, 2011: 1-10. (in Chinese)
    [22] YAO W X, LI ZH R, ZHAO Q,et al. A new anisotropic diffuse radiation model[J].Energy Conversion and Management, 2015, 95: 304-313.doi:10.1016/j.enconman.2015.01.016
    [23] Cimel. Multiband photometer CE318-T User’s Manual (rev. January 2018) [EB/OL]. (2021-8-16).https://aeronet.gsfc.nasa.gov/new_web/Documents/CE318_T_Photometer_UserManual_V4.3.pdf.
    [24] FA T, XIE W Y, WANG Y R,et al. Development of an all-sky imaging system for cloud cover assessment[J].Applied Optics, 2019, 58(20): 5516-5524.doi:10.1364/AO.58.005516
    [25] ESTELLES V, CAMPANELLI M, UTRILLAS M P,et al. Comparison of AERONET and SKYRAD4.2 inversion products retrieved from a Cimel CE318 sunphotometer[J].Atmospheric Measurement Techniques, 2012, 5(3): 569-579.doi:10.5194/amt-5-569-2012
    [26] BARRETO Á, CUEVAS E, GRANADOS-MUÑOZ M J,et al. The new sun-sky-lunar Cimel CE318-T multiband photometer – A comprehensive performance evaluation[J].Atmospheric Measurement Techniques, 2016, 9(2): 631-654.doi:10.5194/amt-9-631-2016
  • 加载中
图(8)/ 表(11)
计量
  • 文章访问数:1421
  • HTML全文浏览量:446
  • PDF下载量:87
  • 被引次数:0
出版历程
  • 收稿日期:2021-06-25
  • 修回日期:2021-07-21
  • 网络出版日期:2021-08-21
  • 刊出日期:2021-09-18

目录

    /

      返回文章
      返回
        Baidu
        map