主权项 |
一种基于广域相量测量系统的功率缺额估计方法,其特征是,它包括以下步骤:1)基于惯量中心坐标的功率缺额估计 <mrow> <msub> <mi>ω</mi> <mi>COI</mi> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>H</mi> <mi>i</mi> </msub> <mo>·</mo> <msub> <mi>ω</mi> <mi>i</mi> </msub> </mrow> <mrow> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>H</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>1,2</mn> <mo>.</mo> <mo>.</mo> <mo>.</mo> <mi>n</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>其中ωi为第i台发电机的频率,单位为标幺值,Hi为第i台发电机的转动惯量,ωCOI为系统惯量中心频率,单位为标幺值;系统的惯量中心频率反映了电力系统作为一个整体在遭受有功缺额冲击后由系统的整体阻尼特性决定的减速或加速过程,将发电机转速n维受扰轨迹经过惯量中心的变换求出系统总体的受扰轨迹即系统惯量中心频率,既而通过扰动初始时刻惯量中心频率的变化率来估计系统在计及初始的频率调节效应的缺额量; <mrow> <mfrac> <mi>dΔω</mi> <mi>dt</mi> </mfrac> <msub> <mo>|</mo> <mrow> <mi>t</mi> <mo>=</mo> <mn>0</mn> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>P</mi> <mi>STEP</mi> </msub> <msub> <mrow> <mn>2</mn> <mi>H</mi> </mrow> <mi>eq</mi> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>其中: <mrow> <msub> <mi>H</mi> <mi>eq</mi> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>n</mi> </mrow> </munderover> <msub> <mi>H</mi> <mi>i</mi> </msub> <mo>·</mo> <msub> <mi>S</mi> <mi>i</mi> </msub> </mrow> <mrow> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>n</mi> </mrow> </munderover> <msub> <mi>S</mi> <mi>i</mi> </msub> </mrow> </mfrac> </mrow>其中Hi为第i台发电机的转动惯量,Si为第i台发电机的额定视在功率,Heq为系统等值转动惯量;PSTEP为不计及负荷电压突变影响缺额估计量,单位为标幺值;上式中△ω及PSTEP均为标幺值,为了得到有名值通过如下公式折算得到: <mrow> <mi>w</mi> <mo>=</mo> <mfrac> <msub> <mi>w</mi> <mi>rad</mi> </msub> <msub> <mi>w</mi> <mi>N</mi> </msub> </mfrac> <mo>=</mo> <mfrac> <msub> <mi>f</mi> <mi>Hz</mi> </msub> <msub> <mi>f</mi> <mi>N</mi> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>其中ωrad,ωN分别为以弧度为单位的系统频率和系统额定频率;fHz,fN为系统频率和系统额定频率,单位Hz; <mrow> <msub> <mi>P</mi> <mi>STEP</mi> </msub> <mo>=</mo> <mfrac> <msub> <mi>P</mi> <mi>def</mi> </msub> <mrow> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>S</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>Pdef为不计及负荷电压突变影响缺额估计量,单位为MW,Si为系统内第i台机的额定容量;由(2)、(4)式得: <mrow> <msub> <mi>P</mi> <mi>def</mi> </msub> <mo>=</mo> <msub> <mrow> <mn>2</mn> <mi>H</mi> </mrow> <mi>eq</mi> </msub> <mfrac> <mi>dΔω</mi> <mi>dt</mi> </mfrac> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>S</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>上式中有关扰动初始时刻的频率变化率通过数值算法计算得来而数值算法如下式所示: <mrow> <mfrac> <mrow> <mi>df</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> <mi>dt</mi> </mfrac> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>)</mo> </mrow> <msub> <mi>T</mi> <mi>i</mi> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>2)频率及其变化率跟踪算法利用最速跟踪微分器求取输入信号的跟踪信号及微分信号,其具体的离散形式表示为: <mrow> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>hx</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>+</mo> </mtd> </mtr> <mtr> <mtd> <mi>hfh</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>k</mi> <mo>)</mo> </mrow> <mo>,</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <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> <mi>r</mi> <mo>,</mo> <mi>h</mi> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>u(k)为第k时刻输入信号,r,h分别为跟踪快慢和滤波效果参数;式中fh函数的表达式如右式所示: <mrow> <mi>fh</mi> <mo>=</mo> <mo>-</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>r</mi> <mo>·</mo> <mfrac> <mi>a</mi> <mi>δ</mi> </mfrac> <mo>,</mo> </mtd> <mtd> <mo>|</mo> <mi>a</mi> <mo>|</mo> <mo>≤</mo> <mi>δ</mi> </mtd> </mtr> <mtr> <mtd> <mi>rsign</mi> <mrow> <mo>(</mo> <mi>a</mi> <mo>)</mo> </mrow> <mo>,</mo> </mtd> <mtd> <mo>|</mo> <mi>a</mi> <mo>|</mo> <mo>></mo> <mi>δ</mi> </mtd> </mtr> </mtable> </mfenced> </mrow>式中 <mrow> <mi>a</mi> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>+</mo> <mfrac> <mrow> <msub> <mi>a</mi> <mn>0</mn> </msub> <mo>-</mo> <mi>δ</mi> </mrow> <mn>2</mn> </mfrac> <mi>sign</mi> <mrow> <mo>(</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>,</mo> </mtd> <mtd> <mo>|</mo> <mi>y</mi> <mo>|</mo> <mo>></mo> <msub> <mi>δ</mi> <mn>0</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>+</mo> <mfrac> <mi>y</mi> <mi>h</mi> </mfrac> <mo>,</mo> </mtd> <mtd> <mo>|</mo> <mi>y</mi> <mo>|</mo> <mo>≤</mo> <msub> <mi>δ</mi> <mn>0</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mi>y</mi> <mo>=</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>-</mo> <mi>u</mi> <mo>+</mo> <msub> <mi>hx</mi> <mn>2</mn> </msub> <mo>,</mo> </mrow>δ0=δh,δ=rh, <mrow> <msub> <mi>a</mi> <mn>0</mn> </msub> <mo>=</mo> <msqrt> <msup> <mi>δ</mi> <mn>2</mn> </msup> <mo>+</mo> <mn>8</mn> <mi>r</mi> <mo>|</mo> <mi>y</mi> <mo>|</mo> </msqrt> </mrow>r,h分别称为“速度因子”和“滤波因子”r决定了跟踪速度的快慢,其值越大跟踪速度也越快但是其相应的误差也越大;h作为“滤波因子”起到了抑制误差的作用,调整参数可以满足跟踪速度与精度的要求;3)计及负荷电压突变影响缺额量的估计考虑到电压及频率影响因素的100%恒阻抗模型表示如下形式: <mrow> <msub> <mi>P</mi> <mi>L</mi> </msub> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>m</mi> </mrow> </munderover> <msub> <mi>P</mi> <mrow> <mi>L</mi> <mn>0</mn> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>·</mo> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>U</mi> <mi>i</mi> </msub> <msub> <mi>U</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>i</mi> </mrow> </msub> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>Q</mi> <mi>L</mi> </msub> <mo>=</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mrow> <mi>i</mi> <mo>=</mo> <mi>m</mi> </mrow> </munderover> <msub> <mi>Q</mi> <mrow> <mi>L</mi> <mn>0</mn> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>·</mo> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>U</mi> <mi>i</mi> </msub> <msub> <mi>U</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>i</mi> </mrow> </msub> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>其中PL,QL分别为扰动后某时刻系统总的负荷有功和无功功率,PL0,i,QL0,i为扰动前第i个负荷的初始有功功率和无功功率,Ui为扰动后某时刻第i个负荷的电压值,U0,i为扰动前第i个负荷的电压值;系统实际的功率缺额量为: <mrow> <msub> <mi>P</mi> <mi>shed</mi> </msub> <mo>=</mo> <msub> <mi>P</mi> <mi>G</mi> </msub> <mo>-</mo> <msub> <mi>P</mi> <mrow> <mi>L</mi> <mn>0</mn> </mrow> </msub> <mo>=</mo> <msub> <mi>P</mi> <mi>def</mi> </msub> <mo>+</mo> <munderover> <mi>Σ</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>P</mi> <mrow> <mi>L</mi> <mn>0</mn> <mo>,</mo> <mi>i</mi> </mrow> </msub> <mo>·</mo> <mo>[</mo> <msup> <mrow> <mo>(</mo> <mfrac> <msub> <mi>U</mi> <mi>i</mi> </msub> <msub> <mi>U</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>i</mi> </mrow> </msub> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mn>1</mn> <mo>]</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>式中PL0,i为第i个负荷的有功功率,PL0为扰动前总的负荷的初始有功功率,PG为扰动初始时刻系统所有发电机的输出的有功功率,Pshed为计及电压影响因素的实际缺额量。 |