Vibroacoustic analysis of a partially submerged cylindrical shell-liquid coupling system is a typical acoustic-structure interaction problem in an acoustic half-space. Generally, the vibration and acoustic solutions of this problem almost depend on numerical method in previous studies. However, whether from the aspect of verifying the numerical methods or from the aspect of revealing the vibroacoustic mechanism of the acoustic-structure interaction system, the development of an analytical or semi-analytical method is indispensable. In this study, a semi-analytical method is proposed to address the vibroacoustic response of a horizontal cylindrical shell partially immersed in water The acoustic coordinate system is established on the free surface, and the sine series function is introduced to express the sound pressure to meet the pressure release boundary condition on the free surface automatically. Then based on the two-dimensional Flgge shell theory, the motion equation of the shell-liquid coupling system is established with using the center of the shell as the origin of the coordinate system. By using the Galerkin method, the velocity continuous condition on the fluid-structure interface is established. Then the relation matrix between the acoustic pressure amplitude and the shell displacement amplitude is derived after converting the sound pressure from different coordinate systems into that in the same coordinate system by using the geometry relationship of two coordinate systems. Then the vibration and sound radiation of this coupling system are predicted by solving the coupled matrix equations, and the coupled free vibration can also be addressed by solving the characteristic equation after setting the determinant of the coefficient matrix to be zero. To verify the accuracy and reliability of the present method, the coupled natural frequencies, the vibration response and the distribution of the sound pressure obtained from the present method are compared with those obtained from the finite element method, showing that they are in good agreement with each other. Then the shapes of the first four order modes for the coupled system are presented and compared with those of the system submerged in unbounded fluid field, and the couple between different modes, and the couple between the symmetric and antisymmetric modes are observed due to the effect of the free surface. The characteristics of the root mean square velocity of the shell with different immersed depths are discussed, which reveals that the peaks of response curves shift to the lower frequencies with increasing immersed depth. The characteristics of far field sound pressure directivity are presented and explained by the image method in detail. This study provides a novel method to analytically predict the vibroacoustic response of an elastic structure partially coupling with the fluid field when bounded in a sound half space.