Systematic DFT calculations of model coal-pyrrole derivatives substituted by different functional groups are carried out. The N-H bond dissociation energies (N-H BDEs) and H-transfer activation energies (H-TAEs) of pyrrole derivatives are fully evaluated to elucidate the effect of the type of substituents and their position on the molecular reactivity. The results indicate that compounds substituted with electron-donating groups (EDGs) are more prone to pyrolysis while those substituted with electron-withdrawing groups (EWGs) are difficult to pyrolyze. Furthermore, quantitative structure-property relationship (QSPR) models for N-H BDEs and H-TAEs about pyrrole derivatives are built with multiple linear regression (MLR) and support vector machine (SVM). The final results show that the SVM-QSPR model has better fitness, prediction, and robustness, while the MLR-QSPR model can express the physical meaning better. The effects of functional groups on pyrolysis are clarified by the models presented in this paper, which will support the optimization of ultra-low NOx combustion.