JP7258757B2 - 水素を含みうるガスの燃焼特性を推定するための方法 - Google Patents

水素を含みうるガスの燃焼特性を推定するための方法 Download PDF

Info

Publication number
JP7258757B2
JP7258757B2 JP2019535340A JP2019535340A JP7258757B2 JP 7258757 B2 JP7258757 B2 JP 7258757B2 JP 2019535340 A JP2019535340 A JP 2019535340A JP 2019535340 A JP2019535340 A JP 2019535340A JP 7258757 B2 JP7258757 B2 JP 7258757B2
Authority
JP
Japan
Prior art keywords
fuel gas
measured
gas
fuel
module
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.)
Active
Application number
JP2019535340A
Other languages
English (en)
Other versions
JP2020514714A (ja
Inventor
ウルリアック，マテュー
カペラ，サンドラ
ラントワン，ローラン
マンジュー，ナウシャド
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Engie SA
Original Assignee
Engie SA
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Engie SA filed Critical Engie SA
Publication of JP2020514714A publication Critical patent/JP2020514714A/ja
Application granted granted Critical
Publication of JP7258757B2 publication Critical patent/JP7258757B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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
    • G01N25/22—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 on combustion or catalytic oxidation, e.g. of components of gas mixtures
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10—Analysis or design of chemical reactions, syntheses or processes
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Physics & Mathematics (AREA)
  • General Chemical & Material Sciences (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Oil, Petroleum & Natural Gas (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Combustion & Propulsion (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)
  • Regulation And Control Of Combustion (AREA)
  • Feeding And Controlling Fuel (AREA)
  • Measuring Volume Flow (AREA)

Description

本発明は、燃料ガスのグループの一部を形成する燃料ガスの燃焼特性を推定する一般分野に関し、少なくとも１つの特性は、ウォッベ指数又は高位発熱量である。
本発明の特定の用途は、水素を含有する燃料ガス、例えば、体積で０％から２０％の範囲にある水素の量を含有する燃料ガスの燃焼特性を推定することにある。
再生可能資源は、断続的かつエネルギーの必要量と相関関係がないように使用することがあるため、再生可能エネルギー源を使用する際には、エネルギーを貯蔵するための手段の使用が必要となることがある。このことは、特に、太陽光発電、又は、風力タービンを用いた発電の場合に適用される。
再生可能エネルギー源からのエネルギーを貯蔵するための手段として、水の電気分解によって得られる水素を使用することが提案されている。
この方法で生成される水素は、次いで貯蔵庫、例えば、地下に貯蔵でき、又は、天然ガス供給網に導入することもできる。得られるガスは、したがって、０％から２０％の範囲で水素を含みうる。
比重の平方根に対する燃料ガスの発熱量の比率であるウォッベ指数（キロワット時毎ノルマル立方メートル（ｋＷｈ／Ｎｍ3））、及び高位発熱量（ｋＷｈ／Ｎｍ3）は、家庭的利用であろうと産業的利用であろうと、ガスバーナーの燃焼の設定に影響を与える２つの主な燃料の特徴的な数値である。
ガスの用途の種類、及び使用される測定手段（流量が吸引発生手段によって測定されるかどうか）に応じて、これら２つのパラメータのうちの一方又は他方が考慮に入れられる。
一例として、一定の供給圧力で供給されるインジェクタ（例えば、クックトップ又はリボンバーナー）の場合、供給される出力は、圧力、インジェクタの貫通部、水頭損失、及び燃料のウォッベ指数に依存する。
質量又は体積の流量計が利用できる適用例では、供給される出力は、測定された流速及び発熱量に依存する。
特に、バーナーによって供給される熱量は、ウォッベ指数又は高位発熱量（ＰＣＳ）のいずれかに比例する。同様に、ヨーロッパで流通している、水素を含まない天然ガスの場合、化学量論での燃焼を得るために必要な空気の量は、燃料ガスのこれらの特性のうちの一方又は他方に依存する。
天然ガスの輸送及び流通のためのヨーロッパのネットワークは、ますます網目状になり、様々な供給元によって供給されている。そのため、天然ガスの特性（中でもウォッベ指数及びＰＣＳ）は、したがって、ネットワークの任意の点で、無視できない程度で（±５％以上）時間とともに変化しうることに留意するべきである。
そのため、ガラス製造、セラミック、発電、石灰、及び冶金産業で行われるような産業的過程は、これらの変化に対して敏感である。結果として、燃焼を最適化するために、燃焼を調節するための具体的な解決策を利用することが必要である。この種類の調節を実行するために、ウォッベ指数又は高位発熱量を測定することが可能である。
天然ガスの場合、単純に、堅牢に、そして正確にこれらのパラメータを測定するための技術的な解決策がある。しかし、水素が燃料中に存在するときには、そのような解決策は不十分である。水素の体積が１％を超えると、これらのパラメータの測定における誤差が大幅に増加し、測定が再現可能となる保証はない。
さらに、産業用の機器、さらには家庭用の機器でさえ、非常に小さい割合であっても、水素の添加に対して敏感でありうる。
結果として、ガラス製造又は冶金工程のような高感度での使用だけでなく、ネットワークの操作者が利用者に供給されるガス燃料の燃焼特性を局所的に管理できるように、これらの新規な燃料の主要な燃焼パラメータを継続的な方法で定量できる新規な測定装置が必要である。
ウォッベ指数及び発熱量を測定又は推定するための方法が提案されている。
具体的に、おそらくゼロでない割合の水素を含む天然ガスのウォッベ指数及び／又は発熱量を測定するために使用されうる装置の４つのグループがある。
・熱量測定（発熱量用）：測定は、制御された量のガスを燃焼させることを含む。その後、放出されるエネルギーは、（ガス又は水の）熱交換器を横切る入力／出力の温度差によって定量される。精度は１０％のオーダー（すなわち、１．１ｋＷｈ／Ｎｍ3）である。しかし、各測定は約１０分かかる。その技術は、次第に使われなくなっている。
・燃焼測定（ウォッベ指数用）：測定は、空気／ガスの混合物を燃焼させることを含む。燃焼生成物における酸素の含有量は、ジルコニアプローブによって測定される。残留している酸素の含有量は、燃焼性指数と相関しており、それ自体は、天然ガスに関して、ウォッベ指数と相関している。この種類の装置は、約１５，０００ユーロかかる。この種類の機器は、水素の含有量が低い混合物について作動する（体積で５％を超える水素の含有量に対しては、ウォッベ指数は、もはや燃焼性指数と直線的に比例しない）。この技術の欠点は、５％のオーダー（すなわち、０．７５ｋＷｈ／Ｎｍ3）と精度が低いこと、及び相当なメンテナンスを要すること（混合物の燃焼が起こるオーブンの高温による劣化）である。
・気相クロマトグラフィー（ウォッベ指数用及び発熱量用）：この技術は、気体混合物から分子を分離するのに役立つ。ウォッベ指数及び発熱量は、ガスの組成に関する情報に基づいて計算されうる。高出力での燃料の燃焼機器（例えば、ガスタービン）の場合、気相クロマトグラフィーは、一般に、熱量測定及び燃焼測定に取って代わる。これらの測定装置の精度はよりよく、０．５％のオーダー（すなわち、ウォッベ指数に関して０．０８ｋＷｈ／Ｎｍ3）である。機器が水素分子を分離できる特定のセンサを備えているならば、この種類の機器は天然ガスと水素との混合物についても作動しうる。気相クロマトグラフィーの主な欠点は、昨今の改善が考慮されたとしても、費用である（費用はおそらく２０，０００ユーロから５０，０００ユーロの範囲にありうる）。さらに、最善の状況でも、気相クロマトグラフィー機器の応答時間は１分のオーダーである。したがって、測定されるガスと燃料の燃焼機器（例えば、バーナー）によって実際に使用されているガスとの間で不一致がありうる。
・相関装置（ウォッベ指数用又は高位発熱量用）：ウォッベ指数又は高位発熱量と相関する１つ又は複数の物理的な数値がそのような装置で測定される。コンピュータによって実行される相関関係は、次に、ウォッベ指数又はＰＣＳを推定するのに使用される。この種類の装置は、ヨーロッパで流通している種類の天然ガスに関して、１％以下の精度を得ることを可能にする。この種類の技術の利点は、応答速度（瞬時）、費用（１０，０００ユーロから２０，０００ユーロの範囲）、堅牢性、及びメンテナンスの削減である。それにもかかわらず、市販の装置は、天然ガスと水素との混合物について作動しない。実際の天然ガス（純粋なメタンではない）が体積で水素を数パーセント含むと、誤差は劇的に増加する。
最新技術に含まれる文献EP 1 346 215には、相関関係によってウォッベ指数を測定するための装置が記載されている。
文献US 4 384 792も知られており、これには、相関関係によってガス燃料のウォッベ指数を測定し調節するための装置、及びウォッベ指数を調節することが記載されている。
文献US 6 244 097も知られており、これには、相関関係によってガス燃料の発熱量を測定するための装置が開示されている。
最後に、既知の文献DE 4 118 781には、相関関係によってガス燃料の発熱量及びウォッベ指数を測定するための装置が開示されている。
それらの文献に開示されている解決策は、相関装置の欠点を有する。特に、それらの精度は非常に低いので、調査対象の燃料ガスが水素を含むと、それらは使用に適さない。
さらに、上記したように、天然ガスと水素との混合物について作動できる他の装置は、特に高価であり、このことは、それらの用途を、主要な産業グループに設置されうる種類の高出力の機器に限定する。
本発明は、特に、これらの欠点の少なくともいくつかを軽減することを目的とする。
本発明は、燃料ガスのグループに属する燃料ガスの少なくとも１つの燃焼特性を推定する方法を提案することによってこの要求を満たし、前記少なくとも１つの特性は、ウォッベ指数又は高位発熱量であり、前記方法は、
・前記燃料ガスの少なくとも２つの流動特性を測定することと、
・前記燃料ガスに含まれている水素の含有量ＸH2を測定することと、
・前記少なくとも１つの特性ΞGN/H2は、次の経験的なアフィン関係を用いて推定されることと、
を含む。
Figure 0007258757000001
ここで、
α，β，及びγは、前記燃料ガスのグループに関する所定の係数であり、
Ｙは、前記燃料ガスの前記少なくとも２つの流動特性の測定値から生成される前記燃料ガスの物理特性を表す変数である。
この方法は、したがって、相関型の方法に関連し、この方法は、したがって、簡単かつ安価に実施されうる。この方法により、結果を極めて迅速に得ることが可能となる。
本発明は、従来技術の方法とは異なり、燃料ガスにおける水素の含有量を考慮に入れている。
本発明は、０％から２０％の範囲の水素の体積含有量を有する燃料ガスに適用できることに留意するべきである。
本発明者らは、高位発熱量及びウォッベ指数の両方を推定するために同じ形式の経験的なアフィン関係を使用することが可能であることを見出した。
参考として、測定されているのが高位発熱量であるかウォッベ指数であるかに応じて、燃料ガスの様々な異なる流動特性が測定されうる。言い換えれば、変数Ｙは、測定されているのが高位発熱量であるかウォッベ指数であるかに応じて異なる形式を有していてもよい。
さらに、係数α，β，及びγは、定数である必要はなく、それらは、測定されているのが高位発熱量であるかウォッベ指数であるかに応じて、異なっていてもよい。
燃料ガスの流動特性は、流動手段において測定された特性であること、及び、それらは、流量、温度、又は圧力を含む群から選択されてもよいことに留意するべきである。
変数Ｙは、ガスの物理特性を表す。参考として、変数Ｙは、燃料ガスの粘度、その比熱容量（単位重量当たりの）、又はその密度などの１つ又は複数の物理特性と関連付けられていてもよい。変数Ｙがこれらの物理特性を表し、又は変数Ｙがこれらの物理特性に関連付けられていることは、これらの物理特性を関連付けている数学関数の形式でそれを記述できることを意味する。これらの特性は変数Ｙを生成するために使用されるので、変数Ｙは、物理特性を関連付けている数学関数の形式、及び、燃料ガスの流動特性を関連付けている数学関数の形式の両方で記述されていてもよい。特に、流動特性の測定の特定の種類に関して、複数のそのような測定を燃料ガスの流れの物理特性と関連付けることが可能である。
したがって、
流動特性Ｕ1の測定に関してＵmes,1の記号を使用し、
流動特性Ｕ2の測定に関してＵmes,2の記号を使用し、
流動特性Ｕ3の測定に関してＵmes,3の記号を使用することによって、以下が与えられる。
Figure 0007258757000002
特定の実施では、前記係数α，β，及びγは、入力として前記水素の含有量ＸH2に関する前記測定値を有し、出力として前記係数α，β，及びγを提供する表から読み取られた係数である。
これにより、推定されている特性のより正確な推定を行うことが可能になる。
特定の実施では、前記表は、前記係数α，β，及びγを、前記水素の含有量ＸH2に関して１％の幅を有する複数の数値範囲と関連付ける。
言い換えれば、表は、１％のオーダーの粒度を有する。
本発明の発明者らは、これにより、許容できるレベルの複雑さで満足できる精度が得られることを見出した。
特定の実施では、前記係数α，β，及びγの値は、前記物理特性を表す前記Ｙの値及び前記燃焼特性が知られている、前記燃料ガスのグループの既知のガスに関するデータセットから得られる。
「既知のガス」の語句は、組成が知られているガス、例えば、決定された化学的性質を持つ様々な成分の割合が知られているガスを包含するべく使用される。当業者は、そのようなガスの燃焼特性を決定する方法を知っており、特に、ゼロでない量の水素を含む既知のガスを利用できる。
特定の実施では、前記燃料ガスのグループの既知のガスに関する前記データセットから、燃焼特性と物理特性を表すＹに関する値とをランダムに生成することを含む。
特定の実施では、前記ウォッベ指数及び前記高位発熱量は、２つの経験的なアフィン関係を用いて推定される。
したがって、この方法を１回実施すれば、これらの燃焼特性の両方を推定することが可能である。
特定の実施では、前記方法は、前記推定されたウォッベ指数及び前記推定された高位発熱量から前記燃料ガスの密度を推定することをさらに含む。
特定の実施では、前記方法は、前記燃料ガスの前記燃焼特性を調節すること、又は、燃料ガスの前記燃焼特性と、推定された空気の化学量論的な体積もしくは前記推定された特性に対応する推定された燃焼性指数とを調節することを含む。
特定の実施では、前記燃料ガスの前記少なくとも１つの燃焼特性は、前記ウォッベ指数ＩＷGN/H2を含み、前記燃料ガスの少なくとも２つの流動特性の前記測定は、流体狭窄部（例えば、オリフィス又はマイクロノズル）を通る音速流の状態にある（すなわち、音速以上の速度で）前記燃料ガスの質量流量を測定することを含み、前記測定は、前記狭窄部の上流で測定された絶対圧力及び前記狭窄部の上流で測定された絶対温度で行われ、
前記方法は、前記狭窄部の上流で測定された基準ガスの絶対圧力及び前記狭窄部の上流で測定された基準ガスの絶対温度で、前記流体狭窄部を通る音速流の状態にある基準ガス（例えば、メタン）の質量流量の測定が行われる校正手順をさらに含み、
前記ウォッベ指数ＩＷGN/H2を推定するために使用される前記経験的なアフィン関係は、次の式で記述される。
Figure 0007258757000003
Figure 0007258757000004
ここで、
Ｑmes,2は、測定された前記燃料ガスの前記質量流量であり、
Ｐmesは、測定された前記燃料ガスの前記絶対圧力であり、
Ｔmesは、測定された前記燃料ガスの前記絶対温度であり、
Ｑrefは、測定された前記基準ガスの前記質量流量であり、
Ｐrefは、測定された前記基準ガスの前記絶対圧力であり、
Ｔrefは、測定された前記基準ガスの前記絶対温度であり、
Ｄ，Ｅ，及びＦは、前記燃料ガスのグループに関する所定の係数であり、それぞれ、係数α，β，及びγに対応する。
特定の実施では、前記燃料ガスの前記ウォッベ指数が推定され、前記方法は、前記燃料ガスの前記密度を測定すること、及び前記推定されたウォッベ指数と測定された前記ガスの密度とから前記高位発熱量を推定すること、をさらに含む。
特定の実施では、前記燃料ガスの前記少なくとも１つの燃焼特性は、前記高位発熱量ＰＣＳGN/H2を含み、前記燃料ガスの少なくとも２つの流動特性の前記測定は、
・圧力降下を生じさせる装置を通る層流の状態にある前記燃料ガスの前記質量流量を測定し、前記測定は、前記燃料ガスの粘度及び基準ガス（例えば、メタン）の粘度に依存すること、
・圧力降下を生じさせる前記装置の下流で、熱式質量流量計を用いて前記燃料ガスの前記質量流量を測定し、前記測定は、前記燃料ガスの比熱容量及び基準ガスの熱容量に依存すること、
を含み、
前記高位発熱量ＰＣＳGN/H2を推定するために使用される前記経験的なアフィン関係は、次の式で記述される。
Figure 0007258757000005
Figure 0007258757000006
ここで、
Ｚは、変数Ｙに対応する変数であり、
Ｑmes,1は、測定された圧力降下を生じさせる装置を通る層流の状態にある前記燃料ガスの前記質量流量であり、
Ｑmes,2は、圧力降下を生じさせる前記装置の下流で測定された前記燃料ガスの前記質量流量であり、
Ａ，Ｂ，及びＣは、前記燃料ガスのグループに関する所定の係数であり、それぞれ、前記係数α，β，及びγに対応する。
例えば、基準ガスの粘度に関する値及び基準ガスの熱容量に関する値を得るために、基準ガスを用いて校正工程を実行することが可能である。
この校正工程は、特に、水素を含有するガスのような、純粋ではないガスに関して、正確な測定を行うことを可能にする。
本発明は、また、燃料ガスのグループに属する燃料ガスの少なくとも１つの燃焼特性を推定するための装置を提供し、前記少なくとも１つの特性は、ウォッベ指数又は高位発熱量であり、前記装置は、
・前記燃料ガスの少なくとも２つの流動特性を測定するための少なくとも２つのモジュールと、
・前記燃料ガスに含まれている水素の含有量ＸH2を測定するためのモジュールと、
・次の経験的なアフィン関係を用いて前記少なくとも１つの特性ΞGN/H2を推定するように構成されたモジュールと、
を含む。
Figure 0007258757000007
ここで、
α，β，及びγは、前記燃料ガスのグループに関する所定の係数であり、
Ｙは、前記燃料ガスの前記少なくとも２つの流動特性の測定値から生成される前記燃料ガスの物理特性を表す変数である。
この装置は、上記したように、本方法の実施の全てを実行するように構成されていてもよい。
この装置は、燃料ガスの流動特性が測定されうるように、燃料ガスが流れる流動部材を含んでいてもよいことに留意するべきである。
特定の実施形態では、前記燃料ガスの前記少なくとも１つの燃焼特性は、前記ウォッベ指数ＩＷGN/H2を含み、前記装置は、
・前記燃料ガスの流れを受け入れるための入口、
・基準ガスの流れを受け入れるための入口、
・前記燃料ガスの前記流れ又は前記基準ガスの前記流れを管に導くためのセレクタ及び案内モジュール、
・流体狭窄部、及び、
・前記流体狭窄部を通る音速流の状態にある前記燃料ガスの質量流量を測定するためのモジュールであって、前記狭窄部の上流で絶対圧力を測定するためのサブモジュール、及び前記狭窄部の上流で絶対温度を測定するためのサブモジュール、を含むモジュール、
を備え、
前記ウォッベ指数ＩＷGN/H2を推定するために使用される前記経験的なアフィン関係は、次の式で記述される。
Figure 0007258757000008
Figure 0007258757000009
ここで、
Ｑmes,2は、測定された前記燃料ガスの前記質量流量であり、
Ｐmesは、測定された前記燃料ガスの前記絶対圧力であり、
Ｔmesは、測定された前記燃料ガスの前記絶対温度であり、
Ｑrefは、測定された前記基準ガスの前記質量流量であり、
Ｐrefは、測定された前記基準ガスの前記絶対圧力であり、
Ｔrefは、測定された前記基準ガスの前記絶対温度であり、
Ｄ，Ｅ，及びＦは、前記燃料ガスのグループに関する所定の係数であり、それぞれ、前記係数α，β，及びγに対応する。
特定の実施形態では、前記装置は、前記燃料ガスの前記ウォッベ指数を推定するのに適しており、前記装置は、前記燃料ガスの密度を測定するためのモジュールをさらに備え、前記少なくとも１つの特性を推定するために構成されている前記モジュールは、さらに、前記推定されたウォッベ指数及び前記ガスの測定された前記密度から、前記高位熱容量を推定するように構成されている。
詳細には、高位発熱量は、測定された密度の二乗を推定されたウォッベ指数に乗ずることによって得られる。
参考として、スイスの供給業者ＴＲＡＦＡＧからの参照番号８７７４を有する密度センサを使用することが可能であることに留意するべきである。
特定の実施形態では、前記燃料ガスの前記少なくとも１つの燃焼特性は、前記高位発熱量ＰＣＳGN/H2を含み、前記装置は、
・前記燃料ガスの流れを受け入れるための入口、
・圧力降下を生じさせる装置を通る層流の状態にある前記燃料ガスの前記質量流量を測定するためのモジュールであって、前記測定は、前記燃料ガスの粘度及び基準ガスの粘度に依存するモジュール、及び、
・圧力降下を生じさせる前記装置の下流で、熱式質量流量計を用いて前記燃料ガスの前記質量流量を測定するためのモジュールであって、前記測定は、前記燃料ガスの比熱容量及び基準ガスの熱容量に依存するモジュール、
を備え、
前記高位発熱量ＰＣＳGN/H2を推定するために使用される前記経験的なアフィン関係は、次の式で記述される。
Figure 0007258757000010
Figure 0007258757000011
ここで、
Ｚは、変数Ｙに対応する変数であり、
Ｑmes,1は、測定された圧力降下を生じさせる装置を通る層流の状態にある前記燃料ガスの前記質量流量であり、
Ｑmes,2は、圧力降下を生じさせる前記装置の下流で測定された前記燃料ガスの前記質量流量であり、
Ａ，Ｂ，及びＣは、前記燃料ガスのグループに関する所定の係数であり、それぞれ、前記係数α，β，及びγに対応する。
特定の実施形態では、前記経験的なアフィン関係を用いて、前記少なくとも１つの特性ΞGN/H2を推定するために構成されている前記モジュールは、さらに、空気の化学量論的な体積又は燃焼性指数を推定するように構成されている。
特定の実施形態では、前記装置は、前記燃料ガスの前記燃焼特性を調節するためのモジュール、又は、燃料ガスの前記燃焼特性と、推定された空気の化学量論的な体積もしくは前記推定された特性に対する推定された燃焼性指数とを調節するためのモジュールをさらに備える。
この調節モジュールは、特に、例えば、空気のような追加のガスを導入するためのアクチュエータを備えていてもよい。調節は、閉ループ制御で行われてもよい。
本発明の他の特徴及び利点は、添付の図面を参照して行われる以下の説明から明らかであり、添付の図面は、非限定的な特徴を有する例を示す。
図１は、燃焼特性を推定する方法の一例の工程を示す図である。 図２は、燃焼特性を推定するための装置の一例を示す図である。 図３は、ウォッベ指数を推定する方法の工程の図である。 図４は、ウォッベ指数を推定するための装置の一例の図である。 図５は、高位発熱量を推定する方法の工程の図である。 図６は、高位発熱量を推定するための装置の一例を示す図である。 図７は、ウォッベ指数及び高位発熱量を推定するための装置の一例を示す図である。 図８は、ウォッベ指数を調節するための装置の一例の図である。
以下、燃料ガスのグループに属する燃料ガスの燃焼特性を推定するための方法及び装置について説明する。特性は、ウォッベ指数又は高位発熱量であってもよい。
本発明は、単一の特性を推定するいずれかの方法に限定されず、ウォッベ指数及び高位発熱量を同時に推定することを含んでいてもよい。
図１は、燃料ガスの燃焼特性を推定する方法の工程を示す図である。
この方法は、特に、ゼロでない量の水素を含みうるガス、特に、０％から２０％の範囲で水素を含む燃料ガスに適している。
第１の工程Ｅ０１において、燃料ガスの少なくとも２つの流動特性が測定される。測定対象の燃料ガスは、次に流動部材を流れ、温度、圧力、又は流量などの流動特性が、例えば、センサを用いて測定される。
第２の工程Ｅ０２において、流動ガスにおける水素の含有量が、測定され、ＸH2で記述される。本発明の実施及び実施形態の全てにおいて、水素の含有量は、モル分率又は体積分率であってもよいことに留意するべきである。圧力が十分に低ければ、理想気体の法則が適用されると仮定してもよく、モル分率及び体積分率が同じ値を有すると仮定してもよい。本発明の適用において、水素の含有量は、モル分率又は体積分率と同じ値を有する。この工程は、工程Ｅ０１と同時に実施されてもよく、又は工程Ｅ０１の前に実施されてもよく、又は工程Ｅ０１の後に実施されてもよい。
第３の工程Ｅ０３において、ΞGN/H2で記述される特性は、次の経験的なアフィン関係を用いて推定される。
Figure 0007258757000012
α，β，及びγは、燃料ガスのグループに関する所定の係数である。Ｙは、燃料ガスの少なくとも２つの流動特性の測定値から生成される燃料ガスの物理特性を表す変数である。
燃料ガスの少なくとも２つの流動特性の測定値は、特に、変数Ｙが、粘度、比熱容量（単位重量当たり）、又は密度などの燃料ガスの物理特性を表すことができるように、また、燃料ガスの流動特性の測定値に応じて表現されうるように、選択される。
したがって、流動特性Ｕ1の測定に関してＵmes,1の記号を使用し、流動特性Ｕ2の測定に関してＵmes,2の記号を使用し、流動特性Ｕ3の測定に関してＵmes,3の記号を使用することによって、以下が与えられる。
Figure 0007258757000013
図２は、燃料ガスのグループに属する燃料ガスの少なくとも１つの燃焼特性を推定するための装置の実施形態を示す。
本装置は、特に、図１を参照して説明した種類の方法を実行できる。
この例では、装置１は、例えば、ウォッベ指数又は高位発熱量を推定することが要求される燃料ガスの流れを受け入れるように構成されている流動部材２を備える。流動部材は、燃料ガスの流れを受け入れるための入口２ａと、出口２ｂとを有する。
装置１は、また、燃料ガスの流動特性を測定するための２つのモジュール３ａ及び３ｂを有する。
一例として、モジュール３ａ及び３ｂは、両方とも、又はそれぞれ、流量、圧力、又は温度を測定してもよい。
具体的には、モジュール３ａ及び３ｂは、流量、圧力、又は温度を測定するための従来型のセンサである。
装置１は、水素の含有量ＸH2を測定するためのモジュール４、例えば、水素のモル百分率を送信するセンサを含む。
モジュール３ａ，３ｂ，及び４は、それらのそれぞれの測定の結果を推定モジュール５に伝えるように、推定モジュール５と通信する。推定モジュール５は、プロセッサ及びメモリ（図示しない）を有するコンピュータであってもよい。
モジュール５は、次の経験的なアフィン関係を用いて、少なくとも１つの特性ΞGN/H2を推定するように構成されている。
Figure 0007258757000014
α，β，及びγは、燃料ガスのグループに関する所定の係数である。Ｙは、モジュール３ａ及び３ｂによって測定された燃料ガスの少なくとも２つの流動特性の値から生成される燃料ガスの物理特性を表す変数である。
この目的のために、モジュール５は、メモリに格納された、所定の係数α，β，及びγについて取りうる値を有していてもよい。モジュール５は、さらに、推定を実行することを可能にするために、メモリに格納されたコンピュータプログラムの命令を有していてもよい。
一例として、コンピュータプログラムは、モジュール３ａ及び３ｂによって測定された値から変数Ｙを計算するための命令と、上記の定義された関数を用いて特性ΞGN/H2を計算するための命令と、を含んでいてもよい。
任意には、装置１は、少なくとも１つの特性を調節するために、推定モジュール５によって制御されるアクチュエータ６を含んでいてもよい。一例として、アクチュエータ６は圧縮空気のインジェクタであってもよい。
図３及び４を参照して、燃料ガスのウォッベ指数ＩＷGN/H2を推定する方法及び実施形態が説明される。
図３は、燃料ガスのウォッベ指数ＩＷGN/H2を推定する方法の工程を示す。
この方法は、流体狭窄部（例えば、オリフィス又はマイクロノズル）を通る音速流の状態にある基準ガス（例えば、メタン）の質量流量の測定が、測定された基準ガスの絶対圧力及び測定された基準ガスの絶対温度で行われる校正の第１工程Ｅ１１を含む。
第２の工程Ｅ１２では、流体狭窄部を通って音速流の状態にある燃料ガスの質量流量の測定が行われる。この測定は、狭窄部の上流で測定された絶対圧力及び狭窄部の上流で測定された絶対温度で行われる。
工程Ｅ１３では、水素の含有量ＸH2が測定される。
工程Ｅ１１からＥ１３は、任意の可能な順番で実行されてもよい。特に、Ｅ１２及びＥ１３は、同時に実行されてもよい。
工程Ｅ１４では、ウォッベ指数ＩＷGN/H2は、以下に記述された経験的なアフィン関係によって推定される。
Figure 0007258757000015
Figure 0007258757000016
Ｑmes,2は、測定された燃料ガスの質量流量である。Ｐmesは、測定された燃料ガスの絶対圧力である。Ｔmesは、測定された燃料ガスの絶対温度である。Ｑrefは、測定された基準ガスの質量流量である。Ｐrefは、測定された基準ガスの絶対圧力である。Ｔrefは、測定された基準ガスの絶対温度である。Ｄ，Ｅ，及びＦは、燃料ガスのグループに関する所定の係数であり、それぞれ、係数α，β，及びγに対応する。
工程Ｅ１１及びＥ１２を実行するために、形状が知られているオリフィス又はマイクロノズルなどの流体狭窄部が使用される。
流体狭窄部のいずれかの側で、次の測定が行われる。流体狭窄部の上流で絶対温度及び絶対圧力が測定される。流体狭窄部の下流で質量流量が測定される（例えば、熱式質量流量計を用いて）。
具体的に、それは、装置を通るＱvnで記述される標準体積流量を２つの異なる方法で測定することに対応する。
特に、流れが音速であるため、以下が適用される。
Figure 0007258757000017
ｋは、オリフィスの形状の特徴によって定まる定数である。ｐ及びＴは、オリフィスの上流におけるガスの絶対圧力及び絶対温度である。ｄは、ガスの密度である。
さらに、熱式質量流量計を使用することによって、以下が得られる。
Figure 0007258757000018
Ｑvn,mes,2は、工程Ｅ１２で、熱式質量流量計によって測定された測定値である。Ｃは、空気の特性と比較されたガスの物理特性（例えば、その比熱容量、その粘度、その熱伝導率）の差を考慮した補正係数である。
既知の組成のガスに対して、以下が適用される。
Figure 0007258757000019
Ｘiは、成分ｉの体積分率である。Ｃiは、成分ｉに関する補正係数である。この係数は、熱式質量流量計に関連付けられたテーブルから読み取られうる。
校正工程Ｅ１１は、正確に知られている組成の基準ガス（好ましくは、純粋なメタン）を装置に供給する工程である。
工程Ｅ１１の校正手順で得られた値を使用する２つの流量の関係の間の方程式は、したがって、次のように記述されうる。
Figure 0007258757000020
そして、燃料ガスの測定が行われる工程Ｅ１２の測定手順において、以下が適用される。
Figure 0007258757000021
Ｅ１４での測定を用いて、水素の含有量ＸH2（例えば、体積値で）も判明する。
上記の方程式において、唯一判明していない２つのパラメータは、Ｃ及び密度ｄであることに留意するべきである。したがって、Ｃ及び密度ｄに関連付けられた変数Ｙを表すために、パラメータｋを除くことが可能である。
具体的に、以下が適用される。
Figure 0007258757000022
Ｙは、燃料ガスの物理特性を表す変数である。測定値からＹを書き換えることが可能である。これらの測定値（Ｑmes,2，ｐ，及びＴ）は、測定された流動特性Ｕ1に関してＵmes,1と記述され、測定された流動特性Ｕ2に関してＵmes,2と記述され、及び、測定された流動特性Ｕ3に関してＵmes,3と記述されうる。したがって、以下が与えられる。
Figure 0007258757000023
Ｕmes,1＝Ｐmes、
Ｕmes,2＝Ｑmes,2、及び、
Ｕmes,3＝Ｔmesである。
最終的に、以下が適用される。
Figure 0007258757000024
変数Ｙのこの定義は、その変数Ｙを有する経験的なアフィン関係において使用されうる。そして、それは、以下の形式で記述される。
Figure 0007258757000025
Ｄ，Ｅ，及びＦは、燃料ガスのグループに関する所定の係数であり、それぞれ、図１及び２を参照して説明した係数α，β，及びγに対応する。
係数Ｄ，Ｅ，及びＦは、既知のガスの組成、例えば、国又は地域のネットワークにおける天然ガスの組成から得られうる。一例として、ヨーロッパで流通し、当業者によく知られている高い発熱量を有するガスに関する既知の組成を使用することが可能である。
これらの既知の組成から、様々な化合物のモル分率に関する限定値を規定することが可能である。例えば、ガス中の化学種Ｋのモル分率に関する表記ＸKを使用することによって、以下の種類の既知の組成を持つことが可能である。
Figure 0007258757000026
今日のガス組成には、水素がないことに留意するべきである。
これらの考えうる範囲は、ウォッベ指数（又は高位発熱量でさえ）を決定することが可能であるガスの組成を、ランダムに生成させることを可能にする。したがって、１０，０００種類のガスが生成されうる。
ランダムに生成された各ガスに対して、０％から２０％のモル分率の範囲のランダムな量の水素を加えることが可能である。また、これらのランダムに生成されたガスに関して、ウォッベ指数を決定することが可能である。
ランダムに生成された組成から、変数Ｙに関連する値を推定することも可能である。最小二乗法を適用することによって、係数Ｄ，Ｅ，及びＦに関する値を得ることが可能である。
Figure 0007258757000027
これらの結果は、基準ガスとしてメタンを使用することによって得られた。
０％から２０％の範囲にある水素のモル分率に関して、ウォッベ指数の推定における誤差は、常に、１．４％未満であるとともに、誤差は、ランダムに得られた１０，０００種類のガスの９８％において、１％未満であることに留意するべきである。
あるいは、係数が定数ではなく、モル分率としての水素の含有量に依存する相関関係を使用することが可能である。
経験的なアフィン関係は次のように書き換えられてもよい。
Figure 0007258757000028
この例では、Ｄ及びＥは、水素の含有量の関数である（方程式は、ＸH2に依存する２つの係数のみを示すように書き換えられている。それにもかかわらず、３つの係数を有する方程式を記述することが可能である）。
同様に、ランダムなガス組成を生成することによって、本発明者らは、幅１％の水素の濃度範囲にわたる区分において一定である関数Ｄ及びＥを使用することが可能であることを見出した。
例えば、２％から３％の水素の範囲で、Ｄ及びＥは次の値をとる。
Figure 0007258757000029
ウォッベ指数の推定における誤差は、したがって、１０，０００種類のランダムガスの９５％について、常に１％未満であることに留意するべきである。
図４は、図３を参照して説明された方法を実施できる、ウォッベ指数を推定するための装置を示す。
この装置１０は、燃料ガスＧ１を受け入れるための入口１１、及び、基準ガスＧ２（典型的にはメタン）を受け入れるための入口１２を有する。
装置１０は、また、燃料ガスの流れ又は基準ガスの流れを管の入口１３に導くためのセレクタ及び案内モジュールを有する。具体的には、セレクタ及び案内モジュールは、弁３１及び弁３５を含む。
弁３１及び３５は、燃料ガス又は基準ガスのいずれかが管に流入するように、各端子５２及び５１を介してモジュール３６によって制御される。
この例では、管は、その入口１３から始まって上流から下流へ：
・端子４１を介してモジュール３６に接続された絶対圧力センサ６０、
・端子４２を介してモジュール３６に接続された絶対温度センサ５０、
・流体狭窄部３２（例えば、オリフィス又はマイクロノズル）、
・端子４３を介してモジュール３６に接続された熱式質量流量センサ３３、及び、
・端子４４を介してモジュール３６に接続された、水素のモル分率を計測するためのセンサ３４、
を備えている。
モジュール３６は、端子４１から４４で受信した信号からウォッベ指数を推定するために、上記したような経験的なアフィン関係を適用しうる。
示された例では、モジュール３６は、２つの出力端子５４及び５３を介して外部モジュール３７、例えば、調節を実施するためのアクチュエータ、又はディスプレイと通信する。図示しない変形例では、モジュール３７は、装置１０に組み込まれている。
図５及び６を参照して、燃料ガスの高位発熱量を推定する方法及び実施形態を説明する。
図５は、燃料ガスの高位発熱量ＰＣＳGN/H2を推定する方法の工程を示す。
工程Ｅ２２において、圧力降下を生じさせる装置を通る層流において、燃料ガスの質量流量が測定される。この測定は、燃料ガスの粘度及び基準ガスの粘度に依存する。
工程Ｅ２３において、圧力降下を生じさせる装置の下流で、燃料ガスの質量流量が熱式質量流量計を用いて測定される。測定は、燃料ガスの比熱容量及び基準ガスの熱容量に依存する。
工程Ｅ２２及びＥ２３を実行する前又は後に実行されうる工程Ｅ２１では、工程Ｅ２２及びＥ２３を実行することに対応する校正工程が実行される。ただし、校正工程は、基準ガス（例えば、メタン）を用いて実行される。
工程２４もまた、モル分率として水素の含有量ＸH2を測定するために実施される。
最後に、工程Ｅ２５において、高位発熱量が、次の式を有する経験的なアフィン関係を用いて推定される。
Figure 0007258757000030
Figure 0007258757000031
Ｚは、図１を参照して説明した変数Ｙに対応する変数である。Ｑmes,1は、測定された圧力降下を生じさせる装置を通る層流の状態にある燃料ガスの質量流量である。Ｑmes,2は、圧力降下を生じさせる装置の下流で測定された燃料ガスの質量流量である。Ａ，Ｂ，及びＣは、燃料ガスのグループに関する所定の係数であり、それぞれ、係数α，β，及びγに対応する。
高位発熱量（ＰＣＳ）に関する次の関係は、先行文献DE 4 118 781から知られている。
Figure 0007258757000032
α及びβは所定の係数である。ρは、燃料ガスの密度である。μは、燃料ガスの粘度である。Ｃpは、燃料ガスの熱容量である。
次の無次元変数Ｚは、次のように定義されうる。
Figure 0007258757000033
ここで、下付きのＧＮは、燃料ガスを示し、ｒｅｆは基準ガスを示す。
高位発熱量に関する関係は、次のように書き換えることができる。
Figure 0007258757000034
Ｚを測定するために、２つの流量の測定が使用される。一例として、Ｑmes,1は、測定された圧力降下を生じさせる装置を通る層流の状態にある燃料ガスの質量流量である。ポアズイユの法則を使用することによって、この装置では、層流要素を通って、ガスの流れの圧力を降下させることが知られている。
測定は、粘度に依存するので、特定の測定装置に関し、それは上流と下流との圧力差ΔＰmesに依存する形式を持ち、Ｑmes,1を送信する流量計は、次式を用いて体積流量を決定しうる。
Figure 0007258757000035
ここで、Ｋ1は、幾何学的定数である。
この例では、圧力センサ及び温度プローブが流量計に組み込まれており、それらによれば、以下のように、質量流量Ｑmes,1（又は通常の温度及び圧力の条件に換算された体積流量）を直接得ることができる。
Figure 0007258757000036
修正された体積流量Ｑmes,1の測定は、燃料ガスの粘度に依存する。
流量Ｑmes,2は、熱式質量流量計を使用して測定される。質量流量Ｑmes,2の測定は、燃料ガスの比熱容量に依存する。
この方法では、質量流量の２つの推定が得られ、各推定は、実際の通常の体積流量（Ｑvn）に対するそれぞれの誤差を示す。以下が適用される。
Figure 0007258757000037
Ｑmes,1について、誤差は、混合ガスの実際の粘度と、基準ガス（この例ではメタン）の粘度との差に関連する。
Ｑmes,2について、誤差は、混合ガスの実際の比熱容量を密度に乗じた値と、基準ガス（この例ではメタン）の同量を密度に乗じた値との間の差に関連する。
２つの測定された流量の比率をとることによって、変数Ｚが決定され、文献DE 4 118 781の相関関係は次のようになる。
Figure 0007258757000038
Ｕmes,1＝Ｑmes,1、及び
Ｕmes,2＝Ｑmes,2で記述される。
以下が適用される。
Figure 0007258757000039
それにもかかわらず、本発明者らは、水素が燃料ガス中に存在すると、この関係は適用できないことを見出した。具体的には、測定誤差が非常に大きくなる。水素は、アルカンのようなガスの特性とは大きく異なる特性（粘度、比熱容量）を有する。
その欠点を軽減するために、測定された水素の含有量ＸH2が使用され、以下の方程式が提案される。
Figure 0007258757000040
Ａ，Ｂ，及びＣは、燃料ガスのグループに関する所定の係数であり、それぞれ、図１及び２を参照して説明した係数α，β，及びγに対応する。
係数Ａ，Ｂ，及びＣは、既知のガスの組成、例えば、国又は地域のネットワークにおける天然ガスの組成から得られうる。一例として、ヨーロッパで流通し、当業者によく知られている高い発熱量を有するガスに関する既知の組成を使用することが可能である。
これらの既知の組成から、様々な化合物のモル分率に関する限定値を規定することが可能である。例えば、ガス中の化学種Ｋのモル分率に関する表記ＸKを使用することによって、以下の種類の既知の組成を持つことが可能である。
Figure 0007258757000041
今日のガス組成には、水素がないことに留意するべきである。
これらの考えうる範囲は、高位発熱量を決定することが可能であるガスの組成を、ランダムに生成させることを可能にする。したがって、１０，０００種類のガスが生成されうる。
ランダムに生成された各ガスに対して、０％から２０％のモル分率の範囲のランダムな量の水素を加えることが可能である。ここでも同様に、これらのランダムに生成されたガスの高位発熱量を決定することが可能である。
ランダムに生成された組成から、変数Ｚに関連する値を推定することも可能である。最小二乗法を適用することによって、係数Ａ，Ｂ，及びＣに関する値を得ることが可能である。
Figure 0007258757000042
あるいは、係数が定数ではなく、モル分率としての水素の含有量に依存する相関関係を使用することが可能である。
経験的なアフィン関係は次のように書き換えられてもよい。
Figure 0007258757000043
この例では、Ａ及びＢは、水素の含有量の関数である（方程式は、ＸH2に依存する２つの係数のみを示すように書き換えられている。それにもかかわらず、３つの係数を有する等式を記述することが可能である）。
同様に、ランダムなガス組成を生成することによって、本発明者らは、幅１％の水素の濃度範囲にわたる区分において一定である関数Ａ及びＢを使用することが可能であることを見出した。
例えば、１９％から２０％の水素の範囲で、Ａ及びＢは次の値をとる。
Figure 0007258757000044
図６は、図５を参照して説明された方法を実施するのに適している、高位発熱量を推定するための装置を示す。
この装置１００は、燃料ガスＧ１を受け入れるための入口１１１、及び、基準ガスＧ２（典型的にはメタン）を受け入れるための入口１１２を有する。
装置１００は、また、燃料ガスの流れ又は基準ガスの流れを管の入口１１３に導くためのセレクタ及び案内モジュールを有する。具体的には、セレクタ及び案内モジュールは、弁１３１及び弁１３５を含む。
弁１３１及び１３５は、燃料ガス又は基準ガスのいずれかが管に流入するように、各端子１５２及び１５１を介してモジュール１３６によって制御される。
この例では、管は、その入口１１３から始まって上流から下流へ：
・圧力降下を生じさせる装置を通る層流において質量流量を計測するためのセンサ１３２であって、測定は、燃料ガスの粘度及び基準ガスの粘度に依存し、端子１４２を介してモジュール１３６に接続されているセンサ１３２、
・熱式質量流量を計測するためのセンサ１３３であって、測定は、測定ガスの比熱容量及び基準ガスの比熱容量に依存し、端子１４３を介してモジュール１３６に接続されているセンサ１３３、及び、
・水素のモル分率を計測し、端子１４４を介してモジュール１３６に接続されたセンサ１３４、
を備えている。
モジュール１３６は、端子１４２から１４４で受信した信号から高位発熱量を推定するために、上記したような経験的なアフィン関係を適用しうる。
示された例では、モジュール１３６は、２つの出力端子１５４及び１５３を介して外部モジュール１３７、例えば、調節を実施するためのアクチュエータ、又はディスプレイと通信する。図示しない変形例では、モジュール１３７は、装置１００に組み込まれている。
図７は、２つの経験的なアフィン関係を用いて高位発熱量及びウォッベ指数を同時に決定できる装置の例を示す。
この装置は、圧縮ガスのための２つの入口を有する。第１の入口１０１１は、ガスのグループ（例えば、フランスの規格ＮＦ ＥＮ ４３７に準拠した第２のグループにおける燃料ガス）に属しており、かつ水素を含んでいる燃料ガスＧ１を受け入れるための入口である。装置１０００は、校正段階で使用される基準ガスＧ０のための別の入口１０１２を有する。
装置１０００は、また、燃料ガスＧ１の流れ又は基準ガスＧ０の流れを管の入口１０１３に導くためのセレクタ及び案内モジュールを有する。セレクタ及び案内モジュールは、弁１０３１及び弁１０３５を含む。
弁１０３１及び１０３５は、燃料ガス又は基準ガスのいずれかが管に流入するように、各端子１０５２及び１０５１を介してモジュール１０３６によって制御される。
管の入口１０１３から始まって、上流から下流へ、装置は以下を備える。
・端子１０４１を介してモジュール１０３６に接続された絶対圧力センサ１００６。このセンサは、燃料ガスについて測定された値Ｐmes及び基準ガスについて測定された値Ｐrefを送信する。
・端子１０４２を介してモジュール１０３６に接続された絶対温度センサ１００５。このセンサは、燃料ガスについて測定された値Ｔmes及び基準ガスについて測定された値Ｔrefを送信する。
・流体狭窄部１０３２（例えば、オリフィス又はマイクロノズル）。
・端子１０４３を介してモジュール１０３６に接続され、圧力降下を生じさせる装置を通る層流における質量流量を計測するためのセンサ。このセンサは、値Ｑmes,1を送信する。
・端子１０４６を介してモジュール１０３６に接続された熱式質量流量センサ１０３３。このセンサは、燃料ガスについて測定された値Ｑmes,2及び基準ガスについて測定された値Ｑrefを送信する。
・水素のモル分率を計測し、端子１０４７を介してモジュール３６に接続されたセンサ１０３４。
・排出口１０３９。
この例では、モジュール１０３６は、次の２つの変数を計算しうる。
Figure 0007258757000045
モジュール１０３６は、したがって、ウォッベ指数ＩＷGN/H2、高位発熱量ＰＣＳGN/H2、空気の化学量論的な体積ＶａGN/H2、燃焼性指数ＢGN/H2、及び混合された燃料ガスの密度ｄGN/H2を決定するために次の式を使用しうる。
Figure 0007258757000046
具体的に、ウォッベ指数ＩＷ及び高位発熱量ＰＣＳに関する値から、空気の化学量論的な体積Ｖａ、燃焼性指数Ｂ、及び燃料ガス混合物の密度を決定することが可能である。
より正確には、ヨーロッパで通常流通している天然ガス（現在、水素が含まれていない）に対して、以下が適用される。
Figure 0007258757000047
本発明者らは、水素の存在下では、値が１．１６２から逸脱するであろうことを見出した。
調査対象の全ての天然ガスに対して０．３％未満の誤差で、空気の化学量論的な体積に対する高位発熱量の比率と、（既知の）水素の含有量とを関連付けることが可能である。
言い換えれば、１％未満の誤差で水素の含有量及び推定された高位発熱量が判明することによって、空気の化学量論的な体積が判明することが可能となり、したがって、開ループ制御にて燃焼を調節することが可能となる。
同様に、燃焼性指数に対するウォッベ指数の比率は、水素の含有量に比例して変化する。このことは、また、開ループ制御を実行することを可能にする。
１０，０００種類のガスが使用される場合、燃焼性指数の推定における誤差は、ランダムな組成を含む１０，０００種類のガスの９９．５％において１％未満であることに留意するべきである。
ウォッベ指数の定義によれば、ウォッベ指数及び高位発熱量の両方が判明すれば、ガスの密度を決定することが可能である。
装置１０００は、端子１０５４でウォッベ指数の信号を、端子１０５５で高位発熱量の信号を、端子１０５６で密度の信号を、端子１０５７で空気の化学量論的な体積の信号を、及び端子１０５８で燃焼性指数の信号を送信することに留意するべきである。これらの信号は、調査対象の燃料ガス、すなわち、ガスＧ１に関連する。
図８は、図４を参照して説明した装置１０を利用する調節の適用例を示す。この装置は、管の２つの点、入口点２０１と出口点２０２との間に接続される。装置１０は、センサ群２０５によって行われた測定に応じて、圧縮空気の供給モジュール２０３を制御する。
この装置は、特に、ヨーロッパで流通しているガスに追加量の水素を添加した類似したガスについて、ウォッベ指数を継続的に調節できる。特に、測定が継続的又はリアルタイムで行われるため、調節を実行することが可能になる。これに対し、クロマトグラフの種類の装置でそのようなことは可能ではない。
ウォッベ指数の測定の正確性は、１％のオーダーである。
所望の調節を得るために、装置１０のウォッベ指数に関する設定値を修正することが可能であることに留意するべきである。
さらに、校正工程は、メタンを使用している間、自動的に実行されうる。開始段階もまた、自動的に、すなわち、操作者が介入することなく実行されうる。
さらに、特性の１つ（例えば、ウォッベ指数）について閉ループ制御が実行され、導入される圧縮空気の流量について開ループ制御が実行されるといった、組み合わせの制御も可能である。これにより、特性に関する設定値がより厳密に満たされるように、特性の変動を特に正確な方法で考慮に入れることが可能となる。
そのような装置は、毎時１５０リットル（Ｌ／ｈ）未満の速度で天然ガスを消費してもよいことに留意するべきである。
上記の方法及び実施形態は、熱力学的データ（Ｃｐ、粘度）を、体積分率又はモル分率での水素の含有量の測定と組み合わせることを可能にする。適切な相関関係を使用することによって、ウォッベ指数及び／又は高位発熱量を計算することが可能となる。
気相クロマトグラフを使用することと比較して、以下が適用されることに留意するべきである。
・推定は、ほとんど瞬時に行われる（応答時間は、５秒未満でありうる）。
・推定は、安価である。
既存の相関方法と比較して、水素が存在するにもかかわらず、良好な正確性が得られるとともに、良好な測定の堅牢性も得られる。
最後に、燃焼を利用する装置（燃焼計又は熱量計）と比較して、本発明は、良好な堅牢性、実施の容易さ、並びに、メンテナンスが削減されること、及びメンテナンスが容易に実施できることを提供する。

Claims (13)

  1. 燃料ガスのグループに属する燃料ガスの少なくとも１つの燃焼特性を推定する方法であって、前記少なくとも１つの燃焼特性は、ウォッベ指数又は高位発熱量であり、前記方法は、
    ・前記燃料ガスの少なくとも２つの流動特性を測定すること（Ｅ０１）と、
    ・前記燃料ガスに含まれている水素の含有量ＸH2を測定すること（Ｅ０２）と、
    ・前記少なくとも１つの燃焼特性ΞGN/H2は、次の経験的なアフィン関係を用いて推定されること（Ｅ０３）と、
    を含み、
    Figure 0007258757000048
     ここで、
    α，β，及びγは、前記燃料ガスのグループに関する所定の係数であり、
    Ｙは、前記燃料ガスの前記少なくとも２つの流動特性の測定値から生成される前記燃料ガスの物理特性を表す変数であり、
    前記燃料ガスの前記少なくとも１つの燃焼特性が、前記ウォッベ指数ＩＷ GN/H2 を含む場合、前記燃料ガスの少なくとも２つの流動特性の前記測定は、流体狭窄部を通る音速流の状態にある前記燃料ガスの質量流量を測定すること（Ｅ１２）を含み、前記測定は、前記流体狭窄部の上流で測定された絶対圧力及び前記流体狭窄部の上流で測定された絶対温度で行われ、
    前記方法は、測定された基準ガスの絶対圧力及び測定された基準ガスの絶対温度で、前記流体狭窄部を通る音速流の状態にある基準ガスの質量流量の測定が行われる校正手順（Ｅ１１）をさらに含み、
    前記ウォッベ指数ＩＷ GN/H2 を推定する（Ｅ１４）ために使用される前記経験的なアフィン関係は、次の式で記述され、
    Figure 0007258757000049
    Figure 0007258757000050
     ここで、
    Ｑ mes,2 は、測定された前記燃料ガスの前記質量流量であり、
    Ｐ mes は、測定された前記燃料ガスの前記絶対圧力であり、
    Ｔ mes は、測定された前記燃料ガスの前記絶対温度であり、
    Ｑ ref は、測定された前記基準ガスの前記質量流量であり、
    Ｐ ref は、測定された前記基準ガスの前記絶対圧力であり、
    Ｔ ref は、測定された前記基準ガスの前記絶対温度であり、
    Ｄ，Ｅ，及びＦは、前記燃料ガスのグループに関する所定の係数であり、それぞれ、前記係数α，β，及びγに対応し、
    前記燃料ガスの前記少なくとも１つの燃焼特性が、前記高位発熱量ＰＣＳ GN/H2 を含む場合、
    前記燃料ガスの少なくとも２つの流動特性の前記測定は、
    ・圧力降下を生じさせる装置を通る層流の状態にある前記燃料ガスの質量流量を測定（Ｅ２２）し、前記測定は、前記燃料ガスの粘度及び基準ガスの粘度に依存すること、
    ・圧力降下を生じさせる前記装置の下流で、熱式質量流量計を用いて前記燃料ガスの前記質量流量を測定（Ｅ２３）し、前記測定は、前記燃料ガスの比熱容量及び前記基準ガスの前記粘度に依存すること、
    を含み、
    前記高位発熱量ＰＣＳ GN/H2 を推定する（Ｅ２５）ために使用される前記経験的なアフィン関係は、次の式で記述される、方法。
    Figure 0007258757000051
    Figure 0007258757000052
     ここで、
    Ｚは、変数Ｙに対応する変数であり、
    Ｑ mes,1 は、測定された圧力降下を生じさせる装置を通る層流の状態にある前記燃料ガスの前記質量流量であり、
    Ｑ mes,2 は、圧力降下を生じさせる前記装置の下流で測定された前記燃料ガスの前記質量流量であり、
    Ａ，Ｂ，及びＣは、前記燃料ガスのグループに関する所定の係数であり、それぞれ、前記係数α，β，及びγに対応する。
  2. 前記係数α，β，及びγは、入力として前記水素の含有量ＸH2に関する前記測定値を有し、出力として前記係数α，β，及びγを提供する表から読み取られた係数である、請求項１に記載の方法。
  3. 前記表は、前記係数α，β，及びγを、前記水素の含有量ＸH2の数値範囲と関連付けており、前記数値範囲は、１％の幅を有する、請求項２に記載の方法。
  4. 前記係数α，β，及びγの値は、前記Ｙの値及び前記燃焼特性が知られている、前記燃料ガスのグループの既知のガスに関するデータセットから得られる、請求項１から３のいずれか１項に記載の方法。
  5. 前記燃料ガスのグループの既知のガスに関する前記データセットから、ランダムに生成されたガスの、燃焼特性と前記Ｙに関する値とを関連付ける関数における、係数を得ることを含む、請求項４に記載の方法。
  6. 前記ウォッベ指数及び前記高位発熱量は、２つの前記経験的なアフィン関係を用いて推定される、請求項１から５のいずれか１項に記載の方法。
  7. 前記推定されたウォッベ指数及び前記推定された高位発熱量から前記燃料ガスの密度を推定することをさらに含む、請求項６に記載の方法。
  8. 前記燃料ガスの前記燃焼特性を調節すること、又は、燃料ガスの前記燃焼特性と、推定された空気の化学量論的な体積もしくは前記推定された燃焼特性に対応する推定された燃焼性指数とを調節することを含む、請求項１から７のいずれか１項に記載の方法。
  9. 前記燃料ガスの前記ウォッベ指数が推定され、前記方法は、前記燃料ガスの密度を測定すること、及び前記推定されたウォッベ指数と測定された前記燃料ガスの密度とから前記高位発熱量を推定すること、をさらに含む、請求項１から８のいずれか１項に記載の方法。
  10. 燃料ガスのグループに属する燃料ガスの少なくとも１つの燃焼特性を推定するための装置であって、前記少なくとも１つの燃焼特性は、ウォッベ指数又は高位発熱量であり、前記装置は、
    ・前記燃料ガスの少なくとも２つの流動特性を測定するための少なくとも２つのモジュール（３ａ，３ｂ）と、
    ・前記燃料ガスに含まれている水素の含有量ＸH2を測定するためのモジュール（４）と、
    ・次の経験的なアフィン関係を用いて前記少なくとも１つの燃焼特性ΞGN/H2を推定するように構成されたモジュール（５）と、
    を含み、
    Figure 0007258757000053
     ここで、
    α，β，及びγは、前記燃料ガスのグループに関する所定の係数であり、
    Ｙは、前記燃料ガスの前記少なくとも２つの流動特性の測定値から生成される前記燃料ガスの物理特性を表す変数であり、
    前記燃料ガスの前記少なくとも１つの燃焼特性が、前記ウォッベ指数ＩＷ GN/H2 を含む場合、前記装置は、
    ・前記燃料ガスの流れを受け入れるための入口（１１）、
    ・基準ガスの流れを受け入れるための入口（１２）、
    ・前記燃料ガスの前記流れ又は前記基準ガスの前記流れを管（１３）に導くためのセレクタ及び案内モジュール（３１，３５）、
    ・流体狭窄部（３２）、及び、
    ・前記流体狭窄部を通る音速流の状態にある前記燃料ガスの質量流量を測定するためのモジュール（３３）であって、前記流体狭窄部の上流で絶対圧力を測定するためのサブモジュール、及び前記流体狭窄部の上流で絶対温度を測定するためのサブモジュール、を含むモジュール（３３）、
    を備え、
    前記ウォッベ指数ＩＷ GN/H2 を推定するために使用される前記経験的なアフィン関係は、次の式で記述され、
    Figure 0007258757000054
    Figure 0007258757000055
     ここで、
    Ｑ mes,2 は、測定された前記燃料ガスの前記質量流量であり、
    Ｐ mes は、測定された前記燃料ガスの前記絶対圧力であり、
    Ｔ mes は、測定された前記燃料ガスの前記絶対温度であり、
    Ｑ ref は、測定された前記基準ガスの前記質量流量であり、
    Ｐ ref は、測定された前記基準ガスの前記絶対圧力であり、
    Ｔ ref は、測定された前記基準ガスの前記絶対温度であり、
    Ｄ，Ｅ，及びＦは、前記燃料ガスのグループに関する所定の係数であり、それぞれ、前記係数α，β，及びγに対応し、
    前記燃料ガスの前記少なくとも１つの燃焼特性が、前記高位発熱量ＰＣＳ GN/H2 を含む場合、前記装置は、
    ・前記燃料ガスの流れを受け入れるための入口（１１１）、
    ・圧力降下を生じさせる装置を通る層流の状態にある前記燃料ガスの質量流量を測定するためのモジュール（１３２）であって、前記測定は、前記燃料ガスの粘度及び基準ガスの粘度に依存するモジュール（１３２）、及び、
    ・圧力降下を生じさせる前記装置の下流で、熱式質量流量計を用いて前記燃料ガスの前記質量流量を測定するためのモジュール（１３３）であって、前記測定は、前記燃料ガスの比熱容量及び基準ガスの熱容量に依存するモジュール（１３３）、
    を備え、
    前記高位発熱量ＰＣＳ GN/H2 を推定するために使用される前記経験的なアフィン関係は、次の式で記述される、装置。
    Figure 0007258757000056
    Figure 0007258757000057
     ここで、
    Ｚは、変数Ｙに対応する変数であり、
    Ｑ mes,1 は、測定された圧力降下を生じさせる装置を通る層流の状態にある前記燃料ガスの前記質量流量であり、
    Ｑ mes,2 は、圧力降下を生じさせる前記装置の下流で測定された前記燃料ガスの前記質量流量であり、
    Ａ，Ｂ，及びＣは、前記燃料ガスのグループに関する所定の係数であり、それぞれ、前記係数α，β，及びγに対応する。
  11. 前記装置は、前記燃料ガスの前記ウォッベ指数を推定するのに適しており、前記装置は、前記燃料ガスの密度を測定するためのモジュールをさらに備え、前記少なくとも１つの燃焼特性を推定するために構成されている前記モジュール（５）は、さらに、前記推定されたウォッベ指数及び前記燃料ガスの測定された前記密度から、前記高位発熱量を推定するように構成されている、請求項１０に記載の装置。
  12. 前記経験的なアフィン関係を用いて、前記少なくとも１つの燃焼特性ΞGN/H2を推定するために構成されている前記モジュールは、さらに、空気の化学量論的な体積又は燃焼性指数を推定するように構成されている、請求項１０又は１１に記載の装置。
  13. 前記燃料ガスの前記燃焼特性を調節するためのモジュール、又は、燃料ガスの前記燃焼特性と、推定された空気の化学量論的な体積もしくは前記推定された燃焼特性に対応する推定された燃焼性指数とを調節するためのモジュールをさらに備える、請求項１０から１２のいずれか１項に記載の装置。
JP2019535340A 2016-12-28 2017-12-15 水素を含みうるガスの燃焼特性を推定するための方法 Active JP7258757B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
FR1663468A FR3061302B1 (fr) 2016-12-28 2016-12-28 Procede d'estimation d'une caracteristique de combustion d'un gaz pouvant contenir du dihydrogene
FR1663468 2016-12-28
PCT/FR2017/053613 WO2018122490A1 (fr) 2016-12-28 2017-12-15 Procede d'estimation d'une caracteristique de combustion d'un gaz pouvant contenir du dihydrogene

Publications (2)

Publication Number Publication Date
JP2020514714A JP2020514714A (ja) 2020-05-21
JP7258757B2 true JP7258757B2 (ja) 2023-04-17

Family

ID=58228289

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019535340A Active JP7258757B2 (ja) 2016-12-28 2017-12-15 水素を含みうるガスの燃焼特性を推定するための方法

Country Status (13)

Country Link
US (1) US11592430B2 (ja)
EP (1) EP3563152B1 (ja)
JP (1) JP7258757B2 (ja)
KR (1) KR102478378B1 (ja)
CN (1) CN110312933B (ja)
BR (1) BR112019013088B1 (ja)
CA (1) CA3048290A1 (ja)
DK (1) DK3563152T3 (ja)
ES (1) ES2912536T3 (ja)
FR (1) FR3061302B1 (ja)
PL (1) PL3563152T3 (ja)
PT (1) PT3563152T (ja)
WO (1) WO2018122490A1 (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3109998B1 (fr) * 2020-05-06 2022-05-13 Grtgaz Procédé de détermination d’un modèle de calcul de la concentration d’un composé injecté dans un mélange fluide
JP2022072451A (ja) * 2020-10-29 2022-05-17 三菱重工業株式会社 燃料ガスの組成分析装置及び組成分析方法、並びに、この組成分析装置を備える原動機制御装置及びこの組成分析方法を含む原動機制御方法
US11391221B2 (en) * 2020-12-16 2022-07-19 Hamilton Sundstrand Corporation Mass flow metering method and system
WO2022230120A1 (ja) * 2021-04-28 2022-11-03 株式会社日立製作所 水素供給システム
DE102022205583A1 (de) * 2022-06-01 2023-12-07 Robert Bosch Gesellschaft mit beschränkter Haftung Verfahren zur Messung eines Massendurchflusses eines gasförmigen Mediums

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040062290A1 (en) 2000-12-26 2004-04-01 Remy Cordier Method and device for the evaluating a fuel gas wobbe index
DE102008029553B3 (de) 2008-06-21 2009-11-26 Elster Gmbh Verfahren zur Bestimmung der Gasqualität eines zumindest teilweise mit Biogas oder aufbereitetem Biogas versetzten Probengasen
JP2016111004A (ja) 2014-11-28 2016-06-20 パナソニックＩｐマネジメント株式会社 燃料電池システムの運転方法および燃料電池システムにおける燃料組成推定方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS53129990U (ja) * 1977-03-24 1978-10-16
DE2928739C2 (de) 1979-07-17 1981-03-19 Ruhrgas Ag, 4300 Essen Verfahren und Vorrichtung zur verbrennungslosen Messung und/oder Regelung der Wärmemengenzufuhr zu Gasverbrauchseinrichtungen
US5122746A (en) * 1990-06-12 1992-06-16 Gas Research Institute Hydrocarbon gas measurements using nuclear magnetic resonance
DE4118781C2 (de) 1991-06-07 2002-10-24 Rmg Mestechnik Gmbh Wobbezahlmesser
DE19736528A1 (de) 1997-08-22 1999-02-25 Ruhrgas Ag Verfahren zur verbrennungslosen Messung des Brennwertes von Brenngas
FR2910967B1 (fr) 2006-12-28 2009-04-03 Commissariat Energie Atomique Dispositif et procede de mesure continue de la concentration en goudrons dans un flux gazeux
EP2679659B1 (en) 2012-06-29 2016-08-24 Global Gateways Lux HoldCo S.A. Method and plant for production of a fuel gas from waste
CN105807027B (zh) * 2016-03-16 2018-11-09 新奥科技发展有限公司 气体能量计量方法及装置
CN106096763B (zh) 2016-06-03 2019-08-09 河北省电力建设调整试验所 一种预测运行机组的燃煤发热值的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040062290A1 (en) 2000-12-26 2004-04-01 Remy Cordier Method and device for the evaluating a fuel gas wobbe index
DE102008029553B3 (de) 2008-06-21 2009-11-26 Elster Gmbh Verfahren zur Bestimmung der Gasqualität eines zumindest teilweise mit Biogas oder aufbereitetem Biogas versetzten Probengasen
JP2016111004A (ja) 2014-11-28 2016-06-20 パナソニックＩｐマネジメント株式会社 燃料電池システムの運転方法および燃料電池システムにおける燃料組成推定方法

Also Published As

Publication number Publication date
CN110312933A (zh) 2019-10-08
CN110312933B (zh) 2022-04-05
KR102478378B1 (ko) 2022-12-16
FR3061302B1 (fr) 2019-05-31
EP3563152B1 (fr) 2022-02-16
DK3563152T3 (da) 2022-05-23
WO2018122490A1 (fr) 2018-07-05
PT3563152T (pt) 2022-05-13
KR20190102008A (ko) 2019-09-02
ES2912536T3 (es) 2022-05-26
EP3563152A1 (fr) 2019-11-06
US20200232962A1 (en) 2020-07-23
PL3563152T3 (pl) 2022-06-20
US11592430B2 (en) 2023-02-28
JP2020514714A (ja) 2020-05-21
BR112019013088B1 (pt) 2022-10-18
FR3061302A1 (fr) 2018-06-29
BR112019013088A2 (pt) 2019-12-17
CA3048290A1 (en) 2018-07-05

Similar Documents

Publication Publication Date Title
JP7258757B2 (ja) 水素を含みうるガスの燃焼特性を推定するための方法
US5486107A (en) Determination of fuel characteristics
RU2690099C2 (ru) Способ и измерительное устройство для определения удельных параметров для свойства газа
JP5647083B2 (ja) 原料濃度検出機構を備えた原料気化供給装置
JP4798774B2 (ja) 混合ガス供給システム
US5401162A (en) Microbridge-based combustion control
EP2241810B1 (en) Flow rate control device
US20170261480A1 (en) Method and measuring apparatus for determining physical properties of gas
CN101182926B (zh) 一种在线炉燃料流量控制装置
EP2241811B1 (en) Fuel supply device
JP2001013097A (ja) 可燃性ガスの火力の評価並びに調節のための装置と方法
US6893152B2 (en) Method and device for evaluating a fuel gas wobbe index
JP2012140896A (ja) ガスタービン燃料制御方法およびガスタービン燃料制御装置
US3783684A (en) Fuel gas flow-meter corrector equipment for gases having variable characteristics
JPS61159143A (ja) 天然ガスの発熱量調整装置
JP2534418B2 (ja) 熱量計
Meyer et al. Performance and validation tests on the NIST hybrid humidity generator
JP2019500586A (ja) 圧力値が可変な液体容器を有する冷却装置の試験バイパス
Meyer et al. Design and performance of the new NIST hybrid humidity generator
KR101172538B1 (ko) 연료 조성 변동시의 연료유량 정밀 측정 장치 및 그 방법
JP2023130036A (ja) 推定装置
JP2012041851A (ja) 熱量調整装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20201106

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210818

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210907

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20211206

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220307

A02 Decision of refusal

Free format text: JAPANESE INTERMEDIATE CODE: A02

Effective date: 20220802

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20221202

C60 Trial request (containing other claim documents, opposition documents)

Free format text: JAPANESE INTERMEDIATE CODE: C60

Effective date: 20221202

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20221205

A911 Transfer to examiner for re-examination before appeal (zenchi)

Free format text: JAPANESE INTERMEDIATE CODE: A911

Effective date: 20221228

C21 Notice of transfer of a case for reconsideration by examiners before appeal proceedings

Free format text: JAPANESE INTERMEDIATE CODE: C21

Effective date: 20230110

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: 20230307

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20230405

R150 Certificate of patent or registration of utility model

Ref document number: 7258757

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150