Field-Scale Modeling of Local Capillary Trapping During CO2 Injection into a Saline Aquifer

Local capillary trapping is the small-scale (10-2 to 10+1 m) CO2 trapping that is caused by the capillary pressure heterogeneity. The benefit of LCT, applied specially to CO2 sequestration, is that saturation of stored CO2 is larger than the residual gas, yet these CO2 are not susceptible to leakage through failed seals. Thus quantifying the extent of local capillary trapping is valuable in design and risk assessment of geologic storage projects. Modeling local capillary trapping is computationally expensive and may even be intractable using a conventional reservoir simulator. In this paper, we propose a novel method to model local capillary trapping by combining geologic criteria and connectivity analysis. The connectivity analysis originally developed for characterizing well-to-reservoir connectivity is adapted to this problem by means of a newly defined edge weight property between neighboring grid blocks, which accounts for the multiphase flow properties, injection rate, and gravity effect. Then the connectivity is estimated from shortest path algorithm to predict the CO2 migration behavior and plume shape during injection. A geologic criteria algorithm is developed to estimate the potential local capillary traps based only on the entry capillary pressure field. The latter is correlated to a geostatistical realization of permeability field. The extended connectivity analysis shows a good match of CO2 plume computed by the full-physics simulation. We then incorporate it into the geologic algorithm to quantify the amount of LCT structures identified within the entry capillary pressure field that can be filled during CO2 injection. Several simulations are conducted in the reservoirs with different level of heterogeneity (measured by the Dykstra-Parsons coefficient) under various injection scenarios. We find that there exists a threshold Dykstra-Parsons coefficient, below which low injection rate gives rise to more LCT; whereas higher injection rate increases LCT in heterogeneous reservoirs. Both the geologic algorithm and connectivity analysis are very fast; therefore, the integrated methodology can be used as a quick tool to estimate local capillary trapping. It can also be used as a potential complement to the full-physics simulation to evaluate safe storage capacity.