(57)【特許請求の範囲】
【請求項１】
複数の光散乱媒質の層により構成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置であって、
前記観測対象に、前記任意の層における水の吸収係数の温度変化が小さい特定波長の光を照射する照射手段と、
前記光を照射することにより前記観測対象より放射される複数種の後方散乱光から前記任意の層より放射される後方散乱光を時間分解計測して選択する光散乱媒質層選択手段と、
前記任意の層から放射される後方散乱光を受光する受光手段と、
前記受光手段が受光した光の強度を取得する光強度取得手段と、
前記光強度取得手段が取得した光強度に基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手段と、
前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手段と、を備え、
前記光を短時間パルス光とし、さらに、
前記観測対象に対して照射する前記短時間パルス光の、前記複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、
前記観測対象に対して照射する前記短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段と、
前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手段と、
前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手段とを備え、
前記光強度取得手段は、前記任意の層の複数の時刻ｔ１〜ｔｍにおける光強度を取得し、
前記光吸収係数算出手段は、前記任意の層の光吸収係数を、下記の式（１）
【数１】
0005674093.tif 000022
（但し、Ｉ（ｔ）は前記受光手段が時刻ｔにて受光した光強度、Ｎ（ｔ）は前記短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度、Ｌｉ（ｔ）は前記複数の光散乱媒質の各々の層における伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長、μｉは第ｉ層の光吸収係数である）から算出することを特徴とする濃度定量装置。
【請求項２】
前記光強度取得手段が光強度を取得する複数の時刻は、前記複数の光散乱媒質の各々の層の伝搬光路長分布のピーク時間を含むことを特徴とする請求項１記載の濃度定量装置。
【請求項３】
複数の光散乱媒質の層により構成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置であって、
前記観測対象に、前記任意の層における水の吸収係数の温度変化が小さい特定波長の光を照射する照射手段と、
前記光を照射することにより前記観測対象より放射される複数種の後方散乱光から前記任意の層より放射される後方散乱光を時間分解計測して選択する光散乱媒質層選択手段と、
前記任意の層から放射される後方散乱光を受光する受光手段と、
前記受光手段が受光した光の強度を取得する光強度取得手段と、
前記光強度取得手段が取得した光強度に基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手段と、
前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手段と、を備え、
前記光を短時間パルス光とし、さらに、
前記観測対象に対して照射する前記短時間パルス光の、前記複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、
前記観測対象に対して照射する前記短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段と、
前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手段と、
前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手段とを備え、
前記光強度取得手段は、所定の時刻から少なくとも所定の時間τの間の光強度の時間変化を取得し、
前記光吸収係数算出手段は、前記任意の層の光吸収係数を、下記の式（２）
【数２】
0005674093.tif 000023
（但し、Ｉ（ｔ）は前記受光手段が時刻ｔにて受光した光強度、Ｎ（ｔ）は前記短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度、Ｌｉ（ｔ）は前記複数の光散乱媒質の層各々の層における伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長、ｎは前記観測対象となる層の数、μｉは第ｉ層の光吸収係数である）から算出することを特徴とする濃度定量装置。
【請求項４】
前記濃度算出手段は、前記任意の層における前記目的成分の濃度を、下記の式（３）
【数３】
0005674093.tif 000024
（但し、μａは前記任意の層である第ａ層における光吸収係数、ｇｊは前記観測対象を構成する第ｊ成分のモル濃度、εｊは第ｊ成分の光吸収係数、ｐは前記観測対象を構成する主成分の個数、ｑは前記特定波長の種類数である）から算出することを特徴とする請求項１ないし３のいずれか１項記載の濃度定量装置。
【請求項５】
複数の光散乱媒質の層により構成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量方法であって、
照射手段により、前記観測対象に、前記任意の層における水の吸収係数の温度変化が小さい特定波長の光を照射し、
次いで、光散乱媒質層選択手段により、前記光を照射することにより前記観測対象より放射される複数種の後方散乱光から前記任意の層より放射される後方散乱光を時間分解計測して選択し、
次いで、受光手段により、前記任意の層から放射される後方散乱光を受光し、
次いで、光強度取得手段により、前記受光手段が受光した光の強度を取得し、
次いで、光吸収係数算出手段により、前記光強度取得手段が取得した光強度に基づいて、前記任意の層の光吸収係数を算出し、
次いで、濃度算出手段により、前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する、ことを含み、
前記光を短時間パルス光とし、さらに、
光路長分布記憶手段により、前記観測対象に対して照射する前記短時間パルス光の、前記複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶し、
時間分解波形記憶手段により、前記観測対象に対して照射する前記短時間パルス光の時間分解波形のモデルを記憶し、
光路長取得手段により、前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得し、
光強度モデル取得手段により、前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得しており、
前記光強度取得手段は、前記任意の層の複数の時刻ｔ１〜ｔｍにおける光強度を取得し、
前記光吸収係数算出手段は、前記任意の層の光吸収係数を、下記の式（１）
【数４】
0005674093.tif 000025
（但し、Ｉ（ｔ）は前記受光手段が時刻ｔにて受光した光強度、Ｎ（ｔ）は前記短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度、Ｌｉ（ｔ）は前記複数の光散乱媒質の各々の層における伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長、μｉは第ｉ層の光吸収係数である）から算出することを特徴とする濃度定量方法。
【請求項６】
複数の光散乱媒質の層により構成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量方法であって、
照射手段により、前記観測対象に、前記任意の層における水の吸収係数の温度変化が小さい特定波長の光を照射し、
次いで、光散乱媒質層選択手段により、前記光を照射することにより前記観測対象より放射される複数種の後方散乱光から前記任意の層より放射される後方散乱光を時間分解計測して選択し、
次いで、受光手段により、前記任意の層から放射される後方散乱光を受光し、
次いで、光強度取得手段により、前記受光手段が受光した光の強度を取得し、
次いで、光吸収係数算出手段により、前記光強度取得手段が取得した光強度に基づいて、前記任意の層の光吸収係数を算出し、
次いで、濃度算出手段により、前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する、ことを含み、
前記光を短時間パルス光とし、さらに、
光路長分布記憶手段により、前記観測対象に対して照射する前記短時間パルス光の、前記複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶し、
時間分解波形記憶手段により、前記観測対象に対して照射する前記短時間パルス光の時間分解波形のモデルを記憶し、
光路長取得手段により、前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得し、
光強度モデル取得手段により、前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得しており、
前記光強度取得手段は、所定の時刻から少なくとも所定の時間τの間の光強度の時間変化を取得し、
前記光吸収係数算出手段は、前記任意の層の光吸収係数を、下記の式（２）
【数５】
0005674093.tif 000026
（但し、Ｉ（ｔ）は前記受光手段が時刻ｔにて受光した光強度、Ｎ（ｔ）は前記短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度、Ｌｉ（ｔ）は前記複数の光散乱媒質の層各々の層における伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長、ｎは前記観測対象となる層の数、μｉは第ｉ層の光吸収係数である）から算出することを特徴とする濃度定量方法。
【請求項７】
複数の光散乱媒質の層により構成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置のコンピュータに、
前記観測対象に、前記任意の層における水の吸収係数の温度変化が小さい特定波長の光を照射する照射手順、
前記光を照射することにより前記観測対象より放射される複数種の後方散乱光から前記任意の層より放射される後方散乱光を時間分解計測して選択する光散乱媒質層選択手順、
前記任意の層から放射される後方散乱光を受光する受光手順、
前記受光手順にて受光した光の強度を取得する光強度取得手順、
前記光強度取得手順にて取得した光強度に基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手順、
前記光吸収係数算出手順にて算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手順、を実行させ、
前記光は短時間パルス光であり、さらに、
前記観測対象に対して照射する前記短時間パルス光の、前記複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手順、
前記観測対象に対して照射する前記短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手順、
前記光路長分布記憶手順で記憶した前記伝搬光路長分布のモデルから、前記伝搬光路長分布のモデルの所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手順、
前記時間分解波形記憶手順で記憶した前記短時間パルス光の時間分解波形のモデルから、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手順、を実行させ、
前記光強度取得手順では、前記任意の層の複数の時刻ｔ１〜ｔｍにおける光強度を取得し、
前記光吸収係数算出手順では、前記任意の層の光吸収係数を、下記の式（１）
【数６】
0005674093.tif 000027
（但し、Ｉ（ｔ）は前記受光手順において時刻ｔにて受光した光強度、Ｎ（ｔ）は前記短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度、Ｌｉ（ｔ）は前記複数の光散乱媒質の各々の層における伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長、μｉは第ｉ層の光吸収係数である）から算出することを特徴とするプログラム。
【請求項８】
複数の光散乱媒質の層により構成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置のコンピュータに、
前記観測対象に、前記任意の層における水の吸収係数の温度変化が小さい特定波長の光を照射する照射手順、
前記光を照射することにより前記観測対象より放射される複数種の後方散乱光から前記任意の層より放射される後方散乱光を時間分解計測して選択する光散乱媒質層選択手順、
前記任意の層から放射される後方散乱光を受光する受光手順、
前記受光手順にて受光した光の強度を取得する光強度取得手順、
前記光強度取得手順にて取得した光強度に基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手順、
前記光吸収係数算出手順にて算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手順、を実行させ、
前記光は短時間パルス光であり、さらに、
前記観測対象に対して照射する前記短時間パルス光の、前記複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手順、
前記観測対象に対して照射する前記短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手順、
前記光路長分布記憶手順で記憶した前記伝搬光路長分布のモデルから、前記伝搬光路長分布のモデルの所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手順、
前記時間分解波形記憶手順で記憶した前記短時間パルス光の時間分解波形のモデルから、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手順、を実行させ、
前記光強度取得手順では、所定の時刻から少なくとも所定の時間τの間の光強度の時間変化を取得し、
前記光吸収係数算出手順では、前記任意の層の光吸収係数を、下記の式（２）
【数７】
0005674093.tif 000028
（但し、Ｉ（ｔ）は前記受光手順において時刻ｔにて受光した光強度、Ｎ（ｔ）は前記短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度、Ｌｉ（ｔ）は前記複数の光散乱媒質の層各々の層における伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長、ｎは前記観測対象となる層の数、μｉは第ｉ層の光吸収係数である）から算出することを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の光散乱媒質の層により構成される観測対象のうち、任意の層における目的成分の濃度を、非侵襲的に定量する濃度定量装置及び濃度定量方法並びにプログラムに関するものである。
【背景技術】
【０００２】
近年、我が国は飽食の時代にあって、糖尿病の患者が毎年増加し続けている。そのために、糖尿病性腎炎の患者も毎年増加し続けることとなり、その結果、慢性腎不全の患者も毎年１万人もの増加を続け、患者数は２８万人を超えるようになってきている。
一方、高齢化社会の到来により、予防医学に対する要求の高まりを受けて、個人における代謝量管理の重要性が急速に増大している。中でも、血糖値計測は、糖尿病のごく初期段階での糖代謝の反応を評価することができ、糖尿病の早期診断に基づく早期治療を可能にしている。
【０００３】
従来、血糖値の測定は、腕あるいは指先等の静脈から採血を行い、この血液中のグルコースに対する酵素活性を測定することで行っている。しかし、このような血糖値の測定方法では、採血が煩雑であり、しかも採血に痛みを伴い、さらには感染症の危険性を伴う等の様々な問題がある。
また、血糖値を連続的に測定する方法としては、静脈に注射針を刺した状態で連続的に血糖値相応のグルコースの定量を行う機器が米国にて開発されており、現在臨床試験中である。しかし、静脈に注射針を刺したままにしているために、血糖値の測定中に針が抜ける危険性や感染症の危険性がある。
そこで、採血無しに頻繁に血糖値を測定することができ、しかも感染症の危険性が無い血糖値の測定装置の開発が求められている。さらには、簡単にかつ常時装着可能であり、小型化可能な血糖値の測定装置の開発が求められている。
【０００４】
そこで、血糖値の測定装置に分子吸光の原理を用いた分光分析装置を適用することにより、非侵襲的に血糖値を測定する装置が提案されている（例えば、特許文献１参照）。
この装置は、皮膚に近赤外の連続光を照射し、その光吸収量からグルコースの濃度を算出する装置である。具体的には、予めグルコース濃度と照射する近赤外光の波長と光の吸収量との関係を示す検量線を作成しておき、皮膚に近赤外の連続光を照射し、この皮膚からの戻り光をモノクロメーター等を用いてある波長域を走査し、その波長域の各波長に対する光の吸収量を求め、この各波長における光の吸収量と検量線とを比較することで、血液中のグルコース濃度、すなわち血糖値を算出している。
【０００５】
一般に、水溶液や含水率の高い試料の近赤外分光分析を行う場合、それらのスペクトルは、水のスペクトルと同様、温度変化に伴うスペクトルのシフト等の変動が大きい。したがって、近赤外分光を用いて定量分析をする場合、水溶液や試料の温度の影響を無視することができない。
そこで、生体表面近傍の組織中のグルコース濃度を近赤外領域における光の吸収を利用して測定する場合に、近赤外光送受用の光ファイババンドルのプロ−ブ先端の測定面と生体の表面近傍組織との接触部分の温度を、ヒータ及び表面温度検知手段を用いて一定にする装置も提案されている（特許文献２参照）。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特許第３９３１６３８号公報
【特許文献２】特開２００１−２９９７２７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかしながら、従来の近赤外光の吸収量から血液中あるいは生体の表面近傍組織のグルコース濃度を測定する方法においては、血液中あるいは生体中に含まれる水の近赤外光に対する吸収係数の温度変化率が大きく、血液中あるいは生体中のグルコース濃度を精度よく測定することが難しいという問題点があった。
例えば、近赤外光送受用の光ファイババンドルのプロ−ブ先端の測定面と生体の表面近傍組織との接触部分の温度を、ヒータ及び表面温度検知手段を用いて一定にすれば、確かに、接触部分の温度は一定になるが、生体自体の温度が変化した場合、この生体の表面近傍組織の温度も変化して近赤外光に対する吸収係数が変化してしまい、やはり、生体中のグルコース濃度を精度よく測定することは難しい。
【０００８】
本発明は、上記の課題を解決するためになされたものであって、皮膚等の被測定物である観測対象自体の温度が変化した場合においても、この観測対象に含まれる水の近赤外光に対する吸収係数の温度変化率を小さくすることで、この観測対象に含まれる目的成分の濃度を、非侵襲的に精度良く測定することが可能な濃度定量装置及び濃度定量方法並びにプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
上記の課題を解決するために、本発明は以下の濃度定量装置及び濃度定量方法並びにプログラムを採用した。
すなわち、本発明の濃度定量装置は、複数の光散乱媒質の層により構成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置であって、前記観測対象に、前記任意の層における水の吸収係数の温度変化が小さい特定波長の光を照射する照射手段と、前記光を照射することにより前記観測対象より放射される複数種の後方散乱光から前記任意の層より放射される後方散乱光を選択する光散乱媒質層選択手段と、前記任意の層から放射される後方散乱光を受光する受光手段と、前記受光手段が受光した光の強度を取得する光強度取得手段と、前記光強度取得手段が取得した光強度に基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手段と、前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手段と、を備えてなることを特徴とする。
【００１０】
本発明の濃度定量装置では、照射手段により、前記観測対象に、前記任意の層における水の吸収係数の温度変化が小さい特定波長の光を照射し、光散乱媒質層選択手段により、前記光を照射することにより前記観測対象より放射される複数種の後方散乱光から前記任意の層より放射される後方散乱光を選択し、受光手段により、前記任意の層から放射される後方散乱光を受光する。
このように、観測対象に照射される光を、水の吸収係数の温度変化が小さい特定波長の光とすることで、任意の層から放射される後方散乱光における水の影響を小さくすることができ、この後方散乱光を基に算出される観測対象の任意の層における目的成分の濃度においても、水の影響を小さくすることができる。したがって、目的成分の濃度における水の影響を低減することができ、目的成分の濃度を、非侵襲的に精度良く測定することができる。
【００１１】
本発明の濃度定量装置は、前記光を短時間パルス光とし、さらに、前記観測対象に対して照射する前記短時間パルス光の、前記複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、前記観測対象に対して照射する前記短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段と、前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手段と、前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手段とを備え、
前記光強度取得手段は、前記任意の層の複数の時刻ｔ１〜ｔｍにおける光強度を取得し、
前記光吸収係数算出手段は、前記任意の層の光吸収係数を、下記の式（１）
【数１】
0005674093.tif 000002
（但し、Ｉ（ｔ）は前記受光手段が時刻ｔにて受光した光強度、Ｎ（ｔ）は前記短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度、Ｌｉ（ｔ）は前記複数の光散乱媒質の各々の層における伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長、μｉは第ｉ層の光吸収係数である）
から算出することを特徴とする。
【００１２】
本発明の濃度定量装置では、光強度取得手段が、任意の層の複数の時刻ｔ１〜ｔｍにおける光強度を取得し、光吸収係数算出手段が、任意の層の光吸収係数を、上記の式（１）から算出する。
このように、後方散乱光を時間分解計測することで、任意の層以外の層からの後方散乱光をノイズとして低減することができ、目的成分の濃度における水の影響を低減することができる。したがって、目的成分の濃度をさらに精度良く測定することができる。
【００１３】
本発明の濃度定量装置は、前記光強度取得手段が光強度を取得する複数の時刻は、前記複数の光散乱媒質の各々の層の伝搬光路長分布のピーク時間を含むことを特徴とする。
【００１４】
本発明の濃度定量装置では、光強度取得手段が光強度を取得する複数の時刻が、複数の光散乱媒質の各々の層の伝搬光路長分布のピーク時間を含むことにより、観測対象中の複数の層から任意の層を効率的に選択することができる。したがって、任意の層における目的成分の濃度をさらに精度良く測定することができる。
【００１５】
本発明の濃度定量装置は、前記光を短時間パルス光とし、さらに、前記観測対象に対して照射する前記短時間パルス光の、前記複数の光散乱媒質の層の各々の層における伝搬光路長分布のモデルを記憶する光路長分布記憶手段と、前記観測対象に対して照射する前記短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶手段と、前記光路長分布記憶手段から、前記伝搬光路長分布のモデルの前記所定の時刻における、前記複数の光散乱媒質の層の各々の層の光路長を取得する光路長取得手段と、前記時間分解波形記憶手段から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光の強度を取得する光強度モデル取得手段とを備え、
前記光強度取得手段は、所定の時刻から少なくとも所定の時刻τの間の光強度を取得し、
前記光吸収係数算出手段は、前記任意の層の光吸収係数を、下記の式（２）
【数２】
0005674093.tif 000003
（但し、Ｉ（ｔ）は前記受光手段が時刻ｔにて受光した光強度、Ｎ（ｔ）は前記短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度、Ｌｉ（ｔ）は前記複数の光散乱媒質の層各々の層における伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長、ｎは前記観測対象となる層の数、μｉは第ｉ層の光吸収係数である）
から算出することを特徴とする。
【００１６】
本発明の濃度定量装置では、光強度取得手段が、所定の時刻から少なくとも所定の時刻τの間の光強度の時間変化を取得し、光吸収係数算出手段が、任意の層の光吸収係数を、上記の式（２）から算出する。
このように、後方散乱光を時間分解計測することで、任意の層以外の層からの後方散乱光をノイズとして低減することができ、目的成分の濃度における水の影響を低減することができる。したがって、目的成分の濃度をさらに精度良く測定することができる。
【００１７】
本発明の濃度定量装置は、前記濃度算出手段は、前記任意の層における前記目的成分の濃度を、下記の式（３）
【数３】
0005674093.tif 000004
（但し、μａは前記任意の層である第ａ層における光吸収係数、ｇｊは前記観測対象を構成する第ｊ成分のモル濃度、εｊは第ｊ成分の光吸収係数、ｐは前記観測対象を構成する主成分の個数、ｑは前記特定波長の種類数である）
から算出することを特徴とする。
【００１８】
本発明の濃度定量装置では、濃度算出手段が、任意の層における目的成分の濃度を、上記の式（３）から算出する。
このように、時間分解計測した後方散乱光を用いて任意の層における目的成分の濃度を算出することで、目的成分の濃度における水の影響を低減することができる。したがって、目的成分の濃度をさらに精度良く測定することができる。
【００１９】
本発明の濃度定量方法は、複数の光散乱媒質の層により構成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量方法であって、
照射手段により、前記観測対象に、前記任意の層における水の吸収係数の温度変化が小さい特定波長の光を照射し、次いで、光散乱媒質層選択手段により、前記光を照射することにより前記観測対象より放射される複数種の後方散乱光から前記任意の層より放射される後方散乱光を選択し、次いで、受光手段により、前記任意の層から放射される後方散乱光を受光し、次いで、光強度取得手段により、前記受光手段が受光した光の強度を取得し、次いで、光吸収係数算出手段により、前記光強度取得手段が取得した光強度に基づいて、前記任意の層の光吸収係数を算出し、次いで、濃度算出手段により、前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する、ことを特徴とする。
【００２０】
本発明の濃度定量方法では、照射手段により、前記観測対象に、前記任意の層における水の吸収係数の温度変化が小さい特定波長の光を照射し、次いで、光散乱媒質層選択手段により、前記光を照射することにより前記観測対象より放射される複数種の後方散乱光から前記任意の層より放射される後方散乱光を選択し、次いで、受光手段により、前記任意の層から放射された後方散乱光を受光する。
このように、観測対象に照射される光を、任意の層における水の吸収係数の温度変化が小さい特定波長の光とすることで、任意の層から放射された後方散乱光における水の影響を小さくすることができ、この後方散乱光を基に算出された観測対象の任意の層における目的成分の濃度においても、水の影響を小さくすることができる。したがって、目的成分の濃度における水の影響を低減することができ、目的成分の濃度を、非侵襲的に精度良くかつ効率良く測定することができる。
【００２１】
本発明のプログラムは、複数の光散乱媒質の層により構成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置のコンピュータに、前記観測対象に、前記任意の層における水の吸収係数の温度変化が小さい特定波長の光を照射する照射手順、前記光を照射することにより前記観測対象より放射される複数種の後方散乱光から前記任意の層より放射される後方散乱光を選択する光散乱媒質層選択手順、前記任意の層から放射される後方散乱光を受光する受光手順、前記受光手段が受光した光の強度を取得する光強度取得手順、前記光強度取得手段が取得した光強度に基づいて、前記任意の層の光吸収係数を算出する光吸収係数算出手順、前記光吸収係数算出手段が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手順、を実行させることを特徴とする。
【００２２】
本発明のプログラムでは、濃度定量装置のコンピュータに、前記観測対象に、前記任意の層における水の吸収係数の温度変化が小さい特定波長の光を照射する照射手順、前記光を照射することにより前記観測対象より放射される複数種の後方散乱光から前記任意の層より放射される後方散乱光を選択する光散乱媒質層選択手順、前記任意の層から放射された後方散乱光を受光する受光手順、を順次実行させる。
このように、任意の層における水の吸収係数の温度変化が小さい特定波長の光を照射する照射手順を実行することで、任意の層から放射された後方散乱光における水の影響を小さくすることができ、この後方散乱光を基に算出された観測対象の任意の層における目的成分の濃度においても、水の影響を小さくすることができる。したがって、目的成分の濃度における水の影響を低減することができ、目的成分の濃度を、非侵襲的に精度良くかつ効率良く測定することができる。
【図面の簡単な説明】
【００２３】
【図１】本発明の第１の実施形態の血糖値測定装置の構成を示す概略ブロック図である。
【図２】皮膚の断面構造を示す模式図である。
【図３】シミュレーション部が算出した各層の伝搬光路長分布を示す図である。
【図４】シミュレーション部が算出した時間分解波形を示す図である。
【図５】水による光吸収波長特性を示す図である。
【図６】皮膚の主成分の吸収スペクトルを示す図である。
【図７】皮膚の皮下組織、真皮層及び表皮層各々に照射される光の波長と吸収係数との関係を示す図である。
【図８】水の吸光度スペクトルの温度依存性を示す図である。
【図９】水の吸光度スペクトルの差の温度依存性を示す図である。
【図１０】グルコース水溶液の吸光度スペクトルの一例を示す図である。
【図１１】本発明の第１の実施形態の血糖値を測定する手順を示すフローチャートである。
【図１２】本発明の第２の実施形態の血糖値測定装置の構成を示す概略ブロック図である。
【図１３】本発明の第２の実施形態の血糖値を測定する手順を示すフローチャートである。
【発明を実施するための形態】
【００２４】
本発明の濃度定量装置及び濃度定量方法並びにプログラムを実施するための形態について説明する。
本発明では、濃度定量装置として血糖値測定装置を、観測対象として人の手のひらの皮膚を、目的成分としてグルコースを、特定波長の光として特定波長の短時間パルス光を、それぞれ例に取り説明する。
【００２５】
［第１の実施形態］
図１は、本発明の第１の実施形態の血糖値測定装置の構成を示す概略ブロック図である。
この血糖値測定装置１は、手のひら等の皮膚（観測対象）を構成する複数層のうちの真皮層（任意の層）に含まれるグルコース（目的成分）の濃度を非侵襲にて定量する装置であり、シミュレーション部２と、光路長分布記憶部（光路長分布記憶手段）３と、時間分解波形記憶部（時間分解波形記憶手段）４と、照射部（照射手段）５と、導光部６と、光散乱媒質層選択部（光散乱媒質層選択手段）７と、受光部（受光手段）８と、光強度取得部（光強度取得手段）９と、光路長取得部（光路長取得手段）１０と、無吸収時光強度取得部（光強度モデル取得手段）１１と、光吸収係数算出部（光吸収係数算出手段）１２と、濃度算出部（濃度算出手段）１３とを備えている。
【００２６】
シミュレーション部２は、光吸収係数がゼロの皮膚モデルに対して光を照射するシミュレーションを行う。
光路長分布記憶部３は、皮膚に対して照射する短時間パルス光の、この皮膚を構成する各々の層における伝搬光路長分布のモデルを記憶する。ここでは、光吸収係数がゼロの皮膚モデルの伝搬光路長分布を記憶する。
時間分解波形記憶部４は、皮膚に対して照射する短時間パルス光の時間分解波形のモデルを記憶する。ここでは、光吸収係数がゼロの皮膚モデルの時間分解波形を記憶する。
照射部５は、皮膚に、この皮膚を構成する真皮層における水の吸収係数の温度変化が小さい特定波長λｋ（λ１、λ２…）の短時間パルス光を皮膚に対して照射する。ここで、短時間パルス光とは、パルス幅が１００ｐｓｅｃ程度かそれ以下のパルス光を意味する。なお、短時間パルス光として０．１ｐｓｅｃから数ｐｓｅｃの範囲のパルス幅を持つパルス光を用いても良い。
【００２７】
導光部６は、この特定波長λｋの短時間パルス光を皮膚に照射することにより、この皮膚から放射される複数種の後方散乱光を集光し、光散乱媒質層選択部７へ導光する。
光散乱媒質層選択部７は、導光部６により集光されかつ導光された皮膚から放射される複数種の後方散乱光から、真皮層により放射される後方散乱光を入出射間距離と生体内の深さ方向の到達距離との関係により選択する。
受光部８は、短時間パルス光が皮膚によって後方散乱した光を受光する。
【００２８】
光強度取得部９は、受光部８が受光した真皮層から放射される後方散乱光の異なる複数の時刻の受光強度を取得する。
ここで、この複数の時刻は、皮膚（観測対象）を構成する各々の層の伝搬光路長分布のピーク時間を含むことが好ましい。
このように、各々の層の伝搬光路長分布のピーク時間を含むことで、皮膚の複数の層から任意の層、例えば真皮層を効率的に選択することができる。
【００２９】
光路長取得部１０は、光路長分布記憶部３から、伝搬光路長分布のモデルの所定の時刻における、皮膚の各々の層の光路長を取得する。ここでは、光路長分布記憶部３からある時刻における光路長を取得する。
無吸収時光強度取得部１１は、時間分解波形記憶部４から、短時間パルス光の時間分解波形のモデルの所定の時刻における光の強度を取得する。ここでは、時間分解波形記憶部４からある時刻における光強度を取得する。
光吸収係数算出部１２は、特定波長λｋの短時間パルス光を照射した皮膚の真皮層における光吸収係数を算出する。
【００３０】
この光吸収係数算出部１２では、皮膚における任意の層の光吸収係数を、下記の式（４）
【数４】
0005674093.tif 000005
（但し、Ｉ（ｔ）は受光部８が時刻ｔにて受光した光強度、Ｎ（ｔ）は特定波長λｋの短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度、Ｌｉ（ｔ）は皮膚の各々の層における伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長、μｉは第ｉ層の光吸収係数である）
から算出する。
ここで、第１層は表皮層、第２層は真皮層、第３層は皮下組織を示し、μ１は表皮層の光吸収係数、μ２は真皮層の光吸収係数、μ３は皮下組織の光吸収係数を示す。
【００３１】
濃度算出部１３は、真皮層における光吸収係数から、真皮層に含まれるグルコースの濃度を算出する。
この濃度算出部１３では、皮膚の任意の層におけるグルコースの濃度を、下記の式（５）
【数５】
0005674093.tif 000006
（但し、μａは皮膚の任意の層である第ａ層における光吸収係数、ｇｊは皮膚を構成する第ｊ成分のモル濃度、εｊは第ｊ成分の光吸収係数、ｐは皮膚を構成する主成分の個数、ｑは特定波長λｋの種類数である）
から算出する。
ここで、第１層は表皮層、第２層は真皮層、第３層は皮下組織を示し、μ１は表皮層の光吸収係数、μ２は真皮層の光吸収係数、μ３は皮下組織の光吸収係数を示す。
【００３２】
この血糖値測定装置１では、照射部５は、皮膚２１に、この皮膚２１を構成する真皮層２３における水の吸収係数の温度変化が小さい特定波長λｋ、例えば、λ１＝１４４５ｎｍ、λ２＝１７８２ｎｍの波長の短時間パルス光を照射する。導光部６は、皮膚２１から放射される複数種の後方散乱光を集光する。光散乱媒質層選択部７は、皮膚２１から放射される複数種の後方散乱光から、真皮層２３により放射される後方散乱光を入出射間距離と生体内の深さ方向の到達距離との関係により選択する。受光部８は、真皮層２３から放射される後方散乱光をより多く含む信号光を受光する。
【００３３】
さらに、光強度取得部９は、時刻ｔにおいて受光部８が受光した真皮層２３から放射される後方散乱光の光強度を取得する。
一方、光路長取得部１０は、光路長分布記憶部３から、皮膚モデルにおける伝搬光路長分布の時刻ｔにおける皮膚の各層の光路長を取得し、無吸収時光強度取得部１１は、時間分解波形記憶部４から、皮膚モデルにおける短時間パルス光の時間分解波形の時刻ｔにおける光の強度を取得する。
【００３４】
次いで、光吸収係数算出部１２は、光強度取得部９が取得した光強度と、光路長取得部１０が取得した皮膚の各層の光路長と、無吸収時光強度取得部１１が取得した光強度とに基づいて、皮膚の真皮層の光吸収係数を算出する。次いで、濃度算出部１３は、光吸収係数算出部１２が算出した光吸収係数に基づいて、真皮層２３に含まれるグルコースの濃度を、上記の式（５）に基づき算出する。
これにより、真皮層以外の層によるノイズの影響を軽減して、真皮層に含まれるグルコースの濃度を算出することができる。
【００３５】
以上により、真皮層２３から放射される後方散乱光における水の影響を小さくすることができ、この後方散乱光を基に算出される真皮層２３に含まれるグルコースの濃度についても、水の影響を小さくすることができる。したがって、グルコース濃度における水の影響を低減することができ、真皮層２３に含まれるグルコース濃度を、非侵襲的に精度良く測定することができる。
【００３６】
次に、血糖値測定装置１の動作を説明する。
血糖値測定装置１は、血糖値を測定する前に、予め皮膚モデルの各層における伝搬光路長分布と時間分解波形とを算出しておく必要がある。
図２は、人の皮膚組織の断面を示す模式図であり、皮膚２１は、概ね水を２０％程度含み、残部が蛋白質からなる厚み０．３ｍｍ程度の表皮層２２と、表皮層２２下に形成され、概ね水を６０％程度、蛋白質、脂質及びグルコースを含有する厚み１．２ｍｍ程度の真皮層（任意の層）２３と、真皮層２３下に形成され、概ね脂質を９０％以上含み、残部が水からなる厚み３．０ｍｍ程度の皮下組織２４とにより構成されている。
【００３７】
ここで、皮膚モデルの伝搬光路長分布及び時間分解波形の算出方法を説明する。
初めに、シミュレーション部２は、皮膚モデルを生成する。皮膚モデルの生成は、皮膚の各層の光散乱係数、光吸収係数及び厚みを決定することで行う。ここで、皮膚の各層の散乱係数及び厚みは、個体による差が少ないので、予めサンプルを取ることなどによって決定すると良い。
また、ここで用いる皮膚モデルの光吸収係数はゼロとする。その理由は、この皮膚モデルを用いて光吸収量を算出するからである。
【００３８】
シミュレーション部２は、皮膚モデルを生成すると、この皮膚モデルに光を照射するシミュレーションを行う。このとき、照射部５の位置と受光部８の位置との間の距離を決定しておく必要がある。シミュレーションは、モンテカルロ法を用いて行うと良い。モンテカルロ法によるシミュレーションは、例えば以下のように行われる。
【００３９】
まず、シミュレーション部２は、照射する光のモデルを光子（光束）とし、この光子を皮膚モデルに照射する計算を行う。皮膚モデルに照射された光子は、皮膚モデル内を移動する。このとき、光子は、次に進む点までの距離Ｌ及び方向θを乱数Ｒによって決定する。シミュレーション部２は、光子が次に進む点までの距離Ｌの計算を、式（６）により行う。
【数６】
0005674093.tif 000007
ただし、μｓは、皮膚モデルの第ｓ層（表皮層、真皮層、皮下組織層の何れか）の散乱係数を示す。
【００４０】
また、シミュレーション部２は、光子が次に進む点までの方向θの計算を、式（７）により行う。
【数７】
0005674093.tif 000008
ただし、ｇは、散乱角度の余弦（ｃｏｓ）の平均である非等方性パラメータを示し、皮膚の非等方性パラメータは、略０．９である。
【００４１】
シミュレーション部２は、上記式（６）及び式（７）の計算を単位時間毎に繰り返すことにより、照射部５から受光部８までの光子の移動経路を算出することができる。シミュレーション部２は、複数の光子について移動距離の算出を行う。例えば、シミュレーション部２は、１０８個の光子について移動距離を算出する。
【００４２】
図３は、シミュレーション部が算出した各層の伝搬光路長分布を示す図である。
図３では、横軸を光子の照射からの経過時間とし、縦軸を光路長の対数表示としている。
シミュレーション部２は、受光部８に到達した光子の各々の移動経路を、移動経路が通過する層毎に分類する。そして、シミュレーション部２は、単位時間毎に到達した光子の移動経路の平均長を分類された層毎に算出することで、図３に示すような皮膚の各層の伝搬光路長分布を算出する。
【００４３】
図４は、シミュレーション部が算出した時間分解波形を示す図である。
図４では、横軸を光子の照射からの経過時間とし、縦軸を受光部８が検出した光子数としている。
シミュレーション部２は、単位時間毎に受光部８に到達した光子の個数を算出することで、図４に示すような皮膚モデルの時間分解波形を算出する。
上述したような処理により、シミュレーション部２は、複数の波長に対して、皮膚モデルの伝搬光路長分布及び時間分解波形を算出する。このとき、シミュレーション部２は、皮膚の主成分（水、たんぱく質、脂質、グルコース等）の吸収スペクトルの差が大きい波長について伝搬光路長分布及び時間分解波形を算出すると良い。
【００４４】
図５は、水による光吸収波長特性を示す図である（久保宇市著、「医用レーザ入門」、第１版、オーム社、昭和６０年６月２５日発行、第７０頁、ＩＳＢＮ４−２７４−０３０６５−２）。
図５では、横軸を照射する光の波長（μｍ）とし、縦軸を照射する光の皮膚への浸透深さ（ｃｍ）とし、水に向かって光を入射した場合、入射時の光強度が１／１０に減少するまでに進む浸透深さを赤外域の各波長の光に対して示している。
例えば、３．０μｍ付近の波長帯域の光では、浸透深さが２×１０−３ｃｍ程度と浅く、水に吸収され易いことが分かる。
【００４５】
図６は、皮膚の主成分の吸収スペクトルを示すグラフである。この図６では、横軸を照射する光の波長とし、縦軸を吸収係数としている。
図６によれば、グルコースの吸収係数は波長が１６００ｎｍのときに極大となり、水の吸収係数は波長が１４５０ｎｍのときに極大となることがわかる。
したがって、シミュレーション部２は、例えば１４００ｎｍ、１４５０ｎｍ、１５００ｎｍ、１６００ｎｍ、１６８０ｎｍ、１７２０ｎｍ、１７４０ｎｍというように皮膚の主成分の吸収スペクトルの差が大きい波長について伝搬光路長分布及び時間分解波形を算出すると良い。
【００４６】
シミュレーション部２は、複数の波長に対する皮膚モデルの伝搬光路長分布及び時間分解波形を算出すると、伝搬光路長分布の情報を光路長分布記憶部３に記憶させ、時間分解波形の情報を時間分解波形記憶部４に記憶させる。
【００４７】
図７は、皮膚２１の表皮層２２、真皮層２３及び皮下組織２４各々に照射される光の波長と吸収係数との関係を示す図であり、図中、Ａは表皮層２２の吸収係数を、Ｂは真皮層２３の吸収係数を、Ｃは皮下組織２４の吸収係数を、それぞれ示している。
この図によれば、真皮層２３の吸収スペクトルには、波長１４５０ｎｍ付近に極大値があり、その吸収係数値は水の吸収係数の６０％程度の値なので、真皮層２３の吸収の６０％は水分によるものと考えられる。また、表皮層２２の吸収スペクトルにおいても、波長１４５０ｎｍ付近に真皮層２３の１／３程度の大きさの極大値があり、その吸収係数値は水の吸収係数の２０％程度の値なので、表皮層２２の吸収の２０％は水分によるものと考えられる。一方、皮下組織２４の吸収スペクトルでは、波長１４５０ｎｍ付近に真皮層２３の１／１０程度の大きさの極大値しかなく、その吸収係数値は水の吸収係数の数％程度の値なので、皮下組織２４の吸収の数％程度が水分によるものと考えられる。
【００４８】
以上により、真皮層２３の吸収の６０％は水分によるものと考えられ、また、表皮層２２の吸収の２０％は水分によるものと考えられるが、皮下組織２４では、その吸収の数％程度が水分によるものと考えられる。したがって、皮膚から血糖値を非侵襲的に測定するには、測定対象としてグルコースを含んでいる真皮層２３を選択し、この真皮層２３に含まれるグルコース量を測定すればよいことが分かる。
【００４９】
ところで、水の吸収係数には温度依存性があることが知られている。
図８は、水の吸光度スペクトルの温度依存性を示す図であり、図中、Ａは４１℃における水の吸光度スペクトル、Ｂは２１℃における水の吸光度スペクトルである。
ここでは、セル長が０．５ｍｍの光学セルを用い、光学セルホルダとして温調ユニットタイプのものを用い、恒温循環槽を用いて±０．１℃の範囲で温度調節を行い、紫外可視近赤外分光光度計 Ｌａｍｂｄａ ９００Ｓ（パーキンエルマー社製）を用いて４１℃及び２１℃各々における水の吸光度スペクトルを測定した。
図８によれば、水の吸光度スペクトルの極大値は、２１℃では波長１４５０ｎｍ付近にあり、温度が２１℃より高くなるにしたがって、極大値が１４５０ｎｍより短波長側にシフトすることが分かる。
【００５０】
そこで、２１℃における水の吸光度スペクトルを基準として、２１℃以外の温度における水の吸光度スペクトルの基準との差を求め、この差が「０」となる点の波長を求めれば、この波長が温度変化の影響を受けない波長となる。
図９は、水の吸光度スペクトルの差の温度依存性を示す図であり、図中、Ａは２５℃における水の吸光度スペクトルと２１℃における水の吸光度スペクトルとの差を、Ｂは３１℃における水の吸光度スペクトルと２１℃における水の吸光度スペクトルとの差を、Ｃは３７℃における水の吸光度スペクトルと２１℃における水の吸光度スペクトルとの差を、Ｄは４１℃における水の吸光度スペクトルと２１℃における水の吸光度スペクトルとの差を、それぞれ示している。
【００５１】
図９によれば、各温度差における水の吸光度スペクトルの差の極大値は、温度差が大きくなるにしたがって短波長側にシフトし、また、差の極小値は、温度差が大きくなるにしたがって長波長側にシフトしているが、水の吸光度スペクトルの差が「０」となる点Ｐ及びＱでは、温度差の大小にかかわらず波長が一定であることが分かる。この場合の点Ｐにおける波長は１４４５ｎｍであり、点Ｑにおける波長は１７８２ｎｍである。したがって、１４４５ｎｍ、１７８２ｎｍのいずれかの波長を特定波長λｋとした光を用いて試料の吸光度スペクトルを測定すれば、得られた吸光度スペクトルには温度変化の影響が無いことになる。
【００５２】
一例として、１４４５ｎｍ、１７８２ｎｍの２種類の波長について、血糖値を１００ｍｇ／ｄｌ、測定精度を±５％としたときの温度変化に対して安定とみなせる吸光度及び吸収係数の範囲を実験により求めたところ、１４４５ｎｍの波長では１．３５±０．００００１７５、±０．００００４／ｍｍ、１７８２ｎｍの波長では０．３９±０．００００１７５、±０．００００４／ｍｍであった。したがって、これらの吸光度及び吸収係数の範囲から波長の範囲を求めると、１４４５±０．０２５ｎｍ、１７８２±０．１ｎｍとなった。
その結果、血糖値を１００ｍｇ／ｄｌ、測定精度を±５％としたときの温度変化に対して安定とみなせる波長範囲は、１４４５±０．０２５ｎｍ、１７８２±０．１ｎｍであることが分かった。
【００５３】
以上により、真皮層における水の吸収係数の温度変化が小さい特定波長λｋの短時間パルス光を用いて真皮層１３に含まれるグルコース量を測定すれば、皮膚から血糖値を非侵襲的にて測定することができる。
【００５４】
図１０は、グルコース水溶液の吸光度スペクトルの一例を示す図であり、図中、Ａは参照側を蒸留水（２１．５℃）として測定した９．４ｇ／ｄｌの高濃度のグルコース水溶液の吸光度スペクトルの測定値を、Ｂは同グルコース水溶液の吸光度スペクトルの測定値を温度補正及び体積補正した補正値を、それぞれ示している。
【００５５】
このグルコース水溶液では、グルコース濃度を正常値の約１００倍である９．４ｇ／ｄｌとしたので、このときの体積増加は約６％である。
また、この濃度では、グルコースと水の体積比率が６：１００と大きく、無視できない。このように、試料側の水の体積が参照側の水の体積と比べて減少しているので、水の吸光度の大きい１４００〜１５００ｎｍ付近と１９００ｎｍ以上の波長領域では、吸光度スペクトル差が大きく負になっている。この体積減少は、セル長１ｍｍに対して０．０５７ｍｍの減少に相当している。そこで、グルコース水溶液の吸光度スペクトルの測定値に、この体積減少及び温度の補正を行うと、図１０中、Ｂに示す補正値となり、固体のグルコースの吸光度スペクトルに近似したものとなる。
このように、グルコース水溶液の吸光度スペクトルの測定値を温度補正及び体積補正することにより、グルコース単体の吸光度スペクトルに近いものが得られる。
【００５６】
次に、この血糖値測定装置１を用いて血糖値を測定する手順について、図１１に基づき説明する。
まず、被測定者が血糖値測定装置１を手首等の皮膚に当て、測定開始スイッチ（図示せず）の押下等により血糖値測定装置１を動作させる。
ここでは、照射部５が、皮膚２１に対して、この皮膚２１を構成する真皮層２３における水の吸収係数の温度変化が小さい特定波長λｋの短時間パルス光を照射する（ステップＳ１）。
この特定波長λｋとしては、例えば、図９にて得られた１４４５ｎｍ、１７８２ｎｍのいずれか、または双方の波長を用いることが好ましい。
【００５７】
次いで、導光部６により、皮膚２１から放射される複数種の後方散乱光、すなわち皮下組織２２、真皮層２３及び表皮層２４各々から放射される後方散乱光を集光し、光散乱媒質層選択部７へ導光する。
光散乱媒質層選択部７では、導光部６により集光されかつ導光された皮下組織２２、真皮層２３及び表皮層２４各々から放射される後方散乱光から、真皮層２３により放射される後方散乱光をより多く含む信号光を選択する（ステップＳ２）。
【００５８】
次いで、受光部８により、真皮層２３から放射される単位時間毎の後方散乱光を受光する（ステップＳ３）。このとき、受光部８では、照射開始からの単位時間毎（例えば、１ピコ秒毎の時刻ｔ１〜ｔｍ）の受光強度を内部メモリに記録しておく。
【００５９】
この受光部８が受光を完了したことを光強度取得部９に知らせると、この光強度取得部９では、真皮層２３から放射される後方散乱光の異なる時刻の受光強度を取得する（ステップＳ４）。すなわち、複数の時刻ｔ１〜ｔｍ各々における後方散乱光の光強度を取得する。
ここで、光強度取得部９が光強度を取得する時刻ｔ１〜ｔｍは、真皮層２３から放射される後方散乱光のピークとなる時刻を含むことが好ましい。すなわち、照射部５が特定波長λ１の短時間パルス光を照射した時刻に、真皮層２３の光路長が極大となる時間を加算した時刻とすることが好ましい。
【００６０】
次いで、光吸収係数算出部１２では、光強度取得部９にて取得した真皮層２３から放射される後方散乱光の異なる時刻の受光強度、すなわち、複数の時刻ｔ１〜ｔｍ各々における後方散乱光の光強度を基に、真皮層２３の光吸収係数を、下記の式（８）
【数８】
0005674093.tif 000009
（但し、Ｉ（ｔ）は受光部５が時刻ｔにて受光した光強度、Ｎ（ｔ）は特定波長λｋの短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度、Ｌｉ（ｔ）は皮膚の各々の層における伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長、μｉは第ｉ層の光吸収係数である）
から算出する（ステップＳ５）。
ここでは、第１層は表皮層、第２層は真皮層、第３層は皮下組織を示し、μ１は表皮層の光吸収係数、μ２は真皮層の光吸収係数、μ３は皮下組織の光吸収係数を示す。
【００６１】
次いで、特定波長λｋの種類数と同じ数の特定波長λｋに対して光吸収係数を算出したか否かを判定し、算出したならば、次の手順に進み、算出していなかったならば、特定波長λｋの短時間パルス光を照射する（ステップＳ１）以降の手順を再度行う（ステップＳ６）。
【００６２】
ここで、特定波長λｋの種類数と同じ数の特定波長λｋに対して光吸収係数を算出したならば、濃度算出部１３では、光吸収係数算出部１２が算出した真皮層２３の光吸収係数μを基に、真皮層２３に含まれるグルコースの濃度を算出する。
ここでは、例えば、単一波長での主要成分の光吸収係数と、各層の光吸収係数との関係から、グルコースの濃度を算出する場合について説明する。
下記の式（９）は、皮膚の光吸収係数が、水、グルコース、蛋白質、脂質及びその他の成分各々の波長に関する関数とその係数との積の和であることを示している。
【数９】
0005674093.tif 000010
【００６３】
式（１０）は、式（９）で示した皮膚の光吸収係数より、皮膚の表皮、真皮、皮下組織の各層の光吸収係数は、水、グルコース、蛋白質、脂質及びその他の成分各々の波長に関する光吸収係数と濃度との積の和であることを示している。
【数１０】
0005674093.tif 000011
【００６４】
式（１１）は、上記の式（１０）を行列で表した式であり、式（１１）中、「Ｃ」は、皮膚の表皮、真皮、皮下組織の各層毎の水、グルコース、蛋白質、脂質及びその他の成分各々の濃度を示す係数行列である。
【数１１】
0005674093.tif 000012
【００６５】
この式（１１）を変形することにより、「Ｃ」を解く連立一次方程式である式（１２）が導き出せる。
【数１２】
0005674093.tif 000013
【００６６】
この式（１２）を行列で表した式が、式（１３）である。
【数１３】
0005674093.tif 000014
【００６７】
ここでは、真皮（Ｌ２）の光吸収係数について、複数の波長（λ）における皮膚の主成分吸収係数とその濃度の関係について、式（１４）で示される一次連立方程式を立てる。
【数１４】
0005674093.tif 000015
【００６８】
この式（１４）は、式（１５）と変形することができる。
【数１５】
0005674093.tif 000016
【００６９】
この式（１５）を変形することにより、「Ｃ」を解く連立一次方程式である式（１６）を導くことができる。
式（１６）中、εＷは水の光吸収係数、εＧはグルコースの光吸収係数、εＰは蛋白質の光吸収係数、εＬは脂質の光吸収係数、μＬ１は表皮の光吸収係数、μＬ２は真皮の光吸収係数、μＬ３は皮下組織の光吸収係数である。なお、皮膚主成分の光吸収係数、皮膚の層はさらに追加されても良い。
【数１６】
0005674093.tif 000017
このように、真皮（Ｌ２）に係わる係数行列「Ｃ」は、式（１６）から算出することができる。
【００７０】
さらに、ここでは、複数波長での差分での説明を行う。
ここでは、真皮層２３に含まれるグルコースの濃度を下記の式（１７）
【数１７】
0005674093.tif 000018
（但し、μａは皮膚の任意の層である第ａ層における光吸収係数、ｇｊは皮膚を構成する第ｊ成分のモル濃度、εｊは第ｊ成分の光吸収係数、ｐは皮膚を構成する主成分の個数、ｑは短時間パルス光の種類数である）
から算出する（ステップＳ７）。
【００７１】
上述の血糖値測定装置１は、コンピュータシステムを内蔵しており、上述した各ステップの処理動作は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されている。そこで、このプログラムをコンピュータが読み出して実行することにより、上記の処理動作を行うことができる。
ここで、コンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等が挙げられる。
また、このコンピュータプログラムを通信回線によりコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしてもよい。
【００７２】
また、上記プログラムは、上記の各ステップの一部を実現するためのものであってもよい。
さらに、上述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であってもよい。
【００７３】
以上説明したように、本実施形態によれば、皮膚に照射される短時間パルス光を、水の吸収係数の温度変化が小さい特定波長の短時間パルス光とすることで、真皮層から放射される後方散乱光における水の影響を小さくすることができ、この後方散乱光を基に算出される皮膚の真皮層におけるグルコースの濃度においても、水の影響を小さくすることができる。したがって、グルコースの濃度における水の影響を低減することができ、真皮層に含まれるグルコースの濃度を、非侵襲的に精度良く測定することができる。
【００７４】
［第２の実施形態］
図１２は、本発明の第２の実施形態の血糖値測定装置の構成を示す概略ブロック図であり、本実施形態の血糖値測定装置３１が第１の実施形態の血糖値測定装置１と異なる点は、光強度取得部９及び光吸収係数算出部１２を、これらとは異なる機能を有する光強度取得部（光強度取得手段）３２及び光吸収係数算出部（光吸収係数算出手段）３３に替えた点である。
【００７５】
光強度取得部３２は、受光部８が受光した真皮層から放射される後方散乱光の所定の時刻から少なくとも所定の時刻τの間の光強度の時間変化を取得する。
光吸収係数算出部３３は、特定波長λｋの短時間パルス光を照射した皮膚の真皮層における光吸収係数を算出する。
【００７６】
この光吸収係数算出部３３では、皮膚における任意の層の光吸収係数を、下記の式（１８）
【数１８】
0005674093.tif 000019
（但し、Ｉ（ｔ）は受光部５が時刻ｔにて受光した光強度、Ｎ（ｔ）は特定波長λｋの短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度、Ｌｉ（ｔ）は皮膚の各々の層における伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長、ｎは皮膚の観測対象となる層の数、μｉは第ｉ層の光吸収係数である）
から算出する。
ここで、第１層は表皮層、第２層は真皮層、第３層は皮下組織を示し、μ１は表皮層の光吸収係数、μ２は真皮層の光吸収係数、μ３は皮下組織の光吸収係数を示す。
【００７７】
次に、この血糖値測定装置３１を用いて血糖値を測定する手順について、図１３に基づき説明する。
この手順では、光散乱媒質層選択部７が皮下組織２２、真皮層２３及び表皮層２４各々から放射される後方散乱光から、真皮層２３により放射される後方散乱光を選択する（ステップＳ２）手順までが図１１に示す手順と同一であるから、説明を省略する。
【００７８】
この後方散乱光を選択した後、受光部８により、真皮層２３から放射される所定の時刻τの間の後方散乱光を受光する（ステップＳ１１）。このとき、受光部８では、照射開始から少なくとも所定の時刻τの間の受光強度の時間変化を内部メモリに記録しておく。
次いで、この受光部８が受光を完了したことを光強度取得部３２に知らせると、この光強度取得部３２では、真皮層２３から放射される後方散乱光の照射開始から少なくとも所定の時刻τの間の受光強度の時間変化を取得する（ステップＳ１２）。
【００７９】
次いで、光吸収係数算出部３３では、光強度取得部３２にて取得した真皮層２３から放射される後方散乱光の照射開始から少なくとも所定の時刻τの間の受光強度の時間変化を基に、真皮層２３の光吸収係数を、下記の式（１９）
【数１９】
0005674093.tif 000020
（但し、Ｉ（ｔ）は受光部５が時刻ｔにて受光した光強度、Ｎ（ｔ）は特定波長λｋの短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度、Ｌｉ（ｔ）は皮膚の各々の層における伝搬光路長分布のモデルの時刻ｔにおける第ｉ層の光路長、ｎは皮膚の観測対象となる層の数、μｉは第ｉ層の光吸収係数である）
から算出する（ステップＳ１３）。
【００８０】
次いで、特定波長λｋの種類数と同じ数の特定波長λｋに対して光吸収係数を算出したか否かを判定し、算出したならば、次の手順に進み、算出していなかったならば、特定波長λｋの短時間パルス光を照射する（ステップＳ１）以降の手順を再度行う（ステップＳ１４）。
【００８１】
ここで、特定波長λｋの種類数と同じ数の特定波長λｋに対して光吸収係数を算出したならば、濃度算出部１３では、光吸収係数算出部３３が算出した真皮層２３の光吸収係数μを基に、真皮層２３に含まれるグルコースの濃度を、下記の式（２０）
【数２０】
0005674093.tif 000021
（但し、μａは皮膚の任意の層である第ａ層における光吸収係数、ｇｊは皮膚を構成する第ｊ成分のモル濃度、εｊは第ｊ成分の光吸収係数、ｐは皮膚を構成する主成分の個数、ｑは特定波長λｋの種類数である）
から算出する（ステップＳ１５）。
【００８２】
本実施形態においても、第１の実施形態と同様に、真皮層から放射される後方散乱光における水の影響を小さくすることができ、この後方散乱光を基に算出される皮膚の真皮層におけるグルコースの濃度においても、水の影響を小さくすることができる。したがって、グルコースの濃度における水の影響を低減することができ、真皮層に含まれるグルコースの濃度を、非侵襲的に精度良く測定することができる。
【００８３】
以上、本発明の各実施形態について、図面を参照して説明してきたが、具体的な構成は上述のものに限られることはなく、本発明の要旨を逸脱しない範囲内において様々な設計変更等が可能である。
例えば、上記の各実施形態では、濃度定量装置として血糖値測定装置を、観測対象として人の手のひらの皮膚を、目的成分としてグルコースを、特定波長の光として特定波長の短時間パルス光を、それぞれ取ることで、皮膚の真皮層に含まれるグルコースの濃度を測定する場合について説明したが、これに限らず、濃度定量方法を、複数の光散乱媒質の層から形成される観測対象の任意の層における目的成分の濃度を定量する他の装置に用いてもよく、特定波長の短時間パルス光を、特定波長の連続光に替えてもよい。
例えば、携帯型の皮膚主成分の濃度測定装置に適用した場合、皮膚疾患の検査や診断や治療に有効利用することが可能である。
【符号の説明】
【００８４】
１…血糖値測定装置（濃度定量装置）、３…光路長分布記憶部（光路長分布記憶手段）、４…時間分解波形記憶部（時間分解波形記憶手段）、５…照射部（照射手段）、７…光散乱媒質層選択部（光散乱媒質層選択手段）、８…受光部（受光手段）、１０…光路長取得部（光路長取得手段）、１１…無吸収時光強度取得部（光強度モデル取得手段）、１２…光吸収係数算出部（光吸収係数算出手段）、１３…濃度算出部（濃度算出手段）、２１…皮膚（観測対象）、２３…真皮層（任意の層）、３１…血糖値測定装置（濃度定量装置）、３２…光強度取得部（光強度取得手段）、３３…光吸収係数算出部（光吸収係数算出手段）、Ｓ１〜Ｓ６、Ｓ１１〜Ｓ１４  ステップ
【図１】
0005674093.tif 000029
【図２】
0005674093.tif 000030
【図３】
0005674093.tif 000031
【図４】
0005674093.tif 000032
【図５】
0005674093.tif 000033
【図６】
0005674093.tif 000034
【図７】
0005674093.tif 000035
【図８】
0005674093.tif 000036
【図９】
0005674093.tif 000037
【図１０】
0005674093.tif 000038
【図１１】
0005674093.tif 000039
【図１２】
0005674093.tif 000040
【図１３】
0005674093.tif 000041