3D Integrated Structural, Facies and Petrophysical Static Modeling Approach for Complex Sandstone Reservoirs: A Case Study from the Coniacian–Santonian Matulla Formation, July Oilfield, Gulf of Suez, Egypt

The Coniacian–Santonian Matulla Formation is one of the important reservoirs in the July oilfield, Gulf of Suez Basin. However, this formation is characterized by uncertainty due to the complexity of reservoir architecture, various lithologies, lateral facies variations and heterogeneous reservoir quality. These reservoir challenges, in turn, affect the effectiveness of further exploitation of this reservoir along the Gulf of Suez Basin. In this work, we conduct an integrated study using multidisciplinary datasets and techniques to determine the precise structural, petrophysical, and facies characteristics of the Matulla Formation and predict their complex geometry in 3D space. To complete this study, 30 2D seismic sections, five digital well logs, and core samples of 75 ft (ft = 0.3048 m) length were used to build 3D models for the Matulla reservoir. The 3D structural model shows strong lateral variation in thickness of the Matulla Formation with NW–SE, NE–SW and N–S fault directions. According to the 3D facies model, shale beds dominate the Matulla Formation, followed by sandstone, carbonate, and siltstone beds. The petrophysical model demonstrates the Matulla reservoir's ability to store and produce oil; its upper and lower zones have good quality reservoir, whereas its middle zone is a poor quality reservoir. The most promising areas for hydrocarbon accumulation and production via the Matulla reservoir are located in the central, southeast, and southwest sectors of the oilfield. In this approach, we combined multiple datasets and used the most likely parameters calibrated by core measurements to improve the reservoir modeling of the complex Matulla reservoir. In addition, we reduced many of the common uncertainties associated with the static modeling process, which can be applied elsewhere to gain better understanding of a complex reservoir.


INTRODUCTION
The current research developments within the energy sector are seeking better exploitation of subsurface resources through enhancement of resource prediction and efficient production. This can be achieved by accurate modeling of the rock physical and fluid characteristics in 3D space in order to improve production efficiency and resource recovery (Radwan, 2022;Wu & Xu, 2004, 2005Yan-Lin et al., 2011). Three-dimensional (3D) modeling of hydrocarbon reservoirs is the key factor that controls hydrocarbon reservoir evaluation and simulation for better exploitation and successful development. Although, 3D geological modeling can be done in oil and gas reservoirs easily using various available software, the modeling accuracy still presents a big challenge that has a great impact on the effective development of hydrocarbon reservoirs (Abdel-Fattah et al., 2010, 2018Bilodeau et al., 2002Bryant & Flint, 1993Radwan, 2022). Multi-source data can be integrated into the reservoir's 3D geological property models. The geological model allows attributes to be assigned to each cell, and the quantitative reservoir models are often cell-dependent. Such models place far greater challenges on the geologist than conventional models do because geology requires a complete explanation at any point within a reservoir's 3D volume. Multi-integrated properties can be used to define a 3D geological model, which is critical for connecting data from a wellbore to a 3D geological model in order to create efficient reservoir simulations Abu Amarah et al., 2019;Cosentino, 2005;Deutsch, 2002;Doyen, 2007;Qadri et al., 2019;Radwan, 2022). The main advantage of 3D modeling approaches is the ability to model complex reservoirs with variable lithologies and reservoir heterogeneities, and so accurate parameters and the best data integration are required for better reservoir modeling.
The Coniacian-Santonian deposits of the Matulla Formation are part of the Nezzazat Group, and they are distinguished by various lithologies, lateral facies and thickness changes, and reservoir heterogeneities across the Gulf of Suez basin ( Fig. 1) (Alsharhan, 2003;El-Sabbagh et al., 2011;Kassem et al., 2020;Shehata et al., 2021). In Gabal Nezzazat, Sinai, the Matulla Formation is a well-exposed siliciclastic/carbonate facies (El-Azabi & El-Araby, 2007;Said, 1990). The Matulla reservoir is the third most important reservoir in the July oilfield ( Fig. 1), following the Rudies and Nubia Sandstone reservoirs (EGPC, 1996). The structural and stratigraphic framework of the exposed Matulla Formation was established by El-Azabi and El-Araby (2007), Hasouba et al. (1992) and Marttila and El Bahr (1994). Due to reservoir heterogeneities and lateral facies changes in the Matulla sandstone reservoir, major exploration problems were encountered during its exploitation in the basin, posing a challenge and increasing exploration risks. In this regard, precise subsurface area modeling is required to solve the reservoir architecture in the July oilfield.
The ultimate goal of this study was to provide a new set of potentially optimistic reservoir targets for modeling hydrocarbon accumulation and reservoir characteristics in the complex and heterogeneous Matulla reservoir, using precise 3D modeling with low uncertainties. To overcome the associated facies changes and reservoir complexity of Matulla in the subsurface of the July oilfield, the reservoir architecture, petrophysical parameters, structural characteristics, and facies distribution must be understood and recognized. Simulation of complex reservoir property was used in this work to populate the grid cells of a 3D model with integrated constructed (structural modeling), discrete (facies modeling), and continuous properties (petrophysical modeling). Therefore, we utilized available data from 2D seismic reflection sections, wireline logs, and core sample analysis to define the main structural characteristics of the study area, recognize the distribution of the thickness/facies of the Matulla Formation reservoir, and evaluate its petrophysical perspective. Moreover, we applied best-fit parameters in the reservoir simulation to decrease final geological model uncertainties.

GEOLOGIC SETTING OF THE STUDY AREA
The Gulf of Suez rift is an extended graben shaped by Oligocene rifting as a northern continuation of the Red Sea (Garfunkel & Bartov, 1977;McClay et al., 1998;. The studied July oilfield is located in the Belayim province (NE-tilted fault blocks and SW-dipping normal faults), and it is bounded to the north by latitudes 28°13'N and 28°18'N and to the east by longitudes 33°11'E and 33°17'E ( Fig. 1) (EGPC, 1996;Gaafar et al., 2009). Day and Hoffman (1997) stated that the July oilfield is a highly faulted, complex structure formed by multiple reservoirs within 46 fault-blocks. Within the oilfield, compartmentalization results in different fluids, production, and depletion characteristics for each reservoir compartment.
An asymmetric NE-dipping tilted fault-block is connected to the SW by a SW-dipping fault, and characterized by large displacement normal faults that form the main part of the July oilfield. The major fault block is subdivided into northern and southern blocks, which have slightly different structural histories and trapping mechanisms, by an ENE-trending cross-fault. Initially, in the Oligocene, transverse movement on the cross-fault occurred, causing uplift of the southern block and erosion prior to Miocene deposition. Another transverse movement occurred in the Miocene (Upper Rudeis), causing uplift of the northern block and truncation of the Lower Rudeis reservoir sandstone on the crest (Brown, 1980). The Gulf of Suez stratigraphy can be divided into three main tectonic stratigraphic successions (Colletta et al., 1988). Some authors describe these successions as lithostratigraphic units (Alsharhan, 2003;Said, 1990). They include pre-rift (pre-Miocene) interval (e.g., Radwan et al., 2020b, syn-rift (Oligocene-Miocene) interval (e.g., Kassem et al., 2021aKassem et al., , 2021bRadwan & Sen, 2021a, 2021bShehata et al., 2021), and post-rift (post-Miocene) interval (e.g., Darwish and El Araby, 1993;Radwan, 2020a;2021a, 2021b. The thickness, areal distribution, lithology, hydrocarbon importance, and depositional environments of these units vary (e.g., Abdelghany et al., 2021;Alsharhan, 2003;Elkhawaga et al., 2021;Radwan et al., 2019aRadwan et al., , 2019bRadwan et al., , 2019c. In the July oilfield, the sedimentary sequence ranges in age from Paleozoic to Recent (EGPC, 1996). The general stratigraphic section of the July oilfield is shown in Figure 2. According to EGPC (1996), the main reservoirs of the July oilfield are Miocene clastic deposits (Rudeis and Nukhul Formations), pre-Miocene Cretaceous Nezzazat Group, and Nubia Formation. The studied deposits of the Matulla Formation were alternating shallow marine shales, carbonates and sandstones of Upper Cretaceous age (Coniacian-Santonian). The Coniacian-Santonian Matulla Formation is well formed in the southern and central areas of the basin, and it can be divided into three depositional units (M1-M3) (El Sharawy & Nabawy, 2015). An unconformity separates the intermediate unit (M2) and the upper unit (M3) (El Sharawy & Nabawy, 2015).

MATERIALS AND METHODS
The Gulf of Suez Petroleum Company (GUP-CO) provided the datasets for the current study as part of an agreement with the Egyptian General Petroleum Corporation (EGPC). Thirty 2D seismic sections were represented by 15 NE-SW direction sections (in-line seismic), 11 NW-SE direction sections (cross-line seismic) and four N-S direction sections (randomly) (Fig. 1). Digital well logs for the July-10, 20A, 29A, 82, and 93 wells were also provided, which include caliper, gamma ray, resistivity, PEF, neutron, density, and sonic datasets. The July-93 well was the only well that was cored within the Matulla Formation section in the central area of the July oilfield. The flowchart of the processing steps used in this work is provided in Figure 3.
During the reservoir 3D modeling process, multi-simulation methods were proposed to understand and model the reservoir. The sequential Gaussian simulation (SGS) method is used widely in reservoir modeling applications (i.e., Petrelä software) (Hu & Le Ravalec-Dupin, 2005;Isaaks, 1991;Radwan, 2022). The ability to simulate complex structures is the most significant advantage of 3D modeling techniques. The modeling process is divided into several stages beginning with seismic, structural, facies, and petrophysical modeling, then ending with the final interpreted model as the following.

Seismic Interpretation
The interpretation of seismic data is an integral part of the hydrocarbon development strategy of the oilfield. It performs a variety of tasks in 3D seismic interpretation, such as analyzing the parameters that govern reservoir architecture (John et al., 2010;Qadri et al., 2017Qadri et al., , 2020. The seismic interpretation, which was based primarily on seismic data provided by the company, included the following.

Performing Seismic to Well Tie
The seismic data are estimated in time domain, and are compared with the well data, which are measured in depth. In this regard, the primary goal of well-to-seismic ties is the correlation of horizon tops that are found in a well with specific seismic reflections. In this study, the data were used to tie well data to seismic are simple because the seismic in depth domain facilitated the connection between wells and seismic lines.

Selecting Studied Horizons on Seismic Sections
Seismic reflection data from every in-line and cross line were interpreted to overcome the data quality issue and retain the precision of horizon and fault interpretation. The interpreted seismic section in different directions is passed near to all wells. All lines intersect in the seismic grid. In this study, we focused on the Matulla Formation, which is the study target.

Interpretation of Faults Affecting These Horizons
Fault movements affect hydrocarbon traps and pathways in addition to deforming sedimentary basins (Gartrell et al., 2004;. Thus, fault interpretation is an important part of a structural framework model. Based on the discontinuity in reflection and the abrupt termination of reflection occurrences, the faults present in the study area were established.

Constructing Depth Structure Contour Maps
In the present study, depth structure contour maps were constructed for the Matulla Formation by using Petrel software to illustrate the subsurface structural configuration of the investigated area. The interpreted seismic data were the main container for the 3D geologic modeling (facies and their petrophysical characteristics) according to Radwan (2022). During the 3D structural modeling, the seismic interpretations were the main input data for the Matulla Formation.

Structural Modeling
The second stage of seismic data analysis was the creation of a structural model. The fundamental benefit of a 3D model is that it allows the interpreter to observe and evaluate structural data by displaying a cross section across the model in any direction (Fagin, 1991;Radwan, 2022). A structural model of a reservoir provides borders for the facies and petrophysical models that define the rock qualities, forming the geometrical framework for the 3D grid. The development of a structural model not only allows us to test the model through structural sections, but it also gives us a vision for new well trajectories Qadri et al., 2017;Radwan, 2022), volumetric calculations, and reservoir simulation grids (Holden et al., 2003).
Three major processes were required to complete the structural modeling (i.e., fault modeling, pillar gridding, and horizon generation) Radwan, 2022). The constructed faults were modeled first, whereby this process defined the faults in the geological model and formed the basis for generating the 3D grid, and defined the shape of the faults being modeled by generating the key pillars. Before going to the next step of structural modeling, the key pillars needed some adjustments. A modeled fault must never cross another fault without being connected; however, there were a few cases of crossing faults. The key pillars of the crossing faults were connected in the 3D window, then the pillar gridding process, which involved creating structural grids from the fault model, was carried out. This procedure created the model's skeleton framework. The skeleton grid was made up of three sections: a top, a middle, and a base skeleton grid. Each grid was tied to the essential pillars top, middle, and bottom points. The pillar b Figure 2. Generalized stratigraphic column of the July oilfield (Alsharhan & Salah, 1995;EGPC 1996).
grid parameters of the study area included a 3D grid increment of 50 m by 50 m and a gridding boundary area limited to depth surfaces. The final step was to create horizons boundaries that denoted the adding of the model's various horizons, which were represented through the Matulla reservoir boundaries (top and bottom) (Abdel-fattah et al., 2010; Abdelmaksoud, 2017; Ali et al., 2019b;Avseth et al., 2005;Cannon, 2018;Cosentino, 2005;Radwan, 2022).

Petrography and Facies Analysis
The petrography is important in well logging because it provides evidence of the mineralogical composition that can affect log reaction (Abudeif et al., 2016a(Abudeif et al., , 2016b(Abudeif et al., , 2018Radwan et al., 2020c;Mutebi et al., 2021). Polarizing microscopy of thin sections was used for petrographic analysis of the rock samples. The lithological interpretation of the reservoir was based on the study of thin sections that were prepared from core samples. The cores samples cover the depth interval of 11,413-11,488 ft 1 in the July-93 well. These data were integrated to evaluate the Matulla Formation and to delineate its facies distribution. The various lithologies/facies properties were used as input in the created 3D structural system via facies modeling processes to obtain a 3D facies model for the Matulla reservoir.

Petrophysical Analysis
The petrophysical evaluation was performed through study of the various well logging curves (Fig. 4). The goal of petrophysical analysis was to determine the key parameters that were needed for reservoir characterization and assessment of hydrocarbon potentiality. In this study, wireline logging, composite logs, and internal geological reports were used to evaluate and delineate the reservoir petrophysical properties. The critical parameters of the Matulla Formation (i.e., effective porosities, shale volume, water and hydrocarbon saturations) were determined using industry standard techniques. The water saturation estimation in the contaminated Matulla sandstone with varied percentage of clays across the reservoir was done using the Indonesian equation of Poupon and Leveaux, (1971). Subtracting Sw from 100 percent saturation yields the hydrocarbon saturation (Shc). Water resistivity (Rw) equals 0.013 X m according to Pickett plots (i.e., resistivity-porosity relationship).
The net reservoir and pay thicknesses were calculated using standard cut-off, where the percentages of porosity, shale volume, and water saturation were 10%, 35%, and 50%, respectively. The results of the petrophysical analysis were used as inputs into the created 3D structural system via petrophysical modeling to obtain a 3D petrophysical model for the reservoir.

Facies Modeling
The basic goal of reservoir description was to identify reservoir heterogeneities that affect the amount, position, accessibility, and movement of fluids across the reservoir. Because standard data collecting methods produce sparse sampling, 3D description at the appropriate resolution was generally problematic. In order to ÔÔpredict'' geological, geophysical, and petrophysical properties between sampling points (e.g., wells and seismic lines) and define probable variabilities, geostatistical approaches were used (Radwan, 2022;Seifert & Jensen, 1999).
To build the 3D facies model, the allocated values were loaded into Petrelä software in a distinct code for the lithology distribution, and then dispersed in two directions (vertically and horizontally) across the grid cells to fill the full 3D grid between the wells (Abdel-Fattah et al., 2018; God- will & Waburoko, 2016; Jika et al., 2020;Okoli et al., 2021;Radwan, 2022). In addition, the 3D facies model was applied using the stochastic sequential indicator simulation (SIS) method (Deutsch & Journel, 1998;Journel, 1982;Radwan, 2022;Remy et al., 2009). This algorithm is the most widely used for discrete or categorical variable data such as facies, whereby upscaled facies values and allocated variograms are the most important factors (Deutsch & Journel, 1998;Radwan, 2022).
When using stochastic modeling, SIS, a number of uncertainties can arise. In this simulation approach, a grid-node that has no simulated value (an unknown lithotype) is selected randomly and neighboring node points are identified and simulated with reference to the original conditioning dataset.
Later, another grid-node is selected randomly, and the variable is simulated using the newly generated CDF (cumulative distribution function), which is now increased in size by one value. In this way, each node is simulated (Deutsch & Journel, 1998;Seifert & Jensen, 1999;Okoli et al., 202;Radwan, 2022).

Petrophysical Modeling
The main petrophysical parameters models were constructed based on the IPTM software's interpretation of the results (Cannon, 2018). Because the available data for deterministic modeling were insufficient, a stochastic method was employed instead. Only the vertical well profiles were used as bases for the 3D spreading of the petrophysical properties. For the Matulla reservoir petrophysical parameters, the normal score transformation was used. As a statistical approach that works with the amount of available data, namely the SGS algo- rithm, was used (Pyrcz & Deutsch, 2014;Radwan, 2022). The major inputs for the petrophysical modeling included upscaling of well logs, input distribution, and variogram construction. The vertical range was set to twice the zone's thickness. To ensure that the transformation applied to the data was reasonable (valid), a statistical quality check was done between both of the upscaled and modeled petrophysical attributes to ensure that the upscaled and modeled petrophysical characteristics' standard deviations and variances were nearly identical Radwan, 2022).

Depth Structure Contour Map of the Matulla Formation
After completing the seismic interpretation process for all lines, the seismic to well tie was performed. The studied horizons were then selected on seismic sections, and the faults affecting them were interpreted (Fig. 5a, b, c). The preceding steps made it simple to create depth structure contour map of the Matulla Formation. The reflection map (Fig. 6) at the top of the Matulla Formation shows strong spatial variability of the horizon in the study area, with a general deepening to the southeast and northwest reaching maximum value of À 15,000 ft TVD SS (true vertical depth subsea), at seismic shotpoint (S.P.) 181 on seismic line 2072. In addition, the seismic line 2134 at S.P. 41 showed topog-  raphy shallowing to the west and southwest sections of the region, with minimum value of À 8200 ft TVD SS . The main faults delineated on seismic sections showed a clear NW-SE-trending fault (Red Sea System), which primarily dips southwest with some faults in the southern part dipping northeast and forming horsts. The other major NE-SW and N-S faults are thrown northwestward. As shown in Figure 5, the Matulla Formation structures contain faults and horsts structural closures, which were also noticeable in the majority of the study area.

Isopach Map of the Matulla Formation
The Matulla Formation reservoir's isopach map was created by subtracting the Wata interval's depth map from the Matulla interval's depth map (see stratigraphic position in Fig. 2). All of the wells investigated (July-10, July-20A, July-29A, July-82, and July-93) penetrated the Matulla Formation in this region. The Matulla Formation isopach map (Fig. 7) shows that the thickness of the Matulla Formation increases in the southeast and west sections of the July oilfield. However, the thickness decreases in the central and north eastern sections of the July oilfield. The basinal area is directed to the southeast and west parts of the July oilfield and has a maximum thickness of 1100 ft, while the central and northeast parts of the July oilfield are on the platform and have minimum thickness of less than 100 ft.

3D Structural Model
A precise structural model of the Matulla Formation was required for the explanation and knowledge of the structural pattern of the reservoir in a 3D outline. By optimizing the information from seismic data, the goal was to create a solid structural model using the available wells (Fig. 8). The July oilfield is a highly faulted, complex system that generates 46 fault blocks from four reservoirs (Day & Hoffman, 1997), some of which are isolated and have different fluid properties as well as performance and depletion characteristics (Baack & Latif, 1995). The 3D structural framework illustrates complex structures that affect the July area, with normal faulting serving as the primary faulting system. According to the 3D structural model, normal faults with NW-SE trends represent the dominant fault system in the July oilfield. At the uppermost sandstone of the Matulla Formation, the crest of the structure occurs at $ 9000 ft TVD SS in the central area of the main part of the oilfield (present study), whereas the OWC is at À 10,800 ft TVD SS , giving an oil column of $ 500 ft (Marttila & El Bahr, 1994). Moreover, the model indicated that the horizon depths of the Matulla Formation reservoir were greatest in the eastern and southern east parts of the July oilfield and decreased gradually toward the center and margins of the basin.

Petrographic Analysis and Lithological Features
Evaluation of sedimentological and well log analysis is the key for better reservoir discrimination and better reservoir quality prediction (Radwan, 2021a;Qadri et al., 2021). The Matulla Formation petrographic studies and microfacies analysis were carried out on samples from well July-93, based on the cumulative core interval from 11,413 ft to 11,488 ft (75 ft of total thickness). The stratigraphic section of the Matulla reservoir was composed of fine-to very fine-grained clastic sediments, according to thin section sample analysis, core samples, and well logging analysis. The dominant lithologies were sandstones, shales and limestones accompanied by minor siltstone. The Matulla shale had a dark gray color, and the shale cuttings had moderate hardness and sub-blocky shape. Moreover, the Matulla shales were non-calcareous to slightly calcareous. In terms of hardness and crystallinity, the Matulla limestones were moderately hard and cryptocrystalline. The Matulla limestones were, in contrast, highly argillaceous, indicating the origin of these limestones gray colors. The petrographic characterization of the Matulla sandstones is discussed below (Fig. 9).

Grain Framework
According to thin section analysis, quartz grains were the most abundant mineral constituents of the Matulla sandstone, accounting for more than 70% of the individual samples. Many samples contained monocrystalline and polycrystalline quartz grains. orthoclase accounting for about 5% of the feldspar grains. Accessory minerals averaged between 1 and 2%, and included glauconite, zircon, pyrite and mica. Glauconite may reach a higher percentage in some samples, forming green sandstone in this case. According to the sandstone classification, the Matulla sandstone can be classified sub-lithic arenite sandstone (Pettijohn & Potter, 2012). Variable porosities, primarily of the inter-particle and intraparticle types, were found between the grains, limiting the ability of the reservoir to store oil.

Cement and Matrix
The upper part of the core interval was cemented with calcite, dolomite rhombs, and glauconitic grains, while the lower part was cemented with argillaceous, dolomite, and glauconitic materi-als. However, silica cement is uncommon. Silt-size matrix is dominant in the Matulla reservoir.

Texture and Fabric
Quartz grains were sub-angular to sub-rounded, with the latter being more prevalent in the succession. The Matulla sandstone reservoirs in the July field had grain sizes ranging from fine to very fine sand. In most samples, the well sorted texture was dominant; however, some samples were moderately sorted.

Diagenetic Features
Many processes have reduced the reservoir potential of different facies types, according to diagenesis studies of the Matulla sandstones. In decreasing order of importance, detrital clays, carbonate cementation (calcite and dolomite), compaction, and kaolinite formation were observed. The pore system and reservoir potential both increased due to two factors occurring during diagenesis: dissolution and fracturing. Some of the samples exhibited minor dissolution of detrital grains (primarily feldspars) and partial dissolution of carbonate cements and clays.

Petrophysical Parameters of the Matulla Formation
The petrophysical and hydrocarbon characteristics of the Matulla Formation were based mainly on well logging analysis of the five studied wells, which were distributed across the study area. The average values of the various derived petrophysical parameters of the examined wells of the July oilfield are summarized in Table 1. The petrophysical data logs of selected wells from the July oilfield are depicted in Figures 10 and 11.
The Matulla Formation in the July oilfield was encountered in the depth range 9766-10,300 ft in July 10, 10,665-11,760 ft in July 20A, 9350-10,690 ft in July 29A, 11,086-11,462 ft in July 82, and 11,027-11,611 ft in July 93. The lateral and vertical changes in petrophysical properties were thought to be directly connected to lateral and vertical facies changes across the Matulla reservoir. The petrophysical analyses of the wells under consideration showed good reservoir properties and hydrocarbon fluid entrapment in the Matulla Formation. The lithosaturation cross-plot revealed a complex lithologic nature that varied between sandstone, carbonate and shale with silt intercalations at different levels. Relatively low GR, high resistivity, high density, and low neutron porosity logs were recorded against carbonate lithology intercalated with shale zones. Along the shale zones, high GR, low resistivity, low density, and high neutron porosity values were observed. However, the sandstone had low gamma ray readings, porosity ranging from low to intermediate, To distinguish the entire Matulla Formation, it was divided into three zones: upper, middle, and lower (see well correlations in Fig. 4). Pore fluid discrimination is critical in petrophysical evaluation (Agbasi et al., 2018;Radwan et al., 2020c). Oil was the fluid type in the Matulla Formation, and the majority of the source rocks in the July oilfield were oil prone (EGPC, 1996;Day & Hoffman, 1997). Figure 12a and 12b display the oil water contact at depth of À 10,600 ft TVD SS , providing a large column of oil in the Matulla reservoir. According to the petrophysical investigation, the Matulla sandstone had excellent reservoir characteristics and that it can store and transfer hydrocarbons. The statistical analyses showed a good range of total and effective porosity, as well as permeability, implying that the Matulla Formation's reservoir quality was good. The extreme lateral and vertical changes in lithology may be linked to the significant lateral and vertical changes in petrophysical properties of the Matulla reservoir.

3D Facies Model
The initial stage in facies modeling is to create discrete facies that are upscaled (Fig. 13) to gridbased model cells (Abd El-Gawad et al., 2019;Radwan, 2022;Tyson, 2007). In terms of lithology, the Matulla Formation displayed various lithologies due to the facies changes from sandstone and shales to carbonate facies. In comparison to sandstone, the shale fraction had a larger proportion of the probability distribution, followed by limestone and siltstone in decreasing order of abundance. The outcomes of the data analysis were employed for facies modeling. Based on the combined sedimentological and petrophysical analysis, the following were the proportions of the various lithologies used: 23% sandstone, 9.6% siltstone, 39.3% shale, 21.7% limestone, 3.8% calcareous shale, and 2.6% dolomite (Fig.14). The interpreted facies of mixed carbonate/siliciclastic facies of the Matulla Formation consisted of different lithologies that have been deposited from lower intertidal to slightly deep subtidal in a fluctuating depositional regime. El-Azabi and El-Araby (2007) interpreted various types of foreshore sandstone environments and shallow subtidal siltstone for the Matulla Formation; however, they also described deep subtidal shale/claystone and shallow water lime-mudstone. Cross sections of the Matulla Formation reservoir were generated in addition to the facies model (Fig. 15) to highlight the lateral and vertical facies changes (Fig. 16). Facies variation from continental fluvial to marine facies with tidal flood plain conditions in between reflects sea level fluctuations through the Matulla Formation (El-Azabi & El-Araby, 2007). The rapid facies changes across the studied wells reflect the changes in depositional conditions.

3D Petrophysical Model
The interpreted petrophysical characteristics of the Interactive Petrophysics Program were upscaled and then modeled using standard procedures. In the process of petrophysical modeling, the SGS algorithm was the utilized statistical approach (Ali et al., 2019a(Ali et al., , 2019bRadwan, 2022). Given the amount of available data, this approach is most likely the best fit in this case, where it had the advantages of simulating continuous variables such as permeability and porosity (data after being normalized) (Fig. 14).
The parameter values of unsampled intervals were determined using the equations derived from data regression during SGS, and then all values were forced to satisfy the normal distribution by normal scoring. The fundamental simulation parameters are illustrated in Figure 17. In order to characterize the 3D distribution of these specific attributes within the July oilfield, several cross sections were derived from these models (Fig. 18).
The  relatively high values of water saturation and shale volume, which were decreasing in the central and southeast sectors. These results explain the increase in effective porosity and permeability in the lower zone, with slight increase in the north and central directions. This lower zone had relatively low water saturation and shale volume. We deduced that the upper zone had overall good to very good reservoir quality (high porosity, low shale volume, and low water saturations), especially in the middle interval; also, the lower zone  was of good quality reservoir. However, the middle zone had poor reservoir quality (low porosity, high shale volume, and high water saturations). Finally, according to the integrated 3d model, the most favorable sectors for hydrocarbon accumulation and production through the Matulla reservoir were located in the central, southeast, and southwest sectors of the July oilfield. Therefore, it is recommended that future development activities be concentrated on these areas.

Volumetric Calculations
For a realistic resource estimation, volumetric calculations were employed in conjunction with geological reservoir models (Attia et al., 2015;Cannon, 2018;El Khadragy et al., 2017;Radwan et al., , 2022. We used the following formula to compute the volume of hydrocarbons (Dean, 2008): where OOIP is original oil in place, A is study area in acres, h is net pay thickness in feet, / is porosity in %, SW is water saturation in % and Bo is the oil formation volume factor. The volumetric estimations of oil within the Matulla reservoir were performed using a 3D grid-based geocellular model developed in this work ( Table 2). The OIIP at the July oilfield was assessed to be 715 million stock tank barrels (MSTB).

Model Uncertainties
Deriving of rock and petrophysical parameters from well logs and cores requires some caution (Radwan et al., 2020c). As with all geological interpretation processes in the energy sector, there are uncertainties that can accompany geological, seismic, and petrophysical interpretations. Uncertainties in interpretation have an impact on resource evaluation and modeling. Uncertainties plague reservoir modeling frequently due to restricted data availability, data quality, difficulty of reservoir design in predicting flow in a specific reservoir, and modeling process uncertainties (Radwan, 2022;Singh et al., 2009). Therefore, it is critical to examine the input data for quality, quantity, and complexities at various scales, as well as to revisit basic assumptions of a modeling workflow's appropriateness in the light of static reservoir uncertainties.
Moreover, the software required integrated resources in order to reposition the base case model with a reasonable degree of certainty for more accurate project assessment, risk analysis, and decision-making. Certain techniques have recently been developed to solve these issues. While great progress has been made in the areas of gridding and upscaling in recent years, significant practical obstacles in reservoir modeling still remain. In the not-too-distant future, effective and accurate depiction of reservoir shape and heterogeneity will come first on the list. This is because, as technology improves, the amount of data acquired from producing oilfields grows, geology understanding improves, and reservoir descriptions are becoming more detailed. Therefore, achieving a balance between explicit representation of reservoir geometry and heterogeneity by gridding and correct capture of their impacts on multi-phase flows via upscaling is often an ongoing challenge in reservoir modeling.
We attempted to perform precise interpretations and extract the best fit input parameters in this work in order to achieve a more precise static model for the Matulla Formation in the July oilfield. Determination, organization, and applications of structural uncertainty describe the key sources of uncertainty in determining seismic quality and its restrictions for improved reservoir characterization. General geologic expertise, regional geology, and some previous work skill from neighboring oilfields were used to fill in the gaps. The reservoir geological information and restricted data are the key sources of uncertainty in lithofacies modeling. As a result, the spatial distribution of facies or flow units within the three-dimensional structural system is highly uncertain. In this work, we tried to reduce uncertainties in the facies model; therefore, we applied integrated facies analysis, petrographic description, and well logging interpretation across the studied wells. To achieve the best lithofacies interpretation along the whole area, we used SIS.
By using calibrated well-log data and core measurements, the petrophysical models were simulated to estimate the level of uncertainty in the petrophysical parameter distribution. Stochastic modeling is the process of creating synthetic geologic architectures and/or petrophysical property distributions based on quantitative and qualitative data. The availability of core data from the entire Matulla Formation allows for precise determination of the petrophysical parameters as critical input parameters. Volumetrics are derived from a number of reservoir factors that are often interpreted or predicted based on limited data and were subject to uncertainty. Bulk, pore, hydrocarbon pore volume, and OIIP are volumetric measurements used in a hydrocarbon resource appraisal. Uncertainty studies of these volumetric variables necessitate a thorough examination of data uncertainties in associated input variables as well as implication uncertainties when covering restricted data over the entire oilfield. In this work, we used the most likely fit parameters, which were calibrated using core data and production data. We provided a summary of reservoir modeling issues, focusing on the complex geometry of the complex Matulla reservoir and diverse reservoir parameters based on the available data and the most likely extracted reservoir parameters. Adding more data from more exploratory and development wells targeting the Matulla reservoir can result in more refined and highly accurate reservoir models that accurately represent the subsurface Matulla reservoir.

CONCLUSIONS
The present work dealt mainly with the interpretation of geological and geophysical data to evaluate the structure, lithofacies and petrophysical characters as well as the volumetric parameters of the Matulla reservoir and their 3D distribution in the July oilfield. In this study, we used integrated datasets to assess the distribution and quality of the complex Matulla reservoir, as well as to reduce uncertainty during oilfield development in July. The following were the main findings of this study.
The contour lines of depth-structure contour maps showed the spatial variability of the top and base of the Matulla reservoir across the July oilfield with a general deep topography to the southeast and northwest reaching maximum value of 15,000 ft. However, the Matulla contour lines displayed shallow topography to the west and southwest sections of the region, which was mapped to minimum value of 8200 ft. The main faults that were delineated on the seismic sections showed a clear NW-SE-trending fault (Red Sea System), which dips mainly to the SW, with some NE-dipping faults in the southern part forming horsts. The other major NE-SW and N-S fault patterns are thrown northwest ward. The isopach map showed that in the southeast and west sections the Matulla Formation, thickness increased and reached maximum thickness of 1100 ft, and decreased in the central and northeastern portions of the region platform reaching minimum thickness of less than 100 ft. Petrographic investigations of accessible core samples revealed that the Matulla sandstone reservoir samples were made up mostly of colorless-light gray sublithic arenite sandstone, which was medium-to fine-grained, sub-angular to sub-rounded, and moderately sorted. The Matulla sandstone is cemented by calcite, dolomite rhombs, argillaceous, and glauconitic minerals. The sandstone facies displayed good visual porosity, which reflects the Matulla sandstone's ability to store oil. The 3D facies model indicates that the Matulla Formation is made up mainly of 23% sandstone, 9.6% siltstone, 39.3% shale, 21.7% limestone, 2.6% dolomite and 3.8% calcareous shale. According to the petrophysical model, the Matulla Formation's upper and lower zones are good quality reservoir, while the middle zone is a poor quality reservoir. It was also deduced that lateral and vertical changes in petrophysical properties were linked to lateral and vertical changes in lithology across the complex Matulla reservoir. To better exploit the Matulla reservoir, the most favorable sectors in the Matulla reser-voir for accumulation and production of hydrocarbons are likely at the central, southeast and southwest parts of the studied oilfield. We provided an overview of reservoir modeling issues, focusing on the complex geometry of the complex Matulla reservoir and the diverse reservoir parameters based on available data and the most likely extracted reservoir parameters. We combined multiple datasets and used the most likely parameters calibrated by core measurements to improve modeling of the complex Matulla reservoir and reduced the common uncertainties associated with the static modeling process.

OPEN ACCESS
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecom mons.org/licenses/by/4.0/.