Probabilistic analysis of tunnels: A hybrid polynomial correlated function expansion based approach