Natural gas physical sensible heat calculation method in programmable gas turbine generator set heat rate calculation

Abstract

The invention discloses a natural gas physical sensible heat calculation method in programmable gas turbine generator set heat rate calculation. The natural gas physical sensible heat calculation method comprises the following steps that the natural gas temperature is collected; measuring natural gas components to obtain a gas molar component matrix; setting a natural gas combustion reference temperature under the test condition; establishing a gas constant matrix; calculating a single-gas specific heat capacity function matrix based on the gas constant matrix; converting the gas molar component matrix into a gas mass component matrix, and obtaining a mixed gas specific heat capacity function based on the single gas specific heat capacity function matrix and the gas mass component matrix; integrating the mixed gas specific heat capacity function to obtain natural gas physical sensible heat, wherein the integral upper limit is natural gas temperature and the integral lower limit is natural gas combustion reference temperature; the natural gas physical sensible heat is used for calculating the heat consumption rate of the gas turbine generator set. The method reduces the calculation error of the physical sensible heat of the natural gas, and improves the calculation accuracy of the heat consumption rate of the gas turbine generator set.

Description

Technical Field

[0001] This invention relates to the field of performance testing of gas turbine generator sets, and more particularly to a programmable method for calculating the physical sensible heat of natural gas in the calculation of the heat rate of gas turbine generator sets. Background Technology

[0002] Gas turbine generator sets (gas-steam combined cycle units) are cleaner, more efficient, and more flexible than traditional coal-fired power units, and they account for a certain proportion of installed power capacity. Heat consumption rate is an important energy consumption indicator for these units.

[0003] In the calculation of heat rate during the performance test of gas turbine generator sets, it is necessary to calculate the change in physical sensible heat of natural gas compared to the test boundary conditions. A rougher calculation method is to simply determine the specific heat capacity of natural gas by looking up a table, and then calculate the change in physical sensible heat based on the temperature difference. A more accurate calculation method is to look up the specific heat capacity of the corresponding gas components in the API technical manual, then calculate the overall specific heat capacity, and finally calculate the change in physical sensible heat based on the temperature difference.

[0004] When using the above methods, the specific heat capacity of the gas is often calculated at the average temperature, ignoring the effect of temperature on the specific heat capacity of natural gas. At the same time, the fitting formula in the API technical manual is not the most accurate formula for calculating specific heat capacity, considering the convenience of calculation. A more accurate formula would lead to a more large and complex calculation. Summary of the Invention

[0005] To address the shortcomings of existing technologies, this invention adopts a more accurate formula for calculating the specific heat capacity of gas derived from the national standard GB / T 30491.1, and establishes a programmable method for calculating the physical sensible heat of natural gas in the calculation of the heat rate of gas turbine generator sets. By integrating and summing, the influence of temperature on specific heat capacity is incorporated into the calculation, while solving the problem of the cumbersome nature of manual calculation and improving the accuracy and convenience of related calculations.

[0006] To achieve the above objectives, the present invention adopts the following technical solution:

[0007] A method for calculating the physical sensible heat of natural gas in the calculation of the heat rate of a programmable gas turbine generator set includes the following steps:

[0008] Natural gas temperature was collected.

[0009] The composition of natural gas was measured to obtain a gas molar composition matrix.

[0010] Set the reference temperature for natural gas combustion under the test conditions;

[0011] Establish the gas constant matrix;

[0012] Calculation of the specific heat capacity function matrix of a single gas based on the gas constant matrix;

[0013] The gas molar composition matrix is ​​converted into a gas mass composition matrix, and the specific heat capacity function of the mixed gas is obtained based on the single gas specific heat capacity function matrix and the gas mass composition matrix.

[0014] The physical sensible heat of natural gas is obtained by integrating the specific heat capacity function of the mixed gas. The upper limit of the integration is the temperature of natural gas, and the lower limit of the integration is the reference temperature for natural gas combustion. The physical sensible heat of natural gas is used to calculate the heat rate of the gas turbine generator set.

[0015] To optimize the above technical solution, the specific measures also include:

[0016] Furthermore, the specific temperature of the natural gas collected is:

[0017] The gas turbine generator set heat rate test was conducted under stable operating conditions for 20 minutes, with automatically collected data at 30-second intervals. The natural gas temperature was taken as the average value (T) of the temperature data collected during the test. t .

[0018] Furthermore, the measurement of natural gas components specifically involves:

[0019] m natural gas samples were collected, and the natural gas components were measured using a gas chromatograph equipped with a calibrated standard gas source. The molar content of n gas components was obtained for each natural gas sample.

[0020] Furthermore, the gas molar composition matrix is ​​calculated as follows:

[0021] The molar content of the i-th gaseous component is calculated using the following formula:

[0022]

[0023] In equation (1), Z i The molar content of the i-th gaseous component is expressed in mol / mol. This represents the molar content of the i-th gas component in the j-th natural gas sample; n is the total number of all gas components in the natural gas component report; and m is the number of natural gas samples.

[0024] Twenty-one gaseous components were detected, namely CH4, C2H6, C3H8, and C4H. 10 i-C4H 10 C5H 12 i-C5H 12 C6H 14 C7H 16 C8H 18 C9H 20 C 10H 22 N2, CO2, H2S, H2O, O2, H2, He, Ar and CO;

[0025] The molar content of each gas component forms a 21-row × 1-column data matrix as follows:

[0026]

[0027] In the formula, Z represents the gas molar composition matrix.

[0028] Furthermore, the reference temperature for natural gas combustion is set to 0°C, 15°C, 15.5°C, 20°C, or 25°C.

[0029] Furthermore, the establishment of the gas constant matrix specifically involves:

[0030] Gas constant matrix C all The structure is 21 rows × 10 columns, as follows:

[0031]

[0032] In equation (3), C i_1 ~C i_9 M represents the constant for calculating the specific heat capacity of the i-th gas. i Let be the molar mass of the i-th gas.

[0033] Furthermore, the calculation of the single-gas specific heat capacity function matrix based on the gas constant matrix specifically involves:

[0034] The gas constant matrix C all Substitute the specific heat capacity constant and the molar mass of each gas into the following formula for calculating the specific heat capacity of gases:

[0035] C pi ={C i_1 +C i_2 ×[C i_3 / T / sinh(C i_2 / T)] 2 +C i_3 ×[C i_5 / T / cosh(C i_5 / T)] 2 +C i_6 ×[C i_7 / T / sinh(C i_7 / T)] 2 +C i_8 ×[C i_9 / T / cosh(C i_9 / T)] 2}×R / M i (4)

[0036] In the formula, C pi C represents the specific heat capacity of the i-th gas, expressed in kJ / (kg·K); i_1 ~C i_9 M represents the constant for calculating the specific heat capacity of the i-th gas. i R is the molar mass of the i-th gas, in kg / kmol; R is the molar gas constant, in kJ / (kmol·K); T is the Kelvin temperature of the natural gas, in K.

[0037] The specific heat capacity function C of 21 gases pi The 21-row × 1-column matrix of specific heat capacity functions for a single gas is as follows:

[0038]

[0039] In the formula, C p This is the specific heat capacity function matrix for a single gas.

[0040] Furthermore, the conversion of the gas molar composition matrix into a gas mass composition matrix specifically involves:

[0041]

[0042] In the formula, G i Z represents the mass of the i-th gas, and Z(i,1) represents the element in the i-th row and 1-th column of the gas molar composition matrix Z, that is, the molar content of the i-th gas component Z. i C all (i, 10) represents the gas constant matrix C. all The element in the i-th row and 10-th column, i.e., the gas molar mass M of the i-th gas component. i C all (:,10) represents the gas constant matrix C. all The vector consisting of the elements in the 10th column, where the superscript T indicates transpose, Z represents the gas molar composition matrix, and G represents the gas mass composition matrix.

[0043] Furthermore, the specific method for obtaining the specific heat capacity function of the mixed gas based on the single-gas specific heat capacity function matrix and the gas mass component matrix is ​​as follows:

[0044] C pnat =G T ×C p (8)

[0045] In the formula, C pnat Let G represent the specific heat capacity function of the gas mixture, and C represent the mass composition matrix of the gas. p This is a single-gas specific heat capacity function matrix, with the superscript T indicating transpose.

[0046] Furthermore, the method of integrating the specific heat capacity function of the mixed gas to obtain the physical sensible heat of natural gas specifically involves:

[0047]

[0048] In equation (9), Q represents the sensible heat of natural gas used in the heat rate calculation, in kJ / kg, and T represents the sensible heat of natural gas. t T represents the average value of the natural gas temperature data collected during the test. s This represents the reference temperature for natural gas combustion, where T is the integral variable, i.e., the Kelvin temperature of the natural gas; C pnat This is a function of the specific heat capacity of the gas mixture;

[0049] For integral calculation, select the integral calculation tool from the program function library;

[0050] Alternatively, an approximate calculation can be performed using Simpson's formula, with the following formula divided into 6 sub-intervals:

[0051]

[0052] In the formula, C pnat (·) represents the specific heat capacity function of the mixed gas.

[0053] The beneficial effects of this invention are as follows: Based on the standard gas specific heat capacity calculation formula and integral calculation method, this invention provides a method for calculating the physical sensible heat of natural gas in the heat consumption calculation of gas turbine generator sets that can be used for programmed calculation. This avoids the tedium of manual calculation, reduces the error of natural gas physical sensible heat calculation, and improves the accuracy of gas turbine generator set heat consumption rate calculation. Attached Figure Description

[0054] Figure 1 This is an overall flowchart of a programmable gas turbine generator set heat rate calculation method for calculating the physical sensible heat of natural gas. Detailed Implementation

[0055] The technical solutions in the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings and examples of MATLAB compilers. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. All other embodiments obtained by those skilled in the art based on the embodiments of this application without creative effort are within the scope of protection of this application.

[0056] This invention proposes a method for calculating the physical sensible heat of natural gas in the calculation of the heat rate of a programmable gas turbine generator set, such as... Figure 1 As shown, it includes the following steps:

[0057] Natural gas temperature was collected.

[0058] The composition of natural gas was measured to obtain a gas molar composition matrix.

[0059] Set the reference temperature for natural gas combustion under the test conditions;

[0060] Establish the gas constant matrix;

[0061] Calculation of the specific heat capacity function matrix of a single gas based on the gas constant matrix;

[0062] The gas molar composition matrix is ​​converted into a gas mass composition matrix, and the specific heat capacity function of the mixed gas is obtained based on the single gas specific heat capacity function matrix and the gas mass composition matrix.

[0063] The physical sensible heat of natural gas is obtained by integrating the specific heat capacity function of the mixed gas. The upper limit of the integration is the temperature of natural gas, and the lower limit of the integration is the reference temperature for natural gas combustion. The physical sensible heat of natural gas is used to calculate the heat rate of the gas turbine generator set.

[0064] The specific process is as follows:

[0065] First, performance tests were conducted on the gas turbine generator set or gas-steam combined cycle unit according to relevant standards. The unit was kept under stable operating conditions for 20 minutes. Natural gas temperature data was recorded using an online data acquisition system, with one set every 30 seconds, for a total of 40 sets of data. The average of the 40 sets of data was calculated to obtain the natural gas temperature data T. t =30℃ =303.15K. Create a variable in the program and input it as follows.

[0066] Tt0 = input('Please enter the natural gas test temperature (unit: °C):'); % Unit: °C

[0067] Tt = Tt0 + 273.15; % Unit: K

[0068] During the test, an online gas chromatograph equipped with a calibrated standard gas source was used to measure the gas components, obtaining 40 sets of data. Alternatively, natural gas samples were taken within 5 minutes before the test, within 5–15 minutes after the test began, and within 5 minutes after the test ended. These samples were then sent to a professional testing institution for natural gas component analysis within 12 hours of sampling, obtaining 3 sets of data. The average value of each component was calculated using the following formula.

[0069]

[0070] In equation (1), Z i The molar content of the i-th gas component is expressed in mol / mol; the superscript j indicates the j-th component data; n represents the total number of all gas components in the natural gas component report; m represents the number of natural gas samples; if online gas chromatograph is used to collect data, m = 40; if sampling and testing are used to collect data, m = 3.

[0071] Create variables in the program and input the natural gas components, as follows.

[0072] Z = zeros(21,1);

[0073] gasname={'CH4','C2H6','C3H8','C4H 10 ','i-C4H 10 ','C5H 12 ','i-C5H 12 ','C6H 14 ','C7H 16 ','C8H 18 ','C9H 20 ','C 10 H 22 ','N2','CO2','H2S','H2O','O2','H2','He','Ar','CO'};

[0074] for i = 1:1:21

[0075] str = {'Please enter the gas', gasname{i}, 'Components (unit: mol / mol):'};

[0076] Z(i,1) = input(str);

[0077] end

[0078] The matrix Z, which stores data for 21 gas components, is obtained.

[0079]

[0080] Subscripts 1-21 correspond to the gases CH4, C2H6, C3H8, and C4H. 10 (n-Butane), i-C4H 10 (isobutane), C5H 12 (n-pentane), i-C5H 12 (isopentane), C6H 14 C7H 16 C8H 18 C9H 20 C 10 H 22 , N2, CO2, H2S, H2O, O2, H2, He, Ar, CO.

[0081] Natural gas combustion reference temperature T sThe temperature should be selected according to the test protocol. If the protocol does not specify a temperature, the following commonly used combustion reference temperatures can be selected: 0℃, 15℃, 15.5℃ (60℉), 20℃, and 25℃. Establish and input the variables in the program as follows. Obtain the natural gas combustion reference temperature T. s =15℃=288.15K.

[0082] Ts0 = input('Please enter the reference temperature for natural gas combustion (unit: °C):'); % Unit: °C

[0083] Ts = Ts0 + 273.15; % Unit: K

[0084] The gas constants are obtained by looking up the table, and a gas constant matrix is ​​established as follows.

[0085]

[0086] The gas constant matrix is ​​established in the program as follows.

[0087] CH4=[4.00088 0.76315 820.659 0.00460 178.410 8.74432 1062.820 -4.4691090.530 16.043];

[0088] C2H6=[4.00263 4.33939 559.314 1.23722 223.284 13.19740 1031.380 -6.02 1071.290 30.07];

[0089] C3H8=[4.02939 6.60569 479.856 3.19700 200.893 19.19210 0955.312 -8.373 1027.29044.097];

[0090] C4H10=[4.33944 9.44893 468.270 6.89406 183.636 24.46180 1914.10014.78 903.18558.124];

[0091] iC4H10=[4.06714 8.97575 438.270 5.25156 198.018 25.14230 1905.02016.14 893.76558.124];

[0092] C5H12=[4 8.95043 178.670 21.83600 840.538 33.40320 1774.250 0 072.151];

[0093] iC5H12=[4 11.76180 292.503 20.11010 910.237 33.16880 1919.370 0 072.151];

[0094] C6H14=[4 11.69770 182.326 26.81420 859.207 38.61640 1826.590 0 086.178];

[0095] C7H16=[4 13.72660 169.789 30.47070 836.195 43.55610 1760.460 0 0100.205];

[0096] C8H18=[4 15.68650 158.922 33.80290 815.064 48.17310 1693.070 0 0114.232];

[0097] C9H20=[4 18.02410 156.854 38.12350 814.882 53.34150 1693.790 0 0128.259];

[0098] C10H22=[4 21.00690 164.947 43.49310 836.264 58.36570 1750.240 0 0142.286];

[0099] N2=[3.50031 0.13732 662.738 -0.14660 680.562 0.90066 1740.060 0 028.013];

[0100] CO2=[3.50002 2.04452 919.306 -1.06044 865.070 2.03366 483.553000.014 341.109 44.01];

[0101] H2S=[4 3.11942 1833.63 1.00243 847.181 0 1 0 0 34.076];

[0102] H2O=[4.00392 0.01059 268.795 0.98763 1141.410 3.06904 2507.370 0 018.015];

[0103] O2=[3.50146 1.07558 2235.71 1.01334 1116.690 0 1 0 0 31.999];

[0104] H2=[2.47906 0.95806 228.734 0.45444 326.843 1.56039 1651.710 -1.3761671.690 2.016];

[0105] He=[2.5 0.00000 1 0 0 0 1 0 0 4.0026];

[0106] Ar=[2.5 0.00000 1 0 0 0 1 0 0 39.948];

[0107] CO=[3.50055 1.02865 1550.450 0.00493 704.525 0 1 0 0 28.01];

[0108] Call=[CH4;

[0109] The formula for calculating the specific heat capacity of a gas, derived from national standards GB / T 30491.1-2014 and ISO 20765-1:2005, is as follows:

[0110] C p ={C1+C2×[C3 / T / sinh(C3 / T)] 2 +C4×[C5 / T / cosh(C5 / T)] 2 +C6×[C7 / T / sinh(C7 / T)] 2 +C8×[C9 / T / cosh(C9 / T)] 2}×R / M(4)

[0111] The calculation formula is established in the program as follows.

[0112] syms Cp C1 C2 C3 C4 C5 C6 C7 C8 C9 RM;

[0113] Cp=(C1+C2*((C3 / T) / sinh(C3 / T))^2+C4*((C5 / T) / cosh(C5 / T))^2+C6*((C7 / T) / sinh(C7 / T))^2+C8*((C9 / T) / cosh(C9 / T))^2)*(R / M);

[0114] R = 8.31451; % Unit: kJ / (kmol·K)

[0115] To calculate the specific heat capacity of a single-component gas, substitute the gas constant matrix corresponding to equation (3) into equation (4) to obtain the specific heat capacity function C for all gases. pi (T), forming a single-gas specific heat capacity function matrix C p (T) is as follows.

[0116]

[0117] The specific heat capacity function matrix of a single gas is established in the program as follows.

[0118] CpCH4=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],CH4);

[0119] CpC2H6=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],C2H6);

[0120] CpC3H8=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],C3H8);

[0121] CpC4H10=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],C4H10);

[0122] CpiC4H10=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],iC4H10);

[0123] CpC5H12=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],C5H12);

[0124] CpiC5H12=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],iC5H12);

[0125] CpC6H14=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],C6H14);

[0126] CpC7H16=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],C7H16);

[0127] CpC8H18=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],C8H18);

[0128] CpC9H20=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],C9H20);

[0129] CpC10H22=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],C10H22);

[0130] CpN2=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],N2);

[0131] CpCO2=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],CO2);

[0132] CpH2S=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],H2S);

[0133] CpH2O=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],H2O);

[0134] CpO2=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],O2);

[0135] CpH2=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],H2);

[0136] CpHe=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],He);

[0137] CpAr=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],Ar);

[0138] CpCO=subs(Cp,[C1 C2 C3 C4 C5 C6 C7 C8 C9 M],CO);

[0139] Cpall={CpCH4 CpC2H6 CpC3H8 CpC4H10 CpiC4H10 CpC5H12 CpiC5H12CpC6H14CpC7H16 CpC8H18 CpC9H20 CpC10H22 CpN2 CpCO2 CpH2S CpH2O CpO2 CpH2 CpHeCpAr CpCO}';

[0140] The specific process for calculating the specific heat capacity function of natural gas is as follows.

[0141] Step 1: Convert the gas molar component Z to the mass component G. The calculation formula is as follows.

[0142]

[0143] In the formula, G i Z represents the mass of the i-th gas, and Z(i,1) represents the element in the i-th row and 1-th column of the gas molar composition matrix Z, that is, the molar content of the i-th gas component Z. i C all (i, 10) represents the gas constant matrix C. all The element in the i-th row and 10-th column, i.e., the gas molar mass M of the i-th gas component. i C all (:,10) represents the gas constant matrix C. all The vector consisting of the elements in the 10th column, where the superscript T indicates transpose, Z represents the gas molar composition matrix, and G represents the gas mass composition matrix.

[0144] Step 2: Calculate the specific heat capacity function of the gas mixture. The calculation formula is as follows.

[0145] C pnat =G T ×C p (8)

[0146] The first step of the calculation in the program is as follows.

[0147] G = zeros(21,1);

[0148] for i = 1:1:21

[0149] G(i,1)=Z(i,1)*Call(i,10) / (Call(:,10)'*Z);

[0150] end

[0151] The second step of the calculation in the program is as follows.

[0152] syms Cp1 Cp2;

[0153] Cp1=eval(Cpall{1,1}*G(1,1));

[0154] for i = 2:1:21

[0155] Cp2=eval(Cp1+eval(Cpall{i,1}*G(i,1)));

[0156] Cp1 = Cp2;

[0157] end

[0158] The physical sensible heat of natural gas is calculated using an integral method, and the formula is as follows.

[0159]

[0160] In equation (9), Q is the sensible heat of natural gas involved in the heat rate calculation, in kJ / kg.

[0161] The physical sensible heat of natural gas is calculated in the program as follows.

[0162] Q = eval(int(Cp1,T,Ts,Tt));

[0163] disp(Q);

[0164] After the program completes execution, the specific heat capacity of natural gas under any temperature condition can be calculated via commands:

[0165] Tr = input('Please enter the natural gas temperature (200K~1000K) to query the specific heat capacity:');

[0166] Cp = eval(subs(Cp1,T,Tr));

[0167] disp(Cp);

[0168] The physical sensible heat is then used for the next step of heat rate calculation.

[0169] The final calculation results are shown below.

[0170] Input variables:

[0171] Natural gas test temperature: T t =30℃

[0172] Natural gas combustion reference temperature: T s =15℃

[0173] The composition of natural gas is shown in Table 1:

[0174] Table 1

[0175]

[0176]

[0177] Output result:

[0178] The sensible heat of natural gas is Q = 32.2512 kJ / kg.

[0179] You can view the specific heat capacity of natural gas at any temperature from 200K to 1000K, for example:

[0180] Cp=eval(subs(Cp1,T,293.15));

[0181] The specific heat capacity at 293.15 K is 2.1437 kJ / (kg·K).

[0182] You can view the function for calculating the specific heat capacity of natural gas with respect to temperature T under the current gas composition:

[0183]

[0184] As can be seen, the calculation formula is extremely complex, and this invention effectively reduces the tediousness of the related calculations.

[0185] In addition to using the integration tools in the program's function library, Simpson's formula can also be used for approximate calculation. The calculation formula, divided into 6 sub-intervals, is as follows:

[0186]

[0187] In the formula, C pnat (·) represents the specific heat capacity function of the mixed gas.

[0188] Those skilled in the art will recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed in this application can be implemented in electronic hardware or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0189] The above are merely preferred embodiments of the present invention. The scope of protection of the present invention is not limited to the above embodiments. All technical solutions falling within the scope of the present invention's concept are within the scope of protection of the present invention. It should be noted that for those skilled in the art, any improvements and modifications made without departing from the principles of the present invention should be considered within the scope of protection of the present invention.

Claims

1. A method for calculating the physical sensible heat of natural gas in the calculation of the heat rate of a programmable gas turbine generator set, characterized in that, Includes the following steps: Natural gas temperature was collected. The composition of natural gas was measured to obtain a gas molar composition matrix. Set the reference temperature for natural gas combustion under the test conditions; Establish the gas constant matrix; Calculation of the specific heat capacity function matrix of a single gas based on the gas constant matrix; The gas molar composition matrix is ​​converted into a gas mass composition matrix, and the specific heat capacity function of the mixed gas is obtained based on the single gas specific heat capacity function matrix and the gas mass composition matrix. The physical sensible heat of natural gas is obtained by integrating the specific heat capacity function of the mixed gas. The upper limit of the integration is the temperature of natural gas, and the lower limit of the integration is the reference temperature for natural gas combustion. The physical sensible heat of natural gas is used to calculate the heat rate of the gas turbine generator set.

2. The method for calculating the physical sensible heat of natural gas in the programmable gas turbine generator set heat rate calculation as described in claim 1, characterized in that, The specific temperature of the natural gas collected is: The gas turbine generator set heat rate test was conducted under stable operating conditions for 20 minutes, with automatically collected data at 30-second intervals. The natural gas temperature was taken as the average value (T) of the temperature data collected during the test. t .

3. The method for calculating the physical sensible heat of natural gas in the calculation of the heat rate of a programmable gas turbine generator set as described in claim 1, characterized in that, The measurement of natural gas components specifically involves: m natural gas samples were collected, and the natural gas components were measured using a gas chromatograph equipped with a calibrated standard gas source. The molar content of n gas components was obtained for each natural gas sample.

4. The method for calculating the physical sensible heat of natural gas in the programmable gas turbine generator set heat rate calculation as described in claim 1, characterized in that, The gas molar composition matrix is ​​calculated as follows: The molar content of the i-th gaseous component is calculated using the following formula: In equation (1), Z i The molar content of the i-th gaseous component is expressed in mol / mol. This represents the molar content of the i-th gas component in the j-th natural gas sample; n is the total number of all gas components in the natural gas component report; and m is the number of natural gas samples. Twenty-one gaseous components were detected, namely CH4, C2H6, C3H8, and C4H. 10 i-C4H 10 C5H 12 i-C5H 12 C6H 14 C7H 16 C8H 18 C9H 20 C 10 H 22 N2, CO2, H2S, H2O, O2, H2, He, Ar and CO; The molar content of each gas component forms a 21-row × 1-column data matrix as follows: In the formula, Z represents the gas molar composition matrix.

5. The method for calculating the physical sensible heat of natural gas in the calculation of the heat rate of a programmable gas turbine generator set as described in claim 1, characterized in that, The reference temperature for natural gas combustion is set to 0℃, 15℃, 15.5℃, 20℃, or 25℃.

6. The method for calculating the physical sensible heat of natural gas in the programmable gas turbine generator set heat rate calculation as described in claim 1, characterized in that, The establishment of the gas constant matrix is ​​specifically as follows: Gas constant matrix C all The structure is 21 rows × 10 columns, as follows: In equation (3), C i_1 ~C i_9 M represents the constant for calculating the specific heat capacity of the i-th gas. i Let be the molar mass of the i-th gas.

7. The method for calculating the physical sensible heat of natural gas in the programmable gas turbine generator set heat rate calculation as described in claim 6, characterized in that, The specific calculation of the single-gas specific heat capacity function matrix based on the gas constant matrix is ​​as follows: The gas constant matrix C all Substitute the specific heat capacity constant and the molar mass of each gas into the following formula for calculating the specific heat capacity of gases: C pi ={C i_1 +C i_2 ×[C i_3 / T / sinh(C i_3 / T)] 2 +C i_4 ×[C i_5 / T / cosh(C i_5 / T)] 2 +C i_6 ×[C i_7 / T / sinh(C i_7 / T)] 2 +C i_8 ×[C i_9 / T / cosh(C i_9 / T)] 2 }×R / M i (4) In the formula, C pi C represents the specific heat capacity of the i-th gas, expressed in kJ / (kg·K); i_1 ~C i_9 M represents the constant for calculating the specific heat capacity of the i-th gas. i R is the molar mass of the i-th gas, in kg / kmol; R is the molar gas constant, in kJ / (kmol·K); T is the Kelvin temperature of the natural gas, in K. The specific heat capacity function C of 21 gases pi The 21-row × 1-column matrix of specific heat capacity functions for a single gas is as follows: In the formula, C p This is the specific heat capacity function matrix for a single gas.

8. The method for calculating the physical sensible heat of natural gas in the calculation of the heat rate of a programmable gas turbine generator set as described in claim 1, characterized in that, The specific steps for converting the gas molar composition matrix into a gas mass composition matrix are as follows: In the formula, G i Z represents the mass of the i-th gas, and Z(i,1) represents the element in the i-th row and 1-th column of the gas molar composition matrix Z, that is, the molar content of the i-th gas component Z. i C all (i, 10) represents the gas constant matrix C. all The element in the i-th row and 10-th column, i.e., the gas molar mass M of the i-th gas component. i C all (:,10) represents the gas constant matrix C. all The vector consisting of the elements in the 10th column, where the superscript T indicates transpose, Z represents the gas molar composition matrix, and G represents the gas mass composition matrix.

9. The method for calculating the physical sensible heat of natural gas in the calculation of the heat rate of a programmable gas turbine generator set as described in claim 1, characterized in that, The specific method for obtaining the specific heat capacity function of the mixed gas based on the single-gas specific heat capacity function matrix and the gas mass component matrix is ​​as follows: C pnat =G T ×C p (8) In the formula, C pnat Let G represent the specific heat capacity function of the gas mixture, and C represent the mass composition matrix of the gas. p This is a single-gas specific heat capacity function matrix, with the superscript T indicating transpose.

10. The method for calculating the physical sensible heat of natural gas in the calculation of the heat rate of a programmable gas turbine generator set as described in claim 1, characterized in that, The specific method for obtaining the physical sensible heat of natural gas by integrating the specific heat capacity function of the mixed gas is as follows: In equation (9), Q represents the sensible heat of natural gas used in the heat rate calculation, in kJ / kg, and T represents the sensible heat of natural gas. t T represents the average value of the natural gas temperature data collected during the test. s This represents the reference temperature for natural gas combustion, where T is the integral variable, i.e., the Kelvin temperature of the natural gas; C pnat Let be the specific heat capacity function of the gas mixture; For integral calculation, select the integral calculation tool from the program function library; Alternatively, an approximate calculation can be performed using Simpson's formula, with the following formula divided into 6 sub-intervals: In the formula, C pnat (·) represents the specific heat capacity function of the gas mixture.

Images