Relating the ground motion intensity measure (IM) and the structural engineering demand parameter is a crucial step in the performance-based earthquake engineering framework. This study investigates the selection of IM for development of probabilistic seismic demand model of urban shield tunnels subjected to earthquake ground motions in liquefiable and non-liquefiable soils. Nonlinear dynamic effective stress analyses are conducted to develop a database of the intensity measures and structural seismic responses exposed to ground shaking and soil liquefaction. Two advanced soil constitutive models (i.e., Pressure DependMultiYield03 and PressureIndependMultiYield for liquefiable and non-liquefiable soils, respectively) are employed to capture the nonlinear behavior. A suite of 23ground motion intensity measures is selected and assessed based on the evaluation criteria of correlation, efficiency, practicality and proficiency. Eventually, the multi-level fuzzy comprehensive evaluation method is employed to comprehensively consider the four evaluation criteria and establish the optimal ground motion IM suitable for probabilistic seismic demand analysis of shield tunnel structures. The obtained results show that the sustained maximum acceleration is the optimal IM for evaluating the structural seismic response, followed by the peak ground acceleration in both liquefiable and non-liquefiable soils. Peak pseudo velocity spectrum, displacement square integral and Housner spectral intensity are found to be not suitable for the probabilistic seismic demand analysis of shield tunnel structures.