A wellbore and a combined reservoir system are essential for the management of subsurface fluid resources and the geological storage of CO 2 . But the interaction between wellbore and reservoir flow is often neglected in studies of the combined system. A 2D radial model, considering the interaction of wellbore and reservoir flow was developed to investigate its impact on CO 2 geological sequestration. The mass, energy and momentum equations for the wellbore and reservoir were solved using T2Well/ECO2N. The gas flow rate of the reservoir and wellbore are predicted, and the impact of interaction between wellbore and reservoir flow on the CO 2 plume distribution and evolution was investigated. Furthermore, the influence of the CO 2 injection rate, reservoir properties and salinity on the distribution of wellbore and reservoir flow was also explored. Interaction between the wellbore and reservoir flows determines the distribution of the reservoir gas flow rate which combined with layer thickness and porosity controls the horizontal distribution and evolution of the CO 2 plume. The CO 2 wellhead injection rate and reservoir properties (including lateral transmissivity, permeability) are vital factors influencing wellbore and reservoir flows. However, reservoir salinity has little effect on the interaction between the wellbore flow and the reservoir flow, but increased reservoir salinity can accelerate the horizontal migration of CO 2 . The results of this study may help to change the widely held opinion that the distribution of the injected CO 2 among the individual layers is simply proportional to their transmissivity, and thereby enhance our understanding of CO 2 evolution beneath the surface and provide theoretical support for safe and potential geological storage of CO 2 .
List of symbols KThe index for components, k = 1 (H 2 O), 2 (salt), 3 (CO 2 ), and 4 (energy) M kThe accumulation terms of the components and energy k, kg m -3 q k Source/sink terms for mass or energy components, kg m -3 s -1 F k The mass or energy transport terms along the borehole due to an advective process, W m -1 X k b Mass fraction of component k in fluid phase b,The Darcy's velocity in phase b, m s -1 AThe well cross-sectional area, m 2 z Distance along-wellbore coordinate (can be vertical, inclined, or horizontal), m U bThe internal energy of phase b per unit mass, J kg -1The momentum of phase b per unit mass, J kg -1 k Thermal conductivity, W K -1 m -1 h bThe specific enthalpy of phase b, J kg -1 gGravitational acceleration constant, m s -2 hThe angle of inclination of the wellbore, dimensionlessThe wellbore heat loss/gain per unit length of wellbore, kg m -3 s -1 T Temperature, C, K TTime, s C 0The profile parameter (or distribution coefficient), dimensionless u d Drift velocity, m s -1 q mThe density of the gas-liquid mixture, kg m -3 u mThe mixture velocity (velocity of the mixture mass center), m s -1 q m * The profile-adjusted average density, kg m -3 u inThe flow in mass flow rate of wellbore block, kg s -1 u outThe flow out mass ...