发明名称 一种基于需求响应的配电网供电可靠性评估方法
摘要 本发明涉及一种基于需求响应的配电网供电可靠性评估方法。其方法步骤为:1建立基于分时电价的需求响应模型;2建立改进的基于需求响应的负荷优化模型;3对优化后的负荷曲线,利用二分法进行优化分级,求最佳分级数;4建立改进的基于二分法的聚类负荷分析可靠性评估算法;为计算方便,将聚类负荷区间变为单调递增的形势。本发明有效的解决了基于需求响应的配电网供电可靠性评估问题。
申请公布号 CN105405061A 申请公布日期 2016.03.16
申请号 CN201510224141.6 申请日期 2015.05.06
申请人 国家电网公司;国网河北省电力公司经济技术研究院 发明人 王颖;王涛;袁建普;康伟;任志刚;李树水;张章;张倩茅;朋;荆志朋;张丽洁;邵华;凌云鹏;贺春光;胡珀;马国真;韩文源;谢晓琳;秦伟;张欣悦;翟广心;张昭旭;孟繁华;彭冲
分类号 G06Q50/06(2012.01)I 主分类号 G06Q50/06(2012.01)I
代理机构 石家庄新世纪专利商标事务所有限公司 13100 代理人 张杰
主权项 一种基于需求响应的配电网供电可靠性评估方法,其特征在于包括如下步骤:第一步:建立基于分时电价的需求响应模型;(1) 时段划分作为峰谷分时电价的定价基础,峰谷时段划分必须能够正确的反映出实际负荷曲线的峰谷特性,并有效的与电力需求相衔接;峰谷时段的划分大致分为三步:首先,按照运行状况运用模糊聚类方法将系统运行划分为峰时段<i>T<sub>on</sub></i>、谷时段<i>T<sub>off</sub></i>、平时段<i>T<sub>mid</sub></i>三个基础时段;然后,通过对比各时段区间分界时刻点电力用户负荷曲线比重结构特征的变化情况,评估各时刻点的电力用户分时电价的响应程度;最后,融合各时刻点需求相应评估结果,修正调整基础时段划分区间;(2) 价格弹性系数电力价格弹性系数是电力需求量变动对价格变动的反应程度,即电力需求量的变化率与电力售价变化率的比值;如式(1)所示:<img file="dest_path_665737dest_path_image001.GIF" wi="82" he="83" />(1)式中,<img file="dest_path_520561dest_path_image002.GIF" wi="28" he="22" />代表电力需求改变量,<img file="dest_path_237981dest_path_image003.GIF" wi="17" he="22" />代表电力需求量,<img file="dest_path_69409dest_path_image004.GIF" wi="26" he="18" />代表用户电价改变量,<img file="dest_path_883781dest_path_image005.GIF" wi="17" he="18" />代表用户电价;执行峰谷分时电价后,用户在某一时段的电力需求不仅受当前电价水平的影响,而且还要受到其他时段的电价水平的影响;用户的电量在各时段之间是有替代效应的;因此根据峰、平、谷三时段的电量分配和电价设置,形成峰谷电力价格弹性矩阵,如下式(2)所示:<img file="dest_path_643927dest_path_image006.GIF" wi="303" he="204" />(2)式中,<img file="dest_path_848643dest_path_image007.GIF" wi="49" he="56" />、<img file="dest_path_47543dest_path_image008.GIF" wi="59" he="59" />和<img file="dest_path_654105dest_path_image009.GIF" wi="52" he="62" />分别表示峰、平、谷三个时段的电力电量改变量;<img file="dest_path_83687dest_path_image010.GIF" wi="62" he="77" />、<img file="dest_path_838016dest_path_image011.GIF" wi="55" he="60" />和<img file="dest_path_512711dest_path_image012.GIF" wi="49" he="61" />分别表示峰、平、谷三个时段的用户电价改变量;在峰谷电力弹性矩阵中,由于不同时段的电力商品之间为互补关系,因此,弹性矩阵中自弹性系数为正值,而交叉弹性系数为负值;各弹性系数数值的大小利用多元回归分析原理,通过统计大量历史数据值求得;(3) 引入分时电价后的负荷修正量分时电价的实施不仅使用户降低自身的用电量,最重要的还是将电力负荷从高电价时段转移到低电价时段,起到削峰填谷的作用,电力电量改变量如上式(2)所示,进而得到相应的负荷的改变量,如下式(3)所示:<img file="dest_path_973780dest_path_image013.GIF" wi="366" he="239" />(3)式中,<img file="dest_path_403624dest_path_image014.GIF" wi="26" he="25" /><img file="dest_path_582933dest_path_image015.GIF" wi="122" he="24" />为相应的小时数;<img file="dest_path_294274dest_path_image016.GIF" wi="39" he="26" />和<img file="dest_path_937745dest_path_image017.GIF" wi="34" he="24" />分别为实施峰谷分时电价前后的第<i>t</i>个小时的负荷;第二步:建立改进的基于需求响应的负荷优化模型;(1) 优化目标需建立一个考虑电价差的分时电价优化模型;建立一个考虑峰负荷值最小及峰谷负电量荷差值最小的多目标优化函数:峰负荷最小值为<img file="dest_path_210595dest_path_image018.GIF" wi="205" he="47" />(4)峰谷负荷差最小值为<img file="dest_path_877199dest_path_image019.GIF" wi="431" he="44" />(5)其中,<i>Q</i>是用户的用电量,<i>t</i>是该时段的持续时间,<img file="dest_path_956014dest_path_image020.GIF" wi="14" he="22" />代表谷时段电价和平时段电价的差值,<img file="dest_path_391674dest_path_image021.GIF" wi="16" he="18" />代表峰时段电价和平时段电价的差值;(2) 约束条件约束条件考虑电网及和电力市场的基本约束条件,主要包括供电方收益、用户侧受益、用电总量不变,谷成本低于边际成本;<img file="dest_path_599540dest_path_image022.GIF" wi="16" he="21" />供电方收益:实施需求响应后,必将起到削峰填谷的作用,这会大大降低供电方的发电成本;实施需求响应前后供电方收益关系如下式所示:<img file="dest_path_550178dest_path_image023.GIF" wi="167" he="43" />(6)式中,<img file="dest_path_370367dest_path_image024.GIF" wi="64" he="59" />表示利益转让系数,表示因为削峰作用引起的供电减少而接受的利益转移百分比;<img file="dest_path_722851dest_path_image025.GIF" wi="32" he="25" />表示分时电价实施后的供电方收益;<img file="dest_path_603082dest_path_image026.GIF" wi="17" he="25" />表示实施分时电价前的供电收益;<img file="dest_path_978700dest_path_image027.GIF" wi="16" he="21" />用户侧受益:由于需求响应的实施效果很大一部分取决于用户的参与情况,为此,需给予用户一定的补偿才能有效保证需求响应的效果,即用户参与需求响应后的平均价格是至少不能增加的:<img file="dest_path_399317dest_path_image028.GIF" wi="84" he="32" />(7)式中,<img file="dest_path_331542dest_path_image029.GIF" wi="48" he="35" />表示用户参与需求响应后的用户平均价格;<img file="dest_path_320358dest_path_image030.GIF" wi="21" he="25" />表示未实施需求响应是的用户平均价格;<img file="dest_path_183271dest_path_image031.GIF" wi="16" he="21" />用电总量不变为了方便计算,假设用户实行峰谷分时电价之后前后的用电量是保持不变的;<img file="dest_path_843798dest_path_image032.GIF" wi="204" he="38" />(8)式中,<img file="dest_path_905295dest_path_image033.GIF" wi="47" he="42" />,<img file="dest_path_127329dest_path_image034.GIF" wi="60" he="40" />和<img file="dest_path_477538dest_path_image035.GIF" wi="47" he="37" />分别表示峰、平、谷时段的用电量;<img file="dest_path_505537dest_path_image036.GIF" wi="21" he="25" />表示实施需求响应前的用户用电量;<img file="dest_path_624803dest_path_image037.GIF" wi="16" he="21" />边际成本约束:在谷时段时,电价不应该低于谷时段供电的边际成本:<img file="dest_path_250694dest_path_image038.GIF" wi="87" he="36" />(9)式中,<img file="dest_path_88200dest_path_image039.GIF" wi="37" he="33" />表示谷时段时的电价;<img file="dest_path_654311dest_path_image040.GIF" wi="45" he="37" />表示谷时段供电的边际成本;<img file="dest_path_893662dest_path_image041.GIF" wi="16" he="21" />价格序关系约束:峰时段的电价应该高于平时段的电价,平时段的电价应该高于谷时段的电价,即:<img file="dest_path_254236dest_path_image042.GIF" wi="128" he="33" />(10)式中,<img file="dest_path_579038dest_path_image043.GIF" wi="26" he="25" />,<img file="dest_path_385058dest_path_image044.GIF" wi="37" he="25" />和<img file="dest_path_275654dest_path_image045.GIF" wi="29" he="26" />分别表示峰、平、谷时段的电价;该优化问题的求解是采用模糊优化理论对峰谷分时电价模型进行连续优化,综合考虑各个单目标最优解和多目标满意解之间的相互关系,寻求得到使各个目标都处于较优状态的最客观的峰谷分时电价;第三步:对优化后的负荷曲线,利用二分法进行优化分级,求最佳分级数;(1) 模型的显著性水平最优的分级、聚类方案应该达到如下的效果:各个负荷级之间差值较大、同一负荷级内各负荷点之间的差值较小;根据概率统计理论,用显著性水平来表征分级、聚类方案的好坏;模型的显著性水平越高,说明分级、聚类的效果越好;定义如下指标:模型的显著性水平<img file="dest_path_10392dest_path_image046.GIF" wi="22" he="22" />:<img file="dest_path_556911dest_path_image047.GIF" wi="54" he="48" />(11)其中,<img file="dest_path_668086dest_path_image048.GIF" wi="23" he="27" />为负荷级间距离的最小值;<img file="dest_path_678767dest_path_image049.GIF" wi="26" he="27" />为同一负荷级内各个负荷点的标准差的最大值;<img file="dest_path_82942dest_path_image048.GIF" wi="23" he="27" />、<img file="dest_path_116757dest_path_image049.GIF" wi="26" he="27" />根据式(12)、(13)计算得到;假设年负荷曲线的<i>N</i>个小时负荷点被分为<img file="dest_path_766044dest_path_image050.GIF" wi="24" he="21" />个负荷级,第<i>i</i>负荷级的负荷值为<img file="dest_path_631232dest_path_image051.GIF" wi="41" he="41" />,负荷点数目<img file="dest_path_707772dest_path_image052.GIF" wi="49" he="36" />、负荷点集合为<img file="dest_path_727418dest_path_image053.GIF" wi="31" he="44" />,则:<img file="dest_path_977134dest_path_image054.GIF" wi="240" he="39" />(12)其中,<img file="dest_path_900091dest_path_image055.GIF" wi="122" he="32" />,它表示第<img file="dest_path_147533dest_path_image056.GIF" wi="33" he="21" />负荷级与第<img file="dest_path_155940dest_path_image057.GIF" wi="12" he="21" />负荷级之间的距离,<img file="dest_path_209347dest_path_image058.GIF" wi="26" he="27" />取其中的最小值;<img file="dest_path_485345dest_path_image059.GIF" wi="250" he="42" />(13)其中,<img file="dest_path_903688dest_path_image060.GIF" wi="253" he="67" />,它表示第<img file="dest_path_196129dest_path_image061.GIF" wi="14" he="24" />负荷级的级内标准差,<img file="dest_path_990910dest_path_image062.GIF" wi="23" he="25" />取其中的最大值;(2) 模型的最佳分级数分级、聚类的效果用模型的显著性水平来评价,把模型显著性水平最高的分级、聚类方案作为最终负荷模型,对应的分级数为模型的最佳分级数;为了快速地确定最佳分级数<img file="dest_path_622879dest_path_image063.GIF" wi="22" he="21" />,本文采用二分法的原理进行查找:假设年负荷曲线的N个小时负荷点的最小分级数为<img file="dest_path_274440dest_path_image064.GIF" wi="22" he="21" />,最大分级数为<img file="dest_path_755975dest_path_image065.GIF" wi="25" he="21" />;显然最佳分级数<img file="dest_path_88867dest_path_image066.GIF" wi="129" he="37" />;假如N个小时负荷点各不相同,不难看出<img file="dest_path_575344dest_path_image064.GIF" wi="22" he="21" />和<img file="dest_path_663385dest_path_image065.GIF" wi="25" he="21" />的初始值分别为1和N;当分级数取中间值<img file="dest_path_868102dest_path_image067.GIF" wi="151" he="33" />时,计算模型的显著性水平,如果显著性水平<img file="dest_path_503220dest_path_image068.GIF" wi="57" he="32" />,比较<img file="dest_path_906520dest_path_image069.GIF" wi="64" he="35" />、<img file="dest_path_103146dest_path_image070.GIF" wi="52" he="28" />和<img file="dest_path_795158dest_path_image071.GIF" wi="56" he="29" />之间的关系,若<img file="dest_path_532170dest_path_image072.GIF" wi="131" he="31" />且<img file="dest_path_993238dest_path_image073.GIF" wi="139" he="33" />则停止查找,<img file="dest_path_859301dest_path_image074.GIF" wi="54" he="20" />;若模型的显著性水平处于递增的趋势即<img file="dest_path_835347dest_path_image075.GIF" wi="245" he="37" />,表明最佳分级数<img file="dest_path_313733dest_path_image076.GIF" wi="121" he="32" />,在<img file="dest_path_894887dest_path_image077.GIF" wi="103" he="44" />中继续查找;若模型的显著性水平处于递减的趋势即<img file="dest_path_167737dest_path_image078.GIF" wi="213" he="34" />,表明最佳分级数<img file="dest_path_332877dest_path_image079.GIF" wi="109" he="31" />,在<img file="dest_path_411691dest_path_image080.GIF" wi="89" he="41" />中继续查找,直到找到最佳分级数为止;第四步:建立改进的基于二分法的聚类负荷分析可靠性评估算法;对于考虑负荷变化的配电网供电可靠性评估算法常采用聚类分析法;其计算过程如下:安排计及负荷曲线的可靠性评估算法中应按照每枚举一个故障时,就计算出其在不同负荷水平下的潮流约束;在考虑利用聚类负荷曲线<img file="dest_path_847351dest_path_image081.GIF" wi="177" he="34" />情况下,共轭负荷矩的二次潮流估算方法;假设已在第<i>k</i>个负荷水平<img file="dest_path_291102dest_path_image082.GIF" wi="35" he="25" />下计算节点m的电压,按照共轭负荷矩定义,得如下式子:<img file="dest_path_507320dest_path_image083.GIF" wi="356" he="89" />(14)式中,<img file="dest_path_327508dest_path_image084.GIF" wi="42" he="25" />表示在负荷水平<img file="dest_path_116211dest_path_image085.GIF" wi="37" he="27" />下,节点<i>m</i>在联络开关闭合前的共轭负荷矩;当全网络的负荷水平按相同比例增加为<img file="dest_path_58759dest_path_image086.GIF" wi="56" he="25" />(<img file="dest_path_434377dest_path_image087.GIF" wi="111" he="25" />,C为常数)时,有共轭负荷矩可线性叠加的特点,这时的共轭负荷矩表示为<img file="dest_path_792677dest_path_image088.GIF" wi="126" he="38" />;这时节点<i>m</i>的电压表示为:<img file="dest_path_265246dest_path_image089.GIF" wi="211" he="99" />(15)对于配电网有:<img file="dest_path_50800dest_path_image090.GIF" wi="66" he="74" />〈〈<img file="dest_path_412249dest_path_image091.GIF" wi="35" he="50" />,且<img file="dest_path_636557dest_path_image092.GIF" wi="53" he="50" />,若<img file="dest_path_901316dest_path_image093.GIF" wi="68" he="17" />(<img file="dest_path_123350dest_path_image094.GIF" wi="33" he="17" />)时,认为<img file="dest_path_535877dest_path_image095.GIF" wi="108" he="53" />,有:<img file="dest_path_235979dest_path_image096.GIF" wi="153" he="72" />(16)用式(16)计算在联络开关闭合前,与联络线直接相连的各节点在下一个负荷水平的电压,并用<img file="dest_path_853780dest_path_image097.GIF" wi="41" he="25" />实现送端馈线潮流估计和受端馈线实现下一个负荷水平的潮流估计和负荷削减策略的确定;为计算方便,将聚类负荷区间变为单调递增的形势。
地址 100031 北京市西城区长安西街甲86号