发明名称 一种基于广域相量测量系统的功率缺额估计方法
摘要 本发明是一种基于广域相量测量系统的功率缺额估计方法,其特点是,包括基于惯量中心坐标的功率缺额估计、频率及其变化率跟踪算法和计及负荷电压突变影响缺额量的估计步骤,将最速跟踪微分器应用于电力系统频率及其变化率的跟踪与计算中,该方法在传统功率缺额的估算的基础上同时计及了电压变化对功率缺额估计的影响,具有科学合理,简便易行,计算准确,工程应用价值高等优点。且通过仿真算例分析验证了本发明在处理低频减载问题时具有较高的工程实际应用价值。
申请公布号 CN102508061B 申请公布日期 2014.02.12
申请号 CN201110323674.1 申请日期 2011.10.22
申请人 东北电力大学;中国大唐集团科学技术研究院有限公司 发明人 蔡国伟;邢亮;王春生;李振新;田世杰;孔令国;章昊;刘铖;李鹏飞;姜又滔
分类号 G01R31/00(2006.01)I 主分类号 G01R31/00(2006.01)I
代理机构 吉林市达利专利事务所 22102 代理人 陈传林
主权项 一种基于广域相量测量系统的功率缺额估计方法,其特征是,它包括以下步骤:1)基于惯量中心坐标的功率缺额估计 <mrow> <msub> <mi>&omega;</mi> <mi>COI</mi> </msub> <mo>=</mo> <mfrac> <mrow> <munderover> <mi>&Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>H</mi> <mi>i</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>&omega;</mi> <mi>i</mi> </msub> </mrow> <mrow> <munderover> <mi>&Sigma;</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&Delta;&omega;</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>&Sigma;</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>&CenterDot;</mo> <msub> <mi>S</mi> <mi>i</mi> </msub> </mrow> <mrow> <munderover> <mi>&Sigma;</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>&Sigma;</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&Delta;&omega;</mi> <mi>dt</mi> </mfrac> <munderover> <mi>&Sigma;</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>&CenterDot;</mo> <mfrac> <mi>a</mi> <mi>&delta;</mi> </mfrac> <mo>,</mo> </mtd> <mtd> <mo>|</mo> <mi>a</mi> <mo>|</mo> <mo>&le;</mo> <mi>&delta;</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>&delta;</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>&delta;</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>&delta;</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>&le;</mo> <msub> <mi>&delta;</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>&delta;</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>&Sigma;</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>&CenterDot;</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>&Sigma;</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>&CenterDot;</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>&Sigma;</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>&CenterDot;</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为计及电压影响因素的实际缺额量。
地址 132012 吉林省吉林市船营区长春路169号