Numerical simulation of CO2 leakage from a geologic disposal reservoir including transitions from super- to sub-critical conditions, and boiling of liquid of CO2

The critical point of CO{sub 2} is at temperature and pressure conditions of T{sub crit} = 31.04 C, P{sub crit} = 73.82 bar. At lower (subcritical) temperatures and/or pressures, CO{sub 2} can exist in two different phase states, a liquid and a gaseous state, as well as in two-phase mixtures of these states. Disposal of CO{sub 2} into brine formations would be made at supercritical pressures. However, CO{sub 2} escaping from the storage reservoir may migrate upwards towards regions with lower temperatures and pressures, where CO{sub 2} would be in subcritical conditions. An assessment of the fate of leaking CO{sub 2} requires a capability to model not only supercritical but also subcritical CO{sub 2}, as well as phase changes between liquid and gaseous CO{sub 2} in sub-critical conditions. We have developed a methodology for numerically simulating the behavior of water-CO{sub 2} mixtures in permeable media under conditions that may include liquid, gaseous, and supercritical CO{sub 2}. This has been applied to simulations of leakage from a deep storage reservoir in which a rising CO{sub 2} plume undergoes transitions from supercritical to subcritical conditions. We find strong cooling effects when liquid CO{sub 2} rises to elevations where it begins to boil and evolve a gaseous CO{sub 2} phase. A three-phase zone forms (aqueous - liquid - gas), which over time becomes several hundred meters thick as decreasing temperatures permit liquid CO{sub 2} to advance to shallower elevations. Fluid mobilities are reduced in the three-phase region from phase interference effects. This impedes CO{sub 2} upflow, causes the plume to spread out laterally, and gives rise to dispersed CO{sub 2} discharge at the land surface. Our simulation suggests that temperatures along a CO{sub 2} leakage path may decline to levels low enough so that solid water ice and CO{sub 2} hydrate phases may be formed