This paper proposes a novel two-dimensional (2D) semi-analytical method for calculating ground vibrations from a tunnel situated in a homogeneous half-space with an irregular surface. The circular tunnel is conceptualized as an elastic solid, while the soil is modelled as an elastic, isotropic, and homogeneous half-space with an irregular surface. A virtual horizontal interface is introduced to divide the soil domain into an irregular region with an arbitrary-shaped surface and a half-space with a circular hollow. The wavefield scattered by the irregular surface is simulated by the boundary integral equation. Through the application of the discrete wavenumber approach and expansion of wavefield on the irregular surface, the boundary integral equation is transformed into a set of simultaneous matrix equations containing unknown expansion coefficients. By utilizing the transformation between cylindrical and plane waves, the boundary conditions on soil-tunnel interface are satisfied, yielding a solution for a harmonic point load acting on the tunnel in a half-space with an irregular surface. The proposed method can simulate topographies with arbitrary shapes and only requires the discretization of the irregular part of the ground surface, which provides a highly efficient tool to investigate the propagation characteristics of train-induced vibrations in complex topographies. The proposed method is verified through comparisons with the existing analytical methods and the finite element method. The dynamic responses of a tunnel embedded in a half-space with a Gaussian-shaped surface are meticulously investigated. A case study involving a site for the Shenzhen metro line 5 is performed. The numerical results demonstrate that the presence of irregular topography alters the distribution of ground vibrations, leading to significant differences in responses within the irregular region and certain disparities in the regular region. The effect of the location, size, and undulating direction of the irregular topography on ground vibrations is frequency-dependent.