(57)【要約】  基準波長を用いて成分濃度を測定するための方法は、被測定成分濃度の変化に対して光の強度がインセンシティブである波長を被測定成分の基準波長として確定するステップと、基準波長及びその他の測定波長におけるスペクトルを同時に測定して取得するステップと、基準波長におけるスペクトルを内部基準とし、その他の測定波長におけるスペクトルに対してデータ処理を行うことによって、被測定成分の特異性情報を含む特徴スペクトルを求めるステップと、特徴スペクトルと被測定成分濃度との数学補正モデルを作るステップと、数学補正モデルに基づいて被測定成分の濃度を取得するステップと、を含む。【選択図】  図４
2011522264.tif 000002
【特許請求の範囲】
【請求項１】
サンプル中の成分濃度を測定するための方法であって、
１）被測定成分濃度の変化に対して光強度がインセンシティブである波長を、被測定成分の基準波長と確定するステップと、
２）基準波長及びその他の測定波長におけるスペクトルを同時に測定して取得するステップと、
３）基準波長におけるスペクトルを内部基準として、その他の測定波長におけるスペクトルに対してデータ処理を行うことによって、被測定成分の特異性情報を含む特徴スペクトルを求めるステップと、
４）特徴スペクトルと被測定成分濃度との数学補正モデルを作るステップと、
５）数学補正モデルに基づいて被測定成分の濃度を取得するステップと、
を含むことを特徴とする方法。
【請求項２】
前記ステップ１）において、
吸収係数及び散乱係数と、被測定成分の濃度との関係に基づいて、各波長における出射光強度が被測定成分濃度の変化に従って変化する関係を解析し、光強度の変化率がゼロである波長を被測定成分の基準波長として確定することを特徴とする請求項１に記載の方法。
【請求項３】
純吸収サンプルに対して、
【数１】
2011522264.tif 000023
散乱サンプルに対して、
【数２】
2011522264.tif 000024
を用いて被測定成分の基準波長λｒを確定する（ただし、αｗ、αｇはそれぞれ水分子および被測定成分分子の吸光係数を示し、ｋは被測定成分と水との間の置換係数を示し、ｇは散乱サンプルの異方性因子を示し、ｎｓは散乱サンプルにおける散乱粒子の屈折率を示し、ｎｗは水の屈折率を示し、ｒは散乱粒子の半径を示し、ρｓは単位体積のサンプル中の散乱粒子の個数（ｎｍ−３）を示す）ことを特徴とする請求項２に記載の方法。
【請求項４】
前記ステップ１）において、
サンプル中の被測定成分の濃度を変更させ、変更後の濃度に対応するスペクトルを測定することによって、異なる濃度のサンプル間の差分スペクトルを求め、
差分スペクトルの曲線がゼロの値に交差する箇所に対応する波長を、被測定成分の基準波長として確定することを特徴とする請求項１に記載の方法。
【請求項５】
ブドウ糖水溶液において、ブドウ糖の基準波長が１５２５ｎｍであることを特徴とする請求項１〜４の何れか一項に記載の方法。
【請求項６】
光強度が被測定成分の濃度の変化に対して最もセンシティブである波長を、最適な測定波長として選択することを特徴とする請求項１に記載の方法。
【請求項７】
前記ステップ３）において、純吸収サンプルに対して、被測定成分の濃度変化の情報を含む特異性情報は、
【数３】
2011522264.tif 000025
であり、散乱サンプルに対して、被測定成分の濃度変化の情報を含む特異性情報は、
【数４】
2011522264.tif 000026
であるように、被測定成分の特異性情報を取得する（ただし、Ｉｓ(λ)、Ｉｂ(λ)はそれぞれ基準波長における被測定サンプルの出射光強度及びバックグラウンドサンプルの出射光強度を示し、Ｉｓ(λ)、Ｉｂ(λ)はそれぞれその他の測定波長における被測定サンプルの出射光強度及びバックグラウンドサンプルの出射光強度を示し、Ｉ０(ｔｂ，λｒ)、Ｉ０(ｔｂ，λ)は、それぞれ、バックグラウンドサンプルを測定するときの、基準波長及びその他の測定波長における入射光強度を示し、ＫΔｃはサンプル中の測定対象成分のΔｃの濃度変化による光路長の変化のパーセント係数を示している）ことを特徴とする請求項１に記載の方法。
【請求項８】
ステップ３）を実行した後に、
二つのサンプル間のスペクトル差信号を、温度変化による差信号と被測定成分の濃度変化による差信号との和で示し、そのうち、被測定成分の基準波長において被測定成分濃度による差信号がゼロであるステップと、
基準波長における光強度の差信号に基づいて、サンプルを測定する時の温度変化を算出してサンプルの実際の測定温度を取得するステップと、
その他の測定波長における光強度に対して温度の補正を行い、補正グループ温度に対応する測定スペクトルを取得するステップと、
を実行することによって、サンプルのスペクトルに対して温度の補正を行うことを特徴とする請求項１に記載の方法。
【請求項９】
部分最小二乗法または純信号解析法を用いて数学補正モデルを作ることを特徴とする請求項１に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、近赤外スペクトルを用いた物質の成分濃度の測定方法に関するものである。
【背景技術】
【０００２】
近赤外スペクトル法とは、主に近赤外吸収スペクトル及び化学的な計量法に基づいて、人体の成分（人体の血糖濃度等を含む）の測定を行うことに用いられる方法である。近赤外領域においては、体液や軟組織（ｓｏｆｔ  ｔｉｓｓｕｅ）が比較的透明なので、近赤外光の透過率が高い。また、化学的な計量法が進歩し、それが近赤外スペクトルの検出技術へ応用されるに伴い、近赤外スペクトル法を用いる定量解析の感度、精度及び信頼性は大きく向上されることになった。近赤外スペクトル法は、人体の成分を損傷させずに検出する技術において最も将来性のある技術であると注目されている。成功した先行例の一つとして、近赤外スペクトル法を用いて、人体を損傷させることなく血液内の酸素飽和度の検出に成功したことが挙げられる。
【０００３】
近赤外スペクトル（波長範囲が７８０ｎｍ〜２５００ｎｍ）は、主に、Ｃ−Ｈ、Ｏ−Ｈ、Ｎ−Ｈ等の化学結合（ｃｈｅｍｉｃａｌ  ｂｏｎｄ）を含む化合物が、中赤外領域での基本周波数振動のダブル周波数吸収及び複合周波数吸収によって得られたものである。Ｈの化学結合を含む有機物及びそれと結合する無機物等のサンプルは、その成分の含有量が変更されると、サンプルの近赤外スペクトルの特徴も変化することになる。スペクトル変化の特徴に基づいて、被測定成分の濃度変化を分析することができる。しかしながら、図１に示すように、人体の血液及び組織液においては、各種成分の近赤外スペクトルのオーバーラップが顕著に起こる。また、スペクトルを測定する過程に、設備機器におけるドリフトの影響、温度の変化及びバックグラウンドの干渉などの要因による影響を受ける。そのため、単一の波長データを用いて補正曲線を作成する場合は、大きな誤差が生じ、高い測定精度を得ることは困難である。
【０００４】
一般的に、近赤外スペクトル法は、多変量解析技術と結びつけて物質の成分に対する定量測定を行うことに用いられている。具体的な方法として、図２に示すように、先ず、測定しようとする成分、即ち測定対象成分の濃度が段階的に変化する一連の補正グループのサンプルを設け、各サンプルに対応する近赤外スペクトルを測定し；次に、多変量解析技術を用いて、サンプルのスペクトルマトリックス及び濃度マトリックスに対して回帰分析を行い、補正モデルを作ると共に補正モデルの回帰係数を求め；次に、未知の濃度である測定対象サンプルに対して、その近赤外スペクトルを測定し、補正モデルに基づいてスペクトルマトリックスを分解し、そして、補正モデルの回帰係数を用いて逆計算を行って測定対象成分の濃度を得る方法である。なお、部分最小二乗法（ＰＬＳ：Ｐａｒｔｉａｌ  Ｌｅａｓｔ  Ｓｑｕａｒｅ）は、現在、近赤外スペクトル解析において最も多く使用されかつ最も効果のある多変量モデル作成法である。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、近赤外スペクトル法は、人体の成分を損傷せずに測定を行う技術によって実現できる最適な方法の一つであるが、その実現のための技術は非常に複雑である。ここで、人体を損傷させずに行う血糖濃度測定を例に挙げて説明する。先ず、組織自体の光学特性が非常に複雑であるので、近赤外光が体内で動態変化する人体組織を透過する際に非線形的な散乱が生じることになり、このため、血糖の信号を抽出することが非常に困難である。また、人の体内における血糖の含有量は低く、血糖濃度の変化によるスペクトル信号は非常に弱い。そして、ブドウ糖自体の吸収が、人体内におけるその他の成分の吸収と非常にオーバーラップする。人体内の水分、蛋白質、脂肪などの物質による強い吸収、体温の変化及びその他の生理条件の変化による影響は、干渉の主な要因となっている。近赤外スペクトル法を用いて人体の血糖濃度を測定するために、オーバーラップが複雑に生じていたバックグラウンドから、微弱な糖濃度の変化による特異性情報を抽出することが必要になる。また、情報の抽出過程で、人体の各種の生理的要因による制約及び影響を受けてしまう。
【課題を解決するための手段】
【０００６】
本発明は、近赤外スペクトル測定の過程において、被測定成分以外のその他の要因による干渉を排除し、成分濃度の測定精度及び感度を向上させるため、基準波長を利用した濃度の測定を実現する方法を提案する。被測定成分の濃度が異なるという条件下で得られたサンプルの近赤外吸収スペクトルの特徴に基づいて、当該成分の基準波長を決定し、さらに、基準波長における信号変化の特徴に基づいて、基準波長法を用いて全帯域のスペクトル信号に対してデータ処理を行うことにより、被測定成分の特徴スペクトル信号を求める。そして、この特徴スペクトル信号のマトリックスと被測定成分の濃度マトリックスとの間の数学補正モデルをたてて、被測定成分の濃度を測定する。
【図面の簡単な説明】
【０００７】
【図１】図１は、血液成分の近赤外吸収スペクトルを示す図である。
【図２】図２は、近赤外スペクトル法による定量解析の原理を示すブロック図である。
【図３】図３は、ブドウ糖と水との吸光係数の比率を示す図である。
【図４】図４は、二成分サンプルの実験におけるブドウ糖の基準波長特性を示す図である。
【図５】図５は、三成分サンプルの実験におけるブドウ糖の基準波長特性を示す図である。
【図６】図６は、四成分サンプルの実験におけるブドウ糖の基準波長特性を示す図である。
【図７】図７は、ブドウ糖水溶液の実験における基準波長の温度特性を示す図である。
【図８】図８は、異なる温度下でのブドウ糖水溶液の置換吸光係数の曲線を示す図である。
【図９】図９は、温度補正の前における、異なる温度でのブドウ糖水溶液の吸光度差の値の信号曲線（図において、Ｔｇはブドウ糖水溶液の測定時のサンプル温度を示し、Ｔｗは純水バックグラウンドの測定時のサンプル温度である）を示す図である。
【図１０】図１０は、温度補正の後における、異なる温度でのブドウ糖の特徴的な信号曲線（図において、Ｔｇはブドウ糖水溶液の測定時のサンプル温度を示し、Ｔｗは純水バックグラウンドの測定時のサンプル温度である）を示す図である。
【発明を実施するための形態】
【０００８】
（被測定成分の基準波長の取得）
本発明において、まず、基準波長の存在特性を取得して証明することと、異なる被測定成分に対する基準波長の値を取得することとを解決する必要がある。そして、これを基に、相関するデータの前処理モデル（Ｐｒｅｔｒｅａｔｍｅｎｔ  Ｍｏｄｅｌ）を作ることによって、濃度信号を有効的に抽出することが実現される。本発明は、理論推理により基準波長を証明する方法と、実験により基準波長を取得する方法とを提案した。
【０００９】
理論推理によりその方法を導く際、被測定成分と水またはその他の溶剤との体積置換原理に基づいて、当該被測定成分の濃度変化がサンプル吸収係数及び散乱係数に与える影響を算出し、その後、スペクトルが成分濃度に従う変化関係を推理し又は擬似算出し、成分濃度の変化に従って光強度がゼロになった波長を、当該被測定成分の基準波長として定義する。本発明においては、ブドウ糖の測定を例として説明し、基準波長の存在を理論上から証明する。
【００１０】
先ず、純吸収サンプルにおける基準波長の存在特性について解析する。ブドウ糖水溶液のサンプルについて、ブドウ糖が加えられた後、ブドウ糖分子の濃度が増大するにつれ、ブドウ糖の固有吸収は増大していく。他方、ブドウ糖分子は、一部の水分子と置換することによって、水分子の濃度が小さくなるにつれ、水の固有吸収が小さくなっていく。そのため、ブドウ糖水溶液のサンプル溶液における吸光度と、純水のバックグラウンド溶液における吸光度が、それぞれ、下記の式（１）と式（２）で示されることができる。すなわち、
【数１】
2011522264.tif 000003
ここで、Ａｂ、Ａｓは、それぞれ、純水のバックグラウンド溶液の吸光度とブドウ糖水溶液の吸光度を示し、Ｃｗは、純水の質量濃度を示し、ΔＣｇは、サンプル中のブドウ糖の質量濃度変化を示し、ΔＣｗは、サンプル中のブドウ糖の濃度変化がΔＣｇであるに応じて水の質量濃度の変化を示す。また、αｗ及びαｇは、それぞれ水分子とブドウ糖分子の吸光係数を示す。図１に、その波長特性を示す。
【００１１】
サンプル及びバックグラウンドの吸光度に対して差分処理を行い、サンプルとバックグラウンドの測定条件が一致した場合、ブドウ糖の濃度変化による特徴信号ＡＧは、下記の式（３）で示される。
【数２】
2011522264.tif 000004
ブドウ糖の基準波長において、サンプルの吸光度はブドウ糖の濃度の変化に従って変化せず、即ち、ＡＧ＝０である。そのため、式（３）から、下記の式（４）が得られる。
【数３】
2011522264.tif 000005
【００１２】
既知の研究成果によれば、温度が２０℃である場合、純水のモル濃度はＣｗ０＝５５．４Ｍである。そして、ブドウ糖濃度が１ｍＭ増加するにつれて、蒸留水分子のモル濃度は０．０１１１％減少する。従って、式（４）に基づいてさらに下記の式（５）が導き出される。
【数４】
2011522264.tif 000006
即ち、ブドウ糖分子と水分子の吸光係数の比率が０．６２１７である場合、対応する波長はブドウ糖測定の基準波長となる。図４から分かる通り、波長１５２５ｎｍにおいて、ブドウ糖分子と水分子の吸光係数の比率が、およそ０．６２１７である。従って、ブドウ糖を測定する際、基準波長は波長１５２５ｎｍである。
【００１３】
また、人体組織等の複雑に散乱するサンプルについて、ブドウ糖の濃度変化は、サンプルの吸収係数（吸収係数は吸光係数と質量濃度との積である）及び散乱係数の変化を共に引き起こす。これら二つの総合的な作用により、出射光の強度が変化する。従って、基準波長は、吸収効果と散乱効果とが相互に相殺した波長である。即ち、
【数５】
2011522264.tif 000007
我々が前期に行った研究によれば、以下の式が導き出された。
【数６】
2011522264.tif 000008
ここで、Ｉ(λ)、Ｉ0(λ)はそれぞれ出射光強度及び入射光強度を示し、μα、μｓはそれぞれ散乱サンプルの吸収係数及び散乱係数を示し、ｇは散乱サンプルの異方性因子を示し、ｌはサンプル中における出射光の伝送の経路長を示す。また、ｎｓは散乱サンプルにおける散乱粒子の屈折率を示し、ｎｗは水の屈折率を示している。ｒは散乱粒子の半径を示し、ρｓは単位体積あたりのサンプル中における散乱粒子の個数（ｎｍ−３）を示す。
【００１４】
式（６）〜（１０）によれば、散乱サンプル中において、基準波長λｒの存在条件は、下記の式（１１）で示される。
【数７】
2011522264.tif 000009
【００１５】
上述した理論推理法において、サンプル中における各組成成分のモル吸光係数及び各成分間の置換関係式を確実に知る必要がある。より複雑なサンプルに対して、本発明は、実験に基づいて有効的に基準波長を取得する方法を提供する。
【００１６】
実験による確定法は具体的に以下のステップを含む。即ち、（１）被測定成分の濃度が異なる条件下でのサンプルのスペクトルIs1(λ)、Is2(λ)、Is3(λ)、・・・と、被測定成分の濃度がゼロである場合のバックグラウンドサンプルのスペクトルIｂ1(λ)）、Iｂ2(λ)、Iｂ3(λ)、・・・とを測定するステップ、（２）各被測定サンプルと、対応するバックグラウンドサンプルの光強度とに対して差分処理を行うことによって、濃度の異なるサンプルの差分スペクトル、即ち被測定成分の濃度情報が含まれた特徴スペクトル信号Ａ濃度１(λ)、Ａ濃度２(λ)、Ａ濃度３(λ)、・・・を取得するステップ；（３）被測定成分の差分スペクトルの特徴を解析し、各差分スペクトル曲線がゼロの値と交差する波長を被測定成分の基準波長とするステップとが含まれることを特徴とする。
【００１７】
本発明は、二成分実験（ブドウ糖＋水）、三成分実験（アルブミン＋ブドウ糖＋水）と、四成分実験（ヘモグロビン＋アルブミン＋ブドウ糖＋水）などの一連のサンプルをそれぞれ用いて、簡単なものから複雑なものへ実験を行うことによって、基準波長の存在特性を検証した。実験結果は、それぞれ図４、図５及び図６に示されている。各回の実験から得られた差分吸光度の曲線は、波長１５２５ｎｍにおいてゼロの値で交差している。即ち、純吸収水溶液サンプル中におけるブドウ糖の基準波長は１５２５ｎｍである。
【００１８】
また、本発明は、二成分ブドウ糖水溶液を用いて、実験により基準波長の位置の温度特性を解析した。ここで、それぞれ、３０度、３４度と３８℃のサンプル温度での差分吸光度の曲線特性を解析した。図７に示すように、異なる温度下での差分吸光度は、いずれも波長１５２５ｎｍの所でゼロの値に交差している。即ち、ブドウ糖の基準波長の位置は、温度の変化に従ってシフトせず、固定の１５２５ｎｍとなっている。
【００１９】
（基準波長法に基づいて被測定成分の特異性信号の抽出）
人体の成分に対して、微細な損傷かまたは損傷させずに検出を行う際、異なる時間における測定スペクトル間には、機器システムのドリフト、その他の人体成分の変化および人体の生理状態の変化などのバックグラウンドノイズが存在するので、被測定成分の測定精度に影響を与える。本発明は、基準波長法を用いて、スペクトルのデータに対して前処理を行うことにより、これらのノイズによる干渉を有効的に排除することができる。基準波長法を用いてデータ処理を行う方法の核心となる内容は、測定過程において、基準波長において得られる光強度値または吸光度値を変動パラメーターとし、このパラメーターの値が、サンプル中に被測定成分以外その他の原因によるノイズの干渉を反映したが被測定成分の濃度変化とは無関係であるということである。基準波長における信号に基づいて、その他の波長におけるスペクトル信号に対して差分処理を行うことにより、各種のノイズの干渉を排除することができ、被測定成分から信号をピックアップする特異性が増強される。
【００２０】
純吸収サンプル条件下での出射光の強度は、完全にランベルト‐ベール法則（Ｌａｍｂｅｒｔ‐Ｂｅｅｒ  ｌａｗ）に従うことになっている。ここで、被測定成分濃度の変化をΔｃとする時、光源のドリフト及び測定環境のランダムな干渉からの影響によって、サンプル測定時の入射光強度が、バックグラウンド物質の測定時の入射光強度に対し、一定の比例ずれηが存在し、且つ、バックグラウンドとサンプルのスペクトル測定時にその他の成分濃度の変化も存在している。よって、バックグラウンドおよびサンプルの受光強度は、それぞれ、以下の式（１２）、（１３）で示される。
【数８】
2011522264.tif 000010
【００２１】
式（１２）、(１３)において、Ｉs、Ｉｂは、それぞれ、サンプル及びバックグラウンドの受光強度を示しており、Ｉ0（ｔｂ，λ）は、バックグラウンド測定時の入射光強度であり、μａはバックグラウンド溶液の吸収係数を示したものであり、Δμa_Δｃは、被測定成分濃度変化がΔｃである時に引き起こされたサンプル吸収係数の変化を示したものであり、Δμａ_Δotherは、その他の干渉成分の変化によって引き起こされた吸収係数の変化を示したものであり、ｌは、測定した光路長を示している。
【００２２】
また、測定対象成分の濃度特異性信号は、下記の式（１４）で示すことができる。
【数９】
2011522264.tif 000011
なお、式（１４）において、ＡＧ(λ)＝Δμa_Δc（λ）・ｌ，ＡΔother(λ)＝Δμa_Δother(λ)・ｌである。
【００２３】
基準波長λrにおいて、受光強度は被測定成分の濃度変化に対してインセンシティブである。そのため、基準波長λrにおいて生じた光強度の変化が、完全にバックグラウンドノイズによるものだと考えられる。したがって、被測定成分の濃度変化による溶液の吸光度の変化は、ゼロとなる。即ち、
ＡＧ(λ)＝０         ・・・（１５）
式（１４）および式（１５）によって、基準波長におけるバックグラウンドノイズは、下記の式（１６）で示される。
【数１０】
2011522264.tif 000012
【００２４】
サンプルとバックグラウンドとの間で、その他の成分の濃度変化が比較的小さい時、ＡΔother(λ)≒ＡΔother(λr)に近似することができる。さらに、純吸収サンプルに対して、基準波長に基づいて被測定成分の特異性の信号を抽出する式は、下記の式（１７）で示される。
【数１１】
2011522264.tif 000013
【００２５】
人体の組織等の複雑な散乱サンプルに対しては、吸収効果および散乱効果による影響を同時に考慮しなければならない。補正されたランベルト・ベール法則に基づき、サンプルに入射される入射光の強度の変動を考慮した場合に、ある濃度におけるバックグラウンド及びサンプルからの出射光の強度は、波長λにおいてそれぞれ下記の式（１８）、（１９）で示される。
【数１２】
2011522264.tif 000014
ここで、μeffは、バックグラウンドサンプルの有効減衰係数を示しており、Δμeff-Δcは、被測定成分濃度のΔｃの増加によるサンプルの有効減衰係数の変化を示しており、leffは、光が散乱サンプルを通過する際の有効光路長を示しており、ＫΔｃは、散乱サンプルにおける被測定成分濃度のΔｃの増加による光路長の変化のパーセントを示している。
【００２６】
式（１８）と式（１９）から、下記の式（２０）で示された被測定成分の濃度特異性信号を得ることができる。
【数１３】
2011522264.tif 000015
基準波長が存在するという条件によって、基準波長λrにおいてはＡＧ(λr)＝０である。即ち、
【数１４】
2011522264.tif 000016
である。
【００２７】
式（２１）を式（２０）に代入すると、人体の組織等の複雑な散乱サンプルから基準波長に基づいて、被測定成分の特異性信号を抽出する下記の式（２２）が求められる。
【数１５】
2011522264.tif 000017
【００２８】
本発明に提案された基準波長のデータ処理方法によって、基準波長と測定波長におけるサンプルおよびバックグラウンドのそれぞれのスペクトル信号を測定した上で、式（１７）又は（２２）によって被測定成分の特異性信号を求めることができるので、システムのドリフト及び各種のバックグラウンド変化による干渉を有効的に排除でき、被測定成分の測定精度を向上させることができる。この方法は、血糖濃度などの人の体内にある微弱な成分の検出に適用できる。
【００２９】
（スペクトルの温度補正方法）
近赤外スペクトル法を用いて成分濃度を検出する時に、測定されるスペクトルはサンプル温度の変動によって顕著なずれが生じるので、成分濃度の測定精度への影響が生じてしまう。本発明は、基準波長の特性に基づいて、サンプルのスペクトルに対するリアルタイム温度補正を行う方法を提案した。
【００３０】
ブドウ糖の測定を例として説明する。ここで、ブドウ糖と蒸留水との体積置換係数をｋと定義する。即ち、
【数１６】
2011522264.tif 000018
また、ブドウ糖水溶液の置換吸光係数は、下記の式（２４）のように示されている。
【数１７】
2011522264.tif 000019
式（２）及び式（２４）において、αｗ、αｇ、ｃｗ及びｋのうちのいずれの値も温度に関連している。
【００３１】
異なる温度下で、ブドウ糖水溶液と、相応する蒸留水とのスペクトルを測定して、異なる温度下でのブドウ糖水溶液の吸光係数αTtotal(λ)を算出する。図８に示すように、αTtotal(λ)は温度の変化に従って変更しないことである。
【００３２】
バックグラウンド溶液の純水温度をＴ１とし、ブドウ糖水溶液の温度をＴ２とすると、二者間における吸光度の差は下記の式（２５）のように示される。
【数１８】
2011522264.tif 000020
ここで、
【数１９】
2011522264.tif 000021
と定義する。
【００３３】
異なる温度下でのサンプルの吸光度の差は、温度の変化による吸光度差と、糖濃度の変化による吸光度差との和として示すことができる。即ち、
【数２０】
2011522264.tif 000022
ここで、ΔＡＧ(λ)は、測定中に抽出すべきブドウ糖の特性信号であり、ΔＡＴ(λ)は、温度変化によって引き起こされた干渉信号を示している。
【００３４】
基準波長の定義によれば、基準波長においてΔＡＧ（λｒ）＝０、従って、
ΔＡ(λr)＝ΔＡＴ(λr) ・・・（２７）
となる。
【００３５】
このため、基準波長における吸光度の変化は、完全に蒸留水の吸収係数が温度につれて変化することによるものである。故に、基準波長における吸光度の差の値によれば、ブドウ糖溶液サンプルの測定時の温度と、純水のバックグラウンドの測定時の温度との実際温度の差を算出することができる。
【００３６】
実験を行うことにより、異なる温度の下での蒸留水の吸光係数及び質量濃度を確定することができる。このため、基準波長解析法を用いてブドウ糖溶液サンプルと純水バックグラウンドとの測定温度を取得した場合には、対応するΔＡＴ(λ)の値を算出することができる。従って、式（２６）に基づいて、温度の干渉が排除された後のブドウ糖の特徴信号を以下のように取得できる。
ΔＡＧ(λ)＝ΔＡ(λ)−ΔＡＴ(λ)・・・（２８）
ここで、上述スペクトルに対する温度補正方法について、実験により検証を行った。具体的には、３０℃の温度下で、純水バックグラウンドのスペクトルを測定し、その後、３０℃、３４℃及び３８℃の温度下で、30mg／dLのブドウ糖溶液サンプルのスペクトルを測定した。そして、式（２５）に基づいて、直接的にサンプルとバックグラウンドとの吸光度の差を算出した。
【００３７】
図９に示された通り、サンプルの温度変化の影響により、同一のブドウ糖濃度を有するサンプルが、異なる温度下での吸光度差の値の情報には明らかな違いが見られる。このようなスペクトルを用いるモデルの解析を行うと、必ず温度変化の干渉が存在するため、ブドウ糖濃度の測定精度に影響されることになる。これに対して、基準波長を用いたスペクトル温度補正方法に基づいて、上述したスペクトルを数学的な処理を行うことによって、ブドウ糖の特徴的なスペクトル信号を取得した。図１０に示された通り、異なる温度下でのスペクトル信号は完全に重なり合っている。従って、本発明に提案された、基準波長の技術に基づいてスペクトルの温度補正を行う方法によれば、スペクトルに対する温度による干渉を効率良く排除することができるため、近赤外スペクトル成分の濃度検出技術の測定精度を高めることができる。
【００３８】
以上で述べられた内容は、本発明の実施例に限られるものであり、本発明によって提案されたシステム又は方法に対する限定ではない。本発明の保護される範囲は、添付した請求の範囲及びその等価物により限定されるべきものである。当業者は、本発明の範囲および主旨を逸脱することなしに本発明に種々の補正、切り替え及び変更を行うことが可能であることを理解できる。
【図１】
2011522264.tif 000027
【図２】
2011522264.tif 000028
【図３】
2011522264.tif 000029
【図４】
2011522264.tif 000030
【図５】
2011522264.tif 000031
【図６】
2011522264.tif 000032
【図７】
2011522264.tif 000033
【図８】
2011522264.tif 000034
【図９】
2011522264.tif 000035
【図１０】
2011522264.tif 000036
【国際調査報告】
2011522264.tif 000037
2011522264.tif 000038
2011522264.tif 000039
2011522264.tif 000040