JP7496551B2 - 深部体温推定装置、深部体温推定方法および深部体温推定プログラム - Google Patents

深部体温推定装置、深部体温推定方法および深部体温推定プログラム Download PDF

Info

Publication number
JP7496551B2
JP7496551B2 JP2021007047A JP2021007047A JP7496551B2 JP 7496551 B2 JP7496551 B2 JP 7496551B2 JP 2021007047 A JP2021007047 A JP 2021007047A JP 2021007047 A JP2021007047 A JP 2021007047A JP 7496551 B2 JP7496551 B2 JP 7496551B2
Authority
JP
Japan
Prior art keywords
subject
temperature
deep
skin
amount
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2021007047A
Other languages
English (en)
Other versions
JP2022111550A (ja
Inventor
優生 橋本
啓 桑原
浩芳 都甲
和彦 高河原
隆子 石原
晃正 平田
涼太 上松
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Telegraph and Telephone Corp
Nagoya Institute of Technology NUC
Original Assignee
Nippon Telegraph and Telephone Corp
Nagoya Institute of Technology NUC
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 Nippon Telegraph and Telephone Corp, Nagoya Institute of Technology NUC filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2021007047A priority Critical patent/JP7496551B2/ja
Publication of JP2022111550A publication Critical patent/JP2022111550A/ja
Application granted granted Critical
Publication of JP7496551B2 publication Critical patent/JP7496551B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Temperature Or Quantity Of Heat (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

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＞
本発明は、上記課題を解決するためになされたもので、測定対象者の深部体温を従来よりも高精度に推定することができる深部体温推定装置、深部体温推定方法および深部体温推定プログラムを提供することを目的とする。
本発明の深部体温推定装置は、測定対象者の心拍数または脈拍数を測定するように構成された第１測定部と、前記測定対象者の近傍の温度を測定するように構成された第２測定部と、前記測定対象者の近傍の湿度を測定するように構成された第３測定部と、前記測定対象者の心拍数または脈拍数と前記測定対象者の近傍の温度と前記測定対象者の近傍の湿度とに基づいて、前記測定対象者の第１部位の深部に流入出する熱量と前記第１部位の皮膚に流入出する熱量と前記測定対象者の第２部位の深部に流入出する熱量と前記第２部位の皮膚に流入出する熱量とを算出するように構成された熱量算出部と、前記熱量に基づいて、前記測定対象者の皮膚温度および深部温度を算出するように構成された温度算出部とを備え、前記第１部位は前記測定対象者の体幹部で、前記第２部位は前記測定対象者の四肢部、または前記第１部位は前記測定対象者の上半身で、前記第２部位は前記測定対象者の下半身であり、前記熱量算出部は、前記測定対象者の心拍数または脈拍数に基づいて、前記測定対象者の運動による前記第１部位と前記第２部位のそれぞれの発熱量を算出するように構成された発熱量算出部と、前記測定対象者の近傍の温度と前記温度算出部によって直前に算出された皮膚温度とに基づいて、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける皮膚と外気間の第１熱交換量を算出するように構成された第１熱交換量算出部と、前記測定対象者の心拍数または脈拍数と前記温度算出部によって直前に算出された皮膚温度と深部温度とに基づいて、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける深部と皮膚間の第２熱交換量を算出するように構成された第２熱交換量算出部と、前記温度算出部によって直前に算出された皮膚温度と深部温度とに基づいて、前記測定対象者の深部における前記第１部位と前記第２部位間の第３熱交換量を算出するように構成された第３熱交換量算出部とを含むことを特徴とするものである。
また、本発明の深部体温推定装置の１構成例において、前記第１熱交換量算出部は、前記測定対象者の皮膚が衣服で覆われている割合を考慮して、前記第１熱交換量を算出することを特徴とするものである。
また、本発明の深部体温推定装置の１構成例において、前記熱量算出部は、前記測定対象者の近傍の温度と前記測定対象者の近傍の湿度と前記温度算出部によって直前に算出された皮膚温度と深部温度とに基づいて、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける皮膚蒸散量を算出するように構成された皮膚蒸散量算出部をさらに含むことを特徴とするものである。
また、本発明の深部体温推定装置の１構成例において、前記熱量算出部は、外気と前記測定対象者の深部間の第４熱交換量を設定するように構成された第４熱交換量設定部と、前記温度算出部によって直前に算出された皮膚温度と深部温度とに基づいて、前記測定対象者の第１部位皮膚層の代謝量と第１部位深部層の代謝量と第２部位皮膚層の代謝量と第２部位深部層の代謝量を算出するように構成された代謝量算出部とをさらに含むことを特徴とするものである。
また、本発明の深部体温推定装置の１構成例において、前記温度算出部は、前記測定対象者の前記第１部位と前記第２部位の発熱量Ｑ１，Ｕ，Ｑ１，Ｌ、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける皮膚と外気間の第１熱交換量Ｑ３，Ｕ，Ｑ３，Ｌ、 前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける深部と皮膚間の第２熱交換量Ｑ６，Ｕ，Ｑ６，Ｌ、前記測定対象者の深部における前記第１部位と前記第２部位間の第３熱交換量Ｑ７、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける皮膚蒸散量Ｑ４，Ｕ，Ｑ４，Ｌ、外気と前記測定対象者の深部間の第４熱交換量Ｑ５、前記測定対象者の第１部位皮膚層の代謝量Ｑ２，ＵＳ、第１部位深部層の代謝量Ｑ２，ＵＣ、第２部位皮膚層の代謝量Ｑ２，ＬＳ、第２部位深部層の代謝量Ｑ２，ＬＣ、前記測定対象者の第１部位皮膚層の熱容量ＷＣＵＳ、第１部位深部層の熱容量ＷＣＵＣ、第２部位皮膚層の熱容量ＷＣＬＳ、第２部位深部層の熱容量ＷＣＬＣ、直前に算出した第１部位皮膚層の温度の推定値ＴＵＳ［ｔ］、第１部位深部層の温度の推定値ＴＵＣ［ｔ］、第２部位皮膚層の温度の推定値ＴＬＳ［ｔ］、第２部位深部層の温度の推定値ＴＬＣ［ｔ］に基づいて、時刻（ｔ＋Δｔ）における第１部位皮膚層の温度の推定値ＴＵＳ［ｔ＋Δｔ］、第１部位深部層の温度の推定値ＴＵＣ［ｔ＋Δｔ］、第２部位皮膚層の温度の推定値ＴＬＳ［ｔ＋Δｔ］、第２部位深部層の温度の推定値ＴＬＣ［ｔ＋Δｔ］を算出することを特徴とするものである。
また、本発明の深部体温推定方法は、測定対象者の心拍数または脈拍数を測定する第１ステップと、前記測定対象者の近傍の温度を測定する第２ステップと、前記測定対象者の近傍の湿度を測定する第３ステップと、前記測定対象者の心拍数または脈拍数と前記測定対象者の近傍の温度と前記測定対象者の近傍の湿度とに基づいて、前記測定対象者の第１部位の深部に流入出する熱量と前記第１部位の皮膚に流入出する熱量と前記測定対象者の第２部位の深部に流入出する熱量と前記第２部位の皮膚に流入出する熱量とを算出する第４ステップと、前記熱量に基づいて、前記測定対象者の皮膚温度および深部温度を算出する第５ステップとを含み、前記第１部位は前記測定対象者の体幹部で、前記第２部位は前記測定対象者の四肢部、または前記第１部位は前記測定対象者の上半身で、前記第２部位は前記測定対象者の下半身であり、前記第４ステップは、前記測定対象者の心拍数または脈拍数に基づいて、前記測定対象者の運動による前記第１部位と前記第２部位のそれぞれの発熱量を算出するステップと、前記測定対象者の近傍の温度と前記第５ステップによって直前に算出された皮膚温度とに基づいて、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける皮膚と外気間の第１熱交換量を算出するステップと、前記測定対象者の心拍数または脈拍数と前記第５ステップによって直前に算出された皮膚温度と深部温度とに基づいて、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける深部と皮膚間の第２熱交換量を算出するステップと、前記第５ステップによって直前に算出された皮膚温度と深部温度とに基づいて、前記測定対象者の深部における前記第１部位と前記第２部位間の第３熱交換量を算出するステップとを含むことを特徴とするものである。
また、本発明の深部体温推定プログラムは、前記の各ステップをコンピュータに実行させることを特徴とするものである。
本発明によれば、測定対象者の心拍数と測定対象者近傍の温度だけでなく、測定対象者近傍の湿度を測定して、測定対象者の第１部位と第２部位の深部と皮膚とに流入出する熱量を算出し、算出した熱量から測定対象者の深部体温を推定するので、従来よりも高精度に深部体温を推定することができる。
図１は、本発明の実施例に係る深部体温推定装置の構成を示すブロック図である。 図２は、測定対象者の身体の２部位２層モデルと、測定対象者の各部位・各層に流入出する熱量を示す図である。 図３は、本発明の実施例に係る深部体温推定装置の動作を説明するフローチャートである。 図４は、本発明の実施例に係る深部体温推定装置を実現するコンピュータの構成例を示すブロック図である。
以下、本発明の実施例について図面を参照して説明する。図１は本発明の実施例に係る深部体温推定装置の構成を示すブロック図である。深部体温推定装置は、第１測定部１と、第２測定部２と、第３測定部３と、発熱量算出部４と、代謝量算出部５と、熱伝達・熱放射量算出部６（第１熱交換量算出部）と、皮膚蒸散量算出部７と、呼気蒸散量設定部８（第４熱交換量設定部）と、熱交換量算出部９（第２熱交換量算出部）と、熱交換量算出部１０（第３熱交換量算出部）と、温度算出部１１とを備えている。
発熱量算出部４と代謝量算出部５と熱伝達・熱放射量算出部６と皮膚蒸散量算出部７と呼気蒸散量設定部８と熱交換量算出部９と熱交換量算出部１０とは、熱量算出部１２を構成している。
第１測定部１は、測定対象者の心拍数または脈拍数を測定する。第１測定部１は、例えば、測定対象者の心電位を計測するウエア型あるいはベルト型の心電計と、心電計が計測した心電位から心拍数を算出する算出部とから構成される。あるいは、第１測定部１は、測定対象者の脈波を計測する腕時計型あるいはイヤホン型の脈波計と、脈波計が計測した脈波から心拍数（脈拍数）を算出する算出部とから構成される。
第２測定部２は、測定対象者近傍の温度（測定対象者の雰囲気の気温）を測定する。第２測定部２は、例えば、温度計から構成される。あるいは、第２測定部２は、外部の気象システムから測定対象者の近傍の気象データを取得するようにしてもよい。
第３測定部３は、測定対象者近傍の湿度（測定対象者の雰囲気の湿度）を測定する。第３測定部３は、例えば、湿度計から構成される。あるいは、第３測定部３は、第２測定部２と同様に、外部の気象システムから測定対象者の近傍の気象データを取得するようにしてもよい。
第１測定部１と第２測定部２と第３測定部３とは、一体型の測定器でもよいし、個別の測定器でもよい。
熱量算出部１２は、測定対象者の心拍数と測定対象者近傍の温度と測定対象者近傍の湿度とに基づいて、測定対象者の体幹部（第１部位）と四肢部（第２部位）の深部と皮膚とに流入出する熱量を算出する。
温度算出部１１は、熱量算出部１２によって算出された熱量に基づいて、測定対象者の皮膚温度および深部温度を算出する。
［深部体温推定方法：時間ステップごとの深部体温変化の計算方法］
本発明は、図２に示すように、測定対象者の身体が体幹部Ｕと四肢部Ｌの２部位で構成され、体幹部Ｕと四肢部Ｌの２部位がそれぞれ深部層Ｃと皮膚層Ｓの２層を有するものとみなす。すなわち、測定対象者の身体は、体幹部深部層ＵＣと、体幹部皮膚層ＵＳと、四肢部深部層ＬＣと、四肢部皮膚層ＬＳとからなる。本発明は、取得する温湿度と心拍数情報とを基に、測定対象者の各部・各層に流入出する、時々刻々と変化する熱量を算出し、熱量の変化から、測定対象者の各部位・各層の温度変化を推定する。
本実施例では、センサデータを取得する時間ステップごとの深部体温変化の計算方法を説明する。本実施例では、上述のとおり、測定対象者の身体を体幹部と四肢部の２部位で構成されるものとみなした場合の計算例を記載しているが、体幹部を上半身に置き替え、四肢部を下半身に置き替えてもよい。つまり、測定対象者の身体を第１部位と第２部位の任意の２部位で構成されるものとしてよい。
式（１）～式（４）に、それぞれ測定対象者の体幹部皮膚層ＵＳの温度変化、体幹部深部層ＵＣの温度変化、四肢部皮膚層ＬＳの温度変化、四肢部深部層ＬＣの温度変化を推定する式の例を示す。ＴUS［ｔ］は時刻ｔにおける体幹部皮膚層ＵＳの温度［℃］、ＴUC［ｔ］は時刻ｔにおける体幹部深部層ＵＣの温度［℃］、ＴLS［ｔ］は時刻ｔにおける四肢部皮膚層ＬＳの温度［℃］、ＴLC［ｔ］は時刻ｔにおける四肢部深部層ＬＣの温度［℃］である。
Figure 0007496551000001
Figure 0007496551000002
Figure 0007496551000003
Figure 0007496551000004
Ｑ1,Uは測定対象者の運動による体幹部Ｕの発熱量［Ｗ］、Ｑ1,Lは運動による四肢部Ｌの発熱量［Ｗ］である。Ｑ2,USは測定対象者の体幹部皮膚層ＵＳの代謝量［Ｗ］、Ｑ2,UCは体幹部深部層ＵＣの代謝量［Ｗ］、Ｑ2,LSは四肢部皮膚層ＬＳの代謝量［Ｗ］、Ｑ2,LCは四肢部深部層ＬＣの代謝量［Ｗ］である。Ｑ3,U（第１熱交換量）は測定対象者の体幹部Ｕにおける皮膚と外気間の熱伝達・熱放射量［Ｗ］、Ｑ3,L（第１熱交換量）は四肢部Ｌにおける皮膚と外気間の熱伝達・熱放射量［Ｗ］である。
Ｑ4,Uは測定対象者の体幹部Ｕにおける皮膚蒸散量、Ｑ4,Lは四肢部Ｌにおける皮膚蒸散量である。Ｑ5（第４熱交換量）は測定対象者の呼気蒸散量［Ｗ］である。Ｑ6,U（第２熱交換量）は測定対象者の体幹部Ｕにおける深部と皮膚間の熱交換量［Ｗ］、Ｑ6,L（第２熱交換量）は四肢部Ｌにおける深部と皮膚間の熱交換量［Ｗ］である。Ｑ7（第３熱交換量）は測定対象者の深部における体幹と四肢間の熱交換量［Ｗ］である。
ＷＣUSは測定対象者の体幹部皮膚層ＵＳの熱容量［Ｊ／℃］、ＷＣUCは体幹部深部層ＵＣの熱容量［Ｊ／℃］、ＷＣLSは四肢部皮膚層ＬＳの熱容量［Ｊ／℃］、ＷＣLCは四肢部深部層ＬＣの熱容量［Ｊ／℃］である。Δｔは計算ステップ時間であり、例えば１［ｓ］以下とする。
また、時刻ｔ＋Δｔにおける平均皮膚温Ｔsk［ｔ＋Δｔ］、深部体温Ｔ［ｔ＋Δｔ］は式（５）、式（６）のようになる。
Figure 0007496551000005
Figure 0007496551000006
ｓｆ＿ｃｏｎｆ＿ＵＳは測定対象者の体表面全体に対する体幹部Ｕの表面積の割合［％］、ｓｆ＿ｃｏｎｆ＿ＬＳは体表面全体に対する四肢部Ｌの表面積の割合［％］である。熱容量ＷＣUS，ＷＣUC，ＷＣLS，ＷＣLC、割合ｓｆ＿ｃｏｎｆ＿ＵＳ，ｓｆ＿ｃｏｎｆ＿ＬＳは、それぞれ既知の値であり、実際に即した値を用いればよく、計算の実施において適宜設定すればよい。
図３は本実施例の深部体温推定装置の動作を説明するフローチャートである。第１測定部１は、測定対象者の心拍数を測定する（図３ステップＳ１００）。第２測定部２は、測定対象者近傍の温度（測定対象者の雰囲気の気温）を測定する（図３ステップＳ１０１）。第３測定部３は、測定対象者近傍の湿度（測定対象者の雰囲気の湿度）を測定する（図３ステップＳ１０２）。
［運動による発熱量Ｑ1の算出］
次に、発熱量算出部４は、第１測定部１によって測定された心拍数に基づいて、測定対象者の運動による体幹部Ｕの深部層における発熱量Ｑ1,U［Ｗ］、運動による四肢部Ｌの深部層における発熱量Ｑ1,L［Ｗ］をそれぞれ式（７）、式（８）のように算出する（図３ステップＳ１０３）。
Figure 0007496551000007
Figure 0007496551000008
式（７）、式（８）は、文献「“健康づくりのための運動指針２００６ ～生活習慣病予防のために～”，厚生労働省，２００６年，＜https://www.mhlw.go.jp/shingi/2006/07/dl/s0719-3c.pdf＞」に開示されている。ｗｅｉｇｈｔは測定対象者の重量［ｋｇ］、ｅｘ＿ｃｏｎｆ＿Ｕは測定対象者の全身に対する体幹部Ｕの筋肉量の割合［％］、ｅｘ＿ｃｏｎｆ＿Ｌは全身に対する四肢部Ｌの筋肉量の割合［％］である。重量ｗｅｉｇｈｔ、割合ｅｘ＿ｃｏｎｆ＿Ｕ，ｅｘ＿ｃｏｎｆ＿Ｌは、それぞれ既知の値であり、実際に即した値を用いればよく、計算の実施において適宜設定すればよい。
また、ＭＥＴｓ［ｔ］はメッツ数を表している。発熱量算出部４は、第１測定部１が取得した時刻ｔにおける心拍数ＨＲ［ｔ］［ｂｐｍ］、安静時心拍数ＨＲrest、最大心拍数ＨＲmaxを用いて、式（９）または式（１０）のいずれか一方の式によりＭＥＴｓ［ｔ］を算出すればよい。
Figure 0007496551000009
Figure 0007496551000010
式（１０）は、文献「J.R.Wicks，et al.，“HR Index-A Simple Method for the Prediction of Oxygen Uptake”，Medicine and Science in Sports and Exercise，2011」に開示されている。
安静時心拍数ＨＲrest、最大心拍数ＨＲmaxは、それぞれ過去の測定で得られた既知の値として、実際に即した値を用いればよく、計算の実施において適宜設定すればよい。また、式（７）、式（８）で用いるＭＥＴｓ［ｔ］に関しては、瞬時値だけでなく、時間平均値（例えば、６分間の平均値）を用いてもよい。
［代謝量Ｑ2の算出］
次に、代謝量算出部５は、温度算出部１１がΔｔ前に算出した、時刻ｔにおける測定対象者の平均皮膚温の推定値Ｔsk［ｔ］［℃］と深部体温の推定値Ｔ［ｔ］［℃］とに基づいて、測定対象者の体幹部皮膚層ＵＳの代謝量Ｑ2,US［Ｗ］、体幹部深部層ＵＣの代謝量Ｑ2,UC［Ｗ］、四肢部皮膚層ＬＳの代謝量Ｑ2,LS［Ｗ］、四肢部深部層ＬＣの代謝量Ｑ2,LC［Ｗ］をそれぞれ式（１１）～式（１４）のように算出する（図３ステップＳ１０４）。
Figure 0007496551000011
Figure 0007496551000012
Figure 0007496551000013
Figure 0007496551000014
式（１１）～式（１４）は、文献「Ronald J Spiegel，“A Review of Numerical Models for Predicting the Energy Deposition and Resultant Thermal Response of Humans Exposed to Electromagnetic Fields”，IEEE Transactions on Microwave Theory and Techniques，Volume32，Issue8，1984」に開示されている。Ａskinは代謝に関する定数、ｖｏｌｕｍｅ＿ＵＳは測定対象者の体幹部皮膚層ＵＳの体積［ｍ3］、ｖｏｌｕｍｅ＿ＬＳは四肢部皮膚層ＬＳの体積［ｍ3］、ｗｅｉｇｈｔ＿Ｕは測定対象者の体幹部Ｕの重量［ｋｇ］、ｗｅｉｇｈｔ＿Ｌは四肢部Ｌの重量［ｋｇ］である。定数Ａskin、体積ｖｏｌｕｍｅ＿ＵＳ，ｖｏｌｕｍｅ＿ＬＳ、重量ｗｅｉｇｈｔ＿Ｕ，ｗｅｉｇｈｔ＿Ｌは、それぞれ既知の値であり、実際に即した値を用いればよく、計算の実施において適宜設定すればよい。
Ｔsk［０］は平均皮膚温Ｔsk［ｔ］の演算時の初期値、Ｔ［０］は深部体温Ｔ［ｔ］の演算時の初期値である。本実施例では、体幹部皮膚層ＵＳの温度の初期値ＴUS［０］、体幹部深部層ＵＣの温度の初期値ＴUC［０］、四肢部皮膚層ＬＳの温度の初期値ＴLS［０］、四肢部深部層ＬＣの温度の初期値ＴLC［０］として、別の温度測定装置によって測定した実測値または一般的な安静時の既知の値などの実際に即した値を用いる。したがって、時刻ｔ＝０の初回の代謝量Ｑ2の算出では、Ｔsk［ｔ］＝Ｔsk［０］である。
同様に、深部体温Ｔ［ｔ］の初期値Ｔ［０］としては、別の温度測定装置によって測定した実測値または一般的な安静時の既知の値などの実際に即した値を用いる。したがって、時刻ｔ＝０の初回の代謝量Ｑ2の算出では、Ｔ［ｔ］＝Ｔ［０］である。
また、代謝量算出部５は、式（１２）、式（１４）のＭを式（１５）のように算出する。
Figure 0007496551000015
式（１５）は、文献「AA Ganpule，et al.，“Interindividual variability in sleeping metabolic rate in Japanese subjects”，European Journal of Clinical Nutrition，volume 61，2007」に開示されている。ｗｅｉｇｈｔ＿Ｃは測定対象者の深部の重量［ｋｇ］、ｈｅｉｇｈｔは測定対象者の身長［ｃｍ］、ａｇｅは測定対象者の年齢である。ｓｅｘｃｏｅｆは測定対象者が男性の場合に０．５４７３、測定対象者が女性の場合に０．５４７３×２となる定数である。ａｃｔｉｖｉｔｙ＿ｌｅｖｅｌは身体活動レベル、Ａcoefは代謝調整用のパラメータである。重量ｗｅｉｇｈｔ＿Ｃ、身長ｈｅｉｇｈｔ、年齢ａｇｅ、定数ｓｅｘｃｏｅｆ、ａｃｔｉｖｉｔｙ＿ｌｅｖｅｌ、パラメータＡcoefは、それぞれ既知の値であり、実際に即した値を用いればよく、計算の実施において適宜設定すればよい。
ａｃｔｉｖｉｔｙ＿ｌｅｖｅｌの設定例としては、測定対象者の生活の大部分が座位で、静的な活動が中心の場合は１．５程度とすればよい。測定対象者が座位中心の仕事をしているが、職場内での移動や立位での作業、接客等を含む場合、あるいは通勤、買物、家事、軽いスポーツ等のいずれかを含む場合は、ａｃｔｉｖｉｔｙ＿ｌｅｖｅｌを１．７５程度とすればよい。測定対象者が移動や立位の多い仕事への従事者である場合、あるいはスポーツなど余暇における活発な運動習慣を持っている場合は、ａｃｔｉｖｉｔｙ＿ｌｅｖｅｌを２．０程度とすればよい。このように、測定対象者の運動状況に応じてａｃｔｉｖｉｔｙ＿ｌｅｖｅｌを適宜設定すればよい。
［皮膚と外気間の熱伝達・熱放射量Ｑ3の算出］
次に、熱伝達・熱放射量算出部６は、第２測定部２によって測定された測定対象者近傍の温度Ｔa［ｔ］［℃］と、温度算出部１１がΔｔ前に算出した、時刻ｔにおける時刻ｔにおける測定対象者の平均皮膚温の推定値Ｔsk［ｔ］［℃］と体幹部皮膚層ＵＳの温度の推定値ＴUS［ｔ］［℃］と四肢部皮膚層ＬＳの温度の推定値ＴLS［ｔ］［℃］とに基づいて、測定対象者の体幹部Ｕにおける皮膚と外気間の熱伝達・熱放射量Ｑ3,U［Ｗ］、四肢部Ｌにおける皮膚と外気間の熱伝達・熱放射量Ｑ3,L［Ｗ］を算出する（図３ステップＳ１０５）。第２測定部２が測定対象者の衣服外の気温Ｔa［ｔ］［℃］を測定する場合、熱伝達・熱放射量Ｑ3,U，Ｑ3,Lは式（１６）、式（１７）のようになる。
Figure 0007496551000016
Figure 0007496551000017
割合ｓｆ＿ｃｏｎｆ＿ＵＳ，ｓｆ＿ｃｏｎｆ＿ＬＳについては上記のとおりである。ｆｃｌ＿ＵＳは体幹部Ｕの衣服による伝熱効率を表す定数、ｆｃｌ＿ＬＳは四肢部Ｌの衣服による伝熱効率を表す定数、ｃｏｖｅｒａｇｅは四肢部Ｌを衣服が覆う割合［％］である。割合ｓｆ＿ｃｏｎｆ＿ＵＳ，ｓｆ＿ｃｏｎｆ＿ＬＳ、定数ｆｃｌ＿ＵＳ，ｆｃｌ＿ＬＳ、割合ｃｏｖｅｒａｇｅは、それぞれ既知の値であり、実際に即した値を用いればよい。定数ｆｃｌ＿ＵＳ，ｆｃｌ＿ＬＳ、割合ｃｏｖｅｒａｇｅについては、測定対象者の服装に応じて適宜設定すればよい。
ＨＳは測定対象者の空気との熱交換係数［Ｗ／℃］である。熱伝達・熱放射量算出部６は、熱交換係数ＨＳ［Ｗ／℃］を式（１８）により算出する。
Figure 0007496551000018
ｓｆは測定対象者の体表面積［ｍ2］である。熱伝達・熱放射量算出部６は、体表面積ｓｆ［ｍ2］を測定対象者の重量ｗｅｉｇｈｔ［ｋｇ］と身長ｈｅｉｇｈｔ［ｍ］とから算出することができる。体表面積ｓｆ［ｍ2］の算出式としては、ＤｕＢｏｉｓの式や藤本の式等の推定式がある。
Ｈcmは対流熱伝達係数［Ｗ／℃／ｍ2］、Ｈrは放射熱伝達係数［Ｗ／℃／ｍ2］である。熱伝達・熱放射量算出部６は、対流熱伝達係数Ｈcm［Ｗ／℃／ｍ2］、放射熱伝達係数Ｈr［Ｗ／℃／ｍ2］を、それぞれ式（１９）、式（２０）により算出する。
Figure 0007496551000019
Figure 0007496551000020
式（１９）、式（２０）は、文献「D.Fiala，et al.，“A computer model of human thermoregulation for a wide range of environmental conditions：the passive system”，Journal of Applied Physiology，1985」に開示されている。Ｖairは風速［ｍ／ｓ］である。熱伝達・熱放射量算出部６は、風速Ｖair［ｍ／ｓ］として、風速計等で測定された実測値を用いてもよいし、文献等で開示されている衣服内の風速の既知の値を用いてもよい。
熱伝達・熱放射量算出部６は、第２測定部２が測定対象者の衣服内の温度Ｔa［ｔ］［℃］を測定する場合、熱伝達・熱放射量Ｑ3,U，Ｑ3,Lを式（２１）、式（２２）により算出する。
Figure 0007496551000021
Figure 0007496551000022
［皮膚蒸散量Ｑ4の算出］
次に、皮膚蒸散量算出部７は、第２測定部２によって測定された測定対象者近傍の温度Ｔa［ｔ］［℃］と、第３測定部３によって測定された測定対象者近傍の相対湿度ｈｕｍｉｄｉｔｙ［ｔ］と、温度算出部１１がΔｔ前に算出した、時刻ｔにおける平均皮膚温の推定値Ｔsk［ｔ］［℃］と深部体温の推定値Ｔ［ｔ］［℃］とに基づいて、測定対象者の体幹部Ｕにおける皮膚蒸散量Ｑ4,U［Ｗ］、四肢部Ｌにおける皮膚蒸散量Ｑ4,L［Ｗ］をそれぞれ式（２３）、式（２４）により算出する（図３ステップＳ１０６）。
Figure 0007496551000023
Figure 0007496551000024
ｓｗ＿ｃｏｎｆ＿ＵＳは測定対象者の体表面全体に対する体幹部Ｕの発汗量の割合［％］、ｓｗ＿ｃｏｎｆ＿ＬＳは体表面全体に対する四肢部Ｌの発汗量の割合［％］である。割合ｓｗ＿ｃｏｎｆ＿ＵＳ，ｓｗ＿ｃｏｎｆ＿ＬＳは、それぞれ既知の値であり、実際に即した値を用いればよく、計算の実施において適宜設定すればよい。ｓｗは全身の皮膚蒸散量［Ｗ］である。皮膚蒸散量算出部７は、皮膚蒸散量ｓｗ［Ｗ］を式（２５）により算出することができる。
Figure 0007496551000025
ｍｉｎ（Ｅ，Ｅmax）はＥとＥmaxのうちいずれか小さい方を採用することを意味する。Ｅは、測定対象者の皮膚における不感蒸散と有感蒸散の和［Ｗ］である。皮膚蒸散量算出部７は、不感蒸散と有感蒸散の和Ｅ［Ｗ］を式（２６）により算出することができる。
Figure 0007496551000026
ＰＩは測定対象者の皮膚における不感蒸散［Ｗ］、Ｑevは水の蒸発熱［Ｊ／ｇ］である。不感蒸散ＰＩ［Ｗ］、蒸発熱Ｑev［Ｊ／ｇ］は、それぞれ既知の値であり、実際に即した値を用いればよく、計算の実施において適宜設定すればよい。ｓｗｒａｔｅは有感蒸散［ｇ／ｍｉｎ］である。皮膚蒸散量算出部７は、有感蒸散ｓｗｒａｔｅ［ｇ／ｍｉｎ］を式（２７）により算出することができる。
Figure 0007496551000027
式（２７）は、文献「D.Fiala，et al.，“Computer prediction of human thermoregulatory and temperature responses to a wide range of environmental conditions”，International Journal of Biometeorology，volume 45，2001」に開示されている。ａｉｊ，ｂｉｊ（ｉ＝１，２，ｊ＝０，１）は、発汗係数である。発汗係数ａｉｊ，ｂｉｊは、それぞれ既知の値であり、測定対象者の汗のかき易さに応じて実際に即した値を用いればよく、計算の実施において適宜設定すればよい。具体的には、ｌｏｗ（汗をかき難い）、ｎｏｒｍａｌ（普通）、ｈｉｇｈ（汗をかき易い）の３水準に応じて、式（２８）に示すように発汗係数ａｉｊ，ｂｉｊを設定すればよい。
Figure 0007496551000028
一方、Ｅmaxは最大蒸発熱［Ｗ］である。皮膚蒸散量算出部７は、第２測定部２が測定対象者の衣服外の気温Ｔa［ｔ］［℃］を測定し、第３測定部３が測定対象者の衣服外の相対湿度ｈｕｍｉｄｉｔｙ［ｔ］を測定する場合、最大蒸発熱Ｅmax［Ｗ］を式（２９）により算出することができる。
Figure 0007496551000029
また、皮膚蒸散量算出部７は、第２測定部２が測定対象者の衣服内の温度Ｔa［ｔ］［℃］を測定し、第３測定部３が測定対象者の衣服内の相対湿度ｈｕｍｉｄｉｔｙ［ｔ］を測定する場合、最大蒸発熱Ｅmax［Ｗ］を式（３０）により算出することができる。
Figure 0007496551000030
割合ｓｗ＿ｃｏｎｆ＿ＵＳ，ｓｗ＿ｃｏｎｆ＿ＬＳ，ｃｏｖｅｒａｇｅについては上記のとおりである。ｆｐｃｌ＿ＵＳは測定対象者の体幹部Ｕにおける衣服による伝熱効率を表す定数、ｆｐｃｌ＿ＬＳは四肢部Ｌにおける衣服による伝熱効率を表す定数、Ｅmax＿coefは最大蒸発熱に関する定数である。定数ｆｐｃｌ＿ＵＳ，ｆｐｃｌ＿ＬＳ，Ｅmax＿coefは、それぞれ既知の値であり、実際に即した値を用いればよく、計算の実施において適宜設定すればよい。
Ｈcは空気の風速に依存した対流による熱移動［Ｗ・ｍ2／℃］である。皮膚蒸散量算出部７は、熱移動Ｈc［Ｗ・ｍ2／℃］を式（３１）により算出することができる。
Figure 0007496551000031
式（２９）～式（３１）は、文献「I.Laakso,et al.，“Dominant factors affecting temperature rise in simulations of human thermoregulation during RF exposure”，Physics in Medicine and Biology，Volume 56，2011」に開示されている。Ｐsは測定対象者の皮膚層での飽和水蒸気圧［ｋＰａ］である。皮膚蒸散量算出部７は、飽和水蒸気圧Ｐs［ｋＰａ］を式（３２）により算出することができる。
Figure 0007496551000032
Ｐaは湿度を計測している雰囲気中での飽和水蒸気圧［ｋＰａ］である。皮膚蒸散量算出部７は、飽和水蒸気圧Ｐa［ｋＰａ］を式（３３）により算出することができる。
Figure 0007496551000033
［呼気蒸散量Ｑ5の決定］
次に、呼気蒸散量設定部８は、測定対象者の呼気による蒸散量Ｑ5［Ｗ］を式（３４）のように設定する（図３ステップＳ１０７）。呼気蒸散量Ｑ5［Ｗ］は、外気と測定対象者の深部間の熱交換量に相当する。
Figure 0007496551000034
［深部と皮膚間の熱交換量Ｑ6の算出］
次に、熱交換量算出部９は、第１測定部１によって測定された心拍数ＨＲ［ｔ］［ｂｐｍ］と、温度算出部１１がΔｔ前に算出した、時刻ｔにおける体幹部皮膚層ＵＳの温度の推定値ＴUS［ｔ］［℃］と体幹部深部層ＵＣの温度の推定値ＴUC［ｔ］［℃］と四肢部皮膚層ＬＳの温度の推定値ＴLS［ｔ］［℃］と四肢部深部層ＬＣの温度の推定値ＴLC［ｔ］［℃］と平均皮膚温の推定値Ｔsk［ｔ］［℃］と深部体温の推定値Ｔ［ｔ］［℃］とに基づいて、測定対象者の体幹部Ｕにおける深部と皮膚間の熱交換量Ｑ6,U［Ｗ］、四肢部Ｌにおける深部と皮膚間の熱交換量Ｑ6,L［Ｗ］をそれぞれ式（３５）、式（３６）により算出する（図３ステップＳ１０８）。
Figure 0007496551000035
Figure 0007496551000036
割合ｓｆ＿ｃｏｎｆ＿ＵＳ，ｓｆ＿ｃｏｎｆ＿ＬＳについては上記のとおりである。ｈｘは測定対象者の皮膚と深部間の熱交換係数である。熱交換量算出部９は、熱交換係数ｈｘを式（３７）により算出することができる。
Figure 0007496551000037
ＭＥＴｓ［ｔ］については上記のとおりである。ａ，ｂ，ｅ，ｈｘ０，ｈｘ１，ｈｘ＿ｍａｘは熱交換係数に関わるパラメータである。パラメータａ，ｂ，ｅ，ｈｘ０，ｈｘ１，ｈｘ＿ｍａｘは、それぞれ既知の値であり、実際に即した値を用いればよく、実験を通じて適宜設定すればよい。
［体幹部と四肢部間の熱交換量Ｑ7の算出］
次に、熱交換量算出部１０は、温度算出部１１がΔｔ前に算出した、時刻ｔにおける体幹部深部層ＵＣの温度の推定値ＴUC［ｔ］［℃］と四肢部深部層ＬＣの温度ＴLC［ｔ］［℃］と平均皮膚温の推定値Ｔsk［ｔ］［℃］と深部体温の推定値Ｔ［ｔ］［℃］とに基づいて、測定対象者の深部における体幹部Ｕと四肢部Ｌ間の熱交換量Ｑ7［Ｗ］を式（３８）により算出する（図３ステップＳ１０９）。
Figure 0007496551000038
ｈｃｃは測定対象者の深部における体幹部Ｕと四肢部Ｌの熱交換係数である。熱交換量算出部１０は、熱交換係数ｈｃｃを式（３９）により算出することができる。
Figure 0007496551000039
ｆ，ｈｃｃ０は熱交換係数に関わるパラメータである。パラメータｆ，ｈｃｃ０は、それぞれ既知の値であり、実際に即した値を用いればよく、実験を通じて適宜設定すればよい。ｈｃｃ＿Ｔは熱交換係数ｈｃｃの温度寄与分、ｈｃｃ＿Ｍは熱交換係数ｈｃｃのＭＥＴｓ寄与分である。熱交換量算出部１０は、熱交換係数ｈｃｃの温度寄与分ｈｃｃ＿Ｔを式（４０）により算出することができる。
Figure 0007496551000040
ｈｃｃ＿Ｔｍａｘはｈｃｃ＿Ｔの規定の上限値である。また、熱交換量算出部１０は、熱交換係数ｈｃｃのＭＥＴｓ寄与分ｈｃｃ＿Ｍを式（４１）により算出することができる。
Figure 0007496551000041
ｈｃｃ＿Ｍｍａｘはｈｃｃ＿Ｍの規定の上限値である。ｈｃｃ１は熱交換係数に関わるパラメータである。上限値ｈｃｃ＿Ｔｍａｘ，ｈｃｃ＿Ｍｍａｘ、パラメータｈｃｃ１は、それぞれ既知の値であり、実際に即した値を用いればよく、実験を通じて適宜設定すればよい。ａｖｅＭＥＴｓ［ｔ］は、メッツ数の時間平均値（例えば、６分間の平均値）である。パラメータａ，ｂについては上記のとおりである。
こうして、発熱量算出部４と代謝量算出部５と熱伝達・熱放射量算出部６と皮膚蒸散量算出部７と呼気蒸散量設定部８と熱交換量算出部９と熱交換量算出部１０とにより、各熱量Ｑ1,U，Ｑ1,L，Ｑ2,US，Ｑ2,UC，Ｑ2,LS，Ｑ2,LC，Ｑ3,U，Ｑ3,L，Ｑ4,U，Ｑ4,L，Ｑ5，Ｑ6,U，Ｑ6,L，Ｑ7を算出することができる。
温度算出部１１は、時刻ｔにおける測定対象者の体幹部皮膚層ＵＳの温度の推定値ＴUS［ｔ］［℃］と、体幹部皮膚層ＵＳの代謝量Ｑ2,US［Ｗ］と、体幹部Ｕにおける皮膚と外気間の熱伝達・熱放射量Ｑ3,U［Ｗ］と、体幹部Ｕにおける皮膚蒸散量Ｑ4,U［Ｗ］と、体幹部Ｕにおける深部と皮膚間の熱交換量Ｑ6,U［Ｗ］とに基づいて、Δｔ後の体幹部皮膚層ＵＳの温度の推定値ＴUS［ｔ＋Δｔ］［℃］を式（１）により算出する（図３ステップＳ１１０）。
温度算出部１１は、時刻ｔにおける測定対象者の体幹部深部層ＵＣの温度の推定値ＴUC［ｔ］［℃］と、体幹部Ｕの深部層における発熱量Ｑ1,U［Ｗ］と、体幹部深部層ＵＣの代謝量Ｑ2,UC［Ｗ］と、呼気蒸散量Ｑ5［Ｗ］と、体幹部Ｕにおける深部と皮膚間の熱交換量Ｑ6,U［Ｗ］と、深部における体幹部Ｕと四肢部Ｌ間の熱交換量Ｑ7［Ｗ］とに基づいて、Δｔ後の体幹部深部層ＵＣの温度の推定値ＴUC［ｔ＋Δｔ］［℃］を式（２）により算出する（図３ステップＳ１１１）。
温度算出部１１は、時刻ｔにおける測定対象者の四肢部皮膚層ＬＳの温度の推定値ＴLS［ｔ］［℃］と、四肢部皮膚層ＬＳの代謝量Ｑ2,LS［Ｗ］と、四肢部Ｌにおける皮膚と外気間の熱伝達・熱放射量Ｑ3,L［Ｗ］と、四肢部Ｌにおける皮膚蒸散量Ｑ4,L［Ｗ］と、四肢部Ｌにおける深部と皮膚間の熱交換量Ｑ6,L［Ｗ］とに基づいて、Δｔ後の四肢部皮膚層ＬＳの温度の推定値ＴLS［ｔ＋Δｔ］［℃］を式（３）により算出する（図３ステップＳ１１２）。
温度算出部１１は、時刻ｔにおける測定対象者の四肢部深部層ＬＣの温度の推定値ＴLC［ｔ］［℃］と、四肢部Ｌの深部層における発熱量Ｑ1,L［Ｗ］と、四肢部深部層ＬＣの代謝量Ｑ2,LC［Ｗ］と、四肢部Ｌにおける深部と皮膚間の熱交換量Ｑ6,L［Ｗ］と、深部における体幹部Ｕと四肢部Ｌ間の熱交換量Ｑ7［Ｗ］とに基づいて、Δｔ後の四肢部深部層ＬＣの温度の推定値ＴLC［ｔ＋Δｔ］［℃］を式（４）により算出する（図３ステップＳ１１３）。
こうして、温度の推定値ＴUS［ｔ＋Δｔ］，ＴUC［ｔ＋Δｔ］，ＴLS［ｔ＋Δｔ］，ＴLC［ｔ＋Δｔ］をそれぞれ逐次計算することができる。
そして、温度算出部１１は、Δｔ後の平均皮膚温の推定値Ｔsk［ｔ＋Δｔ］［℃］を式（５）により算出する（図３ステップＳ１１４）。また、温度算出部１１は、Δｔ後の深部体温の推定値Ｔ［ｔ＋Δｔ］［℃］を式（６）により算出する（図３ステップＳ１１５）。
深部体温推定装置は、以上のステップＳ１００～Ｓ１１５の処理を周期Δｔ毎に行う。Δｔ後の次の計算では、直前の回で算出したＴUS［ｔ＋Δｔ］，ＴUC［ｔ＋Δｔ］，ＴLS［ｔ＋Δｔ］，ＴLC［ｔ＋Δｔ］，Ｔsk［ｔ＋Δｔ］，Ｔ［ｔ＋Δｔ］を、それぞれＴUS［ｔ］，ＴUC［ｔ］，ＴLS［ｔ］，ＴLC［ｔ］，Ｔsk［ｔ］，Ｔ［ｔ］として、ステップＳ１００～Ｓ１１５の処理を行うようにすればよい。
本実施例では、測定対象者の心拍数と測定対象者近傍の温度だけでなく、測定対象者近傍の湿度を測定して、測定対象者の各部位・各層に流入出する熱量を算出し、算出した熱量から測定対象者の深部体温を推定するので、従来よりも高精度に深部体温を推定することができる。
なお、本実施例では、測定対象者の心拍数を用いる場合について説明しているが、心拍数の代わりに脈拍数を用いてもよい。
本実施例の発熱量算出部４と代謝量算出部５と熱伝達・熱放射量算出部６と皮膚蒸散量算出部７と呼気蒸散量設定部８と熱交換量算出部９と熱交換量算出部１０と温度算出部１１とは、ＣＰＵ（Central Processing Unit）、記憶装置及びインタフェースを備えたコンピュータと、これらのハードウェア資源を制御するプログラムによって実現することができる。このコンピュータの構成例を図４に示す。
コンピュータは、ＣＰＵ２００と、記憶装置２０１と、インタフェース装置（Ｉ／Ｆ）２０２とを備えている。Ｉ／Ｆ２０２には、第１測定部１と第２測定部２と第３測定部３などが接続される。このようなコンピュータにおいて、本発明の深部体温推定方法を実現させるための深部体温推定プログラムは記憶装置２０１に格納される。ＣＰＵ２００は、記憶装置２０１に格納されたプログラムに従って本実施例で説明した処理を実行する。プログラムをネットワークを通して提供することも可能である。
また、発熱量算出部４と代謝量算出部５と熱伝達・熱放射量算出部６と皮膚蒸散量算出部７と呼気蒸散量設定部８と熱交換量算出部９と熱交換量算出部１０と温度算出部１１とを、複数のコンピュータ機器に分散させるようにしてもよい。
本発明は、人体の深部体温を推定する技術に適用することができる。
１… 第１測定部、２…第２測定部、３…第３測定部、４…発熱量算出部、５…代謝量算出部、６…熱伝達・熱放射量算出部、７…皮膚蒸散量算出部、８…呼気蒸散量設定部、９…熱交換量算出部、１０…熱交換量算出部、１１…温度算出部、１２…熱量算出部。

Claims (7)

  1. 測定対象者の心拍数または脈拍数を測定するように構成された第１測定部と、
    前記測定対象者の近傍の温度を測定するように構成された第２測定部と、
    前記測定対象者の近傍の湿度を測定するように構成された第３測定部と、
    前記測定対象者の心拍数または脈拍数と前記測定対象者の近傍の温度と前記測定対象者の近傍の湿度とに基づいて、前記測定対象者の第１部位の深部に流入出する熱量と前記第１部位の皮膚に流入出する熱量と前記測定対象者の第２部位の深部に流入出する熱量と前記第２部位の皮膚に流入出する熱量とを算出するように構成された熱量算出部と、
    前記熱量に基づいて、前記測定対象者の皮膚温度および深部温度を算出するように構成された温度算出部とを備え、
    前記第１部位は前記測定対象者の体幹部で、前記第２部位は前記測定対象者の四肢部、または前記第１部位は前記測定対象者の上半身で、前記第２部位は前記測定対象者の下半身であり、
    前記熱量算出部は、
    前記測定対象者の心拍数または脈拍数に基づいて、前記測定対象者の運動による前記第１部位と前記第２部位のそれぞれの発熱量を算出するように構成された発熱量算出部と、
    前記測定対象者の近傍の温度と前記温度算出部によって直前に算出された皮膚温度とに基づいて、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける皮膚と外気間の第１熱交換量を算出するように構成された第１熱交換量算出部と、
    前記測定対象者の心拍数または脈拍数と前記温度算出部によって直前に算出された皮膚温度と深部温度とに基づいて、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける深部と皮膚間の第２熱交換量を算出するように構成された第２熱交換量算出部と、
    前記温度算出部によって直前に算出された皮膚温度と深部温度とに基づいて、前記測定対象者の深部における前記第１部位と前記第２部位間の第３熱交換量を算出するように構成された第３熱交換量算出部とを含むことを特徴とする深部体温推定装置。
  2. 請求項１記載の深部体温推定装置において、
    前記第１熱交換量算出部は、前記測定対象者の皮膚が衣服で覆われている割合を考慮して、前記第１熱交換量を算出することを特徴とする深部体温推定装置。
  3. 請求項１または２記載の深部体温推定装置において、
    前記熱量算出部は、前記測定対象者の近傍の温度と前記測定対象者の近傍の湿度と前記温度算出部によって直前に算出された皮膚温度と深部温度とに基づいて、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける皮膚蒸散量を算出するように構成された皮膚蒸散量算出部をさらに含むことを特徴とする深部体温推定装置。
  4. 請求項３記載の深部体温推定装置において、
    前記熱量算出部は、外気と前記測定対象者の深部間の第４熱交換量を設定するように構成された第４熱交換量設定部と、
    前記温度算出部によって直前に算出された皮膚温度と深部温度とに基づいて、前記測定対象者の第１部位皮膚層の代謝量と第１部位深部層の代謝量と第２部位皮膚層の代謝量と第２部位深部層の代謝量を算出するように構成された代謝量算出部とをさらに含むことを特徴とする深部体温推定装置。
  5. 請求項４記載の深部体温推定装置において、
    前記温度算出部は、
    前記測定対象者の前記第１部位と前記第２部位の発熱量Ｑ１，Ｕ，Ｑ１，Ｌ、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける皮膚と外気間の第１熱交換量Ｑ３，Ｕ，Ｑ３，Ｌ、 前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける深部と皮膚間の第２熱交換量Ｑ６，Ｕ，Ｑ６，Ｌ、前記測定対象者の深部における前記第１部位と前記第２部位間の第３熱交換量Ｑ７、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける皮膚蒸散量Ｑ４，Ｕ，Ｑ４，Ｌ、外気と前記測定対象者の深部間の第４熱交換量Ｑ５、前記測定対象者の第１部位皮膚層の代謝量Ｑ２，ＵＳ、第１部位深部層の代謝量Ｑ２，ＵＣ、第２部位皮膚層の代謝量Ｑ２，ＬＳ、第２部位深部層の代謝量Ｑ２，ＬＣ、前記測定対象者の第１部位皮膚層の熱容量ＷＣＵＳ、第１部位深部層の熱容量ＷＣＵＣ、第２部位皮膚層の熱容量ＷＣＬＳ、第２部位深部層の熱容量ＷＣＬＣ、直前に算出した第１部位皮膚層の温度の推定値ＴＵＳ［ｔ］、第１部位深部層の温度の推定値ＴＵＣ［ｔ］、第２部位皮膚層の温度の推定値ＴＬＳ［ｔ］、第２部位深部層の温度の推定値ＴＬＣ［ｔ］に基づいて、
    Figure 0007496551000042
     により、時刻（ｔ＋Δｔ）における第１部位皮膚層の温度の推定値ＴＵＳ［ｔ＋Δｔ］、第１部位深部層の温度の推定値ＴＵＣ［ｔ＋Δｔ］、第２部位皮膚層の温度の推定値ＴＬＳ［ｔ＋Δｔ］、第２部位深部層の温度の推定値ＴＬＣ［ｔ＋Δｔ］を算出することを特徴とする深部体温推定装置。
  6. 測定対象者の心拍数または脈拍数を測定する第１ステップと、
    前記測定対象者の近傍の温度を測定する第２ステップと、
    前記測定対象者の近傍の湿度を測定する第３ステップと、
    前記測定対象者の心拍数または脈拍数と前記測定対象者の近傍の温度と前記測定対象者の近傍の湿度とに基づいて、前記測定対象者の第１部位の深部に流入出する熱量と前記第１部位の皮膚に流入出する熱量と前記測定対象者の第２部位の深部に流入出する熱量と前記第２部位の皮膚に流入出する熱量とを算出する第４ステップと、
    前記熱量に基づいて、前記測定対象者の皮膚温度および深部温度を算出する第５ステップとを含み、
    前記第１部位は前記測定対象者の体幹部で、前記第２部位は前記測定対象者の四肢部、または前記第１部位は前記測定対象者の上半身で、前記第２部位は前記測定対象者の下半身であり、
    前記第４ステップは、
    前記測定対象者の心拍数または脈拍数に基づいて、前記測定対象者の運動による前記第１部位と前記第２部位のそれぞれの発熱量を算出するステップと、
    前記測定対象者の近傍の温度と前記第５ステップによって直前に算出された皮膚温度とに基づいて、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける皮膚と外気間の第１熱交換量を算出するステップと、
    前記測定対象者の心拍数または脈拍数と前記第５ステップによって直前に算出された皮膚温度と深部温度とに基づいて、前記測定対象者の前記第１部位と前記第２部位のそれぞれにおける深部と皮膚間の第２熱交換量を算出するステップと、
    前記第５ステップによって直前に算出された皮膚温度と深部温度とに基づいて、前記測定対象者の深部における前記第１部位と前記第２部位間の第３熱交換量を算出するステップとを含むことを特徴とする深部体温推定方法。
  7. 請求項６記載の各ステップをコンピュータに実行させることを特徴とする深部体温推定プログラム。
JP2021007047A 2021-01-20 2021-01-20 深部体温推定装置、深部体温推定方法および深部体温推定プログラム Active JP7496551B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021007047A JP7496551B2 (ja) 2021-01-20 2021-01-20 深部体温推定装置、深部体温推定方法および深部体温推定プログラム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2021007047A JP7496551B2 (ja) 2021-01-20 2021-01-20 深部体温推定装置、深部体温推定方法および深部体温推定プログラム

Publications (2)

Publication Number Publication Date
JP2022111550A JP2022111550A (ja) 2022-08-01
JP7496551B2 true JP7496551B2 (ja) 2024-06-07

Family

ID=82655505

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021007047A Active JP7496551B2 (ja) 2021-01-20 2021-01-20 深部体温推定装置、深部体温推定方法および深部体温推定プログラム

Country Status (1)

Country Link
JP (1) JP7496551B2 (ja)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140180027A1 (en) 2012-12-20 2014-06-26 U.S. Government, As Represented By The Secretary Of The Army Estimation of Human Core Temperature based on Heart Rate System and Method
JP2017217224A (ja) 2016-06-08 2017-12-14 国立大学法人大阪大学 深部体温推定装置、その方法及びプログラム
JP2020065823A (ja) 2018-10-26 2020-04-30 日本電信電話株式会社 深部体温推定装置
JP2020162952A (ja) 2019-03-29 2020-10-08 学校法人産業医科大学 深部体温推定装置および方法並びにプログラム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140180027A1 (en) 2012-12-20 2014-06-26 U.S. Government, As Represented By The Secretary Of The Army Estimation of Human Core Temperature based on Heart Rate System and Method
JP2017217224A (ja) 2016-06-08 2017-12-14 国立大学法人大阪大学 深部体温推定装置、その方法及びプログラム
JP2020065823A (ja) 2018-10-26 2020-04-30 日本電信電話株式会社 深部体温推定装置
JP2020162952A (ja) 2019-03-29 2020-10-08 学校法人産業医科大学 深部体温推定装置および方法並びにプログラム

Also Published As

Publication number Publication date
JP2022111550A (ja) 2022-08-01

Similar Documents

Publication Publication Date Title
Tanda Skin temperature measurements by infrared thermography during running exercise
Niedermann et al. Prediction of human core body temperature using non-invasive measurement methods
WO2017213011A1 (ja) 深部体温推定装置、その方法及びプログラム
EP2070470B1 (en) Electronic device, arrangement, and method of estimating fluid loss
US9597038B2 (en) Apparatus for monitoring a subject during exercise
CN110175732B (zh) 一种低温作业环境下人体冷应激精细分级评估系统
Hirata et al. Computational estimation of body temperature and sweating in the aged during passive heat exposure
US20180153463A1 (en) Wrist device and arrangement for measuring and testing performance
JP7026905B2 (ja) 深部体温推定装置
Vargas et al. Thermal behavior alleviates thermal discomfort during steady-state exercise without affecting whole body heat loss
Wu et al. A multi-segmented human bioheat model for cold and extremely cold exposures
Fakir et al. Prediction of individual thermal sensation from exhaled breath temperature using a smart face mask
Hirata et al. Body core temperature estimation using new compartment model with vital data from wearable devices
JP7496551B2 (ja) 深部体温推定装置、深部体温推定方法および深部体温推定プログラム
Wolff et al. Heat exchanges in wet suits
CN106650262B (zh) 基于个性化热调节模型的实时内核温度估计方法及系统
JP7572679B2 (ja) 深部体温推定装置、深部体温推定方法および深部体温推定プログラム
Kamiya et al. Different thermoregulatory responses of people from tropical and temperate zones: A computational study
Cadarette et al. Cross validation of USARIEM heat strain prediction models
JP7576265B2 (ja) 深部体温推定装置、深部体温推定方法および深部体温推定プログラム
Veselá et al. Effect of local skin blood flow during light and medium activities on local skin temperature predictions
CN116942092A (zh) 基于多节点睡眠热调节人体模型的皮肤温度预测方法
JP2022111549A (ja) 深部体温推定装置、深部体温推定方法および深部体温推定プログラム
JP2022111548A (ja) 深部体温推定装置、深部体温推定方法および深部体温推定プログラム
JP7312410B2 (ja) 深部体温推定装置

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20210121

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20230512

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20231031

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20231212

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240213

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20240514

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240517

R150 Certificate of patent or registration of utility model

Ref document number: 7496551

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150