A hybrid numerical model for coupled hydro-mechanical analysis during CO 2 injection into heterogeneous unconventional reservoirs

This paper presents a hybrid numerical model for simulating compressible multiphase multicomponent transport in deformable, heterogeneous unconventional reservoirs. A fully implicit theoretical framework for coupled hydromechanical behaviour is established, in which the flash calculation is avoided. To preserve the local mass conversation, an element based finite volume method is used to spatially discretize mass conversation equations. Linear momentum equilibrium equation is discretized with finite element method. Two numerical approaches are coupled through sharing the same shape functions. An implicit mid -interval backwards difference algorithm is used to handle the temporal derivative. Reliability and efficiency of the developed model have been examined by verification tests against available analytical or numerical solutions. The coupled hydromechanical processes in heterogeneous reservoirs were simulated. The results of numerical experiments show that oscillatory problem for saturation in heterogeneous media can be eliminated, demonstrating better performance of developed model for heterogeneous reservoir simulations. Reservoir performance is significantly shaped by its heterogeneity through controlling fluid flow and time for CO 2 - CH 4 displacement. Simulation results show that the time required for completion of CO 2 - CH 4 displacement increases by about 250 days when considering heterogeneity effect. In conclusion, this study offers a valuable tool for facilitating the development of CO 2 -enhanced unconventional gas recovery.