专利类型:发明专利
语 言:中文
申 请 号:CN201310077949.7
申 请 日:20130312
申 请 人:重庆大学
申请人地址:400044 重庆市沙坪坝区沙正街174号
公 开 日:20150916
公 开 号:CN103106314B
代 理 人:胡正顺
代理机构:重庆大学专利中心 50201
摘 要:一种太阳能光伏电源输出功率的时序概率建模方法,属于电力系统电源输出功率的概率建模技术领域。本发明方法利用计算机,通过程序,首先输入光伏电源输出功率及日功率输出起止时刻的实测数据;再基于非参数核密度估计理论估计光伏电源日功率输出起止时刻的概率密度函数,并产生随机样本;然后,基于非参数核密度估计和条件概率理论依次估计光伏电源各时刻输出功率的概率密度函数,并产生随机样本。本发明方法具有方法简单,实用性强,估计的准确度高,通用性强,便于推广应用的特点。本发明可广泛用于电力系统中建立光伏电源输出功率的时序概率模型。
主 权 项:一种太阳能光伏电源输出功率的时序概率建模方法,其特征在于所述方法的具体步骤如下:(1)输入实测数据输入太阳能光伏电源,即光伏电源m天每天t个整点时刻输出功率的实测数据为pki?(k=1,2,…,m,i=1,2,…,t);m天的日功率输出开始时刻ts的实测数据为Ts=[Ts1,?Ts2,..,?Tsm];m天的日功率输出停止时刻te的实测数据为Te=[Te1,?Te2,..,?Tem];(2)产生光伏电源日功率输出开始和停止时刻的随机样本第(1)步完成后,根据第(1)步输入的实测数据,基于非参数核密度估计理论估计光伏电源日功率输出的开始时刻ts和停止时刻te的概率密度函数,并产生随机样本tss与tes,具体步骤如下:1)估计ts和te的概率密度函数根据第(1)步输入的光伏电源日功率输出的开始时刻ts的实测数据Ts和日功率输出的停止时刻te的实测数据Te,采用内点法求解公式(1)所示的优化问题,确定非参数核密度估计的带宽参数hs和he,公式(1)为:(1)式中:n为光伏电源日功率输出开始、停止时刻的实测数据数,exp表示指数函数,Tsi,Tei分别为第i天光伏电源日功率输出开始、停止时刻的实测数据,Tsj,Tej分别为第j天光伏电源日功率输出开始、停止时刻的实测数据,hs和he为非参数核密度估计的带宽参数;然后,基于非参数核密度估计理论估计光伏电源日功率输出的开始时刻ts和停止时刻te的概率密度函数f(ts,?te),计算公式为:
(2)式中:n为光伏电源日功率输出开始、停止时刻的实测数据数,exp表示指数函数,Tsi,Tei分别为第i天光伏电源日功率输出开始、停止时刻的实测数据,hs和he为非参数核密度估计的带宽参数;2)计算ts和te的取值区间第(2)——1)步完成后,根据光伏电源日功率输出的开始时刻ts的实测数据Ts和日功率输出的停止时刻te的实测数据Te,求取光伏电源日功率输出开始时刻ts的取值区间[ats,?bts]以及停止时刻te的取值区间[ate,?bte],计算公式为:
(3)
(4)
(5)
(6)式中:ats,bts分别为光伏电源日功率输出开始时刻ts的取值上、下限,ate,bte分别为光伏电源日功率输出停止时刻te的取值上、下限,n为光伏电源日功率输出开始、停止时刻的实测数据数,Ts1,Ts2,Tsn分别为第1、2、n天光伏电源日功率输出开始时刻的实测数据,Te1,Te2,Ten分别为第1、2、n天光伏电源日功率输出停止时刻的实测数据;3)计算ts和te的概率密度函数最大值第(2)——2)步完成后,将光伏电源日功率输出的开始时刻ts的实测数据Ts和日功率输出的停止时刻te的实测数据Te依次代入公式(2),求取f(ts,?te)在各实测数据处的函数值f(Ts1,?Te1),f(Ts2,?Te2),…,f(Tsn,?Ten),从而计算ts和te的概率密度函数最大值ftmax,计算公式为:
(7)式中:Ts1,Ts2,Tsn分别为第1、2、n天光伏电源日功率输出开始时刻的实测数据,Te1,Te2,Ten分别为第1、2、n天光伏电源日功率输出停止时刻的实测数据,f(Ts1,?Te1)为Ts1和Te1的概率密度函数值,f(Ts2,?Te2)为Ts2和Te2的概率密度函数值,f(Tsn,?Ten)为Tsn和Ten的概率密度函数值;4)产生ts和te的随机样本第(2)——3)步完成后,利用计算机,在区间[0,?1]内产生服从均匀分布的随机数rt及随机向量Rt=[Rts,?Rte],根据公式(8)和(9)计算随机样本ets和ete,并按照四舍五入的原则对ets和ete取整,公式(8)、(9)为:
(8)
(9)然后,根据公式(2)计算ets和ete的概率密度函数值f(ets,?ete),当满足公式(10)所示条件时,将ets作为光伏电源日功率输出开始时刻ts的随机样本tss,将ete作为光伏电源日功率输出停止时刻te的随机样本tes,并令tss=ets,tes=ete;否则,利用计算机,重新在[0,1]区间内产生随机数rt及随机向量Rt,并计算ets、ete和f(ets,?ete),直至公式(10)所示条件满足为止,公式(10)为:
(10)式中:ftmax为ts和te的概率密度函数最大值;(3)产生光伏电源各时刻输出功率的随机样本第(2)步完成后,设光伏电源在一天各时刻的输出功率分别为p1,?p2,…,?pt,随机样本Ps=[ps1,?ps2,…,?pst],初始化循环变量i,并令i=tss,在时段1~tss以及tes~t内光伏电源的输出功率均为0,那么,接下来需要估计时段tss~tes内光伏电源各时刻输出功率的概率密度函数并产生随机样本,具体步骤如下:1)估计时刻i光伏电源输出功率的概率密度函数根据光伏电源第i、i?1个时刻的输出功率pi、pi?1的实测数据,采用内点法求解公式(11)和(12)所示的优化问题,确定非参数核密度估计的带宽参数hi,?hi?1和hp,公式(11)、(12)为:
(11)
(12)式中:hi,hi?1,hp为非参数核密度估计的带宽参数,exp表示指数函数,n表示有n天的光伏电源输出功率实测数据,pxi、pxi分别为光伏电源第i个时刻的输出功率pi在第x、y天的实测数据,px,i?1、py,i?1分别为光伏电源第i?1个时刻的输出功率pi?1在第x、y天的实测数据;然后,基于非参数核密度估计和条件概率理论估计pi的概率密度函数f(pi),计算公式为:
(13)式中:n表示有n天的光伏电源输出功率实测数据,exp表示指数函数,ps,i?1为光伏电源第i?1个时刻的输出功率pi?1的样本值,pj,i?1为光伏电源第i?1个时刻的输出功率pi?1在第j天的实测数据,pji为光伏电源第i个时刻的输出功率pi在第j天的实测数据,hi,?hi?1,?hp为非参数核密度估计的带宽参数;2)计算pi的取值区间第(3)——1)步完成后,根据第i个时刻光伏电源输出功率pi的实测数据,计算其取值区间[api,?bpi],计算公式为:
(14)
(15)式中:api,bpi分别为第i个时刻光伏电源输出功率pi的取值上、下限,p1i,p2i,pni分别为光伏电源第i个时刻的输出功率pi在第1、2、n天的实测数据;3)计算pi的概率密度函数最大值第(3)——2)步完成后,根据pi和pi?1的实测数据,利用公式(13),求取f(pi)在各实测数据处的函数值f(p1i),f(p2i),…,f(pni),并利用公式(16)计算pi的概率密度函数最大值fpimax,公式(16)为:
(16)式中:p1i,p2i,pni分别为光伏电源第i个时刻的输出功率pi在第1、2、n天的实测数据,f(p1i),f(p2i),f(pni)分别为p1i,p2i,pni的概率密度函数值;4)产生pi的随机样本第(3)——3)步完成后,利用计算机,在[0,?1]区间内产生服从均匀分布的随机数rpi、Rpi,根据公式(17)计算随机样本epi,公式(17)为:
(17)然后,根据公式(13)计算epi的概率密度函数值f(epi),当满足公式(18)所示条件时,将epi作为光伏电源时刻i输出功率的随机样本psi,并令psi=epi;否则,利用计算机,重新在[0,1]区间内产生随机数rpi、Rpi,并计算epi和f(epi),直至公式(18)所示条件满足为止,公式(18)为:
(18)式中:fpimax为epi的概率密度函数最大值;最后,当i<tes时,令i=i+1,并返回第(3)——1)步,如此循环,直至i大于等于tes为止;否则,停止计算,输出光伏电源输出功率的随机样本Ps=[ps1,?ps2,…,?pst]。
关 键 词:
法律状态:授权
IPC专利分类号:G06F17/50(2006.01)I