发明名称 一种交直流混合配电网的电压稳定分岔分析方法
摘要 本发明公开了一种交直流混合配电网的电压稳定分岔分析方法,该方法将交直流混合配电网等效为三部分,基于微分代数方程设计相应的控制器,包括电压环、电流环和稳定器;通过设定系统参数和控制参数,实现基于有功功率参考值变化的电压稳定分岔分析;并可基于所得到的分岔图,确定出电压稳定的分岔拐点和极限功率,从而选出系统的稳定器最优参数,找出系统稳定裕度最大时对应的带宽数值,并可实现系统的渐进稳定分析、静态稳定分析、动态稳定分析和暂态稳定分析,得出稳定器参数与系统极限功率之间的对应关系,最终设定稳定器的反馈控制信号,控制系统运行在最大的稳定裕度范围内。
申请公布号 CN106026100A 申请公布日期 2016.10.12
申请号 CN201610590091.8 申请日期 2016.07.19
申请人 河南理工大学 发明人 韦延方;曾志辉;王晓卫;王帅;郑征;孙岩洲;肖记军
分类号 H02J3/02(2006.01)I 主分类号 H02J3/02(2006.01)I
代理机构 代理人
主权项 一种交直流混合配电网的电压稳定分岔分析方法,其特征在于包括以下步骤:步骤1:将交直流混合配电网等效为三部分,分别为:buck型dc‑dc变换器单元、LC滤波器单元、boost型dc‑dc变换器单元;其中,buck型dc‑dc变换器单元与LC滤波器单元通过电缆线连接,用于向直流母线提供稳定的电压;LC滤波器单元与boost型dc‑dc变换器单元通过电缆线连接,作为boost型dc‑dc变换器单元的输入滤波器,同时作为电缆线的漏感;步骤2:确定交直流混合配电网的微分代数方程,其计算公式为:<maths num="0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mfrac><mrow><msub><mi>di</mi><msub><mi>L</mi><mn>1</mn></msub></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><msub><mi>L</mi><mn>1</mn></msub></mfrac><mrow><mo>(</mo><msub><mi>u</mi><mi>s</mi></msub><mo>&CenterDot;</mo><msub><mi>V</mi><mi>e</mi></msub><mo>-</mo><msub><mi>V</mi><mi>s</mi></msub><mo>-</mo><msub><mi>R</mi><mn>1</mn></msub><mo>&CenterDot;</mo><msub><mi>i</mi><msub><mi>L</mi><mn>1</mn></msub></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><msub><mi>dV</mi><mi>s</mi></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><mi>C</mi></mfrac><mrow><mo>(</mo><msub><mi>i</mi><msub><mi>L</mi><mn>1</mn></msub></msub><mo>-</mo><msub><mi>i</mi><mrow><mi>d</mi><mi>c</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><msub><mi>di</mi><mrow><mi>d</mi><mi>c</mi></mrow></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><msub><mi>L</mi><mi>f</mi></msub></mfrac><mrow><mo>(</mo><msub><mi>V</mi><mi>s</mi></msub><mo>-</mo><msub><mi>V</mi><mrow><mi>c</mi><mi>f</mi></mrow></msub><mo>-</mo><msub><mi>R</mi><mi>f</mi></msub><mo>&CenterDot;</mo><msub><mi>i</mi><mrow><mi>d</mi><mi>c</mi></mrow></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><msub><mi>dV</mi><mrow><mi>c</mi><mi>f</mi></mrow></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><msub><mi>C</mi><mi>f</mi></msub></mfrac><mrow><mo>(</mo><msub><mi>i</mi><mrow><mi>d</mi><mi>c</mi></mrow></msub><mo>-</mo><msub><mi>i</mi><msub><mi>L</mi><mn>2</mn></msub></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><msub><mi>di</mi><msub><mi>L</mi><mn>2</mn></msub></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><msub><mi>L</mi><mn>2</mn></msub></mfrac><mrow><mo>(</mo><msub><mi>V</mi><mrow><mi>c</mi><mi>f</mi></mrow></msub><mo>-</mo><msub><mi>R</mi><mn>2</mn></msub><mo>&CenterDot;</mo><msub><mi>i</mi><msub><mi>L</mi><mn>2</mn></msub></msub><mo>-</mo><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>u</mi><mi>l</mi></msub></mrow><mo>)</mo><msub><mi>V</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mfrac><mrow><msub><mi>dV</mi><mn>0</mn></msub></mrow><mrow><mi>d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><msub><mi>C</mi><mn>0</mn></msub></mfrac><mrow><mo>(</mo><mo>(</mo><mrow><mn>1</mn><mo>-</mo><msub><mi>u</mi><mi>l</mi></msub></mrow><mo>)</mo><msub><mi>i</mi><msub><mi>L</mi><mn>2</mn></msub></msub><mo>-</mo><mfrac><msub><mi>V</mi><mn>0</mn></msub><msub><mi>R</mi><mn>0</mn></msub></mfrac><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FSA0000132717350000011.GIF" wi="789" he="987" /></maths>其中,t表示时间变量,d表示求导;u<sub>s</sub>表示buck型dc‑dc变换器单元中电力电子开关器件的脉冲指令;u<sub>l</sub>表示boost型dc‑dc变换器单元中电力电子开关器件的脉冲指令;R<sub>f</sub>、L<sub>f</sub>和C<sub>f</sub>分别表示LC滤波器单元的电阻、电感和电容;C表示直流母线的电容;L<sub>1</sub>、R<sub>1</sub>分别表示buck型dc‑dc变换器单元的串联电感和电阻;L<sub>2</sub>、R<sub>2</sub>分别表示boost型dc‑dc变换器单元的电感和电阻;V<sub>e</sub>表示buck型dc‑dc变换器单元的等效电源电压;<img file="FSA0000132717350000023.GIF" wi="154" he="73" />和i<sub>dc</sub>分别表示buck型dc‑dc变换器单元电流、boost型dc‑dc变换器单元电流和LC滤波器单元电流;V<sub>s</sub>、V<sub>cf</sub>和V<sub>0</sub>分别表示直流母线电压、LC滤波器单元电压和负载的输出电压;C<sub>0</sub>、R<sub>0</sub>分别表示负载的电容和电阻;步骤3:设计buck型dc‑dc变换器单元的控制器,包括电压环、电流环和稳定器;其中,电压环采用线性的分数阶PI控制器进行校正;电流环采用间接的滑膜控制方法,其控制语句由等效控制规则生成;步骤4:设计boost型dc‑dc变换器单元的控制器,其由RT‑Lab的数字控制卡实现,用于调整系统输出的有功功率P;其中,有功功率参考值P<sub>ref</sub>的计算公式为:P<sub>ref</sub>=I<sub>ref</sub>·V<sub>ref</sub>其中,I<sub>ref</sub>、V<sub>ref</sub>分别表示系统测量到的电流、电压的参考值;I<sub>ref</sub>由电压环确定;步骤5:确定控制变量的采样速率周期T,其计算公式为:<maths num="0002"><math><![CDATA[<mrow><mi>T</mi><mo>=</mo><mfrac><mn>1</mn><msub><mi>f</mi><mi>s</mi></msub></mfrac></mrow>]]></math><img file="FSA0000132717350000021.GIF" wi="152" he="137" /></maths>其中,f<sub>s</sub>表示变换器的开关频率;步骤6:电压环表达式C<sub>v</sub>(s)的计算公式为:<maths num="0003"><math><![CDATA[<mrow><msub><mi>C</mi><mi>v</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>K</mi><mrow><mi>p</mi><mi>v</mi></mrow></msub><mo>+</mo><mfrac><msub><mi>K</mi><mrow><mi>i</mi><mi>v</mi></mrow></msub><msup><mi>s</mi><mi>&lambda;</mi></msup></mfrac></mrow>]]></math><img file="FSA0000132717350000022.GIF" wi="394" he="123" /></maths>其中,s表示自变量,K<sub>pv</sub>表示比例增益,K<sub>iv</sub>表示积分系数,λ表示积分阶次;步骤7:设计稳定器,其由基于比例校正环节的高通滤波器构成,该高通滤波器在系统电压发生波动时,通过产生占空比信号,阻尼系统的振荡;步骤8:设定系统参数和控制参数,包括直流输入电压、buck型dc‑dc变换器单元的电阻和电感、直流母线的电容、滤波器的电阻、滤波器的电容和电感、boost型dc‑dc变换器单元的电阻和电感、输出电容、参考电压、开关频率、稳定器的截止频率、电压和电流的分数阶比例增益和积分系数;步骤9:基于有功功率参考值P<sub>ref</sub>变化的电压稳定分岔分析;首先设定buck型dc‑dc变换器单元控制器的控制参数,接着初始化稳定器,得出V<sub>s</sub>、V<sub>cf</sub>的分岔图;根据分岔图,并确定出电压稳定的分岔拐点,得出此时的极限功率P<sub>0</sub>;在P<sub>0</sub>点,系统失去电压的稳定性;当有功功率继续增大、以致超过P<sub>0</sub>时,观察极限环振荡的情况,确定是否出现Hopf分岔;步骤10:调整稳定器的参数,观察极限功率点随有功功率参考值P<sub>ref</sub>变化的轨迹,找出最大的极限功率P<sub>0max</sub>;并选定系统出现P<sub>0max</sub>时所对应的稳定器参数为最优参数;步骤11:基于带宽控制的电压稳定分岔分析;首先设定有功功率参考值P<sub>ref</sub>,接着初始化稳定器,通过调整buck型dc‑dc变换器单元控制器的控制参数,设定电流环控制信号的带宽,使其高于电压环控制信号的带宽,基于电流环控制的系统动态特性分析,得出系统的分岔图,最后观察系统稳定裕度随着带宽改变而发生的变化轨迹,找出系统稳定裕度最大时对应的带宽数值;步骤12:渐进稳定性分析,通过计算系统雅可比矩阵的特征根,判断系统的不稳定状态,识别出系统的稳定运行点和不稳定运行点,并观察buck型dc‑dc变换器单元和boost型dc‑dc变换器单元中控制信号的带宽发生变化时,电源经直流网络向负载传输的有功功率变化情况;步骤13:基于稳定器的不同参数设置,观察系统的静态稳定情况、动态稳定情况和暂态稳定情况,得出稳定器参数与系统极限功率之间的对应关系,最终设定稳定器的反馈控制信号,控制系统运行在最大的稳定裕度范围内。
地址 454003 河南省焦作市高新区世纪大道2001号河南理工大学