发明名称 一种基于积分法的锁相环控制方法
摘要 本发明涉及属于电力系统运行控制方法技术领域的一种基于积分法的锁相环控制方法。采用park变换的方法将原信号分解到d轴、q轴上,并且把电压分量分别用正负序相量表示;分别对d轴和q轴电压分量做时间t的积分,交叉运算消除不平衡分量;然后再消去一阶函数;将V<sub>d+</sub>、V<sub>q+</sub>合成为测量电压并与电压参考值V<sub>ref</sub>进行比较,将差值输入到PI运算环节,得到用于锁相环控制的相位信号<img file="DDA0000456142870000011.GIF" wi="32" he="39" />;这样,即得到基于积分法的锁相环控制方法。本发明与现有方法相比,平衡条件或不平衡条件均适用;不会引入时间延迟;如果原信号中有高次谐波时,采用基于微分的锁相环控制方法会放大这些谐波,造成不好的影响,本发明则不会放大这些谐波。在电力系统运行控制领域具有重要积极意义。
申请公布号 CN103825605A 申请公布日期 2014.05.28
申请号 CN201410014794.7 申请日期 2014.01.13
申请人 国家电网公司;华北电力大学;中国电力科学研究院;国网上海市电力公司 发明人 姚蜀军;鲍明然;韩民晓;汪燕
分类号 H03L7/08(2006.01)I 主分类号 H03L7/08(2006.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 薄观玖
主权项 1.一种基于积分法的锁相环控制方法,其特征在于,该方法步骤为:(1)采用park变换将三相不平衡电压从abc坐标变换到dq坐标下,并作正负序分解;将电网三相不平衡电压进行park变换,从abc坐标转换为dq坐标下分量V<sub>d</sub>和V<sub>q</sub>,V<sub>d</sub>和V<sub>q</sub>分别表示电压相量在d轴和q轴上的分量,V<sub>ref</sub>表示电压参考值,<img file="FDA0000456142840000011.GIF" wi="46" he="57" />表示锁相环输出的相位信号,M为调制比,δ为移相角;因为是不平衡量,写成正负序分量相加的形式,在不平衡情况下,存在正序分量V<sub>+</sub>和负序分量V<sub>-</sub>;正序分量V<sub>+</sub>以ω的角速度逆时针旋转,与坐标轴同步,相对静止,与d轴夹角为θ<sub>+</sub>,正序分量V<sub>+</sub>分解得到:<maths num="0001"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>V</mi><mrow><mi>d</mi><mo>+</mo></mrow></msub><mo>=</mo><msub><mi>V</mi><mo>+</mo></msub><mi>cos</mi><msub><mi>&theta;</mi><mo>+</mo></msub></mtd></mtr><mtr><mtd><msub><mi>V</mi><mrow><mi>q</mi><mo>+</mo></mrow></msub><mo>=</mo><msub><mi>V</mi><mo>+</mo></msub><mi>sin</mi><msub><mi>&theta;</mi><mo>+</mo></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo></mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,V<sub>d+</sub>是电压相量在d轴上分量V<sub>d</sub>的正序分量,V<sub>q+</sub>是电压相量在q轴上分量V<sub>q</sub>的正序分量,V<sub>+</sub>是电压相量正序分量,θ<sub>+</sub>是V<sub>+</sub>与d轴的夹角;负序分量V<sub>-</sub>以ω的角速度顺时针旋转,相对于坐标轴是以2ω的角速度顺时针旋转,与d轴夹角为θ<sub>-</sub>,负序分量V<sub>-</sub>分解得到<maths num="0002"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>V</mi><mrow><mi>d</mi><mo>-</mo></mrow></msub><mo>=</mo><msub><mi>V</mi><mo>-</mo></msub><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&omega;t</mi><mo>+</mo><msub><mi>&theta;</mi><mo>-</mo></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>V</mi><mrow><mi>q</mi><mo>-</mo></mrow></msub><mo>=</mo><mo>-</mo><msub><mi>V</mi><mo>-</mo></msub><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&omega;t</mi><mo>+</mo><msub><mi>&theta;</mi><mo>-</mo></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo></mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,V<sub>d-</sub>是电压相量在d轴上分量V<sub>d</sub>的负序分量,V<sub>q-</sub>是电压相量在q轴上分量V<sub>q</sub>的正序分量,V<sub>-</sub>是电压相量负序分量,θ<sub>-</sub>是V<sub>-</sub>与d轴的夹角;由以上两式得到在三相不平衡情况下,电压相量在dq坐标下的表达形式:<maths num="0003"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>v</mi><mi>d</mi></msub><mo>=</mo><msub><mi>v</mi><mrow><mi>d</mi><mo>+</mo></mrow></msub><mo>+</mo><msub><mi>v</mi><mrow><mi>d</mi><mo>-</mo></mrow></msub><mo>=</mo><msub><mi>V</mi><mo>+</mo></msub><mi>cos</mi><msub><mi>&theta;</mi><mo>+</mo></msub><mo>+</mo><msub><mi>V</mi><mo>-</mo></msub><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&omega;t</mi><mo>+</mo><msub><mi>&theta;</mi><mo>-</mo></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>v</mi><mi>q</mi></msub><mo>=</mo><msub><mi>v</mi><mrow><mi>q</mi><mo>+</mo></mrow></msub><mo>+</mo><msub><mi>v</mi><mrow><mi>q</mi><mo>-</mo></mrow></msub><mo>=</mo><msub><mi>V</mi><mo>+</mo></msub><mi>sin</mi><msub><mi>&theta;</mi><mo>+</mo></msub><mo>-</mo><msub><mi>V</mi><mo>-</mo></msub><mi>sin</mi><mrow><mo>(</mo><mn>2</mn><mi>&omega;t</mi><mo>+</mo><msub><mi>&theta;</mi><mo>-</mo></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo></mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>由上式可以看出,在三相不平衡情况下,在d轴和q轴上出现了两倍频的扰动量;(2)对电压相量的d轴和q轴分量分别作积分运算,通过积分方法消去两倍频分量;然后通过比较环节和PI运算环节得到同步相角<img file="FDA0000456142840000015.GIF" wi="73" he="71" />所述通过积分方法消去两倍频分量的方法为:对式1)中电压相量V<sub>d</sub>和V<sub>q</sub>分别对时间t作积分,得<maths num="0004"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mtext>&Integral;</mtext><msub><mi>V</mi><mi>d</mi></msub><mi>dt</mi><mo>=</mo><msub><mi>V</mi><mo>+</mo></msub><mi>cos</mi><msub><mi>&theta;</mi><mo>+</mo></msub><mi>t</mi><mo>+</mo><mfrac><msub><mi>V</mi><mo>-</mo></msub><mrow><mn>2</mn><mi>&omega;</mi></mrow></mfrac><mi>sin</mi><mrow><mo>(</mo><mn>2</mn><mi>&omega;t</mi><mo>+</mo><msub><mi>&theta;</mi><mo>-</mo></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mrow><mo>&Integral;</mo><mi>V</mi></mrow><mi>q</mi></msub><mi>dt</mi><mo>=</mo><msub><mi>V</mi><mo>+</mo></msub><mi>sin</mi><msub><mi>&theta;</mi><mo>+</mo></msub><mi>t</mi><mo>+</mo><mfrac><msub><mi>V</mi><mo>-</mo></msub><mrow><mn>2</mn><mi>&omega;</mi></mrow></mfrac><mi>cos</mi><mrow><mo>(</mo><mn>2</mn><mi>&omega;t</mi><mo>+</mo><msub><mi>&theta;</mi><mo>-</mo></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo></mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>对于d轴分量:V<sub>d</sub>-2ω∫V<sub>q</sub>dt=V<sub>+</sub>cosθ<sub>+</sub>-2ωV<sub>+</sub>sinθ<sub>+</sub>t,    5)结果中存在一阶函数;将其消去,设目标函数为:f(t)=V<sub>+</sub>cosθ<sub>+</sub>-2ωV<sub>+</sub>sinθ<sub>+</sub>t,    6)仿真步长为T,在时间轴上取两点t<sub>1</sub>,t<sub>2</sub>,其中t<sub>2</sub>=t<sub>1</sub>+T,对应函数:f(t<sub>1</sub>)=V<sub>+</sub>cosθ<sub>+</sub>-2ωV<sub>+</sub>sinθ<sub>+</sub>t<sub>1</sub>f(t<sub>2</sub>)=V<sub>+</sub>cosθ<sub>+</sub>-2ωV<sub>+</sub>sinθ<sub>+</sub>t<sub>2</sub>,    7)设斜率为K,则,<maths num="0005"><![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><mi>K</mi><mo>=</mo><mo>-</mo><mn>2</mn><mi>&omega;</mi><msub><mi>V</mi><mo>+</mo></msub><mi>sin</mi><msub><mi>&theta;</mi><mo>+</mo></msub><mo>=</mo><mfrac><mrow><mi>f</mi><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow><mi>T</mi></mfrac></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mrow><mo>-</mo><mn>2</mn><mi>&omega;</mi><msub><mi>V</mi><mo>+</mo></msub><mi>sin</mi><msub><mi>&theta;</mi><mo>+</mo></msub><msub><mi>t</mi><mn>2</mn></msub><mo>+</mo><mn>2</mn><mi>&omega;</mi><msub><mi>V</mi><mo>+</mo></msub><mi>sin</mi><msub><mi>&theta;</mi><mo>+</mo></msub><msub><mi>t</mi><mn>1</mn></msub></mrow><mi>T</mi></mfrac></mtd></mtr></mtable></mfenced><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo></mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>消去一阶函数,得:V<sub>d+</sub>=f(t)-Kt=V<sub>+</sub>cosθ<sub>+</sub>-2ωV<sub>+</sub>sinθ<sub>+</sub>t-Kt=V<sub>+</sub>cosθ<sub>+</sub>,    9);对于q轴分量:V<sub>q</sub>+2ω∫V<sub>d</sub>dt=V<sub>+</sub>sinθ<sub>+</sub>+2ωV<sub>+</sub>cosθ<sub>+</sub>t    10)结果中存在一阶函数,将其消去,设目标函数为:g(t)=V<sub>+</sub>sinθ<sub>+</sub>+2ωV<sub>+</sub>cosθ<sub>+</sub>t,    11)仿真步长为T,在时间轴上取两点t<sub>1</sub>,t<sub>2</sub>,其中t<sub>2</sub>=t<sub>1</sub>+T,对应函数g(t<sub>1</sub>)=V<sub>+</sub>sinθ<sub>+</sub>+2ωV<sub>+</sub>cosθ<sub>+</sub>t<sub>1</sub>g(t<sub>2</sub>)=V<sub>+</sub>sinθ<sub>+</sub>+2ωV<sub>+</sub>cosθ<sub>+</sub>t<sub>2</sub>,    12)设斜率为K,则,<maths num="0006"><![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><mi>K</mi><mo>=</mo><mn>2</mn><mi>&omega;</mi><msub><mi>V</mi><mo>+</mo></msub><mi>cos</mi><msub><mi>&theta;</mi><mo>+</mo></msub><mo>=</mo><mfrac><mrow><mi>g</mi><mrow><mo>(</mo><msub><mi>t</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><mi>g</mi><mrow><mo>(</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow><mi>T</mi></mfrac></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mrow><mn>2</mn><mi>&omega;</mi><msub><mi>V</mi><mo>+</mo></msub><mi>cos</mi><msub><mi>&theta;</mi><mo>+</mo></msub><msub><mi>t</mi><mn>2</mn></msub><mo>-</mo><mn>2</mn><mi>&omega;</mi><msub><mi>V</mi><mo>+</mo></msub><mi>cos</mi><msub><mi>&theta;</mi><mo>+</mo></msub><msub><mi>t</mi><mn>1</mn></msub></mrow><mi>T</mi></mfrac></mtd></mtr></mtable></mfenced><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo></mo><mn>13</mn><mo>)</mo></mrow></mrow></math>]]></maths>消去一阶函数,得:V<sub>q+</sub>=g(t)-Kt=V<sub>+</sub>sinθ<sub>+</sub>+2ωV<sub>+</sub>cosθ<sub>+</sub>t-Kt=V<sub>+</sub>sinθ<sub>+</sub>,    14);上式中正弦、余弦分量分别表示电压相量在q轴和d轴上的投影,也即是在q轴和d轴上的分量;将式9)、14)得到的V<sub>d+</sub>和V<sub>q+</sub>合成为测量电压,并与电压参考值V<sub>ref</sub>通过比较环节进行比较,得到差值,将差值输入到PI运算环节消除误差,得到相位信号<img file="FDA0000456142840000032.GIF" wi="70" he="73" />所述比较环节为将测量值与参考值相比较得出差值;所述PI运算环节即比例积分运算环节;这样,即得到基于积分法的锁相环控制方法。
地址 100031 北京市西城区西长安街86号