JP2021194031A - 深部体温推定方法および装置 - Google Patents

Abstract

【課題】より容易に深部温度が推定できるようにする方法を提供する。【解決手段】第２演算部１０５は、時刻ｔにおいて第１演算部１０４が求めた熱量、第２測定部１０２が測定した温度と第３測定部１０３が測定した湿度とから求められる測定対象者の雰囲気と測定対象者の皮膚との間の第１熱交換量、測定対象者の皮膚と深部との間の第２熱交換量、および測定対象者の発汗量より求める汗の蒸発量から、測定対象者の深部体温の変化を推定する。【選択図】図１

Description

本発明は、人体の体内の中心部分の温度である深部体温を推定する深部体温推定方法および装置に関する。

毎年、熱中症により多くの人が救急搬送されており、２０１９年５月から９月の全国における熱中症による救急搬送人員数の累計は、７１，３１７人であった（非特許文献１参照）。熱中症の診療や、熱中症のリスクを測る方法として、深部体温（体内の温度）の測定が有効である（非特許文献２参照）。

２０１９年（５月から９月）の熱中症による救急搬送状況、総務省、［２０２０年４月２１日検索］、（https://www.fdma.go.jp/disaster/heatstroke/items/heatstroke_geppou_2019.pdf）。 熱中症診療ガイドライン、厚生労働省、［２０２０年４月２１日検索］、（https://www.mhlw.go.jp/file/06-Seisakujouhou-10800000-Iseikyoku/heatstroke2015.pdf）。

しかしながら、深部体温を監視するためには、直腸温、膀胱温、食道温などを測定する必要があり、日常生活のなかでこれらの測定を実施することは難しい。

本発明は、以上のような問題点を解消するためになされたものであり、より容易に深部温度が推定できるようにすることを目的とする。

本発明に係る深部体温推定方法は、測定対象者の心拍数を測定する第１ステップと、測定対象者の近傍の温度を測定する第２ステップと、測定対象者の近傍の湿度を測定する第３ステップと、第１ステップで測定した心拍数より測定対象者の運動により発生する熱量を求める第４ステップと、時刻ｔにおいて第４ステップで求めた熱量、第２ステップで測定した温度と第３ステップで測定した湿度とから求められる測定対象者の雰囲気と測定対象者の皮膚との間の第１熱交換量、測定対象者の皮膚と深部との間の第２熱交換量、および測定対象者の発汗量より求める汗の蒸発量から、測定対象者の深部体温の変化を推定する第５ステップとを備える。

また、本発明に係る深部体温推定装置は、測定対象者の心拍数を測定する第１測定部と、測定対象者の近傍の温度を測定する第２測定部と、測定対象者の近傍の湿度を測定する第３測定部と、第１測定部が測定した心拍数より測定対象者の運動により発生する熱量を求める第１演算部と、時刻ｔにおいて第１演算部が求めた熱量、第２測定部が測定した温度と第３測定部が測定した湿度とから求められる測定対象者の雰囲気と測定対象者の皮膚との間の第１熱交換量、測定対象者の皮膚と深部との間の第２熱交換量、および測定対象者の発汗量より求める汗の蒸発量から、測定対象者の深部体温の変化を推定する第２演算部とを備える。

以上に説明したように、本発明によれば、熱量、測定対象者の雰囲気と測定対象者の皮膚との間の第１熱交換量、測定対象者の皮膚と深部との間の第２熱交換量、および測定対象者の汗の蒸発量から、測定対象者の深部体温の変化を推定するので、より容易に深部温度が推定できる。

図１は、本発明の実施の形態に係る深部体温推定装置の構成を示す構成図である。 図２は、本発明の実施の形態に係る深部体温推定方法を説明するフローチャートである。 図３は、本発明の実施の形態に係る深部体温推定方法を実施するためのウエアラブルデバイスの構成を示す構成図である。 図４は、本発明の実施の形態に係る深部体温推定方法を実施するためのウエアラブルデバイスの構成を示す構成図である。

以下、本発明の実施の形態に係る深部体温推定装置について図１を参照して説明する。この深部体温推定装置は、第１測定部１０１、第２測定部１０２、第３測定部１０３、第１演算部１０４、および第２演算部１０５を備える。

第１測定部１０１は、測定対象者の心拍数を測定する。第１測定部１０１は、例えば、被測定者の心臓の心電位を計測する心電計と、心電計が計測した心電位から心拍数を算出する算出部とから構成することができる。また、第１測定部１０１は、ウエア型の心拍計、ベルト型の心拍計、腕時計型の心拍計、イヤホン型の心拍計などから構成することができる。

第２測定部１０２は、測定対象者の近傍の温度を測定する。第２測定部１０２は、例えば温度計から構成することができる。第３測定部１０３は、測定対象者の近傍の湿度を測定する。第３測定部１０３は、例えば、湿度計から構成することができる。また、第２測定部１０２、第３測定部１０３は、第１測定部１０１と一体化した温度計、湿度計から構成することもできる。

第１演算部１０４は、第１測定部１０１が測定した心拍数より測定対象者の運動により発生する熱量を求める。

第２演算部１０５は、時刻ｔにおいて第１演算部１０４が求めた熱量、第２測定部１０２が測定した温度と第３測定部１０３が測定した湿度とから求められる測定対象者の雰囲気と測定対象者の皮膚との間の第１熱交換量、測定対象者の皮膚と深部との間の第２熱交換量、および測定対象者の発汗量より求める汗の蒸発量から、測定対象者の深部体温の変化を推定する。

第１演算部１０４、第２演算部１０５は、例えば、携帯端末装置などの、第１演算部１０４、第２演算部１０５を実現するためのプログラムにより動作する、マイクロプロセッサを備えた小型のコンピュータ機器から構成することができる。

次に、本発明の実施の形態に係る深部体温推定方法について、図２を参照して説明する。まず、第１ステップＳ１０１で、第１測定部１０１が、測定対象者の心拍数を測定する。次に、第２ステップＳ１０２で、第２測定部１０２が、測定対象者の近傍の温度を測定する。次に、第３ステップＳ１０３で、第３測定部１０３が、測定対象者の近傍の湿度を測定する。次に、第４ステップＳ１０４で、第１演算部１０４が、第１ステップＳ１０１で第１測定部１０１が測定した心拍数より、測定対象者の運動により発生する熱量を求める。

次に、第５ステップＳ１０５で、第２演算部１０５が、熱量、測定対象者の雰囲気と測定対象者の皮膚との間の第１熱交換量、測定対象者の皮膚と深部との間の第２熱交換量、および測定対象者の汗の蒸発量から、測定対象者の深部体温の変化を推定する。熱量は、時刻ｔにおいて第４ステップＳ１０４で第１演算部１０４が求める。第１熱交換量は、第２ステップＳ１０２で第２測定部１０２が測定した温度と、第３ステップＳ１０３で第３測定部１０３が測定した湿度と、から求められる。

以下、より詳細に説明する。

人体を表皮層（皮膚）と非表皮層（深部）と簡易化して考え，深部体温は非表皮層と同じであると仮定すると、式（１）、式（２）に示すように、時刻tにおける皮膚層と深部層の熱の流入出に関する方程式はそれぞれ式（１）、式（２）のように記述できる。

式（１）の左辺は、時刻ｔ−Δｔからｔに推移した際の、皮膚の温度上昇により生じる熱量を示している。式（１）の右辺の各項は、時刻ｔ−Δｔからｔに推移した際の、皮膚、深部の温度上昇に起因する、深部より皮膚へ流入する熱量、皮膚の代謝で生じる熱量、皮膚より大気へ流出する熱量、皮膚表面での発汗により生じる蒸発熱を示している。

式（２）の左辺は、時刻ｔ−Δｔからｔに推移した際の、深部の温度上昇により生じる熱量を示している。式（２）の右辺の各項は、時刻ｔ−Δｔからｔに推移した際の、皮膚、深部の温度上昇に起因する、皮膚より深部へ流入する熱量、深部の代謝で生じる熱量、測定対象者の運動により生じる熱量を示している。

Ｔ_i,t[℃]（ｉ＝ｓｋｉｎ，ｂｏｄｙ）は、時刻ｔにおける皮膚、深部の温度を示している。Ｔ_i,tは、各部の温度および、時刻がΔｔに推移した際の各部の熱量を、上述の２つの式より算出する際に使用する。熱量については、後述する。

Ｗ_i［ｋｇ］（ｉ＝ｓｋｉｎまたはｂｏｄｙ）は、皮膚または深部の質量を示している。Ｃ_i［Ｊ／(ｋｇ・℃)］（ｉ＝ｓｋｉｎまたはｂｏｄｙ）は、皮膚または深部の比熱を示している。Ｍ_i［Ｗ］（ｉ＝ｓｋｉｎまたはｂｏｄｙ）は、皮膚または深部の代謝量を示している。各部位の質量と比熱は、情報通信研究機構（ＮＩＣＴ）が開発した「ＴＡＲＯ」や「ＨＡＮＡＫＯ」のような、解剖学的人体モデルや文献値を使用することができる。各部位の代謝量に関しては、後述する。

ｈｘ［Ｗ／℃］は、皮膚と深部との間の熱交換係数であり、時刻ｔにおけるｈｘは、ｈｘ＝ｈｘ₀+ｈｘ₁×２^Aより算出する。なお、Ａ＝２（Ｔ_body,t-Δt−Ｔ_body,0)／１．２＋｛ＥＸ（ｔ）−１）／１．５である。ｈｘ₀およびｈｘ₁は、実験により決定する。例えば、ｈｘ₀＝１７．０［Ｗ／℃］、ｈｘ₁＝２３．５［Ｗ／℃］を使用することができる。

ＥＸ(ｔ)［Ｗ］は，時刻ｔ−Δｔからｔに推移した際の、測定対象者の運動により生じる熱量を示しており，第１演算部１０４が求めた心拍数より算出するが、詳細は後述する。

ＳＷ（ｔ）［Ｗ］は、時刻ｔ−Δｔからｔに推移した際の、発汗により生じる蒸発熱である。ＳＷ（ｔ）は、第３測定部１０３が測定した、湿度ｈｕｍｉｄｉｔｙ（ｔ）および，各部の温度より算出するが詳細は後述する。

Ｔ_air［℃］は、周辺温度を示しており，第２測定部１０２が測定した温度を使用する。式（１）の右辺第３項に示すように，測定対象者の皮膚と雰囲気との間の熱伝達率Ｈ［Ｗ／（ｍ²・℃）］（詳細は後述）を用いて、「Ｈ・（Ｔ_skin,t−Ｔ_air）」を全皮膚面積で積分することにより求めることができる。また、皮膚面積Ｓは、例えば，Ｄｕｂｏｉの式「Ｓ＝身長［ｃｍ］^0.725×体重［ｋｇ］^0.425×７．１８４×１０^(-3)」などの推定式を用いて求めることができる。

上記式を変形すると、以下に示す、式（３），式（４）のようになるので、上述した各定数，時刻ｔ-Δｔにおける測定対象者の深部体温のＴ_body,t-Δt［℃］、測定対象者の皮膚温の初期値Ｔ_skin,t-Δt［℃］，各測定部から算出した心拍数，温度，湿度を用いて、時刻ｔにおける、皮膚温度および深部温度を推定することができる。

「ｈｘ×（Ｔ_body,m-1−Ｔ_skin,m-1）」，「ｈｘ×（Ｔ_skin,m-1−Ｔ_body,m-1）」は、皮膚と深部との間の熱交換（第２熱交換量）を示している。また、「∫_SＨ×（Ｔ_skin,m-1−Ｔ_air）」は、皮膚と雰囲気との熱交換（第１熱交換量）を示している。雰囲気との第１熱交換量、汗の蒸発により失われる熱量は、皮膚温度の計算に用いる。運動により発生する熱量と、代謝によって発生する熱量は、深部体温の計算に用いる。また、皮膚と深部とでの第２熱交換量を、皮膚温度と深部温度のそれぞれの計算に用いる。なお、外気温と熱交換量を計算する項の積分（第２演算部の演算）において、皮膚が外気に露出している部分と露出していない部分とに分けて行い、それぞれの熱伝達率（第１熱交換量、汗の蒸発量）を変化させることもできる。

また、各定数は上記や後述の値に限定したものではなく、実際に即した値を用いることができ、計算実施において適宜設定することができる。各部の温度の初期値は、実測しておくことができる。また、各部の温度の初期値は、上記同様、実際に即した値を用いることができ、計算実施において適宜設定することができる。

次に、第２演算部１０５の計算で用いられる各パラメータについて説明する。

発汗の蒸発で奪われる熱量ＳＷ（ｔ）は、不感蒸散と有感蒸散の和Ｅ［Ｗ］から算出できる。「Ｅ＝Ｑ×（ＩＰ＋ｓｗｒａｔｅ（ｔ））」となる。ここで、ＩＰ［ｇ／ｓ］は不感蒸散、Ｑ［Ｊ／ｇ］は水の蒸発熱、ｓｗｒａｔｅ（ｔ）は、定数α１０、α１１、α２０、α２１、β１０、β１１、β２０、β２１を用い、「ｓｗｒａｔｅ（ｔ）＝［α１１ｔａｎｈ｛β１１（Ｔ_skin,t-Δt−Ｔ_skin,0）−β１０｝＋α１０］×（Ｔ_skin,t-Δt−Ｔ_skin,0）／ｍｉｎｕｔｅ＋［α２１ｔａｎｈ｛β２１（Ｔ_body,t-Δt−Ｔ_body,0）−β２０｝＋α２０］×（Ｔ_body,t-Δt−Ｔ_body,0）／ｍｉｎｕｔｅ」となる。ｍｉｔｕｔｅは分である。

例えば，次に示す値を設定することができる。「α１０＝０．９５，α１１＝０．５５，α２０＝３．８０，α２１＝３．２０，β１０＝０．０９，β１１＝０．５９，β２０＝１．８０，β２１＝２．７０」。

ただし、ＳＷ（ｔ)は、湿潤率ｗｅｔが１を超える場合には最大蒸発熱Ｅ_max［Ｗ］を超えることができない。Ｗｅｔ＝Ｅ／Ｅ_maxで計算できる。Ｅ_maxの計算は空気の動きＶ_air［ｍ／ｓ］に依存する対流による熱移動Ｈ_c［Ｗ×ｍ²／ｋ］、皮膚の飽和水蒸気圧Ｐ_s［ｋＰａ］、空気の飽和水蒸気圧Ｐ_air［ｋＰａ］から計算できる。式は下記の通りである。

・Ｈ_c＝３．０×√（１０×Ｖ_air）

・Ｐ_s＝７．５×（１０^D）
なお、Ｄ＝７．２３−（（１７５０．３／（２７３＋Ｔ_skin,t-Δt−３８．１）））である。

・Ｐ_air＝７．５×（１０^G）
なお、Ｇ＝７．２３−（（１７５０．３／（２７３＋Ｔ_air−３８．１）））である。

・Ｅ_max＝Ｅ_{max_coef}×Ｈ_c×（Ｐ_s−ｈｕｍｉｄｉｔｙ（ｔ−Δｔ）×Ｐａ）×｛ｃｏｖｅｒａｇｅ＋（１−ｃｏｖｅｒａｇｅ）／ｆｐｃｌ｝

なお、Ｅ_{max_coef}は係数、ｃｏｖｅｒａｇｅは衣服の被覆率、ｆｐｃｌは透湿性の係数である。

運動により発生する熱量Ｅｘは、ＭＥＴｓにより計算できる。例えば、「ＭＥＴｓ＝（測定された心拍数［ｂｐｍ］−測定対象者の安静時心拍数［ｂｐｍ］）÷（測定対象者の最大心拍数［ｂｐｍ］−安静時心拍数［ｂｐｍ］）×最大酸素摂取量［ｍｌ］÷３．５［ｍｌ］」によりＭＥＴｓ数を求める。求めたＭＥＴｓ数を用い、「ＥＸ熱量［Ｗ］＝ＭＥＴｓ数×１．０５×体重［ｋｇ］×４１８４÷３６００」により、熱量を求める。最大酸素摂取量は、実測した値を用いることができる。また、推定式（ＶＯ２_max＝１５ＨＲ_max／ＨＲ_rest）より、最大酸素摂取量を推定することもできる。なお、ＨＲ_restは、実測した安静時心拍数である。また、ＨＲ_maxは、実測、あるいは年齢を用いた推定式「最大心拍数＝（２０８−０．７×年齢）ｂｐｍ」による最大心拍数である。

熱伝達率Ｈは、対流による熱伝達Ｈｃと、放射による熱伝達Ｈｒとにより計算できる。「Ｈ＝（Ｈｃ＋Ｈｒ）×ｆｃｌ」。ｆｃｌは、被服係数である。

基礎代謝は、測定対象者の皮膚代謝Ｍ_skin［Ｗ］、測定対象者の深部代謝Ｍ_body［Ｗ］と分けて、時刻ｔにおけるＭ_skinおよびＭ_bodyは、Ｍ_skin＝Ｍ_skin,0×１．１^K、Ｍ_body＝Ｍ_body、0×１．１^Lより計算する。なお、「Ｋ＝（Ｔ_skin,t-Δt−Ｔ_skin,0）」、「Ｌ＝（Ｔ_body,t-Δt−Ｔ_body,0）」である。

ここで、「Ｍ_body,0＝（０．１２３８＋０．０４８１×体重×０．９２＋０．０２３４×ｈｅｉｇｈｔ−０．０１３８×ａｇｅ−０．５４７３×ｓｅｘ_coef）×１００００００／２４／３６００×ａｃｔｉｖｉｔｙ_level」となる。ｓｅｘ_coefは性別ごとの係数であり、男性が１、女性が２となる。ａｃｔｉｖｉｔｙ_levelは身体活動レベルである。また、Ｍ_skin,0＝Ｍ_skin／ＲＯＵ_skin×Ｗ_skinで計算する。ＲＯＵ_skinは、皮膚の密度［ｋｇ／ｍ³］、Ｍ_skinは、皮膚の代謝［Ｗ／ｍ³］、Ｗ_skinは、体重×０．０８で計算する。Ｍ_body,0は、上述した測定対象者の性別、身体活動レベル、体重の他に、身長、年齢などに応じて設定することができる。

以下、実施例を用いて本発明について説明する。

［実施例１］
はじめに、実施例１について、図３を参照して説明する。図３は、実施の形態に係る深部体温推定方法を実施するためのウエアラブルデバイスの構成を示している。このウエアラブルデバイスは、測定対象者２０１が着衣するインナーウエア２０２に装着された測定部２０３を備える。測定部２０３は、アウターウエア２０４とインナーウエア２０２との間に配置される。測定部２０３は、一体化された温湿度センサとウエア型の心電計とから構成されている（参考文献）。測定部２０３により、インナーウエア２０２を着衣している測定対象者２０１の心拍数、および近傍の温度・湿度が測定される。

測定部２０３は、通信機能を内蔵し、携帯端末装置２０５と通信可能とされている。測定部２０３が測定した心拍数、温度、湿度のデータは、携帯端末装置２０５に送信される。携帯端末装置２０５では、インストールされたソフトウエアにより、第１演算部および第２演算部の機能を実現する。

測定部２０３より送信された心拍数、温度、湿度のデータを受信した携帯端末装置２０５は、まず、受信した心拍数のデータより、熱量を求める。次に、携帯端末装置２０５は、受信した温度、湿度のデータより、測定対象者２０１の雰囲気と測定対象者２０１の皮膚との間の第１熱交換量を求める。また、携帯端末装置２０５は、測定対象者２０１の皮膚と深部との間の第２の熱交換量を求める。また、携帯端末装置２０５は、測定対象者２０１の発汗量より求められる汗の蒸発量を求める。これらの各量を用い、携帯端末装置２０５は、測定対象者２０１の深部体温の変化を計算（推定）する。

［実施例２］
次に、実施例２について、図４を参照して説明する。図４は、実施の形態に係る深部体温推定方法を実施するためのウエアラブルデバイスの構成を示している。このウエアラブルデバイスは、測定対象者２０１が着衣するインナーウエア２０２とアウターウエア２０４との間に配置される第１測定部２０３ａを備える。第１測定部２０３ａは、ベルト型の心電計から構成され、ベルト２０６により測定対象者２０１の胸部に装着される。第１測定部２０３ａにより、測定対象者２０１の心拍数が測定される。

また、このウエアラブルデバイスは、測定対象者２０１の頭部に配置される第２測定部２０３ｂを備える。第２測定部２０３ｂは、測定対象者２０１が装着するヘルメット２０７に固定されている。第２測定部２０３ｂにより、測定対象者２０１の近傍の温度・湿度が測定される。

第１測定部２０３ａは、通信機能を内蔵し、携帯端末装置２０５と通信可能とされている。第１測定部２０３ａが測定した心拍数のデータは、携帯端末装置２０５に送信される。また、第２測定部２０３ｂも通信機能を有し、携帯端末装置２０５と通信可能とされている。第２測定部２０３ｂ測定した温度、湿度のデータは、携帯端末装置２０５に送信される。携帯端末装置２０５では、インストールされたソフトウエアにより、第１演算部および第２演算部の機能を実現する。

第１測定部２０３ａより送信された心拍数のデータ、第２測定部２０３ｂより送信された温度、湿度のデータを受信した携帯端末装置２０５は、まず、受信した心拍数のデータより、熱量を求める。次に、携帯端末装置２０５は、受信した温度、湿度のデータより、測定対象者２０１の雰囲気と測定対象者２０１の皮膚との間の第１熱交換量を求める。また、携帯端末装置２０５は、測定対象者２０１の皮膚と深部との間の第２の熱交換量を求める。また、携帯端末装置２０５は、測定対象者２０１の発汗量より求められる汗の蒸発量を求める。これらの各量を用い、携帯端末装置２０５は、測定対象者２０１の深部体温の変化を計算（推定）する。

以上に説明したように、本発明によれば、熱量、測定対象者の雰囲気と測定対象者の皮膚との間の第１熱交換量、測定対象者の皮膚と深部との間の第２熱交換量、および測定対象者の汗の蒸発量から、測定対象者の深部体温の変化を推定するので、より容易に深部温度が推定できるようになる。

なお、本発明は以上に説明した実施の形態に限定されるものではなく、本発明の技術的思想内で、当分野において通常の知識を有する者により、多くの変形および組み合わせが実施可能であることは明白である。

［参考文献］スマートヘルスケアに向けた心電，加速度，温度・湿度の計測を可能にする低電力・小型ウエアラブルセンサを開発 NＴＴ，［２０２０年４月２１日検索］、（https://www.ntt.co.jp/news2019/1911/191108a.html）。

１０１…第１測定部、１０２…第２測定部、１０３…第３測定部、１０４…第１演算部、１０５…第２演算部。

Claims (8)

1. 測定対象者の心拍数を測定する第１ステップと、
   前記測定対象者の近傍の温度を測定する第２ステップと、
   前記測定対象者の近傍の湿度を測定する第３ステップと、
   前記第１ステップで測定した心拍数より前記測定対象者の運動により発生する熱量を求める第４ステップと、
   時刻ｔにおいて前記第４ステップで求めた熱量、前記第２ステップで測定した温度と前記第３ステップで測定した湿度とから求められる前記測定対象者の雰囲気と前記測定対象者の皮膚との間の第１熱交換量、前記測定対象者の皮膚と深部との間の第２熱交換量、および前記測定対象者の発汗量より求める汗の蒸発量から、前記測定対象者の深部体温の変化を推定する第５ステップと
   を備える深部体温推定方法。
2. 測定対象者の心拍数を測定する第１測定部と、
   前記測定対象者の近傍の温度を測定する第２測定部と、
   前記測定対象者の近傍の湿度を測定する第３測定部と、
   前記第１測定部が測定した心拍数より前記測定対象者の運動により発生する熱量を求める第１演算部と、
   時刻ｔにおいて前記第１演算部が求めた熱量、前記第２測定部が測定した温度と前記第３測定部が測定した湿度とから求められる前記測定対象者の雰囲気と前記測定対象者の皮膚との間の第１熱交換量、前記測定対象者の皮膚と深部との間の第２熱交換量、および前記測定対象者の発汗量より求める汗の蒸発量から、前記測定対象者の深部体温の変化を推定する第２演算部と
   を備える深部体温推定装置。
3. 請求項２記載の深部体温推定装置において、
   前記第２演算部は、前記測定対象者の皮膚が雰囲気に露出している部分と露出していない部分とで第１熱交換量、汗の蒸発量を変化させることを特徴とする深部体温推定装置。
4. 請求項２または３記載の深部体温推定装置において、
   前記第２演算部は、時刻ｔにおいて、前記第１演算部が求めた前記測定対象者の運動により発生する熱量ＥＸ［Ｗ］、前記第２測定部が測定した前記測定対象者の近傍の温度Ｔ_a（ｔ）［℃］、時刻ｔにおける前記測定対象者の発汗の蒸発で奪われる熱量ＳＷ（ｔ）［Ｗ］、前記測定対象者の深部体温の初期値Ｔ_body,0［℃］、皮膚温の初期値Ｔ_skin,0［℃］、前記測定対象者の体表面積Ｓ［ｍ²］、前記測定対象者の皮膚／雰囲気間の熱伝達率Ｈ［Ｗ／（ｍ²・℃）］、前記測定対象者の皮膚の比熱Ｃ_skin［Ｊ／（ｋｇ・℃）］、前記測定対象者の深部の比熱Ｃ_body［Ｊ／（ｋｇ・℃）］、前記測定対象者の皮膚の重量Ｗ_skin［ｋｇ］、前記測定対象者の深部の重量Ｗ_body［ｋｇ］、時刻ｔ−Δｔにおいて、前記第２演算部より算出される皮膚／深部間の熱交換係数ｈｘ［Ｗ／℃］、前記測定対象者の皮膚の代謝Ｍ_skin［Ｗ］、深部の代謝Ｍ_body［Ｗ］、皮膚温度Ｔ_skin,t-Δt［℃］、深部温度Ｔ_core,t-Δt［℃］を用い、以下の式より、時刻ｔにおける前記測定対象者の皮膚温度Ｔ_skin,t-Δt［℃］、深部温度Ｔ_core,t-Δt［℃］を推定することを特徴とする深部体温推定装置。
   

   
5. 請求項２または３記載の深部体温推定装置において、
   前記第２演算部は、前記測定対象者の発汗の蒸発で奪われる熱量を、前記第２演算部で算出する不感蒸散と有感蒸散の和が前記第２演算部で算出する最大蒸発熱以下の場合には不感蒸散と有感蒸散の和とし、不感蒸散と有感蒸散の和が最大蒸発熱を上回る場合には最大蒸発熱とすることを特徴とする深部体温推定装置。
6. 請求項５記載の深部体温推定装置において、
   前記第２演算部は、設定されている定数α１０，α１１，α２０，α２１，β１０，β１１，β２０，β２１および分を表すｍｉｎｕｔｅ［ｍｉｎ］用いて、「ｓｗｒａｔｅ（ｔ）＝［α１１ｔａｎｈ｛β１１（Ｔ_skin,t-Δt−Ｔ_skin,0）−β１０｝＋α１０］×（Ｔ_skin,t-Δt−Ｔ_skin,0）／ｍｉｎｕｔｅ＋［α２１ｔａｎｈ｛β２１（Ｔ_body,t-Δt−Ｔ_body,0）−β２０｝＋α２０］×（Ｔ_body,t-Δt−Ｔ_body,0）／ｍｉｎｕｔｅ」より算出されるｓｗｒａｔｅ（ｔ）［ｇ／ｓ］、水の蒸発熱Ｑ［Ｊ／ｇ］、不感蒸散ＩＰ［ｇ／ｓ］を用いて、「Ｅ＝Ｑ×（ＩＰ＋ｓｗｒａｔｅ（ｔ））」より、不感蒸散と有感蒸散の和Ｅ［Ｗ］を求め、
   空気の対流による熱移動Ｈ_c［Ｗ・ｍ²／ｋ］、皮膚の飽和水蒸気圧Ｐ_s［ｋＰａ］、空気の飽和水蒸気圧Ｐ_air［ｋＰａ］、設定されている定数Ｅ_{max_coef}、衣服の被覆率ｃｏｖｅｒａｇｅ、透湿性の係数ｆ_pcl、時刻ｔ−Δｔにおいて第３測定部より測定された湿度ｈｕｍｉｄｉｔｙを用いて、「Ｅ_max＝Ｅ_{max_coef}×Ｈ_c×（Ｐ_s−ｈｕｍｉｄｉｔｙ（ｔ−Δｔ）×Ｐａ）×｛ｃｏｖｅｒａｇｅ＋（１−ｃｏｖｅｒａｇｅ）／ｆｐｃｌ｝」より、最大蒸発熱Ｅ_max［Ｗ］を求めることを特徴とする深部体温推定装置。
7. 請求項４記載の深部体温推定装置において、
   前記第２演算部は、設定されている定数Ｍ_skin,0、Ｍ_body,0を用い、Ｋ＝（Ｔ_skin,t-Δt−Ｔ_skin,0）、Ｌ＝（Ｔ_body,t-Δt−Ｔ_body,0）として、Ｍ_skin＝Ｍ_skin,0×１．１^K、Ｍ_body＝Ｍ_body、0×１．１^Lより、前記測定対象者の皮膚の代謝Ｍ_skin［Ｗ］、深部の代謝Ｍ_body［Ｗ］を求めることを特徴とする深部体温推定装置。
8. 請求項７記載の深部体温推定装置において、
   前記第２演算部は、前記測定対象者の性別、身長、体重、年齢、身体活動レベルのうちの少なくとも一つの身体的特徴に応じて、定数Ｍ_body,0を設定することを特徴とする深部体温推定装置。

Images