发明名称 径向基的农药生产废液焚烧炉炉温最佳化系统及方法
摘要 本发明公开了一种径向基的农药生产废液焚烧炉炉温最佳化系统及方法。该方法使用模糊方程系统对训练样本进行模糊化处理,并采用径向基神经网络作为模糊方程系统的局部方程。在本发明中,信号采集模块依照每次采样的时间间隔,从数据库中采集训练样本数据;训练样本经标准化处理模块处理,作为模糊系统模块的输入,用于软测量建模;模糊系统模块的输出与结果显示模块连接,用于将得到的炉温预报值和使炉温最佳的操作变量值传给DCS系统;模型更新模块,用于按设定的采样时间间隔,采集现场智能仪表信号以更新训练样本集。本发明实现了炉温的实时准确控制、噪声抑制能力强且运算速度快。
申请公布号 CN103488205B 申请公布日期 2015.09.30
申请号 CN201310431704.X 申请日期 2013.09.22
申请人 浙江大学 发明人 刘兴高;李见会;张明明;孙优贤
分类号 G05D23/00(2006.01)I;G05B13/04(2006.01)I 主分类号 G05D23/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 周烽
主权项 一种径向基的农药生产废液焚烧炉炉温最佳化系统,包括焚烧炉、现场智能仪表、DCS系统、数据接口以及上位机,所述的DCS系统包括控制站和数据库;所述现场智能仪表与DCS系统连接,所述DCS系统与上位机连接,其特征在于:所述的上位机包括:标准化处理模块,用于将从DCS数据库输入的模型训练样本进行预处理,对训练样本中心化,即减去样本的平均值,然后对其进行标准化:计算均值:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mrow><mi>T</mi><mi>X</mi></mrow><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>TX</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000759920550000011.GIF" wi="1319" he="142" /></maths>计算方差:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>&sigma;</mi><mi>x</mi><mn>2</mn></msubsup><mo>=</mo><mfrac><mn>1</mn><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mrow><mo>(</mo><msub><mi>TX</mi><mi>i</mi></msub><mo>-</mo><mover><mrow><mi>T</mi><mi>X</mi></mrow><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000759920550000012.GIF" wi="1292" he="142" /></maths>标准化:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>X</mi><mo>=</mo><mfrac><mrow><mi>T</mi><mi>X</mi><mo>-</mo><mover><mrow><mi>T</mi><mi>X</mi></mrow><mo>&OverBar;</mo></mover></mrow><msub><mi>&sigma;</mi><mi>x</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000759920550000013.GIF" wi="1354" he="153" /></maths>其中,TX<sub>i</sub>为第i个训练样本,是从DCS数据库中采集的生产正常时的关键变量、炉温和使炉温最佳化的操作变量的数据,N为训练样本数,<img file="FDA0000759920550000014.GIF" wi="86" he="68" />为训练样本的均值,X为标准化后的训练样本;σ<sub>x</sub>表示训练样本的标准差,σ<sup>2</sup><sub>x</sub>表示训练样本的方差;模糊系统模块,对从标准化处理模块传过来的标准化后的训练样本X,进行模糊化;设模糊系统中有c<sup>*</sup>个模糊群,模糊群k、j的中心分别为v<sub>k</sub>、v<sub>j</sub>,则第i个标准化后的训练样本X<sub>i</sub>对于模糊群k的隶属度μ<sub>ik</sub>为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&mu;</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>=</mo><msup><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>c</mi><mo>*</mo></msup></munderover><msup><mrow><mo>(</mo><mfrac><mrow><mo>|</mo><mo>|</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><msub><mi>v</mi><mi>k</mi></msub><mo>|</mo><mo>|</mo></mrow><mrow><mo>|</mo><mo>|</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><msub><mi>v</mi><mi>j</mi></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>)</mo></mrow><mfrac><mn>2</mn><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></mfrac></msup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000759920550000015.GIF" wi="1659" he="261" /></maths>式中,m为模糊分类过程中需要的分块矩阵指数,取作2;||·||为范数表达式;使用以上隶属度值或者它的变形以获得新的输入矩阵,对于模糊群k,其输入矩阵变形为:Φ<sub>ik</sub>(X<sub>i</sub>,μ<sub>ik</sub>)=[1 func(μ<sub>ik</sub>) X<sub>i</sub>]   (5)其中func(μ<sub>ik</sub>)为隶属度值μ<sub>ik</sub>的变形函数,取<img file="FDA0000759920550000016.GIF" wi="72" he="79" />或exp(μ<sub>ik</sub>),Φ<sub>ik</sub>(X<sub>i</sub>,μ<sub>ik</sub>)表示第i个输入变量X<sub>i</sub>及其模糊群k的隶属度μ<sub>ik</sub>所对应的新的输入矩阵;径向基神经网络作为模糊系统的局部方程,对每个模糊群进行优化拟合;设第k个径向基模糊系统的的输出是:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mover><mi>y</mi><mo>^</mo></mover><mrow><mi>i</mi><mi>k</mi></mrow></msub><mo>=</mo><munder><mo>&Sigma;</mo><mi>l</mi></munder><msub><mi>w</mi><mrow><mi>l</mi><mi>k</mi></mrow></msub><msub><mi>&Phi;</mi><mrow><mi>l</mi><mi>k</mi></mrow></msub><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><msub><mi>C</mi><mrow><mi>l</mi><mi>k</mi></mrow></msub><mo>|</mo><mo>|</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000759920550000021.GIF" wi="1651" he="122" /></maths>其中C<sub>lk</sub>和w<sub>lk</sub>分别是第k个径向基模糊系统第l个结点的中心和输出权值,Φ<sub>lk</sub>(||X<sub>i</sub>‑C<sub>lk</sub>||)是第k个径向基模糊系统第l个结点的径向基函数,由下式确定:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>&Phi;</mi><mrow><mi>l</mi><mi>k</mi></mrow></msub><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><msub><mi>C</mi><mrow><mi>l</mi><mi>k</mi></mrow></msub><mo>|</mo><mo>|</mo><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>X</mi><mi>i</mi></msub><mo>-</mo><msub><mi>C</mi><mrow><mi>l</mi><mi>k</mi></mrow></msub><mo>|</mo><mo>|</mo><mo>)</mo></mrow><mn>2</mn></msup><mrow><mn>2</mn><mo>&times;</mo><msub><mi>&sigma;</mi><mrow><mi>l</mi><mi>k</mi></mrow></msub></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000759920550000022.GIF" wi="1641" he="197" /></maths>其中σ<sub>lk</sub>是对应的径向基函数的高斯宽度,||·||为范数表达式;最后,由反模糊方法中的重心法得到最后的模糊系统的输出:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mover><mi>y</mi><mo>^</mo></mover><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>c</mi><mo>*</mo></msup></msubsup><msub><mi>&mu;</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub><msub><mover><mi>y</mi><mo>^</mo></mover><mrow><mi>i</mi><mi>k</mi></mrow></msub></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><msup><mi>c</mi><mo>*</mo></msup></msubsup><msub><mi>&mu;</mi><mrow><mi>i</mi><mi>k</mi></mrow></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000759920550000023.GIF" wi="1662" he="204" /></maths><img file="FDA0000759920550000024.GIF" wi="52" he="78" />即为对应于第i个标准化后的训练样本X<sub>i</sub>的炉温预报值和使炉温最佳的操作变量值;所述的上位机还包括:模型更新模块,用于按设定的采样时间间隔,采集现场智能仪表信号,将得到的实测炉温与系统预报值比较,如果相对误差大于10%或炉温超出生产正常上下限范围,则将DCS数据库中生产正常时的使炉温最佳的新数据加入训练样本数据,更新软测量模型;结果显示模块,用于将得到的炉温预报值和使炉温最佳的操作变量值传给DCS系统,在DCS的控制站显示,并通过DCS系统和现场总线传递到现场操作站进行显示;信号采集模块,用于依照设定的每次采样的时间间隔,从数据库中采集数据;所述的关键变量包括进入焚烧炉的废液流量、进入焚烧炉的空气流量和进入焚烧炉的燃料流量;所述的操作变量包括进入焚烧炉的空气流量和进入焚烧炉的燃料流量。
地址 310058 浙江省杭州市西湖区余杭塘路866号