(57)【要約】【課題】広い波長範囲の近赤外スペクトル測定を必要とせずに簡便に且つ良好な血糖値測定を行う。【解決手段】生体に近赤外光を照射して生体組織からの拡散反射光あるいは透過光を受光して得られた信号から生体組織中のグルコース濃度を測定するにあたり、グルコース成分の吸収ピーク波長である１６００±４０ｎｍの波長範囲から選択した第一の波長の近赤外光と、生体組織中の散乱係数変化の補正用として１４００ｎｍ±２０ｎｍの波長範囲から選択した第二の波長の２波長の近赤外光とを用いる【選択図】図１
2014018478.tif 000002
【特許請求の範囲】
【請求項１】
生体に近赤外光を照射して生体組織からの拡散反射光あるいは透過光を受光して得られた信号から生体組織中のグルコース濃度を測定するにあたり、グルコース成分の吸収ピーク波長である１６００ｎｍを中心とした±４０ｎｍの波長範囲から選択した第一の波長の近赤外光と、生体組織中の散乱係数変化の補正用として１４００ｎｍを中心とした±２０ｎｍの波長範囲から選択した第二の波長の２波長の近赤外光とを用いることを特徴とする血糖値測定方法。
【請求項２】
脂肪成分の吸収が存在する１６５０ｎｍから１８５０ｎｍの波長範囲から選択した第三の波長の近赤外光を加えた３波長を用いることを特徴とする請求項１記載の血糖値測定方法。
【請求項３】
脂肪成分の吸収ピークである１７２５ｎｍを中心とした±５０ｎｍの波長範囲から選択した第三の波長の近赤外光を加えた３波長を用いて生体組織中グルコース濃度を測定することを特徴とする請求項１または２記載の血糖値測定方法。
【請求項４】
水成分の影響を受ける１３５０ｎｍから１３９０ｎｍの波長範囲から選択した第四の波長の近赤外光を加えた４波長を用いて生体組織中グルコース濃度を測定することを特徴とする請求項１または２記載の血糖値測定方法。
【請求項５】
血糖値の算出用の演算式として、グルコースの吸収スペクトル吸収ピーク波長である１６００ｎｍを中心とした±４０ｎｍの波長範囲から選択した第一の波長の信号から、生体組織中の散乱係数変化の補正用としての１４００ｎｍを中心とした±２０ｎｍの波長範囲から選択する第二の波長の信号を減算し、得られた値に係数をかけるものを用いることを特徴とする請求項１記載の血糖値測定方法。
【請求項６】
血糖値を算出する演算式として、グルコースの吸収スペクトル吸収ピーク波長である１６００ｎｍを中心とした±４０ｎｍの波長範囲から選択した第一の波長の信号から、生体組織中の散乱係数変化の補正用としての１４００ｎｍを中心とした±２０ｎｍの波長範囲から選択する第二の波長信号を減算し、得られた結果に対して、１６５０ｎｍから１８５０ｎｍの波長範囲から選択した第三の波長の信号を用いて脂肪成分の吸収変化に起因する外乱影響の補正を行ない、この補正後の値に係数をかけるものを用いることを特徴とする請求項１または２または５記載の血糖値測定方法。
【請求項７】
血糖値を算出する演算式として、グルコースの吸収スペクトル吸収ピーク波長である１６００ｎｍを中心とした±４０ｎｍの波長範囲から選択した第一の波長の信号から、生体組織中の散乱係数変化の補正用としての１４００ｎｍを中心とした±２０ｎｍの波長範囲から選択する第二の波長信号を減算し、得られた結果に対して、脂肪の吸収スペクトルピークである１７２５ｎｍを中心とした±５０ｎｍの波長範囲から選択した第三の波長の信号を用いて脂肪成分の吸収変化に起因する外乱影響の補正を行ない、補正後の値に係数をかけるものを用いることを特徴とする請求項１または２または３または５または６記載の血糖値測定方法。
【請求項８】
第三の波長もしくは第四の波長の信号の値として、変化を鈍らせた値を用いることを特徴とする請求項２または３または４または６または７記載の血糖値測定方法。
【請求項９】
請求項１〜８のいずれか一項に記載の測定方法を実施する血糖値測定装置であって、生体に近赤外光を照射する発光手段と生体組織からの拡散反射光あるいは透過光を受光する受光手段とを０．３ｍｍ以上２ｍｍ以下の受発光間隔で配置したセンシング手段を備えるとともに、生体に常時接触させた上記センシング手段の受光手段から断続的に得られた吸光度信号を基に血糖値を演算する演算手段と、基準とする時点からの血糖値の相対変化を経時的に表示する表示手段とを備えていることを特徴とする血糖値測定装置。
【請求項１０】
前記センシング手段は生体への装着時に１０ｇ重／ｃｍ2以下の接触圧力で生体に接触するものであることを特徴とする請求項９記載の血糖値測定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、生体組織に照射した近赤外光の拡散反射光あるいは透過光から生体の血糖値（の代用特性としてのグルコース濃度）を測定する血糖値測定方法及び血糖値測定装置に関するものである。
【背景技術】
【０００２】
生体組織に近赤外光を照射し、生体組織内を拡散反射した光あるいは生体組織を透過した透過光を測定して得られる信号やスペクトルを基に生体組織の定性・定量分析を行う場合、生体内の種々の情報を非侵襲的に、試薬なしに、その場で、即時に得ることができる。このために医療分野における多くの用途で注目されており、血中酸素濃度測定については既にパルスオキシメーターとして広く利用されている。
【０００３】
そして血糖値測定への応用については、糖尿病患者の血糖値管理の他に、集中治療室（ＩＣＵ）や手術室で患者の血糖値を適切な範囲に管理する用途などの広い分野で有効性が検討されている。
【０００４】
そして、近赤外光を用いた生体組織の分光分析にあたっては、近赤外領域における各成分スペクトルがブロードで互いに重畳し合う特性があるため、他の成分や状態変化の影響を受けやすい。このために定量・定性分析に際しては多変量解析が用いられることが多い。
【０００５】
多変量解析に際しては数波長の信号を選択して重回帰分析を行なうこともあるが、数百ｎｍの広い波長範囲に及ぶ近赤外スペクトル信号を測定して、主成分回帰分析やＰＬＳ回帰分析のような多変量解析を行うことが多い。特に、生体中のグルコース濃度のように微小信号の変化を捉える必要がある場合、特定の波長を選択して数波長のみで定量化することは難しいと考えられており、このために近赤外スペクトルを用いた手法が多く検討されている（特開２００６−８７９１３号公報（特許文献１）参照）。
【０００６】
しかし、近赤外光の測定には十分な光強度を有する光源や分光手段、受光手段が必要であり、特に第１倍音およびその結合音が観測される１３００〜２５００ｎｍの波長領域においては適切な性能の光源、受光素子の入手が難しく、小型化が困難で、コスト的にも高価になってしまう欠点があった。しかも近赤外スペクトルを多変量解析することにより血糖値を算出するには複雑な演算処理が必要であり、高性能の演算用ソフトウェアや処理能力の高いＣＰＵや大容量メモリが必要なことも小型化や低コスト的が難しい要因となっている。
【０００７】
装置の小型化、低コスト化を行なう上では、広い波長範囲の近赤外スペクトル信号を基に測定するのではなく、数波長の信号から血糖値を測定する技術が有効であるが、それには適切な波長選択を行なう必要があり、このために以前より様々な波長の利用が検討されている。
【０００８】
たとえば、特表平７‐５０５２１５号公報（特許文献２）に開示されたグルコース濃度を測定するための方法及び装置においては、１５４７〜１５７７ｎｍの波長帯域内の第１の近赤外光と、１２９５〜１３０５ｎｍの波長帯域内の第２の近赤外光とを用いることが開示されている。
【０００９】
しかし、上記複数波長の近赤外光を用いても、血糖値（グルコース濃度）の測定において、十分な結果を得ることができない。
【先行技術文献】
【特許文献】
【００１０】
【特許文献１】特開２００６−８７９１３号公報
【特許文献２】特表平７‐５０５２１５号公報
【発明の概要】
【発明が解決しようとする課題】
【００１１】
本発明は上記問題に鑑みてなされたもので、広い波長範囲の近赤外スペクトル測定を必要とせず、簡便に且つ良好な血糖値測定を行うことができる血糖値測定方法及び血糖値測定装置を提供することを課題とする。
【課題を解決するための手段】
【００１２】
本発明にかかる血糖値測定方法は、生体に近赤外光を照射して生体組織からの拡散反射光あるいは透過光を受光して得られた信号から生体組織中のグルコース濃度を測定するにあたり、グルコース成分の吸収ピーク波長である１６００ｎｍを中心とした±４０ｎｍの波長範囲から選択した第一の波長の近赤外光と、生体組織中の散乱係数変化の補正用として１４００ｎｍを中心とした±２０ｎｍの波長範囲から選択した第二の波長の２波長の近赤外光とを用いることに特徴を有している。
【００１３】
脂肪成分の吸収が存在する１６５０ｎｍから１８５０ｎｍの波長範囲から選択した第三の波長の近赤外光を加えた３波長を用いることが好ましく、殊に、脂肪成分の吸収ピークである１７２５ｎｍを中心とした±５０ｎｍの波長範囲から選択した第三の波長の近赤外光を加えた３波長を用いて生体組織中グルコース濃度を測定することが好ましい。
【００１４】
水成分の影響を受ける１３５０ｎｍから１３９０ｎｍの波長範囲から選択した第四の波長の近赤外光を加えた４波長を用いて生体組織中グルコース濃度を測定することも好ましい。
【００１５】
このほか、生じる外乱によっては１４８０nmから１５２０ｎｍの波長範囲から選択した波長を第四の波長の近赤外光を加えた４波長を用いて生体組織中グルコース濃度を測定してもよい。もちろん、上記波長範囲から選択した５波長を用いて生体組織中グルコース濃度を測定してもよい。
【００１６】
血糖値の算出用の演算式としては、グルコースの吸収スペクトル吸収ピーク波長である１６００ｎｍを中心とした±４０ｎｍの波長範囲から選択した第一の波長の信号から、生体組織中の散乱係数変化の補正用としての１４００ｎｍを中心とした±２０ｎｍの波長範囲から選択する第二の波長の信号を減算し、得られた値に係数をかけるものを好適に用いることができる。
【００１７】
血糖値を算出する演算式として、グルコースの吸収スペクトル吸収ピーク波長である１６００ｎｍを中心とした±４０ｎｍの波長範囲から選択した第一の波長の信号から、生体組織中の散乱係数変化の補正用としての１４００ｎｍを中心とした±２０ｎｍの波長範囲から選択する第二の波長信号を減算し、得られた結果に対して、１６５０ｎｍから１８５０ｎｍの波長範囲から選択した第三の波長の信号を用いて脂肪成分の吸収変化に起因する外乱影響の補正を行ない、この補正後の値に係数をかけるものを用いてもよく、特に、血糖値を算出する演算式として、グルコースの吸収スペクトル吸収ピーク波長である１６００ｎｍを中心とした±４０ｎｍの波長範囲から選択した第一の波長の信号から、生体組織中の散乱係数変化の補正用としての１４００ｎｍを中心とした±２０ｎｍの波長範囲から選択する第二の波長信号を減算し、得られた結果に対して、脂肪の吸収スペクトルピークである１７２５ｎｍを中心とした±５０ｎｍの波長範囲から選択した第三の波長の信号を用いて脂肪成分の吸収変化に起因する外乱影響の補正を行ない、補正後の値に係数をかけるものを好適に用いることができる。
【００１８】
第三の波長もしくは第四の波長の信号の値として、変化を鈍らせた値を用いることも好ましい。ただし、１４８０nmから１５２０ｎｍの波長範囲から選択した波長を用いる場合、この波長は発汗等の急峻な状態変化に対応することが目的であるために、生データあるいはその変化を損なわない範囲の処理を施した値を用いることが望ましい。
【００１９】
そして本発明に係る血糖値測定装置は、上記測定方法を実施する血糖値測定装置であって、生体に近赤外光を照射する発光手段と生体組織からの拡散反射光あるいは透過光を受光する受光手段とを０．３ｍｍ以上２ｍｍ以下の受発光間隔で配置したセンシング手段を備えるとともに、生体に常時接触させた上記センシング手段の受光手段から断続的に得られた吸光度信号を基に血糖値を演算する演算手段と、基準とする時点からの血糖値の相対変化を経時的に表示する表示手段とを備えていることに特徴を有している。
【００２０】
前記センシング手段は生体への装着時に１０ｇ重／ｃｍ2以下の接触圧力で生体に接触するものが好適である。
【発明の効果】
【００２１】
本発明によれば、広い波長範囲の近赤外スペクトル測定を必要とせずに複数の波長の拡散反射光または透過光の測定のみで血糖値を演算予測することができるものであり、演算も簡便な演算でよく、近赤外光の受発光手段も簡便なものでよいこともあって、低コストで良好な相関係数を有する血糖値を予測することができる。
【図面の簡単な説明】
【００２２】
【図１】模擬生体における１６００ｎｍと１４００ｎｍの吸光度差とグルコース濃度との関係を示すグラフである。
【図２】同上の散乱係数の影響を示すグラフである。
【図３】同上の波長と相関係数との関係を示すグラフである。
【図４】人体における１６００ｎｍの波長の吸光度の時間変化を示すグラフである。
【図５】同上の人体における１６００ｎｍと１４００ｎｍの吸光度差の時間変化と血糖値の時間変化とを示すグラフである。
【図６】同上の吸光度差から導いた予測血糖値と実測血糖値の時間変化を示すグラフである。
【図７】他例における吸光度差から導いた予測血糖値と実測血糖値の時間変化を示すグラフである。
【図８】同上の差分吸光度と波長との相関を示すグラフである。
【図９】同上の補正を加えた後の予測血糖値と実測血糖値の時間変化を示すグラフである。
【図１０】更に他例にかかる補正を加えた後の予測血糖値と実測血糖値の時間変化を示すグラフである。
【図１１】別の例における吸光度差の時間変化と血糖値の時間変化とを示すグラフである
【図１２】同上の差分吸光度と波長との相関を示すグラフである。
【図１３】同上の補正を加えた後の予測血糖値と実測血糖値の時間変化を示すグラフである。
【図１４】さらに別の例における１６００ｎｍと１４００ｎｍの吸光度差の時間変化と血糖値の時間変化とを示すグラフである。
【図１５】同上の補正を加えた後の吸光度差の時間変化と血糖値の時間変化とを示すグラフである。
【図１６】同上の更なる補正を加えた後の予測血糖値と実測血糖値の時間変化を示すグラフである。
【図１７】(a)は本発明に係る装置の一例を示すブロック図、(b)は同上の測定プローブ先端面の正面図である。
【発明を実施するための形態】
【００２３】
近赤外光で生体組織中のグルコース濃度を測定することを困難にする要因として、生体での外乱がある。近赤外分光法では、通常、測定した近赤外スペクトル（説明変量）を多変量解析することにより外乱の影響を小さくし、目的変量（グルコース濃度変化）に相関する信号を抽出するが、グルコース濃度変化に付随する信号に対して、生体組織では個体差・季節差・日間差・測定部位差等のような大きな外乱が多数存在しており、精度を確保するために必要なＳＮ比が期待できないような場合は、多変量解析手法を用いても良好な関係式を導き出すことは難しい。
【００２４】
上記のような外乱の生体組織への影響を検討するに当り、その前段階として測定条件の制御が容易な模擬生体（ファントム）を用いて生体組織中のグルコース濃度と外乱の関係を明らかにして、その対応策を検討することで近赤外光によるグルコース検出技術の考え方と方向性の検討を行なった。
【００２５】
生体の光学ファントムとして、イントラリピッド（ダイズ油注射液、フレゼニウス  カービ製）を用いた模擬皮膚ファントムを作成し、散乱係数変化等の外乱を人為的に付与することでその光学特性を検討した。その手順を以下に説明する。
【００２６】
模擬皮膚ファントムは以下のように作成した。２０％濃度のイントラリピッドを蒸留水で希釈し、ヒト皮膚の散乱係数（波長１６００ ｎｍにおいて、等価散乱係数（μｓ’）≒ １．５ ｍｍ-1）とほぼ同等の散乱係数（１．３５ｍｍ-1，１．６５ ｍｍ-1の２水準）を有するイントラリピッド溶液を作成した。
【００２７】
次に上記２水準のイントラリピッド溶液にグルコースを溶解し、上記の散乱係数毎に５水準のグルコース濃度（１００ｍｇ／ｄｌ，２００ｍｇ／ｄｌ，３００ｍｇ／ｄｌ，４００ｍｇ／ｄｌ，８００ｍｇ／ｄｌ）を有するイントラリピッド・グルコース溶液を作成した。この１０種類のイントラリピッド・グルコース溶液をホットスターラーで攪拌・過熱し、一定温度（２９±０．１℃）に調整し、測定プローブ先端を溶液中に挿入した状態で近赤外スペクトル測定を行った。
【００２８】
この測定に使用した測定プローブは特許文献１に示されたものと同じで、図１７に示すように、ハロゲンランプ１から発光された近赤外光は熱遮蔽板２、ピンホール３、レンズ４、光ファイババンドル５を介して生体組織（模擬皮膚ファントム）６に入射される。光ファイババンドル５には測定用光ファイバ７の一端とリファレンス用光ファイバ８の一端が接続されている。測定用光ファイバ７の他端は測定用プローブ９に接続されており、リファレンス用光ファイバ８の他端はリファレンス用プローブ１０に接続されている。さらに、測定プローブ９およびリファレンスプローブ１０は光ファイバを介して測定側出射体１１，リファレンス側出射体１２にそれぞれ接続されている。
【００２９】
生体組織６の表面にセンシング手段である測定プローブ９の先端面を接触させて近赤外スペクトル測定を行う時、光源１から光ファイババンドル５に入射した近赤外光は、測定用光ファイバ７内を伝達し、図１７(b)に示すような測定用プローブ９の先端から同心円周上に配置された１２本の発光ファイバ２０より生体組織６の表面に照射される。生体組織６に照射されたこの測定光は生体組織内で拡散反射した後に、拡散反射光の一部が測定プローブ９の先端に配置されている受光ファイバ１９に受光される。受光された光はこの受光側光ファイバ１９を介して、測定側出射体１１から出射される。測定側出射体１１から出射された光は、レンズ１３を通して回折格子１４に入射し、分光された後、受光素子１５において検出される。
【００３０】
受光素子１５で検出された光信号はＡ／Ｄコンバーター１６でＡＤ変換された後、パーソナルコンピュータなどの演算装置１７に入力される。リファレンス測定はセラミック板など基準板１８を反射した光を測定し、これを基準光として行う。すなわち、光源１から光ファイババンドル５に入射した近赤外光はリファレンス用光ファイバ８を通して、リファレンス用プローブ１０の先端から基準板１８の表面に照射される。基準板に照射された光の反射光はリファレンス用プローブ１０の先端に配置された受光光ファイバ１９を介してリファンレス側出射体１２から出射される。上記の測定側出射体１１とレンズ１３の間、及びこのリファンレス側出射体１２とレンズ１３の間にはそれぞれシャッター２１が配置してあり、シャッター２２の開閉によって測定側出射体１１からの光とリファンレス側出射体１２からの光のいずれか一方が選択的に通過するようになっている。
【００３１】
測定プローブ９とリファレンスプローブ１０の端面は図１７(b)に示すように円上に配置された１２本の発光ファイバ２０と中心に配置された１本の受光ファイバ１９で構成されている。発光ファイバ２０と受光ファイバ１９の中心間距離Ｌは０．３ｍｍ以上２ｍｍ以下、好ましくは０．６５ ｍｍである。
【００３２】
この測定プローブを皮膚表面に接触させてスペクトル測定を行うと、入射光ファイバより照射された近赤外光は皮膚組織内を拡散反射し、入射された光の一部が検出用光ファイバに到達し、その光の伝播経路は“バナナ・シェイプ”と呼ばれる経路をとる。この時、発光ファイバ２０と受光ファイバ１９の中心間距離Ｌが０．３ｍｍ以上２ｍｍ以下、特に０．６５ ｍｍであると、表面より表皮、真皮、皮下組織の層状構造を有する皮膚組織における真皮部分のスペクトルを選択的に測定することができることになる。
【００３３】
測定結果をプロットしたものを図２に示す。図２の縦軸は、グルコース成分の特異吸収波長（１５６０〜１６４０ｎｍの波長範囲）の中の１６００ｎｍにおける一定温度（２９±０．１℃）下の吸光度（Ａ．Ｕ．）、横軸は模擬皮膚ファントム中のグルコース濃度（ｍｇ／ｄｌ）である。
【００３４】
図２から明らかなように、等価散乱係数１．３５ ｍｍ-1の測定値は上方の回帰直線付近に、等価散乱係数１．６５ ｍｍ-1の測定値は下方の回帰直線付近に存在する。これは１６００ｎｍにおける吸光度は、散乱係数依存が非常に大きいことを示しており、等価散乱係数が一定で外乱のない状態ではグルコース濃度と波長１６００ｎｍにおける吸光度とは正比例（リニア）の関係が存在するものの、模擬皮膚ファントム中の等価散乱係数が変化すると、グルコース濃度変化に対応する吸光度変化よりもずっと大きな吸光度変化が生ずることになる。つまりは１６００ｎｍの波長の吸光度変化のみでは生体中のグルコース濃度を測定することは実際上無理である。
【００３５】
このために、前記測定で得た１０個の測定スペクトル毎に、グルコースの特異吸収波長である１６００ｎｍの吸光度と、１３８０ｎｍ〜１８５０ｎｍ間の各波長（約２ｎｍ間隔）の吸光度とを取り出して、前者の吸光度と後者の吸光度との吸光度差を夫々算出して、各吸光度差とグルコース濃度との相関係数を求めた。図３は得られた１０個の吸光度差とそれぞれに対応するグルコース濃度の相関係数を波長毎に示したグラフである。得られた吸光度差とグルコース濃度の相関係数は１４００ｎｍに急峻な特異的ピークを有し、そこでの相関係数がほぼ１となっていることが分かる。
【００３６】
また、グルコース分子の特異吸収波長である１６００ｎｍと相関係数がピーク値を有する１４００ｎｍの吸光度差をとり、グルコース濃度との関係をプロットしたものを図１に示す。図１から分かるように、１４００ｎｍにおける吸光度を引くことにより、図２に見られた等価散乱係数の違いによる吸光度の上下分布が解消され、等価散乱係数の値（１．３５ｍｍ-1、１．６５ ｍｍ-1）にかかわらず前記吸光度差とグルコース濃度がほぼ直線状に並ぶことが分かる。
【００３７】
この結果は模擬生体での実験結果ではあるが、生体皮膚に擬して散乱係数を詳細に設計した模擬皮膚ファントムから得られた結果であるために、生体組織内においても散乱係数変化に起因する外乱を解消できることを期待できる。生体組織内での散乱係数変化は水分量変化等の要因により生ずる大きな外乱要因の一つと考えられ、皮膚組織内で生じる散乱係数変化を解消できれば、得られた知見は生体組織中のグルコース濃度測定を行う上で非常に有益なものと考えられる。
【００３８】
また、前記の実験については温度一定の状態で測定を行ったが、同様な実験を散乱係数一定で温度を変化（２７℃と２９℃）させた模擬生体で行なったところ、１４６０ｎｍで温度変化に起因する外乱が散乱係数変化と同様に解消できることが分かった。しかしながら、生体においては皮膚組織の温度管理は比較的容易であることから、以下の生体での検討については皮膚温度を一定とした条件で行なった。
【００３９】
上記実験をふまえて、模擬生体で得られた知見が生体皮膚でも再現できるか否かを下記の要領で確かめた。
【００４０】
被験者（５０歳代男性）に対して経口での糖負荷実験を実施し、その血糖値変化と１６００ｎｍと１４００ｎｍの吸光度差の関係を調べた。糖負荷実験は以下の手順で行なった。
【００４１】
座位の被験者に対して前述の測定プローブを左前腕内側部分に接触圧力が１０ｇ重／ｃｍ2以下で軽く接触する程度に貼付し、５分間隔で近赤外吸光度測定を繰り返した。また、比較データとして近赤外光による血糖値測定のタイミングに合わせ、１５分間隔で採血による血糖値を測定した。また、近赤外スペクトル測定に対応する血糖値測定を行わない１５分間隔の間の５分目、１０分目の２点は直線補間により推定した。
【００４２】
近赤外スペクトル測定開始直後と２時間後の２回、経口による糖負荷を行い、被験者の血糖値を変動させた。糖負荷には液体タイプの液体飲料２００ｍｌ（カロリーメイト缶タイプ  大塚製薬）を用いた。測定は測定開始後２回目の糖負荷後、血糖値が通常の１００ｍｇ／ｄｌ以下で安定するまでの４時間程度実施した。
【００４３】
得られた結果を図４〜図６に示す。図４はグルコース成分の特異吸収波長である１６００ｎｍの吸光度変化を経時的に示したもので、横軸は時間（時：分）、縦軸は吸光度（Ａ．Ｕ．）である。この実験における１６００ｎｍの吸光度は時間と共に単調に増加しており、１６００ｎｍの吸光度と血糖値の相関係数は‐０．２８と低くほぼ無相関で、２つのピークを有する実際の血糖値変化の特徴を全く反映していない。
【００４４】
図５はグルコース成分の特異吸収波長である１６００ｎｍと１４００ｎｍの吸光度差の変化を経時的に示したもので、横軸は時間（時：分）、縦軸は差分吸光度（Ａ．Ｕ．）である。上に示した無処理の１６００ｎｍの吸光度は時間と共に単調に増加した事と比較し、１６００ｎｍと１４００ｎｍとの吸光度差は２つのピークを有する実際の血糖値変化の特徴を強く反映しており、吸光度差と血糖値の相関係数は０．９２である。このことは１６００ｎｍの吸光度と１４００ｎｍの吸光度の差を取ることで、皮膚組織内に生じる散乱係数変化の影響を排除することができ、血糖値の変動に対応する変化を抽出することができたと考えられる。
【００４５】
図５に示した吸光度差を血糖値に換算するには、血糖値１ｍｇ／ｄｌの変動に相当する吸光度（Ａ．Ｕ．）で割り算すればよく、この事例においては０．００００８Ａ．Ｕ．／（ｍｇ／ｄｌ）を用いて血糖値に換算した。これにより相対的な変化ではあるが、血糖値の変動を推定することが可能となった。図６は測定開始から４５分後時に採血により測定した血糖値を基準に、その値からの相対変化として差分吸光度を血糖値に換算したものである。用いた換算係数（０．００００８Ａ．Ｕ．／（ｍｇ／ｄｌ））は皮膚組織を伝播する近赤外光の平均光路長に依存する値であり、測定部位の皮膚状態や測定光の中心波長、半値幅等により変動するもので、この値に限ったものではなく、被験者毎、装置毎に適宜較正する必要がある。ただし、被験者間の差については測定部位を選べば個体差が比較的小さいので代表的な値を固定値として代用することが可能である。
【００４６】
上記の実験から、模擬生体での知見と同様に、グルコース成分の特異吸収波長である１６００ｎｍの吸光度と１４００ｎｍの吸光度の差を取ることで、皮膚組織内で生じる散乱係数変化による吸光度変化をキャンセルして、精度良く血糖値が予測できることが生体においても確認できた。
【００４７】
２つの波長の吸光度のみから血糖値を定量できるために、血糖値算出のための演算負荷を小さくすることができる。また、発光素子として、ＬＥＤ（発光ダイオード）や半導体レーザーを用い、受光素子として常温駆動型のＩｎＧａＡｓ受光素子を利用することで、装置の小型化及び低コスト化が可能である。
【００４８】
また、採血による実測血糖値での校正が必要であるとはいえ、相対的変化としての血糖値の時間的変化は的確に捉えることができるために、図１７に示す演算装置１７におけるモニターに上記変化を表示する場合、集中治療室（ＩＣＵ）や手術室で患者の血糖値の観測を行う場合などにきわめて有用である。
【００４９】
次に、皮膚組織内の散乱係数変化以外の外乱への対応について記述する。前述の実験と同様の装置及び同様の方法で、測定開始から２時間目に経口による糖負荷を一度行うことで被験者（５０歳代男性）の血糖値を変動させた。測定プローブは左前腕内側部分に接触圧力１０ｇ重／ｃｍ2以下で軽く接触する程度に貼付した。
【００５０】
図７は実験で得られたグルコース成分の特異吸収波長である１６００ｎｍの吸光度と１４００ｎｍの吸光度の差の経時変化を示している。グルコースの特異吸収波長である１６００ｎｍの吸光度と１４００ｎｍの吸光度の差をとることで、吸光度差が血糖値の変化に幾分追随しているとは言え、時間経過とともに予測値の誤差が正方向に増加する右上りの傾向が出現しており、結果として相関を低下させている。ここでの両者の相関係数は０．６である。
【００５１】
経時的に値が上昇傾向を有する右上りとなる理由は、１４００ｎｍの吸光度の差を取ることで解消できる皮膚組織内の散乱係数変化以外の外乱が発生しているためと考えられる。そこで測定した近赤外スペクトルの変化を精査したところ、スペクトルの経時的変化として脂肪成分の吸収ピークの増加が特徴的であることが分かった。図８は測定開始時に測定した被験者の皮膚近赤外スペクトルを基準とし、測定終了まで５分毎に測定したスペクトルの変化をプロットしたもので、横軸が波長（ｎｍ）、縦軸が測定スペクトルと開始時のスペクトルの吸光度差（Ａ．Ｕ．）である。
【００５２】
図８に示されるように、時間経過と共に脂肪成分の特異吸収波長である１７２５ｎｍ付近が特徴的に増加している。この脂肪吸収の変化と考えられる現象は皮膚組織内の脂肪濃度の増加とは考えにくいために、皮膚組織内を拡散する近赤外光の光路変化と考えることが合理的である。つまり、測定プローブと皮膚が接触するとその接触面は初期状態においては皮溝、皮丘により粗い状態にあるが、時間経過と共に接触により皮膚表面が平滑化され、皮膚表面における散乱係数が低下したためと考えられる。平滑化されたならば、近赤外光が皮膚組織内に入りやすくなり、脂肪を中心とした皮下組織まで到達する近赤外光が増加し、あたかも脂肪濃度が増加したようなスペクトル変化が生じることとなる。なお、この現象は測定プローブと皮膚組織との接触圧力が大きいと顕著となることから、接触圧力は１０ｇ重／ｃｍ2以下とすることが好ましい。
【００５３】
近赤外領域で観測される吸収スペクトルは、波長２．５μm以上の赤外領域における急峻な吸収ピークとは異なり、裾野が広がったなだらかな吸収ピークを持つ、そのため脂肪成分の吸光度が増加すると１７２５ｎｍの脂肪吸収のピーク値だけでなく、その周辺の波長まで吸光度の増減が影響する。つまり、グルコース濃度の増減を検知している１６００ｎｍにおいても脂肪ピーク成長により見かけの吸光度増加が起こり、上記のような右上りの差分値となると考えられる。したがって、上記の差分値から脂肪成長分に相当する分を除去すれば精度よい推定が可能となる。
【００５４】
脂肪成分成長の影響を除去する為に、ここでは１７２５ｎｍの脂肪ピーク波長の吸光度を用いた。この時、脂肪ピーク波長においてもグルコース変動の影響を受ける為、脂肪ピーク成長の変化を鈍らせたトレンド的な値に加工して脂肪成長の指標にする方が精度よい推定が可能となる。脂肪ピーク成長の変化を鈍らせるには、５分毎に測定する１７２５ｎｍ吸光度の数点の移動平均値を利用したり、数点の吸光度に対する回帰式（１次回帰から数次回帰まで）を利用したりすることが考えられるが、ここでは測定開始時の１７２５ｎｍの吸光度と１４００ｎｍの吸光度差を基準とし、その後５分毎に測定した１７２５ｎｍの吸光度と１４００ｎｍの吸光度との差の変化分の総和を測定回数で割った値を用いた。また、ここでも０．００００８Ａ．Ｕ．／（ｍｇ／ｄｌ）を用いて血糖値に換算した。下の式は血糖値定量に用いた換算式である
血糖値（ｍｇ／ｄｌ）＝（（１６００ｎｍ吸光度‐１４００ｎｍ吸光度）−０．５×（鈍らせた１７２５ｎｍ吸光度と１４００ｎｍの吸光度差の変化））／０．００００８＋定数
結果を図９に示す。なお、本例においては測定開始後１５分の実測血糖値に近赤外光による予測血糖値が一致するように定数項を定めている。脂肪成分の吸収成長を除去することで、図から明らかなように、予測値に右上り傾向が生じる事象が解消され、相関係数０．９５の精度のよい血糖値予測が可能となった。
【００５５】
用いた換算係数（０．００００８Ａ．Ｕ．／（ｍｇ／ｄｌ））は、前述のように、皮膚組織を伝播する近赤外光の平均光路長に依存する値であり、測定部位の皮膚状態や測定光の中心波長、半値幅等により変動するものでこの値に限ったものではなく、被験者毎、装置毎に較正する必要がある。ただし、被験者間の差については測定部位を選べば個体差が比較的小さいので固定値として代用することは可能である。鈍らせた１７２５ｎｍ吸光度と１４００ｎｍの吸光度と差の変化に乗じている係数（−０．５）は実験的に求めた値であり、この値に限定されるものではない。
【００５６】
図７に見られるような吸光度差（グルコース成分の特異吸収波長である１６００ｎｍの吸光度と１４００ｎｍの吸光度の差）の右上りの経時変化は、脂肪成分の吸収ピークである１７２５ｎｍの吸光度を用いることなく、キャンセルすることもできる。
【００５７】
この場合の一例について説明すると、１６００ｎｍにおける吸光度変化は、グルコース濃度変化に起因するものとトレンド変化的なものが合わさったものと解釈することができる。従って、１６００ｎｍの吸光度を加工してトレンドラインを想定できれば１６００ｎｍと１４００ｎｍの２波長での測定が可能となる。
【００５８】
１６００ｎｍの吸光度のトレンドラインは、１６００ｎｍの吸光度を鈍らせることで求めることができる。この鈍らせることについては、５分毎に測定する数点の１６００ｎｍ吸光度を用いて移動平均値を用いたり、数点の吸光度に対する回帰式（１次回帰から数次回帰まで）を用いることができる。なお、ここでは測定開始時の１６００ｎｍの吸光度と１４００ｎｍの吸光度差を基準とし、その後５分毎に測定した１６００ｎｍの吸光度と１４００ｎｍの吸光度差の変化分の総和を測定回数で割った値を用いた。また、この事例においても０．００００８Ａ．Ｕ．／（ｍｇ／ｄｌ）を用いて血糖値に換算した。下の式は血糖値定量に用いた換算式である
血糖値（ｍｇ／ｄｌ）＝（（１６００ｎｍ吸光度‐１４００ｎｍ吸光度）−１．５×（鈍らせた１６００ｎｍ吸光度と１４００ｎｍの吸光度差の変化））／０．００００８＋定数
結果を図１０に示す。定数項は、測定開始後１５分の実測血糖値に、近赤外光による予測血糖値が一致するように定めた。図から明らかなように、脂肪成分の吸収成長の影響を排除して、予測値の右上り傾向を解消することができ、相関係数０．９４の精度のよい血糖値予測が可能となった。用いた換算係数（０．００００８Ａ．Ｕ．／（ｍｇ／ｄｌ））は皮膚組織を伝播する近赤外光の平均光路長に依存する値であり、測定部位の皮膚状態や測定光の中心波長、半値幅等により変動するものでこの値に限ったものではなく、被験者毎、装置毎に較正する必要がある。ただし、被験者間の差については測定部位を選べば個体差が比較的小さいので固定値として代用することは可能である。鈍らせた１６００ｎｍ吸光度と１４００ｎｍの吸光度差の変化にかかる係数（−１．５）は実験的に求めた値であり、これに限定されるものではない。
【００５９】
もっとも、トレンドラインを用いる校正は、測定開始直後から血糖値が単調に増加するような変化に対応することはできないが、実験開始時に血糖値が安定し、その後、健常人の血糖値変動のように１００ｍｇ／ｄｌ前後を基準に増減するような血糖値変動が起こることが予想される用途には適している。従って利用対象を限定する必要がある。
【００６０】
次に、皮膚組織内の散乱係数変化以外の外乱への対処についての他例を説明する。前述と同様の装置及び同様の手法で、測定開始から１時間目に一度、経口による糖負荷を行って被験者（５０歳代男性）の血糖値を変動させた。測定プローブは左前腕内側部分に接触圧力が１０ｇ重／ｃｍ2以下で軽く接触する程度に貼付した。
【００６１】
この実験で得られたグルコース成分の特異吸収波長である１５９８ｎｍの吸光度と１３９８ｎｍの吸光度の差の経時変化を図１１に示す。図に示されるように本実施形態においてはグルコースの特異吸収波長である１５９８ｎｍの吸光度と１３９８ｎｍの吸光度の差をとることで、吸光度差が血糖値の変化に幾分追随しているとは言え、時間経過とともに予測値の誤差が正方向に増加する右上りの傾向が出現していて、結果として相関を低下させている。ここでの両者の相関係数は０．６である。
【００６２】
値が経時的に上昇傾向を有する右上りとなる理由は、前述のように１４００ｎｍの吸光度の差を取ることで解消できる皮膚組織内の散乱係数変化以外の外乱の発生によると考えられる。このために測定開始時に測定した被験者の皮膚近赤外スペクトルを基準とし、測定終了まで５分毎に測定したスペクトルの変化をプロットすることでスペクトル変化の特徴を調べた。図１２に結果を示す。図において、横軸は波長（ｎｍ）、縦軸は測定スペクトルと開始時のスペクトルの吸光度差（Ａ．Ｕ．）である。
【００６３】
図から分かるように、時間経過と共に脂肪成分の特異吸収波長である１７２５ｎｍ付近が特徴的に増加しているほか、水成分の吸収ピークを示す１４５０ｎｍ周辺においても特徴的な変動が見られた。この水成分の吸収の変化と考えられる現象は、皮膚組織、特に皮膚表面の水分量変化の影響と考えられる。つまり、測定プローブと皮膚が接触するとその接触面に汗腺からの水分が蓄積し、角質水分量を上昇させることによる影響である。この実験においては、測定位置での角質水分量が実験開始時３４．７％であったものが終了時には３８．８％まで上昇していた。この表皮層の水による吸収が測定スペクトルの波長全体に影響を与えたものと思われる。
【００６４】
水成分の影響を取除くには１３５０から１５００ｎｍの間の波長、特に１４５０ｎｍの波長を用いることが有効であるが、１４５０ｎｍは吸光度が高いために得られる測定信号が微弱であり、しかも温度影響や生体成分との相互作用の影響を受けやすい。このために信号が大きく且つ水成分の変化を反映しやすい波長範囲として１３５０ｎｍから１３９０ｎｍを選択することが好ましい。このために、脂肪成分成長の影響を除去することを目的として、脂肪成分の吸収として１７２６ｎｍの吸光度を用い、加えて、水成分の影響を取除くために特に適している１３６０から１３８０ｎｍの間の１３７３ｎｍの吸光度を用いた。また、１３７３ｎｍの吸光度についても、直接用いるのではなく、トレンド的な値に加工した上で用いることが好ましいと考えられる。このために、水成分の影響を除去するための演算処理としては、測定開始時の１３７３ｎｍの吸光度と１３９８ｎｍの吸光度との差を基準とし、その後５分毎に測定した１３７３ｎｍの吸光度と１３９８ｎｍの吸光度との差の変化を測定し、その総和を測定回数で割った値、つまりは鈍らせた値を用いた。脂肪成分の吸光度を鈍らせることについては前述の例と同じとした。また、吸光度は０．００００６Ａ．Ｕ．／（ｍｇ／ｄｌ）を用いて血糖値に換算した。下の式は血糖値定量に用いた換算式である。
【００６５】
血糖値（ｍｇ／ｄｌ）＝｛（１５９８ｎｍ吸光度−１３９８ｎｍ吸光度）−０．７×（鈍らせた１７２６ｎｍ吸光度と１３９８ｎｍの吸光度差の変化）＋０．８×（鈍らせた１３７３ｎｍ吸光度と１３９８ｎｍの吸光度差の変化）｝／０．００００６＋定数
結果を図１３に示す。なお、測定開始直後の実測血糖値を近赤外光による予測血糖値と一致するように定数項を定めている。用いた換算係数（０．００００６Ａ．Ｕ．／（ｍｇ／ｄｌ））は皮膚組織を伝播する近赤外光の平均光路長に依存する値であり、測定部位の皮膚状態や測定光の中心波長、半値幅等により変動するものでこの値に限ったものではなく、被験者毎、装置毎に較正する必要がある。ただし、被験者間の差については測定部位を選べば個体差が比較的小さいので固定値として代用することは可能である。鈍らせた１７２５ｎｍ吸光度と１４００ｎｍの吸光度との差の変化等にかかる係数は実験的に求めた値であり、これに限定されるものではない。脂肪成分の吸収成長に加え水成分の影響を除去することで、予測値の右上り傾向が解消され、相関係数０．９６の精度のよい血糖値予測が可能となった。
【００６６】
上述の脂肪成分や水成分の影響の排除に際して用いた「鈍らせ手法」は、本発明においてなんら限定するものではなく、他の手法を用いてもよいのはもちろんである。
【００６７】
次に、皮膚組織内の散乱係数変化以外の外乱への対処についての他例を説明する。前述と同様の装置及び同様の手法で、測定開始から１時間１５分目に一度、経口による糖負荷を行って被験者（５０歳代男性）の血糖値を変動させた。測定プローブは左前腕内側部分に接触圧力が１０ｇ重／ｃｍ2以下で軽く接触する程度に貼付した。
【００６８】
この実験で得られたグルコース成分の特異吸収波長である１５９８ｎｍの吸光度と１３９８ｎｍの吸光度の差の経時変化を図１４に示す。図に示されるように本実施形態においてはグルコースの特異吸収波長である１５９８ｎｍの吸光度と１３９８ｎｍの吸光度の差をとることで、吸光度差が血糖値の変化に幾分追随しているが、１３時３５分と４０分の間に発汗が原因と考えられる差分吸光度の急激な低下が観察される。この事例における両者の相関係数は０．７９である。
【００６９】
上記の差分吸光度の急激な低下は、実験を行った室内温度が２８℃程度で蒸し暑く、糖負荷後に発汗が生じたことが原因と考えられる。発汗があると、光ファイバ端面が接触する部分が水分の影響を受けて差分吸光度が瞬時に変化するためである。このような変化も水変化の代用指標として利用できる波長を用いることで影響を軽減することができる。発汗影響の除去については、１４７０ｎｍから１５２０ｎｍの波長範囲から選択した波長が有効であり、ここでは１５０４ｎｍを用いた。図１４に示す変化から以下の式
吸光度差＝（１５９８ｎｍの吸光度と１３９８ｎｍの吸光度の差）−０．８×（１５０４ｎｍの吸光度と１３９８ｎｍの吸光度の差）
で発汗影響を除去した吸光度差を図１５に示す。
【００７０】
図１５に示すように１３時３５分と４０分の間に発汗が原因と考えられる差分吸光度の急激な低下が解消され、図７に見られるような時間経過とともに予測値の誤差が正方向に増加する右上りの傾向を示す関係にすることができた。
【００７１】
なお、発汗等の急激な変化に対応するために、この例においては鈍らせた値を用いていない。スムージング的なデータ加工を行なう場合は、発汗による急峻な差分吸光度変化を損なわない範囲で行う必要がある。
【００７２】
経時的に値が上昇傾向を有する右上りとなる理由は、前述のように脂肪成分の吸収ピークと影響と考えられることから、この影響を除去するために、１７２５ｎｍの脂肪ピーク波長の吸光度を用いた補正を行なった。この時、脂肪ピーク波長においてもグルコース変動の影響を受けることから、脂肪ピーク成長の変化を鈍らせたトレンド的な値に加工して脂肪成長の指標にする方が精度よい推定が可能となる。
【００７３】
脂肪ピーク成長の変化を鈍らせるには、５分毎に測定する１７２５ｎｍ吸光度の数点の移動平均値を利用したり、数点の吸光度に対する回帰式（１次回帰から数次回帰まで）を利用することが考えられるが、ここでは測定開始時の１７２６ｎｍの吸光度と１３９８ｎｍの吸光度差を基準とし、その後５分毎に測定した１７２６ｎｍの吸光度と１３９８ｎｍの吸光度との差の変化分の総和を測定回数で割った値を用いた。また、ここでも０．００００６Ａ．Ｕ．／（ｍｇ／ｄｌ）を用いて血糖値に換算した。下の式は血糖値定量に用いた換算式である
血糖値（ｍｇ／ｄｌ）＝（（（１５９８ｎｍの吸光度と１３９８ｎｍの吸光度の差）−０．８×（１５０４ｎｍの吸光度と１３９８ｎｍの吸光度の差））−０．４５×（鈍らせた１７２６ｎｍ吸光度と１３９８ｎｍの吸光度差の変化））／０．００００６＋定数
結果を図１６に示す。なお、本例においては測定開始後４０分の実測血糖値に近赤外光による予測血糖値が一致するように定数項を定めている。上記手法により、相関係数０．９１の精度のよい血糖値予測が可能となった。
【００７４】
ところで、グルコース分子の特異吸収波長である１６００ｎｍと相関係数がピーク値を有する１４００ｎｍの吸光度を求めることについては、前述の受発光間隔が２ｍｍ以下となっている測定プローブを好適に用いることができるが、このほか、近赤外光を用いて皮膚組織を測定する系であればどのようなものであってもよい。
【００７５】
たとえば共焦点法と呼ばれる手法を用いてもよい。この共焦点法では、光源に必要な波長の近赤外ＶＣＳＥＬ（垂直共振器面発光レーザー）を光源とし、ＰＤ（フォトダイオード）を光検出器として、奥行き分解能を有する共焦点光学系を作成することから、焦点位置を調整することで皮膚表面近傍の任意位置からの反射光の強度を検出することが可能である。血糖値を測定する場合は、波長１６００±４０ｎｍ、１４００±２０ｎｍ（さらには１７２５ｎｍ±５０ｎｍ、１３７０±１０ｎｍ）の近赤外ＶＣＳＥＬを用い、各波長から得られた信号を上記の各例で示したアルゴリズムを用いて血糖値予測を行なう。
【００７６】
このほか、光干渉断層法（Optical Coherence Tomography;OCT）を用いても可能である。光干渉断層法は光の干渉性（コヒーレンス）を利用して、物体内部の様子を撮像する技術であるが、皮膚表面近傍の一定深さの情報を選択的に取り出すことが可能である。この場合も、波長１６００±４０ｎｍ、１４００±２０ｎｍ（さらには１７２５ｎｍ±５０ｎｍ、１３７０±１０ｎｍ）の低コヒーレンス光源を用いて各波長から得られた信号を上記の各例と同様のアルゴリズムを利用して血糖値予測を行なえばよい。
【図１】
2014018478.tif 000003
【図２】
2014018478.tif 000004
【図３】
2014018478.tif 000005
【図４】
2014018478.tif 000006
【図５】
2014018478.tif 000007
【図６】
2014018478.tif 000008
【図７】
2014018478.tif 000009
【図８】
2014018478.tif 000010
【図９】
2014018478.tif 000011
【図１０】
2014018478.tif 000012
【図１１】
2014018478.tif 000013
【図１２】
2014018478.tif 000014
【図１３】
2014018478.tif 000015
【図１４】
2014018478.tif 000016
【図１５】
2014018478.tif 000017
【図１６】
2014018478.tif 000018
【図１７】
2014018478.tif 000019