"複数の気筒を備える内燃機関である多気筒内燃機関の排気通路に配置されたガスセンサの出力値に基づき、前記多気筒内燃機関の初期異常を検知する多気筒内燃機関の初期異常検知方法であって、前記複数の気筒間の空燃比のバラツキの度合いである気筒間インバランス割合が異なる複数の条件のそれぞれで、前記ガスセンサの出力値の振幅を取得する振幅取得工程と、前記振幅取得工程において取得された複数の前記振幅に基づき、前記気筒間インバランス割合の変化に対する前記ガスセンサの前記出力値の前記振幅の変化の割合を算出する割合算出工程と、前記複数の振幅のうちの少なくとも1つの振幅と、前記割合算出工程において算出された前記変化の割合とに基づき、前記気筒間インバランス割合が、基準となる閾値である基準閾値よりも大きいと判断される場合を、前記多気筒内燃機関に初期異常がある場合として検知する第1異常検知工程とを備えることを特徴とする多気筒内燃機関の初期異常検知方法。前記複数の条件は、前記多気筒内燃機関への燃料の供給条件が異なる第1条件と、第2条件とを少なくとも含み、前記第1条件は、前記多気筒内燃機関の前記空燃比が予め設定された目標値となるように前記多気筒内燃機関の運転条件が制御される場合の前記燃料の前記供給条件であり、前記第2条件は、前記複数の気筒のうちの特定の1つの気筒である特定気筒に対しては、前記第1条件で前記多気筒内燃機関が制御された場合の前記燃料の供給量を特定割合変化させた量の前記燃料を供給する条件であり、且つ、前記複数の気筒のうちの前記特定気筒以外の気筒に対しては、前記第1条件で前記多気筒内燃機関が制御された場合に供給される量の前記燃料を供給する条件であることを特徴とする請求項1に記載の多気筒内燃機関の初期異常検知方法。前記第1異常検知工程は、記憶手段に記憶された、前記変化の割合の基準値である基準割合を取得する基準割合取得工程と、前記基準割合を前記変化の割合で除した値を前記ガスセンサの出力値の補正量として算出する補正量算出工程と、前記補正量を用いて前記ガスセンサの前記出力値を補正した場合の、前記少なくとも1つの振幅を算出する振幅補正工程と、前記振幅補正工程において算出された前記少なくとも1つの振幅に基づき、前記気筒間インバランス割合が前記基準閾値よりも大きいと判断される場合を、前記多気筒内燃機関に初期異常がある場合として検知する第2異常検知工程とを備えることを特徴とする請求項1又は2に記載の多気筒内燃機関の初期異常検知方法。前記第1異常検知工程は、記憶手段に記憶された、前記変化の割合の基準値である基準割合を取得する基準割合取得工程と、前記基準閾値及び前記基準閾値に対応する閾値であって、前記ガスセンサの前記出力値の振幅に関する閾値である対応閾値の少なくともいずれかを、前記変化の割合に基づき補正することによって補正閾値を算出する閾値補正工程と、前記少なくとも1つの振幅と、前記補正閾値とに基づき、前記気筒間インバランス割合が前記基準閾値よりも大きいと判断される場合を、前記多気筒内燃機関に初期異常がある場合として検知する第2異常検知工程とを備えることを特徴とする請求項1又は2に記載の多気筒内燃機関の初期異常検知方法。一般に、同じ構成を有するガスセンサ換言すれば、同一品番のガスセンサであっても、製造誤差及び排気管に対するガスセンサの装着方向等に起因して、ガスセンサの出力値にはバラツキが生じる。一方、同じ構成を有する多気筒内燃機関であっても、製造誤差に起因して、気筒間インバランス割合にはバラツキがある。多気筒内燃機関の初期異常を、気筒間インバランス割合が異常であるか否かの観点から検知する場合には、以下のような問題がある。すなわち、気筒間インバランス割合の検知に用いられるガスセンサのバラツキが大きいほど、気筒間インバランス割合が過大になっているか否かを判断する処理に用いる閾値を設定することが困難である。閾値が適切に設定されない場合、気筒間インバランス割合の異常を精度よく検知できない。本発明は、上記課題を解決するためになされたものであり、ガスセンサの出力値のバラツキを考慮しつつ、気筒間インバランス割合が異常であるか否かの観点から、多気筒内燃機関の初期異常を検知する多気筒内燃機関の初期異常検知方法を提供することを目的とする。上記課題を解決するために、本態様の多気筒内燃機関の初期異常検知方法は、複数の気筒を備える内燃機関である多気筒内燃機関の排気通路に配置されたガスセンサの出力値に基づき、前記多気筒内燃機関の初期異常を検知する多気筒内燃機関の初期異常検知方法であって、前記複数の気筒間の空燃比のバラツキの度合いである気筒間インバランス割合が異なる複数の条件のそれぞれで、前記ガスセンサの出力値の振幅を取得する振幅取得工程と、前記振幅取得工程において取得された複数の前記振幅に基づき、前記気筒間インバランス割合の変化に対する前記ガスセンサの前記出力値の前記振幅の変化の割合を算出する割合算出工程と、前記複数の振幅のうちの少なくとも1つの振幅と、前記割合算出工程において算出された前記変化の割合とに基づき、前記気筒間インバランス割合が、基準となる閾値である基準閾値よりも大きいと判断される場合を、前記多気筒内燃機関に初期異常がある場合として検知する第1異常検知工程とを備えている。インバランス割合の変化に対するガスセンサの出力値の振幅の変化は、原点を通る1次直線に近似可能である。したがって、ガスセンサの出力値の振幅から、気筒間インバランス割合を推定可能である。しかし、インバランス割合の変化に対するガスセンサの出力値の振幅の変化の割合は、同一の構成を有する複数のガスセンサ間でバラツキがある。これに対し、本態様の多気筒内燃機関の初期異常検知方法では、ガスセンサの出力値の振幅に含まれる、複数のガスセンサ間の出力値のバラツキの影響を、ガスセンサ個体の変化の割合を用いて低減させる。ガスセンサ個体の変化の割合とは、ガスセンサと、内燃機関とを用いて実測したインバランス割合の特定量変化に対するガスセンサの出力値の振幅の変化の割合である。したがって、本態様の多気筒内燃機関の初期異常検知方法によれば、ガスセンサの出力値のバラツキを考慮しつつ、気筒間インバランス割合が異常であるか否かの観点から、多気筒内燃機関の初期異常を検知することができる。本態様の多気筒内燃機関の初期異常検知方法において、前記複数の条件は、前記多気筒内燃機関への燃料の供給条件が異なる第1条件と、第2条件とを少なくとも含み、前記第1条件は、前記多気筒内燃機関の前記空燃比が予め設定された目標値となるように前記多気筒内燃機関の運転条件が制御される場合の前記燃料の前記供給条件であり、前記第2条件は、前記複数の気筒のうちの特定の1つの気筒である特定気筒に対しては、前記第1条件で前記多気筒内燃機関が制御された場合の前記燃料の供給量を特定割合変化させた量の前記燃料を供給する条件であり、且つ、前記複数の気筒のうちの前記特定気筒以外の気筒に対しては、前記第1条件で前記多気筒内燃機関が制御された場合に供給される量の前記燃料を供給する条件であってもよい。気筒間インバランス割合が異常であるか否かは、多気筒内燃機関の通常の運転条件を含む条件で検討されることが好ましい。多気筒内燃機関の通常の運転条件は、一般に、多気筒内燃機関の空燃比が予め設定された目標値となるように制御される条件、すなわち、第1条件である。したがって、この場合の多気筒内燃機関の初期異常検知方法は、第1条件と、第2条件とを少なくとも含む複数の条件に基づき算出された、ガスセンサ個体の変化の割合に基づき、気筒間インバランス割合が異常であるか否かを適切に検知することができる。本態様の多気筒内燃機関の初期異常検知方法において、前記第1異常検知工程は、記憶手段に記憶された、前記変化の割合の基準値である基準割合を取得する基準割合取得工程と、前記基準割合を前記変化の割合で除した値を前記ガスセンサの出力値の補正量として算出する補正量算出工程と、前記補正量を用いて前記ガスセンサの前記出力値を補正した場合の、前記少なくとも1つの振幅を算出する振幅補正工程と、前記振幅補正工程において算出された前記少なくとも1つの振幅に基づき、前記気筒間インバランス割合が前記基準閾値よりも大きいと判断される場合を、前記多気筒内燃機関に初期異常がある場合として検知する第2異常検知工程とを備えてもよい。ガスセンサの出力値が補正量算出工程で算出した補正量を用いて補正された場合、ガスセンサの出力値が補正されない場合に比べ、同一の構成を有する複数のガスセンサ間の出力値のバラツキが低減される。したがって、この場合の多気筒内燃機関の初期異常検知方法では、ガスセンサの出力値が補正されない場合に比べ、気筒間インバランス割合が過大になっているかを判断する処理に用いる閾値を設定することが容易である。したがって、多気筒内燃機関の初期異常検知方法は、補正量を用いて補正しない場合の振幅に基づき多気筒内燃機関の初期異常を検知する場合に比べ、多気筒内燃機関の初期異常を検知する精度を向上させることができる。本態様の多気筒内燃機関の初期異常検知方法において、前記第1異常検知工程は、記憶手段に記憶された、前記変化の割合の基準値である基準割合を取得する基準割合取得工程と、前記基準閾値及び前記基準閾値に対応する閾値であって、前記ガスセンサの前記出力値の振幅に関する閾値である対応閾値の少なくともいずれかを、前記変化の割合に基づき補正することによって補正閾値を算出する閾値補正工程と、前記少なくとも1つの振幅と、前記補正閾値とに基づき、前記気筒間インバランス割合が前記基準閾値よりも大きいと判断される場合を、前記多気筒内燃機関に初期異常がある場合として検知する第2異常検知工程とを備えてもよい。同一の構成を有する複数のガスセンサ間に出力値のバラツキがあったとしても、個々のガスセンサの出力値に応じた補正閾値が、ガスセンサ個体の変化の割合に基づき設定される。したがって、この場合の多気筒内燃機関の初期異常検知方法は、基準閾値に基づき多気筒内燃機関の初期異常を検知する場合に比べ、多気筒内燃機関の初期異常を検知する精度を向上させることができる。以下、本発明を具体化したセンサ制御装置の第1及び第2の実施形態について、図面を参照して順に説明する。参照する図面は、本発明が採用し得る技術的特徴を説明するために用いるものであり、記載している装置の構成及びフローチャート等は、それのみに限定する趣旨ではなく、単なる説明例である。以下の説明において、図2の上下方向、左右方向をそれぞれ、ガスセンサ1の左右方向、上下方向として説明する。第1及び第2の実施形態のセンサ制御装置4は、物理的構成と電気的構成が同じである。まず、図1を参照して、多気筒内燃機関以下、単に内燃機関という。100と、内燃機関100に取り付けられるセンサ制御装置4との概略的な構成について説明する。内燃機関100は、直列に配置された4つの気筒111から114を備え、ガソリンを燃料として用いて駆動するエンジンである。内燃機関100には、内燃機関100から排出される排気ガスを車外に放出するための排気管102が接続されている。センサ制御装置4は、ガスセンサ1と、ECU3とを備える。ガスセンサ1は、排気管102の経路上に取り付けられる。ガスセンサ1は、排気管102に流通するガスを検知対象ガスとする全領域空燃比センサ酸素センサである。ECU3は、ガスセンサ1とは離れた位置に配置されており、バッテリ80から電力の供給を受けて駆動する。ガスセンサ1と、ECU3とは、ハーネス91信号線を介して電気的に接続されている。ガスセンサ1は、ECU3によって通電制御されるとともに、ECU3に対して検知対象ガス中の酸素濃度に応じた電気信号を出力する。ECU3は、ガスセンサ1の出力に基づき空燃比を算出し、内燃機関100の空燃比のフィードバック制御を実行する。また、ECU3は、内燃機関100を駆動するための各種制御を実行する。次に、図2を参照して、ガスセンサ1について説明する。ガスセンサ1は、検知素子10と、ヒータ素子40と、ハウジング図示省略とを備える。検知素子10は、固体電解質体11,13と、絶縁基体12,24とを、固体電解質体13,絶縁基体12,固体電解質体11,絶縁基体24の順に積層した構造を有する。固体電解質体11,13と、絶縁基体12,24と、後述する絶縁基体17,18とは、いずれも細長い板状に形成されており、図2ではガスセンサ1の長手方向と直交する断面を模式的に示している。固体電解質体11,13は、イットリアを安定化剤に用いた部分安定化ジルコニアを主成分とする材料によって形成され、酸素イオン伝導性を有する。絶縁基体12,24は、アルミナを主成分とする材料によって形成される。ヒータ素子40は、固体電解質体11,13の早期活性化と、固体電解質体11,13の活性の安定性維持とのために、固体電解質体13に積層されている。ハウジングは、検知素子10と、ヒータ素子40とを内部に保持し、排気管102図1参照に取り付けられる。以下、ガスセンサ1が備える検知素子10と、ヒータ素子40とについて詳述する。まず、図2を参照して、検知素子10の構成を説明する。検知素子10は、検知室23と、拡散律速部15と、酸素ポンプセル27以下、Ipセル27という。と、酸素分圧検知セル28以下、Vsセル28という。と、絶縁基体12,24とを備える。検知室23は、排気管102図1参照内を流通する排気ガスが導入される小空間である。検知室23は、固体電解質体11と、固体電解質体13と、拡散律速部15と、絶縁基体12とによって囲まれている。拡散律速部15は、検知室23の幅方向図1の紙面左右方向の両端にある。拡散律速部15は、多孔質状の部材例えば、アルミナによって形成され、検知室23内に検知対象ガスを導入する際の流入量を規制する。Ip1セル27は、固体電解質体11と、多孔質性の電極19,20とを備える。電極19,20は、Ptを主成分とする材料によって形成される。Ptを主成分とする材料としては、例えば、Ptと、Pt合金と、Pt及びセラミックスを含むサーメットとが挙げられる。電極20は、固体電解質体11の面のうち、検知室23側の面に設けられている。電極19は、固体電解質体11の面のうち、検知室23側の面とは反対側の面に設けられている。すなわち、検知素子10の積層方向において、一対の電極19,20は、固体電解質体11を挟むように配置されている。固体電解質体11の上面には、絶縁基体24が積層されている。絶縁基体24は、電極19の上部となる位置に開口29を有し、開口29には保護層25が設けられている。保護層25は、セラミックス例えば、アルミナからなる多孔質性の部材によって形成され、電極19が検知対象ガスに含まれるシリコン等の被毒成分によって劣化しないように、電極19の上面を覆っている。Ip1セル27は、電極19,20間に電流が供給されることで、電極19が接する雰囲気検知素子10の外部の雰囲気と電極20が接する雰囲気検知室23内の雰囲気との間で、酸素の汲み出し及び汲み入れいわゆる酸素ポンピングを行う。Vsセル28は、固体電解質体13と、多孔質性の電極21,22とを備える。固体電解質体13は、検知室23を挟んで固体電解質体11と対向するように配置されている。電極21は、固体電解質体13の面のうち、検知室23側の面に設けられている。電極22は、固体電解質体13の面のうち、検知室23側の面とは反対側の面に形成されている。すなわち、検知素子10の積層方向において、一対の電極21,22は、固体電解質体13を挟むように配置されている。電極21,22は、上述のPtを主成分とする材料によって形成される。Vsセル28は、主として、固体電解質体13によって隔てられた2つの雰囲気電極21と接する検知室23内の雰囲気と、電極22と接する雰囲気間の酸素分圧差に応じて起電力を発生する。また、電極22は、絶縁基体17によって排気管102図1参照を流通するガスと接触しないように閉塞されている。そして、電極22は、詳細は後述するが、検知室23内の酸素濃度の検知のための基準となる酸素濃度を維持する酸素基準電極として機能する。次に、図2を参照して、ヒータ素子40の構成を説明する。ヒータ素子40は、検知素子10特に、固体電解質体11,13を加熱して活性化させる。ヒータ素子40は、発熱抵抗体26と、絶縁基体17,18とを備える。発熱抵抗体26は、白金を主成分とする材料によって構成され、絶縁基体17と絶縁基体18との間に挟まれている。絶縁基体17,18は、アルミナを主成分とする材料によって形成される。次に、図2を参照して、ECU3の構成について説明する。ECU3は、自動車に搭載された内燃機関100と、ガスセンサ1とを制御する。図2のように、ECU3は、マイクロコンピュータ9と、センサ制御回路としてのASIC特定用途向け集積回路30と、ヒータ電圧供給回路50とを構成主体としている。マイクロコンピュータ9は、公知の構成のCPU6と、ROM7と、RAM8と、EEPROM5とを搭載したマイコンチップである。ROM7には、CPU6が実行する各種制御プログラムと、制御プログラム実行時に参照される各種パラメータとが記憶されている。図2のように、ASIC30は、Ip検知回路31と、Ip駆動回路32と、抵抗検知回路33と、電圧出力回路34と、基準電圧比較回路35と、Icp供給回路36とを備える。Ip検知回路31は、Ipセル27の電極19,20間に流れる電流Ipを電圧変換し、変換した電圧を検知信号としてマイクロコンピュータ9に出力する。抵抗検知回路33は、定期的に、予め規定された値の電流をVsセル28に通電し、その通電に応答して得られる電圧変化量電圧Vsの変化量を検知するための回路である。抵抗検知回路33によって検知された電圧Vsの変化量を示す値は、マイクロコンピュータ9に出力される。マイクロコンピュータ9では、抵抗検知回路33から出力された値と、ROM7に記憶されている電圧Vsの変化量とVsセル28のインピーダンスRiとが予め関連付けられたテーブルとに基づいて、Vsセル28のインピーダンスRiが求められる。Vsセル28のインピーダンスRiは、Vsセル28の温度、すなわち、検知素子10全体の温度と相関があり、マイクロコンピュータ9は、Vsセル28のインピーダンスRiに基づいて、ガスセンサ1検知素子10の温度を検知する。電圧出力回路34は、Vsセル28の電極21,22間に生ずる起電力Vsを検知する。基準電圧比較回路35は、予め定められた基準電圧と、電圧出力回路34において検知された起電力Vsとの比較を行い、比較結果をIp駆動回路32に出力する。Ip駆動回路32は、基準電圧比較回路35から出力された比較結果に基づき、Ipセル27の電極19,20間に供給する電流Ipの大きさ及び向きを制御する。Icp供給回路36は、Vsセル28の電極22から電極21へ向かって流れる微小電流Icpを供給する。ヒータ電圧供給回路50は、CPU6から指示に応じて、発熱抵抗体26の両端に印加される電圧Vhを公知のPI制御のもと生成したり、発熱抵抗体26の両端に一定の電圧例えば、12Vを印加したりすることで、発熱抵抗体26を発熱させる。図2に示すように、ECU3には、内燃機関100を制御するための各種センサと、各種装置とが電気的に接続されている。各種センサには、例えば、内燃機関100のクランク角クランク軸の回転状態であるクランク角を検知するクランク角度センサ104と、冷却水の温度を検知する水温センサ106とがある。各種装置には、例えば、ECU3から指示された量の燃料を、ECU3から指示された気筒に供給する燃料噴射装置103がある。ECU3には、さらに、内燃機関100の初期異常を報知する故障警告灯105が電気的に接続されている。次に、図2を参照して、ガスセンサ1を用いて検知対象ガスの酸素濃度排気ガスの空燃比を検知する動作について簡単に説明する。なお、検知対象ガスの酸素濃度を検知する際には、基準電圧比較回路35で比較対象となる基準電圧例えば450mVが設定される。まず、Icp供給回路36は、Vsセル28の電極22から固体電解質体13を介して電極21に向けて微小電流Icpを供給する。この通電により、検知対象ガス中の酸素は、酸素イオンとなって、固体電解質体13を介して電極21側から電極22側に移動する汲み込まれる。Vsセル28に電流Icpが供給されると、酸素イオンは電極21側から電極22側に移動して電極22自身及び周囲に酸素が蓄積され、検知対象ガス中の酸素濃度に応じた起電力Vsを発生させるための基準となる酸素濃度雰囲気が生成される。電圧出力回路34は、両電極21,22間の起電力Vsを検知し、検知した起電力Vsを基準電圧比較回路35に出力する。基準電圧比較回路35は、起電力Vsと基準電圧とを比較し、比較結果をIp駆動回路32に出力する。Ip駆動回路32は、基準電圧比較回路35による比較結果に基づいて、起電力Vsが基準電圧となるように、Ipセル27の電極19,20間に供給する電流Ipの大きさ及び向きを制御する。これにより、Ipセル27による検知室23内への酸素の汲み入れ又は検知室23からの酸素の汲み出しが行われる。なお、検知室23内に流入した排気ガスの空燃比が理論空燃比よりもリッチであった場合、Ipセル27において外部から検知室23内に酸素を汲み入れるように、電極19,20間に供給される電流Ipの大きさや向きが制御される。一方、検知室23内に流入した排気ガスの空燃比が理論空燃比よりもリーンであった場合、Ipセル27において検知室23から外部へ酸素を汲み出すように、電極19,20間に供給される電流Ipの大きさ及び向きが制御される。このときの電流Ipは、Ip検知回路31において電圧変換され、電圧変換された信号は、検知信号としてマイクロコンピュータ9に出力される。マイクロコンピュータ9CPU6は、検知信号に基づいて検知対象ガス中に含まれる酸素濃度対応値、ひいては空燃比を算出する。次に、図3を参照して、第1の実施形態の内燃機関100の初期異常検知処理の原理について概説する。気筒間インバランス割合以下、単にインバランス割合ともいう。が増減する理由には、上述のように、内燃機関での燃焼に関わる種々の要因が挙げられる。本実施形態では特に、燃料供給量に着目し、インバランス割合を、燃料供給量ズレを起こしている1つの気筒の燃料供給量Qiと、燃料供給量ズレを起こしていないその他の気筒の燃料供給量Qsとの差を、燃料供給量Qsで除した値の百分率であると定義する。インバランス割合が0であり、且つ、ガスセンサ1の出力値のバラツキが0である場合、内燃機関100の1燃焼サイクル本実施形態では、クランク角720度サイクルでの空燃比の変動は0である。インバランス割合が大きくなると、インバランス割合が小さい場合に比べ、1燃焼サイクルでの空燃比の変動は大きくなる。すなわち、インバランス割合に対するガスセンサ1の出力値の振幅は、図3の直線202のように、原点を通る一次直線に近似される。同じ構造を有する複数のガスセンサ1であっても、製造誤差及び排気管102に対するガスセンサ1の装着方向等に起因して、ガスセンサ1の出力値にバラツキが生じることがある。例えば、同じ構造を有する複数のガスセンサ1間には、個体差があることがある。また、排気管102に対するガスセンサ1の装着方向が複数のガスセンサ1間で異なると、ガスセンサ1のハウジング図示省略の外部から検知素子10周囲に排気ガスが到達するのに要する時間排気ガスの到達速度も複数のガスセンサ1間で異なることがある。このため、複数のガスセンサ1間でガス応答性にバラツキが生じることがある。ガス応答性とは、検知対象ガス例えば、排気ガス中の特定ガス例えば、酸素の濃度が変化してから、実際にガスセンサによってその変化が検知されるまでの時間に関する指標である。ガス応答性が高いガスセンサほど、ガス応答性が低いガスセンサに比べ、上記時間が短い。同じインバランス割合を有する内燃機関100から排出される排気ガスを検知対象ガスとした場合、一般に、ガス応答性が高いガスセンサ1では、ガス応答性が低いガスセンサ1に比べ、ガスセンサ1の出力値の振幅は大きくなる。このため、ガスセンサ1個体の変化の割合を表す直線の傾きは、ガスセンサ1のガス応答性に応じて、複数のガスセンサ1間で異なる。具体的には、図3のように、ガス応答性が高いガスセンサ1に対応する直線201の傾きは、ガス応答性が低いガスセンサ1に対応する直線203の傾きに比べ急である。したがって、複数のガスセンサ1において、出力値の振幅が同じであったとしても、複数のガスセンサ1個体間の出力値のバラツキに起因して、各出力値の振幅から推定されるインバランス割合は異なることがある。一方、内燃機関100の空燃比が予め設定された目標値となるように制御される場合以下、通常制御時ともいう。のインバランス割合は、製造誤差等に起因して、0よりも大きい内燃機関100個体毎に異なる値を示す。しかし、ECU3は、通常制御時のインバランス割合の内燃機関100個体間のバラツキを認識できない。このため、以下のような問題が生じることがある。複数の内燃機関100において同じ目標値を用いて通常制御が実行された場合の平均的なインバランス割合がJであるとする。ECU3は、点204及び点206のように、通常制御時におけるインバランス割合がJからズレてしまっていたとしても、直接そのズレ量を認識することができない。ズレ量が許容量ではない場合は、内燃機関100に初期異常がある場合であり、内燃機関100を備える製品が出荷される前に検知されることが好ましい。そこで、本実施形態の初期異常検知処理では、上述のようなガスセンサ1個体間の出力値のバラツキと、内燃機関100のインバランス割合の初期バラツキとを考慮して、以下のように内燃機関100の初期異常を検知する。まず、本実施形態の初期異常検知処理では、ガスセンサ1個体の変化の割合と、ROM7に記憶された基準割合とに基づき、ガスセンサ1の出力値の補正量が算出される。より具体的には、ガスセンサ1個体の変化の割合は、インバランス割合の異なる2つの条件を意図的に作り出し、各条件で得られる出力値の振幅に基づき算出される。本実施形態の初期異常検知処理では、2つの条件として、インバランス割合が互いに異なる第1条件と、第2条件とが設定される。インバランス割合が異なる条件は、複数の気筒のうちの特定の気筒である特定気筒に供給される燃料の供給条件を意図的に変えることによって設定される。特定気筒は、複数の気筒のうちの予め任意に選択された1つの気筒であればよく、例えば、4つの気筒111から114のうち、気筒111が特定気筒として特定される。第1条件は、内燃機関100の空燃比が予め設定された目標値となるように内燃機関100の運転条件が制御される場合の燃料の供給条件である。第1条件では、通常、すべての気筒に対して同じ供給量の燃料が供給されるように燃料噴射装置103が制御されるが、実際には各気筒に供給される供給量にバラツキがあり得る。第2条件は、特定気筒に対しては、第1条件で内燃機関100が制御された場合の燃料の供給量を特定割合変化させた量の燃料を供給する条件であり、複数の気筒のうちの特定気筒以外の気筒に対しては、第1条件で内燃機関100が制御された場合に供給される量の燃料を供給する条件である。ROM7に記憶された基準割合は、基準となる変化の割合である。基準割合は、例えば、複数のガスセンサ1のそれぞれの変化の割合を考慮して予め設定されればよい。本実施形態の基準割合は、同じ構成を有する複数のガスセンサ1と、テスト用の内燃機関100とを用いて、各ガスセンサ1について算出した上記変化の割合の平均値である。基準割合は、例えば、図3の直線202の傾きで示される。ガスセンサ1の出力値の補正量は、ガスセンサ1の出力値のバラツキを低減させるために設定される。本実施形態の補正量は、基準割合をガスセンサ1個体の変化の割合で除した値である。本実施形態の初期異常検知処理では、ガスセンサ1の出力値と、補正量との積を算出することによって、ガスセンサ1の出力値を補正する。前述のように、インバランス割合に対するガスセンサ1の出力値の振幅は、原点を通る一次直線に近似可能であるので、補正後のガスセンサ1の出力値は、基準割合を有するガスセンサ1の出力値とほぼ等しい。すなわち、上述のように算出した補正量を用いてガスセンサ1の出力値が補正されることによって、補正前のガスセンサ1の出力値に比べ、ガスセンサ1の出力値のバラツキが低減される。次に、本実施形態の初期異常検知処理では、第1条件で内燃機関を運転させた場合のガスセンサ1の出力値を、補正量を用いて補正した場合の振幅が、線分208で表す閾値TH2よりも大きい場合を、インバランス割合に初期異常がある場合とする。第1条件は、上述の通常制御時の内燃機関100の燃料供給条件である。閾値TH2は、基準閾値TH1に対応する閾値であって、ガスセンサ1の出力値の振幅に関する閾値である。基準閾値TH1は、インバランス割合に異常があるか否かの基準となる、インバランス割合に関する閾値である。閾値TH2は、基準閾値TH1と、基準割合と、補正後のガスセンサ1の出力値のバラツキとを含む条件を考慮して設定される。本実施形態では、閾値TH2を表す線分208は、基準閾値TH1を表す線分209と、直線202との交点を通る。以下、閾値TH2を対応閾値TH2ともいう。本実施形態では、補正後の振幅が、対応閾値TH2よりも大きいか否かが判断されることによって、間接的に、通常制御時のインバランス割合が、基準閾値TH1よりも大きいか否かが判断される。次に、図3から図6を参照して、第1の実施形態の初期異常検知処理について説明する。第1の実施形態の初期異常検知処理は、例えば、内燃機関100と、センサ制御装置4とを備える自動車の製造検査時または出荷検査時に、作業者によって開始の指示が入力された場合に起動される。図4のフローチャートに示す各処理を実行させるプログラムは、図2のROM7に記憶されており、CPU6によって実行される。具体例として、インバランス割合に対するガスセンサ1の出力値の振幅が図3の直線201の関係となるガスセンサ1以下、ガスセンサAともいう。及び内燃機関100以下、内燃機関Aともいう。の組み合わせと、図3の直線203の関係となるガスセンサ1以下、ガスセンサBともいう。及び内燃機関100以下、内燃機関Bともいう。の組み合わせとのそれぞれで、初期異常検知処理が実行された場合について説明する。図4に示すように、初期異常検知処理が起動された場合、まず、作業者によって内燃機関100を始動させる指示が入力されたか否かが判断されるS5。内燃機関100を始動させる指示が入力されていない場合S5:NO、CPU6は指示が入力されるまで待機する。内燃機関100を始動させる指示が入力された場合S5:YES、内燃機関100及びガスセンサ1が始動されるS10。次に、作業者によって第1条件での内燃機関100の運転を開始させる指示が入力されたか否かが判断されるS15。第1条件での内燃機関100の運転を開始させる指示が入力されていない場合S15:NO、CPU6は指示が入力されるまで待機する。第1条件での内燃機関100の運転を開始させる指示が入力された場合S15:YES、内燃機関100の運転条件に第1条件が設定され、第1条件での内燃機関100の運転が開始されるS20。ステップS20では、例えば、燃料噴射装置103が制御され、内燃機関100が備える4つの気筒のそれぞれに第1条件に従って燃料が供給される。第1条件における、燃料供給条件以外の条件は、基準割合を求めた条件と同じ条件であることが好ましい。第1条件における、燃料供給条件以外の条件としては、例えば、内燃機関100の冷却水の温度が80°C以上になる条件が挙げられる。冷却水の温度は、水温センサ106によって検知される。次に、内燃機関100の冷却水の温度が80°C以上となり、且つ、ガスセンサ1が活性化して、内燃機関100の空燃比を検出可能な状態になった後、振幅取得処理が実行されるS30。ステップS30では、第1条件で内燃機関100が制御された場合の、ガスセンサ1の出力値の振幅AM1が算出される。図5を参照して、振幅取得処理の詳細を説明する。図5に示すように、振幅取得処理ではまず、所定長さのサンプリング期間中、所定のサンプリング周期でガスセンサ1の出力値と、内燃機関100のクランク角度とが取得され、取得された出力値は、出力値の取得時間及びクランク角度と対応付けられてRAM8に記憶されるS32。サンプリング期間及びサンプリング周期は適宜設定されればよい。本実施形態では、サンプリング期間を40secとし、サンプリング周期を1msecとする。クランク角度は、クランク角度センサ104から出力される。ガスセンサ1の出力値の取得時間は、少なくとも、ガスセンサ1の出力値の取得順序を示す。N番目Nは自然数に取得されたガスセンサ1の出力値をIpNとし、IpNと対応付けられるクランク角度をCANとする。サンプリング期間、サンプリング周期、及び出力値の取得時間は、例えば、図示しないタイマ回路によって管理される。ガスセンサ1の出力値IpNは、必要に応じてローパスフィルタLow-passfilterに供され、ノイズが除去される。次に、ガスセンサ1の出力値の移動平均が算出され、算出された移動平均はRAM8に記憶されるS34。ステップS34では、ガスセンサ1の出力値の移動平均として、順次、1燃焼サイクル単位期間に含まれるガスセンサ1の出力値を対象にして、当該単位期間の出力値の平均値が算出される。1燃焼サイクル単位期間にM個Mは自然数のガスセンサ1の出力値が含まれる場合、N番目の移動平均IpvNは、式1によって表される。IpvN=IpNIpN1IpN2IpNM-1/M式1IpvNは、クランク角度CANと対応付けられてRAM8に記憶される。次に、ステップS34で算出された移動平均を用いて、内燃機関100の燃焼要因のバラツキに起因すると推定される、ガスセンサ1の出力値のバラツキが補正されるS36。内燃機関100の燃焼要因のバラツキとしては、例えば、吸気装置によって吸入される空気量と、気筒間に分配される空気量と、燃料噴射装置103によって供給される燃料の量とのそれぞれの、燃焼サイクル毎のバラツキが挙げられる。ステップS36では、例えば、前述のようにN番目の移動平均IpvNが算出された場合、補正後のN番目のガスセンサ1の出力値CNは式2に従って算出される。CN=IpNL-IpvN式2Lは、Mよりも小さい自然数である。IpNLは、N番目のガスセンサ1の出力値IpNが取得されてから、1/2燃焼サイクル後に取得されるガスセンサ1の出力値である。すなわち、クランク角度CANと、クランク角度CANLとは、約360度異なる。補正後のN番目のガスセンサ1の出力値CNは、クランク角度CANLと対応付けられる。次に、補正後のN番目のガスセンサ1の出力値CNについて、燃焼サイクル毎に振幅が算出され、算出された振幅はRAM8に記憶されるS38。図6に示すように、ガスセンサ1の出力値の経時変化は、1燃焼サイクルを1周期として変動する波形を示す。ステップS38では、1燃焼サイクル毎に、出力値の最大値と最小値とが取得される。図6において、クランク角度が0度である場合を1燃焼サイクルの開始点とし、クランク角度が720度である場合を終了点とした場合、矢印231から233で示す期間のそれぞれが、1燃焼サイクルに対応する。矢印231に示す期間では、最大値として出力値211が、最小値として出力値221がそれぞれ取得される。矢印232に示す期間では、最大値として出力値212が、最小値として出力値222がそれぞれ取得される。矢印233に示す期間では、最大値として出力値213が、最小値として出力値223がそれぞれ取得される。以下、同様である。次に、同一の燃焼サイクルに含まれる最大値と最小値との差が燃焼サイクル毎のガスセンサ1の出力値の振幅として算出される。次に、ステップS38で算出された燃焼サイクル毎の出力値の振幅の平均値が算出され、算出された平均値は、現在の内燃機関100の運転条件第1条件に対応する振幅AMn図4のステップ30で実行される振幅処理では、n=1としてRAM8に記憶されるS40。ガスセンサ1がガスセンサAである場合、振幅AM1として図3の点204で示される振幅が取得され、ガスセンサ1がガスセンサBである場合、振幅AM1として点206で示される振幅が取得されたものとする。前述のように、同じ条件で内燃機関100が制御された場合にも、内燃機関100個体間で実際のインバランス割合が異なるため、この時点では、AM1に基づき第1条件に対応するインバランス割合を推定することはできない。振幅取得処理は以上で終了する。図4の初期異常検知処理の説明に戻り、ステップS30の次に、第1条件に対応する振幅が取得されたか否かが判断されるS45。第1条件に対応する振幅が取得されていない場合S45:NO、CPU6は、振幅が取得されるまで待機する。第1条件に対応する振幅が取得された場合S45:YES、作業者によって第2条件での内燃機関100の運転を開始させる指示が入力されたか否かが判断されるS48。第2条件での内燃機関100の運転を開始させる指示が入力されていない場合S48:NO、CPU6は指示が入力されるまで待機する。第2条件での内燃機関100の運転を開始させる指示が入力された場合S48:YES、内燃機関100の運転条件に第2条件が設定され、第2条件での内燃機関100の運転が開始されるS50。ステップS50では、特定気筒例えば、気筒111に対しては、第1条件における燃料の供給量以下、第1供給量FS1ともいう。を特定割合変化させた量以下、第2供給量FS2ともいう。の燃料を供給するための指示が燃料噴射装置103に出力される。特定割合の変化は、予め任意に設定された変化であればよく、例えば、10の増加である。複数の気筒のうちの特定気筒以外の気筒例えば、気筒112から114に対しては、第1条件における燃料の供給量を供給するための指示が燃料噴射装置103に出力される。第1供給量FS1及び第2供給量FS2は、内燃機関100の燃焼要因のバラツキに起因して、実際に気筒に供給される量とは等しいとは限らない。複数の気筒のそれぞれに同じ量の燃料を供給するための指示が出力されたとしても、実際に各気筒に供給される量が等しいとは限らない。第2条件における、燃料供給量以外の条件は、基準割合を求めた条件と同じ条件であることが好ましい。次に、振幅取得処理が実行され、第2条件で内燃機関100が制御された場合の、ガスセンサ1の出力値の振幅AM2が算出されるS60。ステップS60の処理は、基本的にステップS30と同様であるので説明を省略する。ガスセンサ1がガスセンサAである場合、第2条件に対応する振幅AM2として図3の点205で示される振幅が取得され、ガスセンサ1がガスセンサBである場合、振幅AM2として点207で示される振幅が取得されたものとする。この時点では、第2条件に対応するインバランス割合は推定できない。次に、インバランス割合の変化に対するガスセンサ1の出力値の振幅の変化の割合傾きが算出され、算出された変化の割合はRAM8及びEEPROM5に記憶されるS70。変化の割合は、振幅AM1と、振幅AM2とを用いた、式3に基づき算出される。変化の割合=AM2-AM1/第2条件でのインバランス割合-第1条件でのインバランス割合式3式3において、分母は、燃料の供給量の特定割合変化させた量に実質的に起因するインバランス割合の差に該当する。したがって、式3の分母は、特定割合K例えば、10に等しい。ガスセンサ1個体の変化の割合を算出することによって、ガスセンサ1に固有なインバランス割合に対するガスセンサ1の出力値の振幅を表す一次直線式が得られる。したがって、この時点では、ガスセンサ1の出力値の振幅と、直線式とに基づき、出力値の振幅に対応するインバランス割合を推定可能である。次に、ROM7に記憶されている基準割合が取得され、取得された基準割合はRAM8に記憶されるS80。本実施形態の基準割合は、上述のように複数のガスセンサ1を用いて予め設定され、ROM7に記憶されている。次に、基準割合を、ガスセンサ1個体の変化の割合で除した値が補正量として算出され、算出された補正量はRAM8に記憶されるS90。次に、第1条件で内燃機関100が制御された場合のガスセンサ1の出力値が補正量を用いて補正された場合の振幅AM1 ́が算出され、算出された振幅AM1 ́はRAM8に記憶されるS100。ステップS100では、振幅AM1と、補正量との積が振幅AM1 ́として算出されてもよい。また、ステップS100では、第1条件で内燃機関100が制御された場合のガスセンサ1の出力値と、補正量との積が算出された後、補正後の出力値の振幅が図5のステップS38と同様に算出されてもよい。図3において、ガスセンサ1がガスセンサAである場合には、点204で表される値が、例えば、点244で表される値に補正される。ガスセンサ1がガスセンサBである場合には、点206で表される値が、例えば、点246で表される値に補正される。次に、振幅AM1 ́が、対応閾値TH2よりも大きいか否かが判断されるS110。点246が表す振幅AM1 ́は、線分208で表される対応閾値TH2よりも大きいS110:YES。したがって、内燃機関Bに初期異常があることを報知するために、故障警告灯105が点灯された後S120、初期異常検知処理は終了する。作業者は、故障警告灯105を確認することによって、内燃機関Bの初期異常を認識可能である。点244が表す振幅AM1 ́は、対応閾値TH2以下であるS110:NO。この場合、内燃機関Aはインバランス割合の観点から正常であるため、故障警告灯105は点灯されず、初期異常検知処理は終了する。以上のように、CPU6は、第1の実施形態の初期異常検知処理を実行する。[評価試験]第1の実施形態の初期異常検知処理に従って算出される補正量を用いてガスセンサ1の出力値を補正した場合に、ガスセンサ1の出力値のバラツキが低減されるか否かを確認する評価試験を行った。評価試験方法は以下の通りである。上記の同じ構成を有する12個のガスセンサ1のそれぞれを、直列4気筒のテストエンジン図示省略の排気管図示省略の経路上に取りつけた。排気管には、テストエンジンの各気筒から排出される排気ガスが流通する。テストエンジンが備える4つの気筒をそれぞれ1から4で表す。テストエンジンの運転条件は、以下の条件とした。試験1では、テストエンジンの回転数を1000rpm、インテークマニフォールド図示省略の圧力を-40kPa、テストエンジンの冷却水の温度を80°C以上とした。試験2では、テストエンジンの回転数を1500rpm、インテークマニフォールドの圧力を-40kPa、テストエンジンの冷却水の温度を80°C以上とした。試験1及び2における各気筒に対する燃料の供給条件に関して、以下のように、第1条件及び第2条件を設定した。第1条件では、1で表される気筒に対しては、通常制御時の燃料供給条件における燃料の供給量の120の量を供給するための指示が燃料噴射装置103に出力される。また、第1条件では、2から4のそれぞれで表される気筒に対しては、通常制御時の燃料供給条件における燃料の供給量を供給するための指示が燃料噴射装置103に出力される。1で表される気筒と、2から4のそれぞれで表される気筒とで燃料の供給量を変えているのは、以下の理由による。すなわち、複数の気筒のそれぞれに対して、同じ第1供給量の燃料を供給する指示が出力されても、内燃機関100の製造誤差等に起因して、実際に各気筒に供給される燃料にはバラツキがあり得る。したがって、現実には第1条件に対応するインバランス割合は0ではないため、現実の条件に近づけるため、上述のように第1条件を設定した。一方、第2条件では、1で表される気筒に対しては、通常制御時の燃料供給条件における燃料の供給量の130の量を供給するための指示が燃料噴射装置103に出力される。また、第2条件では、2から4のそれぞれで表される気筒に対しては、通常制御時の燃料供給条件における燃料の供給量を供給するための指示が燃料噴射装置103に出力される。したがって、評価試験における特定割合は8.3である。試験1及び2の各条件で、第1の実施形態の初期異常検知処理を行い、12個のガスセンサ1のそれぞれについてガスセンサ1個体の変化の割合を算出した。12個の変化の割合の平均値を基準割合として求めた。補正量を用いて補正した場合のガスセンサ1の出力値の振幅に基づき、インバランス割合に対するガスセンサ1の出力値の振幅を再計算した。試験1の結果を図7に示し、試験2の結果を図8に示す。図7及び図8において、補正前のインバランス割合に対するガスセンサ1の出力値の振幅を、黒い四角及び実線で示し、補正後のインバランス割合に対するガスセンサ1の出力値の振幅を、白丸及び破線で示す。図7及び図8に示すように、テストエンジンの回転数が異なる試験1及び試験2の条件ではいずれも、直線の分布範囲は、実線で示す補正前に対応する直線に比べ、破線で示す補正後に対応する直線の方が狭かった。より具体的には、試験1及び試験2において、第1条件における振幅の最大値と平均値との差を平均値で除した値の百分率Δと、第1条件における振幅の最小値と平均値との差を平均値で除した値の百分率-Δとはそれぞれ次の値であった。試験1では、補正前のΔは22.8であり、-Δは-21.0であったのに対し、補正後のΔは5.3であり、-Δは-6.5であった。試験2では、補正前のΔは28.1であり、-Δは-28.6であったのに対し、補正後のΔは12.6であり、-Δは-13.1であった。以上より、第1の実施形態の初期異常検知処理に従って算出される補正量でガスセンサ1の出力値を補正した場合に、ガスセンサ1の出力値のバラツキが低減されることが確認された。上記初期異常検知処理において、図4のステップS30及びステップS60のそれぞれは、本発明の振幅取得工程に相当する。ステップS70は、本発明の割合算出工程に相当する。ROM7は、本発明の記憶手段に相当する。ステップS80は、本発明の基準割合取得工程に相当する。ステップS90は、本発明の補正量算出工程に相当する。ステップS100は、本発明の振幅補正工程に相当する。ステップS110と、ステップS120とは、本発明の第2異常検知工程に相当する。ステップS80と、ステップS90と、ステップS100と、ステップS110と、ステップS120とは、本発明の第1異常検知工程に相当する。第1の実施形態の初期異常検知処理によれば、下記の効果が得られる。ガスセンサの出力値の振幅に含まれる、複数のガスセンサ間の出力値のバラツキの影響を、基準割合と、ガスセンサ個体の変化の割合とを用いて低減させる。具体的には、ガスセンサ1の出力値をステップS90で算出した補正量を用いて補正した後のガスセンサの出力値は、基準割合を示すガスセンサの出力値とほぼ一致する。例えば、図3において、点204で示されるガスセンサAの出力値の振幅AM1と、点206で示されるガスセンサBの出力値の振幅AM1は、ほぼ同じ値である。内燃機関Aの第1条件におけるインバランス割合はTH1よりも小さい正常範囲の値であり、内燃機関Bの第1条件におけるインバランス割合はTH1よりも大きい異常範囲の値であるとする。振幅AM1に基づき内燃機関100に初期異常があるか否かを判断する場合、ガスセンサAの出力値の振幅AM1と、ガスセンサBの出力値の振幅AM1との僅かな差を区別する閾値を設定しなければならない。これに対し、振幅AM1 ́に基づき内燃機関100が正常であるか否かを判断する場合、振幅AM1に基づき内燃機関100が正常であるか否かを判断する場合に比べ、ガスセンサAの点244で示される振幅AM1 ́と、ガスセンサBの点246で示される振幅AM1 ́との差が大きい。このため、振幅AM1 ́に基づき内燃機関100が正常であるか否かを判断する場合、振幅AM1に基づき内燃機関100が正常であるか否かを判断する場合に比べ、対応閾値TH2の設定が容易である。すなわち、第1の実施形態の初期異常検知処理によれば、ガスセンサ1の出力値のバラツキを考慮して、インバランス割合に初期異常があるか否かの観点から、内燃機関100の初期異常を検知することができる。次に、図9及び図10を参照して、第2の実施形態の初期異常検知処理について説明する。第1の実施形態では、ガスセンサ1個体の変化の割合と、基準割合とから、ガスセンサ1の出力値の補正量を算出していた。第2の実施形態では、ガスセンサ1個体の変化の割合と、基準割合とに基づき、対応閾値TH2を補正する。第2の実施形態の初期異常検知処理は、例えば、内燃機関100と、センサ制御装置4とを備える自動車の製造検査または出荷検査時に、開始の指示が入力された場合に起動される。図9のフローチャートに示す各処理を実行させるプログラムは、図2のROM7に記憶されており、CPU6によって実行される。以下、図10に示す第1の実施形態と同様の具体例を用いて、第2の実施形態の初期異常検知処理を説明する。図9において、図4に示す第1の実施形態の初期異常検知処理と同様の処理を行うステップには、同じステップ番号を付与している。図9に示すように、第2の実施形態の初期異常検知処理は、ステップS80に代えてステップS85が実行され、ステップS90及びステップS100に代えてステップS95が実行され、ステップS110に代えてステップS115が実行される点で、第1の実施形態の初期異常検知処理とは異なる。第1の実施形態と同様な処理については説明を省略し、以下、第1の実施形態と異なる、ステップS85、ステップS95及びステップS115の処理について説明する。ステップS85では、ROM7に記憶されている基準割合と、対応閾値TH2とが取得される。ステップS95では、ステップS70で算出された変化の割合と、ステップS85で取得された基準割合とに基づき、対応閾値TH2が補正され、補正後の対応閾値TH2は、補正閾値TH2 ́としてRAM8に記憶されるS95。ステップS95では、例えば、式4に従って、補正閾値TH2 ́が算出される。補正閾値TH2 ́=対応閾値TH2ステップS70算出された変化の割合/ステップS85で取得された基準割合式4線分208で表す対応閾値TH2及び補正閾値TH2 ́は、それぞれ基準閾値TH1に対応している。ステップS95の処理によって、ガスセンサAに対する補正閾値TH2 ́に、例えば、図10の線分308で示される値が設定されたとする。ガスセンサBに対する補正閾値TH2 ́に、例えば、線分318で示される値が設定されたとする。本実施形態では、線分308は、線分209と、直線201との交点を通る。線分318は、線分209と、直線203との交点を通る。ステップS115では、振幅AM1が、補正閾値TH2 ́よりも大きいか否かが判断される。点204で示されるガスセンサAの出力値の振幅AM1は、線分308で表されるガスセンサAに対する補正閾値TH2 ́以下であるS115:NO。この場合、内燃機関Aはインバランス割合の観点から正常であるため、故障警告灯105が点灯されずに、初期異常検知処理は終了する。点206で示されるガスセンサBの出力値の振幅AM1は、線分318で表されるガスセンサBに対する補正閾値TH2 ́よりも大きいS115:YES。したがって、内燃機関Bに初期異常があることを報知するために、故障警告灯105が点灯された後S120、第2の実施形態の初期異常検知処理は終了する。上記第2の実施形態の初期異常検知処理において、ステップS85は、本発明の基準割合取得工程に相当する。ステップS95は、本発明の閾値補正に相当する。ステップS115と、ステップS120とは、本発明の第2異常検知工程に相当する。ステップS85と、ステップS95と、ステップS115と、ステップS120とは、本発明の第1異常検知工程に相当する。第2の実施形態の初期異常検知処理によれば、下記の効果が得られる。第2の実施形態の初期異常検知処理では、内燃機関100が正常であるか否かを判断するための補正閾値TH2 ́が、対応閾値TH2に基づきガスセンサ1毎に設定される。したがって、振幅AM1と補正閾値TH2 ́とを比較して内燃機関100が正常であるか否かを判断する場合、振幅AM1と対応閾値TH2とを比較して内燃機関100が正常であるか否かを判断する場合に比べ、対応閾値TH2基準閾値TH1の設定が容易である。すなわち、第1の実施形態の初期異常検知処理によれば、ガスセンサ1の出力値のバラツキを考慮して、インバランス割合に初期異常があるか否かの観点から、内燃機関100の初期異常を検知することができる。本発明は、以上詳述した実施形態に限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変更が適宜加えられてもよい。例えば、以下の1、2の変形が適宜加えられてもよい。1本発明は空燃比を検知可能な種々のガスセンサに適用可能である。初期異常の検知対象となる多気筒内燃機関の用途、気筒数及び形式は、適宜変更されてよい。上記実施形態では、ECU3において、初期異常検知処理が実行されていたが、初期異常検知処理の一部又は全部は、ECU3とは別途設けられた制御装置例えば、製造検査用のコンピュータにおいて実行されてもよい。ECU3は、ASIC30を備えていたが、ECU3と、ASIC30は別々に設けられてもよい。2初期異常検知処理の各処理は、必要に応じて適宜変更されてもよい。例えば、特徴的な部分が組み合わされて実行されてもよい。具体的には、以下の2-1及び2-2の変形が加えられてもよい。2-1ガスセンサ1の出力値の振幅を算出する処理は適宜変更されてよい。ガスセンサ1の出力値のノイズを除去するための処理は、公知の各種手法が適用されてもよいし、適宜省略されてもよい。他の例では、振幅の算出に用いるサンプリング期間及びサンプリング周期は適宜変更されてよい。ガスセンサ1個体の変化の割合を算出するための内燃機関100の運転条件は、インバランス割合が異なる複数の条件であればよい。例えば、インバランス割合が異なる条件の数は、上記実施形態の2つに限定されない。同様に、内燃機関100の初期異常を検知する処理に用いる、ガスセンサ1の出力値の振幅の数は、1つ以上であればよい。2-2第1の実施形態の図4のステップS110では、振幅AM1を補正して求めた振幅AM1 ́と、対応閾値TH2とが比較されていたが、これに限定されない。よより具体的には、例えば、内燃機関100に初期異常があるか否かは、ガスセンサ1の出力値の振幅と、ガスセンサ1個体の変化の割合と、インバランス割合の異常があるか否かとの対応を記憶したマップを用いて判断されてもよい。この場合、例えば、次のように内燃機関100に初期異常があるか否かが判断されればよい。まず、図4のステップS5からステップS70までの処理が実行され、ガスセンサ1個体の変化の割合が算出される。ステップS80、ステップS90及びステップS100は省略される。ステップS110では、ガスセンサ1個体の変化の割合と、ステップS30で算出された出力値の振幅と、マップに記憶された対応関係とに基づき、インバランス割合が異常であるか否かが判断されればよい。他の例では、例えば、内燃機関100に初期異常があるか否かは、ガスセンサ1の出力値の振幅と、ガスセンサ1個体の変化の割合とから推定されるインバランス割合と基準閾値とを比較して判断されてもよい。この場合、例えば、次のように内燃機関100に初期異常があるか否かが判断されればよい。図4のステップS5からステップS70までの処理が実行され、ガスセンサ1個体の変化の割合が算出される。ステップS80、ステップS90及びステップS100の処理に代えて、ガスセンサ1個体の変化の割合と、ステップS30で算出された出力値の振幅とから、第1条件で内燃機関100が駆動された場合のインバランス割合が推定される。ステップS110では、推定されたインバランス割合と、基準閾値とが比較される。推定されるインバランス割合が基準閾値よりも大きい場合にインバランス割合が異常であると判断されればよい。他の例では、補正閾値は、基準閾値及対応閾値の少なくともいずれかを、変化の割合に基づき補正した値であればよい。したがって、例えば、内燃機関100に初期異常があるか否かは、ガスセンサ1の出力値の振幅と、基準割合とから推定されるインバランス割合と、ガスセンサ1個体の変化の割合に基づき補正された基準閾値補正閾値とを比較して判断されてもよい。この場合、例えば、次のように内燃機関100に初期異常があるか否かが判断されればよい。図9に示す初期異常検知処理において、ステップS5からステップS70までの処理が実行され、ステップS85では基準割合及び基準閾値TH1が取得される。ステップS95では、ステップS70で取得されたガスセンサ1個体の変化の割合と、ステップS85で取得された基準割合とに基づき、基準閾値TH1が補正され、補正された基準閾値TH1は、補正閾値TH1 ́としてRAM8に記憶される。図10に例示する具体例では、例えば、線分308と、直線202との交点を通り縦軸に平行な線分309で表される値がガスセンサAの補正閾値TH1 ́として取得される。ステップS115では、ステップS30で算出された出力値の振幅と、基準割合とから推定されるインバランス割合と、補正閾値TH1 ́とを比較される。推定されるインバランス割合が補正閾値TH1 ́よりも大きい場合に、インバランス割合の観点から内燃機関100に初期異常があると判断されればよい。図10に示すガスセンサAについての具体例では、点304で表されるインバランス割合と、線分309で表される補正閾値TH1 ́とが比較され、インバランス割合の観点から内燃機関100に初期異常はないと判断される。"