特許 7632266 バイオ燃料混合比率推定装置およびバイオ燃料混合比率推定方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2025-02-10

(45)【発行日】2025-02-19

(54)【発明の名称】バイオ燃料混合比率推定装置およびバイオ燃料混合比率推定方法

(51)【国際特許分類】

   F02D 45/00 20060101AFI20250212BHJP

   F02D 19/08 20060101ALI20250212BHJP

【ＦＩ】

F02D45/00 369

F02D19/08 Z

F02D45/00 362

F02D45/00 372

【請求項の数】 9

(21)【出願番号】P 2021209513

(22)【出願日】2021-12-23

(65)【公開番号】P2023094189

(43)【公開日】2023-07-05

【審査請求日】2024-03-14

(73)【特許権者】

【識別番号】000003218

【氏名又は名称】株式会社豊田自動織機

(74)【代理人】

【識別番号】110001195

【氏名又は名称】弁理士法人深見特許事務所

(72)【発明者】

【氏名】東 泰弘

(72)【発明者】

【氏名】大坪 康彦

【審査官】戸田 耕太郎

(56)【参考文献】

【文献】特開２００９－１３８５５６（ＪＰ，Ａ）

【文献】特開２０１１－２３１７４７（ＪＰ，Ａ）

【文献】特開２００７－１２６９９８（ＪＰ，Ａ）

【文献】米国特許出願公開第２００７／０１６３５４２（ＵＳ，Ａ１）

【文献】特開２００９－２０３９２４（ＪＰ，Ａ）

【文献】特開２０１１－１４９２９７（ＪＰ，Ａ）

【文献】特開２００８－６９７１３（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｆ０２Ｄ ４５／００

Ｆ０２Ｄ １９／０８

(57)【特許請求の範囲】

【請求項1】

内燃機関に供給される燃料に混合されたバイオ燃料の混合比率を推定する、バイオ燃料混合比率推定装置であって、
前記内燃機関のクランク角を検出するクランク角センサと、
前記クランク角に基づいてクランク角加速度を算出するクランク角加速度算出部と、
前記クランク角加速度のうち正の前記クランク角加速度を積分した、角加速度積分値を算出する積分値算出部と、
基準燃料における前記角加速度積分値である基準値を記憶する基準値記憶部と、
前記積分値算出部で算出した前記角加速度積分値と前記基準値との差である積分値偏差を算出する偏差算出部と、
前記積分値偏差に基づいて、発熱量差を算出する発熱量差算出部と、
前記発熱量差に基づいて、前記混合比率を算出する混合比率算出部と、を備えるバイオ燃料混合比率推定装置。

【請求項2】

前記内燃機関の回転速度と燃料噴射量と前記積分値偏差とをパラメータとした発熱量差マップを備え、
前記発熱量差算出部は、前記発熱量差マップを用いて前記発熱量差を算出する、請求項１に記載のバイオ燃料混合比率推定装置。

【請求項3】

前記内燃機関の燃料噴射量と前記発熱量差をパラメータとした混合比率マップを備え、
前記混合比率算出部は、前記混合比率マップを用いて前記混合比率を算出する、請求項１または請求項２に記載のバイオ燃料混合比率推定装置。

【請求項4】

前記内燃機関の回転速度の変化量が所定範囲内にあるか否かを判定する定常状態判定部を備え、
前記定常状態判定部で前記変化量が前記所定範囲内にあると判定されたとき、前記積分値算出部が前記角加速度積分値を算出する、請求項１から請求項３のいずれか一項に記載のバイオ燃料混合比率推定装置。

【請求項5】

燃料タンクへの給油が行われた後、最初に前記内燃機関が運転されたとき、前記積分値算出部が前記角加速度積分値を算出する、請求項１から請求項４のいずれか一項に記載のバイオ燃料混合比率推定装置。

【請求項6】

前記内燃機関の筒内における空気過剰率が１より大きい所定値以上のとき、前記積分値算出部が前記角加速度積分値を算出する、請求項１から請求項５のいずれか一項に記載のバイオ燃料混合比率推定装置。

【請求項7】

前記基準燃料は、バイオ燃料が混合されていない燃料であり、
前記基準値記憶部は、バイオ燃料が混合されていない燃料における前記角加速度積分値である基準値を記憶している、請求項１から請求項６のいずれか一項に記載のバイオ燃料混合比率推定装置。

【請求項8】

内燃機関に供給される燃料に混合されたバイオ燃料の混合比率を推定する、バイオ燃料混合比率推定方法であって、
前記内燃機関のクランク角からクランク角加速度を算出するステップと、
正のクランク角加速度を積分し、角加速度積分値を算出するステップと、
基準燃料における前記角加速度積分値である基準値を読み出すステップと、
算出した前記角加速度積分値と読み出した前記基準値との差である積分値偏差を算出するステップと、
前記積分値偏差を用いて発熱量差を算出するステップと、
前記発熱量差から前記混合比率を算出するステップと、を含むバイオ燃料混合比率推定方法。

【請求項9】

前記発熱量差を算出するステップは、前記内燃機関の回転速度と燃料噴射量と前記積分値偏差とから、前記発熱量差を算出し、
前記混合比率を算出するステップは、前記燃料噴射量と前記発熱量差とから、前記混合比率を算出する、請求項８に記載のバイオ燃料混合比率推定方法。

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、バイオ燃料混合比率推定装置およびバイオ燃料混合比率推定方法に関する。

【背景技術】

【0002】

脱炭素社会に向けて、カーボンニュートラルで環境負荷の少ないバイオ燃料の使用が拡大している。たとえば、菜種油や廃食用油を原料としたバイオディーゼル燃料を軽油に混合した、バイオ燃料混合軽油がディーゼルエンジンの燃料として用いられている。

【0003】

バイオディーゼル燃料は、軽油に比べて発熱量が小さい。このため、軽油に混合されたバイオディーゼル燃料の比率（混合比率）に応じて、燃料噴射量や燃料噴射時期等の制御を行うことが望ましい。

【0004】

排気ガス中の煤等のＰＭ（Particulate Matter）を捕集するＤＰＦ（Diesel Particulate Filter）を備える場合には、ＤＰＦに堆積するＰＭの量を燃料噴射量等の運転条件から推定し、ＤＰＦの推定ＰＭ堆積量が所定量に達する毎に、ＤＰＦの再生制御を行う。ＤＰＦの再生制御は、燃料のポスト噴射を実行することにより、未燃焼燃料をＤＰＦの上流に設置した酸化触媒等で燃焼することにより行われる。バイオディーゼル燃料の混合比率が大きいほど、燃焼による煤の発生量が少なくなる。このため、軽油の使用を前提に適合を行った場合、バイオディーゼル燃料の混合比率が大きい燃料が使用されると、実際のＰＭ堆積量が所定量に達する前に推定ＰＭ堆積量が所定量に達して、余剰な再生制御が実行される。バイオディーゼル燃料は、軽油に比べて揮発し難いため、ポスト噴射に起因したオイル希釈が進行する可能性がある。このため、余剰な再生制御を抑制して、燃費の向上を図るとともにオイル希釈の進行を抑制することが望まれ、バイオディーゼル燃料の混合比率に応じて、ＰＭ堆積量を推定する際の制御定数（適合定数）を変更することが好ましい。

【0005】

特開２０１１－１４９２９７号公報（特許文献１）には、燃料圧力センサ（筒内圧センサ）を用いて検出した筒内圧から熱発生量を求め、バイオディーゼル燃料の混合比率を特定することが開示されている。

【先行技術文献】

【特許文献】

【0006】

【文献】特開２０１１－１４９２９７号公報

【発明の概要】

【発明が解決しようとする課題】

【0007】

特許文献１では、筒内圧から熱発生量を求め、バイオディーゼル燃料の混合比率を特定しており、筒内圧センサを備えない内燃機関（エンジン）においては、混合比率を特定することができない。しかし、一般的な電子制御式の内燃機関では、内燃機関の回転速度を検出するため、クランク角センサを備えている。

【0008】

本開示は、クランク角センサを用いて、バイオ燃料の混合比率を推定することを目的とする。

【課題を解決するための手段】

【0009】

本開示のバイオ燃料混合比率推定装置は、内燃機関に供給される燃料に混合されたバイオ燃料の混合比率を推定するものである。バイオ燃料混合比率推定装置は、内燃機関のクランク角を検出するクランク角センサと、クランク角に基づいてクランク角加速度を算出するクランク角加速度算出部と、クランク角加速度のうち正のクランク角加速度を積分した角加速度積分値を算出する積分値算出部と、基準燃料における角加速度積分値である基準値を記憶する基準値記憶部と、積分値算出部で算出した角加速度積分値と基準値との差である積分値偏差を算出する偏差算出部と、積分値偏差に基づいて発熱量差を算出する発熱量差算出部と、発熱量差に基づいて混合比率を算出する混合比率算出部と、を備える。

【0010】

この構成によれば、クランク角加速度算出部は、クランク角センサで検出したクランク角に基づいて、クランク角加速度を算出する。積分値算出部は、クランク角加速度のうち正のクランク角加速度を積分した角加速度積分値を算出する。基準値記憶部は、基準燃料における角加速度積分値である基準値を記憶している。偏差算出部は、積分値算出部で算出した角加速度積分値と基準値記憶部から読み出した基準値との差である積分値偏差を算出する。発熱量差算出部は、積分値偏差に基づいて発熱量差を算出する。混合比率算出部は、発熱量差に基づいて混合比率を算出する。

【0011】

内燃機関の筒内で燃焼する燃料の発熱量が小さくなると、出力トルクが低下するので、クランク角加速度のうち正のクランク角加速度を積分した角加速度積分値が小さくなる。積分値算出部で算出した角加速度積分値は、燃料の発熱量を表すパラメータと見做すことができる。積分値算出部で算出した角加速度積分値と基準値記憶部から読み出した基準値との差である積分値偏差は、内燃機関に供給された燃料と基準燃料との発熱量差に相当するパラメータである。発熱量差算出部は、偏差算出部で算出した積分値偏差から発熱量差を算出する。この発熱量差に基づいて、混合比率算出部は、混合比率を算出することができる。したがって、クランク角センサで検出したクランク角から、バイオ燃料の混合比率を推定することができる。

【0012】

好ましくは、バイオ燃料混合比率推定装置は、内燃機関の回転速度と燃料噴射量と積分値偏差とをパラメータとした発熱量差マップを備え、発熱量差算出部は、発熱量差マップを用いて発熱量差を算出するようにしてもよい。

【0013】

また、バイオ燃料混合比率推定装置は、燃料噴射量と発熱量差をパラメータとした混合比率マップを備え、混合比率算出部は、混合比率マップを用いて混合比率を算出するようにしてもよい。

【0014】

これらの構成によれば、回転速度や燃料噴射量が異なる運転状態においても、精度よく混合比率を算出することが可能になる。

【0015】

好ましくは、バイオ燃料混合比率推定装置は、回転速度の変化量が所定範囲内にあるか否かを判定する定常状態判定部を備え、定常状態判定部で変化量が所定範囲内にあると判定されたとき、積分値算出部が角加速度積分値を算出するようにしてもよい。

【0016】

この構成によれば、角加速度積分値の変動が少ない運転状態において混合比率を算出するので、混合比率の算出精度を高めることができる。

【0017】

好ましくは、バイオ燃料混合比率推定装置は、燃料タンクへの給油が行われた後、最初に内燃機関が運転されたとき、積分値算出部が角加速度積分値を算出するようにしてもよい。

【0018】

給油が行われると、次回の給油時まで、内燃機関に供給される燃料の混合比率が変化することはない。この構成によれば、燃料タンクへの給油が行われた後、最初に内燃機関が運転された際の角加速度積分値に基づいて、混合比率を算出するので、混合比率を算出するための負荷を低減することができる。

【0019】

好ましくは、内燃機関の筒内における空気過剰率が１より大きい所定値以上のとき、積分値算出部が角加速度積分値を算出するようにしてもよい。

【0020】

筒内における空気過剰率が小さいときには、筒内の燃料が完全燃焼し難くなる。特に、空気過剰率が１より小さいときには、燃料の燃焼に必要な空気が不足する。このため、空気過剰率が小さい場合には、精度よく発熱量差を求めることができない。この構成によれば、空気過剰率が１より大きい所定値以上のとき、角加速度積分値を算出するので、精度よく混合比率を算出することができる。

【0021】

好ましくは、基準燃料はバイオ燃料が混合されていない燃料であり、基準値記憶部は、バイオ燃料が混合されていない燃料における角加速度積分値である基準値を記憶していてよい。

【0022】

この構成によれば、通常の（バイオ燃料が混合されていない）燃料を用いて基準値を設定（適合）することができるので、基準値の設定が比較的容易になる。

【0023】

本開示のバイオ燃料混合比率推定方法は、内燃機関に供給される燃料に混合されたバイオ燃料の混合比率を推定する。バイオ燃料混合比率推定方法は、内燃機関のクランク角からクランク角加速度を算出するステップと、正のクランク角加速度を積分し角加速度積分値を算出するステップと、基準燃料における角加速度積分値である基準値を読み出すステップと、算出した角加速度積分値と読み出した基準値との差である積分値偏差を算出するステップと、積分値偏差を用いて発熱量差を算出するステップと、発熱量差から混合比率を算出するステップと、を含む。

【0024】

この構成によれば、内燃機関のクランク角からクランク角加速度を算出し、正のクランク角加速度を積分した角加速度積分値と基準値の偏差である積分値偏差を算出する。そして、積分値偏差を用いて算出した発熱量差から混合比率を算出する。したがって、内燃機関のクランク角から、バイオ燃料の混合比率を推定することができる。

【0025】

好ましくは、発熱量差を算出するステップは、内燃機関の回転速度と燃料噴射量と積分値偏差とから発熱量差を算出し、混合比率を算出するステップは、燃料噴射量と発熱量差とから混合比率を算出するようにしてもよい。

【0026】

この構成によれば、回転速度や燃料噴射量が異なる運転状態においても、精度よく混合比率を算出することが可能になる。

【発明の効果】

【0027】

本開示によれば、クランク角センサを用いて、バイオ燃料の混合比率を推定することができる。

【図面の簡単な説明】

【0028】

【図1】本実施の形態に係るエンジン１の全体構成図である。

【図2】エンジンＥＣＵ２００に構成されたバイオ燃料混合比率推定部１００の機能ブロックを示す図である。

【図3】クランク角加速度αの一例を示した図である。

【図4】基準値記憶部１０７に記憶されている基準値ΣαＲのマップの例を示した図である。

【図5】発熱量差マップ１０９の一例を示す図である。

【図6】混合比率マップ１１２の一例を示す図である。

【図7】エンジンＥＣＵ２００で実行される、混合比率推定制御の処理を示すフローチャートである。

【図8】エンジンＥＣＵで実行される、混合比率推定制御の割り込み要求処理を示すフローチャートである。

【図9】変形例２において、エンジンＥＣＵ２００で実行される、混合比率推定制御の処理を示すフローチャートである。

【発明を実施するための形態】

【0029】

以下、本開示の実施の形態について、図面を参照しながら詳細に説明する。なお、図中同一または相当部分には同一符号を付してその説明は繰返さない。

【0030】

図１は、本実施の形態に係るエンジン１の全体構成図である。エンジン１は、圧縮自己着火式内燃機関（ディーゼルエンジン）であり、エンジン本体１０のシリンダ（気筒）１２に形成された燃焼室に、燃料噴射弁（インジェクター）１４から燃料を噴射し、圧縮自己着火を行う内燃機関である。本実施の形態において、エンジン１は４気筒である。エンジン１の吸気通路２０には、エアクリーナ２２、インタークーラ２４、および吸気絞り弁（電子制御スロットル）２６が設けられており、エアクリーナ２２で異物が除去された新気（空気）は、ターボ過給機３０のコンプレッサ３２で過給（圧縮）され、インタークーラ２４で冷却されて、吸気マニホールド２８に供給され、吸気ポートから各燃焼室に供給される。

【0031】

燃焼室から排出される排気（排気ガス）は、排気マニホールド５０に集められ、排気通路５２を介して、外気に放出される。また、排気の一部は、ＥＧＲ（Exhaust Gas Recirculation）通路６０を介して、吸気マニホールド２８に還流される。ＥＧＲ通路６０には、ＥＧＲクーラ６２とＥＧＲ弁６４が設けられる。

【0032】

排気通路５２には、上流側から、ターボ過給機３０のタービン３４、酸化触媒７０、ＤＰＦ（Diesel Particulate Filter）７２が設けられている。ＤＰＦ７２は、排気中の微粒子（ＰＭ）を捕集し、捕集した微粒子を適宜燃焼除去することにより浄化する。また、図示しないが、ＤＰＦ７２の下流に、ＮＯｘ吸蔵還元触媒（ＮＳＲ（NOx Storage-Reduction）触媒）が設けられる。なお、ＮＯｘ吸蔵還元触媒に代えて、尿素添加弁および選択還元触媒を設けてもよい。

【0033】

燃料タンク４０には、燃料が貯留されている。燃料タンク４０内の燃料は、高圧燃料ポンプ４１によって、燃料通路４２を介してコモンレール４３に圧送される。コモンレール４３に蓄えられた高圧の燃料が、インジェクター１４から燃焼室（筒内）に噴射される。

【0034】

エンジン１は、さらに、エンジンＥＣＵ２００（Electronic Control Unit）を備える。エンジンＥＣＵ２００は、ＣＰＵ（Central Processing Unit）２２０、処理プログラム等を記憶するＲＯＭ（Read Only Memory）およびデータを一時的に記憶するＲＡＭ（Random Access Memory）等からなるメモリ２４０、各種信号を入出力するための入出力ポート（図示せず）等を含み、メモリ２４０に記憶された情報、各種センサからの情報に基づいて所定の演算処理を実行し、インジェクター１４、吸気絞り弁２６、高圧燃料ポンプ４１等を制御する。

【0035】

各種センサとしては、たとえば、クランク角センサ１２１、アクセルペダルセンサ１２２、吸入空気量センサ１２３、燃料レベルゲージ１２４、等である。クランク角センサ１２１は、エンジン１のクランク角θを検出する。アクセルペダルセンサ１２２は、ユーザによるアクセルペダル操作量（以下「アクセル開度」ともいう）ＡＰを検出する。吸入空気量センサ１２３は、吸入空気量Ｇａを検出する。燃料レベルゲージ１２４は、燃料タンク４０に貯留された燃料量を検出する。

【0036】

上記のように構成されたエンジン１では、たとえば、菜種油や廃食用油を原料としたバイオディーゼル燃料を軽油に混合した、バイオ燃料混合軽油を燃料として用いられることがある。バイオディーゼル燃料は、軽油に比べて発熱量が小さい。また、バイオディーゼル燃料は、燃焼による煤の発生量が少なくなる。このため、軽油に混合されたバイオディーゼル燃料の比率（混合比率）に応じて、燃料噴射量や燃料噴射時期等の制御、ＤＰＦのＰＭ堆積量の推定等を行うことが望ましい。本実施の形態では、クランク角センサ１２１で検出したクランク角θを用いて、バイオディーゼル燃料の混合比率を算出する。

【0037】

図２は、エンジンＥＣＵ２００に構成されたバイオ燃料混合比率推定部１００の機能ブロックを示す図である。バイオ燃料混合比率推定部１００は、本開示のバイオ燃料混合比率推定装置に相当する。

【0038】

回転速度算出部（ＮＥ算出部）１０１は、クランク角センサ１２１で検出したクランク角θに基づいて、エンジン１の回転速度（回転数（ｒｐｍ））ＮＥを算出する。燃料噴射量算出部（Ｑｆ算出部）１０２は、回転速度ＮＥとアクセル開度ＡＰに基づいて、燃料噴射量Ｑｆ（１噴射当たりの噴射量（ｍｍ^３／ｓｔ））を算出する。なお、算出された燃料噴射量Ｑｆは、インジェクター駆動部に送信され、燃料噴射量Ｑｆに相当する燃料が、インジェクター１４から噴射される。

【0039】

定常状態判定部１０３は、回転速度ＮＥの変化量ΔＮＥが所定範囲内にあるか否かを判定する。たとえば、「｜回転速度ＮＥの前回値－回転速度ＮＥの今回値｜＜閾値Ｔ１」であるとき、定常状態であると判定する。

【0040】

クランク角加速度算出部（α算出部）１０４は、クランク角センサ１２１で検出したクランク角θに基づいて、クランク角加速度αを算出する。たとえば、クランク角θからクランク角速度ωを算出し、クランク角速度ωを時間微分することにより、クランク角加速度αを算出する。

【0041】

積分値算出部（Σα算出部）１０５は、定常状態判定部１０３が定常状態であると判定しているとき、正のクランク角加速度αを積分し、角加速度積分値Σαを算出する。図３は、クランク角加速度αを示した図である。なお、本実施の形態において、エンジン１は４気筒／４サイクルエンジンであり、１８０°ＣＡ（１８０度クランク角）毎に、燃焼行程が行われる。図３において、実線は、バイオ燃料混合軽油を使用した場合におけるクランク角加速度αを示しており、積分値算出部１０５において算出される角加速度積分値Σαは、図３における斜線部の面積に相当する。本実施の形態では、１０回の燃焼行程（燃焼＋膨張行程）における平均値を、角加速度積分値Σαとして算出しているが、所定時間における燃焼行程の平均値であってもよい。

【0042】

偏差算出部（ＤΣα算出部）１０６は、角加速度積分値Σαと基準値ΣαＲとの差である積分値偏差ＤΣαを算出する。基準値ΣαＲは、基準燃料における角加速度積分値であり、基準値記憶部（ΣαＲ記憶部）１０７に記憶されている。図４は、基準値記憶部１０７に記憶されている基準値ΣαＲのマップの例を示した図である。なお、基準値記憶部１０７は、メモリ２４０内に構成される。本実施の形態において、基準燃料は、（バイオディーゼル燃料が混合されていない）軽油である。基準値ΣαＲは、基準燃料を使用してエンジン１を作動した場合における、正のクランク角加速度αを積分した値である。軽油の発熱量はバイオ燃料混合軽油より発熱量が大きいので、基準値ΣαＲは、図３に示す破線で囲まれた面積に相当し、予め実験等によって求められる。図４に示すように、本実施の形態では、基準値ΣαＲは、燃料噴射量Ｑｆと回転速度ＮＥをパラメータとしたマップとして基準値記憶部１０７に記憶されているが、燃料噴射量Ｑｆのみをパラメータとしたマップであってもよい。

【0043】

偏差算出部１０６は、燃料噴射量Ｑｆと回転速度ＮＥに基づいて、図４のマップから基準値ΣαＲを読み出す。そして、偏差算出部１０６は、角加速度積分値Σαと基準値ΣαＲの差である積分値偏差ＤΣαを算出する。本実施の形態では、基準値ΣαＲから角加速度積分値Σαを減算することにより、積分値偏差ＤΣα（＝ΣαＲ－Σα）を算出する。

【0044】

発熱量差算出部（Ｃｄ算出部）１０８は、発熱量差マップ（Ｃｄマップ）１０９を用いて発熱量差Ｃｄを算出する。図５は、発熱量差マップ１０９の一例を示す図である。発熱量差マップ１０９は、所謂、三次元マップであり、回転速度ＮＥ、燃料噴射量Ｑｆ、および、積分値偏差ＤΣαをパラメータとしたマップである。本実施の形態では、図５に示すように、燃料噴射量Ｑｆと積分値偏差ＤΣαをパラメータとした二次元マップが、回転速度ＮＥ毎に、メモリ２４０へ記憶されている。発熱量差算出部１０８は、回転速度ＮＥ、燃料噴射量Ｑｆ、および、積分値偏差ＤΣαから、発熱量差マップ１０９を用いて発熱量差Ｃｄを算出する。なお、発熱量差マップ１０９に無い点（数値）に関しては、線形補間によって発熱量差Ｃｄを算出するようにしてよい。発熱量差マップ１０９は、シミュレーションや実験等によって予め設定される。

【0045】

混合比率算出部（Ｍｒ算出部）１１１は、混合比率マップ（Ｍｒマップ）１１２を用いて混合比率Ｍｒを算出する。図６は、混合比率マップ１１２の一例を示す図である。混合比率マップ１１２は、燃料噴射量Ｑｆと発熱量差Ｃｄをパラメータとしたマップである。混合比率算出部１１１は、燃料噴射量Ｑｆおよび発熱量差Ｃｄから、混合比率マップ１１２を用いて混合比率Ｍｒを算出する。なお、混合比率は、軽油に混合したバイオディーゼル燃料の割合であり、たとえば、「％」で表すことができる。混合比率マップ１１２は、シミュレーションや実験等によって予め設定される。

【0046】

定数変更部１１３は、混合比率算出部１１１で算出した混合比率Ｍｒを用いて、燃料噴射量Ｑｆ、燃料噴射時期、ＰＭ堆積量推定制御の制御定数（適合定数）を変更する。

【0047】

図７は、エンジンＥＣＵ２００で実行される、混合比率推定制御の処理を示すフローチャートである。このフローチャートは、エンジン１の作動中に、所定期間毎に割り込み処理される。ステップ（以下、ステップを「Ｓ」と略す）１０では、定常状態判定部１０３で、エンジン１の運転状態が定常状態であるか否かを判定する。「｜回転速度ＮＥの前回値－回転速度ＮＥの今回値｜＜閾値Ｔ１」であれば、定常状態であり肯定判定され、Ｓ１１へ進む。「｜回転速度ＮＥの前回値－回転速度ＮＥの今回値｜≧閾値Ｔ１」であれば、否定判定され、今回のルーチンを終了する。

【0048】

Ｓ１１において、積分値算出部１０５は、正のクランク角加速度αを積分し、角加速度積分値Σαを算出したあと、Ｓ１２へ進む。なお、クランク角加速度αは、図示しないルーチン（たとえば、回転速度ＮＥ算出ルーチン）において、クランク角加速度算出部１０４で、クランク角センサ１２１で検出したクランク角θからクランク角速度ωを算出し、クランク角速度ωを時間微分することにより算出される。

【0049】

積分値算出部（Σα算出部）１０５は、定常状態判定部１０３が定常状態であると判定しているとき、正のクランク角加速度αを積分し、角加速度積分値Σαを算出する。

【0050】

Ｓ１２において、偏差算出部１０６は、基準値記憶部１０７（メモリ２４０）に記憶されている基準値ΣαＲを読み出し、角加速度積分値Σαと基準値ΣαＲとの差である積分値偏差ＤΣαを算出する。

【0051】

続くＳ１３で、発熱量差算出部１０８は、発熱量差マップ１０９を用いて、回転速度ＮＥと燃料噴射量Ｑｆと積分値偏差ＤΣαとから発熱量差Ｃｄを算出し、Ｓ１４へ進む。

【0052】

Ｓ１４において、混合比率算出部１１１は、燃料噴射量Ｑｆおよび発熱量差Ｃｄから、混合比率マップ１１２を用いて混合比率Ｍｒを算出し、今回のルーチンを終了する。

【0053】

エンジン１の燃焼室（筒内）で燃焼する燃料の発熱量が小さくなると、出力トルクが低下するので、クランク角加速度αのうち正のクランク角加速度αを積分した角加速度積分値Σα（図３における斜線部の面積）が小さくなる。積分値算出部１０５で算出した角加速度積分値Σαは、燃料の発熱量を表すパラメータと見做すことができる。積分値算出部１０５で算出した角加速度積分値Σαと基準値記憶部１０７から読み出した基準値ΣαＲとの差である積分値偏差ＤΣαは、エンジン１に供給された燃料と基準燃料との発熱量差に相当するパラメータである。発熱量差算出部１０８は、発熱量差マップ１０９を用いて、偏差算出部１０６で算出した積分値偏差ＤΣαから発熱量差Ｃｄを算出する。この発熱量差Ｃｄに基づいて、混合比率算出部１１１は、混合比率マップ１１２を用いて、混合比率Ｍｒを算出することができる。

【0054】

上記の実施形態によれば、クランク角センサ１２１で検出したクランクθから、バイオ燃料混合軽油の混合比率Ｍｒを推定することができる。上記の実施形態において、クランク角センサ１２１は、回転速度ＮＥを検出（算出）するために、エンジン１に設けられたものであり、新たに燃焼圧センサ（筒内圧センサ）等を設けることなく、バイオ燃料混合軽油の混合比率Ｍｒを推定することができる。

【0055】

（変形例１）
燃料タンク４０への給油が行われると、次回の給油時まで、エンジン１に供給される燃料（バイオ燃料混合軽油）の混合比率が変化することはない。したがって、燃料タンク４０に給油が行われたあと、最初にエンジン１の運転が行われた際に、図７の混合比率推定制御の処理を実行することによって、エンジンＥＣＵ２００の処理負荷（演算負荷）を低減することができる。

【0056】

上記実施の形態では、図７の混合比率推定制御の処理（フローチャート）は、エンジン１の作動中に所定期間毎に割り込み処理していた。変形例１では、混合比率推定制御は割り込み要求によって処理される。図８は、エンジンＥＣＵ２００で実行される、混合比率推定制御の割り込み要求処理を示すフローチャートである。このフローチャートは、エンジン１が始動されたときに実行される。たとえば、イグニッションスイッチ（始動スイッチ）がＯＮされ、エンジン１が始動されると、Ｓ２０において、燃料タンク４０への給油が行われたか否かを判定する。本実施の形態では、燃料レベルゲージ１２４で検出した、燃料タンク４０の燃料量を用いて、給油が行われたか否かを判定する。たとえば、前回、エンジン１が停止したとき（イグニッションスイッチがＯＦＦされたとき）の燃料量よりも、今回のエンジン１始動時の燃料量の方が多いとき、給油が行われたと判定する。なお、エンジン１の停止から今回のエンジン１の始動までの間に、燃料リッド（フェールリッド）が開かれた履歴がある場合、給油が行われたと判定してよい。Ｓ２０において、燃料タンク４０への給油が行われたと判定されるとＳ２１へ進む。燃料タンク４０への給油が行われていないと判定されると、今回のルーチンを終了する。

【0057】

Ｓ２１では、エンジン１の始動後、所定期間が経過したか否かを判定する。所定期間は、給給油前に燃料タンク４０に貯留されていた燃料と給油されたバイオ燃料混合軽油が、混ざり合うのに十分な時間とされる。エンジン１が車両に搭載されている場合、走行時の振動によって、燃料タンク４０内の燃料の攪拌が期待できるので、所定期間は、エンジン１の始動後の数分間であってよい。エンジン１の始動後、所定期間が経過するまで、Ｓ２１が繰り返し処理され、エンジン１の始動後、所定期間が経過すると、Ｓ２１で肯定判定されＳ２２へ進む。

【0058】

Ｓ２２では、混合比率推定制御の割り込み要求を行い、今回のルーチンを終了する。この変形例１では、Ｓ２２の割り込み要求が行われると、図７の混合比率推定制御の処理が実行される。

【0059】

この変形例１によれば、燃料タンク４０への給油が行われた後、最初にエンジン１が運転された際の角加速度積分値Σαに基づいて、混合比率Ｍｒを算出するので、混合比率Ｍｒを算出するための演算負荷を低減することができる。

【0060】

（変形例２）
エンジン１の燃焼室（筒内）における空気過剰率が小さいときには、燃焼室の燃料が完全燃焼し難くなる。特に、空気過剰率が１より小さいときには、燃料の燃焼に必要な空気が不足する。このため、空気過剰率が小さい場合には、バイオ燃料混合軽油と基準燃料の発熱量差Ｃｄを精度よく求めることができない。

【0061】

たとえば、ＮＳＲ触媒に吸蔵されたＮＯｘを還元浄化するとき、吸気絞り弁２６を閉じるとともに燃料噴射量Ｑｆを増量することにより、吸気過剰率を１より小さくして、還元雰囲気を生成する。この場合、燃焼室に供給された燃料の一部は、ＮＯｘ還元用の未燃燃料として排出されるため、角加速度積分値Σαを元いて発熱量差Ｃｄを精度よく求めることができない。また、空気過剰率が１より大きい場合であっても、空気過剰率が小さい場合には、完全燃焼し難く、角加速度積分値Σαを元いて発熱量差Ｃｄを精度よく求めることができない。

【0062】

変形例２では、燃料が良好に完全燃焼する空気過剰率以上のとき、角加速度積分値Σαを用いて混合比率Ｍｒを算出することにより、精度よく混合比率Ｍｒを算出する。図９は、変形例２において、エンジンＥＣＵ２００で実行される、混合比率推定制御の処理を示すフローチャートである。このフローチャートは、エンジン１の作動中に、所定期間毎に割り込み処理される。このフローチャートは、上記の実施形態における図７のフローチャートに対して、Ｓ４０を追加したものであり、Ｓ１０～Ｓ１４は図７のフローチャートと同様であるので、その説明を省略する。

【0063】

定常状態（｜回転速度ＮＥの前回値－回転速度ＮＥの今回値｜＜閾値Ｔ１）であり、Ｓ１０で肯定判定されると、Ｓ４０では、空気過剰率λが所定値Ａより大きいか否かを判定する。空気過剰率λは、たとえば、吸入空気量センサ１２３で検出した吸入空気量Ｇａと燃料噴射量Ｑｆから求める。また、排気通路５２にλセンサ（全域空燃比センサ）を設けて、空気過剰率λを検出するようにしてもよい。空気過剰率λが所定値Ａ（たとえば、１．２）より小さいときには、否定判定され今回のルーチンを終了する。空気過剰率λが所定値Ａより大きいときには、肯定判定されＳ１１へ進む。

【0064】

この変形例２によれば、空気過剰率が１より大きい所定値以上のとき、角加速度積分値Σαを算出し混合比率Ｍｒを算出するので、精度よく混合比率Ｍｒを推定することができる。

【0065】

上記の実施の形態では、基準燃料として軽油を用いていた。しかし、基準燃料として、所定の混合比率を有するバイオディーゼル燃料混合軽油を、基準燃料としてもよい。たとえば、５％のバイオディーゼル燃料を混合した、所謂、Ｂ５と表記されたバイオディーゼル燃料混合軽油を基準燃料とし、基準値ΣαＲ（図４）、発熱量差マップ（図５）、混合比率マップ（図６）を設定し、混合比率を推定するようにしてもよい。

【0066】

上記の実施の形態では、エンジン１としてディーゼルエンジンを説明したが、エンジン１はガソリンエンジン（火花点火式エンジン）であってもよい。ガソリンエンジンでは、ガソリンに、バイオエタノール、あるいは、バイオＥＴＢＥ（ethyl tertiary-butyl ether）を混合した、バイオ燃料混合ガソリンが用いられる。バイオエタノール、バイオＥＴＢＥの発熱量は、ガソリンに比較して小さく、上記の実施の形態と同様に、バイオ燃料混合ガソリンにおいても、バイオ燃料（バイオエタノール、バイオＥＴＢＥ）の混合比率を推定することができる。

【0067】

今回開示された実施の形態は、すべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した実施の形態の説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味及び範囲内でのすべての変更が含まれることが意図される。

【符号の説明】

【0068】

１ エンジン、１０ エンジン本体、１２ シリンダ（気筒）、１４ インジェクター、２０ 吸気通路、２２ エアクリーナ、２４ インタークーラ、２６ 吸気絞り弁、２８ 給気マニホールド、３０ ターボ過給機、３２ コンプレッサ、３４ タービン、４０ 燃料タンク、４１ 高圧燃料ポンプ、４２ 燃料通路、４３ コモンレール、５０ 排気マニホールド、５２ 排気通路、６０ ＥＲＧ通路、６２ ＥＧＲクーラ、６４ ＥＧＲ弁、７０ 酸化触媒、７２ ＤＰＦ、１００ バイオ燃料混合比率推定部、１０１ 回転速度算出部、１０２ 燃料噴射量算出部、１０３ 定常状態判定部、 １０４ クランク角加速度算出部、１０５ 積分値算出部、１０６ 偏差算出部、１０７ 基準値記憶部、１０８ 発熱量差算出部、１０９ 発熱量差マップ、１１１ 混合比率算出部、１１２ 混合比率マップ、１２１ クランク角センサ、１２２ アクセルベダルセンサ、１２３ 吸入空気量センサ、１２４ 燃料レベルゲージ、２００ エンジンＥＣＵ、２１０ ＣＰＵ、２２０ メモリ。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】