(57)【特許請求の範囲】
【請求項１】
各需要者側に設置され少なくともガス流量を計測する計測装置と、前記計測装置から少なくともデータを受信するセンターとからなるガス使用状況判断システムであって、
前記計測装置は、
流路内のガス圧力に応じた信号を出力する圧力センサ、及び、流路内のガス流量に応じた信号を出力する流量センサの少なくとも一方からなる計測センサと、
前記計測センサにより出力された信号に基づく計測値データを送信する送信手段と、を備え、
前記センターは、
前記送信手段により送信された計測値データを受信する受信手段と、
前記受信手段により受信された計測値データの所定以上の変化時からの微小時間中に得られる振動波形からガス漏れ及び使用ガス器具の少なくとも一方を判断する判断手段と、を備える
ことを特徴とするガス使用状況判断システム。
【請求項２】
前記計測装置は、
前記計測センサにより出力された信号の所定以上の変化時にトリガ信号を発生させるトリガ信号発生手段と、
前記トリガ信号発生手段によりトリガ信号が発生された場合に、前記計測センサのサンプリング時間を通常のサンプリング時間よりも短縮するサンプリング時間調整手段と、
をさらに備えることを特徴とする請求項１に記載のガス使用状況判断システム。
【請求項３】
前記送信手段は、前記微小時間経過後に、前記計測センサにより出力された信号に基づく計測値データを、前記微小時間分一括して送信する
ことを特徴とする請求項２に記載のガス使用状況判断システム。
【請求項４】
前記サンプリング時間調整手段は、前記微小時間だけサンプリング時間を短縮した後、サンプリング時間の短縮状態を解除する
ことを特徴とする請求項２又は請求項３のいずれかに記載のガス使用状況判断システム。
【請求項５】
前記センターは、
前記受信手段により受信された前記計測値データに基づく振動波形と、所定の振動波形との類似度推移を算出する類似度推移算出手段をさらに備え、
前記判断手段は、前記類似度推移算出手段により算出された類似度推移に基づいて、ガス漏れ及び使用ガス器具の少なくとも一方を判断する
ことを特徴とする請求項１から請求項４のいずれか１項に記載のガス使用状況判断システム。
【請求項６】
前記センターは、
前記受信手段により受信された前記計測値データに基づく振動波形を解析して、周波数と振幅との相関を示すスペクトルデータを算出する解析手段をさらに備え、
前記判断手段は、前記解析手段により算出されたスペクトルデータに基づいて、ガス漏れ及び使用ガス器具の少なくとも一方を判断する
ことを特徴とする請求項１から請求項４のいずれか１項に記載のガス使用状況判断システム。
【請求項７】
前記判断手段は、ガス漏れ及び使用が開始されたガス器具の双方を判断可能であって、使用が開始されたガス器具を判断するのに先立って、ガス漏れ判断の処理を実行する
ことを特徴とする請求項１から請求項６のいずれか１項に記載のガス使用状況判断システム。
【請求項８】
前記センターは、
リコール対象となるガス器具の情報を記憶したリコール器具記憶手段と、
リコール対象となるガス器具に対する保安処理を実行する保安手段と、をさらに備え、
前記判断手段は、少なくとも使用ガス器具を判断すると共に、前記リコール器具記憶手段により記憶されたリコール対象となるガス器具の情報に基づいて、使用ガス器具がリコール対象となるガス器具であるかを判断し、
前記保安手段は、前記判断手段により使用ガス器具がリコール対象のガス器具であると判断された場合、当該リコール対象のガス器具に対する保安処理を実行させる
ことを特徴とする請求項１から請求項７のいずれか１項に記載のガス器具判断システム。
【請求項９】
各需要者側に設置され少なくともガス流量を計測する計測装置と、前記計測装置から少なくともデータを受信するセンターとからなるガス使用状況判断システムのガス使用状況判断方法であって、
前記計測装置は、
流路内のガス圧力に応じた信号を出力する圧力センサ、及び、流路内のガス流量に応じた信号を出力する流量センサの少なくとも一方からなる計測センサにより出力された信号に基づく計測値データを送信する送信工程を有し、
前記センターは、
前記送信工程において送信された計測値データを受信する受信工程と、
前記受信工程において受信された計測値データの所定以上の変化時からの微小時間中に得られる振動波形からガス漏れ及び使用ガス器具の少なくとも一方を判断する判断工程と、を有する
ことを特徴とするガス使用状況判断システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ガス使用状況判断システム及びガス使用状況判断方法に関する。
【背景技術】
【０００２】
従来、使用ガス器具やガス漏れを判断するガス使用状況判断装置が提案されている。このガス使用状況判断装置では、流路内を流れるガスの流量値パターンから、使用ガス器具を判断する（例えば特許文献１〜３参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００９−１０９２５６号公報
【特許文献２】特開２００９−１４６８６号公報
【特許文献３】特開２００７−１９９００２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかし、従来装置では、使用ガス器具を判断するにあたり、ガス使用状況判断装置側で判断しなければならず、判断にあたってガス使用状況判断装置の電力消費量が多くなってしまう。特に、ガス使用状況判断装置がガスメータである場合には、消費電力の増大は１０年の駆動期間に影響を与えるため、顕著な問題となってしまう。
【０００５】
本発明はこのような従来の課題を解決するためになされたものであり、その発明の目的とするところは、消費電力を抑えつつ、ガス漏れ及び使用ガス器具の少なくとも一方を判断することが可能なガス使用状況判断システム及びガス使用状況判断方法を提供することにある。
【課題を解決するための手段】
【０００６】
本発明のガス使用状況判断システムは、各需要者側に設置され少なくともガス流量を計測する計測装置と、前記計測装置から少なくともデータを受信するセンターとからなるガス使用状況判断システムであって、前記計測装置は、流路内のガス圧力に応じた信号を出力する圧力センサ、及び、流路内のガス流量に応じた信号を出力する流量センサの少なくとも一方からなる計測センサと、前記計測センサにより出力された信号に基づく計測値データを送信する送信手段と、を備え、前記センターは、前記送信手段により送信された計測値データを受信する受信手段と、前記受信手段により受信された計測値データの所定以上の変化時からの微小時間中に得られる振動波形からガス漏れ及び使用ガス器具の少なくとも一方を判断する判断手段と、を備えることを特徴とする。
【０００７】
このガス使用状況判断システムによれば、計測装置は計測値データを送信し、センターは、受信した計測値データの所定以上の変化時からの微小時間中に得られる振動波形からガス漏れ及び使用ガス器具の少なくとも一方を判断する。ここで、本件出願人らは、ガス漏れ発生時や、ガス器具の使用開始時及び終了時など、計測値データの所定以上の変化時から微小時間経過するまでの波形に、ガス漏れやガス器具毎に固有の振動を示すことを見出した。このため、所定以上の変化時からの微小時間中に得られる振動波形からガス漏れ及び使用ガス器具の少なくとも一方を判断することができる。特に、計測装置は計測センサにより出力された信号に基づく計測値データを送信するだけでよく、判断についての演算を行う必要がない。このため、基本的に電池駆動となる計測装置の消費電力を抑えることができる。従って、消費電力を抑えつつ、ガス漏れ及び使用ガス器具の少なくとも一方を判断することができる。
【０００８】
また、本発明のガス使用状況判断システムにおいて、前記計測装置は、前記計測センサにより出力された信号の所定以上の変化時にトリガ信号を発生させるトリガ信号発生手段と、前記トリガ信号発生手段によりトリガ信号が発生された場合に、前記計測センサのサンプリング時間を通常のサンプリング時間よりも短縮するサンプリング時間調整手段と、をさらに備えることが好ましい。
【０００９】
このガス使用状況判断システムによれば、計測センサにより出力された信号の所定以上の変化時にトリガ信号を発生し、トリガ信号が発生された場合に、計測センサのサンプリング時間を通常のサンプリング時間よりも短縮する。このため、微小時間における波形の振動を適切にとらえることができ、適切に計測することができる。
【００１０】
また、本発明のガス使用状況判断システムにおいて、前記送信手段は、前記微小時間経過後に、前記計測センサにより出力された信号に基づく計測値データを、前記微小時間分一括して送信することが好ましい。
【００１１】
このガス使用状況判断システムによれば、微小時間経過後に、計測センサにより出力された信号に基づく計測値データを、微小時間分一括して送信するため、計測の都度に送信する場合と比較して計測装置による送信回数の増加を抑えることができ、一層消費電力を抑えることができる。
【００１２】
また、本発明のガス使用状況判断システムにおいて、前記サンプリング時間調整手段は、前記微小時間だけサンプリング時間を短縮した後、サンプリング時間の短縮状態を解除することが好ましい。
【００１３】
このガス使用状況判断システムによれば、微小時間だけサンプリング時間を短縮した後、サンプリング時間の短縮状態を解除するため、ガス漏れや使用ガス器具の判断に必要なる振動波形が得られた後は、詳細なサンプリングを必要とせず、そのような場合にサンプリング時間の短縮状態を解除することで、一層消費電力を抑えることができる。
【００１６】
また、本発明のガス使用状況判断システムにおいて、前記センターは、前記受信手段により受信された前記計測値データに基づく振動波形と所定の振動波形との類似度推移を算出する類似度推移算出手段をさらに備え、前記判断手段は、前記類似度推移算出手段により算出された類似度推移に基づいて、ガス漏れ及び使用ガス器具の少なくとも一方を判断することが好ましい。
【００１７】
このガス使用状況判断システムによれば、計測値データに基づく振動波形と所定の振動波形との類似度推移を算出し、算出された類似度推移に基づいて、ガス漏れ及び使用ガス器具の少なくとも一方を判断する。ここで、所定の振動波形が特定のガス器具の振動波形であるとする。そして、計測センサから得られた信号に基づく振動波形が特定のガス器具の振動波形である場合、算出される類似度は高くなる傾向にあり、類似度推移についても全体的に高い値を示す傾向にある。また、計測センサから得られた信号に基づく振動波形が他のガス器具の振動波形である場合には、類似度は高くならず、類似度推移についても高くならないが、類似度推移はそのガス器具特有の推移を示す。ガス漏れ時も同様にガス漏れ特有の推移を示す。従って、上記のように類似度推移に基づいて判断を実行することで、振動波形の特徴をとらえて、ガス漏れ及び使用ガス器具の少なくとも一方を判断することができる。なお、センターは、所定の振動波形を生成してもよいし、予めデータテーブルとして所定の振動波形を記憶しておいてもよいし、これらの組み合わせであってもよい。
【００１８】
また、本発明のガス使用状況判断システムにおいて、前記センターは、前記受信手段により受信された前記計測値データに基づく振動波形を解析して、周波数と振幅との相関を示すスペクトルデータを算出する解析手段をさらに備え、前記判断手段は、前記解析手段により算出されたスペクトルデータに基づいて、ガス漏れ及び使用ガス器具の少なくとも一方を判断することが好ましい。
【００１９】
このガス使用状況判断システムによれば、計測値データに基づく振動波形を解析して、周波数と振幅との相関を示すスペクトルデータを算出し、算出されたスペクトルデータに基づいて、ガス漏れ及び使用ガス器具の少なくとも一方を判断する。ここで、ガス器具使用時及び終了時やガス漏れ時における振動波形それぞれには、圧力や流量の波形の周波数及び振幅に特徴があらわれる。よって、波形を解析して周波数と振幅との相関を示すスペクトルデータを得ると共に、このスペクトルデータに基づいて判断を実行することで、ガス漏れ及び使用ガス器具の少なくとも一方を判断することができる。
【００２０】
また、本発明のガス使用状況判断システムにおいて、前記判断手段は、ガス漏れ及び使用が開始されたガス器具の双方を判断可能であって、使用が開始されたガス器具を判断するのに先立って、ガス漏れ判断の処理を実行することが好ましい。
【００２１】
このガス使用状況判断システムによれば、ガス漏れ及び使用が開始されたガス器具の双方を判断可能であって、使用が開始されたガス器具を判断するのに先立って、ガス漏れ判断の処理を実行する。このため、緊急性の高いガス漏れ判断を優先して行うこととなり、安全性の向上を図ることができる。
【００２２】
また、本発明のガス使用状況判断システムにおいて、前記センターは、リコール対象となるガス器具の情報を記憶したリコール器具記憶手段と、リコール対象となるガス器具に対する保安処理を実行する保安手段と、をさらに備え、前記判断手段は、少なくとも使用ガス器具を判断すると共に、前記リコール器具記憶手段により記憶されたリコール対象となるガス器具の情報に基づいて、使用ガス器具がリコール対象となるガス器具であるかを判断し、前記保安手段は、前記判断手段により使用ガス器具がリコール対象のガス器具であると判断された場合、当該リコール対象のガス器具に対する保安処理を実行させることが好ましい。
【００２３】
このガス使用状況判断システムによれば、センターはリコール対象となるガス器具の情報を記憶し、使用ガス器具がリコール対象のガス器具である場合、当該リコール対象のガス器具に対する保安処理を実行させる。このため、リコール品の回収や使用禁止を行うことができ、安全性を向上させることができる。加えて、リコール品をセンター側で登録しておけばよく、各需要者側において登録等する必要がないことから、リコール品に対する保安処理についての工数を低減することができる。
【００２４】
また、本発明のガス使用状況判断方法は、各需要者側に設置され少なくともガス流量を計測する計測装置と、前記計測装置から少なくともデータを受信するセンターとからなるガス使用状況判断システムのガス使用状況判断方法であって、前記計測装置は、流路内のガス圧力に応じた信号を出力する圧力センサ、及び、流路内のガス流量に応じた信号を出力する流量センサの少なくとも一方からなる計測センサにより出力された信号に基づく計測値データを送信する送信工程を有し、前記センターは、前記送信工程において送信された計測値データを受信する受信工程と、前記受信工程において受信された計測値データの所定以上の変化時からの微小時間中に得られる振動波形からガス漏れ及び使用ガス器具の少なくとも一方を判断する判断工程と、を有することを特徴とする。
【００２５】
このガス使用状況判断方法によれば、計測装置は計測値データを送信し、センターは、受信した計測値データの所定以上の変化時からの微小時間中に得られる振動波形からガス漏れ及び使用ガス器具の少なくとも一方を判断する。ここで、本件出願人らは、ガス漏れ発生時や、ガス器具の使用開始時及び終了時など、計測値データの所定以上の変化時から微小時間経過するまでの波形に、ガス漏れやガス器具毎に固有の振動を示すことを見出した。このため、所定以上の変化時からの微小時間中に得られる振動波形からガス漏れ及び使用ガス器具の少なくとも一方を判断することができる。特に、計測装置は計測センサにより出力された信号に基づく計測値データを送信するだけでよく、判断についての演算を行う必要がない。このため、基本的に電池駆動となる計測装置の消費電力を抑えることができる。従って、消費電力を抑えつつ、ガス漏れ及び使用ガス器具の少なくとも一方を判断することができる。
【発明の効果】
【００２６】
本発明によれば、消費電力を抑えつつ、ガス漏れ及び使用ガス器具の少なくとも一方を判断することが可能なガス使用状況判断システム及びガス使用状況判断方法を提供することができる。
【図面の簡単な説明】
【００２７】
【図１】本発明の実施形態に係るガス使用状況判断システムの構成図である。
【図２】図１に示したガスメータの詳細を示す構成図である。
【図３】ガス器具の使用過程における流量推移の一例を示す図である。
【図４】図１に示した管理センターの詳細を示す構成図である。
【図５】図４に示した生成部により生成されるガス漏れ振動波形の概略を示す図である。
【図６】ガス漏れ時における圧力変化を示す図である。
【図７】ガス漏れ時における連続ＮＣＣを示すグラフである。
【図８】ガス器具使用開始時における圧力変化を示すグラフであって、（ａ）はガステーブル使用開始時における圧力変化を示し、（ｂ）は小型湯沸器使用開始時における圧力変化を示し、（ｃ）は給湯器使用開始時における圧力変化を示している。
【図９】図４に示した類似度推移算出部により算出される連続ＮＣＣを示すグラフであって、（ａ）はガステーブル使用開始時における連続ＮＣＣを示し、（ｂ）は小型湯沸器使用開始時における連続ＮＣＣを示し、（ｃ）は給湯器使用開始時における連続ＮＣＣを示している。
【図１０】ガス器具使用終了時における圧力変化を示すグラフであって、（ａ）はガステーブル使用終了時における圧力変化を示し、（ｂ）は小型湯沸器使用終了時における圧力変化を示し、（ｃ）は給湯器使用終了時における圧力変化を示している。
【図１１】図４に示した類似度推移算出部により算出される連続ＮＣＣを示すグラフであって、（ａ）はガステーブル使用終了時における連続ＮＣＣを示し、（ｂ）は小型湯沸器使用終了時における連続ＮＣＣを示し、（ｃ）は給湯器使用終了時における連続ＮＣＣを示している。
【図１２】本実施形態に係るガスメータの動作の一例を示すフローチャートである。
【図１３】本実施形態に係る管理センターの動作の一例を示すメインフローチャートである。
【図１４】図１３に示したガス漏れ／開始ガス器具判断処理（Ｓ１３）の詳細を示すフローチャートである。
【図１５】図１３に示した終了ガス器具判断処理（Ｓ１７）の詳細を示すフローチャートである。
【図１６】第２実施形態に係る管理センターの詳細を示す構成図である。
【図１７】図１６に示した解析部により算出されるスペクトルデータを示すグラフであって、ガス漏れが発生したときの圧力波形をフーリエ変換して得られるスペクトルデータを示すグラフである。
【図１８】図１６に示した解析部により算出されるスペクトルデータを示すグラフであって、（ａ）はガステーブル使用開始時におけるスペクトルデータを示し、（ｂ）は小型湯沸器使用開始時におけるスペクトルデータを示し、（ｃ）は給湯器使用開始時におけるスペクトルデータを示している。
【図１９】図１６に示した解析部により算出されるスペクトルデータを示すグラフであって、（ａ）はガステーブル使用終了時におけるスペクトルデータを示し、（ｂ）は小型湯沸器使用終了時におけるスペクトルデータを示し、（ｃ）は給湯器使用終了時におけるスペクトルデータを示している。
【図２０】図１３に示したガス漏れ／開始ガス器具判断処理（Ｓ１３）の詳細を示す第２のフローチャートである。
【図２１】図１３に示した終了ガス器具判断処理（Ｓ１７）の詳細を示す第２のフローチャートである。
【発明を実施するための形態】
【００２８】
以下、本発明の好適な実施形態を図面に基づいて説明する。図１は、本発明の実施形態に係るガス使用状況判断システムの構成図である。ガス使用状況判断システム１は、ガスストーブ、ファンヒータ、小型湯沸器、給湯器、床暖房及びガステーブルなどの各ガス器具１０に燃料ガスを供給するものであって、複数のガス器具１０と、ガス供給元の調整器２０と、配管３１，３２と、ガスメータ（計測装置）４０と、管理センター（センター）５０とを備えている。なお、図１に示す例では、ガスメータ４０を計測装置の一例として挙げるが、計測装置はガスメータ４０に限るものではない。
【００２９】
調整器２０は上流からの燃料ガスを所定圧力に調整して第１配管３１に流すものである。第１配管３１は、調整器２０とガスメータ４０とを接続するものである。第２配管３２はガスメータ４０とガス器具１０とを接続する配管である。ガスメータ４０は、各家庭に設置され、少なくともガス流量を計測して積算流量を表示するものである。このようなガス供給システム１では、ガスメータ４０内に第１配管３１及び第２配管３２とつながる流路が形成されており、調整器２０を通じて流れてきた燃料ガスは第１配管３１からガスメータ４０、及び第２配管３２を通じてガス器具１０に到達し、ガス器具１０において燃焼されることとなる。
【００３０】
管理センター５０は、ガスメータ４０からデータを受信するものであって、受信したデータに基づいて、ガス漏れや各家庭において使用ガス器具１０を判断するものである。なお、本実施形態においてガスメータ４０は家庭に設置されるものとして説明しているが、ガスメータ４０は一般の家庭に設置されるものに限らず、料理を行う店舗など、ガス使用を必要とする需要者であれば、他の場所に設置されていてもよい。
【００３１】
このようなシステム１では、ガスメータ４０で計測された圧力や流量のデータを管理センター５０側に送信し、管理センター５０側においてガス漏れ及び使用ガス器具１０の少なくとも一方を判断するようになっている。このとき、管理センター５０では、圧力や流量の所定以上の変化時から微小時間（最大で２秒）中に得られた圧力や流量のデータに基づいて、ガス漏れ及び使用ガス器具１０の少なくとも一方を判断する。使用ガス器具１０とは、使用が開始したガス器具１０、及び使用が終了したガス器具１０の少なくとも一方を含む概念である。
【００３２】
ここで、本件出願人らは、ガス器具１０の使用開始時や終了時、及びガス漏れ発生時など、圧力や流量の所定以上の変化時から微小時間経過するまでの波形に、ガス漏れやガス器具１０毎に固有の振動を示すことを見出した。このため、管理センター５０は、所定以上の変化時からの微小時間中に得られる振動波形に基づいて、ガス漏れ及び使用ガス器具の少なくとも一方を判断する。特に、本実施形態では管理センター５０側で判断を行うこととなり、ガスメータ４０での演算量を減らして、ガスメータ４０における消費電力を抑えている。
【００３３】
図２は、図１に示したガスメータ４０の詳細を示す構成図である。図２に示すようにガスメータ４０は、圧力センサ（計測センサ）４１と、流量センサ（計測センサ）４２と、制御部４３と、送受信部（送信手段）４４と、トリガ信号発生部（トリガ信号発生手段）４５とを有している。
【００３４】
圧力センサ４１は、ガスメータ４０の流路内におけるガス圧力に応じた計測値の信号を出力するものであって、ピエゾ抵抗式や静電容量式などのセンサによって構成される。流量センサ４２は、ガスメータ４０の流路内におけるガス流量に応じた計測値の信号を出力するものであって、超音波センサやフローセンサなどで構成される。
【００３５】
制御部４３は、ガスメータ４０の全体を制御するものであって、センサ４１，４２により出力された信号に基づいて圧力値や流量値を算出したり、流量の積算表示の制御を実行したりするものである。この制御部４３はマイコンによって構成される。
【００３６】
送受信部４４は、計測値データを送信するものである。ここで、計測値データとは、制御部４３によって算出された圧力値や流量値の情報であってもよいし、センサ４１，４２によって出力されたそのままの信号の情報であってもよい。また、送受信部４４は、管理センター５０から遮断弁を遮断したり復帰させたりするデータを受信する。このような送受信部４４は例えば通信インターフェイスにより構成される。
【００３７】
トリガ信号発生部４５は、センサ４１，４２により出力された信号の所定以上の変化時にトリガ信号を発生させるものである。このトリガ信号は、制御部４３及び各センサ４１，４２に送信される。このようなトリガ信号発生部４５は、例えば微分回路を含んで構成されており、微分回路により所定以上の変化を検出する。
【００３８】
具体的にトリガ信号発生部４５は、ガス器具１０が使用を開始され又はガス漏れが発生したときに、流量が流れ且つ圧力が低下するときの変化を所定以上の変化としてとらえ、第１トリガ信号を出力する。さらに、トリガ信号発生部４５は、ガス器具１０の使用が終了したときに、流量が低下し且つ圧力が上昇するときの変化を所定以上の変化としてとらえ、第２トリガ信号を出力する。なお、トリガ信号発生部４５は、上記２種類のトリガ信号を区別可能に発生させる。
【００３９】
図３は、ガス器具の使用過程における流量推移の一例を示す図である。ガス器具１０の使用過程においては例えば図３に示すような流量が流れる。まず、時刻ｔ１において第１のガス器具１０（例えばガステーブル）が使用開始されたとする。このとき、流量値はＦ１を示す。そして、時刻ｔ２において第２のガス器具１０（例えば給湯器）が使用開始されたとすると、流量値はＦ３（＝Ｆ１＋Ｆ２）を示す。
【００４０】
その後、時刻ｔ３において第１のガス器具１０の使用が終了したとすると、流量値は第２のガス器具１０のみの流量であるＦ２を示す。次いで、時刻ｔ４において第２のガス器具１０についても使用が終了したとすると、流量値は「０」を示す。
【００４１】
トリガ信号発生部４５は、図３で示す時刻ｔ１及び時刻ｔ２のタイミングで第１トリガ信号を発生させる。また、時刻ｔ３及び時刻ｔ４のタイミングで第２トリガ信号を発生させる。
【００４２】
再度、図２を参照する。制御部４３は、サンプリング時間調整部（サンプリング時間調整手段）４３ａを備えている。サンプリング時調整部４３ａは、トリガ信号発生部４５によりトリガ信号が発生された場合に、センサ４１，４２のサンプリング時間を通常のサンプリング時間（例えば流量では２秒、圧力では１０秒）よりも短縮するものである。この際、サンプリング時間調整部４３ａはサンプリング時間を１マイクロ秒に短縮する。これにより、ガスメータ４０は、微小時間における波形の振動を適切にとらえるようにしている。
【００４３】
また、各センサ４１，４２は、トリガ信号が入力されると、高速サンプリングにあわせて信号を出力することとなる。加えて、トリガ信号発生前においてセンサ４１，４２の駆動電流を通常電流よりも小さくしておき、トリガ信号発生後に駆動電流を通常電流に戻すようにしてもよい。
【００４４】
また、サンプリング時間調整部４３ａは、微小時間だけサンプリング時間を短縮した後に、サンプリング時間の短縮状態を解除する。ガス漏れや使用ガス器具１０の判断には、微小時間中の振動波形が必要であり、その後の波形は必要がない。このため、微小時間経過後には通常のサンプリング時間の短縮状態を解除することで、消費電力を抑えるようにしている。
【００４５】
さらに、サンプリング時間調整部４３ａは、トリガ信号発生部４５によりトリガ信号が発生されるまでは、センサ４１，４２のサンプリング時間を通常のサンプリング時間よりも長くしておく。例えば、サンプリング時間調整部４３ａは、流量について１０秒程度に長くし、圧力について３０秒程度に長くする。ここで、トリガ信号が発生していない場合とは、流量や圧力に所定以上の変化がなく、流量や圧力が安定している場合といえる。すなわち、計測の必要性が少ない場合といえる。よって、このような場合に、サンプリング時間調整部４３ａは、サンプリング時間を長くしておくことで、消費電力を軽減させるようにしている。
【００４６】
また、送受信部４４は、微小時間の経過後に、微小時間分一括して計測値データを送信する。これにより、送信回数を抑えて一層消費電力を抑えるようにしている。なお、送信部４４は、微小時間の経過後に、微小時間分一括して計測値データを送信する場合に限らず、微小時間中に１回、微小時間経過後に残りデータを１回など、複数回の送信を行ってもよい。加えて、送受信部４４は、常時送信を行うようになっていてもよい。
【００４７】
さらに、送受信部４４は、第１トリガ信号と第２トリガ信号とのいずれのトリガ信号が発生したかを示す情報についても送信する。これにより、管理センター５０側にいずれのタイミングのデータかを示し、管理センター５０側でのガス漏れ判断、及び、使用ガス器具１０の判断について処理負担の軽減を図るようにしている。
【００４８】
図４は、図１に示した管理センター５０の詳細を示す構成図である。図４に示すように管理センター５０は、送受信部（受信手段）５１と、制御部５２と、記憶部５３とを有している。
【００４９】
送受信部５１は、ガスメータ４０の送受信部４４から送信された計測値データを受信するものである。記憶部５３は、各種判断処理に必要となる情報を記憶したものであって、類似度推移パターン記憶部５３ａと、リコール器具記憶部５３ｂとを備えている。類似度推移パターン記憶部５３ａはガス漏れ及び使用ガス器具１０の判断に必要となる類似度推移のパターンを記憶しており、リコール器具記憶部５３ｂはリコール対象となるガス器具１０の情報（後述のメーカ名及び型番など）を記憶している。なお、類似度推移とは、類似度を連続的に求めて得られるものである。
【００５０】
制御部５２は、送受信部５１により受信された計測値データに基づいて各種処理を実行するものであって、判断部（判断手段）５２ａと、生成部５２ｂと、類似度推移算出部（類似度推移算出手段）５２ｃと、保安部（保安手段）５２ｄとを備えている。
【００５１】
判断部５２ａは、送受信部５１により受信された計測値データの所定以上の変化時から微小時間中に得られる振動波形から、ガス漏れ及び使用ガス器具１０の少なくとも一方を判断するものである。この判断部５２ａは、生成部５２ｂ及び類似度推移算出部５２ｃの演算結果に基づいて、判断を実行する。なお、以下の本実施形態において判断部５２ａは、ガス漏れ及び使用ガス器具１０の双方を判断するものとするが、いずれか一方のみを判断するものであってもよい。
【００５２】
生成部５２ｂは、所定の振動波形を生成するものである。この生成部５２ｂによって生成される所定の振動波形は、後の処理において、ガスメータ４０からの計測値データからなる振動波形と比較される。
【００５３】
類似度推移算出部５２ｃは、ガスメータ４０からの計測値データからなる振動波形と、生成部５２ｂによって生成された所定の振動波形との類似度推移を算出するものである。そして、判断部５２ａは、類似度推移算出部５２ｃによって算出された類似度推移に基づいて、ガス漏れ及び使用ガス器具１０を判断することとなる。
【００５４】
判断の概略を説明すると、例えば生成された所定の振動波形がガス漏れ発生時の振動波形であって、実際にガス漏れが発生したとする。この場合、実際にセンサ４１，４２によって計測される振動波形と生成された振動波形とは、類似度が高くなり、類似度推移についても高くなる。このため、判断部５２ａは、ガス漏れが発生したと判断できる。
【００５５】
また、生成された所定の振動波形がガス漏れ発生時の振動波形であって、実際にガステーブルの使用が開始されたとする。この場合、実際にセンサ４１，４２によって計測される振動波形と生成された振動波形とは、類似度が高くならず、類似度推移についても高くならない。このため、判断部５２ａは、ガス漏れの発生でないと判断する。この際、類似度推移には、ガス漏れ発生時とガステーブルの使用開始時との相違が生じる。この相違は、ガス器具１０毎に異なる。例えば、ガス漏れ発生時とガステーブルの使用開始時との相違と、ガス漏れ発生時と給湯器の使用開始時との相違とは異なっており、判断部５２ａは、相違の状態からガステーブルの使用が開始したと判断できる。
【００５６】
なお、上記は、生成部５２ｂが所定の振動波形としてガス漏れ発生時の振動波形を生成した場合を例に説明したが、これに限らず、生成部５２ｂは、特定のガス器具１０の使用が開始されたときの振動波形を生成するようにしてもよいし、他のガス器具１０の使用が終了したときの振動波形を生成するようにしてもよい。すなわち、生成部５２ｂは、基準となる振動波形を生成すればよく、その振動波形は何の振動波形であってもよい。
【００５７】
加えて、生成部５２ｂは、全ガス器具１０の振動波形（複数の振動波形）を生成し、類似度推移算出部５２ｃは、生成された全ガス器具１０の振動波形と、ガスメータ４０からの計測値データに基づく振動波形との類似度推移を算出して、類似度推移が最も高いガス器具１０を使用ガス器具１０と判断してもよい。また、どの類似度推移についても高くない場合には、ガス漏れの発生と判断してもよい。
【００５８】
次に、ガス漏れ及び使用ガス器具１０の判断についてより詳細に説明する。まず、判断部５２ａは、時刻ｔ１及び時刻ｔ２に示したようなガス流量の増加時には、ガス漏れが発生しているか否かを判断し、ガス漏れが発生していないと判断できる場合に、使用が開始されたガス器具１０を判断する。また、判断部５２ａは、時刻ｔ３及び時刻ｔ４に示したようなガス流量の減少時には、ガス漏れの判断を行わず、使用が終了したガス器具１０を判断する。
【００５９】
なお、以下の実施形態において生成部５２ｂは、基準となる所定の振動波形として、ガス漏れ時における振動波形を生成する場合を例に説明する。さらに、以下の実施形態では、ガス圧力の振動波形に基づいてガス漏れ及び使用ガス器具１０を判断する場合について説明するが、圧力と流量とには一定の相関があるため、流量に基づいてガス漏れ及び使用ガス器具１０を判断してもよい。
【００６０】
図５は、図４に示した生成部５２ｂにより生成されるガス漏れ振動波形の概略を示す図である。図５に示すように、生成部５２ｂは、圧力が時間の経過と共に低下しながら振動するガス漏れ振動波形を生成する。このガス漏れ振動波形は、減衰振動の周波数、ゲイン、及び減衰比を含む２次遅れのステップ応答の式に基づいて生成された波形である。上記したように、本件発明者らは、ガス漏れ発生直後の微小時間において圧力や流量の計測値に振動が発生することについて見出した。このため、生成部５２ｂは、２次遅れのステップ応答の式に基づいてガス漏れ振動波形を生成する。
【００６１】
再度、図４を参照する。類似度推移算出部５２ｃは、受信した計測値データに基づく微小時間中の振動波形と、生成部５２ｂに生成されたガス漏れ振動波形との類似度推移を算出するものである。なお、類似度推移とは、本実施形態において連続的な正規相互相関（NCC：Normalized Cross Correlation）をいう。より具体的には、以下の式（１）により類似度ＲＮＣＣが求められる。類似度推移算出部５２ｃは、この式（１）による類似度ＲＮＣＣの算出を連続的に行うことにより、類似度推移（以下、連続ＮＣＣという）を求める。
【数１】
0005436175.tif 000002
【００６２】
次に、図６を参照してガス漏れ時における圧力変化を説明する。図６は、ガス漏れ時における圧力変化を示す図である。図６に示すように、ガス漏れ発生時には、圧力が低下しつつ振動する波形を示すこととなる。この波形は、図５に示したように生成部５２ｂにより生成されたガス漏れ振動波形と相関が高い。このため、類似度推移の代表値は高い値を示すこととなり、判断部５２ａはガス漏れが発生したと判断することとなる。
【００６３】
図７は、ガス漏れ時における連続ＮＣＣを示すグラフである。なお、図７において実線と破線は、各家庭における配管状態の相違、ガス漏れ箇所の相違、及び、ガス漏れ流量の相違などの条件が異なる場合の連続ＮＣＣを示している。
【００６４】
図７に示すように、ガス漏れ時において圧力変化の発生直後（時刻０秒付近）における連続ＮＣＣは、「０．７」から「０．８」程度の値を示す。しかし、時刻０．０２５秒以降について連続ＮＣＣは「０．９」以上の値を示す。
【００６５】
このように、ガス漏れ発生時において連続ＮＣＣは「０．９」以上の値を示すことから、判断部５２ａは、類似度推移算出部５２ｃにより算出された類似度推移の代表値が閾値（例えば「０．９」）以上である場合に、ガス漏れが発生していると判断する。ここで、代表値とは、類似度全体又は類似度全体のうち特定期間の平均値であってもよいし、圧力や流量の変化が発生してから、ある特定の時刻における類似度であってもよいし、他の値であってもよい。
【００６６】
なお、類似度推移パターン記憶部５３ａに、図７に示したような連続ＮＣＣのパターンを記憶しておき、判断部５２ａは、記憶された連続ＮＣＣのパターンと、算出された連続ＮＣＣとを比較し、両者が近い場合にガス漏れと判断してもよい。
【００６７】
一方、図８に示すようにガス器具使用時には、ガス漏れ時と異なる振動波形を示すこととなる。図８は、ガス器具使用開始時における圧力変化を示すグラフであって、（ａ）はガステーブル使用開始時における圧力変化を示し、（ｂ）は小型湯沸器使用開始時における圧力変化を示し、（ｃ）は給湯器使用開始時における圧力変化を示している。
【００６８】
図８（ａ）に示すように、ガステーブルの使用開始時には圧力が２．９ｋＰａ程度で滑らかに振動する圧力波形が得られる。また、図８（ｂ）に示すように、小型湯沸器の使用終了時には圧力が２．９３ｋＰａを基準にして０．１ｋＰａ強振動する圧力波形が得られる。さらに、図８（ｃ）に示すように、給湯器の使用終了時には圧力が２．９３ｋＰａを基準にして小型湯沸器よりもやや粗い振動を示す圧力波形が得られる。
【００６９】
図９は、図４に示した類似度推移算出部５２ｃにより算出される連続ＮＣＣを示すグラフであって、（ａ）はガステーブル使用開始時における連続ＮＣＣを示し、（ｂ）は小型湯沸器使用開始時における連続ＮＣＣを示し、（ｃ）は給湯器使用開始時における連続ＮＣＣを示している。
【００７０】
ガステーブルの使用が開始した場合、図８（ａ）の振動波形が得られ、類似度推移算出部５２ｃにより算出される連続ＮＣＣは図９（ａ）に示すようになる。すなわち、連続ＮＣＣは、初期的に「１．０」程度を示し、その後「０．２」を下回り、約０．０４秒において「０．９５」まで復帰する。そして、連続ＮＣＣは、約０．１秒において「０．５」程度となり、その後「０．６５」付近までゆっくりと上昇する。
【００７１】
また、小型湯沸器の使用が開始した場合、図８（ｂ）の振動波形が得られ、類似度推移算出部５２ｃにより算出される連続ＮＣＣは図９（ｂ）に示すようになる。すなわち、連続ＮＣＣは、初期的に「１．０」程度を示し、その後「０．２」を下回り、約０．０４秒において「０．９」まで復帰する。そして、連続ＮＣＣは、再度「０．４」程度まで低下し、その後、「０．７」付近までゆっくりと上昇する。
【００７２】
また、給湯器の使用が開始した場合、図８（ｃ）の振動波形が得られ、類似度推移算出部５２ｃにより算出される連続ＮＣＣは図９（ｃ）に示すようになる。すなわち、連続ＮＣＣは、初期的に「０．８」弱を示し、その後「０．２」を下回り、約０．０２秒において「０．７」まで復帰する。そして、連続ＮＣＣは、「０．６」程度まで低下し、次いで「０．７」程度まで復帰する。その後、連続ＮＣＣは再び「０．５」程度まで低下した後に、約０．１秒において「０．６」弱となる。以後、連続ＮＣＣは「０．６５」付近までゆっくりと上昇していく。
【００７３】
このようにガス器具１０の使用開始時において、連続ＮＣＣは大半の期間で「０．９」以上を示さない。このため、判断部５２ａは、連続ＮＣＣの代表値が閾値以上でない場合、ガス器具１０の使用開始であると判断する。
【００７４】
また、連続ＮＣＣはガス器具１０毎に異なっている。このため、判断部５２ａは、このような連続ＮＣＣのパターンから使用が開始したガス器具１０を判断する。具体的には図４に示す類似度推移パターン記憶部５３ａに、各ガス器具１０の連続ＮＣＣのパターンを記憶させておく。すなわち、類似度推移パターン記憶部５３ａは、ガステーブルについて図９（ａ）に示したような連続ＮＣＣのパターンを記憶し、小型湯沸器について図９（ｂ）に示したような連続ＮＣＣのパターンを記憶し、給湯器について図９（ｃ）に示したような連続ＮＣＣのパターンを記憶している。そして、判断部５２ａは、類似度推移パターン記憶部５３ａにより記憶された連続ＮＣＣデータのうち、算出された連続ＮＣＣと最も近い連続ＮＣＣデータが示すガス器具１０の使用が開始したと判断する。
【００７５】
図１０は、ガス器具使用終了時における圧力変化を示すグラフであって、（ａ）はガステーブル使用終了時における圧力変化を示し、（ｂ）は小型湯沸器使用終了時における圧力変化を示し、（ｃ）は給湯器使用終了時における圧力変化を示している。
【００７６】
図１０（ａ）に示すように、ガステーブルの使用終了時には圧力が２．８５ｋＰａ程度で滑らかに振動する圧力波形が得られる。また、図１０（ｂ）に示すように、小型湯沸器の使用終了時には圧力が２．８５ｋＰａを基準にして０．１ｋＰａ程度振動する圧力波形が得られる。さらに、図１０（ｃ）に示すように、給湯器の使用終了時には圧力が２．８８ｋＰａを基準にしてガステーブルよりもやや粗い振動を示す圧力波形が得られる。
【００７７】
図１１は、図４に示した類似度推移算出部５２ｃにより算出される連続ＮＣＣを示すグラフであって、（ａ）はガステーブル使用終了時における連続ＮＣＣを示し、（ｂ）は小型湯沸器使用終了時における連続ＮＣＣを示し、（ｃ）は給湯器使用終了時における連続ＮＣＣを示している。
【００７８】
ガステーブルの使用が終了した場合、図１０（ａ）の振動波形が得られ、類似度推移算出部５２ｃにより算出される連続ＮＣＣは図１１（ａ）に示すようになる。すなわち、連続ＮＣＣは、初期的に「０．９」程度を示し、その後「０．２」を下回り、約０．０３秒において「０．８」まで復帰する。そして、連続ＮＣＣは、約０．１秒において「０．６」程度となり、その後「０．６」付近を維持する。
【００７９】
また、小型湯沸器の使用が終了した場合、図１０（ｂ）の振動波形が得られ、類似度推移算出部５２ｃにより算出される連続ＮＣＣは図１１（ｂ）に示すようになる。すなわち、連続ＮＣＣは、初期的に「０．９」程度を示し、その後「０．２」を下回り、約０．０１秒において「０．８」まで復帰する。そして、連続ＮＣＣは、再度「０．３」程度まで低下し、次いで「０．６」程度まで復帰する。その後、連続ＮＣＣは小さな振動を繰り返しながら約０．１秒において「０．６」程度となる。次に、連続ＮＣＣは「０．７」付近までゆっくりと上昇する。
【００８０】
また、給湯器の使用が終了した場合、図１０（ｃ）の振動波形が得られ、類似度推移算出部５２ｃにより算出される連続ＮＣＣは図１１（ｃ）に示すようになる。すなわち、連続ＮＣＣは、初期的に「０．９」程度を示し、その後「０．２」を下回り、約０．０２秒において「０．６」まで復帰する。そして、連続ＮＣＣは、再度「０．４５」程度まで低下し、次いで「０．６」程度まで復帰する。その後、連続ＮＣＣは再び「０．４５」程度まで低下した後に、約０．１秒において「０．６」程度となる。次に、連続ＮＣＣは「０．７」付近までゆっくりと上昇する。
【００８１】
このように、ガス器具１０の使用終了時においても連続ＮＣＣはガス器具１０毎に異なり、判断部５２ａは、このような連続ＮＣＣのパターンから使用が終了したガス器具１０を判断する。具体的には図４に示す類似度推移パターン記憶部５３ａに、各ガス器具１０の連続ＮＣＣのパターンを記憶させておく。すなわち、類似度推移パターン記憶部５３ａは、ガステーブルについて図１１（ａ）に示したような連続ＮＣＣのパターンを記憶し、小型湯沸器について図１１（ｂ）に示したような連続ＮＣＣのパターンを記憶し、給湯器について図１１（ｃ）に示したような連続ＮＣＣのパターンを記憶している。そして、判断部５２ａは、類似度推移パターン記憶部５３ａにより記憶された連続ＮＣＣデータのうち、算出された連続ＮＣＣと最も近い連続ＮＣＣデータが示すガス器具１０の使用が終了したと判断する。
【００８２】
次に、生成部５２ｂによって生成される所定の振動波形の生成手法について説明する。生成部５２ｂは、以下のようにしてガス漏れ振動波形を生成する。まず、生成部５２ｂは以下の式（２）を記憶している。
【数２】
0005436175.tif 000003
【００８３】
ここで、ｙ（ｔ）は圧力の変化量を示し、Ｋはゲインを示し、ωｄは減衰振動の周波数を示し、ζは減衰比を示している。特に、ゲインＫ、減衰振動の周波数ωｄ、及び減衰比ζは、圧力センサ４１によって実際に計測された波形から求められるものである。次に、これらの算出方法について図６を参照して説明する。
【００８４】
生成部５２ｂは、以下の式（３）から、減衰振動の周波数ωｄを算出する。
【数３】
0005436175.tif 000004
【００８５】
ここで、Ｔｐは行き過ぎ時間であり、図６で示すように、圧力変化発生時から最初の極値Ｖ１（極小値Ｖ１）までの時間をいう。生成部５２ｂは、計測値データから最初の極値Ｖ１が確認されると、行き過ぎ時間Ｔｐを求め、式（３）から減衰振動の周波数ωｄを算出する。
【００８６】
なお、減衰振動の周波数ωｄは、式（３）から求める場合に限らず、圧力変化発生時から２つ目の極値Ｍ（極大点Ｍ）や、３つ目の極値Ｖ２（極小点Ｖ２）に基づいて算出してもよい。
【００８７】
次に、生成部５２ｂは、以下の式（４）から、ゲインＫを算出する。
【数４】
0005436175.tif 000005
【００８８】
このような式であるため、生成部５２ｂは、計測値データから極値Ｖ１，Ｍ，Ｖ２が確認されると、式（４）からゲインＫを算出する。
【００８９】
なお、図６から明らかなように、ゲインＫは圧力変化発生前の圧力値と圧力変化発生後の圧力値との差分によっても求めることができる。従って、生成部５２ｂは、圧力変化が発生して圧力値が略一定値となったとき（図６では時刻０．４秒）に、差分からゲインＫを求めてもよい。さらに、類似度推移算出部５２ｃは、圧力変化発生時から４つ目以降の極値を加味してゲインＫを算出してもよい。
【００９０】
次いで、生成部５２ｂは、以下の式（５）から、減衰比ζを算出する。
【数５】
0005436175.tif 000006
【００９１】
ここで、δは対数減衰率であり、ｍは周期数である。式（５）の場合、周期数ｍは「０．５」となる。
【００９２】
このような式であるため、生成部５２ｂは、計測値データから極値Ｖ１，Ｍが確認されると、式（５）から減衰比ζを算出する。
【００９３】
以上のように、生成部５２ｂは、ゲインＫ、減衰振動の周波数ωｄ、及び減衰比ζを算出し、式（２）より振動波形の式を求める。そして、類似度推移算出部５２ｃは、求めた式と、計測値データ（圧力波形）とから、式（１）に従って連続ＮＣＣを求めることとなる。
【００９４】
ここで、生成部５２ｂは、減衰振動の周波数ωｄ、及び減衰比ζを以下のようにして算出するようにしてもよい。すなわち、図６に示す振動波形は、ガス漏れ時の流量に依存する傾向にある。このため、生成部５２ｂは、流量値のみを変数に含む式を予め記憶し、この式に流量値を代入して、減衰振動の周波数ωｄ、及び減衰比ζを求めるようにしてもよい。
【００９５】
具体的に生成部５２ｂは、以下の式（６）から減衰振動の周波数ωｄ、及び減衰比ζを求める。
【数６】
0005436175.tif 000007
【００９６】
ここで、Ｌは流量値であり、ａ１，ａ２，ｂ１，ｂ２は定数である。このように、式（６）から求めることで演算量を減らして、算出処理の簡素化を図るようにしてもよい。なお、流量と圧力には一定の相関がある。このため、式（６）に代えて圧力値のみを変数に含む式を記憶し、この式から減衰振動の周波数ωｄ、及び減衰比ζを求めるようにしてもよい。
【００９７】
さらに、この場合、生成部５２ｂは、ゲインＫについて式（４）から算出することなく、圧力変化発生前の圧力値と圧力変化発生後の圧力値との差分によっても求めることが望ましい。これにより、一層演算量を減らすことができるからである。
【００９８】
再度、図４を参照する。保安部５２ｄは、リコール対象となるガス器具１０に対する保安処理を実行するものである。ここで、上記した判断部５２ａでは、連続ＮＣＣに基づくことにより、ガステーブルや給湯器などのガス器具１０の種類のみならず、メーカ名や型番についても判断することができる。よって、リコール器具記憶部５３ｂにリコール対象となるガス器具１０のメーカ名及び型番などを記憶させ、使用ガス器具１０が判断されたときにリコール器具記憶部５３ｂの記憶内容を参照することにより、リコール対象となるガス器具１０であるか否かを判断することができる。
【００９９】
なお、保安部５２ｄによる保安処理とは、ガスメータ４０側に遮断弁を弁閉する信号を送信してガスの使用を禁止したり、外部業者に連絡を行い作業員を派遣したりする処理である。また、保安部５２ｇは、ガス漏れ発生時においてガスメータ４０側に遮断弁を弁閉する信号を送信してガスの使用を禁止してもよい。
【０１００】
なお、リコール器具記憶部５３ｂは、リコール対象となるガス器具１０のメーカ名及び型番のみを記憶する場合に限らず、さらに、リコール対象となるガス器具１０の波形データ、類似度推移を記憶しておき、メーカ名、型番、波形データ及び類似度推移を１グループとしてまとめて記憶していてもよい。また、判断部５２ａは、リコール対象となるガス器具１０を判断するにあたり、メーカ名、型番、波形データ及び類似度推移のうち、どの項目をインデックス（ＩＤ）として呼び出してもよい。特に、波形データや類似度推移を呼び出した場合、判断部５２ａは、波形データ同士や類似度推移同士を比較し、リコール対象となるガス器具１０を判断することで、使用ガス器具１０のメーカ名や型番を判断しなくともよい。
【０１０１】
次に、フローチャートを参照して、本実施形態に係るガスメータ４０の動作を説明する。図１２は、本実施形態に係るガスメータ４０の動作の一例を示すフローチャートである。
【０１０２】
まず、図１２に示すように、制御部４３はトリガ信号が発生したか否かを判断する（Ｓ１）。トリガ信号が発生していないと判断した場合（Ｓ１：ＮＯ）、発生したと判断されるまで、この処理が繰り返される。一方、トリガ信号が発生したと判断した場合（Ｓ１：ＹＥＳ）、サンプリング時間調整部４３ａは、サンプリング時間を通常のサンプリング時間よりも短縮し、短縮されたサンプリング時間で圧力を計測する（Ｓ２）。その後、制御部４３は、微小時間経過したか否かを判断する（Ｓ３）。微小時間経過していないと判断した場合（Ｓ３：ＮＯ）、処理はステップＳ２に移行する。なお、ステップＳ２では、圧力のサンプリング時間を短縮しているが、圧力のサンプリング時間に代えて、流量のサンプリング時間を短縮するようにしてもよい。
【０１０３】
微小時間経過したと判断した場合（Ｓ３：ＹＥＳ）、送受信部４４は、ステップＳ２において計測された計測値データと、第１，第２トリガ信号のいずれの信号が発生されたかを示す情報とを一括で送信する（Ｓ４）。その後、図１２に示す処理は終了する。
【０１０４】
なお、図１２に示す処理の終了時においてサンプリング時間調整部４３ａは、サンプリング時間を通常のサンプリング時間に戻したり、通常のサンプリング時間よりも長く設定したりする。例えば、ステップＳ１において第１トリガ信号が発生された場合、流量等の計測の必要があることから、サンプリング時間調整部４３ａは、サンプリング時間を通常のサンプリング時間にする。一方、ステップＳ１において第２トリガ信号が発生された場合、流量値がゼロになっていれば、サンプリング時間を通常のサンプリング時間よりも長く設定し、流量値がゼロを超える値であれば、サンプリング時間を通常のサンプリング時間にする。このように、サンプリング時間調整部４３ａは、微小時間経過後に、サンプリング時間の短縮状態を解除する。
【０１０５】
図１３は、本実施形態に係る管理センター５０の動作の一例を示すメインフローチャートである。図１３に示すように、制御部５２は計測値データを受信したか否かを判断する（Ｓ１１）。計測値データを受信していないと判断した場合（Ｓ１１：ＮＯ）、受信したと判断されるまで、この処理が繰り返される。
【０１０６】
計測値データを受信したと判断した場合（Ｓ１１：ＹＥＳ）、制御部５２は、計測値データと共に一括送信されたトリガ信号を示す情報に従って、第１トリガ信号が発生されていたか否かを判断する（Ｓ１２）。第１トリガ信号が発生されていなかったと判断した場合（Ｓ１２：ＮＯ）、第１トリガ信号が発生されていたと判断されるまで、この処理が繰り返される。
【０１０７】
一方、第１トリガ信号が発生されていたと判断した場合（Ｓ１２：ＹＥＳ）、制御部５２は、ガス漏れ／開始ガス器具判断処理を実行する（Ｓ１３）。第１トリガ信号が発生したということは、ガス流量が増加したりガス圧力が減少したりした場合であり、ガス漏れの発生かガス器具１０の使用であると判断できるため、制御部５２は、ガス漏れ／開始ガス器具判断処理を実行する。
【０１０８】
その後、保安部５２ｄは、ステップＳ１３の処理においてガス漏れが発生していたか否かを判断する（Ｓ１４）。ガス漏れが発生していたと判断した場合（Ｓ１４：ＹＥＳ）、保安部５２ｄは保安処理を実行し（Ｓ１５）、図１３に示す処理は終了する。
【０１０９】
また、ガス漏れが発生していなかったと判断した場合（Ｓ１４：ＮＯ）、制御部５２は、計測値データと共に一括送信されたトリガ信号を示す情報に従って、第２トリガ信号が発生されていたか否かを判断する（Ｓ１６）。第２トリガ信号が発生されていなかったと判断した場合（Ｓ１６：ＮＯ）、処理はステップＳ１４に移行する。
【０１１０】
一方、第２トリガ信号が発生されていたと判断した場合（Ｓ１６：ＹＥＳ）、制御部５２は、終了ガス器具判断処理を実行する（Ｓ１７）。第２トリガ信号が発生したということは、ガス流量が減少したりガス圧力が上昇したりした場合であり、ガス器具１０の使用終了時であると判断できるため、ガス漏れの発生ではないといえる。よって、制御部５２は、終了ガス器具判断処理を実行する。
【０１１１】
そして、保安部５２ｄは、リコールガス器具１０の使用があったか否かを判断する（Ｓ１８）。このとき、保安部５２ｄは、ステップＳ１３におけるガス漏れ／開始ガス器具判断処理、及び、ステップＳ１７における終了ガス器具判断処理の少なくとも一方においてリコールガス器具１０の使用があったか否かを判断する。リコールガス器具１０の使用があったと判断した場合（Ｓ１８：ＹＥＳ）、保安部５２ｄは保安処理を実行し（Ｓ１５）、図１３に示す処理は終了する。一方、リコールガス器具１０の使用でなかったと判断した場合（Ｓ１８：ＮＯ）、保安部５２ｄは保安処理を実行せず、図１３に示す処理は終了する。
【０１１２】
図１４は、図１３に示したガス漏れ／開始ガス器具判断処理（Ｓ１３）の詳細を示すフローチャートである。図１４に示すように、まず、生成部５２ｂは、微小時間中に得られた振動波形から、減衰振動の周波数ωｄ、ゲインＫ、及び減衰比ζを決定する（Ｓ２１）。このとき、生成部５２ｂは、減衰振動の周波数ωｄ、ゲインＫ、及び減衰比ζを式（３）〜式（５）に基づいて算出してもよいし、式（６）から求めてもよい。
【０１１３】
次に、生成部５２ｂは、ステップＳ２１により決定された減衰振動の周波数ωｄ、ゲインＫ、及び減衰比ζから、２次遅れのステップ応答の式に基づいてガス漏れ振動波形を生成する（Ｓ２２）。このとき、生成部５２ｂは、ステップＳ２１により決定された減衰振動の周波数ωｄ、ゲインＫ、及び減衰比ζを式（２）に代入することにより、ガス漏れ振動波形を生成する。
【０１１４】
そして、判断部５２ａは、ステップＳ２２において生成されたガス漏れ振動波形と、受信した計測値データからなる振動波形とに基づいて、式（１）から連続ＮＣＣを算出する（Ｓ２３）。
【０１１５】
次に、判断部５２は、連続ＮＣＣの代表値を決定し、代表値が閾値以上であるか否かを判断する（Ｓ２４）。代表値が閾値以上であると判断した場合（Ｓ２４：ＹＥＳ）、判断部５２は、ガス漏れが発生していると判断する（Ｓ２６）。その後、図１４に示す処理は終了する。
【０１１６】
代表値が閾値以上でないと判断した場合（Ｓ２４：ＮＯ）、判断部５２ａは、類似度推移パターン記憶部５３ａからガス器具１０毎の類似度推移データを読み出す（Ｓ２６）。次いで、判断部５２ａは、ステップＳ２６にて読み出したガス器具１０毎の連続ＮＣＣデータのうち、ステップＳ２３において算出した連続ＮＣＣと最も近いものを特定し、使用が開始したガス器具１０を判断する（Ｓ２７）。その後、図１４に示す処理は終了する。
【０１１７】
なお、図１４に示す処理では、ステップＳ２７において使用が開始したガス器具１０を判断するのに先立って、ステップＳ２４においてガス漏れを判断することにより、迅速性を必要とするガス漏れの判断を優先し、安全性の向上を図っている。
【０１１８】
図１５は、図１３に示した終了ガス器具判断処理（Ｓ１７）の詳細を示すフローチャートである。図１５に示すステップＳ３１〜Ｓ３３において、図１４に示したステップＳ２１〜Ｓ２３と同様の処理が実行される。
【０１１９】
その後、図１５に示すステップＳ３４，Ｓ３５において、図１４に示したステップＳ２６，Ｓ２７と同様の処理が実行される。そして、図１５に示す処理は終了する。
【０１２０】
このようにして、本実施形態に係るガス使用状況判断システム１及びガス使用状況判断方法によれば、ガスメータ４０は計測値データを送信し、管理センター５０は、受信した計測値データの所定以上の変化時からの微小時間中に得られる振動波形からガス漏れ及び使用ガス器具１０の少なくとも一方を判断する。ここで、本件出願人らは、ガス漏れ発生時や、ガス器具１０の使用開始時及び終了時など、計測値データの所定以上の変化時から微小時間経過するまでの波形に、ガス漏れやガス器具１０毎に固有の振動を示すことを見出した。このため、所定以上の変化時からの微小時間中に得られる振動波形からガス漏れ及び使用ガス器具１０の少なくとも一方を判断することができる。特に、ガスメータ４０はセンサ４１，４２により出力された信号に基づく計測値データを送信するだけでよく、判断についての演算を行う必要がない。このため、基本的に電池駆動となるガスメータ４０の消費電力を抑えることができる。従って、消費電力を抑えつつ、ガス漏れ及び使用ガス器具１０の少なくとも一方を判断することができる。
【０１２１】
また、センサ４１，４２により出力された信号の所定以上の変化時にトリガ信号を発生し、トリガ信号が発生された場合に、センサ４１，４２のサンプリング時間を通常のサンプリング時間よりも短縮する。このため、微小時間における波形の振動を適切にとらえることができ、適切に計測することができる。
【０１２２】
また、微小時間経過後に、センサ４１，４２により出力された信号に基づく計測値データを、微小時間分一括して送信するため、計測の都度に送信する場合と比較してガスメータ４０による送信回数の増加を抑えることができ、一層消費電力を抑えることができる。
【０１２３】
また、微小時間だけサンプリング時間を短縮した後、サンプリング時間の短縮状態を解除するため、ガス漏れや使用ガス器具１０の判断に必要なる振動波形が得られた後は、詳細なサンプリングを必要とせず、そのような場合にサンプリング時間の短縮状態を解除することで、一層消費電力を抑えることができる。
【０１２４】
また、トリガ信号が発生されるまでは、センサ４１，４２のサンプリング時間を通常のサンプリング時間よりも長くする。ここで、トリガ信号が発生していない場合とは、流量や圧力に所定以上の変化がなく、流量や圧力が安定している場合といえる。すなわち、計測の必要性が少ない場合といえる。よって、このような場合に、サンプリング時間を長くしておくことで、一層消費電力を軽減させることができる。
【０１２５】
また、計測値データに基づく振動波形と所定の振動波形との類似度推移を算出し、算出された類似度推移に基づいて、ガス漏れ及び使用ガス器具１０の少なくとも一方を判断する。ここで、所定の振動波形が特定のガス器具１０の振動波形であるとする。そして、センサ４１，４２から得られた信号に基づく振動波形が特定のガス器具１０の振動波形である場合、算出される類似度は高くなる傾向にあり、類似度推移についても全体的に高い値を示す傾向にある。また、センサ４１，４２から得られた信号に基づく振動波形が他のガス器具１０の振動波形である場合には、類似度は高くならず、類似度推移についても高くならないが、類似度推移はそのガス器具１０特有の推移を示す。ガス漏れ時も同様にガス漏れ特有の推移を示す。従って、上記のように類似度推移に基づいて判断を実行することで、振動波形の特徴をとらえて、ガス漏れ及び使用ガス器具１０の少なくとも一方を判断することができる。なお、上記実施形態では所定の振動波形を生成しているが、これに限らず、予め管理センター５０のハードディスクなどにデータテーブルとして記憶されていてもよい。また、これらの組み合わせであってもよい。
【０１２６】
また、ガス漏れ及び使用が開始されたガス器具１０の双方を判断可能であって、使用が開始されたガス器具１０を判断するのに先立って、ガス漏れ判断の処理を実行する。このため、緊急性の高いガス漏れ判断を優先して行うこととなり、安全性の向上を図ることができる。
【０１２７】
また、管理センター５０はリコール対象となるガス器具１０の情報を記憶し、使用ガス器具１０がリコール対象のガス器具１０である場合、当該リコール対象のガス器具１０に対する保安処理を実行させる。このため、リコール品の回収や使用禁止を行うことができ、安全性を向上させることができる。加えて、リコール品を管理センター５０側で登録しておけばよく、各需要者側において登録等する必要がないことから、リコール品に対する保安処理についての工数を低減することができる。
【０１２８】
次に、本発明の第２実施形態を説明する。第２実施形態に係るガス使用状況判断システム１及びガス使用状況判断方法は、第１実施形態のものと同様であるが、構成及び処理内容が一部異なっている。以下、第１実施形態との相違点を説明する。
【０１２９】
図１６は、第２実施形態に係る管理センター５０の詳細を示す構成図である。図１６に示すように、第２実施形態に係る管理センター５０は、生成部５２ｂ及び類似度推移算出部５２ｃに代えて、解析部（解析手段）５２ｅを備えている。
【０１３０】
解析部５２ｅは、計測値データに基づく振動波形を解析して、周波数と振幅との相関を示すスペクトルデータを算出するものである。具体的に本実施形態に係る解析部５２ｅは、振動波形をフーリエ変換することにより、スペクトルデータを算出する。なお、解析部５２ｅはフーリエ変換によりスペクトルデータを算出する場合に限らず、他の方法によってスペクトルデータを算出するようにしてもよい。
【０１３１】
図１７は、図１６に示した解析部５２ｅにより算出されるスペクトルデータを示すグラフであって、ガス漏れが発生したときの圧力波形をフーリエ変換して得られるスペクトルデータを示すグラフである。図１７に示すように、ガス漏れが発生した場合、得られる圧力波形には２０Ｈｚ以上の周波数成分が殆ど含まれていない。なお、６０Ｈｚ付近において存在するピークは、商用電源によるノイズであると考えられる。
【０１３２】
図１８は、図１６に示した解析部５２ｅにより算出されるスペクトルデータを示すグラフであって、（ａ）はガステーブル使用開始時におけるスペクトルデータを示し、（ｂ）は小型湯沸器使用開始時におけるスペクトルデータを示し、（ｃ）は給湯器使用開始時におけるスペクトルデータを示している。
【０１３３】
図１８（ａ）に示すように、ガステーブルの使用が開始した場合、得られる圧力波形には３０Ｈｚ以下の周波数成分が多く、特に１０〜２０Ｈｚ付近において大きな振幅を示す傾向がある。また、図１８（ｂ）に示すように、小型湯沸器の使用が開始した場合、圧力波形は１５０Ｈｚまでの圧力成分を含んでおり、特に３０Ｈｚ程度では非常に大きな振幅を示す傾向がある。さらに、図１８（ｃ）に示すように、給湯器の使用が開始した場合、圧力波形は１８０Ｈｚまでの圧力成分を含んでおり、特に２０Ｈｚ程度では非常に大きな振幅を示す傾向がある。
【０１３４】
図１９は、図１６に示した解析部５２ｅにより算出されるスペクトルデータを示すグラフであって、（ａ）はガステーブル使用終了時におけるスペクトルデータを示し、（ｂ）は小型湯沸器使用終了時におけるスペクトルデータを示し、（ｃ）は給湯器使用終了時におけるスペクトルデータを示している。
【０１３５】
図１９（ａ）に示すように、ガステーブルの使用が終了した場合、得られる圧力波形には３０Ｈｚ以下の周波数成分が多く、特に１０〜２０Ｈｚ付近において大きな振幅を示す傾向がある。また、図１９（ｂ）に示すように、小型湯沸器の使用が終了した場合、圧力波形は１５０Ｈｚまでの圧力成分を含んでおり、特に９０Ｈｚ程度で大きな振幅を示す傾向がある。さらに、図１９（ｃ）に示すように、給湯器の使用が終了した場合、３０Ｈｚ程度でやや大きな振幅を示す程度であり、その他の周波数成分を殆ど含まない傾向がある、なお、５０Ｈｚ付近において存在するピークは、商用電源によるノイズであると考えられる。
【０１３６】
再度、図１６を参照する。記憶部５３は、類似度推移パターン記憶部５３ａに代えて、スペクトルデータ記憶部５３ｃを記憶している。このスペクトルデータ記憶部５３ｃは、図１７〜図１９に示したようなスペクトルデータを記憶している。そして、判断部５２ａは、このスペクトルデータに基づいて、ガス漏れ及び使用ガス器具１０を判断する。
【０１３７】
すなわち、解析部５２ｅは計測値データに基づく振動波形をフーリエ変換してスペクトルデータを算出する。判断部５２ａは、解析部５２ｅにより算出されたスペクトルデータと、スペクトルデータ記憶部５３ｃにより記憶されたスペクトルデータとを比較し、類似度が最も高いスペクトルデータを特定し、ガス漏れの発生や使用ガス器具１０について判断する。ここで、類似度とは、上記したＮＣＣであってもよいし、他の手法により算出された類似度であってもよい。
【０１３８】
なお、第１実施形態においてリコール器具記憶部５３ｂは、少なくともリコール対象となるガス器具１０のメーカ名及び型番を記憶し、望ましくは波形データ及び類似度推移を記憶していた。第２実施形態ではスペクトルデータによりガス器具１０を判断するため、リコール器具記憶部５３ｂは、類似度推移に代えてスペクトルデータを記憶することとなる。
【０１３９】
次に、フローチャートを参照して、第２実施形態に係るガスメータ４０の動作を説明する。なお、第２実施形態においてガスメータ４０のフロー及び管理センター５０のメインフローは第１実施形態と同様であるため、説明を省略する。
【０１４０】
図２０は、図１３に示したガス漏れ／開始ガス器具判断処理（Ｓ１３）の詳細を示す第２のフローチャートである。図２０に示すように、まず、解析部５２ｅは、ガスメータ４０から受信した計測値データに基づく振動波形をフーリエ変換し、スペクトルデータを算出する（Ｓ４１）。その後、判断部５２ａは、ガス漏れのスペクトルデータを読み出し（Ｓ４２）、読み出したガス漏れのスペクトルデータと、ステップＳ４１にて算出したスペクトルデータとの類似度を算出する（Ｓ４３）。
【０１４１】
次に、判断部５２ａは、ステップＳ４３にて算出した類似度が特定値以上であるか否かを判断する（Ｓ４４）。ステップＳ４３にて算出した類似度が特定値以上であると判断した場合（Ｓ４４：ＹＥＳ）、判断部５２ａはガス漏れが発生したと判断する（Ｓ４５）。そして、図２０に示す処理は終了する。
【０１４２】
ところで、ステップＳ４３にて算出した類似度が特定値以上でないと判断した場合（Ｓ４４：ＮＯ）、判断部５２ａは、ガス器具１０毎のスペクトルデータを読み出し（Ｓ４６）、読み出したガス器具１０毎のスペクトルデータと、ステップＳ４１にて算出したスペクトルデータとの類似度を算出する（Ｓ４７）。
【０１４３】
その後、判断部５２ａは、類似度が最大となったスペクトルデータが示す種類のガス器具１０の使用が開始したと判断する（Ｓ４８）。そして、図２０に示す処理は終了する。
【０１４４】
図２１は、図１３に示した終了ガス器具判断処理（Ｓ１７）の詳細を示す第２のフローチャートである。図２１に示すステップＳ５１において、図２０に示したステップＳ４１と同様の処理が実行される。
【０１４５】
その後、図２１に示すステップＳ５２〜Ｓ５４において、図２０に示したステップＳ４６〜Ｓ４８と同様の処理が実行される。そして、図２１に示す処理は終了する。
【０１４６】
このようにして、第２実施形態に係るガス使用状況判断システム１及びガス使用状況判断方法によれば、第１実施形態と同様に、消費電力を抑えつつ、ガス漏れ及び使用ガス器具１０の少なくとも一方を判断することができる。
【０１４７】
さらに、第２実施形態によれば、計測値データに基づく振動波形を解析して、周波数と振幅との相関を示すスペクトルデータを算出し、算出されたスペクトルデータに基づいて、ガス漏れ及び使用ガス器具１０の少なくとも一方を判断する。ここで、ガス器具使用時及び終了時やガス漏れ時における振動波形それぞれには、圧力や流量の波形の周波数及び振幅に特徴があらわれる。よって、波形を解析して周波数と振幅との相関を示すスペクトルデータを得ると共に、このスペクトルデータに基づいて判断を実行することで、ガス漏れ及び使用ガス器具の少なくとも一方を判断することができる。
【０１４８】
以上、実施形態に基づき本発明を説明したが、本発明は上記実施形態に限られるものではなく、本発明の趣旨を逸脱しない範囲で、変更を加えてもよい。
【０１４９】
また、本実施形態においてガスメータ４０を計測装置の一例としているが、これに限らず、計測装置はガスメータ４０でなくともよい。さらには、既存のガスメータに送信機など機能を追加して、計測装置を構成してもよい。
【０１５０】
さらに、第１実施形態において類似度推移を式（１）により算出しているが、これに限らず、他の方法で類似度推移を算出するようにしてもよい。
【０１５１】
また、第１実施形態において判断部５２ａは、類似度推移パターン記憶部５３ａに記憶された連続ＮＣＣデータのうち、判断部５２ａにより算出された連続ＮＣＣと近いものが存在しない場合、類似度推移パターン記憶部５３ａに記憶された連続ＮＣＣデータが示すガス器具１０に不足があると判断してもよい。
【０１５２】
また、本実施形態では燃料ガスをＬＰガスとする場合の例について説明したが、これに限らず、都市ガスの場合にも適用可能である。
【０１５３】
また、本実施形態において微小時間を最大で２秒（望ましくは１秒以内）としているが、２秒よりも長い時間であってもよい。
【０１５４】
また、本実施形態において類似度推移パターン記憶部５３ａは、ガス器具１０毎の連続ＮＣＣデータを記憶している。この連続ＮＣＣデータは、１つのガス器具１０に対して１つだけ記憶されていてもよいし、１つのガス器具１０に対して複数記憶されていてもよい。例えば、給湯器では給湯器内の水温によって連続ＮＣＣが異なってくる。この場合、類似度推移パターン記憶部５３ａに記憶される連続ＮＣＣデータが１つだけであると、給湯器の水温に応じて使用ガス器具１０の判断を誤ってしまう可能性がある。そこで、このようなガス器具１０に対しては複数の連続ＮＣＣデータを記憶しておくことが望ましい。これにより、より精度良く使用ガス器具１０を判断することができるからである。
【０１５５】
また、第２実施形態において解析部５２ｅは、スペクトルデータの全周波数域で類似度を算出しているが、これに限らず一部の周波数域のみで類似度を算出してもよい。例えば、給湯器の使用終了時では１００Ｈｚ以上の周波数域においてもスペクトルデータに大きな振幅が得られるという特徴があるため、１００Ｈｚ以上の周波数域についてスペクトルデータの類似度を算出することによっても使用が終了したガス器具１０を特定することができる。このように、一部の周波数域のみで類似度を算出して演算量を減らすこともできる。
【０１５６】
また、上記実施形態では、使用が開始したガス器具１０、使用が終了したガス器具１０、及び、ガス漏れについて、連続ＮＣＣを求めたり、スペクトルデータを求めたりすることで、判断している。しかし、これに限らず、例えば、図６、図８及び図１０に示すような微小時間における波形を直接記憶しておき、波形同士の類似度などから、使用が開始したガス器具１０、使用が終了したガス器具１０、及び、ガス漏れを判断するようにしてもよい。さらには、波形の特定点など波形の直接の特徴から使用が開始したガス器具１０、使用が終了したガス器具１０、及び、ガス漏れを判断するようにしてもよい。
【０１５７】
また、上記実施形態において送受信部４４は圧力と流量との双方の計測値データを送信しているが、これに限らず、いずれか一方のみの計測値データを送信するようにしてもよい。
【０１５８】
さらに、上記実施形態においてセンターは、所定の振動波形を生成しているが、これに限らず、予めデータテーブルとして所定の振動波形を記憶しておいてもよいし、これらを組み合わせてもよい。
【符号の説明】
【０１５９】
１…ガス使用状況判断システム
１０…ガス器具
２０…調整器
３１…第１配管
３２…第２配管
４０…ガスメータ（ガス使用状況判断システム）
４１…圧力センサ
４２…流量センサ
４３…制御部
４３ａ…サンプリング時間調整部（サンプリング時間調整手段）
４４…送受信部（送信手段）
４５…トリガ信号発生部（トリガ信号発生手段）
４３ｈ…解析部（解析手段）
５０…管理センター（センター）
５１…送受信部（受信手段）
５２…制御部
５２ａ…判断部（判断手段）
５２ｂ…生成部
５２ｃ…類似度推移算出部（類似度推移算出手段）
５２ｄ…保安部（保安手段）
５２ｅ…解析部（解析手段）
５３…記憶部
５３ａ…類似度推移パターン記憶部
５３ｂ…リコール器具記憶部（リコール器具記憶手段）
５３ｃ…スペクトルデータ記憶部
【図１】
0005436175.tif 000008
【図２】
0005436175.tif 000009
【図３】
0005436175.tif 000010
【図４】
0005436175.tif 000011
【図５】
0005436175.tif 000012
【図６】
0005436175.tif 000013
【図７】
0005436175.tif 000014
【図８】
0005436175.tif 000015
【図９】
0005436175.tif 000016
【図１０】
0005436175.tif 000017
【図１１】
0005436175.tif 000018
【図１２】
0005436175.tif 000019
【図１３】
0005436175.tif 000020
【図１４】
0005436175.tif 000021
【図１５】
0005436175.tif 000022
【図１６】
0005436175.tif 000023
【図１７】
0005436175.tif 000024
【図１８】
0005436175.tif 000025
【図１９】
0005436175.tif 000026
【図２０】
0005436175.tif 000027
【図２１】
0005436175.tif 000028