This study presents a comprehensive analysis of rock properties for a selected group of six carbonate reservoir rock samples revealing complex structures at various length scales. Experimental laboratory methods as well as image analysis techniques were conducted in this study to characterize the macro- and micro-pores in mud- and grain-dominated limestones samples from the Upper Jurassic Arab Formation (Arab D member). Mercury Injection Capillary Pressure (MICP), porosimeter, and permeameter lab measurements were employed to assess the pore network heterogeneity and complexity. In addition, a multiscale rock imaging approach was implemented to detect grain and pore phases at several length scales using Thin Sections (TS), Scanning Electron Microscopy (SEM), Focused Ion Beam Scanning Electron Microscopy (FIB-SEM), as well as 3D X-ray Computed Tomography (CT), and micro-computed tomography images (MCT). Furthermore, the multifractal analysis method was applied on the MICP and FIB-SEM to characterize quantitatively the heterogeneity of the pores in the studied samples. Heterogeneous samples 3R, 4M, 5W, and 6M display the highest non-uniformity degree Δα values, falling within the range of [1.21, 1.39] based on FIB-SEM images. Samples 1G, 2R, 3R, and 5W exhibit more heterogeneous pore structure, with Δα values ranging from 0.73 to 1.49 based on the MICP experiments. The results and findings confirm the effectiveness of multifractal parameters Δα and the asymmetry degree in the vertical axis Δf(α) in quantifying and characterizing rock heterogeneity.