发明名称 一种电网连锁故障的抑制方法
摘要 一种电网连锁故障的抑制方法,属于电力系统安全防护技术领域。该方法通过计算出对过载支路影响最大的节点,确定出切机切负荷量范围,并对其进行多次优化,按照每一次优化后的切机切负荷量对电网网络实施切机和切负荷操作,并检测下一时刻过载支路的电流值是否在安全范围内,以此方法连续循环进行切机切负荷量优化和实施切负荷操作,直至过载支路的电流值达到安全范围或者切机切负荷量优化次数超过最大次数停止,该法可在支路过载时优化切机切负荷量,降低损失,同时有效抑制连锁故障的发生。
申请公布号 CN104466953A 申请公布日期 2015.03.25
申请号 CN201410747329.4 申请日期 2014.12.09
申请人 东北大学;国网辽宁省电力有限公司电力科学研究院;国家电网公司 发明人 张化光;高凯;杨珺;张艳军;王馨漪;孙秋野;刘鑫蕊;郝松;冯健;杨东升;汪刚
分类号 H02J3/00(2006.01)I;H02H7/28(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 沈阳东大知识产权代理有限公司 21109 代理人 梁焱
主权项 一种电网连锁故障的抑制方法,其特征在于:包括如下步骤:步骤1:确定电网网络的拓扑结构并统计网络节点总数b和网络支路总数s、确定因故障而被切除的网络支路、发生潮流转移的网络支路;步骤2:建立网络支路电流的节点重要性贡献矩阵并计算该矩阵;将该电网网络中的发电机、负荷都作为节点注入电流表示,则网络支路电流的节点重要性贡献矩阵为:<img file="FDA0000628313440000011.GIF" wi="1422" he="342" />式中,R(γ)=YA<sup>T</sup>B<sup>‑1</sup>为b*s阶矩阵,Y为网络的支路导纳矩阵,A为网络的节点关联矩阵,B为网络的节点导纳矩阵;w<sub>j</sub>为节点的权重系数,j=1,2,…,s,c<sub>j</sub>为节点的代价系数;矩阵I<sub>m</sub>(λ)表示的是支路电流与节点注入电流的关系重要度,即节点注入电流对支路电流影响的重要程度,矩阵I<sub>m</sub>(λ)中第一行元素λ<sub>1j</sub>、第k行元素λ<sub>kj</sub>、第b行元素λ<sub>bj</sub>分别表示的是节点1至节点s的注入电流分别对第一条支路电流影响的重要程度、对第k条支路电流影响的重要程度、对第b条支路电流影响的重要程度;步骤3:根据步骤2中网络支路电流的节点重要性贡献矩阵的计算结果给出的节点对发生潮流转移的支路电流影响的重要程度,确定需要进行切机的节点和切负荷的节点;假设发生潮流转移的支路为第k条支路,则将矩阵I<sub>m</sub>(λ)中第k行元素λ<sub>kj</sub>按照数值大小进行排序后:从小到大选取n个相应的节点l<sub>i</sub>(i=1,...,n),即负荷节点,作为待切负荷节点;同时,从大到小选取与待切负荷节点相同个数的节点g<sub>i</sub>(i=1,...,n),即发电机节点,作为待切机节点;步骤4:确定切机切负荷量范围;在电网网络中切机量与切负荷量应相等,且满足以下公式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>&Delta;</mi><msub><mi>I</mi><mi>k</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>&Delta;</mi><msub><mi>P</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>&gamma;</mi><msub><mi>kg</mi><mi>i</mi></msub></msub><mo>|</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>j</mi><mi>tan</mi><msub><mi>&theta;</mi><msub><mi>g</mi><mi>i</mi></msub></msub><mo>)</mo></mrow><mo>/</mo><msubsup><mover><mi>U</mi><mo>.</mo></mover><msub><mi>g</mi><mi>i</mi></msub><mo>*</mo></msubsup><mo>|</mo><mo>-</mo><msub><mi>&gamma;</mi><msub><mi>kl</mi><mi>i</mi></msub></msub><mo>|</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>j</mi><mi>tan</mi><msub><mi>&theta;</mi><msub><mi>l</mi><mi>i</mi></msub></msub><mo>)</mo></mrow><mo>/</mo><msubsup><mover><mi>U</mi><mo>.</mo></mover><msub><mi>l</mi><mi>i</mi></msub><mo>*</mo></msubsup><mo>|</mo><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>P</mi><mi>i</mi></msub><mo>=</mo><mi>&Delta;</mi><msub><mi>P</mi><msub><mi>l</mi><mi>i</mi></msub></msub><mo>=</mo><mi>&Delta;</mi><msub><mi>P</mi><msub><mi>g</mi><mi>i</mi></msub></msub><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>P</mi><mi>i</mi></msub><mo>&le;</mo><msup><msub><mi>P</mi><msub><mi>l</mi><mi>i</mi></msub></msub><mn>0</mn></msup><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>P</mi><mi>i</mi></msub><mo>&le;</mo><msup><msub><mi>P</mi><msub><mi>g</mi><mi>i</mi></msub></msub><mn>0</mn></msup><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000628313440000012.GIF" wi="1409" he="342" /></maths>式中,<img file="FDA0000628313440000013.GIF" wi="331" he="89" />ΔI<sub>k</sub>为第k条支路上电流变化量的有效值,I<sub>k</sub>为第k条支路的电流有效值;I<sub>kset</sub>为第k条支路的保护启动所对应的最小电流值,S(S<1)为保证支路不发生连锁跳闸的可靠系数;<img file="FDA0000628313440000027.GIF" wi="74" he="60" />为R(γ)矩阵中第k行l<sub>i</sub>列,表示第k条支路的电流与负荷节点l<sub>i</sub>的注入电流之间的数值关系;<img file="FDA0000628313440000028.GIF" wi="71" he="60" />为R(γ)矩阵中第k行g<sub>i</sub>列,表示第k条支路的电流与发电机节点g<sub>i</sub>的注入电流之间的数值关系;<img file="FDA0000628313440000029.GIF" wi="54" he="71" />与<img file="FDA00006283134400000210.GIF" wi="52" he="71" />分别表示负荷节点l<sub>i</sub>的负荷功率因数角和发电机节点g<sub>i</sub>的发电机功率因数角;<img file="FDA0000628313440000021.GIF" wi="67" he="80" />与<img file="FDA0000628313440000022.GIF" wi="66" he="80" />表示负荷节点l<sub>i</sub>与发电机节点g<sub>i</sub>的电压向量的共轭;<img file="FDA0000628313440000023.GIF" wi="201" he="76" />分别为负荷节点l<sub>i</sub>、发电机节点g<sub>i</sub>的切除量;全部容量<img file="FDA0000628313440000024.GIF" wi="54" he="83" />为负荷节点l<sub>i</sub>的全部容量,<img file="FDA0000628313440000025.GIF" wi="78" he="78" />为发电机节点g<sub>i</sub>的全部容量;步骤5:对步骤4得到的切机切负荷量范围进行优化,得到切机切负荷量的优化值;采用公式(2)优化目标函数对步骤4得到的切机切负荷量范围进行优化;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><munder><mi>min</mi><mrow><mi>&Delta;</mi><msub><mi>P</mi><mi>i</mi></msub></mrow></munder><msub><mi>J</mi><mi>T</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>L</mi><mrow><mo>(</mo><mi>I</mi><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mi>&Delta;</mi><msub><mi>P</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mi>F</mi><mrow><mo>(</mo><mi>I</mi><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000628313440000026.GIF" wi="1920" he="122" /></maths>s.t.I(t+1)=f(I(t),ΔP<sub>1</sub>(t),...,ΔP<sub>n</sub>(t))ΔP<sub>i</sub>(t+1)∈[0,ΔP<sub>imax</sub>](i=1,...,n)式中,J为优化目标函数,I(t)和I(t+1)分别为第k条支路在t时刻的电流值和进行切机切负荷后t+1时刻的电流值;L(I,ΔP<sub>i</sub>)表示切机切负荷量优化控制的预定目标,一般取为范数形式:L(I,ΔP<sub>i</sub>)=I<sup>T</sup>QI+ΔP<sub>i</sub><sup>T</sup>RΔP<sub>i</sub>,其中,R和Q是正定对称加权矩阵;f表示控制量与输入量的函数关系;ΔP<sub>i</sub>(t+1)表示的是t+1时刻切机切负荷量;ΔP<sub>imax</sub>表示最大允许切除量;F(I)为终端惩罚函数,表示时刻t的状态与设定值的偏差,一般形式为F(I)=I<sup>T</sup>PI,其中,P为终端加权矩阵;步骤6:按照步骤5得到的切机切负荷量优化值对当前时刻的电网网络实施切机和切负荷;步骤7:检测下一时刻发生潮流转移的网络支路的电流值并记录该电流值,并判断其电流是否进入安全范围内,否,则转至步骤8,是,则转至步骤10;步骤8:判断切机切负荷量优化次数是否达到事先设定的次数,否,则转至步骤2,是,则转至步骤9;步骤9:输出发生潮流转移的网络支路的电流值;步骤10:结束。
地址 110819 辽宁省沈阳市和平区文化路3号巷11号