We report Monte Carlo simulations of the solvation pressure between two planar surfaces, which represent the interface of spherical silica nanoparticles in supercritical carbon dioxide. Carbon dioxide (CO2) was modeled as an atomistic dumbbell or a spherical Lennard-Jones particle. The interaction between CO2 molecules and silica surfaces was characterized by the standard Steele potential with energetic heterogeneities representing the hydrogen bonds. The parameters for the solid-fluid interaction potentials were obtained by fitting our simulations to the experimental isotherms of CO2 sorption on mesoporous siliceous materials. We studied the dependence of the solvation force on the distance between planar silica surfaces at T = 318 K, at equilibrium bulk pressures P-bulk ranging from 69 to 200 atm. At 69 atm, we observed a long-range attraction between the two surfaces, and it vanished when the pressure was increased to 102 and then 200 atm. The results obtained with different fluid models were consistent with each other. According to our observations, energetic heterogeneities of the surface have negligible influence on the solvation pressure. Using the Derjaguin approximation, we calculated the solvation forces between spherical silica nanoparticles in supercritical CO2 from the solvation pressures between the planar surfaces.