Natural Gas Hydrate Prevention & Inhibitor Sizing

Gas hydrates are crystalline, ice-like compounds formed when water and light natural gas molecules (methane, ethane, propane, and isobutane) contact one another at elevated pressures and moderate temperatures.

Because hydrates form at temperatures well above the freezing point of water (often 10°C to 20°C under typical transmission pipeline pressures), preventing hydrate crystallization is the single most critical flow assurance task in upstream and midstream operations.

1. Clathrate Crystal Structures

Hydrates are non-stoichiometric clathrate inclusion compounds. Water forms hydrogen-bonded host lattices consisting of three primary crystalline structures: * Structure I (sI): Small spherical cages that trap small molecules such as Methane ($CH_4$) and Ethane ($C_2H_6$). * Structure II (sII): Larger cages that trap Propane ($C_3H_8$) and Isobutane ($i-C_4H_{10}$). These typically exhibit the highest hydrate formation temperatures. * Structure H (sH): Hexagonal cages capable of trapping heavy hydrocarbons like methylcyclopentane when paired with methane.

2. Predicting the Hydrate Formation Curve

For sweet natural gas, the hydrate formation temperature ($T_{hyd}$) is calculated from pipeline pressure ($P$) and gas specific gravity ($\gamma_g$) using the GPSA / Motiee formulation:

$$T_{hyd} (^\circ C) = -238.29 + 78.997 \log_{10}(P_{kPa}) - 5.352 (\log_{10} P_{kPa})^2 + (34.904 + 0.5527 \log_{10} P_{kPa}) \gamma_g$$

Where: * $P_{kPa}$ is the absolute system pressure in kilopascals ($1 \text{ bar} = 100 \text{ kPa}$). * $\gamma_g$ is the gas specific gravity relative to air ($MW_{gas} / 28.964$).

3. Thermodynamic Inhibition (The Hammerschmidt Equation)

Thermodynamic inhibitors depress the chemical activity of liquid water, shifting the hydrate equilibrium curve toward lower temperatures and higher pressures.

The required concentration ($W$) of inhibitor in the free liquid water phase is governed by the Hammerschmidt equation:

$$\Delta T = \frac{K_H \cdot W}{M \cdot (100 - W)}$$

Inverting this equation allows engineers to directly calculate the required inhibitor concentration ($W$, wt%):

$$W = \frac{100 \cdot M \cdot \Delta T}{K_H + M \cdot \Delta T}$$

Where: * $\Delta T$ is the required hydrate temperature depression ($T_{hyd} - T_{operating} + \text{Safety Margin}$). * $M$ is the molecular weight of the inhibitor (Methanol = $32.04$, MEG = $62.07$, DEG = $106.12$). * $K_H$ is the empirical constant: * For Celsius: $K_H = 1297$ (Methanol), $1500$ (MEG), $2222$ (DEG). * For Fahrenheit: $K_H = 2335$ (Methanol), $2700$ (MEG), $4000$ (DEG).

4. Total Injection Rate Calculation

The total injection rate must satisfy two requirements: 1. Liquid Phase Dosage: $\dot{m}{liq} = \dot{m}$ 2. } \cdot \frac{W}{100 - WVaporization Loss (Methanol): Because methanol has significant vapor pressure, a fraction vaporizes directly into the natural gas stream, requiring an additional vapor loss allowance ($\approx 0.5 - 1.5 \text{ lb MeOH/MMSCF}$ per 100 psi of line pressure).