(57)【特許請求の範囲】
【請求項１】
互いに異なる波長の複数の光を同一の周波数で且つ異なる位相の信号によりそれぞれ強度変調して被測定物に照射する光照射手段と、
前記複数の強度変調光のうち少なくとも１つの光の強度を変化させる光強度制御手段と、
前記被測定物の温度を制御する温度制御手段と、
光照射によって前記被測定物から発生する光音響信号を検出して電気信号を出力する光音響信号検出手段と、
前記光強度制御手段が変化させた光の強度を測定する光強度測定手段と、
前記被測定物の温度が異なる２つの時点について、それぞれ前記電気信号の強度が最低となったときの前記光強度測定手段の測定結果から得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記光音響信号の共振モードである音響モードを推定する音響モード推定手段と、
この音響モード推定手段が推定した音響モードと、前記被測定物に含まれる測定対象の成分に対する濃度測定感度と音響モードとの既知の関係に基づいて、前記測定対象の成分に対する濃度測定感度の補正係数を決定する補正係数決定手段と、
前記光強度測定手段の測定結果から得られる光強度の変化量と前記補正係数に基づいて前記測定対象の成分の濃度を導出する濃度導出手段とを備え、
前記光照射手段は、異なる２つの波長の光を前記被測定物に同時に照射し、
前記音響モード推定手段は、第１の波長の光の強度を一定にして第２の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第２の波長の光の強度の測定結果と、前記温度制御手段により前記被測定物の温度を変化させた後に前記第１の波長の光の強度を一定にして前記第２の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第２の波長の光の強度の測定結果とから得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記音響モードを推定することを特徴とする成分濃度測定装置。
【請求項２】
請求項１記載の成分濃度測定装置において、
前記被測定物に照射される複数の光は、これら複数の光の波長に対する、前記被測定物に含まれる測定対象の成分以外の成分の光吸収係数が略等しいことを特徴とする成分濃度測定装置。
【請求項３】
請求項１または２記載の成分濃度測定装置において、
前記光照射手段は、異なる波長の複数の光をそれぞれ平行光にして前記被測定物の略同一領域に同軸で照射し、
異なる波長の複数の光のビーム径は略等しいことを特徴とする成分濃度測定装置。
【請求項４】
請求項１乃至３のいずれか１項に記載の成分濃度測定装置において、
さらに、前記被測定物を収容し、前記光音響信号を増幅するための共振器を備えることを特徴とする成分濃度測定装置。
【請求項５】
請求項４記載の成分濃度測定装置において、
前記共振器は、前記光照射手段から照射される光の光軸に対して略直交する面を２つ持ち、これらの面は互いに平行平面となっていることを特徴とする成分濃度測定装置。
【請求項６】
互いに異なる波長の複数の光を同一の周波数で且つ異なる位相の信号によりそれぞれ強度変調して被測定物に照射する光照射ステップと、
前記複数の強度変調光のうち少なくとも１つの光の強度を変化させる光強度制御ステップと、
光照射によって前記被測定物から発生する光音響信号を検出して電気信号を出力する光音響信号検出ステップと、
前記光強度制御ステップで変化させた光の強度を測定する光強度測定ステップと、
前記被測定物の温度が異なる２つの時点について、それぞれ前記電気信号の強度が最低となったときの前記光強度測定ステップの測定結果から得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記光音響信号の共振モードである音響モードを推定する音響モード推定ステップと、
この音響モード推定ステップで推定した音響モードと、前記被測定物に含まれる測定対象の成分に対する濃度測定感度と音響モードとの既知の関係に基づいて、前記測定対象の成分に対する濃度測定感度の補正係数を決定する補正係数決定ステップと、
前記光強度測定ステップの測定結果から得られる光強度の変化量と前記補正係数に基づいて前記測定対象の成分の濃度を導出する濃度導出ステップとを含み、
前記光照射ステップは、異なる２つの波長の光を前記被測定物に同時に照射するステップを含み、
前記音響モード推定ステップは、第１の波長の光の強度を一定にして第２の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第２の波長の光の強度の測定結果と、温度制御手段により前記被測定物の温度を変化させた後に前記第１の波長の光の強度を一定にして前記第２の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第２の波長の光の強度の測定結果とから得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記音響モードを推定するステップを含むことを特徴とする成分濃度測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、成分濃度測定装置、特に、血液成分、例えばグルコースを測定対象とした測定装置に関するものである。
【背景技術】
【０００２】
高齢化により、生活習慣病に対する対応が大きな課題になりつつあり、血液の成分に対する注目が集まっている。しかし、現在までに開発された技術は、採取したサンプルを試薬と反応させたり、専用のセンサに吸着させたりといった１つのセンサで１つの成分を検出するようになっており、マルチなセンシング技術が注目されている。そこで、電磁波を照射したときに、測定対象とする血液成分（例えば、血糖値の場合はグルコース分子）に吸収されることを利用して、透過あるいは反射した電磁波を測定する方法が行われている。
【０００３】
しかし、グルコースと電磁波との相互作用は小さく、また生体に安全に照射し得る電磁波の強度には制限があり、さらに生体が電磁波に対して散乱体であるために、生体の血糖値測定においては、十分な効果を挙げるに至っていない。
そこで、電磁波を照射し局所的に加熱して熱膨張を起こして発生する音波を観測する、光音響法が注目されている。
【０００４】
光音響法には、パルス（pulse）法と連続波（continuous-wave、以下ＣＷとする）法の二つの方式がある。しかし、従来のパルス法やＣＷ法では、数回にわたる血漿中のグルコース濃度測定中に、グルコース濃度以外の他の血漿中パラメータ（例えば体温や、他の成分の濃度等）も変わる可能性が高いので、グルコース選択性が悪く、正確なグルコース濃度を得ることが難しいという問題点があった。
【０００５】
そこで、ＣＷ法に属する方法として、カムー  セルジュ（S.Camou）らにより光パワーバランスシフト（Optical power balance shift：ＯＰＢＳ）法が開発されている（特許文献１−３参照）。
【０００６】
このＯＰＢＳ法を実施するための従来の成分濃度測定装置の構成例を図１６に示す。成分濃度測定装置は、レーザ光を放射するレーザダイオード１−１，１−２と、レーザダイオード１−１，１−２を駆動するレーザドライバ２と、レーザダイオード１−１，１−２から放射されたレーザ光を導く光ファイバ３−１，３−２と、レーザダイオード１−１，１−２から放射されたレーザ光を合波する光カプラ４と、光カプラ４によって合波されたレーザ光を導く光ファイバ５と、光音響効果によって被測定物１１から発生する光音響信号を検出し、音圧に比例した電気信号に変換する音響センサ６と、音響センサ６から出力された電気信号を増幅する増幅器７と、参照信号を発生するファンクションジェネレータ８と、増幅器７の出力信号とファンクションジェネレータ８から出力された参照信号とを入力として、増幅器７の出力信号から所望の周波数の測定信号を検出するロックインアンプ９と、ファンクションジェネレータ８およびロックインアンプ９を制御すると共に、ロックインアンプ９が検出した測定信号を処理して特定の成分濃度を導出するコンピュータからなる情報処理装置１０とから構成される。
【０００７】
図１６に示した成分濃度測定装置を用いた従来のＯＰＢＳ法では、まずレーザドライバ２は、ファンクションジェネレータ８から出力される参照信号に応じて矩形波の駆動電流をレーザダイオード１−１（第１の光源）に供給することにより、レーザダイオード１−１から放射される光を強度変調する。レーザダイオード１−１から放射された強度変調光は、光ファイバ３−１によって導かれ、光カプラ４および光ファイバ５を介して被測定物１１に照射される（図１７ステップＳ１００）。
【０００８】
音響センサ６は、被測定物１１から発生する光音響信号を検出し、増幅器７は、音響センサ８から出力された電気信号を増幅する。ロックインアンプ９は、増幅器７の出力に含まれる信号のうち、ファンクションジェネレータ８から出力される参照信号によって決まる周波数の測定信号を検出する。
【０００９】
情報処理装置１０は、ファンクションジェネレータ８が発生する参照信号の周波数を変化させることにより、レーザドライバ２からレーザダイオード１−１に供給される駆動電流の周波数を変化させ、光変調周波数を漸次変化させると共に、ロックインアンプ９が検出する測定信号の周波数（光変調周波数と同一の周波数）を漸次変化させる光変調周波数掃引を行う（図１７ステップＳ１０１）。こうして、変調周波数と測定信号の振幅（音圧）との関係を取得し（図１７ステップＳ１０２）、測定信号の振幅が最大の変調周波数であるピーク変調周波数を探索する（図１７ステップＳ１０３）。
【００１０】
次に、レーザダイオード１−１を再び動作させて光を出力させ（図１７ステップＳ１０４）、続いてレーザダイオード１−２（第２の光源）を動作させて光を出力させる（図１７ステップＳ１０５）。２つのレーザダイオード１−１，１−２から放射される光の波長は異なる。
【００１１】
情報処理装置１０は、ファンクションジェネレータ８が発生する２つの参照信号の周波数を変化させることにより、レーザドライバ２からレーザダイオード１−１，１−２に供給される駆動電流の周波数を変化させ、光変調周波数を上記のピーク変調周波数に設定すると共に、ロックインアンプ９が検出する測定信号の周波数をピーク変調周波数に設定する。このとき、レーザドライバ２から、同一周波数で逆位相の矩形波の駆動電流をレーザダイオード１−１，１−２に供給させることにより、レーザダイオード１−１，１−２から放射される光を同一周波数（ピーク変調周波数）で逆位相の信号によりそれぞれ強度変調する（図１７ステップＳ１０６）。
【００１２】
レーザダイオード１−１，１−２から放射された強度変調光は、それぞれ光ファイバ３−１，３−２によって導かれ、光カプラ４によって合波され、さらに光ファイバ５によって導かれ、被測定物１１に照射される。情報処理装置１０は、レーザドライバ２からレーザダイオード１−１に供給される駆動電流の大きさを所定値にすることにより、レーザダイオード１−１から放射される光の強度を固定し（図１７ステップＳ１０７）、レーザドライバ２からレーザダイオード１−２に供給される駆動電流の大きさを変化させることにより、レーザダイオード１−２から放射される光の強度を漸次変化させる光強度掃引を行い、測定信号の振幅が最小になるようにレーザダイオード１−２から放射される光の強度を調節する（図１７ステップＳ１０８）。
【００１３】
そして、一定時間経過後（図１７ステップＳ１０９）、再びレーザダイオード１−１，１−２を動作させ、レーザダイオード１−１から放射される光の強度を固定し（図１７ステップＳ１１０）、測定信号の振幅が最小になるようにレーザダイオード１−２から放射される光の強度を調節する（図１７ステップＳ１１１）。そして、情報処理装置１０は、ステップＳ１０８で得られた光強度とステップＳ１１１で得られた光強度から、光強度の変化量を計算する（図１７ステップＳ１１２）。
【００１４】
このように、従来のＯＰＢＳ法では、光波長が異なり位相差が１８０°の２つの光ビームを被測定物に対して、同一光出力口から照射し、２つの光ビームの強度を増減させながら、光音響信号の振幅が最小の箇所の位相の変曲点を探し、この結果から血液中に溶解している分子濃度を求める。具体的には、２つの光ビームのうち一方の光ビームの強度を変えながら、光音響信号の振幅が最小となる光強度を探すことで、光強度の変化量より被測定物内の特定成分（例えばグルコース）の濃度の正確な測定を行う。
【先行技術文献】
【特許文献】
【００１５】
【特許文献１】特開２０１４−５０５６３号公報
【特許文献２】特開２０１３−１０６８７４号公報
【特許文献３】特開２０１２−１７９２１２号公報
【発明の概要】
【発明が解決しようとする課題】
【００１６】
被測定物を収容するケースが円筒形、すなわち光音響エネルギーが閉じ込められる空間が円筒空洞である場合の、この円筒空洞内での共振を例にとった音波の共振モードを図１８に示す。ここで、Ｒは円筒の半径、Ｌは円筒の長さ、ｊは半径方向（Radial）のモードの番号、ｍは円周方向（Azimuthal）のモードの番号、ｑは長手方向（Longitudinal）のモードの番号である。
【００１７】
変調周波数と音圧との関係を図１９に示す。図１９中の（１０１）、（２００）、（１０２）、（３００）は共振モード（ｊｍｑ）を示しており、例えば（１０１）は半径方向のモードの番号ｊが１、円周方向のモードの番号ｍが０、長手方向のモードの番号ｑが１であることを示している。このように、円筒形の共振器内では、音波は様々なモードで共振し、それらが複雑に干渉し合った音波の分布が生じている。
【００１８】
ＣＷ光を用いた従来の光音響法においては、光の照射により生じる音波の分布と定在波の音波の干渉によって測定値が大きく変わってしまうという問題があった。すなわち、被測定物に含まれる対象の成分の濃度の検量線の傾きが数倍のオーダーで変わってしまうという問題があった。
【００１９】
本発明は、上記課題を解決するためになされたもので、被測定物に含まれる測定対象の成分に対する濃度測定感度の音響モード依存性を補償することができる成分濃度測定装置および方法を提供することを目的とする。
【課題を解決するための手段】
【００２０】
本発明の成分濃度測定装置は、互いに異なる波長の複数の光を同一の周波数で且つ異なる位相の信号によりそれぞれ強度変調して被測定物に照射する光照射手段と、前記複数の強度変調光のうち少なくとも１つの光の強度を変化させる光強度制御手段と、前記被測定物の温度を制御する温度制御手段と、光照射によって前記被測定物から発生する光音響信号を検出して電気信号を出力する光音響信号検出手段と、前記光強度制御手段が変化させた光の強度を測定する光強度測定手段と、前記被測定物の温度が異なる２つの時点について、それぞれ前記電気信号の強度が最低となったときの前記光強度測定手段の測定結果から得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記光音響信号の共振モードである音響モードを推定する音響モード推定手段と、この音響モード推定手段が推定した音響モードと、前記被測定物に含まれる測定対象の成分に対する濃度測定感度と音響モードとの既知の関係に基づいて、前記測定対象の成分に対する濃度測定感度の補正係数を決定する補正係数決定手段と、前記光強度測定手段の測定結果から得られる光強度の変化量と前記補正係数に基づいて前記測定対象の成分の濃度を導出する濃度導出手段とを備え、前記光照射手段は、異なる２つの波長の光を前記被測定物に同時に照射し、前記音響モード推定手段は、第１の波長の光の強度を一定にして第２の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第２の波長の光の強度の測定結果と、前記温度制御手段により前記被測定物の温度を変化させた後に前記第１の波長の光の強度を一定にして前記第２の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第２の波長の光の強度の測定結果とから得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記音響モードを推定することを特徴とするものである。
【００２３】
また、本発明の成分濃度測定装置の１構成例において、前記被測定物に照射される複数の光は、これら複数の光の波長に対する、前記被測定物に含まれる測定対象の成分以外の成分の光吸収係数が略等しいことを特徴とするものである。
また、本発明の成分濃度測定装置の１構成例において、前記光照射手段は、異なる波長の複数の光をそれぞれ平行光にして前記被測定物の略同一領域に同軸で照射し、異なる波長の複数の光のビーム径は略等しいことを特徴とするものである。
また、本発明の成分濃度測定装置の１構成例は、さらに、前記被測定物を収容し、前記光音響信号を増幅するための共振器を備えることを特徴とするものである。
また、本発明の成分濃度測定装置の１構成例において、前記共振器は、前記光照射手段から照射される光の光軸に対して略直交する面を２つ持ち、これらの面は互いに平行平面となっていることを特徴とするものである。
【００２４】
また、本発明の成分濃度測定方法は、互いに異なる波長の複数の光を同一の周波数で且つ異なる位相の信号によりそれぞれ強度変調して被測定物に照射する光照射ステップと、前記複数の強度変調光のうち少なくとも１つの光の強度を変化させる光強度制御ステップと、光照射によって前記被測定物から発生する光音響信号を検出して電気信号を出力する光音響信号検出ステップと、前記光強度制御ステップで変化させた光の強度を測定する光強度測定ステップと、前記被測定物の温度が異なる２つの時点について、それぞれ前記電気信号の強度が最低となったときの前記光強度測定ステップの測定結果から得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記光音響信号の共振モードである音響モードを推定する音響モード推定ステップと、この音響モード推定ステップで推定した音響モードと、前記被測定物に含まれる測定対象の成分に対する濃度測定感度と音響モードとの既知の関係に基づいて、前記測定対象の成分に対する濃度測定感度の補正係数を決定する補正係数決定ステップと、前記光強度測定ステップの測定結果から得られる光強度の変化量と前記補正係数に基づいて前記測定対象の成分の濃度を導出する濃度導出ステップとを含み、前記光照射ステップは、異なる２つの波長の光を前記被測定物に同時に照射するステップを含み、前記音響モード推定ステップは、第１の波長の光の強度を一定にして第２の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第２の波長の光の強度の測定結果と、温度制御手段により前記被測定物の温度を変化させた後に前記第１の波長の光の強度を一定にして前記第２の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第２の波長の光の強度の測定結果とから得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記音響モードを推定するステップを含むことを特徴とするものである。
【発明の効果】
【００２５】
本発明によれば、互いに異なる波長の複数の光を同一の周波数で且つ異なる位相の信号によりそれぞれ強度変調して被測定物に照射する光照射手段と、複数の強度変調光のうち少なくとも１つの光の強度を変化させる光強度制御手段と、被測定物の温度を制御する温度制御手段と、光照射によって被測定物から発生する光音響信号を検出して電気信号を出力する光音響信号検出手段と、光強度制御手段が変化させた光の強度を測定する光強度測定手段と、被測定物の温度が異なる２つの時点について、それぞれ電気信号の強度が最低となったときの光強度測定手段の測定結果から得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、光音響信号の共振モードである音響モードを推定する音響モード推定手段と、音響モード推定手段が推定した音響モードと、被測定物に含まれる測定対象の成分に対する濃度測定感度と音響モードとの既知の関係に基づいて、測定対象の成分に対する濃度測定感度の補正係数を決定する補正係数決定手段とを設けることにより、被測定物内の測定対象の成分に対する濃度測定感度の、音響モードによる変化を補償することができ、測定感度を向上させることができる。
【００２６】
また、本発明では、被測定物を収容する共振器を設けることにより、光照射によって被測定物から発生する光音響信号を増幅することができる。
【００２８】
また、本発明では、第１の波長の光の強度を一定にして第２の波長の光の強度を変化させて電気信号の強度が最低となったときの第２の波長の光の強度の測定結果と、温度制御手段により被測定物の温度を変化させた後に第１の波長の光の強度を一定にして第２の波長の光の強度を変化させて電気信号の強度が最低となったときの第２の波長の光の強度の測定結果とから得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、音響モードを推定することができる。
【００２９】
また、本発明では、光照射手段が、異なる波長の複数の光をそれぞれ平行光にして被測定物の略同一領域に同軸で照射することにより、音波の干渉の影響を低減することができる。
【図面の簡単な説明】
【００３０】
【図１】グルコース濃度の応答特性を示す図である。
【図２】変調周波数と音波強度の関係を示す図である。
【図３】本発明の参考例に係る成分濃度測定装置の構成を示すブロック図である。
【図４】本発明の参考例に係る成分濃度測定装置の情報処理装置の構成を示すブロック図である。
【図５】本発明の参考例に係る成分濃度測定装置の動作を説明するフローチャートである。
【図６】本発明の参考例における変調周波数と音波強度との関係の１例を示す図である。
【図７】本発明の参考例における光強度変化の例を示す図である。
【図８】本発明の参考例における音響モードと光強度の変化量の比との関係の１例を示す図である。
【図９】本発明の参考例におけるグルコース濃度感度と音響モードとの関係の１例を示す図である。
【図１０】本発明の実施の形態に係る成分濃度測定装置の構成を示すブロック図である。
【図１１】本発明の実施の形態に係る成分濃度測定装置の情報処理装置の構成を示すブロック図である。
【図１２】本発明の実施の形態に係る成分濃度測定装置の動作を説明するフローチャートである。
【図１３】本発明の実施の形態における光強度変化の例を示す図である。
【図１４】本発明の実施の形態における音響モードと光強度の変化量の比との関係の１例を示す図である。
【図１５】音波の反射率と音響インピーダンスとの関係を示す図である。
【図１６】従来の成分濃度測定装置の構成を示すブロック図である。
【図１７】従来の光パワーバランスシフト法を説明するフローチャートである。
【図１８】円筒空洞内での音波の共振モードを説明する図である。
【図１９】変調周波数と音圧との関係の１例を示す図である。
【発明を実施するための形態】
【００３１】
［発明の原理］
本発明では、コリメートされた光を用いることで音波の干渉の影響を低減する。被測定物の形状が自由に変更できる場合は、音波の波長の半分に対して径の方が細い円筒形状の被測定物を用いることで、コリメート光と同等の効果を得ることができる。コリメートされた光は、強度分布として以下の式（１）のように表すことができる。
【００３２】
【数１】
0006606018.tif 000002
【００３３】
また、光音響効果で生じる音波の大きさＳは、以下の式（２）のように表すことができる。
【００３４】
【数２】
0006606018.tif 000003
【００３５】
式（１）、式（２）において、Ｃは定数、αは特定の波長の光に対する被測定物の光吸収係数、ｑは音響モード、Ｆ（α，ｑ）は光吸収係数αと音響モードｑの関数、ωは変調周波数、ｗは光ビーム径、ｚは被測定物表面からの深さ、ｒは光軸から放射方向の距離、Ｉは光強度である。
光吸収係数αは式（３）のように記述することができる。
【００３６】
【数３】
0006606018.tif 000004
【００３７】
式（３）において、ΔＣは被測定物内の特定成分の濃度、ΔＴは被測定物の温度、βは特定成分のモル吸光度、γは被測定物の温度吸光度である。なお、α0は初期状態の被測定物の光吸収係数を表している。
【００３８】
たとえば、円筒形の共振器内においては、音響波のモードは以下のような式（４）、式（５）で与えられる。音響波のモードには、長手方向、円周方向、半径方向のモードがある。
【００３９】
【数４】
0006606018.tif 000005
【００４０】
【数５】
0006606018.tif 000006
【００４１】
式（４）、式（５）において、ｃは被測定物内の音速、ｑは光軸方向の音響モード、Ｒは円筒の半径、Ｌは円筒の長さ、αm,nは被測定物壁面で径方向の音圧変化が０という境界条件を満たすｎ次の根である。
【００４２】
例えば、コリメートされた光を被測定物に照射する場合、ｑは円筒軸方向の定在波モードに相当する（図１８）。また、関数Ｆは、この場合、近似的に式（６）のように記述することができるため、モードごとに離散的な値をとる。
【００４３】
【数６】
0006606018.tif 000007
【００４４】
ここで、光吸収係数αi、関数Ｆiのｉは光の波長に対応する添え字である。平行平板間の音波の共振においても、円筒内の音波の共振と同様に扱うことができる。ここで、波長が異なり、１８０°位相がずれた２つの光を被測定物に照射する場合に発生する音波の大きさＳは、それぞれの光により生ずる音波の重ね合わせとして次のように表すことができる。
【００４５】
【数７】
0006606018.tif 000008
【００４６】
Ｓ1，Ｓ2は２つの音波の大きさである。ＯＰＢＳ法では、音波の信号強度が最少、すなわち理想的にゼロになる状態を考えることと同等であることから、Ｓ1＝Ｓ2となる。
被測定物内の特定の成分がΔＣだけ濃度変化したとすると、以下の式（８）が成立する。
【００４７】
【数８】
0006606018.tif 000009
【００４８】
Ｉ1は第１の光の強度、Ｉ2は第２の光の強度、ΔＩ2は被測定物内の特定成分がΔＣだけ濃度変化したときの第２の光の強度の変化量である。式（８）から、当該成分濃度変化と光強度変化の関係は式（９）、式（１０）のように近似できる。
【００４９】
【数９】
0006606018.tif 000010
【００５０】
【数１０】
0006606018.tif 000011
【００５１】
式（９）、式（１０）より音響モードｑおよび背景光吸収が分かれば、濃度変化の感度が求められる。被測定物をほぼ水と仮定し、グルコース濃度が変化したときの応答特性の例を図１に示す。図１によれば、グルコース濃度に比例して光強度が変化していくこと、音響モードｑにより感度が変化することが示される。また、音響モードｑが変化すると、感度が大きく変化してしまうことが分かる。
【００５２】
そこで、音響モードの補正を行う必要がある。この補正は、測定に用いる波長の光とは別の第３の波長の光あるいは温度に対する吸収スペクトルのデータセットを用意し、実測データに対する数値的補正を行うことで実現することができる。式（９）、式（１０）および既知の吸収係数差を用いて音波の共振モードの干渉誤差を補正する後述の参考例の方法と本発明の方法の２つについて説明する。
【００５３】
［参考例の方法］
第１の波長の光を被測定物に照射し、光の変調周波数を任意の範囲で掃引することで図２に示すような変調周波数と音圧との関係を取得する。この図２のピークから被測定物内での音波の共振周波数を選択することができる。
【００５４】
次に前述のとおり、光波長が異なり位相差が１８０°の２つの光ビームを被測定物に対して、同一の光出力口から照射する。このときの光変調周波数は上記で選択した変調周波数である。第１の波長の光の強度を一定に固定し、第２の波長の光の強度を増減させながら、生じる音波強度が極小となる、第２の波長の光の強度Ｉ0を測定する。
【００５５】
次に、第２の波長の光を被測定物に対する光吸収係数が既知の第３の波長の光に切り替えて、上記と同様に第３の波長の光の強度を増減させながら、生じる音波強度が極小となる、第３の波長の光の強度Ｉ1を測定する。求めた光強度Ｉ0，Ｉ1より光強度Ｉ0に対する光強度の変化量（Ｉ1−Ｉ0）の比ΔＩ／Ｉは式（１１）のようになる。
【００５６】
【数１１】
0006606018.tif 000012
【００５７】
ここで、十分コリメートされた光を用いた場合、上式の関係は各音響モードｑに対して次式で与えられる。
【００５８】
【数１２】
0006606018.tif 000013
【００５９】
【数１３】
0006606018.tif 000014
【００６０】
関数Ｆi，光吸収係数αiのｉは光の波長に対応する添え字である。Ｆ1は第１の波長の光に関する光吸収係数αと音響モードｑの関数、Ｆ2は第２の波長の光に関する光吸収係数αと音響モードｑの関数、Ｆ1’は第３の波長を用いるときの第１の波長の光に関する光吸収係数αと音響モードｑの関数、Ｆ2’は第３の波長の光に関する光吸収係数αと音響モードｑの関数、α2は第２の波長の光に対する被測定物の光吸収係数である。
【００６１】
よって、式（１１）と式（１２）の応答を比較し、最も濃度測定感度が近しい応答を探し音響モードｑを求めることができる。音響モードｑを用いて式（９）より現在の変調周波数における濃度測定感度Ｇｑが推測できる。この濃度測定感度Ｇｑを音響モードｑ＝０のときの濃度測定感度Ｇ０と同じになるように定数係数Ａを決定すると、補償された濃度との関係は次式で表される。
【００６２】
【数１４】
0006606018.tif 000015
【００６３】
［本発明の方法］
参考例の方法と同様に、第１の波長の光を被測定物に照射し、光の変調周波数を任意の範囲で掃引することで図２に示すような変調周波数と音圧との関係を取得する。この図２のピークから被測定物内での音波の共振周波数を選択することができる。
【００６４】
次に前述のとおり、光波長が異なり位相差が１８０°の２つの光ビームを被測定物に対して、同一光出力口から照射する。このときの変調周波数は上記で選択した変調周波数である。第１の波長の光の強度を一定に固定し、第２の波長の光の強度を増減させながら、生じる音波強度が極小となる、第２の波長の光の強度Ｉ0を測定する。
【００６５】
次に、温度制御器を用いて被測定物の温度をΔＴだけ変化させ、上記と同様に第２の波長の光の強度を増減させながら、生じる音波強度が極小となる、第２の波長の光の強度Ｉ1を測定する。求めた光強度Ｉ0，Ｉ1より光強度Ｉ0に対する光強度の変化量（Ｉ1−Ｉ0）の比ΔＩ／Ｉは式（１５）のようになる。
【００６６】
【数１５】
0006606018.tif 000016
【００６７】
被測定物の温度変化に対する光吸収係数変化は、スペクトルデータセットより光吸収係数の温度依存性γで与えた場合、次式で温度に対する応答特性が与えられる。
【００６８】
【数１６】
0006606018.tif 000017
【００６９】
【数１７】
0006606018.tif 000018
【００７０】
γ1は第１の波長の光に対する被測定物の光吸収係数α1の温度依存性、γ2は第２の波長の光に対する被測定物の光吸収係数α2の温度依存性である。
【００７１】
よって、式（１５）と式（１６）の応答を比較し、最も感度が近しい応答を探し音響モードｑを求めることができる。音響モードｑを用いて式（９）より現在の変調周波数における濃度測定感度Ｇｑが推測できる。この濃度測定感度Ｇｑを音響モードｑ＝０のときの濃度測定感度Ｇ０と同じになるように定数係数Ａを決定すると、補償された濃度との関係は次式で表される。
【００７２】
【数１８】
0006606018.tif 000019
【００７３】
光音響法の音波共振モードの補償方法についてグルコース水溶液の場合を例にとって下記に示す。
【００７４】
［参考例］
以下、本発明の参考例について図面を参照して説明する。図３は本発明の参考例に係る成分濃度測定装置の構成を示すブロック図である。本参考例の成分濃度測定装置は、レーザ光を放射するレーザダイオード１−１，１−２，１−３と、レーザダイオード１−１，１−２，１−３を駆動するレーザドライバ２と、レーザダイオード１−１，１−２，１−３から放射されたレーザ光を導く光ファイバ３−１，３−２，３−３と、レーザダイオード１−１，１−２，１−３から放射されたレーザ光を合波する光カプラ４と、光カプラ４によって合波されたレーザ光を導く光ファイバ５と、光音響効果によって被測定物１１から発生する光音響信号を検出し、音圧に比例した電気信号に変換する光音響信号検出手段となる音響センサ６と、音響センサ６から出力された電気信号を増幅する増幅器７と、参照信号を発生するファンクションジェネレータ８と、増幅器７の出力信号とファンクションジェネレータ８から出力された参照信号とを入力として、増幅器７の出力信号から所望の周波数の測定信号を検出するロックインアンプ９と、レーザドライバ２とファンクションジェネレータ８とロックインアンプ９とを制御すると共に、ロックインアンプ９が検出した測定信号を処理して特定の成分濃度を導出するコンピュータからなる情報処理装置１０とから構成される。
【００７５】
レーザダイオード１−１，１−２，１−３とレーザドライバ２とファンクションジェネレータ８とは、光照射手段を構成している。情報処理装置１０とファンクションジェネレータ８とは、光強度制御手段を構成している。
【００７６】
レーザダイオード１−１，１−２，１−３の例としては、例えば分布帰還型半導体レーザ（ＤＦＢ−ＬＤ）等がある。音響センサ６の例としては、圧電センサを用いるマイクロホンがある。
なお、被測定物１１は例えば円筒形の光音響セル（不図示）内に導入され、光音響セルに設けられたガラス製の光学窓を通してレーザ光が照射されるようになっている。
【００７７】
図４は情報処理装置１０の構成を示すブロック図である。情報処理装置１０は、ファンクションジェネレータ８を制御するジェネレータ制御部１００と、ファンクションジェネレータ８を介して光の強度を制御する光強度制御部１０１と、レーザダイオード１−１〜１−３から放射された光の強度を測定する光強度測定手段となる光強度測定部１０２と、強度を変化させる光の波長が異なる２つの時点または被測定物１１の温度が異なる２つの時点について、それぞれ測定信号の強度が最低となったときの光強度の測定結果から得られる光強度の変化量を算出する光強度変化量算出部１０３と、強度を変化させる光の波長が異なる２つの時点または被測定物１１の温度が異なる２つの時点について、それぞれ測定信号の強度が最低となったときの光強度の測定結果から得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、光音響信号の共振モードである音響モードを推定する音響モード推定部１０４と、音響モード推定部１０４が推定した音響モードと、被測定物１１に含まれる特定成分（例えばグルコース）に対する濃度測定感度と音響モードとの既知の関係に基づいて、特定成分に対する濃度測定感度の補正係数を決定する補正係数決定部１０５と、被測定物１１内の特定成分の濃度を導出する濃度導出部１０６と、情報記憶のための記憶部１０７とを有する。
【００７８】
次に、本参考例の成分濃度測定装置の動作を図５のフローチャートを参照して説明する。初めに、レーザダイオード１−１（第１の光源）のみを動作させる。レーザドライバ２から駆動電流が供給されると、レーザダイオード１−１はレーザ光（連続光）を放射する。このとき、レーザドライバ２から矩形波の駆動電流が供給されることにより、レーザダイオード１−１は強度変調光を放射する。レーザダイオード１−１から放射される光の波長（第１の波長）は例えば１３８２ｎｍである。レーザダイオード１−１から放射された強度変調光は、光ファイバ３−１によって導かれ、光カプラ４および光ファイバ５を介して被測定物１１に照射される（図５ステップＳ２００）。
【００７９】
音響センサ６は、被測定物１１から発生する光音響信号を検出し、増幅器７は、音響センサ８から出力された電気信号を増幅する。ロックインアンプ９は、増幅器７の出力に含まれる信号のうち、ファンクションジェネレータ８から出力される参照信号によって決まる周波数の測定信号を検出する。
【００８０】
情報処理装置１０のジェネレータ制御部１００は、ファンクションジェネレータ８が発生する参照信号の周波数を変化させることにより、レーザドライバ２からレーザダイオード１−１に供給される駆動電流の周波数を変化させ、光変調周波数を漸次変化させると共に、ロックインアンプ１１が検出する測定信号の周波数（光変調周波数と同一の周波数）を漸次変化させる光変調周波数掃引を行う（図５ステップＳ２０１）。
【００８１】
こうして、変調周波数と測定信号の振幅（音圧）との関係を取得し（図５ステップＳ２０２）、測定信号の振幅が最大の変調周波数であるピーク変調周波数を探索する（図５ステップＳ２０３）。本参考例において得られた変調周波数と音波強度との関係を図６に示す。
【００８２】
次に、レーザダイオード１−１を再び動作させて連続光を出力させ（図５ステップＳ２０４）、続いてレーザダイオード１−２（第２の光源）を動作させて連続光を出力させる（図５ステップＳ２０５）。レーザダイオード１−１から放射される光の波長は例えば１３８２ｎｍ、レーザダイオード１−２から放射される光の波長（第２の波長）は例えば１６１０ｎｍであり、２つのレーザダイオード１−１，１−２から放射される光の波長は異なる。
【００８３】
情報処理装置１０のジェネレータ制御部１００は、ファンクションジェネレータ８が発生する２つの参照信号の周波数を変化させることにより、レーザドライバ２からレーザダイオード１−１，１−２に供給される駆動電流の周波数を変化させ、光変調周波数を上記のピーク変調周波数に設定すると共に、ロックインアンプ９が検出する測定信号の周波数をピーク変調周波数に設定する。このとき、レーザドライバ２から、同一周波数で逆位相の矩形波の駆動電流をレーザダイオード１−１，１−２に供給させることにより、レーザダイオード１−１，１−２から放射される光を同一周波数（ピーク変調周波数）で逆位相の信号によりそれぞれ強度変調する（図５ステップＳ２０６）。したがって、本参考例では、互いに異なる波長の２波のレーザ光を同一周波数で逆位相の信号によりそれぞれ強度変調して２つの強度変調光を生成することになる。
【００８４】
レーザダイオード１−１，１−２から放射された強度変調光は、それぞれ光ファイバ３−１，３−２によって導かれ、光カプラ４によって合波され、さらに光ファイバ５によって導かれ、被測定物１１に照射される。
【００８５】
情報処理装置１０の光強度制御部１０１は、ファンクションジェネレータ８を制御して、レーザダイオード１−１を駆動するための参照信号の電圧を所定値に設定することにより、レーザドライバ２からレーザダイオード１−１に供給される駆動電流の大きさを一定にし、レーザダイオード１−１から放射される光の強度を一定にする（図５ステップＳ２０７）。
【００８６】
また、情報処理装置１０の光強度制御部１０１は、ファンクションジェネレータ８を制御して、レーザダイオード１−２を駆動するための参照信号の電圧を変化させることにより、レーザドライバ２からレーザダイオード１−２に供給される駆動電流の大きさを変化させ、レーザダイオード１−２から放射される光の強度を漸次変化させる光強度掃引を行い、ロックインアンプ９から出力される測定信号の振幅が最小になるようにレーザダイオード１−２から放射される光の強度を調節する（図５ステップＳ２０８）。
【００８７】
情報処理装置１０の光強度測定部１０２は、レーザダイオード１−２から放射される光の強度を測定する。光強度測定部１０２は、ファンクションジェネレータ８から出力されている、レーザダイオード１−２の駆動用の参照信号の電圧を取得する。情報処理装置１０の記憶部１０７には、参照信号の電圧とレーザダイオード１−２から放射される光の強度との関係を示すキャリブレーションデータが予め記憶されている。キャリブレーションデータは、電圧と強度の実測により予め求めておくことができる。光強度測定部１０２は、このようなキャリブレーションデータを参照し、取得した参照信号の電圧を、レーザダイオード１−２から放射された光の強度に換算する。光強度測定部１０２は、測定信号の振幅が最小となるときの、レーザダイオード１−２から放射される光の強度Ｉ0を記憶部１０７に格納する（図５ステップＳ２０９）。
【００８８】
次に、レーザダイオード１−２の動作を停止させて、レーザダイオード１−３（第３の光源）を動作させて連続光を出力させる（図５ステップＳ２１０）。レーザダイオード１−３から放射される光の波長（第３の波長）は例えば１６５０ｎｍであり、２つのレーザダイオード１−１，１−３から放射される光の波長は異なる。
【００８９】
情報処理装置１０のジェネレータ制御部１００は、ファンクションジェネレータ８が発生する２つの参照信号の周波数を変化させることにより、レーザドライバ２からレーザダイオード１−１，１−３に供給される駆動電流の周波数を変化させ、光変調周波数を上記のピーク変調周波数に設定すると共に、ロックインアンプ９が検出する測定信号の周波数をピーク変調周波数に設定する。このとき、レーザドライバ２から、同一周波数で逆位相の矩形波の駆動電流をレーザダイオード１−１，１−３に供給させることにより、レーザダイオード１−１，１−３から放射される光を同一周波数（ピーク変調周波数）で逆位相の信号によりそれぞれ強度変調する（図５ステップＳ２１１）。
【００９０】
情報処理装置１０の光強度制御部１０１は、ファンクションジェネレータ８を制御して、レーザダイオード１−３を駆動するための参照信号の電圧を変化させることにより、レーザドライバ２からレーザダイオード１−３に供給される駆動電流の大きさを変化させ、レーザダイオード１−３から放射される光の強度を漸次変化させる光強度掃引を行い、ロックインアンプ９から出力される測定信号の振幅が最小になるようにレーザダイオード１−３から放射される光の強度を調節する（図５ステップＳ２１２）。上記のとおり、レーザダイオード１−１から放射される光の強度は一定である。
【００９１】
情報処理装置１０の光強度測定部１０２は、レーザダイオード１−３から放射される光の強度を測定する。光強度測定部１０２は、ファンクションジェネレータ８から出力されている、レーザダイオード１−３の駆動用の参照信号の電圧を取得する。上記と同様に、情報処理装置１０の記憶部１０７には、参照信号の電圧とレーザダイオード１−３から放射される光の強度との関係を示すキャリブレーションデータが予め記憶されている。光強度測定部１０２は、このようなキャリブレーションデータを参照し、取得した参照信号の電圧を、レーザダイオード１−３から放射された光の強度に換算する。光強度測定部１０２は、測定信号の振幅が最小となるときの、レーザダイオード１−３から放射される光の強度Ｉ1を記憶部１０７に格納する（図５ステップＳ２１３）。
【００９２】
情報処理装置１０の光強度変化量算出部１０３は、ステップＳ２０９，Ｓ２１３で求めた光強度Ｉ0，Ｉ1より光強度Ｉ0に対する光強度の変化量ΔＩ＝（Ｉ1−Ｉ0）の比ΔＩ／Ｉを、式（１１）により計算する（図５ステップＳ２１４）。
【００９３】
レーザダイオード１−３から放射される第３の波長の光に対する被測定物１１の光吸収係数がレーザダイオード１−２から放射される第２の波長の光に対する被測定物１１の光吸収係数よりも０．０２／ｍｍ大きい場合の光強度の変化量の比ΔＩ／Ｉを図７に示す。図７によると、ΔＩ／Ｉ＝４．７％という結果が得られる。
【００９４】
次に、情報処理装置１０の音響モード推定部１０４は、表１のような所定のデータセットを用いて、光強度の変化量の比ΔＩ2／Ｉ2を式（１２）、式（１３）により音響モードｑ毎に計算する（図５ステップＳ２１５）。データセットは、第１、第２、第３の波長と、これらの波長の光に対する特定成分（例えばグルコース）の光吸収係数との関係を表すものである。
【００９５】
【表１】
0006606018.tif 000020
【００９６】
関数Ｆ1は第１の波長の光に関する光吸収係数α1＝０．３０／ｍｍと音響モードｑから計算することができ、関数Ｆ2は第２の波長の光に関する光吸収係数α2＝０．２８／ｍｍと音響モードｑから計算することができる。Δα2は光吸収係数α1に対する光吸収係数α2の変化量である。被測定物１１の長さＬは既知の値である。
【００９７】
ステップＳ２１５の計算の結果から、図８に示すような音響モードｑと光強度の変化量の比ΔＩ2／Ｉ2との関係が得られる。情報処理装置１０の音響モード推定部１０４は、ステップＳ２１４で計算された光強度の変化量の比ΔＩ／Ｉと、ステップＳ２１５で計算した光強度の変化量の比ΔＩ2／Ｉ2とが等しくなる音響モードｑの値を、今回のピーク変調周波数での音響モードの値として決定する（図５ステップＳ２１６）。本参考例の例では、ｑ＝２であることが推測できる。
【００９８】
情報処理装置１０の補正係数決定部１０５は、音響モード推定部１０４が推定した音響モードｑの値に対応する、被測定物１１内の特定成分（本参考例の例ではグルコース）の濃度の測定感度Ｇｑを、濃度測定感度と音響モードｑとの既知の関係（図９）から求め、この濃度測定感度Ｇｑに対する音響モードｑ＝０のときの濃度測定感度Ｇ０の比Ｇ０／Ｇｑを補正係数Ａとして決定する（図５ステップＳ２１７）。
【００９９】
図９によると、音響モードｑ＝２のときの濃度測定感度Ｇｑは音響モードｑ＝０のときの濃度測定感度Ｇ０より約５倍大きいので、補正係数決定部１０５は、補正係数Ａ＝１／５と求めることができる。
【０１００】
情報処理装置１０の濃度導出部１０６は、ステップＳ２１４で計算された光強度の変化量の比ΔＩ／Ｉ（音響モードｑ＝２のときの光強度の変化量の比ΔＩ2／Ｉ2）と、補正係数決定部１０５が決定した補正係数Ａとから、測定開始時から任意の時間経過後の時点における被測定物１１内の特定成分の濃度の変化量ΔＣ［％］を式（１４）により計算する（図５ステップＳ２１８）。
【０１０１】
最後に、濃度導出部１０６は、ステップＳ２１８で計算した特定成分の濃度の変化量ΔＣ［％］と特定成分の既知の参照濃度とから、任意の時間経過後の時点における特定成分の濃度を計算する（図５ステップＳ２１９）。なお、参照濃度は、被測定物１１に対して標準的な血糖測定法を実施することにより得ることができる。標準的な血糖測定法を実施するには、血糖測定器の本体に、グルコースセンサーを差し込み、針を専用の機械（または本体）にセットして、指などから採血し、グルコースセンサーに血を吸収させる。標準的な血糖測定法は、既知濃度のグルコース液を標準校正液として機械動作確認用に用いる。初期動作時に機械が正常に動いているかを確認したり、血糖値が異常値にあるか（正常に機械が動作しているか）を確認したりするときに用いる。このような標準的な血糖測定法を、測定開始時の例えばピーク変調周波数の選択時に実施しておけばよい。
【０１０２】
以上のようにして、本参考例では、被測定物１１内の特定成分に対する濃度測定感度の、音響モードによる変化を補償することができ、測定感度を向上させることができる。
【０１０３】
［実施の形態］
次に、本発明の実施の形態について説明する。図１０は本発明の実施の形態に係る成分濃度測定装置の構成を示すブロック図であり、図３と同一の構成には同一の符号を付してある。本実施の形態の成分濃度測定装置は、レーザダイオード１−１，１−２と、レーザドライバ２と、光ファイバ３−１，３−２と、光カプラ４と、光ファイバ５と、音響センサ６と、増幅器７と、ファンクションジェネレータ８と、ロックインアンプ９と、情報処理装置１０ａと、被測定物１１の温度を変化させる温度制御器１２とから構成される。
【０１０４】
図１１は本実施の形態の情報処理装置１０ａの構成を示すブロック図である。情報処理装置１０ａは、ジェネレータ制御部１００と、光強度制御部１０１と、光強度測定部１０２と、光強度変化量算出部１０３と、音響モード推定部１０４と、補正係数決定部１０５と、濃度導出部１０６と、記憶部１０７と、被測定物１１の温度を制御する温度制御部１０８とを有する。
【０１０５】
次に、本実施の形態の成分濃度測定装置の動作を図１２のフローチャートを参照して説明する。初めに、レーザダイオード１−１（第１の光源）のみを動作させる。レーザドライバ２から駆動電流が供給されると、レーザダイオード１−１はレーザ光（連続光）を放射する。このとき、レーザドライバ２から矩形波の駆動電流が供給されることにより、レーザダイオード１−１は強度変調光を放射する。参考例と同様に、レーザダイオード１−１から放射される光の波長（第１の波長）は例えば１３８２ｎｍである。レーザダイオード１−１から放射された強度変調光は、光ファイバ３−１によって導かれ、光カプラ４および光ファイバ５を介して被測定物１１に照射される（図１２ステップＳ３００）。
【０１０６】
音響センサ６は、被測定物１１から発生する光音響信号を検出し、増幅器７は、音響センサ８から出力された電気信号を増幅する。ロックインアンプ９は、増幅器７の出力に含まれる信号のうち、ファンクションジェネレータ８から出力される参照信号によって決まる周波数の測定信号を検出する。
【０１０７】
情報処理装置１０ａのジェネレータ制御部１００は、ファンクションジェネレータ８が発生する参照信号の周波数を変化させることにより、レーザドライバ２からレーザダイオード１−１に供給される駆動電流の周波数を変化させ、光変調周波数を漸次変化させると共に、ロックインアンプ１１が検出する測定信号の周波数を漸次変化させる光変調周波数掃引を行う（図１２ステップＳ３０１）。こうして、変調周波数と測定信号の振幅（音圧）との関係を取得し（図１２ステップＳ３０２）、測定信号の振幅が最大の変調周波数であるピーク変調周波数を探索する（図１２ステップＳ３０３）。
【０１０８】
次に、レーザダイオード１−１を再び動作させて連続光を出力させ（図１２ステップＳ３０４）、続いてレーザダイオード１−２（第２の光源）を動作させて連続光を出力させる（図１２ステップＳ３０５）。レーザダイオード１−１から放射される光の波長は例えば１３８２ｎｍ、レーザダイオード１−２から放射される光の波長（第２の波長）は例えば１６１０ｎｍである。
【０１０９】
情報処理装置１０ａのジェネレータ制御部１００は、ファンクションジェネレータ８が発生する２つの参照信号の周波数を変化させることにより、レーザドライバ２からレーザダイオード１−１，１−２に供給される駆動電流の周波数を変化させ、光変調周波数を上記のピーク変調周波数に設定すると共に、ロックインアンプ９が検出する測定信号の周波数をピーク変調周波数に設定する。このとき、レーザドライバ２から、同一周波数で逆位相の矩形波の駆動電流をレーザダイオード１−１，１−２に供給させることにより、レーザダイオード１−１，１−２から放射される光を同一周波数（ピーク変調周波数）で逆位相の信号によりそれぞれ強度変調する（図１２ステップＳ３０６）。
【０１１０】
レーザダイオード１−１，１−２から放射された強度変調光は、それぞれ光ファイバ３−１，３−２によって導かれ、光カプラ４によって合波され、さらに光ファイバ５によって導かれ、被測定物１１に照射される。
【０１１１】
情報処理装置１０ａの光強度制御部１０１は、ファンクションジェネレータ８を制御して、レーザダイオード１−１を駆動するための参照信号の電圧を所定値に設定することにより、レーザドライバ２からレーザダイオード１−１に供給される駆動電流の大きさを一定にし、レーザダイオード１−１から放射される光の強度を一定にする（図１２ステップＳ３０７）。
【０１１２】
また、光強度制御部１０１は、ファンクションジェネレータ８を制御して、レーザダイオード１−２を駆動するための参照信号の電圧を変化させることにより、レーザドライバ２からレーザダイオード１−２に供給される駆動電流の大きさを変化させ、レーザダイオード１−２から放射される光の強度を漸次変化させる光強度掃引を行い、ロックインアンプ９から出力される測定信号の振幅が最小になるようにレーザダイオード１−２から放射される光の強度を調節する（図１２ステップＳ３０８）。
【０１１３】
情報処理装置１０ａの光強度測定部１０２は、参考例と同様に、レーザダイオード１−２から放射される光の強度を測定する。そして、光強度測定部１０２は、測定信号の振幅が最小となるときの、レーザダイオード１−２から放射される光の強度Ｉ0を記憶部１０７に格納する（図１２ステップＳ３０９）。
【０１１４】
次に、情報処理装置１０ａの温度制御部１０８は、温度制御器１２を制御して被測定物１１の温度をΔＴ（例えば１℃）だけ上げる（図１２ステップＳ３１０）。
【０１１５】
情報処理装置１０ａの光強度制御部１０１は、ファンクションジェネレータ８を制御して、レーザダイオード１−２を駆動するための参照信号の電圧を変化させることにより、レーザドライバ２からレーザダイオード１−２に供給される駆動電流の大きさを変化させ、レーザダイオード１−２から放射される光の強度を漸次変化させる光強度掃引を行い、ロックインアンプ９から出力される測定信号の振幅が最小になるようにレーザダイオード１−２から放射される光の強度を調節する（図１２ステップＳ３１１）。上記のとおり、レーザダイオード１−１から放射される光の強度は一定である。
【０１１６】
情報処理装置１０ａの光強度測定部１０２は、レーザダイオード１−２から放射される光の強度をステップＳ３０９と同様に測定する。そして、光強度測定部１０２は、測定信号の振幅が最小となるときの、レーザダイオード１−２から放射される光の強度Ｉ1を記憶部１０７に格納する（図１２ステップＳ３１２）。
【０１１７】
情報処理装置１０ａの光強度変化量算出部１０３は、ステップＳ３０９，Ｓ３１２で求めた光強度Ｉ0，Ｉ1より光強度Ｉ0に対する光強度の変化量ΔＩ＝（Ｉ1−Ｉ0）の比ΔＩ／Ｉを、式（１５）により計算する（図１２ステップＳ３１３）。被測定物１１の温度を１℃上げた場合の温度変化前後の光強度の変化量の比ΔＩ／Ｉを図１３に示す。図１３によると、ΔＩ／Ｉ＝１．８％という結果が得られる。
【０１１８】
次に、情報処理装置１０ａの音響モード推定部１０４は、表２のような所定のデータセットを用いて、光強度の変化量の比ΔＩ2／Ｉ2を式（１６）、式（１７）により音響モードｑ毎に計算する（図１２ステップＳ３１４）。データセットは、第１、第２、第３の波長と、これらの波長の光に対する特定成分（例えばグルコース）の光吸収係数と、これら光吸収係数の温度依存性との関係を表すものである。
【０１１９】
【表２】
0006606018.tif 000021
【０１２０】
γ1は第１の波長の光に対する被測定物の光吸収係数α1の温度依存性（温度吸光度）、γ2は第２の波長の光に対する被測定物の光吸収係数α2の温度依存性である。関数Ｆ1は第１の波長の光に関する光吸収係数α1＝０．３０／ｍｍと音響モードｑから計算することができ、関数Ｆ2は第２の波長の光に関する光吸収係数α2＝０．２８／ｍｍと音響モードｑから計算することができる。被測定物１１の長さＬは既知の値である。
【０１２１】
ステップＳ３１４の計算の結果から、図１４に示すような音響モードｑと光強度の変化量の比ΔＩ2／Ｉ2との関係が得られる。情報処理装置１０ａの音響モード推定部１０４は、ステップＳ３１３で計算された光強度の変化量の比ΔＩ／Ｉと、ステップＳ３１４で計算した光強度の変化量の比ΔＩ2／Ｉ2とが等しくなる音響モードｑの値を、今回のピーク変調周波数での音響モードの値として決定する（図１２ステップＳ３１５）。本実施の形態の例では、ｑ＝１であることが推測できる。
【０１２２】
情報処理装置１０ａの補正係数決定部１０５は、音響モード推定部１０４が推定した音響モードｑの値に対応する、被測定物１１内の特定成分（本実施の形態の例ではグルコース）の濃度の測定感度Ｇｑを、濃度測定感度と音響モードｑとの既知の関係（図９）から求め、この濃度測定感度Ｇｑに対する音響モードｑ＝０のときの濃度測定感度Ｇ０の比Ｇ０／Ｇｑを補正係数Ａとして決定する（図１２ステップＳ３１６）。
【０１２３】
図９によると、音響モードｑ＝１のときの濃度測定感度Ｇｑは音響モードｑ＝０のときの濃度測定感度Ｇ０より約３倍大きいので、補正係数決定部１０５は、補正係数Ａ＝１／３と求めることができる。
【０１２４】
情報処理装置１０ａの濃度導出部１０６は、ステップＳ３１３で計算された光強度の変化量の比ΔＩ／Ｉ（音響モードｑ＝１のときの光強度の変化量の比ΔＩ2／Ｉ2）と、補正係数決定部１０５が決定した補正係数Ａとから、測定開始時から任意の時間経過後の時点における被測定物１１内の特定成分の濃度の変化量ΔＣ［％］を式（１８）により計算する（図１２ステップＳ３１７）。
【０１２５】
最後に、濃度導出部１０６は、ステップＳ３１７で計算した特定成分の濃度の変化量ΔＣ［％］と特定成分の既知の参照濃度とから、任意の時間経過後の時点における特定成分の濃度を計算する（図１２ステップＳ３１８）。
こうして、本実施の形態では、上記の参考例と同様の効果を得ることができる。
【０１２６】
音波の反射伝搬は、音響インピーダンスを用いて評価することができ、例えば異種材料界面における反射率Ｒおよび透過率Ｔは下式で表すことができる。
【０１２７】
【数１９】
0006606018.tif 000022
【０１２８】
【数２０】
0006606018.tif 000023
【０１２９】
Ｚ1は第１の材料の音響インピーダンス、Ｚ2は第１の材料と接する第２の材料の音響インピーダンスである。音響インピーダンスＺは式（２１）で表すことができる。
【０１３０】
【数２１】
0006606018.tif 000024
【０１３１】
ρは材料の密度、Ｃは材料内での音速である。音響インピーダンスの差によって、反射、透過する音波が決定される。図１５に音波の反射率と音響インピーダンスとの関係を示す。被測定物１１と共振器材料界面との反射率が０．８以上であることが共振器（被測定物１１を収容する光音響セル）として機能することの必要条件である。この条件を満たす材料として、たとえばガラスや鉄や銅をはじめとした各種金属がある。照射する光により生じる音波の周波数は、式（２２）のように記述することができる。共振器は、照射される光の光軸に対して略直交する面を２つ持ち、これらの２つの面は互いに平行平面となっていることが好ましい。
【０１３２】
【数２２】
0006606018.tif 000025
【０１３３】
λは音波の波長、Ｃは被測定物１１内での音速、ｆは光の変調周波数である。本発明の光音響を用いた測定方法においては、光の変調周波数を１００ｋＨｚ−１ＭＨｚとした場合、この時の音波の波長は１ｍｍ弱から１２ｍｍとなる。たとえば、軸方向２次のモードを得るためには、共振長はおよそ波長程度に設定する必要がある。
【０１３４】
なお、参考例および実施の形態において、被測定物１１に照射される複数の光は、これら複数の光の波長に対する、被測定物１１に含まれる測定対象の成分以外の成分の光吸収係数が略等しいことが好ましい。
【０１３５】
また、参考例および実施の形態において、異なる波長の複数の光を光学系（不図示）によりそれぞれ平行光にして被測定物１１の略同一領域に同軸で照射することが好ましい。また、このとき異なる波長の複数の光のビーム径は略等しいことが好ましい。
【０１３６】
参考例および実施の形態の情報処理装置１０，１０ａは、例えばＣＰＵ（Central Processing Unit）、記憶装置およびインタフェースを備えたコンピュータとこれらのハードウェア資源を制御するプログラムによって実現することができる。ＣＰＵは、記憶装置に格納されたプログラムに従って参考例および実施の形態で説明した処理を実行する。
【産業上の利用可能性】
【０１３７】
本発明は、血液グルコース、アルブミン等の成分の濃度をモニターする技術に適用することができる。
【符号の説明】
【０１３８】
１−１〜１−３…レーザダイオード、２…レーザドライバ、３−１〜３−３，５…光ファイバ、４…光カプラ、６…音響センサ、７…増幅器、８…ファンクションジェネレータ、９…ロックインアンプ、１０，１０ａ…情報処理装置、１１…被測定物、１２…温度制御器、１００…ジェネレータ制御部、１０１…光強度制御部、１０２…光強度測定部、１０３…光強度変化量算出部、１０４…音響モード推定部、１０５…補正係数決定部、１０６…濃度導出部、１０７…記憶部、１０８…温度制御部。
【図１】
0006606018.tif 000026
【図２】
0006606018.tif 000027
【図３】
0006606018.tif 000028
【図４】
0006606018.tif 000029
【図５】
0006606018.tif 000030
【図６】
0006606018.tif 000031
【図７】
0006606018.tif 000032
【図８】
0006606018.tif 000033
【図９】
0006606018.tif 000034
【図１０】
0006606018.tif 000035
【図１１】
0006606018.tif 000036
【図１２】
0006606018.tif 000037
【図１３】
0006606018.tif 000038
【図１４】
0006606018.tif 000039
【図１５】
0006606018.tif 000040
【図１６】
0006606018.tif 000041
【図１７】
0006606018.tif 000042
【図１８】
0006606018.tif 000043
【図１９】
0006606018.tif 000044