Mean-field theory for the conserved charges in Richardson-Gaudin models