(57)【要約】  本発明の一態様のガス測定装置は、レーザ光照射部（１）によるレーザ光を光共振器（４）に入射させ、レーザ光の波長を目的成分の吸収ピークを含む所定の波長範囲で走査しつつ光検出部（５）によりＣＲＤＳ法による測定を実施するように各部を制御する測定制御部（６）と、所定の波長範囲内の各波長毎に得られたデータに基づいて、吸収スペクトルを作成するスペクトル作成部（７２）と、吸収スペクトルにおける目的成分の吸収ピークの波形形状を多項式で近似し、該多項式における所定の次数の項の係数を取得する演算処理部（７３）と、予め作成された、所定の次数の項の係数と吸収強度との対応関係を示す参照情報に基づいて、前記係数から吸収強度を導出する結果導出部（７４）と、を備える。
2020059452.tif 000002
【特許請求の範囲】
【請求項１】
波長可変であるレーザ光照射部と、一対の高反射ミラー及び被測定ガスが収容される測定セルを含み、前記レーザ光照射部から発して該測定セル内に導入されたレーザ光を共振させる光共振器と、該光共振器から取り出されたレーザ光を検出する光検出部と、を具備し、キャビティリングダウン吸収分光法により被測定ガス中の目的成分の濃度を求めるガス測定装置において、
前記レーザ光照射部によるレーザ光を前記光共振器に入射させ、該レーザ光の波長を目的成分の吸収ピークを含む所定の波長範囲で走査しつつ前記光検出部によりキャビティリングダウン吸収分光法による測定を実施するように各部を制御する測定制御部と、
前記測定制御部の制御の下で所定の波長範囲内の各波長毎に得られたデータに基づいて、吸収スペクトルを作成するスペクトル作成部と、
前記吸収スペクトルにおける前記目的成分の吸収ピークの波形形状を多項式で近似し、該多項式における所定の次数の項の係数を取得する演算処理部と、
予め作成された、前記所定の次数の項の係数と吸収強度との対応関係を示す参照情報に基づいて、前記演算処理部で得られた係数から吸収強度を導出する結果導出部と、
を備えるガス測定装置。
【請求項２】
前記多項式は２次多項式であり、前記所定の次数の項は２次の項である、請求項１に記載のガス測定装置。
【請求項３】
前記測定制御部は前記多項式の２次の項の係数を算出するのに必要十分な波長範囲の波長走査を行う、請求項２に記載のガス測定装置。
【請求項４】
前記必要十分な波長範囲は目的成分の吸収ピークのピーク幅の略１／２に相当する波長範囲である、請求項３に記載のガス測定装置。
【請求項５】
キャビティリングダウン吸収分光法により被測定ガス中の目的成分の濃度を求めるガス測定方法であって、
一対の高反射ミラー、及び被測定ガスが収容される測定セル、を含む光共振器に、レーザ光を入射させ、該レーザ光の波長を前記目的成分の吸収ピークを含む所定の波長範囲で走査しつつ検出することにより、キャビティリングダウン吸収分光法による測定を実施する測定工程と、
前記測定工程において所定の波長範囲内の各波長毎に得られたデータに基づいて、吸収スペクトルを作成するスペクトル作成工程と、
前記吸収スペクトルにおける前記目的成分の吸収ピークの波形形状を多項式で近似し、該多項式における所定の次数の項の係数を取得する演算工程と、
予め作成された、前記所定の次数の項の係数と吸収強度との対応関係を示す参照情報に基づいて、前記係数から吸収強度を導出する結果導出工程と、
を実行するガス測定方法。
【請求項６】
前記多項式は２次多項式であり、前記所定の次数の項は２次の項である、請求項５に記載のガス測定方法。
【請求項７】
前記測定制御部は前記多項式の２次の項の係数を算出するのに必要十分な波長範囲の波長走査を行う、請求項６に記載のガス測定方法。
【請求項８】
前記必要十分な波長範囲は目的成分の吸収ピークのピーク幅の略１／２に相当する波長範囲である、請求項７に記載のガス測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、レーザ光に対する吸収を利用して、被測定ガス中の特定成分の濃度を測定するガス測定装置及びガス測定方法に関し、さらに詳しくは、キャビティリングダウン吸収分光法を利用したガス測定装置及びガス測定方法に関する。
【背景技術】
【０００２】
被測定ガス中の特定成分の濃度を測定する手法としてレーザ吸収分光法（Laser Absorption Spectroscopy）が広く利用されている。レーザ吸収分光法には、直接レーザ吸収分光法（Direct Laser Absorption Spectroscopy）、波長変調レーザ吸収分光法（Laser-Wavelength-Modulation Absorption Spectroscopy）、キャビティリングダウン吸収分光法（Cavity Ring-down Absorption Spectroscopy）などの方式がある（特許文献１など参照）。【０００３】
直接レーザ吸収分光法や波長変調レーザ吸収分光法では、一般に、被測定ガスを収容した測定セルにレーザ光を照射し、該測定セル内を通過する過程でガスの成分により吸収を受けたレーザ光の強度を光検出器で検出する。検出感度を高めるには光吸収のための実効光路長が長いことが望ましいが、これらレーザ吸収分光法では実効光路長は主として測定セルのサイズに依存する。これに対し、キャビティリングダウン吸収分光法（以下、慣用に従って「ＣＲＤＳ」と称す）では、光共振器を用いて光吸収のための実効光路長を長くすることにより、検出感度を大幅に改善することができる（非特許文献１など参照）。【０００４】
ＣＲＤＳ装置における基本的な測定原理について説明する。図４は、一般的なＣＲＤＳ装置における光路を中心とする概略構成図である。
図４において、レーザ光源部１から射出された所定波長のレーザ光は、光スイッチ３を通して、被測定ガスが収容されている測定セル４０に導入される。筒状である測定セル４０の両端には一対の高反射率（ごく僅かに光が透過する）のミラー４７、４８が対向して配置されており、測定セル４０、及びミラー４７、４８は光共振器４を構成する。この光共振器４は例えばレーザ装置等で一般的に用いられているものと同様のファブリペロー共振器であり、共振し得る光の波長（周波数）は共振条件に応じて決まっている。なお、光共振器４は、２枚のミラーを対向して配置した構成の共振器でなく、３枚以上のミラーで構成されるリング型の共振器であってもよい。【０００５】
光共振器４において共振し得る周波数は一般にモード周波数と呼ばれる。図５に示すように、モード周波数は所定の周波数間隔で存在し、光共振器４に導入されたレーザ光の周波数がこのモード周波数と一致しない場合には、該光共振器４内に光のパワーは蓄積されない。一方、レーザ光源部１でのレーザ光の発振周波数がモード周波数と一致するように調整されると、光共振器４内に光のパワーが蓄積される。【０００６】
ＣＲＤＳ装置では、光のパワーが光共振器４内に十分に蓄積されたあと、該光共振器４へ入射するレーザ光を光スイッチ３によって急峻に遮断する。すると、その直前に光共振器４内に蓄積されていた光は一対のミラー４７、４８の間を多数回（実際には数千〜数万回）往復し、その間、測定セル４０内に封入されている被測定ガス中の成分による吸収によって光は徐々に減衰していく。その際に、光共振器４の一方のミラー４８を経て外部へと漏れ出る一部の光の減衰の状態を光検出器５によって繰り返し検出する。この光検出器５により検出したデータに基づいて光の減衰の時定数（リングダウン時間）を求めることで、そのときのレーザ光の周波数における被測定ガス中の目的成分の吸収係数を算出することができる。そして、その吸収係数から目的成分の絶対濃度を求めることができる。また、レーザ光源部１におけるレーザ光の発振周波数を走査しながら同様の測定を繰り返すことにより、被測定ガス中の目的成分の吸収スペクトルを得ることもできる。【０００７】
被測定ガス中の目的成分の吸収係数αを求めるには、通常、次の(1)式が用いられる（特許文献２等参照）。
α＝１／ｃ｛（１／τ）−（１／τ0）｝      …(1)
ここで、ｃは光速、τは測定セル４０内に被測定ガスが収容されているときのリングダウン時間、τ0は測定セル４０内に被測定ガスが収容されていない（例えば真空状態である）とき或いは被測定ガス中の成分による吸収が全く無視できるときの基準状態のリングダウン時間である。また、目的成分（吸収物質）の吸収係数α、単位体積当たりの目的成分分子の数密度ｎ、及び、目的成分による吸収断面積σの関係は次の(2)式のようになる。
α＝ｎσ      …(2)【０００８】
したがって、(1)式、(2)式を用いて、リングダウン時間τ、τ0から、吸収断面積が既知である成分についての絶対濃度を計算することができる。ＣＲＤＳ装置では、光共振器４を用いることで、光が被測定ガス中を透過する距離を伸ばしているため、リングダウン時間τ、τ0の差が大きくなる。それによって、微量な目的成分によるごく僅かな光吸収も検出することができ、他の方式のレーザ吸収分光法に比べて高い検出感度を実現することができる。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】国際公開特許第２０１４／１０６９４０号パンフレット
【特許文献２】特開２０１１−１１９５４１号公報
【非特許文献】
【００１０】
【非特許文献１】橋口幸治、「ガス中微量水分の高効率な計測技術に関する調査研究」、産業技術総合研究所、産総研計量標準報告、2015年10月、Vol.9、No.2、pp.185-205
【発明の概要】
【発明が解決しようとする課題】
【００１１】
上述したようにＣＲＤＳ装置では、非常に高い感度で以て被測定ガス中の濃度の吸収係数を求めることができ、ガス濃度の超高感度測定にしばしば用いられる。ガス濃度を定量的に算出する際には、吸収スペクトルにおいてベースラインから吸収ピークのトップまでの高さが必要である。しかしながら、ＣＲＤＳ装置では装置の状態によって測定値にドリフトが生じ、これがリングダウン時間を測定する度にリングダウン時間が僅かに変化する原因となることがある。測定値のドリフトには主として次のような要因がある。
（１）光共振器４を構成するミラー４７、４８に被測定ガス中の物質が付着することによる、該ミラー４７、４８の実効的な反射率の低下。
（２）上記ミラー４７、４８の微小な変動や該ミラー４７、４８への入射光位置の微小な変動による該ミラー４７、４８の実効的な反射率の変化、及び光共振器長（つまりは片道の光路長）の実効的な変化。
（３）周囲温度の変化によって生じる熱膨張・収縮による光共振器長の変化。【００１２】
図６は、上記のような要因によって、リングダウン時間τ0を測定した時点からリングダウン時間τを測定する時点までの間に光共振器４のミラー４７、４８の実効的な反射率がＲからＲ’に低下した場合における吸収スペクトルの変化の一例を示す図である。このように、反射率が低下すると吸収係数が大きくなる方向（図６では上方向）に吸収ピークがシフトしてしまう。そのため、単に吸収ピークの位置（波数）における吸収係数の値を求めるだけでは成分濃度を正確に算出することはできず、吸収ピークの位置におけるベースラインからピークトップまでの高さの値を取得して吸収係数を算出する必要がある。【００１３】
特に、生体に由来する試料についての同位体比測定では、数週間或いはそれ以上の長い期間に亘って継続的にサンプルを採取し、サンプルに含まれる同位体比の長期間に亘る変化を測定することもある。そのため、上述したような測定値のドリフトの影響がより顕著になる可能性がある。【００１４】
測定値のドリフトの影響を軽減するために一般的には、定期的に又は測定の直前若しくは直後に成分濃度が管理された標準ガスを測定し、その測定結果に基づく校正を行うという方法が採られている。しかしながら、被測定ガスのほかに標準ガスを測定しなければならず、測定のスループットが低下するという問題がある。また、被測定ガスの測定時点と標準ガスの測定時点との間に生じるような、比較的短い時間で起こるドリフトには対応できないという問題もある。【００１５】
本発明は上記課題に鑑みて成されたものであり、その目的とするところは、ＣＲＤＳ法を用いたガス測定装置において、測定のスループットを低下させることなく、比較的短い時間で生じる変動も含めた測定値のドリフトの影響を軽減して、高い精度で目的成分の濃度を算出することができるガス測定装置及びガス測定方法を提供することである。
【課題を解決するための手段】
【００１６】
本出願人は特許文献１において、上記直接レーザ吸収分光法においてレーザ光の波長を走査しつつ測定を行うことで得られた吸収スペクトルに対し数学的な演算処理を実行することにより、波長変調レーザ吸収分光法における変調周波数の高調波信号のプロファイルに相当する曲線を求め、この曲線から被測定ガスの濃度や圧力、温度などを算出する方法を提案している。この方法は、レーザ光の波長変調を行わない直接レーザ吸収分光法での測定結果から、レーザ光の波長変調を行う波長変調レーザ吸収分光法に近い、高精度の測定結果を得ることを目的としたものであり、被測定ガスに照射されるレーザ光の強度の変化を考慮した正規化処理を行うことで、その光強度に依存しないロバストなガス計測を可能としている。【００１７】
一方、ＣＲＤＳ法では、その測定原理から、測定精度は被測定ガスに照射されるレーザ光の強度に依存しないものの、上述したように光共振器を構成するミラーの反射率低下等の要因が測定値のドリフトを引き起こす。ＣＲＤＳ法では実際にレーザ光の波長を変調させた測定を行うことは想定されないが、本発明者は、上記のような要因による測定値のドリフトの影響を軽減するために特許文献１に記載の方法をＣＲＤＳ法に適用することに想到し、その有効性を検証するとともに、ＣＲＤＳ法に適用する際の処理内容やその手順の最適化を検討し、本発明を完成させるに至った。【００１８】
即ち、上記課題を解決するために成された本発明に係るガス測定装置の一態様は、波長可変であるレーザ光照射部と、一対の高反射ミラー及び被測定ガスが収容される測定セルを含み、前記レーザ光照射部から発して該測定セル内に導入されたレーザ光を共振させる光共振器と、該光共振器から取り出されたレーザ光を検出する光検出部と、を具備し、キャビティリングダウン吸収分光法により被測定ガス中の目的成分の濃度を求めるガス測定装置において、
前記レーザ光照射部によるレーザ光を前記光共振器に入射させ、該レーザ光の波長を目的成分の吸収ピークを含む所定の波長範囲で走査しつつ前記光検出部によりキャビティリングダウン吸収分光法による測定を実施するように各部を制御する測定制御部と、
前記測定制御部の制御の下で所定の波長範囲内の各波長毎に得られたデータに基づいて、吸収スペクトルを作成するスペクトル作成部と、
前記吸収スペクトルにおける前記目的成分の吸収ピークの波形形状を多項式で近似し、該多項式における所定の次数の項の係数を取得する演算処理部と、
予め作成された、前記所定の次数の項の係数と吸収強度との対応関係を示す参照情報に基づいて、前記演算処理部で得られた係数から吸収強度を導出する結果導出部と、
を備えるものである。【００１９】
また上記課題を解決するために成された本発明に係るガス測定方法の一態様は、キャビティリングダウン吸収分光法により被測定ガス中の目的成分の濃度を求めるガス測定方法であって、
一対の高反射ミラー、及び被測定ガスが収容される測定セル、を含む光共振器に、レーザ光を入射させ、該レーザ光の波長を前記目的成分の吸収ピークを含む所定の波長範囲で走査しつつ検出することにより、キャビティリングダウン吸収分光法による測定を実施する測定工程と、
前記測定工程において所定の波長範囲内の各波長毎に得られたデータに基づいて、吸収スペクトルを作成するスペクトル作成工程と、
前記吸収スペクトルにおける前記目的成分の吸収ピークの波形形状を多項式で近似し、該多項式における所定の次数の項の係数を取得する演算工程と、
予め作成された、前記所定の次数の項の係数と吸収強度との対応関係を示す参照情報に基づいて、前記係数から吸収強度を導出する結果導出工程と、
を実行するものである。【００２０】
本発明の上記態様のガス測定装置において、測定制御部は、レーザ光照射部等の各部を制御することで、被測定ガスに対し、所定の波長範囲内の各波長においてＣＲＤＳ法による測定を実施する。スペクトル作成部は、異なる波長毎の測定により得られたデータに基づいてリングダウン時間を求め、この実測によるリングダウン時間と予め取得されている基準状態のリングダウン時間とに基づいて吸収係数を算出する。そして、波長毎の吸収係数から、目的成分の吸収ピークの波長を含む所定の波長範囲における吸収係数の変化を示す吸収スペクトルを作成する。【００２１】
なお、波長は波数の逆数であるから、ここでいう「波長」を「波数」で置き換え可能であることは当然である。ここで、波長走査の範囲は、濃度を求める対象である目的成分の既知の吸収波長を中心として所定の範囲とすればよく、その所定の範囲は想定される吸収ピークの幅に応じて決められる値とするとよい。【００２２】
演算処理部は、得られた吸収スペクトルにおける目的成分の吸収ピークのピークトップを含む所定の波長範囲における波形形状を、予め決められた次数の多項式で近似する。典型的には、吸収ピークの波形形状、つまりは各波長において測定の結果得られた吸収係数の値、を２次以上の高次の多項式で近似するフィッティング処理を実行する。特許文献１に記載の方法では、直接レーザ吸収分光法で得られた吸収スペクトルに基づいて、波長変調レーザ吸収分光法における変調周波数の高調波信号（２ｆ信号、１ｆ信号）のプロファイルに相当する曲線を作成するために、吸収ピーク全体を含む広い波長範囲に亘り、波長の各点における所定の波長幅の範囲毎にそれぞれ多項式での近似を実施している。これに対し、本発明では、吸収ピークのピークトップを含む所定の波長範囲内の波形を近似できさえすればよいから、多項式での近似は吸収ピークのピークトップ付近の一つの波長範囲に対してのみ実行すればよい。【００２３】
なお、フィッティング処理の際には、例えば最小二乗法など既知の様々な手法を利用することができる。また、通常、多項式としては２次多項式を用いれば十分であるが、さらに高次の多項式を用いてもよい。【００２４】
上述したようなミラーの反射率の変化や光共振器長の変化などの要因により測定値にドリフトが生じた場合、図６に示したように吸収スペクトルが変化するが、その変化は、バックグラウンドによるベースラインの変動が主であり、目的成分による吸収ピークの波形形状自体の変化は殆どない。そのため、上述したように実測結果である吸収ピークの波形に対するフィッティングにより求めた近似多項式の２次又はそれ以上の偶数次の項の係数の値における上記ドリフトの影響は殆どなく、あったとしても無視できる程度である。また上述したように、ＣＲＤＳ法においてリングダウン時間の測定値はレーザ光の強度に依存しないので、近似多項式の各項の係数にはレーザ光の強度の変動や相違の影響も現れず、特許文献１に記載の方法で実施されているような正規化処理は不要である。こうしたことから、近似多項式の偶数次の項の係数は測定値のドリフトの影響を除外した目的成分による吸収係数を反映している。そこで演算処理部は、近似多項式における所定の次数、例えば２次の項の係数を取得する。【００２５】
上記の所定の次数の項の係数は被測定ガスの温度や圧力が一定、つまり同じ条件であれば、目的成分の種類とその濃度に依存し、測定値のドリフト等には依存しないものとみなせる。そのため、目的の成分について、係数と吸収強度との対応関係を示す参照情報を実験結果等に基づいて予め作成しておき、それを例えばメモリなどに格納しておく。結果導出部はその参照情報を利用して、演算処理部により得られた係数から吸収強度を導出する。ここで吸収強度は吸収係数、又は吸収係数と吸収断面積から求まる絶対濃度のいずれでもよい。即ち、結果導出部では、参照情報をルックアップテーブルとして係数値から即座に吸収係数や絶対濃度を求めることができる。【００２６】
上述したように、本発明では、上記多項式は２次多項式であり、上記所定の次数の項は２次の項とするとよい。【００２７】
この場合、上記測定制御部は、上記多項式の２次の項の係数を算出するのに必要十分な波長範囲の波長走査を行うようにすればよい。【００２８】
特許文献１に記載の方法では吸収ピークの全体を含み、さらにその外側の広い波長範囲までの吸収スペクトルを測定によって取得する必要がある。これに対し本発明では、上述したように、吸収ピークのピークトップ付近の波形を良好に近似できさえすればよいから、一般的には、上記必要十分な波長範囲は、目的成分の吸収ピークの波長を中心としてピーク幅の略１／２に相当する波長範囲とすればよい。【００２９】
これにより、本発明において、測定が必要である波長の点数は例えば吸収ピーク全体をカバーするように測定する場合に比べて、かなり少なくて済む。それによって、測定に要する時間を短縮して測定のスループットの向上を図ることができる。
【発明の効果】
【００３０】
本発明によれば、被測定ガスのほかに標準ガスを測定する手間が不要であるので、測定値のドリフトの影響を軽減する際にも測定のスループットの低下を回避することができる。また、長期間における測定値のドリフトだけでなく比較的短い時間での測定値の変動の影響も高い精度で除去することができるので、目的成分の濃度を高い精度で算出することができる。
【図面の簡単な説明】
【００３１】
【図１】本発明の一実施形態であるＣＲＤＳ装置の要部の構成図。
【図２】本実施形態のＣＲＤＳ装置において目的成分の濃度を求める際の測定及び処理の手順の一例を示すフローチャート。
【図３】本実施形態のＣＲＤＳ装置における目的成分の濃度算出方法の原理を説明するための図。
【図４】一般的なＣＲＤＳ装置の概略構成図。
【図５】光共振器でのモード周波数とレーザ光の発振周波数との関係を示す概略図。
【図６】リングダウン時間τ0を測定した時点からリングダウン時間τを測定する時点までの間に、光共振器のミラーの実効的な反射率が低下した場合における、吸収スペクトルの変化の一例を示す図。
【発明を実施するための形態】
【００３２】
以下、本発明の一実施形態であるガス測定装置及びガス測定方法について、添付図面を参照して説明する。図１は、本実施形態のガス測定装置の概略構成図である。
【００３３】
本実施形態のガス測定装置において、測定系の構成は図４に示した一般的なＣＲＤＳ装置と同じであり、レーザ光源部１、レーザ駆動部２、光スイッチ３、光共振器４、及び光検出器５、を備える。光共振器４は、被測定ガスであるサンプルガスが収容される略円筒状の測定セル４０と、該測定セル４０の両端に対向して配置された一対の高反射率のミラー４７、４８と、を含む。測定セル４０にはガス導入管４１とガス排出管４３とが接続され、ガス導入管４１には導入バルブ４２が設けられ、ガス排出管４３には排出バルブ４４が設けられている。
【００３４】
入力部８が接続された測定制御部６は、後述する測定やデータ処理を実行するためにレーザ駆動部２等の各部を制御するものである。測定制御部６に含まれる図示しないメモリには、測定対象である成分の種類などに対応する吸収ピーク位置（波長）や波長走査範囲などの情報が予め格納されている。また、光検出器５による検出信号が入力されるデータ処理部７は、機能的なブロックとして、リングダウン時間算出部７１、スペクトル作成部７２、フィッティング処理部７３、濃度算出部７４、及び、参照情報記憶部７５などを備える。また、データ処理部７に接続された出力部９は例えば表示モニタなどである。
【００３５】
本実施形態のガス測定装置において、被測定ガス中の目的成分の濃度を求める際の測定動作及び処理動作、つまり本装置で実施されるガス測定方法を、図１に加え図２を参照して説明する。図２は、このガス測定装置を用いて目的成分の濃度を求める際の、測定及び処理の手順の一例を示すフローチャートである。
ユーザは事前に、目的成分の種類などを入力部８から入力する。
【００３６】
測定制御部６は、排出バルブ４４を閉鎖した状態で導入バルブ４２を開放し、測定セル４０内に被測定ガスを導入する。図示しない圧力センサにより検出される圧力が所定値になったならば、導入バルブ４２を閉鎖して測定セル４０内に被測定ガスを充満させる。そして、測定制御部６は、事前に指定された目的成分に対応した波長範囲の情報を取得し、レーザ駆動部２を通してレーザ光源部１で発生するレーザ光の波長をその波長範囲内で順次走査しつつ、各波長においてＣＲＤＳ法による測定を実施してリングダウン時間を測定する（ステップＳ１：測定工程）。
【００３７】
即ち、所定の波長範囲内の各波長において、レーザ光源部１は測定セル４０中の被測定ガスにレーザ光を照射し、光スイッチ３は所定のタイミングでレーザ光を遮断する。リングダウン時間算出部７１は、レーザ光を遮断する直前から所定の時間が経過するまでの間、光検出器５により検出されたデータを収集する。そして、そのデータに基づいて、波長毎にリングダウン時間τを算出する。
【００３８】
次に、スペクトル作成部７２は、波長毎に、実測データに基づくリングダウン時間τと、参照情報記憶部７５に保存されている基準状態のリングダウン時間τ0とに基づいて、吸収係数を計算する。この吸収係数の計算方法は従来と同じであり、例えば上記(1)、(2)式を用いればよい。そして、波長毎に算出した吸収係数の値を集めることで、所定の波長範囲における吸収スペクトルを求める（ステップＳ２：スペクトル作成工程）。
【００３９】
このときに得られる吸収スペクトルには、上述した様々な要因による測定値のドリフトの影響が含まれる可能性がある。そこで、次のような手順で、こうしたドリフトを含む測定値の変動を除去して、正確な濃度を算出する。まず、測定値の変動の影響を除去する原理を説明する。
【００４０】
いま、基準状態のリングダウン時間τ0を測定したときと、被測定ガスに対するリングダウン時間τを測定したときとで、光共振器４のミラーの反射率が相違する（つまりはドリフトしている）場合を想定する。リングダウン時間τ0の測定時におけるミラーの反射率をＲ、リングダウン時間τの測定時におけるミラーの反射率をＲ’とすると、測定結果から算出されるＣＲＤ信号Ｓ(ν)は次の(3)式となる。
Ｓ(ν)＝（１／ｃ）｛（１／τ）−（１／τ0）｝＝（１／ｃ）｛（ｃ［（１−Ｒ’）＋α(ν)Ｌ］／Ｌ）−（ｃ［１−Ｒ］／Ｌ）｝＝α(ν)＋（Ｒ−Ｒ’）／Ｌ      …(3)
ここで、α(ν)はミラーの反射率に変動がない場合の吸収係数である。
【００４１】
上記ＣＲＤ信号Ｓ(ν)を、波数軸を横軸として描いたものが吸収スペクトルである。この吸収スペクトルに対し、特許文献１に記載の方法の適用を検討する。いま、図３（ａ）に示す吸収スペクトル（但し、ここでは横軸は中心波数との波数差Δν軸である）に対し、特許文献１に記載の多項式近似の処理を適用して、波長変調レーザ吸収分光法における変調周波数の２ｆ信号（２次微分）、１ｆ信号（１次微分）、定数項（０次微分）のプロファイルに相当する曲線を求めると、図３（ｂ）、（ｃ）、（ｄ）に示すようになる。
【００４２】
但し、特許文献１に記載の方法では、吸収スペクトルにおいて、波長変調幅に相当する波長幅の範囲内毎にそれぞれ多項式近似を行っているが、本実施形態における方法では、図３（ｂ）等に示されるような曲線自体を求める必要はなく、求める必要があるのはΔν＝０における２ｆ信号の値のみである。その理由は次のとおりである。
【００４３】
上記(3)式におけるＳ(ν)とα(ν)とがそれぞれ多項式で近似が可能であるとすれば、それぞれ次の式で表すことができる。
Ｓ(ν)＝ｂ0'＋ｂ1'(ν−<ν>)＋ｂ2'（ν−<ν>）2＋ｂ3'（ν−<ν>）3＋・・・      …(4)
α(ν)＝ｂ0＋ｂ1（ν−<ν>）＋ｂ2（ν−<ν>）2＋ｂ3（ν−<ν>）3＋・・・      …(5)
ここで、<ν>は吸収スペクトルを取得する波長範囲の中心の波長、つまり、本実施形態ではΔν＝０の位置である。
【００４４】
(4)式と(5)式の各項の係数を比較すると、以下のようになる。
ｂ0'≒ｂ0＋（Ｒ−Ｒ’）／Ｌ      …(6)
ｂ1'≒ｂ1      …(7)
ｂ2'≒ｂ2      …(8)
【００４５】
(8)式は、多項式の２次の項の係数ｂ2がミラーの反射率の変動に依存しないことを意味している。なお、(7)式によれば、数式上では１次の項の係数ｂ1もミラーの反射率の変動に依存しないことになる。しかしながら、図３（ｃ）を見れば明らかであるように、１次を含め奇数次の項の値は、目的成分の吸収ピークのピークトップの波長において０になる。そのため、目的成分による吸収強度を求める際には、近似多項式の奇数次の項を利用することはできず、２次以上の偶数次の項の係数のみが利用可能である。
【００４６】
また、特許文献１に記載の方法はレーザ直接吸収分光法により取得した吸収スペクトルを処理対象としており、レーザ直接吸収分光法では、被測定ガスに照射されるレーザ光の強度が変動するとそれが吸収スペクトルに現れる。これを回避するために、多項式近似によって得られたｂ2信号をｂ1信号やｂ0信号で除する正規化処理を実行することで、光強度の変動に依存せず透過特性にのみ依存する値を得るようにしている。これに対し、ＣＲＤＳ法では原理上、被測定ガスに入射する光の光強度の変動は測定値に影響しない。そのため、上記のような正規化処理を行うことなく、多項式近似により求まる２次項の係数をそのまま用いることができる。
【００４７】
上記説明では、基準状態のリングダウン時間τ0測定時と被測定ガスに対するリングダウン時間τ測定時とでミラーの反射率が相違する場合を想定していたが、上述したように、実効的な光共振器長が変化する場合や、光共振器長とミラーの反射率との両方が変化する場合でも同様に、吸収スペクトルについての近似多項式における２次以上の偶数次の項の係数を利用して、上記変動要因の影響の少ない信号を得ることができる。
【００４８】
いま、共振器長にドリフトがあり、リングダウン時間τ0測定時における光共振器長をＬ、リングダウン時間τ測定時における光共振器長をＬ’とすると、測定結果から算出されるＣＲＤ信号Ｓ(ν)は次の(9)式となる。
Ｓ(ν)＝（１／ｃ）｛（１／τ）−（１／τ0）｝＝（１／ｃ）｛（ｃ［（１−Ｒ）＋α(ν)Ｌ’］／Ｌ’）−（ｃ［１−Ｒ］／Ｌ）｝＝α(ν)＋｛（１−Ｒ）（Ｌ−Ｌ’）／Ｌ・Ｌ’｝      …(9)
【００４９】
また、光共振器４内のミラー４７、４８の反射率と共振器長との両方にドリフトがある場合、測定結果から算出されるＣＲＤ信号Ｓ(ν)は次の(10)式となる。
Ｓ(ν)＝（１／ｃ）｛（１／τ）−（１／τ0）｝＝（１／ｃ）｛（ｃ［（１−Ｒ’）＋α(ν)Ｌ’］／Ｌ’）−（ｃ［１−Ｒ］／Ｌ）｝＝α(ν)＋｛［（１−Ｒ’）Ｌ−（１−Ｒ）Ｌ’］／Ｌ・Ｌ’｝      …(10)
これらの式は(3)式と同様の式であるから、この場合にも、Ｓ(ν)の近似多項式の２次項の係数ｂ2がそれらの変動に依存しないことを意味していることが分かる。
【００５０】
なお、ミラーの反射率及び光共振器長のドリフトにレーザ光の波長依存性があるとすると、測定により得られるＣＲＤ信号Ｓ(ν)から算出した２ｆ信号にも誤差が生じ、(8)式が成り立たなくなるおそれがある。しかしながら、共振器長の変動についてはレーザ光の波長に依存する要因はないので、その波長依存性は無視できる。一方、ミラーに被測定ガス中の物質が付着することによる反射率の低下に関しては、波長依存性が考えられるものの、ここで使用される１cm-1以下の波数範囲では波長依存性は無視できると考えられる。したがって、ここでは上記変動要因の影響の少ない信号を得るために、吸収スペクトルについての近似多項式における２次以上の偶数次の項の係数を利用すればよい。
【００５１】
図２に示すフローチャートに戻り、処理手順を説明する。本実施形態のガス測定装置では、フィッティング処理部７３が、吸収スペクトル上の吸収ピークの波形に対して２次多項式でフィッティングするように該多項式の各項の係数を決める（ステップＳ３：演算工程）。フィッティング処理の際には、例えば最小二乗法など、既知の手法を利用することができる。吸収ピークのピークトップ付近の形状を適切に近似できればよいので、多項式としては２次の多項式を用いればよく、吸収ピークのピーク幅の１／２程度の波長範囲のみについてフィッティングを行えばよい。したがって、図３（ａ）に示した吸収ピークでは、図中に点線で囲んだ範囲のピーク波形についてフィッティングを行えばよい。つまり、ステップＳ１における測定もその波長範囲のみについて実施すれば十分である。そして、フィッティング処理部７３は、近似多項式の中の２次項の係数を取得する（ステップＳ４：演算工程）。
【００５２】
上述したように、２次項の係数は吸収強度（吸収係数）に直接対応している。そこで、目的成分（及び測定対象となる可能性のある他の成分）についての２次項の係数と吸収強度との関係を予備実験等により求めておき、この関係を例えばテーブル形式で参照情報記憶部７５に格納しておく。なお、この参照情報記憶部７５に格納される情報はユーザ自身が作成することもできるが、本装置のメーカが作成しておくものとすることもできる。また、吸収強度は被測定ガスの温度や圧力に依存するので、参照情報記憶部７５に格納される情報は、所定の温度及び所定の圧力の下での２次項の係数と吸収強度との関係であり、同じ所定の温度、所定の圧力の下で被測定ガスの測定を実施するものとする。
【００５３】
濃度算出部７４は、実測に基づいて得られた２次項の係数を、参照情報記憶部７５に格納されている情報に照合して、対応する吸収強度を取得する（ステップＳ５：結果導出工程）。そして、その吸収強度から被測定ガス中の目的成分の絶対濃度を算出し（ステップＳ６）、その結果を出力部９に出力して表示する。
【００５４】
以上のようにして本実施形態のガス測定装置では、測定値の短期的な変動を含めたドリフトの影響を軽減して、精度の高い濃度を算出することができる。また、本実施形態のガス測定装置では、吸収ピークのピークトップを中心とする比較的狭い波長範囲についてのみ測定を実施すればよいため、測定時間を短縮することができる。測定時間を短縮した分だけ測定のスループットを向上させる、或いは、その代わりに同じ被測定ガスに対する測定の繰り返し回数を増やし測定結果を積算することで、測定精度を向上させることも可能である。
【００５５】
なお、上記実施形態は本発明の一例であり、本発明の趣旨の範囲で適宜に変形や修正、追加などを行っても、本願特許請求の範囲に包含されることは明らかである。
【符号の説明】
【００５６】
１…レーザ光源部
２…レーザ駆動部
３…光スイッチ
４…光共振器
４０…測定セル
４１…ガス導入管
４２…導入バルブ
４３…ガス排出管
４４…排出バルブ
４７、４８…ミラー
５…光検出器
６…測定制御部
７…データ処理部
７１…リングダウン時間算出部
７２…スペクトル作成部
７３…フィッティング処理部
７４…濃度算出部
７５…参照情報記憶部
８…入力部
９…出力部
【図１】
2020059452.tif 000003
【図２】
2020059452.tif 000004
【図３】
2020059452.tif 000005
【図４】
2020059452.tif 000006
【図５】
2020059452.tif 000007
【図６】
2020059452.tif 000008
【国際調査報告】
2020059452.tif 000009
2020059452.tif 000010
2020059452.tif 000011
2020059452.tif 000012