发明名称 适用于SCR静止变频器的抽水蓄能机组低速位置检测方法
摘要 适用于SCR静止变频器的抽水蓄能机组低速位置检测方法。本发明涉及一种适用于SCR静止变频器的抽水蓄能机组低速位置检测方法。现有的低速阶段位置检测方法无法在抽水蓄能机组机端电压幅值低采样电压谐波含量大的情况下准确计算出转子的位置。一种适用于SCR静止变频器的抽水蓄能机组低速位置检测方法,通过对抽水蓄能机组定子三相绕组机端电压检测判断何时启动;待启动后将AD转换器采样的三相电压变换到两相静止坐标系下电压;经过改进型磁链观测器求出两相静止坐标系下第k时刻的定子绕组磁链值;再利用反正切求出转子位置角度。本发明应用于抽水蓄能机组低速位置检测。
申请公布号 CN104393809A 申请公布日期 2015.03.04
申请号 CN201410680857.2 申请日期 2014.11.24
申请人 哈尔滨工业大学 发明人 段建东;毕恺韬;孙力;国海峰;赵克;吴凤江;安群涛
分类号 H02P21/13(2006.01)I 主分类号 H02P21/13(2006.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 一种适用于SCR静止变频器的抽水蓄能机组低速位置检测方法,其特征在于:所述抽水蓄能机组低速位置检测方法通过以下步骤实现:步骤一、判断抽水蓄能机组定子通电电压值:利用AD转换器对抽水蓄能机组机端电压幅值连续采样,得N组机端电压值数据,将N组机端电压值数据等分成两组,分别对两组机端电压值数据求取平均值,将两组机端电压值数据的平均值进行对比,将较高的机端电压值数据的平均值确定为机端电压峰值,当机端电压峰值大于或等于2%的额定电压峰值时,进行第二步的PARK变换,否则继续进行机端电压幅值连续采样;步骤二、对AD转换器获得的N组机端电压值数据进行PARK变换:将机组定子三相相电压的电压e<sub>a</sub>(k)、电压e<sub>b</sub>(k)及电压e<sub>c</sub>(k)利用公式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>e</mi><mi>&alpha;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>e</mi><mi>&beta;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>=</mo><msqrt><mfrac><mn>2</mn><mn>3</mn></mfrac></msqrt><mfenced open='[' close=']'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mo>-</mo><mn>1</mn><mo>/</mo><mn>2</mn></mtd><mtd><mo>-</mo><mn>1</mn><mo>/</mo><mn>2</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msqrt><mn>3</mn></msqrt><mo>/</mo><mn>2</mn></mtd><mtd><mo>-</mo><msqrt><mn>3</mn></msqrt><mo>/</mo><mn>2</mn></mtd></mtr></mtable></mfenced><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>e</mi><mi>A</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>e</mi><mi>B</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>e</mi><mi>C</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000615029540000011.GIF" wi="1975" he="267" /></maths>转换成为α、β两相坐标系下电压e<sub>α</sub>(k)、e<sub>β</sub>(k);其中,k为获取电压值的时刻;步骤三、对步骤二获得的两相坐标系下电压e<sub>α</sub>(k)、e<sub>β</sub>(k)滤除高频谐波:利用由公式:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>e</mi><mi>&alpha;</mi><mo>'</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>&omega;</mi><mrow><mi>c</mi><mn>1</mn></mrow></msub><mrow><mi>s</mi><mo>+</mo><msub><mi>&omega;</mi><mrow><mi>c</mi><mn>1</mn></mrow></msub></mrow></mfrac><msub><mi>e</mi><mi>&alpha;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000615029540000012.GIF" wi="1693" he="156" /></maths>和<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>e</mi><mi>&beta;</mi><mo>'</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>&omega;</mi><mrow><mi>c</mi><mn>1</mn></mrow></msub><mrow><mi>s</mi><mo>+</mo><msub><mi>&omega;</mi><mrow><mi>c</mi><mn>1</mn></mrow></msub></mrow></mfrac><msub><mi>e</mi><mi>&beta;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000615029540000013.GIF" wi="1921" he="154" /></maths>形成的改进型磁链观测器分别计算出各个时刻下机组定子三相绕组在α、β两相坐标系下的磁链值ψ<sub>α</sub>(k)和ψ<sub>β</sub>(k);式中,e'<sub>α</sub>(k)为两相坐标系下电压e<sub>α</sub>(k)经过低通滤波器后的电压值,e'<sub>β</sub>(k)为两相坐标系下电压e<sub>β</sub>(k)经过低通滤波器后的电压值,s为拉普拉斯算子,ω<sub>c1</sub>为前级低通滤波器的截止频率;步骤四、对步骤三得到的低通滤波器后的电压值e'<sub>α</sub>(k)和e'<sub>β</sub>(k)进行电压补偿:利用公式:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>e</mi><mi>&alpha;</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>e</mi><mi>&alpha;</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mfrac><msub><mi>&omega;</mi><mrow><mi>c</mi><mn>2</mn></mrow></msub><mrow><msub><mi>&omega;</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><msubsup><mi>e</mi><mi>&beta;</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000615029540000014.GIF" wi="1763" he="154" /></maths>和<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>e</mi><mi>&beta;</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>e</mi><mi>&beta;</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mfrac><msub><mi>&omega;</mi><mrow><mi>c</mi><mn>2</mn></mrow></msub><mrow><msub><mi>&omega;</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><msubsup><mi>e</mi><mi>&alpha;</mi><mo>&prime;</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000615029540000015.GIF" wi="1945" he="154" /></maths>来求取补偿电压值;式中:<img file="FDA0000615029540000021.GIF" wi="125" he="79" />为经过补偿环节后获取的电压值,<img file="FDA0000615029540000022.GIF" wi="136" he="94" />为经过补偿环节后获取的电压值,ω<sub>c2</sub>为后级低通滤波器的截止频率,ω<sub>e</sub>为前一时刻的机组定子机组角频率;步骤五、对步骤四获得的补偿环节后的电压值<img file="FDA0000615029540000023.GIF" wi="122" he="90" />和<img file="FDA0000615029540000024.GIF" wi="149" he="105" />利用公式:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>&Psi;</mi><mi>&alpha;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>s</mi><mo>+</mo><msub><mi>&omega;</mi><mrow><mi>c</mi><mn>2</mn></mrow></msub></mrow></mfrac><msubsup><mi>e</mi><mi>&alpha;</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000615029540000025.GIF" wi="1992" he="165" /></maths>和<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>&Psi;</mi><mi>&beta;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>s</mi><mo>+</mo><msub><mi>&omega;</mi><mrow><mi>c</mi><mn>2</mn></mrow></msub></mrow></mfrac><msubsup><mi>e</mi><mi>&beta;</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000615029540000026.GIF" wi="1947" he="162" /></maths>求取第k时刻机组定子在静止坐标系下的磁链;步骤六、对步骤五获得的磁链值通过反正切方法,即利用公式:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>&theta;</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>arctan</mi><mrow><mo>(</mo><mfrac><mrow><msub><mi>&psi;</mi><mi>&beta;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>&psi;</mi><mi>&alpha;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000615029540000027.GIF" wi="1985" he="171" /></maths>求取第k时刻转子位置角度θ<sub>e</sub>(k);步骤七、通过步骤六获得转子位置角度θ<sub>e</sub>(k)通过公式:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>&omega;</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>d&theta;</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac><mo>=</mo><mfrac><mrow><msub><mi>&psi;</mi><mi>&alpha;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>e</mi><mi>&beta;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>&psi;</mi><mi>&beta;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><msub><mi>e</mi><mi>&alpha;</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><msubsup><mi>&psi;</mi><mi>&alpha;</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&psi;</mi><mi>&beta;</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000615029540000028.GIF" wi="1974" he="177" /></maths>进一步获得求取第k时刻低速转子转速,即求取第k时刻机组角频率ω<sub>e</sub>(k)。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号