Research on the Nowcasting Method of Solar Energy Resources Base on FY-4A Ground Radiation Products

  • Dong YE ,
  • Jinqiu LIANG ,
  • Yanbo SHEN ,
  • Chuanhui WANG ,
  • Rui CHANG
Expand
  • 1. Public Meteorological Service Center of China Meteorological Administration,Beijing 100081,China
    2. Public Meteorological Service Center of Anhui,Hefei 230031,Anhui,China
    3. Datong Meteorological Bureau of Shanxi Province,Datong 037010,Shanxi,China
    4. National Climate Center,Beijing 100081,China

Received date: 2021-07-05

  Revised date: 2021-11-30

  Online published: 2022-09-08

Cite this article

Dong YE , Jinqiu LIANG , Yanbo SHEN , Chuanhui WANG , Rui CHANG . Research on the Nowcasting Method of Solar Energy Resources Base on FY-4A Ground Radiation Products[J]. Plateau Meteorology, 2022 , 41(4) : 1030 -1040 . DOI: 10.7522/j.issn.1000-0534.2021.00110

1 引言

在对可持续能源的需求日益增长的背景下, 太阳能是满足这一新需求的最有希望的解决方案之一(石广玉, 2007申彦波等, 2013Ayet and Tandeo, 2018)。但由于太阳辐射具有间歇性、 波动性和周期性的特点, 导致太阳能发电功率输出的不稳定性, 对电力系统的安全运行、 调度和控制造成日益严重的影响(王长贵和王斯成, 2009丁宇宇等, 2014Wan et al, 2014)。因此, 准确预测地表太阳辐射可以为电力调度提供重要的决策支持, 有效降低单纯依靠光伏资源的电力系统运行成本, 产生更大的经济效益和社会效益(Luis et al, 2019Sergio et al, 2020)。随着光伏发电产业的迅速发展, 装机规模越来越大, 电网公司明确要求太阳能发电企业要提供太阳辐射和发电功率的短期(3~7天)和短临(0~1 h或4 h)预报, 为电网调度提供依据(中国电力企业联合会, 2013)。非常短期的地表太阳辐照度预测, 也称为“短临预报”(0~4 h前)可以提供高频范围内光伏功率变化的信息, 主要用于微电网的暂态分析和控制策略研究(Yang et al, 2020)。本文将重点关注辐射短临预报问题。
当前国内外太阳能资源短临预报方法根据使用资料的不同, 主要有三类, 一是基于地面辐射观测数据直接外推预报, 这类方法多采用统计法, 根据现在的天气状况建立回归方程统计后延关系从而获得水平面总辐射(程序等, 2011江滢等, 2013崔杨等, 2018); 二是基于地基云图建模反演预报, 这类方法基于全天空成像仪的云图观测, 通过对云图的识别处理以及云团运动预测, 建立云和辐射的关系, 最终获得水平面总辐射预报(丁宇宇等, 2014高晓清等, 2018蒋俊霞等, 2020)。这两类方法都必须依赖于现场观测, 并且长期观测后, 仪器误差可能导致不准确; 另外虽然地基云图的时空分辨率很高, 但观测全天空成像仪的范围是有限的(例如, TSI-880的有效检测范围是5 km), 对于提前3 h的预报, 云团可能已经移出了成像仪的检测范围。因此全天空成像仪适合0~15 min或者0~1 h的预报。由于全天空成像仪检测的局限性, 需要考虑大范围的0~3 h预报。第三类方法是直接基于卫星云图的辐射短临预报方法(Miller et al, 2018Ayet and Tandeo, 2018), 卫星云图具有较大的空间尺度、 较高的空间分辨率的特点, 可以实现大范围的云状态观测和辐射预报信息, 但是从卫星云图到辐射的预报有赖于云到辐射的模型, 尽管云是影响辐射的主要因子, 但是影响辐射不仅只有云的影响, 还有气溶胶, 水汽等的影响(沈钟平和张华, 2009尹青等, 2011Zhang et al, 2013张凤霞等, 2016张华等, 20172019邱宗旭等, 2019)。
FY-4A是我国于2016年12月11日发射的世界上先进的新一代静止卫星(张志清等, 2017), 在沙尘监测、 云分析等方面获得了较好的应用(崔林丽等, 2020黄守友和徐国强, 2020姜红等, 2021)。其反演的地表太阳入射辐射是FY-4A的L2定量反演产品之一, 考虑了云、 气溶胶、 水汽含量、 地表反照率、 地表海拔高度等影响参数。该产品输出的空间分辨率为4 km, 时间分辨率可到15 min(一般为1 h, 另外每3 h除整点外, 在整点前后15 min各有一次观测)(梁静秋等, 2021)。
山西省是我国重要的综合能源基地, 光伏发电发展迅速, 据国家能源局数据统计, 截至2019年底, 山西省光伏发电累计装机容量为1.088×107 kW, 占全国的5.3%。本文基于国家级辐射观测站山西太原气象站的分钟辐射观测数据和台站所在位置的基于FY-4A地面辐射产品研究建立太阳能资源短临预报方法。

2 资料来源与方法介绍

2.1 数据资料

2.1.1 地面气象站总辐射观测数据

太原气象台(37.6167°N, 112.5833°E, 海拔776.3 m)地处山西中部, 位于太原市区南部, 属于城乡接合部, 2013年1月1日迁站到此, 环境评分为84.2分, 该站属暖温带亚湿润区(吴林栋, 2020)。根据统计, 近30年(1989 -2018年)该站的年平均日照时数为2460.1 h, 年平均降水量为431.2 mm。太原气象站为国家级辐射观测站, 目前采用TBQ-2型总辐射表观测总辐照度, 其获得的太阳总辐照度数据为逐分钟观测数据(梁进秋等, 2020)。本文所用辐射数据为该站2019年全年分钟时刻的1 min平均总辐照度, 数据来源于山西省气象信息中心。文中所用时间均为北京时。

2.1.2  FY-4A卫星总辐射数据

本文所用辐射数据为2019年1 -12月的太原气象站对应的FY-4A 像素点的逐时总辐射数据, 从FY-4A圆盘图中匹配距离最近的像素点(37.6248°N, 112.5637°E), 提取像素点2019年1 -12月整点时刻瞬时总辐照度, 数据来源于国家卫星气象中心。FY-4A地表太阳入射辐射产品反演算法采用离散纵标法辐射传输软件包, 受地球曲率的影响, 当太阳天顶角大于70°时, 当前所采用的辐射传输软件包中的平面平行模式已不适用(吴北婴, 1999), 因此, 在反演过程中, 将太阳天顶角的边界值设为70°, 超过该边界值时, 采用地外太阳辐射乘以晴空指数作为反演总辐照度输出。经过综合统计, 太原站2019年同时满足有卫星辐射、 地面观测总辐射的小时数为4377个时次。总体而言, 太原站年平均整点总辐照度FY-4A比气象站大, 相对误差在19.1%~73.5%(梁进秋等, 2020)。

2.2 太阳能资源短临预报方法

2.2.1 总体思路

基于FY-4A地面辐射产品的太阳能资源短临预报方法由多个环节组成, 包括数据插补、 短临预报、 滚动订正, 具体流程见图1。卫星地面辐射产品输入为逐小时数据, 首先将小时数据插补成15 min数据, 然后基于均生函数法利用过去2 h数据预报未来4 h数据, 时间分辨率为15 min, 再对预报数据进行订正, 最后获得订正后的预报总辐射。该工作逐15 min滚动预报和订正一次。最后对结果进行统计检验。
图1 基于FY-4A地面辐射产品的太阳能资源短临预报方法流程图

Fig.1 Flow chart of nowcasting method of solar energy resource based on FY-4A ground radiation products

2.2.2 插补方案

卫星反演辐照度数据一般为小时数据, 在每3 h(02:00, 05:00, 08:00)除整点外, 在整点前后15 min各有一次观测, 进行短临预报的数据需要15 min的数据, 在进行短临预报前, 对无卫星数据的时次先进行数据插补, 本文采取线性插补方案, 将小时数据插补到15 min数据。

2.2.3 均生函数短临预报方法

江滢等(2013)提到考虑到晴空指数是一个相对稳定的物理量, 一般来说, 数分钟或数十分钟内变化较小, 因此, 利用晴空指数作为短临预报方法的基础物理量, 以尽量减少天气突变对分钟级太阳能预报的影响。
本文通过《太阳能资源评估方法(GB/T 37526-2019)》(全国气候与气候变化标准化技术委员会, 2019)中附录B计算瞬时地外水平面太阳辐照度EI; 然后利用EI将观测数据和预报数据转换成晴空指数:
K = G H I E I
式中: K为晴空指数; GHI水平面总辐射辐照度; EI是瞬时地外水平面太阳辐照度, 指到达地球大气上界的太阳辐射能量, 是一个不受大气影响, 随时间而变得辐射量。对白天时段的晴空指数利用常蕊等(2013)的短临预报方案进行短临外推预报, 对获得的预报晴空指数, 再利用式(1), 输入预报时次的EI, 获得预报的GHI
本文参考常蕊等(2013)的基于均生函数的短临预报方案, 该方案分三块: 计算序列的延拓序列、 利用灰色关联度挑选显著周期和利用多元线性回归进行建模预报。该方法既可以进行多步预报, 又可以较好地预报极值。
K _ _ l ( i ) = 1 n l j = 0 n l - 1 K ( i + j l ) ,    i = 1 , ,   l ;   1 l m
ζ i ( k ) = m i n i m i n k | K 0 ( k ) - K i ( k ) | + P m a x i m a x k | K 0 ( k ) - K i ( k ) | | K 0 ( k ) - K i ( k ) | + P m a x i m a x k | K 0 ( k ) - K i ( k ) |
γ i = 1 n k = 1 n ζ i ( k )
K ̂ ( n + q ) = a 0 + i = 1 k a i f i ( n + q ) q = 1,2 , , 16
(5)
式(2)为均生函数的构造, 本文中设置n为8; K n为预报起报时刻前2 h逐15 min 8个时次的晴空指数时间序列; K l ¯为构造的均值生成函数, 将均生函数延长到整个数轴, 可以获得周期性的均生函数延拓序列。式(3)为灰色关联系数计算, 式(4)为灰色关联度计算, 可以获得均生函数和预报序列之间的关联度大小; 通过灰色关联分析, 可以挑选出优势周期, 利用多元线性回归方法可以建立预报模型, 即为式(5); 设模型用作q步预报, 将入选的均生函数作q步外延, 本文中设置q最大为16, 则 K ̂为获得的起报时刻后4 h逐15 min 16个时次的晴空指数; 公式详细意义可参见常蕊等(2013)。值得一提的是, 基于均生函数的短临预报方法是一种数学统计方法, 其预报精度只和历史观测序列本身有关, 在本文中, 即和卫星反演辐照度的精度有关, 在沈钟平和张华(2009)以及吴北婴(1999)的研究中指出, 到达地面的太阳辐射受云、 气溶胶、 吸收气体、 地表、 太阳入射和卫星观测角度等多种参数的综合影响, 经过离散纵标法辐射传输软件包的反演计算时已考虑云、 气溶胶、 水汽等的影响。

2.2.4 订正方法

梁进秋等(2020)的数据统计分析, 太原站的卫星辐射反演值比地面观测值整体偏高, 说明卫星反演辐射存在系统性的偏差, 为解决这一问题, 本文在短临预报基础上, 对预报辐射进行订正。用当前时刻的卫星反演值与地面观测值的偏差订正预报时刻的预报值, 并随着时间推移一直滚动订正, 见式(6)式(7)。先预报再订正的方案, 首先保持了风四产品对辐射连续性观测、 对天气波动变化能较好把握的优点, 然后实时订正, 又消除了一部分风四产品的系统性误差, 这是该方法的优势所在。
δ = G H I S A T - G H I O B S
G H I f s t ' ' q = G H I f s t q - δ ,    q = 1,2 , , 16
式中: G H I S A T为预报当前时刻的卫星反演总辐射; G H I O B S为预报当前时刻的地面观测辐射; δ为两者的偏差; G H I f s t q为预报的总辐射; q为预报步长, 本文q最大为16步; G H I f s t ' ' ( q )为订正后的预报总辐射。

2.3 检验方法

本文做的短临预报方案为逐15 min滚动预报未来4 h的总辐射, 时间分辨率为15 min, 即每次预报未来16个时间步长, 每个时间步长有短临预报总辐射、 订正后的短临预报总辐射和对应的观测总辐射, 试验时间为2019年1月1日至12月31日。
本文设计3种检验手段: (1)典型天气个例分析。挑选3种典型天气(晴天、 云天、 阴天)的个例进行分析, 检验本方法在不同天气条件下的预报能力。(2)不考虑时间变化的16个预报时间步长综合统计分析。每次短临预报, 都有一条16个预报时间步长的预报曲线, 不考虑每条预报曲线的时间变化, 只统计每个预报时间步长的统计结果, 比如将所有预报曲线的第1个预报时间步长的结果进行统计, 如此类推, 检验该方法对预报未来4 h辐射的能力。(3)整点预报时刻的全年综合统计分析。以每天11:00起报为检验对象, 预报未来4 h的结果, 整点预报时刻为12:00、 13:00、 14:00和15:00, 对2019年每天11:00起报的4 h整点预报结果进行统计分析。
文中涉及的检验用到以下指标。
平均相对误差:
A R E = 1 N i = 1 N F G H I i - O G H I i O G H I i × 100 %
均方根误差:
R M S E = i = 1 N ( F G H I i - O G H I i ) 2 N
相关系数:
r = i = N ( F G H I i - F G H I ¯ ) ( O G H I i - O G H I ¯ ) i = N ( F G H I i - F G H I ¯ ) 2 i = N ( O G H I i - O G H I ¯ ) 2
式中: O G H I为地面观测总辐照度; O G H I ¯为地面观测总辐照度的平均; F G H I为预报总辐照度; F G H I ¯为预报总辐照度的平均; N为统计样本数。

3 结果分析

3.1 不同天气下个例分析

为检验试验效果, 挑选3种典型天气(晴天、 云天、 阴天)的个例进行分析, 因中午时分辐射最大, 挑选以11:00为预报时间, 预报未来4 h作为研究对象。本文以小时云量(100为全云)为指标, 分析3种特殊天气条件下短临预报能力(程兴宏等, 2014): (1)11:00 -15:00每个时次的云量≤10, 定义为晴天; (2)11:00 -15:00每个时次的云量在40~60间, 定义为多云; (3)11:00 -15:00每个时次的云量≥90, 定义为阴天。经挑选, 太原站的晴天典型日为3月23日, 多云典型日为5月17日, 阴天典型日为6月5日。
从太原站晴天个例[表1图2(a), 预报起始时间为2019年3月23日11:00, 预报未来4 h]可以看出, 晴天辐射的预报曲线相对比较光滑, 和地外太阳辐射的趋势比较一致, 订正前的短临预报总辐射整体比观测偏低, 最大偏差120 W·m-2左右, 经过订正后的短临预报总辐射得到了有效的改进, 每个时刻都得到了改进, 最大偏差缩小到40 W·m-2左右。晴天经过订正后的短临预报辐射, 太原站得到了有效的改进。
表1 太原站的3种典型日的订正前后不同预报时长的平均总辐射与观测总辐射对比

Table 1 Comparison of mean global radiation and observed global radiation with different forecast duration before and after correction of three typical days at Taiyuan station

北 京 时

观测值

/(W·m-2

晴天(3月23日) 多云天(5月17日) 阴天(6月5日)
订正前 订正后

观测值

/(W·m-2

订正前 订正后 观测值 /(W·m-2 订正前 订正后
预报值 /(W·m-2 相对误差 /% 预报值 /(W·m-2 相对误差 /% 预报值 /(W·m-2 相对误差 /% 预报值 /(W·m-2 相对 误差 /% 预报值 /(W·m-2 相对 误差 /% 预报值 /(W·m-2 相对 误差 /%
11:00 826 750.8 -9.1 803.0 -2.8 860 948.6 10.3 870.9 1.3 109 442.3 305.8 157.9 44.9
11:15 848 770.8 -9.1 823.0 -2.9 878 959.7 9.3 882.0 0.5 110 457.6 316.0 173.2 57.5
11:30 856 791.1 -7.6 843.3 -1.5 920 979.9 6.5 902.2 -1.9 146 456.6 212.7 172.2 17.9
11:45 885 803.7 -9.2 855.9 -3.3 917 988.4 7.8 910.7 -0.7 151 435.8 188.6 151.4 0.3
12:00 890 809.8 -9.0 862.0 -3.1 934 994.7 6.5 917.0 -1.8 132 499.4 278.3 215.0 62.9
12:15 901 824.0 -8.5 876.2 -2.8 914 999.8 9.4 922.1 0.9 92 440.8 379.1 156.4 70.0
12:30 912 822.2 -9.8 874.4 -4.1 914 1003.5 9.8 925.8 1.3 114 468.2 310.7 183.8 61.2
12:45 903 821.8 -9.0 874.0 -3.2 947 994.4 5.0 916.7 -3.2 130 474.8 265.2 190.4 46.5
13:00 904 821.6 -9.1 873.8 -3.3 928 994.7 7.2 917.0 -1.2 100 464.4 364.4 180.0 80.0
13:15 894 813.4 -9.0 865.6 -3.2 923 983.0 6.5 905.3 -1.9 76 434.7 472.0 150.3 97.8
13:30 875 798.7 -8.7 850.9 -2.8 906 969.3 7.0 891.6 -1.6 83 488.4 488.4 204.0 145.8
13:45 867 792.0 -8.7 844.2 -2.6 893 954.3 6.9 876.6 -1.8 83 422.7 409.3 138.3 66.6
14:00 840 769.8 -8.4 822.0 -2.1 857 938.0 9.5 860.3 0.4 102 440.0 331.4 155.6 52.5
14:15 814 749.2 -8.0 801.4 -1.5 855 909.9 6.4 832.2 -2.7 111 437.2 293.9 152.8 37.7
14:30 776 728.5 -6.1 780.7 0.6 835 890.5 6.6 812.8 -2.7 150 418.8 179.2 134.4 -10.4
14:45 740 700.8 -5.3 753.0 1.8 789 860.3 9.0 782.6 -0.8 145 383.6 164.6 99.2 -31.6
图 2 太原站典型天气下的订正前后总辐射对比

Fig.2 Comparison of global radiation before and after correction in typical weather of Taiyuan station

从太原站多云天个例[表1图2(b), 预报起始时间为2019年5月17日11:00, 预报未来4 h]可以看出, 观测总辐射受云的遮挡影响比较大, 呈现较明显的波动性, 短临预报曲线也呈一定波动性, 并且整体大于观测总辐射, 整体偏大约70 W·m-2, 经过订正后的短临预报总辐射, 有效改进了预报效果, 和观测整体比较吻合, 平均绝对误差下降到14 W·m-2。多云天经过订正后的短临预报辐射, 太原站得到了有效的改进, 订正后最大误差不到14 W·m-2
从太原站阴天个例[表1图2(c), 预报时间为2019年6月5日11:00, 预报未来4 h]可以看出, 阴天的观测总辐射被极大削弱, 最大在150 W·m-2, 而且也呈现明显的波动性, 受天气影响极大, 短临预报总辐射整体偏大, 存在一个系统性偏差, 约333 W·m-2。经过订正后的短临预报总辐射, 有效改进了预报效果, 和观测整体比较吻合, 平均绝对误差下降到56 W·m-2。阴天经过订正后的短临预报辐射, 太原站得到了有效的改进。
从对太原站的3种典型天气的个例分析可以看出, 短临预报的总辐射和观测总辐射存在偏差, 经过订正后的短临预报总辐射明显改善了辐射预报。

3.2 不考虑时间变化的16个预报时间步长综合统计分析

从太原站2019年订正前后卫星预报总辐射和观测总辐射的不同预报时间步长的平均总辐射、 平均相对误差、 均方根误差和相关系数(表2)及其变化曲线(图3)可以看出: 不同预报时间步长的总辐射都比较稳定, 呈逐步下降的趋势。太原站短临预报的总辐射从第1个预报时间步长的656 W·m-2下降到第16个预报时间步长的565 W·m-2, 订正后的短临预报平均总辐射从第1个预报时间步长的569 W·m-2下降到第16个预报时间步长的491 W·m-2, 观测平均总辐射从第1个预报时间步长的569 W·m-2下降到第16个预报时间步长的467 W·m-2。订正后的短临预报总辐射和观测总辐射在每个预报时间步长都吻合得非常好, 有效地改善了短临预报的总辐射。
表2 太原站2019年订正前后不同预报步长的平均总辐射与观测总辐射对比

Table 2 Comparison of average global radiation and observed global radiation of different prediction step before and after the correction of Taiyuan station in 2019

预报时间 步长(间隔 15 min) 观测值 /(W·m-2 订正前 订正后
预报值 /(W·m-2 平均相对 误差/% 均方根误差 /(W·m-2 相关 系数 预报值 /(W·m-2 平均相对 误差/% 均方根误差 /(W·m-2 相关 系数
1 569.7 656.0 44.7 166.4 0.84 569.5 21.3 132.2 0.87
2 570.7 654.5 44.3 165.2 0.84 568.3 23.7 143.4 0.85
3 569.6 653.7 45.3 169.9 0.82 567.9 26.6 153.9 0.82
4 568.2 654.4 46.1 173.8 0.81 569.8 28.6 158.2 0.81
5 565.7 646.5 48.1 180.3 0.79 563.1 32.4 170.3 0.78
6 560.9 645.9 47.8 176.5 0.81 563.6 32.5 167.1 0.79
7 554.2 639.3 50.7 185.8 0.79 558.3 35.7 176.5 0.77
8 546.4 630.4 51.2 184.5 0.8 551.0 37.2 177.7 0.77
9 537.5 622.6 53 187.2 0.8 544.4 39.8 182.5 0.76
10 527.3 617.2 54.6 188.8 0.8 540.7 41.6 181.4 0.77
11 517.7 601.4 58.2 195.2 0.78 526.0 45.8 191 0.75
12 506.7 598.4 58.9 191.2 0.81 524.0 46.3 183.6 0.77
13 496.9 588.4 62.2 198.8 0.78 514.7 49.9 192.3 0.75
14 486.1 576.4 63.4 196.5 0.79 502.4 50.9 190.7 0.75
15 477.0 568.9 64.9 196.9 0.78 495.2 52.8 190 0.74
16 467.3 565.4 66.4 198.8 0.78 491.6 54.3 189.9 0.74
图3 太原站2019年订正前后预报总辐射和观测总辐射的不同预报步长的统计

Fig.3 Statistical chart of forecast and observation global radiation before and after correction in 2019

随着预报时间步长的加大, 短临预报总辐射的平均相对误差也在加大, 但总体比较平缓。太原站从预报第1时间步长的44.7%增大到第16个时间步长的66.4%, 订正后的短临预报总辐射的平均相对误差也得到极大的改善, 和订正前的平均相对误差随预报时间步长的加大而缓慢加大趋势一致, 从第一时间步长的21.3%增大到第16个时间步长的54.3% 。每个时间步长都得到有效改善。
不同预报时间步长的短临预报和订正后的短临预报都比较稳定, 太原站短临预报的均方根误差随着预报时间步长的加大而加大, 但总体比较平缓, 从预报第1时间步长的166.4 W·m-2增到第16个时间步长的198.8 W·m-2, 订正后的短临预报总辐射的均方根误差得到了改善, 在每个预报时间步长都得到了改善, 并且也是随着预报时间步长的加大而加大, 从第1个时间步长的132.2 W·m-2增加到第16个时间步长的189.9 W·m-2
短临预报的不同时间步长的相关系数逐步下降趋势; 太原站从0.84逐步降到0.78。订正后的短临预报的不同时长的相关系数有同样的逐步下降趋势, 从0.86下降到0.71。在前面2个时间步长, 订正后的相关系数大于未订正的, 第3~4个时间步长相关系数订正前后一致, 第5~16个时间步长, 订正后的相关系数略低于订正前的。
从上述的各种统计分析可以看出, 订正后的短临预报结果能有效改善短临预报的结果, 有效地提升了辐射预报的能力。

3.3 整点预报时刻的全年综合统计分析

从2019年太原站总辐射不同预报时刻(12:00、 13:00、 14:00、 15:00)的预报值、 订正值和观测值以及预报值和观测值之差、 订正值和观测值之差的逐日变化曲线(图4)可以看出, 以12:00为例, 总体来看, 总辐射的预报值大于观测值, 订正后总辐射与观测值较为接近, 偏差减小。13:00, 14:00和15:00有类似的表现。并且2019年太原站不同预报时刻的订正前后的预报总辐射与观测总辐射的相关系数、 均方根误差、 平均相对误差(图5)可看出, 相关系数随着预报时刻增加而减小, 4个预报时次的相关系数均大于0.72, 订正后前3个时次的相关系数都有提高, 12:00从0.84提高到0.88, 13:00从0.78提高到0.82, 14:00从0.77提高到0.79, 15:00略有下降, 从0.74下降到0.72。均方根误差随着预报时刻增加而增加, 4个预报时次的均方根误差均在200 W·m-2以下, 订正后均方根误差都有明显改进, 12:00从162.7 W·m-2下降到129.3 W·m-2, 13:00从196.8 W·m-2下降到165.8 W·m-2, 14: 00从194.0 W·m-2下降到168.3 W·m-2, 15:00从197.1 W·m-2下降到180.9 W·m-2。平均相对误差随着预报时间增加而增加, 4个预报时次的平均相对误差在80%内, 订正后平均相对误差有显著减少, 12:00从46.0%下降到24.0%, 13:00从60.2%下降到33.7%, 14:00从64.4%下降到37.3%, 15:00从76.8%下降到46.5%。
图4 2019年太原站总辐射不同预报时刻的预报值、 订正值和观测值(a, c, e, g)以及预报值和观测值之差、 订正值和观测值之差(b, d, f, h)的逐日变化曲线

Fig.4 The daily variation curves of the predicted, the corrected and the observed value (a, c, e, g), the difference between the predicted and the observed value, the difference between the corrected and the observed value (b, d, f, h) at different forecast time of the global radiation of Taiyuan station in 2019

图5 2019年太原站不同预报时刻的订正前后的预报总辐射与观测总辐射的相关系数(a)、 均方根误差(b)、 平均相对误差(c)

Fig.5 Correlation coefficient (a), RMSE (b) and average absolute relative error (c) of global radiation before and after correction at different forecast time of Taiyuan station in 2019

4 结论和讨论

以山西太原国家级气象观测站为例, 建立基于均生函数法的短临预报方法, 输入过去2 h逐15 min的卫星数据, 获得未来4 h的逐15 min预报总辐射, 并在此基础上滚动订正获得订正后的短临预报辐射结果。主要结论如下:
(1) FY-4A是新一代静止气象卫星, 高时空分辨率的特点为短临预报提供了可能。经过检验, FY-4A地面辐射相对于地面观测系统性偏高, 需要订正后才能应用。本文先预报再订正, 首先保持了风四产品对辐射连续性观测、 对天气波动变化能较好把握的优点, 然后利用前一时刻的误差进行订正, 又消除了一部分风四产品的系统性误差——这是该方法的优势所在。
(2) 从个例分析看, 基于卫星辐射产品的短临预报能抓住天气变化带来的影响, 经过订正的短临预报辐射能得到有效改进。晴天情况下, 太原站从订正前的最大偏差120 W·m-2左右改进到订正后的40 W·m-2左右。多云天情况下, 平均相对误差下降到14 W·m-2。阴天情况下, 平均相对误差下降到56 W·m-2
(3) 从不考虑时间变化的16个预报时间步长的综合统计分析看, 订正后的短临预报总辐射和观测总辐射在每个预报时间步长都吻合得非常好, 相对误差从44.7%~66.4%下降到21.3%~54.3%, 均方根误差从166.4~198.8 W·m-2下降到132.2~189.9 W·m-2, 每个时间步长都得到有效改善, 订正后的相关系数总体比订正前有所下降。
(4) 从整点预报时刻的全年综合统计分析看, 总辐射的预报值大于观测值, 订正后总辐射与观测值较为接近, 偏差减小。订正后的相关系数总体比订正前提高, 从0.77~0.84提高到0.79~0.88, 仅在15:00略有下降。均方根误差从162.7~197.1 W·m-2下降到129.3~180.9 W·m-2, 平均相对误差从46.0%~76.8%下降到24.0%~46.5%。
本方法的优点, 首先从空间上, 基于风四卫星辐射数据, 在全圆盘任意点上都能进行预报, 其次从时间上, 卫星数据时间连续性好, 可做任意光伏电站任意时刻的短临预报。
尽管本方法获得以上定性和定量结论, 但本文仍存在以下几个问题需要讨论, 并值得进一步深入研究: (1)当前的订正方法依赖于当地的气象站数据, 在没有地面实测辐射数据的情况下, 需建立更适宜的订正方法。(2)风四产品本身有一定滞后性, 通常滞后20 min, 导致预报时长第一小时本身效果好, 但在实际预报应用中来不及使用, 在实际应用中, 用第二个小时的结果。(3)本文所建的方法, 存在几方面的不确定性, 即: 本方法用于预报的历史逐15 min数据序列, 由小时数据插值获得, 存在误差; 基于均生函数的短临预报方法为纯统计的数学方法, 只与输入变量有关。从沈钟平和张华(2009)以及Zhang et al (2013)的研究中可以看出辐射受到云、 气溶胶、 水汽、 太阳天顶角等的影响, 虽然在卫星反演时已考虑, 但仍存在一定的不确定性; 另外一个不确定性, 在做订正时, 假定未来4 h系统误差不变, 这也是一个误差来源。今后将结合数值模式和地面全天空成像仪的观测资料进行预报和订正的改进。(4)本方法仅在太原站进行了检验, 未来将在更多的观测站点进行应用。
因此, 下一步有必要在本文研究的基础上, 对不同地区, 进一步考虑其天气气候特征以及云、 气溶胶、 水汽、 太阳天顶角等影响因子, 建立更适用于本地的辐射短临预报方法。

References

null
Ayet A Tandeo P2018.Nowcasting solar irradiance using an analog method and geostationary satellite images[J].Solar Energy, 164: 301-315.
null
Luis H C Sara G S Victor A G2019.A review of photovoltaic systems: design, operation and maintenance[J].Solar Energy, 188: 426-440.
null
Miller S D Rogers M A Haynes J M, et al, 2018.Short-term solar irradiance forecasting via satellite/model coupling [J].Solar Energy, 168: 102-117.
null
Sergio H A Andres S G Miguel R P, et al, 2020.A solar altitude angle model for efficient solar energy predictions[J].Sensors, 20: 1391.
null
Wan C Xu Z Pinson P, et al, 2014.Probabilistic forecasting of wind power generation using extreme learning machine[J].IEEE Transactions on Power Systems29(3): 1033-1044.
null
Yang L W Gao X Q Hua J J, et al, 2020.Very short-term surface solar irradiance forecasting based on FengYun-4 geostationary satellite[J].Sensors, 20: 2606.
null
Zhang H Yin Q Teruyuki N, et al, 2013.Influence of changes in solar radiation on surface temperature in China[J].Acta Meteorologica Sinica27(1): 87-97.DOI: 10.1007/s13351-013-0109-8 .
null
常蕊, 朱蓉, 柳艳香, 等, 2013.基于均生函数的风电场风速短临预报模型[J].气象39(2): 226-233.
null
程兴宏, 刘瑞霞, 申彦波, 等, 2014.基于卫星资料同化和LAPS-WRF模式系统的云天太阳辐射数值模拟改进方法[J].大气科学38(3): 577-589.
null
程序, 周海, 王知嘉, 等, 2011.晴空工况光伏超短期功率预报方法[J].可再生能源29(6): 142-145.
null
崔林丽, 郭巍, 葛伟强, 等, 2020.FY?4A卫星云顶参数精度检验及台风应用研究[J].高原气象39(1): 196-203.DOI: 10.7522/j.issn.1000-0534.2019.00065 .
null
崔杨, 成驰, 陈正洪, 2018.考虑大气透明度修正的光伏发电超短期功率预测方法[J].水电能源科学36(9): 201-204.
null
丁宇宇, 丁杰, 周海, 等, 2014.基于全天空成像仪的光伏电站水平面总辐射预报[J].中国电机工程学报34(1): 50-56.
null
高晓清, 蒋俊霞, 汪宁渤, 等, 2018.基于TSI地基云图的云量反演方法研究[J].太阳能学报39(1): 95-100.
null
黄守友, 徐国强, 2020.FY4A的LMIE闪电数据对云信息初始化的影响及数值试验[J].高原气象39(2): 378-392.DOI: 10. 7522/j.issn.1000-0534.2019.00110 .
null
江滢, 申彦波, 党军, 2013.两种太阳总辐射分钟级预报方法地比较[J].气象与环境学报29(3): 85-91.
null
姜红, 何清, 曾晓青, 等, 2021.基于随机森林和卷积神经网络的FY-4A号卫星沙尘监测研究[J].高原气象40(3): 680-689.DOI: 10.7522/j.issn.1000-0534.2020.00060 .
null
蒋俊霞, 高晓清, 吕清泉, 等, 2020.基于地基云图的云跟踪与太阳辐照度超短期预报方法研究[J].太阳能学报41(5): 351-358.
null
梁进秋, 申彦波, 胡丽琴, 等, 2020.FY-4A地表太阳入射辐射产品在山西高原的适用性研究[J].气象46(12): 1575-1585.
null
邱宗旭, 杨琳, 李兴荣, 等, 2019.从辐射强度变化分析深圳太阳辐射气象条件[J].太阳能学报40(3): 724-732.
null
全国气候与气候变化标准化技术委员会, 2019. 太阳能资源评估方法: GB/T 37526-2019 [S].北京: 中国标准出版社, 1-24.
null
申彦波, 程兴宏, 赵东, 2013.气象与太阳能生态城[M].北京: 中国建筑工业出版社, 16-25.
null
沈钟平, 张华, 2009.影响地面太阳辐射及其谱分布的因子分析[J].太阳能学报30(10): 1209-1215.
null
石广玉, 2007.大气辐射学[M].北京: 科学出版社.
null
王长贵, 王斯成, 2009.太阳能光伏发电应用技术 (第二版)[M].北京: 化学工业出版社.
null
吴北婴, 1999.大气辐射传输实用算法[M].北京: 气象出版社.
null
吴林栋, 2020.地面气象观测站探测环境调查评估报告书(太原站)[R].北京: 中国气象局.
null
尹青, 张华, 何金海, 2011.近48年华东地区地面太阳总辐射变化特征和影响因子分析[J].大气与环境光学学报6(1): 37-46。
null
张凤霞, 李正强, 李凯涛, 等, 2016.京津唐地区气溶胶直接辐射强迫的遥感观测[J].中国科学院大学学报33(2): 155-161.
null
张华, 荆现文, 彭杰, 等, 2019.云辐射与气候[M].北京: 气象出版社.
null
张华, 王志立, 赵树云, 等, 2017.大气气溶胶及其气候效应[M].北京: 气象出版社.
null
张志清, 陆风, 方翔, 等, 2017.FY-4卫星应用和发展[J].上海航天34(4): 8-19.
null
中国电力企业联合会, 2013. 光伏发电站功率预测系统技术要求: NB/T 32011-2013 [S].北京: 中国电力出版社, 1-8.
Outlines

/