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　23　ground　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.　©　2024　Tongji　University