(57)【特許請求の範囲】
【請求項１】
燃焼システムの燃焼に用いる燃料に含まれている複数種類の分子構造の各々の混合割合を取得する取得部（Ｓ１１）と、
前記取得部により取得された前記混合割合に基づき、燃焼に伴う煤成分の生じやすさを表わした煤生成指数を算出する煤算出部（Ｓ１２）と、
前記取得部により取得された前記混合割合に基づき、燃焼に伴い生じた可溶性有機成分の付着しやすさを表わした付着指数を算出する付着指数算出部（Ｓ１３）と、
前記煤算出部により算出された前記煤生成指数、および前記付着指数算出部により算出された前記付着指数に基づき、前記燃焼システムの所定部位に付着する可溶性有機成分の堆積量を推定する堆積量推定部（Ｓ１４）と、
を備え、
前記燃料に含まれる成分のうち、燃焼前に分解して重合することでアロマ類成分を形成する成分をアロマ類可変成分と呼ぶ場合において、
前記煤算出部は、
前記取得部により取得された複数種類の前記混合割合のうちアロマ類成分の混合割合が多いほど、前記煤成分が生じやすいことを表わす値に前記煤生成指数を算出し、
なおかつ、前記取得部により取得された複数種類の前記混合割合のうち前記アロマ類可変成分の混合割合が多いほど、前記煤成分が生じやすいことを表わす値に前記煤生成指数を算出し、
前記煤算出部の算出に用いられる前記アロマ類可変成分には、ナフテン、側鎖パラフィンおよび直鎖パラフィンが含まれているデポジット推定装置。
【請求項２】
前記付着指数算出部は、
前記取得部により取得された複数種類の前記混合割合が、燃料の揮発性が低くなるような値の組み合わせであるほど、前記可溶性有機成分が付着しやすいことを表わす値に前記付着指数を算出する請求項１に記載のデポジット推定装置。
【請求項３】
前記付着指数算出部は、
前記取得部により取得された複数種類の前記混合割合が、燃料の平均炭素数が多くなるような値の組み合わせであるほど、前記可溶性有機成分が付着しやすいことを表わす値に前記付着指数を算出する請求項１または２に記載のデポジット推定装置。
【請求項４】
前記付着指数算出部は、
前記取得部により取得された複数種類の前記混合割合が、燃料の動粘度が高くなるような値の組み合わせであるほど、前記可溶性有機成分が付着しやすいことを表わす値に前記付着指数を算出する請求項１〜３のいずれか１つに記載のデポジット推定装置。
【請求項５】
燃焼システムの燃焼に用いる燃料に含まれている複数種類の分子構造の各々の混合割合を取得する取得部（Ｓ１１）と、
前記取得部により取得された前記混合割合に基づき、燃焼に伴う煤成分の生じやすさを表わした煤生成指数を算出する煤算出部（Ｓ１２）と、
前記取得部により取得された前記混合割合に基づき、燃焼に伴い生じた可溶性有機成分の付着しやすさを表わした付着指数を算出する付着指数算出部（Ｓ１３）と、
前記煤算出部により算出された前記煤生成指数、および前記付着指数算出部により算出された前記付着指数に基づき、前記燃焼システムの所定部位に付着する可溶性有機成分の堆積量を推定する堆積量推定部（Ｓ１４）と、
前記堆積量推定部により推定された前記堆積量に応じて、前記堆積量を低減させるように前記燃焼システムの作動を制御する制御部（Ｓ１６）と、
を備え、
前記燃料に含まれる成分のうち、燃焼前に分解して重合することでアロマ類成分を形成する成分をアロマ類可変成分と呼ぶ場合において、
前記煤算出部は、
前記取得部により取得された複数種類の前記混合割合のうちアロマ類成分の混合割合が多いほど、前記煤成分が生じやすいことを表わす値に前記煤生成指数を算出し、
なおかつ、前記取得部により取得された複数種類の前記混合割合のうち前記アロマ類可変成分の混合割合が多いほど、前記煤成分が生じやすいことを表わす値に前記煤生成指数を算出し、
前記煤算出部の算出に用いられる前記アロマ類可変成分には、ナフテン、側鎖パラフィンおよび直鎖パラフィンが含まれている燃焼システム制御装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、燃焼システムの所定部位に付着する可溶性有機成分の堆積量を推定するデポジット推定装置に関する。
【背景技術】
【０００２】
燃焼システムでの燃焼に伴い生成される可溶性有機成分（つまりＳＯＦ成分）は粘着性が高い。そのため、燃焼システムのうち排ガスに晒される部分にＳＯＦ成分が付着して堆積していき、燃焼システムの作動不良を招くことが懸念される。このような作動不良を未然に防ぐためには、ＳＯＦ成分の堆積量（つまりデポジット量）が所定量に達したタイミングで、デポジットを低減させることを要する。例えば、内燃機関の停止後に、ＳＯＦ成分が付着しているバルブを開閉作動させてＳＯＦ成分を振り落とす制御を実施したり、デポジットを燃焼除去させたり、排ガス中のＳＯＦ成分量を少なくさせるように燃焼状態を制御することを要する。
【０００３】
このような制御を必要最小限の頻度で実施するためには、デポジット量を精度良く推定することが重要になる。例えば特許文献１には、燃料噴射弁の噴孔周りに堆積していくＳＯＦ成分の量（つまりデポジット量）を、燃料噴射弁からの燃料噴射量、噴孔の雰囲気温度、圧力、排ガス中のＮＯｘ濃度等に基づき推定する技術が開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１０−１１１２９３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、どのような燃料が用いられるかよって、ＳＯＦ成分の発生量や粘性が異なってくる。例えば、高粘性のＳＯＦ成分が発生するような燃料が用いられた場合、ＳＯＦ成分が付着しやすくなるので、デポジット量が多くなる。そして、特許文献１に記載のデポジット量の推定手法では、どのような燃料が用いられているかについては考慮されていないため、推定精度が悪い。
【０００６】
本発明は、上記問題を鑑みてなされたもので、その目的は、デポジット量を高精度で推定可能なデポジット推定装置および燃焼システム制御装置を提供することにある。
【課題を解決するための手段】
【０００７】
ここに開示される発明は上記目的を達成するために以下の技術的手段を採用する。なお、特許請求の範囲およびこの項に記載した括弧内の符号は、後述する実施形態に記載の具体的手段との対応関係を示すものであって、発明の技術的範囲を限定するものではない。
【０００８】
開示される第１の発明は、
燃焼システムの燃焼に用いる燃料に含まれている複数種類の分子構造の各々の混合割合を取得する取得部（Ｓ１１）と、
取得部により取得された混合割合に基づき、燃焼に伴う煤成分の生じやすさを表わした煤生成指数を算出する煤算出部（Ｓ１２）と、
取得部により取得された混合割合に基づき、燃焼に伴い生じた可溶性有機成分の付着しやすさを表わした付着指数を算出する付着指数算出部（Ｓ１３）と、
煤算出部により算出された煤生成指数、および付着指数算出部により算出された付着指数に基づき、燃焼システムの所定部位に付着する可溶性有機成分の堆積量を推定する堆積量推定部（Ｓ１４）と、
を備え、
燃料に含まれる成分のうち、燃焼前に分解して重合することでアロマ類成分を形成する成分をアロマ類可変成分と呼ぶ場合において、
煤算出部は、
取得部により取得された複数種類の混合割合のうちアロマ類成分の混合割合が多いほど、煤成分が生じやすいことを表わす値に煤生成指数を算出し、
なおかつ、取得部により取得された複数種類の混合割合のうちアロマ類可変成分の混合割合が多いほど、煤成分が生じやすいことを表わす値に煤生成指数を算出し、
煤算出部の算出に用いられるアロマ類可変成分には、ナフテン、側鎖パラフィンおよび直鎖パラフィンが含まれているデポジット推定装置である。
【０００９】
また、開示される第２の発明は、
燃焼システムの燃焼に用いる燃料に含まれている複数種類の分子構造の各々の混合割合を取得する取得部（Ｓ１１）と、
取得部により取得された混合割合に基づき、燃焼に伴う煤成分の生じやすさを表わした煤生成指数を算出する煤算出部（Ｓ１２）と、
取得部により取得された混合割合に基づき、燃焼に伴い生じた可溶性有機成分の付着しやすさを表わした付着指数を算出する付着指数算出部（Ｓ１３）と、
煤算出部により算出された煤生成指数、および付着指数算出部により算出された付着指数に基づき、燃焼システムの所定部位に付着する可溶性有機成分の堆積量を推定する堆積量推定部（Ｓ１４）と、
堆積量推定部により推定された堆積量に応じて、堆積量を低減させるように燃焼システムの作動を制御する制御部（Ｓ１６）と、
を備え、
燃料に含まれる成分のうち、燃焼前に分解して重合することでアロマ類成分を形成する成分をアロマ類可変成分と呼ぶ場合において、
煤算出部は、
取得部により取得された複数種類の混合割合のうちアロマ類成分の混合割合が多いほど、煤成分が生じやすいことを表わす値に煤生成指数を算出し、
なおかつ、取得部により取得された複数種類の混合割合のうちアロマ類可変成分の混合割合が多いほど、煤成分が生じやすいことを表わす値に煤生成指数を算出し、
煤算出部の算出に用いられるアロマ類可変成分には、ナフテン、側鎖パラフィンおよび直鎖パラフィンが含まれている燃焼システム制御装置である。
【００１０】
さて、燃焼システムの排ガス中に含まれる微粒子成分（つまりＰＭ）は煤を主成分としているが、煤のままの状態では粘着性を有していないドライ状態である。このようなドライ状態の煤が、排ガス中に含まれる未燃燃料や潤滑油に取り込まれたり、煤前駆体である多環アロマ成分が燃え残ったりすると、粘着性を有した状態、つまりＳＯＦ成分と呼ばれる可溶性有機成分になる。このＳＯＦ成分が付着して堆積してデポジットとなる。したがって、燃焼に伴い煤成分が生じやすい燃料であるほど、ＳＯＦ成分が多くなるのでデポジット量は多くなる。また、ＳＯＦ成分の粘性が高くなるような燃料であるほど、ＳＯＦ成分が付着堆積しやすくなるのでデポジット量は多くなる。つまり、用いられている燃料について、煤成分が生じやすい燃料であるか否かの情報（つまり煤生成指数）、およびＳＯＦ成分が高粘性となる燃料であるか否かの情報（つまり付着指数）が得られれば、デポジット量を高精度で推定できる筈である。
【００１１】
そして、「煤生成指数および付着指数は、燃料に含まれている複数種類の分子構造の各々の混合割合から推測可能である」との知見を本発明者らは得ている。例えば、煤成分は、多数の直鎖や側鎖をもつパラフィン類成分やナフテン類成分が熱分解やラジカルによる分解を経て重合してアロマ類成分に変化し、そのアロマ類成分が重合、縮合により積層して形成されたものである。そのため、アロマ類成分や、上述の如くアロマ類成分に変化し得る成分（以下、アロマ類可変成分と記載）の混合割合が多い燃料であるほど、煤成分が生じやすい燃料、つまり煤生成指数が高い燃料であると言える。例えば、アロマ類成分の中でも炭素数が多いアロマ類成分の混合割合が多い燃料であるほど、ＳＯＦ成分の揮発性が低くなるため、ＳＯＦ成分の粘性が高くなりやすい燃料、つまり付着指数が高い燃料であると言える。
【００１２】
これらの知見に基づき、上記第１および第２の発明によれば、複数種類の分子構造の各々の混合割合に基づき煤生成指数を算出する。また、上記混合割合に基づき付着指数を算出する。そして、このように算出された両指数に基づきＳＯＦ成分の堆積量（つまりデポジット量）を推定する。そのため、デポジット量を高精度で推定できる。
【図面の簡単な説明】
【００１３】
【図１】本発明の第１実施形態に係る燃焼システム制御装置と、その装置が適用される内燃機関の燃焼システムを説明する図。
【図２】着火遅れ時間の説明図。
【図３】複数の着火遅れ時間、燃えやすさを表わす燃焼環境値の組み合わせである燃焼条件、および各種成分の混合量の関係を説明する図。
【図４】筒内酸素濃度に起因して生じる着火遅れ時間の変化を表す特性線と、燃料の分子構造種との関係を示す図。
【図５】筒内温度に起因して生じる着火遅れ時間の変化を表す特性線と、燃料の分子構造種との関係を示す図。
【図６】着火遅れ時間に基づき特定される特性線と、分子構造種の混合割合との関係を示す図。
【図７】図１に示すマイコンの処理フローであって、デポジット量を推定し、その推定結果に基づき燃焼システムの作動を制御する手順を示すフローチャート。
【図８】第１実施形態において、煤生成指数Ｘを算出する行列式を説明する図。
【図９】第１実施形態において、付着指数Ｙを算出する行列式を説明する図。
【図１０】第１実施形態において、煤生成指数Ｘおよび付着指数Ｙと、デポジット量Ｍとの関係を示す図。
【図１１】本発明の第３実施形態において、煤生成指数Ｘおよび付着指数Ｙと、デポジット量Ｍとの関係を示す図。
【発明を実施するための形態】
【００１４】
以下、図面を参照しながら発明を実施するための複数の形態を説明する。各形態において、先行する形態で説明した事項に対応する部分には同一の参照符号を付して重複する説明を省略する場合がある。各形態において、構成の一部のみを説明している場合は、構成の他の部分については先行して説明した他の形態を参照し適用することができる。
【００１５】
（第１実施形態）
本実施形態に係る燃焼システム制御装置は、図１に示す電子制御装置（つまりＥＣＵ８０）により提供される。ＥＣＵ８０は、マイクロコンピュータ（つまりマイコン８０ａ）や、図示しない入力処理回路および出力処理回路等を備える。マイコン８０ａは、図示しない中央処理装置（つまりＣＰＵ）およびメモリ８０ｂを備える。メモリ８０ｂに記憶された所定のプログラムをＣＰＵが実行することで、マイコン８０ａは、燃焼システムが備える燃料噴射弁１５、燃料ポンプ１５ｐ、ＥＧＲバルブ１７ａ、調温バルブ１７ｄ、および過給調圧機器２６等の作動を制御する。これらの制御により、燃焼システムが備える内燃機関１０での燃焼状態は、所望の状態に制御される。燃焼システムおよびＥＣＵ８０は車両に搭載されたものであり、当該車両は、内燃機関１０の出力を駆動源として走行する。
【００１６】
内燃機関１０は、シリンダブロック１１、シリンダヘッド１２およびピストン１３等を備える。シリンダヘッド１２には、吸気バルブ１４ｉｎ、排気バルブ１４ｅｘ、燃料噴射弁１５および筒内圧センサ２１が取り付けられている。コモンレール１５ｃ等の燃料通路を形成する部分または燃料タンクには、燃料の密度を検出する密度センサ２７、および燃料の動粘度を検出する動粘度センサ２８が取り付けられている。
【００１７】
燃料ポンプ１５ｐは、燃料タンク内の燃料をコモンレール１５ｃへ圧送する。ＥＣＵ８０が燃料ポンプ１５ｐの作動を制御することで、コモンレール１５ｃ内の燃料は、目標圧力Ｐｔｒｇに維持された状態でコモンレール１５ｃに蓄えられる。コモンレール１５ｃは、蓄圧された燃料を各気筒の燃料噴射弁１５へ分配する。燃料噴射弁１５から噴射された燃料は、燃焼室１１ａで吸気と混合して混合気を形成し、混合気はピストン１３により圧縮されて自着火する。要するに、内燃機関１０は圧縮自着火式のディーゼルエンジンであり、燃料には軽油が用いられている。
【００１８】
燃料噴射弁１５は、電磁アクチュエータおよび弁体をボデー内部に収容して構成されている。電磁アクチュエータへの通電をＥＣＵ８０がオンさせると、電磁アクチュエータの電磁吸引力により図示しない背圧室のリーク通路が開弁し、背圧低下に伴い弁体が開弁作動し、ボデーに形成されている噴孔が開弁されて噴孔から燃料が噴射される。上記通電をオフさせると、弁体が閉弁作動して燃料噴射が停止される。
【００１９】
シリンダヘッド１２に形成されている吸気ポート１２ｉｎおよび排気ポート１２ｅｘには、吸気管１６ｉｎおよび排気管１６ｅｘが接続されている。吸気管１６ｉｎおよび排気管１６ｅｘにはＥＧＲ管１７が接続されており、排気の一部であるＥＧＲガスが、ＥＧＲ管１７を通じて吸気管１６ｉｎへ流入（つまり還流）する。ＥＧＲ管１７にはＥＧＲバルブ１７ａが取り付けられている。ＥＣＵ８０がＥＧＲバルブ１７ａの作動を制御することで、ＥＧＲ管１７の開度が制御され、ＥＧＲガスの流量が制御される。
【００２０】
さらに、ＥＧＲ管１７のうちＥＧＲバルブ１７ａの上流部分には、ＥＧＲガスを冷却するＥＧＲクーラ１７ｂ、バイパス管１７ｃおよび調温バルブ１７ｄが取り付けられている。バイパス管１７ｃは、ＥＧＲガスがＥＧＲクーラ１７ｂをバイパスするバイパス流路を形成する。調温バルブ１７ｄは、バイパス流路の開度を調整することで、ＥＧＲクーラ１７ｂを流れるＥＧＲガスと、バイパス流路を流れるＥＧＲガスとの割合を調整し、ひいては、吸気管１６ｉｎへ流入するＥＧＲガスの温度を調整する。ここで、吸気ポート１２ｉｎへ流入する吸気には、吸気管１６ｉｎから流入する外部空気（つまり新気）およびＥＧＲガスが含まれる。したがって、調温バルブ１７ｄによりＥＧＲガスの温度を調整することは、吸気ポート１２ｉｎへ流入する吸気の温度であるインマニ温度を調整することに相当する。
【００２１】
燃焼システムは図示しない過給機を備える。過給機は、排気管１６ｅｘに取り付けられるタービン、および吸気管１６ｉｎに取り付けられるコンプレッサを有する。排気の流速エネルギによりタービンが回転すると、タービンの回転力によりコンプレッサが回転し、コンプレッサにより新気が圧縮つまり過給される。先述した過給調圧機器２６は、タービンの容量を変化させる機器であり、ＥＣＵ８０が過給調圧機器２６の作動を制御することで、タービン容量が調整され、これにより、コンプレッサによる過給圧が制御される。
【００２２】
ＥＣＵ８０には、筒内圧センサ２１、酸素濃度センサ２２、レール圧センサ２３、クランク角センサ２４およびアクセルペダルセンサ２５等、各種センサによる検出信号が入力される。
【００２３】
筒内圧センサ２１は、燃焼室１１ａの圧力（つまり筒内圧）に応じた検出信号を出力する。筒内圧センサ２１は、圧力検出素子に加えて温度検出素子２１ａを有しており、燃焼室１１ａの温度（つまり筒内温度）に応じた検出信号も出力する。酸素濃度センサ２２は、吸気管１６ｉｎに取り付けられ、吸気中の酸素濃度に応じた検出信号を出力する。検出対象となる吸気は、新気とＥＧＲガスが混合したものである。レール圧センサ２３はコモンレール１５ｃに取り付けられており、蓄圧されている燃料の圧力（つまりレール圧）に応じた検出信号を出力する。クランク角センサ２４は、ピストン１３により回転駆動するクランク軸の回転速度であって、単位時間あたりのクランク軸の回転数（つまりエンジン回転数）に応じた検出信号を出力する。アクセルペダルセンサ２５は、車両運転者により踏み込み操作されるアクセルペダルの踏込量（つまりエンジン負荷）に応じた検出信号を出力する。
【００２４】
ＥＣＵ８０は、これらの検出信号に基づき、燃料噴射弁１５、燃料ポンプ１５ｐ、ＥＧＲバルブ１７ａ、調温バルブ１７ｄおよび過給調圧機器２６の作動を制御する。これにより、燃料の噴射開始時期、噴射量、噴射圧、ＥＧＲガス流量、インマニ温度および過給圧が制御される。
【００２５】
燃料噴射弁１５の作動を制御している時のマイコン８０ａは、燃料の噴射開始時期、噴射量、および多段噴射に係る噴射段数を制御する噴射制御部８３として機能する。燃料ポンプ１５ｐの作動を制御している時のマイコン８０ａは、噴射圧を制御する燃圧制御部８４として機能する。ＥＧＲバルブ１７ａの作動を制御している時のマイコン８０ａは、ＥＧＲガス流量を制御するＥＧＲ制御部８５として機能する。調温バルブ１７ｄの作動を制御している時のマイコン８０ａは、インマニ温度を制御するインマニ温度制御部８７として機能する。過給調圧機器２６の作動を制御している時のマイコン８０ａは、過給圧を制御する過給圧制御部８６として機能する。
【００２６】
マイコン８０ａは、燃焼に関する物理量の検出値（つまり燃焼特性値）を取得する燃焼特性取得部８１としても機能する。本実施形態に係る燃焼特性値とは、図２に示す着火遅れ時間ＴＤのことである。図２の上段は、マイコン８０ａから出力されるパルス信号を示す。パルス信号にしたがって燃料噴射弁１５への通電が制御される。具体的には、パルスオンのｔ１時点で通電が開始され、パルスオン期間Ｔｑに通電オンが継続される。要するに、パルスオンのタイミングにより噴射開始時期が制御される。また、パルスオン期間Ｔｑにより噴射期間が制御され、ひいては噴射量が制御される。
【００２７】
図２の中段は、パルス信号にしたがって弁体が開弁作動および閉弁作動した結果生じる、噴孔からの燃料の噴射状態の変化を示す。具体的には、単位時間あたりに噴射される燃料の噴射量（つまり噴射率）の変化を示す。図示されるように、通電開始のｔ１時点から、実際に噴射が開始されるｔ２時点までにはタイムラグが存在する。また、通電終了時点から実際に噴射が停止されるまでにもタイムラグが存在する。実際に噴射が為されている期間Ｔｑ１は、パルスオン期間Ｔｑで制御される。
【００２８】
図２の下段は、噴射された燃料の、燃焼室１１ａでの燃焼状態の変化を示す。具体的には、噴射された燃料と吸気の混合気が自着火燃焼することに伴い生じる、単位時間あたりの熱量（つまり熱発生率）の変化を示す。図示されるように、噴射開始のｔ２時点から、実際に燃焼が開始されるｔ３時点までにはタイムラグが存在する。本実施形態では、通電開始のｔ１時点から燃焼開始のｔ３時点までの時間を着火遅れ時間ＴＤと定義する。
【００２９】
燃焼特性取得部８１は、筒内圧センサ２１で検出される筒内圧の変化に基づき、燃焼開始のｔ３時点を推定する。具体的には、ピストン１３が上死点に達してからクランク角が所定量だけ回転する期間において、筒内圧が急上昇した時期を燃焼開始時期（つまりｔ３時点）と推定する。この推定結果に基づき、着火遅れ時間ＴＤは燃焼特性取得部８１により算出される。さらに燃焼特性取得部８１は、燃焼時の各種状態（つまり燃焼条件）を、燃焼毎に取得する。具体的には、筒内圧、筒内温度、吸気酸素濃度、噴射圧力および混合気流速の少なくとも１つを、燃焼環境値として取得する。
【００３０】
これらの燃焼環境値は、燃料の燃えやすさを表わすパラメータであり、燃焼直前での筒内圧が高いほど、燃焼直前での筒内温度が高いほど、吸気酸素濃度が高いほど、噴射圧力が高いほど、混合気流速が速いほど、混合気が自着火しやすく燃えやすいと言える。燃焼直前での筒内圧および筒内温度として、例えば、燃料噴射弁１５への通電を開始するｔ１時点で検出された値を用いればよい。筒内圧は筒内圧センサ２１により検出され、筒内温度は温度検出素子２１ａにより検出され、吸気酸素濃度は酸素濃度センサ２２により検出され、噴射圧力はレール圧センサ２３により検出される。混合気流速は、燃焼直前における燃焼室１１ａ内での混合気の流速である。この流速は、上記エンジン回転数が速いほど速くなるので、エンジン回転数に基づき算出される。燃焼特性取得部８１は、取得した着火遅れ時間ＴＤを、その燃焼に係る上記燃焼環境値の組み合わせ（つまり燃焼条件）と関連付けてメモリ８０ｂに記憶させる。
【００３１】
マイコン８０ａは、異なる燃焼条件で検出された複数の燃焼特性値に基づき、燃料に含まれている各種成分の混合割合を推定する、混合割合推定部８２としても機能する。例えば、異なる燃焼条件毎の着火遅れ時間ＴＤを図３に示す行列式に代入することで、各種成分の混合量を算出する。なお、算出された各々の混合量を総量で除算することで、各種成分の混合割合が算出される。
【００３２】
図３の左辺にある行列は、ｘ行１列であり、この行列が有する数値は、各種成分の混合量を表わす。各種成分とは、分子構造の種類の違いにより分類される成分である。分子構造の種類には、直鎖パラフィン類、側鎖パラフィン類、ナフテン類およびアロマ類が含まれている。
【００３３】
右辺の左側にある行列は、ｘ行ｙ列でありこの行列が有する数値は、予め実施した試験に基づき定められた定数である。右辺の右側にある行列は、ｙ行１列でありこの行列が有する数値は、燃焼特性取得部８１により取得された着火遅れ時間ＴＤである。例えば、１行１列目の数値は、燃焼環境値の所定の組み合わせからなる燃焼条件ｉの時に取得された着火遅れ時間ＴＤ（条件ｉ）であり、２行１列目の数値は、燃焼条件ｊの時に取得された着火遅れ時間ＴＤ（条件ｊ）である。燃焼条件ｉと燃焼条件ｊとでは、全ての燃焼環境値が異なる値に設定されている。以下の説明では、燃焼条件ｉに係る筒内圧、筒内温度、吸気酸素濃度および噴射圧力を、Ｐ（条件ｉ）、Ｔ（条件ｉ）、Ｏ２（条件ｉ）、Ｐｃ（条件ｉ）とする。燃焼条件ｊに係る筒内圧、筒内温度、吸気酸素濃度および噴射圧力を、Ｐ（条件ｊ）、Ｔ（条件ｊ）、Ｏ２（条件ｊ）、Ｐｃ（条件ｊ）とする。
【００３４】
次に、図４、図５および図６を用いて、図３の行列式に燃焼条件毎の着火遅れ時間ＴＤを代入することで各分子構造種の混合量が算出できる理屈を説明する。
【００３５】
図４に示すように、燃焼に係る混合気に含まれる酸素の濃度（つまり筒内酸素濃度）が高いほど自着火しやすくなるので、着火遅れ時間ＴＤが短くなる。図中の３本の実線（１）（２）（３）は、筒内酸素濃度と着火遅れ時間ＴＤとの関係を示す特性線である。但し、この特性線は燃料に応じて異なる。厳密には、燃料に含まれている各々の分子構造種の混合割合に応じて特性線は異なる。したがって、筒内酸素濃度がＯ２（条件ｉ）の場合の着火遅れ時間ＴＤを検出すれば、いずれの分子構造種であるかを推測できる。特に、筒内酸素濃度がＯ２（条件ｉ）の場合とＯ２（条件ｊ）の場合とで着火遅れ時間ＴＤを比較すれば、より高精度で混合割合を推定できる。
【００３６】
同様にして、図５に示すように、筒内温度が高いほど自着火しやすくなるので、着火遅れ時間ＴＤが短くなる。図中の３本の実線（１）（２）（３）は、筒内温度と着火遅れ時間ＴＤとの関係を示す特性線である。但し、この特性線は燃料に応じて異なる。厳密には、燃料に含まれている各々の分子構造種の混合割合に応じて異なる。したがって、筒内温度がＢ１の場合の着火遅れ時間ＴＤを検出すれば、いずれの分子構造種であるかを推測できる。特に、筒内温度がＴ（条件ｉ）の場合とＴ（条件ｊ）の場合とで着火遅れ時間ＴＤを比較すれば、より高精度で混合割合を推定できる。
【００３７】
同様に噴射圧が高ければ、酸素を取り込みやすく自着火しやすくなるので、着火遅れ時間ＴＤが短くなる。厳密には、燃料に含まれている各々の分子構造種の混合割合に応じて感度が異なる。したがって、噴射圧が異なる場合の着火遅れ時間ＴＤを検出すれば、より高精度で混合割合を推定できる。
【００３８】
また、筒内酸素濃度に係る特性線（図４参照）に対する影響度の高い分子構造種と、筒内温度に係る特性線（図５参照）に対する影響度の高い分子構造種とは異なる。このように、複数の燃焼条件の各々に係る特性線に対して影響度の高い分子構造種は異なる。したがって、複数の燃焼環境値の組み合わせ（つまり燃焼条件）を異なる値にして取得された着火遅れ時間ＴＤの組み合わせに基づけば、例えば図６の如くいずれの分子構造種の混合割合が多いのかを高精度で推定できる。なお、以下の説明では筒内酸素濃度を第１燃焼環境値、筒内温度を第２燃焼環境値と呼び、第１燃焼環境値に係る特性線を第１特性線、第２燃焼環境値に係る特性線を第２特性線と呼ぶ。
【００３９】
図６に例示する分子構造種Ａは、第１燃焼環境値としての筒内酸素濃度に係る特性線（以下、第１特性線と呼ぶ）に対する影響度が高い分子構造種である。また、分子構造種Ｂは、第２燃焼環境値としての筒内温度に係る特性線（以下、第２特性線と呼ぶ）に対する影響度が高い分子構造種であり、分子構造種Ｃは、第３燃焼環境値に係る第３特性線に対する影響度が高い分子構造種である。第１燃焼環境値の変化に対して着火遅れ時間ＴＤの変化が大きく現れるほど、分子構造種Ａが多く混合していると言える。同様にして、第２燃焼環境値の変化に対して着火遅れ時間ＴＤの変化が大きく現れるほど分子構造種Ｂが多く混合しており、第３燃焼環境値の変化に対して着火遅れ時間ＴＤの変化が大きく現れるほど分子構造種Ｃが多く混合していると言える。したがって、異なる燃料（１）（２）（３）の各々に対し、分子構造種Ａ、Ｂ、Ｃの混合割合を推定できる。
【００４０】
次に、燃焼特性取得部８１が実行するプログラムの処理について説明する。この処理は、以下に説明するパイロット噴射が指令される毎に実行される。１燃焼サイクル中に同一の燃料噴射弁１５から複数回噴射（つまり多段噴射）させるように噴射制御する場合がある。これら複数回の噴射のうち、最も噴射量が多く設定された噴射をメイン噴射と呼び、その直前の噴射をパイロット噴射と呼ぶ。
【００４１】
先ず、燃焼特性取得部８１は、上述した通り筒内圧センサ２１の検出値に基づき燃焼開始のｔ３時点を推定して、パイロット噴射に係る着火遅れ時間ＴＤを算出する。次に、複数の燃焼環境値の組み合わせ（つまり燃焼条件）と関連付けて、着火遅れ時間ＴＤをメモリ８０ｂに記憶させる。
【００４２】
具体的には、各燃焼環境値が取り得る数値範囲を複数の領域に区分けしておき、複数の燃焼環境値の領域の組み合わせ予め設定しておく。例えば図３に示す着火遅れ時間ＴＤ（条件ｉ）は、Ｐ（条件ｉ）、Ｔ（条件ｉ）、Ｏ２（条件ｉ）、Ｐｃ（条件ｉ）の領域の組み合わせ時に取得された着火遅れ時間ＴＤを表わす。同様に、着火遅れ時間ＴＤ（条件ｊ）は、Ｐ（条件ｊ）、Ｔ（条件ｊ）、Ｏ２（条件ｊ）、Ｐｃ（条件ｊ）の領域の組み合わせ時に取得された着火遅れ時間ＴＤを表わす。
【００４３】
なお、ユーザが給油することに起因して、燃料タンクに貯留されている燃料に別の燃料が混合した可能性が高い場合に、分子構造種の混合割合が変化したとみなし、推定されていた混合量の値をリセットする。例えば、内燃機関１０の運転停止時に、燃料タンクの燃料残量を検出するセンサにより燃料残量の増大が検出された場合にリセットする。
【００４４】
燃焼特性取得部８１は、着火遅れ時間ＴＤを図３の行列式に代入して、分子構造種毎の混合量を算出する。なお、サンプリング数、つまり行列式の右辺右側の行列の行数に応じて、定数を表わす行列の列数を変更する。或いは、取得されていない着火遅れ時間ＴＤについては、予め設定しておいたノミナル値を着火遅れ時間ＴＤの行列に代入する。このように算出された分子構造種毎の混合量に基づき、分子構造種毎の混合割合を算出する。
【００４５】
マイコン８０ａは、分子構造種毎の混合割合に基づき、燃焼システムの所定部位に付着するＳＯＦ成分の堆積量（つまりデポジット量）を推定する堆積量推定部８８としても機能する。デポジット量Ｍの推定手法については図７〜図１０を用いて後に詳述する。可溶性有機成分であるＳＯＦ成分が付着する上記所定部位の具体例としては、ＥＧＲバルブ１７ａ、ＥＧＲクーラ１７ｂ、調温バルブ１７ｄ、燃料噴射弁１５の噴孔周り、吸気バルブ１４ｉｎ、排気バルブ１４ｅｘ等が挙げられる。要するに、上記所定部位とは、燃焼システムのうち排ガスに晒される部分のことである。
【００４６】
先述した通り、マイコン８０ａは、噴射制御部８３、燃圧制御部８４、ＥＧＲ制御部８５、過給圧制御部８６およびインマニ温度制御部８７としても機能する。これらの制御部は、エンジン回転数、エンジン負荷およびエンジン冷却水温度等に基づき目標値を設定し、制御対象が目標値となるようにフィードバック制御する。或いは、目標値に対応する内容でオープン制御する。なお、特許請求の範囲に記載の「燃焼システム」は、内燃機関１０を備えるとともに上記制御対象を備えて構成されている。
【００４７】
噴射制御部８３は、噴射開始時期、噴射量および噴射段数が目標値となるように図２のパルス信号を設定することで、噴射開始時期、噴射量および噴射段数を制御（つまり噴射制御）する。上記噴射段数とは、先述した多段噴射に係る噴射回数のことである。具体的には、上記目標値に対応するパルス信号のオン時間（つまり通電時間）およびパルスオン立ち上がり時期（つまり通電開始時期）を、マップ上に予め記憶させておく。そして、目標値に対応する通電時間および通電開始時期をマップから取得してパルス信号を設定する。
【００４８】
また、噴射により得られた出力トルクや、ＮＯｘ量およびスモーク量等のエミッション状態値を記憶しておく。そして、次回以降の噴射において、エンジン回転数およびエンジン負荷等に基づき目標値を設定するにあたり、上述の如く記憶された値に基づき、目標値を補正する。要するに、実際の出力トルクやエミッション状態値と、所望する出力トルクやエミッション状態値との偏差をゼロにするよう、目標値を補正してフィードバック制御する。
【００４９】
燃圧制御部８４は、燃料ポンプ１５ｐに吸入される燃料の流量を制御する調量弁の作動を制御する。具体的には、レール圧センサ２３で検出された実レール圧と目標圧力Ｐｔｒｇ（つまり目標値）との偏差に基づき、調量弁の作動をフィードバック制御する。その結果、燃料ポンプ１５ｐによる単位時間当りの吐出量が制御され、実レール圧が目標値となるように制御（つまり燃圧制御）される。
【００５０】
ＥＧＲ制御部８５は、エンジン回転数およびエンジン負荷等に基づき、ＥＧＲ量の目標値を設定する。そして、この目標値に基づき、ＥＧＲバルブ１７ａのバルブ開度を制御（つまりＥＧＲ制御）してＥＧＲ量を制御する。過給圧制御部８６は、エンジン回転数およびエンジン負荷等に基づき、過給圧の目標値を設定する。そして、この目標値に基づき、過給調圧機器２６の作動を制御（つまり過給圧制御）して過給圧を制御する。インマニ温度制御部８７は、外気温度、エンジン回転数およびエンジン負荷等に基づき、インマニ温度の目標値を設定する。そして、この目標値に基づき、調温バルブ１７ｄのバルブ開度を制御（つまりインマニ温度制御）してインマニ温度を制御する。
【００５１】
さらに、上述した各種の制御部により設定される目標値は、混合割合に応じて推定されるデポジット量Ｍに応じて、後述するデポジット低減制御により変更される。この補正をマイコン８０ａが実行する処理手順について、図７を用いて以下に説明する。この処理は、内燃機関１０の運転期間中、所定周期で繰返し実行される。
【００５２】
先ず、図７のステップＳ１０において、燃焼室１１ａで燃焼が生じる直前における燃焼条件、つまり先述した各種の燃焼環境値の各々を取得する。例えば、筒内圧、筒内温度、吸気酸素濃度、噴射圧力および混合気流速の少なくとも１つを、燃焼環境値として取得する。
【００５３】
続くステップＳ１１では、混合割合推定部８２により推定された混合割合を取得する。つまり、図３の左辺に示す分子構造種の各々についての混合割合を取得する。続くステップＳ１２では、ステップＳ１１で取得された混合割合に基づき、燃焼に伴う煤成分の生じやすさを表わした煤生成指数Ｘを算出する。例えば、燃料の単位量当りに含まれる分子構造種毎の混合量（つまり混合割合）を図８に示す行列式に代入することで、煤生成指数Ｘを算出する。例えば、図８に示す行列式に分子構造種毎の混合割合を代入して、燃焼環境値毎の煤生成指数Ｘ００・・・Ｘｘ０を算出する。図８の右辺の左側にある行列は、ｘ行ｙ列であり、この行列が有する数値ｂ００、ｂ０１・・・ｂｘｙは、予め実施した試験に基づき燃焼環境値毎に定められた定数である。右辺の右側にある行列は、ｙ行１列である。算出されたＸベクトルのうち、燃焼環境値に応じた値を最終的な煤生成指数Ｘとする。これらの数値は、混合割合推定部８２により推定された値である。
【００５４】
ここで、異なる燃料のうち、セタン価等の燃料性状が同等の燃料であっても、その燃料に含まれている各種成分の混合割合が異なれば、煤成分の発生のしやすさ（つまり発生度合い）は異なってくる。本実施形態では、煤発生度合いを表した指標を煤生成指数Ｘと呼び、煤生成指数Ｘの値が大きいほど、煤成分発生度合いが大きい。燃料に含まれる分子構造種の中には、煤生成指数Ｘに大きく影響する成分と、あまり影響しない成分とが存在する。このような影響度合いを鑑みて、分子構造種の混合割合に基づき煤生成指数Ｘは算出される。
【００５５】
先述した通り、排気中に含まれるＰＭの主成分は煤であり、煤は、多数のアロマ類成分が熱分解やラジカルによる分解を経て重合し、積層して形成されたものである。この重合反応は、アロマ類成分を含んだ燃料が高温の環境に晒されることに起因して生じる。したがって、煤は、燃焼室１１ａへ噴射された燃料から燃焼の直前に生成される。但し、生成された煤の殆どは、生成直後に燃焼室１１ａで燃焼して消失する。そして、燃え残った煤が燃焼室１１ａから排出される。このように排出された煤が、排気スモーク中のＰＭの主成分である。上記煤生成指数Ｘは、正確には、燃焼室１１ａで燃焼直前に存在している煤の増大しやすさを表わす。煤生成指数Ｘが高い燃料であるほど、燃焼直前に存在している煤の量が多くなるので、燃え残る煤の量が多くなる。
【００５６】
さて、多数の直鎖や側鎖をもつパラフィン類成分やナフテン類成分が、熱分解やラジカルによる分解を経て重合してアロマ類成分に変化する場合がある。このようにアロマ類成分に変化し得る成分をアロマ類可変成分と呼ぶ。そして、アロマ類可変成分の変化により生じたアロマ類成分や、もともと燃料に含まれていたアロマ類成分が、重合、縮合により積層して煤成分を形成する。特にこの重合反応は、アロマ類成分を含んだ燃料が高温の環境に晒されることに起因して生じる。したがって、煤成分は、燃焼室１１ａへ噴射された燃料から燃焼の直前に生成される。したがって、ステップＳ１１で取得した分子構造種毎の混合割合のうち、アロマ類成分の混合割合が多いほど、煤生成指数Ｘが高くなる。また、先述したアロマ類可変成分は、燃焼直前にアロマ類成分に変化し得るので、ステップＳ１１で取得した分子構造種毎の混合割合のうち、アロマ類可変成分の混合割合が多いほど、煤生成指数Ｘが高くなる。
【００５７】
これらの知見に鑑みて、ステップＳ１２では、アロマ類成分およびアロマ類可変成分の混合割合が多いほど、煤生成指数Ｘを大きい値に推定する。詳細には、アロマ類成分の煤生成指数Ｘに対する影響度合を表した重み付け係数は、アロマ類可変成分の煤生成指数Ｘに対する重み付け係数よりも大きく設定されている。
【００５８】
アロマ類可変成分の中でも、アロマ類成分に変化しやすい成分であるほど重み付け係数が大きく設定されている。例えば、アロマ類可変成分の具体例として、ナフテン類成分、側鎖パラフィン類成分および直鎖パラフィン類成分等が挙げられる。そして、ナフテン類成分、側鎖パラフィン類成分および直鎖パラフィン類成分の順に、アロマ類成分に変化しやすいので、この順に上記重み付け係数が大きく設定されている。
【００５９】
ナフテン類成分の中でも、環状構造を２つ以上有する構造のナフテン類成分はアロマ類成分に変化しやすい。そのため、環状構造を２つ以上有する構造のナフテン類成分は、２つ未満のナフテン類成分に比べて上記重み付け係数が大きく設定されている。
【００６０】
側鎖パラフィン類成分の中でも、燃料に含まれている複数種類の成分の平均炭素数よりも炭素数が少ない構造の側鎖パラフィン類成分は、アロマ類成分に変化しやすい。そのため、平均炭素数未満の炭素数を有する側鎖パラフィン類成分は、平均炭素数以上の側鎖パラフィン類成分に比べて上記重み付け係数が大きく設定されている。
【００６１】
図８の行列式への代入に係る分子構造の種類には、直鎖パラフィン類、側鎖パラフィン類およびナフテン類等のアロマ類可変成分と、アロマ類とが含まれている。ナフテン類成分は、環状構造を２つ以上有する構造のナフテン類と、環状構造が２つ未満のナフテン類とに区別して代入される。ナフテン類成分の中でも、環状構造を２つ以上有する構造のナフテン類成分は、特にアロマ類成分に変化しやすい。そのため、環状構造を２つ以上有する構造のナフテン類成分は、２つ未満のナフテン類成分に比べて上記重み付け係数が大きく設定されている。なお、環状構造が２つ未満のナフテン類については、２つ以上のナフテン類に比べてアロマ類に変化しにくいので、行列式への代入を廃止してもよい。
【００６２】
側鎖パラフィン類成分は、炭素数が少ない構造の側鎖パラフィン類と、炭素数が多い構造の側鎖パラフィン類とに区別して代入される。具体的には、燃料に含まれている複数種類の成分の平均炭素数を算出し、その平均炭素数よりも該当する側鎖パラフィン類の炭素数が少ないか否かで上記区別を行う。側鎖パラフィン類成分の中でも、燃料に含まれている複数種類の成分の平均炭素数よりも炭素数が少ない構造の側鎖パラフィン類成分は、特にアロマ類成分に変化しやすい。そのため、平均炭素数未満の炭素数を有する側鎖パラフィン類成分は、平均炭素数以上の側鎖パラフィン類成分に比べて上記重み付け係数が大きく設定されている。なお、炭素数が多い構造の側鎖パラフィン類については、少ない構造の側鎖パラフィン類に比べてアロマ類に変化しにくいので、行列式への代入を廃止してもよい。
【００６３】
図７の説明に戻り、続くステップＳ１３では、ステップＳ１１で取得された混合割合に基づき、燃焼に伴い生じたＳＯＦ成分の付着しやすさを表わした付着指数Ｙを算出する。例えば、燃料の単位量当りに含まれる分子構造種毎の混合量（つまり混合割合）を図９に示す行列式に代入することで、付着指数Ｙを算出する。例えば、図９に示す行列式に分子構造種毎の混合割合を代入して、付着指数Ｙを算出する。図９の右辺の左側にある行列は、１行ｙ列であり、例えばｃ００、ｃ０１・・・ｃ０ｙという数値を有する行列になっている。これらの数値ｃ００、ｃ０１・・・ｃ０ｙは、予め実施した試験に基づき定められた定数である。右辺の右側にある行列は、ｙ行１列であり、この行列が有する数値は、混合割合推定部８２により推定された値である。
【００６４】
ここで、異なる燃料のうち、セタン価等の燃料性状が同等の燃料であっても、その燃料に含まれている各種成分の混合割合が異なれば、ＳＯＦ成分の付着しやすさ（つまり付着度合い）は異なってくる。本実施形態では、付着度合いを表した指標を付着指数Ｙと呼び、付着指数Ｙの値が大きいほど、ＳＯＦ成分の付着度合いが大きい。燃料に含まれる分子構造種の中には、付着指数Ｙに大きく影響する成分と、あまり影響しない成分とが存在する。このような影響度合いを鑑みて、分子構造種の混合割合に基づき付着指数Ｙは算出される。
【００６５】
具体的には、燃料が揮発しやすい性状であるほどＳＯＦ成分の粘性が高くなる。より厳密には、ＳＯＦ成分が揮発しやすい性状であるほどＳＯＦ成分の粘性が高くなる。そして、ＳＯＦ成分の粘性が高くなると付着指数Ｙは高くなり、デポジット量Ｍは増大しやすくなる。
【００６６】
各種成分の混合割合に基づけば、分子構造種の平均炭素数を算出できる。その平均炭素数が多いほど、沸点が高く、揮発性が低い蒸留性状の燃料であるとみなすことができ、例えば、燃料の５０％が蒸発する温度、つまり蒸留性状Ｔ５０を平均炭素数から推定できる。そして、推定した平均炭素数が少ないほど揮発しやすい性状の燃料であるとみなして付着指数Ｙを小さい値にする。
【００６７】
さらに、分子構造種によって、ＳＯＦ成分の粘性への影響度合いは異なる。例えば、多環アロマ、単環アロマ、多環ナフテン、直鎖パラフィン系、側鎖パラフィン系の順にＳＯＦ成分の粘性への影響度合いは大きいので、この順に重み付け係数が大きく設定されている。要するに、分子構造種毎の混合割合と付着指数Ｙとは相関があるので、混合割合から付着指数Ｙを算出できる。
【００６８】
続くステップＳ１４では、ステップＳ１２で算出された煤生成指数ＸおよびステップＳ１３で算出された付着指数Ｙに基づき、デポジット量Ｍを算出する。具体的には、内燃機関１０の運転時間が所定時間経過する毎に、煤生成指数Ｘおよび付着指数Ｙに基づき算出される所定時間毎のデポジット量（つまり単位デポジット量）を積算して、デポジット量Ｍの値を更新していく。このように積算するにあたり、ステップＳ１０で取得した燃焼条件の履歴に応じて積算する値を変更させてもよい。例えば、ＥＧＲバルブ１７ａに付着するデポジット量に関し、単位時間当りにＥＧＲ管１７を通過するＥＧＲ量が多いほど、単位デポジット量を増量させるように変更して積算する。或いは、燃料噴射弁１５やＥＧＲバルブ１７ａに付着するデポジット量に関し、筒内温度が低温であるほど揮発量が少ないとみなして、単位デポジット量を増量させるように変更して積算する。或いは、酸素濃度が少ない燃焼条件で燃焼させた場合にはＳＯＦ成分の発生量が少なくなるので、上述した単位デポジット量を少なくするように補正して、積算してもよい。
【００６９】
図１０の横軸は煤生成指数Ｘ、縦軸は付着指数Ｙであり、両指数の値が大きいほど、図中の矢印に示すようにデポジット量Ｍは多くなる。したがって、図１０に示すデポジット量Ｍと両指数との関係を予め試験等により取得してマイコン８０ａにマップ等の状態で記憶させておき、ステップＳ１４において、上記マップを参照して両指数からデポジット量Ｍを算出すればよい。
【００７０】
図１０中の境界線Ｌ１は、煤が生成される下限範囲を示す。この境界線Ｌ１よりも両指数が小さい側の範囲では、単位デポジット量はゼロであるとみなす。そして、境界線Ｌ１よりも両指数が大きい側の範囲において、煤生成指数Ｘの値が大きいほど、また、付着指数Ｙの値が大きいほど、デポジット量Ｍは多い値に算出される。要するに、両指数が大きいほどデポジット量Ｍは多くなる。そして、煤生成指数Ｘの値が大きくても、付着指数Ｙの値が小さければデポジット量Ｍは少なくなり、付着指数Ｙの値が大きくても、煤生成指数Ｘの値が小さければデポジット量Ｍは少なくなる。境界線Ｌ１よりも両指数が大きい側の範囲において、Ｚ＝ａ・Ｘ・Ｙの演算式に基づき単位デポジット量Ｚを算出する。なお演算式中の「ａ」は、先述した燃焼条件の履歴、ＥＧＲ量、筒内温度等の環境条件に応じて設定される係数である。
【００７１】
例えば、アロマ類成分の混合割合が多い燃料は煤生成指数Ｘが大きい。そして、アロマ類成分のうち、炭素数の多いアロマ類成分は炭素数の少ないアロマ類成分に比べて付着指数Ｙが大きい。つまり、炭素数の多いアロマ類成分が多いほど、煤生成指数Ｘおよび付着指数Ｙの両方が大きくなり、デポジット量Ｍは多くなる傾向にある。具体的には、燃料に含まれている複数種類の成分の平均炭素数よりも炭素数が多いアロマ類成分が燃料に多く含まれているほど、デポジット量Ｍは多くなる傾向にある。
【００７２】
例えば、炭素数の多い直鎖パラフィン類成分が燃料に多く含まれているほど揮発しにくくなり粘性が高くなるので、煤生成指数Ｘはそれほど大きくはならないが、付着指数Ｙが大きくなり、デポジット量Ｍが多くなる傾向にある。
【００７３】
続くステップＳ１５では、デポジット量Ｍが、予め記憶させておいた所定量ＴＨ未満であるか否かを判定する。デポジット量Ｍが所定量ＴＨ未満であると判定された場合、図７の処理を終了し、噴射制御部８３、燃圧制御部８４、ＥＧＲ制御部８５、過給圧制御部８６およびインマニ温度制御部８７による先述した制御（つまり通常制御）がそのまま継続される。
【００７４】
一方、デポジット量Ｍが所定量ＴＨ未満でないと判定された場合、続くステップＳ１６において、デポジット量Ｍを低減させるように、以下に説明するデポジット低減制御を実行する。例えば内燃機関１０を停止させた直後に、ＥＧＲバルブ１７ａを開閉作動させる。これにより、ＥＧＲバルブ１７ａに付着したデポジットが振り落とされて、デポジット量の低減が図られる。或いは、内燃機関１０を停止させた直後にＥＧＲバルブ１７ａの開閉作動を常時実行させる場合において、その開閉作動の回数を増量する。
【００７５】
或いは、噴射制御部８３、燃圧制御部８４、ＥＧＲ制御部８５、過給圧制御部８６およびインマニ温度制御部８７の少なくとも１つにおいて、通常制御に係る各種制御量の目標値を、煤成分が低減する側に補正する。例えば、ＥＧＲ制御部８５に係るＥＧＲ量の目標値を低下させて、実ＥＧＲ量を減少させる。或いは、インマニ温度制御部８７に係るインマニ温度の目標値を低下させて、実インマニ温度を低下させる。これによれば、例えばＥＧＲクーラ１７ｂの製品寿命を長くできる。
【００７６】
続くステップＳ１７では、分子構造種の混合割合の情報である燃料情報と、デポジット低減制御の履歴である制御履歴とをマイコン８０ａに記憶させる。例えば、給油する毎に変化する分子構造種の混合割合を記録するとともに、その記録と関連付けて制御履歴を記録する。
【００７７】
なお、ステップＳ１１の処理を実行している時のマイコン８０ａは、特許請求の範囲に記載の「取得部」に相当する。ステップＳ１２、Ｓ１３の処理を実行している時のマイコン８０ａは、それぞれ、特許請求の範囲に記載の「煤算出部」「付着指数算出部」に相当する。ステップＳ１４の処理を実行している時のマイコン８０ａは、特許請求の範囲に記載の「堆積量推定部」に相当する。ステップＳ１６、Ｓ１７の処理を実行している時のマイコン８０ａは、特許請求の範囲に記載の「制御部」に相当する。そして、マイコン８０ａを備えるＥＣＵ８０により、デポジット推定装置が提供される。
【００７８】
以上に説明した通り、本実施形態では、ステップＳ１１、Ｓ１２、Ｓ１３、Ｓ１４による取得部、煤算出部、付着指数算出部および堆積量推定部を備える。取得部は、燃料に含まれている複数種類の分子構造の各々の混合割合を取得する。煤算出部は、取得部により取得された混合割合に基づき、燃焼に伴う煤成分の生じやすさを表わした煤生成指数Ｘを算出する。付着指数算出部は、取得部により取得された混合割合に基づき、燃焼に伴い生じたＳＯＦ成分の付着しやすさを表わした付着指数Ｙを算出する。堆積量推定部は、煤生成指数Ｘおよび付着指数Ｙに基づき、燃焼システムの所定部位に付着するＳＯＦ成分の堆積量を推定する。
【００７９】
このように、本実施形態によれば、取得部、煤算出部および付着指数算出部を備えるので、複数種類の分子構造の各々の混合割合に基づき、煤生成指数Ｘおよび付着指数Ｙの算出が実現可能になる。その上で、本実施形態では堆積量推定部を備えるので、デポジット量Ｍを高精度で推定できる。
【００８０】
さらに本実施形態では、ステップＳ１３による付着指数算出部は、複数種類の分子構造の各々の混合割合が、燃料の揮発性が低くなるような値の組み合わせであるほど付着指数Ｙを大きい値に算出する。また、上記混合割合が、燃料の平均炭素数が多くなるような値の組み合わせであるほど付着指数Ｙを大きい値に算出する。
【００８１】
ここで、上記混合割合と燃料の平均炭素数とは相関があり、平均炭素数と蒸留性状（つまり揮発性）とについても相関があり、燃料の揮発性が低いほどＳＯＦ成分の粘着性が高くなる、との知見を本発明者らは得ている。したがって、上記混合割合の組み合わせが、燃料の揮発性が低くなる値や平均炭素数が多くなる値である場合に、付着指数Ｙを大きい値にする本実施形態によれば、付着指数Ｙを精度良く推定でき、ひいてはデポジット量Ｍを高精度で推定できる。
【００８２】
また、上記混合割合と動粘度とは相関がある。よって、上記混合割合が、燃料の動粘度が高くなるような値の組み合わせであるほど付着指数Ｙを大きい値に算出する本実施形態によれば、付着指数Ｙを精度良く推定でき、ひいてはデポジット量Ｍを高精度で推定できる。
【００８３】
さらに本実施形態では、ステップＳ１２による煤算出部は、燃料に含まれるアロマ類成分の混合割合が多いほど煤生成指数Ｘを大きい値に算出する。煤成分は、多数の直鎖や側鎖をもつパラフィン類成分やナフテン類成分が分解を経て重合することや、アロマ類成分が重合、縮合により多環化することで形成されるものである。したがって、アロマ類成分の混合割合が多いほど煤生成指数Ｘを大きくする本実施形態によれば、デポジット量Ｍを高精度で推定できる。なお、上記分解には熱分解やラジカルによる分解等があり、厳密には、熱分解が生じた後に、ラジカルによる分解が生じる。
【００８４】
ここで、燃焼室１１ａへ噴射された燃焼前の燃料は、高温環境に晒されることに起因して分子構造が変化する。その変化の１つに、以下に説明するアロマ類可変成分が、熱分解やラジカルによる分解を経た後、重合してアロマ類成分へ変化することが挙げられる。アロマ類可変成分の具体例としてはナフテン類やパラフィン類等が挙げられる。アロマ類は不飽和結合を有した環状構造であるが、このような構造にアロマ類可変成分は変化する。
【００８５】
例えば、ナフテン類は環状構造であるものの不飽和結合を有していない。このようなナフテン類であっても、以下に説明するようにアロマ類に変化する可能性がある。すなわち、熱分解等により原子同士の結合が部分的に切れ、さらに水素引き抜き反応により水素が引き抜かれることでその切れた箇所が別の箇所に結合し、その結果、不飽和結合を有した環状構造、つまりアロマ類に変化する可能性がある。また、パラフィン類は環状構造を有していないが、同様に分解して重合することで、不飽和結合を有した環状構造、つまりアロマ類に変化する可能性がある。
【００８６】
さて、燃焼室１１ａでは、燃焼直前にアロマ類成分が重合して煤成分を形成し、その煤成分の大半が燃焼により消失する。そして、煤成分が未燃燃料や潤滑油に取り込まれたり、煤前駆体である多環アロマ成分が燃え残ったりするとＳＯＦ成分となる。したがって、燃料に含まれているアロマ類成分が多いほどＳＯＦ成分は多くなる。
【００８７】
しかし、上述したように、アロマ類可変成分は燃焼直前にアロマ類成分に変化し得るので、常温の状態ではアロマ類成分が少ない燃料であっても、燃焼直前にはアロマ類成分が多くなっている場合がある。このことは、燃料に含まれているアロマ類成分量が同じであっても、アロマ類可変成分量が異なればＳＯＦ成分量、つまりデポジット量Ｍは異なってくることを意味する。
【００８８】
以上の知見に基づき、本実施形態では、ステップＳ１２による煤算出部において、燃料に含まれるアロマ類可変成分の混合割合が多いほど煤生成指数Ｘを大きい値に算出する。そのため、燃焼前に生じる燃料の分子構造変化をも考慮して煤生成指数Ｘが推定されるので、デポジット量Ｍを高精度で推定できる。
【００８９】
さらに本実施形態では、煤生成指数Ｘの推定に用いるアロマ類可変成分には、ナフテン類成分が少なくとも含まれている。各種のアロマ類可変成分の中でも特にナフテン類成分はアロマ類成分に変化しやすい。したがって、煤生成指数Ｘの推定に用いるアロマ類可変成分量にナフテン類成分量を含ませる本実施形態によれば、煤生成指数Ｘの推定精度を向上できる。
【００９０】
さらに本実施形態では、煤生成指数Ｘの推定に用いるナフテン類成分には、環状構造を２つ以上有する構造のナフテン類成分が少なくとも含まれている。ナフテン類成分の中でも特に環状構造を２つ以上有する構造のナフテン類成分は、アロマ類成分に変化しやすい。したがって、煤生成指数Ｘ推定に用いるアロマ類可変成分量に、環状構造を２つ以上有する構造のナフテン類成分を含ませる本実施形態によれば、煤生成指数Ｘの推定精度を向上できる。
【００９１】
さらに本実施形態では、煤生成指数Ｘの推定に用いるアロマ類可変成分には、側鎖パラフィン類成分が少なくとも含まれている。各種のアロマ類可変成分の中でも特にナフテン類成分はアロマ類成分に変化しやすい。したがって、煤生成指数Ｘの推定に用いるアロマ類可変成分量に側鎖パラフィン類成分量を含ませる本実施形態によれば、煤生成指数Ｘの推定精度を向上できる。
【００９２】
さらに本実施形態では、煤生成指数Ｘの推定に用いる側鎖パラフィン類成分には、燃料に含まれている複数種類の成分の平均炭素数よりも炭素数が少ない構造の側鎖パラフィン類成分が少なくとも含まれている。側鎖パラフィン類成分の中でも特に炭素数が少ない構造の側鎖パラフィン類成分は、アロマ類成分に変化しやすい。したがって、煤生成指数Ｘの推定に用いるアロマ類可変成分量に、平均炭素数よりも炭素数が少ない構造の側鎖パラフィン類成分を含ませる本実施形態によれば、デポジット量Ｍの推定精度を向上できる。
【００９３】
さらに本実施形態では、堆積量推定部により推定された堆積量つまりデポジット量Ｍに応じて、堆積量を低減させるように燃焼システムの作動を制御（つまり低減制御）する制御部を備える。これによれば、高精度で推定されたデポジット量Ｍに基づき低減制御を実行するので、低減制御の過不足を抑制できる。
【００９４】
さらに本実施形態では、燃焼特性取得部８１および混合割合推定部８２を備える。燃焼特性取得部８１は、内燃機関１０の燃焼に関する物理量の検出値を燃焼特性値として取得する。混合割合推定部８２は、異なる燃焼条件で検出された複数の燃焼特性値に基づき、燃料に含まれている各種成分の混合割合を推定する。
【００９５】
ここで、全く同じ燃料を燃焼させても、その時の筒内圧や筒内温度等の燃焼条件が異なれば、着火遅れ時間や熱発生量等の燃焼特性値は異なってくる。例えば、図４の燃料（１）は、筒内酸素濃度が多いといった燃焼条件であるほど、着火遅れ時間ＴＤ（燃焼特性値）は短くなる。そして、燃焼条件の変化に対する燃焼特性値の変化の度合い、つまり図４の実線に示す特性線は、分子構造種の混合割合が互いに異なる燃料（１）（２）（３）の各々で、異なってくる。この点を鑑みた本実施形態では、異なる燃焼条件で検出された複数の着火遅れ時間ＴＤ（燃焼特性値）に基づき、燃料に含まれている分子構造種の混合割合を推定するので、燃料の性状をより正確に把握できるようになる。
【００９６】
さらに本実施形態では、燃焼条件は、複数種類の燃焼環境値の組み合わせにより特定される条件である。つまり、複数種類の燃焼環境値各々について、燃焼環境値の値が異なる燃焼時の燃焼特性値を取得する。これによれば、同一種類の燃焼環境値についてその燃焼環境値の値が異なる燃焼時の燃焼特性値を取得し、それらの燃焼条件および燃焼特性値に基づき混合割合を推定する場合に比べて、混合割合を高精度で推定できる。
【００９７】
さらに本実施形態では、燃焼条件に係る複数種類の燃焼環境値には、筒内圧、筒内温度、吸気酸素濃度および燃料噴射圧力の少なくとも１つが含まれている。これらの燃焼環境値は、燃焼状態に与える影響が大きいので、これらの条件が異なる燃焼時の燃焼特性値を用いて混合割合を推定する本実施形態によれば、混合割合を精度良く推定できる。
【００９８】
さらに本実施形態では、燃焼特性値は、燃料噴射を指令してから自着火するまで着火遅れ時間ＴＤである。着火遅れ時間ＴＤは、各種成分の混合割合の影響を大きく受けるので、着火遅れ時間ＴＤに基づき混合割合を推定する本実施形態によれば、混合割合を精度良く推定できる。
【００９９】
さらに本実施形態では、燃焼特性取得部８１は、メイン噴射の前に噴射（パイロット噴射）された燃料の燃焼に関する燃焼特性値を取得する。メイン噴射の燃料が燃焼すると、筒内温度が高くなるので、メイン噴射後の燃料が燃焼しやすくなる。そのため、燃料の混合割合の違いに起因した燃焼特性値の変化が現れにくくなる。これに対し、メイン噴射の前に噴射（パイロット噴射）された燃料は、メイン燃焼の影響を受けないので、混合割合の違いに起因した燃焼特性値の変化が現れやすくなる。よって、燃焼特性値に基づき混合割合を推定するにあたり、その推定精度を向上できる。
【０１００】
（第２実施形態）
上記第１実施形態では、混合割合推定部８２が、複数の燃焼特性値に基づき各種成分の混合割合を推定している。これに対し本実施形態では、燃料の一般性状をセンサ（つまり性状センサ）で検出し、その検出結果に基づき上記混合割合を推定する。
【０１０１】
上記性状センサの具体例としては、密度センサ２７、および動粘度センサ２８等が挙げられる。密度センサ２７は、例えば固有振動周期測定法に基づいて燃料の密度を検出する。動粘度センサ２８は、例えば細管粘度計や、細線加熱法に基づく動粘度計であり、燃料タンク内の燃料の動粘度を検出する。なお、密度センサ２７及び動粘度センサ２８は、ヒータを備えており、ヒータにより所定温度に燃料を加熱した状態で燃料の密度及び動粘度をそれぞれ検出する。
【０１０２】
ここで、本発明者らは、燃料の特定の性状パラメータ、つまり先述した中間パラメータが、燃料組成に含まれる各分子構造の物理量に相関があること、各性状パラメータについては、性状パラメータの種別ごとに分子構造に対する感度が異なることに着目した。つまり、燃料において分子構造が異なると分子間の結合力、構造による立体障害や相互作用などが相違する。また、燃料には複数種の分子構造が含まれ、その混合割合もまちまちである。この場合、分子構造ごとに性状パラメータに寄与する感度が異なると考えられるため、分子構造量に依存して性状パラメータの値が変化する。
【０１０３】
そこで本発明者らは、性状パラメータと分子構造とについて相関式を構築した。この相関式は、複数の性状パラメータに対する複数の分子構造量の依存度を示す感度係数を用い、複数の分子構造量に感度係数を反映することで複数の性状パラメータを導出する性状算出モデルの演算式である。相関式において、上記性状センサにより検出された値を性状パラメータの値として入力することで、燃料組成に含まれる分子構造量の算出が可能となる。
【０１０４】
また、低位発熱量は、燃料の動粘度及び密度と相関があることから、その相関を示すマップや演算式を用いることで、動粘度及び密度に基づいて算出することが可能である。このようにして算出された低位発熱量を、相関式に入力する性状パラメータとしてもよい。
【０１０５】
また、燃料に含まれている水素量と炭素量との比（つまりＨＣ比）は、低位発熱量と相関があることから、その相関を示すマップや演算式を用いることで、低位発熱量に基づいてＨＣ比を算出することが可能である。このようにして算出されたＨＣ比を、相関式に入力する性状パラメータとしてもよい。その他、性状パラメータとして、セタン価や、蒸留性状に関するパラメータを用いることも可能である。
【０１０６】
以上により、本実施形態によれば、燃料の性状を示す複数の性状パラメータを取得する。そして、複数の性状パラメータと燃料における複数の分子構造量との相関を定義した相関データを用い、取得した複数の性状パラメータの取得値に基づいて複数の分子構造量、つまり分子構造種毎の混合割合を推定する。そのため、筒内圧センサ２１の検出値を用いること無く、性状センサの検出値を用いて、デポジット量Ｍの推定に用いる分子構造種の混合割合または中間パラメータを取得できる。
【０１０７】
（第３実施形態）
上記第１実施形態では、煤生成指数Ｘおよび付着指数Ｙに基づきデポジット量Ｍを算出するにあたり、図１０に示すように、煤が生成される下限範囲の境界線を１本の境界線Ｌ１で定義している。これに対し、本実施形態では、図１１に示すように、煤が生成される下限範囲を４本の境界線Ｌ１、Ｌ２、Ｌ３、Ｌ４で定義している。境界線Ｌ１は図１０の境界線Ｌ１と同じである。境界線Ｌ２は付着指数Ｙの下限値を示し、煤生成指数Ｘの値に拘らずに設定された値である。境界線Ｌ３は煤生成指数Ｘの下限値を示し、付着指数Ｙの値に拘らずに設定された値である。ここで、煤生成指数Ｘが小さく、かつ、付着指数Ｙが大きくなる燃料は存在し得ない。境界線Ｌ４は、このような存在し得ない領域を下限範囲の境界として設定するものである。
【０１０８】
以上により、本実施形態によれば、４種類の技術的意味のある境界線Ｌ１、Ｌ２、Ｌ３、Ｌ４により、デポジット量Ｍの下限範囲を設定するので、煤生成指数Ｘおよび付着指数Ｙに基づきデポジット量Ｍを算出するにあたり、その算出精度を向上できる。
【０１０９】
（第４実施形態）
上記第１実施形態では、図７のステップＳ１３における付着指数算出部において、分子構造種毎の混合割合に基づき付着指数Ｙを算出する。これに対し本実施形態では、動粘度センサ２８により検出された検出値に基づき付着指数Ｙを算出する。検出された動粘度が高いほど、ＳＯＦ成分は付着しやすくなるので付着指数Ｙを高い値に算出する。例えば、図９の演算式に替えて、動粘度を変数とした演算式に、動粘度センサ２８による検出値を代入して付着指数Ｙを算出する。なお、煤生成指数Ｘについては第１実施形態と同様にして混合割合に基づき算出する。そして、両指数からデポジット量Ｍを算出する手法についても、第１実施形態と同様である。
【０１１０】
以上により、本実施形態によれば、ステップＳ１３における付着指数算出部において、動粘度センサ２８で検出された燃料の動粘度が高いほど、付着指数Ｙを高い値に算出する。動粘度と付着指数Ｙとの相関は高いので、本実施形態によっても、上記第１実施形態と同様にして付着指数Ｙを精度良く算出でき、ひいてはデポジット量Ｍを精度良く算出できる。
【０１１１】
（他の実施形態）
以上、発明の好ましい実施形態について説明したが、発明は上述した実施形態に何ら制限されることなく、以下に例示するように種々変形して実施することが可能である。各実施形態で具体的に組合せが可能であることを明示している部分同士の組合せばかりではなく、特に組合せに支障が生じなければ、明示してなくとも実施形態同士を部分的に組み合せることも可能である。
【０１１２】
図９に示す実施形態では、分子構造種毎の混合割合を演算式に代入して付着指数Ｙを算出する。これに対し、分子構造種毎の混合割合から蒸留性状Ｔ５０や動粘度等の中間パラメータを推定し、その推定値を演算式に代入して付着指数Ｙを算出するように演算式を設定してもよい。
【０１１３】
上記第４実施形態では、動粘度センサ２８の検出値に基づき付着指数Ｙを算出しているが、密度センサ２７等の他のセンサにより検出された燃料性状に基づき付着指数Ｙを算出してもよい。或いは、分子構造種毎の混合割合と動粘度とは相関があることに着目し、混合割合に基づき動粘度を推定し、その推定値に基づき付着指数Ｙを算出してもよい。
【０１１４】
図２に示す上記実施形態では、通電開始のｔ１時点から燃焼開始のｔ３時点までの時間を着火遅れ時間ＴＤと定義している。これに対し、噴射開始のｔ２時点から燃焼開始のｔ３時点までの時間を着火遅れ時間ＴＤと定義してもよい。噴射開始のｔ２時点は、噴射開始に伴いレール圧等の燃圧に変化が生じた時期を検出し、その検出時期に基づき推定すればよい。
【０１１５】
図１に示す燃焼特性取得部８１は、燃焼に関する物理量の検出値（つまり燃焼特性値）として、着火遅れ時間ＴＤを取得している。これに対し、熱発生率の変化を表わす波形や、該当する燃料の燃焼で発生した熱量（熱発生量）等を燃焼特性値として取得してもよい。また、着火遅れ時間ＴＤ、熱発生率の波形、および熱発生量等、複数種類の燃焼特性値に基づき、各種成分の混合割合を推定してもよい。例えば、図３の右辺左側の行列（定数）を、複数種類の燃焼特性値に対応した値に設定しておき、図３の右辺右側の行列に、複数種類の燃焼特性値を代入して混合割合を推定する。
【０１１６】
図３の例では、複数の着火遅れ時間ＴＤの各々について、全ての燃焼環境値が異なるように燃焼条件が設定されている。つまり、燃焼環境値の所定の組み合わせからなる燃焼条件ｉ、ｊ、ｋ、ｌ（図３参照）の各々について、筒内圧は全て異なる値Ｐ（条件ｉ）、Ｐ（条件ｊ）、Ｐ（条件ｋ）、Ｐ（条件ｌ）に設定されている。同様に、筒内温度Ｔ、吸気酸素濃度Ｏ２および噴射圧力Ｐｃも全て異なる値に設定されている。これに対し、異なる燃焼条件の各々において、少なくとも１つの燃焼環境値の値が異なっていればよい。例えば燃焼条件ｉ、ｊの各々において、筒内温度Ｔ、吸気酸素濃度Ｏ２および噴射圧力Ｐｃを同じ値に設定し、筒内圧だけを異なる値Ｐ（条件ｉ）、Ｐ（条件ｊ）に設定してもよい。
【０１１７】
上述した実施形態では、メイン噴射の直前に噴射（パイロット噴射）された燃料の燃焼に関する燃焼特性値を取得している。これに対し、メイン噴射の後に噴射された燃料の燃焼に関する燃焼特性値を取得してもよい。メイン噴射後の噴射の具体的例として、アフター噴射やポスト噴射が挙げられる。また、メイン噴射の前に複数回噴射する多段噴射を実施する場合には、初回に噴射された燃料の燃焼に関する燃焼特性値を取得すれば、メイン燃焼の影響を大きく受けずに済むので望ましい。
【０１１８】
上述した実施形態では、筒内圧センサ２１の検出値に基づき燃焼特性値を取得している。これに対し、筒内圧センサ２１を備えていない構成において、回転角センサの回転変動（つまり微分値）に基づき燃焼特性値を推定してもよい。例えば、パイロット燃焼に起因して微分値が既定の閾値を超えた時期をパイロット着火時期として推定できる。また、微分値の大きさからパイロット燃焼量を推定できる。
【０１１９】
図１に示す実施形態では、筒内温度は温度検出素子２１ａにより検出されているが、筒内圧センサ２１により検出された筒内圧に基づき推定してもよい。具体的には、筒内温度を、筒内圧力、シリンダ容積、シリンダ内のガス重量、ガス定数から演算して推定する。
【０１２０】
図７に示す制御では、ステップＳ１４による堆積量推定部が推定したデポジット量Ｍに応じて、デポジット量Ｍを低減させるように燃焼システムの作動を制御するデポジット低減制御をステップＳ１６で実施している。これに対し、ステップＳ１６による制御部を廃止してもよい。この場合、デポジット量Ｍが所定量ＴＨ以上であると判定された場合、ステップＳ１７による燃料情報等の記録を実施するとともに、警告音や表示により運転者に異常を報知することが望ましい。
【０１２１】
ＥＣＵ８０（燃焼システム制御装置）が提供する手段および／または機能は、実体的な記憶媒体に記録されたソフトウェアおよびそれを実行するコンピュータ、ソフトウェアのみ、ハードウェアのみ、あるいはそれらの組合せによって提供することができる。例えば、燃焼システム制御装置がハードウェアである回路によって提供される場合、それは多数の論理回路を含むデジタル回路、またはアナログ回路によって提供することができる。
【符号の説明】
【０１２２】
８０…ＥＣＵ（デポジット推定装置）、Ｓ１１…取得部、Ｓ１２…煤算出部、Ｓ１３…付着指数算出部、Ｓ１４…堆積量推定部。
【図１】
0006436064.tif 000002
【図２】
0006436064.tif 000003
【図３】
0006436064.tif 000004
【図４】
0006436064.tif 000005
【図５】
0006436064.tif 000006
【図６】
0006436064.tif 000007
【図７】
0006436064.tif 000008
【図８】
0006436064.tif 000009
【図９】
0006436064.tif 000010
【図１０】
0006436064.tif 000011
【図１１】
0006436064.tif 000012