This study focuses on the heterogeneity of the middle Miocene syn-rift Belayim nullipore (reefal) marine sequences in the Gulf of Suez and its impacts on reservoir quality. The sequences consist of coralline algal reef limestones with a highly complex dual-porosity system of primary and secondary porosities of widely varying percentages. To achieve a precise mathematical modeling of these reservoir sequences, a workflow protocol was applied to separate these sequences into a number of hydraulic flow units (HFUs) and reservoir rock types (RRTs). This has been achieved by conducting a conventional core analysis on the nullipore marine sequence. To illustrate the heterogeneity of the nullipore reservoir, the Dykstra-Parsons coefficient (V) has been estimated (V = 0.91), indicating an extremely heterogeneous reservoir. A slight to high anisotropy (λk) has been assigned for the studied nullipore sequences. A stratigraphic modified Lorenz plot (SMLP) was applied to define the optimum number of HFUs and barriers/baffles in each of the studied wells. Integrating the permeability-porosity, reservoir quality index-normalized porosity index (RQI-NPI) and the RQI-flow zone indicator (RQI-FZI) plots, the discrete rock types (DRT) and the R35 techniques enable the discrimination of the reservoir sequences into 4 RRTs/HFUs. The RRT4 packstone samples are characterized by the best reservoir properties (moderate permeability anisotropy, with a good-to-fair reservoir quality index), whereas the RRT1 mudstone samples have the lowest flow and storage capacities, as well as the tightest reservoir quality.