Lawrence Berkeley National Laboratory

Numerical modeling of Geologic Carbon Sequestration in permeable reservoirs initially containing hydrocarbons is conducted using the multi-phase, multi-component thermohydrologic simulator TOGA (TOUGH Oil, Gas, Aqueous; TOUGH stands for Transport Of Unsaturated Groundwater and Heat), to determine how phase and composition of the original fluids influence the extent of the zone where upward fugitive fluid migration could potentially occur, denoted Rf . The area within Rf comprises regions of substantially elevated pressure and free-phase CO2 saturation, where a breach in reservoir sealing capacity would lead to upward fugitive fluid migration. The model examines the conditions within the storage reservoir that could lead to fugitive flow, but does not model the fugitive flow itself. A one-dimensional radial model of the storage reservoir is used, and three initial phase conditions are considered: single-phase aqueous, two-phase gas-aqueous, and three-phase oil-gas-aqueous. Components that may be present are H2O, CO2, CH4, C4H10, and C10H22. The most important factors controlling Rf are (1) the initial gas-phase saturation within the reservoir, and (2) the lateral extent of multi-phase initial conditions, particularly CO2. The composition of liquid and gas phases has a secondary effect. The impact of reservoir depth, thickness, injection rate, and hydrologic properties are also briefly examined, with thickness (or equivalently injection rate) having the biggest effect. These results can help to understand important trends in potential response of CO2-EOR fields being considered for dedicated CO2 storage.