发明名称 电网自动电压控制系统间的协调量设计方法及其实施方法
摘要 本发明涉及电力系统自动电压控制技术领域,具体涉及一种电网自动电压控制系统间的协调量设计方法及其实施方法,所述协调量为电压区域控制偏差,由不同区域、不同层次的自动电压控制系统间的边界节点电压和外区域注入本区域无功的加权和构成。本发明通过该协调量能够判断无功扰动发生在区域内还是区域外,并衡量本区域控制动作对相邻区域无功影响的方向和大小,每个控制区只负责控制本区域内的无功扰动,使考核公平,由此明确区域控制责任、减小耦合,达到协调控制的目的。
申请公布号 CN101944743B 申请公布日期 2012.12.26
申请号 CN201010291701.7 申请日期 2010.09.21
申请人 东北电网有限公司;清华大学 发明人 王钢;梅生伟;苏辛一;张雪敏;陶家琪;何光宇;贾伟;高德宾;刘家庆;孟令愚;侯凯元;杨宁;李泽宇;邵广惠;徐兴伟;马新;彭晓洁;郗郅;夏德明;岳涵
分类号 H02J3/00(2006.01)I;H02J3/04(2006.01)I;H02J3/16(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 北京路浩知识产权代理有限公司 11002 代理人 王莹
主权项 1.一种电网自动电压控制系统间的协调量设计方法,其特征在于,按如下步骤设计电压区域控制偏差VACE,所述VACE为不同区域、不同层次的自动电压控制系统间的由边界节点电压和外区域注入本区域无功的加权和构成的协调量:S1、对于一电压控制区域S,列出电压对无功功率的灵敏度方程;S2、根据所述灵敏度方程形成支路串联电纳矩阵<img file="FDA00001824052800011.GIF" wi="231" he="51" />支路并联电纳矩阵<img file="FDA00001824052800012.GIF" wi="203" he="39" />和边界节点并联电纳矩阵<img file="FDA00001824052800013.GIF" wi="234" he="54" />所述支路串联电纳矩阵为各支路的串联电纳组成的矩阵;支路并联电纳矩阵为各支路的并联电纳组成的矩阵;边界节点并联电纳矩阵为边界节点间的并联电纳组成的矩阵;R表示实数矩阵,其上标N<sub>B</sub>、N<sub>S</sub>分别表示行数和列数;S3、根据步骤S1和S2的结果计算VACE的常数矩阵;S4、步骤S1、S2和S3的结果计算VACE;S5、计算VACE对区域S内发电机电压的灵敏度;S6、将VACE控制约束条件加入本区域无功电压优化模型,以进行本区域的目标优化;所述步骤S1中的灵敏度方程为:<maths num="0001"><![CDATA[<math><mrow><mo>-</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>B</mi><mi>DD</mi></msub></mtd><mtd><msub><mi>B</mi><mi>DG</mi></msub></mtd><mtd><msub><mi>B</mi><mi>DB</mi></msub></mtd><mtd></mtd></mtr><mtr><mtd><msub><mi>B</mi><mi>GD</mi></msub></mtd><mtd><msub><mi>B</mi><mi>GG</mi></msub></mtd><mtd><msub><mi>B</mi><mi>GB</mi></msub></mtd><mtd></mtd></mtr><mtr><mtd><msub><mi>B</mi><mi>BD</mi></msub></mtd><mtd><msub><mi>B</mi><mi>BG</mi></msub></mtd><mtd><msub><mi>B</mi><mi>BB</mi></msub></mtd><mtd><msub><mi>B</mi><mi>BE</mi></msub></mtd></mtr><mtr><mtd></mtd><mtd></mtd><mtd><msub><mi>B</mi><mi>EB</mi></msub></mtd><mtd><msub><mi>B</mi><mi>EE</mi></msub></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Delta;</mi><msub><mi>V</mi><mi>D</mi></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>V</mi><mi>G</mi></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>V</mi><mi>B</mi></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>V</mi><mi>E</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mi>&Delta;</mi><msub><mi>Q</mi><mi>D</mi></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>Q</mi><mi>G</mi></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>Q</mi><mi>B</mi></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>Q</mi><mi>E</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,B代表节点的电纳矩阵,ΔV代表节点的电压幅值变化量组成的向量,ΔQ代表节点的外区域注入本区域无功变化量组成的向量;下标B代表区域S侧的边界节点,下标E代表区域S以外网络的节点;下标D代表区域S内除了边界节点以外的所有PQ节点;下标G代表代表区域S内除了边界节点以外的所有PV节点和松弛节点,区域S的节点数目为N<sub>S</sub>,边界节点B的数目为N<sub>B</sub>;B<sub>DD</sub>表示D类节点到D类节点之间的电纳;B<sub>DG</sub>表示D类节点到G类节点之间的电纳;B<sub>DB</sub>、B<sub>GD</sub>、B<sub>GG</sub>、B<sub>GB</sub>、B<sub>BD</sub>、B<sub>BG</sub>、B<sub>BB</sub>、B<sub>BE</sub>、B<sub>EB</sub>、B<sub>EE</sub>的含义依次类推;所述步骤S2中的3个矩阵的元素取值由下式(2)至(4)给出:<img file="FDA00001824052800021.GIF" wi="1591" he="149" /><img file="FDA00001824052800022.GIF" wi="1591" he="149" /><maths num="0002"><![CDATA[<math><mrow><msubsup><mover><mi>B</mi><mo>~</mo></mover><mi>ij</mi><mi>Z</mi></msubsup><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>b</mi><mi>z</mi></msub></mtd><mtd><mi>i</mi><mo>=</mo><mi>j</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>i</mi><mo>&NotEqual;</mo><mi>j</mi></mtd></mtr></mtable></mfenced><mo>&ForAll;</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>&Element;</mo><mi>B</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,b<sub>l</sub>是节点b所连支路的等效串联电纳;b<sub>c</sub>是节点b所连支路的等效并联电纳;b<sub>z</sub>是边界节点所连接的用于进行无功补偿的并联无功补偿装置的电纳;所述步骤S3中的VACE的常数矩阵的计算公式如下:<maths num="0003"><![CDATA[<math><mrow><msub><mi>K</mi><mi>QS</mi></msub><mo>=</mo><msub><mi>B</mi><mi>BD</mi></msub><msup><msub><mi>B</mi><mi>DD</mi></msub><mrow><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>B</mi><mi>DB</mi></msub><mo>+</mo><msubsup><mover><mi>B</mi><mo>~</mo></mover><mi>B</mi><mi>L</mi></msubsup><mo>-</mo><mi>diag</mi><mrow><mo>(</mo><msup><mover><mi>B</mi><mo>~</mo></mover><mi>L</mi></msup><mi>e</mi><mo>)</mo></mrow><mo>-</mo><mn>2</mn><mi>diag</mi><mrow><mo>(</mo><msup><mover><mi>B</mi><mo>~</mo></mover><mi>C</mi></msup><mi>e</mi><mo>)</mo></mrow><mo>-</mo><mn>2</mn><msup><mover><mi>B</mi><mo>~</mo></mover><mi>Z</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,<img file="FDA00001824052800025.GIF" wi="398" he="68" />的所有元素都是1,<img file="FDA00001824052800026.GIF" wi="45" he="49" />是<img file="FDA00001824052800027.GIF" wi="44" he="39" />中边界节点所对应的列组成的矩阵,diag表示求对角矩阵;所述步骤S4中的VACE的计算公式如下:<maths num="0004"><![CDATA[<math><mrow><msub><mi>I</mi><mi>VACE</mi></msub><mo>=</mo><mrow><mo>(</mo><msub><mover><mi>Q</mi><mo>~</mo></mover><mi>B</mi></msub><mo>-</mo><msub><mover><mi>Q</mi><mo>~</mo></mover><mi>Bref</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>E</mi><mi>S</mi></msub><mrow><mo>(</mo><msub><mi>V</mi><mi>B</mi></msub><mo>-</mo><msub><mi>V</mi><mi>Bref</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,E<sub>S</sub>=-K<sub>QS</sub>,V<sub>Bref</sub>和<img file="FDA00001824052800029.GIF" wi="81" he="57" />分别是边界节点电压和外区域注入本区域无功的设定值组成的向量,V<sub>B</sub>和<img file="FDA000018240528000210.GIF" wi="49" he="52" />分别是边界节点电压和外区域注入本区域无功的当前值组成的向量;所述步骤S5中VACE对区域S内发电机电压的灵敏度的计算公式如下:<maths num="0005"><![CDATA[<math><mrow><msub><mi>C</mi><mi>GVACE</mi></msub><mo>=</mo><mfrac><mrow><mi>&Delta;</mi><msup><msub><mi>I</mi><mi>VACE</mi></msub><mrow><mo>(</mo><mi>S</mi><mo>)</mo></mrow></msup></mrow><mrow><mi>&Delta;</mi><msub><mi>V</mi><mi>G</mi></msub></mrow></mfrac><mo>=</mo><msup><mover><mi>B</mi><mo>~</mo></mover><mi>L</mi></msup><msub><mi>C</mi><mi>GS</mi></msub><mo>-</mo><msup><mover><mi>B</mi><mo>~</mo></mover><mi>L</mi></msup><msub><mi>C</mi><mi>BS</mi></msub><msub><mi>C</mi><mi>GB</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,C<sub>GS</sub>是区域S内所有节点电压对发电机节点电压的灵敏度,C<sub>GB</sub>是区域S内所有边界节点电压对发电机节点电压的灵敏度,C<sub>BS</sub>是区域S内所有节点电压对边界节点电压的灵敏度;ΔI<sub>VACE</sub><sup>(S)</sup>是区域S内由发电机节点电压V<sub>G</sub>的变化量ΔV<sub>G</sub>所引起的VACE的变化量;步骤S6中,将如下VACE控制约束条件加入本区域无功电压优化模型,以进行本区域的目标优化:-δ≤I<sub>VACE</sub>+C<sub>GVACE</sub>ΔV<sub>G</sub>≤δ                    (8)其中,δ为VACE允许的变化范围;所述优化模型中可调节的变量为无功功率或电压,故称为无功电压优化模型;步骤S5中三个灵敏度矩阵的计算公式如下:<maths num="0006"><![CDATA[<math><mrow><msub><mi>C</mi><mi>GB</mi></msub><mo>=</mo><mo>-</mo><msup><mrow><mo>(</mo><msub><mi>B</mi><mi>BB</mi></msub><mo>-</mo><msub><mi>B</mi><mi>BE</mi></msub><msubsup><mi>B</mi><mi>EE</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>B</mi><mi>EB</mi></msub><mo>-</mo><msub><mi>B</mi><mi>BD</mi></msub><msubsup><mi>B</mi><mi>DD</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>B</mi><mi>DB</mi></msub><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msub><mi>B</mi><mi>BG</mi></msub><mo>-</mo><msub><mi>B</mi><mi>BD</mi></msub><msubsup><mi>B</mi><mi>DD</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>B</mi><mi>DG</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><msub><mi>C</mi><mi>GS</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mo>-</mo><msubsup><mi>B</mi><mi>DD</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>B</mi><mi>DG</mi></msub><mo>+</mo><msub><mi>B</mi><mi>DB</mi></msub><msub><mi>C</mi><mi>GB</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>I</mi></mtd></mtr><mtr><mtd><msub><mi>C</mi><mi>GB</mi></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msub><mi>C</mi><mi>BS</mi></msub><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mo>-</mo><msubsup><mi>B</mi><mi>DD</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>B</mi><mi>DB</mi></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mi>I</mi></mtd></mtr></mtable></mfenced><mo>.</mo></mrow></math>]]></maths>
地址 110181 辽宁省沈阳市浑南新区营盘北街1号