发明名称 一种晶闸管控制变压器式可控电抗器优化控制方法
摘要 本发明公开了一种晶闸管控制变压器式可控电抗器优化控制方法,包括下述步骤:1)检测负荷三相功率;2)计算电流谐波畸变率;3)判断谐波是否越限;4)构造电能质量指标的评分函数;5)利用层次分析法确定指标的权重;6)以电能质量最优为目的,确定三相触发角。本发明方法能够在不增加设备的情况下,依靠可控电抗器晶闸管触发角的优化选择,降低电抗器产生的谐波,改善线路的电能质量。
申请公布号 CN103683304A 申请公布日期 2014.03.26
申请号 CN201310659680.3 申请日期 2013.12.09
申请人 国家电网公司;江苏省电力公司;江苏省电力公司电力科学研究院;东南大学 发明人 袁晓冬;王元凯;顾伟;李群;张帅;柏晶晶
分类号 H02J3/18(2006.01)I;H02J3/01(2006.01)I 主分类号 H02J3/18(2006.01)I
代理机构 南京纵横知识产权代理有限公司 32224 代理人 董建林
主权项 1.一种晶闸管控制变压器式可控电抗器优化控制方法,其特征在于,该控制方法包括下述步骤:1)检测可控电抗器未接入时负荷三相有功功率P<sub>L(ab)</sub>、P<sub>L(bc)</sub>、P<sub>L(ca)</sub>、三相无功功率Q<sub>L(ab)</sub>、Q<sub>L(bc)</sub>、Q<sub>L(ca)</sub>;2)计算三相电流谐波畸变率,具体的过程为:21)利用步骤1)检测的负荷三相有功功率和三相无功功率P<sub>L(ab)</sub>、P<sub>L(bc)</sub>、P<sub>L(ca)</sub>、Q<sub>L(ab)</sub>、Q<sub>L(bc)</sub>、Q<sub>L(ca)</sub>,求解以下方程组,得到线路功率因数<img file="FDA0000432841400000014.GIF" wi="153" he="60" />且三相电流不平衡度ε=0时的三相触发角α<sub>a0</sub>、α<sub>b0</sub>、α<sub>c0</sub>:<![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><msub><mi>P</mi><mo>+</mo></msub><msqrt><msup><mrow><mo>(</mo><msub><mi>P</mi><mo>+</mo></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Q</mi><mo>+</mo></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac><mo>=</mo><mn>1</mn></mtd></mtr><mtr><mtd><mfrac><msqrt><msup><mrow><mo>(</mo><msub><mi>P</mi><mo>-</mo></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Q</mi><mo>-</mo></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><msqrt><msup><mrow><mo>(</mo><msub><mi>P</mi><mo>+</mo></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>Q</mi><mo>+</mo></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced></math>]]></maths>其中P<sub>+</sub>=P<sub>L(ab)</sub>+P<sub>L(bc)</sub>+P<sub>L(ca)</sub>为正序有功功率,Q<sub>+</sub>=(Q<sub>L(ab)</sub>+Q<sub>c(ab)</sub>)+(Q<sub>L(bc)</sub>+Q<sub>c(bc)</sub>)+(Q<sub>L(ca)</sub>+Q<sub>c(ca)</sub>)为正序无功功率,<![CDATA[<math><mrow><msub><mi>P</mi><mo>-</mo></msub><mo>=</mo><mn>0.5</mn><mo>&times;</mo><mo>&lsqb;</mo><msub><mi>P</mi><mrow><mi>L</mi><mrow><mo>(</mo><mi>ab</mi><mo>)</mo></mrow></mrow></msub><mo>-</mo><mn>2</mn><msub><mi>P</mi><mrow><mi>L</mi><mrow><mo>(</mo><mi>bc</mi><mo>)</mo></mrow></mrow></msub><mo>+</mo><msub><mi>P</mi><mrow><mi>L</mi><mrow><mo>(</mo><mi>ca</mi><mo>)</mo></mrow></mrow></msub><mo>-</mo><msqrt><mn>3</mn></msqrt><mrow><mo>(</mo><msub><mi>Q</mi><mrow><mi>L</mi><mrow><mo>(</mo><mi>ca</mi><mo>)</mo></mrow></mrow></msub><mo>+</mo><msub><mi>Q</mi><mrow><mi>c</mi><mrow><mo>(</mo><mi>ca</mi><mo>)</mo></mrow></mrow></msub><mo>-</mo><msub><mi>Q</mi><mrow><mi>L</mi><mrow><mo>(</mo><mi>ab</mi><mo>)</mo></mrow></mrow></msub><mo>-</mo><msub><mi>Q</mi><mrow><mi>c</mi><mrow><mo>(</mo><mi>ab</mi><mo>)</mo></mrow></mrow></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></math>]]></maths>为负序有功功率,<![CDATA[<math><mrow><msub><mi>Q</mi><mo>-</mo></msub><mo>=</mo><mn>0.5</mn><mo>&times;</mo><mo>&lsqb;</mo><msub><mi>Q</mi><mrow><mi>L</mi><mrow><mo>(</mo><mi>ab</mi><mo>)</mo></mrow></mrow></msub><mo>+</mo><msub><mi>Q</mi><mrow><mi>c</mi><mrow><mo>(</mo><mi>ab</mi><mo>)</mo></mrow></mrow></msub><mo>-</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>Q</mi><mrow><mi>L</mi><mrow><mo>(</mo><mi>bc</mi><mo>)</mo></mrow></mrow></msub><mo>+</mo><msub><mi>Q</mi><mrow><mi>c</mi><mrow><mo>(</mo><mi>bc</mi><mo>)</mo></mrow></mrow></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>Q</mi><mrow><mi>L</mi><mrow><mo>(</mo><mi>ca</mi><mo>)</mo></mrow></mrow></msub><mo>+</mo><msub><mi>Q</mi><mrow><mi>c</mi><mrow><mo>(</mo><mi>ca</mi><mo>)</mo></mrow></mrow></msub><mo>+</mo><msqrt><mn>3</mn></msqrt><mrow><mo>(</mo><msub><mi>P</mi><mrow><mi>L</mi><mrow><mo>(</mo><mi>ca</mi><mo>)</mo></mrow></mrow></msub><mo>-</mo><msub><mi>P</mi><mrow><mi>L</mi><mrow><mo>(</mo><mi>ab</mi><mo>)</mo></mrow></mrow></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow></math>]]></maths>为负序无功功率;Q<sub>c(ab)</sub>=U<sup>2</sup>×(sin2α<sub>a0</sub>-2α<sub>a0</sub>+2π)/πX,Q<sub>c(bc)</sub>=U<sup>2</sup>×(sin2α<sub>b0</sub>-2α<sub>b0</sub>+2π)/πX,Q<sub>c(ca)</sub>=U<sup>2</sup>×(sin2α<sub>c0</sub>-2α<sub>c0</sub>+2π)/πX分别为a、b、c三相补偿的无功功率,其中U为网点线电压有效值,X为变压器漏抗;22)按照下式计算当触发角为α<sub>a0</sub>、α<sub>b0</sub>、α<sub>c0</sub>时的谐波畸变率THD<sub>0</sub>:<![CDATA[<math><mrow><msub><mi>THD</mi><mn>0</mn></msub><mo>=</mo><mrow><mo>(</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>h</mi><mo>=</mo><mn>3</mn></mrow><mn>11</mn></munderover><msubsup><mi>I</mi><mi>ha</mi><mn>2</mn></msubsup></msqrt><mo>/</mo><msub><mi>I</mi><mi>fa</mi></msub><mo>+</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>h</mi><mo>=</mo><mn>3</mn></mrow><mn>11</mn></munderover><msubsup><mi>I</mi><mi>hb</mi><mn>2</mn></msubsup></msqrt><mo>/</mo><msub><mi>I</mi><mi>fb</mi></msub><mo>+</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>h</mi><mo>=</mo><mn>3</mn></mrow><mn>11</mn></munderover><msubsup><mi>I</mi><mi>hc</mi><mn>2</mn></msubsup></msqrt><mo>/</mo><msub><mi>I</mi><mi>fc</mi></msub><mo>)</mo></mrow><mo>/</mo><mn>3</mn></mrow></math>]]></maths>其中I<sub>fa</sub>、I<sub>fb</sub>、I<sub>fc</sub>分别为流入可控电抗器a、b、c三相的基波电流有效值,I<sub>ha</sub>、I<sub>hb</sub>、I<sub>hc</sub>分别为流入可控电抗器a、b、c三相的谐波电流有效值,下标h代表谐波次数;基波电流有效值I<sub>f</sub>具体表达式为:<![CDATA[<math><mrow><msub><mi>I</mi><mi>f</mi></msub><mo>=</mo><mfrac><mi>U</mi><mi>&pi;X</mi></mfrac><msup><mrow><mo>&lsqb;</mo><msubsup><mi>A</mi><mi>f</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>B</mi><mi>f</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>A</mi><mi>f</mi></msub><msub><mi>B</mi><mi>f</mi></msub><mo>&rsqb;</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup></mrow></math>]]></maths>A<sub>f</sub>=sin(2α)-2α+2πB<sub>f</sub>=sin(2β)-2β+2π其中α、β为可控电抗器a、b、c三相中其中两相的触发角,A<sub>f</sub>为α触发角对应相的基波电流系数,B<sub>f</sub>为β触发角对应相的基波电流系数,对流入可控电抗器a相的基波电流I<sub>fa</sub>:α=α<sub>a0</sub>对应a相,β=α<sub>c0</sub>对应c相;对流入可控电抗器中b相的基波电流I<sub>f</sub>b:α=α<sub>b0</sub>对应b相,β=α<sub>a0</sub>对应a相;对流入可控电抗器中c相的基波电流I<sub>fc</sub>:α=α<sub>c0</sub>对应c相,β=α<sub>b0</sub>对应b相;谐波电流有效值I<sub>h</sub>具体表达式为:<img file="FDA0000432841400000021.GIF" wi="1040" he="447" /><![CDATA[<math><mrow><msub><mi>A</mi><mi>h</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mi>h</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mi>sin</mi><mrow><mo>(</mo><mi>h</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mi>&alpha;</mi><mo>+</mo><mfrac><mn>1</mn><mrow><mi>h</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mi>sin</mi><mrow><mo>(</mo><mi>h</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>&alpha;</mi><mo>-</mo><mfrac><mn>2</mn><mi>h</mi></mfrac><mi>cos</mi><mi></mi><mi>&alpha;</mi><mi>sin</mi><mrow><mo>(</mo><mi>h&alpha;</mi><mo>)</mo></mrow></mrow></math>]]></maths><![CDATA[<math><mrow><msub><mi>B</mi><mi>h</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mi>h</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mi>sin</mi><mrow><mo>(</mo><mi>h</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mi>&beta;</mi><mo>+</mo><mfrac><mn>1</mn><mrow><mi>h</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mi>sin</mi><mrow><mo>(</mo><mi>h</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>&beta;</mi><mo>-</mo><mfrac><mn>2</mn><mi>h</mi></mfrac><mi>cos</mi><mi></mi><mi>&beta;</mi><mi>sin</mi><mrow><mo>(</mo><mi>h&beta;</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中k为正整数,h为谐波次数,A<sub>h</sub>为α触发角对应相的谐波电流系数,B<sub>h</sub>为β触发角对应相的谐波电流系数,对流入可控电抗器a相的谐波电流I<sub>ha</sub>:α=α<sub>a0</sub>对应a相,β=α<sub>c0</sub>对应c相;对流入可控电抗器b相的谐波电流I<sub>hb</sub>:α=α<sub>b0</sub>对应b相,β=α<sub>a0</sub>对应a相;对流入可控电抗器c相的电流谐波I<sub>hc</sub>:α=α<sub>c0</sub>对应c相,β=α<sub>b0</sub>对应b相;3)用户根据对应的国标或实际情况,分别设定线路功率因数<img file="FDA0000432841400000024.GIF" wi="97" he="48" />、三相电流不平衡度ε、三相电流谐波畸变率THD的阈值,若步骤22)得到的谐波畸变率THD<sub>0</sub>小于电流谐波畸变率THD的阈值L<sub>THD</sub>,则将触发角α<sub>a0</sub>、α<sub>b0</sub>、α<sub>c0</sub>作为电抗器晶闸管触发角输出后结束优化控制方法流程,否则进入步骤4);4)构造电能质量指标的评分函数S<sub>THD</sub>、S<sub>ε</sub>、<img file="FDA0000432841400000025.GIF" wi="94" he="58" />;评分函数S<sub>THD</sub>、S<sub>ε</sub>、<img file="FDA0000432841400000026.GIF" wi="101" he="63" />分别为三相电流谐波畸变率、三相电流不平衡度、线路功率因数三项电能质量指标值到指标评分值的函数关系,其评分值从0到100,用来分别表征三项电能质量指标的优劣情况;5)利用层次分析法分别确定三项电能质量指标的权重ω<sub>THD</sub>、ω<sub>ε</sub>、<img file="FDA0000432841400000027.GIF" wi="105" he="49" />;三项电能质量指标的权重ω<sub>THD</sub>、ω<sub>ε</sub>、<img file="FDA0000432841400000028.GIF" wi="109" he="54" />表示三项指标在电能质量整体评价中的相对重要程度;6)以电能质量最优为目的,利用所述步骤4)得到的电能质量指标的评分函数,以及步骤5)得到的三项电能质量指标的权重,来构造目标函数<img file="FDA00004328414000000210.GIF" wi="372" he="61" /><img file="FDA0000432841400000029.GIF" wi="404" he="66" />然后以步骤21)计算的三相触发角α<sub>a0</sub>、α<sub>b0</sub>、α<sub>c0</sub>作为初值,利用牛顿法确定各相触发角α<sub>a</sub>、α<sub>b</sub>、α<sub>c</sub>并输出。
地址 100031 北京市西城区西长安街86号