JP2018127921A - Engine control device - Google Patents
Engine control device Download PDFInfo
- Publication number
- JP2018127921A JP2018127921A JP2017020560A JP2017020560A JP2018127921A JP 2018127921 A JP2018127921 A JP 2018127921A JP 2017020560 A JP2017020560 A JP 2017020560A JP 2017020560 A JP2017020560 A JP 2017020560A JP 2018127921 A JP2018127921 A JP 2018127921A
- Authority
- JP
- Japan
- Prior art keywords
- sox
- amount
- temperature
- adsorption
- distribution
- 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
- 238000003795 desorption Methods 0.000 claims abstract description 50
- 239000000779 smoke Substances 0.000 claims abstract description 14
- 238000001179 sorption measurement Methods 0.000 claims description 166
- 238000009826 distribution Methods 0.000 claims description 155
- 238000000746 purification Methods 0.000 claims description 51
- 229920006395 saturated elastomer Polymers 0.000 claims description 48
- 238000004364 calculation method Methods 0.000 claims description 36
- 238000006243 chemical reaction Methods 0.000 claims description 15
- QAOWNCQODCNURD-UHFFFAOYSA-L Sulfate Chemical compound [O-]S([O-])(=O)=O QAOWNCQODCNURD-UHFFFAOYSA-L 0.000 claims description 8
- 230000007423 decrease Effects 0.000 claims description 6
- 239000010419 fine particle Substances 0.000 claims description 2
- 239000003054 catalyst Substances 0.000 abstract description 26
- 230000008929 regeneration Effects 0.000 abstract description 14
- 238000011069 regeneration method Methods 0.000 abstract description 14
- 230000001629 suppression Effects 0.000 abstract description 4
- 229910052815 sulfur oxide Inorganic materials 0.000 description 228
- 239000000446 fuel Substances 0.000 description 43
- 238000010586 diagram Methods 0.000 description 21
- 238000000034 method Methods 0.000 description 21
- 239000000463 material Substances 0.000 description 20
- 239000011248 coating agent Substances 0.000 description 18
- 238000000576 coating method Methods 0.000 description 18
- 238000002347 injection Methods 0.000 description 13
- 239000007924 injection Substances 0.000 description 13
- 230000006870 function Effects 0.000 description 12
- 229910000510 noble metal Inorganic materials 0.000 description 9
- 238000007254 oxidation reaction Methods 0.000 description 8
- 229930195733 hydrocarbon Natural products 0.000 description 7
- 150000002430 hydrocarbons Chemical class 0.000 description 7
- 230000000630 rising effect Effects 0.000 description 7
- 230000003647 oxidation Effects 0.000 description 6
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 4
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 4
- 229910002091 carbon monoxide Inorganic materials 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 239000010410 layer Substances 0.000 description 4
- 229910052717 sulfur Inorganic materials 0.000 description 4
- 239000011593 sulfur Substances 0.000 description 4
- 239000002344 surface layer Substances 0.000 description 4
- 230000005484 gravity Effects 0.000 description 3
- 239000003921 oil Substances 0.000 description 3
- YREOLPGEVLLKMB-UHFFFAOYSA-N 3-methylpyridin-1-ium-2-amine bromide hydrate Chemical compound O.[Br-].Cc1ccc[nH+]c1N YREOLPGEVLLKMB-UHFFFAOYSA-N 0.000 description 2
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 230000001590 oxidative effect Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 238000000137 annealing Methods 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 239000010687 lubricating oil Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 229910052763 palladium Inorganic materials 0.000 description 1
- 229910052697 platinum Inorganic materials 0.000 description 1
- 239000010970 precious metal Substances 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- XTQHKBHJIVJGKJ-UHFFFAOYSA-N sulfur monoxide Chemical class S=O XTQHKBHJIVJGKJ-UHFFFAOYSA-N 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Exhaust Gas After Treatment (AREA)
- Processes For Solid Components From Exhaust (AREA)
Abstract
Description
本発明は、エンジン制御装置に関し、詳しくは、排気管に触媒を備えるディーゼルエンジンの制御装置に関する。 The present invention relates to an engine control device, and more particularly to a control device for a diesel engine provided with a catalyst in an exhaust pipe.
従来、ディーゼルエンジンの排気管に設けられる触媒に吸着した硫黄酸化物(SO2またはSO3をいい、以下においてこれらを区別しない場合には「SOx」と総称する。)を、定期的に脱離させる昇温制御を行うことが知られている。昇温制御に関連する文献として、例えば特開2013−029038号公報が挙げられる。この公報には、触媒に蓄積したSOxの蓄積量を推定し、推定した蓄積量が要求放出量に到達したときに、当該触媒の床温を上昇させて500〜550℃に制御する技術が開示されている。この公報によれば、床温が500℃未満であると触媒からSOxが放出されず、床温が500〜550℃であると触媒からSOxが低濃度で放出され、床温が600℃を超えると触媒からSOxが高濃度で放出されるという特性が触媒にあるとしている。従って、このような特性を有する触媒の床温を500〜550℃に制御すれば、触媒からSOxを低濃度で脱離させてその機能を回復させることができる。また、低濃度で脱離させたSOxから白煙が発生するのを抑制することもできる。 Conventionally, sulfur oxides adsorbed on a catalyst provided in an exhaust pipe of a diesel engine (referred to as SO 2 or SO 3 , hereinafter collectively referred to as “SOx” unless they are distinguished below) are periodically desorbed. It is known to perform temperature rise control. For example, Japanese Patent Application Laid-Open No. 2013-029038 is cited as a document related to temperature rise control. This publication discloses a technique for estimating the accumulated amount of SOx accumulated in the catalyst and, when the estimated accumulated amount reaches the required release amount, increasing the bed temperature of the catalyst and controlling it to 500 to 550 ° C. Has been. According to this publication, when the bed temperature is less than 500 ° C., SOx is not released from the catalyst, and when the bed temperature is 500 to 550 ° C., SOx is released from the catalyst at a low concentration, and the bed temperature exceeds 600 ° C. The catalyst has the characteristic that SOx is released at a high concentration from the catalyst. Therefore, if the bed temperature of the catalyst having such characteristics is controlled to 500 to 550 ° C., the function can be recovered by desorbing SOx from the catalyst at a low concentration. It is also possible to suppress the generation of white smoke from SOx desorbed at a low concentration.
上記公報の技術は、触媒の床温に対する触媒からのSOx放出量の特性に着目している。しかしながら、触媒から脱離するSOx量は、触媒の床温だけでなく床温の変化量の影響も受ける。このため、上記公報の技術では、触媒の温度変化が大きい条件において、触媒から脱離するSOx量を精度よく推定できないおそれがある。 The technique of the above publication pays attention to the characteristics of the SOx release amount from the catalyst with respect to the bed temperature of the catalyst. However, the amount of SOx desorbed from the catalyst is affected not only by the bed temperature of the catalyst but also by the amount of change in the bed temperature. For this reason, in the technique of the above publication, there is a possibility that the amount of SOx desorbed from the catalyst cannot be accurately estimated under conditions where the temperature change of the catalyst is large.
本発明は、上述のような課題に鑑みてなされたもので、触媒の昇温制御において、触媒から脱離するSOx量を精度よく推定することにより、SOxの脱離による触媒の機能の回復と、脱離したSOxに起因した白煙の発生の抑制とを高水準で両立させることのできるエンジン制御装置を提供することを目的とする。 The present invention has been made in view of the above-described problems. In the temperature rise control of the catalyst, the amount of SOx desorbed from the catalyst is accurately estimated, thereby recovering the function of the catalyst by desorbing SOx. An object of the present invention is to provide an engine control device that can achieve a high level of suppression of the generation of white smoke caused by desorbed SOx.
本発明は、ディーゼルエンジンの排気管に設けられる浄化装置の温度を、前記浄化装置からSOxが脱離する温度域の目標温度まで上昇させる昇温制御を実行するエンジン制御装置を対象としている。エンジン制御装置は、所定の制御周期ごとに前記浄化装置の温度の代表値である代表温度を取得する温度取得手段と、前記浄化装置に流入するSOx量を流入SOx量として前記制御周期ごとに推定する流入SOx量推定手段と、前記流入SOx量と前記代表温度とを用いて、前記浄化装置から新たに脱離するSOx量を新規脱離SOx量として前記制御周期ごとに推定する新規脱離SOx量推定手段と、前記新規脱離SOx量を用いて、前記浄化装置の温度上昇中の各温度において前記浄化装置に最終的に吸着するSOx量を前記代表温度に関連付けたグラフとして表される最終吸着SOx分布を前記制御周期ごとに推定する最終吸着SOx分布推定手段と、前記最終吸着SOx分布と許容脱離SO3量とを用いて、前記目標温度を前記制御周期ごとに算出する目標温度算出手段と、を備えている。前記新規脱離SOx量推定手段は、前記代表温度の過去値から今回値への変化度合を表す補正値を算出する算出手段と、前記代表温度に前記補正値を反映させる補正を行うことにより補正後代表温度を算出する補正手段と、を含み、前記流入SOx量と前記補正後代表温度を用いて前記新規脱離SOx量を推定するように構成されている。また、前記補正手段は、前記変化度合の前記補正後代表温度への反映度合が、前記代表温度が大きくなるに従って単調非増加で小さくなるように構成されている。 The present invention is directed to an engine control device that performs temperature increase control for increasing the temperature of a purification device provided in an exhaust pipe of a diesel engine to a target temperature in a temperature range where SOx is desorbed from the purification device. The engine control device estimates a representative temperature that is a representative value of the temperature of the purification device every predetermined control cycle, and estimates an SOx amount flowing into the purification device as an inflow SOx amount for each control cycle. New desorbed SOx that estimates the newly desorbed SOx amount as the new desorbed SOx amount for each control cycle using the inflow SOx amount estimating means that performs the inflow SOx amount and the representative temperature Using the amount estimation means and the new desorbed SOx amount, the SOx amount finally adsorbed to the purification device at each temperature during the temperature rise of the purification device is represented as a graph related to the representative temperature. and final suction SOx distribution estimating means for estimating a suction SOx distribution for each of the control period, using said final suction SOx distribution and allowable desorption SO 3 content, the system of the target temperature And a target temperature calculation means for calculating for each cycle. The new desorption SOx amount estimation means is corrected by calculating a correction value that represents a degree of change of the representative temperature from the past value to the current value, and correcting the representative temperature to reflect the correction value. Correction means for calculating a post representative temperature, and configured to estimate the new desorption SOx quantity using the inflow SOx quantity and the post-correction representative temperature. The correction means is configured such that the degree of reflection of the degree of change on the corrected representative temperature decreases as the representative temperature increases, monotonically non-increasing.
本発明のエンジン制御装置において、前記補正手段は、前記代表温度が大きくなるに従って単調非増加で小さくなる補正係数を算出し、前記補正値に前記補正係数を乗じた値を前記代表温度に加算することにより前記補正後代表温度を算出するように構成されていてもよい。 In the engine control apparatus of the present invention, the correction means calculates a correction coefficient that decreases monotonically and does not increase as the representative temperature increases, and adds a value obtained by multiplying the correction value by the correction coefficient to the representative temperature. Thus, the corrected representative temperature may be calculated.
また、本発明のエンジン制御装置において、前記算出手段は、前記代表温度の今回値と前回値の差分値を前記補正値として算出するように構成されていてもよい。 In the engine control device of the present invention, the calculation unit may be configured to calculate a difference value between the current value and the previous value of the representative temperature as the correction value.
また、本発明のエンジン制御装置は、前記浄化装置の温度上昇中の各温度において前記浄化装置に吸着するSOx量を前記代表温度に関連付けたグラフとして表される吸着SOx分布と、前記浄化装置の温度上昇中の各温度において前記浄化装置に吸着するSOx最大量を前記代表温度に関連付けたグラフとして表される飽和SOx分布と、を用いて、前記浄化装置におけるSOx飽和率を前記制御周期ごとに推定するSOx飽和率推定手段と、前記流入SOx量と前記SOx飽和率とを用いて、前記浄化装置に流入して前記浄化装置に新たに吸着するSOx量を新規吸着SOx量として前記制御周期ごとに推定する新規吸着SOx量推定手段と、前記新規吸着SOx量を用いて、前記浄化装置に新たなSOxが吸着した後の前記吸着SOx分布を吸着後SOx分布として前記制御周期ごとに推定する吸着後SOx分布推定手段と、を備えることができる。この場合、前記新規脱離SOx量推定手段は、前記吸着後SOx分布と前記代表温度とを用いて、前記新規脱離SOx量を前記制御周期ごとに推定するように構成されてもよい。また、本発明のエンジン制御装置は、前記新規吸着SOx量を用いて、前記浄化装置に流入して前記浄化装置に吸着することなくすり抜けるSOx量をすり抜けSOx量として前記制御周期ごとに推定するすり抜けSOx量推定手段と、前記浄化装置においてSO3に転化するSO2の転化率と前記代表温度との関係を表した転化率マップと、今回の推定サイクルにおける前記代表温度と、前記すり抜けSOx量と、を用いて、前記浄化装置にSOxの状態で流入して前記浄化装置に吸着することなくすり抜けてSO3の状態で排出されるSO3量をすり抜けSO3量として前記制御周期ごとに推定するすり抜けSO3量推定手段と、サルフェート白煙に関する制約に相当する前記浄化装置の下流におけるSO3量と前記すり抜けSO3量とを用いて、前記浄化装置から脱離することが許されるSO3量を前記許容脱離SO3量として前記制御周期ごとに算出する許容脱離SO3量算出手段と、をさらに備えるように構成されていてもよい。 Further, the engine control device of the present invention includes an adsorption SOx distribution represented as a graph in which the amount of SOx adsorbed to the purification device at each temperature during the temperature increase of the purification device is associated with the representative temperature, The saturation SOx distribution represented as a graph in which the maximum amount of SOx adsorbed to the purification device at each temperature during the temperature increase is related to the representative temperature, and the SOx saturation rate in the purification device for each control cycle. Using the estimated SOx saturation rate estimation means, the inflow SOx amount and the SOx saturation rate, the SOx amount that flows into the purification device and is newly adsorbed to the purification device is defined as a new adsorption SOx amount for each control period. The new adsorption SOx amount estimating means for estimating the new adsorption SOx amount and the adsorption SOx after the new SOx is adsorbed to the purification device using the new adsorption SOx amount A post-adsorption SOx distribution estimating means for estimating for each of the control period as SOx distribution after adsorption fabric may comprise. In this case, the new desorption SOx amount estimation means may be configured to estimate the new desorption SOx amount for each control cycle using the post-adsorption SOx distribution and the representative temperature. Further, the engine control device of the present invention uses the new adsorption SOx amount to estimate the amount of SOx that flows into the purification device and slips through without being adsorbed by the purification device as the slipping SOx amount for each control cycle. SOx amount estimation means, a conversion rate map representing the relationship between the conversion rate of SO 2 converted to SO 3 in the purifier and the representative temperature, the representative temperature in the current estimation cycle, and the slip-by SOx amount , is used to estimate for each of the control period as the purifier slip through SO 3 amount discharged in the form of SO 3 to pass through without being adsorbed on the purifier flows into a state of SOx to SO 3 amount and SO 3 amount estimating means slipping, said the SO 3 content in the downstream of the purification device corresponding to constraints on sulfate white smoke slipped SO With the amount, and the allowable desorbed SO 3 amount calculating means for calculating an SO 3 amount that is allowed to detach from the purification device to each of the control period as the allowable desorbed SO 3 amount, further comprising as the It may be configured.
さらに、本発明のエンジン制御装置において、前記温度取得手段は、前記排気管における前記浄化装置の下流側に流れたガスの温度を前記代表温度として取得するように構成されていてもよい。 Furthermore, in the engine control apparatus of the present invention, the temperature acquisition means may be configured to acquire the temperature of the gas flowing downstream of the purification device in the exhaust pipe as the representative temperature.
また、本発明のエンジン制御装置において、前記浄化装置は、前記排気管を流れる微粒子を捕集するフィルタを含んでいてもよい。この場合、前記昇温制御を、前記フィルタに捕集された微粒子量の推定値が除去要求量に到達したときに開始することとしてもよい。 In the engine control device of the present invention, the purification device may include a filter that collects particulates flowing through the exhaust pipe. In this case, the temperature increase control may be started when the estimated value of the amount of fine particles collected by the filter reaches the removal request amount.
本発明によれば、浄化装置の代表温度を用いて、浄化装置から新たに脱離する脱離SOx量が推定される。この際、代表温度の過去値から今回値への変化度合を表す補正値を代表温度に反映させる補正が行なわれ、補正された補正後代表温度を用いて新規脱離SOx量が推定される。ここで、代表温度の補正は、代表温度が大きくなるに従い補正値の代表温度への反映度合が単調非増加で小さくなるように行なわれる。脱離SOx量は浄化装置の代表温度だけでなく、代表温度の変化度合の影響も受ける。具体的には、浄化装置の低温域では代表温度の変化度合の影響を大きく受けるのに対し、浄化装置の高温域では低温域ほどに影響を受けない。このため、本発明によれば、浄化装置の温度が高温であっても脱離SOx量を精度よく推定することができるので、触媒の昇温制御において、SOxの脱離による触媒の機能の回復と、このSOxに起因した白煙の発生の抑制とを高水準で両立させることができる。 According to the present invention, the amount of desorbed SOx newly desorbed from the purification device is estimated using the representative temperature of the purification device. At this time, a correction value that reflects the degree of change of the representative temperature from the past value to the current value is reflected in the representative temperature, and a new desorbed SOx amount is estimated using the corrected representative temperature after correction. Here, the correction of the representative temperature is performed so that the degree of reflection of the correction value to the representative temperature decreases with a monotonous non-increase as the representative temperature increases. The amount of desorbed SOx is affected not only by the representative temperature of the purification device but also by the degree of change in the representative temperature. Specifically, the temperature of the purification device is greatly affected by the degree of change in the representative temperature, whereas the temperature of the purification device is not affected as much as the low temperature region. Therefore, according to the present invention, since the amount of desorbed SOx can be accurately estimated even when the temperature of the purifier is high, recovery of the function of the catalyst by desorption of SOx can be performed in the temperature rise control of the catalyst. And the suppression of the generation of white smoke due to SOx can be achieved at a high level.
以下、図面を参照して本発明の実施の形態について説明する。ただし、以下に示す実施の形態において各要素の個数、数量、量、範囲等の数に言及した場合、特に明示した場合や原理的に明らかにその数に特定される場合を除いて、その言及した数に、この発明が限定されるものではない。また、以下に示す実施の形態において説明する構造やステップ等は、特に明示した場合や明らかに原理的にそれに特定される場合を除いて、この発明に必ずしも必須のものではない。 Embodiments of the present invention will be described below with reference to the drawings. However, in the embodiment shown below, when referring to the number of each element, quantity, quantity, range, etc., unless otherwise specified or clearly specified in principle, the reference However, the present invention is not limited to these numbers. Further, the structures, steps, and the like described in the embodiments below are not necessarily essential to the present invention unless otherwise specified or clearly specified in principle.
実施の形態1.
本発明の実施の形態1について図を参照して説明する。
[システム構成の説明]
図1は、本発明の実施の形態のシステム構成を示す図である。図1に示すシステムは、車両に搭載されるディーゼルエンジン10(以下単に「エンジン10」ともいう。)を備えている。エンジン10の各気筒には、燃料としての軽油を噴射するインジェクタ12が設けられている。なお、図1に描かれるエンジン10は直列4気筒エンジンであるが、エンジン10の気筒数および気筒配列は特に限定されない。また、図1には、4つのインジェクタ12のうちの1つが描かれている。
[Description of system configuration]
FIG. 1 is a diagram showing a system configuration according to an embodiment of the present invention. The system shown in FIG. 1 includes a diesel engine 10 (hereinafter also simply referred to as “
エンジン10の排気マニホールド14には、ターボチャージャ16の排気タービン16aの入口が接続されている。排気タービン16aは吸気管18に設けられたコンプレッサ16bに連結されている。コンプレッサ16bは、排気タービン16aの回転により駆動して吸気を過給する。排気タービン16aの出口には排気管20が接続されている。排気管20には排気浄化装置22が設けられている。排気浄化装置22は、DOC(Diesel Oxidation Catalyst)22aと、DPF(Diesel Particulate Filter)22bと、を備えている。DOC22aは、排気中の炭化水素(HC)や一酸化炭素(CO)を酸化して、水(H2O)や二酸化炭素(CO2)に転化させる機能を有する触媒である。DPF22bは、排気中に含まれる微粒子(PM)を捕集するフィルタである。排気浄化装置22の上流には、インジェクタ12と共通の燃料を排気管20に添加する燃料添加弁24が設けられている。
The
図1に示すシステムは、制御装置としてのECU(Electronic Control Unit)30を備えている。ECU30は、RAM(Random Access Memory)、ROM(Read Only Memory)、CPU(Central Processing Unit)等を備えている。ECU30は、車両に搭載された各種センサの信号を取り込み処理する。各種センサには、吸気管18の入口付近に設けられたエアフローメータ32、DOC22aの出口温度を検出する温度センサ34、DPF22bの上下流における圧力差を検出する差圧センサ36が含まれている。ECU30は、取り込んだ各センサの信号を処理して所定の制御プログラムに従って各種アクチュエータを操作する。ECU30によって操作されるアクチュエータには、上述したインジェクタ12、燃料添加弁24が含まれている。
The system shown in FIG. 1 includes an ECU (Electronic Control Unit) 30 as a control device. The
[DPF22bの再生制御]
本実施の形態では、ECU30によるエンジン制御として、DPF22bの昇温制御(以下「PM再生制御」ともいう。)が行われる。PM再生制御は、DPF22bで捕集したPMの推定値が除去要求量に到達したときに、燃料添加弁24から燃料を添加する制御である。例えば差圧センサ36で検出した圧力差が所定値に到達したときに、PMの推定値が除去要求量に到達したと判断することができる。燃料添加弁24から燃料を添加することでDOC22aにおいて添加燃料を酸化し、この酸化反応熱によりDPF22bの床温を600℃以上まで上昇させる。これにより、DPF22bで捕集したPMを燃焼除去できるので、DPF22bの捕集機能を回復させることができる。なお、DPF22bの床温を600℃以上に上昇させるための燃料添加弁24からの添加燃料量(以下「DPF用燃料量)は、DPFの床温と関連付けたマップに基づき決定されるものとする。このようなマップは例えばECU30のROMに記憶させておき、DPF22bの実際の床温に応じて適宜読み出すことができる。
[Regeneration control of DPF 22b]
In the present embodiment, as the engine control by the
[PM再生制御における問題点]
ところで、ディーゼルエンジンの燃料や潤滑油には一般に硫黄が含まれており、燃料の燃焼に伴ってこうした硫黄からSOxが生成する。本実施の形態においても同様で、エンジン10での燃料の燃焼に伴ってSOxが生成する。生成したSOxはエンジン10から排出されて排気浄化装置22に流入し、主にDOC22aに吸着する。但し、DOC22aの床温が高くなると、ここに吸着していたSOxが脱離し始める。DOC22aの組成等により多少変動するものの、PM再生制御を行う温度域ではDOC22aからSOxが脱離して下流側に放出される。
[Problems in PM regeneration control]
By the way, the fuel and lubricating oil of diesel engines generally contain sulfur, and SOx is generated from such sulfur as the fuel burns. The same applies to the present embodiment, and SOx is generated with the combustion of fuel in the
DOC22aにおけるSOxの吸着と脱離について、図2を参照しながら説明する。図2は、DOC22aにおけるSOxの吸着と脱離を説明するための図である。この図に示すように、DOC22aは、基材(不図示)の表面を覆うコート材22cと、貴金属22d(Pt,Pd等)とを備えている。貴金属22dはコート材22cに分散担持されており、HCやCOを酸化する際の活性点となる。但し、排気中のSO2が貴金属22dに吸着し、または、排気中のSO3がコート材22cに吸着する。貴金属22dに吸着したSO2の一部は貴金属22dから脱離して排気中に戻り、または、貴金属22d上で酸化されてSO3となり、SO3の状態でコート材22cに吸着する。即ち、貴金属22dにはSO2が吸着し、コート材22cには排気由来のSO3とSO2由来のSO3が吸着する。何れにせよ、SOxが吸着することで、DOC22aにおけるHC等の酸化機能が阻害されることになる。
The adsorption and desorption of SOx in the
上述した2経路によりコート材22cに吸着したSO3は、コート材22cの床温が高くなることで脱離する。また、コート材22cの床温が高くなることで貴金属22d上でのSO2からSO3への転化が促進されるので、このようなSO3もコート材22cから脱離する。従って、PM再生制御を行うことで、上述したDPF22bの捕集機能だけでなく、DOC22aにおけるHC等の酸化機能を回復させることもできる。ところが、図2に示すように、コート材22cから脱離したSO3が排気管20に存在するH2Oと反応することでH2SO4が発生する。そして、このH2SO4の濃度が一定濃度を超えると視認可能な白煙(サルフェート白煙)となるので、エンジン10を搭載した車両の商品価値を損ねてしまうおそれがある。
The SO 3 adsorbed on the
[本実施の形態の特徴]
DOC22aの下流における排気中のH2SO4の濃度が高くなり過ぎないように燃料添加弁24から燃料を添加すれば、PM再生制御中のサルフェート白煙の発生を抑制できる。そこで本実施の形態では、DOC22aの下流におけるSO3の濃度がサルフェート白煙に関する制約を満たすように、PM再生制御中のDOC22aの床温の目標温度(以下「目標床温Ttrg」ともいう。)を算出し、目標床温Ttrgに基づいて、燃料添加弁24から添加する燃料量(以下「制約充足用燃料量」ともいう。)を算出することとしている。なお、このような制約SO3濃度(DOC22aの下流におけるSO3濃度の上限値)は、例えばECU30のROMに記憶させておくことができる。制約充足用燃料量よりもDPF用燃料量の方が多い場合に、DPF用燃料量ではなく制約充足用燃料量を採用することで、サルフェート白煙に関する制約を満たしつつ、DOC22aにおけるHC等の酸化機能を回復させることができる。
[Features of this embodiment]
If fuel is added from the
[目標床温Ttrgの算出ロジック]
図3は、目標床温Ttrgを算出するためのロジックを示す機能ブロック図であり、これはECU30により実現される。この図に示すように、ECU30は、流入SOx量推定部M1と、SOx飽和率推定部M2と、新規吸着SOx量およびすり抜けSOx量推定部M3と、吸着後SOx分布推定部M4と、新規脱離SOx量推定部M5と、最終吸着SOx分布推定部M6と、すり抜けSO3量推定部M7と、許容脱離SO3量算出部M8と、白煙抑制目標床温算出部M9と、を備えており、これらの要素M1〜M9によって制御周期ごとに(具体的にはエンジン10の燃焼サイクルごとに)目標床温Ttrgを算出するものとする。なお、以下の説明においては、要素M1〜M9を簡略化するものとし、例えば流入SOx量推定部M1を「推定部M1」ともいう。
[Calculation logic of target bed temperature Ttrg]
FIG. 3 is a functional block diagram showing logic for calculating the target bed temperature Ttrg, and this is realized by the
推定部M1は、DOC22aに流入するSOxの量(以下「流入SOx量」ともいう。)を推定する。なお、本明細書でいう「DOC22aに流入するSOx」には、エンジン10で生成し、ここから排出されてDOC22aに流入するSOxだけでなく、燃料添加弁24から添加された燃料のDOC22aでの酸化反応に伴って生成し、DOC22a上を流れるSOxも含まれるものとする。
The estimation unit M1 estimates the amount of SOx flowing into the
推定部M1は、具体的に、インジェクタ12からの噴射燃料量(筒内噴射量)および燃料添加弁24からの添加燃料量(排気添加量)を変数とする次式(1)により、第t番目のサイクルにおける流入SOx量を推定する。なお、式(1)の燃料S濃度は、燃料中の硫黄濃度であり、燃料供給系に別途設けた硫黄濃度センサの検出値を用いてもよく、設定値を用いてもよい。
流入SOx量(排気添加量(t),筒内噴射量(t))[μg/s]=流入燃料量(排気添加量(t),筒内噴射量(t))[g/s]×燃料S濃度[ppm] ・・・(1)
Specifically, the estimation unit M1 calculates the t th by the following equation (1) using the injected fuel amount (in-cylinder injection amount) from the
Inflow SOx amount (exhaust addition amount (t), in-cylinder injection amount (t)) [μg / s] = inflow fuel amount (exhaust addition amount (t), in-cylinder injection amount (t)) [g / s] × Fuel S concentration [ppm] (1)
式(1)の流入燃料量(排気添加量(t),筒内噴射量(t))は、「DOC22aに流入するSOx」の由来となった燃料の第t番目のサイクルにおける量であり、燃料の比重(軽油比重)を用いて次式(2)により算出される。
流入燃料量(排気添加量(t),筒内噴射量(t))[g/s]=(排気添加量(t)[g/s]÷1000×軽油比重[g/cm3]+筒内噴射量(t)[g/s]) ・・・(2)
なお、以下の説明においては、流入SOx量(排気添加量(t),筒内噴射量(t))を流入SOx量(t)ともいう。また、流入燃料量(排気添加量(t),筒内噴射量(t))を流入燃料量(t)ともいう。
The inflow fuel amount (exhaust gas addition amount (t), in-cylinder injection amount (t)) in equation (1) is the amount of fuel in the t-th cycle that originated from “SOx flowing into
Inflow fuel amount (exhaust addition amount (t), in-cylinder injection amount (t)) [g / s] = (exhaust addition amount (t) [g / s] ÷ 1000 × light oil specific gravity [g / cm3] + in-cylinder Injection amount (t) [g / s]) (2)
In the following description, the inflow SOx amount (exhaust addition amount (t), in-cylinder injection amount (t)) is also referred to as inflow SOx amount (t). The inflow fuel amount (exhaust gas addition amount (t), in-cylinder injection amount (t)) is also referred to as inflow fuel amount (t).
推定部M2は、DOC22aにおけるSOxの飽和率(以下「SOx飽和率」ともいう。)を推定する。SOx飽和率の推定には、DOC22aの床温上昇中の各床温においてDOC22aに吸着するSOxの量(以下「吸着SOx量」ともいう。)を、DOC22aの床温に関連付けたグラフとして表される分布(以下「吸着SOx分布」ともいう。)と、DOC22aの床温上昇中の各床温においてDOC22aに吸着するSOxの最大量(以下「飽和SOx量」ともいう。)を、DOC22aの床温に関連付けたグラフとして表される分布(以下「飽和SOx分布」ともいう。)と、が用いられる。先ず、吸着SOx分布と飽和SOx分布について、SO3を例として図4を参照しながら説明する。
The estimation unit M2 estimates the SOx saturation rate (hereinafter also referred to as “SOx saturation rate”) in the
図4は、吸着SOx分布と飽和SOx分布を説明するための図である。図4に「吸着SO3量」として示したデータは、次の手法により収集したものである。具体的には先ず、図4に「現在温度」として示した床温において、十分な量のSOxをDOC22aに吸着させる。続いて、DOC22aの床温上昇中の各床温においてDOC22aから脱離したSO3の量を、上昇速度を一定とした条件のもとで測定する。そして、この脱離SO3量をDOC22aの床温に関連付けてグラフを作成する。これにより、脱離SO3量を表した分布(以下「脱離SO3分布」ともいう。)を得ることができる。これと同様の手法により、DOC22aの床温上昇中の各床温においてDOC22aから脱離したSO2の量を、DOC22aの床温に関連付けたグラフ(以下「脱離SO2分布」ともいう。)を得ることもできる。なお、DOC22aから脱離するSO3については、これをセンサで直接的に測定してもよいし、SOxまたはSO2を検出するセンサを用いて両者を測定し、これらの差から算出してもよい(SO3=SOx−SO2)。
FIG. 4 is a diagram for explaining the adsorption SOx distribution and the saturated SOx distribution. The data shown as “adsorbed SO 3 amount” in FIG. 4 is collected by the following method. Specifically, first, a sufficient amount of SOx is adsorbed on the
ここで、DOC22aの床温上昇中にDOC22aから脱離するSO3は、実際には図4に「現在温度」として示した床温でDOC22aに吸着させたSO3である。しかし、ある床温でDOC22aから脱離するSO3は、その床温に至るまでDOC22aに吸着し続けることのできたSO3であり、更に言えば、その床温でDOC22aに吸着することができたSO3であると考えることもできる。このような仮定に基づいて、上述した脱離SO3分布の縦軸を、DOC22aの床温上昇中の各床温においてDOC22aに吸着するSO3の量に置き換えると、図4に示した「吸着SO3量」のデータのグラフ、即ち、吸着SO3分布を得ることができる。そして、これと同様の手法によって、吸着SO2分布を得ることもできる。
Here, SO 3 desorbed from DOC22a in increased bed temperature of DOC22a is SO 3 which actually adsorbed to DOC22a in bed temperature shown as "current temperature" in FIG. However, SO 3 desorbed from DOC22a at some bed temperature, the bed is SO 3 which could continue to adsorb to DOC22a up to temperature and, more could be adsorbed to DOC22a in that bed temperature It can also be considered SO 3 . Based on this assumption, if the vertical axis of the desorption SO 3 distribution described above is replaced with the amount of SO 3 adsorbed on the
また、図4に「飽和SO3量」として示したデータは、「吸着SO3量」のデータと同様の手法により収集したものである。この「飽和SO3量」のデータは、具体的に、上昇速度を極低速とした条件のもと、DOC22aの床温上昇中の各床温(例えば5℃間隔)においてDOC22aから脱離したSO3の量に相当している。DOC22aの床温の上昇速度が極低速であることから、この「飽和SO3量」のデータは、DOC22aから脱離したSO3の量の最大値であると考えることができる。また、この最大値に対しては、上述した仮定を適用することができる。即ち、ある床温でDOC22aから脱離するSO3の最大量は、その床温でDOC22aに吸着することができたSO3の最大量に等しいと考えることができる。このような仮定に基づいて上述した脱離SO3分布の縦軸を上記SO3の最大量に置き換えると、図4に示した「飽和SO3量」のデータのグラフ、即ち、飽和SO3分布を得ることができる。そして、これと同様の手法によって、飽和SO2分布を得ることもできる。
Further, the data shown as “saturated SO 3 amount” in FIG. 4 is collected by the same method as the “adsorption SO 3 amount” data. Specifically, this “saturated SO 3 amount” data is obtained from the SO 2 desorbed from the
推定部M2は、第t番目のサイクルにおけるDOC22aの現在床温T2を変数とする次式(3)により、第t番目のサイクルにおけるSOx飽和率(T2(t),t)を推定する。なお、現在床温T2は、DOC22aの現在の床温の代表温度として、例えば温度センサ34の検出値を使用することができる。
SOx飽和率(T2(t),t)=1−(総吸着余裕量(T2(t),t)/総飽和量(T2(t),t)) ・・・(3)
The estimation unit M2 estimates the SOx saturation rate (T 2 (t), t) in the t-th cycle by the following equation (3) using the current bed temperature T 2 of the
SOx saturation rate (T 2 (t), t) = 1− (total adsorption margin amount (T 2 (t), t) / total saturation amount (T 2 (t), t)) (3)
式(3)のSOx飽和率(T2(t),t)の算出過程は次のとおりである。先ず、DOC22aの床温上昇中における床温T1と、現在床温T2とを変数とする次式(4)および(5)により、第t番目のサイクルにおける飽和SO2分布(T1,T2(t),t)および飽和SO3分布(T1,T2(t),t)をそれぞれ算出する。
飽和SO2分布(T1,T2(t),t)[μg/℃]=基準飽和SO2分布×床温補正SO2map(T2(t))[μg/℃] ・・・(4)
飽和SO3分布(T1,T2(t),t)[μg/℃]=基準飽和SO3分布×床温補正SO3map(T2(t))[μg/℃] ・・・(5)
The calculation process of the SOx saturation rate (T 2 (t), t) in Expression (3) is as follows. First, the saturated SO 2 distribution (T 1 , T 1 , T 2 ) in the t-th cycle is calculated by the following equations (4) and (5) using the bed temperature T 1 during the rise of the bed temperature of the
Saturated SO 2 distribution (T 1 , T 2 (t), t) [μg / ° C.] = Standard saturated SO 2 distribution × bed temperature corrected SO 2 map (T 2 (t)) [μg / ° C.] 4)
Saturated SO 3 distribution (T 1 , T 2 (t), t) [μg / ° C.] = Reference saturated SO 3 distribution × bed temperature corrected SO 3 map (T 2 (t)) [μg / ° C.] 5)
式(4)の基準飽和SO2分布は、十分な量のSOxをDOC22aに吸着させるときの床温(図4の「現在温度」)を基準床温(例えば上述した吸着限界量が最大となる300℃付近の床温)として作成した飽和SO2分布である。式(5)の基準飽和SO3分布もこれと同様である。式(4)の床温補正SO2map(T2(t))は、基準飽和SO3分布を現在床温T2の飽和SO2分布に変換するための補正値を定めたマップである。式(5)の床温補正SO3map(T2(t))もこれと同様である。このような基準飽和SOx分布と補正マップは、例えばECU30のROMに記憶させておくことができ、現在床温T2に応じて適宜読み出すことができる。
In the reference saturated SO 2 distribution of the equation (4), the bed temperature (“current temperature” in FIG. 4) when a sufficient amount of SOx is adsorbed by the
基準飽和SOx分布と補正後の飽和SOx分布の関係を、SO2を例として図5を参照しながら説明する。図5は、基準飽和SOx分布と補正後の飽和SOx分布の関係を説明するための図である。なお、この図の横軸のTLおよびTHは、DOC22aの床温上昇中にDOC22aからSO2が脱離し始める温度(下限温度)と、DOC22aからSO2が脱離し終わる温度(上限温度)に、それぞれ相当している。この図に示す3種類の分布の違いは、現在床温T2にある。即ち、現在床温T2が基準温度と等しい場合は、補正後の飽和SO2分布の形状が基準飽和SO2分布の形状と一致する(中央)。一方、現在床温T2が基準温度よりも低い場合(左方)や、現在床温T2が基準温度よりも高い場合(右方)は、補正後の飽和SO2分布の形状が基準飽和SO2分布の形状と一致しなくなる。なお、現在床温T2が基準温度よりも高い場合(右方)には、補正後の飽和SO2分布の形状が、現在床温T2よりも低温側のデータが欠落しているような形状となる。この理由は、現在床温T2よりも低温側では、本来であればこの床温域においてDOC22aに吸着し続けることのできたはずのSOxが、DOC22aから既に脱離していると考えられるためである。
The relationship between the reference saturated SOx distribution and the corrected saturated SOx distribution will be described with reference to FIG. 5 using SO 2 as an example. FIG. 5 is a diagram for explaining the relationship between the reference saturated SOx distribution and the corrected saturated SOx distribution. Note that TL and TH on the horizontal axis in this figure are the temperature at which SO 2 begins to desorb from the
続いて、式(4)により算出した飽和SO2分布(T1,T2(t),t)を次式(6)に代入して、第t番目のサイクルにおける総飽和SO2量(T2(t),t)を算出する。また、式(5)により算出した飽和SO3分布(T1,T2(t),t)を次式(7)に代入して、第t番目のサイクルにおける総飽和SO3量を算出する。 Subsequently, the saturated SO 2 distribution (T 1 , T 2 (t), t) calculated by the equation (4) is substituted into the following equation (6), and the total saturated SO 2 amount (T 2 (t), t) are calculated. Further, the saturated SO 3 distribution (T 1 , T 2 (t), t) calculated by the equation (5) is substituted into the following equation (7) to calculate the total saturated SO 3 amount in the t-th cycle. .
総飽和SO2量(T2(t),t)および総飽和SO3量(T2(t),t)を算出したら、次式(8)にこれらを代入して、第t番目のサイクルにおける総飽和量(T2(t),t)を算出する。
総飽和量(T2(t),t)=総飽和SO2量(T2(t),t)+総飽和SO3量(T2(t),t) ・・・(8)
なお、以下の説明においては、総飽和SO2量(T2(t),t)を単に総飽和SO2量(t)ともいう。また、総飽和SO3量(T2(t),t)を単に総飽和SO3量(t)ともいう。また、総飽和量(T2(t),t)を単に総飽和量(t)ともいう。
When the total saturated SO 2 amount (T 2 (t), t) and the total saturated SO 3 amount (T 2 (t), t) are calculated, these are substituted into the following equation (8), and the t-th cycle The total saturation amount at (T 2 (t), t) is calculated.
Total saturation amount (T 2 (t), t) = Total saturation SO 2 amount (T 2 (t), t) + Total saturation SO 3 amount (T 2 (t), t) (8)
In the following description, the total saturated SO 2 amount (T 2 (t), t) is also simply referred to as the total saturated SO 2 amount (t). The total saturated SO 3 amount (T 2 (t), t) is also simply referred to as the total saturated SO 3 amount (t). The total saturation amount (T 2 (t), t) is also simply referred to as the total saturation amount (t).
式(8)により総飽和量(t)を算出したら、飽和SO2分布(T1,T2(t),t)と、推定部M6で推定した第t番目のサイクルにおける最終吸着SO2分布(T1,t)とを次式(9)に代入して、第t番目のサイクルにおける吸着余裕SO2分布(T1,T2(t),t)を算出する。また、飽和SO3分布(T1,T2(t),t)と、推定部M6で推定した第t番目のサイクルにおける最終吸着SO3分布(T1,t)とを次式(10)に代入して、第t番目のサイクルにおける吸着余裕SO3分布(T1,T2(t),t)を算出する。
吸着余裕SO2分布(T1,T2(t),t)[μg/℃]=max{飽和SO2分布(T1,T2(t),t)[μg/℃]−最終吸着SO2分布(T1,t)[μg/℃],0} ・・・(9)
吸着余裕SO3分布(T1,T2(t),t)[μg/℃]=max{飽和SO3分布(T1,T2(t),t)[μg/℃]−最終吸着SO3分布(T1,t)[μg/℃],0} ・・・(10)
When the total saturation amount (t) is calculated by the equation (8), the saturated SO 2 distribution (T 1 , T 2 (t), t) and the final adsorption SO 2 distribution in the t-th cycle estimated by the estimation unit M6. Substituting (T 1 , t) into the following equation (9), the adsorption margin SO 2 distribution (T 1 , T 2 (t), t) in the t-th cycle is calculated. Further, the saturated SO 3 distribution (T 1 , T 2 (t), t) and the final adsorption SO 3 distribution (T 1 , t) in the t-th cycle estimated by the estimation unit M6 are expressed by the following equation (10). And the adsorption margin SO 3 distribution (T 1 , T 2 (t), t) in the t-th cycle is calculated.
Adsorption margin SO 2 distribution (T 1 , T 2 (t), t) [μg / ° C.] = Max {saturated SO 2 distribution (T 1 , T 2 (t), t) [μg / ° C.] − Final adsorption SO 2 distribution (T 1 , t) [μg / ° C.], 0} (9)
Adsorption margin SO 3 distribution (T 1 , T 2 (t), t) [μg / ° C.] = Max {saturated SO 3 distribution (T 1 , T 2 (t), t) [μg / ° C.] − Final adsorption SO 3 distribution (T 1 , t) [μg / ° C.], 0} (10)
続いて、式(9)により算出した吸着余裕SO2分布(T1,T2(t),t)を次式(11)に代入して、第t番目のサイクルにおける総吸着余裕SO2量(T2(t),t)を算出する。また、式(10)により算出した吸着余裕SO3分布(T1,T2(t),t)を次式(12)に代入して、第t番目のサイクルにおける総吸着余裕SO3量(T2(t),t)を算出する。 Subsequently, the adsorption margin SO 2 distribution (T 1 , T 2 (t), t) calculated by the equation (9) is substituted into the following equation (11), and the total adsorption margin SO 2 amount in the t-th cycle is calculated. (T 2 (t), t) is calculated. Further, the adsorption margin SO 3 distribution (T 1 , T 2 (t), t) calculated by the equation (10) is substituted into the following equation (12), and the total adsorption margin SO 3 amount in the t-th cycle ( T 2 (t), t) is calculated.
なお、以下の説明においては、吸着余裕SO2分布(T1,T2(t),t)を単に吸着余裕SO2分布(t)ともいう。また、吸着余裕SO3分布(T1,T2(t),t)を単に吸着余裕SO3分布(t)ともいう。また、総吸着余裕SO2量(T2(t),t)を単に総吸着余裕SO2量(t)ともいう。また、総吸着余裕SO3量(T2(t),t)を単に総吸着余裕SO3量(t)ともいう。 In the following description, the adsorption margin SO 2 distribution (T 1 , T 2 (t), t) is also simply referred to as an adsorption margin SO 2 distribution (t). The adsorption margin SO 3 distribution (T 1 , T 2 (t), t) is also simply referred to as an adsorption margin SO 3 distribution (t). The total adsorption margin SO 2 amount (T 2 (t), t) is also simply referred to as the total adsorption margin SO 2 amount (t). The total adsorption margin SO 3 amount (T 2 (t), t) is also simply referred to as the total adsorption margin SO 3 amount (t).
図6を参照して、総吸着余裕SO2量(t)を説明する。図6は、総吸着余裕SO2量を説明するための図である。なお、総吸着余裕SO3量(t)についてはこれと同様である。この図に示すように、総吸着余裕SO2量(t)は、飽和SO2分布から、飽和SO2分布と吸着SO2分布の重複部分を除いた面積として表すことができる。なお、この図の右方の分布に領域Aとして示すように、DOC22aの床温上昇中の各床温においてDOC22aに吸着するSO2の量、即ち、吸着SO2量がその最大量、即ち、飽和SO2量を上回る場合は、DOC22aが飽和していると考えられることから、総吸着余裕SO2量(t)の算出から除外される。また、この右方の分布において、現在床温T2よりも低温側のデータが欠落している理由については、図5で説明した通りである。
The total adsorption margin SO 2 amount (t) will be described with reference to FIG. FIG. 6 is a diagram for explaining the total adsorption margin SO 2 amount. The total adsorption margin SO 3 amount (t) is the same as this. As shown in this figure, the total adsorption margin SO 2 amount (t) may be a saturated SO 2 distribution, expressed as the area excluding the overlapping portions of the adsorption SO 2 distribution and saturated SO 2 distribution. In addition, as shown as a region A in the distribution on the right side of this figure, the amount of SO 2 adsorbed on the
そして、総吸着余裕SO2量(t)および総吸着余裕SO3量(t)を算出したら、これらを次式(13)に代入して、第t番目のサイクルにおける総吸着余裕量(T2(t),t)を算出する。
総吸着余裕量(T2(t),t)=総吸着余裕SO2量(t)+総吸着余裕SO3量(t) ・・・(13)
更に、式(8)により算出した総飽和量(t)と、式(13)により算出した総吸着余裕量(t)とを式(3)に代入すれば、飽和率(T2(t),t)を算出できる。なお、以下の説明においては、飽和率(T2(t),t)を単に飽和率(t)ともいう。
When the total adsorption margin SO 2 amount (t) and the total adsorption margin SO 3 amount (t) are calculated, these are substituted into the following equation (13) to obtain the total adsorption margin amount (T 2 ) in the t-th cycle. (T), t) are calculated.
Total adsorption margin (T 2 (t), t) = Total adsorption margin SO 2 amount (t) + Total adsorption margin SO 3 amount (t) (13)
Furthermore, if the total saturation amount (t) calculated by the equation (8) and the total adsorption margin amount (t) calculated by the equation (13) are substituted into the equation (3), the saturation rate (T 2 (t) , T) can be calculated. In the following description, the saturation rate (T 2 (t), t) is also simply referred to as the saturation rate (t).
図3に戻り目標床温Ttrgの算出ロジックの説明を続ける。推定部M3は、「DOC22aに流入するSOx」であってDOC22aに新たに吸着するSOxの量(以下「新規吸着SOx量」ともいう。)、および、「DOC22aに流入するSOx」であってDOC22aに吸着することなくすり抜けるSOxの量(以下「すり抜けSOx量」ともいう。)を推定する。先ず、新規吸着SOx量とすり抜けSOx量の関係について、図7を参照して説明する。図7は、新規吸着SOx量とすり抜けSOx量の関係を説明するための図である。この図に矢印で示すように、新規吸着SOx量とすり抜けSOx量の和が、流入SOx量に等しくなる。この理由は、「DOC22aに流入するSOx」のうちの一部がDOC22aに吸着し、残りがDOC22aに吸着することなくすり抜けるためである。
Returning to FIG. 3, the description of the calculation logic of the target bed temperature Ttrg will be continued. The estimation unit M3 is “SOx flowing into the
推定部M3は、具体的に、推定部M1で推定した流入SOx量(t)と、推定部M2で推定した飽和率(t)とを変数とする次式(14)により新規吸着SOx量を推定し、次式(15)によりすり抜けSOx量を推定する。
新規吸着SOx量(流入SOx量(t),飽和率(t))[μg/s]=流入SOx量(t)×吸着率map(飽和率(t)) ・・・(14)
すり抜けSOx量(流入SOx量(t),飽和率(t))[μg/s]=流入SOx量(t)×{1−吸着率map(飽和率(t))} ・・・(15)
なお、以下の説明においては、新規吸着SOx量(流入SOx量(t),飽和率(t))を単に新規吸着SOx量(t)ともいう。また、すり抜けSOx量(流入SOx量(t),飽和率(t))を単にすり抜けSOx量(t)ともいう。
Specifically, the estimation unit M3 calculates the new adsorption SOx amount by the following equation (14) using the inflow SOx amount (t) estimated by the estimation unit M1 and the saturation rate (t) estimated by the estimation unit M2 as variables. The estimated SOx amount is estimated by the following equation (15).
New adsorption SOx amount (inflow SOx amount (t), saturation rate (t)) [μg / s] = inflow SOx amount (t) × adsorption rate map (saturation rate (t)) (14)
Slip-through SOx amount (inflow SOx amount (t), saturation rate (t)) [μg / s] = inflow SOx amount (t) × {1−adsorption rate map (saturation rate (t))} (15)
In the following description, the new adsorption SOx amount (inflow SOx amount (t), saturation rate (t)) is also simply referred to as a new adsorption SOx amount (t). The slip-through SOx amount (inflow SOx amount (t), saturation rate (t)) is also simply referred to as a slip-through SOx amount (t).
式(14)および(15)の吸着率mapは、第t番目のサイクルにおいて「DOC22aに流入するSOx」のうち、DOC22aに吸着するSOxの割合(即ち、吸着率)が、飽和率(t)によって変わるという特性に基づいて作成されたマップである。図8は、吸着率mapを説明するための図である。この図に示すように、吸着率mapの特性は、飽和率(t)が低い領域では吸着率が高く、飽和率(t)が高くなるにつれて吸着率が徐々に低下するというものである。なお、このようなマップは、例えばECU30のROMに記憶させておくことができ、現在床温T2に応じて適宜読み出すことができる。
The adsorption rate map of the formulas (14) and (15) is the saturation rate (t) in the ratio of SOx adsorbed to the
図3に戻り、推定部M4は、推定部M3で推定した新規吸着SOx量を吸着SOx分布に反映させて吸着後SOx分布を推定する。吸着後SOx分布について、SO2を例とした図9を参照しながら説明する。図9は、吸着後SOx分布を説明するための図である。この図に示すように、吸着後SO2分布は、前回サイクル(例えば、第t−1番目のサイクル)における最終吸着SO2分布に、今回サイクル(例えば、第t番目のサイクル)においてDOC22aに新たに吸着するSO2の量を表した分布(以下「新規吸着SO2分布」ともいう。)を加算することで推定される。
Returning to FIG. 3, the estimation unit M4 estimates the post-adsorption SOx distribution by reflecting the amount of new adsorption SOx estimated by the estimation unit M3 in the adsorption SOx distribution. The SOx distribution after adsorption will be described with reference to FIG. 9 using SO 2 as an example. FIG. 9 is a diagram for explaining the post-adsorption SOx distribution. As shown in this figure, the post-adsorption SO 2 distribution is updated to the final adsorption SO 2 distribution in the previous cycle (for example, the (t−1) th cycle), and to the
推定部M4は、具体的には先ず、新規吸着SOx量(t)、総吸着余裕量(t)および吸着余裕SO2分布(t)を変数とする次式(16)により、第t番目のサイクルにおける新規吸着SO2分布を算出する。新規吸着SO2分布と同様に、推定部M4は、DOC22aに新たに吸着するSO3の量を表した分布(以下「新規吸着SO3分布」ともいう。)を、次式(17)により算出する。なお、吸着余裕SO2分布(t)と総吸着余裕量(t)には、推定部M2において算出されたものが使用される。
新規吸着SO2分布(新規吸着SOx量(t),吸着余裕SO2分布(t),総吸着余裕量(t))[μg/℃]=吸着余裕SO2分布(t)[μg/℃]×{新規吸着SOx量(t)/総吸着余裕量(t)} ・・・(16)
新規吸着SO3分布(新規吸着SOx量(t),吸着余裕SO3分布(t),総吸着余裕量(t))[μg/℃]=吸着余裕SO3分布(t)[μg/℃]×{新規吸着SOx量(t)/総吸着余裕量(t)} ・・・(17)
なお、以下の説明においては、新規吸着SO2分布(新規吸着SOx量(t),吸着余裕SO2分布(t),総吸着余裕量(t))を、単に新規吸着SO2分布(t)ともいう。また、新規吸着SO3分布(新規吸着SOx量(t),吸着余裕SO3分布(t),総吸着余裕量(t))を、単に新規吸着SO3分布(t)ともいう。
Specifically, the estimation unit M4 firstly calculates the t-th number by the following equation (16) using the new adsorption SOx amount (t), the total adsorption margin amount (t), and the adsorption margin SO 2 distribution (t) as variables. Calculate the new adsorption SO 2 distribution in the cycle. Similar to the new adsorption SO 2 distribution, the estimation unit M4 calculates a distribution representing the amount of SO 3 newly adsorbed on the
New adsorption SO 2 distribution (new adsorption SOx amount (t), adsorption margin SO 2 distribution (t), total adsorption margin amount (t)) [μg / ° C.] = Adsorption margin SO 2 distribution (t) [μg / ° C.] × {New adsorption SOx amount (t) / Total adsorption margin amount (t)} (16)
New adsorption SO 3 distribution (new adsorption SOx amount (t), adsorption margin SO 3 distribution (t), total adsorption margin amount (t)) [μg / ° C.] = Adsorption margin SO 3 distribution (t) [μg / ° C.] × {New adsorption SOx amount (t) / Total adsorption margin amount (t)} (17)
In the following description, the new adsorption SO 2 distribution (new adsorption SOx amount (t), adsorption margin SO 2 distribution (t), total adsorption margin amount (t)) is simply referred to as new adsorption SO 2 distribution (t). Also called. The new adsorption SO 3 distribution (new adsorption SOx amount (t), adsorption margin SO 3 distribution (t), total adsorption margin amount (t)) is also simply referred to as new adsorption SO 3 distribution (t).
推定部M4は、続いて、算出した新規吸着SO2分布と、第t−1番目のサイクルにおける最終吸着SO2分布(t−1)とを次式(18)に代入して、吸着後SO2分布を算出する。また、算出した新規吸着SO3分布と、第t−1番目のサイクルにおいて推定部M6で推定した吸着SO3分布(t−1)とを次式(19)に代入して、吸着後SO3分布を算出する。
吸着後SO2分布(t)[μg/℃]=最終吸着SO2分布(t−1)[μg/℃]+新規吸着SO2分布(t)[μg/℃] ・・・(18)
吸着後SO3分布(t)[μg/℃]=最終吸着SO3分布(t−1)[μg/℃]+新規吸着SO3分布(t)[μg/℃] ・・・(19)
Subsequently, the estimation unit M4 substitutes the calculated new adsorption SO 2 distribution and the final adsorption SO 2 distribution (t−1) in the (t−1) -th cycle into the following equation (18) to obtain the post-adsorption SO 2 Two distributions are calculated. Further, the calculated new adsorption SO 3 distribution and the adsorption SO 3 distribution (t−1) estimated by the estimation unit M6 in the (t−1) -th cycle are substituted into the following equation (19), and the post-adsorption SO 3 Calculate the distribution.
Post-adsorption SO 2 distribution (t) [μg / ° C.] = Final adsorption SO 2 distribution (t−1) [μg / ° C.] + New adsorption SO 2 distribution (t) [μg / ° C.] (18)
Post-adsorption SO 3 distribution (t) [μg / ° C.] = Final adsorption SO 3 distribution (t−1) [μg / ° C.] + New adsorption SO 3 distribution (t) [μg / ° C.] (19)
図3に戻り、推定部M5は、推定部M4で推定した吸着後SOx分布に基づいて、DOC22aから新たに脱離するSOxの量(以下「新規脱離SOx量」ともいう。)を推定する。
Returning to FIG. 3, the estimation unit M5 estimates the amount of SOx newly desorbed from the
推定部M5は、具体的には先ず、DOC22aから脱離することのできるSOxの総量(以下「脱離可能総SOx量」ともいう。)を推定する。脱離可能総SOx量について、SO2を例とした図10を参照しながら説明する。図10は、脱離可能総SOx量を説明するための図である。なお、この図の横軸のTLおよびTHは、上述した下限温度および上限温度に、それぞれ相当している。この図に示すように、脱離可能総SOx量は、現在床温T2よりも低温側で、尚且つ、下限温度TLよりも高温側の吸着後SOx分布の面積に相当する。
Specifically, the estimation unit M5 first estimates the total amount of SOx that can be desorbed from the
DOC22aから脱離することのできるSO2の総量、即ち、脱離可能総SO2量は、現在床温T2を変数とする次式(20)により算出される。DOC22aから脱離することのできるSO3の総量、即ち、脱離可能総SO3量は、現在床温T2を変数とする次式(21)により算出される。
The total amount of SO 2 that can be desorbed from DOC22a, i.e., removable total SO 2 amount is calculated by the current equation of the bed temperature T 2 and a variable (20). The total amount of SO 3 that can be desorbed from the
推定部M5は、算出した脱離可能総SO2量を次式(22)に代入して、第t番目のサイクルにおいてDOC22aから新たに脱離するSO2の量、即ち、新規脱離SO2量を算出する。また、算出した脱離可能総SO3量を次式(23)に代入して、第t番目のサイクルにおいてDOC22aから新たに脱離するSO3の量、即ち、新規脱離SO3量を算出する。なお、式(22)および(23)の脱離率には設定値が使用され、例えばECU30のROMに記憶させておくことができる。
新規脱離SO2量(T2(t),t)[μg]=脱離可能総SO2量[μg]×脱離率 ・・・(22)
新規脱離SO3量(T2(t),t)[μg]=脱離可能総SO3量[μg]×脱離率 ・・・(23)
The estimation unit M5 substitutes the calculated total desorbable SO 2 amount into the following equation (22), and the amount of SO 2 newly desorbed from the
New desorbed SO 2 amount (T 2 (t), t) [μg] = total desorbable SO 2 amount [μg] × desorption rate (22)
New desorption SO 3 amount (T 2 (t), t) [μg] = total desorbable SO 3 amount [μg] × desorption rate (23)
次に、現在床温T2を用いた新規脱離SOx量の推定手法の問題点について説明する。上述した手法では、脱離可能総SOx量は、吸着後SOx分布における下限温度TLと現在床温T2の間の面積として表される。しかしながら、図10に示す吸着後SOx分布は、床温の上昇速度を一定とした条件を前提にしている。このため、床温の上昇速度が変化する条件においては、脱離可能総SOx量を精度よく推定できないおそれがある。 Next, a description will be given problems of estimation method of the current new desorption SOx amount using the bed temperature T 2. In the method described above, removable total SOx amount is expressed as the area between the current and the lower limit temperature TL in the adsorption after SOx distribution floor temperature T 2. However, the post-adsorption SOx distribution shown in FIG. 10 is based on the condition that the rate of increase in bed temperature is constant. For this reason, there is a possibility that the detachable total SOx amount cannot be accurately estimated under the condition that the rising speed of the bed temperature changes.
そこで、実施の形態1のシステムでは、現在床温T2の値だけでなく現在床温T2への変化度合も考慮した新規脱離SOx量の推定手法を採用している。以下、新規脱離SOx量の推定手法の具体的な内容について詳細に説明する。 Therefore, in the system of the first embodiment adopts a method of estimating the current new desorption SOx amount change degree also considered to a value just not current bed temperature T 2 of the floor temperature T 2. Hereinafter, specific contents of the new desorption SOx amount estimation method will be described in detail.
DOC22aに吸着しているSOxには、DOC22aのコート材22cの表層付近に吸着しているものとコート材22cの奥層付近に吸着しているものとがある。コート材22cの表層付近に吸着しているSOxは比較的低いエネルギで脱離するのに対して、コート材22cの奥層付近に吸着しているSOxは、脱離するために比較的高いエネルギが必要となる。このような知見を前提に吸着後SOx分布を考察すると、脱離時床温の低温側では主としてコート材22cの表層付近に吸着しているSOxが脱離可能となり、脱離時床温の高温側では主としてコート材22cの奥層付近に吸着しているSOxが脱離可能となると考えられる。
The SOx adsorbed on the
ここで、DOC22aの内部の温度分布は、流入するガスからの受熱によって複雑に変化する。より詳しくは、DOC22aのコート材22cの表層は、流入するガスに直接晒される位置であるため、流入するガスの温度変化に対する感度が高い。これに対して、DOC22aのコート材22cの奥層は、流入するガスに直接晒されない位置であるため、流入するガスの温度変化に対する感度が低い。このため、吸着後SOx分布を用いて推定された脱離可能SOx量は、現在床温T2が低温側の値であるほど、実際の脱離可能SOx量から乖離すると考えられる。そこで、温度変化の影響を現在床温T2に反映させる補正を施すことが考えられるが、現在床温T2の低温側を基準に一律に補正を施すと、現在床温T2の高温側において補正が過剰になるおそれがある。
Here, the temperature distribution inside the
そこで、実施の形態1のシステムでは、以下に示す制御ロジックを用いて、脱離可能総SOx量の推定に用いる現在床温T2を補正することとしている。図11は、補正後現在床温T2´を算出するための制御ロジックを示す機能ブロック図である。この図に示す補正部M30は、推定部M3の一部として構成される機能ブロックであり、ECU30により実現される。なお、以下の説明では、補正前の現在床温T2を補正前現在床温T2と称し、補正後の現在床温T2を補正後現在床温T2´と称することとする。
Therefore, in the system of the first embodiment, using the control logic described below, it is set to be currently corrected bed temperature T 2 to be used to estimate the removable total SOx amount. FIG. 11 is a functional block diagram illustrating a control logic for calculating the corrected current bed temperature T 2 ′. The correction unit M30 shown in this figure is a functional block configured as a part of the estimation unit M3, and is realized by the
この図に示すように、補正部M30は、温度変化量算出部M31と、温度補正係数算出部M32と、下限ガード部M33と、補正値算出部M34と、補正後現在床温算出部M35と、を備えている。算出部M31には、サイクル毎に補正前現在床温T2が入力される。入力される補正前現在床温T2は、DOC22aの代表温度として、温度センサ34によって検出されるDOC22aの出口温度が用いられる。なお、補正前現在床温T2は、DOC22aの内部の温度を直接検出した値を用いてもよいし、また、公知の手法を用いて推定されたDOC22aの床温を用いてもよい。算出部M31では、第t番目のサイクルの現在床温T2(t)と第(t−1)番目のサイクルの現在床温T2(t−1)の差分値、つまり現在床温T2の今回値と前回値の差分値が温度変化量ΔTとして算出される。なお、温度変化量ΔTは、補正前現在床温T2の過去値から今回値への変化度合を表す値であれば上記の差分値に限定されない。例えば、補正前現在床温T2に重畳するノイズの影響を排除するために、現在床温T2(t−1)と現在床温T2(t−2)との平均値を算出し、現在床温T2(t)と当該平均値との差分値を温度変化量ΔTとしてもよい。また、他の公知のなまし処理を温度変化量ΔTの算出に適用してもよい。
As shown in this figure, the correction unit M30 includes a temperature change amount calculation unit M31, a temperature correction coefficient calculation unit M32, a lower limit guard unit M33, a correction value calculation unit M34, and a corrected current bed temperature calculation unit M35. It is equipped with. The calculation unit M31 is before correction for each cycle currently bed temperature T 2 is input. Uncorrected current floor temperature T 2 which is input, as a representative temperature of DOC22a, outlet temperature of DOC22a detected by the
算出された温度変化量ΔTは、ガード部M33に入力される。ガード部M33では、固定値ゼロと温度変化量ΔTとの最大値選択によって、温度変化量ΔTが下限値ゼロでガードされる。 The calculated temperature change amount ΔT is input to the guard unit M33. In the guard unit M33, the temperature change amount ΔT is guarded at the lower limit value zero by selecting the maximum value between the fixed value zero and the temperature change amount ΔT.
算出部M32では、補正前現在床温T2の入力を受けて温度補正係数Kを算出する。図12は、補正前現在床温T2に温度補正係数Kを関連付けたマップの一例を示している。温度補正係数Kは、例えば図12に示すマップから算出することができる。このマップでは、補正前現在床温T2が400℃より小さい領域では、温度補正係数Kが一定値(正値)に設定され、500℃より大きい領域では、温度補正係数Kがゼロに設定されている。また、このマップでは、補正前現在床温T2が400℃から500℃の間において、補正前現在床温T2が大きいほど温度補正係数Kが小さくなるように設定されている。これにより、温度補正係数Kは、補正前現在床温T2が大きくなるに従い単調非増加で小さくなるように変化する。なお、温度補正係数Kの算出は、図12に示すマップを用いる方法に限られない。すなわち、補正前現在床温T2が大きくなるに従い単調非増加で小さくなる傾向が規定された他のマップを用いる構成でもよい。 The calculation unit M32, receives the uncorrected current inputs bed temperature T 2 calculates the temperature correction coefficient K. Figure 12 shows an example of a map that associates a temperature correction coefficient K before correction current to the bed temperature T 2. The temperature correction coefficient K can be calculated from, for example, the map shown in FIG. In this map, the pre-correction at the current floor temperature T 2 is 400 ° C. smaller area is the temperature correction coefficient K is set to a constant value (positive value), the 500 ° C. larger area, the temperature correction coefficient K is set to zero ing. Further, in this map, the uncorrected current floor temperature T 2 is in between 500 ° C. from 400 ° C., as the temperature correction coefficient K current is large bed temperature T 2 before correction is set to be smaller. Accordingly, the temperature correction coefficient K is changed to be smaller monotonic nonincreasing accordance before correction current floor temperature T 2 becomes large. The calculation of the temperature correction coefficient K is not limited to the method using the map shown in FIG. Namely, before correction may be configured to use other maps smaller tendency defined by monotone nonincreasing accordance currently bed temperature T 2 becomes large.
算出部M34には、温度補正係数Kと温度変化量ΔTが入力される。算出部M34では、温度変化量ΔTに温度補正係数Kを乗算することにより、補正値(K×ΔT)が算出される。補正値(K×ΔT)は、補正前現在床温T2が大きいほど、つまり温度補正係数Kが小さいほど小さな値となり、温度補正係数Kがゼロになると補正値(K×ΔT)もゼロとなる。 A temperature correction coefficient K and a temperature change amount ΔT are input to the calculation unit M34. The calculation unit M34 calculates a correction value (K × ΔT) by multiplying the temperature change amount ΔT by the temperature correction coefficient K. Correction value (K × [Delta] T), the more correction larger before the current bed temperature T 2, i.e. becomes smaller the value temperature correction coefficient K is small, a correction value temperature correction coefficient K is zero (K × [Delta] T) is also zero Become.
算出部M35には、補正前現在床温T2と補正値(K×ΔT)が入力される。算出部M35では、補正前現在床温T2に補正値(K×ΔT)を加算することによって補正後現在床温T2´が算出される。これにより、補正前現在床温T2が小さい場合には温度変化量ΔTの反映度合の大きい補正後現在床温T2´が算出され、補正前現在床温T2が大きい場合には温度変化量ΔTの反映度合が小さい補正後現在床温T2´が算出される。このように、補正部M30の処理によれば、補正前現在床温T2の大きさに応じて温度変化量ΔTの反映度合を変化させることができる。 The calculation unit M35 are uncorrected current bed temperature T 2 and the correction value (K × ΔT) are inputted. The calculation unit M35 calculates the corrected current bed temperature T 2 ′ by adding the correction value (K × ΔT) to the current bed temperature T 2 before correction. As a result, when the current bed temperature T 2 before correction is small, the corrected current bed temperature T 2 ′ having a large degree of reflection of the temperature change ΔT is calculated, and when the current bed temperature T 2 before correction is large, the temperature change. A corrected current bed temperature T 2 ′ with a small reflection degree of the amount ΔT is calculated. Thus, according to the process of the correction unit M30, it is possible to change the reflection degree of the temperature change ΔT in accordance with the magnitude of the correction before the current floor temperature T 2.
なお、補正後現在床温T2´を算出する制御ロジックは、上記の機能ブロックの構成に限られない。すなわち、補正前現在床温T2に温度変化量ΔTを反映させる補正を行うことにより補正後現在床温T2´を算出する制御ロジックにおいて、温度変化量ΔTの補正前現在床温T2への反映度合が補正前現在床温T2が大きくなるに従って単調非増加で小さくなるように構成されていればよい。 The control logic for calculating the corrected current bed temperature T 2 ′ is not limited to the configuration of the functional block described above. Namely, before correction in the control logic for calculating a corrected current floor temperature T 2 'by performing the current compensation to reflect the temperature variation ΔT in bed temperature T 2, the temperature variation ΔT uncorrected Currently bed temperature T 2 the reflection degree may be composed to be smaller monotonic nonincreasing according to the current floor temperature T 2 increases before correction.
実施の形態1のシステムでは、上式(20)及び(21)において補正後の現在床温T2が用いられる。これにより、現在床温の低い温度域では温度変化量ΔTの現在床温への反映度合が高められるので、推定された脱離可能SOx量を実際の脱離可能SOx量に近づけることができる。また、現在床温の低い温度域では温度変化量ΔTの現在床温への反映度合が弱められるので、現在床温が過剰に補正されることによって推定された脱離可能SOx量が実際の脱離可能SOx量から乖離することを防ぐことができる。 In the system of the first embodiment, the current floor temperature T 2 of the corrected in the above equation (20) and (21) are used. Thereby, in the temperature range where the current bed temperature is low, the degree of reflection of the temperature change amount ΔT on the current bed temperature is increased, and thus the estimated detachable SOx amount can be brought close to the actual detachable SOx amount. In addition, in the temperature range where the current bed temperature is low, the degree of reflection of the temperature change amount ΔT on the current bed temperature is weakened, so that the amount of detachable SOx estimated by excessively correcting the current bed temperature is the actual desorption. Deviation from the detachable SOx amount can be prevented.
上式(20)及び(21)において算出された脱離可能SOx量は、上式(22)及び(23)による新規脱離SO2量の算出に使用される。これにより、現在床温の変化度合によらず新規脱離SO2量を精度よく推定することが可能となる。 Eliminable SOx amount calculated in the above equation (20) and (21) are used for the calculation of the above equation (22) and new desorbed SO 2 amount of (23). This makes it possible to accurately estimate the amount of new desorbed SO 2 regardless of the degree of change in the current bed temperature.
図3に戻り、推定部M6は、推定部M5で推定した新規脱離SOx量を吸着後SOx分布に反映させて、最終吸着SOx分布を推定する。 Returning to FIG. 3, the estimation unit M6 reflects the newly desorbed SOx amount estimated by the estimation unit M5 in the post-adsorption SOx distribution, and estimates the final adsorption SOx distribution.
推定部M6は具体的に、推定部M5で推定した新規脱離SOx量の分だけSOxが脱離し、吸着後SOx分布の形状が変形すると仮定して、最終吸着SOx分布(脱離後SOx分布)を推定する。最終吸着SOx分布と吸着後SOx分布の関係について、SO2を例とした図13を参照しながら説明する。図13は、最終吸着SOx分布と吸着後SOx分布の関係を説明するための図である。なお、この図の横軸のTLおよびTHは、上述した下限温度および上限温度に、それぞれ相当している。この図に示すように、吸着後SO2分布の下限温度TLからの積分値が新規脱離SO2量に一致するときの吸着後SO2分布の面積、即ち、下限温度TLから温度TdSO2までの面積を、吸着後SOx分布から削った後に残る分布が、最終吸着SO2分布となる。 Specifically, the estimation unit M6 assumes that the SOx is desorbed by the amount of the new desorption SOx estimated by the estimation unit M5 and the shape of the post-adsorption SOx distribution is deformed, and the final adsorption SOx distribution (post-desorption SOx distribution). ). The relationship between the final adsorption SOx distribution and the post-adsorption SOx distribution will be described with reference to FIG. 13 using SO 2 as an example. FIG. 13 is a diagram for explaining the relationship between the final adsorption SOx distribution and the post-adsorption SOx distribution. Note that TL and TH on the horizontal axis in this figure correspond to the above-described lower limit temperature and upper limit temperature, respectively. As shown in this figure, the area of the post-adsorption SO 2 distribution when the integrated value from the lower limit temperature TL of the post-adsorption SO 2 distribution matches the newly desorbed SO 2 amount, that is, from the lower limit temperature TL to the temperature Td SO 2. The distribution remaining after removing the area of the SOx distribution after adsorption is the final adsorption SO 2 distribution.
図13の温度TdSO2が床温T1を上回る場合には、DOC22aからSO2が全て脱離していることを意味する。これを考慮すると、第t番目のサイクルにおける最終吸着SO2分布は床温T1を変数とする次式(24)で表され、第t番目のサイクルにおける最終吸着SO3分布は次式(25)で表されることになる。なお、式(25)の温度TdSO3は、吸着後SO3分布の下限温度TLからの積分値が新規脱離SO3量に一致するときの床温T1に相当している。 When the above temperature Td SO2 floor temperature T 1 of the FIG. 13 means that the SO 2 is all desorbed from DOC22a. Considering this, the final adsorption SO 2 distribution in the t-th cycle is expressed by the following equation (24) using the bed temperature T 1 as a variable, and the final adsorption SO 3 distribution in the t-th cycle is expressed by the following equation (25 ). The temperature Td SO3 of formula (25) corresponds to the bed temperature T 1 of the when the integrated value from the lower limit temperature TL of SO 3 distribution after adsorption to match the new desorption SO 3 amount.
なお、新規脱離SO2量と温度TdSO2の関係は次式(26)で表すことができ、新規脱離SO3量と温度TdSO3の関係は次式(27)で表すことができる。
図3に戻り、推定部M7は、上述したすり抜けSOxのうちDOC22aからSO3の状態で排出されるSOxの量(以下「すり抜けSO3量」ともいう。)を推定する。
Returning to FIG. 3, the estimation unit M7 estimates the amount of SOx discharged from the
図2で説明したように、DOC22aでは貴金属22dに吸着したSO2の一部がSO3に転化する。この転化がすり抜けSOx中のSO2にも起こると仮定して、推定部M7では、すり抜け量と現在床温T2を変数とする次式(28)により、第t番目のサイクルにおけるすり抜けSO3量を推定する。なお、すり抜けSOxのうちDOC22aからSO2の状態で排出されるSOxの量(以下「すり抜けSO2量」ともいう。)は、次式(29)で表すことができる。
すり抜けSO3量(すり抜け量(t),T2(t))[μg/s]=すり抜けSOx量(t)×SO3転化率map(T2(t)) ・・・(28)
すり抜けSO2量(すり抜け量(t),T2(t))[μg/s]=すり抜けSOx量(t)×{1−SO3転化率map(T2(t))} ・・・(29)
As described with reference to FIG. 2, in the
Amount of slip-through SO 3 (amount of slip-through (t), T 2 (t)) [μg / s] = amount of slip-through SOx (t) × SO 3 conversion map (T 2 (t)) (28)
Amount of slip-through SO 2 (amount of slip-through (t), T 2 (t)) [μg / s] = amount of slip-through SOx (t) × {1-SO 3 conversion rate map (T 2 (t))} 29)
式(28)および(29)のSO3転化率map(T2(t))は、第t番目のサイクルにおいて「DOC22aに流入するSOx」のうち、DOC22aからSO3の状態で排出されるSOxの割合(即ち、SO3転化率)が、DOC22aの現在床温T2によって変わるという特性に基づいて作成されたマップである。図14は、SO3転化率mapを説明するための図である。SO3転化率mapの特性は、例えば図14に示すように、現在床温T2がある温度域Bにある場合はSO3転化率が高くなり、この温度域Bよりも低温側や高温側では、SO2からSO3への転化が起こり難くなるというものである。このようなマップは、例えばECU30のROMに記憶させておくことができ、現在床温T2に応じて適宜読み出すことができる。
The SO 3 conversion rate map (T 2 (t)) of the equations (28) and (29) is the SOx discharged in the state of SO 3 from the
図3に戻り、算出部M8は、DOC22aの床温上昇中にDOC22aから脱離してもよいSO3の量(以下「許容脱離SO3量」ともいう。)を算出する。許容脱離SO3量について、図15を参照して説明する。図15は、許容脱離SO3量を説明するための図である。この図に示す制約SO3量は、サルフェート白煙に関する制約に相当しており、この図においてはすり抜けSO3量と許容脱離SO3量の和が制約SO3量と等しくなっている。すり抜けSO3量と許容脱離SO3量の和は、DOC22aの下流におけるSO3の量であることから、この和の値が制約SO3量よりも小さい値であれば制約が満たされることになる。
Returning to FIG. 3, the calculation unit M8 calculates the amount of SO 3 that may be desorbed from the
制約SO3量は、第t番目のサイクルにおけるエンジン10の排気流量(ガス流量)を変数とする次式(30)により算出することができる。なお、エンジン10の排気流量は、例えばエアフローメータ32の検出値を用いることができる。
制約SO3量(ガス流量(t))[μg/s]=制約SO3濃度[ppm]×ガス流量(t)[g/s]÷空気の平均モル質量×SO3分子量 ・・・(30)
The amount of constraint SO 3 can be calculated by the following equation (30) using the exhaust flow rate (gas flow rate) of the
Restricted SO 3 amount (gas flow rate (t)) [μg / s] = restricted SO 3 concentration [ppm] × gas flow rate (t) [g / s] ÷ average molar mass of air × SO 3 molecular weight (30 )
従って、この制約SO3量とすり抜けSO3量を変数とする次式(31)を許容脱離SO3量が満たせば、制約が満たされることになる。
許容脱離SO3量(制約SO3量(ガス流量(t)),すり抜けSO3量(t))[μg/s]≦制約SO3量(ガス流量(t))[μg/s]−すり抜けSO3量(t)[μg/s] ・・・(31)
なお、以下の説明においては、許容脱離SO3量(制約SO3量(ガス流量(t)),すり抜けSO3量(t))を単に許容脱離SO3量(t)ともいう。
Therefore, the constraint is satisfied if the allowable desorption SO 3 amount satisfies the following equation (31) using the constraint SO 3 amount and the slip-through SO 3 amount as variables.
Allowable desorption SO 3 amount (restricted SO 3 amount (gas flow rate (t)), slip-through SO 3 amount (t)) [μg / s] ≦ restricted SO 3 amount (gas flow rate (t)) [μg / s] − Slip-through SO 3 amount (t) [μg / s] (31)
In the following description, the allowable desorption SO 3 amount (restricted SO 3 amount (gas flow rate (t)), slip-through SO 3 amount (t)) is also simply referred to as allowable desorption SO 3 amount (t).
図3に戻り、算出部M9は、PM再生制御中のサルフェート白煙の発生を抑制するための第t番目のサイクルにおける目標床温Ttrgを算出する。目標床温Ttrgについて、図16を参照しながら説明する。図16は、目標床温Ttrgを説明するための図である。なお、この図の横軸のTLおよびTHは、上述した下限温度および上限温度に、それぞれ相当している。この図に示すように、最終吸着SO3分布の低温側からの積分値に脱離率を乗じた値が、算出部M8で算出した許容脱離SO3量に一致するときの床温T1が目標床温Ttrgに相当する。 Returning to FIG. 3, the calculation unit M9 calculates the target bed temperature Ttrg in the t-th cycle for suppressing the generation of sulfate white smoke during PM regeneration control. The target bed temperature Ttrg will be described with reference to FIG. FIG. 16 is a diagram for explaining the target bed temperature Ttrg. Note that TL and TH on the horizontal axis in this figure correspond to the above-described lower limit temperature and upper limit temperature, respectively. As shown in this figure, the bed temperature T 1 when the value obtained by multiplying the integral value from the low temperature side of the final adsorption SO 3 distribution by the desorption rate matches the allowable desorption SO 3 amount calculated by the calculation unit M8. Corresponds to the target bed temperature Ttrg.
なお、第t番目のサイクルにおける許容脱離SO3量と目標床温Ttrgの関係は次式(32)で表すことができる。式(32)の脱離率には設定値が使用され、例えばECU30のROMに記憶させておくことができる。
The relationship between the allowable desorption SO 3 amount and the target bed temperature Ttrg in the t-th cycle can be expressed by the following equation (32). A set value is used as the desorption rate in the equation (32), and can be stored in the ROM of the
図17を参照して、本実施の形態による効果を説明する。図17は、現在床温の温度補正の効果を説明するためのタイムチャートである。この図の1段目のチャートは、補正前の現在床温T2と補正後の現在床温T2´の時間変化をそれぞれ示している。また、2段目のチャートは、補正前の現在床温T2を用いた場合の新規脱離SO3量の推定値と実測値をそれぞれ示している。これらのチャートに示すように、現在床温T2の高温側の領域では現在床温T2に対する補正度合が小さくなっている。これにより、現在床温T2の高温側の領域において過剰な補正を抑えて新規脱離SO3量の推定値を実測値に近づけることができる。 With reference to FIG. 17, the effect by this Embodiment is demonstrated. FIG. 17 is a time chart for explaining the effect of temperature correction of the current bed temperature. The first chart in this figure shows the time variation of the current bed temperature T 2 before correction and the current bed temperature T 2 ′ after correction. The second chart shows an estimated value and a measured value of the new desorbed SO 3 amount when the current bed temperature T 2 before correction is used. As shown in these charts, in the region of the current high-temperature side of the bed temperature T 2 is a correction degree current for bed temperature T 2 is smaller. This makes it possible to close the currently measured value of the estimated value of the new desorbed SO 3 content by suppressing excessive correction in the high temperature side region of the floor temperature T 2.
なお、上述した実施の形態においては、推定部M1が本発明の「流入SOx量推定手段」に相当し、推定部M2が本発明の「SOx飽和率推定手段」に相当し、推定部M3が本発明の「新規吸着SOx量推定手段」および「すり抜けSOx量推定手段」に相当し、推定部M4が本発明の「吸着後SOx分布推定手段」に相当し、推定部M5が本発明の「新規脱離SOx量推定手段」に相当し、推定部M6が本発明の「最終吸着SOx分布推定手段」に相当し、推定部M7が本発明の「すり抜けSO3量推定手段」に相当し、算出部M8が本発明の「許容脱離SO3量算出手段」に相当し、算出部M9が本発明の「目標温度算出手段」に相当している。 In the embodiment described above, the estimation unit M1 corresponds to the “inflow SOx amount estimation unit” of the present invention, the estimation unit M2 corresponds to the “SOx saturation rate estimation unit” of the present invention, and the estimation unit M3 The estimation unit M4 corresponds to the “new adsorption SOx amount estimation unit” and the “passing SOx amount estimation unit” of the present invention, the estimation unit M4 corresponds to the “post-adsorption SOx distribution estimation unit” of the present invention, and the estimation unit M5 of the present invention “ The estimation unit M6 corresponds to the “final adsorption SOx distribution estimation unit” of the present invention, and the estimation unit M7 corresponds to the “passing SO 3 amount estimation unit” of the present invention; The calculation unit M8 corresponds to “allowable desorption SO 3 amount calculation unit” of the present invention, and the calculation unit M9 corresponds to “target temperature calculation unit” of the present invention.
また、上述した実施の形態においては、現在床温T2及び補正前現在床温T2が本発明の「代表温度」に相当し、補正後現在床温T2が本発明の「補正後代表温度」に相当し、温度補正係数Kが本発明の「補正係数」に相当し、温度変化量ΔTが本発明の「補正値」に相当し、算出部M31が本発明の「算出手段」に相当し、算出部M35が本発明の「補正手段」に相当している。 Further, in the embodiment described above, the current corresponds to the "representative temperature" bed temperature T 2 and the uncorrected current floor temperature T 2 the present invention, the corrected current floor temperature T 2 "corrected representative of the present invention The temperature correction coefficient K corresponds to the “correction coefficient” of the present invention, the temperature change amount ΔT corresponds to the “correction value” of the present invention, and the calculation unit M31 serves as the “calculation means” of the present invention. The calculation unit M35 corresponds to the “correction unit” of the present invention.
ところで、上述した実施の形態では、燃料添加弁24からの燃料の添加によりPM再生制御を行った。しかしこのPM再生制御を、インジェクタ12からの燃料の噴射(具体的には、メイン噴射よりも後のサブ噴射(例えばポスト噴射))により行ってもよい。この場合は、式(1)の排気添加量を、インジェクタ12からのサブ噴射量に置き換えればよい。
また、上述した実施の形態では、PM再生制御中を例としてDOC22aの床温の目標温度を算出した。しかし、DOC22aからSOxを脱離させる制御をPM再生制御と併せて行うような場合に、この脱離制御中に上述した手法によりDOC22aの床温の目標温度を算出してもよい。このように、上述した目標温度の算出手法は、DOC22aからSOxが脱離する温度域までDOC22aの床温を上昇させる制御一般に適用することができる。
また、上述した実施の形態では、DOC22aとDPF22bを備える排気浄化装置22を例として説明した。しかし、DOC22aにおけるHC等の酸化機能をDPF22bに付与して、排気浄化装置22からDOC22aを省略してもよい。この場合は、上述した目標温度の算出手法を、酸化機能が付与されたDPF22bに適用することで上述した実施の形態と同様の効果を得ることができる。
また、上述した実施の形態ではエンジン10がターボチャージャ16を備えるとしたが、エンジン10がターボチャージャ16を備えていなくてもよい。即ち、上述した目標温度の算出手法は、非過給ディーゼルエンジンのシステムにも適用できる。
In the above-described embodiment, PM regeneration control is performed by adding fuel from the
In the above-described embodiment, the target temperature of the bed temperature of the
Moreover, in embodiment mentioned above, the exhaust
In the above-described embodiment, the
10 ディーゼルエンジン
12 インジェクタ
20 排気管
22 排気浄化装置
22a DOC
22b DPF
22c コート材
22d 貴金属
24 燃料添加弁
30 ECU
DESCRIPTION OF
22b DPF
Claims (7)
所定の制御周期ごとに前記浄化装置の温度の代表値である代表温度を取得する温度取得手段と、
前記浄化装置に流入するSOx量を流入SOx量として前記制御周期ごとに推定する流入SOx量推定手段と、
前記流入SOx量と前記代表温度とを用いて、前記浄化装置から新たに脱離するSOx量を新規脱離SOx量として前記制御周期ごとに推定する新規脱離SOx量推定手段と、
前記新規脱離SOx量を用いて、前記浄化装置の温度上昇中の各温度において前記浄化装置に最終的に吸着するSOx量を前記代表温度に関連付けたグラフとして表される最終吸着SOx分布を前記制御周期ごとに推定する最終吸着SOx分布推定手段と、
前記最終吸着SOx分布と許容脱離SO3量とを用いて、前記目標温度を前記制御周期ごとに算出する目標温度算出手段と、を備え、
前記新規脱離SOx量推定手段は、
前記代表温度の過去値から今回値への変化度合を表す補正値を算出する算出手段と、
前記代表温度に前記補正値を反映させる補正を行うことにより補正後代表温度を算出する補正手段と、を含み、前記流入SOx量と前記補正後代表温度を用いて前記新規脱離SOx量を推定するように構成され、
前記補正手段は、前記変化度合の前記補正後代表温度への反映度合が、前記代表温度が大きくなるに従って単調非増加で小さくなるように構成されていることを特徴とするエンジン制御装置。 An engine control device that executes a temperature rise control for increasing a temperature of a purification device provided in an exhaust pipe of a diesel engine to a target temperature in a temperature range where SOx is desorbed from the purification device,
Temperature acquisition means for acquiring a representative temperature, which is a representative value of the temperature of the purification device, for each predetermined control cycle;
An inflow SOx amount estimating means for estimating an SOx amount flowing into the purifier as an inflow SOx amount for each control period;
Using the inflow SOx amount and the representative temperature, a new desorption SOx amount estimation means for estimating a SOx amount newly desorbed from the purification apparatus as a new desorption SOx amount for each control period;
Using the new desorbed SOx amount, the final adsorption SOx distribution represented as a graph in which the SOx amount finally adsorbed to the purification device at each temperature during the temperature increase of the purification device is related to the representative temperature is expressed as the graph. A final adsorption SOx distribution estimating means for estimating each control cycle;
Using the final adsorption SOx distribution and the allowable desorption SO 3 amount, a target temperature calculation means for calculating the target temperature for each control cycle,
The new desorbed SOx amount estimating means includes:
Calculation means for calculating a correction value representing the degree of change from the past value of the representative temperature to the current value;
Correction means for calculating a corrected representative temperature by performing correction that reflects the correction value on the representative temperature, and estimating the new desorbed SOx amount using the inflow SOx amount and the corrected representative temperature. Configured to
The engine control apparatus according to claim 1, wherein the correction means is configured such that the degree of reflection of the degree of change in the corrected representative temperature decreases in a monotonically non-increasing manner as the representative temperature increases.
前記代表温度が大きくなるに従って単調非増加で小さくなる補正係数を算出し、
前記補正値に前記補正係数を乗じた値を前記代表温度に加算することにより前記補正後代表温度を算出する
ように構成されていることを特徴とする請求項1に記載のエンジン制御装置。 The correction means includes
Calculating a correction coefficient that decreases monotonically without increasing as the representative temperature increases,
The engine control device according to claim 1, wherein the corrected representative temperature is calculated by adding a value obtained by multiplying the correction value by the correction coefficient to the representative temperature.
前記流入SOx量と前記SOx飽和率とを用いて、前記浄化装置に流入して前記浄化装置に新たに吸着するSOx量を新規吸着SOx量として前記制御周期ごとに推定する新規吸着SOx量推定手段と、
前記新規吸着SOx量を用いて、前記浄化装置に新たなSOxが吸着した後の前記吸着SOx分布を吸着後SOx分布として前記制御周期ごとに推定する吸着後SOx分布推定手段と、を備え、
前記新規脱離SOx量推定手段は、前記吸着後SOx分布と前記代表温度とを用いて、前記新規脱離SOx量を前記制御周期ごとに推定するように構成されていることを特徴とする請求項1乃至3の何れか1項に記載のエンジン制御装置。 The adsorption SOx distribution represented as a graph in which the amount of SOx adsorbed to the purification device at each temperature during the temperature increase of the purification device is related to the representative temperature, and the purification device at each temperature during the temperature increase of the purification device. SOx saturation rate estimating means for estimating the SOx saturation rate in the purification apparatus for each control cycle using a saturated SOx distribution represented as a graph in which the maximum amount of SOx adsorbed on the representative temperature is related to the representative temperature;
Using the inflow SOx amount and the SOx saturation rate, a new adsorption SOx amount estimation means that estimates the SOx amount that flows into the purification device and is newly adsorbed to the purification device as a new adsorption SOx amount for each control cycle. When,
A post-adsorption SOx distribution estimation unit that estimates the adsorption SOx distribution after the new SOx is adsorbed on the purification device as the post-adsorption SOx distribution for each control cycle using the new adsorption SOx amount;
The new desorption SOx amount estimation means is configured to estimate the new desorption SOx amount for each control cycle using the post-adsorption SOx distribution and the representative temperature. Item 4. The engine control device according to any one of Items 1 to 3.
前記浄化装置においてSO3に転化するSO2の転化率と前記代表温度との関係を表した転化率マップと、今回の推定サイクルにおける前記代表温度と、前記すり抜けSOx量と、を用いて、前記浄化装置にSOxの状態で流入して前記浄化装置に吸着することなくすり抜けてSO3の状態で排出されるSO3量をすり抜けSO3量として前記制御周期ごとに推定するすり抜けSO3量推定手段と、
サルフェート白煙に関する制約に相当する前記浄化装置の下流におけるSO3量と前記すり抜けSO3量とを用いて、前記浄化装置から脱離することが許されるSO3量を前記許容脱離SO3量として前記制御周期ごとに算出する許容脱離SO3量算出手段と、
をさらに備えることを特徴とする請求項4に記載のエンジン制御装置。 By using the new amount of adsorbed SOx, the amount of SOx passing through the purification device and slipping through without being adsorbed by the purification device is estimated as a slip-through SOx amount for each control cycle;
Using the conversion rate map representing the relationship between the conversion rate of SO 2 converted to SO 3 and the representative temperature in the purification device, the representative temperature in the current estimation cycle, and the slip-through SOx amount, slipped SO 3 amount estimating means for estimating for each of the control period as the SO 3 content slipped through SO 3 amount discharged in the form of SO 3 to pass through without the purifying device and flows in the form of SOx adsorbed to the purifier When,
Using the SO 3 content slipping said the SO 3 content in the downstream of the purification device corresponding to constraints on sulfate white smoke, the purifying device wherein the SO 3 amount that is allowed to detach from the allowable desorption SO 3 weight An allowable desorption SO 3 amount calculating means for calculating each control cycle as
The engine control device according to claim 4, further comprising:
前記目標温度まで上昇させる制御を、前記フィルタに捕集された微粒子量の推定値が除去要求量に到達したときに開始することを特徴とする請求項1乃至6の何れか1項に記載のエンジン制御装置。 The purification device includes a filter that collects particulates flowing through the exhaust pipe,
The control for increasing the temperature to the target temperature is started when an estimated value of the amount of fine particles collected by the filter reaches a removal request amount. Engine control device.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017020560A JP6547779B2 (en) | 2017-02-07 | 2017-02-07 | Engine control unit |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017020560A JP6547779B2 (en) | 2017-02-07 | 2017-02-07 | Engine control unit |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2018127921A true JP2018127921A (en) | 2018-08-16 |
JP6547779B2 JP6547779B2 (en) | 2019-07-24 |
Family
ID=63173837
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017020560A Expired - Fee Related JP6547779B2 (en) | 2017-02-07 | 2017-02-07 | Engine control unit |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6547779B2 (en) |
-
2017
- 2017-02-07 JP JP2017020560A patent/JP6547779B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP6547779B2 (en) | 2019-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7707826B2 (en) | System for controlling triggering of adsorber regeneration | |
JP6582409B2 (en) | Exhaust purification system | |
JP2004076589A (en) | Filter control method and device | |
US20130298529A1 (en) | System amd method for controlling an after-treatment component of a compression-ignition engine | |
JP4574460B2 (en) | Exhaust gas purification device for internal combustion engine | |
CN107407182B (en) | Exhaust gas purification device and control method thereof | |
JP2016133064A (en) | Exhaust emission control system | |
JP2016188604A (en) | Exhaust emission control device | |
JP6319282B2 (en) | Engine control device | |
JP2004132358A (en) | Filter control device | |
JP2004293413A (en) | Exhaust emission control system | |
JP6547779B2 (en) | Engine control unit | |
JP5912494B2 (en) | Diesel engine exhaust purification system | |
JP6418014B2 (en) | Exhaust purification system | |
WO2016117568A1 (en) | EXHAUST GAS PURIFICATION SYSTEM, AND NOx PURIFICATION CAPACITY RESTORATION METHOD | |
JP6806025B2 (en) | Engine control unit | |
JP6183659B2 (en) | Exhaust gas purification device | |
CN110945218B (en) | Exhaust gas purification system | |
JP7122873B2 (en) | Exhaust gas temperature control method and exhaust gas purification device | |
JP2015017516A (en) | Exhaust emission control device for internal combustion engine | |
JP2016118135A (en) | Exhaust emission control system | |
JP6424618B2 (en) | Exhaust purification system | |
JP2016180383A (en) | Catalyst temperature estimation device | |
JP2016125375A (en) | Exhaust emission control system | |
JP2016133062A (en) | Exhaust emission control system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20180827 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20190528 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20190529 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190610 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 6547779 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |
|
LAPS | Cancellation because of no payment of annual fees |