脱硝触媒の劣化再生予測プログラム及びそれを用いた脱硝システム
【課題】脱硝触媒の劣化及び再生を総括する反応モデルを用いて、脱硝触媒の劣化再生予測プログラムを提供する。
【解決手段】コンピュータによって、アンモニア系還元剤のパラメータ設定を行い、脱硝触媒に関する条件設定を行い、脱硝触媒でのアンモニアの反応量及び窒素酸化物の反応量を算出し、脱硝触媒での三酸化硫黄の生成量を算出し、脱硝触媒での酸性硫安の増減量を算出し、脱硝触媒の性能劣化量を算出し、算出結果を出力する脱硝触媒の劣化再生予測プログラムとする。
【解決手段】コンピュータによって、アンモニア系還元剤のパラメータ設定を行い、脱硝触媒に関する条件設定を行い、脱硝触媒でのアンモニアの反応量及び窒素酸化物の反応量を算出し、脱硝触媒での三酸化硫黄の生成量を算出し、脱硝触媒での酸性硫安の増減量を算出し、脱硝触媒の性能劣化量を算出し、算出結果を出力する脱硝触媒の劣化再生予測プログラムとする。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、脱硝触媒の劣化再生予測プログラム及びそれを用いた脱硝システムに関する。
【背景技術】
【０００２】
国際海事機関（ＩＭＯ）によるＮＯｘ規制は現在の一次規制から強化されることが決定している。２０１１年には二次規制、２０１６年には三次規制となる予定だが、三次規制になるとエンジン本体だけではなく、後処理装置の導入によるＮＯｘ低減も必要になると考えられている。尿素水を用いた選択触媒還元（ＳＣＲ：Selective Catalytic Reduction）はＮＯｘ低減のための有力な後処理技術の一つである。ＳＣＲのメリットとしては、浄化率が高いこと、ＳＯｘを含んだ排ガス中でも利用ができること等が挙げられる。そのため、舶用ディーゼルエンジンへの適用が検討されているが、触媒性能の劣化やアンモニア（ＮＨ３）スリップ等の問題点もある。
【０００３】
三次規制は指定された海域内でのみ適用され、それ以外の海域では二次規制が適用されるが、この海域は基本的に港湾や輻輳海域等の船舶交通量の比較的大きく、陸地に近い場所に設定されると考えられる。その時、ＳＣＲはこの海域を航行する限られた時間のみで運用され、航行中の大部分はＳＣＲを作動させないことが想定できる。そうなると、ＳＣＲの運用に関する自由度は上がり、例えば劣化した触媒の再生をＮＯｘ還元の必要のない海域航行中に行うようなことも可能になるためＳＣＲの寿命や劣化再生方法について調べることが重要となってくる。しかし、これらの検討の手段として実験を行おうとすると非常に長い時間がかかることや温度や排ガス成分等のパラメータが非常に多いことから、数値解析による予測が有効な手段になると考えられる。
【０００４】
ＳＣＲを用いた脱硝反応に関して基礎的な実験・解析から数値解析まで，数多くの研究が行われている。
【０００５】
例えば、窒素酸化物を含有する排ガス中に還元剤であるアンモニアを投入する際の制御について開示されている（特許文献１，２等）。特許文献１では、触媒ユニット、エンジンの作動に関するパラメータを用いて制御が行われる。特許文献２では、窒素酸化物（ＮＯｘ）浄化特性の多項式を用いてアンモニア吸着量の目標値を設定する方法が採用されている。
【０００６】
また、脱硝システムにおける脱硝触媒の劣化を診断する診断装置が開示されている（特許文献３）。ここでは、アンモニアが漏れるまで脱硝触媒へ供給された還元剤の総量に基づいて脱硝触媒が劣化したか否かを判定している。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特表平０８−５０９７９５号公報
【特許文献２】特開２００８−１９６３４０号公報
【特許文献３】特開２００９−１２７４９６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかし、脱硝触媒の劣化を含めた計算を行っている例は少なく、特に、舶用ディーゼルエンジンを想定した比較的低排気温度で、高硫黄含有燃料を用いた場合を対象にして検討している例はない。さらに、脱硝触媒の劣化及び再生を総括する反応モデルを用いて脱硝触媒の劣化再生予測を行う技術はない。
【０００９】
本発明は、上記の課題を鑑み、脱硝触媒の劣化及び再生を総括する反応モデルを用いて、脱硝触媒の劣化再生予測プログラム及びそれを用いた脱硝システムを提供することを目的とする。
【課題を解決するための手段】
【００１０】
請求項１に対応した脱硝触媒の劣化再生予測プログラムは、コンピュータによって、燃焼機器からの排ガスにアンモニア系還元剤を添加して脱硝触媒で窒素酸化物の分解を行う脱硝システムにおける前記脱硝触媒の劣化と再生を予測するプログラムであって、前記アンモニア系還元剤のパラメータ設定を行うステップ１と、前記脱硝触媒に関する条件設定を行うステップ２と、前記脱硝触媒でのアンモニアの反応量及び窒素酸化物の反応量を算出するステップ３と、前記脱硝触媒での三酸化硫黄の生成量を算出するステップ４と、前記脱硝触媒での酸性硫安の増減量を算出するステップ５と、前記脱硝触媒の性能劣化量を算出するステップ６と、算出結果を出力するステップ７と、を備えることを特徴とする。
【００１１】
ここで、前記アンモニア系還元剤は、尿素水が好適であり、前記ステップ１における前記アンモニア系還元剤のパラメータは、前記尿素水の分解したアンモニアの濃度を含むことが好適である。これにより、還元剤がアンモニアである場合の脱硝システムにおける脱硝触媒の劣化又は再生を適切に予測することが可能となる。なお、アンモニア系還元剤としては、尿素水や尿素以外に、アンモニアや炭酸アンモニウム、その他アンモニア前駆体等を使用することができる。
【００１２】
また、前記ステップ２における条件は、前記脱硝触媒の温度及び体積と、前記排ガスの成分、流量及び温度と、であることが好適である。これにより、脱硝触媒のパラメータに応じて脱硝システムにおける脱硝触媒の劣化又は再生を適切に予測することが可能となる。
【００１３】
また、前記ステップ３は、Ｒｉｄｅａｌ−Ｅｌｅｙ機構を用いた脱硝反応速度式に基づいたアンモニアと窒素酸化物との反応量を算出する脱硝モデルを用いることが好適である。Ｒｉｄｅａｌ−Ｅｌｅｙ機構を用いた脱硝反応速度式を適用することにより、脱硝システムにおける脱硝触媒の劣化又は再生の予測の確度を高めることができる。
【００１４】
また、前記ステップ４は、非定常かつ、排ガス成分の影響を考慮した速度式に基づいた三酸化硫黄の生成量を算出する二酸化硫黄酸化反応モデルを用いることが好適である。これにより、三酸化硫黄の生成量に応じた脱硝システムにおける脱硝触媒の劣化又は再生を適切に予測することが可能となる。
【００１５】
また、前記ステップ５は、アンモニアと三酸化硫黄と水蒸気より酸性硫安を生じる平衡式の平衡定数を雰囲気ガス温度並びにアンモニア、三酸化硫黄及び水蒸気の分圧に基づいて定めた酸性硫安の増減量を算出する酸性硫安析出分解モデルを用いることを特徴とする。これにより、酸性硫安の増減量を考慮して脱硝システムにおける脱硝触媒の劣化又は再生を適切に予測することが可能となる。
【００１６】
また、前記酸性硫安析出分解モデルは、ケルビン方程式で補正されたものであることが好適である。これにより、酸性硫安の析出及び分解を正確に予測することができ、脱硝システムにおける脱硝触媒の劣化又は再生の予測の確度を高めることが可能となる。
【００１７】
また、前記ステップ６は、酸性硫安の増減量と劣化関数との関係を実験結果から求め、以下の式で算出される反応速度定数に基づいて行われることが好適である。これにより、脱硝触媒の劣化関数を考慮して、脱硝システムにおける脱硝触媒の劣化又は再生の予測の確度を高めることが可能となる。
ｋｃ＝ＡＮＯｅｘｐ（−ＥＮＯ／ＲＴｃａｔ）×Φ
ｋｃ：反応速度定数［ｍ３Ｎ／ｍ２ｈ］
ＡＮＯ、ＥＮＯ：ＮＯの反応に関する頻度因子［ｍ３Ｎ／ｍ２ｈ］
Ｒ：一般ガス定数［Ｊ／ｍｏｌＫ］
Ｔｃａｔ：触媒表面温度［Ｋ］
Φ：劣化関数
【００１８】
また、前記脱硝触媒の出口におけるアンモニア量と窒素酸化物量を算出するステップ８をさらに備えることが好適である。これにより、脱硝触媒の劣化又は再生を考慮して、脱硝処理におけるアンモニア量と窒素酸化物量の変化を把握することができる。
【００１９】
また、前記脱硝触媒の温度が所定の温度に達したか否かを判定するステップ９をさらに備えることが好適である。これにより、脱硝触媒の温度に応じた脱硝触媒の劣化又は再生をより正確に予測することが可能となる。
【００２０】
また、前記所定の温度を更新するステップ１０をさらに備えることが好適である。これにより、脱硝触媒の複数温度における脱硝触媒の劣化又は再生をより正確に予測することが可能となる。
【００２１】
また、前記脱硝システムの運転時間が所定の時間に達したか否かを判断するステップ１１をさらに備えることが好適である。これにより、脱硝システムの運転時間を考慮した脱硝触媒の劣化又は再生の予測を行うことが可能となる。
【００２２】
また、前記ステップ７は、時間と触媒性能劣化量の関係、時間と前記脱硝触媒の出口におけるアンモニア量の関係の少なくとも一方を出力することが好適である。これにより、脱硝システムの運転時間に応じた脱硝触媒の劣化又は再生を考慮して、触媒の性能劣化量、脱硝触媒の出口におけるアンモニア量を把握することができる。
【００２３】
請求項１４に対応した脱硝システムは、脱硝触媒の劣化再生予測プログラムを用いて上記脱硝触媒の劣化及び再生の少なくとも一方の予測を行うことが好適である。これにより、脱硝触媒の劣化及び再生の少なくとも一方の予測が可能な脱硝システムを実現することができる。
【００２４】
また、前記脱硝触媒の劣化再生予測プログラムの出力に基づいて、前記脱硝システムの操作量を変更し、前記ステップ１、２及び１０の少なくとも１つの設定を変更することが好適である。これにより、実際の脱硝システムにおける操作量に応じた脱硝触媒の劣化又は再生を考慮し、状況に応じた脱硝システムの制御を行うことが可能となる。
【発明の効果】
【００２５】
本発明の脱硝触媒の劣化再生予測プログラムによれば、アンモニア系還元剤を添加して脱硝触媒で窒素酸化物の分解を行う際の脱硝反応をより正確に予測することができる。また、アンモニア系還元剤のパラメータや脱硝触媒に関する条件を変えた場合の脱硝触媒の性能劣化量の算出が容易にできる。
これにより、例えば、脱硝触媒の出口におけるアンモニア濃度、窒素酸化物の算出や劣化に至る運転時間の算出が容易となる。また、劣化した脱硝触媒の触媒温度を更新した場合の、脱硝触媒の再生の予測を正確に行うことができる。
また、本発明の脱硝触媒の劣化再生予測プログラムを用いた脱硝システムによれば、脱硝触媒の劣化及び再生の少なくとも一方の予測ができ、脱硝システムの的確な運転に寄与することができる。例えば、脱硝触媒の劣化レベルの表示や、運転条件を変更した場合の窒素酸化物やアンモニアの排出状態の予測、脱硝触媒の交換や再生のタイミングの告知、再生温度と再生時間の予測等が容易となる。
また、劣化再生プログラムの出力に基づいて脱硝システムの操作量を変更した場合は、状況に応じ脱硝システムを自動的に運転することも可能となる。
【図面の簡単な説明】
【００２６】
【図１】本発明の実施の形態における脱硝触媒の劣化再生予測システムの構成を示す図である。
【図２】本発明の実施の形態における脱硝触媒の劣化再生予測処理のフローチャートを示す図である。
【図３】脱硝反応を説明する図である。
【図４】酸性硫安付着量と劣化関数Φとの関係を示す図である。
【図５】濃度比αに対する脱硝率と触媒温度との関係を示す図である。
【図６】ＳＶ，温度変更時の脱硝率と濃度比αとの関係を示す図である。
【図７】脱硝率の時間変化の実験値及び計算値との比較を示す図である。
【図８】硫黄成分と酸性硫安析出上限温度との関係を示す図である。
【図９】空気過剰率による脱硝率の時間変化を示す図である。
【図１０】ＮＯ，ＳＯ２濃度が脱硝率に与える影響を示す図である。
【図１１】脱硝触媒の出口ＮＯｘ濃度に及ぼす濃度比αの効果を示す図である。
【図１２】ＮＨ３スリップ濃度の時間変化を示す図である。
【発明を実施するための形態】
【００２７】
＜脱硝触媒の劣化再生予測処理＞
本発明の実施の形態における脱硝触媒の劣化再生予測処理は、コンピュータシステム１００を用いて実行される。コンピュータシステム１００は、図１に示すように、処理部１０、記憶部１２、入力部１４及び出力部１６を含んで構成される。
【００２８】
処理部１０は、ＣＰＵ等で構成され、記憶部１２に予め記憶された脱硝触媒の劣化再生予測プログラムを実行することにより脱硝触媒の劣化再生予測処理の演算を行う。記憶部１２は、脱硝触媒の劣化再生予測プログラム及びそれに関する各種パラメータを格納及び保持する。記憶部１２は、半導体メモリ、ハードディスク等の記憶装置を含んで構成され、処理部１０から読み出し及び書き込み可能とされる。入力部１４は、脱硝触媒の劣化再生予測処理に必要なデータ等を入力するために用いられる。入力部１４は、キーボード、ポインティングデバイス等の入力装置を含んで構成される。出力部１６は、入力部１４から入力されたデータの表示、脱硝触媒の劣化再生予測処理で得られる中間処理結果や最終処理結果の表示等を行う。出力部１６は、ディスプレイ、プリンタ等の出力装置を含んで構成される。
【００２９】
以下、本発明の実施の形態における脱硝触媒の劣化再生予測処理について説明する。脱硝触媒の劣化再生予測処理は、処理部１０によって記憶部１２に記憶されている脱硝触媒の劣化再生予測プログラムを実行することによって、図２に示すフローチャートに沿って行われる。
【００３０】
本実施の形態では、脱硝触媒内の反応のみを考慮した一次元の簡易的な反応律速モデルを用いている。すなわち、脱硝触媒中の各物質は流れ直角方向に十分拡散し、流れ方向のみに変化することを仮定している。また、一般的に脱硝触媒には還元剤として尿素水が導入されるが、尿素水は加水分解後、アンモニア（ＮＨ３）ガスとして触媒に導入されるとした。また窒素酸化物（ＮＯｘ）については一酸化窒素（ＮＯ）のみを仮定した。考慮する物質は、アンモニア（ＮＨ３）、一酸化窒素（ＮＯ）、二酸化硫黄（ＳＯ２）、三酸化硫黄（ＳＯ３）、水（Ｈ２Ｏ）、窒素（Ｎ２）、酸素（Ｏ２）及び硫酸水素アンモニウム（ＮＨ４ＨＳＯ４：以下、酸性硫安と称する）の８種類とする。
【００３１】
本実施の形態における脱硝触媒の劣化再生予測処理では、各物質の流量は脱硝触媒を流れ方向に分割した各ブロックに存在できる時間毎に分割され、そのブロックに各物質の存在する流量（体積）に対して触媒反応の計算を行う。すなわち、ある時刻ステップ（ｔ）における各ブロックにおける脱硝触媒による脱硝反応、二硫化硫黄の酸化反応、酸性硫安の生成分解反応及び酸性硫安による脱硝触媒の劣化・再生反応の計算を行う。算出された各物質は、酸性硫安以外はそのブロックに何も残さず、次の時刻ステップ（ｔ＋１）には次のブロックに輸送されるものとする。そして、算出された各物質の濃度等を用いて、次の時刻ステップ（ｔ＋１）における各ブロックにおける触媒反応の計算を行う。このような処理を繰り返し、脱硝触媒の劣化再生予測処理を実行する。
【００３２】
なお、各ブロックでの滞留時間は、雰囲気ガス流量［ｍ３Ｎ／ｈ］を触媒体積［ｍ３］で除した空塔速度（ＳＶ：Space Velocity）の逆数として定義した。ただし、二次元の反応律速モデルに拡張することも可能である。
【００３３】
ステップＳ１０では、脱硝触媒の劣化再生予測処理で用いられるパラメータの初期設定を行う。ここでのパラメータの設定には、還元剤及び窒素酸化物に関するパラメータの設定が含まれる。上記のように還元剤として尿素水を加水分解したアンモニア（ＮＨ３）ガスが使用されるものとして、還元剤に関するパラメータは、反応速度定数ＫＮＨ３［ｍ３Ｎ／ｍ２ｈ］、アンモニアの体積濃度ＣＮＨ３［ｍ３／ｍ３］、アンモニア吸着平衡定数ＫＮＨ３等の後述の反応モデルで使用されるパラメータである。また、窒素酸化物に関するパラメータは、反応速度定数ＫＮＯ［ｍ３Ｎ／ｍ２ｈ］、窒素酸化物の体積濃度ＣＮＯ［ｍ３／ｍ３］、窒素酸化物の反応に関する頻度因子ＡＮＯ［ｍ３Ｎ／ｍ２ｈ］及び活性化エネルギＥＮＯ［Ｊ／ｍｏｌ］等の後述の反応モデルで使用されるパラメータである。また、二酸化硫黄に関する反応速度定数ｋＳ、酸素に関する反応次数αＯ２、水に関する反応次数αＷ、一酸化窒素及びアンモニアに対する係数ｋＮＯ，ｋＮＨ３、二硫化硫黄に関する頻度因子ＡＳＯ２，活性化エネルギＥＳＯ２が設定される。このようにアンモニア系還元剤のパラメータや脱硝触媒に関する条件を設定することによって、これらの設定を変更した場合の脱硝触媒の性能劣化量の算出が容易にできる。なお、使用するアンモニア系還元剤の種類により、使用するパラメータは異なる。また、尿素水にしても、尿素の含有形態や供給方法、パラメータの表現方法等により、各種のパラメータやその組み合わせがあり得る。
【００３４】
例えば、反応速度定数ＫＮＨ３は後述する実験結果と一致するように０．０５×１０６に設定することが好適である。一酸化窒素（ＮＯ）とアンモニア（ＮＨ３）の濃度比α（ＣＮＨ３／ＣＮＯ）＝１のときの実験結果より、頻度因子ＡＮＯ及び活性化エネルギＥＮＯは、それぞれ８８９８ｍ３Ｎ／ｍ２ｈ及び２３．８７×１０３Ｊ／ｍｏｌとすることが好適である。また、窒素酸化物に関するパラメータは、反応次数αＯ２は文献(E.Tronconi,et.al.,Ind.Eng.Chem.Res.38(1999)2593)から設定することができ、その他は実験結果等から表１のように設定することができる。
【表１】
【００３５】
ステップＳ１２では、脱硝触媒に関するパラメータの設定を行う。ここでは、脱硝触媒の比表面積比Ａｐ［ｍ３／ｍ３］、脱硝触媒の表面温度Ｔｃａｔ［Ｋ］、脱硝触媒の予測処理を行う温度範囲等が設定される。脱硝触媒の表面温度Ｔｃａｔは初期設定であり、脱硝触媒の予測処理を行う温度範囲や脱硝触媒の再生処理温度に応じてステップ２６において変更される。なお、脱硝触媒に関するパラメータの設定は、使用する脱硝触媒により、温度範囲や初期設定される表面温度等は異なる。また、空塔速度や周囲温度等を設定してもよい。
【００３６】
ステップＳ１４では、脱硝触媒でのアンモニアの反応量及び窒素酸化物の反応量を算出する。Ｖ２Ｏ５−ＴｉＯ２系触媒における尿素水加水分解後のアンモニアによる代表的な脱硝反応の化学反応式は、化学式１のように表される。
【化１】
【００３７】
その他に二酸化窒素（ＮＯ２）と一酸化窒素（ＮＯ）が等モルで、化学式１よりも低温、高速で起こる反応もあるが、ディーゼルエンジン排ガス中の窒素酸化物（ＮＯｘ）のほとんどは一酸化窒素（ＮＯ）であるので、酸化触媒による二酸化窒素（ＮＯ２）の生成を考慮せず、ここでは化学式１の反応のみを考える。この反応が触媒表面で起こっているとすると、図３の概念図に示すように、脱硝反応はアンモニア（ＮＨ３）の触媒活性点への吸着、吸着したアンモニア（ＮＨ３）と雰囲気に拡散している一酸化窒素（ＮＯ）との反応、生成物の活性点からの脱離と進行すると考えられる。
【００３８】
脱硝反応速度式は、アンモニア（ＮＨ３）の影響を考慮する必要があるため、脱硝機構をモデル化する手法としてＲｉｄｅａｌ−Ｅｌｅｙ機構を用いた式（１）及び（２）により記述される。
【数１】
【数２】
【００３９】
この反応モデルにより、脱硝触媒でのアンモニアの反応量及び窒素酸化物の反応量を算出する。すなわち、式（１）及び（２）から窒素酸化物（ＮＯ）の体積濃度ＣＮＯの時間的な変化を算出することができる。化学式１から窒素酸化物（ＮＯ）の体積濃度ＣＮＯの時間的な変化に対応したアンモニア（ＮＨ３）の体積濃度ＣＮＨ３の時間的な変化も算出することができる。
【００４０】
ここで、Φは脱硝触媒の劣化関数である。劣化関数Φは、脱硝触媒の劣化後の反応速度定数ｋｃと初期反応速度定数ｋｃ０の比（ｋｃ／ｋｃ０）と定義する。劣化関数Φについては後述する。劣化関数Φの初期値は予め設定しておけばよく、例えば、脱硝触媒がまったく劣化していない状態では劣化関数Φの値は１とすればよい。なお、アンモニアの反応量及び窒素酸化物の反応量の算出は、上記したＲｉｄｅａｌ−Ｅｌｅｙ機構を用いた式（１）及び（２）以外の簡略式や厳密式、また、パラメータの設定や表現方法を変えた各種の式に基づくことができる。
【００４１】
ステップＳ１６では、脱硝触媒での三酸化硫黄の生成量を算出する。燃料中に硫黄が存在する場合には、エンジン内の燃焼によって硫黄が酸化され二酸化硫黄（ＳＯ２）になり、更に二酸化硫黄（ＳＯ２）が触媒上で酸化され三酸化硫黄（ＳＯ３）になる。二酸化硫黄（ＳＯ２）の酸化反応は、化学式２で表される。
【化２】
【００４２】
ここで、二酸化硫黄（ＳＯ２）の酸化反応は、Ｖ２Ｏ５−ＴｉＯ２系触媒において窒素酸化物の脱硝反応と同時に起こり、脱硝反応に比べて非常に遅いことが知られている。この二酸化硫黄（ＳＯ２）から三酸化硫黄（ＳＯ３）への転換率は数％程度である。しかしながら、酸化した三酸化硫黄（ＳＯ３）は後述するようにアンモニア（ＮＨ３）及び水（Ｈ２Ｏ）との反応により酸性硫安になると触媒表面に付着して活性吸着点を覆うため、脱硝反応が阻害される。そのため、劣化を予測するためには二酸化硫黄（ＳＯ２）の酸化反応を予測する。
【００４３】
ここでは、二酸化硫黄（ＳＯ２）の反応モデルとして非定常での計算ができる速度式として式（３）及び（４）を用いる。なお、三酸化硫黄の生成量の算出は式（３），（４）に限定されるものでなく、設定パラメータや表現方法を変えた式、簡略式や厳密式等、各種の式に基づくことができる。
【数３】
【数４】
【００４４】
ステップＳ１８では、脱硝触媒での酸性硫安の増減量を算出する。ここで、舶用ディーゼルエンジンのような硫黄を含む燃料を用いるシステムの場合、ＳＣＲによる脱硝を行おうとすると、未反応のアンモニア（ＮＨ３）と排ガス中の三酸化硫黄（ＳＯ３）及び水（Ｈ２Ｏ）の反応によって生成される化合物が触媒上に付着し、触媒の活性を低下させることが問題となる。この触媒の劣化の有無及び進行速度は、温度によっても変わってくる。ガス中にアンモニア（ＮＨ３）、三酸化硫黄（ＳＯ３）及び水（Ｈ２Ｏ）が共存している場合に生成される化合物は、硫酸アンモニウム（（ＮＨ４）２ＳＯ４）、酸性硫安等が考えられるが、ここではその化合物を酸性硫安のみであるとして演算する。
【化３】
【数５】
【数６】
【００４５】
ここで、Ｋは平衡定数、ｐはアンモニア（ＮＨ３）、三酸化硫黄（ＳＯ３）及び水（Ｈ２Ｏ）の分圧、Ｔａは雰囲気ガス温度を示している。式（５）で示される分圧が式（６）の温度による平衡定数より大きい場合に析出が生じる。ここで、表面張力と密度は文献等から得られる値とし、細孔半径は実測値をＲｏｓｉｎ−Ｒａｍｍｌｅｒの分布関数で近似した値とすることが好適である。
【００４６】
しかし、この式は気体中での反応を示している。それを無数の細孔のある触媒中で適用するため、式（５）をケルビン方程式によって修正する。この式は、細孔内では気体中の飽和蒸気圧ＰＯよりも低い圧力Ｐによって凝縮が起こることを示す。なお、式（６）の平衡定数Ｋが式（７）の飽和蒸気圧ＰＯに相当する。
【数７】
ここで、σは表面張力、Ｍａｂ及びρａｂは酸性硫安の分子量及び密度、ｒｐｏｒｅは細孔半径である。
【００４７】
ここでは、二酸化硫黄（ＳＯ２）の反応モデルを用いて得られた二酸化硫黄（ＳＯ２）の体積濃度ＣＳＯ２の時間変化から三酸化硫黄（ＳＯ３）の体積濃度ＣＳＯ３の時間変化を求め、その算出結果を式（５）〜（７）に導入することにより酸性硫安の増加・減少量を求める。なお、酸性硫安の増減量の算出は、式（５）〜（７）に限定されるものではなく、設定パラメータや表現方法を変えた式、簡略式や厳密式等、各種の式に基づくことができる。
【００４８】
ステップＳ２０では、脱硝触媒の性能劣化再生量を算出する。脱硝触媒の劣化・再生は、脱硝触媒上に酸性硫安が付着することによって生ずる。すなわち、脱硝触媒の性能劣化再生量は、酸性硫安付着による触媒性能の劣化関数Φの値を式（２）の一酸化窒素（ＮＯ）の反応速度定数に乗算し、窒素酸化物（ＮＯ）の濃度の時間変化の減少量（増加量）として算出される。
【００４９】
ここで、劣化関数Φの値は、酸性硫安の付着量の関係から求めることができる。ここでは、図４に示す実験結果から劣化関数Φの値を求めるものとする。すなわち、横軸の触媒１ｇあたりの硫酸イオン付着量［ｍｍｏｌ／ｇｃａｔ］を酸性硫安付着量と置き換え、それに対応する劣化関数Φの値を決定する。
【００５０】
具体的には、図４中の実線で示すように酸性硫安が０．５ｍｍｏｌ／ｇｃａｔ以下の時には５次の近似式で表し、０．５ｍｍｏｌ／ｇｃａｔ以上の時には１次の近似式で表すことが好適である。なお、酸性硫安は二酸化硫黄（ＳＯ２）の酸化反応に影響を与えると考えられるが、付着量と劣化関数Φの関係が明らかでないため、ここでは影響を考慮しないものとする。なお、脱硝触媒の性能劣化再生量の算出に当たっては、他の実験結果に基づいた劣化関数や他の近似式を用いてもよい。また、脱硝触媒の性能劣化量のみを求めることも可能である。
【００５１】
ステップＳ２２では、脱硝触媒の出口でのアンモニア流出量（体積濃度ＣＮＨ３）及び一酸化窒素（ＮＯ）流出量（体積濃度ＣＮＯ）の算出が行われる。すなわち、ステップＳ２０で求められた劣化関数Φを式（２）に導入し、式（１）及び（２）の関係式から脱硝触媒の出口でのアンモニア流出量及び一酸化窒素（ＮＯ）流出量を算出する。
【００５２】
ステップＳ２４では、脱硝触媒の温度を変化させるか否かを判定する。本実施の形態では、脱硝触媒の温度に応じた脱硝触媒の劣化・再生予測を行うために、設定された脱硝触媒の温度範囲での予測処理が終了したか否かを判定する。脱硝触媒の温度範囲及び変更ステップ等は、エンジンの運転条件の変更や再生処理の条件に応じて予め設定しておけばよい。
【００５３】
脱硝触媒の劣化・再生予測を行う脱硝触媒の温度が存在すれば、ステップＳ２６において次の脱硝触媒の温度に設定し、ステップＳ１４に処理を戻し、新たに設定された脱硝触媒の温度での演算を行う。例えば、脱硝触媒の温度を初期値から更新ステップ温度毎に増加させ、所定の演算最高温度に到達するまで脱硝触媒の温度を更新する。一方、必要な脱硝触媒の温度のすべてについて演算が終了していれば、ステップＳ２８に処理を移行させる。なお、温度を上げて脱硝触媒の再生を行う時にも、脱硝触媒の温度を更新する。再生の判断は、運転時間や脱硝触媒出口におけるアンモニア流出量や一酸化窒素の流出量等から可能である。
【００５４】
ステップＳ２８では、脱硝反応の運転時間、すなわち脱硝触媒の劣化・再生予測を行う時間が経過したか否かを判定する。運転時間が経過していなければ、ステップＳ１４に処理を戻し、脱硝触媒の劣化・再生予測の演算を継続する。運転時間が経過していれば、脱硝触媒の劣化・再生予測の演算を終え、ステップＳ３０へ処理を移行させる。
【００５５】
ステップＳ３０では、脱硝触媒の劣化・再生予測の算出結果を出力する。脱硝触媒の劣化・再生予測の算出結果は、脱硝触媒の劣化再生状況（劣化関数Φの値、窒素酸化物（ＮＯ）の濃度の時間変化の減少量（増加量））、脱硝触媒の出口におけるアンモニア流量（濃度）や一酸化窒素流量（濃度）とすることができる。これらの算出結果は、脱硝反応の運転時間、脱硝触媒の劣化・再生予測を行う時間の関数として出力してもよい。このように、例えば、脱硝触媒の出口におけるアンモニア濃度、窒素酸化物の算出や劣化に至る運転時間の算出が容易となる。また、劣化した脱硝触媒の触媒温度を更新した場合の、脱硝触媒の再生の予測を正確に行うことができる。なお、算出結果の出力は、上記したものに限られるものでなく、また、処理途中の結果の出力を含めてもよく、出力の形式も各種の方法を採り得る。
【００５６】
＜脱硝触媒の劣化再生予測処理の検証及び評価＞
まず、式（１）の脱硝反応式のみを用いた解析を行い、モデル式の特徴を確認した。図５にＮＨ３／ＮＯ濃度比αによる脱硝率の温度影響について示す。ＫＮＨ３を０．０５×１０６に設定することによって、記号で示している実験結果と計算結果はほぼ一致した。
【００５７】
濃度比αが１ではアンモニア（ＮＨ３）と一酸化窒素（ＮＯ）が過不足なく反応するだけアンモニア（ＮＨ３）を導入していることになるが、温度が低下するに従って脱硝率は低下する。アンモニア（ＮＨ３）スリップを考慮すると、３００℃以下では反応しきれなかったアンモニア（ＮＨ３）が流出することになるので、濃度比α＝１の条件下で脱硝を行うとすると温度を３００℃以上にする必要がある。また、濃度比αが低下すると脱硝率も低下するが、温度に対して脱硝率の変化が起きなくなる。
【００５８】
そこで、最低雰囲気ガス温度を２５０℃に設定し、濃度比αを０．８程度に維持できるような触媒入口での一酸化窒素（ＮＯ）濃度とすることで、雰囲気ガス温度に対して濃度比αを変更する必要がなくなり、還元剤の投入制御性が向上すると考えられる。
【００５９】
次に、触媒中での滞留時間の影響を調べるために空塔速度（ＳＶ）を変更した時の脱硝率の影響を求める。図６にＳＶ変更時の濃度比αと脱硝率の関係を示す。温度は２００、２５０℃とした。図５と同様に温度が高い方が脱硝率は大きくなる。ＳＶが小さくなると滞留時間が長くなるため、反応時間が確保でき脱硝率が向上する。また、ＳＶが小さくなる程、濃度比αと脱硝率の線形関係が維持されることが確認できた。どの条件においても濃度比αを大きくしても脱硝率が変わらなくなる領域が存在する。この領域はアンモニア（ＮＨ３）スリップが生じていることになるため濃度比と脱硝率が線形関係にある範囲において還元剤の噴射制御を行うことが重要となる。
【００６０】
Ｖ２Ｏ５−ＴｉＯ２系触媒では触媒上で二酸化硫黄（ＳＯ２）が三酸化硫黄（ＳＯ３）に酸化されるため、化学式３に示されるように酸性硫安が生成され触媒性能が劣化する。本モデルでの触媒性能劣化予測の確認のために、１００時間劣化試験との比較を行った。図７に長時間運転時の脱硝率の変化の一例を示す。温度は、２００，２３０，２５０℃、二酸化硫黄（ＳＯ２）の濃度ＣＳＯ２は３０，８０ｐｐｍとした。
【００６１】
雰囲気ガス成分、温度等の様々な条件によって析出の有無は変わってくるが、どの条件でも実験結果とほぼ一致している。これより劣化モデルは本触媒での現象を再現できることが確認できた。なお、脱硝触媒の温度Ｔｃａｔ＝２５０℃、二酸化硫黄（ＳＯ２）の濃度ＣＳＯ２は８０ｐｐｍでの結果は劣化していないように見えるが、後述の図８の酸性硫安析出上限温度に示すように、性能劣化までには至っていないが酸性硫安の析出は起こっていることを確認した。
【００６２】
１００時間での劣化予測結果に見られるように、性能劣化が生じてしまうと長時間の運転中に徐々に脱硝率が低下してしまう。そのため、酸性硫安が生じるときの温度を調べることが重要になる。
【００６３】
ここではマン社（ＭＡＮ）の資料（佐々木，Ｋ．Ａａｂｏ、日マリ学誌、４３−３（２００８）３８２）にある燃料中硫黄含有率に対する酸性硫安析出上限温度を比較対照にして、本モデルにおける上限温度について確認した。
【００６４】
図８には、本モデルとマン社（ＭＡＮ）の資料の酸性硫安析出上限温度を示している。図に示す線以下の温度で酸性硫安は析出する。本計算における式（７）の細孔半径は、Ｒｏｓｉｎ−Ｒａｍｍｌｅｒ分布関数の１％通過径０．４ｎｍとして、析出しやすい半径を選択した。マン社（ＭＡＮ）の資料は燃料中硫黄含有率との関係を示しているが、硫黄は触媒には雰囲気ガス中の硫黄酸化物として導入されるため、硫黄酸化物を二酸化硫黄（ＳＯ２）、燃料を化学式（ＣａＨｂＳｃ：ａ＝１３，ｂ＝２８）で表される物質とし、空気過剰率をλ、完全燃焼を仮定して以下の化学式４から、燃料中の質量比である硫黄含有率と生成物のモル比で表される二酸化硫黄（ＳＯ２）濃度の関係を求めた。
【化４】
【００６５】
ただし、両者の関係はλによって変化するので、ここではλ＝２．５とした（図８の横軸に併せて示した）。およそ燃料中の硫黄分が３．３％で排気中の二酸化硫黄（ＳＯ２）の濃度が約８００ｐｐｍになることが確認できた。硫黄成分が少なくなると、酸性硫安は析出しにくくなる傾向がある。酸性硫安析出は式（５）のように硫黄成分の他にアンモニア（ＮＨ３）、水（Ｈ２Ｏ）の濃度によって影響を受けるが、二酸化硫黄（ＳＯ２）の濃度が８００ｐｐｍ、λ＝１において、アンモニア（ＮＨ３）の濃度を１５００，１０００，５００ｐｐｍに変更しても、上限温度は３２５，３２３，３１９℃であり、ほとんど影響しなかった。
【００６６】
次に、触媒性能劣化に及ぼすガス濃度の影響について確認した。図５及び図７に示したモデル確認のために行った実験での水（Ｈ２Ｏ）及び酸素（Ｏ２）の濃度はそれぞれ１０％及び１３％であった。また、ディーゼルエンジンの排ガス中の水（Ｈ２Ｏ）及び酸素（Ｏ２）の濃度は、化学式４の生成物のモル比から空気過剰率λを２．０〜４．０の範囲内の値とするとそれぞれ７．０７〜３．６０％及び１０．１〜１５．４％となる。そこで空気過剰率λを変更して水（Ｈ２Ｏ）及び酸素（Ｏ２）の濃度の脱硝率に与える影響について調べた。
【００６７】
図９に空気過剰率λが劣化速度に与える効果について示す。温度は２５０℃、ＳＶは１１０００ｈ−１とした。空気過剰率λの増加に従って劣化速度は減少した。式（７）の細孔半径に分布関数を使っているため、計算ステップ毎に酸性硫安析出によって消費されるアンモニア（ＮＨ３）の濃度が異なるものとなり、脱硝開始後３時間付近まで脱硝率が振動した。水（Ｈ２Ｏ）及び酸素（Ｏ２）の濃度が性能劣化に与える影響は、式（３）の二酸化硫黄（ＳＯ２）の酸化率の水（Ｈ２Ｏ）及び酸素（Ｏ２）の濃度と式（５）の酸性硫安平衡定数に関する水（Ｈ２Ｏ）の分圧によってモデルに組み込まれている。空気過剰率λの増加にともなって水（Ｈ２Ｏ）の濃度の減少、酸素（Ｏ２）濃度の増加が生じたが、式（３）内の反応次数は酸素（Ｏ２）よりも水（Ｈ２Ｏ）の方が大きいために二酸化硫黄（ＳＯ２）の酸化率は低下した。水（Ｈ２Ｏ）の濃度に加えて、三酸化硫黄（ＳＯ３）濃度も小さくなるため平衡定数も小さくなり、酸性硫安の生成量も減少した。
【００６８】
この結果は、燃料中の硫黄含有率を変えなくても過給機性能の向上による空気量の増加や燃料消費率の低減によって空気過剰率を増加させられれば、触媒性能の劣化を起こりづらくすることができることを予測している。
【００６９】
次に、二酸化硫黄（ＳＯ２）の濃度と一酸化窒素（ＮＯ）の濃度との変更による脱硝率の変化について調べた。図１０に脱硝開始後８時間での脱硝率の変化を示す。濃度比α＝１に固定したため、一酸化窒素（ＮＯ）の濃度と同時にアンモニア（ＮＨ３）の濃度も変更した。
【００７０】
一酸化窒素（ＮＯ）及びアンモニア（ＮＨ３）の濃度が大きくなると脱硝性能が劣化しやすいことがわかった。また、二酸化硫黄（ＳＯ２）の濃度だけではなく、一酸化窒素（ＮＯ）もしくはアンモニア（ＮＨ３）の濃度を低下させても劣化し難くなることがわかった。例えば、一酸化窒素（ＮＯ）の濃度を２００ｐｐｍにできれば二酸化硫黄（ＳＯ２）の濃度に依存せず、脱硝率は一定になった。また、一酸化窒素（ＮＯ）の濃度が小さくなると脱硝率が低下しているが、これは性能劣化の影響ではなく式（１）で表されるモデル自体の特性である。
【００７１】
一酸化窒素（ＮＯ）及びアンモニア（ＮＨ３）の濃度の低下によって触媒性能が改善することは、燃料性状の改善のみならず、エンジン本体による窒素酸化物（ＮＯｘ）低減によっても触媒寿命の長期化を計れることを示している。これより、エンジン本体での窒素酸化物（ＮＯｘ）の低減効果は、還元剤の減少によるコスト面での利点だけではなく、脱硝率、触媒劣化の面からも有効であることがわかった。
【００７２】
次に、一酸化窒素（ＮＯ）及びアンモニア（ＮＨ３）の濃度比αとアンモニア（ＮＨ３）スリップの関係について確認した。アンモニア（ＮＨ３）スリップは、式（１）の脱硝反応に用いられず、かつ触媒に吸着しきれなくなったアンモニア（ＮＨ３）が流出することで生じる。性能劣化が生じない場合には、温度を上げ、濃度比α≦１になるように還元剤噴射制御を行えばアンモニア（ＮＨ３）スリップは生じない。
【００７３】
本モデルにおいて、アンモニア（ＮＨ３）に関係する式（１）は、アンモニア（ＮＨ３）の吸着平衡後の反応速度を表しているため、反応しきれなくなったアンモニア（ＮＨ３）は触媒に吸着されることはなく、流出する。アンモニア（ＮＨ３）スリップは投入されたアンモニア（ＮＨ３）量に依存することが考えられるため、性能劣化していないときに脱硝触媒の出口における一酸化窒素（ＮＯ）の濃度が３００ｐｐｍになるように一酸化窒素（ＮＯ）の濃度と濃度比αを変更し、アンモニア（ＮＨ３）スリップに対する濃度比αの影響を調べた。濃度比αは、０．５，０．６６，０．７５，０．８，０．８２５，０．８５の６種類とした。
【００７４】
図１１は、脱硝開始から８時間後の濃度比αに対する脱硝触媒の出口における一酸化窒素（ＮＯ）の濃度の変化を示す。また、図１２は、アンモニア（ＮＨ３）スリップ濃度の時間変化を示す。
【００７５】
どの条件下においても脱硝開始時には出口における一酸化窒素（ＮＯ）の濃度は３００ｐｐｍであったが、濃度比αが大きい程、性能劣化の進行が早くなった。また濃度比α≧０．８になると急激に脱硝率が悪化した。濃度比α≧０．８で濃度比αが大きい程性能劣化しやすく、アンモニア（ＮＨ３）が消費できなくなるためにアンモニア（ＮＨ３）スリップが大きくなったものと考えられる。濃度比α＝０．８５では５時間以上この条件で運転できず、再生運転が必要となるといえる。濃度比α＜０．８ではほぼアンモニア（ＮＨ３）スリップが生じないのは、触媒の持つ脱硝性能以下で運転しているため、若干の触媒の性能劣化が生じてもアンモニア（ＮＨ３）を消費して脱硝性能を維持できるからと考えられる。これより、アンモニア（ＮＨ３）スリップの観点からもアンモニア（ＮＨ３）の投入量を減少できるようにエンジンでの窒素酸化物（ＮＯｘ）の低減対策を行う必要があることがわかった。
【００７６】
このようにＶ２Ｏ５−ＴｉＯ２系触媒の脱硝反応を再現できるモデルを用いることによって、硫黄含有燃料を使用したときの各種ガスの成分（濃度、濃度変化）、一酸化窒素（ＮＯ）及びアンモニア（ＮＨ３）の濃度比αに対する脱硝率、アンモニア（ＮＨ３）スリップへの影響について予測できる。
【００７７】
また、本実施の形態における脱硝触媒の劣化再生予測プログラムを用いた予測処理は、実際の脱硝システムに適用することができる。例えば、舶用ディーゼルエンジンに対する脱硝システムにおいて、実際に脱硝触媒に導入される尿素水の分解したアンモニア（ＮＨ３）や一酸化窒素（ＮＯ）及び脱硝触媒に関するパラメータを用いて脱硝触媒の劣化や再生を予測することができる。さらに、本発明の脱硝触媒の劣化再生予測プログラムを用いた脱硝システムによれば、脱硝触媒の劣化及び再生の少なくとも一方の予測ができ、脱硝システムの的確な運転に寄与することができる。例えば、予測結果に応じて脱硝システムの各種の操作量を変更し、操作量の変更に伴ってステップＳ１０、Ｓ１２及びＳ２６における少なくとも１つのパラメータの設定・変更を行い、予測を繰り返すものとしてもよい。また、例えば、脱硝触媒の劣化レベルの表示や、運転条件を変更した場合の窒素酸化物やアンモニアの排出状態の予測、脱硝触媒の交換や再生のタイミングの告知、再生温度と再生時間の予測等が容易となる。また、劣化再生プログラムの出力に基づいて脱硝システムの操作量を変更した場合は、状況に応じ脱硝システムを自動的に運転することも可能となる。
【産業上の利用可能性】
【００７８】
以上の実施形態は、舶用ディーゼルエンジンのみならず、他の内燃機関や燃焼機器の脱硝触媒システムにも適用することができる。
【技術分野】
【０００１】
本発明は、脱硝触媒の劣化再生予測プログラム及びそれを用いた脱硝システムに関する。
【背景技術】
【０００２】
国際海事機関（ＩＭＯ）によるＮＯｘ規制は現在の一次規制から強化されることが決定している。２０１１年には二次規制、２０１６年には三次規制となる予定だが、三次規制になるとエンジン本体だけではなく、後処理装置の導入によるＮＯｘ低減も必要になると考えられている。尿素水を用いた選択触媒還元（ＳＣＲ：Selective Catalytic Reduction）はＮＯｘ低減のための有力な後処理技術の一つである。ＳＣＲのメリットとしては、浄化率が高いこと、ＳＯｘを含んだ排ガス中でも利用ができること等が挙げられる。そのため、舶用ディーゼルエンジンへの適用が検討されているが、触媒性能の劣化やアンモニア（ＮＨ３）スリップ等の問題点もある。
【０００３】
三次規制は指定された海域内でのみ適用され、それ以外の海域では二次規制が適用されるが、この海域は基本的に港湾や輻輳海域等の船舶交通量の比較的大きく、陸地に近い場所に設定されると考えられる。その時、ＳＣＲはこの海域を航行する限られた時間のみで運用され、航行中の大部分はＳＣＲを作動させないことが想定できる。そうなると、ＳＣＲの運用に関する自由度は上がり、例えば劣化した触媒の再生をＮＯｘ還元の必要のない海域航行中に行うようなことも可能になるためＳＣＲの寿命や劣化再生方法について調べることが重要となってくる。しかし、これらの検討の手段として実験を行おうとすると非常に長い時間がかかることや温度や排ガス成分等のパラメータが非常に多いことから、数値解析による予測が有効な手段になると考えられる。
【０００４】
ＳＣＲを用いた脱硝反応に関して基礎的な実験・解析から数値解析まで，数多くの研究が行われている。
【０００５】
例えば、窒素酸化物を含有する排ガス中に還元剤であるアンモニアを投入する際の制御について開示されている（特許文献１，２等）。特許文献１では、触媒ユニット、エンジンの作動に関するパラメータを用いて制御が行われる。特許文献２では、窒素酸化物（ＮＯｘ）浄化特性の多項式を用いてアンモニア吸着量の目標値を設定する方法が採用されている。
【０００６】
また、脱硝システムにおける脱硝触媒の劣化を診断する診断装置が開示されている（特許文献３）。ここでは、アンモニアが漏れるまで脱硝触媒へ供給された還元剤の総量に基づいて脱硝触媒が劣化したか否かを判定している。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特表平０８−５０９７９５号公報
【特許文献２】特開２００８−１９６３４０号公報
【特許文献３】特開２００９−１２７４９６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかし、脱硝触媒の劣化を含めた計算を行っている例は少なく、特に、舶用ディーゼルエンジンを想定した比較的低排気温度で、高硫黄含有燃料を用いた場合を対象にして検討している例はない。さらに、脱硝触媒の劣化及び再生を総括する反応モデルを用いて脱硝触媒の劣化再生予測を行う技術はない。
【０００９】
本発明は、上記の課題を鑑み、脱硝触媒の劣化及び再生を総括する反応モデルを用いて、脱硝触媒の劣化再生予測プログラム及びそれを用いた脱硝システムを提供することを目的とする。
【課題を解決するための手段】
【００１０】
請求項１に対応した脱硝触媒の劣化再生予測プログラムは、コンピュータによって、燃焼機器からの排ガスにアンモニア系還元剤を添加して脱硝触媒で窒素酸化物の分解を行う脱硝システムにおける前記脱硝触媒の劣化と再生を予測するプログラムであって、前記アンモニア系還元剤のパラメータ設定を行うステップ１と、前記脱硝触媒に関する条件設定を行うステップ２と、前記脱硝触媒でのアンモニアの反応量及び窒素酸化物の反応量を算出するステップ３と、前記脱硝触媒での三酸化硫黄の生成量を算出するステップ４と、前記脱硝触媒での酸性硫安の増減量を算出するステップ５と、前記脱硝触媒の性能劣化量を算出するステップ６と、算出結果を出力するステップ７と、を備えることを特徴とする。
【００１１】
ここで、前記アンモニア系還元剤は、尿素水が好適であり、前記ステップ１における前記アンモニア系還元剤のパラメータは、前記尿素水の分解したアンモニアの濃度を含むことが好適である。これにより、還元剤がアンモニアである場合の脱硝システムにおける脱硝触媒の劣化又は再生を適切に予測することが可能となる。なお、アンモニア系還元剤としては、尿素水や尿素以外に、アンモニアや炭酸アンモニウム、その他アンモニア前駆体等を使用することができる。
【００１２】
また、前記ステップ２における条件は、前記脱硝触媒の温度及び体積と、前記排ガスの成分、流量及び温度と、であることが好適である。これにより、脱硝触媒のパラメータに応じて脱硝システムにおける脱硝触媒の劣化又は再生を適切に予測することが可能となる。
【００１３】
また、前記ステップ３は、Ｒｉｄｅａｌ−Ｅｌｅｙ機構を用いた脱硝反応速度式に基づいたアンモニアと窒素酸化物との反応量を算出する脱硝モデルを用いることが好適である。Ｒｉｄｅａｌ−Ｅｌｅｙ機構を用いた脱硝反応速度式を適用することにより、脱硝システムにおける脱硝触媒の劣化又は再生の予測の確度を高めることができる。
【００１４】
また、前記ステップ４は、非定常かつ、排ガス成分の影響を考慮した速度式に基づいた三酸化硫黄の生成量を算出する二酸化硫黄酸化反応モデルを用いることが好適である。これにより、三酸化硫黄の生成量に応じた脱硝システムにおける脱硝触媒の劣化又は再生を適切に予測することが可能となる。
【００１５】
また、前記ステップ５は、アンモニアと三酸化硫黄と水蒸気より酸性硫安を生じる平衡式の平衡定数を雰囲気ガス温度並びにアンモニア、三酸化硫黄及び水蒸気の分圧に基づいて定めた酸性硫安の増減量を算出する酸性硫安析出分解モデルを用いることを特徴とする。これにより、酸性硫安の増減量を考慮して脱硝システムにおける脱硝触媒の劣化又は再生を適切に予測することが可能となる。
【００１６】
また、前記酸性硫安析出分解モデルは、ケルビン方程式で補正されたものであることが好適である。これにより、酸性硫安の析出及び分解を正確に予測することができ、脱硝システムにおける脱硝触媒の劣化又は再生の予測の確度を高めることが可能となる。
【００１７】
また、前記ステップ６は、酸性硫安の増減量と劣化関数との関係を実験結果から求め、以下の式で算出される反応速度定数に基づいて行われることが好適である。これにより、脱硝触媒の劣化関数を考慮して、脱硝システムにおける脱硝触媒の劣化又は再生の予測の確度を高めることが可能となる。
ｋｃ＝ＡＮＯｅｘｐ（−ＥＮＯ／ＲＴｃａｔ）×Φ
ｋｃ：反応速度定数［ｍ３Ｎ／ｍ２ｈ］
ＡＮＯ、ＥＮＯ：ＮＯの反応に関する頻度因子［ｍ３Ｎ／ｍ２ｈ］
Ｒ：一般ガス定数［Ｊ／ｍｏｌＫ］
Ｔｃａｔ：触媒表面温度［Ｋ］
Φ：劣化関数
【００１８】
また、前記脱硝触媒の出口におけるアンモニア量と窒素酸化物量を算出するステップ８をさらに備えることが好適である。これにより、脱硝触媒の劣化又は再生を考慮して、脱硝処理におけるアンモニア量と窒素酸化物量の変化を把握することができる。
【００１９】
また、前記脱硝触媒の温度が所定の温度に達したか否かを判定するステップ９をさらに備えることが好適である。これにより、脱硝触媒の温度に応じた脱硝触媒の劣化又は再生をより正確に予測することが可能となる。
【００２０】
また、前記所定の温度を更新するステップ１０をさらに備えることが好適である。これにより、脱硝触媒の複数温度における脱硝触媒の劣化又は再生をより正確に予測することが可能となる。
【００２１】
また、前記脱硝システムの運転時間が所定の時間に達したか否かを判断するステップ１１をさらに備えることが好適である。これにより、脱硝システムの運転時間を考慮した脱硝触媒の劣化又は再生の予測を行うことが可能となる。
【００２２】
また、前記ステップ７は、時間と触媒性能劣化量の関係、時間と前記脱硝触媒の出口におけるアンモニア量の関係の少なくとも一方を出力することが好適である。これにより、脱硝システムの運転時間に応じた脱硝触媒の劣化又は再生を考慮して、触媒の性能劣化量、脱硝触媒の出口におけるアンモニア量を把握することができる。
【００２３】
請求項１４に対応した脱硝システムは、脱硝触媒の劣化再生予測プログラムを用いて上記脱硝触媒の劣化及び再生の少なくとも一方の予測を行うことが好適である。これにより、脱硝触媒の劣化及び再生の少なくとも一方の予測が可能な脱硝システムを実現することができる。
【００２４】
また、前記脱硝触媒の劣化再生予測プログラムの出力に基づいて、前記脱硝システムの操作量を変更し、前記ステップ１、２及び１０の少なくとも１つの設定を変更することが好適である。これにより、実際の脱硝システムにおける操作量に応じた脱硝触媒の劣化又は再生を考慮し、状況に応じた脱硝システムの制御を行うことが可能となる。
【発明の効果】
【００２５】
本発明の脱硝触媒の劣化再生予測プログラムによれば、アンモニア系還元剤を添加して脱硝触媒で窒素酸化物の分解を行う際の脱硝反応をより正確に予測することができる。また、アンモニア系還元剤のパラメータや脱硝触媒に関する条件を変えた場合の脱硝触媒の性能劣化量の算出が容易にできる。
これにより、例えば、脱硝触媒の出口におけるアンモニア濃度、窒素酸化物の算出や劣化に至る運転時間の算出が容易となる。また、劣化した脱硝触媒の触媒温度を更新した場合の、脱硝触媒の再生の予測を正確に行うことができる。
また、本発明の脱硝触媒の劣化再生予測プログラムを用いた脱硝システムによれば、脱硝触媒の劣化及び再生の少なくとも一方の予測ができ、脱硝システムの的確な運転に寄与することができる。例えば、脱硝触媒の劣化レベルの表示や、運転条件を変更した場合の窒素酸化物やアンモニアの排出状態の予測、脱硝触媒の交換や再生のタイミングの告知、再生温度と再生時間の予測等が容易となる。
また、劣化再生プログラムの出力に基づいて脱硝システムの操作量を変更した場合は、状況に応じ脱硝システムを自動的に運転することも可能となる。
【図面の簡単な説明】
【００２６】
【図１】本発明の実施の形態における脱硝触媒の劣化再生予測システムの構成を示す図である。
【図２】本発明の実施の形態における脱硝触媒の劣化再生予測処理のフローチャートを示す図である。
【図３】脱硝反応を説明する図である。
【図４】酸性硫安付着量と劣化関数Φとの関係を示す図である。
【図５】濃度比αに対する脱硝率と触媒温度との関係を示す図である。
【図６】ＳＶ，温度変更時の脱硝率と濃度比αとの関係を示す図である。
【図７】脱硝率の時間変化の実験値及び計算値との比較を示す図である。
【図８】硫黄成分と酸性硫安析出上限温度との関係を示す図である。
【図９】空気過剰率による脱硝率の時間変化を示す図である。
【図１０】ＮＯ，ＳＯ２濃度が脱硝率に与える影響を示す図である。
【図１１】脱硝触媒の出口ＮＯｘ濃度に及ぼす濃度比αの効果を示す図である。
【図１２】ＮＨ３スリップ濃度の時間変化を示す図である。
【発明を実施するための形態】
【００２７】
＜脱硝触媒の劣化再生予測処理＞
本発明の実施の形態における脱硝触媒の劣化再生予測処理は、コンピュータシステム１００を用いて実行される。コンピュータシステム１００は、図１に示すように、処理部１０、記憶部１２、入力部１４及び出力部１６を含んで構成される。
【００２８】
処理部１０は、ＣＰＵ等で構成され、記憶部１２に予め記憶された脱硝触媒の劣化再生予測プログラムを実行することにより脱硝触媒の劣化再生予測処理の演算を行う。記憶部１２は、脱硝触媒の劣化再生予測プログラム及びそれに関する各種パラメータを格納及び保持する。記憶部１２は、半導体メモリ、ハードディスク等の記憶装置を含んで構成され、処理部１０から読み出し及び書き込み可能とされる。入力部１４は、脱硝触媒の劣化再生予測処理に必要なデータ等を入力するために用いられる。入力部１４は、キーボード、ポインティングデバイス等の入力装置を含んで構成される。出力部１６は、入力部１４から入力されたデータの表示、脱硝触媒の劣化再生予測処理で得られる中間処理結果や最終処理結果の表示等を行う。出力部１６は、ディスプレイ、プリンタ等の出力装置を含んで構成される。
【００２９】
以下、本発明の実施の形態における脱硝触媒の劣化再生予測処理について説明する。脱硝触媒の劣化再生予測処理は、処理部１０によって記憶部１２に記憶されている脱硝触媒の劣化再生予測プログラムを実行することによって、図２に示すフローチャートに沿って行われる。
【００３０】
本実施の形態では、脱硝触媒内の反応のみを考慮した一次元の簡易的な反応律速モデルを用いている。すなわち、脱硝触媒中の各物質は流れ直角方向に十分拡散し、流れ方向のみに変化することを仮定している。また、一般的に脱硝触媒には還元剤として尿素水が導入されるが、尿素水は加水分解後、アンモニア（ＮＨ３）ガスとして触媒に導入されるとした。また窒素酸化物（ＮＯｘ）については一酸化窒素（ＮＯ）のみを仮定した。考慮する物質は、アンモニア（ＮＨ３）、一酸化窒素（ＮＯ）、二酸化硫黄（ＳＯ２）、三酸化硫黄（ＳＯ３）、水（Ｈ２Ｏ）、窒素（Ｎ２）、酸素（Ｏ２）及び硫酸水素アンモニウム（ＮＨ４ＨＳＯ４：以下、酸性硫安と称する）の８種類とする。
【００３１】
本実施の形態における脱硝触媒の劣化再生予測処理では、各物質の流量は脱硝触媒を流れ方向に分割した各ブロックに存在できる時間毎に分割され、そのブロックに各物質の存在する流量（体積）に対して触媒反応の計算を行う。すなわち、ある時刻ステップ（ｔ）における各ブロックにおける脱硝触媒による脱硝反応、二硫化硫黄の酸化反応、酸性硫安の生成分解反応及び酸性硫安による脱硝触媒の劣化・再生反応の計算を行う。算出された各物質は、酸性硫安以外はそのブロックに何も残さず、次の時刻ステップ（ｔ＋１）には次のブロックに輸送されるものとする。そして、算出された各物質の濃度等を用いて、次の時刻ステップ（ｔ＋１）における各ブロックにおける触媒反応の計算を行う。このような処理を繰り返し、脱硝触媒の劣化再生予測処理を実行する。
【００３２】
なお、各ブロックでの滞留時間は、雰囲気ガス流量［ｍ３Ｎ／ｈ］を触媒体積［ｍ３］で除した空塔速度（ＳＶ：Space Velocity）の逆数として定義した。ただし、二次元の反応律速モデルに拡張することも可能である。
【００３３】
ステップＳ１０では、脱硝触媒の劣化再生予測処理で用いられるパラメータの初期設定を行う。ここでのパラメータの設定には、還元剤及び窒素酸化物に関するパラメータの設定が含まれる。上記のように還元剤として尿素水を加水分解したアンモニア（ＮＨ３）ガスが使用されるものとして、還元剤に関するパラメータは、反応速度定数ＫＮＨ３［ｍ３Ｎ／ｍ２ｈ］、アンモニアの体積濃度ＣＮＨ３［ｍ３／ｍ３］、アンモニア吸着平衡定数ＫＮＨ３等の後述の反応モデルで使用されるパラメータである。また、窒素酸化物に関するパラメータは、反応速度定数ＫＮＯ［ｍ３Ｎ／ｍ２ｈ］、窒素酸化物の体積濃度ＣＮＯ［ｍ３／ｍ３］、窒素酸化物の反応に関する頻度因子ＡＮＯ［ｍ３Ｎ／ｍ２ｈ］及び活性化エネルギＥＮＯ［Ｊ／ｍｏｌ］等の後述の反応モデルで使用されるパラメータである。また、二酸化硫黄に関する反応速度定数ｋＳ、酸素に関する反応次数αＯ２、水に関する反応次数αＷ、一酸化窒素及びアンモニアに対する係数ｋＮＯ，ｋＮＨ３、二硫化硫黄に関する頻度因子ＡＳＯ２，活性化エネルギＥＳＯ２が設定される。このようにアンモニア系還元剤のパラメータや脱硝触媒に関する条件を設定することによって、これらの設定を変更した場合の脱硝触媒の性能劣化量の算出が容易にできる。なお、使用するアンモニア系還元剤の種類により、使用するパラメータは異なる。また、尿素水にしても、尿素の含有形態や供給方法、パラメータの表現方法等により、各種のパラメータやその組み合わせがあり得る。
【００３４】
例えば、反応速度定数ＫＮＨ３は後述する実験結果と一致するように０．０５×１０６に設定することが好適である。一酸化窒素（ＮＯ）とアンモニア（ＮＨ３）の濃度比α（ＣＮＨ３／ＣＮＯ）＝１のときの実験結果より、頻度因子ＡＮＯ及び活性化エネルギＥＮＯは、それぞれ８８９８ｍ３Ｎ／ｍ２ｈ及び２３．８７×１０３Ｊ／ｍｏｌとすることが好適である。また、窒素酸化物に関するパラメータは、反応次数αＯ２は文献(E.Tronconi,et.al.,Ind.Eng.Chem.Res.38(1999)2593)から設定することができ、その他は実験結果等から表１のように設定することができる。
【表１】
【００３５】
ステップＳ１２では、脱硝触媒に関するパラメータの設定を行う。ここでは、脱硝触媒の比表面積比Ａｐ［ｍ３／ｍ３］、脱硝触媒の表面温度Ｔｃａｔ［Ｋ］、脱硝触媒の予測処理を行う温度範囲等が設定される。脱硝触媒の表面温度Ｔｃａｔは初期設定であり、脱硝触媒の予測処理を行う温度範囲や脱硝触媒の再生処理温度に応じてステップ２６において変更される。なお、脱硝触媒に関するパラメータの設定は、使用する脱硝触媒により、温度範囲や初期設定される表面温度等は異なる。また、空塔速度や周囲温度等を設定してもよい。
【００３６】
ステップＳ１４では、脱硝触媒でのアンモニアの反応量及び窒素酸化物の反応量を算出する。Ｖ２Ｏ５−ＴｉＯ２系触媒における尿素水加水分解後のアンモニアによる代表的な脱硝反応の化学反応式は、化学式１のように表される。
【化１】
【００３７】
その他に二酸化窒素（ＮＯ２）と一酸化窒素（ＮＯ）が等モルで、化学式１よりも低温、高速で起こる反応もあるが、ディーゼルエンジン排ガス中の窒素酸化物（ＮＯｘ）のほとんどは一酸化窒素（ＮＯ）であるので、酸化触媒による二酸化窒素（ＮＯ２）の生成を考慮せず、ここでは化学式１の反応のみを考える。この反応が触媒表面で起こっているとすると、図３の概念図に示すように、脱硝反応はアンモニア（ＮＨ３）の触媒活性点への吸着、吸着したアンモニア（ＮＨ３）と雰囲気に拡散している一酸化窒素（ＮＯ）との反応、生成物の活性点からの脱離と進行すると考えられる。
【００３８】
脱硝反応速度式は、アンモニア（ＮＨ３）の影響を考慮する必要があるため、脱硝機構をモデル化する手法としてＲｉｄｅａｌ−Ｅｌｅｙ機構を用いた式（１）及び（２）により記述される。
【数１】
【数２】
【００３９】
この反応モデルにより、脱硝触媒でのアンモニアの反応量及び窒素酸化物の反応量を算出する。すなわち、式（１）及び（２）から窒素酸化物（ＮＯ）の体積濃度ＣＮＯの時間的な変化を算出することができる。化学式１から窒素酸化物（ＮＯ）の体積濃度ＣＮＯの時間的な変化に対応したアンモニア（ＮＨ３）の体積濃度ＣＮＨ３の時間的な変化も算出することができる。
【００４０】
ここで、Φは脱硝触媒の劣化関数である。劣化関数Φは、脱硝触媒の劣化後の反応速度定数ｋｃと初期反応速度定数ｋｃ０の比（ｋｃ／ｋｃ０）と定義する。劣化関数Φについては後述する。劣化関数Φの初期値は予め設定しておけばよく、例えば、脱硝触媒がまったく劣化していない状態では劣化関数Φの値は１とすればよい。なお、アンモニアの反応量及び窒素酸化物の反応量の算出は、上記したＲｉｄｅａｌ−Ｅｌｅｙ機構を用いた式（１）及び（２）以外の簡略式や厳密式、また、パラメータの設定や表現方法を変えた各種の式に基づくことができる。
【００４１】
ステップＳ１６では、脱硝触媒での三酸化硫黄の生成量を算出する。燃料中に硫黄が存在する場合には、エンジン内の燃焼によって硫黄が酸化され二酸化硫黄（ＳＯ２）になり、更に二酸化硫黄（ＳＯ２）が触媒上で酸化され三酸化硫黄（ＳＯ３）になる。二酸化硫黄（ＳＯ２）の酸化反応は、化学式２で表される。
【化２】
【００４２】
ここで、二酸化硫黄（ＳＯ２）の酸化反応は、Ｖ２Ｏ５−ＴｉＯ２系触媒において窒素酸化物の脱硝反応と同時に起こり、脱硝反応に比べて非常に遅いことが知られている。この二酸化硫黄（ＳＯ２）から三酸化硫黄（ＳＯ３）への転換率は数％程度である。しかしながら、酸化した三酸化硫黄（ＳＯ３）は後述するようにアンモニア（ＮＨ３）及び水（Ｈ２Ｏ）との反応により酸性硫安になると触媒表面に付着して活性吸着点を覆うため、脱硝反応が阻害される。そのため、劣化を予測するためには二酸化硫黄（ＳＯ２）の酸化反応を予測する。
【００４３】
ここでは、二酸化硫黄（ＳＯ２）の反応モデルとして非定常での計算ができる速度式として式（３）及び（４）を用いる。なお、三酸化硫黄の生成量の算出は式（３），（４）に限定されるものでなく、設定パラメータや表現方法を変えた式、簡略式や厳密式等、各種の式に基づくことができる。
【数３】
【数４】
【００４４】
ステップＳ１８では、脱硝触媒での酸性硫安の増減量を算出する。ここで、舶用ディーゼルエンジンのような硫黄を含む燃料を用いるシステムの場合、ＳＣＲによる脱硝を行おうとすると、未反応のアンモニア（ＮＨ３）と排ガス中の三酸化硫黄（ＳＯ３）及び水（Ｈ２Ｏ）の反応によって生成される化合物が触媒上に付着し、触媒の活性を低下させることが問題となる。この触媒の劣化の有無及び進行速度は、温度によっても変わってくる。ガス中にアンモニア（ＮＨ３）、三酸化硫黄（ＳＯ３）及び水（Ｈ２Ｏ）が共存している場合に生成される化合物は、硫酸アンモニウム（（ＮＨ４）２ＳＯ４）、酸性硫安等が考えられるが、ここではその化合物を酸性硫安のみであるとして演算する。
【化３】
【数５】
【数６】
【００４５】
ここで、Ｋは平衡定数、ｐはアンモニア（ＮＨ３）、三酸化硫黄（ＳＯ３）及び水（Ｈ２Ｏ）の分圧、Ｔａは雰囲気ガス温度を示している。式（５）で示される分圧が式（６）の温度による平衡定数より大きい場合に析出が生じる。ここで、表面張力と密度は文献等から得られる値とし、細孔半径は実測値をＲｏｓｉｎ−Ｒａｍｍｌｅｒの分布関数で近似した値とすることが好適である。
【００４６】
しかし、この式は気体中での反応を示している。それを無数の細孔のある触媒中で適用するため、式（５）をケルビン方程式によって修正する。この式は、細孔内では気体中の飽和蒸気圧ＰＯよりも低い圧力Ｐによって凝縮が起こることを示す。なお、式（６）の平衡定数Ｋが式（７）の飽和蒸気圧ＰＯに相当する。
【数７】
ここで、σは表面張力、Ｍａｂ及びρａｂは酸性硫安の分子量及び密度、ｒｐｏｒｅは細孔半径である。
【００４７】
ここでは、二酸化硫黄（ＳＯ２）の反応モデルを用いて得られた二酸化硫黄（ＳＯ２）の体積濃度ＣＳＯ２の時間変化から三酸化硫黄（ＳＯ３）の体積濃度ＣＳＯ３の時間変化を求め、その算出結果を式（５）〜（７）に導入することにより酸性硫安の増加・減少量を求める。なお、酸性硫安の増減量の算出は、式（５）〜（７）に限定されるものではなく、設定パラメータや表現方法を変えた式、簡略式や厳密式等、各種の式に基づくことができる。
【００４８】
ステップＳ２０では、脱硝触媒の性能劣化再生量を算出する。脱硝触媒の劣化・再生は、脱硝触媒上に酸性硫安が付着することによって生ずる。すなわち、脱硝触媒の性能劣化再生量は、酸性硫安付着による触媒性能の劣化関数Φの値を式（２）の一酸化窒素（ＮＯ）の反応速度定数に乗算し、窒素酸化物（ＮＯ）の濃度の時間変化の減少量（増加量）として算出される。
【００４９】
ここで、劣化関数Φの値は、酸性硫安の付着量の関係から求めることができる。ここでは、図４に示す実験結果から劣化関数Φの値を求めるものとする。すなわち、横軸の触媒１ｇあたりの硫酸イオン付着量［ｍｍｏｌ／ｇｃａｔ］を酸性硫安付着量と置き換え、それに対応する劣化関数Φの値を決定する。
【００５０】
具体的には、図４中の実線で示すように酸性硫安が０．５ｍｍｏｌ／ｇｃａｔ以下の時には５次の近似式で表し、０．５ｍｍｏｌ／ｇｃａｔ以上の時には１次の近似式で表すことが好適である。なお、酸性硫安は二酸化硫黄（ＳＯ２）の酸化反応に影響を与えると考えられるが、付着量と劣化関数Φの関係が明らかでないため、ここでは影響を考慮しないものとする。なお、脱硝触媒の性能劣化再生量の算出に当たっては、他の実験結果に基づいた劣化関数や他の近似式を用いてもよい。また、脱硝触媒の性能劣化量のみを求めることも可能である。
【００５１】
ステップＳ２２では、脱硝触媒の出口でのアンモニア流出量（体積濃度ＣＮＨ３）及び一酸化窒素（ＮＯ）流出量（体積濃度ＣＮＯ）の算出が行われる。すなわち、ステップＳ２０で求められた劣化関数Φを式（２）に導入し、式（１）及び（２）の関係式から脱硝触媒の出口でのアンモニア流出量及び一酸化窒素（ＮＯ）流出量を算出する。
【００５２】
ステップＳ２４では、脱硝触媒の温度を変化させるか否かを判定する。本実施の形態では、脱硝触媒の温度に応じた脱硝触媒の劣化・再生予測を行うために、設定された脱硝触媒の温度範囲での予測処理が終了したか否かを判定する。脱硝触媒の温度範囲及び変更ステップ等は、エンジンの運転条件の変更や再生処理の条件に応じて予め設定しておけばよい。
【００５３】
脱硝触媒の劣化・再生予測を行う脱硝触媒の温度が存在すれば、ステップＳ２６において次の脱硝触媒の温度に設定し、ステップＳ１４に処理を戻し、新たに設定された脱硝触媒の温度での演算を行う。例えば、脱硝触媒の温度を初期値から更新ステップ温度毎に増加させ、所定の演算最高温度に到達するまで脱硝触媒の温度を更新する。一方、必要な脱硝触媒の温度のすべてについて演算が終了していれば、ステップＳ２８に処理を移行させる。なお、温度を上げて脱硝触媒の再生を行う時にも、脱硝触媒の温度を更新する。再生の判断は、運転時間や脱硝触媒出口におけるアンモニア流出量や一酸化窒素の流出量等から可能である。
【００５４】
ステップＳ２８では、脱硝反応の運転時間、すなわち脱硝触媒の劣化・再生予測を行う時間が経過したか否かを判定する。運転時間が経過していなければ、ステップＳ１４に処理を戻し、脱硝触媒の劣化・再生予測の演算を継続する。運転時間が経過していれば、脱硝触媒の劣化・再生予測の演算を終え、ステップＳ３０へ処理を移行させる。
【００５５】
ステップＳ３０では、脱硝触媒の劣化・再生予測の算出結果を出力する。脱硝触媒の劣化・再生予測の算出結果は、脱硝触媒の劣化再生状況（劣化関数Φの値、窒素酸化物（ＮＯ）の濃度の時間変化の減少量（増加量））、脱硝触媒の出口におけるアンモニア流量（濃度）や一酸化窒素流量（濃度）とすることができる。これらの算出結果は、脱硝反応の運転時間、脱硝触媒の劣化・再生予測を行う時間の関数として出力してもよい。このように、例えば、脱硝触媒の出口におけるアンモニア濃度、窒素酸化物の算出や劣化に至る運転時間の算出が容易となる。また、劣化した脱硝触媒の触媒温度を更新した場合の、脱硝触媒の再生の予測を正確に行うことができる。なお、算出結果の出力は、上記したものに限られるものでなく、また、処理途中の結果の出力を含めてもよく、出力の形式も各種の方法を採り得る。
【００５６】
＜脱硝触媒の劣化再生予測処理の検証及び評価＞
まず、式（１）の脱硝反応式のみを用いた解析を行い、モデル式の特徴を確認した。図５にＮＨ３／ＮＯ濃度比αによる脱硝率の温度影響について示す。ＫＮＨ３を０．０５×１０６に設定することによって、記号で示している実験結果と計算結果はほぼ一致した。
【００５７】
濃度比αが１ではアンモニア（ＮＨ３）と一酸化窒素（ＮＯ）が過不足なく反応するだけアンモニア（ＮＨ３）を導入していることになるが、温度が低下するに従って脱硝率は低下する。アンモニア（ＮＨ３）スリップを考慮すると、３００℃以下では反応しきれなかったアンモニア（ＮＨ３）が流出することになるので、濃度比α＝１の条件下で脱硝を行うとすると温度を３００℃以上にする必要がある。また、濃度比αが低下すると脱硝率も低下するが、温度に対して脱硝率の変化が起きなくなる。
【００５８】
そこで、最低雰囲気ガス温度を２５０℃に設定し、濃度比αを０．８程度に維持できるような触媒入口での一酸化窒素（ＮＯ）濃度とすることで、雰囲気ガス温度に対して濃度比αを変更する必要がなくなり、還元剤の投入制御性が向上すると考えられる。
【００５９】
次に、触媒中での滞留時間の影響を調べるために空塔速度（ＳＶ）を変更した時の脱硝率の影響を求める。図６にＳＶ変更時の濃度比αと脱硝率の関係を示す。温度は２００、２５０℃とした。図５と同様に温度が高い方が脱硝率は大きくなる。ＳＶが小さくなると滞留時間が長くなるため、反応時間が確保でき脱硝率が向上する。また、ＳＶが小さくなる程、濃度比αと脱硝率の線形関係が維持されることが確認できた。どの条件においても濃度比αを大きくしても脱硝率が変わらなくなる領域が存在する。この領域はアンモニア（ＮＨ３）スリップが生じていることになるため濃度比と脱硝率が線形関係にある範囲において還元剤の噴射制御を行うことが重要となる。
【００６０】
Ｖ２Ｏ５−ＴｉＯ２系触媒では触媒上で二酸化硫黄（ＳＯ２）が三酸化硫黄（ＳＯ３）に酸化されるため、化学式３に示されるように酸性硫安が生成され触媒性能が劣化する。本モデルでの触媒性能劣化予測の確認のために、１００時間劣化試験との比較を行った。図７に長時間運転時の脱硝率の変化の一例を示す。温度は、２００，２３０，２５０℃、二酸化硫黄（ＳＯ２）の濃度ＣＳＯ２は３０，８０ｐｐｍとした。
【００６１】
雰囲気ガス成分、温度等の様々な条件によって析出の有無は変わってくるが、どの条件でも実験結果とほぼ一致している。これより劣化モデルは本触媒での現象を再現できることが確認できた。なお、脱硝触媒の温度Ｔｃａｔ＝２５０℃、二酸化硫黄（ＳＯ２）の濃度ＣＳＯ２は８０ｐｐｍでの結果は劣化していないように見えるが、後述の図８の酸性硫安析出上限温度に示すように、性能劣化までには至っていないが酸性硫安の析出は起こっていることを確認した。
【００６２】
１００時間での劣化予測結果に見られるように、性能劣化が生じてしまうと長時間の運転中に徐々に脱硝率が低下してしまう。そのため、酸性硫安が生じるときの温度を調べることが重要になる。
【００６３】
ここではマン社（ＭＡＮ）の資料（佐々木，Ｋ．Ａａｂｏ、日マリ学誌、４３−３（２００８）３８２）にある燃料中硫黄含有率に対する酸性硫安析出上限温度を比較対照にして、本モデルにおける上限温度について確認した。
【００６４】
図８には、本モデルとマン社（ＭＡＮ）の資料の酸性硫安析出上限温度を示している。図に示す線以下の温度で酸性硫安は析出する。本計算における式（７）の細孔半径は、Ｒｏｓｉｎ−Ｒａｍｍｌｅｒ分布関数の１％通過径０．４ｎｍとして、析出しやすい半径を選択した。マン社（ＭＡＮ）の資料は燃料中硫黄含有率との関係を示しているが、硫黄は触媒には雰囲気ガス中の硫黄酸化物として導入されるため、硫黄酸化物を二酸化硫黄（ＳＯ２）、燃料を化学式（ＣａＨｂＳｃ：ａ＝１３，ｂ＝２８）で表される物質とし、空気過剰率をλ、完全燃焼を仮定して以下の化学式４から、燃料中の質量比である硫黄含有率と生成物のモル比で表される二酸化硫黄（ＳＯ２）濃度の関係を求めた。
【化４】
【００６５】
ただし、両者の関係はλによって変化するので、ここではλ＝２．５とした（図８の横軸に併せて示した）。およそ燃料中の硫黄分が３．３％で排気中の二酸化硫黄（ＳＯ２）の濃度が約８００ｐｐｍになることが確認できた。硫黄成分が少なくなると、酸性硫安は析出しにくくなる傾向がある。酸性硫安析出は式（５）のように硫黄成分の他にアンモニア（ＮＨ３）、水（Ｈ２Ｏ）の濃度によって影響を受けるが、二酸化硫黄（ＳＯ２）の濃度が８００ｐｐｍ、λ＝１において、アンモニア（ＮＨ３）の濃度を１５００，１０００，５００ｐｐｍに変更しても、上限温度は３２５，３２３，３１９℃であり、ほとんど影響しなかった。
【００６６】
次に、触媒性能劣化に及ぼすガス濃度の影響について確認した。図５及び図７に示したモデル確認のために行った実験での水（Ｈ２Ｏ）及び酸素（Ｏ２）の濃度はそれぞれ１０％及び１３％であった。また、ディーゼルエンジンの排ガス中の水（Ｈ２Ｏ）及び酸素（Ｏ２）の濃度は、化学式４の生成物のモル比から空気過剰率λを２．０〜４．０の範囲内の値とするとそれぞれ７．０７〜３．６０％及び１０．１〜１５．４％となる。そこで空気過剰率λを変更して水（Ｈ２Ｏ）及び酸素（Ｏ２）の濃度の脱硝率に与える影響について調べた。
【００６７】
図９に空気過剰率λが劣化速度に与える効果について示す。温度は２５０℃、ＳＶは１１０００ｈ−１とした。空気過剰率λの増加に従って劣化速度は減少した。式（７）の細孔半径に分布関数を使っているため、計算ステップ毎に酸性硫安析出によって消費されるアンモニア（ＮＨ３）の濃度が異なるものとなり、脱硝開始後３時間付近まで脱硝率が振動した。水（Ｈ２Ｏ）及び酸素（Ｏ２）の濃度が性能劣化に与える影響は、式（３）の二酸化硫黄（ＳＯ２）の酸化率の水（Ｈ２Ｏ）及び酸素（Ｏ２）の濃度と式（５）の酸性硫安平衡定数に関する水（Ｈ２Ｏ）の分圧によってモデルに組み込まれている。空気過剰率λの増加にともなって水（Ｈ２Ｏ）の濃度の減少、酸素（Ｏ２）濃度の増加が生じたが、式（３）内の反応次数は酸素（Ｏ２）よりも水（Ｈ２Ｏ）の方が大きいために二酸化硫黄（ＳＯ２）の酸化率は低下した。水（Ｈ２Ｏ）の濃度に加えて、三酸化硫黄（ＳＯ３）濃度も小さくなるため平衡定数も小さくなり、酸性硫安の生成量も減少した。
【００６８】
この結果は、燃料中の硫黄含有率を変えなくても過給機性能の向上による空気量の増加や燃料消費率の低減によって空気過剰率を増加させられれば、触媒性能の劣化を起こりづらくすることができることを予測している。
【００６９】
次に、二酸化硫黄（ＳＯ２）の濃度と一酸化窒素（ＮＯ）の濃度との変更による脱硝率の変化について調べた。図１０に脱硝開始後８時間での脱硝率の変化を示す。濃度比α＝１に固定したため、一酸化窒素（ＮＯ）の濃度と同時にアンモニア（ＮＨ３）の濃度も変更した。
【００７０】
一酸化窒素（ＮＯ）及びアンモニア（ＮＨ３）の濃度が大きくなると脱硝性能が劣化しやすいことがわかった。また、二酸化硫黄（ＳＯ２）の濃度だけではなく、一酸化窒素（ＮＯ）もしくはアンモニア（ＮＨ３）の濃度を低下させても劣化し難くなることがわかった。例えば、一酸化窒素（ＮＯ）の濃度を２００ｐｐｍにできれば二酸化硫黄（ＳＯ２）の濃度に依存せず、脱硝率は一定になった。また、一酸化窒素（ＮＯ）の濃度が小さくなると脱硝率が低下しているが、これは性能劣化の影響ではなく式（１）で表されるモデル自体の特性である。
【００７１】
一酸化窒素（ＮＯ）及びアンモニア（ＮＨ３）の濃度の低下によって触媒性能が改善することは、燃料性状の改善のみならず、エンジン本体による窒素酸化物（ＮＯｘ）低減によっても触媒寿命の長期化を計れることを示している。これより、エンジン本体での窒素酸化物（ＮＯｘ）の低減効果は、還元剤の減少によるコスト面での利点だけではなく、脱硝率、触媒劣化の面からも有効であることがわかった。
【００７２】
次に、一酸化窒素（ＮＯ）及びアンモニア（ＮＨ３）の濃度比αとアンモニア（ＮＨ３）スリップの関係について確認した。アンモニア（ＮＨ３）スリップは、式（１）の脱硝反応に用いられず、かつ触媒に吸着しきれなくなったアンモニア（ＮＨ３）が流出することで生じる。性能劣化が生じない場合には、温度を上げ、濃度比α≦１になるように還元剤噴射制御を行えばアンモニア（ＮＨ３）スリップは生じない。
【００７３】
本モデルにおいて、アンモニア（ＮＨ３）に関係する式（１）は、アンモニア（ＮＨ３）の吸着平衡後の反応速度を表しているため、反応しきれなくなったアンモニア（ＮＨ３）は触媒に吸着されることはなく、流出する。アンモニア（ＮＨ３）スリップは投入されたアンモニア（ＮＨ３）量に依存することが考えられるため、性能劣化していないときに脱硝触媒の出口における一酸化窒素（ＮＯ）の濃度が３００ｐｐｍになるように一酸化窒素（ＮＯ）の濃度と濃度比αを変更し、アンモニア（ＮＨ３）スリップに対する濃度比αの影響を調べた。濃度比αは、０．５，０．６６，０．７５，０．８，０．８２５，０．８５の６種類とした。
【００７４】
図１１は、脱硝開始から８時間後の濃度比αに対する脱硝触媒の出口における一酸化窒素（ＮＯ）の濃度の変化を示す。また、図１２は、アンモニア（ＮＨ３）スリップ濃度の時間変化を示す。
【００７５】
どの条件下においても脱硝開始時には出口における一酸化窒素（ＮＯ）の濃度は３００ｐｐｍであったが、濃度比αが大きい程、性能劣化の進行が早くなった。また濃度比α≧０．８になると急激に脱硝率が悪化した。濃度比α≧０．８で濃度比αが大きい程性能劣化しやすく、アンモニア（ＮＨ３）が消費できなくなるためにアンモニア（ＮＨ３）スリップが大きくなったものと考えられる。濃度比α＝０．８５では５時間以上この条件で運転できず、再生運転が必要となるといえる。濃度比α＜０．８ではほぼアンモニア（ＮＨ３）スリップが生じないのは、触媒の持つ脱硝性能以下で運転しているため、若干の触媒の性能劣化が生じてもアンモニア（ＮＨ３）を消費して脱硝性能を維持できるからと考えられる。これより、アンモニア（ＮＨ３）スリップの観点からもアンモニア（ＮＨ３）の投入量を減少できるようにエンジンでの窒素酸化物（ＮＯｘ）の低減対策を行う必要があることがわかった。
【００７６】
このようにＶ２Ｏ５−ＴｉＯ２系触媒の脱硝反応を再現できるモデルを用いることによって、硫黄含有燃料を使用したときの各種ガスの成分（濃度、濃度変化）、一酸化窒素（ＮＯ）及びアンモニア（ＮＨ３）の濃度比αに対する脱硝率、アンモニア（ＮＨ３）スリップへの影響について予測できる。
【００７７】
また、本実施の形態における脱硝触媒の劣化再生予測プログラムを用いた予測処理は、実際の脱硝システムに適用することができる。例えば、舶用ディーゼルエンジンに対する脱硝システムにおいて、実際に脱硝触媒に導入される尿素水の分解したアンモニア（ＮＨ３）や一酸化窒素（ＮＯ）及び脱硝触媒に関するパラメータを用いて脱硝触媒の劣化や再生を予測することができる。さらに、本発明の脱硝触媒の劣化再生予測プログラムを用いた脱硝システムによれば、脱硝触媒の劣化及び再生の少なくとも一方の予測ができ、脱硝システムの的確な運転に寄与することができる。例えば、予測結果に応じて脱硝システムの各種の操作量を変更し、操作量の変更に伴ってステップＳ１０、Ｓ１２及びＳ２６における少なくとも１つのパラメータの設定・変更を行い、予測を繰り返すものとしてもよい。また、例えば、脱硝触媒の劣化レベルの表示や、運転条件を変更した場合の窒素酸化物やアンモニアの排出状態の予測、脱硝触媒の交換や再生のタイミングの告知、再生温度と再生時間の予測等が容易となる。また、劣化再生プログラムの出力に基づいて脱硝システムの操作量を変更した場合は、状況に応じ脱硝システムを自動的に運転することも可能となる。
【産業上の利用可能性】
【００７８】
以上の実施形態は、舶用ディーゼルエンジンのみならず、他の内燃機関や燃焼機器の脱硝触媒システムにも適用することができる。
【特許請求の範囲】
【請求項１】
コンピュータによって、燃焼機器からの排ガスにアンモニア系還元剤を添加して脱硝触媒で窒素酸化物の分解を行う脱硝システムにおける前記脱硝触媒の劣化と再生を予測するプログラムであって、
前記アンモニア系還元剤のパラメータ設定を行うステップ１と、
前記脱硝触媒に関する条件設定を行うステップ２と、
前記脱硝触媒でのアンモニアの反応量及び窒素酸化物の反応量を算出するステップ３と、
前記脱硝触媒での三酸化硫黄の生成量を算出するステップ４と、
前記脱硝触媒での酸性硫安の増減量を算出するステップ５と、
前記脱硝触媒の性能劣化量を算出するステップ６と、
算出結果を出力するステップ７と、
を備えることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項２】
請求項１に記載の脱硝触媒の劣化再生予測プログラムであって、
前記アンモニア系還元剤は、尿素水であり、
前記ステップ１における前記アンモニア系還元剤のパラメータは、前記尿素水の分解したアンモニアの濃度を含むことを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項３】
請求項２に記載の脱硝触媒の劣化再生予測プログラムであって、
前記ステップ２における条件は、前記脱硝触媒の温度及び体積と、前記排ガスの成分、流量及び温度と、であることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項４】
請求項１〜３のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記ステップ３は、Ｒｉｄｅａｌ−Ｅｌｅｙ機構を用いた脱硝反応速度式に基づいたアンモニアと窒素酸化物との反応量を算出する脱硝モデルを用いることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項５】
請求項１〜４のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記ステップ４は、非定常かつ、排ガス成分の影響を考慮した速度式に基づいた三酸化硫黄の生成量を算出する二酸化硫黄酸化反応モデルを用いることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項６】
請求項１〜５のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記ステップ５は、アンモニアと三酸化硫黄と水蒸気より酸性硫安を生じる平衡式の平衡定数を雰囲気ガス温度並びにアンモニア、三酸化硫黄及び水蒸気の分圧に基づいて定めた酸性硫安の増減量を算出する酸性硫安析出分解モデルを用いることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項７】
請求項６に記載の脱硝触媒の劣化再生予測プログラムであって、
前記酸性硫安析出分解モデルは、ケルビン方程式で補正されたものであることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項８】
請求項１〜７のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記ステップ６は、酸性硫安の増減量と劣化関数との関係を実験結果から求め、以下の式で算出される反応速度定数に基づいて行われることを特徴とする脱硝触媒の劣化再生予測プログラム。
ｋｃ＝ＡＮＯｅｘｐ（−ＥＮＯ／ＲＴｃａｔ）×Φ
ｋｃ：反応速度定数［ｍ３Ｎ／ｍ２ｈ］
ＡＮＯ、ＥＮＯ：ＮＯの反応に関する頻度因子［ｍ３Ｎ／ｍ２ｈ］
Ｒ：一般ガス定数［Ｊ／ｍｏｌＫ］
Ｔｃａｔ：触媒表面温度［Ｋ］
Φ：劣化関数
【請求項９】
請求項１〜８のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記脱硝触媒の出口におけるアンモニア量と窒素酸化物量を算出するステップ８をさらに備えることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項１０】
請求項１〜９のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記脱硝触媒の温度が所定の温度に達したか否かを判定するステップ９をさらに備えることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項１１】
請求項１〜１０のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記所定の温度を更新するステップ１０をさらに備えることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項１２】
請求項１〜１１のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記脱硝システムの運転時間が所定の時間に達したか否かを判断するステップ１１をさらに備えることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項１３】
請求項１〜１２のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記ステップ７は、時間と触媒性能劣化量の関係、時間と前記脱硝触媒の出口におけるアンモニア量の関係の少なくとも一方を出力することを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項１４】
燃焼機器からの排ガスにアンモニア系還元剤を添加して脱硝触媒で窒素酸化物の分解を行う脱硝システムにおいて、請求項１〜１３のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムを用いて前記脱硝触媒の劣化及び再生の少なくとも一方の予測を行うことを特徴とする脱硝システム。
【請求項１５】
請求項１４に記載の脱硝システムであって、
前記脱硝触媒の劣化再生予測プログラムの出力に基づいて、前記脱硝システムの操作量を変更し、前記ステップ１、２及び１０の少なくとも１つの設定を変更することを特徴とする脱硝システム。
【請求項１】
コンピュータによって、燃焼機器からの排ガスにアンモニア系還元剤を添加して脱硝触媒で窒素酸化物の分解を行う脱硝システムにおける前記脱硝触媒の劣化と再生を予測するプログラムであって、
前記アンモニア系還元剤のパラメータ設定を行うステップ１と、
前記脱硝触媒に関する条件設定を行うステップ２と、
前記脱硝触媒でのアンモニアの反応量及び窒素酸化物の反応量を算出するステップ３と、
前記脱硝触媒での三酸化硫黄の生成量を算出するステップ４と、
前記脱硝触媒での酸性硫安の増減量を算出するステップ５と、
前記脱硝触媒の性能劣化量を算出するステップ６と、
算出結果を出力するステップ７と、
を備えることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項２】
請求項１に記載の脱硝触媒の劣化再生予測プログラムであって、
前記アンモニア系還元剤は、尿素水であり、
前記ステップ１における前記アンモニア系還元剤のパラメータは、前記尿素水の分解したアンモニアの濃度を含むことを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項３】
請求項２に記載の脱硝触媒の劣化再生予測プログラムであって、
前記ステップ２における条件は、前記脱硝触媒の温度及び体積と、前記排ガスの成分、流量及び温度と、であることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項４】
請求項１〜３のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記ステップ３は、Ｒｉｄｅａｌ−Ｅｌｅｙ機構を用いた脱硝反応速度式に基づいたアンモニアと窒素酸化物との反応量を算出する脱硝モデルを用いることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項５】
請求項１〜４のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記ステップ４は、非定常かつ、排ガス成分の影響を考慮した速度式に基づいた三酸化硫黄の生成量を算出する二酸化硫黄酸化反応モデルを用いることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項６】
請求項１〜５のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記ステップ５は、アンモニアと三酸化硫黄と水蒸気より酸性硫安を生じる平衡式の平衡定数を雰囲気ガス温度並びにアンモニア、三酸化硫黄及び水蒸気の分圧に基づいて定めた酸性硫安の増減量を算出する酸性硫安析出分解モデルを用いることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項７】
請求項６に記載の脱硝触媒の劣化再生予測プログラムであって、
前記酸性硫安析出分解モデルは、ケルビン方程式で補正されたものであることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項８】
請求項１〜７のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記ステップ６は、酸性硫安の増減量と劣化関数との関係を実験結果から求め、以下の式で算出される反応速度定数に基づいて行われることを特徴とする脱硝触媒の劣化再生予測プログラム。
ｋｃ＝ＡＮＯｅｘｐ（−ＥＮＯ／ＲＴｃａｔ）×Φ
ｋｃ：反応速度定数［ｍ３Ｎ／ｍ２ｈ］
ＡＮＯ、ＥＮＯ：ＮＯの反応に関する頻度因子［ｍ３Ｎ／ｍ２ｈ］
Ｒ：一般ガス定数［Ｊ／ｍｏｌＫ］
Ｔｃａｔ：触媒表面温度［Ｋ］
Φ：劣化関数
【請求項９】
請求項１〜８のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記脱硝触媒の出口におけるアンモニア量と窒素酸化物量を算出するステップ８をさらに備えることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項１０】
請求項１〜９のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記脱硝触媒の温度が所定の温度に達したか否かを判定するステップ９をさらに備えることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項１１】
請求項１〜１０のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記所定の温度を更新するステップ１０をさらに備えることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項１２】
請求項１〜１１のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記脱硝システムの運転時間が所定の時間に達したか否かを判断するステップ１１をさらに備えることを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項１３】
請求項１〜１２のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムであって、
前記ステップ７は、時間と触媒性能劣化量の関係、時間と前記脱硝触媒の出口におけるアンモニア量の関係の少なくとも一方を出力することを特徴とする脱硝触媒の劣化再生予測プログラム。
【請求項１４】
燃焼機器からの排ガスにアンモニア系還元剤を添加して脱硝触媒で窒素酸化物の分解を行う脱硝システムにおいて、請求項１〜１３のいずれか１項に記載の脱硝触媒の劣化再生予測プログラムを用いて前記脱硝触媒の劣化及び再生の少なくとも一方の予測を行うことを特徴とする脱硝システム。
【請求項１５】
請求項１４に記載の脱硝システムであって、
前記脱硝触媒の劣化再生予測プログラムの出力に基づいて、前記脱硝システムの操作量を変更し、前記ステップ１、２及び１０の少なくとも１つの設定を変更することを特徴とする脱硝システム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【公開番号】特開２０１３−１５０２８（Ｐ２０１３−１５０２８Ａ）
【公開日】平成２５年１月２４日（２０１３．１．２４）
【国際特許分類】
【出願番号】特願２０１１−１４６４０１（Ｐ２０１１−１４６４０１）
【出願日】平成２３年６月３０日（２０１１．６．３０）
【新規性喪失の例外の表示】特許法第３０条第１項適用申請有り 『マリンエンジニアリング 第４６巻 第１号』 ２０１１年１月１日 社団法人 日本マリンエンジニアリング学会 発行
【出願人】（５０１２０４５２５）独立行政法人海上技術安全研究所 (185)
【出願人】（５９５０５５１６２）社団法人日本舶用工業会 (25)
【出願人】（５０３１１６８９９）新潟原動機株式会社 (61)
【出願人】（３９００３３０４２）ダイハツディーゼル株式会社 (43)
【出願人】（０００００５９０２）三井造船株式会社 (1,723)
【Ｆターム（参考）】
【公開日】平成２５年１月２４日（２０１３．１．２４）
【国際特許分類】
【出願日】平成２３年６月３０日（２０１１．６．３０）
【新規性喪失の例外の表示】特許法第３０条第１項適用申請有り 『マリンエンジニアリング 第４６巻 第１号』 ２０１１年１月１日 社団法人 日本マリンエンジニアリング学会 発行
【出願人】（５０１２０４５２５）独立行政法人海上技術安全研究所 (185)
【出願人】（５９５０５５１６２）社団法人日本舶用工業会 (25)
【出願人】（５０３１１６８９９）新潟原動機株式会社 (61)
【出願人】（３９００３３０４２）ダイハツディーゼル株式会社 (43)
【出願人】（０００００５９０２）三井造船株式会社 (1,723)
【Ｆターム（参考）】
[ Back to top ]