发明名称 一种钢铁企业蒸汽产生量与消耗量认证的方法
摘要 本发明涉及一种应用于钢铁企业不同工序蒸汽产生量与消耗量认证的方法。本方法对三种不同工序建立蒸汽流量认证模型:对稳定的蒸汽生产源,建立蒸汽流量与锅炉生产变量之间的关系模型;对回收余热蒸汽源的工序,建立影响每台设备蒸汽产生量的因素与蒸汽流量的关系模型;对消耗蒸汽的工序(用户),建立影响各环节蒸汽用量与蒸汽流量的关系模型。本方法规定了建立蒸汽流量认证模型的具体实施过程,并提出了建立模型库,用检索、修正、更新相结合的方法应用蒸汽流量认证模型。本发明提供了钢铁工序蒸汽即时流量与累积流量的实时认证方法,提高了蒸汽计量的可靠性,利于解决钢铁企业蒸汽计量因精度低、停止计量等引起的输差和计量纠纷问题,实现能源考核和节能降耗。
申请公布号 CN103134570B 申请公布日期 2014.11.26
申请号 CN201110388180.1 申请日期 2011.11.29
申请人 中国科学院沈阳自动化研究所 发明人 罗先喜;苑明哲;王宏
分类号 G01F25/00(2006.01)I 主分类号 G01F25/00(2006.01)I
代理机构 沈阳科苑专利商标代理有限公司 21002 代理人 周秀梅;许宗富
主权项 一种钢铁企业蒸汽产生量与消耗量认证的方法,其特征在于,建立蒸汽流量认证模型数据库,将蒸汽流量认证模型按照季节和生产状态存入蒸汽流量认证模型数据库,实施以下步骤:读取来自能源管理系统EMS和制造执行系统MES与蒸汽流量认证模型有关的实时数据;在蒸汽流量认证模型数据库中检索调入符合当前季节和生产状态的蒸汽流量认证模型,所述蒸汽流量认证模型包括稳定蒸汽源产汽量蒸汽流量认证模型、余热回收汽源产汽量蒸汽流量认证模型和消耗蒸汽工序用汽量蒸汽流量认证模型;用制造执行系统MES数据和蒸汽流量认证模型推算蒸汽流量的认证值,所述认证值为蒸汽的即时流量和累积流量;比较蒸汽流量的测量值与认证值,如果相对偏差小于一定范围时,记录测量值为当前流量;如果相对偏差超过一定范围时,判断流量计是否损坏;如果流量计损坏,记录认证值为当前流量,检修仪表;如果流量计没有损坏,修正模型,并将修正模型发生的季节、时间及生产状态信息存入蒸汽流量认证模型数据库;所述稳定蒸汽源产汽量蒸汽流量认证模型的建立方法为:第一步,检查确认用于检测的仪表工作正常;第二步,按照设计确定单台蒸汽锅炉的产汽量的基值F<sub>b0</sub>,及锅炉稳定生产流量为F<sub>b0</sub>时对应的过程变量的值,即初始值向量,记为:x<sub>0</sub>=(F<sub>w0</sub>,h<sub>0</sub>,p<sub>0</sub>,P<sub>0</sub>,T<sub>0</sub>,χ<sub>0</sub>)第三步,在蒸汽锅炉生产状态不同于第二步状态且稳定运行时,记录一组锅炉运行数据(F<sub>b1</sub>,F<sub>w1</sub>,h<sub>1</sub>,p<sub>1</sub>,P<sub>1</sub>,T<sub>1</sub>,χ<sub>1</sub>),每个数据不超过初始值向量对应项10%,形成样本,并将其处理为:x<sub>1</sub>=(F<sub>w1</sub>,h<sub>1</sub>,p<sub>1</sub>,P<sub>1</sub>,T<sub>1</sub>,χ<sub>1</sub>)‑x<sub>0</sub>y<sub>1</sub>=F<sub>b1</sub>‑F<sub>b0</sub>同理类推,得到第2,…,m组数据,(x<sub>2</sub>,y<sub>2</sub>),....(x<sub>m</sub>,y<sub>m</sub>),令X=(x<sub>1</sub>,x<sub>2</sub>...,x<sub>m</sub>)<sup>T</sup>,Y=(y<sub>1</sub>,y<sub>2</sub>...,y<sub>m</sub>)<sup>T</sup>;第四步,利用公式θ=(XX<sup>T</sup>)<sup>‑1</sup>X<sup>T</sup>Y         (7)求蒸汽流量认证模型参数θ;第五步,得到锅炉蒸汽产生量的静态估算模型认证值F<sub>b</sub>(t)=xθ+F<sub>b0</sub>        (8)和锅炉生产蒸汽流量的累积值的认证值<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>Q</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>kT</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><msub><mi>F</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>dt</mi><mo>&ap;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>k</mi></munderover><msub><mi>F</mi><mi>b</mi></msub><mrow><mo>(</mo><mi>k</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><msub><mi>T</mi><mi>s</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000569087150000021.GIF" wi="1315" he="142" /></maths>其中,T<sub>s</sub>为数据采集周期,k为正整数;所述余热回收汽源产汽量蒸汽流量认证模型的建立方法为:第一步,按照工序的特点,将该工序分为N台内部相互独立的产汽设备,确定每台设备基准蒸汽回收率F<sub>ri0</sub>和生产工序设计的单位时间产量;令i=1;第二步,确定第i台设备的蒸汽回收系数q<sub>αi</sub>;当该设备属连续生产设备时q<sub>αi</sub>取1;当该设备属间歇式周期生产设备时,则对一个周期T=aT<sub>s</sub>内蒸汽回收系数按每间隔T<sub>s</sub>采样一次,得到一个序列,则该序列为:<maths num="0002" id="cmaths0002"><img file="FDA0000569087150000022.GIF" wi="1332" he="141" /></maths>式中q<sub>αij</sub>(j=0,1,...,a)为该序列中第j+1次采样值;n<sub>0</sub>T<sub>s</sub>为该序列坐标;δ(n<sub>0</sub>T<sub>s</sub>‑jT<sub>s</sub>)为单位脉冲函数,取多组采样序列对应值的平均值<img file="FDA0000569087150000023.GIF" wi="79" he="76" />作为模型序列的对应值q<sub>αij</sub>;第三步,得到第i台设备在kT<sub>s</sub>产生蒸汽的即时流量的认证值F<sub>i</sub>(kT<sub>s</sub>)=λ<sub>i</sub>(kT<sub>s</sub>)S<sub>1i</sub>(kT<sub>s</sub>)r<sub>i</sub>(kT<sub>s</sub>)q<sub>αi</sub>(kT<sub>s</sub>‑T<sub>0</sub>)F<sub>ri0</sub>        (17)其中,T<sub>s</sub>为数据采集周期,k为正整数,T<sub>0</sub>为该设备在本次生产周期开始的时刻,若为连续生产设备,q<sub>αi</sub>(kT<sub>s</sub>‑T<sub>0</sub>)项取1;第四步,i=i+1转第二步,直至i=N;第五步,得到蒸汽的累积流量的认证值<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>Q</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>kT</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><msub><mi>F</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>dt</mi><mo>&ap;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>k</mi></munderover><msub><mi>F</mi><mi>r</mi></msub><mrow><mo>(</mo><mi>iT</mi><mo>)</mo></mrow><mi>T</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000569087150000031.GIF" wi="1428" he="142" /></maths>第六步,确认仪表正常时,用实际测量值与认证值对照,如相对偏差超过10%时,则返回第一步;否则确认模型;所述消耗蒸汽工序用汽量蒸汽流量认证模型的建立方法为:第一步,将该工序分为M个用汽环节,确定每个用汽环节的基本蒸汽用量Fc0i和该环节设计的单位时间产量;令i=1;第二步,确定第i个环节的蒸汽消耗系数F<sub>αi</sub>:当该生产环节为连续型时F<sub>αi</sub>取1;当该生产环节为间歇式周期生产时,对一个周期T=bT<sub>s</sub>内蒸汽消耗系数按每间隔T<sub>s</sub>采样一次,得到一个序列,则该序列为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>F</mi><mi>&alpha;i</mi></msub><mrow><mo>(</mo><msub><mi>n</mi><mn>0</mn></msub><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>b</mi></munderover><msub><mi>F</mi><mi>&alpha;ij</mi></msub><mi>&delta;</mi><mrow><mo>(</mo><msub><mi>n</mi><mn>0</mn></msub><msub><mi>T</mi><mi>s</mi></msub><mo>-</mo><mi>j</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>22</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000569087150000032.GIF" wi="1336" he="148" /></maths>式中q<sub>αij</sub>(j=0,1,...,b)为该序列中第j+1次采样值;n<sub>0</sub>T<sub>s</sub>为本序列坐标;δ(n<sub>0</sub>T<sub>s</sub>‑jT<sub>s</sub>)为单位脉冲函数,取多组采样序列的对应值的平均值<img file="FDA0000569087150000033.GIF" wi="89" he="76" />作为模型序列的对应值F<sub>αij</sub>;第三步,确定第i环节在kT<sub>s</sub>消耗蒸汽的即时流量认证值F<sub>ci</sub>(kT<sub>s</sub>)=S<sub>2i</sub>(kT<sub>s</sub>)s<sub>i</sub>(kT<sub>s</sub>)F<sub>αi</sub>(kT<sub>s</sub>‑T<sub>1</sub>)F<sub>0i</sub>     (23)式中T<sub>1</sub>为具有周期性生产特征环节该生产环节在本次生产周期开始的时刻,若为连续生产环节,F<sub>αi</sub>(kT<sub>s</sub>‑T<sub>1</sub>)项取1,T<sub>s</sub>为数据采集周期,k为正整数;第四步,i=i+1转第二步,直至i=M;第五步,得到蒸汽的累积流量的认证值<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>Q</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>Q</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>k</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>t</mi></msubsup><msub><mi>F</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>dt</mi><mo>&ap;</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mi>k</mi></munderover><msub><mi>F</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>i</mi><msub><mi>T</mi><mi>s</mi></msub><mo>)</mo></mrow><msub><mi>T</mi><mi>s</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>21</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000569087150000041.GIF" wi="1373" he="142" /></maths>其中,T<sub>s</sub>为数据采集周期,k为正整数;第六步,确认仪表正常时,用实际测量值与认证值对照,如相对偏差超过10%,则返回第一步;否则确认模型;所述与蒸汽流量认证有关的实时数据包括:制造执行系统MES数据:稳定蒸汽源产汽量蒸汽流量认证模型所需的锅炉锅炉加水流量F<sub>w</sub>、锅炉液位h、锅炉汽包的压力p、加热功率P、锅炉蒸汽出口温度T和减温减压加水系数χ;余热回收汽源产汽量蒸汽流量认证模型所需的蒸汽回收设备的运转状态S<sub>1</sub>、生产设备负荷系数r、第i台设备的基准蒸汽回收率F<sub>ri0</sub>、蒸汽回收系数q<sub>αi</sub>和设备附加系数λ;消耗蒸汽工序用汽量蒸汽流量认证模型所需的第i个用汽环节的蒸汽使用状态S<sub>2i</sub>、基本蒸汽用量F<sub>c0i</sub>、季节因子s<sub>i</sub>和蒸汽消耗系数F<sub>αi</sub>;能源管理系统EMS数据:与三种认证模型相对应的蒸汽的即时流量和累积流量的测量值。
地址 110016 辽宁省沈阳市东陵区南塔街114号