JP6049930B1 - ガス分析方法およびガス測定装置 - Google Patents
ガス分析方法およびガス測定装置 Download PDFInfo
- Publication number
- JP6049930B1 JP6049930B1 JP2016060397A JP2016060397A JP6049930B1 JP 6049930 B1 JP6049930 B1 JP 6049930B1 JP 2016060397 A JP2016060397 A JP 2016060397A JP 2016060397 A JP2016060397 A JP 2016060397A JP 6049930 B1 JP6049930 B1 JP 6049930B1
- Authority
- JP
- Japan
- Prior art keywords
- gas
- wave number
- sample
- signal
- absorption
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000004868 gas analysis Methods 0.000 title claims abstract description 24
- 238000010521 absorption reaction Methods 0.000 claims abstract description 86
- 238000005259 measurement Methods 0.000 claims description 11
- 238000001228 spectrum Methods 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims 1
- 239000007789 gas Substances 0.000 description 43
- 238000000862 absorption spectrum Methods 0.000 description 9
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- CSCPPACGZOOCGX-UHFFFAOYSA-N Acetone Chemical compound CC(C)=O CSCPPACGZOOCGX-UHFFFAOYSA-N 0.000 description 2
- 238000002835 absorbance Methods 0.000 description 2
- 229910002092 carbon dioxide Inorganic materials 0.000 description 2
- 239000001569 carbon dioxide Substances 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- QGZKDVFQNNGYKY-UHFFFAOYSA-O Ammonium Chemical compound [NH4+] QGZKDVFQNNGYKY-UHFFFAOYSA-O 0.000 description 1
- RWSOTUBLDIXVET-UHFFFAOYSA-N Dihydrogen sulfide Chemical compound S RWSOTUBLDIXVET-UHFFFAOYSA-N 0.000 description 1
- 230000002902 bimodal effect Effects 0.000 description 1
- 239000002803 fossil fuel Substances 0.000 description 1
- 229910000037 hydrogen sulfide Inorganic materials 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- -1 isotopes) Chemical compound 0.000 description 1
- 230000000155 isotopic effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3504—Investigating 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/39—Investigating 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)
- 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)
- Optics & Photonics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
【課題】吸収率を高い精度で測定可能なガス分析方法を提供する。【解決手段】ガス分析方法は、サンプルガスに含まれる対象ガスの少なくとも１つの吸収線を含む波数帯域内で赤外線の波数を変化させ、前記サンプルガスを透過した前記赤外線の強度の変化を測定してサンプル信号とする第１ステップと、吸収率が低い区間のサンプル信号から暫定リファレンス信号を推定する第２ステップと、前記サンプル信号から前記暫定リファレンス信号を除算して暫定の吸収率の分布を算出する第３ステップと、前記暫定の吸収率が所定値以下となる波数区間を決定し、前記波数区間のサンプル信号に基づいてリファレンス信号を再推定する第４ステップと、前記サンプル信号から再推定された前記リファレンス信号を除算して吸収率の分布を算出する第５ステップと、を有する。【選択図】図３
Description
本発明の実施形態は、ガス分析方法およびガス測定装置に関する。
ガスは、赤外線に対して固有の吸収スペクトルを有する。
サンプルガスを透過した赤外線強度を測定することにより、ガスの吸収係数やガス濃度を算出できる。
波長が３μｍ以上の赤外線は、量子カスケードレーザから放出可能である。但し、レーザ光は、出力変動や波長変動を生じる。基準信号を用いて毎回測定値を補正すると、特性変動による誤差を低減することができる。
環境測定のように長期間にわたってガス分析を行う場合、測定の度に基準信号を測定することが困難なことがある。
赤外線吸収率を高い精度で測定可能なガス分析方法を提供する。
実施形態のガス分析方法は、サンプルガスに含まれる対象ガスの少なくとも１つの吸収線を含む波数帯域内で赤外線の波数を変化させ、前記サンプルガスを透過した前記赤外線の強度の変化を測定してサンプル信号とする第１ステップと、吸収率が低い区間のサンプル信号から暫定リファレンス信号を推定する第２ステップと、前記サンプル信号から前記暫定リファレンス信号を除算して暫定の吸収率の分布を算出する第３ステップと、前記暫定の吸収率が所定値以下となる波数区間を決定し、前記波数区間のサンプル信号に基づいてリファレンス信号を再推定する第４ステップと、前記サンプル信号から再推定された前記リファレンス信号を除算して吸収率の分布を算出する第５ステップと、を有する。
以下に、本発明の各実施の形態について図面を参照しつつ説明する。
図１は、第１の実施形態に係るガス分析方法に用いるガス分析装置を例示する模式図である。
図１に表したように、ガス分析装置１１０は、セル部２０と、光源部３０と、検出部４０と、制御部４５と、を含む。
図１は、第１の実施形態に係るガス分析方法に用いるガス分析装置を例示する模式図である。
図１に表したように、ガス分析装置１１０は、セル部２０と、光源部３０と、検出部４０と、制御部４５と、を含む。
セル部２０には、サンプルガス５０が導入される。すなわち、セル部２０に設けられる空間２３ｓにサンプルガス５０が導入される。サンプルガス５０は、二酸化炭素（同位体を含む）、硫化水素、アセトン、アンモニウム、などを含む。
光源部３０は、空間２３ｓに赤外線３０Ｌを入射させる。また、光源部３０は、駆動部３０ｂにより赤外線３０Ｌの波長を、変化させる。
この例では、光源部３０は、ＱＣＬ（Quantum Cascade Laser:量子カスケードレーザ)３０ａを含む。ＱＣＬ３０ａとして、例えば、分布帰還（ＤＦＢ：Distributed Feed-Back）型ＱＣＬが用いられる。ＱＣＬ３０ａは、端面発光型でも、面発光型でもよい。
検出部４０は、例えば、空間２３ｓにサンプル５０が導入された状態において空間２３ｓを通過した赤外線３０Ｌを検出する。検出部４０は、空間２３ｓを通過した透過光の強度Ｉを検出する。検出部４０には、赤外領域に感度を有する素子が用いられ、たとえば、サーモパイルまたは半導体センサ素子（例えばＩｎＡｓＳｂなどとされる）。
入射光の強度をＩ０、吸光度をＡ、透過率をＴＲ、吸収係数（ｃｍ−１）をα、セル部２０の光路長をＬとするとき、吸光度Ａは、式（１）で表される。
モル濃度をεとするとき、吸収係数（ｃｍ−１）は、式（２）で表される。
また、サンプルガスのモル濃度ｃは、式（３）で表される。
制御部４５は、測定された透過光強度Ｉに基づいて、吸収係数α、吸収率Ｂ、モル濃度ｃ、ＨＩＴＲＡＮデータベースに基づいた吸収スペクトルの計算などを行う。
図２（ａ）〜図２（ｃ）は赤外線の波数を吸収線にチューニングする作用を説明するグラフ図である。すなわち、図２（ａ）はＱＣＬの動作電流の波形図、図２（ｂ）は波長の変化を示すグラフ図、図２（ｃ）は検出器出力を示すグラフ図、である。
図２（ａ）において、縦軸はＱＣＬ３０ａの電流Ｊ、横軸は時間ｔ、を表す。図２（ｂ）において、縦軸は波長、横軸は時間ｔ、を表す。なお、波数は、波長の逆数である。図２（ｃ）において、縦軸は検出器出力Ｓｇ、横軸は時間ｔ、である。ＱＣＬ３０ａに供給される電流Ｊは、時間ｔの増加とともに単調に変化する。図２（ａ）に表すように電流を三角波とすると、波長λは１周期内で時間とともに増加し（図２（ｂ））時間ｔ１で赤外線の波長はλ１となり、時間ｔ２で赤外線の波長はλ２（＞λ１）となる。なお、図２において、右半分は、対象ガスを含むサンプルガスをセルに導入した場合であり、左半分は、対象ガスを含まない基準ガスを導入した場合である。
図２（ａ）において、縦軸はＱＣＬ３０ａの電流Ｊ、横軸は時間ｔ、を表す。図２（ｂ）において、縦軸は波長、横軸は時間ｔ、を表す。なお、波数は、波長の逆数である。図２（ｃ）において、縦軸は検出器出力Ｓｇ、横軸は時間ｔ、である。ＱＣＬ３０ａに供給される電流Ｊは、時間ｔの増加とともに単調に変化する。図２（ａ）に表すように電流を三角波とすると、波長λは１周期内で時間とともに増加し（図２（ｂ））時間ｔ１で赤外線の波長はλ１となり、時間ｔ２で赤外線の波長はλ２（＞λ１）となる。なお、図２において、右半分は、対象ガスを含むサンプルガスをセルに導入した場合であり、左半分は、対象ガスを含まない基準ガスを導入した場合である。
たとえば、波長λ１が第１の吸収線に対応し、波長λ２が第２の吸収線に対応するものとする。それぞれの吸収線の近傍では、対象ガスを含むサンプルガスがセルに導入されている右半分では、赤外線がサンプルガス５０に吸収されるので、サンプルガス５０を透過した赤外線の検出器出力は対象ガスを含まない基準ガスがセルに導入されている左半分に比べて減少する。なお、１つのサンプルガス５０に対して、三角波を繰り返し発生し、検出器出力の平均値を求めると、測定精度を高めることができる。
光源部３０において、波長λの変化は、短時間（例えば１００ｍｓ程度以下）で行われる。例えば、赤外線の波長λの中心値を、４．３４５マイクロメートル以上４．３８４マイクロメータ以下とすると、二酸化炭素の吸収線を含む。この場合、波数帯域は、２２８１ｃｍ−１以上２３０１ｃｍ−１以下となる。
ガス分析の際には、式（１）に用いる（Ｉ/Ｉ０）がこの吸収の有無における検出器出力の比に一致する。そのため、サンプルガスを透過した検出器出力である「サンプル信号」に対して、吸収が無い場合の「リファレンス信号」が必要である。従来技術ではこのリファレンス信号を上記の通り、対象ガスを含まない基準ガスの測定によって求めていたが、本発明ではサンプル信号から求めるため、基準ガスでの測定を不要にする。
図３は、実施形態にかかるガス分析方法のフロー図である。
また、図４（ａ）はサンプル信号の検出器出力、図４（ｂ）は暫定の吸収率、図４（ｃ）は決定された波数区間、図４（ｄ）は決定されたリファレンス推定信号、を表すグラフ図である。
図４（ａ）において、縦軸は検出器出力（Ｖ）、横軸は時間（ｓ）、である。図４（ｂ）において、縦軸は吸収率、横軸は時間（ｓ）、である。図４（ｃ）において、縦軸は吸収率、横軸は時間（ｓ）、である。
また、図４（ａ）はサンプル信号の検出器出力、図４（ｂ）は暫定の吸収率、図４（ｃ）は決定された波数区間、図４（ｄ）は決定されたリファレンス推定信号、を表すグラフ図である。
図４（ａ）において、縦軸は検出器出力（Ｖ）、横軸は時間（ｓ）、である。図４（ｂ）において、縦軸は吸収率、横軸は時間（ｓ）、である。図４（ｃ）において、縦軸は吸収率、横軸は時間（ｓ）、である。
サンプルガス５０に含まれる少なくとも１つの対象ガスの吸収線を含む波数帯域内で赤外線の波数を変化させ、サンプルガス５０を透過した赤外線の透過光強度Ｉの変化を測定してサンプル信号とする（第１ステップＳ１００および図４（ａ））。
吸収率が常に低い複数の区間（レーザ光の波長変動を考慮しても吸収率が常に低くなる時間範囲、「低い」とは、たとえば５×１０−５以下を意味するものとする）のサンプル信号のレベルの例えば２次曲線による内挿から暫定リファレンス信号を推定する（第２ステップＳ１０２）。
サンプル信号から暫定リファレンス信号を除算し、暫定吸収率を算出する（第３ステップＳ１０４および図４（ｂ））。
所定の吸収率（たとえば、５×１０−５）以下の複数の波数区間を決定し、その波数区間のサンプル信号のレベルの例えば２次曲線による内挿からリファレンス信号を再推定する（第４ステップＳ１０６および図４（ｄ））。
サンプル信号から再推定されたリファレンス信号を除算し、吸収率を算出する（第５ステップＳ１０８）。
次に測定するサンプルガスがあれば第１ステップＳ１００へ戻り、なければ終了する（第６ステップＳ１１０）。
図５は、比較例にかかるガス分析方法を説明するグラフ図である。
縦軸は検出器出力（Ｖ）、横軸は時間ｔ、である。比較例では、リファレンス信号を内挿により推定する波数区間Ｆ１１、Ｆ１２などは、時間軸上で固定されている。他方、ＱＣＬの特性は、経時変化、ガスセル内外の環境（温度や圧力など）変化により、吸収のある範囲が時間軸上で、図５に表すように変動する。このため、波数帯域Ｆ１１、Ｆ１２を時間軸上の範囲で固定すると、サンプル信号Ａ１ではリファレンス信号が正しく推定されるが、サンプル信号Ａ２では、推定されたリファレンス信号が本来のリファレンス推定信号よりも過小となり吸収率Ｂ算出における誤差が増大する。
縦軸は検出器出力（Ｖ）、横軸は時間ｔ、である。比較例では、リファレンス信号を内挿により推定する波数区間Ｆ１１、Ｆ１２などは、時間軸上で固定されている。他方、ＱＣＬの特性は、経時変化、ガスセル内外の環境（温度や圧力など）変化により、吸収のある範囲が時間軸上で、図５に表すように変動する。このため、波数帯域Ｆ１１、Ｆ１２を時間軸上の範囲で固定すると、サンプル信号Ａ１ではリファレンス信号が正しく推定されるが、サンプル信号Ａ２では、推定されたリファレンス信号が本来のリファレンス推定信号よりも過小となり吸収率Ｂ算出における誤差が増大する。
これに対して、第１の実施形態にかかるガス分析方法では、吸収率Ｂが所定の値（たとえば５×１０−５以下）よりも低くなる波数区間を、内挿によりリファレンス信号を推定する波数区間としてサンプルガス毎に決定する。このため、内挿範囲に吸収率Ｂが高い波数区間が含まれないので、リファレンス信号の推定精度が損なわれず、吸収率Ｂの測定精度を高めることができる。
次に、波数区間内に同位体元素により生じる少なくとも２つの吸収線がある場合を考える。第２の実施形態にかかるガス分析方法は、双峰性の吸収スペクトル曲線をフィッティングにより２つのスペクトルに分離し、２つの同位体比を算出する。
図６（ａ）〜（ｃ）は同位体比の測定原理を説明するグラフ図である。すなわち、図６（ａ）は時間に対する検出器出力依存性を表すグラフ図、図６（ｂ）は波数に対する吸収係数依存性を表すグラフ図、図６（ｃ）は波数に対するフィッティング後の吸収係数依存性を表すグラフ図、である。
図６（ａ）において、縦軸は検出器出力、横軸は時間、である。図６（ｂ）において、縦軸は吸収係数、横軸は波数（ｃｍ−１）、である。図６（ｃ）において、縦軸は吸収係数、横軸は波数（ｃｍ−１）、である。
図６（ａ）において、縦軸は検出器出力、横軸は時間、である。図６（ｂ）において、縦軸は吸収係数、横軸は波数（ｃｍ−１）、である。図６（ｃ）において、縦軸は吸収係数、横軸は波数（ｃｍ−１）、である。
なお、本図において、同位体は、１３ＣＯ２と１２ＣＯ２とする。図６（ａ）に表すように、第１の実施形態にかかるガス分析方法により、リファレンス推定信号を算出する。図６（ｂ）に表すように、検出器出力の強度分布からリファレンス推定信号の強度分布を除算し、吸収係数分布を算出する。図６（ｃ）に表すように、理論値とのフィッティングにより、１３ＣＯ２と１２ＣＯ２との吸収係数スペクトルを分離する。
次に、吸収スペクトルのフィッティングを行うために、ガス吸収モデルについて説明する。ガス吸収モデルでは、ＨＩＴＲＡＮデータベースの値を使用し、吸収スペクトルを計算する。また、計算は、各吸収線（３１９０本）ごとに、圧力依存性・温度依存性を考慮して行う。
吸収係数αｉは、式（４）により算出される。
それぞれの吸収線強度Ｓ（Ｔ）は、式（５）により算出される。
スペクトル形状関数ｆ（ν、Ｔ、ｐ）は、式（６）により算出される。
図７（ａ）は圧力をパラメータとした波数に対する吸収係数依存性、図７（ｂ）は温度をパラメータとした波数に対する吸収係数依存性、を表すグラフ図である。
図７（ａ）、（ｂ） において、縦軸は吸収係数（ｃｍ−１）、横軸は波数（ｃｍ−１）、である。
図７（ａ）、（ｂ） において、縦軸は吸収係数（ｃｍ−１）、横軸は波数（ｃｍ−１）、である。
図７（ａ）に表すよう、圧力を０．４ａｔｍ、０．５ａｔｍ、０．６ａｔｍと変化したとき吸収係数がピークとなる波数の変動は少ないが、吸収係数の極小値は変動する。また、図７（ｂ）に表すように、温度を２９６Ｋ、３００Ｋ、３０４Ｋと変化したとき吸収係数がピークとなる波数の変動は少ないが、吸収係数の極大値が変動する。このようなガス吸収モデルを用いて、測定値から求められた吸収スペクトルと、理論値とをフィッティングさせることにより、たとえば、２つの同位体の吸収スペクトルを分離して、同位体比の精度を高めることができる。
図８は、フィッティング範囲の決定方法の一例を説明するグラフ図である。
吸収率Ｂのスペクトル曲線のフィッティング範囲は、１３ＣＯ２の吸収線および１２ＣＯ２の吸収線に対して、たとえば、それぞれのピーク強度の２分の１以上かつ信号対雑音比ＳＮＲが１０以上となる領域とすることができる。
吸収率Ｂのスペクトル曲線のフィッティング範囲は、１３ＣＯ２の吸収線および１２ＣＯ２の吸収線に対して、たとえば、それぞれのピーク強度の２分の１以上かつ信号対雑音比ＳＮＲが１０以上となる領域とすることができる。
ＣＯ２の同位体比を（存在比＝１３ＣＯ２／１２ＣＯ２）を％で表すと、小数点以下の桁数が多くなる。このため。天然の同位体比からの乖離度をパーミル（‰）で表すデルタ記法を用いると都合がよい。天然の同位体比を、（１３ＣＯ２／１２ＣＯ２）ref＝１．１２３７２％とする。
式（７）は、デルタ記法による表示式を表す。
たとえば、環境ガス分析を行う場合、１３ＣＯ２と１２ＣＯ２との濃度比をδ１３Ｃ（‰）で表すことができる。対流圏の自由大気のδ１３Ｃは−８‰などであり、化石燃料の濃度比δ１３Ｃは−２８‰などであり、土壌のδ１３Ｃは−２５．７‰などとなるように、ＣＯ２発生源による濃度比δ１３Ｃが異なる。このため、濃度比δ１３Ｃを測定することにより、ＣＯ２の発生源を特定することが可能である。
同様に、火山活動レベルにより濃度比δ１３Ｃの値が異なる。このため、火口近傍で濃度比δ１３Ｃを測定することにより、噴火予知に活用できる。
すなわち、第２の実施形態では、２つの同位体を含む吸収係数の測定値を理論スペクトルとフィッティングさせることにより、同位体比を高精度に求めることができる。
図９は、第３の実施形態にかかるリファレンス信号推定方法を説明するグラフ図である。
第３の実施形態では、図９に表すようにサンプル信号中に２つの同位体の吸収スペクトルがある場合に、それぞれの吸収波数におけるリファレンス信号を、それぞれの両側にある、低吸収区間から個別に内挿する。第１の実施形態では、サンプル信号中の低吸収区間全てから例えば２次曲線一本によるリファレンス推定を行っているため、それぞれの吸収波数においてリファレンス信号の誤差が大きくなるが、本実施形態ではその誤差を抑制することができる。
第３の実施形態では、図９に表すようにサンプル信号中に２つの同位体の吸収スペクトルがある場合に、それぞれの吸収波数におけるリファレンス信号を、それぞれの両側にある、低吸収区間から個別に内挿する。第１の実施形態では、サンプル信号中の低吸収区間全てから例えば２次曲線一本によるリファレンス推定を行っているため、それぞれの吸収波数においてリファレンス信号の誤差が大きくなるが、本実施形態ではその誤差を抑制することができる。
図１０は、第４の実施形態にかかるリファレンス信号推定方法を説明するグラフ図である。
第４の実施形態では、図１０に表すようにサンプル信号中に複数の対象ガスの吸収スペクトルがある場合に、そのピーク時間に波数を対応させて、例えばＨＩＴＲＡＮデータベース等から決定される理論上の低吸収な波数区間のみを細分化して、それらの内挿によりリファレンス信号を推定する。第１の実施形態では、内挿区間に理論上の低吸収ではない範囲が含まれているのに対して、本実施形態ではより低吸収な区間のみから内挿するため、リファレンス信号の推定精度が向上する。
第４の実施形態では、図１０に表すようにサンプル信号中に複数の対象ガスの吸収スペクトルがある場合に、そのピーク時間に波数を対応させて、例えばＨＩＴＲＡＮデータベース等から決定される理論上の低吸収な波数区間のみを細分化して、それらの内挿によりリファレンス信号を推定する。第１の実施形態では、内挿区間に理論上の低吸収ではない範囲が含まれているのに対して、本実施形態ではより低吸収な区間のみから内挿するため、リファレンス信号の推定精度が向上する。
本発明の実施形態によれば、 吸収率を高い精度で測定可能なガス分析方法が提供される。また、透過光強度の測定値を補正するリファレンス推定信号は、吸収率の強度分布に基づいてそれぞれのサンプルガス毎に推定される。このため、基準ガスを用いたリファレンス信号を測定する必要がなく、分析装置を簡素な構成とできる。
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。
３０ 光源、３０ａ ＱＣＬ、３０Ｌ 赤外線、４０ 検出器、５０ サンプルガス、Ｓ 吸収線強度、Ｂ 吸収率、α、αi 吸収係数、Ｆ１、Ｆ２、Ｆ３ 波数区間、ＳＮＲ 信号対雑音比
Claims (9)
- サンプルガスに含まれる対象ガスの少なくとも１つの吸収線を含む波数帯域内で赤外線の波数を変化させ、前記サンプルガスを透過した前記赤外線の強度の変化を測定してサンプル信号とする第１ステップと、
吸収率が低い区間のサンプル信号から暫定リファレンス信号を推定する第２ステップと、
前記サンプル信号から前記暫定リファレンス信号を除算して暫定の吸収率の分布を算出する第３ステップと、
前記暫定の吸収率が所定値以下となる波数区間を決定し、前記波数区間のサンプル信号に基づいてリファレンス信号を再推定する第４ステップと、
前記サンプル信号から再推定された前記リファレンス信号を除算して吸収率の分布を算出する第５ステップと、
を備えたガス分析方法。 - 時間とともに波数が単調に変化する、請求項１記載のガス分析方法。
- 前記波数帯域内に同位体元素により生じる少なくとも２つの吸収線がある場合、前記吸収率をフィッティングにより少なくとも２つのスペクトルに分離する、請求項１または２に記載のガス分析方法。
- 前記吸収率のスペクトル曲線のフィッティング範囲は、それそれの吸収線に対して、それぞれのピーク強度のおよそ２分の１以上かつ信号対雑音比が約１０以上となるスペクトル強度領域とする請求項３記載のガス分析方法。
- 前記赤外線は、量子カスケードレーザから放出される、請求項１〜４のいずれか１つに記載のガス分析方法。
- サンプルガスに含まれる対象ガスの少なくとも１つの吸収線を含む波長帯域内で赤外線の波数を変化させる光源部と、
前記サンプルガスが導入されるセル部と、
前記サンプルガスを透過した前記赤外線の強度を測定する検出部と、
前記検出部により測定された前記赤外線の前記強度に基づいて、前記赤外線の吸収率の分布を算出する制御部であって、吸収率が低い区間のサンプル信号から暫定リファレンス信号を推定し、前記サンプル信号から前記暫定リファレンス信号を除算して暫定の吸収率の分布を算出し、前記暫定の吸収率が所定値以下となる波数区間を決定し、前記波数区間のサンプル信号に基づいてリファレンス信号を再推定し、前記サンプル信号から再推定された前記リファレンス信号を除算して吸収率の分布を算出する、制御部と、
を備えたガス測定装置。 - 前記光源部は、前記赤外線を放出可能な量子カスケードレーザを有する請求項６記載のガス測定装置。
- 前記光源部は、時間とともに波数を単調に変化させる請求項６または７に記載のガス測定装置。
- 前記制御部は、前記波数帯域内に同位体元素により生じる少なくとも２つの吸収線がある場合、前記吸収率をフィッティングにより少なくとも２つのスペクトルに分離する、請求項６〜８のいずかれ１つに記載のガス測定装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016060397A JP6049930B1 (ja) | 2016-03-24 | 2016-03-24 | ガス分析方法およびガス測定装置 |
| PCT/JP2016/075858 WO2017163452A1 (ja) | 2016-03-24 | 2016-09-02 | ガス分析方法およびガス測定装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016060397A JP6049930B1 (ja) | 2016-03-24 | 2016-03-24 | ガス分析方法およびガス測定装置 |
Related Child Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2016226443A Division JP6293855B2 (ja) | 2016-11-22 | 2016-11-22 | 物質分析方法および分析装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP6049930B1 true JP6049930B1 (ja) | 2016-12-21 |
| JP2017173162A JP2017173162A (ja) | 2017-09-28 |
Family
ID=57572466
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2016060397A Active JP6049930B1 (ja) | 2016-03-24 | 2016-03-24 | ガス分析方法およびガス測定装置 |
Country Status (2)
| Country | Link |
|---|---|
| JP (1) | JP6049930B1 (ja) |
| WO (1) | WO2017163452A1 (ja) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7334502B2 (ja) * | 2018-09-19 | 2023-08-29 | 富士電機株式会社 | レーザ式ガス分析計 |
| WO2020059100A1 (ja) * | 2018-09-21 | 2020-03-26 | 大塚電子株式会社 | 測定装置、および測定方法 |
| JP7147870B2 (ja) * | 2018-11-20 | 2022-10-05 | 株式会社島津製作所 | ガス測定装置及びガス測定方法 |
| EP3884249A4 (en) * | 2018-11-21 | 2022-08-03 | Thermo Fisher Scientific Inc. | SYSTEM AND METHOD FOR RAPID AND ACCURATE MEASUREMENT OF TRACE GAS |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001059772A (ja) * | 1999-08-25 | 2001-03-06 | Jasco Corp | ベースライン補正方法 |
| JP2001343324A (ja) * | 2000-06-01 | 2001-12-14 | Advantest Corp | 赤外線吸光スペクトルのベースライン補正方法及びそのプログラム記録媒体 |
| JP2003232732A (ja) * | 2002-02-06 | 2003-08-22 | Shimadzu Corp | 同位体ガス測定装置 |
| JP2010060387A (ja) * | 2008-09-02 | 2010-03-18 | Jasco Corp | ベースライン設定方法 |
| JP2011013126A (ja) * | 2009-07-03 | 2011-01-20 | Shimadzu Corp | ガス濃度測定装置 |
-
2016
- 2016-03-24 JP JP2016060397A patent/JP6049930B1/ja active Active
- 2016-09-02 WO PCT/JP2016/075858 patent/WO2017163452A1/ja active Application Filing
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2001059772A (ja) * | 1999-08-25 | 2001-03-06 | Jasco Corp | ベースライン補正方法 |
| JP2001343324A (ja) * | 2000-06-01 | 2001-12-14 | Advantest Corp | 赤外線吸光スペクトルのベースライン補正方法及びそのプログラム記録媒体 |
| JP2003232732A (ja) * | 2002-02-06 | 2003-08-22 | Shimadzu Corp | 同位体ガス測定装置 |
| JP2010060387A (ja) * | 2008-09-02 | 2010-03-18 | Jasco Corp | ベースライン設定方法 |
| JP2011013126A (ja) * | 2009-07-03 | 2011-01-20 | Shimadzu Corp | ガス濃度測定装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2017163452A1 (ja) | 2017-09-28 |
| JP2017173162A (ja) | 2017-09-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6049930B1 (ja) | ガス分析方法およびガス測定装置 | |
| US9618391B2 (en) | Collisional broadening compensation using real or near-real time validation in spectroscopic analyzers | |
| EP2597456B1 (en) | Laser gas analyzer | |
| US11346774B2 (en) | Infra red spectrometer | |
| Kluczynski et al. | Detection of acetylene impurities in ethylene and polyethylene manufacturing processes using tunable diode laser spectroscopy in the 3-μm range | |
| Phillips et al. | Real-time trace gas sensing of fluorocarbons using a swept-wavelength external cavity quantum cascade laser | |
| JP6293855B2 (ja) | 物質分析方法および分析装置 | |
| JP6981817B2 (ja) | 分光分析装置及び分光分析方法 | |
| JP2016035385A (ja) | レーザガス分析装置 | |
| TWI732264B (zh) | 具有用於干擾弱化之交錯的資料擷取之共振腔環路衰減光譜術 | |
| Tátrai et al. | Dual-channel photoacoustic hygrometer for airborne measurements: background, calibration, laboratory and in-flight intercomparison tests | |
| KR20230112634A (ko) | 가스 분석 장치 및 가스 분석 방법 | |
| Song et al. | A four-channel-based mid-infrared methane sensor system using novel optical/electrical dual-domain self-adaptive denoising algorithm | |
| JP5170034B2 (ja) | ガス分析装置 | |
| JP6530669B2 (ja) | ガス濃度測定装置 | |
| CN114878515A (zh) | 一种可调谐二极管激光吸收信号时间-频率转换方法 | |
| CN112858206A (zh) | 一种基于可调谐fpi的中红外气体测量方法及装置 | |
| JP6201551B2 (ja) | ガス分析装置 | |
| JP4150760B2 (ja) | アイソトポマー存在比測定方法 | |
| Bozóki | Airborne application of a photoacoustic instrument | |
| Pushkarsky et al. | Performance characteristics of a compact widely tunable external cavity quantum cascade laser | |
| JP6314605B2 (ja) | ガス分析装置 | |
| JP6571476B2 (ja) | ガス濃度測定装置 | |
| Ivantsov et al. | Study of Optic Gas Sensor Calibration Problems | |
| JP2023080041A (ja) | 濃度測定装置および濃度測定方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20161003 |
|
| TRDD | Decision of grant or rejection written | ||
| A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20161014 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20161024 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20161122 |
|
| R151 | Written notification of patent or utility model registration |
Ref document number: 6049930 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |