Geological heterogeneities essentially affect the dynamics of a CO2 plume in subsurface environments. Previously we showed how the dynamics of a CO2 plume is influenced by the multi-scale stratal architecture in deep saline reservoirs. The results strongly suggest that representing small-scale features is critical to understanding capillary trapping processes. Here we present the result of simulation of CO2 trapping using two different conventional approaches, i.e. Brooks-Corey and van Genuchten, for the capillary pressure curves. We showed that capillary trapping and dissolution rates are very different for the Brooks-Corey and van Genuchten approaches when heterogeneity and hysteresis are both represented.