Estimating Plume Volume for Geologic Storage of CO 2 in Saline Aquifers Christine Doughty Earth Sciences Division Lawrence Berkeley National Laboratory July 2008 Typically, when a new subsurface flow and transport problem is first being considered, very simple models with a minimal number of parameters are used to get a rough idea of how the system will evolve. For a hydrogeologist considering the spreading of a contaminant plume in an aquifer, the aquifer thickness, porosity, and permeability might be enough to get started. If the plume is buoyant, aquifer dip comes into play. If regional groundwater flow is significant or there are nearby wells pumping, these features need to be included. Generally, the required parameters tend to be known from pre-existing studies, are parameters that people working in the field are familiar with, and represent features that are easy to explain to potential funding agencies, regulators, stakeholders, and the public. The situation for geologic storage of carbon dioxide (CO 2 ) in saline aquifers is quite different. It is certainly desirable to do preliminary modeling in advance of any field work since geologic storage of CO 2 is a novel concept that few people have much experience with or intuition about. But the parameters that control CO 2 plume behavior are a little more daunting to assemble and explain than those for a groundwater flow problem. Even the most basic question of how much volume a given mass of injected CO 2 will occupy in the subsurface is non-trivial. However, with a number of simplifying assumptions, some preliminary estimates can be made, as described below. To make efficient use of the subsurface storage volume available, CO 2 density should be large, which means choosing a storage formation at depths below about 800 m, where pressure and temperature conditions are above the critical point of CO 2 (P = 73.8 bars, T = 31 o C). Then CO 2 will exist primarily as a free-phase supercritical fluid, while some CO 2 will dissolve into the aqueous phase. A mass balance for CO 2 may be written as M = V + V where M is the total mass of CO 2 injected, φ is the porosity of the storage formation, S g is the saturation of free-phase CO 2 (that is, the fraction of pore space filled with free-phase CO 2 ), S l = 1 – S g is the saturation of the aqueous phase (water plus dissolved salt plus dissolved CO 2 ), ρ g and ρ l are densities of the CO 2 and aqueous phases, respectively, X l is the mass fraction of CO 2 dissolved in the aqueous phase, and V is plume volume. Angle brackets represent a spatial average over the plume. The first term of Equation (1) represents the mass of free-phase CO 2 and the second term the mass of dissolved CO 2 . Assuming that φ, S , ρ, and X l are not correlated enables the angle brackets to be dropped,
Read full abstract