CN107168049A - 光伏阵列输出特性曲线实时获取方法 - Google Patents
光伏阵列输出特性曲线实时获取方法 Download PDFInfo
- Publication number
- CN107168049A CN107168049A CN201710356894.1A CN201710356894A CN107168049A CN 107168049 A CN107168049 A CN 107168049A CN 201710356894 A CN201710356894 A CN 201710356894A CN 107168049 A CN107168049 A CN 107168049A
- Authority
- CN
- China
- Prior art keywords
- photovoltaic array
- photovoltaic
- image
- shade
- array
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 230000005855 radiation Effects 0.000 claims abstract description 25
- 238000001514 detection method Methods 0.000 claims abstract description 22
- 238000005286 illumination Methods 0.000 claims abstract description 18
- 238000004891 communication Methods 0.000 claims abstract description 14
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 13
- 239000000523 sample Substances 0.000 claims abstract description 6
- 238000005094 computer simulation Methods 0.000 claims abstract description 5
- 238000009434 installation Methods 0.000 claims abstract description 4
- 238000012545 processing Methods 0.000 claims abstract description 4
- 230000011218 segmentation Effects 0.000 claims abstract description 4
- 230000001360 synchronised effect Effects 0.000 claims abstract description 4
- 108010074864 Factor XI Proteins 0.000 claims description 6
- 230000001154 acute effect Effects 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000004458 analytical method Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 230000012447 hatching Effects 0.000 claims description 3
- 238000003709 image segmentation Methods 0.000 claims description 3
- 239000004575 stone Substances 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 2
- 238000003708 edge detection Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 4
- 238000005457 optimization Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 241001270131 Agaricus moelleri Species 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/0205—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system
- G05B13/024—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Image Analysis (AREA)
Abstract
一种光伏阵列输出特性曲线实时获取方法，首先综合考虑传感器个数以及辐照度检测效率，确定光伏阵列中辐射传感器的安装位置；然后构建光伏阵列输出特性曲线实时获取系统框架，并将辐照度数据与同步采集的阵列运行图像标记时刻后通过通信网络发送至光伏发电系统控制中心；针对浅色地面与深色光伏板之间的灰度值以及正常光照与遮荫状态下组件的灰度值均存在差别这一特征，提出了基于迭代法阈值分割、LOG算子边缘检测等数字图像处理技术的遮荫组件识别算法，用以识别遮荫组件；最后结合光伏组件参数、阵列结构以及获取的辐照度、识别的遮荫组件，利用计算机仿真绘制局部遮荫时光伏阵列的输出特性曲线。
Description
技术领域
本发明涉及光伏发电系统中光伏阵列遭遇局部阴影时的辐照度测量技术、阴影检测技术以及光伏阵列输出特性曲线实时获取技术，属于光伏发电系统运行条件检测领域。
背景技术
随着传统能源的枯竭与大众环保意识的提高，无污染、分布广且易获取的太阳能日益受到重视。太阳能光伏发电是对太阳能有效利用的主要形式，并成为了继风力发电之后主要的新能源发电方式。近年来光伏组件价格的降低以及政府补贴政策的出台，光伏发电装机容量飞速增长，其中包括大容量的荒漠光伏电站、分散式的中小容量建筑光伏系统。
光伏阵列通常由大量光伏电池以一定的串并联结构组成，以此获取所需的光伏输出电压和功率。为防止出现热斑现象与支路电能倒送现象，光伏阵列还装设有旁路二极管与防逆二极管。采用集中输出控制的光伏阵列受均匀光照时，光伏阵列的输出功率-电压(P-U)特性曲线呈单峰状，传统最大功率点跟踪(MPPT)方法可以很容易跟踪到最大功率点(MPP)；但当阵列受局部遮荫时，局部阴影不仅会削弱光伏阵列潜在的最大功率输出能力，二极管的存在还会导致光伏阵列输出特性复杂化、多峰化，给MPPT控制、重构优化、发电功率预测等带来很大困难。
局部阴影条件下，光伏阵列中各个光伏电池所接收到的太阳光辐照度可能存在差异。考虑到局部阴影的千变万化，理论上可以通过密集布置辐射传感器 的方法，实现对光伏阵列中各个位置辐照度的测量。但大量的传感器意味着建设成本高和系统复杂化，无法实现广泛应用。若对局部阴影形状准确度测量要求越高，传感器的数量将成倍增加。
发明内容
本发明要克服现有技术在局部遮荫条件下光伏阵列输出特性曲线实时获取困难的问题，利用辐射传感器与遮荫组件识别算法获取光伏阵列上的辐照度数据与阴影分布信息，并结合光伏组件参数、阵列结构，利用计算机仿真绘制局部遮荫时光伏阵列的输出特性曲线。本发明所需的传感器数量非常少，但是可以快速、准确地绘制局部阴影条件下光伏阵列的输出特性曲线。
一种光伏阵列输出特性曲线实时获取方法，其特征在于具体实现该方法的步骤如下：
步骤1：综合考虑传感器个数以及辐照度检测效率，确定光伏阵列中辐射传感器的安装位置；具体实现方法如下：
步骤11：对于共有n块光伏组件的光伏发电系统，根据阵列规模，设计m种定点安装辐射传感器的方案，分别记为方案1、方案2、….、方案m，且记方案j(j＝1,2,…,m)中传感器的使用个数为aj；
步骤12：通过排列组合，构造i(i＝1,2,…,n)块光伏组件受遮荫时光伏阵列所有可能的阴影分布类型，记录阴影分布类型数Yi，并计算阴影分布类型总数Y＝ΣYi；
步骤13：采用步骤11中的方案j时，光伏阵列遍历步骤12中总数为Y的阴影分布类型，统计光照区与阴影区的辐照度能同时被检测到的阴影类型数Xj；
步骤14：计算步骤11中方案j的适用率ξj＝Xj/Y；
步骤15：综合考虑传感器个数，定义方案j的辐照度检测效率ρj＝ξj/aj，选取最大辐照度检测效率max(ρj)对应的方案作为最优方案；
步骤16：按照步骤15选取的最优方案在光伏阵列中布置辐射传感器；
步骤2：构建光伏阵列输出特性曲线实时获取系统框架；所述的光伏阵列输出特性曲线实时获取系统框架包括检测光伏阵列接收入射光强度的辐射传感器(1)、采集光伏阵列运行图像的摄像头(2)、传输光伏阵列图像与辐照度数据的通信网络(3)、实现A/D转换与遮荫组件识别并绘制光伏阵列输出特性曲线的光伏发电系统控制中心(4)、光伏阵列(5)，其中辐射传感器(1)、摄像头(2)的输出端与通信网络(3)的输入端连接，通信网络(3)的输出端与光伏发电系统控制中心(4)的输入端相连；摄像头(2)安装在光伏阵列(5)的正前方；辐射传感器(1)按照步骤15所述的辐照度检测效率最高的最优方案来布置；光伏阵列(5)采取包括串-并联结构(SP)、桥型结构(BL)、全连接结构(TCT)在内的所有结构；
步骤3：利用辐射传感器检测光伏阵列中光照区与阴影区的辐照度数据并标记测量时刻，与此同时，利用摄像头同步采集光伏阵列运行图像，将辐照度数据与光伏阵列图像通过通信网络传送到光伏发电系统控制中心；
步骤4：针对浅色地面与深色光伏板之间的灰度值以及正常光照与遮荫状态下组件的灰度值均存在差别这一特征，提出光伏阵列受局部遮荫时的遮荫组件识别算法，识别光伏阵列图像中的遮荫组件；遮荫组件识别算法具体实现方法如下：
步骤41：对光伏阵列图像进行预处理；
步骤411：将图像大小为N1×N2的光伏阵列图像转换为阵列灰度图f(x,y)，其中(x,y)为图像中各像素点的坐标，f(x,y)表示点(x,y)的灰度值；
步骤412：采用迭代法求取阈值K，将灰度值小于等于K的像素点设置为黑，灰度值大于K的像素点设置为白，即将K作为阈值对灰度图f(x,y)进行二值化，得二值图像；
步骤413：去除地面砖块之间的黑色栅线与光伏组件上的白色栅线：设置检测宽度W，逐行遍历步骤412所得的二值图像中横坐标从1+W到N2-W、纵坐标从1到N1范围内所有的像素点(x,y)，若点(x,y)为黑且点(x-W,y)与点(x+W,y)均为白，则认为点(x,y)是黑色栅线点，并将y行中横坐标从x-W到x+W的点设置为白，否则判断下一个像素点，根据同一原理，逐行去除白色栅线点后，再逐列去除图像中的黑、白栅线点；
步骤414：使用中值滤波器过滤去栅线后光伏阵列图像中较多的点状噪声；
步骤42：识别光伏阵列图像中的光伏组件；
步骤421：采用高斯拉普拉斯(LOG)算子检测经过步骤414滤波后的光伏阵列图像中灰度值剧烈变化的边界，产生了n+a+b个亮对象；
步骤422：去除n+a+b个亮对象中由地面石子及阴影产生的a个小连通亮对象；
步骤423：去除经过步骤422处理后剩余的n+b个亮对象中b个非封闭的亮对象：对含有n+b个亮对象的图像使用元素0标记非亮对象区域，元素i(i＝1,2,…,n+b)标记各个亮对象，测量并记录使用元素i标记的亮对象的质心位置Zi与周长Ci，计算Zi在上下左右四个方向与最近亮对象之间的距离，四个距离分别记为Ai、Bi、Ci、Di，计算Li＝Ai+Bi、Wi＝Ci+Di并计算C0i＝2× (Li+Wi)，判断Ci/C0i是否大于0.5且小于1.5，若是则保留该亮对象，若不是则去除该亮对象；
步骤424：对n个亮对象的封闭区域填充白，之外的其他区域填充黑，得组件区域图像I(x,y)，其中白色区域表示组件区域，黑色区域表示背景区域；
步骤43：识别光伏组件中受遮荫的组件；
步骤431：将步骤424所得组件区域图像I(x,y)中的黑色区域记为区域P，并将步骤411所得阵列灰度图f(x,y)中区域P内所有像素点设置为白；
步骤432：使用元素k(k＝1,2,…,n)标记各个组件，计算第k个组件的灰度均值Mk，将Mk作为新的灰度值填充对应的第k个组件区域，得灰度填充图像；
步骤433：计算步骤432的灰度填充图像中n个组件区域的灰度均值M，将M作为阈值对灰度填充图像进行图像分割；
步骤434：采用高斯拉普拉斯算子检测经过步骤433图像分割后的光伏阵列图像中灰度值剧烈变化的边界，并产生了c个亮对象；
步骤435：对c个亮对象的封闭区域填充白，之外的其他区域填充黑，得遮荫组件区域图像G(x,y)，其中白色区域表示遮荫组件区域，黑色区域包括背景区域与受正常光照的组件区域；
步骤5：根据光伏阵列光照区与阴影区的辐照度数据与识别的遮荫组件，结合光伏阵列结构，分析光伏阵列各个组件上的辐照度；
步骤6：结合组件辐照度、光伏组件参数以及光伏阵列结构，利用计算机仿真快速绘制该时刻光伏阵列的输出特性曲线。
本发明的优点是：无需测量各组件或各支路的电流与电压，只需采集少量的辐照度数据与阵列图像，通过简单的遮荫组件识别算法即可绘制光伏阵列的输出特性曲线，不仅有助于降低最大功率点跟踪(MPPT)方法的算法要求，提高MPPT跟踪准确度，而且可以为光伏阵列重构优化提供必要的全局信息，实现光伏发电系统的高效运行。
附图说明
图1是本发明的备选传感器布置方案示意图，其中·表示辐射传感器(1)；点线封闭区域表示光伏阵列(5)；图1a是方案1对角布置的示意图；图1b是方案2三角布置的示意图；图1c是方案3四角布置的示意图；图1d是方案4四边加中心布置的示意图；图1e是方案5六点布置的示意图。
图2是本发明的光伏阵列系统框架图；其中·表示辐射传感器(1)；△表示摄像头(2)；┄表示通信网络(3)；□表示光伏发电系统控制中心(4)；点线封闭区域表示光伏阵列(5)。
图3是本发明的遮荫组件识别算法流程图。
图4是实施案例遮荫组件识别过程的示意图，其中图4a是实施例的光伏阵列灰度图；图4b是组件区域识别结果；图4c是遮荫组件区域识别结果。
图5是图像采集时刻光伏阵列输出特性曲线，其中图5a是I-U特性曲线；图5b是P-U特性曲线。
具体实施方式
下面结合实施例及附图，对本发明作进一步的详细说明，但本发明的实施方式不限于此。实施例采用规模为3×3的光伏发电系统。
步骤一：确定光伏阵列中辐射传感器的安装位置；
a)对于共有9块光伏组件的光伏发电系统，备选的传感器布置方案共5 种：对角布置、三角布置、四角布置、四边加中心布置、六点布置，5种方案分别记为方案1、方案2、方案3、方案4、方案5，传感器布置方案示意图如图1所示，且a1＝2，a2＝3，a3＝4，a4＝5，a5＝6；
b)当1块光伏组件受遮荫时，排列组合后得此时阵列所有可能的阴影分布类型遮荫2块光伏组件时遮荫3块光伏组件时遮荫4块光伏组件时遮荫5块光伏组件时遮荫6块光伏组件时遮荫7块光伏组件时遮荫8块光伏组件时遮荫9块光伏组件时计算阴影分布类型总数Y＝511；
c)采用步骤a)中的方案1时，光伏阵列遍历步骤b)中总数为511的阴影分布类型，统计光照区与阴影区的辐照度能同时被检测到的阴影类型数X1＝145，采用方案2时X2＝385，采用方案3时X3＝449，采用方案4时X4＝481，采用方案5时X5＝497；
d)计算方案1的适用率ξ1＝X1/Y＝0.2838，方案2的适用率ξ2＝X2/Y＝0.7534，方案3的适用率ξ3＝X3/Y＝0.8787，方案4的适用率ξ4＝X4/Y＝0.9413，方案5的适用率ξ5＝X5/Y＝0.9726；
e)综合考虑传感器个数后，计算方案1的辐照度检测效率ρ1＝ξ1/a1＝0.1419，方案2中ρ2＝ξ2/a2＝0.2511，方案3中ρ3＝ξ3/a3＝0.2197，方案4中ρ4＝ξ4/a4＝0.1883，方案5中ρ5＝ξ5/a5＝0.1621，选取最大辐照度检测效率max(ρj)对应的方案2作为最优方案；
f)按照步骤e)选取的最优方案在光伏阵列中布置辐射传感器；
步骤二：构建光伏阵列输出特性曲线实时获取系统框架；实施例使用的光伏阵列输出特性曲线快速实时获取系统框架如图2所示，包括检测光伏阵列接 收入射光强度的辐射传感器(1)、采集光伏阵列运行图像的摄像头(2)、传输光伏阵列图像与辐照度数据的通信网络(3)、实现A/D转换与遮荫组件识别并绘制光伏阵列输出特性曲线的光伏发电系统控制中心(4)、光伏阵列(5)，其中辐射传感器(1)、摄像头(2)的输出端与通信网络(3)的输入端连接，通信网络(3)的输出端与光伏发电系统控制中心(4)的输入端相连，且辐射传感器(1)按照步骤一中辐照度检测效率最高的最优方案来布置，摄像头(2)安装在光伏阵列(5)的正前方，光伏阵列(5)包括3×3块光伏组件，每一块光伏组件并联一个旁路二极管，每3块光伏组件串联成串后，再串上一个防逆二极管，最后将3串组件串并联，形成SP连接结构；
步骤三：利用辐射传感器检测光伏阵列光照区与阴影区的辐照度数据并标记测量时刻，与此同时，利用摄像头同步采集光伏阵列运行图像，将辐照度数据与光伏阵列图像通过通信网络传送到光伏发电系统控制中心，本实施例中局部阴影处辐照度较均匀，辐射传感器测得光照处与阴影处的辐照度分别为1037W/m2和240W/m2；
步骤四：受局部遮荫的光伏组件图像中，浅色地面与深色光伏板之间的灰度值以及正常光照与遮荫状态下组件的灰度值均存在差别，根据这一特征，提出了流程图如图3所示的适用于光伏阵列的遮荫组件识别算法，算法步骤如下：
步骤1：由于摄像头本身的固有特性以及图片拍摄时无法避免的环境干扰，首先对光伏阵列图像进行预处理：
a)将图像大小为685×514的光伏阵列图像转换成如图4a所示的阵列灰度图；
b)采用迭代法求得阈值K为112，将灰度值小于等于112的像素点设置为黑，灰度值大于112的像素点设置为白，即将112作为阈值对图4a 所示的灰度图进行二值化，得二值图像；
c)去除地面砖块之间的黑色栅线与光伏组件上的白色栅线：设置检测宽度W＝3，逐行遍历步骤b)所得的二值图像中横坐标从4到511、纵坐标从1到685范围内所有的像素点(x,y)，若点(x,y)为黑且点(x-3,y)与点(x+3,y)均为白，则认为点(x,y)是黑色栅线点，并将y行中横坐标从x-3到x+3的点设置为白，否则判断下一个像素点，根据同一原理，逐行去除白色栅线点后，再逐列去除图像中的黑、白栅线点；
d)使用中值滤波器过滤去栅线后光伏阵列图像中较多的点状噪声；
步骤2：识别光伏阵列图像中的光伏组件：
a)采用高斯拉普拉斯(LOG)算子检测经过步骤1预处理后的光伏阵列图像中灰度值剧烈变化的边界，并产生了98个亮对象；
b)去除98个亮对象中由地面石子及阴影产生的88个小连通亮对象；
c)去除经过步骤b)处理后剩余的10个亮对象中1个非封闭的亮对象：对含有10个亮对象的图像使用元素0标记非亮对象区域，元素i(i＝1,2,…,10)标记各个亮对象，测量并记录使用元素i标记的亮对象的质心位置Zi与周长Ci，计算Zi在上下左右四个方向与最近亮对象之间的距离，四个距离分别记为Ai、Bi、Ci、Di，计算Li＝Ai+Bi、Wi＝Ci+Di并计算C0i＝2×(Li+Wi)，判断Ci/C0i是否大于0.5且小于1.5，若是则保留该亮对象，若不是则去除该亮对象；
d)对9个亮对象的封闭区域填充白，之外的其他区域填充黑，得如图4b所示的组件区域识别结果，其中白色区域表示组件区域，黑色区域表示背景区域；
步骤3：识别光伏组件中受遮荫的组件：
a)将图4b中的黑色区域记为区域P，并将图4a中区域P内所有像素点设置为白；
b)使用元素k(k＝1,2,…9)标记9个组件，计算第k个组件的灰度均值Mk，将Mk作为新的灰度值填充对应的第k个组件区域，得灰度填充图像；
c)计算步骤b)的灰度填充图像中9个组件区域的灰度均值M＝52，将灰度值52作为阈值对灰度填充图像进行图像分割；
d)采用高斯拉普拉斯算子检测经过步骤c)图像分割后的光伏阵列图像中灰度值剧烈变化的边界，并产生了3个亮对象；
e)对3个亮对象的封闭区域填充白，之外的其他区域填充黑，得如图4c所示的遮荫组件识别结果，其中白色区域表示遮荫组件区域，黑色区域包括背景区域与受正常光照的组件区域；
对比图4a与图4c可见，利用本发明的遮荫组件识别算法可以获得较好的识别效果；
步骤五：根据步骤三获取的光伏阵列光照区与阴影区的辐照度数据与步骤四识别的遮荫组件，结合光伏阵列结构，分析可得光伏阵列各个组件上的辐照度；
步骤六：结合组件辐照度、光伏组件参数以及光伏阵列结构，利用计算机仿真快速绘制该时刻光伏阵列的输出特性曲线如图5所示。
本实施例无需测量各组件或各支路的电流与电压，只需采集少量的辐照度数据与阵列图像，通过简单的遮荫组件识别算法即可绘制光伏阵列的输出特性曲线，不仅有助于光伏发电系统最大输出功率的跟踪与发电量预测，而且可以为光伏阵列重构优化提供必要的全局信息，实现光伏发电系统的高效运行。
Claims (1)
1.一种光伏阵列输出特性曲线实时获取方法，步骤如下：
步骤1：综合考虑传感器个数以及辐照度检测效率，确定光伏阵列中辐射传感器的安装位置；具体实现方法如下：
步骤11：对于共有n块光伏组件的光伏发电系统，根据阵列规模，设计m种定点安装辐射传感器的方案，分别记为方案1、方案2、….、方案m，且记方案j(j＝1,2,…,m)中传感器的使用个数为aj；
步骤12：通过排列组合，构造i(i＝1,2,…,n)块光伏组件受遮荫时光伏阵列所有可能的阴影分布类型，记录阴影分布类型数Yi，并计算阴影分布类型总数Y＝∑Yi；
步骤13：采用步骤11中的方案j时，光伏阵列遍历步骤12中总数为Y的阴影分布类型，统计光照区与阴影区的辐照度能同时被检测到的阴影类型数Xj；
步骤14：计算步骤11中方案j的适用率ξj＝Xj/Y；
步骤15：综合考虑传感器个数，定义方案j的辐照度检测效率ρj＝ξj/aj，选取最大辐照度检测效率max(ρj)对应的方案作为最优方案；
步骤16：按照步骤15选取的最优方案在光伏阵列中布置辐射传感器；
步骤2：构建光伏阵列输出特性曲线实时获取系统框架；所述的光伏阵列输出特性曲线实时获取系统框架包括检测光伏阵列接收入射光强度的辐射传感器(1)、采集光伏阵列运行图像的摄像头(2)、传输光伏阵列图像与辐照度数据的通信网络(3)、实现A/D转换与遮荫组件识别并绘制光伏阵列输出特性曲线的光伏发电系统控制中心(4)、光伏阵列(5)，其中辐射传感器(1)、摄像头(2)的输出端与通信网络(3)的输入端连接，通信网络(3)的输出端与光伏发电系统控制中心(4)的输入端相连；摄像头(2)安装在光伏阵列(5)的正前方；辐射传感器(1)按照步骤15所述的辐照度检测效率最高的最优方案来布置；光伏阵列(5)采取包括串-并联结构(SP)、桥型结构(BL)、全连接结构(TCT)在内的所有结构；
步骤3：利用辐射传感器检测光伏阵列中光照区与阴影区的辐照度数据并标记测量时刻，与此同时，利用摄像头同步采集光伏阵列运行图像，将辐照度数据与光伏阵列图像通过通信网络传送到光伏发电系统控制中心；
步骤4：针对浅色地面与深色光伏板之间的灰度值以及正常光照与遮荫状态下组件的灰度值均存在差别这一特征，提出光伏阵列受局部遮荫时的遮荫组件识别算法，识别光伏阵列图像中的遮荫组件；遮荫组件识别算法具体实现方法如下：
步骤41：对光伏阵列图像进行预处理；
步骤411：将图像大小为N1×N2的光伏阵列图像转换为阵列灰度图f(x,y)，其中(x,y)为图像中各像素点的坐标，f(x,y)表示点(x,y)的灰度值；
步骤412：采用迭代法求取阈值K，将灰度值小于等于K的像素点设置为黑，灰度值大于K的像素点设置为白，即将K作为阈值对灰度图f(x,y)进行二值化，得二值图像；
步骤413：去除地面砖块之间的黑色栅线与光伏组件上的白色栅线：设置检测宽度W，逐行遍历步骤412所得的二值图像中横坐标从1+W到N2-W、纵坐标从1到N1范围内所有的像素点(x,y)，若点(x,y)为黑且点(x-W,y)与点(x+W,y)均为白，则认为点(x,y)是黑色栅线点，并将y行中横坐标从x-W到x+W的点设置为白，否则判断下一个像素点，根据同一原理，逐行去除白色栅线点后，再逐列去除图像中的黑、白栅线点；
步骤414：使用中值滤波器过滤去栅线后光伏阵列图像中较多的点状噪声；
步骤42：识别光伏阵列图像中的光伏组件；
步骤421：采用高斯拉普拉斯(LOG)算子检测经过步骤414滤波后的光伏阵列图像中灰度值剧烈变化的边界，并产生了n+a+b个亮对象；
步骤422：去除n+a+b个亮对象中由地面石子及阴影产生的a个小连通亮对象；
步骤423：去除经过步骤422处理后剩余的n+b个亮对象中b个非封闭的亮对象：对含有n+b个亮对象的图像使用元素0标记非亮对象区域，元素i(i＝1,2,…,n+b)标记各个亮对象，测量并记录使用元素i标记的亮对象的质心位置Zi与周长Ci，计算Zi在上下左右四个方向与最近亮对象之间的距离，四个距离分别记为Ai、Bi、Ci、Di，计算Li＝Ai+Bi、Wi＝Ci+Di并计算C0i＝2×(Li+Wi)，判断Ci/C0i是否大于0.5且小于1.5，若是则保留该亮对象，若不是则去除该亮对象；
步骤424：对n个亮对象的封闭区域填充白，之外的其他区域填充黑，得组件区域图像I(x,y)，其中白色区域表示组件区域，黑色区域表示背景区域；
步骤43：识别光伏组件中受遮荫的组件；
步骤431：将步骤424所得组件区域图像I(x,y)中的黑色区域记为区域P，并将步骤411所得阵列灰度图f(x,y)中区域P内所有像素点设置为白；
步骤432：使用元素k(k＝1,2,…,n)标记各个组件，计算第k个组件的灰度均值Mk，将Mk作为新的灰度值填充对应的第k个组件区域，得灰度填充图像；
步骤433：计算步骤432所得的灰度填充图像中n个组件区域的灰度均值M，将M作为阈值对灰度填充图像进行图像分割；
步骤434：采用高斯拉普拉斯算子检测经过步骤433图像分割后的光伏阵列图像中灰度值剧烈变化的边界，并产生了c个亮对象；
步骤435：对c个亮对象的封闭区域填充白，之外的其他区域填充黑，得遮荫组件区域图像G(x,y)，其中白色区域表示遮荫组件区域，黑色区域包括背景区域与受正常光照的组件区域；
步骤5：根据光伏阵列光照区与阴影区的辐照度数据与识别的遮荫组件，结合光伏阵列结构，分析光伏阵列各个组件上的辐照度；
步骤6：结合组件辐照度、光伏组件参数以及光伏阵列结构，利用计算机仿真快速绘制该时刻光伏阵列的输出特性曲线。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710356894.1A CN107168049B (zh) | 2017-05-19 | 2017-05-19 | 光伏阵列输出特性曲线实时获取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710356894.1A CN107168049B (zh) | 2017-05-19 | 2017-05-19 | 光伏阵列输出特性曲线实时获取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107168049A true CN107168049A (zh) | 2017-09-15 |
| CN107168049B CN107168049B (zh) | 2019-11-29 |
Family
ID=59815759
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710356894.1A Active CN107168049B (zh) | 2017-05-19 | 2017-05-19 | 光伏阵列输出特性曲线实时获取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107168049B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107742171A (zh) * | 2017-10-31 | 2018-02-27 | 浙江工业大学 | 基于移动阴影图像识别的光伏电站发电功率预测方法 |
| CN108008633A (zh) * | 2017-12-11 | 2018-05-08 | 甘肃省电力公司风电技术中心 | 包含多种天气变化的辐照度与光伏组件坐标关联关系建立方法 |
| CN108197371A (zh) * | 2017-12-28 | 2018-06-22 | 中南大学 | 基于bim的固定光伏阵列连线优化 |
| CN108389216A (zh) * | 2018-02-06 | 2018-08-10 | 西安交通大学 | 面向在线铁谱图像磨粒识别的局部自适应阈值分割方法 |
| CN111911841A (zh) * | 2020-08-17 | 2020-11-10 | 长春理工大学 | 一种用于准直式太阳模拟器的辐照均匀度自动调节系统 |
| CN114115418A (zh) * | 2021-11-15 | 2022-03-01 | 华能新能源股份有限公司 | 光伏系统最大功率点递阶跟踪方法与装置 |
| CN117200971A (zh) * | 2023-09-11 | 2023-12-08 | 重庆千信新能源有限公司 | 基于区块链的光伏组件新能源电力数据传输方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102609951A (zh) * | 2012-03-29 | 2012-07-25 | 重庆大学 | 光伏电池板遮荫部分的检测方法 |
| CN103472884A (zh) * | 2013-08-14 | 2013-12-25 | 浙江工业大学 | 一种光伏阵列全局最大功率点跟踪方法 |
| CN103475267A (zh) * | 2013-08-14 | 2013-12-25 | 浙江工业大学 | 一种提高光伏阵列发电效率的方法 |
| CN103712685A (zh) * | 2013-09-29 | 2014-04-09 | 浙江工业大学 | 光伏阵列辐照度测量辨识方法 |
| CN105261021A (zh) * | 2015-10-19 | 2016-01-20 | 浙江宇视科技有限公司 | 去除前景检测结果阴影的方法及装置 |
-
2017
- 2017-05-19 CN CN201710356894.1A patent/CN107168049B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102609951A (zh) * | 2012-03-29 | 2012-07-25 | 重庆大学 | 光伏电池板遮荫部分的检测方法 |
| CN103472884A (zh) * | 2013-08-14 | 2013-12-25 | 浙江工业大学 | 一种光伏阵列全局最大功率点跟踪方法 |
| CN103475267A (zh) * | 2013-08-14 | 2013-12-25 | 浙江工业大学 | 一种提高光伏阵列发电效率的方法 |
| CN103712685A (zh) * | 2013-09-29 | 2014-04-09 | 浙江工业大学 | 光伏阵列辐照度测量辨识方法 |
| CN105261021A (zh) * | 2015-10-19 | 2016-01-20 | 浙江宇视科技有限公司 | 去除前景检测结果阴影的方法及装置 |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107742171A (zh) * | 2017-10-31 | 2018-02-27 | 浙江工业大学 | 基于移动阴影图像识别的光伏电站发电功率预测方法 |
| CN107742171B (zh) * | 2017-10-31 | 2020-08-21 | 浙江工业大学 | 基于移动阴影图像识别的光伏电站发电功率预测方法 |
| CN108008633A (zh) * | 2017-12-11 | 2018-05-08 | 甘肃省电力公司风电技术中心 | 包含多种天气变化的辐照度与光伏组件坐标关联关系建立方法 |
| CN108197371A (zh) * | 2017-12-28 | 2018-06-22 | 中南大学 | 基于bim的固定光伏阵列连线优化 |
| CN108389216A (zh) * | 2018-02-06 | 2018-08-10 | 西安交通大学 | 面向在线铁谱图像磨粒识别的局部自适应阈值分割方法 |
| CN108389216B (zh) * | 2018-02-06 | 2020-06-26 | 西安交通大学 | 面向在线铁谱图像磨粒识别的局部自适应阈值分割方法 |
| CN111911841A (zh) * | 2020-08-17 | 2020-11-10 | 长春理工大学 | 一种用于准直式太阳模拟器的辐照均匀度自动调节系统 |
| CN111911841B (zh) * | 2020-08-17 | 2024-01-19 | 长春理工大学 | 一种用于准直式太阳模拟器的辐照均匀度自动调节系统 |
| CN114115418A (zh) * | 2021-11-15 | 2022-03-01 | 华能新能源股份有限公司 | 光伏系统最大功率点递阶跟踪方法与装置 |
| CN117200971A (zh) * | 2023-09-11 | 2023-12-08 | 重庆千信新能源有限公司 | 基于区块链的光伏组件新能源电力数据传输方法 |
| CN117200971B (zh) * | 2023-09-11 | 2024-03-26 | 重庆千信新能源有限公司 | 基于区块链的光伏组件新能源电力数据传输方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107168049B (zh) | 2019-11-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107168049B (zh) | 光伏阵列输出特性曲线实时获取方法 | |
| Li et al. | Visible defects detection based on UAV‐based inspection in large‐scale photovoltaic systems | |
| CN108347221B (zh) | 一种双面光伏电池板的跟踪控制方法及系统 | |
| Jazayeri et al. | Adaptive photovoltaic array reconfiguration based on real cloud patterns to mitigate effects of non-uniform spatial irradiance profiles | |
| US8076625B2 (en) | Field level tracker controller | |
| CN104272128B (zh) | 太阳光发电监视方法以及在该方法中使用的太阳光发电监视系统 | |
| CN107742171A (zh) | 基于移动阴影图像识别的光伏电站发电功率预测方法 | |
| CN110487397B (zh) | 一种光伏阵列辐照度传感器的布置方法 | |
| CN105915179A (zh) | 一种晶圆及太阳能电池的光致载流子辐射锁相成像检测方法与系统 | |
| CN104166999B (zh) | 一种基于地基云图强度分层的云团提取方法 | |
| CN103712685B (zh) | 光伏阵列辐照度测量辨识方法 | |
| CN108008633A (zh) | 包含多种天气变化的辐照度与光伏组件坐标关联关系建立方法 | |
| CN110533230A (zh) | 一种光伏阵列辐照度传感器的位置优化方法 | |
| CN107133685A (zh) | 光伏发电系统发电量的预测方法及预测系统 | |
| CN110210060A (zh) | 太阳能光伏板表面积灰程度的预测方法 | |
| CN106897998A (zh) | 太阳能直射辐射强度信息预测方法和系统 | |
| CN112446863A (zh) | 基于图像处理的光伏阵列热斑检测方法和检测系统 | |
| Guerriero et al. | Automatic edge identification for accurate analysis of thermographic images of solar panels | |
| CN111860592A (zh) | 一种少样本条件下太阳能电池片缺陷分类检测方法 | |
| CN116800172B (zh) | 一种光电转换方法及储能系统 | |
| CN112182714A (zh) | 顾及飞行员视线和天气条件的建筑太阳能潜力计算方法 | |
| CN113160311A (zh) | 太阳能电池板调节方法及调节装置、追日系统和存储设备 | |
| CN116192005A (zh) | 基于微时空尺度辐照预测的光伏跟踪方法、系统及介质 | |
| Ye et al. | Research on PV array output characteristics based on shadow image recognition | |
| CN111596381B (zh) | 一种双辐照计推测直射辐照比例的方法 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20231102 Address after: 310000 room 506, building C, No. 3, Weiye Road, Puyan street, Binjiang District, Hangzhou City, Zhejiang Province Patentee after: HANGZHOU YUNMU TECHNOLOGY CO.,LTD. Address before: 310014 science and Technology Office of Zhejiang University of technology, No. 18 Chaowang Road, Xiacheng District, Hangzhou, Zhejiang Patentee before: JIANG University OF TECHNOLOGY |