JP2015148445A - 波長走査式分析装置及び方法 - Google Patents

波長走査式分析装置及び方法 Download PDF

Info

Publication number
JP2015148445A
JP2015148445A JP2014019750A JP2014019750A JP2015148445A JP 2015148445 A JP2015148445 A JP 2015148445A JP 2014019750 A JP2014019750 A JP 2014019750A JP 2014019750 A JP2014019750 A JP 2014019750A JP 2015148445 A JP2015148445 A JP 2015148445A
Authority
JP
Japan
Prior art keywords
component
harmonic
light
measurement target
sample
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
JP2014019750A
Other languages
English (en)
Other versions
JP6262554B2 (ja
Inventor
藤原 真人
Masato Fujiwara
真人 藤原
栄富 盧
Rongfu Lu
栄富 盧
敦夫 渡邉
Atsuo Watanabe
敦夫 渡邉
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.)
Nikkiso Co Ltd
Original Assignee
Nikkiso Co Ltd
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 Nikkiso Co Ltd filed Critical Nikkiso Co Ltd
Priority to JP2014019750A priority Critical patent/JP6262554B2/ja
Priority to PCT/JP2015/052802 priority patent/WO2015119068A1/ja
Publication of JP2015148445A publication Critical patent/JP2015148445A/ja
Application granted granted Critical
Publication of JP6262554B2 publication Critical patent/JP6262554B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/359Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using near infrared light
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/02Details
    • G01J3/06Scanning arrangements arrangements for order-selection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/12Generating the spectrum; Monochromators
    • G01J3/26Generating the spectrum; Monochromators using multiple reflection, e.g. Fabry-Perot interferometer, variable interference filters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J3/42Absorption spectrometry; Double beam spectrometry; Flicker spectrometry; Reflection spectrometry
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/3577Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing liquids, e.g. polluted water
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N2021/3595Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using FTIR
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/39Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using tunable lasers

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

【課題】測定妨害成分の影響を確実に排除し、測定対象成分の測定精度を向上させ、かつ複数の測定対象成分を同時に測定する。
【解決手段】干渉フィルタ12で所定の走査周波数で光源からの光を特定の走査帯域において走査し、試料液に照射する。光学セル18を透過した光を光検出器20で電気信号に変換する。処理装置22は、信号を高速フーリエ変換し、演算して得られた複数の高調波成分の強度と、予め測定されてメモリに記憶された、複数種類の測定対象成分の各々と複数の高調波成分の強度との対応関係を用いて試料中に含有される複数種類の測定対象成分を同時に定量する。
【選択図】図4

Description

本発明は、波長走査式分析装置に関し、測定妨害成分を含む試料中の測定対象成分を光吸収により測定する技術に関する。
従来から、測定妨害成分を含む試料中の測定対象成分、例えば測定妨害成分としての水や電解質を含む試料中のグルコースを測定する場合等において、光吸収を利用する分光分析測定が知られている。
特許文献1には、光吸収の波長が測定対象成分の光吸収の最大吸収波長の近辺に極値を持たないが、その裾が重なっている測定妨害成分を含有する試料中の測定対象成分を光吸収により計測するための装置であって、光源と、最大透過波長を周期的に走査させる手段を持つ干渉フィルタと、検出器と、ゼロクロスディテクタとからなり、前記干渉フィルタの周期的な走査波長を、前記測定対象成分の最大吸収波長を中心とし、かつ前記測定妨害成分の光吸収波長の極値を含まない狭い範囲で走査し、前記光源からの光を前記干渉フィルタに透過させ、その前または後ろに置かれた前記試料に透過または反射させたのち、検出器に入れ電気信号として検出し、電気信号の周期的に変化する成分の信号の全周期及び半周期(各ゼロクロス点間の時間)をゼロクロスディテクタにより求め、全周期に対する半周期の変化量から前記成分を定量するようにした干渉フィルタ透過波長走査式光度計が開示されている。
特許第3532870号
上記の従来技術は、電気信号のゼロクロス点間の時間が測定対象成分の濃度と相関関係にあることを利用して測定対象成分を定量するものであるが、例えば水を含む試料中のグルコースを測定する場合において、水の吸収が小さいスペクトル領域ではグルコースの吸収も小さいため測定することは困難であり、逆に、グルコースの吸収がある領域では水の吸収も大きいのでグルコースのスペクトル成分が水のスペクトル成分に埋もれやすいため同様に測定が困難であり、測定妨害成分を確実に排除することが困難である問題がある。また、測定対象成分が複数存在する場合においても、これら複数の測定対象成分を同時に測定することもできない問題がある。
本発明の目的は、測定妨害成分の影響を確実に排除し、測定対象成分の測定精度を向上させるとともに、複数の測定対象成分を同時に測定することができる分析装置及び方法を提供することにある。
本発明は、試料中の測定対象成分を光吸収により計測する波長走査式分析装置であって、光を所定の走査周波数で走査して分析光として試料に照射する走査手段と、前記試料を透過ないし反射した光を受光して信号に変換する光検出手段と、前記光検出手段からの信号をフーリエ変換して前記走査周波数の高調波成分を演算し、演算して得られた複数の高調波成分の強度に基づき前記試料中に含有される前記測定対象成分を定量する処理手段と、予め測定された、前記測定対象成分と複数の高調波成分の強度との対応関係を記憶する記憶手段とを備え、前記処理手段は、演算して得られた複数の高調波成分の強度と前記対応関係を用いて前記試料中に含有される前記測定対象成分を定量することを特徴とする。本発明において、測定対象成分は、1種類あるいは複数種類のいずれでもよい。
本発明の他の実施形態では、前記測定対象成分がS個(S≧2)の場合、S個の高調波成分の強度に基づいて前記対応関係を用いて前記試料中に含有される前記測定対象成分を定量することを特徴とする。前記測定対象成分が1個の場合、複数の高調波成分の強度に基づいて前記対応関係を用いて高調波成分毎の定量値を算出し、複数の定量値を平均することにより、前記試料中に含有される前記測定対象成分を定量してもよい。
本発明のさらに他の実施形態では、前記分析光は、近赤外光または赤外光であることを特徴とする。
また、本発明は、試料中の測定対象成分を光吸収により計測する波長走査式分析方法であって、光を所定の走査周波数で走査して分析光として試料に照射する走査ステップと、前記試料を透過ないし反射した光を受光して信号に変換する光検出ステップと、前記光検出ステップで得られた信号をフーリエ変換して前記走査周波数の高調波成分を演算し、演算して得られた複数の高調波成分の強度に基づき前記試料中に含有される前記測定対象成分を定量する処理ステップとを備え、前記処理ステップでは、演算して得られた複数の高調波成分の強度と、予め測定されてメモリに記憶された、前記測定対象成分と複数の高調波成分の強度との対応関係を用いて前記試料中に含有される前記測定対象成分を定量することを特徴とする。本発明においても、測定対象成分は、1種類あるいは複数種類のいずれでもよい。
本発明によれば、測定妨害成分の影響を確実に排除し、測定対象成分の測定精度を向上させるとともに、単一あるいは複数の測定対象成分を同時に測定することができる。
電気信号の波形と走査周波数及び高調波成分との関係を示す図である。 吸収スペクトルの形状と高調波成分との関係を示す説明図である。 グルコース及び水の吸収スペクトル図である。 実施形態の構成ブロック図である。 実施形態の光学系の斜視図である。 干渉フィルタの角度と波長との関係を示すグラフ図である。 実施形態の処理フローチャートである。
以下、図面に基づき本発明の実施形態について説明する。
<基本原理>
まず、本実施形態の基本原理について説明する。
図1(a)に、既述した従来技術のように、光源からの光を干渉フィルタに透過させ、干渉フィルタの周期的な走査波長として、測定対象成分の最大吸収波長を中心とし、かつ測定妨害成分の光吸収波長の極値を含まない狭い範囲で走査し、干渉フィルタの後ろに置かれた試料に透過させたのち、検出器に入れ電気信号として検出し、電気信号の周期的に変化する成分の信号を取り出した場合の電気信号波形を示す。図において、波形100は水の信号波形であり、波形200はグルコースの信号波形である。検出器では、これら2つの信号が合成した信号波形として検出され、図1(b)に示す信号波形となる。図において、破線は水の信号波形であり、実線は図1(a)に示す2つの信号波形を合成した信号波形であり、検出器で検出される信号波形である。
従来においては、この実線で示す信号波形のひずみ、すなわちゼロクロス点間の時間がグルコース濃度と相関があることに着目してゼロクロス点間の時間を測定することでグルコースを定量しているが、信号波形のひずみは種々の影響を受けるため正確にこれを測定することは困難である。
そこで、本実施形態では、実線で示す信号波形のゼロクロス点間の時間に着目するのではなく、信号波形を高速フーリエ変換(FFT)し、信号波形の高調波成分(2倍波、3倍波、4倍波、・・・)に着目し、高調波成分の強度がグルコース濃度と相関があることに着目し、高調波成分の強度を測定することでグルコース濃度を定量するものである。
図1(c)に、図1(b)の信号波形を高速フーリエ変換して得られるスペクトルを示す。走査周波数、2倍波、3倍波において強度が極大となるが、本実施形態ではこれらの極値に着目する。走査周波数の強度(1倍波の強度)が測定妨害成分の影響を受けていても、3倍波の強度が測定妨害成分の影響をほとんど受けていない場合、この3倍波の強度に着目してグルコース濃度を定量する。
図2に、試料の吸収スペクトルの形状と、走査の1周期(=1走査)で得られる電気信号との対応関係を示す。左側に示すスペクトル形状をもつ物質を図のように波長走査した場合、透過ないし反射した光の強度(=光検出器からの信号)はそれぞれ右側に示す波形のようになる。物質のスペクトルの変曲点の数が、信号の波の数に反映し、信号を高速フーリエ変換したときに異なった高調波成分として検出される。0変曲点の場合は走査周波数に対応し、1変曲点の場合は走査周波数の2倍波に対応し、2変曲点の場合は走査周波数の3倍波に対応し、3変曲点の場合は走査周波数の4倍波に対応する。一般に、変曲点の数がn個(n=0,1,2,3、・・・)の場合は走査周波数の(n+1)倍波に対応するといえる。
走査帯域において、吸収スペクトルが単調減少または単調増加のスペクトルである場合、得られる電気信号には変曲点はなく、これが走査周波数に対応するものとすると、高調波成分は相対的に小さい。
走査帯域において、吸収スペクトルが1つの変曲点を有するスペクトルである場合、高調波成分のうち2倍波成分が相対的に大きくなり、2倍波成分の強度をグルコースの特徴的な強度としてグルコース濃度を高精度に定量し得る。
走査帯域において、吸収スペクトルが2つの変曲点を有するスペクトルである場合、高調波成分のうち3倍波成分が相対的に大きくなり、3倍波成分の強度をグルコースの特徴的な強度としてグルコース濃度を高精度に定量し得る。
走査帯域において、吸収スペクトルが3つの変曲点を有するスペクトルである場合、高調波成分のうち4倍波成分が相対的に大きくなり、4倍波成分の強度をグルコースを特徴的な強度としてグルコース濃度を高精度に定量し得る。
このようにグルコースの吸収スペクトルの波形に応じて、高速フーリエ変換したときの2倍波、3倍波、4倍波・・・の高調波成分のうち特徴的な高調波成分が定まり、この特徴的な高調波成分の強度からグルコースを定量し得る。言い換えれば、グルコースを測定対象成分として測定する際に、測定妨害成分である水その他の成分に対してその吸収スペクトルの波形が大きく異なる波形となる走査帯域を選択し、そのスペクトルに応じた高調波成分を選択してその強度を測定することで、測定妨害成分である水その他の影響を確実に排除しつつグルコース濃度を測定し得ることになる。例えば、走査帯域として、その走査帯域で水のスペクトルは単調減少し、他方、グルコースのスペクトルが1つの変曲点を有するような場合、走査して得られた電気信号を高速フーリエ変換して2倍波成分を測定し、2倍波成分の強度からグルコース濃度を定量できる。また、走査帯域として、その走査帯域で水のスペクトルは単調減少し、他方、グルコースのスペクトルが2つの変曲点を有するような場合、走査して得られた電気信号を高速フーリエ変換して3倍波成分を測定し、3倍波成分の強度からグルコース濃度を定量できる。ここで、走査帯域として、その走査帯域で水のスペクトルは1つの変曲点を有し、他方、グルコースのスペクトルが2つの変曲点を有するような場合においても、走査1周期で得られた電気信号を高速フーリエ変換して3倍波成分を測定し、3倍波成分の強度からグルコース濃度を定量できる。特定の高調波成分に着目することで、測定妨害物質である水のスペクトルが単調に減少していなくてもグルコース濃度を高精度に定量できる点に留意されたい。
従来のように、得られた電気信号のゼロクロス点間の時間を求める方法では、単調減少スペクトル(変曲点が0)と他の変曲点を有するスペクトルを区別するだけの情報しか得られないところ、本実施形態では、
変曲点が0と変曲点が1
変曲点が0と変曲点が2
変曲点が0と変曲点が3
変曲点が1と変曲点が2
変曲点が1と変曲点が3
変曲点が2と変曲点が3
等のスペクトルを区別することができる。このことは、例えば測定妨害物質の変曲点が0だけでなく、測定妨害物質の変曲点が1であったとしても、変曲点が2、あるいは変曲点が3のスペクトルを用いて測定妨害成分の影響を排除しつつ測定対象成分を定量できることを意味する。本実施形態の分析装置は、単に走査して得られる電気信号の高調波成分を用いて定量するのではなく、測定妨害成分と測定対象成分の組み合わせに応じて着目すべき高調波成分を特定し、特定した高調波成分の強度を用いて定量するものである。従って、測定対象成分が同じであっても、測定妨害成分が変われば着目すべき特徴の高調波成分も変化する。
このように、グルコース等の測定対象成分の吸収スペクトルと水等の測定妨害成分の吸収スペクトルの波形が顕著に相違するような帯域を走査帯域として設定し、かつ、その帯域における測定対象成分のスペクトルに応じた高速フーリエ変換の高調波成分の強度を用いて測定妨害成分の影響を排除しつつ測定対象成分の濃度を定量するものである。従って、どの走査帯域を設定するか、及びどの高調波成分を用いるかは、測定妨害成分及び測定対象成分の吸収スペクトルの相対的な波形の相違に応じて決定される。
図3に、水及びグルコースの吸収スペクトルを示す。図において、Aで示す帯域では水の吸収スペクトルは単調に増加する波形であるのに対し、グルコースの吸収スペクトルは1つの変曲点を有するような波形である。従って、Aで示す帯域を走査帯域として設定し、かつ、走査して得られる電気信号を高速フーリエ変換して高調波成分のうち2倍波成分の強度を測定すれば、2倍波成分をほとんど有しない水の影響を排除してグルコースを定量できる。また、Bで示す帯域では水の吸収スペクトルは単調に増加する波形であるのに対し、グルコースの吸収スペクトルは3つの変曲点を有するような波形である。従って、Bで示す帯域を走査帯域として設定し、かつ、走査して得られる電気信号を高速フーリエ変換して高調波成分のうち4倍波成分の強度を測定すれば、4倍波成分をほとんど有しない水の影響を排除してグルコースを定量できる。
なお、水以外の測定妨害物質がAで示す帯域においてグルコースと同様に1つの変曲点を有するような場合、Aで示す帯域ではなくBで示す帯域を走査帯域として設定し、4倍波成分の強度を測定すればよい。本実施形態では、任意の高調波成分を用いることができるため、測定妨害物質の吸収スペクトルに応じて適応的に走査帯域を設定することができる。このことは、たとえ測定妨害成分の種類が相対的に多くても、測定対象成分を定量するために走査帯域を設定する際の選択肢が多いことを意味し、測定対象成分の定量化を容易化し得ることを意味する。
以上が基本原理であるが、測定対象成分あるいは測定妨害物質によっては、適応的に走査帯域を設定することが困難な場合があり、特に試料中に複数の測定対象成分が含まれる場合に困難となる。複数の測定対象成分が含まれる場合でも、複数回の測定を行い、各測定毎に測定対象成分を1つに限定することで走査帯域を設定することが可能であるが、複数の測定が必要となるため最終的な結果を得るまでに時間を要してしまう。
そこで、本実施形態では、複数の測定対象成分が含まれる場合において、その測定結果は各測定対象成分の高調波成分の線形和となることに鑑み、予め複数の測定対象成分を含む校正テーブルを作成しておき、この校正テーブルを用いて複数の測定対象成分を同時に定量する。校正テーブルは、複数の高調波成分における、複数の測定対象成分の強度を規定したものであり、測定対象成分を含む水溶液を用いて高調波成分を測定してその強度を規格化あるいは正規化したものが好適である。ここで、「規格化」あるいは「正規化」とは、水のみの強度の影響を除去し、かつ、単位濃度あたりの強度であることを意味する。
複数の測定対象成分が含まれる場合、そのうちの特定の測定対象成分を定量する際には他の測定対象成分は当該特定の測定対象成分にとって測定妨害成分として機能することになるが、本実施形態では、これを排除するのではなく、むしろ測定妨害成分も測定の対象と把握してその強度を測定し、これを用いて複数の測定対象成分を同時に定量するものといえる。
<具体的構成>
次に、本実施形態における波長走査式分析装置の具体的な構成について説明する。
図4に、波長走査式分析装置の構成ブロック図を示す。波長走査式分析装置は、白熱電球10、干渉フィルタ12、ステッピングモータ14、モータドライバ16、試料液を入れるための光学セル18、光検出器20、処理装置22及び表示部30を含んで構成される。
白熱電球10は、分析光を干渉フィルタ12に向けて照射する。
干渉フィルタ12は、走査手段として機能し、分析光の光軸に対して傾斜角が変化するように回転振動され、その最大透過波長がシフトする。干渉フィルタ12はステッピングモータ14及びモータドライバ16により回転軸を中心として往復回転駆動される。走査周波数は、例えば5Hzに設定される。
試料液を入れるための光学セル18は、分析光の光軸上に配置される。試料液には、測定対象成分及び測定妨害成分が含まれる。測定対象成分には、例えばグルコース、尿素、乳酸等の複数種類が含まれる。
光検出器20は、分析光の光軸上であって光学セル18の後ろに配置され、光学セル18を透過した光を受光し、電気信号に変換して処理装置22に出力する。光検出器20は、例えば4秒間透過光を検出する。
処理装置22は、コンピュータ等で構成され、AD変換部24と、フーリエ変換(FFT)部26と、濃度計算部28及びメモリ29を含む。
AD変換部24は、光検出器20からの電気信号(アナログ電気信号)を所定のサンプリング周波数でデジタル信号に変換してフーリエ変換部26に出力する。AD変換部24は、例えば10kHzでサンプリングし、40000点のデジタルデータを出力する。
フーリエ変換部26は、デジタル信号を高速フーリエ変換して濃度計算部28に出力する。
濃度計算部28は、高速フーリエ変換で得られた複数の高調波成分の強度から測定対象成分の濃度を算出して表示部30に出力する。濃度計算部28は、予め複数の高調波成分の強度と複数の測定対象成分の濃度との相関関係を校正テーブル29aとしてメモリ29に記憶しており、該メモリ29を参照することで、フーリエ変換部26で得られた高調波成分の強度に対応する複数の測定対象成分を同時に算出する。
図5に、波長走査式分析装置の光学系斜視図を示す。光源10、集光レンズ11、干渉フィルタ12、光学セル18、集光レンズ19及び光検出器20は光軸上に配置される。光源10からの光は集光レンズ11を介して干渉フィルタ12に入射する。
干渉フィルタ12は、ステッピングモータ14(図5では不図示)により軸を中心に所定の周波数で往復回転するように駆動される。干渉フィルタ12は、バンドパスフィルタであり、ピーク透過波長を中心として一定の波長帯域の光のみを透過するとともに、その回転駆動によりピーク透過波長がシフトする。干渉フィルタ12の回転駆動により最大透過波長が順次シフトした光は光学セル18に照射され、光学セル18を透過して集光レンズ19により集光されて光検出器20に入射する。
光検出器20は、入射した光を電気信号に変換して処理装置22に出力する。
図6に、干渉フィルタ12の回転角度と干渉フィルタ12を透過する光の波長との関係を示す。図において、横軸は波長、縦軸は透過率(%)である。角度は、干渉フィルタ12の法線と光軸とのなす角度を示す。従って、角度=0°は、光軸に対して干渉フィルタ12の面が垂直(光軸と干渉フィルタ12の法線が平行)であることを示す。図に示すように、角度が0°、5°、10°、・・・、40°、45°と変化するに従い、干渉フィルタ12を透過する光のピーク透過波長は短波長側に順次シフトしていく。従って、干渉フィルタ12を所定角度範囲で往復回転する、例えば25°から35°まで一定の方向に回転し、その後、35°から25°まで逆方向に回転することで、1周期分の走査が行われる。
干渉フィルタ12をどの範囲で回転させるか、つまりどの帯域を走査帯域として設定するかは、既述したように、測定妨害成分及び測定対象成分の吸収スペクトルの相対的な波形の相違に応じて決定される。測定対象成分が複数存在する場合、走査帯域としては、複数の測定対象成分の吸収スペクトルを全て考慮して決定される。このため、単一の測定対象成分の場合よりも走査帯域は広がることになる。
また、処理装置22の濃度計算部28は、フーリエ変換部26で算出された複数の高調波成分と、予め測定されメモリ29に記憶された校正テーブル29aを用いて複数の測定対象成分の濃度を同時に算出する。
図7に、本実施形態における濃度算出処理フローチャートを示す。まず、測定に先立って、測定対象成分を含む水溶液の信号を取得する(S101)。複数の測定対象成分を例えばグルコース、尿素、乳酸、アスコルビン酸、ビリルビンとすると、グルコースを含む水溶液、尿素を含む水溶液、乳酸を含む水溶液、アスコルビン酸を含む水溶液、ビリルビンを含む水溶液をそれぞれ用意し、各水溶液を測定して1倍波、2倍波、3倍波、4倍波、及び5倍波の高調波成分の強度を算出する。
表1に、このようにして算出された各高調波成分の強度の一例を示す。
Figure 2015148445
なお、表1には、参考のため水のみの高調波成分も併せて示す。
次に、処理装置22は、各測定対象成分と水との差分をその濃度で除算することで、校正テーブル29aを作成し、メモリ29に記憶する(S102)。校正テーブル29aの一例を表2に示す。
Figure 2015148445
表2において、グルコースに関しては、1倍波は、
(26134−25134)/100=10.0
として算出される。2倍波は、
(6157−6425)/100=−2.68
として算出される。また、尿素に関しては、1倍波は、
(25391−25134)/50=5.14
として算出される。
以下、同様であり、ビリルビンに関しては、5倍波は、
(2592−3315)/10=−72.3
として算出される。各測定対象成分と水との差分をその濃度で除算することで、規格化された強度値が得られる。
表2のような校正テーブル29aを作成してメモリ29に記憶した後、測定すべき試料液を測定し(S103)、高調波成分の強度を算出する(S104)。すなわち、AD変換部24は、光検出器20からの電気信号(アナログ電気信号)を所定のサンプリング周波数でデジタル信号に変換してフーリエ変換部26に出力し、フーリエ変換部26は、デジタル信号を高速フーリエ変換して濃度計算部28に出力する。濃度計算部28は、高調波成分として、1倍波、2倍波、3倍波、4倍波、及び5倍波の信号強度を算出する(S104)。次に、算出したこれらの強度と校正テーブル29aを用いて複数の測定対象成分のそれぞれの濃度を算出する(S105)。
具体的には、高調波成分の強度が、例えば表3に示す値であるとする。
Figure 2015148445
これらの高調波成分の強度は、各測定対象成分の高調波成分の線形和として得られたものであるから、以下の5元1次連立方程式が成り立つ。
P1=26521
=25134+C1×(10.0)+C2×(5.14)+C3×(117)+
C4×(21.7)+C5×(−3.30)
P2=8513
=6425+C1×(−2.68)+C2×(28.6)+C3×(−7.10)
+C4×(−14.1)+C5×(−88.0)
P3=3674
=3542+C1×(−4.34)+C2×(−0.78)+C3×
(−13.2)+C4×(214)+C5×(12.3)
P4=25
=12+C1×(1.90)+C2×(36.0)+C3×(4.30)+C4
×(6.90)+C5×(−1.10)
P5=3027
=3315+C1×(2.28)×C2×(−2.00)+C3×(334)+
C4×(60.3)+C5×(−72.3)
ここで、C1,C2,C3,C4,及びC5は、試料中に含まれるグルコース、尿素、乳酸、アスコルビン酸、ビリルビンのそれぞれの濃度である。従って、濃度計算部28は、上記の5元1次連立方程式を解くことで、C1,C2,C3,C4,及びC5を算出し、これにより複数の測定対象成分の濃度を同時に算出することができる。
本実施形態において、走査手段として干渉フィルタ12を用いて光の最大透過波長をシフトさせているが、必ずしも干渉フィルタ12に限定されるものではなく、波長を走査し得る任意の光学機器、例えば分光器や波長可変レーザを用いることもできる。また、使用する光についても、赤外、近赤外、可視光、紫外のいずれも用いることができる。
また、本実施形態において、干渉フィルタ12を往復回転させているが、一般的に波長走査は必ずしも往復である必要はなく一方向の走査で電気信号を検出してもよい。
また、本実施形態では、試料液を透過した光を光検出器20で検出して電気信号を得ているが、試料液を反射した光を光検出器20で検出して電気信号を得てもよい。要するに、測定対象成分及び測定妨害成分の吸収スペクトルが反映されたスペクトルの電気信号を得ればよい。
また、本実施形態において、フーリエ変換として、高速フーリエ変換(FFT)を用いる場合を例に説明したが、これに限るものではなく、その他のフーリエ変換を用いても良い。
さらに、本実施形態では、複数種類の測定対象成分として、グルコース、尿素、乳酸、アスコルビン酸、ビリルビンを例示したが、これに限定されることはなく、他の測定対象成分としてもよい。一般に、測定対象成分の数をSとすると、校正テーブル29aにおける高調波成分の数もSだけ必要である。各測定対象成分の濃度をCi(iはS個の自然数)、校正テーブル29aにおけるi倍波の各測定対象成分の値、つまり測定対象成分のみの水溶液のi倍波の強度と水のみのi倍波の強度の差分を水溶液の濃度で除算した値(規格化された値)をTi、水のみの値をWiとし、濃度未知の試料液を測定して、i倍波の信号強度Piが得られたものとすると、
Pi=Wi+Σ(Ci×Ti)
但し、ΣはS個の自然数の加算
のS元1次連立方程式を解いてCiを算出すればよい。
ここで、上記式において、iは必ずしも連続した自然数である必要はないことに留意されたい。すなわち、表2では、i=1,2,3,4,5であるが、複数の測定対象成分の種類によっては、i=2,3,5,7,9あるいは、i=1,2,4,5,9等を用いてもよい。処理装置22は、予め個々の測定対象成分とその測定に用いる高調波成分との対応関係をメモリ29に記憶しておき、ユーザがS個の測定対象成分を入力した場合に、測定に用いるべきS個の高調波成分を自動的に選択する構成としてもよい。例えば、ユーザが測定対象成分として、グルコース、尿素、アスコルビン酸を入力した場合、測定に用いるべき3個の高調波成分として2倍波、3倍波、4倍波を自動選択する等である。
本実施形態では、測定対象成分が複数種類の場合について説明したが、本発明は測定対象成分が1種類の場合でも適用可能であり、測定対象が1種類の場合、例えば以下の3つの方法のいずれかにより定量することができる。
方法1:
算出された5つの結果を均等に取り扱って平均値を結果とする。
方法2:
より望ましいのは、高調波成分に重みを付けて平均を算出する方法である。例えば、グルコースの1mg/dL当たりのキャリブレーション値が次のような値であった場合、最大の3倍波を重み1とし、それ以外はキャリブレーション値に応じて重みを計算する。
1倍波 12 ・・・重み10%
2倍波 30 ・・・重み25%
3倍波 120 ・・・重み100%
4倍波 100 ・・・重み83%
5倍波 6 ・・・重み5%
3倍波の結果を基準とし、重み83%の分だけ4倍波の結果に近づけ、さらに重み25%の分だけ2倍波の結果に近づける等の補正を繰り返すことで最終的な結果を算出する。
方法3:
方法2のような重み補正を行うことなく、最大の3倍波の結果だけを用いて算出する。
上記実施形態の説明において、測定対象成分がS個(S≧2)のとき、1〜S倍波の高調波成分の強度に基づいて試料中に含有される測定対象成分を定量する場合を例に説明したが、本発明はこれに限定されるものではなく、例えば、2〜S+1倍波の高調波成分を用いてもよいし、3〜S+2倍波の高調波成分を用いてもよい。
さらに、本実施形態において、校正テーブルを用いる算出方法の一例について説明したが、本発明はこの算出方法に限定されるものではなく、他の算出方法でもよい。
10 光源(白熱電球)、12 干渉フィルタ、14 ステッピングモータ、16 モータドライバ、18 光学セル、20 光検出器、22 処理装置、24 AD変換部、26 高速フーリエ変換(FFT)部、28 濃度計算部、29 メモリ、29a 校正テーブル、30 表示部。

Claims (9)

  1. 試料中の測定対象成分を光吸収により計測する波長走査式分析装置であって、
    光を所定の走査周波数で走査して分析光として試料に照射する走査手段と、
    前記試料を透過ないし反射した光を受光して信号に変換する光検出手段と、
    前記光検出手段からの信号をフーリエ変換して前記走査周波数の高調波成分を演算し、演算して得られた複数の高調波成分の強度に基づき前記試料中に含有される前記測定対象成分を定量する処理手段と、
    予め測定された、前記測定対象成分と複数の高調波成分の強度との対応関係を記憶する記憶手段と、
    を備え、前記処理手段は、演算して得られた複数の高調波成分の強度と前記対応関係を用いて前記試料中に含有される前記測定対象成分を定量することを特徴とする波長走査式分析装置。
  2. 前記測定対象成分は複数種類であることを特徴とする請求項1記載の波長走査式分析装置。
  3. 前記測定対象成分がS個(S≧2)の場合、S個の高調波成分の強度に基づいて前記対応関係を用いて前記試料中に含有される前記測定対象成分を定量することを特徴とする請求項2記載の波長走査式分析装置。
  4. 前記S個の高調波成分は、1〜S倍波の高調波成分であることを特徴とする請求項3記載の波長走査式分析装置。
  5. 前記測定対象成分が1個の場合、複数の高調波成分の強度に基づいて前記対応関係を用いて高調波成分毎の定量値を算出し、複数の定量値を平均することにより、前記試料中に含有される前記測定対象成分を定量することを特徴とする請求項1記載の波長走査式分析装置。
  6. 前記高調波成分毎の定量値に対して重み付けをして平均することにより前記試料中に含有される前記測定対象成分を定量することを特徴とする請求項5記載の波長走査式分析装置。
  7. 前記分析光は、近赤外光または赤外光であることを特徴とする請求項1〜6のいずれかに記載の波長走査式分析装置。
  8. 試料中の測定対象成分を光吸収により計測する波長走査式分析方法であって、
    光を所定の走査周波数で走査して分析光として試料に照射する走査ステップと、
    前記試料を透過ないし反射した光を受光して信号に変換する光検出ステップと、
    前記光検出ステップで得られた信号をフーリエ変換して前記走査周波数の高調波成分を演算し、演算して得られた複数の高調波成分の強度に基づき前記試料中に含有される前記測定対象成分を定量する処理ステップと、
    を備え、前記処理ステップでは、演算して得られた複数の高調波成分の強度と、予め測定されてメモリに記憶された、前記測定対象成分と複数の高調波成分の強度との対応関係を用いて前記試料中に含有される前記測定対象成分を定量することを特徴とする波長走査式分析方法。
  9. 前記測定対象成分は複数種類であることを特徴とする請求項8記載の波長走査式分析方法。
JP2014019750A 2014-02-04 2014-02-04 波長走査式分析装置及び方法 Active JP6262554B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2014019750A JP6262554B2 (ja) 2014-02-04 2014-02-04 波長走査式分析装置及び方法
PCT/JP2015/052802 WO2015119068A1 (ja) 2014-02-04 2015-02-02 波長走査式分析装置及び方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2014019750A JP6262554B2 (ja) 2014-02-04 2014-02-04 波長走査式分析装置及び方法

Publications (2)

Publication Number Publication Date
JP2015148445A true JP2015148445A (ja) 2015-08-20
JP6262554B2 JP6262554B2 (ja) 2018-01-17

Family

ID=53777873

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014019750A Active JP6262554B2 (ja) 2014-02-04 2014-02-04 波長走査式分析装置及び方法

Country Status (2)

Country Link
JP (1) JP6262554B2 (ja)
WO (1) WO2015119068A1 (ja)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3532870B2 (ja) * 2000-02-28 2004-05-31 敦夫 渡邉 干渉フィルタ透過波長走査式光度計
WO2006029848A1 (de) * 2004-09-14 2006-03-23 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V Vorrichtung zum messen mindestens einer gaskomponente
JP2010164375A (ja) * 2009-01-14 2010-07-29 Mitsubishi Heavy Ind Ltd 濃度計測装置及び濃度計測方法
JP2013134239A (ja) * 2011-12-27 2013-07-08 Horiba Ltd ガス計測装置およびガス計測装置における状態量の測定方法。

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3532870B2 (ja) * 2000-02-28 2004-05-31 敦夫 渡邉 干渉フィルタ透過波長走査式光度計
WO2006029848A1 (de) * 2004-09-14 2006-03-23 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V Vorrichtung zum messen mindestens einer gaskomponente
JP2008513736A (ja) * 2004-09-14 2008-05-01 フラウンホーファー・ゲゼルシャフト ツア フェルデルンク デア アンゲヴァンテン フォルシュンク エー.ファウ. 1種類以上のガス成分を測定する装置
JP2010164375A (ja) * 2009-01-14 2010-07-29 Mitsubishi Heavy Ind Ltd 濃度計測装置及び濃度計測方法
JP2013134239A (ja) * 2011-12-27 2013-07-08 Horiba Ltd ガス計測装置およびガス計測装置における状態量の測定方法。

Also Published As

Publication number Publication date
JP6262554B2 (ja) 2018-01-17
WO2015119068A1 (ja) 2015-08-13

Similar Documents

Publication Publication Date Title
JP5358426B2 (ja) 精製油劣化度測定装置及び精製油劣化度測定方法
JP2001174404A (ja) 光断層像計測装置および計測方法
KR102099230B1 (ko) 형광분석에 의한 체외 검출 및/또는 정량화를 위한 시스템
CN110160976B (zh) 基于二次谐波干涉包络补偿的气体浓度检测方法和系统
WO2018184262A1 (zh) 用于激光诱导击穿光谱采集的阶梯光谱仪动态校正方法
US10561348B2 (en) Device and method for spectrum analysis and device for blood glucose measurement
WO2015156684A2 (en) Signal harmonic error cancellation method and apparatus
KR20170052256A (ko) 라만 산란을 이용한 물질의 농도 측정 장치 및 방법
JP6454211B2 (ja) 検体分析装置、血液凝固分析装置、検体分析方法、及びコンピュータプログラム
WO2001063249A1 (fr) Photometre a balayage en longueur d'onde d'emission a filtre antiparasite
CN111157115B (zh) 一种水下布里渊散射光谱获取方法及装置
JP6262554B2 (ja) 波長走査式分析装置及び方法
WO2018151150A1 (ja) 散乱体濃度計測装置及びその方法
JP6310249B2 (ja) 波長走査式分析装置及び方法
EP3162283B1 (en) Blood measuring apparatus using spectroscope
TW201124712A (en) A method for measuring the property of the suspended solids and other substances in the fluid
CN103954594B (zh) 不同光电倍增管电压下三维光谱数据的峰值换算方法
JP3532870B2 (ja) 干渉フィルタ透過波長走査式光度計
RU2413930C1 (ru) Способ определения оптических характеристик однородного рассеивающего вещества
CN117849017A (zh) 一种基于拉曼位移的硝酸盐浓度检测方法
CN111157116B (zh) 一种水下布里渊散射光谱测试系统
CN114636688B (zh) 模型校正方法、光谱设备、计算机设备和存储介质
Gributs et al. Multiresolution analysis for quantification of optical properties in scattering media using pulsed photon time-of-flight measurements
JP2007534955A (ja) 自己相関関数を決定する方法
CN114813559A (zh) 一种多维度的电化学信息与多光谱信息结合测量方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160908

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170801

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170906

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: 20171205

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20171214

R150 Certificate of patent or registration of utility model

Ref document number: 6262554

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250