Multiple analytical methods and Monte Carlo simulations were performed to evaluate neutron penetra-tion in straight and curved labyrinths. Factors studied included variations in beam losses of off-axis point source, on-axis point source, and line source. For the straight labyrinth, it was found that the analytical expressions neglect the dose rate platform appearing at the bend of the labyrinth, and the agreement between analytical methods and Monte Carlo estimation was related to the type of neutron source term. For the curved labyrinth, the neutron attenuation length obtained under different conditions was nearly identical and appeared to be in quite good accord with the empirical formula calculation. Moreover, the neutron energy spectra along the centerline distance of the labyrinth were also analyzed. In the first leg, differences in beam loss led to variance in the distribution of spectra, while in the second and subsequent legs, the spectra were similar, where the main contributors were thermal neu-trons. This work is valuable for practical design of the labyrinths in the accelerator facilities.