Several laboratory experiments demonstrated that different water compositions cause rocks to change from oil- to water-wet state. Although it is a consensus that wettability alteration is the main recovery mechanism, modeling the underlying mechanism is still a major challenge. Our main goal is to improve and validate a physically based model to predict contact angles from zeta-potential measurements. We propose a new mass-action formulation for surface complexation model (SCM) that includes the energy interaction effect between two close surfaces (PS). Currently, most SCMs consider rock and oil as isolated surfaces (IS). Thus, we hypothesize that, as rock and oil surface approach each other, PS model produce a better description of electrostatic distribution. Additionally, we develop a method of determining SCM equilibrium constants to fit several zeta-potential measurements for different ion concentrations (Na+, Ca2+, Mg2+, SO42- and H+). Finally, we estimate contact angles using disjoining pressure calculations and compare them with ones reported in the literature. From a SCM set of reactions available in the literature, we validate the developed IS model against PHREEQC (a reference simulator for geochemical reactions). For the PS case, the system of equations’ solution is very close to IS approach when the interaction between surfaces are negligible (wide spacing between surfaces). Regarding zeta-potential prediction for calcite-brine system, we argue that Na+ might not be an indifferent ion as suggested previously. Our simulation results indicate that, besides the renowned potential-determining ions, sodium adsorption on calcite can play an important role in electrostatic interactions, switching surface charge polarity. Thus, we only achieve a successful fit of zeta-potential measurements when Na+ is considered in the SCM reactions. Finally, contact angle estimation using the PS model and disjoining pressure theory provide good predictions of seven different cases reported in the literature. We validate our method on a total of 66 and 163 contact angle and zeta-potential measurements, respectively. The present work is a novel approach to represent how electrostatic interactions among rock, brine and oil modify the rock surface charge and the rock wetting state.