发明名称 光伏并网逆变器自适应准PRD控制方法
摘要 光伏并网逆变器自适应准PRD控制方法,涉及光伏并网逆变器控制技术领域。本发明是为了解决现有光伏并网逆变器对电网阻抗变化适应性差,和弱电网下光伏并网系统的稳定性差的问题。本发明所述的通过扰动电流注入法向电网中注入特定频率的扰动电流,然后基于闭环极点参数配置的原则得到电流控制器参数与电网阻抗的关系式,当电网电阻R<sub>g</sub>和电网电感L<sub>g</sub>发生变化时,实现对控制器参数K<sub>p</sub>、K<sub>r</sub>、K<sub>d</sub>的自适应调节,以适应不断变化的电网阻抗,提高了光伏并网逆变器控制系统的适应性与稳定性。它可用于对光伏并网逆变器的控制。
申请公布号 CN104810859A 申请公布日期 2015.07.29
申请号 CN201510278698.8 申请日期 2015.05.27
申请人 哈尔滨工业大学 发明人 王卫;刘桂花;刘鸿鹏;曹小娇;吴辉;徐殿国
分类号 H02J3/38(2006.01)I;H02J3/01(2006.01)I 主分类号 H02J3/38(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 岳泉清
主权项 光伏并网逆变器自适应准PRD控制方法,其特征在于,它包括以下步骤:步骤一、通过逆变器以扰动电流注入法向电网中注入设定频率的扰动电流,然后利用检测元件检测出并网点处的电压和电流响应,再结合信号处理技术,分离出所注入特定次谐波分量的模值,根据该模值获得等效的电网电阻R<sub>g</sub>和电网电感L<sub>g</sub>,步骤二、将逆变器桥臂输出的正弦脉宽调制电压等效为一个电压源,并根据步骤一中获得的电网电阻R<sub>g</sub>和电网电感L<sub>g</sub>,建立LC型单相光伏并网逆变器并网电流控制回路的受控对象的传递函数G<sub>o</sub>(s),对该传递函数G<sub>o</sub>(s)进行离散化,得到相应的离散化表达式G<sub>o</sub>(z<sup>‑1</sup>),步骤三、采用准PRD控制器作为电流内环的控制器,准PRD控制器的传递函数表达式为<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>G</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>K</mi><mi>p</mi></msub><mo>+</mo><mfrac><mrow><mn>2</mn><msub><mi>K</mi><mi>r</mi></msub><msub><mi>&omega;</mi><mi>c</mi></msub><mi>s</mi></mrow><mrow><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><msub><mi>&omega;</mi><mi>c</mi></msub><mi>s</mi><mo>+</mo><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup></mrow></mfrac><mo>+</mo><mfrac><mrow><msub><mi>K</mi><mi>d</mi></msub><mi>s</mi></mrow><mrow><mn>1</mn><mo>+</mo><mi>&tau;s</mi></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000725812120000011.GIF" wi="769" he="147" /></maths>其中,K<sub>p</sub>为准PRD控制器比例系数,K<sub>r</sub>为准PRD控制器谐振系数,K<sub>d</sub>为准PRD控制器微分系数,ω<sub>0</sub>为基波对应的角频率,ω<sub>c</sub>为截止角频率,<img file="FDA0000725812120000012.GIF" wi="138" he="132" />为低通滤波环节,τ为滤波环节中的延迟时间;对准PRD控制器的传递函数G<sub>c</sub>(s)进行离散化,获得其离散化表达式G<sub>c</sub>(z<sup>‑1</sup>),对离散化表达式G<sub>c</sub>(z<sup>‑1</sup>)进行变换得到差分方程:u(k)=b<sub>0</sub>e(k)+b<sub>1</sub>e(k‑1)+b<sub>2</sub>e(k‑2)+b<sub>3</sub>e(k‑3)‑a<sub>1</sub>u(k‑1)‑a<sub>2</sub>u(k‑2)‑a<sub>3</sub>u(k‑3),其中,k为采样步数,u(k)为第k个采样周期内,准PRD数字控制器的输出量,e(k)为控制器的输入量,也即给定电流与实际电流的差值,a<sub>i</sub>和b<sub>j</sub>为系统参数,i=1,2,3;j=1,2,3;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>a</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><mn>3</mn><mi>T</mi><mo>+</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>+</mo><mrow><mo>(</mo><mi>T</mi><mo>-</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mo>-</mo><mn>4</mn><mrow><mo>(</mo><mi>T</mi><mo>+</mo><mn>6</mn><mi>&tau;</mi><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><mi>T</mi><mo>+</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mrow><mo>(</mo><mn>4</mn><mo>+</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mo>+</mo><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000725812120000013.GIF" wi="1001" he="172" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>a</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><mn>3</mn><mi>T</mi><mo>-</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>-</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>+</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mo>-</mo><mn>4</mn><mrow><mo>(</mo><mi>T</mi><mo>-</mo><mn>6</mn><mi>&tau;</mi><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><mi>T</mi><mo>+</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mrow><mo>(</mo><mn>4</mn><mo>+</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mo>+</mo><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000725812120000014.GIF" wi="1007" he="164" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>a</mi><mn>3</mn></msub><mo>=</mo><mfrac><mrow><mrow><mo>(</mo><mi>T</mi><mo>-</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mrow><mo>(</mo><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>-</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mo>+</mo><mn>4</mn><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><mi>T</mi><mo>+</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mrow><mo>(</mo><mn>4</mn><mo>+</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mo>+</mo><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000725812120000015.GIF" wi="681" he="176" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>b</mi><mn>0</mn></msub><mo>=</mo><msub><mi>K</mi><mi>p</mi></msub><mo>+</mo><mfrac><mrow><mn>4</mn><msub><mi>K</mi><mi>r</mi></msub><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub></mrow><mrow><mn>4</mn><mo>+</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mo>+</mo><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>2</mn><msub><mi>K</mi><mi>d</mi></msub></mrow><mrow><mi>T</mi><mo>+</mo><mn>2</mn><mi>&tau;</mi></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000725812120000021.GIF" wi="754" he="147" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>b</mi><mn>1</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>K</mi><mi>p</mi></msub><mo>[</mo><mrow><mo>(</mo><mn>3</mn><mi>T</mi><mo>+</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>+</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>-</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mo>-</mo><mn>4</mn><mrow><mo>(</mo><mi>T</mi><mo>+</mo><mn>6</mn><mi>&tau;</mi><mo>)</mo></mrow><mo>]</mo><mo>+</mo><msub><mi>K</mi><mi>r</mi></msub><mo>&CenterDot;</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>-</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mo>+</mo><mn>2</mn><msub><mi>K</mi><mi>d</mi></msub><mrow><mo>(</mo><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>-</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mo>-</mo><mn>12</mn><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><mi>T</mi><mo>+</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mrow><mo>(</mo><mn>4</mn><mo>+</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mo>+</mo><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000725812120000022.GIF" wi="1945" he="170" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>b</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>K</mi><mi>p</mi></msub><mo>[</mo><mrow><mo>(</mo><mn>3</mn><mi>T</mi><mo>-</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>-</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>+</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mo>-</mo><mn>4</mn><mrow><mo>(</mo><mi>T</mi><mo>-</mo><mn>6</mn><mi>&tau;</mi><mo>)</mo></mrow><mo>]</mo><mo>-</mo><msub><mi>K</mi><mi>r</mi></msub><mo>&CenterDot;</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>+</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mo>+</mo><mn>2</mn><msub><mi>K</mi><mi>d</mi></msub><mrow><mo>(</mo><mrow><mn>12</mn><mo>-</mo></mrow><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mo>-</mo><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><mi>T</mi><mo>+</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mrow><mo>(</mo><mn>4</mn><mo>+</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mo>+</mo><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000725812120000023.GIF" wi="1944" he="168" /></maths><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>b</mi><mn>3</mn></msub><mo>=</mo><mfrac><mrow><msub><mi>K</mi><mi>p</mi></msub><mo>[</mo><mrow><mo>(</mo><mi>T</mi><mo>-</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>-</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>-</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mo>+</mo><mn>4</mn><mrow><mo>(</mo><mi>T</mi><mo>-</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mo>]</mo><mo>-</mo><msub><mi>K</mi><mi>r</mi></msub><mo>&CenterDot;</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>-</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mo>-</mo><mn>2</mn><msub><mi>K</mi><mi>d</mi></msub><mrow><mo>(</mo><mn>4</mn><mo>-</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mo>+</mo><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>)</mo></mrow></mrow><mrow><mrow><mo>(</mo><mi>T</mi><mo>+</mo><mn>2</mn><mi>&tau;</mi><mo>)</mo></mrow><mrow><mo>(</mo><mn>4</mn><mo>+</mo><mn>4</mn><mi>T</mi><msub><mi>&omega;</mi><mi>c</mi></msub><mo>+</mo><msup><mi>T</mi><mn>2</mn></msup><msubsup><mi>&omega;</mi><mn>0</mn><mn>2</mn></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000725812120000024.GIF" wi="1869" he="167" /></maths>步骤四、根据准PRD控制器理想闭环极点的位置,获得理想的闭环特征方程A<sub>m</sub>(z<sup>‑1</sup>),根据步骤二得到的受控对象传递函数G<sub>o</sub>(z<sup>‑1</sup>)和步骤三得到的准PRD控制器传递函数G<sub>c</sub>(z<sup>‑1</sup>),获得系统实际的闭环特征方程A<sub>c</sub>(z<sup>‑1</sup>),令A<sub>c</sub>(z<sup>‑1</sup>)=A<sub>m</sub>(z<sup>‑1</sup>),获得准PRD控制器的控制参数K<sub>p</sub>、K<sub>r</sub>、K<sub>d</sub>与电网电阻R<sub>g</sub>和电网电感L<sub>g</sub>的关系式,当电网电阻R<sub>g</sub>和电网电感L<sub>g</sub>发生变化时,逆变器根据步骤三中的差分方程实现准PRD控制器的输出,以实现自适应准PRD控制器对光伏并网逆变器的控制。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号