发明名称 一种矿井通风系统的联合优调优控方法
摘要 本发明提供一种矿井通风系统的联合优调优控方法,建立风机-系统联合调节非线性规划模型,根据矿井通风系统拓扑结构、根据指定的需风风道的需风要求和系统可调参数,适时计算出通风系统的最优调节方案,即:在最大限度地满足需风要求、并使总功率消耗最小、各种控制参数可实施的目标下,计算出风机控制参数和调节风道的调节参数;并发布调节和控制命令完成矿井通风系统的在线最优调节,从而实现风机和风网的联合调节,保证矿井通风系统的按时按需低功耗最优供风。
申请公布号 CN104948215A 申请公布日期 2015.09.30
申请号 CN201510256440.8 申请日期 2015.05.19
申请人 卢新明;尹红 发明人 卢新明;尹红
分类号 E21F1/00(2006.01)I;E21F1/04(2006.01)I;E21F1/10(2006.01)I;F04D27/00(2006.01)I 主分类号 E21F1/00(2006.01)I
代理机构 北京市盛峰律师事务所 11337 代理人 席小东
主权项 一种矿井通风系统的联合优调优控方法,其特征在于,包括以下步骤:步骤1,根据矿井通风系统调控目的,将风道划分为四类,分别为Ⅰ类风道、Ⅱ类风道、Ⅲ类风道和Ⅳ类风道;其中,Ⅰ类风道、Ⅱ类风道和Ⅲ类风道为控制风道;Ⅳ类风道为非控制风道;具体的,Ⅰ类风道为风机所在风道;Ⅱ类风道为需风风道;Ⅲ类风道为安装有风阻调节装置的调节风道;然后,建立完整的矿井通风系统网络拓扑结构图并编号,形成用于优化调控计算的通风系统网络图;步骤2,设所述通风系统共有n条风道,m个独立回路;则对所述通风系统网络图的n条风道由1开始进行编号,由此得到m行n列的基本回路矩阵<img file="FDA0000720127370000016.GIF" wi="247" he="88" />其中,b<sub>kj</sub>为独立回路矩阵系数;此处,令E代表所有风道的集合;令I代表独立回路对应的风道集;令E<sub>1</sub>代表Ⅰ类风道的集合;令E<sub>2</sub>代表Ⅱ类风道的集合;令E<sub>3</sub>代表Ⅲ类风道的集合;步骤3,对于Ⅲ类风道集合E<sub>3</sub>中的任意一个调节风道,记为调节风道j,将调节风道j的总风阻R<sub>j</sub>分解为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>R</mi><mi>j</mi></msub><mo>=</mo><msubsup><mi>R</mi><mi>j</mi><mrow><mi>F</mi><mn>0</mn></mrow></msubsup><mo>+</mo><msubsup><mi>R</mi><mi>j</mi><mrow><mi>T</mi><mn>0</mn></mrow></msubsup><mo>+</mo><msubsup><mi>R</mi><mi>j</mi><mi>T</mi></msubsup></mrow>]]></math><img file="FDA0000720127370000011.GIF" wi="430" he="89" /></maths>其中:<img file="FDA0000720127370000012.GIF" wi="94" he="92" />—调节风道j的基础风阻,为已知值;<img file="FDA0000720127370000013.GIF" wi="87" he="84" />—调节风道j的风阻调节装置不可开启部分产生的局部固定风阻,为已知值;<img file="FDA0000720127370000014.GIF" wi="69" he="84" />—调节风道j的风阻调节装置可开启部分产生的风阻增量,为控制参数;步骤4,对每一台可调节的风机i,建立以下三个风机多参数性能曲线:1)风压‑(叶片角度,转速,风量)曲线:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>h</mi><mi>i</mi><mi>f</mi></msubsup><mo>=</mo><msubsup><mi>h</mi><mi>i</mi><mi>f</mi></msubsup><mrow><mo>(</mo><msub><mi>n</mi><mi>i</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>,</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>A</mi><mi>h</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>)</mo></mrow><msubsup><mi>q</mi><mi>i</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>n</mi><mi>i</mi></msub><msub><mi>B</mi><mi>h</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mi>q</mi><mi>i</mi></msub><mo>+</mo><msubsup><mi>n</mi><mi>i</mi><mn>2</mn></msubsup><msub><mi>C</mi><mi>h</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000720127370000015.GIF" wi="1083" he="85" /></maths>其中:<img file="FDA0000720127370000021.GIF" wi="72" he="78" />—风机i的风压;n<sub>i</sub>—风机i的转速;q<sub>i</sub>—风机i的风量;θ<sub>i</sub>—风机i的动叶角度;A<sub>h</sub>,B<sub>h</sub>和C<sub>h</sub>—风机i的转速为1时,风机i的风压‑风量特性曲线系数;2)效率‑(叶片角度,转速,风量)曲线<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>&eta;</mi><mi>i</mi><mi>f</mi></msubsup><mo>=</mo><msubsup><mi>&eta;</mi><mi>i</mi><mi>f</mi></msubsup><mrow><mo>(</mo><msub><mi>n</mi><mi>i</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>,</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>n</mi><mi>i</mi><mrow><mo>-</mo><mn>2</mn></mrow></msubsup><msub><mi>A</mi><mi>&eta;</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>)</mo></mrow><msubsup><mi>q</mi><mi>i</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>n</mi><mi>i</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><msub><mi>B</mi><mi>&eta;</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>)</mo></mrow><msub><mi>q</mi><mi>i</mi></msub><mo>+</mo><msub><mi>C</mi><mi>&eta;</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000720127370000022.GIF" wi="1132" he="91" /></maths>其中,<img file="FDA0000720127370000023.GIF" wi="66" he="78" />—风机i的效率3)功率‑(叶片角度,转速,风量)曲线<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>P</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>n</mi><mi>i</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>,</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>C</mi><mi>i</mi><mi>f</mi></msubsup><msub><mi>q</mi><mi>i</mi></msub><msubsup><mi>h</mi><mi>i</mi><mi>f</mi></msubsup><mrow><mo>(</mo><msub><mi>n</mi><mi>i</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>,</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>/</mo><msubsup><mi>&eta;</mi><mi>i</mi><mi>f</mi></msubsup><mrow><mo>(</mo><msub><mi>n</mi><mi>i</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>,</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000720127370000024.GIF" wi="920" he="85" /></maths>其中,P<sub>i</sub>(n<sub>i</sub>,θ<sub>i</sub>,q<sub>i</sub>)—风机i的功率;<img file="FDA0000720127370000025.GIF" wi="73" he="85" />—风机i的功率系数,与量刚有关;对于上述三个风机多参数性能曲线:当安装的风机为变频调速风机,即:仅具有调速功能、不具有动叶角度调节功能时,风机i的动叶角度θ<sub>i</sub>为固定常数值,转速n<sub>i</sub>为变量;当安装的风机为仅具有动叶角度调节功能、不具有调速功能时,风机i的转速n<sub>i</sub>为固定常数值、风机i的动叶角度θ<sub>i</sub>为变量;当安装的风机同时具有动叶角度调节功能和调速功能时,风机i的动叶角度θ<sub>i</sub>和风机i的转速n<sub>i</sub>均为变量;步骤5,在步骤1、2、3、4的基础上,建立如下的风机‑系统联合调节非线性规划模型:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>min</mi><mi>F</mi><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>i</mi><mo>&Element;</mo><msub><mi>E</mi><mn>1</mn></msub></mrow></munder><msubsup><mi>C</mi><mi>i</mi><mi>f</mi></msubsup><msub><mi>q</mi><mi>i</mi></msub><msubsup><mi>h</mi><mi>i</mi><mi>f</mi></msubsup><mrow><mo>(</mo><msub><mi>n</mi><mi>i</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>,</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>/</mo><msubsup><mi>&eta;</mi><mi>i</mi><mi>f</mi></msubsup><mrow><mo>(</mo><msub><mi>n</mi><mi>i</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>,</mo><msub><mi>q</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><mi>M</mi><munder><mi>&Sigma;</mi><mrow><mi>i</mi><mo>&Element;</mo><msub><mi>E</mi><mn>2</mn></msub></mrow></munder><msup><mrow><mo>(</mo><msub><mi>q</mi><mi>i</mi></msub><mo>-</mo><msubsup><mi>q</mi><mi>i</mi><mn>0</mn></msubsup><mo>)</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000720127370000026.GIF" wi="1223" he="132" /></maths>约束条件为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><msub><mi>E</mi><mn>3</mn></msub></mrow></munder><msub><mi>b</mi><mi>kj</mi></msub><mo>[</mo><mrow><mo>(</mo><msubsup><mi>R</mi><mi>j</mi><mrow><mi>F</mi><mn>0</mn></mrow></msubsup><mo>+</mo><msubsup><mi>R</mi><mi>j</mi><mrow><mi>T</mi><mn>0</mn></mrow></msubsup><mo>+</mo><msub><mi>y</mi><mi>j</mi></msub><mo>)</mo></mrow><msup><mrow><mo>|</mo><msub><mi>q</mi><mi>j</mi></msub><mo>|</mo></mrow><mrow><msub><mi>&sigma;</mi><mi>j</mi></msub><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>q</mi><mi>j</mi></msub><mo>-</mo><msubsup><mi>h</mi><mi>j</mi><mi>f</mi></msubsup><mrow><mo>(</mo><msub><mi>n</mi><mi>j</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>j</mi></msub><mo>,</mo><msub><mi>q</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>h</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo><mo>+</mo></mtd></mtr><mtr><mtd><munder><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&Element;</mo><mi>E</mi><mo>-</mo><msub><mi>E</mi><mn>3</mn></msub></mrow></munder><msub><mi>b</mi><mi>kj</mi></msub><mo>[</mo><msubsup><mi>R</mi><mi>j</mi><mrow><mi>F</mi><mn>0</mn></mrow></msubsup><mo>+</mo><msubsup><mi>R</mi><mi>j</mi><mrow><mi>T</mi><mn>0</mn></mrow></msubsup><msup><mrow><mo>|</mo><msub><mi>q</mi><mi>j</mi></msub><mo>|</mo></mrow><mrow><msub><mi>&sigma;</mi><mi>j</mi></msub><mo>-</mo><mn>1</mn></mrow></msup><msub><mi>q</mi><mi>j</mi></msub><mo>-</mo><msubsup><mi>h</mi><mi>j</mi><mi>f</mi></msubsup><mrow><mo>(</mo><msub><mi>n</mi><mi>j</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>j</mi></msub><mo>,</mo><msub><mi>q</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>h</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow><mo>]</mo><mo>=</mo><msub><mi>&gamma;</mi><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>&Element;</mo><mi>I</mi></mtd></mtr><mtr><mtd><msub><mi>q</mi><mi>j</mi></msub><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>k</mi><mo>&Element;</mo><mi>I</mi></mrow></munder><msub><mi>b</mi><mi>kj</mi></msub><msub><mi>q</mi><mi>k</mi></msub><mo>+</mo><msubsup><mi>q</mi><mi>j</mi><mn>0</mn></msubsup><mo>,</mo><mi>j</mi><mo>&Element;</mo><mi>E</mi><mo>-</mo><mi>I</mi></mtd></mtr><mtr><mtd><mn>0</mn><mo>&le;</mo><msub><mi>y</mi><mi>j</mi></msub><mo>&le;</mo><msup><msub><mi>y</mi><mi>j</mi></msub><mn>1</mn></msup><mo>,</mo><mi>j</mi><mo>&Element;</mo><msub><mi>E</mi><mn>3</mn></msub></mtd></mtr><mtr><mtd><msubsup><mi>&theta;</mi><mi>j</mi><mn>0</mn></msubsup><mo>&le;</mo><msub><mi>&theta;</mi><mi>j</mi></msub><mo>&le;</mo><msubsup><mi>&theta;</mi><mi>j</mi><mn>1</mn></msubsup><mo>,</mo><mi>j</mi><mo>&Element;</mo><msub><mi>E</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msup><msub><mi>n</mi><mi>j</mi></msub><mn>0</mn></msup><mo>&le;</mo><msub><mi>n</mi><mi>j</mi></msub><mo>&le;</mo><msup><msub><mi>n</mi><mi>j</mi></msub><mn>1</mn></msup><mo>,</mo><mi>j</mi><mo>&Element;</mo><msub><mi>E</mi><mn>1</mn></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000720127370000031.GIF" wi="1141" he="702" /></maths>其中:F—目标函数;以控制时期整个通风系统动力消耗最小、各需风风道的需风量与调节风量差平方和最小为目标;M—惩罚因子;<img file="FDA0000720127370000032.GIF" wi="58" he="78" />—风道i的需风量;γ<sub>k</sub>—由基本回路的风压平衡方程推导出来的右侧常数项,通常为0;b<sub>kj</sub>—独立回路矩阵系数;σ<sub>j</sub>—风道j中风流的流态指数;h<sub>t</sub>(j)—风道j的自然风压;y<sub>j</sub>—风道j的风阻调节增量,与<img file="FDA0000720127370000033.GIF" wi="63" he="85" />含义相同,此处代表未知控制量;y<sub>j</sub><sup>1</sup>—风道j的风阻调节上限;<img file="FDA0000720127370000034.GIF" wi="58" he="85" />和<img file="FDA0000720127370000035.GIF" wi="51" he="83" />—风机j的动叶角度的调节下限和上限;n<sub>j</sub><sup>0</sup>和n<sub>j</sub><sup>1</sup>—风机j转速调节的下限和上限;在上述风机‑系统联合调节非线性规划模型中,各调节风道的风阻调节增量y为控制参数;调节风机的动叶角度θ和/或调节风机转速n为控制参数,此处,如果风机为变频调速风机,则仅转速n为控制参数;如果风机为动叶角度调节风机,则仅动叶角度θ为控制参数;如果风机同时具有动叶角度调节功能和调速功能时,则转速n和动叶角度θ均为控制参数;除上述控制参数之外,其他所有参数均为已知固定值,通过求解该非线性规划模型,计算出通风系统的最优调节方案,即:在最大限度地满足需风要求、并使总功率消耗最小、各种控制参数可实施的情况下,计算得到各调节风道的控制参数,包括:风阻调节增量y,还包括,调节风机的动叶角度θ和/或调节风机转速n。
地址 271000 山东省泰安市高新技术开发区高创中心
您可能感兴趣的专利