(57)【要約】【課題】簡易な方法で空調システム内の各空調機器の使用エネルギー量を推定することができる空調システム内機器の使用エネルギー量推定方法、及び、この方法が実装された使用エネルギー量推定装置を提供する。【解決手段】使用エネルギー量推定装置１０で実行される使用エネルギー量推定方法は、所定の期間における、空調システムの総電力量と、空調機器毎の設定温度、当該空調機器が制御する空間の室温、及び、当該空調機器の稼働状況と、空調以外の機器の稼働状況と、からなる運転データを測定する運転データ測定ステップと、空調機器毎の設定温度、室温及び稼働状況、並びに、空調以外の機器の稼働状況を入力とし、総電力量を出力とするモデルを作成するモデル作成ステップと、このモデルに基づいて、空調機器毎の設定温度、室温及び稼働状況から、当該空調機器毎の使用電力量の推定値を算出する推定値算出ステップと、を有する。【選択図】図３
2012017930.tif 000002
【特許請求の範囲】
【請求項１】
空調機器と、空調以外の機器とを有する空調システムにおいて、前記空調機器毎の使用エネルギー量を推定する使用エネルギー量推定方法であって、
所定の期間における、前記空調システムの総エネルギー量と、前記空調機器毎の設定温度、当該空調機器が制御する空間の室温、及び、当該空調機器の稼働状況と、前記空調以外の機器の稼働状況と、からなる運転データを測定する運転データ測定ステップと、
前記空調機器毎の前記設定温度、前記室温及び前記稼働状況、並びに、前記空調以外の機器の前記稼働状況を入力とし、前記総エネルギー量を出力とするモデルを作成するモデル作成ステップと、
前記モデルに基づいて、前記空調機器毎の前記設定温度、前記室温及び前記稼働状況から、当該空調機器毎の使用エネルギー量の推定値を算出する推定値算出ステップと、を有することを特徴する空調システム内機器の使用エネルギー量推定方法。
【請求項２】
前記空調機器は、室外機及び室内機から構成され、
前記モデル作成ステップで作成される前記モデルは、前記総エネルギー量をＷとし、前記空調システムを構成するｍ台の前記空調機器のうち、第ｉ番目の前記空調機器の室内機の前記稼働状況をＯ１i、前記設定温度をＳi、前記室温をＴiとし、前記空調システムを構成するｎ組の前記空調以外の機器のうち、第ｊ番目の前記空調以外の機器の稼働状況をＯ２jとし、モデル係数をα１i，α２j，βiとしたとき、次式
Ｗ  ＝  Σα１i・Ｏ１i  ＋  Σα２j・Ｏ２j  ＋  Σβi・（Ｔi−Ｓi）・Ｏ１i
但し、ｉ＝１…ｍ，ｊ＝１…ｎ
で表され、
前記推定値算出ステップで算出される前記空調機器毎の使用エネルギー量の推定値は、第ｓ番目から第ｔ番目の前記室内機が接続された前記室外機の推定エネルギー量の推定値をＷｐとしたとき、次式
Ｗｐ  ＝  Σ（α１k・Ｏ１k  ＋  βk・（Ｔk−Ｓk）・Ｏ１k）
但し、ｋ＝ｓ…ｔ
で表されることを特徴とする請求項１に記載の空調システム内機器の使用エネルギー量推定方法。
【請求項３】
前記運転データ測定ステップにおいて、前記所定の期間内における、前記空調機器及び前記空調以外の機器の運転状況を、適宜変化させることを特徴とする請求項１または２に記載の空調システム内機器の使用エネルギー量推定方法。
【請求項４】
前記空調機器が、１台の室外機に２台以上の室内機が接続された構成のときに、
前記運転データ測定ステップにおいて、前記室内機毎に、前記設定温度、前記室温及び前記稼働状況を測定し、且つ、前記モデル作成ステップにおいて、前記当該室内機毎の前記設定温度、前記室温及び前記稼働状況を入力とする前記モデルを作成することにより、前記推定値算出ステップにおいて、当該空調機器における前記室内機毎の前記使用エネルギー量の推定値に対する寄与度を算出するように構成されたことを特徴とする請求項１〜３いずれか一項に記載の空調システム内機器の使用エネルギー量推定方法。
【請求項５】
前記空調機器に運転モードが設けられているときに、当該運転モード毎に前記運転データ測定ステップ、前記モデル作成ステップ及び前記推定値算出ステップを実行することを特徴する請求項１〜４いずれか一項に記載の空調システム内機器の使用エネルギー量推定方法。
【請求項６】
前記空調システムで測定された前記総エネルギー量の実測値をＷとし、
前記モデル作成ステップで作成された前記空調システムのモデルにより求められる前記総エネルギー量の推定値をＷＴとし、前記モデルにより求められる前記空調機器毎の前記使用エネルギー量の推定値をＷｐとしたとき、
前記推定値算出ステップは、補正された前記空調機器毎の前記使用エネルギー量の推定値Ｗｐ′を、次式
Ｗｐ′  ＝  Ｗｐ・Ｗ／ＷＴ
により算出することを特徴とする請求項１〜５いずれか一項に記載の空調システム内機器の使用エネルギー量推定方法。
【請求項７】
前記空調システムを構成する機器の定義情報が記憶されたデータ関係定義情報ファイルと、
前記運転データが記憶された運転データファイルと、
前記空調機器毎の使用エネルギー量の推定値が記憶される推定結果ファイルと、
前記データ関係情報ファイルから前記定義情報を読み出し、前記運転データファイルから前記運転データを読み出して、請求項１〜６いずれか一項に記載の使用エネルギー量推定方法により、前記空調システムに含まれる前記空調機器毎の前記使用エネルギー量の推定値を算出して前記推定結果ファイルに記憶させるモデル作成・推定演算部と、を有することを特徴する使用エネルギー量推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、空調システム内機器の使用エネルギー量推定方法、及び、この方法が実装された使用エネルギー量推定装置に関する。
【背景技術】
【０００２】
店舗やオフィス等の施設において、省エネを目的とした使用エネルギー、例えば電力の管理を行う場合には、その内部に設置されている消費電力の大きい多数の機器に関して、個別に使用電力量を把握する必要がある。しかし、このような電力量を計測するためには、電力量計測機器を設置するためのコストがかかることなどから、個別の電力量計測ではなく、配電盤で全体の電力量計測が行われたり、あるいは、系統別の電力量計測が行われたりしているのが一般的である。このため、特許文献１にあるように総電力量を各機器の定格電力と運転時間とを用いて推定する方式や、特許文献２にあるように総電力量を各機器の最低必要電力と運転時間とを用いて按分する方式や、特許文献３にあるように空調室外機に電力量測定装置を付加する方式が提案されている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００４−１９４４２９号公報
【特許文献２】特開平５−１９６６４０号公報
【特許文献３】国際公開第ＷＯ２００７／０３２０６５号
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、特許文献３に開示されている方式では、空調機器に電力量計測装置を付加する必要があるため低コストでの個別機器の電力量把握が難しいという課題があった。また、特許文献１または２に開示されている方式では、空調機器で消費される電力量が、その空調機器での設定温度と制御対象温度との差によって変動することが考慮されておらず、定格出力や最低必要電力量を考慮しただけであるので、正確に個別の電力量が推定できないという課題があった。
【０００５】
本発明は、このような課題に鑑みてなされたものであり、簡易な方法で空調システム内の各空調機器の使用エネルギー量を推定することができる空調システム内機器の使用エネルギー量推定方法、及び、この方法が実装された使用エネルギー量推定装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
前記課題を解決するために、本発明に係る空調システム内機器の使用エネルギー量推定方法は、空調機器と、空調以外の機器とを有する空調システムにおいて、空調機器毎の使用エネルギー量を推定する使用エネルギー量推定方法であって、所定の期間における、空調システムの総エネルギー量と、空調機器毎の設定温度、当該空調機器が制御する空間の室温、及び、当該空調機器の稼働状況と、空調以外の機器の稼働状況と、からなる運転データを測定する運転データ測定ステップと、空調機器毎の設定温度、室温及び稼働状況、並びに、空調以外の機器の稼働状況を入力とし、総エネルギー量を出力とするモデルを作成するモデル作成ステップと、このモデルに基づいて、空調機器毎の設定温度、室温及び稼働状況から、当該空調機器毎の使用エネルギー量の推定値を算出する推定値算出ステップと、を有することを特徴する。
【０００７】
このような本発明に係る空調システム内機器の使用エネルギー量推定方法において、空調機器は、室外機及び室内機から構成され、モデル作成ステップで作成されるモデルは、例えば、総エネルギー量をＷとし、空調システムを構成するｍ台の空調機器のうち、第ｉ番目の空調機器の室内機の稼働状況をＯ１i、設定温度をＳi、室温をＴiとし、空調システムを構成するｎ組の空調以外の機器のうち、第ｊ番目の空調以外の機器の稼働状況をＯ２jとし、モデル係数をα１i，α２j，βiとしたとき、次式
Ｗ＝Σα１i・Ｏ１i  ＋  Σα２j・Ｏ２j  ＋  Σβi・（Ｔi−Ｓi）・Ｏ１i
但し、ｉ＝１…ｍ，ｊ＝１…ｎ
で表され、さらに、推定値算出ステップで算出される空調機器毎の使用エネルギー量の推定値は、第ｓ番目から第ｔ番目の室内機が接続された室外機の推定エネルギー量の推定値をＷｐとしたとき、次式
Ｗｐ＝Σ（α１k・Ｏ１k  ＋  βk・（Ｔk−Ｓk）・Ｏ１k）
但し、ｋ＝ｓ…ｔ
で表すことができる。
【０００８】
また、このような本発明に係る空調システム内機器の使用エネルギー量推定方法は、運転データ測定ステップにおいて、所定の期間内における、空調機器及び空調以外の機器の運転状況を、適宜変化させることが好ましい。
【０００９】
また、このような本発明に係る空調システム内機器の使用エネルギー量推定方法は、空調機器が、１台の室外機に２台以上の室内機が接続された構成のときに、運転データ測定ステップにおいて、室内機毎に、設定温度、室温及び稼働状況を測定し、且つ、モデル作成ステップにおいて、当該室内機毎の設定温度、室温及び稼働状況を入力とするモデルを作成することになる。このことより、推定値算出ステップにおいて、当該空調機器における室内機毎の使用エネルギー量の推定値に対する寄与度を算出するように構成されることが可能となる。
【００１０】
また、このような本発明に係る空調システム内機器の使用エネルギー量推定方法は、空調機器に冷房・暖房などの運転モードが設けられているときに、当該運転モード毎に運転データ測定ステップ、モデル作成ステップ及び推定値算出ステップを実行することが好ましい。
【００１１】
さらに、このような本発明に係る空調システム内機器の使用エネルギー量推定方法は、空調システムで測定された総エネルギー量の実測値をＷとし、モデル作成ステップで作成された空調システムのモデルにより求められる総エネルギー量の推定値をＷＴとし、モデルにより求められる空調機器毎の使用エネルギー量の推定値をＷｐとしたとき、推定値算出ステップは、補正された空調機器毎の使用エネルギー量の推定値Ｗｐ′を、次式
Ｗｐ′  ＝  Ｗｐ・Ｗ／ＷＴ
により算出することが好ましい。
【００１２】
また、本発明に係る使用エネルギー量推定装置は、空調システムを構成する機器の定義情報が記憶されたデータ関係定義情報ファイルと、運転データが記憶された運転データファイルと、空調機器毎の使用電力量の推定値が記憶される推定結果ファイルと、データ関係定義情報ファイルから定義情報を読み出し、運転データファイルから運転データを読み出して、上述の使用エネルギー量推定方法のいずれかにより、空調システムに含まれる空調機器毎の使用エネルギー量の推定値を算出して推定結果ファイルに記憶させるモデル作成・推定演算部と、を有することを特徴とする。
【発明の効果】
【００１３】
本発明に係る空調システム内機器の使用エネルギー量推定方法、及び、この方法が実装された使用エネルギー量推定装置を以上のような構成とすると、簡易な方法で空調システム内の各空調機器の使用エネルギー量を精度良く推定することができる。
【図面の簡単な説明】
【００１４】
【図１】空調システムを構成する機器の一例を示す説明図である。
【図２】使用エネルギー量推定装置の構成を示すブロック図である。
【図３】上記使用エネルギー量推定装置における処理を示すフローチャートである。
【発明を実施するための形態】
【００１５】
以下、本発明の好ましい実施形態について図面を参照して説明する。ここでは、本発明に係る空調システム内機器の使用エネルギー量推定方法を、ビル空調管理システムに適用した場合について説明する。また、以降の説明では、電力により稼働する空調機器及び空調以外の機器から構成される場合を示す（そのため、「エネルギー量」を「電力量」として説明する）。
【００１６】
図１に示すように、ビル空調管理システム内には、複数の空調機器と空調以外の機器（例えば、照明等）とが存在する。また、空調機器は、室外機と室内機とから構成され、１台の室外機には１台以上の室内機が接続されている。このようなビル空調管理システムに対して、本発明に係る使用エネルギー量推定方法は、所定の時間内において、空調機器については、室内機毎に、その室内機の稼働状況（発停情報）、設定温度、室温を計測し、空調機器以外の機器については稼働状況（発停情報）を計測し、さらに、これらの空調機器及び空調以外の機器が消費した総電力量を測定して、モデル化を行うことにより、各空調機器の使用電力量を推定するものである。
【００１７】
それでは、図１に示すように、３台の室外機Ａ〜Ｃと、それに接続された４台の室内機Ａ，Ｂ，Ｃ１，Ｃ２（室外機Ｃに２台の室内機Ｃ１，Ｃ２が接続されている）と、空調以外の機器Ｄと、を有する空調システムにおける使用エネルギー量推定方法について説明する。この空調システムにおいて、一定時間内に測定された総電力量をＷとし、同一時間内の空調室外機Ａに対応する室内機Ａの空調制御設定温度の平均値をＳａ、同一時間内におけるこの室内機Ａの空調制御対象である室温の平均値をＴａ、同一時間内の室内機Ａの稼働状況（発停情報であって、稼働しているときを「１」、停止しているときを「０」とする）をＯａとする。同様に、空調室外機Ｂ，Ｃ及びその室内機Ｂ，Ｃ１，Ｃ２についても各変数を定義する。なお、室外機Ｃには２台の室内機Ｃ１，Ｃ２が接続されているが、それぞれの室内機Ｃ１，Ｃ２にデータを区分して扱うものとする。さらに、この空調システムには空調以外の機器Ｄが存在し、その稼働状況をＯｄとする。
【００１８】
このような空調システムにおいて、まず、ある期間において空調運転を実施し、上記運転データ（総電力量Ｗ、各室内機の設定温度Ｓｘ、室温Ｔｘ，発停情報Ｏｘ（ｘ＝Ａ，Ｂ，Ｃ１，Ｃ２）及び空調以外の機器の発停情報Ｏｄ）を測定し、記録する。そして、このようにして測定された運転データに基づき、総電力量Ｗを、稼働状況、設定温度、室温を入力値として表現するモデルを作成する。ここで用いるモデルは、ニューラルネットワークモデルのような非線形モデルでも良いが、ここでは、次式（１）に示されるような線形モデルを例として説明する。なお、以下の式（１）において、α１，α２，・・・，β４はモデル係数である。
【００１９】
Ｗ＝α１・Ｏａ＋α２・Ｏｂ＋α３・Ｏｃ１＋α４・Ｏｃ２＋α５・Ｏｄ
＋β１・（Ｔａ−Ｓａ）・Ｏａ  ＋  β２・（Ｔｂ−Ｓｂ）・Ｏｂ
＋β３・（Ｔｃ１−Ｓｃ１）・Ｏｃ１
＋β４・（Ｔｃ２−Ｓｃ２）・Ｏｃ２                                  （１）
【００２０】
この式（１）に示すモデルは、総電力量Ｗが、各機器の稼働状況に加えて、空調機器については、設定温度と室温との差に最も影響を受けるという考え方に基づいている。なお、測定される総電力量Ｗに含まれる、この空調システム内に設置されている空調以外の機器Ｄ（例えば、照明等）については、稼働状況のみをモデルに含めて計算している。
【００２１】
そして、上述のようにして測定された運転データを用いてフィッティング（重回帰分析）を行い、式（１）におけるモデル係数（α１，α２，・・・，β４）を決定する。なお、適切なモデルを作成するために、上述のモデル構築用の運転データ計測時に、各機器の運転状況に、適宜ランダムもしくは、個々の動きの一部が他の動きと相関がないような変動を与えることが望ましい。
【００２２】
以上のようにして空調システムのモデルが作成されると（運転データによりフィッティングされると）、モデル係数を用いて各空調機器の推定電力量を次式（２）〜（６）により算出する。ここで、Ｗａ，Ｗｂ，Ｗｃは、それぞれ、室外機Ａ，Ｂ，Ｃの推定電力量であり、Ｗｃ１，Ｗｃ２は、室外機Ｃにおける室内機Ｃ１，Ｃ２の電力量寄与度合いを表す値である。
【００２３】
Ｗａ    ＝  α１・Ｏａ    ＋  β１・（Ｔａ−Ｓａ）・Ｏａ                （２）
Ｗｂ    ＝  α２・Ｏｂ    ＋  β２・（Ｔｂ−Ｓｂ）・Ｏｂ                （３）
Ｗｃ１  ＝  α３・Ｏｃ１  ＋  β３・（Ｔｃ１−Ｓｃ１）・Ｏｃ１          （４）
Ｗｃ２  ＝  α４・Ｏｃ２  ＋  β４・（Ｔｃ２−Ｓｃ２）・Ｏｃ２          （５）
Ｗｃ    ＝  Ｗｃ１  ＋  Ｗｃ２                                          （６）
【００２４】
このように、本発明に係る空調システム内機器の使用エネルギー量推定方法によると、空調システム内の総電力量と、このシステム内にある各空調機器の稼働状況（発停情報）、設定温度、室内温度とを測定することによりモデルを作成することができるため、このモデルを用いることにより、簡易な方法で、この空調システム内の各空調機器の使用電力量を推定することができる。すなわち、上記のようにして作成されたモデルに対して、使用電力量を推定したい機器の設定温度、室温及び稼働状況を入力することで、使用電力量の推定値を求めることができる。なお、モデル作成においては、空調システム内に含まれる空調以外の機器についても、稼働状況を測定して変数として加えることにより、総電力量への影響を考慮することができるため、空調機器の使用電力量の推定値の精度を向上させることができる。
【００２５】
ところで、本発明に係る空調システム内機器の使用エネルギー量推定方法は、入力値として現在より先の時刻の予測値を用いた場合でも、既に測定済みの過去値を用いた場合にでも適用できるが、過去値を用いた場合で、総電力量Ｗも既知の場合は、その既知の総電力量と、次式（７）で求められる総電力量の推定値ＷＴとを用いて、推定値Ｗａ〜Ｗｃを次式（８）〜（１０）に示すように総電力量が一致する値に補正した補正値Ｗａ′〜Ｗｃ′を推定値としても良い。このような補正をすることにより、推定値ＷＴを実測値Ｗに合わせて、各空調機器の使用電力量の推定値の精度を向上させることができる。
【００２６】
ＷＴ＝α１・Ｏａ＋α２・Ｏｂ＋α３・Ｏｃ１＋α４・Ｏｃ２＋α５・Ｏｄ
＋β１・（Ｔａ−Ｓａ）・Ｏａ  ＋  β２・（Ｔｂ−Ｓｂ）・Ｏｂ
＋β３・（Ｔｃ１−Ｓｃ１）・Ｏｃ１
＋β４・（Ｔｃ２−Ｓｃ２）・Ｏｃ２                                （７）
Ｗａ′  ＝  Ｗａ・Ｗ／ＷＴ                                              （８）
Ｗｂ′  ＝  Ｗｂ・Ｗ／ＷＴ                                              （９）
Ｗｃ′  ＝  Ｗｃ・Ｗ／ＷＴ                                              （１０）
【００２７】
ここで、一般に、空調設備では暖房、冷房、送風など、複数種類の運転モードで運転がなされる場合がある。そのような場合に、本発明に係る空調システム内機器の使用エネルギー量推定方法を適用する場合には、運転モードを考慮したモデル、例えば、運転モード毎にモデルを用意するなどをすることにより対応可能である。
【００２８】
また、モデル作成時に、線形モデルを採用した場合は、入力データに多重共線性があるとモデル係数が容易に特定できないことがある。そのような場合には、線形モデルとしてＰＬＳ（部分最小二乗法）を用いてモデル化することで、多重共線性の問題を回避してモデルを作成することができる。
【００２９】
さらに、連続運転対象での推定を継続して行う場合には、モデル精度を高めるために、継続的に測定された運転データに基づき、モデルを逐次的に更新することも有効である。
【００３０】
なお、以上の説明では、使用電力量の推定を行う空調システムが、３台の室外機及び４台の室内機からなる空調機器と、１組の空調以外の機器とから構成されている場合について説明したが、空調機器等の台数はこの例に限定されない。例えば、ｍ台の室内機と、ｎ組の空調以外の機器から構成される空調システムについては、次式（１１）によりモデル化される。ここで、Ｗは総電力量であり、Ｏ１iは第ｉ番目の室内機の発停情報であり、Ｓiは第ｉ番目の室内機の設定温度であり、Ｔiは第ｉ番目の室内機の空調制御対象である室温であり、Ｏ２jは第ｊ番目の空調以外の機器の発停情報であり、α１i，α２j，βiはモデル係数である。
【００３１】
Ｗ＝Σα１i・Ｏ１i＋Σα２j・Ｏ２j＋Σβi・（Ｔi−Ｓi）・Ｏ１i        （１１）
但し、ｉ＝１…ｍ，ｊ＝１…ｎ
【００３２】
また、第ｓ番目から第ｔ番目までの室内機が接続された室外機の推定電力量Ｗｐは、次式（１２）で求められる。
【００３３】
Ｗｐ＝Σ（α１k・Ｏ１k＋βk・（Ｔk−Ｓk）・Ｏ１k）                    （１２）
但し、ｋ＝ｓ…ｔ
【００３４】
さらに、上述の式（１１）で求められる総電力量の実績値Ｗ及び推定値ＷＴにより、各室外機の推定電力量を補正したＷｐ′は、次式（１３）で表される。
【００３５】
Ｗｐ′＝Ｗｐ・Ｗ／ＷＴ                                                （１３）
【００３６】
次に、このような空調システム内機器の使用エネルギー量推定方法が実装された使用エネルギー量推定装置の構成について図２及び図３を用いて説明する。この使用エネルギー量推定装置１０は、上述の使用エネルギー量推定方法が実装されたモデル作成・推定演算部１１と、空調システム内機器の構成等の定義情報が記憶されたデータ関係定義情報ファイル１２と、測定された運転データが記憶された運転データファイル１３と、個別の空調機器の推定電力量が記憶される推定結果ファイル１４と、を有して構成される。なお、モデル作成・推定演算部１１は、例えば、コンピュータが有するＣＰＵで処理されるプログラムとして実装される。また、データ関係定義情報ファイル１２、運転データファイル１３及び推定結果ファイル１４は、この使用エネルギー量推定装置１０が有するハードディスク等の記憶装置に記憶されている。また、データ関係定義情報ファイル１２及び運転データファイル１３には、それぞれ、予めデータが記憶されているものとする。
【００３７】
このような構成の使用エネルギー量推定装置１０において、モデル作成・推定演算部１１は、まず、データ関係定義情報ファイル１２から、空調システムを構成する機器の構成情報を読み出す（ステップＳ１００）。また、モデル作成・推定演算部１１は、運転データファイル１３から、所定の期間の運転データを読み出す（ステップＳ１１０）。そして、これらの構成情報及び運転データから、上述の式（１）に基づいてフィッティングを行い、モデル係数を求めることで、対象とする空調システムの総電力量推定モデルを作成する（ステップＳ１２０）。さらに、このようにして作成された総電力量推定モデルに基づいて、上述の式（２）〜（６）により、この空調システムを構成する空調機器の個別の推定電力量を算出し（ステップＳ１３０）、その結果を推定結果ファイル１４に書き出す（ステップＳ１４０）。
【００３８】
なお、以上の説明では、電力により稼働する空調機器の使用電力量を推定する方法について説明したが、例えば、ガスにより稼働する空調機器に本発明に係る推定方法を用いても良いことは言うまでもない。
【符号の説明】
【００３９】
１０  使用エネルギー量推定装置
１１  モデル作成・推定演算部
１２  データ関係定義情報ファイル
１３  運転データファイル
１４  推定結果ファイル
【図１】
2012017930.tif 000003
【図２】
2012017930.tif 000004
【図３】
2012017930.tif 000005