(57)【要約】波数可変レーザ光源１１３と、波長可変レーザ光源から出射されて測定対象ガスを通過したレーザ光の強度を検出する光検出器１１４と、レーザ光を所定の波数範囲で繰り返し掃引するように波数可変レーザ光源１１３に駆動電流を供給するレーザ駆動手段１１２とを備えたガス吸収分光装置において、測定対象ガスの圧力の値又は該圧力に同期して変化する値である、圧力関連値を取得する圧力関連値取得手段１１７と、前記圧力関連値に応じて前記掃引を行う波数範囲を変化させるようレーザ駆動手段１１２を制御する制御手段１３１とを設ける。これにより、測定対象ガスの圧力が変化し且つ高速応答性が求められる状況下においても、低圧から高圧までの広い圧力範囲で高精度な計測を行うことが可能となる。
2019159581.tif 000002
【特許請求の範囲】
【請求項１】
波数可変レーザ光源と、
前記波数可変レーザ光源から出射されて測定対象ガスを通過したレーザ光の強度を検出する光検出器と、
前記レーザ光を所定の波数範囲で繰り返し掃引するように前記波数可変レーザ光源に駆動電流を供給するレーザ駆動手段と、
前記測定対象ガスの圧力の値又は該圧力に同期して変化する値である、圧力関連値を取得する圧力関連値取得手段と、
前記圧力関連値に応じて前記掃引を行う波数範囲を変化させるよう前記レーザ駆動手段を制御する制御手段と、
を有することを特徴とするガス吸収分光装置。
【請求項２】
それぞれ掃引波数幅が異なる掃引波形が規定された複数のテーブルが、前記圧力関連値と対応付けて記憶されたテーブル記憶手段、
を更に有し、
前記制御手段が、前記複数のテーブルのうち前記圧力関連値取得手段から取得された前記圧力関連値に対応するテーブルを、前記テーブル記憶手段から読み出し、該テーブルに従って前記レーザ駆動手段を制御することを特徴とする請求項１に記載のガス吸収分光装置。
【請求項３】
前記光検出器で検出された光強度の変化の曲線を、波数の各点において所定の波数幅の範囲内で近似多項式により近似する多項式近似部と、
前記各点の近似多項式の各項の係数に基づき、前記曲線のゼロ次を含むｎ次微分曲線を作成する微分曲線作成部と、
前記ゼロ次を含むｎ次微分曲線に基づき前記測定対象ガスに含まれる特定ガスの温度、濃度、及び分圧のうちの少なくとも１つを決定する物理量決定手段と、
を更に有し、前記多項式近似部が、前記圧力関連値取得手段で取得された前記圧力関連値に応じて、前記近似を行う波数幅を変化させることを特徴とする請求項１に記載のガス吸収分光装置。
【請求項４】
前記レーザ駆動手段が、前記駆動電流を所定の変調振幅且つ所定の変調周波数で変調するものであって、
前記光検出器による検出信号から前記変調周波数の成分又は前記変調周波数の高調波成分を抽出する復調手段、
を更に有し、
前記制御手段が、更に、前記圧力関連値に応じて前記変調振幅を変化させるよう前記レーザ駆動手段を制御することを特徴とする請求項１に記載のガス吸収分光装置。
【請求項５】
前記測定対象ガスがエンジンの燃焼室内のガスであって、
前記圧力関連値が、前記エンジンのクランク角であることを特徴とする請求項１に記載のガス吸収分光装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、測定対象ガスのレーザ光吸収スペクトルに基づき、該測定対象ガスに含まれる特定ガスの濃度、温度、又は分圧等を測定するガス吸収分光装置に関する。
【背景技術】
【０００２】
レーザを用いたガス吸収分光法としては主に以下の２つの方法が知られている。
(1) ＤＬＡＳ（Direct Laser Absorption Spectroscopy、レーザ直接吸収分光法）
(2) ＷＭＳ（Wavelength Modulated Spectroscopy、波長変調分光法）【０００３】
ＤＬＡＳでは、レーザ光を測定対象ガスに照射し、光検出器によりレーザ光を測定する。ここで、ガスに照射するレーザ光の波長（波数）を特定の値に固定してガスでの吸収を測定する方法と、レーザ光の波長を掃引し、ガスの吸収スペクトルを測定する方法がある。前者の場合には、レーザ光の波長は特定ガスの吸収波長に固定し、その波長での吸光度を測定する。波長を掃引する場合には、レーザ光の波長を特定ガスの吸収波長を含む範囲で変化させてそのスペクトルを測定し、ガスによる吸収ピークの面積を測定する。【０００４】
ＷＭＳは波長掃引型のＤＬＡＳに似ているが、波長掃引に加えて、掃引周期より十分短い周期（すなわち十分高い周波数。ここではｆとする。）で正弦波状に波長を変調する。検出器では、周波数ｆの高調波（一般的に用いられるのは２次高調波）を検出することでＤＬＡＳよりも高い感度で吸収を測定することができる。ＷＭＳでは得られた吸収スペクトルの強度から容易にガス濃度が算出される。【０００５】
特に、産業用のガス吸収分光法としては、感度に優れたＷＭＳが好適とされている。但し、ＷＭＳには以下の理由から、高速測定において精度のよいガス測定が困難であるという問題があった。
１．高速測定を行うためには、掃引周期を短くすると共に高い周波数での波長変調が必要となる。しかしながら、波長可変型のレーザとして一般的に最も普及している注入電流制御型波長可変ダイオードレーザを使用した場合、変調周波数を高めると注入電流に対する波長変化率が低下し、十分な波長変調幅（modulation depth）が得られない。
２．特にＭＨｚを超えるような高速な変調に対しては波長変調幅を正確に測定することは難しく、高速測定において正確な波長変調幅が確定できない。そのため、測定結果から算出されるガス濃度・温度等の情報の不確定性が高くなる。【０００６】
上記の問題を解決するため、本発明者の一部は、特許文献１において新規のガス吸収分光法を提案している（以下これを「改良型ＷＭＳ」とよぶ）。この改良型ＷＭＳでは、上記従来のＷＭＳのようなレーザ光の変調は行わず、波長掃引型のＤＬＡＳと同様に特定ガスの吸収線を含む所定の波長範囲でレーザ光の波長掃引を行う。この光は、測定対象ガスを通過した後、光検出器により受光され、その強度変化が検出される。波長掃引を行う波長範囲は、予め、特定ガスの吸収線の波長を含むように設定されるため、光検出器により検出された光のスペクトルプロファイル（光強度の変化の曲線）には、特定ガスに固有の吸収線の波長を中心とする吸収ピークが現れる。改良型ＷＭＳでは、この吸収ピークを含むスペクトルプロファイルに対し、ＷＭＳ処理に類似の数学的演算を行う。具体的には、各波長ポイントを中心に、ＷＭＳの波長変調幅に相当する区間のスペクトルプロファイルに対してｎ次多項式近似を行い、フーリエ変換の原理に基づいてｎ次多項式の係数を用いて、ＷＭＳ信号振幅を再現する。その原理は次の通りである。【０００７】
一般にＷＭＳ処理では、同期検波して得られるｎ次の高調波のスペクトルプロファイルは、近似的に吸収スペクトルをｎ階微分した波形になることが知られている（非特許文献１：Equation 8）。したがって、波長掃引で得られたスペクトルをｎ階微分すれば、ｎ次の同期検波に相当するスペクトルが得られると考えられる。しかし、ｎ次微分をすると計測データのノイズの影響が大きくなり実用上問題がある。そのため、改良型ＷＭＳでは、高調波信号を求めたい波長を中心としたある範囲に対してｎ次多項式近似を行う。得られる多項式の係数はＷＭＳ処理により得られる高調波信号となる。このとき多項式近似を行う範囲がＷＭＳ処理での変調振幅に相当する。この近似多項式の次数は、高い方がより精度の高い近似を行うことができるが、一般的には１次又は２次多項式近似で十分である。また、ガス吸収以外の光遮断等の光量変化補正処理も行う。【０００８】
このような改良型ＷＭＳでは、光源においては数100kHz以下の波長掃引のみしか行わないため、光源の注入電流に対する発振波長は正確に決定される。そして、その波長情報に基づいて数学的演算によりＷＭＳ処理を行うため、光源駆動電源・光源自体の非線形性の影響も受けず、正確な波長変調幅での高次同期検出が可能となる。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】国際公開第２０１４／１０６９４０号
【非特許文献】
【００１０】
【非特許文献１】Reid, J. and Labrie, D., "Second-harmonic detection with tunable diode lasers-comparison of experiment and theory," Appl. Phys. B 26, 203-210 (1981).
【非特許文献２】福里 克彦, 池田 裕二, 中島 健,「半導体レーザ分光システムを用いたCO2ガスの計測 (第2報)」,日本機械学会論文集 B編, 2002, 68, 2901-2907.
【非特許文献３】G. B. Rieker, J. B. Jeffries, and R. K. Hanson, "Calibration-free wavelength modulation spectroscopy for measurements of gas temperature and concentration in harsh environments," Appl. Opt., submitted 2009.
【発明の概要】
【発明が解決しようとする課題】
【００１１】
ところで、ガス吸収分光法で取得される吸収線の中心波長やピーク幅は、測定対象ガスの圧力値に依存する。一例として、ＨＩＴＲＡＮデータベースから取得した、Ｈ２Ｏ（濃度 1.2％、光路長 1.0cm、温度 400K）の波長1348.4nm（波数7416.0cm-1）付近における吸収線の圧力特性を図１１に示す。同図から明らかなように、ガス圧が高くなるにつれて、吸収線ピークの中心位置がシフトし、ピーク幅も大きく広がっていく。【００１２】
こうした吸収線の圧力特性により、従来のガス吸収分光法では、測定対象ガスの圧力が変化し且つ高速応答性が求められる状況において、低圧又は高圧のいずれかにおける測定精度が悪化するという問題があった。具体的には、例えば、高圧下で広がった吸収ピークを確実に測定するためには、掃引波数（波長）範囲を広く設定する必要があるが、そのようにすると、検出信号のサンプリングレートの限界により、低圧下における幅狭のピークを高いＳ／Ｎで測定することができなくなる。一方、低圧下における幅狭のピークを高いＳ／Ｎで測定できるように掃引波数範囲を狭く設定した場合、高圧下で広がった吸収ピークについては、その一部分しか測定できなくなってＳ／Ｎが悪化したり、ピークシフトによりピークの中心が掃引波数範囲から外れてしまい、そもそも所望する物理量（例えば、特定ガスの濃度等）を算出できなくなったりするおそれがある。【００１３】
本発明は上記の点に鑑みて成されたものであり、その目的とするところは、測定対象ガスの圧力が変化し且つ高速応答性が求められる状況において、低圧から高圧までの広い圧力範囲で高精度な計測を行うことができるガス吸収分光装置を提供することにある。
【課題を解決するための手段】
【００１４】
上記課題を解決するために成された本発明に係るガス吸収分光装置は、
波数可変レーザ光源と、
前記波数可変レーザ光源から出射されて測定対象ガスを通過したレーザ光の強度を検出する光検出器と、
前記レーザ光を所定の波数範囲で繰り返し掃引するように前記波数可変レーザ光源に駆動電流を供給するレーザ駆動手段と、
前記測定対象ガスの圧力の値又は該圧力に同期して変化する値である、圧力関連値を取得する圧力関連値取得手段と、
前記圧力関連値に応じて前記掃引を行う波数範囲を変化させるよう前記レーザ駆動手段を制御する制御手段と、
を有することを特徴としている。【００１５】
なお、ここで言う「波数」は「波長」と一義的に対応するものであり、「波長」を用いて同様の構成を組み立てることももちろん可能である。【００１６】
上記構成によれば、圧力関連値取得手段によって取得される圧力関連値に基づき、測定対象ガスの圧力が大きくなるにつれて掃引波数範囲が広くなるようにレーザ駆動手段を制御することにより、測定対象ガスの圧力が高圧から低圧まで大きく変化する状況下においても、常にそのときのガス圧に応じた適切な波数範囲で波数掃引を行うことができる。そのため、低圧下及び高圧下のいずれにおいても高いＳ／Ｎで測定を行うことができる。【００１７】
なお、本発明に係るガス吸収分光装置及び方法は、例えば、自動車産業において、ガス濃度をはじめ温度、圧力の非接触、高速測定に適用可能であり、その他、プラント炉内の燃焼ガスのような高温・高圧環境におけるガス計測など多岐の分野において応用可能である。ここで例えば、ピストン式の内燃機関又は外燃機関の燃焼室内のガスを前記測定対象ガスとする場合、該燃焼室内のガス圧を圧力センサで直接測定して得られた値を前記圧力関連値とするほか、例えば、燃焼室内のガス圧に同期して変化する値であるクランク角を前記圧力関連値として利用することができる。【００１８】
上記本発明に係るガス吸収分光装置は、
それぞれ掃引波数幅が異なる掃引波形が規定された複数のテーブルが、前記圧力関連値と対応付けて記憶されたテーブル記憶手段、
を更に有し、
前記制御手段が、前記複数のテーブルのうち前記圧力関連値取得手段から取得された前記圧力関連値に対応するテーブルを、前記テーブル記憶手段から読み出し、該テーブルに従って前記レーザ駆動手段を制御するものとすることができる。【００１９】
本発明は、例えば、上述の改良型ＷＭＳによる測定を行うガス吸収分光装置に適用することができる。改良型ＷＭＳでは、レーザ光の変調を行う代わりに、光検出器によって検出された光強度の変化の曲線（スペクトルプロファイル）を、波数の各点においてＷＭＳの変調波数幅に相当する波数幅で多項式近似する。なお、「ＷＭＳの変調波数幅」は、吸収線の半値半幅の２．２倍とした場合に最大のＳ／Ｎが得られることが知られているが、上述の通り、吸収線のピーク幅は測定対象ガスの圧力に応じて変化する。そこで、本発明を改良型ＷＭＳによる測定を行うガス吸収分光装置に適用する場合には、測定対象ガスの圧力変化に応じて、掃引波数幅だけでなく前記多項式近似を行う際の波数幅も変化させることが望ましい。【００２０】
すなわち、上記本発明に係るガス吸収分光装置は、
前記光検出器で検出された光強度の変化の曲線を、波数の各点において所定の波数幅の範囲内で近似多項式により近似する多項式近似部と、
前記各点の近似多項式の各項の係数に基づき、前記曲線のゼロ次を含むｎ次微分曲線を作成する微分曲線作成部と、
前記ゼロ次を含むｎ次微分曲線に基づき前記測定対象ガスに含まれる特定ガスの温度、濃度、及び分圧のうちの少なくとも１つを決定する物理量決定手段と、
を更に有し、前記多項式近似部が、前記圧力関連値取得手段で取得された前記圧力関連値に応じて、前記近似を行う波数幅を変化させるものとすることができる。【００２１】
ここで「特定ガス」とは、測定者等によって決められる任意の成分であって、例えば酸素、水蒸気、二酸化炭素、又は一酸化炭素等である。【００２２】
また、本発明は、上述のＷＭＳによる測定を行うガス吸収分光装置に適用することもできる。ＷＭＳでは、レーザの発振波数が、所定の周波数で変調される。このとき、最大のＳ／Ｎが得られる変調波数幅は測定対象ガスの圧力に応じて変化する。そのため、本発明をＷＭＳによる測定を行うガス吸収分光装置に適用する場合には、測定対象ガスの圧力変化に応じて、掃引波数幅だけでなく変調波数幅も変化させることが望ましい。【００２３】
すなわち、上記本発明に係るガス吸収分光装置は、
前記レーザ駆動手段が、前記駆動電流を所定の変調振幅且つ所定の変調周波数で変調するものであって、
前記光検出器による検出信号から前記変調周波数の成分又は前記変調周波数の高調波成分を抽出する復調手段、
を更に有し、
前記制御手段が、更に、前記圧力関連値に応じて前記変調振幅を変化させるよう前記レーザ駆動手段を制御するものとすることができる。
【発明の効果】
【００２４】
以上の通り、本発明に係るガス吸収分光装置によれば、測定対象ガスの圧力が変化し且つ高速応答性が求められる状況においても、低圧から高圧までの広い圧力範囲で高精度な計測を行うことが可能となる。
【図面の簡単な説明】
【００２５】
【図１】本発明の一実施形態に係るガス吸収分光装置の概略構成図。
【図２】同実施形態における制御部の動作を示すフローチャート。
【図３】本発明の別の実施形態に係るガス吸収分光装置の概略構成図。
【図４】同実施形態における解析部の動作を示すフローチャート。
【図５】本発明の更に別の実施形態に係るガス吸収分光装置の概略構成図。
【図６】同実施形態における制御部の動作を示すフローチャート。
【図７】エンジンの燃焼室内のガスを測定対象ガスとする場合の構成例を示す図。
【図８】改良型ＷＭＳにおいてスペクトルプロファイルを多項式で表す方法を模式的に示す説明図。
【図９】ＷＭＳにおけるレーザ駆動信号を示す波形図であって、(a)が掃引信号の波形であり、(b)が変調信号の波形である。
【図１０】ＷＭＳにおけるレーザ出力波形を示す図。
【図１１】Ｈ２Ｏの波長1348.4nm（波数7416.0cm-1）の吸収線の圧力特性を示す図。
【発明を実施するための形態】
【００２６】
＜実施形態１＞
以下、本発明を実施するための形態について図面を参照しつつ説明する。本発明の一実施形態によるガス吸収分光装置の概略構成を図１に示す。このガス吸収分光装置は、改良型ＷＭＳ（すなわち特許文献１に記載の方法）による測定を行うものであり、測定対象ガスが通過する（又は測定対象ガスが収容された）ガスセル１１０と、ガスセル１１０を挟んで対向配置されたレーザ光源１１３及び光検出器１１４と、レーザ光源１１３に駆動電流を注入するレーザ駆動部１１２と、レーザ駆動部１１２に所定の掃引信号を入力する掃引信号発生部１１１と、前記各部の制御及び光検出器１１４からの出力の解析を行う制御／解析部１２０と、を備えている。更に、ガスセル１１０には、本実施形態における特徴的構成の１つである圧力センサ１１７が設けられている。
【００２７】
制御／解析部１２０は、制御部１３０、解析部１４０、及び記憶部１５０を含んでいる。制御部１３０は、掃引信号発生部１１１を制御するレーザ制御部１３１を有し、解析部１４０は、光検出器１１４の検出信号を処理するための多項式近似部１４１及び微分曲線作成部１４２と、該処理後の信号から所望の物理量を算出する物理量決定部１４３とを備えている。また、記憶部１５０にはテーブル記憶部１５１が設けられている（詳細は後述する）。
【００２８】
制御／解析部１２０の機能は、ＣＰＵ、メモリ、及び大容量記憶媒体（ハードディスク等）などを具備するコンピュータによって実現される。該コンピュータは、ガス吸収分光装置の本体に内蔵された専用コンピュータであってもよいが、典型的には、パーソナルコンピュータ等を使用する。前記コンピュータには、予め所定のプログラムがインストールされており、ＣＰＵがこのプログラムを実行することによって、本システムにおける制御部１３０及び解析部１４０の機能がソフトウェア的に具現化される。また、記憶部１５０の機能は前記大容量記憶媒体によって実現される。
【００２９】
レーザ光源１１３は、波数が可変のものであり、その波数が、特定ガスの吸収線の波数を含む所定の波数範囲で掃引される。波数可変のレーザ光源１１３として用いられる半導体レーザダイオード（以下「ＬＤ」と略す）の発振波数は、注入される電流（注入電流）の大きさに依存するから、レーザ光の波数掃引は該注入電流を掃引することにより実施される。具体的には、レーザ制御部１３１の制御の下に、掃引信号発生部１１１が鋸歯状パターンを有する信号（図１参照）を発生し、これをレーザ駆動部１１２に送出する。レーザ駆動部１１２は、該信号に従って鋸歯状に変化する注入電流をＬＤから成るレーザ光源１１３に注入する。これにより、レーザ光源１１３の発振波数が、所定の波数範囲で繰り返し掃引されることとなる。
【００３０】
まず、本実施形態によるガス吸収分光装置の基本的な動作について説明する。本実施形態のガス吸収分光装置によって測定対象ガス中の特定ガスの濃度、温度、又は分圧等を測定する際には、制御部１３０の制御の下に、レーザ光源１１３より所定の最高波数のレーザ光が放射され、その波数が順次変化されて最低波数まで掃引される。なお、上述のＷＭＳでは、波数掃引に加えて、掃引周期より十分短い周期で正弦波状に波数を変調するが、本実施形態に係る装置は、改良型のＷＭＳによる測定を行うものであるため、このような変調は行わない。レーザ光源１１３からの光はガスセル１１０中の測定対象ガスを通過し、その際に、特定ガスの吸収線の波数において吸収を受ける。測定対象ガスを通過したレーザ光は、光検出器１１４でその強度が検出される。光検出器１１４から出力される、光強度を表す電気信号は、アンプ１１５を経てＡ／Ｄ変換器１１６でデジタルサンプリングされ、解析部１４０に送られる。この電気信号の変化がスペクトルプロファイルとなる。解析部１４０は、このスペクトルプロファイルを表すデータに基づき、所定の数学的演算を行う。
【００３１】
解析部１４０が行う、前記数学的演算について説明する。ここでは、まず、上記スペクトルプロファイルの、波数軸の各点νを中心とする幅２ａ’の範囲［ν−ａ’＜ν＜ν＋ａ’］について以下の式（１）
【数１】
2019159581.tif 000003
で示される多項式で表されると考える（なお、本願の本文中では、電子出願の制約により、νに付された上線を下線で表現している）。これを図８に模式的に示す。
式（１）のｎ次微分を求めると
【数２】
2019159581.tif 000004
となる。ここで、一般にはＷＭＳ処理で同期検波して得られるｎ次の高調波のスペクトルプロファイルは、近似的に次式で示されることが知られている（非特許文献１：Equation 8）。
【数３】
2019159581.tif 000005
したがって、式(２)、(３)より
【数４】
2019159581.tif 000006
となる。したがって、上記スペクトルプロファイルにおいて波数νに対するＷＭＳ信号を算出するためには、［ν−ａ’＜ν＜ν＋ａ’］の波数の範囲を最小二乗法等によりフィッティングし、係数ｂ０、ｂ１、ｂ２、ｂ３…を求める。νを逐次変化させてフィッティングにより求めた係数ｂ１とｂ２のプロファイルが、１ｆと２ｆのＷＭＳプロファイルに相当するものとなる。なお、フィッティングの範囲を表すａ’はＷＭＳの変調振幅ａ（すなわち変調波数幅）に相当する値となる。
【００３２】
具体的には、多項式近似部１４１が、上記スペクトルプロファイルの波数νを中心とする幅２ａ’の範囲［ν−ａ’＜ν＜ν＋ａ’］を最小二乗法等によってフィッティングして係数ｂ０、ｂ１、ｂ２、ｂ３…を求める。
【００３３】
続いて、微分曲線作成部１４２が、前記の波数νを逐次変化させてフィッティングを行うことにより求められた係数ｂ０、ｂ１、ｂ２、ｂ３…を、波数νに対してプロットすることにより各係数のプロファイル（すなわちゼロ次を含む高次微分曲線）を作成する。ここで、係数ｂ１のプロファイルと係数ｂ２のプロファイルが、それぞれＷＭＳにおける一次同期検波プロファイルと二次同期検波プロファイルに相当するものとなる。
【００３４】
続いて、上記処理によって作成された（ゼロ次を含む）高次微分曲線に基づき、物理量決定部１４３が、測定対象ガス中の特定ガスの濃度、分圧、又は温度などを算出する。例えば、特定ガスの濃度はゼロ次微分曲線の吸収ピークの面積から算出することができる。また、２次微分曲線のピーク高さからも特定ガスの濃度を算出することができる。測定対象ガスの分圧Ｐは、ゼロ次微分曲線の吸収ピークの半値幅αＬと次のような関係を有することが知られている。
【数５】
2019159581.tif 000007
ここで、αＬ０は気圧Ｐ０、温度Ｔ０におけるローレンツ広がりの半値半幅、Ｐ０は基準時の測定対象ガスの圧力、Ｔは測定時の測定対象ガスの温度、Ｔ０は基準時の測定対象ガスの温度、γはローレンツ幅の温度依存性を表す定数である。この式により、特定ガスの分圧を求めることができる。また、特定ガスの温度については、温度が変化するに伴い、２つの吸収ピークの大きさの比が変化することが知られており、その関係を用いることにより、測定対象ガスの温度を検出することができる（非特許文献２）。
【００３５】
次に、透過光強度の正規化処理について述べる。ガス吸収分光法では、ガスセルで使用している光学部品の汚れや劣悪環境下での振動による光軸変化に伴い光強度が変化することが実用上の課題の１つとなっている。したがって、光強度の補正処理が必要となるが、補正方法の１つとして、同期検波された２ｆ信号を１ｆ信号で除算する正規化処理（非特許文献３）が知られている。しかし、この方法では、レーザ光を変調する必要がある上、同期検波回路も１ｆ用と２ｆ用の２種類を用意する必要がある。
【００３６】
一方、上記のような多項式近似を用いたＷＭＳ相当処理では、レーザ光の変調や同期検波回路も不要であり、近似をする際に同時に１ｆ及び２ｆの検波信号が算出できるため、非常に簡便に正規化処理が可能である。以下に詳細を示す。
【００３７】
測定対象ガスに入射する光強度をＩ０とすれば、検出される光強度はＳ（ν）＝ＧＩ０τ（ν）となる。Ｇは各光学部品による光強度の低下（及び変動）と検出した光強度に対する電気的ゲインを示す。よって実際の装置では、数学的演算によるＷＭＳ処理をＳ（ν）に対して適用し、次式のようになる。
【数６】
2019159581.tif 000008
したがって、ここで得られる係数は
【数７】
2019159581.tif 000009
となる。よって、光強度の変動に依存せず透過スペクトルのみに依存する値を得るためには、次のようにｂ２’（２ｆ信号）をｂ１’（１ｆ信号）又はｂ０’で除算すればよい。
【数８】
2019159581.tif 000010
なお、吸光が少ない場合はｂ０〜１（ｂ０が１に近い）であるため、式（８ｂ）で示すような近似が可能である。以上より、光強度に依存しないロバストなガス計測が可能となる。
【００３８】
続いて、本実施形態のガス吸収分光装置における特徴的な動作について説明する。上述のように、特定ガスの吸収線のピーク幅は測定対象ガスの圧力に応じて変化する。そこで、本実施形態に係るガス吸収分光装置では、高速応答性が求められる状況において、高圧下における幅広の吸収ピークも低圧化における幅狭の吸収ピークも適切に検出できるよう、圧力センサ１１７によって検出される測定対象ガスの圧力が高いときには掃引波数範囲を相対的に広くし、前記圧力が低いときには掃引波数範囲を相対的に狭くするようレーザ光源１１３を制御する。
【００３９】
上記のような制御を実現するため、テーブル記憶部１５１には、測定対象ガスの圧力範囲毎にそれぞれ掃引波数幅が異なる掃引波形が規定された複数のテーブルが記憶されている。これら複数のテーブル（以下、これをテーブルセットとよぶ）は、ガス吸収分光装置で測定する特定ガスの種類毎に用意し、例えば、ユーザが測定開始前に特定ガスの種類を設定した時点で、該特定ガスの種類に応じたテーブルセットが自動的に選択される。
【００４０】
測定対象ガスの測定中に、圧力センサ１１７から出力される、測定対象ガスの圧力を表す電気信号は、アンプ１１８を経てＡ／Ｄ変換器１１９でデジタルサンプリングされ、制御部１３０に送られる。以下、制御部１３０の動作について図２のフローチャートを参照しつつ説明する。
【００４１】
制御部１３０が前記測定対象ガスの圧力（ガス圧）を表す信号を受け取ると（ステップＳ１１）、制御部１３０に設けられたレーザ制御部１３１によって、該圧力に対応するテーブルがテーブル記憶部１５１から読み出される（ステップＳ１２）。そして、該テーブルに記述された掃引波形の情報（例えば、波数掃引の各時刻における注入電流の値）に基づいて、レーザ制御部１３１が、掃引信号発生部１１１を制御する（ステップＳ１３）。これにより、圧力センサ１１７で検出されたガス圧に応じた適切な掃引波数幅によるレーザ光の波数掃引が実行される。このような処理（ステップＳ１１〜Ｓ１３）は、ユーザから制御／解析部１２０に対して測定終了の指示が入力されるまで（すなわち、ステップＳ１４でＹｅｓになるまで）一定の時間間隔で繰り返し実行される。
【００４２】
以上の処理により、本実施形態に係るガス吸収分光装置によれば、測定対象ガスの圧力が変化し且つ高速応答性が求められる状況下でも、その時々のガス圧に応じた適切な掃引波数幅でのレーザ光の波数掃引が実行されるため、常に高いＳ／Ｎで測定を行うことができる。
【００４３】
本発明の効果を確認するために行ったシミュレーションについて説明する。ここでは、Ａ／Ｄ変換器１１６に入力される光検出器１１４からの電圧信号が1348nmのH2O吸収プロファイルがのった鋸歯状の信号であって、その鋸歯状の電圧値が1〜4V、ノイズ幅が500uVrmsであると仮定した。具体的には、HITRAN2012に基づく1348nmの吸収線にそれぞれ異なるホワイトノイズを足し合わせて成る32種類の信号を用意し、Ａ／Ｄ変換器１１６が一掃引において４００点のデジタル値を取得すると仮定した場合に、該32種類のデータから算出される「二次係数／ゼロ次係数」の平均をＳ（＝Signal）、その分散をＮ（＝Noise）としてＳ／Ｎの計算を行った。その結果、低圧（1atm）において、変調幅を0.15cm-1とし、掃引幅2cm-1で波数掃引を行った場合は、Ｓ／Ｎ＝２０であった。一方、低圧（1atm）において、変調幅を同じく0.15cm-1とし、掃引幅を0.7cm-1とした場合、Ｓ／Ｎ＝３３であった。このことから、低圧下においては掃引幅を相対的に狭くすることによって、より高いＳ／Ｎが達成されることが確かめられた。
【００４４】
＜実施形態２＞
以下、本発明に係るガス吸収分光装置の他の実施形態について、図３及び図４を参照しつつ説明する。図３は、本実施形態に係るガス吸収分光装置の概略構成を示す図であり、図４は、本実施形態に係るガス吸収分光装置における解析部の動作を示すフローチャートである。本実施形態において、図１で説明したものと同一又は対応する構成要素については下二桁が共通する符号を付し、適宜説明を省略する。
【００４５】
本実施形態に係るガス吸収分光装置は、測定対象ガスの圧力変化に応じて、レーザ光の掃引波数幅だけでなく、解析部２４０にてスペクトルプロファイルを多項式近似する際のフィッティング幅をも変化させるものとなっている。そのため、本実施形態に係るガス吸収分光装置では、圧力センサ２１７からの検出信号が制御部２３０だけでなく解析部２４０にも入力される。また、テーブル記憶部２５１には、実施形態１と同様に、測定対象ガスの圧力範囲毎に複数のテーブル（以下、これをテーブルセットとよぶ）が記憶されており、各テーブルには、その圧力範囲での測定に適した掃引波形の情報に加えて、その圧力範囲での測定で得られたスペクトルプロファイルの多項式近似に好適なフィッティング幅の情報が記載されている。このようなテーブルセットは、ガス吸収分光装置で測定する特定ガスの種類毎に用意し、例えば、ユーザが測定開始前に特定ガスの種類を設定した時点で、該特定ガスの種類に応じたテーブルセットが自動的に選択される。
【００４６】
本実施形態において、測定対象ガスの測定中に圧力センサ２１７から出力される、測定対象ガスの圧力を表す電気信号は、アンプ２１８を経てＡ／Ｄ変換器２１９でデジタルサンプリングされた上で、制御部２３０及び解析部２４０へと送られる。以下、このときの解析部２４０の動作について図４のフローチャートを参照しつつ説明する（なお、制御部２３０の動作は、図２のフローチャートで示したものと同様であるため説明を省略する）。
【００４７】
前記測定対象ガスの圧力（ガス圧）を表す信号は、解析部２４０に設けられた多項式近似部２４１に送られる（ステップＳ２１）。該信号を受け取った多項式近似部２４１は、該ガス圧に対応するテーブルをテーブル記憶部２５１から読み出し（ステップＳ２２）、該テーブルに記載されたフィッティング幅でスペクトルプロファイルの多項式近似を行う（ステップＳ２３）。そして、その結果に基づいて、微分曲線作成部２４２による微分曲線の作成（ステップＳ２４）、及び物理量決定部２４３における特定ガスの濃度、温度、又は分圧等の算出（ステップＳ２５）が行われる。このような処理（ステップＳ２１〜Ｓ２５）は、ユーザから測定終了の指示が入力されるまで（すなわち、ステップＳ２６でＹｅｓになるまで）一定の時間間隔で繰り返し実行される。
【００４８】
上述の通り、改良型ＷＭＳでは、レーザ光の変調を行う代わりに、光検出器によって検出された光強度の変化の曲線（スペクトルプロファイル）を、波数の各点においてＷＭＳの変調波数幅に相当する波数幅の範囲内で多項式近似する。このとき、最大のＳ／Ｎが得られるＷＭＳの変調波数幅は、測定対象ガスの圧力に応じて変化する。したがって、上記のように圧力センサ２１７で取得されるガス圧に応じて、前記多項式近似を行う波数幅（フィッティング幅）を随時変化させることにより、測定対象ガスの圧力が変化し且つ高速応答性が求められる状況下でも、常に高いＳ／Ｎを達成することができる。
【００４９】
＜実施形態３＞
以下、本発明に係るガス吸収分光装置の更に別の実施形態について、図５及び図６を参照しつつ説明する。図５は、本実施形態に係るガス吸収分光装置の概略構成を示す図であり、図６は、本実施形態に係るガス吸収分光装置における制御部の動作を示すフローチャートである。本実施形態において、図１で説明したものと同一又は対応する構成要素については下二桁が共通する符号を付し、適宜説明を省略する。
【００５０】
このガス吸収分光装置は、ＷＭＳによる測定を行うものであり、実施形態１及び２と同様の鋸歯状の掃引信号を発生する掃引信号発生部３１１に加えて、前記鋸歯状の掃引信号よりも周波数の高い正弦波である変調信号を発生する変調信号発生部３６１を備えている。制御部３３０の制御の下に掃引信号発生部３１１及び変調信号発生部３６１でそれぞれ生成された掃引信号（図９(a)）及び変調信号（図９(b)）は、加算器３６２において可算され、レーザ駆動部３１２へ送出される。レーザ駆動部３１２は、該信号にしたがって変調された鋸歯状の注入電流をレーザ光源３１３に注入する。これにより、レーザ光源３１３からの出力は、図１０に示すように、所定の掃引波数幅Ａで繰り返し掃引され且つ所定の変調波数幅ａで変調されたものとなる。
【００５１】
制御／解析部３２０は、制御部３３０、解析部３４０、及び記憶部３５０を有している。制御部３３０には、レーザ制御部３３１が設けられ、解析部３４０には、光検出器３１４の検出信号を復調する復調部３４４と、該復調後の信号から所望の物理量を算出する物理量決定部３４５が設けられている。記憶部３５０には、測定対象ガスの圧力範囲毎に複数のテーブル（テーブルセット）を記憶したテーブル記憶部３５１が設けられており、各テーブルには、その圧力範囲での測定に適した掃引波形及び変調波形の情報が記載されている。上記のようなテーブルセットは、ガス吸収分光装置で測定する特定ガスの種類毎に用意され、例えば、ユーザが測定開始前に特定ガスの種類を設定した時点で、該特定ガスの種類に応じたテーブルセットが自動的に選択される。
【００５２】
これらの制御部３３０及び解析部３４０は、ＣＰＵ、メモリ、及び大容量記憶装置等を備えたコンピュータに、専用のソフトウェアをインストールして実行することによってソフトウェア的に実現される機能手段である。また、記憶部３５０の機能は前記大容量記憶媒体によって実現される。
【００５３】
本実施形態において、測定対象ガスの測定中に圧力センサ３１７から出力される、測定対象ガスの圧力を表す電気信号は、アンプ３１８を経てＡ／Ｄ変換器３１９でデジタルサンプリングされ、制御部３３０に送られる。以下、制御部３３０の動作について図６のフローチャートを参照しつつ説明する。
【００５４】
制御部３３０が前記測定対象ガスの圧力（ガス圧）を表す信号を受け取ると（ステップＳ３１）、制御部３３０のレーザ制御部３３１によって、前記圧力に対応するテーブルがテーブル記憶部３５１から読み出される（ステップＳ３２）。そして、該テーブルに記述された掃引波形及び変調波形の情報に従って、レーザ制御部３３１が、掃引信号発生部３１１及び変調信号発生部３６１を制御する（ステップＳ３３）。これにより、圧力センサ３１７で検出されたガス圧に応じた適切な掃引波数幅によるレーザ光の波数掃引、及び該ガス圧に応じた適切な変調波数幅によるレーザ光の変調が実行される。このような処理（ステップＳ３１〜Ｓ３３）は、ユーザから測定終了の指示が入力されるまで（すなわち、ステップＳ３４でＹｅｓになるまで）一定の時間間隔で繰り返し実行される。
【００５５】
このように、本実施形態に係るガス吸収分光装置では、圧力センサ３１７で取得されるガス圧に応じてレーザ光の掃引波数幅及び変調波数幅が随時変化させることにより、測定対象ガスの圧力が変化し且つ高速応答性が求められる状況下でも、常に高いＳ／Ｎを達成することが可能となる。
【００５６】
なお、本実施形態では、ガス圧に応じて掃引波数幅及び変調波数幅を変化させるものとしたが、掃引波数幅又は変調波数幅のいずれか一方のみをガス圧に応じて変化させるようにしてもよい。
【００５７】
以上、本発明を実施するための形態について具体例を挙げて説明を行ったが、本発明は上記実施形態に限定されるものではなく、本発明の趣旨の範囲で適宜変更が許容される。
【００５８】
例えば、実施形態１〜３では、ガスセル内を通過する（又はガスセル内に収容された）ガスを測定対象ガスとしたが、これに代えて、内燃機関や外燃機関の燃焼室内のガスを測定対象ガスとしてもよい。この場合、燃焼室内に圧力センサを配置し、該圧力センサで取得したガス圧の値に応じて、測定対象ガスに照射するレーザ光の掃引波数幅、該レーザ光の変調振幅、又はスペクトルプロファイルのフィッティング幅（以下、これらを総称して「掃引波数幅等」とよぶ）を変化させる構成としてもよいが、前記内燃機関又は外燃機関がピストン式のものである場合は、このような圧力センサは設けず、ピストンの上下動に伴って回転するクランクシャフトの回転角度（すなわちクランク角）に応じて掃引波数幅等を変化させる構成とすることもできる。この場合、クランク角が本発明における圧力関連値に相当する。
【００５９】
図７に、本発明に係るガス吸収分光装置をピストン式のエンジンに適用する場合の構成例を示す。エンジン４７０は、シリンダ４７１と、シリンダ４７１内で摺動可能なピストン４７２を有し、シリンダ４７１の内部空間とピストン４７２によって燃焼室４７３が形成されている。シリンダ４７１の周面にはレンズ等から成る２つの光学窓４７４ａ、４７４ｂが対向配置されており、レーザ光源４１３から出射したレーザ光は、一方の光学窓４７４ａから燃焼室４７３に進入した後、他方の光学窓４７４ｂから外に出て光検出器４１４によって受光される。ピストン４７２は、コンロッド４７５を介してクランクシャフト４７６に連結されており、クランクシャフト４７６には、外周に複数の突起４７７ａを備えたタイミングロータ４７７が挿通されている。タイミングロータ４７７はクランクシャフト４７６の回転に伴って回転するよう構成されており、タイミングロータ４７７の外周近傍にはクランク角センサ４７８が配置されている。クランク角センサ４７８は、タイミングロータ４７７の突起４７７ａを電磁気的又は光学的に検知してパルス状の信号（クランク信号）を出力する。
【００６０】
クランク角センサ４７８からのクランク信号（及び光検出器４１４からの検出信号）は、図示しない制御／解析部に送られる。なお、タイミングロータ４７７の外周の一部には、突起４７７ａが設けられていない領域（欠歯部）４７７ｂがあるため、これに対応して前記クランク信号にもパルスが出現しない時間帯が周期的に現れる。そのため、制御／解析部では、この時間帯を基準としてクランク信号上のパルスをカウントすることにより、現在のクランク角を特定することができる。このクランク角は燃焼室４７３内のガス圧と同期して変化するため、該クランク角に基づいて燃焼室４７３内のガス圧を推定することができる。
【００６１】
なお、制御／解析部としては、例えば、実施形態１〜３における制御／解析部１２０、２２０、又は３２０のいずれかと同様の構成を採用することができる。但し、いずれの場合も、圧力センサ１１７、２１７、又は３１７からの検出信号に代えて、クランク角センサ４７８からの信号（クランク信号）を利用するものとし、該クランク信号に基づいて、掃引波数幅等を変化させるものとする。この場合、上記のようにクランク信号のカウント数から燃焼室４７３内のガス圧を推定し、該ガス圧に応じて掃引波数幅等を変化させるようにしてもよいが、予め、各カウント数における適切な掃引波数幅等を記載した複数のテーブルをテーブル記憶部１５１、２５１、又は３５１に記憶させておくようにしてもよい。この場合、測定中の各時点におけるクランク信号から求められたカウント数に応じたテーブルを、レーザ制御部１３１、２３１、又は３３１がテーブル記憶部１５１、２５１、又は３５１から読み出し、該テーブルに基づいて掃引波数幅等を制御する。
【符号の説明】
【００６２】
１１１、２１１、３１１…掃引信号発生部
１１２、２１２、３１２…レーザ駆動部
１１３、２１３、３１３、４１３…レーザ光源
１１４、２１４、３１４、４１４…光検出器
１１７、２１７、３１７…圧力センサ
１２０、２２０、３２０…制御／解析部
１３０、２３０、３３０…制御部
１３１、２３１、３３１…レーザ制御部
１４０、２４０、３４０…解析部
１４１、２４１…多項式近似部
１４２、２４２…微分曲線作成部
１４３、２４３、３４５…物理量決定部
１５０、２５０、３５０…記憶部
１５１、２５１、３５１…テーブル記憶部
１１０、２１０、３１０…ガスセル
３４４…復調部
３６１…変調信号発生部
３６２…加算器
４７０…エンジン
４７１…シリンダ
４７２…ピストン
４７３…燃焼室
４７４ａ、４７４ｂ…光学窓
４７６…クランクシャフト
４７７…タイミングロータ
４７７ａ…突起
４７８…クランク角センサ
【図１】
2019159581.tif 000011
【図２】
2019159581.tif 000012
【図３】
2019159581.tif 000013
【図４】
2019159581.tif 000014
【図５】
2019159581.tif 000015
【図６】
2019159581.tif 000016
【図７】
2019159581.tif 000017
【図８】
2019159581.tif 000018
【図９】
2019159581.tif 000019
【図１０】
2019159581.tif 000020
【図１１】
2019159581.tif 000021
【手続補正書】
【提出日】令和2年6月26日(2020.6.26)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
波数可変レーザ光源と、
前記波数可変レーザ光源から出射されて測定対象ガスを通過したレーザ光の強度を検出する光検出器と、
前記レーザ光を所定の波数範囲で繰り返し掃引するように前記波数可変レーザ光源に駆動電流を供給するレーザ駆動手段と、
前記測定対象ガスの圧力の値又は該圧力に同期して変化する値である、圧力関連値を取得する圧力関連値取得手段と、
前記圧力関連値に基づき、前記圧力が大きくなるにつれて前記掃引を行う波数範囲が広くなるように前記レーザ駆動手段を制御する制御手段と、
を有することを特徴とするガス吸収分光装置。
【請求項２】
それぞれ掃引波数幅が異なる掃引波形が規定された複数のテーブルが、前記圧力関連値と対応付けて記憶されたテーブル記憶手段、
を更に有し、
前記制御手段が、前記複数のテーブルのうち前記圧力関連値取得手段から取得された前記圧力関連値に対応するテーブルを、前記テーブル記憶手段から読み出し、該テーブルに従って前記レーザ駆動手段を制御することを特徴とする請求項１に記載のガス吸収分光装置。
【請求項３】
前記光検出器で検出された光強度の変化の曲線を、波数の各点において所定の波数幅の範囲内で近似多項式により近似する多項式近似部と、
前記各点の近似多項式の各項の係数に基づき、前記曲線のゼロ次を含むｎ次微分曲線を作成する微分曲線作成部と、
前記ゼロ次を含むｎ次微分曲線に基づき前記測定対象ガスに含まれる特定ガスの温度、濃度、及び分圧のうちの少なくとも１つを決定する物理量決定手段と、
を更に有し、前記多項式近似部が、前記圧力関連値取得手段で取得された前記圧力関連値に応じて、前記近似を行う波数幅を変化させることを特徴とする請求項１に記載のガス吸収分光装置。
【請求項４】
前記レーザ駆動手段が、前記駆動電流を所定の変調振幅且つ所定の変調周波数で変調するものであって、
前記光検出器による検出信号から前記変調周波数の成分又は前記変調周波数の高調波成分を抽出する復調手段、
を更に有し、
前記制御手段が、更に、前記圧力関連値に応じて前記変調振幅を変化させるよう前記レーザ駆動手段を制御することを特徴とする請求項１に記載のガス吸収分光装置。
【請求項５】
前記測定対象ガスがエンジンの燃焼室内のガスであって、
前記圧力関連値が、前記エンジンのクランク角であることを特徴とする請求項１に記載のガス吸収分光装置。
【国際調査報告】
2019159581.tif 000022
2019159581.tif 000023
2019159581.tif 000024
2019159581.tif 000025