发明名称 一种选矿工业数质量平衡计算方法
摘要 本发明涉及一种选矿工业数质量平衡计算方法,其步骤:1)依据实际选矿工业流程总输入质量与总输出质量相等列出质量平衡方程;2)依据先验知识作为约束条件,将质量平衡方程转化为非线性优化问题;3)求解步骤2)中的非线性优化方程,得到的可测量<img file="DDA0000549549080000011.GIF" wi="59" he="73" />作为可测量的修正值,得到的所求量<img file="DDA0000549549080000012.GIF" wi="56" he="75" />作为所求量的推算值,进而解决了数质量平衡问题。本发明由于充分利用先验知识,通过先验知识将单纯的质量平衡方程改写为非线性优化问题进行求解,从而对现场测量的误差有非常大的容忍性,即使在测量存在较大误差的情况下,也能够自动修复不合理的测量值,推算出合理可信的结果,可广泛在选矿工业领域中应用。
申请公布号 CN104123475A 申请公布日期 2014.10.29
申请号 CN201410383948.X 申请日期 2014.08.06
申请人 清华大学 发明人 王焕钢;王正;徐文立;周俊武;徐宁;王庆凯
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 北京纪凯知识产权代理有限公司 11245 代理人 徐宁;孙楠
主权项 一种选矿工业数质量平衡计算方法,其包括以下步骤:1)依据实际选矿工业流程总输入质量与总输出质量相等列出质量平衡方程<img file="FDA0000549549050000011.GIF" wi="241" he="70" />对第k个子过程,有:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>g</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>,</mo><msub><mover><mi>y</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mn>0</mn><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi></mrow>]]></math><img file="FDA0000549549050000012.GIF" wi="551" he="70" /></maths>式中,所有的下标k均代表该质量平衡方程针对的是第k个子过程;<img file="FDA0000549549050000013.GIF" wi="60" he="71" />代表通过能够测量得到的数质量,简称为可测量;<img file="FDA0000549549050000014.GIF" wi="61" he="70" />代表不能测量但希望知道的数质量,简称为所求量;n代表工业流程中子过程的总个数,同时也是质量平衡方程的个数;2)依据先验知识作为约束条件,将质量平衡方程转化为非线性优化问题:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><munder><mi>min</mi><mrow><msub><mover><mi>x</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>,</mo><msub><mover><mi>y</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi></mrow></munder><mi>f</mi><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>,</mo><msub><mover><mi>y</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&alpha;</mi><mi>k</mi></msub><msup><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>-</mo><msub><mover><mi>X</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>&beta;</mi><mi>k</mi></msub><msub><mi>g</mi><mi>k</mi></msub><msup><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>,</mo><msub><mover><mi>y</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000549549050000015.GIF" wi="1394" he="138" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>p</mi><mo>.</mo><mi>t</mi><mo>.</mo><mi></mi><msub><mi>p</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>,</mo><msub><mover><mi>y</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>m</mi><mi>p</mi></msub></mrow>]]></math><img file="FDA0000549549050000016.GIF" wi="923" he="70" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>q</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>,</mo><msub><mover><mi>y</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>&lt;</mo><mn>0</mn><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>m</mi><mi>q</mi></msub></mrow>]]></math><img file="FDA0000549549050000017.GIF" wi="815" he="80" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mover><mi>a</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>&le;</mo><msub><mover><mi>x</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>&le;</mo><msub><mover><mi>b</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi></mrow>]]></math><img file="FDA0000549549050000018.GIF" wi="561" he="84" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mover><mi>c</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>&le;</mo><msub><mover><mi>y</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>&le;</mo><msub><mover><mi>d</mi><mo>&RightArrow;</mo></mover><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi></mrow>]]></math><img file="FDA0000549549050000019.GIF" wi="561" he="87" /></maths>其中,<img file="FDA00005495490500000110.GIF" wi="412" he="65" />代表优化的目标函数,其自变量是所有可测量<img file="FDA00005495490500000111.GIF" wi="54" he="70" />和所有所求量<img file="FDA00005495490500000112.GIF" wi="78" he="71" />n代表质量平衡方程的个数;α<sub>k</sub>、β<sub>k</sub>代表系数;<img file="FDA00005495490500000113.GIF" wi="69" he="82" />表示可测量<img file="FDA00005495490500000114.GIF" wi="46" he="71" />的实测值;<img file="FDA00005495490500000115.GIF" wi="440" he="74" />代表等式形式的先验知识;<img file="FDA00005495490500000116.GIF" wi="436" he="74" />代表不等式形式的先验知识;m<sub>p</sub>代表等式形式先验知识的个数;m<sub>q</sub>代表不等式形式先验知识的个数;<img file="FDA00005495490500000117.GIF" wi="63" he="74" />代表可测量的经验下限;<img file="FDA00005495490500000118.GIF" wi="58" he="87" />代表可测量的经验上限;<img file="FDA00005495490500000119.GIF" wi="54" he="74" />代表所求量的经验下限;<img file="FDA00005495490500000120.GIF" wi="62" he="83" />代表所求量的经验上限;3)求解步骤2)中的非线性优化方程,得到的可测量<img file="FDA00005495490500000122.GIF" wi="58" he="68" />作为可测量的修正值,得到的所求量<img file="FDA00005495490500000121.GIF" wi="58" he="70" />作为所求量的推算值,进而解决了数质量平衡问题。
地址 100084 北京市海淀区100084信箱82分箱清华大学专利办公室