Abstract: The water surface wave radiation problem caused by multiple cylinders oscillating with identical frequency was solved in frequency domain by the boundary element method using simple Green's function in the inner water region combined with the eigenfunction expansions in the outer water region. The numerical method is suitable to the situation of constant depth of outer regions and complicated boundary conditions of inner region, while the oscillating modes, motion amplitudes and phases of the cylinders may be different from one another. The second order potential and hydrodynamic forces acting on ecah cylinder were evaluated completely by perturbation method. Compared with the case of single oscillating cylinder, hydrodynamic interference phenomena, such as wave resonance and negative added mass, of the radiation problem due to the oscillatory motions of multiple cylinders are identified which is of engineering importance to the design of moorings and other facilities involving multiple structures.