Natural gas hydrates represents a huge source of energy. At the same time substantial leakages of natural gas from hydrates contributes significantly to climate changes. One of the most important reasons for these natural gas fluxes is leakage of seawater in to the hydrates from seafloor, through fracture systems. Hydrate dissociates if surrounding seawater is less than hydrate stability limit. Another interesting aspect of natural gas hydrates is the potential for safe CO2 storage. These different aspects of hydrates in natural sediments put demands on thermodynamic models. In addition to accurate description of pressure temperature hydrate stability there also a need to describe hydrate dissociation in concentration gradients towards surrounding water or surrounding gas as two examples. In this work we present new experimental data and an extensive thermodynamic model for hydrate. In contrast to conventional thermodynamic models for hydrate the model is consistent since all thermodynamic properties are derived from the Gibbs free energy. In this work we examine mixtures of CH4, C2H6, N2, CO2 from the China Sea and some synthetic mixtures, using this model. Maximum CO2 content in these mixtures are 60 mol% and the rest is dominated by CH4. Agreement between experimental data and model calculations are generally good and average deviations are below 5.5% for all the systems and conditions examined. Another aspect of the model is the ability for incorporation of effects of mineral surfaces. Specifically it is illustrated that adsorption of water on rust dominates liquid water drop out from gas as compared to water dew-point. Production of natural gas with such high CO2 content requires a strategy for CO2 separation and storage. It is proposed that the CH4 is separated from the C2H6, CO2 and N2 and cracked to H2 and CO2 using steam. Thermodynamic analysis indicates a significant potential for safe CO2 storage in natural gas hydrate and H2 as the only export product.