CH4 hydrate was modeled by considering and as hydrate formers.
Considered competing hydrate phase transitions via Gibb's free energy minimization.
Calculated Gibb's free energy using a non-equilibrium thermodynamic package.
Hydrate dissociation and reformation were modeled based on thermodynamics variables.
Modified code was tested by the depressurization method on a hypothetical model.