Buoyancy-Dominated Multiphase Flow and Its Effect on Geological Sequestration of CO2
- Steven L. Bryant (U. of Texas at Austin) | Srivatsan Lakshminarasimhan (U. of Texas at Austin) | Gary A. Pope (U. of Texas at Austin)
- Document ID
- Society of Petroleum Engineers
- SPE Journal
- Publication Date
- December 2008
- Document Type
- Journal Paper
- 447 - 454
- 2008. Society of Petroleum Engineers
- 5.5 Reservoir Simulation, 1.8.5 Phase Trapping, 5.1.4 Petrology, 5.3.4 Integration of geomechanics in models, 5.3.2 Multiphase Flow, 4.3.4 Scale, 5.4.2 Gas Injection Methods, 4.1.2 Separation and Treating, 2.4.3 Sand/Solids Control, 5.1 Reservoir Characterisation, 5.1.5 Geologic Modeling, 5.4 Enhanced Recovery, 5.1.2 Faults and Fracture Characterisation, 1.2.3 Rock properties, 6.5.7 Climate Change, 5.1.1 Exploration, Development, Structural Geology, 6.5.3 Waste Management
- 2 in the last 30 days
- 1,586 since 2007
- Show more detail
- View rights & permissions
|SPE Member Price:||USD 5.00|
|SPE Non-Member Price:||USD 35.00|
We have previously proposed the "inject low and let rise" strategy of storing CO2 in deep saline aquifers. The idea is to maximize the amount of CO2 stored in immobile forms by letting CO2 rise toward the top seal of the aquifer but not reach it. The distance that the CO2 rises depends on the uniformity of the displacement front. In this paper, we address the question of whether the intrinsic instability of a buoyancy-driven immiscible displacement leads to fingering. Fingers could reach the top seal of the aquifer, leading to an accumulation of CO2 at large saturations. We study the mechanisms governing this type of displacement in a series of fine-grid numerical simulations. Each simulation begins with a finite volume of CO2 placed at large saturation at the bottom of a 2D aquifer. Only buoyancy forces drive the displacement. Boundaries are closed, so CO2 rises and brine falls as the simulation proceeds. Several fine-scale geostatistical realizations of permeability are considered, and the effects of capillary pressure, anisotropy, and dip angle are examined. In these simulations, buoyant instability has very little effect on the uniformity of the displacement front. Instead, the CO2 rises along preferential flow paths that are the consequence of spatially heterogeneous rock properties (permeability, drainage capillary pressure curve, and anisotropy). Capillary pressure broadens the lateral extent of the flow paths. If the formation beds are not horizontal, capillary pressure and anisotropy can cause the CO2 to move predominantly along the bedding plane rather than vertically. Accurate assessment of CO2 migration after injection ends will therefore require accurate characterization of the spatial correlation of permeability in the target formation and of the capillary pressure and relative permeability curves.
Storing CO2 in deep saline aquifers will be a key technology if society elects to limit the amount of greenhouse gases entering the atmosphere. Large-scale (106 tonnes of CO2 per year) examples of this type of storage are underway at Sleipner and In Salah, and more are planned (IPCC 2008). Effective mitigation of CO2 emissions will require many more projects of this type, storing on the order of 109 tonnes per year (Pacala and Socolow 2004). In terms of volumetric flow rates through wellbores, this rate of storage is of the same magnitude as the current global rate of oil production. Thus inexpensive, reliable methods of ensuring that stored CO2 remains in place will be essential.
CO2 can be stored in an aquifer in four modes: as a bulk phase within a structural trap, as a residual phase trapped by capillary forces, as aqueous species dissolved in brine, and as a precipitated mineral. The latter three forms of storage are permanent in the sense that the CO2 will remain in the aquifer at least as long as the residence time of water in the aquifer. On the other hand, CO2 held in a structural trap at large saturations (above residual) is potentially mobile in that it will remain trapped only as long as the seal remains intact. Storage methods that reduce the amount of potentially mobile CO2 correspondingly reduce the risk of leakage over the long term.
The inject-low-and-let-rise strategy is one such method (Kumar et al. 2005; Ozah et al. 2005). Under typical storage conditions, CO2 is less dense than brine. If CO2 is injected only into the lower part of an aquifer, then, after injection ends, the CO2 will continue to migrate, driven only by buoyancy. As CO2 rises into the upper part of the aquifer, it will leave behind a residual phase trapped by capillary forces. The permanency of residual phase trapping is the main motivation for this approach, but an additional benefit is that vertical movement toward the top seal is also retarded. By choosing the volume injected, one can, in principle, prevent the CO2 from reaching the top of the aquifer. The distance that the CO2 rises depends on the uniformity of the displacement front and on the saturation of CO2 behind the front. In this paper, we discuss factors that control the former feature. We will report on the latter in future publications.
Coarse-grid simulations suggest that CO2 will rise in a compact plume having a smooth outline. As the grid is refined, the shape of plume becomes more uneven. Can this loss of uniformity be attributed to the intrinsically unstable character of buoyancy-driven immiscible flow? In analogy with immiscible displacements that exhibit viscous instability, we might anticipate the emergence of fingers as the CO2 rises. Such fingers conceivably could reach the top seal of the aquifer quickly, even when the volume of stored CO2 is insufficient to allow a uniform displacement to reach the top. This could lead to an accumulation of potentially mobile CO2, the very situation the inject-low-and-let-rise strategy seeks to avoid. Thus, it is important to assess the extent to which gravity fingers develop under typical storage conditions for a range of target formations.
Some aspects of this problem are familiar from the long experience of gas and CO2 injection into oil reservoirs (Stalkup 1983). In gas-injection processes, the competition between viscous forces and buoyancy leads to gravity override. The larger mobility of the gas phase also leads to viscous fingering. We will see that some factors that govern gas-injection displacements also influence the situation of interest here--that is, when injection has ended and the only driving force is buoyancy. On one hand, this is not surprising. On the other, it should not be taken for granted because there has been relatively little examination of the buoyancy-dominated dynamics. The key question is whether the absence of competing forces allows the intrinsic instability of a buoyant displacement to dominate the shape of the plume.
The idealized initial condition for our simulations is an approximation of the situation commonly observed at the end of the injection period in simulations of the inject-low-and-let-rise strategy. The simplification allows us to attribute differences in behaviors unequivocally to differences in petrophysical properties and to the physics of buoyant flow. The understanding thus obtained will provide insight into the post-injection behavior when the injection period is simulated more realistically.
|File Size||7 MB||Number of Pages||8|
Bachu, S., Gunter, W.D., and Perkins, E.H. 1994. Aquifer disposal ofCO2: hydrodynamic and mineral trapping. Energy Conversion andManagement 35 (4): 269-279. doi:10.1016/0196-8904(94)90060-4.
Chang, Y.-B., Lim, M.T., Pope, G.A., and Sepehrnoori, K. 1994. CO2 Flow Patterns Under MultiphaseFlow: Heterogeneous Field-Scale Conditions. SPERE 9 (3):208-216. SPE-22654-PA. doi: 10.2118/22654-PA.
Ennis-King, J. and Paterson, L. 2005. Role of Convective Mixing in theLong-Term Storage of Carbon Dioxide in Deep Saline Formations. SPEJ10 (3): 349-356. SPE-84344-PA. doi: 10.2118/84344-PA.
Gunter, W.D., Wiwchar, B., and Perkins, E.H. 1997. Aquifer disposal ofCO2-rich greenhouse gases: Extension of the time scale of experimentfor CO2 -sequestering reactions by geochemical modeling.Mineralogy and Petrology 59 (1-2): 121-140.doi:10.1007/BF01163065.
Holtz, H.M. 2002. Residual GasSaturation to Aquifer Influx: A Calculation Method for 3-D Computer ReservoirModel Construction. Paper SPE 75502 presented at the SPE Gas TechnologySymposium, Calgary, 30 April-2 May. doi: 10.2118/75502-MS.
House, N.J., Faulder, D.D., Olson, G.L., and Fanchi, J.R. 2003. Simulation Study of CO2 Sequestrationin a North Sea Formation. Paper SPE 81202 presented at the SPE/EPA/DOEExploration and Production Environmental Conference, San Antonio, Texas, USA,10-12 March. doi: 10.2118/81202-MS.
Hovorka, S.D., Romero, M.L., Trevino, R.H., Warne, A.G., Ambrose, W.A.,Knox, P.R., and Tremblay, T.A. 2001. Phase II: Optimal Geological Environmentsfor Carbon Dioxide Disposal in Brine-Bearing Formations (Aquifers) in theUnited States. Final report, NETL Contract No. DE-AC26-98FT40417, USDOE/University of Texas at Austin (BEG), Austin, Texas (July 2001).
IPCC Special Report on Carbon Dioxide Capture and Storage, Chap. 5. http://www.ipcc.ch/pdf/special-reports/srccs/srccs_wholereport.pdf.Downloaded 14 October 2008.
Kumar, A., Ozah, R., Noh, M., Pope, G.A., Bryant, S., Sepehrnoori, K., andLake, L.W. 2005. ReservoirSimulation of CO2 Storage in Deep Saline Aquifers. SPEJ 10(3): 336-348. SPE-89343-PA. doi: 10.2118/89343-PA.
Nghiem, L., Sammon, P., Grabenstetter, J., and Ohkuma, H. 2004. Modeling CO2 Storage in Aquifers Witha Fully-Coupled Geochemical EOS Compositional Simulator. Paper SPE 89474presented at the SPE/DOE Symposium on Improved Oil Recovery, Tulsa, 17-21April. doi: 10.2118/89474-MS.
Noh, M., Lake, L.W., Bryant, S.L., and Araque-Martinez, A. 2007. Implicationsof Coupling Fractional Flow and Geochemistry for CO2 Injection in Aquifers.SPEREE 10 (4): 406-414. SPE-89341-PA. doi: 10.2118/89341-PA.
Ozah, R., Lakshminarasimhan, S., Pope, G.A., Sepehrnoori, K., and Bryant,S.L. 2005. Numerical Simulation ofthe Storage of Pure CO2 and CO2-H2S Gas Mixtures in Deep Saline Aquifers.Paper SPE 97255 presented at the SPE Annual Technical Conference andExhibition, Dallas, 9-12 October. doi: 10.2118/97255-MS.
Pacala, S. and Socolow, R. 2004. Stabilization wedges: Solvingthe climate problem for the next 50 years with current technologies.Science 305 (5686): 968-972. doi:10.1126/science.1100103.
Pruess, K., Xu, T., Apps, J., and Garcia, J. 2003. Numerical Modeling of AquiferDisposal of CO2. SPEJ 8 (1): 49-60. SPE-83695-PA. doi:10.2118/83695-PA
Stalkup, F.I. 1983. Miscible Displacement. Monograph Series, SPE,Richardson, Texas 8: 31-70.