Impact of Grid Non-Orthogonality on Co2 Migration and Pressure Prediction Using Tpfa and Mpfa Methods

Abstract

Accurate simulation of multiphase flow is essential for predicting pressure buildup and CO2 plume migration in geological carbon storage. The two-point flux approximation (TPFA) method is widely used in reservoir simulation due to its simplicity, but it is only consistent on K-orthogonal grids. Realistic subsurface models, however, often involve non-orthogonal grids caused by geological features such as faults and fractures, where TPFA may introduce significant errors. This study systematically quantifies the numerical error associated with using TPFA on non-K-orthogonal grids in the context of CO2 storage in saline aquifers. A series of numerical experiments are conducted to compare the TPFA and multi-point flux approximation (MPFA) methods, with high-resolution reference solutions used as benchmarks. Results show that grid non-orthogonality can impact the accuracy of both TPFA and MPFA solutions with larger degree of grid inclination leading to larger solution errors in general. The magnitude of solution errors of MPFA method is significantly smaller than that of TPFA in terms of both pressure solution and CO2 saturation and MPFA method is much more robust against grid non-orthogonality effect. In particular, the TPFA method can produce substantial deviations in CO₂ saturation distributions when grid non-orthogonality is present, indicating the necessity of more advanced discretization methods such as MPFA for modelling the CO2 plume migration behavior more accurately. However, our results also show that for highly heterogenous formations, the permeability/porosity heterogeneity plays a predominant role and the effect of grid non-orthogonality becomes less pronounced. These findings highlight the importance of selecting appropriate discretization methods for geologically complex reservoirs and offer practical guidance on the trade-offs between computational efficiency and physical accuracy in carbon storage modeling.

Keywords: geological CO2 storage, numerical simulation, Non-orthogonality, TPFA, MPFA