A multi-bubble model is developed to study the morphology of a finite array of plasma density depletions (bubbles) in the context of equatorial F-region irregularities during spread F. The Pedersen current conservation equation with quasi-neutrality is solved analytically by using an electrostatic analogy. The solution is exact with no a priori assumption regarding the separation distance. A two-bubble system with a piecewise constant density profile is first analyzed, and the technique is then applied to multi-bubble systems to calculate the polarization electric field and the rise velocities. It is shown that the influence of the neighboring bubbles is relatively short ranged and that a small number of bubbles can adequately model the essential physics in a large array of bubbles. For moderately short separation distances, it is found that the E¿B rise velocity is substantially reduced in comparison with the singel-bubble case and that the rise velocity is strongly sheared resulting in deformation of the contours. It is found that these results compare favorably with observational evidence. The implications of the new morphological results on the stability and dynamical behavior of the bubbles are discussed. The theory is also applied to multi-phase density enhancement such as one might encounter in plasma cloud striation fingers.