Structural and thermodynamic properties of the restrictive primitive model +1:-1 electrolyte solution adsorbed in a disordered charged media were studied by means of the Grand Canonical Monte Carlo simulation and the replica Ornstein-Zernike theory. Disordered media (adsorbent, matrix) was represented by a distribution of negatively charged hard spheres frozen in a particular equilibrium distribution. The annealed counterions and co-ions were assumed to be distributed within the nanoporous adsorbent in thermodynamic equilibrium with an external reservoir of the same electrolyte. In accordance with the primitive model of electrolyte solutions, the solvent was treated as a dielectric continuum. The simulations were performed for a set of model parameters, varying the net charge of the matrix (i.e., concentrations of matrix ions) and of annealed electrolyte, in addition to the dielectric constant of the invading solution. The concentration of adsorbed electrolyte was found to be lower than the corresponding concentration of the equilibrium bulk solution. This electrolyte "exclusion" depends strongly on the dielectric constant of the invading solution, as also on concentrations of all components. The most important parameter is the net charge of the matrix. Interestingly, the electrolyte rejection decreases with increasing Bjerrum length for the range of parameters studied here. The latter finding can be ascribed to strong inter-ionic correlation in cases where the Bjerumm length is high enough. To a minor extent, the adsorption also depends on the spacial distribution of fixed charges in adsorbent material. The replica Ornstein-Zernike theory was modified to cater for this model and tested against the computer simulations. For the range of parameters explored in this work, the agreement between the two methods is very good. These calculations were also compared with the results of the classical Donnan theory for electrolyte exclusion.