JP6983244B2 - 分析装置、分析方法、分析装置用プログラム及び分析用学習装置 - Google Patents

分析装置、分析方法、分析装置用プログラム及び分析用学習装置 Download PDF

Info

Publication number
JP6983244B2
JP6983244B2 JP2019535138A JP2019535138A JP6983244B2 JP 6983244 B2 JP6983244 B2 JP 6983244B2 JP 2019535138 A JP2019535138 A JP 2019535138A JP 2019535138 A JP2019535138 A JP 2019535138A JP 6983244 B2 JP6983244 B2 JP 6983244B2
Authority
JP
Japan
Prior art keywords
correlation
total
reference sample
sample
data
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
JP2019535138A
Other languages
English (en)
Other versions
JPWO2019031331A1 (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
Publication of JPWO2019031331A1 publication Critical patent/JPWO2019031331A1/ja
Application granted granted Critical
Publication of JP6983244B2 publication Critical patent/JP6983244B2/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/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/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/129Using chemometrical methods
    • G01N2201/1293Using chemometrical methods resolving multicomponent spectra

Landscapes

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

Description

本発明は、測定試料に光を照射して得られるスペクトルデータに基づいて当該測定試料を分析する分析装置に関するものである。
従来、例えば自動車の排ガスに含まれるハイドロカーボン全体(THC)の濃度や量の測定には、FID分析装置や、FTIRなどの分光分析装置が用いられている。
しかしながら、FID分析装置は、分析精度において優れてはいるものの、助燃ガスとして水素ガス(H)やその濃度調整などに用いられるヘリウムガス(He)を供給しなければならないので、取り扱いの難しさやランニングコストの増大などの問題がある。
他方、FTIR分析装置は、より簡便に取り扱えてランニングコストも安価に済むという利点はあるものの、分析精度に難がある。すなわち、FTIR分析装置においては、光スペクトルから各ハイドロカーボン(HC)の濃度をそれぞれ個別に求め、次にそれ等を重みづけして足し合わせるという2段階の演算を行っており、各HCの濃度測定で生じ得る誤差に、重みづけ係数の設定において生じ得る誤差が重畳されるので、測定精度を向上させることが非常に難しい。
例えば、自動車の排ガスのTHC濃度は一般に100ppmC以下であるが、FTIR分光分析装置で得られる各HCの濃度精度は精々1〜10ppm程度であり、これらの情報から十分な精度で測定することは困難である。また、FTIR分光分析装置で濃度定量される各HC成分が、全種類のHC成分を網羅しているかどうかを保証できないため、これがTHCの測定精度を高められない要因ともなっている。
その他、例えば、特許文献1においては、個別測定を回避すべく、測定試料に含まれる一部の複数種のHCについては、一つの化合物として単一成分のスペクトルが生じるとみなし、そのスペクトルからその総量を一括して定量測定するようにしたものもある。
特開平4−265842号公報
そこで、本発明は、FTIR分光分析装置などの分析装置において、測定精度を飛躍的に向上させることをその主たる所期課題とするものである。
すなわち本発明に係る分析装置は、測定試料に光を照射して得られるスペクトルデータに基づいて当該測定試料を分析するものであって、所定複数成分の合計分析値が予め求められている参照試料の当該合計分析値に基づいて、測定試料のスペクトルデータから当該測定試料における前記複数成分の合計分析値を算出する合計分析値算出部を備えていることを特徴とする。
また、この合計分析値算出部は、測定試料に光を照射して得られるスペクトルデータに基づいて当該測定試料を分析するものであって、所定複数成分の合計分析値が予め求められている参照試料のスペクトルデータに基づいて、測定試料のスペクトルデータから当該測定試料における前記複数成分の合計分析値を算出するものでもよい。
なお、スペクトルデータとは、測定試料(又は参照試料)を透過、反射又は散乱した光のスペクトルデータの他、測定試料(又は参照試料)で吸収された光のスペクトルデータ(吸光スペクトルデータ)や、測定試料(又は参照試料)に含まれる干渉成分による影響を軽減又は取り除いた補正吸光スペクトルデータをも含むものである。
また、合計分析値とは、複数成分の濃度の合計値、質量の合計値など、複数各成分の物理量の合計値のことである。
このようなものであれば、合計分析値及びスペクトルデータが既知の参照試料に係るデータを基準として、測定試料に含まれる複数成分の合計分析値を算出することに加え、測定試料のスペクトルデータから直接的に合計分析値を算出することができるので、従来のように、スペクトルデータから個別成分の分析値を算出するといった途中演算での誤差を排除することができ、分析精度の飛躍的な向上が可能となる。
また、助燃ガス(H)などの取扱いが難しいガスを使用する必要はないので、簡便性やランニングコストの点にも寄与し得る。
より具体的には、前記合計分析値算出部が、参照試料のスペクトルデータと当該参照試料の合計分析値との相関を示すデータを格納している相関データ格納部と、測定試料のスペクトルデータに前記相関データ格納部の相関を適用して、当該測定試料における前記複数成分の合計分析値を算出する算出本体部とを備えたものを挙げることができる。
この場合、前記合計分析値算出部が、参照試料のスペクトルデータと当該参照試料の合計分析値とに基づいて、スペクトルデータと合計分析値との相関を算出し、前記相関データ格納部に格納するように構成しておけばよい。
分析精度の向上にさらに資するためには、参照試料の合計分析値を精度よく測定できる他の分析装置によって取得することが好ましい。
すなわち、他の方式の分析装置で測定された参照試料の合計分析値を示すデータを受け付ける受付部をさらに備え、前記相関算出部が、前記受付部で受け付けた合計分析値と、本分析装置によって得られた前記参照試料のスペクトルデータとに基づいて前記相関を算出するものであることが好ましい。
スペクトルデータから合計分析値への相関を直接求めるには、機械学習や深層学習を利用することが好ましい。
この場合、算出された相関は教師データが多いほど精度が上がる。したがって、前記相関算出部が、互いに異なる複数の参照試料における各合計分析値と各スペクトルデータとから前記相関関係を算出するように構成されたものが望ましい。
また、使用を経るにつれ相関精度が向上するようにするには、新たな参照試料に係る合計分析値及びそのスペクトルデータを取得した場合、前記相関算出部が、新たな参照試料に係る合計分析値及びそのスペクトルデータを加味して前記相関関係を更新するように構成しておけばよい。
相関精度を向上させる他の態様としては、前記相関算出部が、参照試料の温度及び/又は圧力をもパラメータとして前記相関を算出するものであり、前記算出本体部が、測定試料の温度及び/又は圧力をパラメータとして特定される前記相関に基づいて当該測定試料の合計分析値を算出するものであることが望ましい。
複数成分の分析精度をさらに向上させるためには、前記スペクトルデータが、前記複数成分に干渉する他成分の影響を軽減又は除去したものであることが好ましい。
前記複数成分から1以上の任意の成分を差し引いた成分の合計分析値を算出するには、測定試料に光を照射して得られるスペクトルデータに基づいて、当該測定試料に含まれる個々の単一成分を分析する主分析部をさらに備えていればよい。
例えばHCの場合、分子量の小さいHCについては濃度などの分析値を、従来の例えばFTIR分析装置でも精度よく求めることができるが、分子量の大きなHCになると、スペクトルデータの分離が困難になり、分析精度が悪化する。そのため、分子量の小さなHCはFTIR分析装置で求めておき、そこから分子量の小さなHCのスペクトルデータを求めて元のスペクトルデータから減算して補正し、補正したスペクトルから分子量の大きなHCの分析を行ってもよい。
本発明の効果が顕著に奏される具体的態様としては、前記測定試料又は参照試料が自動車の排ガスであり、分析対象となる成分がハイドロカーボン(HC)を挙げることができる。また、所定複数成分の合計分析値としてTHC濃度を挙げることができる。
また、本分析装置はFTIR方式のものが好ましく、THCを分析する場合には、他の方式の分析装置は、FID分析装置であることが好ましい。
参照試料だけを用いて相関のみを算出する機能に特化させた分析用学習装置も、本発明の一つである。
その場合は、所定複数成分の合計分析値が予め求められている参照試料に光を照射して得られるスペクトルデータを受け付ける受付部と、複数の互いに異なる参照試料の合計分析値及びスペクトルデータを格納している参照試料データ格納部と、前記参照試料データ格納部を参照し、各参照試料のスペクトルデータと合計分析値との間に共通する相関を算出する相関算出部とを備えているものが好適である。
このように構成した本発明によれば、合計分析値及びスペクトルデータが既知の参照試料に係るデータを利用して、測定試料を分析するうえ、従来のように、スペクトルデータから個別成分の分析値を算出するといった途中演算での誤差を排除することができるので、分析精度の飛躍的な向上が可能となる。また、助燃ガス(H)などの取扱いが難しいガスを使用する必要はないので、簡便性やランニングコストの点にも寄与し得る。
本発明の一実施形態における分析装置を含む排ガス測定システムの全体図である。 同実施形態における分析装置の全体を示す模式図である。 同実施形態における演算処理装置の機能ブロック図である。 同実施形態における分析装置の動作を示すフローチャートである。 他の実施形態における演算処理装置の機能ブロック図である。
100・・・分析装置
51・・・主分析部
52・・・合計分析値算出部
521・・・相関算出部
522・・・算出本体部
53・・・受付部
以下に本発明の一実施形態に係る分析装置100について図面を参照して説明する。
本実施形態の分析装置100は、排ガス測定システム200の一部を構成するものである。この排ガス測定システム200は、図1に示すように、シャシダイナモ300と、FID分析装置400と、本分析装置100とを具備している。
しかして、本分析装置100は、図2に示すように、光源1、干渉計(分光部)2、試料セル3、光検出器4、演算処理装置5等を具備した、いわゆるFTIRと称されるフーリエ変換型赤外分光分析装置である。この分析装置100(以下、区別のため、FTIR分析装置100ということもある。)は、測定試料である排ガスのTHC濃度(又は量)を測定する排ガス分析装置として使用される。
光源1は、ブロードなスペクトルを有する光(多数の波数の光を含む連続光)を射出するものであり、例えばタングステン・ヨウ素ランプや、高輝度セラミック光源が用いられる。
干渉計2は、同図に示すように、1枚のハーフミラー(ビームスプリッタ)21、固定鏡22及び移動鏡23を具備した、いわゆるマイケルソン干渉計を利用したものである。この干渉計2に入射した前記光源1からの光は、前記ハーフミラー21によって反射光と透過光に分割される。一方の光は固定鏡22で反射され、もう一方は移動鏡23で反射されて、再びハーフミラー21に戻り、合成されて、この干渉計2から射出される。
試料セル3は、測定試料である排ガスが導入される透明セルであり、前記干渉計2から出た光が、該試料セル3内の測定試料を透過して前記光検出器4に導かれるようにしてある。
光検出器4は、ここでは、いわゆるMCT光検出器4と称されるものである。
演算処理装置5は、バッファ、増幅器などを有したアナログ電気回路と、CPU、メモリ、DSPなどを有したデジタル電気回路と、それらの間に介在するA/Dコンバータ等を有したものである。
該演算処理装置5は、前記メモリに格納した所定プログラムにしたがってCPUやその周辺機器が協動することにより、図3に示すように、試料を透過した光のスペクトルを示す透過光スペクトルデータを前記光検出器4の出力値から算出し、該透過光スペクトルデータから吸光スペクトルデータを算出して測定試料に含まれる種々の成分を特定し、かつそれぞれの成分の濃度(又は量)算出する主分析部51としての機能を発揮する。
この主分析部51は、スペクトルデータ生成部511と、個別成分分析部512とを具備している。
移動鏡23を進退させ、試料を透過した光強度を移動鏡23の位置を横軸にとって観測すると、単波数の光の場合、干渉によって光強度はサインカーブを描く。一方、試料を透過した実際の光は連続光であるから、前記サインカーブは波数毎に異なるから、実際の光強度は、各波数の描くサインカーブの重ね合わせとなり、干渉パターン(インターフェログラム)は波束の形となる。
前記スペクトルデータ生成部511は、移動鏡23の位置を例えば図示しないHeNeレーザなどの測距計(図示しない)によって求めるとともに、移動鏡23の各位置における光強度を光検出器4によって求め、これらから得られる干渉パターンを高速フーリエ変換(FFT)することによって、各波数成分を横軸とした透過光スペクトルデータに変換する。そして、例えば試料セルが空の状態で予め測定しておいた透過光スペクトルデータに基づいて、測定試料の透過光スペクトルデータを吸光スペクトルデータにさらに変換する。
前記個別成分分析部512は、例えば吸光スペクトルデータの各ピーク位置(波数)及びその高さから測定試料に含まれる種々の成分を特定し、かつそれぞれの成分の濃度(又は量)を算出する。
しかして、この実施形態では、図3に示すように、測定試料である排ガスのTHC濃度(又は量)を精度よく測定できるように、前記演算処理装置5に、受付部53、合計分析値算出部52等としての機能をさらに付与してある。
前記受付部53は、FID分析装置400で測定された排ガスのTHC濃度を受け付けるものである。このFID分析装置400で測定されTHC濃度が既知となった排ガスを、以下では参照試料という。
ところで、この参照試料は、FID分析装置400とともに本FTIR分析装置100にも導入され、その吸光スペクトルデータが前記主分析部51によっても取得されるので、前記受付部53は、主分析部51で算出される中間情報である当該参照試料の吸光スペクトルデータをも受け付け、これをFID分析装置400によって測定された参照試料のTHC濃度と紐づけて参照試料データとし、これをメモリの所定領域に設定された参照試料データ格納部D1に格納する。
さらに、この実施形態においては、本システムに設けられた図示しないセンサやオペレータによる入力により、参照試料の温度及び圧力を少なくとも含む周辺状況データを取得できるようになっていて、前記受付部53は、当該参照試料の周辺状況データを取得して前記参照試料データに付帯させ、参照試料データ格納部D1に格納するように構成してある。
前記合計分析値算出部52は、前記参照試料データを教師データとして、測定試料(排ガス)の吸光スペクトルデータから当該測定試料におけるTHCの濃度を算出するものであり、より具体的には、相関算出部521と算出本体部522とを備えている。なお、このTHCが請求項でいう複数成分にあたり、THC濃度が請求項でいう合計分析値にあたる。
前記相関算出部521は、参照試料データ格納部D1に格納されている複数の参照試料データを参照し、それら参照試料データに共通する、吸光スペクトルデータとTHC濃度との相関を機械学習(又は人工知能、ディープラーニング等)によって算出するものである。算出された相関を示す相関データは、メモリの所定領域に設定された相関データ格納部D2に格納される。
なお、この相関算出部521においては、参照試料データが追加されるたびに学習を繰り返し、相関を更新するので、参照試料データが多いほど、その相関精度は向上していく。
また、この実施形態での相関算出部521は、参照試料の周辺状況データをもパラメータとして前記相関を算出する、すなわち、相関が参照試料の温度や圧力などに応じても変わるように構成してあるが、相関の算出に当たって周辺状況データを参照しなくともよい。
算出本体部522は、前記相関算出部521が算出した相関を測定試料のスペクトルデータに当てはめて測定試料のTHC濃度を算出するものである。このとき、受付部53においては測定試料の周辺状況データを取得するようにしてあるので、当該算出本体部522は、THC濃度の算出にあたり、測定試料の周辺状況データに対応する相関を適用するようにしている。
次にかかる構成の排ガス測定システム200の動作を図4を参照して説明する。
まず学習させる。そのために、シャシダイナモ300上で自動車を運転し、その排ガスである参照試料をFID分析装置400及びFTIR分析装置100に導く。なお、シャシダイナモ300上での自動車の運転に限らず、エンジンダイナモに接続されたエンジンを運転してもよいし、駆動系ダイナモにトランスミッション等の駆動系部品を運転してもよい。
すると、FID分析装置400ではTHC濃度が測定され、一方、FTIR分析装置100では、前記主分析部51によって当該排ガスの吸光スペクトルデータが測定される。
この実施形態では、FID分析装置400及びFTIR分析装置100は、例えば、一定タイミング(数msec〜数sec)ごとに同期して排ガス測定(分析)を行う。
各測定の都度、前記受付部53は、主分析部51で算出された排ガスの吸光スペクトルデータ及びFID分析装置400で分析された当該排ガスのTHC濃度を取得し、参照試料データとして参照試料データ格納部D1に次々格納する(ステップS1)。このとき、前記受付部53は、当該排ガスの温度及び圧力を取得して前記参照試料データに付帯させ、参照試料データ格納部D1に格納する。
しかして、自動車のエンジン状態は、運転開始からの時間経過やエンジン回転数の変動にともなって種々変化し、それに伴って排ガスの状態(成分や圧力、温度)も逐次変化していくため、上述した逐次測定によって、少なくともTHC濃度が互いに異なる複数の参照試料のデータが得られる。
次に、前記相関算出部521が、参照試料データ格納部D1に格納されている多くの参照試料データを参照し、それら参照試料データに共通する、吸光スペクトルデータとTHC濃度との相関を機械学習によって算出し(ステップS2)、その結果得られた相関を前記相関データ格納部D2に格納する(ステップS3)。
これで学習が終了する。
このように学習が終了した後、FTIR分析装置100だけを用いた実際のTHC濃度測定を行うことができることとなる。このTHC濃度測定にあたっては、試験対象となる他の自動車をシャシダイナモ300に設置して運転し、その排ガスをFTIR分析装置100に導く。FID分析装置400には、排ガス導入管の図示しないバルブを閉じるなどして排ガスを導入しないようにするとともに、FID分析装置400の動作も停止させておく。
FTIR分析装置100では、主分析部51が、排ガスの吸光スペクトルデータを取得する(ステップS4)。すると前記合計分析値算出部52(算出本体部522)が、前記相関データ格納部D2に格納されている相関を吸光スペクトルデータに適用してTHC濃度を算出する(ステップS5)。
しかして、このように構成した本分析装置100によれば、参照試料のスペクトルデータとTHC濃度の相関を予め求め、その相関に基づいて、測定試料のスペクトルデータから直接的にTHC濃度を算出するようにしており、従来のように、スペクトルデータから個別成分の分析値を算出するといった途中演算を介在させることがないので、途中演算で発生し得る誤差や情報の脱落を排除でき、分析精度を飛躍的に向上させることができる。
また、相関さえ予め求めておけば、THC濃度測定においてFID分析装置を用いる必要はないので、助燃ガス(H)や、その濃度調整に用いられる不活性ガス(He)などの取扱いが難しいガスを使用する必要がなくなり、簡便性やランニングコストの点にも寄与し得る。
なお、本発明は前記実施形態に限られるものではない。
例えば、THCの吸光度は、測定対象に水分やその他の干渉成分が含まれているから、それらの影響を低減ないし排除した補正吸光スペクトルに基づいてTHC濃度を算出するようにしてもよい。このようにすれば、分析精度がさらに向上する。
また他の実施形態では、吸光スペクトルデータとTHC濃度との相関の算出にあたって、参照試料データが示すTHC濃度に応じて吸光スペクトルデータとTHC濃度との相関を算出してもよい。
具体的には、相関算出部521は、THC濃度を複数(例えば2以上)の濃度区分に分けて、この濃度区分毎に吸光スペクトルデータとTHC濃度の相関を算出し、その相関データを相関データ格納部D2に格納してよい。ここで算出本体部522は、測定試料の吸光スペクトルデータを受け付けると、当該吸光スペクトルデータの面積等に基づいて、相関データ格納部D2に格納された複数の相関データから当該吸光スペクトルデータに適した1つの相関データを選択し、当該選択した相関データを測定試料のスペクトルデータに当てはめて、測定試料のTHC濃度を算出するようにしてもよい。
このようなものであれば、相関算出部521が濃度区分毎に吸光スペクトルデータとTHC濃度の相関を算出するので、より精度の高い相関データを得ることができる。そして、算出本体部522が、このような精度の高い相関データを用いて測定試料のTHC濃度を算出するので、測定精度をより向上することができる。
なお、THC濃度はより多くの濃度区分に分けられることが好ましい。このようなものであれば、より一層精度の高い相関データを得ることができる。
また他の実施形態では、吸光スペクトルデータとTHC濃度との相関の算出にあたって、参照試料に含まれる水の濃度をパラメータとしてもよい。
具体的には、受付部53は、個別成分分析部512が算出した参照試料に含まれる水(HO)の濃度を受け付け、前記参照試料データに付帯させて試料データ格納部D1に格納するように構成されていてもよい。相関算出部521は、参照試料データ格納部D1に格納されている複数の参照試料データとこれに付帯する水の濃度を参照し、当該水濃度をパラメータとして、吸光スペクトルデータとTHC濃度との相関を算出するようにしてもよい。算出本体部522は、THC濃度の算出にあたり、測定試料の水濃度に対応する相関を適用するようにしてよい。
FTIRによるTHC濃度の測定においては、排ガスに含まれる水の干渉影響及び共存影響が問題となることがあるが、このようなものであれば、水の濃度をパラメータとする相関を用いること、測定試料のTHC濃度の算出における水の影響を低減でき、測定精度をより向上することができる。
また、前記実施形態では、スペクトルデータをTHC濃度に直接的に結び付ける、機械学習により算出された相関を利用していたが、機械学習によらず、スペクトルデータとTHC濃度との関係を表す演算式を算出するようにしてもよい。
その場合は、例えば、スペクトルデータの値とTHC濃度とを直接的に結び付けるものではなく、主分析部51によってスペクトルデータから求められた個別HC成分それぞれの分析値(ここでは濃度)とTHC濃度との関係を表す演算式(相関の一種でもある。)を用いればよい。
前記演算式の具体例は下記のとおりである。
THC(R)=a・HC1(R)+a・HC2(R)+・・・(1)
ここで、THC(R)は、FID分析装置400によって測定された参照試料のTHC濃度、HC1(R)、HC2(R)、・・・は、主分析部51によって算出された参照試料における各ハイドロカーボンの濃度、a、a、・・・は重み付け係数である。
この場合は、複数の参照試料を分析して多変量解析などを施し、上記式(1)を満たす最適な重み付け係数a、a、・・・を定める。これが相関算出部521の機能である。
次に、測定試料を分析して算出された各ハイドロカーボンの濃度HC1(R)、HC2(R)、・・・を下記式(2)に代入して、測定試料のTHC濃度を算出する。これが、合計分析値算出部52の機能である。
THC(M)=a・HC1(M)+a・HC2(M)+・・・(2)
ここで、THC(M)は、合計分析値算出部52によって算出された測定試料のTHC濃度、HC1(M)、HC2(M)、・・・は、主分析部51によって算出された測定試料における各ハイドロカーボンの濃度である。
なお、上記例では、相関を機械学習にはよらず、多変量解析を用いて求めている。
すなわち、相関の算出には、機械学習のみならず、上述したような理論的に求められる数式を用いてもよい。
また、機械学習と理論的数式を組み合わせることも可能である。例えば、質量の小さなHCに関しては主分析部を用いて(つまり理論式によって)個別の濃度を求めておき、一方、質量の大きなHCについては、前記相関算出部及び算出本体部によって(つまり機械学習によって)トータル濃度を算出し、それらを合算してTHC濃度を求めるようにしてもよい。
さらに、相関の算出にあたって、参照試料の温度、圧力などのように、試料の物理的属性に係る値のみならず、他の周辺状況データをパラメータとしてもよい。例えば、エンジンの燃焼情報(過給、EGR、リッチ/ストイキ/リーン、層流、均一流、直噴、ポート噴射などに係る情報)、エンジンのヘッド形状、点火時期、触媒の構成、燃料の酸素量、無機ガス成分、Soot濃度、SOF濃度、エンジン型式、エンジン回転数、負荷情報、ホットスタート、コールドスタート、酸素濃度、触媒温度、ギア比等をパラメータとして加味してもよい。
これとは逆に、前記周辺状況データの全部または一部を相関算出のパラメータとはせず、分析装置で算出測定されたTHC濃度(合計分析値)に影響を強く及ぼす(つまり関連の高い)周辺状況データを抽出するような構成にしても構わない。
このようにすれば、THC濃度と関連度の高い周辺状況データ、すなわち設計パラメータを把握することができるので、自動車メーカーや触媒メーカーに対し、設計開発支援システムとして提供することが可能になる。
相関に用いる周囲状況データによっては、例えば、エンジンの型式が異なった場合には、都度、再学習をしなければならなくなる場合もある。その点でいえば、加味する周辺状況データを少なくし、スペクトルデータとTHC濃度(合計分析値)との相関を取ることができれば、より汎用性は高くなる。具体的には、相関算出には、周辺状況データを一切加味しなくともよいし、試料自身の物理状態(例えば、圧力、温度、屈折率、粘度など)については相関算出に加味し、その他の外部属性(例えば、エンジンの型式、点火時期など)については相関算出には加味しないといった態様も考えられる。
相関の算出にあたって周辺状況として参照試料の温度をパラメータとする場合、シャシダイナモ300上の自動車のテールパイプの出口温度をセンサ等により測定し、当該測定した温度をパラメータにすることが好ましい。
学習・分析に用いるスペクトルデータの範囲は、分析対象成分が含まれる波数範囲のみでもよいし、それを超えた所定範囲まで拡げてもよい。また干渉成分の波数範囲を除くようにしてもよい。
具体的には、学習・分析に用いるスペクトルデータの範囲を2800cm−1以上3200cm−1以下にしてもよい。
スペクトルデータの範囲がこのようなものであれば、分析対象成分であるHCの波数範囲を含むとともに、干渉成分である水の波数範囲(約3400cm−1以上)が除かれるので、測定試料のTHC濃度の算出における水の影響を低減でき、測定精度をより向上することができる。
分子量の小さいHCについては各濃度を主分析部51で求めておき、分子量の大きな1又は複数のHCの合計濃度を、THC濃度から差し引いて算出してもよい。
さらに、本分析装置から個別成分分析部の機能を排除するとともに参照試料のみを導入するようにして、相関算出だけを行う分析用学習装置を構成してもよい。この分析用学習装置によって求められた相関は、他のFTIR分析装置で利用することができる。
前記実施形態の分析装置100は、参照試料のスペクトルデータとTHC濃度の相関を自ら算出するものであったがこれに限らない。他の実施形態の分析装置100は、相関算出だけを行う他の分析用学習装置により予め算出された相関を用い、その相関に基づいて測定試料のスペクトルデータから直接的にTHC濃度を算出するようにしてもよい。
具体的には分析装置100は、図5に示すように、演算処理装置5が参照試料データ格納部D1及び相関算出部521としての機能を有さないように構成されてもよい。ここでは受付部53が、他の分析用学習装置が予め算出した相関を示す相関データ(すなわち、学習済みデータ)を、ネットワーク等を介して受付けるとともに相関データ格納部D2に予め格納しておいてよい。そして算出本体部522が、相関データ格納部D2に予め格納されている相関データを測定試料の吸光スペクトルデータに当てはめて測定試料のTHC濃度を算出するようにしてよい。
なお受付部53は、所定の期間毎に他の分析用学習装置から新しい相関データを受け付けて、相関データ格納部D2に格納されている相関データを定期的に更新するようにしてよい。
参照試料及び測定対象試料は、前記実施形態では排ガスであったが、大気やその他のガス、あるいは液体でも構わない。
また、参照試料は、測定試料と同種のものである必要はなく、例えば窒素などの主成分に分析対象となる複数成分を混合して生成した標準ガス等を用いてもよい。この場合は、複数成分の合計分析値が既知なので、参照試料の複数成分を分析するための他の分析装置を用いる必要はない。
さらに、分析対象となる成分は、HCに限られず、窒素酸化物やダイオキシンなど、その他の成分でもよい。
分析装置としては、光を測定試料に照射してそのスペクトルから分析をするものであればよく、NDIRにも適用可能であるし、分光分析装置以外の、例えば散乱型粒子径分布測定装置などにも適用可能である。また、本発明は、自動車の排ガスの分析に限らず、船舶、航空機、農業用機械、工作機械等の内燃機関の排ガスを分析することも可能である。
その他、本発明は前記実施形態に限られず、その趣旨を逸脱しない範囲で種々の変形が可能であるのは言うまでもない。
本発明によれば、合計分析値及びスペクトルデータが既知の参照試料に係るデータを利用して、測定試料を分析するうえ、従来のように、スペクトルデータから個別成分の分析値を算出するといった途中演算での誤差を排除することができるので、分析精度の飛躍的な向上が可能となる。また、助燃ガス(H)などの取扱いが難しいガスを使用する必要はないので、簡便性やランニングコストの点にも寄与し得る。

Claims (21)

  1. ガスである測定試料に光を照射して得られるスペクトルデータに基づいて当該測定試料を分析するものであって、
    所定複数成分の合計分析値が予め求められているガスである参照試料の当該合計分析値に基づいて、測定試料のスペクトルデータから当該測定試料における前記複数成分の合計分析値を算出する合計分析値算出部を備え
    前記合計分析値算出部が、
    参照試料に光を照射して得られたスペクトルデータと当該参照試料の合計分析値との相関であって、機械学習により算出された相関を示すデータを格納している相関データ格納部と、
    測定試料のスペクトルデータに前記相関データ格納部に格納されている相関を適用して、当該測定試料における前記複数成分の合計分析値を算出する算出本体部とを備えたものである分析装置
  2. 前記合計分析値算出部が、
    参照試料のスペクトルデータと当該参照試料の合計分析値とに基づいて、スペクトルデータと合計分析値との相関を算出し、前記相関データ格納部に格納する相関算出部をさらに備えたものである請求項記載の分析装置。
  3. 他の方式の分析装置で測定された参照試料の合計分析値を示すデータを受け付ける受付部をさらに備え、
    前記相関算出部が、前記受付部で受け付けた合計分析値と、本分析装置によって得られた前記参照試料のスペクトルデータとに基づいて前記相関を算出するものである請求項記載の分析装置。
  4. 前記相関算出部が、複数の参照試料における各合計分析値と各スペクトルデータとから前記相関を算出するものであって、新たな参照試料に係る合計分析値及びそのスペクトルデータを取得した場合に、新たな参照試料に係る合計分析値及びそのスペクトルデータを加味して前記相関を更新するものである請求項2又は3記載の分析装置。
  5. 前記相関算出部が、参照試料の温度及び/又は圧力をもパラメータとして前記相関を算出するものであり、
    前記算出本体部が、測定試料の温度及び/又は圧力をパラメータとして特定される前記相関に基づいて当該測定試料の合計分析値を算出するものである請求項2乃至4いずれか記載の分析装置。
  6. 前記相関算出部が、参照試料の温度及び/又は圧力に加えて、エンジン又は触媒に関連する周辺状況データをもパラメータとして相関を算出し、
    前記算出本体部が、測定試料の温度及び/又は圧力に加えて、エンジン又は触媒に関連する周辺状況データをパラメータとして特定される前記相関に基づいて当該測定試料の合計分析値を算出する請求項記載の分析装置。
  7. 前記スペクトルデータが、前記複数成分に干渉する他成分の影響を軽減又は除去したものである請求項1乃至6のいずれか記載の分析装置。
  8. 測定試料に光を照射して得られるスペクトルデータに基づいて、当該測定試料に含まれる個々の単一成分を分析する主分析部をさらに備えている請求項1乃至7のいずれか記載の分析装置。
  9. 前記測定試料又は参照試料が自動車の排ガスであり、分析対象となる成分がハイドロカーボンである請求項1乃至8のいずれか記載の分析装置。
  10. 前記所定複数成分の合計分析値が、前記排ガスに含まれるハイドロカーボン全体の濃度である請求項記載の分析装置。
  11. FTIR方式のものである請求項1乃至10いずれか記載の分析装置。
  12. 他の方式の分析装置がFID分析装置である請求項11記載の分析装置。
  13. ガスである測定試料に光を照射して得られるスペクトルデータに基づいて当該測定試料を分析する方法であって、
    所定複数成分の合計分析値が予め求められているガスである参照試料の当該合計分析値に基づいて、測定試料のスペクトルデータから当該測定試料における前記複数成分の合計分析値を算出する合計分析値算出ステップを備え、
    当該分析値算出ステップにおいて、参照試料に光を照射して得られたスペクトルデータと当該参照試料の合計分析値との相関であって、機械学習により算出された相関に、測定試料のスペクトルデータを適用して、当該測定試料における前記複数成分の合計分析値を算出することを特徴とする分析方法。
  14. ガスである測定試料に光を照射して得られるスペクトルデータに基づいて当該測定試料を分析する分析装置に搭載されるプログラムであって、
    所定複数成分の合計分析値が予め求められているガスである参照試料の当該合計分析値に基づいて、測定試料のスペクトルデータから当該測定試料における前記複数成分の合計分析値を算出する合計分析値算出部としての機能を分析装置に発揮させるものであり、
    前記分析値算出部が、参照試料に光を照射して得られたスペクトルデータと当該参照試料の合計分析値との相関であって、機械学習により算出された相関に、測定試料のスペクトルデータを適用して、当該測定試料における前記複数成分の合計分析値を算出することを特徴とする分析装置用プログラム。
  15. 所定複数成分の合計分析値が予め求められているガスである参照試料に光を照射して得られるスペクトルデータを受け付ける受付部と、
    複数の互いに異なる参照試料の合計分析値及びスペクトルデータを格納している参照試料データ格納部と、
    前記参照試料データ格納部を参照し、各参照試料のスペクトルデータと合計分析値との間に共通する相関を機械学習により算出する相関算出部とを備えていることを特徴とする分析用学習装置。
  16. 前記受付部が、互いに異なる方式の分析装置により得られたものである、前記スペクトルデータとこれに対応する前記参照試料の前記合計分析値とを受け付け、これらを前記参照試料データ格納部に格納するように構成され、
    前記相関算出部が、当該異なる方式の分析装置によって得られた前記スペクトルデータと前記合計分析値とに基づいて前記相関を算出するものである請求項15に記載の分析用学習装置。
  17. 前記相関算出部が、複数の参照試料における各合計分析値と各スペクトルデータとから前記相関を算出するものであって、新たな参照試料に係る合計分析値及びそのスペクトルデータを取得した場合に、新たな参照試料に係る合計分析値及びそのスペクトルデータを加味して前記相関を更新するものである請求項15又は16に記載の分析用学習装置。
  18. 前記相関算出部が、前記参照試料の温度及び/又は圧力をもパラメータとして前記相関を算出するものである請求項15〜17のいずれか一項に記載の分析用学習装置。
  19. ガスである参照試料に光を照射して得られるスペクトルデータと、当該参照試料に含まれる所定複数成分の合計分析値との相関を学習する学習方法であって、
    所定複数成分の合計分析値が予め求められているガスである参照試料に光を照射して得られるスペクトルデータを受け付け、これを参照試料データ格納部に格納し、
    前記参照試料データ格納部に格納されている複数の互いに異なる参照試料の合計分析値及びスペクトルデータを参照して、各参照試料のスペクトルデータと合計分析値との間に共通する相関を機械学習により算出する学習方法。
  20. 所定複数成分の合計分析値が予め求められているガスである参照試料に光を照射して得られるスペクトルデータを受け付ける受付部と、
    複数の互いに異なる参照試料の合計分析値及びスペクトルデータを格納している参照試料データ格納部と、
    前記参照試料データ格納部を参照し、各参照試料のスペクトルデータと合計分析値との間に共通する相関を機械学習により算出する相関算出部としての機能をコンピュータに発揮させることを特徴とする分析用学習装置用プログラム。
  21. 測定試料に光を照射して得られるスペクトルデータに基づいて当該測定試料を分析するものであって、
    所定複数成分の合計分析値が予め求められている参照試料の当該合計分析値に基づいて、測定試料のスペクトルデータから当該測定試料における前記複数成分の合計分析値を算出する合計分析値算出部を備え、
    前記合計分析値算出部が、
    参照試料のスペクトルデータと当該参照試料の合計分析値との相関を示すデータを格納している相関データ格納部と、
    測定試料のスペクトルデータに前記相関データ格納部の相関を適用して、当該測定試料における前記複数成分の合計分析値を算出する算出本体部と、
    参照試料のスペクトルデータと当該参照試料の合計分析値とに基づいて、スペクトルデータと合計分析値との相関を算出し、前記相関データ格納部に格納する相関算出部を備え、
    前記相関算出部が、複数の参照試料における各合計分析値と各スペクトルデータとから前記相関を算出するものであって、新たな参照試料に係る合計分析値及びそのスペクトルデータを取得した場合に、新たな参照試料に係る合計分析値及びそのスペクトルデータを加味して前記相関を更新するものである分析装置。
JP2019535138A 2017-08-07 2018-08-01 分析装置、分析方法、分析装置用プログラム及び分析用学習装置 Active JP6983244B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2017152787 2017-08-07
JP2017152787 2017-08-07
PCT/JP2018/028812 WO2019031331A1 (ja) 2017-08-07 2018-08-01 分析装置、分析方法、分析装置用プログラム及び分析用学習装置

Publications (2)

Publication Number Publication Date
JPWO2019031331A1 JPWO2019031331A1 (ja) 2020-06-25
JP6983244B2 true JP6983244B2 (ja) 2021-12-17

Family

ID=65271392

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019535138A Active JP6983244B2 (ja) 2017-08-07 2018-08-01 分析装置、分析方法、分析装置用プログラム及び分析用学習装置

Country Status (4)

Country Link
US (1) US11448588B2 (ja)
EP (1) EP3667295A4 (ja)
JP (1) JP6983244B2 (ja)
WO (1) WO2019031331A1 (ja)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11081220B2 (en) 2018-02-02 2021-08-03 Carefusion 303, Inc. System and method for dispensing medication
US11222721B2 (en) 2018-05-04 2022-01-11 Carefusion 303, Inc. Peer community based anomalous behavior detection
CA3125744A1 (en) 2019-01-07 2020-07-16 Carefusion 303, Inc. Machine learning based safety controller
US11984212B2 (en) 2019-01-10 2024-05-14 Carefusion 303, Inc. System for monitoring dose pattern and patient response
EP3912166A1 (en) 2019-01-18 2021-11-24 CareFusion 303, Inc. Medication tracking system
JP6964114B2 (ja) * 2019-06-25 2021-11-10 株式会社堀場製作所 分析装置、分析方法及び分析装置用プログラム
WO2021005900A1 (ja) 2019-07-05 2021-01-14 株式会社堀場製作所 試料ガス分析装置、試料ガス分析方法及び試料ガス分析用プログラム
US20230296501A1 (en) 2020-07-13 2023-09-21 Horiba, Ltd. Analysis device, analysis method, program for analysis device, learning device for analysis, learning method for analysis, and program for learning device for analysis
WO2023127262A1 (ja) * 2021-12-28 2023-07-06 株式会社堀場製作所 機械学習装置、排ガス分析装置、機械学習方法、排ガス分析方法、機械学習プログラム、及び、排ガス分析プログラム

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0265842A (ja) 1988-08-31 1990-03-06 Shimadzu Corp Mriマルチエコー撮像法
MY107650A (en) * 1990-10-12 1996-05-30 Exxon Res & Engineering Company Method of estimating property and / or composition data of a test sample
JP3099898B2 (ja) 1991-02-21 2000-10-16 株式会社堀場製作所 Ftirを用いた多成分定量分析方法
JP3261264B2 (ja) * 1994-07-13 2002-02-25 株式会社堀場製作所 多成分水溶液の分析方法およびその分析装置
JPH09101259A (ja) 1995-10-08 1997-04-15 Horiba Ltd 分光光度計を用いた複数成分定量方法
JP3780646B2 (ja) 1997-07-18 2006-05-31 富士電機システムズ株式会社 下水中の複数成分定量方法とその装置
JP4079404B2 (ja) * 1999-06-04 2008-04-23 株式会社堀場製作所 Ftir法による多成分ガス分析方法
JP3754581B2 (ja) * 1999-08-31 2006-03-15 株式会社堀場製作所 多成分有機溶液の分析方法
JP2003014627A (ja) 2001-06-28 2003-01-15 Ishikawajima Harima Heavy Ind Co Ltd 多成分分析計の検量線作成方法
JP3899004B2 (ja) 2002-09-27 2007-03-28 株式会社堀場製作所 車載型hc測定装置
US6865472B2 (en) * 2002-09-27 2005-03-08 Horiba Ltd. Vehicle-installed exhaust gas analyzing apparatus
JP4505598B2 (ja) * 2004-02-09 2010-07-21 独立行政法人農業・食品産業技術総合研究機構 茶葉中に含まれる化学成分の定量方法
US9057718B2 (en) * 2005-10-06 2015-06-16 Jp3 Measurement, Llc Optical determination and reporting of hydrocarbon properties
EP2032966B1 (fr) 2006-06-07 2012-12-19 Valorex Methode de determination de la qualite nutritionnelle des lipides du lait
JP2010210360A (ja) 2009-03-10 2010-09-24 Naraken Chusho Kigyo Sien Center クズイソフラボン量の測定方法
JP2011122950A (ja) 2009-12-11 2011-06-23 Naraken Chusho Kigyo Sien Center イソチオシアネート含量の測定方法
JP5809961B2 (ja) 2011-12-22 2015-11-11 株式会社堀場製作所 試料ガス分析装置及び試料ガス分析装置用プログラム
EP3161452B1 (en) 2014-06-30 2020-05-20 MKS Instruments, Inc. Systems, methods, and apparatus for optical hydrocarbon gas composition monitoring
JP6967387B2 (ja) * 2017-07-14 2021-11-17 株式会社堀場製作所 ガス分析装置、ガス分析装置用プログラム、及びガス分析方法

Also Published As

Publication number Publication date
EP3667295A1 (en) 2020-06-17
US11448588B2 (en) 2022-09-20
JPWO2019031331A1 (ja) 2020-06-25
US20210140877A1 (en) 2021-05-13
WO2019031331A1 (ja) 2019-02-14
EP3667295A4 (en) 2021-05-05

Similar Documents

Publication Publication Date Title
JP6983244B2 (ja) 分析装置、分析方法、分析装置用プログラム及び分析用学習装置
US8862445B2 (en) Selecting spectral elements and components for optical analysis systems
US7248370B2 (en) Method to reduce background noise in a spectrum
US7251037B2 (en) Method to reduce background noise in a spectrum
JP2008539417A (ja) 分析物の混合物における分析物の量を求める分光法
CN103175786B (zh) 样品气体分析装置与样品气体分析装置用方法
US11099124B2 (en) Gas analysis apparatus, program for gas analysis apparatus, and gas analysis method
JP6964114B2 (ja) 分析装置、分析方法及び分析装置用プログラム
WO2023127262A1 (ja) 機械学習装置、排ガス分析装置、機械学習方法、排ガス分析方法、機械学習プログラム、及び、排ガス分析プログラム
JP2741376B2 (ja) 分光分析における多成分分析方法
US9146195B2 (en) Systems and methods for radiance efficiency measurement
JP7335727B2 (ja) 排ガス分析装置、排ガス分析方法、排ガス分析用プログラム及び機械学習装置
Ye et al. Graphical characterization of infrared absorption spectroscopic gas sensor using symmetrized dot pattern
WO2022014126A1 (ja) 分析装置、分析方法、分析装置用プログラム、分析用学習装置、分析用学習方法及び分析用学習装置用プログラム
US20210318280A1 (en) System and computer-implemented method for extrapolating calibration spectra
JP7461948B2 (ja) 試料ガス分析装置、試料ガス分析方法及び試料ガス分析用プログラム
JP7014701B2 (ja) 光学分析装置、並びに光学分析装置に用いられる機械学習装置及びその方法
NL2003432A (nl) Werkwijze voor het bepalen van concentratie-, druk- en temperatuurprofielen in willekeurige, bij voorkeur gasvormige media.
JPH04265839A (ja) 分光分析における多成分分析方法
JPH04262237A (ja) フーリエ変換赤外分光計を用いたエンジン排ガス分析装置
Jaeger Robust on-line engine exhaust methane measurements using wavelength modulation spectroscopy
Branco et al. Synthetic stellar spectra to study multiple populations in globular clusters: an extended grid and the effects on the integrated light
Schröder et al. Emission Measurement with Quantum Cascade Laser
JPH04265841A (ja) 自動車排気ガスの定量方法
JPH04262240A (ja) 分光分析における多成分定量分析方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210720

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210831

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20211027

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20211122

R150 Certificate of patent or registration of utility model

Ref document number: 6983244

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150