A robust and highly efficient phase boundary method for determining the thermodynamic equilibrium conditions of bulk gas hydrate systems