CN109960873A - 一种城市固废焚烧过程二噁英排放浓度软测量方法 - Google Patents
一种城市固废焚烧过程二噁英排放浓度软测量方法 Download PDFInfo
- Publication number
- CN109960873A CN109960873A CN201910224790.4A CN201910224790A CN109960873A CN 109960873 A CN109960873 A CN 109960873A CN 201910224790 A CN201910224790 A CN 201910224790A CN 109960873 A CN109960873 A CN 109960873A
- Authority
- CN
- China
- Prior art keywords
- potential
- indicate
- potential feature
- feature
- dxn
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 69
- 230000008569 process Effects 0.000 title claims abstract description 43
- 238000000691 measurement method Methods 0.000 title claims abstract description 9
- 239000002910 solid waste Substances 0.000 title claims description 11
- HGUFODBRKLSHSI-UHFFFAOYSA-N 2,3,7,8-tetrachloro-dibenzo-p-dioxin Chemical compound O1C2=CC(Cl)=C(Cl)C=C2OC2=C1C=C(Cl)C(Cl)=C2 HGUFODBRKLSHSI-UHFFFAOYSA-N 0.000 title claims 5
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 21
- 230000007246 mechanism Effects 0.000 claims abstract description 15
- 238000005457 optimization Methods 0.000 claims abstract description 13
- 238000005259 measurement Methods 0.000 claims abstract description 12
- 238000004458 analytical method Methods 0.000 claims abstract description 9
- 230000003044 adaptive effect Effects 0.000 claims description 16
- 239000011159 matrix material Substances 0.000 claims description 14
- 238000013507 mapping Methods 0.000 claims description 10
- 238000000605 extraction Methods 0.000 claims description 8
- 238000000465 moulding Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 2
- 239000000284 extract Substances 0.000 claims description 2
- 230000017105 transposition Effects 0.000 claims description 2
- KVGZZAHHUNAVKZ-UHFFFAOYSA-N 1,4-Dioxin Chemical compound O1C=COC=C1 KVGZZAHHUNAVKZ-UHFFFAOYSA-N 0.000 description 65
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 21
- 239000003546 flue gas Substances 0.000 description 21
- 238000001514 detection method Methods 0.000 description 8
- 238000012549 training Methods 0.000 description 7
- 239000003517 fume Substances 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 239000007789 gas Substances 0.000 description 5
- 238000010521 absorption reaction Methods 0.000 description 4
- 238000013528 artificial neural network Methods 0.000 description 4
- 238000010276 construction Methods 0.000 description 4
- 239000003344 environmental pollutant Substances 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 231100000719 pollutant Toxicity 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000000007 visual effect Effects 0.000 description 4
- 239000000428 dust Substances 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000003062 neural network model Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 235000008733 Citrus aurantifolia Nutrition 0.000 description 1
- 241001269238 Data Species 0.000 description 1
- RWSOTUBLDIXVET-UHFFFAOYSA-N Dihydrogen sulfide Chemical compound S RWSOTUBLDIXVET-UHFFFAOYSA-N 0.000 description 1
- 238000012952 Resampling Methods 0.000 description 1
- 235000011941 Tilia x europaea Nutrition 0.000 description 1
- 229910052934 alunite Inorganic materials 0.000 description 1
- 239000010424 alunite Substances 0.000 description 1
- 239000002956 ash Substances 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 231100000693 bioaccumulation Toxicity 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 235000019504 cigarettes Nutrition 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229910001385 heavy metal Inorganic materials 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 239000004571 lime Substances 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000002893 slag Substances 0.000 description 1
- 239000000779 smoke Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 231100000167 toxic agent Toxicity 0.000 description 1
- 239000003440 toxic substance Substances 0.000 description 1
- KPZTWMNLAFDTGF-UHFFFAOYSA-D trialuminum;potassium;hexahydroxide;disulfate Chemical compound [OH-].[OH-].[OH-].[OH-].[OH-].[OH-].[Al+3].[Al+3].[Al+3].[K+].[O-]S([O-])(=O)=O.[O-]S([O-])(=O)=O KPZTWMNLAFDTGF-UHFFFAOYSA-D 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/30—Administration of product recycling or disposal
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/0265—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion
- G05B13/027—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion using neural networks only
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/048—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators using a predictor
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
- G06N20/10—Machine learning using kernel methods, e.g. support vector machines [SVM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
- G06N20/20—Ensemble learning
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02W—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO WASTEWATER TREATMENT OR WASTE MANAGEMENT
- Y02W90/00—Enabling technologies or technologies with a potential or indirect contribution to greenhouse gas [GHG] emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Automation & Control Theory (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Economics (AREA)
- Sustainable Development (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Feedback Control In General (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种基于多源潜在特征选择性集成(SEN)建模的DXN排放浓度软测量方法。首先,将MSWI过程数据依据工业流程划分为不同来源的子系统,采用主元分析(PCA)分别提取其潜在特征,并依据经验预设的主元贡献率阈值进行多源潜在特征初选;接着,采用互信息(MI)度量初选的潜在特征与DXN间的相关性,自适应确定潜在特征再选的上下限及阈值;最后,基于再选潜在特征,采用具有超参数自适应选择机制的最小二乘‑支撑向量机(LS‑SVM)算法,建立针对不同子系统的DXN排放浓度子模型,采用基于分支定界(BB)和预测误差信息熵加权算法的策略优化选择子模型和计算权系数,构建DXN排放浓度SEN软测量模型。
Description
本发明由科学技术部国家重点研发计划(No:2018YFC1900801)和国家自然科学基金(No:61573364,61873009)资助。
技术领域
本发明属于固废焚烧技术领域,尤其涉及一种城市固废焚烧过程二噁英排放浓度软测量方法。
背景技术
基于运行优化控制策略降低复杂工业过程的能源消耗和污染排放,是国内外众多流程工业企业所面临的急需解决的难题[1,2,3]。焚烧是进行城市固废(MSW)处理的主要技术手段[4]。对于发展中国家的MSW焚烧(MSWI)企业,最为紧迫的问题是如何控制焚烧造成的污染排放[5,6],其中急需控制得是:造成焚烧建厂“邻避效应”、具有生物累积效应等特点的剧毒物质二噁英(DXN)的排放浓度[7,8,9]。MSWI企业最主要关注如何基于优化的运行参数实现最小化的DXN排放[10]。目前,除采用先进的尾气处理装置外,工业过程普遍采用如下间接策略控制DXN排放,即“3T1E”准则[11],即:焚烧炉内高于850℃的温度(T)、超过2秒(T)的烟气停留时间、较大的湍流程度(T)和合适的过量空气系数(E)。当前MSWI企业不能进行以降低DXN排放为直接目标的运行优化和反馈控制,其主要原因:一是DXN排放浓度的机理模型难以构建[12],二是以月或季为周期的离线直接检测不能提供实时反馈的DXN排放浓度[13]。近年来的研究热点是基于指示物/关联物对DXN排放进行在线间接检测[14,15],但这些方法固有的设备复杂、造价昂贵、检测滞后等原因导致其仍难以用于MSWI过程的运行优化和反馈控制。
数据驱动软测量技术可用于需要离线化验的难以检测参数(如本发明中的二噁英)的在线估计[16,17]。针对MSWI过程,文献[18,19]通过采用多年前欧美研究机构针对不同类型的焚烧炉采集的小样本数据,基于线性回归、人工神经网络(ANN)等算法构建DXN排放浓度软测量模型。近年来,文献[20]结合相关性分析、主元分析(PCA)和ANN等算法,基于MSWI过程数据构建DXN排放预测模型;但是,ANN并不适用于构建DXN浓度排放模型,主要原因在于其固有的易陷入局部最小、易过拟合和对小样本数据建模泛化性能差等缺点。采用具有合适超参数的支撑向量机(SVM)算法可有效地用于小样本数据建模[21]。针对SVM需求解二次规划的问题,最小二乘-支撑向量机(LS-SVM)通过求解线性等式予以克服;模型超参数可以通过单目标或多目标优化求解算法获得[22,23,24],但这些方法耗时且只能得到次优解[25]。因此,目前研究缺少对超参数进行有效选择的自适应机制。
通常,MSWI过程包括固废储运、固废焚烧、蒸汽发电和烟气处理等阶段组成的多个子系统,所涉及的过程变量达到数百维,与DXN的产生、吸收和再合成等机理过程在不同程度上具有相关性。文献[26]指出,模型输入维数的增加和低价值训练样本的增加使得完备数量的训练样本难以获得。在模式识别领域,通常认为训练样本的数量与特征之比应该为2、5或10。文献[27]定义了维数约简后的训练样本与约简特征之比,认为该值应满足构建鲁棒预测学习模型的要求。因此,针对具有小样本高维特性的DXN排放浓度建模数据,进行维数约简是非常必要的。基于非监督特征提取的方法获得的特征,虽蕴含原始高维变量中的主要变化却可能导致所提取特征与被预测参数无关。能够提取高维数据蕴含变化的主元分析(PCA),是目前工业过程在难以检测参数软测量中最为常用的潜在特征提取方法[28],但贡献率低的主元会导致预测稳定性差[29]。
从另外一个视角,针对MSWI过程的不同子系统所提取的特征可视为源于多个视图的多源信息。理论和经验分析表明,面向多源信息采用选择性集成(SEN)学习机制构建的软测量模型具有更佳的预测稳定性和鲁棒性,其中子模型间的差异性尤为重要[29]。文献[30]综述了集成子模型多样性的构造策略,指出训练样本重采样包括划分训练样本(样本空间)、划分或变换特征变量(特征空间)等,在模型预测性能上基于特征空间的集成构造策略优于基于多分类器的构造策略。针对小样本多源高维谱数据,汤等人提出基于选择性融合多源特征和多工况样本的SEN潜结构映射模型[31,32];文献[33]和[32]提出了基于随机采样样本空间的SEN神经网络模型和潜结构映射模型,文献[34]提出基于子空间的集成学习通用框架,文献[35]提出了在特征子空间内随机采样样本空间的面向多尺度机械信号的双层SEN潜结构映射模型,文献[36]提出的SEN神经网络模型分别构建优化视角的候选子模型和优化选择集成子模型及其权重,但上述方法未进行模型参数自适应机制的研究。综上可知,以度量和再选后的非监督潜在特征为输入,基于自适应超参数选择机制的SEN-LSSVM建模策略及其在DXN排放浓度软测量中的研究未见报道。
发明内容
综上,本发明提供一种基于多源潜在特征SEN的MSWI过程DXN排放浓度软测量方法,首先,采用潜在特征提取与初选模块将MSWI过程数据依据工业流程划分为不同来源的子系统,采用主元分析(PCA)分别提取其潜在特征并依据经验预设的主元贡献率阈值进行多源潜在特征初选。接着,潜在特征度量与再选模块采用互信息(MI)度量初选的潜在特征与DXN间的相关性,自适应确定潜在特征再选的上下限及阈值。最后,自适应选择性集成建模模块基于再选潜在特征,采用具有超参数自适应选择机制的最小二乘-支撑向量机(LS-SVM)算法,建立针对不同子系统的DXN排放浓度子模型,采用基于分支定界(BB)和预测误差信息熵加权算法的策略优化选择子模型和计算权系数,构建DXN排放浓度SEN软测量模型。
附图说明
图1基于DXN视角的MSWI过程描述;
图2基于潜在特征SEN建模的DXN排放浓度软测量策略;
图3不同子系统的前6个PC的累积贡献率;
图4(a)Incinerator子系统初选潜在特征与DXN间的MI值;
图4(b)Boiler子系统初选潜在特征与DXN间的MI值;
图4(c)Flue gas子系统初选潜在特征与DXN间的MI值;
图4(d)Steam子系统初选潜在特征与DXN间的MI值;
图4(e)Stack子系统初选潜在特征与DXN间的MI值;
图4(f)Common子系统初选潜在特征与DXN间的MI值;
图4(g)MSWI系统初选潜在特征与DXN间的MI值;
图5(a)焚烧子模型超参数自适应寻优的第1次和第2次的曲线;
图5(b)锅炉子模型超参数自适应寻优的第1次和第2次的曲线;
图5(c)烟气处理子模型超参数自适应寻优的第1次和第2次的曲线;
图5(d)蒸汽发电子模型超参数自适应寻优的第1次和第2次的曲线;
图5(e)烟气排放子模型超参数自适应寻优的第1次和第2次的曲线;
图5(f)公用工程子模型超参数自适应寻优的第1次和第2次的曲线;
图5(g)MSWI全流程子模型超参数自适应寻优的第1次和第2次的曲线。
具体实施方试
面向DXN排放过程的MSWI描述
MSWI的主要设备包括焚烧炉、移动炉排、废锅和尾气处理等设备,其中:焚烧炉将MSW转化为残渣、灰尘、烟气与热量,位于焚烧炉底部的移动炉排促使MSW有效和完全燃烧,废锅产生的蒸汽用于推动汽轮机产生电力,烟气中的灰尘和污染物通过尾气处理设备净化后排入大气。其过程如图1所示。
由图1可知,从污染排放视角,MSWI过程包含了DXN产生、吸收和排放共3个阶段,其分别包含在标记为G1、G2和G3的烟气中。显然这些不同阶段的烟气中所包含的DXN浓度具有差异性。通常,为保证焚烧炉内的有害物质能够有效和完全分解,烟气温度至少应该达到850℃并保持2秒以上。用于焚烧MSW的一次风从炉排底部喷入的同时又对炉排进行冷却,并通过引入湍流和保证过量氧的供应使得二次风能够辅助进行烟气的完全燃烧。在烟气冷却过程中,进行焚烧矿渣和废锅底灰的处置和收集,同时排出烟气G1。活性炭和石灰被注入反应器,用于移除酸性气体和吸收DXN及一些重金属,烟气然后再进入袋式过滤器。在反应器和袋式过滤器内产生的飞灰被注入混涅设备,同时产生烟气G2。引风机将烟气G2吸入到烟囱,进而作为烟气G3排放至空气中,包含HCL、SO2、NOx和HF等多种能够实时在线检测的污染物浓度,以及具有长周期、高成本等特点需离线化验浓度的DXN。由上述描述可知,DXN排放浓度与MSWI过程不同阶段的易检测过程变量均具有相关性。
由图1可知,MSWI过程可以分为焚烧、锅炉、烟气处理、蒸汽发电、烟气排放、公共工程辅助共6个子系统。对于DXN排放浓度建模,本发明将多个子系统视为多源信息。
建模策略
本发明中,模型输入数据X∈RN×M包括N个样本(行)和M个变量(列),其源于MSWI过程的不同子系统。将来自第ith个子系统的建模数据表示为即存在如下关系,
其中,I表示子系统个数,Mi表示第ith个子系统包含的变量个数。相应的,输出数据包括N个样本(行),其来源于离线化验的DXN排放浓度检测数据。显然,输入/输出数据在时间尺度上具有较大的差异性:过程变量以秒为单位在DCS系统采集与存储,DXN排放浓度以月/季为周期离线化验获得,故存在N<<M。
依据上述情况,本发明提一种基于潜在特征SEN建模的DXN排放浓度软测量方法,包括潜在特征提取与初选模块、潜在特征度量与再选模块、自适应选择性集成建模模块,如图2所示。
在图2中,表示从第ith个子系统所采集的全部过程变量;表示针对第ith个子系统的全部过程变量采用PCA提取的数量为的全部潜在特征;表示针对第ith个子系统的全部潜在特征依据设定阈值θContri选择的数量为的初选潜在特征;表示针对第ith个子系统的初选潜在特征与DXN进行互信息(MI)度量后,基于阈值θMI选择的数量为的再选潜在特征;和表示为第ith个基于LS-SVM的子模型所选择的核参数和正则化参数,即超参数对,本发明将其记为表示第ith个子模型的预测输出;y和表示DXN排放浓度软测量模型的真值和预测输出。
上述模块的功能是:
(1)潜在特征提取与初选模块:采用PCA提取从不同子系统所采集的过程变量的全部潜在特征,基于依据经验设定的潜在特征贡献率阈值获得多源初选潜在特征,其目的是防止较小贡献率的潜在特征造成模型预测性能不稳定性。
(2)潜在特征度量与再选模块:采用互信息(MI)度量不同子系统的初选潜在特征与DXN排放浓度间的关系,并结合基于软测量模型预测性能自适应确定的阈值获得再选潜在特征,其目的是使得所选的多源潜在特征与DXN排放浓度间具有较好的映射关系。
(3)自适应选择性集成建模模块:采用超参数自适应选择策略构建面向不同子系统的具有最佳预测性能的子模型,结合分支定界(BB)和预测误差信息熵加权算法自适应地选择子模型和计算其加权系数,其目的是选择具有较好冗余与互补关系的潜在特征构建子模型进行融合,以提高SEN软测量模型的预测性能。
潜在特征提取与初选模块
以第ith个子系统为例,首先采用PCA提取高维输入过程变量的潜在特征。将输入数据Xi进行零均值1方差的标准化后,将其分解为,
其中,和表示第个主元(PC)的得分和负载向量,T表示转置,表示对第ith个子系统所提取的潜在特征数量,其计算公式如下,
基于上述表达式,从数据Xi所提取的全部潜在特征可表示为,
其中,表示得分矩阵,是数据Xi在负载矩阵Pi方向上的正交映射;Pi采用如下公式表示,
其中,
因此,从数据Xi所提取的潜在特征可表示为,
其中,
进一步,全部潜在特征可表示为,
研究表明,采用贡献率较小的潜在变量建模会导致模型预测性能的不稳定。此处,将与第个负载向量对应的特征向量记为相应的第个潜在特征的贡献率采用如下公式计算,
将依据经验选择的阈值记为θContri,其默认取值为1。采用如下规则对全部潜在特征进行初次选择,
其中,表示第个潜在特征是否被选中的标记值,其值为1表示该潜在特征被初次选中。
因此,将针对第ith个子系统的初选潜在特征表示为,
进一步,全部初选潜在特征ZFeSe1st可表示为,
潜在特征度量与再选模块
上一步骤中所获得的初选潜在特征是采用非监督方式提取,并且同一子系统所包含的特征是相互独立的,但未考虑这些特征与DXN排放浓度间的相关性,即贡献率高的潜在特征并不一定与DXN间的相关性强。仍以第ith个子系统为例,将每个初选潜在特征与DXN排放浓度间的互信息(MI)值标记为采用如下公式计算,
其中,和pprob(y)表示和y的边际概率密度;表示联合概率密度;表示条件熵,H(y)表示信息熵。
依据软测量模型的预测性能自适应确定阈值。阈值的上限值下限值及固定步长采用如下公式计算,
其中,函数max(·)和min(·)分别表示取最大值和最小值;表示依据经验确定的候选阈值的数量,其默认值为10。
将选定的阈值记为θContri,其值在和之间以DXN软测量模型的预测性能为准则进行自适应选择。
采用如下规则对初选的潜在特征进行再次选择,
其中,表示第个潜在特征是否被选中的标记值,其值为1表示该潜在特征被再次选中。
进一步,将针对第ith个子系统的再选潜在特征表示为,
因此,全部再选潜在特征ZFeSe2nd可表示为,
自适应选择性集成建模模块
以第ith个子系统为例,描述基于再选潜在特征和模型超参数对构建DXN排放浓度子模型的过程。
首先,将再选潜在特征通过映射变换到高维特征空间,然后求解如下优化问题,
其中,wi表示权重系数,bi表示偏置,是第nth个样本的预测误差。
采用拉格朗日方法,可得如下公式,
其中,表示拉格朗日算子向量,表示预测误差向量。
对上述公式进行求解,
将所采用的核函数表示如下,
进一步,将LS-SVM问题转换为求解以下线性等式系统,
通过求解上述公式,得到βi和bi。
进而,基于LS-SVM构建的DXN排放浓度子模型可表示为,
上述排放浓度子模型的超参数自适应选择机制采用下述的两步法实现:
第1步,采用网格搜索策略以子模型的预测性能为目标函数,在候选超参数矩阵Mpara中自适应选择初始超参数对超参数矩阵Mpara如下所示,
其中,k=1,…,K,K表示候选核参数的数量;r=1,…,R,R表示候选惩罚参数的数量;表示由第kth个核参数和第rth个惩罚参数组成的超参数对,也超参数矩阵Mpara中的第jth个参数对,即存在j=1,…,J,J=K×R表示超参数矩阵Mpara中的全部超参数对的个数。因此,对于初次采用网格搜索策略所选择的超参数对是矩阵Mpara中一个元素,即存在
第2步,基于上述方法选择的采用如下公式获得新的候选超参数集合,
其中,和表示新的候选超参数集合,分别对应核参数向量和惩罚参数向量;Nker和Nreg表示依据经验设定的新的超参数的数量;和为依据经验设定的超参数收缩和扩放因子,其默认值均为10。
通过再次采用网格搜索策略自适应获得第ith个子模型的超参数对
对全部子系统执行上述过程,子模型预测输出的集合可表示为,
其中,表示第ith个子模型。
结合基于分支定界(BB)的最优化选择算法和基于误差的信息熵加权算法对上述子模型进行自适应的优化选择和计算加权系数。在给定候选子模型和加权算法后,最佳子模型选择与加权类似最优特征选择[29]。面向有限数量的候选子模型,通过多次运行优化和加权算法,可获得集成尺寸为2到(I-1)的SEN模型,最后排序这些优化的SEN模型并将预测性能最佳的作为最终DXN软测量模型。
假定最终DXN软测量模型的集成尺寸为Isel,其预测输出值可由下式计算:
其中,表示第iselth个优化选择的子模型,和表示其对应的加权系数及预测值,和表示子模型的超参数和输入特征。
对比式(29)可知,存在如下关系,
利用子模型的预测值和真值,采用基于预测误差信息墒的加权算法得到,如下所示,
其中,
其中,和表示第nth个样本基于第iselth个优化选择的子模型的预测值和相对预测误差,表示第iselth个优化选择的子模型的预测误差的信息熵。
应用研究
本发明中的建模数据源于北京某基于炉排炉的MSWI焚烧企业的1#炉,涵盖了2012~2018年所记录可用的DXN排放浓度检测样本,其数量为39个;相应的输入变量的维数为286维(包含了MSWI过程的全部过程变量)。可见,输入特征数量远远超过建模样本数量,进行维数约简非常有必要。本发明将建模数据等分为两部分,分别用于训练和测试。
本发明中,将6个子系统焚烧、锅炉、烟气处理、蒸汽发电、烟气排放和公用工程辅助分别的标记为Incinerator、Boiler、Flue gas、Steam、Stack和Common。为表示焚烧过程变量的整体变化特性,本发明将包含全部变量的MSWI系统也作为一个特殊的子系统进行分析和建模。因此,本发明共包含7个子系统。
潜在特征提取与初选结果
采用PCA提取的7个子系统的前6个潜在特征的累积贡献率如图3所示。
图3表明,前6个PC的贡献率达到了80%,并且不同子系统的潜在变量所提取的贡献率具有差异性。
基于单个PC的贡献率不小于1%的准则,初选潜在特征的主元个数及其贡献率如表1所示。
表1初选潜在特征的主元个数及其贡献率
由表1可知,不同子系统所初选的潜在特征的数量为13、6、9、8、6、12和15。由于PCA属于非监督特征提取方法,这些所提取特征仅是描述了输入数据的变化,其与DXN间的映射关系需要进行进一步的度量。
潜在特征度量与再选结果
采用MI方法度量针对不同子系统所提取的初选潜在特征与DXN间的映射关系,如图4所示。
图4表明:(1)全部子系统所初选的能够表征所在子系统过程变量最大变化的第1个潜在特征的MI值是最小的,表明了这些潜在特征与DXN排放浓度间的相关性较弱;(2)除第1个潜在特征外,其他潜在特征在表征过程变量的贡献率上虽然是逐渐下降的,但在表征MI值上却没有明显的规律可循;(3)从机理上分析,与DXN的产生、吸收和排放最为相关的是Incinerator、Flue gas和Stack子系统,但这些子系统在所提取的潜在特征的MI值与其他子系统的差别较小。可见,仅仅是基于MI值进行决策存在局限性。不同子系统的初选潜在特征的MI值的最大值和最小值如表2所示。
表2不同子系统的初选潜在特征的MI值的最大值和最小值集合
表2表明:(1)针对最大值集合:最大值源于理论上与DXN排放并无直接关系的公共辅助(Common)子系统,值为0.8613,其是否合理有待于进一步的结合模型预测结果进行验证;排在第2位的是焚烧(Incinerator)子系统,其值为0.8559,该潜在变量理论上与DXN的产生相关,是比较合理的;(2)针对最小值集合:最小值源于焚烧(MSWI)子系统,仅为0.4429,表明针对不同的子系统进行单独分析还是比较必要的;最大值源于烟气排放(Stack)子系统,其为0.7182,由于其他的排放气体与DXN间是存在相关性的,此值也是较为合理的。
由表2可知,MI阈值的上限为0.7882,下限为0.7182,步长为0.006999。结合阈值的上下限和步长,最终确定的阈值为0.7882,再选的潜在特征的数量和MI值如表3所示。
表3再选潜在特征数量和MI值
自适应选择集成建模
本发明中,将候选正则化参数与核参数的集合分别预先选择为{0.0001,0.001,0.01,0.1,1,10,100,1000,2000,4000,6000,8000,10000,20000,40000,60000,80000,160000}和{0.0001,0.001,0.01,0.1,1,10,100,1000,1600,3200,6400,12800,25600,51200,102400}。
结合上文可知,焚烧、锅炉、烟气处理、蒸汽发电、烟气排放、公用工程、MSWI全流程子模型的输入特征数量分别为5、2、1、3、2、6和1。采用网格搜索方法进行超参数自适应寻优的第1次和第2次的曲线如图5所示。
基于上述结果,上述子模型自适应选择的超参数对分别为{109,109}、{10000,25.75}、{5.950,0.0595}、{30.70,2.080}、{5.950,0.5950}、{1520800,22816}和{1362400,158.5},对应的测试数据的均方根误差(RMSE)分别为0.01676、0.02302、0.01348、0.01943、0.01475、0.02261和0.02375。
采用基于BB和预测误差信息熵加权算法的寻优和加权策略,在集成尺寸为2~6时所构建的SEN模型的测试误差分别为0.01345、0.01332、0.01401、0.01460和0.01560。最终DXN软测量模型的集成尺寸为3,其选择的子模型所对应的子系统为烟气处理、烟气排放和焚烧,理论上这3个子系统与DXN的吸收、排放和生成相关。从本发明结果看,验证了所以算法的有效性,也表明了数据的可用性。
比较结果
与通常采用PLS单模型、PCA-LSSVM单模型以及不同加权方法的比较如表4所示。
表4不同建模方法的结果统计
表1表明,采用全部过程变量构建的基于PLS和PCA-LSSVM的DXN单模型的预测性能均弱于本发明所提的SEN建模方法,说明采用基于源特征构建SEN模型的策略是有效;同时,集成全部子模型的方法出采用PLS加权具有强于其他集成全部子模型的EN方法,表明PLS算法在消除子模型的共线性方面较好;此外,SEN模型选择的子模型所对应的子系统均与DXN的产生、吸收和排放机理相关,表明了建模数据的可用性和算法的有效性。
本发明基于北京某焚烧企业的工业过程数据,采用基于PCA和先验知识的潜在特征提取和初选、基于MI和先验知识的初选潜在特征度量和选择和面向再选潜在特征的自适应SEN建模机制,提出了基于多源潜在特征SEN建模的DXN排放浓度软测量,仿真验证了所提方法的有效性。在结合软测量模型的预测出差自适应调整贡献率阈值、MI阈值、超参数和SEN模型结构等方面还有待于深入研究。此外,结合DXN排放机理的分析也有待于深入进行。
参考文献
[1]柴天佑。复杂工业过程运行优化与反馈控制.自动化学报,39(11)(2013)1744-1757.
[2]T.Y.Chai,J.L.Ding,G.Yu,H.Wang.Integrated optimization for theautomation systems of mineral processing[J].IEEE Transactions on AutomationScience&Engineering,11(4)(2014)965-982.
[3]T.Y.Chai,S.J.Qin,H.Wang.Optimal operational control for complexindustrial processes[J].Annu.Rev.Control,38(1)(2014)81-92.
[4]Arafat H A,Jijakli K,Ahsan A,2015.Environmental performance andenergy recovery potential of five processes for municipal solid wastetreatment[J].Journal of Cleaner Production,2015,105:233-240
[5]Yuanan H,Hefa C,Shu T.The growing importance of waste-to-energy(WTE)incineration in China's anthropogenic mercury emissions:Emissioninventories and reduction strategies[J].Renewable and Sustainable EnergyReviews,2018,97:119-137.
[6]Huang T,Zhou L,Liu L,Xia M.Ultrasound-enhanced electrokineticremediation for removal of Zn,Pb,Cu and Cd in municipal solid wasteincineration fly ashes[J].Waste Management,2018,75:226-235.
[7]P.H.Jones,J.Degerlache,E.Marti,G.Mischer,M.C.Scherrer,W.J.Bontinck,H.J.Niessen,The global exposure of man to dioxins-a perspectiveon industrial-waste incineration[J],Chemosphere 26(1993)1491-1497.
[8]Li X,Zhang C,Li Y,Zhi Q.The Status of Municipal Solid WasteIncineration(MSWI)in China and its Clean Development[J].Energy Procedia,2016,104:498-503.
[9]Phillips K,Longhurst P J,Wagland S T.Assessing the perception andreality of arguments against thermal waste treatment plants in terms ofproperty prices[J].Waste Management.2014,34(1):219-225.
[10]Zhang H J,Ni Y W,Chen J P,Zhang Q.Influence of variation in theoperating conditions on PCDD/F distribution in a full-scale MSW incinerator[J].Chemosphere,2008,70(4):721-730.
[11]Mukherjee A,Debnath B,Ghosh S K.A Review on Technologies ofRemoval of Dioxins and Furans from Incinerator Flue Gas[J].ProcediaEnvironmental Sciences,2016,35:528-540.
[12]B.R.Stanmore,Modeling the formation of PCDD/F in solid wasteincinerators,Chemosphere 47(2002)565-773.
[13]乔俊飞,郭子豪,汤健.面向城市固废焚烧过程的二噁英排放浓度检测方法综述[J].自动化学报,2019(在审).
[14]M.Pandelova,D.Lenoir,and K.W.Schramm,Correlation between PCDD/F,PCB and PCBz in coal/waste combustion Influence of various inhibitors,Chemosphere 62(2006)1196-1205.
[15]B.K.Gullett,L.Oudejans,D.Tabor,A.Touati,and S.Ryan,Near-real-timecombustion monitoring for PCDD/PCDF indicators by GC–REMPI–TOFMS[J][J],Environmental Engineering Science 46(2012)923-928.
[16]W.Wang,T.Y.Chai,W.Yu,Modeling component concentrations of sodiumaluminate solution via hammerstein recurrent neural networks[J],IEEETransactions on Control Systems Technology,20(2012)971-982.
[17]J.Tang,T.Y.Chai,W.Yu,L.J.Zhao,Modeling load parameters of ballmill in grinding process based on selective ensemble multisensor information[J],IEEE Transactions on Automation Science&Engineering,10(2013)726-740.
[18]Chang N B,Huang S H.Statistical modelling for the prediction andcontrol of PCDDs and PCDFs emissions from municipal solid waste incinerators[J].Waste Management&Research,1995,13,379-400.
[19]Chang N B,Chen W C.Prediction of PCDDs/PCDFs emissions frommunicipal incinerators by genetic programming and neural network modeling[J].Waste Management&Research,2000,18(4)41-351.
[20]Bunsan S,Chen W Y,Chen H W,Chuang Y H,Grisdanurak N.Modeling thedioxin emission of a municipal solid waste incinerator using neural networks[J].Chemosphere,2013,92:258-264.
[21]T.A.F.Gomes,R.B.C.Prudêncio,C.Soares,A.L.D.Rossi,A.Carvalho,Combining meta-learning and search techniques to select parameters forsupport vector machines[J].Neurocomputing 75(1)(2012)3-13.
[22]J.Tang,T.Y.Chai,W.Yu and L.J.Zhao,Feature extraction andselection based on vibration spectrum with application to estimate the loadparameters of ball mill in grinding process[J],Control Engineering Practice,20(10)(2012)991-1004.
[23]C.Soares,.A hybrid meta-learning architecture for multi-objectiveoptimization of svm parameters[J].Neurocomputing,143(143)(2014)27-43.
[24]G.Yu,T.Y.Chai,X.C.Luo.Multiobjective production planningoptimization using hybrid evolutionary algorithms for mineral processing[J].IEEE Transact.Evolut.Comput.15(4)(2011)487-514.
[25]S.Yin,J.Yin,Tuning kernel parameters for svm based on expectedsquare distance ratio[J].Information Sciences,370-371(2016)92-102.
[26]Li D C,Liu C W.Extending attribute information for small data setclassication.IEEE Transactions on Knowledge and Data Engineering,2010,24(3):452-464
[27]汤健,乔俊飞,柴天佑,等.基于虚拟样本生成技术的多组分机械信号建模[J].自动化学报,2018,44(9):1569-1589.
[28]Jian Tang,Zhuo liu,Jian Zhang,Zhiwei Wu,Tianyou Chai,WenYu.Kernel latent feature adaptive extraction and selection method for multi-component non-stationary signal of industrial mechanical device[J],Neurocomputing,2016,216(C):296-309.
[29]汤健,田福庆,贾美英,李东。基于频谱数据驱动的旋转机械设备负荷软测量[M],国防工业出版社,2015年6月,北京
[30]Gavin Brown,Jeremy Wyatt,Rachel Harris,Xin Yao.Diversity creationmethods:a survey and categorisation[J].Information Fusion 6(2005)5–20
[31]Tang Jian,Chai Tian-you,Yu Wen,Zhao Li-jie.Modeling LoadParameters of Ball Mill in Grinding Process Based on Selective EnsembleMultisensor Information[J].IEEE Tranzaction on Automation Science andEngineering,2013,10(3):726-740.
[32]Jian Tang,Tianyou Chai,Wen Yu,Zhuo Liu,Xiaojie Zhou.A Comparativestudy that measures ball mill load parameters through different single-scaleand multi-scale frequency spectra-based approaches[J],IEEE Transactions onIndustrial Informatics,2016,12(6):2008-2019.
[33]Zhou ZH,Wu J and Tang W,Ensembling neural networks:many could bebetter than all,”Artificial Intelligence[J],2002,137(1-2):239-263.
[34]Ma G,Wu L,Wang Y.A general subspace ensemble learning frameworkvia totally-corrective boosting and tensor-based and local patch-basedextensions for gait recognition[J].Pattern Recognition,2017.
[35]Jian Tang,Junfei Qiao,ZhiWei Wub,Tianyou Chai,Jian Zhang,WenYu.Vibration and acoustic frequency spectra for industrial process modelingusing selective fusion multi-condition samples and multi-sourcefeatures.Mechanical Systems and Signal Processing,99,1-27,2018.
[36]S.Soares,C.H.Antunes,R.Araújo,“Comparison of a genetic algorithmand simulated annealing for automatic neural network ensemble development,”Neurocomputing,Vol.21,No.9,pp.498-511,May.2013.
Claims (4)
1.一种城市固废焚烧过程二噁英排放浓度软测量方法,其特征在于,首先,采用潜在特征提取与初选模块将MSWI过程数据依据工业流程划分为不同来源的子系统,采用主元分析(PCA)分别提取其潜在特征,并依据经验预设的主元贡献率阈值进行多源潜在特征初选;接着,潜在特征度量与再选模块采用互信息(MI)度量初选的潜在特征与DXN间的相关性,自适应确定潜在特征再选的上下限及阈值;最后,自适应选择性集成建模模块基于再选潜在特征,采用具有超参数自适应选择机制的最小二乘-支撑向量机(LS-SVM)算法,建立针对不同子系统的DXN排放浓度子模型,采用基于分支定界(BB)和预测误差信息熵加权算法的策略优化选择子模型和计算权系数,构建DXN排放浓度SEN软测量模型。
2.如权利要求1所述的城市固废焚烧过程二噁英排放浓度软测量方法,其特征在于,潜在特征提取与初选模块的工作过程如下:以第ith个子系统为例,首先采用PCA提取高维输入过程变量的潜在特征。将输入数据Xi进行零均值1方差的标准化后,将其分解为,
其中,和表示第个主元(PC)的得分和负载向量,T表示转置,表示对第ith个子系统所提取的潜在特征数量,其计算公式如下,
基于上述表达式,从数据Xi所提取的全部潜在特征可表示为,
其中,表示得分矩阵,是数据Xi在负载矩阵Pi方向上的正交映射;Pi采用如下公式表示,
其中,
因此,从数据Xi所提取的潜在特征可表示为,
其中,
进一步,全部潜在特征可表示为,
将与第个负载向量对应的特征向量记为相应的第个潜在特征的贡献率采用如下公式计算,
将依据经验选择的阈值记为θContri,其默认取值为1,采用如下规则对全部潜在特征进行初次选择,
其中,表示第个潜在特征是否被选中的标记值,其值为1表示该潜在特征被初次选中,
因此,将针对第ith个子系统的初选潜在特征表示为,
进一步,全部初选潜在特征ZFeSe1st可表示为,
3.如权利要求2所述的城市固废焚烧过程二噁英排放浓度软测量方法,其特征在于,潜在特征度量与再选模块的工作过程如下:以第ith个子系统为例,将每个初选潜在特征与DXN排放浓度间的互信息(MI)值标记为采用如下公式计算,
其中,和pprob(y)表示和y的边际概率密度;表示联合概率密度;表示条件熵,H(y)表示信息熵,
依据软测量模型的预测性能自适应确定阈值。阈值的上限值下限值及固定步长采用如下公式计算,
其中,函数max(·)和min(·)分别表示取最大值和最小值;表示依据经验确定的候选阈值的数量,其默认值为10,
将选定的阈值记为θContri,其值在和之间以DXN软测量模型的预测性能为准则进行自适应选择,
采用如下规则对初选的潜在特征进行再次选择,
其中,表示第个潜在特征是否被选中的标记值,其值为1表示该潜在特征被再次选中,
进一步,将针对第ith个子系统的再选潜在特征表示为,
因此,全部再选潜在特征ZFeSe2nd可表示为,
4.如权利要求3所述的城市固废焚烧过程二噁英排放浓度软测量方法,其特征在于,自适应选择性集成建模模块的工作过程如下:
以第ith个子系统为例,首先,将再选潜在特征通过映射变换到高维特征空间,然后求解如下优化问题,
其中,wi表示权重系数,bi表示偏置,是第nth个样本的预测误差,
采用拉格朗日方法,可得如下公式,
其中,表示拉格朗日算子向量,表示预测误差向量。
对上述公式进行求解,
将所采用的核函数表示如下,
进一步,将LS-SVM问题转换为求解以下线性等式系统,
通过求解上述公式,得到βi和bi,
进而,基于LS-SVM构建的DXN排放浓度子模型可表示为,
上述排放浓度子模型的超参数自适应选择机制采用下述的两步法实现:
第1步,采用网格搜索策略以子模型的预测性能为目标函数,在候选超参数矩阵Mpara中自适应选择初始超参数对超参数矩阵Mpara如下所示,
其中,k=1,…,K,K表示候选核参数的数量;r=1,…,R,R表示候选惩罚参数的数量;表示由第kth个核参数和第rth个惩罚参数组成的超参数对,也超参数矩阵Mpara中的第jth个参数对,即存在j=1,…,J,J=K×R表示超参数矩阵Mpara中的全部超参数对的个数,因此,对于初次采用网格搜索策略所选择的超参数对是矩阵Mpara中一个元素,即存在
第2步,基于上述方法选择的采用如下公式获得新的候选超参数集合,
其中,和表示新的候选超参数集合,分别对应核参数向量和惩罚参数向量;Nker和Nreg表示依据经验设定的新的超参数的数量;和为依据经验设定的超参数收缩和扩放因子,其默认值均为10,
通过再次采用网格搜索策略自适应获得第ith个子模型的超参数对
对全部子系统执行上述过程,子模型预测输出的集合可表示为,
其中,fi(·)表示第ith个子模型,
结合基于分支定界(BB)的最优化选择算法和基于误差的信息熵加权算法对上述子模型进行自适应的优化选择和计算加权系数。在给定候选子模型和加权算法后,最佳子模型选择与加权类似最优特征选择。面向有限数量的候选子模型,通过多次运行优化和加权算法,可获得集成尺寸为2到(I-1)的SEN模型,最后排序这些优化的SEN模型并将预测性能最佳的作为最终DXN软测量模型,
假定最终DXN软测量模型的集成尺寸为Isel,其预测输出值可由下式计算:
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910224790.4A CN109960873B (zh) | 2019-03-24 | 2019-03-24 | 一种城市固废焚烧过程二噁英排放浓度软测量方法 |
US16/967,408 US12002014B2 (en) | 2019-03-24 | 2019-12-02 | Soft measurement method for dioxin emission concentration in municipal solid waste incineration process |
PCT/CN2019/122326 WO2020192166A1 (zh) | 2019-03-24 | 2019-12-02 | 一种城市固废焚烧过程二噁英排放浓度软测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910224790.4A CN109960873B (zh) | 2019-03-24 | 2019-03-24 | 一种城市固废焚烧过程二噁英排放浓度软测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109960873A true CN109960873A (zh) | 2019-07-02 |
CN109960873B CN109960873B (zh) | 2021-09-10 |
Family
ID=67024801
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910224790.4A Active CN109960873B (zh) | 2019-03-24 | 2019-03-24 | 一种城市固废焚烧过程二噁英排放浓度软测量方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN109960873B (zh) |
WO (1) | WO2020192166A1 (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020192166A1 (zh) * | 2019-03-24 | 2020-10-01 | 北京工业大学 | 一种城市固废焚烧过程二噁英排放浓度软测量方法 |
WO2023083009A1 (zh) * | 2021-11-11 | 2023-05-19 | 北京工业大学 | 一种面向mswi过程的多入多出回路控制半实物仿真实验平台 |
WO2023138140A1 (zh) * | 2022-01-19 | 2023-07-27 | 北京工业大学 | 基于宽度混合森林回归的mswi过程二噁英排放软测量方法 |
WO2023165635A1 (zh) * | 2022-03-04 | 2023-09-07 | 北京工业大学 | 基于残差拟合机制简化深度森林回归的炉排炉mswi过程二噁英排放软测量方法 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113189891B (zh) * | 2020-11-08 | 2023-01-13 | 北京工业大学 | 基于双向安全隔离的城市固废焚烧过程半物理仿真平台 |
CN112541242B (zh) * | 2020-12-23 | 2024-01-12 | 燕山大学 | 一种弹性提升城市电与气综合能源系统扩展规划方法 |
CN114661754B (zh) * | 2022-01-20 | 2024-05-03 | 福建中兴电子科技有限公司 | 一种基于分数引导正则化网络的水质污染无监督预警方法 |
CN114518737B (zh) * | 2022-02-28 | 2024-04-16 | 北京工业大学 | 面向mswi过程多模态数据驱动预测模型的仿实时验证平台 |
CN115358152B (zh) * | 2022-08-26 | 2023-04-07 | 绍兴市再生能源发展有限公司 | 垃圾焚烧气体控制与反馈调节系统和方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105910111A (zh) * | 2016-06-15 | 2016-08-31 | 重庆盎瑞悦科技有限公司 | 一种具有废气二次燃烧干燥点火装置的固废焚烧系统及处理方法 |
CN107944173A (zh) * | 2017-12-05 | 2018-04-20 | 北京工业大学 | 一种基于选择性集成最小二乘支撑向量机的二噁英软测量系统 |
CN108062566A (zh) * | 2017-12-13 | 2018-05-22 | 北京工业大学 | 一种基于多核潜在特征提取的智能集成软测量方法 |
CN108549792A (zh) * | 2018-06-27 | 2018-09-18 | 北京工业大学 | 一种基于潜结构映射算法的固废焚烧过程二噁英排放浓度软测量方法 |
US20180311674A1 (en) * | 2017-04-26 | 2018-11-01 | Adr Technology B.V. | Method and Apparatus for Liberating Particles from Moist MSWI Ash |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8002871B2 (en) * | 2008-02-01 | 2011-08-23 | Honeywell International Inc. | Methods and apparatus for an oxygen furnace quality control system |
CN103455635A (zh) * | 2013-09-24 | 2013-12-18 | 华北电力大学 | 基于最小二乘支持向量机集成的热工过程软测量建模方法 |
CN109960873B (zh) * | 2019-03-24 | 2021-09-10 | 北京工业大学 | 一种城市固废焚烧过程二噁英排放浓度软测量方法 |
-
2019
- 2019-03-24 CN CN201910224790.4A patent/CN109960873B/zh active Active
- 2019-12-02 WO PCT/CN2019/122326 patent/WO2020192166A1/zh active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105910111A (zh) * | 2016-06-15 | 2016-08-31 | 重庆盎瑞悦科技有限公司 | 一种具有废气二次燃烧干燥点火装置的固废焚烧系统及处理方法 |
US20180311674A1 (en) * | 2017-04-26 | 2018-11-01 | Adr Technology B.V. | Method and Apparatus for Liberating Particles from Moist MSWI Ash |
CN107944173A (zh) * | 2017-12-05 | 2018-04-20 | 北京工业大学 | 一种基于选择性集成最小二乘支撑向量机的二噁英软测量系统 |
CN108062566A (zh) * | 2017-12-13 | 2018-05-22 | 北京工业大学 | 一种基于多核潜在特征提取的智能集成软测量方法 |
CN108549792A (zh) * | 2018-06-27 | 2018-09-18 | 北京工业大学 | 一种基于潜结构映射算法的固废焚烧过程二噁英排放浓度软测量方法 |
Non-Patent Citations (6)
Title |
---|
JIAN TANG等: "Dioxin Emission Concentration Soft Measuring Method Based on Selective Ensemble Least Square Support Vector Machine Algorithm", 《PROCEEDINGS OF THE 37TH CHINESE CONTROL CONFERENCE》 * |
汤健等: "《基于频谱数据驱动的旋转机械设备负荷软测量》", 30 June 2015, 国防工业出版社 * |
汤健等: "基于主元分析和互信息维数约简策略的网络人侵异常检测", 《NETINFO SECURITY》 * |
汤健等: "基于虚拟样本生成技术的多组分机械信号建模", 《自动化学报》 * |
汤健等: "基于选择性集成核学习算法的固废焚烧过程二英排放浓度软测量", 《化工学报》 * |
龚文: "基于SVM&CBR的E-Learning情绪教学研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020192166A1 (zh) * | 2019-03-24 | 2020-10-01 | 北京工业大学 | 一种城市固废焚烧过程二噁英排放浓度软测量方法 |
WO2023083009A1 (zh) * | 2021-11-11 | 2023-05-19 | 北京工业大学 | 一种面向mswi过程的多入多出回路控制半实物仿真实验平台 |
WO2023138140A1 (zh) * | 2022-01-19 | 2023-07-27 | 北京工业大学 | 基于宽度混合森林回归的mswi过程二噁英排放软测量方法 |
WO2023165635A1 (zh) * | 2022-03-04 | 2023-09-07 | 北京工业大学 | 基于残差拟合机制简化深度森林回归的炉排炉mswi过程二噁英排放软测量方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2020192166A1 (zh) | 2020-10-01 |
CN109960873B (zh) | 2021-09-10 |
US20210233039A1 (en) | 2021-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109960873A (zh) | 一种城市固废焚烧过程二噁英排放浓度软测量方法 | |
Rahat et al. | Data-driven multi-objective optimisation of coal-fired boiler combustion systems | |
Rasheed et al. | Life cycle assessment of a cleaner supercritical coal-fired power plant | |
CN108549792B (zh) | 一种基于潜结构映射算法的固废焚烧过程二噁英排放浓度软测量方法 | |
Nikula et al. | Data-driven framework for boiler performance monitoring | |
Sahu et al. | Targeting for optimal grid-wide deployment of carbon capture and storage (CCS) technology | |
Dirik | Prediction of NOx emissions from gas turbines of a combined cycle power plant using an ANFIS model optimized by GA | |
CN107944173B (zh) | 一种基于选择性集成最小二乘支撑向量机的二噁英软测量系统 | |
Seijo et al. | Modeling and multi-objective optimization of a complex CHP process | |
CN111461355A (zh) | 基于随机森林的二噁英排放浓度迁移学习预测方法 | |
CN110135057A (zh) | 基于多层特征选择的固废焚烧过程二噁英排放浓度软测量方法 | |
NO329798B1 (no) | System og fremgangsmate for empirisk ensemblebasert virtuell sensing av svevestovpartikler | |
Dindarloo et al. | Prediction of the unburned carbon content of fly ash in coal-fired power plants | |
US8019446B2 (en) | Control loop for regulating a combustion process | |
Braun et al. | A neuro-genetic approach for modeling and optimizing a complex cogeneration process | |
CN109978011A (zh) | 一种城市固废焚烧过程二噁英排放浓度预测系统 | |
US11475186B2 (en) | Apparatus and method for deriving boiler combustion model | |
Wen et al. | Dioxin emission prediction from a full-scale municipal solid waste incinerator: Deep learning model in time-series input | |
Weimann et al. | A thermodynamic-based mixed-integer linear model of post-combustion carbon capture for reliable use in energy system optimisation | |
Su et al. | Temporal validation of life cycle greenhouse gas emissions of energy systems in China | |
Chongwatpol et al. | Applying analytics in the energy industry: A case study of heat rate and opacity prediction in a coal-fired power plant | |
Desell et al. | Long term predictions of coal fired power plant data using evolved recurrent neural networks | |
Tang et al. | Dioxin Emission Concentration Soft Measuring Method Based on Selective Ensemble Least Square Support Vector Machine Algorithm | |
Chen et al. | Gas emission prediction for environmental sustainability via heterogeneous data sources correlation with support vector regression | |
Tang et al. | Dioxin soft measuring method in municipal solid waste incineration based on virtual sample generation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |