发明名称 一种基于反求法的变频供水系统运行效率在线检测方法
摘要 本发明提供了一种基于反求法的变频供水系统运行效率在线检测方法,包括:建立输出功率、频率扰动量和水压变化量之间的数学模型及约束条件,形成供水系统运行效率在线检测需要的数学模型;在稳态工况下进行频率小信号扰动,采样管网压力值,根据输出功率数学模型反求出稳态时轴输出功率,进而求得变频供水系统效率。本发明可实现在线输出功率检测,无需流量传感器,节约系统安装调试所需时间和成本,使得系统结构更加简单,系统成本更低。本发明可有效保护电机和变频器低频运行引起的效率低下故障,提高系统的寿命和可靠性,为水泵电机安全、高效运行提供可靠保证。
申请公布号 CN103487187A 申请公布日期 2014.01.01
申请号 CN201310410409.6 申请日期 2013.09.10
申请人 温州大学 发明人 彭志辉;李峰平;付培红;曹宇;赵军平;周余庆
分类号 G01L3/26(2006.01)I 主分类号 G01L3/26(2006.01)I
代理机构 北京中北知识产权代理有限公司 11253 代理人 程春生
主权项 1.一种基于反求法的变频供水系统运行效率检测方法,其特征在于,包括以下步骤:(1)以采样周期T<sub>s</sub>为间隔对供水系统管网的水压值、逆变电路输出频率和输入功率进行采样,将第一次采样值标记为p(1)、f(1)和P<sub>in</sub>(1);当前采样次数为k,令k=1;(2)建立由M个元素构成的水压值数组{p(i)}、逆变电路输出频率数组{f(i)},以及输入功率数组{P<sub>in</sub>(i)};其中i={k-M+1,k-M+2,...k},M为预先设定的大于1的正整数,k为当前采样次数;p(i)|<sub>i<=0</sub>=0,f(i)|<sub>i<=0</sub>=0,P<sub>in</sub>(i)|<sub>i<=0</sub>=0;(3)判断变频供水系统是否处于稳定状态;如果是,进入步骤(4);否则,变频供水系统处于不稳定状态,转入步骤(16);(4)求解水压平均值<img file="FDA0000379811510000011.GIF" wi="392" he="141" />逆变电路输出频率平均值<maths num="0001"><![CDATA[<math><mrow><mover><mi>F</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>k</mi><mo>-</mo><mi>M</mi><mo>+</mo><mn>1</mn></mrow><mi>k</mi></munderover><mi>f</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></math>]]></maths>和输入功率平均值<maths num="0002"><![CDATA[<math><mrow><mover><msub><mi>P</mi><mi>in</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>M</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mi>k</mi><mo>-</mo><mi>M</mi><mo>+</mo><mn>1</mn></mrow><mi>k</mi></munderover><msub><mi>P</mi><mi>in</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>(5)令n=1,<img file="FDA0000379811510000014.GIF" wi="218" he="140" />T<sub>d</sub>为预先定义的观测时间长度;(6)标记当前时刻为t=0时刻,给逆变电路输出频率一个固定的任意扰动ΔF<sub>n</sub>;(7)令m=1;(8)判断mT<sub>s</sub>>T<sub>d</sub>是否成立,如果成立,则转入步骤(11);否则,在t=mT<sub>s</sub>时刻,采样管网压力值p<sub>n</sub>(m),计算得到<img file="FDA0000379811510000015.GIF" wi="469" he="92" />(9)判断<img file="FDA0000379811510000021.GIF" wi="354" he="163" />是否成立,如果不成立,标记N<sup>n</sup>=m,转入步骤(12);否则,将Δp<sub>n</sub>(m)、<img file="FDA0000379811510000022.GIF" wi="190" he="83" />ΔF<sub>n</sub>、T<sub>b</sub>、ρ、g、P<sub>b</sub>、V<sub>b</sub>、T和t=mT<sub>s</sub>代入公式<maths num="0003"><![CDATA[<math><mrow><mfrac><mrow><mi>&Delta;</mi><msub><mi>p</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow></mrow><mover><mi>P</mi><mo>&OverBar;</mo></mover></mfrac><mo>=</mo><mfrac><mrow><mo>(</mo><mn>2</mn><mover><mi>F</mi><mo>&OverBar;</mo></mover><mo>&times;</mo><mi>&Delta;</mi><msub><mi>F</mi><mi>n</mi></msub><mo>&times;</mo><mi>&Delta;</mi><msup><msub><mi>F</mi><mi>n</mi></msub><mn>2</mn></msup><mo>)</mo></mrow><msup><mover><mi>F</mi><mo>&OverBar;</mo></mover><mn>2</mn></msup></mfrac><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mrow><msubsup><mi>P</mi><mi>out</mi><mi>n</mi></msubsup><mo>[</mo><mi>m</mi><mo>]</mo><msub><mi>T</mi><mi>b</mi></msub></mrow><mrow><mi>&rho;g</mi><msub><mi>P</mi><mi>b</mi></msub><msub><mi>V</mi><mi>b</mi></msub><mi>T</mi></mrow></mfrac><mi>t</mi></mrow></msup><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>求解得出<img file="FDA0000379811510000024.GIF" wi="185" he="75" />其中,P<sub>b</sub>为供水系统气压罐额定压力值,V<sub>b</sub>为供水系统气压罐气室额定体积,T<sub>b</sub>为供水系统气压罐额定温度;T为环境温度,ρ为液体密度;g为重力加速度;(10)更新变量,令m=m+1;返回步骤(8);(11)计算平均值<img file="FDA0000379811510000025.GIF" wi="386" he="149" />转入步骤(13);(12)计算平均值<img file="FDA0000379811510000026.GIF" wi="403" he="149" />进入步骤(13);(13)更新变量,令n=n+1;判断n>5是否成立,如果是,进入步骤(14);否则,转入步骤(6);(14)计算标准差<img file="FDA0000379811510000027.GIF" wi="632" he="229" />判断σ<sub>Pout</sub><0.5是否成立;如果是,则进入步骤(15);否则,转入步骤(16)。(15)令<img file="FDA0000379811510000028.GIF" wi="315" he="141" />计算系统效率<img file="FDA0000379811510000029.GIF" wi="350" he="169" />判断η<=η<sub>min</sub>是否满足;其中η<sub>min</sub>为预先设定的低效率阈值;如果是,则说明系统处于低效率工作状态,逆变器输出关闭,退出;否则,进入步骤(16);(16)令k=k+1;在本次采样周期结束后,进行下一次采样,并标记水压值、逆变电路输出频率和输入功率的采样值为p(k)、f(k)和P<sub>in</sub>(k);返回步骤(2)。
地址 325000 浙江省温州市茶山高教园区温州大学机电工程学院