JP2022111548A - Deep part body temperature estimation device, deep part body temperature estimation method, and deep part body temperature estimation program - Google Patents

Deep part body temperature estimation device, deep part body temperature estimation method, and deep part body temperature estimation program Download PDF

Info

Publication number
JP2022111548A
JP2022111548A JP2021007044A JP2021007044A JP2022111548A JP 2022111548 A JP2022111548 A JP 2022111548A JP 2021007044 A JP2021007044 A JP 2021007044A JP 2021007044 A JP2021007044 A JP 2021007044A JP 2022111548 A JP2022111548 A JP 2022111548A
Authority
JP
Japan
Prior art keywords
temperature
deep
measurement subject
trunk
heat exchange
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.)
Pending
Application number
JP2021007044A
Other languages
Japanese (ja)
Inventor
優生 橋本
Yuki Hashimoto
啓 桑原
Hiroshi Kuwabara
浩芳 都甲
Hiroyoshi Toko
和彦 高河原
Kazuhiko Takagahara
隆子 石原
Takako Ishihara
晃正 平田
Akimasa Hirata
涼太 上松
Ryota Uematsu
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 JP2021007044A priority Critical patent/JP2022111548A/en
Publication of JP2022111548A publication Critical patent/JP2022111548A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

To accurately estimate a deep part body temperature of a person to be measured.SOLUTION: A deep part body temperature estimation device includes: a first measuring unit 1 for measuring a heart rate or a pulse rate of a person to be measured; a second measuring unit 2 for measuring a skin temperature of the person to be measured; an amount-of-heat calculation unit 9 for calculating an amount of heat that flows into and outflows from deep parts of a first body part and a second body part, and the skin of the person to be measured on the basis of the heart rate or the pulse rate and the skin temperature of the person to be measured; and a temperature calculation unit 8 for calculating a deep part temperature of the person to be measured on the basis of the amount of heat.SELECTED DRAWING: Figure 1

Description

本発明は、測定対象者の深部体温を推定する深部体温推定装置、深部体温推定方法および深部体温推定プログラムに関するものである。 The present invention relates to a core body temperature estimation device, a core body temperature estimation method, and a core body temperature estimation program for estimating the core body temperature of a person to be measured.

近年、熱中症により多くの人が救急搬送されており、2019年5月から9月の全国における熱中症による救急搬送者数の累計は、71317人であった(非特許文献1参照)。熱中症の診療や熱中症のリスク軽減のために、深部体温の測定が有効であるといわれている(非特許文献2参照)。 In recent years, many people have been transported by ambulance due to heat stroke, and the total number of people transported by ambulance due to heat stroke nationwide from May to September 2019 was 71,317 (see Non-Patent Document 1). It is said that measurement of core body temperature is effective for diagnosis of heat stroke and reduction of the risk of heat stroke (see Non-Patent Document 2).

深部体温を監視するためには、直腸温、膀胱温、食道温などを測定する必要があり、日常生活のなかでこれらの測定を実施することは難しい。そこで、深部体温を推定する方法として、測定対象者の心拍数、雰囲気の気温から深部温度を簡易的に推定する方法が提案されている(特許文献1参照)。 In order to monitor core body temperature, it is necessary to measure rectal temperature, bladder temperature, esophageal temperature, etc., and it is difficult to carry out these measurements in daily life. Therefore, as a method of estimating the core body temperature, a method of simply estimating the core body temperature from the heart rate of the person to be measured and the air temperature of the atmosphere has been proposed (see Patent Document 1).

しかしながら、特許文献1に開示された方法では、深部体温の推定精度が低いという課題があった。 However, the method disclosed in Patent Literature 1 has a problem that the accuracy of estimating core body temperature is low.

特開2020-65823号公報JP 2020-65823 A

“2019年(5月から9月)の熱中症による救急搬送状況”,総務省,[2020年3月31日検索],<https://www.fdma.go.jp/disaster/heatstroke/items/heatstroke_geppou_2019.pdf>“2019 (May to September) emergency transportation due to heatstroke”, Ministry of Internal Affairs and Communications, [searched March 31, 2020], <https://www.fdma.go.jp/disaster/heatstroke/items /heatstroke_geppou_2019.pdf> “熱中症診療ガイドライン 2015”,一般社団法人日本救急医学会,2015年,[2019年3月31日検索],<https://www.mhlw.go.jp/file/06-Seisakujouhou-10800000-Iseikyoku/heatstroke2015.pdf>“Heat Stroke Treatment Guidelines 2015”, Japanese Association for Acute Medicine, 2015, [searched March 31, 2019], <https://www.mhlw.go.jp/file/06-Seisakujouhou-10800000- Iseikyoku/heatstroke2015.pdf>

本発明は、上記課題を解決するためになされたもので、測定対象者の深部体温を従来よりも高精度に推定することができる深部体温推定装置、深部体温推定方法および深部体温推定プログラムを提供することを目的とする。 The present invention has been made to solve the above problems, and provides a core body temperature estimation device, a core body temperature estimation method, and a core body temperature estimation program capable of estimating the core body temperature of a person to be measured with higher accuracy than ever before. intended to

本発明の深部体温推定装置は、測定対象者の心拍数または脈拍数を測定するように構成された第1測定部と、前記測定対象者の皮膚温度を測定するように構成された第2測定部と、前記測定対象者の心拍数または脈拍数と皮膚温度とに基づいて、前記測定対象者の第1部位と第2部位の深部と皮膚とに流入出する熱量を算出するように構成された熱量算出部と、前記熱量に基づいて、前記測定対象者の深部温度を算出するように構成された温度算出部とを備えることを特徴とするものである。 A core body temperature estimating device of the present invention includes a first measuring unit configured to measure the heart rate or pulse rate of a measurement subject, and a second measurement unit configured to measure the skin temperature of the measurement subject. and the heart rate or pulse rate of the measurement subject and the skin temperature of the measurement subject. and a temperature calculation unit configured to calculate the deep temperature of the person to be measured based on the heat quantity.

また、本発明の深部体温推定装置の1構成例において、前記第1部位は、前記測定対象者の体幹部であり、前記第2部位は、前記測定対象者の四肢部であり、前記熱量算出部は、前記測定対象者の心拍数または脈拍数に基づいて、前記測定対象者の運動による前記体幹部と前記四肢部のそれぞれの発熱量を算出するように構成された発熱量算出部と、前記測定対象者の心拍数または脈拍数と皮膚温度と前記温度算出部によって直前に算出された深部温度とに基づいて、前記測定対象者の前記体幹部と前記四肢部のそれぞれにおける深部と皮膚間の第1熱交換量を算出するように構成された第1熱交換量算出部と、前記測定対象者の心拍数または脈拍数と皮膚温度と前記温度算出部によって直前に算出された深部温度とに基づいて、前記測定対象者の深部における前記体幹部と前記四肢部間の第2熱交換量を算出するように構成された第2熱交換量算出部と、前記温度算出部によって直前に算出された深部温度に基づいて、前記測定対象者の体幹部深部層の代謝量と四肢部深部層の代謝量を算出するように構成された代謝量算出部とを含むことを特徴とするものである。 Further, in one configuration example of the core body temperature estimating device of the present invention, the first part is the trunk of the measurement subject, the second part is the extremities of the measurement subject, and the calorie calculation is performed. a calorific value calculation unit configured to calculate calorific value of each of the trunk and extremities due to exercise of the measurement subject based on the heart rate or pulse rate of the measurement subject; Based on the heart rate or pulse rate of the measurement subject, the skin temperature, and the deep temperature calculated immediately before by the temperature calculation unit, the distance between the deep part and the skin in each of the trunk and the extremities of the measurement subject a first heat exchange amount calculation unit configured to calculate the first heat exchange amount of the subject, the heart rate or pulse rate of the measurement subject, the skin temperature, and the deep temperature calculated immediately before by the temperature calculation unit calculated immediately before by a second heat exchange amount calculation unit configured to calculate a second heat exchange amount between the trunk and the extremities in the deep part of the measurement subject based on and a metabolic rate calculation unit configured to calculate the metabolic rates of the deep layers of the trunk and the deep layers of the extremities of the person to be measured based on the obtained deep temperature. be.

また、本発明の深部体温推定装置の1構成例において、前記第2測定部は、前記体幹部と前記四肢部のそれぞれの皮膚温度を測定し、前記第1熱交換量算出部は、前記第2測定部の測定結果を用いて前記測定対象者の前記体幹部と前記四肢部のそれぞれにおける深部と皮膚間の第1熱交換量を算出し、前記第2熱交換量算出部は、前記第2測定部の測定結果を用いて前記測定対象者の深部における前記体幹部と前記四肢部間の第2熱交換量を算出することを特徴とするものである。
また、本発明の深部体温推定装置の1構成例において、前記熱量算出部は、外気と前記測定対象者の深部間の第3熱交換量を設定するように構成された第3熱交換量設定部をさらに含むことを特徴とするものである。
Further, in one configuration example of the core body temperature estimating device of the present invention, the second measuring unit measures the skin temperature of each of the trunk and the extremities, and the first heat exchange amount calculating unit measures the skin temperature of the first 2 Calculate a first heat exchange amount between the deep part and the skin in each of the trunk and the extremities of the measurement subject using the measurement results of the measurement unit, and the second heat exchange amount calculation unit The method is characterized in that a second heat exchange amount between the trunk and the extremities in the deep part of the person to be measured is calculated using the measurement results of two measurement units.
Further, in one configuration example of the core body temperature estimating device of the present invention, the heat amount calculation unit is configured to set a third heat exchange amount between outside air and the deep part of the measurement subject. It is characterized by further comprising a part.

また、本発明の深部体温推定装置の1構成例において、前記温度算出部は、前記測定対象者の前記体幹部と前記四肢部の発熱量Q1,U,Q1,L、前記測定対象者の前記体幹部と前記四肢部のそれぞれにおける深部と皮膚間の第1熱交換量Q4,U,Q4,L、前記測定対象者の深部における前記体幹部と前記四肢部間の第2熱交換量Q5、外気と前記測定対象者の深部間の第3熱交換量Q3、前記測定対象者の体幹部深部層の代謝量Q2,U、四肢部深部層の代謝量Q2,L、前記測定対象者の体幹部深部層の熱容量WCUC、四肢部深部層の熱容量WCLC、直前に算出した体幹部深部層の温度の推定値TUC[t]、四肢部深部層の温度の推定値TLC[t]に基づいて、時刻(t+Δt)における体幹部深部層の温度の推定値TUC[t+Δt]、四肢部深部層の温度の推定値TLC[t+Δt]を算出することを特徴とするものである。 In one configuration example of the apparatus for estimating core body temperature according to the present invention, the temperature calculator may calculate heat generation amounts Q 1,U and Q 1,L of the trunk and extremities of the measurement subject, First heat exchange amounts Q 4,U , Q 4,L between the deep part and the skin in the trunk and the extremities, respectively, and the second heat between the trunk and the extremities in the deep part of the measurement subject exchange amount Q 5 , third heat exchange amount Q 3 between the outside air and the deep part of the subject, metabolic rate Q 2,U in the deep layer of the trunk of the subject, metabolic rate Q 2 , in the deep layer of the extremities L , heat capacity WC UC of the deep layer of the trunk of the subject, heat capacity WC LC of the deep layer of the extremities, estimated value T UC [t] of the temperature of the deep part of the trunk calculated immediately before, temperature of the deep layers of the extremities Based on the estimated value T LC [t] of , calculate the estimated value T UC [t + Δt] of the temperature of the deep layer of the trunk and the estimated value T LC [t + Δt] of the temperature of the deep layer of the extremities at the time (t + Δt) It is characterized by

また、本発明の深部体温推定方法は、測定対象者の心拍数または脈拍数を測定する第1ステップと、前記測定対象者の皮膚温度を測定する第2ステップと、前記測定対象者の心拍数または脈拍数と皮膚温度とに基づいて、前記測定対象者の第1部位と第2部位の深部と皮膚とに流入出する熱量を算出する第3ステップと、前記熱量に基づいて、前記測定対象者の深部温度を算出する第4ステップとを含むことを特徴とするものである。
また、本発明の深部体温推定プログラムは、前記の各ステップをコンピュータに実行させることを特徴とするものである。
Further, the core body temperature estimation method of the present invention includes a first step of measuring the heart rate or pulse rate of a measurement subject, a second step of measuring the skin temperature of the measurement subject, and a heart rate of the measurement subject. Alternatively, a third step of calculating the amount of heat flowing into and out of the deep part and the skin of the first and second parts of the measurement subject based on the pulse rate and the skin temperature; and a fourth step of calculating the deep temperature of the person.
Further, a core body temperature estimating program of the present invention is characterized by causing a computer to execute each of the steps described above.

本発明によれば、測定対象者の心拍数だけでなく、測定対象者の皮膚温度を測定して、測定対象者の第1部位と第2部位の深部と皮膚とに流入出する熱量を算出し、算出した熱量から測定対象者の深部体温を推定するので、従来よりも高精度に深部体温を推定することができる。 According to the present invention, not only the heart rate of the person to be measured but also the skin temperature of the person to be measured is measured, and the amount of heat flowing into and out of the deep parts and skin of the first and second parts of the person to be measured is calculated. Then, the core body temperature of the person to be measured is estimated from the calculated amount of heat, so the core body temperature can be estimated with higher accuracy than in the past.

図1は、本発明の実施例に係る深部体温推定装置の構成を示すブロック図である。FIG. 1 is a block diagram showing the configuration of a core body temperature estimation device according to an embodiment of the present invention. 図2は、測定対象者の身体の2部位2層モデルと、測定対象者の各部位・各層に流入出する熱量を示す図である。FIG. 2 is a diagram showing a two-part, two-layer model of the subject's body and the amount of heat flowing into and out of each part and each layer of the subject. 図3は、本発明の実施例に係る深部体温推定装置の動作を説明するフローチャートである。FIG. 3 is a flow chart for explaining the operation of the core body temperature estimating device according to the embodiment of the present invention. 図4は、本発明の実施例に係る深部体温推定装置を実現するコンピュータの構成例を示すブロック図である。FIG. 4 is a block diagram showing a configuration example of a computer that realizes a core body temperature estimation device according to an embodiment of the present invention.

以下、本発明の実施例について図面を参照して説明する。図1は本発明の実施例に係る深部体温推定装置の構成を示すブロック図である。深部体温推定装置は、第1測定部1と、第2測定部2と、発熱量算出部3と、代謝量算出部4と、呼気蒸散量設定部5(第3熱交換量設定部)と、熱交換量算出部6(第1熱交換量算出部)と、熱交換量算出部7(第2熱交換量算出部)と、温度算出部8とを備えている。
発熱量算出部3と代謝量算出部4と呼気蒸散量設定部5と熱交換量算出部6と熱交換量算出部7とは、熱量算出部9を構成している。
Hereinafter, embodiments of the present invention will be described with reference to the drawings. FIG. 1 is a block diagram showing the configuration of a core body temperature estimation device according to an embodiment of the present invention. The core body temperature estimating device includes a first measuring unit 1, a second measuring unit 2, a calorific value calculating unit 3, a metabolic rate calculating unit 4, and an exhaled transpiration amount setting unit 5 (third heat exchange amount setting unit). , a heat exchange amount calculator 6 (first heat exchange amount calculator), a heat exchange amount calculator 7 (second heat exchange amount calculator), and a temperature calculator 8 .
The calorific value calculator 3 , the metabolic rate calculator 4 , the respiratory transpiration rate setter 5 , the heat exchange rate calculator 6 , and the heat exchange rate calculator 7 constitute a heat rate calculator 9 .

第1測定部1は、測定対象者の心拍数または脈拍数を測定する。第1測定部1は、例えば、測定対象者の心電位を計測するウエア型あるいはベルト型の心電計と、心電計が計測した心電位から心拍数を算出する算出部とから構成される。あるいは、第1測定部1は、測定対象者の脈波を計測する腕時計型あるいはイヤホン型の脈波計と、脈波計が計測した脈波から心拍数(脈拍数)を算出する算出部とから構成される。 The first measurement unit 1 measures the heart rate or pulse rate of the person to be measured. The first measurement unit 1 is composed of, for example, a wear-type or belt-type electrocardiograph that measures the electrocardiogram of the person to be measured, and a calculation unit that calculates the heart rate from the electrocardiogram measured by the electrocardiograph. . Alternatively, the first measurement unit 1 includes a wristwatch-type or earphone-type pulse wave meter that measures the pulse wave of the person to be measured, and a calculation unit that calculates the heart rate (pulse rate) from the pulse wave measured by the pulse wave meter. consists of

第2測定部2は、測定対象者の皮膚温度を測定する。第2測定部2は、例えば温度計から構成される。第2測定部2は、第1測定部と一体化した温度計であってもよい。また、第2測定部2を複数の温度計から構成し、測定対象者の複数個所の皮膚温度を測定する構成とすることもできる。 The second measuring unit 2 measures the skin temperature of the person to be measured. The second measuring unit 2 is composed of, for example, a thermometer. The second measuring section 2 may be a thermometer integrated with the first measuring section. Further, the second measuring unit 2 may be composed of a plurality of thermometers and configured to measure the skin temperature at a plurality of locations of the person to be measured.

また、第2測定部2を赤外線サーモグラフィから構成してもよい。赤外線サーモグラフィは、測定対象者から放射された赤外線放射エネルギーを検出して、見かけの温度に変換し、温度分布を示す画像データとして計測する。第2測定部を赤外線サーモグラフィから構成することで、測定対象者の所定の領域の皮膚温度を測定することができる。
第1測定部1と第2測定部2とは、一体型の測定器でもよいし、個別の測定器でもよい。
In addition, the second measurement section 2 may be configured by an infrared thermography. Infrared thermography detects infrared radiant energy emitted from a person to be measured, converts it into apparent temperature, and measures it as image data showing temperature distribution. By configuring the second measurement unit with an infrared thermography, it is possible to measure the skin temperature of a predetermined region of the measurement subject.
The first measuring unit 1 and the second measuring unit 2 may be integrated measuring instruments or separate measuring instruments.

熱量算出部9は、測定対象者の心拍数と皮膚温度とに基づいて、測定対象者の体幹部(第1部位)と四肢部(第2部位)の深部と皮膚とに流入出する熱量を算出する。
温度算出部8は、熱量算出部9によって算出された熱量に基づいて、測定対象者の深部温度を算出する。
Based on the heart rate and skin temperature of the measurement subject, the calorie calculation unit 9 calculates the amount of heat flowing into and out of the deep parts of the trunk (first part) and extremities (second part) of the measurement subject and the skin. calculate.
The temperature calculator 8 calculates the deep temperature of the person to be measured based on the heat quantity calculated by the heat quantity calculator 9 .

[深部体温推定方法:時間ステップごとの深部体温変化の計算方法]
本発明は、図2に示すように、測定対象者の身体が体幹部Uと四肢部Lの2部位で構成され、体幹部Uと四肢部Lの2部位がそれぞれ深部層Cと皮膚層Sの2層を有するものとみなす。すなわち、測定対象者の身体は、体幹部深部層UCと、体幹部皮膚層USと、四肢部深部層LCと、四肢部皮膚層LSとからなる。本発明は、取得する皮膚温度と心拍数情報とを基に、測定対象者の各部位・各層に流入出する、時々刻々と変化する熱量を算出し、熱量の変化から、測定対象者の各部位・各層の温度変化を推定する。
[Method for Estimating Core Body Temperature: Method for Calculating Changes in Core Body Temperature at Each Time Step]
In the present invention, as shown in FIG. 2, the body of the person to be measured is composed of two parts, the trunk U and the extremities L, and the two parts of the trunk U and the extremities L are the deep layer C and the skin layer S, respectively. is considered to have two layers of That is, the body of the measurement subject is composed of a trunk deep layer UC, a trunk skin layer US, a limb deep layer LC, and a limb skin layer LS. The present invention, based on the acquired skin temperature and heart rate information, calculates the amount of heat that flows into and out of each part and layer of the measurement subject, which changes from time to time, and from the change in the heat amount, each of the measurement subjects Estimate the temperature change of each part and each layer.

本実施例では、センサデータを取得する時間ステップごとの深部体温変化の計算方法を説明する。本実施例では、上述のとおり、測定対象者の身体を体幹部と四肢部の2部位で構成されるものとみなした場合の計算例を記載しているが、体幹部を上半身に置き替え、四肢部を下半身に置き替えてもよい。つまり、測定対象者の身体を第1部位と第2部位の任意の2部位で構成されるものとしてよい。 In this embodiment, a method for calculating changes in core body temperature for each time step of acquiring sensor data will be described. In this embodiment, as described above, a calculation example is described in which the body of the person to be measured is considered to be composed of two parts, the trunk and the limbs, but the trunk is replaced with the upper body, The limbs may be replaced with the lower body. That is, the body of the person to be measured may be composed of arbitrary two parts, the first part and the second part.

式(1)、式(2)に、それぞれ測定対象者の体幹部深部層UCの温度変化、四肢部深部層LCの温度変化を推定する式の例を示す。TUC[t]は時刻tにおける体幹部深部層UCの温度[℃]、TLC[t]は時刻tにおける四肢部深部層LCの温度[℃]である。 Equations (1) and (2) show examples of equations for estimating the temperature change in the core deep layer UC and the temperature change in the extremities deep layer LC, respectively. T UC [t] is the temperature [° C.] of the core deep layer UC at time t, and T LC [t] is the temperature [° C.] of the extremity deep layer LC at time t.

Figure 2022111548000002
Figure 2022111548000002

Figure 2022111548000003
Figure 2022111548000003

1,Uは測定対象者の運動による体幹部Uの発熱量[W]、Q1,Lは運動による四肢部Lの発熱量[W]である。Q2,Uは測定対象者の体幹部深部層UCの代謝量[W]、Q2,Lは四肢部深部層LCの代謝量[W]である。Q3(第3熱交換量)は測定対象者の呼気蒸散量[W]である。Q4,U(第1熱交換量)は測定対象者の体幹部Uにおける深部と皮膚間の熱交換量[W]、Q4,L(第1熱交換量)は四肢部Lにおける深部と皮膚間の熱交換量[W]である。Q5(第2熱交換量)は測定対象者の深部における体幹部Uと四肢部L間の熱交換量[W]である。 Q 1,U is the calorific value [W] of the trunk U due to the exercise of the measurement subject, and Q 1,L is the calorific value [W] of the extremities L due to the exercise. Q 2,U is the metabolic rate [W] of the body trunk deep layer UC of the measurement subject, and Q 2,L is the metabolic rate [W] of the extremity deep layer LC. Q 3 (third heat exchange amount) is the expiratory transpiration amount [W] of the subject. Q 4 ,U (first heat exchange amount) is the amount of heat exchange [W] between the deep part of the trunk U of the subject and the skin; It is the amount of heat exchange [W] between skins. Q 5 (second heat exchange amount) is the heat exchange amount [W] between the trunk U and the extremities L in the deep part of the person to be measured.

WCUCは測定対象者の体幹部深部層UCの熱容量[J/℃]、WCLCは四肢部深部層LCの熱容量[J/℃]である。Δtは計算ステップ時間であり、例えば1[s]以下とする。また、時刻tにおける平均皮膚温Tsk[t]は式(3)のようになる。時刻t+Δtにおける深部体温T[t+Δt]は式(4)のようになる。 WC UC is the heat capacity [J/°C] of the core deep layer UC of the subject, and WC LC is the heat capacity of the extremity deep layer LC [J/°C]. Δt is a calculation step time, and is set to 1 [s] or less, for example. Also, the average skin temperature T sk [t] at time t is given by equation (3). The core body temperature T[t+Δt] at time t+Δt is given by equation (4).

Figure 2022111548000004
Figure 2022111548000004

Figure 2022111548000005
Figure 2022111548000005

sf_conf_USは測定対象者の体表面全体に対する体幹部Uの表面積の割合[%]、sf_conf_LSは体表面全体に対する四肢部Lの表面積の割合[%]である。熱容量WCUC,WCLC、割合sf_conf_US,sf_conf_LSは、それぞれ既知の値であり、実際に即した値を用いればよく、計算の実施において適宜設定すればよい。TUS[t]は時刻tにおける体幹部皮膚層USの温度[℃]、TLS[t]は時刻tにおける四肢部皮膚層LSの温度[℃]である。TUS[t],TLS[t]は第2測定部2によって測定される。 sf_conf_US is the ratio [%] of the surface area of the trunk U to the entire body surface of the subject, and sf_conf_LS is the ratio [%] of the surface area of the extremities L to the entire body surface. The heat capacities WC UC , WC LC and the ratios sf_conf_US, sf_conf_LS are known values, respectively, and may be set as appropriate in the execution of calculations. T US [t] is the temperature [° C.] of the trunk skin layer US at time t, and T LS [t] is the temperature [° C.] of the limb skin layer LS at time t. T US [t] and T LS [t] are measured by the second measuring section 2 .

図3は本実施例の深部体温推定装置の動作を説明するフローチャートである。第1測定部1は、測定対象者の心拍数を測定する(図3ステップS100)。第2測定部2は、測定対象者の皮膚温度を測定する(図3ステップS101)。このとき、第2測定部2は、少なくとも体幹部Uと四肢部Lの2箇所の皮膚温度を測定する。 FIG. 3 is a flow chart for explaining the operation of the core body temperature estimating device of this embodiment. The first measurement unit 1 measures the heart rate of the person to be measured (step S100 in FIG. 3). The second measuring unit 2 measures the skin temperature of the person to be measured (step S101 in FIG. 3). At this time, the second measuring unit 2 measures the skin temperature at least at two points, the trunk U and the extremities L. As shown in FIG.

[平均皮膚温Tskの算出]
温度算出部8は、第2測定部2によって測定された体幹部皮膚層USの温度TUS[t][℃]と四肢部皮膚層LSの温度TLS[t][℃]とに基づいて、時刻tにおける平均皮膚温Tsk[t]を式(3)のように算出する(図3ステップS102)。
[Calculation of average skin temperature Tsk ]
The temperature calculator 8 calculates the temperature T US [t] [°C] of the trunk skin layer US and the temperature T LS [t] [°C] of the limb skin layer LS measured by the second measurement unit 2. , the average skin temperature T sk [t] at time t is calculated as shown in equation (3) (step S102 in FIG. 3).

[運動による発熱量Q1の算出]
次に、発熱量算出部3は、第1測定部1によって測定された心拍数に基づいて、測定対象者の運動による体幹部Uの深部層における発熱量Q1,U[W]、運動による四肢部Lの深部層における発熱量Q1,L[W]をそれぞれ式(5)、式(6)のように算出する(図3ステップS103)。
[Calculation of calorific value Q 1 due to exercise]
Next, based on the heart rate measured by the first measurement unit 1, the calorific value calculation unit 3 calculates the calorific value Q 1,U [W] in the deep layer of the trunk U due to the exercise of the measurement subject, The calorific value Q 1,L [W] in the deep layer of the extremities L is calculated by the equations (5) and (6) respectively (step S103 in FIG. 3).

Figure 2022111548000006
Figure 2022111548000006

Figure 2022111548000007
Figure 2022111548000007

式(5)、式(6)は、文献「“健康づくりのための運動指針2006 ~生活習慣病予防のために~”,厚生労働省,2006年,<https://www.mhlw.go.jp/shingi/2006/07/dl/s0719-3c.pdf>」に開示されている。weightは測定対象者の重量[kg]、ex_conf_Uは測定対象者の全身に対する体幹部Uの筋肉量の割合[%]、ex_conf_Lは全身に対する四肢部Lの筋肉量の割合[%]である。重量weight、割合ex_conf_U,ex_conf_Lは、それぞれ既知の値であり、実際に即した値を用いればよく、計算の実施において適宜設定すればよい。 Equations (5) and (6) are obtained from the document ““Exercise guidelines for health promotion 2006-For prevention of lifestyle-related diseases-”, Ministry of Health, Labor and Welfare, 2006, <https://www.mhlw.go. jp/shingi/2006/07/dl/s0719-3c.pdf>”. weight is the weight [kg] of the person to be measured, ex_conf_U is the ratio [%] of the muscle mass of the trunk U to the whole body of the person to be measured, and ex_conf_L is the ratio [%] of the muscle mass of the extremities L to the whole body. The weight weight and the ratios ex_conf_U and ex_conf_L are known values, and may be set according to the actual situation, and may be appropriately set in the execution of the calculation.

また、METs[t]はメッツ数を表している。発熱量算出部3は、第1測定部1が取得した時刻tにおける心拍数HR[t][bpm]、安静時心拍数HRrest、最大心拍数HRmaxを用いて、式(7)または式(8)のいずれか一方の式によりMETs[t]を算出すればよい。 METs[t] represents the number of METs. The calorific value calculation unit 3 uses the heart rate HR [t] [bpm], the resting heart rate HR rest , and the maximum heart rate HR max at the time t acquired by the first measurement unit 1 to calculate Equation (7) or Equation METs[t] may be calculated by any one of the equations (8).

Figure 2022111548000008
Figure 2022111548000008

Figure 2022111548000009
Figure 2022111548000009

式(8)は、文献「J.R.Wicks,et al.,“HR Index-A Simple Method for the Prediction of Oxygen Uptake”,Medicine and Science in Sports and Exercise,2011」に開示されている。
安静時心拍数HRrest、最大心拍数HRmaxは、それぞれ過去の測定で得られた既知の値として、実際に即した値を用いればよく、計算の実施において適宜設定すればよい。また、式(5)、式(6)で用いるMETs[t]に関しては、瞬時値だけでなく、時間平均値(例えば、6分間の平均値)を用いてもよい。
Equation (8) is disclosed in the document “JR Wicks, et al., “HR Index-A Simple Method for the Prediction of Oxygen Uptake”, Medicine and Science in Sports and Exercise, 2011”.
The resting heart rate HR rest and the maximum heart rate HR max may each be known values obtained from past measurements, and may be set as appropriate for actual calculations. Also, for METs[t] used in equations (5) and (6), not only the instantaneous value but also the time average value (for example, the average value for 6 minutes) may be used.

[代謝量Q2の算出]
次に、代謝量算出部4は、温度算出部8がΔt前に算出した、時刻tにおける深部体温の推定値T[t][℃]に基づいて、測定対象者の体幹部深部層UCの代謝量Q2,U[W]、四肢部深部層LCの代謝量Q2,L[W]をそれぞれ式(9)、式(10)のように算出する(図3ステップS104)。
[Calculation of metabolic rate Q2 ]
Next, the metabolic rate calculation unit 4 calculates the core body temperature UC of the measurement subject based on the estimated value T [t] [° C.] of the core body temperature at time t, which was calculated by the temperature calculation unit 8 before Δt. The metabolic rate Q 2,U [W] and the metabolic rate Q 2,L [W] in the deep layers of the extremities LC are calculated by the equations (9) and (10), respectively (step S104 in FIG. 3).

Figure 2022111548000010
Figure 2022111548000010

Figure 2022111548000011
Figure 2022111548000011

式(9)、式(10)は、文献「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」に開示されている。Askinは代謝に関する定数、volume_USは測定対象者の体幹部皮膚層USの体積[m3]、volume_LSは四肢部皮膚層LSの体積[m3]、weight_Uは測定対象者の体幹部Uの重量[kg]、weight_Lは四肢部Lの重量[kg]である。定数Askin、体積volume_US,volume_LS、重量weight_U,weight_Lは、それぞれ既知の値であり、実際に即した値を用いればよく、計算の実施において適宜設定すればよい。 Formulas (9) and (10) are obtained from the document "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, Volume 32, Issue 8, 1984. A skin is a constant related to metabolism, volume_US is the volume of the trunk skin layer US of the measurement subject [m 3 ], volume_LS is the volume of the extremity skin layer LS [m 3 ], and weight_U is the weight of the trunk U of the measurement subject. [kg], weight_L is the weight [kg] of the extremities L. The constant A skin , the volumes volume_US and volume_LS, and the weights weight_U and weight_L are known values, and may be set as appropriate in the execution of the calculation, using values that are practical.

T[0]は深部体温T[t]の演算時の初期値である。本実施例では、深部体温T[t]の初期値T[0]として、別の温度測定装置によって測定した実測値または一般的な安静時の既知の値などの実際に即した値を用いる。時刻t=0の初回の代謝量Q2の算出では、T[t]=T[0]である。
また、代謝量算出部4は、式(9)、式(10)のMを式(11)のように算出する。
T[0] is the initial value when the core body temperature T[t] is calculated. In this embodiment, as the initial value T[0] of the core body temperature T[t], a realistic value such as an actual value measured by another temperature measuring device or a known value at general rest is used. In the initial calculation of the metabolic rate Q 2 at time t=0, T[t]=T[0].
Also, the metabolic rate calculation unit 4 calculates M in the formulas (9) and (10) as in the formula (11).

Figure 2022111548000012
Figure 2022111548000012

式(11)は、文献「AA Ganpule,et al.,“Interindividual variability in sleeping metabolic rate in Japanese subjects”,European Journal of Clinical Nutrition,volume 61,2007」に開示されている。weight_Cは測定対象者の深部の重量[kg]、heightは測定対象者の身長[cm]、ageは測定対象者の年齢である。sexcoefは測定対象者が男性の場合に0.5473、測定対象者が女性の場合に0.5473×2となる定数である。activity_levelは身体活動レベル、Acoefは代謝調整用のパラメータである。重量weight_C、身長height、年齢age、定数sexcoef、activity_level、パラメータAcoefは、それぞれ既知の値であり、実際に即した値を用いればよく、計算の実施において適宜設定すればよい。 Formula (11) is disclosed in the document “AA Ganpule, et al., “Interindividual variability in sleeping metabolic rate in Japanese subjects”, European Journal of Clinical Nutrition, volume 61, 2007”. weight_C is the deep weight [kg] of the measurement subject, height is the height [cm] of the measurement subject, and age is the age of the measurement subject. sexcoef is a constant that is 0.5473 when the measurement subject is male and 0.5473×2 when the measurement subject is female. activity_level is a physical activity level, and A coef is a parameter for metabolic adjustment. The weight weight_C, height height, age age, constant sexcoef, activity_level, and parameter A coef are known values, and may be set appropriately in the execution of the calculation, using values that match the actual situation.

activity_levelの設定例としては、測定対象者の生活の大部分が座位で、静的な活動が中心の場合は1.5程度とすればよい。測定対象者が座位中心の仕事をしているが、職場内での移動や立位での作業、接客等を含む場合、あるいは通勤、買物、家事、軽いスポーツ等のいずれかを含む場合は、activity_levelを1.75程度とすればよい。測定対象者が移動や立位の多い仕事への従事者である場合、あるいはスポーツなど余暇における活発な運動習慣を持っている場合は、activity_levelを2.0程度とすればよい。このように、測定対象者の運動状況に応じてactivity_levelを適宜設定すればよい。 As an example of setting activity_level, if the person to be measured spends most of his/her life in a sitting position and mainly engages in static activities, the activity_level may be set to about 1.5. If the person to be measured works mainly in a sitting position, but includes movement within the workplace, work in a standing position, customer service, etc., or includes any of commuting, shopping, housework, light sports, etc., Activity_level should be about 1.75. If the person to be measured is engaged in a job that requires a lot of movement or standing, or if he/she has a habit of active exercise in leisure such as sports, activity_level should be set to about 2.0. Thus, the activity_level may be appropriately set according to the exercise status of the person to be measured.

[呼気蒸散量Q3の設定]
次に、呼気蒸散量設定部5は、測定対象者の呼気による蒸散量Q3[W]を式(12)のように設定する(図3ステップS105)。呼気蒸散量Q3[W]は、外気と測定対象者の深部間の熱交換量に相当する。
[Setting of expiratory transpiration Q3 ]
Next, the exhaled transpiration setting unit 5 sets the exhaled transpiration Q 3 [W] of the person to be measured according to the equation (12) (step S105 in FIG. 3). The expiratory transpiration amount Q 3 [W] corresponds to the amount of heat exchange between the outside air and the deep part of the subject.

Figure 2022111548000013
Figure 2022111548000013

[深部と皮膚間の熱交換量Q4の算出]
次に、熱交換量算出部6は、第1測定部1によって測定された心拍数HR[t][bpm]と、第2測定部2によって測定された体幹部皮膚層USの温度TUS[t][℃]と四肢部皮膚層LSの温度TLS[t][℃]と、温度算出部8が算出した平均皮膚温Tsk[t][℃]と、温度算出部8がΔt前に算出した、時刻tにおける体幹部深部層UCの温度の推定値TUC[t][℃]と四肢部深部層LCの温度の推定値TLC[t][℃]と深部体温の推定値T[t][℃]とに基づいて、測定対象者の体幹部Uにおける深部と皮膚間の熱交換量Q4,U[W]、四肢部Lにおける深部と皮膚間の熱交換量Q4,L[W]をそれぞれ式(13)、式(14)により算出する(図3ステップS106)。
[Calculation of heat exchange amount Q 4 between deep part and skin]
Next, the heat exchange amount calculation unit 6 calculates the heart rate HR [t] [bpm] measured by the first measurement unit 1 and the temperature T US [ t] [°C], the temperature T LS [t] [°C] of the limb skin layer LS, the average skin temperature T sk [t] [°C] calculated by the temperature calculator 8, and the temperature calculated by the temperature calculator 8 before Δt Estimated temperature T UC [t] [°C] of core deep layer UC, estimated temperature T LC [t] [°C] of limb deep layer LC, and estimated core body temperature at time t Based on T [t] [°C], the amount of heat exchange Q 4,U [W] between the deep part of the trunk U of the measurement subject and the skin, and the amount of heat exchange Q 4 between the deep part and the skin of the extremities L , L [W] are calculated by the equations (13) and (14), respectively (step S106 in FIG. 3).

Figure 2022111548000014
Figure 2022111548000014

Figure 2022111548000015
Figure 2022111548000015

割合sf_conf_US,sf_conf_LSについては上記のとおりである。hxは測定対象者の皮膚と深部間の熱交換係数である。熱交換量算出部6は、熱交換係数hxを式(15)により算出することができる。 The ratios sf_conf_US and sf_conf_LS are as described above. hx is the heat exchange coefficient between the skin of the person to be measured and the deep part. The heat exchange amount calculation unit 6 can calculate the heat exchange coefficient hx by Equation (15).

Figure 2022111548000016
Figure 2022111548000016

min(X,Y)はXとYのうちいずれか小さい方を採用することを意味する。METs[t]については上記のとおりである。a,b,e,hx0,hx1,hx_maxは熱交換係数に関わるパラメータである。パラメータa,b,e,hx0,hx1,hx_maxは、それぞれ既知の値であり、実際に即した値を用いればよく、実験を通じて適宜設定すればよい。 min(X, Y) means to adopt the smaller one of X and Y. METs[t] is as described above. a, b, e, hx0, hx1, and hx_max are parameters related to heat exchange coefficients. The parameters a, b, e, hx0, hx1, and hx_max are known values, and may be set appropriately through experiments.

sk[0]は平均皮膚温Tsk[t]の演算時の初期値である。本実施例では、体幹部深部層UCの温度の初期値TUC[0]、四肢部深部層LCの温度の初期値TLC[0]として、別の温度測定装置によって測定した実測値または一般的な安静時の既知の値などの実際に即した値を用いる。時刻t=0の初回の熱交換量Q4の算出では、Tsk[t]=Tsk[0]である。 T sk [0] is the initial value when calculating the average skin temperature T sk [t]. In this embodiment, as the initial value T UC [0] of the temperature of the core deep layer UC and the initial value T LC [0] of the temperature of the extremity deep layer LC, measured values measured by another temperature measuring device or general Use realistic values, such as known resting values. In the first calculation of the heat exchange amount Q 4 at time t=0, T sk [t]=T sk [0].

[体幹部と四肢部間の熱交換量Q5の算出]
次に、熱交換量算出部7は、第1測定部1によって測定された心拍数HR[t][bpm]と、温度算出部8が算出した平均皮膚温Tsk[t][℃]と、温度算出部8がΔt前に算出した、時刻tにおける体幹部深部層UCの温度の推定値TUC[t][℃]と四肢部深部層LCの温度TLC[t][℃]と深部体温の推定値T[t][℃]とに基づいて、測定対象者の深部における体幹部Uと四肢部L間の熱交換量Q5[W]を式(16)により算出する(図3ステップS107)。
[Calculation of heat exchange amount Q 5 between trunk and extremities]
Next, the heat exchange amount calculation unit 7 calculates the heart rate HR [t] [bpm] measured by the first measurement unit 1, the average skin temperature T sk [t] [°C] calculated by the temperature calculation unit 8, and , the estimated temperature value T UC [t] [°C] of the core deep layer UC at time t and the temperature T LC [t] [°C] of the extremity deep layer LC at time t calculated by the temperature calculator 8 before Δt, and Based on the estimated core body temperature T [t] [° C.], the amount of heat exchange Q 5 [W] between the trunk U and the extremities L in the deep part of the subject is calculated by equation (16) (Fig. 3 step S107).

Figure 2022111548000017
Figure 2022111548000017

hccは測定対象者の深部における体幹部Uと四肢部Lの熱交換係数である。熱交換量算出部7は、熱交換係数hccを式(17)により算出することができる。 hcc is the heat exchange coefficient between the trunk U and the extremities L in the deep part of the person to be measured. The heat exchange amount calculation unit 7 can calculate the heat exchange coefficient hcc by Equation (17).

Figure 2022111548000018
Figure 2022111548000018

f,hcc0は熱交換係数に関わるパラメータである。パラメータf,hcc0は、それぞれ既知の値であり、実際に即した値を用いればよく、実験を通じて適宜設定すればよい。hcc_Tは熱交換係数hccの温度寄与分、hcc_Mは熱交換係数hccのMETs寄与分である。熱交換量算出部7は、熱交換係数hccの温度寄与分hcc_Tを式(18)により算出することができる。 f and hcc0 are parameters related to the heat exchange coefficient. The parameters f and hcc0 are known values, and may be set appropriately through experiments. hcc_T is the temperature contribution of the heat exchange coefficient hcc, and hcc_M is the METs contribution of the heat exchange coefficient hcc. The heat exchange amount calculation unit 7 can calculate the temperature contribution hcc_T of the heat exchange coefficient hcc by Equation (18).

Figure 2022111548000019
Figure 2022111548000019

hcc_Tmaxはhcc_Tの規定の上限値である。また、熱交換量算出部7は、熱交換係数hccのMETs寄与分hcc_Mを式(19)により算出することができる。 hcc_Tmax is the specified upper limit of hcc_T. Further, the heat exchange amount calculation unit 7 can calculate the METs contribution hcc_M of the heat exchange coefficient hcc by using Equation (19).

Figure 2022111548000020
Figure 2022111548000020

hcc_Mmaxはhcc_Mの規定の上限値である。hcc1は熱交換係数に関わるパラメータである。上限値hcc_Tmax,hcc_Mmax、パラメータhcc1は、それぞれ既知の値であり、実際に即した値を用いればよく、実験を通じて適宜設定すればよい。aveMETs[t]は、メッツ数の時間平均値(例えば、6分間の平均値)である。パラメータa,bについては上記のとおりである。 hcc_Mmax is the specified upper limit of hcc_M. hcc1 is a parameter related to the heat exchange coefficient. The upper limit values hcc_Tmax, hcc_Mmax, and the parameter hcc1 are known values, and may be set appropriately through experiments. aveMETs[t] is the time average value of the METs number (for example, the average value for 6 minutes). Parameters a and b are as described above.

こうして、発熱量算出部3と代謝量算出部4と呼気蒸散量設定部5と熱交換量算出部6と熱交換量算出部7とにより、各熱量Q1,U,Q1,L,Q2,U,Q2,L,Q3,Q4,U,Q4,L,Q5を算出することができる。 In this way, the calorific value calculator 3, the metabolic rate calculator 4, the exhaled transpiration rate setting unit 5, the heat exchange rate calculator 6, and the heat exchange rate calculator 7 calculate the respective heat amounts Q 1,U , Q 1,L , Q 2,U , Q2 ,L , Q3 , Q4 ,U , Q4 ,L , and Q5 can be calculated.

温度算出部8は、時刻tにおける測定対象者の体幹部深部層UCの温度の推定値TUC[t][℃]と、体幹部Uの深部層における発熱量Q1,U[W]と、体幹部深部層UCの代謝量Q2,U[W]と、呼気蒸散量Q3[W]と、体幹部Uにおける深部と皮膚間の熱交換量Q4,U[W]と、深部における体幹部Uと四肢部L間の熱交換量Q5[W]とに基づいて、Δt後の体幹部深部層UCの温度の推定値TUC[t+Δt][℃]を式(1)により算出する(図3ステップS108)。 The temperature calculator 8 calculates an estimated value T UC [t] [° C ] of the temperature of the deep layer UC of the torso of the person to be measured at time t, , the metabolic rate Q 2,U [W] in the core deep layer UC, the expiratory transpiration rate Q 3 [W], the heat exchange rate Q 4,U [W] between the deep region and the skin in the core U, and the deep region Based on the heat exchange amount Q 5 [W] between the trunk U and the extremities L at the time, the estimated value T UC [t + Δt] [°C] of the temperature of the trunk deep layer UC after Δt is calculated by the formula (1) Calculate (step S108 in FIG. 3).

温度算出部8は、時刻tにおける測定対象者の四肢部深部層LCの温度の推定値TLC[t][℃]と、四肢部Lの深部層における発熱量Q1,L[W]と、四肢部深部層LCの代謝量Q2,L[W]と、四肢部Lにおける深部と皮膚間の熱交換量Q4,L[W]と、深部における体幹部Uと四肢部L間の熱交換量Q5[W]とに基づいて、Δt後の四肢部深部層LCの温度の推定値TLC[t+Δt][℃]を式(2)により算出する(図3ステップS109)。 The temperature calculator 8 calculates an estimated value T LC [t] [°C] of the temperature of the deep layers LC of the limbs L of the measurement subject at time t, and the amount of heat generated Q 1,L [W] in the deep layers of the limbs L , the metabolic rate Q 2,L [W] of the limb deep layer LC, the heat exchange rate Q 4,L [W] between the deep part and the skin in the limb L, and the deep part between the trunk U and the limb L Based on the amount of heat exchange Q 5 [W], an estimated value T LC [t+Δt] [° C.] of the temperature of the extremity deep layer LC after Δt is calculated by Equation (2) (step S109 in FIG. 3).

こうして、温度の推定値TUC[t+Δt],TLC[t+Δt]をそれぞれ逐次計算することができる。
そして、温度算出部8は、Δt後の深部体温の推定値T[t+Δt][℃]を式(4)により算出する(図3ステップS110)。
Thus, temperature estimates T UC [t+Δt] and T LC [t+Δt] can be calculated sequentially.
Then, the temperature calculator 8 calculates the estimated value T[t+Δt][° C.] of the core body temperature after Δt using the equation (4) (step S110 in FIG. 3).

深部体温推定装置は、以上のステップS100~S110の処理を周期Δt毎に行う。Δt後の次の計算では、直前の回で算出したTUC[t+Δt],TLC[t+Δt],T[t+Δt]を、それぞれTUC[t],TLC[t],T[t]として、ステップS100~S110の処理を行うようにすればよい。 The core body temperature estimating apparatus performs the above steps S100 to S110 every cycle Δt. In the next calculation after Δt, T UC [t + Δt], T LC [t + Δt], and T [t + Δt] calculated in the previous calculation are used as T UC [t], T LC [t], and T [t], respectively. , steps S100 to S110 may be performed.

本実施例では、測定対象者の心拍数だけでなく、測定対象者の皮膚温度を測定して、測定対象者の各部位・各層に流入出する熱量を算出し、算出した熱量から測定対象者の深部体温を推定するので、従来よりも高精度に深部体温を推定することができる。 In this embodiment, not only the heart rate of the person to be measured but also the skin temperature of the person to be measured is measured, the amount of heat flowing into and out of each part and layer of the person to be measured is calculated, and the calculated amount of heat is used by the person to be measured. Since the core body temperature is estimated, it is possible to estimate the core body temperature with higher accuracy than before.

なお、本実施例では、測定対象者の心拍数を用いる場合について説明しているが、心拍数の代わりに脈拍数を用いてもよい。 In addition, although the present embodiment describes the case of using the heart rate of the person to be measured, the pulse rate may be used instead of the heart rate.

本実施例の発熱量算出部3と代謝量算出部4と呼気蒸散量設定部5と熱交換量算出部6と熱交換量算出部7と温度算出部8とは、CPU(Central Processing Unit)、記憶装置及びインタフェースを備えたコンピュータと、これらのハードウェア資源を制御するプログラムによって実現することができる。このコンピュータの構成例を図4に示す。 The calorific value calculator 3, the metabolic rate calculator 4, the exhaled transpiration rate setter 5, the heat exchange rate calculator 6, the heat exchange rate calculator 7, and the temperature calculator 8 of this embodiment are CPU (Central Processing Unit). , a computer having a storage device and an interface, and a program that controls these hardware resources. A configuration example of this computer is shown in FIG.

コンピュータは、CPU200と、記憶装置201と、インタフェース装置(I/F)202とを備えている。I/F202には、第1測定部1と第2測定部2などが接続される。このようなコンピュータにおいて、本発明の深部体温推定方法を実現させるための深部体温推定プログラムは記憶装置201に格納される。CPU200は、記憶装置201に格納されたプログラムに従って本実施例で説明した処理を実行する。プログラムをネットワークを通して提供することも可能である。 The computer comprises a CPU 200 , a storage device 201 and an interface device (I/F) 202 . The I/F 202 is connected with the first measurement unit 1, the second measurement unit 2, and the like. In such a computer, a core body temperature estimation program for implementing the core body temperature estimation method of the present invention is stored in the storage device 201 . The CPU 200 executes the processing described in this embodiment according to the programs stored in the storage device 201 . It is also possible to provide the program through a network.

また、発熱量算出部3と代謝量算出部4と呼気蒸散量設定部5と熱交換量算出部6と熱交換量算出部7と温度算出部8とを、複数のコンピュータ機器に分散させるようにしてもよい。 Further, the calorific value calculation unit 3, the metabolic amount calculation unit 4, the exhalation amount setting unit 5, the heat exchange amount calculation unit 6, the heat exchange amount calculation unit 7, and the temperature calculation unit 8 are distributed among a plurality of computer devices. can be

本発明は、人体の深部体温を推定する技術に適用することができる。 INDUSTRIAL APPLICABILITY The present invention can be applied to techniques for estimating core body temperature of the human body.

1…第1測定部、2…第2測定部、3…発熱量算出部、4…代謝量算出部、5…呼気蒸散量設定部、6,7…熱交換量算出部、8…温度算出部、9…熱量算出部。 DESCRIPTION OF SYMBOLS 1... 1st measurement part 2... 2nd measurement part 3... Calorific value calculation part 4... Metabolic amount calculation part 5... Evaporation amount setting part 6, 7... Heat exchange amount calculation part 8... Temperature calculation Part 9: Calorie calculation part.

Claims (7)

測定対象者の心拍数または脈拍数を測定するように構成された第1測定部と、
前記測定対象者の皮膚温度を測定するように構成された第2測定部と、
前記測定対象者の心拍数または脈拍数と皮膚温度とに基づいて、前記測定対象者の第1部位と第2部位の深部と皮膚とに流入出する熱量を算出するように構成された熱量算出部と、
前記熱量に基づいて、前記測定対象者の深部温度を算出するように構成された温度算出部とを備えることを特徴とする深部体温推定装置。
a first measuring unit configured to measure the heart rate or pulse rate of the person to be measured;
a second measuring unit configured to measure the skin temperature of the measurement subject;
A calorie calculation configured to calculate an amount of heat flowing into and out of deep parts and skin of a first region and a second region of the measurement subject based on the measurement subject's heart rate or pulse rate and skin temperature. Department and
A core body temperature estimating device, comprising: a temperature calculator configured to calculate a core body temperature of the measurement subject based on the amount of heat.
請求項1記載の深部体温推定装置において、
前記第1部位は、前記測定対象者の体幹部であり、
前記第2部位は、前記測定対象者の四肢部であり、
前記熱量算出部は、
前記測定対象者の心拍数または脈拍数に基づいて、前記測定対象者の運動による前記体幹部と前記四肢部のそれぞれの発熱量を算出するように構成された発熱量算出部と、
前記測定対象者の心拍数または脈拍数と皮膚温度と前記温度算出部によって直前に算出された深部温度とに基づいて、前記測定対象者の前記体幹部と前記四肢部のそれぞれにおける深部と皮膚間の第1熱交換量を算出するように構成された第1熱交換量算出部と、
前記測定対象者の心拍数または脈拍数と皮膚温度と前記温度算出部によって直前に算出された深部温度とに基づいて、前記測定対象者の深部における前記体幹部と前記四肢部間の第2熱交換量を算出するように構成された第2熱交換量算出部と、
前記温度算出部によって直前に算出された深部温度に基づいて、前記測定対象者の体幹部深部層の代謝量と四肢部深部層の代謝量を算出するように構成された代謝量算出部とを含むことを特徴とする深部体温推定装置。
The core body temperature estimating device according to claim 1,
The first part is the trunk of the measurement subject,
The second part is the extremities of the measurement subject,
The heat quantity calculation unit
a calorific value calculation unit configured to calculate calorific value of each of the trunk and extremities due to exercise of the measurement subject based on the heart rate or pulse rate of the measurement subject;
Based on the heart rate or pulse rate of the measurement subject, the skin temperature, and the deep temperature calculated immediately before by the temperature calculation unit, the distance between the deep part and the skin in each of the trunk and the extremities of the measurement subject A first heat exchange amount calculation unit configured to calculate the first heat exchange amount of
Based on the heart rate or pulse rate of the measurement subject, the skin temperature, and the deep temperature calculated immediately before by the temperature calculation unit, a second heat between the trunk and the extremities in the deep part of the measurement subject a second heat exchange amount calculator configured to calculate an exchange amount;
a metabolic rate calculator configured to calculate the metabolic rates of the deep layers of the trunk and the deep layers of the extremities of the measurement subject based on the deep temperature calculated immediately before by the temperature calculator; A core body temperature estimation device comprising:
請求項2記載の深部体温推定装置において、
前記第2測定部は、前記体幹部と前記四肢部のそれぞれの皮膚温度を測定し、
前記第1熱交換量算出部は、前記第2測定部の測定結果を用いて前記測定対象者の前記体幹部と前記四肢部のそれぞれにおける深部と皮膚間の第1熱交換量を算出し、
前記第2熱交換量算出部は、前記第2測定部の測定結果を用いて前記測定対象者の深部における前記体幹部と前記四肢部間の第2熱交換量を算出することを特徴とする深部体温推定装置。
The core body temperature estimating device according to claim 2,
The second measuring unit measures the skin temperature of each of the trunk and the extremities,
The first heat exchange amount calculation unit uses the measurement result of the second measurement unit to calculate a first heat exchange amount between the deep part and the skin in each of the trunk and the extremities of the measurement subject,
The second heat exchange amount calculation unit calculates a second heat exchange amount between the trunk and the extremities in the deep part of the measurement subject using the measurement result of the second measurement unit. Core body temperature estimator.
請求項2または3記載の深部体温推定装置において、
前記熱量算出部は、外気と前記測定対象者の深部間の第3熱交換量を設定するように構成された第3熱交換量設定部をさらに含むことを特徴とする深部体温推定装置。
The core body temperature estimating device according to claim 2 or 3,
The deep body temperature estimating apparatus, wherein the heat amount calculating unit further includes a third heat exchange amount setting unit configured to set a third heat exchange amount between the outside air and the deep part of the measurement subject.
請求項4記載の深部体温推定装置において、
前記温度算出部は、
前記測定対象者の前記体幹部と前記四肢部の発熱量Q1,U,Q1,L、前記測定対象者の前記体幹部と前記四肢部のそれぞれにおける深部と皮膚間の第1熱交換量Q4,U,Q4,L、前記測定対象者の深部における前記体幹部と前記四肢部間の第2熱交換量Q5、外気と前記測定対象者の深部間の第3熱交換量Q3、前記測定対象者の体幹部深部層の代謝量Q2,U、四肢部深部層の代謝量Q2,L、前記測定対象者の体幹部深部層の熱容量WCUC、四肢部深部層の熱容量WCLC、直前に算出した体幹部深部層の温度の推定値TUC[t]、四肢部深部層の温度の推定値TLC[t]に基づいて、
Figure 2022111548000021
により、時刻(t+Δt)における体幹部深部層の温度の推定値TUC[t+Δt]、四肢部深部層の温度の推定値TLC[t+Δt]を算出することを特徴とする深部体温推定装置。
The core body temperature estimating device according to claim 4,
The temperature calculation unit
Heat generation amounts Q 1,U , Q 1,L of the trunk and extremities of the measurement subject, first heat exchange amounts between the deep part and the skin of the trunk and extremities of the measurement subject, respectively Q 4,U , Q 4,L , a second amount of heat exchange Q 5 between the trunk and extremities in the deep part of the person to be measured, and a third amount of heat exchange Q between the outside air and the deep part of the person to be measured 3 , the metabolic rate Q 2,U in the deep layers of the trunk of the measurement subject, the metabolic rate Q 2,L in the deep layers of the limbs, the heat capacity WC UC of the deep layers of the limbs, the deep layers of the limbs Based on the heat capacity WC LC , the estimated temperature value T UC [t] of the core deep layer calculated immediately before, and the estimated temperature value T LC [t] of the extremity deep layer,
Figure 2022111548000021
A deep body temperature estimating device characterized by calculating an estimated value T UC [t+Δt] of the temperature of the deep layer of the trunk and an estimated value T LC [t+Δt] of the temperature of the deep layer of the extremities at time (t+Δt).
測定対象者の心拍数または脈拍数を測定する第1ステップと、
前記測定対象者の皮膚温度を測定する第2ステップと、
前記測定対象者の心拍数または脈拍数と皮膚温度とに基づいて、前記測定対象者の第1部位と第2部位の深部と皮膚とに流入出する熱量を算出する第3ステップと、
前記熱量に基づいて、前記測定対象者の深部温度を算出する第4ステップとを含むことを特徴とする深部体温推定方法。
a first step of measuring the heart rate or pulse rate of the person to be measured;
a second step of measuring the skin temperature of the measurement subject;
a third step of calculating the amount of heat flowing into and out of the first and second deep parts and the skin of the measurement subject based on the heart rate or pulse rate of the measurement subject and the skin temperature;
and a fourth step of calculating the core temperature of the person to be measured based on the amount of heat.
請求項6記載の各ステップをコンピュータに実行させることを特徴とする深部体温推定プログラム。 7. A deep body temperature estimating program that causes a computer to execute each step according to claim 6.
JP2021007044A 2021-01-20 2021-01-20 Deep part body temperature estimation device, deep part body temperature estimation method, and deep part body temperature estimation program Pending JP2022111548A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021007044A JP2022111548A (en) 2021-01-20 2021-01-20 Deep part body temperature estimation device, deep part body temperature estimation method, and deep part body temperature estimation program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2021007044A JP2022111548A (en) 2021-01-20 2021-01-20 Deep part body temperature estimation device, deep part body temperature estimation method, and deep part body temperature estimation program

Publications (1)

Publication Number Publication Date
JP2022111548A true JP2022111548A (en) 2022-08-01

Family

ID=82655502

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021007044A Pending JP2022111548A (en) 2021-01-20 2021-01-20 Deep part body temperature estimation device, deep part body temperature estimation method, and deep part body temperature estimation program

Country Status (1)

Country Link
JP (1) JP2022111548A (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016028662A (en) * 2014-07-25 2016-03-03 船井電機株式会社 Sleep evaluation device and sleep evaluation method
JP2017217224A (en) * 2016-06-08 2017-12-14 国立大学法人大阪大学 Deep body temperature estimation device, method thereof, and program
US20190134396A1 (en) * 2016-04-29 2019-05-09 Lifelens Technologies, Llc Monitoring and management of physiologic parameters of a subject
JP2020065823A (en) * 2018-10-26 2020-04-30 日本電信電話株式会社 Deep temperature estimation device

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2016028662A (en) * 2014-07-25 2016-03-03 船井電機株式会社 Sleep evaluation device and sleep evaluation method
US20190134396A1 (en) * 2016-04-29 2019-05-09 Lifelens Technologies, Llc Monitoring and management of physiologic parameters of a subject
JP2017217224A (en) * 2016-06-08 2017-12-14 国立大学法人大阪大学 Deep body temperature estimation device, method thereof, and program
JP2020065823A (en) * 2018-10-26 2020-04-30 日本電信電話株式会社 Deep temperature estimation device

Similar Documents

Publication Publication Date Title
US12076122B2 (en) Pacing templates for performance optimization
Tanda Skin temperature measurements by infrared thermography during running exercise
US10971271B2 (en) Method and system for personalized blood flow modeling based on wearable sensor networks
US9326735B2 (en) Method and apparatus for monitoring cardiac output of a patient in a home environment
WO2017213011A1 (en) Deep body temperature estimation device, and method and program therefor
JP7026905B2 (en) Deep body temperature estimator
US9597038B2 (en) Apparatus for monitoring a subject during exercise
US20220125388A1 (en) Method for estimating core body temperature
EP2818105B1 (en) Apparatus and method for the mobile determination of a physiological stress threshold value
Selvaraj et al. Fully disposable wireless patch sensor for continuous remote patient monitoring
JP7201712B2 (en) Method and apparatus for estimating trends in blood pressure surrogate values
JP2020162952A (en) Deep body temperature estimation device, method, and program
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
Wolff et al. Heat exchanges in wet suits
JP2022111548A (en) Deep part body temperature estimation device, deep part body temperature estimation method, and deep part body temperature estimation program
US20210127989A1 (en) Method of estimation of the hemodynamic status and capacity of effort from patients with cardiovascular and pulmonar diseases
Kamiya et al. Different thermoregulatory responses of people from tropical and temperate zones: A computational study
JP2022111550A (en) Deep part body temperature estimation device, deep part body temperature estimation method, and deep part body temperature estimation program
JP2022111549A (en) Deep part body temperature estimation device, deep part body temperature estimation method, and deep part body temperature estimation program
JP4873792B2 (en) Calorie consumption measuring device
WO2023073750A1 (en) Dehydration symptom estimation device, estimation method, and program
JP7312410B2 (en) Deep body temperature estimator
JP2022045457A (en) Deep body temperature estimation method, device, and program
WO2022180773A1 (en) Notification device, notification method, and program

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

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20231121

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240122

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20240423

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20240624

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20240917