James L. Palandri | USGS Science

USGS2012

Sequestration of non-pure carbon dioxide streams in iron oxyhydroxide-containing saline repositories

Iron oxyhydroxide, goethite (α-FeOOH), was evaluated as a potential formation mineral reactant for trapping CO 2 in a mineral phase such as siderite (FeCO 3 ), when a mixture of CO 2 -SO 2 flue gas is injected into a saline aquifer. Two thermodynamic simulations were conducted, equilibrating a CO 2 -SO 2 fluid mixture with a NaCl-brine and Fe-rich rocks at 150 °C and 300 bar. The modeling studies evaluated mineral and fluid composition at equilibrium and the influence of pH buffering in the system. Results show siderite precipitates both in the buffered and unbuffered system; however, the presence of an alkaline pH buffer enhances the stability of the carbonate. Based on the model, an experiment was designed to compare with thermodynamic predictions. A CO 2 -SO 2 gas mixture was reacted in 150 ml of NaCl-NaOH brine containing 10 g of goethite at 150 °C and 300 bar for 24 days. Mineralogical and brine chemistry confirmed siderite as the predominant reaction product in the system. Seventy-six mg of CO 2 are sequestered in siderite per 10 g of goethite.

International Journal of Greenhouse Gas Control↗

USGS2004

A compilation of rate parameters of water-mineral interaction kinetics for application to geochemical modeling

Geochemical reaction path modeling is useful for rapidly assessing the extent of water-aqueous-gas interactions both in natural systems and in industrial processes. Modeling of some systems, such as those at low temperature with relatively high hydrologic flow rates, or those perturbed by the subsurface injection of industrial waste such as CO2 or H2S, must account for the relatively slow kinetics of mineral-gas-water interactions. We have therefore compiled parameters conforming to a general Arrhenius-type rate equation, for over 70 minerals, including phases from all the major classes of silicates, most carbonates, and many other non-silicates. The compiled dissolution rate constants range from -0.21 log moles m-2 s-1 for halite, to -17.44 log moles m-2 s-1 for kyanite, for conditions far from equilibrium, at 25 ?C, and pH near neutral. These data have been added to a computer code that simulates an infinitely well-stirred batch reactor, allowing computation of mass transfer as a function of time. Actual equilibration rates are expected to be much slower than those predicted by the selected computer code, primarily because actual geochemical processes commonly involve flow through porous or fractured media, wherein the development of concentration gradients in the aqueous phase near mineral surfaces, which results in decreased absolute chemical affinity and slower reaction rates. Further differences between observed and computed reaction rates may occur because of variables beyond the scope of most geochemical simulators, such as variation in grain size, aquifer heterogeneity, preferred fluid flow paths, primary and secondary mineral coatings, and secondary minerals that may lead to decreased porosity and clogged pore throats.

Open-File Report↗