.[1] Most simulators for subsurface flow of water, gas, and oil phases use empirical correlations, such as Henry's law, for the CO 2 composition in the aqueous phase, and equations of state (EOS) that do not represent the polar interactions between CO 2 and water. Widely used simulators are also based on lowest-order finite difference methods and suffer from numerical dispersion and grid sensitivity. They may not capture the viscous and gravitational fingering that can negatively affect hydrocarbon (HC) recovery, or aid carbon sequestration in aquifers. We present a three-phase compositional model based on higher-order finite element methods and incorporate rigorous and efficient three-phase-split computations for either three HC phases or water-oil-gas systems. For HC phases, we use the Peng-Robinson EOS. We allow solubility of CO 2 in water and adopt a new cubic-plus-association (CPA) EOS, which accounts for cross association between H 2 O and CO 2 molecules, and association between H 2 O molecules. The CPA-EOS is highly accurate over a broad range of pressures and temperatures. The main novelty of this work is the formulation of a reservoir simulator with new EOS-based unique three-phase-split computations, which satisfy both the equalities of fugacities in all three phases and the global minimum of Gibbs free energy. We provide five examples that demonstrate twice the convergence rate of our method compared with a finite difference approach, and compare with experimental data and other simulators. The examples consider gravitational fingering during CO 2 sequestration in aquifers, viscous fingering in water-alternating-gas injection, and full compositional modeling of three HC phases.Citation: Moortgat, J., Z. Li, and A. Firoozabadi (2012), Three-phase compositional modeling of CO 2 injection by higher-order finite element methods with CPA equation of state for aqueous phase, Water Resour. Res., 48, W12511,