WO2023053673A1

WO2023053673A1 - 空気調和機及び空気調和システム

Info

Publication number: WO2023053673A1
Application number: PCT/JP2022/027912
Authority: WO
Inventors: 聡司奥村
Original assignee: 株式会社富士通ゼネラル
Priority date: 2021-09-29
Filing date: 2022-07-15
Publication date: 2023-04-06
Also published as: CN117980670A; EP4411289A1; AU2022357654A1; JP2023049949A; JP7380663B2

Abstract

空気調和機は、圧縮機、室外熱交換器及び膨張弁を有する室外機と、室内熱交換器を有する室内機とを有する。空気調和機は、室外機と室内機とが冷媒配管で接続されて形成される冷媒回路を有し、室内熱交換器を圧縮機において圧縮される冷媒の凝縮器として、かつ、室外熱交換器を前記室内熱交換器において凝縮される冷媒の蒸発器として機能させる暖房運転を少なくとも行うことができる。空気調和機は、少なくとも暖房運転における空気調和機の運転状態量を用いて冷媒回路に残存する冷媒量を推定する推定部を有する。推定部は、冷媒回路に残存する冷媒量の範囲に対応させた異なる複数の推定モデルを含み、複数の推定モデルのうち少なくとも一つは、運転状態量として室内側熱交換器の出口における冷媒の過冷却度を用いる。その結果、所望のタイミングでかつ冷媒の残存量に左右されないで冷媒量を判定できる。

Description

空気調和機及び空気調和システム

　本発明は、空気調和機の冷媒回路内に充填されている冷媒量の不足量（又は残存量）を推定する機能、特に、熱源側ユニット（以下、室外機ともいう）と利用側ユニット（以下、室内機ともいう）とが冷媒連絡配管を介して接続されたセパレートタイプの空気調和機の冷媒回路内に充填されている冷媒量の不足量（又は残存量）を推定する機能を持つ空気調和機または空気調和システムに関する。

　冷媒回路で検知できる運転状態量を用いて冷媒量の適否を判定する空気調和機が提案されている。特許文献１では、例えば、利用側ユニットを冷房運転し利用側熱交換器出口の過熱度が正値（利用側熱交換器出口のガス冷媒が過熱状態）になるように利用側膨張弁を制御しつつ利用側熱交換器の蒸発圧力が所定値になるように圧縮機の運転容量を制御する冷媒量判定運転モード（以下、デフォルト状態ともいう）で熱源側熱交換器出口の過冷却度を用いて冷媒量の適否を判定している。

特開２００６－２３０７２号公報

　空気調和機では、過冷却度等の運転状態量を用いて冷媒量の適否を判定する場合、冷媒回路を上述したデフォルト状態にする必要がある。そして、冷媒量の適否を判定する時点での過冷却度を規定量の冷媒が充填された直後の過冷却度と比較する。過冷却度の値を比較した結果、冷媒量の適否を判定する時点での過冷却度が小さくなっている場合は、冷媒量が少ない状態（適切ではない状態）であると判定する。

　しかしながら、冷媒量の適否を判定する時点での外気温や室内温度などの外部環境が、季節や日射量等の違いによる影響で、規定量の冷媒が充填された直後の外部環境に一致するとは限らない。このため、冷媒量の適否を判定したいときに冷媒回路の状態をデフォルト状態に合わせる（例えば蒸発圧力を所定値に合わせる）ことが難しく、所望のタイミングでデフォルト状態に合わせることができずに冷媒量の適否を判定できない場合があった。また、上述した冷媒回路では、冷媒回路に残存する冷媒量によっても冷媒回路で検知できる運転状態量が変化してしまうため、十分な量の冷媒が充填されている場合とそうでない場合とで冷媒回路の状態が異なる。このため、熱源側熱交換器出口の過冷却度の大小により冷媒量の適否を判定する特許文献１の方法では、冷媒回路における冷媒の不足状態は判定できても、正確に冷媒の不足量（又は残存量）を判定することはできなかった。

　また、利用側ユニットを冷房運転する特許文献１の方法では、利用側ユニットを暖房運転した状態で過冷却度を用いて冷媒の不足量（又は残存量）を判定することもできない。

　本発明ではこのような問題に鑑み、利用側ユニットを暖房運転した状態でも、冷媒の残存量に左右されないで冷媒量の不足量（又は残存量）を判定できる空気調和機を提供することを目的とする。

　一つの態様の空気調和機は、圧縮機、室外熱交換器及び膨張弁を有する室外機と、室内熱交換器を有する室内機とを有する。空気調和機は、室外機と室内機とが冷媒配管で接続されて形成される冷媒回路を有し、室内熱交換器を圧縮機において圧縮される冷媒の凝縮器として、かつ、室外熱交換器を前記室内熱交換器において凝縮される冷媒の蒸発器として機能させる暖房運転を少なくとも行うことができる。空気調和機は、少なくとも暖房運転における空気調和機の運転状態量を用いて冷媒回路に残存する冷媒量を推定する推定部を有する。推定部は、冷媒回路に残存する冷媒量の範囲に対応させた異なる複数の推定モデルを含み、複数の推定モデルのうち少なくとも一つは、運転状態量として室内側熱交換器の出口における冷媒の過冷却度を用いる。

　一つの側面として、所望のタイミングで、かつ、冷媒の残存量に左右されないで冷媒量を判定できる。

図１は、本実施例の空気調和機の一例を示す説明図である。図２は、室外機及び室内機の一例を示す説明図である。図３Ａは、室外機の室外機制御部の一例を示すブロック図である。図３Ｂは、室内機の室内機制御部の一例を示すブロック図である。図４は、集中コントローラ内の制御回路の一例を示すブロック図である。図５は、空気調和機の冷媒変化の状態を示すモリエル線図である。図６Ａは、第１の冷房用推定モデルによる推定結果と第２の冷房用推定モデルによる推定結果との間をシグモイド曲線で補間しなかった場合の一例を示す説明図である。図６Ｂは、第１の冷房用推定モデルによる推定結果と第２の冷房用推定モデルによる推定結果との間のシグモイド曲線で補間した場合の一例を示す説明図である。図７Ａは、第１の暖房用推定モデルによる推定結果と第２の暖房用推定モデルによる推定結果との間をシグモイド曲線で補間しなかった場合の一例を示す説明図である。図７Ｂは、第１の暖房用推定モデルによる推定結果と第２の暖房用推定モデルによる推定結果との間のシグモイド曲線で補間した場合の一例を示す説明図である。図８は、センサ値編集処理の一例を示す説明図である。図９は、推定処理に関わる制御回路の処理動作の一例を示すフローチャートである。図１０は、重回帰分析処理に関わる制御回路の処理動作の一例を示すフローチャートである。図１１は、冷房運転時の室外熱交換機における冷媒出口側の冷媒過冷却度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図１２は、冷房運転時の吸入温度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図１３は、暖房運転時の室外機膨張弁の開度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図１４は、暖房運転時の室内機３の過冷却度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図１５は、吸入過熱度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図１６Ａは、暖房運転時の室外機膨張弁の開度のみを使用した第３の暖房用推定モデルの冷媒不足率毎の推定値の精度の関係を示す説明図である。図１６Ｂは、暖房運転時の室外機膨張弁の開度及び室内過冷却度等を使用した第３の暖房用推定モデルの冷媒不足率毎の推定値の精度の関係を示す説明図である。図１７は、実施例２の空気調和システムの一例を示す説明図である。

　以下、図面に基づいて、本願の開示する空気調和機等の実施例を詳細に説明する。尚、本実施例により、開示技術が限定されるものではない。また、以下に示す各実施例は、矛盾を起こさない範囲で適宜変形しても良い。

＜空気調和機の構成＞
　図１は、本実施例の空気調和機１の一例を示す説明図である。図１に示す空気調和機１は、１台の室外機２と、Ｎ台の室内機３と、各室内機３を個別に制御する個別制御手段としての個別コントローラ（不図示）と、室外機２及び室内機３の状態（例えば後述する運転情報など）の表示と制御を行う集中制御手段としての集中コントローラ７とを有する（Ｎは２上の自然数）。室外機２は、液管４及びガス管５で並列に各室内機３と接続する。そして、室外機２と室内機３とが液管４及びガス管５等の冷媒配管で接続することで、空気調和機１の冷媒回路６が形成されている。室内機３は個別コントローラによる使用者からの操作指示を受け付けて室内機３毎に空調運転を実行する。集中コントローラ７は、室外機２及び室内機３を含む空気調和機本体１Ａの状態を表示するモニタ部８０と、空気調和機本体１Ａを制御する制御回路７０を有する。

＜室外機の構成＞
　図２は、室外機２およびＮ台の室内機３の一例を示す説明図である。室外機２は、圧縮機１１と、四方弁１２と、室外熱交換器１３と、室外機膨張弁１４と、第１の閉鎖弁１５と、第２の閉鎖弁１６と、アキュムレータ１７と、室外機ファン１８と、室外機制御部１９とを有する。これら圧縮機１１、四方弁１２、室外熱交換器１３、室外機膨張弁１４、第１の閉鎖弁１５、第２の閉鎖弁１６及びアキュムレータ１７を用いて、以下で詳述する各冷媒配管で相互に接続されて冷媒回路６の一部を成す室外側冷媒回路を形成する。

　圧縮機１１は、例えば、インバータにより回転数が制御される図示しないモータの駆動に応じて、運転容量を可変できる高圧容器型の能力可変型圧縮機である。圧縮機１１は、その冷媒吐出側と四方弁１２の第１のポート１２Ａとの間を吐出管２１で接続している。また、圧縮機１１は、その冷媒吸入側とアキュムレータ１７の冷媒流出側との間を吸入管２２で接続している。

　四方弁１２は、冷媒回路６における冷媒の流れる方向を切替えるための弁であって、第１～第４のポート１２Ａ～１２Ｄを備えている。第１のポート１２Ａは、圧縮機１１の冷媒吐出側との間を吐出管２１で接続している。第２のポート１２Ｂは、室外熱交換器１３の一方の冷媒出入口との間を室外冷媒管２３で接続している。第３のポート１２Ｃは、アキュムレータ１７の冷媒流入側との間を室外冷媒管２６で接続している。そして、第４のポート１２Ｄは、第２の閉鎖弁１６との間を室外ガス管２４で接続している。

　室外熱交換器１３は、冷媒と、室外機ファン１８の回転により室外機２の内部に取り込まれた外気とを熱交換させる。室外熱交換器１３は、その一方の冷媒出入口と四方弁１２の第２のポート１２Ｂとの間を室外冷媒管２６で接続している。室外熱交換器１３は、その他方の冷媒出入口と第１の閉鎖弁１５との間を室外液管２５で接続している。室外熱交換器１３は、空気調和機１が冷房運転を行う場合に凝縮器として機能し、空気調和機１が暖房運転を行う場合に蒸発器として機能する。

　室外機膨張弁１４は、室外液管２５に設けられており、図示しないパルスモータで駆動する電子膨張弁である。室外機膨張弁１４は、パルスモータに与えられるパルス数に応じて開度が調整されることで、室外熱交換器１３に流入する冷媒量、又は、室外熱交換器１３から流出する冷媒量を調整するものである。室外機膨張弁１４の開度は、空気調和機１が暖房運転を行っている場合、圧縮機１１の冷媒吸入側の冷媒過熱度が目標吸入冷媒過熱度となるように調整される。また、室外機膨張弁１４の開度は、空気調和機１が冷房運転を行っている場合、全開とされる。

　アキュムレータ１７は、その冷媒流入側と四方弁１２の第３のポート１２Ｃとの間を室外冷媒管２６で接続している。更に、アキュムレータ１７は、その冷媒流出側と圧縮機１１の冷媒流入側との間を吸入管２２で接続している。アキュムレータ１７は、室外冷媒管２６からアキュムレータ１７の内部に流入した冷媒をガス冷媒と液冷媒とに分離し、ガス冷媒のみを圧縮機１１に吸入させる。

　室外機ファン１８は、樹脂材で形成されており、室外熱交換器１３の近傍に配置されている。室外機ファン１８は、図示しないファンモータの回転に応じて、図示しない吸込口から室外機２の内部へ外気を取り込み、室外熱交換器１３において冷媒と熱交換した外気を図示しない吹出口から室外機２の外部へ放出する。

　また、室外機２には、複数のセンサが配置されている。吐出管２１には、圧縮機１１から吐出される冷媒の圧力、すなわち吐出圧力を検出する吐出圧センサ３１と、圧縮機１１から吐出された冷媒の温度、すなわち吐出温度を検出する吐出温度センサ３２とが配置されている。室外冷媒管２６のアキュムレータ１７の冷媒流入口近傍には、圧縮機１１に吸入される冷媒の圧力である吸入圧力を検出する吸入圧力センサ３３と、圧縮機１１に吸入される冷媒の温度を検出する吸入温度センサ３４とが配置されている。

　室外熱交換器１３と室外機膨張弁１４との間の室外液管２５には、室外熱交換器１３に流入する冷媒の温度、又は、室外熱交換器１３から流出する冷媒の温度を検出するための冷媒温度センサ３５が配置されている。そして、室外機２の図示しない吸込口付近には、室外機２の内部に流入する外気の温度、すなわち外気温度を検出する外気温度センサ３６が配置されている。

　図３Ａは、室外機２の室外機制御部１９の一例を示すブロック図である。図３Ａに示す室外機制御部１９は、室外側検出部１９Ａと、室外側記憶部１９Ｂと、室外側制御部１９Ｃとを有する。室外側検出部１９Ａは、運転状態量のうち室外機２側の運転状態量である室外側運転状態量を検出する。室外側検出部１９Ａは、室外機２の各センサである。室外側記憶部１９Ｂは、室外側検出部１９Ａで検出した室外側検出結果を記憶する。室外側検出結果は、室外機２の各センサの検出結果及び各センサの検出時刻を含む。室外側制御部１９Ｃは、室外機２の各部の動作を制御する。室外側制御部１９Ｃは、室外側記憶部１９Ｂに記憶中の室外側検出結果を集中コントローラ７に転送する場合、室外側検出結果の内、直前の検出時刻のセンサ値に対して変化がある場合にのみ、そのときのセンサ値を室外側検出結果として集中コントローラ７に転送する。また、室外側制御部１９Ｃは、直前の検出時刻のセンサ値に対して変化がない場合、室外側検出結果の集中コントローラ７へ転送しない。

＜室内機の構成＞
　図２に示すように、室内機３は、室内熱交換器５１と、室内機膨張弁５２と、液管接続部５３と、ガス管接続部５４と、室内機ファン５５と、室内機制御部６５とを有する。これら室内熱交換器５１、室内機膨張弁５２、液管接続部５３及びガス管接続部５４は、後述する各冷媒配管で相互に接続されて、冷媒回路６の一部を成す室内機冷媒回路を構成する。

　室内熱交換器５１は、冷媒と、室内機ファン５５の回転により図示しない吸込口から室内機３の内部に取り込まれた室内空気とを熱交換させる。室内熱交換器５１は、その一方の冷媒出入口と液管接続部５３との間を室内液管５６で接続している。また、室内熱交換器５１は、その他方の冷媒出入口とガス管接続部５４との間を室内ガス管５７で接続している。室内熱交換器５１は、空気調和機１が暖房運転を行う場合、凝縮器として機能する。これに対して、室内熱交換器５１は、空気調和機１が冷房運転を行う場合、蒸発器として機能する。

　室内機膨張弁５２は、室内液管５６に設けられており、電子膨張弁である。室内熱交換器５１が蒸発器として機能する場合、すなわち、室内機３が冷房運転を行う場合、室内機膨張弁５２の開度は、室内熱交換器５１の冷媒出口（ガス管接続部５４側）での冷媒過熱度が目標冷媒過熱度となるように調整される。また、室内熱交換器５１が凝縮器として機能する場合、すなわち室内機３が暖房運転を行う場合、室内機膨張弁５２の開度は、室内熱交換器５１の冷媒出口（液管接続部５３側）での冷媒過冷却度が目標冷媒過冷却度となるように調整される。ここで、目標冷媒過熱度や目標冷媒過冷却度とは、室内機３で十分な冷房能力あるいは暖房能力を発揮するのに必要な冷媒過熱度および冷媒過冷却度である。

　室内機ファン５５は、樹脂材で形成されており、室内熱交換器５１の近傍に配置されている。室内機ファン５５は、図示しないファンモータによって回転することで、図示しない吸込口から室内機３の内部に室内空気を取り込み、室内熱交換器５１において冷媒と熱交換した室内空気を図示しない吹出口から室内へ放出する。

　室内機３には各種のセンサが設けられている。室内液管５６には、室内熱交換器５１と室内機膨張弁５２との間に、室内熱交換器５１に流入する冷媒の温度、又は室内熱交換器５１から流出する冷媒の温度を検出する液側冷媒温度センサ６１が配置されている。室内ガス管５７には、室内熱交換器５１から流出又は室内熱交換器５１に流入する冷媒の温度を検出するガス側温度センサ６２が配置されている。室内機３の図示しない吸込口付近には、室内機３の内部に流入する室内空気の温度、すなわち吸込温度を検出する吸込温度センサ６３が配置されている。

　図３Ｂは、室内機３の室内機制御部６５の一例を示すブロック図である。図３Ｂに示す室内機制御部６５は、室内側検出部６５Ａと、室内側記憶部６５Ｂと、室内側制御部６５Ｃとを有する。室内側検出部６５Ａは、運転状態量のうち室内機３側の運転状態量である室内側運転状態量を検出する。室内側検出部６５Ａは、室内機３内の各センサである。室内側記憶部６５Ｂは、室内側検出部６５Ａで検出した室内側検出結果を記憶する。室内側検出結果は、室内機３内の各センサの検出結果及び各センサの検出時刻を含む。室内側制御部６５Ｃは、個別コントローラ（不図示）から使用者の運転指示を受ける。運転指示を受けた室内側制御部６５Ｃは、指示内容に応じて室内機３の各部の動作を制御する。また、室内側制御部６５Ｃは、室内側記憶部６５Ｂに記憶中の室内側検出結果を、室外機制御部１９を介して集中コントローラ７に転送する。この場合、室内側制御部６５Ｃは、室内側検出結果の内、直前の検出時刻のセンサ値に対して変化がある場合にのみ、そのときのセンサ値を室内側検出結果として集中コントローラ７に転送する。また、室内側制御部６５Ｃは、直前の検出時刻のセンサ値に対して変化がない場合、室外側検出結果を集中コントローラ７へ転送しない。

＜冷媒回路の動作＞
　次に、本実施形態における空気調和機１の空調運転時の冷媒回路６における冷媒の流れや各部の動作について説明する。尚、図１における矢印は暖房運転時の冷媒の流れを示している。

　空気調和機１が暖房運転を行う場合、四方弁１２は、第１のポート１２Ａと第４のポート１２Ｄとが連通し、第２のポート１２Ｂと第３のポート１２Ｃとが連通するように切替えている。これにより、冷媒回路６は、各室内熱交換器５１が凝縮器として機能し、室外熱交換器１３が蒸発器として機能する暖房サイクルとなる。尚、説明の便宜上、暖房運転時の冷媒の流れは、図２に示す実線矢印で表記する。

　冷媒回路６が上記の状態で圧縮機１１が駆動すると、圧縮機１１から吐出された冷媒は、吐出管２１を流れて四方弁１２に流入し、四方弁１２から室外ガス管２４を流れて、第２の閉鎖弁１６を介してガス管５へと流入する。ガス管５を流れる冷媒は、各ガス管接続部５４を介して各室内機３に分流する。各室内機３に流入した冷媒は、各室内ガス管５７を流れて各室内熱交換器５１に流入する。各室内熱交換器５１に流入した冷媒は、各室内機ファン５５の回転により各室内機３の内部に取り込まれた室内空気との間で熱交換することで凝縮する。つまり、各室内熱交換器５１が凝縮器として機能し、各室内熱交換器５１で冷媒によって加熱された室内空気が図示しない吹出口から室内に吹き出されることで、各室内機３が設置された室内の暖房が行われる。

　各室内熱交換器５１から各室内液管５６に流入した冷媒は、各室内熱交換器５１の冷媒出口側での冷媒過冷却度が目標冷媒過冷却度となるように開度が調整された各室内機膨張弁５２を通過して減圧される。ここで、目標冷媒過冷却度は、各室内機３で要求される冷房能力に基づいて定められるものである。

　各室内機膨張弁５２で減圧された冷媒は、各室内液管５６から各液管接続部５３を介して液管４に流出する。液管４で合流した冷媒は、第１の閉鎖弁１５を介して室外機２に流入する。室外機２の第１の閉鎖弁１５に流入した冷媒は、室外液管２５を流れ、室外機膨張弁１４を通過して減圧される。室外機膨張弁１４で減圧された冷媒は、室外液管２５を流れて室外熱交換器１３に流入し、室外機ファン１８の回転によって室外機２の図示しない吸込口から流入した外気と熱交換を行って蒸発する。室外熱交換器１３から室外冷媒管２６へと流出した冷媒は、四方弁１２、室外冷媒管２６、アキュムレータ１７及び吸入管２２の順に流入し、圧縮機１１に吸入されて再び圧縮され、四方弁１２の第１のポート１２Ａ及び第４のポート１２Ｄ経由で室外ガス管２４に流出する。

　また、空気調和機１が冷房運転を行う場合、四方弁１２は、第１のポート１２Ａと第２のポート１２Ｂとが連通し、第３のポート１２Ｃと第４のポート１２Ｄとが連通するように切替えている。これにより、冷媒回路６は、各室内熱交換器５１が蒸発器として機能し、室外熱交換器１３が凝縮器として機能する冷房サイクルとなる。尚、説明の便宜上、冷房運転時の冷媒の流れは、図２に示す破線矢印で表記する。

　冷媒回路６の状態で圧縮機１１が駆動すると、圧縮機１１から吐出された冷媒は、吐出管２１を流れて四方弁１２に流入し、四方弁１２から室外冷媒管２６を流れて、室外熱交換器１３に流入する。室外熱交換器１３に流入した冷媒は、室外機ファン１８の回転により室外機２の内部に取り込まれた室外空気との間で熱交換することで凝縮する。つまり、室外熱交換器１３が凝縮器として機能し、室外熱交換器１３で冷媒によって加熱された室内空気が図示しない吹出口から室外に吹き出す。

　室外熱交換器１３から室外液管２５へと流入した冷媒は、開度が全開とされている室外機膨張弁１４を通過して減圧される。室外機膨張弁１４で減圧された冷媒は、第１の閉鎖弁１５を介して液管４を流れて各室内機３に分流する。各室内機３に流入した冷媒は、各液管接続部５３を通じて室内液管５６を流れて室内熱交換器５１の冷媒出口で冷媒過冷却度が目標冷媒過冷却度となる開度に調整された室内機膨張弁５２を通過して減圧される。室内機膨張弁５２で減圧された冷媒は、室内液管５６を流れて室内熱交換器５１に流入し、室内機ファン５５の回転によって室内機３の図示しない吸入口から流入した室内空気と熱交換を行って蒸発する。つまり、各室内熱交換器５１が蒸発器として機能し、各室内熱交換器５１で冷媒によって冷却された室内空気が図示しない吹出口から室内に吹き出されることで、各室内機３が設置された室内の冷房が行われる。

　室内熱交換器５１からガス管接続部５４を介してガス管５へ流れる冷媒は、室外機２の第２の閉鎖弁１６を介して室外ガス管２４に流れて四方弁１２の第４のポート１２Ｄに流入する。四方弁１２の第４のポート１２Ｄに流入した冷媒は、第３のポート１２Ｃからアキュムレータ１７の冷媒流入側に流入する。アキュムレータ１７の冷媒流入側から流入した冷媒は、吸入管２２を介して流入し、圧縮機１１に吸入されて再び圧縮されることになる。

　暖房運転時に、室内熱交換器５１は圧縮機１１において圧縮される冷媒の凝縮器として、かつ、室外熱交換器１３は室内熱交換器５１において凝縮される冷媒の蒸発器として機能する。

＜集中コントローラ内の制御回路＞
　集中コントローラ７内の制御回路７０は、空気調和機１全体を制御する。図４は、集中コントローラ７内の制御回路７０の一例を示すブロック図である。制御回路７０は、取得部７１と、通信部７２と、記憶部７３と、制御部７４とを有する。取得部７１は、前述した各種センサのセンサ値を取得する。取得部７１は、室外機２内の吐出圧センサ３１、吐出温度センサ３２、吸入圧力センサ３３、吸込温度センサ６３、冷媒温度センサ３５及び外気温度センサ３６のセンサ値を取得する。更に、取得部７１は、各室内機３の液側冷媒温度センサ６１、ガス側温度センサ６２及び吸込温度センサ６３のセンサ値を取得する。

　通信部７２は、室外機２や各室内機３の通信部と通信する通信インタフェースである。記憶部７３は、例えば、フラッシュメモリであって、室外機２の制御プログラムや各種センサからの検出信号に対応した検出値等の運転状態量、室外機２の運転情報（例えば、運転／停止等の情報、圧縮機１１や室外機ファン１８の駆動状態等を含む）、各室内機３から送信される運転情報（例えば、運転／停止等の情報、冷房／暖房等の運転モード等を含む）、室外機２の定格能力及び各室内機３の要求能力、などを記憶する。

　本実施例では、記憶部７３は冷媒回路６に残存する冷媒量を推定する推定モデルを記憶している。本実施例では、冷媒回路６に残存する冷媒量として、例えば相対的な冷媒量を用いている。具体的には、本実施例の記憶部７３は冷媒回路６の冷媒不足率（冷媒回路６に規定量の冷媒が充填されている状態を冷媒充填率１００％の状態としたとき、この規定量からの減少分を指す。以下、同様）を推定する推定モデルを記憶している。記憶部７３が記憶する推定モデルは、例えば冷媒不足率が低い範囲（残存する冷媒量が多い範囲）に対応させた第１の冷房用推定モデル７３Ａを含む。また、記憶部７３が記憶する推定モデルは、例えば冷媒不足率が高い範囲（残存する冷媒量が少ない範囲）に対応させた第２の冷房用推定モデル７３Ｂを含む。また、記憶部７３が記憶する推定モデルは、例えば前記第１の冷房用推定モデル７３Ａと第２の冷房用推定モデル７３Ｂとを組み合わせた第３の冷房用推定モデル７３Ｃを含む。また、記憶部７３が記憶する推定モデルは、例えば冷媒不足率が低い範囲（残存する冷媒量が多い範囲）に対応させた第１の暖房用推定モデル７３Ｄを含む。また、記憶部７３が記憶する推定モデルは、例えば冷媒不足率が高い範囲（残存する冷媒量が少ない範囲）に対応させた第２の暖房用推定モデル７３Ｅを含む。また、記憶部７３が記憶する推定モデルは、例えば、第１の暖房用推定モデル７３Ｄと第２の暖房用推定モデル７３Ｅとを組み合わせた第３の暖房用推定モデル７３Ｆを含む。

　制御部７４は、通信部７２を介して各種センサでの検出値を定期的（例えば、３０秒毎）に取り込み、各室内機３から送信される運転情報を含む信号が通信部７２を介して入力される。制御部７４は、これら入力された各種情報に基づいて、室外機膨張弁１４の開度調整や圧縮機１１の駆動制御を行う。更に、制御部７４は、上述した各推定モデルを用いて冷媒不足率を推定する推定部７４Ａを有する。

　推定部７４Ａは、冷媒回路６内の冷媒不足率の範囲に応じて異なる複数の推定モデルを用いて、例えば、暖房運転における空気調和機本体１Ａの運転状態量を用いて冷媒回路６に残存する冷媒量を推定する。推定部７４Ａは、室内機３のうち、少なくとも２台以上の室内機３の室内熱交換器５１を冷媒の凝縮器として機能させる場合に、凝縮器として機能する室内熱交換器５１の出口における冷媒の過冷却度を用いて、推定モデルで冷媒量を推定する。

　図５は、空気調和機１の冷凍サイクルを示すモリエル線図である。空気調和機１の冷房運転時は、室外熱交換器１３が凝縮器として機能し、室内熱交換器５１が蒸発器として機能する。また、空気調和機１の暖房運転時は、室外熱交換器１３が蒸発器として機能し、室内熱交換器５１が凝縮器として機能する。

　圧縮機１１は、蒸発器から流入する低温低圧のガス冷媒を圧縮して高温高圧のガス冷媒（図５の点Ｂの状態になった冷媒）を吐出する。尚、圧縮機１１が吐出するガス冷媒の温度が吐出温度であり、吐出温度は、吐出温度センサ３２で検出する。

　凝縮器は、圧縮機１１からの高温高圧のガス冷媒を空気と熱交換して凝縮させる。この際、凝縮器では、潜熱変化によってガス冷媒が全て液冷媒となった後は顕熱変化によって液冷媒の温度が低下して過冷却状態となる（図５の点Ｃの状態）。尚、ガス冷媒が潜熱変化で液冷媒へと変化している際の温度が高圧飽和温度であり、高圧飽和温度は吐出圧力センサ３１で検出した圧力値（図５に「ＨＰＳ」と表記している圧力値Ｐ２）に相当する温度である。凝縮器の出口における過冷却状態となっている冷媒の温度が熱交出口温度であり、空気調和機１の冷房運転時における熱交出口温度は、冷媒温度センサ３５で検出する。

　膨張弁は、凝縮器から流出した低温高圧の冷媒を減圧して、ガスと液とが混合した気液二相冷媒（図５の点Ｄの状態になった冷媒）となる。

　蒸発器は、流入した気液二相冷媒を空気と熱交換して蒸発させる。この際、蒸発器では、潜熱変化によって気液二相冷媒が全てガス冷媒となった後は顕熱変化によってガス冷媒の温度が上昇して過熱状態（図５の点Ａの状態）となり、圧縮機１１に吸入される。尚、液冷媒が潜熱変化でガス冷媒へと変化している際の温度が低圧飽和温度である。低圧飽和温度は、吸入圧力センサ３３で検出した圧力値（図５に「ＬＰＳ」と表記している圧力値Ｐ１）に相当する温度である。また、蒸発器で過熱されて圧縮機１１に吸入される冷媒の温度が吸入温度である。吸入温度は、吸入温度センサ３４で検出する。

　なお、凝縮器から流出する際に過冷却状態となっている冷媒の冷媒過冷却度は、高圧飽和温度から凝縮器として機能している熱交換器の冷媒出口における冷媒温度（上述した熱交出口温度）を減じて算出できる。また、蒸発器から流出する際に過熱状態となっている冷媒の冷媒過熱度は、低圧飽和温度から吸入温度を減じて算出できる。

＜推定モデルの構成＞
　推定モデルは、複数の運転状態量の内、任意の運転状態量（特徴量）を用いて、例えば回帰分析法の一種である重回帰分析法を用いて生成されている。複数のシミュレーション結果（数値計算により冷媒回路を再現して、残存する冷媒量の変化（例えば冷媒不足率０％、冷媒不足率１０％、冷媒不足率２０％・・・）に対して運転状態量がどのような値となるかを計算した結果）を重回帰分析法で分析した結果、複数の回帰式が得られた。この回帰式のうち、Ｐ値（生成した推定モデルの精度に運転状態量が与える影響度合いを示す値（所定の重みパラメータ））が小さく、かつ、補正値Ｒ２（生成した推定モデルの精度を示す値）が０．９以上１．０以下の間のできるだけ大きい値となる回帰式を推定モデルとして使用する。ここで、Ｐ値および補正値Ｒ２は、重回帰分析法で推定モデルを生成する際に、当該推定モデルの精度に関わる値であり、Ｐ値が小さいほど、また、補正値Ｒ２が１．０に近い値であるほど、生成された推定モデルの精度が高くなる。

　その結果、冷房時の冷媒不足率が０～３０％の場合では、例えば、運転状態量として冷媒過冷却度、外気温度、高圧飽和温度及び圧縮機１１の回転数を特徴量とする回帰式を推定モデルとする。冷房時の冷媒不足率が４０～７０％の場合では、例えば、吸入温度、外気温度及び圧縮機１１の回転数といった運転状態量を特徴量とする回帰式を推定モデルとする。

　暖房時の冷媒不足率が０～２０％の場合では、例えば、運転状態量として室外機膨張弁１４の開度、室内機３の過冷却度、圧縮機１１の回転数、を特徴量とする回帰式を推定モデルとする。本実施例の室外機２には搭載されていないが、サブクール熱交換器（以下、ＳＣ熱交換器ともいう）を搭載する場合には、運転状態量としてＳＣ熱交出口温度を特徴量としても良い。尚、室内機３の過冷却度は、暖房運転時に凝縮器として機能する室内熱交換器５１から流出する冷媒の冷媒過冷却度である。室内機３の過冷却度は、（室外機２の高圧飽和温度（圧縮機１１の吐出圧力センサ３１で検出した圧力値を温度変換した値）－室内熱交換器５１の熱交出口温度（液側冷媒温度センサ６１の検出温度））で算出する。ここで、室内機３の過冷却度は外気温や室内温度などの外的要因の影響も受けるため、外的要因を反映した運転状態量（外気温度、室内温度）を特徴量に含めれば、冷媒不足率の推定精度を高めることができる。

　また、暖房時の冷媒不足率が３０％～７０％の場合では、例えば、運転状態量として吸入冷媒過熱度（吸入温度から低圧飽和温度を減じて求められる）、室外機膨張弁１４の開度を特徴量とする回帰式を推定モデルとする。

　本実施例の推定モデルは、後述する６つの推定モデル（第１の冷房用推定モデル７３Ａと、第２の冷房用推定モデル７３Ｂと、第３の冷房用推定モデル７３Ｃと、第１の暖房用推定モデル７３Ｄと、第２の暖房用推定モデル７３Ｅと、第３の暖房用推定モデル７３Ｆ）を含む。本実施例では、これら各推定モデルは、後述するシミュレーション結果を用いて生成される。なお、これら推定モデルは、本実施例のように予め空気調和機１に保持（例えば、集中コントローラ７の記憶部７３に記憶）されていても良いし、または空気調和機１と接続されるサーバ１２０に保持されていても良い。

　第１の冷房用推定モデル７３Ａは、例えば冷媒不足率が０％～３０％のような低い範囲（残存する冷媒量が多い範囲（第１の範囲））の場合に有効な推定モデルであって、冷媒不足率を高精度に推定できる第１の回帰式である。第１の回帰式は、例えば、（α１×冷媒過冷却度）＋（α２×外気温度）＋（α３×高圧飽和温度）＋（α４×圧縮機１１の回転数）＋α５である。係数α１～α５は、推定モデル生成の際に決定されるものとする。制御部７４は、第１の回帰式に、取得部７１にて取得された現在の冷媒過冷却度、外気温度、高圧飽和温度及び圧縮機１１の回転数を代入することで、現時点での冷媒回路６の冷媒不足率を算出する。尚、冷媒過冷却度、外気温度、高圧飽和温度及び圧縮機１１の回転数を代入する理由は、第１の冷房用推定モデル７３Ａの生成時に使用した特徴量を使用するためである。冷媒過冷却度は、例えば、（高圧飽和温度－熱交出口温度）で算出できる。外気温度は、外気温度センサ３６で検出する。高圧飽和温度は、吐出圧力センサ３１で検出した圧力値を温度変換した値である。圧縮機１１の回転数は、圧縮機１１の図示しない回転数センサで検出する。

　第２の冷房用推定モデル７３Ｂは、例えば、冷媒不足率が４０％～７０％のような高い範囲（残存する冷媒量が少ない範囲（第２の範囲））の場合に有効な推定モデルであって、冷媒不足率を高精度に推定できる第２の回帰式である。第２の回帰式は、例えば、（α１１×吸入温度）＋（α１２×外気温度）＋（α１３×圧縮機１１の回転数）＋α１４である。係数α１１～α１４は、推定モデル生成の際に決定されるものとする。制御部７４は、第２の回帰式に、取得部７１にて取得された現在の吸入温度、外気温度及び圧縮機１１の回転数を代入することで、現時点での冷媒回路６の冷媒不足率を算出する。尚、吸入温度、外気温度及び圧縮機１１の回転数を代入する理由は、第２の冷房用推定モデル７３Ｂの生成時に使用した特徴量を使用するためである。吸入温度は、吸入温度センサ３４で検出する。外気温度は、外気温度センサ３６で検出する。圧縮機１１の回転数は、圧縮機１１の図示しない回転数センサで検出する。

　ところで、前述したように、第１の回帰式で求めることができる冷媒不足率は０％～３０％であり、第２の回帰式で求めることができる冷媒不足率は４０％～７０％である。この場合、冷媒不足率が３０％～４０％である場合は、第１の回帰式を用いると冷媒不足率は３０％と算出され、第２の回帰式を用いると冷媒不足率は４０％と算出される。つまり、冷媒不足率が３０％～４０％である場合に、冷媒不足率が３０％以下での寄与度の高い冷媒過冷却度、冷媒不足率が４０％以上での寄与度の高い吸入温度の何れも変化が小さく、有効な推定モデルを生成できない。従って、第１の回帰式あるいは第２の回帰式を用いると、図６Ａに示すようにどちらのモデルを使用するのかによって冷媒不足率が異なる。

　上述した第１の冷房用推定モデル７３Ａ及び第２の冷房用推定モデル７３Ｂは、冷媒回路６に残存する冷媒量に応じて切り換えて利用することができる。例えば、空気調和機１を設置した直後であれば、冷媒不足率はほぼゼロであると推定できるため、第１の冷房用推定モデル７３Ａを用いる。そして、第１の冷房用推定モデル７３Ａによって冷媒不足率が高まってきたことが確認された場合には、推定モデルを第２の冷房用推定モデル７３Ｂに切り替える。上記推定モデルの切り換えは、空気調和機１の制御部が行うこともできるし、手動で行うようにしてもよい。

　しかし、次に説明する第３の冷房用推定モデル７３Ｃを用いることで、推定モデルの切り換えを不要にすることができる。

　第３の冷房用推定モデル７３Ｃは、上記のような第１の回帰式あるいは第２の回帰式のいずれを使用しても冷媒不足率を推定できない範囲も含めて、冷媒不足率が０％～７０％の範囲をカバーできる冷房時冷媒不足率算出式である。第３の冷房用推定モデル７３Ｃは第１の冷房用推定モデル７３Ａと第２の冷房用推定モデル７３Ｂとを組み合わせて生成する。具体的には、図６Ｂに示すように、第３の冷房用推定モデル７３Ｃ（冷房時冷媒不足率算出式）は、第１の冷房用推定モデル７３Ａ（第１の回帰式）の推定結果である冷媒不足率と第２の冷房用推定モデル７３Ｂ（第２の回帰式）の推定結果である冷媒不足率との間を、シグモイド係数を使用したシグモイド曲線で連続的につなぐものである。より具体的には、冷房時冷媒不足率算出式は、（シグモイド係数×第１の回帰式で求めた冷媒不足率）＋（（１－シグモイド係数）×第２の回帰式で求めた冷媒不足率）である。制御部７４は、第１の回帰式および第２の回帰式に取得部７１にて取得された現在の運転状態量を代入してそれぞれ算出された冷媒不足率を冷房時冷媒不足率算出式に代入して、現時点での冷媒回路６の冷媒不足率を算出する。

　ここで、シグモイド係数の算出は、運転状態量のいずれかを用いる。本実施例では、サブクールが０となると第１の回帰式による結果がほぼ一定となってしまうことを考慮し、サブクールが５℃のときに、シグモイド係数が０．５となる計算式とした。

　ｐ＝１／（１＋ｅｘｐ-（ｓｃ－５））
　　ｐ:シグモイド係数
　　ｓｃ:サブクール値

　このようにシグモイド係数を決定して第３の冷房用推定モデル７３Ｃに用いることで、冷媒不足率が０％～３０％、つまり、冷媒不足率が第１の範囲であるときは、第３の冷房用推定モデル７３Ｃによる推定値において第１の冷房用推定モデル７３Ａの推定値が支配的となり、また、冷媒不足率が４０％～７０％、つまり、冷媒不足率が第２の範囲であるときは、第３の冷房用推定モデル７３Ｃによる推定値において第２の冷房用推定モデル７３Ｂの推定値が支配的となる。

　なお、シグモイド係数の算出は上述した方法に限らず、実際の冷媒不足率が３０％以上であるとき、つまり、実際の冷媒不足率が第１の範囲でないときは、第３の冷房用推定モデル７３Ｃによる推定値において第２の冷房用推定モデル７３Ｂの推定値が支配的となるように、また、実際の冷媒不足率が４０％以下であるとき、つまり、実際の冷媒不足率が第２の範囲でないときは、第３の冷房用推定モデル７３Ｃによる推定値において第１の冷房用推定モデル７３Ａの推定値が支配的となるように、シグモイド係数を決定すればよい。

　第１の暖房用推定モデル７３Ｄは、冷媒不足率が０％～２０％（残存する冷媒量が多い範囲（第３の範囲））の場合に有効な推定モデルであって、冷媒不足率を高精度に推定できる第４の回帰式である。第４の回帰式は、例えば、（α３１×室外機膨張弁１４の開度）＋（α３２×室内機３の過冷却度）＋（α３３×圧縮機１１の回転数）＋α３４である。係数α３１～α３４は、推定モデル生成の際に決定されるものとする。制御部７４は、第４の回帰式に、取得部７１にて取得された現在の室外機膨張弁１４の開度、室内機３の過冷却度圧縮機１１の回転数を代入することで、冷媒不足率を算出する。尚、室外機膨張弁１４の開度、室内機３の過冷却度、圧縮機１１の回転数を代入する理由は、暖房運転時における室外機膨張弁１４の開度、および室内機３の過冷却度が、冷媒不足量が少ない場合（例えば第３の範囲）における冷媒量の変化に影響を受ける運転状態量であり、圧縮機１１の回転数が、稼働している室内機台数に影響を受ける運転状態量であるからである。第１の暖房用推定モデル７３Ｄの生成時にこれらの運転状態量（特徴量）を使用する。室外機膨張弁１４の開度は、図示しないセンサで検出する。圧縮機１１の回転数は、圧縮機１１の図示しない回転数センサで検出する。なお、圧縮機１１の回転数は、室外側制御部から取得してもよい。室内機３の過冷却度は、例えば、（室外機２の高圧飽和温度－液側冷媒温度センサ６１の検出温度）で算出する。ここで、室内機３の過冷却度は外気温や室内温度などの外的要因の影響も受けるため、外的要因（外気温や室内温度など）を反映した運転状態量（外気温度、室内温度）を特徴量に含めれば、冷媒不足率の検知精度を高めることができる。例えば、外的要因を考慮した推定モデル（第４´の回帰式）は、（α３１´×室外機膨張弁１４の開度）＋（α３２´×室内機３の過冷却度）＋（α３３´×外気温度）＋（α３４´×ＳＣ熱交出口温度）＋（α３５´×圧縮機１１の回転数）＋（α３６´×室内温度）＋α３７´）となる。係数α３１´～α３７´は、推定モデル生成の際に決定されるものとする。外気温度は、外気温度センサ３６で検出する。室内温度は、図示しない室内温度センサで検出する。

　第２の暖房用推定モデル７３Ｅは、冷媒不足率が３０％～７０％（残存する冷媒量が少ない範囲（第４の範囲））の場合に有効な推定モデルであって、冷媒不足率を高精度に推定できる第５の回帰式である。第５の回帰式は、例えば、（α４１×吸入過熱度）＋（α４２×室外機膨張弁１４の開度）＋α４３である。係数α４１～α４３は、推定モデル生成の際に決定されるものとする。制御部７４は、第５の回帰式に、取得部７１にて取得された現在の吸入過熱度、室外機膨張弁１４の開度を代入することで、現時点での冷媒回路６の冷媒不足率を算出する。尚、吸入過熱度、室外機膨張弁１４の開度を代入する理由は、暖房運転時における吸入過熱度、室外機膨張弁１４の開度が、冷媒不足量が多い場合（例えば第４の範囲）における冷媒量の変化に影響を受ける運転状態量だからであり、第２の暖房用推定モデル７３Ｅの生成時に使用した特徴量を使用するためである。吸入過熱度は、例えば、（吸入温度（吸入温度センサ３４の検出値）－低圧飽和温度（吸入圧力センサ３３で検出した圧力値に相当する温度））で算出できる。室外機膨張弁１４の開度は、図示しないセンサで検出する。

　また、前述したように、第４の回帰式で求める冷媒不足率は例えば０％～２０％であり、第５の回帰式で求める冷媒不足率は例えば３０％～７０％である。この場合、冷媒不足率が２０％～３０％の範囲にある空気調和機１に第４の回帰式を用いると冷媒不足率は２０％と算出される。また、同様の空気調和機１に第５の回帰式を用いると冷媒不足率は３０％と算出される。つまり、空気調和機１の冷媒不足率が２０％～３０％の範囲にある場合は、運転状態量（冷媒不足量（冷媒不足率）が低い場合に冷媒量の変化に影響を受ける運転状態量（室外機膨張弁１４の開度及び室内機３の過冷却度）、冷媒不足量が多い（冷媒不足率が高い）場合に冷媒量の変化に影響を受ける運転状態量（室外機膨張弁１４の開度及び吸入過熱度））の何れも変化が小さく、２０％～３０％の間における冷媒不足率の変化を推定することが難しい。従って、第４の回帰式あるいは第５の回帰式を独立に用いると、空気調和機１の冷媒不足率が２０％～３０％の範囲では、図７Ａに示すようにどちらのモデルを使用するのかによって冷媒不足率が異なる点に注意が必要である。

　上述した第１の暖房用推定モデル７３Ｄ及び第２の暖房用推定モデル７３Ｅは、冷媒回路６に残存する冷媒量に応じて切り換えて利用することができる。例えば、空気調和機１を設置した直後であれば、冷媒不足率はほぼゼロであると推定できるため、第１の暖房用推定モデル７３Ｄを用いることができる。そして、第１の暖房用推定モデル７３Ｄによって冷媒不足率が高まってきたことが確認された場合には、推定モデルを第２の暖房用推定モデル７３Ｅに切り換える。上記推定モデルの切り換えは、空気調和機１の制御部が行うこともできるし、手動で行うようにしてもよい。

　しかし、次に説明する第３の暖房用推定モデル７３Ｆを用いることで、推定モデルの切り換えを不要にすることができる。第３の暖房用推定モデル７３Ｆは、上記のような第４の回帰式あるいは第５の回帰式のいずれを使用しても冷媒不足率を推定できない範囲も含めて、冷媒不足率が０％～７０％の範囲をカバーできる暖房時冷媒不足率算出式である。第３の暖房用推定モデル７３Ｆは第１の暖房用推定モデル７３Ｄと第２の暖房用推定モデル７３Ｅとを組み合わせて生成する。具体的には、図７Ｂに示すように、第３の暖房用推定モデル７３Ｆ（暖房時冷媒不足率算出式）は、第１の暖房用推定モデル７３Ｄ（第４の回帰式）の推定結果である冷媒不足率と第２の暖房用推定モデル７３Ｅ（第５の回帰式）の推定結果である冷媒不足率との間を、シグモイド係数を使用したシグモイド曲線で連続的に繋ぐものである。より具体的には、暖房時冷媒不足率算出式は、（シグモイド係数×第５の回帰式で求めた冷媒不足率）＋（（１－シグモイド係数）×第４の回帰式で求めた冷媒不足率）である。制御部７４は、第４の回帰式および第５の回帰式に取得部７１にて取得された現在の運転状態量を代入してそれぞれ算出された冷媒不足率を暖房時冷媒不足率算出式に代入して、現時点での冷媒回路６の冷媒不足率を算出する。

　ここで、シグモイド係数の算出は、冷房運転時と同様に運転状態量のいずれかを用いる。本実施例では、室外機膨張弁１４の開度を用いてシグモイド係数ｐを算出する。室外機膨張弁１４の開度は暖房運転時の冷媒不足率を推定する第４の回帰式と第５の回帰式とのいずれかに用いられる運転状態量である。例えばシグモイド係数ｐは、室外機膨張弁１４の開度Ｄが全閉の場合にＤ＝０、全開の場合にＤ＝１００として、下記計算式から算出される。下記計算式は、室外機膨張弁１４の開度が全開となると第４の回帰式による結果がほぼ一定となってしまうことを考慮し、室外機膨張弁１４の開度が９０のときに、シグモイド係数ｐが０．５となる計算式とした。

　ｐ＝１／（１＋ｅｘｐ-（Ｄ－９０））
　　ｐ:シグモイド係数
　　Ｄ:　室外機膨張弁１４の開度

　このようにシグモイド係数を決定して第３の暖房用推定モデル７３Ｆに用いることで、冷媒不足率が０％～２０％、つまり、冷媒不足率が第３の範囲であるときは、第３の暖房用推定モデル７３Ｆによる推定値において第１の暖房用推定モデル７３Ｄの推定値が支配的となり、また、冷媒不足率が３０％～７０％、つまり、冷媒不足率が第４の範囲であるときは、第３の暖房用推定モデル７３Ｆによる推定値において第２の暖房用推定モデル７３Ｅの推定値が支配的となる。

　なお、シグモイド係数の算出は上述した方法に限らず、実際の冷媒不足率が２０％以上であるとき、つまり、実際の冷媒不足率が第３の範囲でないときは、第３の暖房用推定モデル７３Ｆによる推定値において第２の暖房用推定モデル７３Ｅの推定値が支配的となるように、また、実際の冷媒不足率が３０％以下であるとき、つまり、実際の冷媒不足率が第４の範囲でないときは、第３の暖房用推定モデル７３Ｆによる推定値において第１の暖房用推定モデル７３Ｄの推定値が支配的となるように、シグモイド係数を決定すればよい。

　以上に説明したように、冷房運転時は、冷媒不足率に応じた回帰式（第１の回帰式、第２の回帰式）を用いて冷媒不足率を推定できる。また、第１の回帰式と第２の回帰式とを含んだ冷房時冷媒不足率算出式を使用して冷媒不足率を推定しても良い。第１の回帰式と第２の回帰式とを使い分ける場合には、例えば冷房時の冷媒過冷却度が第１の閾値（図６Ａ及び図６ＢのＴｖ１）より大きい値である場合は、第１の回帰式を選択する。また、冷房時の冷媒過冷却度が第１の閾値以下である場合は、第２の回帰式を選択する。冷房時の冷媒過冷却度が第１の閾値付近の値である場合は、いずれの回帰式を用いるかで冷媒不足率の推定値が不連続に変化する。一方、第１の回帰式と第２の回帰式とを含んだ冷房時冷媒不足率算出式を使用する場合には上記のような切り換えが不要になる。また、第１の回帰式と第２の回帰式とを含んだ冷房時冷媒不足率算出式を選択すれば、冷媒過冷却度が第１の閾値付近にある場合でも、冷房時の冷媒不足率の変化を連続的に推定できる。

　また、暖房運転時は、冷媒不足率に応じた回帰式（第４の回帰式、第５の回帰式）を用いて冷媒不足率を推定できる。また、第４の回帰式と第５の回帰式とを含んだ暖房時冷媒不足率算出式を使用して冷媒不足率を推定しても良い。第４の回帰式と第５の回帰式とを使い分ける場合には、例えば暖房時の室外機膨張弁１４の開度が第２の閾値未満（図７Ａ及び図７ＢのＴｖ２）の場合に第４の回帰式を選択する。また、暖房時の室外機膨張弁１４の開度が第２の閾値以上である場合は、第５の回帰式を選択する。暖房時の室外機膨張弁１４の開度が第２の閾値付近の値である場合は、いずれの回帰式を用いるかで冷媒不足率の推定値が不連続に変化する。一方、第４の回帰式と第５の回帰式とを含んだ暖房時冷媒不足率算出式を使用する場合には上記のような切り換えが不要になる。また、第４の回帰式と第５の回帰式とを含んだ暖房時冷媒不足率算出式を選択すれば、室外機膨張弁１４の開度が第２の閾値付近にある場合でも、暖房時の冷媒不足率の変化を連続的に推定できる。

＜推定処理の動作＞
　図９は、推定処理に関わる制御回路７０の処理動作の一例を示すフローチャートである。尚、制御回路７０は、事前に生成された第１の冷房用推定モデル７３Ａ、第２の冷房用推定モデル７３Ｂ、第３の冷房用推定モデル７３Ｃ、第１の暖房用推定モデル７３Ｄ、第２の暖房用推定モデル７３Ｅ、第３の暖房用推定モデル７３Ｆを保持しているものとする。図９において制御回路７０内の制御部７４は、取得部７１を通じて運転状態量を運転データとして収集する（ステップＳ１１）。制御部７４は、収集した運転データから任意の運転状態量を抽出するデータフィルタリング処理を実行する（ステップＳ１２）。制御部７４は、データクレンジング処理を実行する（ステップＳ１３）。制御部７４内の推定部７４Ａは、各回帰式又は各冷媒不足率算出式を用いて、現時点の冷媒回路６の冷媒不足率を算出し（ステップＳ１４）、図９に示す処理動作を終了する。

　データフィルタリング処理は、複数の運転状態量の全てを使用するのではなく、所定フィルタ条件に基づき、複数の運転状態量の内、冷媒不足率を算出するのに必要な一部の運転状態量のみを抽出する。生成された推定モデルの各回帰式や各冷媒不足率算出式に、データフィルタリング処理を行った（異常値や突出値を除いた）運転状態量を代入することで、より正確に冷媒不足率を推定できる。

　所定のフィルタ条件は、第１のフィルタ条件と、第２のフィルタ条件と、第３のフィルタ条件とを有する。第１のフィルタ条件は、例えば、空気調和機１の全運転モード共通に抽出するデータのフィルタ条件である。第２のフィルタ条件は、冷房運転時に抽出するデータのフィルタ条件である。第３のフィルタ条件は、暖房運転時に抽出するデータのフィルタ条件である。

　第１のフィルタ条件は、例えば、圧縮機１１の駆動状態、運転モードの識別、特殊運転の排除、取得した値における欠損値の排除、各回帰式の生成に際し与える影響の大きい運転状態量について変化量が大きい値の除外（変化量が小さい値の選択）等である。圧縮機１１の駆動状態は、冷媒不足率を推定するために必要な運転状態量である。冷媒不足率を推定するための条件として、圧縮機１１が安定して運転すること（冷媒回路６内の冷媒の循環量が安定していること）が必要である。従って、圧縮機１１の起動時等の過渡期（冷媒回路６内の冷媒の循環量が不安定な状態）に検出された運転状態量を除外する必要があり、そのような運転状態量を除外するためにデータフィルタリング処理が設けられる。

　運転モードの識別とは、冷房運転時及び暖房運転時に取得した運転状態量のみを抽出するためのフィルタ条件である。従って、除湿運転時や送風運転時に取得した運転状態量は除外される。特殊運転の排除とは、例えば、油回収運転や除霜運転といった冷房運転時や暖房運転時と比べて冷媒回路６の状態が大きく異なる特殊運転時に取得した運転状態量を除外するフィルタ条件である。欠損値の排除とは、冷媒不足率の判定に使用する運転状態量に欠損値があった場合、当該運転状態量を用いて各回帰式を生成すれば精度が落ちる可能性があるため、欠損値を含む運転状態量を除外するフィルタ条件である。

　各回帰式や各冷媒不足率算出式に代入する運転状態量について変化量が小さい値の選択とは、空気調和機１の運転状態が安定している状態（冷媒回路６内の冷媒の循環量が安定している状態）での運転状態量のみを抽出するフィルタ条件であり、各回帰式や各冷媒不足率算出式による推定精度を上げるために必要な条件である。尚、冷媒不足率を推定する際に推定精度に大きな影響を与える運転状態量とは、例えば、冷房運転時の冷媒不足率が低い場合（例えば０～３０％の場合）に使用する冷媒過冷却度、冷房運転時の冷媒不足率が高い場合（例えば４０～７０％の場合）に使用する吸入温度、暖房運転時の冷媒不足率が低い場合（例えば０％～２０％の場合）に使用する室内機３の過冷却度や、暖房運転時の冷媒不足率が高い場合（例えば３０％～７０％の場合）に使用する吸入過熱度等である。

　第２のフィルタ条件には、例えば、熱交出口温度の排除、サブクールの異常、吐出温度の異常等がある。

　熱交出口温度の排除は、外気温度センサ３６と熱交出口温度センサ３５とが近い場所に配置されていることにより、冷房運転時に熱交出口温度センサ３５で検出した熱交出口温度が外気温度センサ３６で検出した外気温度より低くなることがないことを考慮したフィルタ条件であり、外気温度より低い熱交出口温度を除外するフィルタ条件である。

　サブクール異常は、冷房負荷が極端に大きいあるいは小さいことに起因して異常に高いあるいは以上に低い冷媒過冷却度検出されたときにこれを除外するフィルタ条件である。吐出温度の異常は、冷房負荷が小さいことに起因して圧縮機１１に吸入される冷媒量が減少する所謂ガス欠状態時に検出した吐出温度を除外するフィルタ条件である。

　第３のフィルタ条件は、例えば、吐出温度の異常等である。暖房運転時に暖房負荷の大きさに起因して吐出温度が高くなって吐出温度保護制御が実行されると、例えば、圧縮機１１の回転数を低下させることで吐出温度が低下するため、このときに検出した吐出温度を除外するフィルタ条件である。

　データクレンジング処理は、取得した全て運転状態量を冷媒不足率の推定に使用するのではなく、誤った推定を行うおそれがある運転状態量を除外するための処理である。具体的には、取得した運転状態量を平滑化してノイズ抑制やデータ数制限等がある。データの平滑化によるノイズ抑制とは、該当区間の平均値を算出し、各モデルにおいて例えば冷媒過冷却度、吸入温度、吸入冷媒過熱度の移動平均をとることで、ノイズを抑える処理である。データ数制限とは、例えば、データ数が少ないものは信頼性が低いため排除する処理である。例えば、１日分の入力データをフィルタリング処理して残ったデータ数がＸ個以上であれば冷媒不足率の推定に使用、それより少なければ、その日のデータはすべて使用しない。つまり、データクレンジング処理では、推定モデルの各回帰式や各冷媒不足率算出式に異常値や突出値を除いた運転状態量を代入することで、より正確に冷媒不足率を推定できる。

＜センサ値編集処理＞
　また、室内機３のセンサで検出した運転状態量を推定モデルに使用する場合には様々な課題がある。例えば、室外機２に対して室内機３が複数台接続されている場合には、複数の室内機３の内、運転中の室内機３と、停止中の室内機３とが混在する場合がある。このため、この点を考慮して各室内機３のセンサで検出した運転状態量を用いて推定モデルを使用することになる。

　さらに、室内機３の過冷却度を用いて第１の暖房推定モデル７３Ｄを使用する場合、室内機３の過冷却度は、室内機３の液側冷媒温度センサ６１の検出温度及び、室外機２の高圧飽和温度を用いて計算することになる。尚、室外機２の高圧飽和温度は、室外機２内の吐出圧力センサ３１のセンサ値で換算した値である。

　しかしながら、運転中の各室内機３の液側冷媒温度センサ６１の検出温度及び室外機２の吐出圧力センサ３１のセンサ値は、室内温度や室外温度などに影響をうけて変動する。この場合、室内機３の過冷却度を正確に計算するためには、各センサ値（室内機３の液側冷媒温度センサ６１の検出温度及び室外機２の吐出圧力センサ３１の圧力値）の検出時刻がなるべく近いセンサ値（以下、同一検出時刻付近のセンサ値ともいう）を用いる必要がある。従って、同一検出時刻付近の液側冷媒温度センサ６１の検出温度及び吐出圧力センサ３１の圧力値を得る仕組みが必要となる。

　そこで、本実施例では、同一検出時刻付近の室外機２の吐出圧力センサ３１のセンサ値と室内機３の液側冷媒温度センサ６１のセンサ値とを対応付けて取得するセンサ値編集処理が必要となる。

　図８は、センサ値編集処理の一例を示す説明図である。図８に示すセンサ値編集処理は、例えば集中コントローラ７の制御回路７０が実行する処理である。尚、説明の便宜上、複数の室内機３の内、運転中の室内機３を３台として説明する。運転中の室内機３は、例えば、「室内機“＃１”」、「室内機“＃２”」及び「室内機“＃３”」と表記する。

　室外側制御部１９Ｃは、室外側記憶部１９Ｂに記憶中の室外側検出結果を集中コントローラ７に転送する。また、室内側制御部６５Ｃは、室内側記憶部６５Ｂに記憶中の室内側検出結果を、室外機制御部１９を介して集中コントローラ７に転送する。各室内機３や室外機２から集中コントローラ７への検出結果の転送は、検出結果（センサ値）に変化がある場合にのみ行われる。例えば、室外側制御部１９Ｃや室内側制御部６５Ｃが前回検出結果と今回検出結果とを比較して、変化があれば（例えば運転モードが変わる、運転オンオフが変わる、センサの温度が変わるなど）検出結果を集中コントローラ７に転送する。集中コントローラ７へ転送される検出結果には、各室内機３や室外機２において検出結果の変化が検出された時刻（検出時刻）が紐づけられる。

　本実施例では、集中コントローラ７の制御回路７０内の取得部７１は、室外機２からは吐出圧力センサ３１で検出したセンサ値及び、その検出時刻を取得する場合を例示する。また、各室内機３からは液側冷媒温度センサ６１で検出したセンサ値及び、その検出時刻を取得する場合を例示する。図８に示す左図はセンサ値編集処理前のセンサ値、右図はセンサ値編集処理後のセンサ値である。

　制御回路７０内の制御部７４は、図８の左図に示すように、検出時刻と検出時刻毎の各センサ値とを取得し、順次記憶する。各室内機３や室外機２は、通信のトラヒックを減らすために、検出結果が変化したときに、検出結果を集中コントローラ７に転送する。このため、集中コントローラ７が各室内機３や室外機２から検出結果を取得する間隔は、それぞれ不規則となる。従って、図８中の「センサ値変化」は、直前のセンサ値と比較して検出時刻のセンサ値に変化がある場合（集中コントローラ７が各室内機３や室外機２から検出結果を取得したこと）を示している。また、図８中の「変化なし」は直前のセンサ値と比較して検出時刻のセンサ値に変化がない（集中コントローラ７が各室内機３や室外機２から検出結果を取得していない）場合を示している。制御部７４は、左図の記憶内容を参照し、検出時刻毎の室外機２及び室内機３の各センサ値を認識できる。

　制御部７４は、左図の検出時刻毎の室外機２及び各室内機３のセンサ値に基づき、空気調和機１全体（ユニット）におけるデータセットを作成する。データセットは、一定の時間間隔（例えば５分間隔）で刻まれた時刻（例えば図８中の「加工後の時刻」のこと。以下、代表時刻ともいう。）と、当該代表時刻に紐づけられる各代表センサ値と、で構成される。例えば、所定の代表時刻から５分後の代表時刻までの間に検出されたセンサ値を代表時刻付近のセンサ値と判断し、代表時刻付近毎のセンサ値を順次編集してデータセットを作成し、記憶する。

　例えば、代表時刻「０：０５」のセンサ値（代表センサ値）を決定する場合を説明する。時刻「０：００」～時刻「０：０５」の間において、制御部７４は検出時刻「０：００」、「０：０１」及び「０：０３」でのセンサ値を取得している。一方、時刻「０：０２」、「０：０４」ではいずれのセンサ値にも変化がないため、制御部７４はセンサ値を取得していない。従って、制御部７４は、代表時刻「０：０５」付近の検出時刻「０：００」、「０：０１」及び「０：０３」の各センサ値を使用して代表センサ値を決定する。例えば制御部７４は、代表時刻「０：０５」の室外機２の代表センサ値を決定する場合、まず検出時刻「０：００」、「０：０１」及び「０：０３」の室外機２のセンサ値内に「センサ値変化」があるか否かを判定する。次に、制御部７４は、「センサ値変化」がある場合に検出時刻のセンサ値の内、例えば最も早い時刻の「０：００」の「センサ値変化」のセンサ値を代表時刻「０：０５」の室外機２の代表センサ値と決定する。同様に、制御部７４は、代表時刻「０：０５」の「室内機＃１」の代表センサ値を決定する場合、まず検出時刻「０：００」、「０：０１」及び「０：０３」の「室内機＃１」のセンサ値内に「センサ値変化」があるか否かを判定する。次に、制御部７４は、「センサ値変化」がある場合に検出時刻のセンサ値の内、例えば最も早い時刻の「センサ値変化」のセンサ値を代表時刻「０：０５」の「室内機＃１」の代表センサ値と決定する。同様に、制御部７４は、代表時刻「０：０５」の「室内機＃２」の代表センサ値を決定する場合、まず検出時刻「０：００」、「０：０１」及び「０：０３」の「室内機＃２」のセンサ値内に「センサ値変化」があるか否かを判定する。次に、制御部７４は、「センサ値変化」がある場合に検出時刻のセンサ値の内、例えば最も早い時刻の「センサ値変化」のセンサ値を代表時刻「０：０５」の「室内機＃２」の代表センサ値と決定する。同様に、制御部７４は、代表時刻「０：０５」の「室内機＃３」の代表センサ値を決定する場合、まず検出時刻「０：００」、「０：０１」及び「０：０３」の「室内機＃３」のセンサ値内に「センサ値変化」があるか否かを判定する。次に、制御部７４は、「センサ値変化」がある場合に検出時刻のセンサ値の内、例えば最も早い時刻の「センサ値変化」のセンサ値を代表時刻「０：０５」の「室内機＃３」の代表センサ値と決定する。

　制御部７４は、例えば、代表時刻「０：１０」の代表センサ値を決定する場合について説明する。制御部７４は、代表時刻「０：１０」付近の検出時刻「０：０６」及び「０：０９」の各センサ値を使用して代表センサ値を決定する。例えば、制御部７４は、例えば、代表時刻「０：１０」の室外機２の代表センサ値を決定する場合、先ず、検出時刻「０：０６」及び「０：０９」の室外機２のセンサ値内に「センサ値変化」があるか否かを判定する。次に、制御部７４は、「センサ値変化」がないため、「変化なし」のセンサ値として、直前の代表時刻「０：０５」の室外機２のセンサ値を「前のセンサ値」として代表時刻「０：１０」の室外機２の代表センサ値と決定する。同様に、制御部７４は、代表時刻「０：１０」の「室内機＃１」の代表センサ値を決定する場合、代表時刻「０：１０」付近の検出時刻「０：０６」及び「０：０９」の「室内機＃１」のセンサ値内に「センサ値変化」があるか否かを判定する。次に、制御部７４は、「センサ値変化」がないため、「変化なし」のセンサ値として、直前の代表時刻「０：０５」の「室内機＃１」の代表センサ値を「前のセンサ値」として代表時刻「０：１０」の「室内機＃１」の代表センサ値と決定する。また、制御部７４は、例えば、代表時刻「０：１０」の「室内機＃２」の代表センサ値を決定する場合、代表時刻「０：１０」付近の検出時刻「０：０６」及び「０：０９」の「室内機＃２」のセンサ値内に「センサ値変化」があるか否かを判定する。次に、制御部７４は、「センサ値変化」がある場合、検出時刻のセンサ値の内、例えば最も早い時刻の「センサ値変化」のセンサ値を代表時刻「０：１０」の「室内機＃２」の代表センサ値と決定する。また、制御部７４は、代表時刻「０：１０」の「室内機＃３」の代表センサ値を決定する場合、代表時刻「０：１０」付近の検出時刻「０：０６」及び「０：０９」の「室内機＃３」のセンサ値内に「センサ値変化」があるか否かを判定する。次に、制御部７４は、「センサ値変化」がある場合に検出時刻のセンサ値の内、例えば最も早い時刻の「センサ値変化」のセンサ値を代表時刻「０：１０」の「室内機＃３」の代表センサ値と決定する。

　制御部７４は、例えば、代表時刻「０：１５」の代表センサ値を決定する場合を説明する。時刻「０：１１」～時刻「０：１５」の間では、いずれのセンサ値にも変化がないため、制御部７４はセンサ値を取得していない。従って、代表時刻「０：１５」付近の各センサ値がない。この場合には、直前の代表時刻「０：１０」の各代表センサ値を、代表時刻「０：１５」の代表センサ値と決定する。

　例えば、代表時刻「０：３０」の代表センサ値を決定する場合について説明する。制御部７４は、代表時刻「０：３０」付近の検出時刻「０：２７」及び「０：２８」の各センサ値を使用して代表センサ値を決定する。例えば、制御部７４は、代表時刻「０：３０」の室外機２の代表センサ値を決定する場合、先ず、検出時刻「０：２７」及び「０：２８」の室外機２の代表センサ値内に「センサ値変化」があるか否かを判定する。次に、制御部７４は、「センサ値変化」がないため、「変化なし」のセンサ値として、直前の代表時刻「０：２５」の室外機２の代表センサ値を「前のセンサ値」として代表時刻「０：３０」の室外機２の代表センサ値と決定する。

　また、制御部７４は、例えば、代表時刻「０：３０」の「室内機＃１」の代表センサ値を決定する場合、代表時刻「０：３０」付近の検出時刻「０：２７」及び「０：２８」の「室内機＃１」のセンサ値内に「センサ値変化」があるか否かを判定する。次に、制御部７４は、「センサ値変化」がないため、「変化なし」のセンサ値として、直前の代表時刻「０：２５」の「室内機＃１」の代表センサ値を「前のセンサ値」として代表時刻「０：３０」の「室内機＃１」の代表センサ値と決定する。また、制御部７４は、例えば、代表時刻「０：３０」の「室内機＃２」の代表センサ値を決定する場合、代表時刻「０：３０」付近の検出時刻「０：２７」及び「０：２８」の「室内機＃２」のセンサ値内に「センサ値変化」があるか否かを判定する。制御部７４は、「センサ値変化」がある場合、例えば最も早い時刻の「センサ値変化」のセンサ値を代表時刻「０：３０」の「室内機＃２」の代表センサ値と決定する。また、制御部７４は、例えば、代表時刻「０：３０」の「室内機＃３」の代表センサ値を決定する場合、代表時刻「０：３０」付近の検出時刻「０：２７」及び「０：２８」の「室内機＃３」のセンサ値内に「センサ値変化」があるか否かを判定する。制御部７４は、「センサ値変化」がある場合、例えば最も早い時刻の「センサ値変化」のセンサ値を代表時刻「０：３０」の「室内機＃３」の代表センサ値と決定する。

　制御部７４は、代表時刻毎に室外機２及び各室内機３のセンサ値を編集し、編集後の室外機２及び各室内機３のセンサ値を代表センサ値として記憶する。尚、制御部７４は、編集後の室内機２及び室内機３のセンサ値以外の不要なセンサ値を記憶部から消去する。このようにして、空気調和機１の運転データが収集される。

　なお、収集された運転データは、図９に示すデータフィルタリング処理、データクレンジング処理が施された後、冷媒不足率の算出に利用される。

　例えば第３の暖房用推定モデル７３Ｆを用いて冷媒不足率を算出する場合について説明する。この場合、代表時刻における室外機２及び室内機３の代表センサ値を参照し、各室内機＃１、＃２及び＃３の代表センサ値の平均値を用いて過冷却度等が算出される。例えば制御部７４は、室外機２のセンサ値を吐出圧力センサ３１のセンサ値、室内機３のセンサ値を液側冷媒温度センサ６１のセンサ値とした場合、代表時刻付近の吐出圧力センサ３１のセンサ値及び、液側冷媒温度センサ６１のセンサ値を参照して各室内機３における代表センサ値を得る。そして、制御部７４は、代表時刻の吐出圧力センサ３１の代表センサ値に基づき、代表時刻の高圧飽和温度を算出するため、代表時刻の吐出圧力センサ３１の代表センサ値及び、各室内機３における液側冷媒温度センサ６１の代表センサ値の平均値に基づき、代表時刻毎の室内機３の過冷却度を算出する。そして、制御部７４は、算出した代表時刻の室内機３の過冷却度等と第３の暖房用推定モデル７３Ｆを用いて、代表時刻における冷媒回路６の冷媒不足率を算出できる。

＜回帰式の生成方法＞
　次に第１～第６の回帰式の生成に使用する特徴量について説明する。第１～第３の回帰式を使用する冷房運転時では、重回帰分析法により第１～第６の回帰式の生成を行う際に使用する特徴量として、例えば、冷媒過冷却度、外気温度、高圧飽和温度、圧縮機１１の回転数、吸入温度等の各運転状態量を用いる。そして、これら各運転状態量は、シミュレーションにより得た結果を使用する。また、第４～第６の回帰式を使用する暖房運転時では、重回帰分析の特徴量として、例えば、室内機３の過冷却度、室内温度、吸入過熱度、外気温度、圧縮機１１の回転数、室外機膨張弁１４の開度等の各運転状態量を用いる。そして、これら各運転状態量は、シミュレーションにより得た結果を使用する。

　具体的には、空気調和機１の設計段階で、一例として室内機３が４台運転している場合に外気温度を異ならせてシミュレーションを行い、特徴量と冷媒不足率との関係をシミュレーション毎に取得する。シミュレーションを行う際の条件としては、例えば、外気温度を２０℃、２５℃、３０℃、３５℃及び４０℃と変化させる。なお、シミュレーションを行うに際しては、外気温度の他のパラメータを加えてもよく、例えば、室内機３の運転台数を１～４台と異ならせてもよい。

　図１１は、冷房運転時の室外熱交換機における冷媒出口側の冷媒過冷却度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図１１に示す冷媒過冷却度は、冷媒不足率が０％～３０％までは右肩下がりで減少し、冷媒不足率が３０％から６０％までは変化なしとなっている。つまり、冷房運転時に冷媒不足率０～３０％である場合は、冷媒回路６における冷媒量の不足が冷媒過冷却度の値に大きな影響を与えるということである。なお、図１１において冷媒不足率が６０％以上であるときの冷媒過冷却度がマイナスの値となっているが、実際は冷媒過冷却度が０℃未満とはならないため、これはシミュレーションでのみ現れる値である。従って、冷媒不足率が６０％以上であるときの冷媒過冷却度は、回帰式の生成に使用しない。

　図１２は、冷房運転時の吸入温度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図１２に示す吸入温度は、冷媒不足率が４０～７０％のとき増加する傾向にある。つまり、冷房運転時の冷媒不足率が４０～７０％である場合は、冷媒回路６における冷媒量の不足が吸入温度の値に大きな影響を与えるということである。なお、図１２において冷媒不足率が７０％以上であるときの吸入温度はほとんど変化しないため、これ以上の冷媒不足率を吸入温度で推測するのは難しい。従って、冷媒不足率が７０％以上であるときの吸入温度は、回帰式の生成に使用しない。

　図１３は、暖房運転時の室外機膨張弁１４の開度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図１３に示す室外機膨張弁１４の開度は、冷媒不足率が０～２０％の場合に変化するのに対し、冷媒不足率が２０％を超えると、室外機膨張弁１４の開度の変化が概ね無くなる。つまり、暖房運転時の冷媒不足率が０～２０％である場合は、冷媒回路６における冷媒量の不足が室外機膨張弁１４の開度に大きな影響を与えるということである。なお、上述したように、冷媒不足率が２０％を超えると、室外機膨張弁１４の開度の変化が概ね無くなる。従って、冷媒不足率が２０％を超えたときの室外機膨張弁１４の開度は、回帰式の生成に使用しない。

　図１４は、暖房運転時の室内機３の過冷却度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図１４に示す室内機３の過冷却度は、冷媒不足率が０～３５％の場合に変化するのに対し、冷媒不足率が３５％を超えると、過冷却度の変化が概ね無くなる。つまり、暖房運転時の冷媒不足率が低い領域（例えば０～２０％）では、冷媒回路６における冷媒量の不足が室内機３の過冷却度に大きな影響を与えるということである。なお、上述したように、冷媒不足率が３５％を超えると、室内機３の過冷却度の変化が概ね無くなる。

　図１５は、吸入過熱度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図１５に示す吸入過熱度は、冷媒不足率が増加するときに吸入過熱度が大きくなる傾向にあり、冷媒不足率が３０％を超えると吸入過熱度が大きく上昇する。つまり、暖房運転時の冷媒不足率が高い領域（例えば３０％）では、冷媒回路６における冷媒量の不足が吸入過熱度に大きな影響を与えるということである。なお、図１５において冷媒不足率が３０％より小さいときの吸入過熱度の変化が緩やかであるため、これ以下の冷媒不足率を吸入過熱度で精度良く推測するのは難しい。従って、本実施例では、冷媒不足率が３０％より小さい時の吸入過熱度は、回帰式の生成に使用しない。

　次に、第１の暖房用推定モデルの運転状態量として暖房運転時の室外機膨張弁１４の開度のみを使用した場合における、第３の暖房用推定モデルの冷媒不足率毎の推定値の精度について説明する。図１６Ａは、第１の暖房用推定モデルの運転状態量として暖房運転時の室外機膨張弁１４の開度のみを使用した第３の暖房用推定モデルの冷媒不足率毎の推定値の精度の関係を示す説明図である。

　例えば、室外機膨張弁１４の開度のみを使用した第１の暖房用推定モデルの冷媒不足率０％～２０％の推定値の補正Ｒ２は０．２９である。尚、補正Ｒ２が“１”に近似する程、推定値の精度が高くなることを示す。第１の暖房用推定モデルの運転状態量として暖房運転時の室外機膨張弁１４の開度のみを使用した場合、第３の暖房用推定モデルでは、図１６Ａに示すように、冷媒不足率が０％～２０％の場合、冷媒不足率毎の推定値が理想値Ｘから離れやすく、推定値の精度が低下している。

　一方、第１の暖房用推定モデルの運転状態量として暖房運転時の室内機３の過冷却度のみを使用した場合、第１の暖房用推定モデルにおける冷媒不足率０％～２０％の推定値の補正Ｒ２は０．５１である。従って、室外機膨張弁１４の開度のみを使用した場合よりも、室内機３の過冷却度を使用した方が、第１の暖房用推定モデルの推定値の精度が高くなる。また、暖房運転時の室内機３の過冷却度に加えて圧縮機１１の回転数を使用すると、第１の暖房用推定モデルにおける冷媒不足率０％～２０％の推定値の補正Ｒ２は０．８０となり、推定値の精度がさらに高くなる。

　本実施例の暖房運転時の第３の暖房用推定モデル７３Ｆでは、第１の暖房用推定モデルの運転状態量として室内機３の過冷却度及び圧縮機１１の回転数に加え、室外機膨張弁１４の開度を使用する。特に、室内機３の過冷却度は、図１４に示すように、冷媒不足率０～２０％のときに大きく変動する。第３の暖房用推定モデル７３Ｆは、冷媒不足率が低い範囲の場合に運転状態量として室内機３の過冷却度も考慮することで、冷媒不足率の変化の検出精度を高めることができる。

　本実施例で用いる第１の暖房用推定モデル７３Ｄは、室内機３の過冷却度と圧縮機１１の回転数に加えて、室外機膨張弁１４の開度を運転状態量として使用するため、冷媒不足率０％～２０％時の推定値の補正Ｒ２は０．８２である。図１６Ｂは、本実施例の第３の暖房運転モデル７３Ｆの冷媒不足率毎の推定値の精度の関係を示す説明図である。本実施例の第３の暖房用推定モデル７３Ｆでは、図１６Ｂに示すように、冷媒不足率が０％～２０％の場合、冷媒不足率毎の推定値が理想値Ｘに近く、冷媒回路６に残存する冷媒量の推定精度が高くなる。なお、上述したように、室内機３の過冷却度は外気温や室内温度などの外的要因の影響も受けるため、外的要因（外気温や室内温度など）を反映した運転状態量（外気温度、室内温度）を特徴量に含めれば、冷媒不足率の検知精度を高めることができる。室内機３の過冷却度、圧縮機１１の回転数、及び室外機膨張弁１４の開度に加えて、さらに外気温度、室内温度を運転状態量として含む場合には、冷媒不足率０％～２０％時における第１の暖房用推定モデル７３Ｄの推定値の補正Ｒ２は０．９２である。

＜実施例１の効果＞
　実施例１の空気調和機１では、暖房運転時の冷媒不足率が低い範囲での冷媒不足率推定モデルである第４の回帰式を生成する場合は、室内機３の過冷却度を用いる。その結果、冷媒不足率が低い範囲（例えば０％～２０％のとき）において、冷媒不足率に応じて値が大きく変動する室内機３の過冷却度を使用するため、冷媒不足率が低い範囲でも暖房運転時に冷媒不足率の変化を安定して推定できる。

　空気調和機１は、暖房運転時の冷媒不足率が高い範囲での冷媒不足率推定モデルである第５の回帰式を生成する場合は、運転状態量として、圧縮機１１の吸入過熱度、室外機膨張弁１４の開度を用いて回帰分析法で生成される。その結果、冷媒不足率が高い範囲において、暖房運転時に冷媒不足率の変化を安定して推定できる。

　空気調和機１は、冷房用推定モデルと、冷房運転時の現在の運転状態量とを用いて、冷房運転時の冷媒不足率を推定すると共に、暖房用推定モデルと、暖房運転時の現在の運転状態量とを用いて、暖房運転時の冷媒不足率を推定する。その結果、運転状態毎に異なる推定モデルを使用することで、冷媒不足率を高精度に推定できる。

　空気調和機１は、第１の暖房用推定モデル７３Ｄと第２の暖房用推定モデル７３Ｅとをシグモイド曲線で繋いだ第３の暖房用推定モデル７３Ｆに現在の運転状態量を代入することで、暖房運転時の冷媒不足率を高精度に推定できる。

　第１の暖房用推定モデル７３Ｄは、運転状態量として室外機膨張弁１４の開度及び室内機３の過冷却度を用いて、冷媒不足率を推定する。その結果、空気調和機１は、暖房運転時に冷媒不足率を高精度に推定できる。

　第２の暖房用推定モデル７３Ｅは、運転状態量として圧縮機１１の吸入過熱度を用いて、冷媒不足率を推定する。その結果、空気調和機１は、暖房運転時に冷媒不足率を高精度に推定できる。

　第３の暖房用推定モデル７３Ｆは、第１の暖房用推定モデル７３Ｄの推定結果と第２の暖房用推定モデル７３Ｅの推定結果との間をシグモイド曲線で補間する。その結果、暖房運転時の冷媒不足率が０～７０％の範囲で、正確な冷媒不足率を推定できる。

　重回帰分析処理において、データフィルタリング処理及びデータクレンジング処理後の現在の運転状態量（センサ値）を推定モデルの各回帰式に代入する。本実施例では、推定モデルの各回帰式の生成は、シミュレーションで得た特徴量を用いており、シミュレーションで得た特徴量には異常な値や他と比べて突出して大きいあるいは小さい値は含まれていない。このような、異常値や突出値を含まない特徴量を用いて生成された推定モデルの各回帰式や各冷媒不足率算出式に、データフィルタリング処理及びデータクレンジング処理を行って異常値や突出値を除いた運転状態量を代入することで、より正確に冷媒不足率を推定できる。

　尚、以上に説明した実施例では、空気調和機１の設計段階で各運転状態量のシミュレーション結果を求め、学習機能を有するサーバなどの端末にシミュレーション結果を学習させて得られた推定モデルを制御回路７０が予め記憶している場合を例示した。これに代えて、空気調和機１との間を通信網１１０で接続するサーバ１２０が存在し、このサーバ１２０が第１～第６の回帰式を生成して空気調和機１に送信するようにしてもよい。この実施の形態につき、以下に説明する。

＜空気調和システムの構成＞
　図１７は、実施例２の空気調和システム１００の一例を示す説明図である。尚、実施例１の空気調和機１と同一の構成には同一符号を付すことで、その重複する構成及び動作の説明については省略する。図１７に示す空気調和システム１００は、空気調和機本体１Ａと、集中コントローラ７と、通信網１１０と、サーバ１２０とを有する。空気調和機本体１Ａは、圧縮機１１、室外熱交換器１３及び室外機膨張弁１４を有する室外機２と、室内熱交換器５１を有する室内機３とを有する。空気調和機本体１Ａは、室外機２と室内機３とが液管４及びガス管５等の冷媒配管で接続されて構成する冷媒回路６を備え、当該冷媒回路６に所定量の冷媒が充填される。集中コントローラ７は、空気調和機本体１Ａと、通信網１１０との間を通信で接続する。集中コントローラ７は、室外機２及び室内機３を含む空気調和機本体１Ａの状態を表示するモニタ部８０と、空気調和機本体１Ａ全体を制御する制御回路７０を有する。

　サーバ１２０は、推定部１２１と、送信部１２２とを有する。推定部１２１は、冷媒回路６に充填される冷媒の冷媒不足率の推定に関わる運転状態量を用いて重回帰分析法を用いて生成された推定モデルを使用して、冷媒不足率を推定する。尚、推定モデルは、例えば、実施例１で説明した第１の冷房用推定モデル７３Ａ、第２の冷房用推定モデル７３Ｂ、第３の冷房用推定モデル７３Ｃ、第１の暖房用推定モデル７３Ｄ、第２の暖房用推定モデル７３Ｅ及び第３の暖房用推定モデル７３Ｆを有する。送信部１２２は、推定部１２１にて推定した推定結果を通信網１１０経由で集中コントローラ７に送信する。集中コンロトローラ７内の制御回路７０は、受信した推定結果を用いて、使用者に対して空気調和機１の冷媒回路６における冷媒不足率を表示する。

＜実施例２の効果＞
　実施例２のサーバ１２０は、現在の運転状態量を用いて、冷媒不足率を推定する。その結果、使用者は集中コントローラ７を介して空気調和機１の冷媒不足率を確認することができる。

　また、本実施例では、冷媒回路６に残存する冷媒量を表すものとして相対的な冷媒量を推定する場合を説明した。具体的には、冷媒回路６に冷媒を充填した際の充填量（初期値）に対する、冷媒回路６から外部に漏洩した冷媒量の割合である冷媒不足率を推定して提供する場合を説明した。しかし、本発明はこれに限られるものではなく、推定した冷媒不足率に初期値を乗じて、冷媒回路６から外部に漏洩した冷媒量を提供するようにしてもよい。また、冷媒回路６から外部に漏洩した絶対的な冷媒量あるいは冷媒回路６に残留する絶対的な冷媒量を推定するようにしてもよい。冷媒回路６から外部に漏洩した絶対的な冷媒量あるいは冷媒回路６に残留する絶対的な冷媒量を推定する場合は、ここまでに説明した各運転状態量に加えて、室外熱交換器１３および各室内熱交換器１の容積や液管４の容積を考慮すればよい。

＜変形例＞
　尚、本実施例では、例えば、第１の暖房用推定モデル７３Ｄの推定結果と第２の暖房用推定モデル７３Ｅの推定結果との間をシグモイド係数で補間する場合を例示したが、シグモイド係数に限定されるものではなく、例えば、線形補間等の補間方法を使用しても良く、適宜変更可能である。

　なお、本実施例では、推定モデルは事前に生成されていたものを利用していた。しかし、サーバ１２０で生成するようにしても良い。例えば、冷媒回路６に充填される冷媒の冷媒不足率の推定に関わる運転状態量と、冷媒量を計測する計測器からの計測結果とを用いて重回帰分析法を使用して、サーバ１２０が冷媒不足率を推定する推定モデルを生成してもよい。また、本実施例では、重回帰分析法を用いて各推定モデルを生成する場合を例示したが、一般の回帰分析法を行える機械学習手法のSVR(Support　Vector　Regression)、NN(Neural　Network)などを用いて推定モデルを生成しても良い。その際、特徴量選択に当たっては重回帰分析法で用いたＰ値や補正値Ｒ２の代わりに、推定モデルの精度が向上するよう特徴量を選択する一般の手法（Forward　Feature　Selection法、Backward　feature　Eliminationなど）を使えばよい。

　尚、センサ値編集処理を実行する制御部７４は、代表時刻付近のセンサ値が複数ある場合に「センサ値変化」があるか否かを判定し、「センサ値変化」がある場合、最も早い時刻の「センサ値変化」のセンサ値を代表センサ値と決定する場合を例示した。しかしながら、最も早い時刻の「センサ値変化」に限定されるものではなく、例えば、「センサ値変化」のセンサ値の平均値や、最も遅い時刻の「センサ値」としても良く、適宜変更可能である。

　制御部７４は、代表時刻付近の検出時刻の各センサ値が取得できていない場合、直前の代表時刻の代表センサ値を、代表時刻の代表センサ値と決定する場合を例示した。しかしながら、制御部７４は、直前の代表時刻の代表センサ値に限定されるものではなく、例えば、直前の「センサ値変化」のセンサ値を使用しても良く、適宜変更可能である。

　また、図示した各部の各構成要素は、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各部の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況等に応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。

　更に、各装置で行われる各種処理機能は、ＣＰＵ（Central　Processing　Unit）（又はＭＰＵ（Micro　Processing　Unit）、ＭＣＵ（Micro　Controller　Unit）等のマイクロ・コンピュータ）上で、その全部又は任意の一部を実行するようにしても良い。また、各種処理機能は、ＣＰＵ（又はＭＰＵ、ＭＣＵ等のマイクロ・コンピュータ）で解析実行するプログラム上、又はワイヤードロジックによるハードウェア上で、その全部又は任意の一部を実行するようにしても良いことは言うまでもない。

　また、以上に説明した各実施例では、冷媒回路６に残存する冷媒量を示す指標として、例えば相対的な冷媒量である冷媒不足率を用いた。冷媒不足率は、冷媒回路６に規定量の冷媒が充填されている状態を冷媒充填率１００％としたとき、この規定量からの冷媒の減少率を意味する。しかし、冷媒回路６に残存する冷媒量を示す指標として、冷媒の減少率に代えて冷媒の充填率を用いても良い。また、冷媒の不足率や充填率を表す際の基準量（規定量）を予め定められた冷媒量としていたが、これに代えて冷媒回路に実際に充填した冷媒量を基準量（規定量）とみなしても良い。この場合は例えば、冷媒回路６に実際に充填された冷媒量が、予め定められる規定量より少ない（又は多い）場合でも、この冷媒量を１００％とすることができる。このように実際に充填された冷媒量を基準量とすることで、冷媒回路毎の冷媒不足率をより正確に推定できる。さらに、冷媒回路６に残存する冷媒量を示す指標として、相対的な指標（割合）に代えて絶対的な指標である冷媒量を用いてもよい。

　１　空気調和機
　２　室外機
　３　室内機
　４　液管
　５　ガス管
　１１　圧縮機
　１２　四方弁
　１３　室外熱交換器
　１４　室外機膨張弁
　１９　室外機制御部
　１９Ａ　室外側検出部
　１９Ｂ　室外側記憶部
　１９Ｃ　室外側制御部
　５１　室内熱交換器
　６５　室内機制御部
　６５Ａ　室内側検出部
　６５Ｂ　室内側記憶部
　６５Ｃ　室内側制御部
　７１　取得部
　７３Ｄ　第１の暖房用推定モデル
　７３Ｅ　第２の暖房用推定モデル
　７３Ｆ　第３の暖房用推定モデル
　７４　制御部
　７４Ａ　推定部

Claims

　圧縮機、室外熱交換器及び膨張弁を有する室外機と、室内熱交換器を有する室内機とを有し、前記室外機と前記室内機とが冷媒配管で接続されて形成される冷媒回路を有し、前記室内熱交換器を前記圧縮機において圧縮される冷媒の凝縮器として、かつ、前記室外熱交換器を前記室内熱交換器において凝縮される冷媒の蒸発器として機能させる暖房運転を少なくとも行うことが可能な空気調和機であって、
　前記空気調和機は、
　少なくとも前記暖房運転における前記空気調和機の運転状態量を用いて前記冷媒回路に残存する冷媒量を推定する推定部を有し、
　前記推定部は、
　前記冷媒回路に残存する冷媒量の範囲に対応させた異なる複数の推定モデルを含み、前記複数の推定モデルのうち一つは、前記運転状態量として前記室内側熱交換器の出口における冷媒の過冷却度を用いることを特徴とする空気調和機。
　前記複数の推定モデルのうち、前記冷媒回路に残存する冷媒量が多い範囲に対応させた推定モデルを第１の推定モデルとし、前記冷媒回路に残存する冷媒量が少ない範囲に対応させた推定モデルを第２の推定モデルとしたとき、前記第１の推定モデルは前記運転状態量として前記室内側熱交換器の出口における冷媒の過冷却度を用いることを特徴とする、請求項１に記載の空気調和機。
　前記第１の推定モデルと前記第２の推定モデルとで構成される推定モデルを第３の推定モデルとしたとき、前記推定部は、前記第３の推定モデルを含む、ことを特徴とする、請求項２に記載の空気調和機。
　前記室内機は、複数台設置されており、
　前記推定部は、
　前記室内機のうち、少なくとも２台以上の室内機の室内熱交換器を、前記圧縮機において圧縮される冷媒の凝縮器として機能させる場合に、前記凝縮器として機能する室内熱交換器の出口における冷媒の過冷却度を用いて、前記推定モデルで冷媒量を推定することを特徴とする請求項１～３の何れか一つに記載の空気調和機。
　前記推定部は、
　前記２台以上の室内機のそれぞれの室内熱交換器の出口における冷媒の温度の平均値に基づく過冷却度を用いて、前記冷媒量を推定することを特徴とする請求項４に記載の空気調和機。
　前記室内機には、
　前記室内機の各部の動作を制御する室内側制御部と、前記運転状態量のうち室内機側の運転状態量である室内側運転状態量を検出する室内側検出部と、前記室内側検出部で検出した室内側検出結果を記憶する室内側記憶部と、を備え、
　前記室外機には、
　前記室外機の各部の動作を制御する室外側制御部と、前記運転状態量のうち室外機側の運転状態量である室外側運転状態量を検出する室外側検出部と、前記室外側検出部で検出した室外側検出結果を記憶する室外側記憶部と、を備え、
　前記室内側制御部は、
　前記室内側検出結果を、検出時刻と紐づけて前記室内側記憶部に格納し、
　前記室外側制御部は、
　前記室外側検出結果を、検出時刻と紐づけて前記室内側記憶部に格納する、
ことを特徴とする請求項１～５の何れか一つに記載の空気調和機。
　前記室内機及び前記室外機の状態を表示する集中制御手段を備え、
　前記集中制御手段は、
　制御部と、記憶部とを備え、
　前記記憶部は、
　検出時刻と紐づけられた前記室内側検出結果と、検出時刻と紐づけられた前記室外側検出結果を記憶し、
　前記制御部は、
　前記室内側検出結果の検出時刻と前記室外側検出結果の検出時刻とが所定範囲内にあるとき、前記室内側検出結果と前記室外側検出結果を新たな時刻とを紐づけて前記記憶部に格納することを特徴とする請求項６に記載の空気調和機。
　前記室内側検出部は、
　前記室内熱交換器の出口における冷媒の温度を前記室内側検出結果として検出する第１のセンサを含み、
　前記室内側検出部は、
　前記室外熱交換器の高圧飽和温度を前記室外側検出結果として検出する第２のセンサを含み、
　前記推定部は、
　前記室内側検出結果の検出時刻と前記室外側検出結果の検出時刻とが所定範囲内にある前記室内側検出結果及び前記室外側検出結果を用いて算出した前記過冷却度を用いて前記冷媒量を推定することを特徴とする請求項６～７の何れか一つに記載の空気調和機。
　前記推定部は、
　前記室内機のうち、少なくとも２台以上の室内機の室内熱交換器を、前記圧縮機において圧縮される前記冷媒の凝縮器として機能させる場合に、前記２台以上の室内機のそれぞれの前記第１のセンサが検出した検出結果の平均値に基づく過冷却度を用いて前記冷媒量を推定することを特徴とする請求項８に記載の空気調和機。
　前記推定モデルは、
　前記冷媒回路に残存する冷媒量として、前記冷媒回路から減少した冷媒の割合を示す冷媒不足率を推定する推定モデルであることを特徴とする請求項１～９の何れか一つに記載の空気調和機。
　圧縮機、室外熱交換器及び膨張弁を有する室外機と、室内熱交換器を有する室内機とを有し、前記室外機と前記室内機とが冷媒配管で接続されて形成される冷媒回路を有し、前記室内熱交換器を前記圧縮機において圧縮される冷媒の凝縮器として、かつ、前記室外熱交換器を前記室内熱交換器において凝縮される冷媒の蒸発器として機能させる暖房運転を少なくとも行うことが可能な空気調和機と、前記空気調和機と通信で接続するサーバとを有する空気調和システムであって、
　前記サーバは、
　少なくとも前記暖房運転における前記空気調和機の運転状態量を用いて前記冷媒回路に残存する冷媒量を推定する推定部を有し、
　前記推定部は、
　前記冷媒回路に残存する冷媒量の範囲に対応させた異なる複数の推定モデルを含み、前記複数の推定モデルのうち一つは、前記暖房運転における前記空気調和機の運転状態量として前記室内側熱交換器の出口における冷媒の過冷却度を用いることを特徴とする空気調和システム。
　前記複数の推定モデルのうち、前記冷媒回路に残存する冷媒量が多い範囲に対応させた推定モデルを第１の推定モデルとし、前記冷媒回路に残存する冷媒量が少ない範囲に対応させた推定モデルを第２の推定モデルとしたとき、前記第１の推定モデルは前記運転状態量として前記室内側熱交換器の出口における冷媒の過冷却度を用いることを特徴とする、請求項１１に記載の空気調和システム。
　前記第１の推定モデルと前記第２の推定モデルとで構成される推定モデルを第３の推定モデルとしたとき、前記推定部は、前記第３の推定モデルを含む、ことを特徴とする、請求項１２に記載の空気調和システム。
　前記室内機は、複数台設置されており、
　前記推定部は、
　前記室内機のうち、少なくとも２台以上の室内機の室内熱交換器を、前記圧縮機において圧縮される冷媒の凝縮器として機能させる場合に、前記凝縮器として機能する室内熱交換器の出口における冷媒の過冷却度を用いて、前記推定モデルで冷媒量を推定することを特徴とする請求項１１～１３の何れか一つに記載の空気調和システム。
　前記室内機及び前記室外機の状態を表示する集中制御手段を備え、前記空気調和機と前記サーバとが、前記集中制御手段を介して通信で接続される請求項１１～１４の何れか一つに記載の空気調和システム。
　前記推定部は、
　前記２台以上の室内機のそれぞれの室内熱交換器の出口における冷媒の温度の平均値に基づく過冷却度を用いて、前記冷媒量を推定することを特徴とする請求項１４に記載の空気調和システム。
　前記室内機には、
　前記室内機の各部の動作を制御する室内側制御部と、前記運転状態量のうち室内機側の運転状態量である室内側運転状態量を検出する室内側検出部と、前記室内側検出部で検出した室内側検出結果を記憶する室内側記憶部と、を備え、
　前記室外機には、
　前記室外機の各部の動作を制御する室外側制御部と、前記運転状態量のうち室外機側の運転状態量である室外側運転状態量を検出する室外側検出部と、前記室外側検出部で検出した室外側検出結果を記憶する室外側記憶部と、を備え、
　前記室内側制御部は、
　前記室内側検出結果を、検出時刻と紐づけて前記室内側記憶部に格納し、
　前記室外側制御部は、
　前記室外側検出結果を、検出時刻と紐づけて前記室内側記憶部に格納する、
ことを特徴とする請求項１１～１６の何れか一つに記載の空気調和システム。
　前記集中制御手段は
　制御部と、記憶部とを備え、
　前記記憶部は、
　検出時刻と紐づけられた前記室内側検出結果と、検出時刻と紐づけられた前記室外側検出結果とを記憶し、
　前記制御部は、
　前記室内側検出結果の検出時刻と前記室外側検出結果の検出時刻とが所定範囲内にあるとき、前記室内側検出結果と前記室外側検出結果とを新たな時刻と紐づけて前記記憶部に格納することを特徴とする請求項１７に記載の空気調和システム。
　前記室内側検出部は、
　前記室内熱交換器の出口における冷媒の温度を前記室内側検出結果として検出する第１のセンサを含み、
　前記室内側検出部は、
　前記室外熱交換器の高圧飽和温度を前記室外側検出結果として検出する第２のセンサを含み、
　前記推定部は、
　前記室内側検出結果の検出時刻と前記室外側検出結果の検出時刻とが所定範囲内にある前記室内側検出結果及び前記室外側検出結果を用いて算出した前記過冷却度を用いて前記冷媒量を推定することを特徴とする請求項１７又は１８に記載の空気調和システム。
　前記推定部は、
　前記室内機のうち、少なくとも２台以上の室内機の室内熱交換器を、前記圧縮機において圧縮される前記冷媒の凝縮器として機能させる場合に、前記２台以上の室内機のそれぞれの前記第１のセンサが検出した検出結果の平均値に基づく過冷却度を用いて前記冷媒量を推定することを特徴とする請求項１９に記載の空気調和システム。
　前記推定モデルは、
　前記冷媒回路に残存する冷媒量として前記冷媒回路から漏洩した冷媒の割合を示す冷媒不足率を推定する推定モデルであることを特徴とする請求項１１～２０の何れか一つに記載の空気調和システム。