Natural Gas Quality, Energy & Wobbe Index

The standard reference framework for calculating calorific values, density, and Wobbe index is established by ISO 6976:2016 (International Organization for Standardization) and GPA 2172-09 / API MPMS Chapter 14.5. The calculation begins with the normalized molar composition of the hydrocarbon mixture.

Molar Mass & Ideal Specific Gravity

The mixture molecular weight (\(M_{mix}\)) is the mole-fraction-weighted sum of individual pure component molecular weights:

\[ M_{mix} = \sum_{j=1}^{N} x_j \cdot M_j \]

The ideal relative density (ideal specific gravity, \(G_{ideal}\)) is the ratio of mixture molar mass to that of standardized dry air (\(M_{air} = 28.9626 \text{ kg/kmol}\) per ISO 6976):

\[ G_{ideal} = \frac{M_{mix}}{M_{air}} \]

Real-Gas Compressibility Factor (\(Z\)) via Summation Factors

At metering conditions (e.g., 15°C, 101.325 kPa), natural gas deviates slightly from ideal gas behavior due to intermolecular van der Waals forces. Rather than relying on computationally heavy cubic equations of state, ISO 6976 mandates the rigorous summation factor method, derived from the second virial coefficient:

\[ \sqrt{1 - Z(T, P)} = \sum_{j=1}^{N} x_j \cdot s_j(T, P) \]

\[ Z(T, P) = 1 - \left( \sum_{j=1}^{N} x_j \cdot s_j(T, P) \right)^2 \]

Where \(s_j(T, P)\) is the dimensionless summation factor for pure component \(j\) at metering temperature \(T\) and pressure \(P\). For dry sweet natural gas at 15°C and atmospheric pressure, \(Z\) typically ranges from 0.9965 to 0.9982. The real relative density (\(G_{real}\)) is then computed by incorporating the compressibility of both the gas mixture and dry air (\(Z_{air} = 0.99958\)):

\[ G_{real} = G_{ideal} \cdot \frac{Z_{air}(T, P)}{Z_{mix}(T, P)} \]

Superior (Gross) vs. Inferior (Net) Heating Value

Combustion energy is calculated from standard molar enthalpies of combustion (\(\Delta H_c^\circ\)):

• Gross Heating Value (HHV / Superior): Assumes that all water produced during hydrocarbon combustion is completely condensed to liquid state at reference temperature, releasing its latent heat of vaporization (\(\Delta H_{vap} \approx 44.01 \text{ kJ/mol}\)).
• Net Heating Value (LHV / Inferior): Assumes that all water produced remains entirely in vapor phase in the flue exhaust, representing the net heat physically recoverable in non-condensing industrial boilers, furnaces, and gas turbines.

The real volumetric heating value is obtained by dividing the molar calorific value by the real molar volume \(V_m(T, P) = \frac{Z \cdot R \cdot T}{P}\):

\[ HHV_{real} = \frac{HHV_{ideal}}{Z_{mix}} = \frac{\sum x_j \cdot H_{c, j}^\circ}{V_{ideal} \cdot Z_{mix}} \quad \left[\frac{\text{MJ}}{\text{Sm}^3} \text{ or } \frac{\text{BTU}}{\text{scf}}\right] \]

A frequent and dangerous misconception in gas utility and power plant operations is assuming that two gas streams with the same Gross Heating Value (HHV) will burn identically in an industrial burner or aeroderivative gas turbine. In reality, fuel gas delivery through a fixed-geometry burner nozzle is controlled by Bernoulli's orifice flow equation:

\[ Q_{flow} = C_d \cdot A_{orifice} \cdot \sqrt{\frac{2 \Delta P}{\rho_{gas}}} \propto \frac{1}{\sqrt{G_{real}}} \]

Because gas flow rate is inversely proportional to the square root of relative density, the total thermal energy delivered through the burner orifice per unit time is the product of volumetric flow and volumetric heating value:

\[ \dot{Q}_{thermal} = Q_{flow} \cdot HHV \propto \frac{HHV}{\sqrt{G_{real}}} = \mathbf{W_G} \]

This critical ratio is the Wobbe Index (\(W_G\) for Gross, \(W_N\) for Net). If two distinct gas compositions possess the exact same Wobbe Index, the burner nozzle will deliver the identical heat release rate at identical manifold supply pressure, requiring zero mechanical adjustment to the gas train or control valves.

Combustion Failure Modes from Wobbe Variations

Consider a pipeline custody transfer metering skid receiving a typical natural gas stream at a standard volume flow of \(10,000 \text{ Sm}^3\text{/h}\) (measured at 15°C, 101.325 kPa):

Step-by-Step Numerical Walkthrough:

1. Mixture Molar Mass: \(M_{mix} = 0.90(16.043) + 0.05(30.070) + 0.02(44.097) + 0.015(44.010) + 0.015(28.013) = \mathbf{18.904 \text{ g/mol}}\).
2. Ideal Specific Gravity: \(G_{ideal} = \frac{18.904}{28.9626} = \mathbf{0.6527}\).
3. Summation Factor & Z: \(\sum x_j s_j = 0.90(0.0447) + 0.05(0.0922) + 0.02(0.1338) + 0.015(0.0755) + 0.015(0.0210) = 0.04896\).
   \(Z_{mix} = 1 - (0.04896)^2 = \mathbf{0.9976}\).
4. Real Specific Gravity: \(G_{real} = 0.6527 \cdot \frac{0.99958}{0.9976} = \mathbf{0.6540}\).
5. Gross Heating Value (Real): \(HHV_{ideal} = \frac{0.90(890.63) + 0.05(1560.69) + 0.02(2219.17)}{23.6409} = 39.083 \text{ MJ/m}^3\).
   \(HHV_{real} = \frac{39.083}{0.9976} = \mathbf{39.177 \text{ MJ/Sm}^3}\) (equivalent to \(\mathbf{1051.5 \text{ BTU/scf}}\)).
6. Gross Wobbe Index: \(W_G = \frac{39.177}{\sqrt{0.6540}} = \mathbf{48.44 \text{ MJ/Sm}^3}\) (\(\mathbf{1300.2 \text{ BTU/scf}}\)).
7. Daily Energy Flow Rate: \(\dot{E} = 10,000 \frac{\text{Sm}^3}{\text{h}} \cdot 39.177 \frac{\text{MJ}}{\text{Sm}^3} \cdot 24 \frac{\text{h}}{\text{d}} = 9,402,480 \text{ MJ/d} = \mathbf{9,402.5 \text{ GJ/d}}\) (\(\mathbf{8,911.8 \text{ MMBTU/day}}\)).
8. Continuous Thermal Power: \(P_{th} = \frac{9,402.5 \times 10^9 \text{ J}}{86,400 \text{ s}} = \mathbf{108.83 \text{ MW}_{th}}\).