説明

発熱量算出式作成システム、発熱量算出式の作成方法、発熱量測定システム、及び発熱量の測定方法

【課題】ガスの発熱量を容易に測定可能な発熱量測定システムを提供する。
【解決手段】ガスが注入されるチャンバ101と、チャンバ101に配置された発熱素子及びガスの温度を計測する測温素子を含むマイクロチップ8と、発熱素子の温度及び測温素子の温度が等しい場合の測温素子の抵抗値に基づいて、発熱素子の抵抗値のドリフトを補正する補正部301と、ガスの計測された温度に基づいて、発熱素子を複数の発熱温度で発熱させたときのガスの放熱係数の計測値を計測する計測機構10と、複数の発熱温度に対する放熱係数を独立変数とし、発熱量を従属変数とする発熱量算出式を保存する式記憶装置402と、発熱量算出式の放熱係数の独立変数にガスの放熱係数の計測値を代入し、ガスの発熱量の値を算出する発熱量算出モジュール305と、を備える、発熱量測定システムを提供する。

【発明の詳細な説明】
【技術分野】
【0001】
本発明はガス検査技術に係り、発熱量算出式作成システム、発熱量算出式の作成方法、発熱量測定システム、及び発熱量の測定方法に関する。
【背景技術】
【0002】
従来、混合ガスの発熱量を求める際には、高価なガスクロマトグラフィ装置等を用いて混合ガスの成分を分析する必要がある。また、混合ガスの熱伝導率及び混合ガスにおける音速を測定することにより、混合ガスに含まれるメタン(CH4)、プロパン(C38)、窒素(N2)、及び炭酸ガス(CO2)の成分比率を算出し、混合ガスの発熱量を求める方法も提案されている(例えば、特許文献1参照。)。
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特表2004−514138号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
しかし、特許文献1に開示された方法は、熱伝導率を測定するためのセンサの他に、音速を測定するための高価な音速センサが必要である。そこで、本発明は、ガスの発熱量を容易に測定可能な発熱量算出式作成システム、発熱量算出式の作成方法、発熱量測定システム、及び発熱量の測定方法を提供することを目的の一つとする。
【課題を解決するための手段】
【0005】
本発明の態様によれば、(a)複数の混合ガスが注入される容器と、(b)容器に配置された発熱素子と、(c)容器に配置され、複数の混合ガスのそれぞれの温度を計測する測温素子と、(d)発熱素子の温度及び測温素子の温度が等しい場合の測温素子の抵抗値に基づいて、発熱素子の抵抗値のドリフトを補正する補正部と、(e)複数の混合ガスのそれぞれの計測された温度に基づいて、発熱素子を複数の発熱温度で発熱させたときの複数の混合ガスのそれぞれの放熱係数の値を計測する計測機構と、(f)複数の混合ガスの既知の発熱量の値と、複数の発熱温度に対する放熱係数の計測値と、に基づいて、複数の発熱温度に対する放熱係数を独立変数とし、発熱量を従属変数とする発熱量算出式を作成する式作成モジュールと、を備える、発熱量算出式作成システムが提供される。
【0006】
本発明の態様によれば、(a)複数の混合ガスを準備することと、(b)複数の混合ガスのそれぞれの温度を測温素子で計測することと、(c)発熱素子の温度及び測温素子の温度が等しい場合の測温素子の抵抗値に基づいて、発熱素子の抵抗値のドリフトを補正することと、(d)複数の混合ガスのそれぞれの計測された温度に基づいて、発熱素子を複数の発熱温度で発熱させたときの複数の混合ガスのそれぞれの放熱係数の値を計測することと、(e)複数の混合ガスの既知の発熱量の値と、複数の発熱温度に対する放熱係数の計測値と、に基づいて、複数の発熱温度に対する放熱係数を独立変数とし、発熱量を従属変数とする発熱量算出式を作成することと、を含む、発熱量算出式の作成方法が提供される。
【0007】
本発明の態様によれば、(a)発熱量が未知の計測対象混合ガスが注入される容器と、(b)容器に配置された発熱素子と、(c)容器に配置され、計測対象混合ガスの温度を計測する測温素子と、(e)発熱素子の温度及び測温素子の温度が等しい場合の測温素子の抵抗値に基づいて、発熱素子の抵抗値のドリフトを補正する補正部と、(f)計測対象混合ガスの計測された温度に基づいて、発熱素子を複数の発熱温度で発熱させたときの計測対象混合ガスの放熱係数の値を計測する計測機構と、(g)複数の発熱温度に対する放熱係数を独立変数とし、発熱量を従属変数とする発熱量算出式を保存する式記憶装置と、(h)発熱量算出式の放熱係数の独立変数に計測対象混合ガスの放熱係数の計測値を代入し、計測対象混合ガスの発熱量の値を算出する発熱量算出モジュールと、を備える、発熱量測定システムが提供される。
【0008】
本発明の態様によれば、(a)発熱量が未知の計測対象混合ガスを用意することと、(b)計測対象混合ガスの温度を測温素子で計測することと、(c)発熱素子の温度及び測温素子の温度が等しい場合の測温素子の抵抗値に基づいて、発熱素子の抵抗値のドリフトを補正することと、(d)計測対象混合ガスの計測された温度に基づいて、発熱素子を複数の発熱温度で発熱させたときの計測対象混合ガスの放熱係数の値を計測することと、(e)複数の発熱温度に対する放熱係数を独立変数とし、発熱量を従属変数とする発熱量算出式を用意することと、(f)発熱量算出式の放熱係数の独立変数に計測対象混合ガスの放熱係数の計測値を代入し、計測対象混合ガスの発熱量の値を算出することと、を含む、発熱量の測定方法が提供される。
【0009】
本発明の態様によれば、(a)ガスが注入される容器と、(b)容器に配置された発熱素子と、(c)容器に配置され、ガスの温度を計測する測温素子と、(d)発熱素子の温度及び測温素子の温度が等しい場合の測温素子の抵抗値に基づいて、発熱素子の抵抗値のドリフトを補正する補正部と、(e)ガスの計測された温度に基づいて、発熱素子を発熱させたときの計測対象混合ガスの放熱係数の値を計測する計測機構と、を備える、物性計測システムが提供される。
【0010】
本発明の態様によれば、(a)ガスを準備することと、(b)ガスの温度を測温素子で計測することと、(c)発熱素子の温度及び測温素子の温度が等しい場合の測温素子の抵抗値に基づいて、発熱素子の抵抗値のドリフトを補正することと、(d)ガスの計測された温度に基づいて、発熱素子を発熱させたときのガスの放熱係数の値を計測することと、を含む、物性計測方法が提供される。
【発明の効果】
【0011】
本発明によれば、ガスの発熱量を容易に測定可能な発熱量算出式作成システム、発熱量算出式の作成方法、発熱量測定システム、及び発熱量の測定方法を提供可能である。
【図面の簡単な説明】
【0012】
【図1】本発明の第1の実施の形態に係るマイクロチップの斜視図である。
【図2】本発明の第1の実施の形態に係るマイクロチップの図1のII−II方向から見た断面図である。
【図3】本発明の第1の実施の形態に係る発熱素子に関する回路図である。
【図4】本発明の第1の実施の形態に係る測温素子に関する回路図である。
【図5】本発明の第1の実施の形態に係る発熱素子の温度と、ガスの放熱係数の関係を示すグラフである。
【図6】本発明の第1の実施の形態に係るガス物性値測定システムの第1の模式図である。
【図7】本発明の第1の実施の形態に係るガス物性値測定システムの第2の模式図である。
【図8】本発明の第1の実施の形態に係る発熱量算出式の作成方法を示すフローチャートである。
【図9】本発明の第2の実施の形態に係るガス物性値測定システムを示す模式図である。
【図10】本発明の第2の実施の形態に係る熱伝導率と放熱係数との関係を示すグラフである。
【図11】本発明の第3の実施の形態に係るガス物性値測定システムを示す模式図である。
【図12】本発明の第3の実施の形態に係るガスの濃度と放熱係数との関係を示すグラフである。
【図13】本発明の第4の実施の形態に係るガス物性値測定システムを示す模式図である。
【図14】本発明の第4の実施の形態に係る発熱量の測定方法を示すフローチャートである。
【図15】本発明の実施の形態の実施例1に係るサンプル混合ガスの組成と、発熱量と、を示す表である。
【図16】本発明の実施の形態の実施例1に係るサンプル混合ガスの算出された発熱量と、真の発熱量と、を示すグラフである。
【図17】本発明の実施の形態の実施例1に係るサンプル混合ガスの真の発熱量と、算出された発熱量と、の関係を示すグラフである。
【発明を実施するための形態】
【0013】
以下に本発明の実施の形態を説明する。以下の図面の記載において、同一又は類似の部分には同一又は類似の符号で表している。但し、図面は模式的なものである。したがって、具体的な寸法等は以下の説明を照らし合わせて判断するべきものである。また、図面相互間においても互いの寸法の関係や比率が異なる部分が含まれていることは勿論である。
【0014】
(第1の実施の形態)
まず、斜視図である図1、及びII−II方向から見た断面図である図2を参照して、第1の実施の形態に係るガス物性値測定システムに用いられるマイクロチップ8について説明する。マイクロチップ8は、キャビティ66が設けられた基板60、及び基板60上にキャビティ66を覆うように配置された絶縁膜65を備える。基板60の厚みは、例えば0.5mmである。また、基板60の縦横の寸法は、例えばそれぞれ1.5mm程度である。絶縁膜65のキャビティ66を覆う部分は、断熱性のダイアフラムをなしている。さらにマイクロチップ8は、絶縁膜65のダイアフラムの部分に設けられた発熱素子61と、発熱素子61を挟むように絶縁膜65のダイアフラムの部分に設けられた第1の測温素子62及び第2の測温素子63と、基板60上に設けられた保温素子64と、を備える。
【0015】
発熱素子61は、キャビティ66を覆う絶縁膜65のダイアフラムの部分の中心に配置されている。発熱素子61は、例えば抵抗器であり、電力を与えられて発熱し、発熱素子61に接する雰囲気ガスを加熱する。第1の測温素子62及び第2の測温素子63は、例えば抵抗器であり、発熱素子61が発熱する前の雰囲気ガスのガス温度を検出する。なお、第1の測温素子62及び第2の測温素子63のいずれかのみを用いてガス温度を検出してもよい。あるいは、第1の測温素子62が検出したガス温度と、第2の測温素子63が検出したガス温度と、の平均値を、ガス温度として採用してもよい。以下においては、第1の測温素子62が検出したガス温度を採用する例を説明するが、これに限定されない。
【0016】
保温素子64は、例えば抵抗器であり、電力を与えられて発熱し、基板60の温度を一定に保つ。基板60の材料としては、シリコン(Si)等が使用可能である。絶縁膜65の材料としては、酸化ケイ素(SiO2)等が使用可能である。キャビティ66は、異方性エッチング等により形成される。また発熱素子61、第1の測温素子62、第2の測温素子63、及び保温素子64のそれぞれの材料には白金(Pt)等が使用可能であり、リソグラフィ法等により形成可能である。また、発熱素子61、第1の測温素子62、及び第2の測温素子63は、同一の部材からなっていてもよい。
【0017】
マイクロチップ8は、マイクロチップ8の底面に配置された断熱部材18を介して、雰囲気ガスが充填されるチャンバ等の容器に固定される。断熱部材18を介してマイクロチップ8をチャンバ等に固定することにより、マイクロチップ8の温度が、チャンバ等の内壁の温度変動の影響を受けにくくなる。ガラス等からなる断熱部材18の熱伝導率は、例えば1.0W/(m・K)以下である。
【0018】
図3に示すように、発熱素子61の一端には、例えば、オペアンプ170のマイナス(−)入力端子が電気的に接続され、他端は接地される。また、オペアンプ170の−入力端子及び出力端子と並列に、抵抗素子160が接続される。オペアンプ170のプラス(+)入力端子は、直列に接続された抵抗素子162と抵抗素子163との間、直列に接続された抵抗素子163と抵抗素子164との間、直列に接続された抵抗素子164と抵抗素子165との間、又は直列に接続された抵抗素子165と抵抗素子166との間に電気的に接続される。各抵抗素子162−166の抵抗値を適当に定めることにより、例えば5.0Vの電圧Vinを抵抗素子162の一端に印加すると、抵抗素子163と抵抗素子162との間には、例えば3.4Vの電圧VL3が生じる。また、抵抗素子164と抵抗素子163との間には例えば2.4Vの電圧VL2が生じ、抵抗素子165と抵抗素子164との間には例えば1.5Vの電圧VL1が生じる。抵抗素子166と抵抗素子165との間には、例えば0.2Vの電圧VL0が生じる。
【0019】
抵抗素子162及び抵抗素子163の間と、オペアンプの−入力端子との間には、スイッチSW1が設けられており、抵抗素子163及び抵抗素子164の間と、オペアンプの−入力端子との間には、スイッチSW2が設けられている。また、抵抗素子164及び抵抗素子165の間と、オペアンプの−入力端子との間には、スイッチSW3が設けられており、抵抗素子165の接地端子と、オペアンプの−入力端子との間には、スイッチSW4が設けられている。
【0020】
オペアンプ170の+入力端子に3.4Vの電圧VL3を印加する場合、スイッチSW1のみが通電され、スイッチSW2,SW3,SW4は切断される。オペアンプ170の+入力端子に2.4Vの電圧VL2を印加する場合、スイッチSW2のみが通電され、スイッチSW1,SW3,SW4は切断される。オペアンプ170の+入力端子に1.5Vの電圧VL1を印加する場合、スイッチSW3のみが通電され、スイッチSW1,SW2,SW4は切断される。オペアンプ170の+入力端子に0.2Vの電圧VL0を印加する場合、スイッチSW4のみが通電され、スイッチSW1,SW2,SW3は切断される。
【0021】
したがって、スイッチSW1,SW2,SW3,SW4の開閉によって、オペアンプ170の−入力端子に4段階の電圧のいずれかを印加可能である。そのため、スイッチSW1,SW2,SW3,SW4の開閉によって、発熱素子61の温度を定める印加電圧を4段階に設定可能である。ここで、オペアンプ170の+入力端子に0.2Vの電圧VL0を印加した場合の発熱素子61の温度をTH0、オペアンプ170の+入力端子に1.5Vの電圧VL1を印加した場合の発熱素子61の温度をTH1とする。また、オペアンプ170の+入力端子に2.4Vの電圧VL2を印加した場合の発熱素子61の温度をTH2、オペアンプ170の+入力端子に3.4Vの電圧VL3を印加した場合の発熱素子61の温度をTH3とする。
【0022】
オペアンプ170の+入力端子に0.2Vのような弱い電圧VL0を印加した場合、発熱素子61の温度TH0は、オペアンプ170の+入力端子に電圧を印加しなかった場合と比較してほとんど上昇しない。そのため、発熱素子61の温度TH0は、雰囲気ガスの温度と同じになるか、又は近似する。これに対し、オペアンプ170の+入力端子に電圧VL1,VL2,VL3,を印加した場合、発熱素子61の温度TH1,TH2,TH3,は、オペアンプ170の+入力端子に電圧を印加しなかった場合と比較して上昇し、雰囲気ガスの温度よりも高くなる。
【0023】
図1及び図2に示す第1の測温素子62は、図4に示すように、例えば抵抗ブリッジ回路の一部をなしている。抵抗ブリッジ回路は、第1の測温素子62と直列に接続された抵抗素子181と、第1の測温素子62及び抵抗素子181と並列に接続された抵抗素子182,183を備える。抵抗ブリッジ回路には、オペアンプ171が接続されている。図1及び図2に示す第2の測温素子63も、例えば抵抗ブリッジ回路の一部をなしている。
【0024】
発熱素子61の抵抗値は、発熱素子61の温度によって変化する。発熱素子61の温度THと、発熱素子61の抵抗値RHの関係は、下記(1)式で与えられる。
RH = RH_STD×[1+α(TH-TH_STD) + β(TH-TH_STD)2] ・・・(1)
ここで、TH_STDは発熱素子61の標準温度を表し、例えば20℃である。RH_STDは標準温度TH_STDにおける予め計測された発熱素子61の抵抗値を表す。αは1次の抵抗温度係数を表す。βは2次の抵抗温度係数を表す。
【0025】
発熱素子61の抵抗値RHは、発熱素子61の駆動電力PHと、発熱素子61の通電電流IHから、下記(2)式で与えられる。
RH = PH / IH2 ・・・(2)
あるいは発熱素子61の抵抗値RHは、発熱素子61にかかる電圧VHと、発熱素子61の通電電流IHから、下記(3)式で与えられる。
RH = VH / IH ・・・(3)
【0026】
ここで、発熱素子61の温度THは、発熱素子61と雰囲気ガスの間が熱的に平衡になったときに安定する。なお、熱的に平衡な状態とは、発熱素子61の発熱と、発熱素子61から雰囲気ガスへの放熱とが釣り合っている状態をいう。平衡状態において、下記(4)式に示すように、発熱素子61の駆動電力PHを、発熱素子61の温度THと雰囲気ガスの温度TIとの差ΔTHで割ることにより、雰囲気ガスの放熱係数MIが得られる。なお、放熱係数MIの単位は、例えばW/℃である。
MI = PH / (TH - TI)
= PH /ΔTH ・・・(4)
【0027】
上記(1)式より、発熱素子61の温度THは下記(5)式で与えられる。
TH = (1 / 2β)×[-α+ [α2 - 4β(1 - RH / RH_STD)]1/2] + TH_STD ・・・(5)
したがって、発熱素子61の温度THと雰囲気ガスの温度TIとの差ΔTHは、下記(6)式で与えられる。
ΔTH = (1 / 2β)×[-α+ [α2 - 4β(1 - RH / RH_STD)]1/2] + TH_STD - TI ・・・(6)
発熱素子61の通電電流IHと、駆動電力PH又は電圧VHは計測可能であるため、上記(2)式又は(3)式から発熱素子61の抵抗値RHを算出可能である。また、雰囲気ガスの温度TIは、図1に示す第1の測温素子62で測定可能である。したがって、図1及び図2に示すマイクロチップ8を用いて、下記(7)式から雰囲気ガスの放熱係数MIが算出可能である。
MI = PH /ΔTH
= PH / [(1 / 2β)×[-α+ [α2 - 4β(1 - RH / RH_STD)]1/2] + TH_STD - TI] ・・・(7)
【0028】
ただし、発熱素子61は高温で発熱するよう高電圧をかけられるため、エレクトロマイグレーションが生じ得る。そのため、発熱素子61の温度が標準温度TH_STDの場合の発熱素子61の抵抗値RH_STDは、経年変化等により予め測定した値からドリフトする場合がある。よって、予め測定した抵抗値RH_STDを用いて(7)式から放熱係数MIを算出すると、誤差が生じる場合がある。
【0029】
エレクトロマイグレーションが生じた後、発熱素子61の温度が標準温度TH_STDの場合の発熱素子61の抵抗値をRH_STD_Dとして、発熱素子61の抵抗値RH_Dは下記(8)式で与えられる。
RH_D = RH_STD_D×[1+α(TH-TH_STD) + β(TH-TH_STD)2] ・・・(8)
【0030】
ここで、第1の測温素子62の温度TIと、第1の測温素子62の抵抗値RIの関係は、下記(9)式で与えられる。
RI = RI_STD×[1+α(TI-TI_STD) + β(TI-TI_STD)2] ・・・(9)
I_STDは第1の測温素子62の標準温度を表し、例えば20℃である。RI_STDは標準温度TI_STDにおける予め計測された第1の測温素子62の抵抗値を表す。上記(9)式より、第1の測温素子62の温度TIは下記(10)式で与えられる。
TI = (1 / 2β)×[-α+ [α2 - 4β(1 - RI / RI_STD)]1/2] + TI_STD ・・・(10)
【0031】
第1の測温素子62は、第1の測温素子62の温度TIが雰囲気ガスの温度TIと同じとみなせるよう、自己発熱しない程度の弱い電圧をかけられて使用される。そのため、発熱素子61と異なり、第1の測温素子62には、高温で発熱するよう高電圧をかけられる機会は少ない。したがって、第1の測温素子62の温度が標準温度TI_STDの場合の第1の測温素子62の抵抗値RI_STDは、予め測定した値から経年変化等によりドリフトする場合が少なく、一定とみなし得る。
【0032】
ここで、第1の測温素子62の抵抗値RIに対する、発熱素子61にエレクトロマイグレーションが生じる前の発熱素子61の抵抗値RHの比γ1は、上記(1)及び(9)式より、下記(11)式で与えられる。
γ1 = RH / RI
= [RH_STD×[1+α(TH-TH_STD) + β(TH-TH_STD)2]] / [RI_STD×[1+α(TI-TI_STD) + β(TI-TI_STD)2]] ・・・(11)
【0033】
発熱素子61及び第1の測温素子62に、発熱素子61及び第1の測温素子62が自己発熱しない程度の弱い同一の電圧をかけた場合、発熱素子61及び第1の測温素子62の温度は、雰囲気ガスの温度とほぼ同じままとなる。そのため、自己発熱しない程度の弱い同一の駆動電力PH0を与えた場合、下記(12)式に示すように、発熱素子61の温度THと、第1の測温素子62の温度TIとは、等しい。
TH = TI ・・・(12)
また、発熱素子61の標準温度TH_STDと、第1の測温素子62の標準温度TI_STDとが、共に20℃等、同じ温度である場合、上記(11)式から下記(13)式が導かれる。
γ1 = RH_STD / RI_STD ・・・(13)
【0034】
次に、第1の測温素子62の抵抗値RIに対する、発熱素子61にエレクトロマイグレーションが生じた後の発熱素子61の抵抗値RH_Dの比γ2は、上記(8)及び(9)式より、下記(14)式で与えられる。
γ2 = RH_D / RI
= [RH_STD_D×[1+α(TH-TH_STD) + β(TH-TH_STD)2]] / [RI_STD×[1+α(TI-TI_STD) + β(TI-TI_STD)2]] ・・・(14)
また、上記(12)及び(13)式より、(14)式から下記(15)式が導かれる。
γ2 = RH_STD_D / RI_STD ・・・(15)
【0035】
下記(16)式に示すように、比γ2を比γ1で割ることにより、ドリフト発生時の発熱素子61の抵抗値の初期値からの変化率Hを算出することが可能である。
H =γ2 1 = RH_D / (RI×γ1) = RH_STD_D / (RI_STD×γ1) = RH_STD_D / RH_STD ・・・(16)
【0036】
発熱素子61にエレクトロマイグレーションが生じた場合の雰囲気ガスの放熱係数MIは、下記(17)式で与えられる。
MI = PH / [(1 / 2β)×[-α+ [α2 - 4β(1 - RH_D / RH_STD_D)]1/2] + TH_STD - TI] ・・・(17)
上記(16)式より、(17)式から下記(18)式が導かれる。
MI = PH / [(1 / 2β)×[-α+ [α2 - 4β(1 - RH_D / (H×RH_STD))]1/2] + TH_STD - TI] ・・・(18)
(18)式を用いることにより、発熱素子61にエレクトロマイグレーションが生じた場合も、雰囲気ガスの放熱係数MIを正確に算出することが可能となる。
【0037】
なお、保温素子64で基板60の温度を一定に保つことにより、発熱素子61が発熱する前のマイクロチップ8の近傍の雰囲気ガスの温度が、基板60の一定の温度と近似する。そのため、発熱素子61が発熱する前の雰囲気ガスの温度の変動が抑制される。温度変動が一度抑制された雰囲気ガスを発熱素子61でさらに加熱することにより、より高い精度で放熱係数MIを算出することが可能となる。
【0038】
ここで、雰囲気ガスが混合ガスであり、混合ガスが、ガスA、ガスB、ガスC、及びガスDの4種類のガス成分からなっていると仮定する。ガスAの体積率VA、ガスBの体積率VB、ガスCの体積率VC、及びガスDの体積率VDの総和は、下記(19)式で与えられるように、1である。
VA+VB+VC+VD=1 ・・・(19)
【0039】
また、ガスAの単位体積当たりの発熱量をKA、ガスBの単位体積当たりの発熱量をKB、ガスCの単位体積当たりの発熱量をKC、ガスDの単位体積当たりの発熱量をKDとすると、混合ガスの単位体積当たりの発熱量Qは、各ガス成分の体積率に、各ガス成分の単位体積当たりの発熱量を乗じたものの総和で与えられる。したがって、混合ガスの単位体積当たりの発熱量Qは、下記(20)式で与えられる。なお、単位体積当たりの発熱量の単位は、例えばMJ/m3である。
Q = KA×VA+ KB×VB+ KC×VC+KD×VD ・・・(20)
【0040】
また、ガスAの放熱係数をMA、ガスBの放熱係数をMB、ガスCの放熱係数をMC、ガスDの放熱係数をMDとすると、混合ガスの放熱係数MIは、各ガス成分の体積率に、各ガス成分の放熱係数を乗じたものの総和で与えられる。したがって、混合ガスの放熱係数MIは、下記(21)式で与えられる。
MI = MA×VA+ MB×VB+ MC×VC+MD×VD ・・・(21)
【0041】
さらに、ガスの放熱係数は発熱素子61の温度THに依存するので、混合ガスの放熱係数MIは、発熱素子61の温度THの関数として、下記(22)式で与えられる。
MI (TH)= MA(TH)×VA+ MB(TH)×VB+ MC(TH)×VC+MD(TH)×VD ・・・(22)
【0042】
したがって、発熱素子61の温度がTH1のときの混合ガスの放熱係数MI1(TH1)は下記(23)式で与えられる。また、発熱素子61の温度がTH2のときの混合ガスの放熱係数MI2(TH2)は下記(24)式で与えられ、発熱素子61の温度がTH3のときの混合ガスの放熱係数MI3(TH3)は下記(25)式で与えられる。
MI1 (TH1)= MA(TH1)×VA+ MB(TH1)×VB+ MC(TH1)×VC+MD(TH1)×VD ・・・(23)
MI2 (TH2)= MA(TH2)×VA+ MB(TH2)×VB+ MC(TH2)×VC+MD(TH2)×VD ・・・(24)
MI3 (TH3)= MA(TH3)×VA+ MB(TH3)×VB+ MC(TH3)×VC+MD(TH3)×VD ・・・(25)
【0043】
ここで、発熱素子61の温度THに対して各ガス成分の放熱係数MA(TH),MB(TH),MC(TH),MD(TH)が非線形性を有する場合、上記(23)乃至(25)式は、線形独立な関係を有する。また、発熱素子61の温度THに対して各ガス成分の放熱係数MA(TH),MB(TH),MC(TH),MD(TH)が線形性を有する場合でも、発熱素子61の温度THに対する各ガス成分の放熱係数MA(TH),MB(TH),MC(TH),MD(TH)の変化率が異なる場合は、上記(23)乃至(25)式は、線形独立な関係を有する。さらに、(23)乃至(25)式が線形独立な関係を有する場合、(19)及び(23)乃至(25)式は線形独立な関係を有する。
【0044】
図5は、天然ガスに含まれるメタン(CH4)、プロパン(C38)、窒素(N2)、及び二酸化炭素(CO2)の放熱係数と、発熱抵抗体である発熱素子61の温度との関係を示すグラフである。発熱素子61の温度に対して、メタン(CH4)、プロパン(C38)、窒素(N2)、及び二酸化炭素(CO2)のそれぞれのガス成分の放熱係数は線形性を有する。しかし、発熱素子61の温度に対する放熱係数の変化率は、メタン(CH4)、プロパン(C38)、窒素(N2)、及び二酸化炭素(CO2)のそれぞれで異なる。したがって、混合ガスを構成するガス成分がメタン(CH4)、プロパン(C38)、窒素(N2)、及び二酸化炭素(CO2)であるである場合、上記(23)乃至(25)式は、線形独立な関係を有する。
【0045】
(23)乃至(25)式中の各ガス成分の放熱係数MA(TH1),MB(TH1),MC(TH1),MD(TH1),MA(TH2),MB(TH2),MC(TH2),MD(TH2),MA(TH3),MB(TH3),MC(TH3),MD(TH3)の値は、計測等により予め得ることが可能である。したがって、(19)及び(23)乃至(25)式の連立方程式を解くと、ガスAの体積率VA、ガスBの体積率VB、ガスCの体積率VC、及びガスDの体積率VDのそれぞれが、下記(26)乃至(29)式に示すように、混合ガスの放熱係数MI1(TH1),MI2(TH2),MI3(TH3)の関数として与えられる。なお、下記(26)乃至(29)式において、nを自然数として、fnは関数を表す記号である。
VA=f1[MI1 (TH1), MI2 (TH2), MI3 (TH3)] ・・・(26)
VB=f2[MI1 (TH1), MI2 (TH2), MI3 (TH3)] ・・・(27)
VC=f3[MI1 (TH1), MI2 (TH2), MI3 (TH3)] ・・・(28)
VD=f4[MI1 (TH1), MI2 (TH2), MI3 (TH3)] ・・・(29)
【0046】
ここで、上記(20)式に(26)乃至(29)式を代入することにより、下記(30)式が得られる。
Q = KA×VA+ KB×VB+ KC×VC+KD×VD
= KA×f1[MI1 (TH1), MI2 (TH2), MI3 (TH3)]
+ KB×f2[MI1 (TH1), MI2 (TH2), MI3 (TH3)]
+ KC×f3[MI1 (TH1), MI2 (TH2), MI3 (TH3)]
+ KD×f4[MI1 (TH1), MI2 (TH2), MI3 (TH3)] ・・・(30)
【0047】
上記(30)式に示すように、混合ガスの単位体積当たりの発熱量Qは、発熱素子61の温度がTH1,TH2,TH3である場合の混合ガスの放熱係数MI1(TH1),MI2(TH2),MI3(TH3)を変数とする方程式で与えられる。したがって、混合ガスの発熱量Qは、gを関数を表す記号として、下記(31)式で与えられる。
Q = g[MI1 (TH1), MI2 (TH2), MI3 (TH3)] ・・・(31)
【0048】
よって、ガスA、ガスB、ガスC、及びガスDからなる混合ガスについて、予め上記(31)式を得れば、ガスAの体積率VA、ガスBの体積率VB、ガスCの体積率VC、及びガスDの体積率VDが未知の検査対象混合ガスの単位体積当たりの発熱量Qを容易に算出可能であることを、発明者らは見出した。具体的には、発熱素子61の発熱温度がTH1,TH2,TH3である場合の検査対象混合ガスの放熱係数MI1(TH1),MI2(TH2),MI3(TH3)を計測し、(31)式に代入することにより、検査対象混合ガスの発熱量Qを一意に求めることが可能となる。
【0049】
なお、混合ガスのガス成分は、4種類に限定されることはない。例えば、混合ガスがn種類のガス成分からなる場合、まず、下記(32)式で与えられる、発熱素子61の少なくともn−1種類の温度TH1,TH2,TH3,・・・,THn-1に対する混合ガスの放熱係数MI1(TH1),MI2(TH2),MI3(TH3),・・・,MIn-1(THn-1)を変数とする方程式を予め取得する。そして、発熱素子61のn−1種類の温度TH1,TH2,TH3,・・・,THn-1に対する、n種類のガス成分のそれぞれの体積率が未知の検査対象混合ガスの放熱係数MI1(TH1),MI2(TH2),MI3(TH3),・・・,MIn-1(THn-1)を計測し、(32)式に代入することにより、検査対象混合ガスの単位体積当たりの発熱量Qを一意に求めることが可能となる。
Q = g[MI1 (TH1), MI2 (TH2), MI3 (TH3), ・・・, MIn-1 (THn-1)] ・・・(32)
【0050】
ただし、混合ガスが、ガス成分としてメタン(CH4)、プロパン(C38)に加えて、jを自然数として、メタン(CH4)とプロパン(C38)以外のアルカン(Cj2j+2)を含む場合、メタン(CH4)とプロパン(C38)以外のアルカン(Cj2j+2)を、メタン(CH4)とプロパン(C38)の混合物とみなしても、(32)式の算出には影響しない。例えば、エタン(C26)、ブタン(C410)、ペンタン(C512)、ヘキサン(C614)を、下記(33)乃至(36)式に示すように、それぞれ所定の係数を掛けられたメタン(CH4)とプロパン(C38)の混合物とみなして(32)式を算出してもかまわない。
C2H6 = 0.5 CH4 + 0.5 C3H8 ・・・(33)
C4H10 = -0.5 CH4 + 1.5 C3H8 ・・・(34)
C5H12 = -1.0 CH4 + 2.0 C3H8 ・・・(35)
C6H14 = -1.5 CH4 + 2.5 C3H8 ・・・(36)
【0051】
したがって、zを自然数として、n種類のガス成分からなる混合ガスが、ガス成分としてメタン(CH4)、プロパン(C38)に加えて、メタン(CH4)とプロパン(C38)以外のz種類のアルカン(Cj2j+2)を含む場合、少なくともn−z−1種類の温度における混合ガスの放熱係数MIを変数とする方程式を求めてもよい。
【0052】
なお、(32)式の算出に用いられた混合ガスのガス成分の種類と、単位体積当たりの発熱量Qが未知の検査対象混合ガスのガス成分の種類が同じ場合に、検査対象混合ガスの発熱量Qの算出に(32)式を利用可能であることはもちろんである。さらに、検査対象混合ガスがn種類より少ない種類のガス成分からなり、かつ、n種類より少ない種類のガス成分が、(32)式の算出に用いられた混合ガスに含まれている場合も、(32)式を利用可能である。例えば、(32)式の算出に用いられた混合ガスが、メタン(CH4)、プロパン(C38)、窒素(N2)、及び二酸化炭素(CO2)の4種類のガス成分を含む場合、検査対象混合ガスが、窒素(N2)を含まず、メタン(CH4)、プロパン(C38)、及び二酸化炭素(CO2)の3種類のガス成分のみを含む場合も、検査対象混合ガスの発熱量Qの算出に(32)式を利用可能である。
【0053】
さらに、(32)式の算出に用いられた混合ガスが、ガス成分としてメタン(CH4)とプロパン(C38)を含む場合、検査対象混合ガスが、(32)式の算出に用いられた混合ガスに含まれていないアルカン(Cj2j+2)を含んでいても、(32)式を利用可能である。これは、上述したように、メタン(CH4)とプロパン(C38)以外のアルカン(Cj2j+2)を、メタン(CH4)とプロパン(C38)の混合物とみなしても、(32)式を用いた単位体積当たりの発熱量Qの算出に影響しないためである。
【0054】
ここで、図6に示す第1の実施の形態に係るガス物性値測定システム20は、発熱量Qの値が既知のサンプル混合ガスが充填されるチャンバ101と、図1及び図2に示す発熱素子61、及び第1の測温素子62を用いて、サンプル混合ガスの複数の放熱係数MIの値を計測する図6に示す計測機構10と、発熱素子61の温度及び第1の測温素子62の温度が等しい場合の第1の測温素子62の抵抗値に基づいて、発熱素子61の抵抗値のドリフトを補正する補正部301と、を備える。さらに、ガス物性値測定システムは、サンプル混合ガスの既知の発熱量Qの値、及びサンプル混合ガスの複数の放熱係数MIの計測値に基づいて、発熱素子61の複数の温度に対するガスの放熱係数MIを独立変数とし、ガスの発熱量Qを従属変数とする発熱量算出式を作成する式作成モジュール303を備える。なお、サンプル混合ガスは、複数種類のガス成分を含む。
【0055】
計測機構10は、サンプル混合ガスが注入されるチャンバ101内に配置された、図1及び図2を用いて説明したマイクロチップ8を備える。マイクロチップ8は、断熱部材18を介してチャンバ101内に配置されている。チャンバ101には、サンプル混合ガスをチャンバ101に送るための流路102と、サンプル混合ガスをチャンバ101から外部に排出するための流路103と、が接続されている。
【0056】
それぞれ発熱量Qが異なる4種類のサンプル混合ガスが使用される場合、図7に示すように、第1のサンプル混合ガスを貯蔵する第1のガスボンベ50A、第2のサンプル混合ガスを貯蔵する第2のガスボンベ50B、第3のサンプル混合ガスを貯蔵する第3のガスボンベ50C、及び第4のサンプル混合ガスを貯蔵する第4のガスボンベ50Dが用意される。第1のガスボンベ50Aには、流路91Aを介して、第1のガスボンベ50Aから例えば0.2MPa等の低圧に調節された第1のサンプル混合ガスを得るための第1のガス圧調節器31Aが接続されている。また、第1のガス圧調節器31Aには、流路92Aを介して、第1の流量制御装置32Aが接続されている。第1の流量制御装置32Aは、流路92A及び流路102を介してガス物性値測定システム20に送られる第1のサンプル混合ガスの流量を制御する。
【0057】
第2のガスボンベ50Bには、流路91Bを介して、第2のガス圧調節器31Bが接続されている。また、第2のガス圧調節器31Bには、流路92Bを介して、第2の流量制御装置32Bが接続されている。第2の流量制御装置32Bは、流路92B,93,102を介してガス物性値測定システム20に送られる第2のサンプル混合ガスの流量を制御する。
【0058】
第3のガスボンベ50Cには、流路91Cを介して、第3のガス圧調節器31Cが接続されている。また、第3のガス圧調節器31Cには、流路92Cを介して、第3の流量制御装置32Cが接続されている。第3の流量制御装置32Cは、流路92C,93,102を介してガス物性値測定システム20に送られる第3のサンプル混合ガスの流量を制御する。
【0059】
第4のガスボンベ50Dには、流路91Dを介して、第4のガス圧調節器31Dが接続されている。また、第4のガス圧調節器31Dには、流路92Dを介して、第4の流量制御装置32Dが接続されている。第4の流量制御装置32Dは、流路92D,93,102を介してガス物性値測定システム20に送られる第4のサンプル混合ガスの流量を制御する。
【0060】
第1乃至第4のサンプル混合ガスのそれぞれは、例えば天然ガスである。第1乃至第4のサンプル混合ガスのそれぞれは、例えばメタン(CH4)、プロパン(C38)、窒素(N2)、及び二酸化炭素(CO2)の4種類のガス成分を含む。
【0061】
図6に示す第1の実施の形態に係るガス物性値測定システム20は、補正情報記憶装置401をさらに備える。補正情報記憶装置401は、上記(13)式で与えられる、第1の測温素子62の抵抗値RIに対する、発熱素子61にエレクトロマイグレーションが生じる前の発熱素子61の抵抗値RHの比γ1の予め取得された値を保存する。
【0062】
図1及び図2に示す発熱素子61及び第1の測温素子62は、それぞれ図6に示す駆動回路310から、自己発熱しない程度の弱い駆動電力PH0を与えられる。補正部301は、上記(2)式より、自己発熱しない程度の駆動電力PH0を与えられた場合の発熱素子61の抵抗値RH_Dと、自己発熱しない程度の駆動電力PH0を与えられた場合の第1の測温素子62の抵抗値RIと、を算出する。また、補正部301は、補正情報記憶装置401から、比γ1の予め取得された値を読み出す。さらに補正部301は、上記(16)式に示すように、発熱素子61の抵抗値RH_Dを、第1の測温素子62の抵抗値RIと比γ1の予め取得された値との積で割って、ドリフト発生時の発熱素子61の抵抗値の初期値からの変化率Hを算出する。またさらに、補正部301は、標準温度TH_STDにおける予め計測された発熱素子61の抵抗値RH_STDと、変化率Hと、の積をとり、標準温度TH_STDにおける発熱素子61の補正された抵抗値H×RH_STDを算出する。
【0063】
次に、チャンバ101に第1のサンプル混合ガスが充填された後、図1及び図2に示すマイクロチップ8の第1の測温素子62は、第1のサンプル混合ガスの温度TIを検出する。その後、図1及び図2に示す発熱素子61は、図6に示す駆動回路310から駆動電力PH1,PH2,PH3を与えられる。駆動電力PH1,PH2,PH3を与えられた場合、発熱素子61は、例えば、100℃の温度TH1、150℃の温度TH2、及び200℃の温度TH3で発熱する。
【0064】
図6に示すチャンバ101から第1のサンプル混合ガスが除去された後、第2乃至第4のサンプル混合ガスがチャンバ101に順次充填される。第2乃至第4のサンプル混合ガスのそれぞれがチャンバ101に充填された後、図1及び図2に示すマイクロチップ8の第1の測温素子62は、第2乃至第4のサンプル混合ガスのそれぞれの温度TIを検出する。また、発熱素子61は、100℃の温度TH1、150℃の温度TH2、及び200℃の温度TH3で発熱する。
【0065】
なお、それぞれのサンプル混合ガスがn種類のガス成分を含む場合、マイクロチップ8の図1及び図2に示す発熱素子61は、少なくともn−1種類の異なる温度で発熱させられる。ただし、上述したように、メタン(CH4)及びプロパン(C38)以外のアルカン(Cj2j+2)は、メタン(CH4)及びプロパン(C38)の混合物とみなしうる。したがって、zを自然数として、n種類のガス成分からなるサンプル混合ガスが、ガス成分としてメタン(CH4)及びプロパン(C38)に加えてz種類のアルカン(Cj2j+2)を含む場合は、発熱素子61は、少なくともn−z−1種類の異なる温度で発熱させられる。
【0066】
図6に示す計測機構10は、マイクロチップ8に接続された放熱係数算出モジュール302を備える。放熱係数算出モジュール302は、上記(18)式を用いて、図1及び図2に示す発熱素子61の温度がTH1(ここでは100℃)の場合の発熱素子61の抵抗値RH1_D、標準温度TH_STDにおける発熱素子61の補正された抵抗値H×RH_STD、及び第1のサンプル混合ガスの温度TIに基づいて、温度TH1が100℃の発熱素子61と熱的に平衡な第1のサンプル混合ガスの放熱係数MI1の値を算出する。計測機構10は、温度TH1が100℃の発熱素子61と熱的に平衡な第2乃至第4のサンプル混合ガスのそれぞれの放熱係数MI1の値も算出する。
【0067】
また、図6に示す放熱係数算出モジュール302は、上記(18)式を用いて、図1及び図2に示す発熱素子61の温度がTH2(ここでは150℃)の場合の発熱素子61の抵抗値RH2、標準温度TH_STDにおける発熱素子61の補正された抵抗値H×RH_STD、及び第1のサンプル混合ガスの温度TIに基づいて、温度TH2が150℃の発熱素子61と熱的に平衡な第1のサンプル混合ガスの放熱係数MI2の値を算出する。計測機構10は、温度TH2が150℃の発熱素子61と熱的に平衡な第2乃至第4のサンプル混合ガスのそれぞれの放熱係数MI2の値も算出する。
【0068】
さらに、図6に示す放熱係数算出モジュール302は、上記(18)式を用いて、図1及び図2に示す発熱素子61の温度がTH3(ここでは200℃)の場合の発熱素子61の抵抗値RH3、標準温度TH_STDにおける発熱素子61の補正された抵抗値H×RH_STD、及び第1のサンプル混合ガスの温度TIに基づいて、温度TH3が200℃の発熱素子61と熱的に平衡な第1のサンプル混合ガスの放熱係数MI3の値を算出する。計測機構10は、温度TH3が200℃の発熱素子61と熱的に平衡な第2乃至第4のサンプル混合ガスのそれぞれの放熱係数MI3の値も算出する。
【0069】
図6に示すガス物性値測定システム20は、CPU300に接続された放熱係数記憶装置402をさらに備える。放熱係数算出モジュール302は、算出した放熱係数MI1,MI2,MI3の値を放熱係数記憶装置402に保存する。
【0070】
式作成モジュール303は、例えば第1乃至第4のサンプル混合ガスのそれぞれの既知の発熱量Qの値と、発熱素子61の温度が100℃の場合のガスの放熱係数MI1の値と、発熱素子61の温度が150℃の場合のガスの放熱係数MI2の値と、発熱素子61の温度が200℃の場合のガスの放熱係数MI3の値と、を収集する。さらに式作成モジュール303は、収集した発熱量Q及び放熱係数MI1,MI2,MI3の値に基づいて、多変量解析により、発熱素子61の温度が100℃の場合の放熱係数MI1、発熱素子61の温度が150℃の場合の放熱係数MI2、及び発熱素子61の温度が200℃の場合の放熱係数MI3を独立変数とし、発熱量Qを従属変数とする発熱量算出式を算出する。
【0071】
なお、「多変量解析」とは、A. J Smola及びB. Scholkopf著の「A Tutorial on Support Vector Regression」(NeuroCOLT Technical Report (NC−TR−98−030)、1998年)に開示されているサポートベクトル回帰、重回帰分析、及び特開平5−141999号公報に開示されているファジィ数量化理論II類等を含む。また、放熱係数算出モジュール302及び式作成モジュール303は、中央演算処理装置(CPU)300に含まれている。
【0072】
ガス物性値測定システム20は、CPU300に接続された式記憶装置403をさらに備える。式記憶装置403は、式作成モジュール303が作成した発熱量算出式を保存する。さらにCPU300には、入力装置312及び出力装置313が接続される。入力装置312としては、例えばキーボード、及びマウス等のポインティングデバイス等が使用可能である。出力装置313には液晶ディスプレイ、モニタ等の画像表示装置、及びプリンタ等が使用可能である。
【0073】
次に、図8に示すフローチャートを用いて第1の実施の形態に係る発熱量算出式の作成方法について説明する。
【0074】
(a)ステップS100で、図7に示す第2乃至第4の流量制御装置32B−32Dの弁を閉じたまま、第1の流量制御装置32Aの弁を開き、図6に示すチャンバ101内に第1のサンプル混合ガスを導入する。ステップS101で、補正部301は、標準温度TH_STDにおける発熱素子61の補正された抵抗値H×RH_STDを算出する。次に、第1の測温素子62は、第1のサンプル混合ガスの温度TIを検出する。その後、図6に示す駆動回路310は、図1及び図2に示す発熱素子61に駆動電力PH1を与え、発熱素子61を100℃で発熱させる。さらに、図6に示す放熱係数算出モジュール302は、上記(18)式を用いて、発熱素子61の温度が100℃の場合の第1のサンプル混合ガスの放熱係数MI1の値を算出する。また、放熱係数算出モジュール302は、発熱素子61の温度が100℃の場合の第1のサンプル混合ガスの放熱係数MI1の値を放熱係数記憶装置402に保存する。
【0075】
(b)ステップS102で、駆動回路310は、図1及び図2に示す発熱素子61の温度の切り替えが完了したか否か判定する。温度150℃及び温度200℃への切り替えが完了していない場合には、ステップS101に戻り、図6に示す駆動回路310は、図1及び図2に示す発熱素子61を150℃で発熱させる。図6に示す放熱係数算出モジュール302は、上記(18)式を用いて、発熱素子61の温度が150℃の場合の第1のサンプル混合ガスの放熱係数MI2の値を算出し、放熱係数記憶装置402に保存する。
【0076】
(c)再びステップS102で、図1及び図2に示す発熱素子61の温度の切り替えが完了したか否か判定する。温度200℃への切り替えが完了していない場合には、ステップS101に戻り、図6に示す駆動回路310は、図1及び図2に示す発熱素子61を200℃で発熱させる。図6に示す放熱係数算出モジュール302は、上記(18)式を用いて、発熱素子61の温度が200℃の場合の第1のサンプル混合ガスの放熱係数MI3の値を算出し、放熱係数記憶装置402に保存する。
【0077】
(d)発熱素子61の温度の切り替えが完了した場合には、ステップS102からステップS103に進む。ステップS103で、サンプル混合ガスの切り替えが完了したか否かを判定する。第2乃至第4のサンプル混合ガスへの切り替えが完了していない場合には、ステップS100に戻る。ステップS100で、図7に示す第1の流量制御装置32Aを閉じ、第3乃至第4の流量制御装置32C−32Dの弁を閉じたまま第2の流量制御装置32Bの弁を開き、図6に示すチャンバ101内に第2のサンプル混合ガスを導入する。
【0078】
(e)第1のサンプル混合ガスと同様に、ステップS101乃至ステップS102のループが繰り返される。放熱係数算出モジュール302は、発熱素子61の温度が100℃の場合の第2のサンプル混合ガスの放熱係数MI1の値、発熱素子61の温度が150℃の場合の第2のサンプル混合ガスの放熱係数MI2の値、及び発熱素子61の温度が200℃の場合の第2のサンプル混合ガスの放熱係数MI3の値を算出する。さらに放熱係数算出モジュール302は、算出した第2のサンプル混合ガスの放熱係数MI1,MI2,MI3の値を放熱係数記憶装置402に保存する。
【0079】
(f)その後、ステップS100乃至ステップS103のループが繰り返される。これにより、発熱素子61の温度が100℃、150℃、200℃の場合の第3のサンプル混合ガスの放熱係数MI1,MI2,MI3の値と、発熱素子61の温度が100℃、150℃、200℃の場合の第4のサンプル混合ガスの放熱係数MI1,MI2,MI3の値とが、放熱係数記憶装置402に保存される。ステップS104で、入力装置312から式作成モジュール303に、第1のサンプル混合ガスの既知の発熱量Qの値、第2のサンプル混合ガスの既知の発熱量Qの値、第3のサンプル混合ガスの既知の発熱量Qの値、及び第4のサンプル混合ガスの既知の発熱量Qの値を入力する。また、式作成モジュール303は、放熱係数記憶装置402から、発熱素子61の温度が100℃、150℃、200℃の場合の第1乃至第4のサンプル混合ガスのそれぞれの放熱係数MI1,MI2,MI3の値を読み出す。
【0080】
(g)ステップS105で、第1乃至第4のサンプル混合ガスの発熱量Qの値と、発熱素子61の温度が100℃、150℃、200℃の場合の第1乃至第4のサンプル混合ガスの放熱係数MI1,MI2,MI3の値と、に基づいて、式作成モジュール303は、重回帰分析を行う。重回帰分析により、式作成モジュール303は、発熱素子61の温度が100℃の場合の放熱係数MI1、発熱素子61の温度が150℃の場合の放熱係数MI2、及び発熱素子61の温度が200℃の場合の放熱係数MI3を独立変数とし、発熱量Qを従属変数とする発熱量算出式を算出する。その後、ステップS106で、式作成モジュール303は作成した発熱量算出式を式記憶装置403に保存し、第1の実施の形態に係る発熱量算出式の作成方法が終了する。
【0081】
以上示したように、第1の実施の形態に係る発熱量算出式の作成方法によれば、計測対象混合ガスの発熱量Qの値を一意に算出可能な発熱量算出式を作成することが可能となる。
【0082】
(第2の実施の形態)
図9に示すように、第2の実施の形態に係るガス物性値測定システム20のCPU300には、熱伝導率記憶装置411が接続されている。ここで、図10は、発熱素子に2mA、2.5mA、及び3mAの電流を流した際の、混合ガスの放熱係数MIと、熱伝導率と、の関係を示す。図10に示すように、混合ガスの放熱係数MIと熱伝導率は一般に比例関係にある。そこで、図9に示す熱伝導率記憶装置411は、チャンバ101に導入されるガスの放熱係数MIと熱伝導率との対応関係を、近似式あるいはテーブル等で予め保存する。
【0083】
第2の実施の形態に係るCPU300は、熱伝導率算出モジュール322をさらに含む。熱伝導率算出モジュール322は、放熱係数記憶装置402から放熱係数MIの値を読み出し、熱伝導率記憶装置411からガスの放熱係数MIと熱伝導率との対応関係を読み出す。さらに熱伝導率算出モジュール322は、ガスの放熱係数MIの値と、ガスの放熱係数MI及び熱伝導率の対応関係とに基づいて、チャンバ101に導入されたガスの熱伝導率を算出する。
【0084】
第2の実施の形態に係るガス物性値測定システム20のその他の構成要素は、第1の実施の形態と同様であるので、説明は省略する。第2の実施の形態に係るガス物性値測定システム20によれば、放熱係数MIに基づいて、ガスの正確な熱伝導率の値を算出することが可能となる。
【0085】
(第3の実施の形態)
図11に示すように、第3の実施の形態に係るガス物性値測定システム20のCPU300には、濃度記憶装置412がさらに接続されている。ここで、図12は、ガス温度TIが0℃、20℃、及び40℃のときのプロパンガスの放熱係数MIと、濃度と、の関係を示す。図12に示すように、ガスの放熱係数MIと、ガスの濃度とは、一般に比例関係にある。そこで、図11に示す濃度記憶装置412は、チャンバ101に導入されるガスの放熱係数MIと、濃度と、の対応関係を、近似式あるいはテーブル等で予め保存する。
【0086】
第3の実施の形態に係るCPU300は、濃度算出モジュール323をさらに含む。濃度算出モジュール323は、放熱係数記憶装置402から放熱係数MIの値を読み出し、濃度記憶装置412からガスの放熱係数MIと、濃度と、の対応関係を読み出す。さらに濃度算出モジュール323は、ガスの放熱係数MIの値と、ガスの放熱係数MI及び濃度の対応関係とに基づいて、チャンバ101に導入されたガスの濃度を算出する。
【0087】
第3の実施の形態に係るガス物性値測定システム20のその他の構成要素は、第1の実施の形態と同様であるので、説明は省略する。第3の実施の形態に係るガス物性値測定システム20によれば、ガスの放熱係数MIに基づいて、ガスの濃度の正確な値を算出することが可能となる。
【0088】
(第4の実施の形態)
図13に示すように、第4の実施の形態に係るガス物性値測定システム21は、発熱量Qの値が未知の計測対象混合ガスが充填されるチャンバ101と、図1及び図2に示す発熱素子61及び第1の測温素子62を用いて、計測対象混合ガスの複数の放熱係数MIの値を計測する図13に示す計測機構10と、を備える。さらに、ガス物性値測定システム21は、発熱素子61の複数の発熱温度に対するガスの放熱係数MIを独立変数とし、発熱量Qを従属変数とする発熱量算出式を保存する式記憶装置403と、発熱量算出式の発熱素子61の複数の発熱温度に対するガスの放熱係数MIの独立変数に、発熱素子61の複数の発熱温度に対する計測対象混合ガスの放熱係数MIの測定値を代入し、計測対象混合ガスの発熱量Qの値を算出する発熱量算出モジュール305と、を備える。
【0089】
式記憶装置403は、第1の実施の形態で説明した発熱量算出式を保存する。ここでは、例として、発熱量算出式の作成のために、メタン(CH4)、プロパン(C38)、窒素(N2)、及び二酸化炭素(CO2)を含む天然ガスがサンプル混合ガスとして使用された場合を説明する。また、発熱量算出式は、発熱素子61の温度TH1が100℃の場合のガスの放熱係数MI1と、発熱素子61の温度TH2が150℃の場合のガスの放熱係数MI2と、発熱素子61の温度TH3が200℃の場合のガスの放熱係数MI3と、を独立変数としているものとする。
【0090】
第4の実施の形態においては、例えば、未知の体積率でメタン(CH4)、プロパン(C38)、窒素(N2)、及び二酸化炭素(CO2)を含む、発熱量Qが未知の天然ガスが、計測対象混合ガスとして、チャンバ101に導入される。発熱素子61及び第1の測温素子62は、それぞれ図6に示す駆動回路310から、自己発熱しない程度の弱い駆動電力PH0を与えられる。補正部301は、上記(2)式より、自己発熱しない程度の駆動電力PH0を与えられた場合の発熱素子61の抵抗値RH_Dと、自己発熱しない程度の駆動電力PH0を与えられた場合の第1の測温素子62の抵抗値RIと、を算出する。また、補正部301は、補正情報記憶装置401から、比γ1の予め取得された値を読み出す。さらに補正部301は、上記(16)式に示すように、発熱素子61の抵抗値RH_Dを、第1の測温素子62の抵抗値RIと比γ1の予め取得された値との積で割って、ドリフト発生時の発熱素子61の抵抗値の初期値からの変化率Hを算出する。またさらに、補正部301は、標準温度TH_STDにおける予め計測された発熱素子61の抵抗値RH_STDと、変化率Hと、の積をとり、標準温度TH_STDにおける発熱素子61の補正された抵抗値H×RH_STDを算出する。
【0091】
次に、図1及び図2に示すマイクロチップ8の第1の測温素子62は、計測対象混合ガスの温度TIを検出する。また、発熱素子61は、図6に示す駆動回路310から駆動電力PH1,PH2,PH3を与えられる。駆動電力PH1,PH2,PH3を与えられた場合、発熱素子61は、例えば、100℃の温度TH1、150℃の温度TH2、及び200℃の温度TH3で発熱する。
【0092】
図13に示す放熱係数算出モジュール302は、上記(18)式を用いて、図1及び図2に示す発熱素子61の温度がTH1(ここでは100℃)の場合の発熱素子61の抵抗値RH1_D、標準温度TH_STDにおける発熱素子61の補正された抵抗値H×RH_STD、及び計測対象混合ガスの温度TIに基づいて、温度100℃で発熱する発熱素子61と熱的に平衡な計測対象混合ガスの放熱係数MI1の値を算出する。また、放熱係数算出モジュール302は、温度150℃で発熱する発熱素子61と熱的に平衡な計測対象混合ガスの放熱係数MI2の値、及び温度200℃で発熱する発熱素子61と熱的に平衡な計測対象混合ガスの放熱係数MI3の値を算出する。放熱係数算出モジュール302は、算出した放熱係数MI1,MI2,MI3の値を放熱係数記憶装置402に保存する。
【0093】
発熱量算出モジュール305は、発熱量算出式のガスの放熱係数MI1,MI2,MI3の独立変数に、計測対象混合ガスの放熱係数MI1,MI2,MI3の測定値を代入し、計測対象混合ガスの発熱量Qの値を算出する。CPU300には、発熱量記憶装置404がさらに接続されている。発熱量記憶装置404は、発熱量算出モジュール305が算出した計測対象混合ガスの発熱量Qの値を保存する。第4の実施の形態に係るガス物性値測定システム21のその他の構成要件は、図6で説明した第1の実施の形態に係るガス物性値測定システム20と同様であるので、説明は省略する。
【0094】
次に、図14に示すフローチャートを用いて、第4の実施の形態に係る発熱量の測定方法について説明する。
【0095】
(a)ステップS200で、図13に示すチャンバ101内に計測対象混合ガスを導入する。ステップS201で、補正部301は、標準温度TH_STDにおける発熱素子61の補正された抵抗値H×RH_STDを算出する。次に、第1の測温素子62は、計測対象混合ガスの温度TIを検出する。その後、図13に示す駆動回路310は、図1及び図2に示す発熱素子61に駆動電力PH1を与え、発熱素子61を100℃で発熱させる。さらに、図13に示す放熱係数算出モジュール302は、上記(18)式を用いて、発熱素子61の温度が100℃の場合の計測対象混合ガスの放熱係数MI1の値を算出する。また、放熱係数算出モジュール302は、発熱素子61の温度が100℃の場合の計測対象混合ガスの放熱係数MI1の値を放熱係数記憶装置402に保存する。
【0096】
(b)ステップS202で、図13に示す駆動回路310は、図1及び図2に示す発熱素子61の温度の切り替えが完了したか否か判定する。温度150℃及び温度200℃への切り替えが完了していない場合には、ステップS201に戻り、図13に示す駆動回路310は、図1及び図2に示す発熱素子61を150℃に発熱させる。図13に示す放熱係数算出モジュール302は、上記(18)式を用いて、発熱素子61の温度が150℃の場合の計測対象混合ガスの放熱係数MI2の値を算出し、放熱係数記憶装置402に保存する。
【0097】
(c)再びステップS202で、図1及び図2に示す発熱素子61の温度の切り替えが完了したか否か判定する。温度200℃への切り替えが完了していない場合には、ステップS201に戻り、図13に示す駆動回路310は、図1及び図2に示す発熱素子61を200℃に発熱させる。図13に示す放熱係数算出モジュール302は、上記(18)式を用いて、発熱素子61の温度が200℃の場合の計測対象混合ガスの放熱係数MI3の値を算出し、放熱係数記憶装置402に保存する。
【0098】
(d)発熱素子61の温度の切り替えが完了した場合には、ステップS202からステップS203に進む。ステップS203で、図13に示す発熱量算出モジュール305は、式記憶装置403から、発熱素子61の温度が100℃、150℃、及び200℃の場合のガスの放熱係数MI1,MI2,MI3を独立変数とする発熱量算出式を読み出す。また、発熱量算出モジュール305は、放熱係数記憶装置402から、発熱素子61の温度が100℃、150℃、及び200℃の場合の計測対象混合ガスの放熱係数MI1,MI2,MI3の計測値を読み出す。
【0099】
(e)ステップS204で、発熱量算出モジュール305は、発熱量算出式の放熱係数MI1,MI2,MI3の独立変数に計測対象混合ガスの放熱係数MI1,MI2,MI3の計測値を代入し、計測対象混合ガスの発熱量Qの値を算出する。その後、発熱量算出モジュール305は、算出した発熱量Qの値を発熱量記憶装置404に保存し、第4の実施の形態に係る発熱量の測定方法を終了する。
【0100】
以上説明した第4の実施の形態に係る発熱量算出方法によれば、高価なガスクロマトグラフィ装置や音速センサを用いることなく、計測対象混合ガスの放熱係数MI1,MI2,MI3の計測値から、計測対象混合ガスの混合ガスの発熱量Qの値を測定することが可能となる。
【0101】
天然ガスは、産出ガス田によって炭化水素の成分比率が異なる。また、天然ガスには、炭化水素の他に、窒素(N2)や炭酸ガス(CO2)等が含まれる。そのため、産出ガス田によって、天然ガスに含まれるガス成分の体積率は異なり、ガス成分の種類が既知であっても、天然ガスの発熱量Qは未知であることが多い。また、同一のガス田由来の天然ガスであっても、発熱量Qが常に一定であるとは限らず、採取時期によって変化することもある。
【0102】
従来、天然ガスの使用料金を徴収する際には、天然ガスの使用発熱量Qでなく、使用体積に応じて課金する方法がとられている。しかし、天然ガスは由来する産出ガス田によって発熱量Qが異なるため、使用体積に課金するのは公平でない。これに対し、第4の実施の形態に係る発熱量算出方法を用いれば、ガス成分の種類が既知であるが、ガス成分の体積率が未知であるために発熱量Qが未知の天然ガス等の混合ガスの発熱量Qを、簡易に算出することが可能となる。そのため、公平な使用料金を徴収することが可能となる。
【0103】
また、ガラス加工品の製造業においては、ガラスを加熱加工する際、加工精度を一定に保つために、一定の発熱量Qを有する天然ガスが供給されることが望まれている。そのためには、複数のガス田由来の天然ガスのそれぞれの発熱量Qを正確に把握し、総ての天然ガスの発熱量Qが同一になるよう調整した上で、ガラスの加熱加工工程に天然ガスを供給することが検討されている。これに対し、第4の実施の形態に係る発熱量算出方法を用いれば、複数のガス田由来の天然ガスのそれぞれ発熱量Qを正確には把握することが可能となるため、ガラスの加熱加工精度を一定に保つことが可能となる。
【0104】
さらに、第4の実施の形態に係る発熱量算出方法によれば、天然ガス等の混合ガスの正確な発熱量Qを容易に知ることが可能となるため、混合ガスを燃焼させる場合に必要な空気量を適切に設定することが可能となる。そのため、無駄な二酸化炭素(CO2)の排出量を削減することも可能となる。
【0105】
(実施例1)
まず、図15に示すように発熱量Qの値が既知の28種類のサンプル混合ガスを用意した。28種類のサンプル混合ガスのそれぞれは、ガス成分としてメタン(CH4)、エタン(C26)、プロパン(C38)、ブタン(C410)、窒素(N2)、及び二酸化炭素(CO2)のいずれか又は全部を含んでいた。例えば、No.7のサンプル混合ガスは、90vol%のメタン、3vol%のエタン、1vol%のプロパン、1vol%のブタン、4vol%の窒素、及び1vol%の二酸化炭素を含んでいた。また、No.8のサンプル混合ガスは、85vol%のメタン、10vol%のエタン、3vol%のプロパン、及び2vol%のブタンを含み、窒素及び二酸化炭素を含んでいなかった。また、No.9のサンプル混合ガスは、85vol%のメタン、8vol%のエタン、2vol%のプロパン、1vol%のブタン、2vol%の窒素、及び2vol%の二酸化炭素を含んでいた。
【0106】
次に、28種類のサンプル混合ガスのそれぞれの放熱係数MIの値を、発熱素子の温度を100℃、150℃、及び200℃に設定して計測した。なお、例えばNo.7のサンプル混合ガスは6種類のガス成分を含んでいるが、上述したように、エタン(C26)とブタン(C410)は、メタン(CH4)とプロパン(C38)の混合物とみなしうるので、放熱係数MIの値を3種類の温度で計測しても問題ない。その後、28種類のサンプル混合ガスの発熱量Qの値と、計測された放熱係数MIの値に基づいて、サポートベクトル回帰により、放熱係数MIを独立変数とし、発熱量Qを従属変数とする、発熱量Qを算出するための1次方程式、2次方程式、及び3次方程式を作成した。
【0107】
発熱量Qを算出するための1次方程式を作成する際には、キャリブレーション・ポイントは、3乃至5個を目安に、適宜決定できる。作成された1次方程式は下記(37)式で与えられた。28種類のサンプル混合ガスの発熱量Qを(37)式で算出し、真の発熱量Qと比較したところ、最大誤差は2.1%であった。
Q = 39.91 - 20.59×MI (100℃) - 0.89×MI (150℃) + 19.73×MI (200℃) ・・・(37)
【0108】
発熱量Qを算出するための2次方程式を作成する際には、キャリブレーション・ポイントは、8乃至9個を目安に、適宜決定できる。28種類のサンプル混合ガスの発熱量Qを作成された2次方程式で算出し、真の発熱量Qと比較したところ、最大誤差は1.2乃至1.4%であった。
【0109】
発熱量Qを算出するための3次方程式を作成する際には、キャリブレーション・ポイントは、10乃至14個を目安に、適宜決定できる。28種類のサンプル混合ガスの発熱量Qを作成された3次方程式で算出し、真の発熱量Qと比較したところ、最大誤差は1.2%未満であった。図16及び図17に示すように、10個のキャリブレーション・ポイントを取って作成された3次方程式で算出された発熱量Qは、真の発熱量Qに良好に近似した。
【符号の説明】
【0110】
8 マイクロチップ
10 計測機構
18 断熱部材
20,21 ガス物性値測定システム
31A,31B,31C,31D ガス圧調節器
32A,32B,32C,32D 流量制御装置
50A,50B,50C,50D ガスボンベ
60 基板
61 発熱素子
62 第1の測温素子
63 第2の測温素子
64 保温素子
65 絶縁膜
66 キャビティ
91A,91B,91C,91D,92A,92B,92C,92D,93,102,103 流路
101 チャンバ
160,161,162,163,164,165,166,181,182,183 抵抗素子
170,171 オペアンプ
301 補正部
302 放熱係数算出モジュール
303 式作成モジュール
305 発熱量算出モジュール
310 駆動回路
312 入力装置
313 出力装置
322 熱伝導率算出モジュール
323 濃度算出モジュール
401 補正情報記憶装置
402 放熱係数記憶装置
403 式記憶装置
404 発熱量記憶装置
411 熱伝導率記憶装置
412 濃度記憶装置

【特許請求の範囲】
【請求項1】
複数の混合ガスが注入される容器と、
前記容器に配置された発熱素子と、
前記容器に配置され、前記複数の混合ガスのそれぞれの温度を計測する測温素子と、
前記発熱素子の温度及び前記測温素子の温度が等しい場合の前記測温素子の抵抗値に基づいて、前記発熱素子の抵抗値のドリフトを補正する補正部と、
前記複数の混合ガスのそれぞれの計測された温度に基づいて、前記発熱素子を複数の発熱温度で発熱させたときの前記複数の混合ガスのそれぞれの放熱係数の値を計測する計測機構と、
前記複数の混合ガスの既知の発熱量の値と、前記複数の発熱温度に対する前記放熱係数の計測値と、に基づいて、前記複数の発熱温度に対する放熱係数を独立変数とし、前記発熱量を従属変数とする発熱量算出式を作成する式作成モジュールと、
を備える、発熱量算出式作成システム。
【請求項2】
前記発熱素子及び前記測温素子が同一の部材からなる、請求項1の発熱量算出式作成システム。
【請求項3】
前記発熱素子の複数の発熱温度の数が、少なくとも、前記複数の混合ガスのそれぞれが含むガス成分の数から1を引いた数である、請求項1又は2に記載の発熱量算出式作成システム。
【請求項4】
前記式作成モジュールが、サポートベクトル回帰を用いて前記発熱量算出式を作成する、請求項1乃至3のいずれか1項に記載の発熱量算出式作成システム。
【請求項5】
前記複数の混合ガスのそれぞれが天然ガスである、請求項1乃至4のいずれか1項に記載の発熱量算出式作成システム。
【請求項6】
複数の混合ガスを準備することと、
前記複数の混合ガスのそれぞれの温度を測温素子で計測することと、
前記発熱素子の温度及び前記測温素子の温度が等しい場合の前記測温素子の抵抗値に基づいて、前記発熱素子の抵抗値のドリフトを補正することと、
前記複数の混合ガスのそれぞれの計測された温度に基づいて、発熱素子を複数の発熱温度で発熱させたときの前記複数の混合ガスのそれぞれの放熱係数の値を計測することと、
前記複数の混合ガスの既知の発熱量の値と、前記複数の発熱温度に対する前記放熱係数の計測値と、に基づいて、前記複数の発熱温度に対する放熱係数を独立変数とし、前記発熱量を従属変数とする発熱量算出式を作成することと、
を含む、発熱量算出式の作成方法。
【請求項7】
前記発熱素子及び前記測温素子が同一の部材からなる、請求項6の発熱量算出式の作成方法。
【請求項8】
前記複数の発熱温度の数が、少なくとも、前記複数の混合ガスのそれぞれが含むガス成分の数から1を引いた数である、請求項6又は7に記載の発熱量算出式の作成方法。
【請求項9】
前記発熱量算出式を作成することにおいて、サポートベクトル回帰が用いられる、請求項6乃至8のいずれか1項に記載の発熱量算出式の作成方法。
【請求項10】
前記複数の混合ガスのそれぞれが天然ガスである、請求項6乃至9のいずれか1項に記載の発熱量算出式の作成方法。
【請求項11】
発熱量が未知の計測対象混合ガスが注入される容器と、
前記容器に配置された発熱素子と、
前記容器に配置され、前記計測対象混合ガスの温度を計測する測温素子と、
前記発熱素子の温度及び前記測温素子の温度が等しい場合の前記測温素子の抵抗値に基づいて、前記発熱素子の抵抗値のドリフトを補正する補正部と、
前記計測対象混合ガスの計測された温度に基づいて、前記発熱素子を複数の発熱温度で発熱させたときの前記計測対象混合ガスの放熱係数の計測値を計測する計測機構と、
前記複数の発熱温度に対する前記放熱係数を独立変数とし、前記発熱量を従属変数とする発熱量算出式を保存する式記憶装置と、
前記発熱量算出式の前記放熱係数の独立変数に前記計測対象混合ガスの前記放熱係数の計測値を代入し、前記計測対象混合ガスの発熱量の値を算出する発熱量算出モジュールと、
を備える、発熱量測定システム。
【請求項12】
前記発熱素子及び前記測温素子が同一の部材からなる、請求項11の発熱量測定システム。
【請求項13】
前記複数の発熱温度の数が、少なくとも、前記計測対象混合ガスに含まれる複数種類のガス成分の数から1を引いた数である、請求項11又は12に記載の発熱量測定システム。
【請求項14】
複数種類のガス成分を含む複数のサンプル混合ガスの発熱量の値と、前記複数の発熱温度に対する前記複数のサンプル混合ガスの放熱係数の値とに基づいて、前記複数の発熱温度に対する放熱係数を独立変数とし、前記発熱量を従属変数とする発熱量算出式が作成された、請求項11乃至13のいずれか1項に記載の発熱量測定システム。
【請求項15】
前記発熱量算出式を作成するために、サポートベクトル回帰が用いられた、請求項14に記載の発熱量測定システム。
【請求項16】
前記複数のサンプル混合ガスのそれぞれが天然ガスである、請求項14又は15に記載の発熱量測定システム。
【請求項17】
前記計測対象混合ガスが天然ガスである、請求項11乃至16のいずれか1項に記載の発熱量測定システム。
【請求項18】
発熱量が未知の計測対象混合ガスを用意することと、
前記計測対象混合ガスの温度を測温素子で計測することと、
前記発熱素子の温度及び前記測温素子の温度が等しい場合の前記測温素子の抵抗値に基づいて、前記発熱素子の抵抗値のドリフトを補正することと、
前記計測対象混合ガスの計測された温度に基づいて、発熱素子を複数の発熱温度で発熱させたときの前記計測対象混合ガスの放熱係数の計測値を計測することと、
前記複数の発熱温度に対する前記放熱係数を独立変数とし、前記発熱量を従属変数とする発熱量算出式を用意することと、
前記発熱量算出式の前記放熱係数の独立変数に前記計測対象混合ガスの前記放熱係数の計測値を代入し、前記計測対象混合ガスの発熱量の値を算出することと、
を含む、発熱量の測定方法。
【請求項19】
前記発熱素子及び前記測温素子が同一の部材からなる、請求項18の発熱量の測定方法。
【請求項20】
前記複数の温度の数が、少なくとも、前記計測対象混合ガスに含まれる複数種類のガス成分の数から1を引いた数である、請求項18又は19に記載の発熱量の測定方法。
【請求項21】
複数種類のガス成分を含む複数のサンプル混合ガスの発熱量の値と、前記複数の発熱温度に対する前記複数のサンプル混合ガスの放熱係数の値とに基づいて、前記複数の発熱温度に対する放熱係数を独立変数とし、前記発熱量を従属変数とする発熱量算出式が作成された、請求項18乃至20のいずれか1項に記載の発熱量の測定方法。
【請求項22】
前記発熱量算出式を作成するために、サポートベクトル回帰が用いられた、請求項21に記載の発熱量の測定方法。
【請求項23】
前記複数のサンプル混合ガスのそれぞれが天然ガスである、請求項21又は22に記載の発熱量の測定方法。
【請求項24】
前記計測対象混合ガスが天然ガスである、請求項18乃至23のいずれか1項に記載の発熱量の測定方法。
【請求項25】
ガスが注入される容器と、
前記容器に配置された発熱素子と、
前記容器に配置され、前記ガスの温度を計測する測温素子と、
前記発熱素子の温度及び前記測温素子の温度が等しい場合の前記測温素子の抵抗値に基づいて、前記発熱素子の抵抗値のドリフトを補正する補正部と、
前記ガスの計測された温度に基づいて、前記発熱素子を発熱させたときの前記計測対象混合ガスの放熱係数の値を計測する計測機構と、
を備える、物性計測システム。
【請求項26】
前記計測された放熱係数の値に基づき、前記ガスの熱伝導率を算出する熱伝導率算出モジュールを更に備える請求項25に記載の物性計測システム。
【請求項27】
前記計測された放熱係数の値に基づき、前記ガスの濃度を算出する濃度算出モジュールを更に備える請求項26に記載の物性計測システム。
【請求項28】
ガスを準備することと、
前記ガスの温度を測温素子で計測することと、
発熱素子の温度及び前記測温素子の温度が等しい場合の前記測温素子の抵抗値に基づいて、前記発熱素子の抵抗値のドリフトを補正することと、
前記ガスの計測された温度に基づいて、前記発熱素子を発熱させたときの前記ガスの放熱係数の値を計測することと、
を含む、物性計測方法。
【請求項29】
前記計測された放熱係数の値に基づき、前記ガスの熱伝導率を算出することを更に含む請求項28に記載の物性計測方法。
【請求項30】
前記計測された放熱係数の値に基づき、前記ガスの濃度を算出することを更に含む請求項29に記載の物性計測方法。

【図1】
image rotate

【図2】
image rotate

【図3】
image rotate

【図4】
image rotate

【図5】
image rotate

【図6】
image rotate

【図7】
image rotate

【図8】
image rotate

【図9】
image rotate

【図10】
image rotate

【図11】
image rotate

【図12】
image rotate

【図13】
image rotate

【図14】
image rotate

【図15】
image rotate

【図16】
image rotate

【図17】
image rotate