发明名称 一种硫化过程嵌入式碳排放监控与检测系统
摘要 一种硫化过程嵌入式碳排放监控与检测系统,包括能耗传感器、能耗采集单元、碳排放监控单元和嵌入式碳排放处理单元;所述能耗采集单元通过网络获得所述能耗传感器采集的能耗数据,并将所述能耗数据传送给所述碳排放监控单元和嵌入式碳排放处理单元进行处理;所述碳排放监控单元根据异常检测模型进行异常告警处理;所述嵌入式碳排放处理单元包括碳排放获得单元、碳排放优化识别单元和检测单元;所述碳排放获得单元包括修正单元和处理单元,本发明在碳排放监控与检测系统架构等方面有较大突破,同时对于提高企业节能管理水平,加大节能技术改造,减轻环境污染,缓解能源瓶颈制约,实现的节约发展、清洁发展和可持续发展具有十分重要的战略意义和现实意义。
申请公布号 CN103425099A 申请公布日期 2013.12.04
申请号 CN201310299342.3 申请日期 2013.07.16
申请人 广东工业大学 发明人 杨海东;杨春
分类号 G05B19/418(2006.01)I 主分类号 G05B19/418(2006.01)I
代理机构 北京品源专利代理有限公司 11332 代理人 马晓亚
主权项 1.一种硫化过程嵌入式碳排放监控与检测系统,其特征在于:包括能耗传感器、能耗采集单元、碳排放监控单元和嵌入式碳排放中央处理单元;所述能耗传感器包括供电系统电表、空调系统电表、锅炉系统温度和压力表、供水系统水表、空压系统压力表和制氮系统压力表;所述能耗采集单元通过网络获得所述能耗传感器采集的能耗数据,并将所述能耗数据传送给所述碳排放监控单元和嵌入式碳排放中央处理单元进行处理;所述碳排放监控单元以发现能耗异常为目标,将所述能耗数据进行监控和分析,根据异常检测模型进行异常告警处理,或者直接对可以控制的能耗设备进行启停处理;所述嵌入式碳排放中央处理单元包括碳排放获得单元、碳排放优化识别单元和检测单元;所述碳排放获得单元包括修正单元和处理单元;所述碳排放优化识别单元设置有分别保存碳排放识别优化程序和碳排放识别优化数据的存储器,所述存储器采用各自的总线、拥有独立的地址,可以独立访问;所述碳排放识别优化程序包括划分单元、提取单元和识别单元,所述划分单元用于将分布式生产过程中碳排放分为直接碳排放和间接碳排放;所述直接碳排放为制造单元碳排放,包括生产过程中机器耗费的电能、水和蒸汽;所述间接碳排放为仪表盘、照明系统和冷却系统碳排放;所述提取单元将实时采集的能源数据按照分钟、小时、日、月和年的时间和表计、生产线、车间和区域的空间进行汇总计算,形成分析所需的能耗用量、标准煤、费用和碳排放的能耗数据,获取生产过程中的排班、设备或车间产出、质量检测等数据,进行单产能耗、单位面积能耗、设备能效、车间能效能耗指标的计算;所述识别单元根据能效优化运行的需要,对各能耗数据进行分析处理,包括能耗平衡分析、制造能效分析和节能效果分析;所述检测单元包括电子键盘和显示触摸屏;所述电子键盘用于实现对碳排放进行分布式、全面的监控,通过以太网或者CAN总线与远程信息系统进行远程数据和命令通信,根据监控中心发出的不同指令,对设备及时执行正确的操作;所述显示触摸屏根据碳排放管理需要,进行选择性的显示,监测人员通过所述显示触摸屏实时查看各项碳排放指标;所述异常检测模型使用面向流程的数据集成模型和数据回归模型来集成碳排放的不同的因素,用于碳排放异常分析;所述数据集成模型按工序为单位进行采集和计量,获得表征能源效率的能流率EFR<sub>vul</sub>,所述能流率EFR<sub>vul</sub>为理论能耗TE和直接能耗DE的比;DE是制造过程中的实际能源消耗;TE是制造过程需要的最小能源:<img file="FDA00003519356400021.GIF" wi="324" he="150" />其中0≤EFR<sub>vul</sub>≤1;TE是执行制造过程需要的最小能源消耗,其包括物理过程的能源消耗TE<sub>p</sub>和化学过程的能源消耗TE<sub>c</sub>,即其表达式如下:TE=TE<sub>p</sub>+TE<sub>c</sub>;在硫化过程中,气体物理能量转换的第一阶段常温条件下的气体压缩所消耗的能量大小如下表达式计算:<maths num="0001"><![CDATA[<math><mrow><mi>W</mi><mo>=</mo><munderover><mo>&Integral;</mo><msub><mi>v</mi><mn>2</mn></msub><msub><mi>v</mi><mn>1</mn></msub></munderover><mi>P</mi><mo>&CenterDot;</mo><mi>dV</mi><mo>=</mo><munderover><mo>&Integral;</mo><msub><mi>v</mi><mn>2</mn></msub><msub><mi>v</mi><mn>1</mn></msub></munderover><mfrac><mi>nRT</mi><mi>V</mi></mfrac><mo>&CenterDot;</mo><mi>dV</mi><mo>=</mo><mi>nRT</mi><mi>ln</mi><mfrac><msup><mi>p</mi><mo>*</mo></msup><msub><mi>p</mi><mn>1</mn></msub></mfrac><mo>;</mo></mrow></math>]]></maths>在第二阶段用于加热材料的能源使用量计算如下:Q=c·m·(t<sub>2</sub>-t<sub>1</sub>);上式中的c是气体比热容,m是气体的质量;物理能耗大小为:<maths num="0002"><![CDATA[<math><mrow><msub><mi>TE</mi><mi>p</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>{</mo><msub><mi>W</mi><mi>i</mi></msub><mo>+</mo><msub><mi>Q</mi><mi>i</mi></msub><mo>}</mo><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>{</mo><msub><mi>nRT</mi><mi>il</mi></msub><mi>ln</mi><mfrac><msup><msub><mi>p</mi><mi>i</mi></msub><mo>*</mo></msup><msub><mi>p</mi><mi>il</mi></msub></mfrac><mo>+</mo><msub><mi>c</mi><mi>i</mi></msub><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>t</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>}</mo><mo>;</mo></mrow></math>]]></maths>化学反应过程中TE<sub>c</sub>的计算表达式如下:TE<sub>c</sub>=m·q上式中的q是单位物质吸热/放热量,m参加反应的橡胶质量,取q=-43×10<sup>3</sup>J·kg<sup>-1</sup>;面向设备DE的计算:DE通过面向设备公式进行计算,其能够对所述实际能源消耗进行测量,包括主要设备的能源消耗PE和辅助设备的能源消耗AE;其表达式如下:DE=PE+AEPE是主要设备使用的电能和煤炭量,其包括硫化、压缩、氮气的制造和水蒸气的制造;AE是电力能源的消耗,其包括照明、通风排气和空调;<maths num="0003"><![CDATA[<math><mrow><mi>DE</mi><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>m</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>Q</mi><mi>i</mi></msub></mrow></math>]]></maths>上式中的m<sub>i</sub>(i=1,2)是电能和煤炭消耗的质量,Q<sub>i</sub>(i=1,2)分别是他们能量转化率;所述数据回归模型是通过利用统计模型和人工神经网络为所述能源效率建立适当的回归模型,在预测区间内能源消耗的实际值来判断碳排放异常,使ER<sub>vul</sub>可以预先预测;硫化过程的能源效率的回归模型是y<sub>t</sub>=f(X<sub>t</sub>),其中y<sub>t</sub>是ER<sub>vul</sub>当前的预测值,f(·)是回归函数,X<sub>t</sub>是影响硫化能效的一个矩阵,它在硫化过程初始时获得,包括操作者爱好、环境、轮胎坯料和硫化机器;f(·)从历史数据库中提取训练样本;根据无偏性标准和最小方差标准,绝对平均误差M<sub>e</sub>和绝对误差S<sub>e</sub>的标准差用于估算预测精度,其定义为如下式:<maths num="0004"><![CDATA[<math><mrow><msub><mi>M</mi><mi>e</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>|</mo><msub><mover><mi>y</mi><mo>~</mo></mover><mi>t</mi></msub><mo>-</mo><msub><mi>y</mi><mi>t</mi></msub><mo>|</mo><mo>;</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>S</mi><mi>e</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><mo>|</mo><msub><mi>e</mi><mi>t</mi></msub><mo>|</mo><mo>-</mo><msub><mi>M</mi><mi>e</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt><mo>;</mo></mrow></math>]]></maths>其中<img file="FDA00003519356400041.GIF" wi="57" he="71" />是预测值,y<sub>t</sub>是ER<sub>vul</sub>观察值,<img file="FDA00003519356400042.GIF" wi="262" he="73" />n是训练样本数;最小二乘标准,MinΣe<sub>t</sub><sup>2</sup>作为目标函数;函数f(·)只给出ER<sub>vul</sub>预测点,诊断异常事件需要一个置信区间,在此区间ER<sub>vul</sub>从正常处理过程产生预计会下降,如果观察值不在置信区间所决定的范围内,那么能耗预警将会触发;对于第(n+1)处理阶段,ER<sub>vul</sub>的预测值由<img file="FDA00003519356400044.GIF" wi="90" he="71" />表示,而观测值则用y<sub>n+1</sub>表示,假设<img file="FDA00003519356400045.GIF" wi="96" he="74" />和y<sub>n+1</sub>与正态随机分布相同,<img file="FDA00003519356400046.GIF" wi="354" he="74" />而e<sub>n+1</sub>同样也是正态分布,其期望E(e<sub>n+1</sub>)=0或e<sub>n+1</sub>~N(0,δ<sup>2</sup>);标准差的估计表达式:<maths num="0006"><![CDATA[<math><mrow><msup><mover><mi>&delta;</mi><mo>~</mo></mover><mn>2</mn></msup><mo>=</mo><mfrac><mn>1</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><msub><mi>e</mi><mi>t</mi></msub><mn>2</mn></msup><mo>;</mo></mrow></math>]]></maths>上式中n是训练样本数量,建立统计z为如下表达式:<maths num="0007"><![CDATA[<math><mrow><mi>z</mi><mo>=</mo><mfrac><mrow><msub><mi>e</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><mi>E</mi><mrow><mo>(</mo><msub><mi>e</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow><mrow><mi>D</mi><mrow><mo>(</mo><msub><mi>e</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>y</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mover><mi>y</mi><mo>~</mo></mover><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow><msup><mover><mi>&delta;</mi><mo>~</mo></mover><mn>2</mn></msup></mfrac><mo>~</mo><mi>N</mi><mrow><mo>(</mo><mn>0,1</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>给定一个显著水平α,置信水平为1-α,则y<sub>n+1</sub>的置信区间为<maths num="0008"><![CDATA[<math><mrow><msub><mover><mi>y</mi><mo>~</mo></mover><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>Z</mi><mrow><mi>&alpha;</mi><mo>/</mo><mn>2</mn></mrow></msub><mover><mi>&delta;</mi><mo>~</mo></mover><mo>&le;</mo><msub><mi>y</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>&le;</mo><msub><mover><mi>y</mi><mo>~</mo></mover><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>Z</mi><mrow><mi>&alpha;</mi><mo>/</mo><mn>2</mn></mrow></msub><mover><mi>&delta;</mi><mo>~</mo></mover><mo>.</mo></mrow></math>]]></maths>
地址 510006 广东省广州市番禺区大学城外环西路100号