发明名称 基于场景削减的配电网抗负荷波动的无功功率控制方法
摘要 本发明涉及一种基于场景削减的配电网抗负荷波动的无功功率控制方法,属于电力系统优化运行控制领域。首先根据电能计量数据得到三相不确定负荷有功功率预测误差的区间以及概率密度函数,将其分为多个离散区间,以各区间的中间值代表不确定负荷节点三相有功功率预测误差值,面积代表与误差值相对应的概率,通过同步回代缩减法对视在功率组合进行削减,将所有削减后的组合代入到配电网无功优化控制模型中,求解适用于所有组合的最优无功功率出力决策值。本方法考虑了无量测负荷功率的不确定性,并减小负荷组合的规模,缩减了无功优化控制的求解时间,使得优化求解的无功功率出力决策值对负荷的预测误差具有适应性,有利于配电网的电压优化控制。
申请公布号 CN104410078B 申请公布日期 2016.08.17
申请号 CN201410598663.8 申请日期 2014.10.29
申请人 国网山东省电力公司潍坊供电公司;国家电网公司;清华大学 发明人 吴琥;吴文传;马金亮;张伯明;张帆;王永杰;杨晶晶;郭庆来
分类号 H02J3/18(2006.01)I 主分类号 H02J3/18(2006.01)I
代理机构 北京清亦华知识产权代理事务所(普通合伙) 11201 代理人 罗文群
主权项 一种基于场景削减的配电网抗负荷波动的无功功率控制方法,其特征在于该方法包括以下步骤:(1)设配电网中的不确定负荷节点数为N,并设配电网中不确定负荷节点的实际有功功率预测区间为<img file="FDA0000933230050000011.GIF" wi="955" he="87" />其中<img file="FDA0000933230050000012.GIF" wi="198" he="79" />分别为N个不确定负荷节点中实际有功功率预测区间的下限值和上限值,从配电网的电能计量历史数据中得到不确定负荷节点有功功率的预测期望值分别为<img file="FDA0000933230050000013.GIF" wi="467" he="78" />不确定负荷节点预测有功功率的标准差分别为σ<sub>L1</sub>,σ<sub>L2</sub>,…,σ<sub>Ln</sub>,…,σ<sub>LN</sub>,配电网的不确定负荷节点的三相视在功率S<sub>L1</sub>,S<sub>L2</sub>,…,S<sub>Ln</sub>,…,S<sub>LN</sub>分别为:<maths num="0001"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><msub><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>L</mi><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>jP</mi><mrow><mi>L</mi><mn>1</mn></mrow></msub><mo>*</mo><msub><mi>tan&theta;</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>S</mi><mrow><mi>L</mi><mn>2</mn></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>L</mi><mn>2</mn></mrow></msub><mo>+</mo><msub><mi>jP</mi><mrow><mi>L</mi><mn>2</mn></mrow></msub><mo>*</mo><msub><mi>tan&theta;</mi><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mrow><msub><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>L</mi><mi>n</mi></mrow></msub><mo>+</mo><msub><mi>jP</mi><mrow><mi>L</mi><mi>n</mi></mrow></msub><mo>*</mo><msub><mi>tan&theta;</mi><mi>n</mi></msub></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mrow><msub><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow></msub><mo>=</mo><msub><mi>P</mi><mrow><mi>L</mi><mi>N</mi></mrow></msub><mo>+</mo><msub><mi>jP</mi><mrow><mi>L</mi><mi>N</mi></mrow></msub><mo>*</mo><msub><mi>tan&theta;</mi><mi>N</mi></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000933230050000014.GIF" wi="494" he="440" /></maths>其中,P<sub>L1</sub>,P<sub>L2</sub>,…,P<sub>Ln</sub>,…,P<sub>LN</sub>分别为各不确定负荷节点三相有功功率的实际值,tanθ<sub>1</sub>,tanθ<sub>2</sub>,…,tanθ<sub>n</sub>,…,tanθ<sub>N</sub>分别为各不确定负荷节点三相负荷的功率因数,从与不确定负荷节点相连的设备获取;将上述各不确定负荷节点实际有功功率的预测区间<img file="FDA0000933230050000015.GIF" wi="187" he="93" />转变为[‑1,1]的形式如下:<img file="FDA0000933230050000016.GIF" wi="1254" he="141" />其中,P<sub>Ln</sub>为N个不确定负荷节点中第n个不确定负荷节点的实际有功功率,与[‑1,1]区间相对应的不确定负荷节点实际有功功率的等效标准差<img file="FDA0000933230050000017.GIF" wi="70" he="79" />为:<img file="FDA0000933230050000018.GIF" wi="854" he="199" />(2)利用核密度估计方法,根据配电网电能计量装置的电量数据,得到各不确定负荷节点的三相有功功率的预测误差分布,求得各不确定负荷节点的三相有功功率预测误差等效区间的正态分布函数为:<img file="FDA0000933230050000021.GIF" wi="732" he="230" />其中,x<sub>n</sub>为N个不确定负荷节点中第n个不确定负荷节点的三相有功功率的预测误差等效区间内的实际值,μ<sub>n</sub>为N个不确定负荷节点中第n个不确定负荷节点的三相有功功率的预测误差等效区间内的期望值,取值为0,<img file="FDA0000933230050000022.GIF" wi="134" he="109" />为N个不确定负荷节点中第n个三相有功功率的预测误差等效区间内的方差;(3)在[‑1,1]区间,将上述正态分布函数分成M个区间,每一区间的中间值<img file="FDA0000933230050000023.GIF" wi="123" he="63" />为该不确定负荷节点三相有功功率预测误差值,每一区间的面积为该三相有功功率预测误差值出现的第n个不确定负荷节点的第m个区间的初始概率,记为<img file="FDA0000933230050000024.GIF" wi="132" he="92" />对M个区间的初始概率进行归一化处理,得到每一区间的面积为该三相有功功率预测误差值出现的概率:<img file="FDA0000933230050000025.GIF" wi="1086" he="222" />其中,<img file="FDA0000933230050000026.GIF" wi="101" he="71" />为第n个不确定负荷节点的第m个区间的概率;(4)根据上述由正态分布函数M个区间的中间值<img file="FDA0000933230050000027.GIF" wi="123" he="63" />得到的不确定负荷节点三相有功功率预测误差值,通过下式得到不确定负荷节点的M×N个三相有功功率的可能实际值<img file="FDA0000933230050000028.GIF" wi="67" he="61" />为:<img file="FDA0000933230050000029.GIF" wi="1220" he="143" />通过下式得到不确定负荷节点的M×N个三相视在功率的可能实际值<img file="FDA00009332300500000210.GIF" wi="67" he="63" />为:<img file="FDA00009332300500000211.GIF" wi="1661" he="262" />其中,<img file="FDA00009332300500000212.GIF" wi="62" he="61" />为第n个不确定负荷节点第m个区间的三相视在功率可能实际值的相对误差对应的比例,<img file="FDA00009332300500000213.GIF" wi="79" he="79" />为不确定负荷节点的三相视在功率的预测期望值,<img file="FDA00009332300500000214.GIF" wi="484" he="80" />tanθ<sub>n</sub>为第n个不确定负荷节点三相负荷的功率因数,从与不确定负荷节点相连的设备获取,j为虚数单位;(5)将N个不确定负荷节点的M个三相视在功率的可能实际值进行排列组合,得到第一组合序列:<maths num="0002"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>A</mi><mo>=</mo><mo>{</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mn>1</mn></msubsup><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mn>2</mn></msubsup><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mi>m</mi></msubsup><mo>;</mo><mn>...</mn><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mi>M</mi></msubsup><mo>;</mo></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mrow><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mn>1</mn></msubsup><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mn>1</mn></msubsup><mo>;</mo><mn>...</mn><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mn>1</mn></msubsup><mo>;</mo><mn>...</mn><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mn>1</mn></msubsup><mo>;</mo></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mn>2</mn></msubsup><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mn>2</mn></msubsup><mo>;</mo><mn>...</mn><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mn>2</mn></msubsup><mo>;</mo><mn>...</mn><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mn>2</mn></msubsup><mo>;</mo></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mrow><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mi>m</mi></msubsup><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mi>m</mi></msubsup><mo>;</mo><mn>...</mn><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mi>m</mi></msubsup><mo>;</mo><mn>...</mn><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mi>m</mi></msubsup><mo>;</mo></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mrow><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mi>M</mi></msubsup><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mi>M</mi></msubsup><mo>;</mo><mn>...</mn><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mi>M</mi></msubsup><mo>;</mo><mn>...</mn><mo>;</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mn>1</mn></mrow><mn>1</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>n</mi></mrow><mn>2</mn></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>S</mi><mrow><mi>L</mi><mi>N</mi></mrow><mi>M</mi></msubsup><mo>;</mo><mo>}</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>M</mi><mo>,</mo><mi>n</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>N</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000933230050000031.GIF" wi="1741" he="758" /></maths>通过下式,计算得到上述第一组合序列中与每个组合相对应的三相视在功率的可能实际值出现的概率P<sub>b</sub>为:<maths num="0003"><math><![CDATA[<mrow><msub><mi>P</mi><mi>b</mi></msub><mo>=</mo><munderover><mo>&Pi;</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>P</mi><mrow><mi>l</mi><mi>e</mi><mi>v</mi><mi>e</mi><mi>l</mi><mo>-</mo><mi>m</mi></mrow><mi>n</mi></msubsup></mrow>]]></math><img file="FDA0000933230050000032.GIF" wi="260" he="127" /></maths>其中,<img file="FDA0000933230050000033.GIF" wi="102" he="70" />为第n个不确定负荷节点的第m个区间的概率,∏为连乘计算;对与上述组合序列中的各组合相对应的不确定负荷节点三相视在功率取模值,将N个不确定负荷节点的M个三相视在功率的模值进行排列组合,得到第二组合序列为:B=|A| m=1,2,…,M,n=1,2,…,N(6)采用同步回代缩减法,对上述第二组合序列进行削减,得到削减后的组合序列,包括以下步骤:(6‑1)从上述第二组合序列中删除两个三相视在功率模值为最大和最小的组合,将删除最大和最小功率模值后的组合序列作为初始集合,记为D,将两个功率模值为最大和最小的组合组成一个第三序列组合序列,记为J;(6‑2)分别计算上述初始集合D中所有任意组合对ε<sub>b</sub>和ε<sub>b′</sub>之间的距离:DT<sub>b,b′</sub>=||ε<sub>b</sub>‑ε<sub>b′</sub>||<sub>2</sub>,b=1,2,…,N<sub>s</sub>,b′=1,2,…,N<sub>s</sub>其中,N<sub>s</sub>为初始集合D中组合数目,ε<sub>b</sub>为初始集合D中的任意一个组合,ε<sub>b′</sub>为初始集合D中的所有组合,定义组合对ε<sub>b</sub>和ε<sub>b′</sub>中的任意一个组合为主体组合时,另一个组合则为配对组合,ε<sub>b</sub>和ε<sub>b′</sub>均为向量,||·||<sub>2</sub>表示向量的2范数,即向量中各个元素平方和再开根号的值,DT<sub>b,b′</sub>表示组合ε<sub>b</sub>与ε<sub>b′</sub>之间的距离;(6‑3)利用步骤(6‑2)的方法,计算初始集合D中的本体组合ε<sub>k</sub>与所有配对组合ε<sub>b′</sub>之间的距离DT<sub>k,b′</sub>,取所有距离DT<sub>k,b′</sub>中的最小值,与该最小值相对应的配对组合ε<sub>r</sub>为初始集合D中距离本体组合ε<sub>k</sub>最近的组合,即DT<sub>k,r</sub>=minDT<sub>k,b′</sub>,其中b′,k∈D,b′≠k;(6‑4)根据步骤(5)得到的每个组合的三相视在功率的可能实际值出现的概率P<sub>s</sub>和上述最小距离DT<sub>k,b′</sub>,计算概率距离PD<sub>k,r</sub>:PD<sub>k,r</sub>=P<sub>k</sub>×DT<sub>k,r</sub>,k∈D,重复本步骤,分别计算得到初始集合D中各组合的概率距离,取所有概率距离中的最小值,根据该最小值,得到初始集合D中与该最小值相对应的组合对ε<sub>u</sub>、ε<sub>v</sub>,并设定ε<sub>u</sub>为主体组合,ε<sub>v</sub>为配对组合;(6‑5)从步骤(6‑2)初始集合D中削减上述主体组合ε<sub>u</sub>,并将ε<sub>u</sub>添加到第三组合序列J中,将主体组合ε<sub>u</sub>的概率P<sub>u</sub>加到距离最近的配对组合ε<sub>v</sub>的概率P<sub>v</sub>上:D=D‑{ε<sub>u</sub>}J=J+{ε<sub>u</sub>}P<sub>v</sub>=P<sub>v</sub>+P<sub>u</sub>(6‑6)设定一个概率距离的阈值PD<sub>min</sub>,将步骤(6‑4)中的概率距离最小值与设定的概率距离阈值进行比较,若概率距离最小值大于或等于设定的概率距离阈值,则重复步骤(6‑2)~(6‑5),若概率距离最小值小于设定的概率距离阈值,则得到削减后的最终集合U,该最终集合U中组合的概率值记为Pr<sub>w</sub>;(7)设削减后最终集合U中的组合数为N<sub>r</sub>,任意组合w的概率为Pr<sub>w</sub>,则最终集合U中的所有组合组成的组合序列以矩阵Q形式的表达式为:Q=[|S<sub>L1</sub>|,|S<sub>L2</sub>|,…,|S<sub>Ln</sub>|,…,|S<sub>LN</sub>|]上式中,|S<sub>Ln</sub>|为削减后最终集合U中第n个不确定负荷节点的三相视在功率模值组成的列向量,该列向量的维度为N<sub>r</sub>×1,Q为N<sub>r</sub>×N的矩阵,矩阵Q中的行代表N个不确定负荷节点三相视在功率模值的组合,矩阵Q中的列代表与N<sub>r</sub>个组合相对应的不确定负荷节点视在功率的模值;将上述最终矩阵Q中的第n个不确定负荷节点的三相视在功率模值组成的列向量与N个不确定负荷中第n个负荷的三相视在功率预测期望值的模值<img file="FDA0000933230050000051.GIF" wi="97" he="102" />的比值记为T:即<img file="FDA0000933230050000052.GIF" wi="1318" he="111" />其中,<img file="FDA0000933230050000053.GIF" wi="45" he="61" />为N个不确定负荷中第n个不确定负荷节点三相视在功率总和的模值与预测期望值模值的比值;由于配电网不确定负荷节点的A、B、C三相视在功率预测的相对误差值对应的比例与三相总和视在功率预测的相对误差值对应的比例相同,得到第n个不确定负荷节点的A、B、C三相视在功率模值与A、B、C三相视在功率预测期望值模值比值的三个列向量<img file="FDA0000933230050000054.GIF" wi="91" he="62" /><img file="FDA0000933230050000055.GIF" wi="184" he="63" />分别为:<maths num="0004"><math><![CDATA[<mrow><msubsup><mi>t</mi><mrow><mi>L</mi><mi>A</mi></mrow><mi>n</mi></msubsup><mo>=</mo><msubsup><mi>t</mi><mrow><mi>L</mi><mi>B</mi></mrow><mi>n</mi></msubsup><mo>=</mo><msubsup><mi>t</mi><mrow><mi>L</mi><mi>C</mi></mrow><mi>n</mi></msubsup><mo>=</mo><msubsup><mi>t</mi><mi>L</mi><mi>n</mi></msubsup></mrow>]]></math><img file="FDA0000933230050000056.GIF" wi="374" he="63" /></maths>并得到第n个不确定负荷节点的A、B、C三相视在功率的实际值<img file="FDA0000933230050000057.GIF" wi="67" he="62" />为:<img file="FDA0000933230050000058.GIF" wi="1037" he="111" />其中,<img file="FDA0000933230050000059.GIF" wi="69" he="95" />和<img file="FDA00009332300500000510.GIF" wi="85" he="93" />分别代表第n个不确定负荷节点的A、B、C三相有功功率和无功功率的期望值,<img file="FDA00009332300500000511.GIF" wi="68" he="70" />和<img file="FDA00009332300500000512.GIF" wi="77" he="70" />为第n个不确定负荷节点的A、B、C三相有功功率和无功功率的实际值;进而得到配电网中各不确定负荷节点的A、B、C实际运行的有功功率和无功功率;(8)将上述得到的配电网中各不确定负荷节点的A、B、C实际运行的有功功率和无功功率用于配电网无功功率优化控制模型,计算得到配电网中各无功设备的无功功率,实现配电网抗负荷波动下的无功功率控制。
地址 261021 山东省潍坊市东风西街425号