JP2012202739A - Heating value measuring system and heating value measuring method - Google Patents

Abstract

PROBLEM TO BE SOLVED: To provide a heating value measuring system capable of easily and accurately measuring a heating value of gas.SOLUTION: A heating value measuring system comprises: a container 101 which has a heater element, where a plurality of different voltages can be applied thereto, as well as a pressure sensor and is injected with gas; a formula storage device 402 which stores a heating value calculation formula with electric signals from the heater element, when energized by each of the plurality of the voltages, as well as the electric signal from the pressure sensor as independent variables and the heating value as a dependent variable; and a heating value calculation section 305 which calculates a measurement of the heating value of the gas injected into the container 101 by substituting the electric signals from the heater element as well as the pressure sensor for the independent variables in the heating value calculation formula.

Description

  The present invention relates to a gas inspection technique, and relates to a calorific value measuring system and a calorific value measuring method.

When determining the calorific value of natural gas, it is necessary to analyze the component ratio of natural gas using an expensive gas chromatography device or the like. In addition, by measuring the thermal conductivity of natural gas and the speed of sound in a natural gas atmosphere, methane (CH ₄ ), propane (C ₃ H ₈ ), nitrogen (N ₂ ), and carbon dioxide (CO ₂ ) contained in natural gas are measured. A method of calculating the component ratio of ₂ ) and obtaining the calorific value of natural gas has also been proposed (see, for example, Patent Document 1).

JP-T-2004-514138

  However, the method disclosed in Patent Document 1 requires an expensive sound speed sensor for measuring sound speed in addition to the sensor for measuring thermal conductivity. Therefore, an object of the present invention is to provide a calorific value measuring system and a calorific value measuring method capable of easily and accurately measuring the calorific value of a gas.

  According to the aspect of the present invention, (a) a heating element to which a plurality of different voltages can be applied, and a pressure sensor are provided, and a container into which gas is injected, and (b) a case where a plurality of different voltages are respectively applied. An equation storage device for storing a calorific value calculation formula having the electric signal from the heating element and the electric signal from the pressure sensor as independent variables and the calorific value as a dependent variable, and (c) an independent variable of the calorific value calculation formula, A calorific value measurement system is provided, comprising: a calorific value calculation unit that calculates a calorific value of gas injected into a container by substituting measured values of electrical signals from the heating element and the pressure sensor.

  In addition, according to the aspect of the present invention, (a) injecting a gas into a container provided with a heating element and a pressure sensor to which a plurality of different voltages can be applied, and (b) a plurality of different voltages can be respectively applied. Preparing a calorific value calculation formula with the electric signal from the heating element and the electric signal from the pressure sensor as independent variables and the calorific value as a dependent variable, and (c) as an independent variable of the calorific value calculation formula Substituting the measured value of the electrical signal from the heating element and the pressure sensor to calculate the measured value of the calorific value of the gas injected into the container, a method of measuring the calorific value is provided.

  Moreover, according to the aspect of the present invention, (a) a measurement unit that measures a measured value of the heat dissipation coefficient or thermal conductivity of gas and a measured value of the pressure of gas, and (b) a heat dissipation coefficient or thermal conductivity. And (c) a measured value of the heat dissipation coefficient or thermal conductivity of the gas, a measured value of the pressure of the gas, and a correlation. And a calorific value calculation unit for calculating a measured value of the calorific value of the gas.

  Moreover, according to the aspect of the present invention, (a) a measurement value of the heat dissipation coefficient or thermal conductivity of the gas and a measurement value of the gas pressure are measured, and (b) the heat dissipation coefficient or the thermal conductivity. Providing a correlation between the pressure and the calorific value, and (c) a measured value of the heat dissipation coefficient or thermal conductivity of the gas, a measured value of the pressure of the gas, and a correlation. And a calorific value measurement method comprising: calculating a calorific value measurement value of the calorific value.

  ADVANTAGE OF THE INVENTION According to this invention, the calorific value measuring system and calorific value measuring method which can measure the calorific value of gas easily and correctly can be provided.

1 is a perspective view of a microchip according to an embodiment of the present invention. It is sectional drawing seen from the II-II direction of FIG. 1 of the microchip which concerns on embodiment of this invention. It is a circuit diagram regarding the heat generating element which concerns on embodiment of this invention. It is a circuit diagram regarding the temperature measuring element which concerns on embodiment of this invention. It is a graph which shows the relationship between the temperature of the heat generating element which concerns on embodiment of this invention, and the thermal radiation coefficient of gas. It is the 1st schematic diagram of the calorific value measuring system concerning an embodiment of the invention. It is the 2nd schematic diagram of the calorific value measuring system concerning an embodiment of the invention. It is a flowchart which shows the preparation method of the emitted-heat amount calculation formula which concerns on embodiment of this invention. It is a flowchart which shows the measuring method of the emitted-heat amount which concerns on embodiment of this invention. It is a graph which shows the relationship between the thermal conductivity of the gas which concerns on other embodiment of this invention, and a thermal radiation coefficient. It is a graph which shows the error from the true value of the calorific value calculated of the sample mixed gas concerning the example of the embodiment of the invention. It is a graph which shows the difference | error from the true value of the calorific value calculated of the sample mixed gas which concerns on the comparative example of embodiment of this invention.

  Embodiments of the present invention will be described below. In the following description of the drawings, the same or similar parts are denoted by the same or similar reference numerals. However, the drawings are schematic. Therefore, specific dimensions and the like should be determined in light of the following description. Moreover, it is a matter of course that portions having different dimensional relationships and ratios are included between the drawings.

  First, the microchip 8 used in the calorific value measurement system according to the embodiment will be described with reference to FIG. 1 which is a perspective view and FIG. 2 which is a cross-sectional view seen from the II-II direction. The microchip 8 includes a substrate 60 provided with a cavity 66, and an insulating film 65 disposed on the substrate 60 so as to cover the cavity 66. The thickness of the substrate 60 is, for example, 0.5 mm. The vertical and horizontal dimensions of the substrate 60 are, for example, about 1.5 mm. A portion of the insulating film 65 covering the cavity 66 forms a heat insulating diaphragm. Furthermore, the microchip 8 includes a heating element 61 provided in the diaphragm portion of the insulating film 65, a first temperature measuring element 62 and a second temperature measuring element 62 provided in the diaphragm portion of the insulating film 65 so as to sandwich the heating element 61. The temperature measuring element 63 and a heat retaining element 64 provided on the substrate 60 are provided.

  The heating element 61 is disposed at the center of the diaphragm portion of the insulating film 65 that covers the cavity 66. The heating element 61 is a resistor, for example, and generates heat when power is applied to heat the atmospheric gas in contact with the heating element 61. The first temperature measuring element 62 and the second temperature measuring element 63 are electronic elements such as passive elements such as resistors, for example, and output electric signals depending on the gas temperature of the atmospheric gas. In the following, an example in which the output signal of the first temperature measuring element 62 is used will be described. However, the present invention is not limited to this. For example, the output signal of the first temperature measuring element 62 and the output of the second temperature measuring element 63 are described. You may utilize the average value of a signal as an output signal of a temperature measuring element.

The heat retaining element 64 is a resistor, for example, and generates heat when supplied with electric power, and keeps the temperature of the substrate 60 constant. As a material of the substrate 60, silicon (Si) or the like can be used. As a material of the insulating film 65, silicon oxide (SiO ₂ ) or the like can be used. The cavity 66 is formed by anisotropic etching or the like. Further, platinum (Pt) or the like can be used as the material of the heating element 61, the first temperature measuring element 62, the second temperature measuring element 63, and the heat retaining element 64, and can be formed by a lithography method or the like. . Further, the heating element 61, the first temperature measuring element 62, and the second temperature measuring element 63 may be made of the same member.

  The microchip 8 is fixed to a container such as a chamber filled with atmospheric gas via a heat insulating member 18 disposed on the bottom surface of the microchip 8. By fixing the microchip 8 to the chamber or the like via the heat insulating member 18, the temperature of the microchip 8 becomes less susceptible to the temperature fluctuation of the inner wall of the chamber or the like. The heat conductivity of the heat insulating member 18 made of glass or the like is, for example, 1.0 W / (m · K) or less.

As shown in FIG. 3, for example, a positive input terminal of an operational amplifier 170 is electrically connected to one end of the heating element 61, and the other end is grounded. A resistance element 161 is connected in parallel with the + input terminal and the output terminal of the operational amplifier 170. The negative input terminal of the operational amplifier 170 is a power supply, between the resistor element 162 and the resistor element 163 connected in series, between the resistor element 163 and the resistor element 164 connected in series, and the resistor element 164 connected in series. And the resistance element 165 or electrically connected to the ground terminal of the resistance element 165. By appropriately determining the resistance value of each of the resistance elements 162 to 165, for example, when a voltage Vin of 3.2 V is applied to one end of the resistance element 162, there is, for example, 2.8 V between the resistance element 163 and the resistance element 162. Voltage V _L3 is generated. Further, for example, a voltage V _{L2 of} 2.2 V is generated between the resistance element 164 and the resistance element 163, and a voltage V _{L1 of} 1.5 V, for example, is generated between the resistance element 165 and the resistance element 164.

  A switch SW1 is provided between the power supply and the negative input terminal of the operational amplifier, and a switch SW2 is provided between the resistive element 162 and the resistive element 163 and between the negative input terminal of the operational amplifier. A switch SW3 is provided between the resistance element 163 and the resistance element 164 and between the negative input terminal of the operational amplifier. Further, a switch SW4 is provided between the resistance element 164 and the resistance element 165 and between the negative input terminal of the operational amplifier, and between the ground terminal of the resistive element 165 and the negative input terminal of the operational amplifier. , A switch SW5 is provided.

When a voltage V _in of 3.2 V is applied to the negative input terminal of the operational amplifier 170, only the switch SW1 is energized and the switches SW2, SW3, SW4, SW5 are disconnected. Of the operational amplifier 170 - When applying the voltage V _L3 of 2.8V to the input terminal, only the switch SW2 is turned on and the switches SW1, SW3, SW4, SW5 are disconnected. When a voltage V _L2 of 2.2 V is applied to the negative input terminal of the operational amplifier 170, only the switch SW3 is energized and the switches SW1, SW2, SW4, SW5 are disconnected. When a voltage V _L1 of 1.5 V is applied to the negative input terminal of the operational amplifier 170, only the switch SW4 is energized and the switches SW1, SW2, SW3, SW5 are disconnected. When a voltage V _L0 of 0 V is applied to the negative input terminal of the operational amplifier 170, only the switch SW5 is energized and the switches SW1, SW2, SW3, SW4 are disconnected. Therefore, either 0V or four levels of voltages can be applied to the negative input terminal of the operational amplifier 170 by opening and closing the switches SW1, SW2, SW3, SW4, and SW5. Therefore, the applied voltage that determines the heat generation temperature of the heat generating element 61 can be set in four stages by opening and closing the switches SW1, SW2, SW3, SW4, and SW5.

Here, the temperature of the heating element 61 when a voltage V _L1 of 1.5 V is applied to the + input terminal of the operational amplifier 170 is T _H1 . Further, the temperature of the heating element 61 when a voltage V _L2 of 2.2 V is applied to the + input terminal of the operational amplifier 170 is T _H2 , and the heat generation when a voltage V _L3 of 2.8 V is applied to the + input terminal of the operational amplifier 170. The temperature of the element 61 is T _H3 .

  As shown in FIG. 4, for example, a negative input terminal of an operational amplifier 270 is electrically connected to one end of the first temperature measuring element 62, and the other end is grounded. In addition, a resistance element 261 is connected in parallel with the negative input terminal and the output terminal of the operational amplifier 270. The positive input terminal of the operational amplifier 270 is electrically connected between the resistance element 264 and the resistance element 265 connected in series. Thereby, a weak voltage of about 0.3 V is applied to the first temperature measuring element 62.

The resistance value of the heating element 61 shown in FIGS. 1 and 2 varies depending on the temperature of the heating element 61. The relationship between the temperature _TH of the heating element 61 and the resistance value R _H of the heating element 61 is given by the following equation (1).
R _H = R _H _ _STD × [1 + α _H (T _H -T _H _ _STD ) + β _H (T _H -T _H _ _STD ) ² ] ・ ・ ・ (1)
Here, T _{H — STD} represents the standard temperature of the heating element 61, and is 20 ° C., for example. R _{H — STD} represents the resistance value of the heating element 61 measured in advance at the standard temperature T _{H — STD} . α _H represents a first-order resistance temperature coefficient. β _H represents a second-order resistance temperature coefficient.

The resistance value R _H of the heating element 61 is given by the following equation (2) from the drive power P _H of the heating element 61 and the energization current I _H of the heating element 61.
R _H = P _H / I _H ² (2)
Alternatively, the resistance value R _H of the heating element 61 is given by the following equation (3) from the voltage V _H applied to the heating element 61 and the energization current I _H of the heating element 61.
R _H = V _H / I _H (3)

Here, the temperature _TH of the heating element 61 is stabilized when the heating element 61 and the ambient gas are in thermal equilibrium. The thermally balanced state refers to a state in which the heat generation of the heating element 61 and the heat dissipation from the heating element 61 to the atmospheric gas are balanced. (4) below, as shown in the expression, the driving power P _H of the heater element 61 at equilibrium, by dividing the difference [Delta] T _H between the temperature T _I of the temperature T _H and the ambient gas of the heater element 61, the atmospheric gas A heat dissipation coefficient M _I is obtained. The unit of the heat dissipation coefficient M _I is, for example, W / ° C.
M _I = P _H / (T _H -T _I )
= P _H / ΔT _H (4)

From the above equation (1), the temperature T _H of the heater element 61 is given by the following equation (5).
_{T H = (1 / 2β H} ) × [-α H + [α H 2 - 4β H (1 - R H / R H _ STD)] 1/2] + T H _ STD ···(Five)
Therefore, the difference ΔT _H between the temperature T _H of the heating element 61 and the temperature T _{I of the} atmospheric gas is given by the following equation (6).
_{ΔT H = (1 / 2β H} ) × [-α H + [α H 2 - 4β H (1 - R H / R H _ STD)] 1/2] + T H _ STD - T I ··· ( 6)

The temperature T _{I of the} atmospheric gas approximates the temperature T _I of the first temperature measuring element 62 to which power that does not generate heat is given. The relationship between the temperature T _I of the first temperature measuring element 62 and the resistance value R _{I of} the first temperature measuring element 62 is given by the following equation (7).
R _I = R _I _ _STD × [1 + α _I (T _I -T _I _ _STD ) + β _I (T _I -T _I _ _STD ) ² ] (7)
T _{I — STD} represents the standard temperature of the first temperature measuring element 62 and is, for example, 20 ° C. R _{I — STD} represents the resistance value of the first temperature measuring element 62 measured in advance at the standard temperature T _{I — STD} . α _I represents the first-order resistance temperature coefficient. β _I represents a second-order resistance temperature coefficient. From the above equation (7), the temperature T _I of the first temperature measuring element 62 is given by the following equation (8).
_{T I = (1 / 2β I} ) × [-α I + [α I 2 - 4β I (1 - R I / R I _ STD)] 1/2] + T I _ STD ... (8)

Therefore, the radiation coefficient M _I of the atmospheric gas is given by the following equation (9).
M _I = P _H / ΔT _{H
= P H / [(1 /} 2β H) [- α H + [α H 2 -4β H (1-R H / R H _ STD)] 1/2] + T H _ STD - (1 / 2β I ) [-α _I + [α _I ² -4β _I (1-R _I / R _I _ _STD )] ^1/2 ] -T _I _ _STD ]
^{_{= (R H I H 2)}} / [(1 / 2β H) [- α H + [α H 2 -4β H (1-R H / R H _ STD)] 1/2] + T H _ STD - (1 / 2β _I ) [-α _I + [α _I ² -4β _I (1-R _I / R _I _ _STD )] ^1/2 ] -T _I _ _STD ] ... (9)

Since the energizing current I _H and the driving power P _H or voltage V _H of the heating element 61 can be measured, the resistance value R _H of the heating element 61 can be calculated from the above equation (2) or (3). Similarly, the resistance value R _{I of} the first temperature measuring element 62 can also be calculated. Therefore, by using the microchip 8, the heat radiation coefficient M _{I of the} atmospheric gas can be calculated from the above equation (9).

Note that, by keeping the temperature of the substrate 60 constant by the heat retaining element 64, the temperature of the ambient gas in the vicinity of the microchip 8 before the heat generating element 61 generates heat approximates the constant temperature of the substrate 60. For this reason, fluctuations in the temperature of the atmospheric gas before the heat generating element 61 generates heat are suppressed. By temperature variation further heated by the heating element 61 to the atmosphere gas is once suppressed, it is possible to calculate the radiation coefficient M _I with high accuracy.

Here, it is assumed that the atmospheric gas is a mixed gas, and the mixed gas is composed of four types of gas components: gas A, gas B, gas C, and gas D. The sum of the volume ratio V _A of gas _A , the volume ratio V _B of gas _B , the volume ratio V _C of gas _C , and the volume ratio V _D of gas D is 1, as given by the following equation (10). .
V _A + V _B + V _C + V _D = 1 (10)

The calorific value per unit volume of gas A is K _A , the calorific value per unit volume of gas B is K _B , the calorific value per unit volume of gas C is K _C , and the calorific value per unit volume of gas D is Is K _D , the calorific value Q per unit volume of the mixed gas is given by the sum of the volume ratio of each gas component multiplied by the calorific value per unit volume of each gas component. Therefore, the calorific value Q per unit volume of the mixed gas is given by the following equation (11). The unit of the calorific value per unit volume is, for example, MJ / m ³ .
Q = K _A × V _A + K _B × V _B + K _C × V _C + K _D × V _D ... (11)

Moreover, the radiation coefficient M _A gas A, the radiation coefficient of gas B M _B, when the radiation coefficient of gas C M _C, the radiation coefficient of the gas D and M _D, the radiation coefficient M _I of the mixed gas, the It is given as the sum of the volume fraction of the gas component multiplied by the heat dissipation coefficient of each gas component. Therefore, the heat dissipation coefficient M _I of the mixed gas is given by the following equation (12).
M _I = M _A × V _A + M _B × V _B + M _C × V _C + M _D × V _D ... (12)

Further, since the radiation coefficient of gas depends on the temperature T _H of the heater element 61, the radiation coefficient M _I of the mixed gas, as a function of the temperature T _H of the heater element 61 is given by the following equation (13).
M _I (T _H ) = M _A (T _H ) × V _A + M _B (T _H ) × V _B + M _C (T _H ) × V _C + M _D (T _H ) × V _D・ ・ ・ ( 13)

Therefore, the heat release coefficient M _I1 (T _H1 ) of the mixed gas when the temperature of the heating element 61 is T _H1 is given by the following equation (14). The heat dissipation coefficient M _I2 (T _H2 ) of the mixed gas when the temperature of the heating element 61 is T _H2 is given by the following equation (15), and the heat dissipation coefficient of the mixed gas when the temperature of the heating element 61 is T _H3 : M _I3 (T _H3 ) is given by the following equation (16).
M _I1 (T _H1 ) = M _A (T _H1 ) × V _A + M _B (T _H1 ) × V _B + M _C (T _H1 ) × V _C + M _D (T _H1 ) × V _D・ ・ ・ ( 14)
M _I2 (T _H2 ) = M _A (T _H2 ) × V _A + M _B (T _H2 ) × V _B + M _C (T _H2 ) × V _C + M _D (T _H2 ) × V _D・ ・ ・ ( 15)
M _I3 (T _H3 ) = M _A (T _H3 ) × V _A + M _B (T _H3 ) × V _B + M _C (T _H3 ) × V _C + M _D (T _H3 ) × V _D・ ・ ・ ( 16)

Here, the heat dissipation coefficients M _A (T _H ), M _B (T _H ), M _C (T _H ), and M _D (T _H ) of each gas component are non-linear with respect to the temperature T _H of the heating element 61. If so, the above equations (14) to (16) have a linearly independent relationship. Further, the heat dissipation coefficients M _A (T _H ), M _B (T _H ), M _C (T _H ), and M _D (T _H ) of each gas component have linearity with respect to the temperature T _H of the heating element 61. Even in the case where the rate of change of the heat radiation coefficient M _A (T _H ), M _B (T _H ), M _C (T _H ), M _D (T _H ) of each gas component with respect to the temperature T _H of the heating element 61 is different. The above formulas (14) to (16) have a linearly independent relationship. Further, when the equations (14) to (16) have a linearly independent relationship, the equations (10) and (14) to (16) have a linearly independent relationship.

FIG. 5 shows the heat release coefficients of methane (CH ₄ ), propane (C ₃ H ₈ ), nitrogen (N ₂ ), and carbon dioxide (CO ₂ ) contained in natural gas, and the heating element 61 that is a heating resistor. It is a graph which shows the relationship between temperature. The heat release coefficient of each gas component of methane (CH ₄ ), propane (C ₃ H ₈ ), nitrogen (N ₂ ), and carbon dioxide (CO ₂ ) has linearity with respect to the temperature of the heating element 61. However, the rate of change of the heat dissipation coefficient with respect to the temperature of the heating element 61 is different for each of methane (CH ₄ ), propane (C ₃ H ₈ ), nitrogen (N ₂ ), and carbon dioxide (CO ₂ ). Therefore, when the gas components constituting the mixed gas are methane (CH ₄ ), propane (C ₃ H ₈ ), nitrogen (N ₂ ), and carbon dioxide (CO ₂ ), the above (14) to (16 ) Has a linearly independent relationship.

The heat release coefficients M _A (T _H1 ), M _B (T _H1 ), M _C (T _H1 ), M _D (T _H1 ), M _A (T _H2 ) of each gas component in the equations (14) to (16) , M _B (T _H2 ), M _C (T _H2 ), M _D (T _H2 ), M _A (T _H3 ), M _B (T _H3 ), M _C (T _H3 ), M _D (T _H3 ) The value can be obtained in advance by measurement or the like. Therefore, when the simultaneous equations (10) and (14) to (16) are solved, the volume ratio V _{A of} the gas _A , the volume ratio V _{B of} the gas _B , the volume ratio V _{C of} the gas _C , and the volume of the gas D Each of the rates V _D is given as a function of the heat release coefficient M _I1 (T _H1 ), M _I2 (T _H2 ), M _I3 (T _H3 ) of the mixed gas as shown in the following equations (17) to (20). It is done. In the following equations (17) to (20), n is a natural number and f _n is a symbol representing a function.
V _A = f ₁ [M _I1 (T _H1 ), M _I2 (T _H2 ), M _I3 (T _H3 )] (17)
V _B = f ₂ [M _I1 (T _H1 ), M _I2 (T _H2 ), M _I3 (T _H3 )] (18)
V _C = f ₃ [M _I1 (T _H1 ), M _I2 (T _H2 ), M _I3 (T _H3 )] (19)
V _D = f ₄ [M _I1 (T _H1 ), M _I2 (T _H2 ), M _I3 (T _H3 )] (20)

Here, the following equation (21) is obtained by substituting the equations (17) to (20) into the above equation (11).
Q = K _A × V _A + K _B × V _B + K _C × V _C + K _D × V _D
= K _A × f ₁ [M _I1 (T _H1 ), M _I2 (T _H2 ), M _I3 (T _H3 )]
+ K _B × f ₂ [M _I1 (T _H1 ), M _I2 (T _H2 ), M _I3 (T _H3 )]
+ K _C × f ₃ [M _I1 (T _H1 ), M _I2 (T _H2 ), M _I3 (T _H3 )]
+ K _D × f ₄ [M _I1 (T _H1 ), M _I2 (T _H2 ), M _I3 (T _H3 )] (21)

As shown in the above equation (21), the calorific value Q per unit volume of the mixed gas is the heat dissipation coefficient M _I1 (T _H1 ) of the mixed gas when the temperature of the heating element 61 is T _H1 , T _H2 , T _H3. ), M _I2 (T _H2 ), M _I3 (T _H3 ) are given as equations. Therefore, the calorific value Q of the mixed gas is given by the following equation (22), where g is a symbol representing a function.
Q = g [M _I1 (T _H1 ), M _I2 (T _H2 ), M _I3 (T _H3 )] (22)

Therefore, if the above equation (22) is obtained in advance for a mixed gas composed of gas A, gas B, gas C, and gas D, the volume ratio V _A of gas _A , the volume ratio V _B of gas _B , and the gas C The inventor has found that the calorific value Q per unit volume of the measurement target mixed gas whose volume ratio V _C and the volume ratio V _{D of the} gas D are unknown can be easily calculated. Specifically, the radiation coefficient M _I1 (T _H1) of the measurement target mixed gas when the heat producing temperature is T _H1, T _H2, T _H3 of the heater element _{61, M I2 (T H2)} , M I3 (T H3 ) Are respectively measured and substituted into the equation (22), the calorific value Q of the measurement target mixed gas can be uniquely obtained.

Further, thermal characteristics such as the amount of heat generated by the gas, the heat dissipation coefficient, and the thermal conductivity depend on the pressure of the gas. Therefore, by adding an independent variable of the pressure Ps of the measurement target mixed gas to the equation of the calorific value Q given by the above formula (22), as shown in the following formula (23), the calculation accuracy of the calorific value Q is increased. .
Q = g [M _I1 (T _H1 ), M _I2 (T _H2 ), M _I3 (T _H3 ), Ps] (23)

Additionally, the radiation coefficient M _I of the mixed gas, as shown in Equation (9), dependent on the resistance value R _H of the heater element 61, the resistance value R _I of the first temperature measurement element 62, the. Therefore, the present inventor has determined that the heating value Q per unit volume of the mixed gas is the heating element 61 when the temperature of the heating element 61 is T _H1 , T _H2 , T _H3 as shown in the following equation (24). The resistance values R _H1 (T _H1 ), R _H2 (T _H2 ), R _H3 (T _H3 ), the resistance value R _{I of} the first temperature measuring element 62 in contact with the mixed gas, and the pressure Ps are defined as variables. I found that the equation is also given.
Q = g [R _H1 (T _H1 ), R _H2 (T _H2 ), R _H3 (T _H3 ), R _I , Ps] (24)

Therefore, the resistance value R _H1 of the heater element 61 when the heat producing temperature of the heater element 61 that comes into contact with the measurement target mixed gas is _{T H1, T H2, T H3} (T H1), R H2 (T H2), R H3 ( T _H3 ), the resistance value R _{I of} the first temperature measuring element 62 in contact with the measurement target mixed gas, and the pressure Ps of the measurement target mixed gas are also measured and substituted into the equation (24). The calorific value Q of the target mixed gas can be uniquely determined.

Further, the calorific value Q per unit volume of the mixed gas is, as shown in the following formula (25), the energization current I _H1 of the heating element 61 when the temperature of the heating element 61 is T _H1 , T _H2 , T _H3. (T _H1 ), I _H2 (T _H2 ), I _H3 (T _H3 ), an energization current I _{I of} the first temperature measuring element 62 in contact with the mixed gas, and a pressure Ps of the mixed gas as variables But given.
Q = g [I _H1 (T _H1 ), I _H2 (T _H2 ), I _H3 (T _H3 ), I _I , Ps] (25)

Alternatively, the calorific value Q per unit volume of the mixed gas is the voltage V _H1 (V _H1 applied to the heating element 61 when the temperature of the heating element 61 is T _H1 , T _H2 , T _H3 , as shown in the following equation (26). T _H1 ), V _H2 (T _H2 ), V _H3 (T _H3 ), the voltage V _I applied to the first temperature measuring element 62 in contact with the mixed gas, and the pressure Ps of the mixed gas are variables. Given.
Q = g [V _H1 (T _H1 ), V _H2 (T _H2 ), V _H3 (T _H3 ), V _I , Ps] (26)

Alternatively, the calorific value Q per unit volume of the mixed gas is an analog connected to the heating element 61 when the temperature of the heating element 61 is T _H1 , T _H2 , T _H3 as shown in the following equation (27). _-The output signal AD _H1 (T _H1 ), AD _H2 (T _H2 ), AD _H3 (T _H3 ) of the digital conversion circuit (hereinafter referred to as “A / D conversion circuit”) and the first measurement in contact with the mixed gas. An equation having the output signal AD _{I of} the A / D conversion circuit connected to the temperature element 62 and the pressure Ps of the mixed gas as variables is also given. For example, when the A / D conversion circuit is a double integration type, the output signal of the A / D conversion circuit is a count value.
Q = g [AD _H1 (T _H1 ), AD _H2 (T _H2 ), AD _H3 (T _H3 ), AD _I , Ps] (27)

Therefore, the calorific value Q per unit volume of the mixed gas is an electric signal from the heating element 61 when the heating temperature of the heating element 61 is T _H1 , T _H2 , T _H3 as shown in the following equation (28). _{S H1 (T H1), S} H2 (T H2), S H3 and (T _H3), and the electric signal S _I from the first temperature measurement element 62 that comes into contact with the mixed gas, the pressure Ps of the mixed gas, the variable Is given by the equation
Q = g [S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3 ), S _I , Ps] (28)

The pressure Ps of the mixed gas is measured using a pressure sensor. The pressure sensor includes a strain gauge made of, for example, an electric resistance element. The strain gauge is deformed by pressure and its electric resistance changes. Therefore, the output voltage of the pressure sensor or the output signal of the A / D conversion circuit connected to the pressure sensor correlates with the pressure Ps of the mixed gas. Therefore, the calorific value Q per unit volume of the mixed gas is expressed by the electric signals S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3 ) from the heating element 61 as shown in the following equation (29). ), And the electric signal S _I from the first temperature measuring element 62 and the electric signal S _P from the pressure sensor are also given as equations.
Q = g [S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3 ), S _I , S _P ] (29)

Further, if the temperature of the mixed gas is constant, the electric signal S _I from the first temperature measuring element 62 is a constant. In this case, the calorific value Q per unit volume of the mixed gas is expressed by the electric signals S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T and _H3), and an electric signal S _P output from the pressure sensor, also given in the equation as a variable.
Q = g [S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3 ), S _P ] (30)

The gas components of the mixed gas are not limited to four types. For example, when the mixed gas is composed of n types of gas components, first, at least n−1 types of heat generation temperatures T _H1 , T _H2 , T _{H3 ,.} , S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3 ),..., S _Hn-1 (T _Hn-1 ), and the _first contact with the mixed gas. and an electric signal S _I from the first temperature measurement element 62, obtained in advance an equation for the electric signal S _P output from the pressure sensor, a variable. Then, the heating element 61 in contact with the measurement target mixed gas whose volume ratio of each of the n types of gas components at the n-1 types of heat generation temperatures T _H1 , T _H2 , T _{H3 ,.} electrical signals from the S _H1 (T _H1), the contact with the _{S H2 (T H2), S} H3 (T H3), ···, and the value of _{S Hn-1 (T Hn-} 1), the measurement target mixed gas A measurement object is obtained by measuring the value of the electric signal S _I from one temperature measuring element 62 and the value of the electric signal S _P from the pressure sensor in contact with the measurement target mixed gas and substituting them into the equation (31). The calorific value Q per unit volume of the mixed gas can be uniquely determined.
Q = g [S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3 ), ..., S _Hn-1 (T _Hn-1 ), S _I , S _P ] ... 31)

However, the mixed gas, methane (CH ₄₎ as a gas component in addition to the propane (C ₃ H _8), a j is a natural number, methane (CH ₄₎ and other than propane (C ₃ H ₈₎ alkane (C _j H _{2j + 2} ), alkanes other than methane (CH ₄ ) and propane (C ₃ H ₈ ) (C _j H _{2j + 2} ), and mixtures of methane (CH ₄ ) and propane (C ₃ H ₈ ) This does not affect the calculation of equation (31). For example, ethane (C ₂ H ₆ ), butane (C ₄ H ₁₀ ), pentane (C ₅ H ₁₂ ), and hexane (C ₆ H ₁₄ ) are respectively represented by the following formulas (32) to (35): The equation (31) may be calculated by regarding the mixture as methane (CH ₄ ) and propane (C ₃ H ₈ ) multiplied by a predetermined coefficient.
C ₂ H ₆ = 0.5 CH ₄ + 0.5 C ₃ H ₈ ... (32)
C ₄ H ₁₀ = -0.5 CH ₄ + 1.5 C ₃ H ₈ ... (33)
C ₅ H ₁₂ = -1.0 CH ₄ + 2.0 C ₃ H ₈ ... (34)
C ₆ H ₁₄ = -1.5 CH ₄ + 2.5 C ₃ H ₈ ... (35)

Accordingly, the z as a natural number, a mixed gas consisting of n kinds of gas components methane (CH ₄₎ as a gas component in addition to the propane (C ₃ H _8), methane (CH ₄₎ and propane (C ₃ H ₈ ) if it contains non-z type of alkane _{(C j H 2j + 2)} , and an electric signal S _H from the at least n-z-1 kinds heating element 61 in the heat generating temperature of the first temperature measurement element 62 and an electric signal S _I, and the electric signal S _P output from the pressure sensor, may be an equation for the variable.

When the type of the gas component of the mixed gas used in the calculation of the equation (31) is the same as the type of the gas component of the measurement target mixed gas whose calorific value Q per unit volume is unknown, the measurement target mixed gas Of course, equation (31) can be used to calculate the calorific value Q. Further, when the measurement target mixed gas is composed of less than n kinds of gas components and less than n kinds of gas components are included in the mixed gas used in the calculation of the equation (31), Equation (31) can be used. For example, the mixed gas used in the calculation of the equation (31) has four types of gas components, methane (CH ₄ ), propane (C ₃ H ₈ ), nitrogen (N ₂ ), and carbon dioxide (CO ₂ ). When it is included, the measurement target mixed gas may not contain nitrogen (N ₂ ) but may contain only three types of gas components, methane (CH ₄ ), propane (C ₃ H ₈ ), and carbon dioxide (CO ₂ ). The equation (31) can be used to calculate the calorific value Q of the measurement target mixed gas.

Furthermore, when the mixed gas used for the calculation of the equation (31) includes methane (CH ₄ ) and propane (C ₃ H ₈ ) as gas components, the measurement target mixed gas is used for the calculation of the equation (31). Even if alkane (C _j H _{2j + 2} ) not contained in the mixed gas is included, the formula (31) can be used. This is because, as described above, the methane (CH ₄₎ and propane (C ₃ H ₈₎ other than the alkane _{(C j H 2j + 2)} , a mixture of methane (CH ₄₎ and propane (C ₃ H ₈₎ Even if it considers, it is because it does not affect calculation of the emitted-heat amount Q per unit volume using (31) Formula.

Here, the calorific value measurement system 20 according to the embodiment shown in FIG. 6 is disposed in the chamber 101, which is a container into which each of a plurality of sample mixed gases is injected, and the first shown in FIG. A microchip 8 including a temperature measuring element 62 and a heating element 61 to which a plurality of different voltages are applied to generate heat at a plurality of heating temperatures _TH , and a pressure sensor 201 are provided. Furthermore, the calorific value measurement system 20 shown in FIG. 6 includes the value of the electric signal S _I from the first temperature measuring element 62 depending on the temperature T _I of each of the plurality of sample mixed gases and the plurality of exothermic temperatures T _H. comprising the value of the electric signal S _H from the heater element 61 in each, and the value of the electric signal S _P output from the pressure sensor 201, a measurement unit 301 that measures. Furthermore, the calorific value measurement system 20, the known calorific value Q of a plurality of mixed gases, the value of the electric signal S _I from the first temperature measuring element 62, electricity from the heat generating element 61 at a plurality of heat producing temperatures the value of the signal S _H, and based on the value of the electric signal S _P output from the pressure sensor 201, the electric signal S _I from the first temperature measuring element 62, the electric signal from the heater element 61 at a plurality of heat producing temperatures T _H S _H, and the electric signal S _P as independent variables from the pressure sensor 201, and a formula creation module 302 that creates a calorific value calculation formula to the calorific value Q of the gas as the dependent variable. The sample mixed gas includes a plurality of types of gas components.

  The microchip 8 is disposed in the chamber 101 via the heat insulating member 18. The chamber 101 is connected to a flow path 102 for sending the sample mixed gas to the chamber 101 and a flow path 103 for discharging the sample mixed gas from the chamber 101 to the outside.

  For example, a gauge pressure sensor or an absolute pressure sensor can be used as the pressure sensor 201 for measuring the pressure of the gas inside the chamber 101. The pressure sensor 201 includes a pressure sensitive element. As the pressure sensitive element, for example, a semiconductor diaphragm type, a capacitance type, an elastic diaphragm type, a piezoelectric type, and a vibration type can be used.

  When four types of sample mixed gases with different calorific values Q are used, as shown in FIG. 7, the first gas cylinder 50A for storing the first sample mixed gas and the second sample mixed gas for storing the second sample mixed gas are used. A second gas cylinder 50B, a third gas cylinder 50C for storing a third sample mixed gas, and a fourth gas cylinder 50D for storing a fourth sample mixed gas are prepared. The first gas cylinder 50A is connected to a first gas pressure adjuster 31A that adjusts the pressure of the first sample mixed gas via a flow path 91A. In addition, a first flow rate control device 32A is connected to the first gas pressure regulator 31A via a flow path 92A. The first flow rate control device 32A controls the flow rate of the first sample mixed gas sent to the calorific value measurement system 20 via the flow path 92A and the flow path 102.

  A second gas pressure regulator 31B is connected to the second gas cylinder 50B via a flow path 91B. The second flow rate controller 32B is connected to the second gas pressure regulator 31B via a flow path 92B. The second flow rate control device 32B controls the flow rate of the second sample mixed gas sent to the calorific value measurement system 20 via the flow paths 92B, 93, 102.

  A third gas pressure regulator 31C is connected to the third gas cylinder 50C via a flow path 91C. In addition, a third flow rate control device 32C is connected to the third gas pressure regulator 31C via a flow path 92C. The third flow rate control device 32C controls the flow rate of the third sample mixed gas sent to the calorific value measurement system 20 via the flow paths 92C, 93, 102.

  A fourth gas pressure regulator 31D is connected to the fourth gas cylinder 50D via a flow path 91D. In addition, a fourth flow rate control device 32D is connected to the fourth gas pressure regulator 31D via a flow path 92D. The fourth flow rate control device 32D controls the flow rate of the fourth sample mixed gas sent to the calorific value measurement system 20 via the flow paths 92D, 93, and 102.

Each of the first to fourth sample mixed gases is natural gas having a different calorific value, for example. Each of the first to fourth sample mixed gases has different ratios of four types of gas components, for example, methane (CH ₄ ), propane (C ₃ H ₈ ), nitrogen (N ₂ ), and carbon dioxide (CO ₂ ). Including.

After the first sample mixed gas is filled into the chamber 101 shown in FIG. 6, the pressure sensor 201 outputs an electric signal S _P which depends on the pressure of the first sample mixed gas. The first temperature measurement element 62 of the microchip 8 shown in Figs. 1 and 2 outputs an electric signal S _I that is dependent on the temperature of the first sample mixed gas. The heating element 61 is supplied with drive powers P _H1 , P _H2 , and P _H3 from the drive circuit 303 shown in FIG. When the driving powers P _H1 , P _H2 , and P _H3 are given, the heating element 61 in contact with the first sample mixed gas has, for example, a temperature T _{H1 of} 100 ° C., a temperature T _{H2 of} 150 ° C., and a temperature of 200 ° C. generates heat T _H3, the electric signal S _H1 at the heat producing temperature T _H1 (T _H1), an electric signal S _H2 at the heat producing temperature T _H2 (T _H2), and the electric signal S _H3 at the heat producing temperature T _H3 of (T _H3) output To do.

After the first sample mixed gas is removed from the chamber 101, the chamber 101 is sequentially filled with the second to fourth sample mixed gases. After the second sample mixed gas is filled into the chamber 101, the pressure sensor 201 outputs an electric signal S _P which depends on the pressure of the second sample mixed gas. The first temperature measurement element 62 of the microchip 8 shown in Figs. 1 and 2 outputs an electric signal S _I that is dependent on the temperature of the second sample mixed gas. Heater element 61 coming into contact with the second sample mixed gas, electric signals S _H1 (T _H1) at the heat producing temperature T _H1, the electric signal S at the heat producing temperature T _{H2 H2} (T _H2), and the electric signal S at the heat producing temperature T _{H3 H3} (T _H3 ) is output.

After the third sample mixed gas is filled into the chamber 101 shown in FIG. 6, the pressure sensor 201 outputs an electric signal S _P which depends on the pressure of the third sample mixed gas. The first temperature measurement element 62 of the microchip 8 shown in Figs. 1 and 2 outputs an electric signal S _I that is dependent on the temperature of the third sample mixed gas. The third sample mixed in contact with the gas heating element 61, the electric signal S _H1 at the heat producing temperature T _H1 (T _H1), an electric signal S _H2 at the heat producing temperature T _H2 (T _H2), and an electrical signal at a heat producing temperature T _H3 S _H3 (T _H3 ) is output.

After the fourth sample mixed gas is filled into the chamber 101 shown in FIG. 6, the pressure sensor 201 outputs an electric signal S _P which depends on the pressure of the fourth sample mixed gas. The first temperature measurement element 62 of the microchip 8 shown in Figs. 1 and 2 outputs an electric signal S _I that is dependent on the temperature of the fourth sample mixed gas. The heater element 61 in contact with the fourth sample mixed gas, an electrical signal at a heat producing temperature T the electrical signal at _H1 S _H1 (T _H1), electric signals S _H2 at the heat producing temperature T _H2 (T _H2), and heating temperature T _H3 S _H3 (T _H3 ) is output.

In addition, when each sample mixed gas contains n types of gas components, the heating element 61 shown in FIGS. 1 and 2 of the microchip 8 is caused to generate heat at at least n−1 types of different temperatures. However, as described above, methane (CH ₄₎ and propane (C ₃ H ₈₎ other than the alkane (C _j H _{2j + 2)} is regarded as a mixture of methane (CH ₄₎ and propane (C ₃ H ₈₎ sell. Accordingly, a sample mixed gas composed of n kinds of gas components, where z is a natural number, is added to methane (CH ₄ ) and propane (C ₃ H ₈ ) as gas components, and z kinds of alkanes (C _j H _{2j + 2} ). In the case where the heat generating element 61 is included, the heat generating element 61 generates heat at at least nz-1 different temperatures.

As shown in FIG. 6, the microchip 8 and the pressure sensor 201 are connected to a central processing unit (CPU) 300 including a measurement unit 301 via an A / D conversion circuit 304. An electrical signal storage device 401 is connected to the CPU 300. Measuring unit 301, the value of the electric signal S _I from the first temperature measuring element 62, the electric signal S _H1 (T _H1) at a heat producing temperature T _H1, from the heater element 61, the electric signal S _H2 at the heat producing temperature T _H2, (T _H2) storage, and a value of the electric signal S _H3 (T _H3) at the heat producing temperature T _H3, and the value of the electric signal S _P output from the pressure sensor 201, is measured, the measured value into an electric signal storage device 401 To do.

The electric signal S _I from the first temperature measuring element 62 includes a resistance value R _I of the first temperature measuring element 62, an energization current I _I of the first temperature measuring element 62, and the first temperature measuring element 62. Either the voltage V _I or the output signal AD _I of the A / D conversion circuit 304 connected to the first temperature measuring element 62 may be used. Similarly, the electric signal S _H from the heater element 61, the resistance value R _H of the heater element 61, current I _H flowing in the heater element 61, which is connected to the voltage V _H and the heating element 61, according to the heating elements 61 it may be any of the output signal AD _H of the a / D conversion circuit 304. Also, the electric signal S _P output from the pressure sensor 201, for example, the resistance value of the strain gauge pressure sensor 201 is provided, current applied strain gauges, voltage across the strain gauge, and the strain gauges connected to the A / D converter circuit Any of the output signals 304 may be used.

Equation generating unit 302 included in the CPU300, for example each of the known calorific value Q of the first to fourth sample mixed gases, the plurality of measured values for the electric signals S _I from the first temperature measurement element 62 And a plurality of measured values of the electric signals S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3 ) from the heating element 61 and a plurality of measured values of the electric signal S _P from the pressure sensor 201. And collect. Wherein creating unit 302 further collected calorific value Q, the electric signals S _I, electric signals S _H, and based on the value of the electric signal S _P, by multivariate analysis, the electrical signal from the first temperature measuring element 62 S _I , electric signals S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3 ) from the heating element 61 and the electric signal S _P from the pressure sensor 201 are independent variables, and the amount of heat generated by the gas A calorific value calculation formula with Q as a dependent variable is calculated.

  Note that “multivariate analysis” refers to A.I. J Smol and B.M. In support vector regression, multiple regression analysis, and Japanese Patent Application Laid-Open No. 5-141999, disclosed in Scholkopf's “A Tutor on Support Vector Regression” (NeuroCOLt Technical Report (NC-TR-98-030), 1998). Includes the disclosed fuzzy quantification theory class II.

  The calorific value measurement system 20 further includes a formula storage device 402 connected to the CPU 300. The formula storage device 402 stores the calorific value calculation formula created by the formula creation unit 302. Further, an input device 312 and an output device 313 are connected to the CPU 300. As the input device 312, for example, a keyboard and a pointing device such as a mouse can be used. As the output device 313, an image display device such as a liquid crystal display and a monitor, a printer, and the like can be used.

  Next, a method for creating a calorific value calculation formula using the calorific value measurement system 20 according to the embodiment will be described using the flowchart shown in FIG. In the following, an example of creating a calorific value calculation formula by setting the atmospheric pressure inside the chamber 101 to atmospheric pressure, 5 kPa, 20 kPa, and 30 kPa will be described.

(A) In step S100, with the valves of the second to fourth flow control devices 32B-32D shown in FIG. 7 closed, the valve of the first flow control device 32A is opened and placed in the chamber 101 shown in FIG. A first sample mixed gas is introduced. In step S101, the pressure inside the chamber 101 is made equal to the atmospheric pressure. Measuring portion 301 measures the value of the electric signal S _P representing the pressure from the pressure sensor 201, the electric signal storage device 401. The measurement unit 301, the value of the electric signal S _I from the first temperature measurement element 62 in contact with the first sample mixed gas is measured, in the electric signal storage device 401. Next, the drive circuit 303 _applies drive power P _H1 to the heating element 61 shown in FIGS. 1 and 2, and causes the heating element 61 to generate heat at 100 ° C. The measurement unit 301 shown in FIG. 6 stores the value of the electrical signal S _H1 (T _H1 ) from the heating element 61 that generates heat at 100 ° C. in the electrical signal storage device 401.

(B) In step S102, the measurement unit 301 determines whether switching of the atmospheric pressure inside the chamber 101 is completed. If switching to 5 kPa, 20 kPa, and 30 kPa has not been completed, the process returns to step S101, and the pressure inside the chamber 101 is set to 5 kPa. Further, the measuring unit 301 has a value of the electric signal S _P from the pressure sensor 201 at 5 kPa, a value of the electric signal S _I from the first temperature measuring element 62, and an electric signal S from the heating element 61 that generates heat at 100 ° C. The value of _H1 (T _H1 ) is stored in the electrical signal storage device 401.

(C) In step S102 again, the measurement unit 301 determines whether or not the switching of the atmospheric pressure inside the chamber 101 has been completed. If switching to 20 kPa and 30 kPa has not been completed, the process returns to step S101, and the pressure inside the chamber 101 is set to 20 kPa. Further, the measurement unit 301 has a value of the electric signal S _P from the pressure sensor 201 at 20 kPa, a value of the electric signal S _I from the first temperature measuring element 62, and an electric signal S from the heating element 61 that generates heat at 100 ° C. The value of _H1 (T _H1 ) is stored in the electrical signal storage device 401.

(D) In step S102 again, the measurement unit 301 determines whether or not the switching of the atmospheric pressure inside the chamber 101 has been completed. If the switching to 30 kPa has not been completed, the process returns to step S101, and the pressure inside the chamber 101 is set to 30 kPa. Further, the measuring unit 301 has a value of the electric signal S _P from the pressure sensor 201 at 30 kPa, a value of the electric signal S _I from the first temperature measuring element 62, and an electric signal S from the heating element 61 that generates heat at 100 ° C. The value of _H1 (T _H1 ) is stored in the electrical signal storage device 401.

(E) When switching of the atmospheric pressure in the chamber 101 is completed, the process proceeds from step S102 to step S103. In step S103, the drive circuit 303 determines whether or not the temperature switching of the heating element 61 shown in FIGS. 1 and 2 has been completed. If the switching to the temperature of 150 ° C. and the temperature of 200 ° C. has not been completed, the process returns to step S101, and the drive circuit 303 shown in FIG. 6 causes the heating element 61 shown in FIGS. Thereafter, repeated loop of steps S101 and S102, the measurement unit 301 shown in FIG. 6, the atmospheric pressure, 5 kPa, 20 kPa, and the value of the electric signal S _P output from the pressure sensor 201 at each of 30 kPa, the first temperature measurement The value of the electric signal S _I from the element 62 and the value of the electric signal S _H1 (T _H1 ) from the heating element 61 that generates heat at 150 ° C. are stored in the electric signal storage device 401.

(F) In step S103 again, the drive circuit 303 determines whether or not the temperature switching of the heating element 61 shown in FIGS. 1 and 2 has been completed. If the switching to the temperature of 200 ° C. has not been completed, the process returns to step S101, and the drive circuit 303 shown in FIG. 6 causes the heating element 61 shown in FIGS. 1 and 2 to generate heat at 200 ° C. Thereafter, repeated loop of steps S101 and S102, the measurement unit 301 shown in FIG. 6, the atmospheric pressure, 5 kPa, 20 kPa, and the value of the electric signal S _P output from the pressure sensor 201 at each of 30 kPa, the first temperature measurement The value of the electric signal S _I from the element 62 and the value of the electric signal S _H1 (T _H1 ) from the heating element 61 that generates heat at 200 ° C. are stored in the electric signal storage device 401.

  (G) When the switching of the temperature of the heating element 61 is completed, the process proceeds from step S103 to step S104. In step S104, it is determined whether switching of the sample mixed gas is completed. If switching to the second to fourth sample mixed gases has not been completed, the process returns to step S100. In step S100, the first flow control device 32A shown in FIG. 7 is closed, the valves of the second flow control device 32B are opened while the valves of the third to fourth flow control devices 32C-32D are closed, and FIG. A second sample mixed gas is introduced into the chamber 101 shown in FIG.

(H) Similar to the first sample mixed gas, the loop of step S101 to step S103 is repeated. The measurement unit 301 has the electric signal S _P from the pressure sensor 201 in contact with the second sample mixed gas at atmospheric pressure, 5 kPa, 20 kPa, and 30 kPa, and the electric signal S _I from the first temperature measuring element 62. And the values of the electric signals S _H1 (T _H1 ), S _H2 (T _H2 ), and S _H3 (T _H3 ) from the heating element 61 that generates heat at 100 ° C., 150 ° C., and 200 ° C. are stored in the electric signal storage. Save to device 401.

(I) Thereafter, the loop from step S100 to step S104 is repeated. Thus, the measurement unit 301, atmospheric pressure, 5 kPa, 20 kPa, and the value of the electric signal S _P output from the pressure sensor 201 in contact with the third sample mixed gas at the respective 30 kPa, electricity from the first temperature measuring element 62 The value of the signal S _{I and} the values of the electrical signals S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3 ) from the heating element 61 that generates heat at 100 ° C., 150 ° C., and 200 ° C. Save in the electrical signal storage device 401. The measurement unit 301, atmospheric pressure, 5 kPa, 20 kPa, and the value of the electric signal S _P output from the pressure sensor 201 in contact with the fourth sample mixed gases at each of 30 kPa, the electrical signal from the first temperature measuring element 62 The values of S _I and the electric signals S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3 ) from the heating elements 61 that generate heat at 100 ° C., 150 ° C. and 200 ° C. Save in the signal storage device 401.

(J) In step S105, the input device 312 sends the formula generating unit 302 the known calorific value Q of the first sample mixed gas, the known calorific value Q of the second sample mixed gas, A known calorific value Q of the sample mixed gas and a known calorific value Q of the fourth sample mixed gas are input. Further, the formula creating unit 302 receives a plurality of measured values of the electric signal S _I from the first temperature measuring element 62 and the electric signals S _H1 (T _H1 ), S from the heating element 61 from the electric signal storage device 401. A plurality of measured values of _H2 (T _H2 ) and S _H3 (T _H3 ) and a plurality of measured values of the electric signal S _P from the pressure sensor 201 are read out.

(K) in step S106, the calorific values Q of the first to fourth sample mixed gases, a plurality of measurements of the electric signal S _I from the first temperature measuring element 62, electricity from the heating elements 61 signal _{S H1 (T H1), S} H2 (T H2), and a plurality of measured values of S _H3 (T _H3), and a plurality of measured values of the electric signal S _P output from the pressure sensor 201, based on the formula creation Unit 302 performs multiple regression analysis. Based on the multiple regression analysis, the formula creating unit 302 generates an electrical signal S _I from the first temperature measuring element 62, and electrical signals S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3), and an electric signal S _P as independent variables from the pressure sensor 201, calculates the heat generation amount Q of the gas calorific value calculation formula to dependent variable. Thereafter, in step S107, the formula creation unit 302 stores the created calorific value calculation formula in the formula storage device 402, and the calorific value calculation formula creation method according to the embodiment ends.

  As described above, according to the method for creating a calorific value calculation formula using the calorific value measurement system 20 according to the embodiment, a calorific value calculation formula capable of uniquely calculating the value of the calorific value Q of the gas is created. It becomes possible.

Next, the function of the calorific value measurement system 20 according to the embodiment when measuring the calorific value Q of the measurement target mixed gas whose calorific value Q is unknown will be described. For example, natural gas with unknown calorific value Q including methane (CH ₄ ), propane (C ₃ H ₈ ), nitrogen (N ₂ ), carbon dioxide (CO ₂ ), etc. with unknown volume ratio. A gas is introduced into the chamber 101. The pressure sensor 201 outputs an electric signal S _P which depends on the pressure of the measurement target mixed gas. The first temperature measuring element of the microchip 8 shown in Figs. 1 and 2 62 outputs an electric signal S _I that is dependent on the temperature of the measurement target mixed gas. Next, the heating element 61 is supplied with driving powers P _H1 , P _H2 , and P _H3 from the driving circuit 303 shown in FIG. When the driving powers P _H1 , P _H2 , and P _H3 are given, the heating element 61 in contact with the measurement target mixed gas has, for example, a temperature T _{H1 of} 100 ° C., a temperature T _{H2 of} 150 ° C., and a temperature T _{H3 of} 200 ° C. _Generates an electric signal S _H1 (T _H1 ) at a heat generation temperature T _H1, an electric signal S _H2 (T _H2 ) at a heat generation temperature T _H2 , and an electric signal S _H3 (T _H3 ) at a heat generation temperature T _H3 .

Measuring unit 301 shown in FIG. 6 is in contact with the measurement target mixed gas, the value of the electric signal S _P output from the pressure sensor 201 which depends on the pressure of the measurement target mixed gas, the temperature dependent T _I of the measurement target mixed gas the value of the electric signal S _I from the first temperature measurement element 62, the electric signal S _H1 at a heat producing temperature T _H1, from the heater element 61 in contact with the measurement target mixed gas (T _H1), an electric signal S _H2 at the heat producing temperature T _H2, (T _H2 ) and the value of the electric signal S _H3 (T _H3 ) at the heat generation temperature T _H3 are measured, and the measured value is stored in the electric signal storage device 401.

As described above, the expression storage device 402 includes the electric signal S _I from the first temperature measuring element 62, the electric signal S _H1 (T _H1 ) from the heating element 61 having a heat generation temperature T _H1 of 100 ° C., and heat generation. An electrical signal S _H2 (T _H2 ) from the heating element 61 having a temperature T _H2 of 150 ° C., an electrical signal S _H3 (T _H3 ) from the heating element 61 having a heating temperature T _H3 of 200 ° C., and a pressure sensor 201 and an electric signal S _P, and the independent variable, and saves the calorific value calculation formula to the calorific value Q of the gas as the dependent variable.

The calorific value measurement system 20 according to the embodiment further includes a calorific value calculation unit 305. Calorific value calculation section 305, the independent variable of the electric signal S _I from the first temperature measuring element 62 of the calorific value calculation formula, the independent variable of the electric signals S _H from the heater element 61, and the electrical signal from the pressure sensor 201 the independent variables of S _P, measured values of the electric signals S _I from the first temperature measuring element 62, measured values of the electric signals S _H from the heater element 61, and the measured value of the electric signal S _P output from the pressure sensor 201 Are respectively substituted, and the measured value of the calorific value Q of the measurement target mixed gas injected into the chamber 101 is calculated. A heat generation amount storage device 403 is further connected to the CPU 300. The calorific value storage device 403 stores the calorific value Q of the measurement target mixed gas calculated by the calorific value calculation unit 305.

  Next, a calorific value measurement method using the calorific value measurement system 20 according to the embodiment will be described using the flowchart shown in FIG.

(A) In step S200, the measurement target mixed gas is introduced into the chamber 101 shown in FIG. In step S201, the measurement unit 301, the value of the electric signal S _P output from the pressure sensor 201 in contact with the measurement target mixed gas was measured, in the electric signal storage device 401. The measurement unit 301, the value of the electric signal S _I from the first temperature measurement element 62 coming into contact with the measurement target mixed gas was measured, in the electric signal storage device 401. Next, the drive circuit 303 _applies drive power P _H1 to the heating element 61 shown in FIGS. 1 and 2, and causes the heating element 61 to generate heat at 100 ° C. The measurement unit 301 shown in FIG. 6 stores the value of the electrical signal S _H1 (T _H1 ) from the heating element 61 that contacts the measurement target mixed gas and generates heat at 100 ° C. in the electrical signal storage device 401.

(B) In step S202, the drive circuit 303 shown in FIG. 6 determines whether or not the temperature switching of the heating element 61 shown in FIGS. 1 and 2 is completed. When the switching to the temperature of 150 ° C. and the temperature of 200 ° C. is not completed, the process returns to step S201, and the drive circuit 303 gives the drive power P _H2 to the heating element 61 shown in FIGS. Is heated at 150 ° C. 6 stores the value of the electric signal S _H2 (T _H2 ) from the heating element 61 that contacts the measurement target mixed gas and generates heat at 150 ° C. in the electric signal storage device 401.

(C) In step S202 again, it is determined whether or not the temperature switching of the heating element 61 shown in FIGS. 1 and 2 has been completed. If the switching to the temperature of 200 ° C. has not been completed, the process returns to step S201, and the drive circuit 303 _{supplies the} drive power P _H3 to the heating element 61 shown in FIGS. Causes fever. 6 stores the value of the electrical signal S _H3 (T _H3 ) from the heating element 61 that is in contact with the measurement target mixed gas and generates heat at 200 ° C. in the electrical signal storage device 401.

(D) When the switching of the temperature of the heating element 61 is completed, the process proceeds from step S202 to step S203. In step S203, the calorific value calculation unit 305 shown in FIG. 6 receives the electric signal S _I from the first temperature measuring element 62 and the electric signals S _H1 (T _H1 ), S from the first temperature measuring element 62 from the formula storage device 402. _{H2 (T H2), S H3} (T H3), and an electric signal S _P as independent variables from the pressure sensor 201 reads the calorific value Q of the gas calorific value calculation formula to dependent variable. Further, the calorific value calculation unit 305 from the electric signal storage device 401, the measured value of the electric signal S _I from the first temperature measurement element 62 coming into contact with the measurement target mixed gas, from the heater element 61 in contact with the measurement target mixed gas electric signals _{S H1 (T H1), S} H2 (T H2), the measured value of S _H3 (T _H3), and reads out the measured value of the electric signal S _P output from the pressure sensor 201.

(E) In step S204, the calorific value calculation unit 305 determines the electric signal S _I , electric signals S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3 ), and electric signals of the calorific value calculation formula. the independent variables of S _P, respectively substitutes the measured values to calculate the calorific value Q of the measurement target mixed gas. Thereafter, the calorific value calculation unit 305 stores the calculated calorific value Q in the calorific value storage device 403, and ends the calorific value measurement method according to the embodiment.

According to the calorific value calculation method according to the embodiment described above, the value of the electric signal S _I from the first temperature measuring element 62 in contact with the measurement target mixed gas without using an expensive gas chromatography device or sound velocity sensor. And the values of the electrical signals S _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3 ) from the heating element 61 in contact with the measurement target mixed gas, and the pressure sensor 201 in contact with the measurement target mixed gas. the value of the electric signal S _P from, it is possible to measure the calorific value Q of the measurement target mixed gas.

Natural gas has a different component ratio of hydrocarbons depending on the gas field. Natural gas includes nitrogen (N ₂ ), carbon dioxide (CO ₂ ) and the like in addition to hydrocarbons. Therefore, the volume ratio of the gas component contained in the natural gas differs depending on the output gas field, and the calorific value Q of the natural gas is often unknown even if the type of the gas component is known. Moreover, even if it is the natural gas derived from the same gas field, the calorific value Q is not always constant and may change depending on the sampling time.

  Conventionally, when collecting the usage fee of natural gas, a method has been used in which charging is performed according to the volume used, not the calorific value Q used for natural gas. However, since the calorific value Q differs depending on the production gas field from which natural gas is derived, it is not fair to charge the volume used. On the other hand, if the calorific value calculation method according to the embodiment is used, the type of gas component is known, but since the volume ratio of the gas component is unknown, the mixed gas such as natural gas whose calorific value Q is unknown The calorific value Q can be easily calculated. Therefore, it becomes possible to collect a fair usage fee.

  Also, it is not environmentally preferable to throw away the gas whose calorific value Q is measured into the atmosphere. Therefore, when measuring the calorific value Q of the gas in the gas pipe, a calorific value measuring system is provided in the gas pipe itself, or a calorific value measuring system is provided in the bypass path of the gas pipe, and the gas whose calorific value Q is measured is supplied to the gas pipe. It is preferable to return. At this time, the atmospheric pressure in the gas pipe can vary significantly. On the other hand, according to the calorific value measurement system according to the embodiment, by including an independent variable of pressure in the calorific value calculation formula, it becomes possible to suppress a calorific value calculation error due to fluctuations in atmospheric pressure. When the same pressure sensor 201 is used for both the creation of the calorific value calculation formula and the measurement of the calorific value, the pressure sensor 201 does not need to have a correction circuit. This is because it is possible to suppress the calculation error of the calorific value due to the fluctuation of the atmospheric pressure, as long as the electrical signal output from the pressure sensor 201 can be measured according to the pressure without necessarily measuring the accurate pressure value. is there.

(Other embodiments)
As described above, the present invention has been described according to the embodiment. However, it should not be understood that the description and drawings constituting a part of this disclosure limit the present invention. From this disclosure, various alternative embodiments, embodiments, and operation techniques should be apparent to those skilled in the art. For example, in the embodiment, the expression storage device 402 shown in FIG. 6 includes the electric signal from the pressure sensor 201, the electric signal from the first temperature measuring element 62 shown in FIG. 1, and the heating element 61 at a plurality of heating temperatures. An example in which the calorific value calculation formula with the calorific value Q as a dependent variable is stored as an independent variable has been described.

On the other hand, as described in the above equation (23), the calorific value Q of the gas is the heat radiation of the gas when the gas pressure Ps and the temperature of the heating element 61 are T _H1 , T _H2 , and T _H3 , respectively. An equation having variables M _I1 (T _H1 ), M _I2 (T _H2 ), and M _I3 (T _H3 ) as variables is also given. Therefore, the formula storage device 402 shown in FIG. 6 stores a calorific value calculation formula having the gas pressure and the heat dissipation coefficient of the gas at a plurality of heat generation temperatures of the heating element 61 as independent variables and the calorific value Q as a dependent variable. May be. In this case, the measurement unit 301 measures the measured value of the heat dissipation coefficient of the gas injected into the chamber 101 by causing the heating element 61 to generate heat at a plurality of heating temperatures. The fact that the heat dissipation coefficient of gas can be measured using the microchip 8 is as described in the above equation (9). The calorific value calculation unit 305 substitutes the measured value of the gas pressure and the measured value of the gas heat dissipation coefficient into the independent variable of the calorific value calculation formula stored in the formula storage device 402, and generates the calorific value Q of the gas. Calculate the measured value.

  Next, FIG. 10 shows the relationship between the heat dissipation coefficient of the mixed gas and the thermal conductivity when currents of 2 mA, 2.5 mA, and 3 mA are passed through the heating resistor. As shown in FIG. 10, the heat dissipation coefficient of the mixed gas and the thermal conductivity are generally in a proportional relationship. Therefore, the formula storage device 402 shown in FIG. 6 stores a calorific value calculation formula having the gas pressure and the thermal conductivity of the gas at a plurality of heating temperatures of the heating element 61 as independent variables and the calorific value Q as a dependent variable. May be. In this case, the measurement unit 301 measures the measured value of the thermal conductivity of the gas injected into the chamber 101 by causing the heating element 61 to generate heat at a plurality of heating temperatures. The calorific value calculation unit 305 substitutes the measured value of the gas pressure and the measured value of the thermal conductivity of the gas into the independent variables of the calorific value calculation formula stored in the formula storage device 402, and generates the calorific value of the gas. Calculate the measured value of Q.

  Thus, it should be understood that the present invention includes various embodiments and the like not described herein.

(Example)
First, 12 types of sample mixed gases with known calorific values Q were prepared. Each of the twelve sample mixed gases includes methane (CH ₄ ), ethane (C ₂ H ₆ ), propane (C ₃ H ₈ ), butane (C ₄ H ₁₀ ), nitrogen (N ₂ ), and It contained any or all of carbon dioxide (CO ₂ ). For example, one sample gas mixture contained 90 vol% methane, 3 vol% ethane, 1 vol% propane, 1 vol% butane, 4 vol% nitrogen, and 1 vol% carbon dioxide. Also, some sample gas mixtures contained 85 vol% methane, 10 vol% ethane, 3 vol% propane, and 2 vol% butane, and did not contain nitrogen and carbon dioxide. One sample gas mixture also contained 85 vol% methane, 8 vol% ethane, 2 vol% propane, 1 vol% butane, 2 vol% nitrogen, and 2 vol% carbon dioxide.

Next, using each of the 12 kinds of sample mixed gases, a plurality of measured values of the electric signal S _P from the pressure sensor and the electric signal S _H1 (T _H1 ) from the heating element to which four different voltages are applied. , S _H2 (T _H2 ), S _H3 (T _H3 ), and a plurality of measured values of S _H4 (T _H4 ). Thereafter, the known calorific value Q of the 12 kinds of sample mixed gas, a plurality of measured values of the electric signal S _P from the pressure sensor, and the electric signals S _H1 (T _H1 ), S _H2 ( Based on a plurality of measured values of T _H2 ), S _H3 (T _H3 ), and S _H4 (T _H4 ), an electric signal S _P from the pressure sensor and an electric signal S from the heating element are obtained by support vector regression. An equation for calculating the calorific value Q, where _H1 (T _H1 ), S _H2 (T _H2 ), S _H3 (T _H3 ), and S _H4 (T _H4 ) are independent variables and the calorific value Q is a dependent variable, Created.

  Using the created equation, the calorific value Q of each of the 12 sample mixed gases was calculated and compared with the true calorific value Q. As shown in FIG. 11, the error was within ± 1%. .

(Comparative example)
Using the same 12 kinds of sample mixed gases as in the embodiment, the electric signals S _H1 (T _H1 ), S _H2 (T _H2 ), S from the heating elements are not included without including the independent variable of the electric signal S _P from the pressure sensor. An equation for calculating the calorific value Q was created, with _H3 (T _H3 ) and S _H4 (T _H4 ) as independent variables and the calorific value Q as a dependent variable. When the calorific value Q of each of the 12 sample mixed gases is calculated using an equation that does not include the independent variable of the electrical signal S _P from the created pressure sensor, and compared with the true calorific value Q, FIG. As shown, the error was within ± 2%.

8 Microchip 18 Heat insulation member 20 Heat generation amount measurement system 31A, 31B, 31C, 31D Gas pressure regulators 32A, 32B, 32C, 32D Flow control devices 50A, 50B, 50C, 50D Gas cylinder 60 Substrate 61 Heating element 62 First measurement Temperature element 63 Second temperature measuring element 64 Thermal insulation element 65 Insulating film 66 Cavity 91A, 91B, 91C, 91D, 92A, 92B, 92C, 92D, 93, 102, 103 Channel 101 Chamber 161, 162, 163, 164 165, 261, 264, 265 Resistance element 170, 270 Operational amplifier 201 Pressure sensor 301 Measurement unit 302 Formula creation unit 303 Drive circuit 304 A / D conversion circuit 305 Heat generation amount calculation unit 312 Input device 313 Output device 401 Electric signal storage device 402 Formula Storage device 403 Heat generation amount storage device

According to the aspect of the present invention, (a) a heating element to which a plurality of different voltages can be applied, a pressure sensor , and a temperature measuring element are provided , and a container into which gas is injected, and (b) a plurality of different voltages are respectively added. The heat signal with the heat generation amount as a dependent variable, with the electrical signal from the heating element, the electrical signal from the pressure sensor, and the electrical signal from the temperature measuring element depending on the temperature of the gas injected into the container as independent variables A formula storage device that stores the quantity calculation formula, and (c) the measured value of the electric signal from the heating element , pressure sensor , and temperature measurement element is substituted into the independent variable of the calorific value calculation formula , and injected into the container A calorific value measuring system is provided, comprising: a calorific value calculating unit that calculates a measured value of the calorific value of the gas.

According to the aspect of the present invention, (a) injecting gas into a container provided with a heating element, a pressure sensor , and a temperature measuring element to which a plurality of different voltages can be applied ; and (b) a plurality of different voltages. The electrical signal from the heating element, the pressure sensor , and the electrical signal from the temperature measuring element depending on the temperature of the gas injected into the container are set as independent variables, and the amount of generated heat is a dependent variable. (C) Substituting the measured values of the electrical signals from the heating element , pressure sensor , and temperature measuring element into the independent variables of the heating value calculation formula , and injecting them into the container And a calorific value measurement method comprising: calculating a calorific value measurement value of the gas.

Moreover, according to the aspect of the present invention, (a) a measurement unit that measures a measured value of a heat dissipation coefficient or thermal conductivity of a gas, a measured value of a gas pressure, and a measured value of a gas temperature ; b) a storage device that stores the correlation between the heat dissipation coefficient or thermal conductivity, pressure, temperature, and calorific value; (c) the measured value of the gas heat dissipation coefficient or thermal conductivity, and the gas pressure. A calorific value measurement system is provided that includes a calorific value calculation unit that calculates a measured value of the calorific value of gas based on the measured value, the measured value of gas temperature, and the correlation.

Further, according to the aspect of the present invention, (a) measuring a measured value of a heat dissipation coefficient or thermal conductivity of a gas, a measured value of a gas pressure, and a measured value of a gas temperature; ) Prepare correlation between heat dissipation coefficient or thermal conductivity, pressure, temperature and calorific value, (c) measured value of gas heat dissipation coefficient or thermal conductivity, measured value of gas pressure And calculating a measurement value of the calorific value of the gas based on the measured value of the gas temperature and the correlation.

Claims (24)

A heating element to which a plurality of different voltages can be applied, and a pressure sensor, and a container into which gas is injected,
A storage device for storing a calorific value calculation formula having the electric signal from the heating element and the electric signal from the pressure sensor when the plurality of different voltages are applied as independent variables and the calorific value as a dependent variable; ,
A calorific value calculation unit that calculates a measured value of the calorific value of the gas injected into the container by substituting the measured value of the electrical signal from the heating element and the pressure sensor into the independent variable of the calorific value calculation formula; ,
A calorific value measuring system. The container is further provided with a temperature measuring element,
The calorific value calculation formula further includes, as an independent variable, an electrical signal from the temperature measuring element that depends on the temperature of the gas injected into the container,
The calorific value calculation unit further substitutes the measured value of the electrical signal from the temperature measuring element into the independent variable of the calorific value calculation formula to calculate the measured value of the calorific value of the gas injected into the container. ,
The calorific value measurement system according to claim 1.   The calorific value measurement system according to claim 1 or 2, wherein the number of the plurality of different voltages is at least one obtained by subtracting 1 from the number of the plurality of types of gas components contained in the gas.   The calorific value based on the value of the calorific value of the plurality of sample mixed gases containing plural kinds of gas components and the value of the electric signal from the heating element and the pressure sensor in contact with each of the plural sample mixed gases The calorific value measurement system according to any one of claims 1 to 3, wherein a calculation formula is created.   The calorific value measurement system according to claim 4, wherein support vector regression is used to create the calorific value calculation formula.   The calorific value measurement system according to claim 4 or 5, wherein each of the plurality of sample mixed gases is natural gas.   The calorific value measurement system according to any one of claims 1 to 6, wherein the gas is natural gas. Injecting gas into a container provided with a heat generating element to which a plurality of different voltages can be applied, and a pressure sensor;
Preparing a calorific value calculation formula with the electric signal from the heating element and the electric signal from the pressure sensor when the plurality of different voltages are respectively applied as independent variables and the calorific value as a dependent variable;
Substituting the measured value of the electrical signal from the heating element and the pressure sensor into the independent variable of the calorific value calculation formula, calculating the measured value of the calorific value of the gas injected into the container;
Method for measuring calorific value, including The container is further provided with a temperature measuring element,
The calorific value calculation formula further includes, as an independent variable, an electrical signal from the temperature measuring element that depends on the temperature of the gas injected into the container,
In calculating the calorific value measurement value, the measured value of the electrical signal from the temperature measuring element is further substituted into the independent variable of the calorific value calculation formula,
The calorific value measuring method according to claim 8.   The calorific value measurement method according to claim 8 or 9, wherein the number of the plurality of different voltages is at least one obtained by subtracting 1 from the number of the plurality of types of gas components contained in the gas.   The calorific value based on the value of the calorific value of the plurality of sample mixed gases containing plural kinds of gas components and the value of the electric signal from the heating element and the pressure sensor in contact with each of the plural sample mixed gases The calorific value measurement method according to any one of claims 8 to 10, wherein a calculation formula is created.   The calorific value measurement method according to claim 11, wherein support vector regression is used to create the calorific value calculation formula.   The calorific value measurement method according to claim 11 or 12, wherein each of the plurality of sample mixed gases is natural gas.   The calorific value measuring method according to any one of claims 8 to 13, wherein the gas is natural gas. A measurement unit for measuring the measured value of the heat dissipation coefficient or thermal conductivity of the gas, and the measured value of the pressure of the gas;
A storage device that stores a correlation between the heat dissipation coefficient or thermal conductivity, the pressure, and the heat generation amount;
A calorific value calculation unit that calculates a measured value of the calorific value of the gas based on the measured value of the heat dissipation coefficient or thermal conductivity of the gas, the measured value of the gas pressure, and the correlation;
A calorific value measuring system.   The value of the calorific value of a plurality of sample mixed gases containing a plurality of types of gas components, the value of the heat dissipation coefficient or thermal conductivity of each of the plurality of sample mixed gases, and the value of the pressure of each of the plurality of sample mixed gases The calorific value measurement system according to claim 15, wherein a calorific value calculation formula representing the correlation is created based on   The calorific value measurement system according to claim 16, wherein support vector regression is used to create the calorific value calculation formula.   The calorific value measurement system according to claim 16 or 17, wherein each of the plurality of sample mixed gases is natural gas.   The calorific value measurement system according to any one of claims 15 to 18, wherein the gas is natural gas. Measuring the measured value of the heat dissipation coefficient or thermal conductivity of the gas and the measured value of the pressure of the gas;
Preparing a correlation between the heat dissipation coefficient or thermal conductivity, the pressure, and the amount of heat generated;
Calculating a measured value of the calorific value of the gas based on a measured value of the heat dissipation coefficient or thermal conductivity of the gas, a measured value of the pressure of the gas, and the correlation;
Method for measuring calorific value, including   The value of the calorific value of a plurality of sample mixed gases containing a plurality of types of gas components, the value of the heat dissipation coefficient or thermal conductivity of each of the plurality of sample mixed gases, and the value of the pressure of each of the plurality of sample mixed gases The calorific value calculation method according to claim 20, wherein a calorific value calculation formula representing the correlation is created based on   The calorific value measurement method according to claim 21, wherein support vector regression is used to create the calorific value calculation formula.   The calorific value measurement method according to claim 21 or 22, wherein each of the plurality of sample mixed gases is natural gas.   The calorific value measurement method according to any one of claims 20 to 23, wherein the gas is natural gas.