(57)【特許請求の範囲】
【請求項１】
複数の層により構成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置であって、
前記観測対象に短時間パルス光を照射する照射部と、
前記短時間パルス光の照射により前記観測対象から後方散乱される光を受光する受光部と、
前記観測対象に対して照射する短時間パルス光の、前記複数の層の各々の層における光路長分布のモデルを記憶する光路長分布記憶部と、
前記光路長分布のモデルの所定の時刻における、前記複数の層の各々の層の光路長を取得する光路長取得部と、
前記観測対象に対して照射する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶部と、
時間分解波形のモデルの前記所定の時刻における光強度を取得する光強度モデル取得部と、
前記受光部が受光した前記光の強度を取得する光強度取得部と、
前記光強度取得部が取得した前記光の強度の光強度分布と、前記光路長取得部が取得した前記複数の層の各々の層の光路長と、前記光強度モデル取得部が取得した光強度モデルとに基づいて、前記光強度分布から前記任意の層の光強度分布に対応する領域の時間の範囲である積分区間を算出する積分区間算出部と、
前記積分区間算出部が算出した前記積分区間を変化させて前記任意の層における目的成分の光吸収係数を算出する光吸収係数算出部と、
前記光吸収係数算出部が取得した前記目的成分の光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出部と、
を備えてなることを特徴とする濃度定量装置。
【請求項２】
前記光吸収係数算出部は、
前記積分区間算出部が算出した積分区間に基づいて前記任意の層の光吸収係数を算出する光吸収係数算出部と、
前記光吸収係数算出部が算出した前記任意の層の光吸収係数に基づき前記積分区間における光吸収係数分布を算出し記憶する光吸収係数分布記憶部と、
前記光吸収係数分布記憶部が取得した光吸収係数分布に基づき、前記任意の層における前記目的成分の光吸収係数を取得する光吸収係数取得部と、
を備えてなることを特徴とする請求項１記載の濃度定量装置。
【請求項３】
前記光強度取得部は、前記積分区間内の前記観測対象の層の数ｎ以上となる複数の時刻ｔ１〜ｔｍにおける光強度を取得し（但し、ｎは１以上の自然数、ｍはｎ以上の自然数）、
前記光吸収係数算出部は、
自然対数を示すｌｎ（・）、前記受光部が前記積分区間内の時刻ｔにおいて受光した光強度を示すＩ（ｔ）、前記短時間パルス光の時間分解波形のモデルの前記積分区間内の時刻ｔにおける光強度を示すＮ（ｔ）、前記光路長分布のモデルの前記積分区間内の時刻ｔにおける第ｉ層の光路長を示すＬｉ（ｔ）、第ｉ層の光吸収係数を示すμｉを用いて、
【数１】
0005521199.tif 000017
から任意の層の光吸収係数を算出する、
ことを特徴とする請求項２記載の濃度定量装置。
【請求項４】
前記光強度取得部が光強度を取得する複数の時刻は、前記複数の層の各々の層の光路長分布のピーク時間を含むことを特徴とする請求項３記載の濃度定量装置。
【請求項５】
前記光強度取得部は、前記積分区間内の所定の時刻から少なくとも所定の時間τの間の光強度を取得し、
前記光吸収係数算出部は、
自然対数を示すｌｎ（・）、前記受光部が前記積分区間内の時刻ｔにおいて受光した光強度を示すＩ（ｔ）、前記短時間パルス光の時間分解波形のモデルの前記積分区間内の時刻ｔにおける光強度を示すＮ（ｔ）、前記光路長分布のモデルの前記積分区間内の時刻ｔにおける第ｉ層の光路長を示すＬｉ（ｔ）、前記観測対象の層の数を示すｎ、第ｉ層の光吸収係数を示すμｉを用いて、
【数２】
0005521199.tif 000018
から任意の層の光吸収係数を算出する、
ことを特徴とする請求項２記載の濃度定量装置。
【請求項６】
前記照射部は、複数の波長１〜ｑの光を照射し、
前記光吸収係数算出部は、前記任意の層における光吸収係数を前記照射部が照射した複数の波長毎に算出し、
前記濃度算出部は、
前記任意の層である第ａ層における波長ｉの光吸収係数を示すμａ（ｉ）、前記観測対象を形成する第ｊ成分のモル濃度を示すｇｊ、第ｊ成分の波長ｉに対する光吸収係数を示すεｊ（ｉ）、前記観測対象を形成する主成分の個数を示すｐ、照射部が照射する波長の種類数を示すｑを用いて、
【数３】
0005521199.tif 000019
から前記任意の層における前記目的成分の濃度を算出する、
ことを特徴とする請求項２ないし５のいずれか１項記載の濃度定量装置。
【請求項７】
前記照射部が照射する複数の光は、前記目的成分の光吸収係数が大きくなる波長の光を含むことを特徴とする請求項６記載の濃度定量装置。
【請求項８】
前記照射部が照射する複数の光は、前記観測対象を構成する主成分の各々の成分の吸収スペクトル分布の直交性が高くなる波長の光を含むことを特徴とする請求項６記載の濃度定量装置。
【請求項９】
複数の層により構成される観測対象に短時間パルス光を照射する照射部と、前記短時間パルス光の照射により前記観測対象から後方散乱される光を受光する受光部と、前記観測対象に対して照射する短時間パルス光の、前記複数の層の各々の層における光路長分布のモデルを記憶する光路長分布記憶部と、前記観測対象に対して照射する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶部とを備え、前記観測対象のうち任意の層における目的成分の濃度を定量する濃度定量装置を用いた濃度定量方法であって、
前記照射部により、前記観測対象に短時間パルス光を照射し、
前記受光部により、前記短時間パルス光の照射により前記観測対象から後方散乱される光を受光し、
光強度取得部により、前記照射部が短時間パルス光を照射した時刻以降の所定の時刻にて前記受光部が受光した光の強度を取得し、
光路長取得部により、前記光路長分布記憶部から、前記光路長分布のモデルの前記所定の時刻における前記複数の層の各々の層の光路長を取得し、
光強度モデル取得部により、前記時間分解波形記憶部から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光強度を取得し、
積分区間算出部により、前記光強度取得部が取得した前記光強度分布と、前記光路長分布記憶部が記憶した前記複数の層の各々の層の光路長と、前記時間分解波形記憶部が記憶した前記光強度とに基づいて、前記光強度分布から前記任意の層の光強度分布に対応する領域の時間の範囲である積分区間を算出し、
光吸収係数算出部により、前記積分区間算出部が算出した積分区間を変化させて前記任意の層における目的成分の光吸収係数を算出し、
濃度算出部により、前記光吸収係数算出部が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する、
ことを特徴とする濃度定量方法。
【請求項１０】
複数の層により構成される観測対象に短時間パルス光を照射する照射部と、前記短時間パルス光の照射により前記観測対象から後方散乱される光を受光する受光部と、前記観測対象に対して照射する短時間パルス光の、前記複数の層の各々の層における光路長分布のモデルを記憶する光路長分布記憶部と、前記観測対象に対して照射する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶部とを備え、前記観測対象のうち任意の層における目的成分の濃度を定量する濃度定量装置のコンピューターに、
前記観測対象に前記短時間パルス光を照射する照射手順、
前記短時間パルス光の照射により前記観測対象から後方散乱される光を受光する受光手順、
前記短時間パルス光を照射した時刻以降の所定の時刻にて受光した光の強度を取得する光強度取得手順、
前記光路長分布記憶部から、前記光路長分布のモデルの前記所定の時刻における前記複数の層の各々の層の光路長を取得する光路長取得手順、
前記時間分解波形記憶部から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光強度を取得する光強度モデル取得手順、
前記光強度取得手順により取得された前記光強度分布と、前記光路長取得手順により取得された前記複数の層の各々の層の光路長と、前記光強度モデル取得手順により取得された前記光強度とに基づいて、前記光強度分布から前記任意の層の光強度分布に対応する領域の時間の範囲である積分区間を算出する積分区間算出手順、
前記積分区間算出手順により算出された積分区間を変化させて前記任意の層における目的成分の光吸収係数を算出する光吸収係数算出手順、
前記光吸収係数算出手順により算出された光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手順、
を実行させることを特徴とするプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の光散乱媒質の層により構成される観測対象のうち、任意の層における目的成分の濃度を、非侵襲的にかつ精度良く定量する濃度定量装置及び濃度定量方法並びにプログラムに関するものである。
【背景技術】
【０００２】
近年、我が国は飽食の時代にあって、糖尿病の患者が毎年増加し続けている。そのために、糖尿病性腎炎の患者も毎年増加し続けることとなり、その結果、慢性腎不全の患者も毎年１万人もの増加を続け、患者数は２８万人を超えるようになってきている。
一方、高齢化社会の到来により、予防医学に対する要求の高まりを受けて、個人における代謝量管理の重要性が急速に増大している。中でも、血糖値測定は、食前や食後の血糖値を測定することで糖代謝の反応が分かることが知られており、糖尿病のごく初期段階での糖代謝の反応を評価することで、糖尿病の早期診断に基づく早期治療が可能になる。
【０００３】
従来、血糖値の測定は、腕あるいは指先等の静脈から採血を行い、この血液中のグルコースに対する酵素活性を測定することで行っている。しかし、このような血糖値の測定方法では、採血が煩雑であり、しかも採血に痛みを伴い、さらには感染症の危険性を伴う等の様々な問題がある。
また、血糖値を連続的に測定する方法としては、静脈に注射針を刺した状態で連続的に血糖値相応のグルコースの定量を行う機器が米国にて開発されており、現在臨床試験中である。しかし、静脈に注射針を刺したままにしているために、血糖値の測定中に針が抜ける危険性や感染症の危険性がある。
そこで、採血無しに頻繁に血糖値を測定することができ、しかも感染症の危険性が無い血糖値の測定装置の開発が求められている。さらには、簡単にかつ常時装着可能であり、小型化可能な血糖値の測定装置の開発が求められている。
【０００４】
近赤外の連続光を用いて非侵襲的に血糖値を測定する装置としては、分子吸光の原理を用いた一般的な分光分析測定の原理を適用した装置が提案されている（例えば、特許文献１、２参照）。
この装置は、皮膚の赤外スペクトルを用いて生体成分濃度の定量をおこなう場合に、皮下脂肪の影響を受けて生体成分濃度の定量に誤差が生じることに対応したもので、より具体的には、皮膚に近赤外の連続光を照射し、その光吸収量からグルコースの濃度を算出する装置である。
この装置では、予めグルコース濃度と照射する近赤外光の波長と光の吸収量との関係を示す検量線を作成しておき、皮膚に近赤外の連続光を照射し、この皮膚からの戻り光をモノクロメーター等を用いてある波長域を走査し、その波長域の各波長に対する光の吸収量を求め、この各波長における光の吸収量と検量線とを比較することで、血液中のグルコース濃度、すなわち血糖値を算出している。
【０００５】
また、１７００ｎｍ〜１８００ｎｍの波長範囲から選択した皮下脂肪の特異吸収波長での吸光度から、皮膚の性状の分類を行い、「皮膚厚さ」の代用特性として検量式を選択している。
さらには、予備的に近赤外の受光部と発光部との間隔を６５０μｍとして推定した「皮膚厚さ」を１．２ｍｍ以上、１．２ｍｍ未満のいずれかに判断し、受光部と発光部との間隔を６５０μｍ、３００μｍのいずれかに選択した後に検量式を選択している。
【０００６】
一方、近赤外光を用いた生体診断としては、例えば、時間分解計測法を用いた生体組織イメージングにより皮膚主成分における近赤外光の吸収量を測定し、この吸収量を基に皮膚主成分の各割合、例えば、血糖相応のグルコース濃度を求める方法が知られている。
この皮膚主成分の吸収量には波長依存性があるので、通常、予め皮膚主成分の定量に影響を及ぼす変動要因を多変量解析で複数の割合で変化させた複数のスペクトラムを作製しておき、皮膚主成分における近赤外光の吸収量の測定結果のスペクトルを上記の複数のスペクトラムと比較し、これらのスペクトラムから一致するスペクトラムを選ぶことにより、皮膚主成分の各割合を推定する方法が採られている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特許第３９３１６３８号公報
【特許文献２】特許第３９９４５８８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、従来の近赤外の連続光を用いた非侵襲的に血糖値を測定する装置では、特定深さを通過する経路の光の吸収量のみを測定することができず、したがって、特定深さの皮膚主成分における血糖相応のグルコース濃度を精度よく定量することができないという問題点があった。
また、特許文献１の装置では、皮膚表面から皮下脂肪までの深さを「皮膚厚さ」として、皮下脂肪の特異吸収波長での吸光度から皮膚の性状を分類すること、例えば、皮膚表面から皮下脂肪までの深さを「皮膚厚さ」として代用することには、（１）皮膚の真皮と皮下組織の境界は、皮膚の表面からの深さとして均一では無いこと、（２）真皮には脂肪を分泌する汗腺があって脂肪分泌物を蓄えていること、（３）皮下脂肪の特異吸収波長での吸光度から皮膚の性状の分類を行う場合、真皮の細胞及び間質液には脂肪が含まれているので、真皮と皮下脂肪との区別が難しい、等の理由により問題点があった。
【０００９】
一般に、皮膚の赤外スペクトルを用いて生体成分濃度の定量を行う場合、受光部と発光部との間隔によって定まるバナナシェイプ特性により、皮膚内での光路の皮膚表面からの深さが概ね推定される。例えば、受光部と発光部との間隔を６５０μｍとすれば、光路の皮膚表面からの深さは３２５μｍと推定され、また、受光部と発光部との間隔を３００μｍとすれば、光路の皮膚表面からの深さは１５０μｍと推定される。
しかしながら、特許文献１の装置では、上記の理由等により、皮膚の赤外スペクトルを用いて生体成分濃度の定量を行う部位を特定することができず、したがって、真皮中で間質成分の一つとしてグルコースが存在している網状層（Stratum reticulare）を特定部位として、この特定部位を透過する光路での吸光度を選択的に測定することはできない。
【００１０】
本発明は、上記の課題を解決するためになされたものであって、複数の層により構成される観測対象のうち、定量を行う特定部位である任意の層からの後方散乱光を他の層からの光と区別することで、任意の層における目的成分の濃度を、非侵襲的にかつ精度良く定量することができる濃度定量装置及び濃度定量方法並びにプログラムを提供することを目的とする。
【課題を解決するための手段】
【００１１】
上記の課題を解決するために、本発明は以下の濃度定量装置及び濃度定量方法並びにプログラムを採用した。
すなわち、本発明の濃度定量装置は、複数の層により構成される観測対象のうち、任意の層における目的成分の濃度を定量する濃度定量装置であって、前記観測対象に短時間パルス光を照射する照射部と、前記短時間パルス光の照射により前記観測対象から後方散乱される光を受光する受光部と、前記観測対象に対して照射する短時間パルス光の、前記複数の層の各々の層における光路長分布のモデルを記憶する光路長分布記憶部と、前記光路長分布のモデルの所定の時刻における、前記複数の層の各々の層の光路長を取得する光路長取得部と、前記観測対象に対して照射する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶部と、時間分解波形のモデルの前記所定の時刻における光強度を取得する光強度モデル取得部と、前記受光部が受光した前記光の強度を取得する光強度取得部と、前記光強度取得部が取得した前記光の強度の光強度分布と、前記光路長取得部が取得した前記複数の層の各々の層の光路長と、前記光強度モデル取得部が取得した光強度モデルとに基づいて、前記光強度分布から前記任意の層の光強度分布に対応する領域の時間の範囲である積分区間を算出する積分区間算出部と、前記積分区間算出部が算出した前記積分区間を変化させて前記任意の層における目的成分の光吸収係数を算出する光吸収係数算出部と、前記光吸収係数算出部が取得した前記目的成分の光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出部と、を備えてなることを特徴とする。
【００１２】
本発明の濃度定量装置では、照射部により、前記観測対象に前記短時間パルス光を照射し、受光部により、前記短時間パルス光の照射により前記観測対象から後方散乱される光を受光し、光強度取得部により、前記受光部が受光した光の強度を取得し、積分区間算出部により、前記光強度取得部が取得した前記光強度分布と、前記光路長分布記憶部が取得した前記複数の層の各々の層の光路長と、前記時間分解波形記憶部が取得した前記光強度とに基づいて、前記光強度分布から前記任意の層の光強度に対応する領域の時間の範囲である積分区間を算出する。
このように、積分区間算出部により算出された積分区間を基に、受光部が受光した光強度から前記積分区間に対応する時間帯の光強度を取得することにより、定量を行う特定部位である任意の層からの光を他の層からの光と区別して測定することができ、任意の層からの光に対する他の層からの光の影響を低減することができる。したがって、任意の層における目的成分の光の吸収量、すなわち目的成分の濃度を精度良く測定することができ、その結果、任意の層における目的成分の濃度を、非侵襲的にかつ精度良く定量することができる。
【００１３】
本発明の濃度定量装置は、前記光吸収係数算出部は、前記積分区間算出部が算出した積分区間に基づいて前記任意の層の光吸収係数を算出する光吸収係数算出部と、前記光吸収係数算出部が算出した前記任意の層の光吸収係数に基づき前記積分区間における光吸収係数分布を算出し記憶する光吸収係数分布記憶部と、前記光吸収係数分布記憶部が取得した光吸収係数分布に基づき、前記任意の層における前記目的成分の光吸収係数を取得する光吸収係数取得部と、を備えてなることを特徴とする。
【００１４】
本発明の濃度定量装置では、光吸収係数算出部により、積分区間算出部が算出した積分区間に基づいて任意の層の光吸収係数を算出し、光吸収係数分布記憶部により、光吸収係数算出部が算出した任意の層の光吸収係数に基づき積分区間における光吸収係数分布を算出し記憶する。
このように、積分区間に基づいて算出した任意の層の光吸収係数に基づき、積分区間における光吸収係数分布を算出することで、任意の層における目的成分の光吸収係数、すなわち目的成分の濃度をさらに精度良く測定することができる。
【００１５】
本発明の濃度定量装置は、前記光強度取得部は、前記積分区間内の前記観測対象の層の数ｎ以上となる複数の時刻ｔ１〜ｔｍにおける光強度を取得し（但し、ｎは１以上の自然数、ｍはｎ以上の自然数）、前記光吸収係数算出部は、自然対数を示すｌｎ（・）、前記受光部が前記積分区間内の時刻ｔにおいて受光した光強度を示すＩ（ｔ）、前記短時間パルス光の時間分解波形のモデルの前記積分区間内の時刻ｔにおける光強度を示すＮ（ｔ）、前記光路長分布のモデルの前記積分区間内の時刻ｔにおける第ｉ層の光路長を示すＬｉ（ｔ）、第ｉ層の光吸収係数を示すμｉを用いて、
【数１】
0005521199.tif 000002
から任意の層の光吸収係数を算出する、ことを特徴とする。
【００１６】
本発明の濃度定量装置では、光強度取得部が、積分区間内の観測対象の層の数ｎ以上となる複数の時刻ｔ１〜ｔｍにおける光強度を取得し、光吸収係数算出部が、任意の層の光吸収係数を、上記の式（１）から算出する。
このように、後方散乱光を時間分解計測することで、任意の層以外の層からの後方散乱光をノイズとして低減することができ、目的成分の濃度における任意の層以外の層からの影響を低減することができる。したがって、目的成分の濃度をさらに精度良く測定することができる。
【００１７】
本発明の濃度定量装置は、前記光強度取得部が光強度を取得する複数の時刻は、前記複数の層の各々の層の光路長分布のピーク時間を含むことを特徴とする。
【００１８】
本発明の濃度定量装置では、光強度取得部が光強度を取得する複数の時刻が、複数の層の各々の層の光路長分布のピーク時間を含むことにより、観測対象中の複数の層から任意の層を効率的に選択することができる。したがって、任意の層における目的成分の濃度をさらに精度良く測定することができる。
【００１９】
本発明の濃度定量装置は、前記光強度取得部は、前記積分区間内の所定の時刻から少なくとも所定の時間τの間の光強度を取得し、前記光吸収係数算出部は、自然対数を示すｌｎ（・）、前記受光部が前記積分区間内の時刻ｔにおいて受光した光強度を示すＩ（ｔ）、前記短時間パルス光の時間分解波形のモデルの前記積分区間内の時刻ｔにおける光強度を示すＮ（ｔ）、前記光路長分布のモデルの前記積分区間内の時刻ｔにおける第ｉ層の光路長を示すＬｉ（ｔ）、前記観測対象の層の数を示すｎ、第ｉ層の光吸収係数を示すμｉを用いて、
【数２】
0005521199.tif 000003
から任意の層の光吸収係数を算出する、ことを特徴とする。
【００２０】
本発明の濃度定量装置では、光強度取得部が、積分区間内の所定の時刻から少なくとも所定の時刻τの間の光強度を取得し、光吸収係数算出部が、任意の層の光吸収係数を、上記の式（２）から算出する。
このように、後方散乱光を時間分解計測することで、任意の層以外の層からの後方散乱光をノイズとして低減することができ、目的成分の濃度における任意の層以外の層からの影響を低減することができる。したがって、目的成分の濃度をさらに精度良く測定することができる。
【００２１】
本発明の濃度定量装置は、前記照射部は、複数の波長１〜ｑの光を照射し、前記光吸収係数算出部は、前記任意の層における光吸収係数を前記照射部が照射した複数の波長毎に算出し、前記濃度算出部は、前記任意の層である第ａ層における波長ｉの光吸収係数を示すμａ（ｉ）、前記観測対象を形成する第ｊ成分のモル濃度を示すｇｊ、第ｊ成分の波長ｉに対する光吸収係数を示すεｊ（ｉ）、前記観測対象を形成する主成分の個数を示すｐ、照射部が照射する波長の種類数を示すｑを用いて、
【数３】
0005521199.tif 000004
から前記任意の層における前記目的成分の濃度を算出する、ことを特徴とする。
【００２２】
本発明の濃度定量装置では、照射部が、複数の波長１〜ｑの光を照射し、光吸収係数算出部が、任意の層における光吸収係数を照射部が照射した複数の波長毎に算出し、濃度算出部が、任意の層における目的成分の濃度を上記の式（３）から算出する。
このように、後方散乱光を時間分解計測することで、任意の層以外の層からの後方散乱光をノイズとして低減することができ、目的成分の濃度における任意の層以外の層からの影響を低減することができる。したがって、目的成分の濃度をさらに精度良く測定することができる。
【００２３】
本発明の濃度定量装置は、前記照射部が照射する複数の光は、前記目的成分の光吸収係数が大きくなる波長の光を含むことを特徴とする。
【００２４】
本発明の濃度定量装置は、前記照射部が照射する複数の光は、前記観測対象を構成する主成分の各々の成分の吸収スペクトル分布の直交性が高くなる波長の光を含むことを特徴とする。
【００２５】
本発明の濃度定量方法は、複数の層により構成される観測対象に短時間パルス光を照射する照射部と、前記短時間パルス光の照射により前記観測対象から後方散乱される光を受光する受光部と、前記観測対象に対して照射する短時間パルス光の、前記複数の層の各々の層における光路長分布のモデルを記憶する光路長分布記憶部と、前記観測対象に対して照射する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶部とを備え、前記観測対象のうち任意の層における目的成分の濃度を定量する濃度定量装置を用いた濃度定量方法であって、前記照射部により、前記観測対象に短時間パルス光を照射し、前記受光部により、前記短時間パルス光の照射により前記観測対象から後方散乱される光を受光し、光強度取得部により、前記照射部が短時間パルス光を照射した時刻以降の所定の時刻にて前記受光部が受光した光の強度を取得し、光路長取得部により、前記光路長分布記憶部から、前記光路長分布のモデルの前記所定の時刻における前記複数の層の各々の層の光路長を取得し、光強度モデル取得部により、前記時間分解波形記憶部から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光強度を取得し、積分区間算出部により、前記光強度取得部が取得した前記光強度分布と、前記光路長分布取得部が取得した前記複数の層の各々の層の光路長と、前記時間分解波形取得部が取得した前記光強度とに基づいて、前記光強度分布から前記任意の層の光強度分布に対応する領域の時間の範囲である積分区間を算出し、光吸収係数算出部により、前記積分区間算出部が算出した積分区間を変化させて前記任意の層における目的成分の光吸収係数を算出し、濃度算出部により、前記光吸収係数算出部が算出した光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する、ことを特徴とする。
【００２６】
本発明の濃度定量方法では、照射部により、観測対象に短時間パルス光を照射し、受光部により、短時間パルス光の照射により観測対象から後方散乱される光を受光し、光強度取得部により、照射部が短時間パルス光を照射した時刻以降の所定の時刻にて受光部が受光した光の強度を取得し、積分区間算出部により、前記光強度取得部が取得した前記光強度分布と、前記光路長分布取得部が取得した前記複数の層の各々の層の光路長と、前記時間分解波形取得部が取得した前記光強度とに基づいて、前記光強度分布から前記任意の層の光強度に対応する領域の積分区間を算出する。
【００２７】
このように、積分区間算出部により、光強度分布から任意の層の光強度に対応する領域の積分区間を算出することで、この積分区間算出部により算出された積分区間を基に、受光部が受光した光強度から前記積分区間に対応する時間帯の光強度を取得することができ、任意の層からの光を他の層からの光と区別して測定することができ、任意の層からの光に対する他の層からの光の影響を低減することができる。したがって、任意の層における目的成分の光の吸収量、すなわち目的成分の濃度を精度良く測定することができ、その結果、任意の層における目的成分の濃度を、非侵襲的に精度良くかつ効率よく定量することができる。
【００２８】
本発明のプログラムは、複数の層により構成される観測対象に短時間パルス光を照射する照射部と、前記短時間パルス光の照射により前記観測対象から後方散乱される光を受光する受光部と、前記観測対象に対して照射する短時間パルス光の、前記複数の層の各々の層における光路長分布のモデルを記憶する光路長分布記憶部と、前記観測対象に対して照射する短時間パルス光の時間分解波形のモデルを記憶する時間分解波形記憶部とを備え、前記観測対象のうち任意の層における目的成分の濃度を定量する濃度定量装置のコンピューターに、前記観測対象に前記短時間パルス光を照射する照射手順、前記短時間パルス光の照射により前記観測対象から後方散乱される光を受光する受光手順、前記短時間パルス光を照射した時刻以降の所定の時刻にて受光した光の強度を取得する光強度取得手順、前記光路長分布記憶部から、前記光路長分布のモデルの前記所定の時刻における前記複数の層の各々の層の光路長を取得する光路長取得手順、前記時間分解波形記憶部から、前記短時間パルス光の時間分解波形のモデルの前記所定の時刻における光強度を取得する光強度モデル取得手順、前記光強度取得手順により取得された前記光強度分布と、前記光路長取得手順により取得された前記複数の層の各々の層の光路長と、前記光強度モデル取得手順により取得された前記光強度とに基づいて、前記光強度分布から前記任意の層の光強度分布に対応する領域の時間の範囲である積分区間を算出する積分区間算出手順、前記積分区間算出手順により算出された積分区間を変化させて前記任意の層における目的成分の光吸収係数を算出する光吸収係数算出手順、前記光吸収係数算出手順により算出された光吸収係数に基づいて、前記任意の層における前記目的成分の濃度を算出する濃度算出手順、を実行させることを特徴とする。
【００２９】
本発明のプログラムでは、濃度定量装置のコンピューターに、観測対象に短時間パルス光を照射する照射手順、短時間パルス光の照射により観測対象から後方散乱される光を受光する受光手順、短時間パルス光を照射した時刻以降の所定の時刻にて受光した光の強度を取得する光強度取得手順、を実行させた後、光強度取得部により取得された光強度分布と、光路長分布記憶部により取得された複数の層の各々の層の光路長と、時間分解波形記憶部により取得された光強度とに基づいて、光強度分布から任意の層の光強度に対応する領域の積分区間を算出する積分区間算出手順、積分区間算出手順により算出された積分区間に基づいて任意の層における目的成分の光吸収係数を算出する光吸収係数算出手順、光吸収係数算出手順により算出された光吸収係数に基づいて、任意の層における目的成分の濃度を算出する濃度算出手順、を順次実行させる。
【００３０】
このように、光強度取得手順により取得された光強度と、光路長取得手順により取得された複数の層の各々の層の光路長と、光強度モデル取得手順により取得された光強度とに基づいて、光強度分布から任意の層の光強度に対応する領域の積分区間を算出する積分区間算出手順を実行することで、この積分区間算出部により算出された積分区間を基に、受光部が受光した光の強度から前記積分区間に対応する時間帯の光の強度のみを取得することができ、任意の層からの光に対する他の層からの光の影響を低減することができる。したがって、任意の層における目的成分の光の吸収量、すなわち目的成分の濃度を精度良く測定することができ、その結果、任意の層における目的成分の濃度を、非侵襲的に精度良くかつ効率よく定量することができる。
【図面の簡単な説明】
【００３１】
【図１】本発明の第１の実施形態の血糖値測定装置の構成を示す概略ブロック図である。
【図２】人の皮膚組織の断面を示す模式図である。
【図３】シミュレーション部が算出した各層の光路長分布を示す図である。
【図４】シミュレーション部が算出した時間分解波形を示す図である。
【図５】皮膚の主成分の吸収スペクトルを示す図である。
【図６】本発明の第１の実施形態の血糖値測定装置が血糖値を測定する動作を示すフローチャートである。
【図７】本発明の第２の実施形態の血糖値測定装置が血糖値を測定する動作を示すフローチャートである。
【図８】本発明の第３の実施形態の血糖値測定装置の構成を示す概略ブロック図である。
【図９】本発明の第４の実施形態の血糖値測定装置の構成を示す概略ブロック図である。
【図１０】シミュレーション部で得た無吸収時光強度Ｎ（ｔ）の時間分解波形を示す図である。
【図１１】シミュレーション部で得た吸収時光強度Ｉ（ｔ）の時間分解波形をリニアスケールで示した図である。
【図１２】シミュレーション部で得た吸収時光強度Ｉ（ｔ）の時間分解波形をログスケールで示した図である。
【図１３】本発明の第４の実施形態の血糖値測定装置が血糖値を測定する動作を示すフローチャートである。
【図１４】本発明の第４の実施形態の血糖値測定装置が血糖値を測定する動作を示すフローチャートである。
【図１５】真皮層における光吸収係数と積分区間との関係を示す図である。
【図１６】本発明の第５の実施形態の血糖値測定装置が血糖値を測定する動作を示すフローチャートである。
【図１７】本発明の第５の実施形態の血糖値測定装置が血糖値を測定する動作を示すフローチャートである。
【図１８】本発明の第６の実施形態の血糖値測定装置の構成を示す概略ブロック図である。
【図１９】本発明の第７の実施形態の血糖値測定装置の構成を示す概略ブロック図である。
【発明を実施するための形態】
【００３２】
本発明の濃度定量装置及び濃度定量方法並びにプログラムを実施するための形態について説明する。
本発明では、濃度定量装置として血糖値測定装置を、観測対象として人の手のひらの皮膚を、目的成分としてグルコースを、それぞれ例に取り説明する。
【００３３】
［第１の実施形態］
図１は、本発明の第１の実施形態の血糖値測定装置の構成を示す概略ブロック図である。
この血糖値測定装置１００は、シミュレーション部１０１、光路長分布記憶部１０２、時間分解波形記憶部１０３、照射部１０４、受光部１０５、計測光強度取得部１０６、光路長取得部１０７、無吸収時光強度取得部（光強度モデル取得部）１０８、光吸収係数算出部１０９、濃度算出部１１０を備えている。
この血糖値測定装置１００は、皮膚（観測対象）の真皮層（任意の層）に含まれるグルコース（目的成分）の濃度を測定する。
【００３４】
シミュレーション部１０１は、光吸収係数がゼロの皮膚モデルに対して光を照射するシミュレーションを行う。
光路長分布記憶部１０２は、光吸収係数がゼロの皮膚モデルの光路長分布を記憶する。
時間分解波形記憶部１０３は、光吸収係数がゼロの皮膚モデルの時間分解波形を記憶する。
照射部１０４は、皮膚に対して短時間パルス光を照射する。この照射部１０４が照射する複数の短時間パルス光は、皮膚を構成する主成分の各々の成分の吸収スペクトル分布の直交性が高くなる波長の光、すなわち、皮膚を構成する主成分の各々の成分のうち、ある主成分における特定成分の吸収スペクトルの極大値が他の成分の吸収スペクトルの極大値と大きく異なる波長の光を含んでいる。
【００３５】
受光部１０５は、短時間パルス光が皮膚によって後方散乱した光を受光する。
ここで、短時間パルス光とは、パルス幅の時間が照射部１０４から受光部１０５へ光が空気中を直接伝搬する時間よりも短いパルス光のことであり、例えば、パルス光の半値幅が０．１ｐｓ〜１０ｐｓ、２つのパルス光の間の時間間隔が１ｐｓ〜１００ｐｓのパルス光のことである。
また、光路長分布とは、光（光子）の移動経路の長さ（光路長）を当該光（光子）が受光部１０５に到達するまでの時間を基に分布関数として表したものである。
【００３６】
計測光強度取得部１０６は、受光部１０５が受光した光のある時刻における光強度を取得する。
光路長取得部１０７は、光路長分布記憶部１０２からある時刻における光路長を取得する。
無吸収時光強度取得部１０８は、時間分解波形記憶部１０３からある時刻における光の吸収係数をゼロ（零）としたときの光強度モデルを取得する。
光吸収係数算出部１０９は、短時間パルス光を照射した皮膚の真皮層における光吸収係数を算出する。
濃度算出部１１０は、真皮層に含まれるグルコースの濃度を算出する。
【００３７】
そして、血糖値測定装置１００において、照射部１０４は、皮膚に短時間パルス光を照射し、受光部１０５は、短時間パルス光が皮膚によって後方散乱した光を受光する。
ここで、観測対象である人の皮膚組織の構造について説明する。
図２は、人の皮膚組織の断面を示す模式図であり、皮膚３１は、表皮層３２と、真皮層（任意の層）３３と、皮下組織３４の３層により構成されている。
表皮層３２は、最も外側にある厚み０．２ｍｍ〜０．３ｍｍの薄い層で、概ね水を６０％程度、蛋白質、脂質及びグルコースを含有する層であり、角質層、顆粒層、有棘層、底層等を含む。
真皮層３３は、表皮層３２下に形成される厚み０．５ｍｍ〜２ｍｍの層で、概ね水を６０％程度、蛋白質、脂質及びグルコースを含有する層であり、この真皮層３３内には神経、毛根、皮脂腺、汗腺、毛包、血管、リンパ管等が存在する。
皮下組織３４は、真皮層３３下に形成される厚み１〜３ｍｍの層で、大部分が概ね脂質を９０％以上含み、残部が水からなる皮下脂肪でできている。
【００３８】
真皮層３３内には毛細血管等が発達しており、血中グルコースに応じた物質移動が速やかに起こり、血中グルコース濃度（血糖値）に対して真皮層３３中のグルコース濃度も追随して変化すると考えられている。そこで、この血糖値測定装置１００では、照射部１０４及び受光部１０５を所定の入出射間距離Ｗをおいて皮膚３１の表面に密着させ、この密着状態で照射部１０４から皮膚３１の表面に光を照射し、この光が皮膚３１内の組織により反射され、この反射光が照射部１０４及び受光部１０５に向かって散乱する光（後方散乱した光）を受光部１０５で検出する。
【００３９】
次に、計測光強度取得部１０６は、時刻ｔにおいて受光部１０５が受光した光の強度を取得する。次に、光路長取得部１０７は、光路長分布記憶部１０２から、皮膚モデルにおける光路長分布の時刻ｔにおける皮膚の各層の光路長を取得し、無吸収時光強度取得部１０８は、時間分解波形記憶部１０３から、皮膚モデルにおける短時間パルス光の時間分解波形の時刻ｔにおける光の強度を取得する。
【００４０】
次に、光吸収係数算出部１０９は、計測光強度取得部１０６が取得した光強度と光路長取得部１０７が取得した皮膚の各層の光路長と無吸収時光強度取得部１０８が取得した光強度とに基づいて、皮膚の真皮層の光吸収係数を算出し、濃度算出部１１０は、光吸収係数算出部１０９が算出した光吸収係数に基づいて、真皮層におけるグルコースの濃度を算出する。
【００４１】
この血糖値測定装置１００では、皮膚３１の特定深さの層である真皮層３３の光吸収係数を用いてグルコース濃度を算出するので、真皮層３３以外の層によるノイズの影響を軽減することができ、特定深さの層である真皮層３３に含まれるグルコースの濃度を算出することができる。したがって、特定深さの層に含まれるグルコースの光の吸収量、すなわちグルコースの濃度を、非侵襲的に精度良く測定することができる。
【００４２】
次に、血糖値測定装置１００の動作を説明する。
血糖値測定装置１００は、血糖値を測定する前に、予め皮膚モデルの各層における光路長分布と時間分解波形とを算出しておく必要がある。
まず、皮膚モデルの光路長分布及び時間分解波形の算出方法を説明する。
初めに、シミュレーション部１０１は、皮膚モデルを生成する。皮膚モデルの生成は、皮膚の各層の光散乱係数、光吸収係数及び厚みを決定することで行う。ここで、皮膚の部分を特定すれば、この特定された皮膚の部分における各層の散乱係数及び厚みは、個体による差が少ないので、予めサンプルを取ることなどによって決定するとよい。なお、表皮層３２の厚みは略０．３ｍｍ、真皮層３３の厚みは略１．２ｍｍ、皮下組織３４の厚みは略３．０ｍｍである。
また、ここで用いる皮膚モデルの光吸収係数はゼロとする。その理由は、この皮膚モデルを用いて光吸収量を算出するからである。
【００４３】
シミュレーション部１０１は、皮膚モデルを生成すると、この皮膚モデルに光を照射するシミュレーションを行う。このとき、照射部１０４の位置と受光部１０５の位置との間の距離Ｗを決定しておく必要がある。シミュレーションは、例えば、モンテカルロ法を用いて行う。
このモンテカルロ法によるシミュレーションは、例えば以下のように行われる。
【００４４】
まず、シミュレーション部１０１は、照射する光のモデルを光子（光束）とし、この光子を皮膚モデルに照射する計算を行う。皮膚モデルに照射された光子は、皮膚モデル内を移動する。このとき、光子は、次に進む点までの距離Ｌ及び方向θを乱数Ｒによって決定する。シミュレーション部１０１は、光子が次に進む点までの距離Ｌの計算を、式（４）により行う。
【００４５】
【数４】
0005521199.tif 000005
【００４６】
但し、ｌｎ（Ａ）はＡの自然対数を示し、μｓは、皮膚モデルの第ｓ層（表皮層、真皮層、皮下組織層の何れか）の散乱係数を示す。
また、シミュレーション部１０１は、光子が次に進む点までの方向θの計算を、式（５）により行う。
【００４７】
【数５】
0005521199.tif 000006
【００４８】
但し、ｇは、散乱角度の余弦（ｃｏｓ）の平均である非等方性パラメーターを示し、皮膚の非等方性パラメーターは、略０．９である。
シミュレーション部１０１は、上記式（４）、式（５）の計算を単位時間毎に繰り返すことにより、照射部１０４から受光部１０５までの光子の移動経路を算出することができる。シミュレーション部１０１は、複数の光子について移動距離の算出を行う。例えば、シミュレーション部１０１は、１０８個の光子について移動距離を算出する。
【００４９】
図３は、シミュレーション部１０１が算出した各層の光路長分布を示す図である。
図３では、横軸を光子の照射からの経過時間とし、縦軸を光路長の対数表示としている。
シミュレーション部１０１は、受光部１０５に到達した光子の各々の移動経路を、移動経路が通過する層毎に分類する。そして、シミュレーション部１０１は、単位時間毎に到達した光子の移動経路の平均長を分類された層毎に算出することで、図３に示すような皮膚の各層の光路長分布を算出する。
【００５０】
また、シミュレーション部１０１は、単位時間毎に受光部１０４に到達した光子の個数を算出することで、図４に示すような皮膚モデルの時間分解波形を算出する。
図４は、シミュレーション部１０１が算出した時間分解波形を示す図である。図４では、横軸を光子の照射からの経過時間とし、縦軸を受光部８が検出した光子数としている。
上述したような処理により、シミュレーション部１０１は、複数の波長に対して、皮膚モデルの光路長分布及び時間分解波形を算出する。このとき、シミュレーション部１０１は、皮膚の主成分（水、たんぱく質、脂質、グルコース等）の吸収スペクトルの直交性高くなる波長の光、すなわち皮膚の主成分（水、たんぱく質、脂質、グルコース等）の各々の成分のうち、ある主成分における特定成分の吸収スペクトルの極大値が他の成分の吸収スペクトルの極大値と大きく異なる波長の光について、光路長分布及び時間分解波形を算出すると良い。
【００５１】
図５は、皮膚の主成分の吸収スペクトルを示すグラフである。この図５では、横軸を照射する光の波長とし、縦軸を吸収係数としている。
図５によれば、グルコースの吸収係数は、波長が１６００ｎｍのときに極大となり、水の吸収係数は波長が１４５０ｎｍのときに極大となることがわかる。
そのため、シミュレーション部１０１は、例えば１４５０ｎｍ、１６００ｎｍというように、皮膚を構成する主成分の各々の成分の吸収スペクトル分布の直交性が高くなる波長の光、すなわち、皮膚を構成する主成分の各々の成分のうち、ある主成分における特定成分の吸収スペクトルの極大値が他の成分の吸収スペクトルの極大値と大きく異なる波長の光について、光路長分布及び時間分解波形を算出すると良い。
【００５２】
シミュレーション部１０１が複数の波長に対する皮膚モデルの光路長分布及び時間分解波形を算出すると、光路長分布記憶部１０２は、算出された光路長分布の情報を記憶し、時間分解波形記憶部１０３は、算出された時間分解波形の情報を記憶する。なお、シミュレーション部１０１は血糖測定装置１００に含まれない構成としても良い。この場合、外部の装置にてシミュレーションを行った結果を光路長分布記憶部１０２と時間分解波形記憶部１０３に記憶することで血糖測定装置１００が血糖値測定を行う。
【００５３】
次に、この血糖値測定装置１００が血糖値を測定する動作について、図６に基づき説明する。
まず、ユーザー（被測定者）が血糖値測定装置１００を手首等の皮膚に当て、測定開始スイッチ（図示せず）の押下等により血糖値測定装置１００を動作させると、照射部１０４は、皮膚３１に対して波長λ１の短時間パルス光を照射する（ステップＳ１）。
【００５４】
ここで、波長λ１は、シミュレーション部１０１が光路長分布及び時間分解波形を算出した複数の波長の中の１つである。
例えば、皮膚を構成する主成分のうち、ある主成分における特定成分の光吸収係数が他の成分の光吸収係数より大きくなる波長の光、すなわち、特定成分の光吸収係数の極小値が他の成分の光吸収係数の極小値と大きく異なる波長の光について光路長分布及び時間分解波形を算出すると良い。
【００５５】
照射部１０４が短時間パルス光を照射すると、受光部１０５は、照射部１０４から照射され皮膚３１によって後方散乱された光を受光する（ステップＳ２）。
このとき、受光部１０５は、照射開始からの単位時間毎（例えば、１ピコ秒毎）の受光強度を内部メモリー（図示せず）に記録しておく。
【００５６】
次いで、受光部１０５が受光を完了すると、計測光強度取得部１０６は、内部メモリーに格納されている、異なる時刻ｔにおける受光強度Ｉ（ｔ）を皮膚の層の数と同じ数だけ取得する（ステップＳ３）。
図６のフローチャートでは、皮膚の３つの層について４種類の波長を用いて濃度測定を行う場合について説明する。
すなわち、計測光強度取得部１０６は、皮膚の３つの層について４種類の波長を用いて濃度測定を行う場合には、３つの異なる時刻ｔ１〜ｔ３における受光強度Ｉ（ｔ１）〜Ｉ（ｔ３）を取得する。ここで、皮膚の層の数と同じ数だけ受光強度を取得する理由は、後述する処理において、皮膚の各層の吸収係数を連立方程式によって算出するためである。
【００５７】
また、計測光強度取得部が光強度を取得する時刻ｔ１〜ｔ３は、皮膚の各層の光路長分布のピークとなる時刻であると良い。すなわち、照射部１０４が短時間パルス光を照射した時刻に、図３に示すグラフにおいて皮膚の各層の光路長が極大となる時間を加算した時刻の光強度をそれぞれ取得すると良い。
【００５８】
計測光強度取得部１０６が受光強度Ｉ（ｔ１）〜Ｉ（ｔ３）を取得すると、光路長取得部１０７は、光路長分布記憶部１０２が記憶する波長λ１の光路長分布から、時刻ｔ１〜ｔ３における皮膚の各層の光路長Ｌ１（ｔ１）〜Ｌ１（ｔ３）、Ｌ２（ｔ１）〜Ｌ２（ｔ３）、Ｌ３（ｔ１）〜Ｌ３（ｔ３）を取得する（ステップＳ４）。
また、計測光強度取得部１０６が受光強度Ｉ（ｔ１）〜Ｉ（ｔ３）を取得すると、無吸収時光強度取得部１０８は、時間分解波形記憶部１０３が記憶する波長λ１の時間分解波形から、時刻ｔ１〜ｔ３における検出光子数Ｎ（ｔ１）〜Ｎ（ｔ３）を取得する（ステップＳ５）。
【００５９】
光路長取得部１０７が皮膚の各層の光路長を取得し、無吸収時光強度取得部１０８が検出光子数を取得すると、光吸収係数算出部１０９は、式（６）に基づいて、皮膚の各層の光吸収係数μ１〜μ３を算出する（ステップＳ６）。ここで、光吸収係数μ１は、表皮層の光吸収係数を示し、光吸収係数μ２は、真皮層の光吸収係数を示し、光吸収係数μ３は、皮下組織の光吸収係数を示す。
【００６０】
【数６】
0005521199.tif 000007
【００６１】
但し、ｌｎ（Ａ）はＡの自然対数を示し、Ｎ（ｔ）は特定波長λｋの短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度を示す。また、Ｉｉｎは、照射部７が照射した短時間パルス光の光強度を示す。また、Ｎｉｎは、シミュレーション部２が照射のシミュレーションを行った光子の個数を示す。
【００６２】
光吸収係数算出部１０９が皮膚の各層の光吸収係数μ１〜μ３を算出すると、光吸収係数算出部１０９は、皮膚の主成分の種類数と同じ数の波長に対して光吸収係数μ１〜μ３を算出したか否かを判定する（ステップＳ７）。
本実施形態では、皮膚の主成分を水、たんぱく質、脂質、グルコースの４種類として血糖値の測定を行うため、光吸収係数算出部１０９は、４種類の波長λ１〜λ４に対して光吸収係数μ１〜μ３を算出したか否かを判定する。ここで、波長λ１〜λ４は、シミュレーション部１０１が光路長分布及び時間分解波形を算出した複数の波長の中から選択する。
【００６３】
光吸収係数算出部１０９が、光吸収係数μ１〜μ３を算出していない波長λ１〜λ４があると判定した場合（ステップＳ７：ＮＯ）、ステップＳ１に戻り、まだ光吸収係数μ１〜μ３を算出していない波長λ１〜λ４の光吸収係数μ１〜μ３の算出を行う。
他方、光吸収係数算出部１０９が、波長λ１〜λ４の光吸収係数μ１〜μ３を算出していると判定した場合（ステップＳ７：ＹＥＳ）、濃度算出部１１０は、式（７）に基づいて真皮層に含まれるグルコースの濃度を算出する（ステップＳ８）。
【００６４】
【数７】
0005521199.tif 000008
【００６５】
但し、μ２（１）〜μ２（４）は、真皮層における波長λ１〜λ４の光吸収係数を示す。また、ｇ１〜ｇ４は、真皮層におけるそれぞれ皮膚の主成分である水、たんぱく質、脂質、グルコースのモル濃度を示す。また、ε１（１）〜ε１（４）は、波長λ１〜λ４に対する水のモル吸光係数を示し、ε２（１）〜ε２（４）は、波長λ１〜λ４に対するたんぱく質のモル吸光係数を示し、ε３（１）〜ε３（４）は、波長λ１〜λ４に対する脂質のモル吸光係数を示し、ε４（１）〜ε４（４）は、波長λ１〜λ４に対するグルコースのモル吸光係数を示す。
つまり、式（７）のｇ４を算出することで、真皮層に含まれるグルコースのモル濃度を求めることができる。
【００６６】
ここで、式（７）によりグルコースのモル濃度を求めることができる理由を説明する。皮膚の散乱係数の波長依存性は小さいので、検出光子数Ｎ（ｔ）及び光路長Ｌｎ（ｔ）の波長に対する変化は無視することができる。また、ベア・ランベルト（Beer-Lambert）の法則により、吸光度＝モル吸光係数×モル濃度で表すことができる。これにより、２波長で得られた時間分解計測より、検出光子数Ｎ（ｔ）を消去することで、真皮層において得られる吸収係数差と皮膚を形成する各成分のモル吸光係数との関係式を示す式（７）を導くことができる。
【００６７】
本実施形態によれば、短時間パルス光を照射し、所定の時刻において受光した光の強度に基づいてグルコースの濃度を定量する。これにより、所定の時刻において受光した光の中から、真皮層の吸収係数を選択的に算出することができる。これにより、特定の皮膚の層のグルコースの濃度を算出することができ、他の層によるノイズの影響を軽減し、精度の高い血糖値を算出することが可能になる。
【００６８】
［第２の実施形態］
本発明の第２の実施形態について詳しく説明する。
第２の実施形態における血糖値測定装置の構成は、第１の実施形態の血糖値測定装置１００の構成と同一であり、計測光強度取得部１０６、光路長取得部１０７、無吸収時光強度取得部１０８、光吸収係数算出部１０９の動作が異なる。
【００６９】
図７は、血糖値測定装置が血糖値を測定する動作を示すフローチャートである。
まず、血糖値測定装置１００を動作させると、照射部１０４は、皮膚に対して波長λ１の短時間パルス光を照射する（ステップＳ１１）。ここで、波長λ１は、シミュレーション部１０１が光路長分布及び時間分解波形を算出した複数の波長の中の１つである。
【００７０】
照射部１０４が短時間パルス光を照射すると、受光部１０５は、照射部１０４から照射され、皮膚によって後方散乱した光を受光する（ステップＳ１２）。このとき、受光部１０５は、照射開始からの単位時間毎（例えば、１ピコ秒毎）の受光強度を内部メモリーに登録しておく。
受光部１０５が受光を完了すると、計測光強度取得部１０６は、内部メモリーに格納されている受光強度から、ある時刻から時間τの間の受光強度の時間分布を取得する（ステップＳ１３）。
【００７１】
計測光強度取得部１０６が、時間τの間の受光強度の時間分布を取得すると、光路長取得部１０７は、光路長分布記憶部１０２が記憶する波長λ１の光路長分布から、ある時刻から時間τの間の皮膚の各層の光路長を取得する。
図７のフローチャートでは、皮膚の３つの層について４種類の波長を用いて濃度測定を行う場合について説明する。
すなわち、計測光強度取得部１０６は、皮膚の３つの層について４種類の波長を用いて濃度測定を行う場合には、ある時刻から時間τの間の皮膚の各層の光路長Ｌ１〜Ｌ３を取得する（ステップＳ１４）。
【００７２】
また、計測光強度取得部１０６が、時間τの間の受光強度の時間分布を取得すると、無吸収時光強度取得部１０８は、時間分解波形記憶部１０３が記憶する波長λ１の時間分解波形から、ある時刻から時間τの間の検出光子数を取得する（ステップＳ１５）。
【００７３】
光路長取得部１０７が皮膚の各層の光路長を取得し、無吸収時光強度取得部１０８が検出光子数を取得すると、光吸収係数算出部１０９は、式（８）に基づいて、皮膚の各層の光吸収係数μ１〜μ３を算出する（ステップＳ１６）。ここで、光吸収係数μ１は、表皮層の光吸収係数を示し、光吸収係数μ２は、真皮層の光吸収係数を示し、光吸収係数μ３は、皮下組織の光吸収係数を示す。
【００７４】
【数８】
0005521199.tif 000009
【００７５】
但し、ｌｎ（Ａ）は、Ａの自然対数を示す。また、Ｉ（ｔ）は、時刻ｔにおける受光部１０５の受光強度を示し、Ｉｉｎは、照射部１０４が照射した短時間パルス光の光強度を示す。また、Ｎ（ｔ）は、時間分解波形の時刻ｔにおける検出光子数を示し、Ｎｉｎは、シミュレーション部１０１が照射のシミュレーションを行った光子の個数を示す。また、Ｌ１（ｔ）〜Ｌ３（ｔ）は、時刻ｔにおける皮膚の各層の光路長を示す。
【００７６】
光吸収係数算出部１０９が皮膚の各層の光吸収係数μ１〜μ３を算出すると、光吸収係数算出部１０９は、皮膚の主成分の種類数と同じ数の波長に対して光吸収係数μ１〜μ３を算出したか否かを判定する（ステップＳ１７）。本実施形態では、皮膚の主成分を水、たんぱく質、脂質、グルコースの４種類として血糖値の測定を行うので、光吸収係数算出部１０９は、４種類の波長λ１〜λ４に対して光吸収係数μ１〜μ３を算出したか否かを判定する。ここで、波長λ１〜λ４は、シミュレーション部１０１が光路長分布及び時間分解波形を算出した複数の波長の中から選択する。
【００７７】
ここで、光吸収係数算出部１０９が光吸収係数μ１〜μ３を算出していない波長λ１〜λ４があると判定した場合（ステップＳ１７：ＮＯ）、ステップＳ１に戻り、まだ光吸収係数μ１〜μ３を算出していない波長λ１〜λ４の光吸収係数μ１〜μ３の算出を行う。
他方、光吸収係数算出部１０９が波長λ１〜λ４の光吸収係数μ１〜μ３を算出していると判定した場合（ステップＳ１７：ＹＥＳ）、濃度算出部１１０は、上述した式（７）に基づいて真皮層に含まれるグルコースの濃度を算出する（ステップＳ１８）。
【００７８】
本実施形態によれば、吸収係数μ１〜μ３を、時間τの間の光路長の積分値によって算出する。これにより、計測した受光強度Ｉ（ｔ）に含まれている誤差による吸収係数μ１〜μ３の算出結果に対する影響を少なくすることができる。
【００７９】
第１の実施形態及び第２の実施形態について、図面を参照して説明してきたが、具体的な構成は上述のものに限られることはなく、この発明の要旨を逸脱しない範囲内において様々な設計変更等をすることが可能である。
例えば、第１の実施形態及び第２の実施形態では、濃度定量方法を血糖値測定装置１００に実装し、皮膚の真皮層に含まれるグルコースの濃度を測定する場合を説明したが、これに限られず、濃度定量方法を、複数の層から形成される観測対象の任意の層における目的成分の濃度を定量する他の装置に用いても良い。
【００８０】
また、第１の実施形態及び第２の実施形態では、血糖値測定装置１００がシミュレーション部１０１と、光路長分布記憶部１０２と、時間分解波形記憶部１０３とを備えた構成としたが、シミュレーション部１０１におけるシミュレーションの結果を、光路長分布記憶部１０２及び時間分解波形記憶部１０３に記憶させた構成とすれば、シミュレーション部１０１を別途設ける必要がなくなる。
【００８１】
上述の血糖値測定装置１００は、内部にコンピューターシステムを備えた構成としてもよい。そして、上述した各処理部の動作または各処理部の一部の動作は、プログラムの形式でコンピューター読み取り可能な記録媒体に記憶されており、このプログラムをコンピューターが読み出して実行することによって、上記処理が行われる。ここで、コンピューターが読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリー等をいう。また、このコンピュータープログラムを通信回線によってコンピューターに配信し、この配信を受けたコンピューターが当該プログラムを実行するようにしても良い。
【００８２】
また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。
さらに、前述した機能をコンピューターシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であっても良い。
【００８３】
［第３の実施形態］
図８は、本発明の第３の実施形態の血糖値測定装置の構成を示す概略ブロック図である。
本実施形態の血糖値測定装置２００が、第１の実施形態の血糖値測定装置１００と異なる点は、本実施形態の血糖値測定装置２００が、照射部１０４と、受光部１０５と、コンピューターシステム２０１とを備え、コンピューターシステム２０１は、記憶部２０２と、制御部（ＣＰＵ）２０３とを備え、記憶部２０２は、シミュレーション部１０１におけるシミュレーションの結果を記憶している光路長分布記憶部１０２及び時間分解波形記憶部１０３の機能を実行し、制御部（ＣＰＵ）２０３は、計測光強度取得部１０６、光路長取得部１０７、無吸収時光強度取得部１０８、光吸収係数算出部１０９及び濃度算出部１１０の機能を実行する点である。
【００８４】
本実施形態の血糖値測定装置２００においても、第１の実施形態の血糖値測定装置１００と同様の作用・効果を奏することができる。
【００８５】
［第４の実施形態］
図９は、本発明の第４の実施形態の血糖値測定装置の構成を示す概略ブロック図である。
この血糖値測定装置３００は、手のひら等の皮膚（観測対象）を構成する複数層のうちの真皮層（任意の層）に含まれるグルコース（目的成分）の濃度を非侵襲にて定量する装置であり、シミュレーション部３０１と、光路長分布記憶部３０２と、時間分解波形記憶部３０３と、照射部３０４と、受光部３０５と、計測光強度取得部（光強度取得部）３０６と、光路長取得部３０７と、無吸収時光強度取得部（光強度モデル取得部）３０８と、積分区間算出部３０９と、光吸収係数算出部３１０と、光吸収係数分布記憶部３１１と、光吸収係数取得部３１２と、濃度算出部３１３とを備えている。
【００８６】
シミュレーション部３０１は、光吸収係数がゼロの皮膚モデルに対して光を照射するシミュレーションを行う。
光路長分布記憶部３０２は、皮膚に対して照射する短時間パルス光の、この皮膚を構成する各々の層における光路長分布のモデルを記憶する。ここでは、光吸収係数がゼロの皮膚モデルの光路長分布を記憶する。
ここで、短時間パルス光とは、パルス幅の時間が照射部３０４から受光部３０５へ光が空気中を直接伝搬する時間よりも短いパルス光のことであり、例えば、パルス光の半値幅が０．１ｐｓ〜１０ｐｓ、２つのパルス光の間の時間間隔が１ｐｓ〜１００ｐｓのパルス光のことである。
また、光路長分布とは、光（光子）の移動経路の長さ（光路長）を当該光（光子）が受光部３０５に到達するまでの時間を基に分布関数として表したものである。
【００８７】
光路長取得部３０７は、光路長分布記憶部３０２から、光路長分布のモデルの所定の時刻における、皮膚の各々の層の光路長を取得する。ここでは、光路長分布記憶部３０２からある時刻における光路長を取得する。ここでいう光路長とは、照射部３０４から照射された短時間パルス光が皮膚内に侵入し、この皮膚内にて散乱されて受光部３０５により検出されるまでの光の路の長さのことであり、後述するように、照射部３０４と受光部３０５との距離を設定することにより、皮膚の各々の層の光路長が推定される。
【００８８】
時間分解波形記憶部３０３は、皮膚に対して照射する短時間パルス光の時間分解波形のモデルを記憶する。ここでは、光吸収係数がゼロの皮膚モデルの時間分解波形を記憶する。ここで、短時間パルス光の時間分解波形とは、受光部３０５にて受光した光（光子）の強度を、この短時間パルス光の照射時からの経過時間を基に分布関数として表したものである。
無吸収時光強度取得部３０８は、時間分解波形記憶部３０３から、短時間パルス光の時間分解波形のモデルの所定の時刻における光の吸収係数をゼロ（零）とした時（無吸収時）の光強度を取得する。ここでは、時間分解波形記憶部３０３からある時刻における無吸収時光強度を取得する。
【００８９】
照射部３０４は、皮膚に対して短時間パルス光を照射する。この光は皮膚内の組織によって散乱され、皮膚内に拡散する。拡散した光の一部は、受光部３０５に到達する（後方散乱光）。この受光部３０５に到達した後方散乱光が皮膚内を伝搬してきた経路は、図２に示すようにバナナシェイプの経路となる。この照射部３０４が照射する複数の短時間パルス光は、皮膚を構成する主成分の各々の成分の吸収スペクトル分布の直交性が高くなる波長の光、すなわち、皮膚を構成する主成分の各々の成分のうち、ある主成分における特定成分の吸収スペクトルの極大値が他の成分の吸収スペクトルの極大値と大きく異なる波長の光を含んでいる。
受光部３０５は、短時間パルス光が皮膚によって後方散乱した光を受光する。この受光部３０５は、受光強度を記録する内部メモリー（図示せず）を備えている。なお、この内部メモリーは、受光部３０５に電気的に接続する外部メモリーに代えた構成としてもよい。
【００９０】
この照射部３０４と受光部３０５との入出射間距離Ｗと皮膚内に侵入する光Ｒの侵入深さとの間には、一定の関係がある。そこで、照射部３０４と受光部３０５との入出射間距離Ｗを規定することにより、皮膚内に侵入する光Ｒの侵入深さも一義的に決定されることとなる。例えば、入出射間距離Ｗを１０ｍｍとすると、光Ｒの侵入深さは１０ｍｍとなり、入出射間距離Ｗを０．８ｍｍとすると、光Ｒの侵入深さは０．８ｍｍとなる。
【００９１】
受光部３０５が後方散乱光を受光すると、計測光強度取得部３０６は、受光部３０５が受光した光の積分区間内の皮膚の各々の層の数ｎ以上となる複数の時刻ｔ１〜ｔｍにおける光強度を取得する（但し、ｎは１以上の自然数、ｍはｎ以上の自然数）。計測光強度取得部３０６が光強度を取得する複数の時刻は、皮膚の各々の層の光路長分布のピーク時間を含んでいる。
積分区間算出部３０９は、光路長取得部３０７が取得した光路長分布のモデルの皮膚の各々の層の光路長と、無吸収時光強度取得部３０８が取得した短時間パルス光の時間分解波形のモデルの無吸収時光強度と、計測光強度取得部３０６が取得した受光部３０５が受光した光強度分布とに基づいて、前記光強度分布から皮膚の各々の層の光強度に対応する領域の積分区間を算出する。
【００９２】
ここで、積分区間とは、光強度分布における任意の層の光強度に対応する領域の時間幅のことであり、開始時刻と、終了時刻と、増分時間とにより決定することができる。
例えば、（１）後方散乱した光を受光する受光部３０５の出力する光強度が計測光強度取得部３０６の最小検出感度を超えて検出された時刻から最小検出感度と等しい光強度で検出された時刻までの時間、（２）シミュレーション部３０１で得られる無吸収時光強度を記憶している時間分解波形記憶部３０３から取得した無吸収時光強度の時間特性、（３）皮膚表面に接する受光部３０５と照射部３０４との間隔、（４）シミュレーション部３０１に与える皮膚モデルのサイズ及び光学特性（散乱係数、吸収係数、非等方性パラメータ−、または屈折率）を用いて、積分区間の開始時刻、終了時刻、増分時間を決定する。
【００９３】
光吸収係数算出部３１０は、積分区間算出部３０９が算出した積分区間に基づいて皮膚の各層の光吸収係数を算出する。
ここでは、積分区間算出部３０９によって定めた積分区間での皮膚の各層の光吸収係数及び推定誤差率を求め、積分区間に対する皮膚の各層の光吸収係数及び推定誤差率の分布を算出する。
この光吸収係数算出部３１０では、皮膚における各層の光吸収係数を、下記の式（９）から算出する。
【００９４】
【数９】
0005521199.tif 000010
【００９５】
但し、ｌｎ（Ａ）はＡの自然対数、Ｉ（ｔ）は受光部３０５が時刻ｔにて受光した光強度、Ｎ（ｔ）は特定波長λｋの短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度、Ｌｉ（ｔ）は皮膚の各々の層における光路長分布のモデルの時刻ｔにおける第ｉ層の光路長、μｉは第ｉ層の光吸収係数を示す。
ここで、第１層は表皮層、第２層は真皮層、第３層は皮下組織を示し、μ１は表皮層の光吸収係数、μ２は真皮層の光吸収係数、μ３は皮下組織の光吸収係数を示す。
【００９６】
光吸収係数分布記憶部３１１は、光吸収係数算出部３１０が算出した積分区間に対する任意の層の光吸収係数及び推定誤差率の分布を記憶する。
光吸収係数取得部３１２は、光吸収係数分布記憶部３１１から取得した積分区間に対する任意の層の光吸収係数及び推定誤差率の分布と、積分区間の変化に対する光吸収係数変動率の範囲等の基準とを用いて、皮膚の表面からの特定深さにおける層の皮膚主成分における血糖相応のグルコース濃度に基づく光吸収係数を取得する。
【００９７】
濃度算出部３１３は、光吸収係数取得部３１２が取得した皮膚の表面からの特定深さにおける層の皮膚主成分における血糖相応のグルコース濃度に基づく光吸収係数から、特定深さの層に含まれるグルコースの濃度を算出する。
この濃度算出部３１３では、皮膚の任意の層におけるグルコースの濃度を、下記の式（１０）から算出する。
【００９８】
【数１０】
0005521199.tif 000011
【００９９】
但し、μａは皮膚の任意の層である第ａ層における光吸収係数、ｇｊは皮膚を構成する第ｊ成分のモル濃度、εｊは第ｊ成分の光吸収係数、ｐは皮膚を構成する主成分の個数、ｑは特定波長λｋの種類数を示す。
ここで、第１層は表皮層、第２層は真皮層、第３層は皮下組織を示し、μ１は表皮層の光吸収係数、μ２は真皮層の光吸収係数、μ３は皮下組織の光吸収係数を示す。
【０１００】
そして、この血糖値測定装置３００では、照射部３０４は、皮膚に短時間パルス光を照射し、受光部３０５は、短時間パルス光が皮膚により後方散乱した光を受光し、計測光強度取得部３０６は、時刻ｔにおいて受光部３０５が受光した光強度を取得する。
次いで、光路長取得部５は、光路長分布記憶部４から、皮膚モデルにおける光路長分布の時刻ｔにおける皮膚の各層の光路長を取得し、無吸収時光強度取得部３０８は、時間分解波形記憶部３０３から、皮膚モデルにおける短時間パルス光の時間分解波形の時刻ｔにおける光の強度を取得する。
【０１０１】
次いで、積分区間算出部３０９は、光路長取得部３０７が取得した光路長分布のモデルの皮膚の各々の層の光路長、例えば、光路長分布記憶部３０２から取得した光路長と、無吸収時光強度取得部３０８が取得した短時間パルス光の時間分解波形のモデルにおける光強度、例えば、時間分解波形記憶部３０３から取得した光強度と、計測光強度取得部３０６が取得した受光部３０５が受光した光の強度分布とに基づいて、前記光の強度分布から任意の層の光強度に対応する領域の積分区間を算出する。
【０１０２】
例えば、（１）後方散乱した光を受光する受光部３０５の出力する光強度が計測光強度取得部３０６の最小検出感度を超えて検出された時刻から最小検出感度と等しい光強度で検出された時刻までの時間、（２）シミュレーション部３０１で得られる無吸収時光強度を記憶している時間分解波形記憶部３０３から取得した無吸収時光強度の時間特性、（３）皮膚表面に接する受光部３０５と照射部３０４との間隔、（４）シミュレーション部３０１に与える皮膚モデルのサイズ及び光学特性（散乱係数、吸収係数、非等方性パラメーター、または屈折率）を用いて、積分区間の開始時刻、終了時刻、増分時間を決定する。
【０１０３】
次いで、光吸収係数算出部３１０は、積分区間算出部３０９によって定めた積分区間での任意の層の光吸収係数及び推定誤差率を求め、積分区間に対する任意の層の光吸収係数及び推定誤差率の分布を算出する。
次いで、光吸収係数分布記憶部３１１は、光吸収係数算出部３１０が算出した積分区間に対する任意の層の光吸収係数及び推定誤差率の分布を記憶する。
次いで、光吸収係数取得部３１２は、光吸収係数分布記憶部３１１から取得した積分区間に対する任意の層の光吸収係数及び推定誤差率の分布と、積分区間の変化に対する光吸収係数変動率の範囲等の基準とを用いて、皮膚の表面からの特定深さにおける層の皮膚主成分における血糖相応のグルコース濃度に基づく光吸収係数を取得する。
【０１０４】
次いで、濃度算出部３１３は、光吸収係数取得部３１２が取得した皮膚の表面からの特定深さにおける層の皮膚主成分における血糖相応のグルコース濃度に基づく光吸収係数に基づいて、皮膚の表面からの特定深さにおける層に含まれるグルコースの濃度を、上記の式（１０）に基づき算出する。
【０１０５】
以上により、この血糖値測定装置３００では、皮膚の特定深さの層である真皮層の光吸収係数を用いてグルコース濃度を算出するので、真皮層以外の層によるノイズの影響を軽減することができ、特定深さの層である真皮層に含まれるグルコースの濃度を精度良く算出することができる。
【０１０６】
また、積分区間算出部３０９が算出した積分区間に基づいて、この積分区間に対応する時間帯の特定深さの層である真皮層に含まれるグルコースの光吸収係数を取得するので、真皮層からの光に対する他の層からの光の影響を低減することができる。したがって、特定深さの層に含まれるグルコースの光の吸収量、すなわちグルコースの濃度を、非侵襲的に精度良く測定することができる。
【０１０７】
次に、血糖値測定装置３００の動作を説明する。
血糖値測定装置３００は、血糖値を測定する前に、予め皮膚モデルの各層における光路長分布と時間分解波形とを算出しておく必要がある。
ここで、皮膚モデルの光路長分布及び時間分解波形の算出方法を説明する。
初めに、シミュレーション部３０１は、皮膚モデルを生成する。皮膚モデルの生成は、皮膚の各層の光散乱係数、光吸収係数及び厚みを決定することで行う。ここで、皮膚の部分を特定すれば、この特定された皮膚の部分における各層の散乱係数及び厚みは、個体による差が少ないので、予めサンプルを取ることなどによって決定するとよい。なお、表皮層３２の厚みは略０．３ｍｍ、真皮層３３の厚みは略１．２ｍｍ、皮下組織３４の厚みは略３．０ｍｍである。
また、ここで用いる皮膚モデルの光吸収係数はゼロとする。その理由は、この皮膚モデルを用いて光吸収量を算出するからである。
【０１０８】
シミュレーション部３０１は、皮膚モデルを生成すると、この皮膚モデルに光を照射するシミュレーションを行う。このとき、照射部３０４の位置と受光部３０５の位置との間の距離を決定しておく必要がある。シミュレーションは、例えば、モンテカルロ法を用いて行う。
このモンテカルロ法によるシミュレーションは、例えば以下のように行われる。
【０１０９】
まず、シミュレーション部３０１は、照射する光のモデルを光子（光束）とし、この光子を皮膚モデルに照射する計算を行う。皮膚モデルに照射された光子は、皮膚モデル内を移動する。このとき、光子は、次に進む点までの距離Ｌ及び方向θを乱数Ｒによって決定する。シミュレーション部３０１は、光子が次に進む点までの距離Ｌの計算を、式（１１）により行う。
【０１１０】
【数１１】
0005521199.tif 000012
【０１１１】
但し、ｌｎ（Ａ）はＡの自然対数を示し、μｓは、皮膚モデルの第ｓ層（表皮層、真皮層、皮下組織層の何れか）の散乱係数を示す。
また、シミュレーション部３０１は、光子が次に進む点までの方向θの計算を、式（１２）により行う。
【０１１２】
【数１２】
0005521199.tif 000013
【０１１３】
ただし、ｇは、散乱角度の余弦（ｃｏｓ）の平均である非等方性パラメーターを示し、皮膚の非等方性パラメーターは、略０．９である。
シミュレーション部３０１は、上記式（１１）及び式（１２）の計算を単位時間毎に繰り返すことにより、照射部３０４から受光部３０５までの光子の移動経路を算出することができる。シミュレーション部３０１は、複数の光子について移動距離の算出を行う。例えば、シミュレーション部３０１は、１０８個の光子について移動距離を算出する。
【０１１４】
シミュレーション部３０１は、受光部３０５に到達した光子の各々の移動経路を、移動経路が通過する層毎に分類する。そして、シミュレーション部３０１は、単位時間毎に到達した光子の移動経路の平均長を分類された層毎に算出することで、既に説明した図３に示すような皮膚の各層の光路長分布を算出する。
また、シミュレーション部３０１は、単位時間毎に受光部３０５に到達した光子の個数を算出することで、既に説明した図４に示すような皮膚モデルの時間分解波形を算出する。
【０１１５】
図１０は、シミュレーション部３０１で得た無吸収時光強度（受光光子数と等しい）Ｎ（ｔ）の時間分解波形を示す図であり、図１１は、シミュレーション部３０１で得た吸収時光強度Ｉ（ｔ）の時間分解波形をリニアスケールで示した図であり、図１２は、同時間分解波形をログスケールで示した図である。図１１及び図１２では、表皮層、真皮層、皮下組織の三層の皮膚モデルを用いており、真皮の吸収係数に対して表皮の吸収係数を２５％から１５０％まで変化させている。図１１及び図１２中、Ａは２５％を、Ｂは５０％を、Ｃは７５％を、Ｄは１００％を、Ｅは１２５％を、Ｆは１５０％を、それぞれ示している。
図１０〜図１２によれば、３０ｐｓ程度で光子数または光強度が極大となる値（ピーク値）に対する振幅が５％以下となっており、真値に近いことが分かる。
また、図１１及び図１２によれば、５０ｐｓ程度では、真皮層の吸収係数に対して表皮の吸収係数を２５％とした場合、ピーク値に対して３桁程度の減少となっていることが分かる。
【０１１６】
上述したような処理により、シミュレーション部３０１は、複数の波長に対して、皮膚モデルの光路長分布及び時間分解波形を算出する。このとき、シミュレーション部３０１は、皮膚の主成分（水、たんぱく質、脂質、グルコース等）の吸収スペクトルの差が大きい波長について光路長分布及び時間分解波形を算出すると良い。
【０１１７】
皮膚の主成分の吸収スペクトルについては、既に図５にて説明しているとおり、グルコースの吸収係数は、波長が１６００ｎｍのときに極大となり、水の吸収係数は波長が１４５０ｎｍのときに極大となる。そのために、シミュレーション部３０１は、例えば１４５０ｎｍ、１６００ｎｍというように、皮膚を構成する主成分の各々の成分の吸収スペクトル分布の直交性が高くなる波長の光、すなわち、皮膚を構成する主成分の各々の成分のうち、ある主成分における特定成分の吸収スペクトルの極大値が他の成分の吸収スペクトルの極大値と大きく異なる波長の光について光路長分布及び時間分解波形を算出すると良い。
【０１１８】
シミュレーション部３０１が複数の波長に対する皮膚モデルの光路長分布及び時間分解波形を算出すると、光路長分布記憶部３０２は、算出された光路長分布の情報を記憶し、時間分解波形記憶部３０３は、算出された時間分解波形の情報を記憶する。
【０１１９】
次に、この血糖値測定装置３００が血糖値を測定する動作について、図１３及び図１４に基づき説明する。
まず、被測定者が血糖値測定装置３００を手首等の皮膚に当て、測定開始スイッチ（図示せず）の押下等により血糖値測定装置３００を動作させると、照射部３０４は、皮膚に対して波長λｋの短時間パルス光を照射する（ステップＳ２１）。
この波長λｋとしては、例えば、シミュレーション部３０１が光路長分布及び時間分解波形を算出した複数の波長の中の１つが好ましい。
例えば、皮膚を構成する主成分のうち、ある主成分における特定成分の光吸収係数が他の成分の光吸収係数より大きくなる波長の光、すなわち、特定成分の光吸収係数の極小値が他の成分の光吸収係数の極小値と大きく異なる波長の光について光路長分布及び時間分解波形を算出すると良い。
【０１２０】
照射部３０４が波長λｋの短時間パルス光を照射すると、受光部３０５は、照射部３０４から照射され皮膚によって後方散乱された光を受光する（ステップＳ２２）。
このとき、受光部３０５は、照射開始からの単位時間毎（例えば、１ピコ秒毎の時刻ｔ１〜ｔｍ）の受光強度を、内部メモリー（図示せず）に記録しておく。
【０１２１】
次いで、積分区間を変化させて真皮層の光吸収係数を算出する（処理Ａ；ステップＳ２３）。
このステップＳ２３は、図１４に示す動作により行う。
まず、積分区間算出部３０９により、（１）後方散乱した光を受光する受光部３０５の出力する光強度が計測光強度取得部３０６の最小検出感度を超えて検出された時刻から最小検出感度と等しい光強度で検出された時刻までの時間、（２）シミュレーション部３０１で得られる無吸収時光強度を記憶している時間分解波形記憶部３０３から取得した無吸収時光強度の時間特性、（３）皮膚表面に接する受光部３０５と照射部３０４との間隔、（４）シミュレーション部３０１に与える皮膚モデルのサイズ及び光学特性（散乱係数、吸収係数、非等方性パラメーター、または屈折率）を用いて、積分区間を算出する。より具体的には、積分区間の開始時刻、終了時刻、増分時間を算出する（ステップＳ３１）。
【０１２２】
図１５は、真皮層における光吸収係数と積分区間との関係を示す図である。図１５では、表皮層、真皮層、皮下組織の三層の皮膚モデルを用いており、真皮層の吸収係数に対して表皮層の吸収係数を２５％から１５０％まで変化させている。図１５中、Ａは２５％を、Ｂは５０％を、Ｃは７５％を、Ｄは１００％を、Ｅは１２５％を、Ｆは１５０％を、それぞれ示している。
【０１２３】
図１５によれば、１４ｐｓ〜３４ｐｓの範囲では、表皮層と真皮層の光吸収係数の値が近づくと、予め与えられた真値の光吸収係数の値である０．５５／ｍｍに近い値で算出されていることが分かる。また、３４ｐｓ〜５４ｐｓの範囲では、増加してピークを持つ傾向があることが分かる。これにより、積分区間が広くなるにしたがって、真皮層の光吸収係数の値はある特定の範囲（図１５では、０．５４〜０．５６／ｍｍ）に収斂していくことが分かる。
【０１２４】
次いで、受光部３０５が受光を完了すると、計測光強度取得部３０６は、内部メモリーに記録されている受光強度から、ある時刻ｔにおける受光強度を皮膚の層の数と同じ数だけ取得する（ステップＳ３２）。
例えば、皮膚の３つの層について４種類の波長を用いて濃度測定を行う場合には、３つの異なる時刻ｔ１〜ｔ３における受光強度Ｉ（ｔ１）〜Ｉ（ｔ３）を取得する。ここで、皮膚の層の数と同じ数だけ受光強度を取得する理由は、後述する処理において、皮膚の各層の吸収係数を連立方程式によって算出するためである。
【０１２５】
また、計測光強度取得部３０６が光強度を取得する時刻ｔ１〜ｔ３は、皮膚の各層の光路長分布のピークとなる時刻であると良い。すなわち、照射部３０４が短時間パルス光を照射した時刻に、既に説明した図３において皮膚の各層の光路長が極大となる時間を加算した時刻の光強度をそれぞれ取得すると良い。
【０１２６】
計測光強度取得部３０６が、受光強度Ｉ（ｔ１）〜Ｉ（ｔ３）を取得すると、光路長取得部３０７は、光路長分布記憶部３０２が記憶する波長λ１の光路長分布から、時刻ｔ１〜ｔ３における皮膚の各層の光路長Ｌ１（ｔ１）〜Ｌ１（ｔ３）、Ｌ２（ｔ１）〜Ｌ２（ｔ３）、Ｌ３（ｔ１）〜Ｌ３（ｔ３）を取得する（ステップＳ３３）。
また、計測光強度取得部３０６が、受光強度Ｉ（ｔ１）〜Ｉ（ｔ３）を取得すると、無吸収時光強度取得部３０８は、時間分解波形記憶部３０３が記憶する波長λ１の時間分解波形から、短時間パルス光の時間分解波形のモデルの所定の時刻における光強度、例えば、時刻ｔ１〜ｔ３における検出光子数（無吸収時光強度）Ｎ（ｔ１）〜Ｎ（ｔ３）を取得する（ステップＳ３４）。
【０１２７】
光路長取得部３０７が皮膚の各層の光路長を取得し、無吸収時光強度取得部３０８が検出光子数（無吸収時光強度）Ｎ（ｔ１）〜Ｎ（ｔ３）を取得すると、光吸収係数算出部３１０は、式（１３）に基づいて、積分区間算出部３０９が算出した、ある積分区間での皮膚の各層の光吸収係数μ１〜μ３を算出する（ステップＳ３５）。ここで、光吸収係数μ１は、表皮層の光吸収係数を示し、光吸収係数μ２は、真皮層の光吸収係数を示し、光吸収係数μ３は、皮下組織層の光吸収係数を示す。
【０１２８】
【数１３】
0005521199.tif 000014
【０１２９】
但し、ｌｎ（Ａ）はＡの自然対数を示し、Ｎ（ｔ）は特定波長λｋの短時間パルス光の時間分解波形のモデルの時刻ｔにおける光強度を示す。また、Ｉｉｎは、照射部３０４が照射した短時間パルス光の光強度を示す。また、Ｎｉｎは、シミュレーション部３０１が照射のシミュレーションを行った光子の個数を示す。
【０１３０】
光吸収係数算出部３１０がある積分区間での皮膚の各層の光吸収係数μ１〜μ３を算出すると、光吸収係数分布記憶部３１１は、光吸収係数算出部３１０が算出した、ある積分区間での皮膚の各層の光吸収係数μ１〜μ３を記憶する（ステップＳ３６）。
【０１３１】
光吸収係数算出部３１０がある積分区間での皮膚の各層の光吸収係数μ１〜μ３を算出すると、光吸収係数算出部３１０は、設定した積分区間での真皮層３３の光吸収係数を算出したか否かを判断する（ステップＳ３７）。
本実施形態では、皮膚の主成分を水、たんぱく質、脂質、グルコースの４種類として血糖値の測定を行うので、光吸収係数算出部３１０は、４種類の波長λ１〜λ４に対して光吸収係数μ１〜μ３を算出したか否かを判定する。ここで、波長λ１〜λ４は、シミュレーション部３０１が光路長分布及び時間分解波形を算出した複数の波長の中から選出する。
【０１３２】
ここで、光吸収係数算出部３１０が設定した積分区間での真皮層の光吸収係数μ１〜μ３に算出しなかった光吸収係数があると判断した場合（ステップＳ３７：ＮＯ）、再度、ある時刻における受光強度の取得（ステップＳ３２）に戻り、まだ算出していない真皮層の光吸収係数を算出し、再度、設定した積分区間での真皮層の光吸収係数の算出の可否の判断（ステップＳ３７）を行う。
一方、光吸収係数算出部３１０が設定した積分区間での真皮層の光吸収係数μ１〜μ３を算出したと判断した場合（ステップＳ３７：ＹＥＳ）、真皮層の光吸収係数分布から光吸収係数を取得する（ステップ３８）。
【０１３３】
光吸収係数取得部３１２は、皮膚の主成分の種類数に対応した波長数の光吸収係数を算出したか否かを判断する（ステップＳ２４）。
ここで、光吸収係数取得部３１２が皮膚の主成分の種類数に対応した波長数の光吸収係数を算出していないと判断した場合（ステップＳ２４：ＮＯ）、短時間パルス光の照射（ステップＳ２１）に戻り、まだ算出していない皮膚の主成分の種類数に対応した波長数の光吸収係数を算出し、再度、光吸収係数の算出の可否の判断（ステップＳ２４）を行う。
【０１３４】
一方、光吸収係数取得部３１２が皮膚の主成分の種類数に対応した波長数の光吸収係数を算出したと判断した場合（ステップＳ２４：ＹＥＳ）、濃度算出部３１３は真皮層に含まれるグルコースの濃度を算出する（ステップＳ２５）。
濃度算出部３１３は、真皮層におけるグルコースの濃度を、下記の式（１４）により算出する。
【０１３５】
【数１４】
0005521199.tif 000015
【０１３６】
但し、μ２（１）〜μ２（４）は、真皮層における波長λ１〜λ４の光吸収係数を示す。また、ｇ１〜ｇ４は、真皮層におけるそれぞれ皮膚の主成分である水、たんぱく質、脂質、グルコースのモル濃度を示す。また、ε１（１）〜ε１（４）は、波長λ１〜λ４に対する水のモル吸光係数を示し、ε２（１）〜ε２（４）は、波長λ１〜λ４に対するたんぱく質のモル吸光係数を示し、ε３（１）〜ε３（４）は、波長λ１〜λ４に対する脂質のモル吸光係数を示し、ε４（１）〜ε４（４）は、波長λ１〜λ４に対するグルコースのモル吸光係数を示す。
つまり、式（１４）のｇ４を算出することで、真皮層に含まれるグルコースのモル濃度を求めることができる。
【０１３７】
ここで、式（１４）によりグルコースのモル濃度を求めることができる理由を説明する。皮膚の散乱係数の波長依存性は小さいので、検出光子数Ｎ（ｔ）及び光路長Ｌｎ（ｔ）の波長に対する変化は無視することができる。また、ベア・ランベルト（Beer-Lambert）の法則により、吸光度＝モル吸光係数×モル濃度で表すことができる。これにより、２波長で得られた時間分解計測より、検出光子数Ｎ（ｔ）を消去することで、真皮層において得られる吸収係数差と皮膚を形成する各成分のモル吸光係数との関係式を示す式（１４）を導くことができる。
【０１３８】
上述の血糖値測定装置３００は、内部にコンピューターシステム（図示略）を有しており、上述した各ステップの処理動作は、プログラムの形式でコンピューター読み取り可能な記録媒体に記憶されている。そこで、このプログラムをコンピューターが読み出して実行することにより、上記の処理動作を行うことができる。
ここで、コンピューター読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリー等が挙げられる。
また、このコンピュータープログラムを通信回線によりコンピューターに配信し、この配信を受けたコンピューターが当該プログラムを実行するようにしてもよい。
【０１３９】
また、上記プログラムは、上記の機能の一部を実現するためのものであってもよい。
さらに、上述した機能をコンピューターシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であってもよい。
【０１４０】
以上説明したように、本実施形態によれば、積分区間算出部３０９により、光路長取得部３０７が取得した光路長分布のモデルの皮膚の各々の層の光路長と、無吸収時光強度取得部３０８が取得した短時間パルス光の時間分解波形のモデルの無吸収時光強度と、計測光強度取得部３０６が取得した受光部３０５が受光した光強度分布とに基づいて、前記光強度分布から真皮層の光強度に対応する領域の積分区間を算出するので、積分区間算出部３０９により算出された積分区間を基に、受光部３０５が受光した光の強度から前記積分区間に対応する時間帯の光の強度を取得することにより、真皮層からの光を他の層からの光と区別して測定することができ、真皮層からの光に対する他の層からの光の影響を低減することができる。したがって、真皮層におけるグルコースの濃度を精度良く測定することができ、その結果、真皮層におけるグルコースの濃度を、非侵襲的にかつ精度良く定量することができる。
また、積分区間を可変させることにより、真皮層におけるグルコースの濃度の測定精度を高めることができる。
【０１４１】
なお、本実施形態では、シミュレーション部３０１が、光吸収係数がゼロの皮膚モデルに対して光を照射するシミュレーションを行うこととしたが、シミュレーション部３０１が行った光吸収係数がゼロの皮膚モデルに対して光を照射するシミュレーションの結果を、光路長分布記憶部３０２及び時間分解波形記憶部３０３に記憶させておけば、シミュレーション部３０１を備えなくとも、本実施形態と同様の作用・効果を奏することができる。
【０１４２】
［第５の実施形態］
図１６及び図１７は、本発明の第５の実施形態の血糖値測定装置（濃度定量装置）が血糖値を測定する動作を示すフローチャートである。
本実施形態の血糖値測定装置は、第３の実施形態の血糖値測定装置３００と同一の構成であり、光路長取得部３０７、無吸収時光強度取得部３０８、計測光強度取得部３０６、光吸収係数算出部３１０の動作が異なる。
【０１４３】
次に、本実施形態の血糖値測定装置が血糖値を測定する動作について説明する。
まず、被測定者が血糖値測定装置３００を手首等の皮膚に当て、測定開始スイッチ（図示せず）の押下等により血糖値測定装置３００を動作させると、照射部３０４は、皮膚３１に対して波長λｋの短時間パルス光を照射する（ステップＳ４１）。
この波長λｋとしては、例えば、シミュレーション部３０１が光路長分布及び時間分解波形を算出した複数の波長の中の１つが好ましい。
【０１４４】
照射部３０４が波長λｋの短時間パルス光を照射すると、受光部３０５は、照射部３０４から照射され皮膚３１によって後方散乱された光を受光する（ステップＳ４２）。
このとき、受光部３０５は、照射開始からの単位時間毎（例えば、１ピコ秒毎の時刻ｔ１〜ｔｍ）の受光強度を、内部メモリー（図示せず）に記録しておく。
【０１４５】
次いで、積分区間を変化させて真皮層の光吸収係数を算出する（処理Ｂ：ステップＳ４３）。
このステップＳ４３は、図１７に示す動作により行う。
まず、積分区間算出部３０９により、（１）後方散乱した光を受光する受光部３０５の出力する光強度が計測光強度取得部３０６の最小検出感度を超えて検出された時刻から最小検出感度と等しい光強度で検出された時刻までの時間、（２）シミュレーション部３０１で得られる無吸収時光強度を記憶している時間分解波形記憶部３０３から取得した無吸収時光強度の時間特性、（３）皮膚表面に接する受光部３０５と照射部３０４との間隔、（４）シミュレーション部３０１に与える皮膚モデルのサイズ及び光学特性（散乱係数、吸収係数、非等方性パラメーター、または屈折率）を用いて、積分区間を算出する。より具体的には、積分区間の開始時刻、終了時刻、増分時間を算出する（ステップＳ５１）。
【０１４６】
次いで、受光部３０５が受光を完了すると、計測光強度取得部３０６は、受光部３０５の内部メモリーに記録されている受光強度から、ある時刻から時間τの間の受光強度の時間分布を取得する（ステップＳ５２）。
例えば、皮膚の３つの層について４種類の波長を用いて濃度測定を行う場合には、３つの異なる時間τ１〜τ３における受光強度の時間分布を取得する。
【０１４７】
計測光強度取得部３０６が、時間τの間の受光強度の時間分布を取得すると、光路長取得部３０７は、光路長分布記憶部３０２が記憶する波長λ１の光路長分布から、ある時刻から時間τの間の皮膚の各層の光路長、例えば、光路長Ｌ１（ｔ１）〜Ｌ１（ｔ３）、Ｌ２（ｔ１）〜Ｌ２（ｔ３）、Ｌ３（ｔ１）〜Ｌ３（ｔ３）を取得する（ステップＳ５３）。
また、計測光強度取得部３０６が、時間τの間の受光強度の時間分布を取得すると、無吸収時光強度取得部３０８は、時間分解波形記憶部３０３が記憶する波長λ１の時間分解波形から、ある時刻から時間τの間の無吸収時光強度、例えば、ある時刻から時間τの間における検出光子数（無吸収時光強度）Ｎ（ｔ１）〜Ｎ（ｔ３）を取得する（ステップＳ５４）。
【０１４８】
光路長取得部３０７が皮膚の各層の光路長を取得し、無吸収時光強度取得部３０８が検出光子数（無吸収時光強度）Ｎ（ｔ１）〜Ｎ（ｔ３）を取得すると、光吸収係数算出部３１０は、式（１５）に基づいて、積分区間算出部３０９が算出した、ある積分区間での皮膚の各層の光吸収係数μ１〜μ３を算出する（ステップＳ５５）。ここで、光吸収係数μ１は、表皮層の光吸収係数を示し、光吸収係数μ２は、真皮層の光吸収係数を示し、光吸収係数μ３は、皮下組織層の光吸収係数を示す。
【０１４９】
【数１５】
0005521199.tif 000016
【０１５０】
但し、ｌｎ（Ａ）はＡの自然対数を示す。また、Ｉ（ｔ）は、時刻ｔにおける受光部３０５の受光強度を示し、Ｉｉｎは、照射部３０４が照射した短時間パルス光の光強度を示す。また、Ｎ（ｔ）は、時間分解波形の時刻ｔにおける検出光子数を示し、Ｎｉｎは、シミュレーション部３０１が照射のシミュレーションを行った光子の個数を示す。また、Ｌ１（ｔ）〜Ｌ３（ｔ）は、時刻ｔにおける皮膚の各層の光路長を示す。
【０１５１】
光吸収係数算出部３１０がある積分区間での皮膚の各層の光吸収係数μ１〜μ３を算出すると、光吸収係数分布記憶部３１１は、光吸収係数算出部３１０が算出した、ある積分区間での皮膚の各層の光吸収係数μ１〜μ３を記憶する（ステップＳ５６）。
【０１５２】
光吸収係数算出部３１０がある積分区間での皮膚の各層の光吸収係数μ１〜μ３を算出すると、光吸収係数算出部３１０は、設定した積分区間での真皮層の光吸収係数を算出したか否かを判断する（ステップＳ５７）。
本実施形態では、皮膚の主成分を水、たんぱく質、脂質、グルコースの４種類として血糖値の測定を行うので、光吸収係数算出部３１０は、４種類の波長λ１〜λ４に対して光吸収係数μ１〜μ３を算出したか否かを判定する。ここで、波長λ１〜λ４は、シミュレーション部３０１が光路長分布及び時間分解波形を算出した複数の波長の中から選出する。
【０１５３】
ここで、光吸収係数算出部３１０が設定した積分区間での真皮層の光吸収係数μ１〜μ３に算出しなかった光吸収係数があると判断した場合（ステップＳ５７：ＮＯ）、再度、ある時刻における受光強度の取得（ステップＳ５２）に戻り、まだ算出していない真皮層の光吸収係数を算出し、再度、設定した積分区間での真皮層の光吸収係数の算出の可否の判断（ステップＳ５７）を行う。
一方、光吸収係数算出部３１０が設定した積分区間での真皮層の光吸収係数μ１〜μ３を算出したと判断した場合（ステップＳ５７：ＹＥＳ）、真皮層の光吸収係数分布から光吸収係数を取得する（ステップ５８）。
【０１５４】
光吸収係数取得部３１２は、皮膚の主成分の種類数に対応した波長数の光吸収係数を算出したか否かを判断する（ステップＳ４４）。
ここで、光吸収係数取得部３１２が皮膚の主成分の種類数に対応した波長数の光吸収係数を算出していないと判断した場合（ステップＳ４４：ＮＯ）、短時間パルス光の照射（ステップＳ４１）に戻り、まだ算出していない皮膚の主成分の種類数に対応した波長数の光吸収係数を算出し、再度、光吸収係数の算出の可否の判断（ステップＳ４４）を行う。
【０１５５】
一方、光吸収係数取得部３１２が皮膚の主成分の種類数に対応した波長数の光吸収係数を算出したと判断した場合（ステップＳ４４：ＹＥＳ）、濃度算出部３１３は、上記の式（１４）に基づいて、真皮層に含まれるグルコースの濃度を算出する（ステップＳ４５）。
【０１５６】
このように、本実施形態によれば、吸収係数μ１〜μ３を、時間τの間の光路長の積分値によって算出する。これにより、計測した受光強度Ｉ（ｔ）に含まれている誤差による吸収係数μ１〜μ３の算出結果に対する影響を少なくすることができる。
【０１５７】
［第６の実施形態］
図１８は、本発明の第６の実施形態の血糖値測定装置（濃度定量装置）の構成を示す概略ブロック図である。
本実施形態の血糖値測定装置４００が、第４の実施形態の血糖値測定装置３００と異なる点は、この血糖値測定装置４００が、照射部３０４と、受光部３０５と、コンピューターシステム４０１とを備え、コンピューターシステム４０１は記憶部４０２と、制御部（ＣＰＵ）４０３とを備え、記憶部４０２は、シミュレーション部３０１におけるシミュレーションの結果を記憶している光路長分布記憶部３０２及び時間分解波形記憶部３０３の機能を実行し、制御部（ＣＰＵ）４０３は、光路長取得部３０７、無吸収時光強度取得部３０８、計測光強度取得部３０６、積分区間算出部３０９、光吸収係数算出部３１０、光吸収係数分布記憶部３１１、光吸収係数取得部３１２及び濃度算出部３１３の機能を実行する点である。
【０１５８】
本実施形態の血糖値測定装置４００においても、第１の実施形態の血糖値測定装置３００と同様の作用・効果を奏することができる。
【０１５９】
［第７の実施形態］
図１９は、本発明の第７の実施形態の血糖値測定装置（濃度定量装置）の構成を示す概略ブロック図であり、本実施形態の血糖値測定装置５００が第４の実施形態の血糖値測定装置３００と異なる点は、時間分解波形記憶部３０３と積分区間算出部３０９との間に、無吸収時光強度取得部３０８と並行に吸収時光強度取得部（光強度モデル取得部）５０１を設けた点である。
【０１６０】
吸収時光強度取得部５０１は、シミュレーション部３０１で得られる無吸収時光強度と同様に得られる吸収時光強度を記憶している時間分解波形記憶部３０３から、短時間パルス光の時間分解波形のモデルの所定の時刻における吸収時光強度を取得する。ここでは、時間分解波形記憶部３０３からある時刻における吸収時の光強度を取得する。
【０１６１】
積分区間算出部３０９は、光路長取得部３０７が取得した光路長分布のモデルの皮膚の各々の層の光路長と、無吸収時光強度取得部３０８が取得した短時間パルス光の時間分解波形のモデルの無吸収時光強度と、吸収時光強度取得部５０１が取得した短時間パルス光の時間分解波形のモデルの吸収時光強度と、計測光強度取得部３０６が取得した受光部３０５が受光した光強度分布とに基づいて、前記光強度分布から任意の層の光強度に対応する領域の積分区間を算出する。
【０１６２】
この血糖値測定装置５００を用いて血糖値を測定する場合、積分区間を変化させて真皮層の光吸収係数を算出する。
この場合、積分区間算出部３０９は、（１）後方散乱した光を受光する受光部３０５の出力する光強度が計測光強度取得部３０６の最小検出感度を超えて検出された時刻から最小検出感度と等しい光強度で検出された時刻までの時間、（２）シミュレーション部３０１で得られる無吸収時光強度を記憶している時間分解波形記憶部３０３から取得した無吸収時光強度の時間特性、（３）シミュレーション部３０１で得られる吸収時光強度を記憶している時間分解波形記憶部３０３から取得した吸収時光強度の時間特性、（４）皮膚表面に接する受光部３０５と照射部３０４との間隔、（５）シミュレーション部３０１に与える皮膚モデルのサイズ及び光学特性（散乱係数、吸収係数、非等方性パラメーター、または屈折率）を用いて、積分区間を算出する。より具体的には、積分区間の開始時刻、終了時刻、増分時間を算出する。
【０１６３】
なお、この積分区間算出部３０９が算出した積分区間に基づいて、光吸収係数算出部３１０〜濃度算出部３１３が特定深さの層に含まれるグルコースの濃度を算出する方法及び手順については、第４の実施形態の血糖値測定装置３００と全く同様であるから、説明を省略する。
【０１６４】
本実施形態においても、第４の実施形態と同様に、真皮層からの光を他の層からの光と区別して測定することができ、真皮層からの光に対する他の層からの光の影響を低減することができる。したがって、真皮層におけるグルコースの濃度を精度良く測定することができ、その結果、真皮層におけるグルコースの濃度を、非侵襲的にかつ精度良く定量することができる。
また、積分区間を可変させることにより、真皮層におけるグルコースの濃度の測定精度を高めることができる。
【０１６５】
以上、本発明の各実施形態について、図面を参照して説明してきたが、具体的な構成は上述のものに限られることはなく、本発明の要旨を逸脱しない範囲内において様々な設計変更等が可能である。
例えば、上記の各実施形態では、濃度定量装置として血糖値測定装置を、観測対象として人の手のひらの皮膚を、目的成分としてグルコースを、パルス光として短時間パルス光を、それぞれ取ることで、皮膚の真皮層に含まれるグルコースの濃度を測定する場合について説明したが、これに限らず、濃度定量方法を、複数の光散乱媒質の層から形成される観測対象の任意の層における目的成分の濃度を定量する他の装置に用いてもよく、特定波長の短時間パルス光を、特定波長の連続光に替えてもよい。
例えば、携帯型の皮膚主成分の濃度測定装置に適用した場合、皮膚疾患の検査や診断や治療に有効利用することが可能である。
【符号の説明】
【０１６６】
３００…血糖値測定装置（濃度定量装置）、３０２…光路長分布記憶部、３０３…時間分解波形記憶部、３０４…照射部、３０５…受光部、３０６…計測光強度取得部（光強度取得部）、３０７…光路長取得部、３０８…無吸収時光強度取得部（光強度モデル取得部）、３０９…積分区間算出部、３１０…光吸収係数算出部、３１１…光吸収係数分布記憶部、３１２…光吸収係数取得部、３１３…濃度算出部、３１…皮膚（観測対象）、３３…真皮層（任意の層）、４００…血糖値測定装置（濃度定量装置）、５００…血糖値測定装置（濃度定量装置）、５０１…吸収時光強度取得部（光強度モデル取得部）、Ｓ２１〜Ｓ２５、Ｓ３１〜Ｓ３８、Ｓ４１〜Ｓ４５、Ｓ５１〜Ｓ５８…ステップ
【図１】
0005521199.tif 000020
【図２】
0005521199.tif 000021
【図３】
0005521199.tif 000022
【図４】
0005521199.tif 000023
【図５】
0005521199.tif 000024
【図６】
0005521199.tif 000025
【図７】
0005521199.tif 000026
【図８】
0005521199.tif 000027
【図９】
0005521199.tif 000028
【図１０】
0005521199.tif 000029
【図１１】
0005521199.tif 000030
【図１２】
0005521199.tif 000031
【図１３】
0005521199.tif 000032
【図１４】
0005521199.tif 000033
【図１５】
0005521199.tif 000034
【図１６】
0005521199.tif 000035
【図１７】
0005521199.tif 000036
【図１８】
0005521199.tif 000037
【図１９】
0005521199.tif 000038