The first part proposed a two-phase flow-reservoir coupling model

in order to apply this model to evaluate the stability of reservoir and fluid transportation

the commercial algorithm

Comsol Multiphysics

is used and the corresponding numerical modeling method is proposed. The determination procedure of model parameters is firstly presented and the simulation results are validated by the laboratory tests. The density and viscosity of CO2 are determined according to the in situ condition of temperature and pressure

the Van Genuchten parameters of rock layer are from the literature. Based on the triaxial compression tests

effective stress coefficient tests and permeability tests

the model parameters involved in mechanical model and coupling relations are identified. We then perform a case study. The rock layer is situated at depth between 680 to 700 m

and the width is 100 m. The evolution of injection pressure under different injection rates is studied

the distributions of pore pressure

vertical strain and damage variable are presented; the CO2 transportation is also analyzed. The present work offers theoretical basis for evaluating the stability of reservoir and fluid transportation.