The operation and maintenance of metro systems play a crucial role in urban development, with a focus on ensuring the serviceability and safety of metro tunnels. Accurate evaluation of the condition of these tunnel states requires investigating the complex interaction of multiple factors that impact their safety state. This study developed a hybrid model that integrates pair copula constructions (PCCs) and Bayesian networks (BN) to assess the safety state of metro tunnels, considering complex dependencies among these factors. First, key performance indicators (KPIs) were selected to assess tunnel safety, based on six failure modes. Then, an improved copula-based PC algorithm was employed to learn the non-Gaussian bayesian network model, which eliminated the normality assumption in the marginal densities of the KPIs. Finally, a combination of forward reasoning and GM(1,1) was utilized to predict the safety state of the metro tunnel in time series, facilitating the formulation of effective operation and maintenance plans. Furthermore, the proposed approach was applied to a real case study of the Wuhan metro system to demonstrate its effectiveness and applicability. The results highlighted the significant influence of some KPIs, such as convergence deformation, dislocation displacement, and stripping area, on metro tunnel safety. These KPIs emerged as key factors requiring focused attention in the operation and maintenance of metro tunnels.