(57)【特許請求の範囲】
【請求項１】
室外機を一台または複数台と、室内機を一台または複数台設け、前記室外機と前記室内機とを配管接続して閉回路をなし、前記閉回路の中に冷媒を封入し、
前記室内機においては、室内空気と熱交換を行う室内熱交換器と冷媒の流量を調節する室内膨張弁を順次配管接続し、
前記室外機においては、圧縮機と第一室外熱交換器及び室外膨張弁を配管接続してヒートポンプ回路を形成すると共に、前記圧縮機の低圧配管と前記室内膨張弁と第二室外熱交換器を配管接続して冷媒自然循環回路を形成し、
空調負荷に応じて前記圧縮機を用いる蒸気圧縮式ヒートポンプ運転と、前記冷媒自然循環回路を用いる冷媒自然循環式サイクル運転とを切換えて運転する制御演算装置を備えた空気調和機において、
前記制御演算装置は、空調負荷を検出する空調場熱負荷同定部と、蒸気圧縮式ヒートポンプの運転能力を検知するヒートポンプ能力検知部と、冷媒自然循環式サイクルの運転能力を予測する自然サイクル能力予測部と、冷媒自然循環式サイクルの運転中の運転能力を検知する自然サイクル能力検知部と、前記空調場熱負荷同定部で検出された空調負荷と冷媒自然循環式サイクルの運転能力とを比較する判定部を設け、
前記制御演算装置は、前記判定部で前記空調場熱負荷同定部で検知された空調負荷より冷媒自然循環式サイクルの運転能力の予測値が大きいと判定したとき冷媒自然循環式サイクル運転を起動し、この起動後に前記空調場熱負荷同定部で空調負荷を再度検知するとともに、前記自然サイクル能力検知部で冷媒自然循環式サイクル運転中の運転能力を検知し、前記再検出された空調負荷より大きいと前記判定部で判定したとき冷媒自然循環式サイクル運転を継続するように構成されたことを特徴とする空気調和機。
【請求項２】
請求項１に記載の空気調和機において、前記制御演算装置は、前記判定部で冷媒自然循環式サイクルの運転中の運転能力が予測値より大きいと判定したとき、前記空調場熱負荷同定部で空調負荷を再度検知することを特徴とする空気調和機。
【請求項３】
請求項１または２に記載の空気調和機において、前記制御演算装置は、さらに、前記判定部で冷媒自然循環式サイクルの運転中の運転能力が予測値より小さいと判定したとき冷媒自然循環式サイクルの運転状況を改善する運転状況改善部を備え、この運転状況改善部での改善の後、再度冷媒自然循環式サイクル運転を継続するように構成されたことを特徴とする空気調和機。
【請求項４】
請求項１に記載の空気調和機において、前記制御演算装置は空調する部屋の熱容量並びに熱通過係数を少なくとも要因とする外部空調負荷に関わるパラメータ及び内部空調負荷を、複数の検知による計測量を元に時々刻々同定し、経年変化や状態変化に追随して切り替え運転することを特徴とする空気調和機。
【請求項５】
請求項１に記載の空気調和機において、前記制御演算装置は空気温度条件と前記冷媒自然循環式サイクルの運転能力の関係を表す係数を、複数の検知による計測量を元に時々刻々同定し、経年変化や状態変化に追随して切り替え運転することを特徴とする空気調和機。
【請求項６】
室外機を一台または複数台と、室内機を一台または複数台設け、前記室外機と前記室内機とを配管接続して閉回路をなし、前記閉回路の中に冷媒を封入し、
前記室内機においては、室内空気と熱交換を行う室内熱交換器と冷媒の流量を調節する室内膨張弁を順次配管接続し、
前記室外機においては、圧縮機と第一室外熱交換器及び室外膨張弁を配管接続してヒートポンプ回路を形成すると共に、前記圧縮機の低圧配管と前記室内膨張弁と第二室外熱交換器を配管接続して冷媒自然循環回路を形成し、
制御演算装置により、空調負荷に応じて前記圧縮機を用いるヒートポンプ運転と、前記冷媒自然循環回路を用いる冷媒自然循環式サイクル運転とを切換えて運転する空気調和機の運転方法において、
前記制御演算装置は、空調負荷を検出する空調場熱負荷同定部と、ヒートポンプ運転の能力を検知するヒートポンプ能力検知部と、冷媒自然循環式サイクルの運転能力を予測する自然サイクル能力予測部と、冷媒自然循環式サイクルの運転中の運転能力を検知する自然サイクル能力検知部と、前記空調場熱負荷同定部で検出された空調負荷と冷媒自然循環式サイクルの運転能力とを比較する判定部を設け、
前記制御演算装置により、検知された空調負荷より冷媒自然循環式サイクルの運転能力の予測値が大きいと判定したとき冷媒自然循環式サイクル運転を起動し、この起動後に空調負荷を再度検知するとともに冷媒自然循環式サイクル運転中の運転能力を検知し、前記再検出された空調負荷より冷媒自然循環式サイクル運転中の運転能力が大きいと判定したとき冷媒自然循環式サイクル運転を継続することを特徴とする空気調和機の運転方法。
【請求項７】
請求項６に記載の空気調和機の運転方法において、前記制御演算装置により、冷媒自然循環式サイクルの運転中の運転能力が予測値より大きいと判定したとき、空調負荷を再度検知することを特徴とする空気調和機の運転方法。
【請求項８】
請求項６または７に記載の空気調和機の運転方法において、前記制御演算装置により、さらに、冷媒自然循環式サイクルの運転中の運転能力が予測値より小さいと判定したとき冷媒自然循環式サイクルの運転状況を改善し、この改善の後再度冷媒自然循環式サイクル運転を継続することを特徴とする空気調和機の運転方法。
【請求項９】
請求項６に記載の空気調和機の運転方法において、前記制御演算装置により空調する部屋の熱容量並びに熱通過係数を少なくとも要因とする外部空調負荷に関わるパラメータ及び内部空調負荷を、複数の検知による計測量を元に時々刻々同定し、経年変化や状態変化に追随して切り替え運転することを特徴とする空気調和機の運転方法。
【請求項１０】
請求項６に記載の空気調和機の運転方法において、前記制御演算装置により空気温度条件と前記冷媒自然循環式サイクルの運転能力の関係を表す係数を、複数の検知による計測量を元に時々刻々同定し、経年変化や状態変化に追随して切り替え運転することを特徴とする空気調和機の運転方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、蒸気圧縮式ヒートポンプと冷媒自然循環式サイクルを併用する空気調和機およびその運転方法に関する。
【背景技術】
【０００２】
近年、社会の省エネルギの推進が要望されている。この要望より、近年製造される空気調和機には省エネルギが実現出来る様な様々な工夫がなされている。空気調和機の実使用としては、種々の研究があるが、定格負荷の１０〜２０％の低負荷にて運転されていると言う報告もある。従って、低負荷に対する運転を向上させる事で、実際の運転効率向上が見込める。低負荷に対する運転を向上させる手段として、特許文献１や、特許文献２にある様に、圧縮機の起動停止（発停）頻度を低減させる技術がある。また駆動エネルギが不必要な冷媒自然循環式サイクルを併用し、蒸気圧縮式ヒートポンプから切り換えて、省エネルギを図る方法がある。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００２−６１９２５号公報
【特許文献２】特開２００９−２２８９２２号公報
【特許文献３】特開２００５−２５７１９７号公報
【特許文献４】特開平１１−２８７５２４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
上記特許文献１のものは、圧縮機の運転時間が短くなると、冷房運転時の目標値を上げることにより、目標値と実際の蒸発温度との偏差を小さくする様に構成されており、上記目標変更回路は、前回の変更量を考慮して今回の変更量を導出するため、調整がなされるためには、何度も圧縮機の発停が必要となり、その分、省エネルギ性の向上分は制限される。
【０００５】
また特許文献２のものは、対象を冷蔵庫としているが、圧縮機が停止した際、圧縮機の平均回転数（以下、周波数と呼ぶ）を算出し、それが最適周波数下限より小さいか、或いは最適周波数上限より大きいかを判定し、低圧圧力の目標値を変更する。しかし、これも圧縮機が発停する必要があり、効率向上分が制限される。
【０００６】
これらの技術は、現地の最適値に調整されていなかったパラメータを調整する技術形態の一つであるが、空気調和機の能力と負荷を合致させると言う考え方ではないため、例えば日中、太陽が昇って室外気温が上昇したり、日射が増えて負荷環境が変化していくと、それを前もって補正したり、素早く追従する事は出来ない。また空気調和機の運転調整に留まっているため、効率を決定する重要要素である圧縮機運転範囲の限界が空気調和機の限界となってしまう。
【０００７】
また自然循環式サイクルを併用するものとして、特許文献３がある。特許文献３は主に判定後の運転制御法に関するものであるが、自然循環式サイクルが使用出来ると判定する方式として、主に室外温度と室内温度を比較して、室外温度が室内温度より（一定値以上）低いと判断した時に、蒸気圧縮式ヒートポンプから自然循環に切り替える方式を取っている。
【０００８】
しかし、確かに室外温度が室内温度より低ければ自然循環式サイクルが作動出来るが、実際の運転では、自然循環式サイクルの能力が低く、空調負荷よりも小さい場合、室内温度が上昇し、著しく快適性に欠ける運転となってしまう。対策として、室外温度と室内温度により、予め自然循環式サイクルの能力を計算して記憶しておき、自然循環式サイクルに切り替わった際の能力を予測する方法があるが、空調負荷よりも大きいかどうかが判定出来ないため、室内温度を快適に保持する事は出来ない。
【０００９】
また特許文献４の様に、室内負荷を室内外温度センサにより推定し、自然循環式サイクルの能力に応じて蒸気圧縮式ヒートポンプと切り換える方式もある。この方式は、室内負荷を温度を元に仮定し、前回の運転状態と比較している。ある検知された温度状態に対し、前回蒸気圧縮式ヒートポンプで運転していた場合は今回も蒸気圧縮式ヒートポンプとし、前回自然循環式サイクルで運転していた場合には、自然循環式サイクルと判定する。
【００１０】
しかし、実際の空調負荷は、例えば冷房負荷の場合、外気温が上昇した際には増加し、オフィス機器等が多く稼働すれば増加する。また室内の人間が増えた場合や日照量や換気量が増えた場合も増加する。この様に、空調負荷は運転中に時々刻々と変化するものであって、一定では無く、前回の運転状態は参考とはなるが、現在の空調負荷を決める事は出来ない。
【００１１】
そこで本発明は、蒸気圧縮式ヒートポンプから自然循環式サイクルを切り替える際、時々刻々と変化する空調負荷を複数回の同定により検知し、現在の空調負荷を把握した上で、自然循環式サイクルでの能力と比べ、自然循環式サイクルにて空調負荷を賄えると判断した時に、自然循環式サイクルに切り替える空気調和機および空気調和機の運転方法を提供する。
【課題を解決するための手段】
【００１２】
上記目的を達成するために、本発明は、室外機を一台または複数台と、室内機を一台または複数台設け、前記室外機と前記室内機とを配管接続して閉回路をなし、前記閉回路の中に冷媒を封入し、
前記室内機においては、室内空気と熱交換を行う室内熱交換器と冷媒の流量を調節する室内膨張弁を順次配管接続し、
前記室外機においては、圧縮機と第一室外熱交換器及び室外膨張弁を配管接続してヒートポンプ回路を形成すると共に、前記圧縮機の低圧配管と前記室内膨張弁と第二室外熱交換器を配管接続して冷媒自然循環回路を形成し、
空調負荷に応じて前記圧縮機を用いるヒートポンプ運転と、前記冷媒自然循環回路を用いる冷媒自然循環式サイクル運転とを切換えて運転する制御演算装置を備えた空気調和機において、
前記制御演算装置は、空調負荷を検出する空調場熱負荷同定部と、ヒートポンプ運転の能力を検知するヒートポンプ能力検知部と、冷媒自然循環式サイクルの運転能力を予測する自然サイクル能力予測部と、冷媒自然循環式サイクルの運転中の運転能力を検知する自然サイクル能力検知部と、前記空調場熱負荷同定部で検出された空調負荷と冷媒自然循環式サイクルの運転能力とを比較する判定部を設け、
前記制御演算装置は、前記判定部で前記空調場熱負荷同定部で検知された空調負荷より冷媒自然循環式サイクルの運転能力の予測値が大きいと判定したとき冷媒自然循環式サイクル運転を起動し、この起動後に前記空調場熱負荷同定部で空調負荷を再度検知するとともに、前記自然サイクル能力検知部で冷媒自然循環式サイクル運転中の運転能力を検知し、前記再検出された空調負荷より大きいと前記判定部で判定したとき冷媒自然循環式サイクル運転を継続するように構成されたことを特徴とする。
【００１３】
また、上記に記載の空気調和機において、前記制御演算装置は、前記判定部で冷媒自然循環式サイクルの運転中の運転能力が予測値より大きいと判定したとき、前記空調場熱負荷同定部で空調負荷を再度検知することを特徴とする。
【００１４】
また、上記に記載の空気調和機において、前記制御演算装置は、さらに、前記判定部で冷媒自然循環式サイクルの運転中の運転能力が予測値より小さいと判定したとき冷媒自然循環式サイクルの運転状況を改善する運転状況改善部を備え、この運転状況改善部での改善の後、再度冷媒自然循環式サイクル運転を継続するように構成されたことを特徴とする。
【００１５】
また、上記に記載の空気調和機において、前記制御演算装置は、空調する部屋の熱容量並びに熱通過係数を少なくとも要因とする外部空調負荷に関わるパラメータ及び内部空調負荷を、複数の検知による計測量を元に時々刻々同定し、経年変化や状態変化に追随して切り替え運転することを特徴とする。
【００１６】
また、上記に記載の空気調和機において、前記制御演算装置は、空気温度条件と前記冷媒自然循環式サイクルの運転能力の関係を表す係数を、複数の検知による計測量を元に時々刻々同定し、経年変化や状態変化に追随して切り替え運転することを特徴とする。
【００１７】
上記目的を達成するために、本発明は、室外機を一台または複数台と、室内機を一台または複数台設け、前記室外機と前記室内機とを配管接続して閉回路をなし、前記閉回路の中に冷媒を封入し、
前記室内機においては、室内空気と熱交換を行う室内熱交換器と冷媒の流量を調節する室内膨張弁を順次配管接続し、
前記室外機においては、圧縮機と第一室外熱交換器及び室外膨張弁を配管接続してヒートポンプ回路を形成すると共に、前記圧縮機の低圧配管と前記室内膨張弁と第二室外熱交換器を配管接続して冷媒自然循環回路を形成し、
制御演算装置により、空調負荷に応じて前記圧縮機を用いるヒートポンプ運転と、前記冷媒自然循環回路を用いる冷媒自然循環式サイクル運転とを切換えて運転する空気調和機の運転方法において、
前記制御演算装置は、空調負荷を検出する空調場熱負荷同定部と、ヒートポンプ運転の能力を検知するヒートポンプ能力検知部と、冷媒自然循環式サイクルの運転能力を予測する自然サイクル能力予測部と、冷媒自然循環式サイクルの運転中の運転能力を検知する自然サイクル能力検知部と、前記空調場熱負荷同定部で検出された空調負荷と自然循環サイクルの運転能力とを比較する判定部を設け、
前記制御演算装置により、検知された空調負荷より冷媒自然循環式サイクルの運転能力の予測値が大きいと判定したとき冷媒自然循環式サイクル運転を起動し、この起動後に空調負荷を再度検知するとともに冷媒自然循環式サイクル運転中の運転能力を検知し、前記再検出された空調負荷より冷媒自然循環式サイクル運転中の運転能力が大きいと判定したとき冷媒自然循環式サイクル運転を継続することを特徴とする。
【発明の効果】
【００１８】
本発明によれば、通常運転中に部屋の熱容量並びに熱通過係数等の外部空調負荷に関わるパラメータや内部空調負荷を時々刻々同定し、部屋の空調負荷に合った能力を発揮する事で圧縮機の頻繁な発停を防止し、冷媒自然循環式サイクルの空調能力が空調負荷を十分賄える場合、冷媒自然循環式サイクルを利用して圧縮機を停止させる事で、圧縮機の頻繁な発停を防止し、効率の良い運転を行う事が出来る。
【図面の簡単な説明】
【００１９】
【図１】本発明実施例の冷凍サイクル系統図である。
【図２】図1の制御演算装置の構成図である。
【図３Ａ】実施例１の動作を全熱負荷で示す第1フローチャートである。
【図３Ｂ】同じく第２フローチャートである。
【図４Ａ】実施例２の動作を顕熱負荷で示す第1フローチャートである。
【図４Ｂ】同じく第２フローチャートである。
【図５】本発明実施例の空調全熱内部負荷と同負荷の推定値を示すグラフである。
【図６】同じく外部負荷パラメータａの推定値を表すグラフである。
【図７】同じく外部負荷パラメータｂの推定値を表すグラフである。
【図８】同じくパラメータｃの推定値を表すグラフである。
【図９】同じくパラメータｄの推定値を表すグラフである。
【発明を実施するための形態】
【００２０】
以下、本発明の実施の形態１を説明する。図１は、空気調和機と制御演算装置の構成を示す冷凍サイクル系統図であり、１台或いは複数台の室外機２０と、１台或いは複数台の室内機３０を有し、室外機２０と室内機３０を配管２８、２９で接続して閉回路をなし、その閉回路の中に冷媒を封入している。
【００２１】
室外機２０は、１台或いは複数台の、運転回転数（以下周波数という）可変、或いは周波数固定の圧縮機（蒸気圧縮式ヒートポンプ）２１と、第一室外熱交換器２２１と、第二室外熱交換器２２２と、これらの室外熱交換器２２１、２２２の冷媒流量を調整する室外膨張弁２７とを配管接続すると共に、室外熱交換器２２１、２２２に送風する室外ファン２３を備えている。
【００２２】
室内機３０は、室内空気と熱交換を行う室内熱交換器３１とその室内熱交換器の冷媒流量を調整する室内膨張弁３３を順次配管すると共に、室内熱交換器３１に送風する室内ファン３２が設けられている。また室内機３０は空気調和の対象となる部屋に配置してある。
【００２３】
室外機２０は、四方弁２５を有し、アキュムレータ２４と液受容器２６を備えているが、これらは無くても使用可能である。また室外機２０及び室内機３０の各液側及びガス側を、各液側管路２８、ガス側管路２９で接続して閉回路となし、その閉回路の中に冷媒が封入してある。そして、圧縮機２１と第一室外熱交換器２２１及び室外膨張弁３６を配管接続してヒートポンプ回路を形成すると共に、圧縮機２１の低圧配管と前記室内膨張弁３３と第二室外熱交換器２２２を配管接続して冷媒自然循環回路を形成している。
【００２４】
さらに、室外機２０には、室外温度を検知する室外温度検知器４４、圧縮機冷媒吸入温度検知器４０、及び圧縮機吐出温度検知器４１、圧縮機冷媒吸入圧力を検知する圧縮機吸入圧力検知器４５、圧縮機冷媒吐出圧力を検知する圧縮機吐出圧力検知器４６が設けられ、圧縮機の周波数を操作するインバータ圧縮機周波数操作器２１、室外ファン２３の送風能力を操作する室外ファン送風能力操作器３５、室外膨張弁２７開度を操作する室外膨張弁開度操作器３６、蒸気圧縮式ヒートポンプ２１と冷媒自然循環式サイクルの第二室外熱交換器２２２を変えた方が良い時に、使用熱交換器を変更する第一の開閉弁５４と第二の開閉弁５５がそれぞれ設けられている。
【００２５】
前記室外機２０には、圧縮機２１、第一室外熱交換器２２１、室外膨張弁２７、室内膨張弁３３、および室内熱交換器３１を配管接続してヒートポンプ回路が形成される。また、圧縮機２１の低圧配管、第二室外熱交換器２２２、室内膨張弁３３、および室内熱交換器３１を配管接続して冷媒自然循環回路が形成される。
【００２６】
室内機３０には、室内温度を検知する室内機吸込温度検知器４８、その部屋への吹出空気温度を検知する室内機吹出温度検知器４９、室内ファン３２の送風能力を操作する室内ファン送風能力操作器３８、室内膨張弁３３の開度を操作する室内膨張弁開度操作器３９、予め与えられた設定値を記憶或いは使用者が好みの熱環境を設定するための室内温度設定器５６を有している。さらに、本実施例の空気調和機の動作手順を演算して制御する制御演算装置５７が接続されている。
【００２７】
上記構成において、前記制御演算装置５７により、空調負荷に応じて前記圧縮機２１を用いるヒートポンプ運転と、前記冷媒自然循環回路を用いる冷媒自然循環式サイクル運転とを切換えて運転がなされる。
【００２８】
図２は、前記制御演算装置５７の構成図である。制御演算装置５７は、この装置全体を制御する制御部１００と、空調する場所の空調熱負荷を検出する空調場熱負荷同定部１０１と、ヒートポンプ２１の運転能力を検知するヒートポンプ能力検知部１０２と、外気温度等から冷媒自然循環回路による冷媒自然循環式サイクルの運転能力を予測する自然サイクル能力予測部１０４と、冷媒自然循環式サイクルの運転中の運転能力を検知する自然サイクル能力検知部１０３と、前記空調場熱負荷同定部１０１で検出された空調負荷と冷媒自然循環式サイクルの運転能力の予測値、または運転能力とを比較する判定部１０６を備えている。判定部１０６は、冷媒自然循環式サイクルの運転中の運転能力と運転能力の予測値とも比較し、制御演算装置５７は判定部１０６が冷媒自然循環式サイクルの運転中の運転能力が小さいと判定したとき、冷媒自然循環式サイクルの運転状況を改善する運転状況改善部１０５を備えている。また、１０７は上記各検知部をまとめて示し、１０８は上記各操作部をまとめて示している。
【００２９】
制御演算装置５７は、判定部１０６が空調場熱負荷同定部１０１で検知された空調負荷より冷媒自然循環式サイクルの運転能力の予測値が大きいと判定したとき、冷媒自然循環式サイクル運転を起動し、この起動後に空調場熱負荷同定部１０１で空調負荷を再度検知するとともに、自然サイクル能力検知部１０３で冷媒自然循環式サイクル運転中の運転能力を検知し、前記再検出された空調負荷より大きいと判定部１０６で判定したとき、冷媒自然循環式サイクル運転を継続するように構成される。
【００３０】
また、制御演算装置５７は、前記判定部１０６が冷媒自然循環式サイクルの運転中の運転能力が予測値より大きいと判定したとき、前記空調場熱負荷同定部１０１で空調負荷を再度検知する。制御演算装置５７は、空調する部屋の熱容量並びに熱通過係数を少なくとも要因とする外部空調負荷に関わるパラメータ及び内部空調負荷を、複数の検知による計測量を元に時々刻々同定し、経年変化や状態変化に追随して切り替え運転する。さらに、制御演算装置５７は、空気温度条件と前記冷媒自然循環式サイクルの運転能力の関係を表す係数を、複数の検知による計測量を元に時々刻々同定し、経年変化や状態変化に追随して切り替え運転する。
【００３１】
次に、本実施例による空気調和機の動作について説明する。初めに、室外機３０は、室内機２０よりも高い位置に設置されているものとする。理由は、自然循環システムを稼動させるために、室外機３０から室内機２０に冷媒を流すために必要であり、近年の小型、中型ビルは、屋上に室外機を設置する場合が多く、室外機もそう言った使い方に合わせて設計されている場合が多いため、本施工に適合していると言える。
【００３２】
一般的に、空気調和機の能力（ヒートポンプの運転能力）は、圧縮機２１の運転周波数を変更する事によって変化し、能力範囲は機種によって違いはあるが、定格を１００％とすると３０％〜１２０％程度である。近年は、ヒートアイランド現象により、夏場に高い気温が発生する事が多く、また、空気調和機が施工される時点では、テナントの負荷が不明なため、空気調和機の容量選定としては、安全側を見越して、やや空調負荷に対して大きな容量を選定する傾向がある。そのため、夏場の酷暑時期以外は、空調負荷に対して空調能力が大きく、空気調和機容量の１０％〜２０％で運転される事が多くなっている。
【００３３】
空気調和機は、特にヒートポンプ（蒸気圧縮式ヒートオンプ）を使用するものは、５０％程度の能力出力時に最も効率が高くなるが、それ以下になると、圧縮機自体の効率低下や、ファン動力の占める割合が多くなる理由等により、効率が低下する。加えて、能力制御の不完全性により、不必要な空調能力を出力し、不要なサーモオフ、サーモオンが頻繁に行われる。この頻繁なサーモオフ、サーモオンによる断続運転は、空気調和機としては効率を低下させる。この圧縮機発停による効率低下は、例えばＪＩＳ  Ｂ８６１６にも効率低下係数ＣＤとして表現されている。この様に、空気調和機は、小さい負荷にて運転した場合、ハードウェア（ヒートポンプ自体の特性）による効率低下、ソフトウェア（運転制御）による効率低下があるため、大きく効率が低下する。
【００３４】
そこで、先ずは運転制御の改善による効率向上について記す。
【００３５】
運転制御の改善は、例えば特許第０３３８８３０５号にある様に、部屋の熱容量、熱通過率等のパラメータをシステム同定によって同定し、圧縮機周波数、能力、空調負荷の関係を明確にする事で、過度な能力を抑制して圧縮機の発停を防止し、不快感、効率低下を防止出来る。ここでは空調負荷として顕熱を例にとり、その要因を室外温度、熱容量と熱通過率を基本としたモデルで表現しているが、実際には、日射や侵入空気、人体発熱や機器発熱の様な内部発熱がある。また、直ちに熱負荷になるものと、時間遅れや振幅減少を伴って熱負荷になるものがあり、これらを一つ一つの熱負荷を分解して検知する事は困難である。
【００３６】
そこで、上記の通り単純なモデルで使用するが、このパラメータは時々刻々同定するため、値も時々刻々変化する。実際の負荷要因は、例えば日射によるものは、１時間単位でゆっくりと変化するが、パラメータが同定されるために要する時間は、２０〜３０ステップ程度なので、同定ステップを１０秒毎に行えば、３〜５分程度である。従って、同定パラメータは、実際の負荷の変動を追随する事が出来るので、同定値はゆっくり変動する事になる。ここで重要な事は、負荷要因の変動に合わせて、同定パラメータが追随し、変化する熱負荷環境に適応していると言う事である。従って、熱負荷が高い午前中にも、熱負荷が低い夜間でも、負荷パラメータに追随して、効率低下防止機能が働いている。
【００３７】
しかし、空調負荷が、空気調和機の最低能力を下回る場合、断続運転をせざるを得ない。これはハードウェアの改善を必要とする。
【００３８】
空気調和機の圧縮機の効率改善は必要であるが、それとは別に、近年、事務所やサーバルーム等、室外温度が室内温度よりも下がっても内部負荷が大きいため、空調負荷が発生する環境において、冷媒自然循環式サイクルを取り入れる事例が増えている。ここで冷媒自然循環式サイクルとは、室外機にて、室外温度で冷媒を凝縮、液化させ、室内機との高低差を利用して、液冷媒を室内機に搬送し、室内温度で冷媒を蒸発、気化(室内冷房の場合)させ、ガス化した冷媒を室外機まで搬送させる循環方式であり、高低差利用による搬送動力不要なシステムである。
【００３９】
冷媒自然循環式サイクルは搬送動力が不要なため、空気調和機の効率は大きく向上する。しかし欠点は、高低差が必要な事、室内温度が室外温度より低い場合で無いと能力を発揮しない事、また発生能力は、室内温度と室外温度の差に比例するため、空調負荷と逆の関係にある事である。
【００４０】
そのため、冷媒自然循環式サイクルに切り替える判断を間違えると、空調能力よりも空調負荷が大きい場合、室温が高くなって快適性に欠ける運転となったり、また室温上昇を防止させるために再度圧縮機を起動させる必要が出たりし、切り替えしない方が省エネルギとなる場合すらある。
【００４１】
そこで次の様な方法にて、無駄な切り替えを防止する。先ず、空調負荷の顕熱に関するモデルを下記数式（１）の通りとする。
【００４２】
【数１】
0005591214.tif 000002
【００４３】
ここでｔは時刻、Ｔｉ（ｔ）は室温、Ｔｏ（ｔ）は外気温、ＱａＳ（ｔ）は空気調和機による顕熱能力、ＬｉＳ（ｔ）は内部負荷による顕熱負荷、Ｃは室内家具や室内空気と壁表面を合わせた熱容量、Ｋは外気導入量、放熱量を合わせた熱通過係数を表す。空気調和機としては、圧縮機を用いる蒸気圧縮式ヒートポンプや自然循環式サイクルどちらでも成り立つ。
【００４４】
実際には、室内空気と内壁、室外空気と外壁が熱交換をし、壁内部は熱伝導により熱交換を行うため、数式（１）よりも高次の方程式となり、加えて空間的にも広がりがあるため、分布定数系による表現をしなければならない。しかし次数の低次元化による誤差、集中定数系による誤差は本発明の本質とは別問題であるので、ここでは数式（１）にて十分精度良くモデル化されている事を前提に話を進める。
【００４５】
また数式（１）は顕熱負荷だけであるため、潜熱負荷に関するモデルを作成し、それについて考慮するとより良い効果が得られる。但し顕熱負荷と同じ考え方であるため、ここでは省略する。
【００４６】
数式（１）を、一度数式（２）の通り書き直す。
【００４７】
【数２】
0005591214.tif 000003
【００４８】
数式（２）をサンプル時間Δｔにて離散化すると、下記の通りとなる。
【００４９】
【数３】
0005591214.tif 000004
【００５０】
【数４】
0005591214.tif 000005
【００５１】
ここでｋはサンプリングを表す番号である。係数ｃとｄはこのモデル（２）では同じ値であるが、敢えて別に扱う。次に、下記ベクトルを定義する。
【００５２】
【数５】
0005591214.tif 000006
【００５３】
【数６】
0005591214.tif 000007
【００５４】
ここで（・）’は転置を表し、
【００５５】
【数２４】
0005591214.tif 000008
【００５６】
はそれぞれａ〜ｄの推定値を表す。また推定ベクトル
【００５７】
【数２５】
0005591214.tif 000009
【００５８】
の初期値は既知、或いは０とする。この推定ベクトルは、下記数式に示す値
【００５９】
【数７】
0005591214.tif 000010
【００６０】
を最小にしたい場合、下記数式に示す逐次計算の形で求められる。
【００６１】
【数８】
0005591214.tif 000011
【００６２】
【数９】
0005591214.tif 000012
【００６３】
ここで、Ｐは共分散行列と呼ばれるもので、上記の例では４×４次元、初期値は既知、或いは０を代入する。
【００６４】
上記数式（７）を最小にする意味合いは、得られたデータが最もモデルである数式（３）に合致する様にパラメータを見出す事である。これはシステム同定であり、元々は観測雑音やシステム雑音が印加するシステムのモデルパラメータを同定する事を目的に発展した。
【００６５】
さてここで、全ての状態量は観測可能として扱っているが、空気調和機の顕熱能力ＱａＳはρＡを空気密度、Ｇを室内ファン風量、ｃｐを空気の比熱、ＴｉＥを室内機吹出温度とすると
【００６６】
【数１０】
0005591214.tif 000013
【００６７】
にて観測出来るが、室内内部顕熱負荷ＬｉＳ（ｋ）は観測出来ない。本質的に、システム同定は、未知パラメータを同定するためには、状態量は観測可能でなければならない。そこで、ここでは数式（２）と（３）で分かる通り、ｃとｄが同じ値である事を利用する。
【００６８】
先ず、ＬｉＳ（ｋ）の初期値として、常識的な初期値（物理的に取りうる値）を入力する。次に、そのＬｉＳ（ｋ）を用いて（４）〜（８）に沿ってシステム同定を行う。その際、正確では無いがａ〜ｄが求まるので、
【００６９】
【数１１】
0005591214.tif 000014
【００７０】
【数１２】
0005591214.tif 000015
【００７１】
として、室内内部顕熱負荷ＬｉＳ（ｋ）を求める事が出来る。以上より、パラメータａ〜ｄ及びＬｉＳ（ｋ）を推定する事が出来た。ａ、ｂが求まると、Ｃ、Ｋの値は、数式（４）を用いて下記の通り求まる。
【００７２】
【数１３】
0005591214.tif 000016
【００７３】
次に、空調顕熱負荷は、外部顕熱負荷Ｋ｛Ｔｉ（ｋ）−Ｔｏ（ｋ）｝と内部顕熱負荷
ＬｉＳ（ｋ）を足したものになるので、この値を
【００７４】
【数１４】
0005591214.tif 000017
【００７５】
とおく。ここでＬＳ（ｋ）は顕熱空調負荷である。
【００７６】
以上より、顕熱だけであるが、現在の空調負荷ＬＳ（ｋ）が推定出来た。以上の過程を潜熱負荷に対しても同様に行うが、潜熱負荷は外気と室内の絶対湿度差と換気量より決まる外部潜熱、人間の呼気等の内部潜熱を足し合わせたものとなり、基本的な考えは顕熱負荷と同じなので省略する。同じ手順で潜熱空調負荷ＬＬ（ｋ）を推定すると、全熱空調負荷Ｌ（ｋ）は
【００７７】
【数１５】
0005591214.tif 000018
【００７８】
にて計算出来る。
【００７９】
次に冷媒自然循環式サイクルについて検討する。冷媒自然循環式サイクルより出力可能な全熱空調能力をＱＮＡｃとし、ＱＮＡｃを外気温度、室内温度、冷媒自然循環式サイクルの配管長や室内外高低差に応じて予め計算で予測しておき、
【００８０】
【数１６】
0005591214.tif 000019
【００８１】
となった場合、つまり自然循環式サイクルにて空調負荷を賄えると判断出来る場合に、蒸気圧縮式ヒートポンプの圧縮機を停止させ、自然循環式サイクルに切り替える。また、切り替え後も、空調能力と室内温度の変化を数式（１）〜（１５）の手順で、現時刻の空調負荷を推定し続ける事が出来る。
【００８２】
次に、逆に蒸気圧縮式ヒートポンプに切り替える方法について説明する。
自然循環式サイクルは室内外ファンにより能力を調整出来るが、その範囲は小さく限られている。そこで室内温度が上昇すれば、その時点で能力不足と認識し、蒸気圧縮式ヒートポンプに切り替えると判断する手法もあるが、能力が不足する理由として、自然循環式サイクルとして発揮出来る能力が元々不足している場合と、本来は能力が発揮出来るにも拘らず、蒸発器に液冷媒が過剰に滞留したり、或いは冷媒量が不足している等、運転状態が悪いために能力不足になっている場合の二種類がある。そこで、先に説明した通り、自然循環方式の際にも空調負荷が推定出来るため、空調負荷推定値を用いる。
【００８３】
現在の空調能力は観測出来るため、それをＱＮＡ（ｋ）とし、熱交換器の仕様、外気温、室温、配管長、高低差、ファン風量等の条件によって予め予測された能力ＱＮＡｃ（ｋ）と比較する。その際、
【００８４】
【数１７】
0005591214.tif 000020
【００８５】
の場合は、現在の運転状態が悪い状態にあると判断し、その悪化条件を見つけて対処し、再度冷媒自然循環式サイクルを継続する。この悪化条件を見つけて対処する具体例としては、蒸発器に液冷媒が過剰に滞留したり、或いは冷媒量が不足している等、運転状態が悪いために能力不足になっている場合が多いので、蒸気圧縮式ヒートポンプを強制的に駆動して、液冷媒を循環させる。
【００８６】
【数１８】
0005591214.tif 000021
【００８７】
の場合は、問題無く運転されていると考える。
【００８８】
また推定空調負荷と現在発揮している能力が推定空調負荷Ｌ（ｋ）に対し
【００８９】
【数１９】
0005591214.tif 000022
【００９０】
であり、かつ室温をＴｉ（ｋ）、設定温度をＴｓｅｔとすると
【００９１】
【数２０】
0005591214.tif 000023
【００９２】
の場合は、現在室内温度が設定温度に近づいている最中か、或いは窓を開ける等、一時的に空調負荷が急変して室内温度が上昇している状態と考え、そのまま冷媒自然循環式サイクルを継続する。
【００９３】
但し、室温が許容値を超えた場合、
【００９４】
【数２１】
0005591214.tif 000024
【００９５】
の時は不快な状態を避けるため、冷媒自然循環式サイクルから蒸気圧縮式ヒートポンプに切り替える。
【００９６】
また
【００９７】
【数２２】
0005591214.tif 000025
【００９８】
となった場合は、空調負荷が大きいため、冷媒自然循環式サイクルでは設定温度に保持出来ないのは確実なので、蒸気圧縮式ヒートポンプに切り替える。
【００９９】
以上が蒸気圧縮式ヒートポンプから冷媒自然循環式サイクルに、また冷媒自然循環式サイクルから蒸気圧縮式ヒートポンプに切り替えるアルゴリズムである。
【０１００】
実施例１は、空調負荷を全熱負荷として記述したが、室内温度のみを制御する場合、顕熱能力と顕熱負荷だけを扱っても良い。図３Ａ、図３Ｂは、実施例１を示すフローチャートであり、以下動作を説明する。
【０１０１】
ステップ１で運転が開始されるとステップ２で空調機が起動し、ステップ３で。蒸気圧縮式ヒートポンプ２１が起動する。ステップ４でヒートポンプ能力検知部１０２によってヒートポンプの全熱能力が検知され、次のステップ５で空調場熱負荷同定部１０１で空調する場所の全熱負荷Ｌが同定される。
【０１０２】
次のステップ６では、前記の全熱負荷Ｌ（ｋ）より冷媒自然循環式サイクル運転の予め予測された全熱能力ＱＮＡｃ（ｋ）（予測値）が大きいか否かを判定部１０６によって判定される（数式１６参照）。全熱能力ＱＮＡｃ（ｋ）は、熱交換器の仕様、外気温、室温、配管長、高低差、ファン風量等の条件によって、予め自然サイクル能力予測部１０４で予測される。判定の結果、不成立（Ｆ）であればステップ４に戻り、成立（Ｔ）であればステップ７でヒートポンプ２１を停止し、ステップ８で冷媒自然循環式サイクルの運転を起動する。
【０１０３】
次のステップ９では、運転中の冷媒自然循環式サイクルの全熱能力ＱＮＡが自然サイクル能力検知部１０３で検知される。続いて、ステップ１０では、上記で検知された全熱能力ＱＮＡが予測された全熱能力ＱＮＡｃより小さいか判定される（前記数式１７参照）。Ｑｔは誤差許容値である。判定の結果、成立（Ｔ）であれば、冷媒自然循環式サイクルの運転状態が悪いために全熱能力ＱＮＡが不足していると判断して、ステップ１８で自然循環サイクル運転状況の改善動作を行う。具体的には前述のように蒸気圧縮式ヒートポンプを強制的に駆動して、液冷媒を循環させる。判定の結果、不成立（Ｆ）であれば、全熱能力ＱＮＡが十分大きいので（数式１８参照）、問題無く冷媒自然循環式サイクルが運転されていると判断し、ステップ１１で再度空調場の全熱負荷Ｌを同定する。
【０１０４】
ステップ１２では、上記で得られた全熱負荷Ｌとステップ９で検知された全熱能力ＱＮＡとが判定部１０６で比較され、全熱能力ＱＮＡが全熱負荷Ｌより大きいか否かが判定される。判定の結果、全熱能力ＱＮＡが全熱負荷Ｌより大きいとき（Ｔ）、ステップ１３に移り、小さいとき（Ｆ）ステップ１７に移る。ステップ１７では自然循環サイクルの運転を停止し、ステップ３に戻ってヒートポンプを起動する。
【０１０５】
ステップ１３では、空調されている室温Ｔｉと設定温度Ｔｓｅｔを比較し、室温Ｔｉが設定温度に第１許容値Ｔｔ１を加えた値を超えたか否か判定する。室温Ｔｉが許容値を超えたとき（前記数式２１参照）、ステップ１７に移って、自然循環サイクルの運転を停止し、ステップ３に戻ってヒートポンプを起動する。室温Ｔｉが許容値を超えてないとき、さらにステップ１４で設定温度Ｔｓｅｔを超えたか否か判定される。判定の結果、超えていれば（Ｔ）ステップ９に戻り、超えていなければ（Ｆ）ステップ１５に移る。
【０１０６】
ステップ１５では、室温Ｔｉが設定温度から第２許容値Ｔｔ２を引いた値より小さいか否を判定する。判定の結果、室温Ｔｉが小さい場合（Ｔ）、室温Ｔｉが十分低いと判断してステップ１６で空調運転を停止する。判定の結果、室温Ｔｉが大きい場合（Ｆ）はステップ９に戻り、冷媒自然循環式サイクルの運転を継続しつつ、冷媒自然循環式サイクルの全熱能力ＱＮＡが自然サイクル能力検知部１０３で検知される。
【０１０７】
図３Ａ、図３Ｂのフローチャートでは、蒸気圧縮式ヒートポンプ運転から冷媒自然循環式サイクルの運転に切り替える際に、冷媒自然循環式サイクル起動の前後のステップ５と１１で空調する場所の全熱負荷Ｌを求めることで、時々刻々と変化する空調負荷を常に正確に求めている。また、上記の各全熱負荷は、全熱能力との比較直前に求められるので、常に新しい全熱負荷が求められる。そして、上記で得た各全熱負荷Ｌと、冷媒自然循環式サイクルの全熱能力の予測値と現在値とをそれぞれ比較して、冷媒自然循環式サイクルの全熱能力で全熱負荷Ｌを賄えると判断したときに、冷媒自然循環式サイクルを起動および運転継続している。
【０１０８】
したがって、冷媒自然循環式サイクルの起動後に同サイクルの全熱能力の不足により、圧縮機の運転に切り換わることが少なくなり、圧縮機の発停の頻発を防止することができる。
【０１０９】
また、図３Ａ、図３Ｂのフローチャートでは、冷媒自然循環式サイクル起動後のステップ９で、運転中の冷媒自然循環式サイクルの全熱能力を検知し、ステップ１０で、上記で検知された全熱能力が予め予測された全熱能力より小さいか判定され、小さい場合はステップ１８で、冷媒自然循環式サイクル運転の改善動作を行い、再度、ステップ９で運転中の冷媒自然循環式サイクルの全熱能力を検知して上記動作を繰返す。
【０１１０】
したがって、冷媒自然循環式サイクルの起動後に同サイクルの運転中の実際の全熱能力を把握し、これに基いて運転継続を決めているので、冷媒自然循環式サイクルの能力不足で運転に切り換わることが少なくなり、圧縮機の発停の頻発を防止することができる。また、冷媒自然循環式サイクルの起動後に同サイクルの全熱能力の不足が判明しても、運転の改善動作により運転中の全熱能力を回復することができ、圧縮機の運転に切り換わることが少なくなる。
【０１１１】
図４Ａ、図４Ｂは、上記で説明した実施例１の全熱負荷を顕熱負荷に置き換えて示した実施例２のフローチャートである。動作は、図３Ａ、図３Ｂのフローチャートと同じであるため、説明を省略する。
【０１１２】
図５〜図９は、空調負荷を表すパラメータを推定している計算結果である。図５で、８１は空調全熱内部負荷Ｌｉ（ｋ）で、８２は空調全熱内部負荷推定値Ｌｉを示し、ここでは、内部負荷量Ｌｉは、周期３６００ステップのゆっくりとした正弦波と、確率分布として正規分布を示す白色雑音を加えた変動をしていると仮定した。図５〜図９の通り、各パラメータａ〜ｄの推定値は比較的安定して推定出来ており、内部負荷Ｌｉは時間が掛かっているが、±５％以内の推定が出来ている。８３はパラメータａ推定値、８４はパラメータｂ推定値、８５はパラメータｃ推定値、８６はパラメータｄ推定値である。
【０１１３】
最後に、年間の省電力、効率を計算する。ＪＩＳ  Ｂ８６１６の指標とは別に、年間効率を擬似的ではあるが、より正確に表す下記の値Ｅを次の数式で定義する。
【０１１４】
【数２３】
0005591214.tif 000026
【０１１５】
ここでＱ（Ｔｊ）は室外温度Ｔｊにおける空気調和機の発生する能力、Ｐ（Ｔｊ）空気調和機の消費電力、ｎｊはＴｊが発生する発生時間である。分子は年間に発生する能力の総和であり、分母は年間に発生する消費電力の総和である。これは、ある室外温度Ｔｊが発生した際、空調負荷が一意に決まり、それと同じ能力と消費電力を空気調和機が発生させ、室内温度を設定温度に保って安定している、と言う事を前提にしている。実際の空調負荷は、前述した通り、室外温度だけに依存せず、日射やその他の影響を受け、また時々刻々変化していくので、Ｅは必ずしも実使用上の空気調和機の年間効率を表す指標とはならないが、室外温度Ｔｊ毎に能力と消費電力を測定し、指標に入れ込む事で、ＪＩＳ  Ｂ８６１６より、より正確に年間効率を表していると言える。
【０１１６】
数式（２３）を元に、従来機と、本発明の空気調和機の年間効率Ｅを計算にて比較したところ、従来の空気調和機に対して、本発明の空気調和機は４％向上する事が確認された。この効果は冷媒自然循環式サイクルに依るものであるが、外気温度発生率、空調期間はＪＩＳ  Ｂ８６１６に準拠しており、東京の気候であるため、更に気温が低くて内部負荷が高い空調場では、効果が大きくなると期待される。上記の例で、発生能力が同じと仮定すると、年間の消費電力も従来機に比べて、本発明の空気調和機は４％低減出来る。
【符号の説明】
【０１１７】
１…運転開始ステップ、２…空調機起動ステップ、３…蒸気圧縮式ヒートポンプ起動ステップ、４…蒸気圧縮式ヒートポンプ全熱能力検知ステップ、５…空調場全熱負荷同定ステップ、６…自然循環式サイクル予測全熱能力対空調全熱負荷判定ステップ、７…蒸気圧縮式ヒートポンプ停止ステップ、８…自然循環式サイクル起動ステップ、９…自然循環式サイクル全熱能力検知ステップ、１０…自然循環式サイクル運転状況判定ステップ、１１…空調全熱負荷同定ステップ、１２…自然循環式サイクル全熱能力対空調全熱負荷判定ステップ、１３…快適性判定ステップ、１４…室温状態判定ステップ、１５…サーモオフ判定ステップ、１６…空調機サーモオフ実施ステップ、１７…自然循環式サイクル停止ステップ、１８…自然循環式サイクル運転状況改善ステップ、２０…室外機、２１…圧縮機、２２１…第一の室外熱交換器、２２２…第二の室外熱交換器２３…室外ファン、２４…アキュムレータ、２５…四方弁、２６…受容器、２７…室外膨張弁、２８…液側配管、２９…ガス側配管、３０…室内機、３１…室内熱交換器、３２…室内ファン、３３…室内膨張弁、３４…インバータ圧縮機周波数操作器、３５…室外ファン送風能力操作器、３６…室外膨張弁開度操作器、３７…四方弁操作器、３８…室内ファン送風能力操作器、３９…室内膨張弁開度操作器、４０…圧縮機吸入温度検知器、４１…圧縮機吐出温度検知器、４２…過冷却熱交換器出口温度検知器、４３…室外熱交換器出口温度検知器、４４…室外温度検知器、４５…圧縮機吸入圧力検知器、４６…圧縮機吐出圧力検知器、４７…室外湿度検知器、４８…室内機吸込温度検知器、４９…室内機吹出温度検知器、５０…室内機液管温度検知器、５１…室内機ガス温度検知器、５２…室内機吸込湿度検知器、５３…室内機吹出湿度検知器、５４…第一の開閉弁、５５…第二の開閉弁、５６…室内温度設定器、５７…制御演算装置、６１…運転開始ステップ、６２…空調機起動ステップ、６３…蒸気圧縮式ヒートポンプ起動ステップ、６４…蒸気圧縮式ヒートポンプ顕熱能力検知ステップ、６５…空調場顕熱負荷同定ステップ、６６…自然循環式サイクル予測顕熱能力対空調顕熱負荷判定ステップ、６７…蒸気圧縮式ヒートポンプ停止ステップ、６８…自然循環式サイクル起動ステップ、６９…自然循環式サイクル顕熱能力検知ステップ、７０…自然循環式サイクル運転状況判定ステップ、７１…空調顕熱負荷同定ステップ、７２…自然循環式サイクル顕熱能力対空調顕熱負荷判定ステップ、７３…快適性判定ステップ、７４…室温状態判定ステップ、７５…サーモオフ判定ステップ、７６…空調機サーモオフ実施ステップ、７７…自然循環式サイクル停止ステップ、７８…自然循環式サイクル運転状況改善ステップ、８１…空調全熱内部負荷、８２…空調全熱内部負荷推定値、８３…パラーメータａ推定値、８４…パラメータｂ推定値、８５…パラメータｃ推定値、８６…パラメータｄ推定値、１００…制御部、１０１…空調場熱負荷同定部、１０２…ヒートポンプ能力検知部、１０３…自然サイクル能力検知部、１０４…自然サイクル能力予測部、１０５…運転状況改善部、１０６…判定部。
【図１】
0005591214.tif 000027
【図２】
0005591214.tif 000028
【図３Ａ】
0005591214.tif 000029
【図３Ｂ】
0005591214.tif 000030
【図４Ａ】
0005591214.tif 000031
【図４Ｂ】
0005591214.tif 000032
【図５】
0005591214.tif 000033
【図６】
0005591214.tif 000034
【図７】
0005591214.tif 000035
【図８】
0005591214.tif 000036
【図９】
0005591214.tif 000037