Information on the sulfur classes present in petroleum is a key factor in determining the value of refined products and processing behavior in the refinery. A large part of the sulfur present is included in polycyclic aromatic sulfur heterocycles (PASHs), which in turn are difficult to desulfurize. Furthermore, some PASHs are potentially more mutagenic and carcinogenic than polycyclic aromatic hydrocarbons, PAHs. All of this calls for improved methods for the identification and quantification of individual sulfur species. Recent advances in analytical techniques such as comprehensive two-dimensional gas chromatography (GC×GC) have enabled the identification of many individual sulfur species. However, full identification of individual components, particularly in virgin oil fractions, is still out of reach as standards for numerous compounds are unavailable. In this work, a method for accurately predicting retention times in GC×GC using a QSRR (quantitative structure retention relationship) method was very helpful for the identification of individual sulfur compounds. Retention times for 89 saturated, aromatic, and polyaromatic sulfur-containing heterocyclic compounds were determined using two-dimensional gas chromatography. These retention data were correlated with molecular descriptors generated with CODESSA software. Two independent QSRR relationships were derived for the primary as well as the secondary retention characteristics. The predictive ability of the relationships was tested by using both independent sets of compounds and a cross-validation technique. When the corresponding chemical standards are unavailable, the equations developed for predicting retention times can be used to identify unknown chromatographic peaks by matching their retention times with those of sulfur compounds of known molecular structure.