To get deep insight into the variability and heterogeneity of pore size distributions (PSDs) in different tectonically deformed coals (TDCs) (mean maximum vitrinite reflectance R-o,R-max ranging from 1.720% to 1.857%) collected from Hancheng Mine, Weibei coalfield, the multifractal analysis using generalized dimensions D-q was employed to study PSDs of TDCs based on mercury intrusion porosimetry. The mercury intrusion data under pressure over 10 MPa was corrected in combination with gas adsorption data. The results show that, for the test coals, tectonic deformation (TD) mainly increases seepage-porosity while maceral composition (MC) is the main control on adsorption-porosity. The variations of D-q versus q show that the PSDs of TDCs exhibit multifractal behavior, however, differ in their multifractality suggested by the extracted parameters from PSDs including information dimension D-1, the Hurst exponent H, the width of the right side D-0-D-10 and the left side D-(10)-D-0 of D-q spectra. TD and MC have great influence on the variability and heterogeneity in the inner distribution of seepage-porosity and adsorption-porosity respectively, as demonstrated by the change of D-1, H and D-0-D-10 for moment q > 0 and D-(10)-D-0 for moment q < 0. Both D-1 and H increase with TD but D-0-D-10 shows a reverse trend. Consequently, TD leads to narrower distribution with higher fluctuation, lower pore connectivity and greater complexity in the distribution of seepage-pores, which may be due to the presence and uneven distribution of microfracture and inter-granular pores formed through deformation process. Increasing vitrinite and decreasing inertinite contribute to a heterogeneous distribution of adsorption-porosity, as observed by the high correlation between D-(10)-D-0 and MC. The aforementioned results confirm that the multifractal approach is useful to characterize the internal heterogeneity of PSDs and to detailedly distinguish between PSDs of TDCs. (C) 2015 Elsevier B.V. All rights reserved.