Although deep saline aquifers are found in all sedimentary basins and provide very large storage capacities, a little is known about them because they are rarely a target for the exploration. Furthermore, nearly all the experiments and simulations made for CO2 sequestration in deep saline aquifers are related to the sandstone formations. The aim of this study is to create a predictive model to estimate the CO2 storage capacity of the deep saline carbonate aquifers since a little is known about them. To create a predictive model, the variables which affect the CO2 storage capacity and their ranges are determined from published literature data. They are rock properties (porosity, permeability, vertical to horizontal permeability ratio), fluid properties (irreducible water saturation, gas permeability end point, Corey water and gas coefficients), reaction properties (forward and backward reaction rates) and reservoir properties (depth, pressure gradient, temperature gradient, formation dip angle, salinity), diffusion coefficient and Kozeny-Carman Coefficient. Other parameters such as pore volume compressibility and density of brine are calculated from correlations found in literature. To cover all possibilities, Latin Hypercube Space Filling Design is used to construct 100 simulation cases and CMG STARS is used for simulation runs. By using least squares method, a linear correlation is found to calculate CO2 storage capacity of the deep saline carbonate aquifers with a correlation coefficient 0.81 by using variables found from literature and simulation results. Numerical dispersion effects have been considered by increasing the grid dimensions. It has been found that correlation coefficient decreased to 0.77 when the grid size was increased from 250 ft to 750 ft. The sensitivity analysis shows that the most important parameter that affects CO2 storage capacity is depth since the pressure difference between formation pressure and fracture pressure increases with depth. Also, CO2 storage mechanisms are investigated at the end of 300 years of simulation. Most of the gas (up to 90%) injected into formation dissolves into the formation water and negligible amount of CO2 reacts with carbonate. This result is consistent with sensitivity analysis results since the variables affecting the solubility of CO2 in brine have greater affect on storage capacity of aquifers. Dimensionless linear and nonlinear predictive models are constructed to estimate the CO2 storage capacity of all deep saline carbonate aquifers and it is found that the best dimensionless predictive model is linear one independent of bulk volume of the aquifer.

Geological sequestration of CO₂ in deep saline reservoirs is one of the ways to reduce its continuous emission into the atmosphere to mitigate the greenhouse effect. The effectiveness of any CO₂ sequestration operation depends on pore volume and the sequestration efficiency of the reservoir. Sequestration efficiency is defined here as the maximum storage with minimum risk of leakage to the overlying formations or to the surface. This can be characterized using three risk parameters i) the time the plume takes to reach the top seal; ii) maximum lateral extent of the plume and iii) the percentage of mobile CO₂ present at any time. The selection among prospective saline reservoirs can be expedited by developing some semi-analytical correlations for these risk parameters which can be used in place of reservoir simulation study for each and every saline reservoir. Such correlations can reduce the cost and time for commissioning a geological site for CO₂ sequestration. To develop such correlations, a database has been created from a large number of compositional reservoir simulations for different elementary reservoir parameters including porosity, permeability, permeability anisotropy, reservoir depth, thickness, dip, perforation interval and constant pressure far boundary condition. This database is used to formulate different correlations that relate the sequestration efficiency to reservoir properties and operating conditions. The various elementary reservoir parameters are grouped together to generate different variants of gravity number used in the correlations. We update a previously reported correlation for time to hit the top seal and develop new correlations for other two parameters using the newly created database. A correlation for percentage of trapped CO₂ is also developed using a previously created similar database. We find that normalizing all risk parameters with their respective characteristic values yields reasonable correlations with different variants of gravity number. All correlations confirm the physics behind plume movement in a reservoir. The correlations reproduce almost all simulation results within a factor of two, and this is adequate for rapid ranking or screening of prospective storage reservoirs. CO₂ injection in saline reservoirs on the scale of tens of millions of tonnes may result in fracturing, fault activation and leakage of brine along conductive pathways. Critical contour of overpressure (CoP) is a convenient proxy to determine the risk associated with pressure buildup at different location and time in the reservoir. The location of this contour varies depending on the target aquifer properties (porosity, permeability etc.) and the geology (presence and conductivity of faults). The CoP location also depends on relative permeability, and we extend the three-region injection model to derive analytical expressions for a specific CoP as a function of time. We consider two boundary conditions at the aquifer drainage radius, constant pressure or an infinite aquifer. The model provides a quick tool for estimating pressure profiles. Such tools are valuable for screening and ranking sequestration targets. Relative permeability curves measured on samples from seven potential storage formations are used to illustrate the effect on the CoPs. In the case of a constant pressure boundary and constant rate injection scenario, the CoP for small overpressures is time-invariant and independent of relative permeability. Depending on the relative values of overall mobilities of two-phase region and of brine region, the risk due to a critical CoP which lies in the two-phase region can either increase or decrease with time. In contrast, the risk due to a CoP in the drying region always decreases with time. The assumption of constant pressure boundaries is optimistic in the sense that CoPs extend the least distance from the injection well. We extend the analytical model to infinite-acting aquifers to get a more widely applicable estimate of risk. An analytical expression for pressure profile is developed by adapting water influx models from traditional reservoir engineering to the "three-region" saturation distribution. For infinite-acting boundary condition, the CoP trends depend on same factors as in the constant pressure case, and also depend upon the rate of change of aquifer boundary pressure with time. Commercial reservoir simulators are used to verify the analytical model for the constant pressure boundary condition. The CoP trends from the analytical solution and simulation results show a good match. To achieve safe and secure CO₂ storage in underground reservoirs several state and national government agencies are working to develop regulatory frameworks to estimate various risks associated with CO₂ injection in saline aquifers. Certification Framework (CF), developed by Oldenburg et al (2007) is a similar kind of regulatory approach to certify the safety and effectiveness of geologic carbon sequestration sites. CF is a simple risk assessment approach for evaluating CO₂ and brine leakage risk associated only with subsurface processes and excludes compression, transportation, and injection-well leakage risk. Certification framework is applied to several reservoirs in different geologic settings. These include In Salah CO₂ storage project Krechba, Algeria, Aquistore CO₂ storage project Saskatchewan, Canada and WESTCARB CO₂ storage project, Solano County, California. Compositional reservoir simulations in CMG-GEM are performed for CO₂ injection in each storage reservoir to predict pressure build up risk and CO₂ leakage risk. CO₂ leakage risk is also estimated using the catalog of pre-computed reservoir simulation results. Post combustion CO₂ capture is required to restrict the continuous increase of carbon content in the atmosphere. Coal fired electricity generating stations are the dominant players contributing to the continuous emissions of CO₂ into the atmosphere. U.S. government has planned to install post combustion CO₂ capture facility in many coal fired power plants including W.A. Parish electricity generating station in south Texas. Installing a CO₂ capture facility in a coal fired power plant increases the capital cost of installation and operating cost to regenerate the turbine solvent (steam or natural gas) to maintain the stripper power requirement. If a coal-fired power plant with CO₂ capture is situated over a viable source for geothermal heat, it may be desirable to use this heat source in the stripper. Geothermal brine can be used to replace steam or natural gas which in turn reduces the operating cost of the CO₂ capture facility. High temperature brine can be produced from the underground geothermal brine reservoir and can be injected back to the reservoir after the heat from the hot brine is extracted. This will maintain the reservoir pressure and provide a long-term supply of hot brine to the stripper. Simulations were performed to supply CO₂ capture facility equivalent to 60 MWe electric unit to capture 90% of the incoming CO₂ in WA Parish electricity generating station. A reservoir simulation study in CMG-GEM is performed to evaluate the feasibility to recycle the required geothermal brine for 30 years time. This pilot study is scaled up to 15 times of the original capacity to generate 900 MWe stripping system to capture CO₂ at surface. / text

In spite of the large number of research works on carbon capture and sequestration (CCS), the migration and behaviour of CO2 in the subsurface (i. e. strata below the earth's surface) still needs further understanding and investigations with the aim of encouraging the governmental policy makers to adopt CCS technology as one of the most viable means to tackle the global warming threats. In this research work, a series of numerical simulations has been carried out using STOMP-CO2 simulation code to determine the flow behaviour and ultimate fate of the injected supercritical carbon dioxide (scCO2) into saline aquifers in medium terms of storage (i. e. few thousand years). The characteristics of the employed simulator, including the mathematical algorithm, governing equations, equations of states and phase equilibria calculations are explained in details.

Of all the strategies to reduce carbon emissions, carbon dioxide (CO₂) geological sequestration is an immediately available option for removing large amounts of the gas from the atmosphere. However, our understanding of the transition behavior between Forchheimer and Darcy flow through porous media during CO₂ injection is currently very limited. In addition, the kinetic mass transfer of SO₂ and CO₂ from CO₂ stream to the saline and the fully coupling between the changes of porosity and permeability and multiphase flow are two significant dimensions to investigate the brine acidification and the induced porosity and permeability changes due to SO₂ co-injection with CO₂. Therefore, this dissertation develops a multiphase flow, contaminant transport and geochemical model which includes the kinetic mass transfer of SO₂ into deep saline aquifers and obtains the critical Forchheimer number for both water and CO₂ by using the experimental data in the literature. The critical Forchheimer numbers and the multiphase flow model are first applied to analyze the application problem involving the injection of CO₂ into deep saline aquifers. The results show that the Forchheimer effect would result in higher displacement efficiency with a magnitude of more than 50% in the Forchheimer regime than that for Darcy flow, which could increase the storage capacity for the same injection rate and volume of a site. Another merit for the incorporation of Forchheimer effect is that more CO₂ would be accumulated in the lower half of the domain and lower pressure would be imposed on the lower boundary of the cap-rock. However, as a price for the advantages mentioned above, the injection pressure required in Forchheimer flow would be higher than that for Darcy flow. The fluid flow and contaminant transport and geochemical model is then applied to analyze the brine acidification and induced porosity and permeability changes due to SO₂ co-injection. The results show that the co-injection of SO₂ with CO₂ would lead to a substantially acid zone near the injecting well and it is important to include the kinetic dissolution of SO₂ from the CO₂ stream to the water phase into the simulation models, otherwise considerable errors would be introduced for the equilibrium assumption. This study provides a useful tool for future analysis and comprehension of multiphase Darcy-Forchheimer flow and brine acidification of CO₂ injection into deep saline aquifers. Results from this dissertation have practical use for scientists and engineers concerned with the description of flow behavior, and transport and fate of SO₂ during SO₂ co-injection with CO₂ in deep saline aquifers.