JP6670706B2 - 発熱量測定装置および方法 - Google Patents
発熱量測定装置および方法 Download PDFInfo
- Publication number
- JP6670706B2 JP6670706B2 JP2016156190A JP2016156190A JP6670706B2 JP 6670706 B2 JP6670706 B2 JP 6670706B2 JP 2016156190 A JP2016156190 A JP 2016156190A JP 2016156190 A JP2016156190 A JP 2016156190A JP 6670706 B2 JP6670706 B2 JP 6670706B2
- Authority
- JP
- Japan
- Prior art keywords
- temperature
- index value
- value
- change rate
- heat
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 14
- 239000007789 gas Substances 0.000 claims description 83
- 230000008859 change Effects 0.000 claims description 64
- 239000002737 fuel gas Substances 0.000 claims description 43
- 238000004364 calculation method Methods 0.000 claims description 31
- 238000005259 measurement Methods 0.000 claims description 20
- 230000005540 biological transmission Effects 0.000 claims description 10
- 238000010438 heat treatment Methods 0.000 claims description 10
- 238000000691 measurement method Methods 0.000 claims description 5
- 238000012546 transfer Methods 0.000 claims description 5
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 16
- 238000010586 diagram Methods 0.000 description 8
- 239000003345 natural gas Substances 0.000 description 6
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 4
- ATUOYWHBWRKTHZ-UHFFFAOYSA-N Propane Chemical compound CCC ATUOYWHBWRKTHZ-UHFFFAOYSA-N 0.000 description 4
- 230000017525 heat dissipation Effects 0.000 description 4
- 238000000491 multivariate analysis Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- MYMOFIZGZYHOMD-UHFFFAOYSA-N Dioxygen Chemical compound O=O MYMOFIZGZYHOMD-UHFFFAOYSA-N 0.000 description 2
- OTMSDBZUPAUEDD-UHFFFAOYSA-N Ethane Chemical compound CC OTMSDBZUPAUEDD-UHFFFAOYSA-N 0.000 description 2
- 239000001273 butane Substances 0.000 description 2
- 239000001569 carbon dioxide Substances 0.000 description 2
- 229910002092 carbon dioxide Inorganic materials 0.000 description 2
- 229910001873 dinitrogen Inorganic materials 0.000 description 2
- 229910001882 dioxygen Inorganic materials 0.000 description 2
- IJDNQMDRQITEOD-UHFFFAOYSA-N n-butane Chemical compound CCCC IJDNQMDRQITEOD-UHFFFAOYSA-N 0.000 description 2
- OFBQJSOFQDEBGM-UHFFFAOYSA-N n-pentane Natural products CCCCC OFBQJSOFQDEBGM-UHFFFAOYSA-N 0.000 description 2
- 239000003209 petroleum derivative Substances 0.000 description 2
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000001294 propane Substances 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 239000002994 raw material Substances 0.000 description 2
- 238000007707 calorimetry Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000020169 heat generation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052697 platinum Inorganic materials 0.000 description 1
- 238000002076 thermal analysis method Methods 0.000 description 1
- 238000004861 thermometry Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/18—Investigating or analyzing materials by the use of thermal means by investigating thermal conductivity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/22—Fuels; Explosives
- G01N33/225—Gaseous fuels, e.g. natural gas
Landscapes
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Chemical & Material Sciences (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Description
本発明は、天然ガスを主原料とする都市ガスなどの発熱量を測定する発熱量測定装置および方法に関する。
例えば都市ガスは、天然ガス（ＮＧ:Natural Gas）を原料としてこれに石油ガス（ＰＧ:Petroleum Gas）を混ぜて所定の熱量（発熱量）に調整して供給されている。このような燃料ガスの熱量測定では、ガスを燃焼させて発熱量を計測する最も基本的な方法と、ガス密度と発熱量の関係を利用して発熱量を計測する方式がある。また、熱伝導率を測定することで、測定対象のガス発熱量を求める方法も提案されている（特許文献１参照）。
日本熱測定学会 編、「熱量測定・熱分析ハンドブック」、丸善出版株式会社、１０４頁、２０１０年発行。
しかしながら、上述した技術では、燃料ガスの組成比が大きく変化すると、熱伝導率の測定結果より求められるガス発熱量の精度が低下するという問題があった。
本発明は、以上のような問題点を解消するためになされたものであり、組成比が大きく変化してもより正確にガス発熱量が求められるようにすることを目的とする。
本発明に係る発熱量測定装置は、３つの異なる第１温度、第２温度、第３温度で、測定対象の燃料ガスの熱伝導率指標値および熱拡散率指標値を求める計測部と、計測部が計測した第１温度および第２温度の間の熱伝導率指標値の温度変化率κ1、第２温度および第３温度の間の熱伝導率指標値の温度変化率κ2、第１温度および第２温度の間の熱拡散率指標値の温度変化率α1、第２温度および第３温度の間の熱拡散率指標値の温度変化率α2を求める変化率算出部と、κ1、κ2、α1、α2を説明変数とし、発熱量を目的変数とする発熱量算出式により燃料ガスの発熱量を算出する発熱量算出部とを備え、発熱量算出式は、発熱量が既知の４つ以上の燃料ガスの各々について３つの異なる第１温度、第２温度、第３温度で熱伝導率指標値および熱拡散率指標値を求め、求めた各々について第１温度および第２温度の間の熱伝導率指標値の温度変化率κ1、第２温度および第３温度の間の熱伝導率指標値の温度変化率κ2、第１温度および第２温度の間の熱拡散率指標値の温度変化率α1、第２温度および第３温度の間の熱拡散率指標値の温度変化率α2を求め、求めた各々の燃料ガスにおけるκ1、κ2、α1、α2を説明変数とし、各々の既知の燃料ガスの発熱量を目的変数とする関係より得られた回帰式である。
上記発熱量測定装置において、計測部は、送熱部および受熱部を備え、送熱部からの電気信号および受熱部からの電気信号より熱伝導率指標値および熱拡散率指標値を求めればよい。
上記発熱量測定装置において、発熱量算出式は、成分ガスの組成比が異なる発熱量が既知の４つ以上の燃料ガスをサンプルガスとして用意し、用意した各々のサンプルガスの温度に依存する送熱部および受熱部からの電気信号の値を取得し、サンプルガスが接する送熱部を複数の発熱温度で発熱させ、各々の発熱温度において送熱部および受熱部からの電気信号の値を取得し、送熱部および受熱部からの電気信号の値から求めた熱伝導率指標値および熱拡散率指標値より各指標値の温度変化率κ1、κ2、α1、α2を求め、求めたκ1、κ2、α1、α2を説明変数とし、サンプルガスの発熱量を目的変数とすることで作成された回帰式である。
また、本発明に係る発熱量測定方法は、３つの異なる第１温度、第２温度、第３温度で、測定対象の燃料ガスの熱伝導率指標値および熱拡散率指標値を求める第１ステップと、第１温度および第２温度の間の熱伝導率指標値の温度変化率κ1、第２温度および第３温度の間の熱伝導率指標値の温度変化率κ2、第１温度および第２温度の間の熱拡散率指標値の温度変化率α1、第２温度および第３温度の間の熱拡散率指標値の温度変化率α2を求める第２ステップと、κ1、κ2、α1、α2を説明変数とし、発熱量を目的変数とする発熱量算出式により燃料ガスの発熱量を算出する第３ステップとを備え、発熱量算出式は、発熱量が既知の４つ以上の燃料ガスの各々について３つの異なる第１温度、第２温度、第３温度で熱伝導率指標値および熱拡散率指標値を求め、求めた各々について第１温度および第２温度の間の熱伝導率指標値の温度変化率κ1、第２温度および第３温度の間の熱伝導率指標値の温度変化率κ2、第１温度および第２温度の間の熱拡散率指標値の温度変化率α1、第２温度および第３温度の間の熱拡散率指標値の温度変化率α2を求め、求めた各々の燃料ガスにおけるκ1、κ2、α1、α2を説明変数とし、各々の既知の燃料ガスの発熱量を目的変数とする関係より得られた回帰式である。
上記発熱量測定方法において、熱伝導率指標値および熱拡散率指標値は、送熱部および受熱部を備える計測部を用いた燃料ガスの測定により、送熱部からの電気信号および受熱部からの電気信号により求めればよい。
上記発熱量測定方法において、成分ガスの組成比が異なる、発熱量が既知の４つ以上の燃料ガスをサンプルガスとして用意し、用意した各々のサンプルガスの温度に依存する送熱部および受熱部からの電気信号の値を取得し、サンプルガスが接する送熱部を複数の発熱温度で発熱させ、各々の発熱温度において送熱部および受熱部からの電気信号の値を取得し、送熱部および受熱部からの電気信号の値から求めた熱伝導率指標値および熱拡散率指標値より各指標値の温度変化率κ1、κ2、α1、α2を求め、求めたκ1、κ2、α1、α2を説明変数とし、サンプルガスの発熱量を目的変数として発熱量算出式を作成する。
以上説明したことにより、本発明によれば、組成比が大きく変化してもより正確にガス発熱量が求められるという優れた効果が得られる。
以下、本発明の実施の形態について図を参照して説明する。図１は、本発明の実施の形態における発熱量測定装置の構成を示す構成図である。この装置は、計測部１０１、変化率算出部１０２、発熱量算出部１０３、および記憶部１０４を備える。
計測部１０１は、３つの異なる第１温度、第２温度、第３温度で、測定対象の燃料ガスの熱伝導率指標値および熱拡散率指標値を求める。
変化率算出部１０２は、計測部１０１が計測した第１温度および第２温度の間の熱伝導率指標値の温度変化率κ1、第２温度および第３温度の間の熱伝導率指標値の温度変化率κ2、第１温度および第２温度の間の熱拡散率指標値の温度変化率α1、第２温度および第３温度の間の熱拡散率指標値の温度変化率α2を求める。
発熱量算出部１０３は、κ1、κ2、α1、α2を説明変数とし、発熱量を目的変数とする発熱量算出式により燃料ガスの発熱量を算出する。発熱量算出式は、記憶部１０４に記憶されている。
ここで、発熱量算出式は、発熱量が既知の４つ以上の燃料ガスの各々について、３つの異なる第１温度、第２温度、第３温度で熱伝導率指標値および熱拡散率指標値を求め、求めた各々について第１温度および第２温度の間の熱伝導率指標値の温度変化率κ1、第２温度および第３温度の間の熱伝導率指標値の温度変化率κ2、第１温度および第２温度の間の熱拡散率指標値の温度変化率α1、第２温度および第３温度の間の熱拡散率指標値の温度変化率α2を求め、求めた各々の燃料ガスにおけるκ1、κ2、α1、α2を説明変数とし、各々の既知の燃料ガスの発熱量を目的変数とする関係より得られた、以下に例示する回帰式である。なお、以下の式において、ｃiは係数である。また、次数は、例えば、ｎ１＋ｎ２＋ｎ３＋ｎ４＝０〜３とされていればよい。
次に、本発明の実施の形態における発熱量測定装置の動作例（発熱量測定方法）について、図２のフローチャートを用いて説明する。
まず、ステップＳ１０１で、計測部１０１が、第１温度Ｔ1における熱伝導率指標値ｋ1を取得する。次に、ステップＳ１０２で、計測部１０１が、第１温度Ｔ1における熱拡散率指標値ａ1を取得する。上述した計測部１０１における熱伝導率指標値および熱拡散率指標値の取得を、各々異なる第１温度Ｔ1，第２温度Ｔ2、第３温度Ｔ3で実施し（ステップＳ１０１〜ステップＳ１０３）、熱伝導率指標値ｋ1、熱伝導率指標値ｋ2、熱伝導率指標値ｋ3、熱拡散率指標値ａ1、熱拡散率指標値ａ2、熱拡散率指標値ａ3を得る。
次に、ステップＳ１０４で、変化率算出部１０２が、以下に示すように、第１温度Ｔ1および第２温度Ｔ2の間の熱伝導率指標値の温度変化率κ1、第２温度Ｔ2および第３温度Ｔ3の間の熱伝導率指標値の温度変化率κ2を求める。
次に、ステップＳ１０５で、変化率算出部１０２が、以下に示すように、第１温度Ｔ1および第２温度Ｔ2の間の熱拡散率指標値の温度変化率α1、第２温度Ｔ2および第３温度Ｔ3の間の熱拡散率指標値の温度変化率α2を求める。
次に、ステップＳ１０６で、発熱量算出部１０３が、κ1、κ2、α1、α2を説明変数とし、発熱量を目的変数とする発熱量算出式により燃料ガスの発熱量を算出する。前述したように、発熱量算出式は、発熱量が既知の４つ以上の燃料ガスの各々について、３つの異なる第１温度Ｔ1、第２温度Ｔ2、第３温度Ｔ3で熱伝導率指標値および熱拡散率指標値を求め、求めた各々について第１温度Ｔ1および第２温度Ｔ2の間の熱伝導率指標値の温度変化率κ1、第２温度Ｔ2および第３温度Ｔ3の間の熱伝導率指標値の温度変化率κ2、第１温度Ｔ1および第２温度Ｔ2の間の熱拡散率指標値の温度変化率α1、第２温度Ｔ2および第３温度Ｔ3の間の熱拡散率指標値の温度変化率α2を求め、求めた各々の燃料ガスにおけるκ1、κ2、α1、α2を説明変数とし、各々の既知の燃料ガスの熱量を目的変数とする関係より得られた回帰式である。
上述した実施の形態によれば、測定対象におけるガス組成が大きく変化しても、より正確にガス発熱量が求められるようになる。なお、実施の形態における発熱量測定装置は、ＣＰＵ（Central Processing Unit；中央演算処理装置）と主記憶装置と外部記憶装置とネットワーク接続装置となどを備えたコンピュータ機器であり、主記憶装置に展開されたプログラムによりＣＰＵが動作することで、上述した各機能が実現される。また、各機能は、複数のコンピュータ機器に分散させるようにしてもよい。
以下、より詳細に説明する。まず、単一成分のガスの発熱量は一意に特定できるため、ガス種（成分）およびこの割合が既知の場合、混合ガスの発熱量はガス組成から算出できる。ここで、天然ガスの主な各成分について、温度と熱伝導率指標値との関係を図３に示し、温度と熱拡散率指標値の関係を図４に示す。
ところで、特許文献１の技術では、発熱量の測定対象のガスにおける各ガス成分の割合を求めること無く、混合ガスの発熱量を求めている。この技術では、まず、複数の混合ガスの放熱係数または熱伝導率指標値の値を複数の温度で計測し、複数の混合ガスの既知の発熱量の値と、複数の温度で計測された放熱係数または熱伝導率指標値の値とに基づいて、複数の温度における放熱係数または熱伝導率指標値を説明変数とし、発熱量を目的変数とする発熱量算出式を作成しておく。
この発熱量算出式を用い、発熱量が未知の計測対象混合ガスの発熱量を測定する。まず、発熱量が未知の計測対象混合ガスの放熱係数または熱伝導率指標値の値を複数の温度で計測し、発熱量算出式の複数の温度における放熱係数または熱伝導率指標値の説明変数に、計測対象混合ガスの複数の温度で計測された放熱係数または熱伝導率指標値の値を代入し、計測対象混合ガスの発熱量の値を算出する。
しかしながら、各ガスの成分の割合（組成比）が大きく変動する場合、上述した従来の技術では、算出される発熱量の精度が得られていないという問題があった。
このような更なる測定精度の向上が求められているなかで、発明者らによる鋭意の検討の結果、温度が変化したときに得られる変化率を用いることで、高い精度で発熱量が求められることに想到した。温度変化による熱伝導率指標値および熱拡散率指標値の変化率が、測定精度の影響などを受けにくく、かつ発熱量推定の精度向上に有用であることを突き止めた。また、各測定値を説明変数とし、発熱量を目的変数として発熱量算出式を作成するので、真の熱伝導率および熱拡散率が測定される必要はなく、真の熱伝導率および熱拡散率が再現性よく反映された指標（値）が得られればよく、測定に対する制約が軽減されることも重要な特徴である。
次に、熱伝導率指標値および熱拡散率指標値の測定について説明する。これらの測定では、以下に示すように、送熱部および受熱部を備える計測部を用い、送熱部からの電気信号および受熱部からの電気信号により熱伝導率指標値および熱拡散率指標値を求める。
まず、熱伝導率指標値の測定について説明する。熱伝導率指標値の測定では、例えば、図５に示すように、送熱部２０１および受熱部２０２を用いて周囲のガス２１１の測定を行う。送熱部２０１および受熱部２０２は、例えば白金などより構成された測温抵抗体素子である。ガス２１１は、例えば、所定の処理装置の計測室内に収容され、気流および気圧などが安定している状態とする。
送熱部２０１、受熱部２０２、周辺のガス２１１は、温度Ｔ0で定常状態になっているものとする。この状態では、送熱部２０１および受熱部２０２の電気抵抗値は、温度Ｔ0相当値ｒbase＝ｒ0で一致している。
この状態より、例えば，ｔ1の時点で送熱部２０１の電気抵抗値が温度Ｔ1相当値ｒdrive＝ｒ1となるように送熱部２０１へ電流を印加すると、図６に示すように、受熱部２０２の温度（電気抵抗値）はｔ2の時点で定常状態となる。図６において、実線が送熱部２０１の電気抵抗値の変化を示し、点線が受熱部２０２の電気抵抗値の変化を示す。このとき、送熱部２０１と受熱部２０２間のガス２１１には、温度勾配が形成される。
特許文献１に記載されているように、送熱部２０１を温度Ｔ1となるように駆動した時の消費電力と、送熱部２０１と受熱部２０２との間の温度差から放熱係数を求めることができ、放熱係数は温度Ｔ1におけるガスの熱伝導率指標値と相関があることが知られている。従って、温度Ｔ１で送熱部２０１を駆動した時に得られる放熱係数から、上記相関関係を利用して導出される値を，温度Ｔ1における熱伝導率指標値ｋ1とする。
次に、熱拡散率指標値の測定について説明する。送熱部２０１、受熱部２０２、およびガス２１１は、温度Ｔ1で定常状態になっている場合、送熱部２０１と受熱部２０２の電気抵抗値は、温度Ｔ1相当値ｒbase＝ｒ1で一致している。
この状態より、例えばｔ1のタイミングで送熱部２０１の温度がΔＴだけ上昇するように電流を印加する（温度Ｔ1＋ΔＴ相当の電気抵抗値をｒdriveとする）。
送熱部２０１によってガス２１１を加熱して受熱部２０２でガス２１１の温度を観測すると、受熱部２０２で観測される温度上昇の過渡応答とガス２１１の熱拡散率指標値は相関があることは一般的に知られている（非特許文献１参照）。過渡応答を評価するパラメータとしては、受熱部２０２の温度上昇の過渡応答の時定数や、送熱部２０１が発熱開始してから受熱部２０２の温度がある一定値上昇するまでの時間、送熱部２０１が発熱開始してからある一定時間経過後の受熱部２０２の温度などが挙げられる。
従って、送熱部２０１、受熱部２０２、およびガス２１１が温度Ｔ１で定常状態となった後、送熱部２０１の温度をΔＴだけ上昇させたときの受熱部２０２の温度（電気抵抗値）の過渡応答のパラメータ（時定数など上記に列挙される指標）から相関関係を利用して導出される値を、温度Ｔ1における熱拡散率指標値ａ1とする。
次に、温度変化率について説明する。上述したように、温度Ｔ1における熱伝導率指標値ｋ1，熱拡散率指標値ａ1と同様にすることで、温度Ｔ2における熱伝導率指標値ｋ2，熱拡散率指標値ａ2、温度Ｔ3における熱伝導率指標値ｋ3，熱拡散率指標値ａ3を求める。このようにして求めた各値をもとに、前述した式（１）および式（２）よりκ1、κ2、α1、α2を求める。
次に、発熱量算出式の作成について説明する。発熱量算出式の作成では、例えば、４つ以上の発熱量既知のガスを用いて各々のκ1、κ2、α1、α2を求める。例えば、酸素ガス、窒素ガス、二酸化炭素ガス、メタンガス、エタンガス、プロパンガス、ブタンガスなどを成分として各成分の割合を変えた複数のサンプルガスを作製し、作製した複数のサンプルガスについて、各々のκ1、κ2、α1、α2を求める。上述した各サンプルガスについては、発熱量が既知である。
上述した各サンプルガスについて、組成（表１）と各組成に対する測定結果（表２）について以下に示す。表中の数値は体積％である。
各々のガスについて、κ1、κ2、α1、α2を説明変数とし、各々の既知の燃料ガスの熱量を目的変数とする複数の関係より、回帰式を多変量解析によって作成する。回帰式作成に用いる多変量解析の手法としては、例えば、サポートベクトル回帰が挙げられる。このようにして作成した回帰式を発熱量算出式として、記憶部１０４に記憶しておけばよい。
上述では、１０種類のガスを用い、κ1、κ2、α1、α2を説明変数とし、各々の既知の燃料ガスの熱量を目的変数とする１０個の関係より回帰式を得るようにしたが、これに限るものではない、より多くの種類のガスを対象として回帰式を得ることで、得られる回帰式（発熱量算出式）を用いて求める発熱量の統計的信頼性がより向上する。
実施の形態によれば、例えば特許文献１の技術に比較して、より高い精度で発熱量が求められるようになる。例えば、メタンガスの組成比が６８〜１００％，エタンガスの組成比が０〜２４％，プロパンガスの組成比が０〜６％、ブタンガスの組成比が０〜３％、酸素ガスの組成比が０〜０．２％、窒素ガスの組成比が０〜１５％、二酸化炭素ガスの組成比が０〜９％の範囲で変化する２１７２種類の混合ガス（燃料ガス）について、各々発熱量を求めることを検討する。
実施の形態においては、組成比が大きく変化している上述した２１７２種類の混合ガスのうち３４３種類の混合ガスをサンプルガスとし、各々のサンプルガスについて、κ1、κ2、α1、α2を説明変数とし、各々のサンプルガスの熱量を目的変数とする複数の関係より、多変量解析によって回帰式を作成する。
比較例として、特許文献１に示されているように、５つの異なる温度で、各燃料ガスの熱伝導率指標値ｋ1，ｋ2，ｋ3，ｋ4，ｋ5を求め、求めた熱伝導率指標値ｋ1，ｋ2，ｋ3，ｋ4，ｋ5を説明変数とし、各々のサンプルガスの熱量を目的変数とする複数の関係より、多変量解析によって回帰式を作製する。
上述した実施の形態による回帰式で上記２１７２種類の混合ガス（燃料ガス）について、各々発熱量を求めた結果における精度（推定誤差率）のばらつきは、図７に示すように、小さいものとなった。これに対し、比較例による回帰式で上記２１７２種類の混合ガス（燃料ガス）について、各々発熱量を求めた結果におけるばらつきは、図８に示すように、大きいものとなった。
以上に説明したように、本発明によれば、熱伝導率指標値の温度変化率および熱拡散率指標値の温度変化率を説明変数とし、燃料ガスの発熱量を目的変数とする発熱量算出式で発熱量を求めるようにしたので、組成比が大きく変化してもより正確にガス発熱量が求められるようになる。
なお、本発明は以上に説明した実施の形態に限定されるものではなく、本発明の技術的思想内で、当分野において通常の知識を有する者により、多くの変形および組み合わせが実施可能であることは明白である。
１０１…計測部、１０２…変化率算出部、１０３…発熱量算出部、１０４…記憶部。
Claims (6)
- ３つの異なる第１温度、第２温度、第３温度で、測定対象の燃料ガスの熱伝導率指標値および熱拡散率指標値を求める第１ステップと、
第１温度および第２温度の間の熱伝導率指標値の温度変化率κ1、第２温度および第３温度の間の熱伝導率指標値の温度変化率κ2、第１温度および第２温度の間の熱拡散率指標値の温度変化率α1、第２温度および第３温度の間の熱拡散率指標値の温度変化率α2を求める第２ステップと、
κ1、κ2、α1、α2を説明変数とし、発熱量を目的変数とする発熱量算出式により前記燃料ガスの発熱量を算出する第３ステップと
を備え、
前記発熱量算出式は、
発熱量が既知の４つ以上の燃料ガスの各々について３つの異なる第１温度、第２温度、第３温度で熱伝導率指標値および熱拡散率指標値を求め、求めた各々について第１温度および第２温度の間の熱伝導率指標値の温度変化率κ1、第２温度および第３温度の間の熱伝導率指標値の温度変化率κ2、第１温度および第２温度の間の熱拡散率指標値の温度変化率α1、第２温度および第３温度の間の熱拡散率指標値の温度変化率α2を求め、求めた各々の燃料ガスにおけるκ1、κ2、α1、α2を説明変数とし、各々の既知の燃料ガスの発熱量を目的変数とする関係より得られた回帰式である
ことを特徴とする発熱量測定方法。 - 請求項１記載の発熱量測定方法において、
前記熱伝導率指標値および前記熱拡散率指標値は、送熱部および受熱部を備える計測部を用いた前記燃料ガスの測定により、前記送熱部からの電気信号および前記受熱部からの電気信号より求める
ことを特徴とする発熱量測定方法。 - 請求項２記載の発熱量測定方法において、
成分ガスの組成比が異なる、発熱量が既知の４つ以上の燃料ガスをサンプルガスとして用意し、
用意した各々の前記サンプルガスの温度に依存する前記送熱部および前記受熱部からの電気信号の値を取得し、
前記サンプルガスが接する前記送熱部を複数の発熱温度で発熱させ、
前記各々の発熱温度において前記送熱部および前記受熱部からの電気信号の値を取得し、
前記送熱部および前記受熱部からの電気信号の値から求めた熱伝導率指標値および熱拡散率指標値より各指標値の温度変化率κ1、κ2、α1、α2を求め、
求めたκ1、κ2、α1、α2を説明変数とし、前記サンプルガスの発熱量を目的変数として前記発熱量算出式を作成する
ことを特徴とする発熱量測定方法。 - ３つの異なる第１温度、第２温度、第３温度で、測定対象の燃料ガスの熱伝導率指標値および熱拡散率指標値を求める計測部と、
前記計測部が計測した第１温度および第２温度の間の熱伝導率指標値の温度変化率κ1、第２温度および第３温度の間の熱伝導率指標値の温度変化率κ2、第１温度および第２温度の間の熱拡散率指標値の温度変化率α1、第２温度および第３温度の間の熱拡散率指標値の温度変化率α2を求める変化率算出部と、
κ1、κ2、α1、α2を説明変数とし、発熱量を目的変数とする発熱量算出式により前記燃料ガスの発熱量を算出する発熱量算出部と
を備え、
前記発熱量算出式は、
発熱量が既知の４つ以上の燃料ガスの各々について３つの異なる第１温度、第２温度、第３温度で熱伝導率指標値および熱拡散率指標値を求め、求めた各々について第１温度および第２温度の間の熱伝導率指標値の温度変化率κ1、第２温度および第３温度の間の熱伝導率指標値の温度変化率κ2、第１温度および第２温度の間の熱拡散率指標値の温度変化率α1、第２温度および第３温度の間の熱拡散率指標値の温度変化率α2を求め、求めた各々の燃料ガスにおけるκ1、κ2、α1、α2を説明変数とし、各々の既知の燃料ガスの発熱量を目的変数とする関係より得られた回帰式である
ことを特徴とする発熱量測定装置。 - 請求項４記載の発熱量測定装置において、
前記計測部は、送熱部および受熱部を備え、前記送熱部からの電気信号および前記受熱部からの電気信号より前記熱伝導率指標値および熱拡散率指標値を求める
ことを特徴とする発熱量測定装置。 - 請求項５記載の発熱量測定装置において、
前記発熱量算出式は、
成分ガスの組成比が異なる発熱量が既知の４つ以上の燃料ガスをサンプルガスとして用意し、
用意した各々の前記サンプルガスの温度に依存する前記送熱部および前記受熱部からの電気信号の値を取得し、
前記サンプルガスが接する前記送熱部を複数の発熱温度で発熱させ、
前記各々の発熱温度において前記送熱部および前記受熱部からの電気信号の値を取得し、
前記送熱部および前記受熱部からの電気信号の値から求めた熱伝導率指標値および熱拡散率指標値より各指標値の温度変化率κ1、κ2、α1、α2を求め、
求めたκ1、κ2、α1、α2を説明変数とし、前記サンプルガスの発熱量を目的変数とすることで作成された回帰式である
ことを特徴とする発熱量測定装置。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016156190A JP6670706B2 (ja) | 2016-08-09 | 2016-08-09 | 発熱量測定装置および方法 |
| CN201780046516.1A CN109564178B (zh) | 2016-08-09 | 2017-07-24 | 发热量测定装置以及方法 |
| US16/321,652 US11029270B2 (en) | 2016-08-09 | 2017-07-24 | Apparatus and method for measuring calorific value |
| EP17839200.7A EP3499223A4 (en) | 2016-08-09 | 2017-07-24 | DEVICE AND METHOD FOR CALORIE VALUE MEASUREMENT |
| PCT/JP2017/026611 WO2018030127A1 (ja) | 2016-08-09 | 2017-07-24 | 発熱量測定装置および方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016156190A JP6670706B2 (ja) | 2016-08-09 | 2016-08-09 | 発熱量測定装置および方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2018025420A JP2018025420A (ja) | 2018-02-15 |
| JP6670706B2 true JP6670706B2 (ja) | 2020-03-25 |
Family
ID=61161967
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2016156190A Expired - Fee Related JP6670706B2 (ja) | 2016-08-09 | 2016-08-09 | 発熱量測定装置および方法 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US11029270B2 (ja) |
| EP (1) | EP3499223A4 (ja) |
| JP (1) | JP6670706B2 (ja) |
| CN (1) | CN109564178B (ja) |
| WO (1) | WO2018030127A1 (ja) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2020148684A (ja) * | 2019-03-14 | 2020-09-17 | オムロン株式会社 | 流量測定装置、流量測定装置を備えたガスメータ及び、ガスメータのための流量測定装置ユニット |
| WO2021256356A1 (ja) * | 2020-06-15 | 2021-12-23 | 理研計器株式会社 | 熱量測定装置および熱量算出方法 |
| CN113848231B (zh) * | 2020-06-28 | 2024-03-08 | 宝山钢铁股份有限公司 | 基于炼焦煤热解过程中热扩散率的结焦性判断方法 |
| CN115790908B (zh) * | 2023-02-08 | 2023-06-23 | 成都千嘉科技股份有限公司 | 基于热量计量的天然气计量方法和装置 |
Family Cites Families (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU3055992A (en) * | 1991-10-23 | 1993-05-21 | Niagara Mohawk Power Corporation | On-line combustionless measurement of gaseous fuels fed to gas consumption devices |
| EP0554095A3 (en) * | 1992-01-30 | 1994-12-14 | Honeywell Inc | Determination of fuel characteristics |
| CA2126481C (en) * | 1994-06-22 | 2001-03-27 | Andreas Mandelis | Non-contact photothermal method for measuring thermal diffusivity and electronic defect properties of solids |
| WO2007036983A1 (ja) * | 2005-09-27 | 2007-04-05 | Yamatake Corporation | 熱伝導率測定方法および装置、並びにガス成分比率測定装置 |
| US20110185789A1 (en) | 2008-10-01 | 2011-08-04 | Yamatake Corporation | Calorific value calculation formula generating system, calorific value calculation formula generating method, calorific value calculating system, and calorific value calculating method |
| JP5344958B2 (ja) * | 2009-03-12 | 2013-11-20 | アズビル株式会社 | 発熱量算出式作成システム、発熱量算出式の作成方法、発熱量算出システム、及び発熱量の算出方法 |
| JP5420456B2 (ja) * | 2010-03-17 | 2014-02-19 | アズビル株式会社 | 発熱量算出式作成システム、発熱量算出式の作成方法、発熱量測定システム、及び発熱量の測定方法 |
| JP5421832B2 (ja) * | 2010-03-25 | 2014-02-19 | アズビル株式会社 | 発熱量算出式作成システム、発熱量算出式の作成方法、発熱量測定システム、及び発熱量の測定方法 |
| JP5335722B2 (ja) * | 2010-03-26 | 2013-11-06 | アズビル株式会社 | 発熱量算出式作成システム、発熱量算出式の作成方法、発熱量測定システム、及び発熱量の測定方法 |
| JP5335727B2 (ja) * | 2010-03-29 | 2013-11-06 | アズビル株式会社 | 発熱量算出式作成システム、発熱量算出式の作成方法、発熱量測定システム、及び発熱量の測定方法 |
| JP5335728B2 (ja) * | 2010-03-29 | 2013-11-06 | アズビル株式会社 | 発熱量算出式作成システム、発熱量算出式の作成方法、発熱量測定システム、及び発熱量の測定方法 |
| JP5534193B2 (ja) * | 2010-04-20 | 2014-06-25 | アズビル株式会社 | 温度拡散率測定システム及び流量測定システム |
| JP5759780B2 (ja) * | 2011-05-09 | 2015-08-05 | アズビル株式会社 | 発熱量測定システム及び発熱量の測定方法 |
| JP5781968B2 (ja) * | 2012-03-27 | 2015-09-24 | アズビル株式会社 | 発熱量測定システム及び発熱量の測定方法 |
| US9354220B2 (en) * | 2013-09-27 | 2016-05-31 | Caterpillar Inc. | Engine system having fuel quality sensor |
| DE102014115566A1 (de) * | 2013-11-19 | 2015-05-21 | Endress + Hauser Flowtec Ag | Messgerät und Verfahren zur Bestimmung eines korrigierten Massedurchflusses und Verwendungen des Messgerätes |
| JP2016057198A (ja) * | 2014-09-10 | 2016-04-21 | セイコーエプソン株式会社 | 熱流計測装置及び代謝計測装置 |
| CN105606786B (zh) * | 2014-11-14 | 2019-12-24 | Mems股份公司 | 用于确定燃气品质的特定量值的方法和测量装置 |
| JP5813260B1 (ja) * | 2015-02-25 | 2015-11-17 | トーソー株式会社 | コード結合具 |
| CN105807027B (zh) * | 2016-03-16 | 2018-11-09 | 新奥科技发展有限公司 | 气体能量计量方法及装置 |
-
2016
- 2016-08-09 JP JP2016156190A patent/JP6670706B2/ja not_active Expired - Fee Related
-
2017
- 2017-07-24 WO PCT/JP2017/026611 patent/WO2018030127A1/ja unknown
- 2017-07-24 CN CN201780046516.1A patent/CN109564178B/zh active Active
- 2017-07-24 EP EP17839200.7A patent/EP3499223A4/en not_active Ceased
- 2017-07-24 US US16/321,652 patent/US11029270B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| EP3499223A4 (en) | 2020-03-04 |
| US20200049636A1 (en) | 2020-02-13 |
| US11029270B2 (en) | 2021-06-08 |
| CN109564178B (zh) | 2021-10-22 |
| EP3499223A1 (en) | 2019-06-19 |
| WO2018030127A1 (ja) | 2018-02-15 |
| CN109564178A (zh) | 2019-04-02 |
| JP2018025420A (ja) | 2018-02-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6670706B2 (ja) | 発熱量測定装置および方法 | |
| Kubicár et al. | A step-wise method for measuring thermophysical parameters of materials | |
| CN101464422B (zh) | 一种固体材料的导热系数测量仪 | |
| RU2343466C1 (ru) | Способ определения теплопроводности материалов | |
| CN105973504B (zh) | 一种77k～90k温区热电偶标定系统及方法 | |
| Malinarič et al. | Comparison of the new plane source method to the step wise transient method for thermal conductivity and diffusivity measurement | |
| Haloua et al. | Uncertainty analysis of theoretical methods for adiabatic temperature rise determination in calorimetry | |
| Beckmüller et al. | A fundamental equation of state for the calculation of thermodynamic properties of n-octane | |
| Sairanen et al. | Validation of a calibration set-up for radiosondes to fulfil GRUAN requirements | |
| CN101464423B (zh) | 一种固体材料的导热系数的测量装置 | |
| JPH0875688A (ja) | 熱量計 | |
| Haloua et al. | Traceable measurement and uncertainty analysis of the gross calorific value of methane determined by isoperibolic calorimetry | |
| JP2017156278A (ja) | 流速の評価方法、及び評価システム | |
| Boersma et al. | Gas composition sensor for natural gas and biogas | |
| JP2017026385A (ja) | 熱伝導率測定装置、熱伝導率測定方法、及び真空度評価装置 | |
| JP4504320B2 (ja) | 熱量測定方法 | |
| Pokhodun | Current status and prospects for development of thermometry | |
| JP2017194323A (ja) | 熱伝導率計 | |
| RU2755330C1 (ru) | Способ измерения теплопроводности | |
| Scandurra et al. | Fast Johnson noise thermometry using a temperature dependent sensor | |
| Kubičár et al. | „Transient methods for the measurement of thermophysical properties “ | |
| Berberian | Note: Low temperature calibration of a commercial differential scanning calorimeter | |
| JP2017142130A (ja) | 窒素濃度測定システム及び窒素濃度の測定方法 | |
| RU2550306C1 (ru) | Способ измерения объемной концентрации водорода | |
| Yang et al. | Modelling the calibration of the industrial platinum-resistance thermometers according to GUM |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190311 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20200218 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200302 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6670706 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |