Two near-field methods, namely the integral method and differential method, were presented for giving second-order mean drift forces and moments between two fixed submerged bodies in regular waves. For the integral method, with a series of mathematical manipulations, second-order drift forces and moments could be easily expressed by distributed sources which could be calculated by source distribution techniques with the assumption that the amplitude of ship motions are small on the basis of the linear 3D frequency theory. For the differential method, drift forces and moments could be expressed by the derivative of velocity potential with respect to space coordinate. Because two bodies would behave as a single body while the clearance is very large, the numerical results of one sphere in such case were given and compared with analytical results of a single sphere which does not involve the effect of free surface. When submerged depth becomes enough large, a good agreement can be reached. Then the integral method was used to predict the second-order drift forces and moments of two submerged spheres and spheroids with a small lateral separation distance in waves compared with the numerical results obtained by the differential method and they agree well. By comparison, it indicates the interaction effects between two submerged bodies have a profound influence on the drift forces and moments. In this paper, the forward speed effect on submerged spheres was also considered.