Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP7014701B2 - 光学分析装置、並びに光学分析装置に用いられる機械学習装置及びその方法 - Google Patents
[go: Go Back, main page]

JP7014701B2 - 光学分析装置、並びに光学分析装置に用いられる機械学習装置及びその方法 - Google Patents

光学分析装置、並びに光学分析装置に用いられる機械学習装置及びその方法 Download PDF

Info

Publication number
JP7014701B2
JP7014701B2 JP2018234772A JP2018234772A JP7014701B2 JP 7014701 B2 JP7014701 B2 JP 7014701B2 JP 2018234772 A JP2018234772 A JP 2018234772A JP 2018234772 A JP2018234772 A JP 2018234772A JP 7014701 B2 JP7014701 B2 JP 7014701B2
Authority
JP
Japan
Prior art keywords
output value
photodetector
light
correction
interferogram
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.)
Active
Application number
JP2018234772A
Other languages
English (en)
Other versions
JP2020094990A (ja
Inventor
嘉健 安藤
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.)
Horiba Ltd
Original Assignee
Horiba 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 Horiba Ltd filed Critical Horiba Ltd
Priority to JP2018234772A priority Critical patent/JP7014701B2/ja
Priority to EP19194274.7A priority patent/EP3667294A1/en
Priority to US16/564,861 priority patent/US20200191707A1/en
Publication of JP2020094990A publication Critical patent/JP2020094990A/ja
Application granted granted Critical
Publication of JP7014701B2 publication Critical patent/JP7014701B2/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
    • 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/3504Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis
    • 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/255Details, e.g. use of specially adapted sources, lighting or optical systems
    • 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/27Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
    • G01N21/274Calibration, base line adjustment, drift correction
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/14Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N20/00Machine learning
    • 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
    • G01N2201/00Features of devices classified in G01N21/00
    • G01N2201/12Circuits of general importance; Signal processing
    • G01N2201/121Correction signals
    • G01N2201/1211Correction signals for temperature

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Chemical & Material Sciences (AREA)
  • Biochemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Software Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Spectrometry And Color Measurement (AREA)

Description

本発明は、FTIRなどの光学分析装置、並びに光学分析装置に用いられる機械学習装置及びその方法に関するものである。
FTIRなどの分光分析器は、試料を経た光のスペクトルに基づいて試料成分を分析する。この試料を経た光の強度を検出するには、MCT光検出器などの光検出器が用いられるが、当然のことながら、分析に先立ってその出力値と入射光の強度との関係(検量線)を求めておく必要がある。
そこで、従来は、所定の基準光を光検出器に入射してそのときの光検出器の出力値を取得し、検量線を作成する。
例えばMCT光検出器を例にとった場合、入射光の強度が所定以下では、図1に示すように、その出力値と入射光強度との関係は線形(比例関係)となる。したがって、この領域であれば、検量線は、オフセットとスパンを求めればよいので比較的容易に作成できる。
しかしながら、入射光強度がそれ以上になると、同図に示すように、非線形な関係となり、正確な検量線を求めるためには、異なる強度の多数の基準光を用意し、それぞれの出力値を測定しなければならなくなる。具体的には、1ポイントの測定に10分程度かかり、これを複数ポイント(例えば10ポイント以上)に亘って行わなければならないので、検量線の作成だけで数時間という多大な工数(時間)を要する場合があるうえ、それなりの設備も必要になるという問題が生じる。また、前記工数や設備の問題から、メンテナンスや修理時に、製品納入先での検量線の点検や再作成が難しいという問題も生じる。そして、これは、MCT光検出器のみならず、他の種類の光検出器にも共通する問題である。
特開2017-26506号公報
このような問題を解決すべく、本出願人は、特許文献1に示すように、光路上に光学フィルタを介在させた場合とそうでない場合との双方で、光源から射出された光を光検出器に入射させ、該入射光に含まれる複数の所定波数の光のそれぞれに対する光検出器の出力値である第1出力値及び第2出力値を取得し、所定波数毎の第1出力値及び第2出力値の比と、光学フィルタの光学特性とをパラメータとして、光検出器の出力値から入射光の強度を算出するための演算式(検量線)を求める方法を考えている。
しかしながら、この出力補正方法は、1つの光学フィルタを用いた1点補正であり、光検出器に入射する光量が変わった場合には、それとは異なる光量により得られた演算式から入射光の強度を算出することになってしまい、測定誤差の要因となる。なお、光学特性の異なる複数の光学フィルタを用いて異なる光量毎に演算式を求めることも考えられるが、装置構成が複雑になってしまう。
そこで本発明は、上述した課題を解決すべくなされたものであり、光学分析装置に用いられる光検出器の出力補正を、複数の光量を用いて、かつ、簡単な構成で行うことをその主たる課題とするものである。
すなわち本発明に係る光学分析装置は、光を検出して出力値を得る光検出器の非線形性を補正する光学分析装置であって、前記光を照射する光源と、前記光源から照射された光が入力される測定セルと、前記測定セルから出力される光を検出して出力値を得る前記光検出器と、前記光検出器の補正前の出力値と、前記光検出器の線形領域における出力値又は前記測定セルに導入されたサンプルの既知濃度とを用いて、前記光検出器の出力値を線形補正するための補正用モデルを機械学習する学習部と、前記学習部により得られた補正用モデルを保存する保存部とを備えることを特徴とする。
このようなものであれば、機械学習アルゴリズムにより光検出器の出力値を線形補正するための補正用モデルを学習しているので、従来の補正式演算アルゴリズムを不要とし、光検出器の出力補正を、複数の光量を用いて、かつ、簡単な構成で行うことができる。
前記光学分析装置が、前記光検出器が検出する前の光の強度を変動させる変動部をさらに備えている場合には、前記学習部は、前記光検出器の補正前の出力値に含まれる前記変動部による変動成分と、前記光検出器の線形領域における出力値とに基づいて、前記光検出器の出力値を線形補正するための補正用モデルを学習することが望ましい。
前記光学分析装置は、前記測定セルを複数の温度に温調する温調部をさらに備えている場合には、前記測定セルから出力される光は、温調された前記測定セルから発生し、その温度によって光量が変化する光と、前記光源から照射されて前記測定セルを通過した光とを含むものとなる。この場合、前記学習部は、前記測定セルを複数の温度に温調させたときのそれぞれにおいて前記光検出器により検出された複数の補正前の出力値を用いて、前記補正用モデルを学習することが望ましい。
この構成であれば、従来のように1点補正では無く多点補正することになり、頑健な補正用モデルを作成することができる。
測定セルの赤外輻射はその温度変化に伴って変化し、これによって、光検出器に入射する光量が変化することになる。図2には、測定セルが25.9℃の場合と、測定セルが164.1℃の場合とにおける光検出器からの出力値を示している。なお、この出力値は、光源と測定セルとの間に干渉器、又は、光学チョッパや分光器等の入射光強度を変調させる光変調器を設けた構成により得られたものであり、図2は干渉計の移動鏡位置にしたがって光検出器の出力値が変化する様子を示している。また、両温度において、測定セル内のサンプルは同じであり、温度以外は同条件である。図2から分かるように、測定セルが25.9℃の場合と、測定セルが164.1℃の場合とで、光検出器の出力において測定セルからの赤外輻射に依存する直流オフセット成分(DC成分)が変化してしまう。一方、光検出器の出力特性が線形の場合には、干渉計による変動成分(非DC成分)は、いずれの温度の場合であっても等しくなるはずであるが、実際は164.1℃の場合では非DC成分の振幅が小さくなっている。これは、光検出器の出力特性が非線形であることに起因する。
補正用モデルの頑健性を向上させるためには、前記学習部は、さらに、前記光の強度をゼロとしたときの前記光検出器の出力値、前記光検出器の出力値に含まれる前記測定セルの赤外輻射によるオフセット成分、又は、前記測定セルの温度の少なくとも1つを取得して、前記補正用モデルを学習することが望ましい。
変動部の具体的な実施の態様としては、前記光源及び前記測定セルの間に設けられ、前記光源から照射された光を反射させる固定ミラーと、前記光源から照射された光を往復移動しながら反射する移動ミラーとを備えた干渉計であるか、又は光を通過又は遮断させて前記光の強度を変動させる光変調器であることが考えられる。
また、変動部としては、前記測定セルに導入するサンプルの濃度を変動させ、前記光検出器が検出する前の光の強度を変動させる構成としても良い。
具体的に前記学習部は、強化学習アルゴリズムを用いて補正用モデルを学習するものであり、前記光検出器の補正前の出力値が入力されて、前記補正前の出力値を線形補正するための補正用モデルを演算する強化学習部と、前記補正用モデルを、前記光検出器の線形領域における出力値又は前記測定セルに導入されたサンプルの既知濃度を用いて評価するモデル評価部とを備え、前記強化学習部は、前記モデル評価部の評価結果に基づく報酬を用いて前記補正用モデルを更新することが望ましい。
前記モデル評価部は、前記補正用モデルを用いて補正された補正後の出力値と前記線形領域における出力値との間の時間成分の一致度、周波数成分の一致度又は分光特性の一致度の少なくとも1つにより評価することが望ましい。
学習部は、強化学習アルゴリズムの他に、教師あり学習アルゴリズムを用いたものであっても良い。この場合、前記学習部は、教師データとして前記検出器の線形領域における出力値又はサンプルの既知濃度を用いて前記補正用モデルを学習することが望ましい。
また本発明に係る機械学習装置は、光を検出して出力値を得る光検出器の非線形性を補正する光学分析装置に用いられる機械学習装置であって、前記光検出器の補正前の出力値と、前記光検出器の線形領域における出力値又は前記光学分析装置に導入されたサンプルの既知濃度とを用いて、前記光検出器の出力値を線形補正するための補正用モデルを学習する学習部と、前記学習部により得られた補正用モデルを保存する保存部とを備えることを特徴とする。
さらに本発明に係る機械学習方法は、光を検出して出力値を得る光検出器の非線形性を補正する光検出器を備える光学分析装置に用いられる機械学習方法であって、前記光検出器の補正前の出力値と、前記光検出器の線形領域における出力値又は前記測定セルに導入されたサンプルの既知濃度とを用いて、前記光検出器の出力値を線形補正するための補正用モデルを学習することを特徴とする。
その上本発明に係る光学分析装置は、光を検出して出力値を得る光検出器の非線形性を補正する光学分析装置であって、前記学習部により生成された補正用モデルを用いて、前記光検出器の出力値を線形補正する補正部と、前記補正部により補正された補正後の出力値を用いて測定対象成分の濃度を算出する演算部とを備えることを特徴とする。
このように構成した本発明によれば、光学分析装置に用いられる検出器の出力補正を、複数の光量を用いて、かつ、簡単な構成で行うことができる。
MCT光検出器の入射光強度と出力値との関係を示す特性図である。 セル温度25.9℃、164.1℃におけるMCT光検出器の出力値を示す図である。 本発明の一実施形態における光学分析装置の全体を示す模式図である。 同実施形態の学習部を示す機能ブロック図である。 変形実施形態の機械学習方法の一例を示す図である。
以下、本発明に係る光学分析装置の一実施形態について図面を参照して説明する。
<装置構成>
本実施形態の光学分析装置100は、いわゆるFTIRと称されるフーリエ変換型赤外分光分析器100であり、図3に示すように、光源1、干渉計2、測定セル3、光検出器4及び演算処理装置5を具備したものである。
光源1は、ブロードなスペクトルを有する光(多数の波数の光を含む連続光)を射出するものであり、例えばタングステン・ヨウ素ランプや、高輝度セラミック光源が用いられる。
干渉計2は、光検出器4が検出する前の光の強度を変動させる変動部であり、同図に示すように、1枚のハーフミラー(ビームスプリッタ)21、固定鏡22及び移動鏡23を具備した、いわゆるマイケルソン干渉計を利用したものである。この干渉計2に入射した光源1からの光は、ハーフミラー21によって反射光と透過光に分割される。一方の光は固定鏡22で反射され、もう一方は移動鏡23で反射されて、再びハーフミラー21に戻り、合成されて、この干渉計2から射出される。
測定セル3は、測定対象となるガス(以下、サンプルガスともいう。)が導入されるものであり、干渉計2から出た光が、測定セル3内のサンプルガスを透過して光検出器4に導かれるようにしてある。この測定セル3は、例えばアルミニウムなどの金属製の筒状をなすセル本体と、セル本体の両端壁に設けられた石英ガラス等の透光性材料からなる窓部材とから構成されている。
光検出器4は、測定セル3から出力される光を検出するものであり、入射光強度が所定範囲内においてその出力値と線形となり、その他の範囲で非線形となるような光検出器である。ここでは、いわゆるMCT光検出器4と称されるものである。
演算処理装置5は、バッファ、増幅器などを有したアナログ電気回路と、CPU、メモリ、DSPなどを有したデジタル電気回路と、それらの間に介在するA/Dコンバータ等を有したものである。この演算処理装置5は、前記メモリに格納した所定プログラムにしたがってCPUやその周辺機器が協働することにより、サンプルガスを透過した光のスペクトルを光検出器4の出力値から算出するとともに、この光スペクトルから各波数の光の吸光度を求めてサンプルガスを分析する分析部51としての機能を発揮する。
この分析部51は、光スペクトルを以下のようにして算出する。
移動鏡23を進退させ、サンプルガスを透過した光強度を移動鏡23の位置を横軸にとって観測すると、単波数の光の場合、干渉によって光強度はサインカーブを描く。一方、サンプルガスを透過した実際の光は多数の波数を有する光であり、前記サインカーブは波数毎に異なるから、実際の光強度は、各波数の描くサインカーブの重ね合わせとなり、干渉パターン(インターフェログラム)は波束の形となる。
分析部51は、移動鏡23の位置を例えば図示しないHeNeレーザなどの測距計(図示しない)によって求めるとともに、移動鏡23の各位置における光強度を光検出器4によって求め、これらから得られる干渉パターンを高速フーリエ変換(FFT)することによって、各波数成分を横軸とし、光強度信号を縦軸とした光スペクトルに変換する。
このとき、サンプルガスを透過した光、つまり光検出器4に入射した光の強度は、光検出器4の出力値に基づいて算出されるが、そのためには、光検出器4の出力値と入射光強度との関係(検量線)を予め求めて、メモリに記憶させておく必要がある。なお、この実施形態における検量線とは、光検出器4の出力値を入射光強度に変換するための演算式であるが、テーブルやマップでも構わない。
本実施形態では、測定セル3を多段階に温調できるように構成されており、光学分析装置100は、測定セル3を温調するための温調機構6と、当該温調機構6を制御する温調機構制御部52とを有している。これら温調機構6及び温調機構制御部52により温調部が構成される。
温調機構6は、測定セル3の周囲に設けられた例えば電熱線などのヒータである。また、温調機構制御部52は、ヒータ6に電力を供給する電源回路を制御することによって、ヒータ6による測定セル3の加熱温度を制御するものである。例えば、温調機構制御部52は、測定セル3の温度が25.9℃(低温時)と164.1℃(高温時)となるように制御する。
<機械学習機能について>
しかして本実施形態の演算処理装置5は、光検出器4の出力値を線形補正するための補正用モデルを機械学習する機能を有している。具体的に演算処理装置5は、同一の試料を測定した場合に、インターフェログラムにおける直流オフセット成分(DC成分)の大きさに関わらず、干渉計2による変動成分(AC成分)を等しくなるように補正する補正用モデルを機械学習する機能を有する。
具体的に演算処理装置5の分析部51は、図3に示すように、測定対象成分の濃度を算出する演算部51aの他に、光検出器4の出力値を線形補正する補正用モデルを学習する学習部51bと、学習部51bにより得られた補正用モデルを保存する保存部51cと、補正用モデルを用いて、光検出器4の出力値を線形補正する補正部51dとを有している。
学習部51bは、強化学習アルゴリズムを用いて補正用モデルを機械学習するものであり、光検出器4の補正前の出力値(例えばインターフェログラム)が入力されて、補正前の出力値を線形補正するための補正用モデルを演算する強化学習部53と、補正用モデルを、光検出器4の線形領域における出力値を用いて評価するモデル評価部54と、モデル評価部54の評価結果に応じて予め定められた報酬を算出する報酬算出部55とを備えている。
強化学習部53は、測定セル3を複数の温度に変化させたときのそれぞれの出力値を取得して、補正用モデルを学習するものである。
また、強化学習部53は、モデル評価部54の評価結果に基づいて報酬算出部55により得られた報酬を用いて、補正用モデルを更新する。ここで、補正用モデルを更新することには、補正用モデルに含まれる各種パラメータを変更することを含む。
モデル評価部54は、「補正用モデルを用いて補正された補正後の出力値」と「光検出器4の線形領域における出力値」との間の時間成分の一致度、周波数成分の一致度又は分光特性の一致度の少なくとも1つにより評価する。時間成分の一致度とは、例えば、時間を横軸とし、光強度信号を縦軸としたグラフ波形の一致度である。周波数成分の一致度とは、例えば、光強度信号をフーリエ変換して得られた周波数成分を用いたものであり、波数成分を横軸とし、光強度信号を縦軸としたグラフ波形の一致度である。分光特性の一致度とは、例えば、波長成分を横軸とし、光強度信号を縦軸としたグラフ波形の一致度である。
なお、光検出器4の線形領域における出力値は、光源1からの入射光強度が弱く、入射光強度と出力値との関係が線形となる所定範囲内のものを用いてもよいし、光源1と測定セル3との間にフィルタ(図示せず)を設置して、光源1から入射光強度を弱めて、入射光強度と出力値との関係が線形となるものを用いてもよい。
報酬算出部55は、モデル評価部54により得られた「補正用モデルを用いて補正された補正後の出力値」と「光検出器4の線形領域における出力値」との一致度に応じた報酬を算出するものである。ここで、一致度に応じた報酬は予め定められて記憶されている。
<強化学習方法の一例>
次に、強化学習における補正用モデルの算出方法の一例について説明する。
本実施形態のFTIRのように変動部として干渉計2を有する光学分析装置100の場合には、干渉計2を通して光源1の光を測定セル3に照射する。これにより、光検出器4の出力値(インターフェログラム)がDC成分及びAC成分の両方を含むことになる。
一方、干渉計2を有さない例えばNDIR(非分散形赤外分光分析器)等の分析装置の場合には、図5に示すように、光源1と測定セル3との間に変動部として光学チョッパ等の光変調器7を配置し、当該光変調器7により、光検出器4の出力値(インターフェログラム)がDC成分及びAC成分の両方を含むようにする。ここで、光学チョッパ7に形成された複数の光通過孔はそれぞれ開口径が異なり、光通過孔によって通過する光量が異なるように構成されていてもよい。なお、上記何れの場合であっても、測定セル3にサンプルガスや基準ガスを封入してもよいし、封入していなくても良い。
この状態で、温調機構制御部52により温調機構6の温度を例えば、25.9℃(低温時)から164.1℃(高温時)となるように制御する。この間の各温度において光検出器4は出力値を出力する。このように測定セル3の温度を変更させることによって、光検出器4の出力値に含まれるDC成分の大きさを変化させる。
各温度における補正前の出力値(測定セル3からの赤外輻射によってDC成分が変化している。)を強化学習部53に入力する。そうすると、強化学習部53が補正用モデルを算出して、モデル評価部54が補正用モデルを評価して、報酬算出部55が報酬を算出する。そして、強化学習部53は、報酬算出部55の報酬が最大化するように補正用モデルを更新する。
上記により得られた補正用モデルは、保存部51cに保存される。そして、補正部は、保存部51cに保存された補正用モデルを用いて、光検出器4の補正前の出力値を線形補正する。この線形補正された出力値が演算部51aに入力され、演算部51aは、線形補正された出力値を用いて測定対象成分の濃度を算出する。
<本実施形態の効果>
このように構成した本実施形態によれば、機械学習アルゴリズムにより光検出器4の出力値を線形補正するための補正用モデルを学習しているので、従来の補正式演算アルゴリズムを不要とし、光検出器4の出力補正を、複数の光量を用いて、かつ、簡単な構成で行うことができる。
<その他の変形実施形態>
なお、本発明は前記実施形態に限られるものではない。
例えば、強化学習部53への入力信号として、光検出器4の補正前の出力値の他に、入射光強度をゼロ(光検出器への入射光量をゼロ)としたときの光検出器4の出力値、光検出器4の出力値(インターフェログラム)から分離された直流オフセット成分の値、又は、測定セル3の温度の少なくとも1つを用いても良い。この場合、強化学習部53は、光検出器4の補正前の出力値とともに、それらの値を用いて、補正用モデルを学習する。
モデル評価部54は、光検出器4の線形領域における出力値を用いる他に、補正用モデルを評価するためのサンプルの既知濃度を用いて補正用モデルを評価するようにしても良い。例えば、補正用モデルを用いて補正された補正後の出力値により得られた算出濃度と評価用の正しい濃度(既知濃度)とを比較することによって補正用モデルを評価しても良い。
前記実施形態では、光学分析装置100が機械学習機能を有する構成であったが、光学分析装置100とは別体に機械学習装置を設けて、当該機械学習装置により補正用モデルを強化学習し、当該補正用モデルを光学分析装置100に記憶させても良い。この場合、機械学習装置には、学習部51b及び保存部51cが設けられる。一方、光学分析装置100には、保存部51c及び補正部51dが設けられることになる。
前記実施形態では、光検出器の出力値におけるAC成分を形成する変動部として干渉計2又は光学チョッパ等の光変調器7を用いていたが、干渉計2や光変調器7を用いること無く、変動部は、測定セル3に導入するサンプルの濃度を変化させて光検出器4が検出する前の光の強度を変動させることによって、光検出器4の出力値におけるAC成分を形成するものであっても良い。
前記実施形態の学習部51bは、強化学習アルゴリズムにより機械学習を行うものであったが、教師あり学習アルゴリズムにより機械学習を行うものであっても良い。
この場合、学習部51bは、教師データ(学習データセット)として、光検出器4の線形領域における出力値又は測定セル3に導入されたサンプルの既知濃度を用いて、補正用モデルを学習する。具体的には、教師データとして、補正前の出力値と、光検出器4の線形領域における出力値の時間成分、周波数成分又は分光特性の少なくとも1つとからなるデータセットが入力される。このデータセットにより、学習部51bは、補正用モデルを学習する。
光検出器4は、MCT光検出器に限られない。入射光強度が所定範囲内においてその出力値と線形となり、その他の範囲で非線形となるような光検出器に本発明を適用することができる。
また、本発明は、フーリエ変換型赤外分光分析器に限られず、他のタイプの分光分析器や分光器(前記実施形態で言えば光スペクトルを求めるまでの構成のもの)に適用して同様の効果を奏し得るものである。
その他、本発明は前記実施形態に限られること無く、その趣旨を逸脱しない範囲で種々の変形が可能である。
100・・・光学分析装置
1 ・・・光源
3 ・・・測定セル
4 ・・・光検出器
5 ・・・演算処理装置
51 ・・・分析部
52 ・・・温調機構制御部(温調部)
51a・・・演算部
51b・・・学習部
51c・・・保存部
51d・・・補正部
7 ・・・温調機構(温調部)

Claims (10)

  1. 光を検出して出力値を得る光検出器の非線形性を補正する光学分析装置であって、
    前記光を照射する光源と、
    前記光源から照射された光が入力される測定セルと、
    前記測定セルから出力される光を検出して出力値を得る前記光検出器と、
    前記光検出器が検出する前の光の強度を変動させる変動部と、
    前記光検出器の補正前の出力値となるインターフェログラムと、前記光検出器の線形領域における出力値となるインターフェロフラム又は前記測定セルに導入されたサンプルの既知濃度とを用いて、前記光検出器の出力値となるインターフェログラムに含まれる変動成分を補正して前記光検出器の出力値を線形補正するための補正用モデルを機械学習する学習部と、
    前記学習部により得られた補正用モデルを保存する保存部とを備える、光学分析装置。
  2. 前記光学分析装置は、前記測定セルを複数の温度に温調する温調部をさらに備え、
    前記測定セルから出力される光は、温調された前記測定セルから発生し、その温度によって光量が変化する光と、前記光源から照射されて前記測定セルを通過した光とを含むものであり、
    前記学習部は、前記測定セルを複数の温度に温調させたときのそれぞれにおいて前記光検出器により検出された複数の補正前の出力値となるインターフェログラムを用いて、前記補正用モデルを学習する、請求項記載の光学分析装置。
  3. 前記学習部は、さらに、前記光の強度をゼロとしたときの前記光検出器の出力値となるインターフェログラム、前記光検出器の出力値に含まれる前記測定セルの赤外輻射によるオフセット成分、又は、前記測定セルの温度の少なくとも1つを取得して、前記補正用モデルを学習する、請求項1又は2に記載の光学分析装置。
  4. 前記変動部は、前記光源及び前記測定セルの間に設けられ、前記光源から照射された光を反射させる固定ミラーと、前記光源から照射された光を往復移動しながら反射する移動ミラーとを備えた干渉計であるか、又は光を通過又は遮断させて前記光の強度を変動させる光変調器である、請求項乃至の何れか一項に記載の光学分析装置。
  5. 前記学習部は、
    前記光検出器の補正前の出力値となるインターフェログラムが入力されて、前記補正前の出力値となるインターフェログラムに含まれる変動成分を補正して前記光検出器の出力値を線形補正を線形補正するための補正用モデルを演算する強化学習部と、
    前記補正用モデルを、前記光検出器の線形領域における出力値となるインターフェログラム又は前記測定セルに導入されたサンプルの既知濃度を用いて評価するモデル評価部とを備え、
    前記強化学習部は、前記モデル評価部の評価結果に基づく報酬を用いて前記補正用モデルを更新する、請求項1乃至の何れか一項に記載の光学分析装置。
  6. 前記モデル評価部は、前記補正用モデルを用いて補正された補正後の出力値となるインターフェログラムと前記線形領域における出力値となるインターフェログラムとの間の時間成分の信号波形の一致度、周波数成分の信号波形の一致度又は分光特性の信号波形の一致度の少なくとも1つにより評価する、請求項記載の光学分析装置。
  7. 前記光学分析装置は、前記光源からの光をインターフェログラムに変調する干渉計と、前記光検出器の出力値となるインターフェログラムをフーリエ変換し、このフーリエ変換によって得られるスペクトルに基づいて測定対象成分の濃度を算出する演算部とをさらに備え、フーリエ変換型分光分析を行うものである、請求項1乃至の何れか一項に記載の光学分析装置。
  8. 光を検出して出力値を得る光検出器の非線形性を補正する光学分析装置に用いられる補正用モデルを生成する機械学習装置であって、
    前記光検出器の補正前の出力値となるインターフェログラムと、前記光検出器の線形領域における出力値となるインターフェログラム又は前記光学分析装置に導入されたサンプルの既知濃度とを用いて、前記光検出器の出力値となるインターフェログラムに含まれる変動成分を補正して前記光検出器の出力値を線形補正するための補正用モデルを学習する学習部と、
    前記学習部により得られた補正用モデルを保存する保存部とを備える、機械学習装置。
  9. 光を検出して出力値を得る光検出器の非線形性を補正する光検出器を備える光学分析装置に用いられる補正用モデルを生成する機械学習方法であって、
    前記光検出器の補正前の出力値となるインターフェログラムと、前記光検出器の線形領域における出力値となるインターフェログラム又は測定セルに導入されたサンプルの既知濃度とを用いて、前記光検出器の出力値となるインターフェログラムに含まれる変動成分を補正して前記光検出器の出力値を線形補正するための補正用モデルを学習する、機械学習方法。
  10. 光を検出して出力値を得る光検出器の非線形性を補正する光学分析装置であって、
    請求項記載の機械学習装置により生成された補正用モデルを用いて、前記光検出器の出力値を線形補正する補正部と、
    前記補正部により補正された補正後の出力値を用いて測定対象成分の濃度を算出する演算部とを備える光学分析装置。
JP2018234772A 2018-12-14 2018-12-14 光学分析装置、並びに光学分析装置に用いられる機械学習装置及びその方法 Active JP7014701B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2018234772A JP7014701B2 (ja) 2018-12-14 2018-12-14 光学分析装置、並びに光学分析装置に用いられる機械学習装置及びその方法
EP19194274.7A EP3667294A1 (en) 2018-12-14 2019-08-29 Optical analysis device, and machine learning device used for optical analysis device, and method for the same
US16/564,861 US20200191707A1 (en) 2018-12-14 2019-09-09 Optical analysis device, and machine learning device used for optical analysis device, and method for the same

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018234772A JP7014701B2 (ja) 2018-12-14 2018-12-14 光学分析装置、並びに光学分析装置に用いられる機械学習装置及びその方法

Publications (2)

Publication Number Publication Date
JP2020094990A JP2020094990A (ja) 2020-06-18
JP7014701B2 true JP7014701B2 (ja) 2022-02-01

Family

ID=67809273

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018234772A Active JP7014701B2 (ja) 2018-12-14 2018-12-14 光学分析装置、並びに光学分析装置に用いられる機械学習装置及びその方法

Country Status (3)

Country Link
US (1) US20200191707A1 (ja)
EP (1) EP3667294A1 (ja)
JP (1) JP7014701B2 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023127262A1 (ja) * 2021-12-28 2023-07-06 株式会社堀場製作所 機械学習装置、排ガス分析装置、機械学習方法、排ガス分析方法、機械学習プログラム、及び、排ガス分析プログラム

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017026506A (ja) 2015-07-24 2017-02-02 株式会社堀場製作所 分光分析器に用いられる光検出器の出力補正方法
JP2018191073A (ja) 2017-04-28 2018-11-29 東海光学株式会社 着色材料の配合量の算出方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4927269A (en) * 1989-01-31 1990-05-22 Bruke Analytische Messtechnik Gmbh Correction of non-linearities in detectors in fourier transform spectroscopy
JPH11304705A (ja) * 1998-04-23 1999-11-05 Nippon Soken Inc 赤外線吸収方式のガス濃度計測方法および装置
ITPD20130095A1 (it) * 2013-04-15 2014-10-16 Cnr Consiglio Naz Delle Ric Erche Apparato per la misura di concentrazione di gas in tracce mediante la spettroscopia scar
US9535053B1 (en) * 2016-02-04 2017-01-03 Nova Biomedical Corporation Analyte system and method for determining hemoglobin parameters in whole blood

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017026506A (ja) 2015-07-24 2017-02-02 株式会社堀場製作所 分光分析器に用いられる光検出器の出力補正方法
JP2018191073A (ja) 2017-04-28 2018-11-29 東海光学株式会社 着色材料の配合量の算出方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Evalution of calibration datra in capillary electrophresis using artificial neural networks to increase preecision of analysis,JOUNAL OF CHROMATOGRAPHY A,ELSEVIER,2002年12月06日,vol.979,no.1-2,59-67

Also Published As

Publication number Publication date
JP2020094990A (ja) 2020-06-18
US20200191707A1 (en) 2020-06-18
EP3667294A1 (en) 2020-06-17

Similar Documents

Publication Publication Date Title
JP5983779B2 (ja) ガス吸収分光装置及びガス吸収分光方法
JP6981817B2 (ja) 分光分析装置及び分光分析方法
JP6863831B2 (ja) 光検出器の出力補正用演算式の算出方法、及び光検出器の出力補正方法
EP3218695B1 (en) Target analyte detection and quantification in sample gases with complex background compositions
EP3332231B1 (en) Determination and correction of frequency registration deviations for quantitative spectroscopy
EP3133380B1 (en) Photodetector output correction method used for spectroscopic analyzer or spectroscope, spectroscopic analyzer or spectroscope using this method and program for spectroscopic analyzer or spectroscope instructing this method
EP3332230B1 (en) Reconstruction of frequency registration deviations for quantitative spectroscopy
RU2400715C2 (ru) Способ градуировки спектрометра
US10101208B2 (en) Apparatus and method for measuring the concentration of trace gases by SCAR spectroscopy
CN108226064B (zh) 分析装置、计算机可读存储介质和分析方法
JP7135608B2 (ja) ガス吸収分光装置、及びガス吸収分光方法
JP2017026506A5 (ja)
US20180149586A1 (en) Gas analysis apparatus and gas analysis method
CN109520941B (zh) 在线光谱测量仪器的响应函数校正方法
CN106053021B (zh) 一种分布反馈式激光器时频响应曲线的确定方法
RU2571185C2 (ru) Способ компенсации дрейфа амплитуды в спектрометре и спектрометр, осуществляющий указанный способ
JP7014701B2 (ja) 光学分析装置、並びに光学分析装置に用いられる機械学習装置及びその方法
JP2018169203A (ja) ガス吸収分光装置及びガス吸収分光法
US9970867B2 (en) Method of determining the concentration of a gas component and spectrometer therefor
JP7440866B2 (ja) レーザ式ガス分析装置
JP2024504351A (ja) ガスセンサの校正方法およびその校正を用いたガスの測定方法
US20060238768A1 (en) Fourier transform infrared spectrometer
CN120558906A (zh) 基于双谱线的免标定激光吸收光谱信号处理方法
CN120390870A (zh) 光谱分析装置、光谱分析方法、光谱分析程序、学习装置、学习方法和学习程序
WO2006115490A1 (en) Fourier transform infrared spectrometer

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20201216

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210930

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20211102

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211227

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220120

R150 Certificate of patent or registration of utility model

Ref document number: 7014701

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