(57)【要約】      （修正有）【課題】エネルギー消費量の基準値からの差分を機器ごとに算出する省エネルギー要素抽出装置を提供する。【解決手段】本発明の省エネルギー要素抽出装置１は、系統に属する要素のエネルギー消費量の実測値を要素ごとに取得し、系統に関する所定の前提値を入力としエネルギー消費量の基準値を出力とする基準モデル３４から、基準値を要素ごとに取得する基準値推定部２１と、取得した実測値と取得した基準値との差分を要素ごとに算出する解析部２２とを備える。【選択図】図２
2019086999.tif 000002
【特許請求の範囲】
【請求項１】
系統に属する要素のエネルギー消費量の実測値を前記要素ごとに取得し、
前記系統に関する所定の前提値を入力とし前記エネルギー消費量の基準値を出力とする基準モデルから、前記基準値を前記要素ごとに取得する基準値推定部と、
前記取得した実測値と前記取得した基準値との差分を前記要素ごとに算出する解析部と、
を備えることを特徴とする省エネルギー要素抽出装置。
【請求項２】
前記基準値推定部は、
前記実測値を前記要素ごとに取得できない場合、前記系統のエネルギー消費量の実測値を取得し、
前記取得した系統のエネルギー消費量の実測値を、所定の方法で前記要素のそれぞれに按分し、
前記解析部は、
前記按分した実測値と前記取得した基準値との差分を前記要素ごとに算出すること、
を特徴とする請求項１に記載の省エネルギー要素抽出装置。
【請求項３】
前記基準値推定部は、
前記実測値を前記要素ごとに取得できない場合、前記系統のエネルギー消費量の実測値を取得し、
前記要素のそれぞれのエネルギー消費量を増減させる事象の組合せが発生する場合における前記要素のそれぞれのエネルギー消費量の基準値を前記基準モデルから前記組合せごとに取得し、
当該取得した要素のそれぞれのエネルギー消費量の基準値を合計することによって前記系統のエネルギー消費量の基準値の候補を前記組合せごとに算出し、
前記解析部は、
前記基準値の候補から、前記取得した系統のエネルギー消費量の実測値に所定の基準を満たす程度に近似するものを抽出し、
前記抽出した候補を前記組合せが発生する確率の順に並び替えること、
を特徴とする請求項２に記載の省エネルギー要素抽出装置。
【請求項４】
前記解析部は、
前記基準モデルを使用しない場合は、前記系統又は前記要素のエネルギー消費量の過去の実測値と現時点の実測値との差分を、前記要素ごとに算出すること、
を特徴とする請求項３に記載の省エネルギー要素抽出装置。
【請求項５】
前記解析部は、
同種の要素が複数存在する場合、同種の要素のうち、前記エネルギー消費量が所定の基準を満たす程度に大きいものを抽出すること、
を特徴とする請求項４に記載の省エネルギー要素抽出装置。
【請求項６】
前記実測値、前記基準値及び前記差分を画面表示する表示処理部を備えること、
を特徴とする請求項５に記載の省エネルギー要素抽出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、省エネルギー要素抽出装置に関する。
【背景技術】
【０００２】
例えばオフィスビルには大規模な空調設備が存在する。空調設備は、冷凍機、冷却塔、空調機、冷却水ポンプ、冷水ポンプ等の機器、及び、それらの間を接続する配管等を備える。そして、これらの機器の動力源は、電力であることが多い。例えば冷凍機の圧縮機の動力は、電動機であり、冷却塔のファンの動力も電動機である。これらの機器は、例えば室温がある設定値になるように制御される。これらの機器の電動機は、当然のことながら電気エネルギーを消費する。空調設備の規模が大きくなると、消費する電気エネルギー量も膨大になる。すると、単に室内の快適さの観点だけではなく、電気エネルギーの節約、地球環境の維持等の観点からも、空調設備が適正に制御されることが必要になる。
【０００３】
特許文献１の建物省エネルギー評価監視装置は、建物内に設置された設備のエネルギー消費量を実測する一方、予測モデルを使用して、設備が理想的に制御された場合のエネルギー消費量を予測する。そして、当該建物省エネルギー評価監視装置は、実測されたエネルギー消費量を予測されたエネルギー使用量と比較することによって、設備の制御状態の不良を発見し、その原因を追及する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００３−２１６７１５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
実際の空調設備の制御においては、多く存在する機器のうち、どの機器がエネルギーを必要以上に多く消費しているかを知ることが肝要である。そして、それを知って初めて、その機器の制御を見直し、エネルギー消費量を削減することが可能になる。しかしながら、特許文献１の建物省エネルギー評価監視装置は、多くの機器を一体として評価・監視するに留まり、エネルギーを必要以上に多く消費している機器を機器群のうちから抽出するという発想を欠いている。
【０００６】
特許文献１の建物省エネルギー評価監視装置を、単体の機器に使用することも可能である。しかしながら、エネルギー消費量を実測するセンサ（電力量計等）が、単体の機器に対して設置されておらず、複数の機器を含む系統に対して設置されている場合も多い。このような場合、特許文献１の建物省エネルギー評価監視装置を使用しても、エネルギーを必要以上に多く消費している機器を発見することはできない。
そこで、本発明は、エネルギー消費量の基準値からの差分を機器ごとに算出することを目的とする。
【課題を解決するための手段】
【０００７】
本発明の省エネルギー要素抽出装置は、系統に属する要素のエネルギー消費量の実測値を要素ごとに取得し、系統に関する所定の前提値を入力としエネルギー消費量の基準値を出力とする基準モデルから、基準値を要素ごとに取得する基準値推定部と、取得した実測値と取得した基準値との差分を要素ごとに算出する解析部と、を備えることを特徴とする。
その他の手段については、発明を実施するための形態のなかで説明する。
【発明の効果】
【０００８】
本発明によれば、エネルギー消費量の基準値からの差分を機器ごとに算出することができる。
【図面の簡単な説明】
【０００９】
【図１】本実施形態の環境を説明する図である。
【図２】省エネルギー要素抽出装置等の構成を説明する図である。
【図３】機器特性情報を説明する図である。
【図４】（ａ）及び（ｂ）は、運転履歴情報を説明する図である。
【図５】シミュレーション結果情報を説明する図である。
【図６】基準モデルを説明する図である。
【図７】エネルギー増減事象情報を説明する図である。
【図８】第１の処理手順のフローチャートである。
【図９】第２の処理手順のフローチャートである。
【図１０】第３の処理手順のフローチャートである。
【図１１】（ａ）及び（ｂ）は、電力量計の配置を説明する図である。
【図１２】（ａ）及び（ｂ）は、電力量計の配置を説明する図である。
【図１３】（ａ）は、理想基準値及び実測値の表示例である。（ｂ）は、理想基準値と実測値との機器ごとの差分の表示例である。
【図１４】絞り込み前の劣化基準値の表示例である。
【図１５】絞り込み後の劣化基準値の表示例である。
【発明を実施するための形態】
【００１０】
以降、本発明を実施するための形態（“本実施形態”という）を、図等を参照しながら詳細に説明する。本実施形態は、電力及び都市ガスを動力源とする複数の機器から構成される空調系統を例にしている。しかしながら、動力源は、圧縮空気等の他の動力源であってもよい。さらに、複数の機器は、空調系統以外の系統、例えば生産・用役ライン等を構成する機器であってもよい。さらに、ここでの機器は、それ自身に対して消費エネルギーが定義できる任意のもの（“要素”と呼ぶ）を代表している。要素は、物体であってもよいし、生体であってもよい。つまり、本発明は、相互に関連しながら個別にエネルギーを消費する複数の要素を有する系統に対して適用可能である。
【００１１】
（本実施形態の環境）
図１に沿って、本実施形態の環境を説明する。あるオフィスビルの空調設備５１は、冷却塔５２、冷凍機５３、冷却水ポンプ５４、冷水一次ポンプ５５、冷水二次ポンプ５６、冷水往ヘッダ５７、冷水二次ヘッダ５８、冷水還ヘッダ５９及び空調機６０を有する。例えば夏季において、室内６１内に通じる空調機６０は、人間等から発生する熱を空気を介して吸収する。空調機６０は、吸収した熱を冷水に与える。冷水は、空調機６０から冷水還ヘッダ５９に流れる。この流れの動力となっているのは、冷水二次ポンプ５６である。冷水一次ポンプ５５は、冷水を冷水還ヘッダ５９から冷凍機５３に吐出する。
【００１２】
冷凍機５３は、自身が有する圧縮機（図示せず）による冷媒の状態変化（気体／液体）を利用して、冷水の熱を冷却水に与える。冷却水は、冷凍機５３から冷却塔５２に流れる。この流れの動力となっているのは、冷却水ポンプ５４である。冷却塔５２は、ファンで発生させた風（外気）を冷却水と熱交換することにより、冷却水中の熱を大気中に逃がす。冷却水ポンプ５４は、熱を逃がした冷却水を冷凍機５３に吐出する。冷凍機５３において熱を冷媒に与えた冷水は、冷水往ヘッダ５７に流れる。この流れの動力となっているのは、冷水一次ポンプ５５である。冷水二次ポンプ５６は、冷水を、冷水二次ヘッダ５８を経由して空調機６０に吐出する。なお、図１における矢印は、配管及びその内側を流れる冷水及び冷却水の方向を示す。空調機６０と室内６１との間の双方向矢印は、往ダクト及び還ダクトを示す。
【００１３】
配管及びダクトの要所（例えば、冷却塔５２の冷却水出口、冷凍機５３の冷水出口等）に、温度計６２、流量計６３、圧力計６４及び制御弁６５が設置されている。また、温度計（外気温度計）６２は、室内６１及び外気に触れる位置にも設置されている。湿度計(外気湿度計)６６は、外気に触れる位置に設置されている。さらに、冷却塔５２、冷凍機５３、冷却水ポンプ５４、冷水一次ポンプ５５及び冷水二次ポンプ５６に対して動力としての電力及び／又は都市ガスが供給されている。これらの機器の近傍には、当該機器が消費する電力及び／又は都市ガスの量を実測する電力量計及び／又はガス量計（図示せず）が設置されている。電力量計及びガス量計の位置及び個数については詳しく後記する。
【００１４】
（省エネルギー要素抽出装置等の構成）
図２に沿って、省エネルギー要素抽出装置１等の構成を説明する。省エネルギー要素抽出装置１は、一般的なコンピュータであり、オフィスビルの設備監視室等に設置される。省エネルギー要素抽出装置１は、中央制御装置１１、入力装置１２、出力装置１３、主記憶装置１４、補助記憶装置１５及び通信装置１６を備える。これらは、相互にバスで接続されている。補助記憶装置１５は、機器特性情報３１、運転履歴情報３２、シミュレーション結果情報３３、基準モデル３４及びエネルギー増減事象情報３５を格納する。
【００１５】
主記憶装置１４における、基準値推定部２１、解析部２２及び表示処理部２３は、プログラムである。以降の説明において、“○○部は”と動作主体を記した場合、それは、中央制御装置１１が補助記憶装置１５から○○部を読み出し、主記憶装置１４にロードしたうえで○○部の機能（詳細後記）を実現することを意味する。省エネルギー要素抽出装置１は、ネットワーク４を介して、１又は複数の現地監視装置２と接続され得る。
【００１６】
機器３は、例えば、図１における冷却塔５２、冷凍機５３、冷却水ポンプ５４、冷水一次ポンプ５５、冷水二次ポンプ５６等である。請求項における“要素”には、ここでの機器３が相当する。センサ４１は、図１における温度計６２、流量計６３、圧力計６４、制御弁６５、電力量計及びガス量計のうち、機器３に付属するものである。センサ４１は、付属する機器から、温度、流量、圧力、弁開度、消費電力量、消費ガス量等を実測する。そしてセンサ４１は、実測したデータを現地監視装置２に送信する。１つの機器３に対して複数のセンサ４１が対応していてもよい。センサ４１ｂは、室内に設置されている温度計６２、並びに、建物の外部に設置されている温度計６２及び湿度計６６である。センサ４１ｂもまた、外気の温度等を実測し、実測したデータを現地監視装置２に送信する。
【００１７】
現地監視装置２もまた、一般的なコンピュータであり、例えば、機器３の近傍に設置される。現地監視装置２は、中央制御装置、入力装置（図示せず）、出力装置４２、主記憶装置４３、補助記憶装置４４及び通信装置４５を備える。これらは、相互にバスで接続されている。主記憶装置４３におけるデータ取得部４６は、センサ４１及び４１ｂからセンサ値を取得し、省エネルギー要素抽出装置１に送信する処理の動作主体となるプログラムである。補助記憶装置４４には運転履歴情報３２が記憶されている。この運転履歴情報３２は、省エネルギー要素抽出装置１の補助記憶装置１５に記憶されている運転履歴情報３２のうち、実測されたデータの送信元である機器３に関する部分である。本実施形態では、現地監視装置２及び機器３は、１対１に対応している。しかしながら、１つの現地監視装置２に対して複数の機器３が対応していてもよい。
【００１８】
（機器特性情報）
図３に沿って、機器特性情報３１を説明する。機器特性情報３１においては、機器ＩＤ（Identifier）欄１０１に記憶された機器ＩＤに関連付けて、機器種類欄１０２には機器種類が、系統欄１０３には系統が、諸元１欄１０４から諸元８欄１１１までには各諸元値が、上流圧力損失係数欄１１２には上流圧力損失係数が、定格電力欄１１３には定格電力が記憶されている。
【００１９】
機器ＩＤ欄１０１の機器ＩＤは、機器３を一意に特定する識別子である。
機器種類欄１０２の機器種類は、機器３の種類であり、ここでは、“冷凍機”、“冷却塔”、“空調機”、“冷水ポンプ”（冷水一次ポンプ及び冷水二次ポンプを含む）、並びに、“冷却水ポンプ”のうちの何れかである。
系統欄１０３の系統は、ある機器が他の機器とともに属する機器の群である。図１において熱の遷移サイクルを説明したが、同じ系統に属する機器は、同じ熱の遷移サイクルに属している。
【００２０】
諸元１欄１０４から諸元８欄１１１までの各諸元値は、機器の能力、特性を示す設計上の数値である。各諸元値が具体的に何を示すかは、機器種類によって異なる。例えばポンプの諸元値は、内径、最大流量、Ｐ−Ｑ特性等である。冷却塔の諸元値は、設計上の冷却能力、風量、冷却水往還温度、外気湿球温度、エンタルピー総括熱伝達係数等である。
上流圧力損失係数欄１１２の上流圧力損失係数は、当該機器の直ぐ上流における配管又はダクトの圧力損失係数である。上流圧力損失係数に代替して下流圧力損失係数が記憶されてもよい。
定格電力欄１１３の定格電力は、機器の定格電力である。
各欄の“＃”は、異なる値を省略的に示している。“−”は、その欄が空欄（該当なし）であることを示している。
【００２１】
（運転履歴情報）
図４（ａ）及び図４（ｂ）に沿って、運転履歴情報３２を説明する。図４（ａ）の運転履歴情報３２においては、機器ＩＤ欄１２１に記憶された機器ＩＤに関連付けて、機器種類欄１２２には機器種類が、取得時点欄１２３には取得時点が、一次エネルギー消費量欄１２４には一次エネルギー消費量が、運転状態欄１２５には運転状態値が、負荷条件欄１２６には負荷条件値が、異常欄１２７には異常フラグが、保守欄１２８には保守フラグが記憶されている。
【００２２】
機器ＩＤ欄１２１の機器ＩＤは、図３の機器ＩＤと同じである。
機器種類欄１２２の機器種類は、図３の機器種類と同じである。
取得時点欄１２３の取得時点は、当該レコード（行）のデータが取得された時点の年月日時分秒である。
【００２３】
一次エネルギー消費量欄１２４の一次エネルギー消費量は、機器３が消費する一次エネルギーの消費量である。機器３の動力源としては、電力、都市ガス等様々なものが存在する。これらのそれぞれの消費量にそれぞれの換算係数を乗算した結果を合計したものが一次エネルギー消費量である。一次エネルギー消費量は、センサ４１が実測した値に基づいて省エネルギー算出装置１によって算出される。なお、ここでの各レコードの一次エネルギー消費量は、取得時点における瞬間値であってもよいし、取得時点を含む期間における積算値であってもよい。一次エネルギー消費量が瞬間値である場合は、その瞬間値は、取得時点の前後の複数の時点における瞬間値の平均値であってもよい。
【００２４】
運転状態欄１２５の運転状態値は、機器の運転状態を示す値である。運転状態値は、機器種類によって異なる。例えば、冷凍機の運転状態値の１つは、冷水出口温度１２５ａであってもよい。冷却塔の運転状態値の１つは、冷却水出口温度１２５ｃであってもよい。運転状態値は、センサ４１が実測した値である。なお、ここでの各レコードの運転状態値は、取得時点における瞬間値であってもよいし、取得時点の前後の複数の時点における瞬間値の平均値であってもよい。
【００２５】
負荷条件欄１２６の負荷条件値は、機器がその能力を発揮する前提となる負荷を算出するための値であり、ここでは、外気温度１２６ａ、外気湿度１２６ｂ及び室内温度１２６ｃである。一般に、室内温度は、空調機に付属するセンサ４１によって実測される。室内温度は、人間が快適に活動できる目標値に制御される。そして、図４（ａ）の例では、外気温度及び外気湿度も、空調機に付属する（空調機の近辺の建物外壁に設置される）センサ４１によって実測される。
【００２６】
異常欄１２７の異常フラグは、事後的に運転状態値を分析した結果、機器に異常が認められたことを示すフラグである。ここでの異常フラグは、その異常が一次エネルギー消費量を正常状態に比して増加させることを示す“増加”、又は、その異常が一次エネルギー消費量を正常状態に比して減少させることを示す“減少”の何れかである。機器が劣化すると、一般的にはその一次エネルギー消費量は増加する。しかしながら、機器の種類及び制御方法によっては、機器が劣化すると、その一次エネルギー消費量が減少する場合もある。例えば、ポンプの回転速度が一定に維持されている場合、ポンプの劣化や配管の詰り等の圧力損失上昇により、ポンプの一次エネルギー消費量は減少する。
【００２７】
保守欄１２８の保守フラグは、事後的に保守員が機器を保守したことを示す“○”である。
各欄の“＃”、“♭”及び“＃♭”は、異なる値を省略的に示している。“−”は、その欄が空欄（該当なし）であることを示している。なお、“＃”、“♭”及び“＃♭”の区別については後記する。
【００２８】
図４（ｂ）もまた運転履歴情報３２である。図４（ａ）においては一次エネルギー消費量が機器ごとに記憶されている一方、図４（ｂ）においては、一次エネルギー消費量が機器ごとには記憶されておらず、機器が属する系統全体についての数値として記憶されている。なぜならば、１個の電力量計等が１個の系統に対して設置されているからである。本実施形態の省エネルギー要素抽出装置１は、まさにこのような場合に適している。
【００２９】
（シミュレーション結果情報）
図５に沿って、シミュレーション結果情報３３を説明する。図５のシミュレーション結果情報３３の構成は、図４（ａ）の運転履歴情報３２の構成から、異常欄１２７及び保守欄１２８を削除したものである。省エネルギー要素抽出装置１は、基準モデル３４（図２）を使用して機器の一次エネルギー消費量を推定する。より具体的には、省エネルギー要素抽出装置１は、基準モデル３４に対して前提条件を入力する。すると、基準モデル３４は、一次エネルギー消費量を機器ごとに出力する。つまり、省エネルギー要素抽出装置１は、基準モデル３４に入力する前提条件を様々に変化させることによって、一次エネルギー消費量を機器ごとにシミュレーションし、その結果をシミュレーション結果情報３３として出力することができる。なお、本実施形態において、“基準値”という語は、基準モデル３４が出力する一次エネルギー消費量を意味する。
【００３０】
（基準モデル）
図６に沿って、基準モデル３４を説明する。基準モデル３４は、複数の数式（等式、不等式、条件式等）の集合である。個々の数式は、複数の変数を含む。基準モデル３４の一部の変数に具体的な値（前提値）が入力（代入）されると、基準モデル３４は、他の変数を未知数とする連立方程式を解き、解を出力する。入力値及び出力値の例は、以下の通りである。
【００３１】
[入力値]
・機器の諸元値（図３の欄１０４〜欄１１１）
・機器の上流圧力損失係数（図３の欄１１２）
・機器の定格電力（図３の欄１１３）
・外気温度、外気湿度及び室内温度（図４（ａ）又は図４（ｂ）の欄１２６ａ〜欄１２６ｃ）
[出力値]
・機器の一次エネルギー消費量（図５の欄１２４）
[入力値又は出力値の何れかになる値]
・運転状態値（図４（ａ）又は図４（ｂ）の欄１２５ａ〜欄１２５ｆ）
【００３２】
図３、図４（ａ）、図４（ｂ）及び図５において、“＃”は、基準モデル３４に対する入力値である。“♭”は基準モデル３４からの出力値である。“＃♭”は、入力値又は出力値の何れかになる値である。例えば、運転状態値（欄１２５）のうち、冷水出口温度１２５ａ及び冷水入口温度１２５ｂが入力値であるとする。このとき、運転状態値のうち、冷却水出口温度１２５ｃ、冷却水入口温度１２５ｄ、送風温度１２５ｅ及び冷水流量１２５ｆが出力値となる。
【００３３】
（エネルギー増減事象の組合せ）
ある機器が劣化すると、その機器は、理想的な一次エネルギー消費量で稼働できなくなる。前記のように、劣化の進行に伴い、一次エネルギー消費量が理想的な基準値に比して増加するか減少するかは、制御方法、機器の種類、機器の諸元等により異なる。ここでは、単純化のために、冷凍機及び冷却塔が劣化すると、それらの一次エネルギー消費量は増加し、冷水ポンプ及び冷却水ポンプが劣化すると、それらの一次エネルギー消費量は減少するものとする。
【００３４】
いま、“冷凍機の一次エネルギー消費量の増加”を短縮して“冷凍機増”と表現する。同様に、“冷却塔の一次エネルギー消費量の増加”、“冷却水ポンプの一次エネルギー消費量の減少”及び “冷水ポンプの一次エネルギー消費量の減少”をそれぞれ“冷却塔増”、“冷却水ポンプ減”及び“冷水ポンプ減”と表現する。そして、“冷凍機増”、“冷却塔増”、“冷却水ポンプ減”及び“冷水ポンプ減”のそれぞれをエネルギー増減事象と呼ぶ。
【００３５】
（発生確率）
機器が属する系統が劣化すると、前記した４つのエネルギー増減事象のうちの１つが単独で発生し、又は、４つのエネルギー増減事象のうちの複数が同時に発生する。これらの組合せの数は、４Ｃ１＋４Ｃ２＋４Ｃ３＋４Ｃ４＝１４である。そして、１４個の組合せのそれぞれに対して発生確率が定義され得る。省エネルギー要素抽出装置１は、例えば以下の算出方法１〜４の何れかで発生確率を算出することができる。
【００３６】
〈算出方法１〉保守員による保守結果を利用する。
例えば、１０日間の評価対象期間のうち、“冷凍機増”が発生していた期間が“２日”であった場合、発生確率は、２／１０＝０.２である。“冷凍機増”及び“冷却塔増”が同時に発生していた期間が“１日”であった場合、発生確率は、１／１０＝０.１である
〈算出方法２〉メーカが作成した設計値を利用する。
機器のメーカは、機器の不具合が発生する確率を公開している場合がある。
【００３７】
〈算出方法３〉ベイズの定理を適用して条件付確率を算出する。
例えば、いま、エネルギー増減事象“Ａ”が発生する確率Ｐ（Ａ）、及び、エネルギー増減事象“Ａ”が発生したときエネルギー増減事象“Ｂ”が発生する条件付確率Ｐ（Ｂ｜Ａ）が既知であるとする。すると、Ｐ（Ａ∩Ｂ）＝Ｐ（Ａ）×Ｐ（Ｂ｜Ａ）が成立するので、Ｐ（Ａ∩Ｂ）も既知となる。Ｐ（Ａ∩Ｂ）は、エネルギー増減事象“Ａ”及び“Ｂ”が同時に発生する確率である。
〈算出方法４〉ベイジアンネットワークを学習する。
算出方法３は、１つのエネルギー増減事象が、相互の因果関係を有することなく独立に発生することを前提としている。例えば“冷凍機増”が発生すると必ず“冷却塔増”が発生する（両者は独立ではない）ことがわかっているとする。この場合、もし教師データがあればベイジアンネットワークを学習することによって１４個の発生確率を算出してもよい。
【００３８】
（エネルギー増減事象情報）
図７に沿って、エネルギー増減事象情報３５を説明する。エネルギー増減事象情報３５においては、組合せＩＤ欄１３１に記憶された組合せＩＤに関連付けて、エネルギー増減事象の組合せ欄１３２にはエネルギー増減事象の組合せが、発生確率欄１３３には発生確率が記憶されている。
【００３９】
組合せＩＤ欄１３１の組合せＩＤは、前記したエネルギー増減事象の組合せを一意に特定する識別子である。
エネルギー増減事象の組合せ欄１３２のエネルギー増減事象の組合せは、組合せを構成するエネルギー増減事象１３２ａ〜１３２ｄを示す“○”の組合せである。
発生確率欄１３３の発生確率は、前記したエネルギー増減事象の組合せの発生確率である。“＃＃”は、異なる値（０以上１以下）を省略的に示している。
【００４０】
（エネルギー増減事象の一般化）
例えば“冷凍機増”というエネルギー増減事象の粒度をさらに細かくしていくと、以下のようなエネルギー増減事象を得ることができる。
・特定の冷凍機の一次エネルギー消費量が増加すること
・特定の冷凍機の一次エネルギー消費量が特定の値だけ増加すること
・特定の冷凍機の一次エネルギー消費量が特定の範囲に入ること
・特定の冷凍機が特定の症状を示すこと
【００４１】
結局、エネルギー増減事象は、“任意の機器が任意の異常状態にあること”を指す。そして、エネルギー増減事象のそれぞれに対し、その発生確率が定義され得る。また。エネルギー増減事象のそれぞれを引き起こす機器の諸元値（劣化諸元値）が決定され得る。さらに、基準モデルは、劣化諸元値を入力値として受け付けると、機器が劣化した場合を想定した一次エネルギー消費量を機器ごとに出力できる。
【００４２】
（処理手順）
図８〜図１０は、処理手順のフローチャートである。フローチャートは３種類存在し、それぞれ第１の処理手順、第２の処理手順、及び、第３の処理手順と呼ばれる。これら３つの処理手順は、それぞれ別々に実行され得るが、処理手順ごとに前提が異なる。なお、以降では、電力量計及びガス量計を代表して単に“電力量計”という。
【００４３】
〈第１の処理手順における前提〉
・１個の機器に対して１個の電力量計（Ｗｈ）が配置されている（図１１（ａ）参照）。・基本モデル３４は、各機器が理想的な状態で運転される場合の基準値（“理想基準値”と呼ぶ）を機器ごとに出力する。
【００４４】
〈第２の処理手順における前提〉
・１個の系統に対して１個の電力量計（Ｗｈ）が配置されている（図１１（ｂ）参照）。なお、このように電力量計が配置されている場合が、“実測値を要素ごとに取得できない場合”に相当する。
・省エネルギー要素抽出装置１は、所定の方法（詳細後記）で、系統全体が消費する一次エネルギー消費量を当該系統に属する機器のそれぞれに按分することができる。
・基本モデル３４は、理想基準値を機器ごとに出力する。
【００４５】
〈第３の処理手順における前提〉
・１個の系統に対して１個の電力量計（Ｗｈ）が配置されている（図１１（ｂ）参照）。
・基本モデル３４は、各機器が劣化した状態で運転される場合の基準値（“劣化基準値”と呼ぶ）を機器ごとに出力する。
【００４６】
以降では、図８〜図１０に沿って、第１の処理手順、第２の処理手順及び第３の処理手順をこの順に説明する。その過程で、図１１〜図１５を適宜参照する。
【００４７】
（第１の処理手順）
図８は、第１の処理手順のフローチャートである。
ステップＳ２０１において、省エネルギー要素抽出装置１の基準値推定部２１は、診断対象期間及び診断対象系統を受け付ける。具体的には、基準値推定部２１は、ユーザが入力装置１２を介して、診断対象期間の始期の年月日及び終期の年月日、並びに、診断対象系統を入力するのを受け付ける。このとき、基準値推定部２１は、ユーザの入力を受け付けるまでもなく、直近の１週間、直近の１月等予め指定された期間が診断対象期間として入力されたと看做してもよい。
【００４８】
ステップＳ２０２において、基準値推定部２１は、運転履歴情報３２（図４（ａ））を取得する。具体的には、基準値推定部２１は、運転履歴情報３２から、診断対象系統に属する各機器の診断対象期間のレコード（行）をすべて取得する。
【００４９】
ステップＳ２０３において、基準値推定部２１は、一次エネルギー消費量の実測値を機器ごとに取得する。具体的には、基準値推定部２１は、ステップＳ２０２において取得したレコードから一次エネルギー消費量を取得する。このとき、基準値推定部２１は、実測値として、１日ごとの平均値を算出してもよいし、診断対象期間の積算値（合計値）を算出してもよい。基準値推定部２１は、１日ごとの平均値を算出する場合、最大値及び最小値を“はずれ値”として無視したうえで、平均値を算出してもよいし、１日の値を前後の所定の日数の移動平均値としてもよい。
【００５０】
ステップＳ２０４において、基準値推定部２１は、基準モデル３４を読み出す。具体的には、基準値推定部２１は、補助記憶装置１５から基準モデル３４を読み出す。
【００５１】
ステップＳ２０５において、基準値推定部２１は、一次エネルギー消費量の理想基準値を機器ごとに算出する。具体的には、第１に、基準値推定部２１は、ステップＳ２０４において読み出した基準モデル３４に対して以下の値を入力（代入）する。
・機器の諸元値（図３の欄１０４〜欄１１１）
・機器の上流圧力損失係数（図３の欄１１２）
・機器の定格電力（図３の欄１１３）
・外気温度、外気湿度及び室内温度（図４（ａ）の欄１２６ａ〜欄１２６ｃ）
・運転状態値（図４（ａ）の欄１２５ａ〜欄１２５ｆ）のうち予めユーザが入力値として指定したもの
なお、ここで入力される機器の諸元値及び機器の上流圧力損失係数は、機器が劣化していない理想状態にある場合の値である。
【００５２】
第２に、基準値推定部２１は、基準モデル３４が以下の値を機器ごとに出力するのを受け付ける。
・機器の一次エネルギー消費量（理想基準値、図５の欄１２４）
・運転状態値（図５の欄１２５ａ〜欄１２５ｆ）のうち予めユーザが出力値として指定したもの
なお、基準値推定部２１は、実測値にあわせ、理想基準値として、１日ごとの平均値を算出してもよいし、診断対象期間の積算値（合計値）を算出してもよい。
第３に、基準値推定部２１は、シミュレーション結果情報３３（図５）を作成し、補助記憶装置１５に記憶する。
【００５３】
ステップＳ２０６において、省エネルギー要素抽出装置１の解析部２２は、一次エネルギー消費量の実測値と理想基準値との差分を算出する。具体的には、解析部２２は、ステップＳ２０３において取得した一次エネルギー消費量の実測値から、ステップＳ２０５の“第２”において受け付けた一次エネルギー消費量（理想基準値）を減算した差分を機器ごとに算出する。
【００５４】
ステップＳ２０７において、省エネルギー要素抽出装置１の表示処理部２３は、実測値、理想基準値及び差分を機器ごとに表示する。具体的には、第１に、表示処理部２３は、出力装置１３及び／又は現地監視装置２の出力装置４２に図１３（ａ）のグラフを表示する。図１３（ａ）の左側の棒グラフは、ステップＳ２０５の“第２”において受け付けた一次エネルギー消費量（理想基準値）を機器（冷凍機、冷却塔、冷水ポンプ及び冷却水ポンプ）ごとに積み上げたものである。図１３（ａ）の右側の棒グラフは、ステップＳ２０３において取得した一次エネルギー消費量の実測値を機器ごとに積み上げたものである。
【００５５】
第２に、表示処理部２３は、出力装置１２及び／又は現地監視装置２の出力装置４２に図１３（ｂ）のグラフを表示する。図１３（ｂ）の棒グラフは、ステップＳ２０６において算出した機器ごとの差分を、その値が大きい順に左から並べたものである。なお、冷水ポンプの差分は“０”であるので、その棒グラフは現れていない。
【００５６】
ステップＳ２０８において、省エネルギー要素抽出装置１の解析部２２は、同種の機器のうち、一次エネルギー消費量の実測値が最大であるものを抽出する。例えば、冷水ポンプが複数存在する場合、解析部２２は、ステップＳ２０３において取得した冷水ポンプごとの一次エネルギー消費量の実測値のうち、最大のものを抽出する。そして、解析部２２は、その冷水ポンプの機器ＩＤ、一次エネルギー消費量の実測値、理想基準値、及び、差分を出力装置１３に表示する。より一般的には、解析部２２は、所定の基準を満たす程度にその一次エネルギー消費量の実測値が大きい機器を抽出してもよい。例えば解析部２２は、“一次エネルギー消費量の実測値が大きい順に○個の機器”、“一次エネルギー消費量の実測値が所定の数値以上である機器”等を抽出してもよい。
その後、第１の処理手順を終了する。
【００５７】
（第１の処理手順の変形例）
前記の例では、基準値推定部２１は、ステップＳ２０５において基準モデル３４を使用している。しかしながら、基準値推定部２１は、所定の入力値を縦横の軸とするテーブルデータ（マトリクス）を使用して、理想基準値を取得してもよい。このようなテーブルデータは機器ごとに補助記憶装置１５に記憶され、例えば、縦軸の目盛は、外気比エンタルピー又は外気湿球温度であり、横軸の目盛は、負荷である。マトリクスの交点のセルには、理想基準値（出力値）が記憶されている。
【００５８】
（第２の処理手順）
図９は、第２の処理手順のフローチャートである。図９のステップＳ２０１及びＳ２０２の処理内容は、図８のステップＳ２０１及びＳ２０２の処理内容と同じである。省エネルギー要素抽出装置１の基準値推定部２１は、ステップＳ２０２の直後にステップＳ３０１に進む。
【００５９】
ステップＳ３０１において、基準値推定部２１は、一次エネルギー消費量の実測値を系統全体として取得する。具体的には、基準値推定部２１は、ステップＳ２０２において取得したレコードから一次エネルギー消費量を取得する。このとき取得される一次エネルギー消費量は、系統全体の一次エネルギー消費量である（図４（ｂ）及び図１１（ｂ）参照）。なお前記と同様、基準値推定部２１は、実測値として、１日ごとの平均値を算出してもよいし、診断対象期間の積算値（合計値）を算出してもよい。
【００６０】
ステップＳ３０２において、基準値推定部２１は、実測値を機器ごとに按分する。具体的には、基準値推定部２１は、ステップＳ３０１において取得した一次エネルギー消費量の実測値を、以下の按分方法１〜３の何れかで機器ごとに按分する。
〈按分方法１〉
・各機器の定格電力（図３の欄１１３）で按分する。
【００６１】
〈按分方法２〉
・系統に属するすべての機器を稼働したうえで、系統全体の一次エネルギー消費量の実測値を取得する。
・ある機器を停止したうえで、系統全体の一次エネルギー消費量の実測値を取得する。
・当該機器の停止前後における系統全体の一次エネルギー消費量の実測値の差分を、当該機器の一次エネルギー消費量の実測値とする。
・各機器の実測値の合計が系統全体の実測値の合計に一致しない場合は、各機器の実測値で、系統全体の実測値を按分する。
なお、電力量計（Ｗｈ）の配置が図１２（ａ）又は図１２（ｂ）のようになっている場合であっても按分方法２は有効である。
【００６２】
〈按分方法３〉
・系統に属するすべての機器を一旦停止する。
・ある機器を稼働したうえで、系統全体の一次エネルギー消費量の実測値を取得する。
・当該機器の稼働前後における系統全体の一次エネルギー消費量の実測値の差分を、当該機器の一次エネルギー消費量の実測値とする。
・各機器の実測値の合計が系統全体の実測値の合計に一致しない場合は、各機器の実測値で、系統全体の実測値を按分する。
なお、電力量計（Ｗｈ）の配置が図１２（ａ）又は図１２（ｂ）のようになっている場合であっても按分方法３は有効である。
【００６３】
その後、基準値推定部２１は、ステップＳ２０４に進む。図９のステップＳ２０４〜Ｓ２０８の処理内容は、図８のステップＳ２０４〜Ｓ２０８の処理内容と同じである。
【００６４】
（第３の処理手順）
図１０は、第３の処理手順のフローチャートである。図１０のステップＳ２０１及びＳ２０２の処理内容は、図８のステップＳ２０１及びＳ２０２の処理内容と同じである。図１０のステップＳ３０１の処理内容は、図９のステップ３０１の処理内容と同じである。図１０のステップＳ２０４の処理内容は、図８のステップＳ２０４の処理内容と同じである。省エネルギー要素抽出装置１の基準値推定部２１は、ステップＳ２０４の直後にステップＳ４０１に進む。
【００６５】
ステップＳ４０１において、基準値推定部２１は、エネルギー増減事象の組合せを作成する。具体的には、第１に、基準値推定部２１は、前記したエネルギー増減事象を作成する。ここで作成されたエネルギー増減事象は、“冷凍機増”、“冷却塔増”、“冷却水ポンプ減”及び“冷水ポンプ減”であるとする。
第２に、基準値推定部２１は、エネルギー増減事象の組合せを作成する。この例では、図７に示すように１４個の組合せが作成されることになる。
【００６６】
ステップＳ４０２において、基準値推定部２１は、エネルギー増減事象の組合せの発生確率を算出する。具体的には、基準値推定部２１は、前記した算出方法１〜４の何れかを使用して、発生確率をエネルギー増減事象の組合せごとに算出する。
【００６７】
ステップＳ４０３において、基準値推定部２１は、エネルギー増減事象情報３５（図７）を作成する。具体的には、第１に、基準値推定部２１は、エネルギー増減事象情報３５のレコードを組合せの数（この例では“１４”）だけ作成する。
第２に、基準値推定部２１は、作成したレコードのエネルギー増減事象の組合せ欄１３２に、ステップＳ４０１の“第２”において作成した組合せを構成するエネルギー増減事象の“○”を記憶する。
第３に、基準値推定部２１は、発生確率欄１３３に、ステップＳ４０２において算出した発生確率を記憶する。
第４に、基準値推定部２１は、組合せＩＤを採番したうえで、組合せＩＤ欄１３１に記憶する。
【００６８】
ステップＳ４０４において、基準値推定部２１は、一次エネルギー消費量の劣化基準値の候補を機器ごとに算出する。具体的には、第１に、基準値推定部２１は、ステップＳ２０４において読み出した基準モデル３４に対して以下の値を入力（代入）する。
・機器の諸元値（図３の欄１０４〜欄１１１）
・機器の上流圧力損失係数（図３の欄１１２）
・機器の定格電力（図３の欄１１３）
・外気温度、外気湿度及び室内温度（図４（ｂ）の欄１２６ａ〜欄１２６ｃ）
・運転状態値（図４（ｂ）の欄１２５ａ〜欄１２５ｆ）のうち予めユーザが入力値として指定したもの
【００６９】
なお、ここで入力される機器の諸元値及び機器の上流圧力損失係数は、機器が劣化状態にある場合の値である。例えば“冷凍機増”というエネルギー増減事象が発生している場合、冷凍機の諸元値又は上流圧力損失係数は、劣化していない場合に比して変化している。また、“冷凍機増”及び“冷却塔増”という２個のエネルギー増減事象が同時に発生している場合、冷凍機及び冷却塔の、諸元値又は上流圧力損失係数は、劣化していない場合に比して変化している。そこで、基準値推定部２１は、この変化した値の組合せを１４通り取得し、基準モデル３４に入力することになる。
【００７０】
第２に、基準値推定部２１は、基準モデル３４が以下の値を機器ごとに出力するのを受け付ける。
・機器の一次エネルギー消費量（劣化基準値、図５の欄１２４）
・運転状態値（図５の欄１２５ａ〜欄１２５ｆ）のうち予めユーザが出力値として指定したもの
なお、基準値推定部２１は、実測値にあわせ、劣化基準値として、１日ごとの平均値を算出してもよいし、診断対象期間の積算値（合計値）を算出してもよい。基準値推定部２１は、１４通りの入力値に対応して、基準モデル３４から１４通りの基準値の出力を受け付けることとなる。１４通りの基準値のそれぞれは、機器ごとの劣化基準値からなる。そして、基準値推定部２１は、シミュレーション結果情報３３（図５）を１４通り作成し、補助記憶装置１５に記憶する。
【００７１】
第３に、基準値推定部２１は、機器ごとの劣化基準値を合計することによって、系統の一次エネルギー消費量の劣化基準値を算出する。この段階で、基準値推定部２１は、系統全体の一次エネルギー消費量の劣化基準値を１４通り保持していることになる。ここで保持される１４通りの劣化基準値を“候補”と呼ぶ。なぜならば、エネルギー消費量の増減事象の組合せの数が大きくなる（例えば数百）場合、ユーザはこれらの“候補”を絞り込む必要があるからである。
【００７２】
なお、ステップＳ４０４の“第３”において、基準値推定部２１は、絞り込む前の劣化基準値の候補を出力装置１３に表示してもよい。図１４は、その表示例である。一番左側の棒グラフは、理想基準値である。それ以外の１４本の棒グラフは、１４通りの劣化基準値の“候補”であり、図７の組合せＩＤ“Ｃ０１〜Ｃ１４”に対応している。何れの棒グラフも、冷凍機、冷却塔、冷水ポンプ及び冷却水ポンプの一次エネルギー消費量を縦に積み上げたものである。
【００７３】
ステップＳ４０５において、省エネルギー要素抽出装置１の解析部２２は、劣化基準値の候補を絞り込む。具体的には、解析部２２は、ステップＳ４０４の“第３”において算出した系統全体の劣化基準値のうちから、ステップＳ３０１において取得した系統全体の一次エネルギー消費量の実測値に所定の基準を満たす程度に近似しているものを抽出する。例えば、系統全体の一次エネルギー消費量の実測値が“１０５”であったとする。そして、所定の基準が、“実測値に対して±２.５の範囲内”であったとする。このとき、解析部２２は、系統全体の一次エネルギー消費量の劣化基準値の候補のうちから、“１０２.５以上１０７.５以下”のものを抽出する。
【００７４】
ステップＳ４０６において、解析部２２は、絞り込んだ劣化基準値を発生確率の大きい順に並び替える。具体的には、解析部２２は、ステップＳ４０５において絞り込んだ候補を、エネルギー増減事象情報３５（図７）の発生確率が大きい順に並べ替える。解析部２２は、絞り込んだ候補を、発生確率が小さい順に並べ替えてもよい。
【００７５】
ステップＳ４０７において、省エネルギー要素抽出装置１の表示処理部２３は、並び替えた劣化基準値を表示する。具体的には、表示処理部２３は、ステップＳ４０６において並べ替えた劣化基準値を出力装置１３に表示する。図１５は、その表示例である。一番左側の棒グラフは、理想基準値である。それ以外の５本の棒グラフは、絞り込まれた後の５通りの劣化基準値である。これらの５本の棒グラフは、エネルギー増減事象の組合せが発生する発生確率の大きい順に左から並んでいる。
【００７６】
ステップＳ４０８において、表示処理部２３は、並び替えた劣化基準値と理想基準値との差分を機器ごとに表示する。具体的には、第１に、表示処理部２３は、並び替えた劣化基準値のうちの任意の１つをユーザが選択するのを受け付ける。いま、ユーザは、左から２番目の“冷凍機増”に対応する劣化基準値を選択したとする。
第２に、表示処理部２３は、“冷凍機増”に対応する劣化基準値と理想基準値との差分を機器ごとに出力装置１３に表示する。図１３（ｂ）は、その表示例である。
【００７７】
ステップＳ４０９において、省エネルギー要素抽出装置１の解析部２２は、同種の機器が複数存在する場合、同種の機器のうち、一次エネルギー消費量の劣化基準値が最大であるものを抽出する。例えば、冷水ポンプが３個存在するとする。このとき、基準値推定部２１は、ステップＳ４０４の“第２”において、冷水ポンプの一次エネルギー消費量の劣化基準値を、３×１４＝４２通り算出していることになる。解析部２２は、これらのうち、最大のものを抽出する。そして、省エネルギー要素抽出装置１の表示処理部２３は、その冷水ポンプの機器ＩＤ、エネルギー増減事象の組合せ、一次エネルギー消費量の劣化基準値、理想基準値、及び、差分を出力装置１３に表示する。
【００７８】
より一般的には、解析部２２は、所定の基準を満たす程度にその一次エネルギー消費量の劣化基準値が大きい機器を抽出してもよい。例えば解析部２２は、“一次エネルギー消費量の劣化基準値が大きい順に○個の機器”、“一次エネルギー消費量の劣化基準値が所定の数値以上である機器”等を抽出してもよい。
その後、第３の処理手順を終了する。
【００７９】
（変形例）
省エネルギー要素抽出装置１は、官公庁等が運営する外部機器から取得した気象データを使用して、外気温度及び外気湿度を取得してもよい。
省エネルギー要素抽出装置１は、常設のセンサ４１に加えて又は代替して、臨時に設置されたセンサからデータを取得してもよい。
冷却水及び冷水は、ブライン等の不凍剤を含んでいてもよい。
系統は、冷水二次ポンプを有さない冷水一次ポンプ方式であってもよい。
省エネルギー要素抽出装置１は、現地監視装置２と一体の構成であってもよい。
【００８０】
省エネルギー要素抽出装置１は、現地監視装置２の出力装置４２にデータを表示する場合、当該現地監視装置２に対応する機器３以外の機器についてのデータを表示しなくてもよい。
省エネルギー要素抽出装置１の解析部２２は、一次エネルギー消費量の実測値と、基準モデル３４が出力した基準値との差分を算出することに代えて、過去の一次エネルギー消費量の実測値と、現時点の一次エネルギー消費量の実測値との差分を算出してもよい。このようにすれば、省エネルギー要素抽出装置１は、基準モデル３４を使用することなく、複数の時点間について、一次エネルギー消費量の実測値の増減を知ることができる。
【００８１】
省エネルギー要素抽出装置１は、一次エネルギー消費量に代えて、以下のような指標を使用することも可能である。
・通年エネルギー消費効率（ＡＰＦ：Annual Performance Factor）
・定格電力／定格消費電力比（ＣＯＰ：Coefficient of Performance）
・搬出熱量／ポンプ消費電力比（ＷＴＦ：Water Transportation Factor）
【００８２】
（診断後の機器制御）
本実施形態において一次エネルギー消費量が増加していると診断された機器（要注意機器）は、その後、診断結果に応じた具体的な制御の対象となる。要注意機器は、例えば台数制御のローテーションから除外されることが好ましい。また、ユーザは、優先的に要注意機器の制御を見直す、又は、保守を行うことが望ましい。
【００８３】
（本実施形態の効果）
本実施形態の省エネルギー要素抽出装置の効果は以下の通りである。
（１）一次エネルギー消費量の実測値と基準値との差分を機器ごとに知ることができる。
（２）実測値が機器ごとに実測できない場合であっても、機器ごとの実測値の目安を知ることができる。
（３）系統全体の一次エネルギー消費量を実測値に近似させるようなエネルギー増減事象の組合せを、発生確率の大きい順に知ることができる。
（４）基準モデルを使用することなく、一次エネルギー消費量を時系列で比較できる。
（５）同種の機器が複数存在する場合、最も早く対策を講じるべき機器を特定できる。
（６）一次エネルギー消費量の実測値と基準値との差分を機器ごとの容易に視認できる。
【００８４】
なお、本発明は前記した実施例に限定されるものではなく、様々な変形例が含まれる。例えば、前記した実施例は、本発明を分かり易く説明するために詳細に説明したものであり、必ずしも説明したすべての構成を備えるものに限定されるものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることも可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能である。
【符号の説明】
【００８５】
１      省エネルギー要素抽出装置
２      現地監視装置
３      機器
４      ネットワーク
１１    中央制御装置
１２    入力装置
１３    出力装置
１４    主記憶装置
１５    補助記憶装置
１６    通信装置
２１    基準値推定部
２２    解析部
２３    表示処理部
３１    機器特性情報
３２    運転履歴情報
３３    シミュレーション結果情報
３４    基準モデル
３５    エネルギー増減事象情報
４１、４１ｂ  センサ
４２    出力装置
４３    主記憶装置
４４    補助記憶装置
４５    通信装置
４６    データ取得部
【図１】
2019086999.tif 000003
【図２】
2019086999.tif 000004
【図３】
2019086999.tif 000005
【図４】
2019086999.tif 000006
【図５】
2019086999.tif 000007
【図６】
2019086999.tif 000008
【図７】
2019086999.tif 000009
【図８】
2019086999.tif 000010
【図９】
2019086999.tif 000011
【図１０】
2019086999.tif 000012
【図１１】
2019086999.tif 000013
【図１２】
2019086999.tif 000014
【図１３】
2019086999.tif 000015
【図１４】
2019086999.tif 000016
【図１５】
2019086999.tif 000017