CN109060618A - 一种光伏发电功率预测装置及方法 - Google Patents
一种光伏发电功率预测装置及方法 Download PDFInfo
- Publication number
- CN109060618A CN109060618A CN201810946791.5A CN201810946791A CN109060618A CN 109060618 A CN109060618 A CN 109060618A CN 201810946791 A CN201810946791 A CN 201810946791A CN 109060618 A CN109060618 A CN 109060618A
- Authority
- CN
- China
- Prior art keywords
- power generation
- wave band
- photovoltaic
- aod
- prediction
- 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.)
- Pending
Links
- 238000010248 power generation Methods 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 title claims abstract description 51
- 230000005855 radiation Effects 0.000 claims abstract description 68
- 238000013528 artificial neural network Methods 0.000 claims abstract description 32
- 238000004891 communication Methods 0.000 claims abstract description 19
- 230000005540 biological transmission Effects 0.000 claims abstract description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 17
- 238000006243 chemical reaction Methods 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims description 14
- 239000000443 aerosol Substances 0.000 claims description 13
- 210000002569 neuron Anatomy 0.000 claims description 11
- 238000012546 transfer Methods 0.000 claims description 9
- 238000012549 training Methods 0.000 claims description 7
- 238000012360 testing method Methods 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 3
- 239000003673 groundwater Substances 0.000 claims description 3
- 230000001360 synchronised effect Effects 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims 1
- 210000005036 nerve Anatomy 0.000 claims 1
- 238000013277 forecasting method Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 102000003712 Complement factor B Human genes 0.000 description 2
- 108090000056 Complement factor B Proteins 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 241001365789 Oenanthe crocata Species 0.000 description 1
- 239000005427 atmospheric aerosol Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 239000008187 granular material Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000003595 mist Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 229910021420 polycrystalline silicon Inorganic materials 0.000 description 1
- 229920005591 polysilicon Polymers 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/06—Investigating concentration of particle suspensions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02S—GENERATION OF ELECTRIC POWER BY CONVERSION OF INFRARED RADIATION, VISIBLE LIGHT OR ULTRAVIOLET LIGHT, e.g. USING PHOTOVOLTAIC [PV] MODULES
- H02S50/00—Monitoring or testing of PV systems, e.g. load balancing or fault identification
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/50—Photovoltaic [PV] energy
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Evolutionary Computation (AREA)
- Computational Linguistics (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- Artificial Intelligence (AREA)
- Dispersion Chemistry (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Photovoltaic Devices (AREA)
Abstract
本发明公开了一种光伏发电功率预测装置及方法,本发明涉及光伏发电技术领域,所述光伏发电功率预测装置包含PM10传感器、PM2.5传感器、湿度传感器、能见度传感器、温度传感器、RS232通信接口、单片机、电源电路、时钟电路、3G通信模块、FPGA芯片、FLASH芯片;所述功率预测方法在获取大量PM2.5浓度‑PM10浓度‑湿度‑能见度‑AOD的实际样本的基础上,利用BP神经网络方法对AOD进行预测,选取AOD与可降水量为主要变化参数,采用双波段大气辐射传输模型计算太阳直接辐射强度与散射辐射强度,进而预测出光伏电池板接收的总辐射强度,最终通过光电转换模型实现对光伏发电功率的准确、实时预测,预测光伏发电的输出功率,具有计算过程清晰,预测精度高的优点。
Description
技术领域
本发明涉及光伏发电技术领域,具体为一种光伏发电功率预测装置及方法。
背景技术
光伏发电功率预测对于电力系统的调度有着重要意义。已有的光伏发电预测方法可分为统计预测方法以及基于大气辐射传输模型的物理预测方法。相对于统计预测方法,物理预测方法具有计算简单、不依赖于历史数据的优点,但其预测精度受模型参数影响较大。
大气气溶胶光学厚度(Aerosol Optical Depth,AOD)是大气辐射传输模型中最重要的参数,用于描述气溶胶对光的衰减作用,其获取方式有卫星观测方式与地基观测方式。卫星观测方式是基于遥感卫星搭载的中分辨率成像光谱仪(MODIS)监测数据对AOD进行反演,受卫星过境时间与数据发布周期影响,反演的实时性较差;地基观测是通过在地面观测站点利用光度计测量数据对AOD进行反演,反演精度较高,但是光度计的标定工作颇为复杂,且观测站点较少;因此,如何实时、准确而又便捷的获取AOD,成为光伏发电物理预测方法的关键问题。
近年来,我国出现了大范围的雾霾天气。雾和霾分别指悬浮于大气中的细微水滴和灰尘颗粒,是大气气溶胶的重要组成部分。2012年,国家环境空气监测网开始实时发布空气质量指数(AirQualityIndex,AQI)以及其中的PM10、PM2.5等6项基本监测数据,至2015年1月2日,已经覆盖367个城市。2016年1月1日起,环保部正式向社会发布空气质量预测预报信息,主要内容包括重点城市未来24小时、48小时空气质量预报,城市AQI范围、空气质量级别及首要污染物等。
受上述结论启发,鉴于大气气溶胶光学厚度(AOD)与PM10浓度、PM2.5浓度等数据具有较高的相关性,且空气质量监测站的站点数量将不断增加,本发明考虑建立AOD的BP神经网络预测模型,进而提出一种光伏发电功率预测装置及方法,以更好的进行光伏发电功率的预测。
发明内容
(一)解决的技术问题
针对现有技术的不足,本发明提供一种光伏发电功率预测装置及方法,以对光伏发电功率进行实时、准确地预测。
(二)技术方案
为实现以上目的,本发明提供一种光伏发电功率预测装置,包含PM10传感器、PM2.5传感器、湿度传感器、能见度传感器、温度传感器、RS232通信接口、单片机、电源电路、时钟电路、3G通信模块、FPGA芯片、FLASH芯片。
优选的,所述PM10传感器、PM2.5传感器、湿度传感器、能见度传感器、温度传感器分别进行PM10浓度、PM2.5浓度、湿度、能见度、气温的采集,并与RS232通信接口的一端连接,所述RS232通信接口另一端与单片机连接,其作用为将传感器采集的数据发送至单片机。
优选的,所述电源电路与单片机及FPGA芯片连接,用于提供电源,优选的,所述时钟电路与单片机和FPGA连接,用于提供同步的时钟信号。
优选的,所述单片机的15个端口直接与FPGA的15个端口相连,采用并行的方式传输数据,其中端口1用作发送使能端,端口2用作读入使能端,端口3用作读入时钟端,其余12个端口为数据口。
优选的,所述FLASH芯片与FPGA连接,用于存放BP神经网络预测模型的配置文件,所述FPGA芯片与单片机相连接,用于进行BP神经网络预测模型的计算,所述3G通信模块与单片机连接,用于光伏发电功率预测结果的远传,所述LED指示灯与单片机连接,用于指示工作状态。
本发明还提供了基于上述光伏发电功率预测装置的预测方法,所述预测方法在获取大量PM2.5浓度-PM10浓度-湿度-能见度-AOD的实际样本的基础上,采用BP神经网络方法对AOD进行预测,选取AOD与可降水量为主要变化参数,采用双波段大气辐射传输模型计算太阳直接辐射强度与散射辐射强度,进而预测出光伏电池板接收的总辐射强度,最终通过光电转换模型实现对光伏发电功率的准确、实时预测。
优选的,所述预测方法包括以下具体步骤进行:
步骤1.以能见度v,湿度h,PM2.5浓度C0,PM10浓度C1为输入,以440nm波长AODτ440与1020nm波长AODτ1020为输出,建立并训练BP神经网络预测模型;
步骤2.将BP神经网络预测模型所预测的440nm波长AOD与1020nm波长分别代表波段1的AODτa1与波段2的AODτa2,采用线性拟合模型计算可降水量w,进而利用双波段大气辐射传输模型计算太阳直接辐射强度与散射辐射强度,从而预测出光伏电池板所接收的总辐射强度;
步骤3.基于光电转换模型,预测光伏发电的输出功率。
优选的,所述步骤1按以下方法进行处理:
a.采集数据构造样本集;
记录某一时刻k下PM2.5浓度C0、PM10浓度C1、能见度v、湿度h和440nm波长的AODτ440与1020nm波长的AODτ1020,得到一个样本(C0(k),C1(k),v(k),h(k),τ440(k),τ1020(k)),通过收集多个不同时刻下的样本,形成样本集{(C0(k),C1(k),v(k),h(k),τ440(k),τ1020(k))};
b.建立三层结构BP神经网络预测模型,其中,输入层神经元节点数为4个,隐含层神经元节点数为5个,输出层神经元节点数为2个,隐含层神经元转移函数使用双曲正切函数,输出层神经元转移函数使用S型函数;
c.随机抽取样本集中的一部分做为训练样本,剩余部分做为检验样本,对所建立的BP神经网络预测模型进行训练,对于第k个时刻所构造的样本(C0(k),C1(k),v(k),h(k),τ440(k),τ1020(k)),取BP神经网络预测模型的第一个输入为v(k),第二个输入为h(k),第三个输入为C0(k),第四个输入为C1(k),第一个输出为τ440(k),第二个输出为τ1020(k);
优选的,所述步骤2按以下方法进行处理:
利用所述光伏发电功率预测装置采集某一环境下PM2.5浓度C0、PM10浓度C1、能见度v、湿度h,作为BP神经网络预测模型的输入,则预测模型的第一个输出为当前环境下440nm波长AODτ440的预测值第二个输出为1020nm波长AODτ1020的预测值
采用线性拟合模型计算可降水量w,即
w=-0.518+0.177·SVP (1)
其中,SVP为地面水汽压,可根据空气相对湿度与环境温度求得:
其中,R为空气相对湿度,es0为0℃时的饱和水汽压,es0=6.1078,经验系数a=7.5,b=237.3,t为环境温度;
所述双波段大气辐射传输模型,其直射辐射强度的计算公式为:
式中,Ebn为总直射辐射强度,i=1、2,Ebni为波段i直射辐射强度,E0ni为波段i大气层外直射辐射强度,且有E0n1=635.4W/m2,E0n2=709.7W/m2,Tai为波段i气溶胶传输衰减系数,计算公式为:
Tai=exp(-maτai) (4)
其中,取且有
ma=[cosZ+0.16851Z0.18198/(95.318-Z)1.9542]-1 (5)
Twi为波段i可降水量衰减系数,计算公式为:
Tw1=(1+h1mw)/(1+h2mw) (6)
TRi为波段i雷利系数,各波段计算公式为
TR2=(1-0.010394m'R)/(1-0.00011042m'R 2) (9)
以上式中参数h1,h2,mw,c1,c2,c3,c4,m'R的取值本专利不作过多说明;
所述双波段大气辐射传输模型,其散射辐射强度的计算公式为:
式中,Ed为总散射辐射强度,i=1、2,Edi为波段i散射辐射强度,Eddi为波段i反向散射因子,Edpi为波段i入射散射辐射因子,Ba为气溶胶前向散射因子,BRi为波段i雷利衰减的前向散射因子,Fi为波段i修正系数,Tasi为波段i气溶胶散射系数,ρgi为波段i地表反照率,ρsi为波段i天空反照率,以上参数本专利不作过多说明;
对于朝向赤道,与水平地面呈β倾角布置的光伏电池板,考虑其面板法线方向与太阳直射辐射方向的夹角,并假定散射辐射呈各向同性,则其接收的总辐射强度为:
其中,Z为太阳天顶角,θh为太阳时角。
优选的,所述步骤3按以下方法进行处理:
所述光电转换模型为:
Ps=ηpvSEx[1-0.005(t0+25)] (12)
式中:Ps为光伏电池板输出功率,ηpv为光伏电池板转换效率,S为光伏电池板的面积,t0为光伏电池板工作温度,由太阳辐射强度EX与气温T进行估算:
t0=T+kEx (13)
其中,气温T由所述光伏发电功率预测装置中温度传感器进行采集,k为温度系数。
(三)有益效果
本发明提供的一种光伏发电功率预测装置及方法,具备以下有益效果:
(1)该光伏发电功率预测装置及方法,利用可实时获取的PM2.5浓度C0、PM10浓度C1、能见度v、湿度h等数据对AOD进行直接预测,与现有的卫星观测法与地基观测法相比,省去了其中的复杂反演过程,并且能保证预测的实时性与准确性。
(2)该光伏发电功率预测装置及方法,基于BP神经网络预测模型所预测的440nm波长AOD与1020nm波长AOD,以及采用线性拟合模型计算的可降水量w,利用双波段大气辐射传输模型计算太阳直接辐射强度与散射辐射强度,最终基于光电转换模型,预测光伏发电的输出功率,具有计算过程清晰,预测精度高的优点。
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
附图说明
图1为本发明光伏发电功率预测装置的总体结构框图;
图2为本发明光伏发电功率预测装置的FPGA与单片机并行连接示意图;
图3为本发明光伏发电功率预测装置的3G通信模块的结构框图;
图4为本发明光伏发电功率预测方法中BP神经网络AOD预测模型;
图5为本发明实施例中AOD预测模型在训练样本上的预测结果;
图6为本发明实施例中AOD预测模型在检验样本上的预测误差;
图7为本发明实施例中光伏发电功率的预测结果。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为实现以上目的,本发明提供一种光伏发电功率预测装置,参阅图1,包含PM10传感器、PM2.5传感器、湿度传感器、能见度传感器、温度传感器、RS232通信接口、单片机、电源电路、时钟电路、3G通信模块、FPGA芯片、FLASH芯片。实施例中,PM2.5传感器选用DS8001-232传感器模块、PM10传感器选用T-PM2510泵吸式激光传感器、湿度传感器选用MHTRS232传感器、能见度传感器选用NWD-10传感器、温度传感器选用IAQ50传感器,单片机选用dsPIC33FJ256GP510单片机,FPGA芯片选用Altera Cyclone IV-EP4CE115F29芯片,3G通信模块选用华为EM770W模块,FLASH芯片选用AT24C02芯片。
PM10传感器、PM2.5传感器、湿度传感器、能见度传感器、温度传感器与RS232通信接口的一端连接,RS232通信接口另一端与单片机连接,用于采集传感器数据。
电源电路与单片机及FPGA芯片连接,用于提供电源,时钟电路与单片机和FPGA连接,用于提供同步的时钟信号。
单片机的15个端口直接与FPGA的15个端口相连,采用并行的方式传输数据,参阅图2,其中端口1用作发送使能端,端口2用作读入使能端,端口3用作读入时钟端,其余12个端口为数据口。
FLASH芯片与FPGA连接,用于存放BP神经网络预测模型的配置文件,FPGA芯片与单片机相连接,用于BP神经网络预测模型的计算。
3G通信模块与单片机连接,用于光伏发电功率预测结果的远传。参阅图3,本发明实施例中的3G通信模块由华为EM770W模块及其外围电路组成,由于EM770W模块为5V电平,而dsPIC33FJ256GP510单片机I/O接口为3.3V电平,故采用74LVC4245芯片进行双向电平转换。
LED指示灯与单片机连接,用于指示工作状态。
本发明所提供的基于上述光伏发电功率预测装置的预测方法,包括以下具体步骤进行:
步骤1.以能见度v,湿度h,PM2.5浓度C0,PM10浓度C1为输入,以440nm波长AODτ440与1020nm波长AODτ1020为输出,建立并训练BP神经网络预测模型;
步骤2.将BP神经网络预测模型所预测的440nm波长AOD与1020nm波长分别代表波段1的AODτa1与波段2的AODτa2,采用线性拟合模型计算可降水量w,进而利用双波段大气辐射传输模型计算太阳直接辐射强度与散射辐射强度,从而预测出光伏电池板所接收的总辐射强度;
步骤3.基于光电转换模型,预测光伏发电的输出功率。
步骤1按以下方法进行处理:
a.采集数据构造样本集;
记录某一时刻k下PM2.5浓度C0、PM10浓度C1、能见度v、湿度h和440nm波长的AODτ440与1020nm波长的AODτ1020,得到一个样本(C0(k),C1(k),v(k),h(k),τ440(k),τ1020(k)),通过收集多个不同时刻下的样本,形成样本集{(C0(k),C1(k),v(k),h(k),τ440(k),τ1020(k))},实施例中,选取北京市进行数据收集,自2015年1月1日至2016年11月30日,共收集样本2141组;
b.建立三层结构BP神经网络预测模型,参阅图4,其中,输入层神经元节点数为4个,隐含层神经元节点数为5个,输出层神经元节点数为2个,隐含层神经元转移函数使用双曲正切函数,输出层神经元转移函数使用S型函数;
c.随机抽取样本集中的一部分做为训练样本,剩余部分做为检验样本,对所建立的BP神经网络预测模型进行训练,对于第k个时刻所构造的样本(C0(k),C1(k),v(k),h(k),τ440(k),τ1020(k)),取BP神经网络预测模型的第一个输入为v(k),第二个输入为h(k),第三个输入为C0(k),第四个输入为C1(k),第一个输出为τ440(k),第二个输出为τ1020(k),实施例中,随机抽取2141组样本中的1971组做为训练样本,其余170组作为检验样本;
实施例中,训练后的BP神经网络预测模型的权值、阈值矩阵如下表1所示,以配置文件的形式存放于光伏发电功率预测装置的Flash芯片中。
表1
实施例中,BP神经网络预测模型在训练样本上的预测结果参阅图5,在检验样本上的预测误差参阅图6,可知,该预测模型具有较高的预测精度。
步骤2按以下方法进行处理:
利用光伏发电功率预测装置采集某一环境下PM2.5浓度C0、PM10浓度C1、能见度v、湿度h,作为BP神经网络预测模型的输入,则预测模型的第一个输出为当前环境下440nm波长AODτ440的预测值第二个输出为1020nm波长AODτ1020的预测值
采用线性拟合模型计算可降水量w,即
w=-0.518+0.177·SVP (1)
其中,SVP为地面水汽压,可根据空气相对湿度与环境温度求得:
其中,R为空气相对湿度,es0为0℃时的饱和水汽压,es0=6.1078,经验系数a=7.5,b=237.3,t为环境温度;
双波段大气辐射传输模型,其直射辐射强度的计算公式为:
式中,Ebn为总直射辐射强度,i=1、2,Ebni为波段i直射辐射强度,E0ni为波段i大气层外直射辐射强度,且有E0n1=635.4W/m2,E0n2=709.7W/m2,Tai为波段i气溶胶传输衰减系数,计算公式为:
Tai=exp(-maτai) (4)
其中,取且有
ma=[cosZ+0.16851Z0.18198/(95.318-Z)1.9542]-1 (5)
Twi为波段i可降水量衰减系数,计算公式为:
Tw1=(1+h1mw)/(1+h2mw) (6)
TRi为波段i雷利系数,各波段计算公式为
TR2=(1-0.010394m'R)/(1-0.00011042m'R 2) (9)
实施例中,以上式中参数h1,h2,mw,c1,c2,c3,c4,m'R的取值由太阳天顶角Z或者可降水量w进行拟合:
mw=[cosZ+0.10648Z0.11423/(93.781-Z)1.9203]-1
h1=w(0.065445+0.00029901w)/(1+1.2728w)
h2=w(0.065687+0.0013218w))/(1+1.2008w)
c1=w(19.566-1.6506w+1.0672w2)/(1+5.4248w+1.6005w2)
c2=w(0.50158-0.14732w+0.047584w2)/(1+1.1811w+1.0699w2)
c3=w(21.286-0.39232w+1.2692w2)/(1+4.8318w+1.412w2)
c4=w(0.70992-0.23155w+0.096514w2)/(1+0.44907w+0.75425w2)
m'R=[cosZ+0.48353Z0.095846/(96.741-Z)1.754]-1
双波段大气辐射传输模型,其散射辐射强度的计算公式为:
式中,Ed为总散射辐射强度,i=1、2,Edi为波段i散射辐射强度,Eddi为波段i反向散射因子,Edpi为波段i入射散射辐射因子,Ba为气溶胶前向散射因子,BRi为波段i雷利衰减的前向散射因子,Fi为波段i修正系数,Tasi为波段i气溶胶散射系数,ρgi为波段i地表反照率,ρsi为波段i天空反照率,实施例中,Ba的取值由太阳天顶角Z进行拟合:
Ba=1-exp(-0.6931-1.8326cosZ)
Tasi=exp(-maωiτai)
其中,取ω1=0.92,ω2=0.84,且有
实施例中,对于波段1,雷利衰减的前向散射因子BR1的取值由下式进行计算:
其中mR由太阳天顶角进行拟合:
mR=[cosZ+0.48353Z0.095846/(96.741-Z)1.754]-1
实施例中,对于波段2,取雷利衰减的前向散射因子BR2=0.5;
实施例中,取F1=F2=1.2,ρg1=ρg2=0.2,ρs1=ρs2=0.1。
对于朝向赤道,与水平地面呈β倾角布置的光伏电池板,考虑其面板法线方向与太阳直射辐射方向的夹角,并假定散射辐射呈各向同性,则其接收的总辐射强度为:
其中,θh为太阳时角。
步骤3按以下方法进行处理:
光电转换模型为:
Ps=ηpvSEx[1-0.005(t0+25)] (12)
式中:Ps为光伏电池板输出功率,ηpv为光伏电池板转换效率,S为光伏电池板的面积,t0为光伏电池板工作温度,由太阳辐射强度EX与气温T进行估算:
t0=T+kEx (13)
其中,气温T由光伏发电功率预测装置中温度传感器进行采集,k为温度系数,实施例中取典型值0.03。
对于上述功率预测方法,通过将步骤2及步骤3写入光伏发电功率预测装置中的单片机,并且由FPGA芯片完成BP神经网络预测模型的计算,即可实现光伏发电功率的预测。
为了验证本发明所提功率预测方法的功效,利用保定市华北电力大学一校区自动化楼楼顶平台布置的1.2kW多晶硅光伏发电系统对开展功率预测实验。该系统的光伏电池板朝向赤道方向布置,与地面倾角呈35°。选取晴朗无云的2016年12月14日(空气质量为优)、12月16日(空气质量为重度污染),对每日9时至15时各整点时刻的光伏发电功率进行滚动预测,并与实测功率数据进行对比,参阅图7。计算14日与16日的光伏发电功率预测的平均相对误差,为12.8%,说明所提光伏发电功率预测方法具有较高的精度。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下。由语句“包括一个......限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素”。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (10)
1.一种光伏发电功率预测装置,其特征在于:包括以下元件:PM10传感器、PM2.5传感器、湿度传感器、能见度传感器、温度传感器、RS232通信接口、单片机、电源电路、时钟电路、3G通信模块、FPGA芯片、FLASH芯片。
2.根据权利要求1所述的一种光伏发电功率预测装置,其特征在于:所述PM10传感器、PM2.5传感器、湿度传感器、能见度传感器、温度传感器分别进行PM10浓度、PM2.5浓度、湿度、能见度、气温的采集,并与RS232通信接口的一端连接,所述RS232通信接口另一端与单片机连接,其作用为将传感器采集的数据发送至单片机。
3.根据权利要求2所述的一种光伏发电功率预测装置,其特征在于:所述所述电源电路与单片机及FPGA芯片连接,用于提供电源,优选的,所述时钟电路与单片机和FPGA连接,用于提供同步的时钟信号。
4.根据权利要求3所述的一种光伏发电功率预测装置,其特征在于:所述单片机的15个端口直接与FPGA的15个端口相连,采用并行的方式传输数据,其中端口1用作发送使能端,端口2用作读入使能端,端口3用作读入时钟端,其余12个端口为数据口。
5.根据权利要求4所述的一种光伏发电功率预测装置,其特征在于:所述FLASH芯片与FPGA连接,用于存放BP神经网络预测模型的配置文件,所述FPGA芯片与单片机相连接,用于进行BP神经网络预测模型的计算,所述3G通信模块与单片机连接,用于光伏发电功率预测结果的远传,所述LED指示灯与单片机连接,用于指示工作状态。
6.一种光伏发电功率预测装置的预测方法,其特征在于:包括以下具体步骤进行:
步骤1.以能见度v,湿度h,PM2.5浓度C0,PM10浓度C1为输入,以440nm波长AODτ440与1020nm波长AODτ1020为输出,建立并训练BP神经网络预测模型;
步骤2.将BP神经网络预测模型所预测的440nm波长与1020nm波长分别代表波段1的AODτa1与波段2的AODτa2,采用线性拟合模型计算可降水量w,进而利用双波段大气辐射传输模型计算太阳直接辐射强度与散射辐射强度,从而预测出光伏电池板所接收的总辐射强度;
步骤3.基于光电转换模型,预测光伏发电的输出功率。
7.根据权利要求6所述的一种光伏发电功率预测装置的预测方法,其特征在于:所述步骤1按以下方法进行处理:
a.采集数据构造样本集;
记录某一时刻k下PM2.5浓度C0、PM10浓度C1、能见度v、湿度h和440nm波长的AODτ440与1020nm波长的AODτ1020,得到一个样本(C0(k),C1(k),v(k),h(k),τ440(k),τ1020(k)),通过收集多个不同时刻下的样本,形成样本集{(C0(k),C1(k),v(k),h(k),τ440(k),τ1020(k))};
b.建立三层结构BP神经网络预测模型,其中,输入层神经元节点数为4个,隐含层神经元节点数为5个,输出层神经元节点数为2个,隐含层神经元转移函数使用双曲正切函数,输出层神经元转移函数使用S型函数;
c.随机抽取样本集中的一部分做为训练样本,剩余部分做为检验样本,对所建立的BP神经网络预测模型进行训练,对于第k个时刻所构造的样本(C0(k),C1(k),v(k),h(k),τ440(k),τ1020(k)),取BP神经网络预测模型的第一个输入为v(k),第二个输入为h(k),第三个输入为C0(k),第四个输入为C1(k),第一个输出为τ440(k),第二个输出为τ1020(k)。
8.根据权利要求6所述的一种光伏发电功率预测装置的预测方法,其特征在于:所述步骤2按以下方法进行处理:
利用所述光伏发电功率预测装置采集某一环境下PM2.5浓度C0、PM10浓度C1、能见度v、湿度h,作为BP神经网络预测模型的输入,则预测模型的第一个输出为当前环境下440nm波长AODτ440的预测值第二个输出为1020nm波长AODτ1020的预测值
采用线性拟合模型计算可降水量w,即
w=-0.518+0.177·SVP (1)
其中,SVP为地面水汽压,可根据空气相对湿度与环境温度求得:
其中,R为空气相对湿度,es0为0℃时的饱和水汽压,es0=6.1078,经验系数a=7.5,b=237.3,t为环境温度;
所述双波段大气辐射传输模型,其直射辐射强度的计算公式为:
式中,Ebn为总直射辐射强度,i=1、2,Ebni为波段i直射辐射强度,E0ni为波段i大气层外直射辐射强度,且有E0n1=635.4W/m2,E0n2=709.7W/m2,Tai为波段i气溶胶传输衰减系数,计算公式为:
Tai=exp(-maτai) (4)
其中,取且有
ma=[cosZ+0.16851Z0.18198/(95.318-Z)1.9542]-1
(5)
Twi为波段i可降水量衰减系数,计算公式为:
Tw1=(1+h1mw)/(1+h2mw) (6)
TRi为波段i雷利系数,各波段计算公式为
以上式中参数h1,h2,mw,c1,c2,c3,c4,m'R的取值本专利不作过多说明;
所述双波段大气辐射传输模型,其散射辐射强度的计算公式为:
式中,Ed为总散射辐射强度,i=1、2,Edi为波段i散射辐射强度,Eddi为波段i反向散射因子,Edpi为波段i入射散射辐射因子,Ba为气溶胶前向散射因子,BRi为波段i雷利衰减的前向散射因子,Fi为波段i修正系数,Tasi为波段i气溶胶散射系数,ρgi为波段i地表反照率,ρsi为波段i天空反照率,以上参数本专利不作过多说明;
对于朝向赤道,与水平地面呈β倾角布置的光伏电池板,考虑其面板法线方向与太阳直射辐射方向的夹角,并假定散射辐射呈各向同性,则其接收的总辐射强度为:
其中,Z为太阳天顶角,θh为太阳时角。
9.根据权利要求6所述的一种光伏发电功率预测装置的预测方法,其特征在于:所述步骤3按以下方法进行处理:
所述光电转换模型为:
Ps=ηpvSEx[1-0.005(t0+25)] (12)
式中:Ps为光伏电池板输出功率,ηpv为光伏电池板转换效率,S为光伏电池板的面积,t0为光伏电池板工作温度,由太阳辐射强度EX与气温T进行估算:
t0=T+kEx (13)
其中,气温T由所述光伏发电功率预测装置中温度传感器进行采集,k为温度系数。
10.根据权利要求6所述的一种光伏发电功率预测装置的预测方法,其特征在于:所述预测方法在获取大量PM2.5浓度-PM10浓度-湿度-能见度-AOD的实际样本的基础上,采用BP神经网络方法对AOD进行预测,选取AOD与可降水量为主要变化参数,采用双波段大气辐射传输模型计算太阳直接辐射强度与散射辐射强度,进而预测出光伏电池板接收的总辐射强度,最终通过光电转换模型实现对光伏发电功率的准确、实时预测。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810946791.5A CN109060618A (zh) | 2018-08-20 | 2018-08-20 | 一种光伏发电功率预测装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810946791.5A CN109060618A (zh) | 2018-08-20 | 2018-08-20 | 一种光伏发电功率预测装置及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109060618A true CN109060618A (zh) | 2018-12-21 |
Family
ID=64686431
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810946791.5A Pending CN109060618A (zh) | 2018-08-20 | 2018-08-20 | 一种光伏发电功率预测装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109060618A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109829572A (zh) * | 2019-01-14 | 2019-05-31 | 国网江苏省电力有限公司苏州供电分公司 | 雷电气候下的光伏发电功率预测方法 |
CN111769904A (zh) * | 2020-06-23 | 2020-10-13 | 电子科技大学 | 一种反向散射通信系统中多反射设备并行传输的检测方法 |
CN111861784A (zh) * | 2020-06-01 | 2020-10-30 | 大唐东北电力试验研究院有限公司 | 基于人工蜂群优化神经网络的光伏发电功率短期预测方法 |
CN117740816A (zh) * | 2023-11-10 | 2024-03-22 | 南京旗云中天科技有限公司 | 提高光伏组件积灰污染比检测精度的方法、传感器、系统 |
CN117879490A (zh) * | 2024-03-13 | 2024-04-12 | 国网山东省电力公司微山县供电公司 | 一种光伏设备发电功率预测方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011158363A1 (ja) * | 2010-06-17 | 2011-12-22 | 株式会社四国総合研究所 | 日射強度予測システムと太陽光発電出力予測システム |
CN103500365A (zh) * | 2013-09-18 | 2014-01-08 | 广州供电局有限公司 | 光伏发电功率预测方法和系统 |
CN105574619A (zh) * | 2016-01-06 | 2016-05-11 | 国家电网公司 | 一种太阳能光伏发电出力预测系统及预测方法 |
CN105631558A (zh) * | 2016-03-22 | 2016-06-01 | 国家电网公司 | 基于相似日的bp神经网络光伏发电系统功率预测方法 |
CN106096246A (zh) * | 2016-06-06 | 2016-11-09 | 华北电力大学(保定) | 基于pm2.5和pm10的气溶胶光学厚度估计方法 |
CN106372718A (zh) * | 2016-08-29 | 2017-02-01 | 华北电力大学(保定) | 雾霾积灰条件下的光伏发电功率输出减少率估计方法 |
CN207036621U (zh) * | 2017-07-07 | 2018-02-23 | 宁夏大学 | 一种光伏电池表面沙尘沉积量测量装置 |
-
2018
- 2018-08-20 CN CN201810946791.5A patent/CN109060618A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011158363A1 (ja) * | 2010-06-17 | 2011-12-22 | 株式会社四国総合研究所 | 日射強度予測システムと太陽光発電出力予測システム |
CN103500365A (zh) * | 2013-09-18 | 2014-01-08 | 广州供电局有限公司 | 光伏发电功率预测方法和系统 |
CN105574619A (zh) * | 2016-01-06 | 2016-05-11 | 国家电网公司 | 一种太阳能光伏发电出力预测系统及预测方法 |
CN105631558A (zh) * | 2016-03-22 | 2016-06-01 | 国家电网公司 | 基于相似日的bp神经网络光伏发电系统功率预测方法 |
CN106096246A (zh) * | 2016-06-06 | 2016-11-09 | 华北电力大学(保定) | 基于pm2.5和pm10的气溶胶光学厚度估计方法 |
CN106372718A (zh) * | 2016-08-29 | 2017-02-01 | 华北电力大学(保定) | 雾霾积灰条件下的光伏发电功率输出减少率估计方法 |
CN207036621U (zh) * | 2017-07-07 | 2018-02-23 | 宁夏大学 | 一种光伏电池表面沙尘沉积量测量装置 |
Non-Patent Citations (1)
Title |
---|
刘卫亮等: "计及雾霾影响因素的光伏发电超短期功率预测", 《中国电机工程学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109829572A (zh) * | 2019-01-14 | 2019-05-31 | 国网江苏省电力有限公司苏州供电分公司 | 雷电气候下的光伏发电功率预测方法 |
CN109829572B (zh) * | 2019-01-14 | 2021-03-12 | 国网江苏省电力有限公司苏州供电分公司 | 雷电气候下的光伏发电功率预测方法 |
CN111861784A (zh) * | 2020-06-01 | 2020-10-30 | 大唐东北电力试验研究院有限公司 | 基于人工蜂群优化神经网络的光伏发电功率短期预测方法 |
CN111769904A (zh) * | 2020-06-23 | 2020-10-13 | 电子科技大学 | 一种反向散射通信系统中多反射设备并行传输的检测方法 |
CN117740816A (zh) * | 2023-11-10 | 2024-03-22 | 南京旗云中天科技有限公司 | 提高光伏组件积灰污染比检测精度的方法、传感器、系统 |
CN117879490A (zh) * | 2024-03-13 | 2024-04-12 | 国网山东省电力公司微山县供电公司 | 一种光伏设备发电功率预测方法及系统 |
CN117879490B (zh) * | 2024-03-13 | 2024-05-31 | 国网山东省电力公司微山县供电公司 | 一种光伏设备发电功率预测方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109060618A (zh) | 一种光伏发电功率预测装置及方法 | |
CN107341566A (zh) | 基于气象参数与太阳能电池板运行状态参数的光伏系统发电功率预测装置及其方法 | |
CN102156180B (zh) | 一种区域能见度监视预测系统及方法 | |
Hicks | Measuring dry deposition: a re-assessment of the state of the art | |
CN102288222A (zh) | 一种基于无线传感器网络的农情信息实时监测平台 | |
CN213302548U (zh) | 智能一体式固态自动遥测气象站系统 | |
CN104777108B (zh) | 一种叶绿素含量的检测装置及方法 | |
CN102568171A (zh) | 一种用于测量田间农作物土壤墒情的可控无线传输仪 | |
CN106405055A (zh) | 一种连续在线测定土壤co2通量的系统及方法 | |
CN106096246B (zh) | 基于pm2.5和pm10的气溶胶光学厚度估计方法 | |
CN204854800U (zh) | 一种水文水质在线监测及预警系统 | |
CN103217507B (zh) | 一种无线智能碳汇监测系统 | |
CN109656977A (zh) | 一种农业面源污染监测大数据管理系统 | |
CN101587195A (zh) | 一种城市热岛强度的测量方法及其测量装置 | |
CN201666871U (zh) | 雾区能见度检测设备 | |
CN206248138U (zh) | 一种水土流失监测系统 | |
CN103234583A (zh) | 包含测光和测风功能的实时监测系统 | |
CN110907318A (zh) | 一种近地面大气总悬浮颗粒物质量浓度遥感物理估算方法 | |
CN208607220U (zh) | 一种智能远程土壤检测设备 | |
CN103247231A (zh) | 一种新型环保led公交站牌 | |
CN207081442U (zh) | 一种建筑装饰装修工程质量验收设备系统 | |
CN105717048A (zh) | 一种基于养殖水体光学特性的水质采集处理系统 | |
CN205404848U (zh) | 一种基于树莓派的多气象传感器数据采集系统 | |
CN103267717A (zh) | 一种基于太阳能供电模式的环保型pm2.5监测设备 | |
CN203178234U (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181221 |
|
RJ01 | Rejection of invention patent application after publication |