ИСТИНА

Two-phase models are widely used for simulation of CO2 storage in saline aquifers. These models support gaseous phase mainly saturated with CO2 and liquid phase mainly saturated with H2O (e.g. TOUGH2 code). The models can be applied to analysis of CO2 storage only in relatively deeply-buried reservoirs where pressure exceeds CO2 critical pressure. At these supercritical reservoir conditions only one supercritical CO2-rich phase appears in aquifer due to CO2 injection. In shallow aquifers where reservoir pressure is less than the critical pressure CO2 can split in two different liquid-like and gas-like phases (e.g. Spycher et al., 2003). Thus a region of three-phase flow of water, liquid and gaseous CO2 can appear near the CO2 injection point. Today there is no widely used and generally accepted numerical model capable of the three-phase flows with two CO2-rich phases. In this work we propose a new hydrodynamic simulator MUFITS (Multiphase Filtration Transport Simulator) for multiphase compositional modeling of CO2-H2O mixture flows in porous media at conditions of interest for carbon sequestration. The simulator is effective both for supercritical flows in a wide range of pressure and temperature and for subcritical three-phase flows of water, liquid CO2 and gaseous CO2 in shallow reservoirs. The distinctive feature of the proposed code lies in the methodology for mixture properties determination. Transport equations and Darcy correlation are solved together with calculation of the entropy maximum that is reached in thermodynamic equilibrium and determines the mixture composition. To define and solve the problem only one function – mixture thermodynamic potential – is required. The potential is determined using a three-parametric generalization of Peng-Robinson equation of state fitted to experimental data (Todheide, Takenouchi, Altunin etc.). We apply MUFITS to simple 1D and 2D test problems of CO2 injection in shallow reservoirs subjected to phase changes between liquid and gaseous CO2. We consider CO2 injection into highly heterogeneous the 10th SPE reservoir. We provide analysis of physical phenomena that have control temperature distribution in the reservoir. The distribution is non-monotonic with regions of high and low temperature. The main phenomena responsible for considerable temperature decline around CO2 injection point is the liquid CO2 evaporation process. We also apply the code to real-scale 3D simulations of CO2 geological storage at supercritical conditions in Sleipner field and Johansen formation (Fig).