发明名称 静态混成自动电压控制最高层经济性调控方法
摘要 本发明属于电网电压无功自动控制技术领域,其特征在于在根据给定的电力系统网络参数和各母线节点的电压、电流、有功和无功功率实测值形成Jacobian矩阵,求出实际的总有功损耗值;再以总有功损耗最小为目标构建由状态变量和控制变量组成的优化目标函数,得到系统总有功损耗最小值;然后以系统的实际的总有功损耗值和优化的最小总有功损耗值形成离散的逻辑事件,在实际总有功损耗值大于优化的最小总有功损耗值W倍时,把此时的求解优化目标函数得到的各发电机节点的电压值、无功补偿器节点的电压值和有载调压变压器的变比值作为控制指令输出,以调整上述节点的参数。本发明在确保电压质量的同时提高运行的经济性。
申请公布号 CN100483888C 申请公布日期 2009.04.29
申请号 CN200710098782.7 申请日期 2007.04.27
申请人 清华大学 发明人 卢强;陶家琪;胡伟;张雪敏;陈颖;梅生伟;高德宾
分类号 H02J3/00(2006.01)I;H02J3/12(2006.01)I;H02J3/16(2006.01)I;H02J3/38(2006.01)I;H02J3/46(2006.01)I 主分类号 H02J3/00(2006.01)I
代理机构 代理人
主权项 1、静态混成自动电压控制最高层经济性调控方法,其特征在于,该方法在静态混成自动电压控制最高层经济性调控计算机中依次按以下步骤实现:步骤(1):初始化设定:①电力系统的网络参数,其中包括输电线路的串联电阻、串联电抗、并联电导和并联电纳;变压器的变比和阻抗;并联在输电线路上的电容器和电抗器的阻抗;电力系统的母线节点的名称;②给定的发电机节点个数α<sub>G</sub>和母线名,无功补偿器节点个数α<sub>S</sub>和母线名,除了发电机节点和无功补偿器节点以外的母线节点的个数α<sub>D</sub>和母线名,以及有载调压变压器个数α<sub>F</sub>和母线名;③电力系统的静态电压经济性指标W,5≥W≥1,给定:电力系统的实时量测数据,其中包括各母线节点的电压、电流、有功功率和无功功率,以及有载调压变压器的变比值;步骤(2):在第k个采样间隔后的当前时刻,根据步骤(1)中所述的实时数据,以及电力系统的网络参数,形成电力系统当前潮流计算使用的Jacobian矩阵,并按照下式计算当前时刻电力系统的实时网络有功功率损耗值S<sub>loss</sub>[k],<maths num="0001"><![CDATA[<math><mrow><msub><mi>S</mi><mi>loss</mi></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>=</mo><mi>Re</mi><mrow><mo>{</mo><mover><mi>V</mi><mo>&CenterDot;</mo></mover><msup><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mi>T</mi></msup><mo>&CenterDot;</mo><mrow><mo>(</mo><mi>Y</mi><mo>&CenterDot;</mo><mover><mi>V</mi><mo>&CenterDot;</mo></mover><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>)</mo></mrow><mo>}</mo></mrow></mrow></math>]]></maths>其中Y为以地为参考节点建立的节点导纳矩阵,符号“T”为矩阵的转置操作符号,<maths num="0002"><![CDATA[<math><mfenced open='' close='' separators=','><mtable><mtr><mtd><mover><mi>V</mi><mo>&CenterDot;</mo></mover><msub><mrow><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>=</mo><mi>V</mi></mrow><mrow><mi>G</mi><mn>1</mn></mrow></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&theta;</mi><mrow><mi>G</mi><mn>1</mn></mrow></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow></msup></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>V</mi><msub><mi>G&alpha;</mi><mi>G</mi></msub></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&theta;</mi><msub><mi>G&alpha;</mi><mi>G</mi></msub></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow></msup></mtd><mtd><msub><mi>V</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&theta;</mi><mrow><mi>S</mi><mn>1</mn></mrow></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow></msup></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>V</mi><msub><mi>S&alpha;</mi><mi>S</mi></msub></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&theta;</mi><mi>S&alpha;S</mi></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow></msup></mtd></mtr></mtable></mfenced></math>]]></maths><maths num="0003"><![CDATA[<math><msup><mrow><mfenced open='' close=']' separators=''><mtable><mtr><mtd><msub><mi>V</mi><mrow><mi>D</mi><mn>1</mn></mrow></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&theta;</mi><mrow><mi>D</mi><mn>1</mn></mrow></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow></msup></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>V</mi><msub><mi>D&alpha;</mi><mi>D</mi></msub></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><msup><mi>e</mi><mrow><mi>j</mi><msub><mi>&theta;</mi><mi>D&alpha;D</mi></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mrow></msup></mtd></mtr></mtable></mfenced></mrow><mi>T</mi></msup></math>]]></maths>其中<img file="C200710098782C00024.GIF" wi="265" he="58" />为第m个发电机节点电压量测值,幅值为V<sub>Gm</sub>[k],相角为θ<sub>Gm</sub>[k],m为发电机节点的序号,m=1,…,α<sub>G</sub>,<img file="C200710098782C00025.GIF" wi="244" he="58" />为第n个无功补偿器节点的电压量测值,幅值为V<sub>Sn</sub>[k],相角为θ<sub>Sn</sub>[k],n为无功补偿器节点的序号,n=1,…,α<sub>S</sub>,<img file="C200710098782C00026.GIF" wi="252" he="57" />为第r个除了发电机节点和无功补偿器节点以外的母线节点的电压量测值,幅值为V<sub>Dr</sub>[k],相角为θ<sub>Dr</sub>[k],r为除了发电机节点和无功补偿器节点以外的母线节点的序号,r=1,…,α<sub>D</sub>,<img file="C200710098782C00027.GIF" wi="96" he="49" />为系统所有母线节点电压量测值组成的列向量;步骤(3):根据步骤(1)得到的当前时刻电力系统的实时数据,根据下式计算最优化潮流条件下电力系统总的线路有功功率损耗的最小值S<sub>eco</sub>[k],S<sub>eco</sub>[k]=J(x[k],u[k])<maths num="0004"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><munder><mi>min</mi><mi>u</mi></munder><mi>J</mi><mrow><mo>(</mo><mi>x</mi><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>,</mo><mi>u</mi><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mi></mi><mi></mi><mi></mi><mi></mi><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>,</mo><mi>u</mi><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mi>h</mi><mrow><mo>(</mo><mi>x</mi><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>,</mo><mi>u</mi><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>)</mo></mrow><mo>&le;</mo><mn>0</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中<maths num="0005"><![CDATA[<math><mrow><mi>u</mi><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>=</mo><msup><mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>V</mi><mrow><mi>G</mi><mn>1</mn></mrow></msub></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>V</mi><msub><mi>G&alpha;</mi><mi>G</mi></msub></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mtd><mtd><msub><mi>V</mi><mrow><mi>S</mi><mn>1</mn></mrow></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>V</mi><msub><mi>S&alpha;</mi><mi>S</mi></msub></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mtd><mtd><msub><mi>t</mi><mn>1</mn></msub><mrow><mo>]</mo><mi>k</mi><mo>]</mo></mrow></mtd><mtd><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo></mtd><mtd><msub><mi>t</mi><msub><mi>&alpha;</mi><mi>F</mi></msub></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow></mtd></mtr></mtable></mfenced></mrow><mi>T</mi></msup></mrow></math>]]></maths>u[k]代表控制变量,包括各发电机节点的电压参考值、各无功补偿器节点的电压参考值和有载调压变压器的变比参考值,x[k]代表状态变量,包括系统中除了控制变量以外各母线节点的电压、电流、有功功率和无功功率值,J(x[k],u[k])为描述系统运行经济性的目标函数,即总的线路损耗最小的优化目标函数f(x[k],u[k])=0为系统潮流方程,h(x[k],u[k])≤0为系统运行时的不等式约束,步骤(4):根据步骤(3)所述的基于总的线路损耗最小的优化目标函数求出:各发电机节点的电压参考值<img file="C200710098782C00033.GIF" wi="166" he="58" />m为发电机节点的序号,m=1,…,α<sub>G</sub>,各无功补偿器节点的电压参考值<img file="C200710098782C00034.GIF" wi="167" he="59" />n为无功补偿器节点的序号,n=1,…,α<sub>S</sub>,和有载调压变压器的变比参考值<img file="C200710098782C00035.GIF" wi="161" he="53" />其中l为有载调压变压器F的序号,l=1,…,α<sub>F</sub>;步骤(5):设定如下的逻辑条件来构成“事件”E<sub>eco</sub>[k]来形成离散事件<maths num="0006"><![CDATA[<math><mrow><msub><mi>E</mi><mi>eco</mi></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>=</mo><mrow><mfenced open='{' close='' separators=','><mtable><mtr><mtd><msub><mi>P</mi><mi>econ</mi></msub><mo>,</mo></mtd><mtd><mi>if</mi></mtd><mtd><msub><mi>S</mi><mi>loss</mi></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>></mo><mi>W</mi><mo>&times;</mo><msub><mi>S</mi><mi>eco</mi></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>,</mo></mtd></mtr><mtr><mtd><msub><mi>P</mi><mi>enon</mi></msub><mo>,</mo></mtd><mtd><mi>if</mi></mtd><mtd><msub><mi>S</mi><mi>loss</mi></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>&le;</mo><mi>W</mi><mo>&times;</mo><msub><mi>S</mi><mi>eco</mi></msub><mrow><mo>[</mo><mi>k</mi><mo>]</mo></mrow><mo>,</mo></mtd></mtr></mtable></mfenced></mrow></mrow></math>]]></maths>其中P<sub>econ</sub>代表当前需要进行电压经济控制,P<sub>enon</sub>代表当前不需要进行电压经济控制,当事件为P<sub>enon</sub>时,则经济性调控环节输出的控制指令为0,当事件为P<sub>econ</sub>时,则经济性调控环节把步骤(4)得到的各控制节点的控制参考量作为控制指令输出。
地址 100084北京市100084信箱82分箱清华大学专利办公室