WO2019026439A1 - 深部体温推定システム、ヒートストレス警報システム、及び、深部体温推定方法 - Google Patents

深部体温推定システム、ヒートストレス警報システム、及び、深部体温推定方法 Download PDF

Info

Publication number
WO2019026439A1
WO2019026439A1 PCT/JP2018/022899 JP2018022899W WO2019026439A1 WO 2019026439 A1 WO2019026439 A1 WO 2019026439A1 JP 2018022899 W JP2018022899 W JP 2018022899W WO 2019026439 A1 WO2019026439 A1 WO 2019026439A1
Authority
WO
WIPO (PCT)
Prior art keywords
temperature
heat
user
body temperature
deep body
Prior art date
Application number
PCT/JP2018/022899
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 帝人株式会社
Priority to JP2019533946A priority Critical patent/JP6826667B2/ja
Priority to EP18841143.3A priority patent/EP3662819B1/en
Priority to US16/635,390 priority patent/US11291373B2/en
Publication of WO2019026439A1 publication Critical patent/WO2019026439A1/ja

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/1118—Determining activity level
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/01—Measuring temperature of body parts ; Diagnostic temperature sensing, e.g. for malignant or inflamed tissue
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/1121—Determining geometric values, e.g. centre of rotation or angular range of movement
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
    • A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
    • A61B5/6802—Sensor mounted on worn items
    • A61B5/6804—Garments; Clothes
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/7405—Details of notification to user or communication with user or patient ; user input means using sound
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/746—Alarms related to a physiological condition, e.g. details of setting alarm thresholds or avoiding false alarms
    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00—Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/02—Alarms for ensuring the safety of persons
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/02—Details of sensors specially adapted for in-vivo measurements
    • A61B2562/0219—Inertial sensors, e.g. accelerometers, gyroscopes, tilt switches
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/02—Details of sensors specially adapted for in-vivo measurements
    • A61B2562/0223—Magnetic field sensors
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/02—Details of sensors specially adapted for in-vivo measurements
    • A61B2562/0247—Pressure sensors
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/02—Details of sensors specially adapted for in-vivo measurements
    • A61B2562/0271—Thermal or temperature sensors

Definitions

  • the present invention relates to a deep body temperature estimation system, a heat stress alarm system, and a deep body temperature estimation method.
  • Patent Document 1 describes a technique for grasping a heat stress state based on the inside temperature of protective equipment detected by a temperature sensor.
  • the temperature sensor detects the inner temperature of the protective equipment on the premise that the heatstroke can be predicted by grasping the inner temperature of the protective equipment.
  • core body temperature can be an index indicating the danger of heat stress felt by the human body, including heat stroke.
  • core body temperature does not directly predict core body temperature, and is not sufficient as a prediction accuracy of heat stroke.
  • the present invention has been made in view of the above problems, and an object thereof is to provide a deep body temperature estimation system, a heat stress alarm system, and a deep body temperature estimation method capable of accurately estimating a deep body temperature of a user. I assume.
  • a deep body temperature estimation system is A temperature sensor that measures the temperature in the clothes, which is the temperature inside the clothes worn by the user; An acceleration sensor for detecting an acceleration applied to the user, The core temperature which is the temperature in the user's body is estimated.
  • It may further include deep body temperature estimation means for estimating the deep body temperature based on the temperature in the clothes measured by the temperature sensor and the acceleration detected by the acceleration sensor.
  • the system further comprises heat generation amount estimation means for estimating a heat generation amount that is a heat amount generated in the body of the user based on the acceleration detected by the acceleration sensor.
  • the deep body temperature estimation means may estimate the deep body temperature based on the temperature in the clothes measured by the temperature sensor and the heat generation amount estimated by the heat generation amount estimation means.
  • the deep body temperature estimating means may estimate the deep body temperature using a mathematical model indicating a correspondence between the temperature in the clothes, the heat generation amount, and the deep body temperature.
  • the mathematical model is A change amount of a first heat amount which is a heat amount of the user's muscle is indicated by a difference between the heat generation amount and a first heat flow amount which is a heat flow amount supplied to the user's blood from the user's muscle.
  • the formula of 1 The difference between the first heat flow and the second heat flow, which is the heat flow supplied from the blood of the user to the inside of the clothes, of the amount of change in the second heat, which is the heat of the user's blood.
  • the first heat flow rate is the amount of change in the first temperature, which is the temperature of the user's muscle, relative to the first heat resistance, which is the heat resistance from the user's muscle to the user's blood.
  • a third equation represented by a ratio of a difference from a second temperature change amount which is the temperature of the blood;
  • the second heat flow rate is defined by the amount of change in the second temperature and the amount of change in the temperature in the clothes with respect to the second heat resistance value that is the heat resistance value from the blood of the user to the inside of the clothes.
  • a fourth equation which is expressed as a percentage of difference A fifth equation indicating the amount of change in the first temperature as a ratio of the first heat amount to the first heat amount that is the heat capacity of the user's muscle;
  • a variational equation of the second temperature is represented by a simultaneous differential equation including a sixth equation indicating a ratio of the second heat amount to a second heat amount which is a heat capacity of blood of the user.
  • the core body temperature estimation means may estimate the second temperature as the core body temperature.
  • the first heat capacity and the second heat capacity are values based on the physical features of the user,
  • the first thermal resistance value and the second thermal resistance value may be predetermined values.
  • the pressure sensor further comprises a pressure sensor for detecting a pressure around the user
  • the heat generation amount estimation means may estimate the heat generation amount based on the acceleration detected by the acceleration sensor and the air pressure detected by the air pressure sensor.
  • An angular velocity sensor that detects an angular velocity applied to the user; And a magnetic sensor for detecting the direction of the magnetic field around the user.
  • the acceleration sensor, the angular velocity sensor, and the magnetic sensor are mutually fixed,
  • the heat generation amount estimation unit may estimate the heat generation amount based on the acceleration detected by the acceleration sensor, the angular velocity detected by the angular velocity sensor, and the direction detected by the magnetic sensor.
  • the information processing apparatus may further include notification means for notifying an abnormality of the user based on the deep body temperature estimated by the deep body temperature estimation means.
  • the clothes may be clothes having a heat shielding effect.
  • the heat shielding effect may satisfy that the HTI 24 measured by the method according to ISO 9151 is 3 seconds or more.
  • the clothing may have a permeability of 1.0 cm 3 / cm 2 / s or less, which is a method according to JIS L 1096 and measured by a method using a Frazier tester.
  • the clothes may be a method in accordance with JIS L 1099, and the moisture permeability value measured by a cup method may be 1000 g / m 2 / h or less.
  • the thermal conductivity of the clothes may be 1 kcal / h / m / ° C. or less.
  • a heat stress alarm system is: Deep body temperature estimation means for estimating a deep body temperature which is a temperature in the body of the user based on an in-garment temperature which is an inner temperature of clothes worn by the user and an acceleration applied to the user; And alarming means for emitting an alarm based on the core body temperature estimated by the core body temperature estimating means.
  • a deep body temperature estimation method is Measure the temperature inside the clothes, which is the temperature inside the clothes worn by the user, Detecting an acceleration applied to the user; Based on the measured temperature in the clothes and the detected acceleration, a deep body temperature, which is a temperature in the user's body, is estimated.
  • the deep body temperature of the user can be accurately estimated.
  • Configuration diagram of deep body temperature estimation system according to an embodiment of the present invention
  • Configuration diagram of a processing apparatus Diagram showing the storage location of each module
  • Explanation of mathematical model A flowchart showing deep body temperature estimation processing performed by the processing device according to the embodiment of the present invention Figure showing basic data of subjects
  • the figure which shows the result data of Example 1 Diagram showing the relationship between measured values and calculated values
  • the figure which shows the result data of Example 2 The figure which shows the result data 1 of Example 3.
  • Diagram showing the relationship between weight and C1 The figure which shows the result data 2 of Example 3.
  • FIG. 1 is a block diagram of a deep body temperature estimation system 1000 according to an embodiment of the present invention.
  • Core body temperature estimation system 1000 is a system for estimating the body temperature of the user's body (hereinafter referred to as “core body temperature”). More specifically, deep body temperature estimation system 1000 measures the amount of heat generated by the operation of the user (hereinafter referred to as “heat generation amount”) and the temperature of the space between the user and the clothes worn by the user (hereinafter referred to as “in-garment temperature It is a system which presumes deep body temperature of a user using a mathematical model using ". And.”.
  • the core temperature estimation system 1000 has a function of emitting an alarm when there is a possibility of heatstroke from the core temperature.
  • the user is a firefighter who carries out a firefighting activity, and is wearing fireproof clothes as clothes.
  • the user may not be a fireman but an operator working in a hot sun.
  • the clothes worn by the user are not limited to fire protection clothes.
  • Thermal insulation which is an index of thermal insulation effect, refers to the time it takes for the fabric used for clothes to be exposed to a specified flame by the method in accordance with ISO 9151 and the temperature rise to reach 24 ° C ( It can be mentioned that the HTI 24) satisfies 3 seconds or more, but is not limited thereto.
  • examples of the clothes having a high heat shielding effect include clothes having a multilayer structure cloth described in JP-A-2014-091307 and JP-A-2011-106069. That is, it is preferable that the garment has a laminated fabric comprising two or more layers of the outer layer and the inner layer.
  • the thickness of each layer of the laminated fabric greatly affects the heat shielding property. For this reason, for example, as described in JP-A-2010-255124, it is preferable that the thickness of the outer layer and the thickness of the inner layer be clothes using a fabric that satisfies the following formula, but is limited thereto I will not.
  • the fabric shape may change from the steady state under flame exposure. For example, an increase in fabric thickness under flame exposure is contemplated.
  • HTI24 is 3 second or more when the heat insulation of the fabric which comprises clothes is measured by the method based on ISO9151. HTI 24 is the time for the temperature rise to reach 24 ° C. when exposed to a defined flame.
  • the limiting oxygen index (LOI) of the fibers constituting the multilayer structure fabric is 21 or more, preferably 24 or more.
  • the limit oxygen index is the oxygen concentration (%) of the atmosphere required to continue combustion, and when it is 21 or more, it means that self-extinguishing does not continue in normal air, and high heat resistance It can be demonstrated.
  • the limiting oxygen index (LOI) is a value measured by JIS L1091 (method E).
  • the fiber whose outermost layer is the limit oxygen index (LOI) 21 or more.
  • the fibers include meta-aramid fiber, para-aramid fiber, polybenzimidazole fiber, polyimide fiber, polyamideimide fiber, polyetherimide fiber, polyarylate fiber, polyparaphenylene benzobisoxazole fiber, novoloid fiber, polychlore Fibers, flame-retardant acrylic fibers, flame-retardant rayon fibers, flame-retardant polyester fibers, flame-retardant cotton fibers, flame-retardant wool fibers and the like can be mentioned.
  • meta-aramid fibers such as polymetaphenylene isophthalamide, para-aramid fibers for the purpose of improving the strength of woven or knitted fabrics, that is, polyparaphenylene terephthalamide, or fibers obtained by copolymerizing the same with a third component
  • polyparaphenylene terephthalamide copolymer copolyparaphenylene 3.4'-oxydiphenylene terephthalamide is exemplified.
  • easily twistable materials such as polyester fibers, polyamide fibers, nylon fibers, acrylic fibers and the like may be used within the range in which the flame retardancy is not impaired.
  • the fibers may be base-laid fibers or post-dyed fibers. Also, if necessary, flame-retardant processing may be imparted after weaving the fabric.
  • long fibers or short fibers may be used. Further, two or more types of the above fibers may be mixed or mixed and used.
  • a meta-aramid fiber and a para-aramid fiber are preferably exemplified in the form of filaments or mixed yarn and used in the form of spun yarn.
  • the spun yarn used may be single ply or double ply.
  • the mixing ratio of the para-aramid fiber is preferably 5% by weight or more based on the total fibers constituting the fabric, but the para-aramid fiber is easily fibrillated, so the mixing ratio is 60% by weight or less It is preferable to keep the
  • the fabric may be used in the form of a woven fabric, a knitted fabric, a non-woven fabric or the like, but a woven fabric is particularly preferred.
  • the woven fabric may be any woven structure such as plain weave, twill weave and satin weave. In the case of woven or knitted fabric, two types of fibers may be used for cross-weaving and inter-knitting.
  • the fabric used in the outermost layer (outer layer) preferably has a basis weight of 140 to 500 g / m 2 , more preferably 160 to 400 g / m 2 , and still more preferably 200 to 400 g / m 2. It is preferred to use. If the basis weight is less than 140 g / m 2 , sufficient heat resistance may not be obtained, while if the basis weight exceeds 500 g / m 2 , the wearing comfort when heat shielding active clothing is obtained May be inhibited.
  • fabric inner layer tensile modulus 80 ⁇ 800cN / dtex, the thermal conductivity of the fabric 6.0W ⁇ m -1 ⁇ k -1 or less, preferably 5.0W ⁇ m -1 ⁇ k
  • the tensile modulus of elasticity of the fiber is preferably 80 to 800 cN / dtex (more preferably 80 to 500 cN / dtex, still more preferably 120 to 500 cN / dtex).
  • the tensile elastic modulus is less than 80 cN / dtex, depending on the movement and posture of the wearer, when it is used as a thermal barrier active clothing, the fibers are stretched at a part of the part, making the fabric thinner and sufficient thermal barrier effect It may not be obtained.
  • the tensile modulus exceeds 800 cN / dtex, a so-called “stiff” wearing feeling may be provided. Although it may be possible to avoid this by using a spun yarn, it is preferable that the tensile modulus is 800 cN / dtex or less in order to exert a sufficient effect.
  • the fabric weight of the fabric is preferably 60 to 500 g / m 2 (more preferably 80 to 400 g / m 2 , still more preferably 100 to 350 g / m 2 ). If the basis weight is less than 60 g / m 2 , transmission of electromagnetic waves may not be sufficiently prevented. On the other hand, when the basis weight is higher than 500 g / m 2 , the tendency to store heat becomes remarkable, which may impair the heat shielding property, and the lightness may also be impaired.
  • the fibers constituting the above-mentioned multilayer structure fabric there is no particular limitation on the fibers constituting the above-mentioned multilayer structure fabric.
  • a metal or carbon kneaded or surface-adhered In order to improve the absorption and reflection of the electromagnetic wave, it is also possible to use a metal or carbon kneaded or surface-adhered.
  • carbon fibers can be used as the above-mentioned fibers, aramid fibers, polybenzimidazole fibers, polyimide fibers, polyamideimide fibers, polyetherimide fibers, polyarylate fibers, polyparaphenylene benzobisoxazole fibers, novoloid fibers, Fibers made of organic polymers such as polychlor fiber, flame-retardant acrylic fiber, flame-retardant rayon fiber, flame-retardant polyester fiber, flame-retardant cotton fiber, flame-retardant wool fiber (hereinafter sometimes referred to as organic polymer fiber) It can be mentioned
  • fine particles of carbon, gold, silver, copper, aluminum, etc. are contained in an organic polymer fiber or attached to the surface of an organic polymer fiber to improve the electromagnetic wave absorptivity and thermal conductivity.
  • carbon or the like may be contained in the organic polymer fiber or applied to the surface as a pigment or paint containing the carbon or the like.
  • the content or adhesion ratio of the fine particles to the total weight of the organic polymer fiber is preferably 0.05 to 60% by weight, more preferably 0.05 to 40% by weight, although it depends on the specific gravity of the fine particles, etc. .
  • the content is preferably 0.05% by weight or more, more preferably 0.05 to 10% by weight, and still more preferably 0.05% to 5% by weight.
  • aluminum fine particles it is preferably 1% by weight or more, more preferably 1 to 20% by weight, and still more preferably 1 to 10% by weight.
  • the number average particle diameter of the fine particles is preferably 10 ⁇ m or less (more preferably 0.01 to 1 ⁇ m). If carbon fiber, metal fiber, etc. satisfy the requirements such as the above-mentioned LOI value and heat conductivity, they can be used as they are without kneading the fine particles.
  • a fabric having a carbon fiber or metal fiber content of preferably 50% by weight or more, more preferably 80% by weight or more, still more preferably 100% by weight can be mentioned as an example. .
  • the multilayer structure fabric it is also possible to dispose, as an intermediate layer, a laminate obtained by laminating and fixing a moisture-permeable waterproof film to a fabric composed of fibers having an LOI value of 25 or more between the outer layer and the inner layer. Thereby, the penetration of water from the outside can be suppressed while maintaining the comfort as the fabric structure, and therefore, it is more suitable as a protective suit for fire brigade personnel performing a fire fighting activity such as discharging water.
  • the basis weight of the intermediate layer to be used is preferably in the range of 50 to 200 g / m 2 . If the basis weight is less than 50 g / m 2 , sufficient heat shielding performance may not be obtained, while if the basis weight exceeds 200 g / m 2 , the wearing comfort when heat shielding active clothing is obtained May be inhibited.
  • This fabric is preferably laminated with a thin film made of polytetrafluoroethylene or the like having moisture permeability and water resistance, whereby moisture permeability and water resistance and chemical resistance are improved, and the perspiration of the wearer is transpirationed. It can promote and reduce the heat stress of the wearer.
  • the basis weight per unit area of the thin film laminated to the above intermediate layer is preferably in the range of 10 to 50 g / m 2 . Even when the thin film is laminated on the intermediate layer fabric, it is preferable that the fabric weight of the intermediate layer subjected to the processing be in the range of 50 to 200 g / m 2 described above.
  • the value measured by JIS L 1096 is 1.0 cm 3 / cm 2 / s or less, and more preferably 0.01 cm 3 / cm 2 / s or less.
  • the value measured by JIS L 1099 is 1000 g / m 2 / h or less, and more preferably 10 g / m 2 / h or less.
  • a low thermal conductivity fabric it is desirable to use a low thermal conductivity fabric. Desirably, it is 1 kcal / h / m / ° C. or less, more preferably 1 cal / h / m / ° C. or less.
  • the deep body temperature estimation system 1000 includes a processing device 100, a temperature sensor 210, an acceleration sensor 220, an air pressure sensor 230, an angular velocity sensor 240, a magnetic sensor 250, a buzzer 300, and a terminal device. And 400.
  • the processing apparatus 100 and various sensors are mutually connected.
  • the processing device 100 and the buzzer 300 are connected to each other, and the processing device 100 and the terminal device 400 are connected to each other via a network (not shown).
  • the processing device 100 is a module that forms the core of the deep body temperature estimation system 1000.
  • the processing device 100 is, for example, a microcomputer or a personal computer.
  • the processing apparatus 100 estimates the deep body temperature of the user from information acquired from various sensors, using a predetermined mathematical model.
  • the processing apparatus 100 determines that the user is highly likely to develop heat stroke, that is, when it is determined that the estimated deep body temperature is abnormal, the processing apparatus 100 issues an alarm by sound, light, or screen display.
  • the processing device 100 may issue an alarm on its own or may issue an alarm to the buzzer 300 or the terminal device 400.
  • the configuration of the processing apparatus 100 will be described below with reference to FIG.
  • the processing device 100 includes a central processing unit (CPU) 11, a read only memory (ROM) 12, a random access memory (RAM) 13, a flash memory 14, a real time clock (RTC) 15, and a touch screen. 16, a first communication interface 17, and a second communication interface 18.
  • CPU central processing unit
  • ROM read only memory
  • RAM random access memory
  • RTC real time clock
  • the CPU 11 controls the overall operation of the processing device 100.
  • the CPU 11 operates according to a program stored in the ROM 12 and uses the RAM 13 as a work area.
  • the ROM 12 stores programs and data for controlling the overall operation of the processing device 100.
  • the RAM 13 functions as a work area of the CPU 11. That is, the CPU 11 temporarily writes a program or data to the RAM 13 and refers to the program or data as appropriate.
  • the function of the core body temperature estimation means and the function of the heat production amount estimation means are realized, for example, by the CPU 11 executing a program stored in the ROM 12 using the RAM 13 as a work area.
  • the flash memory 14 is a non-volatile memory that stores various types of information.
  • the RTC 15 is a device for clocking.
  • the RTC 15 has a built-in battery, for example, and continues clocking even while the processing apparatus 100 is powered off.
  • the RTC 15 includes, for example, an oscillation circuit including a quartz oscillator.
  • the touch screen 16 detects a touch operation made by the user and supplies a signal indicating the result of the detection to the CPU 11.
  • the touch screen 16 also displays an image based on the image signal supplied from the CPU 11.
  • the touch screen 16 functions as a user interface of the processing device 100.
  • the first communication interface 17 is an interface for connecting the processing apparatus 100 to a network (for example, the Internet) not shown.
  • the processing device 100 communicates with the terminal device 400 via a network (not shown).
  • the first communication interface 17 is, for example, an interface for wireless LAN (Local Area Network).
  • the second communication interface 18 is an interface for connecting the processing device 100 to various sensors and the buzzer 300.
  • the second communication interface 18 may be connected to various sensors and the buzzer 300 via a communication cable (not shown), or may be connected to the various sensors and the buzzer 300 by wireless communication.
  • the second communication interface 18 may be an interface for a wireless LAN, or may be an interface compliant with a universal serial bus (USB) or a recommended standard 232 (RS-232).
  • the temperature sensor 210 is a sensor that measures the ambient temperature.
  • the temperature sensor 210 is, for example, disposed inside the clothes to measure the temperature in the clothes.
  • the temperature sensor 210 supplies the processing device 100 with information (hereinafter referred to as “temperature information”) indicating the measured temperature in the clothes.
  • the temperature sensor 210 is, for example, a mercury type thermometer, an alcohol type thermometer, a thermistor type thermometer, an infrared type thermometer, a resistance temperature detector, or a CMOS (Complememyary MOS) temperature sensor.
  • CMOS Compplememyary MOS
  • the acceleration sensor 220 is a sensor that measures the acceleration applied to the acceleration sensor 220.
  • the acceleration sensor 220 is attached to a user or clothes and measures the acceleration applied to the user, that is, the acceleration generated on the acceleration sensor 220 by the user's operation.
  • the acceleration sensor 220 is assumed to be an acceleration sensor that measures acceleration of three axes of the XYZ axes.
  • the acceleration sensor 220 supplies the processing device 100 with information indicating the X-axis acceleration, the Y-axis acceleration, and the Z-axis acceleration (hereinafter referred to as “acceleration information”).
  • the acceleration sensor 220 is, for example, a semiconductor acceleration sensor, a mechanical displacement measurement acceleration sensor, or an optical acceleration sensor.
  • the semiconductor type acceleration sensor is an acceleration sensor using MEMS (Micro Electro Mechanical Systems) technology, and is, for example, a capacitance type acceleration sensor, a piezoresistance type acceleration sensor, or a gas temperature distribution type acceleration sensor.
  • the barometric pressure sensor 230 is a sensor that measures the surrounding barometric pressure.
  • the barometric pressure sensor 230 is, for example, disposed outside the clothes and measures the barometric pressure around the user.
  • the pressure sensor 230 supplies the processing device 100 with information indicating the measured pressure (hereinafter referred to as “pressure information”).
  • the pressure sensor 230 is, for example, an electric pressure sensor using a semiconductor, and is a capacitive pressure sensor or a vibration sensor.
  • the angular velocity sensor 240 is a sensor that measures the angular velocity applied to the angular velocity sensor 240.
  • the angular velocity sensor 240 is attached to, for example, a user or clothes, and measures the angular velocity applied to the user, that is, the angular velocity generated in the angular velocity sensor 240 by the user's operation.
  • the angular velocity sensor 240 is assumed to be an acceleration sensor that measures acceleration of three axes of the XYZ axes.
  • the angular velocity sensor 240 supplies the processing device 100 with information indicating the measured angular velocity (hereinafter referred to as “angular velocity information”).
  • the angular velocity sensor 240 is, for example, a gyro sensor that uses a quartz oscillator or a piezoelectric ceramic element to extract strain generated by Coriolis force as an electrical signal.
  • the magnetic sensor 250 is a sensor that detects the magnitude and direction of the surrounding magnetic field.
  • the magnetic sensor 250 is attached to, for example, a user or clothes, detects geomagnetism (a magnetic field generated by the earth), and detects a direction.
  • the magnetic sensor 250 supplies information (hereinafter referred to as “direction information”) indicating the detected direction to the processing device 100.
  • the magnetic sensor 250 includes, for example, a Hall element, a magnetoresistive element, and a magnetoresistive element.
  • the buzzer 300 issues an alarm according to the control by the processing device 100.
  • This alarm is generated, for example, when the estimated core temperature is higher than a predetermined threshold, or when the estimated core temperature rises above a predetermined value (for example, 1.5 ° C. or more) This is a sound that informs the user that the user may have heatstroke.
  • the buzzer 300 includes an electromagnet and a diaphragm, and emits a sound in accordance with the electrical signal supplied from the processing device 100.
  • the function of the alarm means is realized, for example, by the processing device 100 and the buzzer 300 cooperating with each other.
  • the terminal device 400 is a device that functions as a user interface of the processing device 100.
  • the terminal device 400 issues an alarm by sound or screen display, for example, according to control by the processing device 100.
  • the terminal device 400 includes, for example, a CPU, a ROM, a RAM, a flash memory, and a communication interface (not shown).
  • the terminal device 400 is, for example, a smartphone, a tablet terminal, a smart watch, or a head mounted display.
  • FIG. 3 is a view showing how the processing apparatus 100, various sensors, and the buzzer 300 are worn on the clothes 500.
  • the garment 500 comprises a backing 510 and an outer surface 520, with a space 530 between the backing 510 and the outer surface 520.
  • Garment 500 has space 530, thereby providing space 540 which is the space inside garment 500 (the space between garment 500 and underwear 600 worn by the user), and space 550 which is the space outside garment 500.
  • the outer fabric 520 is made of a highly heat insulating material.
  • the lining 510 is provided with an inner pocket 511, and the outer surface 520 is provided with an outer pocket 521.
  • the processing device 100, various sensors, and the buzzer 300 are housed, for example, in any one of the inner pocket 511 and the outer pocket 521.
  • the temperature sensor 210 is housed in the inner pocket 511 in order to measure the temperature in the clothes.
  • the pressure sensor 230 is housed in the outer pocket 521 because it is preferable to measure the pressure inside the clothes 500 rather than the pressure inside the clothes 500.
  • the buzzer 300 is accommodated in the outer pocket 521 so that an alarm can be easily transmitted to the user.
  • the processing device 100, the acceleration sensor 220, the angular velocity sensor 240, and the magnetic sensor 250 may be housed in either the inner pocket 511 or the outer pocket 521. However, the axis of the acceleration detected by the acceleration sensor 220 and the axis of the angular velocity detected by the angular velocity sensor 240 are determined based on the direction detected by the magnetic sensor 250. Therefore, the acceleration sensor 220, the angular velocity sensor 240, and the magnetic sensor 250 are fixed to each other so that their positional relationships do not change. In the present embodiment, the processing device 100, the temperature sensor 210, the acceleration sensor 220, the angular velocity sensor 240, and the magnetic sensor 250 are fixed to each other so that their positional relationships do not change.
  • a mathematical model used to estimate core body temperature will be described.
  • a mathematical model is expressed by a simultaneous differential equation composed of six equations from Equation (1) to Equation (6) shown below.
  • the first equation corresponds to, for example, equation (1).
  • the second equation corresponds to, for example, equation (2).
  • the third equation corresponds to, for example, equation (3).
  • the fourth equation corresponds to, for example, equation (4).
  • the fifth equation corresponds to, for example, equation (5).
  • the sixth equation corresponds to, for example, equation (6).
  • Formula (1) is a formula which represents the variation
  • the first amount of heat is the amount of heat that the user's muscle has.
  • the first heat flow rate is a heat flow rate supplied to the user's blood from the user's muscle.
  • Q1 is a first heat quantity
  • dQ1 / dt is a variation of the first heat quantity per unit time.
  • the unit time is, for example, 1 minute, 10 seconds, 1 second. Hereinafter, unit time shall mean one minute.
  • F0 is a heat production amount. The amount of heat production will be described later.
  • F1 is a first heat flow rate.
  • the equation (1) indicates that the amount of change in the amount of heat possessed by a muscle is defined by the difference between the amount of heat production that is the heat flow entering the muscle and the amount of heat flow coming from the muscle into the blood.
  • Formula (2) is a formula which represents the variation
  • the second amount of heat is the amount of heat that the user's blood has.
  • the second heat flow rate is the heat flow rate supplied from the user's blood to the inside of the garment 500.
  • Q2 is a second heat quantity
  • dQ2 / dt is a variation of the second heat quantity per unit time.
  • F2 is a second heat flow rate. Equation (2) shows that the amount of change in the amount of heat the blood has is defined as the difference between the heat flow from the muscle into the blood and the heat flow from the blood to the inside of the garment 500. .
  • Expression (3) is an expression representing the first heat flow rate as a ratio of the difference between the first temperature change amount and the second temperature change amount with respect to the first heat resistance value.
  • the first thermal resistance value is a thermal resistance value from the user's muscle to the user's blood, and is a value indicating the difficulty in transferring heat from the muscle to the blood.
  • the first temperature is the temperature of the user's muscle.
  • the second temperature is the temperature of the user's blood.
  • R1 is a first thermal resistance value.
  • T1 is the amount of change of the first temperature per unit time.
  • T2 is the amount of change of the second temperature per unit time.
  • the heat flow from muscle to blood is defined as the ratio of the difference between the change in muscle temperature and the change in blood temperature to the difficulty of heat transfer from muscle to blood It is shown that.
  • the equation (4) is an equation representing the second heat flow rate as a ratio of the difference between the change amount of the second temperature and the change amount of the temperature in the clothes with respect to the second heat resistance value.
  • the second heat resistance value is a heat resistance value from the user's blood to the inside of the clothes 500, and is a value indicating the difficulty of heat transfer from the blood to the inside of the clothes.
  • R2 is a second thermal resistance value.
  • T3 is a change amount of the temperature in the clothes per unit time.
  • the heat flow from the muscle to the blood is defined as the ratio of the difference between the change in muscle temperature and the change in blood temperature to the difficulty of heat transfer from muscle to blood It is shown that.
  • Formula (5) is a formula which shows the variation
  • the first heat capacity is the heat capacity of the user's muscle.
  • C1 is a first heat capacity. Equation (5) indicates that muscle temperature change is defined by the ratio of muscle heat capacity to the amount of heat the muscle has.
  • Formula (6) is a formula which shows the variation
  • the second heat capacity is the heat capacity of the user's blood.
  • C2 is a second heat capacity. Equation (6) shows that the temperature change of blood is defined as the ratio between the heat capacity of blood and the amount of heat the blood has.
  • FIG. 4 shows that when heat of the heat generation amount (F0) is generated in the muscle, the heat of the first heat flow (F1) flows from the muscle to the blood at a rate according to the first heat resistance value (R1)
  • the heat of the second heat flow (F2) flows from the blood toward the inside of the clothes 500 at a speed according to the second heat resistance value (R2).
  • the second temperature is estimated based on the heat generation amount and the temperature in the clothes, and the second temperature is estimated as the core temperature. That is, in the present embodiment, the second temperature, which is the temperature of the user's blood, is regarded as the user's core temperature.
  • F0 is obtained from acceleration information and pressure information
  • T3 is obtained from temperature information
  • R1, R2, C1, and C2 are fixed values
  • Q1, Q2, F1, F2, T1 Let T2 be a variable.
  • the fixed value is, for example, a value when the correlation coefficient between the calculated value waveform indicating the temporal change of the calculated value of the core temperature and the measured value waveform indicating the temporal change of the measured value of the core temperature is closest to 1. It is determined.
  • the fixed value may be calculated by a correlation equation between a value when the correlation coefficient is closest to 1 and a value based on the physical feature of the user, for example, a value based on the weight of the user. It is.
  • the heat production can basically be determined from the acceleration of the user during exercise. However, it can be expected that the heat production amount can be accurately obtained by considering not only the acceleration of the user but also the amount of movement in the height direction of the user. Therefore, in the present embodiment, a method for obtaining the heat generation amount from the acceleration of the user and the movement amount in the height direction of the user will be described. This method is described in detail in Japanese Patent Laid-Open No. 2008-220517, and therefore will be summarized and described here.
  • the heat-production amount can be obtained, for example, from the four equations (7) to (10) shown below.
  • VM indicates the magnitude of the acceleration of the user. Strictly speaking, VM indicates the magnitude of accumulated acceleration, which is the accumulated value of the acceleration imposed on the user over a unit time.
  • the unit of VM is mG.
  • x is an instantaneous value of acceleration in the x-axis direction.
  • y is an instantaneous value of acceleration in the y-axis direction.
  • z is an instantaneous value of acceleration in the z-axis direction.
  • the VM is calculated based on the acceleration information supplied from the acceleration sensor 220 to the processing device 100. Specifically, the VM is obtained by accumulating, over a unit time, instantaneous values of acceleration indicated by acceleration information acquired every predetermined sampling period (for example, 0.1 seconds).
  • a sampling period shall mean 0.1 second.
  • Gravity acceleration is considered not to affect oxygen consumption and heat production. Accordingly, it is assumed that the above-mentioned instantaneous value of acceleration is a value excluding gravity acceleration.
  • the z-axis direction is the vertical direction
  • a value obtained by excluding the gravitational acceleration from the instantaneous value of the acceleration in the z-axis direction is adopted as z in the equation (7).
  • VO 2 is per unit time
  • unit body weight e.g., 1 kg
  • the unit of VO 2 is mL / kg / min.
  • Hu is a cumulative increase per unit time.
  • Hd is the cumulative descent amount per unit time.
  • the unit of Hu and Hd is m / min.
  • Hu and Hd can be calculated from the amount of change in barometric pressure indicated by barometric pressure information acquired for each predetermined sampling cycle.
  • the change in air pressure is used. That is, when the user moves upward, the air pressure decreases, and when the user moves downward, the air pressure increases.
  • EE represents the total energy consumption.
  • the unit of EE is kcal / min. 4.825 is a predetermined constant, and the unit is kcal / L.
  • BW is weight.
  • the unit of BW is kg. 1000 is a predetermined constant, and the unit is mL / L.
  • F0 heat-production amount per unit time.
  • the unit of F0 is kcal / min.
  • 0.8 is a predetermined constant. 0.8 indicates that of the total energy consumption, 20% is consumed by exercise and 80% is converted to heat.
  • the heat generation amount of the user can be calculated from the acceleration applied to the user and the movement amount of the user in the vertical direction. That is, the heat production amount of the user can be obtained from the acceleration information and the pressure information.
  • the heat generation amount of the user can be calculated from the acceleration applied to the user, and the barometric pressure information is unnecessary.
  • the motion of the user can be accurately identified, and the calculation accuracy of the heat generation amount of the user can be improved.
  • the direction information supplied by the magnetic sensor 250 the movement of the user can be accurately identified, and the calculation accuracy of the heat generation amount of the user can be improved.
  • the core body temperature estimation process is performed, for example, in response to the processing apparatus 100 being powered on.
  • the CPU 11 measures an initial value of the temperature in the clothes (step S101). For example, the CPU 11 acquires temperature information from the temperature sensor 210, and stores the acquired temperature information in the flash memory 14 as information indicating the initial value of the temperature in the clothes.
  • step S102 the CPU 11 obtains acceleration information from the acceleration sensor 220 and stores the acceleration information in the flash memory 14.
  • step S102 will measure air pressure, if the process of step S102 is completed (step S103). For example, the CPU 11 acquires pressure information from the pressure sensor 230 and stores the pressure information in the flash memory 14.
  • step S104 determines whether it is an estimation timing of deep body temperature.
  • the estimation timing of the core temperature is a timing that arrives every unit time, for example, a timing that arrives every minute.
  • the CPU 11 determines, based on the information supplied from the RTC 15, whether or not it is an estimation timing of the deep body temperature. If the CPU 11 determines that it is not the estimation timing of the deep body temperature (step S104: NO), the process returns to step S102. Note that the process of step S102 and the process of step S103 are performed every sampling cycle, for example, every 0.1 seconds.
  • Step S104 If CPU11 judges that it is the presumed timing of deep body temperature (Step S104: YES), it will measure the temperature in clothes (Step S105). For example, the CPU 11 acquires temperature information from the temperature sensor 210, and stores the acquired temperature information in the flash memory 14 as information indicating the temperature in the clothes.
  • step S106 the CPU 11 calculates an accumulated acceleration (step S106). For example, the CPU 11 calculates the cumulative acceleration from the acceleration measured in step S102 using the equation (7).
  • the CPU 11 calculates the movement amount in the height direction (step S107). For example, the CPU 11 calculates the amount of movement in the height direction from the atmospheric pressure measured in step S103.
  • the amount of movement in the height direction is a cumulative amount of increase and a cumulative amount of decrease.
  • step S108 the CPU 11 calculates the heat generation amount using the cumulative acceleration calculated in step S106, the equation (8), the equation (9), and the equation (10).
  • the CPU 11 estimates core body temperature (step S109). For example, the CPU 11 estimates the core body temperature using the temperature in the clothes measured in step S105, the heat generation amount calculated in step S108, and the simultaneous differential equations of equations (1) to (6). Although it is the amount of change in core body temperature that is calculated from simultaneous differential equations, core body temperature is calculated by adding the amount of core body temperature change to a predetermined reference value (for example, 37 degrees) for core body temperature. It can be calculated. In addition, estimating core body temperature is basically synonymous with estimating the amount of change in core body temperature.
  • step S110 determines whether or not the deep body temperature is abnormal. For example, if the core temperature exceeds the threshold (for example, 39 degrees), or if the core temperature rises above a predetermined value (for example, 1.5 degrees), the CPU 11 determines the core temperature Is determined to be abnormal. When the CPU 11 determines that the core temperature is not abnormal (step S110: NO), the process returns to step S102.
  • the threshold for example, 39 degrees
  • a predetermined value for example, 1.5 degrees
  • step S110 when determining that the core temperature is abnormal (step S110: YES), the CPU 11 issues an alarm (step S111). For example, the CPU 11 causes the buzzer 300 to output a sound corresponding to the alarm. When the process of step S111 is completed, the CPU 11 ends the deep body temperature estimation process.
  • simultaneous differential equations indicating the correspondence between the temperature in the clothes, which is the temperature in the clothes of the user, and the heat generation amount, which is the heat generated in the user's body, are used.
  • the core temperature is estimated from Therefore, according to the present embodiment, the deep body temperature of the user can be accurately estimated.
  • core body temperature can be estimated accurately. Therefore, the estimated core temperature can also be used as an indicator of the risk of heat stress felt by the human body, including heat stroke.
  • the temperature of the user's blood is estimated to be regarded as the user's core temperature.
  • the temperature of the blood is considered to be closer to core body temperature than, for example, the temperature of the rectum. Therefore, according to the present embodiment, the deep body temperature of the user can be accurately estimated.
  • the temperature to be measured is the temperature inside the clothes. Therefore, for example, measurement of the temperature is easier than when measuring the temperature of the rectum.
  • a probe-type deep body temperature measuring device (not shown) for measuring esophagus temperature as deep body temperature for a subject wearing clothes 500, a temperature sensor 210 for measuring the temperature in the clothes, and an acceleration sensor 220, an atmospheric pressure sensor 230, an angular velocity sensor 240, and a magnetic sensor 250 were mounted.
  • a capsule-type sweat rate measuring device, a mask-type oxygen uptake measuring device, and the like were also used.
  • the subjects were 6 persons indicated by S, H, M, Y, K and U.
  • the probe of the probe type deep body temperature measuring device was attached to the esophagus of the subject.
  • the temperature sensor 210, the acceleration sensor 220, the angular velocity sensor 240, and the magnetic sensor 250 were disposed in the inner pocket 511 provided in the clothes 500.
  • the barometric pressure sensor 230 was disposed in the outer pocket 521 of the garment 500.
  • Clothes 500 worn by the subject are clothes described in Comparative Example 4 of JP-A-2014-091307.
  • This garment 500 is a multilayer structure fabric and is sewn in the form of fire protection clothing (fire protection clothing).
  • polymetaphenylene isophthalamide fiber manufactured by Teijin Ltd., trade name: Cornex
  • copolyparaphenylene-3,4'-oxydiphenylene terephthalamide fiber manufactured by Teijin Techno Products Co., Ltd.
  • a woven fabric composed of heat-resistant fibers (number count: 40/2) made of heat-resistant fibers mixed at a mixing ratio of 90:10 was used to weave a fabric configured as a plain weave ripstop.
  • the basis weight of this outer layer was 380 g / m 2 .
  • polymetaphenylene isophthalamide fiber manufactured by Teijin Ltd., trade name: Cornex
  • copolyparaphenylene-3,4′-oxydiphenylene terephthalamide fiber manufactured by Teijin Techno Products Co., Ltd., trade name: Technora
  • a woven fabric fabric weight: 80 g / m 2 ) woven into a plain weave using spun yarn (count: 40 / ⁇ ) consisting of heat-resistant fibers mixed at a mixing ratio of 95: 5.
  • stacked the moisture-permeable waterproof film made by Japan Gore-Tex company made from fluoro ethylene was employ
  • the heat shield layer (inner layer) is a plain weave using a filament with a total fineness of 1670 dtex made of copolyparaphenylene-3,4'-oxydiphenylene terephthalamide fiber yarn (manufactured by Teijin Techno Products, trade name: Technola).
  • the constructed fabric was woven.
  • the basis weight of the heat shielding layer was 210 g / m 2 .
  • the inclination angle of the treadmill for allowing the subject to walk was set to 5%, and the running speed of the treadmill was gradually increased.
  • the operating speed is 0 km / h for 3 minutes from 0 to 3 minutes from the start of measurement, and 1.5 km / h for 9 minutes from 3 minutes to 12 minutes.
  • the time is 3.0km / h for 9 minutes from 12 minutes to 21 minutes, 4.5km / h for 9 minutes from 21 minutes to 30 minutes, and 6.0km / h for 30 minutes after 30 minutes did.
  • evaluation was terminated when the heart rate reached 180 or higher or the esophagus temperature reached 38.5 degrees.
  • FIG. 6 shows basic data of the subject.
  • Basic data include age, gender, height (cm), body weight (kg), BMI, body surface area (cm 2 ), measurement time (min), post-exercise weight (kg), weight difference (kg).
  • BMI is indicated by weight / height 2 .
  • Body surface area is shown by weight 0.444 x height 0.663 x 88.83 / 10,000.
  • the weight difference is a value obtained by subtracting the weight after exercise from the weight (weight before exercise).
  • the heat production amount was calculated by the method described in the embodiment. However, as the heat generation amount, a method described in Japanese Patent Application Laid-Open No. 2008-220517 or a method described in Japanese Patent No. 3571272 may be adopted.
  • the temperature inside the clothes is the temperature measured by the temperature sensor 210.
  • the result data of Example 1 are shown in FIG. FIG. 8 shows the relationship between measured values and calculated values.
  • the first embodiment is an embodiment in which all parameters are individually calculated.
  • the parameters are four of C1, C2, R1, and R2 in a simultaneous differential equation including six equations from Equation (1) to Equation (6).
  • the measured value of core body temperature rises with the passage of time.
  • core body temperature is indicated by a temperature change amount (core body temperature rise) based on core body temperature at the start of measurement.
  • the mathematical model represented by the above-described simultaneous differential equations is appropriate, the difference between the temporal change in core temperature calculated using the mathematical model and the temporal change in measured core temperature should be small. is there. That is, the correlation coefficient between the calculated value waveform indicating the temporal change of the calculated value of the core temperature and the measured value waveform indicating the temporal change of the measured value of the core temperature should be close to 1.
  • the calculated value of the core temperature is determined to obtain the calculated value waveform while switching between C1, C2, R1, and R2 for each subject. Then, parameters (C1, C2, R1, R2) at which the correlation coefficient between the calculated value waveform and the measured value waveform is closest to 1 were specified.
  • FIG. 7 shows parameters (C1, C2, R1, R2) that maximize the correlation coefficient for each subject.
  • the correlation coefficient is between 1 and 0.7 means that there is a strong positive correlation.
  • a correlation coefficient of between 0.7 and 0.4 means that there is moderate positive correlation.
  • a correlation coefficient of between 0.4 and 0.2 means that there is a weak positive correlation.
  • a correlation coefficient of between 0.2 and -0.2 means that there is almost no correlation.
  • a correlation coefficient between -0.2 and -0.4 means that there is a weak negative correlation.
  • a correlation coefficient between -0.4 and -0.7 means that there is a moderate negative correlation.
  • a correlation coefficient between -0.7 and -1 means that there is a strong negative correlation.
  • Slope is the slope of a straight line approximated by the plotted points when the amount of change in core temperature per minute is plotted with the calculated value as the horizontal axis and the actual value as the vertical axis. Note that this straight line can be determined, for example, by the least squares method. Intercept is the intercept of this straight line. AVG shows an average value. SD indicates standard deviation, that is, variation. SE indicates the standard error, that is, the estimation error of the estimator.
  • Example 2 In Example 2, based on the result data of Example 1, an attempt was made to generalize the parameters. Specifically, in Example 2, the average value of the parameters for each subject determined in Example 1 was employed. The result data of Example 2 are shown in FIG.
  • the correlation coefficient was about 0.97 or more for any subject. That is, according to the second embodiment, it is considered that the mathematical model holds even if values (that is, average values) predetermined in all the parameters (C1, C2, R1, R2) are adopted in the above-described mathematical model.
  • Example 3 In Example 3, based on the result data of Example 1, it was verified whether the parameter was correlated with the weight of the user. Specifically, in Example 3, first, all parameters (C1, C2, R1, R2) were examined to see if there was a correlation with their weight.
  • FIG. 10 shows the result data 1 of the third embodiment.
  • the correlation coefficient is a correlation coefficient between each parameter and body weight.
  • a method of determining the correlation coefficient between each parameter and weight will be described with reference to FIG.
  • FIG. 11 is a diagram showing the relationship between body weight and C1, and is a scatter diagram of each subject when the body weight is taken along the horizontal axis and C1 is taken along the vertical axis.
  • PS is a point specified by the weight of the subject indicated by S (hereinafter simply referred to as “S”. The same applies to other subjects) and C1 after adjustment according to Example 1 of S.
  • Show. PH indicates a point specified by the weight of H and C1 after adjustment of H.
  • PM indicates a point specified by M's weight and C1 after M's adjustment.
  • PY indicates a point specified by the weight of Y and C1 after adjustment of Y.
  • PK indicates a point specified by the weight of K and C1 after adjustment of K.
  • PU shows the point specified by U's weight and C1 after U's adjustment.
  • r which is the correlation coefficient between the weight and C1
  • r is a correlation coefficient between the weight waveform when the weight of all the subjects is arranged in order and the C1 waveform when C1 of all the subjects is arranged in order.
  • r was 0.754 and r 2 was 0.568. Since the correlation coefficient between weight and C1 is 0.754, it is considered that the weight and C1 have a strong positive correlation. In other words, it is considered effective to set C1 using the weight as an index. Therefore, we try to identify the correspondence between weight and C1.
  • L1 represents a straight line approximating PS, PH, PM, PY, PK, and PU. L1 is specified, for example, by the least squares method.
  • the weight is BW
  • C1 0.1234 ⁇ BW + 13.02.
  • FIG. 12 shows the result data 2 of the third embodiment.
  • the correlation coefficient was about 0.97 or more for all subjects. That is, according to the third embodiment, even if C1 and C2 are specified using the weight in the above-described mathematical model, and a predetermined value (that is, an average value) is adopted as R1 and R2, the mathematical operation is performed. It is thought that a model holds.
  • Example 4 In Example 4, based on the result data of Example 1, it was verified whether the parameter was correlated with the body surface area of the user. Specifically, in Example 4, first, all parameters (C1, C2, R1, R2) were examined to see if there was a correlation coefficient with the body surface area.
  • FIG. 13 shows result data 1 of the fourth embodiment.
  • the correlation coefficient is the correlation coefficient between each parameter and the body surface area.
  • the method of determining the correlation coefficient between each parameter and body surface area is basically the same as the method of determining the correlation coefficient between each parameter and body weight.
  • the body surface area is BSA
  • C1 0.0010 ⁇ BSA + 0.7450.
  • FIG. 14 shows the result data 2 of the fourth embodiment.
  • the correlation coefficient was about 0.97 or more for all subjects. That is, according to the fourth embodiment, even if C1 and C2 are specified using BSA in the above-described mathematical model, and a predetermined value (that is, an average value) is adopted as R1 and R2, the mathematical operation is performed. It is thought that a model holds.
  • Example 5 (Example 5)
  • C1 was specified by body weight
  • C2 was specified by body surface area
  • the result data of Example 5 are shown in FIG.
  • the T value is a statistical test quantity.
  • N is the number of data.
  • P value is the two-sided probability of T distribution.
  • C1 is identified using BW
  • C2 is identified using BSA
  • the average value is adopted for R1 and R2
  • the correlation coefficient is about 0.97 or more for any subject. It became. That is, according to the fifth embodiment, in the mathematical model described above, C1 is identified using BW, C2 is identified using BSA, and predetermined values (that is, average values) for R1 and R2 are determined. Even if adopted, it is considered that a mathematical model holds.
  • FIG. 16 shows the correlation coefficient, the SLOPE, and the Intercept under four conditions.
  • the four conditions are a first condition for identifying C1 by body weight and C2 by body surface area, a second condition for identifying C1 by body weight and adopting an average value as C2, and C1 and C2 by body weight.
  • a third condition to be specified and a fourth condition to specify C1 and C2 by body surface area are included.
  • predetermined values that is, average values
  • the first condition corresponds to the fifth embodiment
  • the third condition corresponds to the third embodiment
  • the fourth condition corresponds to the fourth embodiment.
  • C1 and C2 are specified by the value (for example, weight and body surface area) based on the physical feature of the user, and predetermined values (that is, average values) for R1 and R2 are determined. If adopted, it is considered that a mathematical model is established. In addition, considering the correlation coefficient, SLOPE, and Intercept, C1 is specified by the weight, C2 is specified by the body surface area as in the first condition (Example 5), and R1 and R2 are made in advance. It is considered to be most preferable to adopt a defined value (that is, an average value).
  • the temperature in the clothes mentioned above includes the air temperature between the clothes and the human body, the air temperature between the clothes worn in the outermost layer and the clothes inside the clothes, the skin temperature of the human body covered by the clothes, etc. However, it is not limited thereto.
  • the physical characteristics of the user include, but are not limited to, weight, height, BMI, body fat rate, sweat rate, respiratory rate, body surface area, muscle mass and the like.
  • the notification means gives notification, for example, when the deep body temperature exceeds a predetermined threshold or when the degree of increase in the deep body temperature exceeds a predetermined increase limit, but is not limited thereto.
  • the first embodiment describes an example in which many functions of the processing apparatus 100 are realized by software (or firmware), that is, many functions of the processing apparatus 100 are realized by execution of a program by the processor.
  • a function may be realized by hardware.
  • the processing device 100 includes a processing circuit instead of the CPU 11.
  • the processing circuit is configured by, for example, a single circuit, a complex circuit, a programmed processor, a parallel programmed processor, an application specific integrated circuit (ASIC), a field programmable gate array (FPGA), or a combination thereof.
  • the personal computer it is also possible to cause the personal computer to function as the processing device 100 according to the present invention by applying an operation program defining the operation of the processing device 100 according to the present invention to an existing personal computer or information terminal device.
  • the distribution method of such a program is arbitrary, and for example, it is stored by being stored in a computer readable recording medium such as a CD-ROM (Compact Disk Read-Only Memory), a DVD (Digital Versatile Disk), or a memory card. Or may be distributed via a communication network (eg, the Internet).

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Physiology (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Geometry (AREA)
  • Business, Economics & Management (AREA)
  • Emergency Management (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Professional, Industrial, Or Sporting Protective Garments (AREA)
  • Details Of Garments (AREA)
  • Emergency Alarm Devices (AREA)

Abstract

温度センサ（２１０）は、ユーザが着用する衣服の内側の温度である衣服内温度を測定する。加速度センサ（２２０）は、ユーザに加わる加速度を検出する。処理装置（１００）は、温度センサ（２１０）により測定された衣服内温度と加速度センサ（２２０）により検出された加速度とに基づいて、ユーザの体内の温度である深部体温を推定する。

Description

深部体温推定システム、ヒートストレス警報システム、及び、深部体温推定方法
　本発明は、深部体温推定システム、ヒートストレス警報システム、及び、深部体温推定方法に関する。
　消火活動を実施する消防士や炎天下で作業する作業者などのユーザが熱中症になることを防止することは重要である。熱中症を防止するためには、例えば、熱中症を引き起こすヒートストレスを検出し、この検出結果に応じて警報を発する方法が考えられる。
　例えば、特許文献１には、温度センサにより検出された防護装備の内側温度に基づいて、ヒートストレス状況を把握する技術が記載されている。特許文献１に記載された技術では、防護装備の内側温度を把握することで熱中症の予測ができるという前提のもと、温度センサにより防護装備の内側温度を検出している。
特開２００４－３０１８０号公報
　ところで、熱中症と深部体温との間には関係があることは知られている。すなわち、深部体温は熱中症をはじめとする、人体が感じるヒートストレスの危険性を示す指標となりうるものである。しかし、先行技術では深部体温を直接予測するものではなく、熱中症の予測精度としては十分なものではなかった。また、深部体温を測定するためには温度測定プローブを口腔、耳、あるいは肛門より体内に挿入する必要があり、ユーザの作業中に深部体温を測定することは困難であった。
　本発明は、上記問題に鑑みてなされたものであり、ユーザの深部体温を精度良く推定することが可能な深部体温推定システム、ヒートストレス警報システム、及び、深部体温推定方法を提供することを目的とする。
　上記目的を達成するために、本発明の第１の観点に係る深部体温推定システムは、
　ユーザが着用する衣服の内側の温度である衣服内温度を測定する温度センサと、
　前記ユーザに加わる加速度を検出する加速度センサと、を備え、
　前記ユーザの体内の温度である深部体温を推定する。
　前記温度センサにより測定された衣服内温度と前記加速度センサにより検出された加速度とに基づいて、前記深部体温を推定する深部体温推定手段を更に備えていてもよい。
　前記加速度センサにより検出された加速度に基づいて、前記ユーザの体内で産出される熱量である産熱量を推定する産熱量推定手段を更に備え、
　前記深部体温推定手段は、前記温度センサにより測定された衣服内温度と前記産熱量推定手段により推定された産熱量とに基づいて、前記深部体温を推定してもよい。
　前記深部体温推定手段は、前記衣服内温度と前記産熱量と前記深部体温との対応関係を示す数理モデルを用いて、前記深部体温を推定してもよい。
　前記数理モデルは、
　前記ユーザの筋肉が有する熱量である第１の熱量の変化量を、前記産熱量と前記ユーザの筋肉から前記ユーザの血液に供給される熱流量である第１の熱流量との差で示す第１の式と、
　前記ユーザの血液が有する熱量である第２の熱量の変化量を、前記第１の熱流量と前記ユーザの血液から前記衣服の内側に供給される熱流量である第２の熱流量との差で示す第２の式と、
　前記第１の熱流量を、前記ユーザの筋肉から前記ユーザの血液への熱抵抗値である第１の熱抵抗値に対する、前記ユーザの筋肉の温度である第１の温度の変化量と前記ユーザの血液の温度である第２の温度の変化量との差の割合で示す第３の式と、
　前記第２の熱流量を、前記ユーザの血液から前記衣服の内側への熱抵抗値である第２の熱抵抗値に対する、前記第２の温度の変化量と前記衣服内温度の変化量との差の割合で示す第４の式と、
　前記第１の温度の変化量を、前記ユーザの筋肉が有する熱容量である第１の熱容量に対する前記第１の熱量の割合で示す第５の式と、
　前記第２の温度の変化量を、前記ユーザの血液が有する熱容量である第２の熱容量に対する前記第２の熱量の割合で示す第６の式と、を含む連立微分方程式により表され、
　前記深部体温推定手段は、前記第２の温度を前記深部体温として推定してもよい。
　前記第１の熱容量と前記第２の熱容量とは、前記ユーザの身体上の特徴に基づく値であり、
　前記第１の熱抵抗値と前記第２の熱抵抗値とは、予め定められた値であってもよい。
　前記ユーザの周囲の気圧を検出する気圧センサを更に備え、
　前記産熱量推定手段は、前記加速度センサにより検出された加速度と前記気圧センサにより検出された気圧とに基づいて、前記産熱量を推定してもよい。
　前記ユーザに加わる角速度を検出する角速度センサと、
　前記ユーザの周囲の磁場の方向を検出する磁気センサと、を更に備え、
　前記加速度センサと前記角速度センサと前記磁気センサとは、相互に固定されており、
　前記産熱量推定手段は、前記加速度センサにより検出された加速度と前記角速度センサにより検出された角速度と前記磁気センサにより検出された方向とに基づいて、前記産熱量を推定してもよい。
　前記深部体温推定手段により推定された深部体温に基づいて、前記ユーザの異常を報知する報知手段を更に備えてもよい。
　前記衣服が、遮熱効果を有する衣服であってもよい。
　前記遮熱効果が、ＩＳＯ９１５１に準拠した方法で測定されるＨＴＩ２４が３秒以上であることを満たしてもよい。
　前記衣服が、ＪＩＳ　Ｌ　１０９６に準拠した方法であってフラジール形試験機を用いた方法で測定される通気度の値が１．０ｃｍ３／ｃｍ２／ｓ以下であってもよい。
　前記衣服が、ＪＩＳ　Ｌ　１０９９に準拠した方法であってカップ法により測定される透湿度の値が１０００ｇ／ｍ２／ｈ以下であってもよい。
　前記衣服の熱伝導率が、１ｋｃａｌ／ｈ／ｍ／℃以下であってもよい。
　上記目的を達成するために、本発明の第２の観点に係るヒートストレス警報システムは、
　ユーザが着用する衣服の内側の温度である衣服内温度と前記ユーザに加わる加速度とに基づいて、前記ユーザの体内の温度である深部体温を推定する深部体温推定手段と、
　前記深部体温推定手段により推定された深部体温に基づいて、警報を発する警報手段と、を備える。
　上記目的を達成するために、本発明の第３の観点に係る深部体温推定方法は、
　ユーザが着用する衣服の内側の温度である衣服内温度を測定し、
　前記ユーザに加わる加速度を検出し、
　前記測定された衣服内温度と前記検出された加速度とに基づいて、前記ユーザの体内の温度である深部体温を推定する。
　本発明によれば、ユーザの深部体温を精度良く推定することができる。
本発明の実施形態に係る深部体温推定システムの構成図 本発明の実施形態に係る処理装置の構成図 各モジュールの格納場所を示す図 数理モデルの説明図 本発明の実施形態に係る処理装置が実行する深部体温推定処理を示すフローチャート 被験者の基本データを示す図 実施例１の結果データを示す図 実測値と算出値との関係を示す図 実施例２の結果データを示す図 実施例３の結果データ１を示す図 体重とＣ１との関係を示す図 実施例３の結果データ２を示す図 実施例４の結果データ１を示す図 実施例４の結果データ２を示す図 実施例５の結果データを示す図 実施例の最終結果データを示す図
　以下、本発明の実施形態について、図面を参照して説明する。なお、図中において、同一又は対応する部分には、同一の符号を付す。
（実施形態）
　図１は、本発明の実施形態に係る深部体温推定システム１０００の構成図である。深部体温推定システム１０００は、ユーザの体の中心部の体温（以下「深部体温」という。）を推定するシステムである。より詳細には、深部体温推定システム１０００は、ユーザの動作により発生する熱量（以下「産熱量」という。）と、ユーザとユーザが着用する衣服との間の空間の温度（以下「衣服内温度」という。）と、を用いた数理モデルを用いて、ユーザの深部体温を推定するシステムである。
　また、深部体温推定システム１０００は、深部体温を推定する機能に加え、深部体温から熱中症のおそれがある場合に警報を発する機能を有する。本実施形態では、ユーザは、消防活動を実施する消防士であり、衣服として防火服を着用しているものとする。ただし、ユーザは、消防士ではなく、炎天下で作業する作業者であってもよい。また、ユーザが着用する衣服は、防火服に限定されない。ただし、ユーザが着用する衣服は、一般の衣服よりも、遮熱効果が高い衣服であることが望ましい。
　ここで言う遮熱効果の指標である遮熱性とは、例えば衣服に用いられている布帛が、ＩＳＯ９１５１に準拠した方法により、規定の火炎に暴露し、温度上昇が２４℃に達するまでの時間（ＨＴＩ２４）が３秒以上であることを満たすことが挙げられるが、これに限定されない。
　また、遮熱効果が高い衣服としては、例えば特開２０１４－０９１３０７号公報や特開２０１１－１０６０６９号公報に記載された多層構造布帛を有する衣服が挙げられる。すなわち、外層および内層の２層以上からなる積層布帛を有する衣服であることが好適である。積層布帛の各層の厚みは、遮熱性に大きく影響する。このため、例えば、特開２０１０－２５５１２４号公報に記載されているように、表地層の厚みと内層の厚みは、下記式を満足する布帛を用いた衣服であることが好ましいが、これに限定されない。
５．０ｍｍ≧遮熱層厚み（ｍｍ）≧－２９．６×（表地層厚み（ｍｍ））＋１４．１（ｍｍ）
　多層構造布帛においては、火炎暴露下において定常時と布帛形態が変化するものであっても良い。例えば、火炎暴露下において布帛厚みが増加することが考えられる。また、衣服を構成する布帛の遮熱性は、ＩＳＯ９１５１に準拠した方法で測定してＨＴＩ２４が３秒以上であることが好ましい。ＨＴＩ２４は、規定の火炎に暴露した場合において、温度上昇が２４℃に達するまでの時間である。
　また、多層構造布帛においては、上記温度センサを保護する為、難燃性の高い繊維素材であることが好ましい。例えば、多層構造布帛を構成する繊維の限界酸素指数（ＬＯＩ）が２１以上であり、好ましくは２４以上である。限界酸素指数とは燃焼継続するのに必要な雰囲気の酸素濃度（％）であり、２１以上であると通常の空気中では燃焼が継続せずに自己消火することを意味し、高い耐熱性を発揮することができる。ここで、限界酸素指数（ＬＯＩ）は、ＪＩＳ　Ｌ１０９１（Ｅ法）により測定された値である。
　このように、最外層が限界酸素指数（ＬＯＩ）２１以上の繊維を用いることにより、高い耐熱性を発揮することができる。上記繊維としては、例えば、メタ型アラミド繊維、パラ型アラミド繊維、ポリベンゾイミダゾール繊維、ポリイミド繊維、ポリアミドイミド繊維、ポリエーテルイミド繊維、ポリアリレート繊維、ポリパラフェニレンベンゾビスオキサゾール繊維、ノボロイド繊維、ポリクラール繊維、難燃アクリル繊維、難燃レーヨン繊維、難燃ポリエステル繊維、難燃綿繊維、難燃ウール繊維などを挙げることができる。
　特に、ポリメタフェニレンイソフタルアミドなどのメタ系アラミド繊維や、織物や編物強度を向上させる目的でパラ系のアラミド繊維、すなわち、ポリパラフェニレンテレフタルアミド、あるいは、これに第３成分を共重合した繊維などを用いることが有用である。ポリパラフェニレンテレフタルアミド共重合体の一例として、コポリパラフェニレン・３．４’－オキシジフェニレンテレフタルアミドが例示される。ただし、難燃性を阻害しない範囲で、ポリエステル繊維やポリアミド繊維、ナイロン繊維、アクリル繊維などの易撚素材を混用しても構わない。また、繊維は原着繊維であっても、後染め繊維であっても構わない。また、必要に応じて、布帛を製織後に難燃加工を付与しても構わない。
　上述した繊維には、長繊維または短繊維を用いてもよい。また、上記繊維を２種類以上混繊または混紡して用いてもよい。特に、外層に用いられる布帛として、本発明において、メタ系アラミド繊維とパラ系アラミド繊維はフィラメントあるいは混紡して紡績糸の形態で使用するものが好ましく例示される。使用される紡績糸はシングルプライ、ダブルプライであってもよい。該パラ系のアラミド繊維の混合比率としては、布帛を構成する全繊維に対して５重量％以上が好ましいが、パラ系のアラミド繊維は、フィブリル化を起こしやすいため、混合比率を６０重量％以下に抑えることが好ましい。
　上記布帛は、織物、編物、不織布などの形態で用いてもよいが、特に織物が好ましい。また、織物としては、平織、綾織、朱子織など、どのような織組織であってもよい。また、織物、編物では、２種類の繊維を用い、交織、交編させてもよい。なお、最外層（表地層）に用いる布帛は、目付けが、好ましくは１４０～５００ｇ／ｍ２、より好ましくは１６０～４００ｇ／ｍ２、さらに好ましくは２００～４００ｇ／ｍ２の範囲にあるものを使用することが好ましい。上記目付けが１４０ｇ／ｍ２未満の場合には、十分な耐熱性能が得られないおそれがあり、一方該目付けが５００ｇ／ｍ２を超える場合には、遮熱活動服にした場合の着用感が阻害されるおそれがある。
　多層構造布帛においては、内層が、引張弾性率が８０～８００ｃＮ／ｄｔｅｘ、該布帛の熱伝導率が６．０Ｗ・ｍ－１・ｋ－１以下、好ましくは５．０Ｗ・ｍ－１・ｋ－１以下、かつ密度が３．０ｇ／ｃｍ３以下の繊維で構成された布帛であり、８００～３０００ｎｍの波長の電磁波の透過率が１０％以下であり、かつ目付けが６０～５００ｇ／ｍ２であることが好ましい。
　また、繊維の引張弾性率は、８０～８００ｃＮ／ｄｔｅｘ（より好ましくは８０～５００ｃＮ／ｄｔｅｘ、さらに好ましくは１２０～５００ｃＮ／ｄｔｅｘ）であることが好ましい。引張弾性率が８０ｃＮ／ｄｔｅｘ未満の場合、遮熱活動服などとして使用する場合に、着用者の動きや姿勢によっては一部の部位で繊維が伸びて、布帛が薄くなり十分な遮熱効果が得られない場合がある。また、引張弾性率が８００ｃＮ／ｄｔｅｘを超える場合、いわゆる「突っ張る」着用感を与える場合がある。紡績糸を用いることによりこれを避けることも可能な場合もあるが、十分な効果を発揮するには引張弾性率が８００ｃＮ／ｄｔｅｘ以下であることが好ましい。
　多層構造布帛において、布帛の目付けは６０～５００ｇ／ｍ２（より好ましくは８０～４００ｇ／ｍ２、さらに好ましくは１００～３５０ｇ／ｍ２）であることが好ましい。目付けが６０ｇ／ｍ２より低いと電磁波の透過を十分に防ぐことができない場合がある。一方、目付けが５００ｇ／ｍ２より高いと、熱を溜め込む傾向が顕著になり遮熱性を阻害するおそれがあり、軽量性も損なわれるおそれがある。
　上記の多層構造布帛を構成する繊維には、特に限定を受けるものではない。電磁波の吸収や反射を向上させるために、金属やカーボンなどを練り込んだり、表面付着させたりしたものを用いることも可能である。上記繊維としては、炭素繊維を用いることもができるが、アラミド繊維、ポリベンゾイミダゾール繊維、ポリイミド繊維、ポリアミドイミド繊維、ポリエーテルイミド繊維、ポリアリレート繊維、ポリパラフェニレンベンゾビスオキサゾール繊維、ノボロイド繊維、ポリクラール繊維、難燃アクリル繊維、難燃レーヨン繊維、難燃ポリエステル繊維、難燃綿繊維、難燃ウール繊維などの有機高分子からなる繊維（以下、有機高分子繊維と称することがある）を好適に挙げることができる。
　多層構造布帛においては、電磁波吸収率と熱伝導率の改良のため、カーボン、金、銀、銅、アルミニウムなどの微粒子を、有機高分子繊維に含有させたり、有機高分子繊維の表面に付着させたりすることができる。この際、カーボン等は、これを含有する顔料または塗料として有機高分子繊維に含有または表面に付与等してもよい。これらの微粒子の有機高分子繊維の全重量に対する含有率または付着率は、微粒子の比重等にもよるが、好ましくは０．０５～６０重量％、より好ましくは０．０５～４０重量％である。また、カーボン微粒子の場合は、好ましくは０．０５重量％以上、より好ましくは０．０５～１０重量％、さらに好ましくは０．０５重量％以上５重量％未満である。さらに、アルミニウム微粒子の場合は、好ましくは１重量％以上、より好ましくは１～２０重量％、さらに好ましくは１～１０重量％である。
　また、上記微粒子の数平均粒子径は、好ましくは１０μｍ以下（より好ましくは０．０１～１μｍ）である。炭素繊維や金属繊維など、それ自体が上記のＬＯＩ値や熱伝導性等の要件を満たせば、微粒子を練り込んだりすることなく、そのまま用いることもできる。特に、内層を構成する繊維としては、炭素繊維や金属繊維の含有率が好ましくは５０重量％以上、より好ましくは８０重量％以上、さらに好ましくは１００重量％からなる布帛を例として挙げることができる。
　多層構造布帛には、上記の外層と内層の間に、中間層として、ＬＯＩ値が２５以上の繊維からなる布帛に透湿防水性フィルムを積層固着したものを配することも可能である。これにより布帛構造体としての快適性を保持したまま外からの水の浸入を抑えることができるため、放水などが行われる消防活動を行う消防隊員用の防護服としてより好適である。用いる中間層の目付けは、５０～２００ｇ／ｍ２の範囲にあるものを使用することが好ましい。目付が５０ｇ／ｍ２未満の場合には、十分な遮熱性能が得られない恐れがあり、一方該目付けが２００ｇ／ｍ２を超える場合には、遮熱活動服にした場合の着用感が阻害されるおそれがある。
　この布帛は、透湿防水性のあるポリテトラフルオロエチレン等からなる薄膜フィルムをラミネート加工されていることが好ましく、これにより透湿防水性や耐薬品性が向上し、着用者の汗の蒸散を促進することができ、着用者のヒートストレスを減少することができる。上記の中間層にラミネート加工する薄膜フィルムの単位面積あたりの目付けは１０～５０ｇ／ｍ２の範囲とすることが好ましい。なお、このように中間層の布帛に薄膜フィルムをラミネート加工する場合でも、該加工を施した中間層の布帛の目付けが前述した５０～２００ｇ／ｍ２の範囲にあることが好ましい。
　さらに、本実施形態においては通気度の低い生地を用いられていることが望ましい。望ましくはＪＩＳ　Ｌ　１０９６（フラジール形試験機）で測定した値が１．０ｃｍ３／ｃｍ２／ｓ以下であり、さらに望ましくは０．０１ｃｍ３／ｃｍ２／ｓ以下である。
　さらに、本実施形態においては透湿度の低い生地を用いられていることが望ましい。望ましくはＪＩＳ　Ｌ　１０９９（カップ法）で測定した値が１０００ｇ／ｍ２／ｈ以下であり、さらに望ましくは１０ｇ／ｍ２／ｈ以下である。
　さらに、本実施形態においては熱伝導率の低い生地を用いられていることが望ましい。望ましくは１ｋｃａｌ／ｈ／ｍ／℃以下であり、さらに望ましくは１ｃａｌ／ｈ／ｍ／℃以下である。
　図１に示すように、深部体温推定システム１０００は、処理装置１００と、温度センサ２１０と、加速度センサ２２０と、気圧センサ２３０と、角速度センサ２４０と、磁気センサ２５０と、ブザー３００と、端末装置４００と、を備える。処理装置１００と各種のセンサ（温度センサ２１０、加速度センサ２２０、気圧センサ２３０、角速度センサ２４０、磁気センサ２５０）とは相互に接続される。また、処理装置１００とブザー３００とは相互に接続され、処理装置１００と端末装置４００とは図示しないネットワークを介して相互に接続される。
　処理装置１００は、深部体温推定システム１０００の中核をなすモジュールである。処理装置１００は、例えば、マイクロコンピュータ、パーソナルコンピュータである。処理装置１００は、予め定められた数理モデルを用いて、各種のセンサから取得した情報からユーザの深部体温を推定する。また、処理装置１００は、ユーザが熱中症になる可能性が高いと判別した場合、つまり、推定した深部体温が異常であると判別した場合、音、光、画面表示により警報を発する。処理装置１００は、自ら警報を発してもよいし、ブザー３００や端末装置４００に警報を発させてもよい。以下、図２を参照して、処理装置１００の構成について説明する。
　図２に示すように、処理装置１００は、ＣＰＵ（Central Processing Unit）１１、ＲＯＭ（Read Only Memory）１２、ＲＡＭ（Random Access Memory）１３、フラッシュメモリ１４、ＲＴＣ（Real Time Clock）１５、タッチスクリーン１６、第１の通信インターフェース１７、第２の通信インターフェース１８を備える。処理装置１００が備える各構成要素は、バスを介して相互に接続される。
　ＣＰＵ１１は、処理装置１００の全体の動作を制御する。なお、ＣＰＵ１１は、ＲＯＭ１２に格納されているプログラムに従って動作し、ＲＡＭ１３をワークエリアとして使用する。ＲＯＭ１２には、処理装置１００の全体の動作を制御するためのプログラムやデータが記憶される。ＲＡＭ１３は、ＣＰＵ１１のワークエリアとして機能する。つまり、ＣＰＵ１１は、ＲＡＭ１３にプログラムやデータを一時的に書き込み、これらのプログラムやデータを適宜参照する。深部体温推定手段の機能と産熱量推定手段の機能とは、例えば、ＣＰＵ１１がＲＡＭ１３をワークエリアとして使用してＲＯＭ１２に格納されたプログラムを実行することにより実現される。
　フラッシュメモリ１４は、各種の情報を記憶する不揮発性メモリである。ＲＴＣ１５は、計時用のデバイスである。ＲＴＣ１５は、例えば、電池を内蔵し、処理装置１００の電源がオフの間も計時を継続する。ＲＴＣ１５は、例えば、水晶発振子を備える発振回路を備える。タッチスクリーン１６は、ユーザによりなされたタッチ操作を検知し、検知の結果を示す信号をＣＰＵ１１に供給する。また、タッチスクリーン１６は、ＣＰＵ１１から供給された画像信号に基づく画像を表示する。このように、タッチスクリーン１６は、処理装置１００のユーザインターフェースとして機能する。
　第１の通信インターフェース１７は、処理装置１００を、図示しないネットワーク（例えば、インターネット）に接続するためのインターフェースである。処理装置１００は、図示しないネットワークを介して端末装置４００と通信する。第１の通信インターフェース１７は、例えば、無線ＬＡＮ（Local Area Network）用のインターフェースである。
　第２の通信インターフェース１８は、処理装置１００を、各種のセンサやブザー３００と接続するためのインターフェースである。第２の通信インターフェース１８は、図示しない通信ケーブルを介して各種のセンサやブザー３００に接続されてもよいし、無線通信により各種のセンサやブザー３００に接続されてもよい。第２の通信インターフェース１８は、無線ＬＡＮ用のインターフェースであってもよいし、ＵＳＢ（Universal Serial Bus）やＲＳ－２３２（Recommended Standard 232）に準拠したインターフェースであってもよい。
　温度センサ２１０は、周囲の温度を測定するセンサである。温度センサ２１０は、例えば、衣服の内部に配置され、衣服内温度を測定する。温度センサ２１０は、測定した衣服内温度を示す情報（以下「温度情報」という。）を処理装置１００に供給する。温度センサ２１０は、例えば、水銀式温度計、アルコール式温度計、サーミスタ式温度計、赤外線式温度計、測温抵抗体、ＣＭＯＳ（ＣｏｍｐｌｅｍｅｍｙａｒｙＭＯＳ）温度センサである。
　加速度センサ２２０は、加速度センサ２２０に加わる加速度を測定するセンサである。加速度センサ２２０は、ユーザや衣服に装着され、ユーザに加わる加速度、つまり、ユーザの動作により加速度センサ２２０に生じる加速度を測定する。加速度センサ２２０は、ＸＹＺ軸の３軸の加速度を測定する加速度センサであるものとする。加速度センサ２２０は、Ｘ軸の加速度とＹ軸の加速度とＺ軸の加速度とを示す情報（以下「加速度情報」という。）を、処理装置１００に供給する。加速度センサ２２０は、例えば、半導体方式の加速度センサ、機械的変位測定方式の加速度センサ、光学的方式の加速度センサである。半導体方式の加速度センサは、ＭＥＭＳ（Micro Electro Mechanical Systems）技術を用いた加速度センサであり、例えば、静電容量型の加速度センサ、ピエゾ抵抗型の加速度センサ、ガス温度分布型の加速度センサである。
　気圧センサ２３０は、周囲の気圧を測定するセンサである。気圧センサ２３０は、例えば、衣服の外部に配置され、ユーザの周囲の気圧を測定する。気圧センサ２３０は、測定した気圧を示す情報（以下「気圧情報」という。）を処理装置１００に供給する。気圧センサ２３０は、例えば、半導体を用いた電気式気圧センサであり、静電容量式の気圧センサ、振動式のセンサである。
　角速度センサ２４０は、角速度センサ２４０に加わる角速度を測定するセンサである。角速度センサ２４０は、例えば、ユーザや衣服に装着され、ユーザに加わる角速度、つまり、ユーザの動作により角速度センサ２４０に生じる角速度を測定する。角速度センサ２４０は、ＸＹＺ軸の３軸の加速度を測定する加速度センサであるものとする。角速度センサ２４０は、測定した角速度を示す情報（以下「角速度情報」という。）を、処理装置１００に供給する。角速度センサ２４０は、例えば、水晶発振子や圧電セラミック素子を用いて、コリオリの力によって生じるひずみを電気信号として取り出すジャイロセンサである。
　磁気センサ２５０は、周囲の磁場の大きさや方向を検出するセンサである。磁気センサ２５０は、例えば、ユーザや衣服に装着され、地磁気（地球により生じる磁場）を検出し、方角を検出する。磁気センサ２５０は、検出した方角を示す情報（以下「方角情報」という。）を、処理装置１００に供給する。磁気センサ２５０は、例えば、ホール素子、磁気抵抗効果素子、磁気インピーダンス素子を備える。
　ブザー３００は、処理装置１００による制御に従って、警報を発する。この警報は、例えば、推定された深部体温が予め定められた閾値よりも高い場合、あるいは、推定された深部体温が予め定められた値を超えて上昇した（例えば、１．５℃以上）場合、ユーザが熱中症になる可能性があることをユーザに知らせる音である。ブザー３００は、電磁石と振動板とを備え、処理装置１００から供給された電気信号に従って、音を発する。警報手段の機能は、例えば、処理装置１００とブザー３００とが協働することにより実現される。
　端末装置４００は、処理装置１００のユーザインターフェースとして機能する装置である。端末装置４００は、例えば、処理装置１００による制御に従って、音又は画面表示により警報を発する。端末装置４００は、例えば、図示しない、ＣＰＵ、ＲＯＭ、ＲＡＭ、フラッシュメモリ、通信インターフェースを備える。端末装置４００は、例えば、スマートフォン、タブレット端末、スマートウォッチ、ヘッドマウントディスプレイである。
　次に、図３を参照して、処理装置１００、各種のセンサ、ブザー３００の配置例について説明する。図３は、処理装置１００、各種のセンサ、ブザー３００が、衣服５００に装着されている様子を示す図である。図３には、衣服５００が、裏地５１０と表地５２０とを備え、裏地５１０と表地５２０との間に、空間５３０を有する例が示されている。衣服５００は、空間５３０を有することにより、衣服５００の内側の空間（衣服５００とユーザが着用する下着６００との間の空間）である空間５４０と、衣服５００の外側の空間である空間５５０との間の断熱性を高めている。また、表地５２０は、断熱性の高い素材で構成される。
　裏地５１０には、内ポケット５１１が設けられ、表地５２０には、外ポケット５２１が設けられている。処理装置１００、各種のセンサ、ブザー３００は、例えば、内ポケット５１１と外ポケット５２１とのいずれか一方に収容される。温度センサ２１０は、衣服内温度を測定するために、内ポケット５１１に収容される。気圧センサ２３０は、衣服５００の内側の気圧よりも衣服５００の外側の気圧を測定する方が好ましいため、外ポケット５２１に収容される。ブザー３００は、警報がユーザに伝わりやすいように、外ポケット５２１に収容される。
　処理装置１００と加速度センサ２２０と角速度センサ２４０と磁気センサ２５０とは、内ポケット５１１と外ポケット５２１とのいずれに収容されてもよい。ただし、加速度センサ２２０により検出される加速度の軸や角速度センサ２４０により検出される角速度の軸は、磁気センサ２５０により検出された方向を基準として決定される。このため、加速度センサ２２０と角速度センサ２４０と磁気センサ２５０とは、互いの位置関係が変化しないように、互いに固定される。本実施形態では、処理装置１００と温度センサ２１０と加速度センサ２２０と角速度センサ２４０と磁気センサ２５０とが、互いの位置関係が変化しないように、互いに固定されるものとする。
　次に、図４を参照して、深部体温の推定に用いられる数理モデルについて説明する。数理モデルは、以下に示す式（１）から式（６）までの６つの式により構成される連立微分方程式により表現される。第１の式は、例えば、式（１）に対応する。第２の式は、例えば、式（２）に対応する。第３の式は、例えば、式（３）に対応する。第４の式は、例えば、式（４）に対応する。第５の式は、例えば、式（５）に対応する。第６の式は、例えば、式（６）に対応する。
　ｄＱ１／ｄｔ＝Ｆ０－Ｆ１　式（１）
　ｄＱ２／ｄｔ＝Ｆ１－Ｆ２　式（２）
　Ｆ１＝（Ｔ１－Ｔ２）／Ｒ１　式（３）
　Ｆ２＝（Ｔ２－Ｔ３）／Ｒ２　式（４）
　Ｔ１＝Ｑ１／Ｃ１　式（５）
　Ｔ２＝Ｑ２／Ｃ２　式（６）
　式（１）は、第１の熱量の変化量を、産熱量と第１の熱流量との差で表す式である。第１の熱量は、ユーザの筋肉が有する熱量である。第１の熱流量は、ユーザの筋肉からユーザの血液に供給される熱流量である。Ｑ１は、第１の熱量であり、ｄＱ１／ｄｔは、単位時間当たりの第１の熱量の変化量である。単位時間は、例えば、１分間、１０秒間、１秒間である。以下、単位時間は、１分間を意味するものとする。Ｆ０は、産熱量である。産熱量に関しては、後述する。Ｆ１は、第１の熱流量である。式（１）は、筋肉が有する熱量の変化量は、筋肉に入ってくる熱流量である産熱量と、筋肉から血液に出ていく熱流量との差で定義されることを示している。
　式（２）は、第２の熱量の変化量を、第１の熱流量と第２の熱流量との差で表す式である。第２の熱量は、ユーザの血液が有する熱量である。第２の熱流量は、ユーザの血液から衣服５００の内側に供給される熱流量である。Ｑ２は、第２の熱量であり、ｄＱ２／ｄｔは、単位時間当たりの第２の熱量の変化量である。Ｆ２は、第２の熱流量である。式（２）は、血液が有する熱量の変化量は、筋肉から血液に入ってくる熱流量と、血液から衣服５００の内側に出ていく熱流量との差で定義されることを示している。
　式（３）は、第１の熱流量を、第１の熱抵抗値に対する、第１の温度の変化量と第２の温度の変化量との差の割合で表す式である。第１の熱抵抗値は、ユーザの筋肉からユーザの血液への熱抵抗値であり、筋肉から血液への熱の伝わりにくさを示す値である。第１の温度は、ユーザの筋肉の温度である。第２の温度は、ユーザの血液の温度である。Ｒ１は、第１の熱抵抗値である。Ｔ１は、単位時間当たりの第１の温度の変化量である。Ｔ２は、単位時間当たりの第２の温度の変化量である。式（３）は、筋肉から血液に流れる熱流量が、筋肉の温度の変化量と血液の温度の変化量との差と、筋肉から血液への熱の伝わりにくさとの比で定義されることを示している。
　式（４）は、第２の熱流量を、第２の熱抵抗値に対する、第２の温度の変化量と衣服内温度の変化量との差の割合で表す式である。第２の熱抵抗値は、ユーザの血液から衣服５００の内側への熱抵抗値であり、血液から衣服の内側への熱の伝わりにくさを示す値である。Ｒ２は、第２の熱抵抗値である。Ｔ３は、単位時間当たりの衣服内温度の変化量である。式（４）は、筋肉から血液に流れる熱流量が、筋肉の温度の変化量と血液の温度の変化量との差と、筋肉から血液への熱の伝わりにくさとの比で定義されることを示している。
　式（５）は、第１の温度の変化量を、第１の熱容量に対する第１の熱量の割合で示す式である。第１の熱容量は、ユーザの筋肉が有する熱容量である。Ｃ１は、第１の熱容量である。式（５）は、筋肉の温度変化が、筋肉の熱容量と筋肉が有する熱量との比で定義されることを示している。
　式（６）は、第２の温度の変化量を、第２の熱容量に対する第２の熱量の割合で示す式である。第２の熱容量は、ユーザの血液が有する熱容量である。Ｃ２は、第２の熱容量である。式（６）は、血液の温度変化が、血液の熱容量と血液が有する熱量との比で定義されることを示している。
　図４は、筋肉において産熱量（Ｆ０）の熱が発生すると、第１の熱抵抗値（Ｒ１）に応じた速度で、筋肉から血液に向けて第１の熱流量（Ｆ１）の熱が流れ、第２の熱抵抗値（Ｒ２）に応じた速度で、血液から衣服５００の内側に向けて第２の熱流量（Ｆ２）の熱が流れる様子を示している。ここで、産熱量と衣服内温度とに基づいて、第２の温度が推定され、第２の温度が深部体温として推定される。つまり、本実施形態では、ユーザの血液の温度である第２の温度が、ユーザの深部体温と見做される。
　なお、本実施形態では、Ｆ０が加速度情報と気圧情報とから求められ、Ｔ３が温度情報から求められ、Ｒ１、Ｒ２、Ｃ１、Ｃ２が固定値であり、Ｑ１、Ｑ２、Ｆ１、Ｆ２、Ｔ１、Ｔ２が変数であるものとする。上記固定値は、例えば、深部体温の算出値の時間変化を示す算出値波形と、深部体温の実測値の時間変化を示す実測値波形との相関係数が、最も１に近いときの値により決定される。或いは、この固定値は、この相関係数が最も１に近い時の値と、ユーザの身体上の特徴に基づく値、例えば、ユーザの体重などに基づく値との相関式により算出することも可能である。
　ここで、産熱量を求める手法について説明する。産熱量は、基本的に、運動時におけるユーザの加速度から求めることができる。ただし、産熱量は、ユーザの加速度だけではなく、ユーザの高さ方向における移動量も考慮した方が、正確に求めることが期待できる。そこで、本実施形態では、産熱量を、ユーザの加速度とユーザの高さ方向における移動量とから求める手法について説明する。この手法に関しては、特開２００８－２２０５１７号公報に詳細に説明されているため、ここでは、要約して説明する。産熱量は、例えば、以下に示す式（７）から式（１０）の４つの式から求められる。
　ＶＭ＝Σ√（ｘ２＋ｙ２＋ｚ２）　式（７）
　ＶＯ２＝０．０４４×ＶＭ＋１．３６５×Ｈｕ＋０．５５３×Ｈｄ　式（８）
　ＥＥ＝ＶＯ２×４．８２５×ＢＷ／１０００　式（９）
　Ｆ０＝ＥＥ×０．８　式（１０）
　式（７）において、ＶＭは、ユーザの加速度の大きさを示す。厳密には、ＶＭは、単位時間に亘ってユーザに課された加速度の累積値である累積加速度の大きさを示す。ＶＭの単位は、ｍＧである。ｘは、ｘ軸方向における加速度の瞬時値である。ｙは、ｙ軸方向における加速度の瞬時値である。ｚは、ｚ軸方向における加速度の瞬時値である。ＶＭは、加速度センサ２２０から処理装置１００に供給される加速度情報に基づいて算出される。具体的には、ＶＭは、予め定められたサンプリング周期（例えば、０．１秒）毎に取得される加速度情報により示される加速度の瞬時値を、単位時間に亘って累積することにより求められる。以下、サンプリング周期は、０．１秒を意味するものとする。なお、重力加速度は、酸素消費量や産熱量に影響を与えないと考えられる。従って、上述した加速度の瞬時値は、重力加速度を除外した値であるものとする。例えば、ｚ軸方向が鉛直方向である場合、ｚ軸方向における加速度の瞬時値から重力加速度を除外した値が、式（７）においてｚとして採用される。
　式（８）において、ＶＯ２は、単位時間当たり、単位体重（例えば、１ｋｇ）当たりの酸素消費量である。ＶＯ２の単位は、ｍＬ／ｋｇ／ｍｉｎである。Ｈｕは、単位時間当たりの累積上昇量である。Ｈｄは、単位時間当たりの累積下降量である。Ｈｕ、Ｈｄの単位は、ｍ／ｍｉｎである。Ｈｕ、Ｈｄは、予め定められたサンプリング周期毎に取得される気圧情報により示される気圧の変化量から算出可能である。なお、ユーザが高さ方向に移動した場合、気圧が変化することを利用する。つまり、ユーザが上方向に移動した場合、気圧が低下し、ユーザが下方向に移動した場合、気圧が上昇することを利用する。
　式（９）において、ＥＥは、総エネルギー消費量を示す。ＥＥの単位は、ｋｃａｌ／ｍｉｎである。４．８２５は、予め定められた定数であり、単位は、ｋｃａｌ／Ｌである。ＢＷは、体重である。ＢＷの単位は、ｋｇである。１０００は、予め定められた定数であり、単位は、ｍＬ／Ｌである。
　式（１０）において、Ｆ０は、単位時間当たりの産熱量である。Ｆ０の単位は、ｋｃａｌ／ｍｉｎである。０．８は、予め定められた定数である。０．８は、総消費エネルギーのうち、運動による消費分が２０％であり、熱に変換される分が８０％であることを示している。
　以上説明したように、ユーザの産熱量は、ユーザに加わる加速度とユーザの上下方向の移動量とから算出可能である。つまり、ユーザの産熱量は、加速度情報と気圧情報とから求めることができる。なお、ユーザの上下方向の移動量が小さい場合、ユーザの産熱量は、ユーザに加わる加速度から算出可能であり、気圧情報は不要である。
　また、角速度センサ２４０により供給される角速度情報に基づいて、ユーザの動きを正確に特定し、ユーザの産熱量の算出精度を向上させることができる。また、磁気センサ２５０により供給される方角情報に基づいて、ユーザの動きを正確に特定し、ユーザの産熱量の算出精度を向上させることができる。
　次に、図５を参照して、本実施形態に係る処理装置１００が実行する深部体温推定処理について説明する。深部体温推定処理は、例えば、処理装置１００の電源が投入されたことに応答して実行される。
　まず、ＣＰＵ１１は、衣服内温度の初期値を測定する（ステップＳ１０１）。例えば、ＣＰＵ１１は、温度センサ２１０から温度情報を取得し、取得した温度情報を衣服内温度の初期値を示す情報としてフラッシュメモリ１４に記憶する。
　ＣＰＵ１１は、ステップＳ１０１の処理を完了すると、加速度を測定する（ステップＳ１０２）。例えば、ＣＰＵ１１は、加速度センサ２２０から加速度情報を取得し、フラッシュメモリ１４に記憶する。
　ＣＰＵ１１は、ステップＳ１０２の処理を完了すると、気圧を測定する（ステップＳ１０３）。例えば、ＣＰＵ１１は、気圧センサ２３０から気圧情報を取得し、フラッシュメモリ１４に記憶する。
　ＣＰＵ１１は、ステップＳ１０３の処理を完了すると、深部体温の推定タイミングであるか否かを判別する（ステップＳ１０４）。深部体温の推定タイミングは、単位時間毎に到来するタイミングであり、例えば、１分毎に到来するタイミングである。ＣＰＵ１１は、ＲＴＣ１５から供給される情報に基づいて、深部体温の推定タイミングであるか否かを判別する。ＣＰＵ１１は、深部体温の推定タイミングでないと判別すると（ステップＳ１０４：ＮＯ）、ステップＳ１０２に処理を戻す。なお、ステップＳ１０２の処理とステップＳ１０３の処理とは、サンプリング周期毎、例えば、０．１秒毎に実行されるものとする。
　ＣＰＵ１１は、深部体温の推定タイミングであると判別すると（ステップＳ１０４：ＹＥＳ）、衣服内温度を測定する（ステップＳ１０５）。例えば、ＣＰＵ１１は、温度センサ２１０から温度情報を取得し、取得した温度情報を衣服内温度を示す情報としてフラッシュメモリ１４に記憶する。
　ＣＰＵ１１は、ステップＳ１０５の処理を完了すると、累積加速度を算出する（ステップＳ１０６）。例えば、ＣＰＵ１１は、式（７）を用いて、ステップＳ１０２で測定された加速度から累積加速度を算出する。
　ＣＰＵ１１は、ステップＳ１０６の処理を完了すると、高さ方向の移動量を算出する（ステップＳ１０７）。例えば、ＣＰＵ１１は、ステップＳ１０３で測定された気圧から高さ方向の移動量を算出する。高さ方向の移動量は、累積上昇量と累積下降量とである。
　ＣＰＵ１１は、ステップＳ１０７の処理を完了すると、産熱量を算出する（ステップＳ１０８）。例えば、ＣＰＵ１１は、ステップＳ１０６で算出した累積加速度と、式（８）と、式（９）と、式（１０）とを用いて、産熱量を算出する。
　ＣＰＵ１１は、ステップＳ１０８の処理を完了すると、深部体温を推定する（ステップＳ１０９）。例えば、ＣＰＵ１１は、ステップＳ１０５で測定した衣服内温度と、ステップＳ１０８で算出した産熱量と、式（１）から式（６）までの連立微分方程式とを用いて、深部体温を推定する。なお、連立微分方程式から算出されるのは、深部体温の変化量であるが、予め定められた深部体温の基準値（例えば、３７度）に深部体温の変化量を加えることにより、深部体温を算出することができる。なお、深部体温を推定することは、基本的に、深部体温の変化量を推定することと同義である。
　ＣＰＵ１１は、ステップＳ１０９の処理を完了すると、深部体温が異常であるか否かを判別する（ステップＳ１１０）。例えば、ＣＰＵ１１は、深部体温が閾値を超えた場合（例えば、３９度）を超えた場合や、深部体温が予め定められた値を超えて上昇した（例えば、１．５度）場合、深部体温が異常であると判別する。ＣＰＵ１１は、深部体温が異常でないと判別すると（ステップＳ１１０：ＮＯ）、ステップＳ１０２に処理を戻す。
　一方、ＣＰＵ１１は、深部体温が異常であると判別すると（ステップＳ１１０：ＹＥＳ）、警報を発する（ステップＳ１１１）。例えば、ＣＰＵ１１は、ブザー３００から警報を応じた音声を出力させる。ＣＰＵ１１は、ステップＳ１１１の処理を完了すると、深部体温推定処理を終了する。
　本実施形態では、ユーザの衣服内の温度である衣服内温度とユーザの体内で産出される熱量である産熱量との対応関係を示す連立微分方程式が用いられて、衣服内温度と産熱量とから深部体温が推定される。従って、本実施形態によれば、ユーザの深部体温を精度良く推定することができる。
　また、本実施形態によれば、深部体温を精度よく推定することができる。このため、推定された深部体温は、熱中症をはじめとする、人体が感じるヒートストレスの危険性の指標としても使用することができる。
　また、本実施形態では、ユーザの血液の温度が、ユーザの深部体温と見做されて推定される。ここで、血液の温度は、例えば、直腸の温度よりも、深部体温に近いと考えられる。このため、本実施形態によれば、ユーザの深部体温を精度良く推定することができる。また、本実施形態では、測定対象の温度は、衣服内温度である。従って、例えば、直腸の温度を測定する場合に比べ、温度の測定が容易である。
　本実施形態では、６つの式により構成される連立微分方程式により定義される数理モデルが、産熱量と衣服内温度とから深部体温を精度良く求めることができる前提で説明した。以下、実施例１から実施例３では、この数理モデルの妥当性の検証結果や、この数理モデルで用いるパラメータの決定方法について説明する。
（実施例１）
　まず、実施例１における条件について説明する。実施例１では、衣服５００を着用した被験者に、深部体温としての食道温を実測するためのプローブ式深部体温測定装置（図示せず）と、衣服内温度を測定する温度センサ２１０と、加速度センサ２２０と、気圧センサ２３０と、角速度センサ２４０と、磁気センサ２５０とを装着させた。また、カプセル式発汗量測定装置、マスク型酸素摂取量測定装置なども使用した。被験者は、Ｓ、Ｈ、Ｍ、Ｙ、Ｋ、Ｕにより示される６人とした。プローブ式深部体温測定装置のプローブは、被験者の食道に装着された。温度センサ２１０と加速度センサ２２０と角速度センサ２４０と磁気センサ２５０とは、衣服５００が備える内ポケット５１１内に配置された。気圧センサ２３０は、衣服５００が備える外ポケット５２１内に配置された。
　被験者が着用する衣服５００は、特開２０１４－０９１３０７号公報の比較例４に記載された衣服であるものとする。この衣服５００は、多層構造布帛であり、消防用防護服（防火服）の形状に縫合されたものである。具体的には、最外層には、ポリメタフェニレンイソフタルアミド繊維（帝人社製、商標名：コーネックス）と、コポリパラフェニレン・３、４’－オキシジフェニレンテレフタルアミド繊維（帝人テクノプロダクツ社製、商標名：テクノーラ）とを、混合比率が９０：１０となる割合で混合した耐熱繊維からなる紡績糸（番手：４０／２）を用いて、平織リップストップに構成した布帛を製織した。この表地層の目付けは、３８０ｇ／ｍ２であった。
　中間層には、ポリメタフェニレンイソフタルアミド繊維（帝人社製、商標名：コーネックス）と、コポリパラフェニレン・３、４’－オキシジフェニレンテレフタルアミド繊維（帝人テクノプロダクツ社製、商標名：テクノーラ）とを、混合比率が９５：５となる割合で混合した耐熱繊維からなる紡績糸（番手：４０／－）を用いて平織に織成した織布（目付：８０ｇ／ｍ２）に、ポリテトラフルオロエチレン製の透湿防水性フィルム（ジャパンゴアテックス社製）をラミネートしたものを採用した。遮熱層（内層）には、コポリパラフェニレン・３、４’－オキシジフェニレンテレフタルアミド繊維ヤーン（帝人テクノプロダクツ社製、商標名：テクノーラ）からなる総繊度１６７０ｄｔｅｘのフィラメントを用いて、平織に構成した布帛を製織した。遮熱層の目付は２１０ｇ／ｍ２であった。
　被験者を歩行させるトレッドミルの傾斜角度は５％に設定され、トレッドミルの運転速度を、段階的に早くした。具体的には、運転速度は、測定開始からの経過時間が０分から３分までの３分間は０ｋｍ／ｈとし、経過時間が３分から１２分までの９分間は１．５ｋｍ／ｈとし、経過時間が１２分から２１分までの９分間は３．０ｋｍ／ｈとし、経過時間が２１分から３０分までの９分間は４．５ｋｍ／ｈとし、経過時間が３０分以後は６．０ｋｍ／ｈとした。なお、安全のため、心拍数１８０以上、あるいは、食道温３８．５度に達した時点で評価を終了した。
　図６に、被験者の基本データを示す。基本データは、年齢、性別、身長（ｃｍ）、体重（ｋｇ）、ＢＭＩ、体表面積（ｃｍ２）、測定時間（ｍｉｎ）、運動後体重（ｋｇ）、体重差（ｋｇ）を含む。ＢＭＩは、体重／身長２により示される。体表面積は、体重０．４４４×身長０．６６３×８８．８３／１００００により示される。体重差は、体重（運動前体重）から運動後体重を減じた値である。
　産熱量は、実施形態で説明した手法により算出した。ただし、産熱量は、特開２００８－２２０５１７号公報に記載された手法や、特許第３５７１２７２号公報に記載された手法を採用してもよい。衣服内温度は、温度センサ２１０により測定された温度である。図７に、実施例１の結果データを示す。図８に、実測値と算出値との関係を示す。なお、実施例１は、全パラメータを個別に算出する実施例である。パラメータは、式（１）から式（６）までの６つの式を含む連立微分方程式における、Ｃ１、Ｃ２、Ｒ１、Ｒ２の４つである。
　図８に示すように、深部体温の実測値は、時間の経過とともに上昇する。なお、理解を容易にするため、深部体温は、測定開始時の深部体温を基準とした温度変化量（深部体温上昇度）により示している。ここで、上述した連立微分方程式により表される数理モデルが適切である場合、数理モデルを用いて算出された深部体温の時間変化と、実測された深部体温の時間変化との差が小さいはずである。つまり、深部体温の算出値の時間変化を示す算出値波形と、深部体温の実測値の時間変化を示す実測値波形との相関係数が、１に近いはずである。そこで、被験者毎に、Ｃ１、Ｃ２、Ｒ１、Ｒ２のそれぞれを切り替えながら、深部体温の算出値を求めて算出値波形を取得する。そして、算出値波形と実測値波形との相関係数が、最も１に近いときのパラメータ（Ｃ１、Ｃ２、Ｒ１、Ｒ２）を特定した。図７は、被験者毎に相関係数が最大となるパラメータ（Ｃ１、Ｃ２、Ｒ１、Ｒ２）を示している。
　なお、相関係数が１から０．７の間であることは、強い正の相関があることを意味する。相関係数が０．７から０．４の間であることは、中程度の正の相関があることを意味する。相関係数が０．４から０．２の間であることは、弱い正の相関があることを意味する。相関係数が０．２から－０．２の間であることは、ほとんど相関がないことを意味する。相関係数が－０．２から－０．４の間であることは、弱い負の相関があることを意味する。相関係数が－０．４から－０．７の間であることは、中程度の負の相関があることを意味する。相関係数が－０．７から－１の間であることは、強い負の相関があることを意味する。
　Ｓｌｏｐｅは、算出値を横軸、実測値を縦軸として、１分毎の深部体温の変化量をプロットした場合において、プロットされた点により近似される直線の傾きである。なお、この直線は、例えば、最小二乗法により求めることができる。Ｉｎｔｅｒｃｅｐｔは、この直線の切片である。ＡＶＧは、平均値を示す。ＳＤは、標準偏差、つまり、ばらつきを示す。ＳＥは、標準誤差、つまり、推定量の推定誤差を示す。
　図７に示すように、被験者毎に、Ｃ１、Ｃ２、Ｒ１、Ｒ２を調整した場合、どの被験者についても、相関係数がほぼ０．９７以上となった。つまり、実施例１によれば、上述した数理モデルが尤もらしいと考えられる。
（実施例２）
　実施例２では、実施例１の結果データを踏まえて、パラメータの一般化を試みた。具体的には、実施例２では、実施例１で求められた被験者毎のパラメータの平均値を採用した。図９に、実施例２の結果データを示す。
　図９に示すように、全パラメータ（Ｃ１、Ｃ２、Ｒ１、Ｒ２）に平均値を採用した場合、どの被験者についても、相関係数がほぼ０．９７以上となった。つまり、実施例２によれば、上述した数理モデルにおいて、全パラメータ（Ｃ１、Ｃ２、Ｒ１、Ｒ２）に予め定められた値（つまり、平均値）を採用しても、数理モデルが成立すると考えられる。
（実施例３）
　実施例３では、実施例１の結果データを踏まえて、パラメータがユーザの体重と相関関係があるか検証した。具体的には、実施例３では、まず、全パラメータ（Ｃ１、Ｃ２、Ｒ１、Ｒ２）について、体重との間に相関があるか否かを検証した。図１０に、実施例３の結果データ１を示す。
　図１０において、相関係数は、各パラメータと体重との相関係数である。図１１を参照して、各パラメータと体重との相関係数を求める手法について説明する。図１１は、体重とＣ１との関係を示す図であり、体重を横軸、Ｃ１を縦軸にしたときの各被験者の散布図である。図１１において、ＰＳは、Ｓにより示される被験者（以下、単に「Ｓ」という。他の被験者についても同様。）の体重と、Ｓの実施例１による調整後のＣ１とにより特定される点を示す。ＰＨは、Ｈの体重と、Ｈの調整後のＣ１とにより特定される点を示す。ＰＭは、Ｍの体重と、Ｍの調整後のＣ１とにより特定される点を示す。ＰＹは、Ｙの体重と、Ｙの調整後のＣ１とにより特定される点を示す。ＰＫは、Ｋの体重と、Ｋの調整後のＣ１とにより特定される点を示す。ＰＵは、Ｕの体重と、Ｕの調整後のＣ１とにより特定される点を示す。
　ここで、体重とＣ１との相関係数であるｒは、全被験者の体重を順に並べたときの体重波形と、全被験者のＣ１を順に並べたときのＣ１波形との相関係数である。ここで、ｒ＝０．７５４、ｒ２＝０．５６８となった。体重とＣ１との相関係数は、０．７５４であることから、体重とＣ１とは正の相関が強いと考えられる。つまり、体重を指標としてＣ１を設定することは、有効であると考えられる。そこで、体重とＣ１との対応関係を特定することを試みる。まず、Ｌ１は、ＰＳ、ＰＨ、ＰＭ、ＰＹ、ＰＫ、ＰＵを近似する直線を示す。Ｌ１は、例えば、最小二乗法により特定される。ここで、体重をＢＷとすると、Ｃ１＝０．１２３４×ＢＷ＋１３．０６２となる。
　体重とＣ２との関係、体重とＲ１との関係、体重とＲ２との関係についても、体重とＣ１との関係と同様に、相関係数を求めることにより検証することができる。体重とＣ２との相関係数は０．５５５であることから、体重とＣ２とは正の相関がある程度あると考えられる。そこで、体重とＣ２との対応関係を求めると、Ｃ２＝０．０８８７×ＢＷ＋１５．９６３となる。一方、体重とＲ１との相関係数は０．３５６であることから、体重とＲ１とはあまり相関がないと考えられる。同様に、体重とＲ２との相関係数は０．０３９であることから、体重とＲ２とはあまり相関がないと考えられる。
　そこで、実施例３では、Ｃ１とＣ２とは、体重を用いた相関式で求め、Ｒ１とＲ２とは、平均値を採用する。図１２に、実施例３の結果データ２を示す。図１２に示すように、Ｃ１とＣ２とを体重を用いて特定し、Ｒ１とＲ２とに平均値を採用すると、どの被験者についても、相関係数がほぼ０．９７以上となった。つまり、実施例３によれば、上述した数理モデルにおいて、Ｃ１とＣ２とを体重を用いて特定し、Ｒ１とＲ２とに予め定められた値（つまり、平均値）を採用しても、数理モデルが成立すると考えられる。
（実施例４）
　実施例４では、実施例１の結果データを踏まえて、パラメータがユーザの体表面積と相関関係があるか検証した。具体的には、実施例４では、まず、全パラメータ（Ｃ１、Ｃ２、Ｒ１、Ｒ２）について、体表面積との間に相関係数があるか否かを検証した。図１３に、実施例４の結果データ１を示す。
　図１３において、相関係数は、各パラメータと体表面積との相関係数である。各パラメータと体表面積との相関係数を求める手法は、基本的に、各パラメータと体重との相関係数を求める手法と同様である。図１３に示す例では、Ｃ１に関しては、ｒ＝０．８４４、ｒ２＝０．７１２である。体表面積とＣ１との相関係数は、０．８４４であることから、体表面積とＣ１とは正の相関が強いと考えられる。つまり、体表面積を指標としてＣ１を設定することは、有効であると考えられる。そこで、体表面積とＣ１との対応関係を特定することを試みる。体表面積をＢＳＡとすると、Ｃ１＝０．００１０×ＢＳＡ＋０．７４５０となる。
　体表面積とＣ２との関係、体表面積とＲ１との関係、体表面積とＲ２との関係についても、体表面積とＣ１との関係と同様に、相関係数を求めることにより検証することができる。体表面積とＣ２との相関係数は０．６８０であることから、体表面積とＣ２とは正の相関がある程度あると考えられる。そこで、体表面積とＣ２との対応関係を求めると、Ｃ２＝０．０００９×ＢＳＡ＋５．７３６０となる。一方、体表面積とＲ１との相関係数は０．４２９であることから、体表面積とＲ１とはあまり相関がないと考えられる。同様に、体表面積とＲ２との相関係数は０．１２５であることから、体表面積とＲ２とはあまり相関がないと考えられる。
　そこで、実施例４では、Ｃ１とＣ２とは、ＢＳＡを用いた式で求め、Ｒ１とＲ２とは、平均値を採用する。図１４に、実施例４の結果データ２を示す。図１４に示すように、Ｃ１とＣ２とをＢＳＡを用いて特定し、Ｒ１とＲ２とに平均値を採用すると、どの被験者についても、相関係数がほぼ０．９７以上となった。つまり、実施例４によれば、上述した数理モデルにおいて、Ｃ１とＣ２とをＢＳＡを用いて特定し、Ｒ１とＲ２とに予め定められた値（つまり、平均値）を採用しても、数理モデルが成立すると考えられる。
（実施例５）
　実施例５では、実施例３及び実施例４の結果データを踏まえて、Ｃ１を体重により特定し、Ｃ２を体表面積により特定し、Ｒ１及びＲ２として平均値を採用した。つまり、実施例５では、Ｃ１＝０．１２３４×ＢＷ＋１３．０６２、Ｃ２＝０．０００９×ＢＳＡ＋５．７３６０、Ｒ１＝２．４、Ｒ２＝２４．２を採用した。図１５に、実施例５の結果データを示す。
　図１５において、Ｔ値は、統計検定量である。Ｔ値は、Ｔ＝ｒ×（√（Ｎ－２））／（√１－ｒ２）により表現される。なお、Ｎは、データ数である。また、Ｐ値は、Ｔ分布の両側確率である。図１５に示すように、Ｃ１をＢＷを用いて特定し、Ｃ２をＢＳＡを用いて特定し、Ｒ１とＲ２とに平均値を採用すると、どの被験者についても、相関係数がほぼ０．９７以上となった。つまり、実施例５によれば、上述した数理モデルにおいて、Ｃ１をＢＷを用いて特定し、Ｃ２をＢＳＡを用いて特定し、Ｒ１とＲ２とに予め定められた値（つまり、平均値）を採用しても、数理モデルが成立すると考えられる。
　図１６に、実施例１から実施例５に基づく、最終結果データを示す。図１６には、４つの条件において、相関係数と、ＳＬＯＰＥと、Ｉｎｔｅｒｃｅｐｔとを示している。４つの条件は、体重によりＣ１を特定し体表面積によりＣ２を特定する第１の条件と、体重によりＣ１を特定しＣ２として平均値を採用する第２の条件と、体重によりＣ１とＣ２とを特定する第３の条件と、体表面積によりＣ１とＣ２とを特定する第４の条件と、を含む。なお、いずれの条件においても、Ｒ１とＲ２とに予め定められた値（つまり、平均値）を採用している。また、第１の条件は実施例５に対応し、第３の条件は実施例３に対応し、第４の条件は実施例４に対応する。
　図１６に示すように、Ｃ１とＣ２とをユーザの身体上の特徴に基づく値（例えば、体重や体表面積）により特定し、Ｒ１とＲ２とに予め定められた値（つまり、平均値）を採用すれば、数理モデルが成立すると考えられる。なお、相関係数と、ＳＬＯＰＥと、Ｉｎｔｅｒｃｅｐｔとを考慮すると、第１の条件（実施例５）のように、体重によりＣ１を特定し、体表面積によりＣ２を特定し、Ｒ１とＲ２とに予め定められた値（つまり、平均値）を採用すると、最も好ましいと考えられる。
（変形例）
　以上、本発明の実施形態を説明したが、本発明を実施するにあたっては、種々の形態による変形及び応用が可能である。
　本発明において、上記実施形態において説明した構成、機能、動作のどの部分を採用するのかは任意である。また、本発明において、上述した構成、機能、動作のほか、更なる構成、機能、動作が採用されてもよい。また、上述した実施形態は、適宜、自由に組み合わせることができる。また、本発明において採用可能な素材、サイズ、電気的特性などが、上記実施形態において示したものに限定されないことは勿論である。
　上記衣服内温度は、衣服と人体との間の空気温度、最外層に着用している衣服とその内側の衣服との間の空気温度、衣服に覆われている人体の皮膚温度、などが挙げられるがこれに限定されるものではない。
　上記ユーザの身体上の特徴とは体重や身長、ＢＭＩ、体脂肪率、発汗量、呼吸量、体表面積、筋肉量などが挙げられるが、これに限定されない。
　上記報知手段は、例えば、上記深部体温が予め定められた閾値を超える場合、あるいは上記深部体温の上昇度合いが予め定められた上昇限界度合いを超える場合に報知するが、これに限定されない。
　実施形態１では、処理装置１００が有する機能を多くがソフトウェア（又は、ファームウェア）により実現される例、つまり、処理装置１００が有する機能を多くがプロセッサによるプログラムの実行により実現されるについて説明した。本発明において、このような機能は、ハードウェアにより実現されてもよい。この場合、例えば、処理装置１００は、ＣＰＵ１１に代えて、処理回路を備える。この処理回路は、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ（Application Specific Integrated Circuit）、ＦＰＧＡ（Field Programmable Gate Array）、又はこれらの組合せにより構成される。
　本発明に係る処理装置１００の動作を規定する動作プログラムを既存のパーソナルコンピュータや情報端末装置に適用することで、当該パーソナルコンピュータを本発明に係る処理装置１００として機能させることも可能である。また、このようなプログラムの配布方法は任意であり、例えば、ＣＤ－ＲＯＭ（Compact Disk Read-Only Memory）、ＤＶＤ（Digital Versatile Disk）、メモリカードなどのコンピュータ読み取り可能な記録媒体に格納して配布してもよいし、通信ネットワーク（例えば、インターネット）を介して配布してもよい。
　本発明は、本発明の広義の精神と範囲を逸脱することなく、様々な実施形態及び変形が可能とされるものである。また、上述した実施形態は、この発明を説明するためのものであり、本発明の範囲を限定するものではない。すなわち、本発明の範囲は、実施形態ではなく、請求の範囲によって示される。そして、請求の範囲内及びそれと同等の発明の意義の範囲内で施される様々な変形が、この発明の範囲内とみなされる。
　本出願は、２０１７年７月３１日に出願された日本国特許出願特願２０１７－１４７３５５号に基づく。本明細書中に日本国特許出願特願２０１７－１４７３５５号の明細書、特許請求の範囲、図面全体を参照として取り込むものとする。
　１１　ＣＰＵ、１２　ＲＯＭ、１３　ＲＡＭ、１４　フラッシュメモリ、１５　ＲＴＣ、１６　タッチスクリーン、１７　第１の通信インターフェース、１８　第２の通信インターフェース、１００　処理装置、２１０　温度センサ、２２０　加速度センサ、２３０　気圧センサ、２４０　角速度センサ、２５０　磁気センサ、３００　ブザー、４００　端末装置、５００　衣服、５１０　裏地、５１１　内ポケット、５２０　表地、５２１　外ポケット、５３０，５４０，５５０　空間、６００　下着、１０００　深部体温推定システム。

Claims (16)

  1. 　ユーザが着用する衣服の内側の温度である衣服内温度を測定する温度センサと、
    　前記ユーザに加わる加速度を検出する加速度センサと、を備え、
    　前記ユーザの体内の温度である深部体温を推定する、
    　深部体温推定システム。
  2. 　前記温度センサにより測定された衣服内温度と前記加速度センサにより検出された加速度とに基づいて、前記深部体温を推定する深部体温推定手段を更に備える、
    　請求項１に記載の深部体温推定システム。
  3. 　前記加速度センサにより検出された加速度に基づいて、前記ユーザの体内で産出される熱量である産熱量を推定する産熱量推定手段を更に備え、
    　前記深部体温推定手段は、前記温度センサにより測定された衣服内温度と前記産熱量推定手段により推定された産熱量とに基づいて、前記深部体温を推定する、
    　請求項２に記載の深部体温推定システム。
  4. 　前記深部体温推定手段は、前記衣服内温度と前記産熱量と前記深部体温との対応関係を示す数理モデルを用いて、前記深部体温を推定する、
    　請求項３に記載の深部体温推定システム。
  5. 　前記数理モデルは、
    　前記ユーザの筋肉が有する熱量である第１の熱量の変化量を、前記産熱量と前記ユーザの筋肉から前記ユーザの血液に供給される熱流量である第１の熱流量との差で示す第１の式と、
    　前記ユーザの血液が有する熱量である第２の熱量の変化量を、前記第１の熱流量と前記ユーザの血液から前記衣服の内側に供給される熱流量である第２の熱流量との差で示す第２の式と、
    　前記第１の熱流量を、前記ユーザの筋肉から前記ユーザの血液への熱抵抗値である第１の熱抵抗値に対する、前記ユーザの筋肉の温度である第１の温度の変化量と前記ユーザの血液の温度である第２の温度の変化量との差の割合で示す第３の式と、
    　前記第２の熱流量を、前記ユーザの血液から前記衣服の内側への熱抵抗値である第２の熱抵抗値に対する、前記第２の温度の変化量と前記衣服内温度の変化量との差の割合で示す第４の式と、
    　前記第１の温度の変化量を、前記ユーザの筋肉が有する熱容量である第１の熱容量に対する前記第１の熱量の割合で示す第５の式と、
    　前記第２の温度の変化量を、前記ユーザの血液が有する熱容量である第２の熱容量に対する前記第２の熱量の割合で示す第６の式と、を含む連立微分方程式により表され、
    　前記深部体温推定手段は、前記第２の温度を前記深部体温として推定する、
    　請求項４に記載の深部体温推定システム。
  6. 　前記第１の熱容量と前記第２の熱容量とは、前記ユーザの身体上の特徴に基づく値であり、
    　前記第１の熱抵抗値と前記第２の熱抵抗値とは、予め定められた値である、
    　請求項５に記載の深部体温推定システム。
  7. 　前記ユーザの周囲の気圧を検出する気圧センサを更に備え、
    　前記産熱量推定手段は、前記加速度センサにより検出された加速度と前記気圧センサにより検出された気圧とに基づいて、前記産熱量を推定する、
    　請求項３から６の何れか１項に記載の深部体温推定システム。
  8. 　前記ユーザに加わる角速度を検出する角速度センサと、
    　前記ユーザの周囲の磁場の方向を検出する磁気センサと、を更に備え、
    　前記加速度センサと前記角速度センサと前記磁気センサとは、相互に固定されており、
    　前記産熱量推定手段は、前記加速度センサにより検出された加速度と前記角速度センサにより検出された角速度と前記磁気センサにより検出された方向とに基づいて、前記産熱量を推定する、
    　請求項３から７の何れか１項に記載の深部体温推定システム。
  9. 　前記深部体温推定手段により推定された深部体温に基づいて、前記ユーザの異常を報知する報知手段を更に備える、
    　請求項２から８の何れか１項に記載の深部体温推定システム。
  10. 　前記衣服が、遮熱効果を有する衣服である、
    　請求項１から９の何れか１項に記載の深部体温推定システム。
  11. 　前記遮熱効果が、ＩＳＯ９１５１に準拠した方法で測定されるＨＴＩ２４が３秒以上であることを満たす、
    　請求項１０に記載の深部体温推定システム。
  12. 　前記衣服が、ＪＩＳ　Ｌ　１０９６に準拠した方法であってフラジール形試験機を用いた方法で測定される通気度の値が１．０ｃｍ３／ｃｍ２／ｓ以下である、
    　請求項１から１１の何れか１項に記載の深部体温推定システム。
  13. 　前記衣服が、ＪＩＳ　Ｌ　１０９９に準拠した方法であってカップ法により測定される透湿度の値が１０００ｇ／ｍ２／ｈ以下である、
    　請求項１から１２の何れか１項に記載の深部体温推定システム。
  14. 　前記衣服の熱伝導率が、１ｋｃａｌ／ｈ／ｍ／℃以下である、
    　請求項１から１３の何れか１項に記載の深部体温推定システム。
  15. 　ユーザが着用する衣服の内側の温度である衣服内温度と前記ユーザに加わる加速度とに基づいて、前記ユーザの体内の温度である深部体温を推定する深部体温推定手段と、
    　前記深部体温推定手段により推定された深部体温に基づいて、警報を発する警報手段と、を備える、
    　ヒートストレス警報システム。
  16. 　ユーザが着用する衣服の内側の温度である衣服内温度を測定し、
    　前記ユーザに加わる加速度を検出し、
    　前記測定された衣服内温度と前記検出された加速度とに基づいて、前記ユーザの体内の温度である深部体温を推定する、
    　深部体温推定方法。
PCT/JP2018/022899 2017-07-31 2018-06-15 深部体温推定システム、ヒートストレス警報システム、及び、深部体温推定方法 WO2019026439A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2019533946A JP6826667B2 (ja) 2017-07-31 2018-06-15 深部体温推定システム、ヒートストレス警報システム、及び、深部体温推定方法
EP18841143.3A EP3662819B1 (en) 2017-07-31 2018-06-15 Deep body temperature estimation system, and deep body temperature estimation method
US16/635,390 US11291373B2 (en) 2017-07-31 2018-06-15 Deep-body-temperature estimation system, heat stress warning system, and deep-body-temperature estimation method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2017-147355 2017-07-31
JP2017147355 2017-07-31

Publications (1)

Publication Number Publication Date
WO2019026439A1 true WO2019026439A1 (ja) 2019-02-07

Family

ID=65233680

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2018/022899 WO2019026439A1 (ja) 2017-07-31 2018-06-15 深部体温推定システム、ヒートストレス警報システム、及び、深部体温推定方法

Country Status (4)

Country Link
US (1) US11291373B2 (ja)
EP (1) EP3662819B1 (ja)
JP (1) JP6826667B2 (ja)
WO (1) WO2019026439A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2021040872A (ja) * 2019-09-10 2021-03-18 株式会社テクノ・コモンズ 生体信号処理方法およびその装置
JP2021133184A (ja) * 2020-02-28 2021-09-13 公立大学法人公立諏訪東京理科大学 熱中症予防システム
GB2593433A (en) * 2020-02-10 2021-09-29 Prevayl Ltd Electronics arrangement for a wearable article
WO2022013911A1 (ja) * 2020-07-13 2022-01-20 日本電信電話株式会社 温度測定装置、方法およびプログラム

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2584492B (en) 2019-06-07 2021-08-18 Prevayl Ltd Method, garment and system

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS571272B2 (ja) 1976-02-13 1982-01-09
JP2004030180A (ja) 2002-06-25 2004-01-29 Masami Murayama 防護装具のヒートストレス報知システム
JP2008220517A (ja) 2007-03-09 2008-09-25 Npo Jukunen Taiiku Daigaku Research Center 消費カロリー算出方法および携帯用消費カロリー測定装置
JP2010255124A (ja) 2009-04-21 2010-11-11 Teijin Techno Products Ltd 耐熱性防護服
JP2011106069A (ja) 2009-11-19 2011-06-02 Teijin Techno Products Ltd 軽量耐熱防護服
JP2014091307A (ja) 2012-11-06 2014-05-19 Teijin Ltd 遮熱活動服用積層布帛およびそれを用いた遮熱活動服
JP2017147355A (ja) 2016-02-18 2017-08-24 本田技研工業株式会社 電装部品防水構造

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8830068B2 (en) 2010-09-15 2014-09-09 Colorado State University Research Foundation Multi-sensor environmental and physiological monitor system and methods of use
US9823138B2 (en) 2013-08-21 2017-11-21 Petpace Ltd Non-invasive automatic monitoring of pet animal's core temperature
CN109074716A (zh) * 2016-04-19 2018-12-21 帝人株式会社 具备警报系统的物品
JP6675552B2 (ja) * 2016-06-08 2020-04-01 国立大学法人大阪大学 深部体温推定装置、その方法及びプログラム

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS571272B2 (ja) 1976-02-13 1982-01-09
JP2004030180A (ja) 2002-06-25 2004-01-29 Masami Murayama 防護装具のヒートストレス報知システム
JP2008220517A (ja) 2007-03-09 2008-09-25 Npo Jukunen Taiiku Daigaku Research Center 消費カロリー算出方法および携帯用消費カロリー測定装置
JP2010255124A (ja) 2009-04-21 2010-11-11 Teijin Techno Products Ltd 耐熱性防護服
JP2011106069A (ja) 2009-11-19 2011-06-02 Teijin Techno Products Ltd 軽量耐熱防護服
JP2014091307A (ja) 2012-11-06 2014-05-19 Teijin Ltd 遮熱活動服用積層布帛およびそれを用いた遮熱活動服
JP2017147355A (ja) 2016-02-18 2017-08-24 本田技研工業株式会社 電装部品防水構造

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HAMATANI, TAKASHI: "A solar radiation model and parameter calibration in a human thermal model using a wearable sensor (non-official translation)", DICOMO 2015 PROCEEDINGS OF THE SYMPOSIUM OF MULTIMEDIA, DISTRIBUTED, COOPERATIVE, AND MOBILE (IPSJ SYMPOSIUM SERIES) CD-ROM, July 2015 (2015-07-01), pages 381 - 391, XP009519251 *
HAMATANI, TAKASHI: "A study on deep body temperature estimation using wearable sensor (non-official translation)", RESEARCH REPORT OF INFORMATION PROCESSING SOCIETY OF JAPAN, vol. 2014-MBL-072, 27 August 2014 (2014-08-27), pages 1 - 5, XP009519253, ISSN: 0919-6072 *
HAMATI, TAKASHI: "Proposal of deep body temperature estimation method based on body heat model using wearable sensor", IPSJ SIG TECHNICAL REPORTS, ITS (INTELLIGENT TRANSPORT SYSTEMS AND SMART COMMUNITY), vol. 2014-ITS-059, 20 November 2014 (2014-11-20), pages 1 - 5, XP009519257, ISSN: 0919-6072 *
See also references of EP3662819A4

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2021040872A (ja) * 2019-09-10 2021-03-18 株式会社テクノ・コモンズ 生体信号処理方法およびその装置
JP7417919B2 (ja) 2019-09-10 2024-01-19 株式会社テクノ・コモンズ 生体信号処理装置
GB2593433A (en) * 2020-02-10 2021-09-29 Prevayl Ltd Electronics arrangement for a wearable article
GB2593433B (en) * 2020-02-10 2024-01-10 Prevayl Innovations Ltd Electronics arrangement for a wearable article
JP2021133184A (ja) * 2020-02-28 2021-09-13 公立大学法人公立諏訪東京理科大学 熱中症予防システム
JP7417933B2 (ja) 2020-02-28 2024-01-19 公立大学法人公立諏訪東京理科大学 熱中症予防システム
WO2022013911A1 (ja) * 2020-07-13 2022-01-20 日本電信電話株式会社 温度測定装置、方法およびプログラム

Also Published As

Publication number Publication date
JP6826667B2 (ja) 2021-02-03
EP3662819B1 (en) 2022-03-30
US20200367758A1 (en) 2020-11-26
JPWO2019026439A1 (ja) 2020-05-28
US11291373B2 (en) 2022-04-05
EP3662819A1 (en) 2020-06-10
EP3662819A4 (en) 2021-06-30

Similar Documents

Publication Publication Date Title
JP6826667B2 (ja) 深部体温推定システム、ヒートストレス警報システム、及び、深部体温推定方法
US9858794B2 (en) Detecting and notifying of various potential hazards
TWI673035B (zh) 具備警報系統的防護裝備
JP7222158B2 (ja) 頭部装着装置、熱中症予防システム及び水分補給警告システム
EP3447745A1 (en) Article provided with warning system
Caya et al. Warning system for firefighters using e-textile
WO2020184687A1 (ja) 頭部装着装置、熱中症予防システム及び水分補給警告システム
JP7411948B2 (ja) 衣服、熱中症予防システム及び水分補給警告システム
JP7432204B2 (ja) 全身発汗量推定システム及び熱中症予防システム
JP2017192492A (ja) 運動量算出システムを備えた物品および警報システムを備えた物品
Florea et al. Wearable system for heat stress monitoring in firefighting applications
JP2017194781A (ja) 警報システムを備えた物品
Jia et al. Moisture based perspiration level estimation
JP2017187913A (ja) 警報システムを備えた防護装備
Magenes et al. Biosensing and environmental sensing for emergency and protection e-Textiles
Presti et al. Feasibility assessment of an FBG-based soft sensor embedded into a single-use surgical mask for respiratory monitoring
US11672290B2 (en) Protective glove
WO2020184688A1 (ja) 頭部装着装置、熱中症予防システム及び水分補給警告システム
Jo et al. Machine Embroidery Enclosure for Stretchable Fiber Optic Respiration Sensor
JP2021115370A (ja) トレーニング状態評価システム
JP2017187943A (ja) センサーを備えた防護装備
JP2017187942A (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: 18841143

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2019533946

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2018841143

Country of ref document: EP

Effective date: 20200302