Description: Highlights • CO2-methane exchange in a pressure vessel was simulated. • The model uses a detailed description of the kinetics for the CO2-methane exchange and simplifies the transport phenomena. • Irreversible dissociation rate of CH4- and CO2-hydrate in the pressure vessel was estimated as 0.02 and 0.03 mol m−3.s−1. • Formation of CO2-hydrate not only improved the quality of CO2 retention but also enhanced the methane recovery. Carbon dioxide exchange with methane in the clathrate structure has been shown beneficial in laboratory experiments and has been suggested as a field-scale technique for production of natural gas from gas-hydrate bearing sediments. Furthermore, the method is environmentally attractive due to the formation of CO2-hydrate in the sediments, leading to the geosequestration of carbon dioxide. However, the knowledge is still limited on the impact of small-scale heterogeneities on hydrate dissociation kinetics. In the present study, we developed a model for simulating laboratory experiments of carbon dioxide injection into a pressure vessel containing a mixture of gas hydrate and quartz sand. Four experiments at different temperature and pressure conditions were modeled. The model assumes that the contents are ideally mixed and aims to estimate the effective dissociation rate of gas hydrate by matching the model results with the experimental observations. Simulation results indicate that with a marginal offset the model was able to simulate different hydrate dissociation experiments, in particular, those that are performed at high pressures and low temperatures. At low pressures and high temperatures large discrepancies were noticed between the model results and the experimental observations. The mismatches were attributed to the development of extremely heterogeneous flow patterns at pore-scale, where field-scale models usually assume the characteristics to be uniform. Through this modeling study we estimated the irreversible dissociation rate of methane- and CO2-hydrate as 0.02 and 0.03 mol m(-3) s(-1), respectively.

Description: We present a transport-reaction model (TRACTION) specifically designed to account for non-ideal transport effects in the presence of thermodynamic (e.g. salinity or temperature) gradients. The model relies on the most fundamental concept of solute diffusion, which states that the chemical potential gradient (Maxwell’s model) rather than the concentration gradient (Fick’s law) is the driving force for diffusion. In turn, this requires accounting for species interactions by applying Pitzer’s method to derive species chemical potentials and Onsager coefficients instead of using the classical diffusion coefficients. Electrical imbalances arising from varying diffusive fluxes in multicomponent systems, like seawater, are avoided by applying an electrostatic gradient as an additional transport contribution. We apply the model to pore water data derived from the seawater mixing zone at the submarine Mercator mud volcano in the Gulf of Cadiz. Two features are particularly striking at this site: (i) Ascending halite-saturated fluids create strong salinity (NaCl) gradients in the seawater mixing zone that result in marked chemical activity, and thus chemical potential gradients. The model predicts strong transport-driven deviations from the mixing profile derived from the commonly used Fick’s diffusion model, and is capable of matching well with the profile shapes observed in the pore water concentration data. Even better agreement to the observed data is achieved when ion pairs are transported separately. (ii) The formation of authigenic gypsum (several wt%) occurs in the surface sediments, which is typically restricted to evaporitic surface processes. Very little is known about the gypsum paragenesis in the subseafloor and we first present possible controls on gypsum solubility, such as pressure, temperature, and salinity (pTS), as well as the common ion and ion pairing effects. Due to leaching of deep diapiric salt, rising fluids of the MMV are saturated with respect to gypsum (as well as celestite and barite). Several processes that could drive these fluids towards gypsum supersaturation and hence precipitation were postulated and numerically quantified. In line with the varied morphology of the observed gypsum crystals, gypsum paragenesis at the MMV is likely a combination of two temperature-related processes. Gypsum solubility increases with increasing temperature, especially in strong electrolyte solutions and the first mechanism involves the cooling of saturated fluids along the geothermal gradient during their ascent. Secondly, local temperature changes, i.e. cooling during the transition from MMV activity towards dormancy results in the cyclic build-up of gypsum. The model showed that the interpretation of field data can be majorly misguided when ignoring non-ideal effects in extreme diagenetic settings. While at first glance the pore water profiles at the Mercator mud volcano would indicate strong reactive influences in the seawater mixing zone, our model shows that the observed species distributions are in fact primarily transport-controlled. The model results for SO4 are particularly intriguing, as SO4 is shown to diffuse into the sediment along its increasing (!) concentration gradient. Also, a pronounced gypsum saturation peak can be observed in the seawater mixing zone. This peak is not related to the dissolution of gypsum but is simply a result of the non-ideal transport forces acting on the activity profile of SO4 and Ca profiles.

Description: Highlights: • Clay dehydration water expelled from buried sediments drives mud volcanism. • Rise of fluids mediated by crustal-scale strike-slip faults cross-cutting wedge. • On active accretionary wedge, petroleum accumulations were dismantled in Neogene. • 4He enrichment and δ13C-CH4 ~−50‰ in fluids reflect an open hydrocarbon system. • Petroleum pools remain on shallow margin. Microbial gas vented out of active wedge. Abstract: A geochemical study of the composition of hydrocarbon gases and helium isotopes (3He/4He) in fluids from Mud Volcanoes (MVs) located on and out of the active accretionary wedge of the Gulf of Cadiz (GoC) provides information on fluid sources and migrations in deeply buried sediments. The GoC is a tectonically active segment of the Africa-Iberia plate boundary occluded beneath the thick sediments of an accretionary wedge dissected by crustal-scale strike-slip faults. Initially built during the Miocene Gibraltar Arc subduction, the wedge has since developed toward the W-NW in an oblique convergent setting. Interstitial water expelled from clays undergoing diagenesis in buried sediments drives mud volcanism on the wedge, with MVs located along strike-slip faults mediating fluid ascent. The large excess of radiogenic helium (4He) in all GoC fluids agrees with a clay mineral dehydration source of water. Hydrocarbon gases from all deepwater MVs bear methane having similar stable carbon isotope compositions of ~−50‰VPDB whether fluids are highly enriched in methane relative to heavier homologues (C2+) or not (Methane / (Ethane + Propane) ~10 to 10,000). We suggest that methane with −50‰VPDB was largely diffused out of early generating source rocks, and became dissolved in the water expelled by the buried sediments. Consistently, low 3He/4He ratios suggest an open hydrocarbon system: Petroleum accumulations and 3He dissolved in the original sedimentary pore water have mostly escaped into the water column during the major Late Neogene compressional events. At present, some MVs vent CH4-rich fluids from dewatering sediments, while other structures located on active thrusts additionally vent C2+-rich gases generated by active Cretaceous source intervals. By contrast, evaporitic seals preserved petroleum accumulations on the shallow Moroccan Margin, while the westernmost MVs located out of the accretionary wedge vent microbial gas.

Description: Carbon dioxide (CO2) capture and storage (CCS) has been discussed as a potentially significant mitigation option for the ongoing climate warming. Natural CO2 release sites serve as natural laboratories to study subsea CO2 leakage in order to identify suitable analytical methods and numerical models to develop best-practice procedures for the monitoring of subseabed storage sites. We present a new model of bubble (plume) dynamics, advection-dispersion of dissolved CO2, and carbonate chemistry. The focus is on a medium-sized CO2 release from 294 identified small point sources around Panarea Island (South-East Tyrrhenian Sea, Aeolian Islands, Italy) in water depths of about 40–50 m. This study evaluates how multiple CO2 seep sites generate a temporally variable plume of dissolved CO2. The model also allows the overall flow rate of CO2 to be estimated based on field measurements of pH. Simulations indicate a release of ∼6900 t y–1 of CO2 for the investigated area and highlight an important role of seeps located at 〉20 m water depth in the carbon budget of the Panarea offshore gas release system. This new transport-reaction model provides a framework for understanding potential future leaks from CO2 storage sites.

Description: Highlights • CO2 gas bubbles are completely dissolved within 2 m above the seabed. • CO2 is not emitted into the atmosphere but retained in the North Sea. • Dissolved CO2 is rapidly dispersed by tidal currents in the North Sea. • Harmful effects on benthic biota occur in the direct vicinity of the leak. • Monitoring has to be performed at the seabed and close to the leak. Abstract Existing wells pose a risk for the loss of carbon dioxide (CO2) from storage sites, which might compromise the suitability of carbon dioxide removal (CDR) and carbon capture and storage (CCS) technologies as climate change mitigation options. Here, we show results of a controlled CO2 release experiment at the Sleipner CO2 storage site and numerical simulations that evaluate the detectability and environmental consequences of a well leaking CO2 into the Central North Sea (CNS). Our field measurements and numerical results demonstrate that the detectability and impact of a leakage of 〈55 t yr−1 of CO2 would be limited to bottom waters and a small area around the leak, due to rapid CO2 bubble dissolution in seawater within the lower 2 m of the water column and quick dispersion of the dissolved CO2 plume by strong tidal currents. As such, the consequences of a single well leaking CO2 are found to be insignificant in terms of storage performance. Only prolonged leakage along numerous wells might compromise long-term CO2 storage and may adversely affect the local marine ecosystem. Since many abandoned wells leak natural gas into the marine environment, hydrocarbon provinces with a high density of wells may not always be the most suitable areas for CO2 storage.