"ガスに接する測温素子と、前記ガスに接し、複数の発熱温度で発熱する発熱素子と、前記ガスの温度に依存する前記測温素子からの電気信号の値と、前記複数の発熱温度のそれぞれにおける前記発熱素子からの電気信号の値と、を測定する測定部と、前記測温素子からの電気信号及び前記複数の発熱温度における前記発熱素子からの電気信号を表す独立変数と、ガスに含まれるそれぞれの分子を構成する水素原子の数と前記水素原子の原子量と前記それぞれの分子の体積比との積の総和を表す従属変数と、を含む第1の式を保存する式記憶装置と、前記第1の式の独立変数に、前記測温素子からの電気信号の値、及び前記発熱素子からの電気信号の値を代入し、前記ガスに含まれるそれぞれの分子を構成する水素原子の数と前記水素原子の原子量と前記それぞれの分子の体積比との積の総和を表す値を算出する算出部と、前記ガスから抽出された水素を供給される燃料電池と、前記算出された積を表す値に基づき、前記燃料電池への前記水素の供給量を制御する制御装置と、を備える、発電システム。前記式記憶装置が、前記測温素子からの電気信号及び前記複数の発熱温度における前記発熱素子からの電気信号を表す独立変数と、ガスに含まれるそれぞれの分子を構成する炭素原子の数と前記炭素原子の原子量と前記それぞれの分子の体積比との積の総和を表す従属変数と、を含む第2の式を更に保存し、前記算出部が、前記第2の式の独立変数に、前記測温素子からの電気信号の値、及び前記発熱素子からの電気信号の値を代入し、前記ガスに含まれるそれぞれの分子を構成する炭素原子の数と前記炭素原子の原子量と前記それぞれの分子の体積比との積の総和を表す値を算出する、請求項1に記載の発電システム。前記ガスが天然ガス又は都市ガスである、請求項1又は2に記載の発電システム。ガスに接する測温素子と、前記ガスに接し、複数の発熱温度で発熱する発熱素子と、前記ガスの温度に依存する前記測温素子からの電気信号の値と、前記複数の発熱温度のそれぞれにおける前記発熱素子からの電気信号の値と、を測定する測定部と、前記測温素子からの電気信号及び前記複数の発熱温度における前記発熱素子からの電気信号を表す独立変数と、ガスに含まれるそれぞれの分子を構成する原子の数と前記原子の原子量と前記それぞれの分子の体積比との積の総和を表す従属変数と、を含む式を保存する式記憶装置と、前記式の独立変数に、前記測温素子からの電気信号の値、及び前記発熱素子からの電気信号の値を代入し、前記ガスに含まれるそれぞれの分子を構成する原子の数と前記原子の原子量と前記それぞれの分子の体積比との積の総和を表す値を算出する算出部と、を備える、ガス測定システム。前記ガスに含まれるそれぞれの分子を構成する原子の数と、前記原子の原子量と、前記それぞれの分子の体積比と、の積の総和が、前記ガスに含まれるそれぞれの分子を構成する水素原子の数と、前記水素原子の原子量と、前記それぞれの分子の体積比と、の積の総和である、請求項4に記載のガス測定システム。前記ガスに含まれるそれぞれの分子を構成する原子の数と、前記原子の原子量と、前記それぞれの分子の体積比と、の積の総和が、前記ガスに含まれるそれぞれの分子を構成する炭素原子の数と、前記炭素原子の原子量と、前記それぞれの分子の体積比と、の積の総和である、請求項4に記載のガス測定システム。前記ガスに含まれるそれぞれの分子を構成する原子の数と、前記原子の原子量と、前記それぞれの分子の体積比と、の積の総和が、前記ガスに含まれるそれぞれの分子を構成する水素原子の数と、前記水素原子の原子量と、前記それぞれの分子の体積比と、の積の総和と、前記ガスに含まれるそれぞれの分子を構成する炭素原子の数と、前記炭素原子の原子量と、前記それぞれの分子の体積比と、の積の総和と、の和である、請求項4に記載のガス測定システム。前記ガスが天然ガス又は都市ガスである、請求項4乃至7のいずれか1項に記載のガス測定システム。ガスの放熱係数又は熱伝導率の測定値を測定する測定部と、前記放熱係数又は熱伝導率と、ガスに含まれるそれぞれの分子を構成する原子の数と前記原子の原子量と前記それぞれの分子の体積比との積の総和と、の相関関係を保存する記憶装置と、前記ガスの放熱係数又は熱伝導率の測定値と、前記相関関係と、に基づいて、前記ガスに含まれるそれぞれの分子を構成する原子の数と前記原子の原子量と前記それぞれの分子の体積比との積の総和を表す値を算出する算出部と、を備える、ガス測定システム。前記ガスに含まれるそれぞれの分子を構成する原子の数と、前記原子の原子量と、前記それぞれの分子の体積比と、の積の総和が、前記ガスに含まれるそれぞれの分子を構成する水素原子の数と、前記水素原子の原子量と、前記それぞれの分子の体積比と、の積の総和である、請求項9に記載のガス測定システム。前記ガスに含まれるそれぞれの分子を構成する原子の数と、前記原子の原子量と、前記それぞれの分子の体積比と、の積の総和が、前記ガスに含まれるそれぞれの分子を構成する炭素原子の数と、前記炭素原子の原子量と、前記それぞれの分子の体積比と、の積の総和である、請求項9に記載のガス測定システム。前記ガスに含まれるそれぞれの分子を構成する原子の数と、前記原子の原子量と、前記それぞれの分子の体積比と、の積の総和が、前記ガスに含まれるそれぞれの分子を構成する水素原子の数と、前記水素原子の原子量と、前記それぞれの分子の体積比と、の積の総和と、前記ガスに含まれるそれぞれの分子を構成する炭素原子の数と、前記炭素原子の原子量と、前記それぞれの分子の体積比と、の積の総和と、の和である、請求項9に記載のガス測定システム。前記ガスが天然ガス又は都市ガスである、請求項9乃至12のいずれか1項に記載のガス測定システム。上述したように、固体高分子形燃料電池は、燃料として供給された水素を、酸素と化学反応させて発電する。そのため、燃料電池に供給される水素の量を正確に計測可能な技術が求められている。また、燃料電池に限らず、様々な分野において、ガスに含まれる分子を構成する原子の量を正確に計測可能な技術が求められている。そこで、本発明は、ガスに含まれる分子を構成する原子の数を評価可能な発電システム及びガス測定システムを提供することを目的の一つとする。本発明の態様によれば、aガスに接する測温素子と、bガスに接し、複数の発熱温度で発熱する発熱素子と、cガスの温度に依存する測温素子からの電気信号の値と、複数の発熱温度のそれぞれにおける発熱素子からの電気信号の値と、を測定する測定部と、d測温素子からの電気信号及び複数の発熱温度における発熱素子からの電気信号を表す独立変数と、ガスに含まれるそれぞれの分子を構成する水素原子の数と水素原子の原子量とそれぞれの分子の体積比との積の総和を表す従属変数と、を含む第1の式を保存する式記憶装置と、e第1の式の独立変数に、測温素子からの電気信号の値、及び発熱素子からの電気信号の値を代入し、ガスに含まれるそれぞれの分子を構成する水素原子の数と水素原子の原子量とそれぞれの分子の体積比との積の総和を表す値を算出する算出部と、fガスから抽出された水素を供給される燃料電池と、g算出された積を表す値に基づき、燃料電池への水素の供給量を制御する制御装置と、を備える、発電システムが提供される。本発明の態様によれば、aガスに接する測温素子と、bガスに接し、複数の発熱温度で発熱する発熱素子と、cガスの温度に依存する測温素子からの電気信号の値と、複数の発熱温度のそれぞれにおける発熱素子からの電気信号の値と、を測定する測定部と、d測温素子からの電気信号及び複数の発熱温度における発熱素子からの電気信号を表す独立変数と、ガスに含まれるそれぞれの分子を構成する原子の数と原子の原子量とそれぞれの分子の体積比との積の総和を表す従属変数と、を含む式を保存する式記憶装置と、e式の独立変数に、測温素子からの電気信号の値、及び発熱素子からの電気信号の値を代入し、ガスに含まれるそれぞれの分子を構成する原子の数と原子の原子量とそれぞれの分子の体積比との積の総和を表す値を算出する算出部と、を備える、ガス測定システムが提供される。本発明の態様によればaガスの放熱係数又は熱伝導率の測定値を測定する測定部と、b放熱係数又は熱伝導率と、ガスに含まれるそれぞれの分子を構成する原子の数と原子の原子量とそれぞれの分子の体積比との積の総和と、の相関関係を保存する記憶装置と、cガスの放熱係数又は熱伝導率の測定値と、相関関係と、に基づいて、ガスに含まれるそれぞれの分子を構成する原子の数と原子の原子量とそれぞれの分子の体積比との積の総和を表す値を算出する算出部と、を備える、ガス測定システムが提供される。本発明によれば、ガスに含まれる分子を構成する原子の数を評価可能な発電システム及びガス測定システムを提供可能である。以下に本発明の実施の形態を説明する。以下の図面の記載において、同一又は類似の部分には同一又は類似の符号で表している。但し、図面は模式的なものである。したがって、具体的な寸法等は以下の説明を照らし合わせて判断するべきものである。また、図面相互間においても互いの寸法の関係や比率が異なる部分が含まれていることは勿論である。第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と、を備える。ダイアフラムには、複数の孔が設けられている。ダイアフラムに複数の孔を設けることにより、キャビティ66内のガスの置換が速くなる。あるいは、絶縁膜65は、図3、及びIV-IV方向から見た断面図である図4に示すように、キャビティ66をブリッジ状に覆うように、基板60上に配置されてもよい。これによっても、キャビティ66内が露出し、キャビティ66内のガスの置換が速くなる。発熱素子61は、キャビティ66を覆う絶縁膜65のダイアフラムの部分の中心に配置されている。発熱素子61は、例えば抵抗器であり、電力を与えられて発熱し、発熱素子61に接する雰囲気ガスを加熱する。第1の測温素子62及び第2の測温素子63は、例えば抵抗器等の受動素子等の電子素子であり、雰囲気ガスのガス温度に依存した電気信号を出力する。以下においては、第1の測温素子62の出力信号を利用する例を説明するが、これに限定されず、例えば第1の測温素子62の出力信号及び第2の測温素子63の出力信号の平均値を、測温素子の出力信号として利用してもよい。保温素子64は、例えば抵抗器であり、電力を与えられて発熱し、基板60の温度を一定に保つ。基板60の材料としては、シリコンSi等が使用可能である。絶縁膜65の材料としては、酸化ケイ素SiO2等が使用可能である。キャビティ66は、異方性エッチング等により形成される。また発熱素子61、第1の測温素子62、第2の測温素子63、及び保温素子64のそれぞれの材料には白金Pt等が使用可能であり、リソグラフィ法等により形成可能である。また、発熱素子61、第1の測温素子62、及び第2の測温素子63は、同一の部材からなっていてもよい。マイクロチップ8は、マイクロチップ8の底面に配置された断熱部材18を介して、雰囲気ガスが流されるパイプ等に固定される。断熱部材18を介してマイクロチップ8をパイプに固定することにより、マイクロチップ8の温度が、パイプの内壁の温度変動の影響を受けにくくなる。ガラス等からなる断熱部材18の熱伝導率は、例えば1.0W/mK以下である。図5に示すように、発熱素子61の一端には、例えば、オペアンプ170の入力端子が電気的に接続され、他端は接地される。また、オペアンプ170の入力端子及び出力端子と並列に、抵抗素子161が接続される。オペアンプ170の-入力端子は、直列に接続された抵抗素子162と抵抗素子163との間、直列に接続された抵抗素子163と抵抗素子164との間、直列に接続された抵抗素子164と抵抗素子165との間、又は抵抗素子165の接地端子に電気的に接続される。各抵抗素子162-165の抵抗値を適当に定めることにより、例えば5.0Vの電圧Vinを抵抗素子162の一端に印加すると、抵抗素子163と抵抗素子162との間には、例えば2.4Vの電圧VL3が生じる。また、抵抗素子164と抵抗素子163との間には、例えば1.9Vの電圧VL2が生じ、抵抗素子165と抵抗素子164との間には、例えば1.4Vの電圧VL1が生じる。抵抗素子162及び抵抗素子163の間と、オペアンプの-入力端子との間には、スイッチSW1が設けられており、抵抗素子163及び抵抗素子164の間と、オペアンプの-入力端子との間には、スイッチSW2が設けられている。また、抵抗素子164及び抵抗素子165の間と、オペアンプの-入力端子との間には、スイッチSW3が設けられており、抵抗素子165の接地端子と、オペアンプの-入力端子との間には、スイッチSW4が設けられている。オペアンプ170の-入力端子に2.4Vの電圧VL3を印加する場合、スイッチSW1のみが通電され、スイッチSW2,SW3,SW4は切断される。オペアンプ170の-入力端子に1.9Vの電圧VL2を印加する場合、スイッチSW2のみが通電され、スイッチSW1,SW3,SW4は切断される。オペアンプ170の-入力端子に1.4Vの電圧VL1を印加する場合、スイッチSW3のみが通電され、スイッチSW1,SW2,SW4は切断される。オペアンプ170の-入力端子に0Vの電圧VL0を印加する場合、スイッチSW4のみが通電され、スイッチSW1,SW2,SW3は切断される。したがって、スイッチSW1,SW2,SW3,SW4の開閉によって、オペアンプ170の-入力端子に0V又は3段階の電圧のいずれかを印加可能である。そのため、スイッチSW1,SW2,SW3,SW4の開閉によって、発熱素子61の発熱温度を定める印加電圧を3段階に設定可能である。ここで、オペアンプ170の入力端子に1.4Vの電圧VL1を印加した場合の発熱素子61の温度をTH1とする。また、オペアンプ170の入力端子に1.9Vの電圧VL2を印加した場合の発熱素子61の温度をTH2、オペアンプ170の入力端子に2.4Vの電圧VL3を印加した場合の発熱素子61の温度をTH3とする。図6に示すように、第1の測温素子62の一端には、例えば、オペアンプ270の-入力端子が電気的に接続され、他端は接地される。また、オペアンプ270の-入力端子及び出力端子と並列に、抵抗素子261が接続される。オペアンプ270の入力端子は、直列に接続された抵抗素子264と抵抗素子265との間に電気的に接続される。これにより、第1の測温素子62には、0.3V程度の弱い電圧が加えられる。図1及び図2に示す発熱素子61の抵抗値は、発熱素子61の温度によって変化する。発熱素子61の温度THと、発熱素子61の抵抗値RHの関係は、下記1式で与えられる。RHRHSTD×1αHβH2ここで、THSTDは発熱素子61の標準温度を表し、例えば20°Cである。RHSTDは標準温度THSTDにおける予め測定された発熱素子61の抵抗値を表す。αHは1次の抵抗温度係数を表す。βHは2次の抵抗温度係数を表す。発熱素子61の抵抗値RHは、発熱素子61の駆動電力PHと、発熱素子61の通電電流IHから、下記2式で与えられる。RHPHIH2あるいは発熱素子61の抵抗値RHは、発熱素子61にかかる電圧VHと、発熱素子61の通電電流IHから、下記3式で与えられる。RHVHIHここで、発熱素子61の温度THは、発熱素子61と雰囲気ガスの間が熱的に平衡になったときに安定する。なお、熱的に平衡な状態とは、発熱素子61の発熱と、発熱素子61から雰囲気ガスへの放熱とが釣り合っている状態をいう。下記4式に示すように、平衡状態における発熱素子61の駆動電力PHを、発熱素子61の温度THと雰囲気ガスの温度TIとの差ΔTHで割ることにより、雰囲気ガスの放熱係数MIが得られる。なお、放熱係数MIの単位は、例えばW/°Cである。MIPHPHΔTH上記1式より、発熱素子61の温度THは下記5式で与えられる。TH×αHαH24βH12THSTDしたがって、発熱素子61の温度THと雰囲気ガスの温度TIとの差ΔTHは、下記6式で与えられる。ΔTH×αHαH24βH12THSTDTI雰囲気ガスの温度TIは、自己発熱しない程度の電力を与えられる第1の測温素子62の温度TIに近似する。第1の測温素子62の温度TIと、第1の測温素子62の抵抗値RIの関係は、下記7式で与えられる。RIRISTD×1αIβI2TISTDは第1の測温素子62の標準温度を表し、例えば20°Cである。RISTDは標準温度TISTDにおける予め測定された第1の測温素子62の抵抗値を表す。αIは1次の抵抗温度係数を表す。βIは2次の抵抗温度係数を表す。上記7式より、第1の測温素子62の温度TIは下記8式で与えられる。TI×αIαI24βI12TISTDよって、雰囲気ガスの放熱係数MIは、下記9式で与えられる。MIPHΔTHPHαHαH24βH12THSTDαIαI24βI12TISTD発熱素子61の通電電流IHと、駆動電力PH又は電圧VHは測定可能であるため、上記2式又は3式から発熱素子61の抵抗値RHを算出可能である。同様に、第1の測温素子62の抵抗値RIも算出可能である。よって、マイクロチップ8を用いて、上記9式から雰囲気ガスの放熱係数MIが算出可能である。なお、保温素子64で基板60の温度を一定に保つことにより、発熱素子61が発熱する前のマイクロチップ8の近傍の雰囲気ガスの温度が、基板60の一定の温度と近似する。そのため、発熱素子61が発熱する前の雰囲気ガスの温度の変動が抑制される。温度変動が一度抑制された雰囲気ガスを発熱素子61でさらに加熱することにより、より高い精度で放熱係数MIを算出することが可能となる。ここで、雰囲気ガスが混合ガスであり、混合ガスが、ガスA、ガスB、ガスC、及びガスDの4種類のガス成分からなっているとする。ガスAの体積率VA、ガスBの体積率VB、ガスCの体積率VC、及びガスDの体積率VDの総和は、下記10式で与えられるように、1である。VAVBVCVD1ガスAの単位体積当たりの発熱量をKA、ガスBの単位体積当たりの発熱量をKB、ガスCの単位体積当たりの発熱量をKC、ガスDの単位体積当たりの発熱量をKDとすると、混合ガスの単位体積当たりの発熱量Qは、各ガス成分の体積率に、各ガス成分の単位体積当たりの発熱量を乗じたものの総和で与えられる。したがって、混合ガスの単位体積当たりの発熱量Qは、下記11式で与えられる。なお、単位体積当たりの発熱量の単位は、例えばMJ/m3である。QKA×VAKB×VBKC×VCKD×VDまた、ガスAを構成する水素原子による単位体積当たりの発熱量をKAH、ガスBを構成する水素原子による単位体積当たりの発熱量をKBH、ガスCを構成する水素原子による単位体積当たりの発熱量をKCH、ガスDを構成する水素原子による単位体積当たりの発熱量をKDHとすると、混合ガスに含まれる分子を構成する水素原子による単位体積当たりの発熱量QHは、各ガス成分の体積率に、各ガス成分を構成する水素原子による単位体積当たりの発熱量を乗じたものの総和で与えられる。したがって、混合ガスに含まれる分子を構成する水素原子による単位体積当たりの発熱量QHは、下記12式で与えられる。QHKAH×VAKBH×VBKCH×VCKDH×VD次に、ガスAの放熱係数をMA、ガスBの放熱係数をMB、ガスCの放熱係数をMC、ガスDの放熱係数をMDとすると、混合ガスの放熱係数MIは、各ガス成分の体積率に、各ガス成分の放熱係数を乗じたものの総和で与えられる。したがって、混合ガスの放熱係数MIは、下記13式で与えられる。MIMA×VAMB×VBMC×VCMD×VDさらに、ガスの放熱係数は発熱素子61の温度THに依存するので、混合ガスの放熱係数MIは、発熱素子61の温度THの関数として、下記14式で与えられる。MIMA×VAMB×VBMC×VCMD×VDしたがって、発熱素子61の温度がTH1のときの混合ガスの放熱係数MI1TH1は下記15式で与えられる。また、発熱素子61の温度がTH2のときの混合ガスの放熱係数MI2TH2は下記16式で与えられ、発熱素子61の温度がTH3のときの混合ガスの放熱係数MI3TH3は下記17式で与えられる。MI1MA×VAMB×VBMC×VCMD×VDMI2MA×VAMB×VBMC×VCMD×VDMI3MA×VAMB×VBMC×VCMD×VDここで、発熱素子61の温度THに対して各ガス成分の放熱係数MATH,MBTH,MCTH,MDTHが非線形性を有する場合、上記15乃至17式は、線形独立な関係を有する。また、発熱素子61の温度THに対して各ガス成分の放熱係数MATH,MBTH,MCTH,MDTHが線形性を有する場合でも、発熱素子61の温度THに対する各ガス成分の放熱係数MATH,MBTH,MCTH,MDTHの変化率が異なる場合は、上記15乃至17式は、線形独立な関係を有する。さらに、15乃至17式が線形独立な関係を有する場合、10及び15乃至17式は線形独立な関係を有する。図7は、天然ガスや都市ガスに含まれるメタンCH4、プロパンC3H8、窒素N2、及び二酸化炭素CO2の放熱係数と、発熱抵抗体である発熱素子61の温度との関係を示すグラフである。発熱素子61の温度に対して、メタンCH4、プロパンC3H8、窒素N2、及び二酸化炭素CO2のそれぞれのガス成分の放熱係数は線形性を有する。しかし、発熱素子61の温度に対する放熱係数の変化率は、メタンCH4、プロパンC3H8、窒素N2、及び二酸化炭素CO2のそれぞれで異なる。したがって、混合ガスを構成するガス成分がメタンCH4、プロパンC3H8、窒素N2、及び二酸化炭素CO2であるである場合、上記15乃至17式は、線形独立な関係を有する。15乃至17式中の各ガス成分の放熱係数MATH1,MBTH1,MCTH1,MDTH1,MATH2,MBTH2,MCTH2,MDTH2,MATH3,MBTH3,MCTH3,MDTH3の値は、測定等により予め得ることが可能である。したがって、10及び15乃至17式の連立方程式を解くと、ガスAの体積率VA、ガスBの体積率VB、ガスCの体積率VC、及びガスDの体積率VDのそれぞれが、下記18乃至21式に示すように、混合ガスの放熱係数MI1TH1,MI2TH2,MI3TH3の関数として与えられる。なお、下記18乃至21式において、nを自然数として、fnは関数を表す記号である。VAf1MI1MI2MI3VBf2MI1MI2MI3VCf3MI1MI2MI3VDf4MI1MI2MI3ここで、上記11式に18乃至21式を代入することにより、下記22式が得られる。QKA×VAKB×VBKC×VCKD×VDKA×f1MI1MI2MI3KB×f2MI1MI2MI3KC×f3MI1MI2MI3KD×f4MI1MI2MI3上記22式に示すように、混合ガスの単位体積当たりの発熱量Qは、発熱素子61の温度がTH1,TH2,TH3である場合の混合ガスの放熱係数MI1TH1,MI2TH2,MI3TH3を変数とする方程式で与えられる。したがって、混合ガスの発熱量Qは、g1を関数を表す記号として、下記23式で与えられる。Qg1MI1MI2MI3また、上記12式、及び18乃至21式より、混合ガスに含まれる分子を構成する水素原子による単位体積当たりの発熱量QHは、g2を関数を表す記号として、下記24式で与えられる。QHg2MI1MI2MI3よって、ガスA、ガスB、ガスC、及びガスDからなる混合ガスについて、予め上記24式を得れば、ガスAの体積率VA、ガスBの体積率VB、ガスCの体積率VC、及びガスDの体積率VDが未知の測定対象混合ガスに含まれる分子を構成する水素原子による単位体積当たりの発熱量QHを容易に算出可能であることを、発明者は見出した。具体的には、発熱素子61の発熱温度がTH1,TH2,TH3である場合の測定対象混合ガスの放熱係数MI1TH1,MI2TH2,MI3TH3を測定し、24式に代入することにより、測定対象混合ガスに含まれる分子を構成する水素原子による発熱量QHを一意に求めることが可能となる。次に、混合ガスの放熱係数MIは、上記9式に示すように、発熱素子61の抵抗値RHと、第1の測温素子62の抵抗値RIと、に依存する。そこで、本発明者は、混合ガスに含まれる分子を構成する水素原子による単位体積当たりの発熱量QHは、h1を関数を表す記号として、下記25式に示すように、発熱素子61の温度がTH1,TH2,TH3である場合の発熱素子61の抵抗値RH1TH1,RH2TH2,RH3TH3と、混合ガスに接する第1の測温素子62の抵抗値RIと、を変数とする方程式でも与えられることを見出した。QHh1RH1RH2RH3RIよって、測定対象混合ガスに接する発熱素子61の発熱温度がTH1,TH2,TH3である場合の発熱素子61の抵抗値RH1TH1,RH2TH2,RH3TH3と、測定対象混合ガスに接する第1の測温素子62の抵抗値RIを測定し、25式に代入することによっても、測定対象混合ガスに含まれる分子を構成する水素原子による発熱量QHを一意に求めることが可能となる。また、混合ガスに含まれる分子を構成する水素原子による単位体積当たりの発熱量QHは、h2を関数を表す記号として、下記26式に示すように、発熱素子61の温度がTH1,TH2,TH3である場合の発熱素子61の通電電流IH1TH1,IH2TH2,IH3TH3と、混合ガスに接する第1の測温素子62の通電電流IIと、を変数とする方程式でも与えられる。QHh2IH1IH2IH3IIあるいは混合ガスに含まれる分子を構成する水素原子による単位体積当たりの発熱量QHは、h3を関数を表す記号として、下記27式に示すように、発熱素子61の温度がTH1,TH2,TH3である場合の発熱素子61にかかる電圧VH1TH1,VH2TH2,VH3TH3と、混合ガスに接する第1の測温素子62にかかる電圧VIと、を変数とする方程式でも与えられる。QHh3VH1VH2VH3VIまたあるいは混合ガスに含まれる分子を構成する水素原子による単位体積当たりの発熱量QHは、h4を関数を表す記号として、下記28式に示すように、発熱素子61の温度がTH1,TH2,TH3である場合の発熱素子61に接続されたアナログ-デジタル変換回路以下においてA/D変換回路という。の出力信号ADH1TH1,ADH2TH2,ADH3TH3と、混合ガスに接する第1の測温素子62に接続されたA/D変換回路の出力信号ADIと、を変数とする方程式でも与えられる。QHh4ADH1ADH2ADH3ADIよって、混合ガスに含まれる分子を構成する水素原子による単位体積当たりの発熱量QHは、h5を関数を表す記号として、下記29式に示すように、発熱素子61の発熱温度がTH1,TH2,TH3である場合の発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3と、混合ガスに接する第1の測温素子62からの電気信号SIと、を変数とする方程式で与えられる。QHh5SH1SH2SH3SI次に、混合ガスに含まれる分子を構成する水素原子による単位体積当たりの発熱量QHは、混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量例えば、1.00794と、それぞれの分子の混合ガスにおける体積比と、の積の総和と相関する。したがって、上記29式に基づき、混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和GHは、h6を関数を表す記号として、下記30式で与えられる。GHh6SH1SH2SH3SIなお、混合ガスのガス成分は、4種類に限定されることはない。例えば、混合ガスがn種類のガス成分からなる場合、まず、下記31式で与えられる、少なくともn-1種類の発熱温度TH1,TH2,TH3,,THn1における発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3,,SHn1THn1と、混合ガスに接する第1の測温素子62からの電気信号SIと、を変数とする方程式を予め取得する。そして、n-1種類の発熱温度TH1,TH2,TH3,,THn1における、n種類のガス成分のそれぞれの体積率が未知の測定対象混合ガスに接する発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3,,SHn1THn1の値と、測定対象混合ガスに接する第1の測温素子62からの電気信号SIの値と、を測定し、31式に代入することにより、測定対象混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和GHを一意に求めることが可能となる。GHh12SH1SH2SH3SHn1SI例えば、混合ガスがメタンCH4を90体積、エタンC2H6を5体積、プロパンC3H8を1体積、ブタンC4H10を1体積、窒素N2を1体積、二酸化炭素CO2を2含むとする。この場合、メタンを構成する水素原子の数は4であるから、メタンを構成する水素原子の数と、水素原子の原子量と、メタンの体積比と、の積は、4×1.00794×0.9=3.628584である。また、エタンを構成する水素原子の数は6であるから、エタンを構成する水素原子の数と、水素原子の原子量と、エタンの体積比と、の積は、6×1.00794×0.05=0.302382である。また、プロパンを構成する水素原子の数は8であるから、プロパンを構成する水素原子の数と、水素原子の原子量と、プロパンの体積比と、の積は、8×1.00794×0.01=0.0806352である。また、ブタンを構成する水素原子の数は10であるから、ブタンを構成する水素原子の数と、水素原子の原子量と、ブタンの体積比と、の積は、10×1.00794×0.01=0.100794である。窒素及び二酸化炭素を構成する水素原子の数は0である。したがって、メタンを構成する水素原子の数と水素原子の原子量とメタンの体積比との積と、エタンを構成する水素原子の数と水素原子の原子量とエタンの体積比との積と、プロパンを構成する水素原子の数と水素原子の原子量とプロパンの体積比との積と、ブタンを構成する水素原子の数と水素原子の原子量とブタンの体積比との積と、の総和は、3.6285840.3023820.08063520.100794=4.1123952である。ただし、混合ガスが、ガス成分としてメタンCH4、プロパンC3H8に加えて、jを自然数として、メタンCH4とプロパンC3H8以外のアルカンCjH2j2を含む場合、メタンCH4とプロパンC3H8以外のアルカンCjH2j2を、メタンCH4とプロパンC3H8の混合物とみなしても、31式の算出には影響しない。例えば、エタンC2H6、ブタンC4H10、ペンタンC5H12、ヘキサンC6H14を、下記32乃至35式に示すように、それぞれ所定の係数を掛けられたメタンCH4とプロパンC3H8の混合物とみなして31式を算出してもかまわない。C2H605CH405C3H8C4H1005CH415C3H8C5H1210CH420C3H8C6H1415CH425C3H8したがって、zを自然数として、n種類のガス成分からなる混合ガスが、ガス成分としてメタンCH4、プロパンC3H8に加えて、メタンCH4とプロパンC3H8以外のz種類のアルカンCjH2j2を含む場合、少なくともn-z-1種類の発熱温度における発熱素子61からの電気信号SHと、第1の測温素子62からの電気信号SIと、を変数とする方程式を求めてもよい。31式の算出に用いられた混合ガスのガス成分の種類と、測定対象混合ガスのガス成分の種類が同じ場合に、31式を利用可能であることはもちろんである。さらに、測定対象混合ガスがn種類より少ない種類のガス成分からなり、かつ、n種類より少ない種類のガス成分が、31式の算出に用いられた混合ガスに含まれている場合も、31式を利用可能である。例えば、31式の算出に用いられた混合ガスが、メタンCH4、プロパンC3H8、窒素N2、及び二酸化炭素CO2の4種類のガス成分を含む場合、測定対象混合ガスが、窒素N2を含まず、メタンCH4、プロパンC3H8、及び二酸化炭素CO2の3種類のガス成分のみを含む場合も、31式を利用可能である。さらに、31式の算出に用いられた混合ガスが、ガス成分としてメタンCH4とプロパンC3H8を含む場合、測定対象混合ガスが、31式の算出に用いられた混合ガスに含まれていないアルカンCjH2j2を含んでいても、31式を利用可能である。これは、上述したように、メタンCH4とプロパンC3H8以外のアルカンCjH2j2を、メタンCH4とプロパンC3H8の混合物とみなし得るためである。ここで、図8及び図9に示す第1の実施の形態に係るガス測定システム20は、複数のサンプル混合ガスのそれぞれが流れるパイプ101と、パイプ101に配置され、図1に示す第1の測温素子62及び複数の発熱温度THで発熱する発熱素子61を含むマイクロチップ8と、を備える。なお、それぞれのサンプル混合ガスは、複数種類のガス成分を含む。また、図8に示すガス測定システム20は、複数のサンプル混合ガスのそれぞれの温度TIに依存する第1の測温素子62からの電気信号SIの値と、複数の発熱温度THのそれぞれにおける発熱素子61からの電気信号SHの値と、を測定する測定部301を備える。さらに、ガス測定システム20は、複数の混合ガスのそれぞれに含まれる分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の既知の値と、第1の測温素子62からの電気信号SIの値と、複数の発熱温度における発熱素子61からの電気信号の値と、に基づいて、第1の測温素子62からの電気信号SI及び複数の発熱温度THにおける発熱素子61からの電気信号SHを独立変数とし、ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和を従属変数とする式を作成する式作成部302を備える。それぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の値がそれぞれ異なる4種類のサンプル混合ガスが使用される場合、図9に示すように、第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及びパイプ101を介してガス測定システム20に送られる第1のサンプル混合ガスの流量を制御する。第2のガスボンベ50Bには、パイプ91Bを介して、第2のガス圧調節器31Bが接続されている。また、第2のガス圧調節器31Bには、パイプ92Bを介して、第2の流量制御装置32Bが接続されている。第2の流量制御装置32Bは、パイプ92B,93,101を介してガス測定システム20に送られる第2のサンプル混合ガスの流量を制御する。第3のガスボンベ50Cには、パイプ91Cを介して、第3のガス圧調節器31Cが接続されている。また、第3のガス圧調節器31Cには、パイプ92Cを介して、第3の流量制御装置32Cが接続されている。第3の流量制御装置32Cは、パイプ92C,93,101を介してガス測定システム20に送られる第3のサンプル混合ガスの流量を制御する。第4のガスボンベ50Dには、パイプ91Dを介して、第4のガス圧調節器31Dが接続されている。また、第4のガス圧調節器31Dには、パイプ92Dを介して、第4の流量制御装置32Dが接続されている。第4の流量制御装置32Dは、パイプ92D,93,101を介してガス測定システム20に送られる第4のサンプル混合ガスの流量を制御する。第1乃至第4のサンプル混合ガスのそれぞれは、例えば天然ガス又は都市ガスである。第1乃至第4のサンプル混合ガスは、例えばメタンCH4、プロパンC3H8、窒素N2、及び二酸化炭素CO2の4種類のガス成分を、それぞれ異なる体積比で含む。図8に示すパイプ101に第1のサンプル混合ガスが供給された際、図1及び図2に示すマイクロチップ8の第1の測温素子62は、第1のサンプル混合ガスの温度に依存する電気信号SIを出力する。次に、発熱素子61は、図8に示す駆動回路303から駆動電力PH1,PH2,PH3を順次与えられる。駆動電力PH1,PH2,PH3を与えられた場合、第1のサンプル混合ガスに接する発熱素子61は、例えば、100°Cの温度TH1、150°Cの温度TH2、及び200°Cの温度TH3で順次発熱し、発熱温度TH1における電気信号SH1TH1、発熱温度TH2における電気信号SH2TH2、及び発熱温度TH3における電気信号SH3TH3を出力する。パイプ101から第1のサンプル混合ガスが除去された後、第2乃至第4のサンプル混合ガスがパイプ101に順次供給される。パイプ101に第2のサンプル混合ガスが供給された際、図1及び図2に示すマイクロチップ8の第1の測温素子62は、第2のサンプル混合ガスの温度に依存する電気信号SIを出力する。次に、第2のサンプル混合ガスに接する発熱素子61は、発熱温度TH1における電気信号SH1TH1、発熱温度TH2における電気信号SH2TH2、及び発熱温度TH3における電気信号SH3TH3を出力する。図8に示すパイプ101に第3のサンプル混合ガスが供給された際、図1及び図2に示すマイクロチップ8の第1の測温素子62は、第3のサンプル混合ガスの温度に依存する電気信号SIを出力する。次に、第3のサンプル混合ガスに接する発熱素子61は、発熱温度TH1における電気信号SH1TH1、発熱温度TH2における電気信号SH2TH2、及び発熱温度TH3における電気信号SH3TH3を出力する。図8に示すパイプ101に第4のサンプル混合ガスが供給された際、図1及び図2に示すマイクロチップ8の第1の測温素子62は、第4のサンプル混合ガスの温度に依存する電気信号SIを出力する。次に、第4のサンプル混合ガスに接する発熱素子61は、発熱温度TH1における電気信号SH1TH1、発熱温度TH2における電気信号SH2TH2、及び発熱温度TH3における電気信号SH3TH3を出力する。なお、それぞれのサンプル混合ガスがn種類のガス成分を含む場合、マイクロチップ8の図1及び図2に示す発熱素子61は、少なくともn-1種類の異なる温度で発熱させられる。ただし、上述したように、メタンCH4及びプロパンC3H8以外のアルカンCjH2j2は、メタンCH4及びプロパンC3H8の混合物とみなしうる。したがって、zを自然数として、n種類のガス成分からなるサンプル混合ガスが、ガス成分としてメタンCH4及びプロパンC3H8に加えてz種類のアルカンCjH2j2を含む場合は、発熱素子61は、少なくともn-z-1種類の異なる温度で発熱させられる。図8に示すように、マイクロチップ8は、測定部301を含む中央演算処理装置CPU300に接続されている。CPU300には、電気信号記憶装置401が接続されている。測定部301は、第1の測温素子62からの電気信号SIの値と、発熱素子61からの発熱温度TH1における電気信号SH1TH1、発熱温度TH2における電気信号SH2TH2、及び発熱温度TH3における電気信号SH3TH3の値と、を測定し、測定値を電気信号記憶装置401に保存する。なお、第1の測温素子62からの電気信号SIとは、第1の測温素子62の抵抗値RI、第1の測温素子62の通電電流II、第1の測温素子62にかかる電圧VI、及び第1の測温素子62に接続されたA/D変換回路304の出力信号ADIのいずれであってもよい。同様に、発熱素子61からの電気信号SHとは、発熱素子61の抵抗値RH、発熱素子61の通電電流IH、発熱素子61にかかる電圧VH、及び発熱素子61に接続されたA/D変換回路304の出力信号ADHのいずれであってもよい。CPU300に含まれる式作成部302は、例えば第1のサンプル混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の既知の値と、第2のサンプル混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の既知の値と、第3のサンプル混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の既知の値と、第4のサンプル混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の既知の値と、第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の複数の測定値と、を収集する。さらに式作成部302は、収集した積の総和の値と、電気信号SIの値と、電気信号SHの値と、に基づいて、多変量解析により、第1の測温素子62からの電気信号SI及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3を独立変数とし、ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和を従属変数とする式を算出する。なお、多変量解析とは、A.JSmola及びB.Scholkopf著のATutorialonSupportVectorRegressionNeuroCOLTTechnicalReportNC-TR-98-030、1998年に開示されているサポートベクトル回帰、重回帰分析、及び特開平5-141999号公報に開示されているファジィ数量化理論II類等を含む。ガス測定システム20は、CPU300に接続された式記憶装置402をさらに備える。式記憶装置402は、式作成部302が作成した式を保存する。さらにCPU300には、入力装置312及び出力装置313が接続される。入力装置312としては、例えばキーボード、及びマウス等のポインティングデバイス等が使用可能である。出力装置313には液晶ディスプレイ、モニタ等の画像表示装置、及びプリンタ等が使用可能である。次に、図10に示すフローチャートを用いて第1の実施の形態に係る、ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和を算出するための式の作成方法について説明する。aステップS100で、図9に示す第2乃至第4の流量制御装置32B-32Dの弁を閉じたまま、第1の流量制御装置32Aの弁を開き、図8に示すパイプ101内に第1のサンプル混合ガスを導入する。ステップS101で、測定部301は、第1のサンプル混合ガスに接する第1の測温素子62からの電気信号SIの値を測定し、電気信号記憶装置401に保存する。次に、駆動回路303は、図1及び図2に示す発熱素子61に駆動電力PH1を与え、発熱素子61を100°Cで発熱させる。図8に示す測定部301は、100°Cで発熱する発熱素子61からの電気信号SH1TH1の値を、電気信号記憶装置401に保存する。bステップS102で、駆動回路303は、図1及び図2に示す発熱素子61の温度の切り替えが完了したか否か判定する。温度150°C及び温度200°Cへの切り替えが完了していない場合には、ステップS101に戻り、図8に示す駆動回路303は、図1及び図2に示す発熱素子61を150°Cで発熱させる。図8に示す測定部301は、第1のサンプル混合ガスに接し、150°Cで発熱する発熱素子61からの電気信号SH2TH2の値を、電気信号記憶装置401に保存する。c再びステップS102で、図1及び図2に示す発熱素子61の温度の切り替えが完了したか否か判定する。温度200°Cへの切り替えが完了していない場合には、ステップS101に戻り、図8に示す駆動回路303は、図1及び図2に示す発熱素子61を200°Cで発熱させる。図8に示す測定部301は、第1のサンプル混合ガスに接し、200°Cで発熱する発熱素子61からの電気信号SH3TH3の値を、電気信号記憶装置401に保存する。d発熱素子61の温度の切り替えが完了した場合には、ステップS102からステップS103に進む。ステップS103で、サンプル混合ガスの切り替えが完了したか否かを判定する。第2乃至第4のサンプル混合ガスへの切り替えが完了していない場合には、ステップS100に戻る。ステップS100で、図9に示す第1の流量制御装置32Aを閉じ、第3乃至第4の流量制御装置32C-32Dの弁を閉じたまま第2の流量制御装置32Bの弁を開き、図8に示すパイプ101内に第2のサンプル混合ガスを導入する。e第1のサンプル混合ガスと同様に、ステップS101乃至ステップS102のループが繰り返される。測定部301は、第2のサンプル混合ガスに接する第1の測温素子62からの電気信号SIの値を測定し、電気信号記憶装置401に保存する。また、測定部301は、第2のサンプル混合ガスに接し、100°C、150°C、及び200°Cで発熱する発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の値を、電気信号記憶装置401に保存する。fその後、ステップS100乃至ステップS103のループが繰り返される。これにより、パイプ101内に供給された第3のサンプル混合ガスに接する第1の測温素子62からの電気信号SIの値と、第3のサンプル混合ガスに接し、100°C、150°C、及び200°Cで発熱する発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の値と、が電気信号記憶装置401に保存される。また、パイプ101内に供給された第4のサンプル混合ガスに接する第1の測温素子62からの電気信号SIの値と、第4のサンプル混合ガスに接し、100°C、150°C、及び200°Cで発熱する発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の値と、が電気信号記憶装置401に保存される。gステップS104で、入力装置312から式作成部302に、第1のサンプル混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の既知の値と、第2のサンプル混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の既知の値と、第3のサンプル混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の既知の値と、第4のサンプル混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の既知の値と、を入力する。また、式作成部302は、電気信号記憶装置401から、第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の複数の測定値と、を読み出す。hステップS105で、積の総和の既知の値と、第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の複数の測定値と、に基づいて、式作成部302は、重回帰分析を行う。重回帰分析により、式作成部302は、第1の測温素子62からの電気信号SI及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3を独立変数とし、ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和を従属変数とする式を算出する。その後、ステップS106で、式作成部302は作成した式を式記憶装置402に保存し、第1の実施の形態に係る式の作成方法が終了する。以上示したように、本発明の第1の実施の形態によれば、測定対象混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の値を一意に算出可能な式を作成することが可能となる。次に、測定対象混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の値を測定する際の、図8に示す第1の実施の形態に係るガス測定システム20の機能を説明する。例えば未知の体積率でメタンCH4、プロパンC3H8、窒素N2、及び二酸化炭素CO2等を含む天然ガス又は都市ガス等の測定対象混合ガスが、パイプ101に導入される。図1及び図2に示すマイクロチップ8の第1の測温素子62は、測定対象混合ガスの温度に依存する電気信号SIを出力する。次に、発熱素子61は、図8に示す駆動回路303から駆動電力PH1,PH2,PH3を与えられる。駆動電力PH1,PH2,PH3を与えられた場合、測定対象混合ガスに接する発熱素子61は、例えば、100°Cの温度TH1、150°Cの温度TH2、及び200°Cの温度TH3で発熱し、発熱温度TH1における電気信号SH1TH1、発熱温度TH2における電気信号SH2TH2、及び発熱温度TH3における電気信号SH3TH3を出力する。図8に示す測定部301は、パイプ101内に供給された測定対象混合ガスに接し、測定対象混合ガスの温度TIに依存する第1の測温素子62からの電気信号SIの値と、測定対象混合ガスに接する発熱素子61からの発熱温度TH1における電気信号SH1TH1、発熱温度TH2における電気信号SH2TH2、及び発熱温度TH3における電気信号SH3TH3の値と、を測定し、測定値を電気信号記憶装置401に保存する。上述したように、式記憶装置402は、第1の測温素子62からの電気信号SIと、発熱温度TH1が100°Cの発熱素子61からの電気信号SH1TH1と、発熱温度TH2が150°Cの発熱素子61からの電気信号SH2TH2と、発熱温度TH3が200°Cの発熱素子61からの電気信号SH3TH3と、を独立変数とし、ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和を従属変数とする式を保存している。第1の実施の形態に係るガス測定システム20は、さらに、積算出部305を備える。積算出部305は、式記憶装置402に保存されている式の第1の測温素子62からの電気信号SIの独立変数、及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の独立変数に、第1の測温素子62からの電気信号SIの測定値、及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の測定値をそれぞれ代入し、測定対象混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の値を算出する。CPU300には、積記憶装置403がさらに接続されている。積記憶装置403は、積算出部305が算出した積の総和の値を保存する。次に、図11に示すフローチャートを用いて、第1の実施の形態に係るガス測定方法について説明する。aステップS200で、図8に示すパイプ101内に測定対象混合ガスを導入する。ステップS201で、測定部301は、測定対象混合ガスに接する第1の測温素子62からの電気信号SIの値を測定し、電気信号記憶装置401に保存する。次に、駆動回路303は、図1及び図2に示す発熱素子61に駆動電力PH1を与え、発熱素子61を100°Cで発熱させる。図8に示す測定部301は、測定対象混合ガスに接し、100°Cで発熱する発熱素子61からの電気信号SH1TH1の値を、電気信号記憶装置401に保存する。bステップS202で、図8に示す駆動回路303は、図1及び図2に示す発熱素子61の温度の切り替えが完了したか否か判定する。温度150°C及び温度200°Cへの切り替えが完了していない場合には、ステップS201に戻り、駆動回路303は、図1及び図2に示す発熱素子61に駆動電力PH2を与え、発熱素子61を150°Cで発熱させる。図8に示す測定部301は、測定対象混合ガスに接し、150°Cで発熱する発熱素子61からの電気信号SH2TH2の値を、電気信号記憶装置401に保存する。c再びステップS202で、図1及び図2に示す発熱素子61の温度の切り替えが完了したか否か判定する。温度200°Cへの切り替えが完了していない場合には、ステップS201に戻り、駆動回路303は、図1及び図2に示す発熱素子61に駆動電力PH3を与え、発熱素子61を200°Cで発熱させる。図8に示す測定部301は、測定対象混合ガスに接し、200°Cで発熱する発熱素子61からの電気信号SH3TH3の値を、電気信号記憶装置401に保存する。d発熱素子61の温度の切り替えが完了した場合には、ステップS202からステップS203に進む。ステップS203で、図8に示す積算出部305は、式記憶装置402から、第1の測温素子62からの電気信号SI及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3を独立変数とし、ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和を従属変数とする式を読み出す。また、積算出部305は、電気信号記憶装置401から、測定対象混合ガスに接する第1の測温素子62からの電気信号SIの測定値、及び測定対象混合ガスに接する発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の測定値を読み出す。eステップS204で、積算出部305は、式の電気信号SI及び電気信号SH1TH1,SH2TH2,SH3TH3の独立変数に、それぞれ測定値を代入し、測定対象混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の値を算出する。その後、積算出部305は、算出した積の総和の値を積記憶装置403に保存し、第1の実施の形態に係るガス測定方法を終了する。以上説明した本発明の第1の実施の形態によれば、測定対象混合ガスに接する第1の測温素子62からの電気信号SIの値と、測定対象混合ガスに接する発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の値と、から、測定対象混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の値を測定することが可能となる。第1の実施の形態の第1の実施例まず、組成が既知の40種類のサンプル混合ガスを用意した。40種類のサンプル混合ガスのそれぞれは、ガス成分としてメタン、エタン、プロパン、ブタン、ペンタン、ヘキサン、窒素、及び二酸化炭素のいずれか又は全部を含んでいた。次に、40種類のサンプル混合ガスのそれぞれを用いて、図1に示す第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3,SH4TH4,SH5TH5の複数の測定値と、を取得した。その後、40種類のサンプル混合ガスのそれぞれに含まれるそれぞれの分子を構成する水素原子の数と、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の既知の値と、第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3,SH4TH4,SH5TH5の複数の測定値と、に基づいて、サポートベクトル回帰により、第1の測温素子62からの電気信号SI及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3,SH4TH4,SH5TH5を独立変数とし、ガスに含まれるそれぞれの分子を構成する水素原子の数と、水素原子の原子量と、それぞれの分子の体積比と、の積の総和を従属変数とする方程式を作成した。方程式を作成する際には、キャリブレーションポイントは、3乃至5個を目安に、適宜決定できる。作成された方程式を用いて、40種類のサンプル混合ガスのそれぞれに含まれるそれぞれの分子を構成する水素原子の数と、水素原子の原子量と、それぞれの分子の体積比と、の積の総和を算出し、真の値と比較したところ、図12に示すように、誤差は±1.25であった。第1の実施の形態の第2の実施例第1の実施例と同様に、組成が既知の40種類のサンプル混合ガスを用意した。次に、40種類のサンプル混合ガスのそれぞれを用いて、図1に示す第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3,SH4TH4,SH5TH5の複数の測定値と、を取得した。その後、上記32式及び33式に示すように、エタン及びブタンのそれぞれを、メタンとプロパンの混合物とみなして、組成が既知の40種類のサンプル混合ガスのそれぞれに含まれるそれぞれの分子を構成する水素原子の数と、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の値を算出した。例えば、サンプル混合ガスの一つがメタンを90体積、エタンを5体積、プロパンを1体積、ブタンを1体積、窒素を1体積、二酸化炭素を2含むとする。この場合、エタン及びブタンを、メタン及びプロパンに分解すると、メタンの体積は、元のメタンの90体積と、エタンから分解されたメタンの5体積×0.5と、ブタンから分解されたメタンの1体積×-0.5と、の和である92体積とみなしうる。また、プロパンの体積は、エタンから分解されたプロパンの5体積×0.5と、元のプロパンの1体積と、ブタンから分解されたプロパンの1体積×1.5と、の和の5体積とみなしうる。さらに、算出した積の総和の値と、第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3,SH4TH4,SH5TH5の複数の測定値と、に基づいて、サポートベクトル回帰により、第1の測温素子62からの電気信号SI及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3,SH4TH4,SH5TH5を独立変数とし、ガスに含まれるそれぞれの分子を構成する水素原子の数と、水素原子の原子量と、それぞれの分子の体積比と、の積の総和を従属変数とする方程式を作成した。方程式を作成する際には、キャリブレーションポイントは、3乃至5個を目安に、適宜決定できる。作成された方程式を用いて、40種類のサンプル混合ガスのそれぞれに含まれるそれぞれの分子を構成する水素原子の数と、水素原子の原子量と、それぞれの分子の体積比と、の積の総和を算出し、真の値と比較したところ、図13に示すように、誤差は±2.0であった。第2の実施の形態図1乃至図4に示すマイクロチップ8の雰囲気ガスが混合ガスであり、混合ガスが、ガスA、ガスB、ガスC、及びガスDの4種類のガス成分からなっているとする。ガスAに含まれる炭素原子による単位体積当たりの発熱量をKAC、ガスBに含まれる炭素原子による単位体積当たりの発熱量をKBC、ガスCに含まれる炭素原子による単位体積当たりの発熱量をKCC、ガスDに含まれる炭素原子による単位体積当たりの発熱量をKDCとすると、混合ガスに含まれる炭素原子による単位体積当たりの発熱量QCは、各ガス成分の体積率に、各ガス成分に含まれる炭素原子による単位体積当たりの発熱量を乗じたものの総和で与えられる。したがって、混合ガスに含まれる分子を構成する炭素原子による単位体積当たりの発熱量QCは、下記36式で与えられる。QCKAC×VAKBC×VBKCC×VCKDC×VD上記36式、及び18乃至21式より、混合ガスに含まれる分子を構成する炭素原子による単位体積当たりの発熱量QCは、g3を関数を表す記号として、下記37式で与えられる。QCg3MI1MI2MI3また、上述したように、混合ガスの放熱係数MIは、発熱素子61の発熱温度がTH1,TH2,TH3である場合の発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3と、混合ガスに接する第1の測温素子62からの電気信号SIと、に依存する。したがって、混合ガスに含まれる分子を構成する炭素原子による単位体積当たりの発熱量QCは、h13を関数を表す記号として、下記38式でも与えられる。QCh13SH1SH2SH3SIさらに、混合ガスに含まれる分子を構成する炭素原子による単位体積当たりの発熱量QCは、混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと、炭素原子の原子量例えば12.0107と、それぞれの分子の体積比と、の積の総和と相関する。したがって、上記38式に基づき、混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと、炭素原子の原子量と、それぞれの分子の体積比と、の積の総和GCは、h14を関数を表す記号として、下記39式で与えられる。GCh14SH1SH2SH3SIまた、n種類のガス成分からなる混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと、炭素原子の原子量と、それぞれの分子の体積比と、の積の総和GCは、下記40式から求められる。GCh15SH1SH2SH3SHn1SI例えば、混合ガスがメタンCH4を90体積、エタンC2H6を5体積、プロパンC3H8を1体積、ブタンC4H10を1体積、窒素N2を1体積、二酸化炭素CO2を2含むとする。この場合、メタンを構成する炭素原子の数は1であるから、メタンを構成する炭素原子の数と、炭素原子の原子量と、メタンの体積比と、の積は、1×12.0107×0.9=10.80963である。また、エタンを構成する炭素原子の数は2であるから、エタンを構成する炭素原子の数と、炭素原子の原子量と、エタンの体積比と、の積は、2×12.0107×0.05=1.20107である。また、プロパンを構成する炭素原子の数は3であるから、プロパンを構成する炭素原子の数と、炭素原子の原子量と、プロパンの体積比と、の積は、3×12.0107×0.01=0.360321である。また、ブタンを構成する炭素原子の数は4であるから、ブタンを構成する炭素原子の数と、炭素原子の原子量と、ブタンの体積比と、の積は、4×12.0107×0.01=0.480428である。窒素を構成する炭素原子の数は0である。二酸化炭素を構成する炭素原子の数は1であるから、二酸化炭素を構成する炭素原子の数と、炭素原子の原子量と、二酸化炭素の体積比と、の積は、1×12.0107×0.02=0.240214である。したがって、メタンを構成する炭素原子の数と炭素原子の原子量とメタンの体積比との積と、エタンを構成する炭素原子の数と炭素原子の原子量とエタンの体積比との積と、プロパンを構成する炭素原子の数と炭素原子の原子量とプロパンの体積比との積と、ブタンを構成する炭素原子の数と炭素原子の原子量とブタンの体積比との積と、二酸化炭素を構成する炭素原子の数と炭素原子の原子量と二酸化炭素の体積比との積と、の総和は、10.809631.201070.3603210.4804280.240214=13.091663である。なお、エタンC2H6、ブタンC4H10、ペンタンC5H12、ヘキサンC6H14を、上記32乃至35式に示すように、それぞれ所定の係数を掛けられたメタンCH4とプロパンC3H8の混合物とみなして39式を算出してもかまわない。第2の実施の形態において、図8に示す式作成部302は、例えば第1のサンプル混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと、炭素原子の原子量と、それぞれの分子の混合ガスにおける体積比と、の積の総和の既知の値と、第2のサンプル混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと、炭素原子の原子量と、それぞれの分子の混合ガスにおける体積比と、の積の総和の既知の値と、第3のサンプル混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと、炭素原子の原子量と、それぞれの分子の混合ガスにおける体積比と、の積の総和の既知の値と、第4のサンプル混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと、炭素原子の原子量と、それぞれの分子の混合ガスにおける体積比と、の積の総和の既知の値と、第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の複数の測定値と、を収集する。さらに式作成部302は、収集した積の総和の値と、電気信号SIの値と、電気信号SHの値と、に基づいて、多変量解析により、第1の測温素子62からの電気信号SI及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3を独立変数とし、ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと、炭素原子の原子量と、それぞれの分子の体積比と、の積の総和を従属変数とする式を算出する。式作成部302は、作成した式を式記憶装置402に保存する。また、第2の実施の形態において、積算出部305は、式記憶装置402に保存されている式の第1の測温素子62からの電気信号SIの独立変数、及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の独立変数に、第1の測温素子62からの電気信号SIの測定値、及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の測定値をそれぞれ代入し、測定対象混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと、炭素原子の原子量と、それぞれの分子の混合ガスにおける体積比と、の積の総和の値を算出する。積算出部305は、積の総和の値を、積記憶装置403に保存する。以上説明した本発明の第2の実施の形態によれば、測定対象混合ガスに接する第1の測温素子62からの電気信号SIの値と、測定対象混合ガスに接する発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の値と、から、測定対象混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと、炭素原子の原子量と、それぞれの分子の体積比と、の積の総和の値を測定することが可能となる。第2の実施の形態の実施例第1の実施の形態の第1の実施例と同様に、組成が既知の40種類のサンプル混合ガスを用意した。次に、40種類のサンプル混合ガスのそれぞれを用いて、図1に示す第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3,SH4TH4,SH5TH5の複数の測定値と、を取得した。その後、40種類のサンプル混合ガスのそれぞれに含まれるそれぞれの分子を構成する炭素原子の数と、炭素原子の原子量と、それぞれの分子の体積比と、の積の総和の既知の値と、第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3,SH4TH4,SH5TH5の複数の測定値と、に基づいて、サポートベクトル回帰により、第1の測温素子62からの電気信号SI及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3,SH4TH4,SH5TH5を独立変数とし、ガスに含まれるそれぞれの分子を構成する炭素原子の数と、炭素原子の原子量と、それぞれの分子の体積比と、の積の総和を従属変数とする方程式を作成した。方程式を作成する際には、キャリブレーションポイントは、3乃至5個を目安に、適宜決定できる。作成された方程式を用いて、40種類のサンプル混合ガスのそれぞれに含まれるそれぞれの分子を構成する炭素原子の数と、炭素原子の原子量と、それぞれの分子の体積比と、の積の総和を算出し、真の値と比較したところ、図14に示すように、誤差は±1.5であった。第3の実施の形態図1乃至図4に示すマイクロチップ8の雰囲気ガスが混合ガスであり、混合ガスが、ガスA、ガスB、ガスC、及びガスDの4種類のガス成分からなっているとする。ガスAに含まれる炭素原子及び水素原子による単位体積当たりの発熱量をKACH、ガスBに含まれる炭素原子及び水素原子による単位体積当たりの発熱量をKBCH、ガスCに含まれる炭素原子及び水素原子による単位体積当たりの発熱量をKCCH、ガスDに含まれる炭素原子及び水素原子による単位体積当たりの発熱量をKDCHとすると、混合ガスに含まれる炭素原子及び水素原子による単位体積当たりの発熱量QCHは、各ガス成分の体積率に、各ガス成分に含まれる炭素原子及び水素原子による単位体積当たりの発熱量を乗じたものの総和で与えられる。したがって、混合ガスに含まれる炭素原子及び水素原子による単位体積当たりの発熱量QCHは、下記41式で与えられる。QCHKACH×VAKBCH×VBKCCH×VCKDCH×VD上記41式、及び18乃至21式より、混合ガスに含まれる分子を構成する炭素原子及び水素原子による単位体積当たりの発熱量QCHは、g4を関数を表す記号として、下記42式で与えられる。QCHg4MI1MI2MI3また、上述したように、混合ガスの放熱係数MIは、発熱素子61の発熱温度がTH1,TH2,TH3である場合の発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3と、混合ガスに接する第1の測温素子62からの電気信号SIと、に依存する。したがって、混合ガスに含まれる分子を構成する炭素原子及び水素原子による単位体積当たりの発熱量QCHは、h16を関数を表す記号として、下記43式でも与えられる。QCHh16SH1SH2SH3SIさらに、混合ガスに含まれる分子を構成する炭素原子及び水素原子による単位体積当たりの発熱量QCHは、混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと炭素原子の原子量とそれぞれの分子の体積比との積の総和と、混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと水素原子の原子量とそれぞれの分子の体積比との積の総和と、の和と相関する。したがって、上記43式に基づき、混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと炭素原子の原子量とそれぞれの分子の体積比との積の総和と、それぞれの分子を構成する水素原子の数NHと水素原子の原子量とそれぞれの分子の体積比との積の総和と、の和GCHは、h17を関数を表す記号として、下記44式で与えられる。GCHh17SH1SH2SH3SIまた、n種類のガス成分からなる混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと炭素原子の原子量とそれぞれの分子の体積比との積の総和と、それぞれの分子を構成する水素原子の数NHと水素原子の原子量とそれぞれの分子の体積比との積の総和と、の和GCHは、下記45式から求められる。GCHh18SH1SH2SH3SHn1SIなお、エタンC2H6、ブタンC4H10、ペンタンC5H12、ヘキサンC6H14を、上記32乃至35式に示すように、それぞれ所定の係数を掛けられたメタンCH4とプロパンC3H8の混合物とみなして45式を算出してもかまわない。第3の実施の形態において、図8に示す式作成部302は、第1のサンプル混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと炭素原子の原子量とそれぞれの分子の体積比との積の総和の既知の値と、それぞれの分子を構成する水素原子の数NHと水素原子の原子量とそれぞれの分子の体積比との積の総和の既知の値と、を収集する。また、式作成部302は、第2のサンプル混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと炭素原子の原子量とそれぞれの分子の体積比との積の総和の既知の値と、それぞれの分子を構成する水素原子の数NHと水素原子の原子量とそれぞれの分子の体積比との積の総和の既知の値と、を収集する。また、式作成部302は、第3のサンプル混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと炭素原子の原子量とそれぞれの分子の体積比との積の総和の既知の値と、それぞれの分子を構成する水素原子の数NHと水素原子の原子量とそれぞれの分子の体積比との積の総和の既知の値と、を収集する。また、式作成部302は、第4のサンプル混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと炭素原子の原子量とそれぞれの分子の体積比との積の総和の既知の値と、それぞれの分子を構成する水素原子の数NHと水素原子の原子量とそれぞれの分子の体積比との積の総和の既知の値と、を収集する。さらに、式作成部302は、第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の複数の測定値と、を収集する。式作成部302は、収集した積の総和の値と、電気信号SIの値と、電気信号SHの値と、に基づいて、多変量解析により、第1の測温素子62からの電気信号SI及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3を独立変数とし、ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと炭素原子の原子量とそれぞれの分子の体積比との積の総和と、それぞれの分子を構成する水素原子の数NHと水素原子の原子量とそれぞれの分子の体積比との積の総和と、の和を従属変数とする式を算出する。式作成部302は、作成した式を式記憶装置402に保存する。また、第3の実施の形態において、積算出部305は、式記憶装置402に保存されている式の第1の測温素子62からの電気信号SIの独立変数、及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の独立変数に、第1の測温素子62からの電気信号SIの測定値、及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の測定値をそれぞれ代入し、測定対象混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと炭素原子の原子量とそれぞれの分子の体積比との積の総和と、それぞれの分子を構成する水素原子の数NHと水素原子の原子量とそれぞれの分子の体積比との積の総和と、の和の値を算出する。積算出部305は、算出した和の値を、積記憶装置403に保存する。以上説明した本発明の第3の実施の形態によれば、測定対象混合ガスに接する第1の測温素子62からの電気信号SIの値と、測定対象混合ガスに接する発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の値と、から、測定対象混合ガスに含まれるそれぞれの分子を構成する炭素原子の数NCと炭素原子の原子量とそれぞれの分子の体積比との積の総和と、それぞれの分子を構成する水素原子の数NHと水素原子の原子量とそれぞれの分子の体積比との積の総和と、の和の値を測定することが可能となる。第4の実施の形態第4の実施の形態に係る発電システムは、図15に示すように、パイプ101で接続されたガス測定システム20と、流量制御装置501と、改質器502と、シフト器503と、選択酸化器504と、燃料電池505と、を備える。ガス測定システム20は、ガスを供給され、第1の実施の形態で説明したように、パイプ101を流れるガスに含まれるそれぞれの分子を構成する水素原子の数と、水素原子の原子量と、それぞれの分子の体積比と、の積の総和の値を算出する。ガス測定システム20の下流に配置された流量制御装置501は、ガス測定システム20が算出した積の総和の値に基づいて、パイプ101を流れるガスの流量を制御する。例えば、ガス測定システム20が算出した積の総和の値が大きい場合は、燃料電池505に供給可能な水素分子が豊富であるため、パイプ101を流れるガスの流量を減らしてもよい。また、ガス測定システム20が算出した積の総和の値が小さい場合は、燃料電池505に供給可能な水素分子が不足する可能性があるため、パイプ101を流れるガスの流量を増やしてもよい。流量制御装置501の下流に配置された改質器502は、パイプ101を流れるガスに水を加え、水蒸気改質と呼ばれる改質法によって、水素分子を生成する。例えば、ガス中のメタンは、水と反応させられ、一酸化炭素、二酸化炭素、及び水素が生成される。改質器502の下流に配置されたシフト器503は、ガス中の一酸化炭素と水を反応させ、二酸化炭素と水素分子とを生成するシフト反応によって、ガス中の一酸化炭素濃度を低減させる。シフト器503の下流に配置された選択酸化器504は、ガス中に残留する一酸化炭素と、酸素と、を反応させて、二酸化炭素を生成し、ガス中の一酸化炭素濃度をさらに低下させる。シフト器503の下流に配置された燃料電池505は、水素分子を豊富に含み、一酸化炭素濃度が低減されたガスを供給され、発電する。以上説明した第4の実施の形態に係る発電システムによれば、燃料電池505に供給される水素分子の量が予測可能となり、また、燃料電池505に供給される水素分子の量を一定に保つことが可能となる。そのため、燃料電池505を安定に駆動させることが可能となる。なお、ガス測定システム20は、ガスに含まれるそれぞれの分子を構成する炭素原子の数と、炭素原子の原子量と、それぞれの分子の体積比と、の積の総和をさらに算出してもよい。その他の実施の形態上記のように、本発明は実施の形態によって記載したが、この開示の一部をなす記述及び図面はこの発明を限定するものであると理解するべきではない。この開示から当業者には様々な代替実施の形態、実施の形態及び運用技術が明らかになるはずである。例えば、実施の形態では、図8に示す式記憶装置402が、図1に示す第1の測温素子62からの電気信号及び複数の発熱温度における発熱素子61からの電気信号を独立変数とし、ガスに含まれるそれぞれの分子を構成する水素原子の数と水素原子の原子量とそれぞれの分子の体積比との積の総和を従属変数とする式を保存する例を説明した。これに対し、上記24式で説明したように、混合ガスに含まれる分子を構成する水素原子による単位体積当たりの発熱量QHは、発熱素子61の温度がそれぞれTH1,TH2,TH3である場合の気体の放熱係数MI1TH1,MI2TH2,MI3TH3を変数とする方程式でも与えられる。また、混合ガスに含まれる分子を構成する水素原子による単位体積当たりの発熱量QHは、混合ガスに含まれるそれぞれの分子を構成する水素原子の数NHと、水素原子の原子量例えば、1.00794と、それぞれの分子の混合ガスにおける体積比と、の積の総和と相関する。よって、図8に示す式記憶装置402は、発熱素子61の複数の発熱温度における気体の放熱係数を独立変数とし、ガスに含まれるそれぞれの分子を構成する水素原子の数と水素原子の原子量とそれぞれの分子の体積比との積の総和を従属変数とする式等の、放熱係数と、積の総和と、の相関関係を保存してもよい。この場合、測定部301は、パイプ101に注入された気体の放熱係数の測定値を、発熱素子61を複数の発熱温度に発熱させて測定する。なお、気体の放熱係数がマイクロチップ8を用いて測定可能であることは、上記9式で説明したとおりである。積算出部305は、式記憶装置402に保存されている式の独立変数に、気体の放熱係数の測定値を代入して、ガスに含まれるそれぞれの分子を構成する水素原子の数と水素原子の原子量とそれぞれの分子の体積比との積の総和の測定値を算出する。次に、図16は、発熱抵抗体に2mA、2.5mA、及び3mAの電流を流した際の、混合ガスの放熱係数と、熱伝導率と、の関係を示す。図16に示すように、混合ガスの放熱係数と、熱伝導率と、は一般に比例関係にある。したがって、図8に示す式記憶装置402は、発熱素子61の複数の発熱温度における気体の熱伝導率を独立変数とし、ガスに含まれるそれぞれの分子を構成する水素原子の数と水素原子の原子量とそれぞれの分子の体積比との積の総和を従属変数とする式を保存してもよい。この場合、測定部301は、パイプ101に注入された気体の熱伝導率の測定値を、発熱素子61を複数の発熱温度に発熱させて測定する。積算出部305は、式記憶装置402に保存されている式の独立変数に、気体の熱伝導率の測定値を代入して、ガスに含まれるそれぞれの分子を構成する水素原子の数と水素原子の原子量とそれぞれの分子の体積比との積の総和の測定値を算出する。この様に、本発明はここでは記載していない様々な実施の形態等を包含するということを理解すべきである。"