Artificial ground freezing method is widely applied in the construction of metro tunnel and significantly impact the microstructure of soils in artificial frozen-walls. To delve into the pore distribution characteristics of thawed residual soils, Nuclear Magnetic Resonance Imaging (NMRI) measurements were performed to investigate the relaxation time (T₂) spectrums and T₂-weighted images of saturated samples after freezing at different temperatures. The pore volume distributions were determined from T₂ spectrums based on the surface relaxation coefficient ( ${\rho }_2$ ) and the pore structures were visualized by T₂-weighted images. Subsequently, the pore size distribution curves from NMRI were compared and validated by mercury intrusion porosimetry (MIP) tests. According to the results, the peak areas of T₂ spectrums were linearly related to freezing temperatures in a positive manner. Pore volume distribution curves of thawed soils have two peaks, which are the major peaks with diameters of 0.5–20 μm and the secondary peaks with diameters of 20–500 μm. As the freezing temperature drops, the volumes of pores with different diameters all increased. The damage degree of microstructure in thawed soils increases as the temperature drops, according to the visualized pore structure. Besides, NMRI measurements of saturated soils are more accurate to reflect the full diameter range of pores, compared to MIP method.