In this article, we present a classical density functional theory for electrical double layers of spherical macroions that extends the capabilities of conventional approaches by accounting for electrostatic ion correlations, size asymmetry, and excluded volume effects. The approach is based on a recent approximation introduced by Hansen-Goos and Roth for the hard sphere excess free energy of inhomogeneous fluids [J. Chem. Phys. 124, 154506 (2006); J. Phys.: Condens. Matter 18, 8413 (2006)]. It accounts for the proper and efficient description of the effects of ionic asymmetry and solvent excluded volume, especially at high ion concentrations and size asymmetry ratios including those observed in experimental studies. Additionally, we utilize a leading functional Taylor expansion approximation of the ion density profiles. In addition, we use the mean spherical approximation for multi-component charged hard sphere fluids to account for the electrostatic ion correlation effects. These approximations are implemented in our theoretical formulation into a suitable decomposition of the excess free energy which plays a key role in capturing the complex interplay between charge correlations and excluded volume effects. We perform Monte Carlo simulations in various scenarios to validate the proposed approach, obtaining a good compromise between accuracy and computational cost. We use the proposed computational approach to study the effects of ion size, ion size asymmetry, and solvent excluded volume on the ion profiles, integrated charge, mean electrostatic potential, and ionic coordination number around spherical macroions in various electrolyte mixtures. Our results show that both solvent hard sphere diameter and density play a dominant role in the distribution of ions around spherical macroions, mainly for experimental water molarity and size values where the counterion distribution is characterized by a tight binding to the macroion, similar to that predicted by the Stern model. © 2014 AIP Publishing LLC. [http://dx