JP2021004691A - Warmth realization model generation device and method - Google Patents

Warmth realization model generation device and method Download PDF

Info

Publication number
JP2021004691A
JP2021004691A JP2019118176A JP2019118176A JP2021004691A JP 2021004691 A JP2021004691 A JP 2021004691A JP 2019118176 A JP2019118176 A JP 2019118176A JP 2019118176 A JP2019118176 A JP 2019118176A JP 2021004691 A JP2021004691 A JP 2021004691A
Authority
JP
Japan
Prior art keywords
model
information
thermal
learning data
past
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
JP2019118176A
Other languages
Japanese (ja)
Other versions
JP7272877B2 (en
Inventor
眞由美 三浦
Mayumi Miura
眞由美 三浦
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.)
Azbil Corp
Original Assignee
Azbil Corp
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 Azbil Corp filed Critical Azbil Corp
Priority to JP2019118176A priority Critical patent/JP7272877B2/en
Priority to SG10202005263SA priority patent/SG10202005263SA/en
Publication of JP2021004691A publication Critical patent/JP2021004691A/en
Application granted granted Critical
Publication of JP7272877B2 publication Critical patent/JP7272877B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Air Conditioning Control Device (AREA)

Abstract

To bring an estimation result on warmth realization of a resident close to the reality.SOLUTION: A warmth realization model generation device includes: a past model generation unit 15 for generating a past model that models relations between a residence environment index of an evaluation object space and a resident's warmth realization for the evaluation object space; a season-dependent management unit 21 for obtaining a feature amount modification function indicating relations between a model feature amount that determines the shape of a function of the past model and a representative value of outdoor warmth environment based on a plurality of the past models generated over a plurality of learning data periods, and a representative value of the past outdoor warmth environment corresponding to the plurality of learning data periods; and a future model generation unit 22 for determining a future learning data period to be set for a future date in which the resident's warmth realization is desired to be estimated, calculating a model feature amount corresponding to the representative value of the outdoor warmth environment in the learning data period based on the feature amount modification function, and generating a future model modifying the latest past model based on the calculated model feature amount.SELECTED DRAWING: Figure 2

Description

本発明は、未来の日時における評価対象空間に対する居住者の温熱実感を推定するための温熱実感モデルを生成する技術に関するものである。 The present invention relates to a technique for generating a thermal sensation model for estimating a resident's thermal sensation with respect to an evaluation target space at a future date and time.

室内環境の空調制御や運用管理において、室内環境への居住者の満足度合いを精度良く把握する技術が求められている。
居住者自身からの申告や居住者の行動、画像、生理量計測等から、居住者の温熱実感情報を取得し、取得した温熱実感情報と、温熱実感情報発生時の居住者の環境に関する指標(以下、居住環境指標)との関係を学習した学習モデルを構築し、この学習済みモデルを個別温熱実感モデルとして、任意の環境に対する対象居住者の温熱実感を推定する技術がある。居住者の温熱実感情報は、暑い/寒いなどの全身温冷感の尺度や、温熱環境に起因する快適感/満足感など、である。居住環境指標としては、温度、湿度などの物理環境指標、あるいは物理環境指標を利用した指標であるPMV(Predicted Mean Vote)などがある。
In the air conditioning control and operation management of the indoor environment, there is a demand for a technology for accurately grasping the degree of satisfaction of the resident with the indoor environment.
The resident's thermal sensation information is acquired from the resident's own declaration, the resident's behavior, images, physiological amount measurement, etc., and the acquired thermal sensation information and the index related to the resident's environment when the thermal sensation information is generated ( Hereinafter, there is a technique of constructing a learning model that learns the relationship with the living environment index) and using this learned model as an individual thermal feeling model to estimate the thermal feeling of the target resident with respect to an arbitrary environment. Residents' feelings of heat information include a measure of whole body warmth / coldness such as hot / cold, and comfort / satisfaction caused by the thermal environment. Examples of the living environment index include a physical environment index such as temperature and humidity, or a PMV (Predicted Mean Vote) which is an index using the physical environment index.

特許文献1には、居住者の温熱実感情報(特許文献1では温冷感申告)と、居住環境指標(特許文献1では室内温度)との関係を実質的に学習する方法が記載されている。
また、非特許文献1には、ユーザの温熱実感情報(非特許文献1では不快感申告)と、ユーザの居住環境指標(非特許文献1ではPMV)との関係を学習する方法が記載されている。
Patent Document 1 describes a method of substantially learning the relationship between the resident's feeling of warmth (declaration of hot / cold feeling in Patent Document 1) and the living environment index (indoor temperature in Patent Document 1). ..
Further, Non-Patent Document 1 describes a method of learning the relationship between the user's thermal feeling information (declaration of discomfort in Non-Patent Document 1) and the user's living environment index (PMV in Non-Patent Document 1). There is.

居住者の温熱実感は、季節によって異なる着衣や代謝量の変化、睡眠時の自宅寝室の室温など、季節により推移する要素(以下、季節推移要素)の影響を受けて常に変化する(非特許文献2参照)。このため、過去のデータによる学習済みモデルを、学習時以降の温熱実感情報の推定に適用すると推定誤差が発生するという課題があり、推定精度の向上が求められている。 Residents' feelings of heat constantly change under the influence of seasonal factors (hereinafter referred to as seasonal factors) such as seasonal changes in clothing and metabolic rate, and room temperature in the bedroom at home during sleep (non-patent documents). 2). Therefore, there is a problem that an estimation error occurs when the trained model based on the past data is applied to the estimation of the thermal feeling information after the learning time, and improvement of the estimation accuracy is required.

特許第2713532号公報Japanese Patent No. 2713532

田中 規敏他,“既存オフィスのZEB化に関する研究 (第4報)個人の温冷感に配慮した空調制御手法の研究”,空気調和・衛生工学会大会学術講演論文集,第10巻,pp213−216,2017/9Noritoshi Tanaka et al., "Study on ZEB conversion of existing offices (4th report) Research on air-conditioning control method considering individual warmth and cold feeling", Proceedings of the Conference of the Air Conditioning and Sanitary Engineering Society, Vol. 10, pp213 216, 2017/9 久保 博子他,“温冷感と快適感の季節差に関する実験的研究”,人間と生活環境,1巻1号,pp51−57,1994Hiroko Kubo et al., "Experimental Study on Seasonal Differences in Feeling of Warm and Cold", Humans and Living Environment, Vol. 1, No. 1, pp51-57, 1994

本発明は、上記課題を解決するためになされたもので、居住者の温熱実感の推定結果を実態に近づけることができる温熱実感モデル生成装置および方法を提供することを目的とする。 The present invention has been made to solve the above problems, and an object of the present invention is to provide a thermal sensation model generation device and a method capable of bringing the estimation result of the resident's thermal sensation closer to the actual situation.

本発明の温熱実感モデル生成装置は、評価対象空間の居住環境指標と前記評価対象空間に対する居住者の温熱実感との関係をモデル化した第1のモデルの更新タイミングの情報と学習データ期間の情報とを記憶するように構成されたデータ期間管理部と、前記更新タイミングになったときに、前記学習データ期間における前記居住環境指標と前記居住者の温熱実感情報とに基づいて前記第1のモデルを生成するように構成された第1のモデル生成部と、複数の学習データ期間によってそれぞれ生成された複数の前記第1のモデルと前記複数の学習データ期間のそれぞれに対応する過去の屋外温熱環境の代表値とに基づいて、前記第1のモデルの関数の形状を決定付けるモデル特性量と屋外温熱環境の代表値との関係を示す特性量修正関数を求めるように構成された季節依存管理部と、前記居住者の温熱実感を推定したい未来の日時に対して設定されるべき未来の学習データ期間を決定して、この学習データ期間における屋外温熱環境の代表値に対応するモデル特性量を前記特性量修正関数に基づいて算出して、算出したモデル特性量に基づいて最新の前記第1のモデルを修正した第2のモデルを生成するように構成された第2のモデル生成部とを備えることを特徴とするものである。 The thermal sensation model generator of the present invention has information on the update timing and learning data period of the first model that models the relationship between the living environment index of the evaluation target space and the resident's thermal sensation with respect to the evaluation target space. The first model is based on the data period management unit configured to store the data, the living environment index in the learning data period, and the resident's thermal feeling information at the update timing. A first model generator configured to generate, and a past outdoor thermal environment corresponding to each of the plurality of first models and the plurality of training data periods generated by the plurality of training data periods. A seasonally dependent management unit configured to obtain a characteristic quantity correction function indicating the relationship between the model characteristic quantity that determines the shape of the function of the first model and the representative value of the outdoor thermal environment based on the representative value of. Then, the future learning data period to be set for the future date and time for which the resident's feeling of heat is to be estimated is determined, and the model characteristic quantity corresponding to the representative value of the outdoor thermal environment in this learning data period is calculated. It is provided with a second model generation unit configured to generate a second model obtained by calculating based on a characteristic amount correction function and modifying the latest first model based on the calculated model characteristic amount. It is characterized by that.

また、本発明の温熱実感モデル生成装置の1構成例において、前記データ期間管理部は、前記更新タイミングの情報と前記学習データ期間の情報に加えて、オーバーラップ期間の情報を記憶し、前記複数の学習データ期間は、連続する前後の学習データ期間が前記オーバーラップ期間の幅の分だけ重なるように設定されることを特徴とするものである。
また、本発明の温熱実感モデル生成装置の1構成例は、前記評価対象空間が含まれる建物の過去の屋外温熱環境情報を記憶するように構成された屋外温熱環境情報管理部をさらに備え、前記季節依存管理部は、前記複数の学習データ期間のそれぞれに対応する過去の前記屋外温熱環境情報を前記屋外温熱環境情報管理部から取得し、前記複数の学習データ期間のそれぞれに対応する過去の屋外温熱環境の代表値を決定して、前記特性量修正関数を求め、前記第2のモデル生成部は、前記未来の学習データ期間に対応する過去の前記屋外温熱環境情報を前記屋外温熱環境情報管理部から取得して、前記未来の学習データ期間に対応する屋外温熱環境の代表値を決定することを特徴とするものである。
また、本発明の温熱実感モデル生成装置の1構成例において、前記第2のモデル生成部は、前記特性量修正関数に基づいて算出したモデル特性量と最新の前記第1のモデルのモデル特性量とに基づいて、最新の前記第1のモデルを修正することを特徴とするものである。
また、本発明の温熱実感モデル生成装置の1構成例において、前記居住環境指標は、PMVであり、前記居住者の温熱実感は、温熱環境に対する不満足を申告する居住者の割合である。
Further, in one configuration example of the thermal feeling model generation device of the present invention, the data period management unit stores information on the overlap period in addition to the information on the update timing and the information on the learning data period, and the plurality of data periods are stored. The learning data period of the above is characterized in that the learning data periods before and after the continuation are set so as to overlap by the width of the overlapping period.
Further, one configuration example of the thermal experience model generation device of the present invention further includes an outdoor thermal environment information management unit configured to store past outdoor thermal environment information of the building including the evaluation target space. The seasonally dependent management unit acquires the past outdoor thermal environment information corresponding to each of the plurality of learning data periods from the outdoor thermal environment information management unit, and the past outdoor corresponding to each of the plurality of learning data periods. The representative value of the thermal environment is determined, the characteristic quantity correction function is obtained, and the second model generation unit manages the past outdoor thermal environment information corresponding to the future learning data period to the outdoor thermal environment information management. It is characterized in that the representative value of the outdoor thermal environment corresponding to the future learning data period is determined by acquiring from the department.
Further, in one configuration example of the thermal feeling model generation device of the present invention, the second model generation unit has a model characteristic amount calculated based on the characteristic amount correction function and the latest model characteristic amount of the first model. It is characterized in that the latest first model is modified based on the above.
Further, in one configuration example of the thermal sensation model generator of the present invention, the living environment index is PMV, and the thermal sensation of the resident is the ratio of the resident who declares dissatisfaction with the thermal environment.

また、本発明の温熱実感モデル生成方法は、評価対象空間の居住環境指標と前記評価対象空間に対する居住者の温熱実感との関係をモデル化した第1のモデルの更新タイミングの情報と学習データ期間の情報とを記憶するデータ期間管理部を参照し、前記更新タイミングになったときに、前記学習データ期間における前記居住環境指標と前記居住者の温熱実感情報とに基づいて前記第1のモデルを生成する第1のステップと、複数の学習データ期間によってそれぞれ生成された複数の前記第1のモデルと前記複数の学習データ期間のそれぞれに対応する過去の屋外温熱環境の代表値とに基づいて、前記第1のモデルの関数の形状を決定付けるモデル特性量と屋外温熱環境の代表値との関係を示す特性量修正関数を求める第2のステップと、前記居住者の温熱実感を推定したい未来の日時に対して設定されるべき未来の学習データ期間を決定して、この学習データ期間における屋外温熱環境の代表値に対応するモデル特性量を前記特性量修正関数に基づいて算出して、算出したモデル特性量に基づいて最新の前記第1のモデルを修正した第2のモデルを生成する第3のステップとを含むことを特徴とするものである。 Further, in the method for generating a thermal sensation model of the present invention, information on the update timing of the first model that models the relationship between the living environment index of the evaluation target space and the resident's thermal sensation with respect to the evaluation target space and the learning data period. With reference to the data period management unit that stores the information of the above, when the update timing is reached, the first model is created based on the living environment index in the learning data period and the thermal feeling information of the resident. Based on the first step to generate, the plurality of first models generated by the plurality of training data periods, and the representative values of the past outdoor thermal environment corresponding to each of the plurality of training data periods. The second step of finding the characteristic quantity correction function showing the relationship between the model characteristic quantity that determines the shape of the function of the first model and the representative value of the outdoor thermal environment, and the future in which the resident's thermal feeling is to be estimated. The future training data period to be set for the date and time was determined, and the model characteristic quantity corresponding to the representative value of the outdoor thermal environment in this training data period was calculated and calculated based on the characteristic quantity correction function. It is characterized by including a third step of generating a second model obtained by modifying the latest first model based on a model characteristic quantity.

また、本発明の温熱実感モデル生成方法の1構成例において、前記データ期間管理部は、前記更新タイミングの情報と前記学習データ期間の情報に加えて、オーバーラップ期間の情報を記憶し、前記複数の学習データ期間は、連続する前後の学習データ期間が前記オーバーラップ期間の幅の分だけ重なるように設定されることを特徴とするものである。
また、本発明の温熱実感モデル生成方法の1構成例において、前記第2のステップは、前記評価対象空間が含まれる建物の過去の屋外温熱環境情報を記憶する屋外温熱環境情報管理部を参照し、前記複数の学習データ期間のそれぞれに対応する過去の前記屋外温熱環境情報を前記屋外温熱環境情報管理部から取得し、前記複数の学習データ期間のそれぞれに対応する過去の屋外温熱環境の代表値を決定して、前記特性量修正関数を求めるステップを含み、前記第3のステップは、前記未来の学習データ期間に対応する過去の前記屋外温熱環境情報を前記屋外温熱環境情報管理部から取得して、前記未来の学習データ期間に対応する屋外温熱環境の代表値を決定するステップを含むことを特徴とするものである。
また、本発明の温熱実感モデル生成方法の1構成例において、前記第3のステップは、前記特性量修正関数に基づいて算出したモデル特性量と最新の前記第1のモデルのモデル特性量とに基づいて、最新の前記第1のモデルを修正するステップを含むことを特徴とするものである。
Further, in one configuration example of the thermal feeling model generation method of the present invention, the data period management unit stores information on the overlap period in addition to the information on the update timing and the information on the learning data period, and the plurality of data periods are stored. The learning data period of the above is characterized in that the learning data periods before and after the continuation are set so as to overlap by the width of the overlapping period.
Further, in one configuration example of the thermal experience model generation method of the present invention, the second step refers to the outdoor thermal environment information management unit that stores the past outdoor thermal environment information of the building including the evaluation target space. , The past outdoor thermal environment information corresponding to each of the plurality of learning data periods is acquired from the outdoor thermal environment information management unit, and a representative value of the past outdoor thermal environment corresponding to each of the plurality of learning data periods. Including the step of determining the characteristic quantity correction function and obtaining the characteristic quantity correction function, the third step acquires the past outdoor thermal environment information corresponding to the future learning data period from the outdoor thermal environment information management unit. The present invention is characterized by including a step of determining a representative value of the outdoor thermal environment corresponding to the future learning data period.
Further, in one configuration example of the method for generating a thermal feeling model of the present invention, the third step is divided into a model characteristic amount calculated based on the characteristic amount correction function and the latest model characteristic amount of the first model. Based on this, it is characterized by including a step of modifying the latest first model.

本発明によれば、データ期間管理部と第1のモデル生成部と季節依存管理部と第2のモデル生成部とを設け、季節推移が反映された屋外温熱環境状態に基づいて第1のモデルを修正して第2のモデルを生成することにより、季節推移要素による推定誤差を改善することができる。その結果、本発明では、居住者の温熱実感の推定結果を実態に近づけることができる。 According to the present invention, the data period management unit, the first model generation unit, the season-dependent management unit, and the second model generation unit are provided, and the first model is based on the outdoor thermal environment state in which the seasonal transition is reflected. By modifying the above to generate a second model, it is possible to improve the estimation error due to the seasonal transition element. As a result, in the present invention, the estimation result of the resident's feeling of heat can be brought closer to the actual situation.

また、本発明では、更新タイミングの情報と学習データ期間の情報に加えて、オーバーラップ期間の情報をデータ期間管理部に記憶させ、連続する複数の学習データ期間がオーバーラップ期間の幅の分だけ重なるように設定することにより、モデルの信頼性を更に向上させることができる、 Further, in the present invention, in addition to the update timing information and the learning data period information, the overlap period information is stored in the data period management unit, and a plurality of consecutive learning data periods are limited to the width of the overlap period. By setting them to overlap, the reliability of the model can be further improved.

図1は、本発明の第1の実施例に係る空調制御システムの例を示す図である。FIG. 1 is a diagram showing an example of an air conditioning control system according to a first embodiment of the present invention. 図2は、本発明の第1の実施例に係る温熱実感モデル生成装置の構成を示すブロック図である。FIG. 2 is a block diagram showing a configuration of a thermal sensation model generator according to the first embodiment of the present invention. 図3は、本発明の第1の実施例に係る温熱実感モデル生成装置の居住環境情報保持部と温熱実感情報保持部とデータ統合部の動作を説明するフローチャートである。FIG. 3 is a flowchart illustrating the operation of the living environment information holding unit, the heat sensation information holding unit, and the data integration unit of the thermal sensation model generator according to the first embodiment of the present invention. 図4は、本発明の第1の実施例における温冷感申告情報と申告日時と申告者IDの1例を示す図である。FIG. 4 is a diagram showing an example of the warm / cold feeling declaration information, the declaration date / time, and the filer ID in the first embodiment of the present invention. 図5は、本発明の第1の実施例における統合情報の1例を示す図である。FIG. 5 is a diagram showing an example of integrated information in the first embodiment of the present invention. 図6は、本発明の第1の実施例に係る温熱実感モデル生成装置の過去モデル生成部の動作を説明するフローチャートである。FIG. 6 is a flowchart illustrating the operation of the past model generation unit of the thermal experience model generation device according to the first embodiment of the present invention. 図7は、居住環境指標であるPMVと「暑い」という申告者数との関係を示す図である。FIG. 7 is a diagram showing the relationship between PMV, which is a living environment index, and the number of filers who report “hot”. 図8は、本発明の第1の実施例におけるモデリング情報の生成方法を説明するフローチャートである。FIG. 8 is a flowchart illustrating a method of generating modeling information according to the first embodiment of the present invention. 図9は、PMVとPPDとの関係を示す図である。FIG. 9 is a diagram showing the relationship between PMV and PPD. 図10は、不満足度関数の例を示す図である。FIG. 10 is a diagram showing an example of the dissatisfaction function. 図11は、本発明の第1の実施例に係る温熱実感モデル生成装置の季節依存管理部の動作を説明するフローチャートである。FIG. 11 is a flowchart illustrating the operation of the season-dependent management unit of the thermal sensation model generator according to the first embodiment of the present invention. 図12は、本発明の第1の実施例に係る温熱実感モデル生成装置の季節依存管理部の特性量修正関数決定動作を説明するフローチャートである。FIG. 12 is a flowchart illustrating the operation of determining the characteristic quantity correction function of the season-dependent management unit of the thermal feeling model generation device according to the first embodiment of the present invention. 図13は、本発明の第1の実施例におけるモデル特性量を説明する図である。FIG. 13 is a diagram illustrating a model characteristic quantity in the first embodiment of the present invention. 図14は、本発明の第1の実施例における特性量修正関数の決定方法を説明する図である。FIG. 14 is a diagram illustrating a method for determining a characteristic quantity correction function in the first embodiment of the present invention. 図15は、本発明の第1の実施例における別のモデル特性量を説明する図である。FIG. 15 is a diagram illustrating another model characteristic quantity in the first embodiment of the present invention. 図16は、本発明の第1の実施例に係る温熱実感モデル生成装置の未来モデル生成部の動作を説明するフローチャートである。FIG. 16 is a flowchart illustrating the operation of the future model generation unit of the thermal feeling model generation device according to the first embodiment of the present invention. 図17は、本発明の第1、第2の実施例の学習データ期間を説明する図である。FIG. 17 is a diagram illustrating a learning data period of the first and second embodiments of the present invention. 図18は、外気温の期間平均推移値の1例を示す図である。FIG. 18 is a diagram showing an example of a period average transition value of the outside air temperature. 図19は、本発明の第1、第2の実施例に係る温熱実感モデル生成装置を実現するコンピュータの構成例を示すブロック図である。FIG. 19 is a block diagram showing a configuration example of a computer that realizes the thermal sensation model generator according to the first and second embodiments of the present invention.

[発明の原理]
季節推移があると、学習用のデータが収集された時点で温熱実感の実態が変化していることに対応しなければならない。よって、発明者は、過去のデータによる学習済みの個別温熱実感モデル(以下、過去モデル)を、季節推移が反映された屋外温熱環境状態に基づいて修正し、未来モデルとすることに着眼した。季節推移が反映された屋外温熱環境状態で過去モデルを修正することにより、季節変動の影響をモデルに反映させることが可能である。
[Principle of invention]
When there is a seasonal transition, it is necessary to respond to the change in the actual state of thermal sensation when the learning data is collected. Therefore, the inventor focused on modifying the trained individual thermal experience model based on past data (hereinafter referred to as the past model) based on the outdoor thermal environment condition reflecting the seasonal transition to make it a future model. It is possible to reflect the influence of seasonal fluctuations in the model by modifying the past model in the outdoor thermal environment condition that reflects the seasonal transition.

季節推移が反映された屋外温熱環境状態としては、屋外の外気温などがある。例えば、過去モデル生成時の学習データ期間の外気温平均値と、居住者の温熱実感を推定したい未来の平年外気温や外気温予報値を利用して、学習済みの過去モデルから未来モデルを生成する。 Outdoor thermal environment conditions that reflect seasonal changes include outdoor outside air temperature. For example, a future model is generated from a trained past model by using the average outside temperature during the training data period when the past model is generated and the future normal outside temperature or outside temperature forecast value for which the resident's feeling of heat is to be estimated. To do.

さらに発明者は、2日間、3日間などの短期的な変動が含まれる期間のデータを使用すると過去モデルおよび未来モデルの信頼性が保てない点に着目した。季節推移という年単位の変動を考慮するためには、この変動の周期に対応する、ある程度の期間幅のデータをモデルの学習用に確保する必要がある。この期間幅としては1か月程度が適当であるが、さらに、季節推移の連続性を維持するため、信頼できる期間(2週間など)を前期間と後期間でオーバーラップさせるオーバーラップ方式とするのが有効であることに想到した。 Furthermore, the inventor noted that the reliability of past and future models cannot be maintained when using data for periods that include short-term fluctuations such as 2 days and 3 days. In order to take into account the annual fluctuation of seasonal transition, it is necessary to secure data with a certain period width corresponding to the cycle of this fluctuation for model training. A suitable period width is about one month, but in order to maintain the continuity of seasonal transitions, an overlap method is adopted in which a reliable period (such as two weeks) overlaps between the previous period and the later period. I came up with the idea that is effective.

[第1の実施例]
以下、本発明の実施例について図面を参照して説明する。本発明の第1の実施例では、過去モデルから未来モデルを生成する基本手法について説明する。
本実施例は、月単位で各月のモデルを生成する例とし、過去モデルのデータ期間月の平年外気温と、予測したい未来が含まれる未来モデルデータ期間月の平年外気温に基づき、過去モデルのモデルパラメータを特性量修正関数で修正し、特性量修正関数を利用して過去モデルを修正し、未来モデルとする。
[First Example]
Hereinafter, examples of the present invention will be described with reference to the drawings. In the first embodiment of the present invention, a basic method for generating a future model from a past model will be described.
In this example, a model for each month is generated on a monthly basis, and a past model is based on the data period month's normal outside temperature of the past model and the future model data period month's normal outside temperature including the future to be predicted. The model parameter of is modified by the characteristic quantity correction function, and the past model is modified by using the characteristic quantity correction function to make it a future model.

本実施例では、説明の簡単のために、図1に示すように、対象とする1つの空調ゾーンに居住者2名が在席する居住者申告型空調制御システムの例で説明する。居住者申告型空調制御システムについては、例えば文献「大曲 康仁他,“温冷感申告対応空調システムの実証試験”,空気調和・衛生工学会学術講演論文集,第3巻,pp41−44,2016/9」に開示されている。 In this embodiment, for the sake of simplicity, as shown in FIG. 1, an example of a resident declaration type air conditioning control system in which two residents are present in one target air conditioning zone will be described. Regarding the resident-declared air-conditioning control system, for example, the literature "Yasuhito Omagari et al.," Demonstration test of air-conditioning system for hot and cold feeling declaration ", Proceedings of the Society of Air Conditioning and Sanitary Engineering, Vol. 3, pp41-44, 2016 / 9 ”is disclosed.

図1において、100は居住者、101は評価対象空間(居住者の在席空間)、102は温冷感申告を受ける空調制御装置(コントローラ)、103は評価対象空間101の室温を計測する温度センサ、104は評価対象空間101の湿度を計測する湿度センサ、105は室内機、106は室外機である。空調制御装置102は、温度センサ103によって計測される室温が室温設定値と一致し、湿度センサ104によって計測される湿度が湿度設定値と一致するように空調機器(室内機105および室外機106)を制御する。 In FIG. 1, 100 is a resident, 101 is an evaluation target space (resident's seating space), 102 is an air conditioning controller (controller) that receives a hot / cold feeling report, and 103 is a temperature for measuring the room temperature of the evaluation target space 101. The sensor, 104 is a humidity sensor that measures the humidity of the evaluation target space 101, 105 is an indoor unit, and 106 is an outdoor unit. The air conditioner controller 102 is an air conditioner (indoor unit 105 and outdoor unit 106) so that the room temperature measured by the temperature sensor 103 matches the room temperature set value and the humidity measured by the humidity sensor 104 matches the humidity set value. To control.

図2は本実施例の温熱実感モデル生成装置の構成を示すブロック図である。温熱実感モデル生成装置は、過去モデル生成ユニット1と、未来モデル生成ユニット2とから構成される。 FIG. 2 is a block diagram showing a configuration of a thermal experience model generator of this embodiment. The thermal feeling model generation device is composed of a past model generation unit 1 and a future model generation unit 2.

過去モデル生成ユニット1は、データ収集部10と、評価対象空間101の居住環境指標と評価対象空間101に対する居住者の温熱実感との関係をモデル化した個別温熱実感モデルの更新タイミングの情報と学習データ期間の情報とを記憶するデータ期間管理部14と、更新タイミングになったときに、学習データ期間における居住環境指標と居住者の温熱実感情報とに基づいて過去モデル(第1のモデル)を生成する過去モデル生成部15(第1のモデル生成部)と、過去モデル生成部15によって生成された過去モデルを記憶する過去モデル記憶管理部16とを備えている。 The past model generation unit 1 is the data collection unit 10, and information and learning of the update timing of the individual thermal sensation model that models the relationship between the living environment index of the evaluation target space 101 and the resident's thermal sensation with respect to the evaluation target space 101. The past model (first model) is created based on the data period management unit 14 that stores the data period information, and the living environment index and the resident's thermal feeling information in the learning data period when the update timing comes. It includes a past model generation unit 15 (first model generation unit) to be generated, and a past model storage management unit 16 that stores the past model generated by the past model generation unit 15.

データ収集部10は、温熱環境情報から居住環境指標を算出する居住環境情報保持部11と、評価対象空間の居住者から申告された温熱実感情報を取得する温熱実感情報保持部12と、温熱実感情報の発生日時に対応する計測日時の温熱環境情報から算出された居住環境指標を抽出して、温熱実感情報と居住環境指標とを統合した統合情報を生成するデータ統合部13とから構成される。 The data collection unit 10 includes a living environment information holding unit 11 that calculates a living environment index from the thermal environment information, a thermal feeling information holding unit 12 that acquires the thermal feeling information declared by the resident of the evaluation target space, and a thermal feeling. It is composed of a data integration unit 13 that extracts the living environment index calculated from the thermal environment information of the measurement date and time corresponding to the generation date and time of the information and generates integrated information that integrates the thermal feeling information and the living environment index. ..

未来モデル生成ユニット2は、評価対象空間が含まれる建物の過去の屋外温熱環境情報を記憶する屋外温熱環境情報管理部20と、複数の学習データ期間によってそれぞれ生成された複数の過去モデルと複数の学習データ期間のそれぞれに対応する過去の屋外温熱環境の代表値とに基づいて、過去モデルの関数の形状を決定付けるモデル特性量と屋外温熱環境の代表値との関係を示す特性量修正関数を求める季節依存管理部21と、居住者の温熱実感を推定したい未来の日時に対して設定されるべき未来の学習データ期間を決定して、この学習データ期間における屋外温熱環境の代表値に対応するモデル特性量を特性量修正関数に基づいて算出して、算出したモデル特性量に基づいて最新の過去モデルを修正した未来モデル(第2のモデル)を生成する未来モデル生成部22(第2のモデル生成部)とを備えている。 The future model generation unit 2 includes an outdoor thermal environment information management unit 20 that stores past outdoor thermal environment information of the building including the evaluation target space, and a plurality of past models and a plurality of past models generated by a plurality of learning data periods. Based on the representative values of the past outdoor thermal environment corresponding to each of the training data periods, the characteristic quantity correction function showing the relationship between the model characteristic quantity that determines the shape of the function of the past model and the representative value of the outdoor thermal environment Determine the future learning data period to be set for the desired season-dependent management unit 21 and the future date and time for which the resident's feeling of heat is to be estimated, and correspond to the representative value of the outdoor thermal environment in this learning data period. Future model generator 22 (second model) that calculates the model characteristic quantity based on the characteristic quantity correction function and generates a future model (second model) in which the latest past model is modified based on the calculated model characteristic quantity. It has a model generator).

図3は居住環境情報保持部11と温熱実感情報保持部12とデータ統合部13の動作を説明するフローチャートである。
居住環境情報保持部11は、例えば評価対象空間101に設置された環境計測デバイス(温度センサ103や湿度センサ104など)から温熱環境情報(温度計測値や湿度計測値などの温熱環境計測値)を受信して、居住環境指標を算出する。そして、居住環境情報保持部11は、温熱環境計測値の計測日時と居住環境指標とを対応付けて居住環境情報(計測日時、居住環境指標)として記憶する(図3ステップS100)。
FIG. 3 is a flowchart illustrating the operation of the living environment information holding unit 11, the thermal feeling information holding unit 12, and the data integration unit 13.
The living environment information holding unit 11 obtains thermal environment information (thermal environment measurement values such as temperature measurement value and humidity measurement value) from, for example, an environment measurement device (temperature sensor 103, humidity sensor 104, etc.) installed in the evaluation target space 101. Receive and calculate the living environment index. Then, the living environment information holding unit 11 stores the measurement date and time of the thermal environment measurement value and the living environment index as living environment information (measurement date and time, living environment index) in association with each other (step S100 in FIG. 3).

居住環境情報保持部11には、居住環境指標の算出に必要な算出式およびパラメータ値、算出に利用する温熱環境計測値の種類の情報などが、室内環境評価/制御等のソリューションプロバイダや設備管理者により予め設定されている。 In the living environment information holding unit 11, the calculation formula and parameter values necessary for calculating the living environment index, information on the type of thermal environment measurement value used for the calculation, and the like are provided to the solution provider such as indoor environment evaluation / control and facility management. It is preset by the person.

居住環境情報保持部11は、データ統合部13からの情報送信要求に応じて、居住環境情報(計測日時、居住環境指標)をデータ統合部13に送信する(図3ステップS101)。 The living environment information holding unit 11 transmits the living environment information (measurement date and time, living environment index) to the data integration unit 13 in response to the information transmission request from the data integration unit 13 (step S101 in FIG. 3).

居住環境指標は、居住環境の温熱的な環境指標であり、空調環境で一般的に計測される空気温度や湿度といった物理環境指標、あるいは物理環境指標を用いて算出される一般的な指標(例えば作用温度、SET*(Standard new Effective Temperature)、PMVなど)である。本実施例では、居住環境指標を国際標準化されているPMVとし、PMVの演算に通常必要となる6要素:空気温度Ta[℃]、湿度H[%]、平均放射温度Tr[℃]、気流速度v[m/s]、代謝量M[met]、着衣量Icl[clo]については各々、以下の値とする。つまり、空気温度Ta[℃]と湿度H[%]は対象建物の空気温度計測値と湿度計測値計測値とし、平均放射温度Trを空気温度Taと等しい値とし、気流速度vを0.1[m/s]、代謝量Mを1.0[met]、着衣量Iclを0.8[clo]としてPMVを求める例で説明する。これら以外でPMV演算に必要な情報はPMVの算出式と共に国際標準などを参照して設定されているとする。 The living environment index is a thermal environmental index of the living environment, and is a physical environment index such as air temperature and humidity generally measured in an air-conditioned environment, or a general index calculated using the physical environment index (for example). Operating temperature, SET * (Standard new Effective Temperature), PMV, etc.). In this embodiment, the living environment index is PMV, which is internationally standardized, and the six elements normally required for PMV calculation: air temperature Ta [° C], humidity H [%], average radiation temperature Tr [° C], and airflow. The speed v [m / s], the amount of metabolism M [met], and the amount of clothing Icl [clo] are set to the following values, respectively. That is, the air temperature Ta [° C.] and the humidity H [%] are the measured values of the air temperature and the measured humidity of the target building, the average radiation temperature Tr is set to be equal to the air temperature Ta, and the airflow velocity v is 0.1. An example of obtaining PMV will be described with [m / s], a metabolic amount M of 1.0 [met], and a clothing amount of Icl of 0.8 [cl]. Other than these, it is assumed that the information required for PMV calculation is set by referring to the international standard and the like together with the PMV calculation formula.

ここで、PMVの算出方法について説明する。PMVの算出方法は、ISO−7730で国際標準化されているので、この標準の算出方法に従えばよい。PMVは、以下の式(1)〜式(3)を用いて求められる。算出プログラムは、ANSI/ASHRAE Standard 55−2010等においても公開されている。
PMV=Q(M)×L ・・・(1)
Q(M)=0.303exp-0.036M+0.028 ・・・(2)
L=M−(C+R+Ed+Es)−(Cre+Ere) ・・・(3)
Here, a method of calculating PMV will be described. Since the PMV calculation method is internationally standardized by ISO-7730, the calculation method of this standard may be followed. PMV is obtained by using the following formulas (1) to (3). The calculation program is also published in ANSI / ASHRAE Standard 55-2010 and the like.
PMV = Q (M) x L ... (1)
Q (M) = 0.303exp -0.036M + 0.028 ... (2)
L = M- (C + R + Ed + Es)-(Cre + Ere) ... (3)

Mは代謝量[met]、Lは人体の熱負荷L[W/m2]、Cは対流熱損失量[W/m2]、Rは放射熱損失量[W/m2]、Edは不感蒸せつ量[W/m2]、Esは発汗による蒸発熱損失量[W/m2]、Creは呼吸による顕熱損失量[W/m2]、Ereは呼吸による潜熱損失量[W/m2]である。 M is the amount of metabolism [met], L is the heat load of the human body L [W / m 2 ], C is the amount of convection heat loss [W / m 2 ], R is the amount of radiant heat loss [W / m 2 ], and Ed is Insensitive steam amount [W / m 2 ], Es is evaporation heat loss amount due to sweating [W / m 2 ], Cre is manifest heat loss amount due to breathing [W / m 2 ], Ere is latent heat loss amount due to breathing [W] / M 2 ].

以下の式(4)〜式(6)を収束計算し、対流熱損失量Cおよび放射熱損失量Rを算出する。tclは着衣外表面温度[℃]である。
C=fcl×hc×(tcl−Ta) ・・・(4)
R=3.96×10-8×fcl×{(tcl+273.15)4
−(Tr+273.15)4} ・・・(5)
tcl=ts−0.155×Icl×(C+R) ・・・(6)
The following equations (4) to (6) are convergently calculated to calculate the convection heat loss amount C and the radiant heat loss amount R. tcl is the clothing outer surface temperature [° C.].
C = fcl × hc × (tcl-Ta) ・ ・ ・ (4)
R = 3.96 × 10 -8 × fcl × {(tcl + 273.15) 4
-(Tr + 273.15) 4 } ... (5)
tcl = ts-0.155 x Icl x (C + R) ... (6)

式(4)〜式(6)において、Iclは着衣量[clo]である。fclは着衣面積増加係数であり、Icl<0.78の時、fcl=1+1.29×Icl、Icl≧0.78の時、fcl=1.05+0.645×Iclである。Trは平均放射温度[℃]である。tsは平均皮膚温度[℃]であり、次式により求められる。Wは機械的仕事量[W/m2]である(通常は0)。
ts=35.7−0.028×(M−W) ・・・(7)
In formulas (4) to (6), Icl is the amount of clothing [clo]. fcl is a clothing area increase coefficient, and when Icl <0.78, fcl = 1 + 1.29 × Icl, and when Icl ≧ 0.78, fcl = 1.05 + 0.645 × Icl. Tr is the average radiation temperature [° C.]. ts is the average skin temperature [° C.] and is calculated by the following equation. W is the mechanical work load [W / m 2 ] (usually 0).
ts = 35.7-0.028 × (MW) ・ ・ ・ (7)

また、hcは人体の対流熱伝達率[W/(m2・℃)]であり、2.38×|tcl−ta|0.25または12.1×v0.5のうちの大きい方となる。vは気流速度[m/s]である。 Further, hc is the convective heat transfer coefficient [W / (m 2 · ° C.)] of the human body, which is the larger of 2.38 × | tcl-ta | 0.25 or 12.1 × v 0.5 . v is the airflow velocity [m / s].

不感蒸せつ量Edは、次式により得られる。
Ed=3.05×(5.73−0.007×M−pa) ・・・(8)
発汗による蒸発熱損失量Esは、次式により得られる。
Es=0.42×(M−58.15) ・・・(9)
The insensitive steaming amount Ed is obtained by the following equation.
Ed = 3.05 × (5.73-0.007 × M-pa) ・ ・ ・ (8)
The amount of heat of vaporization Es due to sweating is obtained by the following equation.
Es = 0.42 × (M-58.15) ・ ・ ・ (9)

呼吸による顕熱損失量Creは、次式により得られる。
Cre=0.0014×M×(34−Ta) ・・・(10)
呼吸による潜熱損失量Ereは、次式により得られる。paは水蒸気圧である。
Ere=0.0173×M×(5.87−pa) ・・・(11)
The amount of sensible heat loss Cre due to respiration is obtained by the following equation.
Cre = 0.0014 × M × (34-Ta) ・ ・ ・ (10)
The latent heat loss amount Ere due to respiration is obtained by the following equation. pa is the vapor pressure.
Ere = 0.0173 × M × (5.87-pa) ・ ・ ・ (11)

以上により、居住環境情報保持部11は、居住環境指標としてPMVを算出することができる。本実施例では、1分周期で計測された温熱環境計測値が居住環境情報保持部11に送信され、居住環境情報保持部11が1分周期でPMVを算出するものとする。空調の制御装置や建物を管理する中央監視システムなどに周期的に温熱環境計測値が蓄積されている場合は、これら空調の制御装置や中央監視システムから温熱環境計測値を取得すればよい。 From the above, the living environment information holding unit 11 can calculate PMV as a living environment index. In this embodiment, it is assumed that the thermal environment measurement value measured in a 1-minute cycle is transmitted to the living environment information holding unit 11, and the living environment information holding unit 11 calculates the PMV in a 1-minute cycle. If the thermal environment measurement values are periodically accumulated in the air conditioning control device or the central monitoring system that manages the building, the thermal environment measurement values may be acquired from these air conditioning control devices or the central monitoring system.

一方、温熱実感情報保持部12は、評価対象空間101の居住者の温熱実感情報を受信し、この温熱実感情報の発生日時や温熱実感情報を申告した居住者の情報と関連付けて温熱実感情報を記憶する(図3ステップS102)。 On the other hand, the heat feeling information holding unit 12 receives the heat feeling information of the resident of the evaluation target space 101, and associates the heat feeling information with the information of the resident who declared the heat feeling information and the generation date and time of the heat feeling information. Store it (step S102 in FIG. 3).

温熱実感情報保持部12は、データ統合部13からの情報送信要求に応じて、温熱実感情報と発生日時情報と居住者情報とをデータ統合部13に送信する(図3ステップS103)。 The thermal feeling information holding unit 12 transmits the thermal feeling information, the occurrence date / time information, and the resident information to the data integrating unit 13 in response to the information transmission request from the data integration unit 13 (step S103 in FIG. 3).

温熱実感情報は、居住環境指標に対応する居住者の感じ方を示す量であり、暑い/寒いなどの全身温冷感の尺度や、温熱環境に起因する快適感/満足感などの情報を示す。対象空間の不満足者数や、不満足者数から算出される不満足者率などの指標も温熱実感情報に含まれる。温熱実感情報の入力端末としては、スマートフォン、PC(personal computer)などを適宜利用すればよい。 The thermal sensation information is an amount indicating the resident's feeling corresponding to the living environment index, and indicates information such as a scale of whole body warming sensation such as hot / cold and comfort / satisfaction caused by the thermal environment. .. Indicators such as the number of dissatisfied persons in the target space and the dissatisfied person rate calculated from the number of dissatisfied persons are also included in the thermal feeling information. As an input terminal for the heat feeling information, a smartphone, a PC (personal computer), or the like may be appropriately used.

本実施例では、居住者からの温冷感申告情報(暑い/寒い)を温熱実感情報とし、居住者が自由に暑い、あるいは寒いという温冷感申告を入力端末を通じて行うものとする。温冷感申告情報と申告日時(発生日時情報)と申告者ID(居住者情報)の1例を図4に示す。図4の例では、暑いを「hot」、寒いを「cold」で示している。 In this embodiment, the hot / cold feeling report information (hot / cold) from the resident is used as the hot / cold feeling information, and the resident freely makes a hot / cold feeling report through the input terminal. FIG. 4 shows an example of the warm / cold feeling declaration information, the declaration date / time (occurrence date / time information), and the reporter ID (resident information). In the example of FIG. 4, hot is indicated by "hot" and cold is indicated by "cold".

申告者IDは、居住者が温冷感申告と合わせて入力して温熱実感情報保持部12に送信してもよいし、居住者個々人が入力端末を利用する場合は端末IDを申告者IDとして、温熱実感情報保持部12に送信してもよい。
また、申告日時は、入力端末が温熱実感情報保持部12に送信してもよいし、温熱実感情報保持部12が温冷感申告情報を受信した日時を申告日時としてもよい。
The reporter ID may be input by the resident together with the heat / cold feeling report and transmitted to the heat / cold feeling information holding unit 12, or when the individual resident uses the input terminal, the terminal ID is used as the filer ID. , May be transmitted to the thermal feeling information holding unit 12.
Further, the date and time of the declaration may be transmitted by the input terminal to the heat sensation information holding unit 12, or may be the date and time when the heat sensation information holding unit 12 receives the heat sensation report information.

次に、データ統合部13は、居住環境情報保持部11から居住環境指標を受信し、温熱実感情報保持部12から温熱実感情報を受信すると、温熱実感情報の発生日時に対応する計測日時の居住環境指標を抽出して、温熱実感情報と居住環境指標とを統合した統合情報(発生日時情報、居住者ID、温熱実感情報、居住環境指標)を記憶する(図3ステップS104)。 Next, when the data integration unit 13 receives the living environment index from the living environment information holding unit 11 and receives the heat feeling information from the heat feeling information holding unit 12, the data integration unit 13 resides on the measurement date and time corresponding to the generation date and time of the heat feeling information. The environmental index is extracted, and integrated information (occurrence date / time information, resident ID, thermal feeling information, living environment index) that integrates the thermal feeling information and the living environment index is stored (FIG. 3, step S104).

データ統合部13は、過去モデル生成部15からの情報送信要求に応じて、統合情報を過去モデル生成部15に送信する(図3ステップS105)。 The data integration unit 13 transmits the integrated information to the past model generation unit 15 in response to the information transmission request from the past model generation unit 15 (step S105 in FIG. 3).

温熱実感情報発生日時情報(申告日時情報)と居住者ID(申告者ID)と温熱実感情報と居住環境指標(PMV)とからなる統合情報の1例を図5に示す。
温熱実感情報の発生日時に対応する計測日時の居住環境指標を抽出する際に、温熱実感情報の発生日時(申告日時)の分解能が例えば1分単位で、温熱環境計測値の計測日時の分解能が例えば10分単位などというように、温熱実感情報の発生と温熱環境計測値の計測の時間間隔が異なる場合がある。この場合には、予めソリューションプロバイダや設備管理者が抽出ルールを決定しておけばよい。
FIG. 5 shows an example of integrated information consisting of heat feeling information generation date / time information (declaration date / time information), resident ID (reporter ID), heat feeling information, and living environment index (PMV).
When extracting the living environment index of the measurement date and time corresponding to the occurrence date and time of the thermal feeling information, the resolution of the generation date and time (declaration date and time) of the thermal actual feeling information is, for example, 1 minute unit, and the resolution of the measurement date and time of the thermal environment measurement value For example, the time interval between the generation of the thermal feeling information and the measurement of the thermal environment measurement value may be different, such as in units of 10 minutes. In this case, the solution provider or the equipment manager may determine the extraction rule in advance.

例えば温熱実感情報の発生日時(申告日時)の分解能が1分単位で、温熱環境計測値が10分おきに計測される場合、温熱実感情報の発生日時の直近の計測日時、または温熱実感情報の発生日時を10分単位で繰り上げた日時を、温熱実感情報の発生日時に対応する計測日時とする、という抽出ルールを決定しておけばよい。 For example, if the resolution of the date and time of occurrence of the heat sensation information (declaration date and time) is in units of 1 minute and the measured value of the thermal environment is measured every 10 minutes, the latest measurement date and time of the heat sensation information or the heat sensation information It is sufficient to determine the extraction rule that the date and time when the occurrence date and time is advanced in units of 10 minutes is set as the measurement date and time corresponding to the occurrence date and time of the thermal feeling information.

14時50分、15時00分、15時10分、・・・・といったように温熱環境計測値が10分おきに計測される場合、温熱実感情報の発生日時15時03分に対応する計測日時は、直近の15時00分、または15時03分を10分単位で繰り上げた15時10分となる。 When the thermal environment measurement value is measured every 10 minutes, such as 14:50, 15:00, 15:10, etc., the measurement corresponds to the occurrence date and time of the thermal feeling information at 15:03. The date and time will be 15:10, which is the latest 15:00 or 15:03, which is advanced by 10 minutes.

居住環境情報保持部11と温熱実感情報保持部12とデータ統合部13とは、継続的に図3の処理を実行し、データ統合部13に統合情報が蓄積される。 The living environment information holding unit 11, the thermal feeling information holding unit 12, and the data integration unit 13 continuously execute the process of FIG. 3, and the integrated information is accumulated in the data integration unit 13.

次に、データ期間管理部14には、個別温熱実感モデル(過去モデル)の更新タイミングと学習データ期間とを決定するためのデータ期間管理情報が室内環境評価/制御等のソリューションプロバイダや設備管理者によって予め設定されている。
データ期間管理部14は、過去モデル生成部15または未来モデル生成部22からの情報送信要求に応じ、データ期間管理情報を送信する。
Next, in the data period management unit 14, the data period management information for determining the update timing and the learning data period of the individual thermal feeling model (past model) is provided to the solution provider or equipment manager such as indoor environment evaluation / control. It is preset by.
The data period management unit 14 transmits data period management information in response to an information transmission request from the past model generation unit 15 or the future model generation unit 22.

本実施例では、個別温熱実感モデルの更新タイミングの情報と学習データ期間幅の情報とをデータ期間管理情報とする。具体的には、個別温熱実感モデルの更新タイミングを、毎月1日(X月1日、Xは任意の月)の0時00分とし、学習データ期間幅を、前回の更新タイミングから今回の更新タイミングまでの1か月とする。つまり、個別温熱実感モデルの更新タイミングが毎月1日であるので、この個別温熱実感モデルの更新のための学習データ期間の開始日時は前月の1日の0時00分、学習データ期間の終了日時は前月の末日の23時59分となる。 In this embodiment, the update timing information of the individual thermal feeling model and the learning data period width information are used as data period management information. Specifically, the update timing of the individual thermal feeling model is set to 0:00 on the first day of every month (X month 1, X is an arbitrary month), and the learning data period width is updated from the previous update timing to this time. One month until the timing. That is, since the update timing of the individual thermal feeling model is the first day of every month, the start date and time of the learning data period for updating this individual thermal feeling model is 0:00 on the first day of the previous month, and the end date and time of the learning data period. Will be 23:59 on the last day of the previous month.

データ期間管理情報から個別温熱実感モデルの更新タイミングと学習データ期間とを決定できればよいので、更新スタート日時と、その更新スタート日時から遡る学習データの期間幅(14日間毎や168時間毎など)とをデータ期間管理情報として設定してももちろん構わない。 Since it is only necessary to determine the update timing and training data period of the individual thermal feeling model from the data period management information, the update start date and time and the period width of the learning data (every 14 days, every 168 hours, etc.) that goes back from the update start date and time Of course, may be set as data period management information.

図6は過去モデル生成部15の動作を説明するフローチャートである。過去モデル生成部15は、データ期間管理部14から取得したデータ期間管理情報により、個別温熱実感モデルの更新タイミングになったことを認識すると(図6ステップS200においてYES)、今回のモデル更新に対応する学習データ期間の開始日時と終了日時とを決定する(図6ステップS201)。 FIG. 6 is a flowchart illustrating the operation of the past model generation unit 15. When the past model generation unit 15 recognizes that the update timing of the individual thermal feeling model has come based on the data period management information acquired from the data period management unit 14 (YES in step S200 of FIG. 6), it corresponds to the current model update. The start date and time and the end date and time of the learning data period to be performed are determined (FIG. 6, step S201).

過去モデル生成部15は、決定した学習データ期間の統合情報をデータ統合部13から取得し(図6ステップS202)、取得した統合情報を利用して、個別温熱実感モデルを同定するための入出力データ(モデリング情報)を生成する(図6ステップS203)。過去モデル生成部15は、生成したモデリング情報を利用して、個別温熱実感モデルを同定し、この個別温熱実感モデルを最新の過去モデルとして更新する(図6ステップS204)。 The past model generation unit 15 acquires the integrated information of the determined learning data period from the data integration unit 13 (step S202 in FIG. 6), and uses the acquired integrated information to input / output for identifying the individual thermal feeling model. Data (modeling information) is generated (FIG. 6, step S203). The past model generation unit 15 uses the generated modeling information to identify the individual thermal sensation model, and updates the individual thermal sensation model as the latest past model (FIG. 6, step S204).

そして、過去モデル生成部15は、最新の過去モデルのモデリング結果とモデリングに利用した学習データ期間の情報とを過去モデル記憶管理部16に送信する(図6ステップS205)。
過去モデル生成部15は、図6の処理をモデルの更新タイミング毎に行う。
Then, the past model generation unit 15 transmits the modeling result of the latest past model and the information of the learning data period used for modeling to the past model storage management unit 16 (step S205 in FIG. 6).
The past model generation unit 15 performs the process of FIG. 6 at each model update timing.

本実施例では、上記のとおり個別温熱実感モデルの更新タイミングを、毎月1日の0時00分としている。例えば個別温熱実感モデルが2018年3月1日0時00分に更新される場合、学習データ期間は2018年2月1日0時00分から2月28日23時59分までとなる。また、個別温熱実感モデルが2019年8月1日0時00分に更新される場合、学習データ期間は2018年7月1日0時00分から7月31日23時59分までとなる。 In this embodiment, as described above, the update timing of the individual thermal feeling model is set to 0:00 on the first day of every month. For example, when the individual thermal feeling model is updated at 0:00 on March 1, 2018, the learning data period is from 0:00 on February 1, 2018 to 23:59 on February 28, 2018. If the individual thermal experience model is updated at 0:00 on August 1, 2019, the learning data period will be from 0:00 on July 1, 2018 to 23:59 on July 31, 2018.

また、本実施例では、居住環境指標(PMV)と、「暑い」という申告者率PPVhとをモデリング情報とし、予測不満足者率PPD(Predicted Percentage Dissatisfied)の関数を個別温熱実感モデルの基本関数とする。PPDも国際標準化されており、ISO−7730でPMVを用いた演算式が示されている。 Further, in this embodiment, the living environment index (PMV) and the reporter rate PPVh of "hot" are used as modeling information, and the function of the predicted dissatisfaction rate PPD (Predicted Percentage Dissatisfied) is used as the basic function of the individual thermal feeling model. To do. PPD has also been internationally standardized, and ISO-7730 shows an arithmetic expression using PMV.

次に、ステップS203のモデリング情報の生成方法について説明する。上記のとおり、過去モデル生成部15は、統合情報からモデリング情報を生成する。
任意のPMVにおける「暑い」という申告者率PPVhは、任意のPMVにおける「暑い」という申告者数をN、評価対象空間における居住者数をN_allとしたとき、次式により算出される。
PPVh=N/N_all ・・・(12)
Next, a method of generating modeling information in step S203 will be described. As described above, the past model generation unit 15 generates modeling information from the integrated information.
The rate of filers who say "hot" in any PMV PPVh is calculated by the following equation, where N is the number of reporters who say "hot" in any PMV and the number of residents in the evaluation target space is N_all.
PPVh = N / N_all ... (12)

評価対象空間における居住者数N_allは、室内環境評価/制御等のソリューションプロバイダや設備管理者により予め設定されている。
ここで、式(12)における申告者数Nは、温冷感に対する生理的な妥当性を考慮して算出した申告者数としてもよい。この算出方法は、任意のPMVにおいて「暑い」という申告をした人は、その環境よりも更に暑い側のPMVではより暑く感じるという原理に基づくものである。
The number of residents N_all in the evaluation target space is preset by the solution provider for indoor environment evaluation / control and the facility manager.
Here, the number of filers N in the formula (12) may be the number of filers calculated in consideration of the physiological validity for the feeling of warmth and cold. This calculation method is based on the principle that a person who declares "hot" in any PMV feels hotter in the PMV on the hotter side than the environment.

具体的には、各々の申告者が、PMVが低い(暑さが弱い)ときに最初に申告を行ったPMVをラッチ開始点PMVstart(m)(mは正の整数、本実施例ではm=1〜4)とし、ラッチ開始点PMVstart(m)よりも大きいPMV領域では「暑い」という申告の発生状態を維持して人数を算出する。この算出方法によれば、図7に示すように「暑い」という申告者の数は、より暑い環境になるに従って増加し、生理的な妥当性が高くなる。以降ではこの方法を申告ラッチ法と記述し、この申告ラッチ法によりモデリング情報を生成する例を説明する。 Specifically, each filer latches the PMV that was first declared when the PMV is low (the heat is weak) at the latch start point PMV start (m) (m is a positive integer, m = in this embodiment. In the PMV region larger than the latch start point PMV start (m), the number of people is calculated while maintaining the state of occurrence of the declaration of "hot". According to this calculation method, as shown in FIG. 7, the number of reporters who say "hot" increases as the environment becomes hotter, and the physiological validity becomes higher. Hereinafter, this method will be described as the declaration latch method, and an example of generating modeling information by this declaration latch method will be described.

図8に申告ラッチ法を用いて算出した申告者数から申告者率PPVhを算出して、モデリング情報を生成するフローチャートを示す。過去モデル生成部15は、ステップS202で取得した統合情報から「暑い」という申告の情報のみを抽出し(図8ステップS300)、抽出した情報に含まれる申告者IDが異なる申告者の数(「暑い」という申告をした居住者の人数)N_voteを求める(図8ステップS301)。図5に示したデータ例を、学習データ期間の統合情報とすると、「暑い」(hot)という申告をした申告者の数はN_vote=4である。 FIG. 8 shows a flowchart for generating modeling information by calculating the filer rate PPVh from the number of filers calculated by using the report latch method. The past model generation unit 15 extracts only the information of the declaration of "hot" from the integrated information acquired in step S202 (step S300 of FIG. 8), and the number of filers with different filer IDs included in the extracted information ("" The number of residents who declared "hot") N_vote is requested (Fig. 8, step S301). Assuming that the data example shown in FIG. 5 is the integrated information of the learning data period, the number of filers who have declared "hot" is N_vote = 4.

次に、過去モデル生成部15は、申告者を数えるための変数m(mはN_vote以下の正の整数)1に初期化する(図8ステップS302)。
過去モデル生成部15は、m番目の申告者について、「暑い」という申告が行われたときのPMVのうち最小値をラッチ開始点PMVstart(m)として求める(図8ステップS303)。
Next, the past model generation unit 15 initializes the variable m (m is a positive integer less than or equal to N_vote) 1 for counting the filers (step S302 in FIG. 8).
The past model generation unit 15 obtains the minimum value of the PMV when the declaration of "hot" is made for the m-th filer as the latch start point PMV start (m) (step S303 in FIG. 8).

過去モデル生成部15は、m=N_vote、すなわち「暑い」という申告をした全申告者についてステップS303の処理を終えたかどうかを判定し(図8ステップS304)、変数mがN_vote未満の場合には処理を終えていないとして、変数mを1増やす(図8ステップS305)。こうして、「暑い」という申告をした申告者毎にステップS303の処理が行われる。 The past model generation unit 15 determines whether m = N_vote, that is, whether or not the processing of step S303 has been completed for all the filers who have declared “hot” (step S304 in FIG. 8), and if the variable m is less than N_vote, Assuming that the processing is not completed, the variable m is incremented by 1 (step S305 in FIG. 8). In this way, the process of step S303 is performed for each filer who declares "hot".

過去モデル生成部15は、m=N_voteとなり、「暑い」という申告をした全申告者についてステップS303の処理を終えたときに(ステップS304においてYES)、任意のPMVとこれに対応する「暑い」という申告者率PPVhからなるモデリング情報を生成する(図8ステップS306)。 The past model generation unit 15 has m = N_vote, and when the processing of step S303 is completed for all the filers who have declared “hot” (YES in step S304), any PMV and the corresponding “hot” Modeling information consisting of the filer rate PPVh is generated (FIG. 8 step S306).

任意のPMVに対する申告者P(m)の申告の有無の状態V_latch(m)は、PMV<PMVstart(m)のとき、V_latch(m)=0(申告無し)、PMVstart(m)≦PMVのとき、V_latch(m)=1(申告有り)となる。
このとき、任意のPMVに対する「暑い」という申告者率PPVhは次式により得られる。
The state of presence / absence of declaration of the filer P (m) for any PMV When V_latch (m) is PMV <PMVstart (m), V_latch (m) = 0 (no declaration), PMVstart (m) ≤ PMV , V_latch (m) = 1 (with declaration).
At this time, the reporter rate PPVh of "hot" for any PMV is obtained by the following equation.

以上により、「暑い」という申告者率PPVhの算出が終了し、学習データ期間のPMVと、学習データ期間の「暑い」という申告者率PPVhとからなるモデリング情報の生成が終了する。図5のデータ例の場合、図7の縦軸を「暑い」という申告者の人数で除算し、「暑い」という申告者率としたものがモデリング情報となる。モデリング情報は、PMVを一定の刻みで生成しても良いし、PPVhが増加する点のみのデータとするなどとしても構わない。 As a result, the calculation of the reporter rate PPVh of "hot" is completed, and the generation of modeling information including the PMV of the learning data period and the reporter rate PPVh of "hot" in the learning data period is completed. In the case of the data example of FIG. 5, the vertical axis of FIG. 7 is divided by the number of filers who say "hot" to obtain the rate of filers who say "hot", which is the modeling information. The modeling information may be generated at regular intervals of PMV, or may be data only at points where PPVh increases.

次に、ステップS204の個別温熱実感モデルの同定方法について説明する。上記のとおり、本実施例では、予測不満足者率PPDを個別温熱実感モデルの基本関数とする。PPDは、PMVの関数として以下の関係式で示される。
PPD=F(PMV)=100−95×exp(−0.03353×PMV4
−0.2179×PMV2)[%] ・・・(14)
Next, the method of identifying the individual thermal feeling model in step S204 will be described. As described above, in this embodiment, the predicted dissatisfaction rate PPD is used as the basic function of the individual thermal experience model. PPD is represented by the following relational expression as a function of PMV.
PPD = F (PMV) = 100-95 x exp (-0.03353 x PMV 4)
−0.2179 × PMV 2 ) [%] ・ ・ ・ (14)

PMVとPPDとの関係は図9のようになる。個別温熱実感モデルは、PMVとPPDとを対応付けた、図9のようなV字型の分布関数として定義すればよい。あるいは、暑い側と寒い側の不満足度を分離して、不満足関数(個別温熱実感モデル)を定義することも可能である。 The relationship between PMV and PPD is as shown in FIG. The individual thermal feeling model may be defined as a V-shaped distribution function as shown in FIG. 9 in which PMV and PPD are associated with each other. Alternatively, it is also possible to separate the degree of dissatisfaction on the hot side and the cold side and define the dissatisfaction function (individual thermal sensation model).

図10は特開2015−4480号公報に開示された不満足度関数の例を示す図である。「暑い」という不満足度を求める場合には、図10の110で示す不満足度関数を使用すればよく、「寒い」という不満足度を求める場合には、図10の111で示す不満足度関数を使用すればよい。図10のような関数としては、例えばシグモイド関数などがあり、個別温熱実感モデルの関数はこのような一般的なモデルで構わない。 FIG. 10 is a diagram showing an example of a dissatisfaction function disclosed in Japanese Patent Application Laid-Open No. 2015-4480. The dissatisfaction function shown in 110 of FIG. 10 may be used to obtain the dissatisfaction of "hot", and the dissatisfaction function shown in 111 of FIG. 10 may be used to obtain the dissatisfaction of "cold". do it. Examples of the function as shown in FIG. 10 include a sigmoid function, and the function of the individual thermal feeling model may be such a general model.

本実施例では、式(14)の基本関数に居住者の温熱実感を反映させて関数を修正した、以下の式(15)のモデルを個別温熱実感モデル(不満足度モデル)とした例で説明する。この式(15)の個別温熱実感モデルの形状は、図10の110で示した関数と同様のものである。式(15)では、不満足度に相当する式(14)の予測不満足者率PPDを居住者の温熱実感を反映させた「暑い」という申告者率PPVhとしている。 In this embodiment, the basic function of the equation (14) is modified by reflecting the resident's feeling of heat, and the model of the following formula (15) is used as an individual feeling of heat model (dissatisfaction model). To do. The shape of the individual thermal feeling model of the equation (15) is the same as the function shown by 110 in FIG. In the formula (15), the predicted dissatisfaction rate PPD of the formula (14) corresponding to the degree of dissatisfaction is set to the reporter rate PPVh of "hot" reflecting the feeling of heat of the resident.

式(15)において、a,b,c,d,eは探索対象となるモデルのパラメータである。式(16)は、式(15)においてa=b(「暑い」という申告者率PPVhの最小値が0)としたものであり、これを個別温熱実感モデルのモデル関数とする場合、探索対象となるモデルのパラメータはa,c,d,eとなる。 In equation (15), a, b, c, d, and e are parameters of the model to be searched. Equation (16) is the one in which a = b (the minimum value of the reporter rate PPVh of "hot" is 0) in equation (15), and when this is used as the model function of the individual thermal feeling model, the search target. The parameters of the model are a, c, d, and e.

式(15)または式(16)のような個別温熱実感モデルのモデル関数と、探索するモデルのパラメータと、利用する最適化手法とが、室内環境評価/制御等のソリューションプロバイダや設備管理者により、予め設定されている。 The model function of the individual thermal feeling model such as equation (15) or equation (16), the parameters of the model to be searched, and the optimization method to be used are determined by the solution provider or equipment manager such as indoor environment evaluation / control. , Preset.

過去モデル生成部15は、ステップS203で生成したモデリング情報を利用して、ソリューションプロバイダや設備管理者により予め設定された個別温熱実感モデルのモデル関数に対し、同様に予め設定された探索対象のモデルパラメータを探索して個別温熱実感モデルを同定する。パラメータを探索する際に利用する最適化手法としては、最小二乗法やシンプレックス法など汎用的な手法を利用すればよい。
以上により、過去モデルの更新が終了する。
The past model generation unit 15 uses the modeling information generated in step S203 to search the model function of the individual thermal feeling model preset by the solution provider or the equipment manager. Search the parameters to identify the individual thermal sensation model. As the optimization method used when searching for parameters, a general-purpose method such as the least squares method or the simplex method may be used.
With the above, the update of the past model is completed.

過去モデル記憶管理部16は、過去モデル生成部15から受信した過去モデルのモデリング結果(探索したモデルのパラメータ値)と学習データ期間の情報とを記憶し、季節依存管理部21からのモデル情報の送信要求に応じて、記憶している情報を季節依存管理部21に送信する。 The past model storage management unit 16 stores the modeling result (parameter value of the searched model) of the past model received from the past model generation unit 15 and the training data period information, and stores the model information from the season-dependent management unit 21. In response to the transmission request, the stored information is transmitted to the seasonal dependence management unit 21.

屋外温熱環境情報管理部20は、評価対象空間101が含まれる建物の過去の屋外温熱環境情報と未来の屋外温熱環境情報を記憶している。過去の屋外温熱環境情報は、屋外温熱環境状態情報と、屋外温熱環境状態情報を取得した日時の情報とからなる。未来の屋外温熱環境情報は、未来の屋外温熱環境を推定するために必要な情報であり、屋外温熱環境状態情報と、屋外温熱環境状態情報の出現が推定される日時の情報とからなる。 The outdoor thermal environment information management unit 20 stores the past outdoor thermal environment information and the future outdoor thermal environment information of the building including the evaluation target space 101. The past outdoor thermal environment information consists of outdoor thermal environment status information and information on the date and time when the outdoor thermal environment status information was acquired. The future outdoor thermal environment information is information necessary for estimating the future outdoor thermal environment, and includes information on the outdoor thermal environment state information and information on the date and time when the appearance of the outdoor thermal environment state information is estimated.

季節推移を反映する屋外温熱環境状態情報としては、屋外の気温や湿度、気温と湿度とを統合した不快指数などの指標、WBGT(湿球黒球温度)などがある。また、過去の屋外温熱環境情報は、平年のデータ(月平均気温、日平均気温等)などの統計量としてもよいし、屋外に設置された温度センサ等による計測値としてもよい。未来の屋外温熱環境情報は、本質的には未来の学習データ期間に対応する過去の屋外温熱環境情報であり、過去の屋外温熱環境情報が未来で再現するとして過去の屋外温熱環境情報と同様としてもよいし、気象予報などに基づいて推定される屋外温熱環境情報の推定値としてもよい。 As outdoor thermal environment state information that reflects seasonal changes, there are indicators such as outdoor air temperature and humidity, an index such as a discomfort index that integrates temperature and humidity, and WBGT (wet-bulb globe temperature). Further, the past outdoor thermal environment information may be a statistic such as normal year data (monthly average temperature, daily average temperature, etc.), or may be a measured value by a temperature sensor or the like installed outdoors. The future outdoor thermal environment information is essentially the past outdoor thermal environment information corresponding to the future learning data period, and is similar to the past outdoor thermal environment information as the past outdoor thermal environment information is reproduced in the future. It may be used as an estimated value of outdoor thermal environment information estimated based on a weather forecast or the like.

屋外温熱環境情報の取得方法については、室内環境評価/制御等のソリューションプロバイダや設備管理者によって予め設定されている。屋外温熱環境状態情報として月平均気温などの、予め入手可能なデータを用いる場合には、これらのデータを予め屋外温熱環境情報管理部20に設定しておくことで構わない。
本実施例では、平年の各月の月平均外気温を過去および未来の屋外温熱環境状態情報として使用する。
The method of acquiring outdoor thermal environment information is preset by the solution provider such as indoor environment evaluation / control and the equipment manager. When data available in advance such as monthly average temperature is used as the outdoor thermal environment state information, these data may be set in the outdoor thermal environment information management unit 20 in advance.
In this embodiment, the monthly average outside temperature of each month of the normal year is used as past and future outdoor thermal environmental condition information.

次に、季節依存管理部21は、未来モデル生成部22からの要求により、特性量修正関数Gpmを決定し、この特性量修正関数Gpmと過去モデル記憶管理部16に記憶されている過去モデルのモデリング結果(探索したモデルのパラメータ値)とを未来モデル生成部22へ送信する。 Next, the seasonal dependence management unit 21 determines the characteristic quantity correction function Gpm in response to the request from the future model generation unit 22, and the characteristic quantity correction function Gpm and the past model stored in the past model storage management unit 16 The modeling result (parameter value of the searched model) is transmitted to the future model generation unit 22.

図11は季節依存管理部21の動作を説明するフローチャートである。季節依存管理部21は、最新のモデリング結果から(N_ref−1)回前のモデリング結果までのN_ref回(N_refは予め設定された参照モデル数)のモデリング結果と、これらN_ref回のモデリングに利用した学習データ期間の情報とを過去モデル記憶管理部16から取得する(図11ステップS400)。 FIG. 11 is a flowchart illustrating the operation of the seasonal dependence management unit 21. The season-dependent management unit 21 used the modeling results of N_ref times (N_ref is a preset number of reference models) from the latest modeling result to the modeling result of (N_ref-1) times before, and modeling of these N_ref times. Information on the training data period is acquired from the past model storage management unit 16 (step S400 in FIG. 11).

続いて、季節依存管理部21は、N_ref回のモデリングに利用した学習データ期間の各々に対応する過去の屋外温熱環境情報を屋外温熱環境情報管理部20から取得する(図11ステップS401)。 Subsequently, the seasonally dependent management unit 21 acquires the past outdoor thermal environment information corresponding to each of the learning data periods used for the modeling of N_ref times from the outdoor thermal environment information management unit 20 (FIG. 11, step S401).

そして、季節依存管理部21は、取得した過去の屋外温熱環境情報に基づいて、N_ref回のモデリングに利用した学習データ期間に対応する過去の屋外温熱環境代表値Eotを学習データ期間毎に算出する(図11ステップS402)。例えば、過去の屋外温熱環境情報管理部20から取得した屋外温熱環境状態情報が1時間毎に計測された外気温であり、学習データ期間が2018年8月1日0時00分から8月10日23時59分であれば、この学習データ期間の外気温の平均値(平均外気温)を算出して過去の屋外温熱環境代表値Eotとすればよい。 Then, the seasonal dependence management unit 21 calculates the past outdoor thermal environment representative value Eot corresponding to the learning data period used for modeling N_ref times for each learning data period based on the acquired past outdoor thermal environment information. (FIG. 11 step S402). For example, the outdoor thermal environment state information acquired from the past outdoor thermal environment information management unit 20 is the outside air temperature measured every hour, and the learning data period is from 0:00 on August 1, 2018 to August 10, 2018. At 23:59, the average value (average outside air temperature) of the outside air temperature during this learning data period may be calculated and used as the past outdoor thermal environment representative value Eot.

ただし、上記のように、平年の各月の月平均外気温を過去の屋外温熱環境状態情報として使用する場合、季節依存管理部21は、学習データ期間を含む月の月平均外気温をそのまま過去の屋外温熱環境代表値Eotとすればよい。また、複数の月にわたって学習データ期間が設定されている場合には、これら複数の月の月平均外気温の平均値や日数に応じた重み付け平均値を過去の屋外温熱環境代表値Eotとすればよい。 However, as described above, when the monthly average outside temperature of each month of the normal year is used as the past outdoor thermal environment state information, the seasonal dependence management unit 21 keeps the past monthly average outside temperature of the month including the learning data period as it is. The representative value of the outdoor thermal environment of Eot may be used. In addition, when the training data period is set over a plurality of months, the average value of the monthly average outside temperature of these multiple months and the weighted average value according to the number of days can be set as the past outdoor thermal environment representative value Eot. Good.

次に、季節依存管理部21は、過去の屋外温熱環境代表値Eotと過去モデルのモデル特性量SZとの関係を示す特性量修正関数Gpmを求める(図11ステップS403)。モデル特性量SZは、過去モデルの関数の形状を決定付ける量のひとつである。特性量修正関数Gpmが決定されれば、任意の屋外温熱環境代表値Eotに対応するモデル特性量SZが推定可能となる。 Next, the seasonal dependence management unit 21 obtains a characteristic quantity correction function Gpm showing the relationship between the past outdoor thermal environment representative value Eot and the model characteristic quantity SZ of the past model (FIG. 11, step S403). The model characteristic quantity SZ is one of the quantities that determines the shape of the function of the past model. Once the characteristic quantity correction function Gpm is determined, the model characteristic quantity SZ corresponding to an arbitrary outdoor thermal environment representative value Eot can be estimated.

図12は季節依存管理部21の特性量修正関数決定動作を説明するフローチャートである。以下の説明において、iは過去モデルを数えるための変数であり、1〜N_refの整数とする。変数iが小さいほど、新しい過去モデルであることを示している。 FIG. 12 is a flowchart illustrating the operation of determining the characteristic quantity correction function of the seasonal dependence management unit 21. In the following description, i is a variable for counting past models, and is an integer of 1 to N_ref. The smaller the variable i, the newer the past model.

また、図13に示すように、モデル特性量SZを、「暑い」という申告者率PPVhの飽和値とし、特性量修正関数Gpmを、過去の屋外温熱環境代表値Eotとモデル特性量SZとの線形近似の関数とする。つまり、任意の屋外温熱環境代表値Eotが決定されれば、これに対応するモデル特性量SZの推定値が特性量修正関数Gpmによって算出できる。個別温熱実感モデルを式(15)、式(16)に設定した場合、申告者率PPVhの飽和値はモデルの探索パラメータa,b,c,d,eのうちのaに対応する。 Further, as shown in FIG. 13, the model characteristic quantity SZ is set to the saturation value of the reporter rate PPVh of "hot", and the characteristic quantity correction function Gpm is the past outdoor thermal environment representative value Eot and the model characteristic quantity SZ. Let it be a function of linear approximation. That is, if an arbitrary outdoor thermal environment representative value Eot is determined, the estimated value of the model characteristic amount SZ corresponding to this can be calculated by the characteristic amount correction function Gpm. When the individual thermal feeling model is set in the equations (15) and (16), the saturation value of the filer rate PPVh corresponds to a of the model search parameters a, b, c, d, and e.

季節依存管理部21は、ステップS400でN_ref回のモデリング結果Sp(i)、すなわち、i番目の過去モデルに対して探索された探索パラメータ値のセット(本実施例では、探索されたパラメータa(i),b(i),c(i),d(i),e(i)のセット)と、これらN_ref回のモデリングに利用した学習データ期間LP(i)の情報とを過去モデル記憶管理部16から取得した後、変数iを1に初期化する(図12ステップS500)。 The seasonally dependent management unit 21 sets the modeling result Sp (i) of N_ref times in step S400, that is, the search parameter value searched for the i-th past model (in this embodiment, the searched parameter a (in this embodiment). i), b (i), c (i), d (i), e (i) set) and the information of the training data period LP (i) used for modeling these N_ref times are stored in the past model. After acquiring from the part 16, the variable i is initialized to 1 (step S500 in FIG. 12).

季節依存管理部21は、最新からi回目のモデリングに利用した学習データ期間LP(i)に対応する過去の屋外温熱環境代表値Eot(i)をステップS402で算出した値の中から取得する(図12ステップS501)。 The seasonal dependence management unit 21 acquires the past outdoor thermal environment representative value Eot (i) corresponding to the learning data period LP (i) used for the i-th modeling from the latest from the values calculated in step S402 ( FIG. 12 step S501).

季節依存管理部21は、i=N_ref、すなわちN_ref個の全ての学習データ期間LP(i)についてステップS501の処理を終えたかどうかを判定し(図12ステップS502)、変数iがN_ref未満の場合には処理を終えていないとして、変数iを1増やす(図12ステップS503)。こうして、学習データ期間LP(i)毎にステップS501の処理が行われる。 The season-dependent management unit 21 determines whether i = N_ref, that is, whether or not the processing of step S501 has been completed for all N_ref training data period LPs (i) (FIG. 12 step S502), and when the variable i is less than N_ref. The variable i is incremented by 1 (FIG. 12, step S503), assuming that the processing has not been completed. In this way, the process of step S501 is performed for each learning data period LP (i).

季節依存管理部21は、i=N_refとなり、N_ref個の全ての学習データ期間LP(i)についてステップS501の処理を終えたときに(ステップS502においてYES)、個別温熱実感モデルの関数とパラメータ探索値のセットであるモデリング結果Sp(i)より決定するモデル特性量SZ(i)と屋外温熱環境代表値Eot(i)との関係を線形近似した関数を求めて、この関数を特性量修正関数Gpmとする(図12ステップS504)。ここで、モデル特性量SZを「暑い」という申告者率PPVhの飽和値とした本実施例では、SZ(i)=a(i)である。 The season-dependent management unit 21 sets i = N_ref, and when the processing of step S501 is completed for all the training data period LPs (i) of N_ref (YES in step S502), the function and parameter search of the individual thermal feeling model Find a function that linearly approximates the relationship between the model characteristic quantity SZ (i) determined from the modeling result Sp (i), which is a set of values, and the outdoor thermal environment representative value Eot (i), and use this function as the characteristic quantity correction function. Let it be Gpm (step S504 in FIG. 12). Here, in this embodiment in which the model characteristic quantity SZ is the saturation value of the reporter rate PPVh of “hot”, SZ (i) = a (i).

図14(A)、図14(B)は特性量修正関数Gpmの決定方法を説明する図である。図14(A)、図14(B)の例では、N_ref=3とし、最新の過去モデルの関数をMp(1)、最新から1回前の過去モデルの関数をMp(2)、最新から2回前の過去モデルの関数をMp(3)としている。 14 (A) and 14 (B) are diagrams illustrating a method for determining the characteristic quantity correction function Gpm. In the examples of FIGS. 14 (A) and 14 (B), N_ref = 3, the function of the latest past model is Mp (1), the function of the past model one time before the latest is Mp (2), and from the latest. The function of the past model two times ago is Mp (3).

例えば最新の過去モデルのモデリングの際に利用した学習データ期間LP(1)を9月、1回前の過去モデルのモデリングの際に利用した学習データ期間LP(2)を8月、2回前の過去モデルのモデリングの際に利用した学習データ期間LP(3)を7月とすれば、過去の屋外温熱環境代表値Eot(1),Eot(2),Eot(3)は、それぞれ9月、8月、7月の平年外気温である。また、モデル特性量SZ(1),SZ(2),SZ(3)は、それぞれ関数Mp(1),Mp(2),Mp(3)の飽和値(パラメータa(1),a(2),a(3))である。 For example, the training data period LP (1) used when modeling the latest past model is in September, and the training data period LP (2) used when modeling the previous past model is in August and twice before. Assuming that the learning data period LP (3) used in modeling the past model is July, the past outdoor thermal environment representative values Eot (1), Eot (2), and Eot (3) are September, respectively. , August and July are the average outside temperatures. Further, the model characteristic quantities SZ (1), SZ (2), and SZ (3) are saturated values of the functions Mp (1), Mp (2), and Mp (3) (parameters a (1) and a (2), respectively. ), A (3)).

そして、図14(B)に示すように、過去の屋外温熱環境代表値{Eot(1),Eot(2),Eot(3)}とモデル特性量{SZ(1),SZ(2),SZ(3)}との関係を線形近似した関数を特性量修正関数Gpmとする。特性量修正関数Gpmは、任意の月平年外気温から、これに対応するモデル特性量SZ(「暑い」という申告者率PPVhの飽和値)を推定する関数となる。 Then, as shown in FIG. 14 (B), past outdoor thermal environment representative values {Eot (1), Eot (2), Eot (3)} and model characteristic quantities {SZ (1), SZ (2), Let the characteristic quantity correction function Gpm be a function that linearly approximates the relationship with SZ (3)}. The characteristic quantity correction function Gpm is a function for estimating the model characteristic quantity SZ (saturation value of the reporter rate PPVh of "hot") corresponding to the arbitrary monthly outside temperature.

なお、本実施例では、特性量修正関数Gpmを求める手法として線形近似を用いたが、特性量修正関数Gpmは、多項式近似や指数近似などの手法により、温熱実感への季節推移の影響を考慮して適宜決定すればよい。 In this embodiment, linear approximation is used as a method for obtaining the characteristic quantity correction function Gpm, but the characteristic quantity correction function Gpm considers the influence of the seasonal transition on the feeling of heat by a method such as polynomial approximation or exponential approximation. And it may be decided appropriately.

なお、本実施例では、「暑い」という申告者率PPVhの飽和値をモデル特性量SZとして特性量修正関数Gpmを決定したが、モデル特性量SZは過去モデルの関数の形状を決定付ける別の量でも良く、また、モデル特性量SZは複数種類設定してもよい。 In this embodiment, the characteristic quantity correction function Gpm is determined by using the saturation value of the reporter rate PPVh of "hot" as the model characteristic quantity SZ, but the model characteristic quantity SZ is another that determines the shape of the function of the past model. The amount may be used, and a plurality of types of model characteristic quantity SZ may be set.

例えば「暑い」という申告者率PPVhの飽和値をモデル特性量SZ1(図13のSZ)、図15(A)に示すように申告者率PPVhの立上り位置のPMV値をモデル特性量SZ2、図15(B)に示すように申告者率PPVhが飽和値の50%に到達したときのPMV値をモデル特性量SZ3とする。 For example, the saturation value of the filer rate PPVh of "hot" is the model characteristic amount SZ1 (SZ in FIG. 13), and the PMV value of the rising position of the filer rate PPVh is the model characteristic amount SZ2 as shown in FIG. 15 (A). As shown in 15 (B), the PMV value when the filer rate PPVh reaches 50% of the saturation value is defined as the model characteristic quantity SZ3.

図15(A)の例では、モデル特性量SZ2は、個別温熱実感モデルを式(15)、式(16)とした場合、モデル探索パラメータa,b,c,d,eのうちのeに対応する。
したがって、季節依存管理部21は、N_ref回のモデリング結果Sp(i)のそれぞれからモデル特性量SZ2(i)=e(i)を取得し、モデル特性量SZ1の場合と同様に、モデル特性量SZ2(i)=e(i)と過去の屋外温熱環境代表値Eot(i)との関係を線形近似した関数を求めて、この関数を特性量修正関数Gpm_2とすればよい(ステップS504)。
In the example of FIG. 15 (A), the model characteristic quantity SZ2 is set to e of the model search parameters a, b, c, d, and e when the individual thermal feeling models are the equations (15) and (16). Correspond.
Therefore, the season-dependent management unit 21 acquires the model characteristic quantity SZ2 (i) = e (i) from each of the modeling results Sp (i) of N_ref times, and similarly to the case of the model characteristic quantity SZ1, the model characteristic quantity. A function that linearly approximates the relationship between SZ2 (i) = e (i) and the past outdoor thermal environment representative value Eot (i) may be obtained, and this function may be set as the characteristic quantity correction function Gpm_2 (step S504).

また、図15(B)の例では、モデル特性量SZ3は、申告者率PPVhが飽和値の50%に到達したときのPMV値である。したがって、季節依存管理部21は、個別温熱実感モデルのモデル関数と、N_ref回のモデリング結果Sp(i)の各々からa(i)の50%の値a_50(i)を算出し、「暑い」という申告者率PPVhがa_50(i)となるPMV値を、式(16)(または式(15))とモデリング結果Sp(i)とから算出して、算出した結果をSZ3(i)とすればよい。 Further, in the example of FIG. 15B, the model characteristic quantity SZ3 is a PMV value when the filer rate PPVh reaches 50% of the saturation value. Therefore, the season-dependent management unit 21 calculates the value a_50 (i) of 50% of a (i) from each of the model function of the individual thermal feeling model and the modeling result Sp (i) of N_ref times, and is "hot". The PMV value at which the reporter rate PPVh is a_50 (i) is calculated from the equation (16) (or the equation (15)) and the modeling result Sp (i), and the calculated result is referred to as SZ3 (i). Just do it.

そして、季節依存管理部21は、モデル特性量SZ1の場合と同様に、モデル特性量SZ3(i)と過去の屋外温熱環境代表値Eot(i)との関係を線形近似した関数を求めて、この関数を特性量修正関数Gpm_3とすればよい(ステップS504)。上記のモデル特性量SZ1(図13のSZ)から求めた特性量修正関数をGpm_1とする。
こうして、モデル特性量SZの種類毎に特性量修正関数を求めることができる。
Then, the seasonal dependence management unit 21 obtains a function that linearly approximates the relationship between the model characteristic quantity SZ3 (i) and the past outdoor thermal environment representative value Eot (i), as in the case of the model characteristic quantity SZ1. This function may be a characteristic quantity correction function Gpm_3 (step S504). Let Gpm_1 be the characteristic quantity correction function obtained from the model characteristic quantity SZ1 (SZ in FIG. 13).
In this way, the characteristic quantity correction function can be obtained for each type of model characteristic quantity SZ.

最後に、季節依存管理部21は、未来モデル生成部22からの情報送信要求に対して、特性量修正関数Gpm(Gpm_1,Gpm_2,Gpm_3)と最新の過去モデルのモデリング結果とを送信する(図11ステップS404)。 Finally, the seasonal dependence management unit 21 transmits the characteristic quantity correction function Gpm (Gpm_1, Gpm_2, Gpm_3) and the modeling result of the latest past model in response to the information transmission request from the future model generation unit 22 (Fig.). 11 steps S404).

なお、参照モデル数N_ref、過去の屋外温熱環境代表値Eotの決定方法、モデル特性量SZの決定方法、および各特性量に対応する特性量修正関数Gpmの決定方法は、室内環境評価/制御等のソリューションプロバイダや設備管理者によって予め設定されている。 The reference model number N_ref, the method for determining the past outdoor thermal environment representative value Eot, the method for determining the model characteristic quantity SZ, and the method for determining the characteristic quantity correction function Gpm corresponding to each characteristic quantity include indoor environment evaluation / control, etc. It is preset by the solution provider and equipment manager of.

図16は未来モデル生成部22の動作を説明するフローチャートである。未来モデル生成部22は、評価対象空間101の居住者の温熱実感(評価対象空間101の温熱環境に対する不満足を申告する居住者の割合であり、「暑い」という申告者率PPVh)を推定したい未来の日時(推定未来日時Tm)をオペレータからの指示などで取得する(図16ステップS600)。 FIG. 16 is a flowchart illustrating the operation of the future model generation unit 22. The future model generation unit 22 wants to estimate the feeling of heat of the residents of the evaluation target space 101 (the ratio of the residents who report dissatisfaction with the thermal environment of the evaluation target space 101, and the reporter rate PPVh of "hot"). (Estimated future date and time Tm) is acquired by an instruction from the operator or the like (FIG. 16 step S600).

未来モデル生成部22は、データ期間管理部14から個別温熱実感モデルの更新タイミングの情報を取得し、推定未来日時Tmを含む未来の学習データ期間LPFを決定する(図16ステップS601)。例えば、推定未来日時Tmを2019年10月15日とすると、個別温熱実感モデルの更新タイミングが毎月1日の0時00分、学習データ期間幅が1か月であるので、未来モデル生成部22は、未来の学習データ期間LPFを2019年10月と決定する。 The future model generation unit 22 acquires information on the update timing of the individual thermal realization model from the data period management unit 14, and determines the future learning data period LPF including the estimated future date and time Tm (step S601 in FIG. 16). For example, assuming that the estimated future date and time Tm is October 15, 2019, the update timing of the individual thermal realization model is 0:00 on the first day of every month, and the training data period width is one month. Therefore, the future model generation unit 22 Determines the future learning data period LPF as October 2019.

そして、未来モデル生成部22は、屋外温熱環境情報管理部20から、未来の学習データ期間LPFに対応する未来の屋外温熱環境情報を取得し(図16ステップS602)、取得した未来の屋外温熱環境情報に基づいて、未来の学習データ期間LPFに対応する未来の屋外温熱環境代表値Eotfを決定する(図16ステップS603)。例えば、未来の学習データ期間LPFが2019年10月の場合、未来モデル生成部22は、10月の平年月平均外気温を取得して未来の屋外温熱環境代表値Eotfとすればよい。 Then, the future model generation unit 22 acquires the future outdoor thermal environment information corresponding to the future learning data period LPF from the outdoor thermal environment information management unit 20 (step S602 in FIG. 16), and the acquired future outdoor thermal environment. Based on the information, the future outdoor thermal environment representative value Eotf corresponding to the future learning data period LPF is determined (FIG. 16 step S603). For example, when the future learning data period LPF is October 2019, the future model generation unit 22 may acquire the average outside temperature of the average month of October and set it as the future outdoor thermal environment representative value Etf.

次に、未来モデル生成部22は、季節依存管理部21から送信された特性量修正関数Gpm(ここではGpm_1)を利用して、未来の屋外温熱環境代表値Eotfに対応するモデル特性量SZ1=Gpm_1(Eoft)を算出し、算出した値を推定特性量SZestとする(図16ステップS604)。 Next, the future model generation unit 22 uses the characteristic quantity correction function Gpm (here, Gpm_1) transmitted from the seasonal dependence management unit 21 to model characteristic quantity SZ1 = corresponding to the future outdoor thermal environment representative value Eotf. Gpm_1 (Eoft) is calculated, and the calculated value is used as the estimated characteristic quantity SZest (step S604 in FIG. 16).

最後に、未来モデル生成部22は、算出した推定特性量SZestに基づいて最新の過去モデルを修正した個別温熱実感モデル(未来モデル)を生成する(図16ステップS605)。過去モデルを修正する方法は、室内環境評価/制御等のソリューションプロバイダや設備管理者によって予め設定されている。 Finally, the future model generation unit 22 generates an individual thermal realization model (future model) in which the latest past model is modified based on the calculated estimated characteristic quantity SZest (step S605 in FIG. 16). The method of modifying the past model is preset by the solution provider such as indoor environment evaluation / control and the equipment manager.

本実施例では、推定特性量SZestと、最新の過去モデルのモデル特性量SZ1(1)=a(1)との比率をモデル補正パラメータηとする。
η=SZest/SZ1(1) ・・・(17)
In this embodiment, the ratio between the estimated characteristic quantity SZest and the model characteristic quantity SZ1 (1) = a (1) of the latest past model is set as the model correction parameter η.
η = SZest / SZ1 (1) ... (17)

そして、未来モデル生成部22は、モデル補正パラメータηに基づいて、最新の過去モデルの関数Mp(1)を次式のように修正する。 Then, the future model generation unit 22 modifies the function Mp (1) of the latest past model based on the model correction parameter η as follows.

式(18)は、最新の過去モデル(式(16))のパラメータa(1),c(1),d(1),e(1)のうち、パラメータa(1)をモデル補正パラメータηに基づいて修正することを示している。 In the equation (18), among the parameters a (1), c (1), d (1), and e (1) of the latest past model (equation (16)), the parameter a (1) is used as the model correction parameter η. It is shown to be corrected based on.

なお、本実施例の過去モデルの修正では、補正パラメータηによりPPVh軸についての倍率修正のみを行っているが、上記のモデル特性量SZ2,SZ3を用いることにより、申告者率PPVhの立上り位置の修正やPMV軸の倍率修正などを行うようにしてもよい。 In the modification of the past model of this embodiment, only the magnification of the PPVh axis is modified by the correction parameter η, but by using the above model characteristic quantities SZ2 and SZ3, the rising position of the filer rate PPVh The correction or the magnification correction of the PMV axis may be performed.

具体的には、未来モデル生成部22は、季節依存管理部21から送信された特性量修正関数Gpm_2を利用して、未来の屋外温熱環境代表値Eotfに対応するモデル特性量SZ2=Gpm_2(Eoft)を算出し、算出した値を推定特性量SZest_2とする。また、未来モデル生成部22は、特性量修正関数Gpm_3を利用して、未来の屋外温熱環境代表値Eotfに対応するモデル特性量SZ3=Gpm_3(Eoft)を算出し、算出した値を推定特性量SZest_3とする(図16ステップS604)。 Specifically, the future model generation unit 22 uses the characteristic quantity correction function Gpm_2 transmitted from the seasonal dependence management unit 21 to model characteristic quantity SZ2 = Gpm_2 (Eoft) corresponding to the future outdoor thermal environment representative value Eotf. ) Is calculated, and the calculated value is used as the estimated characteristic quantity SZest_2. Further, the future model generation unit 22 calculates the model characteristic quantity SZ3 = Gpm_3 (Eoft) corresponding to the future outdoor thermal environment representative value Eotf by using the characteristic quantity correction function Gpm_3, and estimates the calculated value. It is set to SZest_3 (step S604 in FIG. 16).

そして例えば、推定特性量SZest_2と、最新の過去モデルのモデル特性量SZ2(1)=e(1)との差をモデル補正パラメータδとする。
δ=SZest_2−SZ2(1) ・・・(19)
Then, for example, the difference between the estimated characteristic quantity SZest_2 and the model characteristic quantity SZ2 (1) = e (1) of the latest past model is set as the model correction parameter δ.
δ = SZest_2-SZ2 (1) ... (19)

さらに例えば、推定特性量SZest_3とSZest_2との差(SZest_3−SZest_2)と、最新の過去モデルのモデル特性量SZ3(1)とSZ2(1)との差(SZ3(1)−SZ2(1))との比率をモデル補正パラメータξとする。
ξ=(SZest_3−SZest_2)/(SZ3(1)−SZ2(1))
・・・(20)
Further, for example, the difference between the estimated characteristic quantities SZest_3 and SZest_2 (SZest_3-SZest_2) and the difference between the model characteristic quantities SZ3 (1) and SZ2 (1) of the latest past model (SZ3 (1) -SZ2 (1)). Let the ratio with and be the model correction parameter ξ.
ξ = (SZest_3-SZest_2) / (SZ3 (1) -SZ2 (1))
... (20)

そして、未来モデル生成部22は、モデル補正パラメータη,δ,ξに基づいて、最新の過去モデルの関数Mp(1)を次式のように修正する。 Then, the future model generation unit 22 modifies the function Mp (1) of the latest past model based on the model correction parameters η, δ, ξ as follows.

式(21)のように過去モデルの修正を行えば、モデル補正パラメータδによる申告者率PPVhの立上り位置の修正や、モデル補正パラメータξによるPMV軸の倍率修正が可能である。 By modifying the past model as in equation (21), it is possible to modify the rising position of the filer rate PPVh by the model correction parameter δ and to modify the magnification of the PMV axis by the model correction parameter ξ.

以上により、未来モデル生成部22の処理が終了する。本実施例によれば、季節推移が反映された屋外温熱環境状態に基づいて過去モデルを修正して未来モデルを生成することにより、季節推移要素による推定誤差を改善することができ、未来モデルに推定未来日時の推定PMV値を入力してやれば、推定未来日時における評価対象空間101に対する居住者の温熱実感(評価対象空間101の温熱環境に対する不満足を申告する居住者の割合であり、「暑い」という申告者率PPVh)を推定することが可能となる。その結果、本実施例では、居住者の温熱実感の推定結果を実態に近づけることができる。 As a result, the processing of the future model generation unit 22 is completed. According to this embodiment, by modifying the past model based on the outdoor thermal environment condition reflecting the seasonal transition to generate the future model, the estimation error due to the seasonal transition element can be improved, and the future model can be used. If the estimated PMV value of the estimated future date and time is input, the resident's feeling of warmth with respect to the evaluation target space 101 at the estimated future date and time (the ratio of the resident who declares dissatisfaction with the thermal environment of the evaluation target space 101, which is called "hot". It is possible to estimate the filer rate PPVh). As a result, in this embodiment, the estimation result of the resident's feeling of heat can be brought closer to the actual situation.

[第2の実施例]
次に、本発明の第2の実施例について説明する。本実施例においても、温熱実感モデル生成装置の構成は第1の実施例と同様であるので、図2の符号を用いて説明する。
短期的な変動が含まれる数日程度の期間に対応する学習データを使用すると、モデルの信頼性が保てないため、学習データ(統合情報)として一定の期間のデータを確保する必要がある。第1の実施例で説明したとおり、学習データ期間幅としては1か月程度が適当であるが、さらに季節推移の連続性を維持するため、個々の学習データ期間が重なり合うオーバーラップ期間を設ける。本実施例では、信頼できる一定の期間としてオーバーラップ期間を2週間とした例で説明する。
[Second Example]
Next, a second embodiment of the present invention will be described. Also in this embodiment, the configuration of the thermal sensation model generator is the same as that in the first embodiment, and thus the reference numerals in FIG. 2 will be used for description.
If training data corresponding to a period of several days including short-term fluctuations is used, the reliability of the model cannot be maintained, so it is necessary to secure data for a certain period as training data (integrated information). As described in the first embodiment, the learning data period width is appropriately about one month, but in order to maintain the continuity of the seasonal transition, an overlapping period in which the individual learning data periods overlap is provided. In this embodiment, an example in which the overlap period is set to 2 weeks as a reliable constant period will be described.

本実施例の温熱実感モデル生成装置のデータ期間管理部14には、第1の実施例と同様に個別温熱実感モデル(過去モデル)の更新タイミングと学習データ期間とを決定するためのデータ期間管理情報が室内環境評価/制御等のソリューションプロバイダや設備管理者によって予め設定されている。 The data period management unit 14 of the thermal sensation model generator of this embodiment has the data period management for determining the update timing and the learning data period of the individual thermal sensation model (past model) as in the first embodiment. Information is preset by solution providers such as indoor environment evaluation / control and equipment managers.

第1の実施例と異なる点は、データ期間管理情報が個別温熱実感モデルの更新タイミングの情報と学習データ期間幅の情報の他に、オーバーラップ期間幅Loの情報を含むことである。 The difference from the first embodiment is that the data period management information includes the information of the overlap period width Lo in addition to the information of the update timing of the individual thermal feeling model and the information of the learning data period width.

図17(A)、図17(B)は第1、第2の実施例の学習データ期間を説明する図である。なお、図17(A)、図17(B)では、分かり易い説明とするため、各月を一律に30日として更新タイミングの日にちを示しており、2018年10月1日を最新の更新タイミングとして2018年8月1日の更新タイミングまで遡った各更新タイミングの学習データ期間を示している。また、図17(A)、図17(B)いずれも、更新スタート日時を2018年8月1日0時00分、学習データの期間幅を30日毎とし、図17(B)の例では、本実施例のオーバーラップ期間幅Loを15日としている。図17(A)のような連続する更新タイミングの学習データ期間に対して、オーバーラップ期間幅Loを設定する本実施例の図17(B)では、モデル更新の周期がオーバーラップ期間幅Loに対応して短くなる。 17 (A) and 17 (B) are diagrams for explaining the learning data period of the first and second embodiments. In addition, in FIGS. 17 (A) and 17 (B), the update timing dates are uniformly set to 30 days for each month for easy understanding, and October 1, 2018 is the latest update timing. The learning data period of each update timing that goes back to the update timing of August 1, 2018 is shown. Further, in both FIGS. 17 (A) and 17 (B), the update start date and time is set to 0:00 on August 1, 2018, and the period width of the learning data is set to every 30 days. In the example of FIG. 17 (B), The overlap period width Lo of this embodiment is set to 15 days. In FIG. 17 (B) of the present embodiment in which the overlap period width Lo is set for the learning data period of continuous update timing as shown in FIG. 17 (A), the model update cycle is set to the overlap period width Lo. Correspondingly shortened.

第1の実施例に対応する本説明の例では、図17(A)に示すように、例えば2018年9月1日0時00分に更新される過去モデルMp(2)の学習データ期間SWは、8月1日0時00分を開始タイミングとして、8月30日(ここでは各月を一律に30日としているので8月末日に相当)23時59分を終了タイミングとする30日である。さらに、ひとつまえの過去モデルMp(3)の学習データ期間SWは、7月1日0時00分を開始タイミングとして、7月30日(ここでは各月を一律に30日としているので7月末日に相当)23時59分を終了タイミングとする30日であり、連続する更新タイミングの過去モデルMp(2)とMp(3)の学習データ期間SWは重複しない。 In the example of this description corresponding to the first embodiment, as shown in FIG. 17 (A), for example, the training data period SW of the past model Mp (2) updated at 0:00 on September 1, 2018. Is 30 days with the start timing of 0:00 on August 1st and the end timing of 23:59 on August 30th (here, since each month is uniformly set to 30th, it corresponds to the end of August). is there. Furthermore, the learning data period SW of the previous past model Mp (3) starts at 0:00 on July 1st and is July 30th (here, each month is uniformly 30th, so the end of July). It is 30 days with 23:59 as the end timing (corresponding to the day), and the learning data period SWs of the past models Mp (2) and Mp (3) of the continuous update timing do not overlap.

一方、本実施例では、図17(B)に示すように、例えば2018年9月1日0時00分に更新される過去モデルMp(3)の学習データ期間SWは、8月1日0時00分を開始タイミングとして、8月30日(ここでは各月を一律に30日としているので8月末日に相当)23時59分を終了タイミングとする30日であり、図17(A)の過去モデルMp(2)の学習データ期間SWと同じであるが、図17(B)のひとつまえの過去モデルMp(4)の学習データ期間SWは、7月16日0時00分を開始タイミングとして、図17(B)の過去モデルMp(3)の学習データ期間SWと重複するオーバーラップ期間幅Loを含み(すなわち、8月1日から8月15日の15日間)、8月15日23時59分を終了タイミングとする30日であり、連続する更新タイミングの過去モデルMp(3)とMp(4)の学習データ期間SWはオーバーラップ期間幅Loだけ重複する。 On the other hand, in this embodiment, as shown in FIG. 17B, for example, the training data period SW of the past model Mp (3) updated at 0:00 on September 1, 2018 is 0 on August 1. The start timing is 0:00, and the end timing is August 30 (corresponding to the end of August because each month is uniformly set to 30 days), which is the 30th with the end timing shown in FIG. 17 (A). The training data period SW of the past model Mp (2) is the same as that of the past model Mp (2), but the training data period SW of the previous model Mp (4) in FIG. 17 (B) starts at 0:00 on July 16. The timing includes an overlap period width Lo that overlaps with the training data period SW of the past model Mp (3) of FIG. 17 (B) (that is, 15 days from August 1 to August 15), and August 15 It is 30 days with the end timing at 23:59 on the day, and the learning data period SWs of the past models Mp (3) and Mp (4) of the continuous update timing overlap by the overlap period width Lo.

その他の構成は第1の実施例で説明したとおりである。なお、本実施例の学習データ期間幅としては1か月程度が適当であるが、この時、さらに、季節推移の連続性を維持するため、信頼できる期間幅を前後の学習データ期間でオーバーラップさせることが望ましい。オーバーラップ期間幅Loとしては、2週間以上の長い程度の期間とするのが適当である。 Other configurations are as described in the first embodiment. It should be noted that the learning data period width of this embodiment is appropriately about one month, but at this time, in order to maintain the continuity of the seasonal transition, the reliable period width overlaps with the learning data period before and after. It is desirable to let it. It is appropriate that the overlap period width Lo is a long period of 2 weeks or more.

ある建物で2017年7月から9月にかけて計測した外気温の、2週間平均の推移を図18(A)に、4週間平均の推移を図18(B)に示す。図18(A)、図18(B)の横軸は1週目からカウントした経過週の数である。 The two-week average transition of the outside air temperature measured from July to September 2017 in a certain building is shown in FIG. 18 (A), and the four-week average transition is shown in FIG. 18 (B). The horizontal axis of FIGS. 18A and 18B is the number of elapsed weeks counted from the first week.

図18(A)の例ではオーバーラップ期間を1週間、図18(B)の例ではオーバーラップ期間を2週間としている。図18(A)、図18(B)の比較から明らかなように、オーバーラップ期間を2週間として、4週間の平均値をとると、外気温平均値に短期的な変動が含まれ難いことが分かる。 In the example of FIG. 18 (A), the overlap period is one week, and in the example of FIG. 18 (B), the overlap period is two weeks. As is clear from the comparison of FIGS. 18 (A) and 18 (B), when the overlap period is 2 weeks and the average value of 4 weeks is taken, it is difficult for the average value of the outside air temperature to include short-term fluctuations. I understand.

第1、第2の実施例で説明した温熱実感モデル生成装置は、CPU(Central Processing Unit)、記憶装置及びインターフェースを備えたコンピュータと、これらのハードウェア資源を制御するプログラムによって実現することができる。このコンピュータの構成例を図19に示す。 The thermal sensation model generation device described in the first and second embodiments can be realized by a computer provided with a CPU (Central Processing Unit), a storage device, and an interface, and a program for controlling these hardware resources. .. A configuration example of this computer is shown in FIG.

コンピュータは、CPU200と、記憶装置201と、インターフェース装置(以下、I/Fと略する)202とを備えている。I/F202には、環境計測デバイス等が接続される。このようなコンピュータにおいて、本発明の温熱実感モデル生成方法を実現させるためのプログラムは記憶装置201に格納される。CPU200は、記憶装置201に格納されたプログラムに従って第1、第2の実施例で説明した処理を実行する。 The computer includes a CPU 200, a storage device 201, and an interface device (hereinafter, abbreviated as I / F) 202. An environment measurement device or the like is connected to the I / F 202. In such a computer, a program for realizing the method for generating a thermal sensation model of the present invention is stored in the storage device 201. The CPU 200 executes the processes described in the first and second embodiments according to the program stored in the storage device 201.

本発明は、未来の日時における居住者の温熱実感を推定する技術に適用することができる。 The present invention can be applied to a technique for estimating a resident's feeling of heat at a future date and time.

1…過去モデル生成ユニット、2…未来モデル生成ユニット、10…データ収集部、11…居住環境情報保持部、12…温熱実感情報保持部、13…データ統合部、14…データ期間管理部、15…過去モデル生成部、16…過去モデル記憶管理部、20…屋外温熱環境情報管理部、21…季節依存管理部、22…未来モデル生成部。 1 ... Past model generation unit, 2 ... Future model generation unit, 10 ... Data collection unit, 11 ... Living environment information holding unit, 12 ... Thermal feeling information holding unit, 13 ... Data integration unit, 14 ... Data period management department, 15 ... Past model generation unit, 16 ... Past model memory management department, 20 ... Outdoor thermal environment information management department, 21 ... Seasonal dependence management department, 22 ... Future model generation department.

Claims (10)

評価対象空間の居住環境指標と前記評価対象空間に対する居住者の温熱実感との関係をモデル化した第1のモデルの更新タイミングの情報と学習データ期間の情報とを記憶するように構成されたデータ期間管理部と、
前記更新タイミングになったときに、前記学習データ期間における前記居住環境指標と前記居住者の温熱実感情報とに基づいて前記第1のモデルを生成するように構成された第1のモデル生成部と、
複数の学習データ期間によってそれぞれ生成された複数の前記第1のモデルと前記複数の学習データ期間のそれぞれに対応する過去の屋外温熱環境の代表値とに基づいて、前記第1のモデルの関数の形状を決定付けるモデル特性量と屋外温熱環境の代表値との関係を示す特性量修正関数を求めるように構成された季節依存管理部と、
前記居住者の温熱実感を推定したい未来の日時に対して設定されるべき未来の学習データ期間を決定して、この学習データ期間における屋外温熱環境の代表値に対応するモデル特性量を前記特性量修正関数に基づいて算出して、算出したモデル特性量に基づいて最新の前記第1のモデルを修正した第2のモデルを生成するように構成された第2のモデル生成部とを備えることを特徴とする温熱実感モデル生成装置。
Data configured to store the update timing information and the learning data period information of the first model that models the relationship between the living environment index of the evaluation target space and the resident's feeling of heat with respect to the evaluation target space. With the period management department
When the update timing comes, the first model generation unit configured to generate the first model based on the living environment index in the learning data period and the thermal feeling information of the resident. ,
The function of the first model is based on the plurality of first models generated by the plurality of training data periods and the representative values of the past outdoor thermal environments corresponding to each of the plurality of training data periods. A seasonally dependent management unit configured to obtain a characteristic quantity correction function that indicates the relationship between the model characteristic quantity that determines the shape and the representative value of the outdoor thermal environment.
The future learning data period to be set for the future date and time for which the resident's feeling of heat is to be estimated is determined, and the model characteristic amount corresponding to the representative value of the outdoor thermal environment in this learning data period is defined as the characteristic amount. It is provided with a second model generation unit configured to generate a second model obtained by calculating based on a modification function and modifying the latest first model based on the calculated model characteristic quantity. A featured thermal sensation model generator.
請求項1記載の温熱実感モデル生成装置において、
前記データ期間管理部は、前記更新タイミングの情報と前記学習データ期間の情報に加えて、オーバーラップ期間の情報を記憶し、
前記複数の学習データ期間は、連続する前後の学習データ期間が前記オーバーラップ期間の幅の分だけ重なるように設定されることを特徴とする温熱実感モデル生成装置。
In the thermal feeling model generator according to claim 1,
The data period management unit stores information on the overlap period in addition to the information on the update timing and the information on the learning data period.
The plurality of learning data periods are set so that the learning data periods before and after the continuation overlap by the width of the overlapping period.
請求項1または2記載の温熱実感モデル生成装置において、
前記評価対象空間が含まれる建物の過去の屋外温熱環境情報を記憶するように構成された屋外温熱環境情報管理部をさらに備え、
前記季節依存管理部は、前記複数の学習データ期間のそれぞれに対応する過去の前記屋外温熱環境情報を前記屋外温熱環境情報管理部から取得し、前記複数の学習データ期間のそれぞれに対応する過去の屋外温熱環境の代表値を決定して、前記特性量修正関数を求め、
前記第2のモデル生成部は、前記未来の学習データ期間に対応する過去の前記屋外温熱環境情報を前記屋外温熱環境情報管理部から取得して、前記未来の学習データ期間に対応する屋外温熱環境の代表値を決定することを特徴とする温熱実感モデル生成装置。
In the thermal feeling model generator according to claim 1 or 2.
It is further equipped with an outdoor thermal environment information management unit configured to store past outdoor thermal environment information of the building including the evaluation target space.
The seasonally dependent management unit acquires the past outdoor thermal environment information corresponding to each of the plurality of learning data periods from the outdoor thermal environment information management unit, and the past corresponding to each of the plurality of learning data periods. The representative value of the outdoor thermal environment was determined, and the characteristic quantity correction function was obtained.
The second model generation unit acquires the past outdoor thermal environment information corresponding to the future learning data period from the outdoor thermal environment information management unit, and obtains the outdoor thermal environment corresponding to the future learning data period. A thermal sensation model generator characterized by determining a representative value of.
請求項1乃至3のいずれか1項に記載の温熱実感モデル生成装置において、
前記第2のモデル生成部は、前記特性量修正関数に基づいて算出したモデル特性量と最新の前記第1のモデルのモデル特性量とに基づいて、最新の前記第1のモデルを修正することを特徴とする温熱実感モデル生成装置。
In the thermal feeling model generator according to any one of claims 1 to 3.
The second model generation unit modifies the latest first model based on the model characteristic quantity calculated based on the characteristic quantity correction function and the latest model characteristic quantity of the first model. A thermal sensation model generator featuring.
請求項1乃至4のいずれか1項に記載の温熱実感モデル生成装置において、
前記居住環境指標は、PMVであり、前記居住者の温熱実感は、温熱環境に対する不満足を申告する居住者の割合であることを特徴とする温熱実感モデル生成装置。
In the thermal feeling model generator according to any one of claims 1 to 4.
The living environment index is PMV, and the resident's feeling of heat is the proportion of the resident who declares dissatisfaction with the heating environment.
評価対象空間の居住環境指標と前記評価対象空間に対する居住者の温熱実感との関係をモデル化した第1のモデルの更新タイミングの情報と学習データ期間の情報とを記憶するデータ期間管理部を参照し、前記更新タイミングになったときに、前記学習データ期間における前記居住環境指標と前記居住者の温熱実感情報とに基づいて前記第1のモデルを生成する第1のステップと、
複数の学習データ期間によってそれぞれ生成された複数の前記第1のモデルと前記複数の学習データ期間のそれぞれに対応する過去の屋外温熱環境の代表値とに基づいて、前記第1のモデルの関数の形状を決定付けるモデル特性量と屋外温熱環境の代表値との関係を示す特性量修正関数を求める第2のステップと、
前記居住者の温熱実感を推定したい未来の日時に対して設定されるべき未来の学習データ期間を決定して、この学習データ期間における屋外温熱環境の代表値に対応するモデル特性量を前記特性量修正関数に基づいて算出して、算出したモデル特性量に基づいて最新の前記第1のモデルを修正した第2のモデルを生成する第3のステップとを含むことを特徴とする温熱実感モデル生成方法。
Refer to the data period management unit that stores the update timing information and the learning data period information of the first model that models the relationship between the living environment index of the evaluation target space and the resident's feeling of heat with respect to the evaluation target space. Then, when the update timing is reached, the first step of generating the first model based on the living environment index in the learning data period and the thermal feeling information of the resident,
The function of the first model is based on the plurality of first models generated by the plurality of training data periods and the representative values of the past outdoor thermal environments corresponding to each of the plurality of training data periods. The second step of finding the characteristic quantity correction function showing the relationship between the model characteristic quantity that determines the shape and the representative value of the outdoor thermal environment, and
The future learning data period to be set for the future date and time for which the resident's feeling of heat is to be estimated is determined, and the model characteristic amount corresponding to the representative value of the outdoor thermal environment in this learning data period is defined as the characteristic amount. Thermal sensation model generation, which includes a third step of calculating based on a modification function and generating a second model modified from the latest first model based on the calculated model characteristic quantity. Method.
請求項6記載の温熱実感モデル生成方法において、
前記データ期間管理部は、前記更新タイミングの情報と前記学習データ期間の情報に加えて、オーバーラップ期間の情報を記憶し、
前記複数の学習データ期間は、連続する前後の学習データ期間が前記オーバーラップ期間の幅の分だけ重なるように設定されることを特徴とする温熱実感モデル生成方法。
In the method for generating a thermal feeling model according to claim 6,
The data period management unit stores information on the overlap period in addition to the information on the update timing and the information on the learning data period.
The method for generating a thermal sensation model, wherein the plurality of learning data periods are set so that the learning data periods before and after the continuation overlap by the width of the overlapping period.
請求項6または7記載の温熱実感モデル生成方法において、
前記第2のステップは、前記評価対象空間が含まれる建物の過去の屋外温熱環境情報を記憶する屋外温熱環境情報管理部を参照し、前記複数の学習データ期間のそれぞれに対応する過去の前記屋外温熱環境情報を前記屋外温熱環境情報管理部から取得し、前記複数の学習データ期間のそれぞれに対応する過去の屋外温熱環境の代表値を決定して、前記特性量修正関数を求めるステップを含み、
前記第3のステップは、前記未来の学習データ期間に対応する過去の前記屋外温熱環境情報を前記屋外温熱環境情報管理部から取得して、前記未来の学習データ期間に対応する屋外温熱環境の代表値を決定するステップを含むことを特徴とする温熱実感モデル生成方法。
In the method for generating a thermal sensation model according to claim 6 or 7.
The second step refers to the outdoor thermal environment information management unit that stores the past outdoor thermal environment information of the building including the evaluation target space, and refers to the past outdoor thermal environment corresponding to each of the plurality of learning data periods. Including the step of acquiring the thermal environment information from the outdoor thermal environment information management unit, determining the representative value of the past outdoor thermal environment corresponding to each of the plurality of learning data periods, and obtaining the characteristic quantity correction function.
In the third step, the past outdoor thermal environment information corresponding to the future learning data period is acquired from the outdoor thermal environment information management unit, and the representative of the outdoor thermal environment corresponding to the future learning data period is obtained. A method for generating a thermal sensation model, which comprises a step of determining a value.
請求項6乃至8のいずれか1項に記載の温熱実感モデル生成方法において、
前記第3のステップは、前記特性量修正関数に基づいて算出したモデル特性量と最新の前記第1のモデルのモデル特性量とに基づいて、最新の前記第1のモデルを修正するステップを含むことを特徴とする温熱実感モデル生成方法。
In the method for generating a thermal sensation model according to any one of claims 6 to 8.
The third step includes a step of modifying the latest first model based on the model characteristic quantity calculated based on the characteristic quantity correction function and the latest model characteristic quantity of the first model. A method for generating a thermal sensation model, which is characterized in that.
請求項6乃至9のいずれか1項に記載の温熱実感モデル生成方法において、
前記居住環境指標は、PMVであり、前記居住者の温熱実感は、温熱環境に対する不満足を申告する居住者の割合であることを特徴とする温熱実感モデル生成方法。
In the method for generating a thermal sensation model according to any one of claims 6 to 9.
The method for generating a thermal sensation model, wherein the living environment index is PMV, and the resident's thermal sensation is the proportion of resident who declares dissatisfaction with the thermal environment.
JP2019118176A 2019-06-26 2019-06-26 Thermal sensation model generation device and method Active JP7272877B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2019118176A JP7272877B2 (en) 2019-06-26 2019-06-26 Thermal sensation model generation device and method
SG10202005263SA SG10202005263SA (en) 2019-06-26 2020-06-04 Thermal sensation model generating apparatus and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019118176A JP7272877B2 (en) 2019-06-26 2019-06-26 Thermal sensation model generation device and method

Publications (2)

Publication Number Publication Date
JP2021004691A true JP2021004691A (en) 2021-01-14
JP7272877B2 JP7272877B2 (en) 2023-05-12

Family

ID=74097073

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019118176A Active JP7272877B2 (en) 2019-06-26 2019-06-26 Thermal sensation model generation device and method

Country Status (2)

Country Link
JP (1) JP7272877B2 (en)
SG (1) SG10202005263SA (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06207734A (en) * 1993-01-08 1994-07-26 Matsushita Refrig Co Ltd Air conditioner
JP2010266975A (en) * 2009-05-13 2010-11-25 Sony Corp Learning device and method, data generating device and method, and program
US20120259469A1 (en) * 2009-12-16 2012-10-11 John Ward Hvac control system and method
JP2018062190A (en) * 2016-10-11 2018-04-19 公益財団法人鉄道総合技術研究所 In-vehicle air-conditioning method and system

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104394044B (en) 2014-10-29 2018-02-02 小米科技有限责任公司 The method and apparatus of self-defined smart machine scene mode

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH06207734A (en) * 1993-01-08 1994-07-26 Matsushita Refrig Co Ltd Air conditioner
JP2010266975A (en) * 2009-05-13 2010-11-25 Sony Corp Learning device and method, data generating device and method, and program
US20120259469A1 (en) * 2009-12-16 2012-10-11 John Ward Hvac control system and method
JP2013514510A (en) * 2009-12-16 2013-04-25 コモンウェルス サイエンティフィック アンドインダストリアル リサーチ オーガナイゼーション HVAC control system and method
JP2018062190A (en) * 2016-10-11 2018-04-19 公益財団法人鉄道総合技術研究所 In-vehicle air-conditioning method and system

Also Published As

Publication number Publication date
SG10202005263SA (en) 2021-01-28
JP7272877B2 (en) 2023-05-12

Similar Documents

Publication Publication Date Title
Laftchiev et al. An IoT system to estimate personal thermal comfort
US20160274963A1 (en) Equipment evaluation device, equipment evaluation method and non-transitory computer readable medium
US11808473B2 (en) Action optimization device, method and program
JP6126380B2 (en) Living environment evaluation method and apparatus
CN109899933B (en) Cold and heat feeling report information processing device and method
JP5070307B2 (en) Basic unit calculation system, program for executing the system, and recording medium recording the program
JP6878361B2 (en) Insulation performance diagnostic device
Biyik et al. Cloud-based model predictive building thermostatic controls of commercial buildings: Algorithm and implementation
Ding et al. Influence of occupancy-oriented interior cooling load on building cooling load design
Love Understanding the interactions between occupants, heating systems and building fabric in the context of energy efficient building fabric retrofit in social housing
Wang et al. The virtual in-situ calibration of various physical sensors in air handling units
JP7272877B2 (en) Thermal sensation model generation device and method
EP2945105A1 (en) Device for estimating thermal characteristics of room, program
Zhou et al. A model-driven learning approach for predicting the personalized dynamic thermal comfort in ordinary office environment
Merabet et al. Measuring human comfort for smart building application: experimental set-up using WSN
JP2018155435A (en) Air conditioner control device and air conditioner control method
Tanimoto et al. Cooling load simulation considering actual variation of inhabitants’ behavior for accurate estimation of urban maximum energy requirement
JP7375355B2 (en) Estimation method, estimation device, computer program, and method for generating a trained model
Fraga et al. Optimisation as a tool for gaining insight: an application to the built environment
KR102324304B1 (en) Demand response Airconditioner control system and method based on demander&#39;s comfort
Carroll Index to quantify thermal comfort in homes
CN111928458B (en) Indoor environment temperature control method and device based on actual thermal comfort of user
Varela et al. Predicting historical indoor temperatures from available local weather data
Khamseh et al. How people use their thermostats? high-resolution data analysis of two years of system operations in office rooms
JP2021196616A (en) Management device and management method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220324

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20230206

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: 20230404

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230427

R150 Certificate of patent or registration of utility model

Ref document number: 7272877

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150