Monitoring field-scale CO2 geological utilizations is of paramount importance but challenging due to the complexities in microscopic heterogeneity. In this paper, complex geological characterizations and fluid properties are specifically analysed and a prediction model is developed, which is capable to track the dynamic behaviour of miscible CO2 multiphase flow in microscopically-heterogeneous porous media. More specifically, first, pore-throat sizes and distributions were characterised from the constant-rate mercury injection and the CO2-displacement seepage resistance was evaluated from a capillary bundle model. The differences of seepage resistance caused from the throat changing and coupling diffusion-dissolution effects and viscosity-resistance reductions were specifically studied in the process of continuous miscible CO2 displacement. Accordingly, a comprehensive mathematical model was developed with the time-node analysis method and validated through comparison with the experimental results. The leading edge of CO2 displacement as well as the timing of gas breakthrough and displacement completion are determined to be different with varying throat sizes. It is further found that the gas breakthrough time and the time required to complete the displacement are reduced with increasing throat size and their differences also decrease. Moreover, the interval area of injection and production wells could be characterised as pure CO2, diffusion and pure oil zones, whose positions could be dynamically tracked from the recovery performance. The calculated oil recovery and gas-oil ratio from the developed model share good agreement with experimental results, with deviations of 2.2 and 0.7%. Such the validated mathematical model is then applied to successfully predict the dynamic performance of an actual reservoir (H3).