CN106368772A - Scr***尿素喷射控制方法 - Google Patents
Scr***尿素喷射控制方法 Download PDFInfo
- Publication number
- CN106368772A CN106368772A CN201610962526.7A CN201610962526A CN106368772A CN 106368772 A CN106368772 A CN 106368772A CN 201610962526 A CN201610962526 A CN 201610962526A CN 106368772 A CN106368772 A CN 106368772A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- theta
- sigma
- phi
- scr
- 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
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F01—MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
- F01N—GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR MACHINES OR ENGINES IN GENERAL; GAS-FLOW SILENCERS OR EXHAUST APPARATUS FOR INTERNAL COMBUSTION ENGINES
- F01N3/00—Exhaust or silencing apparatus having means for purifying, rendering innocuous, or otherwise treating exhaust
- F01N3/08—Exhaust or silencing apparatus having means for purifying, rendering innocuous, or otherwise treating exhaust for rendering innocuous
- F01N3/10—Exhaust or silencing apparatus having means for purifying, rendering innocuous, or otherwise treating exhaust for rendering innocuous by thermal or catalytic conversion of noxious components of exhaust
- F01N3/18—Exhaust or silencing apparatus having means for purifying, rendering innocuous, or otherwise treating exhaust for rendering innocuous by thermal or catalytic conversion of noxious components of exhaust characterised by methods of operation; Control
- F01N3/20—Exhaust or silencing apparatus having means for purifying, rendering innocuous, or otherwise treating exhaust for rendering innocuous by thermal or catalytic conversion of noxious components of exhaust characterised by methods of operation; Control specially adapted for catalytic conversion ; Methods of operation or control of catalytic converters
- F01N3/2066—Selective catalytic reduction [SCR]
- F01N3/208—Control of selective catalytic reduction [SCR], e.g. dosing of reducing agent
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- 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
-
- 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
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/12—Improving ICE efficiencies
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Computer Hardware Design (AREA)
- Toxicology (AREA)
- Health & Medical Sciences (AREA)
- Mechanical Engineering (AREA)
- Evolutionary Computation (AREA)
- Combustion & Propulsion (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Exhaust Gas After Treatment (AREA)
- Feedback Control In General (AREA)
Abstract
一种SCR***尿素喷射控制方法,属于柴油机控制技术领域。本发明的目的是从尿素SCR***化学反应机理出发,研究***的偏微分方程模型(PDE),选择合适的逼近方法简化模型,利用滚动优化控制处理时滞和不确定性优势的SCR***尿素喷射控制方法。本发明的步骤:a、尿素SCR***偏微分建模,b、偏微分***预测控制器设计。本发明基于实际的尿素SCR***控制需求,将分布参数***建模与预测控制器设计结合在一起,能够在精确描述***动态性能的同时,解决该***的约束条件下的非线性优化问题。
Description
技术领域
本发明属于柴油机控制技术领域。
背景技术
与汽油机相比,柴油机具有更高的燃油经济性和动力输出。然而,由于其稀燃的特性,柴油机会比汽油机产生更多的NOx有害气体。随着全世界范围内针对NOx排放越来越严格的法规出台,多种为降低NOx的排放后处理***面世了。这些技术就包括尿素选择性催化还原(SCR)***。尿素SCR***工作时则不需额外燃油,而且尿素消耗也相对较低,凭借这些优势,已经在汽车工业界占据一定优势。在我国,目前的实际国情是燃油中硫含量较高,而且许多种排放控制技术推广都受到限制。所以,凭借其对硫的敏感性较低的特性,尿素SCR排放后处理技术在我国的发展更具优势。
尿素SCR技术的基本原理是利用NOx与氨(NH3)之间的氧化还原反应,而所用的氨一般都来源于32.5%的尿素溶液(添蓝溶液)。虽然氨能够还原NOx,但其较高的排放也是对人体有害的,并且有着刺鼻的气味。为实现较高的NOx转化效率,要有充分的氨做为还原剂;但是,这一点反过来会增加氨的逃逸量,这一矛盾成为了尿素SCR***研究面临的主要挑战之一。目前更为普遍的共识是,通过改进尿素喷射控制技术达到上述目标,是一种较便捷且经济的方法。针对尿素SCR控制***的研究,当前较为先进的是一些基于模型的反馈、前馈控制方法。但是,这些模型大多是基于常微分方程(ODE)建立的。SCR***作为一个典型的多物质沿催化器轴向不均匀传播、并发生复杂化学反应的分布参数***,不但具有无限维特性,复杂的多物质轴向传播的环境也使***具有较大的滞后和不确定性,上述基于集中参数的ODE模型的控制方法很难达到理想的控制效果。此外,对于尿素SCR***参数来说,也存在较多的约束条件。例如,***输入(尿素喷射器喷射量)有最大值限制,***输出(NOx与NH3)受到排放法规的限制等。这些带约束的优化控制问题,适合于在模型预测控制的框架内解决。
发明内容
本发明的目的是从尿素SCR***化学反应机理出发,研究***的偏微分方程模型(PDE),选择合适的逼近方法简化模型,利用滚动优化控制处理时滞和不确定性优势的SCR***尿素喷射控制方法。
本发明的步骤:
a、尿素SCR***偏微分建模
基于尿素SCR***各主要物质沿轴向传输特性,主要参考了动态传输反应模型:
磁通密度公式:
Jk=θzCk, (2)
以及一种物质浓度在催化器表面的动态及稳态假设得到的:
建立得到***的一阶双曲型PDE模型:
其中,
进一步整理得到:
其中,
表1和表2分别显示了模型中所有常量和变量的相关定义及参数名义参考值:
表1常量命名法
表2变量命名法
将T,CNOx等看做可测量变量,考虑的状态量仅为和定义状态变量时变参数控制输入变量被控输出为得到***面向控制模型:
其中,
进一步整理成状态空间形式:
上述一阶双曲型***表示为以下状态空间形式:
其中,
b、偏微分***预测控制器设计:
选取Haar小波尺度函数φ(x),将时空变量展开如下形式:
一阶空间偏导数表示为:
采用Eular前向差分法对时间偏导数进行离散:
应用上述方法,同时将时空偏导数替换得到:
对上式进行整理,得到状态变量(zn,i+1|i)的预测值为:
由y=C·X得P步向前预测输出值:
为了快速而平滑地得到***设定值,定义如下参考轨迹:
y*(zn,k+i|k) (19)
预测控制输出为:
定义***误差函数为:
本发明基于实际的尿素SCR***控制需求,将分布参数***建模与预测控制器设计结合在一起,能够在精确描述***动态性能的同时,解决该***的约束条件下的非线性优化问题。其优点:
1、本发明提出的偏微分建模方法能够准确描述SCR***动态。
2、本发明基于实际的尾气排放处理***,将偏微分建模过程与模型预测控制器推导相结合。
3、本发明提出的模型预测控制器,能够处理排放约束的优化控制问题。
附图说明
图1是基于分布式参数建模及预测控制的SCR***尿素喷射控制器设计技术路线图。
具体实施方式
本发明主要技术路线包括:根据尿素SCR***控制需求与控制问题的提炼与描述,建立偏微分模型;参照化学反应机理以及数学逼近方法,推导出偏微分***预测控制器;对优化问题进行归纳,并给出优化问题的求解方式;最后,得到满足排放法规要求的控制优化算法。其的步骤是:
a、尿素SCR***偏微分建模
基于尿素SCR***各主要物质沿轴向传输特性,主要参考了动态传输反应模型:
磁通密度公式:
Jk=θzCk, (2)
以及一种物质浓度在催化器表面的动态及稳态假设得到的:
建立得到***的一阶双曲型PDE模型:
其中,
进一步整理得到:
其中,
表1和表2分别显示了模型中所有常量和变量的相关定义及参数名义参考值:
表1常量命名法
表2变量命名法
为了实现氨覆盖率控制,将T,CNOx等看做可测量变量,考虑的状态量仅为和定义状态变量时变参数控制输入变量被控输出为得到***面向控制模型:
其中,
进一步整理成状态空间形式:
上述一阶双曲型***表示为以下状态空间形式:
其中,
b、偏微分***预测控制器设计:
应用时空分解理论,选择恰当的空间基函数进行逼近。为了降低计算量,本发明选取简单实用的Haar小波尺度函数φ(x),将时空变量展开如下形式:
一阶空间偏导数表示为:
采用Eular前向差分法对时间偏导数进行离散:
应用上述方法,同时将时空偏导数替换得到:
对上式进行整理,得到状态变量(zn,i+1|i)的预测值为:
由y=C·X得P步向前预测输出值:
为了快速而平滑地得到***设定值,定义如下参考轨迹:
y*(zn,k+i|k) (19)
考虑扰动,预测控制输出为:
定义***误差函数为:
上述预测控制求解,是典型的约束非线性二次规划问题,可以采用粒子群优化方法求解。
随着分布参数***理论与预测控制理论发展,已经有学者将两者相结合,并提出了一些理论方法,例如:基于特征线法的鲁棒预测控制算法、基于小波正交逼近的分布参数***预测控制方法等。有学者已经将一些能够处理真实分布参数***约束问题的预测控制算法,应用到如薄膜沉积过程、化工生产过程以及超深井温度场控制等方面。从尿素SCR***化学反应机理出发,研究***的偏微分方程模型(PDE),选择合适的逼近方法简化模型,利用滚动优化控制处理时滞和不确定性的优势,是得到高精度和较强适应能力的控制策略的有效途径。
根据尿素SCR***控制需求与控制问题的提炼与描述,建立偏微分模型。依据该***化学反应特性,选定工业中的常用一阶双曲型***描述流体流动过程。选择Haar小波作为空间基函数,将分布参数***在配置点上转化为集中参数***。将一阶空间偏导数投影到拟Haar小波基上,采用前向Eular法对时间变量离散,得到***的自回归模型。在处理***边界问题时,采用对称延拓的方法抑制边界效应。选择标准二次优化性能指标,将***转化为一个有约束滚动优化问题,设计基于低阶模型的非线性预测控制器,有效处理***的输入-输出约束。最后,采用了粒子群优化方法求解该约束非线性二次规划问题。同时,复杂多变的运行工况和环境给控制带来了大范围的不确定性和时变因素,这就要求指定的控制算法具有自适应的功能,能够随着工况和环境参数变化而进行自调整。因此,选择了自适应滚动优化控制方法优化排放。
仿真在TESIS公司开发的enDYNA发动机精确建模软件平台上进行。选择瞬态标准测试循环ECE,并将氨覆盖率的理想值设定为连续的阶跃信号,验证本发明提出的偏微分模型预测控制器的效果。仿真结果如下图所示。初始阶段氨覆盖率的理想值为0.6;第60秒时刻,氨覆盖率的理想值阶跃到0.8;在第120秒时刻该值阶跃下降到0.7。尿素喷射器约束为0.003mol/s,执行情况在图的下部位置。结果显示,基于尾气排放SCR后处理***的实际情况,本发明设计的控制器能够完成较好的跟踪控制效果。
Claims (1)
1.一种SCR***尿素喷射控制方法,其特征在于:
a、尿素SCR***偏微分建模
基于尿素SCR***各主要物质沿轴向传输特性,主要参考了动态传输反应模型:
磁通密度公式:
Jk=θzCk, (2)
以及一种物质浓度在催化器表面的动态及稳态假设得到的:
建立得到***的一阶双曲型PDE模型:
其中,
进一步整理得到:
其中,
表1和表2分别显示了模型中所有常量和变量的相关定义及参数名义参考值:
表1常量命名法
表2变量命名法
将T,CNOx等看做可测量变量,考虑的状态量仅为和定义状态变量时变参数控制输入变量被控输出为得到***面向控制模型:
其中,
进一步整理成状态空间形式:
上述一阶双曲型***表示为以下状态空间形式:
其中,
b、偏微分***预测控制器设计:
选取Haar小波尺度函数φ(x),将时空变量展开如下形式:
一阶空间偏导数表示为:
采用Eular前向差分法对时间偏导数进行离散:
应用上述方法,同时将时空偏导数替换得到:
对上式进行整理,得到状态变量(zn,i+1|i)的预测值为:
由y=C·X得P步向前预测输出值:
为了快速而平滑地得到***设定值,定义如下参考轨迹:
y*(zn,k+i|k) (19)
预测控制输出为:
定义***误差函数为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610962526.7A CN106368772B (zh) | 2016-11-04 | 2016-11-04 | Scr***尿素喷射控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610962526.7A CN106368772B (zh) | 2016-11-04 | 2016-11-04 | Scr***尿素喷射控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106368772A true CN106368772A (zh) | 2017-02-01 |
CN106368772B CN106368772B (zh) | 2019-03-01 |
Family
ID=57894579
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610962526.7A Active CN106368772B (zh) | 2016-11-04 | 2016-11-04 | Scr***尿素喷射控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106368772B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106762049A (zh) * | 2017-03-14 | 2017-05-31 | 吉林师范大学 | 基于nmpc的两核尿素scr***排放控制方法 |
CN106841531A (zh) * | 2017-03-22 | 2017-06-13 | 吉林大学 | 基于滑膜控制的氨覆盖率非线性观测器设计方法 |
CN109304087A (zh) * | 2018-10-31 | 2019-02-05 | 华中科技大学 | 一种基于脱硝反应动力学方程的电站scr喷氨控制方法 |
CN109424398A (zh) * | 2017-08-28 | 2019-03-05 | 通用汽车环球科技运作有限责任公司 | 内燃机排气***的排放控制*** |
CN109985523A (zh) * | 2019-04-17 | 2019-07-09 | 江苏科技大学 | 一种基于scr的船舶尾气脱硝控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100050614A1 (en) * | 2008-08-28 | 2010-03-04 | Michael Parmentier | System and method for selective catalytic reduction control |
CN101994553A (zh) * | 2009-08-14 | 2011-03-30 | 万国引擎知识产权有限责任公司 | 尿素喷射控制方法 |
CN102155279A (zh) * | 2011-04-19 | 2011-08-17 | 潍柴动力股份有限公司 | 用于控制柴油发动机的尿素喷射***的设备和方法 |
CN103527293A (zh) * | 2013-10-08 | 2014-01-22 | 潍柴动力股份有限公司 | 一种尿素喷射控制方法和控制单元 |
CN105443212A (zh) * | 2015-11-24 | 2016-03-30 | 吉林师范大学 | 一种基于观测器的单传感器双闭环urea-SCR反馈控制方法 |
-
2016
- 2016-11-04 CN CN201610962526.7A patent/CN106368772B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100050614A1 (en) * | 2008-08-28 | 2010-03-04 | Michael Parmentier | System and method for selective catalytic reduction control |
CN101994553A (zh) * | 2009-08-14 | 2011-03-30 | 万国引擎知识产权有限责任公司 | 尿素喷射控制方法 |
CN102155279A (zh) * | 2011-04-19 | 2011-08-17 | 潍柴动力股份有限公司 | 用于控制柴油发动机的尿素喷射***的设备和方法 |
CN103527293A (zh) * | 2013-10-08 | 2014-01-22 | 潍柴动力股份有限公司 | 一种尿素喷射控制方法和控制单元 |
CN105443212A (zh) * | 2015-11-24 | 2016-03-30 | 吉林师范大学 | 一种基于观测器的单传感器双闭环urea-SCR反馈控制方法 |
Non-Patent Citations (4)
Title |
---|
胡云峰,蒋冰晶,宫洵,赵靖华,陈虹: "柴油机尿素SCR***氨覆盖率跟踪控制器设计", 《农业机械学报》 * |
赵靖华,胡云峰,谭振江,高炳钊: "基于降维观测器的氨覆盖率自抗扰反馈控制器设计", 《汽车技术》 * |
赵靖华,胡云峰,高炳钊,陈虹: "基于尿素选择催化还原***的氨覆盖率非线性降维观测器设计", 《吉林大学学报(工学版)》 * |
赵靖华,陈志刚,胡云峰,陈虹: "基于"三步法"的柴油机urea-SCR***控制设计", 《吉林大学学报(工学版)》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106762049A (zh) * | 2017-03-14 | 2017-05-31 | 吉林师范大学 | 基于nmpc的两核尿素scr***排放控制方法 |
CN106841531A (zh) * | 2017-03-22 | 2017-06-13 | 吉林大学 | 基于滑膜控制的氨覆盖率非线性观测器设计方法 |
CN106841531B (zh) * | 2017-03-22 | 2019-04-23 | 吉林大学 | 基于滑膜控制的氨覆盖率非线性观测器设计方法 |
CN109424398A (zh) * | 2017-08-28 | 2019-03-05 | 通用汽车环球科技运作有限责任公司 | 内燃机排气***的排放控制*** |
CN109304087A (zh) * | 2018-10-31 | 2019-02-05 | 华中科技大学 | 一种基于脱硝反应动力学方程的电站scr喷氨控制方法 |
CN109985523A (zh) * | 2019-04-17 | 2019-07-09 | 江苏科技大学 | 一种基于scr的船舶尾气脱硝控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106368772B (zh) | 2019-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106368772A (zh) | Scr***尿素喷射控制方法 | |
CN106837480B (zh) | 一种基于模型的尿素喷射量控制方法及后处理控制*** | |
Chen et al. | Estimation and adaptive nonlinear model predictive control of selective catalytic reduction systems in automotive applications | |
US9650934B2 (en) | Engine and aftertreatment optimization system | |
US10621291B2 (en) | Approach for aftertreatment system modeling and model identification | |
Cloudt et al. | Integrated Emission Management strategy for cost-optimal engine-aftertreatment operation | |
Faghihi et al. | Development of a neural network model for selective catalytic reduction (SCR) catalytic converter and ammonia dosing optimization using multi objective genetic algorithm | |
Zhang et al. | Optimal dosing and sizing optimization for a ground-vehicle diesel-engine two-cell selective catalytic reduction system | |
CN104632323A (zh) | 一种尿素scr氨覆盖率反馈跟踪控制方法 | |
Meisami-Azad et al. | LPV gain-scheduled control of SCR aftertreatment systems | |
Zambrano et al. | Identification of a discrete-time nonlinear Hammerstein-Wiener model for a selective catalytic reduction system | |
CN106762049B (zh) | 基于nmpc的两核尿素scr***排放控制方法 | |
Meisami-Azad et al. | An adaptive control strategy for urea-SCR aftertreatment system | |
Zhao et al. | Sequential optimization of eco-driving taking into account fuel economy and emissions | |
You et al. | A fuzzy logic urea dosage controller design for two-cell selective catalytic reduction systems | |
CN105443212A (zh) | 一种基于观测器的单传感器双闭环urea-SCR反馈控制方法 | |
Chen et al. | A robust ammonia coverage ratio control method for a two-cell selective catalytic reduction system in low temperature operations | |
CN108868976B (zh) | 一种基于apso的双串联scr***尿素喷射规律云计算方法 | |
CN110244565B (zh) | 一种scr***分区控制方法和装置 | |
Zhang et al. | Multi-objective optimization of Fe-based SCR catalyst on the NOx conversion efficiency for a diesel engine based on FGRA-ANN/RF | |
Tayamon et al. | Control of selective catalytic reduction systems using feedback linearisation | |
CN108229091A (zh) | Scr后处理***化学反应动力学模型的搭建方法 | |
CN108386262B (zh) | 一种柴油机串联scr***氨气覆盖率与存储量观测方法 | |
Lim et al. | Gain-scheduling Selective Catalytic Reduction Control in Diesel Engines with Switched H∞ Controllers | |
Nekooei et al. | Hybrid fuzzy logic controller in Matlab/Simulink for controlling AFR of SI engine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |