Chemo-Metrical Prediction of Methane Index for the Natural Gas

Abstract

A new fuel quality sensor enables the use of RNG without removing CO.sub.2. Efficiency and the economy of the process improves significantly if CO.sub.2 separation cost can be avoided. Fuel quality information from the sensor makes it easy to adjust air/fuel ratio as well as efficient combustion inside a combustion engine or boiler.

Claims

1. A fuel quality sensor, comprising: a database of Wobbe index values, methane number values, thermal conductivity values, and sound velocities of a gas mixture; and a predictive model coupled to the database and a sensor to determine parameters of renewable nature gas and to adjust an air/fuel ratio for an engine based on the predictive model.

2. The system of claim 1, further comprising a data set builder for the predictive model that builds the database.

3. The system of claim 1, wherein the predictive model comprises a linear regression predictive model.

4. The system of claim 1, wherein the predictive model is to adjust the air/fuel ratio in real time during operation of the engine.

5. A method for sensing fuel quality, the method comprising: building a database of Wobbe index values, methane number values, thermal conductivity values, and sound velocities of a gas mixture; and applying a predictive model coupled to the database and a sensor to determine parameters of renewable nature gas and to adjust an air/fuel ratio for an engine based on the predictive model.

6. The method of claim 5, wherein the step of building the database comprises using a data set builder for the predictive model.

7. The method of claim 6, wherein the predictive model comprises a linear regression predictive model.

8. The method of claim 5, wherein the predictive model is to adjust the air/fuel ratio in real time during operation of the engine.

9. An onboard engine controller comprising: a non-resident memory containing a database of Wobbe index values, methane number values, thermal conductivity values, and sound velocities of a gas mixture; a gas sensor to provide real time parameters of the fuel input stream; a processor to run a predictive model coupled to the database and the sensor to determine the real-time parameters of the fuel input stream, process air flow data, and provide real-time air/fuel ratio adjustment commands based on the predictive model; a fuel control value to adjust the air/fuel ratio based on the adjustment command; whereby the air/fuel ratio is optimized for an engine based on the specific fuel provided.

10. The controller of claim 9, where the processor further hosts a data set builder for the predictive model that builds the database.

11. The controller of claim 9, where the processor executes a linear regression predictive model.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0018] These and other features, aspects and advantages of the present invention will become better understood with regard to the following description, appended claims, and accompanying figures where:

[0019] FIG. 1 is a high level schematic diagram that illustrates an apparatus for heating value calculation according to one embodiment;

[0020] FIG. 2 is a graph of the predicted Wobbe index verses the actual Wobbe index;

[0021] FIG. 3 is a graph of the predicted methane number verses the actual methane number;

[0022] FIG. 4 graph of the predicted CH.sub.4 composition verses the actual CH.sub.4 composition;

[0023] FIG. 5 is a graph of the predicted C.sub.2H.sub.6 composition verses the actual C.sub.2H.sub.6 composition; and

[0024] FIG. 6 is a graph of the predicted CO.sub.2 composition verses the actual CO.sub.2 composition.

[0025] FIG. 7 is an exemplary embodiment of a fuel system controller employing the current invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0026] For illustrative purposes, the present invention is embodied in the apparatus and method generally shown and described herein with reference to FIG. 1 through FIG. 6. It will be appreciated that the apparatus may vary as to configuration and as to details of the parts, and that the method may vary as to the specific steps and sequence, without departing from the basic concepts as disclosed herein.

Development of a Predictive Model

[0027] In order to develop a predictive model, a data set is developed. The system creates a data set that comprises attributes such as thermal conductivity, sound velocity, temperature, pressure, gas composition, WI, and MN. In this data set, the attributes, such as thermal conductivity, sound velocity, temperature, and pressure are physical properties that depend on the gas composition and are easily measured using widely available sensors or calculated through the gas mixture using software tools. The attributes, WI and MN also depend on the gas composition. Thus, in one embodiment, one method is to first measure the gas composition and calculate the WI and MN using the measured composition, which requires a long and expensive process. However, more preferably, the present system estimates WI, MN, and gas composition based on physical properties such as thermal conductivity, sound velocity, temperature, and pressure of the gas mixture in real-time. The main advantage of the system is that the physical properties can be easily and economically measured directly in the gas using sensors currently available in the market without analyzing the gas composition through the long and expensive process. A detailed process of preparing the data set and developing the predictive model to estimate the WI, MN, and gas composition will now be discussed.

Creating a Data Set

[0028] To create a data set with the right attributes of gas composition, fossil, anaerobic digester gas, and landfill gas are mixed together to get the various combination of gas mixture containing components CH.sub.4, CO.sub.2 and C.sub.2H.sub.6 by 10% incremental from 0% to 100%. Table 1 illustrates the composition of the components in Fossil natural gas, anaerobic digester gas, and landfill gas.

TABLE-US-00001 TABLE 1 Component of Gaseous Fossil Anaerobic Landfill Fuel Mixture Gas Digester Gas Gas Methane (mol %) 97 68 60 Carbon Dioxide (mol %) 0 26 33 Water (mol %) 0 5 6.5 Ethane (mol %) 2 0 0 Other (N2, O2) (mol %) 1 1 0.5 Percent of Volumes of Components in Normalized Gases Methane (mol %) 98 72.3 64.5 Carbon Dioxide (mol %) 0 27.7 35.5 Ethane (mol %) 2 0 0

[0029] For example, a mixture of 40% fossil, 30% anaerobic digester gas and 30% landfill gas contains 78.64% CH.sub.4, 18.96% CO.sub.2, and 2.4% C.sub.2H.sub.6. Using this method, by way of example, and not by way of limitation, sixty-six different combinations of the gas mixture is created for the gas composition. The combination of the gas mixture along with the composition of the components is shown in Table 2.

TABLE-US-00002 TABLE 2 Gaseous Fuel Mixture Composition (%) Fossil Anaerobic Landfill Composition (%) Natural Gas Disaster Gas Gas CH.sub.4 C.sub.2H.sub.6 CO.sub.2 100 0 0 94 6 0 90 10 0 91.83 5.4 2.77 80 0 20 88.1 4.8 7.1 70 20 10 86.71 4.2 9.09 60 30 10 84.54 3.6 11.86 50 30 20 81.59 3 15.41 40 40 20 79.42 2.4 18.18 30 20 50 74.91 1.8 23.29 20 50 30 74.3 1.2 24.5 10 60 30 72.13 0.6 27.27 10 20 70 69.01 0.6 30.39 0 100 0 72.3 0 27.7 0 0 100 64.5 0 35.5

[0030] In addition to the gas composition, pressure and temperature of the fuel are related to the thermal conductivity, sound velocity, WI, and MN. That indicates that thermal conductivity, sound velocity, WI, and MN may be determined by a function of temperature, pressure and composition of the gas mixture. Thus, combinations of temperature and pressure were created. For the combination of pressure, 500 psi intervals in the range between 500 psi and 3000 psi are used. For the combinations of temperature, 20 C. intervals in the range between 20 and 80 C. are used. From 6 different temperature intervals and 6 different pressure intervals, 36 possible combinations can be made. With these 36 combinations of temperature and pressure along, with 66 combinations of the gas mixture, the total of 2376 possible combinations were created.

[0031] The specific gas properties, such as thermal conductivity and sound velocity, are calculated. Thermal conductivity of the gaseous fuel mixture can be calculated by k= x.sub.ik.sub.i where k.sub.i is the thermal conductivity of pure individual gas and x.sub.i is the composition of the component in the mixture. The thermal conductivity data for each components at the specific T and P are taken as described in F. Uribe, E. Mason, J. Kestin, Thermal conductivity of nine polyatomic gases at low density, J. Phys. Chem. (1990), http://scitation.aip.org/content/aip/journal/jpcrd/19/5/10.1063/1.555864 (accessed Nov. 1, 2016), and B. A. Younglove, J. F. Ely, Thermophysical Properties of Fluids. II. Methane, Ethane, Propane, Isobutane, and Normal Butane, J. Phys. Chem. Ref. Data. 16 (1987) 577-798. doi:10.1063/1.555785. The thermal conductivity of the individual gas component in the mixture is calculated by equation 1 and equation 2. The equations are valid for methane when the reduced temperature above 1 and valid for other hydrocarbons at any temperature condition.

[00001]$\begin{array}{cc}{k}_i={10}^{-7}.Math.{\left(14.52.Math..Math.T_r-5.14\right)}^{2/3}.Math.\frac{C_p}{}& \left[\mathrm{equation}.Math..Math.1\right]\\ =T.Math..Math.c^{1/6}.Math.{M^{1/2}\left(\frac{101.325}{p_c}\right)}^{2/3}& \left[\mathrm{equation}.Math..Math.2\right]\end{array}$

where k.sub.i=vapor thermal conductivity of pure components, W/m K; Tr=reduced temperature, T/T.sub.c; T=temperature, K; T.sub.c=critical temperature, K; C.sub.p=heat capacity at constant pressure, J/kmol K; M=molecular weight and p.sub.c=critical pressure, kPa.

[0032] Sound velocity of the mixed gas is calculated by the following equation: Sound velocity,

[00002]$\begin{array}{cc}c=\sqrt{\frac{.Math.\mathrm{ZRT}}{M_g}}& \left[\mathrm{equation}.Math..Math.3\right]\end{array}$

[0033] where

[00003]$=\frac{C_p}{C_v};$

Z=compressibility factor; Constant pressure specific heat, C.sub.p=/(x.sub.iC.sub.pi(T.sub.i)) [J/mol.Math.K].

[0034] Temperature based constant pressure heat capacity for each component is calculated for ideal gas condition.

Constant volume specific heat, C.sub.y=C.sub.pR[J/mol.Math.K]

where the compressibility factor is calculated as a function of pseudo reduced pressure and temperature. Pseudo reduced pressure and temperature are calculated by using Suttons gas gravity method:

[00004]$\begin{array}{cc}{p}_{\mathrm{pc}}=7.56.Math.\mathrm{.8}-131.07.Math.-3.6.Math.2^{}& \left[\mathrm{equation}.Math..Math.4\right]\\ T_{\mathrm{pc}}=169.2+349.5.Math.-74.0.Math.2^{}& \left[\mathrm{equation}.Math..Math.5\right]\\ p_{\mathrm{pr}}=\frac{p}{p_{\mathrm{pc}}}& \left[\mathrm{equation}.Math..Math.6\right]\\ T_{\mathrm{pr}}=\frac{T}{T_{\mathrm{pc}}}& \left[\mathrm{equation}.Math..Math.7\right]\end{array}$

[0035] where relative density/Specific gravity

[00005]$\left(\right)=\frac{\mathrm{Molecular}.Math..Math.\mathrm{weight}.Math..Math.\mathrm{of}.Math..Math.\mathrm{the}.Math..Math.\mathrm{gas}.Math..Math.\left(M_g\right)}{\mathrm{Molecular}.Math..Math.\mathrm{weight}.Math..Math.\mathrm{of}.Math..Math.\mathrm{air}.Math..Math.\left(M_a\right)}$

[0036] Molecular weight of the gas mixture, M.sub.g=/M.sub.ix.sub.i where M.sub.i is the molecular weight of the component and x.sub.i is the composition of the component.

[0037] WI for the gas mixture is calculated by the following equation:

[00006]$\begin{array}{cc}\mathrm{WI}=\frac{H_c}{\sqrt{}}& \left[\mathrm{equation}.Math..Math.8\right]\end{array}$

where H.sub.c is the heating value of the gas mixture at the specific temperature and pressure and .sub.g is the relative density/specific gravity of the mixed gas.

[0038] An Aspen Plus process simulator is used for calculating the heating value of combustible components. A diagram of the apparatus developed by Aspen Plus is shown in FIG. 1. In the model, the combustible components CH.sub.4 and C.sub.2H.sub.6 are burned in presence of stoichiometric air ratio and the product leaves the reactor at the reaction temperature and pressure (at which temperature and pressure Wobbe index is calculated). CO.sub.2 also contributes to increase/decrease the Wobbe index of the mixed gas by changing the volume not by changing the total heat content.

[0039] The methane number of gas mixture is obtained from Cummins Westport's (CWI) web site, http://www.cumminswestport.com/fuel-quality-calculator (accessed Jul. 27, 2017). CWI uses the SAE based methane number calculation (SAE 922359 equation 4). In 2015, CWI's fuel quality calculator switched to a Cummins Proprietary methane number calculation. The Cummins Proprietary calculation provides a more accurate representation of the true MN of the fuel. Table 3 illustrates a portion of the data set including all the attributes for physical properties, WI, MN, and gas composition.

TABLE-US-00003 TABLE 3 Composition Thermal Sound Wobbe Methane Temperature, Pressure, CH.sub.4 C.sub.2H.sub.6 CO.sub.2 Conductivity, Velocity, Index Number T [K] P [kPa] % % % k 10.sup.3 [W/m .Math. K] v [m/s] [MJ/Nm.sup.3] [MN] 253.15 3447 94 6 0 27.36 350.86 54.9 84.5 253.15 3447 91.83 5.4 2.77 27.02 343.9 52.16 88.8 253.15 3447 89.66 4.8 5.54 26.68 336.48 49.53 93.2 253.15 3447 87.49 4.2 8.31 26.33 330.29 47.03 97.6 253.15 3447 85.32 3.6 11.08 25.99 324.43 44.62 102 253.15 3447 83.15 3 13.85 25.64 318.88 42.32 106.4 253.15 3447 80.98 2.4 16.62 25.3 313.6 40.1 110.7 253.15 3447 78.81 1.8 19.39 24.95 307.78 37.97 115 253.15 3447 76.64 1.2 22.16 24.61 303 35.91 119.2 253.15 3447 74.47 0.6 24.93 24.26 298.44 33.92 123.3 253.15 3447 72.3 0 27.7 23.92 294.08 32 127.3 253.15 3447 71.52 0 28.48 23.8 292.74 31.51 128.1 253.15 3447 73.69 0.6 25.71 24.14 297.04 33.41 124.2 253.15 3447 70.74 0 29.26 23.68 291.42 31.03 129 253.15 3447 75.86 1.2 22.94 24.49 301.54 35.38 120.1 253.15 3447 72.91 0.6 26.49 24.03 295.66 32.91 125 253.15 3447 69.96 0 30.04 23.56 290.11 30.55 129.8 253.15 3447 78.03 1.8 20.17 24.83 306.24 37.42 115.9 253.15 3447 75.08 1.2 23.72 24.37 300.09 34.86 121 253.15 3447 72.13 0.6 27.27 23.91 294.3 32.41 125.9 253.15 3447 69.18 0 30.82 23.45 288.83 30.07 130.7 253.15 3447 80.2 2.4 17.4 25.18 311.98 39.53 111.6

Constructing a Predictive Model

[0040] In a statistical modelling, regression analysis is commonly used for estimating the relationships among variables in a data set. Among many techniques for analysing a data set, a multiple regression approach is used to build a predictive model that estimates the WI. MN, and gas composition. Multiple regression allows for multiple dependent and independent variables, often called predictors, as well as multiple coefficients, .

[0041] The general multiple regression equation for a dependent variable, y.sub.i is

y.sub.i=.sub.0+.sub.1x.sub.1i+ . . . +.sub.mx.sub.mi+.sub.i for i=1 . . . n [equation 9]

where x.sub.mi represents independent variables; represents the coefficients that are normally unknown but to be found; represents the error; n represents the number of dependent values (that is equal to the samples); m represents the number of independent variables. Based on this, a multiple regression model is typically presented in the following form:

Y=X+[equation 10]

where Y is a vector representing values of dependant variables; X is a matrix representing values of independent variables. In most regression analysis, the value of is assumed to be randomly distributed. A regression model in general relates or approximates Y to a function of X and , that is, Y=f(X, ). In order to find a functional relational relationship between X and Y, an appropriate number of samples (or measurements) may be used as shown below, where n is number of samples and m is the number of variables. The data collection process for the sample data set used in this modelling is described above.

[00007]$Y=\left[\begin{array}{c}{y}_1\\ y_2\\ \mathrm{.Math.}\\ y_n\end{array}\right],X=\left[\begin{array}{ccccc}1& x_{11}& x_{21}& \mathrm{.Math.}& x_{k.Math..Math.1}\\ 1& x_{12}& x_{22}& \mathrm{.Math.}& x_{k.Math..Math.2}\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}& & \mathrm{.Math.}\\ 1& x_{1.Math..Math.n}& x_{2.Math.n}& \mathrm{.Math.}& x_{\mathrm{kn}}\end{array}\right],=\left[\begin{array}{c}{}_1\\ {}_2\\ \mathrm{.Math.}\\ {}_2\end{array}\right],\mathrm{.Math.}=\left[\begin{array}{c}{\mathrm{.Math.}}_1\\ {\mathrm{.Math.}}_2\\ \mathrm{.Math.}\\ {\mathrm{.Math.}}_2\end{array}\right]$

Y represents a vector consisting of WI, MN, and gas composition of all three components as dependent variables. X is a matrix representing the values of 2376 samples for all the physical properties as independent variables. The goal is to find a vector {tilde over ()} representing coefficients for the model that minimizes the square sum of the error, that is {tilde over ()}=argmin.sub..sub.i=1.sup.n(y.sub.i.sup.Tx.sub.i).sup.2. The coefficients (-vector) are calculated using the Least Squares Equation:

{tilde over ()}=(X.sup.TX).sup.1X.sup.TY [equation 11 ]

With the calculated coefficients {tilde over ()}, the estimated Y, is calculated by =X{tilde over ()}.

[0042] Further, the average error of the model is calculated by the root mean square (rms) between the actual value in the data set and the predicted value by the model as shown by the formula below:

[00008]$\begin{array}{cc}{\mathrm{error}}_{r.Math..Math.m.Math..Math.s}=\sqrt{\frac{{.Math.\mathrm{error}}^2}{2376}}& \left[\mathrm{equation}.Math..Math.12\right]\end{array}$

Results from the System

[0043] The system to predict the WI, MN, and gas composition from the physical properties is developed based on multiple regression using MATLAB.

Wobbe Index Prediction Model

[0044] The model is created base on four different variables: constant (1),

[00009]$\frac{1}{\mathrm{temperature}.Math..Math.\left(T\right)}.Math..Math.\left(x_1\right),.Math.\frac{1}{\mathrm{pressure}.Math..Math.\left(P\right)}.Math..Math.\left(x_2\right),.Math.\frac{1}{\mathrm{thermal}.Math..Math.\mathrm{conductivity}.Math..Math.\left(k\right)}.Math..Math.\left(x_3\right),\mathrm{and}$$\frac{1}{\mathrm{sound}.Math..Math.\mathrm{velocity}.Math..Math.\left(v\right)}.Math..Math.\left(x_4\right).$

Y is the matrix representing the WIs for all the samples in the data set. The coefficients for the model are:

[00010]$=\left[\begin{array}{c}70.76\\ -0.55\\ 0.000225\\ 3635.79\\ 0.05\end{array}\right].$

[0045] FIG. 2 shows the comparison of actual value and predicted value of WI. The trend line is desirable. While the accuracy of the predicted WIs in the range of 30-40 is very high, it is dropped slightly outside the range. The average error, error.sub.rms is 1.02 and the average percentage of error was 2.76%. Table 4 also shows the detailed regression statistics for the WI prediction model. These statistics strongly support the accuracy of the model for WI.

TABLE-US-00004 TABLE 4 R Square 0.98 Standard Error 0.94 Standard Error t Stat P-value Intercept (constant) 0.48 238.35 0.00 Variable 1 0.00 171.64 0.00 (Temperature, T) Variable 2 0.00 22.29 0.00 (Pressure, P) Variable 3 37.38 112.05 0.00 (Thermal Conductivity, k) Variable 4 0.00 24.07 0.00 (Sound Velocity, v)

[0046] The model for MN prediction is developed using four different variables: constant (1),

[00011]$\frac{1}{\mathrm{temperature}.Math..Math.\left(T\right)}.Math..Math.\left(x_1\right),.Math.\frac{1}{\mathrm{pressure}.Math..Math.\left(P\right)}.Math..Math.\left(x_2\right),.Math.\frac{1}{\mathrm{thermal}.Math..Math.\mathrm{conductivity}.Math..Math.\left(k\right)}.Math..Math.\left(x_3\right),\mathrm{and}$$\frac{1}{\mathrm{sound}.Math..Math.\mathrm{velocity}.Math..Math.\left(v\right)}.Math..Math.\left(x_4\right).$

Y is the matrix representing the methane number for all the samples. The coefficients for the model are:

[00012]$=\left[\begin{array}{c}51.71\\ 1.06\\ 0.000402\\ -7240.22\\ -0.10\end{array}\right].$

[0047] FIG. 3 shows the accuracy of the predicted values of methane number. As illustrated in the figure, the accuracy of the predicted MNs is much higher than the predicted WIs. The average error, error.sub.rms was 1.92 and the average percentage of error was 1.65%. Table 5 shows the detailed regression statistics for the MN prediction model.

TABLE-US-00005 TABLE 5 R Square 0.98 Standard Error 1.92 Standard Error t Stat P-value Intercept (constant) 0.98 52.67 0.00 Variable 1 0.01 147.54 0.00 (Temperature, T) Variable 2 0.00 22.66 0.00 (Pressure, P) Variable 3 76.20 95.01 0.00 (Thermal Conductivity, k) Variable 4 0.00 24.47 0.00 (Sound Velocity, v)

Gas Composition Prediction Model

[0048] The model for the gas composition prediction is developed using four different variables: constant (1),

[00013]$\frac{1}{\mathrm{temperature}.Math..Math.\left(T\right)}.Math..Math.\left(x_1\right),.Math.\frac{1}{\mathrm{pressure}.Math..Math.\left(P\right)}.Math..Math.\left(x_2\right),.Math.\frac{1}{\mathrm{thermal}.Math..Math.\mathrm{conductivity}.Math..Math.\left(k\right)}.Math..Math.\left(x_3\right),\mathrm{and}$$\frac{1}{\mathrm{sound}.Math..Math.\mathrm{velocity}.Math..Math.\left(v\right)}.Math..Math.\left(x_4\right).$

Table 6 shows the detailed regression statistics for the Gas Composition prediction model. The predicted data vs. the actual data plot for CH.sub.4, C.sub.2H.sub.6, and CO.sub.2 are presented in FIG. 4, FIG. 5, and FIG. 6, respectively.

TABLE-US-00006 TABLE 6 CH.sub.4 C.sub.2H.sub.6 CO.sub.2 R Square 0.98 0.93 0.98 Standard Error 0.94 0.43 1.16 Standard Error t Stat P-value CH.sub.4 C.sub.2H.sub.6 CO.sub.2 CH.sub.4 C.sub.2H.sub.6 CO.sub.2 CH.sub.4 C.sub.2H.sub.6 CO.sub.2 Intercept (constant) 0.48 0.22 0.59 238.35 43.39 40.99 0.00 0.00 0.00 Variable 1 (Temperature, T) 0.00 0.00 0.00 171.64 79.11 168.40 0.00 0.00 0.00 Variable 2 (Pressure, P) 0.00 0.00 0.00 22.29 14.13 23.30 0.00 0.00 0.00 Variable 3 (Thermal Cond., k) 37.38 17.16 46.15 112.05 50.21 109.40 0.00 0.00 0.00 Variable 4 (Sound Velocity, v) 0.00 0.00 0.00 24.07 15.26 25.17 0.00 0.00 0.00

[0049] The predicted methane composition shown in FIG. 4 shows negligible error in the range where methane composition is in the range of 70 to 85%. The predicted data out of the range has the largest fluctuation. The prediction of methane composition is very closely connected with the WI prediction since WI is mostly contributed by combustible component methane and ethane. Thus, as the ethane concentration is very small, the predicted methane composition and WI in the range is accurate. The errorm.sub.rms for methane composition prediction model is 0.94 and the average percentage of error was 1.22%.

[0050] FIG. 5 shows the predicted and actual ethane composition. As illustrated in the figure, with limited number of data points, the prediction of ethane composition is relatively inaccurate. The error.sub.rms for this prediction was 0.43 and the average percentage of error was 5.51%.

[0051] FIG. 6 shows the predicted CO.sub.2 composition vs actual CO.sub.2 composition. Particularly, the accuracy of the prediction is very high where the CO.sub.2 composition is above 10%. The error.sub.rms for CO.sub.2 composition prediction model was 1.16 and the average percentage of error was 21.59%.

Conclusion

[0052] The development in worldwide alternative fuel use is significantly dependent on the efficient combustion compatibility of the fuel in engine. WI and MN are the fuel quality indicator and the proper information on WI can make the combustion process efficient for the engine. In this research, in order to develop a predictive model that estimate WI, MN, and the composition of the gas, a data set is created that includes temperature, pressure, thermal conductivity, sound velocity, WI, MN, and composition of the gas mixture from fossil natural gas, anaerobic digester gas and landfill gas. The model takes thermal conductivity, sound velocity, temperature and pressure and predict WI, MN, and the gas composition. A database is developed for thermal conductivity, sound velocity, WI and MN of the gas mixture dependent on the temperature, pressure, and composition of the gas mixture. A model for prediction of WI, MN and composition of the components of fossil natural gas blended with anaerobic digester gas and landfill gas is developed.

[0053] The system can predict a WI and an MN of a gas mixture and also efficiently predicts the composition of methane, ethane and CO.sub.2 in the gas mixture. The prediction model, once coupled with a gas sensor, has the potential to make the combustion of alternative fuel more efficient.

[0054] One possible embodiment can be an engine controller that adjusts the air/fuel ratio for the specific input gas for optimum efficiency. For some controllers, the adjustment can be automated and in real-time. An onboard engine controller can comprise: a non-resident memory containing a database of Wobbe index values, methane number values, thermal conductivity values, and sound velocities of a gas mixture; a gas sensor to provide real time parameters of the fuel input stream; and a processor to run a predictive model coupled to the database and the sensor to determine the real-time parameters of the fuel input stream, process air flow input data, and provide real-time air/fuel ratio adjustment commands based on the predictive model; and a fuel control value to adjust the air/fuel ratio based on the adjustment command; whereby the air/fuel ratio is optimized for an engine based on the specific fuel provided. In some controller embodiments, the processor can further host a data set builder for the predictive model to build the database of Wobbe index values, methane number values, thermal conductivity values, and sound velocities of a gas mixture. In some controllers, the processor executes a linear regression predictive model.

[0055] The various embodiments described above are provided by way of illustration only and should not be construed to limit the invention. Those skilled in the art will readily recognize various modifications and changes that may be made to the claimed invention without following the example embodiments and applications illustrated and described herein, and without departing from the true spirit and scope of the claimed invention, which is set forth in the following claims.