CN105443212B - 一种基于观测器的单传感器双闭环urea‑SCR反馈控制方法 - Google Patents
一种基于观测器的单传感器双闭环urea‑SCR反馈控制方法 Download PDFInfo
- Publication number
- CN105443212B CN105443212B CN201510497376.2A CN201510497376A CN105443212B CN 105443212 B CN105443212 B CN 105443212B CN 201510497376 A CN201510497376 A CN 201510497376A CN 105443212 B CN105443212 B CN 105443212B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtr
- mtd
- mover
- 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.)
- Expired - Fee Related
Links
- XSQUKJJJFZCRTK-UHFFFAOYSA-N Urea Chemical compound NC(N)=O XSQUKJJJFZCRTK-UHFFFAOYSA-N 0.000 title claims abstract description 32
- 239000004202 carbamide Substances 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000002347 injection Methods 0.000 claims abstract description 18
- 239000007924 injection Substances 0.000 claims abstract description 18
- 230000007246 mechanism Effects 0.000 claims abstract description 9
- 238000012360 testing method Methods 0.000 claims abstract description 5
- 238000013461 design Methods 0.000 claims abstract description 4
- 238000005259 measurement Methods 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 3
- 230000001568 sexual effect Effects 0.000 claims 1
- 210000002700 urine Anatomy 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 6
- 230000035945 sensitivity Effects 0.000 abstract description 4
- 230000008901 benefit Effects 0.000 abstract description 3
- 238000011217 control strategy Methods 0.000 abstract description 3
- 238000011161 development Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 abstract description 3
- 238000001914 filtration Methods 0.000 abstract description 3
- 230000009977 dual effect Effects 0.000 abstract description 2
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 description 31
- 229910021529 ammonia Inorganic materials 0.000 description 13
- 238000006243 chemical reaction Methods 0.000 description 10
- 229910000069 nitrogen hydride Inorganic materials 0.000 description 5
- 239000003054 catalyst Substances 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000007789 gas Substances 0.000 description 3
- 238000001179 sorption measurement Methods 0.000 description 3
- 239000007921 spray Substances 0.000 description 3
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 238000006555 catalytic reaction Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000006722 reduction reaction Methods 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 239000005864 Sulphur Substances 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000010531 catalytic reduction reaction Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000003795 desorption Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 238000013525 hypothesis research Methods 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005507 spraying Methods 0.000 description 1
- 229910052717 sulfur Inorganic materials 0.000 description 1
- 239000011593 sulfur Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Classifications
-
- 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/40—Engine management systems
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Feedback Control In General (AREA)
- Exhaust Gas After Treatment (AREA)
Abstract
本发明公开了一种基于观测器的单传感器双闭环urea‑SCR反馈控制方法,根据NOx传感器交叉敏感性因子与排气温度的关系map标定试验,基于urea‑SCR系统数学模型,应用非线性卡尔曼滤波观测理论,设计了基于数据/机理的柴油机urea‑SCR系统和参数观测器。以上述两个估计量与各自目标值的差值作为控制器输入,设计了尿素喷射双闭环反馈控制器。本发明主要以该观测器/控制器以及安装在urea‑SCR系统排气端的NOx传感器为主要部件,构造了尿素喷射反馈控制系统;本发明节约硬件成本,双闭环系统可以在单控制输入且权限受限、双输出约束且互为矛盾条件下,达到折中控制效果;观测器设计简单合理,数据/机理观测器可以减少发动机电控单元开发过程中过多的标定过程,控制策略灵活性及通用性。
Description
技术领域
本发明涉及一种柴油机尾气后处理电子反馈控制方法,特别涉及一种基于观测器的单传感器双闭环urea-SCR反馈控制方法。
背景技术
与汽油机相比,柴油机具有更高的燃油经济性和更高的动力输出。然而,柴油机的排放问题对其发展一直是个挑战,特别是NOx排放。研究人员提出了很多降低NOx排放的方法,其中urea-SCR即尿素选择性催化还原排放后处理技术,凭借其低成本、高效率等特点具备一定优势。而且urea-SCR系统对硫的敏感性较低,根据目前我国燃油中硫含量较高的实际国情,该技术更具优势。urea-SCR系统工作的基本原理是利用氨(NH3)还原发动机产生的NOx,这些氨一般都来源于32.5%的尿素溶液(添蓝溶液)。虽然氨能够还原NOx,但其本身也是对人体有害的,并且有着刺鼻的气味,尽管在当前的排放法规中没有明文规定。所以,同时实现较低的NOx排放与氨逃逸成为了urea-SCR系统研究的主要挑战。
针对上述问题,有学者提出了一些硬件改进措施。但在实际应用中,大多数都是以增加成本和所占空间为代价。目前,人们已经普遍认识到,通过改进尿素喷射技术达到上述控制目标,是一种较便捷且经济的方法。当前的一些尿素喷射控制方法,主要集中在反馈控制设计上。其中,一部分反馈控制是基于NOx传感器的。但NOx传感器对氨有交叉敏感性,导致测量信号存在误差。而且,只取NOx反馈信号也难以保证氨排放满足标准。虽然德尔福公司最近推出了NH3传感器,但由于其还未被量产,在实际的工业应用中很少用到。此外,常规的控制算法多数不是基于模型的,开发过程需要大量的标定试验,并且对于参数不确定性等因素抗干扰能力差。
因此,当前urea-SCR系统反馈控制依然面临以下挑战:
1、单控制输入(尿素喷射)、双受限输出(并且NOx与氨排放量互为矛盾);
2、执行器的控制权限受限。即控制器只能增加尿素喷射量,不能减少;
3、传感器测量反馈控制方法存在困难。即NOx传感器测量对氨有交叉敏感性,NH3传感器尚未投入量产;
4、控制策略的易实施性,以及工业应用成本问题;
5、以标定为主的电子控制算法面临耗时耗力、受人为主观因素影响大以及抗干扰能力差等问题。
发明内容
本发明的目的是为了解决现有NOx传感器测量的信号存在误差,以及同时实现较低的NOx排放与氨逃逸的问题,而提供的一种基于观测器的单传感器双闭环urea-SCR反馈控制方法。
本发明根据NOx传感器交叉敏感性因子与排气温度的关系map标定试验,基于urea-SCR系统数学模型,应用非线性卡尔曼滤波观测理论,设计了基于数据/机理的和参数观测器。以上述两个估计量与各自目标值的差值作为控制器输入,设计了尿素喷射双闭环反馈控制器。本发明主要以该观测器/控制器以及安装在urea-SCR系统排气端的NOx传感器为主要部件,构造了尿素喷射反馈控制系统。该方法具体包括以下步骤:
一、采集NOx传感器的信号,传送给观测器/控制器模块;
二、观测器/控制器模块依据NOx传感器的信号,估计出和参数值;
三、观测器/控制器模块利用上述两个估计量与各自目标值的差值作为反馈信号,设计了尿素喷射双闭环反馈控制系统。其中,以信号与其排放指标的差值作为主误差反馈量,以信号与其排放指标的差值作为副误差反馈量;
四、观测器/控制器模块通过CAN总线,将当前状态下计算得到的尿素喷射量传递给尿素计量喷射控制器;
五、尿素计量喷射控制器将信号传递给尿素喷头完成喷射任务。
其中,和数据/机理观测器具体的实现步骤如下:
(1)、建立urea-SCR系统数学模型
其中:
(2)、简化和观测器模型
定义状态变量时变参数控制输入变量得到:
z=x1+Kcsx2, (II)
其中:
进一步整理为连续状态方程和离散观测方程,以及模型噪声W和测量噪声V构成的非线性系统:
(3)、基于urea-SCR系统数学模型设计估计算法并整理形式:
其中,(-)表示测量更新还没有发生,(+)表示测量更新已经发生,tk表示当前时刻,k=0,1…N。
(4)、交叉敏感性因子Kcs与排气温度T之间的关系map标定
根据交叉敏感性因子Kcs与排气温度T之间的数学关系Kcs=eaT,通过全工况的试验数据拟合得到常数a。
上述公式中的常量命名见下表1:
表1常量命名法
上述公式中的变量命名见下表2:
表2变量命名法
本发明的有益效果:
1、节约硬件成本。该方法只用了一个NOx传感器、一个数据/机理观测器单元,即可实现NOx和氨双闭环反馈控制;
2、双闭环系统可以在单控制输入且权限受限、双输出约束且互为矛盾条件下,达到折中控制效果;
3、观测器设计简单合理。利用NOx传感器信号、交叉敏感性因子与排气温度的关系map,即可观测出NOx和氨估计值;
4、数据/机理观测器可以减少发动机电控单元开发过程中过多的标定过程,而且基于机理模型的设计方法可以增强观测器算法的鲁棒性,对于系统外部扰动有较好的抑制作用;
5、控制策略灵活性及通用性。该方法可以根据NOx和氨排放法规的要求调整排放控制目标。
附图说明
图1为本发明实施例增压柴油机urea-SCR系统尿素喷射反馈控制结构示意图。
图2为本发明实施例基于数据/机理观测器的单NOx传感器双闭环反馈控制原理框图。
图3为本发明实施例混合扩展卡尔曼滤波算法的流程图。
具体实施方式
请参阅图1、图2和图3所示,为本发明的实施例,本实施例根据NOx传感器交叉敏感性因子与排气温度的关系map标定试验,基于urea-SCR系统数学模型,应用非线性卡尔曼滤波观测理论,设计了基于数据/机理的和参数观测器。以上述两个估计量与各自目标值的差值作为控制器输入,设计了尿素喷射双闭环反馈控制器。参阅图1,本发明主要以该观测器/控制器以及安装在urea-SCR系统排气端的NOx传感器为主要部件,构造了尿素喷射反馈控制系统。该方法具体包括以下步骤:
一、采集NOx传感器3的信号,传送给观测器/控制器模块5;
二、观测器/控制器模块5依据NOx传感器3的信号,估计出和参数值;
三、观测器/控制器模块5利用上述两个估计量与各自目标值的差值作为反馈信号,设计了尿素喷射双闭环反馈控制系统。其中,以信号与其排放指标的差值作为主误差反馈量,以信号与其排放指标的差值作为副误差反馈量;
四、观测器/控制器模块5通过CAN总线6,将当前状态下计算得到的尿素喷射量传递给尿素计量喷射控制器2;
五、尿素计量喷射控制器2将信号传递给尿素喷头1完成喷射任务。
其中,和数据/机理观测器具体的实现步骤如下:
1、建立urea-SCR系统数学模型
urea-SCR系统中的化学反应主要包括下面几个过程。首先,发动机排出的NOx,in和添蓝溶液产生的NH3,in进入系统。在催化器内,部分NH3被吸附在催化剂表面,变成吸附态的氨(NH3(ads)),其余部分仍为气态氨(NH3(g)),并且这两部分可以相互转化,反应方程式如下:
吸附态的氨能够在催化剂的作用下选择性的与NOx反应生成N2和H2O,当温度高于200摄氏度时,该反应主要按照Eley-Rideal机制进行,反应方程式如下所示:
同时,吸附态的氨也能够被氧化成N2和H2O,反应方程式如下:
经过上述几步化学反应,剩余的NOx(NOx,out)和氨(NH3,out)从urea-SCR系统排出。
根据阿列纽斯方程,研究对化学反应过程的反应速率进行建模。方程(1)中的吸附(rads)和解吸附(rdes)速率由如下公式表示:
方程(2)中的还原速率(rSCR)由如下公式表示:
方程(3)对应的反应速率(rox)由如下公式表示:
需要指出的是,在柴油机废气中氧浓度总是很高的,相比于其他气体成分来说十分充足。所以,大多数的urea-SCR反应速率对氧浓度的变化都不敏感,可以忽略。
研究假设上述催化反应是一个连续搅拌釜式反应器(CSTR),即催化系统的所有状态是均质的。基于上述假设,根据质量平衡以及热量平衡原则,urea-SCR系统模型的常微分方程(ODEs)表示如下:
考虑到urea-SCR系统中,气体成分的流量与浓度之间的关系为
并且,将公式(4),(5),(6)以及(8)代入公式(7)中,得到:
其中:
2、简化和观测器模型
为了实现和观测,将看作可测量变量,考虑的状态量仅为和定义状态变量时变参数控制输入变量基于urea-SCR模型(9),得到系统面向控制模型:
其中:
令
(13)
h(x(k),k)=x1+Kcsx2
进一步整理为连续状态方程和离散观测方程,以及模型噪声W和测量噪声V构成的非线性系统:
E[w(t)]=0,
E[v(t)]=0,
E[w(t)wT(τ)]=Q(t)δ(t-τ), (15)
E[v(t)vT(τ)]=R(t)δ(t-τ),
E[w(t)vT(τ)]=0.
其中x∈Rp,z∈Rm,Q(t)为p×p阶系统噪声方差阵,对任意t≥t0为连续且非负定,R(t)为m×m阶观测噪声方差阵,对任意t≥t0为连续且正定。
3、基于urea-SCR系统数学模型设计估计算法并整理形式:
再对观测模型公式(13)应用混合扩展卡尔曼滤波,其算法流程图如图3所示。扩展卡尔曼滤波中状态初值x(0)、初始误差方差阵p(0)的选取关系到观测性能的好坏以及数值稳定性,模型噪声W和测量噪声V都假设是高斯白噪声。并且,雅可比矩阵:得到的观测器形式如下:
其中,(-)表示测量更新还没有发生,(+)表示测量更新已经发生,tk表示当前时刻,k=0,1…N。
4、交叉敏感性因子Kcs与排气温度T之间的关系map标定
根据交叉敏感性因子Kcs与排气温度T之间的数学关系Kcs=eaT,通过全工况的试验数据拟合得到常数a。
Claims (1)
1.一种基于观测器的单传感器双闭环urea-SCR反馈控制方法,其特征在于:该方法具体包括以下步骤:
一、采集NOx传感器的信号,传送给观测器/控制器模块;
二、观测器/控制器模块依据NOx传感器的信号,估计出和参数值;
三、观测器/控制器模块利用上述两个估计量与各自目标值的差值作为反馈信号,设计了尿素喷射双闭环反馈控制系统,其中,以信号与其排放指标的差值作为主误差反馈量,以信号与其排放指标的差值作为副误差反馈量;
四、观测器/控制器模块通过CAN总线,将当前状态下计算得到的尿素喷射量传递给尿素计量喷射控制器;
五、尿素计量喷射控制器将信号传递给尿素喷头完成喷射任务;
其中,和数据/机理观测器具体的实现步骤如下:
(一)、建立urea-SCR系统数学模型:
<mrow>
<mover>
<mi>T</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>6</mn>
</msub>
<msubsup>
<mi>m</mi>
<mrow>
<mi>E</mi>
<mi>G</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>a</mi>
<mn>7</mn>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>T</mi>
<mn>4</mn>
</msup>
<mo>-</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mi>m</mi>
<mi>b</mi>
</mrow>
<mn>4</mn>
</msubsup>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mover>
<mi>C</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<msub>
<mi>NO</mi>
<mi>x</mi>
</msub>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>n</mi>
<mrow>
<msub>
<mi>NO</mi>
<mi>x</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>C</mi>
<mrow>
<msub>
<mi>NO</mi>
<mi>x</mi>
</msub>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>m</mi>
<mrow>
<mi>E</mi>
<mi>G</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mi>T</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
<msub>
<mi>&Theta;</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>C</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>n</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
<mo>,</mo>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>-</mo>
<msub>
<mi>C</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
</mrow>
</msub>
<mo>&lsqb;</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>m</mi>
<mrow>
<mi>E</mi>
<mi>G</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mi>T</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&Theta;</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>&Theta;</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
</mrow>
</msub>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mi>s</mi>
</msub>
<msub>
<mover>
<mi>&Theta;</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&Theta;</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
</mrow>
</msub>
</mrow>
<mo>)</mo>
<msub>
<mi>C</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
</mrow>
</msub>
<mo>-</mo>
<mo>&lsqb;</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>C</mi>
<mrow>
<msub>
<mi>NO</mi>
<mi>x</mi>
</msub>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<msub>
<mi>&Theta;</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
</mrow>
</msub>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>S</mi>
<mo>,</mo>
<mi>E</mi>
<mi>G</mi>
</mrow>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mi>m</mi>
<mi>b</mi>
</mrow>
</msub>
</mfrac>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>&epsiv;V</mi>
<mi>c</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>c</mi>
<mi>s</mi>
</msub>
<msub>
<mi>S</mi>
<mi>c</mi>
</msub>
<msub>
<mi>&alpha;</mi>
<mrow>
<mi>Pr</mi>
<mi>o</mi>
<mi>b</mi>
</mrow>
</msub>
<msqrt>
<mfrac>
<mrow>
<mi>R</mi>
<mi>T</mi>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;M</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
</mrow>
</msub>
</mrow>
</mfrac>
</msqrt>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>c</mi>
<mi>s</mi>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mi>D</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>D</mi>
<mi>e</mi>
<mi>s</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>R</mi>
<mi>T</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>c</mi>
<mi>s</mi>
</msub>
<msub>
<mi>RTk</mi>
<mrow>
<mi>S</mi>
<mi>C</mi>
<mi>R</mi>
</mrow>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>S</mi>
<mi>C</mi>
<mi>R</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>R</mi>
<mi>T</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>c</mi>
<mi>s</mi>
</msub>
<msub>
<mi>k</mi>
<mrow>
<mi>O</mi>
<mi>x</mi>
</mrow>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>O</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mi>R</mi>
<mi>T</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</msup>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>6</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>c</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>E</mi>
<mi>G</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>c</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<msub>
<mi>m</mi>
<mi>c</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
<msub>
<mi>a</mi>
<mn>7</mn>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&epsiv;</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>d</mi>
<mo>,</mo>
<mi>s</mi>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
<msub>
<mi>&sigma;</mi>
<mrow>
<mi>s</mi>
<mi>b</mi>
</mrow>
</msub>
<msub>
<mi>A</mi>
<mrow>
<mi>r</mi>
<mi>a</mi>
<mi>d</mi>
<mo>,</mo>
<mi>s</mi>
<mi>c</mi>
<mi>r</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>c</mi>
<mrow>
<mi>p</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<msub>
<mi>m</mi>
<mi>c</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
上述公式中的常量命名法见下表1,变量命名法见下表2:
表1 常量命名法
表2 变量命名法
(二)、简化和观测器模型:
定义状态变量时变参数控制输入变量得到:
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
</mrow>
<mrow>
<msub>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mi>B</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
</mrow>
z=x1+Kcsx2, (3)
其中:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>m</mi>
<mrow>
<mi>E</mi>
<mi>G</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mi>T</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>&Theta;</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
</mrow>
</msub>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msubsup>
<mi>m</mi>
<mrow>
<mi>E</mi>
<mi>G</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mi>T</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&Theta;</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>B</mi>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<msub>
<mi>&Theta;</mi>
<mrow>
<msub>
<mi>NH</mi>
<mn>3</mn>
</msub>
</mrow>
</msub>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
令
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>,</mo>
<mi>u</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<msub>
<mi>A</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mi>B</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>K</mi>
<mrow>
<mi>c</mi>
<mi>s</mi>
</mrow>
</msub>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
进一步整理为连续状态方程和离散观测方程,以及模型噪声W和测量噪声V构成的非线性系统:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>,</mo>
<mi>u</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>w</mi>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>z</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>v</mi>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>E</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>E</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>E</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>w</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mo>=</mo>
<mi>Q</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>E</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>v</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mo>=</mo>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>&delta;</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>E</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mi>w</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>v</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中x∈Rp,z∈Rm,Q(t)为p×p阶系统噪声方差阵,对任意t≥t0为连续且非负定,R(t)为m×m阶观测噪声方差阵,对任意t≥t0为连续且正定;
(三)、基于urea-SCR系统数学模型设计估计算法并整理形式:
再对观测模型公式(5)应用混合扩展卡尔曼滤波,扩展卡尔曼滤波中状态初值x(0)、初始误差方差阵p(0)的选取关系到观测性能的好坏以及数值稳定性,模型噪声W和测量噪声V都假设是高斯白噪声;并且,雅可比矩阵:得到的观测器形式如下:
<mrow>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>(</mo>
<mo>+</mo>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>(</mo>
<mo>-</mo>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>K</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>{</mo>
<mi>z</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>(</mo>
<mrow>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mrow>
<mo>(</mo>
<mo>-</mo>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>,</mo>
<msub>
<mi>t</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mo>,</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,(-)表示测量更新还没有发生,(+)表示测量更新已经发生,tk表示当前时刻,k=0,1…N;
(四)、交叉敏感性因子Kcs与排气温度T之间的关系map标定:
根据交叉敏感性因子Kcs与排气温度T之间的数学关系Kcs=eaT,通过全工况的试验数据拟合得到常数a。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510497376.2A CN105443212B (zh) | 2015-11-24 | 2015-11-24 | 一种基于观测器的单传感器双闭环urea‑SCR反馈控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510497376.2A CN105443212B (zh) | 2015-11-24 | 2015-11-24 | 一种基于观测器的单传感器双闭环urea‑SCR反馈控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105443212A CN105443212A (zh) | 2016-03-30 |
CN105443212B true CN105443212B (zh) | 2018-02-13 |
Family
ID=55553781
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510497376.2A Expired - Fee Related CN105443212B (zh) | 2015-11-24 | 2015-11-24 | 一种基于观测器的单传感器双闭环urea‑SCR反馈控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105443212B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106014571B (zh) * | 2016-05-31 | 2018-10-19 | 无锡威孚力达催化净化器有限责任公司 | 发动机原机NOx值的计算方法 |
CN106437956B (zh) * | 2016-09-23 | 2018-10-12 | 上海海事大学 | 一种串联选择性催化还原系统的模糊控制方法 |
CN106368772B (zh) * | 2016-11-04 | 2019-03-01 | 吉林大学 | Scr系统尿素喷射控制方法 |
CN110262254A (zh) * | 2019-07-15 | 2019-09-20 | 吉林师范大学 | 氨覆盖率观测和跟踪整体设计控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101790623A (zh) * | 2007-08-28 | 2010-07-28 | 戴姆勒股份公司 | Scr排气后处理系统的运行方法及诊断方法 |
CN101907012A (zh) * | 2009-03-06 | 2010-12-08 | 德国曼商用车辆股份公司 | 调整还原剂计量添加量进行选择性催化还原的方法 |
CN104314650A (zh) * | 2014-07-29 | 2015-01-28 | 潍柴动力股份有限公司 | 一种scr前馈控制方法和装置 |
CN104405481A (zh) * | 2014-09-29 | 2015-03-11 | 同济大学 | 一种尿素喷射剂量标定方法 |
CN104632323A (zh) * | 2014-12-18 | 2015-05-20 | 吉林大学 | 一种尿素scr氨覆盖率反馈跟踪控制方法 |
-
2015
- 2015-11-24 CN CN201510497376.2A patent/CN105443212B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101790623A (zh) * | 2007-08-28 | 2010-07-28 | 戴姆勒股份公司 | Scr排气后处理系统的运行方法及诊断方法 |
CN101907012A (zh) * | 2009-03-06 | 2010-12-08 | 德国曼商用车辆股份公司 | 调整还原剂计量添加量进行选择性催化还原的方法 |
CN104314650A (zh) * | 2014-07-29 | 2015-01-28 | 潍柴动力股份有限公司 | 一种scr前馈控制方法和装置 |
CN104405481A (zh) * | 2014-09-29 | 2015-03-11 | 同济大学 | 一种尿素喷射剂量标定方法 |
CN104632323A (zh) * | 2014-12-18 | 2015-05-20 | 吉林大学 | 一种尿素scr氨覆盖率反馈跟踪控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105443212A (zh) | 2016-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105443212B (zh) | 一种基于观测器的单传感器双闭环urea‑SCR反馈控制方法 | |
Zhang et al. | Adaptive sliding-mode observer design for a selective catalytic reduction system of ground-vehicle diesel engines | |
Devarakonda et al. | Model-based estimation and control system development in a urea-SCR aftertreatment system | |
CN108087065B (zh) | Scr催化器氨浓度修正方法及scr处理系统 | |
CN103527293B (zh) | 一种尿素喷射控制方法和控制单元 | |
Zhang et al. | Optimal dosing and sizing optimization for a ground-vehicle diesel-engine two-cell selective catalytic reduction system | |
CN108915827B (zh) | 基于SCR化学反应数学模型改善发动机NOx排放量的方法 | |
CN101832167B (zh) | 一种scr催化器中的氨气吸附控制方法 | |
CN201329254Y (zh) | 以排气中NOx为变量的车载SCR计量喷射系统 | |
CN102588055A (zh) | 一种独立式scr计量喷射控制方法及系统 | |
Zhang et al. | Cycle-based optimal NOx emission control of selective catalytic reduction systems with dynamic programming algorithm | |
Zhang et al. | Integrated study of inland-vessel diesel engine two-cell SCR systems with dynamic references | |
CN108595802A (zh) | 基于仿真的scr系统尿素结晶风险评价方法 | |
CN102678239A (zh) | 基于NOx传感器独立的SCR系统控制方法及控制装置 | |
CN109404108A (zh) | 一种用于柴油机NOx排放值的计算方法及装置 | |
CN106368772B (zh) | Scr系统尿素喷射控制方法 | |
CN104405481A (zh) | 一种尿素喷射剂量标定方法 | |
CN106762049A (zh) | 基于nmpc的两核尿素scr系统排放控制方法 | |
DE102011087082B4 (de) | Verfahren zum Betreiben eines SCRF-Katalysatorsystems | |
Ma et al. | Observer-based estimation of aging condition for selective catalytic reduction systems in vehicle applications | |
CN106527143B (zh) | 基于数据驱动预测控制的scr系统尿素喷射控制方法 | |
Zhou et al. | State estimation in the automotive SCR DeNOx process | |
Kang et al. | Research on cross sensitivity of NOx sensor and Adblue injection volume in accordance with the actual situation based on cubature Kalman filter | |
US8808630B2 (en) | Exhaust gas catalytic converter system and method for operating an exhaust gas catalytic converter | |
Zhao et al. | An ammonia coverage ratio observing and tracking controller: stability analysis and simulation evaluation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180213 |