Large diameter shield tunnels traversing liquefiable soil-rock strata are highly susceptible to seismic hazards, as earthquake-induced soil liquefaction significantly reduces soil strength and stiffness. Therefore, it is crucial to accurately assess the seismic performance of these tunnels. This study first establishes a numerical model for tunnel seismic response analysis, considering soil liquefaction, segment nonlinearity, and joint deformation. The validity of the model is affirmed through experimental, theoretical, and additional numerical simulations. The probabilistic seismic demand models are established employing the seismic database consisting of 120 ground motion records. Subsequently, a quantitative selection method for the optimal Intensity Measure (IM) based on fuzzy comprehensive evaluation is proposed, identifying Velocity Spectrum Intensity (VSI) as the most suitable among 29 commonly used IMs, and the IMs related to duration exhibit poor performance. The study then categorizes tunnel damage into three states: minor, moderate, and extensive, using joint opening as the damage measure. Finally, seismic fragility analysis is employed to assess seismic performance of tunnel, and fragility curves derived using VSI and Peak Ground Acceleration (PGA) is compared. The results indicate that PGA, a commonly used IM, significantly underestimates the probability of damage to the tunnel, with a maximum underestimation of 22.4%.
