The Upper Cretaceous Qingshankou Formation's lacustrine shales in the Songliao Basin are among China's most promising shale oil reservoirs. To elucidate their pore and fractal characteristics, a comprehensive set of analyses encompassing total organic carbon (TOC), X-ray diffraction (XRD), and low-temperature N2 adsorption (LTNA), Rock-Eval pyrolysis experiments and two-dimensional nuclear magnetic resonance (2D-NMR) were conducted. Using the Frenkel-Halsey-Hill (FHH) method, fractal dimensions (D) were calculated, and their relationship with pore metrics and shale compositions were explored. Two distinct fractal dimensions, D1 (0 < P/P0 < 0.5) and D2 (0.5 <P/P0 <1.0), were derived from LTNA isotherms via the FHH approach. D1 values fluctuated between 2.5715 and 2.7551 (mean 2.6564), while D2 spanned from 2.3247 to 2.4209 (mean 2.3653). Notably, D1 consistently surpassed D2, signifying that smaller pores exhibit greater homogeneity compared to their larger counterparts. D1 gradually increases with the increase of clay content. A direct correlation was observed between pore volume (PV), specific surface area (SSA), and D (both D1 and D2), whereas the association between average pore diameter (APD) and D was inverse. Both D1 and D2 escalated with diminishing TOC, 2D-NMR solid organic matter (OM), S1 content and 2D-NMR light oil. Intriguingly, D1 showed a stronger association with key pore and "sweet spot" parameters, highlighting its utility in assessing pore structural complexity and shale oil potential. This study illustrates how fractal theory enhances our understanding of pore structures and the shale oil enrichment process for the lacustrine shale.