Method for estimating a combustion characteristic of a gas that may contain dihydrogen
11592430 · 2023-02-28
Assignee
Inventors
- Mathieu Ourliac (Saint Soupplets, FR)
- Sandra Capela (Pantin, FR)
- Laurent Lantoine (Conflans Sainte Honorine, FR)
- Naushad Manjoo (Eaubonne, FR)
Cpc classification
G01N25/22
PHYSICS
International classification
Abstract
A method is provided for estimating at least one combustion characteristic of a fuel gas belonging to a family of fuel gases, where the at least one characteristic includes at least one of a Wobbe index or a higher heating value. The method includes measuring at least two flow properties of the fuel gas and measuring a dihydrogen content X.sub.H.sub.
using an empirical affine relationship of
Here, α, β, and γ are coefficients predetermined for the family of fuel gases, and Y is a variable representative of physical properties of the fuel gas prepared from the measurements of the at least two flow properties of the fuel gas.
Claims
1. A method of estimating at least one combustion characteristic of a fuel gas belonging to a family of fuel gases, the at least one combustion characteristic comprising at least one of a Wobbe index or a higher heating value, the method comprising: measuring at least two flow properties of the fuel gas, the at least two flow properties comprising two or more of: a flow rate, a temperature, or a pressure; measuring a dihydrogen content X.sub.H.sub.
2. The method according to claim 1, wherein the coefficients α, β, and γ are coefficients read from a chart having as input the measurement of the dihydrogen content X.sub.H.sub.
3. The method according to claim 2, wherein the chart associates different values of the coefficients α, β, and γ with different value ranges for the dihydrogen content X.sub.H.sub.
4. The method according to claim 1, wherein values of the coefficients α, β, and γ are obtained from a dataset relating to known gases of the family of fuel gases for which values of Y and combustion characteristics that are representative of the physical properties are known.
5. The method according to claim 4, further comprising: randomly generating combustion characteristics and values for Y that are representative of the physical properties from the dataset relating to the known gases of the family of fuel gases.
6. The method according to claim 1, wherein the Wobbe index and the higher heating value are estimated using two empirical affine relationships.
7. The method according to claim 6, further comprising: estimating the density of the fuel gas from the estimated Wobbe index and from the estimated higher heating value.
8. The method according to claim 1, further comprising: regulating the at least one combustion characteristic of the fuel gas or regulating the at least one combustion characteristic of the fuel gas and an estimated stoichiometric volume of air or an estimated combustibility index corresponding to the at least one combustion characteristic.
9. The method according to claim 1, wherein: the at least one combustion characteristic of the fuel gas comprises the Wobbe index
10. The method according to claim 1, wherein: the Wobbe index of the fuel gas is estimated; and the method further comprises measuring the density of the fuel gas and estimating the higher heating value from the estimated Wobbe index and from the density of the fuel gas.
11. The method according to claim 1, wherein: the at least one combustion characteristic of the fuel gas comprises the higher heating value
12. A device for estimating at least one combustion characteristic of a fuel gas belonging to a family of fuel gases, the at least one combustion characteristic comprising at least one of a Wobbe index or a higher heating value, the device comprising: at least two first sensors configured to measure at least two flow properties of the fuel gas, the at least two flow properties comprising two or more of: a flow rate, a temperature, or a pressure; a second sensor configured to measure a dihydrogen content X.sub.H.sub.
13. The device according to claim 12, wherein: the at least one combustion characteristic of the fuel gas comprises the Wobbe index
14. The device according to claim 12, wherein: the device is configured to estimate the Wobbe index of the fuel gas; the device further comprises a third sensor configured to measure the density of the fuel gas; and the processor is configured to estimate the higher heating value from the estimated Wobbe index and from the density of the fuel gas.
15. The device according to claim 12, wherein: the at least one combustion characteristic of the fuel gas comprises the higher heating value
16. The device according to claim 12, wherein the processor is also configured to estimate a stoichiometric volume of air or a combustibility index.
17. The device according to claim 12, further comprising: an actuator configured to regulate the at least one combustion characteristic of the fuel gas or to regulate the at least one combustion characteristic of the fuel gas and an estimated stoichiometric volume of air or an estimated combustibility index corresponding to the at least one combustion characteristic.
18. A non-transitory computer-readable medium containing instructions that when executed cause a processor to: estimate at least one combustion characteristic of a fuel gas belonging to a family of fuel gases, the at least one combustion characteristic comprising at least one of a Wobbe index or a higher heating value; wherein, to estimate the at least one combustion characteristic of the fuel gas, the instructions when executed cause the processor to: obtain measurements of at least two flow properties of the fuel gas, the at least two flow properties comprising two or more of: a flow rate, a temperature, or a pressure; obtain a measurement of a dihydrogen content X.sub.H.sub.
19. The non-transitory computer-readable medium of claim 18, wherein: the at least one combustion characteristic of the fuel gas comprises the Wobbe index
20. The non-transitory computer-readable medium of claim 18, wherein: the at least one combustion characteristic of the fuel gas comprises the higher heating value
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) Other characteristics and advantages of the present invention appear from the following description made with reference to the accompanying drawings, which show an example having no limiting character.
(2) In the figures:
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
DETAILED DESCRIPTION
(11) There follows a description of a method and a device for estimating a combustion characteristic of a fuel gas belonging to a family of fuel gases. The characteristic may be the Wobbe index or the higher heating value.
(12) The invention is not limited in any way to estimating a single characteristic, and it may comprise estimating the Wobbe index and the higher heating value simultaneously.
(13)
(14) This method is particularly adapted to gases that might contain a non-zero quantity of dihydrogen, and in particular to fuel gases containing dihydrogen in the range 0% to 20%.
(15) In a first step E01, at least two flow properties of the fuel gas are measured. The fuel gas under study is then flowing through a flow member and flow properties such as temperature, pressure, or indeed flow rate are measured, e.g. by means of sensors.
(16) In a second step E02 the dihydrogen content in the flowing gas is measured, and written X.sub.H.sub.
(17) In a third step E03, the characteristic written
(18)
is estimated by means of the following empirical affine relationship:
(19)
where:
α,β, and γ are coefficients predetermined for the family of fuel gases; and
Y is a variable representative of physical properties of said fuel gas prepared from said measured values of said at least two flow properties of said fuel gas.
(20) The measured values of said at least two flow properties of said fuel gas are selected specifically so that the variable Y can be both representative of the physical properties of the fuel gas, such as viscosity, specific heat capacity (per unit weight), or indeed density, and also capable of being expressed as a function of the measured values of the flow properties of the fuel gas.
(21) Thus, by using the notation:
(22) U.sub.mes,1 for the measurement of the flow property U.sub.1,
(23) U.sub.mes,2 for the measurement of the flow property U.sub.2; and
(24) U.sub.mes,3 for the measurement of the flow property U.sub.3, this gives:
Y=f(U.sub.mes,1;U.sub.mes,2:U.sub.mes,3)
(25)
(26) The device is capable in particular of performing a method of the kind described with reference to
(27) In this example, the device 1 comprises a flow member 2 configured to receive a stream of fuel gas for which it is desired to estimate the Wobbe index or the higher heating value, for example. The flow member has an inlet 2a for receiving the stream of fuel gas, and an outlet 2b.
(28) The device 1 also has two modules 3a and 3b for measuring flow properties of the fuel gas.
(29) By way of example, the modules 3a and 3b may both, or each, measure a flow rate, a pressure, or indeed a temperature.
(30) Specifically, the modules 3a and 3b are conventional sensors for measuring a flow rate, a pressure, or indeed a temperature.
(31) The device 1 includes a module 4 for measuring a dihydrogen content X.sub.H.sub.
(32) The modules 3a, 3b, and 4 are in communication with an estimator module 5 so as to communicate the results of their respective measurements to the estimator module 5. The estimator module 5 may be a computer having a processor and a memory (not shown).
(33) The module 5 is configured to estimate said at least one characteristic
(34)
by means of the following empirical affine relationship:
(35)
where:
α, β, and γ are coefficients predetermined for the family of fuel gases; and
Y is a variable representative of physical properties of said fuel gas prepared from said values of said at least two flow properties of said fuel gas, as measured by the modules 3a and 3b.
(36) To this end, the module 5 may have, stored in memory, possible values for the predetermined coefficients α,β, and γ. The module 5 may also have computer program instructions stored in memory to enable it to execute the estimation.
(37) By way of example, the computer program may comprise instructions for calculating the variable Y from the values measured by the modules 3a and 3b, and instructions for calculating the characteristic
(38)
by means of the above-defined function.
(39) Optionally, the device 1 may include an actuator 6 controlled by the estimator module 5 in order to regulate at least said characteristic. By way of example, the actuator 6 may be a compressed air injector.
(40) With reference to
(41)
of a fuel gas is estimated.
(42)
(43)
of a fuel gas.
(44) This method comprises a calibration first step E11, during which a measurement is taken of a mass flow rate of a reference gas (e.g. methane) in sonic flow through a fluid flow constriction (e.g. an orifice or a micro-nozzle), at a measured reference absolute pressure and at a measured reference absolute temperature.
(45) In a second step E12, a measurement is taken of a mass flow rate of a fuel gas in sonic flow through the fluid flow constriction, this measurement being taken at an absolute pressure measured upstream from the constriction and at an absolute temperature measured upstream from the constriction.
(46) In a step E13, the dihydrogen content X.sub.H.sub.
(47) Steps E11 to E13 may be performed in any possible order. In particular, steps E12 and E13 may be performed simultaneously.
(48) In a step E14, the Wobbe index
(49)
is estimated by an empirical affine relationship that is written:
(50)
with:
(51)
where:
Q.sub.mes,2 is the measured mass flow rate of the fuel gas;
p.sub.mes is the measured absolute pressure of the fuel gas;
T.sub.mes is the measured absolute temperature of the fuel gas;
Q.sub.ref is the measured mass flow rate of the reference gas;
p.sub.ref is the measured absolute pressure of the reference gas; and
T.sub.ref is the measured absolute temperature of the reference gas; and
D, E, and F are predetermined coefficients for the fuel gas family and correspond respectively to the coefficients α, β, and γ.
(52) In order to perform the steps E11 and E12, a fluid flow constriction is used, such as an orifice or a micro-nozzle of known shape.
(53) On either side of the fluid flow constriction, the following measurements are taken: upstream from the fluid flow constriction, the absolute temperatures and the absolute pressures are measured; and downstream from the fluid flow constriction, the mass flow rate is measured (e.g. by means of a thermal mass flow meter).
(54) Specifically, that corresponds to measuring the normal volume flow rate, written Q.sub.vn, passing through the device, in two different ways.
(55) Specifically, since the flow is sonic, the following applies:
(56)
with:
k a constant characteristic of the shape of the orifice;
p and T the absolute pressure and the absolute temperature of the gas upstream from the orifice; and
d the density of the gas.
(57) Also, by using a thermal mass flow meter, the following is obtained:
Q.sub.vn=Q.sub.mes,2.Math.C
with:
Q.sub.vn,mes,2 the measurement taken by the thermal mass flow meter in step E12;
C a correction coefficient that takes account of the differences of the physical properties of the gas (e.g. its specific heat capacity, its viscosity, its thermal conductivity) compared with the properties of air.
(58) For a gas of known composition, the following applies:
(59)
with:
X.sub.i the volume fraction of the component i;
C.sub.i the correction coefficient relating to the component i.
This coefficient can be read from a table associated with the thermal mass flow meter.
(60) The calibration step E11 is a step during which the device is fed with a reference gas of composition that is accurately known (preferably pure methane).
(61) The equality between the two flowrate relationships using the values obtained in the calibration procedure of step E11 can thus be written as follows:
(62)
(63) And, in the measuring procedure of step E12 in which measurements are taken on the fuel gas, the following applies:
(64)
(65) By means of the measurement in step E14, the hydrogen content X.sub.H.sub.
(66) It may be observed that in the above equations, the only two parameters that are unknown are C and density d. It is then possible to eliminate the parameter k in order to express the variable Y, associated with C and with the density d.
(67) Specifically, the following applies:
(68)
(69) And, Y is a variable representative of the physical properties of the fuel gas. It is possible to rewrite Y from the measured values. These measured values (Q.sub.mes,2, p, and T) can be written U.sub.mes,1 for the measured flow property U.sub.1, U.sub.mes,2 for the measured flow property U.sub.2, and U.sub.mes,3 for the measured flow property U.sub.3, thus giving:
Y=f(U.sub.mes,1;U.sub.mes,2;U.sub.mes,3)
(70) Where:
(71) U.sub.mes,1=p.sub.mes;
(72) U.sub.mes,2=Q.sub.mes,2; and
(73) U.sub.mes,3=T.sub.mes
(74) Finally, the following applies:
(75)
(76) This definition of the variable Y can be used in an empirical affine relationship having as its variable Y, and that is written in the following form:
(77)
D, E, and F being predetermined coefficients for the fuel gas family and corresponding respectively to the coefficients α, β, and γ described with reference to
(78) The coefficients D, E, and F can be obtained from known gas compositions, e.g. compositions of natural gas in the network of a country or a region. By way of example, it is possible to use known compositions for higher heating value gas distributed in Europe and well known to the person skilled in the art.
(79) From these known compositions, it is possible to define limit values concerning the molar fractions of the various compounds. For example, by using the notation X.sub.K for the molar fraction of the species K in a gas, it is possible to have known compositions of the following type:
0%<X.sub.N2<5.5%
0.5%<X.sub.C2H6<12.5%
0%<X.sub.C4H10<3.5%
0%<X.sub.CO2<1.3%
0%<X.sub.C3H.sub.8<3.5%
0%<X.sub.C5H12<3.5%
X.sub.CH4=1−ΣX.sub.K
(80) It should be observed that there is no dihydrogen in present-day gas compositions.
(81) These possible ranges enable gas compositions to be generated randomly for which it is possible to determine the Wobbe index (or even the higher heating value). 10,000 gases can thus be generated.
(82) For each randomly generated gas, it is possible to add a random dihydrogen quantity lying in the range 0% to 20% molar fraction. Once more, it is possible to determine the Wobbe index for these randomly generated gases.
(83) It is also possible to deduce from the randomly generated compositions the associated values for the variable Y. By applying a least squares method, it is possible to obtain values for the coefficients D, and F:
(84)
(85) These results were obtained by using methane as the reference gas.
(86) It may be observed that for a dihydrogen molar fraction lying in the range 0% to 20%, the error in estimating the Wobbe index is always less than 1.4%, and that it is less than 1% in 98% of the 10,000 gases that were obtained randomly.
(87) Alternatively, it is possible to use a correlation in which the coefficients are not constants but depend on the dihydrogen content as a molar fraction.
(88) The empirical affine relationship may be rewritten as follows:
(89)
(90) In this example, D and E are functions of the dihydrogen content (the equation has been rewritten so as to show only two coefficients that depend on X.sub.H.sub.
(91) In the same manner, by generating random gas compositions, the inventors have observed that it is possible to use functions D and that are piecewise constant over dihydrogen concentration ranges of width 1%.
(92) For example, in the range 2% to 3% dihydrogen, D and E take the following values:
(93) D=−18.181
(94) E=32.75
(95) It has been observed that the error in estimating the Wobbe index is then always less than 1% for 95% of the 10,000 random gases.
(96)
(97) This device 10 has an inlet 11 for receiving a fuel gas G1 and an inlet 12 for receiving a reference gas G2 (typically methane).
(98) The device 10 also has a selector and guide module for bringing the fuel gas stream or the reference gas stream to the inlet of a pipe 13. Specifically, the selector and guide module comprises a valve 31 and a valve 35.
(99) The valves 31 and 35 are controlled by a module 36 via respective terminals 52 and 51 so as to cause either the fuel gas or else the reference gas to flow into the pipe.
(100) In this example, the pipe comprises, from upstream to downstream starting from its inlet 13: An absolute pressure sensor 60 connected to the module 36 via its terminal 41; An absolute temperature sensor 50 connected to the module 36 via its terminal 42; A fluid flow constriction 32 (e.g. an orifice or a micro-nozzle); A thermal mass flow rate sensor 33 connected to the module 36 via its terminal 43; and A sensor 34 for sensing the dihydrogen molar fraction and connected to the module 36 via its terminal 44.
(101) The module 36 can apply an empirical affine relationship such as those described above in order to estimate the Wobbe index from the signals received at the terminals 41 to 44.
(102) In the example shown, the module 36 communicates via two output terminals 54 and 53 with an external module 37, e.g. an actuator for performing regulation, or indeed a display. In a variant that is not shown, the module 37 is integrated in the device 10.
(103) With reference to
(104)
(105)
of a fuel gas.
(106) In a step E22, the mass flow rate of said fuel gas is measured in laminar flow through an apparatus giving rise to a pressure drop, this measurement depending on the viscosity of the fuel gas and on the viscosity of a reference gas.
(107) In a step E23, downstream from said apparatus giving rise to a pressure drop, the mass flow rate of said fuel gas is measured by means of a thermal mass flow meter, the measurement depending on the specific heat capacity of the fuel gas and on the heat capacity of the reference gas.
(108) In a step E21, which may be performed before or after performing the steps E22 and E23, a calibration step is performed that corresponds to performing the steps E22 and E23, but with a reference gas (e.g. methane).
(109) A step E24 is also performed of measuring the dihydrogen content X.sub.H.sub.
(110) Finally, in a step E25, the higher heating value is estimated by means of an empirical affine relationship having the form:
(111)
with:
(112)
where:
Z is a variable corresponding to the variable Y described with reference to
Q.sub.mes,1 is the mass flow rate of the fuel gas in laminar flow through an apparatus giving rise to a measured pressure drop; and
Q.sub.mes,2 is the mass flow rate of the fuel gas measured downstream from said apparatus giving rise to a pressure drop; and
A, B, and C are predetermined coefficients for the fuel gas family and correspond respectively to the coefficients α,β, and γ.
(113) The following relationship for higher heating value (HHV) is known from prior document DE 4 118 781:
(114)
with:
α and β predetermined constants;
ρ the density of the fuel gas;
μ the viscosity of the fuel gas; and
C.sub.P the heat capacity of the fuel gas.
(115) The following dimensionless variable Z can be defined as follows:
(116)
(117) Here, the subscript GN specifies the fuel gas, and ref specifies a reference gas.
(118) The relationship concerning the higher heating value can be rewritten as follows:
HHV.sub.GN=A.Math.Z+B
(119) In order to measure Z, the two flow rate measurements are used. By way of example, Q.sub.mes,1 is the mass flow rate of the fuel gas in laminar flow through an apparatus giving rise to a measured pressure drop. By using Poiseuille's law, it is known that the apparatus subjects the gas stream to a pressure drop through a laminar element.
(120) Since the measurement depends on viscosity, for certain measuring apparatuses, it has a form that depends on the upstream-to-downstream pressure difference Δp.sub.mes and the flowmeter that delivers Q.sub.mes,1 can determine the volume flow rate by means of the following formula:
(121)
(122) Where K.sub.1 is a geometrical constant.
(123) In this example, a pressure sensor and a temperature probe are integrated in the flowmeter and they make it possible to obtain directly the mass flow rate Q.sub.mes,1 (or volume flow rate reduced to normal temperature and pressure conditions), as follows:
(124)
(125) The measurement of the corrected volume flow rate Q.sub.mes,1 depends on the viscosity of the fuel gas.
(126) The flow rate Q.sub.mes,2 is measured using a thermal mass flow meter. The measurement of the mass flow rate Q.sub.mes,2 depends on the specific heat capacity of the fuel gas.
(127) In this way, two estimates are obtained of the mass flow rate, each presenting a respective error relative to the real normal volume flow rate (Q.sub.vn). The following applies:
(128)
(129) For Q.sub.mes,1, the error is associated with the difference between the real viscosity of the mixture gas and the viscosity of the reference gas (methane in this example).
(130) For Q.sub.mes,2, the error is associated with the difference between the density multiplied by the real specific heat capacity of the mixture gas, and by the same quantity for the reference gas (methane in this example).
(131) By taking the ratio of the two measured flow rates, the variable Z is determined and the correlation relationship of Document DE 4 118 781 becomes:
(132)
(133) By writing:
(134) U.sub.mes,1=Q.sub.mes,1; and
(135) U.sub.mes,2=Q.sub.mes,2
(136) The following applies:
(137)
(138) Nevertheless, the inventors have observed that this relationship is not applicable once dihydrogen is present in the fuel gas. Specifically, the measurement error becomes too large, dihydrogen has properties (viscosity, specific heat capacity) that are very different from those of gases such as alkanes.
(139) In order to mitigate that drawback, use is made of the measured dihydrogen content X.sub.H.sub.
(140)
A, B, and C being predetermined coefficients for the fuel gas family and corresponding respectively to the coefficients α, β, and γ described with reference to
(141) The coefficients A, B, and C can be obtained from known gas compositions, e.g. compositions of natural gas in the network of a country or a region. By way of example, it is possible to use known compositions for high heating value gas distributed in Europe and well known to the person skilled in the art.
(142) From these known compositions, it is possible to define limit values concerning the molar fractions of the various compounds. For example, by using the notation XK for the molar fraction of the species K in a gas, it is possible to have known compositions of the following type:
0%<X.sub.N2<5.5%
0.5%<X.sub.C2H6<12.5%
0%<X.sub.C4H10<3.5%
0%<X.sub.CO2<1.3%
0%<X.sub.C3H.sub.8<3.5%
0%<X.sub.C5H12<3.5%
X.sub.CH4=1−ΣX.sub.K
It should be observed that there is no dihydrogen in present-day gas compositions.
(143) These possible ranges make it possible to generate random gas compositions for which it is possible to determine the higher heating value. 10,000 gases can thus be generated.
(144) For each randomly generated gas, it is possible to add a random dihydrogen quantity lying in the range 0% to 20% molar fraction. Here likewise, it is possible to determine the higher heating value of these randomly-generated gases.
(145) It is also possible to deduce from the randomly-generated compositions the associated values for the variable Z. By applying a least squares method, it is possible to obtain values for the coefficients A, B, and C:
A=0.61650996
B=10.428
C=−0.0645996
(146) Alternatively, it is possible to use a correlation in which the coefficients are not constants but depend on the dihydrogen content as a molar fraction.
(147) The empirical affine relationship may be rewritten as follows:
(148)
(149) In this example, A and B are functions of the dihydrogen content (the equation has been rewritten so as to show only two coefficients that depend on X.sub.H.sub.
(150) In the same manner, by generating random gas compositions, the inventors have observed that it is possible to use functions A and B that are piecewise constant over dihydrogen concentration ranges of width 1%.
(151) For example, in the range 19% to 20% dihydrogen, A and B take the following values:
(152) A=0.397217
(153) B=10.189715
(154)
(155) This device 100 has an inlet 111 for receiving a fuel gas G1 and an inlet 112 for receiving a reference gas G2 (typically methane).
(156) The device 100 also has a selector and guide module for bringing the fuel gas stream or the reference gas stream to the inlet of a pipe 113. Specifically, the selector and guide module comprises a valve 131 and a valve 135.
(157) The valves 131 and 135 are controlled by a module 136 via respective terminals 152 and 151 so as to cause either the fuel gas or else the reference gas to flow into the pipe.
(158) In this example, the pipe comprises, from upstream to downstream and starting from its inlet 113: A sensor 132 for sensing mass flow rate in laminar flow through an apparatus that gives rise to a pressure drop, the measurement depending on the viscosity of the fuel gas and on the viscosity of the reference gas, the sensor being connected to the module 136 via its terminal 142; A sensor 133 for sensing thermal mass flow rate, the measurement depending on the specific heat capacity of the measured gas and on the specific heat capacity of the reference gas, the sensor being connected to the module 136 via its terminal 143; and A sensor 134 for sensing the dihydrogen molar fraction and connected to the module 136 via its terminal 144.
(159) The module 136 can apply an empirical affine relationship such as those described above in order to estimate the higher heating value from the signals received at the terminals 142 to 144.
(160) In the example shown, the module 136 communicates via two output terminals 154 and 153 with an external module 137, e.g. an actuator for performing regulation, or indeed a display. In a variant that is not shown, the module 137 is integrated in the device 100.
(161)
(162) This device has two inlets for gas under pressure. A first inlet 1011 is for receiving a fuel gas G1 belonging to a family of gases (e.g. the fuel gas in the second family specified in French standard NF EN 437) and also containing dihydrogen. The device 1000 has another inlet 1012 for a reference gas G0 used for the calibration stages.
(163) The device 1000 also has a selector and guide module for bringing the stream of the fuel gas G1 or of the reference gas G0 to the inlet of a pipe 1013. The selector and guide module comprises a valve 1031 and a valve 1035.
(164) The valves 1031 and 1035 are controlled by a module 1036 via respective terminals 1052 and 1051 so as to cause either the fuel gas or else the reference gas to flow into the pipe.
(165) Starting from the inlet of the pipe 1013, and going from upstream to downstream, there are to be found: An absolute pressure sensor 1006 connected to the module 1036 via its terminal 1041. This sensor delivers the value P.sub.mes measured on the fuel gas and the value P.sub.ref measured on the reference gas. An absolute temperature sensor 1005 connected to the module 1036 via its terminal 1042. This sensor delivers the value T.sub.mes measured on the fuel gas and the value T.sub.ref measured on the reference gas. A fluid flow constriction 1032 (e.g. an orifice or a micro-nozzle); A sensor for sensing of the mass flow rate through a laminar pressure drop connected to the module 1036 via its terminal 1043. This sensor delivers the value Q.sub.mes,1. A thermal mass flow rate sensor 1033 connected to the module 1036 via its terminal 1046. This sensor delivers the value Q.sub.mes,2 measured on the fuel gas and the value Q.sub.ref measured on the reference gas. A sensor 1034 for sensing the dihydrogen molar fraction and connected to the module 36 via its terminal 1047. A vent (39).
(166) In this example, the module 1036 can calculate the following two variables:
(167)
(168) The module 1036 can then use the following formulae for determining the Wobbe index
(169)
the higher heating value
(170)
the stoichiometric volume of air
(171)
combustibility index
(172)
and the density of the fuel gas mixture
(173)
(174) Specifically, it is possible to determine the stoichiometric volume of air Va, the combustibility index B, and the density of the fuel gas mixture from the values for the Wobbe index IW and for the higher heating value HHV.
(175) More precisely, for the natural gases usually distributed in Europe (and presently not containing any dihydrogen, the following applies:
(176)
(177) The inventors have observed that in the presence of dihydrogen, the values may depart from 1.162.
(178) It is possible to associate the ratio of the higher heating value to the stoichiometric volume of air to the (known) dihydrogen content, and to do so with an error of less than 0.3% for all the natural gases that have been studied.
(179) In other words, knowing the dihydrogen content and the estimated higher heating value with an error of less than 1% makes it possible to know the stoichiometric volume of air and thus to be able to regulate combustion with open-loop regulation.
(180) In the same manner, the ratio of the Wobbe index to the combustibility index varies linearly with the dihydrogen content. This also makes it possible to perform open-loop regulation.
(181) It may be observed that the error in estimating the combustibility index is less than 1% in 99.5% of the 10,000 gases of random composition, if 10,000 gases are used.
(182) By definition of the Wobbe index, it is possible to determine the density of the gas, providing both the Wobbe index and the higher heating value are known.
(183) It may be observed that the device 1000 delivers a Wobbe index signal at its terminal 1054, a higher heating value signal at its terminal 1055, a density signal at its terminal 1056, a stoichiometric volume of air signal at its terminal 1057, and a combustibility index signal at its terminal 1058. These signals relate to the fuel gas under study, i.e. to the gas G1.
(184)
(185) This apparatus is capable in particular of continuously regulating the Wobbe index for gases analogous to those distributed in Europe, plus an additional quantity of dihydrogen. In particular, it is because the measurement is continuous or in real time that it becomes possible to perform regulation, whereas that is not possible with apparatuses of the chromatograph type.
(186) The accuracy with which the Wobbe index is measured is of the order of 1%.
(187) It may be observed that it is possible to modify a setpoint value for the Wobbe index in the device 10 in order to obtain the desired regulation.
(188) Furthermore, the calibration steps can be performed automatically while using methane. Starting stages can also be performed automatically, i.e. without operator intervention.
(189) In addition, combined regulation is possible in which closed-loop regulation is performed on one of the characteristics (e.g. the Wobbe index) and open-loop regulation is performed on a flow rate of compressed air to be injected. This makes it possible to take account in particularly accurate manner of variations in the characteristic in order to satisfy more closely a setpoint for the characteristic.
(190) It may be observed that such apparatus may consume natural gas at a rate of less than 150 liters per hour (L/h).
(191) The implementations and embodiments described above enable thermodynamic data (Cp, viscosity) to be coupled with a measurement of dihydrogen content in volume or in molar fraction. Using an appropriate correlation makes it possible to calculate the Wobbe index and/or the higher heating value.
(192) It may be observed that compared with using a gas phase chromatograph, the following applies: estimation is almost instantaneous (with a response time that may be less than 5 seconds (s)); estimation is inexpensive.
(193) Compared with existing correlation methods, good accuracy is obtained even though dihydrogen is present, together with good measurement robustness.
(194) Finally, compared with apparatuses making use of combustion (combustibility meter or calorimeter), the invention provides good robustness, ease of implementation, and maintenance that is reduced and easy to perform.