"気体が流れるパイプと、前記パイプ内を流れる前記気体の流量を制御する流量制御装置と、前記パイプに配置された測温素子と、前記パイプに配置された、複数の発熱温度で発熱する発熱素子と、前記パイプ内を流れている前記気体の温度に依存する前記測温素子からの電気信号の値と、前記複数の発熱温度のそれぞれにおける前記発熱素子からの電気信号の値と、を測定する測定部と、前記測温素子からの電気信号及び前記複数の発熱温度における前記発熱素子からの電気信号を独立変数とし、前記発熱量を従属変数とする発熱量算出式を保存する式記憶装置と、前記発熱量算出式の独立変数に、前記測温素子からの電気信号の値、及び前記前記発熱素子からの電気信号の値を代入し、前記気体の発熱量の値を算出する発熱量算出部と、を備える、発熱量測定システム。前記測温素子及び前記発熱素子が、前記パイプ内に突き出るように配置されている、請求項1に記載の発熱量測定システム。前記測温素子及び前記発熱素子が、前記パイプ内を流れている前記気体に対向するよう、前記パイプ内に斜めに突き出るように配置されている、請求項1に記載の発熱量測定システム。複数種類のガス成分を含む複数のサンプル混合ガスの発熱量の値と、前記複数のサンプル混合ガスのそれぞれに接する前記発熱素子からの電気信号の値とに基づいて、前記発熱量算出式が作成された、請求項1乃至3のいずれか1項に記載の発熱量測定システム。前記発熱量算出式を作成するために、サポートベクトル回帰が用いられた、請求項4に記載の発熱量測定システム。前記気体が天然ガスである、請求項1乃至5のいずれか1項に記載の発熱量測定システム。パイプに気体を流すことと、前記パイプ内を流れている前記気体の温度に依存する測温素子からの電気信号の値を得ることと、前記パイプ内を流れている前記気体に接する発熱素子を複数の発熱温度で発熱させることと、前記複数の発熱温度のそれぞれにおける前記発熱素子からの電気信号の値を得ることと、前記測温素子からの電気信号及び前記複数の発熱温度における前記発熱素子からの電気信号を独立変数とし、前記発熱量を従属変数とする発熱量算出式を用意することと、前記発熱量算出式の独立変数に、前記測温素子からの電気信号の値、及び前記発熱素子からの電気信号の値を代入し、前記気体の発熱量の値を算出することと、を含む、発熱量の測定方法。前記測温素子及び前記発熱素子が、前記パイプ内に突き出るように配置されている、請求項7に記載の発熱量の測定方法。前記測温素子及び前記発熱素子が、前記パイプ内を流れている前記気体に対向するよう、前記パイプ内に斜めに突き出るように配置されている、請求項7に記載の発熱量の測定方法。複数種類のガス成分を含む複数のサンプル混合ガスの発熱量の値と、前記複数のサンプル混合ガスのそれぞれに接する前記発熱素子からの電気信号の値とに基づいて、前記発熱量算出式が作成された、請求項7乃至9のいずれか1項に記載の発熱量の測定方法。前記発熱量算出式を作成するために、サポートベクトル回帰が用いられた、請求項10に記載の発熱量の測定方法。前記気体が天然ガスである、請求項7乃至11のいずれか1項に記載の発熱量の測定方法。気体が流れるパイプと、前記パイプ内を流れる前記気体の流量を制御する流量制御装置と、前記パイプ中を流れる気体の放熱係数又は熱伝導率の測定値を測定する測定部と、前記放熱係数又は熱伝導率と、発熱量と、の相関関係を保存する記憶装置と、前記気体の放熱係数又は熱伝導率の測定値と、前記相関関係と、に基づいて、前記気体の発熱量の測定値を算出する発熱量算出部と、を備える、発熱量測定システム。前記放熱係数又は熱伝導率を測定するマイクロチップが、前記パイプ内に突き出るように配置されている、請求項13に記載の発熱量測定システム。前記放熱係数又は熱伝導率を測定するマイクロチップが、前記パイプ内を流れている前記気体に対向するよう、前記パイプ内に斜めに突き出るように配置されている、請求項13に記載の発熱量測定システム。複数種類のガス成分を含む複数のサンプル混合ガスの発熱量の値と、前記複数のサンプル混合ガスのそれぞれの放熱係数又は熱伝導率と、に基づいて、前記相関関係を表す発熱量算出式が作成された、請求項13乃至15のいずれか1項に記載の発熱量測定システム。前記発熱量算出式を作成するために、サポートベクトル回帰が用いられた、請求項16に記載の発熱量測定システム。前記複数のサンプル混合ガスのそれぞれが天然ガスである、請求項16又は17に記載の発熱量測定システム。前記気体が天然ガスである、請求項15乃至18のいずれか1項に記載の発熱量測定システム。しかし、特許文献1に開示された方法は、熱伝導率を測定するためのセンサの他に、音速を測定するための高価な音速センサが必要である。そこで、本発明は、ガスの発熱量を容易に測定可能な発熱量測定システム及び発熱量の測定方法を提供することを目的の一つとする。本発明の態様によれば、a気体が流れるパイプと、bパイプ内を流れる気体の流量を制御する流量制御装置と、cパイプに配置された測温素子と、dパイプに配置された、複数の発熱温度で発熱する発熱素子と、eパイプ内を流れている気体の温度に依存する測温素子からの電気信号の値と、複数の発熱温度のそれぞれにおける発熱素子からの電気信号の値と、を測定する測定部と、f測温素子からの電気信号及び複数の発熱温度における発熱素子からの電気信号を独立変数とし、発熱量を従属変数とする発熱量算出式を保存する式記憶装置と、g発熱量算出式の測温素子からの電気信号の独立変数、及び発熱素子からの電気信号の独立変数に、測温素子からの電気信号の値、及び発熱素子からの電気信号の値を代入し、気体の発熱量の値を算出する発熱量算出部と、を備える、発熱量測定システムが提供される。本発明の態様によれば、aパイプに気体を流すことと、bパイプ内を流れている気体の温度に依存する測温素子からの電気信号の値を得ることと、cパイプ内を流れている気体に接する発熱素子を複数の発熱温度で発熱させることと、d複数の発熱温度のそれぞれにおける発熱素子からの電気信号の値を得ることと、e測温素子からの電気信号及び複数の発熱温度における発熱素子からの電気信号を独立変数とし、発熱量を従属変数とする発熱量算出式を用意することと、f発熱量算出式の測温素子からの電気信号の独立変数、及び発熱素子からの電気信号の独立変数に、測温素子からの電気信号の値、及び発熱素子からの電気信号の値を代入し、気体の発熱量の値を算出することと、を含む、発熱量の測定方法が提供される。本発明の態様によれば、a気体が流れるパイプと、bパイプ内を流れる気体の流量を制御する流量制御装置と、cパイプ中を流れる気体の放熱係数又は熱伝導率の測定値を測定する測定部と、d放熱係数又は熱伝導率と、発熱量と、の相関関係を保存する記憶装置と、e気体の放熱係数又は熱伝導率の測定値と、相関関係と、に基づいて、気体の発熱量の測定値を算出する発熱量算出部と、を備える、発熱量測定システムが提供される。本発明によれば、ガスの発熱量を容易に測定可能な発熱量測定システム及び発熱量の測定方法を提供可能である。以下に本発明の実施の形態を説明する。以下の図面の記載において、同一又は類似の部分には同一又は類似の符号で表している。但し、図面は模式的なものである。したがって、具体的な寸法等は以下の説明を照らし合わせて判断するべきものである。また、図面相互間においても互いの寸法の関係や比率が異なる部分が含まれていることは勿論である。まず、斜視図である図1、及びII-II方向から見た断面図である図2を参照して、実施の形態に係る発熱量測定システムに用いられるマイクロチップ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の放熱係数をMA、ガスBの放熱係数をMB、ガスCの放熱係数をMC、ガスDの放熱係数をMDとすると、混合ガスの放熱係数MIは、各ガス成分の体積率に、各ガス成分の放熱係数を乗じたものの総和で与えられる。したがって、混合ガスの放熱係数MIは、下記12式で与えられる。MIMA×VAMB×VBMC×VCMD×VDさらに、ガスの放熱係数は発熱素子61の温度THに依存するので、混合ガスの放熱係数MIは、発熱素子61の温度THの関数として、下記13式で与えられる。MIMA×VAMB×VBMC×VCMD×VDしたがって、発熱素子61の温度がTH1のときの混合ガスの放熱係数MI1TH1は下記14式で与えられる。また、発熱素子61の温度がTH2のときの混合ガスの放熱係数MI2TH2は下記15式で与えられ、発熱素子61の温度がTH3のときの混合ガスの放熱係数MI3TH3は下記16式で与えられる。MI1MA×VAMB×VBMC×VCMD×VDMI2MA×VAMB×VBMC×VCMD×VDMI3MA×VAMB×VBMC×VCMD×VDここで、発熱素子61の温度THに対して各ガス成分の放熱係数MATH,MBTH,MCTH,MDTHが非線形性を有する場合、上記14乃至16式は、線形独立な関係を有する。また、発熱素子61の温度THに対して各ガス成分の放熱係数MATH,MBTH,MCTH,MDTHが線形性を有する場合でも、発熱素子61の温度THに対する各ガス成分の放熱係数MATH,MBTH,MCTH,MDTHの変化率が異なる場合は、上記14乃至16式は、線形独立な関係を有する。さらに、14乃至16式が線形独立な関係を有する場合、10及び14乃至16式は線形独立な関係を有する。図7は、天然ガスに含まれるメタンCH4、プロパンC3H8、窒素N2、及び二酸化炭素CO2の放熱係数と、発熱抵抗体である発熱素子61の温度との関係を示すグラフである。発熱素子61の温度に対して、メタンCH4、プロパンC3H8、窒素N2、及び二酸化炭素CO2のそれぞれのガス成分の放熱係数は線形性を有する。しかし、発熱素子61の温度に対する放熱係数の変化率は、メタンCH4、プロパンC3H8、窒素N2、及び二酸化炭素CO2のそれぞれで異なる。したがって、混合ガスを構成するガス成分がメタンCH4、プロパンC3H8、窒素N2、及び二酸化炭素CO2であるである場合、上記14乃至16式は、線形独立な関係を有する。14乃至16式中の各ガス成分の放熱係数MATH1,MBTH1,MCTH1,MDTH1,MATH2,MBTH2,MCTH2,MDTH2,MATH3,MBTH3,MCTH3,MDTH3の値は、測定等により予め得ることが可能である。したがって、10及び14乃至16式の連立方程式を解くと、ガスAの体積率VA、ガスBの体積率VB、ガスCの体積率VC、及びガスDの体積率VDのそれぞれが、下記17乃至20式に示すように、混合ガスの放熱係数MI1TH1,MI2TH2,MI3TH3の関数として与えられる。なお、下記17乃至20式において、nを自然数として、fnは関数を表す記号である。VAf1MI1MI2MI3VBf2MI1MI2MI3VCf3MI1MI2MI3VDf4MI1MI2MI3ここで、上記11式に17乃至20式を代入することにより、下記21式が得られる。QKA×VAKB×VBKC×VCKD×VDKA×f1MI1MI2MI3KB×f2MI1MI2MI3KC×f3MI1MI2MI3KD×f4MI1MI2MI3上記21式に示すように、混合ガスの単位体積当たりの発熱量Qは、発熱素子61の温度がTH1,TH2,TH3である場合の混合ガスの放熱係数MI1TH1,MI2TH2,MI3TH3を変数とする方程式で与えられる。したがって、混合ガスの発熱量Qは、gを関数を表す記号として、下記22式で与えられる。QgMI1MI2MI3よって、ガスA、ガスB、ガスC、及びガスDからなる混合ガスについて、予め上記22式を得れば、ガスAの体積率VA、ガスBの体積率VB、ガスCの体積率VC、及びガスDの体積率VDが未知の測定対象混合ガスの単位体積当たりの発熱量Qを容易に算出可能であることを、発明者は見出した。具体的には、発熱素子61の発熱温度がTH1,TH2,TH3である場合の測定対象混合ガスの放熱係数MI1TH1,MI2TH2,MI3TH3を測定し、22式に代入することにより、測定対象混合ガスの発熱量Qを一意に求めることが可能となる。また、混合ガスの放熱係数MIは、上記9式に示すように、発熱素子61の抵抗値RHと、第1の測温素子62の抵抗値RIと、に依存する。そこで、本発明者は、混合ガスの単位体積当たりの発熱量Qは、下記23式に示すように、発熱素子61の温度がTH1,TH2,TH3である場合の発熱素子61の抵抗値RH1TH1,RH2TH2,RH3TH3と、混合ガスに接する第1の測温素子62の抵抗値RIと、を変数とする方程式でも与えられることを見出した。QgRH1RH2RH3RIよって、測定対象混合ガスに接する発熱素子61の発熱温度がTH1,TH2,TH3である場合の発熱素子61の抵抗値RH1TH1,RH2TH2,RH3TH3と、測定対象混合ガスに接する第1の測温素子62の抵抗値RIを測定し、23式に代入することによっても、測定対象混合ガスの発熱量Qを一意に求めることが可能となる。また、混合ガスの単位体積当たりの発熱量Qは、下記24式に示すように、発熱素子61の温度がTH1,TH2,TH3である場合の発熱素子61の通電電流IH1TH1,IH2TH2,IH3TH3と、混合ガスに接する第1の測温素子62の通電電流IIと、を変数とする方程式でも与えられる。QgIH1IH2IH3IIあるいは混合ガスの単位体積当たりの発熱量Qは、下記25式に示すように、発熱素子61の温度がTH1,TH2,TH3である場合の発熱素子61にかかる電圧VH1TH1,VH2TH2,VH3TH3と、混合ガスに接する第1の測温素子62にかかる電圧VIと、を変数とする方程式でも与えられる。QgVH1VH2VH3VIまたあるいは混合ガスの単位体積当たりの発熱量Qは、下記26式に示すように、発熱素子61の温度がTH1,TH2,TH3である場合の発熱素子61に接続されたアナログ-デジタル変換回路以下においてA/D変換回路という。の出力信号ADH1TH1,ADH2TH2,ADH3TH3と、混合ガスに接する第1の測温素子62に接続されたA/D変換回路の出力信号ADIと、を変数とする方程式でも与えられる。QgADH1ADH2ADH3ADIよって、混合ガスの単位体積当たりの発熱量Qは、下記27式に示すように、発熱素子61の発熱温度がTH1,TH2,TH3である場合の発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3と、混合ガスに接する第1の測温素子62からの電気信号SIと、を変数とする方程式で与えられる。QgSH1SH2SH3SIなお、混合ガスのガス成分は、4種類に限定されることはない。例えば、混合ガスがn種類のガス成分からなる場合、まず、下記28式で与えられる、少なくとも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の値と、を測定し、28式に代入することにより、測定対象混合ガスの単位体積当たりの発熱量Qを一意に求めることが可能となる。QgSH1SH2SH3SHn1SIただし、混合ガスが、ガス成分としてメタンCH4、プロパンC3H8に加えて、jを自然数として、メタンCH4とプロパンC3H8以外のアルカンCjH2j2を含む場合、メタンCH4とプロパンC3H8以外のアルカンCjH2j2を、メタンCH4とプロパンC3H8の混合物とみなしても、28式の算出には影響しない。例えば、エタンC2H6、ブタンC4H10、ペンタンC5H12、ヘキサンC6H14を、下記29乃至32式に示すように、それぞれ所定の係数を掛けられたメタンCH4とプロパンC3H8の混合物とみなして28式を算出してもかまわない。C2H605CH405C3H8C4H1005CH415C3H8C5H1210CH420C3H8C6H1415CH425C3H8したがって、zを自然数として、n種類のガス成分からなる混合ガスが、ガス成分としてメタンCH4、プロパンC3H8に加えて、メタンCH4とプロパンC3H8以外のz種類のアルカンCjH2j2を含む場合、少なくともn-z-1種類の発熱温度における発熱素子61からの電気信号SHと、第1の測温素子62からの電気信号SIと、を変数とする方程式を求めてもよい。なお、28式の算出に用いられた混合ガスのガス成分の種類と、単位体積当たりの発熱量Qが未知の測定対象混合ガスのガス成分の種類が同じ場合に、測定対象混合ガスの発熱量Qの算出に28式を利用可能であることはもちろんである。さらに、測定対象混合ガスがn種類より少ない種類のガス成分からなり、かつ、n種類より少ない種類のガス成分が、28式の算出に用いられた混合ガスに含まれている場合も、28式を利用可能である。例えば、28式の算出に用いられた混合ガスが、メタンCH4、プロパンC3H8、窒素N2、及び二酸化炭素CO2の4種類のガス成分を含む場合、測定対象混合ガスが、窒素N2を含まず、メタンCH4、プロパンC3H8、及び二酸化炭素CO2の3種類のガス成分のみを含む場合も、測定対象混合ガスの発熱量Qの算出に28式を利用可能である。さらに、28式の算出に用いられた混合ガスが、ガス成分としてメタンCH4とプロパンC3H8を含む場合、測定対象混合ガスが、28式の算出に用いられた混合ガスに含まれていないアルカンCjH2j2を含んでいても、28式を利用可能である。これは、上述したように、メタンCH4とプロパンC3H8以外のアルカンCjH2j2を、メタンCH4とプロパンC3H8の混合物とみなしても、28式を用いた単位体積当たりの発熱量Qの算出に影響しないためである。ここで、図8及び図9に示す実施の形態に係る発熱量測定システム20は、複数のサンプル混合ガスのそれぞれが流れるパイプ101と、パイプ101に配置され、図1に示す第1の測温素子62及び複数の発熱温度THで発熱する発熱素子61を含むマイクロチップ8と、を備える。図8に示すように、マイクロチップ8は、断熱部材18を介してパイプ101内に突き出るように配置されている。なお、少なくとも図1に示す第1の測温素子62及び発熱素子61が図8に示すパイプ101内に突き出ていれば、マイクロチップ8は部分的にパイプ101の側壁に埋め込まれていてもよい。マイクロチップ8をパイプ101内に突き出すように配置することにより、図2に示すマイクロチップ8のキャビティ66内のガスの置換が速くなり、マイクロチップ8の応答性が向上する。また、図10に示すように、マイクロチップ8を、パイプ101の側壁に対して斜めに配置してもよい。この場合、図10に示すように、マイクロチップ8を含むパッケージ118全体を、パイプ101の側壁に対して斜めに配置してもよいし、図11に示すように、マイクロチップ8をパッケージ118に対して斜めにダイボンドすることにより、マイクロチップ8をパイプ101の側壁に対して斜めに配置してもよい。マイクロチップ8をパイプ101の側壁に対して斜めに配置することによって、マイクロチップ8がパイプ101内を流れてくるガスに対向するため、図2に示すマイクロチップ8のキャビティ66内のガスが効率良く置換する。図4に示すマイクロチップ8においても同様である。図8に示す発熱量測定システム20は、さらに、複数のサンプル混合ガスのそれぞれの温度TIに依存する第1の測温素子62からの電気信号SIの値と、複数の発熱温度THのそれぞれにおける発熱素子61からの電気信号SHの値と、を測定する測定部301と、複数の混合ガスの既知の発熱量Qの値、第1の測温素子62からの電気信号SIの値、及び複数の発熱温度における発熱素子61からの電気信号の値に基づいて、第1の測温素子62からの電気信号SI及び複数の発熱温度THにおける発熱素子61からの電気信号SHを独立変数とし、発熱量Qを従属変数とする発熱量算出式を作成する式作成部と、を備える。なお、サンプル混合ガスは、複数種類のガス成分を含む。それぞれ発熱量Qが異なる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を与えられた場合、パイプ101内を留まることなく流れる第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乃至第4のサンプル混合ガスのそれぞれの既知の発熱量Qの値と、第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の複数の測定値と、を収集する。さらに式作成部302は、収集した発熱量Q、電気信号SI、及び電気信号SHの値に基づいて、多変量解析により、第1の測温素子62からの電気信号SI及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3を独立変数とし、発熱量Qを従属変数とする発熱量算出式を算出する。なお、多変量解析とは、A.JSmola及びB.Scholkopf著のATutorialonSupportVectorRegressionNeuroCOLTTechnicalReportNC-TR-98-030、1998年に開示されているサポートベクトル回帰、重回帰分析、及び特開平5-141999号公報に開示されているファジィ数量化理論II類等を含む。発熱量測定システム20は、CPU300に接続された式記憶装置402をさらに備える。式記憶装置402は、式作成部302が作成した発熱量算出式を保存する。さらにCPU300には、入力装置312及び出力装置313が接続される。入力装置312としては、例えばキーボード、及びマウス等のポインティングデバイス等が使用可能である。出力装置313には液晶ディスプレイ、モニタ等の画像表示装置、及びプリンタ等が使用可能である。次に、図12に示すフローチャートを用いて実施の形態に係る発熱量算出式の作成方法について説明する。aステップS100で、図9に示す第2乃至第4の流量制御装置32B-32Dの弁を閉じたまま、第1の流量制御装置32Aの弁を開き、図8に示すパイプ101内に第1のサンプル混合ガスを導入する。ステップS101で、測定部301は、パイプ101内を留まることなく流れる第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は、パイプ101内を留まることなく流れる第1のサンプル混合ガスに接し、150°Cで発熱する発熱素子61からの電気信号SH2TH2の値を、電気信号記憶装置401に保存する。c再びステップS102で、図1及び図2に示す発熱素子61の温度の切り替えが完了したか否か判定する。温度200°Cへの切り替えが完了していない場合には、ステップS101に戻り、図8に示す駆動回路303は、図1及び図2に示す発熱素子61を200°Cで発熱させる。図8に示す測定部301は、パイプ101内を留まることなく流れる第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は、パイプ101内を留まることなく流れる第2のサンプル混合ガスに接する第1の測温素子62からの電気信号SIの値を測定し、電気信号記憶装置401に保存する。また、測定部301は、パイプ101内を留まることなく流れる第2のサンプル混合ガスに接し、100°C、150°C、及び200°Cで発熱する発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の値を、電気信号記憶装置401に保存する。fその後、ステップS100乃至ステップS103のループが繰り返される。これにより、パイプ101内を留まることなく流れる第3のサンプル混合ガスに接する第1の測温素子62からの電気信号SIの値と、パイプ101内を留まることなく流れる第3のサンプル混合ガスに接し、100°C、150°C、及び200°Cで発熱する発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の値と、が電気信号記憶装置401に保存される。また、パイプ101内を留まることなく流れる第4のサンプル混合ガスに接する第1の測温素子62からの電気信号SIの値と、パイプ101内を留まることなく流れる第4のサンプル混合ガスに接し、100°C、150°C、及び200°Cで発熱する発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の値と、が電気信号記憶装置401に保存される。gステップS104で、入力装置312から式作成部302に、第1のサンプル混合ガスの既知の発熱量Qの値、第2のサンプル混合ガスの既知の発熱量Qの値、第3のサンプル混合ガスの既知の発熱量Qの値、及び第4のサンプル混合ガスの既知の発熱量Qの値を入力する。また、式作成部302は、電気信号記憶装置401から、第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の複数の測定値と、を読み出す。hステップS105で、第1乃至第4のサンプル混合ガスの発熱量Qの値と、第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の複数の測定値と、に基づいて、式作成部302は、重回帰分析を行う。重回帰分析により、式作成部302は、第1の測温素子62からの電気信号SI及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3を独立変数とし、発熱量Qを従属変数とする発熱量算出式を算出する。その後、ステップS106で、式作成部302は作成した発熱量算出式を式記憶装置402に保存し、実施の形態に係る発熱量算出式の作成方法が終了する。以上示したように、実施の形態に係る発熱量算出式の作成方法によれば、測定対象混合ガスの発熱量Qの値を一意に算出可能な発熱量算出式を作成することが可能となる。次に、発熱量Qが未知の測定対象混合ガスの発熱量Qの値を測定する際の、実施の形態に係る発熱量測定システム20の機能を説明する。例えば未知の体積率でメタンCH4、プロパンC3H8、窒素N2、及び二酸化炭素CO2等を含む、発熱量Qが未知の天然ガス等の測定対象混合ガスが、パイプ101に導入される。図1及び図2に示すマイクロチップ8の第1の測温素子62は、パイプ101内を留まることなく流れる測定対象混合ガスの温度に依存する電気信号SIを出力する。次に、発熱素子61は、図8に示す駆動回路303から駆動電力PH1,PH2,PH3を与えられる。駆動電力PH1,PH2,PH3を与えられた場合、パイプ101内を留まることなく流れる測定対象混合ガスに接する発熱素子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と、を独立変数とし、気体の発熱量Qを従属変数とする発熱量算出式を保存している。実施の形態に係る発熱量測定システム20は、さらに、発熱量算出部305を備える。発熱量算出部305は、発熱量算出式の第1の測温素子62からの電気信号SIの独立変数、及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の独立変数に、第1の測温素子62からの電気信号SIの測定値、及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の測定値をそれぞれ代入し、パイプ101内を留まることなく流れる測定対象混合ガスの発熱量Qの測定値を算出する。CPU300には、発熱量記憶装置403がさらに接続されている。発熱量記憶装置403は、発熱量算出部305が算出した測定対象混合ガスの発熱量Qの値を保存する。次に、図13に示すフローチャートを用いて、実施の形態に係る発熱量の測定方法について説明する。aステップS200で、図8に示すパイプ101内に測定対象混合ガスを導入する。ステップS201で、測定部301は、パイプ101内を留まることなく流れる測定対象混合ガスに接する第1の測温素子62からの電気信号SIの値を測定し、電気信号記憶装置401に保存する。次に、駆動回路303は、図1及び図2に示す発熱素子61に駆動電力PH1を与え、発熱素子61を100°Cで発熱させる。図8に示す測定部301は、パイプ101内を留まることなく流れる測定対象混合ガスに接し、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は、パイプ101内を留まることなく流れる測定対象混合ガスに接し、150°Cで発熱する発熱素子61からの電気信号SH2TH2の値を、電気信号記憶装置401に保存する。c再びステップS202で、図1及び図2に示す発熱素子61の温度の切り替えが完了したか否か判定する。温度200°Cへの切り替えが完了していない場合には、ステップS201に戻り、駆動回路303は、図1及び図2に示す発熱素子61に駆動電力PH3を与え、発熱素子61を200°Cで発熱させる。図8に示す測定部301は、パイプ101内を留まることなく流れる測定対象混合ガスに接し、200°Cで発熱する発熱素子61からの電気信号SH3TH3の値を、電気信号記憶装置401に保存する。d発熱素子61の温度の切り替えが完了した場合には、ステップS202からステップS203に進む。ステップS203で、図8に示す発熱量算出部305は、式記憶装置402から、第1の測温素子62からの電気信号SI及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3を独立変数とし、発熱量Qを従属変数とする発熱量算出式を読み出す。また、発熱量算出部305は、電気信号記憶装置401から、測定対象混合ガスに接する第1の測温素子62からの電気信号SIの測定値、及び測定対象混合ガスに接する発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の測定値を読み出す。eステップS204で、発熱量算出部305は、発熱量算出式の電気信号SI及び電気信号SH1TH1,SH2TH2,SH3TH3の独立変数に、それぞれ測定値を代入し、測定対象混合ガスの発熱量Qの値を算出する。その後、発熱量算出部305は、算出した発熱量Qの値を発熱量記憶装置403に保存し、実施の形態に係る発熱量の測定方法を終了する。以上説明した実施の形態に係る発熱量算出方法によれば、高価なガスクロマトグラフィ装置や音速センサを用いることなく、パイプ101内を留まることなく流れる測定対象混合ガスに接する第1の測温素子62からの電気信号SIの値と、パイプ101内を留まることなく流れる測定対象混合ガスに接する発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の値と、から、測定対象混合ガスの混合ガスの発熱量Qの値を測定することが可能となる。天然ガスは、産出ガス田によって炭化水素の成分比率が異なる。また、天然ガスには、炭化水素の他に、窒素N2や炭酸ガスCO2等が含まれる。そのため、産出ガス田によって、天然ガスに含まれるガス成分の体積率は異なり、ガス成分の種類が既知であっても、天然ガスの発熱量Qは未知であることが多い。また、同一のガス田由来の天然ガスであっても、発熱量Qが常に一定であるとは限らず、採取時期によって変化することもある。従来、天然ガスの使用料金を徴収する際には、天然ガスの使用発熱量Qでなく、使用体積に応じて課金する方法がとられている。しかし、天然ガスは由来する産出ガス田によって発熱量Qが異なるため、使用体積に課金するのは公平でない。これに対し、実施の形態に係る発熱量算出方法を用いれば、ガス成分の種類が既知であるが、ガス成分の体積率が未知であるために発熱量Qが未知の天然ガス等の混合ガスの発熱量Qを、簡易に算出することが可能となる。そのため、公平な使用料金を徴収することが可能となる。また、ガスタービンを駆動する際には、ガスタービンに供給する燃料ガスとしての天然ガスの発熱量Qをタイムラグなく監視することが求められる。燃料ガスの発熱量Qが一定でない場合、燃焼振動等でガスタービンが破損することがあるためである。しかし、従来の熱量計は、応答時間が分単位と長く、ガスタービンに供給する燃料ガスの発熱量Qの制御に適さない。これに対し、実施の形態に係る発熱量測定システムは、秒単位で発熱量を測定することが可能であるため、ガスタービンに供給する燃料ガスの発熱量Qの制御にも適する。さらに、実施の形態に係る発熱量算出方法によれば、天然ガス等の混合ガスの正確な発熱量Qを容易に知ることが可能となるため、混合ガスを燃焼させる場合に必要な空気量を適切に設定することが可能となる。そのため、無駄な二酸化炭素CO2の排出量を削減することも可能となる。実施例1まず、発熱量Qの値が既知の23種類のサンプル混合ガスを用意した。23種類のサンプル混合ガスのそれぞれは、ガス成分としてメタンCH4、エタンC2H6、プロパンC3H8、ブタンC4H10、窒素N2、及び二酸化炭素CO2のいずれか又は全部を含んでいた。例えば、あるサンプル混合ガスは、90volのメタン、3volのエタン、1volのプロパン、1volのブタン、4volの窒素、及び1volの二酸化炭素を含んでいた。また、あるサンプル混合ガスは、85volのメタン、10volのエタン、3volのプロパン、及び2volのブタンを含み、窒素及び二酸化炭素を含んでいなかった。また、あるサンプル混合ガスは、85volのメタン、8volのエタン、2volのプロパン、1volのブタン、2volの窒素、及び2volの二酸化炭素を含んでいた。次に、23種類のサンプル混合ガスのそれぞれを用いて、図8に示す第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の複数の測定値と、を取得した。その後、23種類のサンプル混合ガスの既知の発熱量Qの値と、第1の測温素子62からの電気信号SIの複数の測定値と、発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3の複数の測定値と、に基づいて、サポートベクトル回帰により、第1の測温素子62からの電気信号SI及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3を独立変数とし、発熱量Qを従属変数とする、発熱量Qを算出するための1次方程式、2次方程式、及び3次方程式を作成した。発熱量Qを算出するための1次方程式を作成する際には、キャリブレーションポイントは、3乃至5個を目安に、適宜決定できる。作成された1次方程式は下記33乃至35式で与えられた。23種類のサンプル混合ガスの発熱量Qを33乃至35式で算出し、真の発熱量Qと比較したところ、最大誤差は2.1であった。Q401174×VH1179×VH2289×VH3104×VIQ401238×RH1607×RH2228×RH3114×RIQ401174×ADH1179×ADH2289×ADH3104×ADI発熱量Qを算出するための2次方程式を作成する際には、キャリブレーションポイントは、8乃至9個を目安に、適宜決定できる。23種類のサンプル混合ガスの発熱量Qを作成された2次方程式で算出し、真の発熱量Qと比較したところ、最大誤差は1.2乃至1.4であった。発熱量Qを算出するための3次方程式を作成する際には、キャリブレーションポイントは、10乃至14個を目安に、適宜決定できる。23種類のサンプル混合ガスの発熱量Qを作成された3次方程式で算出し、真の発熱量Qと比較したところ、最大誤差は1.2未満であった。実施例2実施例1で使用したサンプル混合ガスと同様に、発熱量Qの値が既知の23種類のサンプル混合ガスを用意した。ここで、発熱素子61で加熱される前のサンプル混合ガスの温度を、-10°C、5°C、23°C、40°C、及び50°Cに設定した。次に、サポートベクトル回帰により、第1の測温素子62からの電気信号SI及び発熱素子61からの電気信号SH1TH1,SH2TH2,SH3TH3を独立変数とし、発熱量Qを従属変数とする、発熱量Qを算出するための3次方程式を作成した。すると、図14乃至図16に示すように、発熱素子61で加熱される前のサンプル混合ガスの温度にかかわらず、算出される発熱量Qの誤差にばらつきが生じなかった。なお、図14の結果は、電気信号Sとして抵抗Rを用いて得られた。図15の結果は、電気信号Sとして電圧Vを用いて得られた。図16の結果は、電気信号SとしてA/D変換回路304からの出力信号ADを用いて得られた。実施例3図17に示すように、マイクロチップ8をパイプ101内に突き出るように配置した、実施例3に係る発熱量測定システムを用意した。また、図18に示すように、パイプ101の側壁に凹部を設け、マイクロチップ8をパイプ101の側壁の凹部に配置した、比較例に係る発熱量測定システムも用意した。次に、実施例3に係る発熱量測定システムと、比較例に係る発熱量測定システムと、を用いて、同一条件下で、測定対象混合ガスの発熱量を算出した。すると、図19に示すように、実施例3に係る発熱量測定システムは、約8秒で発熱量の算出が可能であった。しかし、図20に示すように、比較例に係る発熱量測定システムは、約16秒で発熱量の算出が可能であった。したがって、マイクロチップ8をパイプ101内に突き出るように配置することにより、マイクロチップ8の応答速度が向上することが示された。また、実施例3に係る発熱量測定システムを用いて、様々な流量で、測定対象混合ガスの発熱量を算出した。結果として、図21に示すように、流量を変化させても、発熱量の算出誤差は0.3以内に収まることが示された。したがって、実施例3に係る発熱量測定システムは、ガスタービンへのガスの供給路等、ガスの流量を変更することが好ましくない環境においても、正確に発熱量を算出可能であることが示された。その他の実施の形態上記のように、本発明は実施の形態によって記載したが、この開示の一部をなす記述及び図面はこの発明を限定するものであると理解するべきではない。この開示から当業者には様々な代替実施の形態、実施の形態及び運用技術が明らかになるはずである。例えば、実施の形態では、図8に示す式記憶装置402が、図1に示す第1の測温素子62からの電気信号及び複数の発熱温度における発熱素子61からの電気信号を独立変数とし、発熱量Qを従属変数とする発熱量算出式を保存する例を説明した。これに対し、上記22式で説明したように、気体の発熱量Qは、発熱素子61の温度がそれぞれTH1,TH2,TH3である場合の気体の放熱係数MI1TH1,MI2TH2,MI3TH3を変数とする方程式でも与えられる。よって、図8に示す式記憶装置402は、発熱素子61の複数の発熱温度における気体の放熱係数を独立変数とし、発熱量Qを従属変数とする発熱量算出式等の、放熱係数と、発熱量Qと、の相関関係を保存してもよい。この場合、測定部301は、パイプ101に注入された気体の放熱係数の測定値を、発熱素子61を複数の発熱温度に発熱させて測定する。なお、気体の放熱係数をマイクロチップ8を用いて測定可能であることは、上記9式で説明したとおりである。発熱量算出部305は、式記憶装置402に保存されている発熱量算出式の独立変数に、気体の放熱係数の測定値を代入して、気体の発熱量Qの測定値を算出する。次に、図22は、発熱抵抗体に2mA、2.5mA、及び3mAの電流を流した際の、混合ガスの放熱係数と、熱伝導率と、の関係を示す。図22に示すように、混合ガスの放熱係数と、熱伝導率と、は一般に比例関係にある。したがって、図8に示す式記憶装置402は、発熱素子61の複数の発熱温度における気体の熱伝導率を独立変数とし、発熱量Qを従属変数とする発熱量算出式を保存してもよい。この場合、測定部301は、パイプ101に注入された気体の熱伝導率の測定値を、発熱素子61を複数の発熱温度に発熱させて測定する。発熱量算出部305は、式記憶装置402に保存されている発熱量算出式の独立変数に、気体の熱伝導率の測定値を代入して、気体の発熱量Qの測定値を算出する。この様に、本発明はここでは記載していない様々な実施の形態等を包含するということを理解すべきである。"