WO2023181778A1 - 推定方法、プログラム、及び推定システム - Google Patents

推定方法、プログラム、及び推定システム Download PDF

Info

Publication number
WO2023181778A1
WO2023181778A1 PCT/JP2023/006667 JP2023006667W WO2023181778A1 WO 2023181778 A1 WO2023181778 A1 WO 2023181778A1 JP 2023006667 W JP2023006667 W JP 2023006667W WO 2023181778 A1 WO2023181778 A1 WO 2023181778A1
Authority
WO
WIPO (PCT)
Prior art keywords
space
carbon dioxide
information
concentration distribution
estimation
Prior art date
Application number
PCT/JP2023/006667
Other languages
English (en)
French (fr)
Inventor
浩史 久保田
将大 鶴居
Original Assignee
パナソニックＩｐマネジメント株式会社
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 パナソニックＩｐマネジメント株式会社 filed Critical パナソニックＩｐマネジメント株式会社
Publication of WO2023181778A1 publication Critical patent/WO2023181778A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Definitions

  • the present disclosure generally relates to an estimation method, program, and estimation system, and more specifically relates to an estimation method, program, and estimation system regarding the concentration of carbon dioxide in an indoor space.
  • Patent Document 1 describes a wide area distribution monitoring device that measures space using a path configuration optical system and calculates the concentration of carbon dioxide in the space using a weight distribution function.
  • the present disclosure has been made in view of the above reasons, and aims to provide an estimation method, program, and estimation system that can estimate the carbon dioxide concentration distribution in an indoor space when a person is present in the indoor space. do.
  • An estimation method includes an analysis step, an extraction step, an approximation step, a position information acquisition step, an estimation step, and a notification control step.
  • a simulation is performed to determine the concentration distribution of carbon dioxide in the virtual space as analysis information.
  • the simulation includes three-dimensional information about the indoor space, first environment information about the environment of the indoor space, second environment information about the outdoor environment, and virtual information in which a person's position is set in the virtual space corresponding to the indoor space. Based on location information.
  • information on a virtual concentration distribution of carbon dioxide in the height direction at an arbitrarily determined location in the virtual space is extracted from the analysis information.
  • an analytical formula that approximates the concentration distribution of carbon dioxide in the height direction is obtained based on the information on the virtual concentration distribution extracted in the extraction step.
  • position information acquisition step position information of the person existing in the real space of the indoor space is acquired.
  • estimation step the actual concentration distribution of carbon dioxide in the real space is estimated based on the analytical formula obtained in the approximation step and the position information obtained in the position information acquisition step.
  • notification control step information regarding the actual concentration distribution estimated in the estimation step is notified.
  • a program according to one aspect of the present disclosure is a program for causing one or more processors to execute the estimation method.
  • An estimation system includes an analysis section, an extraction section, an approximation section, a position information acquisition section, an estimation section, and a notification control section.
  • the analysis unit performs a simulation and obtains the concentration distribution of carbon dioxide in the virtual space as analysis information.
  • the simulation includes three-dimensional information about the indoor space, first environment information about the environment of the indoor space, second environment information about the outdoor environment, and virtual information in which a person's position is set in the virtual space corresponding to the indoor space. Based on location information.
  • the extraction unit extracts information on a virtual concentration distribution of carbon dioxide in a height direction at an arbitrarily determined location in the virtual space from the analysis information.
  • the approximation unit calculates an analytical formula that approximates the concentration distribution of carbon dioxide in the height direction based on the information on the virtual concentration distribution extracted by the extraction unit.
  • the position information acquisition unit acquires position information of the person existing in the real space of the indoor space.
  • the estimation section estimates the actual concentration distribution of carbon dioxide in the real space based on the analytical formula obtained by the approximation section and the position information acquired by the position information acquisition section.
  • the notification control unit notifies information regarding the actual concentration distribution estimated by the estimation unit.
  • FIG. 1 is a schematic diagram showing the configuration of an estimation system according to an embodiment.
  • FIG. 2 is a block diagram showing the configuration of the server according to the above.
  • FIG. 3 is a schematic diagram showing a virtual space corresponding to the indoor space according to the above.
  • FIG. 4 is a schematic diagram showing a display screen displayed by the estimation system according to the above.
  • FIG. 5 is a flowchart showing the operation of the server according to the above.
  • FIG. 6 is a flowchart showing the operation of the server according to the first modification.
  • FIG. 7 is a schematic diagram showing an example of approximating the shape of the carbon dioxide concentration distribution around a person according to the second modification.
  • FIG. 8 is a schematic diagram showing an approximate example of the shape of the carbon dioxide concentration distribution around a person according to the third modification.
  • FIG. 9 is a conceptual diagram for explaining a method for calculating the concentration of carbon dioxide around a person according to the above.
  • orthogonal in the present disclosure includes not only a state where the angle between the two is strictly 90 degrees, but also a state where the two are orthogonal within a certain error range. That is, the angle between two orthogonal angles falls within a certain error range (for example, 10 degrees or less) with respect to 90 degrees.
  • parallel as used in the present disclosure includes not only a state in which the two do not strictly intersect, but also a state in which the two are parallel within a certain error range. For example, "parallel” as used in the present disclosure includes an inclination of one to the other of 10 degrees or less.
  • Estimation system 1 is a system that executes an estimation method.
  • the server 4 included in the estimation system 1 executes the estimation method.
  • the concentration distribution of carbon dioxide (CO2) in the indoor space 20 where the person H1 exists in the facility 2 is estimated.
  • the person H1 is, for example, a resident of the facility 2.
  • “Facilities” in this disclosure include residential facilities used for residential purposes, as well as non-residential facilities such as stores (tenants), offices, welfare facilities, educational facilities, hospitals, and factories. Non-residential facilities also include restaurants, playgrounds, hotels, inns, kindergartens, daycare centers, community centers, etc. That is, the facility 2 may be a residential facility such as an apartment or a non-residential facility such as an office building. Furthermore, the facilities 2 also include facilities in which residential facilities and non-residential facilities coexist, such as stores on lower floors and residential units on upper floors. In this embodiment, it is assumed that the facility 2 is a detached house.
  • a plurality of (three in the example of FIG. 1) beacon terminals 25 are provided in the indoor space 20.
  • the plurality of beacon terminals 25 are included in a local positioning system (LPS) for measuring the position of the person H1.
  • LPS local positioning system
  • the local positioning system of this embodiment is a system that estimates the position of a person H1 who has a mobile terminal 3 based on communication information between a plurality of beacon terminals 25 and the mobile terminal 3 carried by the person H1.
  • mobile phone in the present disclosure includes putting the mobile terminal 3 in a bag or pocket, wearing the mobile terminal 3 by hooking it on a strap, etc., holding the mobile terminal 3 in the hand, etc. may include.
  • the mobile terminal 3 is placed in a pocket of clothing worn by a person H1.
  • the local positioning system of this embodiment includes a mobile terminal 3.
  • the mobile terminal 3 transmits position information including communication information with the plurality of beacon terminals 25 to the server 4, for example, via the communication device 29 and the network NT1.
  • the server 4 of this embodiment includes an analysis section 70, an extraction section 71, an approximation section 72, a position information acquisition section 73, an estimation section 74, and a notification control section 75. .
  • the analysis unit 70 performs analysis processing (analysis step).
  • analysis step a simulation is performed to obtain the carbon dioxide concentration distribution in the virtual space 20a (see FIG. 3) corresponding to the indoor space 20 as analysis information.
  • the simulation uses three-dimensional information about the indoor space 20, first environment information about the environment of the indoor space 20, second environment information about the outdoor environment, and a person H1a in the virtual space 20a corresponding to the indoor space 20 (see FIG. 3).
  • the location is based on virtual location information.
  • the extraction unit 71 performs an extraction process (extraction step).
  • extraction step information on the virtual concentration distribution of carbon dioxide in the height direction at an arbitrarily determined location in the virtual space 20a is extracted from the analysis information.
  • the "height direction" in the present disclosure is a direction parallel to the direction of gravity.
  • the approximation unit 72 performs approximation processing (approximation step).
  • approximation step an analytical formula that approximates the concentration distribution of carbon dioxide in the height direction is determined based on the information on the virtual concentration distribution extracted in the extraction step.
  • the location information acquisition unit 73 performs location information acquisition processing (location information acquisition step).
  • position information acquisition step position information of the person H1 existing in the real space of the indoor space 20 is acquired.
  • position information is acquired from the mobile terminal 3 carried by the person H1 present in the indoor space 20.
  • the estimation unit 74 performs estimation processing (estimation step).
  • estimation step the actual concentration distribution of carbon dioxide in the real space of the indoor space 20 is estimated based on the analytical formula obtained in the approximation step and the position information obtained in the position information acquisition step.
  • the notification control unit 75 performs notification control processing (notification control step).
  • notification control step information regarding the actual concentration distribution estimated in the estimation step is notified.
  • the concentration distribution of carbon dioxide in the indoor space 20 where the person H1 exists is affected by the position of the person H1.
  • the estimation system 1 and the estimation method of the present embodiment by using the position information of the person H1 existing in the indoor space 20 and the analytical formula based on simulation, the indoor space when the person H1 exists in the indoor space 20 is The carbon dioxide concentration distribution in the space 20 can be estimated with high accuracy.
  • the estimation system 1 of this embodiment includes a plurality of equipment provided in a facility 2, a mobile terminal 3 carried by a person H1 present in an indoor space 20, a server 4, and an information terminal. 8.
  • the facility 2 includes an indoor space 20, an air supply duct 23, and an exhaust duct 24.
  • the air supply duct 23 is an air passageway that connects the space outside the facility 2 (outside air) and the indoor space 20 of the facility 2, and is a duct for taking outside air into the indoor space 20.
  • the air supply duct 23 has a first air supply port 231 provided on the indoor space 20 side (ceiling surface 21 of the indoor space 20) and a second air supply port 232 provided on the outside of the facility 2. are doing.
  • the second air supply port 232 is provided with, for example, an air supply fan.
  • the exhaust duct 24 is an air passageway that connects the space outside the facility 2 (outside air) and the indoor space 20 of the facility 2, and is a duct for discharging the air in the indoor space 20 to the outside of the facility 2.
  • the exhaust duct 24 has a first exhaust port 241 provided on the indoor space 20 side (ceiling surface 21 of the indoor space 20) and a second exhaust port 242 provided on the outside of the facility 2.
  • the second exhaust port 242 is provided with, for example, an exhaust fan.
  • the facility 2 is provided with a plurality of equipment included in the estimation system 1.
  • the plurality of equipment includes a plurality of beacon terminals 25, a carbon dioxide sensor (CO2 sensor) 27, and a communication device 29.
  • CO2 sensor carbon dioxide sensor
  • the plurality of beacon terminals 25 are provided, for example, on the ceiling surface 21 of the indoor space 20. In the following description, when each of the plurality of beacon terminals 25 is not distinguished, each of the plurality of beacon terminals 25 may be referred to as a "beacon terminal 25.”
  • the beacon terminal 25 is configured to be able to transmit a beacon signal to the mobile terminal 3 carried by the person H1 and which generates position information based on the beacon signal. That is, the beacon terminal 25 of this embodiment functions as a beacon transmitter. Beacon terminal 25 transmits a beacon signal using a predetermined communication method.
  • the predetermined communication method is, for example, a BLE (Bluetooth (registered trademark) Low Energy) communication method. Note that the predetermined communication method is not limited to BLE, and may be a communication method such as WiFi (registered trademark).
  • Beacon terminal 25 transmits beacon signals at predetermined time intervals and with predetermined transmission power.
  • the beacon signal includes unique information (identification information) specific to each of the plurality of beacon terminals 25.
  • the carbon dioxide sensor 27 is a sensor that detects the concentration of carbon dioxide contained in the air.
  • the carbon dioxide sensor 27 is provided near the first exhaust port 241 of the exhaust duct 24. However, the carbon dioxide sensor 27 only needs to be provided at a location where it can detect the concentration of carbon dioxide in the indoor space 20, and may be provided in a location other than the vicinity of the first exhaust port 241.
  • the carbon dioxide sensor 27 is configured to be able to communicate with the communication device 29.
  • the carbon dioxide sensor 27 transmits information on the detected concentration of carbon dioxide to the communication device 29 .
  • Communicatible means that information can be exchanged directly or indirectly via a network, repeater, etc., using an appropriate communication method such as wired communication or wireless communication.
  • the communication device 29 is connected to a network NT1 such as the Internet.
  • the communication device 29 of this embodiment is, for example, a router.
  • the communication device 29 transmits information on the concentration of carbon dioxide detected by the carbon dioxide sensor 27 to the server 4 via the network NT1.
  • the communication device 29 of this embodiment transmits the position information received from the mobile terminal 3 to the server 4 via the network NT1.
  • the communication device 29 transmits to the server 4 information such as the operating status of whether or not the air supply fan and the exhaust fan are operating, or the operating status of air conditioning equipment such as an air conditioner installed in the indoor space 20. You can also send it.
  • the mobile terminal 3 is a mobile terminal such as a smartphone, a tablet terminal, or a notebook personal computer carried by the person H1.
  • a case where the mobile terminal 3 is a smartphone will be exemplified.
  • the mobile terminal 3 is configured to be able to communicate with a plurality of beacon terminals 25 and a communication device 29.
  • the mobile terminal 3 receives a beacon signal transmitted from the beacon terminal 25. That is, the mobile terminal 3 of this embodiment functions as a beacon receiver. Furthermore, the mobile terminal 3 transmits position information including communication information with the plurality of beacon terminals 25 to the server 4 via the communication device 29 and the like.
  • the communication information includes received signal strength indication (RSSI) of the beacon signal and identification information of the beacon terminal 25 included in the beacon signal.
  • RSSI received signal strength indication
  • the mobile terminal 3 may transmit the position information to the server 4 not through the communication device 29 but through a mobile phone network (carrier network) provided by a communication carrier.
  • the location information may include identification information of the mobile terminal 3.
  • the identification information of the mobile terminal 3 is information including, for example, a MAC (Media Access Control) address, an IP (Internet Protocol) address, or a product number.
  • the information terminal 8 is, for example, a desktop or laptop personal computer.
  • the information terminal 8 is a terminal operated by, for example, a resident of the facility 2 (for example, a person H1) or an employee of a company that provides a service to monitor (manage) the facility 2.
  • the information terminal 8 is configured to be able to communicate with the server 4 via the network NT1.
  • the display unit 81 is, for example, a liquid crystal display or an organic EL (Electro Luminescence) display.
  • the display unit 81 displays a display screen D1 (see FIG. 4) under the control of the notification control unit 75 of the server 4.
  • the server 4 is an information terminal installed, for example, in a company that provides a service for monitoring (managing) the environment of the indoor space 20 of the facility 2.
  • the server 4 includes a communication section 5, a storage section 6, and a control section 7.
  • the server 4 includes, for example, a microcomputer having a processor and memory.
  • the computer system functions as the control unit 7 by the processor executing one or more appropriate programs. That is, the control unit 7 is realized by a computer system having a processor and a memory.
  • Each of the one or more programs may be pre-recorded in a memory, or may be provided via a telecommunications line such as the Internet, or recorded on a non-transitory recording medium such as a memory card.
  • the communication unit 5 is configured to be able to communicate with the communication device 29 provided in the facility 2, the mobile terminal 3 carried by the person H1, and the information terminal 8.
  • the storage unit 6 is a semiconductor memory such as ROM (Read Only Memory), RAM (Random Access Memory), or EEPROM (Electrically Erasable Programmable Read Only Memory). Note that the storage unit 6 is not limited to a semiconductor memory, and may be a hard disk drive or the like.
  • the storage unit 6 of this embodiment stores positional information of each of the plurality of beacon terminals 25, three-dimensional information of the indoor space 20, first environmental information regarding the environment of the indoor space 20, and second environmental information regarding the outdoor environment. and virtual position information that sets the position of the person H1a (see FIG. 3) in the virtual space 20a corresponding to the indoor space 20.
  • the three-dimensional information of the indoor space 20 is, for example, three-dimensional model data such as BIM (Building Information Modeling) data of the facility.
  • the first environment information includes air conditioning equipment information including the performance, installation position, driving conditions, etc. of the air conditioning equipment installed in the indoor space 20, and ventilation equipment information including the performance, installation position, driving conditions, etc. of the ventilation equipment.
  • the second environmental information includes information on outside air sound, humidity, season, and concentration of carbon dioxide in the outside air.
  • the virtual position information includes information such as the position of the person H1a existing in the virtual space 20a, state information indicating the state of the person H1a such as a sitting state or a standing state, and the number of persons H1a existing in the virtual space 20a.
  • the first environment information, the second environment information, and the virtual position information are information used in the simulation performed by the analysis unit 70, and are information set in advance before the analysis unit 70 performs the simulation.
  • the control unit 7 includes an analysis unit 70, an extraction unit 71, an approximation unit 72, a position information acquisition unit 73, an estimation unit 74, a notification control unit 75, an actual measurement value acquisition unit 76, and a difference detection unit 78. , and a correction section 79.
  • the analysis unit 70 performs analysis processing (analysis step).
  • the analysis unit 70 performs an analysis step, for example, before the facility 2 is constructed or before a resident lives in the facility 2.
  • a simulation is performed using, for example, the finite element method, and the concentration distribution of carbon dioxide in the virtual space 20a (see FIG. 3) of the indoor space 20 is obtained as analysis information.
  • the simulation set three-dimensional information about the indoor space 20, first environment information about the environment of the indoor space 20, second environment information about the outdoor environment, and the position of the person H1a in the virtual space 20a corresponding to the indoor space 20. Based on virtual location information.
  • the obtained analysis information is stored in the storage unit 6.
  • the simulation is performed multiple times by changing the settings of the first environment information, the settings of the second environment information, or the settings of the virtual position information, and the analysis results (simulation (Results) information is accumulated.
  • the extraction unit 71 performs extraction processing (extraction step).
  • extraction step information on the virtual concentration distribution of carbon dioxide in the height direction at an arbitrarily determined location in the virtual space 20a is extracted from the analysis information.
  • FIG. 3 is a schematic diagram showing a virtual space 20a corresponding to the indoor space 20.
  • the z-axis in FIG. 3 is an axis along the height direction (gravitational direction) of the indoor space 20.
  • the x-axis and the y-axis are axes that are orthogonal to the z-axis and are orthogonal to each other.
  • the virtual space 20a is formed to have the same shape and size as the indoor space 20.
  • the height from the floor surface 22a to the ceiling surface 21a of the virtual space 20a is the same as the height from the floor surface 22 to the ceiling surface 21 of the indoor space 20.
  • a virtual person H1a corresponding to the person H1 is arranged in the virtual space 20a.
  • the analysis information includes information on the concentration distribution of carbon dioxide in the virtual space 20a.
  • the extraction step for example, information on the virtual concentration distribution of carbon dioxide in the height direction of point P1 on the ceiling surface 21a is extracted from the analysis information.
  • the virtual concentration distribution of carbon dioxide in the height direction of point P1 is the virtual concentration distribution of carbon dioxide in the space between point P1 on the ceiling surface 21a and the floor surface 22a directly below point P1.
  • point P1 is an arbitrary location where no person H1a is present in the height direction.
  • information on a virtual concentration distribution of carbon dioxide in the height direction of point P2 on the ceiling surface 21a is extracted from the analysis information.
  • point P2 is an arbitrary location where the person H1a is present in the height direction.
  • arbitrarily determined locations in the virtual space 20a are selected according to the request of the person H1 (user) who is a resident of the facility 2.
  • the request of the person H1 is, for example, a request regarding in what situation (state) the person H1 wants to use (intends to use) the indoor space 20.
  • the requests of the person H1 include, for example, a request for the position (sitting position) of the person H1 in the indoor space 20, a request for the number of people using the indoor space 20, a request for the number of times of ventilation (scheduled), and the like.
  • the approximation unit 72 performs approximation processing (approximation step).
  • approximation step an analytical formula that approximates the concentration distribution of carbon dioxide in the height direction is determined based on the information on the virtual concentration distribution extracted in the extraction step.
  • the virtual temperature distribution of carbon dioxide in the virtual space 20a can be expressed as a function F(x, y, z).
  • x is the value of the x-axis in FIG. 3
  • y is the value of the y-axis in FIG. 3
  • z is the value of the z-axis in FIG.
  • Vxy(z) representing the equation is an analytical expression that approximates the (virtual) concentration distribution of carbon dioxide in the height direction.
  • facility 2 is equipped with ventilation equipment and air conditioning equipment such as air conditioners that supply and ventilate air, and the air in indoor space 20 is constantly being diffused. After accumulation, it is vented outdoors. That is, in the carbon dioxide concentration distribution in the indoor space 20, the concentration of carbon dioxide near the ceiling surface 21 tends to be higher than the concentration of carbon dioxide near the floor surface 22.
  • the concentration distribution of carbon dioxide in the height direction is as follows. It has the characteristic of being easy to express using analytical formulas. If the concentration distribution of carbon dioxide in the height direction is expressed as Vxy (z), the concentration distribution of carbon dioxide in three dimensions (indoor space 20) can be divided into a two-dimensional (xy plane) distribution and a one-dimensional (z direction) distribution. Can be separated. That is, it becomes possible to express the carbon dioxide concentration distribution in three dimensions in a simplified manner.
  • the location information acquisition unit 73 performs location information acquisition processing (location information acquisition step).
  • position information of the person H1 existing in the real space of the indoor space 20 is acquired.
  • location information is acquired (received) from the mobile terminal 3 carried by the person H1.
  • the position information acquisition unit 73 of this embodiment performs position estimation processing (position estimation step).
  • position estimation step the position of the mobile terminal 3 is estimated based on the position information acquired from the mobile terminal 3.
  • the position estimation step the position of the mobile terminal 3 is estimated by performing three-point positioning, for example, based on the received signal strength of the beacon signal, the identification information of the beacon terminal 25 included in the beacon signal, and the position information of the beacon terminal 25. do.
  • the actual measurement value acquisition unit 76 performs actual measurement value acquisition processing (actual measurement value acquisition step).
  • actual measurement value acquisition step an actual measurement value of the concentration of carbon dioxide at a predetermined location in the real space of the indoor space 20 is acquired.
  • the actual measurement value of the concentration of carbon dioxide at the first exhaust port 241 provided in the actual space of the indoor space 20 is acquired.
  • information on the concentration of carbon dioxide detected by the carbon dioxide sensor 27 provided near the first exhaust port 241 is received from the carbon dioxide sensor 27 via the communication device 29 or the like. Get it by doing that.
  • the difference detection unit 78 performs difference detection processing (difference detection step).
  • difference detection step the difference between the measured value of the concentration of carbon dioxide obtained in the measured value acquisition step and the analytical value of the concentration of carbon dioxide at a predetermined location in the virtual space 20a corresponding to a predetermined location in the real space is detected. To detect.
  • the difference detection step of this embodiment the measured value of the concentration of carbon dioxide obtained in the measured value acquisition step and the first exhaust port 241a (exhaust port) of the virtual space 20a corresponding to the first exhaust port 241 (exhaust port) are used. ) to detect the difference between the analysis value and the concentration of carbon dioxide.
  • an analytical value of the concentration of carbon dioxide at the first exhaust port 241a (predetermined location) of the virtual space 20a is obtained from the analytical information stored in the storage unit 6.
  • the correction unit 79 performs correction processing (correction step).
  • the correction step the analytical formula obtained in the approximation step is corrected based on the difference between the measured value and the analytical value of the concentration of carbon dioxide detected in the difference detection step. For example, if the actual measured value has a higher concentration of carbon dioxide than the analytical value, in the correction step, the analytical formula is corrected so that the result of the analytical formula becomes larger by the difference between the actual measured value and the analytical value. Further, for example, when the concentration of carbon dioxide is higher in the analytical value than in the actual measurement value, in the correction step, the analytical expression is corrected so that the result of the analytical expression becomes smaller by the difference between the actual measurement value and the analytical value.
  • the analytical formula can be corrected based on the difference in the concentration of carbon dioxide at a predetermined location in the indoor space 20 (in this embodiment, the first exhaust port 241). , the estimation accuracy of the carbon dioxide concentration distribution can be improved with a small number of sensors (one in this embodiment).
  • the estimation unit 74 performs estimation processing (estimation step).
  • the estimation step the actual concentration distribution of carbon dioxide in the real space of the indoor space 20 is estimated based on the analytical formula obtained in the approximation step and the position information obtained in the position information acquisition step. More specifically, in the estimation step of this embodiment, the actual amount of carbon dioxide in the real space of the indoor space 20 is estimated based on the analytical formula corrected in the correction step and the position information obtained in the position information acquisition step. Estimate the concentration distribution of
  • the notification control unit 75 performs notification control processing (notification control step).
  • notification control step information regarding the actual concentration distribution of carbon dioxide estimated in the estimation step is notified.
  • notification control step of this embodiment information regarding the actual concentration distribution of carbon dioxide estimated in the estimation step is notified by displaying it on the display unit 81 of the information terminal 8.
  • a display screen D1 including an image G1 showing the actual concentration distribution of carbon dioxide estimated in the estimation step is displayed on the display unit 81.
  • the image G1 a space with a higher concentration of carbon dioxide is displayed in a lighter (white) display color, and a space with a lower concentration of carbon dioxide is displayed in a darker (black) display color.
  • the shading of the display color is shown by dot hatching in five levels from level 1 to level 5 (white).
  • the image G11 included in the image G1 is an image showing the person H1
  • the image G12 is an image showing the ceiling surface 21
  • the image G13 is an image showing the floor surface 22.
  • the information regarding the actual concentration distribution of carbon dioxide may be notified by displaying it on a display section other than the display section 81, such as the display section of the mobile terminal 3, for example. good.
  • the estimation system 1 and the estimation method of the present embodiment since the information regarding the actual concentration distribution of carbon dioxide estimated in the estimation step is displayed on the display unit 81 of the information terminal 8, the user can visually see the carbon dioxide concentration distribution. Information about the actual concentration distribution can be grasped.
  • the estimation system 1 performs analysis processing (analysis step) (S1).
  • the estimation system 1 performs a simulation and obtains the carbon dioxide concentration distribution in the virtual space 20a corresponding to the indoor space 20 as analysis information.
  • the estimation system 1 performs an extraction process (extraction step (S2).
  • the estimation system 1 analyzes information on the virtual concentration distribution of carbon dioxide in the height direction at an arbitrarily determined location in the virtual space 20a. Extract from information.
  • the estimation system 1 performs approximation processing (approximation step) (S3).
  • the estimation system 1 obtains an analytical formula that approximates the concentration distribution of carbon dioxide in the height direction based on the information on the virtual concentration distribution extracted in the extraction step.
  • the estimation system 1 performs actual measurement value acquisition processing (actual value acquisition step) (S4).
  • the estimation system 1 acquires the actual measured value of the concentration of carbon dioxide at the first exhaust port 241 (exhaust port) provided in the real space of the indoor space 20 .
  • the estimation system 1 performs a difference detection process (difference detection step) (S5).
  • the estimation system 1 detects the difference between the measured value and the analyzed value of the concentration of carbon dioxide.
  • the estimation system 1 performs a correction process (correction step) (S6).
  • the estimation system 1 of this embodiment corrects the analytical formula based on the difference in carbon dioxide concentration detected in the difference detection step.
  • the estimation system 1 performs a position information acquisition process (position information acquisition step) (S7).
  • the estimation system 1 acquires position information of the person H1 existing in the real space of the indoor space 20.
  • estimation system 1 performs estimation processing (estimation step) (S8).
  • the estimation system 1 of this embodiment estimates the actual concentration distribution in the indoor space 20 based on the analytical formula corrected in the correction step and the position information acquired in the position information acquisition step.
  • the estimation system 1 performs notification control processing (notification control step) (S9).
  • the estimation system 1 of the present embodiment notifies information regarding the actual concentration distribution estimated in the estimation step by displaying it on the display unit 81 of the information terminal 8.
  • the analysis information may be corrected based on the difference in carbon dioxide concentration detected in the difference detection step.
  • the simulation results obtained in the analysis step are corrected based on the difference in carbon dioxide concentration detected in the difference detection step.
  • the estimation system 1 performs analysis processing (analysis step) (S1).
  • the estimation system 1 performs actual measurement value acquisition processing (actual measurement value acquisition step) (S4), and performs difference detection processing (difference detection step) (S5).
  • the measured value of the concentration of carbon dioxide obtained in the measured value acquisition step is 550 ppm
  • a predetermined location (first exhaust port 241) in the virtual space 20a corresponding to a predetermined location (first exhaust port 241) in the real space is assumed to be 550 ppm.
  • the analytical value of the concentration of carbon dioxide at the port 241a is 500 ppm.
  • the difference in carbon dioxide concentration (actual measurement value ⁇ analytical value) determined in the difference detection step is +50 ppm.
  • simulation result after correction is SA (x, y, z)
  • simulation result before correction is SB (x, y, z)
  • correction term difference in carbon dioxide concentration
  • the estimation system 1 performs a correction process (correction step) (S11).
  • equation (1) is used to correct the simulation results obtained in the analysis step.
  • the estimation system 1 performs an extraction process (extraction step) (S12).
  • extraction step the virtual concentration distribution of carbon dioxide in the height direction at an arbitrarily determined location in the virtual space 20a is calculated from the corrected analysis information including the corrected simulation result SA (x, y, z). Extract information (information on concentration distribution after correction).
  • the estimation system 1 performs approximation processing (approximation step) (S13).
  • an analytical formula that approximates the concentration distribution of carbon dioxide in the height direction is determined based on the information on the virtual concentration distribution extracted in the extraction step (information on the concentration distribution after correction).
  • the analytical expression obtained in the approximation step of the first modification is based on the difference in carbon dioxide concentration detected in the difference detection step.
  • a polynomial such as a linear function or a quadratic function can be used.
  • the estimation system 1 performs a position information acquisition process (position information acquisition step) (S7), performs an estimation process (estimation step) (S8), and performs a notification control process (notification control step) (S9).
  • the analysis information is corrected based on the difference in the concentration of carbon dioxide at a predetermined location in the indoor space 20 (the first exhaust port 241 in the first modification). This makes it possible to improve the estimation accuracy of the carbon dioxide concentration distribution with a small number of sensors (one in this embodiment).
  • the concentration distribution of carbon dioxide in the horizontal direction is uniform.
  • the concentration of carbon dioxide is locally higher than that in other spaces due to exhalation, but since carbon dioxide mainly diffuses upward, the carbon dioxide relatively diffuses horizontally. Carbon is low. That is, even if the concentration of carbon dioxide in the horizontal direction in the space A1 around the person H1 is approximated as uniform, the influence on the concentration distribution of carbon dioxide in the entire indoor space 20 is small, so the concentration distribution of carbon dioxide in the horizontal direction is assumed to be uniform.
  • the shape of the carbon dioxide concentration distribution in the space A1 around the person H1 existing in the indoor space 20 may be approximated by a predetermined three-dimensional shape.
  • the concentration distribution of carbon dioxide in the height direction inside a predetermined three-dimensional shape is based on a function Vxy(z) which is an analytical formula.
  • the concentration distribution of carbon dioxide in the horizontal direction inside the predetermined three-dimensional shape is uniform.
  • the predetermined three-dimensional shape is a truncated cone shape centered on the person H1. More specifically, the predetermined three-dimensional shape is a truncated cone shape centered on the person H1, with the top and bottom surfaces perpendicular to the height direction, and the cross-sectional area along the horizontal direction increasing from the bottom to the top. It is a shape that grows larger.
  • the estimation step is performed by approximating the shape of the carbon dioxide concentration distribution in the space A1 around the person H1 existing in the indoor space 20 with a predetermined three-dimensional shape.
  • the time required for the estimation step can be shortened by simplifying the processing.
  • the predetermined three-dimensional shape of the second modification is a truncated cone shape in which the cross-sectional area along the horizontal direction increases as it approaches the top surface from the bottom surface, so the predetermined three-dimensional shape of the second modification is The shape approximates the distribution of carbon dioxide that accumulates in the ceiling surface 21 and spreads horizontally along the ceiling surface 21.
  • the shape of the carbon dioxide concentration distribution in the space A1 around the person H1 existing in the indoor space 20 is It may be approximated by a cylindrical shape orthogonal to . Since the predetermined three-dimensional shape of the third modification is a cylindrical shape in which the top and bottom surfaces are perpendicular to the height direction, the processing in the estimation step can be more simplified and the time required for the estimation step can be further shortened.
  • the predetermined three-dimensional shape of the third modification is a cylindrical shape centered on the person H1 in a plan view from the height direction (top view). Further, the predetermined three-dimensional shape of the third modification has a cross-sectional radius of 0.5 m in the horizontal direction.
  • the exhaled breath of the person H1 rises within a radius of approximately 0.5 m around the person H1 due to the upward airflow around the person H1.
  • the cylindrical shape with a cross-sectional radius of 0.5 m in the horizontal direction has a shape that approximates the distribution of carbon dioxide rising in the updraft around the person H1.
  • FIG. 9 there may be a person H1 (first person) and a person H2 (second person) who is different from the person H1 in the indoor space 20.
  • the shape of the concentration distribution of carbon dioxide in the space A1 around the person H1 and the shape of the concentration distribution of carbon dioxide in the space A2 around the person H2 are assumed to be such that the top and bottom surfaces are perpendicular to the height direction. It is approximated by a cylindrical shape.
  • the concentration of carbon dioxide in space A2 is added to the concentration of carbon dioxide in space A1.
  • the concentration distribution of carbon dioxide in space A3 where space A1 and space A2 overlap is estimated.
  • the space A1 is a space around the person H1 existing in the real space of the indoor space 20.
  • the space A2 is a space around the person H2 existing in the real space of the indoor space 20.
  • the carbon dioxide concentration distribution in space A3 is estimated by adding the carbon dioxide concentration in space A2 to the carbon dioxide concentration in space A1.
  • the concentration distribution of can be estimated with higher accuracy.
  • the predetermined three-dimensional shape may be a columnar shape with a horizontal cross section having a radius of 1.0 m.
  • a columnar shape with a horizontal cross section having a radius of 1.0 m.
  • the exhaled breath of the person H1 reaches a position further away from the person H1, compared to when the person H1 is not making a sound. Therefore, for example, when a sound sensor, a camera, etc. detect that the person H1 is making a sound, the shape of the carbon dioxide concentration distribution around the person H1 can be changed to a circle whose horizontal cross section is 1.0 m in radius. It may be columnar.
  • the cylindrical shape whose cross section in the horizontal direction has a radius of 1.0 m has a shape that approximates the distribution of carbon dioxide when the person H1 is making a sound.
  • the predetermined three-dimensional shape may be, for example, a barrel-shaped columnar shape in which the vicinity of the face of the person H1 is bulged (the cross section in the horizontal direction is large).
  • the exhaled air of the person H1 is mainly released into the indoor space 20 from the nose and mouth of the person H1's face and spreads in the horizontal direction (xy direction), so the barrel-shaped cylindrical shape is a shape that approximates the actual distribution of carbon dioxide. becomes.
  • the barrel-shaped cylindrical bulging portion has a cross-sectional radius of 1.0 m in the horizontal direction, for example.
  • the function equivalent to the estimation system 1 (server 4) or the estimation method according to the above embodiment may be realized by a (computer) program or a non-temporary recording medium on which the program is recorded. It's okay.
  • a program according to one aspect is a program for causing one or more processors to execute the estimation method described in the above embodiment.
  • the estimation system 1 (server 4) or the execution entity of the estimation method in the present disclosure includes a computer system.
  • a computer system mainly consists of a processor and a memory as hardware.
  • a processor executes a program recorded in the memory of a computer system, the function as an execution entity of the estimation system 1 or the estimation method according to the present disclosure is realized.
  • the program may be pre-recorded in the memory of the computer system, may be provided through a telecommunications line, or may be recorded on a non-transitory storage medium readable by the computer system, such as a memory card, optical disc, hard disk drive, etc. may be provided.
  • a processor in a computer system is comprised of one or more electronic circuits including semiconductor integrated circuits (ICs) or large scale integrated circuits (LSIs).
  • the integrated circuits such as IC or LSI referred to herein have different names depending on the degree of integration, and include integrated circuits called system LSI, VLSI (Very Large Scale Integration), or ULSI (Ultra Large Scale Integration).
  • FPGAs Field-Programmable Gate Arrays
  • logic devices that can reconfigure the connections inside the LSI or reconfigure the circuit sections inside the LSI, may also be used as processors. I can do it.
  • the plurality of electronic circuits may be integrated into one chip, or may be provided in a distributed manner over a plurality of chips.
  • a plurality of chips may be integrated into one device, or may be distributed and provided in a plurality of devices.
  • the computer system herein includes a microcontroller having one or more processors and one or more memories. Therefore, the microcontroller is also composed of one or more electronic circuits including semiconductor integrated circuits or large-scale integrated circuits.
  • estimation system 1 it is not an essential configuration for the estimation system 1 that the multiple functions of the estimation system 1 are concentrated in one housing (server 4), and the components of the estimation system 1 are distributed over multiple housings. may be provided. Furthermore, at least some functions of the estimation system 1, for example, some functions of the server 4, may be realized by a cloud (cloud computing) or the like.
  • the mobile terminal 3 may estimate the position of the mobile terminal 3 (itself) based on the beacon signal. That is, the mobile terminal 3 performs, for example, three-point positioning based on the received signal strength of the beacon signal, the identification information of the beacon terminal 25 included in the beacon signal, and the position information of the beacon terminal 25, and determines the position of the mobile terminal 3. may be estimated. Then, the mobile terminal 3 transmits the estimation result to the server 4 as position information.
  • the beacon terminal 25 is a beacon transmitter and the mobile terminal 3 is a beacon receiver, but the beacon terminal 25 may be a beacon receiver and the mobile terminal 3 may be a beacon transmitter. .
  • the beacon terminal 25 transmits position information to the server 4 via the communication device 29 or the like.
  • the actual measurement value acquisition step for example, information on a plurality of carbon dioxides detected by a plurality of carbon dioxide sensors 27 provided at a plurality of different locations in the indoor space 20 may be obtained.
  • a plurality of analytical values of the concentration of carbon dioxide at a plurality of locations in the virtual space 20a that correspond one-to-one with a plurality of locations in the real space are acquired from the analysis information.
  • a plurality of differences between the plurality of actual measured values acquired in the actual measured value acquisition step and the plurality of analysis values corresponding to the plurality of actual measured values are detected.
  • the analytical formula or analysis information (the simulation result obtained in the analysis step) obtained in the approximation step is corrected based on the plurality of differences detected in the difference detection step.
  • the analytical formula or analysis information obtained in the approximation step may be corrected based on the average of a plurality of differences.
  • the indoor space 20 is divided according to a plurality of locations where a plurality of carbon dioxide sensors 27 are provided, and for each divided space, the actual measurement value of the carbon dioxide sensor 27 provided in the divided space is added.
  • the analytical expression or analytical information obtained in the approximation step may be corrected based on the difference between the corresponding analytical value and the corresponding analytical value.
  • the estimation method includes an analysis step, an extraction step, an approximation step, a position information acquisition step, an estimation step, and a notification control step.
  • a simulation is performed to determine the concentration distribution of carbon dioxide in the virtual space (20a) as analysis information.
  • the simulation includes three-dimensional information about the indoor space (20), first environment information about the environment of the indoor space (20), second environment information about the outdoor environment, and a virtual space (20a) corresponding to the indoor space (20). Based on the virtual position information in which the position of the person (H1) is set.
  • the extraction step information on a virtual concentration distribution of carbon dioxide in the height direction at an arbitrarily determined location in the virtual space (20a) is extracted from the analysis information.
  • an analytical formula that approximates the concentration distribution of carbon dioxide in the height direction is determined based on the information on the virtual concentration distribution extracted in the extraction step.
  • position information acquisition step position information of the person (H1) existing in the real space of the indoor space (20) is acquired.
  • estimation step the actual concentration distribution of carbon dioxide in real space is estimated based on the analytical formula obtained in the approximation step and the position information obtained in the position information acquisition step.
  • notification control step information regarding the actual concentration distribution estimated in the estimation step is notified.
  • the indoor space when the person (H1) is present in the indoor space (20) is The carbon dioxide concentration distribution in the space (20) can be estimated with high accuracy.
  • the estimation method further includes an actual measurement value acquisition step, a difference detection step, and a correction step in the first aspect.
  • the actual measurement value acquisition step an actual measurement value of the concentration of carbon dioxide at the exhaust port (first exhaust port 241) provided in the real space is acquired.
  • the difference detection step a difference is detected between the actual measurement value acquired in the actual measurement value acquisition step and the analytical value of the concentration of carbon dioxide at the exhaust port (first exhaust port 241a) in the virtual space (20a) corresponding to the exhaust port. do.
  • the correction step the analytical formula is corrected based on the difference detected in the difference detection step.
  • the estimation step the actual concentration distribution is estimated based on the analytical formula corrected in the correction step and the position information obtained in the position information acquisition step.
  • the estimation method according to the third aspect is the same as the first aspect, further comprising an actual value acquisition step, a difference detection step, and a correction step.
  • the actual measurement value acquisition step an actual measurement value of the concentration of carbon dioxide at the exhaust port (first exhaust port 241) provided in the real space is acquired.
  • the difference detection step a difference is detected between the actual measurement value acquired in the actual measurement value acquisition step and the analytical value of the concentration of carbon dioxide at the exhaust port (first exhaust port 241a) in the virtual space (20a) corresponding to the exhaust port. do.
  • the correction step the analysis information is corrected based on the difference detected in the difference detection step.
  • the user can visually grasp the information regarding the actual concentration distribution of carbon dioxide. I can do it.
  • the shape of the carbon dioxide concentration distribution around the person (H1) existing in the real space is calculated using a predetermined three-dimensional shape. Approximate by shape.
  • the process in the estimation step is simplified and estimated.
  • the time required for each step can be reduced.
  • the estimation step By adding the carbon dioxide concentration in the second space to the carbon dioxide concentration in the first space, the carbon dioxide concentration distribution in the space (A3) where the first space and the second space overlap is estimated.
  • the first space is a space around a first person (person H1) existing in real space.
  • the second space is a space around the second person (person H2) existing in real space.
  • the carbon dioxide concentration in the second space is added to the carbon dioxide concentration in the first space.
  • the carbon dioxide concentration distribution in the indoor space (20) can be estimated with higher accuracy.
  • the predetermined three-dimensional shape is a cylinder whose bottom surface is orthogonal to the height direction.
  • the processing in the estimation step can be more simplified and the time required for the estimation step can be further shortened.
  • the configurations other than the first aspect are not essential to the estimation method and can be omitted as appropriate.
  • the program according to the eighth aspect is a program for causing one or more processors to execute the estimation method according to any one of the first to seventh aspects.
  • the indoor space when the person (H1) is present in the indoor space (20) is The carbon dioxide concentration distribution in the space (20) can be estimated with high accuracy.
  • the estimation system (1) includes an analysis section (70), an extraction section (71), an approximation section (72), a position information acquisition section (73), an estimation section (74), A notification control unit (75) is provided.
  • the analysis unit (70) performs a simulation and obtains the concentration distribution of carbon dioxide in the virtual space (20a) as analysis information.
  • the simulation includes three-dimensional information about the indoor space (20), first environment information about the environment of the indoor space (20), second environment information about the outdoor environment, and a virtual space (20a) corresponding to the indoor space (20). Based on the virtual position information in which the position of the person (H1) is set.
  • the extraction unit (71) extracts information on a virtual concentration distribution of carbon dioxide in the height direction at an arbitrarily determined location in the virtual space (20a) from the analysis information.
  • the approximation unit (72) obtains an analytical formula that approximates the concentration distribution of carbon dioxide in the height direction based on the information on the virtual concentration distribution extracted by the extraction unit (71).
  • the position information acquisition unit (73) acquires position information of a person (H1) existing in the real space of the indoor space (20).
  • the estimation unit (74) estimates the actual concentration distribution of carbon dioxide in real space based on the analytical formula obtained by the approximation unit (72) and the position information acquired by the position information acquisition unit (73). do.
  • the notification control unit (75) notifies information regarding the actual concentration distribution estimated by the estimation unit (74).
  • the indoor space when the person (H1) is present in the indoor space (20) is The carbon dioxide concentration distribution in the space (20) can be estimated with high accuracy.

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Marketing (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本開示の課題は、室内空間に人が存在する場合の室内空間の二酸化炭素の濃度分布を推定することである。推定方法は、解析ステップと、抽出ステップと、近似ステップと、位置情報取得ステップと、推定ステップと、通知制御ステップと、を有する。解析ステップでは、仮想空間における二酸化炭素の濃度分布を解析情報として求める。抽出ステップでは、高さ方向における二酸化炭素の仮想の濃度分布の情報を、解析情報から抽出する。近似ステップでは、高さ方向における二酸化炭素の濃度分布を近似した解析式を求める。位置情報取得ステップでは、人の位置情報を取得する。推定ステップでは、解析式と、位置情報とに基づいて、実空間における二酸化炭素の実際の濃度分布を推定する。通知制御ステップでは、推定ステップにて推定した実際の濃度分布に関する情報を通知する。

Description

推定方法、プログラム、及び推定システム
　本開示は、一般に推定方法、プログラム、及び推定システムに関し、より詳細には、室内空間の二酸化炭素の濃度に関する推定方法、プログラム、及び推定システムに関する。
　特許文献１には、パス構成光学系により空間の測定を行い、重み分布関数を用いて空間の二酸化炭素の濃度を計算する広域分布監視装置が記載されている。
特開２００３－３４４２７７号公報
　人が空間（室内空間）に存在する場合、人の呼気に二酸化炭素が含まれるため、室内空間における二酸化炭素の濃度に影響がある。
　本開示は上記事由に鑑みてなされており、室内空間に人が存在する場合の室内空間の二酸化炭素の濃度分布を推定することができる推定方法、プログラム、及び推定システムを提供することを目的とする。
　本開示の一態様に係る推定方法は、解析ステップと、抽出ステップと、近似ステップと、位置情報取得ステップと、推定ステップと、通知制御ステップと、を有する。前記解析ステップでは、シミュレーションを行い仮想空間における二酸化炭素の濃度分布を解析情報として求める。前記シミュレーションは、室内空間の３次元情報、前記室内空間の環境に関する第１環境情報、室外の環境に関する第２環境情報、及び、前記室内空間に対応する前記仮想空間において人の位置を設定した仮想位置情報に基づく。前記抽出ステップでは、前記仮想空間における任意に定めた箇所の高さ方向における二酸化炭素の仮想の濃度分布の情報を、前記解析情報から抽出する。前記近似ステップでは、前記抽出ステップにて抽出した前記仮想の濃度分布の情報に基づいて、前記高さ方向における二酸化炭素の濃度分布を近似した解析式を求める。前記位置情報取得ステップでは、前記室内空間の実空間に存在する前記人の位置情報を取得する。前記推定ステップでは、前記近似ステップにて求めた前記解析式と、前記位置情報取得ステップにて取得した前記位置情報とに基づいて、前記実空間における二酸化炭素の実際の濃度分布を推定する。前記通知制御ステップでは、前記推定ステップにて推定した前記実際の濃度分布に関する情報を通知する。
　本開示の一態様に係るプログラムは、前記推定方法を、１以上のプロセッサに実行させるためのプログラムである。
　本開示の一態様に係る推定システムは、解析部と、抽出部と、近似部と、位置情報取得部と、推定部と、通知制御部と、を備える。前記解析部は、シミュレーションを行い仮想空間における二酸化炭素の濃度分布を解析情報として求める。前記シミュレーションは、室内空間の３次元情報、前記室内空間の環境に関する第１環境情報、室外の環境に関する第２環境情報、及び、前記室内空間に対応する前記仮想空間において人の位置を設定した仮想位置情報に基づく。前記抽出部は、前記仮想空間における任意に定めた箇所の高さ方向における二酸化炭素の仮想の濃度分布の情報を、前記解析情報から抽出する。前記近似部は、前記抽出部によって抽出された前記仮想の濃度分布の情報に基づいて、高さ方向における二酸化炭素の濃度分布を近似した解析式を求める。前記位置情報取得部は、前記室内空間の実空間に存在する前記人の位置情報を取得する。前記推定部は、前記近似部によって求められた前記解析式と、前記位置情報取得部によって取得された前記位置情報とに基づいて、前記実空間における二酸化炭素の実際の濃度分布を推定する。前記通知制御部は、前記推定部によって推定された前記実際の濃度分布に関する情報を通知する。
図１は、実施形態に係る推定システムの構成を示す概略図である。 図２は、同上に係るサーバの構成を示すブロック図である。 図３は、同上に係る室内空間に対応する仮想空間を示す概略図である。 図４は、同上に係る推定システムが表示する表示画面を示す概略図である。 図５は、同上に係るサーバの動作を示すフローチャートである。 図６は、第１変形例に係るサーバの動作を示すフローチャートである。 図７は、第２変形例に係る人の周囲における二酸化炭素の濃度分布の形状の近似例を示す概略図である。 図８は、第３変形例に係る人の周囲における二酸化炭素の濃度分布の形状の近似例を示す概略図である。 図９は、同上に係る人の周囲における二酸化炭素の濃度の算出方法を説明するための概念図である。
　以下、本開示に関する好ましい実施形態について図面を参照しつつ詳細に説明する。なお、以下に説明する実施形態において互いに共通する要素には同一符号を付しており、共通する要素についての重複する説明は省略する場合がある。以下の実施形態は、本開示の様々な実施形態の一つに過ぎない。実施形態は、本開示の目的を達成できれば、設計等に応じて種々の変更が可能である。本開示において説明する各図は、模式的な図であり、各図中の各構成要素の大きさ及び厚さのそれぞれの比が、必ずしも実際の寸法比を反映しているとは限らない。
　なお、本開示でいう「直交」は、二者間の角度が厳密に９０度である状態だけでなく、二者がある程度の誤差の範囲内で直交する状態も含む意味である。つまり、直交する二者間の角度は、９０度に対してある程度の誤差（一例として１０度以下）の範囲内に収まる。本開示でいう「平行」についても同様に、厳密に二者が交わらない状態だけでなく、二者がある程度の誤差の範囲内で平行である状態も含む意味である。例えば、本開示でいう「平行」は、一方に対する他方の傾きが１０度以下であることを含む。
　（１）概要
　まず、本実施形態に係る推定システム１及び推定方法の概要について、図１～図３を参照して説明する。推定システム１は、推定方法を実行するシステムである。本実施形態では、推定システム１が有するサーバ４が、推定方法を実行する。
　推定方法では、施設２のうち人Ｈ１が存在する室内空間２０における二酸化炭素（ＣＯ２）の濃度分布を推定する。人Ｈ１は、例えば施設２の住人である。
　本開示でいう「施設」は、居住用途で用いられる住宅施設、並びに店舗（テナント）、オフィス、福祉施設、教育施設、病院及び工場等の非住宅施設を含む。非住宅施設は、飲食店、遊技場、ホテル、旅館、幼稚園、保育所及び公民館等も含む。つまり、施設２は、マンション等の住宅施設であってもよいし、オフィスビル等の非住宅施設であってもよい。さらに、施設２は、例えば、低層階が店舗で高層階が住戸というように、住宅施設と非住宅施設とが混在する態様の施設も含む。本実施形態では、施設２が戸建住宅である場合を想定する。
　室内空間２０には、複数（図１の例では３つ）のビーコン端末２５が設けられている。複数のビーコン端末２５は、人Ｈ１の位置を測定するためのローカル測位システム（ＬＰＳ：Local Positioning System）に含まれる。本実施形態のローカル測位システムは、複数のビーコン端末２５と人Ｈ１が携帯する携帯端末３との通信情報に基づいて、携帯端末３を有する人Ｈ１の位置を推定するシステムである。
　なお、本開示でいう「携帯」は、携帯端末３をかばん又はポケット等に入れたり、携帯端末３をストラップ等に引っ掛けたりして身につけること、及び、携帯端末３を手に持つこと等を含み得る。図１の例では、携帯端末３は、人Ｈ１が着ている衣服のポケットに入れられている。
　本実施形態のローカル測位システムは、携帯端末３を含んでいる。携帯端末３は、複数のビーコン端末２５との通信情報を含む位置情報を、例えば通信装置２９及びネットワークＮＴ１を介して、サーバ４に送信する。
　図２に示すように、本実施形態のサーバ４は、解析部７０と、抽出部７１と、近似部７２と、位置情報取得部７３と、推定部７４と、通知制御部７５と、を備える。
　解析部７０は、解析処理（解析ステップ）を行う。解析ステップでは、シミュレーションを行い、室内空間２０に対応する仮想空間２０ａ（図３参照）における二酸化炭素の濃度分布を解析情報として求める。シミュレーションは、室内空間２０の３次元情報、室内空間２０の環境に関する第１環境情報、室外の環境に関する第２環境情報、及び、室内空間２０に対応する仮想空間２０ａにおいて人Ｈ１ａ（図３参照）の位置を設定した仮想位置情報に基づいている。
　抽出部７１は、抽出処理（抽出ステップ）を行う。抽出ステップでは、仮想空間２０ａにおける任意に定めた箇所の高さ方向における二酸化炭素の仮想の濃度分布の情報を、解析情報から抽出する。なお、本開示でいう「高さ方向」は、重力方向と平行な方向である。
　近似部７２は、近似処理（近似ステップ）を行う。近似ステップでは、抽出ステップにて抽出した仮想の濃度分布の情報に基づいて、高さ方向における二酸化炭素の濃度分布を近似した解析式を求める。
　位置情報取得部７３は、位置情報取得処理（位置情報取得ステップ）を行う。位置情報取得ステップでは、室内空間２０の実空間に存在する人Ｈ１の位置情報を取得する。本実施形態では、室内空間２０に存在する人Ｈ１が携帯する携帯端末３から位置情報を取得する。
　推定部７４は、推定処理（推定ステップ）を行う。推定ステップでは、近似ステップにて求めた解析式と、位置情報取得ステップにて取得した位置情報とに基づいて、室内空間２０の実空間における二酸化炭素の実際の濃度分布を推定する。
　通知制御部７５は、通知制御処理（通知制御ステップ）を行う。通知制御ステップでは、推定ステップにて推定した実際の濃度分布に関する情報を通知する。
　人Ｈ１の呼気には二酸化炭素が含まれるため、人Ｈ１が存在する室内空間２０における二酸化炭素の濃度分布は人Ｈ１の位置に影響される。本実施形態の推定システム１及び推定方法によれば、室内空間２０に存在する人Ｈ１の位置情報と、シミュレーションに基づく解析式とを用いることで、室内空間２０に人Ｈ１が存在する場合の室内空間２０の二酸化炭素の濃度分布を精度よく推定することができる。
　（２）詳細
　以下、本実施形態に係る推定システム１の詳細な構成について、図１～図５を参照して説明する。図１に示すように、本実施形態の推定システム１は、施設２に設けられている複数の設備と、室内空間２０に存在する人Ｈ１が携帯する携帯端末３と、サーバ４と、情報端末８と、を備える。
　（２．１）施設の構成
　まず、施設２の詳細について図１を参照して説明する。図１に示すように、施設２は、室内空間２０と、給気ダクト２３及び排気ダクト２４とを備える。
　給気ダクト２３は、施設２の外の空間（外気）と、施設２の室内空間２０とを繋ぐ空気の通路であり、外気を室内空間２０に取り込むためのダクトである。給気ダクト２３は、室内空間２０側（室内空間２０の天井面２１）に設けられている第１給気口２３１と、施設２の外側に設けられている第２給気口２３２とを有している。第２給気口２３２には、例えば給気用のファンが設けられる。
　排気ダクト２４は、施設２の外の空間（外気）と、施設２の室内空間２０とを繋ぐ空気の通路であり、室内空間２０の空気を施設２の外に排出するためのダクトである。排気ダクト２４は、室内空間２０側（室内空間２０の天井面２１）に設けられている第１排気口２４１と、施設２の外側に設けられている第２排気口２４２とを有している。第２排気口２４２には、例えば排気用のファンが設けられる。
　また、施設２には、推定システム１に含まれる複数の設備が設けられている。複数の設備は、複数のビーコン端末２５、二酸化炭素センサ（ＣＯ２センサ）２７と、通信装置２９と、を含む。
　複数のビーコン端末２５は、例えば室内空間２０の天井面２１に設けられている。以下の説明において、複数のビーコン端末２５の各々を区別しない場合、複数のビーコン端末２５の各々のことを「ビーコン端末２５」と呼ぶことがある。
　ビーコン端末２５は、人Ｈ１が携帯しビーコン信号に基づいて位置情報を生成する携帯端末３に、ビーコン信号を送信可能に構成されている。すなわち、本実施形態のビーコン端末２５は、ビーコン送信機として機能する。ビーコン端末２５は、所定の通信方式でビーコン信号を送信する。所定の通信方式は、例えば、ＢＬＥ（Bluetooth（登録商標） Low Energy）の通信方式である。なお、所定の通信方式はＢＬＥに限定されず、所定の通信方式はＷｉＦｉ（登録商標）等の通信方式でもよい。ビーコン端末２５は、所定の時間間隔及び所定の送信電力でビーコン信号を送信する。ビーコン信号には、複数のビーコン端末２５の各々に固有の固有情報（識別情報）を含んでいる。
　二酸化炭素センサ２７は、空気中に含まれる二酸化炭素の濃度の値を検知するセンサである。二酸化炭素センサ２７は、排気ダクト２４の第１排気口２４１の近傍に設けられている。ただし、二酸化炭素センサ２７は、室内空間２０における二酸化炭素の濃度を検知できる箇所に設けられていればよく、第１排気口２４１の近傍以外に設けられていてもよい。二酸化炭素センサ２７は、通信装置２９と通信可能に構成されている。二酸化炭素センサ２７は、検知した二酸化炭素の濃度の情報を、通信装置２９に送信する。
　なお、本開示でいう「通信可能」とは、有線通信又は無線通信の適宜の通信方式により、直接的、又はネットワーク若しくは中継器等を介して間接的に、情報を授受できることを意味する。
　通信装置２９は、インターネット等のネットワークＮＴ１に接続されている。本実施形態の通信装置２９は、例えばルータである。通信装置２９は、二酸化炭素センサ２７が検知した二酸化炭素の濃度の情報を、ネットワークＮＴ１を介してサーバ４に送信する。また、本実施形態の通信装置２９は、携帯端末３から受信した位置情報を、ネットワークＮＴ１を介してサーバ４に送信する。なお、通信装置２９は、給気用ファン及び排気用ファンが駆動しているか否かの稼働状況、又は、室内空間２０に設けられたエアーコンディショナ等の空調機器の稼働状況等をサーバ４に送信してもよい。
　（２．２）携帯端末の構成
　携帯端末３は、人Ｈ１が携帯するスマートフォン、タブレット端末、又はノート型のパーソナルコンピュータ等の携帯端末である。本実施形態では、携帯端末３がスマートフォンである場合を例示する。
　携帯端末３は、複数のビーコン端末２５と、通信装置２９と通信可能に構成されている。携帯端末３は、ビーコン端末２５から送信されるビーコン信号を受信する。すなわち、本実施形態の携帯端末３は、ビーコン受信機として機能する。また、携帯端末３は、複数のビーコン端末２５との通信情報を含む位置情報を、通信装置２９等を介してサーバ４に送信する。通信情報は、ビーコン信号の受信信号強度（ＲＳＳＩ：Received Signal Strength Indication）と、ビーコン信号に含まれるビーコン端末２５の識別情報と、を含む。なお、携帯端末３は、通信装置２９を介さずに、通信事業者が提供する携帯電話網（キャリア網）を介して、サーバ４に位置情報を送信してもよい。
　位置情報には、携帯端末３の識別情報が含まれていてもよい。携帯端末３の識別情報は、例えばＭＡＣ（Media Access Control）アドレス、ＩＰ（Internet Protocol）アドレス、又は、品番等の情報を含む情報である。
　（２．３）情報端末の構成
　情報端末８は、例えばデスクトップ型又はラップトップ型のパーソナルコンピュータ等である。情報端末８は、例えば、施設２の住人（例えば人Ｈ１）や、施設２を監視（管理）するサービスを行う会社等の従業員によって操作される端末である。情報端末８は、ネットワークＮＴ１を介して、サーバ４と通信可能に構成されている。
　表示部８１は、例えば、液晶ディスプレイ、又は有機ＥＬ（Electro Luminescence）ディスプレイ等である。表示部８１は、サーバ４の通知制御部７５による制御に応じて、表示画面Ｄ１（図４参照）を表示する。
　（２．４）サーバの構成
　サーバ４は、例えば施設２の室内空間２０の環境を監視（管理）するサービスを行う会社等に設置される情報端末である。
　図２に示すように、サーバ４は、通信部５と、記憶部６と、制御部７と、を備える。
　サーバ４は、例えば、プロセッサ及びメモリを有するマイクロコンピュータを備える。プロセッサが適宜の１以上のプログラムを実行することにより、コンピュータシステムが制御部７として機能する。つまり、制御部７は、プロセッサ及びメモリを有するコンピュータシステムで実現されている。１以上のプログラムの各々は、メモリに予め記録されていてもよいし、インターネット等の電気通信回線を通じて、又はメモリカード等の非一時的な記録媒体に記録されて提供されてもよい。
　通信部５は、施設２に設けられた通信装置２９、人Ｈ１が携帯する携帯端末３、及び、情報端末８と通信可能に構成されている。
　記憶部６は、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）、又はＥＥＰＲＯＭ（Electrically Erasable Programmable Read Only Memory）等の半導体メモリである。なお、記憶部６は、半導体メモリに限らず、ハードディスクドライブ等であってもよい。本実施形態の記憶部６は、複数のビーコン端末２５の各々の位置情報と、室内空間２０の３次元情報と、室内空間２０の環境に関する第１環境情報と、室外の環境に関する第２環境情報と、室内空間２０に対応する仮想空間２０ａにおいて人Ｈ１ａ（図３参照）の位置を設定した仮想位置情報と、を記憶している。
　室内空間２０の３次元情報は、例えば、施設のＢＩＭ(Building Information Modeling)データ等の３次元モデルデータである。第１環境情報は、室内空間２０に設けられる空調機器の性能、設置位置及び駆動条件等を含む空調機器情報と、換気機器の性能、設置位置及び駆動条件等を含む換気機器情報と、を含む。第２環境情報は、外気音、湿度、季節、及び外気における二酸化炭素の濃度の情報等を含む。仮想位置情報は、仮想空間２０ａに存在する人Ｈ１ａの位置、座位状態又は起立状態等の人Ｈ１ａの状態を示す状態情報と、仮想空間２０ａに存在する人Ｈ１ａの数等の情報を含む。第１環境情報、第２環境情報、及び仮想位置情報は、解析部７０が行うシミュレーションに用いられる情報であり、解析部７０がシミュレーションを行う前に予め設定される情報である。
　制御部７は、解析部７０と、抽出部７１と、近似部７２と、位置情報取得部７３と、推定部７４と、通知制御部７５と、実測値取得部７６と、差異検出部７８と、補正部７９とを有する。
　上述のように、解析部７０は、解析処理（解析ステップ）を行う。解析部７０は、例えば施設２が建設される前、又は、施設２に住人が住む前等に、解析ステップを行う。解析ステップでは、例えば有限要素法等を用いたシミュレーションを行い、室内空間２０の仮想空間２０ａ（図３参照）における二酸化炭素の濃度分布を解析情報として求める。シミュレーションは、室内空間２０の３次元情報、室内空間２０の環境に関する第１環境情報、室外の環境に関する第２環境情報、及び、室内空間２０に対応する仮想空間２０ａにおいて人Ｈ１ａの位置を設定した仮想位置情報に基づいている。解析ステップでは、求めた解析情報を、記憶部６に記憶させる。なお、本実施形態の解析ステップでは、第１環境情報の設定、第２環境情報の設定、又は、仮想位置情報の設定等を変更して、複数回シミュレーションを行い、解析情報に解析結果（シミュレーション結果）の情報を蓄積させる。
　上述のように、抽出部７１は、抽出処理（抽出ステップ）を行う。抽出ステップでは、仮想空間２０ａにおける任意に定めた箇所の高さ方向における二酸化炭素の仮想の濃度分布の情報を、解析情報から抽出する。
　図３は、室内空間２０に対応する仮想空間２０ａを示す概略図である。図３中のｚ軸は、室内空間２０の高さ方向（重力方向）に沿った軸である。ｘ軸及びｙ軸は、ｚ軸と直交する軸であるとともに、互いに直交する軸である。仮想空間２０ａは、室内空間２０と同じ形状、大きさになるように形成されている。例えば、仮想空間２０ａの床面２２ａから天井面２１ａまでの高さは、室内空間２０の床面２２から天井面２１までの高さと同じである。また仮想空間２０ａには、人Ｈ１に対応する仮想の人Ｈ１ａが配置されている。解析情報は、仮想空間２０ａにおける二酸化炭素の濃度分布の情報を含んでいる。
　抽出ステップでは、例えば、天井面２１ａのポイントＰ１の高さ方向における二酸化炭素の仮想の濃度分布の情報を解析情報から抽出する。なお、ポイントＰ１の高さ方向における二酸化炭素の仮想の濃度分布とは、天井面２１ａのポイントＰ１から、ポイントＰ１の直下となる床面２２ａまでの間の空間における二酸化炭素の仮想の濃度分布である。図３に示すように、ポイントＰ１は、高さ方向に人Ｈ１ａが存在しない任意の箇所である。また、抽出ステップでは、例えば、天井面２１ａのポイントＰ２の高さ方向における二酸化炭素の仮想の濃度分布の情報を解析情報から抽出する。図３に示すように、ポイントＰ２は、高さ方向に人Ｈ１ａが存在する任意の箇所である。
　なお、本実施形態の抽出ステップでは、ポイントＰ１及びポイントＰ２等の仮想空間２０ａにおける任意に定めた箇所は、施設２の住人である人Ｈ１（ユーザ）の要望に応じて選定されている。ここで、人Ｈ１の要望とは、例えば、人Ｈ１が室内空間２０をどのような状況（状態）で使用したい（使用するつもり）なのかという要望である。人Ｈ１の要望は、例えば、室内空間２０における人Ｈ１の位置（座る位置）の要望、室内空間２０を利用する人数の要望、又は、換気回数の要望（予定）等を含む。
　近似部７２は、近似処理（近似ステップ）を行う。近似ステップでは、抽出ステップにて抽出した仮想の濃度分布の情報に基づいて、高さ方向における二酸化炭素の濃度分布を近似した解析式を求める。
　より具体的には、仮想空間２０ａにおける二酸化炭素の仮想の温度分布は、関数Ｆ（ｘ，ｙ，ｚ）として表すことができる。なお、ｘは図３中のｘ軸の値であり、ｙは図３中のｙ軸の値であり、ｚは図３中のｚ軸の値である。ここで、関数Ｆ（ｘ，ｙ，ｚ）を、Ｈ（ｘ，ｙ）・Ｖｘｙ（ｚ）と変数分離したときの、位置（ｘ，ｙ）における高さ方向の二酸化炭素の仮想の濃度分布を表す関数Ｖｘｙ（ｚ）が、高さ方向における二酸化炭素の（仮想の）濃度分布を近似した解析式である。
　人Ｈ１の呼気に含まれる二酸化炭素は、人体から発せられた後、人体の体温に起因する人体周辺の上昇気流に乗って天井面２１付近に蓄積する。このため、短中期的には、天井面２１付近の二酸化炭素の濃度は、床面２２付近の二酸化炭素の濃度より高い傾向がある。また、天井面２１付近に蓄積した二酸化炭素は、天井面２１に沿って水平方向に広がっていく。一方、長期的には二酸化炭素は空気より重いため、床面２２付近に蓄積する性質がある。施設２には給気及び換気を行う換気設備及びエアーコンディショナ等の空調機器が設けられており室内空間２０の空気は常に拡散されているため、二酸化炭素は短中期的に天井面２１付近に蓄積した後、換気により屋外に排気されるようになっている。すなわち、室内空間２０における二酸化炭素の濃度分布は、床面２２付近の二酸化炭素の濃度より、天井面２１付近の二酸化炭素の濃度が高くなる傾向にある。
　室内空間２０における二酸化炭素の濃度分布は、床面２２付近の二酸化炭素の濃度より、天井面２１付近の二酸化炭素の濃度が高くなる傾向にあることから、高さ方向における二酸化炭素の濃度分布は解析式で表現し易いという特徴がある。高さ方向の二酸化炭素の濃度分布をＶｘｙ（ｚ）で表すと、３次元（室内空間２０）における二酸化炭素の濃度分布を２次元（ｘｙ平面）の分布と１次元（ｚ方向）の分布に分離できる。すなわち、３次元における二酸化炭素の濃度分布を簡略して表現することが可能となる。
　上述のように、位置情報取得部７３は、位置情報取得処理（位置情報取得ステップ）を行う。位置情報取得ステップでは、室内空間２０の実空間に存在する人Ｈ１の位置情報を取得する。本実施形態の位置情報取得ステップでは、人Ｈ１が携帯する携帯端末３から位置情報を取得（受信）する。
　本実施形態の位置情報取得部７３は、位置推定処理（位置推定ステップ）を行う。位置推定ステップでは、携帯端末３から取得する位置情報に基づいて、携帯端末３の位置を推定する。位置推定ステップでは、ビーコン信号の受信信号強度と、ビーコン信号に含まれるビーコン端末２５の識別情報と、ビーコン端末２５の位置情報とに基づいて、例えば３点測位を行い携帯端末３の位置を推定する。
　実測値取得部７６は、実測値取得処理（実測値取得ステップ）を行う。実測値取得ステップでは、室内空間２０の実空間の所定の箇所における二酸化炭素の濃度の実測値を取得する。本実施形態の実測値取得ステップでは、室内空間２０の実空間に設けられた第１排気口２４１における二酸化炭素の濃度の実測値を取得する。本実施形態の実測値取得ステップでは、第１排気口２４１の近傍に設けられた二酸化炭素センサ２７が検知した二酸化炭素の濃度の情報を、通信装置２９等を介して二酸化炭素センサ２７から受信することで取得する。
　差異検出部７８は、差異検出処理（差異検出ステップ）を行う。差異検出ステップでは、実測値取得ステップにて取得した二酸化炭素の濃度の実測値と、実空間の所定の箇所に対応する仮想空間２０ａの所定の箇所における二酸化炭素の濃度の解析値との差異を検出する。本実施形態の差異検出ステップでは、実測値取得ステップにて取得した二酸化炭素の濃度の実測値と、第１排気口２４１（排気口）に対応する仮想空間２０ａの第１排気口２４１ａ（排気口）における二酸化炭素の濃度の解析値との差異を検出する。
　本実施形態の差異検出ステップでは、記憶部６に記憶された解析情報から、仮想空間２０ａの第１排気口２４１ａ（所定の箇所）における二酸化炭素の濃度の解析値を取得する。
　補正部７９は、補正処理（補正ステップ）を行う。補正ステップでは、差異検出ステップにて検出した二酸化炭素の濃度の実測値と解析値との差異に基づいて、近似ステップにて求めた解析式を補正する。例えば実測値の方が解析値よりも二酸化炭素の濃度が高い場合、補正ステップでは、解析式の結果が実測値と解析値との差異分だけ大きくなるように、解析式を補正する。また、例えば解析値の方が実測値よりも二酸化炭素の濃度が高い場合、補正ステップでは、解析式の結果が実測値と解析値との差異分だけ小さくなるように、解析式を補正する。
　本実施形態の推定システム１及び推定方法によれば、室内空間２０の所定の箇所（本実施形態では第１排気口２４１）での二酸化炭素の濃度の差異に基づいて解析式を補正することで、少ないセンサ数（本実施形態では１つ）で二酸化炭素の濃度分布の推定精度を向上させることができる。
　上述のように、推定部７４は、推定処理（推定ステップ）を行う。推定ステップでは、近似ステップにて求めた解析式と、位置情報取得ステップにて取得した位置情報とに基づいて、室内空間２０の実空間における二酸化炭素の実際の濃度分布を推定する。より具体的には、本実施形態の推定ステップでは、補正ステップにて補正した解析式と、位置情報取得ステップにて取得した位置情報とに基づいて、室内空間２０の実空間における二酸化炭素の実際の濃度分布を推定する。
　上述のように、通知制御部７５は、通知制御処理（通知制御ステップ）を行う。通知制御ステップでは、推定ステップにて推定した二酸化炭素の実際の濃度分布に関する情報を通知する。本実施形態の通知制御ステップでは、推定ステップにて推定した二酸化炭素の実際の濃度分布に関する情報を、情報端末８の表示部８１に表示させることで通知する。
　図４に示すように、例えば本実施形態の通知制御ステップでは、推定ステップにて推定した二酸化炭素の実際の濃度分布を示す画像Ｇ１を含む表示画面Ｄ１を表示部８１に表示させる。画像Ｇ１では、二酸化炭素の濃度が高い空間ほど淡い（白い）表示色で表示され、二酸化炭素の濃度が低い空間ほど濃い（黒い）表示色で表示される。図４では一例として、表示色の濃淡をレベル１～レベル５（白）の５段階でドットハッチングにより示している。なお、画像Ｇ１に含まれる画像Ｇ１１は人Ｈ１を示す画像であり、画像Ｇ１２は天井面２１を示す画像であり、画像Ｇ１３は床面２２を示す画像である。
　なお、通知制御ステップでは、二酸化炭素の実際の濃度分布に関する情報（表示画面Ｄ１）を、例えば携帯端末３の表示部等の表示部８１とは別の表示部に表示させることで通知してもよい。
　本実施形態の推定システム１及び推定方法によれば、推定ステップにて推定した二酸化炭素の実際の濃度分布に関する情報を情報端末８の表示部８１に表示させるため、ユーザは視覚的に二酸化炭素の実際の濃度分布に関する情報を把握することができる。
　（３）推定システムの動作
　次に、図５を参照して推定システム１の動作について説明する。
　まず、推定システム１は、解析処理（解析ステップ）を行う（Ｓ１）。推定システム１は、シミュレーションを行い、室内空間２０に対応する仮想空間２０ａにおける二酸化炭素の濃度分布を解析情報として求める。
　次に、推定システム１は、抽出処理（抽出ステップを行う（Ｓ２）。推定システム１は、仮想空間２０ａにおける任意に定めた箇所の高さ方向における二酸化炭素の仮想の濃度分布の情報を、解析情報から抽出する。
　次に、推定システム１は、近似処理（近似ステップ）を行う（Ｓ３）。推定システム１は、抽出ステップにて抽出した仮想の濃度分布の情報に基づいて、高さ方向における二酸化炭素の濃度分布を近似した解析式を求める。
　次に、推定システム１は、実測値取得処理（実測値取得ステップ）を行う（Ｓ４）。推定システム１は、室内空間２０の実空間に設けられた第１排気口２４１（排気口）における二酸化炭素の濃度の実測値を取得する。
　次に、推定システム１は、差異検出処理（差異検出ステップ）を行う（Ｓ５）。推定システム１は、二酸化炭素の濃度の実測値と解析値との差異を検出する。
　次に、推定システム１は、補正処理（補正ステップ）を行う（Ｓ６）。本実施形態の推定システム１は、差異検出ステップにて検出した二酸化炭素の濃度の差異に基づいて、解析式を補正する。
　次に、推定システム１は、位置情報取得処理（位置情報取得ステップ）を行う（Ｓ７）。推定システム１は、室内空間２０の実空間に存在する人Ｈ１の位置情報を取得する。
　次に、推定システム１は、推定処理（推定ステップ）を行う（Ｓ８）。本実施形態の推定システム１は、補正ステップにて補正した解析式と、位置情報取得ステップにて取得した位置情報とに基づいて、室内空間２０の実際の濃度分布を推定する。
　次に、推定システム１は、通知制御処理（通知制御ステップ）を行う（Ｓ９）。本実施形態の推定システム１は、推定ステップにて推知した実際の濃度分布に関する情報を、情報端末８の表示部８１に表示させることで通知する。
　図５に示すフローチャートは、一例に過ぎず、処理の順番が適宜変更されてもよいし、処理が適宜追加又は削除されてもよい。
　（４）変形例
　以下、上記実施形態の変形例を列挙する。以下に説明する変形例は、適宜組み合わせて適用可能である。
　（４．１）第１変形例
　補正ステップでは、差異検出ステップにて検出した二酸化炭素の濃度の差異に基づいて、解析情報を補正してもよい。言い換えると、補正ステップでは、差異検出ステップにて検出した二酸化炭素の濃度の差異に基づいて、解析ステップで求めたシミュレーション結果を補正する。
　図６を参照して、第１変形例に係る推定システム１の動作を説明する。まず、推定システム１は、解析処理（解析ステップ）を行う（Ｓ１）。
　次に推定システム１は、実測値取得処理（実測値取得ステップ）を行い（Ｓ４）、差異検出処理（差異検出ステップ）を行う（Ｓ５）。
　ここで、実測値取得ステップにて取得した二酸化炭素の濃度の実測値が５５０ｐｐｍであり、実空間の所定の箇所（第１排気口２４１）に対応する仮想空間２０ａの所定の箇所（第１排気口２４１ａ）における二酸化炭素の濃度の解析値が５００ｐｐｍである場合を想定する。差異検出ステップで求める二酸化炭素の濃度の差異（実測値－解析値）は、＋５０ｐｐｍである。
　補正後のシミュレーション結果をＳＡ（ｘ，ｙ，ｚ）とし、補正前のシミュレーション結果をＳＢ（ｘ，ｙ，ｚ）とし、補正項（二酸化炭素の濃度の差異）をγとすると、補正後のシミュレーション結果は、以下の式（１）で表すことができる。
［数１］
ＳＡ（ｘ，ｙ，ｚ）＝ＳＢ（ｘ，ｙ，ｚ）＋γ　　　式（１）
　推定システム１は、補正処理（補正ステップ）を行う（Ｓ１１）。第１変形例の補正ステップでは、式（１）を用いて、解析ステップで求めたシミュレーション結果を補正する。
　次に、推定システム１は、抽出処理（抽出ステップ）を行う（Ｓ１２）。抽出ステップでは、補正後のシミュレーション結果であるＳＡ（ｘ，ｙ，ｚ）を含む補正後の解析情報から、仮想空間２０ａにおける任意に定めた箇所の高さ方向における二酸化炭素の仮想の濃度分布の情報（補正後の濃度分布の情報）を抽出する。
　次に、推定システム１は、近似処理（近似ステップ）を行う（Ｓ１３）。近似ステップでは、抽出ステップにて抽出した仮想の濃度分布の情報（補正後の濃度分布の情報）に基づいて、高さ方向における二酸化炭素の濃度分布を近似した解析式を求める。第１変形例の近似ステップにて求める解析式は、差異検出ステップにて検出した二酸化炭素の濃度の差異に基づいている。解析式としては、例えば１次関数、２次関数などの多項式を用いることができる。
　そして、推定システム１は、位置情報取得処理（位置情報取得ステップ）を行い（Ｓ７）、推定処理（推定ステップ）を行って（Ｓ８）、通知制御処理（通知制御ステップ）を行う（Ｓ９）。
　第１変形例の推定システム１及び推定方法によれば、室内空間２０の所定の箇所（第１変形例では第１排気口２４１）での二酸化炭素の濃度の差異に基づいて解析情報を補正することで、少ないセンサ数（本実施形態では１つ）で二酸化炭素の濃度分布の推定精度を向上させることができる。
　図６に示すフローチャートは、一例に過ぎず、処理の順番が適宜変更されてもよいし、処理が適宜追加又は削除されてもよい。
　（４．２）第２変形例
　推定ステップにて室内空間２０の実空間における二酸化炭素の実際の濃度分布を推定するとき、関数Ｖｘｙ（ｚ）の水平方向（ｘｙ方向）に広がりをもたせてもよい。つまり、推定ステップにて室内空間２０の実空間における二酸化炭素の実際の濃度分布を推定するとき、関数Ｖｘｙ（ｚ）の水平方向を所定範囲で設定してもよい。
　関数Ｖｘｙ（ｚ）の水平方向を所定範囲で設定する場合、水平方向における二酸化炭素の濃度分布は均一であるとする。人Ｈ１の周辺の空間Ａ１は呼気により局所的に二酸化炭素の濃度が他の空間と比べて高濃度であるが、二酸化炭素は上方へ主に拡散するため、相対的に水平方向へ拡散する二酸化炭素は少ない。すなわち、人Ｈ１の周辺の空間Ａ１における水平方向の二酸化炭素の濃度を均一と近似しても、室内空間２０全体の二酸化炭素の濃度分布への影響が小さいため、水平方向における二酸化炭素の濃度分布を均一であるとしている。
　例えば、推定ステップでは、室内空間２０に存在する人Ｈ１の周囲の空間Ａ１における二酸化炭素の濃度分布の形状を所定の立体形状で近似してもよい。所定の立体形状の内部における高さ方向の二酸化炭素の濃度分布は解析式である関数Ｖｘｙ（ｚ）に基づいている。所定の立体形状の内部における水平方向の二酸化炭素の濃度分布は均一である。
　図７に示すように、所定の立体形状は、人Ｈ１を中心とした円錐台状である。より具体的には、所定の立体形状は、人Ｈ１を中心とした円錐台状であり、上面及び底面が高さ方向と直交し、底面から上面に近付くにつれて水平方向に沿った断面の面積が大きくなる形状である。
　第２変形例の推定システム１及び推定方法によれば、室内空間２０に存在する人Ｈ１の周囲の空間Ａ１における二酸化炭素の濃度分布の形状を所定の立体形状で近似することで、推定ステップでの処理を簡易的にして推定ステップにかかる時間を短縮することができる。また、第２変形例の所定の立体形状は、底面から上面に近付くにつれて水平方向に沿った断面の面積が大きくなる円錐台状なので、第２変形例の所定の立体形状は、天井面２１付近に蓄積し天井面２１に沿って水平方向に広がる二酸化炭素の分布を近似した形状となっている。
　（４．３）第３変形例
　図８に示すように、推定ステップでは、室内空間２０に存在する人Ｈ１の周囲の空間Ａ１における二酸化炭素の濃度分布の形状を、上面及び底面が高さ方向と直交する円柱状で近似してもよい。第３変形例の所定の立体形状は、上面及び底面が高さ方向と直交する円柱状なので、推定ステップでの処理をより簡易的にして推定ステップにかかる時間をより短縮することができる。
　なお、第３変形例の所定の立体形状は、高さ方向からの平面視（上面視）において、人Ｈ１を中心とする円柱状である。また、第３変形例の所定の立体形状は、水平方向における断面の半径が０．５ｍである。人Ｈ１の呼気は、人Ｈ１を中心とした概ね半径０．５ｍの範囲内で、人Ｈ１周辺の上昇気流で上昇する。水平方向における断面の半径が０．５ｍの円柱状は、人Ｈ１周辺の上昇気流で上昇する二酸化炭素の分布を近似した形状となる。
　また、図９に示すように、室内空間２０に、人Ｈ１（第１の人）と人Ｈ１とは別の人Ｈ２（第２の人）が存在することがある。図９の例では、人Ｈ１の周囲の空間Ａ１における二酸化炭素の濃度分布の形状、及び、人Ｈ２の周囲の空間Ａ２における二酸化炭素の濃度分布の形状を、上面及び底面が高さ方向と直交する円柱状で近似している。
　推定ステップでは、空間Ａ１（第１空間）の少なくとも一部と空間Ａ２（第２空間）の少なくとも一部とが重なる場合、空間Ａ１の二酸化炭素の濃度に空間Ａ２の二酸化炭素の濃度を加算することで、空間Ａ１と空間Ａ２とが重なる空間Ａ３における二酸化炭素の濃度分布を推定する。空間Ａ１は、室内空間２０の実空間に存在する人Ｈ１の周囲の空間である。空間Ａ２は、室内空間２０の実空間に存在する人Ｈ２の周囲の空間である。
　空間Ａ１と空間Ａ２とが重なる場合に、空間Ａ１の二酸化炭素の濃度に空間Ａ２の二酸化炭素の濃度を加算して空間Ａ３における二酸化炭素の濃度分布を推定することで、室内空間２０の二酸化炭素の濃度分布をより精度よく推定することができる。
　所定の立体形状は、水平方向における断面が半径１．０ｍの円柱状であってもよい。例えば、人Ｈ１が会話をしている等声を発している場合には、人Ｈ１が声を発していない場合と比較して、人Ｈ１の呼気は人Ｈ１からより離れた位置に到達する。そのため、例えば、音センサ及びカメラ等で人Ｈ１が声を発していることを検知した場合に、人Ｈ１の周囲における二酸化炭素の濃度分布の形状を、水平方向における断面が半径１．０ｍの円柱状としてもよい。水平方向における断面が半径１．０ｍの円柱状は、人Ｈ１が声を発している場合の二酸化炭素の分布を近似した形状となる。
　なお、所定の立体形状は、例えば人Ｈ１の顔付近が膨らんだ（水平方向における断面が大きい）樽型の円柱状であってもよい。人Ｈ１の呼気は主に人Ｈ１の顔の鼻及び口から室内空間２０に放出されて水平方向（ｘｙ方向）に広がるため、樽型の円柱状は、実際の二酸化炭素の分布を近似した形状となる。なお、樽型の円柱状の膨らんだ部分は、例えば水平方向における断面の半径が１．０ｍである。
　（４．４）他の変形例
　上記実施形態に係る推定システム１（サーバ４）又は推定方法と同等の機能は、（コンピュータ）プログラム、又はプログラムを記録した非一時的記録媒体等で具現化されてもよい。一態様に係るプログラムは、上記実施形態で説明した推定方法を、１以上のプロセッサに実行させるためのプログラムである。
　本開示における推定システム１（サーバ４）又は推定方法の実行主体は、コンピュータシステムを含んでいる。コンピュータシステムは、ハードウェアとしてのプロセッサ及びメモリを主構成とする。コンピュータシステムのメモリに記録されたプログラムをプロセッサが実行することによって、本開示における推定システム１又は推定方法の実行主体としての機能が実現される。プログラムは、コンピュータシステムのメモリに予め記録されてもよく、電気通信回線を通じて提供されてもよく、コンピュータシステムで読み取り可能なメモリカード、光学ディスク、ハードディスクドライブ等の非一時的記録媒体に記録されて提供されてもよい。コンピュータシステムのプロセッサは、半導体集積回路（ＩＣ）又は大規模集積回路（ＬＳＩ）を含む１又は複数の電子回路で構成される。ここでいうＩＣ又はＬＳＩ等の集積回路は、集積の度合いによって呼び方が異なっており、システムＬＳＩ、ＶＬＳＩ（Very Large Scale Integration）、又はＵＬＳＩ（Ultra Large Scale Integration）と呼ばれる集積回路を含む。さらに、ＬＳＩの製造後にプログラムされる、ＦＰＧＡ（Field-Programmable Gate Array）、又はＬＳＩ内部の接合関係の再構成若しくはＬＳＩ内部の回路区画の再構成が可能な論理デバイスについても、プロセッサとして採用することができる。複数の電子回路は、１つのチップに集約されていてもよいし、複数のチップに分散して設けられていてもよい。複数のチップは、１つの装置に集約されていてもよいし、複数の装置に分散して設けられていてもよい。ここでいうコンピュータシステムは、１以上のプロセッサ及び１以上のメモリを有するマイクロコントローラを含む。したがって、マイクロコントローラについても、半導体集積回路又は大規模集積回路を含む１又は複数の電子回路で構成される。
　また、推定システム１における複数の機能が、１つの筐体内（サーバ４）に集約されていることは推定システム１に必須の構成ではなく、推定システム１の構成要素は、複数の筐体に分散して設けられていてもよい。さらに、推定システム１の少なくとも一部の機能、例えば、サーバ４の一部の機能がクラウド（クラウドコンピューティング）等によって実現されてもよい。
　携帯端末３がビーコン信号に基づいて携帯端末３（自身）の位置を推定してもよい。すなわち、携帯端末３が、ビーコン信号の受信信号強度と、ビーコン信号に含まれるビーコン端末２５の識別情報と、ビーコン端末２５の位置情報とに基づいて、例えば３点測位を行い携帯端末３の位置を推定してもよい。そして、携帯端末３は、推定結果を位置情報として、サーバ４に送信する。
　上記実施形態では、ビーコン端末２５がビーコン送信機であり携帯端末３がビーコン受信機である場合を例示したが、ビーコン端末２５がビーコン受信機であり携帯端末３がビーコン送信機であってもよい。ビーコン端末２５がビーコン受信機である場合、例えばビーコン端末２５は、位置情報を、通信装置２９等を介してサーバ４に送信する。
　実測値取得ステップでは、例えば、室内空間２０の互いに異なる複数の箇所に設けられる複数の二酸化炭素センサ２７が検知した複数の二酸化炭素の情報を取得してもよい。差異検出ステップでは、実空間の複数の箇所と１対１で対応する仮想空間２０ａの複数の箇所における二酸化炭素の濃度の複数の解析値を解析情報から取得する。そして差異検出ステップでは、実測値取得ステップにて取得した複数の実測値と、複数の実測値に対応する複数の解析値との複数の差異を検出する。補正ステップでは、差異検出ステップにて検出した複数の差異に基づいて、近似ステップにて求めた解析式又は解析情報（解析ステップにて求めたシミュレーション結果）を補正する。補正ステップでは、複数の差異の平均に基づいて、近似ステップにて求めた解析式又は解析情報を補正してもよい。また、補正ステップでは、複数の二酸化炭素センサ２７が設けられた複数の箇所に応じて室内空間２０を区分して、区分した空間毎に、区分した空間に設けられた二酸化炭素センサ２７における実測値と対応する解析値との差異に基づいて、近似ステップにて求めた解析式又は解析情報を補正してもよい。
　（まとめ）
　以上説明したように、第１の態様に係る推定方法は、解析ステップと、抽出ステップと、近似ステップと、位置情報取得ステップと、推定ステップと、通知制御ステップと、を有する。解析ステップでは、シミュレーションを行い仮想空間（２０ａ）における二酸化炭素の濃度分布を解析情報として求める。シミュレーションは、室内空間（２０）の３次元情報、室内空間（２０）の環境に関する第１環境情報、室外の環境に関する第２環境情報、及び、室内空間（２０）に対応する仮想空間（２０ａ）において人（Ｈ１）の位置を設定した仮想位置情報に基づく。抽出ステップでは、仮想空間（２０ａ）における任意に定めた箇所の高さ方向における二酸化炭素の仮想の濃度分布の情報を、解析情報から抽出する。近似ステップでは、抽出ステップにて抽出した仮想の濃度分布の情報に基づいて、高さ方向における二酸化炭素の濃度分布を近似した解析式を求める。位置情報取得ステップでは、室内空間（２０）の実空間に存在する人（Ｈ１）の位置情報を取得する。推定ステップでは、近似ステップにて求めた解析式と、位置情報取得ステップにて取得した位置情報とに基づいて、実空間における二酸化炭素の実際の濃度分布を推定する。通知制御ステップでは、推定ステップにて推定した実際の濃度分布に関する情報を通知する。
　この態様によれば、室内空間（２０）に存在する人（Ｈ１）の位置情報と、シミュレーションに基づく解析式とを用いることで、室内空間（２０）に人（Ｈ１）が存在する場合の室内空間（２０）の二酸化炭素の濃度分布を精度よく推定することができる。
　第２の態様に係る推定方法は、第１の態様において、実測値取得ステップと、差異検出ステップと、補正ステップと、を更に有する。実測値取得ステップでは、実空間に設けられた排気口（第１排気口２４１）における二酸化炭素の濃度の実測値を取得する。差異検出ステップでは、実測値取得ステップにて取得した実測値と、排気口に対応する仮想空間（２０ａ）の排気口（第１排気口２４１ａ）における二酸化炭素の濃度の解析値との差異を検出する。補正ステップでは、差異検出ステップにて検出した差異に基づいて、解析式を補正する。推定ステップでは、補正ステップにて補正した解析式と、位置情報取得ステップにて取得した位置情報とに基づいて、実際の濃度分布を推定する。
　この態様によれば、例えば排気口（第１排気口２４１）等の室内空間（２０）の所定の箇所での二酸化炭素の濃度の差異に基づいて解析式を補正することで、少ないセンサ数で二酸化炭素の濃度分布の推定精度を向上させることができる。
　第３の態様に係る推定方法は、第１の態様において、実測値取得ステップと、差異検出ステップと、補正ステップと、を更に有する。実測値取得ステップでは、実空間に設けられた排気口（第１排気口２４１）における二酸化炭素の濃度の実測値を取得する。差異検出ステップでは、実測値取得ステップにて取得した実測値と、排気口に対応する仮想空間（２０ａ）の排気口（第１排気口２４１ａ）における二酸化炭素の濃度の解析値との差異を検出する。補正ステップでは、差異検出ステップにて検出した差異に基づいて、解析情報を補正する。
　この態様によれば、例えば排気口（第１排気口２４１）等の室内空間（２０）の所定の箇所での二酸化炭素の濃度の差異に基づいて解析情報を補正することで、少ないセンサ数で二酸化炭素の濃度分布の推定精度を向上させることができる。
　第４の態様に係る推定方法では、第１から第３のいずれかの態様において、通知制御ステップでは、推定ステップにて推定した実際の濃度分布に関する情報を表示部（８１）に表示させることで通知する。
　この態様によれば、推定ステップにて推定した二酸化炭素の実際の濃度分布に関する情報を表示部（８１）に表示させるため、ユーザは視覚的に二酸化炭素の実際の濃度分布に関する情報を把握することができる。
　第５の態様に係る推定方法では、第１から第４のいずれかの態様において、推定ステップにて、実空間に存在する人（Ｈ１）の周囲における二酸化炭素の濃度分布の形状を所定の立体形状で近似する。
　この態様によれば、室内空間（２０）に存在する人（Ｈ１）の周囲における二酸化炭素の濃度分布の形状を所定の立体形状で近似することで、推定ステップでの処理を簡易的にして推定ステップにかかる時間を短縮することができる。
　第６の態様に係る推定方法では、第５の態様において、第１空間（空間Ａ１）の少なくとも一部と第２空間（空間Ａ２）の少なくとも一部とが重なる場合、推定ステップにて、第１空間の二酸化炭素の濃度に第２空間の二酸化炭素の濃度を加算することで、第１空間と第２空間とが重なる空間（Ａ３）における二酸化炭素の濃度分布を推定する。第１空間は、実空間に存在する第１の人（人Ｈ１）の周囲における空間である。第２空間は、実空間に存在する第２の人（人Ｈ２）の周囲における空間である。
　この態様によれば、第１空間（空間Ａ１）と第２空間（空間Ａ２）とが重なる場合に、第１空間の二酸化炭素の濃度に第２空間の二酸化炭素の濃度を加算して空間（Ａ３）における二酸化炭素の濃度分布を推定することで、室内空間（２０）の二酸化炭素の濃度分布をより精度よく推定することができる。
　第７の態様に係る推定方法では、第５又は第６の態様において、所定の立体形状は、底面が高さ方向と直交する円柱状である。
　この態様によれば、所定の立体形状は、底面が高さ方向と直交する円柱状なので、推定ステップでの処理をより簡易的にして推定ステップにかかる時間をより短縮することができる。
　第１の態様以外の構成については、推定方法に必須の構成ではなく、適宜省略可能である。
　第８の態様に係るプログラムは、第１から第７のいずれかの態様に係る推定方法を、１以上のプロセッサに実行させるためのプログラムである。
　この態様によれば、室内空間（２０）に存在する人（Ｈ１）の位置情報と、シミュレーションに基づく解析式とを用いることで、室内空間（２０）に人（Ｈ１）が存在する場合の室内空間（２０）の二酸化炭素の濃度分布を精度よく推定することができる。
　第９の態様に係る推定システム（１）は、解析部（７０）と、抽出部（７１）と、近似部（７２）と、位置情報取得部（７３）と、推定部（７４）と、通知制御部（７５）と、を備える。解析部（７０）は、シミュレーションを行い仮想空間（２０ａ）における二酸化炭素の濃度分布を解析情報として求める。シミュレーションは、室内空間（２０）の３次元情報、室内空間（２０）の環境に関する第１環境情報、室外の環境に関する第２環境情報、及び、室内空間（２０）に対応する仮想空間（２０ａ）において人（Ｈ１）の位置を設定した仮想位置情報に基づく。抽出部（７１）は、仮想空間（２０ａ）における任意に定めた箇所の高さ方向における二酸化炭素の仮想の濃度分布の情報を、解析情報から抽出する。近似部（７２）は、抽出部（７１）によって抽出された仮想の濃度分布の情報に基づいて、高さ方向における二酸化炭素の濃度分布を近似した解析式を求める。位置情報取得部（７３）は、室内空間（２０）の実空間に存在する人（Ｈ１）の位置情報を取得する。推定部（７４）は、近似部（７２）によって求められた解析式と、位置情報取得部（７３）によって取得された位置情報とに基づいて、実空間における二酸化炭素の実際の濃度分布を推定する。通知制御部（７５）は、推定部（７４）によって推定された実際の濃度分布に関する情報を通知する。
　この態様によれば、室内空間（２０）に存在する人（Ｈ１）の位置情報と、シミュレーションに基づく解析式とを用いることで、室内空間（２０）に人（Ｈ１）が存在する場合の室内空間（２０）の二酸化炭素の濃度分布を精度よく推定することができる。
１　推定システム
２０　室内空間
２０ａ　仮想空間
２４１　第１排気口（排気口）
２４１ａ　第１排気口（排気口）
７０　解析部
７１　抽出部
７２　近似部
７３　位置情報取得部
７４　推定部
７５　通知制御部
８１　表示部
Ａ１　空間（第１空間）
Ａ２　空間（第２空間）
Ａ３　空間（第１空間と第２空間とが重なる空間）
Ｈ１　人（第１の人）
Ｈ２　人（第２の人）

Claims (9)

  1. 　室内空間の３次元情報、前記室内空間の環境に関する第１環境情報、室外の環境に関する第２環境情報、及び、前記室内空間に対応する仮想空間において人の位置を設定した仮想位置情報に基づいて、シミュレーションを行い前記仮想空間における二酸化炭素の濃度分布を解析情報として求める解析ステップと、
    　前記仮想空間における任意に定めた箇所の高さ方向における二酸化炭素の仮想の濃度分布の情報を、前記解析情報から抽出する抽出ステップと、
    　前記抽出ステップにて抽出した前記仮想の濃度分布の情報に基づいて、前記高さ方向における二酸化炭素の濃度分布を近似した解析式を求める近似ステップと、
    　前記室内空間の実空間に存在する人の位置情報を取得する位置情報取得ステップと、
    　前記近似ステップにて求めた前記解析式と、前記位置情報取得ステップにて取得した前記位置情報とに基づいて、前記実空間における二酸化炭素の実際の濃度分布を推定する推定ステップと、
    　前記推定ステップにて推定した前記実際の濃度分布に関する情報を通知する通知制御ステップと、
    を有する、
    　推定方法。
  2. 　前記実空間に設けられた排気口における二酸化炭素の濃度の実測値を取得する実測値取得ステップと、
    　前記実測値取得ステップにて取得した前記実測値と、前記排気口に対応する前記仮想空間の排気口における二酸化炭素の濃度の解析値との差異を検出する差異検出ステップと、
    　前記差異検出ステップにて検出した前記差異に基づいて、前記解析式を補正する補正ステップと、
    を更に有し、
    　前記推定ステップでは、前記補正ステップにて補正した前記解析式と、前記位置情報取得ステップにて取得した前記位置情報とに基づいて、前記実際の濃度分布を推定する、
    　請求項１に記載の推定方法。
  3. 　前記実空間に設けられた排気口における二酸化炭素の濃度の実測値を取得する実測値取得ステップと、
    　前記実測値取得ステップにて取得した前記実測値と、前記排気口に対応する前記仮想空間の排気口における二酸化炭素の濃度の解析値との差異を検出する差異検出ステップと、
    　前記差異検出ステップにて検出した前記差異に基づいて、前記解析情報を補正する補正ステップと、
    を更に有する、
    　請求項１に記載の推定方法。
  4. 　前記通知制御ステップは、前記推定ステップにて推定した前記実際の濃度分布に関する情報を表示部に表示させることで通知する、
    　請求項１から３のいずれか１項に記載の推定方法。
  5. 　前記推定ステップでは、前記実空間に存在する人の周囲における二酸化炭素の濃度分布の形状を所定の立体形状で近似する、
    　請求項１から３のいずれか１項に記載の推定方法。
  6. 　前記推定ステップでは、前記実空間に存在する第１の人の周囲の空間である第１空間の少なくとも一部と、前記実空間に存在する第２の人の周囲の空間である第２空間の少なくとも一部とが重なる場合、前記第１空間の二酸化炭素の濃度に前記第２空間の二酸化炭素の濃度を加算することで、前記第１空間と前記第２空間とが重なる空間における二酸化炭素の濃度分布を推定する、
    　請求項５に記載の推定方法。
  7. 　前記所定の立体形状は、底面が前記高さ方向と直交する円柱状である、
    　請求項５に記載の推定方法。
  8. 　請求項１から３のいずれか１項に記載の推定方法を、１以上のプロセッサに実行させるための、
    　プログラム。
  9. 　室内空間の３次元情報、前記室内空間の環境に関する第１環境情報、室外の環境に関する第２環境情報、及び、前記室内空間に対応する仮想空間において人の位置を設定した仮想位置情報に基づいて、シミュレーションを行い前記仮想空間における二酸化炭素の濃度分布を解析情報として求める解析部と、
    　前記仮想空間における任意に定めた箇所の高さ方向における二酸化炭素の仮想の濃度分布の情報を、前記解析情報から抽出する抽出部と、
    　前記抽出部によって抽出された前記仮想の濃度分布の情報に基づいて、前記高さ方向における二酸化炭素の濃度分布を近似した解析式を求める近似部と、
    　前記室内空間の実空間に存在する人の位置情報を取得する位置情報取得部と、
    　前記近似部によって求められた前記解析式と、前記位置情報取得部によって取得された前記位置情報とに基づいて、前記実空間における二酸化炭素の実際の濃度分布を推定する推定部と、
    　前記推定部によって推定された前記実際の濃度分布に関する情報を通知する通知制御部と、
    を備える、
    　推定システム。
PCT/JP2023/006667 2022-03-24 2023-02-24 推定方法、プログラム、及び推定システム WO2023181778A1 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2022-048820 2022-03-24
JP2022048820 2022-03-24

Publications (1)

Publication Number Publication Date
WO2023181778A1 true WO2023181778A1 (ja) 2023-09-28

Family

ID=88100516

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2023/006667 WO2023181778A1 (ja) 2022-03-24 2023-02-24 推定方法、プログラム、及び推定システム

Country Status (1)

Country Link
WO (1) WO2023181778A1 (ja)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016185630A1 (ja) * 2015-05-18 2016-11-24 三菱電機株式会社 室内環境モデル作成装置
US20210190360A1 (en) * 2019-04-23 2021-06-24 Lg Electronics Inc. Artificial intelligence device
WO2021235139A1 (ja) * 2020-05-18 2021-11-25 パナソニックＩｐマネジメント株式会社 空調システム、建物及びプログラム

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016185630A1 (ja) * 2015-05-18 2016-11-24 三菱電機株式会社 室内環境モデル作成装置
US20210190360A1 (en) * 2019-04-23 2021-06-24 Lg Electronics Inc. Artificial intelligence device
WO2021235139A1 (ja) * 2020-05-18 2021-11-25 パナソニックＩｐマネジメント株式会社 空調システム、建物及びプログラム

Similar Documents

Publication Publication Date Title
US20220154954A1 (en) System and method for controlling indoor air quality
JP6425854B1 (ja) 大気汚染に対する曝露の監視
JP2007107782A (ja) 空調制御システムおよび空調制御装置
KR20160134107A (ko) 공기질 측정 데이터를 이용한 정보 처리 방법 및 이를 위한 장치, 이를 위한 프로그램을 기록한 컴퓨터 판독 가능한 기록매체
US20190050762A1 (en) Air Quality Prediction System and Method Through Survey, and Air Quality Prediction Server
JP2009150590A (ja) 空調システム
JP7489676B2 (ja) 空調システム、建物及びプログラム
WO2020052167A1 (zh) 空调吹风广角确定方法、装置以及空调
JP2008304124A (ja) 二酸化炭素濃度測定用センサーシステム
JP5377577B2 (ja) データ処理装置及び二酸化炭素濃度測定用センサーシステム
CN205352469U (zh) 基于物联网的室内环境远程监测装置
WO2023181778A1 (ja) 推定方法、プログラム、及び推定システム
CN106030217A (zh) 用于操作空调的方法和设备
JP2012002504A (ja) 二酸化炭素濃度測定用センサーシステム及び二酸化炭素センサー
JP2021006946A (ja) 管理システム、空間設備制御システム、空間設備運用システム、管理方法、プログラム
JP2018078428A (ja) 制御情報出力システム、機器制御システム、機器制御方法、及びプログラム
JP2022151708A (ja) 判定装置、判定システム、判定方法および判定プログラム
JP2022182864A (ja) 密対策システム、密対策方法、及びプログラム
Liu Simulation of Local Climate Control in Shared Offices Based on Occupants Locations and Preferences
US20240028281A1 (en) Display system, display method, and program
KR20190103501A (ko) 사물인터넷에 기반하여 환기장치를 제어하는 관제시스템 및 이를 포함하는 실내 공기질 개선 시스템
JP7557733B2 (ja) 環境制御システム、及び、環境制御方法
CN116583698A (zh) 显示系统、显示方法和程序
JPWO2019244306A1 (ja) 監視システム及び監視方法
JP6474221B2 (ja) 環境診断装置及び環境診断方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 23774377

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2401005665

Country of ref document: TH