(57)【要約】  測定セル（４０）中に被測定ガスを所定圧力となるように導入したあと、レーザ光源部（１）、光スイッチ（３）、光共振器（４）、及び光検出器（５）により、ＣＲＤＳによる所定の波数における測定を行う。次に、測定セル（４０）内の被測定ガスの一部を排出して圧力を下げ、目的成分である14CO2の吸収ピークの波数における測定を行う。高圧力の下での測定時には14CO2による吸収の影響が無視できるため、リングダウン時間から12CO2等のバックグラウンドの濃度が求まる。低圧力の下での測定データから求まるリングダウン時間から算出される吸収係数にはバックグラウンドの影響が含まれるが、高圧力の下で求めたバックグラウンドの濃度から低圧力の下でのバックグラウンドの吸収係数を計算することができる。濃度演算部（７３）はそれを用いてバックグラウンドの影響を除去した14CO2のみの吸収係数を求め、14CO2のみの濃度を算出する。異なる圧力の下で実施された２回の測定の結果から14CO2等の目的成分の精度の高い絶対濃度を得ることができ、従来より測定時間を短縮することができる。
2020105118.tif 000002
【特許請求の範囲】
【請求項１】
キャビティリングダウン吸収分光法により、被測定ガスに含まれる目的成分の濃度を求めるガス測定方法において、
第１圧力の下で前記被測定ガスに対してレーザ光を照射することにより、前記目的成分の吸収ピークの波長についてキャビティリングダウン吸収分光法による測定を実施する第１測定ステップと、
前記第１圧力と異なる第２圧力の下で前記被測定ガスに対してレーザ光を照射することにより、キャビティリングダウン吸収分光法による測定を実施する第２測定ステップと、
前記第１測定ステップの結果と前記第２測定ステップの結果とに対する演算を行うことにより、前記目的成分の濃度を算出する演算ステップと、
を有する、ガス測定方法。
【請求項２】
請求項１に記載のガス測定方法であって、
前記第２測定ステップでは、前記目的成分の吸収ピークの波長とは異なる、該目的成分による吸収の影響を無視できる波長についてキャビティリングダウン吸収分光法による測定を実施し、
前記演算ステップでは、前記第２測定ステップの結果に基づいて前記第２圧力の下での前記被測定ガス中の前記目的成分以外の成分の濃度を推定し、該濃度から、前記第１測定ステップの結果から求まる吸収係数における前記目的成分以外の成分の吸収の寄与を推定し、その影響を除去する演算を行うものである、ガス測定方法。
【請求項３】
請求項１に記載のガス測定方法であって、
前記第２測定ステップでは、前記目的成分の吸収ピークの波長についてキャビティリングダウン吸収分光法による測定を実施し、
前記演算ステップでは、前記第１測定ステップの結果と前記第２測定ステップの結果とに基づく連立方程式を作成し、該連立方程式を解くことで、前記被測定ガス中の前記目的成分の濃度以外の成分による吸収の影響を除去又は軽減した前記目的成分の濃度を算出するものである、ガス測定方法。
【請求項４】
請求項１に記載のガス測定方法であって、
前記被測定ガスはＣＯ2ガスを含み、前記目的成分は該ＣＯ2中の同位体の一つである14ＣＯ2である、ガス測定方法。
【請求項５】
キャビティリングダウン吸収分光法により、被測定ガス中の目的成分の濃度を求めるガス測定装置において、
レーザ光照射部と、
被測定ガスが収容される測定セルを含み、前記レーザ光照射部から発して該測定セル内に導入されたレーザ光を共振させる光共振器と、
該光共振器から取り出されたレーザ光を検出する光検出部と、
前記測定セル中の被測定ガスの圧力を調整する圧力調整部と、
前記測定セル中の被測定ガスに対してキャビティリングダウン吸収分光法による測定を実施するときに前記圧力調整部を制御する制御部と、
該制御部による制御の下で互いに異なる圧力の下で得られた複数の測定結果に対する演算を行うことにより、前記目的成分の濃度を算出する演算処理部と、
を備える、ガス測定装置。
【請求項６】
請求項５に記載のガス測定装置であって、
前記制御部は、
第１圧力の下で前記被測定ガスに対してレーザ光を照射することにより、前記目的成分の吸収ピークの波長についてキャビティリングダウン吸収分光法による測定を実施する第１測定ステップと、
前記第１圧力と異なる第２圧力の下で前記被測定ガスに対してレーザ光を照射することにより、キャビティリングダウン吸収分光法による測定を実施する第２測定ステップと、
を実施するように、前記圧力調整部のほか、前記レーザ光照射部及び前記光検出部を制御する、ガス測定装置。
【請求項７】
請求項６に記載のガス測定装置であって、
前記制御部は、前記第２圧力の下での測定の際に、前記目的成分の吸収ピークの波長とは異なる、該目的成分による吸収の影響を無視できる波長についてキャビティリングダウン吸収分光法による測定を実施し、
前記演算処理部は、前記第２圧力の下での結果に基づいて該第２圧力の下での前記被測定ガス中の前記目的成分以外の成分の濃度を推定し、該濃度から、前記第１圧力の下での結果から求まる吸収係数における前記目的成分以外の成分の吸収の寄与を推定し、その影響を除去する演算を行う、ガス測定装置。
【請求項８】
請求項６に記載のガス測定装置であって、
前記制御部は、前記第２圧力の下での測定の際に、前記目的成分の吸収ピークの波長についてキャビティリングダウン吸収分光法による測定を実施し、
前記演算処理部は、前記第１圧力の下での結果と前記第２圧力の下での結果とに基づく連立方程式を作成し、該連立方程式を解くことで、前記被測定ガス中の前記目的成分の濃度以外の成分による吸収の影響を除去又は軽減した前記目的成分の濃度を算出する、ガス測定装置。
【請求項９】
請求項６に記載のガス測定装置であって、
前記圧力調整部は、前記測定セル中に第２圧力で被測定ガスを封入した状態から、該測定セル中から一部の被測定ガスを外部に強制的に排出することで該測定セル中の被測定ガスの圧力を第１圧力に調整するものである、ガス測定装置。
【請求項１０】
請求項６に記載のガス測定装置であって、
前記圧力調整部は、前記測定セル中に被測定ガスを供給し第１の所定の圧力で被測定ガスを封入した状態から、先に供給されずに残存されていた被測定ガスを該測定セル中に追加供給して封入することで該測定セル中の被測定ガスの圧力を第２の所定の圧力に調整するものである、ガス測定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、レーザ光に対する吸収を利用して、被測定ガス中の特定成分の濃度を測定するガス測定装置及びガス測定方法に関する。
【背景技術】
【０００２】
被測定ガス中の特定成分の濃度を測定する手法として、レーザ吸収分光法が広く利用されている。レーザ吸収分光法には幾つかの手法が知られているが、その一つの手法として、キャビティリングダウン吸収分光法（Cavity Ring-down Absorption Spectroscopy、以下、慣用に従って「ＣＲＤＳ」と称す）がある。ＣＲＤＳは、光共振器を用いて光吸収のための実効光路長を長くすることにより、検出可能な吸光度、つまりはその検出感度を大幅に改善できる手法である（非特許文献１など参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１１−１１９５４１号公報
【特許文献２】特開２０１８−４６５６号公報
【非特許文献】
【０００４】
【非特許文献１】橋口幸治、「ガス中微量水分の高効率な計測技術に関する調査研究」、産業技術総合研究所、産総研計量標準報告、2015年10月、Vol.9、No.2、pp.185-205
【非特許文献２】鈴木彌生子、「キャビティリングダウン分光分析法を用いた果物および野菜の水分の酸素・水素同位体比」、日本水文科学会誌、2016年、第46巻、第2号、pp.157-166
【非特許文献３】パン・ドゥ（Pan Du）ほか２名、「インプルーブド・ピーク・デテクション・イン・マス・スペクトラム・バイ・インコーポレーティング・コンティニュアス・ウェブレット・トランスフォーム-ベースド・パターン・マッチング（Improved peak detection in mass spectrum by incorporating continuous wavelet transform-based pattern matching）」、オックスフォード・ユニバーシティ・プレス（Oxford University Press）、バイオインフォマティックス（Bioimformatics）、2006年、Vol.22、No.17、pp.2059-2065
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明の課題について、図面を参照して具体的に説明する。
図８は、一般的なＣＲＤＳ装置の概略構成図である。図８において、レーザ光源部１から射出された所定波長のレーザ光は、光スイッチ３を通して、被測定ガスが収容されている測定セル４０に導入される。筒状である測定セル４０の両端には一対の高反射率（ごく僅かに光が透過する）のミラー４７、４８が対向して配置されており、測定セル４０、ミラー４７、４８は光共振器４を構成する。この光共振器４は例えばレーザ装置等で一般的に用いられているものと同様のファブリペロー共振器であり、共振し得る光の波長（周波数）は共振条件に応じて決まっている。なお、光共振器４は、２枚のミラーを対向して配置した構成の共振器でなく、３枚以上のミラーで構成されるリング型の共振器であってもよい。【０００６】
光共振器４において共振し得る周波数は一般にモード周波数と呼ばれる。図９に示すように、モード周波数は所定の周波数間隔で存在し、光共振器４に導入されたレーザ光の周波数がこのモード周波数と一致しない場合には、該光共振器４内に光のパワーは蓄積されない。一方、レーザ光源部１でのレーザ光の発振周波数がモード周波数と一致するように調整されると、光共振器４内に光のパワーが蓄積される。【０００７】
ＣＲＤＳ装置では、光のパワーが光共振器４内に十分に蓄積されたあと、該光共振器４へ入射するレーザ光を光スイッチ３によって急峻に遮断する。すると、その直前に光共振器４内に蓄積されていた光は一対のミラー４７、４８の間を多数回（実際には数千〜数万回）往復し、その間、測定セル４０内に封入されている被測定ガス中の成分による吸収によって光は徐々に減衰していく。その際に、光共振器４の一方のミラー４８を経て外部へと漏れ出る一部の光の減衰の状態を光検出器５によって繰り返し検出する。この光検出器５により検出したデータに基づいて光の減衰の時定数（リングダウン時間）を求めることで、そのときのレーザ光の周波数における被測定ガス中の目的成分の吸収係数を算出することができる。そして、その吸収係数から目的成分の絶対濃度を求めることができる。また、レーザ光源部１におけるレーザ光の発振周波数を走査しながら同様の測定を繰り返すことにより、被測定ガス中の目的成分の吸収スペクトルを得ることもできる。【０００８】
被測定ガス中の目的成分の吸収係数αを求めるには、通常、次の(1)式が用いられる（特許文献１等参照）。
α＝１／ｃ｛（１／τ）−（１／τ0）｝      …(1)
ここで、ｃは光速、τは測定セル４０内に被測定ガスが収容されているときのリングダウン時間、τ0は測定セル４０内に被測定ガスが収容されていない（例えば真空状態である）とき或いは被測定ガス中の成分による吸収が全く無視できるときのリングダウン時間である。また、目的成分（吸収物質）の吸収係数α、数密度ｎ、吸収断面積σの関係は次の(2)のようになる。
α＝ｎσ      …(2)【０００９】
したがって、(1)、(2)式を用いてリングダウン時間τ、τ0から、吸収断面積が既知である成分についての絶対濃度を計算することができる。ＣＲＤＳ装置では、光共振器４を用いて光が被測定ガスを透過する実効的な距離を伸ばしているため、リングダウン時間τ、τ0の差が大きくなる。それによって、微量な目的成分によるごく僅かな光吸収も検出することができ、他の方式のレーザ吸収分光法に比べて高い検出感度を実現することができる。【００１０】
上述したようにＣＲＤＳ装置では、非常に高い感度で以て被測定ガス中の成分の濃度を測定することができる。そのため、ＣＲＤＳ装置は、被測定ガス中のＣＯ2やＨ2Ｏの同位体比を高い精度で以て測定する目的でしばしば用いられる（特許文献２等参照）。また、ＣＲＤＳ装置を用いた同位体比測定の応用は、農産物の産地特定等の様々な分野で進んでいる（非特許文献２参照）。【００１１】
しかしながら、例えば炭素の放射性同位体14Ｃを含む14ＣＯ2（天然同位体存在比：１×１０-12）のように、被測定ガス中にごく微量しか含まれない成分の濃度を測定したい場合、その被測定ガスに含まれる別の成分による吸収の影響、つまりはバックグラウンドが無視できない。図１０は、14ＣＯ2吸収線のピーク波長付近における被測定ガスの吸収スペクトルの概略図である。この吸収線ピークのベースラインは実際にはその殆どが、高い濃度で含まれる12ＣＯ2や13ＣＯ2による吸収に由来するもの、つまりはバックグラウンドである。このバックグラウンドを無視してしまうと目的成分の濃度を正確に求めることができない。【００１２】
そこで、高い精度で以て比較的微量である同位体ガスの濃度を測定したい場合、従来、次のようにしてバックグラウンドを除去する作業が行われていた。【００１３】
即ち、レーザ光の発振波長を所定の範囲で変化させながら、測定対象である同位体ガスの吸収ピーク付近の複数の波長でそれぞれＣＲＤＳによる測定を行い、その測定結果から吸収係数をそれぞれ算出する。そして、異なる波長に対してそれぞれ得られた複数の吸収係数に基づき、フォークト（Voigt）関数やローレンツ（Lorentz）関数によるフィッティング処理を行うことで、バックグラウンド（図１０でいえば、14ＣＯ2以外のガス種による吸収）のスペクトル波形を推定する。即ち、図１０ではＡ部やＢ部のスペクトル波形から、ピークのベースラインに相当するＣ部のスペクトル波形を推定する。そして、その推定したスペクトル波形を用いてバックグラウンド除去を行うことで、測定対象である同位体ガスのみの吸収係数を求め、その吸収係数から濃度を算出する。【００１４】
バックグラウンドのスペクトル波形のフィッティングを精度良く行うには、或る程度多くの波長における吸収係数を測定によって求めることが必要である。そのため、一つの被測定ガスに対して多数回の測定を行う必要があり測定時間が長くなる。もちろん、こうしたバックグラウンドのスペクトル波形の推定は、被測定ガスが変わる毎に行う必要がある。【００１５】
測定中には、測定ガス中の一部の成分の吸着によるミラーの実効的な反射率の低下や、ミラーの微小な変動或いは入射光位置の微小な変動による実効的な反射率の変動や共振器長の変動、さらには温度の微小な変動によって生じる光共振器の熱膨張による共振器長等の変動が生じ得る。そのため、測定時間が長引くほど、測定途中で測定状態が変化する可能性が高く又その変化量が大きくなる可能性があり、これにより正確に濃度を求めることができない場合がある。【００１６】
また、以下に説明するように、上述したように推定したベースラインに相当するスペクトル波形が正確でない場合、フィッティングを行ったとしても濃度を正確に求めることができない。図１１（ａ）は、14ＣＯ2の吸収ピーク位置（波数）付近の所定の波数範囲における、炭素の安定同位体12Ｃ、13Ｃをそれぞれ含む12ＣＯ2、13ＣＯ2、及び上記14ＣＯ2についての波数と吸収係数との関係を示す吸収スペクトルを計算により求めたものである。また図１１（ｂ）は、14ＣＯ2の吸収ピーク位置での12ＣＯ2、13ＣＯ2、及び14ＣＯ2による吸収の寄与度合いを計算した結果である。【００１７】
図１１（ａ）から分かるように、この例では、14ＣＯ2の吸収ピークの位置に13ＣＯ2の吸収ピークの裾部が重なってしまっている。一般的には、13ＣＯ2の吸収が単独で存在するスペクトル領域まで測定対象の波数範囲を広げ、14ＣＯ2と13ＣＯ2とのそれぞれに対してスペクトルフィッティングを行い、14ＣＯ2の吸収ピークの位置に重なっている13ＣＯ2の吸収ピークを分離し除去することで14ＣＯ2の単独の吸収ピークを求め、これを評価する。しかしながら、上述したようにスループットを改善するために14ＣＯ2の吸収ピークのみを測定した場合、13ＣＯ2単独の吸収ピークが確認できない。そのため、13ＣＯ2の吸収に対して正しくスペクトルフィッティングが行えず、推定したベースラインに相当するスペクトル波形が正確でなくなってしまう。【００１８】
実際に発振波長を走査しながらＣＲＤＳによる測定を繰り返すことで求まるスペクトルは図１０に実線で示す波形であるが、この波形からベースラインを推定しても、図１１（ａ）に示すような、重なっているピークを反映したベースラインを推定することはできない。そのため、バックグランドを正確に除去することはできず、目的成分の吸収係数を正確に求めることができない。【００１９】
このように吸収ピークの位置が互いに近接した複数の成分が被測定ガス中に存在する場合に、被測定ガスの吸収スペクトルにおいて当該複数の成分に対応する信号同士が重なり合ってしまい、目的成分の吸収係数を正確に求められない場合がある。こうした問題に対し、重なり合った信号同士を画像処理により分離する技術（いわゆる、ピークピッキング）（例えば非特許文献３参照）を適用することも可能であるが、一般的なピークピッキングでは、信号分離の正確性はオペレータの熟練度や技量に依存することが多い。そのため、常に正確な信号分離が行えるとは限らず、測定結果の信頼性や再現性を確保することが難しい。【００２０】
本発明は上述したような課題に鑑みて成されたものであり、その目的は、目的成分の正確な吸収係数や濃度を求めることができるガス測定装置及びガス測定方法を提供することである。
【課題を解決するための手段】
【００２１】
特許文献２等にも記載されているように、ＣＲＤＳにおけるリングダウン特性（光強度の指数関数的な減衰の程度）は、温度や圧力に依存する。そこで本発明者は、圧力によってリングダウン特性、即ち吸収係数が変化することに着目した。何故なら、被測定ガスの温度を変化させると、熱膨張効果によって光共振器の光路長やミラーの反射率が変化してしまい、それによって光共振器のモード周波数やモード線幅が変化してしまうため、安定した測定が困難となる。これに対し、被測定ガスの圧力であれば、比較的容易に且つ正確に変化させることが可能であるからである。そして本発明者はシミュレーション計算などを繰り返し、実用的に変化させることが可能な圧力の範囲で同じ成分による吸収の程度が大幅に変化し得ることを見いだし、本発明を完成させるに至った。【００２２】
即ち、上記課題を解決するために成された本発明に係るガス測定方法は、キャビティリングダウン吸収分光法（ＣＲＤＳ）により、被測定ガスに含まれる目的成分の濃度を求めるガス測定方法において、
第１圧力の下で前記被測定ガスに対してレーザ光を照射することにより、前記目的成分の吸収ピークの波長についてキャビティリングダウン吸収分光法による測定を実施する第１測定ステップと、
前記第１圧力と異なる第２圧力の下で前記被測定ガスに対してレーザ光を照射することにより、キャビティリングダウン吸収分光法による測定を実施する第２測定ステップと、
前記第１測定ステップの結果と前記第２測定ステップの結果とに対する演算を行うことにより、前記目的成分の濃度を算出する演算ステップと、
を有するものである。【００２３】
また上記課題を解決するために成された本発明に係るガス測定装置は、上記本発明に係るガス測定方法を実施するための一つの装置であり、キャビティリングダウン吸収分光法により、被測定ガス中の目的成分の濃度を求めるガス測定装置において、
レーザ光照射部と、
被測定ガスが収容される測定セルを含み、前記レーザ光照射部から発して該測定セル内に導入されたレーザ光を共振させる光共振器と、
該光共振器から取り出されたレーザ光を検出する光検出部と、
前記測定セル中の被測定ガスの圧力を調整する圧力調整部と、
前記測定セル中の被測定ガスに対してキャビティリングダウン吸収分光法による測定を実施するときに前記圧力調整部を制御する制御部と、
該制御部による制御の下で互いに異なる圧力の下で得られた複数の測定結果に対する演算を行うことにより、前記目的成分の濃度を算出する演算処理部と、
を備えるものである。【００２４】
上記本発明に係るガス測定装置において、上記本発明に係るガス測定方法を実施するために、前記制御部は、
第１圧力の下で前記被測定ガスに対してレーザ光を照射することにより、前記目的成分の吸収ピークの波長についてキャビティリングダウン吸収分光法による測定を実施する第１測定ステップと、
前記第１圧力と異なる第２圧力の下で前記被測定ガスに対してレーザ光を照射することにより、キャビティリングダウン吸収分光法による測定を実施する第２測定ステップと、
を実施するように、前記圧力調整部のほか、前記レーザ光照射部及び前記光検出部を制御する構成とすることができる。【００２５】
上述したように、ＣＲＤＳは被測定ガス中の低濃度の成分を高い感度で検出するのに優れた手法である。したがって、本発明において目的成分とは通常、被測定ガス中に比較的低い濃度で含まれる成分であり、典型的には、同じ化学式である同位体中の含有比率の低い同位体、例えば14ＣＯ2、ＤＨＯ（重水）、15ＮＨ3などである。【００２６】
例えばＣＯ2の炭素同位体の一つである14ＣＯ2をＣＲＤＳにより測定する場合、つまり14ＣＯ2が目的成分である場合、被測定ガスの圧力は14ＣＯ2による吸収が他の同位体である12ＣＯ2や13ＣＯ2による吸収に対してできるだけ大きくなるような条件に定められるのが一般的である。この場合、本発明における第１圧力とは、このように14ＣＯ2を測定するのに最適な（又は最適に近い）条件の圧力であり、その圧力の下で目的成分の吸収ピークの波長についてのＣＲＤＳによる測定が実施される。但し、13ＣＯ2による吸収ピークの波長と14ＣＯ2による吸収ピークの波長とはかなり近接しており、上述したように、吸収スペクトルにおいて14ＣＯ2による吸収ピークには13ＣＯ2による吸収ピークが重なる可能性がある。【００２７】
被測定ガスの圧力を14ＣＯ2の測定に最適な圧力条件から変化させると、12ＣＯ2、13ＣＯ2、及び14ＣＯ2のそれぞれの吸収の程度は変化する。圧力や観測する波長によっては、14ＣＯ2による吸収が殆ど無視できる程度になる。また、14ＣＯ2による吸収が無視できない程度に存在する場合でも、14ＣＯ2による吸収の割合は最適な圧力条件のときに比べて大幅に低下する。本発明における第２圧力とは例えば、目的成分である14ＣＯ2による吸収が被測定ガス中に存在する目的成分以外の12ＣＯ2、13ＣＯ2等による吸収に比べて無視できるとき、又は十分に小さいときの圧力である。【００２８】
例えば第２圧力において、目的成分である14ＣＯ2による吸収が目的成分以外の12ＣＯ2、13ＣＯ2等による吸収に比べて無視できる程度に小さい場合、第２測定ステップでの測定では目的成分以外の成分のみによる吸収が反映された結果（リングダウンレート又はリングダウン時間）が得られる。即ち、この測定結果には目的成分による吸収の影響は実質的にないため、そのときの吸収はバックグラウンドであるとみなすことができる。そこで演算ステップでは、被測定ガスの圧力が異なる条件の下で実施された２回のＣＲＤＳによる測定の結果を利用して、バックグラウンドの影響を除去又は軽減するような演算処理を行い、目的成分の濃度を算出する。【００２９】
本発明に係るガス測定方法及びガス測定装置では、第２測定ステップにおける測定で使用されるレーザ光の波長は第１測定ステップにおける測定で使用されるレーザ光の波長と同じでよい。即ち、第１測定ステップと第２測定ステップとで同じ波長のレーザ光を用いてＣＤＲＳ法による測定を実施しても、14ＣＯ2等の低濃度の成分による吸収に対応する信号とこれに重なっている12ＣＯ2、13ＣＯ2等のより高い濃度の成分による吸収に対応する信号とを分離し、低濃度の成分の吸収係数や濃度を算出することができる。【００３０】
ＣＲＤＳ法を利用したガス測定装置においてスペクトルフィッティングを正確に行うために、広範囲に亘ってレーザ光の波長を掃引する場合、広範囲に波長掃引できるレーザ光源が必要になるのは当然であるが、その他に、その波長掃引範囲に対応する高反射ミラーを備えた光共振器が必要である。これに対し本発明では、レーザ光の波長を変えることなく２回の測定を行えばよく、通常、スペクトルフィッティングの正確性のために必要とされる広範囲に波長掃引可能な又は波長切替え可能なレーザ光源やその波長範囲に対応する高反射ミラーを備えた光共振器が不要である。【００３１】
また、被測定ガスの圧力を変化させてバックグラウンドの測定を実施することで、次のような利点もある。【００３２】
上述したようにＣＲＤＳの最大の特徴の一つは検出感度が高いことであるため、目的成分の濃度はかなり低い場合が多く、検出下限に近い状態である場合もしばしばある。そうした場合、目的成分の吸収ピーク位置付近における他の成分による吸収も比較的小さく、フィッティング処理によりバックグラウンドのスペクトル波形を推定する際に用いられる、測定による吸収係数自体があまり正確に得られない可能性もある。そうなると、仮に図１１（ａ）に示したような目的成分以外の成分による吸収ピークが存在していないとしても、バックグラウンドのスペクトル波形が不正確になり、これを用いたバックグラウンド除去の精度が低下して目的成分の濃度の精度も低下する。【００３３】
これに対し、例えば目的成分が14ＣＯ2である場合、被測定ガスの圧力を14ＣＯ2の測定に最適な圧力条件よりも高くすると、14ＣＯ2による吸収に比べて12ＣＯ2や13ＣＯ2による吸収が大きく増加する。そのため、バックグラウンドのレベルが全体的に上がる。その結果、バックグラウンドの測定結果の精度が向上し、バックグラウンド除去を精度良く行うことで目的成分の濃度をより正確に求めることができる。【００３４】
なお、第１圧力及び第２圧力はそれぞれ、目的成分に応じて予め計算により又は実験により決めておけばよい。【００３５】
また、上述したように、第２圧力において、目的成分による吸収が被測定ガス中の目的成分以外の成分による吸収に比べて無視できない場合であっても、圧力によって複数成分による吸収の割合が相違することを利用して、目的成分の吸収係数や濃度を求めることが可能である。【００３６】
即ち、本発明に係るガス測定方法の一つの態様として、
前記第２測定ステップでは、前記目的成分の吸収ピークの波長についてキャビティリングダウン吸収分光法による測定を実施し、
前記演算ステップでは、前記第１測定ステップの結果と前記第２測定ステップの結果とに基づく連立方程式を作成し、該連立方程式を解くことで、前記被測定ガス中の前記目的成分の濃度以外の成分による吸収の影響を除去又は軽減した前記目的成分の濃度を算出するようにしてもよい。【００３７】
また本発明に係るガス測定装置の一つの態様として、
前記制御部は、前記第２圧力の下での測定の際に、前記目的成分の吸収ピークの波長についてキャビティリングダウン吸収分光法による測定を実施し、
前記演算処理部は、前記第１圧力の下での結果と前記第２圧力の下での結果とに基づく連立方程式を作成し、該連立方程式を解くことで、前記被測定ガス中の前記目的成分の濃度以外の成分による吸収の影響を除去又は軽減した前記目的成分の濃度を算出する構成としてもよい。【００３８】
これらの態様では、バックグラウンドに目的成分による吸収の影響も含まれるため、目的成分の濃度や吸収係数、目的成分以外の成分の濃度や吸収係数を未知の値とする連立方程式を解くことで、目的成分のみの濃度を算出する。これにより、目的成分による吸収の影響が完全になくなるような圧力まで被測定ガスの圧力を変更することができない場合であっても、バックグラウンドの影響を適切に除去して目的成分の吸収係数や濃度を得ることができる。【００３９】
また、第２測定ステップでの測定で使用するレーザ光の波長は、第１測定ステップでの測定で使用するレーザ光の波長と必ずしも同じである必要はない。
即ち、本発明に係るガス測定方法の別の態様として、前記第２測定ステップでは、前記目的成分の吸収ピークの波長とは異なる、該目的成分による吸収の影響を無視できる波長についてキャビティリングダウン吸収分光法による測定を実施し、
前記演算ステップでは、前記第２測定ステップの結果に基づいて前記第２圧力の下での前記被測定ガス中の前記目的成分以外の成分の濃度を推定し、該濃度から、前記第１測定ステップの結果から求まる吸収係数における前記目的成分以外の成分の吸収の寄与を推定し、その影響を除去する演算を行うようにするとよい。【００４０】
また本発明に係るガス測定装置の別の態様として、
前記制御部は、前記第２圧力の下での測定の際に、前記目的成分の吸収ピークの波長とは異なる、該目的成分による吸収の影響を無視できる波長についてキャビティリングダウン吸収分光法による測定を実施し、
前記演算処理部は、前記第２圧力の下での結果に基づいて該第２圧力の下での前記被測定ガス中の前記目的成分以外の成分の濃度を推定し、該濃度から、前記第１圧力の下での結果から求まる吸収係数における前記目的成分以外の成分の吸収の寄与を推定し、その影響を除去する演算を行う構成としてもよい。【００４１】
この場合、第２測定ステップでの測定に用いるレーザ光の波長は、目的成分による吸収の影響を無視できることが事前に明らかである適宜の波長に定めればよい。これらの態様では、第１測定ステップと第２測定ステップとで測定に用いるレーザ光の波長を切り替える必要があるものの、第２測定ステップでの測定による測定結果には目的成分による吸収の影響は実質的にない。そのため、目的成分以外の成分のみによる吸収、つまりはバックグラウンドであるとみなすことができるので、上述したような連立方程式を解く場合に比べて、バックグラウンドの除去は容易である。また、目的成分による吸収の影響が完全になくなるような圧力まで被測定ガスの圧力を変更することができない場合であっても、比較的容易にバックグラウンド除去処理を行うことができる。【００４２】
また、上記態様ではレーザ光の波長を切り替える必要があるものの、一般には、目的成分の別の吸収ピークを用いる場合のように遠く離れた波長における吸収係数を測定する必要はなく、狭い波長範囲内で波長を切り替えれば十分である。そのため、ＣＲＤＳ装置における一般的な光源やミラーで問題なく実現することができる。【００４３】
なお、本発明に係るガス測定装置において、前記圧力調整部は、前記測定セル中に第２圧力で被測定ガスを封入した状態から、該測定セル中から一部の被測定ガスを外部に強制的に排出することで該測定セル中の被測定ガスの圧力を第１圧力に調整するものとすることができる。【００４４】
具体的には、測定セルに接続されたガス導入管とガス排出管にそれぞれ設けられた開閉バルブと、ガス排出管を通して測定セル中の被測定ガスを外部へと排出する真空ポンプと、測定セル中のガスの圧力を検出する圧力検出部と、該圧力検出部により圧力をモニタしつつ前記開閉バルブの開閉動作と前記真空ポンプの動作とを制御する圧力制御部と、を含む構成とすることができる。【００４５】
或いは本発明に係るガス測定装置において、前記圧力調整部は、前記測定セル中に被測定ガスを供給し第１の所定の圧力で被測定ガスを封入した状態から、先に供給されずに残存されていた被測定ガスを該測定セル中に追加供給して封入することで該測定セル中の被測定ガスの圧力を第２の所定の圧力に調整するものである構成としてもよい。【００４６】
具体的には、測定セルに接続されたガス導入管とガス排出管にそれぞれ設けられた開閉バルブと、ガス導入管を通して測定セル中に被測定ガスを供給する送給ポンプと、測定セル中のガスの圧力を検出する圧力検出部と、該圧力検出部により圧力をモニタしつつ前記開閉バルブの開閉動作と前記送給ポンプの動作とを制御する圧力制御部と、を含む構成とすることができる。
これら構成によれば、測定セル中の被測定ガスの圧力を容易に、目的とする値に調整することができる。
【発明の効果】
【００４７】
本発明によれば、被測定ガスに対する２回の測定によって、その被測定ガス中の目的成分以外の成分による吸収に起因するバックグラウンドを高い精度で以て除去し、目的成分の正確な濃度を取得することができる。これにより、バックグラウンドのスペクトルを推定するのに必要な多数回の繰り返し測定が不要になるので、測定時間を短縮し測定のスループットを向上させることができる。また、測定時間が短いために、例えば半減期が短い放射性同位体のように目的成分が比較的不安定なものであっても、正確な濃度測定を行うことが可能となる。【００４８】
さらにまた、本発明によれば、測定時間が短いため、ミラーでの測定ガスの一部の吸着による実効的なミラー反射率の低下や、ミラーの微小な変動や入射光位置の微小な変動による実効的な反射率や共振器長の変動、温度の微小な変動によって生じる光共振器の熱膨張による共振器長等の変動といった測定中に生じる測定状態の変化に対する影響を、最小限又はそれに近い状態に抑えることができる。【００４９】
また本発明によれば、目的成分による吸収ピークに他の成分による吸収ピークが重なっているような場合であっても、バックグラウンドを正確に推定してバックグラウンドを除去することで、目的成分の吸収係数や濃度を精度良く求めることができる。
【図面の簡単な説明】
【００５０】
【図１】本発明の一実施例であるＣＲＤＳ装置の要部の構成図。
【図２】圧力が１０１３．２５Paである場合におけるＣＯ2同位体ガスによる吸収特性の計算結果を示す図。
【図３】圧力が１０１３２．５Paである場合におけるＣＯ2同位体ガスによる吸収特性の計算結果を示す図。
【図４】本実施例のＣＲＤＳ装置において目的成分の濃度を求める際の測定及び処理の手順の一例を示すフローチャート。
【図５】本実施例のＣＲＤＳ装置において目的成分の濃度を求める際の測定及び処理の手順の他の例を示すフローチャート。
【図６】圧力が１０１３．２５Paである場合におけるＨ2Ｏ同位体による吸収特性の計算結果を示す図。
【図７】圧力が１０１３２５Paである場合におけるＨ2Ｏ同位体による吸収特性の計算結果を示す図。
【図８】一般的なＣＲＤＳ装置の概略構成図。
【図９】光共振器でのモード周波数とレーザ光の発振周波数との関係を示す概略図。
【図１０】14ＣＯ2吸収線のピーク波長付近における被測定ガスの吸収スペクトルの概略図。
【図１１】ＣＯ2同位体ガスに対する吸収スペクトルの計算結果を示す図（ａ）及び14ＣＯ2の吸収ピーク位置での12ＣＯ2、13ＣＯ2、及び14ＣＯ2による吸収の寄与度合いの計算結果を示す図（ｂ）。
【図１２】被測定ガスの圧力を高くした状態でのＣＯ2同位体ガスに対する吸収スペクトルの計算結果を示す図（ａ）及び14ＣＯ2の吸収ピーク位置での12ＣＯ2、13ＣＯ2、及び14ＣＯ2による吸収の寄与度合いの計算結果を示す図（ｂ）。
【発明を実施するための形態】
【００５１】
以下、本発明に係るガス測定装置の一実施例であるＣＲＤＳ装置及び該装置を用いたガス測定方法について、添付図面を参照して説明する。
まず、すでに図１１及び図１２を用いて、上述した本発明の課題を整理して説明するとともに、本発明におけるバックグラウンド除去の原理について説明する。
【００５２】
図１１（ａ）と図１２（ａ）はいずれも、14ＣＯ2の吸収ピーク位置付近での12ＣＯ2、13ＣＯ2、及び14ＣＯ2についての吸収スペクトルの計算結果であり、図１１（ａ）と図１２（ａ）との相違は想定した被測定ガスの圧力のみである。即ち、圧力が相対的に低い（ここでは０．０３atm）ときには、14ＣＯ2の吸収ピークが明確に観測されるものの、この14ＣＯ2の吸収ピークに裾部が重なっている13ＣＯ2の吸収ピークは単独で確認することができない。そのため、13ＣＯ2の吸収に対して正しくスペクトルフィッティングを行うことはできず、バックグラウンド（ピークのベースライン）を正確に推定することは困難である。
【００５３】
また、一般にガス分子は、回転、並進、及び振動に応じた複数の吸収ピークを有しているため、上述したように或る波長における吸収ピークで信号同士が重なり合っていたとしても、別の波長における別の吸収ピークを用いることで吸収係数を測定できる場合もある。しかしながら、ＣＲＤＳ装置で用いる光源とミラーはそれぞれ限られた波長範囲にのみ有効であり、被測定ガス中の目的成分の複数の吸収ピークの波長には対応できないことがしばしばある。そうした場合、吸収係数を測定する有効な方法が存在しないのが実状である。
【００５４】
従来のＣＲＤＳによる測定では、被測定ガスの圧力や温度を一定とすることが前提であり、14ＣＯ2の濃度を求めたい場合には、図１１（ａ）に示したように14ＣＯ2の吸収ピークができるだけ高くなるような圧力に設定されるのが一般的である。一方、図１２（ａ）に示すように、被測定ガスの圧力を適当な圧力まで高めると、14ＣＯ2の吸収ピークとともに13ＣＯ2の吸収ピークも殆ど消滅する。但し、図１２（ａ）、（ｂ）に示されているように、13ＣＯ2による吸収は無くなったわけではなく、あくまでも吸収ピークが消滅しただけであり、吸収自体は存在している。また、より濃度が高い12ＣＯ2による吸収も存在する。即ち、図１２（ａ）に示すように被測定ガスの圧力を高くした状態では、14ＣＯ2の吸収ピーク位置における吸収は目的成分（14ＣＯ2）以外の成分によるもの、つまりはバックグランドのみであるとみなせる。なお、圧力を高めることで吸収量は増加するため、吸収係数自体が大きく増加していることに注意すべきである。
【００５５】
上述したように被測定ガスの圧力を高めた状態でＣＲＤＳによる測定を実施すると、目的成分以外の成分による吸収を反映した結果（リングダウンレート又はリングダウン時間）が求まるから、その結果から計算される吸収係数に基づいて目的成分以外の成分の絶対濃度を算出することができる。この絶対濃度から、図１１（ａ）に示すような被測定ガスの圧力が相対的に低い状態での14ＣＯ2の吸収ピーク位置におけるベースラインのスペクトルを推算することができる。そして、被測定ガスの圧力が相対的に低い状態でＣＲＤＳによる測定を実施することで得られた測定結果から求まる吸収係数からベースラインを差し引くことで純粋な14ＣＯ2の吸収係数を求め、この吸収係数から目的成分である14ＣＯ2のみの絶対濃度を算出することができる。
【００５６】
なお、被測定ガスの圧力によっては、14ＣＯ2の吸収の度合いを無視できない場合がある。その場合でも後述するように連立方程式を解くことにより、目的成分である14ＣＯ2のみの絶対濃度を算出することができる。また、被測定ガスの圧力を変えるとともに測定に使用するレーザ光の波長も変えることで、目的成分による吸収の影響がないベースラインのスペクトルを反映した測定結果を取得し、これを用いて純粋な14ＣＯ2の吸収係数を求めることも考えられる。これについても後述する。
【００５７】
次に上述した測定原理を用いたＣＲＤＳ装置の一実施例について説明する。図１は本実施例のＣＲＤＳ装置の要部の構成図である。
【００５８】
本実施例のＣＲＤＳ装置は、測定系として、レーザ光源部１、レーザ駆動部２、光スイッチ３、光共振器４、及び光検出器５を備える。光共振器４は、被測定ガスであるサンプルガスが収容される略円筒状の測定セル４０と、該測定セル４０の両端に対向して配置された一対の高反射率のミラー４７、４８と、を含む。測定セル４０にはガス導入管４１とガス排出管４３とが接続され、ガス導入管４１には導入バルブ４２が設けられ、ガス排出管４３には排出バルブ４４と真空ポンプ４５とが設けられている。また、測定セル４０には該セル４０内に収容されているガスの圧力を検出するための圧力センサ４６が付設されている。
【００５９】
制御部６は後述する測定やデータ処理を実行するためにレーザ駆動部２等の各部を制御するものであり、機能的なブロックとして、測定制御部６１、レーザ制御部６２、圧力制御部６３、及び測定パラメータ記憶部６４などを含む。測定パラメータ記憶部６４には、レーザ光の波数（又は波長）や圧力などの測定パラメータが測定対象である成分の種類などに対応して予め格納される。光検出器５による検出信号が入力されるデータ処理部７は、機能的なブロックとして、測定データ格納部７１、リングダウン時間算出部７２、濃度演算部７３、及び演算用既知情報保存部７４などを備える。測定データ格納部７１はアナログ検出信号をデジタル化するアナログデジタル変換器を含む。また、データ処理部７に接続された出力部８は例えば表示モニタなどである。
【００６０】
本実施例のＣＲＤＳ装置において、被測定ガスがＣＯ2であり、目的成分がＣＯ2の同位体の一つである14ＣＯ2である場合を例として具体的な動作を説明する。なお、放射性同位体である炭素14Ｃを含む14ＣＯ2の濃度測定は様々な分野で広く利用されている。
【００６１】
図２（ａ）は、圧力が１０１３．２５Pa（＝０．０１atm）であるときのＣＯ2同位体ガスについてＣＲＤＳによる吸収スペクトルを計算した結果を示す図であり、横軸は光の波数、縦軸は吸収係数である。図２（ｂ）は図２（ａ）中に条件１で示す光の波数（14ＣＯ2の吸収ピークの波数）における吸収の要因の内訳を示す円グラフである。これら計算に際して、温度は２００Ｋ、14ＣＯ2濃度は２×１０-12で14ＣＯ2以外のＣＯ2濃度は０．２であると想定した。また、被測定ガスに含まれる14ＣＯ2以外のＣＯ2同位体は天然同位体存在比で以て測定セル４０に導入されると想定した。
【００６２】
図２（ａ）中に記載されているように、上記圧力の下で、14ＣＯ2の吸収ピークの波数における吸収係数は３．４９×１０-10であると計算される。図２（ｂ）で分かるように、このときの光の吸収は約８２％が14ＣＯ2によるものであるが、残りの約１８％は14ＣＯ2以外のＣＯ2同位体（12ＣＯ2、13ＣＯ2）によるものである。したがって、14ＣＯ2の正確な濃度を算出するには、14ＣＯ2以外のＣＯ2の吸収によるベースラインを差し引くことが必要である。このベースラインのスペクトルを求めるには、14ＣＯ2以外のＣＯ2の濃度が分かればよい。そのためには、14ＣＯ2以外のＣＯ2による吸収係数が測定データから求まればよいが、ここで問題となるのは、14ＣＯ2による吸収ピークが図２（ａ）に示すように急峻に観測されるような、14ＣＯ2に適した測定条件の下では、14ＣＯ2以外のＣＯ2による光の吸収量が小さいためにその吸収量を反映したデータが正確に得られない、という点である。この問題を克服するために本発明では、被測定ガスの圧力を変化させて測定を実施し、その測定結果を利用して14ＣＯ2以外のＣＯ2の吸収によるベースラインのスペクトルを求めている。
【００６３】
一般的なＣＲＤＳでは、圧力を常に一定に維持した状態で被測定ガスに対し測定を実行する。よく知られているように、ガス中の成分の吸収係数は、温度、圧力、光の波長などに依存する。そこで一般に、14ＣＯ2による吸収を測定する際の圧力は、目的成分である14ＣＯ2による吸収ピークの波数において14ＣＯ2の吸収係数と14ＣＯ2以外のＣＯ2の吸収係数との差ができるだけ大きくなる等の条件を満たすような圧力に設定される。何故なら、こうした圧力が14ＣＯ2による吸収ピークを観測するうえで最もＳＮ比が良好な条件と考えられるからである。実際に、この圧力よりも被測定ガスの圧力を高くしていくと、12ＣＯ2や13ＣＯ2による吸収ピークの高さがそれぞれの吸収ピークの波数において大幅に上昇するとともにそのピーク幅も広がる。その結果、バックグラウンドのレベルがかなり高くなり、14ＣＯ2による吸収ピークのＳＮ比は低くなる。
【００６４】
図３（ａ）は、圧力が図２（ａ）の場合の１０倍である１０１３２．５Pa（０．１atm）であるときのＣＯ2同位体ガスについて、ＣＲＤＳによる吸収スペクトルを計算した結果を示す図である。図３（ｂ）は図３（ａ）中に条件２で示す波数（上述の条件１における波数と同一）における吸収の要因の内訳を示す円グラフ、図３（ｃ）は図３（ａ）中に条件３で示す波数（条件１よりも小さい適宜の波数）における吸収の要因の内訳を示す円グラフである。なお、圧力以外の計算条件は図２の場合と同じである。
【００６５】
12ＣＯ2による吸収ピークの波数は図３（ａ）に示したグラフから左に外れる範囲に存在するが、その吸収ピークの高さとピーク幅とは圧力上昇に伴って急激に大きくなる。このときの14ＣＯ2による吸収ピークの高さは図２（ａ）に示すグラフにおける吸収ピークの高さの１０倍以上であるものの、その吸収ピークは12ＣＯ2による吸収ピークのテーリングに殆ど埋もれてしまっている。即ち、このように被測定ガスの圧力を高くするとバックグラウンド全体のレベルがかなり高くなり、14ＣＯ2以外のＣＯ2による吸収を検出し易くなる、又はその検出の精度が向上することが分かる。そこで本発明では、目的成分（14ＣＯ2）の吸収の割合が相対的に大きい圧力条件の下でＣＲＤＳの測定を行うほかに、このようにバックグラウンドが高くなる圧力条件の下で同じ被測定ガスに対するＣＲＤＳの測定を実行する。
【００６６】
図３（ｂ）、（ｃ）に示すように、条件２の波数においては14ＣＯ2による吸収の割合は７％程度残るが、条件３の波数においては14ＣＯ2による吸収の割合は０％である。このバックグラウンドを求めるための相対的に高い圧力の下での測定は条件２又は条件３のいずれか一つで行えばよいが、いずれの条件で測定を実施するのかによって測定結果の処理方法が異なる。なお、条件３の波数は図３（ａ）に示した位置に限るものではなく、14ＣＯ2の吸収ピークがバックグラウンドに埋もれ且つバックグラウンドのレベルが高く（つまりは図３（ａ）で14ＣＯ2の吸収ピークの位置よりも左側であり）、さらにレーザ光の波数の調整可能範囲な範囲でありさえすれば適宜に定めることができる。
【００６７】
［条件３：14ＣＯ2による吸収を無視できる場合］
条件３の波数では、14ＣＯ2による吸収が無視できる程度に12ＣＯ2や13ＣＯ2の吸収に起因するバックグラウンドが大きい。そのため、ＣＲＤＳによる測定結果には14ＣＯ2による吸収の影響が実質的に現れない。そこで、このときに得られる測定データから求まるリングダウン時間に基づいて14ＣＯ2以外のＣＯ2同位体の吸収係数を算出し、その吸収係数から相対的に高い圧力の下での4ＣＯ2以外のＣＯ2同位体の濃度を算出することができる。そして、算出された14ＣＯ2以外のＣＯ2同位体の濃度を用いることで、14ＣＯ2の吸収が大きい、相対的に低い圧力の下での条件１の波数におけるバックグラウンドに相当する吸収係数を算出することができる。そこで、条件１において得られた測定結果から求まる吸収係数からバックグラウンドに相当する吸収係数を差し引くことで14ＣＯ2のみの吸収係数を求め、この吸収係数から目的成分である14ＣＯ2のみの濃度を算出することができる。
【００６８】
［条件２：14ＣＯ2による吸収を無視できない場合］
条件２の波数では14ＣＯ2による吸収が７％程度の割合で存在するため、14ＣＯ2による吸収を無視することができない。この場合には、相対的に低い圧力の下で得られた14ＣＯ2の吸収ピークの波数における測定結果と、相対的に高い圧力の下で得られた条件２の波数における測定結果とに基づいて、14ＣＯ2の濃度とそれ以外のＣＯ2同位体の濃度とをそれぞれ未知の値とする連立方程式を作成してそれを解く。それにより、14ＣＯ2の濃度とそれ以外のＣＯ2同位体の濃度とを算出することができる。或いは、14ＣＯ2の吸収係数とそれ以外のＣＯ2同位体の吸収係数とをそれぞれ未知の値とする連立方程式を作成してそれを解いてもよい。
【００６９】
いずれにしても、上記のようにして算出された14ＣＯ2の濃度には14ＣＯ2以外のＣＯ2同位体による吸収の影響が含まれない又はその影響が無視できる程度であるので、高い精度で以て被測定ガス中の14ＣＯ2の濃度を求めることができる。なお、それぞれの測定を行う際の圧力の条件や使用するレーザ光の波数は測定対象の成分の種類等に応じて予め適切に決めておけばよい。
【００７０】
また、上記説明では12ＣＯ2と13ＣＯ2とを区別せずに14ＣＯ2以外のＣＯ2同位体の濃度を求めていたが、12ＣＯ2、13ＣＯ2、14ＣＯ2の濃度の全てを求めたい場合には、各同位体ガスの吸収の割合が大きくなるそれぞれ異なる圧力条件で以て測定を行い、その三つの測定結果に基づいて各同位体ガスの濃度を算出すればよい。
【００７１】
本実施例のＣＲＤＳ装置における被測定ガス中の目的成分（14ＣＯ2）濃度の測定動作のフローチャートを図４及び図５に示す。
【００７２】
図４は、条件３の下での測定結果を利用してバックグラウンド除去を行う場合の、測定及び処理の手順の一例を示すフローチャートである。なお、ＣＯ2を含む被測定ガスが測定セル４０に導入されていない状態での各圧力の下でのリングダウン時間は予め測定され、演算用既知情報保存部７４に格納されているものとする。また、濃度演算の際に用いられる目的成分の吸収断面積などの先験情報も演算用既知情報保存部７４に格納されているものとする。
【００７３】
まず、制御部６において圧力制御部６３は排出バルブ４４を閉鎖した状態で導入バルブ４２を開放し、測定セル４０内に被測定ガスを導入する。そして、圧力センサ４６により検出される圧力が所定値になったならば導入バルブ４２を閉鎖して測定セル４０内に被測定ガスを充満させる（ステップＳ１）。次いで圧力制御部６３は、排出バルブ４４を開放するとともに真空ポンプ４５を動作させ、測定セル４０内の被測定ガスをガス排出管４３を通して排出し始める。そして、圧力センサ４６により検出される圧力が測定パラメータ記憶部６４に保存されている所定のバックグラウンド（ＢＧ）測定圧Ｐ3にまで下がったならば排出バルブ４４を閉鎖する（ステップＳ２）。これにより、測定セル４０内には圧力Ｐ3の被測定ガスが封入された状態となる。
【００７４】
レーザ制御部６２はレーザ光の波数が予め定められているバックグラウンド（ＢＧ）測定用値ν3になるようにレーザ駆動部２を通してレーザ光源部１を動作させる（ステップＳ３）。そして、測定制御部６１は、レーザ光波数ν3、圧力Ｐ3の下での測定を実行する。即ち、測定セル４０中の被測定ガスにレーザ光を照射し、所定のタイミングで光スイッチ３によりレーザ光を遮断する。そして、レーザ光を遮断する直前から所定の時間が経過するまでの間、光検出器５により検出されたデータを収集する（ステップＳ４）。このときに光検出器５により得られる高圧力の下での測定データは、測定データ格納部７１に一旦保存される。このときの測定データは条件３でのリングダウン時間ｔ3を情報として含むデータである。
【００７５】
そのあと圧力制御部６３は、排出バルブ４４を再び開放するとともに真空ポンプ４５を動作させ、ガス排出管４３を通して測定セル４０内の被測定ガスを外部へ排出し始める。そして、圧力センサ４６により検出される圧力が演算用既知情報保存部７４に保存されている目的成分測定圧Ｐ1にまで下がったならば、排出バルブ４４を閉鎖する（ステップＳ５）。これにより、測定セル４０内には圧力Ｐ3よりも低い圧力Ｐ1の被測定ガスが封入された状態となる。
【００７６】
一方、レーザ制御部６２はレーザ光の波数が目的成分測定用値ν1になるようにレーザ駆動部２を通してレーザ光源部１を動作させる（ステップＳ６）。そして、測定制御部６１はレーザ光波数ν1、圧力Ｐ1の下での測定を実行し、ステップＳ４と同様に、所定期間の測定データを取得する（ステップＳ７）。このときに光検出器５により得られる一連の測定データも測定データ格納部７１に一旦保存される。この低圧力の下での測定データは、図２（ａ）に示した条件１でのリングダウン時間ｔ1を情報として含むデータである。
【００７７】
なお、厳密に言えば、測定セル４０内の圧力を下げるためにステップＳ５では測定セル４０内の被測定ガスの一部を外部へ排出するため、ステップＳ４とＳ７とで測定される被測定ガスは全く同一ではない。しかしながら、測定セル４０内における被測定ガス中の成分の分布は均一であるとみなせるので、ステップＳ４とＳ７とで測定されるガスは同じ被測定ガスであり圧力のみが相違するものとみなすことができる。
【００７８】
データ処理部７においてリングダウン時間算出部７２は測定データ格納部７１に保存した高圧力の下での測定データに基づいてリングダウン時間ｔ3を算出する（ステップＳ８）。そして、濃度演算部７３は、その算出結果と演算用既知情報保存部７４に保存されている条件３の下での被測定ガス無し時のリングダウン時間とに基づいて吸収係数を計算し、さらにその吸収係数から濃度を求める（ステップＳ９）。これらの計算方法は従来と同じであり、例えば上記(1)、(2)式を用いればよい。このときには14ＣＯ2による吸収の影響は無視できるので、ステップＳ８において求まるのは被測定ガス中の14ＣＯ2以外のＣＯ2同位体ガスの濃度である。
【００７９】
引き続いてリングダウン時間算出部７２は、測定データ格納部７１に保存した低圧力の下での測定データに基づいてリングダウン時間ｔ1を算出する（ステップＳ１０）。そして、濃度演算部７３は、その算出結果と演算用既知情報保存部７４に保存されている条件１の下での被測定ガス無し時のリングダウン時間とから吸収係数を計算する（ステップＳ１１）。このときに求まるのは被測定ガス中の14ＣＯ2を含むＣＯ2同位体ガスの吸収係数である。
【００８０】
ステップＳ８では14ＣＯ2を除くＣＯ2同位体ガスの濃度が得られている。そこで、濃度演算部７３はこの濃度から条件１の圧力及びレーザ波数の下での14ＣＯ2を除くＣＯ2同位体ガスによる吸収係数を計算する。この吸収係数がバックグラウンドであるから、14ＣＯ2を含むＣＯ2同位体ガスによる吸収係数から14ＣＯ2を除くＣＯ2同位体による吸収係数を差し引き、条件１における14ＣＯ2による吸収係数を算出し、この吸収係数から14ＣＯ2のみの濃度を算出する（ステップＳ１２）。そして、その結果を出力部８から出力する。
以上のようにして、本実施例のＣＲＤＳ装置では、バックグラウンドが除去された14ＣＯ2のみの正確な濃度の情報をユーザに提供することができる。
【００８１】
次に、条件２の下での測定結果を利用してバックグラウンド除去を行う場合の、測定及び処理の手順の一例を、図５に示すフローチャートに従って説明する。この場合にも、ＣＯ2を含む被測定ガスが測定セル４０に導入されていない状態での各圧力の下でのリングダウン時間や目的成分の吸収断面積などの情報は、演算用既知情報保存部７４に格納されているものとする。
【００８２】
ステップＳ２１、Ｓ２２は上記ステップＳ１、Ｓ２と全く同じ処理であり、これらステップによって、測定セル４０内に圧力Ｐ3の被測定ガスを封入する。レーザ制御部６２はレーザ光の波数が目的成分測定用値ν1になるようにレーザ駆動部２を通してレーザ光源部１を動作させる（ステップＳ２３）。そして、測定制御部６１はレーザ光波数ν1、圧力Ｐ3の下で測定を実行し、測定データを取得する（ステップＳ２４）。この高圧力の下での測定データは、条件２の下でのリングダウン時間ｔ2を情報として含むデータである。
【００８３】
次にステップＳ５と同じステップＳ２５の処理を行うことで、測定セル４０内に封入している被測定ガスの圧力を目的成分測定圧Ｐ1にまで下げる。そしてレーザ波数を目的成分測定用値ν1に維持したまま、測定制御部６１はレーザ光波数ν1、圧力Ｐ1の下での測定を実行し、測定データを取得する（ステップＳ２６）。この低圧力の下での測定データは、図４の例と同様の、条件１の下でのリングダウン時間ｔ1を情報として含むデータである。
【００８４】
データ処理部７においてリングダウン時間算出部７２は、測定データ格納部７１に保存した高圧力の下での測定データに基づいてリングダウン時間ｔ2を算出する（ステップＳ２７）。濃度演算部７３はその算出結果と条件２の下での被測定ガス無し時のリングダウン時間とから吸収係数α2を計算する（ステップＳ２８）。このときには14ＣＯ2による吸収の影響を無視できないので、ここで求まるのは、条件２における被測定ガス中の14ＣＯ2による吸収係数と14ＣＯ2以外のＣＯ2同位体ガスによる吸収係数とを加えた値である。
【００８５】
続いてリングダウン時間算出部７２は、測定データ格納部７１に保存した低圧力の下での測定データに基づいてリングダウン時間ｔ1を算出する（ステップＳ２９）。これはステップＳ１０と同じである。そして、濃度演算部７３は、その算出結果と条件１の下での被測定ガス無し時のリングダウン時間とから吸収係数α1を計算する（ステップＳ３０）。このときに求まるのは、条件１における被測定ガス中の14ＣＯ2による吸収係数と14ＣＯ2以外のＣＯ2同位体による吸収係数とを加えた値である。
【００８６】
ここでは、14ＣＯ2の濃度ｘと14ＣＯ2以外のＣＯ2同位体の濃度ｙとがそれぞれ未知の値である。そこで、条件１の下での測定で得られた吸収係数α1と濃度ｘ、ｙの関係を示す式と、条件２の下での測定で得られた吸収係数α2と濃度ｘ、ｙの関係を示す式とを連立方程式として作成する。そして、濃度演算部７３は、この連立方程式を解くことで14ＣＯ2のみの濃度を算出する（ステップＳ３１）。
以上のようにして、バックグラウンドが除去された14ＣＯ2のみの正確な濃度を算出し、ユーザに提供することができる。
【００８７】
なお、上記の例では、リングダウン時間から吸収係数を求めたあと、所定の計算式に従って14ＣＯ2の濃度を算出していたが、計算式を用いる代わりにデータベースを参照して14ＣＯ2の濃度の導出できるようにしてもよい。即ち、それぞれの測定条件（つまりは条件１、２、３）において、14ＣＯ2とそれ以外のＣＯ2同位体の様々な濃度の組合せに対して観測される吸収係数の値を予め求めてデータベース化しておく。そして、測定データに基づいて吸収係数が求まったならば、その吸収係数を入力としてデータベース検索を行うことで、対応する14ＣＯ2の濃度と14ＣＯ2以外のＣＯ2同位体の濃度とを導出するようにしてもよい。
【００８８】
また上記実施例の説明では、ＣＯ2同位体の一つである14ＣＯ2の濃度を求めていたが、本実施例のＣＲＤＳ装置が被測定ガス中の他の成分の分析にも利用できることは当然である。ここでは、別の応用例として、ＣＯ2同位体の測定と並んで利用価値が高い被測定ガス中のＨ2Ｏ同位体の測定について簡単に述べる。
【００８９】
図６（ａ）は、ガス圧力が１０１３．２５Pa（＝０．０１atm）であるときのＨ2Ｏ同位体であるＤＨＯの吸収ピーク測定を行う場合に得られる、ＣＲＤＳによる吸収スペクトルを計算した結果を示す図である。図６（ｂ）は図６（ａ）中に条件１で示す波数における吸収の要因の内訳を示す円グラフである。一方、図７（ａ）はガス圧力が図６の１００倍である１０１３２５Pa（＝１atm）であるときのＨ2Ｏ同位体であるＤＨＯの吸収ピーク測定を行う場合に得られる、ＣＲＤＳによる吸収スペクトルを計算した結果を示す図である。図７（ｂ）は図７（ａ）において条件１と同じ波数における吸収の要因の内訳を示す円グラフである。
これら計算に際し、温度は３５３Ｋ、被測定ガス中のＨ2Ｏ2の濃度は０．０３であると想定した。また、被測定ガスに含まれるＨ2Ｏ同位体は全て天然同位体存在比で以て測定セル４０に導入されると想定した。
【００９０】
図６に示すように、ガス圧力が１０１３．２５Paであるとき、Ｈ2Ｏの同位体の一つであるＤＨＯの吸収ピークの波数ではＤＨＯの吸収の割合は９３％である。一方、図７に示すように、ガス圧力がその１００倍である高い圧力の下では、吸収の度合いが全体に大きく増加するとともに、ＤＨＯの吸収の割合は２％と大きく減り、殆どがそれ以外のＨ2Ｏ同位体による吸収が占める。したがって、この場合にもＤＨＯの吸収を完全に無視することはできないから、上記例における条件２と同じようにしてバックグラウンドを除去することで、目的とするＤＨＯの正確な絶対濃度を求めることができる。
【００９１】
もちろん、それ以外の様々な同位体ガスにおいて特に濃度が低い同位体の濃度を測定したい場合に、本発明に係る装置や方法が有効なことは明らかである。
【００９２】
なお、上記実施例はいずれも本発明の一例であり、本発明の趣旨の範囲で適宜に変形や修正、追加などを行っても、本願特許請求の範囲に包含されることは明らかである。
【符号の説明】
【００９３】
１…レーザ光源部
２…レーザ駆動部
３…光スイッチ
４…光共振器
４０…測定セル
４１…ガス導入管
４２…導入バルブ
４３…ガス排出管
４４…排出バルブ
４５…真空ポンプ
４６…圧力センサ
４７、４８…ミラー
５…光検出器
６…制御部
６１…測定制御部
６２…レーザ制御部
６３…圧力制御部
６４…測定パラメータ記憶部
７…データ処理部
７１…測定データ格納部
７２…リングダウン時間算出部
７３…濃度演算部
７４…演算用既知情報保存部
８…出力部
【図１】
2020105118.tif 000003
【図２】
2020105118.tif 000004
【図３】
2020105118.tif 000005
【図４】
2020105118.tif 000006
【図５】
2020105118.tif 000007
【図６】
2020105118.tif 000008
【図７】
2020105118.tif 000009
【図８】
2020105118.tif 000010
【図９】
2020105118.tif 000011
【図１０】
2020105118.tif 000012
【図１１】
2020105118.tif 000013
【図１２】
2020105118.tif 000014
【手続補正書】
【提出日】令和3年2月8日(2021.2.8)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００５１
【補正方法】変更
【補正の内容】
【００５１】
以下、本発明に係るガス測定装置の一実施例であるＣＲＤＳ装置及び該装置を用いたガス測定方法について、添付図面を参照して説明する。
まず、図１１及び図１２を用いて、上述した本発明の課題を整理して説明するとともに、本発明におけるバックグラウンド除去の原理について説明する。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００８９
【補正方法】変更
【補正の内容】
【００８９】
図６（ａ）は、ガス圧力が１０１３．２５Pa（＝０．０１atm）であるときのＨ2Ｏ同位体であるＤＨＯの吸収ピーク測定を行う場合に得られる、ＣＲＤＳによる吸収スペクトルを計算した結果を示す図である。図６（ｂ）は図６（ａ）中に下向き矢印で示す波数（ＤＨＯ吸収ピークの位置）における吸収の要因の内訳を示す円グラフである。一方、図７（ａ）はガス圧力が図６の１００倍である１０１３２５Pa（＝１atm）であるときのＨ2Ｏ同位体であるＤＨＯの吸収ピーク測定を行う場合に得られる、ＣＲＤＳによる吸収スペクトルを計算した結果を示す図である。図７（ｂ）は図７（ａ）中に下向き矢印で示す波数（図６（ａ）中の下向き矢印と同じ波数）における吸収の要因の内訳を示す円グラフである。
これら計算に際し、温度は３５３Ｋ、被測定ガス中のＨ2Ｏ2の濃度は０．０３であると想定した。また、被測定ガスに含まれるＨ2Ｏ同位体は全て天然同位体存在比で以て測定セル４０に導入されると想定した。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００９０
【補正方法】変更
【補正の内容】
【００９０】
図６（ｂ）に示すように、ガス圧力が１０１３．２５Paであるとき、Ｈ2Ｏの同位体の一つであるＤＨＯの吸収ピークの波数ではＤＨＯの吸収の割合は９３％である。一方、図７（ｂ）に示すように、ガス圧力がその１００倍である高い圧力の下では、吸収の度合いが全体に大きく増加するとともに、ＤＨＯの吸収の割合は２％と大きく減り、殆どがそれ以外のＨ2Ｏ同位体による吸収が占める。したがって、この場合にもＤＨＯの吸収を完全に無視することはできないから、上記例における条件２と同じようにしてバックグラウンドを除去することで、目的とするＤＨＯの正確な絶対濃度を求めることができる。
【手続補正書】
【提出日】令和3年3月15日(2021.3.15)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項３
【補正方法】変更
【補正の内容】
【請求項３】
請求項１に記載のガス測定方法であって、
前記第２測定ステップでは、前記目的成分の吸収ピークの波長についてキャビティリングダウン吸収分光法による測定を実施し、
前記演算ステップでは、前記第１測定ステップの結果と前記第２測定ステップの結果とに基づく連立方程式を作成し、該連立方程式を解くことで、前記被測定ガス中の前記目的成分以外の成分による吸収の影響を除去又は軽減した前記目的成分の濃度を算出するものである、ガス測定方法。
【手続補正２】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項８
【補正方法】変更
【補正の内容】
【請求項８】
請求項６に記載のガス測定装置であって、
前記制御部は、前記第２圧力の下での測定の際に、前記目的成分の吸収ピークの波長についてキャビティリングダウン吸収分光法による測定を実施し、
前記演算処理部は、前記第１圧力の下での結果と前記第２圧力の下での結果とに基づく連立方程式を作成し、該連立方程式を解くことで、前記被測定ガス中の前記目的成分以外の成分による吸収の影響を除去又は軽減した前記目的成分の濃度を算出する、ガス測定装置。
【手続補正３】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項９
【補正方法】変更
【補正の内容】
【請求項９】
請求項６に記載のガス測定装置であって、
前記圧力調整部は、前記測定セル中に前記第２圧力で被測定ガスを封入した状態から、該測定セル中から一部の被測定ガスを外部に強制的に排出することで該測定セル中の被測定ガスの圧力を前記第１圧力に調整するものである、ガス測定装置。
【手続補正４】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項１０
【補正方法】変更
【補正の内容】
【請求項１０】
請求項６に記載のガス測定装置であって、
前記圧力調整部は、前記測定セル中に被測定ガスを供給し前記第１圧力で該被測定ガスを封入した状態から、先に供給されずに残存されていた被測定ガスを該測定セル中に追加供給して封入することで該測定セル中の被測定ガスの圧力を前記第２圧力に調整するものである、ガス測定装置。
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】００２９
【補正方法】変更
【補正の内容】
【００２９】
本発明に係るガス測定方法及びガス測定装置では、第２測定ステップにおける測定で使用されるレーザ光の波長は第１測定ステップにおける測定で使用されるレーザ光の波長と同じでよい。即ち、第１測定ステップと第２測定ステップとで同じ波長のレーザ光を用いてＣＲＤＳ法による測定を実施しても、14ＣＯ2等の低濃度の成分による吸収に対応する信号とこれに重なっている12ＣＯ2、13ＣＯ2等のより高い濃度の成分による吸収に対応する信号とを分離し、低濃度の成分の吸収係数や濃度を算出することができる。
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】００３６
【補正方法】変更
【補正の内容】
【００３６】
即ち、本発明に係るガス測定方法の一つの態様として、
前記第２測定ステップでは、前記目的成分の吸収ピークの波長についてキャビティリングダウン吸収分光法による測定を実施し、
前記演算ステップでは、前記第１測定ステップの結果と前記第２測定ステップの結果とに基づく連立方程式を作成し、該連立方程式を解くことで、前記被測定ガス中の前記目的成分以外の成分による吸収の影響を除去又は軽減した前記目的成分の濃度を算出するようにしてもよい。
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】００３７
【補正方法】変更
【補正の内容】
【００３７】
また本発明に係るガス測定装置の一つの態様として、
前記制御部は、前記第２圧力の下での測定の際に、前記目的成分の吸収ピークの波長についてキャビティリングダウン吸収分光法による測定を実施し、
前記演算処理部は、前記第１圧力の下での結果と前記第２圧力の下での結果とに基づく連立方程式を作成し、該連立方程式を解くことで、前記被測定ガス中の前記目的成分以外の成分による吸収の影響を除去又は軽減した前記目的成分の濃度を算出する構成としてもよい。
【手続補正８】
【補正対象書類名】明細書
【補正対象項目名】００４５
【補正方法】変更
【補正の内容】
【００４５】
或いは本発明に係るガス測定装置において、前記圧力調整部は、前記測定セル中に被測定ガスを供給し第１圧力で該被測定ガスを封入した状態から、先に供給されずに残存されていた被測定ガスを該測定セル中に追加供給して封入することで該測定セル中の被測定ガスの圧力を第２圧力に調整するものである構成としてもよい。
【手続補正９】
【補正対象書類名】明細書
【補正対象項目名】００６７
【補正方法】変更
【補正の内容】
【００６７】
［条件３：14ＣＯ2による吸収を無視できる場合］
条件３の波数では、14ＣＯ2による吸収が無視できる程度に12ＣＯ2や13ＣＯ2の吸収に起因するバックグラウンドが大きい。そのため、ＣＲＤＳによる測定結果には14ＣＯ2による吸収の影響が実質的に現れない。そこで、このときに得られる測定データから求まるリングダウン時間に基づいて14ＣＯ2以外のＣＯ2同位体の吸収係数を算出し、その吸収係数から相対的に高い圧力の下での14ＣＯ2以外のＣＯ2同位体の濃度を算出することができる。そして、算出された14ＣＯ2以外のＣＯ2同位体の濃度を用いることで、14ＣＯ2の吸収が大きい、相対的に低い圧力の下での条件１の波数におけるバックグラウンドに相当する吸収係数を算出することができる。そこで、条件１において得られた測定結果から求まる吸収係数からバックグラウンドに相当する吸収係数を差し引くことで14ＣＯ2のみの吸収係数を求め、この吸収係数から目的成分である14ＣＯ2のみの濃度を算出することができる。
【手続補正１０】
【補正対象書類名】明細書
【補正対象項目名】００７８
【補正方法】変更
【補正の内容】
【００７８】
データ処理部７においてリングダウン時間算出部７２は測定データ格納部７１に保存した高圧力の下での測定データに基づいてリングダウン時間ｔ3を算出する（ステップＳ８）。そして、濃度演算部７３は、その算出結果と演算用既知情報保存部７４に保存されている条件３の下での被測定ガス無し時のリングダウン時間とに基づいて吸収係数を計算し、さらにその吸収係数から濃度を求める（ステップＳ９）。これらの計算方法は従来と同じであり、例えば上記(1)、(2)式を用いればよい。このときには14ＣＯ2による吸収の影響は無視できるので、ステップＳ９において求まるのは被測定ガス中の14ＣＯ2以外のＣＯ2同位体ガスの濃度である。
【手続補正１１】
【補正対象書類名】明細書
【補正対象項目名】００８０
【補正方法】変更
【補正の内容】
【００８０】
ステップＳ９では14ＣＯ2を除くＣＯ2同位体ガスの濃度が得られている。そこで、濃度演算部７３はこの濃度から条件１の圧力及びレーザ波数の下での14ＣＯ2を除くＣＯ2同位体ガスによる吸収係数を計算する。この吸収係数がバックグラウンドであるから、14ＣＯ2を含むＣＯ2同位体ガスによる吸収係数から14ＣＯ2を除くＣＯ2同位体による吸収係数を差し引き、条件１における14ＣＯ2による吸収係数を算出し、この吸収係数から14ＣＯ2のみの濃度を算出する（ステップＳ１２）。そして、その結果を出力部８から出力する。
以上のようにして、本実施例のＣＲＤＳ装置では、バックグラウンドが除去された14ＣＯ2のみの正確な濃度の情報をユーザに提供することができる。
【手続補正１２】
【補正対象書類名】明細書
【補正対象項目名】００８９
【補正方法】変更
【補正の内容】
【００８９】
図６（ａ）は、ガス圧力が１０１３．２５Pa（＝０．０１atm）であるときのＨ2Ｏ同位体であるＤＨＯの吸収ピーク測定を行う場合に得られる、ＣＲＤＳによる吸収スペクトルを計算した結果を示す図である。図６（ｂ）は図６（ａ）中に下向き矢印で示す波数（ＤＨＯ吸収ピークの位置）における吸収の要因の内訳を示す円グラフである。一方、図７（ａ）はガス圧力が図６の１００倍である１０１３２５Pa（＝１atm）であるときのＨ2Ｏ同位体であるＤＨＯの吸収ピーク測定を行う場合に得られる、ＣＲＤＳによる吸収スペクトルを計算した結果を示す図である。図７（ｂ）は図７（ａ）中に下向き矢印で示す波数（図６（ａ）中の下向き矢印と同じ波数）における吸収の要因の内訳を示す円グラフである。
これら計算に際し、温度は３５３Ｋ、被測定ガス中のＨ2Ｏの濃度は０．０３であると想定した。また、被測定ガスに含まれるＨ2Ｏ同位体は全て天然同位体存在比で以て測定セル４０に導入されると想定した。
【国際調査報告】
2020105118.tif 000015
2020105118.tif 000016
2020105118.tif 000017
2020105118.tif 000018