A new method to investigate the thermodynamic properties of electrolyte solutions is presented on the perturbation theory of fluids, where a mixture of neutral hard spheres and charged hard spheres is used as the reference system. The interaction forces of dispersion, induction, dipole and quadrupole are taken as perturbation potential. Explicit expressions for the Helmholtz free energy and internal energy are deduced and tested with MC data.