The volcanic rocks of the Es3 Formation (the third member of the Eocene Shahejie) in the Liaohe Eastern Sag can be divided into four facies and twelve subfacies. The porosity spectrum, porosity bin, variation coefficient (VC), and porosity width derived from electrical imaging log data were applied to study and characterize the heterogeneity of four facies and nine subfacies, both qualitatively and quantitatively. However, the VC and porosity width cannot be used to quantitatively classify heterogeneity when the VC is small and the porosity width is large. In the present study, the authors propose a new parameter, Pvcd = variation coefficient × porosity width. Based on this parameter, using a combination of porosity spectra, porosity bin features, VC, and porosity width, lithofacies heterogeneity is divided into three categories. The first is weak heterogeneity, which has a Pvcd < 1.1, a VC < 0.15, and a porosity width < 6. The second is moderate heterogeneity, which has a Pvcd 1.1–1.6, a VC 0.15–0.25, and a porosity width 6–9. The third is strong heterogeneity, which has a Pvcd > 1.6, a porosity VC > 0.25, and a porosity width > 9. In these three cases, the porosity spectra mainly display unimodal features, the porosity bins are concentrated, the bimodal features with tails and porosity bins are less concentrated, or the multimodal features with tails and porosity bins are scattered, respectively. Favorable reservoirs of volcanic rocks are controlled by lithofacies. In the study area, the favorable reservoirs appear to be the diatreme subfacies with low or medium heterogeneity, the pyroclastic flow subfacies with low heterogeneity, the compound lava flow subfacies with low or medium heterogeneity, and the outer zone subfacies with strong heterogeneity.