发明名称 基于影响因素筛选的短期负荷预测方法
摘要 本发明提供一种基于影响因素筛选的短期负荷预测方法,用于电力系统短期负荷预测,该方法在对偏最小二乘回归分析的主成分进行有效地提取的基础上,通过对影响负荷因素的重要性的分析,实现对影响负荷因素的筛选,最后将经过筛选的影响因素用于短期负荷预测中;本发明不仅减少了短期负荷预测的工作量,而且保证了负荷预测精度。
申请公布号 CN102073922A 申请公布日期 2011.05.25
申请号 CN201010581653.5 申请日期 2010.12.10
申请人 沈阳工业大学 发明人 卢芸
分类号 G06Q10/00(2006.01)I;G06Q50/00(2006.01)I 主分类号 G06Q10/00(2006.01)I
代理机构 沈阳智龙专利事务所(普通合伙) 21115 代理人 宋铁军
主权项 1.一种基于影响因素筛选的短期负荷预测方法,其特征在于:该方法步骤如下: (1)、建立影响负荷因素筛选模型:如果有<img file="2010105816535100001DEST_PATH_IMAGE002.GIF" wi="17" he="18" />个影响负荷因素<img file="2010105816535100001DEST_PATH_IMAGE004.GIF" wi="84" he="26" />,对它们分别进行了<img file="2010105816535100001DEST_PATH_IMAGE006.GIF" wi="14" he="16" />次采样,得到<img file="822079DEST_PATH_IMAGE006.GIF" wi="14" he="16" />个样本点<img file="2010105816535100001DEST_PATH_IMAGE008.GIF" wi="205" he="26" />则所构成的数据表<img file="2010105816535100001DEST_PATH_IMAGE010.GIF" wi="20" he="18" />可以写成一个<img file="2010105816535100001DEST_PATH_IMAGE012.GIF" wi="38" he="18" />维的矩阵<img file="2010105816535100001DEST_PATH_IMAGE014.GIF" wi="292" he="89" />式中<img file="2010105816535100001DEST_PATH_IMAGE016.GIF" wi="173" he="28" />,<img file="2010105816535100001DEST_PATH_IMAGE018.GIF" wi="80" he="22" />,<img file="2010105816535100001DEST_PATH_IMAGE020.GIF" wi="17" he="25" />被称为第<img file="2010105816535100001DEST_PATH_IMAGE022.GIF" wi="10" he="18" />个样本点;所有的样本点所在的空间被称为样本点空间<img file="2010105816535100001DEST_PATH_IMAGE024.GIF" wi="18" he="18" />;(2)、根据历史数据,对负荷进行主成分的提取:首先将<img file="2010105816535100001DEST_PATH_IMAGE026.GIF" wi="22" he="18" />做标准化处理,<img file="549339DEST_PATH_IMAGE026.GIF" wi="22" he="18" />经标准化处理后的数据矩阵记为<img file="2010105816535100001DEST_PATH_IMAGE028.GIF" wi="145" he="26" />,负荷量<img file="2010105816535100001DEST_PATH_IMAGE030.GIF" wi="17" he="18" />经标准化处理后的数据矩阵记为<img file="2010105816535100001DEST_PATH_IMAGE032.GIF" wi="140" he="26" />;第一步,记<img file="2010105816535100001DEST_PATH_IMAGE034.GIF" wi="14" he="24" />是<img file="2010105816535100001DEST_PATH_IMAGE036.GIF" wi="24" he="25" />的第一个成分,<img file="2010105816535100001DEST_PATH_IMAGE038.GIF" wi="68" he="25" />,<img file="2010105816535100001DEST_PATH_IMAGE040.GIF" wi="21" he="24" />是<img file="851049DEST_PATH_IMAGE036.GIF" wi="24" he="25" />的第一个轴,它是一个单位向量,即<img file="2010105816535100001DEST_PATH_IMAGE042.GIF" wi="53" he="28" />;记<img file="2010105816535100001DEST_PATH_IMAGE044.GIF" wi="18" he="24" />是<img file="2010105816535100001DEST_PATH_IMAGE046.GIF" wi="21" he="25" />的第一个成分,<img file="2010105816535100001DEST_PATH_IMAGE048.GIF" wi="65" he="25" />;<img file="2010105816535100001DEST_PATH_IMAGE050.GIF" wi="17" he="24" />是<img file="682476DEST_PATH_IMAGE046.GIF" wi="21" he="25" />的第一个轴,并且<img file="2010105816535100001DEST_PATH_IMAGE052.GIF" wi="49" he="28" />;如果要求<img file="372215DEST_PATH_IMAGE034.GIF" wi="14" he="24" />,<img file="256994DEST_PATH_IMAGE044.GIF" wi="18" he="24" />能分别很好地代表<img file="2010105816535100001DEST_PATH_IMAGE054.GIF" wi="19" he="17" />与<img file="710978DEST_PATH_IMAGE030.GIF" wi="17" he="18" />中的数据变异信息,根据主成分分析原理,应该有:<img file="2010105816535100001DEST_PATH_IMAGE056.GIF" wi="108" he="49" />另一方面,由于回归建模的需要,又要求<img file="283780DEST_PATH_IMAGE034.GIF" wi="14" he="24" />与<img file="14975DEST_PATH_IMAGE044.GIF" wi="18" he="24" />有最大的解释能力,由典型相关分析的思路,<img file="821388DEST_PATH_IMAGE034.GIF" wi="14" he="24" />与<img file="638035DEST_PATH_IMAGE044.GIF" wi="18" he="24" />的相关度应达到最大值,即:<img file="2010105816535100001DEST_PATH_IMAGE058.GIF" wi="106" he="24" />其中,<img file="2010105816535100001DEST_PATH_IMAGE060.GIF" wi="180" he="49" />为相关系数,<img file="2010105816535100001DEST_PATH_IMAGE062.GIF" wi="73" he="24" />是<img file="2010105816535100001DEST_PATH_IMAGE064.GIF" wi="14" he="24" />与<img file="2010105816535100001DEST_PATH_IMAGE066.GIF" wi="17" he="24" />的协方差,<img file="2010105816535100001DEST_PATH_IMAGE068.GIF" wi="52" he="24" />和<img file="2010105816535100001DEST_PATH_IMAGE070.GIF" wi="56" he="24" />是<img file="185166DEST_PATH_IMAGE064.GIF" wi="14" he="24" />和<img file="770868DEST_PATH_IMAGE066.GIF" wi="17" he="24" />的方差;因此,在偏最小二乘回归中,要求<img file="13762DEST_PATH_IMAGE034.GIF" wi="14" he="24" />与<img file="317704DEST_PATH_IMAGE044.GIF" wi="18" he="24" />的协方差达到最大,即<img file="2010105816535100001DEST_PATH_IMAGE072.GIF" wi="301" he="28" />正规的数学表述是求解下列优化问题,即<img file="2010105816535100001DEST_PATH_IMAGE074.GIF" wi="114" he="84" />因此,将在<img file="2010105816535100001DEST_PATH_IMAGE076.GIF" wi="61" he="30" />和<img file="2010105816535100001DEST_PATH_IMAGE078.GIF" wi="57" he="30" />的约束条件下,去求<img file="2010105816535100001DEST_PATH_IMAGE080.GIF" wi="81" he="25" />的最大值;求得轴<img file="29046DEST_PATH_IMAGE040.GIF" wi="21" he="24" />和<img file="672517DEST_PATH_IMAGE050.GIF" wi="17" he="24" />后,即可得到成分<img file="2010105816535100001DEST_PATH_IMAGE082.GIF" wi="68" he="49" />然后,求<img file="381585DEST_PATH_IMAGE036.GIF" wi="24" he="25" />和<img file="172823DEST_PATH_IMAGE046.GIF" wi="21" he="25" />的三个回归方程<img file="2010105816535100001DEST_PATH_IMAGE084.GIF" wi="100" he="25" /><img file="2010105816535100001DEST_PATH_IMAGE086.GIF" wi="101" he="26" /><img file="2010105816535100001DEST_PATH_IMAGE088.GIF" wi="90" he="25" />式中,回归系数向量是:<img file="2010105816535100001DEST_PATH_IMAGE090.GIF" wi="72" he="52" /><img file="2010105816535100001DEST_PATH_IMAGE092.GIF" wi="70" he="52" /><img file="2010105816535100001DEST_PATH_IMAGE094.GIF" wi="64" he="52" />而<img file="2010105816535100001DEST_PATH_IMAGE096.GIF" wi="22" he="24" />,<img file="2010105816535100001DEST_PATH_IMAGE098.GIF" wi="24" he="25" />,<img file="2010105816535100001DEST_PATH_IMAGE100.GIF" wi="20" he="24" />分别是三个回归方程的残差矩阵;第二步,用残差矩阵<img file="937124DEST_PATH_IMAGE096.GIF" wi="22" he="24" />和<img file="497418DEST_PATH_IMAGE100.GIF" wi="20" he="24" />取代<img file="82114DEST_PATH_IMAGE036.GIF" wi="24" he="25" />和<img file="32753DEST_PATH_IMAGE046.GIF" wi="21" he="25" />,然后,计算得到第二个轴<img file="2010105816535100001DEST_PATH_IMAGE102.GIF" wi="22" he="24" />和<img file="2010105816535100001DEST_PATH_IMAGE104.GIF" wi="18" he="24" />以及第二个成分<img file="2010105816535100001DEST_PATH_IMAGE106.GIF" wi="17" he="24" />,<img file="2010105816535100001DEST_PATH_IMAGE108.GIF" wi="20" he="24" />,有<img file="2010105816535100001DEST_PATH_IMAGE110.GIF" wi="164" he="76" /><img file="600744DEST_PATH_IMAGE102.GIF" wi="22" he="24" />是对应于矩阵<img file="2010105816535100001DEST_PATH_IMAGE112.GIF" wi="69" he="24" />最大特征值<img file="2010105816535100001DEST_PATH_IMAGE114.GIF" wi="21" he="25" />的特征向量,<img file="976116DEST_PATH_IMAGE104.GIF" wi="18" he="24" />是对应于矩阵<img file="2010105816535100001DEST_PATH_IMAGE116.GIF" wi="69" he="24" />最大特征值的特征向量;计算回归系数:<img file="2010105816535100001DEST_PATH_IMAGE118.GIF" wi="73" he="102" />因此,有回归方程:<img file="2010105816535100001DEST_PATH_IMAGE120.GIF" wi="104" he="24" /><img file="2010105816535100001DEST_PATH_IMAGE122.GIF" wi="94" he="24" />如此计算下去,如果<img file="43298DEST_PATH_IMAGE010.GIF" wi="20" he="18" />的秩是<img file="2010105816535100001DEST_PATH_IMAGE124.GIF" wi="18" he="18" />,则会有<img file="2010105816535100001DEST_PATH_IMAGE126.GIF" wi="144" he="25" /><img file="2010105816535100001DEST_PATH_IMAGE128.GIF" wi="166" he="25" />由于<img file="2010105816535100001DEST_PATH_IMAGE130.GIF" wi="57" he="24" />均可以表示成<img file="2010105816535100001DEST_PATH_IMAGE132.GIF" wi="80" he="26" />的线性组合,因此,得到:<img file="2010105816535100001DEST_PATH_IMAGE134.GIF" wi="196" he="28" />,<img file="2010105816535100001DEST_PATH_IMAGE136.GIF" wi="84" he="22" />其中,<img file="2010105816535100001DEST_PATH_IMAGE138.GIF" wi="29" he="25" />是残差矩阵<img file="2010105816535100001DEST_PATH_IMAGE140.GIF" wi="24" he="24" />的第<img file="2010105816535100001DEST_PATH_IMAGE142.GIF" wi="14" he="20" />列;(3)、实现对影响负荷因素的筛选:在偏最小二乘回归分析中,自变量对因变量的解释能力是以变量投影重要性指标VIP来测度的,其定义为:<img file="2010105816535100001DEST_PATH_IMAGE144.GIF" wi="270" he="53" />式中,<img file="2010105816535100001DEST_PATH_IMAGE146.GIF" wi="25" he="26" />是轴<img file="DEST_PATH_IMAGE148.GIF" wi="24" he="25" />的第<img file="DEST_PATH_IMAGE150.GIF" wi="14" he="21" />个分量,用于衡量<img file="DEST_PATH_IMAGE152.GIF" wi="21" he="26" />对构造<img file="DEST_PATH_IMAGE154.GIF" wi="17" he="25" />主成分的贡献大小;<img file="DEST_PATH_IMAGE156.GIF" wi="66" he="25" />和<img file="DEST_PATH_IMAGE158.GIF" wi="106" he="25" />分别代表了<img file="789887DEST_PATH_IMAGE154.GIF" wi="17" he="25" />对<img file="522089DEST_PATH_IMAGE030.GIF" wi="17" he="18" />的解释能力和<img file="DEST_PATH_IMAGE160.GIF" wi="57" he="25" />对<img file="870025DEST_PATH_IMAGE030.GIF" wi="17" he="18" />的累计解释能力,且有:<img file="DEST_PATH_IMAGE162.GIF" wi="294" he="46" /><img file="DEST_PATH_IMAGE164.GIF" wi="202" he="46" />由于影响负荷因素<img file="419692DEST_PATH_IMAGE152.GIF" wi="21" he="26" />对负荷量<img file="157972DEST_PATH_IMAGE030.GIF" wi="17" he="18" />的解释是通过<img file="444597DEST_PATH_IMAGE154.GIF" wi="17" he="25" />来传递的,如果<img file="817679DEST_PATH_IMAGE154.GIF" wi="17" he="25" />对负荷量<img file="164346DEST_PATH_IMAGE030.GIF" wi="17" he="18" />的解释能力很强,而影响负荷因素<img file="389922DEST_PATH_IMAGE152.GIF" wi="21" he="26" />在构造<img file="417921DEST_PATH_IMAGE154.GIF" wi="17" he="25" />时又起到了相当重要的作用,则影响负荷因素<img file="661821DEST_PATH_IMAGE152.GIF" wi="21" he="26" />对负荷量<img file="163078DEST_PATH_IMAGE030.GIF" wi="17" he="18" />的解释能力就被视为很大,进而实现对影响负荷因素的筛选;(4)、定义出影响负荷因素的测定指标:影响系数:在此作出如下定义:影响负荷因素的测定指标用影响系数来表示:<img file="DEST_PATH_IMAGE166.GIF" wi="302" he="53" />(5)、去除影响系数低的影响负荷因素,将影响系数高的影响负荷因素用于神经网络的输入中,实现对影响负荷因素筛选,采用神经网络模型对短期负荷进行预测。
地址 110870 辽宁省沈阳市经济技术开发区沈辽西路111号