CN110273737B - 柴油机后处理系统催化器热失效容错控制方法及其装置 - Google Patents
柴油机后处理系统催化器热失效容错控制方法及其装置 Download PDFInfo
- Publication number
- CN110273737B CN110273737B CN201910285859.4A CN201910285859A CN110273737B CN 110273737 B CN110273737 B CN 110273737B CN 201910285859 A CN201910285859 A CN 201910285859A CN 110273737 B CN110273737 B CN 110273737B
- Authority
- CN
- China
- Prior art keywords
- scr catalyst
- downstream
- concentration
- catalyst
- value
- 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.)
- Active
Links
- 239000003054 catalyst Substances 0.000 title claims abstract description 383
- 238000000034 method Methods 0.000 title claims abstract description 35
- QGZKDVFQNNGYKY-UHFFFAOYSA-N Ammonia Chemical compound N QGZKDVFQNNGYKY-UHFFFAOYSA-N 0.000 claims abstract description 205
- 229910021529 ammonia Inorganic materials 0.000 claims abstract description 102
- 238000011144 upstream manufacturing Methods 0.000 claims abstract description 49
- XSQUKJJJFZCRTK-UHFFFAOYSA-N Urea Chemical compound NC(N)=O XSQUKJJJFZCRTK-UHFFFAOYSA-N 0.000 claims abstract description 46
- 239000004202 carbamide Substances 0.000 claims abstract description 46
- 238000002347 injection Methods 0.000 claims abstract description 43
- 239000007924 injection Substances 0.000 claims abstract description 43
- 238000005259 measurement Methods 0.000 claims abstract description 29
- 238000001914 filtration Methods 0.000 claims abstract description 11
- 238000003745 diagnosis Methods 0.000 claims description 34
- 238000009529 body temperature measurement Methods 0.000 claims description 17
- 238000006243 chemical reaction Methods 0.000 claims description 16
- 229910000069 nitrogen hydride Inorganic materials 0.000 claims description 16
- 230000003044 adaptive effect Effects 0.000 claims description 15
- 230000007246 mechanism Effects 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000010531 catalytic reduction reaction Methods 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims description 3
- 238000003795 desorption Methods 0.000 claims description 3
- 238000007254 oxidation reaction Methods 0.000 claims description 3
- 238000001179 sorption measurement Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 2
- 238000006722 reduction reaction Methods 0.000 claims description 2
- 239000007789 gas Substances 0.000 description 13
- 239000011159 matrix material Substances 0.000 description 7
- 239000013618 particulate matter Substances 0.000 description 5
- 238000011069 regeneration method Methods 0.000 description 5
- 238000011217 control strategy Methods 0.000 description 4
- MWUXSHHQAYIFBG-UHFFFAOYSA-N nitrogen oxide Inorganic materials O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000033228 biological regulation Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000008929 regeneration Effects 0.000 description 3
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 2
- 230000004913 activation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 229910002091 carbon monoxide Inorganic materials 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 206010008428 Chemical poisoning Diseases 0.000 description 1
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 239000010687 lubricating oil Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 238000013021 overheating Methods 0.000 description 1
- 230000003647 oxidation Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 239000004071 soot Substances 0.000 description 1
- 229910052717 sulfur Inorganic materials 0.000 description 1
- 239000011593 sulfur Substances 0.000 description 1
Images
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
- F01N11/00—Monitoring or diagnostic devices for exhaust-gas treatment apparatus, e.g. for catalytic activity
-
- 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
-
- 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
- F01N9/00—Electrical control of exhaust gas treating apparatus
-
- 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
-
- 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
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Combustion & Propulsion (AREA)
- Mechanical Engineering (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Toxicology (AREA)
- Exhaust Gas After Treatment (AREA)
Abstract
本发明公开了一种柴油机后处理系统催化器热失效容错控制方法,包括以下步骤:获取发动机转速和转矩信息,SCR催化器上、下游温度值测量值,SCR催化器下游NOX浓度测量值;预测排气质量流量以及SCR催化器上游NOX浓度估计值;通过扩展卡尔曼滤波算法获取SCR催化器热失效程度信息、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值以及SCR催化器下游NH3浓度估计值;通过基于李雅普诺夫稳定性理论设计得到的模型参考自适应容错控制器计算尿素喷射量。本发明提出的柴油机后处理系统催化器热失效容错控制方法,能快速且准确诊断当前Urea‑SCR催化器热失效程度并进行容错控制。
Description
技术领域
本发明涉及发动机技术领域,尤其涉及一种柴油机后处理系统催化器热失 效容错控制方法及其装置。
背景技术
柴油机与汽油机相比,虽然其一氧化碳(CO)和碳氢化合物(HC)排放较 低,但其氮氧化物(NOX)和颗粒物(PM)的排放量却远超汽油机。Urea-SCR技 术因其良好的燃油经济性以及抗硫能力而成为减少柴油机NOX排放的主流后处 理技术。随着国IV/V排放法规的实施,Urea-SCR系统得以大规模应用。但由于 我国油品、润滑油质量、道路状况、法制监督等都存在一些问题,以及在汽车 行驶过程中,由于过热老化、化学中毒、积垢、机械损坏等原因使得SCR催化 器失效,从而使Urea-SCR系统的转化效率下降,同时还会造成二次污染(氨气 泄漏)。同时,随着国VI排放法规的颁布,NOX和PM排放限值进一步降低, DOC+DPF+SCR(Diesel oxidation catalyst柴油氧化催化器,Diesel particulate filter柴油微粒捕集器,Selective catalytic reduction选择性催化还原) 的集成式后处理系统将被应用于柴油机排放控制。其中,DPF在工作过程中捕集 炭烟颗粒而被堵塞,从而造成排气背压过高、柴油机动力性降低等问题。因此, 系统必须能自行进行PM的清洁,使DPF恢复到原来的工作状态,从而实现DPF 的再生。DPF再生方法通常分为两种:主动再生和被动再生。其中,主动再生方 法是将排气温度加热到500℃以上将PM通过高温燃烧,这将进一步加剧SCR催化器的热失效,而SCR催化器热失效会对Urea-SCR系统的排放控制性能造成不 利影响。
发明内容
本发明的主要目的在于提供一种能快速且准确诊断当前Urea-SCR催化器热 失效程度并进行容错控制的柴油机后处理系统催化器热失效容错控制方法及其 装置。
为实现上述目的,本发明提供一种柴油机后处理系统催化器热失效容错控 制方法,其特征在于,包括以下步骤:
获取发动机转速和转矩信息,从尿素喷射电控单元获取尿素喷射量信号, 根据温度传感器获取SCR催化器上、下游温度值测量值,根据NOX传感器获取SCR 催化器下游NOX浓度测量值;
原排NOX浓度估计模块根据当前发动机转速和转矩信息预测排气质量流量 以及SCR催化器上游NOX浓度估计值;
SCR催化器热失效诊断模块根据SCR催化器上游NOX浓度估计值、排气质量 流量、尿素喷射量、SCR催化器上下游温度测量值和SCR催化器下游NOX浓度测 量值通过扩展卡尔曼滤波算法获取SCR催化器热失效程度信息、SCR催化器氨覆 盖度估计值、SCR催化器下游NOX浓度估计值以及SCR催化器下游NH3浓度估计 值;
SCR催化器热失效容错模块基于SCR催化器上游NOX浓度估计值、排气质量 流量、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值、SCR催化器 下游NH3浓度估计值、SCR催化器上下游温度测量值、SCR催化器下游NOX浓度测 量值以及SCR催化器热失效程度信息,通过基于李雅普诺夫稳定性理论设计得 到的模型参考自适应容错控制器计算尿素喷射量。
优选地,通过扩展卡尔曼滤波算法获取SCR催化器热失效程度信息、SCR催 化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值以及SCR催化器下游NH3浓度估计值的步骤之后还包括:
通过SCR催化器热失效诊断信息有效性判定模块对从SCR催化器热失效诊 断模块获取的SCR催化器热失效程度信息的有效性进行判定后,输出并存储有 效的SCR催化器失效程度信息;
SCR催化器热失效容错模块基于SCR催化器上游NOX浓度估计值、排气质量 流量、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值、SCR催化器 下游NH3浓度估计值、SCR催化器上下游温度测量值、SCR催化器下游NOX浓度测 量值以及从SCR催化器热失效诊断信息有效性判定模块获取的有效的SCR催化 器热失效程度信息,通过基于李雅普诺夫稳定性理论设计得到的模型参考自适 应容错控制器计算尿素喷射量。
优选地,为衡量SCR催化器失效后储氨能力相比SCR催化器失效前储氨能 力减小的程度,SCR催化器热失效程度信息即为SCR催化器失效因子α,其定义 如下:
其中,Cs,aged为催化器老化状态下的最大储氨量,Cs,fresh为催化器健康状态下 的最大储氨量。
优选地,通过SCR催化器热失效诊断信息有效性判定模块对从SCR催化器 热失效诊断模块获取的SCR催化器热失效程度信息的有效性进行判定时,采用 以下方法:
计算最近多个SCR催化器失效因子估计值采样点的方差;
如果这多个个采样点的方差小于阈值,则说明SCR催化器失效诊断信息是 有效的,其可以用于容错控制;否则说明目前SCR催化器失效诊断信息不稳定, 则调用SCR催化器热失效诊断信息有效性判定模块存储的最近且稳定的失效诊 断信息用于容错控制。
优选地,获取SCR催化器热失效程度信息、SCR催化器氨覆盖度估计值、SCR 催化器下游NOX浓度估计值以及SCR催化器下游NH3浓度估计值采用以下方法:
SCR催化器状态向量表示为:
SCR催化器状态方程为:
SCR催化器状态预测模型为:
SCR催化器测量模型表示为:
上式中,nNOx,in为流入催化器的NOX的摩尔流量,为流入催化器的NH3的摩尔流量,为下游NOX的摩尔浓度,为NH3的摩尔浓度,Cs为催化剂 最大储氨量,Sc为1mol表面活性原子面积,αprob为粘着概率,ki为化学反应频 率因子,θ为氨覆盖度,R为气体常数,T为催化剂床温,Ei为反应活化能,其 中,i为ads时表示氨吸附反应,i为des时表示氨脱附反应,i为scr时表示选择 性催化还原反应,i为ox表示氨氧化反应,为氨覆盖度估计值,为下游NOX的摩尔浓度估计值,为NH3的摩尔浓度估计值,mEG为排气质量流量,REG为 排气气体常数,pamb为环境大气压,ε为催化器孔隙率,Vc为催化器体积,cp,c为催化器比热容,cp,EG排气定压比热容,mc为催化器质量,εrad为黑度,σ为气 体辐射常数,MNH3表示NH3的摩尔质量;
通过NOX传感器获取SCR催化器下游NOX浓度测量值,代入SCR催化器测量 模型中,得到SCR催化器下游NOX浓度实际值以及SCR催化器下游NH3浓度实际 值,并根据SCR催化器下游NOX浓度实际值以及SCR催化器下游NH3浓度实际值 计算得到氨覆盖度实际值后,根据SCR催化器状态预测模型和SCR催化器状态 方程,计算得到SCR催化器失效因子、SCR催化器氨覆盖度估计值、SCR催化器 下游NOX浓度估计值以及SCR催化器下游NH3浓度估计值。
优选地,SCR催化器热失效容错模块基于SCR催化器上游NOX浓度估计值、 排气质量流量、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值、 SCR催化器下游NH3浓度估计值、SCR催化器上下游温度测量值、SCR催化器下游 NOX浓度测量值以及SCR催化器热失效程度信息,通过基于李雅普诺夫稳定性理 论设计得到的模型参考自适应容错控制器计算尿素喷射量具体采用以下方法:
模型参考自适应容错控制器中参考模型为:
其中,θref *为氨覆盖度参考值,θref'为初始氨覆盖度参考值,θref为MRAC氨 覆盖度参考模型值;
Urea-SCR系统尿素喷射模型表示如下:
定义:
Urea-SCR系统尿素喷射模型重新表示为:
其中:
优选地,自适应容错控制器的自适应调节机构设计为:
其中:
其中,λ1和λ2均为正常数;
通过SCR催化器热失效诊断信息有效性判定模块获取的SCR催化器失效因 子用于修正模型参考自适应容错控制器中的参考模型以及自适应调节机构,输 入信号同时进入参考模型以及修正模型参考自适应容错控制器中进行计算,其 中参考模型直接得到氨覆盖度参考值,控制器将控制信号传递给受控系统执行 器执行相应指令从而得到氨覆盖度实际值;比较氨覆盖度参考值与氨覆盖度实 际值,两者差值将会传递给自适应调节机构对控制器进行参数进行修正;当氨 覆盖度参考值与氨覆盖度实际值相同时,自适应调节终止。
本发明进一步提出一种基于上述的柴油机后处理系统催化器热失效容错控 制方法的控制装置,包括:
处理器,用于实现各指令;
存储器,用于存储多条指令,所述指令由处理器加载并执行:
获取发动机转速和转矩信息,从尿素喷射电控单元获取尿素喷射量信号, 根据温度传感器获取SCR催化器上、下游温度值测量值,根据NOX传感器获取SCR 催化器下游NOX浓度测量值;
原排NOX浓度估计模块根据当前发动机转速和转矩信息预测排气质量流量 以及SCR催化器上游NOX浓度估计值;
SCR催化器热失效诊断模块根据SCR催化器上游NOX浓度估计值、排气质量 流量、尿素喷射量、SCR催化器上下游温度测量值和SCR催化器下游NOX浓度测 量值通过扩展卡尔曼滤波算法获取SCR催化器热失效程度信息、SCR催化器氨覆 盖度估计值、SCR催化器下游NOX浓度估计值以及SCR催化器下游NH3浓度估计 值;
SCR催化器热失效容错模块基于SCR催化器上游NOX浓度估计值、排气质量 流量、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值、SCR催化器 下游NH3浓度估计值、SCR催化器上下游温度测量值、SCR催化器下游NOX浓度测 量值以及SCR催化器热失效程度信息,通过基于李雅普诺夫稳定性理论设计得 到的模型参考自适应容错控制器计算尿素喷射量。
本发明提出的柴油机后处理系统催化器热失效容错控制方法,具有以下有 益效果:
(1)本发明基于扩展卡尔曼滤波算法构建的SCR催化器热失效诊断模块可 有效抵御传感器测量噪声对SCR催化器失效程度估计的影响,可快速且精确地 诊断SCR催化器热失效程度;
(2)本发明基于李雅普诺夫稳定性理论设计得到的模型参考自适应容错控 制器具有结构简单,易于实现和计算负荷小等特点,可实时地根据SCR催化器 热失效程度自适应调节尿素喷射量,从而使得Urea-SCR系统在SCR催化器老化 失效的情况下仍然能保持良好的排放控制性能;
(3)本发明可广泛应用于SCR催化器失效诊断和容错控制中。
附图说明
图1为本发明柴油机后处理系统催化器热失效容错控制方法中Urea-SCR系 统示意图;
图2为本发明柴油机后处理系统催化器热失效容错控制方法的流程示意图。
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说 明。
具体实施方式
应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定 本发明。
参照图1和图2,本实施例中,一种柴油机后处理系统催化器热失效容错控 制方法,包括以下步骤:
步骤S10,获取发动机转速和转矩信息,从尿素喷射电控单元获取尿素喷射 量信号,根据温度传感器获取SCR催化器上、下游温度值测量值,根据NOX传感 器获取SCR催化器下游NOX浓度测量值;
步骤S20,原排NOX浓度估计模块根据当前发动机转速和转矩信息预测排气 质量流量以及SCR催化器上游NOX浓度估计值;
步骤S30,SCR催化器热失效诊断模块根据SCR催化器上游NOX浓度估计值、 排气质量流量、尿素喷射量、SCR催化器上下游温度测量值和SCR催化器下游 NOX浓度测量值通过扩展卡尔曼滤波算法获取SCR催化器热失效程度信息、SCR 催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值以及SCR催化器下游NH3浓度估计值;
步骤S40,SCR催化器热失效容错模块基于SCR催化器上游NOX浓度估计值、 排气质量流量、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值、 SCR催化器下游NH3浓度估计值、SCR催化器上下游温度测量值、SCR催化器下游 NOX浓度测量值以及SCR催化器热失效程度信息,通过基于李雅普诺夫稳定性理 论设计得到的模型参考自适应容错控制器计算尿素喷射量。
进一步地,步骤S30之后还包括:
步骤S50,通过SCR催化器热失效诊断信息有效性判定模块对从SCR催化器 热失效诊断模块获取的SCR催化器热失效程度信息的有效性进行判定后,输出 并存储有效的SCR催化器失效程度信息;
步骤S60,SCR催化器热失效容错模块基于SCR催化器上游NOX浓度估计值、 排气质量流量、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值、 SCR催化器下游NH3浓度估计值、SCR催化器上下游温度测量值、SCR催化器下游 NOX浓度测量值以及从SCR催化器热失效诊断信息有效性判定模块获取的有效的 SCR催化器热失效程度信息,通过基于李雅普诺夫稳定性理论设计得到的模型参 考自适应容错控制器计算尿素喷射量。
即将步骤S40替换为步骤S50和S60,此时增加了一步判断的过程,即对SCR 催化器热失效程度信息的有效性进行判断,从而使本热失效容错控制方法更加 精确。
具体地,步骤S50中,计算最近多个SCR催化器失效因子估计值采样点的 方差;
如果这多个个采样点的方差小于阈值,则说明SCR催化器失效诊断信息是 有效的,其可以用于容错控制;否则说明目前SCR催化器失效诊断信息不稳定, 则调用SCR催化器热失效诊断信息有效性判定模块存储的最近且稳定的失效诊 断信息用于容错控制。
以下具体说明本柴油机后处理系统催化器热失效容错控制方法的具体过程。
根据SCR催化器内主要化学反应和质量守恒定律,SCR催化器状态方程可描 述如下:
其中,nNOx,in为流入催化器的NOx的摩尔流量(单位mol/s),为流入催 化器的NH3的摩尔流量(单位mol/s),为下游NOX的摩尔浓度(单位mol/m3), 为NH3的摩尔浓度(单位mol/m3),Cs为催化剂最大储氨量(单位单位mol/m3), Sc为1mol表面活性原子面积(单位m2/mol),αprob为粘着概率,ki为化学反应频 率因子,Ei为反应活化能(单位J/mol),θ为氨覆盖度,R为气体常数(单位 8.3145J/molK),T为催化剂床温(单位K)。其中,i代表ads(氨吸附反应)、des(氨 脱附反应)、scr(选择性催化还原反应)和ox(氨氧化反应),为氨覆盖度估计 值,为下游NOX的摩尔浓度估计值,为NH3的摩尔浓度估计值,mEG为排气质量流量。
式中:
其中,REG为排气气体常数(单位J/kgK),pamb为环境大气压(单位Pa), ε为催化器孔隙率,Vc为催化器体积(单位m3),cp,c为催化器比热容(单位J/kgK), cp,EG排气定压比热容(单位J/kgK),mc为催化器质量(单位kg),εrad为黑度, σ为气体辐射常数(单位W/m2K4),MNH3表示NH3的摩尔质量。
根据SCR催化器失效机理,SCR催化器失效将会导致其储氨能力下降。为衡 量SCR催化器失效后储氨能力相比SCR催化器失效前储氨能力减小的程度,SCR 催化器失效因子α定义如下:
其中,Cs,aged为催化器老化状态下的最大储氨量,Cs,fresh为催化器健康状态下 的最大储氨量。
SCR催化器失效状态方程:
非线性系统经过EKF算法处理可表现为以下形式:
其中,f(x,u)为非线性状态函数,h(x)为非线性测量函数。
EKF同样通过预测和更新两个步骤来估计系统的状态参数。在确定了状态参 数初值x(0|0)以及误差协方差矩阵初值P(0|0)后,便可进入EKF的预测步骤。
在预测步骤中,状态量和误差协方差矩阵在k-1时刻对k时刻的预测如下:
x(k|k-1)=f[x(k-1|k-1),u(k)] (6)
P(k|k-1)=F(k)P(k-1|k-1)F(k)T+Q(k) (7)
其中,F为f(x,u)的Jacobian矩阵。
在更新步骤中,EKF会通过卡尔曼增益K(k)对系统状态预测值x(k|k-1)和 预测的误差协方差矩阵P(k|k-1)进行更新,从而得到系统状态参数估计值x(k|k) 以及估计的误差协方差矩阵P(k|k)。假设噪声为零均值高斯白噪声,卡尔曼增 益K(k)、系统状态参数估计值x(k|k)以及估计的误差协方差矩阵P(k|k)可由以 下方程得到:
K(k)=P(k|k-1)H(k)T[H(k)P(k|k-1)H(k)T+R(k)]-1 (8)
x(k|k)=x(k|k-1)+K(k){z(k)-h[x(k|k-1),u(k)]} (9)
P(k|k)=[I-K(k)H(k)]P(k|k-1) (10)
其中,H为h(x)的Jacobian矩阵。
SCR催化器状态向量应表示如下:
其中,
SCR催化器状态方程如下所示:
那么根据式(13),SCR催化器状态预测模型可表示如下:
根据NOx传感器数学模型,SCR催化器测量模型可表示如下:
其中,Kcs为NOX和NH3的交叉感应系数。
在确定SCR催化器状态预测模型和测量模型之后,便可根据式(6)~(10) 进行迭代运算得到SCR催化器状态参数(包括失效因子、SCR催化器氨覆盖度估 计值、SCR催化器下游NOX浓度估计值以及SCR催化器下游NH3浓度估计值这四 个参数)的最优估计。
SCR催化器状态方程用于描述正常催化器工作状态,SCR催化器失效状态方 程用于描述失效时催化器工作状态,SCR催化器状态预测模型用于在设计EKF观 测器时需要确定SCR催化器状态方程与失效因子的状态方程,通过这个模型来 预测四个状态参数。SCR催化器测量模型:在SCR催化器下游安装NOX传感器来 测量NOX浓度,通过该模型得到NOX和NH3的浓度测量值即作为对应实际值后代入 公式(13)和(14)(因NH3传感器成本太高,利用NH3与NOX的交感反应,借由 NOX浓度推得NH3浓度)。
通过NOX传感器获取SCR催化器下游NOX浓度测量值,代入SCR催化器测量 模型中,得到SCR催化器下游NOX浓度实际值以及SCR催化器下游NH3浓度实际 值,并根据SCR催化器下游NOX浓度实际值以及SCR催化器下游NH3浓度实际值 计算得到氨覆盖度实际值后,根据SCR催化器状态预测模型和SCR催化器状态 方程,计算得到SCR催化器失效因子、SCR催化器氨覆盖度估计值、SCR催化器 下游NOX浓度估计值以及SCR催化器下游NH3浓度估计值。
模型参考自适应容错控制器设计方法如下:
氨覆盖度参考值是SCR催化器床温的单值函数,并可表示如下:
为了使SCR催化器失效前后的储氨量保持在同一水平,模型参考自适应容 错控制器参考模型可设计如下:
其中,θref *为氨覆盖度参考值,θref'为初始氨覆盖度参考值,θref为MRAC氨 覆盖度参考模型值;
Urea-SCR系统尿素喷射模型可表示如下:
定义:
故式(21)可重新表示如下:
其中:
MRAC的自适应调节机构可设计如下:
其中:
其中,λ1和λ2均为正常数。
模型参考自适应控制(MRAC)原理:通过SCR催化器热失效诊断信息有效 性判定模块获取的SCR催化器失效因子用于修正模型参考自适应容错控制器中 的参考模型以及自适应调节机构,输入信号同时进入参考模型以及修正模型参 考自适应容错控制器中进行计算,其中参考模型直接得到氨覆盖度参考值,控 制器将控制信号传递给受控系统执行器执行相应指令从而得到氨覆盖度实际值; 比较氨覆盖度参考值与氨覆盖度实际值,两者差值将会传递给自适应调节机构 对控制器进行参数进行修正;当氨覆盖度参考值与氨覆盖度实际值相同时,自 适应调节终止。
公式(19-20)为MRAC参考模型设计,(25-26)为MRAC自适应调节机构设 计,(19-20)由于催化剂的最大储氨量会随着SCR催化器失效程度的增加而减 少,即失效因子α会下降至0到1之间的某一个值,为了使SCR催化器失效前 后的储氨量保持在同一水平,那么氨覆盖度参考值也应该自适应地随着SCR催 化器失效程度的增加而上升,(21-24)为了能使SCR催化器实际氨覆盖度逼近 氨覆盖度参考值,Urea-SCR系统是仅通过控制尿素喷射器喷射适量尿素来实现。
本柴油机后处理系统催化器热失效容错控制方法,为了模拟SCR催化器失 效,引入失效因子来衡量催化器失效前后储氨能力的程度。在SCR催化器状态 方程加入失效因子,得到SCR催化器失效状态方程。然后,基于EKF算法,设 计EKF状态观测器,对SCR催化器状态进行估计,得到失效因子的估计值,即 失效信息。再根据得到的失效信息设计MRAC自适应控制器,从而实现容错控制。
经过测试表明,在采用原尿素喷射控制策略时,随着SCR催化器失效程度 增加,原尿素喷射控制策略下的NOX转化效率会随之下降,原尿素喷射控制策略 下的氨泄漏量也会随之增加。采用本发明的容错控制后,在进行ETC测试循环 时,容错控制下NOX转化效率对催化剂失效并不敏感,容错控制下氨泄漏量相比 原尿素喷射控制策略下的氨泄漏量仅轻微上升。
本发明提出的柴油机后处理系统催化器热失效容错控制方法,具有以下有 益效果:
(1)本发明基于扩展卡尔曼滤波算法构建的SCR催化器热失效诊断模块可 有效抵御传感器测量噪声对SCR催化器失效程度估计的影响,可快速且精确地 诊断SCR催化器热失效程度;
(2)本发明基于李雅普诺夫稳定性理论设计得到的模型参考自适应容错控 制器具有结构简单,易于实现和计算负荷小等特点,可实时地根据SCR催化器 热失效程度自适应调节尿素喷射量,从而使得Urea-SCR系统在SCR催化器老化 失效的情况下仍然能保持良好的排放控制性能;
(3)本发明可广泛应用于SCR催化器失效诊断和容错控制中。
本发明进一步提出一种柴油机后处理系统催化器热失效容错控制方法的控 制装置。
本优选实施例中,柴油机后处理系统催化器热失效容错控制方法的控制装 置,包括:
处理器,用于实现各指令;
存储器,用于存储多条指令,所述指令由处理器加载并执行:
获取发动机转速和转矩信息,从尿素喷射电控单元获取尿素喷射量信号, 根据温度传感器获取SCR催化器上、下游温度值测量值,根据NOX传感器获取SCR 催化器下游NOX浓度测量值;
原排NOX浓度估计模块根据当前发动机转速和转矩信息预测排气质量流量 以及SCR催化器上游NOX浓度估计值;
SCR催化器热失效诊断模块根据SCR催化器上游NOX浓度估计值、排气质量 流量、尿素喷射量、SCR催化器上下游温度测量值和SCR催化器下游NOX浓度测 量值通过扩展卡尔曼滤波算法获取SCR催化器热失效程度信息、SCR催化器氨覆 盖度估计值、SCR催化器下游NOX浓度估计值以及SCR催化器下游NH3浓度估计 值;
SCR催化器热失效容错模块基于SCR催化器上游NOX浓度估计值、排气质量 流量、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值、SCR催化器 下游NH3浓度估计值、SCR催化器上下游温度测量值、SCR催化器下游NOX浓度 测量值以及SCR催化器热失效程度信息,通过基于李雅普诺夫稳定性理论设计 得到的模型参考自适应容错控制器计算尿素喷射量。
控制装置中包括发动机电控单元,尿素喷射电控单元,尿素喷嘴,SCR催化 器上下游温度传感器,SCR催化器下游NOX传感器,尿素喷射电控单元包括柴油 机原排NOX浓度估计模块,柴油机排气质量流量估计模块,SCR催化器失效诊断 模块,SCR催化器失效诊断信息有效性判定模块和SCR催化器失效容错模块。柴 油机原排NOX浓度估计模块,柴油机排气质量流量估计模块,SCR催化器失效诊 断模块,SCR催化器失效诊断信息有效性判定模块和SCR催化器失效容错模块均 集成于处理器中。
以上仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利 用本发明说明书及附图内容所作的等效结构变换,或直接或间接运用在其他相 关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (2)
1.一种柴油机后处理系统催化器热失效容错控制方法,其特征在于,包括以下步骤:
获取发动机转速和转矩信息,从尿素喷射电控单元获取尿素喷射量信号,根据温度传感器获取SCR催化器上、下游温度值测量值,根据NOX传感器获取SCR催化器下游NOX浓度测量值;
原排NOX浓度估计模块根据当前发动机转速和转矩信息预测排气质量流量以及SCR催化器上游NOX浓度估计值;
SCR催化器热失效诊断模块根据SCR催化器上游NOX浓度估计值、排气质量流量、尿素喷射量、SCR催化器上下游温度测量值和SCR催化器下游NOX浓度测量值通过扩展卡尔曼滤波算法获取SCR催化器热失效程度信息、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值以及SCR催化器下游NH3浓度估计值;
SCR催化器热失效容错模块基于SCR催化器上游NOX浓度估计值、排气质量流量、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值、SCR催化器下游NH3浓度估计值、SCR催化器上下游温度测量值、SCR催化器下游NOX浓度测量值以及SCR催化器热失效程度信息,通过基于李雅普诺夫稳定性理论设计得到的模型参考自适应容错控制器计算尿素喷射量;
通过扩展卡尔曼滤波算法获取SCR催化器热失效程度信息、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值以及SCR催化器下游NH3浓度估计值的步骤之后还包括:
通过SCR催化器热失效诊断信息有效性判定模块对从SCR催化器热失效诊断模块获取的SCR催化器热失效程度信息的有效性进行判定后,输出并存储有效的SCR催化器失效程度信息;
SCR催化器热失效容错模块基于SCR催化器上游NOX浓度估计值、排气质量流量、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值、SCR催化器下游NH3浓度估计值、SCR催化器上下游温度测量值、SCR催化器下游NOX浓度测量值以及从SCR催化器热失效诊断信息有效性判定模块获取的有效的SCR催化器热失效程度信息,通过基于李雅普诺夫稳定性理论设计得到的模型参考自适应容错控制器计算尿素喷射量;
为衡量SCR催化器失效后储氨能力相比SCR催化器失效前储氨能力减小的程度,SCR催化器热失效程度信息即为SCR催化器失效因子α,其定义如下:
其中,Cs,aged为催化器老化状态下的最大储氨量,Cs,fresh为催化器健康状态下的最大储氨量;
通过SCR催化器热失效诊断信息有效性判定模块对从SCR催化器热失效诊断模块获取的SCR催化器热失效程度信息的有效性进行判定时,采用以下方法:
计算最近多个SCR催化器失效因子估计值采样点的方差;
如果这多个采样点的方差小于阈值,则说明SCR催化器失效诊断信息是有效的,其可以用于容错控制;否则说明目前SCR催化器失效诊断信息不稳定,则调用SCR催化器热失效诊断信息有效性判定模块存储的最近且稳定的失效诊断信息用于容错控制;
获取SCR催化器热失效程度信息、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值以及SCR催化器下游NH3浓度估计值采用以下方法:
SCR催化器状态向量表示为:
SCR催化器状态方程为:
SCR催化器状态预测模型为:
SCR催化器测量模型表示为:
上式中,nNOx,in为流入催化器的NOX的摩尔流量,为流入催化器的NH3的摩尔流量,为下游NOX的摩尔浓度,为NH3的摩尔浓度,Cs为催化剂最大储氨量,Sc为1mol表面活性原子面积,αprob为粘着概率,ki为化学反应频率因子,θ为氨覆盖度,R为气体常数,T为催化剂床温,Ei为反应活化能,其中,i为ads时表示氨吸附反应,i为des时表示氨脱附反应,i为scr时表示选择性催化还原反应,i为ox表示氨氧化反应,为氨覆盖度估计值,为下游NOX的摩尔浓度估计值,为NH3的摩尔浓度估计值,mEG为排气质量流量,REG为排气气体常数,pamb为环境大气压,ε为催化器孔隙率,Vc为催化器体积,cp,c为催化器比热容,cp,EG排气定压比热容,mc为催化器质量,εrad为黑度,σ为气体辐射常数,MNH3表示NH3的摩尔质量;
通过NOX传感器获取SCR催化器下游NOX浓度测量值,代入SCR催化器测量模型中,得到SCR催化器下游NOX浓度实际值以及SCR催化器下游NH3浓度实际值,并根据SCR催化器下游NOX浓度实际值以及SCR催化器下游NH3浓度实际值计算得到氨覆盖度实际值后,根据SCR催化器状态预测模型和SCR催化器状态方程,计算得到SCR催化器失效因子、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值以及SCR催化器下游NH3浓度估计值;
SCR催化器热失效容错模块基于SCR催化器上游NOX浓度估计值、排气质量流量、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值、SCR催化器下游NH3浓度估计值、SCR催化器上下游温度测量值、SCR催化器下游NOX浓度测量值以及SCR催化器热失效程度信息,通过基于李雅普诺夫稳定性理论设计得到的模型参考自适应容错控制器计算尿素喷射量具体采用以下方法:
模型参考自适应容错控制器中参考模型为:
其中,θref *为氨覆盖度参考值,θref'为初始氨覆盖度参考值,θref为MRAC氨覆盖度参考模型值;
Urea-SCR系统尿素喷射模型表示如下:
定义:
Urea-SCR系统尿素喷射模型重新表示为:
其中:
自适应容错控制器的自适应调节机构设计为:
其中:
其中,λ1和λ2均为正常数;
通过SCR催化器热失效诊断信息有效性判定模块获取的SCR催化器失效因子用于修正模型参考自适应容错控制器中的参考模型以及自适应调节机构,输入信号同时进入参考模型以及修正模型参考自适应容错控制器中进行计算,其中参考模型直接得到氨覆盖度参考值,控制器将控制信号传递给受控系统执行器执行相应指令从而得到氨覆盖度实际值;比较氨覆盖度参考值与氨覆盖度实际值,两者差值将会传递给自适应调节机构对控制器进行参数进行修正;当氨覆盖度参考值与氨覆盖度实际值相同时,自适应调节终止。
2.一种基于权利要求1所述的柴油机后处理系统催化器热失效容错控制方法的控制装置,其特征在于,包括:
处理器,用于实现各指令;
存储器,用于存储多条指令,所述指令由处理器加载并执行:
获取发动机转速和转矩信息,从尿素喷射电控单元获取尿素喷射量信号,根据温度传感器获取SCR催化器上、下游温度值测量值,根据NOX传感器获取SCR催化器下游NOX浓度测量值;
原排NOX浓度估计模块根据当前发动机转速和转矩信息预测排气质量流量以及SCR催化器上游NOX浓度估计值;
SCR催化器热失效诊断模块根据SCR催化器上游NOX浓度估计值、排气质量流量、尿素喷射量、SCR催化器上下游温度测量值和SCR催化器下游NOX浓度测量值通过扩展卡尔曼滤波算法获取SCR催化器热失效程度信息、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值以及SCR催化器下游NH3浓度估计值;
SCR催化器热失效容错模块基于SCR催化器上游NOX浓度估计值、排气质量流量、SCR催化器氨覆盖度估计值、SCR催化器下游NOX浓度估计值、SCR催化器下游NH3浓度估计值、SCR催化器上下游温度测量值、SCR催化器下游NOX浓度测量值以及SCR催化器热失效程度信息,通过基于李雅普诺夫稳定性理论设计得到的模型参考自适应容错控制器计算尿素喷射量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910285859.4A CN110273737B (zh) | 2019-04-10 | 2019-04-10 | 柴油机后处理系统催化器热失效容错控制方法及其装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910285859.4A CN110273737B (zh) | 2019-04-10 | 2019-04-10 | 柴油机后处理系统催化器热失效容错控制方法及其装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110273737A CN110273737A (zh) | 2019-09-24 |
CN110273737B true CN110273737B (zh) | 2021-08-17 |
Family
ID=67959474
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910285859.4A Active CN110273737B (zh) | 2019-04-10 | 2019-04-10 | 柴油机后处理系统催化器热失效容错控制方法及其装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110273737B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110863890A (zh) * | 2019-10-09 | 2020-03-06 | 华东交通大学 | 一种远程诊断柴油车scr系统尿素消耗量合理性的方法 |
CN110645077B (zh) * | 2019-10-15 | 2021-11-02 | 黄冈师范学院 | 柴油机Urea-SCR系统氨喷射速率和状态同步估计方法 |
CN111044907B (zh) * | 2019-12-24 | 2022-06-14 | 苏州正力新能源科技有限公司 | 一种基于微片段数据和电压滤波的soh统计方法 |
CN116357435A (zh) * | 2023-03-30 | 2023-06-30 | 中国汽车技术研究中心有限公司 | 柴油发动机排放的控制方法、设备及介质 |
CN116642562B (zh) * | 2023-07-27 | 2023-10-20 | 黑龙江惠达科技股份有限公司 | 一种植保无人机药液质量测量系统、方法和无人机 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104234802B (zh) * | 2014-07-14 | 2017-01-11 | 浙江大学 | 基于NOx反馈和储氨预测的SCR催化器老化判定方法 |
CN107023367B (zh) * | 2017-03-29 | 2019-04-12 | 北京航空航天大学 | 一种柴油机scr系统氨气输入传感器故障诊断和容错控制方法 |
CN107620620B (zh) * | 2017-09-05 | 2022-12-20 | 天津大学 | 一种实验教学用选择性催化还原发动机尾气系统 |
CN108087065B (zh) * | 2017-11-01 | 2019-08-27 | 联创汽车电子有限公司 | Scr催化器氨浓度修正方法及scr处理系统 |
CN108425730A (zh) * | 2018-02-27 | 2018-08-21 | 西南交通大学 | 柴油机scr后处理系统的obd诊断及处理方法 |
CN108915827B (zh) * | 2018-07-04 | 2020-02-11 | 中国汽车技术研究中心有限公司 | 基于SCR化学反应数学模型改善发动机NOx排放量的方法 |
-
2019
- 2019-04-10 CN CN201910285859.4A patent/CN110273737B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110273737A (zh) | 2019-09-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110273737B (zh) | 柴油机后处理系统催化器热失效容错控制方法及其装置 | |
US11313268B2 (en) | Method of monitoring an SCR catalyst | |
EP2813693B1 (en) | Systems and methods for NOx sensor diagnostics | |
US8596045B2 (en) | On-board-diagnosis method for an exhaust aftertreatment system and on-board-diagnosis system for an exhaust aftertreatment system | |
US8893475B2 (en) | Control system for doser compensation in an SCR system | |
US8051645B2 (en) | Determination of diesel particulate filter load under both transient and steady state drive cycles | |
EP2180157B1 (en) | Diagnostic method for selective catalytic reduction (SCR) exhaust treatment system | |
US8474248B2 (en) | Model based method for selective catalyst reducer urea dosing strategy | |
US8726723B2 (en) | Detection of aftertreatment catalyst degradation | |
US20120137657A1 (en) | Scr closed loop control system | |
CN109931129B (zh) | 用于监测内燃机的排气后处理系统的方法和装置 | |
US10100701B2 (en) | Method for the diagnosis of an exhaust gas aftertreatment system for an internal combustion engine | |
EP2940280B1 (en) | Fuel-cetane-number estimation method and apparatus | |
CN110295978A (zh) | 柴油机scr控制系统温度传感器容错控制方法及其装置 | |
US9261002B2 (en) | Feed-back for closed-loop selective catalytic reduction control | |
JP2020118077A (ja) | 内燃機関の排気浄化装置、及び車両 | |
WO2010096793A2 (en) | Detection of aftertreatment catalyst degradation | |
CN113356983B (zh) | 监测排气后处理装置的方法 | |
US20230141389A1 (en) | Diagnosis device and diagnosis method | |
CN113514169A (zh) | 一种scr系统下游温度传感器可信性故障诊断方法 | |
US20190249586A1 (en) | Abnormality diagnosis apparatus and vehicle | |
US9441519B2 (en) | System variation adaption for feed-forward controller | |
KR102563490B1 (ko) | 배기 정화 장치의 내구수명 향상을 위한 배기 정화 장치의 재생 제어 방법 | |
US20240003284A1 (en) | UREA CRYSTAL DETECTION BY MONITORING A RATIO BETWEEN A NOx SENSOR SIGNAL AND A EQUIVALENCE RATIO | |
Haga et al. | New urea-SCR control system for super clean diesel |
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 |