The objective of this study was to investigate the effect of wettability heterogeneity on pore-scale characteristics of Supercritical (sc) CO2 displacement dynamics and its capillary trapping mechanism during a scCO(2)-brine drainage and imbibition cycle. A multiphase lattice Boltzmann (LB) model was employed to simulate scCO(2)-brine flow in rock samples of Tuscaloosa sandstone taken from the Cranfield CO2 injection site. Using a spectral method, we adopted various wettability fields to generate rock samples containing distributed CO2-wet regions. To gain a better insight into the effect of fractional wettability on scCO(2) displacement patterns during drainage, we quantified the evolution of scCO(2) interface with brine and rock surface for samples with various wettability heterogeneities. In addition, the effect of heterogeneous wettability on the drainage relative permeability and capillary pressure curves has been investigated in this study. According to our results, heterogeneous distribution of CO2-wet regions in the rock leads to more dispersed fluid distribution and, hence, more tortuous flow paths, resulting in higher interfacial area between fluid phases and rock surface at any given scCO(2) saturation. Furthermore, the spatial distribution of wettability controls the scCO(2) entrapment pattern and spatial distribution of residual scCO(2) clusters during brine flooding. In fractional-wet samples, residence of scCO(2) phase in CO2 -wet regions creates more trapped scCO(2) clusters, suppressing the connectivity of the CO2 phase, thus enhancing more residual trapping. Our results imply that the total number of scCO(2) clusters and, as a result, their residual trapping, increases as the fraction of CO2-wet regions becomes larger, leading to a larger surface area of scCO(2) with brine and rock surface, potentially, facilitating the likelihood of long-term dissolution and mineral trapping.