ホーム > 特許ランキング > 株式会社富士通ゼネラル
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-02-13
(45)【発行日】2024-02-21
(54)【発明の名称】空気調和機
(51)【国際特許分類】
F25B 49/02 20060101AFI20240214BHJP
F24F 11/38 20180101ALI20240214BHJP
【ＦＩ】
F25B49/02 520D
F24F11/38
F25B49/02 520H
F25B49/02 520Z
【請求項の数】
12
(21)【出願番号】P 2020058823
(22)【出願日】2020-03-27
(65)【公開番号】P2021156532
(43)【公開日】2021-10-07
【審査請求日】2022-09-29
(73)【特許権者】
【識別番号】000006611
【氏名又は名称】株式会社富士通ゼネラル
(74)【代理人】
【識別番号】110002147
【氏名又は名称】弁理士法人酒井国際特許事務所
(72)【発明者】
【氏名】佐々木 寛
(72)【発明者】
【氏名】奥村 聡司
(72)【発明者】
【氏名】上田 晴康
(72)【発明者】
【氏名】白石 譲熙
(72)【発明者】
【氏名】下谷 亮
【審査官】笹木 俊男
(56)【参考文献】
【文献】特開２０１８－１９４１７８（ＪＰ，Ａ）
【文献】米国特許出願公開第２０１８／０２１６８５９（ＵＳ，Ａ１）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｆ２５Ｂ ４９／０２
Ｆ２４Ｆ １１／００ ～ １１／８９
(57)【特許請求の範囲】
【請求項1】
圧縮機、室外熱交換器及び膨張弁を有する室外機と、室内熱交換器を有する室内機とを有し、前記室外機と前記室内機とが冷媒配管で接続されて形成される冷媒回路を有し、当該冷媒回路に所定量の冷媒が充填される空気調和機であって、前記空気調和機の現在の運転状態量を用いて、前記冷媒回路に残存する冷媒量を推定する、前記空気調和機の運転モード毎に備えた推定モデルは、前記残存する冷媒量の範囲に応じて異なり、
前記推定モデルは、
前記空気調和機内の冷房運転時の前記残存する冷媒量が第１の範囲にある場合に使用する第１の冷房用推定モデルと、
前記冷房運転時の前記残存する冷媒量が前記第１の範囲を超え、かつ、当該第１の範囲とは異なる第２の範囲にある場合に使用する第２の冷房用推定モデルとを有することを特徴とする空気調和機。
【請求項2】
前記推定モデルは、前記空気調和機の現在の運転状態を示す運転状態量の内、少なくとも前記圧縮機の回転数及び外気温度の何れか一つを用いて、前記残存する冷媒量を推定することを特徴とする請求項１に記載の空気調和機。
【請求項3】
前記第１の冷房用推定モデルは、前記運転状態量として、前記室外熱交換器又は前記室内熱交換器の内、凝縮器として機能する当該凝縮器の冷媒過冷却度を用いて、前記残存する冷媒量を推定することを特徴とする請求項１に記載の空気調和機。
【請求項4】
前記第２の冷房用推定モデルは、前記運転状態量として、前記圧縮機の吸入温度を用いて、前記残存する冷媒量を推定することを特徴とする請求項１に記載の空気調和機。
【請求項5】
圧縮機、室外熱交換器及び膨張弁を有する室外機と、室内熱交換器を有する室内機とを有し、前記室外機と前記室内機とが冷媒配管で接続されて形成される冷媒回路を有し、当該冷媒回路に所定量の冷媒が充填される空気調和機であって、
前記空気調和機の現在の運転状態量を用いて、前記冷媒回路に残存する冷媒量を推定する推定モデルを有し、前記推定モデルは前記残存する冷媒量に応じて異なり、
前記推定モデルは、
前記空気調和機内の暖房運転時の前記残存する冷媒量が第３の範囲内の場合に使用する第１の暖房用推定モデルと、
前記暖房運転時の前記残存する冷媒量が前記第３の範囲を超え、かつ、第４の範囲内の場合に使用する第２の暖房用推定モデルと
を有することを特徴とする空気調和機。
【請求項6】
前記第１の暖房用推定モデルは、前記運転状態量として前記膨張弁の開度を用いて、前記冷媒回路に残存する冷媒量を推定することを特徴とする請求項５に記載の空気調和機。
【請求項7】
前記第２の暖房用推定モデルは、前記運転状態量として前記圧縮機の吸入過熱度を用いて、前記残存する冷媒量を推定することを特徴とする請求項５に記載の空気調和機。
【請求項8】
前記推定モデルは、前記残存する冷媒量の推定に使用する前記運転状態量毎に所定の重みパラメータを含むことを特徴とする請求項１～７の何れか一つに記載の空気調和機。
【請求項9】
前記推定モデルは、前記残存する冷媒量の推定に与える影響の大きい前記運転状態量の変化量が、前記空気調和機の運転状態が安定している状態である所定範囲内の運転状態量を用いて前記残存する冷媒量を推定することを特徴とする請求項１～８の何れか一つに記載の空気調和機。
【請求項10】
前記空気調和機の現在の運転状態量を用いて、前記残存する冷媒量として、前記所定量から減少した冷媒の割合である冷媒不足率を推定する推定モデルを有し、前記推定モデルは、前記冷媒不足率に応じて異なることを特徴とする請求項１に記載の空気調和機。
【請求項11】
前記推定モデルは、前記空気調和機内の冷房運転時の冷媒不足率が０％～３０％の範囲内の場合に使用する、第１の冷房用推定モデルと、
前記冷房運転時の冷媒不足率が前記０％～３０％の範囲を超え、かつ、４０％～７０％の範囲内の場合に使用する、第２の冷房用推定モデルと
を有することを特徴とする請求項１０に記載の空気調和機。
【請求項12】
前記推定モデルは、前記空気調和機内の暖房運転時の冷媒不足率が０％～２０％の範囲内の場合に使用する、第１の暖房用推定モデルと、
前記暖房運転時の冷媒不足率が前記０％～２０％の範囲を超え、かつ、３０％～７０％の範囲内の場合に使用する、前記冷媒不足率を推定する第２の暖房用推定モデルと
を有することを特徴とする請求項１０に記載の空気調和機。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、空気調和機に関する。
【背景技術】
【0002】
冷媒回路で検知できる運転状態量を用いて冷媒量の判定を行う空気調和機が提案されている。特許文献１では、例えば、冷房サイクル時の冷媒回路の蒸発器出口の過熱度や蒸発器の圧力を所定の値にした状態（以下、デフォルト状態という）で凝縮器出口の過冷却度を用いて冷媒量を判定している。
【先行技術文献】
【特許文献】
【0003】
【文献】特開２００６－２３０７２号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
空気調和機では、過冷却度等の運転状態量を用いて冷媒量を判定する場合、冷媒回路を上述したデフォルト状態、例えば、冷媒を規定量に充填した直後に運転状態量を用いて冷媒量を求めたときと同じ状態にする必要がある。
【0005】
しかしながら、冷媒量を判定する時点の外気温や室内温度などの外部環境が季節の違いや日射量等による影響で、冷媒を充填した直後の外部環境に一致するとは限らない。このため、冷媒量を判定したいときに冷媒回路の状態をデフォルト状態に合わせることが難しく、所望のタイミングでデフォルト状態に合わせることができずに冷媒量を判定できない場合があった。また、特許文献１に記載の方法では、冷媒量の判定を行うために冷媒回路の状態をデフォルト状態にする必要があるために、これまでに行っていた冷房運転や暖房運転が中断されるので、使用者に不快感を与える恐れがあった。
【0006】
本発明ではこのような問題に鑑み、所望のタイミングで、かつ、現在行っている空調運転を中断せずに冷媒量を判定できる空気調和機を提供することを目的とする。
【課題を解決するための手段】
【0007】
一つの態様の空気調和機は、圧縮機、室外熱交換器及び膨張弁を有する室外機と、室内熱交換器を有する室内機とを有し、前記室外機と前記室内機とが冷媒配管で接続されて形成される冷媒回路を有し、当該冷媒回路に所定量の冷媒が充填される空気調和機であって、前記空気調和機の現在の運転状態量を用いて、前記冷媒回路に残存する冷媒量を推定する推定モデルを有し、前記推定モデルは前記残存する冷媒量に応じて異なる。
【発明の効果】
【0008】
一つの側面として、所望のタイミングで、かつ、現在行っている空調運転を中断せずに冷媒量を判定できる。
【図面の簡単な説明】
【0009】
【発明を実施するための形態】
【0010】
以下、図面に基づいて、本願の開示する空気調和機等の実施例を詳細に説明する。尚、本実施例により、開示技術が限定されるものではない。また、以下に示す各実施例は、矛盾を起こさない範囲で適宜変形しても良い。
【実施例1】
【0011】
＜空気調和機の構成＞
図１は、本実施例の空気調和機１の一例を示す説明図である。図１に示す空気調和機１は、１台の室外機２と、Ｎ台の室内機３とを有する（Ｎは２以上の自然数）。室外機２は、液管４及びガス管５で並列に各室内機３と接続する。そして、室外機２と室内機３とが液管４及びガス管５等の冷媒配管で接続することで、空気調和機１の冷媒回路６が形成されている。
図１は、本実施例の空気調和機１の一例を示す説明図である。図１に示す空気調和機１は、１台の室外機２と、Ｎ台の室内機３とを有する（Ｎは２以上の自然数）。室外機２は、液管４及びガス管５で並列に各室内機３と接続する。そして、室外機２と室内機３とが液管４及びガス管５等の冷媒配管で接続することで、空気調和機１の冷媒回路６が形成されている。
【0012】
＜室外機の構成＞
図２は、室外機２およびＮ台の室内機３の一例を示す説明図である。室外機２は、圧縮機１１と、四方弁１２と、室外熱交換器１３と、室外機膨張弁１４と、第１の閉鎖弁１５と、第２の閉鎖弁１６と、アキュムレータ１７と、室外機ファン１８と、制御手段１９とを有する。これら圧縮機１１、四方弁１２、室外熱交換器１３、室外機膨張弁１４、第１の閉鎖弁１５、第２の閉鎖弁１６及びアキュムレータ１７を用いて、以下で詳述する各冷媒配管で相互に接続されて冷媒回路６の一部を成す室外側冷媒回路を形成する。
図２は、室外機２およびＮ台の室内機３の一例を示す説明図である。室外機２は、圧縮機１１と、四方弁１２と、室外熱交換器１３と、室外機膨張弁１４と、第１の閉鎖弁１５と、第２の閉鎖弁１６と、アキュムレータ１７と、室外機ファン１８と、制御手段１９とを有する。これら圧縮機１１、四方弁１２、室外熱交換器１３、室外機膨張弁１４、第１の閉鎖弁１５、第２の閉鎖弁１６及びアキュムレータ１７を用いて、以下で詳述する各冷媒配管で相互に接続されて冷媒回路６の一部を成す室外側冷媒回路を形成する。
【0013】
圧縮機１１は、例えば、インバータにより回転数が制御される図示しないモータの駆動に応じて、運転容量を可変できる高圧容器型の能力可変型圧縮機である。圧縮機１１は、その冷媒吐出側と四方弁１２の第１のポート１２Ａとの間を吐出管２１で接続している。また、圧縮機１１は、その冷媒吸入側とアキュムレータ１７の冷媒流出側との間を吸入管２２で接続している。
【0014】
四方弁１２は、冷媒回路６における冷媒の流れる方向を切替えるための弁であって、第１～第４のポート１２Ａ～１２Ｄを備えている。第１のポート１２Ａは、圧縮機１１の冷媒吐出側との間を吐出管２１で接続している。第２のポート１２Ｂは、室外熱交換器１３の一方の冷媒出入口との間を室外冷媒管２３で接続している。第３のポート１２Ｃは、アキュムレータ１７の冷媒流入側との間を室外冷媒管２６で接続している。そして、第４のポート１２Ｄは、第２の閉鎖弁１６との間を室外ガス管２４で接続している。
【0015】
室外熱交換器１３は、冷媒と、室外機ファン１８の回転により室外機２の内部に取り込まれた外気とを熱交換させる。室外熱交換器１３は、その一方の冷媒出入口と四方弁１２の第２のポート１２Ｂとの間を室外冷媒管２６で接続している。室外熱交換器１３は、その他方の冷媒出入口と第１の閉鎖弁１５との間を室外液管２５で接続している。室外熱交換器１３は、空気調和機１が冷房運転を行う場合に凝縮器として機能し、空気調和機１が暖房運転を行う場合に蒸発器として機能する。
【0016】
室外機膨張弁１４は、室外液管２５に設けられており、図示しないパルスモータで駆動する電子膨張弁である。室外機膨張弁１４は、パルスモータに与えられるパルス数に応じて開度が調整されることで、室外熱交換器１３に流入する冷媒量、又は、室外熱交換器１３から流出する冷媒量を調整するものである。室外機膨張弁１４の開度は、空気調和機１が暖房運転を行っている場合、圧縮機１１の冷媒吸入側の冷媒過熱度が目標吸入冷媒過熱度となるように調整される。また、室外機膨張弁１４の開度は、空気調和機１が冷房運転を行っている場合、全開とされる。
【0017】
アキュムレータ１７は、その冷媒流入側と四方弁１２の第３のポート１２Ｃとの間を室外冷媒管２６で接続している。更に、アキュムレータ１７は、その冷媒流出側と圧縮機１１の冷媒流入側との間を吸入管２２で接続している。アキュムレータ１７は、室外冷媒管２６からアキュムレータ１７の内部に流入した冷媒をガス冷媒と液冷媒とに分離し、ガス冷媒のみを圧縮機１１に吸入させる。
【0018】
室外機ファン１８は、樹脂材で形成されており、室外熱交換器１３の近傍に配置されている。室外機ファン１８は、図示しないファンモータの回転に応じて、図示しない吸込口から室外機２の内部へ外気を取り込み、室外熱交換器１３において冷媒と熱交換した外気を図示しない吹出口から室外機２の外部へ放出する。
【0019】
また、室外機２には、複数のセンサが配置されている。吐出管２１には、圧縮機１１から吐出される冷媒の圧力、すなわち吐出圧力を検出する吐出圧センサ３１と、圧縮機１１から吐出された冷媒の温度、すなわち吐出温度を検出する吐出温度センサ３２とが配置されている。室外冷媒管２６のアキュムレータ１７の冷媒流入口近傍には、圧縮機１１に吸入される冷媒の圧力である吸入圧力を検出する吸入圧力センサ３３と、圧縮機１１に吸入される冷媒の温度を検出する吸入温度センサ３４とが配置されている。
【0020】
室外熱交換器１３と室外機膨張弁１４との間の室外液管２５には、室外熱交換器１３に流入する冷媒の温度、又は、室外熱交換器１３から流出する冷媒の温度を検出するための冷媒温度センサ３５が配置されている。そして、室外機２の図示しない吸込口付近には、室外機２の内部に流入する外気の温度、すなわち外気温度を検出する外気温度センサ３６が配置されている。
【0021】
制御手段１９は、空気調和機１全体を制御する。図３は、室外機２の制御手段１９の一例を示すブロック図である。制御手段１９は、取得部４１と、通信部４２と、記憶部４３と、制御部４４とを有する。取得部４１は、前述した各種センサのセンサ値を取得する。通信部４２は、各室内機３の通信部と通信する通信インタフェースである。記憶部４３は、例えば、フラッシュメモリであって、室外機２の制御プログラムや各種センサからの検出信号に対応した検出値等の運転状態量、圧縮機１１や室外機ファン１８の駆動状態、各室内機３から送信される運転情報（例えば、運転・停止情報、冷房／暖房等の運転モード等を含む）、室外機２の定格能力及び各室内機３の要求能力、などを記憶する。
【0022】
また、記憶部４３は冷媒回路６に残存する冷媒量を推定する推定モデルを記憶している。本実施例では、冷媒回路６に残存する冷媒量として、例えば相対的な冷媒量を用いている。具体的には、本実施例の記憶部４３は冷媒回路６の冷媒不足率（冷媒が規定量充填されているときを１００％としたとき、この規定量からの減少分を指す。以下、同様）を推定する推定モデルを記憶している。推定モデルは、第１の冷房用推定モデル４３Ａと、第２の冷房用推定モデル４３Ｂと、第３の冷房用推定モデル４３Ｃと、第１の暖房用推定モデル４３Ｄと、第２の暖房用推定モデル４３Ｅと、第３の暖房用推定モデル４３Ｆとを有する。
【0023】
制御部４４は、通信部４２を介して各種センサでの検出値を定期的（例えば、３０秒毎）に取り込み、各室内機３から送信される運転情報を含む信号が通信部４２を介して入力される。制御部４４は、これら入力された各種情報に基づいて、室外機膨張弁１４の開度調整や圧縮機１１の駆動制御を行う。更に、制御部４４は、上述した各推定モデルを用いて冷媒不足率を推定する。
【0024】
＜室内機の構成＞
図２に示すように、室内機３は、室内熱交換器５１と、室内機膨張弁５２と、液管接続部５３と、ガス管接続部５４と、室内機ファン５５とを有する。これら室内熱交換器５１、室内機膨張弁５２、液管接続部５３及びガス管接続部５４は、後述する各冷媒配管で相互に接続されて、冷媒回路６の一部を成す室内機冷媒回路を構成する。
図２に示すように、室内機３は、室内熱交換器５１と、室内機膨張弁５２と、液管接続部５３と、ガス管接続部５４と、室内機ファン５５とを有する。これら室内熱交換器５１、室内機膨張弁５２、液管接続部５３及びガス管接続部５４は、後述する各冷媒配管で相互に接続されて、冷媒回路６の一部を成す室内機冷媒回路を構成する。
【0025】
室内熱交換器５１は、冷媒と、室内機ファン５５の回転により図示しない吸込口から室内機３の内部に取り込まれた室内空気とを熱交換させる。室内熱交換器５１は、その一方の冷媒出入口と液管接続部５３との間を室内液管５６で接続している。また、室内熱交換器５１は、その他方の冷媒出入口とガス管接続部５４との間を室内ガス管５７で接続している。室内熱交換器５１は、空気調和機１が暖房運転を行う場合、凝縮器として機能する。これに対して、室内熱交換器５１は、空気調和機１が冷房運転を行う場合、蒸発器として機能する。
【0026】
室内機膨張弁５２は、室内液管５６に設けられており、電子膨張弁である。室内熱交換器５１が蒸発器として機能する場合、すなわち、室内機３が冷房運転を行う場合、室内機膨張弁５２の開度は、室内熱交換器５１の冷媒出口（ガス管接続部５４側）での冷媒過熱度が目標冷媒過熱度となるように調整される。また、室内熱交換器５１が凝縮器として機能する場合、すなわち室内機３が暖房運転を行う場合、室内機膨張弁５２の開度は、室内熱交換器５１の冷媒出口（液管接続部５３側）での冷媒過冷却度が目標冷媒過冷却度となるように調整される。ここで、目標冷媒過熱度や目標冷媒過冷却度とは、室内機３で十分な冷房能力あるいは暖房能力を発揮するのに必要な冷媒過熱度および冷媒過冷却度である。
【0027】
室内機ファン５５は、樹脂材で形成されており、室内熱交換器５１の近傍に配置されている。室内機ファン５５は、図示しないファンモータによって回転することで、図示しない吸込口から室内機３の内部に室内空気を取り込み、室内熱交換器５１において冷媒と熱交換した室内空気を図示しない吹出口から室内へ放出する。
【0028】
室内機３には各種のセンサが設けられている。室内液管５６には、室内熱交換器５１と室内機膨張弁５２との間に、室内熱交換器５１に流入する冷媒の温度、又は室内熱交換器５１から流出する冷媒の温度を検出する液側冷媒温度センサ６１が配置されている。室内ガス管５７には、室内熱交換器５１から流出又は室内熱交換器５１に流入する冷媒の温度を検出するガス側温度センサ６２が配置されている。室内機３の図示しない吸込口付近には、室内機３の内部に流入する室内空気の温度、すなわち吸込温度を検出する吸込温度センサ６３が配置されている。
【0029】
＜冷媒回路の動作＞
次に、本実施形態における空気調和機１の空調運転時の冷媒回路６における冷媒の流れや各部の動作について説明する。尚、図１における矢印は暖房運転時の冷媒の流れを示している。
次に、本実施形態における空気調和機１の空調運転時の冷媒回路６における冷媒の流れや各部の動作について説明する。尚、図１における矢印は暖房運転時の冷媒の流れを示している。
【0030】
空気調和機１が暖房運転を行う場合、四方弁１２は、第１のポート１２Ａと第４のポート１２Ｄとが連通し、第２のポート１２Ｂと第３のポート１２Ｃとが連通するように切替えている。これにより、冷媒回路６は、各室内熱交換器５１が凝縮器として機能し、室外熱交換器１３が蒸発器として機能する暖房サイクルとなる。尚、説明の便宜上、暖房運転時の冷媒の流れは、図２に示す実線矢印で表記する。
【0031】
冷媒回路６が上記の状態で圧縮機１１が駆動すると、圧縮機１１から吐出された冷媒は、吐出管２１を流れて四方弁１２に流入し、四方弁１２から室外ガス管２４を流れて、第２の閉鎖弁１６を介してガス管５へと流入する。ガス管５を流れる冷媒は、各ガス管接続部５４を介して各室内機３に分流する。各室内機３に流入した冷媒は、各室内ガス管５７を流れて各室内熱交換器５１に流入する。各室内熱交換器５１に流入した冷媒は、各室内機ファン５５の回転により各室内機３の内部に取り込まれた室内空気との間で熱交換することで凝縮する。つまり、各室内熱交換器５１が凝縮器として機能し、各室内熱交換器５１で冷媒によって加熱された室内空気が図示しない吹出口から室内に吹き出されることで、各室内機３が設置された室内の暖房が行われる。
【0032】
各室内熱交換器５１から各室内液管５６に流入した冷媒は、各室内熱交換器５１の冷媒出口側での冷媒過冷却度が目標冷媒過冷却度となるように開度が調整された各室内機膨張弁５２を通過して減圧される。ここで、目標冷媒過冷却度は、各室内機３で要求される冷房能力に基づいて定められるものである。
【0033】
各室内機膨張弁５２で減圧された冷媒は、各室内液管５６から各液管接続部５３を介して液管４に流出する。液管４で合流した冷媒は、第１の閉鎖弁１５を介して室外機２に流入する。室外機２の第１の閉鎖弁１５に流入した冷媒は、室外液管２５を流れ、室外機膨張弁１４を通過して減圧される。室外機膨張弁１４で減圧された冷媒は、室外液管２５を流れて室外熱交換器１３に流入し、室外機ファン１８の回転によって室外機２の図示しない吸込口から流入した外気と熱交換を行って蒸発する。室外熱交換器１３から室外冷媒管２６へと流出した冷媒は、四方弁１２、室外冷媒管２６、アキュムレータ１７及び吸入管２２の順に流入し、圧縮機１１に吸入されて再び圧縮され、四方弁１２の第１のポート１２Ａ及び第４のポート１２Ｄ経由で室外ガス管２４に流出する。
【0034】
また、空気調和機１が冷房運転を行う場合、四方弁１２は、第１のポート１２Ａと第２のポート１２Ｂとが連通し、第３のポート１２Ｃと第４のポート１２Ｄとが連通するように切替えている。これにより、冷媒回路６は、各室内熱交換器５１が蒸発器として機能し、室外熱交換器１３が凝縮器として機能する冷房サイクルとなる。尚、説明の便宜上、冷房運転時の冷媒の流れは、図２に示す破線矢印で表記する。
【0035】
冷媒回路６の状態で圧縮機１１が駆動すると、圧縮機１１から吐出された冷媒は、吐出管２１を流れて四方弁１２に流入し、四方弁１２から室外冷媒管２６を流れて、室外熱交換器１３に流入する。室外熱交換器１３に流入した冷媒は、室外機ファン１８の回転により室外機２の内部に取り込まれた室外空気との間で熱交換することで凝縮する。つまり、室外熱交換器１３が凝縮器として機能し、室外熱交換器１３で冷媒によって加熱された室内空気が図示しない吹出口から室外に吹き出す。
【0036】
室外熱交換器１３から室外液管２５へと流入した冷媒は、開度が全開とされている室外機膨張弁１４を通過して減圧される。室外機膨張弁１４で減圧された冷媒は、第１の閉鎖弁１５を介して液管４を流れて各室内機３に分流する。各室内機３に流入した冷媒は、各液管接続部５３を通じて室内液管５６を流れて室内熱交換器５１の冷媒出口で冷媒過冷却度が目標冷媒過冷却度となる開度に調整された室内機膨張弁５２を通過して減圧される。室内機膨張弁５２で減圧された冷媒は、室内液管５６を流れて室内熱交換器５１に流入し、室内機ファン５５の回転によって室内機３の図示しない吸入口から流入した室内空気と熱交換を行って蒸発する。つまり、各室内熱交換器５１が蒸発器として機能し、各室内熱交換器５１で冷媒によって冷却された室内空気が図示しない吹出口から室内に吹き出されることで、各室内機３が設置された室内の冷房が行われる。
【0037】
室内熱交換器５１からガス管接続部５４を介してガス管５へ流れる冷媒は、室外機２の第２の閉鎖弁１６を介して室外ガス管２４に流れて四方弁１２の第４のポート１２Ｄに流入する。四方弁１２の第４のポート１２Ｄに流入した冷媒は、第３のポート１２Ｃからアキュムレータ１７の冷媒流入側に流入する。アキュムレータ１７の冷媒流入側から流入した冷媒は、吸入管２２を介して流入し、圧縮機１１に吸入されて再び圧縮されることになる。
【0038】
制御手段１９内の取得部４１は、室外機２内の吐出圧センサ３１、吐出温度センサ３２、吸入圧力センサ３３、吸込温度センサ６３、冷媒温度センサ３５及び外気温度センサ３６のセンサ値を取得する。更に、取得部４１は、各室内機３の液側冷媒温度センサ６１、ガス側温度センサ６２及び吸込温度センサ６３のセンサ値を取得する。
【0039】
図４は、空気調和機１の冷凍サイクルを示すモリエル線図である。空気調和機１の冷房運転時は、室外熱交換器１３が凝縮器として機能し、室内熱交換器５１が蒸発器として機能する。また、空気調和機１の暖房運転時は、室外熱交換器１３が蒸発器として機能し、室内熱交換器５１が凝縮器として機能する。
【0040】
圧縮機１１は、蒸発器から流入する低温低圧のガス冷媒を圧縮して高温高圧のガス冷媒（図４の点Ｂの状態になった冷媒）を吐出する。尚、圧縮機１１が吐出するガス冷媒の温度が吐出温度であり、吐出温度は、吐出温度センサ３２で検出する。
【0041】
凝縮器は、圧縮機１１からの高温高圧のガス冷媒を空気と熱交換して凝縮させる。この際、凝縮器では、潜熱変化によってガス冷媒が全て液冷媒となった後は顕熱変化によって液冷媒の温度が低下して過冷却状態となる（図４の点Ｃの状態）。尚、ガス冷媒が潜熱変化で液冷媒へと変化している際の温度が高圧飽和温度であり、凝縮器の出口における過冷却状態となっている冷媒の温度が熱交出口温度である。高圧飽和温度は、吐出圧力センサ３１で検出した圧力値（図４に「ＨＰＳ」と表記している圧力値Ｐ２）に相当する温度である。熱交出口温度は、冷媒温度センサ３５で検出する。
【0042】
膨張弁は、凝縮器から流出した低温高圧の冷媒を減圧して、ガスと液とが混合した気液二相冷媒（図４の点Ｄの状態になった冷媒）となる。
【0043】
蒸発器は、流入した気液二相冷媒を空気と熱交換して蒸発させる。この際、蒸発器では、潜熱変化によって気液二相冷媒が全てガス冷媒となった後は顕熱変化によってガス冷媒の温度が上昇して過熱状態（図４の点Ａの状態）となり、圧縮機１１に吸入される。尚、液冷媒が潜熱変化でガス冷媒へと変化している際の温度が低圧飽和温度である。低圧飽和温度は、吸入圧力センサ３３で検出した圧力値（図４に「ＬＰＳ」と表記している圧力値Ｐ１）に相当する温度である。また、蒸発器で過熱されて圧縮機１１に吸入される冷媒の温度が吸入温度である。吸入温度は、吸入温度センサ３４で検出する。
【0044】
なお、凝縮器から流出する際に過冷却状態となっている冷媒の冷媒過冷却度は、高圧飽和温度から凝縮器として機能している熱交換器の冷媒出口における冷媒温度（上述した熱交出口温度）を減じて算出できる。また、蒸発器から流出する際に過熱状態となっている冷媒の冷媒過熱度は、低圧飽和温度から吸入温度を減じて算出できる。
【0045】
＜推定モデルの構成＞
推定モデルは、複数の運転状態量の内、任意の運転状態量（特徴量）を用いて回帰分析法の一種である重回帰分析法で生成されている。重回帰分析法では、複数のシミュレーション結果（数値計算により冷媒回路を再現して、残存する冷媒量に対して運転状態量がどのような値となるかを計算した結果）から得られた回帰式のうち、Ｐ値（生成した推定モデルの精度に運転状態量が与える影響度合いを示す値（所定の重みパラメータ））が一番小さく、かつ、補正値Ｒ２（生成した推定モデルの精度を示す値）が０．９以上１．０以下の間のできるだけ大きい値となる回帰式を選択して推定モデルとして生成する。ここで、Ｐ値および補正値Ｒ２は、重回帰分析法で推定モデルを生成する際に、当該推定モデルの精度に関わる値であり、Ｐ値が小さいほど、また、補正値Ｒ２が１．０に近い値であるほど、生成された推定モデルの精度が高くなる。その結果、冷房時の冷媒不足率が０～３０％の場合では、例えば、冷媒過冷却度、外気温度、高圧飽和温度及び圧縮機１１の回転数といった運転状態量を特徴量とする。冷房時の冷媒不足率が４０～７０％の場合では、例えば、吸入温度、外気温度及び圧縮機１１の回転数といった運転状態量を特徴量とする。暖房時の冷媒不足率が０～２０％の場合では、例えば、運転状態量として室外機膨張弁１４の開度を特徴量とする。また、暖房時の冷媒不足率が３０％～７０％の場合では、例えば、吸入冷媒過熱度（吸入温度から低圧飽和温度を減じて求められる）、外気温度、圧縮機１１の回転数及び室外機膨張弁１４といった運転状態量を特徴量とする。
推定モデルは、複数の運転状態量の内、任意の運転状態量（特徴量）を用いて回帰分析法の一種である重回帰分析法で生成されている。重回帰分析法では、複数のシミュレーション結果（数値計算により冷媒回路を再現して、残存する冷媒量に対して運転状態量がどのような値となるかを計算した結果）から得られた回帰式のうち、Ｐ値（生成した推定モデルの精度に運転状態量が与える影響度合いを示す値（所定の重みパラメータ））が一番小さく、かつ、補正値Ｒ２（生成した推定モデルの精度を示す値）が０．９以上１．０以下の間のできるだけ大きい値となる回帰式を選択して推定モデルとして生成する。ここで、Ｐ値および補正値Ｒ２は、重回帰分析法で推定モデルを生成する際に、当該推定モデルの精度に関わる値であり、Ｐ値が小さいほど、また、補正値Ｒ２が１．０に近い値であるほど、生成された推定モデルの精度が高くなる。その結果、冷房時の冷媒不足率が０～３０％の場合では、例えば、冷媒過冷却度、外気温度、高圧飽和温度及び圧縮機１１の回転数といった運転状態量を特徴量とする。冷房時の冷媒不足率が４０～７０％の場合では、例えば、吸入温度、外気温度及び圧縮機１１の回転数といった運転状態量を特徴量とする。暖房時の冷媒不足率が０～２０％の場合では、例えば、運転状態量として室外機膨張弁１４の開度を特徴量とする。また、暖房時の冷媒不足率が３０％～７０％の場合では、例えば、吸入冷媒過熱度（吸入温度から低圧飽和温度を減じて求められる）、外気温度、圧縮機１１の回転数及び室外機膨張弁１４といった運転状態量を特徴量とする。
【0046】
推定モデルは、第１の冷房用推定モデル４３Ａと、第２の冷房用推定モデル４３Ｂと、第３の冷房用推定モデル４３Ｃと、第１の暖房用推定モデル４３Ｄと、第２の暖房用推定モデル４３Ｅと、第３の暖房用推定モデル４３Ｆとを有する。本実施例では、これら各推定モデルは、後述するシミュレーション結果を用いて生成されて、予め空気調和機１の制御手段１９に記憶されている。
【0047】
第１の冷房用推定モデル４３Ａは、冷媒不足率が０％～３０％（第１の範囲）の場合に有効な推定モデルであって、冷媒不足率を高精度に推定できる第１の回帰式である。第１の回帰式は、例えば、（α１×冷媒過冷却度）＋（α２×外気温度）＋（α３×高圧飽和温度）＋（α４×圧縮機１１の回転数）＋α５である。係数α１～α５は、推定モデル生成の際に決定されるものとする。制御部４４は、第１の回帰式に、取得部４１にて取得された現在の冷媒過冷却度、外気温度、高圧飽和温度及び圧縮機１１の回転数を代入することで、現時点での冷媒回路６の冷媒不足率を算出する。尚、冷媒過冷却度、外気温度、高圧飽和温度及び圧縮機１１の回転数を代入する理由は、第１の冷房用推定モデル４３Ａの生成時に使用した特徴量を使用するためである。冷媒過冷却度は、例えば、（高圧飽和温度－熱交出口温度）で算出できる。外気温度は、外気温度センサ３６で検出する。高圧飽和温度は、吐出圧力センサ３１で検出した圧力値を温度変換した値である。圧縮機１１の回転数は、圧縮機１１の図示しない回転数センサで検出する。
【0048】
第２の冷房用推定モデル４３Ｂは、冷媒不足率が４０％～７０％（第２の範囲）の場合に有効な推定モデルであって、冷媒不足率を高精度に推定できる第２の回帰式である。第２の回帰式は、例えば、（α１１×吸入温度）＋（α１２×外気温度）＋（α１３×圧縮機１１の回転数）＋α１４である。係数α１１～α１４は、推定モデル生成の際に決定されるものとする。制御部４４は、第２の回帰式に、取得部４１にて取得された現在の吸入温度、外気温度及び圧縮機１１の回転数を代入することで、現時点での冷媒回路６の冷媒不足率を算出する。尚、吸入温度、外気温度及び圧縮機１１の回転数を代入する理由は、第２の冷房用推定モデル４３Ｂの生成時に使用した特徴量を使用するためである。吸入温度は、吸入温度センサ３４で検出する。外気温度は、外気温度センサ３６で検出する。圧縮機１１の回転数は、圧縮機１１の図示しない回転数センサで検出する。
【0049】
ところで、前述したように、第１の回帰式で求めることができる冷媒不足率は０％～３０％であり、第２の回帰式で求めることができる冷媒不足率は４０％～７０％である。この場合、冷媒不足率が３０％～４０％である場合は、第１の回帰式を用いると冷媒不足率は３０％と算出され、第２の回帰式を用いると冷媒不足率は４０％と算出される。つまり、冷媒不足率が３０％～４０％である場合に、冷媒不足率が３０％以下での寄与度の高い冷媒過冷却度、冷媒不足率が４０％以上での寄与度の高い吸入温度の何れも変化が小さく、有効な推定モデルを生成できない。従って、第１の回帰式あるいは第２の回帰式を用いると、図５Ａに示すようにどちらのモデルを使用するのかによって冷媒不足率が大きく異なる。
【0050】
第３の冷房用推定モデル４３Ｃは、上記のような第１の回帰式あるいは第２の回帰式のいずれを使用しても冷媒不足率を推定できない範囲も含めて、冷媒不足率が０％～７０％の範囲をカバーできる冷房時冷媒不足率算出式である。図５Ｂに示すように、冷房時冷媒不足率算出式は、第１の回帰式の推定結果である冷媒不足率と第２の回帰式の推定結果である冷媒不足率との間を、シグモイド係数を使用したシグモイド曲線で連続的につなぐものである。具体的には、冷房時冷媒不足率算出式は、（シグモイド係数×第１の回帰式で求めた冷媒不足率）＋（（１－シグモイド係数）×第２の回帰式で求めた冷媒不足率）である。制御部４４は、第１の回帰式および第２の回帰式に取得部４１にて取得された現在の運転状態量を代入してそれぞれ算出された冷媒不足率を冷房時冷媒不足率算出式に代入して、現時点での冷媒回路６の冷媒不足率を算出する。
【0051】
ここで、シグモイド係数の算出は、運転状態量のいずれかを用いる。本実施例では、サブクールが０となると第１の回帰式による結果がほぼ一定となってしまうことを考慮し、サブクールが５℃のときに、シグモイド係数が０．５となる計算式とした。
【0052】
ｐ＝１／（１＋ｅｘｐ-（ｓｃ－５））
ｐ:シグモイド係数
ｓｃ:サブクール値
ｐ:シグモイド係数
ｓｃ:サブクール値
【0053】
このようにシグモイド係数を決定して第３の冷房用推定モデル４３Ｃに用いることで、冷媒不足率が０％～３０％、つまり、冷媒不足率が第１の範囲であるときは、第３の冷房用推定モデル４３Ｃによる推定値において第１の推定モデル４３Ａの推定値が支配的となり、また、冷媒不足率が４０％～７０％、つまり、冷媒不足率が第２の範囲であるときは、第３の冷房用推定モデル４３Ｃによる推定値において第２の推定モデル４３Ｂの推定値が支配的となる。
【0054】
なお、シグモイド係数の算出は上述した方法に限らず、実際の冷媒不足率が３０％以上であるとき、つまり、実際の冷媒不足率が第１の範囲でないときは、第３の冷房用推定モデル４３Ｃによる推定値において第２の冷房用推定モデル４３Ｂの推定値が支配的となるように、また、実際の冷媒不足率が４０％以下であるとき、つまり、実際の冷媒不足率が第２の範囲でないときは、第３の冷房用推定モデル４３Ｃによる推定値において第１の冷房用推定モデル４３Ａの推定値が支配的となるように、シグモイド係数を決定すればよい。
【0055】
第１の暖房用推定モデル４３Ｄは、冷媒不足率が０％～２０％（第３の範囲）の場合に有効な推定モデルであって、冷媒不足率を高精度に推定できる第４の回帰式である。第４の回帰式は、例えば、（α３１×室外機膨張弁１４の開度）＋α３２である。制御部４４は、第４の回帰式に、取得部４１にて取得された現在の室外機膨張弁１４の開度を代入することで、冷媒不足率を算出する。尚、室外機膨張弁１４の開度を代入する理由は、第１の暖房用推定モデル４３Ｄの生成時に使用した特徴量を使用するためである。
【0056】
第２の暖房用推定モデル４３Ｅは、冷媒不足率が３０％～７０％（第４の範囲）の場合に有効な推定モデルであって、冷媒不足率を高精度に推定できる第５の回帰式である。第５の回帰式は、例えば、（α４１×吸入過熱度）＋（α４２×外気温度）＋（α４３×圧縮機１１の回転数）＋（α４４×室外機膨張弁１４の開度）＋α４５である。係数α４１～α４５は、推定モデル生成の際に決定されるものとする。制御部４４は、第５の回帰式に、取得部４１にて取得された現在の吸入過熱度、外気温度、圧縮機１１の回転数及びメイン側の膨張弁の開度を代入することで、現時点での冷媒回路６の冷媒不足率を算出する。尚、吸入過熱度、外気温度、圧縮機１１の回転数及びメイン側の膨張弁の開度を代入する理由は、第２の暖房用推定モデル４３Ｅの生成時に使用した特徴量を使用するためである。吸入過熱度は、例えば、（吸入温度－低圧飽和温度）で算出できる。外気温度は、外気温度センサ３６で検出する。圧縮機１１の回転数は、圧縮機１１の図示しない回転数センサで検出する。膨張弁の開度は、図示しないセンサで検出する。
【0057】
また、前述したように、第４の回帰式で求めることができる冷媒不足率は０％～２０％であり、第５の回帰式で求めることができる冷媒不足率は３０％～７０％である。この場合、冷媒不足率が２０％～３０％である場合は、第４の回帰式を用いると冷媒不足率は２０％と算出され、第５の回帰式を用いると冷媒不足率は３０％と算出される。つまり、冷媒不足率が２０％～３０％である場合に、冷媒不足率が２０％以下での寄与度の高い室外機膨張弁１４の開度、冷媒不足率が３０％以上での寄与度の高い吸入過熱度の何れも変化が小さく、有効な推定モデルを生成できない。従って、第４の回帰式あるいは第５の回帰式を用いると、図６Ａに示すようにどちらのモデルを使用するのかによって冷媒不足率が大きく異なる。
【0058】
第３の暖房用推定モデル４３Ｆは、上記のような第４の回帰式あるいは第５の回帰式のいずれを使用しても冷媒不足率を推定できない範囲も含めて、冷媒不足率が０％～７０％の範囲をカバーできる暖房時冷媒不足率算出式である。図６Ｂに示すように、暖房時冷媒不足率算出式は、第４の回帰式の推定結果である冷媒不足率と第５の回帰式の推定結果である冷媒不足率との間を、シグモイド係数を使用したシグモイド曲線で連続的に繋ぐものである。具体的には、暖房時冷媒不足率算出式は、（シグモイド係数×第５の回帰式で求めた冷媒不足率）＋（（１－シグモイド係数）×第４の回帰式で求めた冷媒不足率）である。制御部４４は、第４の回帰式および第５の回帰式に取得部４１にて取得された現在の運転状態量を代入してそれぞれ算出された冷媒不足率を暖房時冷媒不足率算出式に代入して、現時点での冷媒回路６の冷媒不足率を算出する。
【0059】
ここで、シグモイド係数の算出は、冷房運転時と同様に運転状態量のいずれかを用いる。本実施例では、室外膨張弁１４の開度を全閉：０／全開：１００としたときに室外膨張弁１４の開度が全開となると第４の回帰式による結果がほぼ一定となってしまうことを考慮し、室外膨張弁１４の開度が９０のときに、シグモイド係数が０．５となる計算式とした。
【0060】
ｐ＝１／（１＋ｅｘｐ-（Ｄ－９０））
ｐ:シグモイド係数
Ｄ: 室外膨張弁１４の開度
ｐ:シグモイド係数
Ｄ: 室外膨張弁１４の開度
【0061】
このようにシグモイド係数を決定して第３の暖房用推定モデル４３Ｆに用いることで、冷媒不足率が０％～２０％、つまり、冷媒不足率が第３の範囲であるときは、第３の暖房用推定モデル４３Ｆによる推定値において第１の暖房用推定モデル４３Ｄの推定値が支配的となり、また、冷媒不足率が３０％～７０％、つまり、冷媒不足率が第４の範囲であるときは、第３の暖房用推定モデル４３Ｆによる推定値において第２の暖房用推定モデル４３Ｅの推定値が支配的となる。
【0062】
なお、シグモイド係数の算出は上述した方法に限らず、実際の冷媒不足率が２０％以上であるとき、つまり、実際の冷媒不足率が第３の範囲でないときは、第３の暖房用推定モデル４３Ｆによる推定値において第２の暖房用推定モデル４３Ｅの推定値が支配的となるように、また、実際の冷媒不足率が３０％以下であるとき、つまり、実際の冷媒不足率が第４の範囲でないときは、第３の暖房用推定モデル４３Ｆによる推定値において第１の暖房用推定モデル４３Ｄの推定値が支配的となるように、シグモイド係数を決定すればよい。
【0063】
以上に説明したように、冷房運転時は、第１の回帰式、第２の回帰式及び冷房時冷媒不足率算出式を使用して冷媒不足率を推定する。冷房時の冷媒過冷却度が第１の閾値（図５のＴｖ１）より大きい値である場合は、第１の回帰式を選択する方が第２の回帰式を選択するより冷媒不足率を精度よく推定できる。また、冷房時の冷媒過冷却度が第１の閾値より小さい値である場合は、第２の回帰式を選択する方が第１の回帰式を選択するより冷媒不足率を精度よく推定できる。そして、冷房時の冷媒過冷却度が第１の閾値付近の値である場合は、いずれの回帰式を用いるかで冷媒不足率の推定値が大きく変わる。そこで、冷房時は、第１の回帰式と第２の回帰式とを含んだ冷房時冷媒不足率算出式を選択する。これにより、冷房時の冷媒不足率を精度よく推定できる。
【0064】
また、暖房運転時は、第４の回帰式、第５の回帰式及び暖房時冷媒不足率算出式を使用して冷媒不足率を推定する。暖房時のメイン膨張弁の開度が第２の閾値未満（図６のＴｖ２）の場合は、第４の回帰式を選択する方が第５の回帰式を選択するよりもりも冷媒不足率を精度よく推定できる。また、暖房時のメイン膨張弁の開度が第２の閾値未満でない場合は、第５の回帰式を選択する方が第４の回帰式を選択するより冷媒不足率を精度よく推定できる。そして、暖房時のメイン膨張弁の開度が第１の閾値付近の値である場合は、いずれの回帰式を用いるかで冷媒不足率の推定値が大きく変わる。そこで、暖房時は、第４の回帰式と第５の回帰式とを含んだ暖房時冷媒不足率算出式を選択する。これにより、暖房時の冷媒不足率を精度よく推定できる。
【0065】
＜推定処理の動作＞
図７は、推定処理に関わる制御手段１９の処理動作の一例を示すフローチャートである。尚、制御手段１９は、事前に生成された第１の冷房用推定モデル４３Ａ、第２の冷房用推定モデル４３Ｂ、第３の冷房用推定モデル４３Ｃ、第１の暖房用推定モデル４３Ｄ、第２の暖房用推定モデル４３Ｅ、第３の暖房用推定モデル４３Ｆを保持しているものとする。図７において制御手段１９内の制御部４４は、取得部４１を通じて運転状態量を運転データとして収集する（ステップＳ１１）。制御部４４は、収集した運転データから任意の運転状態量を抽出するデータフィルタリング処理を実行する（ステップＳ１２）。制御部４４は、データクレンジング処理を実行する（ステップＳ１３）。制御部４４は、各回帰式又は各冷媒不足率算出式を用いて、現時点の冷媒回路６の冷媒不足率を算出し（ステップＳ１４）、図７に示す処理動作を終了する。
図７は、推定処理に関わる制御手段１９の処理動作の一例を示すフローチャートである。尚、制御手段１９は、事前に生成された第１の冷房用推定モデル４３Ａ、第２の冷房用推定モデル４３Ｂ、第３の冷房用推定モデル４３Ｃ、第１の暖房用推定モデル４３Ｄ、第２の暖房用推定モデル４３Ｅ、第３の暖房用推定モデル４３Ｆを保持しているものとする。図７において制御手段１９内の制御部４４は、取得部４１を通じて運転状態量を運転データとして収集する（ステップＳ１１）。制御部４４は、収集した運転データから任意の運転状態量を抽出するデータフィルタリング処理を実行する（ステップＳ１２）。制御部４４は、データクレンジング処理を実行する（ステップＳ１３）。制御部４４は、各回帰式又は各冷媒不足率算出式を用いて、現時点の冷媒回路６の冷媒不足率を算出し（ステップＳ１４）、図７に示す処理動作を終了する。
【0066】
データフィルタリング処理は、複数の運転状態量の全てを使用するのではなく、所定フィルタ条件に基づき、複数の運転状態量の内、冷媒不足率を算出するのに必要な一部の運転状態量のみを抽出する。生成された推定モデルの各回帰式や各冷媒不足率算出式に、データフィルタリング処理を行った（異常値や突出値を除いた）運転状態量を代入することで、より正確に冷媒不足率を推定できる。
【0067】
所定のフィルタ条件は、第１のフィルタ条件と、第２のフィルタ条件と、第３のフィルタ条件とを有する。第１のフィルタ条件は、例えば、空気調和機１の全運転モード共通に抽出するデータのフィルタ条件である。第２のフィルタ条件は、冷房運転時に抽出するデータのフィルタ条件である。第３のフィルタ条件は、暖房運転時に抽出するデータのフィルタ条件である。
【0068】
第１のフィルタ条件は、例えば、圧縮機１１の駆動状態、運転モードの識別、特殊運転の排除、取得した値における欠損値の排除、各回帰式の生成に際し与える影響の大きい運転状態量について変化量が小さい値の選択、等である。圧縮機１１の駆動状態は、圧縮機が安定して運転していることで冷媒回路６に冷媒が循環していないと冷媒不足率を推定できないために判断する必要のある条件であり、圧縮機１１の立ち上がり時等の過渡期に検出した運転状態量を除外するために設けられるフィルタ条件である。
【0069】
運転モードの識別とは、冷房運転時及び暖房運転時に取得した運転状態量のみを抽出するためのフィルタ条件である。従って、除湿運転時や送風運転時に取得した運転状態量は除外される。特殊運転の排除とは、例えば、油回収運転や除霜運転といった冷房運転時や暖房運転時と比べて冷媒回路６の状態が大きく異なる特殊運転時に取得した運転状態量を除外するフィルタ条件である。欠損値の排除とは、冷媒不足率の判定に使用する運転状態量に欠損値があった場合、当該運転状態量を用いて各回帰式を生成すれば精度が落ちる可能性があるため、欠損値を含む運転状態量を除外するフィルタ条件である。
【0070】
各回帰式や各冷媒不足率算出式に代入する運転状態量について変化量が小さい値の選択とは、空気調和機１の運転状態が安定している状態の運転状態量のみを抽出するフィルタ条件であり、各回帰式や各冷媒不足率算出式による推定精度を上げるために必要な条件である。尚、影響の大きい運転状態量とは、例えば、冷房運転時の冷媒不足率が０～３０％の場合に使用する冷媒過冷却度、冷房運転時の冷媒不足率が４０～７０％の場合に使用する吸入温度や、暖房運転時の吸入過熱度等である。
【0071】
第２のフィルタ条件には、例えば、熱交出口温度の排除、サブクールの異常、吐出温度の異常等がある。
【0072】
熱交出口温度の排除は、外気温度センサ３６と熱交出口温度センサ３５とが近い場所に配置されていることにより、冷房運転時に熱交出口温度センサ３５で検出した熱交出口温度が外気温度センサ３６で検出した外気温度より低くなることがないことを考慮したフィルタ条件であり、外気温度より低い熱交出口温度を除外するフィルタ条件である。
【0073】
サブクール異常は、冷房負荷が極端に大きいあるいは小さいことに起因して異常に高いあるいは以上に低い冷媒過冷却度検出されたときにこれを除外するフィルタ条件である。吐出温度の異常は、冷房負荷が小さいことに起因して圧縮機１１に吸入される冷媒量が減少する所謂ガス欠状態時に検出した吐出温度を除外するフィルタ条件である。
【0074】
第３のフィルタ条件は、例えば、吐出温度の異常等である。暖房運転時に暖房負荷の大きさに起因して吐出温度が高くなって吐出温度保護制御が実行されると、例えば、圧縮機１１の回転数を低下させることで吐出温度が低下するため、このときに検出した吐出温度を除外するフィルタ条件である。
【0075】
データクレンジング処理は、取得した全て運転状態量を冷媒不足率の推定に使用するのではなく、誤った推定を行うおそれがある運転状態量を除外するための処理である。具体的には、取得した運転状態量を平滑化してノイズ抑制やデータ数制限等がある。データの平滑化によるノイズ抑制とは、該当区間の平均値を算出し、各モデルにおいて例えば冷媒過冷却度、吸入温度、吸入冷媒過熱度の移動平均をとることで、ノイズを抑える処理である。データ数制限とは、例えば、データ数が少ないものは信頼性が低いため排除する処理である。例えば、１日分の入力データをフィルタリング処理して残ったデータ数がＸ個以上であれば冷媒不足率の推定に使用、それより少なければ、その日のデータはすべて使用しない。つまり、データクレンジング処理では、推定モデルの各回帰式や各冷媒不足率算出式に異常値や突出値を除いた運転状態量を代入することで、より正確に冷媒不足率を推定できる。
【0076】
図８は、重回帰分析処理に関わる制御手段１９の処理動作の一例を示すフローチャートである。重回帰分析処理とは、例えば、データフィルタリング処理及びデータクレンジング処理後の現在の運転状態量（センサ値）を推定モデルの各回帰式や各冷媒不足率算出式に代入することで、現時点の冷媒回路６の冷媒不足率を算出する処理である。図８において制御手段１９内の制御部４４は、現在が冷房運転中であるか否かを判定する（ステップＳ２１）。制御部４４は、現在が冷房運転中の場合（ステップＳ２１：Ｙｅｓ）、第３の冷房用推定モデル４３Ｃを使用して現在の運転状態量を代入し（ステップＳ２２）、図８に示す処理動作を終了する。その結果、制御部４４は、第３の冷房用推定モデル４３Ｃを使用して現時点での冷媒不足率を算出する。
【0077】
制御部４４は、現在が冷房運転中でない場合（ステップＳ２１：Ｎｏ）、すなわち暖房運転中の場合、第３の暖房用推定モデル４３Ｆを使用して現在の運転状態量を代入し（ステップＳ２３）、図８に示す処理動作を終了する。その結果、制御部４４は、第３の暖房用推定モデル４３Ｆを使用して現時点での冷媒不足率を算出する。
【0078】
＜回帰式の生成方法＞
次に第１～第６の回帰式の生成に使用する特徴量について説明する。第１～第３の回帰式を使用する冷房運転時では、重回帰分析法により第１～第６の回帰式の生成を行う際に使用する特徴量として、例えば、冷媒過冷却度、外気温度、高圧飽和温度、圧縮機１１の回転数、吸入温度等の各運転状態量を用いる。そして、これら各運転状態量は、シミュレーションにより得た結果を使用する。また、第４～第６の回帰式を使用する暖房運転時では、重回帰分析の特徴量として、例えば、吸入過熱度、外気温度、圧縮機１１の回転数、室外機膨張弁１４の開度等の各運転状態量を用いる。そして、これら各運転状態量は、シミュレーションにより得た結果を使用する。
次に第１～第６の回帰式の生成に使用する特徴量について説明する。第１～第３の回帰式を使用する冷房運転時では、重回帰分析法により第１～第６の回帰式の生成を行う際に使用する特徴量として、例えば、冷媒過冷却度、外気温度、高圧飽和温度、圧縮機１１の回転数、吸入温度等の各運転状態量を用いる。そして、これら各運転状態量は、シミュレーションにより得た結果を使用する。また、第４～第６の回帰式を使用する暖房運転時では、重回帰分析の特徴量として、例えば、吸入過熱度、外気温度、圧縮機１１の回転数、室外機膨張弁１４の開度等の各運転状態量を用いる。そして、これら各運転状態量は、シミュレーションにより得た結果を使用する。
【0079】
具体的には、空気調和機１の設計段階で、一例として室内機３が４台運転している場合に外気温度を異ならせてシミュレーションを行い、特徴量と冷媒不足率との関係をシミュレーション毎に取得する。シミュレーションを行う際の条件としては、例えば、外気温度を２０℃、２５℃、３０℃、３５℃及び４０℃と変化させる。なお、シミュレーションを行うに際しては、外気温度の他のパラメータを加えてもよく、例えば、室内機３の運転台数を１～４台と異ならせてもよい。
【0080】
図９は、冷房運転時の室外熱交換機における冷媒出口側の冷媒過冷却度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図９に示す冷媒過冷却度は、冷媒不足率が０％～３０％までは右肩下がりで減少し、冷媒不足率が３０％から６０％までは変化なしとなっている。つまり、冷房運転時に冷媒不足率０～３０％である場合は、冷媒回路６における冷媒量の不足が冷媒過冷却度の値に大きな影響を与えるということである。なお、図９において冷媒不足率が６０％以上であるときの冷媒過冷却度がマイナスの値となっているが、実際は冷媒過冷却度が０℃未満とはならないため、これはシミュレーションでのみ現れる値である。従って、冷媒不足率が６０％以上であるときの冷媒過冷却度は、回帰式の生成に使用しない。
【0081】
図１０は、冷房運転時の吸入温度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図１０に示す吸入温度は、冷媒不足率が４０～７０％のとき増加する傾向にある。つまり、冷房運転時の冷媒不足率が４０～７０％である場合は、冷媒回路６における冷媒量の不足が吸入温度の値に大きな影響を与えるということである。なお、図１０において冷媒不足率が７０％以上であるときの吸入温度はほとんど変化しないため、これ以上の冷媒不足率を吸入温度で推測するのは難しい。従って、冷媒不足率が７０％以上であるときの吸入温度は、回帰式の生成に使用しない。
【0082】
図１１は、暖房運転時の室外機膨張弁１４の開度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図１１に示す室外機膨張弁１４の開度は、冷媒不足率が０～２０％の場合に変化するのに対し、冷媒不足率が２０％を超えると、室外機膨張弁１４の開度の変化が概ね無くなる。つまり、暖房運転時の冷媒不足率が０～２０％である場合は、冷媒回路６における冷媒量の不足が室外機膨張弁１４の開度に大きな影響を与えるということである。なお、上述したように、冷媒不足率が２０％を超えると、室外機膨張弁１４の開度の変化が概ね無くなる。従って、冷媒不足率が２０％を超えたときの室外機膨張弁１４の開度は、回帰式の生成に使用しない。
【0083】
図１２は、吸入過熱度と冷媒不足率の関係についてのシミュレーション結果の一例を示す説明図である。図１２に示す吸入過熱度は、冷媒不足率が増加するときに吸入過熱度が大きくなる傾向にあり、冷媒不足率が３０％を超えると吸入過熱度が大きく上昇する。つまり、暖房運転時の冷媒不足率が３０％を超えると、冷媒回路６における冷媒量の不足が吸入過熱度に大きな影響を与えるということである。なお、図１２において冷媒不足率が３０％より小さいときの吸入過熱度の変化が緩やかであるため、これ以下の冷媒不足率を吸入過熱度で推測するのは難しい。従って、冷媒不足率が３０％より小さい時の吸入過熱度は、回帰式の生成に使用しない。
【0084】
＜実施例１の効果＞
実施例１の空気調和機１では、冷媒回路６に充填される冷媒の冷媒不足率の推定に関わる運転状態量を用いて重回帰分析法で生成された推定モデルと、現在の運転状態量とを用いて、冷媒不足率を推定する。推定モデルを生成する際に使用する運転状態量は、前述したように空気調和機１を様々な環境下で運転した場合を想定したシミュレーションによって求められたものであるため、この推定モデルを用いた冷媒不足率の推定は、いかなる環境下で行っても正確な推定が行える。その結果、冷媒回路６をデフォルト状態に整えることなく、現時点の冷媒不足率を推定できる。
実施例１の空気調和機１では、冷媒回路６に充填される冷媒の冷媒不足率の推定に関わる運転状態量を用いて重回帰分析法で生成された推定モデルと、現在の運転状態量とを用いて、冷媒不足率を推定する。推定モデルを生成する際に使用する運転状態量は、前述したように空気調和機１を様々な環境下で運転した場合を想定したシミュレーションによって求められたものであるため、この推定モデルを用いた冷媒不足率の推定は、いかなる環境下で行っても正確な推定が行える。その結果、冷媒回路６をデフォルト状態に整えることなく、現時点の冷媒不足率を推定できる。
【0085】
空気調和機１に搭載される推定モデルは、複数の運転状態量の内、冷媒回路６に充填される冷媒の冷媒不足率の推定に与える影響の大きい運転状態量を用いて回帰分析法で予め生成される。この推定モデルでは、全ての運転状態量を使用するのではなく、推定モデルに与える影響の大きい運転状態量を選択して推定モデルを生成するため、高精度な推定モデルを生成できる。
【0086】
空気調和機１は、影響の大きい同一種の運転状態量の内、運転状態量の変化量が所定範囲内の運転状態量のみを用いて回帰分析法で生成される。例えば、冷房運転時の冷媒不足率の推定モデルである第１の回帰式を生成する場合は、冷媒不足率が０％～３０％のときの冷媒過冷却度を用いる。また、冷房運転時の冷媒不足率の推定モデルである第２の回帰式を生成する場合は、冷媒不足率が４０％～７０％のときの吸入温度を用いる。その結果、安定した高精度の推定モデルを生成できる。
【0087】
空気調和機１は、冷房運転時の影響の大きい運転状態量として、凝縮器として機能する当該凝縮器出口の冷媒過冷却度、高圧飽和温度、低圧飽和温度及び圧縮機１１の吸入温度を用いて回帰分析法で生成される。その結果、冷房運転時の高精度な冷房用推定モデルを生成できる。
【0088】
空気調和機１は、暖房運転時の影響の大きい運転状態量として、圧縮機１１の吸入過熱度及び膨張弁の開度を用いて回帰分析法で生成される。その結果、暖房運転時の高精度な暖房用推定モデルを生成できる。
【0089】
空気調和機１は、冷房用推定モデルと、冷房運転時の現在の運転状態量とを用いて、冷房運転時の冷媒不足率を推定すると共に、暖房用推定モデルと、暖房運転時の現在の運転状態量とを用いて、暖房運転時の冷媒不足率を推定する。その結果、運転状態毎に異なる推定モデルを使用することで、冷媒不足率を高精度に推定できる。
【0090】
空気調和機１は、第１の冷房用推定モデル４３Ａと第２の冷房用推定モデル４３Ｂとをシグモイド曲線で繋いだ第３の冷房用推定モデル４３Ｃに現在の運転状態量を代入することで、冷房運転時の冷媒不足率を高精度に推定できる。
【0091】
空気調和機１は、第１の暖房用推定モデル４３Ｄと第２の暖房用推定モデル４３Ｅとをシグモイド曲線で繋いだ第３の暖房用推定モデル４３Ｆに現在の運転状態量を代入することで、暖房運転時の冷媒不足率を高精度に推定できる。
【0092】
空気調和機１は、最初に冷媒回路６に充填された冷媒の冷媒量からの冷媒不足率に応じて異なる複数の推定モデルを組み合わせた推定モデルを有する。その結果、空気調和機１は、冷媒不足率を正確に推定できる。
【0093】
第１の冷房用推定モデル４３Ａは、運転状態量として、凝縮器として機能する当該凝縮器の冷媒過冷却度を用いて、冷媒不足率を推定する。その結果、空気調和機１は、冷房運転時に冷媒不足率を高精度に推定できる。
【0094】
第２の冷房用推定モデル４３Ｂは、運転状態量として、圧縮機１１の吸入温度を用いて、冷媒不足率を推定する。その結果、空気調和機１は、冷房運転時に冷媒不足率を高精度に推定できる。
【0095】
第１の暖房用推定モデル４３Ｄは、運転状態量として室外機膨張弁１４の開度を用いて、冷媒不足率を推定する。その結果、空気調和機１は、暖房運転時に冷媒不足率を高精度に推定できる。
【0096】
第２の暖房用推定モデル４３Ｅは、運転状態量として圧縮機１１の吸入過熱度を用いて、冷媒不足率を推定する。その結果、空気調和機１は、暖房運転時に冷媒不足率を高精度に推定できる。
【0097】
第３の冷房用推定モデル４３Ｃは、第１の冷房用推定モデル４３Ａの推定結果と第２の冷房用推定モデル４３Ｂの推定結果との間をシグモイド曲線で補間する。その結果、冷房運転時の冷媒不足率が０～７０％の範囲で、正確な冷媒不足率を推定できる。
【0098】
第３の暖房用推定モデル４３Ｆは、第１の暖房用推定モデル４３Ｄの推定結果と第２の暖房用推定モデル４３Ｅの推定結果との間をシグモイド曲線で補間する。その結果、暖房運転時の冷媒不足率が０～７０％の範囲で、正確な冷媒不足率を推定できる。
【0099】
重回帰分析処理において、データフィルタリング処理及びデータクレンジング処理後の現在の運転状態量（センサ値）を推定モデルの各回帰式に代入する。本実施例では、推定モデルの各回帰式の生成は、シミュレーションで得た特徴量を用いており、シミュレーションで得た特徴量には異常な値や他と比べて突出して大きいあるいは小さい値は含まれていない。このような、異常値や突出値を含まない特徴量を用いて生成された推定モデルの各回帰式や各冷媒不足率算出式に、データフィルタリング処理及びデータクレンジング処理を行って異常値や突出値を除いた運転状態量を代入することで、より正確に冷媒不足率を推定できる。
【0100】
尚、以上に説明した実施例では、空気調和機１の設計段階で各運転状態量のシミュレーション結果を求め、学習機能を有するサーバなどの端末にシミュレーション結果を学習させて得られた推定モデルを制御手段１９が予め記憶している場合を例示した。これに代えて、空気調和機１との間を通信網１１０で接続するサーバ１２０が存在し、このサーバ１２０が第１～第６の回帰式を生成して空気調和機１に送信するようにしてもよい。この実施の形態につき、以下に説明する。
【実施例2】
【0101】
＜空気調和システムの構成＞
図１３は、実施例２の空気調和システム１００の一例を示す説明図である。尚、実施例１の空気調和機１と同一の構成には同一符号を付すことで、その重複する構成及び動作の説明については省略する。図１３に示す空気調和システム１００は、空気調和機１と、通信網１１０と、サーバ１２０とを有する。空気調和機１は、圧縮機１１、室外熱交換器１３及び室外機膨張弁１４を有する室外機２と、室内熱交換器５１を有する室内機３とを有する。空気調和機１は、室外機２と室内機３とが液管４及びガス管５等の冷媒配管で接続されて構成する冷媒回路６を備え、当該冷媒回路６に所定量の冷媒が充填される。
図１３は、実施例２の空気調和システム１００の一例を示す説明図である。尚、実施例１の空気調和機１と同一の構成には同一符号を付すことで、その重複する構成及び動作の説明については省略する。図１３に示す空気調和システム１００は、空気調和機１と、通信網１１０と、サーバ１２０とを有する。空気調和機１は、圧縮機１１、室外熱交換器１３及び室外機膨張弁１４を有する室外機２と、室内熱交換器５１を有する室内機３とを有する。空気調和機１は、室外機２と室内機３とが液管４及びガス管５等の冷媒配管で接続されて構成する冷媒回路６を備え、当該冷媒回路６に所定量の冷媒が充填される。
【0102】
サーバ１２０は、生成部１２１と、送信部１２２とを有する。生成部１２１は、冷媒回路６に充填される冷媒の冷媒不足率の推定に関わる運転状態量を用いて重回帰分析法で推定モデルを生成する。尚、推定モデルは、例えば、第１の実施例で説明した第１の冷房用推定モデル４３Ａ、第２の冷房用推定モデル４３Ｂ、第３の冷房用推定モデル４３Ｃ、第１の暖房用推定モデル４３Ｄ、第２の暖房用推定モデル４３Ｅ及び第３の暖房用推定モデル４３Ｆを有する。送信部１２２は、生成部１２１にて生成した各推定モデルを通信網１１０経由で空気調和機１に送信する。空気調和機１内の制御手段１９は、受信した各推定モデルを用いて空気調和機１の冷媒回路６における冷媒不足率を算出する。
【0103】
サーバ１２０内の生成部１２１は、冷媒回路６における冷媒不足率を実測できる空気調和機１の標準機（製造メーカーの試験室などに設置されている）から定期的に冷房運転時の運転状態量を収集し、各推定モデルで推定した冷媒不足率と実測した冷媒不足率との比較結果と収集した運転状態量とを用いて、第１の冷房用推定モデル４３Ａ、第２の冷房用推定モデル４３Ｂ及び第３の冷房用推定モデル４３Ｃを生成あるいは更新する。そして、サーバ１２０内の送信部１２２は、生成あるいは更新した第１の冷房用推定モデル４３Ａ、第２の冷房用推定モデル４３Ｂ及び第３の冷房用推定モデル４３Ｃを空気調和機１に定期的に送信する。なお、実施例１のように、各推定モデルの生成に使用する運転状態量をシミュレーションで得て、生成部１２１がシミュレーションで得た運転状態量を用いて各推定モデルを生成してもよい。
【0104】
サーバ１２０内の生成部１２１は、上述した空気調和機１の標準機から定期的に暖房運転時の運転状態量を収集し、各推定モデルで推定した冷媒不足率と実測した冷媒不足率との比較結果と収集した運転状態量とを用いて、第１の暖房用推定モデル４３Ｄ、第２の暖房用推定モデル４３Ｅ及び第３の暖房用推定モデル４３Ｆを生成する。そして、サーバ１２０内の送信部１２２は、生成した第１の暖房用推定モデル４３Ｄ、第２の暖房用推定モデル４３Ｅ及び第３の暖房用推定モデル４３Ｆを空気調和機１に定期的に送信する。なお、実施例１のように、各推定モデルの生成に使用する運転状態量をシミュレーションで得て、生成部１２１がシミュレーションで得た運転状態量を用いて各推定モデルを生成してもよい。
【0105】
＜実施例２の効果＞
実施例２のサーバ１２０は、冷媒回路６に充填される冷媒の冷媒不足率の推定に関わる運転状態量を用いて重回帰分析法を使用して、冷媒不足率を推定する推定モデルを生成し、生成した推定モデルを空気調和機１に送信する。空気調和機１は、サーバ１２０から受信した推定モデルと、現在の運転状態量とを用いて、冷媒不足率を推定する。その結果、空気調和機１では、高精度な推定モデルを生成できる。
実施例２のサーバ１２０は、冷媒回路６に充填される冷媒の冷媒不足率の推定に関わる運転状態量を用いて重回帰分析法を使用して、冷媒不足率を推定する推定モデルを生成し、生成した推定モデルを空気調和機１に送信する。空気調和機１は、サーバ１２０から受信した推定モデルと、現在の運転状態量とを用いて、冷媒不足率を推定する。その結果、空気調和機１では、高精度な推定モデルを生成できる。
【0106】
尚、実施例１及び２の空気調和機１では、１台の室外機２に対してＮ台の室内機３を接続した場合の冷媒不足率を推定する推定モデルを例示した。これに対し、１台の室外機２と１台の室内機３とが接続した空気調和機１についても、実施例１や実施例２と同様の方法で冷媒不足率を推定できる。上記のような空気調和機１につき、実施例３として以下に説明する。
【実施例3】
【0107】
制御手段は、室外機：室内機が１：１の場合、冷房運転時の現時点の冷媒不足率を推定する第４の冷房用推定モデルと、暖房運転時の現時点の冷媒不足率を推定する第５の暖房用推定モデルとを有する。尚、説明の便宜上、実施例１の空気調和機１と同一の構成には同一符号を付すことで、その重複する構成及び動作の説明については省略する。実施例１の空気調和機１と実施例３の空気調和機１とが異なるところは、室内機３を１台にし、第１～第３の冷房用推定モデルとは異なる運転状態量を用いて生成された第４の冷房用推定モデルを使用し、第１～第３の暖房用推定モデルとは異なる運転状態量を用いて生成された第４の暖房用推定モデルを使用する点にある。
【0108】
第４の冷房用推定モデルは、重回帰分析法により生成された第７の回帰式である。第７の回帰式は、例えば、（α７１×室外熱交温度）－（α７２×外気温度）－（α７３×吐出温度）＋（α７４×圧縮機１１の回転数）－（α７５×膨張弁の開度）＋α７６である。係数α７１～α７５は、推定モデル生成の際に決定されるものとする。制御部４４は、第７の回帰式に、取得部４１にて取得された現在の運転状態量、例えば、室外熱交温度、外気温度、吐出温度、圧縮機１１の回転数及び膨張弁の開度を代入することで、現時点での冷媒不足率を算出する。尚、室外熱交温度、外気温度、吐出温度、圧縮機１１の回転数及び膨張弁の開度を代入する理由は、第４の冷房用推定モデルの生成時に使用した特徴量を使用するためである。尚、室外熱交温度は、冷媒温度センサ３５で検出する。
【0109】
第４の暖房用推定モデルは、重回帰分析法により生成された第８の回帰式である。第８の回帰式は、例えば、（α８１×室内熱交温度）＋（α８２×圧縮機１１の回転数）＋（α８３×外気温度）－（α８４×室外熱交温度）－（α８５×膨張弁の開度）＋α８６である。係数α８１～α８５は、推定モデル生成の際に決定されるものとする。制御部４４は、第８の回帰式に、取得部４１にて取得された現在の運転状態量、例えば、室内熱交温度、圧縮機１１の回転数、外気温度、室外熱交温度、外気温度、吐出温度及び膨張弁の開度を代入することで、現時点での冷媒不足率を算出する。尚、室内熱交温度、圧縮機１１の回転数、外気温度、室外熱交温度、外気温度、吐出温度及び膨張弁の開度を代入する理由は、第４の暖房用推定モデルの生成時に使用した特徴量を使用するためである。尚、暖房時の室内熱交温度は、吐出圧力センサ３１で検出した圧力値から換算できる。
【0110】
また、本実施例では、冷媒回路６に残存する冷媒量を表すものとして相対的な冷媒量を推定する場合を説明した。具体的には、冷媒回路６に冷媒を充填した際の充填量（初期値）に対する、冷媒回路６から外部に漏洩した冷媒量の割合である冷媒不足率を推定して提供する場合を説明した。しかし、本発明はこれに限られるものではなく、推定した冷媒不足率に初期値を乗じて、冷媒回路６から外部に漏洩した冷媒量を提供するようにしてもよい。また、冷媒回路６から外部に漏洩した絶対的な冷媒量あるいは冷媒回路６に残留する絶対的な冷媒量を推定する推定モデルを生成し、この推定モデルによる推定結果を提供するようにしてもよい。冷媒回路６から外部に漏洩した絶対的な冷媒量あるいは冷媒回路６に残留する絶対的な冷媒量を推定する推定モデルを生成する場合は、ここまでに説明した各運転状態量に加えて、室外熱交換器１３および各室内熱交換器１の容積や液管４の容積を考慮すればよい。
【0111】
＜変形例＞
尚、本実施例では、例えば、第１の冷房用推定モデル４３Ａの推定結果と第２の冷房用推定モデル４３Ｂの推定結果との間をシグモイド係数で補間する場合を例示したが、シグモイド係数に限定されるものではなく、例えば、線形補間等の補間方法を使用しても良く、適宜変更可能である。
尚、本実施例では、例えば、第１の冷房用推定モデル４３Ａの推定結果と第２の冷房用推定モデル４３Ｂの推定結果との間をシグモイド係数で補間する場合を例示したが、シグモイド係数に限定されるものではなく、例えば、線形補間等の補間方法を使用しても良く、適宜変更可能である。
【0112】
本実施例では、複数のシミュレーション結果の内、全てのシミュレーション結果を使用するのではなく、一部のシミュレーション結果を使用する。例えば、冷房運転時の冷媒不足率が０～３０％の場合に使用する第１の冷房用推定モデル４３Ａ、冷媒不足率が４０～７０％の場合に使用する第２の冷房用推定モデル４３Ｂ、冷媒不足率３０～４０％の場合に使用する第３の冷房用推定モデル４３Ｃのように個々に分けて生成する。従って、運転状態量をシミュレーションで用意するため、空気調和機１を動作して運転状態量を収集する場合に比較して簡単かつ必要な量の運転状態量を収集できる。
【0113】
本実施例では、推定モデルは、サーバで生成する場合を例示したが、利用者がシミュレーション結果から推定モデルを算出しても良い。また、本実施例では、重回帰分析法を用いて各推定モデルを生成する場合を例示したが、一般の回帰分析法を行える機械学習手法のSVR(Support Vector Regression)、NN(Neural Network)などを用いて推定モデルを生成しても良い。その際、特徴量選択に当たっては重回帰分析法で用いたＰ値や補正値Ｒ２の代わりに、推定モデルの精度が向上するよう特徴量を選択する一般の手法（Forward Feature Selection法、Backward feature Eliminationなど）を使えばよい。
【0114】
また、図示した各部の各構成要素は、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各部の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況等に応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。
【0115】
更に、各装置で行われる各種処理機能は、ＣＰＵ（Central Processing Unit）（又はＭＰＵ（Micro Processing Unit）、ＭＣＵ（Micro Controller Unit）等のマイクロ・コンピュータ）上で、その全部又は任意の一部を実行するようにしても良い。また、各種処理機能は、ＣＰＵ（又はＭＰＵ、ＭＣＵ等のマイクロ・コンピュータ）で解析実行するプログラム上、又はワイヤードロジックによるハードウェア上で、その全部又は任意の一部を実行するようにしても良いことは言うまでもない。
【0116】
また、以上に説明した各実施例では、冷媒不足率を、冷媒が規定量充填されているときを１００％としたとき、この規定量からの減少分とした。これに代えて、冷媒回路６に冷媒を規定量充填した直後に、本実施例に記載した方法で冷媒不足率を推定し、この推定結果を１００％としてもよい。例えば、冷媒回路６に冷媒を規定量充填した直後に推定した冷媒不足率が９０％である場合、つまり、冷媒回路６に充填されている冷媒量が規定量充填より１０％少ないと推定した場合は、この規定量充填より１０％少ない冷媒量と１００％としてもよい。このように１００％とする冷媒量を推定結果に合わせることで、これ以降の冷媒不足率をより正確に推定できる。
【符号の説明】
【0117】
１ 空気調和機
２ 室外機
３ 室内機
４ 液管
５ ガス管
１１ 圧縮機
１２ 四方弁
１３ 室外熱交換器
１４ 室外機膨張弁
４１ 取得部
４３Ａ 第１の冷房用推定モデル
４３Ｂ 第２の冷房用推定モデル
４３Ｃ 第３の冷房用推定モデル
４３Ｄ 第１の暖房用推定モデル
４３Ｅ 第２の暖房用推定モデル
４３Ｆ 第３の暖房用推定モデル
４４ 制御部
５１ 室内熱交換器
２ 室外機
３ 室内機
４ 液管
５ ガス管
１１ 圧縮機
１２ 四方弁
１３ 室外熱交換器
１４ 室外機膨張弁
４１ 取得部
４３Ａ 第１の冷房用推定モデル
４３Ｂ 第２の冷房用推定モデル
４３Ｃ 第３の冷房用推定モデル
４３Ｄ 第１の暖房用推定モデル
４３Ｅ 第２の暖房用推定モデル
４３Ｆ 第３の暖房用推定モデル
４４ 制御部
５１ 室内熱交換器
【図1】
【図2】
【図3】
【図4】
【図5A】
【図5B】
【図6A】
【図6B】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】