CN104281986B - 一种微电网功率预测方法 - Google Patents
一种微电网功率预测方法 Download PDFInfo
- Publication number
- CN104281986B CN104281986B CN201410608857.1A CN201410608857A CN104281986B CN 104281986 B CN104281986 B CN 104281986B CN 201410608857 A CN201410608857 A CN 201410608857A CN 104281986 B CN104281986 B CN 104281986B
- Authority
- CN
- China
- Prior art keywords
- power
- module
- micro
- capacitance sensor
- output
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000005611 electricity Effects 0.000 claims abstract description 34
- 238000004146 energy storage Methods 0.000 claims abstract description 24
- 230000008859 change Effects 0.000 claims abstract description 8
- 210000002569 neuron Anatomy 0.000 claims description 52
- 238000012544 monitoring process Methods 0.000 claims description 45
- 230000007935 neutral effect Effects 0.000 claims description 15
- 230000003203 everyday effect Effects 0.000 claims description 9
- 238000005286 illumination Methods 0.000 claims description 9
- 230000004044 response Effects 0.000 claims description 9
- 238000007726 management method Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000012549 training Methods 0.000 claims description 6
- 238000013500 data storage Methods 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000010248 power generation Methods 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 230000001276 controlling effect Effects 0.000 claims description 3
- 230000006866 deterioration Effects 0.000 claims description 3
- 238000007599 discharging Methods 0.000 claims description 3
- 230000001537 neural effect Effects 0.000 claims description 3
- 238000003062 neural network model Methods 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 3
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 2
- 240000002853 Nelumbo nucifera Species 0.000 claims description 2
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 238000011897 real-time detection Methods 0.000 claims description 2
- 230000001105 regulatory effect Effects 0.000 claims description 2
- 210000005036 nerve Anatomy 0.000 claims 1
- 238000001514 detection method Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000001172 regenerating effect Effects 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 210000000352 storage cell Anatomy 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Entrepreneurship & Innovation (AREA)
- Mathematical Physics (AREA)
- Computing Systems (AREA)
- Molecular Biology (AREA)
- Evolutionary Computation (AREA)
- Development Economics (AREA)
- Data Mining & Analysis (AREA)
- Game Theory and Decision Science (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- General Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明涉及一种微电网功率预测方法,该预测方法可以预测微电网的光伏模块、风电模块的发电功率,预测负载的变化情况,并基于实时检测的蓄电池模块的储能情况和实时获取的大电网的运行情况,来预测微电网的功率供给策略,使得微电网处于安全、经济、令用户满意的运行状态之下。
Description
所属技术领域
本发明涉一种微电网功率预测方法。
背景技术
随着全球性能源危机的不断加大,低碳、清洁的可再生能源已经成为研究热点,目前世界各国均开展了含可再生能源的微电网研究项目。微电网是指多个分布式电源及其相关负载按照一定的拓扑结构组成的网络,并通过静态开关关联至常规电网。
大量的风/光分布式发电以微电网的形式作为一个双向可调度单元从低压配电网接入,形成多个微电网组成的微网群或主动型配电网。然而,由于风力发电、光伏发电等可再生电源的高度随机性,传统负荷、燃料电池等清洁电源的不确定性,以及弹性负荷的可控可调,使得微电网成为公共电网中的不确定负载,对公共电网的稳定与优化带来巨大的冲击。
因此,要保证微电网安全、可靠、经济地运行,实现微电网的推广和应用,就需要对微电网供电管理问题进行研究。与负荷被动参与电网运行的传统电网供电管理相比,微电网的供电管理更期望用户负荷主动参与电网的运行,实现对对电能的积极消费。但随着智能用户数量的增加,用户时变的用电数量和随机的用电时间将给负荷端的能量分配带来扰动。另一方面,微电网的供电端引入了大量可再生、节能环保的新能源,新能源的随机性也给微电网的供电管理带来了扰动的影响。
解决上述问题的一个关键因素在于微电网的功率预测。
发明内容
为解决上述问题,本发明提供一种微电网功率预测方法,该预测方法可以预测微电网的光伏模块、风电模块的发电功率,预测负载的变化情况,并基于实时检测的蓄电池模块的储能情况和实时获取的大电网的运行情况,来预测微电网的功率供给策略,使得微电网处于安全、经济、令用户满意的运行状态之下。
为了实现上述目的,本发明提供一种微电网功率预测方法,该方法基于如下监控***来实现,该监控***包括:
光伏发电监控模块,用于实时监控微电网中的光伏发电模块,并对光伏发电模块的发电功率进行预测;
风电监控模块,用于实时监控风电网中的风电模块,并对风电模块的发电功率进行预测;
蓄电池监控模块,用于实时监控微电网中的蓄电池模块;
负载监控模块,用于实时监控微电网中的负载,并对负载的功率变化情况进行预测;
大电网联络模块,用于实时从大电网调控中心获知大电网的运行情况以及相关调度信息;
中控模块,用于确定微电网的供电策略,并向上述各模块发出指令,以执行该供电策略;
总线模块,用于该监控***的各个模块的通信联络;
该预测方法包括如下步骤:
(1)光伏发电监控模块、风电监控模块实时获取微电网中光伏发电模块和风电模块的运行数据,并存储数据,负载监控模块实时获取负载的负荷变化情况;
(2)根据现有微电网中光伏发电模块和风电模块的运行数据,对未来预定时刻内的微电网中的输出功率进行预测,根据现有微电网中的负载的负荷变化情况,对负载的负荷需求进行预测;
(3)实时检测获取蓄电池模块的储能充放容量,以及大电网的调度信息;
(4)以大电网的调度信息、当前蓄电池储能的充放容量、未来微电网中的输出功率以及对未来负荷需求的变化作为约束条件,建立微电网供电管理的目标函数,来预测微电网在未来一段时间内的功率供给情况。
优选的,在步骤(2)中,基于实时气象数据对光伏发电、风机发电的输出功率进行预测,具体步骤包括:
(201)对温度、光照度、风速进行预测;
(202)然后利用光伏模块和风机出力的数学模型对其输出功率进行预测。
优选的,光伏发电模块的输出功率可以采用如下的经验公式,该公式采用与标准环境下的光照温度进行比对的方法,给出近似参考值,该公式表述如下:
PPV=Pstc×Gw×[1+kr(Tw-Tstc)]/Gstc
式中:Pstc、Tstc、Gstc分别表示标准环境下的最大测试功率、环境参考温度以及光照强度,其中Tstc=300K,Gstc=1000W/m2;Gw、Tw分别表示光伏发电模块工作时的光照强度和温度;kr表示光伏发电模块功率温度系数。
优选的,在步骤(201)中,采用如下方式进行温度预测:
记录样本中的温度数据T1[24]、T2[24]…Tm[24],ΔT1[23]、ΔT2[23]、ΔT3[23]、ΔT4[23]、ΔT5[23];
运行到t时刻监测的温度T[t],ΔT[t-1];
计算温度变化趋势的相似度;
对相似度进行归一化处理;
对t时刻后温度进行预测,得到T[24-t]。
优选的,在步骤(4)中,微电网总功率Pg的约束为:
在非响应调度时段1下,Pg,min≤Pg(l)≤Pg,max,Pg,min为微电网能够从大电网吸收的最大功率,Pg,max为微电网能够向大电网输送功率的最大功率;
在响应调度时段2下,Pg(2)=Pset,Pset为响应调度时段2下要求的联络线功率。
优选的,在步骤(2)中,采用神经网络模型预测负载负荷需求,具体步骤如下:
(211)每一天采集12组负载的有功功率和无功功率,共连续采集8天,这样共有96组数据P(k)和Q(k),k=1,2,…,96;
(212)将96组数据P(k)和Q(k)进行归一化处理,使得 n=1,2,…,96;首先将每一天的12个有功功率P(k)作为一组输入矢量R(m),12个无功功率Q(k)作为一组输入矢量S(m),m=1,2,…,8,m表示神经网络的训练次数;同时预先假设第9天的12个有功功率P′(k)作为预测功率的输出矢量R′,第9天的12个无功功率Q′(k)作为预测功率的输出矢量S′;这样前8天的有功功率输入矢量就为R(1),R(2),R(3),R(4),R(5),R(6),R(7),R(8),第9天预测有功功率的输出矢量为R′;前8天的无功功率输入矢量就为S(1),S(2),S(3),S(4),S(5),S(6),S(7),S(8),第9天预测有功功率的输出矢量为S′;
(213)将8组输入矢量R(m)和S(m)作为神经网络的输入层,隐含层神经元的传递函数采用S型正切函数tansig,输出层神经元的传递函数采用S型对数函数logsig,如图2所示,这样经过8次神经网络训练后,就确定了神经网络中各连接权的权值;;
(214)对于8个有功功率输入矢量R(m),在隐含层神经元有a1=tan sig(IW1R+b1),其中a1为隐含层神经元输出,IW1为隐含层神经元的权值,b1为隐含层神经元的阈值;在输出层神经元有a2=log sig(LW2a1+b2),其中a2为输出层神经元输出,IW2为输出层神经元的权值,b2为输出层神经元的阈值;
(215)对于8个有功功率输入矢量S(m),在隐含层神经元有c1=tan sig(IW1S+b1),其中c1为隐含层神经元输出,IW1为隐含层神经元的权值,b1为隐含层神经元的阈值;在输出层神经元有c2=log sig(LW2c1+b2),其中c2为输出层神经元输出,IW2为输出层神经元的权值,b2为输出层神经元的阈值;
(216)将第8天的输入矢量R(8)和S(8)再次作为神经网络的输入层,此时神经网络中输出的预测功率的输出矢量R′和S′即为第九天的功率预测归一化值,再用反归一化算法,即k=1,2,…,96,输出的矢量值和就是第九天预测功率的12个有功功率P′(k)和12个无功功率Q′(k)。这样以此类推,可以重复上面的步骤利用第二天到第九天的数据预测到第十天的功率,这样后面每一天的功率都可以被预测出来。
优选的,在步骤(4)中,微电网蓄电池容量约束为:
[Socref-Soc(k+1)]=a[Socref-Soc(k)]+ηEs(k)
Socmin≤Soc(k)≤Socmax
其中,Soc(k)为k时刻蓄电池的容量状态,Socref是为保证储能电池工作的可靠性而设置的一个设定值,Es(k)代表储能设备与其它电力设备间流动的电量,储能的物理损耗系数a∈(0,1),η为蓄电池的充放电效率,充电效率记为ηc、放电效率记为ηd,且它们之间满足以下关系:
储能的充放电过程可以看作是一个同时包含连续变量和离散变量的动态过程,这里采用混合逻辑动态模型处理方法,通过引入二进制变量δ(k)来表示储能在当前时刻的工作状态,
Z(k)=δ(k)Es(k)
Z(k)表示当前时刻储能充/放电的电量,则蓄电池的动态特性可描述为:
[Socref-Soc(k+1)]=a[Socref-Soc(k)]+(ηc-ηd)Z(k)+ηdEs(k)
满足以下约束条件:E1δ(k)+E2Z(k)≤E3Es(k)+E4
其中,系数矩阵E1,E2,E3和E4是在将逻辑命题转换为线性不等式时二进制变量和连续变量要满足的线性不等式约束,可通过数学公式的推导得到。
本发明的监控方法具有如下优点:(1)准确预测微电网的功率变化情况;(2)供电策略兼顾大电网调度要求、微电网运行情况和负载的负荷需求,满足用户同时,兼顾了供电可靠性,同时提高了供电效益。
附图说明
图1示出了本发明方法所使用的一种微电网监控***的框图;
图2示出了本发明方法的流程图。
具体实施方式
图1是示出了本发明的一种微电网监控***100,该***100包括:光伏发电监控模块104,用于实时监控微电网200中的光伏发电模块201,并对光伏发电模块201的发电功率进行预测;风电监控模块105,用于实时监控风电网200中的风电模块202,并对风电模块202的发电功率进行预测;蓄电池监控模块106,用于实时监控微电网200中的蓄电池模块203;负载监控模块108,用于实时监控微电网200中的负载204,并对负载204的功率变化情况进行预测;大电网联络模块102,用于实时从大电网300调控中心获知大电网300的运行情况以及相关调度信息;断路器109,用于连接或隔离微电网200与大电网300;并网监控模块103,用于控制断路器109连接或隔离微电网200和大电网300;中控模块107,用于确定微电网200的供电策略,并向上述各模块发出指令,以执行该供电策略;总线模块101,用于该监控***的各个模块的通信联络。
通信模块101,用于上述各个模块之间的通信,所述总线通信模块101通过冗余双CAN总线与其他模块相连。
光伏发电模块201至少包括电压、电流检测设备以及太阳光强度检测设备和温度检测设备。风电监控模块105至少包括风力发电机定压、电流、频率检测设备,以及风速检测设备。风力发电机的输出功率由风力发电机所在地点的风速、风向和自身特征所决定。
基于实时气象数据的光伏发电、风机发电预测,需对温度、光照度、风速进行预测。其中,温度预测可采用以下方法:
样本中的温度数据T1[24]、T2[24]…Tm[24],ΔT1[23]、ΔT2[23]、ΔT3[23]、ΔT4[23]、ΔT5[23];运行到t时刻监测的温度T[t],ΔT[t-1];计算温度变化趋势的相似度;对相似度进行归一化处理;对t时刻后温度进行预测,得到T[24-t]:
此外,对光照度、风速预测时,可按与温度预测类似的方法进行预测,然后利用光伏模块和风机出力的数学模型对其输出功率进行预测。
蓄电池监控模块106至少包括蓄电池端电压、电流检测设备以及温度检测设备。用于实时监控蓄电池模块的充放电过程。通过调节蓄电池的充/放电来存储/补充富余/缺少的能量,蓄电池的能量流动可描述为,
[Socref-Soc(k+1)]=a[Socref-Soc(k)]+ηEs(k)
Socmin≤Soc(k)≤Socmax
其中,Soc(k)为k时刻蓄电池的容量状态,Socref是为保证储能电池工作的可靠性而设置的一个设定值,Es(k)代表储能设备与其它电力设备间流动的电量,储能的物理损耗系数a∈(0,1),η为蓄电池的充放电效率,充电效率记为ηc、放电效率记为ηd,且它们之间满足以下关系:
储能的充放电过程可以看作是一个同时包含连续变量和离散变量的动态过程,这里采用混合逻辑动态模型处理方法,通过引入二进制变量δ(k)来表示储能在当前时刻的工作状态,
Z(k)=δ(k)Es(k)
Z(k)表示当前时刻储能充/放电的电量,则蓄电池的动态特性可描述为:
[Socref-Soc(k+1)]=a[Socref-Soc(k)]+(ηc-ηd)Z(k)+ηdEs(k)
满足以下约束条件:E1δ(k)+E2Z(k)≤E3Es(k)+E4
其中,系数矩阵E1,E2,E3和E4是在将逻辑命题转换为线性不等式时二进制变量和连续变量要满足的线性不等式约束,可通过数学公式的推导得到。
在将逻辑命题转换为线性不等式时二进制变量和连续变量过程中,要满足的线性不等式约束E1δ(k)+E2Z(k)≤E3Es(k)+E4,其中系数矩阵E1,E2,E3和E4分别为:
中控模块107至少包括CPU单元、数据存储单元和显示单元。
大电网联络模块102至少包括无线通信设备。该无线通信设备可以为有线设备或无线设备。
参见附图2,本发明的方法包括如下步骤:
S1.光伏发电监控模块104、风电监控模块实105时获取微电网200中光伏发电模块201和风电模块202的运行数据,并存储数据,负载监控模块108实时获取负载204的负荷变化情况;
S2.根据现有微电网200中光伏发电模块201和风电模块202的运行数据,对未来预定时刻内的微电网200中的输出功率进行预测,根据现有微电网200中的负载204的负荷变化情况,对负载204的负荷需求进行预测;
S3.蓄电池监控模块106实时检测获取蓄电池模块203的储能充放容量,大电网联络模块102实时检测大电网的调度信息;
S4.以大电网300的调度信息、当前蓄电池模块203储能的充放容量、未来微电网中的输出功率以及对未来负荷需求的变化作为约束条件,建立微电网供电管理的目标函数,来预测微电网在未来一段时间内的功率供给情况。
在步骤S2中,基于实时气象数据对光伏发电、风机发电的输出功率进行预测,具体步骤包括:
S201.对温度、光照度、风速进行预测;
S202.然后利用光伏模块和风机出力的数学模型对其输出功率进行预测。其中,温度预测可采用以下方法:
在步骤S201中,采用如下方式进行温度预测:
记录样本中的温度数据T1[24]、T2[24]…Tm[24],ΔT1[23]、ΔT2[23]、ΔT3[23]、ΔT4[23]、ΔT5[23];
运行到t时刻监测的温度T[t],ΔT[t-1];
计算温度变化趋势的相似度;
对相似度进行归一化处理;
对t时刻后温度进行预测,得到T[24-t]。
光伏发电模块的输出功率可以采用如下的经验公式,该公式采用与标准环境下的光照温度进行比对的方法,给出近似参考值。该公式表述如下:
PPV=Pstc×Gw×[1+kr(Tw-Tstc)]/Gstc
式中:Pstc、Tstc、Gstc分别表示标准环境下的最大测试功率、环境参考温度以及光照强度,其中Tstc=300K,Gstc=1000W/m2;Gw、Tw分别表示光伏发电模块工作时的光照强度和温度;kr表示光伏发电模块功率温度系数。
在S2中,采用神经网络模型预测负荷需求,具体步骤如下:
S211.每一天采集12组有功功率和无功功率,共连续采集8天,这样共有96组数据P(k)和Q(k),k=1,2,…,96。
S212.将96组数据P(k)和Q(k)进行归一化处理,使得 n=1,2,…,96;首先将每一天的12个有功功率P(k)作为一组输入矢量R(m),12个无功功率Q(k)作为一组输入矢量S(m),m=1,2,…,8,m表示神经网络的训练次数;同时预先假设第9天的12个有功功率P′(k)作为预测功率的输出矢量R′,第9天的12个无功功率Q′(k)作为预测功率的输出矢量S′;这样前8天的有功功率输入矢量就为R(1),R(2),R(3),R(4),R(5),R(6),R(7),R(8),第9天预测有功功率的输出矢量为R′;前8天的无功功率输入矢量就为S(1),S(2),S(3),S(4),S(5),S(6),S(7),S(8),第9天预测有功功率的输出矢量为S′。
S213.将8组输入矢量R(m)和S(m)作为神经网络的输入层,隐含层神经元的传递函数采用S型正切函数tansig,输出层神经元的传递函数采用S型对数函数logsig,如图2所示,这样经过8次神经网络训练后,就确定了神经网络中各连接权的权值。
S214.对于8个有功功率输入矢量R(m),在隐含层神经元有a1=tan sig(IW1R+b1),其中a1为隐含层神经元输出,IW1为隐含层神经元的权值,b1为隐含层神经元的阈值;在输出层神经元有a2=log sig(LW2a1+b2),其中a2为输出层神经元输出,IW2为输出层神经元的权值,b2为输出层神经元的阈值。
S215.对于8个有功功率输入矢量S(m),在隐含层神经元有c1=tan sig(IW1S+b1),其中c1为隐含层神经元输出,IW1为隐含层神经元的权值,b1为隐含层神经元的阈值;在输出层神经元有c2=log sig(LW2c1+b2),其中c2为输出层神经元输出,IW2为输出层神经元的权值,b2为输出层神经元的阈值。
S216.将第8天的输入矢量R(8)和S(8)再次作为神经网络的输入层,此时神经网络中输出的预测功率的输出矢量R′和S′即为第九天的功率预测归一化值,再用反归一化算法,即k=1,2,…,96,输出的矢量值和就是第九天预测功率的12个有功功率P′(k)和12个无功功率Q′(k)。这样以此类推,可以重复上面的步骤利用第二天到第九天的数据预测到第十天的功率,这样后面每一天的功率都可以被预测出来。
在步骤S4中,微电网总功率Pg的约束为:
在非响应调度时段1下,Pg,min≤Pg(l)≤Pg,max,Pg,min为微电网200能够从大电网300吸收的最大功率,Pg,max为微电网200能够向大电网300输送功率的最大功率;
在响应调度时段2下,Pg(2)=Pset,Pset为响应调度时段2下要求的联络线功率。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,做出若干等同替代或明显变型,而且性能或用途相同,都应当视为属于本发明的保护范围。
Claims (5)
1.一种微电网功率预测方法,该方法基于如下监控***来实现,该监控***包括:
光伏发电监控模块,用于实时监控微电网中的光伏发电模块,并对光伏发电模块的发电功率进行预测;
风电监控模块,用于实时监控风电网中的风电模块,并对风电模块的发电功率进行预测;
蓄电池监控模块,用于实时监控微电网中的蓄电池模块;
负载监控模块,用于实时监控微电网中的负载,并对负载的功率变化情况进行预测;
大电网联络模块,用于实时从大电网调控中心获知大电网的运行情况以及相关调度信息;
中控模块,用于确定微电网的供电策略,并向上述各模块发出指令,以执行该供电策略;
总线模块,用于该监控***的各个模块的通信联络;
其特征在于,该预测方法包括如下步骤:
(1)光伏发电监控模块、风电监控模块实时获取微电网中光伏发电模块和风电模块的运行数据,并存储数据,负载监控模块实时获取负载的负荷变化情况;
(2)根据现有微电网中光伏发电模块和风电模块的运行数据,对未来预定时刻内的微电网中的输出功率进行预测,根据现有微电网中的负载的负荷变化情况,对负载的负荷需求进行预测;
(3)实时检测获取蓄电池模块的储能充放容量,以及大电网的调度信息;
(4)以大电网的调度信息、当前蓄电池储能的充放容量、未来微电网中的输出功率以及对未来负荷需求的变化作为约束条件,建立微电网供电管理的目标函数,来预测微电网在未来一段时间内的功率供给情况;
在步骤(2)中,基于实时气象数据对光伏发电、风机发电的输出功率进行预测,具体步骤包括:
(201)对温度、光照度、风速进行预测;
(202)然后利用光伏模块和风机出力的数学模型对其输出功率进行预测;
光伏发电模块的输出功率可以采用如下的经验公式,该公式采用与标准环境下的光照温度进行比对的方法,给出近似参考值,该公式表述如下:
PPV=Pstc×Gw×[1+kr(Tw-Tstc)]/Gstc
式中:Pstc、Tstc、Gstc分别表示标准环境下的最大测试功率、环境参考温度以及光照强度,其中Tstc=300K,Gstc=1000W/m2;Gw、Tw分别表示光伏发电模块工作时的光照强度和温度;kr表示光伏发电模块功率温度系数。
2.如权利要求1所述的方法,其特征在于,在步骤(201)中,采用如下方式进行温度预测:
记录样本中的温度数据T1[24]、T2[24]…Tm[24],ΔT1[23]、ΔT2[23]、ΔT3[23]、ΔT4[23]、ΔT5[23];
运行到t时刻监测的温度T[t],ΔT[t-1];
计算温度变化趋势的相似度;
对相似度进行归一化处理;
对t时刻后温度进行预测,得到T[24-t]。
3.如权利要求2所述的方法,其特征在于,在步骤(2)中,采用神经网络模型预测负载负荷需求,具体步骤如下:
(211)每一天采集12组负载的有功功率和无功功率,共连续采集8天,这样共有96组数据P(k)和Q(k),k=1,2,…,96;
(212)将96组数据P(k)和Q(k)进行归一化处理,使得 n=1,2,…,96;首先将每一天的12个有功功率P(k)作为一组输入矢量R(m),12个无功功率Q(k)作为一组输入矢量S(m),m=1,2,…,8,m表示神经网络的训练次数;同时预先假设第9天的12个有功功率P′(k)作为预测功率的输出矢量R′,第9天的12个无功功率Q′(k)作为预测功率的输出矢量S′;这样前8天的有功功率输入矢量就为R(1),R(2),R(3),R(4),R(5),R(6),R(7),R(8),第9天预测有功功率的输出矢量为R′;前8天的无功功率输入矢量就为S(1),S(2),S(3),S(4),S(5),S(6),S(7),S(8),第9天预测有功功率的输出矢量为S′;
(213)将8组输入矢量R(m)和S(m)作为神经网络的输入层,隐含层神经元的传递函数采用S型正切函数tansig,输出层神经元的传递函数采用S型对数函数logsig,如图2所示,这样经过8次神经网络训练后,就确定了神经网络中各连接权的权值;
(214)对于8个有功功率输入矢量R(m),在隐含层神经元有a1=tansig(IW1R+b1),其中a1为隐含层神经元输出,IW1为隐含层神经元的权值,b1 为隐含层神经元的阈值;在输出层神经元有a2=log sig(LW2a1+b2),其中a2为输出层神经元输出,IW2为输出层神经元的权值,b2为输出层神经元的阈值;
(215)对于8个有功功率输入矢量S(m),在隐含层神经元有c1=tansig(IW1S+b1),其中c1为隐含层神经元输出,IW1为隐含层神经元的权值,b1为隐含层神经元的阈值;在输出层神经元有c2=log sig(LW2c1+b2),其中c2为输出层神经元输出,IW2为输出层神经元的权值,b2为输出层神经元的阈值;
(216)将第8天的输入矢量R(8)和S(8)再次作为神经网络的输入层,此时神经网络中输出的预测功率的输出矢量R′和S′即为第九天的功率预测归一化值,再用反归一化算法,即k=1,2,…,96,输出的矢量值和就是第九天预测功率的12个有功功率P′(k)和12个无功功率Q′(k),这样以此类推,可以重复上面的步骤利用第二天到第九天的数据预测到第十天的功率,这样后面每一天的功率都可以被预测出来。
4.如权利要求3所述的方法,其特征在于,在步骤(4)中,微电网总功率Pg的约束为:
在非响应调度时段1下,Pg,min≤Pg(l)≤Pg,max,Pg,min为微电网能够从大电网吸收的最大功率,Pg,max为微电网能够向大电网输送功率的最大功率;
在响应调度时段2下,Pg(2)=Pset,Pset为响应调度时段2下要求的联络线功率。
5.如权利要求4所述的方法,其特征在于,在步骤(4)中,微电网蓄电池容量约束为:
[Socref-Soc(k+1)]=a[Socref-Soc(k)]+ηEs(k)
Socmin≤Soc(k)≤Socmax
其中,Soc(k)为k时刻蓄电池的容量状态,Socref是为保证储能电池工作的可靠性而设置的一个设定值,Es(k)代表储能设备与其它电力设备间流动的电量,储能的物理损耗系数a∈(0,1),η为蓄电池的充放电效率,充电效率记为ηc、放电效率记为ηd,且它们之间满足以下关系:
储能的充放电过程可以看作是一个同时包含连续变量和离散变量的动态过程,这里采用混合逻辑动态模型处理方法,通过引入二进制变量δ(k)来表示储能在当前时刻的工作状态,
Z(k)=δ(k)Es(k)
Z(k)表示当前时刻储能充/放电的电量,则蓄电池的动态特性可描述为:
[Socref-Soc(k+1)]=a[Socref-Soc(k)]+(ηc-ηd)Z(k)+ηdEs(k)
满足以下约束条件:E1 δ(k)+E2Z(k)≤E3Es(k)+E4
其中,系数矩阵E1,E2,E3和E4是在将逻辑命题转换为线性不等式时二进制变量和连续变量要满足的线性不等式约束,可通过数学公式的推导得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410608857.1A CN104281986B (zh) | 2014-11-03 | 2014-11-03 | 一种微电网功率预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410608857.1A CN104281986B (zh) | 2014-11-03 | 2014-11-03 | 一种微电网功率预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104281986A CN104281986A (zh) | 2015-01-14 |
CN104281986B true CN104281986B (zh) | 2017-09-26 |
Family
ID=52256837
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410608857.1A Active CN104281986B (zh) | 2014-11-03 | 2014-11-03 | 一种微电网功率预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104281986B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105375533B (zh) * | 2015-09-23 | 2018-02-16 | 许继集团有限公司 | 一种含风、光可再生能源的独立微电网运行模式预测方法 |
CN105186570B (zh) * | 2015-10-19 | 2019-04-12 | 国网北京市电力公司 | 微电网供电控制方法和装置 |
CN105826944B (zh) * | 2016-03-18 | 2018-08-14 | 上海电机学院 | 一种微电网群的功率预测方法和*** |
CN107886445B (zh) * | 2017-11-09 | 2020-12-04 | 华北电力大学 | 一种基于神经元大数据分析的功率调整方法 |
CN109961182A (zh) * | 2019-03-18 | 2019-07-02 | 珠海格力电器股份有限公司 | 供电数据管理方法、装置、计算机设备和存储介质 |
CN114172256B (zh) * | 2021-12-20 | 2022-09-27 | 武汉美格科技股份有限公司 | 一种太阳能发电智能控制装置 |
CN114844114B (zh) * | 2022-06-20 | 2022-11-04 | 北京信云筑科技有限责任公司 | 一种基于分布式光储的配电网远程安全运维方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101414751A (zh) * | 2008-11-20 | 2009-04-22 | 北京方鸿溪科技有限公司 | 风电功率预测***及其方法、电网*** |
CN101630840A (zh) * | 2009-08-12 | 2010-01-20 | 电子科技大学 | 微电网能量智能控制*** |
CN101706335A (zh) * | 2009-11-11 | 2010-05-12 | 华南理工大学 | 一种基于遗传算法优化bp神经网络的风电功率预测方法 |
CN102355057A (zh) * | 2011-09-25 | 2012-02-15 | 国网电力科学研究院 | 微电网***计算机监控方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7715950B2 (en) * | 2007-03-01 | 2010-05-11 | Wisconsin Alumni Research Foundation | Non-inverter based distributed energy resource for use in a dynamic distribution system |
-
2014
- 2014-11-03 CN CN201410608857.1A patent/CN104281986B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101414751A (zh) * | 2008-11-20 | 2009-04-22 | 北京方鸿溪科技有限公司 | 风电功率预测***及其方法、电网*** |
CN101630840A (zh) * | 2009-08-12 | 2010-01-20 | 电子科技大学 | 微电网能量智能控制*** |
CN101706335A (zh) * | 2009-11-11 | 2010-05-12 | 华南理工大学 | 一种基于遗传算法优化bp神经网络的风电功率预测方法 |
CN102355057A (zh) * | 2011-09-25 | 2012-02-15 | 国网电力科学研究院 | 微电网***计算机监控方法 |
Non-Patent Citations (1)
Title |
---|
《风光互补型微电网规划设计典型方案 》;栗晓政 等;《电气时代》;20111010(第10期);94-96 * |
Also Published As
Publication number | Publication date |
---|---|
CN104281986A (zh) | 2015-01-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104281986B (zh) | 一种微电网功率预测方法 | |
CN104319768B (zh) | 一种微电网供电及监控方法 | |
Brekken et al. | Optimal energy storage sizing and control for wind power applications | |
Khatod et al. | Evolutionary programming based optimal placement of renewable distributed generators | |
CN104348188A (zh) | 一种分布式电源运行及监控方法 | |
CN104348189B (zh) | 一种分布式电源*** | |
CN104485681B (zh) | 一种风电场储能***的监控方法 | |
CN104268806A (zh) | 一种微电网功率监控*** | |
CN104318494A (zh) | 一种分布式电源智能监控*** | |
CN105787270B (zh) | 一种输电线路多时间尺度负载能力动态预测方法 | |
CN106875035B (zh) | 一种主动配电网可靠性评估优化方法 | |
Bogaraj et al. | A novel energy management scheme using ANFIS for independent microgrid | |
CN106712075A (zh) | 一种考虑风电并网***安全约束的调峰策略优化方法 | |
Valdés et al. | Operational simulation of wind power plants for electrolytic hydrogen production connected to a distributed electricity generation grid | |
CN114430165A (zh) | 基于深度模型预测的微网群智能协调控制方法和装置 | |
CN113131529A (zh) | 一种计及多种灵活性资源的可再生能源承载力评估的方法 | |
Akbari-Zadeh et al. | Dstatcom allocation in the distribution system considering load uncertainty | |
CN106786789B (zh) | 微网发电量实时控制***及方法 | |
CN109638837B (zh) | 一种面向新能源并网的动态负载管理方法及*** | |
CN114759616B (zh) | 一种考虑电力电子器件特性的微电网鲁棒优化调度方法 | |
Li et al. | Optimal operation of AC/DC hybrid microgrid under spot price mechanism | |
Thorat et al. | Optimization of fuel cost incorporating with wind, solar PV and Electric vehicle energy sources using improved artificial bee colony algorithm | |
CN104505850A (zh) | 一种风电场储能*** | |
CN104485682A (zh) | 一种风电场储能***的监控装置 | |
CN111585305B (zh) | 适用于多能互补联动经济性评价的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |