(57)【要約】      （修正有）【課題】液体表面に発生する泡が要因となって起こる液面検知の誤検知を、より高い精度で検出できるようにする。【解決手段】自動分析装置１０において、特徴量抽出部（演算部７１）が、分注プローブ１ａが下降を開始してから一定時間が経過するまでに、発振回路４が出力した交流信号の発振周波数の時系列データから特徴量を抽出し、抽出した特徴量をデータ系列として出力する。そして、泡接触判定処理部７２Ａが、特徴量のデータ系列の波形と、異常波形モデルとの相互相関を計算し、相互相関の計算結果に基づいて液面の検知が正常に行われたかどうかを判定する。さらに、該判定結果に基づいて、第２制御部８が、分注プローブ１ａの先端部１ｅと容器２内の液面との乖離、及び乖離の要因を判定する。【選択図】図１
2020112380.tif 000002
【特許請求の範囲】
【請求項１】
分注プローブを有し、当該分注プローブの先端部を容器内の液面に移動するとともに、液体を吸引及び吐出するように構成された分注部と、
前記分注プローブと接続された、前記分注プローブの先端部と周辺部との間の静電容量に応じた発振周波数の交流信号を出力する発振回路と、
前記発振回路から出力される前記交流信号の発振周波数に基づいて、前記分注プローブの先端部が前記容器内の液面に接触したかどうかを検知する検出部と、
前記検出部の検知結果に基づき前記分注部の動作を制御する第１の制御部と、
前記分注プローブが下降を開始してから一定時間が経過するまでの間に前記発振回路が出力した前記交流信号の発振周波数の時系列データを一定区間ごとに分け、前記一定区間ごとの特徴量を抽出し、前記抽出した特徴量をデータ系列として出力する特徴量抽出部と、
前記特徴量抽出部から出力された前記特徴量のデータ系列の波形と、前記分注プローブの先端部が前記容器内の液面の泡に接触した場合に観察される波形に基づく異常波形モデルとの相互相関を計算し、前記相互相関の計算結果に基づいて前記液面の検知が正常に行われたかどうかを判定する泡接触判定処理部と、
前記泡接触判定処理部の判定結果に基づいて、前記分注プローブの先端部と前記容器内の液面との乖離、及び前記乖離の要因を判定する第２の制御部と、を備える
自動分析装置。
【請求項２】
前記特徴量は、前記発振周波数の時系列データの一定区間ごとの最小値及び／又は最大値である
請求項１に記載の自動分析装置。
【請求項３】
前記異常波形モデルは、数ミリ秒に対応する区間においてピーク値が連続する台形形状の波形である
請求項２に記載の自動分析装置。
【請求項４】
前記異常波形モデルとの相関関係が計算される前記特徴量のデータ系列の波形は、前記特徴量のデータ系列から、前記分注プローブの先端部が前記容器内の液面に接触したと前記検出部が判定する前の段階で前記特徴量抽出部によって抽出されていた特徴量が削除されたデータ系列の波形である
請求項３に記載の自動分析装置。
【請求項５】
前記発振周波数の時系列データの一定区間ごとの特徴量の微分値を算出し、前記微分値のデータ系列の波形と対応する正常な波形との相互相関を計算し、計算結果から前記液面の検知が正常に行われたかどうかを判定する第１の判定処理部と、
前記発振周波数の時系列データの一定区間ごとの特徴量の微分値を算出し、前記微分値のデータ系列の最大値と閾値とを比較し、比較結果から前記液面の検知が正常に行われたかどうかを判定する第２の判定処理部と、
前記発振周波数の時系列データから抽出された一定区間ごとの特徴量が所定の条件を満たす区間を算出し、前記区間の長さと閾値とを比較し、比較結果から前記液面の検知が正常に行われたかどうかを判定する第３の判定処理部と、
前記発振周波数の時系列データから抽出された一定区間ごとの特徴量としての最大値のデータ系列及び最小値のデータ系列から各データ系列の最大値を算出し、各データ系列の最大値間の差分と閾値とを比較し、比較結果から前記液面の検知が正常に行われたかどうかを判定する第４の判定処理部と、をさらに備え、
前記第２の制御部は、前記泡接触判定処理部、並びに、前記第１の判定処理部から前記第４の判定処理部までの各判定処理部による判定結果のうちの２以上の判定結果の組み合わせにより、前記分注プローブの先端部と前記容器内の液面との乖離の判定、及び、前記乖離の要因を判定する
請求項１〜４のいずれか一項に記載の自動分析装置。
【請求項６】
前記第２の制御部は、前記乖離の要因として、液面に発生する泡、前記分注プローブの先端部と前記容器の内壁との接触又は静電気、若しくは不明を判定する
請求項１〜５のいずれか一項に記載の自動分析装置。
【請求項７】
前記第２の制御部は、前記分注プローブの先端部と前記容器内の液面との乖離の判定結果、及び、前記乖離の要因についての判定結果を表示部に出力するよう制御する
請求項１〜６のいずれかに記載の自動分析装置。
【請求項８】
容器内の液面に移動され当該容器内の液体を吸引及び吐出する分注プローブと接続された発振回路が、前記分注プローブの先端部と周辺部との間の静電容量に応じた発振周波数の交流信号を出力する手順と、
前記分注プローブが下降を開始してから一定時間が経過するまでの間に前記発振回路が出力した前記交流信号の発振周波数の時系列データを一定区間ごとに分け、前記一定区間ごとの特徴量を抽出し、前記抽出した特徴量をデータ系列として出力する手順と、
前記出力された前記特徴量のデータ系列の波形と、前記分注プローブの先端部が前記容器内の液面の泡に接触した場合に観察される波形に基づく異常波形モデルとの相互相関を計算し、前記相互相関の計算結果に基づいて前記液面の検知が正常に行われたかどうかを判定する手順と、
前記判定の結果に基づいて、前記分注プローブの先端部と前記容器内の液面との乖離の判定、及び、前記乖離の要因を判定する手順と、
をコンピューターに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、自動分析装置及びプログラムに関し、特に分注プローブの先端部が容器内の液体の表面（以下「液面」と記す）に対する誤検知を検知する技術に関する。
【背景技術】
【０００２】
臨床検査用の自動分析装置は、生化学検査、免疫検査、輸血検査などさまざまな分野での検査に用いられる。このような自動分析装置では、予め装置上の決められた位置に試料や試薬を格納する容器が配置されている。自動分析装置は、その容器に入っている液体（試料、試薬）を分注プローブとそれに接続されたポンプからなる分注機構を用いて吸引、保持しながらの運搬、そして目的の容器（反応容器など）への吐出を行っている。
【０００３】
この分注機構においては、分注プローブや液体の汚染を避けるため、液面で分注プローブが停止する動作が望ましい。液面位置を検知する方法として現在広く用いられているのが、静電容量方式である。
【０００４】
静電容量方式は、分注プローブの吸引部（先端部）と周辺部（例えば装置筐体のグラウンド）との静電容量をモニタリングし、その変化を検知する方式である。分注プローブの先端部が試料等の液体に触れると静電容量値が変化するため、閾値等を用いてその変化を検知すれば、分注プローブが液面に接触しているか空中にあるかを識別することが可能である。通常は、分注プローブが容器へ降下する際に静電容量のモニタリングを行い、静電容量が閾値を超えて大きく変化した場合、その時点で分注プローブの下降を停止させる。それにより、分注プローブの先端部のみが液面に触れた状態で、分注プローブを保持することができる。
【０００５】
静電容量方式の具体例としては、例えば特許文献１に示されたＣＲ発振回路を構成し、発振周波数の変化から静電容量Ｃの変化を求める方法（特許文献１を参照）の他、類似又は関連する多数の方法（例えば特許文献２〜４を参照）が考案されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開平１１−２７１３２３号公報
【特許文献２】特開平１０−２０６４３２号公報
【特許文献３】特開２０１１−０９４９８５号公報
【特許文献４】特開２００４−０２８６７３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
従来の静電容量方式の液面検知に関する技術的な問題として、試料等の液面と分注プローブ先端が乖離しているにもかかわらず、分注プローブの先端が液面に到達したと判定してしまう事象（以下「液面乖離」と呼ぶ）がある。このとき分注プローブは空中にあるため、目的の試料を吸引できない。したがって対象試料の目的成分の濃度としてゼロ付近の値を、誤って報告してしまうリスクがある。誤った液面の検知の主要因の一つに、液体表面に発生する泡がある。
【０００８】
泡によって起こる問題は、分注プローブが液体の泡の表面に触れた時点で分注プローブが停止してしまうということである。泡の内部は空気であり、分注プローブは空気を吸引してしまう。特許文献１及び２のいずれにも、誤検知の要因が泡であるとを判定することについては記載されていない。特許文献３，４に記載の技術は、液面を検知した高さに基づいているので、泡の検知にはある程度有効と考えられる。一方で、泡が分注プローブと接触したときの泡の状態の変化、ひいては静電容量の変化は多様なバリエーションがあると考えられるため、高さという情報で判断するのは精度の点で疑問がある。
【０００９】
さらに、特許文献３，４に記載された液面の高さによって判断する方法では、泡以外の要因で誤検知が起きる可能性を否定できず、誤検知の要因が泡であると通知することは難しい。ある試料の液面乖離の可能性をユーザーに通知したとき、要因が複数あるため、ユーザーがどのような対処をすべきか判断するのは困難である。異常の可能性がある試料を再検査する際、その判断をするために時間を要することは、検査の遅延を招く。これは最終的に、検査結果を待つ患者の利益を損なうものである。
【００１０】
本発明は、上記の状況を考慮してなされたものであり、液体表面に発生する泡が要因となって起こる液面検知の誤検知を、より高い精度で検出できるようにするものである。
【課題を解決するための手段】
【００１１】
本発明の一態様の自動分析装置は、分注プローブを有し、当該分注プローブの先端部を容器内の液面に移動するとともに、液体を吸引及び吐出するように構成された分注部と、分注プローブと接続された、分注プローブの先端部と周辺部との間の静電容量に応じた発振周波数の交流信号を出力する発振回路と、発振回路から出力される交流信号の発振周波数に基づいて、分注プローブの先端部が容器内の液面に接触したかどうかを検知する検出部と、検出部の検知結果に基づき分注部の動作を制御する第１の制御部と、分注プローブが下降を開始してから一定時間が経過するまでの間に発振回路が出力した交流信号の発振周波数の時系列データを一定区間ごとに分け、一定区間ごとの特徴量を抽出し、抽出した特徴量をデータ系列として出力する特徴量抽出部と、特徴量抽出部から出力された特徴量のデータ系列の波形と、分注プローブの先端部が容器内の液面の泡に接触した場合に観察される波形に基づく異常波形モデルとの相互相関を計算し、相互相関の計算結果に基づいて液面の検知が正常に行われたかどうかを判定する泡接触判定処理部と、泡接触判定処理部の判定結果に基づいて、分注プローブの先端部と容器内の液面との乖離、及び乖離の要因を判定する第２の制御部と、を備える。
【発明の効果】
【００１２】
本発明の少なくとも一態様によれば、泡接触判定処理部の判定結果に基づいて、液体表面に発生する泡が要因となって起こる液面検知の誤検知を、より高い精度で検出できるようになる。
上記した以外の課題、構成及び効果は、以下の実施形態の説明により明らかにされる。
【図面の簡単な説明】
【００１３】
【図１】本発明の一実施形態に係る自動分析装置の構成例を示すブロック図である。
【図２】本発明の一実施形態に係る自動分析装置が備えるコンピューターのハードウェア構成例を示すブロック図である。
【図３】本発明の一実施形態に係る自動分析装置の第１処理部による液面検知時の動作の手順の例を示すフローチャートである。
【図４】本発明の一実施形態に係る液面を正常に検知したときの静電容量波形の例を示すグラフである。
【図５】本発明の一実施形態に係る自動分析装置の第２処理部７による液面乖離判定及び乖離要因判定の手順の例を示すフローチャートである。
【図６】本発明の一実施形態に係る泡接触判定処理、第１の判定処理〜第４の判定処理の各判定結果及び乖離要因の組み合わせを示した液面乖離及び乖離要因判定テーブルの一例を示す図である。
【図７】本発明の一実施形態に係る乖離要因判定結果に対する推奨対処手順が登録された推奨対処手順テーブルの一例を示す図である。
【図８】本発明の一実施形態に係る泡接触判定処理の手順の例を示すフローチャートである。
【図９】本発明の一実施形態に係る静電容量波形を処理したときのイメージを表す図である。
【図１０】本発明の一実施形態に係る正常波形に対するデータ処理及び判定処理を示す図である。
【図１１】本発明の一実施形態に係る分注プローブの先端部が泡に接触して停止したときに起こり得る静電容量波形（山形）の例、データ処理及び判定処理を示す図である。
【図１２】本発明の一実施形態に係る分注プローブの先端部が泡に接触して停止したときに起こり得る静電容量波形（谷形）の例、データ処理及び判定処理を示す図である。
【図１３】本発明の一実施形態に係る実際に泡に対する液面検知を行った場合に観測された波形の実例を示す図である。
【図１４】本発明の一実施形態に係る第１の判定処理の手順の例を示すフローチャートである。
【図１５】本発明の一実施形態に係る静電容量波形を処理したときのイメージを表す図である。
【図１６】本発明の一実施形態に係る第２の判定処理の手順の例を示すフローチャートである。
【図１７】本発明の一実施形態に係る静電容量波形を処理したときのイメージを表す図である。
【図１８】本発明の一実施形態に係る第３の判定処理の手順の例を示すフローチャートである。
【図１９】本発明の一実施形態に係る静電容量波形を処理したときのイメージを表す図である。
【図２０】本発明の一実施形態に係る第４の判定処理の手順の例を示すフローチャートである。
【図２１】本発明の一実施形態に係る静電容量波形を処理したときのイメージを表す図である。
【発明を実施するための形態】
【００１４】
以下、本発明を実施するための形態の例について、添付図面を参照しながら説明する。なお、各図において実質的に同一の機能又は構成を有する構成要素については、同一の符号を付して重複する説明を省略する。
【００１５】
＜自動分析装置の全体構成＞
図１は、一実施形態に係る自動分析装置の構成例を示すブロック図である。
本実施形態に係る自動分析装置１０は、試料や試薬などの液体を分注する分注プローブと周辺部との静電容量をモニタリングすることにより、分注プローブ先端が液面に接触したときの静電容量値の変化を検出し、検出した変化に基づいて液面を検知する機構を備える。
【００１６】
図１に示すように、自動分析装置１０は、分注ユニット１、ＣＲ発振回路４、第１処理部５、記憶部６、第２処理部７、第２制御部８及び表示部９を備える。
【００１７】
分注ユニット１（分注部の一例）は、分注プローブ１ａ、分注プローブ１ａを保持するアーム１ｂ、及び、アーム１ｂを駆動する駆動機構１ｃ等から構成される。分注プローブ１ａは、金属等の導電性部材で構成されており、吸引した液体が収まるよう空洞になっている。分注プローブ１ａの外周面は、液体を吸引及び吐出する先端部１ｅ（吸引部）を除いてシールドで保護されている。
【００１８】
分注ユニット１は、駆動機構１ｃ及びアーム１ｂによって、分注プローブ１ａの水平移動、上下移動、吸引及び吐出の動作が可能に構成される。駆動機構１ｃは、例えばアーム１ｂを駆動する不図示のモーターを有する。分注ユニット１は、駆動機構１ｃ及びアーム１ｂにより分注プローブ１ａの先端部１ｅを容器２の試料３（液体の一例）の液面に移動し、試料３を吸引する。そして、分注ユニット１は、試料３を吸引した分注プローブ１ａを目的の場所（反応容器など）へ移送し、試料３を吐出する。これらの分注ユニット１の動作は、第１処理部５（第１制御部５２）によって制御される。駆動機構１ｃは、ＣＲ発振回路４及び第１処理部５に接続される。
【００１９】
容器２は、ターンテーブル等の保持部材によって保持された状態で移送される。保持部材は、自動分析装置１０の筐体（以下「装置筐体」）のグラウンドに電気的に接地されている。
【００２０】
ＣＲ発振回路４（発振回路の一例）は、抵抗（Ｒ）及びコンデンサ（Ｃ）で構成されるＣＲ回路を用いて帰還する帰還型の発振回路であり、正弦波の交流信号を発生する。ＣＲ発振回路４は、分注プローブ１ａの先端部１ｅと周辺部（例えば保持部材、すなわち装置筐体のグラウンド）との静電容量に応じた発振周波数の交流信号を出力する。ＣＲ発振回路４は、不図示のアナログ−デジタル変換回路を有し、デジタルの交流信号を出力する。なお、発振回路には、ＣＲ発振回路の他にも、ＬＣ発振回路などその他種々の発振回路を適用可能である。
【００２１】
第１処理部５は、検出部５１及び第１制御部５２を有し、ＣＲ発振回路４から出力される交流信号の解析、及び、解析結果に基づく分注ユニット１の制御を行う。
【００２２】
検出部５１は、ＣＲ発振回路４と接続されており、ＣＲ発振回路４から出力される交流信号の発振周波数を監視する。すなわち、検出部５１は、交流信号を取り込み、所定のサンプリング周期で交流信号の発振周波数を検出する。そして、検出部５１は、検出した発振周波数に基づいて、分注プローブ１ａの先端部１ｅが容器２内の液面に接触したかどうかを検知し、液面を検知した場合には液面検知を示す信号を出力する。
【００２３】
第１制御部５２（第１の制御部）は、検出部５１の検出結果に基づいて分注ユニット１の駆動機構１ｃに制御信号を出力し、分注プローブ１ａの動作を制御する。
【００２４】
記憶部６は、分注プローブ１ａが下降を開始してから一定時間が経過するまでにＣＲ発振回路４が出力したデジタルの交流信号のデータを記憶する。
【００２５】
第２処理部７は、第１処理部５による液面検知に誤りがないかどうかを検知する処理を行う。第２処理部７は、演算部７１、泡接触判定処理部７２Ａ、第１液面乖離判定処理部７２Ｂ、第２液面乖離判定処理部７２Ｃ、第３液面乖離判定処理部７２Ｄ及び第４液面乖離判定処理部７２Ｅを有する。
【００２６】
演算部７１（特徴量抽出部の一例）は、記憶部６に記憶されている交流信号のデータから交流信号の発振周波数の時系列データ（静電容量波形に相当）を取得する。そして、その発振周波数の時系列データから特徴量を抽出してデータ系列として出力する。特徴量の抽出は、例えば、発振周波数の時系列データの一定区間ごとに行われる。
【００２７】
泡接触判定処理部７２Ａ、第１液面乖離判定処理部７２Ｂ〜第４液面乖離判定処理部７２Ｅは、演算部７１から出力された特徴量のデータ系列を用いて、容器２内の液面の検知が正常に行われたかどうかを判定する。下記に泡接触判定処理部７２Ａ、第１液面乖離判定処理部７２Ｂ〜第４液面乖離判定処理部７２Ｅの動作の概要を説明するが、各判定処理の詳細については、図８，図１４，図１６，図１８，図２０を参照して後述する。
【００２８】
泡接触判定処理部７２Ａは、演算部７１から出力される、発振周波数の時系列データの一定区間ごとの特徴量（例えば最小値）のデータ系列の波形と、異常波形モデルとの相互相関を計算する。そして、泡接触判定処理部７２Ａは、相互相関の計算結果から液面の検知が正常に行われたかどうかを判定し、判定結果を第２制御部８へ出力する。異常波形モデルは、分注プローブ１ａの先端部１ｅが容器２内の液面の泡に接触した場合に特徴的に観察される波形に基づく、波形モデルである。演算部７１及び泡接触判定処理部７２Ａによる泡接触判定処理では、発振周波数の時系列データによる波形（静電容量波形）の変化の形を見る。
【００２９】
第１液面乖離判定処理部７２Ｂ（第１の判定処理部の一例）は、演算部７１から出力される、発振周波数の時系列データの一定区間ごとの特徴量（例えば最小値）の微分値のデータ系列の波形と、対応する正常な波形との相互相関を計算する。そして、第１液面乖離判定処理部７２Ｂは、計算結果から液面の検知が正常に行われたかどうかを判定し、判定結果を第２制御部８へ出力する。演算部７１及び第１液面乖離判定処理部７２Ｂによる第１の判定処理おいても、発振周波数の時系列データによる波形（静電容量波形）の変化の形を見る。
【００３０】
第２液面乖離判定処理部７２Ｃ（第２の判定処理部の一例）は、演算部７１から出力される、発振周波数の時系列データの一定区間ごとの特徴量（例えば最小値）の微分値のデータ系列の最大値と、閾値とを比較する。そして、第２液面乖離判定処理部７２Ｃは、比較結果から液面の検知が正常に行われたかどうかを判定し、判定結果を第２制御部８へ出力する。演算部７１及び第２液面乖離判定処理部７２Ｃによる第２の判定処理では、発振周波数の時系列データによる波形（静電容量波形）の変化の大きさ（急峻さ）を見る。
【００３１】
第３液面乖離判定処理部７２Ｄ（第３の判定処理部の一例）は、演算部７１から出力される、発振周波数の時系列データから抽出された一定区間ごとの特徴量（例えば最小値）が、所定の条件を満たす区間を算出する。そして、第３液面乖離判定処理部７２Ｄは、その区間の長さと閾値とを比較し、比較結果から液面の検知が正常に行われたかどうかを判定し、判定結果を第２制御部８へ出力する。演算部７１及び第３液面乖離判定処理部７２Ｄによる第３の判定処理では、発振周波数の時系列データによる波形（静電容量波形）の形状を見る。
【００３２】
第４液面乖離判定処理部７２Ｅ（第４の判定処理部の一例）は、演算部７１から出力される、発振周波数の時系列データから抽出された一定区間ごとの特徴量としての最大値のデータ系列及び最小値のデータ系列から、各データ系列の最大値を算出する。そして、第４液面乖離判定処理部７２Ｅは、各データ系列の最大値間の差分と閾値とを比較し、比較結果から液面の検知が正常に行われたかどうかを判定し、判定結果を第２制御部８へ出力する。演算部７１及び第４液面乖離判定処理部７２Ｅによる第４の判定処理では、発振周波数の時系列データによる波形（静電容量波形）の安定性（ノイズの有無など）を見る。
【００３３】
第２制御部８は、泡接触判定処理部７２Ａ、第１液面乖離判定処理部７２Ｂ〜第４液面乖離判定処理部７２Ｅの各判定結果の組み合わせにより、分注プローブ１ａの先端部１ｅと容器２内の液面との乖離、及び、乖離の要因を判定する。第２制御部８の判定結果は、表示部９に出力されて表示部９の画面に表示される。なお、泡との接触が原因となって起こる液面乖離判定については、泡接触判定処理部７２Ａによる判定結果のみに基づいて行ってもよい。
【００３４】
＜コンピューターのハードウェア構成＞
図２は、自動分析装置１０が備えるコンピューターのハードウェア構成例を示すブロック図である。
【００３５】
コンピューター２０は、バス２４にそれぞれ接続されたＣＰＵ（Central Processing Unit）２１、ＲＯＭ（Read Only Memory）２２及びＲＡＭ（Random Access Memory）２３を備える。さらに、コンピューター２０は、表示部２５、操作部２６、不揮発性ストレージ２７及びネットワークインターフェース２８を備える。
【００３６】
ＣＰＵ２１は、本実施形態に係る各機能を実現するソフトウェアのプログラムコードをＲＯＭ２２から読み出して実行する。例えば、ＣＰＵ２１により、第１処理部５、第２処理部７及び第２制御部８の各機能を実現することができる。なお、コンピューター２０は、ＣＰＵ２１の代わりに、ＭＰＵ（Micro-Processing Unit）等の処理装置を備えてもよい。
【００３７】
ＲＡＭ２３には、演算処理の途中に発生した変数やパラメータ等が一時的に書き込まれる。例えばＲＡＭ２３が、図１の記憶部６として、ＣＲ発振回路４から出力されたデジタルの交流信号のデータを記憶してもよい。
【００３８】
表示部２５は図１の表示部９に相当し、コンピューター２０で行われる処理の結果等を表示する。例えば表示部２５は、液晶ディスプレイモニタである。操作部２６には、例えば、キーボード、マウス又はタッチパネル等が用いられ、ユーザーが所定の操作入力、指示を行うことが可能である。
【００３９】
不揮発性ストレージ２７としては、例えば、ＨＤＤ（Hard Disk Drive）、ＳＳＤ（Solid State Drive）、フレキシブルディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、磁気テープ、不揮発性のメモリカード等が用いられる。この不揮発性ストレージ２７には、ＯＳ（Operating System）、各種のパラメータの他に、コンピューター２０を機能させるためのプログラムが記録される。例えば、不揮発性ストレージ２７には、液面乖離判定及び乖離要因判定処理（図５参照）を規定したプログラム、液面乖離及び乖離要因判定テーブル４１（図６参照）、及び、推奨対処手順テーブル４２（図７参照）等が記憶される。
【００４０】
ネットワークインターフェース２８には、例えば、ＮＩＣ（Network Interface Card）等が用いられ、ＬＡＮ等のネットワークＮを介して各装置間で各種のデータを送受信することが可能である。
【００４１】
＜液面検知時の動作＞
図３は、自動分析装置１０の第１処理部５による液面検知時の動作の手順の例を示すフローチャートである。この処理は、図１の第１処理部５、分注ユニット１及びＣＲ発振回路４により実現される。
【００４２】
前提として、第１処理部５の検出部５１は、ＣＲ発振回路４が出力する交流信号の発振周波数を逐次モニタリングしている。まず、第１処理部５の第１制御部５２は、検査対象の容器２が搬送されてくると、設定されたタイミングで分注プローブ１ａの下降動作を開始する（ステップＳ１）。
【００４３】
このとき、検出部５１による発振周波数のモニタリングは継続される（ステップＳ２）。このモニタリングと並行して、検出部５１は、記憶部６に予め設定された、発振周波数に対する閾値を記憶部６から読み出す（ステップＳ３）。そして、検出部５１は、モニタリングした発振周波数が閾値以上であるかどうかを判定する（ステップＳ４）。
【００４４】
もし分注プローブ１ａの先端部１ｅが容器２内の液面に接触した場合には、分注プローブ１ａと周辺部との間の静電容量値は顕著に増加するため、発振周波数は閾値を超える。ステップＳ４で、発振周波数が閾値以上であると判定した場合（ステップＳ４がＹＥＳ判定の場合）、検出部５１は第１制御部５２に対して、容器２内の液面を検知したことを示す信号（液面検知オン信号）を出力する。第１制御部５２は、液面検知オン信号を受信した後、速やかに分注ユニット１の駆動機構１ｃに分注プローブ１ａの下降停止の制御信号を送信する（ステップＳ５）。これにより、駆動機構１ｃは分注プローブ１ａの下降動作を停止する。そして、分注プローブ１ａの先端部１ｅは容器２内の液面に接触した状態となる。
【００４５】
一方、ステップＳ４で、発振周波数は閾値未満であると判定された場合（ステップＳ４がＮＯ判定の場合）、第１制御部５２は、ステップＳ１に移行して分注プローブ１ａの下降動作を継続する。
【００４６】
上述した第１処理部５を介して行われる液面検知の動作の手順は、従来技術を利用したものであり、本発明の一実施形態の前提となるものである。なお、この液面検知の動作の手順は、ハードウェア又はソフトウェアのいずれで実現してもよい。
【００４７】
＜正常な静電容量の波形例＞
図４は、液面を正常に検知したときの静電容量波形の例を示すグラフである。図４では、横軸に時間、縦軸に発振周波数（静電容量値に相関）を示している。検査対象の試料は正常な状態を想定している。
【００４８】
液面検知の動作開始から終了（分注プローブ１ａが停止）の後も、ＣＲ発振回路４は交流信号を出力している。記憶部６は、分注プローブ１ａの下降開始から分注プローブ１ａが液面に触れて停止した状態を経由し一定時間経過するまでの交流信号のデータを、時系列データとして記憶する。このときの静電容量の変化は、図４に示すように交流信号の発振周波数の変化（静電容量波形）として現れる。
【００４９】
分注プローブ１ａの先端部１ｅが周辺部から一定の距離離れている場合には、波形３１の発振周波数は、所定の基準値付近の値もしくは基準値よりも大きい値となる。分注プローブ１ａが下降する間、静電容量はほとんど変化しないか若しくはわずかに変化するだけであり、破線で示す閾値を超えることはない。しかし、分注プローブ１ａの先端部１ｅが液面に接触することで（タイミング３１ａ）、静電容量値は閾値を超え、分注プローブ１ａは停止する（タイミング３１ｂ）。分注プローブ１ａの停止直後、発振周波数は変動しているが、これは液面が揺らいでいるためであると考えられる。その後、分注プローブ１ａの先端部１ｅ及び周辺部の状態（液面接触）に変化はないため、静電容量値はほぼ一定値に収束する。
【００５０】
＜液面乖離判定及び液面乖離要因判定＞
図５は、自動分析装置１０の第２処理部７による液面乖離判定及び乖離要因判定の手順の例を示すフローチャートである。この処理は、図１の第２処理部７及び第２制御部８により実現される。
【００５１】
第２処理部７では、まず、泡接触判定処理部７２Ａ、第１液面乖離判定処理部７２Ｂ〜第４液面乖離判定処理部７２Ｅ（図１参照）のそれぞれが、泡接触判定、第１の判定〜第４の判定の各演算処理を実行する（ステップＳ１１）。ＣＲ発振回路４から出力されて記憶部６（図１）に格納された時系列データは、まず演算部７１によって、短い一定区間のまとまりに分けられ、それらの一定区間ごとに特徴量（最小値、最大値等）が算出及び抽出される。そして、抽出された特徴量は、データ系列として演算部７１から出力される。演算部７１から出力されたデータ系列は、泡接触判定処理部７２Ａ、第１液面乖離判定処理部７２Ｂ〜第４液面乖離判定処理部７２Ｅのそれぞれにおいて処理される。そして、泡接触判定処理部７２Ａ、第１液面乖離判定処理部７２Ｂ〜第４液面乖離判定処理部７２Ｅのそれぞれは、正常（Ｐａｓｓ）又は異常（Ｆａｉｌ）の二値で表される判定結果値を、第２制御部８へ出力する（ステップＳ１２）。
【００５２】
第２制御部８は、泡接触判定処理部７２Ａ、第１液面乖離判定処理部７２Ｂ〜第４液面乖離判定処理部７２Ｅの各判定結果値を用い、その組み合わせを元に液面乖離判定を行い（ステップＳ１３）、液面乖離の判定結果値を出力する（ステップＳ１４）。この液面乖離判定は、後述する図６の液面乖離及び乖離要因判定テーブル４１を参照して行われる。
【００５３】
続いて、第２制御部８は、ステップＳ１４の液面乖離判定結果値が正常又は異常のいずれであるかを判定する（ステップＳ１５）。ステップＳ１５で、正常であると判定された場合（ステップＳ１５の判定が「正常」の場合）、第２処理部７及び第２制御部８は、図５の液面乖離判定及び乖離要因判定の処理を終了する。
【００５４】
一方、液面乖離判定結果値が異常である場合には（ステップＳ１５の判定が「異常」の場合）、第２制御部８は、泡接触判定処理部７２Ａ、第１液面乖離判定処理部７２Ｂ〜第４液面乖離判定処理部７２Ｅの各判定結果値を用い、乖離要因の判定を行う（ステップＳ１６）。そして、第２制御部８は、乖離要因の判定結果値を出力する（ステップＳ１７）。この乖離要因判定は、後述する図６の液面乖離及び乖離要因判定テーブル４１を参照して行われる。
【００５５】
次に、第２制御部８は、分注プローブ１ａの先端部１ｅと液面が乖離しているか否かの判定結果、及び、乖離要因のステータスを、表示部２５の画面に表示することによりユーザーに通知する（ステップＳ１８）。このとき、液面乖離判定が‘Ｆａｉｌ’であれば、ユーザーにその要因に対する推奨対処手順を通知する。本処理が終了すると、第２処理部７及び第２制御部８は、図５の液面乖離判定及び乖離要因判定の処理を終了する。
【００５６】
＜液面乖離及び乖離要因判定テーブル＞
図６は、泡接触判定処理、第１の判定処理〜第４の判定処理の各判定結果及び乖離要因の組み合わせを示した液面乖離及び乖離要因判定テーブルの一例を示す図である。
【００５７】
図６の液面乖離及び乖離要因判定テーブル４１は、「泡接触判定」、「第１の判定」、「第２の判定」、「第３の判定」、「第４の判定」、「液面乖離判定 正常or異常」及び「乖離要因第１候補」の各項目を有する。図６に示す液面乖離及び乖離要因判定テーブル４１では、これらの各判定処理結果の組み合わせごとにパターン分類され、分類されたパターンと液面乖離要因とが対応付けられている。図６に示す液面乖離及び乖離要因判定テーブル４１には、「パターン０」〜「パターン３１」の３２種類のパターンが登録されている。
【００５８】
乖離要因は、「接触」、「泡」、「静電気」の３つに分類される。図６の組み合わせ表は、自動分析装置１０の内部（例えば不揮発性ストレージ２７）に予め記憶されており、液面乖離判定及び乖離要因判定の際に参照される。乖離要因は、自動分析装置１０を用いて実験により得られた結果から、それぞれの組み合わせパターンの発生する割合を求め、その中で最も割合の高いものを最も蓋然性の高い要因として決定する。よって、乖離要因と組み合わせパターンの対応関係は、自動分析装置１０の構成や判定ロジックの閾値によって変化し得る。また、これらの乖離要因はあくまで蓋然性が高いと推定されるものであり、決定事項ではない。よって、他の乖離要因も列挙しつつ、まず優先的に確認するべき要因として表示部２５の画面に表示し、ユーザーに通知する。
【００５９】
例えば、パターン０は、泡接触判定、第１の判定、第２の判定、第３の判定、及び第４の判定の全ての処理結果が‘Ｐａｓｓ’の場合であり、液面乖離判定結果は‘正常’、液面乖離要因は‘該当なし’である。また、パターン１は、泡接触判定の処理の結果が‘Fａｉｌ’であり、第１の判定〜第４の判定の各処理の結果が‘Ｐａｓｓ’である場合であり、液面乖離判定結果は‘異常’、液面乖離要因は‘泡’である。また、パターン３は、泡接触判定の処理の結果が‘Ｆａｉｌ’であり、第１の判定〜第３の判定の各処理の結果が‘Ｐａｓｓ’であり、第４の判定の処理結果が‘Ｆａｉｌ’である場合である。パターン３では、液面乖離判定結果は‘異常’であり、液面乖離要因は‘静電気’である。これらの各乖離要因に関連付けて、推奨対処手順をユーザーに通知するとよい。
【００６０】
＜乖離要因判定結果と推奨対処手順＞
図７は、乖離要因判定結果に対する推奨対処手順が登録された推奨対処手順テーブルの一例を示す図である。
【００６１】
図７の推奨対処手順テーブル４２は、「乖離要因判定結果」及び「モニターに表示する推奨対処手順」の項目を有する。推奨対処手順テーブル４２には、乖離要因の‘接触’に対し、モニターに表示する推奨対処手順として、‘容器が正しくセットされているか確認してください。試料が少ない場合は小容量の容器に移しかえてください。’という文例が示されている。
【００６２】
また、乖離要因の‘泡’に対し、推奨対処手順として‘試料表面の泡を取り除いてください。’という文例が示されている。
【００６３】
また、乖離要因の‘静電気’に対し、推奨対処手順として‘試料容器を湿らせた布で拭いてください’という文例が示されている。
【００６４】
さらに、乖離要因の‘不明’に対し、推奨対処手順として‘接触、泡、静電気が発生していないか確認してください’という文例が示されている。不明と判定される場合は、５つの判定結果があり得ない組み合わせ等の場合である。
【００６５】
＜第２処理部の動作の詳細＞
次に、図５のステップＳ１１のサブルーチンである、第２処理部７による泡接触判定処理、第１の判定処理〜第４の判定処理の各処理の詳細を説明する。
【００６６】
［泡接触判定処理］
まず、泡接触判定処理について説明する。泡接触判定処理は、演算部７１及び泡接触判定処理部７２Ａにより実行される。
図８は、泡接触判定処理の手順の例を示すフローチャートである。
図９は、静電容量波形を処理したときのイメージを表す図である。図９の波形１００のグラフにおいて、横軸は図４と同様に時間を示し、縦軸も図４と同様に発振周波数を示す。
【００６７】
第２処理部７の演算部７１は、まず、記憶部６（図１参照）に記憶されている交流信号のデータから、静電容量波形データ（図９の波形１００）として交流信号の発振周波数の時系列データを取得する（ステップＳ２１）。次に、演算部７１は、この静電容量波形データに対して、測定開始（分注プローブ１ａの下降開始）から終了まで一定間隔の時間領域（一定区間）を設定し、特徴量として各時間領域内での発振周波数の最小値を抽出する（ステップＳ２２）。特徴量としての最小値は、図９において黒丸印で示される。そして、演算部７１は、この複数の最小値のデータを、データ系列Ａ（図９のイメージ１０１）として出力する（ステップＳ２３）。
【００６８】
次に、泡接触判定処理部７２Ａは、演算部７１から最小値のデータ系列Ａを取り込み、最小値のデータ系列Ａから、液面検知前のデータを除外する（ステップＳ２４）。液面検知前のデータとは、分注プローブ１ａの先端部１ｅが容器２内の液面に接触したと検出部５１が判定する前に、演算部７１によって抽出されていた特徴量（最小値）を示す。図９のイメージ１０１において、液面検知前のデータは、特徴量Ｆ１、Ｆ２及びＦ３である。
【００６９】
次に、泡接触判定処理部７２Ａは、不揮発性ストレージ２７等に予め登録されている検索波形データＲ１（異常波形モデルの一例、図９のイメージ１０３）を参照する（ステップＳ２５）。検索波形データＲは、処理フィルタの一例であり、本実施形態では、検索波形データＲ１を十点（−２，−１，０，１，１，１，１，０，−１，−２）から構成される波形とする。この検索波形データＲ１は、分注プローブ１ａの先端部１ｅが容器２内の液面の泡に接触した場合に観察される異常波形との適合度を計算することを意図して設計されたものである。異常波形は、段階的に上昇し、数ミリ秒等の比較的長い時間をかけて山型又は谷型に推移するという特徴を有する。なお、検索波形データＲ１が負の値を含むのは、データ系列の検索波形データＲ１が適用された特徴量の平均値が、相互相関処理の前後で同じ値になるようにする意図によるものである。
【００７０】
検索波形データＲ１の形状は、段階的に上昇し、数ミリ秒に対応する一定の期間等しい値で推移し、その後段階的に下降する形状であれば、どのような形状であってもよい。また、検索波形データＲ１を構成する関数の数も、１０点以上であってもよく、１０点以下であってもよい。また、検索波形データＲ１を構成する関数の各値も、他の値を取ってもよい。検索波形データＲ１を構成する関数の数（波形の時間方向における長さ）は、分注プローブ１ａの先端部１ｅが容器２内の液面に接触してから停止するまでの区間に対して、十分大きいものであることが重要である。
【００７１】
容器２内の液面に発生した泡に分注プローブ１ａの先端部１ｅが接触した場合に取得される静電容量波形データの波形の変動は、１ミリ秒〜１０ミリ秒のオーダーで起こることが実験的に確認されている。したがって、検索波形データＲ１の波形の時間方向における長さを、数ミリ秒程度の長さに設定することにより、液面乖離要因が‘泡’であるか否かを判断するのに十分な区間、静電容量波形データのモニタリングを行うことが可能となる。
【００７２】
次に、泡接触判定処理部７２Ａは、データ系列Ａと検索波形データＲ１とを用いて相互相関演算を行う（ステップＳ２６）。相互相関演算に用いられる相互相関関数Ｃ１（ｎ）は、例えば下記の式（１）によって表される。下記の式（１）において、Ｒ１（１）は、検索波形データＲ１の左端の設定値（関数）、Ｒ１（２）は左端から２番目の設定値、Ｒ１（１０）は右端の設定値である。Ａ（ｎ）は、データ系列Ａを構成する特徴量のうちのｎ番目の特徴量である。
【００７３】
Ｃ１（ｎ）＝Ａ（ｎ）＊Ｒ１（１）＋Ａ（ｎ＋１）＊Ｒ１（２）＋Ａ（ｎ＋２）＊Ｒ１（３）＋Ａ（ｎ＋３）＊Ｒ１（４）＋Ａ（ｎ＋４）＊Ｒ１（５）＋Ａ（ｎ＋５）＊Ｒ１（６）＋Ａ（ｎ＋６）＊Ｒ１（７）＋Ａ（ｎ＋７）＊Ｒ１（８）＋Ａ（ｎ＋８）＊Ｒ１（９）＋Ａ（ｎ＋９）＊Ｒ１（１０）・・・・・式（１）
【００７４】
泡接触判定処理部７２Ａは、上記式（１）におけるｎの値を‘１’から‘（データ系列Ａのポイント数−９）’の範囲で取り、各ｎについて相互相関演算を行う。ここでいうポイントとは、特徴量の抽出点である。ステップＳ２６の相互相関演算の結果、図９のイメージ１０４に示す相互相関関数Ｃ１（ｎ）が得られる。続いて、泡接触判定処理部７２Ａは、相互相関関数Ｃ１（ｎ）の最大値Ｃ１ｍａｘを求める（ステップＳ２７）。
【００７５】
そして、泡接触判定処理部７２Ａは、予めＲＯＭ２２等に格納されていた判定閾値Ｚａ、Ｚｂを参照し（ステップＳ２８）、相互相関関数Ｃ１（ｎ）の最大値Ｃ１ｍａｘの値が判定閾値Ｚａ以下であるか否かを判定する（ステップＳ２９）。
【００７６】
ステップＳ２９において、相互相関関数Ｃ１（ｎ）の最大値Ｃ１ｍａｘの値は判定閾値Ｚａの値より大きいと判定された場合（ステップＳ２９がＮＯ判定の場合）、泡接触判定処理部７２Ａは、液面検知が正常に行われなかった（Ｆａｉｌ）と判定する（ステップＳ３０）。そして、泡接触判定処理部７２Ａは、本フローチャートの処理を終了するとともに、図５のステップＳ１２の処理に移行する。
【００７７】
一方、ステップＳ２９において、相互相関関数Ｃ１（ｎ）の最大値Ｃ１ｍａｘの値が判定閾値Ｚａ以下であると判定された場合（ステップＳ２９がＹＥＳ判定の場合）、泡接触判定処理部７２Ａは、相互相関関数Ｃ１（ｎ）の最大値Ｃ１ｍａｘの値が判定閾値Ｚｂ以上であるか否かを判定する（ステップＳ３１）。ステップＳ３１で、相互相関関数Ｃ１（ｎ）の最大値Ｃ１ｍａｘの値は判定閾値Ｚｂより小さいと判定された場合（ステップＳ３１がＮＯ判定の場合）、泡接触判定処理部７２Ａは、液面検知が正常に行われなかった（Ｆａｉｌ）と判定する（ステップＳ３０）。
【００７８】
一方、ステップＳ３１で、相互相関関数Ｃ１（ｎ）の最大値Ｃ１ｍａｘの値は判定閾値Ｚｂより小さいと判定された場合（ステップＳ３１がＮＯ判定の場合）、泡接触判定処理部７２Ａは、液面検知が正常に行われなかった（Ｆａｉｌ）と判定する（ステップＳ３０）。一方、ステップＳ３１で、相互相関関数の最大値Ｃ１ｍａｘの値は判定閾値Ｚｂ以上であると判定された場合（ステップＳ３１がＹＥＳ判定の場合）、泡接触判定処理部７２Ａは、液面検知は正常に行われた（Ｐａｓｓ）と判定する（ステップＳ３２）。そして、泡接触判定処理部７２Ａは、本フローチャートの処理を終了するとともに、図５のステップＳ１２の処理に移行する。
【００７９】
つまり、泡接触判定処理部７２Ａは、相互相関関数の最大値Ｃ１ｍａｘの値が判定閾値Ｚｂ以上であり、かつ、判定閾値Ｚａ以下である場合に、液面検知は正常に行われた（Ｐａｓｓ）と判定する。相互相関関数の最大値Ｃ１ｍａｘの値が判定閾値Ｚｂより小さい場合、又は、判定閾値Ｚａより大きい場合には、液面検知は正常に行われなかった（Ｆａｉｌ）と判定する。
【００８０】
なお、液面検知が正常に行われなかったと判定した場合には、泡接触判定処理部７２Ａは、判定の結果を、表示部９を介してユーザーへ通知してもよい。
【００８１】
次に、演算部７１及び泡接触判定処理部７２Ａによる静電容量波形データに対するデータ処理について、具体例を挙げて説明する。説明は、液面検知が正常に行われた場合における静電容量波形（以下、「正常波形」と称する）、及び、液面検知が正常に行われなかった場合における静電容量波形（以下、「異常波形」と称する）の両方について行う。
【００８２】
（正常波形に対するデータ処理例）
図１０は、正常波形に対するデータ処理及び判定処理を示す図である。
【００８３】
正常波形の特徴を、図１０の一番上に波形（正常波形）１００として示す。正常波形１００の特徴は、静電容量値（ＣＲ発振回路４の発振周波数として検出）が一過性に上昇し、その後は小さい幅の変動に収まる点である（図４の波形３１も同様）。正常波形１００がこのような特徴を有する第１の理由は、分注プローブ１ａが下降する空気中での静電容量の値と、先端部１ｅが液面に接触している時の静電容量の値とが大きく違うことにある。第２の理由は、先端部１ｅが空気中及び液面に接触した後の各状態における静電容量値が安定していることにある。
【００８４】
正常波形１００は、小さい幅の変動に収まった後、すなわち、液面を検知した後の波形は平坦な形状となる。よって、泡接触判定処理部７２Ａによる相互相関演算の結果は、ほぼ０に近い値となる。それゆえ、相互相関関数Ｃ１（ｎ）の最大値Ｃ１ｍａｘの値も一定となり、最大値Ｃ１ｍａｘは判定閾値Ｚｂ及びＺａにより規定される範囲内に収まる。したがって、泡接触判定処理部７２Ａによる泡接触判定の結果は‘ＰＡＳＳ’となる。
【００８５】
図１１及び図１２は、分注プローブ１ａの先端部１ｅが泡に接触して停止したときに起こり得る静電容量波形、データ処理及び判定処理の例を示す図である。
【００８６】
容器２内には、生体試料や臨床検査用試薬等の液体（試料３）が入っている。生体試料は、その成分により泡立ちが起きることが知られており、臨床検査用試薬は、界面活性剤を含有する場合に泡立ちが発生することが知られている。泡の表面に接触した状態で分注プローブ１ａが停止した場合、停止時の振動により泡と分注プローブ１ａの先端部１ｅとの接触面積が変わることがある。また、分注プローブ１ａの停止時の振動が泡に伝わることにより泡に振動が発生し、泡の形状が変化することがある。
【００８７】
生体試料又は臨床検査用試薬は一般的に電解質成分を含むため、泡が分注プローブ１ａの先端部１ｅに触れることは、分注プローブ１ａと泡とが導通することを意味する。そして、泡との接触面積、又は、泡の形状が変化することは、泡を通して分注プローブ１ａと導通する試料部分の面積が絶えず変化することを意味する。これは、分注プローブ１ａ周縁部（例えば装置筐体のグラウンド）に対して、分注プローブ１ａを含んだ導通物質全体の面積が変化するということであり、二つの相対する導体間の静電容量は、それら導体の面積の変化に伴い変動する。したがって、分注プローブ１ａの先端部１ｅに泡が接触した場合、静電容量値はランダムに増減する。
【００８８】
この場合の静電容量値の変化は、１ミリ秒から１０ミリ秒のオーダーで観察されるものであり、液面接触による静電容量の急激な上昇時における変動と比較すると、緩やかな変動となる。図１１には、静電容量値が山形に変動する例を示し、図１２には、静電容量値が谷形に変動する例を示した。なお、図１３には、本発明の実施形態に近い構成の液面検知機構を用いて、実際に泡に対する液面検知を行った場合に観測された波形の実例を示す。図１３Ａは、山形波形の実例を示す図であり、図１３Ｂは、谷形波形の実例を示す図である。
【００８９】
図１１に示すような山形の変動においては、相互相関関数（検索波形データＲ１）とデータ系列Ａとはよく適合するため、相互相関関数Ｃ１（ｎ）の最大値Ｃ１ｍａｘの値は、液面検知が正常に行われた場合と比較して高くなる。これにより、相互相関関数の最大値Ｃ１ｍａｘの値は閾値Ｚａを超えるため、泡接触判定の結果は‘Ｆａｉｌ’となる。
【００９０】
図１２に示すような谷形の変動においては、相互相関関数（検索波形データＲ１）とデータ系列Ａとは適合しない。しかしながら、谷形は山形を逆転したものであるため、相互相関演算によって演算される値は、山形の波形を基準値から負方向に反転させた数値となる。つまり、相互相関関数の最大値Ｃ１ｍａｘの値は、液面検知が正常に行われた場合と比較して、顕著に低い値となる。これにより、相互相関関数の最大値Ｃ１ｍａｘの値は閾値Ｚｂ未満となり、泡接触判定の結果は‘Ｆａｉｌ’となる。
【００９１】
このように、泡接触判定処理では、演算部７１が、液面検出前から検出後まで適当な時間間隔で静電容量値をサンプリングする。続いて、演算部７１は、サンプリングしたデータ（時系列データ）について、一定数のまとまりごとにそれらの特徴量を抽出し、その特徴量を元にデータ系列Ａを作成及び出力する。続いて、泡接触判定処理部７２Ａは、データ系列Ａについて、異常波形モデルとの相互相関を計算し、相互相関の計算結果に基づいて液面の検知が正常に行われたかどうかを判定する。異常波形モデルは、分注プローブ１ａの先端部１ｅが容器２内の液面の泡に接触した場合に観察される波形に基づき生成された波形モデルである。この泡接触判定処理により、静電容量波形の変化の形が、異常波形モデルに近いかどうかがわかるため、液面乖離要因が泡である場合の誤検出の判定精度を高めることができる。
【００９２】
なお、上述した実施形態では、時系列データＡの特徴量として最小値を抽出する例を説明したが、最大値を抽出してもよい。
【００９３】
［第１の判定処理］
次に、第１の判定処理について説明する。第１の判定処理は、演算部７１及び第１液面乖離判定処理部７２Ｂにより実行される。
図１４は、第１の判定処理の手順の例を示すフローチャートである。
図１５は、静電容量波形を処理したときのイメージを表す図である。図１５の波形１００のグラフにおいて、図４と同様に横軸は時間、縦軸は発振周波数を示す。
【００９４】
第２処理部７の演算部７１は、まず、第２処理部７の演算部７１は、まずステップＳ３１〜Ｓ３３の処理を行い、最小値のデータ系列Ａを作成する。このステップＳ３１〜Ｓ３３の処理は、図８のステップＳ２１〜Ｓ２３の処理と同じであるため、詳細な説明を省略する。
【００９５】
次に、第１液面乖離判定処理部７２Ｂは、演算部７１から最小値のデータ系列Ａを取り込み、データ系列Ａについて微分演算を行う（ステップＳ３４）。すなわち、第１液面乖離判定処理部７２Ｂは、データ系列Ａの各点に対して下記式（２）の関数Ｂ（ｎ）を適用し、データ系列Ａの微分値からなるデータ系列Ｂ（図１５のイメージ１０２）を得る（ステップＳ３５）。下記式（２）において、‘ｎ’は自然数である。
【００９６】
Ｂ（ｎ）＝Ａ（ｎ＋１）−Ａ（ｎ）      ・・・・・式（２）
【００９７】
次に、第１液面乖離判定処理部７２Ｂは、不揮発性ストレージ２７等に予め登録されている検索波形データＲ２（図１５のイメージ１０３）を参照する（ステップＳ３６）。検索波形データＲ２は、処理フィルタの一例でもあり、本実施形態では、検索波形データＲ２を三点（−１，２，−１）から構成される波形とする。この検索波形データＲ２は、数値が一過性に上昇する波形、つまりスパイク状の波形との適合度を計算することを意図して設計されたものである。なお、検索波形データＲ２が負の値を含むのは、データ系列の検索波形データＲ２が適用された特徴量の平均値が、相互相関処理の前後で同じ値になるようにする意図によるものである。
【００９８】
検索波形データＲ２は、スパイク状の波形、つまりある一点とその両脇の点との差が相対的に十分に大きい形状であれば、どのような形状であってもよい。また、検索波形データＲ２を構成する関数の数も、３点以上であってもよく、３点以下であってもよい。また、検索波形データＲ２を構成する関数の各値も、他の値を取ってもよい。
【００９９】
データ系列Ｂと検索波形データＲ２との相互相関関数Ｃ２（ｎ）は、例えば下記の式（３）で表される。下記式（３）において、Ｒ２（１）は、検索波形データＲの３点の左側の設定値、Ｒ２（２）は３点の中央の設定値、Ｒ２（３）は３点の右側の設定値である。
【０１００】
Ｃ２（ｎ）＝Ｂ（ｎ）＊Ｒ２（１）＋Ｂ（ｎ＋１）＊Ｒ２（２）＋Ｂ（ｎ＋２）＊Ｒ２（３）      ・・・・・式（３）
【０１０１】
第１液面乖離判定処理部７２Ｂは、ｎの値を‘１’から‘（データ系列Ｂのポイント数−２）’の範囲で取り、各ｎについて相互相関演算を行う（ステップＳ３７）。ステップＳ３７の相互相関演算の結果、図１５のイメージ１０４に示す相互相関関数Ｃ２（ｎ）が得られる。続いて、第１液面乖離判定処理部７２Ｂは、相互相関関数Ｃ２（ｎ）の最大値Ｃ２ｍａｘを求める（ステップＳ３８）。
【０１０２】
続いて、第１液面乖離判定処理部７２Ｂは、予めＲＯＭ２２等に格納されていた判定閾値Ｚｃを参照し（ステップＳ３９）、相互相関関数の最大値Ｃ２ｍａｘが判定閾値Ｚｃ以上であるか否かを判定する（ステップＳ４０）。
【０１０３】
ステップＳ４０において、相互相関関数の最大値Ｃ２ｍａｘの値が判定閾値Ｚｃの値以上である場合（Ｓ４０がＹＥＳ判定の場合）、第１液面乖離判定処理部７２Ｂは、液面検知が正常に行われた（Ｐａｓｓ）と判定する（ステップＳ４１）。そして、本フローチャートの処理を終了するとともに、図５のステップＳ１２に移行する。
【０１０４】
一方、ステップＳ４０において、相互相関関数の最大値Ｃ２ｍａｘの値は判定閾値Ｚｃより小さいと判定された場合（ステップＳ４０がＮＯ判定の場合）、第１液面乖離判定処理部７２Ｂは、液面検知が正常に行われなかった（Ｆａｉｌ）と判定する（ステップＳ４２）。そして、本フローチャートの処理を終了するとともに、図５のステップＳ１２に移行する。
【０１０５】
ここで、液面検知が正常に行われなかったと判定した場合には、第１液面乖離判定処理部７２Ｂは、判定結果を、表示部９を介してユーザーに通知してもよい。
【０１０６】
このように、第１の判定処理では、液面検出前から検出後まで適当な時間間隔で静電容量値をサンプリングする。続いて、サンプリングしたデータ（時系列データ）について、一定数のまとまりごとにそれらの特徴量を抽出し、その特徴量を元に新たなデータ系列を作成及び記憶する。そして、その新たなデータ系列について、正常な波形のモデル（例えばスパイク状波形）との相互相関を求め、その結果から液面を正しく検知できたか否かを判定する。この第１の判定処理により、静電容量波形の変化の形が、正常な波形のものに近いかどうかがわかる。
【０１０７】
なお、時系列データの特徴量として最小値を抽出する例を説明したが、最大値を抽出してもよい。また、データ系列は、時系列データの特徴量の微分演算により作成される。
【０１０８】
［第２の判定処理］
次に、第２の判定処理について説明する。第２の判定処理は、演算部７１及び第２液面乖離判定処理部７２Ｃにより実行される。
図１６は、第２の判定処理の手順の例を示すフローチャートである。
図１７は、静電容量波形を処理したときのイメージを表す図である。
【０１０９】
第２処理部７の演算部７１は、まず、ステップＳ５１〜Ｓ５３の処理を行い、最小値のデータ系列Ａを作成する。このステップＳ５１〜Ｓ５３の処理は、図８のステップＳ２１〜Ｓ２３の処理と同じであるため、詳細な説明を省略する。
【０１１０】
続いて、第２液面乖離判定処理部７２Ｃは、ステップＳ５４〜Ｓ５５の処理を行い、微分値のデータ系列Ｂを作成する。このステップＳ５４〜Ｓ５５の処理は、図１４のステップＳ３４及びＳ３５の処理と同様であるため、詳細な説明を省略する。
【０１１１】
次に、第２液面乖離判定処理部７２Ｃは、微分値の関数Ｂ（ｎ）の最大値Ｂｍａｘ（図１７のイメージ１０２）を求める（ステップＳ５６）。最大値Ｂｍａｘは、微分値のデータ系列Ｂによる波形の最も大きな傾きの大きさに相当する。そして、第２液面乖離判定処理部７２Ｃは、予めＲＯＭ２２等に格納されていた判定閾値Ｚｄを参照し（ステップＳ５７）、微分値の最大値Ｂｍａｘの値は判定閾値Ｚｄ以上であるか否かを判定する（ステップＳ５８）。
【０１１２】
ステップＳ５８において、微分値の最大値Ｂｍａｘの値は判定閾値Ｚｄの値以上であると判定された場合（ステップＳ５８がＹＥＳ判定の場合）、第２液面乖離判定処理部７２Ｃは、液面検知が正常に行われた（Ｐａｓｓ）と判定する（ステップＳ５９）。そして、本フローチャートの処理を終了するとともに、図５のステップＳ１２に移行する。
【０１１３】
一方、ステップＳ５８において、微分値の最大値Ｂｍａｘの値は判定閾値Ｚｄの値より小さいと判定された場合（ステップＳ５８がＮＯ判定の場合）、第２液面乖離判定処理部７２Ｃは、液面検知が正常に行われなかった（Ｆａｉｌ）と判定する（ステップＳ６０）。そして、本フローチャートの処理を終了するとともに、図５のステップＳ１２に移行する。
【０１１４】
ここで、液面検知が正常に行われなかったと判定された場合には、第２液面乖離判定処理部７２Ｃは、判定結果を、表示部９を介してユーザーに通知してもよい。
【０１１５】
このように、第２の判定処理では、液面検出前から検出後まで適当な時間間隔で静電容量値をサンプリングする。続いて、サンプリングしたデータ（時系列データ）について、一定数のまとまりごとにそれらの特徴量を抽出し、その特徴量を元に新たなデータ系列を作成及び記憶する。そして、その新たなデータ系列の最大値を求め、その結果から液面を正しく検知できたか否かを判定する。この第２の判定処理により、静電容量波形の変化の大きさ（急峻さ）が、正常であるか否かがわかる。第２の判定処理は、主として接触の判定に有効であると考えられるが、最終的な乖離要因は５つの判定結果の組み合わせに基づいて判定される。
【０１１６】
なお、時系列データの特徴量として最小値を抽出する例を説明したが、最大値を抽出してもよい。また、データ系列は、時系列データの特徴量の微分演算により作成される。
【０１１７】
［第３の判定処理］
次に、第３の判定処理について説明する。第３の判定処理は、演算部７１及び第３液面乖離判定処理部７２Ｄにより実行される。
図１８は、第３の判定処理の手順の例を示すフローチャートである。 
図１９は、静電容量波形を処理したときのイメージを表す図である。
【０１１８】
第２処理部７の演算部７１は、まず静電容量波形データ（図１９の波形１１０）に対してステップＳ７１〜Ｓ７３の処理を行い、最小値のデータ系列Ａを作成する。このステップＳ７１〜Ｓ７３の処理は、図８のステップＳ２１〜Ｓ２３の処理と同じであるため、詳細な説明を省略する。
【０１１９】
次に、第３液面乖離判定処理部７２Ｄは、予め記憶部６に格納されていたデータ処理のための閾値Ｗを参照し、データ系列Ａの一番目のポイントから数えて最初に閾値Ｗを超えるｓ番目のポイントを検索する（ステップＳ７４，Ｓ７５）。言い換えると、データ系列Ａにおいて最初に閾値Ｗを超えたポイントをｓとする。ここでポイントとは、特徴量の抽出点である。もし最初の検索でｓが見つからない場合、つまり、データ系列Ａのポイントが一度も閾値Ｗを超えない場合には、ｓをデータ系列Ａの終了ポイント（最終ポイント）とする。
【０１２０】
続いて、第３液面乖離判定処理部７２Ｄは、データ系列Ａのｓ以降のポイントで最初に閾値Ｗ以下となるポイントを検索し、そのポイントをｔとする（ステップＳ７６，Ｓ７７）。もし閾値Ｗ以下となるポイントが存在しない場合、つまりデータ系列Ａのポイントが閾値Ｗを超えて以降、一度も閾値Ｗ以下に下がらない場合には、ｔは終了ポイントに１を加算した値とする。
【０１２１】
次に、第３液面乖離判定処理部７２Ｄは、データ系列Ａの閾値Ｗを超えた区間Ｃの長さを、式（４）により計算する（ステップＳ７８）。区間Ｃの長さは、データ系列Ａの閾値Ｗを超えたポイント数として示される（図１９のイメージ１１１）。
【０１２２】
Ｃ＝ｔ−ｓ＋１      ・・・・・式（４）
【０１２３】
そして、第３液面乖離判定処理部７２Ｄは、予めＲＯＭ２２等に格納されていた判定閾値Ｚｅを参照し（ステップＳ７９）、区間Ｃの長さが判定閾値Ｚｅ以上であるか否かを判定する（ステップＳ８０）。この判定閾値Ｚｅは、分注プローブ１ａの先端部１ｅが液面に接触してから分注プローブ１ａが上昇するまでの時間に基づいて設定される。
【０１２４】
ステップＳ８０において、区間Ｃの長さは判定閾値Ｚｃの値以上であると判定された場合（ステップＳ８０がＹＥＳ判定の場合）、第３液面乖離判定処理部７２Ｄは、液面検知が正常に行われた（Ｐａｓｓ）と判定する（ステップＳ８１）。そして、本フローチャートの処理を終了するとともに、図５のステップＳ１２に移行する。
【０１２５】
一方、ステップＳ８０において、区間Ｃの長さは判定閾値Ｚｃの値より小さいと判定された場合（ステップＳ８０がＮＯ判定の場合）、第３液面乖離判定処理部７２Ｄは、液面検知が正常に行われなかった（Ｆａｉｌ）と判定する（ステップＳ８２）。そして、本フローチャートの処理を終了するとともに、図５のステップＳ１２に移行する。
【０１２６】
ここで、液面検知が正常に行われなかったと判定した場合には、第３液面乖離判定処理部７２Ｄは、判定結果を、表示部９を介してユーザーに通知してもよい。
【０１２７】
このように、第３の判定処理では、液面検出前から検出後まで適当な時間間隔で静電容量値をサンプリングする。続いて、サンプリングしたデータ（時系列データ）について、予め記憶された固有の閾値Ｗを最初に超える（又は達した）ポイントと、それ以降に閾値Ｗ以下へ下がったポイントを抽出し、その二つのポイントの間にあるポイント数を、閾値Ｗを超えた区間Ｃの長さとして求める。そして、その結果から液面を正しく検知できたか否かを判定する。この第３の判定処理により、静電容量波形の形状が正常である（最小値が一定ポイント数以上連続する）か否かがわかる。第３の判定処理は、主として泡の判定に有効であると考えられるが、最終的な乖離要因は５つの判定結果の組み合わせに基づいて判定される。
【０１２８】
なお、時系列データの特徴量として最小値を抽出する例を説明したが、最大値を抽出してもよい。
【０１２９】
［第４の判定処理］
次に、第４の判定処理について説明する。第４の判定処理は、演算部７１及び第４液面乖離判定処理部７２Ｅにより実行される。
図２０は、第４の判定処理の手順の例を示すフローチャートである。
図２１は、静電容量波形を処理したときのイメージを表す図である。
【０１３０】
第２処理部７の演算部７１は、まず記憶部６に記憶されている交流信号のデータから、静電容量波形データ（図２１の波形１００）として交流信号の発振周波数の時系列データを取得する（ステップＳ９１）。次に、演算部７１は、この静電容量波形データに対して、測定開始（分注プローブ１ａの下降開始）から終了まで一定間隔の時間領域（一定区間）を設定し、特徴量として各時間領域内での発振周波数の最大値（図２１に白抜き丸印で示す）及び最小値（図２１に黒丸印で示す）を抽出する（ステップＳ９２）。そして、演算部７１は、は、この最大値及び最小値のデータを、データ系列Ａａ及びデータ系列Ａｂ（図２１のイメージ１２１）として出力する（ステップＳ９３，Ｓ９４）。
【０１３１】
次に、第４液面乖離判定処理部７２Ｅは、最大値のデータ系列Ａａの最大値ＭＡＸ（Ａａ）及び最小値のデータ系列Ａｂの最大値ＭＡＸ（Ａｂ）を演算する（ステップＳ９５）。続いて、第４液面乖離判定処理部７２Ｅは、最大値ＭＡＸ（Ａａ）から最大値ＭＡＸ（Ａｂ）を減算した数値Ｄ（差分）を求める（図２１のイメージ１２２）（ステップＳ９６）。
【０１３２】
原理的に最大値のデータ系列Ａａの最大値ＭＡＸ（Ａａ）は、必ず最小値のデータ系列Ａｂの最大値ＭＡＸ（Ａｂ）以上となるので、差分Ｄは‘０’以上の値となる。
【０１３３】
そして、第４液面乖離判定処理部７２Ｅは、予めＲＯＭ２２等に格納されていた判定閾値Ｚｆを参照し（ステップＳ９７）、差分Ｄが判定閾値Ｚｆ以下であるか否かを判定する（ステップＳ９８）。
【０１３４】
ステップＳ９８において、差分Ｄは判定閾値Ｚｆの値以下であれると判定された場合（ステップＳ９８がＹＥＳ判定の場合）、第４液面乖離判定処理部７２Ｅは、液面検知が正常に行われた（Ｐａｓｓ）と判定する（ステップＳ９９）。そして、本フローチャートの処理を終了するとともに、図５のステップＳ１２に移行する。
【０１３５】
一方、ステップＳ９８において差分Ｄは判定閾値Ｚｆの値より大きいと判定された場合（ステップＳ９８がＮＯ判定の場合）、第４液面乖離判定処理部７２Ｅは、液面検知が正常に行われなかった（Ｆａｉｌ）と判定する（ステップＳ１００）。そして、本フローチャートの処理を終了するとともに、図５のステップＳ１２に移行する。
【０１３６】
ここで、液面検知が正常に行われなかったと判定した場合には、第４液面乖離判定処理部７２Ｅは、判定結果を、表示部９を介してユーザーに通知してもよい。
【０１３７】
このように、第４の判定処理では、液面検出前から検出後まで適当な時間間隔で静電容量値をサンプリングする。続いて、サンプリングしたデータ（時系列データ）について、一定数のまとまりごとにそれらの特徴量として最大値及び最小値を抽出し、最大値からなるデータ系列と最小値からなるデータ系列を作成及び記憶する。そして、この２つのデータ系列の最大値を求め、それらの差を演算し、その結果から液面を正しく検知できたか否かを判定する。この第４の判定処理により、静電容量波形の安定性が高いか否か（ノイズの有無）がわかる。第４の判定処理は、主として静電気の判定に有効であると考えられるが、最終的な乖離要因は５つの判定結果の組み合わせに基づいて判定される。
【０１３８】
＜各種効果＞
上述した実施形態では、演算部７１が、分注プローブ１ａが下降を開始してから一定時間が経過するまでの間に発振回路が出力した、交流信号の発振周波数の時系列データを一定区間ごとに分け、一定区間ごとの特徴量を抽出し、抽出した特徴量をデータ系列Ａとして出力する。また、泡接触判定処理部７２Ａが、演算部７１から出力された特徴量のデータ系列Ａの波形と、分注プローブ１ａの先端部１ｅが容器２内の液面の泡に接触した場合に観察される波形に基づく異常波形モデル（検索波形データＲ１）との相互相関を計算し、相互相関の計算結果に基づいて液面の検知が正常に行われたかどうかを判定する。さらに、第２制御部８が、泡接触判定処理部７２Ａの判定結果値を元に液面乖離判定を行い、液面乖離の判定結果値を出力する。それゆえ、本実施形態によれば、液面乖離要因が泡である場合の誤検出の判定精度を高めることができる。
【０１３９】
また、本実施形態によれば、静電容量方式の液面検知の誤検知に対して、５つの判定処理の判定結果の組み合わせに基づいて、検知結果が誤検知であったことを判定することができる。すなわち、分注プローブ１ａの先端部１ｅが液面に接触していない（液面乖離）にも関わらず検出部５１が液面を検知した場合に、その検知結果が誤りであることをより高い精度で判定することができる。
【０１４０】
また、本実施形態によれば、液面の誤検知に対して、「泡」、「静電気」、「接触」又は「不明」のいずれの要因によって液面乖離が起こったかを、５つの判定処理の判定結果の組み合わせることにより蓋然性高く推定することができる。
【０１４１】
また、本実施形態では、液面検知の誤検知（液面乖離）が発生したことが表示部９を介してユーザーに通知することにより、ユーザーによる試料の誤った成分濃度の報告を防止できる。また、本実施形態によれば、ユーザーは、試料の分析結果が、分注プローブ１ａが試料を適切に吸引していない状態で測定が実施されたものであることを把握できる。
【０１４２】
また、本実施形態によれば、液面検知の誤検知の要因及び試料等に対してどのような処理をすべきか（推奨対処手順）が表示部９に出力される。これにより、ユーザーは、表示部９に表示された内容を参照して、誤検知の要因に対する対処として蓋然性の高い要因に対する対処から実行することができる。それゆえ、本実施形態によれば、ユーザーが誤検知の要因に対する対処にかける時間を短縮できる。つまり、異常試料を再検査するまでの時間を短縮することができる。これにより、自動分析装置１０から早期に誤検知の原因を取り除いて運転を再開することが可能となるので、ユーザーの利便性を高めることができる。
【０１４３】
また、本実施形態では、静電容量方式で用いられる発振回路（例えばＣＲ発振回路４）から出力される、分注プローブ１ａの先端部１ｅと周辺部との間の静電容量に基づく交流信号のデータが記憶部６に記憶される。そして、記憶部６に記憶された交流信号のデータを元に発振周波数（波形）が解析されるため、静電容量値を演算する回路（例えば第１処理部５）に変更を加える必要がなく、別途設置する演算処理部（第２処理部７）で判定を行うことができる。第２処理部７としては、例えばＰＣなどのオペレーション・コンソールを用いることができる。このことは、閾値などの条件をソフトウェア上で簡易に変えられることを意味し、本実施形態に係る自動分析装置１０は柔軟性がある。それゆえ、本実施形態に係る自動分析装置１０は、実用性が高い。
【０１４４】
また、本実施形態によれば、既存の静電容量方式の自動分析装置に、第２処理部７及び第２制御部８を追加するだけで、液面検知の誤検知を高い精度で判定することができる。例えば、第２処理部７及び第２制御部８をソフトウェアで実現した場合には、本実施形態の機能を既存の自動分析装置に簡単に追加することが可能である。
【０１４５】
さらに、第２処理部７の泡接触判定処理部７２Ａ、第１液面乖離判定処理部７２Ｂ〜第４液面乖離判定処理部７２Ｅによる液面乖離判定の処理方法は、それぞれ極簡単な計算式（例えば式（１）〜（４））を用いて実現できるため、処理負荷が小さく、かつ演算処理にかかる時間が短くて済む。
【０１４６】
＜各種変形例＞
上述した実施形態では、泡接触判定処理部７２Ａ、第１液面乖離判定処理部７２Ｂ（第１の判定処理）〜第４液面乖離判定処理部７２Ｅ（第４の判定処理）による５つの判定結果の組み合わせに基づいて、図５の液面乖離判定及び乖離要因判定を行った。しかしながら、本発明はこれに限定されない。
【０１４７】
例えば、第１の判定処理〜第４の判定処理は行わずに、泡接触判定処理のみを行ってもよい。また、第１の判定処理及び第２の判定処理に代えて、泡接触判定処理を行ってもよい。
【０１４８】
また、例えば、泡接触判定処理、第１の判定処理〜第４の判定処理の５つの判定結果のうち２以上の判定結果の組み合わせに基づいて、液面乖離判定及び乖離要因判定を行うようにしてもよい。
【０１４９】
例えば、まず、第１の判定処理〜第４の判定処理を行い、各判定処理による液面乖離判定結果が‘異常’であり、かつ、乖離要因が‘泡’又は‘不明（Ｎ／Ａ）’と判定された場合に、泡接触判定処理を行ってもよい。このような処理を行うことにより、液面乖離要因が泡である場合の誤検出の判定精度を高めることができる。
【０１５０】
また、例えば、判定対象とする乖離要因に応じて、泡接触判定処理、第１の判定処理〜第４の判定処理の中から使用する組み合わせを決定してもよい。例えば、泡及び静電気は要因として特徴的であるので、第３の判定処理と第４の判定処理との組み合わせにより判定してもよい。あるいは、第３の判定処理及び第４の判定処理に他の判定処理を組み合わせて、液面乖離判定及び乖離要因判定をしてもよい。このように使用する判定処理（判定結果）の組み合わせを変えた場合には、それに合わせて液面乖離及び乖離要因判定テーブル４１の内容を変更する。
【０１５１】
また、第１の判定処理及び第２の判定処理はともに、微分演算を行い微分値のデータ系列に基づいて判定を行っており、相互に類似する処理が含まれるため、いずれか一方を使用又は双方を集約してもよい。
【０１５２】
また、第３の判定処理は、静電容量波形の形状を、閾値Ｗに相当する一定時間以上の時間幅で見ているという点で、第１の判定処理から第３の判定処理までとは方向性が異なる。そこで、判定結果の組み合わせに、少なくとも第３の判定処理の判定結果を含むことが望ましい。それにより、判定結果の組み合わせが多様化し、乖離要因の精度が上がる。
【０１５３】
また、上述した実施形態において、特徴量として静電容量値（発振周波数として検出）の時系列データの一定のポイント数（一定区間）ごとに最小値を抽出する代わりに、平均値や最大値を抽出する方法が考えられる。
【０１５４】
また、上述した実施形態において、一定のポイント数ごとの特徴量を抽出する際に、ポイントは必ずしも時間的に等間隔に得られたものである必要はない。つまり、可変かつ適正な時間間隔で都度サンプリングされたデータから、一定のポイント数ごと若しくは一定の時間領域ごとの最小値等を抽出する方法であっても、原理的には検出性能は変わらない（すなわち、判定結果に影響しない）。
【０１５５】
また、上述した実施形態において、誤った液面検知がなされたと判定した場合のユーザーへの通知は、スピーカーにより警報を鳴らす、画面上で警告メッセージを表示する、画面上の測定結果にマークを付加するなどの方法が考えられる。
【０１５６】
さらに、本発明は上述した実施形態例に限られるものではなく、特許請求の範囲に記載した本発明の要旨を逸脱しない限りにおいて、その他種々の応用例、変形例を取り得ることは勿論である。
【０１５７】
例えば、上述した実施形態例は本発明を分かりやすく説明するために装置（自動分析装置）の構成を詳細且つ具体的に説明したものであり、必ずしも説明した全ての構成を備えるものに限定されるものではない。
【０１５８】
また、上記の各構成、機能、処理部等は、それらの一部又は全部を、例えば集積回路で設計する等によりハードウェアで実現してもよい。また、上記の各構成、機能等は、プロセッサがそれぞれの機能を実現するプログラムを解釈し、実行することによりソフトウェアで実現してもよい。各機能を実現するプログラム、テーブル、ファイル等の情報は、メモリやハードディスク、ＳＳＤ（Solid State Drive）等の記録装置、又はＩＣカード、ＳＤカード、ＤＶＤ等の記録媒体に置くことができる。
【０１５９】
また、図１に実線で示す制御線や情報線等は、説明上必要と考えられるものを示しており、製品上必ずしも全ての制御線や情報線を示しているとは限らない。実際には殆ど全ての構成が相互に接続されていると考えてもよい。
【０１６０】
また、本明細書において、時系列的な処理を記述する処理ステップは、記載された順序に沿って時系列的に行われる処理はもちろん、必ずしも時系列的に処理されなくとも、並列的あるいは個別に実行される処理（例えば、並列処理あるいはオブジェクトによる処理）をも含むものである。
【符号の説明】
【０１６１】
１…分注ユニット、  １ａ…分注プローブ、  １ｂ…アーム、  １ｃ…駆動機構、  ２…容器、  ３…試料、  ４…ＣＲ発振回路、  ５…第１処理部、  ６…記憶部、  ７…第２処理部、  ７１…演算部、  ７２Ａ…泡接触判定処理部、７２Ｂ…第１液面乖離判定処理部、７２Ｃ…第２液面乖離判定処理部、７２Ｄ…第３液面乖離判定処理部、７２Ｅ…第４液面乖離判定処理部、  ８…第２制御部、  ９…表示部、  １０…自動分析装置、  ２０…コンピューター、  ２１…ＣＰＵ、  ４１…液面乖離及び乖離要因判定テーブル、  ４２…推奨対処手順テーブル
【図１】
2020112380.tif 000003
【図２】
2020112380.tif 000004
【図３】
2020112380.tif 000005
【図４】
2020112380.tif 000006
【図５】
2020112380.tif 000007
【図６】
2020112380.tif 000008
【図７】
2020112380.tif 000009
【図８】
2020112380.tif 000010
【図９】
2020112380.tif 000011
【図１０】
2020112380.tif 000012
【図１１】
2020112380.tif 000013
【図１２】
2020112380.tif 000014
【図１３】
2020112380.tif 000015
【図１４】
2020112380.tif 000016
【図１５】
2020112380.tif 000017
【図１６】
2020112380.tif 000018
【図１７】
2020112380.tif 000019
【図１８】
2020112380.tif 000020
【図１９】
2020112380.tif 000021
【図２０】
2020112380.tif 000022
【図２１】
2020112380.tif 000023