The invention relates to a chemo-metrical prediction of methane index for natural gas. Specifically, the invention relates to a method and system for predicting the Wobbe index and methane number of a renewable natural gas by the measurement of simple physical properties.
Renewable natural gas (RNG), i.e., natural gas produced from renewable feedstocks (e.g., landfill gas, anaerobic digestion gas, etc.) is an important alternative fuel that can contribute to achieving a number of goals set by local and federal government agencies relating to fossil fuel replacement and greenhouse gas (GHG) emissions reduction in the transportation sector. Most RNG production projects are small to medium scale by nature and comprehensive gas clean-up or upgrading to meet the fuel specifications of pipeline natural gas is often not feasible from a project economic perspective. This results in most RNG resources being wasted (e.g., flaring) or being left unused.
The RNG from a landfill site or anaerobic digester comes with a significant amount of CO2; the CH4 concentration varies from 35 to 70%, whereas CO2 composition varies from 15 to 45%. RNG upgrading to pipeline natural gas requires removal of CO2 from biogas as well as clean-up of impurities, drying, and compression. Among these series of processes, CO2 removal is the most expensive, and it can cost more than $2/mcf.
The typical calorific value of RNG without CO2 removal is around 50 to 60% of equal volume of fossil natural gas, and it varies significantly by project site and season.
Two of the most important parameters for fuel quality of natural gas are Wobbe Index (WI) and Methane Number (MN). WI is the ratio of the fuel calorific value to the square root of its specific gravity or relative density. WI is a critical factor in evaluating interchangeability of fuel gases such as natural gas and liquefied petroleum gas (LPG), and is used to compare the combustion energy output of different composition of fuel gases. Two fuels having identical WI have the identical energy output for a set of given operating condition. Therefore, WI is used in estimating the energy output in a wide variety of equipment and processes that involve natural gas (NG) combustion.
The MN of a gaseous fuel is defined as the methane composition (vol %) combined with hydrogen that makes the same knocking of the gas fuel under specified operating conditions in a knock testing engine. A different engine has a different MN range for suitable operation. MN is an important parameter in measuring engine performance, especially when the fuel is from renewable source with a high possibility of MN variation.
The WI is typically measured using bulky, complex and expensive analyzers. These devices measure the energy value of the fuel by direct calorimetry followed by a separate measurement of density by an optical method. The complex, destructive, and expensive nature of existing WI measurement systems prevent its use for on-site application. Traditionally MN is calculated from the gas composition. The gas composition can be measured by gas chromatography (GC) or a residual gas analyzer (RGA), which is time consuming and also requires a laboratory environment.
There is also a method and apparatus available for measuring the calorific value of a gas based on measuring the thermal conductivity value at two different temperatures and sound velocity values. An on-chip system that measures the energy content of fuel gases, such as natural gas, biogas and hydrogen, has also been proposed. A proposed thermal conductivity microsensor system uses a device that is able to detect the thermal conductivity of the gas mixture with insignificant errors, but the effect of temperature variation has not been studied. Another proposed an approach is to measure the MN and lower heating value of a gas online, which is valid for certain temperature and thermal conductivity of the gas. However, a variation of a wide range of temperature and pressure has not been studied. Studies for determining the gas mixture composition or energy content or thermal conductivity are available for either online or offline systems. Another proposed idea uses a thermal conductivity microsensor to measure the MN online, but no practical application or implementation of the proposed idea has been produced. Another study for designing and implementation of an analytical calculation model for determining the real-time composition of gas mixtures has been performed where the effects caused by the compressibility of the gas mixture and temperature are not taken into account in the model. The thermal conductivity (k) and sound velocity (v) through the gas is dependent on the temperature, pressure, density, and the composition of the gas mixture.
As discussed above, the previous research results show that the WI, MN, thermal conductivity, and sound velocity are dependent on the gas mixture. A limited scale of research also shows a possibility of estimating MN and gas composition using thermal conductivity and sound velocity instead of using the long and expensive process of using GC or RGA.
Given the shortcomings of the prior art, there is a need for a method and system for predicting the Wobbe Index and Methane Number of renewable natural gas by the measurement of simple physical properties. The present invention accomplishes this goal by using a predictive model that can accurately and efficiently estimate the WI, MN, and composition of the gas mixture based on only easily measurable physical properties, such as the thermal conductivity, sound velocity, temperature, and pressure of the gas to predict the quality and composition of gas in real-time.
According to one aspect of a preferred embodiment, the WI is predicted by using a single thermal conductivity value.
According to yet another aspect of a preferred embodiment, a new fuel quality sensor enables the use of RNG without removing CO2. Efficiency and the economy of the process improves significantly if CO2 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.
According to yet another aspect of a preferred embodiment, a method for fuel quality prediction enables the use of renewable natural gas (RNG) without removing CO2. The calorific value of RNG without CO2 removal is around 50 to 60% of fossil natural gas and varies significantly with RNG sources. Due to this fact, onsite measurement of the gas properties is preferable to ensure that expected quality is maintained, and better combustion efficiency is provided by adjusting the air/fuel ratio. Wobbe Index (WI) and Methane Number (MN) are the natural gas quality indicator and these parameters are dependent on the composition of gas mixture. This system and method builds a database for WI, MN, thermal conductivity and sound velocity through the gas mixture as a function of temperature, pressure and composition of the gas mixture. The system and method also develops a model for predicting the WI, MN and composition of the gas mixture. It takes temperature, pressure, thermal conductivity and sound velocity of the gas mixture and uses these variables to predict the WI, MN and composition. The prediction model can be based on MATLAB and may use the linear regression method for optimization.
In some embodiments, the system has been found to predict the WI and MN with an average prediction error of 1.386, and 0.559, respectively. In some embodiments, the system also predicts the composition of methane, ethane and CO2 with an average error of 0.590, 0.279, and 0.509, respectively. The prediction model of the system, combined with a micro-sensor, has the potential to make the combustion of alternative fuel more efficient and ensure the quality of gaseous fuel supplied to an engine.
According to yet another aspect of a preferred embodiment, a database for WI, MN, thermal conductivity, sound velocity of a gas mixture is built; there is a prediction of WI and MN for a gaseous fuel mixture; a prediction of a compositions of the gaseous fuel mixture; and there is a prediction model coupled to a sensor that enables the combustion of gas efficiently.
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:
For illustrative purposes, the present invention is embodied in the apparatus and method generally shown and described herein with reference to
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.
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 CH4, CO2 and C2H6 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.
For example, a mixture of 40% fossil, 30% anaerobic digester gas and 30% landfill gas contains 78.64% CH4, 18.96% CO2, and 2.4% C2H6. 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.
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.
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=Σ xiki where ki is the thermal conductivity of pure individual gas and xi 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.
where ki=vapor thermal conductivity of pure components, W/m K; Tr=reduced temperature, T/Tc; T=temperature, K; Tc=critical temperature, K; Cp=heat capacity at constant pressure, J/kmol K; M=molecular weight and pc=critical pressure, kPa.
Sound velocity of the mixed gas is calculated by the following equation: Sound velocity,
where
Z=compressibility factor; Constant pressure specific heat, Cp=Σ/(xiCpi(Ti)) [J/mol·K].
Temperature based constant pressure heat capacity for each component is calculated for ideal gas condition.
Constant volume specific heat, Cy=Cp−R[J/mol·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:
where relative density/Specific gravity
Molecular weight of the gas mixture, Mg=Σ/Mixi where Mi is the molecular weight of the component and xi is the composition of the component.
WI for the gas mixture is calculated by the following equation:
where Hc is the heating value of the gas mixture at the specific temperature and pressure and g is the relative density/specific gravity of the mixed gas.
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
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.
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, β.
The general multiple regression equation for a dependent variable, yi is
y
i=β0+β1x1i+ . . . +βmxmi+εi for i=1 . . . n [equation 9]
where xmi 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.
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βΣi=1n(yi−βTxi)2. The coefficients (β-vector) are calculated using the Least Squares Equation:
{tilde over (β)}=(XTX)−1XTY [equation 11 ]
With the calculated coefficients {tilde over (β)}, the estimated Y, Ŷ is calculated by Ŷ=X{tilde over (β)}.
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:
The system to predict the WI, MN, and gas composition from the physical properties is developed based on multiple regression using MATLAB.
The model is created base on four different variables: constant (1),
Y is the matrix representing the WIs for all the samples in the data set. The coefficients for the model are:
The model for MN prediction is developed using four different variables: constant (1),
Y is the matrix representing the methane number for all the samples. The coefficients for the model are:
The model for the gas composition prediction is developed using four different variables: constant (1),
Table 6 shows the detailed regression statistics for the Gas Composition prediction model. The predicted data vs. the actual data plot for CH4, C2H6, and CO2 are presented in
The predicted methane composition shown in
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.
The system can predict a WI and an MN of a gas mixture and also efficiently predicts the composition of methane, ethane and CO2 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.
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.
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.
This Application claims the benefit of U.S. Provisional Patent Application No. 62/632,839, titled “Chemo-Metrical Prediction of Methane Index for the Natural Gas,” filed Feb. 20, 2018, the contents of which are incorporated herein by reference in their entirety.
This invention was made with Government support under Grant No. UCR-DOT-306 awarded by U.S. Department of Transportation. The government has certain rights in this invention.
Number | Date | Country | |
---|---|---|---|
62632839 | Feb 2018 | US |