Solvent primitive model for electrolyte solutions in disordered porous matrices of charged species. Replica Ornstein-Zernike theory and grand canonical Monte Carlo simulations