(57)【要約】【課題】  レーザダイオードの発振波長が吸収ピークの周波数上限νｍａｘから周波数下限νｍｉｎまでの範囲をカバーできないときにも、水分の分圧値ＰＨ２Ｏを正確に求めることができるガス分析装置を提供する。【解決手段】  測定対象ガス中の特定ガス量情報を算出する演算部６１ａを備えるガス分析装置１であって、演算部６１ａは、分割測定光の強度変化Ｉ（ν）と参照光の強度変化Ｉ０（ν）とに基づいて、ｌｎ（Ｉ０（ν）／Ｉ（ν））を作成し、最小二乗法を用いてｌｎ（Ｉ０（ν）／Ｉ（ν））から近似式［ｌｎ（Ｉ０（ν）／Ｉ（ν））］を作成した後、［ｌｎ（Ｉ０（ν）／Ｉ（ν））］が最大値となる周波数νをν０ｎｅｗとしたときの周波数ν１ｎｅｗと周波数ν２ｎｅｗとを求め、［ｌｎ（Ｉ０（ν０ｎｅｗ）／Ｉ（ν０ｎｅｗ））］と周波数ν１ｎｅｗと周波数ν２ｎｅｗとに基づいて、測定対象ガス中の特定ガス量情報を算出する。【選択図】図１
2015055486.tif 000002
【特許請求の範囲】
【請求項１】
ガスセル内の測定対象ガスに、特定ガスに吸収される吸収波長を含む測定光を照射する光源部と、
ガスセル内の測定対象ガスに照射する分割測定光と、測定対象ガスに照射しない参照光とに分割する測定光分割部と、
前記測定対象ガス中を通過した分割測定光の強度を受光する分割測定光受光部と、
前記参照光の強度を受光する参照光受光部と、
分割測定光の強度変化Ｉ（ν）と参照光の強度変化Ｉ０（ν）とに基づいて、前記測定対象ガス中の特定ガス量情報を算出する演算部とを備えるガス分析装置であって、
前記演算部は、分割測定光の強度変化Ｉ（ν）と参照光の強度変化Ｉ０（ν）とに基づいて、ｌｎ（Ｉ０（ν）／Ｉ（ν））を作成し、最小二乗法を用いてｌｎ（Ｉ０（ν）／Ｉ（ν））から近似式［ｌｎ（Ｉ０（ν）／Ｉ（ν））］を作成した後、［ｌｎ（Ｉ０（ν）／Ｉ（ν））］が最大値となる周波数νをν０ｎｅｗとしたときに、下記式の関係が成立する周波数ν１ｎｅｗと周波数ν２ｎｅｗとを求め、［ｌｎ（Ｉ０（ν０ｎｅｗ）／Ｉ（ν０ｎｅｗ））］と周波数ν１ｎｅｗと周波数ν２ｎｅｗとに基づいて、前記測定対象ガス中の特定ガス量情報を算出することを特徴とするガス分析装置。
2015055486.tif 000016
【請求項２】
前記演算部は、周波数ν１ｎｅｗと周波数ν２ｎｅｗとを含む分割測定光の強度変化Ｉ（ν）と、周波数ν１ｎｅｗと周波数ν２ｎｅｗとを含む参照光の強度変化Ｉ０（ν）とに基づいて、ｌｎ（Ｉ０（ν）／Ｉ（ν））を作成することを特徴とする請求項１に記載のガス分析装置。
【請求項３】
［ｌｎ（Ｉ０（ν）／Ｉ（ν））］は、下記式で表されることを特徴とする請求項１又は請求項２に記載のガス分析装置。
［ｌｎ（Ｉ０（ν）／Ｉ（ν））］＝１／（ａν２＋ｂν＋ｃ）  ・・・（１３）
ただし、ａ、ｂ、ｃは定数である。
【請求項４】
前記演算部は、［ｌｎ（Ｉ０（ν）／Ｉ（ν））］におけるａ、ｂ、ｃを記憶させ表示させることを特徴とする請求項１〜請求項３のいずれか１項に記載のガス分析装置。
【請求項５】
前記演算部は、Ｌｏｒｅｎｔｚ型の吸収波形プロファイル関数Φ（ν）又はＧａｕｓｓ型の吸収波形プロファイル関数Φ（ν）を下記式に代入することにより、前記測定対象ガス中の特定ガス量情報を算出することを特徴とする請求項１〜請求項４のいずれか１項に記載のガス分析装置。
2015055486.tif 000017
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、レーザ吸収分光法を利用して気体中の特定ガス量情報を測定するガス分析装置に関し、特に、半導体製造装置における真空領域中や煙道中や燃焼プロセス中や自動車測定対象ガス中や燃料電池における流路中等の水蒸気量（特定ガス量情報）を測定するガス分析装置に関する。
【背景技術】
【０００２】
気体中の水蒸気量を測定する方法の一つとして、水分子が特定波長領域（例えば、１．３μｍ帯）の光のみを吸収することを利用した吸収分光法が挙げられる。この吸収分光法は、測定対象のサンプルガスに対し非接触で測定可能であるため、サンプルガスの場を乱さずにサンプルガス中の水蒸気量を計測することができる。
【０００３】
このような吸収分光法の中でも、特に光源に波長可変半導体レーザを利用した「波長可変半導体レーザ吸収分光法」は、シンプルな装置構成で実現することができる。
例えば、「波長可変半導体レーザ吸収分光法」を利用したガス分析装置では、サンプルガスが所定方向に流れている配管（ガスセル）に対して、配管に形成された入射用光学窓と出射用光学窓とを介して、配管を横切って光路（光路長Ｌ）が形成されるようにそれぞれ対向して設けられる波長可変半導体レーザと光検出センサ（受光部）とを追加することが一般的である（例えば、特許文献１参照）。
【０００４】
このようなガス分析装置によれば、波長可変半導体レーザから照射された所定波長のレーザ光は、配管内を通過する過程でサンプルガス中に存在する水分子の遮光作用によってレーザ光の進行が阻害されて、サンプルガス中における水分子の濃度に対応して光検出センサに入射する光量が減少することを利用して、波長可変半導体レーザから放射されたレーザ光の光量に対する光検出センサに入射するレーザ光の光量を測定することによって水分子の濃度が算出される。図６は、ガス分析装置で得られた吸収スペクトルの一例を示すグラフである。縦軸は受光強度Ｉであり、横軸は周波数νである。なお、Ｉ０（ν）は周波数νにおいて水分子の吸収を受けなかった場合の受光強度Ｉであり、非吸収波長の受光強度Ｉに基づいて近似式を作成することで導出されることになる。
【０００５】
ここで、図６に示す吸収スペクトルを用いた演算処理の一例について説明する。Lambert-Beerの法則より下記式（１）が成り立つ。
【０００６】
【数１】
2015055486.tif 000003
【０００７】
なお、Ｉ０（ν）は周波数νにおいて水分子の吸収を受けなかった場合の光強度、Ｉ（ν）は周波数νにおける透過光強度、ｃ（ｍｏｌ／ｃｍ３）は水分子の数密度、Ｌ（ｃｍ）はサンプルガスを通過する光路の長さ、Ｓ（Ｔ）（ｃｍ−１／（ｍｏｌ／ｃｍ−２））は所定の吸収線強度におけるガス温度Ｔの関数である。
【０００８】
ここで、図７は、縦軸をｌｎ（Ｉ０（ν０）／Ｉ（ν０））とし、横軸を周波数νとしたグラフである。よって、式（１）の左辺の値は、図７に示すグラフの面積を求めることで得られる。図７のグラフの面積を求める方法として、長方形近似を一例に挙げると、式（１）の左辺は下記式（２）のように変形することができる。
【０００９】
【数２】
2015055486.tif 000004
【００１０】
なお、νｍａｘは吸収帯（吸収ピーク）の周波数上限、νｍｉｎは吸収帯の周波数下限、ｎは１波形あたりの測定点数である。
【００１１】
一方、式（１）の右辺におけるＳ（Ｔ）に関しては下記式（３）が成り立つ。
【００１２】
【数３】
2015055486.tif 000005
【００１３】
なお、Ｓ０は標準状態での線強度、Ｑ（Ｔ）は分配関数、Ｂ（Ｔ）はボルツマン因子、ＳＥ（Ｔ）は誘導放射の補正式である。
さらに、式（３）の右辺におけるＱ（Ｔ）、Ｂ（Ｔ）、ＳＥ（Ｔ）は、それぞれ下記式（４）、（５）、（６）のように表すことができる。
【００１４】
【数４】
2015055486.tif 000006
【００１５】
なお、Ｓ０、定数ａ〜ｄ、Ｅｌは、ＨＩＴＲＡＮデータベース等から得られることができる。よって、ガス温度値Ｔと光強度変化Ｉ（ν）、Ｉ０（ν）とを得ることができれば、水分子の数密度ｃが算出できることになる。
【００１６】
そして、水分子の数密度ｃと水分の分圧値ＰＨ２Ｏとの関係は、下記式（７）のように表すことができる。
【００１７】
【数５】
2015055486.tif 000007
【００１８】
なお、ｋはボルツマン定数である。これにより、水分の分圧値ＰＨ２Ｏ（水蒸気量情報）を算出することができる。
【００１９】
ここで、図４は、波長可変半導体レーザ吸収分光法を利用したガス分析装置の一例を示す概略構成図であり、図５は、図４に示すガス分析装置のブロック図である。なお、地面に水平な一方向をＸ方向とし、地面に水平でＸ方向と垂直な方向をＹ方向とし、Ｘ方向とＹ方向とに垂直な方向をＺ方向とする。
ガス分析装置１０１は、光源部１１０と、受光部１２０と、圧力値Ｐｔｏｔａｌを測定する圧力センサ３１と、ガス温度値Ｔを測定するガス温度センサ３２と、光源部１１０を制御するレーザ制御部５０と、マイコンやＰＣで構成される制御部１６０とを備える。
【００２０】
このようなガス分析装置１０１は、燃料電池システムへの給排気の各ラインに連結されたサンプル流路７０内を流れるサンプルガス中の水分の分圧値ＰＨ２Ｏを測定するために用いられている。サンプル流路７０はＺ方向に伸びており、サンプル流路７０の側壁には、入射用光学窓７１と、入射用光学窓７１にＸ方向に距離Ｌを空けて対向する出射用光学窓７２とが形成されている。そして、サンプルガスはサンプル流路７０内をＺ方向に流れている。
【００２１】
光源部１１０は、例えば光通信用分布帰還系形（ＤＦＢ：distributed feedback）半導体レーザダイオードが用いられる。レーザダイオードは、近赤外領域の所定波長範囲（位相速度定数Ａ／νｅｎｄ〜位相速度定数Ａ／νｓｔａｒｔ）内においてレーザの発振波長を調整することが可能となっており、レーザ制御部５０からの制御信号によって制御されるようになっている。なお、上記光源部として、波長可変型の半導体レーザであれば赤外光その他の波長領域のいずれの構造のレーザでもよく、また、比較的高価ではあるが、量子カスケードレーザ等が用いてもよい。
そして、レーザダイオードは、入射用光学窓７１からサンプル流路７０内にＸ方向でレーザ光を入射させるように配置されており、レーザ光がサンプルガスに対して照射されるようになっている。
【００２２】
受光部１２０は、光強度を電気信号に変換できるものであればよく、例えばフォトダイオードが用いられる。そして、フォトダイオードは、出射用光学窓７２からサンプル流路７０外にＸ方向で出射されたレーザ光を受光するように配置されており、サンプルガスを通過したレーザ光の強度Ｉ（ν）を受光する。これにより、水分子の吸収スペクトルの中心波長部分のレーザ光の強度Ｉ（ν）と、中心波長部分の両側となる非吸収波長部分のレーザ光の強度Ｉ（ν）とを含むスペクトル波形をフォトダイオードにより取得することで、制御部１６０がＩ０（ν）とＩ（ν）とを算出するようになっている。
【００２３】
圧力センサ３１は、サンプル流路７０内に設置されており、サンプルガスの全圧である圧力値Ｐｔｏｔａｌを所定時間間隔で測定する。また、ガス温度センサ３２も、サンプル流路７０内に設置されており、サンプルガスの温度であるガス温度値Ｔを所定時間間隔で測定する。
【００２４】
レーザ制御部５０は、レーザ電流を制御するレーザ電流制御部５１と、レーザ温度を制御するレーザ温調部５２とにより構成される。
制御部１６０は、ＣＰＵ１６１とメモリ１６２と表示部６３と入力装置６４とを備える。また、ＣＰＵ１６１が処理する機能をブロック化して説明すると、サンプルガス中の分圧値ＰＨ２Ｏを算出する演算部１６１ａを有する。
【００２５】
そして、レーザ光の強度Ｉ（ν）、圧力値Ｐｔｏｔａｌ、ガス温度値ＴをそれぞれＡ／Ｄ変換部１、２、３によってデジタル値に変換し、演算部１６１ａは、このデジタル値からＩ０（ν）とＩ（ν）とを算出して式（１）に当てはめて数密度ｃを得て、得られた数密度ｃを式（７）に当てはめて分圧値ＰＨ２Ｏを得ていく。
なお、演算部１６１ａは、分圧値ＰＨ２Ｏと圧力値Ｐｔｏｔａｌとを用いて、サンプルガス中の水蒸気量（濃度）を演算してもよい。
【先行技術文献】
【特許文献】
【００２６】
【特許文献１】特開２０１０−２３７０７５号公報
【発明の概要】
【発明が解決しようとする課題】
【００２７】
ところで、上述したようなガス分析装置１０１により、図６に示すような吸収スペクトルを得ることになるが、サンプル流路７０内を流れるサンプルガスの圧力値Ｐｔｏｔａｌが高圧であるときや、サンプルガスの濃度が高濃度であるときには、レーザダイオードの発振波長を最大限に調整しても、（νｓｔａｒｔがνｍｉｎより大きくなったりνｅｎｄがνｍａｘより小さくなったりして）吸収スペクトル中の吸収ピークの周波数上限νｍａｘから周波数下限νｍｉｎまでの範囲をカバーできず、その結果、Ｉ０（ν）を正確に算出できなかったり、図７に示すようなグラフの面積を正確に求められないことがあった。
【００２８】
また、図６に示すような吸収スペクトルからＩ０（ν）を算出するのではなく、レーザ光を分割測定光と参照光とに分割することにより参照光の強度Ｉ０（ν）を受光するダブルビーム法を用いた場合、ダブルビーム法では、Ｉ０（ν）を正確に算出することはできても、やはり吸収ピークの周波数上限νｍａｘから周波数下限νｍｉｎまでの範囲をカバーできなければ、図７のグラフの面積を正確に求めることはできなかった。
さらに、ダブルビーム法によりグラフの面積を求めるのではなく、Ｉ（ν）の最大強度値を求め、Ｉ（ν）の最大強度値から水分の分圧値ＰＨ２Ｏを算出することも可能ではあるが、この算出方法では吸収ピークは圧力値Ｐｔｏｔａｌ、ガス温度値Ｔ、ガス濃度といった複数のパラメータの影響を受けて変化するため、精度が低くなるという問題点があった。
【課題を解決するための手段】
【００２９】
本出願人は、レーザダイオードの発振波長が吸収ピークの周波数上限νｍａｘから周波数下限νｍｉｎまでの範囲をカバーできないときにも、水分の分圧値ＰＨ２Ｏを正確に求めることができる方法について検討した。まず、Ｉ０（ν）を正確に算出するために、ダブルビーム法を用いることにした。次に、吸収スペクトル中の吸収ピークについては、サンプルガスの圧力値Ｐｔｏｔａｌの範囲により幾つかの数式で表現できる曲線となることがわかっている。そこで、ｌｎ（Ｉ０（ν）／Ｉ（ν））は、吸収波形プロファイル関数Φ（ν）を用いて下記式（８）のように表すことにした。
【００３０】
【数６】
2015055486.tif 000008
【００３１】
言い換えると、吸収波形プロファイル関数Φ（ν）は、周波数ν全域にわたって積分すると１になる関数である。そして、サンプルガスの圧力値Ｐｔｏｔａｌが大気圧（１０１．３ｋＰａ）以上である場合には、Ｌｏｒｅｎｔｚ型の吸収波形プロファイル関数Φ（ν）をとることが知られており、下記式（９）が成り立つ。なお、図３は、縦軸をｌｎ（Ｉ０（ν０）／Ｉ（ν０））とし、横軸を周波数νとしたグラフであり、ｌｎ（Ｉ０（ν）／Ｉ（ν））が最大値となる周波数νをν０とすると、周波数ν１と周波数ν２とには下記式（１０）の関係が成立する。
【００３２】
【数７】
2015055486.tif 000009
【００３３】
ここで、周波数ν０において式（８）は下記式（１１）のようになる。
【００３４】
【数８】
2015055486.tif 000010
【００３５】
よって、ｌｎ（Ｉ０（ν０）／Ｉ（ν０））と周波数ν１と周波数ν２とがわかれば、式（９）と式（１１）とを用いて水分子の数密度ｃが算出でき、水分の分圧値ＰＨ２Ｏが求められることになる。
【００３６】
一方、サンプルガスの圧力値Ｐｔｏｔａｌが大気圧（１３３Ｐａ）未満である場合には、Ｇａｕｓｓ型の吸収波形プロファイル関数Φ（ν）をとることが知られており、下記式（１２）が成り立つ。
【００３７】
【数９】
2015055486.tif 000011
【００３８】
よって、ｌｎ（Ｉ０（ν０）／Ｉ（ν０））と周波数ν１と周波数ν２とがわかれば、式（１１）と式（１２）とを用いて水分子の数密度ｃが算出でき、水分の分圧値ＰＨ２Ｏが求められることになる。
したがって、本出願人は、ｌｎ（Ｉ０（ν０）／Ｉ（ν０））と周波数ν１と周波数ν２とを用いて水分の分圧値ＰＨ２Ｏを求める方法を見出した。これにより、吸収スペクトル中の吸収ピークの周波数上限νｍａｘから周波数下限νｍｉｎまでの範囲をカバーできない場合に、図７のグラフの面積を求める方法では、吸収ピークの裾野部分の取りこぼしがあったが、本発明ではｌｎ（Ｉ０（ν０）／Ｉ（ν０））と周波数ν１と周波数ν２とを用いているので、水分の分圧値ＰＨ２Ｏを正確に算出することができる。また、吸収ピークの裾野部分が他の吸収ピークと重なっている場合においても、他の吸収ピークの影響を受けることを抑制することができる。
【００３９】
ところで、上述したような方法では、ｌｎ（Ｉ０（ν）／Ｉ（ν））における周波数ν０と周波数ν１と周波数ν２とを求めることが重要となる。しかしながら、電気的なノイズやレーザ特有の光学的な干渉ノイズが発生した場合には、ｌｎ（Ｉ０（ν）／Ｉ（ν））のグラフは図８に示すようになる。図８に示すように吸収の中心波長付近で何らかのノイズが混入した場合、ｌｎ（Ｉ０（ν）／Ｉ（ν））の最大値はノイズによって大きく誤って検出され、さらにその周波数νを吸収の中心波長とみなしてしまうことになる。また、光学的な干渉に由来する揺らぎもみられ、半値幅の算出誤差に大きく影響する。さらに、光学的な干渉ノイズが周波数νに対して定在波となって出現し、時間的な平均化処理で取り除けないノイズとなるため、ｌｎ（Ｉ０（ν）／Ｉ（ν））における周波数ν０と周波数ν１と周波数ν２とを求めることが困難となっている。
【００４０】
ここで、ｌｎ（Ｉ０（ν）／Ｉ（ν））は、吸収波形プロファイル関数Φ（ν）で表されるので、周波数νの２次関数となる。そこで、計測した離散的な測定点数ｎのｌｎ（Ｉ０（ν）／Ｉ（ν））から最小二乗法を用いて近似式［ｌｎ（Ｉ０（ν）／Ｉ（ν））］を作成することにした。実際には、一般的な２次の多項式近似処理を適用するために、ｌｎ（Ｉ０（ν）／Ｉ（ν））の逆数を、最小二乗法を用いて下記式（１３）のように表すことにした。
［ｌｎ（Ｉ（ν）／Ｉ０（ν））］＝ａν２＋ｂν＋ｃ  ・・・（１３）
ただし、ａ、ｂ、ｃは定数である。
【００４１】
そして、式（１３）の逆数を求めると、下記式（１４）のようになる。
［ｌｎ（Ｉ０（ν）／Ｉ（ν））］＝１／（ａν２＋ｂν＋ｃ）  ・・・（１４）
ただし、ａ、ｂ、ｃは定数である。
【００４２】
ここで、図９は、［ｌｎ（Ｉ０（ν）／Ｉ（ν））］のグラフである。［ｌｎ（Ｉ０（ν）／Ｉ（ν））］が最大値となる周波数νをν０ｎｅｗとすると、周波数ν１ｎｅｗと周波数ν２ｎｅｗとには下記式（１５）の関係が成立する。
【００４３】
【数１０】
2015055486.tif 000012
【００４４】
これにより、電気的なノイズ等の影響を低減できた近似式［ｌｎ（Ｉ０（ν）／Ｉ（ν））］における周波数ν０ｎｅｗと周波数ν１ｎｅｗと周波数ν２ｎｅｗとが求められ、その結果、測定対象ガス中の特定ガス量情報を正確に測定することができる。
【００４５】
すなわち、本発明のガス分析装置は、ガスセル内の測定対象ガスに、特定ガスに吸収される吸収波長を含む測定光を照射する光源部と、ガスセル内の測定対象ガスに照射する分割測定光と、測定対象ガスに照射しない参照光とに分割する測定光分割部と、前記測定対象ガス中を通過した分割測定光の強度を受光する分割測定光受光部と、前記参照光の強度を受光する参照光受光部と、分割測定光の強度変化Ｉ（ν）と参照光の強度変化Ｉ０（ν）とに基づいて、前記測定対象ガス中の特定ガス量情報を算出する演算部とを備えるガス分析装置であって、前記演算部は、分割測定光の強度変化Ｉ（ν）と参照光の強度変化Ｉ０（ν）とに基づいて、ｌｎ（Ｉ０（ν）／Ｉ（ν））を作成し、最小二乗法を用いてｌｎ（Ｉ０（ν）／Ｉ（ν））から近似式［ｌｎ（Ｉ０（ν）／Ｉ（ν））］を作成した後、［ｌｎ（Ｉ０（ν）／Ｉ（ν））］が最大値となる周波数νをν０ｎｅｗとしたときに、下記式の関係が成立する周波数ν１ｎｅｗと周波数ν２ｎｅｗとを求め、ｌｎ（Ｉ０（ν０ｎｅｗ）／Ｉ（ν０ｎｅｗ））と周波数ν１ｎｅｗと周波数ν２ｎｅｗとに基づいて、前記測定対象ガス中の特定ガス量情報を算出するようにしている。
【００４６】
【数１１】
2015055486.tif 000013
【発明の効果】
【００４７】
以上のように、本発明のガス分析装置によれば、光源部（レーザ）のもつ波長制御範囲にとらわれずに、測定対象ガスの圧力が高圧であるときや測定対象ガスの濃度が高濃度であるときにも、測定対象ガス中の特定ガス量情報を測定することができる。このとき、近似式［ｌｎ（Ｉ０（ν）／Ｉ（ν））］を作成しているため、電気的なノイズ等の影響を低減でき、適切な周波数ν０ｎｅｗと周波数ν１ｎｅｗと周波数ν２ｎｅｗとを求めているので、測定対象ガス中の特定ガス量情報を正確に測定することができる。
【００４８】
（他の課題を解決するための手段及び効果）
また、上記の発明では、前記演算部は、周波数ν１ｎｅｗと周波数ν２ｎｅｗとを含む分割測定光の強度変化Ｉ（ν）と、周波数ν１ｎｅｗと周波数ν２ｎｅｗとを含む参照光の強度変化Ｉ０（ν）とに基づいて、ｌｎ（Ｉ０（ν）／Ｉ（ν））を作成するようにしてもよい。
【００４９】
また、上記の発明では、［ｌｎ（Ｉ０（ν）／Ｉ（ν））］は、下記式で表されるようにしてもよい。
［ｌｎ（Ｉ０（ν）／Ｉ（ν））］＝１／（ａν２＋ｂν＋ｃ）  ・・・（１３）
ただし、ａ、ｂ、ｃは定数である。
【００５０】
そして、上記の発明では、前記演算部は、［ｌｎ（Ｉ０（ν）／Ｉ（ν））］におけるａ、ｂ、ｃを記憶させて表示させるようにしてもよい。
本発明のガス分析装置によれば、計測時の吸収スペクトル波形情報を多項式近似の係数という極めて少ない数字として変換できるため、演算結果（特定ガス量情報）とともに付加情報としてこれらの情報を保存しておけば、簡便に計測当時の吸光状態を確認することができ、装置の異常判定に役立つ。
【００５１】
さらに、上記の発明では、前記演算部は、Ｌｏｒｅｎｔｚ型の吸収波形プロファイル関数Φ（ν）又はＧａｕｓｓ型の吸収波形プロファイル関数Φ（ν）を下記式に代入することにより、前記測定対象ガス中の特定ガス量情報を算出するようにしてもよい。
【００５２】
【数１２】
2015055486.tif 000014
【図面の簡単な説明】
【００５３】
【図１】本発明の一実施形態であるガス分析装置の一例を示す概略構成図。
【図２】図１に示すガス分析装置のブロック図。
【図３】縦軸をｌｎ（Ｉ０（ν０）／Ｉ（ν０））とし、横軸を周波数νとしたグラフ。
【図４】波長可変半導体レーザ吸収分光法を利用したガス分析装置の一例を示す概略構成図。
【図５】図４に示すガス分析装置のブロック図。
【図６】ガス分析装置で得られた吸収スペクトルの一例を示すグラフ。
【図７】縦軸をｌｎ（Ｉ０（ν０）／Ｉ（ν０））とし、横軸を周波数νとしたグラフ。
【図８】ｌｎ（Ｉ０（ν）／Ｉ（ν））のグラフ。
【図９】［ｌｎ（Ｉ０（ν）／Ｉ（ν））］のグラフ。
【発明を実施するための形態】
【００５４】
以下、本発明の実施形態について図面を用いて説明する。なお、本発明は、以下に説明するような実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で種々の態様が含まれることはいうまでもない。
【００５５】
図１は、本発明の一実施形態であるガス分析装置の一例を示す概略構成図であり、図２は、図１に示すガス分析装置のブロック図である。なお、地面に水平な一方向をＸ方向とし、地面に水平でＸ方向と垂直な方向をＹ方向とし、Ｘ方向とＹ方向とに垂直な方向をＺ方向とする。また、上述した従来のガス分析装置１０１と同様のものについては、同じ符号を付している。
ガス分析装置１は、光源部１０と、分割測定光と参照光とに分割する測定光分割部１５と、分割測定光受光部２０と、参照光受光部２１と、圧力値Ｐｔｏｔａｌを測定する圧力センサ３１と、ガス温度値Ｔを測定するガス温度センサ３２と、光源部１０を制御するレーザ制御部５０と、マイコンやＰＣで構成される制御部６０とを備える。
【００５６】
光源部１０は、例えば光通信用分布帰還系形（ＤＦＢ）半導体レーザダイオードである。レーザダイオードは、近赤外領域の所定波長範囲（位相速度定数Ａ／νｅｎｄ〜位相速度定数Ａ／νｓｔａｒｔ）内においてレーザの発振波長を調整することが可能となっており、レーザ制御部５０からの制御信号によって制御されるようになっている。なお、光源部１０は、水分子（特定ガス）に吸収される（Ａ／ν２ｎｅｗ）以上（Ａ／ν１ｎｅｗ）以下の波長を含むレーザ光（測定光）を照射することができればよく、周波数ν１ｎｅｗから周波数ν２ｎｅｗまでの範囲を高い分解能で測定する点から、下記式（１６）の関係でｍが２に近くなる（Ａ／ν２ｎｅｗ’）以上（Ａ／ν１ｎｅｗ’）以下の波長のレーザ光を照射することが好ましい。
【００５７】
【数１３】
2015055486.tif 000015
【００５８】
測定光分割部１５は、例えばレーザ光を２方向へ分割するハーフミラーであり、レーザダイオードから出射されたレーザ光の一部を透過することでサンプルガスに照射する分割測定光と、レーザ光の残りを反射することでサンプルガスに照射しない参照光とに分割する。これにより、レーザダイオードの発振波長が吸収ピークの周波数上限νｍａｘから周波数下限νｍｉｎまでの範囲をカバーできないときにも、Ｉ０（ν）を正確に取得できることになる。
【００５９】
分割測定光受光部２０は、光強度を電気信号に変換できるものであればよく、例えばフォトダイオードが用いられる。そして、フォトダイオードは、出射用光学窓７２からサンプル流路７０外にＸ方向で出射されたレーザ光（分割測定光）を受光するように配置されており、サンプルガスを通過した強度Ｉ（ν）のレーザ光（分割測定光）を測定点数（測定ｃｈ）ｎで受光する。
参照光受光部２１は、光強度を電気信号に変換できるものであればよく、例えばフォトダイオードが用いられる。そして、フォトダイオードは、測定光分割部１５でＺ方向に反射されたレーザ光（参照光）を受光するように配置されており、サンプルガスを通過しない強度Ｉ０（ν）のレーザ光（参照光）を測定点数（測定ｃｈ）ｎで受光する。つまり、水分子の吸収を受けなかった場合の測定光の強度変化Ｉ０（ν）を取得する。
【００６０】
制御部６０は、ＣＰＵ６１とメモリ６２と表示部６３と入力装置６４とを備える。また、ＣＰＵ６１が処理する機能をブロック化して説明すると、近似式［ｌｎ（Ｉ０（ν）／Ｉ（ν））］を作成してサンプルガス中の分圧値ＰＨ２Ｏを算出する演算部６１ａを有する。さらに、メモリ６２には、Ｌｏｒｅｎｔｚ型の吸収波形プロファイル関数Φ（ν）やＧａｕｓｓ型の吸収波形プロファイル関数Φ（ν）等を予め記憶するプロファイル関数記憶領域６２ａと、作成された近似式［ｌｎ（Ｉ０（ν）／Ｉ（ν））］におけるａ、ｂ、ｃを記憶するための近似式係数記憶領域６２ｂとを有する。
【００６１】
演算部６１ａは、近似式［ｌｎ（Ｉ０（ν）／Ｉ（ν））］を作成する近似式作成部を有する。近似式作成部は、波長に対する分割測定光の強度変化Ｉ（ν）、波長に対する参照光の強度変化Ｉ０（ν）をそれぞれＡ／Ｄ変換部１、４によってデジタル値に変換し、そのデジタル値から、ｌｎ（Ｉ０（ν）／Ｉ（ν））を作成して、ｌｎ（Ｉ０（ν）／Ｉ（ν））から最小二乗法を用いて近似式［ｌｎ（Ｉ０（ν）／Ｉ（ν））］を作成する制御を行う。
【００６２】
具体的には、レーザ制御部５０は、光源部１０のレーザダイオードからＡ／νｅｎｄの波長からＡ／νｓｔａｒｔの波長までのレーザ光を順にｎ点で出射させ、近似式作成部は、分割測定光受光部２０で受光された分割測定光の強度変化Ｉ（ν）と、参照光受光部２１で受光された参照光の強度変化Ｉ０（ν）とを順にｎ点で取得していく。このとき、式（１６）の関係でｍが２に近くなる（Ａ／ν２ｎｅｗ’）以上（Ａ／ν１ｎｅｗ’）以下の波長のレーザ光を照射させるように入力装置６４で設定することが好ましい。
【００６３】
そして、近似式作成部は、各ｎ点でｌｎ（Ｉ０（ν）／Ｉ（ν））を算出し、ｌｎ（Ｉ０（ν）／Ｉ（ν））の逆数ｌｎ（Ｉ（ν）／Ｉ０（ν））から最小二乗法を用いて式（１３）で示す近似２次多項式を作成する。例えば、ａ＝０．００６３、ｂ＝−３．１８６７、ｃ＝４０７．６０６９と算出する。さらに、式（１４）で示す近似式［ｌｎ（Ｉ０（ν）／Ｉ（ν））］を作成する。そして、ａ＝０．００６３、ｂ＝−３．１８６７、ｃ＝４０７．６０６９を近似式係数記憶領域６２ｂに記憶させる。
【００６４】
演算部６１ａは、圧力値Ｐｔｏｔａｌ、ガス温度値ＴをそれぞれＡ／Ｄ変換部２、３によってデジタル値に変換し、［ｌｎ（Ｉ０（ν）／Ｉ（ν））］が最大値となる周波数νをν０ｎｅｗとしたときに、式（１５）の関係が成立する周波数ν１ｎｅｗと周波数ν２ｎｅｗとを求め、［ｌｎ（Ｉ０（ν０ｎｅｗ）／Ｉ（ν０ｎｅｗ））］と周波数ν１ｎｅｗと周波数ν２ｎｅｗとデジタル値とに基づいて、サンプルガスの圧力値Ｐｔｏｔａｌが大気圧（１０１．３ｋＰａ）以上で行われる計測には、式（９）と式（１１）とに当てはめて数密度ｃを得て、数密度ｃを式（７）に当てはめて分圧値ＰＨ２Ｏを得る。また、サンプルガスの圧力値Ｐｔｏｔａｌが大気圧（１３３Ｐａ）未満で行われる計測には、式（１１）と式（１２）とに当てはめて数密度ｃを得て、数密度ｃを式（７）に当てはめて分圧値ＰＨ２Ｏを得る制御を行う。
【００６５】
具体的には、近似式［ｌｎ（Ｉ０（ν）／Ｉ（ν））］が最大値となるｎの周波数νをν０ｎｅｗとする。次に、演算部６１ａは、式（１５）に基づいて周波数ν１ｎｅｗと周波数ν２ｎｅｗとを求める。そして、演算部６１ａは、圧力値Ｐｔｏｔａｌが大気圧（１０１．３ｋＰａ）以上で行われる計測には、ガス温度センサ３２で測定されたガス温度値Ｔを用いて式（９）と式（１１）とに当てはめて数密度ｃを得て、数密度ｃを式（７）に当てはめて分圧値ＰＨ２Ｏを得る。このとき、付加情報として近似式係数記憶領域６２ｂに記憶されているａ＝０．００６３、ｂ＝−３．１８６７、ｃ＝４０７．６０６９を表示してもよい。一方、圧力値Ｐｔｏｔａｌが大気圧（１３３Ｐａ）未満で行われる計測には、ガス温度センサ３２で測定されたガス温度値Ｔを用いて式（１１）と式（１２）とに当てはめて数密度ｃを得て、数密度ｃを式（７）に当てはめて分圧値ＰＨ２Ｏを得る。このとき、付加情報として近似式係数記憶領域６２ｂに記憶されているａ＝０．００６３、ｂ＝−３．１８６７、ｃ＝４０７．６０６９を表示してもよい。
【００６６】
以上のように、本発明のガス分析装置１によれば、レーザダイオードのもつ波長制御範囲にとらわれずに、サンプルガスの圧力値Ｐｔｏｔａｌが高圧であるときやサンプルガスの濃度が高濃度であるときにも、サンプルガス中の分圧値ＰＨ２Ｏ（特定ガス量情報）を測定することができる。このとき、近似式［ｌｎ（Ｉ０（ν）／Ｉ（ν））］を作成しているため、電気的なノイズ等の影響を低減でき、適切な周波数ν０ｎｅｗと周波数ν１ｎｅｗと周波数ν２ｎｅｗとを求めているので、サンプルガス中の分圧値ＰＨ２Ｏを正確に測定することができる。また、グラフの面積を求める演算よりもかなり少ない負荷で分圧値ＰＨ２Ｏが得られる。
【００６７】
＜他の実施形態＞
上述したガス分析装置１においては圧力センサ３１を備える構成としたが、圧力値Ｐｔｏｔａｌが既知である場合や、測定に影響を及ぼすような変化が想定されない場合には、式（９）と式（１１）と式（１２）とには圧力値Ｐｔｏｔａｌが必要なく、大型化する原因となる圧力センサ３１を備えないような構成としてもよく、また、ガス温度値Ｔが既知である場合や、測定に影響を及ぼすような変化が想定されない場合には、ガス温度センサ３２を備えないような構成としてもよい。
【産業上の利用可能性】
【００６８】
本発明は、レーザ吸収分光法を利用して気体中の水蒸気量情報を測定するガス分析装置等に利用することができる。
【符号の説明】
【００６９】
１  ガス分析装置
１０  光源部
１５  測定光分割部
２０  分割測定光受光部
２１  参照光受光部
６１ａ  演算部
７０  サンプル流路（ガスセル）
７１  入射用光学窓
７２  出射用光学窓
【図１】
2015055486.tif 000018
【図２】
2015055486.tif 000019
【図３】
2015055486.tif 000020
【図４】
2015055486.tif 000021
【図５】
2015055486.tif 000022
【図６】
2015055486.tif 000023
【図７】
2015055486.tif 000024
【図８】
2015055486.tif 000025
【図９】
2015055486.tif 000026