Carbon dioxide (CO₂) dissolution is the secondary trapping mechanism enhancing the long-term security of CO₂ in confined geological formations. CO₂ injected into a randomly multilayered formation will preferentially migrate along high permeability layers, increasing CO₂ dissolution efficiency. In this study, sequential Gaussian simulation is adopted to construct the stratified saline formations, and two-phase flow based on MRST is established to illustrate the spatial mobility and distribution of CO₂ migration. The results show that gravity index $G$ and permeability heterogeneity ${\sigma }_Y^2$ are the two predominant factors controlling the spatial mobility and distribution of CO₂ transports. The CO₂ migration shows a totally different spatial mobility under different gravity index and heterogeneity. When the permeability discrepancy is relatively larger, CO₂ preferentially migrates along the horizontal layer without accompanying the vertical migration. For the formation controlled by gravity index, CO₂ migration is governed by supercritical gaseous characteristics. For the medium gravity index, the upward and lateral flow characteristics of the CO₂ plume is determined by gravity index and heterogeneity. When the gravity index is smaller, permeability heterogeneity is the key factor influencing CO₂ plume characteristics. Permeability heterogeneity is the decisive factor in determining final CO₂ dissolution efficiency. This investigation of CO₂ mobility in randomly multilayered reservoirs provides an effective reference for CO₂ storage.