Abstract The diffusiophoretic and electrophoretic motions of two colloidal spheres in the solution of a symmetrically charged electrolyte are analyzed using a method of reflections. The particles are oriented arbitrarily with respect to the electrolyte gradient or the electric field, and they are allowed to differ in radius and in zeta potential. The thickness of the electric double layers surrounding the particles is assumed to be small relative to the radius of each particle and to the gap width between the particles, but the effect of polarization of the mobile ions in the diffuse layer is taken into account. A slip velocity of fluid and normal fluxes of solute ions at the outer edge of the thin double layer are used as the boundary conditions for the fluid phase outside the double layers. The method of reflections is based on an analysis of the electrochemical potential and fluid velocity disturbances produced by a single dielectric sphere placed in an arbitrarily varying electrolyte gradient or electric field. The solution for two-sphere interactions is obtained in expansion form correct to O(r12−7), where r12 is the distance between the particle centers. Our analytical results are found to be in good agreement with the available numerical solutions obtained using a boundary collocation method. On the basis of a model of statistical mechanics, the results of two-sphere interactions are used to analytically determine the first-order effect of the volume fraction of particles of each type on the mean diffusiophoretic and eletrophoretic velocities in a bounded suspension. For a suspension of identical spheres, the mean diffusiophoretic velocity can be decreased or increased as the volume fraction of the particles is increased, while the mean electrophoretic velocity is reduced with the increase in the particle concentration. Generally speaking, the particle interaction effects can be quite significant in typical situations.