濃度定量装置、濃度定量方法及びプログラム
【課題】真皮層以外の層によるノイズの影響を軽減して、真皮層に含まれるグルコースの濃度を算出する。
【解決手段】照射部１０４は、皮膚に短時間パルス光を照射し、受光部１０５は、皮膚によって後方散乱した光を受光し、計測光強度取得部１０６は、時刻ｔにおける受光部１０５の受光強度を取得する。光吸収係数算出部１０９は、計測光強度取得部１０６が取得した光強度と、光路長分布記憶部１０２が記憶する時刻ｔにおける皮膚モデルの各層の光路長と、時間分解波形記憶部１０３が記憶する時刻ｔにおける皮膚モデルの光強度とに基づいて、皮膚の真皮層の光吸収係数を算出し、濃度算出部１１０は、光吸収係数算出部１０９が算出した真皮層の光吸収係数に基づいて、真皮層のグルコース濃度を算出する。受光部１０５は複数用意されており、所望の層で後方散乱が生じるような照射受光間隔を有するものが用いられる。
【解決手段】照射部１０４は、皮膚に短時間パルス光を照射し、受光部１０５は、皮膚によって後方散乱した光を受光し、計測光強度取得部１０６は、時刻ｔにおける受光部１０５の受光強度を取得する。光吸収係数算出部１０９は、計測光強度取得部１０６が取得した光強度と、光路長分布記憶部１０２が記憶する時刻ｔにおける皮膚モデルの各層の光路長と、時間分解波形記憶部１０３が記憶する時刻ｔにおける皮膚モデルの光強度とに基づいて、皮膚の真皮層の光吸収係数を算出し、濃度算出部１１０は、光吸収係数算出部１０９が算出した真皮層の光吸収係数に基づいて、真皮層のグルコース濃度を算出する。受光部１０５は複数用意されており、所望の層で後方散乱が生じるような照射受光間隔を有するものが用いられる。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明のいくつかの様態は、複数の光散乱媒質の層から形成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置、濃度定量方法及びプログラムに関する。
【背景技術】
【０００２】
従来、血糖値の測定は、指先などから採血を行い、血中のグルコースに対する酵素活性を測ることで行っていた。しかし、このような血糖値の測定方法は、指先などから血液を採取して測定しなければならず、採血に手間と痛みを伴うことや、血液を付着させる測定チップが必要なことから、採血を必要としない非侵襲型の血糖値の測定方法が望まれている。
【０００３】
そこで、皮膚に近赤外光を照射し、その光吸収量からグルコースの濃度を求める方法が検討されている（例えば、特許文献１を参照）。具体的には、予めグルコース濃度と照射する光の波長と光の吸収量との関係を示す検量線を作成しておき、モノクロメーター等を用いてある波長域を走査し、その波長域の各波長に対する吸収量を求め、当該波長及び吸収量と検量線とを比較することでグルコース濃度を算出する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第３９３１６３８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、従来の非侵襲血糖値測定方法は、光の入出射間距離を定めることによって、真皮層の近赤外吸収スペクトルを測定するため、測定した吸収スペクトルに、真皮層の吸収スペクトルのみならず表皮層や皮下組織層の吸収スペクトルも含まれ、観測される吸収係数の変化には表皮層や皮下組織層によるノイズが含まれてしまうという問題があった。
本発明は上記の点に鑑みてなされたものであり、その目的は、目的の層以外の層によるノイズの影響を軽減する濃度定量装置、濃度定量方法及びプログラムを提供することにある。
【課題を解決するための手段】
【０００６】
本発明のいくつかの態様は上記の課題を解決するためになされたものであり、複数の光散乱媒質の層から形成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置であって、前記観測対象に短時間パルス光を照射する照射手段と、前記短時間パルス光が前記観測対象によって後方散乱した光を受光する複数の受光部を有する受光手段と、前記複数の受光部のうち、前記短時間パルス光が前記任意の層によって後方散乱した光を受光する特定の受光部を選択する選択手段と、前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記特定の受光手段が受光した光の強度を取得する光強度取得手段と、前記観測対象に対して照射する短時間パルス光の、前記照射手段から前記特定の受光部に至る光の伝搬経路上に配置された複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、前記観測対象に対して照射し前記特定の受光部において受光する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段と、前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手段と、前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手段と、前記光強度取得手段が取得した光強度と前記光路長取得手段が取得した前記複数の光散乱媒質の層の各々の層の光路長と前記光強度モデル取得手段が取得した光強度モデルとに基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手段と、前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手段と、を備えることを特徴とする。
これにより、受光した光の時間分解波形から、任意の層の光吸収計数を選択的に算出することができる。そのため、算出した光吸収計数に基づいて目的成分の濃度を算出することにより、他の層によるノイズの影響を軽減し、精度の高い濃度の定量を行うことができる。
ここで、短時間パルス光とは、パルス幅が１０ｐｓｅｃ程度かそれ以下のパルス光を意味する。短時間パルス光として０．１ｐｓｅｃから１０ｐｓｅｃの範囲のパルス幅を持つパルス光を用いても良い。
【０００７】
また、本発明のいくつかの態様は、前記光強度取得手段は、前記観測対象の層の数ｎ以上となる複数の時刻ｔ１〜ｔｍにおける光強度を取得し、前記光吸収係数算出手段は、下記式（１）に基づいて任意の層の光吸収係数を算出する、ことを特徴とする。
【０００８】
【数１】
【０００９】
但し、ｌｎ（Ａ）は、Ａの自然対数を示し、Ｉ（ｔ）は、前記受光手段が時刻ｔにおいて受光した光強度を示し、Ｎ（ｔ）は、前記短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度を示し、Ｌｉ（ｔ）は、前記伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長を示し、μｉは、第ｉ層の光吸収係数を示す。
【００１０】
また、本発明のいくつかの態様は、前記照射手段は、複数の波長１〜ｑの光を照射し、前記光吸収係数算出手段は、前記任意の層における光吸収係数を前記照射手段が照射した複数の波長毎に算出し、前記濃度算出手段は、下記式（２）に基づいて前記任意の層における前記目的成分の濃度を算出する、ことを特徴とする。
【００１１】
【数２】
【００１２】
但し、μａ（ｉ）は、前記任意の層である第ａ層における波長ｉの光吸収係数を示し、ｇｊは、前記観測対象を形成する第ｊ成分のモル濃度を示し、εｊ（ｉ）は、第ｊ成分の波長ｉに対する光吸収係数を示し、ｐは、前記観測対象を形成する主成分の個数を示し、ｑは、照射手段が照射する波長の種類数を示す。
【００１３】
また、本発明のいくつかの態様は、前記照射手段から前記観測対象に短時間パルス光を照射する位置を照射位置、前記短時間パルス光が前記観測対象によって後方散乱した光が前記観測対象から前記受光部に向けて出射する位置を受光位置、前記照射位置と前記受光位置との間隔を照射受光間隔、としたときに、前記受光手段は、互いに照射受光間隔の異なる複数の受光部を有し、前記選択手段は、前記照射受光間隔に応じて定まる前記観測対象の内部への短時間パルス光の到達深さに基づいて、前記任意の層に短時間パルス光が伝搬可能な照射受光間隔を有する受光部を選択することを特徴とする。
【００１４】
また、本発明のいくつかの態様は、前記照射手段は、短時間パルス光を前記観測対象の表面に伝送する照射用光ファイバーを有し、前記受光部は、前記観測対象によって後方散乱した光を伝送する受光用光ファイバーを有し、前記照射用光ファイバーと前記受光用光ファイバーは、両者の光ファイバーコアの中心間隔を所定の照射受光間隔だけ離間して固定するプローブ装置に装着されており、前記プローブ装置の先端部に露出した前記照射用光ファイバーの先端部と前記受光用光ファイバーの先端部とを前記観測対象の表面に接触させることによって、前記観測対象への短時間パルスの照射処理と、前記短時間パルス光が前記観測対象によって後方散乱した光の受光処理とを行うことを特徴とする。
【００１５】
また、本発明のいくつかの態様は、前記受光部は、互いに等しい照射受光間隔で配置された複数の受光用光ファイバーを有することを特徴とする。
【００１６】
また、本発明のいくつかの態様は、前記受光部は、互いに等しい照射受光間隔で配置された複数の受光用光ファイバーによって伝送された光を同一受光面上に集光する集光素子を有することを特徴とする。
【００１７】
また、本発明のいくつかの態様は、前記照射用光ファイバーと前記受光用光ファイバーは、前記照射手段が照射する前記複数の波長１〜ｑの光の波長分散を補償する分散補償型シングルモード光ファイバーであることを特徴とする。
【００１８】
また、本発明のいくつかの態様は、前記照射用光ファイバーと前記受光用光ファイバーは、前記照射手段が照射する前記複数の波長１〜ｑの光の波長分散に伴う群遅延時間差が、前記複数の光散乱媒質層のうち最も表面側の層の伝搬光路長分布のピークに対応する伝搬時間よりも短いことを特徴とする。
【００１９】
また、本発明のいくつかの態様は、複数の光散乱媒質の層から形成される観測対象に対して短時間パルス光を照射する照射手段と、前記短時間パルス光が前記観測対象によって後方散乱した光を受光する複数の受光部を有する受光手段と、前記複数の受光部のうち、前記短時間パルス光が前記任意の層によって後方散乱した光を受光する特定の受光部を選択する選択手段と、前記観測対象に対して照射する短時間パルス光の、前記照射手段から前記特定の受光部に至る光の伝搬経路上に配置された複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、前記観測対象に対して照射し前記特定の受光部において受光する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段とを備え、前記観測対象のうち任意の層における目的成分の濃度を定量する濃度定量装置を用いた濃度定量方法であって、光強度取得手段は、前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記特定の受光部が受光した光の強度を取得し、光路長取得手段は、前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得し、光強度モデル取得手段は、前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得し、光吸収係数算出手段は、前記光強度取得手段が取得した光強度と前記光路長取得手段が取得した前記複数の光散乱媒質の層の各々の層の光路長と前記光強度モデル取得手段が取得した光強度モデルとに基づいて、前記任意の層の光吸収係数を算出し、濃度算出手段は、前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する、ことを特徴とする。
【００２０】
また、本発明のいくつかの態様は、複数の光散乱媒質の層から形成される観測対象に対して短時間パルス光を照射する照射手段と、前記短時間パルス光が前記観測対象によって後方散乱した光を受光する複数の受光部を有する受光手段と、前記複数の受光部のうち、前記短時間パルス光が前記任意の層によって後方散乱した光を受光する特定の受光部を選択する選択手段と、前記観測対象に対して照射する短時間パルス光の、前記照射手段から前記特定の受光部に至る光の伝搬経路上に配置された複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、前記観測対象に対して照射し前記特定の受光部において受光する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段とを備え、前記観測対象のうち任意の層における目的成分の濃度を定量する濃度定量装置を、前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記特定の受光部が受光した光の強度を取得する光強度取得手段、前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手段、前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手段、前記光強度取得手段が取得した光強度と前記光路長取得手段が取得した前記複数の光散乱媒質の層の各々の層の光路長と前記光強度モデル取得手段が取得した光強度モデルとに基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手段、前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手段、として動作させるためのプログラムである。
【図面の簡単な説明】
【００２１】
【図１】本発明による血糖値測定装置の構成を示す概略ブロック図である。
【図２】シミュレーション部が算出した各層の伝搬光路長分布を示すグラフである。
【図３】シミュレーション部が算出した時間分解波形を示すグラフである。
【図４】皮膚の主成分の吸収スペクトルを示すグラフである。
【図５】皮膚組織の断面図である。
【図６】プローブ装置の断面図である。
【図７】血糖値測定装置が血糖値を測定する動作を示す第１のフローチャートである。
【図８】血糖値測定装置が血糖値を測定する動作を示す第２のフローチャートである。
【図９】光路長取得部が取得する所定時間の光路長を示す図である。
【発明を実施するための形態】
【実施例１】
【００２２】
以下、図面を参照しながら本発明の第１の実施形態について詳しく説明する。図１は、本発明の第１の実施形態による血糖値測定装置の構成を示す概略ブロック図である。
血糖値測定装置１００（濃度定量装置）は、シミュレーション部１０１、光路長分布記憶部１０２（光路長分布記憶手段）、時間分解波形記憶部１０３（時間分解波形記憶手段）、照射部１０４（照射手段）、受光部１０５（受光手段）、計測光強度取得部１０６（光強度取得手段）、光路長取得部１０７（光路長取得手段）、無吸収時光強度取得部１０８（光強度モデル取得手段）、光吸収係数算出部１０９（光吸収係数算出手段）、濃度算出部１１０（濃度算出手段）を備える。
血糖値測定装置１００は、皮膚（観測対象）の真皮層（任意の層）に含まれるグルコース（目的成分）の濃度を測定する。
【００２３】
シミュレーション部１０１は、光吸収係数がゼロの皮膚モデルに対して光を照射するシミュレーションを行う。
光路長分布記憶部１０２は、光吸収係数がゼロの皮膚モデルの伝搬光路長分布を記憶する。
時間分解波形記憶部１０３は、光吸収係数がゼロの皮膚モデルの時間分解波形を記憶する。
照射部１０４は、皮膚に対して短時間パルス光を照射する。
受光部１０５は、短時間パルス光が皮膚によって後方散乱した光を受光する。
計測光強度取得部１０６は、受光部１０５が受光した光のある時刻における光強度を取得する。
光路長取得部１０７は、光路長分布記憶部１０２からある時刻における光路長を取得する。
無吸収時光強度取得部１０８は、時間分解波形記憶部１０３からある時刻における光強度を取得する。
光吸収係数算出部１０９は、短時間パルス光を照射した皮膚の真皮層における光吸収係数を算出する。
濃度算出部１１０は、真皮層に含まれるグルコースの濃度を算出する。
ここで、短時間パルス光とは、パルス幅が１０ｐｓｅｃ程度かそれ以下のパルス光を意味する。短時間パルス光として０．１ｐｓｅｃから１０ｐｓｅｃの範囲のパルス幅を持つパルス光を用いても良い。
【００２４】
そして、血糖値測定装置１００において、照射部１０４は、皮膚に短時間パルス光を照射し、受光部１０５は、短時間パルス光が皮膚によって後方散乱した光を受光し、計測光強度取得部１０６は、時刻ｔにおいて受光部１０５が受光した光の強度を取得する。次に、光路長取得部１０７は、光路長分布記憶部１０２から、皮膚モデルにおける伝搬光路長分布の時刻ｔにおける皮膚の各層の光路長を取得し、無吸収時光強度取得部１０８は、時間分解波形記憶部１０３から、皮膚モデルにおける短時間パルス光の時間分解波形の時刻ｔにおける光の強度を取得する。
【００２５】
次に、光吸収係数算出部１０９は、計測光強度取得部１０６が取得した光強度と光路長取得部１０７が取得した皮膚の各層の光路長と無吸収時光強度取得部１０８が取得した光強度とに基づいて、皮膚の真皮層の光吸収係数を算出し、濃度算出部１１０は、光吸収係数算出部１０９が算出した光吸収係数に基づいて、真皮層におけるグルコースの濃度を算出する。
これにより、真皮層以外の層によるノイズの影響を軽減して、真皮層に含まれるグルコースの濃度を算出することができる。
【００２６】
次に、血糖値測定装置１００の動作を説明する。
血糖値測定装置１００は、血糖値を測定する前に、予め皮膚モデルの各層における伝搬光路長分布と時間分解波形とを算出しておく必要がある。
まず、皮膚モデルの伝搬光路長分布及び時間分解波形の算出方法を説明する。
初めに、シミュレーション部１０１は、皮膚モデルを生成する。皮膚モデルの生成は、皮膚の各層の光散乱係数、光吸収係数及び厚みを決定することで行う。ここで、皮膚の各層の散乱係数及び厚みは、個体による差が少ないため、予めサンプルを取ることなどによって決定すると良い。なお、表皮層の厚みは略０．３ｍｍ、真皮層の厚みは略１．２ｍｍ、皮下組織層の厚みは略３．０ｍｍである。
また、ここで用いる皮膚モデルの光吸収係数はゼロとする。これは、当該皮膚モデルを用いて光吸収量を算出するためである。
【００２７】
シミュレーション部１０１は、皮膚モデルを生成すると、当該皮膚モデルに光を照射するシミュレーションを行う。このとき、照射部１０４の位置と受光部１０５の位置との間の距離を決定しておく必要がある。シミュレーションは、モンテカルロ法を用いて行うと良い。モンテカルロ法によるシミュレーションは、例えば以下のように行われる。
まず、シミュレーション部１０１は、照射する光のモデルを光子（光束）とし、当該光子を皮膚モデルに照射する計算を行う。皮膚モデルに照射された光子は、皮膚モデル内を移動する。このとき、光子は、次に進む点までの距離Ｌ及び方向θを乱数Ｒによって決定する。シミュレーション部１０１は、光子が次に進む点までの距離Ｌの計算を、式（３）によって行う。
【００２８】
【数３】
【００２９】
但し、ｌｎ（Ａ）は、Ａの自然対数を示す。また、μｓは、皮膚モデルの第ｓ層（表皮層、真皮層、皮下組織層の何れか）の散乱係数を示す。
また、シミュレーション部１０１は、光子が次に進む点までの方向θの計算を、式（４）によって行う。
【００３０】
【数４】
【００３１】
但し、ｇは、散乱角度の余弦（ｃｏｓ）の平均である非等方性パラメータを示し、皮膚の非等方性パラメータは、略０．９である。
シミュレーション部１０１は、上記式（３）、式（４）の計算を単位時間毎に繰り返すことにより、照射部１０４から受光部１０５までの光子の移動経路を算出することができる。シミュレーション部１０１は、複数の光子について移動距離の算出を行う。例えば、シミュレーション部１０１は、１０８個の光子について移動距離を算出する。
【００３２】
図２は、シミュレーション部が算出した各層の伝搬光路長分布を示すグラフである。
図２の横軸は光子の照射からの経過時間を示し、縦軸は光路長の対数表示を示している。シミュレーション部１０１は、受光部１０５に到達した光子の各々の移動経路を、移動経路が通過する層毎に分類する。そして、シミュレーション部１０１は、単位時間毎に到達した光子の移動経路の平均長を分類された層毎に算出することで、図２に示すような皮膚の各層の伝搬光路長分布を算出する。
【００３３】
図３は、シミュレーション部が算出した時間分解波形を示すグラフである。
図３の横軸は光子の照射からの経過時間を示し、縦軸は受光部１０５が検出した光子数を示している。シミュレーション部１０１は、単位時間毎に受光部１０５に到達した光子の個数を算出することで、図３に示すような皮膚モデルの時間分解波形を算出する。
上述したような処理により、シミュレーション部１０１は、複数の波長に対して、皮膚モデルの伝搬光路長分布及び時間分解波形を算出する。このとき、シミュレーション部１０１は、皮膚の主成分（水、たんぱく質、脂質、グルコース等）の吸収スペクトルの差が大きくなる波長について伝搬光路長分布及び時間分解波形を算出すると良い。
【００３４】
図４は、皮膚の主成分の吸収スペクトルを示すグラフである。
図４の横軸は照射する光の波長を示し、縦軸は吸収係数を示している。図４を参照すると、グルコースの吸収係数は、波長が１６００ｎｍのときに極大となり、水の吸収係数は、波長が１４５０ｎｍのときに極大となる。そのため、シミュレーション部１０１は、例えば１４００ｎｍ、１４５０ｎｍ、１５００ｎｍ、１６００ｎｍ、１６８０ｎｍ、１７２０ｎｍ、１７４０ｎｍといった皮膚の主成分の吸収スペクトルの差が大きくなる波長について伝搬光路長分布及び時間分解波形を算出すると良い。
【００３５】
シミュレーション部１０１は、複数の波長に対する皮膚モデルの伝搬光路長分布及び時間分解波形を算出すると、伝搬光路長分布の情報を光路長分布記憶部１０２に記憶させ、時間分解波形の情報を時間分解波形記憶部１０３に記憶させる。
【００３６】
ここで、図５に基づいて、観測対象である皮膚組織の構造を説明する。
図５に示すように、人間の皮膚組織は、表皮と、真皮と、皮下組織の３層によって形成されている。表皮は、最も外側にある厚さ０．２〜０．３mmの薄い層で、角質層、顆粒層、有棘層、底層等を含む。真皮は、表皮と皮下組織の間に存在する厚さ０．５〜２ｍｍの層で、真皮内には神経、毛根、皮脂腺や汗腺、毛包、血管、リンパ管が存在する。皮下組織は、真皮の下にある厚さ１〜３ｍｍの層で、大部分が皮下脂肪でできている。
【００３７】
真皮内には毛細血管等が発達しており、血中グルコースに応じた物質移動が速やかに起こり、血中グルコース濃度（血糖値）に対して真皮中のグルコース濃度は追随して変化すると考えられている。そのため、血糖値測定装置１００では、照射部１０４から皮膚組織の表面に光を照射し、皮膚組織内を透過、拡散して後方散乱された光を受光部１０５で検出することによって、真皮中のグルコース濃度を測定している。
【００３８】
照射部１０４は、光源１０４ｂと、光源１０４ｂから射出された光を表皮に向けて伝送する照射用光ファイバー１０４ａとを含む。受光部１０５は、光電変換素子１０５ｂと、皮膚組織内を透過、拡散して表皮から出射した光を光電変換素子１０５ｂに向けて伝送する受光用光ファイバー１０５ａとを含む。照射用光ファイバー１０４ａと受光用光ファイバー１０５ａは、光ファイバーコアの中心間隔が所定の照射受光間隔だけ離間した状態で、互いに光軸（コアの中心軸）が平行となるようにして図示略のプローブ装置に固定されている。
【００３９】
受光用光ファイバー１０５ａに到達する光の経路分布ＲはＵ字型の形状をしている。実際には、図示点線で示したような複数のＵ字型の経路が存在し、それぞれの経路における光の伝搬距離及び伝搬時間は互いに異なっている。皮膚組織におけるこのような経路分布の到達深さは、照射部１０４と受光部１０５との間隔、具体的には、照射用光ファイバー１０４ａと受光用光ファイバー１０５ａのコア中心間の距離である照射受光間隔Ｗによって変化する。したがって、照射部１０４と受光部１０５との照射受光間隔Ｗを適切に設定することにより、主として真皮を透過する光の吸収スペクトルを測定することができる。
【００４０】
図６（ａ）〜図６（ｃ）は、照射用光ファイバー１０４ａと受光用光ファイバー１０５ａの光軸に垂直な面で切った断面図である。なお、図６（ａ）〜図６（ｃ）において、符号Ａで示した光ファイバーは照射用光ファイバー１０４ａであり、符号Ｂ１〜Ｂ４で示した光ファイバーは受光用光ファイバー１０５ａである。
【００４１】
図６（ａ）のプローブ装置１３０は、左端に配置された１つの照射用光ファイバー１０４ａと、その右側に等間隔で配置された複数の受光用光ファイバー１０５ａとを含む。プローブ装置１３０では、照射用光ファイバー１０４ａ及び光源１０４ｂ（図５参照）が本発明の照射手段に該当し、その周囲に配置された複数の受光用光ファイバー１０５ａ及び光電変換素子１０５ｂ（図５参照）が本発明の受光手段に該当する。
【００４２】
照射用光ファイバー１０４ａと受光用光ファイバー１０５ａは、光軸に垂直な方向に一列に配置されている。プローブ装置１３０では、プローブ装置１３０の先端部に露出した照射用光ファイバー１０４ａの先端部と受光用光ファイバー１０５ａの先端部とを皮膚表面に接触させることによって、皮膚組織への短時間パルス光の照射処理と、短時間パルス光が皮膚組織によって後方散乱した光の受光処理とを行うようになっている。
【００４３】
複数の受光用光ファイバー１０５ａには、照射用光ファイバー１０４ａに対して第１の照射受光間隔で配置された第１の受光用光ファイバー１０５ａ（符号Ｂ１で示す）と、その右側に第１の照射受光間隔よりも大きい第２の照射受光間隔で配置された第２の受光用光ファイバー１０５ａ（符号Ｂ２で示す）と、その右側に第２の照射受光間隔よりも大きい第３の照射受光間隔で配置された第３の受光用光ファイバー１０５ａ（符号Ｂ３で示す）と、その右側に第３の照射受光間隔よりも大きい第４の照射受光間隔で配置された第４の受光用光ファイバー１０５ａ（符号Ｂ４で示す）と、が含まれる。
【００４４】
複数の受光用光ファイバー１０５ａは、皮膚組織への到達深さに応じて、１つ又は複数の受光用光ファイバー１０５ａが選択され使用される。例えば、表皮を透過する光のスペクトルを精度良く測定したい場合には、照射受光間隔が最も小さい第１の受光用光ファイバー１０５ａ（符号Ｂ１で示す）を選択し、真皮を透過する光のスペクトルを精度良く測定したい場合には、第２の照射受光間隔で配置される第２の受光用光ファイバー１０５ａ（符号Ｂ２で示す）と第３の照射受光間隔で配置される第３の受光用光ファイバー１０５ａ（符号Ｂ３で示す）とを選択するという方法が用いられる。また、皮膚組織の厚みは個人差があり、例えば大人と子供、男性と女性では、皮膚組織の厚みが異なる場合がある。そのため、同じく真皮の情報を得る場合であっても、表皮及び真皮の厚みが厚いと予想される被検体には、厚みが薄いと予想される被検体よりも照射受光間隔の大きい受光用光ファイバー１０５ａを用いることも考えられる。
【００４５】
図６（ａ）では、４本の受光用光ファイバー１０５ａのみを示したが、受光用光ファイバー１０５ａの数はこれに限らず、３本若しくは５本以上とすることも可能である。受光用光ファイバー１０５ａの数を増やせば、皮膚組織への到達深さを多段階で制御することができる。
【００４６】
図６（ｂ）のプローブ装置１３１は、中央に配置された１つの照射用光ファイバー１０４ａと、その両側に等間隔で配置された複数の受光用光ファイバー１０５ａとを含む。プローブ装置１３１では、照射用光ファイバー１０４ａ及び光源１０４ｂ（図５参照）が本発明の照射手段に該当し、その周囲に配置された複数の受光用光ファイバー１０５ａ及び光電変換素子１０５ｂ（図５参照）が本発明の受光手段に該当する。
【００４７】
照射用光ファイバー１０４ａと受光用光ファイバー１０５ａは、光軸に垂直な方向に一列に配置されている。プローブ装置１３１では、プローブ装置１３１の先端部に露出した照射用光ファイバー１０４ａの先端部と受光用光ファイバー１０５ａの先端部とを皮膚表面に接触させることによって、皮膚組織への短時間パルス光の照射処理と、短時間パルス光が皮膚組織によって後方散乱した光の受光処理とを行うようになっている。
【００４８】
複数の受光用光ファイバー１０５ａには、照射用光ファイバー１０４ａの両側に第１の照射受光間隔で配置された２つの第１の受光用光ファイバー１０５ａ（符号Ｂ１で示す）と、その両側に第１の照射受光間隔よりも大きい第２の照射受光間隔で配置された２つの第２の受光用光ファイバー１０５ａ（符号Ｂ２で示す）と、が含まれる。
【００４９】
複数の受光用光ファイバー１０５ａは、皮膚組織への到達深さに応じて、１つ又は複数の受光用光ファイバー１０５ａが選択され使用される。例えば、表皮を透過する光のスペクトルを精度良く測定したい場合には、照射受光間隔が最も小さい第１の受光用光ファイバー１０５ａ（符号Ｂ１で示す）を選択し、真皮を透過する光のスペクトルを精度良く測定したい場合には、第１の照射受光間隔で配置される第１の受光用光ファイバー１０５ａ（符号Ｂ１で示す）と第２の照射受光間隔で配置される第２の受光用光ファイバー１０５ａ（符号Ｂ２で示す）とを選択するという方法が用いられる。図６（ｂ）のプローブ装置１３１は、図６（ａ）のプローブ装置１３０に比べて、選択できる照射受光間隔が少ないが、１つの照射受光間隔に対して複数の受光用光ファイバー１０５ａが設けられているので、受光強度は大きくなる。そのため、ノイズの少ない測定が可能となる。
【００５０】
図６（ｂ）では、照射受光間隔の異なる２種類の受光用光ファイバー１０５ａを示したが、照射受光間隔は２種類に限らず、３種類以上とすることもできる。例えば、３種類の照射受光間隔で測定を行う場合には、第２の受光用光ファイバー１０５ａの両側に第２の照射受光間隔よりも大きい第３の照射受光間隔で配置された２つの第３の受光用光ファイバー１０５ａを設ければよい。照射受光間隔を４種類以上とする場合も同様である。
【００５１】
図６（ｃ）のプローブ装置１３２は、中央に配置された１つの照射用光ファイバー１０４ａと、その周囲に配置された複数の受光用光ファイバー１０５ａとを含む。プローブ装置１３２では、照射用光ファイバー１０４ａ及び光源１０４ｂ（図５参照）が本発明の照射手段に該当し、その周囲に配置された複数の受光用光ファイバー１０５ａ及び光電変換素子１０５ｂ（図５参照）が本発明の受光手段に該当する。
【００５２】
照射用光ファイバー１０４ａと受光用光ファイバー１０５ａは、光軸に垂直な平面内で最密充填されるように配置されている。すなわち、受光用光ファイバー１０５ａは、照射用光ファイバー１０４ａを中心に６０°の内角を持つように放射方向に並んで配置されている。プローブ装置１３２では、プローブ装置１３２の先端部に露出した照射用光ファイバー１０４ａの先端部と受光用光ファイバー１０５ａの先端部とを皮膚表面に接触させることによって、皮膚組織への短時間パルスの照射処理と、短時間パルス光が皮膚組織によって後方散乱した光の受光処理とを行うようになっている。
【００５３】
複数の受光用光ファイバー１０５ａには、照射用光ファイバー１０４ａの周囲に第１の照射受光間隔で配置された６つの第１の受光用光ファイバー１０５ａ（符号Ｂ１で示す）と、その周囲に第１の照射受光間隔よりも大きい第２の照射受光間隔で配置された６つの第２の受光用光ファイバー１０５ａ（符号Ｂ２で示す）及び第２の照射受光間隔よりも大きい第３の照射受光間隔で配置された６つの第３の受光用光ファイバー１０５ａ（符号Ｂ３で示す）と、が含まれる。
【００５４】
複数の受光用光ファイバー１０５ａは、皮膚組織への到達深さに応じて、１つ又は複数の受光用光ファイバー１０５ａが選択され使用される。図６（ｃ）のプローブ装置１３２は、図６（ｂ）のプローブ装置１３１に比べて、１つの照射受光間隔に対して設けられる受光用光ファイバー１０５ａの数が多いため、よりノイズの少ない測定が可能となる。
【００５５】
なお、図６（ｂ）及び図６（ｃ）では、１つの照射受光間隔に対して複数の受光用光ファイバー１０５ａが設けられている。この場合、複数の受光用光ファイバー１０５ａで伝送される光は、受光用光ファイバー毎に設けられた光電変換素子によって電気信号に変換しても良いし、同一照射受光間隔で配置された複数の受光用光ファイバー１０５ａで伝送した光を集光素子によって同一光電変換素子の同一受光面上に集光し、これらを一括して電気信号に変換してもよい。後者の方法を用いれば、光電変換素子を複数の受光用光ファイバー１０５ａに対して共通化することができるので、部材コストを低減することができる。
【００５６】
受光用光ファイバー毎に１つの光電変換素子を設ける場合には、１つ受光用光ファイバーと１つの光電変換素子によって１つの受光部が構成される。複数の受光用光ファイバーに対して１つの光電変換素子が設けられる場合には、複数の受光用光ファイバーと１つの光電変換素子によって１つの受光部が構成される。
【００５７】
照射用光ファイバー１０４ａと受光用光ファイバー１０５ａは、照射部１０４が照射する複数の波長１〜ｑの光の波長分散を補償する分散補償型シングルモード光ファイバーであることが望ましい。これにより、波長分散による測定誤差を低減することができる。
【００５８】
照射用光ファイバー１０４ａと受光用光ファイバー１０５ａは、分散補償型でないマルチモード光ファイバーを用いることもできる。マルチモード光ファイバーを用いることで、シングルモード光ファイバーを用いる場合に比べて受光面積（コア径）を大きくすることができ、ノイズの少ない測定を行うことが可能となる。ただし、波長分散による測定誤差が大きくなるため、照射用光ファイバー１０４ａと受光用光ファイバー１０５ａにおいては、照射部１０４が照射する複数の波長１〜ｑの光の波長分散に伴う群遅延時間差が、皮膚組織に含まれる複数の光散乱媒質層のうち最も表面側の層である表皮の伝搬光路長分布のピークに対応する伝搬時間よりも短くなるように、各々の光ファイバーの長さが設定されることが望ましい。
【００５９】
例えば、ステップインデックス型マルチモード光ファイバーを用いる場合を考える。ステップインデックス型マルチモード光ファイバの群遅延時間差は下記の式（５）で表せる。
【００６０】
【数５】
【００６１】
但し、Ｔｏは光ファイバの光軸上を進む最も早い光の群遅延時間を示し、Δは光ファイバーのコアとクラッドとの比屈折率差（＝１％）を示し、θcは臨界角、ｎ１はコアの屈折率を示し、Lは光ファイバーの長さを示し、ｃは光速を示す。Ｔｏは図５で示した真皮を通る複数の経路のうち最も短い伝搬距離を有する経路、すなわち、短時間パルス光の到達深さが表皮と真皮の界面付近となる光の経路である。ΔＴmが真皮伝搬光の最短到達時間より短くなる様なマルチモードファイバにすることでシングルモードファイバと比べて信号光の光強度を大きくできる。
【００６２】
例えば、群遅延時間差を１．０ｐｓ以内で、光の強度を受光するための受光用光ファイバー１０５ａの長さは、コアの屈折率を１．５、光速を３×１０８ｍ／ｓとして、式（５）より、２０ｍｍと算出される。
【００６３】
マルチモードファイバ内での光の伝搬を幾何学的に説明すると、光軸に沿って（θ１＝０）進む光がｈだけ進む間に、光軸に対して角度θ１で進む光はｄ＝ｈ／（ｃｏｓθ１）だけ進む。スネルの法則からｃｏｓθ１＝ｎ２／ｎ１（但し、ｎ２はクラッドの屈折率を示し、ｎ１はコアの屈折率を示す）なので、ｄ＝ｈ（ｎ１／ｎ２）になる。
【００６４】
ここで、比屈折率差Δ（＝１−ｎ２／ｎ１）は１％なので、ｎ１／ｎ２＝１．０１である。よってｄ＝１．０１×ｈになる。ｈ＝２０ｍｍだからｄ＝２０．２ｍｍである。つまり、hが２０ｍｍのときのｈとｄの伝搬経路差はｄ−ｈ＝０．２ｍｍである。そして、屈折率ｎの物質中を進む光の速度ｖ＝空気中の光の速度／ｎなので、この場合ｖ＝３×１０８／１．５＝２×１０８ｍ／ｓとなる。よって、伝搬時間差ｔ＝伝搬経路差／光の速度という関係より伝搬経路差＝伝搬時間差ｔ×光の速度＝２０ｍｍになる。
【００６５】
次に、血糖値測定装置１００が血糖値を測定する動作について説明する。
図７は、血糖値測定装置が血糖値を測定する動作を示す第１のフローチャートである。
まず、ユーザが血糖値測定装置１００を皮膚にあてがい、測定開始スイッチ（図示せず）の押下等によって血糖値測定装置１００を動作させると、照射部１０４は、皮膚に対して波長λ１の短時間パルス光を照射する（ステップＳ１）。ここで、波長λ１は、シミュレーション部１０１が伝搬光路長分布及び時間分解波形を算出した複数の波長の中の１つである。
【００６６】
照射部１０４が短時間パルス光を照射すると、受光部１０５は、照射部１０４から照射され、皮膚によって後方散乱した光を受光する（ステップＳ２）。このとき、受光部１０５は、照射開始からの単位時間毎（例えば、１ピコ秒毎）の受光強度を内部メモリに登録しておく。
【００６７】
ここで、対象となる受光部１０５は、受光手段に設けられた複数の受光部の中から選択された特定の受光部である。受光手段に設けられる複数の受光部１０５は、それぞれ受光量に応じた電気信号を受光強度として図示略の内部メモリに記憶している。しかし、それらの受光強度のデータが全て光吸収係数算出部１０９によって取得されるわけではなく、特定の受光部１０５の内部メモリに記憶された受光強度のデータが光吸収係数算出部１０９によって選択的に取得される。すなわち、光吸収係数算出部１０９が本発明の選択手段に該当する。光吸収係数算出部１０９は、ユーザーがメニュー画面に従って直接的又は間接的に選択した受光部を選択しても良いし、予めデフォルトで設定された受光部を選択しても良い。前者の例としては、直接ユーザーが「受光部１」「受光部２」などの表示ボタンを操作することにより選択するものや、ユーザーが被検体の年齢、性別などをメニュー画面で選択し、年齢、性別などに関連付けられた受光部を間接的に選択するものなどが挙げられる。
【００６８】
受光部１０５が受光を完了すると、計測光強度取得部１０６は、光吸収係数算出部１０９によって選択される特定の受光部１０５の内部メモリに格納されている、異なる時刻ｔにおける受光強度Ｉ（ｔ）を皮膚の層の数と同じ数だけ取得する（ステップＳ３）。すなわち、計測光強度取得部１０６は、３つの異なる時刻ｔ１〜ｔ３における受光強度Ｉ（ｔ１）〜Ｉ（ｔ３）を取得する。ここで、皮膚の層の数と同じ数だけ受光強度を取得する理由は、後述する処理において、皮膚の各層の吸収係数を連立方程式によって算出するためである。
【００６９】
また、計測光強度取得部が光強度を取得する時刻ｔ１〜ｔ３は、皮膚の各層の伝搬光路長分布のピークとなる時刻であると良い。すなわち、照射部１０４が短時間パルス光を照射した時刻に、図２に示すグラフにおいて皮膚の各層の光路長が極大となる時間を加算した時刻の光強度をそれぞれ取得すると良い。
【００７０】
計測光強度取得部１０６が、受光強度Ｉ（ｔ１）〜Ｉ（ｔ３）を取得すると、光路長取得部１０７は、光路長分布記憶部１０２が記憶する波長λ１の伝搬光路長分布から、時刻ｔ１〜ｔ３における皮膚の各層の光路長Ｌ１（ｔ１）〜Ｌ１（ｔ３）、Ｌ２（ｔ１）〜Ｌ２（ｔ３）、Ｌ３（ｔ１）〜Ｌ３（ｔ３）を取得する（ステップＳ４）。
また、計測光強度取得部１０６が、受光強度Ｉ（ｔ１）〜Ｉ（ｔ３）を取得すると、無吸収時光強度取得部１０８は、時間分解波形記憶部１０３が記憶する波長λ１の時間分解波形から、時刻ｔ１〜ｔ３における無吸収時光強度Ｎ（ｔ１）〜Ｎ（ｔ３）を取得する（ステップＳ５）。
【００７１】
光路長取得部１０７が皮膚の各層の光路長を取得し、無吸収時光強度取得部１０８が無吸収時光強度を取得すると、光吸収係数算出部１０９は、式（６）に基づいて、皮膚の各層の光吸収係数μ１〜μ３を算出する（ステップＳ６）。ここで、光吸収係数μ１は、表皮層の光吸収係数を示し、光吸収係数μ２は、真皮層の光吸収係数を示し、光吸収係数μ３は、皮下組織層の光吸収係数を示す。
【００７２】
【数６】
【００７３】
但し、ｌｎ（Ａ）は、Ａの自然対数を示す。また、Ｉｉｎは、照射部１０４が照射した短時間パルス光の光強度を示す。また、Ｎｉｎは、シミュレーション部１０１が照射のシミュレーションを行った光子の個数を示す。
光吸収係数算出部１０９が皮膚の各層の光吸収係数μ１〜μ３を算出すると、光吸収係数算出部１０９は、求める成分の種類数αとして、xＣ2≧αを満足する波長の種類数（例えばα=3→x=3、α=4→x=4、α=5→x=4、）に対して光吸収係数μ１〜μ３を算出したか否かを判定する（ステップＳ７）。本実施形態では、皮膚の主成分を水、たんぱく質、脂質、グルコースの４種類として血糖値の測定を行うため、光吸収係数算出部１０９は、４種類の波長λ１〜λ４に対して光吸収係数μ１〜μ３を算出したか否かを判定する。ここで、波長λ１〜λ４は、シミュレーション部１０１が伝搬光路長分布及び時間分解波形を算出した複数の波長の中から選出する。
【００７４】
光吸収係数算出部１０９が、光吸収係数μ１〜μ３を算出していない波長λ１〜λ４があると判定した場合（ステップＳ７：ＮＯ）、ステップＳ１に戻り、まだ光吸収係数μ１〜μ３を算出していない波長λ１〜λ４の光吸収係数μ１〜μ３の算出を行う。
他方、光吸収係数算出部１０９が、波長λ１〜λ４の光吸収係数μ１〜μ３を算出していると判定した場合（ステップＳ７：ＹＥＳ）、濃度算出部１１０は、式（７）に基づいて真皮質に含まれるグルコースの濃度を算出する（ステップＳ８）。
【００７５】
【数７】
【００７６】
但し、μ２（１）〜μ２（４）は、真皮層における波長λ１〜λ４の光吸収係数を示す。また、ｇ１〜ｇ４は、真皮層におけるそれぞれ皮膚の主成分である水、たんぱく質、脂質、グルコースのモル濃度を示す。また、ε１（１）〜ε１（４）は、波長λ１〜λ４に対する水のモル吸光係数を示し、ε２（１）〜ε２（４）は、波長λ１〜λ４に対するたんぱく質のモル吸光係数を示し、ε３（１）〜ε３（４）は、波長λ１〜λ４に対する脂質のモル吸光係数を示し、ε４（１）〜ε４（４）は、波長λ１〜λ４に対するグルコースのモル吸光係数を示す。
つまり、式（７）のｇ４を算出することで、真皮質に含まれるグルコースのモル濃度を求めることができる。
【００７７】
ここで、式（７）によってグルコースのモル濃度を求めることができる理由を説明する。皮膚の散乱係数の波長依存性は小さいため、無吸収時光強度Ｎ（ｔ）及び光路長Ｌｎ（ｔ）の波長に対する変化は無視することができる。また、ベア・ランベルト（Beer-Lambert）の法則により、吸光度＝モル吸光係数×モル濃度で表すことができる。これにより、２波長で得られた時間分解計測より、無吸収時光強度Ｎ（ｔ）を消去することで、真皮層において得られる吸収計数差と皮膚を形成する各成分のモル吸光係数との関係式を示す式（７）を導くことができる。
【００７８】
このように、本実施形態によれば、短時間パルス光を照射し、所定の時刻において受光した光の強度に基づいてグルコースの濃度を定量する。これにより、所定の時刻において受光した光の中から、真皮層の吸収計数を選択的に算出することができる。これにより、特定の皮膚の層のグルコースの濃度を算出することができ、他の層によるノイズの影響を軽減し、精度の高い血糖値を算出することが可能になる。
【００７９】
なお、式（７）では、複数波長での光吸収係数の差分に基づいてグルコースのモル濃度を求めた。しかし、複数波長での光吸収係数の差分によらずに下式（８）、（９）に基づいて直接グルコースのモル濃度を求めることもできる。
【００８０】
式（８）は、各層における光吸収係数μ１〜μ３を示す式である。
【００８１】
【数８】
【００８２】
式（８）において、μ1（ｉ）（ｉ＝１〜４）は、表皮層における波長λｉ（ｉ＝１〜４）の光吸収係数を示し、μ２（ｉ）（ｉ＝１〜４）は、真皮層における波長λｉ（ｉ＝１〜４）の光吸収係数を示し、μ３（ｉ）（ｉ＝１〜４）は、皮下組織層における波長λｉ（ｉ＝１〜４）の光吸収係数を示す。ε１（ｉ）は、波長λｉ（ｉ＝１〜４）に対する水のモル吸光係数を示し、ε２（ｉ）は、波長λｉ（ｉ＝１〜４）に対するたんぱく質のモル吸光係数を示し、ε３（ｉ）は、波長λｉ（ｉ＝１〜４）に対する脂質のモル吸光係数を示し、ε４（ｉ）は、波長λｉ（ｉ＝１〜４）に対するグルコースのモル吸光係数を示す。ｇ１ｊ、ｇ２ｊ、ｇ３ｊ、ｇ４ｊ（ｊ＝１〜３）は、ｊ層におけるそれぞれ皮膚の主成分である水、たんぱく質、脂質、グルコースのモル濃度を示す。ただし、ｊ＝１は表皮層、ｊ＝２は真皮層、ｊ＝３は皮下組織層をそれぞれ示す。
【００８３】
式（９）は、真皮層（ｊ＝２）における水、たんぱく質、脂質、グルコースのモル濃度ｇ１２、ｇ２２、ｇ３２、ｇ４２を求める式である。式（９）のｇ４２を算出することで、真皮質に含まれるグルコースのモル濃度を求めることができる。
【００８４】
【数９】
【実施例２】
【００８５】
次に、本発明の第２の実施形態について詳しく説明する。
第２の実施形態は、第１の実施形態による血糖値測定装置１００と同じ構成であり、計測光強度取得部１０６、光路長取得部１０７、無吸収時光強度取得部１０８、光吸収係数算出部１０９の動作が異なる。
【００８６】
図８は、血糖値測定装置が血糖値を測定する動作を示す第２のフローチャートである。
まず、血糖値測定装置１００を動作させると、照射部１０４は、皮膚に対して波長λ１の短時間パルス光を照射する（ステップＳ１１）。ここで、波長λ１は、シミュレーション部１０１が伝搬光路長分布及び時間分解波形を算出した複数の波長の中の１つである。
【００８７】
照射部１０４が短時間パルス光を照射すると、受光部１０５は、照射部１０４から照射され、皮膚によって後方散乱した光を受光する（ステップＳ１２）。このとき、受光部１０５は、照射開始からの単位時間毎（例えば、１ピコ秒毎）の受光強度を内部メモリに登録しておく。
受光部１０５が受光を完了すると、計測光強度取得部１０６は、受光部１０５の内部メモリに格納されている受光強度から、ある時刻から時刻τまでの間の受光強度の時間分布を取得する（ステップＳ１３）。時刻τまでの光強度の変化としては、例えば図９に示す時間が選択される。すなわち、表皮を透過する光のスペクトルを精度良く測定したい場合には、１ｐｓ〜１０ｐｓの期間τ１を選択し、真皮を透過する光のスペクトルを精度良く測定したい場合には、２ｐｓ〜３０ｐｓの期間τ２を選択し、皮下組織を透過する光のスペクトルを精度良く検出したい場合には、５０ｐｓ〜７０ｐｓの期間τ３を選択する。τ１、τ２、τ３は、互いに時間軸上では重ならず、皮膚組織の表面に近い層を透過する光のスペクトルを検出しようとするほど、測定開始時間及び測定終了時間は時間軸上で早い時期に設定される。
【００８８】
計測光強度取得部１０６が、時刻τまでの間の受光強度の時間分布を取得すると、光路長取得部１０７は、光路長分布記憶部１０２が記憶する波長λ１の伝搬光路長分布から、ある時刻から時刻τまでの間の皮膚の各層の光路長Ｌ１〜Ｌ３を取得する（ステップＳ１４）。
また、計測光強度取得部１０６が、時刻τまでの間の受光強度を取得すると、無吸収時光強度取得部１０８は、時間分解波形記憶部１０３が記憶する波長λ１の時間分解波形から、ある時刻から時刻τまでの間の無吸収時光強度を取得する（ステップＳ１５）。
【００８９】
光路長取得部１０７が皮膚の各層の光路長を取得し、無吸収時光強度取得部１０８が無吸収時光強度を取得すると、光吸収係数算出部１０９は、式（１０）に基づいて、皮膚の各層の光吸収係数μ１〜μ３を算出する（ステップＳ１６）。ここで、光吸収係数μ１は、表皮層の光吸収係数を示し、光吸収係数μ２は、真皮層の光吸収係数を示し、光吸収係数μ３は、皮下組織層の光吸収係数を示す。
【００９０】
【数１０】
【００９１】
但し、ｌｎ（Ａ）は、Ａの自然対数を示す。また、Ｉ（ｔ）は、時刻ｔにおける受光部１０５の受光強度を示し、Ｉｉｎは、照射部１０４が照射した短時間パルス光の光強度を示す。また、Ｎ（ｔ）は、時間分解波形の時刻ｔにおける無吸収時光強度を示し、Ｎｉｎは、シミュレーション部１０１が照射のシミュレーションを行った光子の個数を示す。また、Ｌ１（ｔ）〜Ｌ３（ｔ）は、時刻ｔにおける皮膚の各層の光路長を示す。
【００９２】
光吸収係数算出部１０９は、皮膚の各層の光吸収係数μ１〜μ３を算出すると、光吸収係数算出部１０９が、皮膚の主成分の種類数と同じ数の波長に対して光吸収係数μ１〜μ３を算出したか否かを判定する（ステップＳ１７）。本実施形態では、皮膚の主成分を水、たんぱく質、脂質、グルコースの４種類として血糖値の測定を行うため、光吸収係数算出部１０９は、４種類の波長λ１〜λ４に対して光吸収係数μ１〜μ３を算出したか否かを判定する。ここで、波長λ１〜λ４は、シミュレーション部１０１が伝搬光路長分布及び時間分解波形を算出した複数の波長の中から選出する。
【００９３】
光吸収係数算出部１０９が、光吸収係数μ１〜μ３を算出していない波長λ１〜λ４があると判定した場合（ステップＳ１７：ＮＯ）、ステップＳ１に戻り、まだ光吸収係数μ１〜μ３を算出していない波長λ１〜λ４の光吸収係数μ１〜μ３の算出を行う。
他方、光吸収係数算出部１０９が、波長λ１〜λ４の光吸収係数μ１〜μ３を算出していると判定した場合（ステップＳ１７：ＹＥＳ）、濃度算出部１１０は、上述した式（７）に基づいて真皮質に含まれるグルコースの濃度を算出する（ステップＳ１８）。
【００９４】
このように、本実施形態によれば、吸収係数μ１〜μ３を、時刻τまでの間の光路長の積分値によって算出する。これにより、計測した受光強度Ｉ（ｔ）に含まれている誤差による吸収係数μ１〜μ３の算出結果に対する影響を少なくすることができる。
【００９５】
以上、図面を参照してこの発明の一実施形態について詳しく説明してきたが、具体的な構成は上述のものに限られることはなく、この発明の要旨を逸脱しない範囲内において様々な設計変更等をすることが可能である。
例えば、実施形態では、濃度定量方法を血糖値測定装置１００に実装し、皮膚の真皮層に含まれるグルコースの濃度を測定する場合を説明したが、これに限られず、濃度定量方法を、複数の光散乱媒質の層から形成される観測対象の任意の層における目的成分の濃度を定量する他の装置に用いても良い。
【００９６】
上述の血糖値測定装置１００は内部に、コンピュータシステムを有している。そして、上述した各処理部の動作は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータが読み出して実行することによって、上記処理が行われる。ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしても良い。
【００９７】
また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。
さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であっても良い。
【符号の説明】
【００９８】
１００…血糖値測定装置、１０１…シミュレーション部、１０２…光路長分布記憶部、１０３…時間分解波形記憶部、１０４…照射部（照射手段）、１０４ａ…照射用光ファイバー、１０５…受光部、１０５ａ…受光用光ファイバー、１０６…計測光強度取得部、１０７…光路長取得部、１０８…無吸収時光強度取得部、１０９…光吸収係数算出部（選択手段）、１１０…濃度算出部
【技術分野】
【０００１】
本発明のいくつかの様態は、複数の光散乱媒質の層から形成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置、濃度定量方法及びプログラムに関する。
【背景技術】
【０００２】
従来、血糖値の測定は、指先などから採血を行い、血中のグルコースに対する酵素活性を測ることで行っていた。しかし、このような血糖値の測定方法は、指先などから血液を採取して測定しなければならず、採血に手間と痛みを伴うことや、血液を付着させる測定チップが必要なことから、採血を必要としない非侵襲型の血糖値の測定方法が望まれている。
【０００３】
そこで、皮膚に近赤外光を照射し、その光吸収量からグルコースの濃度を求める方法が検討されている（例えば、特許文献１を参照）。具体的には、予めグルコース濃度と照射する光の波長と光の吸収量との関係を示す検量線を作成しておき、モノクロメーター等を用いてある波長域を走査し、その波長域の各波長に対する吸収量を求め、当該波長及び吸収量と検量線とを比較することでグルコース濃度を算出する。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第３９３１６３８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、従来の非侵襲血糖値測定方法は、光の入出射間距離を定めることによって、真皮層の近赤外吸収スペクトルを測定するため、測定した吸収スペクトルに、真皮層の吸収スペクトルのみならず表皮層や皮下組織層の吸収スペクトルも含まれ、観測される吸収係数の変化には表皮層や皮下組織層によるノイズが含まれてしまうという問題があった。
本発明は上記の点に鑑みてなされたものであり、その目的は、目的の層以外の層によるノイズの影響を軽減する濃度定量装置、濃度定量方法及びプログラムを提供することにある。
【課題を解決するための手段】
【０００６】
本発明のいくつかの態様は上記の課題を解決するためになされたものであり、複数の光散乱媒質の層から形成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置であって、前記観測対象に短時間パルス光を照射する照射手段と、前記短時間パルス光が前記観測対象によって後方散乱した光を受光する複数の受光部を有する受光手段と、前記複数の受光部のうち、前記短時間パルス光が前記任意の層によって後方散乱した光を受光する特定の受光部を選択する選択手段と、前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記特定の受光手段が受光した光の強度を取得する光強度取得手段と、前記観測対象に対して照射する短時間パルス光の、前記照射手段から前記特定の受光部に至る光の伝搬経路上に配置された複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、前記観測対象に対して照射し前記特定の受光部において受光する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段と、前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手段と、前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手段と、前記光強度取得手段が取得した光強度と前記光路長取得手段が取得した前記複数の光散乱媒質の層の各々の層の光路長と前記光強度モデル取得手段が取得した光強度モデルとに基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手段と、前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手段と、を備えることを特徴とする。
これにより、受光した光の時間分解波形から、任意の層の光吸収計数を選択的に算出することができる。そのため、算出した光吸収計数に基づいて目的成分の濃度を算出することにより、他の層によるノイズの影響を軽減し、精度の高い濃度の定量を行うことができる。
ここで、短時間パルス光とは、パルス幅が１０ｐｓｅｃ程度かそれ以下のパルス光を意味する。短時間パルス光として０．１ｐｓｅｃから１０ｐｓｅｃの範囲のパルス幅を持つパルス光を用いても良い。
【０００７】
また、本発明のいくつかの態様は、前記光強度取得手段は、前記観測対象の層の数ｎ以上となる複数の時刻ｔ１〜ｔｍにおける光強度を取得し、前記光吸収係数算出手段は、下記式（１）に基づいて任意の層の光吸収係数を算出する、ことを特徴とする。
【０００８】
【数１】
【０００９】
但し、ｌｎ（Ａ）は、Ａの自然対数を示し、Ｉ（ｔ）は、前記受光手段が時刻ｔにおいて受光した光強度を示し、Ｎ（ｔ）は、前記短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度を示し、Ｌｉ（ｔ）は、前記伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長を示し、μｉは、第ｉ層の光吸収係数を示す。
【００１０】
また、本発明のいくつかの態様は、前記照射手段は、複数の波長１〜ｑの光を照射し、前記光吸収係数算出手段は、前記任意の層における光吸収係数を前記照射手段が照射した複数の波長毎に算出し、前記濃度算出手段は、下記式（２）に基づいて前記任意の層における前記目的成分の濃度を算出する、ことを特徴とする。
【００１１】
【数２】
【００１２】
但し、μａ（ｉ）は、前記任意の層である第ａ層における波長ｉの光吸収係数を示し、ｇｊは、前記観測対象を形成する第ｊ成分のモル濃度を示し、εｊ（ｉ）は、第ｊ成分の波長ｉに対する光吸収係数を示し、ｐは、前記観測対象を形成する主成分の個数を示し、ｑは、照射手段が照射する波長の種類数を示す。
【００１３】
また、本発明のいくつかの態様は、前記照射手段から前記観測対象に短時間パルス光を照射する位置を照射位置、前記短時間パルス光が前記観測対象によって後方散乱した光が前記観測対象から前記受光部に向けて出射する位置を受光位置、前記照射位置と前記受光位置との間隔を照射受光間隔、としたときに、前記受光手段は、互いに照射受光間隔の異なる複数の受光部を有し、前記選択手段は、前記照射受光間隔に応じて定まる前記観測対象の内部への短時間パルス光の到達深さに基づいて、前記任意の層に短時間パルス光が伝搬可能な照射受光間隔を有する受光部を選択することを特徴とする。
【００１４】
また、本発明のいくつかの態様は、前記照射手段は、短時間パルス光を前記観測対象の表面に伝送する照射用光ファイバーを有し、前記受光部は、前記観測対象によって後方散乱した光を伝送する受光用光ファイバーを有し、前記照射用光ファイバーと前記受光用光ファイバーは、両者の光ファイバーコアの中心間隔を所定の照射受光間隔だけ離間して固定するプローブ装置に装着されており、前記プローブ装置の先端部に露出した前記照射用光ファイバーの先端部と前記受光用光ファイバーの先端部とを前記観測対象の表面に接触させることによって、前記観測対象への短時間パルスの照射処理と、前記短時間パルス光が前記観測対象によって後方散乱した光の受光処理とを行うことを特徴とする。
【００１５】
また、本発明のいくつかの態様は、前記受光部は、互いに等しい照射受光間隔で配置された複数の受光用光ファイバーを有することを特徴とする。
【００１６】
また、本発明のいくつかの態様は、前記受光部は、互いに等しい照射受光間隔で配置された複数の受光用光ファイバーによって伝送された光を同一受光面上に集光する集光素子を有することを特徴とする。
【００１７】
また、本発明のいくつかの態様は、前記照射用光ファイバーと前記受光用光ファイバーは、前記照射手段が照射する前記複数の波長１〜ｑの光の波長分散を補償する分散補償型シングルモード光ファイバーであることを特徴とする。
【００１８】
また、本発明のいくつかの態様は、前記照射用光ファイバーと前記受光用光ファイバーは、前記照射手段が照射する前記複数の波長１〜ｑの光の波長分散に伴う群遅延時間差が、前記複数の光散乱媒質層のうち最も表面側の層の伝搬光路長分布のピークに対応する伝搬時間よりも短いことを特徴とする。
【００１９】
また、本発明のいくつかの態様は、複数の光散乱媒質の層から形成される観測対象に対して短時間パルス光を照射する照射手段と、前記短時間パルス光が前記観測対象によって後方散乱した光を受光する複数の受光部を有する受光手段と、前記複数の受光部のうち、前記短時間パルス光が前記任意の層によって後方散乱した光を受光する特定の受光部を選択する選択手段と、前記観測対象に対して照射する短時間パルス光の、前記照射手段から前記特定の受光部に至る光の伝搬経路上に配置された複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、前記観測対象に対して照射し前記特定の受光部において受光する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段とを備え、前記観測対象のうち任意の層における目的成分の濃度を定量する濃度定量装置を用いた濃度定量方法であって、光強度取得手段は、前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記特定の受光部が受光した光の強度を取得し、光路長取得手段は、前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得し、光強度モデル取得手段は、前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得し、光吸収係数算出手段は、前記光強度取得手段が取得した光強度と前記光路長取得手段が取得した前記複数の光散乱媒質の層の各々の層の光路長と前記光強度モデル取得手段が取得した光強度モデルとに基づいて、前記任意の層の光吸収係数を算出し、濃度算出手段は、前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する、ことを特徴とする。
【００２０】
また、本発明のいくつかの態様は、複数の光散乱媒質の層から形成される観測対象に対して短時間パルス光を照射する照射手段と、前記短時間パルス光が前記観測対象によって後方散乱した光を受光する複数の受光部を有する受光手段と、前記複数の受光部のうち、前記短時間パルス光が前記任意の層によって後方散乱した光を受光する特定の受光部を選択する選択手段と、前記観測対象に対して照射する短時間パルス光の、前記照射手段から前記特定の受光部に至る光の伝搬経路上に配置された複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、前記観測対象に対して照射し前記特定の受光部において受光する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段とを備え、前記観測対象のうち任意の層における目的成分の濃度を定量する濃度定量装置を、前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記特定の受光部が受光した光の強度を取得する光強度取得手段、前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手段、前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手段、前記光強度取得手段が取得した光強度と前記光路長取得手段が取得した前記複数の光散乱媒質の層の各々の層の光路長と前記光強度モデル取得手段が取得した光強度モデルとに基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手段、前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手段、として動作させるためのプログラムである。
【図面の簡単な説明】
【００２１】
【図１】本発明による血糖値測定装置の構成を示す概略ブロック図である。
【図２】シミュレーション部が算出した各層の伝搬光路長分布を示すグラフである。
【図３】シミュレーション部が算出した時間分解波形を示すグラフである。
【図４】皮膚の主成分の吸収スペクトルを示すグラフである。
【図５】皮膚組織の断面図である。
【図６】プローブ装置の断面図である。
【図７】血糖値測定装置が血糖値を測定する動作を示す第１のフローチャートである。
【図８】血糖値測定装置が血糖値を測定する動作を示す第２のフローチャートである。
【図９】光路長取得部が取得する所定時間の光路長を示す図である。
【発明を実施するための形態】
【実施例１】
【００２２】
以下、図面を参照しながら本発明の第１の実施形態について詳しく説明する。図１は、本発明の第１の実施形態による血糖値測定装置の構成を示す概略ブロック図である。
血糖値測定装置１００（濃度定量装置）は、シミュレーション部１０１、光路長分布記憶部１０２（光路長分布記憶手段）、時間分解波形記憶部１０３（時間分解波形記憶手段）、照射部１０４（照射手段）、受光部１０５（受光手段）、計測光強度取得部１０６（光強度取得手段）、光路長取得部１０７（光路長取得手段）、無吸収時光強度取得部１０８（光強度モデル取得手段）、光吸収係数算出部１０９（光吸収係数算出手段）、濃度算出部１１０（濃度算出手段）を備える。
血糖値測定装置１００は、皮膚（観測対象）の真皮層（任意の層）に含まれるグルコース（目的成分）の濃度を測定する。
【００２３】
シミュレーション部１０１は、光吸収係数がゼロの皮膚モデルに対して光を照射するシミュレーションを行う。
光路長分布記憶部１０２は、光吸収係数がゼロの皮膚モデルの伝搬光路長分布を記憶する。
時間分解波形記憶部１０３は、光吸収係数がゼロの皮膚モデルの時間分解波形を記憶する。
照射部１０４は、皮膚に対して短時間パルス光を照射する。
受光部１０５は、短時間パルス光が皮膚によって後方散乱した光を受光する。
計測光強度取得部１０６は、受光部１０５が受光した光のある時刻における光強度を取得する。
光路長取得部１０７は、光路長分布記憶部１０２からある時刻における光路長を取得する。
無吸収時光強度取得部１０８は、時間分解波形記憶部１０３からある時刻における光強度を取得する。
光吸収係数算出部１０９は、短時間パルス光を照射した皮膚の真皮層における光吸収係数を算出する。
濃度算出部１１０は、真皮層に含まれるグルコースの濃度を算出する。
ここで、短時間パルス光とは、パルス幅が１０ｐｓｅｃ程度かそれ以下のパルス光を意味する。短時間パルス光として０．１ｐｓｅｃから１０ｐｓｅｃの範囲のパルス幅を持つパルス光を用いても良い。
【００２４】
そして、血糖値測定装置１００において、照射部１０４は、皮膚に短時間パルス光を照射し、受光部１０５は、短時間パルス光が皮膚によって後方散乱した光を受光し、計測光強度取得部１０６は、時刻ｔにおいて受光部１０５が受光した光の強度を取得する。次に、光路長取得部１０７は、光路長分布記憶部１０２から、皮膚モデルにおける伝搬光路長分布の時刻ｔにおける皮膚の各層の光路長を取得し、無吸収時光強度取得部１０８は、時間分解波形記憶部１０３から、皮膚モデルにおける短時間パルス光の時間分解波形の時刻ｔにおける光の強度を取得する。
【００２５】
次に、光吸収係数算出部１０９は、計測光強度取得部１０６が取得した光強度と光路長取得部１０７が取得した皮膚の各層の光路長と無吸収時光強度取得部１０８が取得した光強度とに基づいて、皮膚の真皮層の光吸収係数を算出し、濃度算出部１１０は、光吸収係数算出部１０９が算出した光吸収係数に基づいて、真皮層におけるグルコースの濃度を算出する。
これにより、真皮層以外の層によるノイズの影響を軽減して、真皮層に含まれるグルコースの濃度を算出することができる。
【００２６】
次に、血糖値測定装置１００の動作を説明する。
血糖値測定装置１００は、血糖値を測定する前に、予め皮膚モデルの各層における伝搬光路長分布と時間分解波形とを算出しておく必要がある。
まず、皮膚モデルの伝搬光路長分布及び時間分解波形の算出方法を説明する。
初めに、シミュレーション部１０１は、皮膚モデルを生成する。皮膚モデルの生成は、皮膚の各層の光散乱係数、光吸収係数及び厚みを決定することで行う。ここで、皮膚の各層の散乱係数及び厚みは、個体による差が少ないため、予めサンプルを取ることなどによって決定すると良い。なお、表皮層の厚みは略０．３ｍｍ、真皮層の厚みは略１．２ｍｍ、皮下組織層の厚みは略３．０ｍｍである。
また、ここで用いる皮膚モデルの光吸収係数はゼロとする。これは、当該皮膚モデルを用いて光吸収量を算出するためである。
【００２７】
シミュレーション部１０１は、皮膚モデルを生成すると、当該皮膚モデルに光を照射するシミュレーションを行う。このとき、照射部１０４の位置と受光部１０５の位置との間の距離を決定しておく必要がある。シミュレーションは、モンテカルロ法を用いて行うと良い。モンテカルロ法によるシミュレーションは、例えば以下のように行われる。
まず、シミュレーション部１０１は、照射する光のモデルを光子（光束）とし、当該光子を皮膚モデルに照射する計算を行う。皮膚モデルに照射された光子は、皮膚モデル内を移動する。このとき、光子は、次に進む点までの距離Ｌ及び方向θを乱数Ｒによって決定する。シミュレーション部１０１は、光子が次に進む点までの距離Ｌの計算を、式（３）によって行う。
【００２８】
【数３】
【００２９】
但し、ｌｎ（Ａ）は、Ａの自然対数を示す。また、μｓは、皮膚モデルの第ｓ層（表皮層、真皮層、皮下組織層の何れか）の散乱係数を示す。
また、シミュレーション部１０１は、光子が次に進む点までの方向θの計算を、式（４）によって行う。
【００３０】
【数４】
【００３１】
但し、ｇは、散乱角度の余弦（ｃｏｓ）の平均である非等方性パラメータを示し、皮膚の非等方性パラメータは、略０．９である。
シミュレーション部１０１は、上記式（３）、式（４）の計算を単位時間毎に繰り返すことにより、照射部１０４から受光部１０５までの光子の移動経路を算出することができる。シミュレーション部１０１は、複数の光子について移動距離の算出を行う。例えば、シミュレーション部１０１は、１０８個の光子について移動距離を算出する。
【００３２】
図２は、シミュレーション部が算出した各層の伝搬光路長分布を示すグラフである。
図２の横軸は光子の照射からの経過時間を示し、縦軸は光路長の対数表示を示している。シミュレーション部１０１は、受光部１０５に到達した光子の各々の移動経路を、移動経路が通過する層毎に分類する。そして、シミュレーション部１０１は、単位時間毎に到達した光子の移動経路の平均長を分類された層毎に算出することで、図２に示すような皮膚の各層の伝搬光路長分布を算出する。
【００３３】
図３は、シミュレーション部が算出した時間分解波形を示すグラフである。
図３の横軸は光子の照射からの経過時間を示し、縦軸は受光部１０５が検出した光子数を示している。シミュレーション部１０１は、単位時間毎に受光部１０５に到達した光子の個数を算出することで、図３に示すような皮膚モデルの時間分解波形を算出する。
上述したような処理により、シミュレーション部１０１は、複数の波長に対して、皮膚モデルの伝搬光路長分布及び時間分解波形を算出する。このとき、シミュレーション部１０１は、皮膚の主成分（水、たんぱく質、脂質、グルコース等）の吸収スペクトルの差が大きくなる波長について伝搬光路長分布及び時間分解波形を算出すると良い。
【００３４】
図４は、皮膚の主成分の吸収スペクトルを示すグラフである。
図４の横軸は照射する光の波長を示し、縦軸は吸収係数を示している。図４を参照すると、グルコースの吸収係数は、波長が１６００ｎｍのときに極大となり、水の吸収係数は、波長が１４５０ｎｍのときに極大となる。そのため、シミュレーション部１０１は、例えば１４００ｎｍ、１４５０ｎｍ、１５００ｎｍ、１６００ｎｍ、１６８０ｎｍ、１７２０ｎｍ、１７４０ｎｍといった皮膚の主成分の吸収スペクトルの差が大きくなる波長について伝搬光路長分布及び時間分解波形を算出すると良い。
【００３５】
シミュレーション部１０１は、複数の波長に対する皮膚モデルの伝搬光路長分布及び時間分解波形を算出すると、伝搬光路長分布の情報を光路長分布記憶部１０２に記憶させ、時間分解波形の情報を時間分解波形記憶部１０３に記憶させる。
【００３６】
ここで、図５に基づいて、観測対象である皮膚組織の構造を説明する。
図５に示すように、人間の皮膚組織は、表皮と、真皮と、皮下組織の３層によって形成されている。表皮は、最も外側にある厚さ０．２〜０．３mmの薄い層で、角質層、顆粒層、有棘層、底層等を含む。真皮は、表皮と皮下組織の間に存在する厚さ０．５〜２ｍｍの層で、真皮内には神経、毛根、皮脂腺や汗腺、毛包、血管、リンパ管が存在する。皮下組織は、真皮の下にある厚さ１〜３ｍｍの層で、大部分が皮下脂肪でできている。
【００３７】
真皮内には毛細血管等が発達しており、血中グルコースに応じた物質移動が速やかに起こり、血中グルコース濃度（血糖値）に対して真皮中のグルコース濃度は追随して変化すると考えられている。そのため、血糖値測定装置１００では、照射部１０４から皮膚組織の表面に光を照射し、皮膚組織内を透過、拡散して後方散乱された光を受光部１０５で検出することによって、真皮中のグルコース濃度を測定している。
【００３８】
照射部１０４は、光源１０４ｂと、光源１０４ｂから射出された光を表皮に向けて伝送する照射用光ファイバー１０４ａとを含む。受光部１０５は、光電変換素子１０５ｂと、皮膚組織内を透過、拡散して表皮から出射した光を光電変換素子１０５ｂに向けて伝送する受光用光ファイバー１０５ａとを含む。照射用光ファイバー１０４ａと受光用光ファイバー１０５ａは、光ファイバーコアの中心間隔が所定の照射受光間隔だけ離間した状態で、互いに光軸（コアの中心軸）が平行となるようにして図示略のプローブ装置に固定されている。
【００３９】
受光用光ファイバー１０５ａに到達する光の経路分布ＲはＵ字型の形状をしている。実際には、図示点線で示したような複数のＵ字型の経路が存在し、それぞれの経路における光の伝搬距離及び伝搬時間は互いに異なっている。皮膚組織におけるこのような経路分布の到達深さは、照射部１０４と受光部１０５との間隔、具体的には、照射用光ファイバー１０４ａと受光用光ファイバー１０５ａのコア中心間の距離である照射受光間隔Ｗによって変化する。したがって、照射部１０４と受光部１０５との照射受光間隔Ｗを適切に設定することにより、主として真皮を透過する光の吸収スペクトルを測定することができる。
【００４０】
図６（ａ）〜図６（ｃ）は、照射用光ファイバー１０４ａと受光用光ファイバー１０５ａの光軸に垂直な面で切った断面図である。なお、図６（ａ）〜図６（ｃ）において、符号Ａで示した光ファイバーは照射用光ファイバー１０４ａであり、符号Ｂ１〜Ｂ４で示した光ファイバーは受光用光ファイバー１０５ａである。
【００４１】
図６（ａ）のプローブ装置１３０は、左端に配置された１つの照射用光ファイバー１０４ａと、その右側に等間隔で配置された複数の受光用光ファイバー１０５ａとを含む。プローブ装置１３０では、照射用光ファイバー１０４ａ及び光源１０４ｂ（図５参照）が本発明の照射手段に該当し、その周囲に配置された複数の受光用光ファイバー１０５ａ及び光電変換素子１０５ｂ（図５参照）が本発明の受光手段に該当する。
【００４２】
照射用光ファイバー１０４ａと受光用光ファイバー１０５ａは、光軸に垂直な方向に一列に配置されている。プローブ装置１３０では、プローブ装置１３０の先端部に露出した照射用光ファイバー１０４ａの先端部と受光用光ファイバー１０５ａの先端部とを皮膚表面に接触させることによって、皮膚組織への短時間パルス光の照射処理と、短時間パルス光が皮膚組織によって後方散乱した光の受光処理とを行うようになっている。
【００４３】
複数の受光用光ファイバー１０５ａには、照射用光ファイバー１０４ａに対して第１の照射受光間隔で配置された第１の受光用光ファイバー１０５ａ（符号Ｂ１で示す）と、その右側に第１の照射受光間隔よりも大きい第２の照射受光間隔で配置された第２の受光用光ファイバー１０５ａ（符号Ｂ２で示す）と、その右側に第２の照射受光間隔よりも大きい第３の照射受光間隔で配置された第３の受光用光ファイバー１０５ａ（符号Ｂ３で示す）と、その右側に第３の照射受光間隔よりも大きい第４の照射受光間隔で配置された第４の受光用光ファイバー１０５ａ（符号Ｂ４で示す）と、が含まれる。
【００４４】
複数の受光用光ファイバー１０５ａは、皮膚組織への到達深さに応じて、１つ又は複数の受光用光ファイバー１０５ａが選択され使用される。例えば、表皮を透過する光のスペクトルを精度良く測定したい場合には、照射受光間隔が最も小さい第１の受光用光ファイバー１０５ａ（符号Ｂ１で示す）を選択し、真皮を透過する光のスペクトルを精度良く測定したい場合には、第２の照射受光間隔で配置される第２の受光用光ファイバー１０５ａ（符号Ｂ２で示す）と第３の照射受光間隔で配置される第３の受光用光ファイバー１０５ａ（符号Ｂ３で示す）とを選択するという方法が用いられる。また、皮膚組織の厚みは個人差があり、例えば大人と子供、男性と女性では、皮膚組織の厚みが異なる場合がある。そのため、同じく真皮の情報を得る場合であっても、表皮及び真皮の厚みが厚いと予想される被検体には、厚みが薄いと予想される被検体よりも照射受光間隔の大きい受光用光ファイバー１０５ａを用いることも考えられる。
【００４５】
図６（ａ）では、４本の受光用光ファイバー１０５ａのみを示したが、受光用光ファイバー１０５ａの数はこれに限らず、３本若しくは５本以上とすることも可能である。受光用光ファイバー１０５ａの数を増やせば、皮膚組織への到達深さを多段階で制御することができる。
【００４６】
図６（ｂ）のプローブ装置１３１は、中央に配置された１つの照射用光ファイバー１０４ａと、その両側に等間隔で配置された複数の受光用光ファイバー１０５ａとを含む。プローブ装置１３１では、照射用光ファイバー１０４ａ及び光源１０４ｂ（図５参照）が本発明の照射手段に該当し、その周囲に配置された複数の受光用光ファイバー１０５ａ及び光電変換素子１０５ｂ（図５参照）が本発明の受光手段に該当する。
【００４７】
照射用光ファイバー１０４ａと受光用光ファイバー１０５ａは、光軸に垂直な方向に一列に配置されている。プローブ装置１３１では、プローブ装置１３１の先端部に露出した照射用光ファイバー１０４ａの先端部と受光用光ファイバー１０５ａの先端部とを皮膚表面に接触させることによって、皮膚組織への短時間パルス光の照射処理と、短時間パルス光が皮膚組織によって後方散乱した光の受光処理とを行うようになっている。
【００４８】
複数の受光用光ファイバー１０５ａには、照射用光ファイバー１０４ａの両側に第１の照射受光間隔で配置された２つの第１の受光用光ファイバー１０５ａ（符号Ｂ１で示す）と、その両側に第１の照射受光間隔よりも大きい第２の照射受光間隔で配置された２つの第２の受光用光ファイバー１０５ａ（符号Ｂ２で示す）と、が含まれる。
【００４９】
複数の受光用光ファイバー１０５ａは、皮膚組織への到達深さに応じて、１つ又は複数の受光用光ファイバー１０５ａが選択され使用される。例えば、表皮を透過する光のスペクトルを精度良く測定したい場合には、照射受光間隔が最も小さい第１の受光用光ファイバー１０５ａ（符号Ｂ１で示す）を選択し、真皮を透過する光のスペクトルを精度良く測定したい場合には、第１の照射受光間隔で配置される第１の受光用光ファイバー１０５ａ（符号Ｂ１で示す）と第２の照射受光間隔で配置される第２の受光用光ファイバー１０５ａ（符号Ｂ２で示す）とを選択するという方法が用いられる。図６（ｂ）のプローブ装置１３１は、図６（ａ）のプローブ装置１３０に比べて、選択できる照射受光間隔が少ないが、１つの照射受光間隔に対して複数の受光用光ファイバー１０５ａが設けられているので、受光強度は大きくなる。そのため、ノイズの少ない測定が可能となる。
【００５０】
図６（ｂ）では、照射受光間隔の異なる２種類の受光用光ファイバー１０５ａを示したが、照射受光間隔は２種類に限らず、３種類以上とすることもできる。例えば、３種類の照射受光間隔で測定を行う場合には、第２の受光用光ファイバー１０５ａの両側に第２の照射受光間隔よりも大きい第３の照射受光間隔で配置された２つの第３の受光用光ファイバー１０５ａを設ければよい。照射受光間隔を４種類以上とする場合も同様である。
【００５１】
図６（ｃ）のプローブ装置１３２は、中央に配置された１つの照射用光ファイバー１０４ａと、その周囲に配置された複数の受光用光ファイバー１０５ａとを含む。プローブ装置１３２では、照射用光ファイバー１０４ａ及び光源１０４ｂ（図５参照）が本発明の照射手段に該当し、その周囲に配置された複数の受光用光ファイバー１０５ａ及び光電変換素子１０５ｂ（図５参照）が本発明の受光手段に該当する。
【００５２】
照射用光ファイバー１０４ａと受光用光ファイバー１０５ａは、光軸に垂直な平面内で最密充填されるように配置されている。すなわち、受光用光ファイバー１０５ａは、照射用光ファイバー１０４ａを中心に６０°の内角を持つように放射方向に並んで配置されている。プローブ装置１３２では、プローブ装置１３２の先端部に露出した照射用光ファイバー１０４ａの先端部と受光用光ファイバー１０５ａの先端部とを皮膚表面に接触させることによって、皮膚組織への短時間パルスの照射処理と、短時間パルス光が皮膚組織によって後方散乱した光の受光処理とを行うようになっている。
【００５３】
複数の受光用光ファイバー１０５ａには、照射用光ファイバー１０４ａの周囲に第１の照射受光間隔で配置された６つの第１の受光用光ファイバー１０５ａ（符号Ｂ１で示す）と、その周囲に第１の照射受光間隔よりも大きい第２の照射受光間隔で配置された６つの第２の受光用光ファイバー１０５ａ（符号Ｂ２で示す）及び第２の照射受光間隔よりも大きい第３の照射受光間隔で配置された６つの第３の受光用光ファイバー１０５ａ（符号Ｂ３で示す）と、が含まれる。
【００５４】
複数の受光用光ファイバー１０５ａは、皮膚組織への到達深さに応じて、１つ又は複数の受光用光ファイバー１０５ａが選択され使用される。図６（ｃ）のプローブ装置１３２は、図６（ｂ）のプローブ装置１３１に比べて、１つの照射受光間隔に対して設けられる受光用光ファイバー１０５ａの数が多いため、よりノイズの少ない測定が可能となる。
【００５５】
なお、図６（ｂ）及び図６（ｃ）では、１つの照射受光間隔に対して複数の受光用光ファイバー１０５ａが設けられている。この場合、複数の受光用光ファイバー１０５ａで伝送される光は、受光用光ファイバー毎に設けられた光電変換素子によって電気信号に変換しても良いし、同一照射受光間隔で配置された複数の受光用光ファイバー１０５ａで伝送した光を集光素子によって同一光電変換素子の同一受光面上に集光し、これらを一括して電気信号に変換してもよい。後者の方法を用いれば、光電変換素子を複数の受光用光ファイバー１０５ａに対して共通化することができるので、部材コストを低減することができる。
【００５６】
受光用光ファイバー毎に１つの光電変換素子を設ける場合には、１つ受光用光ファイバーと１つの光電変換素子によって１つの受光部が構成される。複数の受光用光ファイバーに対して１つの光電変換素子が設けられる場合には、複数の受光用光ファイバーと１つの光電変換素子によって１つの受光部が構成される。
【００５７】
照射用光ファイバー１０４ａと受光用光ファイバー１０５ａは、照射部１０４が照射する複数の波長１〜ｑの光の波長分散を補償する分散補償型シングルモード光ファイバーであることが望ましい。これにより、波長分散による測定誤差を低減することができる。
【００５８】
照射用光ファイバー１０４ａと受光用光ファイバー１０５ａは、分散補償型でないマルチモード光ファイバーを用いることもできる。マルチモード光ファイバーを用いることで、シングルモード光ファイバーを用いる場合に比べて受光面積（コア径）を大きくすることができ、ノイズの少ない測定を行うことが可能となる。ただし、波長分散による測定誤差が大きくなるため、照射用光ファイバー１０４ａと受光用光ファイバー１０５ａにおいては、照射部１０４が照射する複数の波長１〜ｑの光の波長分散に伴う群遅延時間差が、皮膚組織に含まれる複数の光散乱媒質層のうち最も表面側の層である表皮の伝搬光路長分布のピークに対応する伝搬時間よりも短くなるように、各々の光ファイバーの長さが設定されることが望ましい。
【００５９】
例えば、ステップインデックス型マルチモード光ファイバーを用いる場合を考える。ステップインデックス型マルチモード光ファイバの群遅延時間差は下記の式（５）で表せる。
【００６０】
【数５】
【００６１】
但し、Ｔｏは光ファイバの光軸上を進む最も早い光の群遅延時間を示し、Δは光ファイバーのコアとクラッドとの比屈折率差（＝１％）を示し、θcは臨界角、ｎ１はコアの屈折率を示し、Lは光ファイバーの長さを示し、ｃは光速を示す。Ｔｏは図５で示した真皮を通る複数の経路のうち最も短い伝搬距離を有する経路、すなわち、短時間パルス光の到達深さが表皮と真皮の界面付近となる光の経路である。ΔＴmが真皮伝搬光の最短到達時間より短くなる様なマルチモードファイバにすることでシングルモードファイバと比べて信号光の光強度を大きくできる。
【００６２】
例えば、群遅延時間差を１．０ｐｓ以内で、光の強度を受光するための受光用光ファイバー１０５ａの長さは、コアの屈折率を１．５、光速を３×１０８ｍ／ｓとして、式（５）より、２０ｍｍと算出される。
【００６３】
マルチモードファイバ内での光の伝搬を幾何学的に説明すると、光軸に沿って（θ１＝０）進む光がｈだけ進む間に、光軸に対して角度θ１で進む光はｄ＝ｈ／（ｃｏｓθ１）だけ進む。スネルの法則からｃｏｓθ１＝ｎ２／ｎ１（但し、ｎ２はクラッドの屈折率を示し、ｎ１はコアの屈折率を示す）なので、ｄ＝ｈ（ｎ１／ｎ２）になる。
【００６４】
ここで、比屈折率差Δ（＝１−ｎ２／ｎ１）は１％なので、ｎ１／ｎ２＝１．０１である。よってｄ＝１．０１×ｈになる。ｈ＝２０ｍｍだからｄ＝２０．２ｍｍである。つまり、hが２０ｍｍのときのｈとｄの伝搬経路差はｄ−ｈ＝０．２ｍｍである。そして、屈折率ｎの物質中を進む光の速度ｖ＝空気中の光の速度／ｎなので、この場合ｖ＝３×１０８／１．５＝２×１０８ｍ／ｓとなる。よって、伝搬時間差ｔ＝伝搬経路差／光の速度という関係より伝搬経路差＝伝搬時間差ｔ×光の速度＝２０ｍｍになる。
【００６５】
次に、血糖値測定装置１００が血糖値を測定する動作について説明する。
図７は、血糖値測定装置が血糖値を測定する動作を示す第１のフローチャートである。
まず、ユーザが血糖値測定装置１００を皮膚にあてがい、測定開始スイッチ（図示せず）の押下等によって血糖値測定装置１００を動作させると、照射部１０４は、皮膚に対して波長λ１の短時間パルス光を照射する（ステップＳ１）。ここで、波長λ１は、シミュレーション部１０１が伝搬光路長分布及び時間分解波形を算出した複数の波長の中の１つである。
【００６６】
照射部１０４が短時間パルス光を照射すると、受光部１０５は、照射部１０４から照射され、皮膚によって後方散乱した光を受光する（ステップＳ２）。このとき、受光部１０５は、照射開始からの単位時間毎（例えば、１ピコ秒毎）の受光強度を内部メモリに登録しておく。
【００６７】
ここで、対象となる受光部１０５は、受光手段に設けられた複数の受光部の中から選択された特定の受光部である。受光手段に設けられる複数の受光部１０５は、それぞれ受光量に応じた電気信号を受光強度として図示略の内部メモリに記憶している。しかし、それらの受光強度のデータが全て光吸収係数算出部１０９によって取得されるわけではなく、特定の受光部１０５の内部メモリに記憶された受光強度のデータが光吸収係数算出部１０９によって選択的に取得される。すなわち、光吸収係数算出部１０９が本発明の選択手段に該当する。光吸収係数算出部１０９は、ユーザーがメニュー画面に従って直接的又は間接的に選択した受光部を選択しても良いし、予めデフォルトで設定された受光部を選択しても良い。前者の例としては、直接ユーザーが「受光部１」「受光部２」などの表示ボタンを操作することにより選択するものや、ユーザーが被検体の年齢、性別などをメニュー画面で選択し、年齢、性別などに関連付けられた受光部を間接的に選択するものなどが挙げられる。
【００６８】
受光部１０５が受光を完了すると、計測光強度取得部１０６は、光吸収係数算出部１０９によって選択される特定の受光部１０５の内部メモリに格納されている、異なる時刻ｔにおける受光強度Ｉ（ｔ）を皮膚の層の数と同じ数だけ取得する（ステップＳ３）。すなわち、計測光強度取得部１０６は、３つの異なる時刻ｔ１〜ｔ３における受光強度Ｉ（ｔ１）〜Ｉ（ｔ３）を取得する。ここで、皮膚の層の数と同じ数だけ受光強度を取得する理由は、後述する処理において、皮膚の各層の吸収係数を連立方程式によって算出するためである。
【００６９】
また、計測光強度取得部が光強度を取得する時刻ｔ１〜ｔ３は、皮膚の各層の伝搬光路長分布のピークとなる時刻であると良い。すなわち、照射部１０４が短時間パルス光を照射した時刻に、図２に示すグラフにおいて皮膚の各層の光路長が極大となる時間を加算した時刻の光強度をそれぞれ取得すると良い。
【００７０】
計測光強度取得部１０６が、受光強度Ｉ（ｔ１）〜Ｉ（ｔ３）を取得すると、光路長取得部１０７は、光路長分布記憶部１０２が記憶する波長λ１の伝搬光路長分布から、時刻ｔ１〜ｔ３における皮膚の各層の光路長Ｌ１（ｔ１）〜Ｌ１（ｔ３）、Ｌ２（ｔ１）〜Ｌ２（ｔ３）、Ｌ３（ｔ１）〜Ｌ３（ｔ３）を取得する（ステップＳ４）。
また、計測光強度取得部１０６が、受光強度Ｉ（ｔ１）〜Ｉ（ｔ３）を取得すると、無吸収時光強度取得部１０８は、時間分解波形記憶部１０３が記憶する波長λ１の時間分解波形から、時刻ｔ１〜ｔ３における無吸収時光強度Ｎ（ｔ１）〜Ｎ（ｔ３）を取得する（ステップＳ５）。
【００７１】
光路長取得部１０７が皮膚の各層の光路長を取得し、無吸収時光強度取得部１０８が無吸収時光強度を取得すると、光吸収係数算出部１０９は、式（６）に基づいて、皮膚の各層の光吸収係数μ１〜μ３を算出する（ステップＳ６）。ここで、光吸収係数μ１は、表皮層の光吸収係数を示し、光吸収係数μ２は、真皮層の光吸収係数を示し、光吸収係数μ３は、皮下組織層の光吸収係数を示す。
【００７２】
【数６】
【００７３】
但し、ｌｎ（Ａ）は、Ａの自然対数を示す。また、Ｉｉｎは、照射部１０４が照射した短時間パルス光の光強度を示す。また、Ｎｉｎは、シミュレーション部１０１が照射のシミュレーションを行った光子の個数を示す。
光吸収係数算出部１０９が皮膚の各層の光吸収係数μ１〜μ３を算出すると、光吸収係数算出部１０９は、求める成分の種類数αとして、xＣ2≧αを満足する波長の種類数（例えばα=3→x=3、α=4→x=4、α=5→x=4、）に対して光吸収係数μ１〜μ３を算出したか否かを判定する（ステップＳ７）。本実施形態では、皮膚の主成分を水、たんぱく質、脂質、グルコースの４種類として血糖値の測定を行うため、光吸収係数算出部１０９は、４種類の波長λ１〜λ４に対して光吸収係数μ１〜μ３を算出したか否かを判定する。ここで、波長λ１〜λ４は、シミュレーション部１０１が伝搬光路長分布及び時間分解波形を算出した複数の波長の中から選出する。
【００７４】
光吸収係数算出部１０９が、光吸収係数μ１〜μ３を算出していない波長λ１〜λ４があると判定した場合（ステップＳ７：ＮＯ）、ステップＳ１に戻り、まだ光吸収係数μ１〜μ３を算出していない波長λ１〜λ４の光吸収係数μ１〜μ３の算出を行う。
他方、光吸収係数算出部１０９が、波長λ１〜λ４の光吸収係数μ１〜μ３を算出していると判定した場合（ステップＳ７：ＹＥＳ）、濃度算出部１１０は、式（７）に基づいて真皮質に含まれるグルコースの濃度を算出する（ステップＳ８）。
【００７５】
【数７】
【００７６】
但し、μ２（１）〜μ２（４）は、真皮層における波長λ１〜λ４の光吸収係数を示す。また、ｇ１〜ｇ４は、真皮層におけるそれぞれ皮膚の主成分である水、たんぱく質、脂質、グルコースのモル濃度を示す。また、ε１（１）〜ε１（４）は、波長λ１〜λ４に対する水のモル吸光係数を示し、ε２（１）〜ε２（４）は、波長λ１〜λ４に対するたんぱく質のモル吸光係数を示し、ε３（１）〜ε３（４）は、波長λ１〜λ４に対する脂質のモル吸光係数を示し、ε４（１）〜ε４（４）は、波長λ１〜λ４に対するグルコースのモル吸光係数を示す。
つまり、式（７）のｇ４を算出することで、真皮質に含まれるグルコースのモル濃度を求めることができる。
【００７７】
ここで、式（７）によってグルコースのモル濃度を求めることができる理由を説明する。皮膚の散乱係数の波長依存性は小さいため、無吸収時光強度Ｎ（ｔ）及び光路長Ｌｎ（ｔ）の波長に対する変化は無視することができる。また、ベア・ランベルト（Beer-Lambert）の法則により、吸光度＝モル吸光係数×モル濃度で表すことができる。これにより、２波長で得られた時間分解計測より、無吸収時光強度Ｎ（ｔ）を消去することで、真皮層において得られる吸収計数差と皮膚を形成する各成分のモル吸光係数との関係式を示す式（７）を導くことができる。
【００７８】
このように、本実施形態によれば、短時間パルス光を照射し、所定の時刻において受光した光の強度に基づいてグルコースの濃度を定量する。これにより、所定の時刻において受光した光の中から、真皮層の吸収計数を選択的に算出することができる。これにより、特定の皮膚の層のグルコースの濃度を算出することができ、他の層によるノイズの影響を軽減し、精度の高い血糖値を算出することが可能になる。
【００７９】
なお、式（７）では、複数波長での光吸収係数の差分に基づいてグルコースのモル濃度を求めた。しかし、複数波長での光吸収係数の差分によらずに下式（８）、（９）に基づいて直接グルコースのモル濃度を求めることもできる。
【００８０】
式（８）は、各層における光吸収係数μ１〜μ３を示す式である。
【００８１】
【数８】
【００８２】
式（８）において、μ1（ｉ）（ｉ＝１〜４）は、表皮層における波長λｉ（ｉ＝１〜４）の光吸収係数を示し、μ２（ｉ）（ｉ＝１〜４）は、真皮層における波長λｉ（ｉ＝１〜４）の光吸収係数を示し、μ３（ｉ）（ｉ＝１〜４）は、皮下組織層における波長λｉ（ｉ＝１〜４）の光吸収係数を示す。ε１（ｉ）は、波長λｉ（ｉ＝１〜４）に対する水のモル吸光係数を示し、ε２（ｉ）は、波長λｉ（ｉ＝１〜４）に対するたんぱく質のモル吸光係数を示し、ε３（ｉ）は、波長λｉ（ｉ＝１〜４）に対する脂質のモル吸光係数を示し、ε４（ｉ）は、波長λｉ（ｉ＝１〜４）に対するグルコースのモル吸光係数を示す。ｇ１ｊ、ｇ２ｊ、ｇ３ｊ、ｇ４ｊ（ｊ＝１〜３）は、ｊ層におけるそれぞれ皮膚の主成分である水、たんぱく質、脂質、グルコースのモル濃度を示す。ただし、ｊ＝１は表皮層、ｊ＝２は真皮層、ｊ＝３は皮下組織層をそれぞれ示す。
【００８３】
式（９）は、真皮層（ｊ＝２）における水、たんぱく質、脂質、グルコースのモル濃度ｇ１２、ｇ２２、ｇ３２、ｇ４２を求める式である。式（９）のｇ４２を算出することで、真皮質に含まれるグルコースのモル濃度を求めることができる。
【００８４】
【数９】
【実施例２】
【００８５】
次に、本発明の第２の実施形態について詳しく説明する。
第２の実施形態は、第１の実施形態による血糖値測定装置１００と同じ構成であり、計測光強度取得部１０６、光路長取得部１０７、無吸収時光強度取得部１０８、光吸収係数算出部１０９の動作が異なる。
【００８６】
図８は、血糖値測定装置が血糖値を測定する動作を示す第２のフローチャートである。
まず、血糖値測定装置１００を動作させると、照射部１０４は、皮膚に対して波長λ１の短時間パルス光を照射する（ステップＳ１１）。ここで、波長λ１は、シミュレーション部１０１が伝搬光路長分布及び時間分解波形を算出した複数の波長の中の１つである。
【００８７】
照射部１０４が短時間パルス光を照射すると、受光部１０５は、照射部１０４から照射され、皮膚によって後方散乱した光を受光する（ステップＳ１２）。このとき、受光部１０５は、照射開始からの単位時間毎（例えば、１ピコ秒毎）の受光強度を内部メモリに登録しておく。
受光部１０５が受光を完了すると、計測光強度取得部１０６は、受光部１０５の内部メモリに格納されている受光強度から、ある時刻から時刻τまでの間の受光強度の時間分布を取得する（ステップＳ１３）。時刻τまでの光強度の変化としては、例えば図９に示す時間が選択される。すなわち、表皮を透過する光のスペクトルを精度良く測定したい場合には、１ｐｓ〜１０ｐｓの期間τ１を選択し、真皮を透過する光のスペクトルを精度良く測定したい場合には、２ｐｓ〜３０ｐｓの期間τ２を選択し、皮下組織を透過する光のスペクトルを精度良く検出したい場合には、５０ｐｓ〜７０ｐｓの期間τ３を選択する。τ１、τ２、τ３は、互いに時間軸上では重ならず、皮膚組織の表面に近い層を透過する光のスペクトルを検出しようとするほど、測定開始時間及び測定終了時間は時間軸上で早い時期に設定される。
【００８８】
計測光強度取得部１０６が、時刻τまでの間の受光強度の時間分布を取得すると、光路長取得部１０７は、光路長分布記憶部１０２が記憶する波長λ１の伝搬光路長分布から、ある時刻から時刻τまでの間の皮膚の各層の光路長Ｌ１〜Ｌ３を取得する（ステップＳ１４）。
また、計測光強度取得部１０６が、時刻τまでの間の受光強度を取得すると、無吸収時光強度取得部１０８は、時間分解波形記憶部１０３が記憶する波長λ１の時間分解波形から、ある時刻から時刻τまでの間の無吸収時光強度を取得する（ステップＳ１５）。
【００８９】
光路長取得部１０７が皮膚の各層の光路長を取得し、無吸収時光強度取得部１０８が無吸収時光強度を取得すると、光吸収係数算出部１０９は、式（１０）に基づいて、皮膚の各層の光吸収係数μ１〜μ３を算出する（ステップＳ１６）。ここで、光吸収係数μ１は、表皮層の光吸収係数を示し、光吸収係数μ２は、真皮層の光吸収係数を示し、光吸収係数μ３は、皮下組織層の光吸収係数を示す。
【００９０】
【数１０】
【００９１】
但し、ｌｎ（Ａ）は、Ａの自然対数を示す。また、Ｉ（ｔ）は、時刻ｔにおける受光部１０５の受光強度を示し、Ｉｉｎは、照射部１０４が照射した短時間パルス光の光強度を示す。また、Ｎ（ｔ）は、時間分解波形の時刻ｔにおける無吸収時光強度を示し、Ｎｉｎは、シミュレーション部１０１が照射のシミュレーションを行った光子の個数を示す。また、Ｌ１（ｔ）〜Ｌ３（ｔ）は、時刻ｔにおける皮膚の各層の光路長を示す。
【００９２】
光吸収係数算出部１０９は、皮膚の各層の光吸収係数μ１〜μ３を算出すると、光吸収係数算出部１０９が、皮膚の主成分の種類数と同じ数の波長に対して光吸収係数μ１〜μ３を算出したか否かを判定する（ステップＳ１７）。本実施形態では、皮膚の主成分を水、たんぱく質、脂質、グルコースの４種類として血糖値の測定を行うため、光吸収係数算出部１０９は、４種類の波長λ１〜λ４に対して光吸収係数μ１〜μ３を算出したか否かを判定する。ここで、波長λ１〜λ４は、シミュレーション部１０１が伝搬光路長分布及び時間分解波形を算出した複数の波長の中から選出する。
【００９３】
光吸収係数算出部１０９が、光吸収係数μ１〜μ３を算出していない波長λ１〜λ４があると判定した場合（ステップＳ１７：ＮＯ）、ステップＳ１に戻り、まだ光吸収係数μ１〜μ３を算出していない波長λ１〜λ４の光吸収係数μ１〜μ３の算出を行う。
他方、光吸収係数算出部１０９が、波長λ１〜λ４の光吸収係数μ１〜μ３を算出していると判定した場合（ステップＳ１７：ＹＥＳ）、濃度算出部１１０は、上述した式（７）に基づいて真皮質に含まれるグルコースの濃度を算出する（ステップＳ１８）。
【００９４】
このように、本実施形態によれば、吸収係数μ１〜μ３を、時刻τまでの間の光路長の積分値によって算出する。これにより、計測した受光強度Ｉ（ｔ）に含まれている誤差による吸収係数μ１〜μ３の算出結果に対する影響を少なくすることができる。
【００９５】
以上、図面を参照してこの発明の一実施形態について詳しく説明してきたが、具体的な構成は上述のものに限られることはなく、この発明の要旨を逸脱しない範囲内において様々な設計変更等をすることが可能である。
例えば、実施形態では、濃度定量方法を血糖値測定装置１００に実装し、皮膚の真皮層に含まれるグルコースの濃度を測定する場合を説明したが、これに限られず、濃度定量方法を、複数の光散乱媒質の層から形成される観測対象の任意の層における目的成分の濃度を定量する他の装置に用いても良い。
【００９６】
上述の血糖値測定装置１００は内部に、コンピュータシステムを有している。そして、上述した各処理部の動作は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータが読み出して実行することによって、上記処理が行われる。ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしても良い。
【００９７】
また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。
さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であっても良い。
【符号の説明】
【００９８】
１００…血糖値測定装置、１０１…シミュレーション部、１０２…光路長分布記憶部、１０３…時間分解波形記憶部、１０４…照射部（照射手段）、１０４ａ…照射用光ファイバー、１０５…受光部、１０５ａ…受光用光ファイバー、１０６…計測光強度取得部、１０７…光路長取得部、１０８…無吸収時光強度取得部、１０９…光吸収係数算出部（選択手段）、１１０…濃度算出部
【特許請求の範囲】
【請求項１】
複数の光散乱媒質の層から形成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置であって、
前記観測対象に短時間パルス光を照射する照射手段と、
前記短時間パルス光が前記観測対象によって後方散乱した光を受光する複数の受光部を有する受光手段と、
前記複数の受光部のうち、前記短時間パルス光が前記任意の層によって後方散乱した光を受光する特定の受光部を選択する選択手段と、
前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記特定の受光手段が受光した光の強度を取得する光強度取得手段と、
前記観測対象に対して照射する短時間パルス光の、前記照射手段から前記特定の受光部に至る光の伝搬経路上に配置された複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、
前記観測対象に対して照射し前記特定の受光部において受光する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段と、
前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手段と、
前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手段と、
前記光強度取得手段が取得した光強度と前記光路長取得手段が取得した前記複数の光散乱媒質の層の各々の層の光路長と前記光強度モデル取得手段が取得した光強度モデルとに基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手段と、
前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手段と、
を備えることを特徴とする濃度定量装置。
【請求項２】
前記光強度取得手段は、前記観測対象の層の数ｎ以上となる複数の時刻ｔ１〜ｔｍにおける光強度を取得し、
前記光吸収係数算出手段は、
自然対数を示すｌｎ、前記受光手段が時刻ｔにおいて受光した光強度を示すＩ（ｔ）、前記短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度を示すＮ（ｔ）、前記伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長を示すＬｉ（ｔ）、第ｉ層の光吸収係数を示すμｉを用いて、
【数１】
から任意の層の光吸収係数を算出する、
ことを特徴とする請求項１に記載の濃度定量装置。
【請求項３】
前記照射手段は、複数の波長１〜ｑの光を照射し、
前記光吸収係数算出手段は、前記任意の層における光吸収係数を前記照射手段が照射した複数の波長毎に算出し、
前記濃度算出手段は、
前記任意の層である第ａ層における波長ｉの光吸収係数を示すμａ（ｉ）、前記観測対象を形成する第ｊ成分のモル濃度を示すｇｊ、第ｊ成分の波長ｉに対する光吸収係数を示すεｊ（ｉ）、前記観測対象を形成する主成分の個数を示すｐ、照射手段が照射する波長の種類数を示すｑを用いて、
【数２】
から前記任意の層における前記目的成分の濃度を算出する、
ことを特徴とする請求項２に記載の濃度定量装置。
【請求項４】
前記照射手段から前記観測対象に短時間パルス光を照射する位置を照射位置、前記短時間パルス光が前記観測対象によって後方散乱した光が前記観測対象から前記受光部に向けて出射する位置を受光位置、前記照射位置と前記受光位置との間隔を照射受光間隔、としたときに、前記受光手段は、互いに照射受光間隔の異なる複数の受光部を有し、前記選択手段は、前記照射受光間隔に応じて定まる前記観測対象の内部への短時間パルス光の到達深さに基づいて、前記任意の層に短時間パルス光が伝搬可能な照射受光間隔を有する受光部を選択することを特徴とする請求項３に記載の濃度定量装置。
【請求項５】
前記照射手段は、短時間パルス光を前記観測対象の表面に伝送する照射用光ファイバーを有し、前記受光部は、前記観測対象によって後方散乱した光を伝送する受光用光ファイバーを有し、前記照射用光ファイバーと前記受光用光ファイバーは、両者の光ファイバーコアの中心間隔を所定の照射受光間隔だけ離間して固定するプローブ装置に装着されており、前記プローブ装置の先端部に露出した前記照射用光ファイバーの先端部と前記受光用光ファイバーの先端部とを前記観測対象の表面に接触させることによって、前記観測対象への短時間パルスの照射処理と、前記短時間パルス光が前記観測対象によって後方散乱した光の受光処理とを行うことを特徴とする請求項４に記載の濃度定量装置。
【請求項６】
前記受光部は、互いに等しい照射受光間隔で配置された複数の受光用光ファイバーを有することを特徴とする請求項５に記載の濃度定量装置。
【請求項７】
前記受光部は、互いに等しい照射受光間隔で配置された複数の受光用光ファイバーによって伝送された光を同一受光面上に集光する集光素子を有することを特徴とする請求項６に記載の濃度定量装置。
【請求項８】
前記照射用光ファイバーと前記受光用光ファイバーは、前記照射手段が照射する前記複数の波長１〜ｑの光の波長分散を補償する分散補償型シングルモード光ファイバーであることを特徴とする請求項５〜７のいずれか1項に記載の濃度定量装置。
【請求項９】
前記照射用光ファイバーと前記受光用光ファイバーは、前記照射手段が照射する前記複数の波長１〜ｑの光の波長分散に伴う群遅延時間差が、前記複数の光散乱媒質層のうち最も表面側の層の伝搬光路長分布のピークに対応する伝搬時間よりも短いことを特徴とする請求項５〜７のいずれか1項に記載の濃度定量装置。
【請求項１０】
複数の光散乱媒質の層から形成される観測対象に対して短時間パルス光を照射する照射手段と、前記短時間パルス光が前記観測対象によって後方散乱した光を受光する複数の受光部を有する受光手段と、前記複数の受光部のうち、前記短時間パルス光が前記任意の層によって後方散乱した光を受光する特定の受光部を選択する選択手段と、前記観測対象に対して照射する短時間パルス光の、前記照射手段から前記特定の受光部に至る光の伝搬経路上に配置された複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、前記観測対象に対して照射し前記特定の受光部において受光する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段とを備え、前記観測対象のうち任意の層における目的成分の濃度を定量する濃度定量装置を用いた濃度定量方法であって、
光強度取得手段は、前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記特定の受光部が受光した光の強度を取得し、
光路長取得手段は、前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得し、
光強度モデル取得手段は、前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得し、
光吸収係数算出手段は、前記光強度取得手段が取得した光強度と前記光路長取得手段が取得した前記複数の光散乱媒質の層の各々の層の光路長と前記光強度モデル取得手段が取得した光強度モデルとに基づいて、前記任意の層の光吸収係数を算出し、
濃度算出手段は、前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する、
ことを特徴とする濃度定量方法。
【請求項１１】
複数の光散乱媒質の層から形成される観測対象に対して短時間パルス光を照射する照射手段と、前記短時間パルス光が前記観測対象によって後方散乱した光を受光する複数の受光部を有する受光手段と、前記複数の受光部のうち、前記短時間パルス光が前記任意の層によって後方散乱した光を受光する特定の受光部を選択する選択手段と、前記観測対象に対して照射する短時間パルス光の、前記照射手段から前記特定の受光部に至る光の伝搬経路上に配置された複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、前記観測対象に対して照射し前記特定の受光部において受光する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段とを備え、前記観測対象のうち任意の層における目的成分の濃度を定量する濃度定量装置を、
前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記特定の受光部が受光した光の強度を取得する光強度取得手段、
前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手段、
前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手段、
前記光強度取得手段が取得した光強度と前記光路長取得手段が取得した前記複数の光散乱媒質の層の各々の層の光路長と前記光強度モデル取得手段が取得した光強度モデルとに基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手段、
前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手段、
として動作させるためのプログラム。
【請求項１】
複数の光散乱媒質の層から形成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置であって、
前記観測対象に短時間パルス光を照射する照射手段と、
前記短時間パルス光が前記観測対象によって後方散乱した光を受光する複数の受光部を有する受光手段と、
前記複数の受光部のうち、前記短時間パルス光が前記任意の層によって後方散乱した光を受光する特定の受光部を選択する選択手段と、
前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記特定の受光手段が受光した光の強度を取得する光強度取得手段と、
前記観測対象に対して照射する短時間パルス光の、前記照射手段から前記特定の受光部に至る光の伝搬経路上に配置された複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、
前記観測対象に対して照射し前記特定の受光部において受光する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段と、
前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手段と、
前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手段と、
前記光強度取得手段が取得した光強度と前記光路長取得手段が取得した前記複数の光散乱媒質の層の各々の層の光路長と前記光強度モデル取得手段が取得した光強度モデルとに基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手段と、
前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手段と、
を備えることを特徴とする濃度定量装置。
【請求項２】
前記光強度取得手段は、前記観測対象の層の数ｎ以上となる複数の時刻ｔ１〜ｔｍにおける光強度を取得し、
前記光吸収係数算出手段は、
自然対数を示すｌｎ、前記受光手段が時刻ｔにおいて受光した光強度を示すＩ（ｔ）、前記短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度を示すＮ（ｔ）、前記伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長を示すＬｉ（ｔ）、第ｉ層の光吸収係数を示すμｉを用いて、
【数１】
から任意の層の光吸収係数を算出する、
ことを特徴とする請求項１に記載の濃度定量装置。
【請求項３】
前記照射手段は、複数の波長１〜ｑの光を照射し、
前記光吸収係数算出手段は、前記任意の層における光吸収係数を前記照射手段が照射した複数の波長毎に算出し、
前記濃度算出手段は、
前記任意の層である第ａ層における波長ｉの光吸収係数を示すμａ（ｉ）、前記観測対象を形成する第ｊ成分のモル濃度を示すｇｊ、第ｊ成分の波長ｉに対する光吸収係数を示すεｊ（ｉ）、前記観測対象を形成する主成分の個数を示すｐ、照射手段が照射する波長の種類数を示すｑを用いて、
【数２】
から前記任意の層における前記目的成分の濃度を算出する、
ことを特徴とする請求項２に記載の濃度定量装置。
【請求項４】
前記照射手段から前記観測対象に短時間パルス光を照射する位置を照射位置、前記短時間パルス光が前記観測対象によって後方散乱した光が前記観測対象から前記受光部に向けて出射する位置を受光位置、前記照射位置と前記受光位置との間隔を照射受光間隔、としたときに、前記受光手段は、互いに照射受光間隔の異なる複数の受光部を有し、前記選択手段は、前記照射受光間隔に応じて定まる前記観測対象の内部への短時間パルス光の到達深さに基づいて、前記任意の層に短時間パルス光が伝搬可能な照射受光間隔を有する受光部を選択することを特徴とする請求項３に記載の濃度定量装置。
【請求項５】
前記照射手段は、短時間パルス光を前記観測対象の表面に伝送する照射用光ファイバーを有し、前記受光部は、前記観測対象によって後方散乱した光を伝送する受光用光ファイバーを有し、前記照射用光ファイバーと前記受光用光ファイバーは、両者の光ファイバーコアの中心間隔を所定の照射受光間隔だけ離間して固定するプローブ装置に装着されており、前記プローブ装置の先端部に露出した前記照射用光ファイバーの先端部と前記受光用光ファイバーの先端部とを前記観測対象の表面に接触させることによって、前記観測対象への短時間パルスの照射処理と、前記短時間パルス光が前記観測対象によって後方散乱した光の受光処理とを行うことを特徴とする請求項４に記載の濃度定量装置。
【請求項６】
前記受光部は、互いに等しい照射受光間隔で配置された複数の受光用光ファイバーを有することを特徴とする請求項５に記載の濃度定量装置。
【請求項７】
前記受光部は、互いに等しい照射受光間隔で配置された複数の受光用光ファイバーによって伝送された光を同一受光面上に集光する集光素子を有することを特徴とする請求項６に記載の濃度定量装置。
【請求項８】
前記照射用光ファイバーと前記受光用光ファイバーは、前記照射手段が照射する前記複数の波長１〜ｑの光の波長分散を補償する分散補償型シングルモード光ファイバーであることを特徴とする請求項５〜７のいずれか1項に記載の濃度定量装置。
【請求項９】
前記照射用光ファイバーと前記受光用光ファイバーは、前記照射手段が照射する前記複数の波長１〜ｑの光の波長分散に伴う群遅延時間差が、前記複数の光散乱媒質層のうち最も表面側の層の伝搬光路長分布のピークに対応する伝搬時間よりも短いことを特徴とする請求項５〜７のいずれか1項に記載の濃度定量装置。
【請求項１０】
複数の光散乱媒質の層から形成される観測対象に対して短時間パルス光を照射する照射手段と、前記短時間パルス光が前記観測対象によって後方散乱した光を受光する複数の受光部を有する受光手段と、前記複数の受光部のうち、前記短時間パルス光が前記任意の層によって後方散乱した光を受光する特定の受光部を選択する選択手段と、前記観測対象に対して照射する短時間パルス光の、前記照射手段から前記特定の受光部に至る光の伝搬経路上に配置された複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、前記観測対象に対して照射し前記特定の受光部において受光する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段とを備え、前記観測対象のうち任意の層における目的成分の濃度を定量する濃度定量装置を用いた濃度定量方法であって、
光強度取得手段は、前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記特定の受光部が受光した光の強度を取得し、
光路長取得手段は、前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得し、
光強度モデル取得手段は、前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得し、
光吸収係数算出手段は、前記光強度取得手段が取得した光強度と前記光路長取得手段が取得した前記複数の光散乱媒質の層の各々の層の光路長と前記光強度モデル取得手段が取得した光強度モデルとに基づいて、前記任意の層の光吸収係数を算出し、
濃度算出手段は、前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する、
ことを特徴とする濃度定量方法。
【請求項１１】
複数の光散乱媒質の層から形成される観測対象に対して短時間パルス光を照射する照射手段と、前記短時間パルス光が前記観測対象によって後方散乱した光を受光する複数の受光部を有する受光手段と、前記複数の受光部のうち、前記短時間パルス光が前記任意の層によって後方散乱した光を受光する特定の受光部を選択する選択手段と、前記観測対象に対して照射する短時間パルス光の、前記照射手段から前記特定の受光部に至る光の伝搬経路上に配置された複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、前記観測対象に対して照射し前記特定の受光部において受光する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段とを備え、前記観測対象のうち任意の層における目的成分の濃度を定量する濃度定量装置を、
前記照射手段が短時間パルス光を照射した時刻以降の所定の時刻において前記特定の受光部が受光した光の強度を取得する光強度取得手段、
前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手段、
前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手段、
前記光強度取得手段が取得した光強度と前記光路長取得手段が取得した前記複数の光散乱媒質の層の各々の層の光路長と前記光強度モデル取得手段が取得した光強度モデルとに基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手段、
前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手段、
として動作させるためのプログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【公開番号】特開２０１２−１９８３３（Ｐ２０１２−１９８３３Ａ）
【公開日】平成２４年２月２日（２０１２．２．２）
【国際特許分類】
【出願番号】特願２０１０−１５８０９８（Ｐ２０１０−１５８０９８）
【出願日】平成２２年７月１２日（２０１０．７．１２）
【出願人】（０００００２３６９）セイコーエプソン株式会社 (51,324)
【出願人】（５０４１７３４７１）国立大学法人北海道大学 (971)
【Ｆターム（参考）】
【公開日】平成２４年２月２日（２０１２．２．２）
【国際特許分類】
【出願日】平成２２年７月１２日（２０１０．７．１２）
【出願人】（０００００２３６９）セイコーエプソン株式会社 (51,324)
【出願人】（５０４１７３４７１）国立大学法人北海道大学 (971)
【Ｆターム（参考）】
[ Back to top ]