(57)【要約】【課題】冷凍機の劣化兆候の有無を簡便に検知することができる冷凍機劣化兆候検知装置を提供する。【解決手段】前記冷温水の出口温度の設定温度を記憶する設定温度記憶部１５ａと、前記冷温水の出口温度を時系列で検出する検出部１４と、検出した前記出口温度を時系列データとして記憶する出口温度記憶部１５ｂと、前記時系列データとして記憶された前記出口温度のうち、前記冷凍機の運転開始から、前記出口温度が前記設定温度に到達するまでの期間における前記出口温度と前記設定温度との乖離量に基づいて前記冷凍機の応答度に関する値を算出する算出部１６ａと、前記応答度に関する値に基づいて前記劣化兆候を判断する判断部１６ｂと、を備える冷凍機劣化兆候検知装置である。前記時系列で検出する時間間隔は、１５分以上とすることができる。【選択図】図２
2014066452.tif 000002
【特許請求の範囲】
【請求項１】
冷温水を供給する冷凍機（５０）の劣化兆候を検知する冷凍機劣化兆候検知装置（１０）であって、
前記冷温水の出口温度の設定温度を記憶する設定温度記憶部（１５ａ）と、
前記冷温水の出口温度を時系列で検出する検出部（１４）と、
検出した前記出口温度を時系列データとして記憶する出口温度記憶部（１５ｂ）と、
前記時系列データとして記憶された前記出口温度のうち、前記冷凍機の運転開始から、前記出口温度が前記設定温度に到達するまでの期間における前記出口温度と前記設定温度との乖離量に基づいて、前記冷凍機の応答度に関する値を算出する算出部（１６ａ）と、
前記応答度に関する値に基づいて前記劣化兆候を判断する判断部（１６ｂ）と、
を備える、冷凍機劣化兆候検知装置。
【請求項２】
前記算出部が前記乖離量の積分値および前記乖離量の積算値の少なくともいずれか１つを算出する、
請求項１に記載の冷凍機劣化兆候検知装置。
【請求項３】
前記時系列で検出する時間間隔は、１５分以上である、
請求項１または請求項２に記載の冷凍機劣化兆候検知装置。
【請求項４】
前記時系列データは、前記冷凍機の運転開始から前記冷凍機の運転終了までの期間に検出された前記出口温度のデータであって、
前記出口温度記憶部は複数の前記時系列データを記憶し、
前記算出部は前記複数の時系列データに基づき前記応答度に関する値を算出する、
請求項１から３のいずれか１項に記載の冷凍機劣化兆候検知装置。
【請求項５】
前記応答度に関する値は、前記複数の時系列データによって得られる前記応答度に関する分布であり、
前記判断部は、前記応答度に関する分布と第１の基準分布とを、Ｆ検定およびｔ検定のうち少なくともいずれか１つを用いて比較し前記劣化兆候を判断する、
請求項４に記載の冷凍機劣化兆候検知装置。
【請求項６】
前記第１の基準分布が、前記判断部が劣化兆候を有さないと判断した前記冷凍機の前記複数の時系列データによって得られる前記応答度に関する分布である、
請求項５に記載の冷凍機劣化兆候検知装置。
【請求項７】
前記判断部が、前記応答度に関する値および前記冷凍機の運転負荷に関する値と、第２の基準分布とのマハラノビス距離を用いて前記劣化兆候を判断する、
請求項４に記載の冷凍機劣化兆候検知装置。
【請求項８】
前記第２の基準分布が、前記判断部が劣化兆候を有さないと判断した前記冷凍機の前記複数の時系列データによって得られる前記応答度に関する値と、前記複数の時系列データを検出した期間と同一期間における前記冷凍機の運転負荷に関する値と、の分布である、
請求項７に記載の冷凍機劣化兆候検知装置。
【請求項９】
前記冷凍機の前記運転負荷に関する値が、前記冷凍機の運転を開始する時の前記出口温度と前記冷凍機に送り込まれる前記冷温水の温度である入口温度との差である、
請求項７または請求項８に記載の冷凍機劣化兆候検知装置。
【請求項１０】
前記冷凍機の前記運転負荷に関する値が、前記冷凍機の運転を開始する時の前記出口温度と前記入口温度との差と前記冷凍機の前記冷温水の流量との乗積である、
請求項７または請求項８に記載の冷凍機劣化兆候検知装置。
【請求項１１】
前記算出部は、前記時系列データとして記憶された前記出口温度のうち、前記冷凍機の運転開始後、前記冷凍機の前記出口温度が前記設定温度に到達した時から前記冷凍機の運転を停止するまでの前記出口温度と前記設定温度との乖離量の絶対値の積分値に基づいて安定度に関する値をさらに算出し、
前記判断部は、さらに、前記安定度に関する値に基づいて前記劣化兆候を判断する、
請求項１から４のいずれか１項に記載の冷凍機劣化兆候検知装置。
【請求項１２】
前記算出部は、前記時系列データとして記憶された前記出口温度のうち、前記冷凍機の運転中の前記出口温度の平均値と前記設定温度との差に基づいて追従度に関する値をさらに算出し、
前記判断部は、さらに、前記追従度に関する値に基づいて前記劣化兆候を判断する、
請求項１から４または請求項１１のいずれか１項に記載の冷凍機劣化兆候検知装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、冷凍機劣化兆候検知装置に関する。
【背景技術】
【０００２】
従来より、例えば、特許文献１（特開平１１−１９３９７７号公報）に開示されているように、冷凍機が吐出する冷温水の出口温度を監視することにより、冷凍機の劣化兆候の有無を検知する冷凍機劣化兆候検知装置が提案されている。この冷凍機劣化兆候検知装置は、冷凍機の運転開始以降からの出口温度の変化を予測するために予測式を用いる。予測式は、運転開始からの経過時間をパラメータとしている。予測式は、冷凍機の劣化兆候の有無を検知する前に作成する。そして、実際の出口温度と予測式に基づいて予測した出口温度とを比較して、冷凍機の劣化兆候の有無を検知するものである。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
ところで、上記特許文献１に係る冷凍機劣化兆候検知装置は、事前に予測式を作成し、予測式に基づいて冷凍機の出口温度を予測する。予測式は、冷凍機に応じて作成する必要がある。従って、予測式を作成する作業は煩雑である。
【０００４】
本発明は上述した点に鑑みてなされたものであり、本発明の課題は、冷凍機の劣化兆候の有無を簡便に検知することができる冷凍機劣化兆候検知装置を提供することにある。
【課題を解決するための手段】
【０００５】
本発明の第１観点に係る冷凍機劣化兆候検知装置は、冷温水を供給する冷凍機の劣化兆候を検知する冷凍機劣化兆候検知装置であって、設定温度記憶部と、検出部と、出口温度記憶部と、算出部と、判断部と、を備える。設定温度記憶部は、冷温水の出口温度の設定温度を記憶する。検出部は、冷温水の出口温度を時系列で検出する。出口温度記憶部は、検出した出口温度を時系列データとして記憶する。算出部は、時系列データとして記憶された出口温度のうち、冷凍機の運転開始から出口温度が設定温度に到達するまでの期間における出口温度と、設定温度との乖離量に基づいて、冷凍機の応答度に関する値を算出する。判断部は、応答度に関する値に基づいて劣化兆候を判断する。
【０００６】
冷凍機劣化兆候検知装置では、冷温水の出口温度の変化を表すデータが時系列データとして検出される。また、時系列データに基づいて応答度に関する値が算出される。そして、応答度に関する値に基づいて劣化兆候の有無が判断される。よって、冷凍機劣化兆候検知装置は、劣化兆候の有無を簡便に判断することが可能である。
【０００７】
本発明の第２観点に係る冷凍機劣化兆候検知装置は、第１観点に係る冷凍機劣化兆候検知装置であって、算出部は、乖離量の積分値および乖離量の積算値の少なくともいずれか１つを算出する。
【０００８】
冷凍機劣化兆候検知装置では、算出部は、出口温度記憶部が記憶した時系列データに対して、演算処理を行い、乖離量の積分値および乖離量の積算値の少なくともいずれか１つが算出される。よって、迅速に冷凍機の劣化兆候の有無を判断することができる。
【０００９】
本発明の第３観点に係る冷凍機劣化兆候検知装置は、第１観点または第２観点に係る冷凍機劣化兆候検知装置であって、検出部の時系列で検出する時間間隔は、１５分以上である。
【００１０】
冷凍機劣化兆候検知装置では、検出部が検出する時間間隔を、１５分以上にすることができる。よって、少ない計算量で応答度に関する値を求めることができる。
【００１１】
本発明の第４観点に係る冷凍機劣化兆候検知装置は、第１観点から第３観点のいずれかに係る冷凍機劣化兆候検知装置であって、時系列データは、冷凍機の運転開始から冷凍機の運転終了までの期間に検出された出口温度のデータである。出口温度記憶部は複数の時系列データを記憶する。算出部は複数の時系列データに基づき応答度に関する値を算出する。
【００１２】
冷凍機劣化兆候検知装置では、出口温度記憶部は、出口温度の時系列データを記憶する。この時系列データは、冷凍機の運転開始から出口温度が設定温度に到達するまでの複数の期間の時系列データである。その複数の期間の時系列データに基づいて、複数の応答度に関する値を算出する。よって、劣化兆候の有無とトレンドとを正確に把握することができる。
【００１３】
本発明の第５観点に係る冷凍機劣化兆候検知装置は、第４観点に係る冷凍機劣化兆候検知装置であって、応答度に関する値は、複数の時系列データによって得られる応答度に関する分布である。判断部は、応答度に関する分布と第１の基準分布とを、Ｆ検定およびｔ検定のうち少なくともいずれか１つを用いて比較し劣化兆候を判断する。
【００１４】
冷凍機劣化兆候検知装置では、算出部は、複数の出口温度の時系列データに基づいて算出された、応答度に関する分布を求める。さらに、判断部は、劣化兆候の有無を判断できる第１の基準分布と、応答度に関する分布とを比較して劣化兆候を判断する。よって、複数のデータによる分布に基づいて劣化兆候の有無を判断するので、劣化兆候のトレンドを把握できる。
【００１５】
本発明の第６観点に係る冷凍機劣化兆候検知装置は、第５観点に係る冷凍機劣化兆候検知装置であって、第１の基準分布は、応答度に関する分布である。応答度に関する分布は、判断部が劣化兆候を有さないと判断した冷凍機の複数の時系列データによって得られる。
【００１６】
冷凍機劣化兆候検知装置では、第１の基準分布を判断部が劣化兆候を有さないと判断した冷凍機の複数の時系列データによって得られる応答度に関する分布とする。判断部は、劣化兆候を有さないと判断された分布に基づいて、劣化兆候に有無を確認する。これにより信頼性の高い判断結果を得ることができる。
【００１７】
本発明の第７観点に係る冷凍機劣化兆候検知装置は、第４観点に係る冷凍機劣化兆候検知装置であって、判断部が、マハラノビス距離を用いて劣化兆候を判断する。マハラノビス距離は、第２の基準分布の相関に基づき、第２の基準分布と応答度に関する値および冷凍機の運転負荷に関する値との統計上の距離である。
【００１８】
冷凍機劣化兆候検知装置では、算出部が応答度に関する値および冷凍機の運転負荷に関する値と、第２の基準分布とのマハラノビス距離を算出する。判断部は算出されたマハラノビス距離から劣化兆候の有無を判断する。よって、劣化兆候の有無を定量的なデータに基づいて判断することができる。
【００１９】
本発明の第８観点に係る冷凍機劣化兆候検知装置は、第７観点に係る冷凍機劣化兆候検知装置であって、第２の基準分布が、判断部が劣化兆候を有さないと判断した冷凍機の複数の時系列データによって得られる応答度に関する値と、複数の時系列データを検出した期間と同一期間における冷凍機の運転負荷に関する値と、の分布である。
【００２０】
冷凍機劣化兆候検知装置では、算出部は、判断部が劣化兆候を有さないと判断した複数の応答度に関する値と、その応答度のデータを検出した時の運転負荷に関する値と、の分布を第２の基準分布をとする。そして、算出部は、劣化兆候の判断対象となる応答度に関する値および冷却機の運転負荷に関する値のデータと第２の基準分布とのマハラノビス距離を求める。判断部は、算出されたマハラノビス距離から劣化兆候の有無を判断する。よって、判断部は、劣化兆候が無いデータを基準とするマハラノビス距離により、劣化兆候の有無を定量的なデータに基づいて、より確実に判断することができる。
【００２１】
本発明の第９観点に係る冷凍機劣化兆候検知装置は、第７観点または第８観点に係る冷凍機劣化兆候検知装置であって、冷凍機の運転負荷に関する値が、冷凍機の運転を開始する時の出口温度と冷凍機に送り込まれる冷温水の温度である入口温度との差である。
【００２２】
冷凍機劣化兆候検知装置では、冷凍機の運転負荷に関する値は、応答度に関する値を算出した時系列データ検出時の冷凍機の運転を開始する時の出口温度と入口温度との差となる。この運転開始時の出口温度と入口温度の差は、迅速に算出できる。よって、算出部は、マハラノビス距離を容易に算出することができる。判断部は、算出部が算出したマハラノビス距離より迅速、且つ定量的に劣化兆候の有無を判断できる。
【００２３】
本発明の第１０観点に係る冷凍機劣化兆候検知装置は、第７観点または第８観点に係る冷凍機劣化兆候検知装置であって、冷凍機の運転負荷に関する値が、冷凍機の運転を開始する時の出口温度と入口温度との差と、冷凍機の冷温水の流量との乗積である。
【００２４】
冷凍機劣化兆候検知装置では、マハラノビス距離を求めるデータのパラメータである冷凍機の運転負荷に関する値は、応答度に関する値を算出した時系列データ検出時の冷凍機の運転を開始する時の出口温度と入口温度との差と、冷温水の流量との乗積となる。この運転開始時の出口温度と入口温度の差と、冷温水の流量の乗積は、冷凍機に係る運転負荷をより正確に表す。よって、算出部は、劣化兆候を判断するのに適切なマハラノビス距離を算出できる。判断部は、算出部が算出したマハラノビス距離より正確、且つ定量的に劣化兆候の有無を判断できる。
【００２５】
本発明の第１１観点に係る冷凍機劣化兆候検知装置は、第１観点から第４観点のいずれか係る冷凍機劣化兆候検知装置であって、算出部は、時系列データとして記憶された出口温度のうち、冷凍機の運転開始後、冷凍機の出口温度が設定温度に到達した時から冷凍機の運転を停止するまでの出口温度と設定温度との乖離量の絶対値の積分値に基づいて安定度に関する値をさらに算出する。判断部は、さらに、安定度に関する値に基づいて劣化兆候を判断する。
【００２６】
冷凍機劣化兆候検知装置では、冷温水の出口温度の変化を表すデータが時系列データとして検出され、この時系列データを用いて算出した安定度に関する値に基づいて、劣化兆候を判断する。よって、容易なデータ処理で、安定度に関する値に基づいて、確実に冷凍機の劣化兆候を検知できる。
【００２７】
本発明の第１２観点に係る冷凍機劣化兆候検知装置は、第１観点から第４観点、または第１１観点のいずれかに係る冷凍機劣化兆候検知装置であって、算出部は、時系列データとして記憶された出口温度のうち、冷凍機の運転中の出口温度の平均値と設定温度との差に基づいて追従度に関する値をさらに算出する。判断部は、さらに、追従度に関する値に基づいて前記劣化兆候を判断する。
【００２８】
冷凍機劣化兆候検知装置では、時系列データのうち冷凍機の運転中のデータの平均を用いて追従度に関する値を算出する。判断部は、追従度に関する値より劣化兆候を判断する。よって、容易なデータ処理で、追従度に関する値に基づいて、確実に冷凍機の劣化兆候の有無を検知できる。
【発明の効果】
【００２９】
本発明の第１観点に係る冷凍機劣化兆候検知装置では、劣化兆候の有無を簡便に判断することが可能である。
【００３０】
本発明の第２観点に係る冷凍機劣化兆候検知装置では、迅速に冷凍機の劣化兆候の有無を判断することができる。
【００３１】
本発明の第３観点に係る冷凍機劣化兆候検知装置では、少ない計算量で応答度に関する値を求めることができる。
【００３２】
本発明の第４観点に係る冷凍機劣化兆候検知装置では、劣化兆候の有無とトレンドとを正確に把握することができる。
【００３３】
本発明の第５観点に係る冷凍機劣化兆候検知装置では、複数のデータによる分布に基づいて劣化兆候の有無を判断するので、劣化兆候のトレンドを把握できる。
【００３４】
本発明の第６観点に係る冷凍機劣化兆候検知装置では、信頼性の高い判断結果を得ることができる。
【００３５】
本発明の第７観点に係る冷凍機劣化兆候検知装置では、劣化兆候の有無を定量的なデータに基づいて判断することができる。
【００３６】
本発明の第８観点に係る冷凍機劣化兆候検知装置では、劣化兆候の有無を定量的なデータに基づいて、より確実に判断することができる。
【００３７】
本発明の第９観点に係る冷凍機劣化兆候検知装置では、迅速、且つ定量的に劣化兆候の有無を判断できる。
【００３８】
本発明の第１０観点に係る冷凍機劣化兆候検知装置では、正確、且つ定量的に劣化兆候の有無を判断できる。
【００３９】
本発明の第１１観点に係る冷凍機劣化兆候検知装置では、容易なデータ処理で、安定度に関する値に基づいて、確実に冷凍機の劣化兆候を検知できる。
【００４０】
本発明の第１２観点に係る冷凍機劣化兆候検知装置では、容易なデータ処理で、追従度に関する値に基づいて、確実に冷凍機の劣化兆候の有無を検知できる。
【図面の簡単な説明】
【００４１】
【図１】本発明の一実施形態に係るエネルギー管理システムの概略構成図である。
【図２】冷凍機劣化兆候検知装置の概略構成図である。
【図３】出口温度記憶部に記憶されるデータである。
【図４】冷凍機劣化兆候検知装置における処理の流れを示すフローチャートである。
【図５】冷凍機の冷温水の出口温度の推移のグラフである。
【図６】（ａ）劣化兆候が観察される冷温水の設定温度と出口温度との乖離量の積分値のグラフである。（ｂ）劣化兆候が観察されない冷温水の設定温度と出口温度との乖離量の積分値のグラフである。
【図７】（ａ）変形例Ａに係る劣化兆候が観察される冷温水の設定温度と出口温度との乖離量の積算値のグラフである。（ｂ）変形例Ａに係る劣化兆候が観察されない冷温水の設定温度と出口温度との乖離量の積算値のグラフである。
【図８】（ａ）変形例Ｂに係る本発明の劣化兆候が観察される冷温水の出口温度を一定の時間間隔で検出しない場合の出口温度と設定温度との乖離量の積分値のグラフである。（ｂ）変形例Ｂに係る劣化兆候が確認されない冷温水の出口温度を一定の時間間隔で検出しない場合の出口温度と設定温度との乖離量の積分値のグラフである。
【図９】変形例Ｃに係る複数の時系列データを説明するためのグラフである。
【図１０】変形例Ｄに係る劣化兆候を診断する応答度に関する値の分布と、第１の基準分布とを示すグラフである。
【図１１】変形例Ｅに係る劣化兆候を診断する応答度に関する値の分布と、判断部が劣化兆候がないと判断した応答度に関する値の分布を示すグラフである。
【図１２】変形例Ｆに係る複数の時系列データ番号と、その時系列データ番号に関連する応答度に関する値および運転負荷に関する値をあらわす表である。
【図１３】変形例Ｆに係るマハラビノスの基準空間となる第２の基準分布と、劣化兆候を診断する応答度に関する値と運転負荷に関する値とを示すグラフである。
【図１４】変形例Ｇに係るマハラノビスの基準空間となる判断部が劣化兆候を有さないと判断した冷凍機の応答度に関する値と運転負荷に関する値と、劣化兆候を診断する応答度に関する値と運転負荷とを示すグラフである。
【図１５】変形例Ｈに係るマハラノビスの基準空間となる判断部が劣化兆候を有さないと判断した冷凍機の応答度に関する値と、冷温水の入口温度と出口温度との差と、劣化兆候を診断する応答度に関する値と、冷温水の入口温度と出口温度との差を示すグラフである。
【図１６】変形例Ｉに係るマハラノビスの基準空間となる判断部が劣化兆候を有さないと判断した冷凍機の応答度に関する値と、冷温水の入口温度と出口温度との差と流量の乗積と、劣化兆候を診断する応答度に関する値と、冷温水の入口温度と出口温度との差と流量の乗積とを示すグラフである。
【図１７】変形例Ｊに係る安定度に関する値を示すグラフである。
【図１８】変形例Ｊに係るマハラノビスの基準空間となる判断部が劣化兆候を有さないと判断した冷凍機の安定度に関する値と、冷温水の入口温度と出口温度との差と流量の乗積と、劣化兆候を診断する安定度に関する値と、冷温水の入口温度と出口温度との差と流量の乗積とを示すグラフである。
【図１９】変形例Ｋに係る追従度に関する値を示すグラフである。
【図２０】変形例Ｋに係るマハラノビスの基準空間となる判断部が劣化兆候を有さないと判断した冷凍機の追従度に関する値と、冷温水の入口温度と出口温度との差と流量の乗積と、劣化兆候を診断する追従度に関する値と、冷温水の入口温度と出口温度との差と流量の乗積とを示すグラフである。
【発明を実施するための形態】
【００４２】
以下、図１を参照しながら、本発明にかかる冷凍機劣化兆候検知装置１０を含むエネルギー管理装置２００について説明する。
【００４３】
（１）エネルギー管理システム
図１にエネルギー管理システム２００を示す。エネルギー管理システム２００は、建物内で消費されるエネルギーを管理するシステムである。エネルギー管理システム２００は、熱源システム１００と冷凍機劣化兆候検知装置１０とを含む。エネルギー管理システム２００は、冷凍機劣化兆候検知装置１０により、熱源システム１００に含まれる冷凍機５０の劣化兆候の有無を検知する。以下、熱源システム１００および冷凍機劣化兆候検知装置１０の構成について説明する。
【００４４】
（１−１）熱源システム
熱源システム１００は、温度調整された空気を空調対象空間に供給するシステムである。熱源システム１００は、冷凍機５０と、１次ポンプ４１と、温度センサ６０と、上流側配管７１と、空調機３２と、下流側配管７２と、を有する。
【００４５】
冷凍機５０は、熱源システム１００内の冷温水を冷却する。１次ポンプ４１は、冷温水を冷凍機５０に送る。温度センサ６０は、冷凍機５０が冷温水を吐出する部分に設けられる。温度センサ６０は、冷凍機５０が吐出する冷温水の温度（出口温度）を測定する。温度センサ６０は、例えば、信号線８０などのネットワークにより冷凍機劣化兆候検知装置１０と接続されている。温度センサ６０は、冷温水の出口温度を測定すると、出口温度に関する信号を冷凍機劣化兆候検知装置１０に送る。上流側配管７１は、冷凍機５０と空調機３２とを繋ぐ配管である。上流側配管７１は、冷凍機５０から吐出された冷温水を空調機３２に通す。空調機３２は、空調対象空間における空気の温度調整を行う。下流側配管７２もまた、空調機３２と冷凍機５０とを繋ぐ配管である。下流側配管７２は、空調機３２によって室内空気と熱交換された冷温水を１次ポンプ４１に通す。
【００４６】
（１−２）冷凍機劣化兆候検知装置１０の構成
冷凍機劣化兆候検知装置１０は、熱源システム１００に設けられた冷凍機５０の劣化兆候を検知する装置である。冷凍機劣化兆候検知装置１０は、上述した温度センサ６０から送られる出口温度に関する信号を受信し、冷凍機５０の劣化兆候の有無を検知する。
【００４７】
以下、図２を参照して、冷凍機劣化兆候検知装置１０を説明する。図２に冷凍機劣化兆候検知装置１０の概略構成図を示す。冷凍機劣化兆候検知装置１０は、熱源システム１００にある、温度センサ６０の出力値を検知する。冷凍機劣化兆候検知装置１０は、入力部１１、出力部１２、時間管理部１３、検出部１４、記憶部１５、および制御部１６を備える。
【００４８】
（１−２−１）入力部
入力部１１は、主として操作ボタン、キーボード、およびマウス等から構成される。入力部１１は、ユーザーから冷凍機５０の出口温度の設定温度を受け付ける。入力部１１が受け付けた設定温度は、後述の設定温度記憶部１５ａに記憶される。
【００４９】
また、入力部１１は、ユーザーからの出口温度の検出の開始要求を受け付ける。入力部１１は、ユーザーからの出口温度の検出要求の終了を受け付ける。
【００５０】
（１−２−２）出力部
出力部１２は、主としてディスプレイから構成される。出力部１２には、後述する記憶部１５に記憶される各種情報を示す管理画面が表示される。また、出力部１２には、後述する算出部１６ａが算出した値、後述する判断部１６ｂが判断した結果等が表示される。
【００５１】
（１−２−３）時間管理部
時間管理部１３は、冷凍機劣化兆候検知装置１０の実行する各種制御の時間管理を行う。例えば、後述する検出部１４が出口温度を検出する間隔を決めるタイマーとなる。
【００５２】
（１−２−４）検出部
検出部１４は、温度センサ６０の出力値を検出する。検出部１４は、前述の時間管理部１３の管理により、温度センサ６０の出力値の検出を１５分以上のサンプリンプ間隔とすることができる。また、検出部１４は、図９に示されるように冷凍機５０の複数回の運転の開始と終了を含む出口温度の変化を検出する。
【００５３】
（１−２−５）記憶部
記憶部１５は、ハードディスク等から構成され、設定温度記憶部１５ａと、出口温度記憶部１５ｂと、を有する。
【００５４】
（１−２−５−１）設定温度記憶部
設定温度記憶部１５ａは、入力部１１がユーザーより受け付けた冷凍機５０の出口温度の設定温度を記憶する。
【００５５】
（１−２−５−２）出口温度記憶部
出口温度記憶部１５ｂは、温度センサ６０が検出した冷凍機５０の出口温度として記憶する。出口温度記憶部１５ｂは、出口温度の記憶を、入力部１１が開始要求を受け付けてから、終了要求を受け付けるまで行う。出口温度は、時系列データとして出口温度記憶部１５ｂに記憶される。出口温度記憶部１５ｂが記憶する時系列データは、例えば図３に示されるように、出口温度を検出した年月日および時刻と、出口温度と、冷凍機５０の運転または停止の状況と、が記憶される。
【００５６】
（１−２−６）制御部
制御部１６は、ＣＰＵ，ＲＯＭ，ＲＡＭ等から構成される。制御部１６は、上述の記憶部１５に記憶されているプログラムを読み出して実行する。制御部１６は、プログラムを実行することにより設定温度記憶部１５ａに記憶された設定温度を使用し、出口温度記憶部１５ｂに記憶された時系列データを解析する。制御部１６は、時系列データの解析結果より冷凍機５０の劣化兆候の有無を判断する。制御部１６は、上述の解析と判断をするために、主として算出部１６ａ、判断部１６ｂとして機能する。
【００５７】
（１−２−６−１）算出部
算出部１６ａは、出口温度記憶部１５ｂに記憶された、図３に示される時系列データから応答度に関する値を算出する。ここで応答度とは、冷凍機の運転を開始してから、設定温度に到達するまでの冷凍機の応答の程度を意味する。算出部１６ａは、応答度に関する値として、乖離量の積分値を算出する。ここで、乖離量とは、検出された出口温度から設定温度記憶部１５ａに記憶されている設定温度を減算した値をいう。以下、図６（ａ）を用いて乖離量の積分値について説明する。
【００５８】
図６に示す例では、記憶されている出口温度の設定温度は７度である。また、冷凍機５０の運転開始の出口温度は２６度である。出口温度は、冷凍機５０の運転開始から７分後に設定温度に到達する。冷凍機５０の運転開始から、設定温度に到達するまでの期間における出口温度と設定温度の乖離量を積分すると、図６（ａ）に示されるＳ１の部分が算出される。本実施形態では、Ｓ１で示す面積を乖離量の積分値という。
【００５９】
（１−２−６−２）判断部
判断部１６ｂは、算出部１６ａが算出した乖離量の積分値に基づいて、冷凍機５０が劣化兆候を有するか否かを判断する。以下、図６（ａ）と図６（ｂ）に示すグラフを参照しながら冷凍機５０が劣化兆候を有するか否かの判断方法について説明する。図６（ａ）は劣化兆候が観察される冷温水の設定温度と出口温度との乖離量の積分値のグラフである。図６（ｂ）は、劣化兆候が観察されない冷温水の設定温度と出口温度との乖離量の積分値のグラフである。図６（ｂ）の乖離量の積分値は、判断部１６ｂが、図６（ａ）の劣化兆候の有無を判断する場合の基準となる。冷凍機５０に劣化兆候があると応答度は悪化する。具体的には、冷凍機５０に劣化兆候があると、冷凍機５０の運転を開始してから、設定温度に到達するまでの時間が長くなる。図6（ａ）では、冷凍機５０の運転開始から７分で出口温度が設定温度に到達している。図６（ｂ）では、冷凍機５０の運転開始から６分で出口温度が設定温度に到達している。よって、図６（ｂ）の乖離量の積分値を示す面積Ｓ２は、図６（ａ）のＳ１よりも小さくなっている。よって、図６（ａ）の方が冷凍機５０の応答度が遅く、図６（ｂ）は、劣化兆候があると判断部１６ｂは、判断できる。
【００６０】
具体的には、判断部１６ｂは、劣化について問題が無い冷凍機５０の時系列データに基づいて、劣化兆候の有無を判断する基準となる乖離量の積分値（基準積分値）を、記憶部１５に記憶させる。次に、算出部１６ａが算出した乖離量の積分値と基準積分値とを比較し、基準積分値より算出した乖離量の積分値が大きい、または基準積分値と算出した積分値の差が所定の値を超えた場合に、乖離量の積分値を求めた冷凍機５０を劣化兆候があると判断する。
【００６１】
（２）冷凍機劣化兆候検知装置で行われる制御処理について
以下、冷凍機劣化兆候検知装置１０における、冷凍機劣化兆候の検知の流れを図４を用いて説明する。
【００６２】
まず、ステップＳ１０１では、入力部１１がユーザーから冷凍機５０の出口温度の設定温度を受付け、設定温度記憶部１５ａに記憶させる。
【００６３】
ステップＳ１０２では、入力部１１がユーザーからの出口温度の検出の開始を受け付けたか否かを判定する。受け付けたと判定される場合は、ステップS１０３に移行する。受け付けていないと判定される場合は、ステップS１０２に移行する。
【００６４】
ステップＳ１０３では、検出部１４が、温度センサ６０が測定した冷凍機５０の出口温度を検出する。
【００６５】
ステップＳ１０４では、検出部１４が、ステップS１０３で検出した冷凍機５０の出口温度を出口温度記憶部１５ｂに記憶する。具体的には、図３に示すように、出口温度と、出口温度を検出した年月日および時刻と、冷凍機５０の運転または停止の状況とを関連付けて記憶する。
【００６６】
なお、ステップＳ１０２において入力部１１が出口温度の検出開始要求を受け付け、ステップＳ１０６において入力部１１が検出終了を受け付ける間に、冷凍機５０が複数回の運転の開始と停止を繰り返す。そうすると、出口温度記憶部１５ｂに記憶された出口温度は、例えば、図５のグラフに示されるように推移する。
【００６７】
ステップＳ１０５では、時間管理部１３が出口温度を検出するための所定の待ち時間を管理する。所定の待ち時間が終了するとステップS１０６に移行する。なお、この所定の待ち時間は、冷凍機５０の出口温度のサンプリング間隔を決める待ち時間となる。この所定のサンプリングの待ち時間は、任意の時間にすることができるが、本実施形態では、１５分以上とすることが好ましい。
【００６８】
ステップＳ１０６では、入力部１１がユーザーからの出口温度の検出の終了を受け付けたか否かを判定する。受け付けたと判定される場合は、ステップＳ１０７へ移行する。受け付けていないと判定される場合は、ステップＳ１０３に移行する。
【００６９】
ステップＳ１０７では、算出部１６ａが応答度に関する値として、乖離量の積分値を算出する。
【００７０】
具体的には、まず出口温度記憶部１５ｂに記憶されている時系列データに基づいて、算出部１６ａが冷凍機５０の運転を開始した時刻を把握する（図３参照）。また、算出部１６ａは、設定温度記憶部１５ａに記憶されている設定温度を参照する。冷凍機５０の運転が開始すると、時間の経過と共に出口温度は、設定温度に近づき、ある時刻で出口温度は、設定温度に到達する。この冷凍機５０の運転開始時刻から、出口温度が設定温度となった時刻までの間において、出口温度から設定温度を減算した値の積分した値を算出する。これは、図６（ａ）のＳ１の部分となる。
【００７１】
ステップＳ１０８では、判断部１６ｂが劣化兆候を判断する。判断部１６ｂは、図６（ａ）のＳ１を求めた後、Ｓ２と比較し、Ｓ１が大きい場合、またはＳ１とＳ２の差が所定の値を超えている場合に冷凍機５０に劣化兆候があると判断する。
【００７２】
ステップＳ１０９では、出力部１２が、判断部１６ｂが判断した結果をディスプレイ等に表示する。
【００７３】
（３）特徴
（３−１）
本実施形態では、冷凍機５０の出口温度を予測する予測式を事前に用意する必要はない。よって、簡易に且つ効率的に冷凍機５０の劣化兆候の有無を判断することができる。
【００７４】
（３−２）
本実施形態では、冷温水の出口温度の変化を表すデータを時系列データとして検出し、記憶する。次に、算出部１６ａは、記憶した時系列データを解析することにより応答度に関する値を算出し、劣化兆候を判断する。よって、冷凍機５０の内部を透視する目的のＸ線照射装置などの他の装置を用意する必要がない。また、冷温水の出口温度を観察する期間が冷凍機５０の運転開始から、出口温度が設定温度に到達するまでの期間なので、処理するデータ量が少ない。さらに、算出部が行うデータ処理は、設定温度と出口温度の乖離量に関する演算処理だけである。よって、冷凍機５０の劣化兆候の有無を、他の装置を使うことなく、少ないデータ量で、容易な演算処理により簡便、確実に検知することができる。
【００７５】
（３−３）
本実施形態では、算出部１６ａは、出口温度記憶部１５ｂが記憶した時系列データに対して、容易な演算処理を行うだけなので、簡便に冷凍機５０の劣化兆候の有無を判断することができる。
【００７６】
（３−４）
本実施形態では、検出部１４が時系列で検出する出口温度のサンプリング間隔を、１５分以上に設定している。従って、検出する時系列のデータの数が、より少なくなる。算出部１６ａの取り扱うデータの量が少なくなるので、より迅速に冷凍機５０の劣化兆候の有無を判断することができる。
【００７７】
（３−５）
本実施形態では、冷凍機５０の音または振動を高いサンプリングレートで採取し、採取した音または振動の周波数成分を解析する場合のように多くのサンプリングデータを取り扱わない。よって、制御部１６は、データ解析に高い計算能力を必要としない。
【００７８】
（４）変形例
以上、本発明の実施形態について図面に基づいて説明したが、具体的な構成は、上記の実施形態に限られるものでなく、発明の要旨を逸脱しない範囲で変更可能である。
【００７９】
以下に本実施例の変形例を示す。なお、複数の変形例は適宜組み合わされてもよい。
【００８０】
（４−１）変形例A
上記実施形態では、応答度に関する値を、冷凍機５０の運転開始時刻から出口温度が設定温度に到達したまでの間（到達必要時間）における、出口温度と設定温度との乖離量を積分した値としていた（図６（ａ）参照）。すなわち、上記実施形態では、到達必要時間における出口温度と設定温度との乖離量を面積として求めた。
【００８１】
ここで、応答度に関する値は、図７（ａ）に示されるように、到達必要時間における、出口温度と設定温度との乖離量を積算した値（乖離量の積算値）であってもよい。言い換えると、乖離量の積算値とは、到達必要時間において、検出部１４によって所定の時間間隔で検出された、それぞれの出口温度と、設定温度との乖離量（符号ａ０からａ８）を積算した値である。以下に、応答度に関する値を、乖離量の積算値とした場合の劣化兆候判断処理について、図７（ａ）と図７（ｂ）とを用いて説明する。
【００８２】
図７（ａ）は、劣化兆候を判断する対象となる冷凍機５０の応答度に関する値を示す図である。一方、図７（ｂ）は、劣化兆候がないと判断された冷凍機５０の応答度に関する値を示す図である。図７（ａ），（ｂ）の符号ａ０〜ａ８，ｂ０〜ｂ７は、それぞれ、到達必要時間に検出部１４によって検出された出口温度と、設定温度との乖離量を示す。劣化兆候がないと判断された冷凍機５０の応答度に関する値は、例えば、新しい冷凍機５０を建物に設置した直後に応答度に関する値（正常積算値）を求めることで可能である。正常積算値は、入力部１１により、あらかじめ記憶部１５に記憶させておく。
【００８３】
判断部１６ｂは、劣化兆候の判断対象となる冷凍機５０について乖離量の積算値（対象積算値）を求めた後、対象積算値と正常積算値とを比較する。判断部１６ｂは、対象積算値が正常積算値より大きい場合、または対象積算値が正常積算値を基準とした正常領域値を外れた場合に冷凍機５０に劣化兆候があると判断する。
【００８４】
乖離量の積算値の算出は、乖離量の積分値の算出よりもさらに演算が容易である。したがって、容易且つ迅速に劣化兆候の有無を判断することができる。
【００８５】
（４−２）変形例Ｂ
上記実施形態では、検出部１４が行う出口温度の検出は一定の時間間隔で行われている。
【００８６】
ここで、図８（ａ）の横軸である時間軸上の出口温度のタイミングが示すように、出口温度の検出は一定の時間間隔でなくても良い。ここで、出口温度の検出を一定の時間間隔で検出しないことを不定期間隔検出いう。以下に乖離量の積算値を不定期間隔検出した場合の劣化兆候判断処理について、図８（ａ）と図８（ｂ）を用いて説明する。
【００８７】
図８（ａ）は、劣化兆候を判断する対象となる冷凍機５０の応答度に関する値を示す図である。一方、図８（ｂ）は、劣化兆候がないと判断された冷凍機５０の応答度に関する値を示す図である。図８（ａ），（ｂ）のＲ１，Ｒ２は、それぞれ、到達必要時間の乖離量の積分値である。
【００８８】
ここで、図８（ａ），（ｂ）のＲ１，Ｒ２は、不定期間隔検出による乖離量の積分値であるが、冷凍機５０の運転開始から、同じタイミングで、出口温度を検出している。よって、Ｒ１とＲ２を比較して、冷凍機５０の劣化兆候の有無を検知することができる。
【００８９】
不定期間隔検出において、劣化兆候がないと判断された冷凍機５０の乖離量の積算値は、例えば、新しい冷凍機５０を建物に設置した直後の乖離量の積算値（不定期間隔検出正常値）を求めることで可能である。不定期間隔検出正常値は、入力部１１により、あらかじめ記憶部１５に記憶させておく。
【００９０】
判断部１６ｂは、劣化兆候の判断対象となる冷凍機５０について不定期間隔検出で算出した乖離量の積算値（不定期間隔検出対象値）を求めた後、不定期間隔検出正常値と不定期間隔検出対象値を比較する。判断部１６ｂは、不定期間隔検出対象値が不定期間隔検出正常値より大きい場合、または不定期間隔検出対象値が不定期間隔検出正常値を基準とした正常領域を外れた場合に冷凍機５０に劣化兆候があると判断する。
【００９１】
冷凍機５０の出口温度の変化については、冷凍機５０の運転開始時など、変化が急な時間がある。そのような時間において、出口温度の検出の時間間隔を短時間する不定期間隔検出を行うことにより、乖離量の積分値の精度を上げることなどが可能である。
【００９２】
（４−３）変形例Ｃ
上記実施形態では、算出部１６ａは、冷凍機５０の到達必要時間の時系列データに基づいて、１つの応答度に関する値を算出した。
【００９３】
ここで、変形例Ｃでは、複数の到達必要時間の時系列データに基づき複数の応答度に関する値を算出する。
【００９４】
図９は、冷凍機運転中ａから冷凍機運転中ｅまで、冷凍機５０が複数回運転されたことを示している。ここで、検出部１４が図９に示されるように、複数回の到達必要時間の出口温度を検出してもよい。そして、出口温度記憶部１５ｂが、上述の期間の出口温度のデータを出口温度記憶部１５ｂに記憶してもよい。次に、算出部１６ａが、それぞれの到達必要時間の出口温度の時系列データに基づいて、複数の応答度に関する値を算出してもよい。さらに、判断部１６ｂが複数の応答度に関する値に基づいて、冷凍機５０の劣化兆候の有無を判断してもよい。
【００９５】
ここで、上述の冷凍機５０の劣化兆候の検知方法とは別の方法であるが、算出部１６ａが複数の応答度に関する値を算出し、その複数の応答度に関する値の平均値を算出しても良い。判断部１６ｂが、この平均値に基づき冷凍機５０の劣化兆候の有無を判断してもよい。
【００９６】
上述の２つの冷凍機５０の劣化兆候の有無の判断方法は、複数の応答度に関する値より、冷凍機５０の劣化兆候の有無を判断するので、劣化兆候検知の信頼性をあげることができる。
【００９７】
（４−４）変形例Ｄ
上記変形例Ｃでは、複数の到達必要時間の時系列データから算出した、複数の応答度に関する値に基づいて、冷凍機５０の劣化兆候の有無を判断部１６ｂが判断した。
【００９８】
ここで、変形例Ｄでは、複数の到達必要時間の時系列データから算出した、複数の応答度に関する値の分布のデータとし、その分布のデータに基づいて、冷凍機５０の劣化兆候の有無を判断してもよい。
【００９９】
以下に応答度に関する値の分布のデータに基づいた劣化兆候判断処理について図１０（ａ）と図１０（ｂ）とを用いて説明する。
【０１００】
ここで、図１０（ｂ）は、劣化兆候の判断対象となる冷凍機５０についての複数の応答度に関する値の分布データである診断分布とする。図１０（ａ）は、第１の基準分布とする。第１の基準分布としては、特に限定されないが、例えば冷凍機５０を設置して始動した時の応答度に関する値の分布、冷凍機５０の過去の定期点検の応答度に関する値の分布、他の新しい複数の冷凍機５０から採取した応答度に関する値の分布などがある。
【０１０１】
判断部１６ｂは、図１０（ｂ）に示される診断分布と、図１０（ａ）に示される第１の基準分布とにＦ検定を行い、等分散であるかを判断する。次に等分散であると判断された場合は、診断分布と第１の基準分布とにｔ検定を行う。Ｆ検定の結果等分散でないと判断された場合は、診断分布と第１の基準分布とにＷｅｌｃｈのｔ検定を行う。この結果、判断部１６ｂが診断分布と第１の基準分布とが統計的に有意差がなく、同じ分布と判断した場合は、冷凍機５０の劣化兆候が無いと判断することができる。判断部１６ｂが統計的に有意差があり、異なる分布と判断した場合は、冷凍機５０の劣化兆候が有ると判断することができる。
【０１０２】
これにより、判断部１６ｂは、複数の応答度の分布同士で、冷凍機５０の劣化兆候を判断するので、劣化兆候の有無と、劣化しているかのトレンドを確実に判断することができる。
【０１０３】
（４−５）変形例Ｅ
上記変形例Ｄでは、診断分布と第１の基準分布に基づいて、冷凍機の劣化兆候の有無を判断した。上記変形例Ｄでは、判断部１６ｂが第１の基準分布を、特に限定していない。
【０１０４】
ここで、変形例Ｅでは、第１の基準分布は、判断部１６ｂが、劣化兆候を有さないと判断した冷凍機５０の複数の時系列データによって得られる応答度に関する値のデータであってもよい。
【０１０５】
以下に応答度に関する値の分布のデータに基づいた劣化兆候判断処理について図１１（ｄ）と図１１（ｃ）とを用いて説明する。
【０１０６】
ここで、図１１（ｄ）は、診断分布である。図１１（ｃ）を判断部１６ｂが劣化兆候がないと判断した冷凍機５０の応答度に関する値の分布（正常分布）である。
【０１０７】
ここで、判断部１６ｂが図１１（ｃ）に示される正常分布と図１１（ｄ）に示される診断分布に対して、変形例Ｄの判断部１６ｂの処理と同じようにＦ検定を行い、等分散であるかを判断する。次に等分散であると判断した場合は、比較する分布と第１の基準分布とにｔ検定を行う。Ｆ検定の結果等分散でないと判断した場合は、比較する分布と第１の基準分布とにＷｅｌｃｈのｔ検定を行う。この結果、判断部１６ｂが比較する分布と第１の基準分布とが統計的に有意差がなく、同じ分布と判断した場合は、冷凍機５０の劣化兆候が無いと判断することができる。判断部１６ｂが統計的に有意差があり、異なる分布と判断した場合は、冷凍機５０の劣化兆候が有ると判断することができる。
【０１０８】
これにより、判断部１６ｂは、判断部１６ｂが劣化兆候がないと判断した分布と、診断分布の分布同士を比較するため、信頼性の高い劣化兆候の有無の判断をすることができる。
【０１０９】
（４−６）変形例Ｆ
上記実施形態では、検出部１４が冷凍機５０の出口温度の時系列データを検出して、算出部１６ａが乖離量の積分値を算出していた。
【０１１０】
ここで変形例Ｆでは、図１２に示すように、検出部１４が、応答度に関する値を求め、冷凍機５０の時系列データを検出する期間における冷凍機の運転負荷に関する値を、応答度に関する値と共にデータとして出口温度記憶部１５ｂに記憶させてもよい。
【０１１１】
運転負荷に関する値を表すデータとしては、特に限定されないが、例えば、圧縮機を持つ冷凍機５０の圧縮機の負荷、冷凍機５０の外気温と出口温度の設定温度との差などがあげられる。
【０１１２】
算出部１６ａは、図１２で出口温度記憶部１５ｂに記憶させた複数の応答度に関する値と運転負荷に関する値のデータを図１３で示すように横軸に応答度に関する値、縦軸に運転負荷に関する値を示す２次元の散布図を作成する。次に、この散布図上にプロットされた複数の応答度に関する値と運転負荷に関する値のデータをマハラノビスの基準空間（第２の基準分布）Ａとする。算出部１６ａは、劣化兆候の有無を判断する冷凍機５０の応答度に関する値と運転負荷に関する値のデータを図１３のａ点とする。そして、算出部１６ａは、ａ点とマハラノビスの基準空間Ａとした第２の基準分布とのマハラノビス距離を算出する。
【０１１３】
判断部１６ｂは、算出部１６ａが算出したマハラノビス距離が記憶部１５に記憶させた所定の値を超えている場合は劣化兆候があると判断することができる。
【０１１４】
これにより、判断部は、マハラノビス距離という客観的な数値により冷凍機５０の劣化兆候の有無を判断することができる。
【０１１５】
（４−７）変形例Ｇ
上記変形例Ｆでは、第２の基準分布となる応答度に関する値と運転負荷に関する値について、特に限定はしていない。
【０１１６】
ここで、変形例Ｇでは、応答度に関する値と運転負荷に関する値を、判断部１６ｂが劣化兆候を有さないと判断した冷凍機５０のそれぞれの値とすることができる。そして、劣化兆候がないと判断された、応答度に関する値と運転負荷に関する値を図１４に示す２次元の散布図上にプロットする。そして、散布図上にプロットされたデータをマハラノビスの基準空間Ｂとすることができる。
【０１１７】
判断部１６ｂが劣化兆候を有さないと判断した冷凍機５０としては、新しい冷凍機５０を設置したときの応答度に関する値と運転負荷に関する値のデータ、他の正常に動作している冷凍機５０の応答度に関する値と運転負荷に関する値のデータなどがあげられる。
【０１１８】
算出部１６ａは、劣化兆候の有無を判断する冷凍機５０の応答度に関する値と運転負荷に関する値のデータを図１４のｂ点とする。そして、算出部１６ａは、図１３のｂ点とマハラノビスの基準空間Ｂとのマハラノビス距離を算出する。
【０１１９】
判断部１６ｂは、算出部１６ａが算出したマハラノビス距離が記憶部１５に記憶させた所定の値を超えている場合は劣化兆候があると判断することができる。
【０１２０】
これにより、判断部１６ｂは、劣化兆候がない応答度に関する値と運転負荷に関する値とで構成されるマハラノビスの基準空間Ｂからのマハラノビス距離という客観的な数値により冷凍機５０の劣化兆候を判断する。よって、より正確に冷凍機５０の劣化兆候の有無を判断することができる。
【０１２１】
（４−８）変形例Ｈ
上記変形例Ｇでは、応答度に関する値と運転負荷に関する値を判断部１６ｂが劣化兆候を有さないと判断した冷凍機５０のそれぞれの値としている。
【０１２２】
ここで変形例Ｈでは、運転負荷に関する値を冷凍機５０の運転開始時に冷凍機５０に送り込まれる冷温水の温度である入口温度と冷凍機５０の出口温度との差（温度差）とすることができる。図１５に示すように、判断部１６ｂが劣化兆候が無いと判断した冷凍機５０の応答度に関する値と温度差とを第２の基準分布とするマハラノビスの基準空間Ｃとした散布図を示す。算出部１６ａは、劣化兆候の有無を判断する冷凍機５０の応答度に関する値と温度差とのデータを図１５のｃ点とする。そして、算出部１６ａは、図１５のｃ点とマハラノビスの基準空間Ｃとのマハラノビス距離を算出する。
【０１２３】
判断部１６ｂは、算出部１６ａが算出したマハラノビス距離が記憶部１５に記憶させた所定の値を超えている場合は劣化兆候があると判断することができる。
【０１２４】
これにより、判断部は、劣化兆候がない応答度に関する値と温度差で構成されるマハラノビスの基準空間Ｃからのマハラノビス距離という冷凍機５０での冷温水の温度の変化による客観的な数値により冷凍機５０の劣化兆候を判断するため、より正確に冷凍機５０の劣化兆候の有無を判断することができる。
【０１２５】
（４−９）変形例Ｉ
上記変形例Ｈでは、運転負荷に関する値を冷凍機５０の運転開始時の温度差としている。
【０１２６】
ここで、変形例Ｉでは、運転負荷に関する値を温度差と冷凍機５０を流れる冷温水との流量との乗積（温度差流量乗積値）とすることができる。図１６に示すように、判断部１６ｂが劣化兆候が無いと判断した冷凍機５０の応答度に関する値と温度差流量乗積値とを第２の基準分布とするマハラノビスの基準空間Ｄとする散布図を示す。算出部１６ａは、図１６の劣化兆候の有無を判断する冷凍機５０の応答度に関する値と温度差流量乗積値とのデータを図１６のｄ点とする。そして、算出部１６ａは図１６のｄ点と、マハラノビスの基準空間Ｄとのマハラノビス距離を算出する。
【０１２７】
判断部１６ｂは、算出部１６ａが算出したマハラノビス距離が記憶部１５に記憶させた所定の値を超えている場合は劣化兆候があると判断することができる。
【０１２８】
これにより、判断部は、劣化兆候がない応答度に関する値と温度差流量乗積値で構成されるマハラノビスの基準空間からのマハラノビス距離という冷凍機５０での冷温水の温度差と冷温水の流量による客観的な数値により冷凍機５０の劣化兆候を判断するため、より正確に冷凍機５０の劣化兆候の有無を判断することができる。
【０１２９】
（４−１０）変形例Ｊ
上記実施形態では、算出部１６ａが冷凍機５０の応答度に関する値である乖離量の積分値を算出する。そして、判断部１６ｂが乖離量の積分値に基づいて、冷凍機５０の劣化兆候の有無を判断する。
【０１３０】
ここで、変形例Ｊでは、算出部１６ａが安定度に関する値を算出し、判断部１６ｂが安定度に関する値に基づいて劣化兆候を判断することもできる。
【０１３１】
ここで、安定度とは冷凍機の運転開始後、出口温度が設定温度となってからの設定温度に対する安定の程度を意味する。
【０１３２】
以下、安定度に関する値について、いくつかの変形例を説明する。
【０１３３】
（４−１０−１）
算出部１６ａは、安定度に関する値として、冷凍機５０の運転開始後、設定温度に到達した時から冷凍機５０を停止するまでの出口温度と設定温度との乖離量の絶対値の積分値（乖離量の絶対値の積分値）を算出する。
【０１３４】
以下、図１７を用いて乖離量の絶対値の積分値を説明する。
【０１３５】
冷凍機５０の運転開始後、出口温度は下がり、設定温度に到達する。この後、冷凍機５０は設定温度を保持する運転となる。実際の出口温度は設定温度を上回ったり、下回ったりする。この出口温度と設定温度との差の絶対値の積分値である、図１７の中で斜線を引いた部分の面積が、乖離量の絶対値の積分値である。ここで、冷凍機５０に劣化兆候があると出口温度は、安定性を失い設定温度に対して、大きく振動する。よって、乖離量の絶対値の積分値は、正常な冷凍機５０の値よりも大きくなる。
【０１３６】
乖離量の絶対値の積分値を用いた、劣化兆候の有無の判断について説明する。判断部１６ｂは、劣化について問題が無い冷凍機５０の時系列データに基づいて、劣化兆候の有無を判断する基準となる乖離量の絶対値の積分値（基準絶対値積分値）を、記憶部１５に記憶させる。次に、算出部１６ａが算出した乖離量の絶対値の積分値と基準絶対値積分値を比較し、基準絶対積分値より算出した乖離量の絶対値の積分値が大きい場合、または、算出した乖離量の絶対値の積分値が基準絶対積分値を基準とした正常領域を外れた場合に、冷凍機５０に劣化兆候があると判断する。
【０１３７】
（４−１０−２）
また別の、算出部１６ａが算出する安定度に関する値として、温度差の標準偏差を算出することができる。温度差の標準偏差とは、冷凍機５０の運転開始後、出口温度が設定温度となった後から、運転を停止までの出口温度の値と設定温度との差を自乗し、その自乗した値の総和の平方根である。
【０１３８】
冷凍機５０に劣化兆候があると、出口温度は、設定温度を、ほぼ中心として大きく振動するため、温度差の標準偏差も大きくなる。ユーザーは、記憶部１５に、劣化について問題が無い冷凍機５０の時系列データに基づく温度差の標準偏差を基準温度差標準偏差として記憶させる。判断部１６ｂは、記憶部１５に記憶された基準温度差標準偏差と、劣化兆候の有無を判断する冷凍機５０の温度差の標準偏差とを比較することにより、冷凍機５０の劣化兆候の有無を判断することができる。
【０１３９】
（４−１０−３）
さらに、別の算出部１６ａが算出する安定度に関する値として、安定度の分布の標準偏差を算出することができる。安定度の分布の標準偏差について説明する。複数回の冷凍機５０の運転により、複数の乖離量の絶対値の積分値を算出することができる。この複数の乖離量の絶対の積分値を、安定度の分布のデータとすることができる。分布のデータはばらつきを有するため、この分布するデータの標準偏差を求めることができる。この標準偏差を安定度の分布の標準偏差とする。
【０１４０】
冷凍機５０に劣化兆候の傾向があると、安定度の分布の標準偏差も大きくなる。ユーザーは、記憶部１５に、劣化について問題が無い冷凍機５０の複数の時系列データに基づく安定度の分布の標準偏差を基準安定度標準偏差として記憶させる。判断部１６ｂは、記憶部１５に記憶された基準安定度標準偏差と、劣化兆候の有無を判断する冷凍機５０の安定度の分布の標準偏差とを比較することにより、冷凍機５０の劣化兆候の有無を判断することができる。
【０１４１】
（４−１０−４）
また、算出部１６ａは、劣化兆候の無い複数の時系列データに基づいて算出した複数の乖離量の絶対値の積分値と、その乖離量の絶対値の積分値を算出した時の冷凍機５０の運転開始時の温度差流量乗積値と、のデータを図１８に示すようにマハラノビスの基準空間Ｅとすることができる。算出部１６ａは、劣化兆候の有無を判断する冷凍機５０の乖離量の絶対値の積分値と温度差流量乗積値のデータを図１８のｅ点とする。そして、算出部１６ａは、図１８のｅ点とマハラノビスの基準空間Ｅとのマハラノビス距離を算出する。
【０１４２】
判断部１６ｂは、算出部１６ａが算出したマハラノビス距離が記憶部１５に記憶させた所定の値よりも大きい場合に、乖離量の積分値を求めた冷凍機５０を劣化兆候があると判断することができる。
【０１４３】
（４−１１）変形例Ｋ
上記実施形態では、算出部１６ａが冷凍機５０の応答度に関する値である乖離量の積分値を算出する。そして、判断部１６ｂが乖離量の積分値に基づいて、冷凍機５０の劣化兆候の有無を判断する。
【０１４４】
ここで、変形例Ｋでは、算出部１６ａが追従度に関する値を算出し、判断部１６ｂが追従度に関する値に基づいて劣化兆候を判断することもできる。
【０１４５】
ここで、追従度とは、運転中の冷凍機の設定温度に対する追従の程度を意味する。
【０１４６】
以下、追従度に関する値について、いくつかの変形例を説明する。
【０１４７】
（４−１１−１）
算出部１６ａは、追従度に関する値として、冷凍機５０の運転中の出口温度の平均値と設定温度との差（平均値との差）を算出する。以下、図１９を用いて平均値との差を説明する。
【０１４８】
冷凍機５０の運転開始から出口温度は低下し、設定温度に到達する。この後、冷凍機５０は設定温度を保つように運転する。算出部１６ａは、運転期間中の出口温度の平均値を求め、求めた平均値から設定温度を減算し、図１９中、設定温度と平均値との差として示されている値を算出する。この値を平均値との差という。
【０１４９】
ここで、冷凍機５０に劣化兆候があると出口温度は、設定温度に到達するのに時間がかかり、設定温度に到達後も設定温度に対して振動する。さらに冷凍機５０が劣化すると、出口温度が設定温度に到達しない場合もある。また、出口温度が設定温度を超えてしまう場合がある。よって、平均値との差は、冷凍機５０に劣化兆候があると大きくなる。
【０１５０】
平均値との差を用いて冷凍機５０の劣化兆候の有無を判断する方法について説明する。ユーザーは、劣化兆候がない冷凍機５０の時系列データに基づいて、劣化兆候の有無を判断する基準となる平均値との差（基準平均値との差）を、記憶部１５に記憶させる。次に、判断部１６ｂは、算出部１６ａが算出した平均値との差と基準平均値との差を比較し、基準平均値との差より平均値との差が大きい場合、または、平均値との差が基準平均値との差を基準とした正常領域値を外れた場合に冷凍機５０に劣化兆候があると判断する。
【０１５１】
これにより、冷凍機５０の運転開始後から冷凍機５０を停止するまでの期間の冷凍機５０の運転状態についても劣化兆候を検知することができる。
【０１５２】
（４−１１−２）
また、算出部１６ａは、劣化兆候のない冷凍機５０の複数の平均値との差と、その平均値との差を算出した時系列データの冷凍機５０の運転開始時の温度差流量乗積値を図２０に示すようにマハラノビスの基準空間Ｆとすることができる。算出部１６ａは、劣化兆候の有無を判断する冷凍機５０の平均値との差と温度差流量乗積値のデータを図２０のｆ点とする。そして、算出部１６ａは図２０中のｆ点とマハラノビスの基準空間Ｆとのマハラノビス距離を算出する。
【０１５３】
判断部１６ｂは、算出部１６ａが算出したマハラノビス距離が記憶部１５に記憶させた所定の値よりも大きい場合に、乖離量の積分値を求めた冷凍機５０を劣化兆候があると判断することができる。
【産業上の利用可能性】
【０１５４】
本発明は、冷媒または冷温水を吐出する熱源機において吐出する冷媒または冷温水の温度より劣化兆候を検知することに、特に有用である。
【符号の説明】
【０１５５】
１０  冷凍機劣化兆候検知装置
３２  空調機
４１  １次ポンプ
５０  冷凍機
６０  温度センサ
７１  上流側配管
７２  下流側配管
８０  信号線
１００  熱源システム１００
２００  エネルギー管理システム
【先行技術文献】
【特許文献】
【０１５６】
【特許文献１】特開平１１−１９３９７７号公報
【図１】
2014066452.tif 000003
【図２】
2014066452.tif 000004
【図３】
2014066452.tif 000005
【図４】
2014066452.tif 000006
【図５】
2014066452.tif 000007
【図６】
2014066452.tif 000008
【図７】
2014066452.tif 000009
【図８】
2014066452.tif 000010
【図９】
2014066452.tif 000011
【図１０】
2014066452.tif 000012
【図１１】
2014066452.tif 000013
【図１２】
2014066452.tif 000014
【図１３】
2014066452.tif 000015
【図１４】
2014066452.tif 000016
【図１５】
2014066452.tif 000017
【図１６】
2014066452.tif 000018
【図１７】
2014066452.tif 000019
【図１８】
2014066452.tif 000020
【図１９】
2014066452.tif 000021
【図２０】
2014066452.tif 000022