WO2014190651A1

WO2014190651A1 - 一种基于地基云图的光伏功率预测方法

Info

Publication number: WO2014190651A1
Application number: PCT/CN2013/084700
Authority: WO
Inventors: 陈志宝; 程序; 周海; 丁杰; 朱想; 王知嘉; 彭佩佩; 江星星; 何洁琼; 孙檬檬
Original assignee: 国家电网公司; 中国电力科学研究院
Priority date: 2013-05-31
Filing date: 2013-09-30
Publication date: 2014-12-04
Also published as: CN103353952A; CN103353952B

Abstract

本发明提供一种基于地基云图的光伏功率预测方法，包括以下步骤：步骤1：云团运动预估；步骤2：对地面辐照度进行预测；步骤3：建立综合光电转换模型；步骤4：进行光伏功率预测。本发明综合运用多种数字图像处理技术，突出问题需求，有很强的针对性和实用性，与传统的统计方法结合为光伏功率精细化预测提供了有效的技术支撑和方法参考。且本发明综合多因素多步建模，方法的描述已涵盖光伏预测领域所有关键环节，有较强的普适性和鲁棒性。本发明为了在使用中结合使用的具体场景和应用对象，有些环节预留了调节参数，操作方便并保证预测效果。

Description

一种基于地基云图的光伏功率预测方法

技术领域

本发明属于新能源光伏发电技术领域，具体涉及一种基于地基云图的光伏功率预测方法。背景技术

近年来，大规模光伏电站接入电网，由于光伏发电输出功率具有随机性和波动性，对电网安全稳定和经济运行造成影响。对光伏电站的输出功率进行准确预测，为电力调度提供重要的决策支持，统筹安排常规电源和光伏发电的协调配合，不仅降低电力***运行成本，而且有利于充分利用光伏资源，获得更大的经济效益和社会效益。

地面辐照度作为光伏电站输出功率的主要影响因素之一，它的不确定性直接导致输出功率的随机性和波动性。而云作为影响太阳地面辐照量的主要气象要素，其生消和移动变化是地面辐照度变化不确定性的主要原因。目前，光伏功率预测方法主要有统计分析或机器学习方法，它们基于历史气象要素数据和光伏电站输出功率数据进行预测，但均未涉及云图数据，对云的遮挡导致的功率突变很难学习导致预测效果欠佳。国外有学者根据气象卫星云图预估云层移动情况进行太阳辐照度预测，但其时空分辨率相对较低,为了满足区域内光伏电站超短期 0-4小时内分钟级功率预测的要求，需要局地高时空分辨率的云图像数据。

近年来，随着电荷耦合器件 CCD等硬件技术的发展和数字图像处理技术的不断完善，很多地基遥感测云仪器研制成功，比如全天空成像仪（total sky imager, TSI)，该仪器能实现对全天空云的自动化观测和实时采集，最小时间分辨率可达 30秒，能很好地实现光伏电站上空云图的实时监测和采集。发明内容

为了克服上述现有技术的不足，本发明提供一种基于地基云图的光伏功率预测方法，利用数字图像处理技术对地基云图进行处理和分析，对未来时刻云团运动和太阳遮挡状况进行预测，相当于将云这一主要气象要素从众多影响辐照度的随机因素中单独分离出来进行建模。同时结合衡量太阳通过大气层时衰减情况的辐射清晰度指数实现未来时刻地面辐照度的预测，最后通过光电转换模型实现光伏功率超短期 1小时内分钟级光伏功率精细化预测。

为了实现上述发明目的，本发明采取如下技术方案：

提供一种基于地基云图的光伏功率预测方法，所述方法包括以下步骤：

步骤 1 : 云团运动预估；

步骤 2: 对地面辐照度进行预测；

步骤 3 : 建立综合光电转换模型；

步骤 4: 进行光伏功率预测。

所述步骤 1包括以下步骤：

步骤 1-1 : 确定图像上太阳位置；

所述步骤 1-1中，用天顶角和方位角来确定图像上太阳位置；

假设观测地的纬度为，赤纬角为，时角为《，天顶角和方位角分别用和; ^表示为：

δ ₌ 23·45。·ή_ηί2_πχ^{284 +} (2) ω = {ΐ— 12、χ15° (3) 其中， ώϊ)^表示一年中的日期序号，为时间；天顶角为和方位角分别为和; ^时，距离图像中心点的距离为 R，于是即可确定某一时刻图像上太阳位置；距离图像中心点的距离 R表示为：

R = -17.59ΘΙ + 118.1.35θ_ζ (4)。步骤 1-2: 对地基云图原始图像进行预处理；

所述步骤 1-2包括以下步骤：

步骤 1-2-1: 对地基云图原始图像进行图像复原；

首先，通过遮光带覆盖像素点范围和镜头支臂的像素点范围确定遮光带和镜头支臂的掩模图像；根据已确定的图像上太阳位置，标定出复原图像对应时刻太阳在图像中的太阳位置点，然后由中心点到太阳位置点至到圆周边沿点确定一条直线，然后根据图像形态学膨胀运算扩展至遮光带宽度，从而确定遮光带覆盖像素点范围；同时由中心点垂直向上至到圆周边沿点确定另一条直线，然后根据图像形态学膨胀运算把该直线扩展至宽度为镜头支臂宽度，从而确定镜头支臂的像素点范围。

然后在地基云图原始图像上，与掩模图像对于遮光带覆盖范围的像素点和镜头支臂范围的像素点的周边像素点进行线性插值，最终生成复原图像 /；

步骤 1-2-2: 对所述复原图像 /进行中值滤波，得到中值滤波图像；

M (i, j) = Med{I{i + r,j + s)} (5) 其中， Μ(ί, 表示中值滤波图像在（, )点的灰度值， /( + r，j' + 表示复原图像 /在

( +r,j'+ 点的灰度值；滤波窗口为 5x5领域时， r，^{-2，-l，0，l，2}_; 滤波窗口为 3x3领域时， r， ^ e {-1,0,1}； Med{I(i + r, j + s)表示随着 r , s取不同值而得到的 25个序列或 9个序列的中值；

步骤 1-2-3: 对中值滤波图像进行锐化；

对中值滤波图像进行锐化后的锐化图像假定为 g(，_y)，则有： g(x,y)=M(x,y)-V²M(x,y) (6) 其中， Μ( ，})表示中值滤波图像，且有 (7) 。

步骤 1-3: 对预处理后的锐化图像 g cy)进行二值化处理，并提取云团质心；所述步骤 1-3包括以下步骤：

步骤 1-3-1 : 对预处理最终得到的锐化图像 g c y)进行二值化处理; 由下式得到 (i, j)点的二值图像：

Binary (i, j) ( 8 )

其中， Bi腿 ry(i, j)表示二值图像 Binary在 (i, j)点的像素值， R(i， j), B(i, j)分别表示锐化图像 g(， 3；)红通道分量和蓝通道分量在 (j, j)点的灰度值； y表示锐化图像 y)红蓝通道分量比的阈值，由如下公式计算：

γ = max—level {RBRJiist) + TH ( 9 ) 其中， RBR— toi表示红蓝分量比图像的灰度直方图， .level表示灰度直方图中频度最大的灰度值， H表示阈值微调量;

步骤 1-3-2: 进行云团检测和标注;

二值图像 V »7中，若像素值为 1且坐标为的像素点为 Ρ，则搜索 Ρ的 8连接点，

8连接点的坐标集为：

对 Ρ点的 8连接点的坐标集对应的像素点中选择像素值为 1的点，并分别对其再次展开 8连接点搜索，如此反复，并且每个像素点不能重复被搜索，直至所有连接点的 8连接点均为已搜索的点或者像素值为 0的点，则结束搜索并标注云团；

云团的标注按照检索的顺序，第一次检索到的团块标记为 1，第二次检索到的团块标记为 2，以此类推；

步骤 1-3-3 : 云团质心提取；

设 ^， Al₂为顶点的矩形框经过云团最小外接矩形的左上角和右下角像素点，若点坐标为 041χ₁ ， Al₂点坐标为（Alx₂, Al;y₂)，云团的质心坐标为（Abe, Aly)，则有:

Al¾+Alx₂

Alx = -

2 ( 10 )

Al^+Aly,

Aly =

2

按照上述云团质心提取的方法，依次提取云团 …和 Am共 m个云团的质心。步骤 1-4: 云团运动预估：

所述步骤 1-4包括以下步骤：

步骤 1-4-1 : 定义云团运动方向和速度；

建立以球面图像中心为坐标原点的极坐标系，球面图像中心的坐标为（ , ），像素点 P 的原始坐标为（ ,_y)，其中 ≠ ; 则像素点/³在极坐标下的极径和极角分别表示为: p = ^(x-x_c)² + (y- y_c)²

Θ = arctan (-— ~~―)

(x_t,y_t)和分别为时刻和 +l时刻某云团的质心坐标；则 t时刻该云团移动速度 V表示为：

v =l(x_t+1-x_t)² + (y_t+1-y_t)² /At (12) 其中 Δί为时刻和 +l时刻之间的采样时间间隔；

云团移动方向定义如下：

(1) 如果 x > X,，且 y_i+1≤y,，则有 ;

(2) 如果; c_i+1 =_¾ ，且 < y,，则有 θ = πΙ1

(3) 如果 χ_ί+1 < ，且 < y_t，则有 Θ = π/2+arctani ^

|¾+1 _

(4) 如果 χ_ί+1 <χ, ，且 y_i+1≥y,，则有 6> = +arctan^^t+1 ~^y^

(5) 如果 χ_ί+1 =χ, ，且 > y,，则有 θ = π/2；

(6) 如果 χ_ί+1 >χ, ，且 > y_t，则有 θ =

步骤 1-4-2: 进行云团运动预估；

假设二值图像 V »7有 W张，每张二值图像 V »7中云团个数有个，第个云团根据第张二值图像 V »7中的质心位置和第 + l张二值图像 V »7中的质心位置所得到的速度和角度分别为 V ^和 ^，其中： = l,2,...， m; = 1,2,...,π-1;

假设云团在相对较短时间内的运动为刚体运动，位于全天空成像仪视野范围内的云团近似匀速直线运动，通过求取 V ^和 ^平均值预估未来第个云团的移动速度和移动方向 θ_; , 具体表示为： (13) 。

所述步骤 2包括以下步骤：

步骤 2-1: 计算学习阶段太阳被云团遮挡时刻对应的辐射清晰度指数；

时刻天顶辐射 H。(0和大气质量 AM (0分别表示为：

H₀ (t) = I_sc (1 + 0.033co5(2^ · dayth 1365))cos0_z ( 14)

AM (t) = ^1229 + ( 6 cose_zf - 614cos0_z (15) 其中， /_sc =1367±5W/m²; 于是：时刻地面理论辐射值表示为：

I_theory(t) = H₀(t)IAM^c(t) (16) 在晴空状况下，令时刻地面理论辐射值等于地面实测辐射值。_fo(0根据式（16) 计算出 C，进而确定时刻地面理论辐射值；

学习阶段太阳被云团遮挡时刻对应的辐射清晰度指数用来衡量太阳通过大气层时云团对辐照度的衰减情况，其表示为：

KT(t) = K) (17) 其中， KT t、为学习阶段太阳被云团遮挡时刻对应的辐射清晰度指数；

步骤 2-2: 进行太阳遮挡预测；

通过判断二值图像上太阳位置处像素点是否为二值图像 V »7中云团的像素点，进而判断此刻太阳是否被云团遮挡，具体表示为：

ίθ, 如果 ( (t)+ex, y_s (t) + ex)点像素值为 0

sisn(t) = i (is

[l, 如果 (t) + (t) + e )点像素值为 1 其中 signW表示时刻太阳的遮挡信息， w'gn(i)=0表示太阳未被云团遮挡， w'gn(i)=l表太阳被云团遮挡；（x_s(0,y_s(0)表示时刻太阳在图像中的位置， ex的取值为 {-2， -1,0,1,2}；步骤 2-3: 对地面辐照度进行预测；假设预测阶段的辐射清晰度指数为 pKT, 其由学***均确定，表示为：

其中， N为学***均的权重系数，且满足 Λ<Λ.₊₁; 7^_;)为学习阶段太阳被云团遮挡时刻对应的第个辐射清晰度指数；

于是，预测：

其中， ρ 表示预测阶段 t时刻的地面理论辐照度。

所述步骤 3包括以下步骤：

步骤 3-1: 建立基础光电转换模型；

通过所读取的历史功率数据和辐射数据绘制散点图，若 p₀(o、 /(0分别为时刻的实测有功功率和实测地面辐照度，并根据散点图确定曲线类型和曲线拟合方程式中的阶数^ 所述曲线拟合方程式为

(0 -a l^q(t) + a «-1 ^q-¹ (t) +… + α₀ (21) 其中 ^,α^，…和 _Ω。均为曲线拟合方程式的系数;

根据上述曲线拟合方程式建立目标函数： ιηϊηΫ(α Ι"(ί)

+■■■ + a_(j-P_(j(t)f (22) 通过最小二乘法确定 i，…和 α。，从而确定时刻的实测有功功率 (0，表示为：

Ρ₀ (0 = f (Ht)) = a_qI^q(t) + a_qJ^q-¹ (0 +… + (23)

/(/(0)为时刻的实测有功功率 (0关于时刻的实测地面辐照度 /(0的函数；步骤 3-2: 结合基础光电转换模型建立综合光电转换模型；

设时刻的组建温度分别为 Γ(ί)，由温度导致的时刻的转换效率≠t、表示为：

77(Ο = 1+0·0035(Γ(Ο- r(t-l》 (24) 由此可得预测阶段时刻的实测有功功率，具体表示为：

其中， P(0表示预测阶段时刻的实测有功功率。

所述步骤 4中，通过建立光伏功率预测模型进行光伏功率预测；

所述光伏功率预测模型为： P{t) = (26)。

与现有技术相比，本发明的有益效果在于：

( 1 )本发明综合运用多种数字图像处理技术，突出问题需求，有很强的针对性，与传统的统计方法结合为光伏功率精细化预测提供了有效的技术支撑和方法参考。

(2)本发明综合多因素多步建模，方法的描述已涵盖光伏预测领域所有关键环节，有较强的普适性和鲁棒性。

(2)本发明采用必要且合理的假设，实现对问题简化的同时又能最大限度地准确描述问题。

( 3 ) 本发明为了在使用中结合使用的具体场景和应用对象，有些环节预留了调节参数，操作方便并保证预测效果。

(4)流程化的处理，详尽的算法描述和形式化表达，为工程实现提供了方便，具有很好的工程实用性。附图说明

图 1 是本发明实施例中某一时刻图像上太阳的位置示意图；

图 2 是图像形态学膨胀运算原理示意图；

图 3 是本发明实施例中掩模图像示意图；

图 4 是本发明实施例中对地基云图原始图像进行图像复原示意图；

图 5 是本发明实施例中 8连接点的连接方式示意图；

图 6 是本发明实施例中云团质心提取示意图；

图 7是本发明实施例中建立以球面图像中心为坐标原点的极坐标系中云团运动预估示意具体实施方式

下面结合附图对本发明作进一步详细说明。

步骤 1 : 云团运动预估；

步骤 2: 对地面辐照度进行预测；

步骤 3: 建立综合光电转换模型；

步骤 4: 进行光伏功率预测。

所述步骤 1包括以下步骤：

步骤 1-1 : 确定图像上太阳位置；

所述步骤 1-1中，用天顶角和方位角来确定图像上太阳位置；

假设观测地的纬度为，赤纬角为，时角为《，天顶角和方位角分别用和; ^表示为:

其巾 _:

S = 23.45^_Sin^,^{284 + dayth}] (2) ω = { — 12、χ15。 (3) 其中， ώϊ)^表示一年中的日期序号，为时间；天顶角为和方位角分别为和; ^时，距离图像中心点的距离为 R，于是即可确定某一时刻图像上太阳位置；由于天顶角在与图像半径有相关性，根据式（4)可将距离图像中心点的距离 R表示为：

R = -17.59ΘΙ + 118.1.35θ_ζ (4)。如图 1， Center点为图像中心点，这样就可以标定出某一时刻图像上太阳的位置。

步骤 1-2: 对地基云图原始图像进行预处理；

所述步骤 1-2包括以下步骤：

步骤 1-2-1: 对地基云图原始图像进行图像复原；

如图 2，图像形态学膨胀运算相当于用的中心在 A的边界走一圈，的边界扫过的区域加上 A本身的区域就是 A膨胀 B的结果。

M(i,j) = Med{I(i + r,j + s)} (5) 其中， M( , )表示中值滤波图像在（, 点的灰度值， /( + r，j' + 表示复原图像 /在

( +r,j'+ 点的灰度值；滤波窗口为 5x5领域时， r， {-2，-1，0，1，2}；滤波窗口为 3x3领域时， r， ^ e {-1,0,1}； Med{I(i + r, j + s)表示随着 r , s取不同值而得到的 25个序列或 9个序列的中值；

步骤 1-2-3: 对中值滤波图像进行锐化；

对中值滤波图像进行锐化后的锐化图像假定为 g( ，_y)，则有： g(x,y)=M(x,y)-V²M(x,y) (6) 其中， Μ( ，})表示中值滤波图像，且有 V²M₍x,) (7)

步骤 1-3: 对预处理后的锐化图像 g( ，})进行二值化处理，并提取云团质心;

所述步骤 1-3包括以下步骤：

步骤 1-3-1: 对预处理最终得到的锐化图像 g cy)进行二值化处理；由下式得到 (i, j)点的二值图像：

Binary (i, j) (8)

其中， Binary(i, j)表示二值图像 Binary在 (i, j)点的像素值， R(i， j), B(i, j)分别表示锐化图像 g(， 3；)红通道分量和蓝通道分量在 (j, j)点的灰度值； y表示锐化图像 y)红蓝通道分量比的阈值，由如下公式计算：

γ = max—level (RBR_hist) + TH (9) 其中， RBR— toi表示红蓝分量比图像的灰度直方图， .level表示灰度直方图中频度最大的灰度值， TH表示阈值微调量;

步骤 1-3-2: 进行云团检测和标注;

二值图像中，若像素值为 1且坐标为的像素点为 Ρ，则搜索 Ρ的 8连接点，其中 8连接点连接方式如图 5所示；

8连接点的坐标集为：

步骤 1-3-3: 云团质心提取；

如图 6，设， Al₂为顶点的矩形框经过云团最小外接矩形的左上角和右下角像素点，若 ^点坐标为 (ΑΙ^,ΑΙ^)， Al₂点坐标为 (A\x₂,A\y₂)，云团的质心坐标为 (A\x,A\y)，则有：

Al¾+Alx₂

Alx = -

2 (10)

Al^+Aly,

A\y =

2 按照上述云团质心提取的方法，依次提取云团 …和 Am共 m个云团的质心步骤 1-4: 云团运动预估：

所述步骤 1-4包括以下步骤：

步骤 1-4-1: 定义云团运动方向和速度；

如图 7，建立以球面图像中心为坐标原点的极坐标系，球面图像中心的坐标为（ , ），像素点 P的原始坐标为（ ,_y)，其中； C≠ ；则像素点 /⁵在极坐标下的极径和极角分别表示为:

p = (x- x_cY + (y- y_cY

(11) θ = arctan (- y-y,

-)

(x_t,y_t)和 _y_i+1)分别为时刻和 +1时刻某云团的质心坐标；则 t时刻该云团移动速度

V表示为：

v = ^ '{x_t+l-x_tf + {y_t+l-y_t) (12) 其中 Δί为时刻和 +l时刻之间的采样时间间隔；

云团移动方向定义如下：

(1) 如果 χ_ί+1 >χ,，且 ₊₁≤ ，则有

(2) 如果 x_t+1 = ，且 < y_t，则有 θ = πΙ1

(3) 如果 χ_ί+1 < χ,，且 < y_t，则有 θ =

(4) 如果 ,₊₁< ,，且 y,₊₁≥y,，则有 6>

(5) 如果 = · 且 yy", > ^ yy_t,，则有 θ = 3^/2； (6) 如果 χ_ί+1 > x,，且 > y_t，则有 θ =

步骤 1-4-2: 进行云团运动预估；

假设二值图像 V »7有 W张，每张二值图像 V »7中云团个数有个，第个云团根据第张二值图像 V »7中的质心位置和第 + l张二值图像 V »7中的质心位置所得到的速度和角度分别为 V 和 ^，其中： j = l,2"..,m'， i = l,2,...,n-l; 假设云团在相对较短时间内的运动为刚体运动，位于全天空成像仪视野范围内的云团近似匀速直线运动，通过求取 V ^和 ^平均值预估未来第个云团的移动速度和移动方向 θ_;, 具体表示为：

所述步骤 2包括以下步骤：

步骤 2-1: 计算学习阶段太阳被云团遮挡时刻对应的辐射清晰度指数;

t时刻天顶辐射 H。(0和大气质量 AM ( 分别表示为：

H₀ (t) = I_sc (1 + 0.033co5(2^- · dayth 1365))cos0_z (14)

AM (t) = ^1229 + ( 6 cose_zf - 614cos0_z (15) 其中， /_sc =1367±5W/m²; 于是：时刻地面理论辐射值 /_i¾∞ry(0表示为：

(t) = H₀(t)IAM (t) (16) 在晴空状况下，令时刻地面理论辐射值 /_ifc。„(0等于地面实测辐射值。_fc(0根据式（16) 计算出 C，进而确定时刻地面理论辐射值. (t)

学习阶段太阳被云团遮挡时刻对应的辐射清晰度指数用来衡] t太阳通过大气层时云团对辐照度的衰减情况，其表示

其中， κ 为学习阶段太阳被云团遮挡时刻对应的辐射清晰度指数；

步骤 2-2: 进行太阳遮挡预测；

0, 如果 ( (t)+ex, y_s (t) + ex)点像素值为 0

sign(t) = (18)

1, 如果 < (t) + ex, y_s (t) + ex)点像素值为 1 其中 signfy表示时刻太阳的遮挡信息， wXi)=0表示太阳未被云团遮挡， w'gn(i)=l表示太阳被云团遮挡；（x_s(0,y_s(0)表示时刻太阳在图像中的位置， ex的取值为 {-2， -1,0,1,2}；步骤 2-3: 对地面辐照度进行预测；

假设预测阶段的辐射清晰度指数为 pKT, 其由学***均确定，表示为：

于是，预测：

其中， ρ 表示预测阶段 t时刻的地面理论辐照度。

所述步骤 3包括以下步骤：

步骤 3-1: 建立基础光电转换模型；

通过所读取的历史功率数据和辐射数据绘制散点图，若 p₀(o、 /(0分别为时刻的实测有功功率和实测地面辐照度，并根据散点图确定曲线类型和曲线拟合方程式中的阶数^ 所述曲线拟合方程式为：

(0 = (0 + _q ^q-¹ (0 +… + α₀ (21) 其中 ^,α^，…和 _Ω。均为曲线拟合方程式的系数；

根据上述曲线拟合方程式建立目标函数： ιηϊηΫ(α Ι^ς(ί)

+ + α_ΰ- Ρ_ΰ(ΐ))² (22) 通过最小二乘法确定 i，…和 α。，从而确定时刻的实测有功功率 (0，表示为：

Ρ₀ (0 = f (Ht)) = a_qI^q(t) + a_qJ^q-¹ (0 +… + (23)

设时刻的组建温度分别为 Γ(ί)，由温度导致的时刻的转换效率 η(Λ表示为： 77(Ο = 1+0·0035(Γ(Ο - T(t - 1》 ( 24) 由此可得预测阶段时刻的实测有功功率，具体表示为：

其中， Ρ(0表示预测阶段时刻的实测有功功率。

所述步骤 4中，通过建立光伏功率预测模型进行光伏功率预测;

所述光伏功率预测模型为：

最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求范围当中。

Claims

1. 一种基于地基云图的光伏功率预测方法，其特征在于：所述方法包括以下步骤：步骤 1: 云团运动预估；

步骤 2: 对地面辐照度进行预测；

步骤 3: 建立综合光电转换模型；

步骤 4: 进行光伏功率预测。

2. 根据权利要求 1所述的基于地基云图的光伏功率预测方法，其特征在于：所述步骤 1 包括以下步骤：

步骤 1-1: 确定图像上太阳位置；

步骤 1-2: 对地基云图原始图像进行预处理；

步骤 1-3: 对预处理后的锐化图像 g( ，})进行二值化处理，并提取云团质心；

步骤 1-4: 云团运动预估。

3. 根据权利要求 2所述的基于地基云图的光伏功率预测方法，其特征在于：所述步骤 1-1中，用天顶角和方位角来确定图像上太阳位置；

假设观测地的纬度为，赤纬角为，时角为《，天顶角和方位角分别用和; ^表示为： θ_ζ = αηχο8( ηφήηδ + coscpcosScosco)

. ,cos5sin >、 (1) y_s = arcsin{ )

S = 23.45。._Si x^{284 + dayth} (2)

365

其中， ώϊ)^表示一年中的日期序号，为时间; 天顶角为和方位角分别为和; ^时，距离图像中心点的距离为 R，于是即可确定某一时刻图像上太阳位置；距离图像中心点的距离 R表示为：

R = -17.59ΘΙ + 118.1.35θ_ζ (4)。

4. 根据权利要求 2所述的基于地基云图的光伏功率预测方法，其特征在于：所述步骤 1-2包括以下步骤：

步骤 1-2-1: 对地基云图原始图像进行图像复原；

首先，通过遮光带覆盖像素点范围和镜头支臂的像素点范围确定遮光带和镜头支臂的掩模图像；然后在地基云图原始图像上，与掩模图像对于遮光带覆盖范围的像素点和镜头支臂范围的像素点的周边像素点进行线性插值，最终生成复原图像 /；

M(i, j) = Med{I(i + r,j + s)} (5) 其中， Μ(ί, 表示中值滤波图像在（, )点的灰度值， /( + r，j' + 表示复原图像 /在 ( +r,j'+ 点的灰度值；滤波窗口为 5x5领域时， r，^{-2，-l，0，l，2}_; 滤波窗口为 3x3领域时， r， ^ e {-1,0,1}； Med{I(i + r, j + s)表示随着 r , s取不同值而得到的 25个序列或 9个序列的中值；

步骤 1-2-3: 对中值滤波图像进行锐化；

5. 根据权利要求 4所述的基于地基云图的光伏功率预测方法，其特征在于：根据已确定的图像上太阳位置，标定出复原图像对应时刻太阳在图像中的太阳位置点，然后由中心点到太阳位置点至到圆周边沿点确定一条直线，然后根据图像形态学膨胀运算扩展至遮光带宽度，从而确定遮光带覆盖像素点范围；

同时由中心点垂直向上至到圆周边沿点确定另一条直线，然后根据图像形态学膨胀运算把该直线扩展至宽度为镜头支臂宽度，从而确定镜头支臂的像素点范围。

6. 根据权利要求 2所述的基于地基云图的光伏功率预测方法，其特征在于：所述步骤 1-3包括以下步骤：

步骤 1-3-1: 对预处理最终得到的锐化图像 g(，} 进行二值化处理；由下式得到 (i， j)点的二值图像：

_D. .、 \1； R(i,j)/B(i,j)>r

Binary (I, j) = < (8)

{0 ; R(i,j)/B(i,j)≤r 其中， Binary(i, j)表示二值图像 Binary在 (i， j)点的像素值， R(i, j), B(i, j)分别表示锐化图像 g( ， })红通道分量和蓝通道分量在 (i, j)点的灰度值； γ表示锐化图像 g(， y)红蓝通道分量比的阈值，由如下公式计算：

γ = max—level {RBRJiist) + TH (9) 其中， RBRJ^i表示红蓝分量比图像的灰度直方图， m«」_eve/表示灰度直方图中频度最大的灰度值， H表示阈值微调量；

步骤 1-3-2: 进行云团检测和标注；

二值图像 V »7中，若像素值为 1且坐标为（i'，J')的像素点为 Ρ，则搜索 Ρ的 8连接点，对 P点的 8连接点的坐标集对应的像素点中选择像素值为 1的点，并分别对其再次展开 8连接点搜索，如此反复，并且每个像素点不能重复被搜索，直至所有连接点的 8连接点均为已搜索的点或者像素值为 0的点，则结束搜索并标注云团；

步骤 1-3-3: 云团质心提取；

设 ^， Al₂为顶点的矩形框经过云团最小外接矩形的左上角和右下角像素点，若点坐标为 041χ₁ ， Al₂点坐标为（Alx₂, Al;y₂)，云团的质心坐标为（Abe, Aly)，则有：

Al¾+Alx₂

Alx = -

2

(10)

Al^+Aly,

Aly =

2

按照上述云团质心提取的方法，依次提取云团 …和 Am共 m个云团的质心。

7. 根据权利要求 2所述的基于地基云图的光伏功率预测方法，其特征在于：所述步骤 1-4包括以下步骤：

步骤 1-4-1: 定义云团运动方向和速度；

建立以球面图像中心为坐标原点的极坐标系，球面图像中心的坐标为（ , ），像素点 ρ 的原始坐标为（ ,_y)，其中 ≠ ; 则像素点/³在极坐标下的极径 p和极角分别表示为：

P X- x_cr + (y- y_c (11) θ = arctan (- y- y, -)

(x_t,y_t)和 _y_i+1)分别为时刻和 +l时刻某云团的质心坐标；则 t时刻该云团移动速度 V表示为：

V = j '(x_t+1 - x_tr + (y_t+1 -y_t) I Δί (12) 其中 Δί为时刻和 +l时刻之间的采样时间间隔；

云团移动方向定义如下：

( 1 ) 如果 x_i+1 > X,，且 < y,，则有 Θ =

(2) 如果 A₊₁ = ，且 < y_t，则有 θ = πΙ1

(3) 如果 χ_ί+1 < χ,，且 < y_t，则有 θ =

(4) 如果 χ_ί+1 < x,，且 > y_t，则有 θ

(5) 如果且 y_t+l > y_t , 则有 θ = 3π/2；

(6) 如果 χ_ί+1 >χ,， Ky_t+l≥y_t, 则有 θ = πΙ2+ατείαη{ ⁺¹ ) ；

Ι

步骤 1-4-2: 进行云团运动预估；

假设云团在相对较短时间内的运动为刚体运动，位于全天空成像仪视野范围内的云团近似匀速直线运动，通过求取 V ^和 ^平均值预估未来第个云团的移动速度和移动方向 θ_;, 具体表示为：

8. 根据权利要求 1所述的基于地基云图的光伏功率预测方法，其特征在于：所述步骤包括以下步骤：

时刻天顶辐射 H。(0和大气质量 AM (0分别表示为：

H₀ (t) = I_sc (1 + 0.033co5(2^- · dayth 1365))cos0_z (14)

AM (t) = 229 + ( 614cos0_z)² -614cos0_z (15) 其中， /_sc =1367±5W/m² ; 于是：时刻地面理论辐射值表示为：

(t) = H₀(t)/AM (t) (16) 在晴空状况下，令时刻地面理论辐射值等于地面实测辐射值。_fo(0根据式（16) 计算出 C，进而确定时刻地面理论辐射值；

步骤 2-2: 进行太阳遮挡预测；

ίθ, 如果 ( (t)+ex, y_s (t) + ex)点像素值为 0

sisn(t) = i (is

[l, 如果 (t) + (t) + e )点像素值为 1 其中 signW表示时刻太阳的遮挡信息， w'gn(i)=0表示太阳未被云团遮挡， w'gn(i)=l表示太阳被云团遮挡；（x_s(0,y_s(0)表示时刻太阳在图像中的位置， ex的取值为 {-2， -1,0,1,2}；步骤 2-3: 对地面辐照度进行预测；

假设预测阶段的辐射清晰度指数为 pKT, 其由学***均确定，表示为： ρΚΤ = ^ ,ΚΤ^) (19) 其中， N为学***均的权重系数，且满足 Λ<Λ.₊₁; 7^_;)为学习阶段太阳被云团遮挡时刻对应的第个辐射清晰度指数；

于是，预测阶段时刻的地面理论辐照度表示为：

,,、

pI(t) (20)

其中， _P 表示预测阶段 t时刻的地面理论辐照度。

9. 根据权利要求 1所述的基于地基云图的光伏功率预测方法，其特征在于：所述步骤 3 包括以下步骤：

步骤 3-1: 建立基础光电转换模型；

通过所读取的历史功率数据和辐射数据绘制散点图，若 P₀(0、 /(0分别为时刻的实测有功功率和实测地面辐照度，并根据散点图确定曲线类型和曲线拟合方程式中的阶数^ 所述曲线拟合方程式为:

Ρ₀ (0 = a_qI^q (t) + a_q_ ^q-¹ ( + ··· + α₀ (21) 其中 ^A ，…和 _Ω。均为曲线拟合方程式的系数；

根据上述曲线拟合方程式建立目标函数： ιηϊηΫ(α Ι"(ί)

+■■■ + a_(j-P_(j(t)f (22) 通过最小二乘法确定，…和 α。，从而确定时刻的实测有功功率 (0，表示为：

Ρ₀ (0 = f (Ht)) = a_qI^q(t) + a_qJ^q-¹ (0 +… + (23)

其中， Ρ(0表示预测阶段时刻的实测有功功率。

10. 根据权利要求 1所述的基于地基云图的光伏功率预测方法，其特征在于：所述步骤 4 中，通过建立光伏功率预测模型进行光伏功率预测；

所述光伏功率预测模型为： P{t) (26)。