JP2020201268A - 居住者の不在期間の確定方法、装置、及びプログラム - Google Patents

居住者の不在期間の確定方法、装置、及びプログラム Download PDF

Info

Publication number
JP2020201268A
JP2020201268A JP2020131315A JP2020131315A JP2020201268A JP 2020201268 A JP2020201268 A JP 2020201268A JP 2020131315 A JP2020131315 A JP 2020131315A JP 2020131315 A JP2020131315 A JP 2020131315A JP 2020201268 A JP2020201268 A JP 2020201268A
Authority
JP
Japan
Prior art keywords
time
measurement
gas
data
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2020131315A
Other languages
English (en)
Other versions
JP6938735B2 (ja
Inventor
篠原　靖志
Yasushi Shinohara
靖志 篠原
剛 上野
Takeshi Ueno
剛 上野
俊一 服部
Shunichi Hattori
俊一 服部
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Central Research Institute of Electric Power Industry
Original Assignee
Central Research Institute of Electric Power Industry
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Central Research Institute of Electric Power Industry filed Critical Central Research Institute of Electric Power Industry
Publication of JP2020201268A publication Critical patent/JP2020201268A/ja
Application granted granted Critical
Publication of JP6938735B2 publication Critical patent/JP6938735B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measuring Volume Flow (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Feeding And Controlling Fuel (AREA)
  • Details Of Flowmeters (AREA)

Abstract

【課題】関連する情報やデータを必要とすることなく、対象住宅や対象住戸のガス総消費量の計測データに基づいて家庭におけるガス消費量を用途別に良好な精度で分解する方法を提供する。【解決手段】ガス配管の主幹における計測によって所定の時間間隔の計測時刻と、計測時刻毎に出力される値であって所定の流量毎に１カウントとして計数する流量カウント数との組み合わせデータとして取得されるガス総消費量の計測データからガス器具が使用されている時間帯の組み合わせデータが連続動作区間のデータとして抽出され、連続動作区間のデータにおける計測時刻の時系列順に付与される計測時刻毎の流量カウント数と流量とを掛け合わせた値が算出され、計測時刻毎の値が用いられて計測時刻毎の主幹におけるガスの流速が算定された上で、ガス総消費量の用途が給湯分、暖房分、及び厨房分に分解されるようにした。【選択図】図１

Description

本発明は、ガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムに関する。さらに詳述すると、本発明は、家庭におけるガスの用途毎の需要量の把握に用いて好適な技術に関する。また、本発明は、居住者の在宅不在の判定方法、判定装置、及び判定プログラムに関する。さらに詳述すると、本発明は、住居や住戸内に人が存在しているか存在していないかの推定に用いて好適な技術に関する。
以下の説明において、単位であることを明確にするために［ ］を用いる場合がある。
家庭における最適なエネルギーの使い方の探索や熱源機器の選択などを行うためには、電気やガスといったエネルギーの種類毎に用途毎の最終需要を把握することが必要である。しかしながら、ガスについては、配管切断を伴う流量計の設置は需要家の負担が大きく手間もコストもかかるため、個別機器の実測調査はほとんど行われていない。また、ガスメータの計測値から用途別（具体的には例えば、「給湯」，「厨房（コンロ）」，及び「暖房（ガスヒータ等）」の別）の消費量に分解する良好な手法は未だ無い。
ガスの配管を切断することなくガス流量を計測する装置として、例えば、ガス配管に設けられている膜式ガスメータに装着され、当該膜式ガスメータ内部の、ガスの流動に応じて周期運動する永久磁石の移動に伴う周期的な磁気変化を磁気センサによって検出し、当該周期的な磁気変化の検出に基づいて所定のガス流量（具体的には、ガスメータの単位計量体積）に対応したパルス出力を流量カウント数として所定時間毎に出力するものがある（例えば、特許文献１，特許文献２）。
上述のようなガス流量計測装置では、家庭におけるガスの消費量が所定時間毎にカウント数として把握され、所定時間は例えば１０秒などであり、１カウントに対応する実際のガス流量は例えば１．２［Ｌ］や１．０［Ｌ］などである。言い換えると、上述のようなガス流量計測装置では、所定の時間毎に、所定の流量毎に１カウントとして計数するカウント数として、ガスの消費量が出力される。
したがって、上述のようなガス流量計測装置では、具体的には例えば、所定時間が１０秒であると共に１カウント当たりの実際の流量が１．２［Ｌ］である場合には、１０秒間隔での、１．２Ｌ刻みのガス消費量が把握される。
また、ガス使用量を用途別に分解する従来の方法としては、例えば、世帯人数，厨房用ガス器具の所有台数，及び住居形態などのガス使用者の属性情報を用い、ガス使用者の世帯状況によってその使用量が異なる厨房用ガス使用量を求め、続いて、暖房を使用しない時期の総ガス使用量が厨房用と給湯用との和であることを利用してこの時期の給湯用ガス使用量を求め、次に、暖房を使用しない時期の給湯用ガス使用量と気温との相関関係を利用して暖房を使用する時期の給湯用ガス使用量を求め、さらに、総ガス使用量から厨房用及び給湯用ガス使用量を減算して暖房用ガス使用量を求めるものがある（特許文献３）。
また、電気に関しては、電気の使用（即ち、電力需要）から生活状態を推定することによって独居高齢者や家族を見守るサービスの実現が検討されたり、住居内の例えば空調等の機器の制御が自動で行われるシステムについて運用効率を改善して節電に繋げるために居住者が在宅であるか不在であるかの判定が行われたりしている例がある。具体的には例えばドアの開閉，人感センサ，或いはＣＯ2センサなどを利用するものがあり、また、住宅内への機器の設置が必要とされない非侵入型モニタリングの手法の一つとして電力需要データを用いるものが開発されている（非特許文献１）。
このような非侵入型モニタリングに関連しては、我が国において、スマートメータと呼ばれる、通信機能を備えて電力の利用状況（即ち、電力需要）をリアルタイムに確認することができる次世代電力計の低圧部門（具体的には、一般家庭）への設置が進められている（非特許文献２）。
特開２００４−１２５６３４号公報 特開２０１４−００６０５８号公報 特開２０００−２２１０６４号公報
Ｔ．Ａ．Ｎｇｕｙｅｎ ａｎｄ Ｍ．Ａｉｅｌｌｏ「Ｅｎｅｒｇｙ ｉｎｔｅｌｌｉｇｅｎｔ ｂｕｉｌｄｉｎｇｓ ｂａｓｅｄ ｏｎ ｕｓｅｒ ａｃｔｉｖｉｔｙ：Ａ ｓｕｒｖｅｙ」，Ｅｎｅｒｇｙ ａｎｄ Ｂｕｉｌｄｉｎｇｓ，Ｖｏｌ．５６，ｐｐ．２４４−２５７，２０１３年 経済産業省スマートメータ制度検討会「スマートメータの導入促進に伴う課題と対応について」，［online］，［平成２８年４月１４日アクセス］，インターネット＜ＵＲＬ：http://www.meti.go.jp/committee/summary/0004668/pdf/015_03_00.pdf＞
特許文献１や特許文献２のガス流量計測装置によって得られる計測データは、所定の時間毎に出力される、所定の流量毎に１カウントとして計数する流量カウント数である。しかしながら、例えばガスファンヒータやガスコンロのガス消費量は１．２［Ｌ／１０秒］に比べて通常は小さいため、１０秒間隔で出力される、実際の流量１．２［Ｌ］を１カウントとして計数する流量カウント数としての計測ではカウント値の変化は散発的（言い換えると、間欠的）になる。したがって、ガスファンヒータやガスコンロの各器具単体でのガス消費量に対応する計測値としての流量カウント数は、これら器具がたとえ継続的に使用されていても０若しくは１の出力となる。このため、時刻毎のガスの配管における実際のガス流量を適切に把握することは困難であり、延いては時刻毎の実際のガス流速を適切に把握することが困難であるという問題がある。
また、特許文献３のガス使用量の用途別分解方法では世帯人数，厨房用ガス器具の所有台数，及び住居形態などのガス使用者の属性情報が必要とされるところ、このような属性情報の入手は困難であり、仮に可能であったとしても多大な手間がかかってしまう。このため、汎用性が高いとは言い難い。加えて、特許文献３のガス使用量の用途別分解方法では気温の情報が更に必要とされるところ、対象世帯が所在する地域毎の、用途分解を行う期間の日毎の気温データの収集は多大な手間がかかるためにこの点においても汎用性が高いとは言い難く、さらに、気温データの地域区分の細かさによっては対象世帯の所在地における実際の気温との乖離が大きくなってしまうことも考えられるために良好な精度が常に確保されるとは言えず信頼性が高いとは言い難い。
そこで、本発明は、関連する情報やデータが必要とされることなく、対象住宅や対象住戸のガス総消費量の計測データであって値の変化の生起について所定の間隔で間欠的に把握され且つ値の変化について所定の量を単位として離散的に把握された計測データ群に基づいて家庭におけるガス消費量を用途別に良好な精度で分解することができるガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムを提供することを目的とする。
また、次世代電力計として設置が進められているスマートメータで自動検針された電力需要値を電力会社（具体的には、送配電部門）へと送るルートは「Ａルート」と呼ばれる。このＡルートから得られる電力需要データ（「Ａルート計測値」と呼ぶ）は、３０分間隔で出力され、且つ、１００ Ｗｈ 刻み（言い換えると、１００ Ｗｈ 単位）の数値である。したがって、Ａルート計測値では、３０分間における実際の電力需要として１００ Ｗｈ に満たない電力需要は、当該時刻（３０分時間帯）の計測値としては切り捨てられ（言い換えると、計上されず）、次の時刻（３０分時間帯）の計測値に繰り越される。
実際の電力需要に対応して出力されるＡルート計測値は、具体的には例えば以下の表１のようになる。なお、表１における「実際の電力需要」は、電力需要の瞬時値が仮に計測されているとした場合の仮想の値であり、３０分間積分値である。
表１に示す例から、３０分間の電力需要が１００ Ｗｈ に満たない時刻区分（３０分時間帯）では実際には電力需要が有るにも拘わらず、Ａルート計測値では電力需要が無いことになったり、或いは、連続する時刻区分（具体的には、１１：００と１１：３０）において実際の電力需要としての差違は小さい（具体的には、約１６ Ｗｈ）にも拘わらず、Ａルート計測値では１００ Ｗｈ もの差違になったりすることが確認される。
このため、Ａルート計測値をそのまま利用して居住者が在宅であるか不在であるかを判定した場合には、適切な判定が行われない虞があり、高い精度で判定を行うことは困難である。特に、３０分間における実際の電力需要が１００ Ｗｈ に満たない時間帯においては、Ａルート計測値としては０ Ｗｈ と１００ Ｗｈ とが交互に出力される状態になり、実際の電力需要に比してその差違が相対的に大きな割合となって誤判定が多くなる。この点に関連しては、省エネルギー機器への切り替えなどにより、エアコンや電子レンジなどの電力消費量が比較的大きい機器を利用しない時間帯における電力需要は減少する傾向にあり、前記誤判定は益々多くなると考えられる。
そこで、本発明は、追加的な特別の機器の設置及び関連する情報やデータが必要とされることなく、例えばスマートメータのＡルート計測値のような、対象住宅や対象住戸の電力需要量の計測データであって値の変化の生起について所定の間隔で間欠的に把握され且つ値の変化の態様について所定の量を単位として離散的に把握された計測データ群に基づいて住居や住戸内に人が存在しているか存在していないかを良好な精度で推定することができる居住者の在宅不在の判定方法、判定装置、及び判定プログラムを提供することを目的とする。
また、短い時間や特定の時刻毎に各々独立して（言い換えると、或る短い時間や特定の時刻のみを取り出して当該の或る短い時間や特定の時刻よりも前の状況を合わせて考慮すること無く）在宅不在を判定すると、居住者が在宅していても、短時間、使用する電気機器が少なかったり電気機器の使用態様に変化が無かったりする状況に対し、不在であると誤った判断をしてしまうことがあるという問題がある。
そこで、本発明は、短時間のみの状況に影響されること無く住居や住戸内に人が存在しているか存在していないかを推定することができる不在期間の確定方法を提供することを目的とする。
本発明の説明ではガス消費量を用途別に分解する際の単位となる家庭（世帯）の住宅や住戸へ差し向けられるガス配管のことを「主幹」と呼ぶと共に当該主幹を流れるガスの流量を「主幹のガス流量」又は「ガス総消費量」と呼び、本発明は当該ガス総消費量を用途別に分解するものとする。なお、住宅内や住戸内において、主幹から枝分かれした屋内配管が給湯用，厨房用，及び暖房用の各ガス器具へと接続される。
そして、本発明のガス消費量の用途分解方法は、ガス配管の主幹における計測によって所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量毎に１カウントとして計数する流量カウント数との組み合わせデータとして取得されたガス総消費量の計測データから、主幹から枝分かれした屋内配管に接続されているガス器具が使用されている時間帯の組み合わせデータが連続動作区間のデータとして抽出され、連続動作区間における計測時刻毎の流量カウント数の累積値と前記所定の流量とを掛け合わせた値を下限とすると共に前記掛け合わせた値に前記所定の流量を加えた値を上限とする領域内に存在し且つ連続動作区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測されたガス総消費量の総和と推定される各計測時刻のガス総消費量の総和とが一致する折れ線が推定され、当該折れ線における各計測時刻間の直線の傾きとして主幹におけるガスの流速が算定された上で、連続動作区間のデータが冬季であるか否か判断される処理と、ガスの流速が所定の範囲である状態が所定の時間に亙って継続しているか否かに基づいて暖房連続動作区間が特定される処理と、連続動作区間のデータにおいて特定された全ての暖房連続動作区間の中でのガスの流速の最小値が暖房ベース流速として設定される処理と、暖房ベース流速と給湯最小流速との和とガスの流速との大きさが比較される処理と、給湯最小流速とガスの流速との大きさが比較される処理と、厨房最小流速と間欠暖房最大流速とのうちの少なくとも一方とガスの流速との大きさが比較される処理とが実行されることにより、ガス総消費量の用途が給湯分，暖房分，及び厨房分に分類されるようにしている。
また、本発明のガス消費量の用途分解装置は、ガス配管の主幹における計測によって所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量毎に１カウントとして計数する流量カウント数との組み合わせデータとして取得されたガス総消費量の計測データから、主幹から枝分かれした屋内配管に接続されているガス器具が使用されている時間帯の組み合わせデータを連続動作区間のデータとして抽出する手段と、連続動作区間における計測時刻毎の流量カウント数の累積値と前記所定の流量とを掛け合わせた値を下限とすると共に前記掛け合わせた値に前記所定の流量を加えた値を上限とする領域内に存在し且つ連続動作区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測されたガス総消費量の総和と推定される各計測時刻のガス総消費量の総和とが一致する折れ線を推定して当該折れ線における各計測時刻間の直線の傾きとして主幹におけるガスの流速を算定する手段と、連続動作区間のデータが冬季であるか否かを判断し、ガスの流速が所定の範囲である状態が所定の時間に亙って継続しているか否かに基づいて暖房連続動作区間を特定すると共に連続動作区間のデータにおいて特定された全ての暖房連続動作区間の中でのガスの流速の最小値を暖房ベース流速として設定し、暖房ベース流速と給湯最小流速との和とガスの流速との大きさを比較し、給湯最小流速とガスの流速との大きさを比較し、さらに、厨房最小流速と間欠暖房最大流速とのうちの少なくとも一方とガスの流速との大きさを比較することにより、ガス総消費量の用途を給湯分，暖房分，及び厨房分に分類する手段とを有するようにしている。
また、本発明のガス消費量の用途分解プログラムは、ガス配管の主幹における計測によって所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量毎に１カウントとして計数する流量カウント数との組み合わせデータとして取得されたガス総消費量の計測データから、主幹から枝分かれした屋内配管に接続されているガス器具が使用されている時間帯の組み合わせデータを連続動作区間のデータとして抽出する処理と、連続動作区間における計測時刻毎の流量カウント数の累積値と前記所定の流量とを掛け合わせた値を下限とすると共に前記掛け合わせた値に前記所定の流量を加えた値を上限とする領域内に存在し且つ連続動作区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測されたガス総消費量の総和と推定される各計測時刻のガス総消費量の総和とが一致する折れ線を推定して当該折れ線における各計測時刻間の直線の傾きとして主幹におけるガスの流速を算定する処理と、連続動作区間のデータが冬季であるか否かを判断し、ガスの流速が所定の範囲である状態が所定の時間に亙って継続しているか否かに基づいて暖房連続動作区間を特定すると共に連続動作区間のデータにおいて特定された全ての暖房連続動作区間の中でのガスの流速の最小値を暖房ベース流速として設定し、暖房ベース流速と給湯最小流速との和とガスの流速との大きさを比較し、給湯最小流速とガスの流速との大きさを比較し、さらに、厨房最小流速と間欠暖房最大流速とのうちの少なくとも一方とガスの流速との大きさを比較することにより、ガス総消費量の用途を給湯分，暖房分，及び厨房分に分類する処理とをコンピュータに行わせるようにしている。
したがって、これらのガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムによると、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量毎に１カウントとして計数する流量カウント数との組み合わせデータ、すなわち、所定の時間間隔且つ所定の流量刻みという値の変化が散発的（言い換えると、間欠的）であり且つ離散的であるデータ、或いは言い換えると、値の変化の生起が時系列において散発的（間欠的）であり且つ値の変化の態様が離散的であるデータに基づいて、計測時刻毎のガスの流速が算定される。
これらのガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムによると、さらに、算定されたガスの流速に基づいてガスの総消費量の用途が給湯分，暖房分，及び厨房分に分類される。
また、本発明のガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムは、ガスの流速が算定された後に、連続動作区間のデータが冬季であるか否か判断される処理がまず行われ、暖房連続動作区間が特定される処理が次に行われ、暖房ベース流速が設定される処理が続いて行われるようにしても良い。この場合には、暖房用ガス器具の使用の有無が最初に判断されると共に暖房用ガス器具が使用されている場合にはガス総消費量から暖房分のガス消費量がまず切り分けられ、その後に給湯分と厨房分とのガス消費量の分類が行われる。
また、本発明の居住者の在宅不在の判定方法は、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量を計測単位として計量する電力需要量との組み合わせデータとして取得された電力需要量の計測データから、所定の時間帯の組み合わせデータが分析対象区間のデータとして抽出され、分析対象区間における計測時刻毎の電力需要量の累積値を下限とすると共に計測時刻毎の電力需要量の累積値に前記所定の電力量を加えた値を上限とする領域内に存在し且つ分析対象区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測された電力需要量の総和と推定される各計測時刻の電力需要量の総和とが一致する折れ線が推定され、当該折れ線における各計測時刻間の直線の傾きとして電力需要速度が算定され、在宅不在の判定の対象である計測時刻以前の所定の期間に含まれる電力需要速度が用いられて前記所定の期間における電力需要速度の統計量が特徴量として算出され、少なくとも前記特徴量が用いられて居住者の在宅不在が判定されるようにしている。
また、本発明の居住者の在宅不在の判定装置は、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量を計測単位として計量する電力需要量との組み合わせデータとして取得された電力需要量の計測データから、所定の時間帯の組み合わせデータを分析対象区間のデータとして抽出し、分析対象区間における計測時刻毎の電力需要量の累積値を下限とすると共に計測時刻毎の電力需要量の累積値に前記所定の電力量を加えた値を上限とする領域内に存在し且つ分析対象区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測された電力需要量の総和と推定される各計測時刻の電力需要量の総和とが一致する折れ線を推定して当該折れ線における各計測時刻間の直線の傾きとして電力需要速度を算定する手段と、在宅不在の判定の対象である計測時刻以前の所定の期間に含まれる電力需要速度を用いて前記所定の期間における電力需要速度の統計量を特徴量として算出すると共に少なくとも前記特徴量を用いて居住者の在宅不在を判定する手段とを有するようにしている。
また、本発明の居住者の在宅不在の判定プログラムは、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量を計測単位として計量する電力需要量との組み合わせデータとして取得された電力需要量の計測データから、所定の時間帯の組み合わせデータを分析対象区間のデータとして抽出し、分析対象区間における計測時刻毎の電力需要量の累積値を下限とすると共に計測時刻毎の電力需要量の累積値に前記所定の電力量を加えた値を上限とする領域内に存在し且つ分析対象区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測された電力需要量の総和と推定される各計測時刻の電力需要量の総和とが一致する折れ線を推定して当該折れ線における各計測時刻間の直線の傾きとして電力需要速度を算定する処理と、在宅不在の判定の対象である計測時刻以前の所定の期間に含まれる電力需要速度を用いて前記所定の期間における電力需要速度の統計量を特徴量として算出すると共に少なくとも前記特徴量を用いて居住者の在宅不在を判定する処理とをコンピュータに行わせるようにしている。
したがって、これらの居住者の在宅不在の判定方法、判定装置、及び判定プログラムによると、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量を計測単位として計量する電力需要量との組み合わせデータ、すなわち、所定の時間間隔且つ所定の電力量刻みという値の変化が散発的（言い換えると、間欠的）であり且つ離散的であるデータ、或いは言い換えると、値の変化の生起が時系列において散発的（間欠的）であり且つ値の変化の態様が離散的であるデータに基づいて、計測時刻毎の電力需要速度が算定される。
これらの居住者の在宅不在の判定方法、判定装置、及び判定プログラムによると、さらに、算定された電力需要速度に基づいて住居や住戸に人が存在しているか存在していないかが推定される。
また、本発明の居住者の在宅不在の判定方法、判定装置、及び判定プログラムは、電力需要量の計測データがスマートメータのＡルートから出力されるデータであるようにしても良い。この場合には、スマートメータのＡルート計測値が、電力需要量の計測データとして用いられる。
また、本発明の不在期間の確定方法は、電力需要量の計測データに基づいて時刻間毎の電力需要速度が算定され、判定対象時刻以前の所定の期間に含まれる電力需要速度が用いられて前記所定の期間における電力需要速度の統計量が特徴量として算出され、少なくとも前記特徴量が用いられて判定対象時刻別の居住者の在宅不在が判定され、所定の時間長さを有して時間軸方向に移動する判定期間に対応する判定対象時刻別の判定結果が全て不在であるときに当該判定期間中は不在であると判定され、不在と判定された一つの判定期間の始まりの時刻を不在開始時刻とすると共に終わりの時刻を不在終了時刻として特定する、または、連続して不在と判定された複数の判定期間のうちの最初の判定期間の始まりの時刻を不在開始時刻とすると共に最後の判定期間の終わりの時刻を不在終了時刻として特定するようにしている。
したがって、この不在期間の確定方法によると、判定対象時刻毎に各々独立して（即ち、当該の判定対象時刻よりも前の時刻における状況を合わせて考慮すること無く）在宅不在を判定する場合のような、居住者が在宅しているにもかかわらず、短時間、使用する電気機器が少なかったり電気機器の使用態様に変化が無かったりする状況での（言い換えると、電力需要速度の変化が小さい状況での）不在であるとの誤った判断をしてしまうことが防止される。
本発明のガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムによれば、値の変化の生起が時系列において散発的（間欠的）であり且つ値の変化の態様が離散的であるデータに基づいて計測時刻毎のガスの流速を算定することができるので、ガスの総消費量の実態の再現性の向上を図ることが可能になる。
本発明のガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムによれば、さらに、実態が良好に再現されたガスの流速に基づいてガスの総消費量の用途を給湯分，暖房分，及び厨房分に分類することができるので、ガスの総消費量の用途分解の良好な精度を確保してガス消費量の用途分解技術としての信頼性の向上を図ることが可能になる。
また、本発明のガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムは、連続動作区間のデータが冬季であるか否かがまず判断されてから暖房連続動作区間が特定された後に暖房ベース流速が設定されようにした場合には、暖房用ガス器具が使用されている場合にガス総消費量から暖房分のガス消費量がまず切り分けられてから給湯分と厨房分とのガス消費量の分類が行われるので、暖房分のガス消費量の分類を適確に行うことが可能になり、延いては、残りの給湯分と厨房分とのガス消費量の分類も適確に行うことが可能になる。
また、本発明の居住者の在宅不在の判定方法、判定装置、及び判定プログラムによれば、値の変化の生起が時系列において散発的（間欠的）であり且つ値の変化の態様が離散的であるデータに基づいて計測時刻毎の電力需要速度を算定することができるので、電力需要の実態の再現性の向上を図ることが可能になる。
本発明の居住者の在宅不在の判定方法、判定装置、及び判定プログラムによれば、さらに、実態が良好に再現された電力需要速度に基づいて住居や住戸に人が存在しているか存在していないかを推定することができるので、居住者の在宅不在の判定の良好な精度を確保して居住者の在宅不在の判定技術としての信頼性の向上を図ることが可能になる。
また、本発明の不在期間の確定方法によれば、短時間、使用する電気機器が少なかったり電気機器の使用態様に変化が無かったりする状況での（言い換えると、電力需要速度の変化が小さい状況での）不在であるとの誤った判断をしてしまうことを防止することができ、不在期間の確定手法としての信頼性の向上を図ることが可能になる。
本発明のガス消費量の用途分解方法の実施形態の一例を示すフローチャートである。 実施形態のガス消費量の用途分解方法をガス消費量の用途分解プログラムを用いて実施する場合の当該プログラムによって実現されるガス消費量の用途分解装置の機能ブロック図である。 実施形態の連続データの作成方法を説明する図であり、計測データとしての流量カウント数に基づいて流速を推定する方法を説明する図である。 ガス消費量の用途分解において用いられるアルゴリズムの一例を説明する用途決定チャートである。 本発明の居住者の在宅不在の判定方法の実施形態の一例を示すフローチャート（Ｓ１乃至Ｓ３）である。 本発明の居住者の在宅不在の判定方法の実施形態の一例を示すフローチャート（Ｓ４乃至Ｓ１１）である。 実施形態の居住者の在宅不在の判定方法を居住者の在宅不在の判定プログラムを用いて実施する場合の当該プログラムによって実現される居住者の在宅不在の判定装置の機能ブロック図である。 実施形態の連続データの作成方法を説明する図であり、計測データとしての電力需要量に基づいて電力需要速度を推定する方法を説明する図である。 真の累積需要曲線としての最短折れ線を推定する方法を説明する図である。 不在開始時刻及び不在終了時刻を特定する方法並びに不在期間を確定する方法の一例を説明する図である。
以下、本発明の構成を図面に示す実施の形態の一例に基づいて詳細に説明する。
図１から図４に、本発明のガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムの実施形態の一例を示す。
本実施形態のガス消費量の用途分解方法は、ガス配管の主幹における計測（Ｓ１）によって所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量ｆu毎に１カウントとして計数する流量カウント数との組み合わせデータとして取得（Ｓ１）されたガス総消費量の計測データから、主幹から枝分かれした屋内配管に接続されているガス器具が使用されている時間帯の組み合わせデータが連続動作区間のデータとして抽出され（Ｓ２）、連続動作区間における計測時刻毎の流量カウント数の累積値と流量ｆuとを掛け合わせた値を下限とすると共に前記掛け合わせた値に流量ｆuを加えた値を上限とする領域内に存在し且つ連続動作区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測されたガス総消費量の総和と推定される各計測時刻のガス総消費量の総和とが一致する折れ線が推定され、当該折れ線における各計測時刻間の直線の傾きとして主幹におけるガスの流速ｙtが算定された上で（Ｓ３）、連続動作区間のデータが冬季であるか否か判断される処理（Ｓ４−１）と、ガスの流速ｙtが所定の範囲である状態が所定の時間に亙って継続しているか否かに基づいて暖房連続動作区間が特定される処理（Ｓ４−２）と、連続動作区間のデータにおいて特定された全ての暖房連続動作区間の中でのガスの流速ｙtの最小値が暖房ベース流速として設定される処理（Ｓ４−３）と、暖房ベース流速と給湯最小流速との和とガスの流速ｙtとの大きさが比較される処理（Ｓ４−４）と、給湯最小流速とガスの流速ｙtとの大きさが比較される処理（Ｓ４−５）と、連続動作区間のデータが夏季であるか否か判断される処理（Ｓ４−６）と、厨房最小流速と間欠暖房最大流速とのうちの少なくとも一方とガスの流速ｙtとの大きさが比較される処理（Ｓ４−７）とが実行されることにより、ガス総消費量の用途が給湯分，暖房分，及び厨房分に分類される（Ｓ４）ようにしている（図１参照）。
また、本実施形態では、Ｓ３の処理において、連続動作区間における計測時刻の時系列順に付与される識別子ｔ（ただし、ｔ＝１，２，３，…，Ｔ）毎の当該計測時刻の識別子ｔに対応する計測時刻における流量カウント数と流量ｆuとを掛け合わせた値ｘtが算出され、計測時刻の識別子ｔ毎の値ｘtが用いられて数式１乃至数式３に示す最適化問題のうちのいずれかが解かれることによって計測時刻の識別子ｔ毎の値ｙtとして主幹におけるガスの流速ｙtが算定されるようにしている。
上記ガス消費量の用途分解方法は、ガス消費量の用途分解装置によって実施され得る。本実施形態のガス消費量の用途分解装置１０は、ガス配管の主幹における計測によって所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量ｆu毎に１カウントとして計数する流量カウント数との組み合わせデータとして取得されたガス総消費量の計測データの入力を受ける手段としてのデータ受部１１ａと、計測データから、主幹から枝分かれした屋内配管に接続されているガス器具が使用されている時間帯の組み合わせデータを連続動作区間のデータとして抽出する手段としての分解区間抽出部１１ｂと、連続動作区間における計測時刻毎の流量カウント数の累積値と流量ｆuとを掛け合わせた値を下限とすると共に前記掛け合わせた値に流量ｆuを加えた値を上限とする領域内に存在し且つ連続動作区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測されたガス総消費量の総和と推定される各計測時刻のガス総消費量の総和とが一致する折れ線を推定して当該折れ線における各計測時刻間の直線の傾きとして主幹におけるガスの流速ｙtを算定する手段としての最短折線作成部１１ｃと、連続動作区間のデータが冬季であるか否かを判断し、ガスの流速ｙtが所定の範囲である状態が所定の時間に亙って継続しているか否かに基づいて暖房連続動作区間を特定すると共に連続動作区間のデータにおいて特定された全ての暖房連続動作区間の中でのガスの流速ｙtの最小値を暖房ベース流速として設定し、暖房ベース流速と給湯最小流速との和とガスの流速ｙtとの大きさを比較し、給湯最小流速とガスの流速ｙtとの大きさを比較し、連続動作区間のデータが夏季であるか否かを判断し、さらに、厨房最小流速と間欠暖房最大流速とのうちの少なくとも一方とガスの流速ｙtとの大きさを比較することにより、ガス総消費量の用途を給湯分，暖房分，及び厨房分に分類する手段としての用途分解部１１ｄとを有する。
また、本実施形態では、最短折線作成部１１ｃが、連続動作区間における計測時刻の時系列順に付与される識別子ｔ（ただし、ｔ＝１，２，３，…，Ｔ）毎の当該計測時刻の識別子ｔに対応する計測時刻における流量カウント数と流量ｆuとを掛け合わせた値ｘtを算出すると共に、計測時刻の識別子ｔ毎の値ｘtを用いて数式１乃至数式３に示す最適化問題のうちのいずれかを解くことによって計測時刻の識別子ｔ毎の値ｙtとして主幹におけるガスの流速ｙtを算定するようにしている。
また、上記ガス消費量の用途分解方法及びガス消費量の用途分解装置は、ガス消費量の用途分解プログラムがコンピュータ上で実行されることによっても実施・実現され得る。ここでは、ガス消費量の用途分解プログラムがコンピュータ上で実行されることによってガス消費量の用途分解装置が実現されると共にガス消費量の用途分解方法が実施される場合を説明する。
本実施形態のガス消費量の用途分解プログラム１７を実行するためのコンピュータ１０（本実施形態では、ガス消費量の用途分解装置１０でもある）の全体構成を図２に示す。このコンピュータ１０（ガス消費量の用途分解装置１０）は制御部１１，記憶部１２，入力部１３，表示部１４，及びメモリ１５を備え、これらが相互にバス等の信号回線によって接続されている。
制御部１１は、記憶部１２に記憶されているガス消費量の用途分解プログラム１７によってコンピュータ１０全体の制御並びにガス総消費量の用途の分解に係る演算を行うものであり、例えばＣＰＵ（中央演算処理装置）である。
記憶部１２は、少なくともデータやプログラムを記憶可能な装置であり、例えばハードディスクである。
入力部１３は、少なくとも作業者の命令や種々の情報を制御部１１に与えるためのインターフェイス（即ち、情報入力の仕組み）であり、例えばキーボードやマウスである。なお、例えばキーボードとマウスとの両方のように複数種類のインターフェイスを入力部１３として有するようにしても良い。
表示部１４は、制御部１１の制御によって文字や図形或いは画像等の描画・表示を行うものであり、例えばディスプレイである。
メモリ１５は、制御部１１が種々の制御や演算を実行する際の作業領域であるメモリ空間となるものであり、例えばＲＡＭ(Ｒandom Ａccess Ｍemory の略)である。
そして、コンピュータ１０（「ガス消費量の用途分解装置１０」と表記する）の制御部１１には、ガス消費量の用途分解プログラム１７が実行されることにより、ガス配管の主幹における計測によって所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量ｆu毎に１カウントとして計数する流量カウント数との組み合わせデータとして取得されたガス総消費量の計測データの入力を受ける処理を行うデータ受部１１ａと、計測データから、主幹から枝分かれした屋内配管に接続されているガス器具が使用されている時間帯の組み合わせデータを連続動作区間のデータとして抽出する処理を行う分解区間抽出部１１ｂと、連続動作区間における計測時刻毎の流量カウント数の累積値と流量ｆuとを掛け合わせた値を下限とすると共に前記掛け合わせた値に流量ｆuを加えた値を上限とする領域内に存在し且つ連続動作区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測されたガス総消費量の総和と推定される各計測時刻のガス総消費量の総和とが一致する折れ線を推定して当該折れ線における各計測時刻間の直線の傾きとして主幹におけるガスの流速ｙtを算定する処理を行う最短折線作成部１１ｃと、連続動作区間のデータが冬季であるか否かを判断し、ガスの流速ｙtが所定の範囲である状態が所定の時間に亙って継続しているか否かに基づいて暖房連続動作区間を特定すると共に連続動作区間のデータにおいて特定された全ての暖房連続動作区間の中でのガスの流速ｙtの最小値を暖房ベース流速として設定し、暖房ベース流速と給湯最小流速との和とガスの流速ｙtとの大きさを比較し、給湯最小流速とガスの流速ｙtとの大きさを比較し、連続動作区間のデータが夏季であるか否かを判断し、さらに、厨房最小流速と間欠暖房最大流速とのうちの少なくとも一方とガスの流速ｙtとの大きさを比較することにより、ガス総消費量の用途を給湯分，暖房分，及び厨房分に分類する処理を行う用途分解部１１ｄとが構成される。
また、本実施形態では、最短折線作成部１１ｃが、連続動作区間における計測時刻の時系列順に付与される識別子ｔ（ただし、ｔ＝１，２，３，…，Ｔ）毎の当該計測時刻の識別子ｔに対応する計測時刻における流量カウント数と流量ｆuとを掛け合わせた値ｘtを算出すると共に、計測時刻の識別子ｔ毎の値ｘtを用いて数式１乃至数式３に示す最適化問題のうちのいずれかを解くことによって計測時刻の識別子ｔ毎の値ｙtとして主幹におけるガスの流速ｙtを算定する処理を行うようにしている。
そして、ガス消費量の用途分解方法の実施として、まず、ガス総消費量の計測が行われて計測データの取得が行われる（Ｓ１）。
ガス総消費量の計測の仕方は、所定の流量毎に１カウントとして計数するカウント数として主幹のガス流量を計測し、所定の時間毎に例えば電気的な信号として前記カウント数を出力するものであれば、特定の方法や機器によるものに限定されるものではなく、前記のような計測を行い得る方法や機器が適宜選択される。
具体的には例えば、あくまで一例として挙げると、膜式ガスメータに装着される前述のガス流量計測装置（例えば、特開２００４−１２５６３４号公報，特開２０１４−６０５８号公報）が利用され得る。あるいは、一定量（例えば、ガスメータの単位計量体積）のガスがガスメータを通過するたびに動作サイクルに起因してガスメータから発生する周期的な動作音や振動を検出してその回数をカウント数として所定の時間毎に出力する機器が利用され得る。
ガス総消費量（主幹のガス流量）の計測における、上述の「所定の流量」としての、１カウント当たりの実際の流量は、例えば、主幹に設けられているガスメータの単位計量体積に従って決定される。具体的には例えば、ガスメータの単位計量体積が１［Ｌ］である場合にはガス流量の計測機器から出力される１カウント毎に１［Ｌ］のガスが実際に流動した（言い換えると、消費された）ことを意味し、或いは、ガスメータの単位計量体積が１．２［Ｌ］である場合にはガス流量の計測機器から出力される１カウント毎に１．２［Ｌ］のガスが実際に流動したことを意味する。
また、ガス総消費量の計測における、上述の「所定の時間」としての、カウント数を出力する時間間隔は、特定の時間［秒，分］に限定されるものではなく、計測に用いられる計測機器の仕様を踏まえた上で、例えばガス消費量に関して想定される変動ピッチや必要とされる推定精度などが考慮されて適当な時間に適宜設定される。具体的には例えば、あくまで一例として挙げると、５〜６０秒程度の範囲で適当な値に設定され得る。
本実施形態では、上述の「所定の流量」が１．２［Ｌ］であると共に上述の「所定の時間」が１０秒である計測データ、すなわち、実際の流量１．２［Ｌ］を１カウントとして計数するカウント数（「流量カウント数」と呼ぶ）が１０秒間隔で出力される計測データがＳ１の処理によって取得されるとする。
そして、計測によって取得されたガス総消費量（主幹のガス流量）の計測値（具体的には、流量カウント数）は、データ受部１１ａを介してガス消費量の用途分解装置１０に入力される。
ガス総消費量の計測値（流量カウント数）は、ガス消費量の用途分解装置１０に、記録媒体やデータサーバを介して入力されるようにしても良いし、ガス消費量の用途分解装置１０と計測機器とがデータや制御指令等の信号の送受信（即ち、出入力）が可能であるように電気的に接続されて入力されるようにしても良い。
記録媒体を介して計測値が入力される場合には、計測機器から出力された流量カウント数が記録媒体に記録され、当該記録媒体がガス消費量の用途分解装置１０の記録媒体接続用端子（図示していない）に差し込まれることによって流量カウント数が入力されるようにしても良い。
また、データサーバを介して計測値が入力される場合には、データサーバがバス等の信号回線によってガス消費量の用途分解装置１０に接続され、計測機器から出力された流量カウント数が前記データサーバにデータファイル等として格納(保存)され、当該データファイル等として保存された流量カウント数が読み込まれるようにしても良い。
また、ガス消費量の用途分解装置１０と計測機器とが電気的に接続されて計測値が入力される場合には、例えば、各々に接続されて敷設されたケーブル等が用いられる有線による信号送受の仕組みを介して信号の送受信が可能であるように電気的に接続されるようにしても良いし、各々に接続された無線信号送受信機が用いられる無線による信号送受の仕組みを介して信号の送受信が可能であるように電気的に接続されるようにしても良い。そして、これら信号送受の仕組みによって流量カウント数が入力されるようにしても良い。
さらに、記録媒体やデータサーバや有線・無線による信号送受の仕組みが組み合わされて用いられて流量カウント数がガス消費量の用途分解装置１０に入力されるようにしても良い。
また、ガス総消費量の各計測値（各流量カウント数）が計測機器から出力された時刻が計測機器に備えられた時計機能によって計測値と一緒に出力されて当該計測値と対応づけられて記録されたり、或いは、計測機器から出力された各計測値が記録媒体やデータサーバやガス消費量の用途分解装置１０に記録されたり入力されたりした時の時刻が当該計測値と対応づけられて記録されたりする。そして、この時刻の間隔［秒，分］が上述の「所定の時間」になる。
ここでの説明においては、各計測値（各流量カウント数）と対応づけられて記録される時刻のことを「計測時刻」と呼び、流量カウント数と計測時刻との組み合わせデータのことを「計測データ」と呼ぶ。
そして、データ受部１１ａにより、計測データがメモリ１５に記憶させられる。
なお、例えば、対象の家庭に関する（言い換えると、対象住宅や対象住戸の主幹における）、あくまで一例として挙げると２４時間に亙る一年分の計測データが、ガス消費量の用途分解のひとかたまりの処理対象としてメモリ１５に記憶させられる。
次に、制御部１１の分解区間抽出部１１ｂにより、Ｓ１の処理によって計測されて取得された計測データからの、用途の分解が行われる時間帯のデータの抽出が行われる（Ｓ２）。
Ｓ２の処理は、Ｓ１の処理においてメモリ１５に記憶された、ガス消費量の用途分解のひとかたまりの処理対象としての計測データの全体に亙って実行される。
本発明はガス総消費量を用途別に分解するものであるので、対象住宅・住戸においてガス器具が使用されてガスが実際に消費されている時間帯が処理対象になり、ガス器具が使用されておらずガスが消費されていない時間帯は処理対象にならない。
そこで、計測データから、何らかのガス器具が使用される（言い換えると、動作する）ことによってガス総消費量が連続して変化（具体的には、増加）している可能性がある時間帯（「連続動作区間」と呼ぶ）の計測データが切り出される。
具体的には、動作休止判定時間［秒，分］に亙って計測データの流量カウント数が０(零)である場合に、ガス器具の使用が途切れたと判定される。動作休止判定時間以上に亙って計測データの流量カウント数が０である時間帯のことを「動作休止区間」と呼ぶ。
そして、時間的に前後する二つの動作休止区間に挟まれる時間帯が連続動作区間として切り出される。したがって、連続動作区間は、計測データの流量カウント数が動作休止判定時間以上に亙って０になることなく継続的に変化している時間帯の計測データであるとも言える。
動作休止判定時間［秒，分］は、特定の値に限られるものではなく、必要とされる推定精度などが考慮されて適当な値に適宜設定される。具体的には例えば、あくまで一例として挙げると、主幹におけるガスの流速が０．２［Ｌ／１０秒］以下である変化は微小な流量として無視しても大きな誤差は生じない（言い換えると、必要な精度は確保される）と判断される場合に、１．２Ｌ／０．２Ｌ＝６（即ち、計測時刻の間隔６個分）に基づいて、動作休止判定時間が１０秒×６＝６０秒に設定されることが考えられる。
なお、動作休止判定時間の値［秒，分］は、ガス消費量の用途分解プログラム１７内に予め規定されるようにしても良いし、Ｓ２の処理が行われる前に入力部１３を介して作業者によって入力されるようにしても良い。
本実施形態では、分解区間抽出部１１ｂにより、Ｓ１の処理においてメモリ１５に記憶された計測データが時系列に従って順次読み込まれ、動作休止判定時間に基づいて動作休止区間が特定されると共に、時間的に前後する二つの動作休止区間に挟まれる時間帯が連続動作区間として抽出される。
そして、分解区間抽出部１１ｂにより、抽出された連続動作区間毎の計測データが、以降のＳ３及びＳ４の処理対象の計測データとしてメモリ１５に記憶させられる。
次に、制御部１１の最短折線作成部１１ｃにより、Ｓ２の処理によって抽出された連続動作区間に関する計測データが用いられて累積需要曲線の作成が行われる（Ｓ３）。
Ｓ３の処理は、Ｓ２の処理において抽出された連続動作区間の各々について実行される。
Ｓ１の処理によって取得される計測データは、１０秒間隔での１．２Ｌ刻みのガス総消費量である。しかしながら、何らかのガス器具が使用されているときは、ガス総消費量が常に増加するという点においてガス総消費量の変化の生起は時系列において連続であり、また、ガスが常に流動するという点においてガス総消費量の変化の態様は漸次的である。
つまり、Ｓ１の処理における計測では、実際には(値の)変化の生起が時系列において連続であると共に(値の)変化の態様が漸次的である現象が、所定の時間間隔且つ所定の値刻みという散発的（言い換えると、間欠的）且つ離散的な変化として把握されていると言える。
そして、本発明では、所定の時間間隔且つ所定の流量刻みという計測値の変化が散発的（言い換えると、間欠的）であり且つ離散的である計測データがそのまま用いられるのではなく、当該散発的（間欠的）且つ離散的な計測データが、計測値の変化の生起が時系列において連続であり且つ計測値の変化の態様が漸次的であるデータに変換された上で用途分解の処理に用いられる。
本発明のガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムにおいて用いられる、計測値の変化の生起が時系列において散発的（間欠的）であり且つ計測値の変化の態様が離散的であるデータに基づいて、計測値の変化の生起が時系列において連続であり且つ計測値の変化の態様が漸次的である累積需要曲線が作成される、連続データの作成方法を、図３を用いて説明する。
なお、連続動作区間では、何らかのガス器具が使用されることによってガス総消費量が常に増加するという点においてガス総消費量の変化の生起が時系列において連続であり、また、何らかのガス器具が使用されることによってガスが常に流動するという点においてガス総消費量の変化の態様が漸次的である。
ここでの説明の例として、図３に示すような、或る連続動作区間に関する計測データ（具体的には、１０秒毎に出力される流量カウント数の累積値；図中の△印）がＳ２の処理によって抽出されたとする。なお、１カウント当たりの実際の流量は１．２［Ｌ］である。
図３の横軸は、１０秒間隔である計測時刻の各時点を、当該の連続動作区間の最初の計測時刻を１とした上で識別子としての連続番号で表したものである。すなわち、横軸における例えば「１０」は、当該の連続動作区間の始まりから１００秒（＝１０×１０秒）だけ時間が経過していることを表す。
また、図３の縦軸は、当該の連続動作区間おける流量カウント数の累積値を表す。すなわち、縦軸における例えば「４」は、当該の連続動作区間の始まりから４．８［Ｌ］（＝４×１．２［Ｌ］）だけガスが消費されていることを表す。
Ｓ１の処理によって取得される計測データは所定の時間間隔且つ所定の流量刻みという値の変化が散発的（間欠的）且つ離散的なデータであるので、すなわち、計測値の変化の生起が時系列において散発的（間欠的）且つ計測値の変化の態様が離散的なデータであるので、データプロットである△印を時系列で結ぶと図３に示すように階段状になる（図３中において符号２１；「下限階段直線２１」と呼ぶ）。
また、計測データの流量カウント数の累積値に計測時刻毎に１を加えたデータのプロットである▽印を時系列で結んだものも、図３に示すように階段状になる（図３中において符号２２；「上限階段直線２２」と呼ぶ）。
そして、下限階段直線２１が実際のガス総消費量に基づく真の累積値の下限（言い換えると、実際のガス総消費量が取り得る累積値の下限）であり、上限階段直線２２が実際のガス総消費量に基づく真の累積値の上限（言い換えると、実際のガス総消費量が取り得る累積値の上限）である。
したがって、実際のガス総消費量の時系列における累積を表す真の累積需要曲線は、下限階段直線２１と上限階段直線２２とに挟まれた領域（「存在可能領域」と呼ぶ）内に存在し、全体としては増加傾向を示す非減少の曲線であると考えられる。
また、家庭において一般的に使用されるガス器具の実際のガス消費の特性を本発明者が分析することにより、ガス器具によって消費されるガスの流量の大きな変化の回数は比較的少なく、ほぼ一定の流速が或る程度の時間に亙って継続する傾向があることが確認された。
このため、真の累積需要曲線としては、存在可能領域内において増加傾向を示す非減少の折れ線のうち、折れ点（言い換えると、変曲点、若しくは、流速変化点）の個数が比較的少ないものが妥当性が高いと考えられる。
そこで、実際のガス総消費量の時系列における累積を表す真の累積需要曲線として、存在可能領域内に存在する、総長（即ち、当該の連続動作区間の始まりから終わりまでに亙る軌跡の長さ）が最短の折れ線（「最短折れ線」と呼ぶ）を推定する。
真の累積需要曲線としての、存在可能領域内に存在する、最短折れ線の推定手法は、特定の方法に限定されるものではなく、特定の領域内に存在して所定の始点から終点へと至る折れ線のうちの総長が最短の折れ線であり且つ各計測時刻において計測された電力需要量の総和と推定される各計測時刻の電力需要量の総和とが一致する折れ線を推定し得る適当な方法が適宜選択される。
最短折れ線の推定手法として、具体的には例えば、あくまで一例として挙げると、当該の連続動作区間において確定済みの折れ線の最後の点を起点とした（未確定の）直線部の許容される傾きの上限値・下限値を順次更新し、全ての直線が存在可能領域を逸脱せざるを得なくなる時点で直線の傾きを決定し、その直線が存在可能領域の境界に最初に接した時点までの折れ線の確定を行い、そして、新たに確定した折れ線の最後の点を新たな起点に設定した上で上述と同様の処理を繰り返すことを基本的な考え方とする方法が用いられ得る。
この推定方法は、計測時刻ｔ（尚、計測単位は３０分であり、したがって計測時刻の間隔は３０分である）における計測流量をｘt［Ｌ／１０秒］とすると共に推定流速をｙt［Ｌ／１０秒］とすると、最適化問題として数式４のように表される。
数式４においてｔは連続動作区間における計測時刻の識別子であってｔ＝１，２，３，…であり、数式４では当該の連続動作区間の範囲として計測時刻の識別子ｔが１からＴまでの期間に亙って処理が実行される。
そして、計測時刻を１刻みの識別子とすることにより、計測の時間間隔である１０秒を一単位として扱うことになり、計測時刻 ｔ−１ から ｔ にかけての増分は計測の時間間隔当たりの増加量であって計測時刻 ｔ−１ から ｔ にかけての速度［Ｌ／１０秒］であり、真の累積需要曲線としては計測時刻 ｔ−１ から ｔ にかけての直線の傾きである。
数式４において、また、計測流量ｘt［Ｌ／１０秒］は、計測時刻ｔにおける流量カウント数（言い換えると、計測時刻 ｔ−１ から ｔ にかけての流量カウント数の増分）と１カウント当たりの実際の流量［Ｌ］とを掛け合わせたものである。
数式４において、さらに、ｆuは、１カウント当たりの実際の流量［Ｌ］である（即ち、本実施形態ではｆu＝１．２である）。なお、ｆuの値［Ｌ］は、ガス消費量の用途分解プログラム１７内に予め規定されるようにしても良いし、Ｓ３の処理が行われる前に入力部１３を介して作業者によって入力されるようにしても良い。
制約条件の第一式は流速の非負性に対応し、第二式は実際のガス総消費量の累積値の範囲制約に対応し、さらに、第三式は当該の連続動作区間における計測値ベースのガス総消費量の総和と推定値ベースのガス総消費量の総和との一致に対応する。
上述の最適化問題を解く方法によると、最悪の場合でも、時系列長（即ち、計測時刻の個数）の２乗のオーダーで真の累積需要曲線に相当する最短折れ線が推定される。
図３に示す計測データのプロット（図中の△印）を結んだ下限階段直線２１に対して数式４で表される最適化問題を適用した結果として、図３において点線で示す、実際のガス総消費量に近いと考えられる推定累積需要曲線２３が得られる。
ここで、上記の最適化問題は計測時刻の識別子ｔ＝１から処理が実行されるところ、ｔ＝１の処理が実行されるためには、ｔ＝１での起点としてｔ＝０（即ち、計測時刻の識別子ｔ＝１に対応する計測時刻の一つ前の計測時刻）からｔ＝１にかけての折れ線の最後の点が特定される必要が、即ち、ｔ＝０からｔ＝１にかけての直線が特定される必要がある。この直線は、具体的には、当該の連続動作区間の最初の計測時刻（即ち、計測時刻の識別子ｔ＝１）の直前の計測時刻として識別子ｔを０(零)とすると共に流量カウント数の累積値を０(零)とする仮想的な原点を置き、当該原点と計測時刻の識別子ｔ＝１における流量カウント数の累積値の値の点とを結ぶことによって特定される。
あるいは、ほぼ一定の流速が或る程度の時間に亙って継続する傾向があるという家庭におけるガス消費の特性に基づいて、真の累積需要曲線としては折れ点の個数が少ないものが妥当性が高いという考えのもと、当該の連続動作区間についての計測時刻の識別子ｔ＝０からｔ＝１にかけての直線の傾きが、識別子ｔ＝１において連続する直線の傾きに一致させられるようにしても良い。
ここで、真の累積需要曲線の上述の推定方法（即ち、数式４で表される最適化問題）は、目的関数部を「距離の和」から「距離の２乗和」に変更した数式５で表される最適化問題、即ち「存在可能領域の任意の曲線のうち各計測時刻の流速（即ち、流量の累積値の計測時刻間差分）の２乗和が最小の曲線」を求める問題と等価である。
数式５のように表される最適化問題では、２乗和の最小化により、各計測時刻の流速が同じ値を持つ傾向が生じるので、折れ線が最適解になる。
真の累積需要曲線の推定方法としての数式４で表される最適化問題は、また、数式５で表される最適化問題の双対問題である数式６で表される最適化問題とも等価である。なお、数式６で表される最適化問題は、全変動正規化問題と呼ばれる凸最適化問題の一種である。なお、数式６において、ｘ0は計測流量ｘt［Ｌ／１０秒］の初期値（即ち、ｔ＝０のときの値）であり、図３に示す例ではｘ0＝０である。
上述の真の累積需要曲線（即ち、最短折れ線）は、一例として以下に示すアルゴリズムが用いられることによって求められ得る。
−真の累積需要曲線の決定アルゴリズム−
《保持・更新情報》
○起点 折れ線
○上限直線 上限需要と領域境界接触点
○下限直線 下限需要と領域境界接触点
《更新内容》
(1) 時刻更新：時刻＝時刻＋１
(2) 下限直線＜上限累積需要値
起点から下限境界接触点までの需要量＝需要下限（確定）
起点＝下限直線の境界接触点
時刻＝起点，上限／下限直線の再設定
(1)へ
(3) 上限直線＞下限累積需要
起点から上限境界接触点までの需要量＝需要上限（確定）
起点＝上限直線の境界接触点
時刻＝起点，上限／下限直線の再設定
(1)へ
(4) 下限直線＜領域下限
下限需要，下限境界設定更新
上限直線＞領域上限
上限需要，上限境界設定更新
(5) 計測時刻の識別子ｔ＝Ｔになるまで(1)へ戻る
本実施形態では、最短折線作成部１１ｃにより、Ｓ２の処理においてメモリ１５に記憶された連続動作区間に関する計測データ（即ち、流量カウント数と計測時刻との組み合わせデータ）が読み込まれる。
そして、メモリ１５から読み込まれた流量カウント数は１カウント当たりの実際の流量である１．２［Ｌ］と掛け合わされて実際のガスの流量［Ｌ］に変換されると共に計測時刻は当該の連続動作区間における計測時刻の識別子として時系列順に１，２，３，…という連続番号に変換される。
続いて、最短折線作成部１１ｃにより、上記変換された計測データが用いられて数式４乃至数式６に表される最適化問題のうちのいずれかが解かれて、計測時刻の識別子ｔ毎のガスの推定流速ｙt［Ｌ／１０秒］が算定される。
そして、最短折線作成部１１ｃにより、Ｓ２の処理において抽出された連続動作区間毎に、算定されたガスの推定流速ｙt［Ｌ／１０秒］の値が計測時刻の識別子ｔの値と対応づけられて（言い換えると、計測時刻の識別子ｔの値とガスの推定流速ｙtの値との組み合わせデータとして）メモリ１５に記憶させられる。
次に、制御部１１の用途分解部１１ｄにより、Ｓ３の処理によって算定された計測時刻別のガスの流速が用いられてガス総消費量の用途毎のガス使用量への分解が行われる（Ｓ４）。
Ｓ４の処理は、Ｓ２の処理において抽出されてＳ３の処理において計測時刻別のガスの流速が算定された連続動作区間の各々について実行される。
具体的には、用途分解部１１ｄにより、Ｓ３の処理においてメモリ１５に記憶された連続動作区間に関する計測時刻の識別子ｔの値（ｔ＝１，２，…，Ｔ）とガスの推定流速ｙtの値との組み合わせデータが読み込まれる。
続いて、本実施形態では、図４に示すアルゴリズムに従ってガス総消費量が用途毎のガス使用量に分解される。
具体的には、まず、当該の連続動作区間が冬季におけるものであるか否かが判断される（Ｓ４−１）。
冬季であるか否かの判断の仕方は、特定の判断基準に限定されるものではなく、冬季とその他の季節とにおける例えば暖房器具の使用有無や給湯器具の温度設定などの差異が顕れる（若しくは、顕れると考えられる）時期などが考慮された上で適当な方法に適宜設定される。具体的には例えば、あくまで一例として挙げると、１０月から翌年４月までの期間に取得された計測データから抽出された連続動作区間のデータは冬季であると判断する方法が採用され得る。
そして、当該の連続動作区間が冬季におけるものである場合には（Ｓ４−１：Ｙｅｓ）Ｓ４−２の処理に進む。
そして、推定流速ｙtが暖房基準流速の最小値Ｖhs-min以上であると共に暖房基準流速の最大値Ｖhs-max以下であり、且つ、前記条件を満たす流速の継続時間（「所定流速継続時間」と呼ぶ）が暖房運転判断時間Ｔh以上であるか否かが判断される（Ｓ４−２）。
Ｓ４−２の判断処理は、流速は必ずしも大きくはないもののそのような流速の継続時間が長い場合には少なくとも暖房用ガス器具が使用されていると考え、暖房用ガス器具が連続して動作している時間帯（「暖房連続動作区間」と呼ぶ）を特定するための処理である。
暖房基準流速の最小値Ｖhs-min及び最大値Ｖhs-max［Ｌ／１０秒］は、暖房用ガス器具が使用されているときのガス消費量として想定される範囲の最小値及び最大値に基づいて設定され、特定の値に限定されるものではなく、例えば、対象住宅・住戸に設置されている（若しくは、設置されていると想定される一般的な）暖房用ガス器具の仕様などが考慮されたり、或いは、実測データが収集されている場合には当該実測データが考慮されたりして適当な値に適宜設定される。具体的には例えば、あくまで一例として挙げると、最小値Ｖhs-minが０．０１〜０．２［Ｌ／１０秒］程度の範囲で適当な値に、最大値Ｖhs-maxが０．５〜１．５［Ｌ／１０秒］程度の範囲で適当な値に、それぞれ設定され得る。
また、暖房運転判断時間Ｔh［秒，分］は、特定の値に限定されるものではなく、例えば暖房器具を使用する場合の最低継続時間として想定される時間などが考慮されて適当な値に適宜設定される。具体的には例えば、あくまで一例として挙げると、５〜１５分程度の範囲で適当な値に設定され得る。
Ｓ４−２の処理は、当該の連続動作区間の全体に亙って実行される。そして、Ｓ４−２の処理の結果として、各連続動作区間について、暖房連続動作区間が一つも特定されないか、または、暖房連続動作区間が一つ以上特定される。
そして、Ｖhs-min≦ｙt≦Ｖhs-max であると共に所定流速継続時間がＴh［分］以上である場合には（Ｓ４−２：Ｙｅｓ）Ｓ４−３の処理に進む。言い換えると、Ｓ４−２の処理において暖房連続動作区間が一つ以上特定された場合にＳ４−３の処理に進む。
そして、暖房ベース流速Ｖhbが設定される（Ｓ４−３）。
暖房ベース流速Ｖhb［Ｌ／１０秒］は、当該の連続動作区間において特定された全ての暖房連続動作区間の中での、推定流速ｙt［Ｌ／１０秒］の最小値に設定される。
したがって、暖房ベース流速Ｖhbは、或る一つの連続動作区間に対して共通に用いられる値として設定され、異なる連続動作区間では値が異なることがあり得る。
次に、推定流速ｙtが暖房ベース流速Ｖhbと給湯最小流速Ｖwとの和以上であるか否かが判断される（Ｓ４−４）。
Ｓ４−４の処理は、Ｓ４−２の処理において特定された暖房連続動作区間それぞれの各計測時刻ｔについて実行される。
Ｓ４−４の判断処理は、当該の連続動作区間において暖房連続動作区間であると特定された時間帯の各計測時刻ｔについて、即ち暖房用ガス器具が使用されている状態での、給湯用ガス器具の使用の有無を判断するための処理である。
給湯最小流速Ｖw［Ｌ／１０秒］は、給湯用ガス器具が使用されているときのガス消費量として想定される範囲の最小値に基づいて設定され、特定の値に限定されるものではなく、例えば、対象住宅・住戸に設置されている（若しくは、設置されていると想定される一般的な）給湯用ガス器具の仕様などが考慮されたり、或いは、実測データが収集されている場合には当該実測データが考慮されたりして適当な値に適宜設定される。具体的には例えば、あくまで一例として挙げると、０．５〜１．５［Ｌ／１０秒］程度の範囲で適当な値に設定され得る。
そして、当該の連続動作区間における暖房連続動作区間内の計測時刻ｔ毎に、ｙt≧Ｖhb＋Ｖwである場合には（Ｓ４−４：Ｙｅｓ）推定流速ｙtから暖房ベース流速Ｖhbを引いたものが給湯分のガス使用量に分類されると共に暖房ベース流速Ｖhbが暖房分のガス使用量に分類され、ｙt≧Ｖhb＋Ｖwでない場合には（Ｓ４−４：Ｎｏ）推定流速ｙtが暖房分のガス使用量に分類される。
一方、Ｓ４−１の処理において、当該の連続動作区間が冬季におけるものでない場合には（Ｓ４−１：Ｎｏ）Ｓ４−５の処理に進む。そして、計測時刻ｔのそれぞれについてＳ４−５乃至Ｓ４−７の処理が実行される。
また、Ｓ４−２の処理において、Ｖhs-min≦ｙt≦Ｖhs-max でない場合、若しくは、所定流速継続時間がＴh［分］以上でない場合にも（Ｓ４−２：Ｎｏ）Ｓ４−５の処理に進む。言い換えると、暖房連続動作区間ではない計測時刻ｔのそれぞれについてＳ４−５乃至Ｓ４−７の処理が実行される。
まず、推定流速ｙtが給湯最小流速Ｖw以上であるか否かが判断される（Ｓ４−５）。
Ｓ４−５の判断処理は、給湯用ガス器具の使用の有無を判断するための処理である。
そして、ｙt≧Ｖwである場合には（Ｓ４−５：Ｙｅｓ）推定流速ｙtが給湯分のガス使用量に分類される。
一方、ｙt≧Ｖwでない場合には（Ｓ４−５：Ｎｏ）Ｓ４−６の処理に進み、当該の連続動作区間が夏季におけるものであるか否かが判断される。
夏季であるか否かの判断（Ｓ４−６）の仕方は、特定の判断基準に限定されるものではなく、夏季とその他の季節とにおける例えば暖房器具の使用有無や給湯器具の温度設定などの差異が顕れる（若しくは、顕れると考えられる）時期などが考慮された上で適当な方法に適宜設定される。具体的には例えば、あくまで一例として挙げると、１１月から翌年３月までの期間に取得された計測データから抽出された連続動作区間のデータは夏季であると判断する方法が採用され得る。
なお、冬季や夏季であるか否かを期間によって判断するようにした場合に、冬季と夏季として指定される期間は、切れ目無く連続しているようにしても良いし、冬季にも夏季にも指定（分類）されない期間があっても良く、また、冬季として指定される期間と夏季として指定される期間とで一部が重複しているようにしても良い。
そして、当該の連続動作区間が夏季におけるものである場合には（Ｓ４−６：Ｙｅｓ）推定流速ｙtが厨房分のガス使用量に分類される。
一方、当該の連続動作区間が夏季におけるものでない場合には（Ｓ４−６：Ｎｏ）Ｓ４−７の処理に進み、推定流速ｙtが厨房最小流速Ｖkよりも大きいか否かが判断される。
Ｓ４−７の判断処理は、厨房用ガス器具の使用の有無を判断するための処理である。
厨房最小流速Ｖk［Ｌ／１０秒］は、厨房用ガス器具が使用されているときのガス消費量として想定される範囲の最小値に基づいて設定され、特定の値に限定されるものではなく、例えば、対象住宅・住戸に設置されている（若しくは、設置されていると想定される一般的な）厨房用ガス器具の仕様などが考慮されたり、或いは、実測データが収集されている場合には当該実測データが考慮されたりして適当な値に適宜設定される。具体的には例えば、あくまで一例として挙げると、０．２〜１．０［Ｌ／１０秒］程度の範囲で適当な値に設定され得る。
ここで、本実施形態のＳ４−７の処理では、厨房用ガス器具のみが使用されている状態と暖房用ガス器具が保温のために間欠的に動作している状態とが区別されるようにしている。
そこで、厨房最小流速Ｖkの値は、暖房用ガス器具が保温のために間欠的に動作しているときのガス消費量として想定される範囲の最大値（「間欠暖房最大流速」と呼ぶ）が合わせて考慮されたり当該間欠暖房最大流速に基づいて設定されたりするようにしても良い。
そして、推定流速ｙtが厨房最小流速Ｖkよりも大きい場合には（Ｓ４−７：Ｙｅｓ）推定流速ｙtが厨房分のガス使用量に分類される。
一方、推定流速ｙtが厨房最小流速Ｖkよりも大きくない場合には（Ｓ４−７：Ｎｏ）推定流速ｙtが暖房分のガス使用量に分類される。
そして、用途分解部１１ｄにより、ガス器具毎のガス使用量の例えば３０分間や１時間など時間，日，週間，月，季節，及び／又は年などの単位での積算値［Ｌ］が算出され、データファイル等に記録されて記憶部１２などに格納(保存)されたり表示部１４に表示されたりする。
なお、Ｓ４の処理において用いられる各閾値などの種々の値は、ガス消費量の用途分解プログラム１７内に予め規定されるようにしても良いし、Ｓ４の処理が行われる前に入力部１３を介して作業者によって入力されるようにしても良い。
そして、制御部１１は、Ｓ１の処理においてメモリ１５に記憶された、ガス消費量の用途分解のひとかたまりの処理対象としての計測データに関する処理を終了する（ＥＮＤ）。
以上のように構成されたガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムによれば、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量ｆu［Ｌ］毎に１カウントとして計数する流量カウント数との組み合わせデータ、すなわち、所定の時間間隔且つ所定の流量ｆu刻みという値の変化が散発的（言い換えると、間欠的）であり且つ離散的であるデータ、或いは言い換えると、値の変化の生起が時系列において散発的（間欠的）であり且つ値の変化の態様が離散的であるデータに基づいて計測時刻毎のガスの流速ｙtを算定することができるので、ガスの総消費量の実態の再現性の向上を図ることが可能になる。
以上のように構成されたガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムによれば、さらに、実態が良好に再現されたガスの流速ｙtに基づいてガスの総消費量の用途を給湯分，暖房分，及び厨房分に分類することができるので、ガスの総消費量の用途分解の良好な精度を確保してガス消費量の用途分解技術としての信頼性の向上を図ることが可能になる。
なお、上述の形態は本発明を実施する際の好適な形態の一例ではあるものの本発明の実施の形態が上述のものに限定されるものではなく、本発明の要旨を逸脱しない範囲において本発明は種々変形実施可能である。
例えば、上述の実施形態ではＳ４の処理においてガス総消費量が用途毎のガス使用量に分解される際に図４に示すアルゴリズムが用いられるようにしているが、用途分解に用いられるアルゴリズムは図４に示すものに限られるものではない。具体的には、本発明におけるガス総消費量の用途分解の要点は、
１）当該の連続動作区間が冬季であるか否かが判断される処理（上述の実施形態ではＳ４−１の処理が該当する）
２）主幹におけるガスの流速(ｙt)が所定の範囲である状態が所定の時間(Ｔh)以上継続しているか否かに基づいて暖房連続動作区間が特定される処理（上述の実施形態ではＳ４−２の処理が該当する）
３）当該の連続動作区間において特定された全ての暖房連続動作区間の中でのガスの流速(ｙt)の最小値が暖房ベース流速(Ｖhb)として設定される処理（上述の実施形態ではＳ４−３の処理が該当する）
４）暖房ベース流速(Ｖhb)と給湯最小流速(Ｖw)との和が考慮される処理（上述の実施形態ではＳ４−４の処理が該当する）
５）給湯最小流速(Ｖw)が考慮される処理（上述の実施形態ではＳ４−５の処理が該当する）、並びに、厨房最小流速(Ｖk)と間欠暖房最大流速とのうちの少なくとも一方が考慮される処理（上述の実施形態ではＳ４−７の処理が該当する）
が含まれることである。
そして、上記１乃至４の処理が実行される順序は図４に示すものに限られるものではない。
また、上記１乃至４として挙げた要点の列挙からも分かるように、当該の連続動作区間が夏季であるか否かが判断される処理（上述の実施形態ではＳ４−６の処理が該当する）は本発明において必須の処理ではない。例えば、上述の実施形態のようにＳ４−１の処理において当該の連続動作区間が冬季であるか否かを期間によって判断するようにした場合の冬季として指定される期間を調整することにより、当該の連続動作区間が夏季であるか否かが判断される処理は含まれないようにしても良い。
また、各処理において用いられる各閾値も、特定の値に限定されるものではなく、適宜設定され得る。
次に、図５乃至図１０に、本発明の居住者の在宅不在の判定方法、判定装置、及び判定プログラムの実施形態の一例を示す。
ここで、本発明における在宅不在の判定とは、住居や住戸内に居住者が存在している状態（即ち、在宅である）と存在していない状態（即ち、不在である）とのどちらであるかを推定することである。
本発明における居住者の在宅不在の判定でも、図３に示される連続データの作成方法の考え方が利用される。そして、所定の時間間隔且つ所定の計量単位で出力される計測値から実際の電力需要量が推定され、当該推定された電力需要量が用いられて居住者の在宅不在の判定が行われる。
具体的には、本実施形態の居住者の在宅不在の判定方法は、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量ｌを計測単位として計量する電力需要量との組み合わせデータとして取得（Ｓ１）された電力需要量の計測データから、所定の時間帯の組み合わせデータが分析対象区間のデータとして抽出され、分析対象区間における計測時刻毎の電力需要量の累積値を下限とすると共に計測時刻毎の電力需要量の累積値に所定の電力量ｌを加えた値を上限とする領域内に存在し且つ分析対象区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測された電力需要量の総和と推定される各計測時刻の電力需要量の総和とが一致する折れ線が推定され（Ｓ２）、当該折れ線における各計測時刻間の直線の傾きとして電力需要速度ｘtが算定され（Ｓ２）、在宅不在の判定の対象である計測時刻以前の所定の期間に含まれる電力需要速度ｘtが用いられて前記所定の期間における電力需要速度ｘtの統計量が特徴量として算出され、少なくとも前記特徴量が用いられて居住者の在宅不在が判定される（Ｓ３）ようにしている（図５参照）。
また、本実施形態の居住者の在宅不在の判定装置は、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量ｌを計測単位として計量する電力需要量との組み合わせデータとして取得された電力需要量の計測データの入力を受ける手段としてのデータ受部３１ａと、計測データから、所定の時間帯の組み合わせデータを分析対象区間のデータとして抽出し、分析対象区間における計測時刻毎の電力需要量の累積値を下限とすると共に計測時刻毎の電力需要量の累積値に所定の電力量ｌを加えた値を上限とする領域内に存在し且つ分析対象区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測された電力需要量の総和と推定される各計測時刻の電力需要量の総和とが一致する折れ線を推定して当該折れ線における各計測時刻間の直線の傾きとして電力需要速度を算定する手段としての最短折線作成部３１ｂと、在宅不在の判定の対象である計測時刻以前の所定の期間に含まれる電力需要速度を用いて前記所定の期間における電力需要速度の統計量を特徴量として算出すると共に少なくとも前記特徴量を用いて居住者の在宅不在を判定する手段としての時刻別判定部３１ｃとを有する。
上記居住者の在宅不在の判定方法及び居住者の在宅不在の判定装置は、居住者の在宅不在の判定プログラムがコンピュータ上で実行されることによっても実施・実現され得る。ここでは、居住者の在宅不在の判定プログラムがコンピュータ上で実行されることによって居住者の在宅不在の判定方法が実施されると共に居住者の在宅不在の判定装置が実現される場合を説明する。
本実施形態の居住者の在宅不在の判定プログラム３７を実行するためのコンピュータ３０（本実施形態では、居住者の在宅不在の判定装置３０でもある）の全体構成を図７に示す。
このコンピュータ３０（居住者の在宅不在の判定装置３０）は制御部３１，記憶部３２，入力部３３，表示部３４，及びメモリ３５を備え、これらが相互にバス等の信号回線によって接続されている。
制御部３１は、記憶部３２に記憶されている居住者の在宅不在の判定プログラム３７に従ってコンピュータ３０全体の制御並びに居住者の在宅不在の判定に係る演算を行うものであり、例えばＣＰＵ（中央演算処理装置）である。
記憶部３２は、少なくともデータやプログラムを記憶可能な装置であり、例えばハードディスクである。
入力部３３は、少なくとも作業者の命令や種々の情報を制御部３１に与えるためのインターフェイス（即ち、情報入力の仕組み）であり、例えばキーボードやマウスである。なお、例えばキーボードとマウスとの両方のように複数種類のインターフェイスを入力部３３として有するようにしても良い。
表示部３４は、制御部３１の制御によって文字や図形或いは画像等の描画・表示を行うものであり、例えばディスプレイである。
メモリ３５は、制御部３１が種々の制御や演算を実行する際の作業領域であるメモリ空間となるものであり、例えばＲＡＭ（Ｒandom Ａccess Ｍemory の略）である。
そして、コンピュータ３０（以下、「居住者の在宅不在の判定装置３０」と表記する）の制御部３１には、居住者の在宅不在の判定プログラム３７が実行されることにより、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量ｌを計測単位として計量する電力需要量との組み合わせデータとして取得された電力需要量の計測データの入力を受ける処理を行うデータ受部３１ａと、計測データから所定の時間帯の組み合わせデータを分析対象区間のデータとして抽出し、さらに、分析対象区間における計測時刻毎の電力需要量の累積値を下限とすると共に計測時刻毎の電力需要量の累積値に所定の電力量ｌを加えた値を上限とする領域内に存在し且つ分析対象区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測された電力需要量の総和と推定される各計測時刻の電力需要量の総和とが一致する折れ線を推定して当該折れ線における各計測時刻間の直線の傾きとして電力需要速度ｘtを算定する処理を行う最短折線作成部３１ｂと、在宅不在の判定の対象である計測時刻以前の所定の期間に含まれる電力需要速度ｘtを用いて所定の期間における電力需要速度ｘtの統計量を特徴量として算出すると共に少なくとも特徴量を用いて居住者の在宅不在を判定する処理を行う時刻別判定部３１ｃとが構成される。
ここで、以下に説明する処理を行う際に用いられる変数についての計算機における処理として、最初に、不在開始時刻，不在終了時刻，及び不在区間指標を表す変数は全て初期化され、具体的には、不在開始時刻及び不在終了時刻を表す変数には例えば時刻としてはあり得ない数値が与えられたり、不在区間指標を表す変数には例えば０(ゼロ)が与えられたりする。
そして、居住者の在宅不在の判定方法の実施として、まず、電力需要量の計測が行われて計測データの取得が行われる（Ｓ１）。
電力需要量の計測の仕方は、在宅不在の判定の対象である住宅や住戸における電力需要量を計測して所定の計量単位（言い換えると、所定の計量の刻み）で所定の時間毎に例えば電気的な信号として前記所定の計量単位での需要量を出力するものであれば、特定の方法や機器によるものに限定されるものではなく、前記のような計測を行い得る方法や機器が適宜選択される。
具体的には例えば、あくまで一例として挙げると、電力需要量の計測の仕方（具体的には、計測機器）としてスマートメータ（例えば、上述の非特許文献２）が利用され得る。
電力需要量の計測における上述の「所定の計量単位」としての、１単位として扱われる（言い換えると、１単位として計量される）実際の電力量は、特定の値［Ｗｈ］に限定されるものではなく、例えば計測に用いられる計測機器の仕様などが踏まえられた上で適当な値に適宜設定される。具体的には例えば、スマートメータのＡルート計測値では、１単位として扱われる実際の電力量は１００ Ｗｈ である。
また、電力需要量の計測における上述の「所定の時間」としての、計測値を出力する時間間隔は、特定の時間［分］に限定されるものではなく、例えば計測に用いられる計測機器の仕様などが踏まえられた上で適当な時間に適宜設定される。具体的には例えば、スマートメータのＡルート計測値では、計測値が出力される時間間隔は３０分である。
本実施形態では、電力需要量の計測においてスマートメータが利用され、上述の「所定の計量単位」が１００ Ｗｈ であると共に上述の「所定の時間」が３０分である計測データ、すなわち、実際の電力量１００ Ｗｈ を１単位として計量する電力需要の値（「需要単位計量値」と呼ぶ）が３０分間隔で出力される計測データが、Ｓ１の処理によって取得されるとする。
そして、計測によって取得された電力需要量の計測値（具体的には、需要単位計量値）は、制御部３１のデータ受部３１ａを介して居住者の在宅不在の判定装置３０へと入力される。
電力需要量の計測値（需要単位計量値）は、居住者の在宅不在の判定装置３０に、各種記録媒体やデータサーバを介して入力されるようにしても良く、或いは、例えばスマートメータのように計測機器が備えている通信機能または計測機器に付加された通信機能や通信機序によって入力されるようにしても良い。
記録媒体を介して計測値が入力される場合には、例えば、計測機器から出力された需要単位計量値が記録媒体に記録され、当該記録媒体が居住者の在宅不在の判定装置３０の記録媒体接続用端子（図示していない）に差し込まれることによって需要単位計量値が入力されるようにすることが考えられる。
また、データサーバを介して計測値が入力される場合には、例えば、データサーバがバス等の信号回線によって居住者の在宅不在の判定装置３０に接続され、計測機器から出力された需要単位計量値が前記データサーバにデータファイル等として格納(保存)され、当該データファイル等として保存された需要単位計量値が読み込まれるようにすることが考えられる。
また、例えばスマートメータのように計測機器が備えている通信機能または計測機器に付加された通信機能や通信機序によって計測値が入力される場合には、需要単位計量値がリアルタイムで（即ち、即時的に）居住者の在宅不在の判定装置３０に入力されることも可能であり、したがってほぼリアルタイムでの（即ち、即時的な）居住者の在宅不在の判定の処理も可能である。
なお、記録媒体，データサーバ，及び通信機能や通信機序が組み合わされて用いられて需要単位計量値が居住者の在宅不在の判定装置３０に入力されるようにしても良い。
また、電力需要量の各計測値（各需要単位計量値）が計測機器から出力された時刻が計測機器に備えられた時計機能によって計測値と一緒に出力されて当該計測値と対応づけられて記録されたり、或いは、計測機器から出力された各計測値が記録媒体，データサーバ，或いは居住者の在宅不在の判定装置３０に記録されたり入力されたりした時の時刻が当該計測値と対応づけられて記録されたりする。なお、この時刻の間隔［分］が上述の「所定の時間」になる。
ここでの説明においては、各計測値（各需要単位計量値）と対応づけられて記録される時刻のことを「計測時刻」と呼び、また、需要単位計量値と計測時刻との組み合わせデータのことを「計測データ」と呼ぶ。
そして、データ受部３１ａにより、計測データがメモリ３５に記憶させられる。
次に、制御部３１の最短折線作成部３１ｂにより、Ｓ１の処理によって計測されて取得された計測データが用いられて累積需要曲線の作成が行われる（Ｓ２）。
Ｓ２の処理は、分析対象区間について実行される。分析対象区間は、少なくとも、対象住宅や対象住戸についての居住者の在宅不在の判定が行われる時刻（「判定対象時刻」と呼ぶ）を含む期間として設定される。
分析対象区間は、特定の時間長さ［時間］に限定されるものではなく、例えば検知の主要な対象として想定される不在期間について想定される時間が考慮されるなどして適当な時間長さに適宜設定される。分析対象区間は、具体的には例えば、あくまで一例として挙げると、６〜３６時間程度の範囲で適当な時間長さに設定され得る。また、分析対象区間は、例えばリアルタイムでの（即ち、即時的な）居住者の在宅不在の判定が行われる場合などには、判定対象時刻の進行に伴って当該の判定対象時刻を基準とするなどして適時・適宜に更新されて（言い換えると、設定し直されて）時間軸方向に移動させられるようにしても良い。
なお、判定対象時刻は、全ての計測時刻の各々が順次に判定対象時刻とされるようにしても良く、或いは、各計測時刻のうちから所定の間隔で選択された計測時刻が順次に判定対象時刻とされるようにしても良い。
Ｓ１の処理によって取得される計測データは、３０分間隔での１００ Ｗｈ 単位（言い換えると、１００ Ｗｈ 刻み）の電力需要量である。しかしながら、何らかの電気機器が使用されているときは、さらに、一般の住居では冷蔵庫など常時稼働且つ所定の間隔（或いは、所定のトリガー）で電力需要が変動する電気機器が通常は使用されていることを考慮すれば、電力（待機電力等を含む）の需要が常に継続しているという点において電力需要量の変化の生起は時系列において連続であり、また、電力が常に消費されているという点において電力需要量の変化の態様は漸次的である。
つまり、Ｓ１の処理における計測では、実際には(値の)変化の生起が時系列において連続であると共に(値の)変化の態様が漸次的である現象が、所定の時間間隔且つ所定の値単位（言い換えると、所定の値刻み，所定の計測単位，或いは所定の計量単位）という散発的（言い換えると、間欠的）且つ離散的な変化として把握されていると言える。
そして、本発明では、所定の時間間隔且つ所定の値単位（言い換えると、所定の値刻み，所定の計測単位，或いは所定の計量単位）という計測値の変化が散発的（言い換えると、間欠的）であり且つ離散的である計測データがそのまま用いられるのではなく、当該散発的（間欠的）且つ離散的な計測データが、計測値の変化の生起が時系列において連続であり且つ計測値の変化の態様が漸次的であるデータに変換された上で在宅不在の判定の処理に用いられる。
本発明の居住者の在宅不在の判定方法、判定装置、及び判定プログラムにおいて用いられる、計測値の変化の生起が時系列において散発的（間欠的）であり且つ計測値の変化の態様が離散的であるデータに基づいて、計測値の変化の生起が時系列において連続であり且つ計測値の変化の態様が漸次的である累積需要曲線が作成される、連続データの作成方法を、図８を用いて説明する。そして、この作成方法の考え方（言い換えると、要点）は、図３を用いて説明した作成方法の考え方（要点）と同様である。
ここでの説明の例として、図８に示すような、或る時間帯の計測データが、具体的には例えばスマートメータのＡルート計測値のように３０分毎に出力される需要単位計量値の累積値であって１００ Ｗｈ 単位の電力需要量の累積値（図中の●印）がＳ１の処理によって取得されて分析対象区間分の計測データとして抽出されたとする。
図８の横軸は、３０分間隔である計測時刻の各時点を、当該の分析対象区間の最初の計測時刻を１とした上で識別子としての連続番号で表したものである。すなわち、横軸における例えば「３」は、当該の分析対象区間の始まりから１．５時間（＝９０分＝３０分×３）だけ時間が経過していることを表す。
また、図３の縦軸は、当該の分析対象区間における需要単位計量値の累積値を表す。すなわち、縦軸における例えば「４００」は、当該の分析対象区間の始まりから４００［Ｗｈ］だけ電力が消費されていることを表す。
Ｓ１の処理によって取得される計測データは所定の時間間隔且つ所定の値単位（値刻み，計測単位，計量単位）という値の変化が散発的（間欠的）且つ離散的なデータであるので、すなわち、計測値の変化の生起が時系列において散発的（間欠的）且つ計測値の変化の態様が離散的なデータであるので、計測データの需要単位計測値の累積値のデータプロットである●印を時系列で結ぶと図８に示すように不規則な昇り階段状になる（図８中において符号４１；「下限階段直線４１」と呼ぶ）。
また、計測データの需要単位計量値の累積値に計測時刻毎に１００ Ｗｈ（即ち、需要単位計量値における計量単位であり値刻みの大きさである）を加えたデータのプロットである▲印を時系列で結んだものも、図８に示すように不規則な昇り階段状になる（図８中において符号４２；「上限階段直線４２」と呼ぶ）。
そして、下限階段直線４１が実際の電力需要量に基づく真の累積値の下限（言い換えると、実際の電力需要量が取り得る累積値の下限）であり、上限階段直線４２が実際の電力需要量に基づく真の累積値の上限（言い換えると、実際の電力需要量が取り得る累積値の上限）である。
したがって、実際の電力需要量の時系列における累積を表す真の累積需要曲線は、下限階段直線４１と上限階段直線４２とに挟まれた領域（「存在可能領域」と呼ぶ）内に存在し、全体としては増加傾向を示す非減少の曲線であると考えられる。
そして、真の累積需要曲線としては、存在可能領域内において増加傾向を示す非減少の折れ線のうち、折れ点（言い換えると、変曲点、若しくは、電力需要速度変化点）の個数が比較的少ないものを想定する。
そこで、実際の電力需要量の時系列における累積を表す真の累積需要曲線として、存在可能領域内に存在する、総長（即ち、当該の分析対象区間の始まりから終わりまでに亙る軌跡の長さ）が最短の折れ線（「最短折れ線」と呼ぶ）を推定する。
真の累積需要曲線としての、存在可能領域内に存在する、最短折れ線の推定手法は、特定の方法に限定されるものではなく、特定の領域内に存在して所定の始点から終点へと至る折れ線のうちの総長が最短の折れ線であり且つ各計測時刻において計測された電力需要量の総和と推定される各計測時刻の電力需要量の総和とが一致する折れ線を推定し得る適当な方法が適宜選択される。
最短折れ線の推定手法として、具体的には例えば、あくまで一例として挙げると、当該の分析対象区間において確定済みの折れ線の最後の点を起点とした（未確定の）直線部の許容される傾きの上限値・下限値を順次更新し、全ての直線が存在可能領域を逸脱せざるを得なくなる時点で直線の傾きを決定し、その直線が存在可能領域の境界に最初に接した時点までの折れ線の確定を行い、そして、新たに確定した折れ線の最後の点を新たな起点に設定した上で上述と同様の処理を繰り返すことを基本的な考え方とする方法が用いられ得る。
この推定方法は、計測時刻ｔ（尚、計測単位は３０分であり、したがって計測時刻の間隔は３０分である）における計側電力量をｘt［Ｗｈ／３０分］とすると共に、推定電力需要速度をｙt［Ｗｈ／３０分］とすると（図９参照）、最適化問題として数式７のように表される。
数式７においてｔは分析対象区間における計測時刻の識別子であってｔ＝１，２，３，…であり、数式７では当該の分析対象区間の範囲として計測時刻の識別子ｔが１からＴまでの期間に亙って処理が実行される。
そして、計測時刻を１刻みの識別子とすることにより、計測の時間間隔である３０分を一単位として扱うことになり、計測時刻 ｔ−１ から ｔ にかけての増分は計測の時間間隔当たりの増加量であって計測時刻 ｔ−１ から ｔ にかけての速度［Ｗｈ／３０分］であり、真の累積需要曲線としては計測時刻 ｔ−１ から ｔ にかけての直線の傾きである。
数式７において、また、計測電力量ｘt［Ｗｈ／３０分］は、計測時刻ｔにおける需要単位計量値（言い換えると、計測時刻 ｔ−１ から ｔ にかけての需要単位計量値の増分）である。
数式７において、さらに、ｌ(エル)は、電力需要［Ｗｈ］の最小単位（即ち、計測単位，計量単位）であり、本実施形態ではｌ＝１００である。
そして、数式７における最小化の対象は、経過時間（即ち、各々１）と各経過時間に関する推定需要（即ち、推定電力需要速度ｙt）とからなる折れ線の距離である。折れ線の距離は、図９に示す例のように、経過時間（即ち、各々１）と当該経過時間に関する推定電力需要速度ｙtとに対する斜辺の長さの総和として求められ、本発明ではこの斜辺の長さの総和が最も小さくなる曲線が推定累積需要曲線とされる。
制約条件の第一式は電力需要速度の非負性に対応し、また、第二式は実際の電力需要量の累積値の範囲制約に対応し、さらに、第三式は当該の分析対象区間における計測値ベースの電力需要量の総和と推定値ベースの電力需要量の総和との一致に対応する。
上述の最適化問題を解く方法によると、最悪の場合でも、時系列長（即ち、計測時刻の個数）の二乗のオーダーで真の累積需要曲線に相当する最短折れ線が推定される。
図８に示す計測データのプロット（図中の●印）を結んだ下限階段直線４１に対して数式７で表される最適化問題を適用した結果として（言い換えると、数式７に示す三つの制約条件を満たす累積需要曲線の全体が）、図８において破線で示す、存在可能領域内の単調増加曲線として、実際の電力需要量に近いと考えられる推定累積需要曲線４３が得られる。
ここで、上記の最適化問題は計測時刻の識別子ｔ＝１から処理が実行されるところ、ｔ＝１の処理が実行されるためには、ｔ＝１での起点としてｔ＝０（即ち、計測時刻の識別子ｔ＝１に対応する計測時刻の一つ前の計測時刻）からｔ＝１にかけての折れ線の最後の点が特定される必要が、すなわち、ｔ＝０からｔ＝１にかけての直線が特定される必要がある。この直線は、具体的には、当該の分析対象区間の最初の計測時刻（即ち、計測時刻の識別子ｔ＝１）の直前の計測時刻として識別子ｔを０(ゼロ)とすると共に需要単位計量値の累積値を０(ゼロ)とする仮想の原点を置き、当該原点と計測時刻の識別子ｔ＝１における需要単位計量値の累積値の値の点とを結ぶことによって特定される。
あるいは、当該の分析対象区間についての計測時刻の識別子ｔ＝０からｔ＝１にかけての直線の傾きが、識別子ｔ＝１において連続する直線の傾きに一致させられるようにしても良い。
ここで、真の累積需要曲線の上述の推定方法（即ち、数式７で表される最適化問題）は、目的関数部を「距離の和」から「距離の２乗和」に変更した数式８で表される最適化問題、すなわち、「存在可能領域の任意の曲線のうち各計測時刻の電力需要速度（即ち、電力需要量の累積値の計測時刻間差分）の２乗和が最小の曲線」を求める問題と等価である。
数式８で表される最適化問題では、２乗和の最小化により、各計測時刻の電力需要速度（電力需要量）が同じ値を持つ傾向が生じるので、折れ線が最適解になる。
真の累積需要曲線の推定方法としての数式７で表される最適化問題は、また、数式８で表される最適化問題の双対問題である数式９で表される最適化問題とも等価である。なお、数式９で表される最適化問題は、全変動正規化問題と呼ばれる凸最適化問題の一種である。なお、数式９において、ｘ0は計測電力量ｘt［Ｗｈ／３０分］の初期値（即ち、ｔ＝０のときの値）であり、図８，図９に示す例ではｘ0＝０である。
上述の真の累積需要曲線（即ち、最短折れ線）は、一例として以下に示すアルゴリズムが用いられることによって求められ得る。
−真の累積需要曲線の決定アルゴリズム−
《保持・更新情報》
○起点 折れ線
○上限直線 上限需要と領域境界接触点
○下限直線 下限需要と領域境界接触点
《更新内容》
(1) 時刻更新：時刻＝時刻＋１
(2) 下限直線＜上限累積需要値
起点から下限境界接触点までの需要量＝需要下限（確定）
起点＝下限直線の境界接触点
時刻＝起点，上限／下限直線の再設定
(1)へ
(3) 上限直線＞下限累積需要
起点から上限境界接触点までの需要量＝需要上限（確定）
起点＝上限直線の境界接触点
時刻＝起点，上限／下限直線の再設定
(1)へ
(4) 下限直線＜領域下限
下限需要，下限境界設定更新
上限直線＞領域上限
上限需要，上限境界設定更新
(5) 計測時刻の識別子ｔ＝Ｔになるまで(1)へ戻る
本実施形態では、最短折線作成部３１ｂにより、Ｓ１の処理においてメモリ３５に記憶された計測データ（即ち、需要単位計量値と計測時刻との組み合わせデータ）のうち、分析対象区間に該当する期間の計測データが時系列に従って順次読み込まれる。
そして、メモリ３５から読み込まれた計測時刻は、当該の分析対象区間における計測時刻の識別子として時系列順に１，２，３…という連続番号に変換される。
続いて、最短折線作成部３１ｂにより、上記の変換が施された計測データが用いられて数式７乃至数式９に表される最適化問題のうちのいずれかが解かれて、計測時刻の識別子ｔ毎の推定電力需要速度ｙt［Ｗｈ／３０分］が算定される。
そして、最短折線作成部３１ｂにより、当該の分析対象区間について、算定された推定電力需要速度ｙt［Ｗｈ／３０分］の値が計測時刻の識別子ｔの値と対応づけられて（言い換えると、計測時刻の識別子ｔの値と推定電力需要速度ｙtの値との組み合わせデータとして）メモリ３５に記憶させられる。
次に、制御部３１の時刻別判定部３１ｃにより、Ｓ２の処理によって算定された計測時刻別の推定電力需要速度が用いられて居住者の在宅不在の時刻別判定指標の決定が行われる（Ｓ３）。
Ｓ３の処理は、Ｓ２の処理が実行された分析対象区間に含まれる判定対象時刻の各々について実行される。なお、判定対象時刻に該当する計測時刻の識別子をｔeと表記する。
本発明では、判定対象時刻における在宅不在の判定を行う際に、当該の判定対象時刻以前の一定期間内の推定電力需要速度のデータが用いられて平均値や最大値などの統計量が特徴量として算出され、当該特徴量が用いられて在宅か不在かの判定が行われる。特徴量の算出対象期間であって、当該の判定対象時刻以前の一定期間のことを「特徴量期間」と呼ぶ。
特徴量期間は、特定の時間長さ［分，時間］に限定されるものではなく、例えば事前の分析結果が考慮されるなどして適当な時間長さに適宜設定される。特徴量期間は、具体的には例えば、あくまで一例として挙げると、１〜３時間程度の範囲で適当な時間長さに設定され得る。
具体的には例えば、或る判定対象時刻（計測時刻の識別子ｔ＝ｔe）における判定を行う場合に、特徴量期間が１時間であるときは、当該の判定対象時刻（計測時刻の識別子ｔ＝ｔe）における推定電力需要速度ｙte［Ｗｈ／３０分］（尚、当該の判定対象時刻以前３０分間の状態が反映されている）と当該の判定対象時刻の一つ前の計測時刻（計測時刻の識別子ｔ＝ｔe−１）における推定電力需要速度ｙte-1［Ｗｈ／３０］（尚、当該の判定対象時刻の６０分前から３０分前までの状態が反映されている）とが用いられてこれら二つの値についての統計量が特徴量として算出される。
特徴量としての統計量は、特定の種類に限定されるものではなく、例えば事前の分析結果が考慮されるなどして適当な種類が取り上げられて適宜設定される。特徴量としては、具体的には例えば、平均値（即ち、特徴量期間に含まれる計測時刻別の推定電力需要速度の平均の値），最大値（即ち、特徴量期間に含まれる計測時刻別の推定電力需要速度のうちの最も大きな値），最小値（即ち、特徴量期間に含まれる計測時刻別の推定電力需要速度のうちの最も小さな値），範囲（即ち、特徴量期間に含まれる計測時刻別の推定電力需要速度のうちの最も大きな値と最も小さな値との差の値），標準偏差（即ち、特徴量期間に含まれる計測時刻別の推定電力需要速度の標準偏差の値）などが挙げられる。
なお、Ｓ３の処理において用いられる特徴量の個数は、特定の数に限定されるものではなく、一個でも良く、或いは、複数個でも良い。
また、特徴量として、電力需要量（具体的には、推定電力需要速度）に関する各種統計量に加え、種々の情報が追加的に用いられるようにしても良い。この場合の情報の種類としては、例えば、世帯の識別子，判定対象時刻（実際の時刻），気温，季節などが挙げられる。
そして、上記のようにして算出されたり設定されたりした特徴量（言い換えると、説明変数）が用いられて、在宅と不在との２クラス分類が行われることによって在宅不在の判定が行われる。
２クラス分類を行う手法としては、例えばランダムフォレストやサポートベクターマシン（ＳＶＭ：Ｓupport Ｖector Ｍachine の略）などの機械学習手法が用いられ得る。
なお、機械学習手法が用いられる際には、時刻別の電力需要量の実測データと当該時刻における実際の在宅不在に関する状態データとの組み合わせデータが正解データ（「教師データ」や「訓練データ」とも呼ばれる）が必要とされる。この正解データは、例えば、実態調査，モニタ調査，アンケート調査などが行われて収集されたり、既存の観測データが利用されて準備されたりする。そして、正解データが機械学習手法における学習に用いられて２クラス分類の判定基準が作成される。なお、２クラス分類に用いられる機械学習手法の手順や判定基準の内容などは、例えば、居住者の在宅不在の判定プログラム３７内に予め規定される。
本実施形態では、時刻別判定部３１ｃにより、Ｓ２の処理においてメモリ３５に記憶された分析対象区間についての計測時刻の識別子ｔの値と推定電力需要速度ｙtの値との組み合わせデータのうち、当該の判定対象時刻についての特徴量期間に該当する（言い換えると、特徴量期間に含まれる）計測時刻の識別子と対応づけられている全ての推定電力需要速度の値が読み込まれる。
そして、時刻別判定部３１ｃにより、読み込まれた推定電力需要速度の値が用いられて、２クラス分類に用いられる特徴量が算出される。
時刻別判定部３１ｃにより、さらに、機械学習手法が用いられて２クラス分類が行われ、当該の判定対象時刻について対象住宅や対象住戸に人が存在しているか存在していないかが決定される。
そして、時刻別判定部３１ｃにより、２クラス分類の結果（即ち、「在宅」と「不在」とのうちのどちらであるかの決定結果）が、在宅不在の時刻別判定指標として判定対象時刻（計測時刻の識別子ｔ＝ｔe）と対応づけられて（言い換えると、計測時刻の識別子の値と在宅不在の時刻別判定指標を表す変数の値との組み合わせデータとして）メモリ３５に記憶させられる。
なお、計算機における処理としては、２クラス分類の結果が、「在宅」である場合には時刻別判定指標を表す変数に「在宅」に対応する数値として適宜選択された適当な数値（具体的には例えば、１）が与えられ、「不在」である場合には時刻別判定指標を表す変数に「不在」に対応する数値として適宜選択された適当な数値（具体的には例えば、０(ゼロ)）が与えられる。
本実施形態では、引き続いて、上述までの処理によって決定された在宅不在の時刻別判定指標が用いられて不在開始時刻及び不在終了時刻が特定されて不在期間が確定されるようにしている（Ｓ４乃至Ｓ１１；図６，図１０）。
具体的にはまず、制御部３１の不在区間判断部３１ｄにより、Ｓ３の処理によって決定された当該の判定対象時刻における在宅不在の時刻別判定指標が用いられて、当該時刻別判定指標が「不在」であるか否かの判断が行われる（Ｓ４）。
本実施形態では、不在区間判断部３１ｄにより、Ｓ３の処理においてメモリ３５に記憶された当該の判定対象時刻の計測時刻の識別子ｔeと対応づけられている在宅不在の時刻別判定指標が読み込まれる。
そして、読み込まれた時刻別判定指標が「不在」である場合には（Ｓ４：Ｙｅｓ）、処理ステップがＳ５の処理へと進められる。
そして、Ｓ５の処理として、不在区間判断部３１ｄにより、後述するＳ８の処理において設定・操作される不在区間指標が「不在」であるか否かの判断が行われる。
そして、不在区間指標が「不在」でない場合には（Ｓ５：Ｎｏ）、不在開始時刻がＳ１の処理の前の初期化された状態のままであるときは、当該の判定対象時刻の計測時刻の識別子ｔ（＝ｔe）の一つ前の計測時刻の識別子ｔ（＝ｔe−１）に対応する実際の時刻が不在開始時刻としてメモリ３５に記憶させられる（Ｓ６）。なお、不在開始時刻がＳ１の処理の前の初期化された状態のままでないときは、不在開始時刻に対しては何らの操作もしない。
次に、制御部３１の在宅不在判定部３１ｅにより、Ｓ３の処理によって決定された在宅不在の時刻別判定指標が用いられて、居住者の在宅不在の判定が行われる（Ｓ７）。
Ｓ３の処理によって決定される在宅不在の時刻別判定指標は、判定対象時刻毎に各々独立して判定された結果であり、居住者が在宅しているものの、短時間、使用する電気機器が少なかったり電気機器の使用態様に変化が無かったりする場合に（言い換えると、電力需要速度の変化が小さい場合に）、不在であるとの判断結果になる虞がある。
このため、本実施形態では、居住者の在宅不在の最終的な判定は、在宅不在の時刻別判定指標が当該の判定対象時刻以前の所定の時間に亙って連続して「不在」であるか否か（言い換えると、Ｓ３の処理における「不在」の決定が所定の時間に亙って継続しているか否か）が判断されることによって行われる。ここで、在宅不在の判定を行うための前記所定の時間のことを「判定期間」と呼ぶ。
判定期間は、特定の時間長さ［時間］（言い換えると、判定対象時刻の個数）に限定されるものではなく、例えば事前の分析結果が考慮されるなどして適当な時間長さに適宜設定される。判定期間は、具体的には例えば、１〜１２時間程度の範囲で適当な時間長さに設定されることが考えられ、４〜８時間程度のうちのいずれかの時間長さに設定されることが好ましく、６時間程度に設定されることが更に好ましい。
本実施形態では、在宅不在判定部３１ｅにより、Ｓ３の処理においてメモリ３５に記憶された計測時刻の識別子ｔの値と在宅不在の時刻別判定指標を表す変数の値との組み合わせデータのうち、当該の判定対象時刻の計測時刻の識別子ｔ（＝ｔe）と対応づけられている在宅不在の時刻別判定指標を表す変数の値が読み込まれると共に、当該の判定対象時刻についての判定期間に該当する（言い換えると、判定期間に含まれる）計測時刻の識別子ｔ（＝ｔe−１，ｔe−２ 等）と対応づけられている全ての在宅不在の時刻別判定指標を表す変数の値が読み込まれる。
そして、在宅不在判定部３１ｅにより、当該の判定対象時刻についての判定期間中の在宅不在の時刻別判定指標が全て「不在」であるか否かが判断される。
そして、判定期間中の在宅不在の時刻別判定指標の中に一つでも「在宅」が有る場合には（Ｓ７：Ｎｏ）、処理ステップがＳ１の処理へと戻される。
他方、判定期間中の在宅不在の時刻別判定指標が全て「不在」である場合には（Ｓ７：Ｙｅｓ）、当該の判定期間の間は居住者は不在であると判定されるということになり、在宅不在判定部３１ｅにより、不在区間指標が「不在」にされた上でメモリ３５に記憶させられると共に当該の判定対象時刻の計測時刻の識別子ｔ（＝ｔe）に対応する実際の時刻が不在終了時刻としてメモリ３５に記憶させられ（Ｓ８）、その上で処理ステップがＳ１の処理へと戻される。なお、不在区間指標が「不在」にされるとは、計算機における処理としては、不在区間指標を表す変数に「不在」に対応する数値として適宜選択された適当な数値、具体的には例えば１が与えられたりすることである。
また、Ｓ５の処理において、不在区間判断部３１ｄにより、不在区間指標が「不在」であるか否かの判断が行われて当該不在区間指標が「不在」である場合には（Ｓ５：Ｙｅｓ）、当該の判定対象時刻の計測時刻の識別子ｔ（＝ｔe）に対応する実際の時刻が不在終了時刻としてメモリ３５に記憶させられ（Ｓ９）、その上で処理ステップがＳ１の処理へと戻される。なお、Ｓ９の処理では、Ｓ８の処理若しくは前回のＳ９の処理において設定された不在終了時刻が更新されることになる。
また、Ｓ４の処理において、不在区間判断部３１ｄにより、当該の判定対象時刻における在宅不在の時刻別判定指標が「不在」であるか否かの判断が行われて当該時刻別判定指標が「不在」でない場合には（Ｓ４：Ｎｏ）、処理ステップがＳ１０の処理へと進められる。
そして、Ｓ１０の処理として、制御部３１の不在期間確定部３１ｆにより、不在区間指標が「不在」であるか否かの判断が行われ、不在区間指標が「不在」である場合には（Ｓ１０：Ｙｅｓ）、上述のＳ６の処理において特定された不在開始時刻から、上述のＳ８の処理若しくはＳ９の処理において特定された不在終了時刻までが、不在期間であると確定される（Ｓ１１）。
そして、対象住宅や対象住戸についての在宅不在の判定の結果として、必要に応じて個々の不在期間を区別するための識別子と対応づけられて、不在開始時刻と不在終了時刻との組み合わせデータが、データファイル等に記録されて記憶部３２などに格納(保存)されたり、表示部３４に表示されたりする。その上で、処理ステップがＳ１の処理へと戻される。この際、不在開始時刻，不在終了時刻，及び不在区間指標は全てリセットされ、計算機における処理としては、Ｓ１の処理の前の初期化された状態に戻される。
他方、不在区間指標が「不在」でない場合には（Ｓ１０：Ｎｏ）、処理ステップがＳ１の処理へと戻される。この際、不在開始時刻がＳ１の処理の前の初期化された状態のままでないときは、不在開始時刻はリセットされ、計算機における処理としては、Ｓ１の処理の前の初期化された状態に戻される。
以上のように構成された居住者の在宅不在の判定方法、判定装置、及び判定プログラムによれば、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量ｌ［Ｗｈ］を計測単位として計量する電力需要量との組み合わせデータ、すなわち、所定の時間間隔且つ所定の電力量ｌ刻みという値の変化が散発的（言い換えると、間欠的）であり且つ離散的であるデータ、或いは言い換えると、値の変化の生起が時系列において散発的（間欠的）であり且つ値の変化の態様が離散的であるデータに基づいて計測時刻毎の電力需要速度ｙtを算定することができるので、電力需要の実態の再現性の向上を図ることが可能なる。
以上のように構成された居住者の在宅不在の判定方法、判定装置、及び判定プログラムによれば、さらに、実態が良好に再現された電力需要速度に基づいて住居や住戸に人が存在しているか存在していないかを推定することができるので、居住者の在宅不在の判定の良好な精度を確保して居住者の在宅不在の判定技術としての信頼性の向上を図ることが可能になる。
また、不在開始時刻から不在終了時刻までの期間は、つまりは対象住宅や対象住戸において人の活動に伴う電気の使用が為されていない期間であるので、対象住宅や対象住戸には人が存在していないと理解したり、或いは、対象住宅や対象住戸に人が存在しているものの例えば病気や怪我などのために動けなくなっているような異常事態が発生している可能性があると理解したりすることができる。したがって、不在開始時刻から不在終了時刻までの時間長さが所定の閾値を超える場合を要安否確認事象として検知し、対象住宅や対象住戸に関して長期不在の有無を確認したり異常事態の発生の有無を確認したりするようにすることも考えられる。すなわち、本発明の居住者の在宅不在の判定手法は、例えば、高齢者等の見守りシステムとしても利用されることが可能である。
なお、上述の形態は本発明を実施する際の好適な形態の一例ではあるものの本発明の実施の形態が上述のものに限定されるものではなく、本発明の要旨を逸脱しない範囲において本発明は種々変形実施可能である。
例えば、上述の実施形態ではＳ３の処理において居住者の在宅不在に関する時刻別判定指標が決定された上で更にＳ４乃至Ｓ１１の処理が実行されて不在開始時刻と不在終了時刻とが特定されるようにしているが、前記Ｓ４乃至Ｓ１１の処理は本発明において必須の構成ではない。すなわち、本発明の要点の一つは前記時刻別判定指標が決定されることであり、前記Ｓ４乃至Ｓ１１の処理は前記時刻別判定指標の利用の仕方の一つであるものの、前記時刻別判定指標の利用の仕方はこれに限定されるものではない。
上記のことも踏まえると、本発明の実施の形態としては、時刻別判定指標が決定されるまで（即ち、上述の実施形態におけるＳ３の処理まで）という態様や、上述の実施形態におけるＳ３の処理に続いて判定期間中の在宅不在の時刻別判定指標が全て「不在」であるか否かが判断される（即ち、上述の実施形態におけるＳ７の処理に相当する処理が実行される。この場合、判定対象時刻以前の判定期間中の在宅不在が判定され、これら処理が判定対象時刻の進行に伴って判定対象時刻毎に繰り返されることによって判定対象時刻以前の判定期間中の在宅不在の情報が更新される）という態様などが挙げられる。
また、上述の実施形態におけるＳ４乃至Ｓ１１の処理の要点は、Ｓ３の処理によって決定された在宅不在の時刻別判定指標を用いて不在期間を確定し得るように、在宅不在の判定を行うための所定の時間としての判定期間が判定対象時刻の進行に伴って時間軸方向に移動させられながら判定期間毎に在宅不在が判定され、不在と判定された一つの判定期間の始まりを不在開始時刻とすると共に終わりを不在終了時刻として特定する、または、連続して不在と判定された複数の判定期間のうちの最初の判定期間の始まりを不在開始時刻とすると共に最後の判定期間の終わりを不在終了時刻として特定することである。したがって、時間軸方向に移動する判定期間毎に在宅不在を判定しながら不在開始時刻と不在終了時刻とを特定し得るものであれば、上述の実施形態におけるＳ４乃至Ｓ１１の処理に対応する具体的な方法・手順は上述の実施形態として説明したもの（図１０を含む）に限定されるものではない。
なお更に付け加えると、所定の時間長さを有すると共に時間軸方向に移動する判定期間を用いることによって不在開始時刻と不在終了時刻とを特定することの利点は、判定対象時刻毎に各々独立して（即ち、当該の判定対象時刻よりも前の時刻における状況を合わせて考慮すること無く）在宅不在を判定する場合のような、居住者が在宅しているにもかかわらず、短時間、使用する電気機器が少なかったり電気機器の使用態様に変化が無かったりする状況での（言い換えると、電力需要速度の変化が小さい状況での）不在であるとの誤った判断を避けることである。この考え方は、上述の実施形態のように値の変化が散発的（言い換えると、間欠的）であり且つ離散的であるデータ、或いは言い換えると、値の変化の生起が時系列において散発的（間欠的）であり且つ値の変化の態様が離散的であるデータに基づいて作成される連続データが用いられる場合のみに適用され得るものではなく、種々の特性のデータが用いられる場合に適用され得るものであり、さらに言えば、判定対象時刻別の在宅不在の判定の仕方は上述の実施形態における方法に限定されるものではない。
例えば、もとより連続であるデータ、言い換えると、値の変化の生起が時系列において連続であり且つ値の変化の態様が漸次的であるデータ、具体的には例えば、あくまで一例として挙げると、１分毎に出力される１ Ｗｈ 単位のデータが用いられる場合にも適用され得る。この場合には、上述の実施形態のように最適化問題を解くこと無く、あくまで一例として挙げると３０分間隔の時刻の識別子ｔ毎の電力需要速度ｙt［Ｗｈ／３０分］や１分間隔の時刻の識別子ｔ毎の電力需要速度ｙt［Ｗｈ／分］、或いは一般化するとｎ分間隔の時刻の識別子ｔ毎の電力需要速度ｙt［Ｗｈ／ｎ分］が得られる。そして、上述の実施形態におけるＳ３以降の処理が同様に行われ得る。
あるいは、値の変化の生起が時系列において散発的（間欠的）であり且つ値の変化の態様が漸次的であるデータ、具体的には例えば、あくまで一例として挙げると、３０分毎に出力される１ Ｗｈ 単位のデータが用いられる場合にも適用され得る。この場合には、上述の実施形態のように最適化問題を解くこと無く、あくまで一例として挙げると、３０分間隔の時刻の識別子ｔ毎の電力需要速度ｙt［Ｗｈ／３０分］が得られる。そして、上述の実施形態におけるＳ３以降の処理が同様に行われ得る。
１０ ガス消費量の用途分解装置
１７ ガス消費量の用途分解プログラム
３０ 居住者の在宅不在の判定装置
３７ 居住者の在宅不在の判定プログラム
本発明は、居住者の不在期間の確定方法、装置、及びプログラムに関する。
以下の説明において、単位であることを明確にするために［ ］を用いる場合がある。
家庭における最適なエネルギーの使い方の探索や熱源機器の選択などを行うためには、電気やガスといったエネルギーの種類毎に用途毎の最終需要を把握することが必要である。しかしながら、ガスについては、配管切断を伴う流量計の設置は需要家の負担が大きく手間もコストもかかるため、個別機器の実測調査はほとんど行われていない。また、ガスメータの計測値から用途別（具体的には例えば、「給湯」，「厨房（コンロ）」，及び「暖房（ガスヒータ等）」の別）の消費量に分解する良好な手法は未だ無い。
ガスの配管を切断することなくガス流量を計測する装置として、例えば、ガス配管に設けられている膜式ガスメータに装着され、当該膜式ガスメータ内部の、ガスの流動に応じて周期運動する永久磁石の移動に伴う周期的な磁気変化を磁気センサによって検出し、当該周期的な磁気変化の検出に基づいて所定のガス流量（具体的には、ガスメータの単位計量体積）に対応したパルス出力を流量カウント数として所定時間毎に出力するものがある（例えば、特許文献１，特許文献２）。
上述のようなガス流量計測装置では、家庭におけるガスの消費量が所定時間毎にカウント数として把握され、所定時間は例えば１０秒などであり、１カウントに対応する実際のガス流量は例えば１．２［Ｌ］や１．０［Ｌ］などである。言い換えると、上述のようなガス流量計測装置では、所定の時間毎に、所定の流量毎に１カウントとして計数するカウント数として、ガスの消費量が出力される。
したがって、上述のようなガス流量計測装置では、具体的には例えば、所定時間が１０秒であると共に１カウント当たりの実際の流量が１．２［Ｌ］である場合には、１０秒間隔での、１．２Ｌ刻みのガス消費量が把握される。
また、ガス使用量を用途別に分解する従来の方法としては、例えば、世帯人数，厨房用ガス器具の所有台数，及び住居形態などのガス使用者の属性情報を用い、ガス使用者の世帯状況によってその使用量が異なる厨房用ガス使用量を求め、続いて、暖房を使用しない時期の総ガス使用量が厨房用と給湯用との和であることを利用してこの時期の給湯用ガス使用量を求め、次に、暖房を使用しない時期の給湯用ガス使用量と気温との相関関係を利用して暖房を使用する時期の給湯用ガス使用量を求め、さらに、総ガス使用量から厨房用及び給湯用ガス使用量を減算して暖房用ガス使用量を求めるものがある（特許文献３）。
また、電気に関しては、電気の使用（即ち、電力需要）から生活状態を推定することによって独居高齢者や家族を見守るサービスの実現が検討されたり、住居内の例えば空調等の機器の制御が自動で行われるシステムについて運用効率を改善して節電に繋げるために居住者が在宅であるか不在であるかの判定が行われたりしている例がある。具体的には例えばドアの開閉，人感センサ，或いはＣＯ2センサなどを利用するものがあり、また、住宅内への機器の設置が必要とされない非侵入型モニタリングの手法の一つとして電力需要データを用いるものが開発されている（非特許文献１）。
このような非侵入型モニタリングに関連しては、我が国において、スマートメータと呼ばれる、通信機能を備えて電力の利用状況（即ち、電力需要）をリアルタイムに確認することができる次世代電力計の低圧部門（具体的には、一般家庭）への設置が進められている（非特許文献２）。
特開２００４−１２５６３４号公報 特開２０１４−００６０５８号公報 特開２０００−２２１０６４号公報
Ｔ．Ａ．Ｎｇｕｙｅｎ ａｎｄ Ｍ．Ａｉｅｌｌｏ「Ｅｎｅｒｇｙ ｉｎｔｅｌｌｉｇｅｎｔ ｂｕｉｌｄｉｎｇｓ ｂａｓｅｄ ｏｎ ｕｓｅｒ ａｃｔｉｖｉｔｙ：Ａ ｓｕｒｖｅｙ」，Ｅｎｅｒｇｙ ａｎｄ Ｂｕｉｌｄｉｎｇｓ，Ｖｏｌ．５６，ｐｐ．２４４−２５７，２０１３年 経済産業省スマートメータ制度検討会「スマートメータの導入促進に伴う課題と対応について」，［online］，［平成２８年４月１４日アクセス］，インターネット＜ＵＲＬ：http://www.meti.go.jp/committee/summary/0004668/pdf/015_03_00.pdf＞
特許文献１や特許文献２のガス流量計測装置によって得られる計測データは、所定の時間毎に出力される、所定の流量毎に１カウントとして計数する流量カウント数である。しかしながら、例えばガスファンヒータやガスコンロのガス消費量は１．２［Ｌ／１０秒］に比べて通常は小さいため、１０秒間隔で出力される、実際の流量１．２［Ｌ］を１カウントとして計数する流量カウント数としての計測ではカウント値の変化は散発的（言い換えると、間欠的）になる。したがって、ガスファンヒータやガスコンロの各器具単体でのガス消費量に対応する計測値としての流量カウント数は、これら器具がたとえ継続的に使用されていても０若しくは１の出力となる。このため、時刻毎のガスの配管における実際のガス流量を適切に把握することは困難であり、延いては時刻毎の実際のガス流速を適切に把握することが困難であるという問題がある。
また、特許文献３のガス使用量の用途別分解方法では世帯人数，厨房用ガス器具の所有台数，及び住居形態などのガス使用者の属性情報が必要とされるところ、このような属性情報の入手は困難であり、仮に可能であったとしても多大な手間がかかってしまう。このため、汎用性が高いとは言い難い。加えて、特許文献３のガス使用量の用途別分解方法では気温の情報が更に必要とされるところ、対象世帯が所在する地域毎の、用途分解を行う期間の日毎の気温データの収集は多大な手間がかかるためにこの点においても汎用性が高いとは言い難く、さらに、気温データの地域区分の細かさによっては対象世帯の所在地における実際の気温との乖離が大きくなってしまうことも考えられるために良好な精度が常に確保されるとは言えず信頼性が高いとは言い難い。
また、次世代電力計として設置が進められているスマートメータで自動検針された電力需要値を電力会社（具体的には、送配電部門）へと送るルートは「Ａルート」と呼ばれる。このＡルートから得られる電力需要データ（「Ａルート計測値」と呼ぶ）は、３０分間隔で出力され、且つ、１００ Ｗｈ 刻み（言い換えると、１００ Ｗｈ 単位）の数値である。したがって、Ａルート計測値では、３０分間における実際の電力需要として１００ Ｗｈ に満たない電力需要は、当該時刻（３０分時間帯）の計測値としては切り捨てられ（言い換えると、計上されず）、次の時刻（３０分時間帯）の計測値に繰り越される。
実際の電力需要に対応して出力されるＡルート計測値は、具体的には例えば以下の表１のようになる。なお、表１における「実際の電力需要」は、電力需要の瞬時値が仮に計測されているとした場合の仮想の値であり、３０分間積分値である。
表１に示す例から、３０分間の電力需要が１００ Ｗｈ に満たない時刻区分（３０分時間帯）では実際には電力需要が有るにも拘わらず、Ａルート計測値では電力需要が無いことになったり、或いは、連続する時刻区分（具体的には、１１：００と１１：３０）において実際の電力需要としての差違は小さい（具体的には、約１６ Ｗｈ）にも拘わらず、Ａルート計測値では１００ Ｗｈ もの差違になったりすることが確認される。
このため、Ａルート計測値をそのまま利用して居住者が在宅であるか不在であるかを判定した場合には、適切な判定が行われない虞があり、高い精度で判定を行うことは困難である。特に、３０分間における実際の電力需要が１００ Ｗｈ に満たない時間帯においては、Ａルート計測値としては０ Ｗｈ と１００ Ｗｈ とが交互に出力される状態になり、実際の電力需要に比してその差違が相対的に大きな割合となって誤判定が多くなる。この点に関連しては、省エネルギー機器への切り替えなどにより、エアコンや電子レンジなどの電力消費量が比較的大きい機器を利用しない時間帯における電力需要は減少する傾向にあり、前記誤判定は益々多くなると考えられる。
そこで、本発明は、追加的な特別の機器の設置及び関連する情報やデータが必要とされることなく、例えばスマートメータのＡルート計測値のような、対象住宅や対象住戸の電力需要量の計測データであって値の変化の生起について所定の間隔で間欠的に把握され且つ値の変化の態様について所定の量を単位として離散的に把握された計測データ群に基づいて住居や住戸内に人が存在しているか存在していないかを良好な精度で推定することができる居住者の不在期間の確定方法、装置、及びプログラムを提供することを目的とする。
また、短い時間や特定の時刻毎に各々独立して（言い換えると、或る短い時間や特定の時刻のみを取り出して当該の或る短い時間や特定の時刻よりも前の状況を合わせて考慮すること無く）在宅不在を判定すると、居住者が在宅していても、短時間、使用する電気機器が少なかったり電気機器の使用態様に変化が無かったりする状況に対し、不在であると誤った判断をしてしまうことがあるという問題がある。
そこで、本発明は、短時間のみの状況に影響されること無く住居や住戸内に人が存在しているか存在していないかを推定することができる不在期間の確定方法を提供することを目的とする。
本発明の居住者の不在期間の確定方法は、電力需要量の計測データに基づいて時刻間毎の電力需要速度を算定し、判定対象時刻以前の所定の期間に含まれる前記電力需要速度を用いて前記所定の期間における前記電力需要速度の統計量を特徴量として算出し、少なくとも前記特徴量を用いて前記判定対象時刻別の居住者の在宅不在を判定し、所定の時間長さを有して時間軸方向に移動する判定期間に対応する前記判定対象時刻別の判定結果が全て不在であるときに当該判定期間中は不在であると判定し、不在と判定された一つの前記判定期間の始まりの時刻を不在開始時刻とすると共に終わりの時刻を不在終了時刻として特定し、または、連続して不在と判定された複数の前記判定期間のうちの最初の前記判定期間の始まりの時刻を不在開始時刻とすると共に最後の前記判定期間の終わりの時刻を不在終了時刻として特定するようにしている。
したがって、この不在期間の確定方法によると、判定対象時刻毎に各々独立して（即ち、当該の判定対象時刻よりも前の時刻における状況を合わせて考慮すること無く）在宅不在を判定する場合のような、居住者が在宅しているにもかかわらず、短時間、使用する電気機器が少なかったり電気機器の使用態様に変化が無かったりする状況での（言い換えると、電力需要速度の変化が小さい状況での）不在であるとの誤った判断をしてしまうことが防止される。
また、本発明の居住者の不在期間の確定方法は、前記電力需要速度が、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量を計測単位として計量する電力需要量との組み合わせデータとして取得された電力需要量の計測データから、所定の時間帯の前記組み合わせデータが分析対象区間のデータとして抽出され、前記分析対象区間における前記計測時刻毎の前記電力需要量の累積値を下限とすると共に前記計測時刻毎の前記電力需要量の累積値に前記所定の電力量を加えた値を上限とする領域内に存在し且つ前記分析対象区間の始点から終点へと至る総長が最短であり尚且つ各前記計測時刻において計測された前記電力需要量の総和と推定される各前記計測時刻の電力需要量の総和とが一致する折れ線が推定され、当該折れ線における各前記計測時刻間の直線の傾きとして算定されるようにしている。
また、本発明の居住者の不在期間の確定方法は、前記電力需要量の計測データがスマートメータのＡルートから出力されるデータである。
また、本発明の居住者の不在期間の確定装置は、電力需要量の計測データに基づいて時刻間毎の電力需要速度を算定する手段と、判定対象時刻以前の所定の期間に含まれる前記電力需要速度が用いられて前記所定の期間における前記電力需要速度の統計量が特徴量として算出され、少なくとも前記特徴量が用いられて前記判定対象時刻別の居住者の在宅不在を判定する手段と、所定の時間長さを有して時間軸方向に移動する判定期間に対応する前記判定対象時刻別の判定結果が全て不在であるときに当該判定期間中は不在であると判定され、不在と判定された一つの前記判定期間の始まりの時刻を不在開始時刻とすると共に終わりの時刻を不在終了時刻として特定する、または、連続して不在と判定された複数の前記判定期間のうちの最初の前記判定期間の始まりの時刻を不在開始時刻とすると共に最後の前記判定期間の終わりの時刻を不在終了時刻として特定する不在期間確定手段とを備えるようにしている。
また、本発明の居住者の不在期間の確定装置は、前記電力需要速度を算定する手段が、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量を計測単位として計量する電力需要量との組み合わせデータとして取得された電力需要量の計測データから、所定の時間帯の前記組み合わせデータを分析対象区間のデータとして抽出し、前記分析対象区間における前記計測時刻毎の前記電力需要量の累積値を下限とすると共に前記計測時刻毎の前記電力需要量の累積値に前記所定の電力量を加えた値を上限とする領域内に存在し且つ前記分析対象区間の始点から終点へと至る総長が最短であり尚且つ各前記計測時刻において計測された前記電力需要量の総和と推定される各前記計測時刻の電力需要量の総和とが一致する折れ線を推定して当該折れ線における各前記計測時刻間の直線の傾きとして電力需要速度を算定するものである。
また、本発明の居住者の不在期間の確定装置は、前記電力需要量の計測データがスマートメータのＡルートから出力されるデータである。
また、本発明の居住者の不在期間の確定プログラムは、電力需要量の計測データに基づいて時刻間毎の電力需要速度を算定する処理と、判定対象時刻以前の所定の期間に含まれる前記電力需要速度が用いられて前記所定の期間における前記電力需要速度の統計量が特徴量として算出され、少なくとも前記特徴量が用いられて前記判定対象時刻別の居住者の在宅不在を判定する処理と、所定の時間長さを有して時間軸方向に移動する判定期間に対応する前記判定対象時刻別の判定結果が全て不在であるときに当該判定期間中は不在であると判定され、不在と判定された一つの前記判定期間の始まりの時刻を不在開始時刻とすると共に終わりの時刻を不在終了時刻として特定する、または、連続して不在と判定された複数の前記判定期間のうちの最初の前記判定期間の始まりの時刻を不在開始時刻とすると共に最後の前記判定期間の終わりの時刻を不在終了時刻として特定する処理とをコンピュータに行わせるようにしている。
また、本発明の居住者の不在期間の確定プログラムは、前記電力需要速度を算定する処理が、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量を計測単位として計量する電力需要量との組み合わせデータとして取得された電力需要量の計測データから、所定の時間帯の前記組み合わせデータを分析対象区間のデータとして抽出し、前記分析対象区間における前記計測時刻毎の前記電力需要量の累積値を下限とすると共に前記計測時刻毎の前記電力需要量の累積値に前記所定の電力量を加えた値を上限とする領域内に存在し且つ前記分析対象区間の始点から終点へと至る総長が最短であり尚且つ各前記計測時刻において計測された前記電力需要量の総和と推定される各前記計測時刻の電力需要量の総和とが一致する折れ線を推定して当該折れ線における各前記計測時刻間の直線の傾きとして電力需要速度を算定するものである。
また、本発明の居住者の不在期間の確定プログラムは、前記電力需要量の計測データがスマートメータのＡルートから出力されるデータである。
本発明の居住者の不在期間の確定方法、装置、及びプログラムによれば、短時間、使用する電気機器が少なかったり電気機器の使用態様に変化が無かったりする状況での（言い換えると、電力需要速度の変化が小さい状況での）不在であるとの誤った判断をしてしまうことを防止することができ、不在期間の確定手法としての信頼性の向上を図ることが可能になる。
ガス消費量の用途分解方法の実施形態の一例を示すフローチャートである。 実施形態のガス消費量の用途分解方法をガス消費量の用途分解プログラムを用いて実施する場合の当該プログラムによって実現されるガス消費量の用途分解装置の機能ブロック図である。 実施形態の連続データの作成方法を説明する図であり、計測データとしての流量カウント数に基づいて流速を推定する方法を説明する図である。 ガス消費量の用途分解において用いられるアルゴリズムの一例を説明する用途決定チャートである。 本発明の居住者の不在期間の確定方法の一実施形態において用いる居住者の在宅不在の判定方法の実施形態の一例を示すフローチャート（Ｓ１乃至Ｓ３）である。 本発明の居住者の不在期間の確定方法の一実施形態を示すフローチャート（Ｓ４乃至Ｓ１１）である。 同実施形態の居住者の不在期間の確定プログラムによって実現される居住者の不在期間の確定装置の機能ブロック図である。 同実施形態の連続データの作成方法を説明する図であり、計測データとしての電力需要量に基づいて電力需要速度を推定する方法を説明する図である。 真の累積需要曲線としての最短折れ線を推定する方法を説明する図である。 本発明の居住者の不在期間の確定方法の実施形態の一例を説明する図であり、不在開始時刻及び不在終了時刻を特定する方法並びに不在期間を確定する方法の一例を説明する。
以下、本発明の構成を図面に示す実施の形態の一例に基づいて詳細に説明する。
図１から図４に、ガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムの実施形態の一例を示す。
本実施形態のガス消費量の用途分解方法は、ガス配管の主幹における計測（Ｓ１）によって所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量ｆu毎に１カウントとして計数する流量カウント数との組み合わせデータとして取得（Ｓ１）されたガス総消費量の計測データから、主幹から枝分かれした屋内配管に接続されているガス器具が使用されている時間帯の組み合わせデータが連続動作区間のデータとして抽出され（Ｓ２）、連続動作区間における計測時刻毎の流量カウント数の累積値と流量ｆuとを掛け合わせた値を下限とすると共に前記掛け合わせた値に流量ｆuを加えた値を上限とする領域内に存在し且つ連続動作区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測されたガス総消費量の総和と推定される各計測時刻のガス総消費量の総和とが一致する折れ線が推定され、当該折れ線における各計測時刻間の直線の傾きとして主幹におけるガスの流速ｙtが算定された上で（Ｓ３）、連続動作区間のデータが冬季であるか否か判断される処理（Ｓ４−１）と、ガスの流速ｙtが所定の範囲である状態が所定の時間に亙って継続しているか否かに基づいて暖房連続動作区間が特定される処理（Ｓ４−２）と、連続動作区間のデータにおいて特定された全ての暖房連続動作区間の中でのガスの流速ｙtの最小値が暖房ベース流速として設定される処理（Ｓ４−３）と、暖房ベース流速と給湯最小流速との和とガスの流速ｙtとの大きさが比較される処理（Ｓ４−４）と、給湯最小流速とガスの流速ｙtとの大きさが比較される処理（Ｓ４−５）と、連続動作区間のデータが夏季であるか否か判断される処理（Ｓ４−６）と、厨房最小流速と間欠暖房最大流速とのうちの少なくとも一方とガスの流速ｙtとの大きさが比較される処理（Ｓ４−７）とが実行されることにより、ガス総消費量の用途が給湯分，暖房分，及び厨房分に分類される（Ｓ４）ようにしている（図１参照）。
また、本実施形態では、Ｓ３の処理において、連続動作区間における計測時刻の時系列順に付与される識別子ｔ（ただし、ｔ＝１，２，３，…，Ｔ）毎の当該計測時刻の識別子ｔに対応する計測時刻における流量カウント数と流量ｆuとを掛け合わせた値ｘtが算出され、計測時刻の識別子ｔ毎の値ｘtが用いられて数式１乃至数式３に示す最適化問題のうちのいずれかが解かれることによって計測時刻の識別子ｔ毎の値ｙtとして主幹におけるガスの流速ｙtが算定されるようにしている。
上記ガス消費量の用途分解方法は、ガス消費量の用途分解装置によって実施され得る。本実施形態のガス消費量の用途分解装置１０は、ガス配管の主幹における計測によって所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量ｆu毎に１カウントとして計数する流量カウント数との組み合わせデータとして取得されたガス総消費量の計測データの入力を受ける手段としてのデータ受部１１ａと、計測データから、主幹から枝分かれした屋内配管に接続されているガス器具が使用されている時間帯の組み合わせデータを連続動作区間のデータとして抽出する手段としての分解区間抽出部１１ｂと、連続動作区間における計測時刻毎の流量カウント数の累積値と流量ｆuとを掛け合わせた値を下限とすると共に前記掛け合わせた値に流量ｆuを加えた値を上限とする領域内に存在し且つ連続動作区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測されたガス総消費量の総和と推定される各計測時刻のガス総消費量の総和とが一致する折れ線を推定して当該折れ線における各計測時刻間の直線の傾きとして主幹におけるガスの流速ｙtを算定する手段としての最短折線作成部１１ｃと、連続動作区間のデータが冬季であるか否かを判断し、ガスの流速ｙtが所定の範囲である状態が所定の時間に亙って継続しているか否かに基づいて暖房連続動作区間を特定すると共に連続動作区間のデータにおいて特定された全ての暖房連続動作区間の中でのガスの流速ｙtの最小値を暖房ベース流速として設定し、暖房ベース流速と給湯最小流速との和とガスの流速ｙtとの大きさを比較し、給湯最小流速とガスの流速ｙtとの大きさを比較し、連続動作区間のデータが夏季であるか否かを判断し、さらに、厨房最小流速と間欠暖房最大流速とのうちの少なくとも一方とガスの流速ｙtとの大きさを比較することにより、ガス総消費量の用途を給湯分，暖房分，及び厨房分に分類する手段としての用途分解部１１ｄとを有する。
また、本実施形態では、最短折線作成部１１ｃが、連続動作区間における計測時刻の時系列順に付与される識別子ｔ（ただし、ｔ＝１，２，３，…，Ｔ）毎の当該計測時刻の識別子ｔに対応する計測時刻における流量カウント数と流量ｆuとを掛け合わせた値ｘtを算出すると共に、計測時刻の識別子ｔ毎の値ｘtを用いて数式１乃至数式３に示す最適化問題のうちのいずれかを解くことによって計測時刻の識別子ｔ毎の値ｙtとして主幹におけるガスの流速ｙtを算定するようにしている。
また、上記ガス消費量の用途分解方法及びガス消費量の用途分解装置は、ガス消費量の用途分解プログラムがコンピュータ上で実行されることによっても実施・実現され得る。ここでは、ガス消費量の用途分解プログラムがコンピュータ上で実行されることによってガス消費量の用途分解装置が実現されると共にガス消費量の用途分解方法が実施される場合を説明する。
本実施形態のガス消費量の用途分解プログラム１７を実行するためのコンピュータ１０（本実施形態では、ガス消費量の用途分解装置１０でもある）の全体構成を図２に示す。このコンピュータ１０（ガス消費量の用途分解装置１０）は制御部１１，記憶部１２，入力部１３，表示部１４，及びメモリ１５を備え、これらが相互にバス等の信号回線によって接続されている。
制御部１１は、記憶部１２に記憶されているガス消費量の用途分解プログラム１７によってコンピュータ１０全体の制御並びにガス総消費量の用途の分解に係る演算を行うものであり、例えばＣＰＵ（中央演算処理装置）である。
記憶部１２は、少なくともデータやプログラムを記憶可能な装置であり、例えばハードディスクである。
入力部１３は、少なくとも作業者の命令や種々の情報を制御部１１に与えるためのインターフェイス（即ち、情報入力の仕組み）であり、例えばキーボードやマウスである。なお、例えばキーボードとマウスとの両方のように複数種類のインターフェイスを入力部１３として有するようにしても良い。
表示部１４は、制御部１１の制御によって文字や図形或いは画像等の描画・表示を行うものであり、例えばディスプレイである。
メモリ１５は、制御部１１が種々の制御や演算を実行する際の作業領域であるメモリ空間となるものであり、例えばＲＡＭ(Ｒandom Ａccess Ｍemory の略)である。
そして、コンピュータ１０（「ガス消費量の用途分解装置１０」と表記する）の制御部１１には、ガス消費量の用途分解プログラム１７が実行されることにより、ガス配管の主幹における計測によって所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量ｆu毎に１カウントとして計数する流量カウント数との組み合わせデータとして取得されたガス総消費量の計測データの入力を受ける処理を行うデータ受部１１ａと、計測データから、主幹から枝分かれした屋内配管に接続されているガス器具が使用されている時間帯の組み合わせデータを連続動作区間のデータとして抽出する処理を行う分解区間抽出部１１ｂと、連続動作区間における計測時刻毎の流量カウント数の累積値と流量ｆuとを掛け合わせた値を下限とすると共に前記掛け合わせた値に流量ｆuを加えた値を上限とする領域内に存在し且つ連続動作区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測されたガス総消費量の総和と推定される各計測時刻のガス総消費量の総和とが一致する折れ線を推定して当該折れ線における各計測時刻間の直線の傾きとして主幹におけるガスの流速ｙtを算定する処理を行う最短折線作成部１１ｃと、連続動作区間のデータが冬季であるか否かを判断し、ガスの流速ｙtが所定の範囲である状態が所定の時間に亙って継続しているか否かに基づいて暖房連続動作区間を特定すると共に連続動作区間のデータにおいて特定された全ての暖房連続動作区間の中でのガスの流速ｙtの最小値を暖房ベース流速として設定し、暖房ベース流速と給湯最小流速との和とガスの流速ｙtとの大きさを比較し、給湯最小流速とガスの流速ｙtとの大きさを比較し、連続動作区間のデータが夏季であるか否かを判断し、さらに、厨房最小流速と間欠暖房最大流速とのうちの少なくとも一方とガスの流速ｙtとの大きさを比較することにより、ガス総消費量の用途を給湯分，暖房分，及び厨房分に分類する処理を行う用途分解部１１ｄとが構成される。
また、本実施形態では、最短折線作成部１１ｃが、連続動作区間における計測時刻の時系列順に付与される識別子ｔ（ただし、ｔ＝１，２，３，…，Ｔ）毎の当該計測時刻の識別子ｔに対応する計測時刻における流量カウント数と流量ｆuとを掛け合わせた値ｘtを算出すると共に、計測時刻の識別子ｔ毎の値ｘtを用いて数式１乃至数式３に示す最適化問題のうちのいずれかを解くことによって計測時刻の識別子ｔ毎の値ｙtとして主幹におけるガスの流速ｙtを算定する処理を行うようにしている。
そして、ガス消費量の用途分解方法の実施として、まず、ガス総消費量の計測が行われて計測データの取得が行われる（Ｓ１）。
ガス総消費量の計測の仕方は、所定の流量毎に１カウントとして計数するカウント数として主幹のガス流量を計測し、所定の時間毎に例えば電気的な信号として前記カウント数を出力するものであれば、特定の方法や機器によるものに限定されるものではなく、前記のような計測を行い得る方法や機器が適宜選択される。
具体的には例えば、あくまで一例として挙げると、膜式ガスメータに装着される前述のガス流量計測装置（例えば、特開２００４−１２５６３４号公報，特開２０１４−６０５８号公報）が利用され得る。あるいは、一定量（例えば、ガスメータの単位計量体積）のガスがガスメータを通過するたびに動作サイクルに起因してガスメータから発生する周期的な動作音や振動を検出してその回数をカウント数として所定の時間毎に出力する機器が利用され得る。
ガス総消費量（主幹のガス流量）の計測における、上述の「所定の流量」としての、１カウント当たりの実際の流量は、例えば、主幹に設けられているガスメータの単位計量体積に従って決定される。具体的には例えば、ガスメータの単位計量体積が１［Ｌ］である場合にはガス流量の計測機器から出力される１カウント毎に１［Ｌ］のガスが実際に流動した（言い換えると、消費された）ことを意味し、或いは、ガスメータの単位計量体積が１．２［Ｌ］である場合にはガス流量の計測機器から出力される１カウント毎に１．２［Ｌ］のガスが実際に流動したことを意味する。
また、ガス総消費量の計測における、上述の「所定の時間」としての、カウント数を出力する時間間隔は、特定の時間［秒，分］に限定されるものではなく、計測に用いられる計測機器の仕様を踏まえた上で、例えばガス消費量に関して想定される変動ピッチや必要とされる推定精度などが考慮されて適当な時間に適宜設定される。具体的には例えば、あくまで一例として挙げると、５〜６０秒程度の範囲で適当な値に設定され得る。
本実施形態では、上述の「所定の流量」が１．２［Ｌ］であると共に上述の「所定の時間」が１０秒である計測データ、すなわち、実際の流量１．２［Ｌ］を１カウントとして計数するカウント数（「流量カウント数」と呼ぶ）が１０秒間隔で出力される計測データがＳ１の処理によって取得されるとする。
そして、計測によって取得されたガス総消費量（主幹のガス流量）の計測値（具体的には、流量カウント数）は、データ受部１１ａを介してガス消費量の用途分解装置１０に入力される。
ガス総消費量の計測値（流量カウント数）は、ガス消費量の用途分解装置１０に、記録媒体やデータサーバを介して入力されるようにしても良いし、ガス消費量の用途分解装置１０と計測機器とがデータや制御指令等の信号の送受信（即ち、出入力）が可能であるように電気的に接続されて入力されるようにしても良い。
記録媒体を介して計測値が入力される場合には、計測機器から出力された流量カウント数が記録媒体に記録され、当該記録媒体がガス消費量の用途分解装置１０の記録媒体接続用端子（図示していない）に差し込まれることによって流量カウント数が入力されるようにしても良い。
また、データサーバを介して計測値が入力される場合には、データサーバがバス等の信号回線によってガス消費量の用途分解装置１０に接続され、計測機器から出力された流量カウント数が前記データサーバにデータファイル等として格納(保存)され、当該データファイル等として保存された流量カウント数が読み込まれるようにしても良い。
また、ガス消費量の用途分解装置１０と計測機器とが電気的に接続されて計測値が入力される場合には、例えば、各々に接続されて敷設されたケーブル等が用いられる有線による信号送受の仕組みを介して信号の送受信が可能であるように電気的に接続されるようにしても良いし、各々に接続された無線信号送受信機が用いられる無線による信号送受の仕組みを介して信号の送受信が可能であるように電気的に接続されるようにしても良い。そして、これら信号送受の仕組みによって流量カウント数が入力されるようにしても良い。
さらに、記録媒体やデータサーバや有線・無線による信号送受の仕組みが組み合わされて用いられて流量カウント数がガス消費量の用途分解装置１０に入力されるようにしても良い。
また、ガス総消費量の各計測値（各流量カウント数）が計測機器から出力された時刻が計測機器に備えられた時計機能によって計測値と一緒に出力されて当該計測値と対応づけられて記録されたり、或いは、計測機器から出力された各計測値が記録媒体やデータサーバやガス消費量の用途分解装置１０に記録されたり入力されたりした時の時刻が当該計測値と対応づけられて記録されたりする。そして、この時刻の間隔［秒，分］が上述の「所定の時間」になる。
ここでの説明においては、各計測値（各流量カウント数）と対応づけられて記録される時刻のことを「計測時刻」と呼び、流量カウント数と計測時刻との組み合わせデータのことを「計測データ」と呼ぶ。
そして、データ受部１１ａにより、計測データがメモリ１５に記憶させられる。
なお、例えば、対象の家庭に関する（言い換えると、対象住宅や対象住戸の主幹における）、あくまで一例として挙げると２４時間に亙る一年分の計測データが、ガス消費量の用途分解のひとかたまりの処理対象としてメモリ１５に記憶させられる。
次に、制御部１１の分解区間抽出部１１ｂにより、Ｓ１の処理によって計測されて取得された計測データからの、用途の分解が行われる時間帯のデータの抽出が行われる（Ｓ２）。
Ｓ２の処理は、Ｓ１の処理においてメモリ１５に記憶された、ガス消費量の用途分解のひとかたまりの処理対象としての計測データの全体に亙って実行される。
ここで説明する手法はガス総消費量を用途別に分解するものであるので、対象住宅・住戸においてガス器具が使用されてガスが実際に消費されている時間帯が処理対象になり、ガス器具が使用されておらずガスが消費されていない時間帯は処理対象にならない。
そこで、計測データから、何らかのガス器具が使用される（言い換えると、動作する）ことによってガス総消費量が連続して変化（具体的には、増加）している可能性がある時間帯（「連続動作区間」と呼ぶ）の計測データが切り出される。
具体的には、動作休止判定時間［秒，分］に亙って計測データの流量カウント数が０(零)である場合に、ガス器具の使用が途切れたと判定される。動作休止判定時間以上に亙って計測データの流量カウント数が０である時間帯のことを「動作休止区間」と呼ぶ。
そして、時間的に前後する二つの動作休止区間に挟まれる時間帯が連続動作区間として切り出される。したがって、連続動作区間は、計測データの流量カウント数が動作休止判定時間以上に亙って０になることなく継続的に変化している時間帯の計測データであるとも言える。
動作休止判定時間［秒，分］は、特定の値に限られるものではなく、必要とされる推定精度などが考慮されて適当な値に適宜設定される。具体的には例えば、あくまで一例として挙げると、主幹におけるガスの流速が０．２［Ｌ／１０秒］以下である変化は微小な流量として無視しても大きな誤差は生じない（言い換えると、必要な精度は確保される）と判断される場合に、１．２Ｌ／０．２Ｌ＝６（即ち、計測時刻の間隔６個分）に基づいて、動作休止判定時間が１０秒×６＝６０秒に設定されることが考えられる。
なお、動作休止判定時間の値［秒，分］は、ガス消費量の用途分解プログラム１７内に予め規定されるようにしても良いし、Ｓ２の処理が行われる前に入力部１３を介して作業者によって入力されるようにしても良い。
本実施形態では、分解区間抽出部１１ｂにより、Ｓ１の処理においてメモリ１５に記憶された計測データが時系列に従って順次読み込まれ、動作休止判定時間に基づいて動作休止区間が特定されると共に、時間的に前後する二つの動作休止区間に挟まれる時間帯が連続動作区間として抽出される。
そして、分解区間抽出部１１ｂにより、抽出された連続動作区間毎の計測データが、以降のＳ３及びＳ４の処理対象の計測データとしてメモリ１５に記憶させられる。
次に、制御部１１の最短折線作成部１１ｃにより、Ｓ２の処理によって抽出された連続動作区間に関する計測データが用いられて累積需要曲線の作成が行われる（Ｓ３）。
Ｓ３の処理は、Ｓ２の処理において抽出された連続動作区間の各々について実行される。
Ｓ１の処理によって取得される計測データは、１０秒間隔での１．２Ｌ刻みのガス総消費量である。しかしながら、何らかのガス器具が使用されているときは、ガス総消費量が常に増加するという点においてガス総消費量の変化の生起は時系列において連続であり、また、ガスが常に流動するという点においてガス総消費量の変化の態様は漸次的である。
つまり、Ｓ１の処理における計測では、実際には(値の)変化の生起が時系列において連続であると共に(値の)変化の態様が漸次的である現象が、所定の時間間隔且つ所定の値刻みという散発的（言い換えると、間欠的）且つ離散的な変化として把握されていると言える。
そして、ここで説明する手法は、所定の時間間隔且つ所定の流量刻みという計測値の変化が散発的（言い換えると、間欠的）であり且つ離散的である計測データがそのまま用いられるのではなく、当該散発的（間欠的）且つ離散的な計測データが、計測値の変化の生起が時系列において連続であり且つ計測値の変化の態様が漸次的であるデータに変換された上で用途分解の処理に用いられる。
本実施形態において用いられる、計測値の変化の生起が時系列において散発的（間欠的）であり且つ計測値の変化の態様が離散的であるデータに基づいて、計測値の変化の生起が時系列において連続であり且つ計測値の変化の態様が漸次的である累積需要曲線が作成される、連続データの作成方法を、図３を用いて説明する。
なお、連続動作区間では、何らかのガス器具が使用されることによってガス総消費量が常に増加するという点においてガス総消費量の変化の生起が時系列において連続であり、また、何らかのガス器具が使用されることによってガスが常に流動するという点においてガス総消費量の変化の態様が漸次的である。
ここでの説明の例として、図３に示すような、或る連続動作区間に関する計測データ（具体的には、１０秒毎に出力される流量カウント数の累積値；図中の△印）がＳ２の処理によって抽出されたとする。なお、１カウント当たりの実際の流量は１．２［Ｌ］である。
図３の横軸は、１０秒間隔である計測時刻の各時点を、当該の連続動作区間の最初の計測時刻を１とした上で識別子としての連続番号で表したものである。すなわち、横軸における例えば「１０」は、当該の連続動作区間の始まりから１００秒（＝１０×１０秒）だけ時間が経過していることを表す。
また、図３の縦軸は、当該の連続動作区間おける流量カウント数の累積値を表す。すなわち、縦軸における例えば「４」は、当該の連続動作区間の始まりから４．８［Ｌ］（＝４×１．２［Ｌ］）だけガスが消費されていることを表す。
Ｓ１の処理によって取得される計測データは所定の時間間隔且つ所定の流量刻みという値の変化が散発的（間欠的）且つ離散的なデータであるので、すなわち、計測値の変化の生起が時系列において散発的（間欠的）且つ計測値の変化の態様が離散的なデータであるので、データプロットである△印を時系列で結ぶと図３に示すように階段状になる（図３中において符号２１；「下限階段直線２１」と呼ぶ）。
また、計測データの流量カウント数の累積値に計測時刻毎に１を加えたデータのプロットである▽印を時系列で結んだものも、図３に示すように階段状になる（図３中において符号２２；「上限階段直線２２」と呼ぶ）。
そして、下限階段直線２１が実際のガス総消費量に基づく真の累積値の下限（言い換えると、実際のガス総消費量が取り得る累積値の下限）であり、上限階段直線２２が実際のガス総消費量に基づく真の累積値の上限（言い換えると、実際のガス総消費量が取り得る累積値の上限）である。
したがって、実際のガス総消費量の時系列における累積を表す真の累積需要曲線は、下限階段直線２１と上限階段直線２２とに挟まれた領域（「存在可能領域」と呼ぶ）内に存在し、全体としては増加傾向を示す非減少の曲線であると考えられる。
また、家庭において一般的に使用されるガス器具の実際のガス消費の特性を本発明者が分析することにより、ガス器具によって消費されるガスの流量の大きな変化の回数は比較的少なく、ほぼ一定の流速が或る程度の時間に亙って継続する傾向があることが確認された。
このため、真の累積需要曲線としては、存在可能領域内において増加傾向を示す非減少の折れ線のうち、折れ点（言い換えると、変曲点、若しくは、流速変化点）の個数が比較的少ないものが妥当性が高いと考えられる。
そこで、実際のガス総消費量の時系列における累積を表す真の累積需要曲線として、存在可能領域内に存在する、総長（即ち、当該の連続動作区間の始まりから終わりまでに亙る軌跡の長さ）が最短の折れ線（「最短折れ線」と呼ぶ）を推定する。
真の累積需要曲線としての、存在可能領域内に存在する、最短折れ線の推定手法は、特定の方法に限定されるものではなく、特定の領域内に存在して所定の始点から終点へと至る折れ線のうちの総長が最短の折れ線であり且つ各計測時刻において計測されたガス総消費量の総和と推定される各計測時刻のガス総消費量の総和とが一致する折れ線を推定し得る適当な方法が適宜選択される。
最短折れ線の推定手法として、具体的には例えば、あくまで一例として挙げると、当該の連続動作区間において確定済みの折れ線の最後の点を起点とした（未確定の）直線部の許容される傾きの上限値・下限値を順次更新し、全ての直線が存在可能領域を逸脱せざるを得なくなる時点で直線の傾きを決定し、その直線が存在可能領域の境界に最初に接した時点までの折れ線の確定を行い、そして、新たに確定した折れ線の最後の点を新たな起点に設定した上で上述と同様の処理を繰り返すことを基本的な考え方とする方法が用いられ得る。
この推定方法は、計測時刻ｔ（尚、計測間隔は１０秒であり、したがって計測時刻の間隔は１０秒である）における計測流量をｘt［Ｌ／１０秒］とすると共に推定流速をｙt［Ｌ／１０秒］とすると、最適化問題として数式４のように表される。
数式４においてｔは連続動作区間における計測時刻の識別子であってｔ＝１，２，３，…であり、数式４では当該の連続動作区間の範囲として計測時刻の識別子ｔが１からＴまでの期間に亙って処理が実行される。
そして、計測時刻を１刻みの識別子とすることにより、計測の時間間隔である１０秒を一単位として扱うことになり、計測時刻 ｔ−１ から ｔ にかけての増分は計測の時間間隔当たりの増加量であって計測時刻 ｔ−１ から ｔ にかけての速度［Ｌ／１０秒］であり、真の累積需要曲線としては計測時刻 ｔ−１ から ｔ にかけての直線の傾きである。
数式４において、また、計測流量ｘt［Ｌ／１０秒］は、計測時刻ｔにおける流量カウント数（言い換えると、計測時刻 ｔ−１ から ｔ にかけての流量カウント数の増分）と１カウント当たりの実際の流量［Ｌ］とを掛け合わせたものである。
数式４において、さらに、ｆuは、１カウント当たりの実際の流量［Ｌ］である（即ち、本実施形態ではｆu＝１．２である）。なお、ｆuの値［Ｌ］は、ガス消費量の用途分解プログラム１７内に予め規定されるようにしても良いし、Ｓ３の処理が行われる前に入力部１３を介して作業者によって入力されるようにしても良い。
制約条件の第一式は流速の非負性に対応し、第二式は実際のガス総消費量の累積値の範囲制約に対応し、さらに、第三式は当該の連続動作区間における計測値ベースのガス総消費量の総和と推定値ベースのガス総消費量の総和との一致に対応する。
上述の最適化問題を解く方法によると、最悪の場合でも、時系列長（即ち、計測時刻の個数）の２乗のオーダーで真の累積需要曲線に相当する最短折れ線が推定される。
図３に示す計測データのプロット（図中の△印）を結んだ下限階段直線２１に対して数式４で表される最適化問題を適用した結果として、図３において点線で示す、実際のガス総消費量に近いと考えられる推定累積需要曲線２３が得られる。
ここで、上記の最適化問題は計測時刻の識別子ｔ＝１から処理が実行されるところ、ｔ＝１の処理が実行されるためには、ｔ＝１での起点としてｔ＝０（即ち、計測時刻の識別子ｔ＝１に対応する計測時刻の一つ前の計測時刻）からｔ＝１にかけての折れ線の最後の点が特定される必要が、即ち、ｔ＝０からｔ＝１にかけての直線が特定される必要がある。この直線は、具体的には、当該の連続動作区間の最初の計測時刻（即ち、計測時刻の識別子ｔ＝１）の直前の計測時刻として識別子ｔを０(零)とすると共に流量カウント数の累積値を０(零)とする仮想的な原点を置き、当該原点と計測時刻の識別子ｔ＝１における流量カウント数の累積値の値の点とを結ぶことによって特定される。
あるいは、ほぼ一定の流速が或る程度の時間に亙って継続する傾向があるという家庭におけるガス消費の特性に基づいて、真の累積需要曲線としては折れ点の個数が少ないものが妥当性が高いという考えのもと、当該の連続動作区間についての計測時刻の識別子ｔ＝０からｔ＝１にかけての直線の傾きが、識別子ｔ＝１において連続する直線の傾きに一致させられるようにしても良い。
ここで、真の累積需要曲線の上述の推定方法（即ち、数式４で表される最適化問題）は、目的関数部を「距離の和」から「距離の２乗和」に変更した数式５で表される最適化問題、即ち「存在可能領域の任意の曲線のうち各計測時刻の流速（即ち、流量の累積値の計測時刻間差分）の２乗和が最小の曲線」を求める問題と等価である。
数式５のように表される最適化問題では、２乗和の最小化により、各計測時刻の流速が同じ値を持つ傾向が生じるので、折れ線が最適解になる。
真の累積需要曲線の推定方法としての数式４で表される最適化問題は、また、数式５で表される最適化問題の双対問題である数式６で表される最適化問題とも等価である。なお、数式６で表される最適化問題は、全変動正規化問題と呼ばれる凸最適化問題の一種である。なお、数式６において、ｘ0は計測流量ｘt［Ｌ／１０秒］の初期値（即ち、ｔ＝０のときの値）であり、図３に示す例ではｘ0＝０である。
上述の真の累積需要曲線（即ち、最短折れ線）は、一例として以下に示すアルゴリズムが用いられることによって求められ得る。
−真の累積需要曲線の決定アルゴリズム−
《保持・更新情報》
○起点 折れ線
○上限直線 上限需要と領域境界接触点
○下限直線 下限需要と領域境界接触点
《更新内容》
(1) 時刻更新：時刻＝時刻＋１
(2) 下限直線＜上限累積需要値
起点から下限境界接触点までの需要量＝需要下限（確定）
起点＝下限直線の境界接触点
時刻＝起点，上限／下限直線の再設定
(1)へ
(3) 上限直線＞下限累積需要
起点から上限境界接触点までの需要量＝需要上限（確定）
起点＝上限直線の境界接触点
時刻＝起点，上限／下限直線の再設定
(1)へ
(4) 下限直線＜領域下限
下限需要，下限境界設定更新
上限直線＞領域上限
上限需要，上限境界設定更新
(5) 計測時刻の識別子ｔ＝Ｔになるまで(1)へ戻る
本実施形態では、最短折線作成部１１ｃにより、Ｓ２の処理においてメモリ１５に記憶された連続動作区間に関する計測データ（即ち、流量カウント数と計測時刻との組み合わせデータ）が読み込まれる。
そして、メモリ１５から読み込まれた流量カウント数は１カウント当たりの実際の流量である１．２［Ｌ］と掛け合わされて実際のガスの流量［Ｌ］に変換されると共に計測時刻は当該の連続動作区間における計測時刻の識別子として時系列順に１，２，３，…という連続番号に変換される。
続いて、最短折線作成部１１ｃにより、上記変換された計測データが用いられて数式４乃至数式６に表される最適化問題のうちのいずれかが解かれて、計測時刻の識別子ｔ毎のガスの推定流速ｙt［Ｌ／１０秒］が算定される。
そして、最短折線作成部１１ｃにより、Ｓ２の処理において抽出された連続動作区間毎に、算定されたガスの推定流速ｙt［Ｌ／１０秒］の値が計測時刻の識別子ｔの値と対応づけられて（言い換えると、計測時刻の識別子ｔの値とガスの推定流速ｙtの値との組み合わせデータとして）メモリ１５に記憶させられる。
次に、制御部１１の用途分解部１１ｄにより、Ｓ３の処理によって算定された計測時刻別のガスの流速が用いられてガス総消費量の用途毎のガス使用量への分解が行われる（Ｓ４）。
Ｓ４の処理は、Ｓ２の処理において抽出されてＳ３の処理において計測時刻別のガスの流速が算定された連続動作区間の各々について実行される。
具体的には、用途分解部１１ｄにより、Ｓ３の処理においてメモリ１５に記憶された連続動作区間に関する計測時刻の識別子ｔの値（ｔ＝１，２，…，Ｔ）とガスの推定流速ｙtの値との組み合わせデータが読み込まれる。
続いて、本実施形態では、図４に示すアルゴリズムに従ってガス総消費量が用途毎のガス使用量に分解される。
具体的には、まず、当該の連続動作区間が冬季におけるものであるか否かが判断される（Ｓ４−１）。
冬季であるか否かの判断の仕方は、特定の判断基準に限定されるものではなく、冬季とその他の季節とにおける例えば暖房器具の使用有無や給湯器具の温度設定などの差異が顕れる（若しくは、顕れると考えられる）時期などが考慮された上で適当な方法に適宜設定される。具体的には例えば、あくまで一例として挙げると、１０月から翌年４月までの期間に取得された計測データから抽出された連続動作区間のデータは冬季であると判断する方法が採用され得る。
そして、当該の連続動作区間が冬季におけるものである場合には（Ｓ４−１：Ｙｅｓ）Ｓ４−２の処理に進む。
そして、推定流速ｙtが暖房基準流速の最小値Ｖhs-min以上であると共に暖房基準流速の最大値Ｖhs-max以下であり、且つ、前記条件を満たす流速の継続時間（「所定流速継続時間」と呼ぶ）が暖房運転判断時間Ｔh以上であるか否かが判断される（Ｓ４−２）。
Ｓ４−２の判断処理は、流速は必ずしも大きくはないもののそのような流速の継続時間が長い場合には少なくとも暖房用ガス器具が使用されていると考え、暖房用ガス器具が連続して動作している時間帯（「暖房連続動作区間」と呼ぶ）を特定するための処理である。
暖房基準流速の最小値Ｖhs-min及び最大値Ｖhs-max［Ｌ／１０秒］は、暖房用ガス器具が使用されているときのガス消費量として想定される範囲の最小値及び最大値に基づいて設定され、特定の値に限定されるものではなく、例えば、対象住宅・住戸に設置されている（若しくは、設置されていると想定される一般的な）暖房用ガス器具の仕様などが考慮されたり、或いは、実測データが収集されている場合には当該実測データが考慮されたりして適当な値に適宜設定される。具体的には例えば、あくまで一例として挙げると、最小値Ｖhs-minが０．０１〜０．２［Ｌ／１０秒］程度の範囲で適当な値に、最大値Ｖhs-maxが０．５〜１．５［Ｌ／１０秒］程度の範囲で適当な値に、それぞれ設定され得る。
また、暖房運転判断時間Ｔh［秒，分］は、特定の値に限定されるものではなく、例えば暖房器具を使用する場合の最低継続時間として想定される時間などが考慮されて適当な値に適宜設定される。具体的には例えば、あくまで一例として挙げると、５〜１５分程度の範囲で適当な値に設定され得る。
Ｓ４−２の処理は、当該の連続動作区間の全体に亙って実行される。そして、Ｓ４−２の処理の結果として、各連続動作区間について、暖房連続動作区間が一つも特定されないか、または、暖房連続動作区間が一つ以上特定される。
そして、Ｖhs-min≦ｙt≦Ｖhs-max であると共に所定流速継続時間がＴh［分］以上である場合には（Ｓ４−２：Ｙｅｓ）Ｓ４−３の処理に進む。言い換えると、Ｓ４−２の処理において暖房連続動作区間が一つ以上特定された場合にＳ４−３の処理に進む。
そして、暖房ベース流速Ｖhbが設定される（Ｓ４−３）。
暖房ベース流速Ｖhb［Ｌ／１０秒］は、当該の連続動作区間において特定された全ての暖房連続動作区間の中での、推定流速ｙt［Ｌ／１０秒］の最小値に設定される。
したがって、暖房ベース流速Ｖhbは、或る一つの連続動作区間に対して共通に用いられる値として設定され、異なる連続動作区間では値が異なることがあり得る。
次に、推定流速ｙtが暖房ベース流速Ｖhbと給湯最小流速Ｖwとの和以上であるか否かが判断される（Ｓ４−４）。
Ｓ４−４の処理は、Ｓ４−２の処理において特定された暖房連続動作区間それぞれの各計測時刻ｔについて実行される。
Ｓ４−４の判断処理は、当該の連続動作区間において暖房連続動作区間であると特定された時間帯の各計測時刻ｔについて、即ち暖房用ガス器具が使用されている状態での、給湯用ガス器具の使用の有無を判断するための処理である。
給湯最小流速Ｖw［Ｌ／１０秒］は、給湯用ガス器具が使用されているときのガス消費量として想定される範囲の最小値に基づいて設定され、特定の値に限定されるものではなく、例えば、対象住宅・住戸に設置されている（若しくは、設置されていると想定される一般的な）給湯用ガス器具の仕様などが考慮されたり、或いは、実測データが収集されている場合には当該実測データが考慮されたりして適当な値に適宜設定される。具体的には例えば、あくまで一例として挙げると、０．５〜１．５［Ｌ／１０秒］程度の範囲で適当な値に設定され得る。
そして、当該の連続動作区間における暖房連続動作区間内の計測時刻ｔ毎に、ｙt≧Ｖhb＋Ｖwである場合には（Ｓ４−４：Ｙｅｓ）推定流速ｙtから暖房ベース流速Ｖhbを引いたものが給湯分のガス使用量に分類されると共に暖房ベース流速Ｖhbが暖房分のガス使用量に分類され、ｙt≧Ｖhb＋Ｖwでない場合には（Ｓ４−４：Ｎｏ）推定流速ｙtが暖房分のガス使用量に分類される。
一方、Ｓ４−１の処理において、当該の連続動作区間が冬季におけるものでない場合には（Ｓ４−１：Ｎｏ）Ｓ４−５の処理に進む。そして、計測時刻ｔのそれぞれについてＳ４−５乃至Ｓ４−７の処理が実行される。
また、Ｓ４−２の処理において、Ｖhs-min≦ｙt≦Ｖhs-max でない場合、若しくは、所定流速継続時間がＴh［分］以上でない場合にも（Ｓ４−２：Ｎｏ）Ｓ４−５の処理に進む。言い換えると、暖房連続動作区間ではない計測時刻ｔのそれぞれについてＳ４−５乃至Ｓ４−７の処理が実行される。
まず、推定流速ｙtが給湯最小流速Ｖw以上であるか否かが判断される（Ｓ４−５）。
Ｓ４−５の判断処理は、給湯用ガス器具の使用の有無を判断するための処理である。
そして、ｙt≧Ｖwである場合には（Ｓ４−５：Ｙｅｓ）推定流速ｙtが給湯分のガス使用量に分類される。
一方、ｙt≧Ｖwでない場合には（Ｓ４−５：Ｎｏ）Ｓ４−６の処理に進み、当該の連続動作区間が夏季におけるものであるか否かが判断される。
夏季であるか否かの判断（Ｓ４−６）の仕方は、特定の判断基準に限定されるものではなく、夏季とその他の季節とにおける例えば暖房器具の使用有無や給湯器具の温度設定などの差異が顕れる（若しくは、顕れると考えられる）時期などが考慮された上で適当な方法に適宜設定される。具体的には例えば、あくまで一例として挙げると、１１月から翌年３月までの期間に取得された計測データから抽出された連続動作区間のデータは夏季であると判断する方法が採用され得る。
なお、冬季や夏季であるか否かを期間によって判断するようにした場合に、冬季と夏季として指定される期間は、切れ目無く連続しているようにしても良いし、冬季にも夏季にも指定（分類）されない期間があっても良く、また、冬季として指定される期間と夏季として指定される期間とで一部が重複しているようにしても良い。
そして、当該の連続動作区間が夏季におけるものである場合には（Ｓ４−６：Ｙｅｓ）推定流速ｙtが厨房分のガス使用量に分類される。
一方、当該の連続動作区間が夏季におけるものでない場合には（Ｓ４−６：Ｎｏ）Ｓ４−７の処理に進み、推定流速ｙtが厨房最小流速Ｖkよりも大きいか否かが判断される。
Ｓ４−７の判断処理は、厨房用ガス器具の使用の有無を判断するための処理である。
厨房最小流速Ｖk［Ｌ／１０秒］は、厨房用ガス器具が使用されているときのガス消費量として想定される範囲の最小値に基づいて設定され、特定の値に限定されるものではなく、例えば、対象住宅・住戸に設置されている（若しくは、設置されていると想定される一般的な）厨房用ガス器具の仕様などが考慮されたり、或いは、実測データが収集されている場合には当該実測データが考慮されたりして適当な値に適宜設定される。具体的には例えば、あくまで一例として挙げると、０．２〜１．０［Ｌ／１０秒］程度の範囲で適当な値に設定され得る。
ここで、本実施形態のＳ４−７の処理では、厨房用ガス器具のみが使用されている状態と暖房用ガス器具が保温のために間欠的に動作している状態とが区別されるようにしている。
そこで、厨房最小流速Ｖkの値は、暖房用ガス器具が保温のために間欠的に動作しているときのガス消費量として想定される範囲の最大値（「間欠暖房最大流速」と呼ぶ）が合わせて考慮されたり当該間欠暖房最大流速に基づいて設定されたりするようにしても良い。
そして、推定流速ｙtが厨房最小流速Ｖkよりも大きい場合には（Ｓ４−７：Ｙｅｓ）推定流速ｙtが厨房分のガス使用量に分類される。
一方、推定流速ｙtが厨房最小流速Ｖkよりも大きくない場合には（Ｓ４−７：Ｎｏ）推定流速ｙtが暖房分のガス使用量に分類される。
そして、用途分解部１１ｄにより、ガス器具毎のガス使用量の例えば３０分間や１時間など時間，日，週間，月，季節，及び／又は年などの単位での積算値［Ｌ］が算出され、データファイル等に記録されて記憶部１２などに格納(保存)されたり表示部１４に表示されたりする。
なお、Ｓ４の処理において用いられる各閾値などの種々の値は、ガス消費量の用途分解プログラム１７内に予め規定されるようにしても良いし、Ｓ４の処理が行われる前に入力部１３を介して作業者によって入力されるようにしても良い。
そして、制御部１１は、Ｓ１の処理においてメモリ１５に記憶された、ガス消費量の用途分解のひとかたまりの処理対象としての計測データに関する処理を終了する（ＥＮＤ）。
以上のように構成されたガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムによれば、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量ｆu［Ｌ］毎に１カウントとして計数する流量カウント数との組み合わせデータ、すなわち、所定の時間間隔且つ所定の流量ｆu刻みという値の変化が散発的（言い換えると、間欠的）であり且つ離散的であるデータ、或いは言い換えると、値の変化の生起が時系列において散発的（間欠的）であり且つ値の変化の態様が離散的であるデータに基づいて計測時刻毎のガスの流速ｙtを算定することができるので、ガスの総消費量の実態の再現性の向上を図ることが可能になる。
以上のように構成されたガス消費量の用途分解方法、用途分解装置、及び用途分解プログラムによれば、さらに、実態が良好に再現されたガスの流速ｙtに基づいてガスの総消費量の用途を給湯分，暖房分，及び厨房分に分類することができるので、ガスの総消費量の用途分解の良好な精度を確保してガス消費量の用途分解技術としての信頼性の向上を図ることが可能になる。
なお、各処理において用いられる各閾値は、特定の値に限定されるものではなく、適宜設定され得る。
次に、図５乃至図１０に、本発明の居住者の不在期間の確定方法、装置、及びプログラムの実施形態の一例を示す。
ここで、本発明の居住者の不在期間の確定方法における在宅不在の判定とは、住居や住戸内に居住者が存在している状態（即ち、在宅である）と存在していない状態（即ち、不在である）とのどちらであるかを推定することである。
居住者の在宅不在の判定では、図３に示される連続データの作成方法の考え方が利用される。そして、所定の時間間隔且つ所定の計量単位で出力される計測値から実際の電力需要量が推定され、当該推定された電力需要量が用いられて居住者の在宅不在の判定ひいては不在期間の確定が行われる。
具体的には、本実施形態の居住者の在宅不在の判定方法は、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量ｌを計測単位として計量する電力需要量との組み合わせデータとして取得（Ｓ１）された電力需要量の計測データから、所定の時間帯の組み合わせデータが分析対象区間のデータとして抽出され、分析対象区間における計測時刻毎の電力需要量の累積値を下限とすると共に計測時刻毎の電力需要量の累積値に所定の電力量ｌを加えた値を上限とする領域内に存在し且つ分析対象区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測された電力需要量の総和と推定される各計測時刻の電力需要量の総和とが一致する折れ線が推定され（Ｓ２）、当該折れ線における各計測時刻間の直線の傾きとして電力需要速度ｘtが算定され（Ｓ２）、在宅不在の判定の対象である計測時刻以前の所定の期間に含まれる電力需要速度ｘtが用いられて前記所定の期間における電力需要速度ｘtの統計量が特徴量として算出され、少なくとも前記特徴量が用いられて居住者の在宅不在が判定される（Ｓ３）ようにしている（図５参照）。
また、本実施形態の居住者の在宅不在の判定装置は、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量ｌを計測単位として計量する電力需要量との組み合わせデータとして取得された電力需要量の計測データの入力を受ける手段としてのデータ受部３１ａと、計測データから、所定の時間帯の組み合わせデータを分析対象区間のデータとして抽出し、分析対象区間における計測時刻毎の電力需要量の累積値を下限とすると共に計測時刻毎の電力需要量の累積値に所定の電力量ｌを加えた値を上限とする領域内に存在し且つ分析対象区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測された電力需要量の総和と推定される各計測時刻の電力需要量の総和とが一致する折れ線を推定して当該折れ線における各計測時刻間の直線の傾きとして電力需要速度を算定する手段としての最短折線作成部３１ｂと、在宅不在の判定の対象である計測時刻以前の所定の期間に含まれる電力需要速度を用いて前記所定の期間における電力需要速度の統計量を特徴量として算出すると共に少なくとも前記特徴量を用いて居住者の在宅不在を判定する手段としての時刻別判定部３１ｃとを有する。
上記居住者の在宅不在の判定方法及び居住者の在宅不在の判定装置は、居住者の在宅不在の判定プログラムがコンピュータ上で実行されることによっても実施・実現され得る。ここでは、居住者の在宅不在の判定プログラムがコンピュータ上で実行されることによって居住者の在宅不在の判定方法が実施されると共に居住者の在宅不在の判定装置が実現される場合を説明する。
本実施形態の居住者の在宅不在の判定及び不在期間の確定プログラム３７を実行するためのコンピュータ３０（本実施形態では、居住者の在宅不在の判定及び不在期間の確定装置３０でもある）の全体構成を図７に示す。
このコンピュータ３０は制御部３１，記憶部３２，入力部３３，表示部３４，及びメモリ３５を備え、これらが相互にバス等の信号回線によって接続されている。
制御部３１は、記憶部３２に記憶されている居住者の在宅不在の判定プログラム３７に従ってコンピュータ３０全体の制御並びに居住者の在宅不在の判定に係る演算を行うものであり、例えばＣＰＵ（中央演算処理装置）である。
記憶部３２は、少なくともデータやプログラムを記憶可能な装置であり、例えばハードディスクである。
入力部３３は、少なくとも作業者の命令や種々の情報を制御部３１に与えるためのインターフェイス（即ち、情報入力の仕組み）であり、例えばキーボードやマウスである。なお、例えばキーボードとマウスとの両方のように複数種類のインターフェイスを入力部３３として有するようにしても良い。
表示部３４は、制御部３１の制御によって文字や図形或いは画像等の描画・表示を行うものであり、例えばディスプレイである。
メモリ３５は、制御部３１が種々の制御や演算を実行する際の作業領域であるメモリ空間となるものであり、例えばＲＡＭ（Ｒandom Ａccess Ｍemory の略）である。
そして、コンピュータ３０（以下、「居住者の在宅不在の判定装置３０」と表記する）の制御部３１には、居住者の在宅不在の判定プログラム３７が実行されることにより、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量ｌを計測単位として計量する電力需要量との組み合わせデータとして取得された電力需要量の計測データの入力を受ける処理を行うデータ受部３１ａと、計測データから所定の時間帯の組み合わせデータを分析対象区間のデータとして抽出し、さらに、分析対象区間における計測時刻毎の電力需要量の累積値を下限とすると共に計測時刻毎の電力需要量の累積値に所定の電力量ｌを加えた値を上限とする領域内に存在し且つ分析対象区間の始点から終点へと至る総長が最短であり尚且つ各計測時刻において計測された電力需要量の総和と推定される各計測時刻の電力需要量の総和とが一致する折れ線を推定して当該折れ線における各計測時刻間の直線の傾きとして電力需要速度ｘtを算定する処理を行う最短折線作成部３１ｂと、在宅不在の判定の対象である計測時刻以前の所定の期間に含まれる電力需要速度ｘtを用いて所定の期間における電力需要速度ｘtの統計量を特徴量として算出すると共に少なくとも特徴量を用いて居住者の在宅不在を判定する処理を行う時刻別判定部３１ｃとが構成される。
ここで、以下に説明する処理を行う際に用いられる変数についての計算機における処理として、最初に、不在開始時刻，不在終了時刻，及び不在区間指標を表す変数は全て初期化され、具体的には、不在開始時刻及び不在終了時刻を表す変数には例えば時刻としてはあり得ない数値が与えられたり、不在区間指標を表す変数には例えば０(ゼロ)が与えられたりする。
そして、居住者の在宅不在の判定方法の実施として、まず、電力需要量の計測が行われて計測データの取得が行われる（Ｓ１）。
電力需要量の計測の仕方は、在宅不在の判定の対象である住宅や住戸における電力需要量を計測して所定の計量単位（言い換えると、所定の計量の刻み）で所定の時間毎に例えば電気的な信号として前記所定の計量単位での需要量を出力するものであれば、特定の方法や機器によるものに限定されるものではなく、前記のような計測を行い得る方法や機器が適宜選択される。
具体的には例えば、あくまで一例として挙げると、電力需要量の計測の仕方（具体的には、計測機器）としてスマートメータ（例えば、上述の非特許文献２）が利用され得る。
電力需要量の計測における上述の「所定の計量単位」としての、１単位として扱われる（言い換えると、１単位として計量される）実際の電力量は、特定の値［Ｗｈ］に限定されるものではなく、例えば計測に用いられる計測機器の仕様などが踏まえられた上で適当な値に適宜設定される。具体的には例えば、スマートメータのＡルート計測値では、１単位として扱われる実際の電力量は１００ Ｗｈ である。
また、電力需要量の計測における上述の「所定の時間」としての、計測値を出力する時間間隔は、特定の時間［分］に限定されるものではなく、例えば計測に用いられる計測機器の仕様などが踏まえられた上で適当な時間に適宜設定される。具体的には例えば、スマートメータのＡルート計測値では、計測値が出力される時間間隔は３０分である。
本実施形態では、電力需要量の計測においてスマートメータが利用され、上述の「所定の計量単位」が１００ Ｗｈ であると共に上述の「所定の時間」が３０分である計測データ、すなわち、実際の電力量１００ Ｗｈ を１単位として計量する電力需要の値（「需要単位計量値」と呼ぶ）が３０分間隔で出力される計測データが、Ｓ１の処理によって取得されるとする。
そして、計測によって取得された電力需要量の計測値（具体的には、需要単位計量値）は、制御部３１のデータ受部３１ａを介して居住者の在宅不在の判定装置３０へと入力される。
電力需要量の計測値（需要単位計量値）は、居住者の在宅不在の判定装置３０に、各種記録媒体やデータサーバを介して入力されるようにしても良く、或いは、例えばスマートメータのように計測機器が備えている通信機能または計測機器に付加された通信機能や通信機序によって入力されるようにしても良い。
記録媒体を介して計測値が入力される場合には、例えば、計測機器から出力された需要単位計量値が記録媒体に記録され、当該記録媒体が居住者の在宅不在の判定装置３０の記録媒体接続用端子（図示していない）に差し込まれることによって需要単位計量値が入力されるようにすることが考えられる。
また、データサーバを介して計測値が入力される場合には、例えば、データサーバがバス等の信号回線によって居住者の在宅不在の判定装置３０に接続され、計測機器から出力された需要単位計量値が前記データサーバにデータファイル等として格納(保存)され、当該データファイル等として保存された需要単位計量値が読み込まれるようにすることが考えられる。
また、例えばスマートメータのように計測機器が備えている通信機能または計測機器に付加された通信機能や通信機序によって計測値が入力される場合には、需要単位計量値がリアルタイムで（即ち、即時的に）居住者の在宅不在の判定装置３０に入力されることも可能であり、したがってほぼリアルタイムでの（即ち、即時的な）居住者の在宅不在の判定の処理も可能である。
なお、記録媒体，データサーバ，及び通信機能や通信機序が組み合わされて用いられて需要単位計量値が居住者の在宅不在の判定装置３０に入力されるようにしても良い。
また、電力需要量の各計測値（各需要単位計量値）が計測機器から出力された時刻が計測機器に備えられた時計機能によって計測値と一緒に出力されて当該計測値と対応づけられて記録されたり、或いは、計測機器から出力された各計測値が記録媒体，データサーバ，或いは居住者の在宅不在の判定装置３０に記録されたり入力されたりした時の時刻が当該計測値と対応づけられて記録されたりする。なお、この時刻の間隔［分］が上述の「所定の時間」になる。
ここでの説明においては、各計測値（各需要単位計量値）と対応づけられて記録される時刻のことを「計測時刻」と呼び、また、需要単位計量値と計測時刻との組み合わせデータのことを「計測データ」と呼ぶ。
そして、データ受部３１ａにより、計測データがメモリ３５に記憶させられる。
次に、制御部３１の最短折線作成部３１ｂにより、Ｓ１の処理によって計測されて取得された計測データが用いられて累積需要曲線の作成が行われる（Ｓ２）。
Ｓ２の処理は、分析対象区間について実行される。分析対象区間は、少なくとも、対象住宅や対象住戸についての居住者の在宅不在の判定が行われる時刻（「判定対象時刻」と呼ぶ）を含む期間として設定される。
分析対象区間は、特定の時間長さ［時間］に限定されるものではなく、例えば検知の主要な対象として想定される不在期間について想定される時間が考慮されるなどして適当な時間長さに適宜設定される。分析対象区間は、具体的には例えば、あくまで一例として挙げると、６〜３６時間程度の範囲で適当な時間長さに設定され得る。また、分析対象区間は、例えばリアルタイムでの（即ち、即時的な）居住者の在宅不在の判定が行われる場合などには、判定対象時刻の進行に伴って当該の判定対象時刻を基準とするなどして適時・適宜に更新されて（言い換えると、設定し直されて）時間軸方向に移動させられるようにしても良い。
なお、判定対象時刻は、全ての計測時刻の各々が順次に判定対象時刻とされるようにしても良く、或いは、各計測時刻のうちから所定の間隔で選択された計測時刻が順次に判定対象時刻とされるようにしても良い。
Ｓ１の処理によって取得される計測データは、３０分間隔での１００ Ｗｈ 単位（言い換えると、１００ Ｗｈ 刻み）の電力需要量である。しかしながら、何らかの電気機器が使用されているときは、さらに、一般の住居では冷蔵庫など常時稼働且つ所定の間隔（或いは、所定のトリガー）で電力需要が変動する電気機器が通常は使用されていることを考慮すれば、電力（待機電力等を含む）の需要が常に継続しているという点において電力需要量の変化の生起は時系列において連続であり、また、電力が常に消費されているという点において電力需要量の変化の態様は漸次的である。
つまり、Ｓ１の処理における計測では、実際には(値の)変化の生起が時系列において連続であると共に(値の)変化の態様が漸次的である現象が、所定の時間間隔且つ所定の値単位（言い換えると、所定の値刻み，所定の計測単位，或いは所定の計量単位）という散発的（言い換えると、間欠的）且つ離散的な変化として把握されていると言える。
そして、本実施形態では、所定の時間間隔且つ所定の値単位（言い換えると、所定の値刻み，所定の計測単位，或いは所定の計量単位）という計測値の変化が散発的（言い換えると、間欠的）であり且つ離散的である計測データがそのまま用いられるのではなく、当該散発的（間欠的）且つ離散的な計測データが、計測値の変化の生起が時系列において連続であり且つ計測値の変化の態様が漸次的であるデータに変換された上で在宅不在の判定の処理に用いられる。
本実施形態の居住者の在宅不在の判定方法、判定装置、及び判定プログラムにおいて用いられる、計測値の変化の生起が時系列において散発的（間欠的）であり且つ計測値の変化の態様が離散的であるデータに基づいて、計測値の変化の生起が時系列において連続であり且つ計測値の変化の態様が漸次的である累積需要曲線が作成される、連続データの作成方法を、図８を用いて説明する。そして、この作成方法の考え方（言い換えると、要点）は、図３を用いて説明した作成方法の考え方（要点）と同様である。
ここでの説明の例として、図８に示すような、或る時間帯の計測データが、具体的には例えばスマートメータのＡルート計測値のように３０分毎に出力される需要単位計量値の累積値であって１００ Ｗｈ 単位の電力需要量の累積値（図中の●印）がＳ１の処理によって取得されて分析対象区間分の計測データとして抽出されたとする。
図８の横軸は、３０分間隔である計測時刻の各時点を、当該の分析対象区間の最初の計測時刻を１とした上で識別子としての連続番号で表したものである。すなわち、横軸における例えば「３」は、当該の分析対象区間の始まりから１．５時間（＝９０分＝３０分×３）だけ時間が経過していることを表す。
また、図３の縦軸は、当該の分析対象区間における需要単位計量値の累積値を表す。すなわち、縦軸における例えば「４００」は、当該の分析対象区間の始まりから４００［Ｗｈ］だけ電力が消費されていることを表す。
Ｓ１の処理によって取得される計測データは所定の時間間隔且つ所定の値単位（値刻み，計測単位，計量単位）という値の変化が散発的（間欠的）且つ離散的なデータであるので、すなわち、計測値の変化の生起が時系列において散発的（間欠的）且つ計測値の変化の態様が離散的なデータであるので、計測データの需要単位計測値の累積値のデータプロットである●印を時系列で結ぶと図８に示すように不規則な昇り階段状になる（図８中において符号４１；「下限階段直線４１」と呼ぶ）。
また、計測データの需要単位計量値の累積値に計測時刻毎に１００ Ｗｈ（即ち、需要単位計量値における計量単位であり値刻みの大きさである）を加えたデータのプロットである▲印を時系列で結んだものも、図８に示すように不規則な昇り階段状になる（図８中において符号４２；「上限階段直線４２」と呼ぶ）。
そして、下限階段直線４１が実際の電力需要量に基づく真の累積値の下限（言い換えると、実際の電力需要量が取り得る累積値の下限）であり、上限階段直線４２が実際の電力需要量に基づく真の累積値の上限（言い換えると、実際の電力需要量が取り得る累積値の上限）である。
したがって、実際の電力需要量の時系列における累積を表す真の累積需要曲線は、下限階段直線４１と上限階段直線４２とに挟まれた領域（「存在可能領域」と呼ぶ）内に存在し、全体としては増加傾向を示す非減少の曲線であると考えられる。
そして、真の累積需要曲線としては、存在可能領域内において増加傾向を示す非減少の折れ線のうち、折れ点（言い換えると、変曲点、若しくは、電力需要速度変化点）の個数が比較的少ないものを想定する。
そこで、実際の電力需要量の時系列における累積を表す真の累積需要曲線として、存在可能領域内に存在する、総長（即ち、当該の分析対象区間の始まりから終わりまでに亙る軌跡の長さ）が最短の折れ線（「最短折れ線」と呼ぶ）を推定する。
真の累積需要曲線としての、存在可能領域内に存在する、最短折れ線の推定手法は、特定の方法に限定されるものではなく、特定の領域内に存在して所定の始点から終点へと至る折れ線のうちの総長が最短の折れ線であり且つ各計測時刻において計測された電力需要量の総和と推定される各計測時刻の電力需要量の総和とが一致する折れ線を推定し得る適当な方法が適宜選択される。
最短折れ線の推定手法として、具体的には例えば、あくまで一例として挙げると、当該の分析対象区間において確定済みの折れ線の最後の点を起点とした（未確定の）直線部の許容される傾きの上限値・下限値を順次更新し、全ての直線が存在可能領域を逸脱せざるを得なくなる時点で直線の傾きを決定し、その直線が存在可能領域の境界に最初に接した時点までの折れ線の確定を行い、そして、新たに確定した折れ線の最後の点を新たな起点に設定した上で上述と同様の処理を繰り返すことを基本的な考え方とする方法が用いられ得る。
この推定方法は、計測時刻ｔ（尚、計測単位は３０分であり、したがって計測時刻の間隔は３０分である）における計側電力量をｘt［Ｗｈ／３０分］とすると共に、推定電力需要速度をｙt［Ｗｈ／３０分］とすると（図９参照）、最適化問題として数式７のように表される。
数式７においてｔは分析対象区間における計測時刻の識別子であってｔ＝１，２，３，…であり、数式７では当該の分析対象区間の範囲として計測時刻の識別子ｔが１からＴまでの期間に亙って処理が実行される。
そして、計測時刻を１刻みの識別子とすることにより、計測の時間間隔である３０分を一単位として扱うことになり、計測時刻 ｔ−１ から ｔ にかけての増分は計測の時間間隔当たりの増加量であって計測時刻 ｔ−１ から ｔ にかけての速度［Ｗｈ／３０分］であり、真の累積需要曲線としては計測時刻 ｔ−１ から ｔ にかけての直線の傾きである。
数式７において、また、計測電力量ｘt［Ｗｈ／３０分］は、計測時刻ｔにおける需要単位計量値（言い換えると、計測時刻 ｔ−１ から ｔ にかけての需要単位計量値の増分）である。
数式７において、さらに、ｌ(エル)は、電力需要［Ｗｈ］の最小単位（即ち、計測単位，計量単位）であり、本実施形態ではｌ＝１００である。
そして、数式７における最小化の対象は、経過時間（即ち、各々１）と各経過時間に関する推定需要（即ち、推定電力需要速度ｙt）とからなる折れ線の距離である。折れ線の距離は、図９に示す例のように、経過時間（即ち、各々１）と当該経過時間に関する推定電力需要速度ｙtとに対する斜辺の長さの総和として求められ、本発明ではこの斜辺の長さの総和が最も小さくなる曲線が推定累積需要曲線とされる。
制約条件の第一式は電力需要速度の非負性に対応し、また、第二式は実際の電力需要量の累積値の範囲制約に対応し、さらに、第三式は当該の分析対象区間における計測値ベースの電力需要量の総和と推定値ベースの電力需要量の総和との一致に対応する。
上述の最適化問題を解く方法によると、最悪の場合でも、時系列長（即ち、計測時刻の個数）の二乗のオーダーで真の累積需要曲線に相当する最短折れ線が推定される。
図８に示す計測データのプロット（図中の●印）を結んだ下限階段直線４１に対して数式７で表される最適化問題を適用した結果として（言い換えると、数式７に示す三つの制約条件を満たす累積需要曲線の全体が）、図８において破線で示す、存在可能領域内の単調増加曲線として、実際の電力需要量に近いと考えられる推定累積需要曲線４３が得られる。
ここで、上記の最適化問題は計測時刻の識別子ｔ＝１から処理が実行されるところ、ｔ＝１の処理が実行されるためには、ｔ＝１での起点としてｔ＝０（即ち、計測時刻の識別子ｔ＝１に対応する計測時刻の一つ前の計測時刻）からｔ＝１にかけての折れ線の最後の点が特定される必要が、すなわち、ｔ＝０からｔ＝１にかけての直線が特定される必要がある。この直線は、具体的には、当該の分析対象区間の最初の計測時刻（即ち、計測時刻の識別子ｔ＝１）の直前の計測時刻として識別子ｔを０(ゼロ)とすると共に需要単位計量値の累積値を０(ゼロ)とする仮想の原点を置き、当該原点と計測時刻の識別子ｔ＝１における需要単位計量値の累積値の値の点とを結ぶことによって特定される。
あるいは、当該の分析対象区間についての計測時刻の識別子ｔ＝０からｔ＝１にかけての直線の傾きが、識別子ｔ＝１において連続する直線の傾きに一致させられるようにしても良い。
ここで、真の累積需要曲線の上述の推定方法（即ち、数式７で表される最適化問題）は、目的関数部を「距離の和」から「距離の２乗和」に変更した数式８で表される最適化問題、すなわち、「存在可能領域の任意の曲線のうち各計測時刻の電力需要速度（即ち、電力需要量の累積値の計測時刻間差分）の２乗和が最小の曲線」を求める問題と等価である。
数式８で表される最適化問題では、２乗和の最小化により、各計測時刻の電力需要速度（電力需要量）が同じ値を持つ傾向が生じるので、折れ線が最適解になる。
真の累積需要曲線の推定方法としての数式７で表される最適化問題は、また、数式８で表される最適化問題の双対問題である数式９で表される最適化問題とも等価である。なお、数式９で表される最適化問題は、全変動正規化問題と呼ばれる凸最適化問題の一種である。なお、数式９において、ｘ0は計測電力量ｘt［Ｗｈ／３０分］の初期値（即ち、ｔ＝０のときの値）であり、図８，図９に示す例ではｘ0＝０である。
上述の真の累積需要曲線（即ち、最短折れ線）は、一例として以下に示すアルゴリズムが用いられることによって求められ得る。
−真の累積需要曲線の決定アルゴリズム−
《保持・更新情報》
○起点 折れ線
○上限直線 上限需要と領域境界接触点
○下限直線 下限需要と領域境界接触点
《更新内容》
(1) 時刻更新：時刻＝時刻＋１
(2) 下限直線＜上限累積需要値
起点から下限境界接触点までの需要量＝需要下限（確定）
起点＝下限直線の境界接触点
時刻＝起点，上限／下限直線の再設定
(1)へ
(3) 上限直線＞下限累積需要
起点から上限境界接触点までの需要量＝需要上限（確定）
起点＝上限直線の境界接触点
時刻＝起点，上限／下限直線の再設定
(1)へ
(4) 下限直線＜領域下限
下限需要，下限境界設定更新
上限直線＞領域上限
上限需要，上限境界設定更新
(5) 計測時刻の識別子ｔ＝Ｔになるまで(1)へ戻る
本実施形態では、最短折線作成部３１ｂにより、Ｓ１の処理においてメモリ３５に記憶された計測データ（即ち、需要単位計量値と計測時刻との組み合わせデータ）のうち、分析対象区間に該当する期間の計測データが時系列に従って順次読み込まれる。
そして、メモリ３５から読み込まれた計測時刻は、当該の分析対象区間における計測時刻の識別子として時系列順に１，２，３…という連続番号に変換される。
続いて、最短折線作成部３１ｂにより、上記の変換が施された計測データが用いられて数式７乃至数式９に表される最適化問題のうちのいずれかが解かれて、計測時刻の識別子ｔ毎の推定電力需要速度ｙt［Ｗｈ／３０分］が算定される。
そして、最短折線作成部３１ｂにより、当該の分析対象区間について、算定された推定電力需要速度ｙt［Ｗｈ／３０分］の値が計測時刻の識別子ｔの値と対応づけられて（言い換えると、計測時刻の識別子ｔの値と推定電力需要速度ｙtの値との組み合わせデータとして）メモリ３５に記憶させられる。
次に、制御部３１の時刻別判定部３１ｃにより、Ｓ２の処理によって算定された計測時刻別の推定電力需要速度が用いられて居住者の在宅不在の時刻別判定指標の決定が行われる（Ｓ３）。
Ｓ３の処理は、Ｓ２の処理が実行された分析対象区間に含まれる判定対象時刻の各々について実行される。なお、判定対象時刻に該当する計測時刻の識別子をｔeと表記する。
本発明では、判定対象時刻における在宅不在の判定を行う際に、当該の判定対象時刻以前の一定期間内の推定電力需要速度のデータが用いられて平均値や最大値などの統計量が特徴量として算出され、当該特徴量が用いられて在宅か不在かの判定が行われる。特徴量の算出対象期間であって、当該の判定対象時刻以前の一定期間のことを「特徴量期間」と呼ぶ。
特徴量期間は、特定の時間長さ［分，時間］に限定されるものではなく、例えば事前の分析結果が考慮されるなどして適当な時間長さに適宜設定される。特徴量期間は、具体的には例えば、あくまで一例として挙げると、１〜３時間程度の範囲で適当な時間長さに設定され得る。
具体的には例えば、或る判定対象時刻（計測時刻の識別子ｔ＝ｔe）における判定を行う場合に、特徴量期間が１時間であるときは、当該の判定対象時刻（計測時刻の識別子ｔ＝ｔe）における推定電力需要速度ｙte［Ｗｈ／３０分］（尚、当該の判定対象時刻以前３０分間の状態が反映されている）と当該の判定対象時刻の一つ前の計測時刻（計測時刻の識別子ｔ＝ｔe−１）における推定電力需要速度ｙte-1［Ｗｈ／３０］（尚、当該の判定対象時刻の６０分前から３０分前までの状態が反映されている）とが用いられてこれら二つの値についての統計量が特徴量として算出される。
特徴量としての統計量は、特定の種類に限定されるものではなく、例えば事前の分析結果が考慮されるなどして適当な種類が取り上げられて適宜設定される。特徴量としては、具体的には例えば、平均値（即ち、特徴量期間に含まれる計測時刻別の推定電力需要速度の平均の値），最大値（即ち、特徴量期間に含まれる計測時刻別の推定電力需要速度のうちの最も大きな値），最小値（即ち、特徴量期間に含まれる計測時刻別の推定電力需要速度のうちの最も小さな値），範囲（即ち、特徴量期間に含まれる計測時刻別の推定電力需要速度のうちの最も大きな値と最も小さな値との差の値），標準偏差（即ち、特徴量期間に含まれる計測時刻別の推定電力需要速度の標準偏差の値）などが挙げられる。
なお、Ｓ３の処理において用いられる特徴量の個数は、特定の数に限定されるものではなく、一個でも良く、或いは、複数個でも良い。
また、特徴量として、電力需要量（具体的には、推定電力需要速度）に関する各種統計量に加え、種々の情報が追加的に用いられるようにしても良い。この場合の情報の種類としては、例えば、世帯の識別子，判定対象時刻（実際の時刻），気温，季節などが挙げられる。
そして、上記のようにして算出されたり設定されたりした特徴量（言い換えると、説明変数）が用いられて、在宅と不在との２クラス分類が行われることによって在宅不在の判定が行われる。
２クラス分類を行う手法としては、例えばランダムフォレストやサポートベクターマシン（ＳＶＭ：Ｓupport Ｖector Ｍachine の略）などの機械学習手法が用いられ得る。
なお、機械学習手法が用いられる際には、時刻別の電力需要量の実測データと当該時刻における実際の在宅不在に関する状態データとの組み合わせデータが正解データ（「教師データ」や「訓練データ」とも呼ばれる）が必要とされる。この正解データは、例えば、実態調査，モニタ調査，アンケート調査などが行われて収集されたり、既存の観測データが利用されて準備されたりする。そして、正解データが機械学習手法における学習に用いられて２クラス分類の判定基準が作成される。なお、２クラス分類に用いられる機械学習手法の手順や判定基準の内容などは、例えば、居住者の在宅不在の判定プログラム３７内に予め規定される。
本実施形態では、時刻別判定部３１ｃにより、Ｓ２の処理においてメモリ３５に記憶された分析対象区間についての計測時刻の識別子ｔの値と推定電力需要速度ｙtの値との組み合わせデータのうち、当該の判定対象時刻についての特徴量期間に該当する（言い換えると、特徴量期間に含まれる）計測時刻の識別子と対応づけられている全ての推定電力需要速度の値が読み込まれる。
そして、時刻別判定部３１ｃにより、読み込まれた推定電力需要速度の値が用いられて、２クラス分類に用いられる特徴量が算出される。
時刻別判定部３１ｃにより、さらに、機械学習手法が用いられて２クラス分類が行われ、当該の判定対象時刻について対象住宅や対象住戸に人が存在しているか存在していないかが決定される。
そして、時刻別判定部３１ｃにより、２クラス分類の結果（即ち、「在宅」と「不在」とのうちのどちらであるかの決定結果）が、在宅不在の時刻別判定指標として判定対象時刻（計測時刻の識別子ｔ＝ｔe）と対応づけられて（言い換えると、計測時刻の識別子の値と在宅不在の時刻別判定指標を表す変数の値との組み合わせデータとして）メモリ３５に記憶させられる。
なお、計算機における処理としては、２クラス分類の結果が、「在宅」である場合には時刻別判定指標を表す変数に「在宅」に対応する数値として適宜選択された適当な数値（具体的には例えば、１）が与えられ、「不在」である場合には時刻別判定指標を表す変数に「不在」に対応する数値として適宜選択された適当な数値（具体的には例えば、０(ゼロ)）が与えられる。
本実施形態では、引き続いて、上述までの処理によって決定された在宅不在の時刻別判定指標が用いられて不在開始時刻及び不在終了時刻が特定されて不在期間が確定されるようにしている（Ｓ４乃至Ｓ１１；図６，図１０）。
具体的にはまず、制御部３１の不在区間判断部３１ｄにより、Ｓ３の処理によって決定された当該の判定対象時刻における在宅不在の時刻別判定指標が用いられて、当該時刻別判定指標が「不在」であるか否かの判断が行われる（Ｓ４）。
本実施形態では、不在区間判断部３１ｄにより、Ｓ３の処理においてメモリ３５に記憶された当該の判定対象時刻の計測時刻の識別子ｔeと対応づけられている在宅不在の時刻別判定指標が読み込まれる。
そして、読み込まれた時刻別判定指標が「不在」である場合には（Ｓ４：Ｙｅｓ）、処理ステップがＳ５の処理へと進められる。
そして、Ｓ５の処理として、不在区間判断部３１ｄにより、後述するＳ８の処理において設定・操作される不在区間指標が「不在」であるか否かの判断が行われる。
そして、不在区間指標が「不在」でない場合には（Ｓ５：Ｎｏ）、不在開始時刻がＳ１の処理の前の初期化された状態のままであるときは、当該の判定対象時刻の計測時刻の識別子ｔ（＝ｔe）の一つ前の計測時刻の識別子ｔ（＝ｔe−１）に対応する実際の時刻が不在開始時刻としてメモリ３５に記憶させられる（Ｓ６）。なお、不在開始時刻がＳ１の処理の前の初期化された状態のままでないときは、不在開始時刻に対しては何らの操作もしない。
次に、制御部３１の在宅不在判定部３１ｅにより、Ｓ３の処理によって決定された在宅不在の時刻別判定指標が用いられて、居住者の在宅不在の判定が行われる（Ｓ７）。
Ｓ３の処理によって決定される在宅不在の時刻別判定指標は、判定対象時刻毎に各々独立して判定された結果であり、居住者が在宅しているものの、短時間、使用する電気機器が少なかったり電気機器の使用態様に変化が無かったりする場合に（言い換えると、電力需要速度の変化が小さい場合に）、不在であるとの判断結果になる虞がある。
このため、本実施形態では、居住者の在宅不在の最終的な判定は、在宅不在の時刻別判定指標が当該の判定対象時刻以前の所定の時間に亙って連続して「不在」であるか否か（言い換えると、Ｓ３の処理における「不在」の決定が所定の時間に亙って継続しているか否か）が判断されることによって行われる。ここで、在宅不在の判定を行うための前記所定の時間のことを「判定期間」と呼ぶ。
判定期間は、特定の時間長さ［時間］（言い換えると、判定対象時刻の個数）に限定されるものではなく、例えば事前の分析結果が考慮されるなどして適当な時間長さに適宜設定される。判定期間は、具体的には例えば、１〜１２時間程度の範囲で適当な時間長さに設定されることが考えられ、４〜８時間程度のうちのいずれかの時間長さに設定されることが好ましく、６時間程度に設定されることが更に好ましい。
本実施形態では、在宅不在判定部３１ｅにより、Ｓ３の処理においてメモリ３５に記憶された計測時刻の識別子ｔの値と在宅不在の時刻別判定指標を表す変数の値との組み合わせデータのうち、当該の判定対象時刻の計測時刻の識別子ｔ（＝ｔe）と対応づけられている在宅不在の時刻別判定指標を表す変数の値が読み込まれると共に、当該の判定対象時刻についての判定期間に該当する（言い換えると、判定期間に含まれる）計測時刻の識別子ｔ（＝ｔe−１，ｔe−２ 等）と対応づけられている全ての在宅不在の時刻別判定指標を表す変数の値が読み込まれる。
そして、在宅不在判定部３１ｅにより、当該の判定対象時刻についての判定期間中の在宅不在の時刻別判定指標が全て「不在」であるか否かが判断される。
そして、判定期間中の在宅不在の時刻別判定指標の中に一つでも「在宅」が有る場合には（Ｓ７：Ｎｏ）、処理ステップがＳ１の処理へと戻される。
他方、判定期間中の在宅不在の時刻別判定指標が全て「不在」である場合には（Ｓ７：Ｙｅｓ）、当該の判定期間の間は居住者は不在であると判定されるということになり、在宅不在判定部３１ｅにより、不在区間指標が「不在」にされた上でメモリ３５に記憶させられると共に当該の判定対象時刻の計測時刻の識別子ｔ（＝ｔe）に対応する実際の時刻が不在終了時刻としてメモリ３５に記憶させられ（Ｓ８）、その上で処理ステップがＳ１の処理へと戻される。なお、不在区間指標が「不在」にされるとは、計算機における処理としては、不在区間指標を表す変数に「不在」に対応する数値として適宜選択された適当な数値、具体的には例えば１が与えられたりすることである。
また、Ｓ５の処理において、不在区間判断部３１ｄにより、不在区間指標が「不在」であるか否かの判断が行われて当該不在区間指標が「不在」である場合には（Ｓ５：Ｙｅｓ）、当該の判定対象時刻の計測時刻の識別子ｔ（＝ｔe）に対応する実際の時刻が不在終了時刻としてメモリ３５に記憶させられ（Ｓ９）、その上で処理ステップがＳ１の処理へと戻される。なお、Ｓ９の処理では、Ｓ８の処理若しくは前回のＳ９の処理において設定された不在終了時刻が更新されることになる。
また、Ｓ４の処理において、不在区間判断部３１ｄにより、当該の判定対象時刻における在宅不在の時刻別判定指標が「不在」であるか否かの判断が行われて当該時刻別判定指標が「不在」でない場合には（Ｓ４：Ｎｏ）、処理ステップがＳ１０の処理へと進められる。
そして、Ｓ１０の処理として、制御部３１の不在期間確定部３１ｆにより、不在区間指標が「不在」であるか否かの判断が行われ、不在区間指標が「不在」である場合には（Ｓ１０：Ｙｅｓ）、上述のＳ６の処理において特定された不在開始時刻から、上述のＳ８の処理若しくはＳ９の処理において特定された不在終了時刻までが、不在期間であると確定される（Ｓ１１）。
そして、対象住宅や対象住戸についての在宅不在の判定の結果として、必要に応じて個々の不在期間を区別するための識別子と対応づけられて、不在開始時刻と不在終了時刻との組み合わせデータが、データファイル等に記録されて記憶部３２などに格納(保存)されたり、表示部３４に表示されたりする。その上で、処理ステップがＳ１の処理へと戻される。この際、不在開始時刻，不在終了時刻，及び不在区間指標は全てリセットされ、計算機における処理としては、Ｓ１の処理の前の初期化された状態に戻される。
他方、不在区間指標が「不在」でない場合には（Ｓ１０：Ｎｏ）、処理ステップがＳ１の処理へと戻される。この際、不在開始時刻がＳ１の処理の前の初期化された状態のままでないときは、不在開始時刻はリセットされ、計算機における処理としては、Ｓ１の処理の前の初期化された状態に戻される。
以上の実施形態にかかる居住者の不在期間の確定方法、確定装置、及び確定プログラムによれば、所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量ｌ［Ｗｈ］を計測単位として計量する電力需要量との組み合わせデータ、すなわち、所定の時間間隔且つ所定の電力量ｌ刻みという値の変化が散発的（言い換えると、間欠的）であり且つ離散的であるデータ、或いは言い換えると、値の変化の生起が時系列において散発的（間欠的）であり且つ値の変化の態様が離散的であるデータに基づいて計測時刻毎の電力需要速度ｙtを算定することができるので、電力需要の実態の再現性の向上を図ることが可能なる。
さらに、実態が良好に再現された電力需要速度に基づいて住居や住戸に人が存在しているか存在していないかを推定することができるので、居住者の在宅不在の判定の良好な精度を確保して居住者の在宅不在の判定技術としての信頼性の向上を図ることが可能になる。
また、不在開始時刻から不在終了時刻までの確定された不在期間は、つまりは対象住宅や対象住戸において人の活動に伴う電気の使用が為されていない期間であるので、対象住宅や対象住戸には人が存在していないと理解したり、或いは、対象住宅や対象住戸に人が存在しているものの例えば病気や怪我などのために動けなくなっているような異常事態が発生している可能性があると理解したりすることができる。したがって、不在開始時刻から不在終了時刻までの時間長さすなわち不在期間が所定の閾値を超える場合を要安否確認事象として検知し、対象住宅や対象住戸に関して長期不在の有無を確認したり異常事態の発生の有無を確認したりするようにすることも考えられる。すなわち、本発明の居住者の不在期間の確定方法、装置、及びプログラムは、例えば、高齢者等の見守りシステムとしても利用されることが可能である。
なお、上述の形態は本発明を実施する際の好適な形態の一例ではあるものの本発明の実施の形態が上述のものに限定されるものではなく、本発明の要旨を逸脱しない範囲において本発明は種々変形実施可能である。
また、上述の実施形態におけるＳ４乃至Ｓ１１の処理の要点は、Ｓ３の処理によって決定された在宅不在の時刻別判定指標を用いて不在期間を確定し得るように、在宅不在の判定を行うための所定の時間としての判定期間が判定対象時刻の進行に伴って時間軸方向に移動させられながら判定期間毎に在宅不在が判定され、不在と判定された一つの判定期間の始まりを不在開始時刻とすると共に終わりを不在終了時刻として特定する、または、連続して不在と判定された複数の判定期間のうちの最初の判定期間の始まりを不在開始時刻とすると共に最後の判定期間の終わりを不在終了時刻として特定することである。したがって、時間軸方向に移動する判定期間毎に在宅不在を判定しながら不在開始時刻と不在終了時刻とを特定し得るものであれば、上述の実施形態におけるＳ４乃至Ｓ１１の処理に対応する具体的な方法・手順は上述の実施形態として説明したもの（図１０を含む）に限定されるものではない。
なお更に付け加えると、所定の時間長さを有すると共に時間軸方向に移動する判定期間を用いることによって不在開始時刻と不在終了時刻とを特定することの利点は、判定対象時刻毎に各々独立して（即ち、当該の判定対象時刻よりも前の時刻における状況を合わせて考慮すること無く）在宅不在を判定する場合のような、居住者が在宅しているにもかかわらず、短時間、使用する電気機器が少なかったり電気機器の使用態様に変化が無かったりする状況での（言い換えると、電力需要速度の変化が小さい状況での）不在であるとの誤った判断を避けることである。この考え方は、上述の実施形態のように値の変化が散発的（言い換えると、間欠的）であり且つ離散的であるデータ、或いは言い換えると、値の変化の生起が時系列において散発的（間欠的）であり且つ値の変化の態様が離散的であるデータに基づいて作成される連続データが用いられる場合のみに適用され得るものではなく、種々の特性のデータが用いられる場合に適用され得るものであり、さらに言えば、判定対象時刻別の在宅不在の判定の仕方は上述の実施形態における方法に限定されるものではない。
例えば、もとより連続であるデータ、言い換えると、値の変化の生起が時系列において連続であり且つ値の変化の態様が漸次的であるデータ、具体的には例えば、あくまで一例として挙げると、１分毎に出力される１ Ｗｈ 単位のデータが用いられる場合にも適用され得る。この場合には、上述の実施形態のように最適化問題を解くこと無く、あくまで一例として挙げると３０分間隔の時刻の識別子ｔ毎の電力需要速度ｙt［Ｗｈ／３０分］や１分間隔の時刻の識別子ｔ毎の電力需要速度ｙt［Ｗｈ／分］、或いは一般化するとｎ分間隔の時刻の識別子ｔ毎の電力需要速度ｙt［Ｗｈ／ｎ分］が得られる。そして、上述の実施形態におけるＳ３以降の処理が同様に行われ得る。
あるいは、値の変化の生起が時系列において散発的（間欠的）であり且つ値の変化の態様が漸次的であるデータ、具体的には例えば、あくまで一例として挙げると、３０分毎に出力される１ Ｗｈ 単位のデータが用いられる場合にも適用され得る。この場合には、上述の実施形態のように最適化問題を解くこと無く、あくまで一例として挙げると、３０分間隔の時刻の識別子ｔ毎の電力需要速度ｙt［Ｗｈ／３０分］が得られる。そして、上述の実施形態におけるＳ３以降の処理が同様に行われ得る。
３０ 居住者の在宅不在の判定並びに不在期間の確定装置
３７ 居住者の在宅不在の判定並びに不在期間の確定プログラム

Claims (13)

  1. ガス配管の主幹における計測によって所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量毎に１カウントとして計数する流量カウント数との組み合わせデータとして取得されたガス総消費量の計測データから、前記主幹から枝分かれした屋内配管に接続されているガス器具が使用されている時間帯の前記組み合わせデータが連続動作区間のデータとして抽出され、前記連続動作区間における前記計測時刻毎の前記流量カウント数の累積値と前記所定の流量とを掛け合わせた値を下限とすると共に前記掛け合わせた値に前記所定の流量を加えた値を上限とする領域内に存在し且つ前記連続動作区間の始点から終点へと至る総長が最短であり尚且つ各前記計測時刻において計測された前記ガス総消費量の総和と推定される各前記計測時刻のガス総消費量の総和とが一致する折れ線が推定され、当該折れ線における各前記計測時刻間の直線の傾きとして前記主幹におけるガスの流速が算定された上で、前記連続動作区間のデータが冬季であるか否か判断される処理と、前記ガスの流速が所定の範囲である状態が所定の時間に亙って継続しているか否かに基づいて暖房連続動作区間が特定される処理と、前記連続動作区間のデータにおいて特定された全ての前記暖房連続動作区間の中での前記ガスの流速の最小値が暖房ベース流速として設定される処理と、前記暖房ベース流速と給湯最小流速との和と前記ガスの流速との大きさが比較される処理と、前記給湯最小流速と前記ガスの流速との大きさが比較される処理と、厨房最小流速と間欠暖房最大流速とのうちの少なくとも一方と前記ガスの流速との大きさが比較される処理とが実行されることにより、前記ガス総消費量の用途が給湯分，暖房分，及び厨房分に分類されることを特徴とするガス消費量の用途分解方法。
  2. 前記ガスの流速が算定された後に、前記連続動作区間のデータが前記冬季であるか否か判断される処理がまず行われ、前記暖房連続動作区間が特定される処理が次に行われ、前記暖房ベース流速が設定される処理が続いて行われることを特徴とする請求項１記載のガス消費量の用途分解方法。
  3. ガス配管の主幹における計測によって所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量毎に１カウントとして計数する流量カウント数との組み合わせデータとして取得されたガス総消費量の計測データから、前記主幹から枝分かれした屋内配管に接続されているガス器具が使用されている時間帯の前記組み合わせデータを連続動作区間のデータとして抽出する手段と、前記連続動作区間における前記計測時刻毎の前記流量カウント数の累積値と前記所定の流量とを掛け合わせた値を下限とすると共に前記掛け合わせた値に前記所定の流量を加えた値を上限とする領域内に存在し且つ前記連続動作区間の始点から終点へと至る総長が最短であり尚且つ各前記計測時刻において計測された前記ガス総消費量の総和と推定される各前記計測時刻のガス総消費量の総和とが一致する折れ線を推定して当該折れ線における各前記計測時刻間の直線の傾きとして前記主幹におけるガスの流速を算定する手段と、前記連続動作区間のデータが冬季であるか否かを判断し、前記ガスの流速が所定の範囲である状態が所定の時間に亙って継続しているか否かに基づいて暖房連続動作区間を特定すると共に前記連続動作区間のデータにおいて特定された全ての前記暖房連続動作区間の中での前記ガスの流速の最小値を暖房ベース流速として設定し、前記暖房ベース流速と給湯最小流速との和と前記ガスの流速との大きさを比較し、前記給湯最小流速と前記ガスの流速との大きさを比較し、さらに、厨房最小流速と間欠暖房最大流速とのうちの少なくとも一方と前記ガスの流速との大きさを比較することにより、前記ガス総消費量の用途を給湯分，暖房分，及び厨房分に分類する手段とを有することを特徴とするガス消費量の用途分解装置。
  4. 前記ガスの流速を算定した後に、前記連続動作区間のデータが前記冬季であるか否かを判断する処理がまず行われ、前記暖房連続動作区間を特定すると共に前記暖房ベース流速を設定する処理が次に行われることを特徴とする請求項３記載のガス消費量の用途分解装置。
  5. ガス配管の主幹における計測によって所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の流量毎に１カウントとして計数する流量カウント数との組み合わせデータとして取得されたガス総消費量の計測データから、前記主幹から枝分かれした屋内配管に接続されているガス器具が使用されている時間帯の前記組み合わせデータを連続動作区間のデータとして抽出する処理と、前記連続動作区間における前記計測時刻毎の前記流量カウント数の累積値と前記所定の流量とを掛け合わせた値を下限とすると共に前記掛け合わせた値に前記所定の流量を加えた値を上限とする領域内に存在し且つ前記連続動作区間の始点から終点へと至る総長が最短であり尚且つ各前記計測時刻において計測された前記ガス総消費量の総和と推定される各前記計測時刻のガス総消費量の総和とが一致する折れ線を推定して当該折れ線における各前記計測時刻間の直線の傾きとして前記主幹におけるガスの流速を算定する処理と、前記連続動作区間のデータが冬季であるか否かを判断し、前記ガスの流速が所定の範囲である状態が所定の時間に亙って継続しているか否かに基づいて暖房連続動作区間を特定すると共に前記連続動作区間のデータにおいて特定された全ての前記暖房連続動作区間の中での前記ガスの流速の最小値を暖房ベース流速として設定し、前記暖房ベース流速と給湯最小流速との和と前記ガスの流速との大きさを比較し、前記給湯最小流速と前記ガスの流速との大きさを比較し、さらに、厨房最小流速と間欠暖房最大流速とのうちの少なくとも一方と前記ガスの流速との大きさを比較することにより、前記ガス総消費量の用途を給湯分，暖房分，及び厨房分に分類する処理とをコンピュータに行わせることを特徴とするガス消費量の用途分解プログラム。
  6. 前記ガスの流速を算定した後に、前記連続動作区間のデータが前記冬季であるか否かを判断する処理がまず行われ、前記暖房連続動作区間を特定すると共に前記暖房ベース流速を設定する処理が次に行われることを特徴とする請求項５記載のガス消費量の用途分解プログラム。
  7. 所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量を計測単位として計量する電力需要量との組み合わせデータとして取得された電力需要量の計測データから、所定の時間帯の前記組み合わせデータが分析対象区間のデータとして抽出され、前記分析対象区間における前記計測時刻毎の前記電力需要量の累積値を下限とすると共に前記計測時刻毎の前記電力需要量の累積値に前記所定の電力量を加えた値を上限とする領域内に存在し且つ前記分析対象区間の始点から終点へと至る総長が最短であり尚且つ各前記計測時刻において計測された前記電力需要量の総和と推定される各前記計測時刻の電力需要量の総和とが一致する折れ線が推定され、当該折れ線における各前記計測時刻間の直線の傾きとして電力需要速度が算定され、在宅不在の判定の対象である計測時刻以前の所定の期間に含まれる前記電力需要速度が用いられて前記所定の期間における前記電力需要速度の統計量が特徴量として算出され、少なくとも前記特徴量が用いられて居住者の在宅不在が判定されることを特徴とする居住者の在宅不在の判定方法。
  8. 前記電力需要量の計測データがスマートメータのＡルートから出力されるデータであることを特徴とする請求項７記載の居住者の在宅不在の判定方法。
  9. 所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量を計測単位として計量する電力需要量との組み合わせデータとして取得された電力需要量の計測データから、所定の時間帯の前記組み合わせデータを分析対象区間のデータとして抽出し、前記分析対象区間における前記計測時刻毎の前記電力需要量の累積値を下限とすると共に前記計測時刻毎の前記電力需要量の累積値に前記所定の電力量を加えた値を上限とする領域内に存在し且つ前記分析対象区間の始点から終点へと至る総長が最短であり尚且つ各前記計測時刻において計測された前記電力需要量の総和と推定される各前記計測時刻の電力需要量の総和とが一致する折れ線を推定して当該折れ線における各前記計測時刻間の直線の傾きとして電力需要速度を算定する手段と、在宅不在の判定の対象である計測時刻以前の所定の期間に含まれる前記電力需要速度を用いて前記所定の期間における前記電力需要速度の統計量を特徴量として算出すると共に少なくとも前記特徴量を用いて居住者の在宅不在を判定する手段とを有することを特徴とする居住者の在宅不在の判定装置。
  10. 前記電力需要量の計測データがスマートメータのＡルートから出力されるデータであることを特徴とする請求項９記載の居住者の在宅不在の判定装置。
  11. 所定の時間間隔の計測時刻と当該計測時刻毎に出力される値であって所定の電力量を計測単位として計量する電力需要量との組み合わせデータとして取得された電力需要量の計測データから、所定の時間帯の前記組み合わせデータを分析対象区間のデータとして抽出し、前記分析対象区間における前記計測時刻毎の前記電力需要量の累積値を下限とすると共に前記計測時刻毎の前記電力需要量の累積値に前記所定の電力量を加えた値を上限とする領域内に存在し且つ前記分析対象区間の始点から終点へと至る総長が最短であり尚且つ各前記計測時刻において計測された前記電力需要量の総和と推定される各前記計測時刻の電力需要量の総和とが一致する折れ線を推定して当該折れ線における各前記計測時刻間の直線の傾きとして電力需要速度を算定する処理と、在宅不在の判定の対象である計測時刻以前の所定の期間に含まれる前記電力需要速度を用いて前記所定の期間における前記電力需要速度の統計量を特徴量として算出すると共に少なくとも前記特徴量を用いて居住者の在宅不在を判定する処理とをコンピュータに行わせることを特徴とする居住者の在宅不在の判定プログラム。
  12. 前記電力需要量の計測データがスマートメータのＡルートから出力されるデータであることを特徴とする請求項１１記載の居住者の在宅不在の判定プログラム。
  13. 電力需要量の計測データに基づいて時刻間毎の電力需要速度が算定され、判定対象時刻以前の所定の期間に含まれる前記電力需要速度が用いられて前記所定の期間における前記電力需要速度の統計量が特徴量として算出され、少なくとも前記特徴量が用いられて前記判定対象時刻別の居住者の在宅不在が判定され、所定の時間長さを有して時間軸方向に移動する判定期間に対応する前記判定対象時刻別の判定結果が全て不在であるときに当該判定期間中は不在であると判定され、不在と判定された一つの前記判定期間の始まりの時刻を不在開始時刻とすると共に終わりの時刻を不在終了時刻として特定する、または、連続して不在と判定された複数の前記判定期間のうちの最初の前記判定期間の始まりの時刻を不在開始時刻とすると共に最後の前記判定期間の終わりの時刻を不在終了時刻として特定することを特徴とする不在期間の確定方法。
JP2020131315A 2015-05-13 2020-08-02 居住者の不在期間の確定方法、装置、及びプログラム Active JP6938735B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2015097846 2015-05-13
JP2015097846 2015-05-13
JP2017032219A JP6765763B2 (ja) 2015-05-13 2017-02-23 居住者の在宅不在の判定装置、及び判定プログラム

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2017032219A Division JP6765763B2 (ja) 2015-05-13 2017-02-23 居住者の在宅不在の判定装置、及び判定プログラム

Publications (2)

Publication Number Publication Date
JP2020201268A true JP2020201268A (ja) 2020-12-17
JP6938735B2 JP6938735B2 (ja) 2021-09-22

Family

ID=57549599

Family Applications (5)

Application Number Title Priority Date Filing Date
JP2016090861A Pending JP2016212104A (ja) 2015-05-13 2016-04-28 ガス消費量の用途分解方法、用途分解装置、及び用途分解プログラム、並びに、居住者の在宅不在の判定方法、判定装置、及び判定プログラム
JP2017032219A Active JP6765763B2 (ja) 2015-05-13 2017-02-23 居住者の在宅不在の判定装置、及び判定プログラム
JP2017032218A Active JP6713709B2 (ja) 2015-05-13 2017-02-23 ガス消費量の用途分解方法、用途分解装置、及び用途分解プログラム
JP2017036604A Pending JP2017201298A (ja) 2015-05-13 2017-02-28 エネルギー需要量の推定方法、推定装置、及び推定プログラム
JP2020131315A Active JP6938735B2 (ja) 2015-05-13 2020-08-02 居住者の不在期間の確定方法、装置、及びプログラム

Family Applications Before (4)

Application Number Title Priority Date Filing Date
JP2016090861A Pending JP2016212104A (ja) 2015-05-13 2016-04-28 ガス消費量の用途分解方法、用途分解装置、及び用途分解プログラム、並びに、居住者の在宅不在の判定方法、判定装置、及び判定プログラム
JP2017032219A Active JP6765763B2 (ja) 2015-05-13 2017-02-23 居住者の在宅不在の判定装置、及び判定プログラム
JP2017032218A Active JP6713709B2 (ja) 2015-05-13 2017-02-23 ガス消費量の用途分解方法、用途分解装置、及び用途分解プログラム
JP2017036604A Pending JP2017201298A (ja) 2015-05-13 2017-02-28 エネルギー需要量の推定方法、推定装置、及び推定プログラム

Country Status (1)

Country Link
JP (5) JP2016212104A (ja)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016212104A (ja) * 2015-05-13 2016-12-15 一般財団法人電力中央研究所 ガス消費量の用途分解方法、用途分解装置、及び用途分解プログラム、並びに、居住者の在宅不在の判定方法、判定装置、及び判定プログラム
JP6134458B1 (ja) * 2017-02-15 2017-05-24 東京瓦斯株式会社 通知システム、サーバ、通知方法
JP7091156B2 (ja) * 2018-06-12 2022-06-27 四国計測工業株式会社 パルス発生装置
JP7462259B2 (ja) * 2020-01-14 2024-04-05 パナソニックＩｐマネジメント株式会社 計測装置、システム、計測方法、及びプログラム
JP2021111222A (ja) * 2020-01-14 2021-08-02 パナソニックＩｐマネジメント株式会社 計測システム、学習装置、システム、計測方法、学習方法、及びプログラム
CN112097315B (zh) * 2020-09-22 2021-12-10 天津大学 适用于学校区域供热系统的分时分区互补间歇供暖方法
JP7122714B2 (ja) * 2020-12-18 2022-08-22 ｂｏｏｏｓｔ ｔｅｃｈｎｏｌｏｇｉｅｓ株式会社 在不在推定システム、電力制御システム、在不在推定方法および在不在推定プログラム

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005174249A (ja) * 2003-12-15 2005-06-30 Central Res Inst Of Electric Power Ind 電力需要家居住者の生活状況推定方法およびシステム
JP2012159896A (ja) * 2011-01-31 2012-08-23 Mitsubishi Electric Corp 監視装置及び監視方法
JP2013131097A (ja) * 2011-12-22 2013-07-04 Nippon Telegr & Teleph Corp <Ntt> 人間活動検知システム、人間活動検知方法および人間活動検知プログラム
JP2014074954A (ja) * 2012-10-02 2014-04-24 Nippon Telegr & Teleph Corp <Ntt> 人間活動の検知方法、及び人間活動の検知システム
JP2014225215A (ja) * 2013-04-16 2014-12-04 一般財団法人電力中央研究所 居住者安否確認システム
JP2015001878A (ja) * 2013-06-17 2015-01-05 一般財団法人電力中央研究所 生活状況の推定方法、推定装置及び推定プログラム

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3199536B2 (ja) * 1993-09-21 2001-08-20 東京瓦斯株式会社 ガス器具別使用量解析装置
JP4785272B2 (ja) * 2001-02-09 2011-10-05 東洋計器株式会社 多機能型ガスカウンタおよび多機能型ガスメータ
JP4170836B2 (ja) * 2002-09-06 2008-10-22 矢崎総業株式会社 流体計数装置、ガスメータ、並びに使用料金算出装置
JP5109636B2 (ja) * 2007-12-14 2012-12-26 パナソニック株式会社 流量計測装置とこれを用いたガス供給システム
JP5285980B2 (ja) * 2008-07-03 2013-09-11 矢崎エナジーシステム株式会社 ガス器具判別装置及びガス器具判別方法
JP5183515B2 (ja) * 2009-02-02 2013-04-17 東洋計器株式会社 ガス分計装置
JP5300761B2 (ja) * 2010-02-25 2013-09-25 株式会社東芝 エネルギー配分計算装置
GB2491109B (en) * 2011-05-18 2014-02-26 Onzo Ltd Identification of a utility consumption event
JP6084456B2 (ja) * 2012-12-14 2017-02-22 関西電力株式会社 ガス使用用途の判別システム
JP5997641B2 (ja) * 2013-03-25 2016-09-28 株式会社東芝 データ処理装置およびその方法、ならびにプログラム
US11079417B2 (en) * 2014-02-25 2021-08-03 Itron, Inc. Detection of electric power diversion
JP6292515B2 (ja) * 2014-04-08 2018-03-14 パナソニックＩｐマネジメント株式会社 需要予測装置、プログラム
JP2016212104A (ja) * 2015-05-13 2016-12-15 一般財団法人電力中央研究所 ガス消費量の用途分解方法、用途分解装置、及び用途分解プログラム、並びに、居住者の在宅不在の判定方法、判定装置、及び判定プログラム

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005174249A (ja) * 2003-12-15 2005-06-30 Central Res Inst Of Electric Power Ind 電力需要家居住者の生活状況推定方法およびシステム
JP2012159896A (ja) * 2011-01-31 2012-08-23 Mitsubishi Electric Corp 監視装置及び監視方法
JP2013131097A (ja) * 2011-12-22 2013-07-04 Nippon Telegr & Teleph Corp <Ntt> 人間活動検知システム、人間活動検知方法および人間活動検知プログラム
JP2014074954A (ja) * 2012-10-02 2014-04-24 Nippon Telegr & Teleph Corp <Ntt> 人間活動の検知方法、及び人間活動の検知システム
JP2014225215A (ja) * 2013-04-16 2014-12-04 一般財団法人電力中央研究所 居住者安否確認システム
JP2015001878A (ja) * 2013-06-17 2015-01-05 一般財団法人電力中央研究所 生活状況の推定方法、推定装置及び推定プログラム

Also Published As

Publication number Publication date
JP2017125856A (ja) 2017-07-20
JP2016212104A (ja) 2016-12-15
JP2017201298A (ja) 2017-11-09
JP6713709B2 (ja) 2020-06-24
JP6765763B2 (ja) 2020-10-07
JP6938735B2 (ja) 2021-09-22
JP2017122572A (ja) 2017-07-13

Similar Documents

Publication Publication Date Title
JP6938735B2 (ja) 居住者の不在期間の確定方法、装置、及びプログラム
Ali et al. A data-driven approach to optimize urban scale energy retrofit decisions for residential buildings
Westermann et al. Unsupervised learning of energy signatures to identify the heating system and building type using smart meter data
Loureiro et al. Water distribution systems flow monitoring and anomalous event detection: A practical approach
Berges et al. Enhancing electricity audits in residential buildings with nonintrusive load monitoring
US20150377935A1 (en) Signal identification methods and systems
EP2595098B1 (en) Method and system for detecting an appliance based on users&#39; feedback information
US20120290230A1 (en) Methods and Apparatuses for Monitoring Energy Consumption and Related Operations
EP3262576A1 (en) Method and system of monitoring appliance usage
WO2017064486A1 (en) System for energy consumption prediction
WO2013090026A2 (en) Automated monitoring for changes in energy consumption patterns
Bacher et al. Non-parametric method for separating domestic hot water heating spikes and space heating
KR20150128163A (ko) 건물 에너지 소비정보 분석 시스템
Kavousian et al. Data-driven benchmarking of building energy efficiency utilizing statistical frontier models
WO2017066435A1 (en) Sensor data analytics and alarm management
JPWO2014126131A1 (ja) 生活行動推定システム、生活行動推定装置、生活行動推定プログラム、及び記録媒体
Zaeri et al. Disaggregation of electricity and heating consumption in commercial buildings with building automation system data
CN111066003A (zh) 对流程行业中的设备进行状态监测的方法、装置及介质
Panagopoulos et al. AdaHeat: A general adaptive intelligent agent for domestic heating control
Hashim et al. False alarm moderation for performance monitoring in industrial water distribution systems
EP2733647A1 (en) On-demand power control system, on-demand power control system program, and computer-readable recording medium recorded with said program
CN107430149A (zh) 具有电源管理器的传感器集线器
JP2005514694A (ja) 取引可能な住宅排出物削減のための測定および検証プロトコル
WO2004029557A1 (en) Energy performance monitoring system
Ioannidis et al. Occupancy-driven facility management and building performance analysis

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200825

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200825

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20210831

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210901

R150 Certificate of patent or registration of utility model

Ref document number: 6938735

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250