(57)【要約】【課題】天然ガスの発熱量値の計測性能を向上させる。【解決手段】カウント値Ｎに応じて定められている閾値ＰＨthをテーブルＴＡから読み出し（Ｓ１０３）、基準ピークハイトＰＨbaseと今回測定されたピークハイトＰＨmeasとの差をΔＰＨとして求め（Ｓ１０４）、ΔＰＨと閾値ＰＨthとを比較する（Ｓ１０５）。ΔＰＨが閾値ＰＨth以下であれば（Ｓ１０５のＹＥＳ）、ＰＨcalc＝（１−α）×ＰＨbase＋α×ＰＨmeasとして、ピークハイトＰＨmeasの補正値を求め（Ｓ１１１）、Ｎ＝Ｎ−１とする（Ｓ１１４）。ΔＰＨが閾値ＰＨthよりも大きければ（Ｓ１０５のＮＯ）、Ｎ＝Ｎ−１とし（Ｓ１１４）、ΔＰＨの算出に用いたピークハイトＰＨmeasを新しい基準ピークハイトＰＨbaseとする（Ｓ１０９）。Ｎ＝０であれば（Ｓ１０６のＹＥＳ）、カウント値Ｎを初期値に再セットし、ΔＰＨの算出に用いたピークハイトＰＨmeasを新しい基準ピークハイトＰＨbaseとする（Ｓ１０９）。【選択図】  図３
2014070949.tif 000002
【特許請求の範囲】
【請求項１】
周期的に、キャリアガスとサンプルガスとを混合し、この混合ガスに含まれる被測定ガス成分を測定周期期間内に予め設定されているゲート時間帯に出現させ、前記測定周期の開始点からピーク値が出現するまでの時間をリテンションタイムとし、このリテンションタイム経過時点における前記被測定ガス成分のクロマトグラムの値と、前記ゲート時間帯の開始点における前記被測定ガス成分のクロマトグラムの値と前記ゲート時間帯の終了点における前記被測定ガス成分のクロマトグラムの値とを結ぶゲート間接線上における前記リテンションタイム経過時点のクロマトグラムの値との差を、前記被測定ガス成分のクロマトグラムのピークハイトＰＨmeasとして測定するガスクロマトグラフにおいて、
測定開始時の前記ピークハイトＰＨmeasを基準ピークハイトＰＨbaseとするとともに予めその値が初期値としてセットされているカウント値Ｎをカウントダウンする第１ステップと、
前記基準ピークハイトＰＨbaseと今回測定された前記ピークハイトＰＨmeasとの差をΔＰＨとして求める第２ステップと、
前記第２ステップで求められた差ΔＰＨと前記カウント値Ｎに応じて定められている閾値ＰＨthとを比較する第３ステップと、
前記第３ステップでの比較結果に基づいて、前記差ΔＰＨが前記閾値ＰＨth以下であった場合、今回測定された前記ピークハイトＰＨmeasの補正値として前記基準ピークハイトＰＨbaseに近づくようなピークハイトＰＨcalcを求めるとともに、前記カウント値Ｎをカウントダウンして、前記第２ステップの処理に戻る第４ステップと、
前記第３ステップでの比較結果に基づいて、前記差ΔＰＨが前記閾値ＰＨthよりも大きかった場合、前記カウント値Ｎが０であるか否かをチェックし、前記カウント値Ｎが０であった場合にはそのカウント値を前記初期値に再セットし、カウント値Ｎが０でなかった場合にはそのカウント値をカウントダウンし、前記差ΔＰＨの算出に用いられた前記ピークハイトＰＨmeasを新しい前記基準ピークハイトＰＨbaseとして、前記第２ステップの処理に戻る第５ステップと
を備えることを特徴とするガスクロマトグラフにおけるピークハイト処理方法。
【請求項２】
請求項１に記載されたガスクロマトグラフにおけるピークハイト処理方法において、
前記第４ステップは、
前記ピークハイトＰＨcalcを、αを所定の補正係数として、ＰＨcalc＝（１−α）×ＰＨbase＋α×ＰＨmeasなる式より求める
ことを特徴とするガスクロマトグラフにおけるピークハイト処理方法。
【請求項３】
請求項１又は２に記載されたガスクロマトグラフにおけるピークハイト処理方法において、
前記閾値ＰＨthは、前記カウント値が減少するほど大きくなるように定められている
ことを特徴とするガスクロマトグラフにおけるピークハイト処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、周期的に、キャリアガスとサンプルガスとを混合し、この混合ガスに含まれる被測定ガス成分を測定周期期間内に予め設定されているゲート時間帯に出現させ、このゲート時間帯に出現する被測定ガス成分のクロマトグラムのピークハイトを求めて、被測定ガス成分の濃度測定を行うガスクロマトグラフにおけるピークハイト処理方法に関するものである。
【背景技術】
【０００２】
従来より、天然ガスの発熱量値などを測定するために、ガスクロマトグラフが用いられている。ガスクロマトグラフを用いて天然ガスの発熱量値を測定する場合、天然ガスの微量サンプルを採取して分析計に取り込み、天然ガス中に含まれる混合物をカラムによって１１成分に分離し、この分離した１１成分の濃度値を測定し、この測定した濃度値より天然ガスの発熱量値を算出する。
【０００３】
このような用途に用いられるガスクロマトグラフにおいては、周期的に、キャリアガス（ヘリウムガス、窒素ガスなど）とサンプルガスとを混合し、この混合ガスに含まれる被測定ガス成分を測定周期期間内に予め設定されているゲート時間帯に出現させ、このゲート時間帯に出現する被測定ガス成分のクロマトグラムのピークハイトを測定し、この測定したピークハイトに基づいて被測定ガス成分の濃度測定を行う。
【０００４】
図１５は被測定ガス成分のピークハイトの測定状況を例示する図である。この例では、測定周期をＴとしてキャリアガスとサンプルガスとを混合し、この混合ガスに含まれる被測定ガス成分をゲート時間帯Ｇに出現させている。この場合、測定周期Ｔの開始点からピーク値が出現するまでの時間をリテンションタイムＲｔとし、このリテンションタイムＲｔ経過時点における被測定ガス成分のクロマトグラムの値ＣＰRtと、ゲート時間帯Ｇの開始点ＧONにおける被測定ガス成分のクロマトグラムの値ＣＰONとゲート時間帯Ｇの終了点ＧOFFにおける被測定ガス成分のクロマトグラムの値ＣＰOFFとを結ぶゲート間接線上におけるリテンションタイムＲｔ経過時点のクロマトグラムの値ＣＰBRtとの差を、被測定ガス成分のクロマトグラムのピークハイトＰＨとして測定し、この測定したピークハイトＰＨ（ＰＨmeas）に基づいて被測定ガス成分の濃度値を算出する（例えば、特許文献１，２参照）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平５−１４９９３６号公報
【特許文献２】特開平８−３５９５９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、図１６に示すように上述したピークハイトＰＨの測定方法によると、ゲート時間帯Ｇの開始点ＧONから終了点ＧOFFまでベースラインＢＬが直線で移動したと仮定してピークハイトＰＨを測定しているので、ベースラインＢＬの変動によってピークハイトＰＨに測定誤差が生じていた。例えば、図１６に点線で示すように、ベースラインＢＬが上方に変動した場合（ベースラインＢＬ＋に変動した場合）、ＰＨ−が真のピークハイトＰＨとして測定されるべきところ、ＰＨmeasがピークハイトＰＨとして測定されてしまう。また、反対に、ベースラインが下方に変動した場合（ベースラインＢＬ−に変動した場合）、ＰＨ＋が真のピークハイトＰＨとして測定されるべきところ、ＰＨmeasがピークハイトＰＨとして測定されてしまう。このように、ベースラインが変動すると測定されるピークハイトＰＨに誤差が生ずる。
【０００７】
ガスクロマトグラフを用いて天然ガスの発熱量値を測定する場合、天然ガスの発熱量値の計測性能は繰り返し性で規定されることが一般的であり、繰り返し性が小さいほど計測性能が良いものとされる。この天然ガスの発熱量値の測定に際し、ベースラインＢＬの変動によって生じるピークハイトＰＨの測定誤差は、発熱量値の繰り返し性を大きくし、発熱量値の計測性能を悪化させる。
【０００８】
本発明は、このような課題を解決するためになされたもので、その目的とするところは、天然ガスの発熱量値の計測性能を向上させることが可能なガスクロマトグラフにおけるピークハイト処理方法を提供することにある。
【課題を解決するための手段】
【０００９】
このような目的を達成するために、本発明に係るガスクロマトグラフにおけるピークハイト処理方法は、周期的に、キャリアガスとサンプルガスとを混合し、この混合ガスに含まれる被測定ガス成分を測定周期期間内に予め設定されているゲート時間帯に出現させ、測定周期の開始点からピーク値が出現するまでの時間をリテンションタイムとし、このリテンションタイム経過時点における被測定ガス成分のクロマトグラムの値と、ゲート時間帯の開始点における被測定ガス成分のクロマトグラムの値とゲート時間帯の終了点における被測定ガス成分のクロマトグラムの値とを結ぶゲート間接線上におけるリテンションタイム経過時点のクロマトグラムの値との差を、被測定ガス成分のクロマトグラムのピークハイトＰＨmeasとして測定するガスクロマトグラフにおいて、 測定開始時のピークハイトＰＨmeasを基準ピークハイトＰＨbaseとするとともに予めその値が初期値としてセットされているカウント値Ｎをカウントダウンする第１ステップと、基準ピークハイトＰＨbaseと今回測定されたピークハイトＰＨmeasとの差をΔＰＨとして求める第２ステップと、第２ステップで求められた差ΔＰＨとカウント値Ｎに応じて定められている閾値ＰＨthとを比較する第３ステップと、第３ステップでの比較結果に基づいて、差ΔＰＨが閾値ＰＨth以下であった場合、今回測定されたピークハイトＰＨmeasの補正値として基準ピークハイトＰＨbaseに近づくようなピークハイトＰＨcalcを求めるとともに、カウント値Ｎをカウントダウンして、第２ステップの処理に戻る第４ステップと、第３ステップでの比較結果に基づいて、差ΔＰＨが閾値ＰＨthよりも大きかった場合、カウント値Ｎが０であるか否かをチェックし、カウント値Ｎが０であった場合にはそのカウント値を初期値に再セットし、カウント値Ｎが０でなかった場合にはそのカウント値をカウントダウンし、差ΔＰＨの算出に用いられたピークハイトＰＨmeasを新しい基準ピークハイトＰＨbaseとして、第２ステップの処理に戻る第５ステップとを備えることを特徴とする。
【００１０】
本発明では、測定開始時のピークハイトＰＨmeasを基準ピークハイトＰＨbaseとするとともに予めその値が初期値としてセットされているカウント値Ｎをカウントダウンする。例えば、カウント値ＮがＮ＝７としてセットされていた場合、カウント値ＮをカウントダウンしてＮ＝６とする。そして、次の測定周期で、この基準ピークハイトＰＨbaseと今回測定されたピークハイトＰＨmeasとの差をΔＰＨとして求め、この差ΔＰＨとカウント値Ｎに応じて定められている閾値ＰＨthとを比較する。
【００１１】
ここで、ΔＰＨが閾値ＰＨth以下であれば、今回測定されたピークハイトＰＨmeasの補正値として、基準ピークハイトＰＨbaseに近づくようなピークハイトＰＨcalcを求める。例えば、ピークハイトＰＨcalcを、αを所定の補正係数として、ＰＨcalc＝（１−α）×ＰＨbase＋α×ＰＨmeasなる式より求め、今回測定されたピークハイトＰＨmeasの補正値とする。この場合、αはプラスの補正係数であっても、マイナスの補正係数であってもよい。
【００１２】
そして、カウント値Ｎをカウントダウンして、次のピークハイトＰＨmeasを測定し、基準ピークハイトＰＨbaseと今回測定されたピークハイトＰＨmeasとの差をΔＰＨとして求め、この差ΔＰＨとカウント値Ｎに応じて定められている閾値ＰＨthとを比較し、閾値ＰＨth以下であれば、上述と同様にして、基準ピークハイトＰＨbaseに近づくようなピークハイトＰＨcalcを求めて、カウント値Ｎをカウントダウンする、という動作を繰り返す。ここで、例えばカウント値Ｎが減少するほど大きくなるように閾値ＰＨthを定めておくものとすれば、カウント値Ｎが減少するにつれて、すなわち測定が安定するにつれて、閾値ＰＨthが大きくなって行く。
【００１３】
このような動作の繰り返し中、基準ピークハイトＰＨbaseと今回測定されたピークハイトＰＨmeasとの差ΔＰＨが閾値ＰＨthよりも大きくなれば、カウント値Ｎが０であるか否かをチェックし、カウント値Ｎが０であった場合にはそのカウント値を初期値（例えば、Ｎ＝７）に再セットし、カウント値Ｎが０でなかった場合にはそのカウント値をカウントダウンし、差ΔＰＨの算出に用いたピークハイトＰＨmeasを新しい基準ピークハイトＰＨbaseとする。そして、次の測定周期で、この新しい基準ピークハイトＰＨbaseと今回測定されたピークハイトＰＨmeasとの差をΔＰＨとして求め、この差ΔＰＨとカウント値Ｎに応じて定められている閾値ＰＨthとを比較し、ΔＰＨが閾値ＰＨth以下であれば、上述と同様にして、基準ピークハイトＰＨbaseに近づくようなピークハイトＰＨcalcを求めて、カウント値Ｎをカウントダウンするという動作を繰り返す。
【発明の効果】
【００１４】
本発明によれば、周期的に、基準ピークハイトＰＨbaseと今回測定されたピークハイトＰＨmeasとの差をΔＰＨとして求め、この差ΔＰＨがカウント値Ｎに応じて定められている閾値ＰＨthよりも大きかった場合、カウント値Ｎが０であるか否かをチェックし、カウント値Ｎが０であった場合にはそのカウント値を初期値に再セットし、カウント値Ｎが０でなかった場合にはそのカウント値をカウントダウンし、差ΔＰＨの算出に用いられたピークハイトＰＨmeasを新しい基準ピークハイトＰＨbaseとする一方、この差ΔＰＨがカウント値Ｎに応じて定められている閾値ＰＨth以下であった場合、今回測定されたピークハイトＰＨmeasの補正値として基準ピークハイトＰＨbaseに近づくようなピークハイトＰＨcalcを求めるとともに、カウント値Ｎをカウントダウンするようにしたので、閾値ＰＨthをカウント値Ｎに応じて自動的に調整しながら、また、ΔＰＨが閾値ＰＨthよりも大きくなれば、ΔＰＨの算出に用いたピークハイトＰＨmeasを新しい基準ピークハイトＰＨbaseとしながら、この基準ピークハイトＰＨbaseに近づくようにピークハイトＰＨmeasが補正されるものとなり、この補正されたピークハイトＰＨmeas（ＰＨcalc）を用いて被測定ガス成分の濃度測定を行うようにすることにより、天然ガスの発熱量値を測定するような場合、発熱量値の繰り返し性を小さくし、発熱量値の計測性能を向上させることが可能となる。
【図面の簡単な説明】
【００１５】
【図１】本発明の実施に用いるガスクロマトグラフの基本構成を示す概略的なブロック図である。
【図２】このガスクロマトグラフにおいてメモリに格納されたカウント値Ｎと閾値ＰＨthとの関係を定めたテーブルＴＡを例示する図である。
【図３】このガスクロマトグラフにおいてＣＰＵが行う特徴的な処理動作を示すフローチャートである。
【図４】このガスクロマトグラフにおける測定開始時のピークハイトＰＨmeasを例示する図である。
【図５】このガスクロマトグラフにおけるカウント値Ｎ＝６の時（測定状態が「サンプル変動（変動期）」の時）の基準ピークハイトＰＨbaseおよび測定されたピークハイトＰＨmeasを例示する図である。
【図６】このガスクロマトグラフにおけるカウント値Ｎ＝５の時（測定状態が「変動→安定（過渡期）」の時）の基準ピークハイトＰＨbaseおよび測定されたピークハイトＰＨmeasを例示する図である。
【図７】このガスクロマトグラフにおけるカウント値Ｎ＝０の時（測定状態が「安定（安定期）」の時）の基準ピークハイトＰＨbaseおよび測定されたピークハイトＰＨmeasを例示する図である。
【図８】測定状態が「安定（安定期）」の時にΔＰＨ＞ＰＨthとなった状態を示す図である。
【図９】測定状態が「安定（安定期）」の時にΔＰＨ＞ＰＨthとなって新しい基準ピークハイトＰＨbaseが定められた状態を示す図である。
【図１０】測定状態が「サンプル変動（変動期）」の時にΔＰＨ＞ＰＨthとなった状態を示す図である。
【図１１】測定状態が「サンプル変動（変動期）」の時にΔＰＨ＞ＰＨthとなって新しい基準ピークハイトＰＨbaseが定められた状態を示す図である。
【図１２】測定状態が「変動→安定（過渡期）」の時にΔＰＨ＞ＰＨthとなった状態を示す図である。
【図１３】測定状態が「変動→安定（過渡期）」の時にΔＰＨ＞ＰＨthとなって新しい基準ピークハイトＰＨbaseが定められた状態を示す図である。
【図１４】実施の形態による効果（発熱量測定結果）を説明する図である。
【図１５】従来のガスクロマトグラフにおける被測定ガス成分のピークハイトの測定状況を例示する図である。
【図１６】ベースラインの変動によって生じるピークハイトの測定誤差を説明する図である。
【発明を実施するための形態】
【００１６】
以下、本発明の実施の形態を図面に基づいて詳細に説明する。図１はこの発明の実施に用いるガスクロマトグラフの基本構成を示す概略的なブロック図である。
【００１７】
図１において、１は恒温槽を形成し所定温度に保持されるアナライザ本体、２はこのアナライザ本体１内に配置されるサンプルバルブ、３はカラム、４は検出器、５は計量管、６はキャリアガス（ヘリウムガス等）ＣＧを所定圧に減圧する減圧弁、７はＣＰＵ、８はメモリである。
【００１８】
このガスクロマトグラフでは、周期的に、サンプルバルブ２の流路を切り替えることにより、計量管５によって分取した測定すべきサンプルガスＳＧをキャリアガスＣＧと混合してカラム３に送り込み、このカラム３内で各ガス成分を分離しながら検出器４へ給送し、この検出器４により各ガス成分の熱伝導度を測定し、この測定した熱伝導度に基づいて各ガス成分の濃度を測定する。この濃度測定はメモリ８に格納されたプログラムに従うＣＰＵ７の処理動作により実行される。
【００１９】
また、メモリ８には、後述するカウント値Ｎと閾値ＰＨthとの関係を定めたテーブルＴＡが格納されている。この例では、図２に示すように、カウント値ＮがＮ＝０である場合、カウント値ＮがＮ＞０である場合、カウント値ＮがＮ＞５である場合について、基準ピークハイトＰＨbaseの大きさを３つに区分し、またゲインを３つに区分し、それぞれについて閾値ＰＨthを定めている。基本的には、カウント値Ｎが減少するほど閾値ＰＨthが大きくなるように、基準ピークハイトＰＨbaseが大きくなるほど閾値ＰＨthが大きくなるように、ゲインが大きくなるほど閾値ＰＨthが大きくなるように、閾値ＰＨthを定めている。
【００２０】
以下、図３に示すフローチャートを参照して、本実施の形態のガスクロマトグラフにおいて、ＣＰＵ７が行う特徴的な処理動作を説明する。
【００２１】
〔測定開始〕
ＣＰＵ７は、カウント値Ｎの初期値を７（Ｎ＝７）として、また基準ピークハイトＰＨbaseの初期値を０（ＰＨbase＝０）として（ステップＳ１０１）、被測定ガス成分のクロマトグラムのピークハイトＰＨmeasの測定を開始し（ステップＳ１０２）、測定開始時のピークハイトＰＨmeasを得る（図４に示す波形Ｉ参照）。
【００２２】
そして、ＣＰＵ７は、その時のカウント値Ｎと基準ピークハイトＰＨbaseとゲインとに応じた閾値ＰＨthをテーブルＴＡから読み出す（ステップＳ１０３）。この時、基準ピークハイトＰＨbaseは０であり、カウント値Ｎは７であるから、ゲインを「０」とした場合、閾値ＰＨthが「８」として読み出される。以下、この例では、ゲインを「０」として説明する。
【００２３】
次に、ＣＰＵ７は、ステップＳ１０２で測定されたピークハイトＰＨmeasと基準ピークハイトＰＨbaseとの差ΔＰＨをΔＰＨ＝｜ＰＨmeas−ＰＨbase｜として求め（ステップＳ１０４）、この求めた差ΔＰＨと閾値ＰＨthとを比較する（ステップＳ１０５）。この場合、ＰＨbase＝０であるから、ΔＰＨはＰＨmeasとして求められ、ＰＨthは「８」であるので、ΔＰＨ＞ＰＨthとなる。
【００２４】
このため、ＣＰＵ７は、ステップＳ１０５でのＮＯに応じてステップＳ１０６へ進み、カウント値Ｎが０でるか否かをチェックする。この場合、Ｎ＝７とされているので、ステップＳ１０７へ進み、カウント値ＮをカウントダウンしてＮ＝６とする。そして、ステップＳ１０２で得た測定開始時のピークハイトＰＨmeasを基準ピークハイトＰＨbaseとするとともに（ステップＳ１０９）、この測定開始時のピークハイトＰＨmeasをＰＨoutとして出力し（ステップＳ１１０，Ｓ１１５）、ステップＳ１０２に戻る。
【００２５】
〔測定状態：サンプル変動（変動期）〕
ＣＰＵ７は、ステップＳ１０２に戻って、次の測定周期のピークハイトＰＨmeasを測定すると（図５に示す波形II参照）、その時のカウント値Ｎと基準ピークハイトＰＨbaseとゲインとに応じた閾値ＰＨthをテーブルＴＡから読み出す（ステップＳ１０３）。この時、例えば、基準ピークハイトＰＨbaseをＰＨbase＞８ｋとすると、カウント値Ｎは６、ゲインは「０」であるから、閾値ＰＨthが「１５」として読み出される。
【００２６】
次に、ＣＰＵ７は、今回測定されたピークハイトＰＨmeasと基準ピークハイトＰＨbaseとの差ΔＰＨをΔＰＨ＝｜ＰＨmeas−ＰＨbase｜として求め（ステップＳ１０４）、この求めた差ΔＰＨと閾値ＰＨthとを比較する（ステップＳ１０５）。
【００２７】
〔ΔＰＨ≦ＰＨthの場合〕
ここで、ΔＰＨがＰＨth以下であれば（ステップＳ１０５のＹＥＳ）、ＣＰＵ７は、αを所定の補正係数とし、下記の（１）式に従ってピークハイトＰＨcalc（ＳＷ処理されたピークハイト）を求め、今回測定されたピークハイトＰＨmeasの補正値とする（ステップＳ１１１）。
ＰＨcalc＝（１−α）×ＰＨbase＋α×ＰＨmeas  ・・・・（１）
【００２８】
この（１）式は、ＰＨcalc＝ＰＨbase＋α（ＰＨmeas−ＰＨbase）を変形して得られた式であり、本実施の形態では補正係数αを例えば＋０．３とする。このような式を用いることによって、今回測定されたピークハイトＰＨmeasの補正値として、基準ピークハイトＰＨbaseに近づくようなピークハイトＰＨcalcが求められる。
【００２９】
例えば、ＰＨbase＝１００、ＰＨmeas＝９０とした場合、ＰＨcalc＝（１−０．３）×１００＋０．３×９０＝７０＋２７＝９７とされる。また、ＰＨbase＝１００、ＰＨmeas＝９５とした場合、ＰＨcalc＝（１−０．３）×１００＋０．３×９５＝７０＋２８．５＝９８．５とされる。
【００３０】
なお、αはプラスの補正係数であっても、マイナスの補正係数であってもよい。例えば、αを−０．３とした場合、上述したＰＨbase＝１００、ＰＨmeas＝９０の例では、ＰＨcalc＝（１＋０．３）×１００−０．３×９０＝１３０−２７＝１０３とされる。また、上述したＰＨbase＝１００、ＰＨmeas＝９５の例では、ＰＨcalc＝（１＋０．３）×１００−０．３×９５＝１３０−２８．５＝１０１．５とされる。
【００３１】
このようにして、ＣＰＵ７は、基準ピークハイトＰＨbaseに近づくようなピークハイトＰＨcalcを求めた後、ＰＨcalcをＰＨoutとし（ステップＳ１１２）、カウント値ＮがＮ＝０であるか否かをチェックする（ステップＳ１１３）。この場合、カウント値ＮはＮ＝６であるので、ステップＳ１１３でのＮＯに応じてステップＳ１１４へ進み、カウント値ＮをカウントダウンしてＮ＝５とし、ＰＨoutを出力して（ステップＳ１１５）、ステップＳ１０２に戻る。
【００３２】
〔測定状態：変動→安定（過渡期）〕
ＣＰＵ７は、ステップＳ１０２に戻って、次の測定周期のピークハイトＰＨmeasを測定すると（図６に示す波形III参照）、その時のカウント値Ｎと基準ピークハイトＰＨbaseとゲインとに応じた閾値ＰＨthをテーブルＴＡから読み出す（ステップＳ１０３）。この時、例えば、基準ピークハイトＰＨbaseをＰＨbase＞８ｋとすると、カウント値Ｎは５、ゲインは「０」であるから、閾値ＰＨthが「３０」として読み出される。
【００３３】
次に、ＣＰＵ７は、今回測定されたピークハイトＰＨmeasと基準ピークハイトＰＨbaseとの差ΔＰＨをΔＰＨ＝｜ＰＨmeas−ＰＨbase｜として求め（ステップＳ１０４）、この求めた差ΔＰＨと閾値ＰＨthとを比較する（ステップＳ１０５）。
【００３４】
〔ΔＰＨ≦ＰＨthの場合〕
ここで、ΔＰＨがＰＨth以下であれば（ステップＳ１０５のＹＥＳ）、ＣＰＵ７は、ステップＳ１１１へ進み、前記（１）式に従って、基準ピークハイトＰＨbaseに近づくようなピークハイトＰＨcalcを求め、今回測定されたピークハイトＰＨmeasの補正値とする。
【００３５】
そして、この求めたピークハイトＰＨcalcをＰＨoutとし（ステップＳ１１２）、カウント値ＮがＮ＝０であるか否かをチェックする（ステップＳ１１３）。この場合、カウント値ＮはＮ＝５であるので、ステップＳ１１３でのＮＯに応じてステップＳ１１４へ進み、カウント値ＮをカウントダウンしてＮ＝４とし、ＰＨoutを出力して（ステップＳ１１５）、ステップＳ１０２に戻る。
【００３６】
ＣＰＵ７は、ΔＰＨ≦ＰＨthである場合、カウント値ＮがＮ＝０となるまで、閾値ＰＨthを「３０」として、ステップＳ１０２〜Ｓ１０５、Ｓ１１１〜Ｓ１１５の処理動作を繰り返す。
【００３７】
〔測定状態：安定（安定期）〕
ＣＰＵ７は、カウント値ＮがＮ＝０となると（ステップＳ１１３のＹＥＳ）、ステップＳ１１４を経ずに直ちにステップＳ１１５へ進む。すなわち、カウント値Ｎのカウントダウンを行わず、Ｎ＝０としたままステップＳ１１５へ進んで、ステップＳ１０２に戻る。
【００３８】
ＣＰＵ７は、ステップＳ１０２に戻って、次の測定周期のピークハイトＰＨmeasを測定すると（図７に示す波形VIII参照）、その時のカウント値Ｎと基準ピークハイトＰＨbaseとゲインとに応じた閾値ＰＨthをテーブルＴＡから読み出す（ステップＳ１０３）。この時、例えば、基準ピークハイトＰＨbaseをＰＨbase＞８ｋとすると、カウント値Ｎは０、ゲインは「０」であるから、閾値ＰＨthが「６０」として読み出される。
【００３９】
次に、ＣＰＵ７は、ステップＳ１０２で測定されたピークハイトＰＨmeasと基準ピークハイトＰＨbaseとの差ΔＰＨをΔＰＨ＝｜ＰＨmeas−ＰＨbase｜として求め（ステップＳ１０４）、この求めた差ΔＰＨと閾値ＰＨthとを比較する（ステップＳ１０５）。
【００４０】
〔ΔＰＨ≦ＰＨthの場合〕
ここで、ΔＰＨがＰＨth以下であれば（ステップＳ１０５のＹＥＳ）、ＣＰＵ７は、ステップＳ１１１へ進み、前記（１）式に従って、基準ピークハイトＰＨbaseに近づくようなピークハイトＰＨcalcを求め、今回測定されたピークハイトＰＨmeasの補正値とする。
【００４１】
そして、この求めたピークハイトＰＨcalcをＰＨoutとし（ステップＳ１１２）、カウント値ＮがＮ＝０であるか否かをチェックする（ステップＳ１１３）。この場合、カウント値ＮはＮ＝０であるので、ステップＳ１１３でのＹＥＳに応じて直ちにステップＳ１１５へ進み、測定されたピークハイトＰＨoutを出力して、ステップＳ１０２に戻る。
【００４２】
ＣＰＵ７は、ΔＰＨ≦ＰＨthである場合、閾値ＰＨthを「６０」として、ステップＳ１０２〜Ｓ１０５、Ｓ１１１〜Ｓ１１３，Ｓ１１５の処理動作を繰り返す。
【００４３】
〔測定状態が「安定（安定期）」の時にΔＰＨ＞ＰＨthとなった場合〕
ＣＰＵ７は、測定状態が「安定（安定期）」の時にΔＰＨ＞ＰＨthとなると（ステップＳ１０５のＮＯ、図８に示す波形VX参照）、カウント値ＮがＮ＝０か否かをチェックする（ステップＳ１０６）。この場合、カウント値ＮはＮ＝０であるので（ステップＳ１０６のＹＥＳ）、初期値７に再セットし（ステップＳ１０８）、すなわちカウント値ＮをリセットしてＮ＝７とし、差ΔＰＨの算出に用いたピークハイトＰＨmeasを新しい基準ピークハイトＰＨbaseとする（ステップＳ１０９、図９に示す波形VX参照）。すなわち、それまでの基準ピークハイトＰＨbaseを破棄し、新しい基準ピークハイトＰＨbaseに変更する（ステップＳ１１０）。そして、ＣＰＵ７は、ＰＨmeasをＰＨoutとして出力し（ステップＳ１１５）、ステップＳ１０２に戻る。
【００４４】
ＣＰＵ７は、ステップＳ１０２に戻って、次の測定周期のピークハイトＰＨmeasを測定すると（図９に示す波形VX+1参照）、その時のカウント値Ｎと基準ピークハイトＰＨbaseとゲインとに応じた閾値ＰＨthをテーブルＴＡから読み出す（ステップＳ１０３）。この時、例えば、基準ピークハイトＰＨbaseをＰＨbase＞８ｋとすると、カウント値Ｎは７、ゲインは「０」であるから、閾値ＰＨthが「１５」として読み出される。すなわち、測定状態が「サンプル変動（変動期）」に戻り、それまで「６０」であった閾値ＰＨthが「１５」に戻される。
【００４５】
次に、ＣＰＵ７は、ステップＳ１０２で測定されたピークハイトＰＨmeasと基準ピークハイトＰＨbaseとの差ΔＰＨをΔＰＨ＝｜ＰＨmeas−ＰＨbase｜として求め（ステップＳ１０４）、この求めた差ΔＰＨと閾値ＰＨthとを比較する（ステップＳ１０５）。
【００４６】
ここで、ΔＰＨがＰＨth以下であれば（ステップＳ１０５のＹＥＳ）、ＣＰＵ７は、ステップＳ１１１へ進み、前記（１）式に従って、基準ピークハイトＰＨbaseに近づくようなピークハイトＰＨcalcを求め、今回測定されたピークハイトＰＨmeasの補正値とする。
【００４７】
そして、この求めたピークハイトＰＨcalcをＰＨoutとし（ステップＳ１１２）、カウント値ＮがＮ＝０であるか否かをチェックする（ステップＳ１１３）。この場合、カウント値ＮはＮ＝７であるので、ステップＳ１１３でのＮＯに応じてステップＳ１１４へ進み、カウント値ＮをカウントダウンしてＮ＝６とし、測定されたピークハイトＰＨoutを出力して（ステップＳ１１５）、ステップＳ１０２に戻る。
【００４８】
〔測定状態が「サンプル変動（変動期）」の時にΔＰＨ＞ＰＨthとなった場合〕
ＣＰＵ７は、測定状態が「サンプル変動（変動期）」の時にΔＰＨ＞ＰＨthとなると（ステップＳ１０５のＮＯ、図１０に示す波形II参照）、カウント値Ｎをチェックする（ステップＳ１０６）。
【００４９】
図１０に示した例の場合、カウント値ＮはＮ＝６であるので、ステップＳ１０６のＮＯに応じてステップＳ１０７へ進み、カウント値ＮをカウントダウンしてＮ＝５とし、今回測定されたピークハイトＰＨmeasを基準ピークハイトＰＨbaseとするとともに（ステップＳ１０９、図１１に示す波形II参照）、今回測定されたピークハイトＰＨmeasをＰＨoutとし（ステップＳ１１０）、ＰＨoutを出力し（ステップＳ１１５）、ステップＳ１０２に戻る。これにより、新しい基準ピークハイトＰＨbaseが定められるとともに、カウント値Ｎのカウントダウンが続けられる。
【００５０】
〔測定状態が「変動→安定（過渡期）」の時にΔＰＨ＞ＰＨthとなった場合〕
ＣＰＵ７は、測定状態が「変動→安定（過渡期）」の時にΔＰＨ＞ＰＨthとなると（ステップＳ１０５のＮＯ、図１２に示す波形III参照）、カウント値Ｎをチェックする（ステップＳ１０６）。
【００５１】
図１２に示した例の場合、カウント値ＮはＮ＝５であるので、ステップＳ１０６のＮＯに応じてステップＳ１０７へ進み、カウント値ＮをカウントダウンしてＮ＝４とし、今回測定されたピークハイトＰＨmeasを基準ピークハイトＰＨbaseとするとともに（ステップＳ１０９、図１３に示す波形III参照）、今回測定されたピークハイトＰＨmeasをＰＨoutとし（ステップＳ１１０）、ＰＨoutを出力し（ステップＳ１１５）、ステップＳ１０２に戻る。これにより、新しい基準ピークハイトＰＨbaseが定められるとともに、カウント値Ｎのカウントダウンが続けられる。
【００５２】
このようにして、本実施の形態では、閾値ＰＨthをカウント値Ｎに応じて自動的に調整しながら、また、ΔＰＨが閾値ＰＨthよりも大きくなれば、ΔＰＨの算出に用いたピークハイトＰＨmeasを新しい基準ピークハイトＰＨbaseとしながら、この基準ピークハイトＰＨbaseに近づくようにピークハイトＰＨmeasが補正されるものとなり、この補正されたピークハイトＰＨmeas（ＰＨcalc）を用いて被測定ガス成分の濃度測定を行うようにすることにより、天然ガスの発熱量値を測定するような場合、発熱量値の繰り返し性を小さくし、発熱量値の計測性能を向上させることが可能となる。
【００５３】
なお、上述した実施の形態では、ピークハイトＰＨcalcを求める際に用いる補正係数αを０．３としたが、０．３に限られるものでないことは言うまでもない。本実施の形態では、０．３〜０．６の間で０．１刻みに変動させて実験を繰り返し、目標とする発熱量値の繰り返し性をクリアできる値として０．３を採用している。
【００５４】
図１４に工場の出荷前検査で使用する濃度の異なる３本のガス（medium，low，high）を用いて発熱量値の測定を実施した結果を示す。同図（ａ）は「high」、同図（ｂ）は「medium」、同図（ｃ）は「low」による発熱量値の測定結果であり、測定結果Ｉ（改良前）は本発明を適用しなかったもの、測定結果II（改良後）は本発明を適用したものである。同図（ｄ）はデータの分析結果であり、「ｒｅｐ％」と表示されている箇所が繰り返し性を示す。
【００５５】
図１４に示した結果からも分かるように、「high」の例では、改良前では繰り返し性が０．０３６％であったものが、改良後では０．０１１％となり、良くなっている。「medium」の例でも、改良前では繰り返し性が０．０２７％であったものが、改良後では０．００８％となり、良くなっている。「low」の例も、改良前では繰り返し性が０．０２３％であったものが、改良後では０．００８％となり、良くなっている。
【００５６】
なお、上述した実施の形態では、ピークハイトＰＨcalcを求める際に用いる補正係数αを０．３としたが、０．３に限られるものでないことは言うまでもない。本実施の形態では、０．３〜０．６の間で０．１刻みに変動させて実験を繰り返し、目標とする発熱量値の繰り返し性を満足できる値として０．３を採用している。
【００５７】
また、上述した実施の形態では、カウント値Ｎの値に応じて測定状態を「サンプル変動（変動期）」、「変動→安定（過渡期）」、「安定（安定期）」の３段階に分けたが、さらに多くの段階に分けてもよく、また２段階で十分な場合もあるし、閾値ＰＨthも製品に応じて自由に変更しても構わない。
【００５８】
〔実施の形態の拡張〕
以上、実施の形態を参照して本発明を説明したが、本発明は上記の実施の形態に限定されるものではない。本発明の構成や詳細には、本発明の技術思想の範囲内で当業者が理解し得る様々な変更をすることができる。
【符号の説明】
【００５９】
１…アナライザ本体、２…サンプルバルブ、３…カラム、４…検出器、５…計量管、６…減圧弁、７…ＣＰＵ、８…メモリ。
【図１】
2014070949.tif 000003
【図２】
2014070949.tif 000004
【図３】
2014070949.tif 000005
【図４】
2014070949.tif 000006
【図５】
2014070949.tif 000007
【図６】
2014070949.tif 000008
【図７】
2014070949.tif 000009
【図８】
2014070949.tif 000010
【図９】
2014070949.tif 000011
【図１０】
2014070949.tif 000012
【図１１】
2014070949.tif 000013
【図１２】
2014070949.tif 000014
【図１３】
2014070949.tif 000015
【図１４】
2014070949.tif 000016
【図１５】
2014070949.tif 000017
【図１６】
2014070949.tif 000018