发明名称 智能电网建设中基于相关向量机的旋转热备用的调度方法
摘要 智能电网建设中基于相关向量机的旋转热备用的调度方法,本发明涉及一种基于相关向量机的旋转热备用的调度方法。为解决含规模化风电等新能源电力系统难以设置旋转热备用以平抑风电并网功率波动问题。将初始化设置结果传递到风功率相关向量机预测系统中;风电场风功率采集模块实时将风电场风功率的测量值进行采集,进行数据预处理后将数据传递给风功率相关向量机预测系统中;风功率相关向量机预测系统接收数据,对未来时刻的风功率进行预测,预测的结果为未来时刻的风功率值和风功率的误差带;将得到的预测值和误差带送入到调度控制器中,预测值即为风电场未来的发电计划,误差带代表的功率范围即为对风电场配备的旋转热备用。本发明用于智能电网建设中。
申请公布号 CN102709926A 申请公布日期 2012.10.03
申请号 CN201210216019.0 申请日期 2012.06.28
申请人 哈尔滨工业大学 发明人 刘金福;苏鹏宇;李照忠;邢媛;万杰;于达仁
分类号 H02J3/24(2006.01)I 主分类号 H02J3/24(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 高媛
主权项 1.一种智能电网建设中基于相关向量机的旋转热备用的调度方法,其特征在于:所述方法是这样实现的:一、初始化;根据当地气候条件以及地理因素,进行当地风资源评估,根据风资源评估结果对计算控制器(3)进行初始化设置,所需设定的参数为预测结果的分布形式,所述分布形式设定为正态分布,将所述初始化设置结果传递到风功率相关向量机预测系统(3)中,初始化完毕;二、数据传输;风电场风功率采集模块(1)实时将风电场风功率的测量值进行采集,经数据预处理后将数据传递给风功率相关向量机预测系统(3)中;三、对未来时刻的风功率进行预测;风功率相关向量机预测系统(3)接收由风电场风功率采集模块(1)传递过来的数据,根据当前时刻的风功率值,对未来时刻的风功率进行预测,预测的结果为未来时刻的风功率值和风功率的误差带;风功率值和风功率的误差带的实现方式如下:集合<img file="FDA00001819471000011.GIF" wi="114" he="60" />由风功率历史数据组成,表示训练集中有n个历史样本,其中每一个历史样本x<sub>i</sub>有10个维数,由连续的10个时刻的风功率值组成,x<sub>i</sub>=[P<sub>i1</sub>,P<sub>i2</sub>,……P<sub>i10</sub>];t<sub>i</sub>为样本x<sub>i</sub>下一时刻的风功率值,即t<sub>i</sub>=P<sub>i11</sub>;对于本模型,即模型输入为x<sub>i</sub>,模型输出为t<sub>i</sub>,即根据当前10个时刻的风功率值,预测下一个时刻的风功率值,n为总的样本个数;将以上<img file="FDA00001819471000012.GIF" wi="114" he="53" />和<img file="FDA00001819471000013.GIF" wi="99" he="50" />作为训练样本集对相关向量机模型进行训练,<img file="FDA00001819471000014.GIF" wi="99" he="50" />为由t<sub>i</sub>组成的集合,由n个样本组成,t<sub>i</sub>与x<sub>i</sub>一一对应;三(一)、相关向量机模型的训练过程为:1)选定核函数K为高斯核函数,<img file="FDA00001819471000015.GIF" wi="563" he="130" />其中x<sub>1</sub>,x<sub>2</sub>为两个长度相等的向量,δ为核宽度,自行设定;2)预测方程为<maths num="0001"><![CDATA[<math><mrow><mi>t</mi><mrow><mo>(</mo><mi>x</mi><mo>;</mo><mi>w</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>w</mi><mi>i</mi></msub><mi>K</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>w</mi><mn>0</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中w=[w<sub>0</sub>,w<sub>1</sub>,……w<sub>n</sub>]表示权系数向量;3)设模型预测值满足正态分布形式,即p(t<sub>i</sub>)=N(t<sub>i</sub>|t(xi;w),σ<sup>2</sup>)(2)其中p表示分布函数,N表示正态分布;上述公式(2)表示预测结果为以公式(1)的结果t(x<sub>i</sub>;w)为均值,σ<sup>2</sup>为方差的正态分布形式;4)公式(1)和公式(2)中需要确定的参数仅为权系数向量w和正态分布的方差σ<sup>2</sup>,相关向量机法假设权系数向量w满足均值为0,方差为<img file="FDA00001819471000021.GIF" wi="44" he="107" />的正态分布,写成公式的形式即<maths num="0002"><![CDATA[<math><mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>w</mi><mi>i</mi></msub><mo>|</mo><msub><mi>&alpha;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>N</mi><mrow><mo>(</mo><msub><mi>w</mi><mi>i</mi></msub><mo>|</mo><mn>0</mn><mo>,</mo><msubsup><mi>&alpha;</mi><mi>i</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow></mrow></math>]]></maths>α=[α<sub>0</sub>,α<sub>1</sub>,......,α<sub>n</sub>]<sup>T</sup><maths num="0003"><![CDATA[<math><mrow><mi>p</mi><mrow><mo>(</mo><mi>w</mi><mo>|</mo><mi>&alpha;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Pi;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><mfrac><msub><mi>&alpha;</mi><mi>i</mi></msub><msqrt><mn>2</mn><mi>&pi;</mi></msqrt></mfrac><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mrow><msub><mi>&alpha;</mi><mi>i</mi></msub><msup><msub><mi>w</mi><mi>i</mi></msub><mn>2</mn></msup></mrow><mn>2</mn></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>经过上述变量替换后需要确定的参数仅为向量α=[α<sub>0</sub>,α<sub>1</sub>,……α<sub>n</sub>]以及方差σ<sup>2</sup>;5)使用数值解法求解α和σ<sup>2</sup>,先给出α与σ<sup>2</sup>猜测值,然后根据下式进行更新:<maths num="0004"><![CDATA[<math><mrow><msubsup><mi>&alpha;</mi><mi>i</mi><mi>NEW</mi></msubsup><mo>=</mo><mfrac><msub><mi>&gamma;</mi><mi>i</mi></msub><msup><msub><mi>&mu;</mi><mi>i</mi></msub><mn>2</mn></msup></mfrac></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msup><mrow><mo>(</mo><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>)</mo></mrow><mi>NEW</mi></msup><mo>=</mo><mfrac><msup><mrow><mo>|</mo><mo>|</mo><mi>t</mi><mo>-</mo><mi>&Phi;&mu;</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mrow><mi>n</mi><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>n</mi></munderover><msub><mi>&gamma;</mi><mi>i</mi></msub></mrow></mfrac></mrow></math>]]></maths>γ<sub>1</sub>=1-α<sub>i</sub>B(i,i)其中γi为中间变量,Φ为用核函数构成的矩阵,具体形式如下所示:<maths num="0006"><![CDATA[<math><mrow><mi>&Phi;</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd><mtd><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>x</mi><mi>n</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>,</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd><mtd><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>,</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>,</mo><msub><mi>x</mi><mi>n</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>n</mi></msub><mo>,</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd><mtd><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>n</mi></msub><mo>,</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mi>K</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>n</mi></msub><mo>,</mo><msub><mi>x</mi><mi>n</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>矩阵A为以向量α为对角线元素的方阵,矩阵B=(σ<sup>-2</sup>Φ<sup>T</sup>Φ+A)<sup>-1</sup>,B(i,i)是B中对角线上的第i项元素,向量μ=σ<sup>-2</sup>BΦ<sup>T</sup>t,T表示的是对矩阵进行转置操作,在足够多的更新之后,大部分的α<sub>i</sub>会趋近无限大,意即对应的w<sub>i</sub>为0;其它的α<sub>i</sub>会稳定趋近有限值,与之对应的x<sub>i</sub>就称为相关向量,在通过上述迭代过程得到了α和σ<sup>2</sup>以后,模型的训练过程结束;三(二)、进入预测环节;给定输入x<sub>*</sub>为当前时刻的风功率序列(包含十个时刻的风功率值),根据公式(1)和公式(2)模型输出风功率预测值t<sub>*</sub>和以及预测值的概率分布,由此得到了预测值和误差带;四、由风功率相关向量机预测系统(3)中得到的预测值和误差带送入到调度控制器(4)中,预测值即为风电场未来的发电计划,误差带代表的功率范围即为对风电场配备的旋转热备用(5),调度控制器(4)将控制旋转热备用(5)按照指令进行发电供用户使用。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号