发明名称 基于非线性鲁棒电力系统稳定器的励磁控制方法
摘要 基于非线性鲁棒电力系统稳定器的励磁控制方法,属于电力系统稳定性控制领域,其特征在于,在多机励磁系统非线性数学模型基础上,用微分几何控制理论与非线性H<SUB>∞</SUB>方法相结合,提出了一种非线性鲁棒电力系统稳定器,其输出信号与原有的自动调压器输出控制信号相加,所得的励磁控制信号,对发电机组的激励系统进行控制,从而避免了外界不确定性因素的影响,解决了电力系统动态呈强非线性和高耦合性的问题,克服了现有技术中忽略了系统固有的非线性特征的缺陷,改善了动态稳定性,最大限度地发电机组的装机容量。
申请公布号 CN101127445A 申请公布日期 2008.02.20
申请号 CN200710119266.8 申请日期 2007.07.19
申请人 清华大学 发明人 卢强;梅生伟;薛建伟;郑少明;彭晓洁
分类号 H02J3/24(2006.01);H02P9/14(2006.01) 主分类号 H02J3/24(2006.01)
代理机构 代理人
主权项 1.基于非线性鲁棒电力系统稳定器励磁控制方法,其特征在于,所述方法是在数字集成电路中依次按照以下步骤完成的:步骤(1)初始化以下发电机组参数并输入所述数字集成电路中:稳态频率值ω<sub>0i</sub>=314.16,历次绕组时间常数I<sub>d0i</sub>′,机组转动惯量H<sub>i</sub>,D轴电抗X<sub>di</sub>、Q轴电抗X<sub>qi</sub>、D轴暂态电抗X<sub>di</sub>′、定子电阻r<sub>i</sub>,其中:i为所述发电机组中第i台发电机的称号;步骤(2)利用电压互感器PT、电流互感器CT实时测得以下个各量:发电机i的三相电压u<sub>ai</sub>,u<sub>bi</sub>,u<sub>ci</sub>,三相电流瞬时值i<sub>ai</sub>,i<sub>bi</sub>,i<sub>ci</sub>,电机转速近似值ω<sub>i</sub>=2πf<sub>i</sub>,系统频率f<sub>i</sub>可以直接测得,再把u<sub>ai</sub>,u<sub>bi</sub>,u<sub>ci</sub>,i<sub>ai</sub>,i<sub>bi</sub>,i<sub>ci</sub>输入到该数字集成电路;步骤(3)根据步骤(2)输入的各个瞬时值计算得以下状态量:线电压有效值V<sub>ti</sub>,电流有效值I<sub>ti</sub>,有功功率P<sub>ei</sub>,无功功率Q<sub>ei</sub>;步骤(4)分别以发电机的额定电压V<sub>B</sub>、额定电流I<sub>B</sub>、额定阻抗Z<sub>B</sub>=V<sub>B</sub>/I<sub>B</sub>为基值,再以发电机组视在功率S<sub>N</sub>以及定子角频率ω<sub>B</sub>=ω<sub>0</sub>=314.1529为基值队步骤(3)得到的各状态量作标幺化处理;步骤(5)在步骤(4)得到的各标幺值且r<sub>1</sub>=0的条件下计算下列各状态量:●Q轴电势:<maths num="0001"><![CDATA[<math><mrow><msub><mi>E</mi><mi>qi</mi></msub><mo>=</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>V</mi><mi>ti</mi></msub><mo>+</mo><mfrac><mrow><msub><mi>Q</mi><mi>ei</mi></msub><msub><mi>x</mi><mi>di</mi></msub></mrow><msub><mi>V</mi><mi>ti</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>P</mi><mi>ei</mi></msub><msub><mi>x</mi><mi>di</mi></msub></mrow><msub><mi>V</mi><mi>ti</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>,</mo></mrow></math>]]></maths>●Q轴暂态电势:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>E</mi><mi>qi</mi><mo>&prime;</mo></msubsup><mo>=</mo><msqrt><msup><mrow><mo>(</mo><msub><mi>V</mi><mi>ti</mi></msub><mo>+</mo><mfrac><mrow><msub><mi>Q</mi><mi>ei</mi></msub><msubsup><mi>x</mi><mi>di</mi><mo>&prime;</mo></msubsup></mrow><msub><mi>V</mi><mi>ti</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>P</mi><mi>ei</mi></msub><msubsup><mi>x</mi><mi>di</mi><mo>&prime;</mo></msubsup></mrow><msub><mi>V</mi><mi>ti</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>,</mo></mrow></math>]]></maths>●D轴电流:<maths num="0003"><![CDATA[<math><mrow><msub><mi>i</mi><mi>di</mi></msub><mo>=</mo><mfrac><msub><mi>c</mi><mi>i</mi></msub><msub><mi>d</mi><mi>i</mi></msub></mfrac><msubsup><mi>E</mi><mi>qi</mi><mo>&prime;</mo></msubsup><msubsup><mi>I</mi><mi>ti</mi><mn>2</mn></msubsup><mo>,</mo></mrow></math>]]></maths>●<maths num="0004"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>a</mi><mi>i</mi></msub><mo>=</mo><msub><mi>Q</mi><mi>ei</mi></msub><mo>+</mo><msubsup><mi>x</mi><mi>di</mi><mo>&prime;</mo></msubsup><msubsup><mi>I</mi><mi>ti</mi><mn>2</mn></msubsup></mtd></mtr><mtr><mtd><msub><mi>b</mi><mi>i</mi></msub><mo>=</mo><msub><mi>P</mi><mi>ei</mi></msub></mtd></mtr><mtr><mtd><msub><mi>c</mi><mi>i</mi></msub><mo>=</mo><msub><mi>Q</mi><mi>ei</mi></msub><mo>+</mo><msub><mi>x</mi><mi>qi</mi></msub><msubsup><mi>I</mi><mi>ti</mi><mn>2</mn></msubsup></mtd></mtr><mtr><mtd><msub><mi>d</mi><mi>i</mi></msub><mo>=</mo><msub><mi>a</mi><mi>i</mi></msub><msub><mi>c</mi><mi>i</mi></msub><mo>+</mo><msubsup><mi>b</mi><mi>i</mi><mn>2</mn></msubsup></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>●Q轴电流:<maths num="0005"><![CDATA[<math><mrow><msub><mi>i</mi><mi>qi</mi></msub><mo>=</mo><mfrac><msub><mi>b</mi><mi>i</mi></msub><msub><mi>d</mi><mi>i</mi></msub></mfrac><msubsup><mi>E</mi><mi>qi</mi><mo>&prime;</mo></msubsup><msubsup><mi>I</mi><mi>ti</mi><mn>2</mn></msubsup><mo>,</mo></mrow></math>]]></maths>步骤(6)按下式计算标幺制下的各反馈量●频率偏差:Δω<sub>i</sub>=(ω<sub>i</sub>-ω<sub>0i</sub>)●功角偏差:<maths num="0006"><![CDATA[<math><mrow><msub><mi>&Delta;&delta;</mi><mi>i</mi></msub><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><msub><mi>&Delta;&omega;</mi><mi>i</mi></msub><mi>d&tau;</mi><mo>,</mo></mrow></math>]]></maths>t为积分周期,dτ为积分间隔,●功率偏差:ΔP<sub>ei</sub>=(P<sub>ei</sub>-P<sub>ei0</sub>),P<sub>ei</sub>由额定功率下相电压、相电流有效值算出,P<sub>ei0</sub>为视在功率,●电压偏差:ΔV<sub>ti</sub>=(V<sub>ti</sub>-V<sub>ti0</sub>),V<sub>ti0</sub>为额定线电压,V<sub>ti</sub>为线电压;步骤(7)计算D轴电流i<sub>di</sub>和Q轴电流i<sub>qi</sub>的微分量<img file="A2007101192660003C2.GIF" wi="39" he="65" />和<img file="A2007101192660003C3.GIF" wi="55" he="71" />:用线性拟合的方法求得;步骤(8)设置Q轴电流i<sub>qi</sub>的最小值,取值区间为[0.1,0.2],按下式计算i<sub>qi</sub>的轻载修正系数k<sub>qi</sub>,以及空载修正系数C<sub>qi</sub>,当i<sub>qi</sub>=i<sub>qimin</sub>,或i<sub>qi</sub>=-i<sub>qimin</sub>时,k<sub>qi</sub>=0,C<sub>qi</sub>=0.1,当i<sub>qi</sub>>i<sub>qimin</sub>,或i<sub>qi</sub><-i<sub>qimin</sub>时,k<sub>qi</sub>=1,C<sub>qi</sub>=1,步骤(9)按照下式计算非线性鲁棒电力系统稳定器NR-PSS的输出V<sub>fiNR-PSS</sub>:<maths num="0007"><![CDATA[<math><mrow><msub><mi>V</mi><mrow><mi>fiNR</mi><mo>-</mo><mi>PSS</mi></mrow></msub><mo>=</mo><msub><mi>E</mi><mi>qi</mi></msub><mo>-</mo><mfrac><msub><mi>T</mi><mrow><mi>d</mi><mn>0</mn><mi>i</mi></mrow></msub><msub><mi>i</mi><mi>qi</mi></msub></mfrac><mo>[</mo><msub><mi>K</mi><mi>qi</mi></msub><msubsup><mi>E</mi><mi>qi</mi><mo>&prime;</mo></msubsup><msub><mover><mi>i</mi><mo>&CenterDot;</mo></mover><mi>qi</mi></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>x</mi><mi>qi</mi></msub><mo>-</mo><msubsup><mi>x</mi><mi>di</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>i</mi><mi>qi</mi></msub><msub><mover><mi>i</mi><mo>&CenterDot;</mo></mover><mi>di</mi></msub><mo>+</mo><msub><mi>K</mi><mi>qi</mi></msub><msub><mi>i</mi><mi>di</mi></msub><msub><mover><mi>i</mi><mo>&CenterDot;</mo></mover><mi>qi</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><mo>+</mo><msub><mi>C</mi><mi>qi</mi></msub><msub><mi>C</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mfrac><msub><mi>HT</mi><mrow><mi>d</mi><mn>0</mn><mi>i</mi></mrow></msub><mrow><msub><mi>&omega;</mi><mn>0</mn></msub><msub><mi>i</mi><mi>qi</mi></msub></mrow></mfrac><mrow><mo>(</mo><msub><mi>k</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><msub><mi>&Delta;&delta;</mi><mi>i</mi></msub><mo>+</mo><msub><mi>k</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><msub><mi>&Delta;&omega;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>k</mi><mrow><mn>3</mn><mi>i</mi></mrow></msub><mfrac><msub><mi>&omega;</mi><mn>0</mn></msub><msub><mi>H</mi><mi>i</mi></msub></mfrac><msub><mi>&Delta;P</mi><mi>ei</mi></msub><mo>)</mo></mrow></mrow></math>]]></maths>其中:反馈系数k<sub>1i</sub>,k<sub>2i</sub>,k<sub>3i</sub>,各自在[0,300]间取值,阻尼调节系数C<sub>1i</sub>,C<sub>1i</sub>=1/m<sub>i</sub>,m<sub>i</sub>在[0.1,10]间取值;步骤(10)按下式计算PID控制的电压调节器AVR的输出V<sub>fiAVR</sub>:<maths num="0009"><![CDATA[<math><mrow><msub><mi>V</mi><mi>fiAVR</mi></msub><mo>=</mo><msub><mrow><mo>-</mo><mi>k</mi></mrow><mi>pi</mi></msub><msub><mi>&Delta;V</mi><mi>ti</mi></msub><mo>-</mo><msub><mi>k</mi><mi>Ii</mi></msub><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><msub><mi>&Delta;V</mi><mi>ti</mi></msub><mi>d&tau;</mi><mo>-</mo><msub><mi>k</mi><mi>di</mi></msub><mi>&Delta;</mi><msub><mover><mi>V</mi><mo>&CenterDot;</mo></mover><mi>ti</mi></msub></mrow></math>]]></maths>其中:ΔV<sub>ti</sub>为机端电压偏差量,k<sub>pi</sub>为比例系数,在[1,300]间取值,k<sub>Ii</sub>为积分系数,在[0,100]间取值,k<sub>di</sub>为微分系数,在[0,100]间取值;步骤(11)按照下式计算用于控制第i台发电机可控硅控制器的励磁信号V<sub>fi</sub>:V<sub>fi</sub>=C<sub>3i</sub>·V<sub>fiAVR</sub>+C<sub>2i</sub>·V<sub>fiNR-PSS</sub>C<sub>1</sub>,其中:C<sub>2i</sub>为NR-PSS的增益系数,取值范围为[0.1,1],C<sub>3i</sub>为AVR的增益系数,取值范围为[0.5,1]。
地址 100084北京市100084-82信箱