JP2017167847A - 運用計画案作成装置、運用計画案作成方法、プログラムおよび運用計画案作成システム - Google Patents

Abstract

【課題】本発明の実施形態に係る運用計画案作成装置は、年月の経過に伴い性能が劣化する設備または機器に対する運用計画案を作成する。【解決手段】本発明の実施形態に係る運用計画案作成装置は、運用対象と類似するとされた計測対象である類似計測対象の計測値に基づき算出された、前記類似計測対象の性能の劣化モデルを取得する取得部と、前記類似計測対象の性能の劣化モデルと、前記運用対象に想定される利用事例とに基づき、前記運用対象の性能の劣化に関するシミュレーションを行うシミュレーション部と、前記シミュレーション結果に基づき、前記運用対象に対して行われる保全作業の実施時期を示す運用計画案を作成する運用計画案作成部と、を備える。【選択図】図１

Description

本発明の実施形態は、運用計画案作成装置、運用計画案作成方法、プログラムおよび運用計画案作成システムに関する。

近年、年月の経過に伴い性能が劣化する設備または機器に対しては、異常の早期発見などを目的として、センサ等による常時監視が行われている。これにより、オンサイトで実施していた従来の保守に比べ異常を迅速に発見し、設備等が故障する前に保全作業を行うことができる。

しかし、異常を検知する度に保全作業を行うと、突発的なコストが発生することになる。また、機器の部品交換を行った後に、別の部品が異常となり、機器全体を交換することになる事態も多々ある。このような事態を回避するためにも、設備、機器または建物全体のライフサイクルを見越した長期の運用計画を立案する必要がある。

運用計画は、設備等の耐用年数またはリース契約などの更新時期などに合わせて作成されるのが一般的である。設備等の劣化の進行具合を高精度に把握しなければ、適切な運用計画を作成することはできないが、設備等の劣化の進行具合は、使用状況、設置場所の環境などにより異なる。これから設置予定の設備等の運用計画案を作成したい場合もある。また、性能の劣化を計測項目から直接算出することができない場合もある。

特開２００９−００３５０２号公報

"Ｃｏｍｍｅｒｃｉａｌ Ｐｒｏｔｏｔｙｐｅ Ｂｕｉｌｄｉｎｇ Ｍｏｄｅｌｓ"、［ｏｎｌｉｎｅ］、［２０１６年３月１日検索］、米国エネルギー省、インターネット＜ＵＲＬ：ｈｔｔｐｓ：／／ｗｗｗ．ｅｎｅｒｇｙｃｏｄｅｓ．ｇｏｖ／ｃｏｍｍｅｒｃｉａｌ−ｐｒｏｔｏｔｙｐｅ−ｂｕｉｌｄｉｎｇ−ｍｏｄｅｌｓ）

本発明の実施形態に係る運用計画案作成装置は、年月の経過に伴い性能が劣化する設備または機器に対する運用計画案を作成する。

本発明の実施形態に係る運用計画案作成装置は、運用対象と類似するとされた計測対象である類似計測対象の計測値に基づき算出された、前記類似計測対象の性能の劣化モデルを取得する取得部と、前記類似計測対象の性能の劣化モデルと、前記運用対象に想定される利用事例とに基づき、前記運用対象の性能の劣化に関するシミュレーションを行うシミュレーション部と、前記シミュレーション結果に基づき、前記運用対象に対して行われる保全作業の実施時期を示す運用計画案を作成する運用計画案作成部と、を備える。

第１の実施形態に係る運用計画案作成装置の概略構成の一例を示すブロック図。 劣化モデル生成部が生成した劣化モデルの一例を示す図。 劣化モデル生成処理のフローチャート。 推定手法としてパーティクルフィルタを用いる場合における劣化モデル生成部の概略構成の一例を示すブロック図。 パーティクルフィルタの処理の内容を示す図。 パーティクルフィルタによる内部パラメータの推定処理のフローチャート。 建物モデル抽出処理のフローチャート。 運用計画案の一例を示す図。 運用計画案の他の一例を示す図。 運用計画案作成処理のフローチャート 第２の実施形態に係る運用計画案作成装置の概略構成の一例を示すブロック図。 要素簡略化の一例を示す図。 直線化の一例を示す図。 分割について説明する図。 分割片の再構成について説明する図。 集約について説明する図。 空間形状加工処理のフローチャート。 空間形状加工部の概略構成の一例を示すブロック図。 方向軸を取得する方法の一例を示す図。 分割線を生成するフローチャート。 簡略区間設定の処理について説明する図。 簡略面積閾値を算出するフローチャート。 要素簡略化処理のフローチャート。 要素簡略化における凹部の簡略化について説明する図。 外周の加工処理のフローチャート。 内部の加工処理のフローチャート。 直線化処理のフローチャート。 直線化における凸部の簡略化について説明する図。 直線化における凹部の簡略化について説明する図。 凹エッジの簡略化について説明する図。 双方簡略化について説明する図。 エッジ部の簡略化のフローチャート。 凹エッジの簡略化のフローチャート。 エッジ部の整形について説明する図。 空間構造加工部の概略構成の一例を示すブロック図。 空間構造加工処理の概略フローチャート。 本実施形態に係る空間情報生成装置を実現したハードウェア構成の一例を示すブロック図。

以下、図面を参照しながら、本発明の実施形態について説明する。

（第１の実施形態）
図１は、第１の実施形態に係る運用計画案作成装置の概略構成の一例を示すブロック図である。第１の実施形態に係る運用計画案作成装置は、運用計画案作成処理部１と、劣化モデル処理部２と、建物モデル処理部３とを備える。

運用計画案作成処理部１は、入力部１１と、取得部１２と、運用計画案作成部１３と、シミュレーション部１４と、運用計画案記憶部１５と、出力部１６と、を備える。

劣化モデル処理部２は、計測データ（センサデータ）管理部２１と、オントロジー管理部２２と、劣化モデル管理部２３とを備える。計測データ管理部２１は、計測データ取得部２１１と、計測データ記憶部２１２とを備える。オントロジー管理部２２は、オントロジー記憶部２２１と、特徴量データ抽出部（利用事例抽出部）２２２と、オントロジーデータ記憶部２２３とを備える。劣化モデル管理部２３は、劣化モデル生成部（パラメータキャリブレータ部）２３１と、オントロジー取得部２３２と、劣化モデル記憶部２３３とを備える。

建物モデル処理部３は、建物データ記憶部３１と、建物モデル抽出部３２と、抽出結果記憶部３３とを備える。

第１の実施形態に係る運用計画案作成装置は、運用対象の運用計画案を作成する。運用対象は、設備または機器（設備等）などであり、経年劣化により性能が劣化するものであればよい。例えば、空調機器、電源機器などが運用対象となり得る。また、運用対象の劣化は、運用対象の使われ方、設置場所の環境などに依存するものとする。

なお、運用対象の使われ方、設置場所の環境など、運用対象の劣化を引き起こす要因となるものを、利用事例と称することとする。

運用計画案とは、運用対象に対して行われる保全作業の実施時期を示すものである。保全作業には、設備等の一部または全部の交換、点検、清掃、および補修、新しいタイプの機器へのリプレースなどの作業が含まれる。なお、設備等ごとに運用計画案を作成するだけでなく、複数の設備等が設置される建物全体の運用計画案を作成してもよい。

運用計画案の保全作業の実施時期は、運用対象の性能の劣化以外に基づいて決定されてもよい。例えば、運用対象にかかるコストなどに基づき決定されてもよい。

運用計画案は、運用対象の利用事例と、運用対象の劣化モデルとに基づき、作成される。劣化モデルは、運用対象などにおける性能の劣化の推移を示す。具体的には、性能に関する所定のパラメータの推移データである。

運用計画案は、さらに運用対象の建物モデルに基づいて作成されてもよい。建物モデルは、運用対象の設置場所のモデルとして用いられる。また、運用対象が空調設備などの場合に、空調設備が空調を行う対象とする空間を建物モデルとしてもよい。建物モデルの違いによっても性能の劣化の推移は異なるからである。

建物モデルは、建物または建物の構成要素の形状および構造を示す。建物の構成要素は、建物内にあるものであれば特に限られものではない。例えば、部屋、廊下、壁、階段、設備、機器などでもよい。運用対象の建物モデルは、運用対象が設置された建物または設置予定の建物の建物モデルとする。

但し、本実施形態の運用計画案は、運用対象に基づく劣化モデルと建物モデルではなく、運用対象と類似する別の対象に基づく劣化モデルと建物モデルを再利用して、運用対象の運用計画案を作成することを想定する。

この別の対象は、計測装置（センサ）などの計測対象であるとする。そして、劣化モデル処理部２が、計測装置による計測データに基づき、計測対象の劣化モデルを生成する。また、運用対象と類似するとは、運用対象と同種の設備等であって、運用対象と属性などが一致または属性の値が所定の閾値以内なものとする。また、属性が一致しなくとも、類似関係であることを示す予め定められた類似関係データに、両属性の関係が登録されている場合には、両属性は類似であるとしてもよい。運用対象の属性は特に限られるものではない。例えば、運用対象の用途、目的、利用方法、利用時間、設置建物、または設置場所などでもよい。

運用対象と類似する別の対象に基づく劣化モデルを用いることにより、センサ等が配置されていない建物の設備等または今後建築される建物に設置予定の設備等の運用計画案も作成することができる。

なお、運用対象の利用事例と建物モデルも、運用対象と類似する別の対象の利用事例と建物モデルを用いてもよい。

なお、本実施形態では、運用計画案作成処理部１と劣化モデル処理部２と建物モデル処理部３とを備える運用計画案作成装置としたが、これらの各部が個別の装置として用意され、データの授受を行うシステムとして構築されてもよい。データの授受は、有線または無線通信にて行われてもよいし、電気信号にて行われてもよい。また、劣化モデル処理部２と建物モデル処理部３がネットワーク上に存在し、クラウドサービスなどとして、劣化モデルと建物モデルを運用計画案作成処理部１に送信してもよい。

運用計画案作成処理部１と劣化モデル処理部２と建物モデル処理部３の内部構成も、個別の装置として用意されてもよい。例えば、計測データ管理部２１が独立の装置として存在し、有線または無線通信にて計測データを取得し、劣化モデル管理装置とオントロジー管理装置に計測データを送信してもよい。

まず、劣化モデル処理部２について説明する。

劣化モデル処理部２の計測データ管理部２１は、設備等の計測対象を計測することにより得られた計測データを収集し管理する。計測対象には、運用対象と同種の設備等が含まれているとする。例えば、運用対象が空調装置の場合、計測対象として空調装置が含まれているとする。運用対象と計測対象は同種であれば、メーカ、型番、設定値といった運用対象の属性は同じでもよいし、異なっていてもよい。

計測データ管理部２１の計測データ取得部２１１は、計測対象自身、計測対象を監視する計測装置（センサ）または計測装置を束ねる計測システムから、通信または電気信号などにより、計測データを収集する。本実施形態おいて、計測対象、計測装置、および計測システムは特に限られるものではない。

計測データは、計測対象または計測装置が計測できるデータであれば何でもよい。例えば、設定値、消費電力、制御信号、エラー等のログなどでもよい。例えば、計測装置が空調設備であれば、部屋の温度、湿度、熱交換器に出入りする水の流量と温度、機器の稼働音などでもよい。計測データに含まれる項目の種類は、１つでも複数でもよい。

計測データは、任意のタイミングで、計測データ取得部２１１がポーリングして取得してもよい。もしくは、運用対象、計測装置または計測システムが、任意のタイミングで、計測データ取得部２１１に送信してもよい。収集された計測データは、計測データ記憶部２１２に送られ、計測データ記憶部２１２に記憶される。

劣化モデル処理部２のオントロジー管理部２２は、オントロジーを管理する。オントロジーとは、概念同士の関係、概念と具体例との関係などを体系化したものである。オントロジーのモデルとしては、下記に説明するＲＤＦ（Ｒｅｓｏｕｒｃｅ Ｄｅｓｃｒｉｐｔｉｏｎ Ｆｒａｍｅｗｏｒｋ）などがあるが、本実施形態において特に限定されるものではない。

例えば、ＲＤＦのモデルでは、主語（ｓｕｂｊｅｃｔ）、述語（ｐｒｅｄｉｃａｔｅ）、目的語（ｏｂｊｅｃｔ）の３つの要素を用いて、リソースを表現する。主語は表現されるリソース自身であり、述語は主語の特徴または主語と目的語との関係を示す。目的語は主語との関係のある物または述語の値を示す。３つの要素の関係を関係情報（トリプル）と称する。一般に、トリプルの集合はＲＤＦグラフと称される。ＲＤＦグラフでは、主語と目的語はノードとして表され、述語はリンクとして表され、全体で１つの知識グラフとして表される。この知識グラフにて、オントロジーは概念同士の関係を表す。

オントロジー管理部２２のオントロジー記憶部２２１（知識グラフ記憶部）は、計測対象に係るオントロジーを記憶し、劣化モデル管理部２３が類似事例を検索する際に利用する。オントロジー記憶部２２１に記憶されるオントロジーは、計測データと、計測対象データ（仕様データ）と、空間データと、特徴量データと、インシデントデータとが相互に関連付けられた、ＲＤＦグラフのような知識グラフとして記憶される。

空間データは、計測対象が設置されている空間に関するデータである。例えば、空間データは、個人宅、商業ビル、工場といった設置されている建物の種類を示すデータでもよい。また、空間データは、計測対象が設置されていているフロア数、部屋番号、部屋における位置などといった設置場所も示すデータでもよい。

計測対象データ（仕様データ）は、計測対象に関するデータである。例えば、計測対象データは、設備等の種類、用途、役割、メーカ名、初期性能、使用条件、想定耐久年数などを示すデータでもよい。また、計測対象に対して行われた保全作業の内容、計測対象に生じた異常報告もしくは故障の記録、当該設備の設置場所のレイアウト変更もしくはテナントの入れ替えなどの計測対象に影響を与えるイベントなど、インシデントの記録が電子化されたデータも含まれる。

特徴量データは、計測データの特徴量を示すデータである。特徴量は、例えば、計測データの値の平均値、最大値、最小値などにしてもよい。または、例えば、所定期間には必ず特定の状態である場合または所定時刻に必ず設定が変化されるといった特徴的な状態または事象（イベント）などにしてもよい。特徴量データは、例えば、特徴量の内容、特徴量の継続時間、特徴量の抽出方法、当該抽出方法に必要な情報、特徴量を表す値などを示すデータでもよい。また、特徴量データは、計測対象の利用事例として用いられてもよい。

インシデントデータは、計測データに含まれる特定の事象（インシデント）に関するデータである。インシデントデータは、例えば、計測対象に対して行われた保全作業の内容でもよい。または計測対象に生じた異常もしくは故障の内容でもよい。または異常等を確認した報告者でもよい。または当該設備の設置場所のレイアウト変更もしくはテナントの入れ替えなど、計測対象に影響を与えるイベントでもよい。インシデントデータも特徴量データは、計測対象の利用事例として用いられてもよい。

オントロジー管理部２２の特徴量データ抽出部２２２は、計測データ記憶部２１２の計測データに基づき、特徴量データまたはインシデントデータの抽出を行う。抽出を行うための情報、例えば、計測対象、対象期間（計測日時）、特徴量の抽出方法、当該抽出方法に必要な情報は、予め与えられるものとする。

特徴量データなどを抽出する方法としては、例えば、対象期間における計測データの平均値または閾値との比較といった統計量に基づき抽出する方法がある。閾値の場合は、閾値を上回る計測データの個数または閾値を下回る計測データの個数を特徴量（頻度）として集計する。またＳＡＸ法と呼ばれる時系列データの近似表現手法を採用して、計測データを文字列表現に変換してもよい。ＳＡＸ法は、指定されたセグメント数で対象期間を分割し、各セグメント内でのデータの平均値を算出した後、指定されたアルファベット数で正規分布の各面積が均等になるように分割し、各分割区間に対して文字列（アルファベット）を割り当てる。ＳＡＸ法を用いるためのセグメント数なども与えられるものとする。

特徴量データ抽出部２２２は、抽出した特徴量データなどにより、オントロジー記憶部２２１に記憶されたオントロジー（知識グラフ）を更新する。

オントロジーにより、計測対象の種類、利用環境、設置場所、ビルの仕様などに関する抽象的な検索キーワードによっても、各データを検出することができる。例えば、「設置場所は、夏はすごく暑い」といった検索キーワードでも、計測データに基づき、オントロジーに関する他のデータの検索が可能である。また、「設置場所は西側上層階」といった検索キーワードでも、空間データに基づき、オントロジーに関する他のデータの検索が可能である。

なお、特徴量データ抽出部２２２は、オントロジー（知識グラフ）生成部として、オントロジーを生成してもよい。オントロジーの生成は、空間データ、計測対象データ、計測データ、特徴データ、インシデントデータをいずれに配置するかが定められた変換フォーマットに基づき作成すればよい。オントロジー（知識グラフ）生成部として、オントロジーを生成する場合は、オントロジーデータ記憶部２２３が、変換フォーマットと空間データと計測対象データとを予め記憶しているものとする。特徴量データ抽出部２２２は、オントロジーデータ記憶部２２３から変換フォーマットと空間データと計測対象データを取得し、計測データ記憶部２１２から計測データを取得し、計測データから特徴量データとインシデントデータを算出した上で、オントロジー（知識グラフ）を一から生成してもよい。

劣化モデル処理部２の劣化モデル管理部２３は、劣化モデルを管理する。本実施形態の劣化モデルは、計測対象の性能を示すパラメータの推移を示す推移データであり、計測対象のオントロジーと対応づけられている。これにより、特徴量データまたは抽象的な検索キーワードなどを用いて、利用事例などが運用対象と類似する計測対象の劣化モデルを検索することができる。

劣化モデル管理部２３の劣化モデル生成部（パラメータキャリブレータ部）２３１は、計測データ記憶部２１２に格納された計測データに基づき、劣化モデルを生成する。

劣化モデル生成部２３１は、所定期間における計測データに基づき、ある時刻における計測対象の所定のパラメータの値を算出する。パラメータの値の算出は周期的に行われる。このように、複数の時刻におけるパラメータの値に基づき、パラメータの推移を示すデータである劣化モデルを生成する。

なお、劣化モデル生成部２３１は、計測データに含まれる計測項目からでは直接算出することができないパラメータに対しても、パラメータの値を推定することにより、劣化モデルを生成してもよい。例えば、計測対象が空調設備の場合、空調の設定温度や部屋の温度などは、計測装置などで計測することができるが、空調設備の冷暖房効率（ＣＯＰ：Ｃｏｅｆｆｉｃｉｅｎｔ ｏｆ ｐｅｒｆｏｒｍａｎｃｅ）は、計測することができず、計測データには含まれない。このような直接計測できない内部パラメータ（非計測パラメータ）も、計測データに基づくシミュレーションなどにより推定を行う。そして、複数の時刻の推定値または確率密度分布に基づき、内部パラメータの劣化モデルが生成される。

なお、計測装置等で計測は可能であるが、実際に計測しておらず、計測データに含まれていないパラメータも内部パラメータとして推定してもよい。

推定方法は、特に限られるものではない。例えば、シミュレイティドアニーリング（Ｓｉｍｕｌａｔｅｄ Ａｎｎｅａｌｉｎｇ：ＳＡ）法などの公知の逐次最適化手法、パーティクルフィルタなどの公知の確率分布推定手法を用いてもよい。また、ビルシミュレーション部１４などの既存のシミュレーション部１４を用いてもよい。算出されたパラメータの推定値は、一意に定めてもよいし、確率密度分布にて表されてもよい。

図２は、劣化モデル生成部２３１が生成した、パラメータの推定値を確率密度分布にて表現する劣化モデルの一例を示す図である。横軸が計測対象の稼働時間を示す。縦軸は計測対象のパラメータの値を示す。劣化モデルは、図２のように内部パラメータの時系列の推移にて表される。図２には３つのグラフが示されているが、一番上の点線のグラフ（最大期待性能）は、推定された最大の性能を示す。一番下の点線のグラフ（最小期待性能）が推定された最小の性能を示す。最大期待性能と最小期待性能の間にある実線のグラフ（平均期待性能）は推定された平均の性能を示す。

図２では、時刻ｔ１とｔ２において、上記３つのグラフの上に正規分布が示されている。これは、時刻ｔ１とｔ２それぞれにおいて、劣化モデル生成部２３１が算出した、パラメータの確率密度分布である。時刻ｔ２における確率密度分布は、時刻ｔ１における確率密度分布よりも確率密度分布のばらつきは大きい。このように、推定される確率密度分布は、一般に時間の経過に伴い大きくなることが多い。

劣化モデル生成部２３１は、各時刻におけるパラメータの確率密度分布を算出して、算出した確率密度分布をつなぎ合わせることにより、推移データを生成する。劣化モデル生成部２３１が生成した推移データは劣化モデル記憶部２３３に記憶される。

劣化モデル管理部２３のオントロジー取得部２３２は、オントロジー記憶部２２１からオントロジーを取得し、劣化モデル記憶部２３３に記憶させる。なお、オントロジーではなく、オントロジー記憶部２２１に記憶されているオントロジーの位置を示す位置情報（リンク）を取得してもよい。

劣化モデル記憶部２３３は、計測対象ごとに、劣化モデル生成部２３１から送られた推移データと、オントロジー取得部２３２から送られた特徴量データを含むオントロジーを検索用のインデックスとして対応付けて記憶する。なお、オントロジーの代わりにオントロジー記憶部２２１に記憶されているオントロジーの位置を示す位置情報（リンク）を対応付けて記憶してもよい。

劣化モデル記憶部２３３は、取得部１２から検索条件を受取り、検索条件に適合した劣化モデルを抽出する。オントロジーは、推移データを抽出する際のインデックスとして用いられる。これにより、オントロジーに含まれる空間データ、計測対象データ、計測データ、特徴量データ、インシデントデータに係る検索キーワードを用いて、運用対象と類似する計測対象の劣化モデルを検索することができる。

例えば、運用計画案処理部１の取得部１２が入力部１１を介して運用対象に想定される利用事例を受け取った場合は、当該利用事例が劣化モデル記憶部２３３に渡され、劣化モデル記憶部２３３は当該利用事例と類似する利用事例を有する計測対象の劣化モデルを取得部１２に渡してもよい。また、劣化モデル記憶部２３３は、運用対象の情報または建物の利用条件などを受け取り、運用対象に類似する計測対象または当該利用条件などに適応する計測対象を検出し、当該計測対象の利用事例と劣化モデルを取得部１２に渡してもよい。

図３は、劣化モデル生成処理のフローチャートである。計測データは既に計測データ記憶部２１２に記憶されていることを想定する。

劣化モデル生成部２３１が、計測データ記憶部２１２から計測データを取得し、内部パラメータの推定処理を実行する（Ｓ１０１）。内部パラメータの推定処理のフローについては後述する。

劣化モデル生成部２３１は、算出した各時刻の内部パラメータの推定値より、推移データである劣化モデルを生成する（Ｓ１０２）。劣化モデルは既に過去に作成された同一対象のものに新たに推定した内部パラメータの推定値を追加し更新してもよい。劣化モデル生成部２３１は、登録周期経過後に、劣化モデル記憶部２３３に劣化モデルを記録する（Ｓ１０３）。

一方、特徴量データ抽出部２２２は、計測データ記憶部２１２から計測データを取得し、計測データから特徴量データを抽出する（Ｓ１０４）。特徴量データ抽出部２２２は、抽出した特徴量データにて、オントロジー記憶部２２１のオントロジーを更新する（Ｓ１０５）。

オントロジー取得部２３２は、定期的にまたはオントロジーが更新されたときに、オントロジー記憶部２２１からオントロジーまたは位置情報を取得する（Ｓ１０６）。そして、オントロジー取得部２３２は、運用ごとに、劣化モデル記憶部２３３の劣化モデルと取得したオントロジーとを対応付ける（Ｓ１０７）。以上が第１の実施形態に係る劣化モデル生成処理のフローである。

なお、このフローチャートは一例であり、これに限られるものではない。例えば、Ｓ１０４とＳ１０５の処理は、Ｓ１０１よりも前に行われても問題は生じない。このように、問題が生じなければ、順番などは入れ替えてもよい。以降に説明されるフローチャートについても同様である。

次に、劣化モデル生成部２３１が行う内部パラメータの推定について説明する。内部パラメータの推定方法としては、ベイズ推定などが用いられる。計測データに基づく計測された状態をＹ、計測されていない状態（推定状態、非計測状態）をＸとすると、状態Ｙに基づき、状態Ｘを推定することは、状態Ｙが起きた場合における状態Ｘの起きる確率（事後確率）Ｐ（Ｘ｜Ｙ）を求めることと同じである。そして、事後確率Ｐ（Ｘ｜Ｙ）はベイズの定理により次式で表される。

ベイズ推定では、上記数式において、Ｘを確率変数とし、Ｘを確率密度関数Ｐにおけるパラメータとみなす。以降、Ｘを推定パラメータと称する。そうすると、Ｐ（Ｘ）は、推定パラメータＸの事前確率密度分布、Ｐ（Ｘ｜Ｙ）は、状態Ｙが計測されたときの推定パラメータＸの事後確率密度分布となる。Ｐ（Ｙ）は状態Ｙが起きる事前確率、Ｐ（Ｙ｜Ｘ）は、パラメータＸの時にＹが得られる事後確率であり、尤度と称される。

さらに、時刻ｔ（ｔは正の実数）における推定パラメータをＸｔ、数式１は次式に置き換えることができる。
Ｙ１：ｔとは、時刻ｔまでに計測されたデータの集合Ｙ＝｛Ｙ１、Ｙ２、・・・Ｙｔ｝を意味する。つまり、Ｐ（Ｘｔ│Ｙ１：ｔ）は、計測開始時刻から現在時刻までの計測値に基づく、推定パラメータＸの確率密度分布を意味する。

なお、確率密度分布の分布形状に着目する場合、Ｐ（Ｙｔ│Ｙ１：ｔ−１）は、Ｘに依存しない定数であるため、無視してもよい。よって、次式で表される。

上記数式３によれは、新たに計測値Ｙｔを得て、尤度Ｐ（Ｙｔ│Ｘｔ）を求めることにより、先の時刻ｔ−１までの計測データから推定した事後確率密度分布Ｐ（Ｘｔ│Ｙ１：ｔ−１）を、現在の時刻までの計測データから推定する事後確率密度分布Ｐ（Ｘｔ│Ｙ１：ｔ）に、逐次更新できることを意味する。したがって、初期時刻ｔ＝０における適当な初期確率密度分布Ｐ（Ｘ０）から始めて、尤度の算出と事後確率密度分布の更新を繰り返すことで、現在時刻の推定パラメータＸの確率密度分布を求めることができる。

このように、事後確率密度分布を求める方法としては、ギブス法、メトロポリス法などを含むマルコフ連鎖モンテカルロ法（ＭＣＭＣ：Ｍａｒｋｏｖ ｃｈａｉｎ Ｍｏｎｔｅ Ｃａｒｌｏ ｍｅｔｈｏｄｓ）、逐次モンテカルロ法の一種であるパーティクル法（パーティクルフィルタ）などを用いればよい。

劣化モデル生成部２３１は、予め定められた上記の手法を用いて、事後確率密度分布を算出する。なお、尤度Ｐ（Ｙｔ│Ｘｔ）はシミュレーションにて求めればよい。シミュレーションを利用する際は、運用計画案作成処理部１のシミュレーション部１４を用いるが、劣化モデル生成部２３１自体がシミュレーション部を含む場合もあり得る。

劣化モデル生成部２３１が、事後確率密度分布を推定する一例として、パーティクルフィルタを推定手法として用いる場合を次に説明する。

パーティクルフィルタは、推定パラメータＸの事後確率密度分布Ｐ（Ｘ｜Ｙ）を、多数のパーティクルを有するパーティクル群の分布で近似する手法である。パーティクルフィルタは、予測、尤度計算、リサンプリング（パーティクル群の分布の更新）を逐次的に繰り返すことによって、現在時刻における推定パラメータＸの事後確率密度分布を算出する。

パーティクルの数は、一般的に１００から１万個程度の範囲で任意に定めるとする。パーティクルの総数が多くなれば、推定精度が向上するが、推定計算に要する時間が長くなる。なお、パーティクル群は、パーティクルの数をｎ（ｎは正の整数）個とすると、Ｐ＝｛ｐ１、ｐ２、・・、ｐｉ・・・ｐｎ｝で表される。ｉは１以上、ｎ以下の整数である。

なお、推定する状態が複数ある場合、推定パラメータＸはｍ（ｍは正の整数）個の構成要素を含むｎ次元ベクトルＸ＝｛ｘ１、ｘ２、・・・ｘｍ｝で表すことができる。例えば、ＣＯＰと、１人あたりの想定発熱量の２つを推定したい場合は、ｘ１をＣＯＰ、ｘ２を１人あたりの想定発熱量とするが、他の情報も含んでいる場合もある。それぞれのパーティクルは、前記の計測値Ｙｔとパーティクルの各構成要素を入力として、乱数と予め定められたモデル式（状態方程式）を用いて、時刻ｔ＋１における各パーティクルの構成要素の予測値と計測予測値Ｙｔ＋１を算出することが可能となる全ての情報を含んでいる。この場合、ｉ番目のパーティクルは、次式で表される。
ｐ_ｉ＝｛ｘ１_ｉ、ｘ２_ｉ、・・・・・，ｘｍ_ｉ、重みｉ｝
重みｉは、後述するリサンプリングの処理にて用いられる数値である。パーティクルの各要素の値および重みは、浮動小数点または整数で表される。

図４は、推定手法としてパーティクルフィルタを用いる場合における劣化モデル生成部２３１の概略構成の一例を示すブロック図である。この場合の劣化モデル生成部２３１は、パーティクル初期設定部２３１１と、シミュレーション制御部２３１２と、パーティクルシミュレーション部２３１３と、パーティクル尤度計算部２３１４と、パーティクル変更演算部２３１５と、合成部２３１６と、を備える。

パーティクル初期設定部２３１１は、初期時刻における各パーティクルの構成要素および重みの初期値を設定する。構成要素の初期値は０、重みの初期値は１を想定しているが、他の値でもよい。

シミュレーション制御部２３１２は、各パーティクルの構成要素および重みの値をパーティクルシミュレーション部２３１３に送り、シミュレーションの実行を指示する。

パーティクルシミュレーション部２３１３は、乱数と予め定められたモデル式（状態方程式）を用いて、時刻ｔ＋１における各パーティクルの構成要素の予測値を計算する。

パーティクル尤度計算部２３１４は、パーティクルシミュレーション部２３１３が算出した時刻ｔ＋１における各パーティクルの予測値と、時刻ｔ＋１における計測データの実測値の差に基づき、尤度を算出する。

尤度の算出方法は、例えば、ガウス分布に基づくノイズが観測値に入ることを仮定し、計測データの計測値と、パーティクルシミュレーション部２３１３の予測値とのユークリッド距離を正規化するなどの方法があるが、特に限定されるものではない。

パーティクル変更演算部２３１５は、パーティクル尤度計算部２３１４が算出した各パーティクルの尤度を、各パーティクルの重み値とし、リサンプリングを行う。リサンプリングとは、重み値に基づき、各パーティクルを複製または消滅させ、新たなパーティクル群を生成することを意味する。なお、消滅させたパーティクルの数だけ、パーティクルは複製されるため、パーティクルの数は一定である。

リサンプリングの方法は、パーティクルｐｉの重みｉを、全てのパーティクルの重みの総和で除算した値（重みｉ／Σ重みｉ）である選択確率Ｒｉに基づき、各パーティクルに対し、複製および消滅を行う。そして、リサンプリング終了後に存在するｎ個のパーティクルを、新しいパーティクルの集合とする。

パーティクル変更演算部２３１５は、新しいパーティクル群の全てのパーティクルの全ての構成要素の値に対し、一定の長さで予め区切られた範囲に含まれるパーティクルの構成要素の値を、当該範囲内の予め定められた値に変更する。これは、パーティクルの個数によって、確率密度分布の値を決定するためである。そして、各パーティクルの重みを１にする。このようにして、時刻ｔ＋１のパーティクル群が生成される。

図５は、パーティクルフィルタの処理の内容を示す図である。横軸は確率変数ｘ１、縦軸は確率密度を表す。

図５（Ａ）は、時刻ｔにおけるパーティクル群の分布を示す。パーティクルが別のパーティクルの上に表示されているのは、便宜上、ｘ１の値が同じパーティクルが複数あることを示す。

図５（Ｂ）は、時刻ｔ＋１におけるパーティクルの分布をシミュレーションにより予測した分布である。

図５（Ｃ）は、尤度のグラフと、パーティクルの重みを色で分類した図である。曲線で示された尤度の大きさに基づき、各パーティクルの重みが決定される。尤度の大小の判断基準は、予め定められているものとする。ここでは、尤度の小さいパーティクルを黒色に、尤度の大きいパーティクルを斜線に、それ以外のパーティクルは白色に示している。

図５（Ｄ）は、リサンプリングの結果を示す。尤度が小さい黒色のパーティクルは消滅し、尤度の大きい斜線のパーティクルは複製されている。なお、複製する個数は、重みによって異なってもよい。例えば、図５（Ｃ）の尤度が最大となるパーティクルは、図５（Ｄ）にて２つ複製されている。

図５（Ｅ）は、時刻ｔ＋１におけるパーティクル群の分布を示す。一定区間内に存在するパーティクルの値を、全て一定値にするという調整により、同じ値のパーティクルが複数存在することになり、時刻ｔ＋１における確率密度分布の形状となる。

この処理を現在の時刻まで繰り返すことで、最終的に現在時刻の事後確率密度分布が求まる。そして、事後確率密度分布の算出処理が周期的に行われることにより、事後確率密度分布の時系列のデータが求まる。

合成部は、各時刻の事後確率密度分布の値を合成して推移データとし劣化モデルを生成する。例えば、合成部は、各時刻の事後確率密度分布の平均値をつなぎ合わせて平均期待性能を生成する。

図６は、パーティクルフィルタによる内部パラメータの推定処理のフローチャートである。本フローは、パーティクルフィルタにより内部パラメータを推定する場合において、図３に示した劣化モデル生成処理のフローのＳ１０１に該当する。

パーティクル初期設定部２３１１は、確率密度分布を生成する推定パラメータに対し、以前に生成したパーティクル群があるかを確認する（Ｓ２０１）。ある場合は、Ｓ２０３の処理に移る。ない場合は、パーティクル初期設定部２３１１は、各パーティクルの初期値を決定する（Ｓ２０２）。パーティクルの数は、予め定められていることを想定しているが、この時にパーティクル初期設定部２３１１が決定してもよい。

シミュレーション制御部２３１２は、全パーティクルの構成要素の値をシミュレーション部１４２３１５に送る（Ｓ２０３）。パーティクルシミュレーション部２３１３は、取得した全パーティクルに対し、シミュレーションを行い、次の時刻における各パーティクルの予測値を算出する（Ｓ２０４）。

パーティクル尤度計算部２３１４は、シミュレーション制御部２３１２から予測値を、計測データ記憶部２１２１０３から計測データを取得し、予測値と計測データに基づき、各パーティクルの尤度を計算する（Ｓ２０５）。

パーティクル尤度計算部２２１４は、リサプリングと各パーティクルの値の調整を行い、新たなパーティクル群を生成する（Ｓ２０６）。生成した新たなパーティクル群が現在時刻のパーティクル群かを確認し（Ｓ２０７）、現在時刻でのパーティクル群でない場合は（Ｓ２０７のＮＯ）、Ｓ２０３の処理に戻る。現時時刻でのパーティクル群である場合は（Ｓ２０７のＹＥＳ）、当該処理は終了し、確率密度分布が内部パラメータの推定値（範囲）となる。

次に建物モデル処理部３について説明する。建物モデル処理部３は、建物モデルを含む、建物に関する様々なデータ（建物データ）を管理する。そして、所定の情報に基づき、建物モデルを抽出および加工することにより、設備運用計画の作成に用いられる建物モデルを生成する。

建物データ管理部の建物データ記憶部３１は、様々な建物の建物データが予め記憶されている。記憶される建物データは、例えば、ＢＩＭモデル（Ｂｕｉｌｄｉｎｇ Ｉｎｆｏｒｍａｔｉｏｎ Ｍｏｄｅｌ）などのＣＡＤデータなどがある。

ＢＩＭモデルなどの建物データには、オブジェクト、そのオブジェクトの属性に関する属性情報（建物属性）、他のオブジェクトとの関係性を表す関係情報などが含まれる。オブジェクトには、建物を構成する、空間、部材（構成物）、設備、機器などを表すオブジェクトなどがある。また、これらのオブジェクトには、頂点の位置座標など形状に関する情報が含まれる。また空間は、床、壁、天井、仮想の区切りなどにより囲まれた空間（室）を表す。扉などで仕切られておらず、空間の境界となる建物部材がない場合でも、仮想の区切りがあるものとしてよい。空間は、平面も立体も含むものとする。建物の一部または構成物は、例えば、窓、柱、階段といったものがある。設備等または機器は、空調、照明、センサ、無線アクセスポイントなど、建物内に存在している設備等であればよい。

属性情報には、例えば、そのオブジェクトの名称、面積、体積、材料、材質、性能、用途、状態、存在する階（フロア）などがある。関係情報には、構造関係、構成関係、および接続関係などがある。

なお、建物データには、加工処理に用いられる情報が含まれていればよく、加工処理に用いられない情報は、含まれていなくともよい。例えば、加工処理に材料の属性が不要であれば、材料の属性の値が空でもよい。また、建物データは、ＢＩＭソフトウェアにより生成されたものでもよいし、空間情報生成装置のために加工または新規作成されたものでもよい。また、ここでは、ＢＩＭモデルを加工することを想定して説明を行うが、ＢＩＭモデルに限らず、必要な情報を含む建物データであればよい。

建物データ管理部の建物モデル抽出部３２は、運用対象が設置される建物と類似する建物を類似建物と判断する。そして、類似建物に係る建物モデルを、運用対象の建物モデルとして、建物データ記憶部３１から取得する。なお、建物モデル抽出部３２が抽出した結果は、取得部１２に渡されてもよいし、抽出結果記憶部３３に記憶されてもよい。

建物が類似するか否かの判定条件は、任意に定めてよい。例えば、建物モデル抽出部３２は、建物内のオブジェクトの属性、形状、または構造のいずれかが一致または類似しているオブジェクトを有する建物を類似建物と判定する。

例えば、建物データが有する属性を比較し、両属性が一致するかを確かめればよい。両属性が一致しなくとも、類似関係であることを示す予め定められた類似関係データに、両属性の関係が登録されている場合には、両属性は類似であると判定してもよい。または、両属性が値で表されている場合に、両属性の値の差分が閾値以下の場合に、両属性は類似であると判定してもよい。

また、例えば、壁、底面などといった平面のオブジェクトの形状に着目し、平面のオブジェクトの一部または全部の形状が一致または相似であるときに、両形状は一致または類似であると判定してもよい。

また、窓もしくは扉などの開口部の向き、または定められた方向軸の向きなどに着目し、これらの向きが一致しているか所定範囲内であるかにより、両構造は一致または類似であると判定してもよい。

その他にも、例えば、形状に関しては、公知の形状判定方法を用いて、両建物が一致または類似であるかを判定してもよい。また、構造に関しては、ＢＩＭＱＬ（Ｂｕｉｌｄｉｎｇ Ｉｎｆｏｒｍａｔｉｏｎ Ｍｏｄｅｌ Ｑｕｅｒｙ Ｌａｎｇｕａｇｅ）などの公知のＢＩＭモデル属性検索方法を用いて、両建物が一致または類似であるかを判定してもよい。また、例えば、建物の情報を意味的な関係にて結び付けた木構造にて表し、木構造の類似度をＴＥＤ（Ｔｒｅｅ Ｅｄｉｔ Ｄｉｓｔａｎｃｅ）により算出する方法などが考えられる

図７は、建物モデル抽出処理のフローチャートである。建物モデル処理部３の建物データ記憶部３１には建物データが既に記録されていることを想定する。

建物モデル抽出部３２は、取得部１２から検索条件を取得する（Ｓ３０１）。建物モデル抽出部３２は建物データ記憶部３１を検索し、検索条件と適合する建物データを有する建物を類似建物と判定し、類似建物の建物モデルを取得する（Ｓ３０２）。建物モデル抽出部３２は、取得した建物モデルを取得部１２に渡す（Ｓ３０３）。建物モデル抽出部３２は、取得した建物モデルを抽出結果記憶部３３に記録してもよい。以上が、建物モデル抽出処理のフローである。

次に運用計画案作成処理部１について説明する。運用計画案作成処理部１は、与えられた情報に基づき、運用計画案の作成に必要な情報を劣化モデル処理部２と建物モデル処理部３から取得した上で、運用計画案を作成する。

入力部１１は、運用計画案に関する情報を受け付ける。例えば、作成される運用計画案の条件として、運用計画案の計画年数、保全作業の実施期限などがある。運用対象に契約期間があり、契約期間前に運用対象を返却しなければならないといった場合には、契約期間前に運用対象を更新する運用計画案を作成する。その他、各保全作業にかかる費用、設備等を交換する場合の新しい設備等の機種候補などがある。

入力部１１は、劣化モデルを取得するための情報を受け付ける。劣化モデルを取得するための情報としては、運用対象の利用事例、または運用対象のオントロジーを利用するためのオントロジーに含まれる空間データ、計測対象データ、計測データ、特徴量データ、インシデントデータに係る検索キーワード情報がある。運用対象の利用事例は、入力部１１から受け取るのではなく、類似する計測対象または建物の利用条件などに基づいて、オントロジー記憶部２２１または劣化モデル記憶部２３３から取得してもよい。

入力部１１は、建物モデルを取得するための情報を受け付ける。建物モデルを取得するための情報としては、例えば、建物の面積、体積、材料、材質、性能、用途、状態などの建物の属性に関する情報がある。

取得部１２は、劣化モデル記憶部２３３から劣化モデルと利用事例を取得する。利用事例に係る情報は、運用対象の使われ方を特定するものであれば特に限られるものではない。例えば、運用対象を空調機器とする場合は、日時別の空調機器のＯＮ／ＯＦＦ時間、設定温度の変化、各部屋の室温、外気温などでもよい。

運用計画案作成部１３は、運用計画案を作成する。運用計画案の基となる運用対象全体の経済性（運用コストと保全コストの和）や快適性などの性能の予測は、シミュレーション部１４が利用事例と、劣化モデルと、建物モデルとに基づいてシミュレーションを行うことにより求まる。運用計画案作成部１３は、利用事例と、劣化モデルと、建物モデルをシミュレーション部１４に設定する。そして、シミュレーションのパラメータとして、保全作業の内容、時期などを変えて、シミュレーションを行なわせる。これにより、保全作業の内容、時期などが異なるシミュレーション結果が生成される。

なお、シミュレーションに用いられる利用事例は、非特許文献１に記載した米国エネルギー省のホームページに開示されているものを例に作成してもよい。

図８は、運用計画案の一例を示す図である。図８では、運用対象に対して行われる保全作業の実施時期が、横軸（時間軸）上の△にて示されている。また、運用計画案を評価する指標として性能が縦軸に示されている。図８のように、運用計画案には、保全作業実施前後の性能の推移が示されている。これにより、保全作業の効果を見ることができる。

図８（Ａ）に示す運用計画案１（プラン１）は、更新時期が早いため、更新時期までにおける性能の劣化は小さい。図８（Ｂ）に示す運用計画案２（プラン２）は、更新時期が遅いため、期待性能のばらつきはあるが、更新時期において性能の劣化は大きい。ゆえに、プラン２では、更新時期の直前に、設備等を利用するユーザ等の不満が高まるといった可能性がある。

運用計画案作成部１３は、図８のような運用計画案を作成し、出力部１６を介して出力する。運用計画案作成部１３は、作成した運用計画案を全て出力してもよい。または、運用計画案作成部１３は、作成した運用計画案のうち、条件を満たす運用計画案または最適と判定される運用計画案を出力してもよい。

例えば、平均期待特性が閾値以下になってはならないという条件の場合において、図８に示したプラン２の平均期待特性が閾値以下であれば、プラン２は出力されなくともよい。例えば、最大期待性能が閾値以上であればよいという条件が入力された場合において、プラン２の最大期待性能が閾値以上であれば、プラン１とプラン２のいずれかまたは両方を出力してもよい。

なお、図８では、保全作業を同種の機器に更新することを想定する。ゆえに、更新後の性能の値は初期値と同じである。また、更新後の使用状況も同じとする。ゆえに、グラフの形状も更新前後において同じである。しかし、保全作業の内容は更新に限られないし、更新後の使用状況も変えてよい。

また、更新により運用対象を異なる機種に交換するとしてもよい。その場合は、運用対象のシミュレーション結果と異なる機種のシミュレーション結果をつなぎ合わせればよい。異なる機種のシミュレーションも、異なる機種の劣化モデルを運用対象と同様に取得すればよい。運用対象を異なる機種に交換した場合の運用計画案では、図８とは異なり、性能指標の値とグラフの形状が変化する。また、更新後の使用条件も変化させてよい。例えば、異なるテナントに設置され、使用条件が変わったと想定してもよい。その場合、異なるテナントに対応する建物モデルと利用事例を用いてシミュレーションを行う。

図８は評価の指標として性能を用いたが、性能以外の指標を用いてもよい。図９は、運用計画案の他の一例を示す図である。図９の例では、累積コストを評価指標とした運用計画案である。累積コストは、今までの運用対象の更新時のコストと運用コストとの総和で表される。図９（Ａ）は、図８（Ａ）で示されたプラン１を示す。図９（Ｂ）は、図８（Ｂ）で示されたプラン２を示す。

図９（Ｂ）に示されたプラン２では、更新時期の直前において、最大期待累積コストが加速度的に増加されている。これは、性能の劣化に伴い消費電力コストなどが増加することを示している。これにより、プラン２の更新後の平均期待累積コストは、プラン１の更新後の平均期待累積コストよりも大きくなる。ゆえに、例えば、図９（Ｂ）の更新時期における平均期待累積コストが最も小さい運用計画案を出力するという条件である場合は、プラン１が出力される。

運用計画案記憶部１５は、運用計画案作成部１３が作成した運用計画案を記憶する。また、ユーザなどから入力部１１を介して検索条件を受け付けて、検索条件に適合する運用計画案を、出力部１６を介して出力してもよい。

図１０は、運用計画案作成処理のフローチャートである。劣化モデルは既に生成され劣化モデル記憶部２３３に記憶されていることを想定する。また、建物モデル処理部３の建物データ記憶部３１には建物データが格納されていることを想定する。

入力部１１が入力情報を受け付ける（Ｓ４０１）。入力部１１は、必要な情報を取得部１２に渡す。取得部１２は、運用対象と運用対象が設置される建物の使用条件などの情報に基づき、劣化モデル処理部２に利用事例を要求する（Ｓ４０２）。なお、ここでは、利用事例をオントロジー記憶部２２１から取得することを想定しているが、取得部１２は、入力部１１を介してユーザまたは他のシステムなどから利用条件を取得してもよい。そのときは、Ｓ４０２の処理は省略される。劣化モデル処理部２は、取得部１２から与えられた情報に適合する利用事例をオントロジー記憶部２２１から抽出し、当該利用事例を取得部１２に渡す（Ｓ４０３）。なお、Ｓ４０２とＳ４０３の処理は、取得部１２とオントロジー記憶部２２１との間で直接行われてもよいし、オントロジー取得部２３２を介して行われてもよい。

取得部１２は、運用対象と取得した利用事例とに基づき、劣化モデル処理部２に劣化モデルを要求する（Ｓ４０４）。劣化モデル処理部２は、取得部１２から与えられた情報に適合する劣化モデルを劣化モデル記憶部２３３から抽出し、当該劣化モデルを取得部１２に渡す（Ｓ４０５）。なお、劣化モデル処理部２はＳ４０３の処理において抽出した利用事例に基づき、Ｓ４０５の処理を連続して行い、利用事例と劣化モデルを取得部１２に一度に渡してもよい。そのときは、Ｓ４０４の処理は省略される。

取得部１２は、運用対象が設置される建物の情報に基づき、建物モデル処理部３に建物モデルを要求する（Ｓ４０６）。建物モデル処理部３は、図７に示した建物モデル抽出処理を行い、建物モデルを取得部１２に渡す（Ｓ４０７）。

取得部１２は、取得した利用事例、劣化モデル、および建物モデルを運用計画案作成部１３に渡す（Ｓ４０８）。運用計画案作成部１３がシミュレーション部１４に利用事例、劣化モデル、建物モデルを設定する（Ｓ４０９）。運用計画案作成部１３が、保全作業の内容、保全作業の時期などのパラメータを変えながら、シミュレーション部１４にシミュレーションを行わせる（Ｓ４１０）。そして、運用計画案作成部１３は、取得したシミュレーション結果に基づき、運用計画案を作成する（Ｓ４１１）。

作成された運用計画案は出力部１６に渡され、出力部１６が運用計画案を出力する（Ｓ４１２）。また、作成された運用計画案は運用計画案記憶部１５に渡され、運用計画案記憶部１５に記憶されてもよい。以上が、運用計画案作成処理のフローである。

以上のように、第１の実施形態によれば、運用対象と類似する計測対象の計測データを用いて運用計画案を作成する。この際、計測データからは直接求めることができない内部パラメータも確率密度分布を用いて推定することにより、性能の劣化を予測でき、保全作業の実施時期が適切な運用計画案を作成することができる。

また、計測対象の計測データと計測対象に係るその他のデータとを体系づけたオントロジーを用いることにより、運用対象と類似する計測対象および計測対象の利用事例を簡易なキーワードでも検出することができる。

また、運用対象が設置される建物と類似する建物の建物モデルを用いることにより、運用対象が設置される建物の詳細情報がない場合または建物が建設中である場合でも、運用計画案を作成することができる。

（第２の実施形態）
第２の実施形態では、シミュレーションに用いられる建物モデルから、不要なデータを取り除くことで建物モデルを簡略化し、シミュレーションの負荷を抑える。例えば、柱等の特定の建物要素、または外気に接している壁といった特定の条件を満たす建物要素を排除してもよい。また、対象空間の外周形状を短絡または直線化してもよい。第１の実施形態と同様な点は、説明を省略する。

図１１は、第２の実施形態に係る運用計画案作成装置の概略構成の一例を示すブロック図である。第２の実施形態は、第１の実施形態に対し、建物モデル処理部３が建物モデル加工部３４をさらに備える。建物モデル加工部３４は、空間形状加工部３４１と、空間構造加工部３４２とを備える。

建物モデル加工部３４は、取得部１２から受け付けたパラメータに基づき、建物モデルを加工して簡略する。取得部１２から受け付けるパラメータとしては、加工する対象、加工する部分または範囲、加工レベル、加工方法などがある。加工レベルは、加工によって失われる面積・体積などの閾値などが考えられる。

建物モデル加工部３４の空間形状加工部３４１は、建物モデルの形状についての加工を行う。形状に関する加工とは、例えば、建物内の部屋などの外周、内周などの形状を簡略化することなどがある。例えば、建物モデルの形状の指定された要素に関する部分、あるいは指定された種類の要素の部分の形状を簡略化する。これにより、平面の当該要素に関する辺数を少なくする。

空間形状加工部３４１は、建物モデル抽出部３２または抽出結果記憶部３３から取得した建物モデルから、建物モデルの一部である平面オブジェクトを取得し、平面オブジェクトの形状を生成する。ここでは、この平面オブジェクトを加工面（基準面）と称する。

空間形状加工部３４１は、生成した加工面の形状から、指定された要素に関する部分、あるいは指定された種類の要素の部分の形状を簡略化する。これにより、加工面の当該要素に関する辺数を少なくする。この簡略化を、ここでは要素簡略化と称する。

また、空間形状加工部３４１は、加工面において、取得した建物モデルと、当該建物モデルに隣接する建物モデルとが接している隣接辺上に存在する、閾値より小さい凸部または凹部を簡略化する。この簡略化を、ここでは直線化と称する。

図１２は、要素簡略化の一例を示す図である。図１２（Ａ）は加工前の加工面を示す図である。図１２（Ｂ）には、この例の指定要素である柱に係る辺は実線で、柱以外の線は点線で示す。図１２（Ｃ）は、簡略化処理の途中を示す図である。図１２（Ｄ）は、加工後の加工面を示す。

加工前の加工面は、外周部に柱による窪み（凹部）と、内部に柱による空き空間が存在する。このような窪み、空間などは、シミュレーション部１４のシミュレーションにおいて不要とされる場合もあり得る。例えば、柱による内部の空き空間の情報は必要だが、外周部の柱による窪みは不要といった場合もあり得る。そのため、空間形状加工部３４１は、指定された省略すべき不要な情報を削除する。

空間形状加工部３４１は、指定要素の柱に関する面と、それ以外の面を区別し、柱に関する面を簡略化する。まず、外周の柱が簡略化され、図１２（Ｃ）では、外周の凹部が消滅した状態となっている。そして、内部に柱による空き空間が簡略化され、図１２（Ｄ）では、柱に関する面が全て削除されている。このようにして、空間形状加工部３４１は、加工面を簡略化する。

図１３は、直線化の一例を示す図である。ここでは、空間の外周に存在する、予め定められた閾値より小さい凸部および凹部を直線化し、オブジェクトが有する情報量を削減する。図１３（Ａ）は直線化処理前の加工面を示す図である。図１３（Ｂ）と図１３（Ｃ）は、直線化処理の途中を示すものであり、図１３（Ｂ）は、凸部および凹部を、予め定められた方法に基づき、簡略化したものである。図１３（Ｃ）は、簡略化された空間と、他の空間との重複部分を示すものである。重複部分について、さらに簡略化処理が行われる。図１３（Ｄ）は、さらなる簡略化後の加工面を示す。このようにして、空間形状加工部３４１は、加工面を直線化する。

空間形状加工部３４１は、要素簡略化と直線化のどちらか一方または両方を行うことにより、不要な情報が排除された簡略化された加工面を生成する。これにより、シミュレーションの処理の負荷を削減することができ、また計算結果の算出までの時間を短くすることができる。空間形状加工部３４１の処理の詳細については後述する。

空間構造加工部３４２は、指定された加工方法に基づき、加工面の分割または集約を行い、建物モデルを簡略化する。ここでは、分割とは、加工面を複数の分割片に分割することを意味する。また、ここでは、集約とは、複数の加工面を１つに合成することを意味する。

図１４は、分割について説明する図である。図１４（Ａ）は、これから簡略化される加工面を示す図である。図１４（Ｂ）は、当該加工面に対し、分割線を引いた図である。図１４（Ｃ）は、生成された分割片を示す図である。図１４（Ａ）に示す加工面の外周に接している黒色の四角は、外周に接する柱を示すものである。空間構造加工部３４２は、例えば柱などの構成要素を基準として分割線を生成する。そして、１つの平面を複数の分割片に分ける。

図１５は、分割片の再構成について説明する図である。図１５（Ａ）は、図１４（Ｃ）で示した図と同じであり、分割片を示す。図１５（Ｂ）は、矢印の末端の分割片が、矢印の先端の分割片に吸収されることを示す。図１５（Ｃ）は、再構成された分割片と、さらなる再構成の方向を矢印にて示す。図１５（Ｄ）は、再構成の結果を示す。分割片の再構成は、このように小さな分割片をなくす。

次に、集約について説明する。図１６は、集約について説明する図である。図１６（Ａ）の実線で囲まれた部分は加工面である。点線は分割線である。灰色で示された加工面は、分割対象に指定されていない加工面であり、白色で示された加工面は、分割対象に指定され、分割片が生成された加工面である。このように、複数の加工面がある場合において、分割対象でない加工面を対象として集約を行う。

空間構造加工部３４２は、加工面の外周の一部が隣接または共有することにより隣接関係にあると言える加工面を取得し、加工面の外周が最長となるように合成する。隣接する複数の加工面を１グループと考えれば、加工面を分割片とみなすことができる。そして、分割片の再構成と同様にすれば、集約を行うことができる。図１６（Ａ）では、白色で示された加工面の上側の３つの加工面を１グループに、白色で示された加工面の下側の２つの加工面を別の１グループとすれば、図１６（Ｂ）のように、集約される。

このようにして、分割または集約が行われることにより、建物モデルが簡略化される。空間構造加工部３４２の処理の詳細については後述する。

次に、空間形状加工処理の詳細について説明する。図１７は、空間形状加工処理のフローチャートである。空間形状加工部３４１は、全ての加工対象の建物モデルに対して処理を行う。まず空間形状加工部３４１は、加工面の形状の生成を行う（Ｓ５０１）。次に空間形状加工部３４１は、加工面の生成後、加工面の方向軸を取得する（Ｓ５０２）。加工面の方向軸は、加工を行う際の基準軸となるものである。

また、空間形状加工部３４１は、簡略区間の設定（Ｓ５０３）および簡略区間における簡略面積閾値を設定する（Ｓ５０４）。簡略区間は、加工面を形成する辺を複数の区間に分割することにより生成された、形状を簡略化する対象の区間である。簡略面積閾値は、空間形状加工部３４１の簡略化により削減される面積の上限値を示す。簡略面積閾値は、簡略化より面積が削減され過ぎるのを防ぐ。

方向軸の取得（Ｓ５０２）は、加工区間および簡略面積閾値の設定（Ｓ５０３、Ｓ５０４）と、並行して行われてもよいし、先または後に行われてもよい。方向軸の取得（Ｓ５０２）と、加工区間および簡略面積閾値の設定（Ｓ５０３、Ｓ５０４）が完了した後に、空間形状加工部３４１は、加工面の形状を簡略化する（Ｓ５０５）。簡略化は、要素簡略化と直線化のどちらか一方または両方でもよい。以上が、空間形状加工処理の概略フローチャートである。

さらに空間形状加工部３４１の詳細について説明する。図１８は、空間形状加工部３４１の概略構成の一例を示すブロック図である。空間形状加工部３４１は、加工面取得部３４１１と、方向軸取得部３４１２と、簡略区間設定部３４１３と、形状簡略部３４１４と、加工程度評価部３４１５と、加工区間情報管理部３４１６とを備える。

加工面取得部３４１１は、加工面の形状を生成する。加工面となる面は、予め定めておいてもよいし、取得部１２などから指定されてもよい。建築分野では、加工面を床面（底面）とすることが多いため、ここでは、加工面を床面として、説明する。

加工面として床面が設定されていた場合、加工面取得部３４１１は、建物モデルの属性情報と関係情報に基づき、床面を検出する。床面を検出後、予め定められた生成方法に基づき、加工面の形状を生成する。生成方法としては、例えば、床面に関する全ての要素の全ての頂点の２次元座標を取得し、各頂点間を結ぶ辺を算出し、最大の閉ループとなる形状を生成するという方法が考えられる。また、別の方法としては、空間を囲い込む側面、例えば壁に関する全ての要素の全ての頂点から、床面に関する頂点のみを抽出し、それらの２次元座標と各頂点間を結ぶ辺とに基づき、最大の閉ループとなる形状を生成する。なお、座標に誤差がある場合などは、壁同士の接続関係を考慮してもよい。

方向軸取得部３４１２は、加工面ごとに方向軸を取得する。図１９は、方向軸を取得する方法の一例を示す図である。方向軸取得部３４１２は、加工面を形成する辺のうち、方向基準として指定された要素に係る辺の向き（べクトル）を取得する。図１９では、指定要素に関する辺を実線で示されている。そして、方向軸取得部３４１２は、指定要素の辺全てにおいて、辺の向き把握した後で、直行する辺の組み合わせがあるかを確認する。直交する辺の組を発見した場合は、その辺の組を方向軸とする。直交する辺の組を複数発見した場合は、方向軸を複数としてもよいし、１つを選択してもよい。

図２０は、分割線を生成するフローチャートである。方向軸取得部３４１２は、加工面の外周を形成する辺の接続関係を取得し（Ｓ６０１）、当該接続関係に基づき、柱などの指定要素の辺が連続する区間を取得する（Ｓ６０２）。連続区間がある場合（Ｓ６０３のＹＥＳ）は、当該連続区間それぞれに対し、分割線の生成を行う。具体的には、指定要素の辺と重なる分割線を生成する（Ｓ６０４）。また、両隣も指定要素である辺を取得する（Ｓ６０５）。この辺は、凹部の窪んだ部分の辺（加工面の外周と接しない辺）を意味する。取得することができたならば（Ｓ６０６のＹＥＳ）、その辺の中点を直交する分割線を生成する（Ｓ６０７）。これにより、連続区間の分割線を生成する。

連続区間がない場合（Ｓ６０３のＮＯ）または全ての連続区間に対する分割線の生成処理（Ｓ６０７）をした後は、両隣が別要素である指定要素の辺を取得する（Ｓ６０８）。取得することができたならば（Ｓ６０９のＹＥＳ）、取得した辺それぞれに対し、辺の中点を直交する分割線を生成する（Ｓ６１０）。該当する辺がない場合（Ｓ６０９のＮＯ）または取得した辺全てに対する分割線の生成処理（Ｓ６１０）をした後は、簡略化した後の外周と直交しない分割線を取得する（Ｓ６１１）。当該分割線がない場合（Ｓ６１２のＮＯ）は処理を終了する。当該分割線を取得した場合（Ｓ６１２のＹＥＳ）は、他の分割線と直交しているかを確認し、直行していない場合（Ｓ６１３のＹＥＳ）は、分割線を削除する（Ｓ６１４）。これにより、方向軸とすることができない不要な分割線を削除することができる。全ての分割線に対し、確認および削除を行ったらば、本フローは終了する。

上記のような予め定められた方法にて方向軸が取得できない場合は、便宜的に、隣接空間の方向軸と合わせる。隣接空間の方向軸も取得できない場合は、探索の範囲を徐々に広げていき、取得可能な空間を見つける。

なお、方向軸を生成する際に、必要となる指定要素は、取得部１２などから指定されればよい。

簡略区間設定部３４１３は、加工面を形成する各辺それぞれに対し、他の空間との隣接関係に基づき、簡略区間を設定（生成）する。

図２１は、簡略区間設定の処理について説明する図である。加工対象である空間Ａが、建物外と、空間Ｂ、Ｃ、およびＤと隣接しているとする。簡略区間設定部３４１３は、対象空間Ａが別の空間と隣接する区間（辺）の両端をそれぞれ区間端に設定する。図２１では、区間端を黒の丸で示す。これにより、隣り合う空間同士の隣接辺の簡略区間が、両隣接空間同士で一致する。同じ辺であっても簡略区間の両端が異なれば、加工結果が異なる場合があり得る。したがって、これにより、各空間それぞれに対して行われた加工処理の結果が、隣接辺において整合性を保つことができる。

そして、簡略区間設定部３４１３は、隣接空間のない区間、つまり建物外に面する辺を取得し、その辺上にある頂点を取得する。そして、取得した各頂点と隣接する２つの区間端とを接続線で結び、２つの接続線が空間内にあるかを確認する。図２１では、空間内にある接続線を１点破線で表示し、空間外にはみ出してしまう接続線を破線で表示している。なお、接続線が、区間端同士を結ぶ線上にある場合も、その接続線は空間内にあるとする。頂点から出ている２つの接続線がともに空間内にある場合、その頂点を空間内頂点とする。図２１では、空間内頂点を白抜きの丸と、内部が斜線で表された丸で示す。頂点から出ている２つの接続線が一方でも空間内にない場合、その頂点を空間外頂点とする。図２１では、空間外頂点を、内部が灰色で表された丸で示す。

そして、空間内頂点のうち、空間内頂点と隣接する２つの区間端とを結ぶ線で囲まれる範囲の面積が最大となる空間内頂点を区間端に追加する。図２１では、内部が斜線で示された丸が、面積が最大となる頂点を示している。区間端に追加された頂点は、簡略化処理により削除されることがなくなる。

簡略区間設定部３４１３は、上記のように区間端を追加した後で、区間端の１つを基点として任意に選び、時計回りに外周を辿り、区間端と区間端との間の区間を簡略区間として設定する。なお、ここでは時計回りとしたが、反時計回りでもよい。なお、以降の説明において行われる処理は、時計回りを前提としており、反時計回りで設定したときは、処理の向きが逆になる。

簡略区間設定部３４１３は、簡略区間ごとに加工区間情報を生成する。加工区間情報は、簡略区間に関する情報と、当該簡略区間に行われた加工処理に関する情報が含まれる。例えば、簡略区間のＩＤ、簡略区間上に存在する頂点のＩＤと位置座標、簡略区間ごとに設定される加工面積閾値、行われた加工処理（加工ステップ）の順番を表す加工ステップ数、各加工ステップにおいて追加または削除された部位の面積、今までの加工ステップにおいて追加または削除された部位の面積の積算値、復元フラグなどが含まれることが考えられる。

復元フラグは、簡略処理によって削除された部位または区間などを、復元するかを判断するためのフラグである。復元対象となる指定要素が削除された場合に、復元フラグの値がｔｒｕｅにされればよい。指定要素は、取得部から取得すればよい。復元対象の指定要素は、前述の省略対象で指定したものの一部でも全部でもよい。

簡略区間設定部３４１３は、算出した簡略区間それぞれに対し、簡略面積閾値を設定する。図２２は、簡略面積閾値を算出するフローチャートである。簡略区間設定部３４１３は、まず加工対象の空間全体の簡略面積閾値ｄ_{ｌｉｍｉｔ} ^ｓを算出する（Ｓ７０１）。簡略面積閾値ｄ_{ｌｉｍｉｔ} ^ｓは、対象空間Ｓの面積と加工割合の積で求められる。

加工割合は、簡略対象とされる凹凸部分の元の面積に対する加除された部分の面積の比である。加工割合の値は、任意に定めてよい。

そして、簡略区間それぞれに対し、各区間の簡略面積閾値を算出する（Ｓ７０２）。とある区間ｊの簡略面積閾値ｄ_{ｌｉｍｉｔ} ^ｓｊとすると、ｄ_{ｌｉｍｉｔ} ^ｓｊは、ｄ_{ｌｉｍｉｔ} ^ｓに対し、区間ｊの長さが加工対象の空間の外周長に占める割合を積算することにより求められる。

次に、簡略区間設定部３４１３は、区間ｊを共有する隣接空間ｓｒにおける区間ｊの簡略面積閾値ｄ_{ｌｉｍｉｔ} ^ｓｒｊと、ｄ_{ｌｉｍｉｔ} ^ｓｊを絶対値で比較する（Ｓ７０３）。ｄ_{ｌｉｍｉｔ} ^ｓｊの絶対値のほうが大きい場合（Ｓ７０４のＹＥＳ）は、ｄ_{ｌｉｍｉｔ} ^ｓｊの値をｄ_{ｌｉｍｉｔ} ^ｓｒｊに置き換える。そうでない場合（Ｓ７０４のＮＯ）は、そのままにする。これにより、区間ｊを有する各空間において、区間ｊの簡略面積閾値が異なるという事態を防ぐことができる。なお、ｄ_{ｌｉｍｉｔ} ^ｓｒｊがまだ算出されていない場合は、ｄ_{ｌｉｍｉｔ} ^ｓｒｊの値を非常に大きな値にして比較してもよいし、比較を省略してもよい。そして、当該簡略区間の加工区間情報の加工面積閾値を更新し（Ｓ７０６）、次の区間の処理に移る。全ての簡略区間で処理が終了すると、本フローは終了する。なお、ここでは、絶対値により比較を行ったが、面積の増減量に対する、負の値から正の値までの許容範囲を定めてもよい。

なお、加工区間情報には、加工ステップごとに、当該加工ステップ時における簡略区間の情報が含まれる。ゆえに、加工区間情報を参照することにより、最後の加工処理後の簡略区間の状態のみならず、各加工ステップにおける状態も参照することができる。

また、簡略区間設定部３４１３は、簡略すべき指定要素が指定されたときは、当該指定要素にかかる面（辺）の形状の一部または全部を、簡略区間として設定してもよい。

形状簡略部３４１４は、対象の加工面に対し、要素簡略化または直線化を行う。要素簡略化および直線化は、いずれか一方のみ行われてもよいし、両方行われてもよい。いずれの処理または両方の処理を行うか否かは、予め定めておいてもよいし、判断基準を定めておいてもよい。判断基準は、例えば、指定された要素の種類、または簡略対象の面積などにすればよい。

要素簡略化の詳細について説明する。図２３は、要素簡略化処理のフローチャートである。形状簡略部３４１４は、外周の加工（Ｓ８０１）または内部の加工（Ｓ８０２）またはその両方を行う。外周の加工と内部の加工については後述する。上記片方または両方の処理を行った後は、これらの処理により削除された指定要素を後で復元するか否かで処理が異なる。

指定要素を後で復元する場合（Ｓ８０３のＹＥＳ）は、指定部位単位で復元するか否かを確認する。指定部位単位で復元する場合（Ｓ８０４のＹＥＳ）は、加工区間情報ごとに復元する指定部位が加工区間情報に含まれているかを確認する。指定部位が加工区間情報に含まれていた場合（Ｓ８０５のＹＥＳ）は、当該部位の復元フラグをｔｒｕｅにする（Ｓ８０６）。これにより、指定された特定の部位だけを復元させることができる。全ての加工区間情報に対し処理を行った場合は、処理を終了する。

指定要素を後で復元しない場合（Ｓ８０３のＮＯ）は、加工した全区間の加工区間情報の変化した面積を積算してｄ_{ｅｌｅｍｅｎｔ} ^ｓを算出する（Ｓ８０７）。算出したｄ_{ｅｌｅｍｅｎｔ} ^ｓの絶対値が上限値を超えた場合（Ｓ８０８のＹＥＳ）には、元に戻す必要があるため、加工した全区間の加工区間情報の復元フラグをｔｒｕｅにし（Ｓ８０９）、処理を終了する。これにより、指定要素の全部位を復元させる。算出したｄ_{ｅｌｅｍｅｎｔ} ^ｓの絶対値が上限値を超えていない場合（Ｓ８０８のＹＥＳ）には、元に戻す必要はないため、処理は終了する。

指定要素を後で復元するが、指定部位単位では復元しない場合（Ｓ８０４のＮＯ）、つまり指定要素の全部位を復元する場合は、加工した全区間の加工区間情報の復元フラグをｔｒｕｅにし（Ｓ８０９）、処理を終了する。これにより、指定要素の全部位を復元させることができる。以上が、要素簡略化処理のフローチャートである。

次に、外周の加工の詳細について説明する。外周の加工は、外周に存在する指定要素に関する面を簡略化することである。簡略化の方法は、簡略化すべき面の形状に応じ、予め定めておけばよい。図２４は、要素簡略化における凹部の簡略化について説明する図である。ｃａｓｅ１から４までの４つのパタンが示されている。なお、これらのパタンは一例であり、これらのパタンに限られるものではない。

図２４（Ａ）に示すｃａｓｅ１では、省略すべき指定要素の辺（実線）と接続されている２辺（点線）を、２辺の交点まで延長させることにより、凹部を簡略化するパタンである。図２４（Ｂ）に示すｃａｓｅ２では、前述の２辺が平行な場合に、省略すべき指定要素の辺と２辺との各接点から等距離にある前述の２辺の垂線と、前述の２辺の延長線とにより、凹部を簡略化するパタンある。図２４（Ｃ）に示すｃａｓｅ３では、前述の２辺の１つを延長した場合に、残りの１つと重なる場合に、前述の２辺の延長線により、凹部を簡略化するパタンである。図２４（Ｄ）に示すｃａｓｅ４では、前述の２辺は平行ではないが、前述の２辺の延長線が交差しない場合に、省略すべき指定要素の辺と２辺との各接点を結ぶ線により、凹部を簡略化するパタンである。

図２５は、外周の加工処理のフローチャートである。形状簡略部３４１４は、簡略区間を形成する辺の接続関係を取得する（Ｓ９０１）。また、指定要素の辺が連続する区間を取得する（Ｓ９０２）。連続する区間が取得できなかった場合（Ｓ９０３のＮＯ）は、次の簡略区間に移る。連続する区間が取得できた場合（Ｓ９０３のＹＥＳ）は、連続区間それぞれに対し、処理を行う。

まず、両端の辺それぞれと隣接する２辺を連続区間方向に延長し、その交点を取得する（Ｓ９０４）。取得できた場合（Ｓ９０５のＹＥＳ）は、連続区間の頂点を取得した交点のみとして簡略化する（Ｓ９０６）。この簡略化は図１１で示したＣａｓｅ１に該当する。

取得できなかった場合（Ｓ９０５のＮＯ）は、両辺のベクトルが同じかを確認し、同じでない場合（Ｓ９０７のＮＯ）は、連続区間の両端を接続し、他の頂点を削除して簡略化する（Ｓ９０８）。この簡略化は図１１で示したＣａｓｅ４に該当する。

両辺のベクトルが同じな場合（Ｓ９０７のＹＥＳ）は、２辺が重なるか否かを確認し、２辺が重なる場合（Ｓ９０９のＮＯ）は、連続区間の全頂点を削除し簡略化する（Ｓ９１０）。この簡略化は図１１で示したＣａｓｅ３に該当する。２辺が重ならない場合（Ｓ９０９のＹＥＳ）には、連続区間両端から等距離の地点を通る２辺と直交する線と、２辺との交点を取得し、連続区間の頂点を取得した交点のみとして簡略化する（Ｓ９１１）。この簡略化は図１１で示したＣａｓｅ２に該当する。これらにより、連続区間を４つの方法のいずれかで簡略化することができる。

上記簡略化の処理を、全ての連続区間にて行い、全ての連続区間に対する処理が完了した後は、形状簡略部３４１４は、当該簡略区間の加工区間情報を更新し（Ｓ９１２）、次の簡略区間に対する処理に移る。なお、加工区間情報の更新とは、加工区間情報を上書きするのではなく、形状簡略部３４１４が行った加工ステップにおいて、加工された結果に関する情報を追加することをいう。したがって、加工区間情報には、加工ステップの前後の情報が含まれる。全ての簡略区間に対して処理を行ったらば、本フローは終了する。

なお、簡略化する連続区間の対象を制限してもよい。例えば、連続区間の両端距離を短絡距離とし、その上限値を定める。そして短絡距離の上限値以下の連続区間を加工対象としておもよい。短絡距離の上限値は、任意に定めてよい。シミュレーション部１４の処理の負荷などに基づき、定めればよい。

次に、内部の加工の詳細について説明する。図２６は、内部の加工処理のフローチャートである。簡略区間設定部３４１３は、外周以外の辺の接続関係を取得し（Ｓ１００１）、取得した接続関係に基づき、指定要素の辺上にある連続かつ閉ループの区間を検索する（Ｓ１００２）。該当の区間が存在しない場合（Ｓ１００３のＮＯ）は、処理は終了する。該当の区間が存在した場合（Ｓ１００３のＹＥＳ）は、当該区間を簡略区間とし、加工区間情報を設定する（Ｓ１００４）。そして、形状簡略部３４１４は、当該区間を削除する（Ｓ１００５）。そして、削除した簡略区間の加工区間情報を更新する（Ｓ１００６）。連続かつ閉ループの区間が他にも存在する場合は、他の区間に対しても処理を行う。全ての連続かつ閉ループの区間に対する処理が完了したときは、本フローは終了となる。なお、簡略区間設定部３４１３と形状簡略部３４１４の処理は分けてもよい。

次に、直線化の詳細について説明する。図２７は、直線化処理のフローチャートである。当該フローは簡略区間それぞれに対して行われる。

形状簡略部３４１４は、加工区間情報の頂点ＩＤのリストから各頂点の向きを取得する（Ｓ１１０１）。頂点の向きとは、簡略区間設定部３４１３が、基点とした区間端から時計回りに外周を辿り簡略区間を設定したときに、当該頂点において、曲がった方向が、時計回りか反時計回りかを意味する。詳細は後述する。

次に、形状簡略部３４１４は、凸部優先処理および凹部優先処理を行う。凸部優先処理は、凸部の簡略化（Ｓ１１０２）、凹部の簡略化（Ｓ１１０３）、エッジ部の簡略化（Ｓ１１０４）の順に処理を行うものである。凹部優先処理は、凹部の簡略化（Ｓ１１０６）、凸部の簡略化（Ｓ１１０７）、エッジ部の簡略化（Ｓ１１０８）の順に処理を行うものである。凸部、凹部、エッジ部については後述する。それぞれの簡略化方法は同じではあるが、凸部の簡略化と凹部の簡略化のいずれかを最初に行うかで処理結果が異なる。そのため、形状簡略部３４１４は、凸部優先処理および凹部優先処理の両方を行う。凸部優先処理または凹部の簡略化の処理は、並列に行われても、別々に行われてもよく、どちらを先に行ってもよい。

凸部優先処理または凹部優先処理の後、形状簡略部３４１４は、加工区間情報に追加する情報があるかを確認する（Ｓ１１０５、Ｓ１１０９）。ある場合（Ｓ１１０５のＮＯ、Ｓ１１０９のＮＯ）は、さらに直線化を行うべき部分が残っている可能性があるため、凸部優先処理または凹部優先処理に戻る（Ｓ１１０２、Ｓ１１０６）。

凸部優先処理と凹部優先処理の両方の処理が完了した後、簡略形状を決定する（Ｓ１１１０）。簡略形状の決定とは、凸部優先処理による加工結果と凹部優先処理による加工結果を比較し、より適した加工結果である方を簡略形状として決定するものである。簡略形状の決定は、加工程度評価部３４１５が行う。詳細は、加工程度評価部３４１５にて説明する。

簡略形状が決定した後は、エッジ部の整形を行う（Ｓ１１１１）。エッジ部の整形とは、方向軸のＸ軸またはＹ軸と平行ではないエッジ部の辺を、Ｘ軸またはＹ軸に平行な線に変更することである。エッジ部の整形処理が完了したときは、次の簡略区間の処理に移る。これを繰り返して、全ての簡略区間に対し処理を行ったとき、直線化処理は終了となる。

次に、凸部および凹部の簡略化について説明する。図２８は、直線化における凸部の簡略化について説明する図である。ここでは、図２８（Ａ）に示された、空間Ａと空間Ｃとに隣接する、頂点（９）と頂点（２０）を区間端とする簡略区間を簡略化する。

凸部は、簡略区間の始端から終端までを辿るときに、簡略区間上にある頂点において、時計回り（ＣＷ：Ｃｌｏｃｋｗｉｓｅ）の向きに曲がる頂点が２つ以上連続し、かつ、反時計回り（ＣＣＷ：Ｃｏｕｎｔｅｒ Ｃｌｏｃｋｗｉｓｅ）の向きに曲がる頂点に挟まれる部分と定義する。図２８（Ｂ）に示す通り、簡略区間上には、区間端を除くと、（１０）から（１９）までの頂点が存在する。各頂点それぞれに、簡略区間の始端（９）から終端（２０）までを辿るときの当該頂点を曲がる向きの矢印が示されている。ここで、頂点（１１）の矢印の向きはＣＣＷである。頂点（１２）および（１３）の矢印の向きはＣＷである。そして、頂点（１４）の矢印の向きはＣＣＷである。ゆえに、ＣＷの向きに曲がる頂点（１２）と（１３）が連続し、かつ頂点（１２）と（１３）は、ＣＷの向きに曲がる頂点（１１）と（１４）に挟まれている。したがって、上記凸部の定義により、頂点（１１）から（１４）までの部分（図２８（Ｃ）の斜線部分）は凸部となる。このようにして、形状簡略部３４１４は、簡略区間上の凸部を認識し、簡略化処理を行う。

簡略化は、凸部の始端と終端を結ぶ線を生成し、始端と終端の間に存在する頂点を削除することとする。凸部の始端は、簡略区間の始端に最も近い頂点であり、凸部の始端は、簡略区間の終端に最も近い頂点である。先ほどの例では、頂点（１１）と（１４）が結ばれ、頂点（１２）と（１３）が削除される。これにより、図２８（Ｄ）に示す形状となる。形状簡略部３４１４は、簡略化後、再度、凸部があるかを確認する。そうすると、頂点（１０）から頂点（１６）までの部分が新たな凸部であると認識できる。先ほどと同様、凹部の始端（１０）から終端（１６）を線で結び、頂点（１１）、（１４）、（１５）を削除する。これにより、図２８（Ｅ）に示す形状となる。この形状は、頂点１８は突出しているものの、凸部の定義に合致しないため、凸部ではない。凸部がなくなったため、凸部の簡略化の処理が終了する。なお、頂点１８のような突出部分、または、逆に空間内部に切り込んだ形状である埋没部分をエッジ部と称する。

また、形状簡略部３４１４は、加工後、簡略化区間の加工区間情報を更新する。凸部を簡略した際は、簡略化した凸部の面積と、今までの簡略化処理により簡略化された凸部の総面積ｄ_凸 ^ｓｊを算出する。

図２９は、直線化における凹部の簡略化について説明する図である。図２９（Ａ）は、図２８（Ｂ）と同じである。凹部は、簡略区間の始端から終端までを辿るときに、簡略区間上にある頂点において、ＣＣＷの向きに曲がる頂点が２つ以上連続し、かつ、ＣＷの向きに曲がる頂点に挟まれる部分と定義する。ゆえに、図２９（Ｂ）、（Ｃ）、および（Ｄ）に示す灰色の部分が凹部である。凹部の簡略化は、対象が凹部なこと以外は、凸部の簡略化と同じである。形状簡略部３４１４は、簡略区間上の凹部を認識し、簡略化処理を繰り返すことで、図２９（Ｅ）に示す簡略結果を得る。図２８（Ｅ）と図２９（Ｅ）から分かるように、凸部の簡略化結果と凹部の簡略化結果は異なる。ゆえに、前述の通り、凸部の簡略化と凹部の簡略化のいずれかを最初に行うかで処理結果が異なる。

次にエッジ部の簡略化について説明する。図２８（Ｅ）のように、凸部または凹部の簡略化を行っても、突出または埋没部分であるエッジ部分が残る場合がある。このような場合に対応するため、形状簡略部３４１４は、予め定められた方法にて、エッジ部を簡略化する。

なお、ここでは、エッジ部を凹エッジと凸エッジの２つとする。凹エッジは、簡略区間の始端から終端までを辿るときに、簡略区間上にある頂点において、ＣＣＷの向きに曲がる頂点が、ＣＷの向きに曲がる頂点に挟まれる部分と定義する。凸エッジは、簡略区間上にある頂点において、ＣＷの向きに曲がる頂点が、ＣＣＷの向きに曲がる頂点に挟まれる部分と定義する。

簡略化の方法は、簡略化すべき部分の形状に応じ、予め定めておけばよい。図３０は、凹エッジの簡略化について説明する図である。ここでは、ｃａｓｅ１から４までの４つのパタンが示されている。なお、これらのパタンは一例であり、これらのパタンに限られるものではない。なお、図３０では、凹エッジが表されているが、凸エッジでもこれらのパタンは同じである。

図３０（Ａ）に示すｃａｓｅ１は、エッジ部に隣接する２辺を延長した際の交点が、当該２辺の線上に存在しないときに、２辺を交点まで延長させることにより、エッジ部を簡略化するパタンである。図３０（Ｂ）に示すｃａｓｅ２は、エッジ部に隣接する２辺を延長した際の交点が、当該２辺のいずれかの線上に存在するときに、当該２辺の一方をその交点まで延長させることにより、エッジ部を簡略化するパタンである。図３０（Ｃ）に示すｃａｓｅ３は、エッジ部に隣接する２辺を延長しても交点がない場合に、当該２辺のうちの１つを延長した線がエッジ部の辺に接触するときに、その延長した線により、エッジ部を簡略化するパタンである。図３０（Ｄ）に示すｃａｓｅ４では、エッジ部に隣接する２辺の１つを延長した場合、他の１辺と重なる場合に、その延長線により、エッジ部を簡略化するパタンである。

また、エッジ部の簡略化では、他の空間との整合性も考慮する。例えば、簡略化した形状が、他の空間との関係により、不適切な場合もあり得る。図３０（Ｅ）のｃａｓｅ０は、不適切な場合の一例である。空間Ｘと空間Ｙとの隣接辺のエッジ部をｃａｓｅ４で簡略化したものである。しかし、このように簡略化すると、空間Ｙと空間Ｚの隣接辺を分断してしまい、整合性がとれなくなる。このように、隣接辺との整合性を考慮して、簡略化したエッジ部を元に戻す場合もある。

また、隣接する空間がある場合において、片方の空間の簡略化処理結果と、他方の空間の簡略化処理結果とが必ずしも一致するとは限らない。そこで、双方簡略化を行う。図３１は、双方簡略化について説明する図である。図３１（Ａ）は、空間Ａに凸部優先処理にて簡略化が行われた結果と、空間Ｃに凹部優先処理にて簡略化が行われた結果を示す。空間Ａと空間Ｃの隣接辺には、エッジ部分がある。図３１（Ｂ）は、空間Ａおよび空間Ｃに凹エッジ簡略化処理が行われた結果を示す。凹エッジ簡略化処理のため、空間Ａ側の突出部分は削除されていない。一方、空間Ｃ側の埋没部分は削除されている。空間Ａと空間Ｃを接合すると、図３１（Ｃ）で示すように重複部分ができる。双方簡略化処理では、この重複部分を削除する。図３１（Ｄ）は、双方簡略化処理後を示す。これにより、簡略化されつつ空間の整合性が取れた形状となる。

図３２は、エッジ部の簡略化のフローチャートである。形状簡略部３４１４は、始めに凹エッジの簡略化を行う（Ｓ１２０１）。そして、隣接空間の有無を確認し、隣接空間がある場合（Ｓ１２０２のＹＥＳ）は、当該隣接空間との双方簡略化を行う。双方簡略化では、エッジ部の簡略化の前に行われた凸部と凹部の簡略化がいずれかが先に行われたかによって処理が異なる。凹部を先に簡略化していた場合（Ｓ１２０３のＮＯ）は、隣接空間は凸部を先に簡略化した結果と比較する（Ｓ１２０４）。逆に凸部を先に簡略化していた（Ｓ１２０３のＹＥＳ）は、隣接空間は凹部を先に簡略化した結果と比較する（Ｓ１２０５）。

隣接空間との比較（Ｓ１２０４，Ｓ１２０５）の結果、重複部分がない場合（Ｓ１２０６のＮＯ）は、今回の処理で簡略化された部分があるとき（Ｓ１２１０のＹＥＳ）のみ、加工区間情報を更新する（Ｓ１２１１）。

隣接空間との比較（Ｓ１２０４，Ｓ１２０５）の結果、重複部分がある場合（Ｓ１２０６のＹＥＳ）は、隣接空間を分断する簡略化結果がないか確認し、分断する簡略化結果があるとき（Ｓ１２０７のＮＯ）は、エッジの簡略化を元に戻す。隣接空間を分断する部分がないとき（Ｓ１２０７のＹＥＳ）または簡略化を元に戻した後（Ｓ１２０８）は、隣接する空間の重複部分を削除する（Ｓ１２０９）。そして、今回の処理で簡略化された部分がある場合（Ｓ１２１０のＹＥＳ）は、加工区間情報を更新する（Ｓ１２１１）。

隣接空間がない場合（Ｓ１２０２のＮＯ）は、凸エッジの簡略化を行う（Ｓ１２１２）。隣接空間がある場合は、当該隣接空間との調整により、凸エッジがなくなるため、凸エッジの簡略化を行う必要はない。しかし、隣接空間がない場合は、凸エッジの簡略化を行う必要がある。凸エッジの簡略化処理（Ｓ１２１２）後、簡略化された凹エッジまたは凸エッジがあった場合（Ｓ１２１０のＹＥＳ）は、簡略区間の加工区間情報を更新する（Ｓ１２１１）。以上が、エッジ部の簡略化のフローである。

次に、凹エッジの簡略化と凸エッジの簡略化について説明する。凹エッジの簡略化と凸エッジの簡略化は、簡略化の対象が凸部であるか凹部であるかの違いしかない。そのため、ここでは、凹エッジの簡略化について説明し、凸部簡略化については省略する。

図３３は、凹エッジの簡略化のフローチャートである。形状簡略部３４１４は、まず凹エッジを取得する（Ｓ１３０１）。凹エッジを取得できなかった場合（Ｓ１３０２のＮＯ）は、処理は終了する。凹エッジを取得できた場合（Ｓ１３０２のＹＥＳ）は、取得した凹エッジそれぞれに対し、処理を行う。

まず、凹エッジの両端の辺それぞれと隣接する２辺を連続区間方向に延長し、延長線を生成する（Ｓ１３０３）。２つの延長線の交点がある場合（Ｓ１３０４のＹＥＳ）は、当該交点が凹エッジ領域内であるかを確認する。凹エッジ領域内でない場合（Ｓ１３０５のＮＯ）は、次の凹エッジの処理に移る。凹エッジ領域内であった場合（Ｓ１３０５のＹＥＳ）は、凹エッジの頂点を取得した交点に変更して簡略化する（Ｓ１３０６）。そして、次の凹エッジの処理に移る。この簡略化は図３０で示したＣａｓｅ１に該当する。

２つの延長線の交点がない場合（Ｓ１３０４のＮＯ）は、他方の隣接辺との交点があるかを確認する。当該交点がある場合（Ｓ１３０７のＹＥＳ）は、凹エッジの頂点を取得した交点に変更して簡略化する（Ｓ１３０６）。そして、次の凹エッジの処理に移る。この簡略化は図３０で示したＣａｓｅ２に該当する。当該交点がない場合（Ｓ１３０７のＮＯ）は、凹エッジの辺と交点があることを確認する。

凹エッジの辺との交点がある場合（Ｓ１３０８のＹＥＳ）は、凹エッジの頂点を凹エッジの辺との交点に変更して簡略化して（Ｓ１３１１）、次の凹エッジの処理に移る。この簡略化は図３０で示したＣａｓｅ３に該当する。凹エッジの辺との交点がない場合は、先ほど生成した延長線同士が重なるかを確認する（Ｓ１３１０）。重なる場合（Ｓ１３１０のＹＥＳ）は、凹エッジの頂点を削除し、その延長線により凹エッジを簡略化して（Ｓ１３１１）、次の凹エッジの処理に移る。この簡略化は図３０で示したＣａｓｅ４に該当する。重ならない場合（Ｓ１３１０のＮＯ）は、簡略化せずに、次の凹エッジの処理に移る。

取得した全ての凹エッジに対する処理が完了すると、本フローは終了する。

次に、エッジ部の整形について説明する。形状簡略部３４１４は、方向軸のＸ軸またはＹ軸と平行ではないエッジ部の辺を、Ｘ軸またはＹ軸に平行な線に変更する。図３４は、エッジ部の整形について説明する図である。図３４（Ａ）は、整形前のエッジ部である。黒の丸はエッジ部の３つの頂点のうちの２つである。この２つの頂点間の辺は、方向軸のＸ軸とＹ軸いずれにも平行ではないため、形状簡略部３４１４は、この辺に対し、整形処理を行う。但し、整形処理を行うのは、対象のエッジ部の辺と接続されている２辺が、方向軸と平行である場合に限られる。なお、この方法の場合、簡略面積に変動はないため、簡略形状を決定した後でも行うことができる。

形状簡略部３４１４は、対象のエッジ部の辺と接続されている２辺がともに、方向軸のＸ軸またはＹ軸と平行である場合に、対象のエッジ部の辺の中点を通り、当該２辺の延長線との垂線を生成する。そして、当該垂線が２辺の延長線と交差する交点（図３４（Ａ）に示す白ぬきの丸）を取得する。そして、取得した２つの交点を接続した線と、各交点まで延長した２辺の延長線とにより、対象のエッジ部の辺を置き換える。図３４（Ｂ）が整形後のエッジ部である。これにより、方向軸のＸ軸またはＹ軸と平行でない加工面の形状を減らすことができる。

加工程度評価部３４１５は、簡略加工の結果が形状加工の制約範囲内であるかを判定する。具体的には、形状簡略部３４１４による直線化において、算出された凸部優先処理による加工結果と凹部優先処理による加工結果を比較し、簡略形状を決定する。但し、凸部優先処理による加工結果と凹部優先処理による加工結果が、簡略区間設定部３４１３が算出した簡略面積閾値を超えている可能性がある。そこで、加工程度評価部３４１５は、加工結果が簡略面積閾値を超えているかを確認し、超えている場合は、加工ステップを１つずつ遡り、遡った加工ステップにおける加工処理の結果が、簡略面積閾値を超えているかを確認する。これにより、加工処理の結果が簡略面積閾値未満である直近の加工ステップと、その加工ステップにおける加工結果を認識することができる。そして、簡略面積閾値未満である凸部優先処理による加工結果と、簡略面積閾値未満である凹部優先処理による加工結果を比較して、簡略形状を決定する。

加工程度評価部３４１５は、加工結果に対する評価値を算出し、当該評価値に基づき、簡略形状を決定する。評価値は、利用目的に応じて任意に定めてよい。例えば、基本軸に基づき、評価値を算出する方法が考えられる。平面の基本軸の方向（ベクトル）と簡略区間の方向（ベクトル）との差分（ずれ）を求め、評価値を差分の逆数にするなどして、差分が小さいほど、評価値が高いようにしてもよい。また、基本軸が複数ある場合は、各基本軸との簡略区間との各差分を算出し、差分の絶対値の総和が小さいほど、評価値が高くなるようにしてもよい。また、簡略化によって加除された面積が小さいほど、評価値が高いようにしてもよい。または、簡略区間に存在する頂点の数が小さいほど、評価値が高いようにしてもよい。また、評価値を算出する方法は１つでもよいし、複数の方法を組み合わせてもよい。複数の方法を組み合わせる場合は、方法ごとに重み付けを行ってもよく、重みは任意に定めてよい。

次に、空間構造加工部３４２の処理の詳細について説明する。図３５は、空間構造加工部３４２の概略構成の一例を示すブロック図である。空間構造加工部３４２は、分割片生成部３４２１と、分割片再構成部３４２２と、分割結果評価部３４２３と、分割片情報管理部３４２４とを備える。

分割片生成部３４２１は、予め指定された指定要素の種類のオブジェクトの位置を分割基準として、加工対象である加工面を分割する線を生成する。そして、分割線で囲まれた領域、または、加工面の形状の輪郭線と分割線で囲まれた領域を、分割片とする。

なお、加工面は、空間形状加工部３４１から取得してもよい。または、空間構造加工部３４２も空間形状加工部３４１の加工面取得部３４１１と同様の部を備え、加工面を生成してもよい。

分割基準となる指定要素は、壁、柱などの建物の構造に関する要素でもよいし、設備等の建物の設備に関する要素でもよい。分割基準および分割方法は、予め定められておいてもよいし、入力部１１と取得部１２を介して指定されてもよい。

分割片再構成部３４２２は、分割片を再構成する。再構成とは、複数の分割片を合成することを意味する。

分割片情報管理部３４２３は、加工された結果を、分割片情報として、管理する。分割片情報は、分割片の生成時に、分割片生成部３４２１により生成されるものである。分割片情報には、分割片に対応付けられるＩＤ、当該分割片が生成された加工ステップ数、分割片に含まれる頂点のＩＤと位置座標、分割片が合成された合成片のリストである合成片ＩＤリスト、隣接する分割片のリストである隣接片ＩＤリスト、元の空間ＩＤ、分割片の形状と重なる簡略区間を表す区間ＩＤリストなどが含まれることが考えられる。

なお、分割片情報には、加工ステップごとに、当該加工ステップ時における分割片の情報が含まれる。ゆにえ、分割片情報を参照することにより、最後の加工処理後の分割片の状態のみならず、各加工ステップにおける状態も参照することができる。

図３６は、空間構造加工処理の概略フローチャートである。空間構造加工部３４２は、初めに空間の分割に関する処理を、分割対象である加工面それぞれに対し、行う。空間の分割に関する処理は、分割線の生成（Ｓ１４０１）と、分割片の生成（Ｓ１４０２）と、分割片の再構成（Ｓ１４０３）の３つの処理からなる。分割線の生成と、分割片の生成は、分割片生成部３４２１が行う。分割片の再構成は、分割片再構成部３４２２が行う。

次に、空間構造加工部３４２は、空間の集約に関する処理を行う。集約は、分割対象以外の加工面を対象に行われる。集約対象がないまたは集約を行わないとする場合（Ｓ１４０４のＮＯ）は、集約処理は省略される。集約対象がある場合（Ｓ１４０４のＹＥＳ）は、空間構造加工部３４２は、まず集約対象であって隣接している加工面をグループ化する（Ｓ１４０５）。そして各グループそれぞれに対し、加工面を合成する（Ｓ１４０６）。これらの集約処理は、分割片再構成部３４２２が行う。

次に、先に説明して図１４を用いて、分割片を生成する方法について説明する。分割片生成部３４２１は、柱の辺と重なるような分割線を生成する。図１４（Ｂ）では、このように生成された分割線を点線で表している。また、分割片生成部３４２１は、柱の辺であって、加工面の外周と接しない辺の中点を通る垂線を生成する。図１４（Ｂ）では、この垂線は破線で表されている。また、このように生成された分割線のうち、加工面の外周とも他の分割線とも直行しない分割線は削除するものとする。そして、図１４（Ｃ）示すように、分割線で囲まれた領域、または加工面の形状の輪郭線と分割線で囲まれた領域が、分割片となる。分割片が生成されたとき、分割片生成部３４２１は、分割片情報を生成する。この分割線の生成方法は、先に説明した空間形状加工部３４１の方向軸取得部３４１２が行う方向軸の取得の方法の１つと同じである。なお、方向軸の取得の方法と異なる方法にて、分割線を生成してもよい。

次に、先に説明して図１５を用いて、分割片を再構成する方法について説明する。分割片再構成部３４２２は、図１５（Ａ）の分割片に対し、合成処理を行う。合成処理は、面積が最小の分割片を、当該分割片に対し基本軸のＸ軸またはＹ軸の方向にて隣接する分割片に合成（吸収）させるものである。図１５（Ｂ）は、Ｘ軸方向に隣接する分割片を合成する場合を示す。複数の分割片と接している場合は、合成させる分割片を任意に選んでもよいが、ここでは、面積の大きいほうに合成するものとする。この合成を、合成により新たに生成される分割片面積が予め指定された閾値を超えない限り繰り返す。これにより、一定値以上の面積を有する分割片のみが残る。次に、先ほどの合成処理とは異なる軸方向に隣接する分割片に対し、同様の合成処理を行う。図１５（Ｃ）は、Ｘ軸方向に隣接する分割片を合成した後に、Ｙ軸方向に隣接する分割片を合成する場合を示す。図１５（Ｂ）では存在した小さな分割片がなくなっていることが分る。図１５（Ｃ）では、さらに、Ｙ軸方向に分割片を合成させて、より大きい分割片を生成する。このようにして、図１５（Ｄ）のようになる。

なお、方向軸の決定方法で説明した通り、方向軸が複数ある場合には、方向軸ごとに分割片の合成を行ってもよい。

なお、Ｘ軸またはＹ軸のいずれかを先に合成するかにより合成の結果は異なる。分割片再構成部３４２２は、このため、Ｘ軸を先に行う合成と、Ｙ軸を先に行う合成の両方を行った上で、各合成結果の評価値を算出する。そして、より良い評価値の合成結果を最終結果とする。算出方法は任意に定めてよい。例えば、生成された分割片の数が少ないほうがよい場合は、分割数に基づき、評価値を算出する。また、生成された分割片の大きさが均一のほうがよい場合は、分割片の面積の標準偏差に基づき、評価値を算出する。また、生成された分割片の大きさができるだけ大きいほうがよい場合は、生成された分割片の面積と、予め定められた分割片の面積の上限値との偏差に基づき、評価値を算出する。なお、評価値を算出する方法は１つでもよいし、複数の方法を組み合わせてもよい。複数の方法を組み合わせる場合は、方法ごとに重み付けを行ってもよく、重みは任意に定めてよい。

また、分割片再構成部３４２２は、最終結果として採用した再構成による分割片に関する分割片情報と、加工区間情報を更新する。以上により、指定要素が省かれた分割片が生成される。

以上のように、第２の実施形態によれば、建物モデルの形状および構造を簡略化することができ、シミュレーション部１４の処理の負荷を減少することができる。

また、上記に説明した実施形態における各処理は、ソフトウェア（プログラム）によって実現することが可能である。よって、上記に説明した実施形態における運用計画案作成装置１は、例えば、汎用のコンピュータ装置を基本ハードウェアとして用い、コンピュータ装置に搭載されたプロセッサにプログラムを実行させることにより実現することが可能である。

図３７は、本発明の一実施形態に係る運用計画案作成装置１を実現したハードウェア構成の一例を示すブロック図である。運用計画案作成装置１は、プロセッサ４１、主記憶装置４２、補助記憶装置４３、ネットワークインタフェース４４、デバイスインタフェース４５、入力装置４６、出力装置４７を備え、これらがバス４８などを介して接続された、コンピュータ装置４として実現できる。

プロセッサ４１が、補助記憶装置４３からプログラムを読み出して、主記憶装置４２に展開して、実行することで、運用計画案作成処理部１、劣化モデル処理部２、建物モデル処理部３の機能を実現することができる。

プロセッサ４１は、コンピュータの制御装置及び演算装置を含む電子回路である。プロセッサ４１は、例えば、汎用目的プロセッサ、中央処理装置（ＣＰＵ）、マイクロプロセッサ、デジタル信号プロセッサ（ＤＳＰ）、コントローラ、マイクロコントローラ、状態マシン、特定用途向け集積回路、フィールドプログラマブルゲートアレイ（ＦＰＧＡ）、プログラム可能論理回路（ＰＬＤ）、及びこれらの組合せを用いることができる。

本実施形態の運用計画案作成装置１は、当該運用計画案作成装置１で実行されるプログラムをコンピュータ装置４に予めインストールすることで実現してもよいし、プログラムをＣＤ−ＲＯＭなどの記憶媒体に記憶して、あるいはネットワークを介して配布して、コンピュータ装置４に適宜インストールすることで実現してもよい。

ネットワークインタフェース４４は、ネットワークに接続するためのインタフェースである。ネットワークインタフェース４４は、既存の無線規格に適合したものを用いればよい。入力部１１、取得部１２、出力部１６は、このネットワークインタフェース４４にて、データの入出力を実現してもよい。ここではネットワークインタフェースを１つのみ示しているが、複数のネットワークインタフェースが搭載されていてもよい。

デバイスインタフェース４５は、外部記憶媒体５などの機器に接続するインタフェースである。外部記憶媒体５は、ＨＤＤ、ＣＤ−Ｒ、ＣＤ−ＲＷ、ＤＶＤ−ＲＡＭ、ＤＶＤ−Ｒ、ＳＡＮ（Ｓｔｏｒａｇｅ ａｒｅａ ｎｅｔｗｏｒｋ）等の任意の記録媒体でよい。各記憶部は、外部記憶媒体５としてデバイスインタフェース４５に接続されてもよい。

主記憶装置４２は、プロセッサ４１が実行する命令、および各種データ等を一時的に記憶するメモリ装置であり、ＤＲＡＭ等の揮発性メモリでも、ＭＲＡＭ等の不揮発性メモリでもよい。補助記憶装置４３は、プログラムやデータ等を永続的に記憶する記憶装置であり、例えば、ＨＤＤまたはＳＳＤ等がある。各記憶部は、主記憶装置４２、補助記憶装置４３として実現されてもよい。

また、運用計画案作成装置１の各部は、プロセッサ４１などを実装している半導体集積回路などの専用のハードウェアにて構成されてもよい。

入力装置４６は、キーボード、マウス、タッチパネル等の入力デバイスを備え、入力部１１の機能を実現する。入力装置４６からの入力デバイスの操作による操作信号はプロセッサ４１に出力される。入力装置４６または出力装置４７は、外部からデバイスインタフェース４５に接続されてもよい。

出力装置４６は、出力部１６の機能を実現する。出力装置４６は、ＬＣＤ（Ｌｉｑｕｉｄ Ｃｒｙｓｔａｌ Ｄｉｓｐｌａｙ）、ＣＲＴ（Ｃａｔｈｏｄｅ Ｒａｙ Ｔｕｂｅ）等の表示ディスプレイでもよい。

上記に、本発明の一実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。

１ 運用計画案作成処理部
１１ 入力部
１２ 取得部
１３ 運用計画案作成部
１４ シミュレーション部
１５ 運用計画案記憶部
１６ 出力部
２ 劣化モデル処理部
２１ 計測データ（センサデータ）管理部
２１１ 計測データ取得部
２１２ 計測データ記憶部
２２ オントロジー管理部
２２１ オントロジー記憶部
２２２ 特徴量データ抽出部（利用事例抽出部）
２２３ オントロジーデータ記憶部
２３ 劣化モデル管理部
２３１ 劣化モデル生成部（パラメータキャリブレータ部）
２３１１ パーティクル初期設定部
２３１２ シミュレーション部制御部
２３１３ パーティクルシミュレーション部
２３１４ パーティクル尤度計算部
２３１５ パーティクル変更演算部
２３１６ 合成部
２３２ 劣化モデル記憶部
２３３ オントロジー取得部
３ 建物モデル処理部
３１ 建物データ記憶部
３２ 建物モデル抽出部
３３ 抽出結果記憶部
３４ 建物モデル加工部
３４１ 空間形状加工部
３４１１ 基準面取得部
３４１２ 方向軸取得部
３４１３ 簡略区間設定部
３４１４ 形状簡略部
３４１５ 加工程度評価部
３４１６ 加工区間情報管理部
３４２ 空間構造加工部
３４２１ 分割片生成部
３４２２ 分割片再構成部
３４２３ 分割結果評価部
３４２４ 分割片情報管理部
４ コンピュータ装置
４１ プロセッサ
４２ 主記憶装置
４３ 補助記憶装置
４４ ネットワークインタフェース
４５ デバイスインタフェース
４６ バス
５ 外部記憶媒体

Claims (11)

1. 運用対象と類似するとされた計測対象である類似計測対象の計測値に基づき算出された、前記類似計測対象の性能の劣化モデルを取得する取得部と、
   前記類似計測対象の性能の劣化モデルと、前記運用対象に想定される利用事例とに基づき、前記運用対象の性能の劣化に関するシミュレーションを行うシミュレーション部と、
   前記シミュレーション結果に基づき、前記運用対象に対して行われる保全作業の実施時期を示す運用計画案を作成する運用計画案作成部と、
   を備える運用計画案作成装置。
2. 前記類似計測対象の性能の劣化モデルは、
   複数の時刻のそれぞれごとに、前記時刻までに計測された、前記類似計測対象の計測値に基づき算出された、前記類似計測対象の性能を表すパラメータの確率密度分布に基づき算出されたものであり
   複数の時刻のそれぞれごとに、前記時刻までに計測された、計測対象の計測値に基づき算出された、計測対象の性能を表すパラメータの確率密度分布を周期的に算出し、算出された複数の確率密度分布に基づき、前記計測対象の性能の劣化モデルを生成する劣化モデル生成部
   をさらに備える請求項１に記載の運用計画案作成装置。
3. 前記計測対象の計測データから、前記計測対象の利用事例を抽出する利用事例抽出部と、
   前記計測対象ごとに、前記計測対象の性能の劣化モデルと利用事例とを対応づけて記憶する劣化モデル記憶部と、
   をさらに備え、
   前記取得部は、前記運用対象に想定される利用事例と類似する利用事例を有する計測対象を前記類似計測対象とすることとし、前記類似計測対象の性能の劣化モデルを取得する
   請求項１または２に記載の運用計画案作成装置。
4. 前記計測データから、前記計測対象の利用事例を抽出する利用事例抽出部と、
   前記計測対象の前記特徴的な利用事例と前記計測対象に関する情報とが体系づけられたオントロジーを記憶するオントロジー記憶部と、
   前記計測対象ごとに、前記計測対象の性能の劣化モデルと前記オントロジーを対応づけて記憶する劣化モデル記憶部と
   をさらに備え、
   前記取得部は、前記運用対象に関する情報に基づき劣化モデル記憶部から選択された、前記類似計測対象の性能の劣化モデルおよび特徴的な利用事例を取得し、
   前記シミュレーション部は、前記運用対象に想定される利用事例として、前記類似計測対象の特徴的な利用事例を用いる
   請求項１から３のいずれか一項に記載の運用計画案作成装置。
5. 建物モデルを含む、建物に関するデータを記憶する建物データ記憶部と、
   前記運用対象が設置される第１建物に関するデータに基づき、前記建物データ記憶部から、前記第１建物と類似する第２建物の建物モデルを抽出する建物モデル抽出部と、
   をさらに備え、
   前記取得部は、前記第１建物の建物モデルとして、前記第２建物の建物モデルを取得し、
   前記シミュレーション部は、前記類似計測対象の性能の劣化モデルと、前記運用対象に想定される利用事例と、前記第１建物の建物モデルとに基づき、前記シミュレーションを行い、
   前記建物データ記憶部が記憶する建物に関するデータには、少なくとも建物内のオブジェクトの属性、形状、または構造に関する情報が含まれ、
   前記建物モデル抽出部は、前記１建物のオブジェクトの属性、形状、または構造の少なくともいずれかが一致または類似しているオブジェクトを有する建物を類似と判断する
   請求項１から４のいずれか一項に記載の運用計画案作成装置。
6. 少なくとも、建物モデルに含まれる平面の外周もしくは内周の形状、または指定された要素に関する部分もしくは部分の形状を、直線化または簡略化することにより、建物モデルの簡略化を行う空間形状加工部
   をさらに備える請求項５に記載の運用計画案作成装置。
7. 少なくとも、建物モデルに含まれる平面の分割または建物モデルに含まれる複数の平面の集約を行うことにより、建物モデルの簡略化を行う空間構造加工部
   をさらに備える請求項５または６に記載の運用計画案作成装置。
8. 前記取得部は、第１運用対象の劣化モデルと、第２運用対象の劣化モデルを取得し、
   前記シミュレーション部は、第１運用対象の劣化モデルに基づく第１シミュレーションと、第２運用対象の劣化モデルに基づく第２シミュレーションとを行い、
   前記運用計画案作成部は、前記第１シミュレーションの結果と前記第２シミュレーションの結果に基づき、前記第１運用対象から前記第２運用対象へ交換する場合の運用計画案を作成する
   請求項１から７のいずれか一項に記載の運用計画案作成装置。
9. 運用対象と類似するとされた計測対象である類似計測対象の計測値に基づき算出された、前記類似計測対象の性能の劣化モデルを取得する第１取得ステップと、
   前記類似計測対象の性能の劣化モデルと、前記運用対象に想定される利用事例とに基づき、前記運用対象の性能の劣化に関するシミュレーションを行うシミュレーションステップと、
   前記シミュレーション結果に基づき、前記運用対象に対して行われる保全作業の実施時期を示す運用計画案を作成する運用計画案作成ステップと、
   をコンピュータが実行する運用計画案作成方法。
10. 運用対象と類似するとされた計測対象である類似計測対象の計測値に基づき算出された、前記類似計測対象の性能の劣化モデルを取得する第１取得ステップと、
    前記類似計測対象の性能の劣化モデルと、前記運用対象に想定される利用事例とに基づき、前記運用対象の性能の劣化に関するシミュレーションを行うシミュレーションステップと、
    前記シミュレーション結果に基づき、前記運用対象に対して行われる保全作業の実施時期を示す運用計画案を作成する運用計画案作成ステップと、
    をコンピュータに実行させるためのプログラム。
11. 計測対象と、第１の通信装置と、第２の通信装置と、第３の通信装置と、を備えた運用計画案作成システムであって、
    前記第１通信装置は、
    前記計測対象の計測値を前記第２通信装置に送り、
    前記第２通信装置は、
    複数の時刻のそれぞれごとに、前記時刻までに計測された、計測対象の計測値に基づき算出された、計測対象の性能を表すパラメータの確率密度分布を周期的に算出し、算出された複数の確率密度分布に基づき、計測対象の性能の劣化モデルを生成する劣化モデル生成部を備え、
    前記第３の通信装置は、
    運用対象と類似するとされた計測対象である類似計測対象の計測値に基づき算出された、前記類似計測対象の性能の劣化モデルを取得する取得部と、
    前記類似計測対象の性能の劣化モデルと、前記運用対象に想定される利用事例とに基づき、前記運用対象の性能の劣化に関するシミュレーションを行うシミュレーション部と、
    前記シミュレーション結果に基づき、前記運用対象に対して行われる保全作業の実施時期を示す運用計画案を作成する運用計画案作成部と、
    を備える運用計画案作成システム。