发明名称 一种碳酸化分解过程原液流量的计量模型
摘要 一种碳酸化分解过程原液流量的计量模型,本发明将碳酸化过程基于汇流的伯努利方程,推导出高位槽液位、原液控制阀门与原液流量的函数关系,同时设计了一个结合一阶动态响应的滤波器对计算结果进行在线修正。本发明提供的原液流量的计量模型,精确度高,能够对流量进行在线测量和修正,为生产氧化铝的连续碳酸化分解控制原液流量提供了依据。
申请公布号 CN101339057A 申请公布日期 2009.01.07
申请号 CN200810032081.8 申请日期 2008.08.15
申请人 中南大学 发明人 桂卫华;胡志坤;阳春华;王晓丽;刘伟;赵东峰
分类号 G01F1/00(2006.01) 主分类号 G01F1/00(2006.01)
代理机构 中南大学专利中心 代理人 胡燕瑜
主权项 1.一种碳酸化分解过程原液流量的计量模型,其特征在于:原液流量Q的计算式为Q=Qa+Qb 其中,Qa=αγπd2ua为第一个高位槽对原液流量的贡献,Qb=αγπd2ub为第二个高位槽对原液流量的贡献,<math><mrow><msub><mi>u</mi><mi>a</mi></msub><mo>=</mo><msqrt><mn>2.02</mn><mrow><mo>(</mo><msub><mi>z</mi><mn>1</mn></msub><mo>-</mo><msub><mi>z</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>+</mo><mn>0.7</mn><msup><msub><mi>u</mi><mi>b</mi></msub><mn>2</mn></msup></msqrt></mrow>表示为第一个高位槽原液的流速,<math><mrow><msub><mi>u</mi><mi>b</mi></msub><mo>=</mo><msqrt><mfrac><mrow><mrow><mo>(</mo><msub><mi>k</mi><mn>1</mn></msub><msub><mi>z</mi><mn>1</mn></msub><mo>+</mo><msub><mi>k</mi><mn>2</mn></msub><msub><mi>z</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><msqrt><msub><mi>k</mi><mn>3</mn></msub><msup><msub><mi>z</mi><mn>1</mn></msub><mn>2</mn></msup><mo>+</mo><msub><mi>k</mi><mn>4</mn></msub><msub><mi>z</mi><mn>1</mn></msub><msub><mi>z</mi><mn>2</mn></msub><mo>+</mo><msub><mi>k</mi><mn>5</mn></msub><msup><msub><mi>z</mi><mn>2</mn></msub><mn>2</mn></msup></msqrt></mrow><msub><mi>k</mi><mn>6</mn></msub></mfrac></msqrt></mrow>表示为第二个高位槽的原液流速,α为阀门系数,当阀门采用螺旋型球形阀时取900,蝶形阀门时取1200,式中,ki,i=1,2...6是中间变量,z1、z2分别是两高位槽的液位,γ为修正因子;中间变量ki,i=1,2...6与损失系数ξc的函数关系是:k1=0.34ξc 2+24.2ξc+279.5;k2=-0.34ξc 2-7.4ξc-189.3k3=-0.4ξc 4+33ξc 3+2459ξc 2+38652ξc+197400k4=-0.9ξc 4+27ξc 3+3480ξc 2-538ξc+294054k5=0.1ξc 4+4.2ξc 3+228ξc 2+6050ξc+66442k6=0.26ξc 2-22ξc-284h为阀门的上升高度,D为阀座的直径,开度<math><mrow><mi>&eta;</mi><mo>=</mo><mfrac><mi>h</mi><mi>D</mi></mfrac><mo>,</mo></mrow>阀门损失系数ξc和开度的函数关系是:当η∈(0,0.25], ξc=-534.8(k-0)+150 (1)当η∈(0.25,0.5],ξc=-24(k-0.25)+16.3 (2)当η∈(0.5,0.75],ξc=-10.48(k-0.5)+10.3 (3)当η∈(0.75,1], ξc=-6.36(k-0.75)+7.68 (4)在计量模型中的修正因子γ,采用低通滤波结合一阶动态响应处理的方法,依靠历史的偏差修正值得到当前偏差的修正值,对检测结果进行在线修正,具体修正方法如下:<math><mrow><msub><mi>E</mi><mi>n</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mo>(</mo><mi>k</mi><mo>=</mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>Q</mi><mi>A</mi><mi>n</mi></msubsup><mo>-</mo><msubsup><mi>Q</mi><mi>E</mi><mi>n</mi></msubsup></mtd><mtd><mrow><mo>(</mo><mi>k</mi><mo>></mo><mn>0</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow><math><mrow><mi>&Delta;</mi><msub><mi>Q</mi><mi>k</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>E</mi><mi>k</mi></msub></mtd><mtd><mrow><mo>(</mo><mi>k</mi><mo>=</mo><mn>0,1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&gamma;</mi><mi>k</mi></msub><mo>)</mo></mrow><mi>&Delta;</mi><msub><mi>Q</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>&gamma;</mi><mi>k</mi></msub><msub><mi>E</mi><mi>k</mi></msub></mtd><mtd><mrow><mo>(</mo><mi>k</mi><mo>></mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow><math><mrow><mi>&Delta;</mi><msub><mi>Q</mi><mi>n</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>E</mi><mi>n</mi></msub></mtd><mtd><mrow><mo>(</mo><mi>k</mi><mo>=</mo><mn>0,1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&gamma;</mi><mi>n</mi></msub><mo>)</mo></mrow><mi>&Delta;</mi><msub><mi>Q</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>&gamma;</mi><mi>n</mi></msub><msub><mi>E</mi><mi>n</mi></msub></mtd><mtd><mrow><mo>(</mo><mi>k</mi><mo>></mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></mrow><math><mrow><mi>&Delta;</mi><msub><mi>E</mi><mi>n</mi></msub><mo>=</mo><mi>&Delta;</mi><msub><mi>Q</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><mrow><mo>(</mo><mi>&Delta;</mi><msub><mi>Q</mi><mi>n</mi></msub><mo>-</mo><mi>&Delta;</mi><msub><mi>Q</mi><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mi>t</mi><mi>&tau;</mi></mfrac></mrow></msup><mo>)</mo></mrow></mrow><math><mrow><msubsup><mi>Q</mi><mi>n</mi><mi>upd</mi></msubsup><mo>=</mo><msubsup><mi>Q</mi><mi>n</mi><mi>cal</mi></msubsup><mo>+</mo><mi>&Delta;</mi><msub><mi>E</mi><mi>n</mi></msub></mrow><math><mrow><msub><mi>&gamma;</mi><mi>n</mi></msub><mo>=</mo><mo>|</mo><mfrac><mrow><msubsup><mi>Q</mi><mi>n</mi><mi>upd</mi></msubsup><mo>-</mo><msubsup><mi>Q</mi><mi>n</mi><mi>cal</mi></msubsup></mrow><msub><mi>E</mi><mi>n</mi></msub></mfrac><mo>|</mo><mo>,</mo><mn>0</mn><mo>&le;</mo><msub><mi>&gamma;</mi><mi>n</mi></msub><mo>&le;</mo><mn>1</mn></mrow>其中,n为修正次数;En为计算值与实际值的偏差;QA n为检测计量模型计算值;QE n为流量测量值;ΔQn为第n次修正的偏差值;γn为修正因子;ΔEn为t时刻的修正偏差;τ为设备的响应时间;Qn upd是t时刻检测计量模型的修正值;Qn cal 是t时刻检测计量模型的计算值。
地址 410083湖南省长沙市麓山南路1号