JP2019207124A - 生体内温度測定装置および生体内温度測定方法 - Google Patents

Abstract

【課題】核心部温度の変動をより正確に把握する。【解決手段】二つの温度センサ（２０ｕ、２０ｓ）を含む熱流束センサ（２０）を用いて生体の核心部温度の時系列データを取得する一方、血流センサ（３０）で前記熱流束センサ（２０）近傍の血流量を測定し、演算回路（４０）によって、血流センサ（３０）で測定された血流量と、予め用意した血流量と生体の核心部温度の変動が表皮温度に反映されるまでの遅延時間との関係とに基づいて、生体の核心部温度の変動が表皮温度に反映されるまでの遅延時間を算出し、前記遅延時間に基づいて前記生体の核心部温度に対応づけられた時刻を補正する。【選択図】 図１

Description

本発明は、生体の核心部温度を測定する生体内温度測定装置および生体内温度測定方法に関する。

生体において、表皮から核心部に向かってある一定の深さを超えると、外気温の変化等に左右されない温度領域が存在する（以下、その部分の温度を「核心部温度」または「深部体温」という。）。核心部温度の変動を計測することは、体内リズムの把握に有用であることが知られている。
核心部温度を測定するにあたり、体内留置といった侵襲的な測定ではなく、経皮的な温度測定法は、手軽で、日常的な体温管理に有用である。

中山昭雄、「新生理科体系 第22巻」、医学書院（1987） 中川慎也 他、「MEMS熱流束センサによるウェアラブル深部体温計の提案」、電気学会論文誌Ｅ、135 巻 (2015) 8 号 p. 343-348

しかしながら、従来の経皮的な体温計測装置では、核心部温度の変動を正確に把握することが困難であった。その理由の一つが、核心部温度の変動が表皮で測定される温度に反映されるまでの時間の変動、すなわち遅延時間の変動である。

図１５に示すように、表皮から核心部温度Ｔｃの温度領域までの深さは、血流量に依存することが知られている（非特許文献１、図５９）。体表部の血流は、真皮層に存在する動静脈吻合（ＡＶＡ：Ａｒｔｅｒｉｏｖｅｎｏｕｓ ａｎａｓｔｏｍｏｓｅｓ）と呼ばれる血管が身体の神経活動によって拡張することで増加する。体表部の血流量が増加すると、核心部の熱エネルギーが血流に乗って表層部へ移動してくるので、表皮から核心部温度の温度領域までのみかけの深さが浅くなる。

体表部の血流量が時々刻々と変化すれば、表皮から核心部温度の温度領域までのみかけの深さが変化して、表皮で測定される温度に核心部温度の変動が反映されるまでの遅延時間も変動する。そのため、経皮的に測定された値が、実際にはどの時刻の核心部温度を反映した値なのか分からず、結果的に、核心部温度の変動を正確に捉えることができなかった。

そこで、本発明は、核心部温度の変動をより正確に把握することができる生体内温度測定装置を提供することを目的とする。

上述した目的を達成するために、本発明に係る生体内温度測定装置は、生体の表皮温度を測定する温度センサ（２０ｓ）と、前記生体の体表部から放出される熱流束の大きさを測定する熱流束センサ（２０）と、前記熱流束センサの近傍の血流量を測定する血流センサ（３０）と、前記熱流束センサの近傍の血流量と前記生体の核心部温度の変動が表皮温度に反映されるまでの遅延時間に関するパラメータとの関係を記憶した記憶部（５０）と、前記温度センサで測定された表皮温度と前記熱流束センサで測定された熱流束の大きさから前記生体の核心部温度を算出するとともに、前記熱流束センサの近傍の血流量と前記記憶部に記憶された前記関係とに基づいて前記遅延時間を算出し、前記遅延時間に基づいて前記生体の核心部温度に対応づけられた時刻を補正するように構成された演算回路（４０）とを備える。

本発明に係る生体内温度測定装置において、前記演算回路（４０）は、前記温度センサで測定された表皮温度と前記熱流束センサで測定された熱流束の大きさから前記生体の核心部温度を算出する第１算出部（４１）と、前記熱流束センサの近傍の血流量と前記記憶部に記憶された前記関係とに基づいて前記遅延時間を算出する第２算出部（４３）と、前記遅延時間に基づいて前記生体の核心部温度に対応づけられた時刻を補正する補正部（４４）とを有していてもよい。

本発明に係る生体内温度測定装置において、前記パラメータは、前記生体の核心部温度の変動に対する表皮温度の変動の時定数または前記遅延時間とすることができる。

本発明に係る生体内温度測定装置は、二以上の前記血流センサ（３０−１、３０−２）を備え、前記演算回路（４０）は、二以上の前記血流センサでそれぞれ測定した血流量の代表値を求め、前記関係を用いて前記血流量の代表値に対応する前記パラメータの値を求めるようにしてもよい。

また、本発明に係る生体内温度測定方法は、生体の表皮温度と前記生体の体表部から放出される熱流束の大きさとを測定するステップ（Ｓ１０）と、前記体表部の血流量を測定するステップと（Ｓ２０）と、測定された前記表皮温度と前記熱流束の大きさとから前記生体の核心部温度を算出するとともに、予め用意した、前記血流量と前記生体の核心部温度の変動が表皮温度に反映されるまでの遅延時間に関するパラメータとの関係と、前記体表部の血流量と、に基づいて前記遅延時間を算出し、前記遅延時間に基づいて前記生体の核心部温度に対応づけられた時刻を補正するステップ（Ｓ４０、Ｓ５０、Ｓ７０）と、を有する。

本発明によれば、生体の体表部の血流量に基づいて算出した遅延時間に基づいて核心部温度の時系列データの時刻を補正することによって、核心部温度の変動をより正確に把握することができる。

図１は、本発明の第１の実施の形態に係る生体内温度測定装置の構成を示す図である。 図２は、温度センサを含む熱流束センサの構成を示す図である。 図３は、第１の実施の形態に係る生体内温度測定装置の構成を示す図である。 図４は、第１の実施の形態に係る生体内温度測定装置における熱流束センサと血流センサとの位置関係を説明する図である。 図５は、生体の皮下組織の温度分布を模式的に説明する図である。 図６は、生体の皮下組織の熱流束を模式的に説明する図である。 図７は、生体の皮下組織の熱流束と血流との関係を模式的に説明する図である。 図８は、血流量と、生体の核心部温度の変動に対する表皮温度の変動の時定数との関係を説明する図である。 図９は、第１の実施の形態に係る生体内温度測定装置の動作を説明するフローチャートである。 図１０は、遅延時間に基づく核心部温度の時系列データの補正について説明する図である。 図１１は、遅延時間に基づくデータの補正によって復元された核心部温度変動の周波数情報について説明する図である。 図１２は、本発明の第２の実施の形態に係る生体内温度測定装置の構成を示す図である。 図１３は、本発明の第２の実施の形態に係る生体内温度測定装置における熱流束センサと血流センサとの位置関係を説明する図である。 図１４は、本発明の第２の実施の形態に係る生体内温度測定装置における熱流束センサと血流センサとの位置関係の変形例を説明する図である。 図１５は、表皮から核心部温度の温度領域までの深さと血流量との関係を示す図である。

以下に本発明の実施の形態について図面を参照して説明する。

［第１の実施の形態］
［生体内温度測定装置の測定原理］
本発明の第１の実施の形態に係る生体内温度測定装置１は、図１に示すように、シート状の基材８０の上に、熱流束センサ２０と、血流センサ３０と、演算回路４０と、メモリ５０と、外部とのＩ／Ｆ回路として機能する通信回路６０と、演算回路４０や通信回路６０等に電力を供給する電池７０とを備えている。

ここで熱流束センサ２０は、単位時間・単位面積当たりの熱の移動を計測するデバイスである。本実施の形態においては、図２に示すように、熱抵抗体２０ｒの上面と下面とにそれぞれ温度センサ２０ｕと温度センサ２０ｓと有する熱流束センサ２０を用い、生体９０の体表部から放出される熱流束の大きさを測定すると同時に、温度センサ２０ｓで表皮温度Ｔｓを測定することとする。温度センサ２０ｕ、２０ｓとしては、例えば、公知のサーミスタや、熱電対を用いたサーモパイル、温度により音速が変わることを利用した超音波温度計、温度により光吸収率が変わることを利用した赤外線温度センサその他の光温度計などを用いることができる。

血流センサ３０は、熱流束センサ２０の近傍に配置され、生体９０の体表部の血流量を測定するデバイスである。このような血流センサ３０としては、例えば、レーザを皮膚に照射することで皮下組織内の血流量を測定するレーザドップラー血流計その他の光学式血流量センサや、超音波血流計を用いることができる。

メモリ５０は、生体の核心部温度の変動が表皮温度に反映されるまでの遅延時間に関するパラメータと体表部の血流量との関係を記憶している。
生体の核心部温度の変動が表皮温度に反映されるまでの遅延時間に関するパラメータは、具体的には、生体の核心部温度の変動に対する表皮温度の変動の時定数や、生体の核心部温度の変動が表皮温度に反映されるまでの遅延時間である。血流量と時定数または遅延時間との関係は、テーブルの形式でメモリ５０に記憶しておけばよいが、関数として記憶されていてもよい。なお、血流量と時定数または遅延時間との関係については後述する。
また、メモリ５０は、熱流束センサ２０で測定した結果から算出した核心部温度の時系列データ、すなわち、核心部温度とその核心部温度を測定した時刻とが互いに関連づけられたデータや、後述するように、その時刻を補正した、補正後の時系列データを記憶する。

演算回路４０は、計時部４２を備え、温度センサ２０ｓで測定された表皮温度と熱流束センサ２０で測定された熱流束の大きさから生体の核心部温度Ｔｃを算出して、核心部温度Ｔｃの時系列データを生成するとともに、熱流束センサ２０近傍の体表部の血流量と、メモリ５０に記憶された遅延時間と血流量との関係とに基づいて、生体の核心部温度Ｔｃの変動が表皮温度に反映されるまでの遅延時間Δｔを算出し、この遅延時間Δｔに基づいて、時系列データを構成する各核心部温度Ｔｃに対応づけられた時刻を補正するように構成されている。

このような演算回路４０は、演算装置とコンピュータ・プログラムによって構成することができる。例えば、演算回路４０は、温度センサ２０ｓで測定された表皮温度Ｔsと熱流束センサ２０で測定された熱流束の大きさから、生体の核心部温度Ｔｃを算出し、計時部４２によって計時された核心部温度Ｔｃの計測時刻とともに、核心部温度Ｔｃの時系列データを生成する第１算出部４１と、熱流束センサ２０近傍の体表部の血流量とメモリ５０に記憶された、血流量と遅延時間との関係とに基づいて、生体の核心部温度Ｔｃの変動が表皮温度Ｔｓに反映されるまでの遅延時間Δｔを算出する第２算出部４３と、遅延時間Δｔに基づいて核心部温度Ｔｃに対応づけられた時刻を補正する補正部４４とから構成される。

通信回路６０は、演算回路４０によって補正された温度の時系列データを外部に出力したり、エラーが発生したときにアラームを出力したりするためのＩ／Ｆ回路である。このような通信回路６０としては、有線でデータ等を出力する場合は、ＵＳＢその他のケーブルが接続できる出力回路となるが、例えば、Ｂｌｕｅｔｏｏｔｈ（登録商標）等に準拠した無線通信回路を用いることもできる。

シート状の基材８０は、熱流束センサ２０、血流センサ３０、演算回路４０、メモリ５０、通信回路６０および電池７０を載置するための土台として機能する他、これらの要素を電気的に接続する図示しない配線を備えている。生体内温度測定装置１を生体の表皮上に載置することを考えると、シート状の基材８０には、変形可能なフレキシブル基板を用いることが望ましい。
また、図３に示すように、シート状の基材８０の一部には開口８２、８３が設けられており、熱流束センサ２０および血流センサ３０は、これらの開口８２、８３からそれぞれ生体の表皮に接するように、基材８０に載置される。

血流センサ３０によって熱流束センサ２０近傍の体表部の血流量を測定するためには、仮に熱流束センサ２０の熱抵抗体２０ｒが例えば円盤状に形成されていたとすると、熱抵抗体２０ｒの直径Ｒの２倍程度の領域が核心部温度Ｔｃの計測に影響を与えるため、図４に示すように、熱流束センサ２０の熱抵抗体２０ｒの直径Ｒの２倍、すなわち平面視で熱抵抗体２０ｒを中心に直径２Ｒ程度の領域内に血流センサ３０を設置する。１個の熱流束センサ２０に対して血流センサ３０を１個または複数個設けることができるが、本実施の形態においては、１個の熱流束センサ２０に対して血流センサ３０を１個設けるものとする。

［生体内温度測定装置の測定原理］
次に、本実施の形態に係る生体内温度測定装置の測定原理について図５ないし図８を参照して説明する。
図５に示すように、生体９０においては、表皮から皮下組織の深さ方向に一定の深さを超えたところに、外気温の変化等に左右されない温度、すなわち、核心部温度の領域が存在し、通常は、核心部温度Ｔｃより表皮温度Ｔｓの方が低く、核心部から表皮に向かって温度勾配が生じている。

一般に、表皮温度Ｔｓの時間温度変化は、図６に示すように、核心部からの熱流束ｑ_skinと空気中への熱流束ｑ_airとから、次の式（１）のように表される。

ここで、ρは皮下組織の密度、Ａは熱流束が貫く面の面積、ｃは皮下組織の熱容量、Ｌは表皮から核心部温度の温度領域までの深さ（距離）である。

また、表皮から空気中への熱流束ｑ_airおよび核心部から表皮への熱流束ｑ_skinは、それぞれ比例定数ｈ（空気中への熱放射係数、熱蒸散係数、熱伝達係数を含む）および皮下組織の熱伝導率ｋを用いて、次の式（２）のように表される。

これを式（１）式に代入すると、次の（３）式のようになる。

この（３）式より、表皮温度Ｔｓは時定数τ＝（ρＬ²ｃ）／（ｈＬ＋ｋ）を持つため、核心部温度Ｔｃが表皮温度Ｔｓに反映されるまでの時間は、表皮から核心部温度の温度領域までの深さＬの２乗に比例することがわかる。

上述したように、表皮から核心部温度の温度領域までのみかけの深さＬは血流量に依存することが知られており（図１５参照。）、生体９０においては、血流によって時々刻々と表皮から核心部温度の温度領域までのみかけの深さが変化している。これをモデル化するため、図７に示すように、血管からの熱流束ｑ_bloodを新たに定義して（１）式に追加し同様に解くと、血流量と時定数τとは図８で表される関係となり、このように時定数τの変化は血流量の変化の関数となる。

以上のように、表皮から核心部温度の温度領域までの深さＬが変化すると時定数τが変化すること、および血流量が変化すると時定数τが変化することから、血流量の変化によって表皮から核心部温度の温度領域までのみかけの深さＬが変化することがわかる。

なお、熱流束センサ２０を用いた場合、２つの温度センサ２０ｕ、２０ｓで計測された空気温度Ｔａおよび表皮温度Ｔｓの差ΔＴと熱抵抗体２０ｒの熱伝導率λと厚さＬ_rから、熱抵抗体２０ｒを貫く熱流束ｑ_rの大きさ（この熱流束ｑ_rは、図６または図７における空気中への熱流束ｑ_airに相当する。）を求めることができる。
核心部から表皮に向かって皮下組織を貫く熱流束の大きさと熱抵抗体２０ｒを貫く熱流束の大きさが等しいと仮定すれば、演算回路４０の第１算出部４１（図１参照。）は、皮下組織の熱抵抗Ｒｘ、熱抵抗体２０ｒの熱抵抗Ｒｒ、表皮温度Ｔｓ、空気温度Ｔａから、次の（４）式に基づいて核心部温度Ｔｃを算出することができる。

［生体内温度測定装置の測定方法］
次に、本実施の形態に係る生体内温度測定装置の動作について、図９を参照して説明する。
なお、メモリ５０には、事前に用意された血流量と時定数τとの関係ｆ（vblood）（図８参照。）を表す関係式f(vblood)またはテーブルが記憶されているものとする。

まず、熱流束センサ２０の２つの温度センサ２０ｕ、２０ｓを用いて温度を測定するとともに（ステップＳ１０）、血流センサ３０を用いて体表部の血流量を測定する（ステップＳ２０）。この動作を複数回繰り返したところで空気温度Ｔａおよび表皮温度Ｔｓの変動が所定の範囲内に収まっているか否かを判断し（ステップＳ３０）、収まっていなければ（ステップＳ３０：Ｎｏ）熱流束は定常状態にないと判断して、ステップＳ１０に戻り、ステップＳ１０〜Ｓ３０を繰り返す。

一方、空気温度Ｔａおよび表皮温度Ｔｓの変動が所定の範囲内に収まっていれば（ステップＳ３０：Ｙｅｓ）、熱流束は定常状態にあると判断して、表皮温度Ｔｓおよび空気温度Ｔａから核心部温度Ｔｃを推定し（ステップＳ４０）、遅延時間Δｔを算出する（ステップＳ５０）。遅延時間Δｔの算出（ステップＳ５０）では、メモリ５０に記憶した血流量と時定数τとの関係ｆ（vblood）を参照して、血流センサ３０によって検出された体表部の血流量に対応した時定数τを求め、遅延時間Δｔを算出する。遅延時間Δｔは、事前に決定した定数αを用いて、下記の（５）式によって求めることができる。

Δｔ＝α×ｆ（vblood） （５）
ここで定数αは、時定数τより求められる周波数応答における位相遅延から算出することができる。

遅延時間Δｔが算出されたら、予め定めた閾値と比較して（ステップＳ６０）、算出された遅延時間Δｔがその閾値より大きい場合（ステップＳ６０：Ｎｏ）には、核心部温度の変化に追従できず、核心部温度Ｔｃはもとより核心部温度Ｔｃの変化を捉えることができないので、核心部温度の推定が不可能であったと判定して、アラームの発報など、所定のエラー処理を実行する（ステップＳ８０）。

一方、算出された遅延時間Δｔがその閾値以下だった場合（ステップＳ６０：Ｙｅｓ）は、その時点で推定された核心部温度Ｔｃの時刻を遅延時間Δｔの分だけ時間的にシフトして、時系列データを補正する（ステップＳ７０）。
以上のステップを終了の指示があるまで（ステップＳ９０：Ｙｅｓ）繰り返す。

以上のようにして補正された核心部温度Ｔｃの時系列データの一例を図１０に示す。図１０において、実線は実際の核心部温度を表し、〇印と□印は、それぞれ推定された補正前の核心部温度と補正後の核心部温度とを表す。血流量の変化に伴う遅延時間Δｔを考慮しない場合（〇印）は、実際の核心部温度Ｔｃと比較して時間軸上でずれが生じている。
これに対し、血流量に応じて算出された遅延時間Δtを用いて補正した時系列データでは、実際の核心部温度Ｔｃとの時間軸上のずれが低減される。したがって、本実施の形態に係る生体内温度測定装置によれば、核心部温度の変動をより正確に把握することができる。

例えば、本実施の形態に係る生体内温度測定装置において算出される補正量、すなわち遅延時間Δｔは、血流量によって変動するため、補正された核心部温度Ｔｃの時系列データからは、図１１に示すように、従来失われていた核心部温度Ｔｃの変動の周波数情報も復元することができる。これによって体内リズムとしての核心部温度Ｔｃの変動をより精確に把握することが出来る。

［第２の実施の形態］
次に図１２および図１３を参照して、本発明の第２の実施の形態について説明する。なお、上述した第１の実施の形態に係る生体内温度測定装置１と共通する構成要素については同一の符号を用い、その詳細な説明は省略する。

第１の実施の形態に係る生体内温度測定装置１は、１個の熱流束センサ２０に対して１個の血流センサ３０を用いたが、第２の実施の形態に係る生体内温度測定装置１ａは、１個の熱流束センサ２０に対して複数の血流センサを備えている。具体的には、第２の実施の形態に係る生体内温度測定装置１ａは、図１２に示すように、シート状の基材８０の上に、２つの温度センサ２０ｕ、２０ｓを含む熱流束センサ２０と、２つの血流センサ３０−１、３０−２と、演算回路４０ａと、メモリ５０と、外部とのＩ／Ｆ回路として機能する通信回路６０と、演算回路４０ａや通信回路６０等に電力を供給する電池７０とを備えている。

本実施の形態に係る生体内温度測定装置１ａは、１個の熱流束センサ２０に対して２個の血流センサ３０−１、３０−２を備えていることに伴い、演算回路４０ａは、第２の算出部４３ａにおいて、血流量の代表値として２個の血流センサ３０−１、３０−２でそれぞれ測定した血流量の平均値を求め、メモリ５０に記憶された関係を用いて、血流量の平均値に対応する時定数または遅延時間Δｔを求めるように構成されている。

２個の血流センサ３０−１、３０−２は、熱流束センサ２０の熱抵抗体２０ｒの近傍に配置される。このとき、熱抵抗体２０ｒの直径Ｒの２倍程度の領域が核心部温度Ｔｃの計測に影響を与えることを考慮すれば、図１３に示すように、熱流束センサ２０の熱抵抗体２０ｒの直径Ｒの２倍、すなわち平面視で熱抵抗体２０ｒを中心に直径２Ｒ程度の領域内に２個の血流センサ３０−１、３０−２を設置する。

なお、本実施の形態においては、２個の血流センサ３０−１、３０−２を、図１３に示すように、平面視で熱流束センサ２０の熱抵抗体２０ｒを挟んで線対象の位置に配置する。

このように複数の血流センサを用いることによって熱流束センサの近傍の血流量を精度よく計測することができるので、遅延時間Δｔの推定精度が上がり、核心部温度の変動をより正確に把握することが可能となる。

なお、本実施の形態においては、２個の血流センサを平面視で熱流束センサ２０の熱抵抗体２０ｒを挟んで線対象の位置に配置しているが、これに限定されるものではなく、例えば、図１４に示すように、２個の血流センサ３０−１、３０−２を、平面視で熱流束センサ２０の熱抵抗体２０ｒの中心を通り互いに直交する２つの線上に位置するように配置してもよい。

また、本実施の形態においては、２個の血流センサ３０−１、３０−２を用いることとしたが、３個以上の血流センサを用いてもよいことは言うまでもない。
また、複数の血流センサでそれぞれ測定した血流量の代表値としては、これらの血流量の平均値を用いることができるが、平均値に代えて、最大値や最小値、または３個以上の血流センサを用いた場合は、これらの加えて中央値等を用いてもよい。

１、１ａ…生体内温度測定装置、２０…熱流束センサ、２０ｕ、２０ｓ…温度センサ、２０ｒ…熱抵抗体、３０…血流センサ、４０…演算回路、５０…メモリ、６０…通信回路、７０…電池、８０…基材、９０…生体。

Claims (5)

1. 生体の表皮温度を測定する温度センサと、
   前記生体の体表部から放出される熱流束の大きさを測定する熱流束センサと、
   前記熱流束センサの近傍の血流量を測定する血流センサと、
   前記熱流束センサの近傍の血流量と前記生体の核心部温度の変動が表皮温度に反映されるまでの遅延時間に関するパラメータとの関係を記憶した記憶部と、
   前記温度センサで測定された表皮温度と前記熱流束センサで測定された熱流束の大きさから前記生体の核心部温度を算出するとともに、前記熱流束センサの近傍の血流量と前記記憶部に記憶された前記関係とに基づいて前記遅延時間を算出し、前記遅延時間に基づいて前記生体の核心部温度に対応づけられた時刻を補正するように構成された演算回路と
   を備える生体内温度測定装置。
2. 請求項１に記載された生体内温度測定装置において、
   前記演算回路は、
   前記温度センサで測定された表皮温度と前記熱流束センサで測定された熱流束の大きさから前記生体の核心部温度を算出する第１算出部と、
   前記熱流束センサの近傍の血流量と前記記憶部に記憶された前記関係とに基づいて前記遅延時間を算出する第２算出部と、
   前記遅延時間に基づいて前記生体の核心部温度に対応づけられた時刻を補正する補正部と
   を有することを特徴とする生体内温度測定装置。
3. 請求項１または２に記載された生体内温度測定装置において、
   前記パラメータは、前記生体の核心部温度の変動に対する表皮温度の変動の時定数または前記遅延時間であることを特徴とする生体内温度測定装置。
4. 請求項１〜３のいずれか１つに記載された生体内温度測定装置において、
   二以上の前記血流センサを備え、
   前記演算回路は、二以上の前記血流センサでそれぞれ測定した血流量の代表値を求め、前記血流量の代表値と前記記憶部に記憶された前記関係とに基づいて前記遅延時間を算出する
   を有することを特徴とする生体内温度測定装置。
5. 生体の表皮温度と前記生体の体表部から放出される熱流束の大きさとを測定するステップと、
   前記体表部の血流量を測定するステップと、
   測定された前記表皮温度と前記熱流束の大きさとから前記生体の核心部温度を算出するとともに、予め用意した、前記血流量と前記生体の核心部温度の変動が表皮温度に反映されるまでの遅延時間に関するパラメータとの関係と、前記体表部の血流量と、に基づいて前記遅延時間を算出し、前記遅延時間に基づいて前記生体の核心部温度に対応づけられた時刻を補正するステップと、を有する生体内温度測定方法。