JP2017217202A - Device and method for measuring constituent concentration - Google Patents

Device and method for measuring constituent concentration Download PDF

Info

Publication number
JP2017217202A
JP2017217202A JP2016113652A JP2016113652A JP2017217202A JP 2017217202 A JP2017217202 A JP 2017217202A JP 2016113652 A JP2016113652 A JP 2016113652A JP 2016113652 A JP2016113652 A JP 2016113652A JP 2017217202 A JP2017217202 A JP 2017217202A
Authority
JP
Japan
Prior art keywords
light
intensity
measured
wavelength
measurement
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2016113652A
Other languages
Japanese (ja)
Other versions
JP6606018B2 (en
Inventor
雄次郎 田中
Yujiro Tanaka
雄次郎 田中
倫子 瀬山
Michiko Seyama
倫子 瀬山
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nippon Telegraph and Telephone Corp
Original Assignee
Nippon Telegraph and Telephone Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nippon Telegraph and Telephone Corp filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2016113652A priority Critical patent/JP6606018B2/en
Publication of JP2017217202A publication Critical patent/JP2017217202A/en
Application granted granted Critical
Publication of JP6606018B2 publication Critical patent/JP6606018B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

PROBLEM TO BE SOLVED: To compensate acoustic mode dependency of concentration measurement sensitivity to a measurement target constituent included in a measurement target object.SOLUTION: Two of laser diodes 1-1 to 1-3 emit light with two waves different in wavelength while modulating intensity of the light by signals with the same frequency and also a different phase. An acoustic sensor 6 detects photoacoustic signals generated from a measurement target object 11. An information processing device 10 estimates acoustic modes with respect to two time points different in wavelength of light changing the intensity or two time points different in temperature of the measurement target object 11 respectively on the basis of a change amount of optical intensity acquired from a result of measurement of the optical intensity when intensity of the measured signal gets lowest and a change amount of optical intensity acquired from a data set, and determines a correction coefficient of concentration measurement sensitivity to a measurement target constituent on the basis of the estimated acoustic mode and the known relationship between the concentration measurement sensitivity to the measurement target constituent included in the measurement target object 11 and the acoustic mode.SELECTED DRAWING: Figure 3

Description

本発明は、成分濃度測定装置、特に、血液成分、例えばグルコースを測定対象とした測定装置に関するものである。   The present invention relates to a component concentration measuring apparatus, and more particularly to a measuring apparatus that uses a blood component such as glucose as a measurement target.

高齢化により、生活習慣病に対する対応が大きな課題になりつつあり、血液の成分に対する注目が集まっている。しかし、現在までに開発された技術は、採取したサンプルを試薬と反応させたり、専用のセンサに吸着させたりといった1つのセンサで1つの成分を検出するようになっており、マルチなセンシング技術が注目されている。そこで、電磁波を照射したときに、測定対象とする血液成分(例えば、血糖値の場合はグルコース分子)に吸収されることを利用して、透過あるいは反射した電磁波を測定する方法が行われている。   With the aging of society, dealing with lifestyle-related diseases is becoming a major issue, and attention is being paid to blood components. However, the technology that has been developed so far is to detect one component with a single sensor, such as reacting a collected sample with a reagent or adsorbing it to a dedicated sensor. Attention has been paid. Therefore, there is a method of measuring transmitted or reflected electromagnetic waves by utilizing absorption by a blood component to be measured (for example, glucose molecules in the case of blood glucose level) when irradiated with electromagnetic waves. .

しかし、グルコースと電磁波との相互作用は小さく、また生体に安全に照射し得る電磁波の強度には制限があり、さらに生体が電磁波に対して散乱体であるために、生体の血糖値測定においては、十分な効果を挙げるに至っていない。
そこで、電磁波を照射し局所的に加熱して熱膨張を起こして発生する音波を観測する、光音響法が注目されている。
However, the interaction between glucose and electromagnetic waves is small, and there is a limit to the intensity of electromagnetic waves that can be safely irradiated to a living body. Furthermore, since living bodies are scatterers against electromagnetic waves, , Has not been effective enough.
Therefore, a photoacoustic method has been attracting attention, in which a sound wave generated by irradiating an electromagnetic wave and locally heating it to cause thermal expansion is observed.

光音響法には、パルス(pulse)法と連続波(continuous-wave、以下CWとする)法の二つの方式がある。しかし、従来のパルス法やCW法では、数回にわたる血漿中のグルコース濃度測定中に、グルコース濃度以外の他の血漿中パラメータ(例えば体温や、他の成分の濃度等)も変わる可能性が高いので、グルコース選択性が悪く、正確なグルコース濃度を得ることが難しいという問題点があった。   There are two types of photoacoustic methods: a pulse method and a continuous-wave (hereinafter referred to as CW) method. However, in the conventional pulse method and CW method, during the measurement of the glucose concentration in plasma several times, there is a high possibility that other plasma parameters other than the glucose concentration (for example, body temperature, concentration of other components, etc.) will also change. Therefore, there is a problem that glucose selectivity is poor and it is difficult to obtain an accurate glucose concentration.

そこで、CW法に属する方法として、カムー セルジュ(S.Camou)らにより光パワーバランスシフト(Optical power balance shift:OPBS)法が開発されている(特許文献1−3参照)。   Therefore, as a method belonging to the CW method, an optical power balance shift (OPPBS) method has been developed by S. Camou et al. (See Patent Documents 1-3).

このOPBS法を実施するための従来の成分濃度測定装置の構成例を図16に示す。成分濃度測定装置は、レーザ光を放射するレーザダイオード1−1,1−2と、レーザダイオード1−1,1−2を駆動するレーザドライバ2と、レーザダイオード1−1,1−2から放射されたレーザ光を導く光ファイバ3−1,3−2と、レーザダイオード1−1,1−2から放射されたレーザ光を合波する光カプラ4と、光カプラ4によって合波されたレーザ光を導く光ファイバ5と、光音響効果によって被測定物11から発生する光音響信号を検出し、音圧に比例した電気信号に変換する音響センサ6と、音響センサ6から出力された電気信号を増幅する増幅器7と、参照信号を発生するファンクションジェネレータ8と、増幅器7の出力信号とファンクションジェネレータ8から出力された参照信号とを入力として、増幅器7の出力信号から所望の周波数の測定信号を検出するロックインアンプ9と、ファンクションジェネレータ8およびロックインアンプ9を制御すると共に、ロックインアンプ9が検出した測定信号を処理して特定の成分濃度を導出するコンピュータからなる情報処理装置10とから構成される。   A configuration example of a conventional component concentration measuring apparatus for carrying out this OPBS method is shown in FIG. The component concentration measurement apparatus includes laser diodes 1-1 and 1-2 that emit laser light, a laser driver 2 that drives laser diodes 1-1 and 1-2, and laser diodes 1-1 and 1-2. Optical fibers 3-1 and 3-2 for guiding the laser light emitted, an optical coupler 4 for combining the laser light emitted from the laser diodes 1-1 and 1-2, and a laser combined by the optical coupler 4. An optical fiber 5 that guides light, an acoustic sensor 6 that detects a photoacoustic signal generated from the object 11 to be measured by the photoacoustic effect, and converts it into an electrical signal proportional to sound pressure, and an electrical signal output from the acoustic sensor 6 The amplifier 7 for amplifying the signal, the function generator 8 for generating a reference signal, the output signal of the amplifier 7 and the reference signal output from the function generator 8 as inputs, and the output of the amplifier 7 The lock-in amplifier 9 that detects a measurement signal of a desired frequency from the signal, the function generator 8 and the lock-in amplifier 9 are controlled, and the measurement signal detected by the lock-in amplifier 9 is processed to derive a specific component concentration. And an information processing apparatus 10 including a computer.

図16に示した成分濃度測定装置を用いた従来のOPBS法では、まずレーザドライバ2は、ファンクションジェネレータ8から出力される参照信号に応じて矩形波の駆動電流をレーザダイオード1−1(第1の光源)に供給することにより、レーザダイオード1−1から放射される光を強度変調する。レーザダイオード1−1から放射された強度変調光は、光ファイバ3−1によって導かれ、光カプラ4および光ファイバ5を介して被測定物11に照射される(図17ステップS100)。   In the conventional OPBS method using the component concentration measuring apparatus shown in FIG. 16, first, the laser driver 2 generates a rectangular-wave drive current according to the reference signal output from the function generator 8 to the laser diode 1-1 (first The light emitted from the laser diode 1-1 is intensity-modulated. The intensity-modulated light emitted from the laser diode 1-1 is guided by the optical fiber 3-1, and is irradiated to the device under test 11 through the optical coupler 4 and the optical fiber 5 (step S100 in FIG. 17).

音響センサ6は、被測定物11から発生する光音響信号を検出し、増幅器7は、音響センサ8から出力された電気信号を増幅する。ロックインアンプ9は、増幅器7の出力に含まれる信号のうち、ファンクションジェネレータ8から出力される参照信号によって決まる周波数の測定信号を検出する。   The acoustic sensor 6 detects a photoacoustic signal generated from the device under test 11, and the amplifier 7 amplifies the electrical signal output from the acoustic sensor 8. The lock-in amplifier 9 detects a measurement signal having a frequency determined by the reference signal output from the function generator 8 among the signals included in the output of the amplifier 7.

情報処理装置10は、ファンクションジェネレータ8が発生する参照信号の周波数を変化させることにより、レーザドライバ2からレーザダイオード1−1に供給される駆動電流の周波数を変化させ、光変調周波数を漸次変化させると共に、ロックインアンプ9が検出する測定信号の周波数(光変調周波数と同一の周波数)を漸次変化させる光変調周波数掃引を行う(図17ステップS101)。こうして、変調周波数と測定信号の振幅(音圧)との関係を取得し(図17ステップS102)、測定信号の振幅が最大の変調周波数であるピーク変調周波数を探索する(図17ステップS103)。   The information processing apparatus 10 changes the frequency of the reference signal generated by the function generator 8, thereby changing the frequency of the drive current supplied from the laser driver 2 to the laser diode 1-1, and gradually changing the light modulation frequency. At the same time, an optical modulation frequency sweep is performed to gradually change the frequency of the measurement signal detected by the lock-in amplifier 9 (the same frequency as the optical modulation frequency) (step S101 in FIG. 17). In this way, the relationship between the modulation frequency and the amplitude (sound pressure) of the measurement signal is acquired (step S102 in FIG. 17), and the peak modulation frequency that is the modulation frequency with the maximum amplitude of the measurement signal is searched (step S103 in FIG. 17).

次に、レーザダイオード1−1を再び動作させて光を出力させ(図17ステップS104)、続いてレーザダイオード1−2(第2の光源)を動作させて光を出力させる(図17ステップS105)。2つのレーザダイオード1−1,1−2から放射される光の波長は異なる。   Next, the laser diode 1-1 is operated again to output light (step S104 in FIG. 17), and then the laser diode 1-2 (second light source) is operated to output light (step S105 in FIG. 17). ). The wavelengths of light emitted from the two laser diodes 1-1 and 1-2 are different.

情報処理装置10は、ファンクションジェネレータ8が発生する2つの参照信号の周波数を変化させることにより、レーザドライバ2からレーザダイオード1−1,1−2に供給される駆動電流の周波数を変化させ、光変調周波数を上記のピーク変調周波数に設定すると共に、ロックインアンプ9が検出する測定信号の周波数をピーク変調周波数に設定する。このとき、レーザドライバ2から、同一周波数で逆位相の矩形波の駆動電流をレーザダイオード1−1,1−2に供給させることにより、レーザダイオード1−1,1−2から放射される光を同一周波数(ピーク変調周波数)で逆位相の信号によりそれぞれ強度変調する(図17ステップS106)。   The information processing apparatus 10 changes the frequency of the drive current supplied from the laser driver 2 to the laser diodes 1-1 and 1-2 by changing the frequency of the two reference signals generated by the function generator 8. The modulation frequency is set to the above peak modulation frequency, and the frequency of the measurement signal detected by the lock-in amplifier 9 is set to the peak modulation frequency. At this time, light emitted from the laser diodes 1-1 and 1-2 is supplied from the laser driver 2 to the laser diodes 1-1 and 1-2 by supplying a rectangular-wave drive current of the same frequency and opposite phase. Intensity modulation is performed with signals having the same frequency (peak modulation frequency) and opposite phases (step S106 in FIG. 17).

レーザダイオード1−1,1−2から放射された強度変調光は、それぞれ光ファイバ3−1,3−2によって導かれ、光カプラ4によって合波され、さらに光ファイバ5によって導かれ、被測定物11に照射される。情報処理装置10は、レーザドライバ2からレーザダイオード1−1に供給される駆動電流の大きさを所定値にすることにより、レーザダイオード1−1から放射される光の強度を固定し(図17ステップS107)、レーザドライバ2からレーザダイオード1−2に供給される駆動電流の大きさを変化させることにより、レーザダイオード1−2から放射される光の強度を漸次変化させる光強度掃引を行い、測定信号の振幅が最小になるようにレーザダイオード1−2から放射される光の強度を調節する(図17ステップS108)。   The intensity-modulated lights emitted from the laser diodes 1-1 and 1-2 are respectively guided by the optical fibers 3-1 and 3-2, combined by the optical coupler 4, and further guided by the optical fiber 5 to be measured. The object 11 is irradiated. The information processing apparatus 10 fixes the intensity of light emitted from the laser diode 1-1 by setting the magnitude of the drive current supplied from the laser driver 2 to the laser diode 1-1 to a predetermined value (FIG. 17). Step S107), by changing the magnitude of the drive current supplied from the laser driver 2 to the laser diode 1-2, performing a light intensity sweep that gradually changes the intensity of the light emitted from the laser diode 1-2, The intensity of light emitted from the laser diode 1-2 is adjusted so that the amplitude of the measurement signal is minimized (step S108 in FIG. 17).

そして、一定時間経過後(図17ステップS109)、再びレーザダイオード1−1,1−2を動作させ、レーザダイオード1−1から放射される光の強度を固定し(図17ステップS110)、測定信号の振幅が最小になるようにレーザダイオード1−2から放射される光の強度を調節する(図17ステップS111)。そして、情報処理装置10は、ステップS108で得られた光強度とステップS111で得られた光強度から、光強度の変化量を計算する(図17ステップS112)。   Then, after a predetermined time has elapsed (step S109 in FIG. 17), the laser diodes 1-1 and 1-2 are operated again, the intensity of light emitted from the laser diode 1-1 is fixed (step S110 in FIG. 17), and measurement is performed. The intensity of light emitted from the laser diode 1-2 is adjusted so that the amplitude of the signal is minimized (step S111 in FIG. 17). Then, the information processing apparatus 10 calculates the amount of change in light intensity from the light intensity obtained in step S108 and the light intensity obtained in step S111 (step S112 in FIG. 17).

このように、従来のOPBS法では、光波長が異なり位相差が180°の2つの光ビームを被測定物に対して、同一光出力口から照射し、2つの光ビームの強度を増減させながら、光音響信号の振幅が最小の箇所の位相の変曲点を探し、この結果から血液中に溶解している分子濃度を求める。具体的には、2つの光ビームのうち一方の光ビームの強度を変えながら、光音響信号の振幅が最小となる光強度を探すことで、光強度の変化量より被測定物内の特定成分(例えばグルコース)の濃度の正確な測定を行う。   As described above, in the conventional OPBS method, two light beams having different light wavelengths and a phase difference of 180 ° are irradiated to the object to be measured from the same light output port, and the intensity of the two light beams is increased or decreased. Then, the phase inflection point of the portion where the amplitude of the photoacoustic signal is minimum is searched, and the concentration of the molecule dissolved in the blood is obtained from the result. Specifically, by searching for the light intensity at which the amplitude of the photoacoustic signal is minimized while changing the intensity of one of the two light beams, the specific component in the object to be measured is determined from the amount of change in the light intensity. Make an accurate measurement of the concentration of (eg glucose).

特開2014−50563号公報JP 2014-50563 A 特開2013−106874号公報JP 2013-106874 A 特開2012−179212号公報JP 2012-179212 A

被測定物を収容するケースが円筒形、すなわち光音響エネルギーが閉じ込められる空間が円筒空洞である場合の、この円筒空洞内での共振を例にとった音波の共振モードを図18に示す。ここで、Rは円筒の半径、Lは円筒の長さ、jは半径方向(Radial)のモードの番号、mは円周方向(Azimuthal)のモードの番号、qは長手方向(Longitudinal)のモードの番号である。   FIG. 18 shows a resonance mode of a sound wave taking a resonance in this cylindrical cavity as an example when the case for accommodating the object to be measured is cylindrical, that is, the space in which the photoacoustic energy is confined is a cylindrical cavity. Where R is the radius of the cylinder, L is the length of the cylinder, j is the number of the radial mode, m is the number of the circumferential mode, and q is the longitudinal mode. Is the number.

変調周波数と音圧との関係を図19に示す。図19中の(101)、(200)、(102)、(300)は共振モード(jmq)を示しており、例えば(101)は半径方向のモードの番号jが1、円周方向のモードの番号mが0、長手方向のモードの番号qが1であることを示している。このように、円筒形の共振器内では、音波は様々なモードで共振し、それらが複雑に干渉し合った音波の分布が生じている。   FIG. 19 shows the relationship between the modulation frequency and the sound pressure. In FIG. 19, (101), (200), (102), and (300) indicate resonance modes (jmq). For example, (101) is a mode number j in the radial direction, and a mode in the circumferential direction. This indicates that the number m is 0 and the mode number q in the longitudinal direction is 1. In this manner, in the cylindrical resonator, the sound wave resonates in various modes, and a sound wave distribution in which they interfere with each other in a complicated manner is generated.

CW光を用いた従来の光音響法においては、光の照射により生じる音波の分布と定在波の音波の干渉によって測定値が大きく変わってしまうという問題があった。すなわち、被測定物に含まれる対象の成分の濃度の検量線の傾きが数倍のオーダーで変わってしまうという問題があった。   In the conventional photoacoustic method using CW light, there is a problem that the measurement value greatly changes due to the distribution of the sound wave generated by the light irradiation and the interference of the sound wave of the standing wave. That is, there is a problem that the slope of the calibration curve of the concentration of the target component contained in the object to be measured changes on the order of several times.

本発明は、上記課題を解決するためになされたもので、被測定物に含まれる測定対象の成分に対する濃度測定感度の音響モード依存性を補償することができる成分濃度測定装置および方法を提供することを目的とする。   The present invention has been made to solve the above problems, and provides a component concentration measuring apparatus and method that can compensate the acoustic mode dependence of the concentration measuring sensitivity for the component to be measured included in the object to be measured. For the purpose.

本発明の成分濃度測定装置は、互いに異なる波長の複数の光を同一の周波数で且つ異なる位相の信号によりそれぞれ強度変調して被測定物に照射する光照射手段と、前記複数の強度変調光のうち少なくとも1つの光の強度を変化させる光強度制御手段と、光照射によって前記被測定物から発生する光音響信号を検出して電気信号を出力する光音響信号検出手段と、前記光強度制御手段が変化させた光の強度を測定する光強度測定手段と、強度を変化させる光の波長が異なる2つの時点または前記被測定物の温度が異なる2つの時点について、それぞれ前記電気信号の強度が最低となったときの前記光強度測定手段の測定結果から得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記光音響信号の共振モードである音響モードを推定する音響モード推定手段と、この音響モード推定手段が推定した音響モードと、前記被測定物に含まれる測定対象の成分に対する濃度測定感度と音響モードとの既知の関係に基づいて、前記測定対象の成分に対する濃度測定感度の補正係数を決定する補正係数決定手段と、前記光強度測定手段の測定結果から得られる光強度の変化量と前記補正係数に基づいて前記測定対象の成分の濃度を導出する濃度導出手段とを備えることを特徴とするものである。   The component concentration measuring apparatus according to the present invention includes a light irradiating unit that irradiates a measured object with a plurality of lights having different wavelengths by using signals having the same frequency and different phases, and irradiating the object to be measured. A light intensity control means for changing the intensity of at least one of the light; a photoacoustic signal detection means for detecting a photoacoustic signal generated from the object to be measured by light irradiation and outputting an electrical signal; and the light intensity control means. The intensity of the electrical signal is the lowest at two time points at which the wavelength of light whose intensity is changed or two time points at which the temperature of the object to be measured is different. Based on the change amount of the light intensity obtained from the measurement result of the light intensity measurement means and the change amount of the light intensity obtained from the known data set, the resonance of the photoacoustic signal An acoustic mode estimation means for estimating an acoustic mode that is a mode, an acoustic mode estimated by the acoustic mode estimation means, and a known relationship between a density measurement sensitivity and an acoustic mode for a component to be measured included in the object to be measured Correction coefficient determining means for determining a correction coefficient of concentration measurement sensitivity for the component to be measured, based on the measurement result based on the amount of change in light intensity obtained from the measurement result of the light intensity measuring means and the correction coefficient And a concentration deriving unit for deriving the concentration of the target component.

また、本発明の成分濃度測定装置の1構成例において、前記光照射手段は、異なる2つの波長の光を前記被測定物に同時に照射し、前記音響モード推定手段は、第1の波長の光の強度を一定にして第2の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第2の波長の光の強度の測定結果と、前記第1の波長の光の強度を一定にして第3の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第3の波長の光の強度の測定結果とから得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記音響モードを推定することを特徴とするものである。   Moreover, in one structural example of the component concentration measuring apparatus of this invention, the said light irradiation means irradiates the said to-be-measured object light of two different wavelengths simultaneously, and the said acoustic mode estimation means has the light of 1st wavelength. The measurement result of the light intensity of the second wavelength when the intensity of the electric signal becomes the minimum by changing the intensity of the light of the second wavelength while keeping the intensity of The intensity of light obtained from the measurement result of the intensity of the light of the third wavelength when the intensity of the electric signal is minimized by changing the intensity of the light of the third wavelength while keeping the intensity of light constant. The acoustic mode is estimated based on the amount of change and the amount of change in light intensity obtained from a known data set.

また、本発明の成分濃度測定装置の1構成例は、さらに、前記被測定物の温度を制御する温度制御手段を備え、前記光照射手段は、異なる2つの波長の光を前記被測定物に同時に照射し、前記音響モード推定手段は、第1の波長の光の強度を一定にして第2の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第2の波長の光の強度の測定結果と、前記温度制御手段により前記被測定物の温度を変化させた後に前記第1の波長の光の強度を一定にして前記第2の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第2の波長の光の強度の測定結果とから得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記音響モードを推定することを特徴とするものである。   Moreover, one structural example of the component concentration measuring apparatus of this invention is further equipped with the temperature control means which controls the temperature of the said to-be-measured object, and the said light irradiation means makes the said to-be-measured object light of two different wavelengths. Irradiating at the same time, the acoustic mode estimating means changes the intensity of the light of the second wavelength while keeping the intensity of the light of the first wavelength constant, so that the intensity of the electric signal becomes the minimum. The measurement result of the light intensity of the second wavelength and the intensity of the light of the second wavelength with the light intensity of the first wavelength constant after changing the temperature of the object to be measured by the temperature control means. The amount of change in the light intensity obtained from the measurement result of the intensity of the light of the second wavelength when the intensity of the electric signal is minimized by changing, and the amount of change in the light intensity obtained from a known data set Based on the above, the acoustic mode is estimated A.

また、本発明の成分濃度測定装置の1構成例において、前記被測定物に照射される複数の光は、これら複数の光の波長に対する、前記被測定物に含まれる測定対象の成分以外の成分の光吸収係数が略等しいことを特徴とするものである。
また、本発明の成分濃度測定装置の1構成例において、前記光照射手段は、異なる波長の複数の光をそれぞれ平行光にして前記被測定物の略同一領域に同軸で照射し、異なる波長の複数の光のビーム径は略等しいことを特徴とするものである。
また、本発明の成分濃度測定装置の1構成例は、さらに、前記被測定物を収容し、前記光音響信号を増幅するための共振器を備えることを特徴とするものである。
また、本発明の成分濃度測定装置の1構成例において、前記共振器は、前記光照射手段から照射される光の光軸に対して略直交する面を2つ持ち、これらの面は互いに平行平面となっていることを特徴とするものである。
Moreover, in one structural example of the component concentration measuring apparatus of this invention, the some light irradiated to the said to-be-measured object is components other than the component of the measuring object contained in the to-be-measured object with respect to the wavelength of these some light Are characterized in that their light absorption coefficients are substantially equal.
Further, in one configuration example of the component concentration measuring apparatus of the present invention, the light irradiating means irradiates a plurality of lights having different wavelengths as parallel lights, coaxially irradiating substantially the same region of the object to be measured, and having different wavelengths. The plurality of light beams have substantially the same beam diameter.
Further, one configuration example of the component concentration measuring apparatus of the present invention is characterized by further comprising a resonator for accommodating the object to be measured and amplifying the photoacoustic signal.
Moreover, in one configuration example of the component concentration measuring apparatus of the present invention, the resonator has two surfaces substantially orthogonal to the optical axis of the light emitted from the light irradiation means, and these surfaces are parallel to each other. It is a flat surface.

また、本発明の成分濃度測定方法は、互いに異なる波長の複数の光を同一の周波数で且つ異なる位相の信号によりそれぞれ強度変調して被測定物に照射する光照射ステップと、前記複数の強度変調光のうち少なくとも1つの光の強度を変化させる光強度制御ステップと、光照射によって前記被測定物から発生する光音響信号を検出して電気信号を出力する光音響信号検出ステップと、前記光強度制御ステップで変化させた光の強度を測定する光強度測定ステップと、強度を変化させる光の波長が異なる2つの時点または前記被測定物の温度が異なる2つの時点について、それぞれ前記電気信号の強度が最低となったときの前記光強度測定ステップの測定結果から得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記光音響信号の共振モードである音響モードを推定する音響モード推定ステップと、この音響モード推定ステップで推定した音響モードと、前記被測定物に含まれる測定対象の成分に対する濃度測定感度と音響モードとの既知の関係に基づいて、前記測定対象の成分に対する濃度測定感度の補正係数を決定する補正係数決定ステップと、前記光強度測定ステップの測定結果から得られる光強度の変化量と前記補正係数に基づいて前記測定対象の成分の濃度を導出する濃度導出ステップとを含むことを特徴とするものである。   The component concentration measurement method of the present invention includes a light irradiation step of irradiating a measured object with a plurality of lights having different wavelengths by using signals having the same frequency and different phases, and the plurality of intensity modulations. A light intensity control step for changing the intensity of at least one of the light; a photoacoustic signal detection step for detecting a photoacoustic signal generated from the object to be measured by light irradiation and outputting an electrical signal; and the light intensity. The intensity of the electric signal at a light intensity measuring step for measuring the intensity of the light changed in the control step and at two time points at which the wavelength of the light at which the intensity is changed is different or at two time points at which the temperature of the object to be measured is different. The amount of change in light intensity obtained from the measurement result of the light intensity measurement step when the value becomes the minimum, and the amount of change in light intensity obtained from a known data set. Therefore, an acoustic mode estimation step for estimating an acoustic mode that is a resonance mode of the photoacoustic signal, an acoustic mode estimated in the acoustic mode estimation step, and a concentration measurement sensitivity for a component to be measured included in the object to be measured A correction coefficient determining step for determining a correction coefficient of concentration measurement sensitivity for the component to be measured based on a known relationship between the sound mode and the acoustic mode; and a change amount of the light intensity obtained from the measurement result of the light intensity measuring step; And a concentration derivation step for deriving the concentration of the component to be measured based on the correction coefficient.

本発明によれば、互いに異なる波長の複数の光を同一の周波数で且つ異なる位相の信号によりそれぞれ強度変調して被測定物に照射する光照射手段と、複数の強度変調光のうち少なくとも1つの光の強度を変化させる光強度制御手段と、光照射によって被測定物から発生する光音響信号を検出して電気信号を出力する光音響信号検出手段と、光強度制御手段が変化させた光の強度を測定する光強度測定手段と、強度を変化させる光の波長が異なる2つの時点または被測定物の温度が異なる2つの時点について、それぞれ電気信号の強度が最低となったときの光強度測定手段の測定結果から得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、光音響信号の共振モードである音響モードを推定する音響モード推定手段と、音響モード推定手段が推定した音響モードと、被測定物に含まれる測定対象の成分に対する濃度測定感度と音響モードとの既知の関係に基づいて、測定対象の成分に対する濃度測定感度の補正係数を決定する補正係数決定手段とを設けることにより、被測定物内の測定対象の成分に対する濃度測定感度の、音響モードによる変化を補償することができ、測定感度を向上させることができる。   According to the present invention, the light irradiating means for irradiating the object to be measured by intensity-modulating a plurality of lights having different wavelengths with signals having the same frequency and different phases, and at least one of the plurality of intensity-modulated lights. A light intensity control means for changing the light intensity, a photoacoustic signal detection means for detecting a photoacoustic signal generated from the object to be measured by light irradiation and outputting an electrical signal, and a light intensity control means for changing the light intensity Light intensity measurement means for measuring intensity, and light intensity measurement at the time when the intensity of the electric signal is the lowest at two time points where the wavelength of light changing the intensity is different or two time points where the temperature of the object to be measured is different An acoustic mode that estimates the acoustic mode, which is the resonance mode of the photoacoustic signal, based on the amount of change in light intensity obtained from the measurement results of the means and the amount of change in light intensity obtained from a known data set. Based on the known relationship between the estimation means, the acoustic mode estimated by the acoustic mode estimation means, and the concentration measurement sensitivity and the acoustic mode for the measurement target component contained in the object to be measured, the concentration measurement sensitivity of the measurement target component By providing the correction coefficient determining means for determining the correction coefficient, it is possible to compensate for the change in the concentration measurement sensitivity with respect to the component to be measured in the object to be measured due to the acoustic mode, and to improve the measurement sensitivity.

また、本発明では、被測定物を収容する共振器を設けることにより、光照射によって被測定物から発生する光音響信号を増幅することができる。   In the present invention, a photoacoustic signal generated from an object to be measured by light irradiation can be amplified by providing a resonator that accommodates the object to be measured.

また、本発明では、第1の波長の光の強度を一定にして第2の波長の光の強度を変化させて電気信号の強度が最低となったときの第2の波長の光の強度の測定結果と、第1の波長の光の強度を一定にして第3の波長の光の強度を変化させて電気信号の強度が最低となったときの第3の波長の光の強度の測定結果とから得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、音響モードを推定することができる。   In the present invention, the intensity of the light of the second wavelength when the intensity of the electric signal becomes the minimum by changing the intensity of the light of the second wavelength while keeping the intensity of the light of the first wavelength constant. The measurement result and the measurement result of the intensity of the light of the third wavelength when the intensity of the light of the third wavelength is changed and the intensity of the light of the third wavelength is changed to make the intensity of the electric signal minimum. The acoustic mode can be estimated based on the change amount of the light intensity obtained from the above and the change amount of the light intensity obtained from the known data set.

また、本発明では、第1の波長の光の強度を一定にして第2の波長の光の強度を変化させて電気信号の強度が最低となったときの第2の波長の光の強度の測定結果と、温度制御手段により被測定物の温度を変化させた後に第1の波長の光の強度を一定にして第2の波長の光の強度を変化させて電気信号の強度が最低となったときの第2の波長の光の強度の測定結果とから得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、音響モードを推定することができる。   In the present invention, the intensity of the light of the second wavelength when the intensity of the electric signal becomes the minimum by changing the intensity of the light of the second wavelength while keeping the intensity of the light of the first wavelength constant. After the measurement result and the temperature of the object to be measured are changed by the temperature control means, the intensity of the light of the first wavelength is made constant by changing the intensity of the light of the second wavelength, and the intensity of the electric signal becomes the minimum. The acoustic mode can be estimated based on the change amount of the light intensity obtained from the measurement result of the light intensity of the second wavelength at that time and the change amount of the light intensity obtained from the known data set.

また、本発明では、光照射手段が、異なる波長の複数の光をそれぞれ平行光にして被測定物の略同一領域に同軸で照射することにより、音波の干渉の影響を低減することができる。   Further, in the present invention, the light irradiating means converts a plurality of lights having different wavelengths into parallel light and irradiates substantially the same region of the object to be measured coaxially, thereby reducing the influence of the interference of the sound wave.

グルコース濃度の応答特性を示す図である。It is a figure which shows the response characteristic of glucose concentration. 変調周波数と音波強度の関係を示す図である。It is a figure which shows the relationship between a modulation frequency and a sound wave intensity. 本発明の第1の実施の形態に係る成分濃度測定装置の構成を示すブロック図である。It is a block diagram which shows the structure of the component concentration measuring apparatus which concerns on the 1st Embodiment of this invention. 本発明の第1の実施の形態に係る成分濃度測定装置の情報処理装置の構成を示すブロック図である。It is a block diagram which shows the structure of the information processing apparatus of the component concentration measuring apparatus which concerns on the 1st Embodiment of this invention. 本発明の第1の実施の形態に係る成分濃度測定装置の動作を説明するフローチャートである。It is a flowchart explaining operation | movement of the component concentration measuring apparatus which concerns on the 1st Embodiment of this invention. 本発明の第1の実施の形態における変調周波数と音波強度との関係の1例を示す図である。It is a figure which shows one example of the relationship between the modulation frequency and sound wave intensity in the 1st Embodiment of this invention. 本発明の第1の実施の形態における光強度変化の例を示す図である。It is a figure which shows the example of the light intensity change in the 1st Embodiment of this invention. 本発明の第1の実施の形態における音響モードと光強度の変化量の比との関係の1例を示す図である。It is a figure which shows one example of the relationship between the acoustic mode in the 1st Embodiment of this invention, and the ratio of the variation | change_quantity of light intensity. 本発明の第1の実施の形態におけるグルコース濃度感度と音響モードとの関係の1例を示す図である。It is a figure which shows one example of the relationship between the glucose concentration sensitivity in 1st Embodiment of this invention, and an acoustic mode. 本発明の第2の実施の形態に係る成分濃度測定装置の構成を示すブロック図である。It is a block diagram which shows the structure of the component concentration measuring apparatus which concerns on the 2nd Embodiment of this invention. 本発明の第2の実施の形態に係る成分濃度測定装置の情報処理装置の構成を示すブロック図である。It is a block diagram which shows the structure of the information processing apparatus of the component concentration measuring apparatus which concerns on the 2nd Embodiment of this invention. 本発明の第2の実施の形態に係る成分濃度測定装置の動作を説明するフローチャートである。It is a flowchart explaining operation | movement of the component concentration measuring apparatus which concerns on the 2nd Embodiment of this invention. 本発明の第2の実施の形態における光強度変化の例を示す図である。It is a figure which shows the example of the light intensity change in the 2nd Embodiment of this invention. 本発明の第2の実施の形態における音響モードと光強度の変化量の比との関係の1例を示す図である。It is a figure which shows one example of the relationship between the acoustic mode in the 2nd Embodiment of this invention, and ratio of the variation | change_quantity of light intensity. 音波の反射率と音響インピーダンスとの関係を示す図である。It is a figure which shows the relationship between the reflectance of a sound wave, and acoustic impedance. 従来の成分濃度測定装置の構成を示すブロック図である。It is a block diagram which shows the structure of the conventional component density | concentration measuring apparatus. 従来の光パワーバランスシフト法を説明するフローチャートである。It is a flowchart explaining the conventional optical power balance shift method. 円筒空洞内での音波の共振モードを説明する図である。It is a figure explaining the resonance mode of the sound wave in a cylindrical cavity. 変調周波数と音圧との関係の1例を示す図である。It is a figure which shows one example of the relationship between a modulation frequency and a sound pressure.

[発明の原理]
本発明では、コリメートされた光を用いることで音波の干渉の影響を低減する。被測定物の形状が自由に変更できる場合は、音波の波長の半分に対して径の方が細い円筒形状の被測定物を用いることで、コリメート光と同等の効果を得ることができる。コリメートされた光は、強度分布として以下の式(1)のように表すことができる。
[Principle of the Invention]
In the present invention, the influence of sound wave interference is reduced by using collimated light. When the shape of the object to be measured can be freely changed, an effect equivalent to that of collimated light can be obtained by using the object to be measured having a cylindrical shape whose diameter is smaller than half the wavelength of the sound wave. The collimated light can be expressed as the following expression (1) as an intensity distribution.

また、光音響効果で生じる音波の大きさSは、以下の式(2)のように表すことができる。   Moreover, the magnitude | size S of the sound wave produced by a photoacoustic effect can be represented like the following formula | equation (2).

式(1)、式(2)において、Cは定数、αは特定の波長の光に対する被測定物の光吸収係数、qは音響モード、F(α,q)は光吸収係数αと音響モードqの関数、ωは変調周波数、wは光ビーム径、zは被測定物表面からの深さ、rは光軸から放射方向の距離、Iは光強度である。
光吸収係数αは式(3)のように記述することができる。
In equations (1) and (2), C is a constant, α is the light absorption coefficient of the object to be measured for light of a specific wavelength, q is the acoustic mode, and F (α, q) is the light absorption coefficient α and the acoustic mode. A function of q, ω is a modulation frequency, w is a light beam diameter, z is a depth from the surface of the object to be measured, r is a distance in the radial direction from the optical axis, and I is a light intensity.
The light absorption coefficient α can be described as in equation (3).

式(3)において、ΔCは被測定物内の特定成分の濃度、ΔTは被測定物の温度、βは特定成分のモル吸光度、γは被測定物の温度吸光度である。なお、α0は初期状態の被測定物の光吸収係数を表している。 In Expression (3), ΔC is the concentration of the specific component in the object to be measured, ΔT is the temperature of the object to be measured, β is the molar absorbance of the specific component, and γ is the temperature absorbance of the object to be measured. Α 0 represents the light absorption coefficient of the object to be measured in the initial state.

たとえば、円筒形の共振器内においては、音響波のモードは以下のような式(4)、式(5)で与えられる。音響波のモードには、長手方向、円周方向、半径方向のモードがある。   For example, in a cylindrical resonator, the acoustic wave mode is given by the following equations (4) and (5). Acoustic wave modes include longitudinal, circumferential, and radial modes.

式(4)、式(5)において、cは被測定物内の音速、qは光軸方向の音響モード、Rは円筒の半径、Lは円筒の長さ、αm,nは被測定物壁面で径方向の音圧変化が0という境界条件を満たすn次の根である。 In equations (4) and (5), c is the speed of sound in the object to be measured, q is the acoustic mode in the optical axis direction, R is the radius of the cylinder, L is the length of the cylinder, and α m, n is the object to be measured. This is the nth-order root that satisfies the boundary condition that the change in sound pressure in the radial direction on the wall surface is zero.

例えば、コリメートされた光を被測定物に照射する場合、qは円筒軸方向の定在波モードに相当する(図18)。また、関数Fは、この場合、近似的に式(6)のように記述することができるため、モードごとに離散的な値をとる。   For example, when the object to be measured is irradiated with collimated light, q corresponds to the standing wave mode in the cylindrical axis direction (FIG. 18). In this case, the function F can be described approximately as shown in the equation (6), and thus takes a discrete value for each mode.

ここで、光吸収係数αi、関数Fiのiは光の波長に対応する添え字である。平行平板間の音波の共振においても、円筒内の音波の共振と同様に扱うことができる。ここで、波長が異なり、180°位相がずれた2つの光を被測定物に照射する場合に発生する音波の大きさSは、それぞれの光により生ずる音波の重ね合わせとして次のように表すことができる。 Here, the light absorption coefficient α i and i of the function F i are subscripts corresponding to the wavelength of light. The resonance of the sound wave between the parallel plates can be handled in the same manner as the resonance of the sound wave in the cylinder. Here, the magnitude S of the sound wave generated when the object to be measured is irradiated with two lights having different wavelengths and having a phase difference of 180 ° is expressed as the superposition of the sound waves generated by the respective lights as follows. Can do.

1,S2は2つの音波の大きさである。OPBS法では、音波の信号強度が最少、すなわち理想的にゼロになる状態を考えることと同等であることから、S1=S2となる。
被測定物内の特定の成分がΔCだけ濃度変化したとすると、以下の式(8)が成立する。
S 1 and S 2 are the magnitudes of two sound waves. In the OPBS method, S 1 = S 2 because it is equivalent to considering a state where the signal intensity of the sound wave is minimum, that is, ideally zero.
Assuming that a specific component in the object to be measured changes in concentration by ΔC, the following equation (8) is established.

1は第1の光の強度、I2は第2の光の強度、ΔI2は被測定物内の特定成分がΔCだけ濃度変化したときの第2の光の強度の変化量である。式(8)から、当該成分濃度変化と光強度変化の関係は式(9)、式(10)のように近似できる。 I 1 is the intensity of the first light, I 2 is the intensity of the second light, and ΔI 2 is the amount of change in the intensity of the second light when the concentration of the specific component in the object to be measured changes by ΔC. From the equation (8), the relationship between the component concentration change and the light intensity change can be approximated as in the equations (9) and (10).

式(9)、式(10)より音響モードqおよび背景光吸収が分かれば、濃度変化の感度が求められる。被測定物をほぼ水と仮定し、グルコース濃度が変化したときの応答特性の例を図1に示す。図1によれば、グルコース濃度に比例して光強度が変化していくこと、音響モードqにより感度が変化することが示される。また、音響モードqが変化すると、感度が大きく変化してしまうことが分かる。   If the acoustic mode q and the background light absorption are known from the equations (9) and (10), the sensitivity of density change can be obtained. FIG. 1 shows an example of response characteristics when the measurement object is almost water and the glucose concentration changes. FIG. 1 shows that the light intensity changes in proportion to the glucose concentration and that the sensitivity changes depending on the acoustic mode q. It can also be seen that the sensitivity changes greatly when the acoustic mode q changes.

そこで、音響モードの補正を行う必要がある。この補正は、測定に用いる波長の光とは別の第3の波長の光あるいは温度に対する吸収スペクトルのデータセットを用意し、実測データに対する数値的補正を行うことで実現することができる。式(9)、式(10)および既知の吸収係数差を用いて音波の共振モードの干渉誤差を補正する本発明の方法を2つ説明する。   Therefore, it is necessary to correct the acoustic mode. This correction can be realized by preparing a data set of an absorption spectrum with respect to light having a third wavelength different from the light having a wavelength used for measurement or temperature, and performing numerical correction on actually measured data. Two methods of the present invention for correcting the interference error of the resonance mode of the sound wave using the equations (9) and (10) and the known absorption coefficient difference will be described.

[第1の補正方法]
第1の波長の光を被測定物に照射し、光の変調周波数を任意の範囲で掃引することで図2に示すような変調周波数と音圧との関係を取得する。この図2のピークから被測定物内での音波の共振周波数を選択することができる。
[First correction method]
The relationship between the modulation frequency and the sound pressure as shown in FIG. 2 is acquired by irradiating the object to be measured with light of the first wavelength and sweeping the modulation frequency of the light in an arbitrary range. The resonance frequency of the sound wave in the object to be measured can be selected from the peak in FIG.

次に前述のとおり、光波長が異なり位相差が180°の2つの光ビームを被測定物に対して、同一の光出力口から照射する。このときの光変調周波数は上記で選択した変調周波数である。第1の波長の光の強度を一定に固定し、第2の波長の光の強度を増減させながら、生じる音波強度が極小となる、第2の波長の光の強度I0を測定する。 Next, as described above, two light beams having different light wavelengths and a phase difference of 180 ° are irradiated to the object to be measured from the same light output port. The light modulation frequency at this time is the modulation frequency selected above. The intensity I 0 of the second wavelength light is measured while the intensity of the generated sound wave is minimized while the intensity of the first wavelength light is fixed and the intensity of the second wavelength light is increased or decreased.

次に、第2の波長の光を被測定物に対する光吸収係数が既知の第3の波長の光に切り替えて、上記と同様に第3の波長の光の強度を増減させながら、生じる音波強度が極小となる、第3の波長の光の強度I1を測定する。求めた光強度I0,I1より光強度I0に対する光強度の変化量(I1−I0)の比ΔI/Iは式(11)のようになる。 Next, by switching the light of the second wavelength to the light of the third wavelength whose light absorption coefficient for the object to be measured is known, the sound wave intensity generated while increasing or decreasing the intensity of the light of the third wavelength in the same manner as described above. Measure the intensity I 1 of the light of the third wavelength that minimizes. The ratio ΔI / I of the change amount (I 1 −I 0 ) of the light intensity with respect to the light intensity I 0 from the obtained light intensity I 0 , I 1 is expressed by the equation (11).

ここで、十分コリメートされた光を用いた場合、上式の関係は各音響モードqに対して次式で与えられる。   Here, when sufficiently collimated light is used, the relationship of the above equation is given by the following equation for each acoustic mode q.

関数Fi,光吸収係数αiのiは光の波長に対応する添え字である。F1は第1の波長の光に関する光吸収係数αと音響モードqの関数、F2は第2の波長の光に関する光吸収係数αと音響モードqの関数、F1’は第3の波長を用いるときの第1の波長の光に関する光吸収係数αと音響モードqの関数、F2’は第3の波長の光に関する光吸収係数αと音響モードqの関数、α2は第2の波長の光に対する被測定物の光吸収係数である。 The function F i and i of the light absorption coefficient α i are subscripts corresponding to the wavelength of light. F 1 is a function of light absorption coefficient α and acoustic mode q for light of the first wavelength, F 2 is a function of light absorption coefficient α and acoustic mode q for light of the second wavelength, and F 1 ′ is a third wavelength. Is a function of the light absorption coefficient α and the acoustic mode q for the first wavelength light, F 2 ′ is a function of the light absorption coefficient α and the acoustic mode q for the third wavelength light, and α 2 is the second function. It is a light absorption coefficient of the object to be measured with respect to light of a wavelength.

よって、式(11)と式(12)の応答を比較し、最も濃度測定感度が近しい応答を探し音響モードqを求めることができる。音響モードqを用いて式(9)より現在の変調周波数における濃度測定感度Gqが推測できる。この濃度測定感度Gqを音響モードq=0のときの濃度測定感度G0と同じになるように定数係数Aを決定すると、補償された濃度との関係は次式で表される。   Therefore, the acoustic mode q can be obtained by comparing the responses of the equations (11) and (12) and searching for the response with the closest density measurement sensitivity. The density measurement sensitivity Gq at the current modulation frequency can be estimated from the equation (9) using the acoustic mode q. When the constant coefficient A is determined so that the density measurement sensitivity Gq is the same as the density measurement sensitivity G0 when the acoustic mode q = 0, the relationship with the compensated density is expressed by the following equation.

[第2の補正方法]
第1の補正方法と同様に、第1の波長の光を被測定物に照射し、光の変調周波数を任意の範囲で掃引することで図2に示すような変調周波数と音圧との関係を取得する。この図2のピークから被測定物内での音波の共振周波数を選択することができる。
[Second correction method]
Similar to the first correction method, the relationship between the modulation frequency and the sound pressure as shown in FIG. 2 is obtained by irradiating the object to be measured with light of the first wavelength and sweeping the modulation frequency of the light in an arbitrary range. To get. The resonance frequency of the sound wave in the object to be measured can be selected from the peak in FIG.

次に前述のとおり、光波長が異なり位相差が180°の2つの光ビームを被測定物に対して、同一光出力口から照射する。このときの変調周波数は上記で選択した変調周波数である。第1の波長の光の強度を一定に固定し、第2の波長の光の強度を増減させながら、生じる音波強度が極小となる、第2の波長の光の強度I0を測定する。 Next, as described above, two light beams having different light wavelengths and a phase difference of 180 ° are irradiated to the object to be measured from the same light output port. The modulation frequency at this time is the modulation frequency selected above. The intensity I 0 of the second wavelength light is measured while the intensity of the generated sound wave is minimized while the intensity of the first wavelength light is fixed and the intensity of the second wavelength light is increased or decreased.

次に、温度制御器を用いて被測定物の温度をΔTだけ変化させ、上記と同様に第2の波長の光の強度を増減させながら、生じる音波強度が極小となる、第2の波長の光の強度I1を測定する。求めた光強度I0,I1より光強度I0に対する光強度の変化量(I1−I0)の比ΔI/Iは式(15)のようになる。 Next, the temperature of the object to be measured is changed by ΔT using a temperature controller, and the intensity of the generated sound wave is minimized while increasing or decreasing the intensity of the second wavelength light in the same manner as described above. The light intensity I 1 is measured. The ratio ΔI / I of the change amount (I 1 −I 0 ) of the light intensity with respect to the light intensity I 0 from the obtained light intensities I 0 and I 1 is expressed by the following equation (15).

被測定物の温度変化に対する光吸収係数変化は、スペクトルデータセットより光吸収係数の温度依存性γで与えた場合、次式で温度に対する応答特性が与えられる。   When the change in the light absorption coefficient with respect to the temperature change of the object to be measured is given by the temperature dependence γ of the light absorption coefficient from the spectral data set, the response characteristic with respect to temperature is given by

γ1は第1の波長の光に対する被測定物の光吸収係数α1の温度依存性、γ2は第2の波長の光に対する被測定物の光吸収係数α2の温度依存性である。 γ 1 is the temperature dependence of the light absorption coefficient α 1 of the object to be measured with respect to the first wavelength light, and γ 2 is the temperature dependence of the light absorption coefficient α 2 of the object to be measured with respect to the light of the second wavelength.

よって、式(15)と式(16)の応答を比較し、最も感度が近しい応答を探し音響モードqを求めることができる。音響モードqを用いて式(9)より現在の変調周波数における濃度測定感度Gqが推測できる。この濃度測定感度Gqを音響モードq=0のときの濃度測定感度G0と同じになるように定数係数Aを決定すると、補償された濃度との関係は次式で表される。   Therefore, the acoustic mode q can be obtained by comparing the responses of the equations (15) and (16) and searching for the response with the closest sensitivity. The density measurement sensitivity Gq at the current modulation frequency can be estimated from the equation (9) using the acoustic mode q. When the constant coefficient A is determined so that the density measurement sensitivity Gq is the same as the density measurement sensitivity G0 when the acoustic mode q = 0, the relationship with the compensated density is expressed by the following equation.

光音響法の音波共振モードの補償方法についてグルコース水溶液の場合を例にとって下記に示す。   A method for compensating the acoustic resonance mode of the photoacoustic method will be described below using a glucose aqueous solution as an example.

[第1の実施の形態]
以下、本発明の実施の形態について図面を参照して説明する。図3は本発明の第1の実施の形態に係る成分濃度測定装置の構成を示すブロック図である。本実施の形態の成分濃度測定装置は、レーザ光を放射するレーザダイオード1−1,1−2,1−3と、レーザダイオード1−1,1−2,1−3を駆動するレーザドライバ2と、レーザダイオード1−1,1−2,1−3から放射されたレーザ光を導く光ファイバ3−1,3−2,3−3と、レーザダイオード1−1,1−2,1−3から放射されたレーザ光を合波する光カプラ4と、光カプラ4によって合波されたレーザ光を導く光ファイバ5と、光音響効果によって被測定物11から発生する光音響信号を検出し、音圧に比例した電気信号に変換する光音響信号検出手段となる音響センサ6と、音響センサ6から出力された電気信号を増幅する増幅器7と、参照信号を発生するファンクションジェネレータ8と、増幅器7の出力信号とファンクションジェネレータ8から出力された参照信号とを入力として、増幅器7の出力信号から所望の周波数の測定信号を検出するロックインアンプ9と、レーザドライバ2とファンクションジェネレータ8とロックインアンプ9とを制御すると共に、ロックインアンプ9が検出した測定信号を処理して特定の成分濃度を導出するコンピュータからなる情報処理装置10とから構成される。
[First Embodiment]
Hereinafter, embodiments of the present invention will be described with reference to the drawings. FIG. 3 is a block diagram showing the configuration of the component concentration measuring apparatus according to the first embodiment of the present invention. The component concentration measuring apparatus according to the present embodiment includes laser diodes 1-1, 1-2, 1-3 that emit laser light, and a laser driver 2 that drives the laser diodes 1-1, 1-2, 1-3. And optical fibers 3-1, 3-2 and 3-3 for guiding the laser light emitted from the laser diodes 1-1, 1-2 and 1-3, and the laser diodes 1-1, 1-2 and 1-1. 3 detects the photoacoustic signal generated from the object to be measured 11 by the photoacoustic effect, the optical coupler 4 that combines the laser light emitted from the optical coupler 4, the optical fiber 5 that guides the laser light combined by the optical coupler 4. , An acoustic sensor 6 serving as a photoacoustic signal detecting means for converting the electrical signal proportional to the sound pressure, an amplifier 7 for amplifying the electrical signal output from the acoustic sensor 6, a function generator 8 for generating a reference signal, and an amplifier 7 output signals and A lock-in amplifier 9 that receives a reference signal output from the function generator 8 and detects a measurement signal of a desired frequency from an output signal of the amplifier 7, a laser driver 2, a function generator 8, and a lock-in amplifier 9. The information processing apparatus 10 includes a computer that controls and processes a measurement signal detected by the lock-in amplifier 9 to derive a specific component concentration.

レーザダイオード1−1,1−2,1−3とレーザドライバ2とファンクションジェネレータ8とは、光照射手段を構成している。情報処理装置10とファンクションジェネレータ8とは、光強度制御手段を構成している。   The laser diodes 1-1, 1-2, 1-3, the laser driver 2, and the function generator 8 constitute light irradiation means. The information processing apparatus 10 and the function generator 8 constitute light intensity control means.

レーザダイオード1−1,1−2,1−3の例としては、例えば分布帰還型半導体レーザ(DFB−LD)等がある。音響センサ6の例としては、圧電センサを用いるマイクロホンがある。
なお、被測定物11は例えば円筒形の光音響セル(不図示)内に導入され、光音響セルに設けられたガラス製の光学窓を通してレーザ光が照射されるようになっている。
Examples of the laser diodes 1-1, 1-2, and 1-3 include, for example, a distributed feedback semiconductor laser (DFB-LD). An example of the acoustic sensor 6 is a microphone using a piezoelectric sensor.
The object to be measured 11 is introduced into, for example, a cylindrical photoacoustic cell (not shown), and laser light is irradiated through a glass optical window provided in the photoacoustic cell.

図4は情報処理装置10の構成を示すブロック図である。情報処理装置10は、ファンクションジェネレータ8を制御するジェネレータ制御部100と、ファンクションジェネレータ8を介して光の強度を制御する光強度制御部101と、レーザダイオード1−1〜1−3から放射された光の強度を測定する光強度測定手段となる光強度測定部102と、強度を変化させる光の波長が異なる2つの時点または被測定物11の温度が異なる2つの時点について、それぞれ測定信号の強度が最低となったときの光強度の測定結果から得られる光強度の変化量を算出する光強度変化量算出部103と、強度を変化させる光の波長が異なる2つの時点または被測定物11の温度が異なる2つの時点について、それぞれ測定信号の強度が最低となったときの光強度の測定結果から得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、光音響信号の共振モードである音響モードを推定する音響モード推定部104と、音響モード推定部104が推定した音響モードと、被測定物11に含まれる特定成分(例えばグルコース)に対する濃度測定感度と音響モードとの既知の関係に基づいて、特定成分に対する濃度測定感度の補正係数を決定する補正係数決定部105と、被測定物11内の特定成分の濃度を導出する濃度導出部106と、情報記憶のための記憶部107とを有する。   FIG. 4 is a block diagram illustrating a configuration of the information processing apparatus 10. The information processing apparatus 10 is radiated from the generator control unit 100 that controls the function generator 8, the light intensity control unit 101 that controls the intensity of light via the function generator 8, and the laser diodes 1-1 to 1-3. The intensity of the measurement signal at each of the light intensity measuring unit 102 serving as a light intensity measuring means for measuring the light intensity and two time points at which the wavelength of the light whose intensity is changed are different or two time points at which the temperature of the object to be measured 11 is different. Of the light intensity change amount calculation unit 103 for calculating the change amount of the light intensity obtained from the measurement result of the light intensity when the light intensity becomes the minimum, and two time points at which the wavelengths of the light whose intensity changes are different or the measurement object 11 For two time points with different temperatures, the amount of change in the light intensity obtained from the measurement result of the light intensity when the intensity of the measurement signal is the lowest, The acoustic mode estimation unit 104 that estimates the acoustic mode that is the resonance mode of the photoacoustic signal based on the amount of change in the light intensity obtained from the data set, the acoustic mode estimated by the acoustic mode estimation unit 104, and the device under test A correction coefficient determining unit 105 that determines a correction coefficient of the concentration measurement sensitivity for the specific component based on a known relationship between the concentration measurement sensitivity for the specific component (for example, glucose) and the acoustic mode included in the measurement object 11, A concentration deriving unit 106 for deriving the concentration of the specific component, and a storage unit 107 for storing information.

次に、本実施の形態の成分濃度測定装置の動作を図5のフローチャートを参照して説明する。初めに、レーザダイオード1−1(第1の光源)のみを動作させる。レーザドライバ2から駆動電流が供給されると、レーザダイオード1−1はレーザ光(連続光)を放射する。このとき、レーザドライバ2から矩形波の駆動電流が供給されることにより、レーザダイオード1−1は強度変調光を放射する。レーザダイオード1−1から放射される光の波長(第1の波長)は例えば1382nmである。レーザダイオード1−1から放射された強度変調光は、光ファイバ3−1によって導かれ、光カプラ4および光ファイバ5を介して被測定物11に照射される(図5ステップS200)。   Next, the operation of the component concentration measuring apparatus according to the present embodiment will be described with reference to the flowchart of FIG. First, only the laser diode 1-1 (first light source) is operated. When a drive current is supplied from the laser driver 2, the laser diode 1-1 emits laser light (continuous light). At this time, when a rectangular-wave drive current is supplied from the laser driver 2, the laser diode 1-1 emits intensity-modulated light. The wavelength (first wavelength) of the light emitted from the laser diode 1-1 is, for example, 1382 nm. The intensity-modulated light radiated from the laser diode 1-1 is guided by the optical fiber 3-1, and irradiated to the device under test 11 through the optical coupler 4 and the optical fiber 5 (step S200 in FIG. 5).

音響センサ6は、被測定物11から発生する光音響信号を検出し、増幅器7は、音響センサ8から出力された電気信号を増幅する。ロックインアンプ9は、増幅器7の出力に含まれる信号のうち、ファンクションジェネレータ8から出力される参照信号によって決まる周波数の測定信号を検出する。   The acoustic sensor 6 detects a photoacoustic signal generated from the device under test 11, and the amplifier 7 amplifies the electrical signal output from the acoustic sensor 8. The lock-in amplifier 9 detects a measurement signal having a frequency determined by the reference signal output from the function generator 8 among the signals included in the output of the amplifier 7.

情報処理装置10のジェネレータ制御部100は、ファンクションジェネレータ8が発生する参照信号の周波数を変化させることにより、レーザドライバ2からレーザダイオード1−1に供給される駆動電流の周波数を変化させ、光変調周波数を漸次変化させると共に、ロックインアンプ11が検出する測定信号の周波数(光変調周波数と同一の周波数)を漸次変化させる光変調周波数掃引を行う(図5ステップS201)。   The generator control unit 100 of the information processing apparatus 10 changes the frequency of the reference current generated by the function generator 8 to change the frequency of the drive current supplied from the laser driver 2 to the laser diode 1-1, thereby performing optical modulation. The optical modulation frequency sweep is performed to gradually change the frequency and gradually change the frequency of the measurement signal detected by the lock-in amplifier 11 (the same frequency as the optical modulation frequency) (step S201 in FIG. 5).

こうして、変調周波数と測定信号の振幅(音圧)との関係を取得し(図5ステップS202)、測定信号の振幅が最大の変調周波数であるピーク変調周波数を探索する(図5ステップS203)。本実施の形態において得られた変調周波数と音波強度との関係を図6に示す。   In this way, the relationship between the modulation frequency and the amplitude (sound pressure) of the measurement signal is acquired (step S202 in FIG. 5), and the peak modulation frequency that is the modulation frequency with the maximum amplitude of the measurement signal is searched (step S203 in FIG. 5). FIG. 6 shows the relationship between the modulation frequency and sound wave intensity obtained in the present embodiment.

次に、レーザダイオード1−1を再び動作させて連続光を出力させ(図5ステップS204)、続いてレーザダイオード1−2(第2の光源)を動作させて連続光を出力させる(図5ステップS205)。レーザダイオード1−1から放射される光の波長は例えば1382nm、レーザダイオード1−2から放射される光の波長(第2の波長)は例えば1610nmであり、2つのレーザダイオード1−1,1−2から放射される光の波長は異なる。   Next, the laser diode 1-1 is operated again to output continuous light (step S204 in FIG. 5), and then the laser diode 1-2 (second light source) is operated to output continuous light (FIG. 5). Step S205). The wavelength of the light emitted from the laser diode 1-1 is, for example, 1382 nm, and the wavelength (second wavelength) of the light emitted from the laser diode 1-2 is, for example, 1610 nm. The two laser diodes 1-1, 1- The wavelengths of light emitted from the two are different.

情報処理装置10のジェネレータ制御部100は、ファンクションジェネレータ8が発生する2つの参照信号の周波数を変化させることにより、レーザドライバ2からレーザダイオード1−1,1−2に供給される駆動電流の周波数を変化させ、光変調周波数を上記のピーク変調周波数に設定すると共に、ロックインアンプ9が検出する測定信号の周波数をピーク変調周波数に設定する。このとき、レーザドライバ2から、同一周波数で逆位相の矩形波の駆動電流をレーザダイオード1−1,1−2に供給させることにより、レーザダイオード1−1,1−2から放射される光を同一周波数(ピーク変調周波数)で逆位相の信号によりそれぞれ強度変調する(図5ステップS206)。したがって、本実施の形態では、互いに異なる波長の2波のレーザ光を同一周波数で逆位相の信号によりそれぞれ強度変調して2つの強度変調光を生成することになる。   The generator control unit 100 of the information processing apparatus 10 changes the frequency of the two reference signals generated by the function generator 8 to change the frequency of the drive current supplied from the laser driver 2 to the laser diodes 1-1 and 1-2. And the optical modulation frequency is set to the peak modulation frequency, and the frequency of the measurement signal detected by the lock-in amplifier 9 is set to the peak modulation frequency. At this time, light emitted from the laser diodes 1-1 and 1-2 is supplied from the laser driver 2 to the laser diodes 1-1 and 1-2 by supplying a rectangular-wave drive current of the same frequency and opposite phase. Intensity modulation is performed with signals having the same frequency (peak modulation frequency) and opposite phases (step S206 in FIG. 5). Therefore, in this embodiment, two intensity-modulated lights are generated by intensity-modulating two laser beams having different wavelengths with signals having the same frequency and opposite phases.

レーザダイオード1−1,1−2から放射された強度変調光は、それぞれ光ファイバ3−1,3−2によって導かれ、光カプラ4によって合波され、さらに光ファイバ5によって導かれ、被測定物11に照射される。   The intensity-modulated lights emitted from the laser diodes 1-1 and 1-2 are respectively guided by the optical fibers 3-1 and 3-2, combined by the optical coupler 4, and further guided by the optical fiber 5 to be measured. The object 11 is irradiated.

情報処理装置10の光強度制御部101は、ファンクションジェネレータ8を制御して、レーザダイオード1−1を駆動するための参照信号の電圧を所定値に設定することにより、レーザドライバ2からレーザダイオード1−1に供給される駆動電流の大きさを一定にし、レーザダイオード1−1から放射される光の強度を一定にする(図5ステップS207)。   The light intensity control unit 101 of the information processing apparatus 10 controls the function generator 8 to set the voltage of the reference signal for driving the laser diode 1-1 to a predetermined value. -1 is made constant, and the intensity of light emitted from the laser diode 1-1 is made constant (step S207 in FIG. 5).

また、情報処理装置10の光強度制御部101は、ファンクションジェネレータ8を制御して、レーザダイオード1−2を駆動するための参照信号の電圧を変化させることにより、レーザドライバ2からレーザダイオード1−2に供給される駆動電流の大きさを変化させ、レーザダイオード1−2から放射される光の強度を漸次変化させる光強度掃引を行い、ロックインアンプ9から出力される測定信号の振幅が最小になるようにレーザダイオード1−2から放射される光の強度を調節する(図5ステップS208)。   Further, the light intensity control unit 101 of the information processing apparatus 10 controls the function generator 8 to change the voltage of the reference signal for driving the laser diode 1-2, so that the laser driver 2 changes the laser diode 1-. The intensity of the measurement signal output from the lock-in amplifier 9 is minimized by performing a light intensity sweep that changes the intensity of the drive current supplied to 2 and gradually changes the intensity of light emitted from the laser diode 1-2. The intensity of light emitted from the laser diode 1-2 is adjusted so as to become (step S208 in FIG. 5).

情報処理装置10の光強度測定部102は、レーザダイオード1−2から放射される光の強度を測定する。光強度測定部102は、ファンクションジェネレータ8から出力されている、レーザダイオード1−2の駆動用の参照信号の電圧を取得する。情報処理装置10の記憶部107には、参照信号の電圧とレーザダイオード1−2から放射される光の強度との関係を示すキャリブレーションデータが予め記憶されている。キャリブレーションデータは、電圧と強度の実測により予め求めておくことができる。光強度測定部102は、このようなキャリブレーションデータを参照し、取得した参照信号の電圧を、レーザダイオード1−2から放射された光の強度に換算する。光強度測定部102は、測定信号の振幅が最小となるときの、レーザダイオード1−2から放射される光の強度I0を記憶部107に格納する(図5ステップS209)。 The light intensity measurement unit 102 of the information processing apparatus 10 measures the intensity of light emitted from the laser diode 1-2. The light intensity measurement unit 102 acquires the voltage of the reference signal for driving the laser diode 1-2 output from the function generator 8. The storage unit 107 of the information processing apparatus 10 stores in advance calibration data indicating the relationship between the voltage of the reference signal and the intensity of light emitted from the laser diode 1-2. Calibration data can be obtained in advance by actual measurement of voltage and intensity. The light intensity measurement unit 102 refers to such calibration data, and converts the voltage of the acquired reference signal into the intensity of light emitted from the laser diode 1-2. The light intensity measurement unit 102 stores the intensity I 0 of the light emitted from the laser diode 1-2 when the amplitude of the measurement signal is minimum in the storage unit 107 (step S209 in FIG. 5).

次に、レーザダイオード1−2の動作を停止させて、レーザダイオード1−3(第3の光源)を動作させて連続光を出力させる(図5ステップS210)。レーザダイオード1−3から放射される光の波長(第3の波長)は例えば1650nmであり、2つのレーザダイオード1−1,1−3から放射される光の波長は異なる。   Next, the operation of the laser diode 1-2 is stopped, and the laser diode 1-3 (third light source) is operated to output continuous light (step S210 in FIG. 5). The wavelength of light emitted from the laser diode 1-3 (third wavelength) is, for example, 1650 nm, and the wavelengths of light emitted from the two laser diodes 1-1 and 1-3 are different.

情報処理装置10のジェネレータ制御部100は、ファンクションジェネレータ8が発生する2つの参照信号の周波数を変化させることにより、レーザドライバ2からレーザダイオード1−1,1−3に供給される駆動電流の周波数を変化させ、光変調周波数を上記のピーク変調周波数に設定すると共に、ロックインアンプ9が検出する測定信号の周波数をピーク変調周波数に設定する。このとき、レーザドライバ2から、同一周波数で逆位相の矩形波の駆動電流をレーザダイオード1−1,1−3に供給させることにより、レーザダイオード1−1,1−3から放射される光を同一周波数(ピーク変調周波数)で逆位相の信号によりそれぞれ強度変調する(図5ステップS211)。   The generator control unit 100 of the information processing apparatus 10 changes the frequency of the two reference signals generated by the function generator 8 to change the frequency of the drive current supplied from the laser driver 2 to the laser diodes 1-1 and 1-3. And the optical modulation frequency is set to the peak modulation frequency, and the frequency of the measurement signal detected by the lock-in amplifier 9 is set to the peak modulation frequency. At this time, light emitted from the laser diodes 1-1 and 1-3 is supplied from the laser driver 2 to the laser diodes 1-1 and 1-3 by supplying a rectangular-wave drive current having the same frequency and opposite phase. Intensity modulation is performed with signals having the same frequency (peak modulation frequency) and opposite phases (step S211 in FIG. 5).

情報処理装置10の光強度制御部101は、ファンクションジェネレータ8を制御して、レーザダイオード1−3を駆動するための参照信号の電圧を変化させることにより、レーザドライバ2からレーザダイオード1−3に供給される駆動電流の大きさを変化させ、レーザダイオード1−3から放射される光の強度を漸次変化させる光強度掃引を行い、ロックインアンプ9から出力される測定信号の振幅が最小になるようにレーザダイオード1−3から放射される光の強度を調節する(図5ステップS212)。上記のとおり、レーザダイオード1−1から放射される光の強度は一定である。   The light intensity control unit 101 of the information processing apparatus 10 controls the function generator 8 to change the voltage of the reference signal for driving the laser diode 1-3 to change the laser driver 2 to the laser diode 1-3. The intensity of the measurement signal output from the lock-in amplifier 9 is minimized by changing the magnitude of the supplied drive current and performing a light intensity sweep that gradually changes the intensity of the light emitted from the laser diode 1-3. In this way, the intensity of light emitted from the laser diode 1-3 is adjusted (step S212 in FIG. 5). As described above, the intensity of light emitted from the laser diode 1-1 is constant.

情報処理装置10の光強度測定部102は、レーザダイオード1−3から放射される光の強度を測定する。光強度測定部102は、ファンクションジェネレータ8から出力されている、レーザダイオード1−3の駆動用の参照信号の電圧を取得する。上記と同様に、情報処理装置10の記憶部107には、参照信号の電圧とレーザダイオード1−3から放射される光の強度との関係を示すキャリブレーションデータが予め記憶されている。光強度測定部102は、このようなキャリブレーションデータを参照し、取得した参照信号の電圧を、レーザダイオード1−3から放射された光の強度に換算する。光強度測定部102は、測定信号の振幅が最小となるときの、レーザダイオード1−3から放射される光の強度I1を記憶部107に格納する(図5ステップS213)。 The light intensity measurement unit 102 of the information processing apparatus 10 measures the intensity of light emitted from the laser diode 1-3. The light intensity measurement unit 102 acquires the voltage of the reference signal for driving the laser diode 1-3 output from the function generator 8. Similarly to the above, the storage unit 107 of the information processing apparatus 10 stores in advance calibration data indicating the relationship between the voltage of the reference signal and the intensity of light emitted from the laser diode 1-3. The light intensity measurement unit 102 refers to such calibration data, and converts the voltage of the acquired reference signal into the intensity of light emitted from the laser diode 1-3. The light intensity measurement unit 102 stores the intensity I 1 of the light emitted from the laser diode 1-3 when the amplitude of the measurement signal is minimum in the storage unit 107 (step S213 in FIG. 5).

情報処理装置10の光強度変化量算出部103は、ステップS209,S213で求めた光強度I0,I1より光強度I0に対する光強度の変化量ΔI=(I1−I0)の比ΔI/Iを、式(11)により計算する(図5ステップS214)。 The light intensity change amount calculation unit 103 of the information processing apparatus 10 compares the light intensity change amount ΔI = (I 1 −I 0 ) with respect to the light intensity I 0 based on the light intensities I 0 and I 1 obtained in steps S209 and S213. ΔI / I is calculated by equation (11) (step S214 in FIG. 5).

レーザダイオード1−3から放射される第3の波長の光に対する被測定物11の光吸収係数がレーザダイオード1−2から放射される第2の波長の光に対する被測定物11の光吸収係数よりも0.02/mm大きい場合の光強度の変化量の比ΔI/Iを図7に示す。図7によると、ΔI/I=4.7%という結果が得られる。   The light absorption coefficient of the device under test 11 for the third wavelength light emitted from the laser diode 1-3 is based on the light absorption coefficient of the device under test 11 for the second wavelength light emitted from the laser diode 1-2. FIG. 7 shows the ratio ΔI / I of the change amount of the light intensity when 0.02 / mm is larger. According to FIG. 7, the result ΔI / I = 4.7% is obtained.

次に、情報処理装置10の音響モード推定部104は、表1のような所定のデータセットを用いて、光強度の変化量の比ΔI2/I2を式(12)、式(13)により音響モードq毎に計算する(図5ステップS215)。データセットは、第1、第2、第3の波長と、これらの波長の光に対する特定成分(例えばグルコース)の光吸収係数との関係を表すものである。 Next, the acoustic mode estimation unit 104 of the information processing apparatus 10 uses the predetermined data set as shown in Table 1 to calculate the light intensity change amount ratio ΔI 2 / I 2 using the equations (12) and (13). Is calculated for each acoustic mode q (step S215 in FIG. 5). The data set represents the relationship between the first, second, and third wavelengths and the light absorption coefficient of a specific component (for example, glucose) with respect to light of these wavelengths.

関数F1は第1の波長の光に関する光吸収係数α1=0.30/mmと音響モードqから計算することができ、関数F2は第2の波長の光に関する光吸収係数α2=0.28/mmと音響モードqから計算することができる。Δα2は光吸収係数α1に対する光吸収係数α2の変化量である。被測定物11の長さLは既知の値である。 The function F 1 can be calculated from the light absorption coefficient α 1 = 0.30 / mm for the first wavelength light and the acoustic mode q, and the function F 2 is the light absorption coefficient α 2 = for the second wavelength light. It can be calculated from 0.28 / mm and the acoustic mode q. Δα 2 is a change amount of the light absorption coefficient α 2 with respect to the light absorption coefficient α 1 . The length L of the DUT 11 is a known value.

ステップS215の計算の結果から、図8に示すような音響モードqと光強度の変化量の比ΔI2/I2との関係が得られる。情報処理装置10の音響モード推定部104は、ステップS214で計算された光強度の変化量の比ΔI/Iと、ステップS215で計算した光強度の変化量の比ΔI2/I2とが等しくなる音響モードqの値を、今回のピーク変調周波数での音響モードの値として決定する(図5ステップS216)。本実施の形態の例では、q=2であることが推測できる。 The relationship between the acoustic mode q and the light intensity change ratio ΔI 2 / I 2 as shown in FIG. 8 is obtained from the calculation result of step S215. The acoustic mode estimation unit 104 of the information processing apparatus 10 has the light intensity change ratio ΔI / I calculated in step S214 equal to the light intensity change ratio ΔI 2 / I 2 calculated in step S215. The value of the acoustic mode q is determined as the value of the acoustic mode at the current peak modulation frequency (step S216 in FIG. 5). In the example of the present embodiment, it can be estimated that q = 2.

情報処理装置10の補正係数決定部105は、音響モード推定部104が推定した音響モードqの値に対応する、被測定物11内の特定成分(本実施の形態の例ではグルコース)の濃度の測定感度Gqを、濃度測定感度と音響モードqとの既知の関係(図9)から求め、この濃度測定感度Gqに対する音響モードq=0のときの濃度測定感度G0の比G0/Gqを補正係数Aとして決定する(図5ステップS217)。   The correction coefficient determination unit 105 of the information processing apparatus 10 has a concentration of a specific component (glucose in the example of the present embodiment) in the DUT 11 corresponding to the value of the acoustic mode q estimated by the acoustic mode estimation unit 104. The measurement sensitivity Gq is obtained from a known relationship between the density measurement sensitivity and the acoustic mode q (FIG. 9), and the ratio G0 / Gq of the density measurement sensitivity G0 when the acoustic mode q = 0 with respect to the density measurement sensitivity Gq is a correction coefficient. A is determined (step S217 in FIG. 5).

図9によると、音響モードq=2のときの濃度測定感度Gqは音響モードq=0のときの濃度測定感度G0より約5倍大きいので、補正係数決定部105は、補正係数A=1/5と求めることができる。   According to FIG. 9, the density measurement sensitivity Gq when the acoustic mode q = 2 is approximately five times larger than the density measurement sensitivity G0 when the acoustic mode q = 0, so that the correction coefficient determination unit 105 performs the correction coefficient A = 1 / 5 can be obtained.

情報処理装置10の濃度導出部106は、ステップS214で計算された光強度の変化量の比ΔI/I(音響モードq=2のときの光強度の変化量の比ΔI2/I2)と、補正係数決定部105が決定した補正係数Aとから、測定開始時から任意の時間経過後の時点における被測定物11内の特定成分の濃度の変化量ΔC[%]を式(14)により計算する(図5ステップS218)。 The concentration deriving unit 106 of the information processing apparatus 10 calculates the light intensity change ratio ΔI / I (the light intensity change ratio ΔI 2 / I 2 when the acoustic mode q = 2) calculated in step S214. Based on the correction coefficient A determined by the correction coefficient determination unit 105, the amount of change ΔC [%] of the concentration of the specific component in the DUT 11 at an arbitrary time after the start of measurement is calculated by the equation (14). Calculation is performed (step S218 in FIG. 5).

最後に、濃度導出部106は、ステップS218で計算した特定成分の濃度の変化量ΔC[%]と特定成分の既知の参照濃度とから、任意の時間経過後の時点における特定成分の濃度を計算する(図5ステップS219)。なお、参照濃度は、被測定物11に対して標準的な血糖測定法を実施することにより得ることができる。標準的な血糖測定法を実施するには、血糖測定器の本体に、グルコースセンサーを差し込み、針を専用の機械(または本体)にセットして、指などから採血し、グルコースセンサーに血を吸収させる。標準的な血糖測定法は、既知濃度のグルコース液を標準校正液として機械動作確認用に用いる。初期動作時に機械が正常に動いているかを確認したり、血糖値が異常値にあるか(正常に機械が動作しているか)を確認したりするときに用いる。このような標準的な血糖測定法を、測定開始時の例えばピーク変調周波数の選択時に実施しておけばよい。   Finally, the concentration deriving unit 106 calculates the concentration of the specific component at an arbitrary time point after the elapse of an arbitrary time from the change amount ΔC [%] of the concentration of the specific component calculated in step S218 and the known reference concentration of the specific component. (Step S219 in FIG. 5). The reference concentration can be obtained by performing a standard blood glucose measurement method on the object to be measured 11. To perform a standard blood glucose measurement method, insert a glucose sensor into the body of the blood glucose meter, set the needle on a dedicated machine (or body), collect blood from a finger, etc., and absorb the blood into the glucose sensor. Let A standard blood glucose measurement method uses a known concentration of glucose solution as a standard calibration solution for machine operation confirmation. It is used to check whether the machine is operating normally during initial operation, or to check whether the blood glucose level is abnormal (ie, whether the machine is operating normally). Such a standard blood glucose measurement method may be performed at the time of starting measurement, for example, when selecting a peak modulation frequency.

以上のようにして、本実施の形態では、被測定物11内の特定成分に対する濃度測定感度の、音響モードによる変化を補償することができ、測定感度を向上させることができる。   As described above, in the present embodiment, it is possible to compensate for the change in the concentration measurement sensitivity with respect to the specific component in the DUT 11 due to the acoustic mode, and to improve the measurement sensitivity.

[第2の実施の形態]
次に、本発明の第2の実施の形態について説明する。図10は本発明の第2の実施の形態に係る成分濃度測定装置の構成を示すブロック図であり、図3と同一の構成には同一の符号を付してある。本実施の形態の成分濃度測定装置は、レーザダイオード1−1,1−2と、レーザドライバ2と、光ファイバ3−1,3−2と、光カプラ4と、光ファイバ5と、音響センサ6と、増幅器7と、ファンクションジェネレータ8と、ロックインアンプ9と、情報処理装置10aと、被測定物11の温度を変化させる温度制御器12とから構成される。
[Second Embodiment]
Next, a second embodiment of the present invention will be described. FIG. 10 is a block diagram showing a configuration of a component concentration measuring apparatus according to the second embodiment of the present invention. The same components as those in FIG. 3 are denoted by the same reference numerals. The component concentration measuring apparatus according to the present embodiment includes laser diodes 1-1 and 1-2, a laser driver 2, optical fibers 3-1 and 3-2, an optical coupler 4, an optical fiber 5, and an acoustic sensor. 6, an amplifier 7, a function generator 8, a lock-in amplifier 9, an information processing device 10 a, and a temperature controller 12 that changes the temperature of the DUT 11.

図11は本実施の形態の情報処理装置10aの構成を示すブロック図である。情報処理装置10aは、ジェネレータ制御部100と、光強度制御部101と、光強度測定部102と、光強度変化量算出部103と、音響モード推定部104と、補正係数決定部105と、濃度導出部106と、記憶部107と、被測定物11の温度を制御する温度制御部108とを有する。   FIG. 11 is a block diagram showing a configuration of the information processing apparatus 10a of the present embodiment. The information processing apparatus 10a includes a generator control unit 100, a light intensity control unit 101, a light intensity measurement unit 102, a light intensity change amount calculation unit 103, an acoustic mode estimation unit 104, a correction coefficient determination unit 105, a concentration A derivation unit 106, a storage unit 107, and a temperature control unit 108 that controls the temperature of the DUT 11 are provided.

次に、本実施の形態の成分濃度測定装置の動作を図12のフローチャートを参照して説明する。初めに、レーザダイオード1−1(第1の光源)のみを動作させる。レーザドライバ2から駆動電流が供給されると、レーザダイオード1−1はレーザ光(連続光)を放射する。このとき、レーザドライバ2から矩形波の駆動電流が供給されることにより、レーザダイオード1−1は強度変調光を放射する。第1の実施の形態と同様に、レーザダイオード1−1から放射される光の波長(第1の波長)は例えば1382nmである。レーザダイオード1−1から放射された強度変調光は、光ファイバ3−1によって導かれ、光カプラ4および光ファイバ5を介して被測定物11に照射される(図12ステップS300)。   Next, the operation of the component concentration measuring apparatus of the present embodiment will be described with reference to the flowchart of FIG. First, only the laser diode 1-1 (first light source) is operated. When a drive current is supplied from the laser driver 2, the laser diode 1-1 emits laser light (continuous light). At this time, when a rectangular-wave drive current is supplied from the laser driver 2, the laser diode 1-1 emits intensity-modulated light. Similar to the first embodiment, the wavelength of light emitted from the laser diode 1-1 (first wavelength) is, for example, 1382 nm. The intensity-modulated light emitted from the laser diode 1-1 is guided by the optical fiber 3-1, and irradiated onto the device under test 11 through the optical coupler 4 and the optical fiber 5 (step S300 in FIG. 12).

音響センサ6は、被測定物11から発生する光音響信号を検出し、増幅器7は、音響センサ8から出力された電気信号を増幅する。ロックインアンプ9は、増幅器7の出力に含まれる信号のうち、ファンクションジェネレータ8から出力される参照信号によって決まる周波数の測定信号を検出する。   The acoustic sensor 6 detects a photoacoustic signal generated from the device under test 11, and the amplifier 7 amplifies the electrical signal output from the acoustic sensor 8. The lock-in amplifier 9 detects a measurement signal having a frequency determined by the reference signal output from the function generator 8 among the signals included in the output of the amplifier 7.

情報処理装置10aのジェネレータ制御部100は、ファンクションジェネレータ8が発生する参照信号の周波数を変化させることにより、レーザドライバ2からレーザダイオード1−1に供給される駆動電流の周波数を変化させ、光変調周波数を漸次変化させると共に、ロックインアンプ11が検出する測定信号の周波数を漸次変化させる光変調周波数掃引を行う(図12ステップS301)。こうして、変調周波数と測定信号の振幅(音圧)との関係を取得し(図12ステップS302)、測定信号の振幅が最大の変調周波数であるピーク変調周波数を探索する(図12ステップS303)。   The generator control unit 100 of the information processing apparatus 10a changes the frequency of the reference current generated by the function generator 8 to change the frequency of the drive current supplied from the laser driver 2 to the laser diode 1-1, thereby performing optical modulation. The optical modulation frequency sweep is performed to gradually change the frequency and gradually change the frequency of the measurement signal detected by the lock-in amplifier 11 (step S301 in FIG. 12). In this way, the relationship between the modulation frequency and the amplitude (sound pressure) of the measurement signal is acquired (step S302 in FIG. 12), and the peak modulation frequency that is the modulation frequency with the maximum amplitude of the measurement signal is searched (step S303 in FIG. 12).

次に、レーザダイオード1−1を再び動作させて連続光を出力させ(図12ステップS304)、続いてレーザダイオード1−2(第2の光源)を動作させて連続光を出力させる(図12ステップS305)。レーザダイオード1−1から放射される光の波長は例えば1382nm、レーザダイオード1−2から放射される光の波長(第2の波長)は例えば1610nmである。   Next, the laser diode 1-1 is operated again to output continuous light (FIG. 12, step S304), and then the laser diode 1-2 (second light source) is operated to output continuous light (FIG. 12). Step S305). The wavelength of the light emitted from the laser diode 1-1 is, for example, 1382 nm, and the wavelength of the light emitted from the laser diode 1-2 (second wavelength) is, for example, 1610 nm.

情報処理装置10aのジェネレータ制御部100は、ファンクションジェネレータ8が発生する2つの参照信号の周波数を変化させることにより、レーザドライバ2からレーザダイオード1−1,1−2に供給される駆動電流の周波数を変化させ、光変調周波数を上記のピーク変調周波数に設定すると共に、ロックインアンプ9が検出する測定信号の周波数をピーク変調周波数に設定する。このとき、レーザドライバ2から、同一周波数で逆位相の矩形波の駆動電流をレーザダイオード1−1,1−2に供給させることにより、レーザダイオード1−1,1−2から放射される光を同一周波数(ピーク変調周波数)で逆位相の信号によりそれぞれ強度変調する(図12ステップS306)。   The generator control unit 100 of the information processing apparatus 10a changes the frequency of the two reference signals generated by the function generator 8 to thereby change the frequency of the drive current supplied from the laser driver 2 to the laser diodes 1-1 and 1-2. And the optical modulation frequency is set to the peak modulation frequency, and the frequency of the measurement signal detected by the lock-in amplifier 9 is set to the peak modulation frequency. At this time, light emitted from the laser diodes 1-1 and 1-2 is supplied from the laser driver 2 to the laser diodes 1-1 and 1-2 by supplying a rectangular-wave drive current of the same frequency and opposite phase. Intensity modulation is performed with signals having the same frequency (peak modulation frequency) and opposite phases (step S306 in FIG. 12).

レーザダイオード1−1,1−2から放射された強度変調光は、それぞれ光ファイバ3−1,3−2によって導かれ、光カプラ4によって合波され、さらに光ファイバ5によって導かれ、被測定物11に照射される。   The intensity-modulated lights emitted from the laser diodes 1-1 and 1-2 are respectively guided by the optical fibers 3-1 and 3-2, combined by the optical coupler 4, and further guided by the optical fiber 5 to be measured. The object 11 is irradiated.

情報処理装置10aの光強度制御部101は、ファンクションジェネレータ8を制御して、レーザダイオード1−1を駆動するための参照信号の電圧を所定値に設定することにより、レーザドライバ2からレーザダイオード1−1に供給される駆動電流の大きさを一定にし、レーザダイオード1−1から放射される光の強度を一定にする(図12ステップS307)。   The light intensity control unit 101 of the information processing apparatus 10a controls the function generator 8 to set the voltage of the reference signal for driving the laser diode 1-1 to a predetermined value, so that the laser driver 1 sends the laser diode 1 to the laser diode 1. -1 is made constant, and the intensity of light emitted from the laser diode 1-1 is made constant (step S307 in FIG. 12).

また、光強度制御部101は、ファンクションジェネレータ8を制御して、レーザダイオード1−2を駆動するための参照信号の電圧を変化させることにより、レーザドライバ2からレーザダイオード1−2に供給される駆動電流の大きさを変化させ、レーザダイオード1−2から放射される光の強度を漸次変化させる光強度掃引を行い、ロックインアンプ9から出力される測定信号の振幅が最小になるようにレーザダイオード1−2から放射される光の強度を調節する(図12ステップS308)。   Further, the light intensity controller 101 is supplied from the laser driver 2 to the laser diode 1-2 by controlling the function generator 8 to change the voltage of the reference signal for driving the laser diode 1-2. The intensity of the measurement signal output from the lock-in amplifier 9 is minimized by performing a light intensity sweep that changes the magnitude of the drive current and gradually changes the intensity of the light emitted from the laser diode 1-2. The intensity of light emitted from the diode 1-2 is adjusted (step S308 in FIG. 12).

情報処理装置10aの光強度測定部102は、第1の実施の形態と同様に、レーザダイオード1−2から放射される光の強度を測定する。そして、光強度測定部102は、測定信号の振幅が最小となるときの、レーザダイオード1−2から放射される光の強度I0を記憶部107に格納する(図12ステップS309)。 The light intensity measurement unit 102 of the information processing apparatus 10a measures the intensity of light emitted from the laser diode 1-2, as in the first embodiment. Then, the light intensity measurement unit 102 stores the intensity I 0 of the light emitted from the laser diode 1-2 when the amplitude of the measurement signal is minimum in the storage unit 107 (step S309 in FIG. 12).

次に、情報処理装置10aの温度制御部108は、温度制御器12を制御して被測定物11の温度をΔT(例えば1℃)だけ上げる(図12ステップS310)。   Next, the temperature control unit 108 of the information processing apparatus 10a controls the temperature controller 12 to increase the temperature of the DUT 11 by ΔT (for example, 1 ° C.) (step S310 in FIG. 12).

情報処理装置10aの光強度制御部101は、ファンクションジェネレータ8を制御して、レーザダイオード1−2を駆動するための参照信号の電圧を変化させることにより、レーザドライバ2からレーザダイオード1−2に供給される駆動電流の大きさを変化させ、レーザダイオード1−2から放射される光の強度を漸次変化させる光強度掃引を行い、ロックインアンプ9から出力される測定信号の振幅が最小になるようにレーザダイオード1−2から放射される光の強度を調節する(図12ステップS311)。上記のとおり、レーザダイオード1−1から放射される光の強度は一定である。   The light intensity control unit 101 of the information processing apparatus 10a controls the function generator 8 to change the voltage of the reference signal for driving the laser diode 1-2, thereby changing the laser driver 2 to the laser diode 1-2. The intensity of the measurement signal output from the lock-in amplifier 9 is minimized by changing the magnitude of the supplied drive current and performing a light intensity sweep that gradually changes the intensity of the light emitted from the laser diode 1-2. In this way, the intensity of the light emitted from the laser diode 1-2 is adjusted (step S311 in FIG. 12). As described above, the intensity of light emitted from the laser diode 1-1 is constant.

情報処理装置10aの光強度測定部102は、レーザダイオード1−2から放射される光の強度をステップS309と同様に測定する。そして、光強度測定部102は、測定信号の振幅が最小となるときの、レーザダイオード1−2から放射される光の強度I1を記憶部107に格納する(図12ステップS312)。 The light intensity measurement unit 102 of the information processing apparatus 10a measures the intensity of light emitted from the laser diode 1-2 in the same manner as in step S309. Then, the light intensity measurement unit 102 stores the intensity I 1 of the light emitted from the laser diode 1-2 when the amplitude of the measurement signal is minimum in the storage unit 107 (step S312 in FIG. 12).

情報処理装置10aの光強度変化量算出部103は、ステップS309,S312で求めた光強度I0,I1より光強度I0に対する光強度の変化量ΔI=(I1−I0)の比ΔI/Iを、式(15)により計算する(図12ステップS313)。被測定物11の温度を1℃上げた場合の温度変化前後の光強度の変化量の比ΔI/Iを図13に示す。図13によると、ΔI/I=1.8%という結果が得られる。 The light intensity change amount calculation unit 103 of the information processing apparatus 10a calculates the ratio of the light intensity change amount ΔI = (I 1 −I 0 ) to the light intensity I 0 from the light intensities I 0 and I 1 obtained in steps S309 and S312. ΔI / I is calculated by equation (15) (step S313 in FIG. 12). FIG. 13 shows the ratio ΔI / I of the amount of change in the light intensity before and after the temperature change when the temperature of the DUT 11 is raised by 1 ° C. According to FIG. 13, a result of ΔI / I = 1.8% is obtained.

次に、情報処理装置10aの音響モード推定部104は、表2のような所定のデータセットを用いて、光強度の変化量の比ΔI2/I2を式(16)、式(17)により音響モードq毎に計算する(図12ステップS314)。データセットは、第1、第2、第3の波長と、これらの波長の光に対する特定成分(例えばグルコース)の光吸収係数と、これら光吸収係数の温度依存性との関係を表すものである。 Next, the acoustic mode estimation unit 104 of the information processing apparatus 10a uses the predetermined data set as shown in Table 2 to calculate the light intensity change amount ratio ΔI 2 / I 2 using the equations (16) and (17). Is calculated for each acoustic mode q (step S314 in FIG. 12). The data set represents the relationship between the first, second, and third wavelengths, the light absorption coefficient of a specific component (for example, glucose) with respect to light of these wavelengths, and the temperature dependence of these light absorption coefficients. .

γ1は第1の波長の光に対する被測定物の光吸収係数α1の温度依存性(温度吸光度)、γ2は第2の波長の光に対する被測定物の光吸収係数α2の温度依存性である。関数F1は第1の波長の光に関する光吸収係数α1=0.30/mmと音響モードqから計算することができ、関数F2は第2の波長の光に関する光吸収係数α2=0.28/mmと音響モードqから計算することができる。被測定物11の長さLは既知の値である。 γ 1 is the temperature dependence of the light absorption coefficient α 1 of the object to be measured for the first wavelength light (temperature absorbance), and γ 2 is the temperature dependence of the light absorption coefficient α 2 of the object to be measured for the second wavelength light. It is sex. The function F 1 can be calculated from the light absorption coefficient α 1 = 0.30 / mm for the first wavelength light and the acoustic mode q, and the function F 2 is the light absorption coefficient α 2 = for the second wavelength light. It can be calculated from 0.28 / mm and the acoustic mode q. The length L of the DUT 11 is a known value.

ステップS314の計算の結果から、図14に示すような音響モードqと光強度の変化量の比ΔI2/I2との関係が得られる。情報処理装置10aの音響モード推定部104は、ステップS313で計算された光強度の変化量の比ΔI/Iと、ステップS314で計算した光強度の変化量の比ΔI2/I2とが等しくなる音響モードqの値を、今回のピーク変調周波数での音響モードの値として決定する(図12ステップS315)。本実施の形態の例では、q=1であることが推測できる。 The relationship between the acoustic mode q and the light intensity change ratio ΔI 2 / I 2 as shown in FIG. 14 is obtained from the calculation result of step S314. The acoustic mode estimation unit 104 of the information processing apparatus 10a has the light intensity change ratio ΔI / I calculated in step S313 equal to the light intensity change ratio ΔI 2 / I 2 calculated in step S314. The value of the acoustic mode q is determined as the value of the acoustic mode at the current peak modulation frequency (step S315 in FIG. 12). In the example of the present embodiment, it can be estimated that q = 1.

情報処理装置10aの補正係数決定部105は、音響モード推定部104が推定した音響モードqの値に対応する、被測定物11内の特定成分(本実施の形態の例ではグルコース)の濃度の測定感度Gqを、濃度測定感度と音響モードqとの既知の関係(図9)から求め、この濃度測定感度Gqに対する音響モードq=0のときの濃度測定感度G0の比G0/Gqを補正係数Aとして決定する(図12ステップS316)。   The correction coefficient determination unit 105 of the information processing device 10a has a concentration of a specific component (glucose in the example of the present embodiment) in the DUT 11 corresponding to the value of the acoustic mode q estimated by the acoustic mode estimation unit 104. The measurement sensitivity Gq is obtained from a known relationship between the density measurement sensitivity and the acoustic mode q (FIG. 9), and the ratio G0 / Gq of the density measurement sensitivity G0 when the acoustic mode q = 0 with respect to the density measurement sensitivity Gq is a correction coefficient. A is determined (step S316 in FIG. 12).

図9によると、音響モードq=1のときの濃度測定感度Gqは音響モードq=0のときの濃度測定感度G0より約3倍大きいので、補正係数決定部105は、補正係数A=1/3と求めることができる。   According to FIG. 9, since the density measurement sensitivity Gq when the acoustic mode q = 1 is approximately three times larger than the density measurement sensitivity G0 when the acoustic mode q = 0, the correction coefficient determination unit 105 determines that the correction coefficient A = 1 / 3 can be obtained.

情報処理装置10aの濃度導出部106は、ステップS313で計算された光強度の変化量の比ΔI/I(音響モードq=1のときの光強度の変化量の比ΔI2/I2)と、補正係数決定部105が決定した補正係数Aとから、測定開始時から任意の時間経過後の時点における被測定物11内の特定成分の濃度の変化量ΔC[%]を式(18)により計算する(図12ステップS317)。 The concentration deriving unit 106 of the information processing apparatus 10a calculates the light intensity change amount ratio ΔI / I (the light intensity change amount ratio ΔI 2 / I 2 when the acoustic mode q = 1) calculated in step S313. Based on the correction coefficient A determined by the correction coefficient determination unit 105, the change amount ΔC [%] of the concentration of the specific component in the DUT 11 at an arbitrary time after the start of the measurement is expressed by Expression (18). Calculation is performed (step S317 in FIG. 12).

最後に、濃度導出部106は、ステップS317で計算した特定成分の濃度の変化量ΔC[%]と特定成分の既知の参照濃度とから、任意の時間経過後の時点における特定成分の濃度を計算する(図12ステップS318)。
こうして、本実施の形態では、上記の第1の実施の形態と同様の効果を得ることができる。
Finally, the concentration deriving unit 106 calculates the concentration of the specific component at an arbitrary time point after the elapse of an arbitrary time from the change amount ΔC [%] of the concentration of the specific component calculated in step S317 and the known reference concentration of the specific component. (Step S318 in FIG. 12).
Thus, in this embodiment, the same effect as that of the first embodiment can be obtained.

音波の反射伝搬は、音響インピーダンスを用いて評価することができ、例えば異種材料界面における反射率Rおよび透過率Tは下式で表すことができる。   The reflected propagation of sound waves can be evaluated using acoustic impedance. For example, the reflectance R and transmittance T at the interface between different materials can be expressed by the following equations.

1は第1の材料の音響インピーダンス、Z2は第1の材料と接する第2の材料の音響インピーダンスである。音響インピーダンスZは式(21)で表すことができる。 Z 1 is the acoustic impedance of the first material, and Z 2 is the acoustic impedance of the second material in contact with the first material. The acoustic impedance Z can be expressed by Equation (21).

ρは材料の密度、Cは材料内での音速である。音響インピーダンスの差によって、反射、透過する音波が決定される。図15に音波の反射率と音響インピーダンスとの関係を示す。被測定物11と共振器材料界面との反射率が0.8以上であることが共振器(被測定物11を収容する光音響セル)として機能することの必要条件である。この条件を満たす材料として、たとえばガラスや鉄や銅をはじめとした各種金属がある。照射する光により生じる音波の周波数は、式(22)のように記述することができる。共振器は、照射される光の光軸に対して略直交する面を2つ持ち、これらの2つの面は互いに平行平面となっていることが好ましい。   ρ is the density of the material, and C is the speed of sound within the material. The sound wave to be reflected and transmitted is determined by the difference in acoustic impedance. FIG. 15 shows the relationship between sound wave reflectance and acoustic impedance. A reflectance between the DUT 11 and the resonator material interface being 0.8 or more is a necessary condition for functioning as a resonator (a photoacoustic cell accommodating the DUT 11). Examples of materials that satisfy this condition include various metals such as glass, iron, and copper. The frequency of the sound wave generated by the irradiating light can be described as in Expression (22). It is preferable that the resonator has two planes substantially orthogonal to the optical axis of the irradiated light, and these two planes are parallel to each other.

λは音波の波長、Cは被測定物11内での音速、fは光の変調周波数である。本発明の光音響を用いた測定方法においては、光の変調周波数を100kHz−1MHzとした場合、この時の音波の波長は1mm弱から12mmとなる。たとえば、軸方向2次のモードを得るためには、共振長はおよそ波長程度に設定する必要がある。   λ is the wavelength of the sound wave, C is the speed of sound in the DUT 11, and f is the modulation frequency of the light. In the measurement method using the photoacoustic of the present invention, when the modulation frequency of light is 100 kHz-1 MHz, the wavelength of the sound wave at this time is less than 1 mm to 12 mm. For example, in order to obtain a secondary mode in the axial direction, the resonance length needs to be set to about the wavelength.

なお、第1、第2の実施の形態において、被測定物11に照射される複数の光は、これら複数の光の波長に対する、被測定物11に含まれる測定対象の成分以外の成分の光吸収係数が略等しいことが好ましい。   In the first and second embodiments, the plurality of lights irradiated on the object to be measured 11 is light of components other than the component to be measured included in the object to be measured 11 with respect to the wavelengths of the plurality of lights. It is preferable that the absorption coefficients are substantially equal.

また、第1、第2の実施の形態において、異なる波長の複数の光を光学系(不図示)によりそれぞれ平行光にして被測定物11の略同一領域に同軸で照射することが好ましい。また、このとき異なる波長の複数の光のビーム径は略等しいことが好ましい。   In the first and second embodiments, it is preferable that a plurality of lights having different wavelengths are converted into parallel lights by an optical system (not shown) and irradiated on the substantially same region of the DUT 11 coaxially. At this time, it is preferable that the beam diameters of a plurality of lights having different wavelengths are substantially equal.

第1、第2の実施の形態の情報処理装置10,10aは、例えばCPU(Central Processing Unit)、記憶装置およびインタフェースを備えたコンピュータとこれらのハードウェア資源を制御するプログラムによって実現することができる。CPUは、記憶装置に格納されたプログラムに従って第1、第2の実施の形態で説明した処理を実行する。   The information processing apparatuses 10 and 10a according to the first and second embodiments can be realized by, for example, a computer having a CPU (Central Processing Unit), a storage device, and an interface, and a program that controls these hardware resources. . The CPU executes the processing described in the first and second embodiments in accordance with a program stored in the storage device.

本発明は、血液グルコース、アルブミン等の成分の濃度をモニターする技術に適用することができる。   The present invention can be applied to a technique for monitoring the concentration of components such as blood glucose and albumin.

1−1〜1−3…レーザダイオード、2…レーザドライバ、3−1〜3−3,5…光ファイバ、4…光カプラ、6…音響センサ、7…増幅器、8…ファンクションジェネレータ、9…ロックインアンプ、10,10a…情報処理装置、11…被測定物、12…温度制御器、100…ジェネレータ制御部、101…光強度制御部、102…光強度測定部、103…光強度変化量算出部、104…音響モード推定部、105…補正係数決定部、106…濃度導出部、107…記憶部、108…温度制御部。   DESCRIPTION OF SYMBOLS 1-1 to 1-3 ... Laser diode, 2 ... Laser driver, 3-1 to 3-3, 5 ... Optical fiber, 4 ... Optical coupler, 6 ... Acoustic sensor, 7 ... Amplifier, 8 ... Function generator, 9 ... Lock-in amplifier, 10, 10a ... Information processing device, 11 ... Object to be measured, 12 ... Temperature controller, 100 ... Generator control unit, 101 ... Light intensity control unit, 102 ... Light intensity measurement unit, 103 ... Light intensity change amount Calculation unit 104 ... acoustic mode estimation unit 105 ... correction coefficient determination unit 106 ... concentration derivation unit 107 ... storage unit 108 ... temperature control unit

Claims (8)

互いに異なる波長の複数の光を同一の周波数で且つ異なる位相の信号によりそれぞれ強度変調して被測定物に照射する光照射手段と、
前記複数の強度変調光のうち少なくとも1つの光の強度を変化させる光強度制御手段と、
光照射によって前記被測定物から発生する光音響信号を検出して電気信号を出力する光音響信号検出手段と、
前記光強度制御手段が変化させた光の強度を測定する光強度測定手段と、
強度を変化させる光の波長が異なる2つの時点または前記被測定物の温度が異なる2つの時点について、それぞれ前記電気信号の強度が最低となったときの前記光強度測定手段の測定結果から得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記光音響信号の共振モードである音響モードを推定する音響モード推定手段と、
この音響モード推定手段が推定した音響モードと、前記被測定物に含まれる測定対象の成分に対する濃度測定感度と音響モードとの既知の関係に基づいて、前記測定対象の成分に対する濃度測定感度の補正係数を決定する補正係数決定手段と、
前記光強度測定手段の測定結果から得られる光強度の変化量と前記補正係数に基づいて前記測定対象の成分の濃度を導出する濃度導出手段とを備えることを特徴とする成分濃度測定装置。
A light irradiation means for irradiating the object to be measured by modulating the intensity of a plurality of lights having different wavelengths with signals of the same frequency and different phases,
Light intensity control means for changing the intensity of at least one of the plurality of intensity-modulated lights;
Photoacoustic signal detection means for detecting a photoacoustic signal generated from the object to be measured by light irradiation and outputting an electrical signal; and
A light intensity measuring means for measuring the intensity of the light changed by the light intensity control means;
It is obtained from the measurement result of the light intensity measuring means when the intensity of the electric signal is the lowest for two time points with different wavelengths of light whose intensity changes or two time points with different temperatures of the object to be measured. An acoustic mode estimation means for estimating an acoustic mode that is a resonance mode of the photoacoustic signal based on a variation in light intensity and a variation in light intensity obtained from a known data set;
Based on the known relationship between the acoustic mode estimated by the acoustic mode estimation means, the concentration measurement sensitivity for the measurement target component contained in the object to be measured, and the acoustic mode, the concentration measurement sensitivity correction for the measurement target component is performed. Correction coefficient determining means for determining a coefficient;
A component concentration measuring apparatus comprising: a concentration deriving unit for deriving a concentration of the component to be measured based on a change amount of light intensity obtained from a measurement result of the light intensity measuring unit and the correction coefficient.
請求項1記載の成分濃度測定装置において、
前記光照射手段は、異なる2つの波長の光を前記被測定物に同時に照射し、
前記音響モード推定手段は、第1の波長の光の強度を一定にして第2の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第2の波長の光の強度の測定結果と、前記第1の波長の光の強度を一定にして第3の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第3の波長の光の強度の測定結果とから得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記音響モードを推定することを特徴とする成分濃度測定装置。
In the component concentration measuring apparatus according to claim 1,
The light irradiating means simultaneously irradiates the object to be measured with light of two different wavelengths,
The acoustic mode estimation means changes the intensity of the light of the second wavelength while keeping the intensity of the light of the first wavelength constant, so that the light of the second wavelength when the intensity of the electrical signal becomes the minimum And the third wavelength when the intensity of the electric signal becomes the minimum by changing the intensity of the light of the third wavelength while keeping the intensity of the light of the first wavelength constant. A component concentration measuring apparatus that estimates the acoustic mode based on a change amount of light intensity obtained from a measurement result of light intensity and a change amount of light intensity obtained from a known data set.
請求項1記載の成分濃度測定装置において、
さらに、前記被測定物の温度を制御する温度制御手段を備え、
前記光照射手段は、異なる2つの波長の光を前記被測定物に同時に照射し、
前記音響モード推定手段は、第1の波長の光の強度を一定にして第2の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第2の波長の光の強度の測定結果と、前記温度制御手段により前記被測定物の温度を変化させた後に前記第1の波長の光の強度を一定にして前記第2の波長の光の強度を変化させて前記電気信号の強度が最低となったときの前記第2の波長の光の強度の測定結果とから得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記音響モードを推定することを特徴とする成分濃度測定装置。
In the component concentration measuring apparatus according to claim 1,
Furthermore, a temperature control means for controlling the temperature of the object to be measured is provided,
The light irradiating means simultaneously irradiates the object to be measured with light of two different wavelengths,
The acoustic mode estimation means changes the intensity of the light of the second wavelength while keeping the intensity of the light of the first wavelength constant, so that the light of the second wavelength when the intensity of the electrical signal becomes the minimum And the intensity of the light of the first wavelength is made constant after changing the temperature of the object to be measured by the temperature control means, and the intensity of the light of the second wavelength is changed. Based on the change amount of the light intensity obtained from the measurement result of the light intensity of the second wavelength when the intensity of the electric signal is the lowest, and the change amount of the light intensity obtained from the known data set, A component concentration measuring apparatus for estimating the acoustic mode.
請求項1乃至3のいずれか1項に記載の成分濃度測定装置において、
前記被測定物に照射される複数の光は、これら複数の光の波長に対する、前記被測定物に含まれる測定対象の成分以外の成分の光吸収係数が略等しいことを特徴とする成分濃度測定装置。
In the component concentration measuring apparatus according to any one of claims 1 to 3,
The plurality of lights irradiated on the object to be measured have a component concentration measurement characterized in that light absorption coefficients of components other than the component to be measured included in the object to be measured are substantially equal to the wavelengths of the plurality of lights. apparatus.
請求項1乃至4のいずれか1項に記載の成分濃度測定装置において、
前記光照射手段は、異なる波長の複数の光をそれぞれ平行光にして前記被測定物の略同一領域に同軸で照射し、
異なる波長の複数の光のビーム径は略等しいことを特徴とする成分濃度測定装置。
In the component concentration measuring apparatus according to any one of claims 1 to 4,
The light irradiating means coaxially irradiates substantially the same region of the object to be measured with each of a plurality of lights having different wavelengths as parallel light,
A component concentration measuring apparatus, wherein the beam diameters of a plurality of lights having different wavelengths are substantially equal.
請求項1乃至5のいずれか1項に記載の成分濃度測定装置において、
さらに、前記被測定物を収容し、前記光音響信号を増幅するための共振器を備えることを特徴とする成分濃度測定装置。
In the component concentration measuring apparatus according to any one of claims 1 to 5,
Furthermore, a component concentration measuring apparatus comprising a resonator for accommodating the object to be measured and amplifying the photoacoustic signal.
請求項6記載の成分濃度測定装置において、
前記共振器は、前記光照射手段から照射される光の光軸に対して略直交する面を2つ持ち、これらの面は互いに平行平面となっていることを特徴とする成分濃度測定装置。
In the component concentration measuring apparatus according to claim 6,
The resonator has two surfaces substantially perpendicular to the optical axis of the light emitted from the light irradiation means, and these surfaces are parallel to each other.
互いに異なる波長の複数の光を同一の周波数で且つ異なる位相の信号によりそれぞれ強度変調して被測定物に照射する光照射ステップと、
前記複数の強度変調光のうち少なくとも1つの光の強度を変化させる光強度制御ステップと、
光照射によって前記被測定物から発生する光音響信号を検出して電気信号を出力する光音響信号検出ステップと、
前記光強度制御ステップで変化させた光の強度を測定する光強度測定ステップと、
強度を変化させる光の波長が異なる2つの時点または前記被測定物の温度が異なる2つの時点について、それぞれ前記電気信号の強度が最低となったときの前記光強度測定ステップの測定結果から得られる光強度の変化量と、既知のデータセットから得られる光強度の変化量に基づいて、前記光音響信号の共振モードである音響モードを推定する音響モード推定ステップと、
この音響モード推定ステップで推定した音響モードと、前記被測定物に含まれる測定対象の成分に対する濃度測定感度と音響モードとの既知の関係に基づいて、前記測定対象の成分に対する濃度測定感度の補正係数を決定する補正係数決定ステップと、
前記光強度測定ステップの測定結果から得られる光強度の変化量と前記補正係数に基づいて前記測定対象の成分の濃度を導出する濃度導出ステップとを含むことを特徴とする成分濃度測定方法。
A light irradiating step of irradiating the object to be measured with intensity modulation of a plurality of lights having different wavelengths with signals of the same frequency and different phases;
A light intensity control step of changing the intensity of at least one of the plurality of intensity-modulated lights;
A photoacoustic signal detection step for detecting a photoacoustic signal generated from the object to be measured by light irradiation and outputting an electrical signal;
A light intensity measurement step for measuring the intensity of the light changed in the light intensity control step;
It is obtained from the measurement result of the light intensity measurement step when the intensity of the electric signal is the lowest for two time points with different wavelengths of light whose intensity changes or two time points with different temperatures of the object to be measured. An acoustic mode estimation step for estimating an acoustic mode, which is a resonance mode of the photoacoustic signal, based on a variation in light intensity and a variation in light intensity obtained from a known data set;
Based on the known relationship between the acoustic mode estimated in this acoustic mode estimation step, the concentration measurement sensitivity for the measurement target component included in the object to be measured, and the acoustic mode, the correction of the concentration measurement sensitivity for the measurement target component is performed. A correction coefficient determination step for determining a coefficient;
A component concentration measurement method, comprising: a concentration derivation step for deriving a concentration of a component to be measured based on a change amount of light intensity obtained from a measurement result of the light intensity measurement step and the correction coefficient.
JP2016113652A 2016-06-07 2016-06-07 Component concentration measuring apparatus and method Active JP6606018B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2016113652A JP6606018B2 (en) 2016-06-07 2016-06-07 Component concentration measuring apparatus and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016113652A JP6606018B2 (en) 2016-06-07 2016-06-07 Component concentration measuring apparatus and method

Publications (2)

Publication Number Publication Date
JP2017217202A true JP2017217202A (en) 2017-12-14
JP6606018B2 JP6606018B2 (en) 2019-11-13

Family

ID=60656912

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016113652A Active JP6606018B2 (en) 2016-06-07 2016-06-07 Component concentration measuring apparatus and method

Country Status (1)

Country Link
JP (1) JP6606018B2 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018171178A (en) * 2017-03-31 2018-11-08 日本電信電話株式会社 Constituent concentration measuring device and method
WO2019181375A1 (en) * 2018-03-23 2019-09-26 日本電信電話株式会社 Component-concentration measuring apparatus and method

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5348002A (en) * 1992-04-23 1994-09-20 Sirraya, Inc. Method and apparatus for material analysis
JP2007229320A (en) * 2006-03-03 2007-09-13 Nippon Telegr & Teleph Corp <Ntt> Component concentration measuring apparatus
JP2007259918A (en) * 2006-03-27 2007-10-11 Nippon Telegr & Teleph Corp <Ntt> Component concentration measuring apparatus
JP2010139510A (en) * 2004-05-06 2010-06-24 Nippon Telegr & Teleph Corp <Ntt> Component concentration measuring device and method of controlling component concentration measuring device
JP2012179212A (en) * 2011-03-01 2012-09-20 Nippon Telegr & Teleph Corp <Ntt> Component concentration measuring method and device
JP2013106874A (en) * 2011-11-24 2013-06-06 Nippon Telegr & Teleph Corp <Ntt> Component concentration measuring method and device
JP2014050563A (en) * 2012-09-07 2014-03-20 Nippon Telegr & Teleph Corp <Ntt> Component concentration measurement method
JP2015031670A (en) * 2013-08-07 2015-02-16 日本電信電話株式会社 Component concentration analyzer and component concentration analytic method

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5348002A (en) * 1992-04-23 1994-09-20 Sirraya, Inc. Method and apparatus for material analysis
JP2010139510A (en) * 2004-05-06 2010-06-24 Nippon Telegr & Teleph Corp <Ntt> Component concentration measuring device and method of controlling component concentration measuring device
JP2007229320A (en) * 2006-03-03 2007-09-13 Nippon Telegr & Teleph Corp <Ntt> Component concentration measuring apparatus
JP2007259918A (en) * 2006-03-27 2007-10-11 Nippon Telegr & Teleph Corp <Ntt> Component concentration measuring apparatus
JP2012179212A (en) * 2011-03-01 2012-09-20 Nippon Telegr & Teleph Corp <Ntt> Component concentration measuring method and device
JP2013106874A (en) * 2011-11-24 2013-06-06 Nippon Telegr & Teleph Corp <Ntt> Component concentration measuring method and device
JP2014050563A (en) * 2012-09-07 2014-03-20 Nippon Telegr & Teleph Corp <Ntt> Component concentration measurement method
JP2015031670A (en) * 2013-08-07 2015-02-16 日本電信電話株式会社 Component concentration analyzer and component concentration analytic method

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018171178A (en) * 2017-03-31 2018-11-08 日本電信電話株式会社 Constituent concentration measuring device and method
WO2019181375A1 (en) * 2018-03-23 2019-09-26 日本電信電話株式会社 Component-concentration measuring apparatus and method
JP2019165982A (en) * 2018-03-23 2019-10-03 日本電信電話株式会社 Component concentration measurement device and method
JP7010103B2 (en) 2018-03-23 2022-01-26 日本電信電話株式会社 Component concentration measuring device and method

Also Published As

Publication number Publication date
JP6606018B2 (en) 2019-11-13

Similar Documents

Publication Publication Date Title
JP5411180B2 (en) Component concentration measuring method and apparatus
JP4742279B2 (en) Temperature measuring apparatus, heat treatment apparatus using the same, and temperature measuring method
JP6606019B2 (en) Component concentration measuring apparatus and method
JP5647092B2 (en) Component concentration measuring method and apparatus
JP6606018B2 (en) Component concentration measuring apparatus and method
KR20120103659A (en) Fluorescence measurement device and fluorescence measurement method
JP5839489B2 (en) Component concentration measurement method
JP2014194344A (en) Method for measurement using terahertz wave
JP6619379B2 (en) Component concentration measuring apparatus and method
WO2012114730A1 (en) Photoacoustic measurement device and photoacoustic signal detection method
JP5947761B2 (en) Component concentration analyzer and component concentration analysis method
JP2008304453A (en) Device and method for measuring temperature
TW201250226A (en) Interferometric biometric sensing apparatus including adjustable coupling and associated methods
JP5336438B2 (en) Component concentration measuring method and apparatus
JP2018013417A (en) Component concentration measuring device and method
WO2019181375A1 (en) Component-concentration measuring apparatus and method
JP2016154607A (en) Constituent concentration measuring apparatus and measuring method
JP2006084392A (en) Online crystal grain diameter measuring device, and measuring method using laser ultrasonic wave
JP6730963B2 (en) Component concentration measuring device and analysis method
JP6404741B2 (en) Component concentration measuring apparatus and measuring method
JP2016158971A (en) Component concentration measuring method
Camou et al. Towards non-invasive and continuous monitoring of blood glucose level based on CW photoacoustics: New concept for selective and sensitive measurements of aqueous glucose
CN214374255U (en) High-precision photoacoustic spectrometer
JP2005300356A (en) Method and apparatus for measuring diameter distribution of crystal grain
JP2016154585A (en) Constituent concentration measuring apparatus and measuring method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180904

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190528

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190611

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190813

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20191015

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20191017

R150 Certificate of patent or registration of utility model

Ref document number: 6606018

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150