Delineating the main structural outlines and the petrophysical properties of the Albian-upper Cretaceous Reservoirs using seismic and well log data, Shushan Basin, north Western Desert, Egypt

This study concerns with delineating and analyzing the subsurface structural setting of Khalda Oil Field at Shushan Basin as a key factor for evaluating the petrophysical properties and the deliverability of the Bahariya Formation. Shushan Basin is considered one of the most prospective basins in the Western Desert. The oil is trapped within a gentle seismically defined anticline that contains a series of coastal sandstone reservoirs interbedded with a neritic and tidal flat mudstone and a few carbonate interbeds of the Albian-Cenomanian Bahariya Formation forming multiple oil-bearing sandstone reservoirs. A set of borehole logging and seismic data were available, tied, processed, and mapped to delineate the predominant structures in the study area. Results show the implication of the NW–SE, NE–SW, and E–W gravity faulting on the Bahariya and the overlying Abu Roash formations. This set of gravity faults is attributed to the Mesozoic tensional stresses. The depth structure contour maps on the upper surfaces of the Lower and Upper Bahariya members reveal the presence of some NE–SW gravity faults that dissected the Khalda Field into some slightly tilted blocks, shallowing in depth to the south as deepening to the north. Besides, the study area has been affected later in the Jurassic by E–W gravity faults during the Atlantic Ocean opening. This faulting trend was changed later into the NE–SW direction as a result of the movement of North Africa against Laurasia. The petrophysical analysis indicated that the Lower Bahariya sandstone reservoir has a good reservoir quality (16 ≤ ∅ ≤ 25%, 53 ≤ Sw ≤ 59%, 6.0 ≤ Vsh ≤ 32.2%, and good net-pay thickness 18.0–38.0 ft). It is revealed that the drilled wells penetrating the crest of the anticline are prospective while that located in the anticline flanks are water-bearing recommending exploration away from the anticline flanks. By applying this workflow, it is possible to explore for the similar subsurface hydrocarbons-bearing sequences in the Western Desert and North Africa in future exploitation plans.


Introduction
The Sedimentary basins of the Western Desert (WD) in Egypt are among the most prolific areas for hydrocarbon exploitation in North Africa. The WD is bounded to the north by the Mediterranean shorelines and extends for 600 km from the Libyan borders in the west to the Nile Valley in the east. The sedimentary basins in the north WD are aged as Mesozoic and Early Tertiary basins, which have been variously described as the Northern Basins or the Coastal Basins. They lie to the north of Qattara Ridge which is known as Sharib-Sheila High which separates these northern basins from the Abu Gharadig Basin in the center of the WD (Fig. 1). These Coastal Basins comprise four discrete Mesozoic basins: Matruh and Shushan Basins to the west, while to east the Dahab/Mireir and Natrun basins are situated (Fig. 1). The hydrocarbons have been trapped in these basins within structural traps that follow the Syrian Arc structural trend. These traps are formed as fault-blocks, and three to four way closure structures (Teama and Nabawy 2016;Shehata et al. 2020Shehata et al. , 2021Nabawy et al. 2022).
The Shushan Basin contains a northward thickening Mesozoic and Tertiary section, which unconformably overlies the thick Paleozoic sediments of Siwa Basin to the west near the Libyan borders (EGPC 1992). It contains two prominent Jurassic to Lower Cretaceous grabens, namely, the NE-SW Shushan Trough, and the triangularshaped Matruh Trough, which is defined by a combination of NE-SW and NNE-SSW faults.
The Cretaceous Bahariya Formation with its dominant clastic composition and regional extension is considered among the most prolific formations in the coastal Basins of the north Western Desert, and in particular the Shushan Basin. Also, it is considered among the highly heterogeneous lithologies in Egypt; it is composed of intermixed medium to coarse-grained sandstones and shales with a few carbonate interbeds (Abu Shady et al. 2010). Therefore, it has been intensively studied in the different coastal basins, but it is still in need to more detailed studies in Shushan Basin and particularly in in Khalda Oil Field (Hassan et al. 2015;Nabawy et al. 2022). However, reservoir characterization based on integrated seismic and petrophysical data was widely by many authors in different areas in all over the world (Emujakporue et al. 2012;Aminzadeh 2012, 2015;Mohamed et al. 2016;Wang et al. 2017;Abdelmaksoud, and Radwan 2022). Thereby, a detailed structural study, based on two-way-time seismic data, is of critical importance to support the future development and exploration plans in the Northern Basins in NW Africa.  (Rossi et al. 2001)

Lithostratigraphy
The Lower Cretaceous lithostratigraphic formations in the studied sequence are represented by Abu Roash, Bahariya, Kharita, Alamein, and Alam El Bueib formations (Fig. 2). Alam El Bueib Formation consists of siltstone, shale, sandstone, dolostone, and limestone intercalations. The sandstone interbeds are prospective as hydrocarbons-bearing zones, whereas the shale beds are considered source rocks. The Aptian Alamein Formation consists of dolostones with some shale and limestone interbeds that have been deposited in a marine transgression (Schlumberger 1995). Upwards, the Albian Kharita Formation consists primarily of fluvial medium to coarse-grained sandstone and friable sand intercalated with thin siltstone and shale layers. It is overlain by the Albian-Cenomanian Bahariya Formation.
The Bahariya Formation, the core of the present study, in Khalda Oil Field is primarily hydrocarbon-bearing sandstone sequences. Upwards the Abu Roash Formation is built of interbedded shale and limestone, intercalated with some thin siltstone, sandstone, and anhydrite interbeds. The Bahariya Formation is the most promising oilbearing reservoirs in Khalda Oil Field. It is primarily composed of sandstone facies accumulated in a transgressive shallow marine coastal depositional setting, with some neritic shales and thin carbonate interbeds. The Bahariya Formation is divided into upper and lower members (Fig. 2) and reaches up to 1013 ft thick in Khalda-1 well (EGPC 1992). Thickness of the Upper Bahariya Member is several hundred feet and is mainly composed of shale and thin beds of limestone, whereas the Lower Bahariya Member is thick and composed of sandstone and shale interbeds. The Upper Bahariya Member in the North Western Desert has been studied sedimentologically in detail by many authors (e.g., Almoqaddam et al. 2018;Kassab et al. 2019;Abdelmaksoud et al. 2019;Nabawy et al. 2022 (Schlumberger 1995) acts as a seal rock that effectively traps the oil in the sandstones of the Lower Bahariya. It overlies unconformably the Middle Jurassic Khatatba Formation which consists of three depositional units: upper shale/siltstone, middle sandstone, and lower shale unit. The lower and upper shale units are considered the prime source rocks for producing the trapped oil in Khalda structure (EGPC 1992).

Structure setting
The Khalda Oil Field occurs at the northern end of the Shushan Trough which subsided rapidly during the Jurassic and Early Cretaceous. The Shushan Trough is mostly defined by a set of NE-SW faults which is mostly recognized on top of the Upper Jurassic Masajid horizon (Fig. 2), and less commonly on top of the Aptian-Alamein dolomite. These rift-bounding faults were apparently sealed shortly after the deposition of the Alamein Formation, so the structural map on top of the Cenomanian Bahariya reservoir does not reveal the presence of these deep-seated grabens (Fig. 3). The detailed morphology of the Shushan Trough is poorly known because seismic data quality below 2.5 s is poor (Abdou 1998). However, the structural framework of the Late Cretaceous-Tertiary sequence is dominated by the ESE-WNW or E-W dextral strike-slip faults, and the SE-NW gravity faults that have been formed during the Upper Cretaceous-Eocene Syrian Arc tectonic event (Fig. 3).
The setting of Khalda structure may have been controlled by deep-seated NE-SW normal faults associated with the Shushan Trough (Fig. 4), which were subjected to NW-SE compression during the 'Syrian Arc' tectonic phase. This same compressive phase is responsible for the anticlinal nature of the field (Fig. 5).
The Khalda structure on top of the Bahariya reservoir is represented by a gentle slope anticline oriented NE-SW of 5.5 km length and 2 km width. The remarkably rectilinear east flank has a pronounced NNE-SSW trend, and is steeper (up to 60°) than the irregular W to NW flank (dips up to 30°). The structure is cut by some normal SE-NW faults and WNW-ESE dextral strike-slip faults (Sehim et al. 1992) (Fig. 6). The faults have thrown the cut beds for more than 50 ft with some sealing implication which is noticeably situated at the southern borders of the Khalda structure (EGPC 1992). The highest point on the structure occurs in the south of the anticline, in the vicinity of Khalda-2 and -4 wells, and is apparently associated with a significant NW-SE cross-fault.  (Sehim et al. 1992)

Data and Methods
In this section, the available borehole logging data and the seismic lines will be listed and the applied methodology and techniques will be also discussed in some details.

Petrophysical data
The formation evaluation and reservoir characterization is based on the petrophysical analysis of this formation (Shehata et al. 2021;Nabawy et al. 2022;Haque et al. 2022;Fallah-Bagtash et al. 2022). Consequently, for the present study, a complete set of borehole logs were released for five wells in Khalda Field (Kh_45, Kh_56, Kh_58, Kh_21x, and Kh-S-1 × wells, Fig. 7) crossing the Bahariya Formation including resistivity logs, neutron, density, PEF, and gamma-ray logs. The petrophysical analysis includes lithologic definition based on a set of X-Y cross-plots for the petrophysical parameters, and the litho-saturation vertical plot, and the contour maps of the petrophysical parameters in the different wells (Baouche and Nabawy 2021;Abdeen et al. 2021). In this regard, the available well log data of the studied wells have been evaluated using the IP3.6 (Interactive Petrophysics); it is a proficient software applied to process and present the petrophysical well logging data as a function of depth and as a set of X-Y plots. The borehole and the environmental conditions were considered during the processing these petrophysical data using the Schlumberger charts (Schlumberger 1995;Nabawy and Shehata 2015;Radwan et al. 2021Radwan et al. , 2022a.

Lithology determination
Different methods were applied simultaneously to determine the lithology including the density-neutron, and the M-N lithology identification cross-plots which were recommended by many authors (Radwan and Nabawy 2022). Estimating the M and N values is based on the apparent matrix density (Δma) and the apparent transit time of rock matrix (Δtma) as follows.
(1) where ΔT f is the transit time of the saturating fluid, ΔT Log is the transit time log reading, ρ Log and ρ f are the density log reading and the fluid density, respectively, Ø Nf is the neutron porosity of the saturating fluid, and Ø NLog is the porosity log reading.

Shale volume determination (Vsh)
Increasing the shale volume (Vsh) causes a reduction in the reservoir quality and the formation resistivity factor (David et al. 2015). The shale volume of the Bahariya Formation has been calculated using different techniques following Asquith and Gibson (1982) and Larionov (1969), and the most actual values were recorded as follows. (

Porosity estimation
For the present study, the effective (∅e) and total porosity (∅ T ) values were calculated using the density-neutron values (Schlumberger 1995) as follows.
where Ø NLog and Ø DLog are the neutron and density porosities.

Water saturation (Sw)
Due to the complexity of the reservoir lithologies and properties, different equations have been widely applied to determine the water saturation (Sw). Among these, the Simandoux's and the Indonesian equations are the most appropriate for the argillaceous clastic sequences (Vsh > 10%) like the Bahariya Formation (Simandoux 1963;Asquith and Krygowski 2004).

Seismic data
For this study, a check-shot survey and a total of 27 seismic profiles ( Fig. 7) covering the Khalda Oil Field and its dominant structure are available for the seismic data processing interpretation. The available data have a seismic bandwidth (10-70 Hz) and a reverse polarity. Also, the borehole logs helped in the well tie process to interpret the geologic and structural elements of Khalda Field. The seismic data were processed using the Halliburton Land-mark and Petrel software to perform the synthetic generation. The seismic data were grouped into stacks of the same orientation which means that the lines of the same stack have an orientation in the same dip direction of the major geologic features. The seismic interpretation workflow starting with collecting the seismic data and checking its quality and then the well tying, the fault and horizon interpretation, and the time to depth conversion process was applied to enable mapping the structural setting of the Khalda structure (Fig. 8).
In more details, the mapped structures and horizons were chosen considering their strength and continuity. Tying has been performed between wells and horizons that intersected them in the seismic profiles. Considering the fact that the present boreholes have been drilled vertically; so it is only require to get a correct check-shot point. The two-way travel time was converted into their equivalent depth points to present it correctly on the seismic section. After interpreting all the seismic profiles, the obtained information was then presented on a base map to construct a subsurface map for the picked horizon. The most obvious type of data to get from the seismic profiles is the two-way travel time (TWT) of the picked surface that represents the mapped geologic horizon. It is the same technique for presenting the formation tops on a 2D map using the borehole log data. Additional useful information is detecting the up and down-thrown intersection points of a fault surface with the mapped horizon. Eventually, the TWT data are converted into depth values and contoured to present a structure contour map. All maps have been presented on a grid of Xinc 50 m x Yinc 50 m and contoured in a convergent interpolation.

Results and interpretation
In this section, the obtained seismic results and interpretation will be discussed. This includes creating the synthetic seismogram to convert the time to depth domain, the seismic interpretation for detecting the dominant fault blocks and mapping the different seismic horizons of the tops of different interesting zones. Also, petrophysical results and interpretation of the borehole data will be discussed using a set of X-Y and litho-saturation plots.

Synthetic seismogram
It is used to correlate between depths of the formation tops, as indicated from the well log data, and the picked horizons in the 2D seismic profiles. It is estimated from the borehole data and correlated with a particular seismic zone. The constructed seismogram is the principal link between the seismic and the borehole data, and it is the main tool that enables geological picking. This will be supported with a vertical sounding profile (VSP) as an alternative for the synthetic seismogram. In case of using combined density and acoustic logs for the studied borehole, the reflection coefficient is convolved with the wavelet to form the synthetic seismic trace. In the present study, a set of synthetic seismograms were constructed for all wells that penetrated or crossed the seismic horizons. In general, tying between the seismic data and these synthetic seismograms were performed at a satisfactory level (Fig. 9). For example, the synthetic seismogram of KHLD-56 well indicates a successful general tying with the seismic data at the studied time zone in particular.

Seismic interpretation
Interpreting the seismic data of the subsurface sequences delineates both the geological and the structural features which are presented on the seismic profiles as extended horizons. On contrary, faults are presented as discontinuities in these horizons (Abd-El Gawad et al. 2015;Abdeen et al. 2021).
In Khalda structure, the Cretaceous section was interpreted on the seismic sections covering the Khalda Oil Field. It is indicated that the Upper Cretaceous Abu Roash and the Lower Cretaceous Alam-El Bueib formations have been affected by normal faulting. Though the implication of folding is minimized, it is traceable in the field. The E-W and NW-SE gravity faults are detectable in the seismic in-lines along the E-W direction (Figs. 10 and 11). These in-line profiles show a gravity faulting with a folded structure along its downthrown side. This is due to the fact that there is an elongated folded structure affected by faulting along this level (Fig. 10). The fold extends along the NE-SW direction while faulting extends in the NW-SE direction which is in accordance with the structural regime of the north WD of Egypt during the Late Cretaceous (Moustafa 2008).
For the present study, the normal faults are primarily represented by steep fault surfaces dipping from 60° Fig. 9 A synthetic seismogram created for the Khalda-56 well 1 3 to 70°. Clustering the gravity faults in groups can create some low and high-standing fault blocks as shown in the rifted crest of the Khalda structure (Fig. 11). These faults are commonly recognized crossing the Miocene and the Pre-Miocene sequences. The dominant structural style trends to the northwest-southeast direction.

Fault blocks
For the present study, the fault blocks which are bounded from the different sides by fault surfaces which are widely presented in the central parts of Khalda structure. The triggering mechanisms of these fault blocks are referred to the Fig. 10 The interpreted TWT seismic section that E-W and pass through the Kh_58 well. It shows a normal faulting with a somewhat folded structure along its downthrown side. The fault trends NW-SE direction with its dip direction toward NE Miocene event. The majority of these blocks are tilted toward the depositional basin of the Bahariya Formation, sometimes are rotated backwards. The top of these blocks are overlain by sandstone alternated with shale sediments. In general, Khalda structure is dissected by some northwest-southeast gravity faults that bound the northward tilted blocks of some half-grabens (Figs. 10 and 11). Also, the N-S direction profiles indicate simplicity of the dominant structures represented by simple normal faults and grabens of slight throws and horizontality of the bedding planes (Fig. 12).

Fig. 11
The interpreted TWT (msec.) seismic section that extends along the E-W direction passing through the Kh_45 well. The interpreted depth structure map on top of the Abu Roash G member is shown to the right down corner of the figure

The subsurface structural maps
Four horizons were mapped for top surfaces of the Lower and Upper Bahariya members, and the Abu Roash F and G members (Figs. 13, 14, 15 and 16). The mapped structure contour surfaces indicate the presence of some large-scale NE-SW major gravity faults trending with a few NW-SE gravity faults. Additional major E-W branching gravity fault lies to S of the Khalda structure. These aforementioned step-gravity faults dissect the main hydrocarbonsbearing structure into some tilted blocks which are deeper to the north and shallower to the south. A three-way Fig. 12 The Interpreted TWT (msec.) seismic section that extends along the S-N direction, and passing through the KH_58 well. Interpreted and un-interpreted TWT (msec.) Seismic Section (XLINE 13,850), extend along the S-N direction, and passing through KH_58 well structural closure has been appeared in some parts due to the intersection of these structural blocks forming some faulted anticlines. These structural closures are caused due to the vertical basement movements.
Abu Roash F (ARF) member The depth structure contour map for the upper surface of the Abu Roash F Member, as derived from the 2D seismic data, indicates the presence of two major fault trends (NE-SW and NW-SE) with more dominance for the NW-SE trend (Fig. 13). These faults are NW-SE gravity faults with its down throw to the east forming step zone faults and a half-graben structure. The structural setting of the ARF Member is structurally high in most parts of the Khalda Field, while highest location extends from the central to the western parts and recording a depth interval 3500-5000 ft subsea. On the other hand, a low structure is presented along the NE direction and extends down to depth 5800 ft at the deepest point of the basin.
Abu Roash G (ARG) member The depth structure contour map for the upper surface of the Abu Roash G Member, as derived from 2D seismic data, shows the dominance of two main fault trends (NE-SW and NW-SE) with more dominance for the NW-SE trend (Fig. 14). The structural configuration of the upper surface of the ARG Member is structurally high area in most parts, recording the smallest depth values ranging from -5700 ft in the central parts to -4500ft to the SW. The low stand areas are situated in the eastern parts of the study area with the deepest site reaching down to -6800 ft depth. The faults are NW-SE gravity faults with the main throw located to the east forming step zone faults and a half-graben structure.
The major faults trend NW-SE directions in the study area with its down throw oriented to the east direction in some major faults and to the west direction in the other faults. These major faults were created with the seismicity initiation in the Khalda Field. The depocenter location of the depositional basin of this Cretaceous sequence is mainly Upper Bahariya member The depth structure map of the top surface of the Upper Bahariya Member, as created from the 2D seismic data, shows that it is structurally high area (Fig. 15). Its lowest depth values are recorded between depths -5500ft in the western, and -7350 ft in the central part, while the highest depth is recorded in the southern parts with depth values down to -8800ft. Most faults are trending NW-SE with some cross elements which directed along the NE-SW trend but with decrease in the faults intensity.
Lower Bahariya member The depth structure contour map on upper surface of the Lower Bahariya Member indicates that the structural high stand areas are represented by the lowest depth values between -5500 ft in the western side, and -7350 ft in the central part, while structurally low areas are located in the southern part of the area recording the highest depth values at -8800ft (Fig. 16). The predominant faulting trend is NW-SE with some cross elements which directed the fault patterns in the NE-SW but with less faulting intensity. The Lower Bahariya Formation is the main producing reservoir in this area.
Therefore, the depth intervals and their increase or decrease directions of the depth structure maps for the different horizons that were abovementioned are tabulated in Table 1. It indicates that the minimum contour values are primarily observed to the southwestern parts of the study

Results
Due to the prospectivity of the Albian-Cenomanian Bahariya Formation as indicated from the seismic study and the well log data, its petrophysical properties and reservoir characteristics were studied in detail. The total gross thickness of the Bahariya Formation including its upper and lower members is estimated as 383-415 ft in the studied wells (Table 2). Due to the structural setting and situation of the Kh_45 and Kh_S-1 × wells in as low stand in a graben between two normal faults as indicated for the seismic profiles (Fig. 12), it is fully saturated with water and indicates not oil prospectivity which primarily migrated in accumulated in the high stand sites in the center of the structure at the crest of the anticline (at kh_56 and kh_58 wells, Fig. 10). The petrophysical analysis of the estimated net-pay thickness, including the shale volume, the net to gross, the effective values, and the water saturation, indicates that the shale volume is less 32.2%, the water saturation is less than 59.0% (for Kh_56, Kh_58 and Kh_21X wells), while the effective porosity is more than 16% up to 25% (Table 2). To estimate the net-pay thickness and based on the well testing data, 10% porosity, 35% shale volume and 70% water saturation were applied as cutoff values. Based on these cutoffs, a total net-pay thickness equals is estimated as 18 to 38 feet in the center of the structure (Table 2). These prospective reservoir properties and the oil saturation resulted in a good Phi*H values (the weighted reservoir pore volume) vary between 1.79 and 3.68. The Bulk volume of water BVW (Phi*Sw) varies between 0.08 and 0.25. The obtained petrophysical and reservoir properties are tabulated in Table 2.

Interpretation
The variation of the petrophysical properties and reservoir parameters of the Bahariya Formation is primarily attributed to its structural complicated setting in the Khalda Field which is controlled by a series of normal faults forming a horst (Fig. 6), graben (Fig. 12), anticline (Fig. 10), and a step-fault zone (Fig. 11). Consequently, to  check the lithology composition and to match the promising reservoir zones a set of X-Y plots and litho-saturation plots are presented and discussed in the following section.

X-Y plots
To determine the lithology of the studied sequence, a set of X-Y plots including the density-neutron, and the M-N plots were presented in accompany with the gamma-ray (GR) values (Fig. 17a, b). The density-neutron plot for the lower and upper members of the Bahariya reservoir sequence in the kH_45 well indicates that all the measured points have porosity higher than the cutoff value (10%) up to 30% and that the readings are shifted toward higher carbonate content with GR fluctuating between 30 and 120 API with most values locates between 60 and 90 API. The API values are higher for the Upper Bahariya Member (Fig. 17a), than that for the Lower Bahariya Member (Fig. 17b). These relatively high GR readings may be attributed to the relatively high content of siltstone in the rock phase and the U content in the water content. Also, the GR values are due to the presence of many shale streaks especially in the Upper Bahariya Member which has been deposited in estuarine to shallow marine environment as referred by many authors (Shehata et al. 2018(Shehata et al. , 2020(Shehata et al. , 2021Nabawy et al. 2022). Nabawy et al. (2022), in their studies on the mineral composition of the Bahariya sequence in Shushan Basin, referred that the Bahariya sequence is relatively rich in K-feldspars (up to 19.2%), glauconite (up to 11.2%), and pyrite content (up to 8.1%) as revealed from the XRD data. This relatively high content of the K-feldspars and accessory minerals causes a decrease in the resistivity values and increase in the GR readings even in the oil-bearing zones . So, it is recommended in the future exploration plan to apply the spectral gamma logging to enable estimating the corrected GR values.
Also, plotting the M-N values on X-Y plot (Fig. 18a, b) for the Upper and Lower Bahariya Members in the Kh_56 well indicates that the measured points are characterized by relatively moderate to high GR readings (almost 30-90 API). It is also indicated that the measured points are shifted toward the higher carbonate content with relatively high Ca +2 and Mg +2 content with more shift for the upper member (Fig. 18b) than the lower member (Fig. 18a). Decreasing the M and N values may be also attributed to the relatively high acoustic and neutron porosity readings for the measured points. The silt and the carbonate content may be thought to play a main role in this behavior.

Litho-saturation plots
The litho-saturation plot of the Kh_45 and the Kh_56 wells indicates different petrophysical properties of both the upper and lower members of the Bahariya Formation with relatively high GR values up to 120 API (Fig. 19a, b). Plotting the petrophysical values of the Kh_45 well indicates relatively higher shale content in the upper member than the lower one where a few thin sand interbeds are recorded in the upper member (Fig. 19a).
In spite of increasing the porosity with increasing the sand content (tracks 6 and 10), the bulk volume of water (BVW) increases which is attributed to the full saturation of this sequence by water, i.e., it is beneath the OWC in the study area. As aforementioned, this may be attributed its low stand position, so it is recommended to avoid exploration to the NE direction where the Kh_45 is situated in a graben structure.
On moving to the SW and central parts of the study area where the anticline crest at the Kh_56 well (Fig. 19b) and the Kh_58 well, the sand content increases in both the upper and the lower members which caused increasing the potentiality of the Bahariya Formation to trap hydrocarbons. So, many movable hydrocarbons-bearing zones are assigned through the different levels of the sequence. Increasing the Table 2 The net-pay reservoir parameters derived from the borehole log data of the Bahariya Formation, Khalda Field where N/G is the net to gross ratio, ∅ eff is the effective porosity, Sw is the water saturation ratios, Vsh is the shale volume, and H is the net-pay thickness  potentiality of the Bahariya Formation in the Kh_56 well increases the net-pay thickness to 38 ft. The BVW data of this well are relatively low causing low water cut production. The relatively low deep resistivity values in both wells though being the Kh_45 well a water-producing while the Kh_56 are a hydrocarbon producing well indicates that the relatively high glauconite and pyrite content is more effective in reducing the resistivity values than that of the saturating interstitial hydrocarbons which tend to increase the resistivity values.

Regional applications
The Albian-Cenomanian deposits are a promising reservoir sequence in the Western Desert, Egypt, and in NE Africa. So, they are studied by many authors during the last decades. The present integrated petrophysical and seismic study is considered an additional study conducted on the highly prolific Albian-Cenomanian Bahariya and Abu Roash reservoirs which are widely distributed in the different promising basins of the Western Desert and produce most of the oil production of the Western Desert (Nabawy and ElHariri 2008;Shalaby et al. 2008;Bosworth et al. 2015;Teama and Nabawy 2016;Gentzis et al. 2018;Nabawy et al. 2022). As a result, the present study on this sequence is extendable to other extensions and subsurface equivalents in Egypt and NE Africa. The low resistivity values of the studied Upper Bahariya Member even in the oil-bearing zones are due of the presence of glauconite, pyrite, and conductive clay minerals which may be considered an ideal case for low resistivity reservoirs and can be applied to other similar case studies in the vicinity. As a result of this low resistivity phenomenon and the relatively high heterogeneity of the Bahariya Formation, more advanced processing techniques such Fig. 19 Litho-saturation cross-plots of the Bahariya Formation, showing the lithology, total and effective, shale volume, porosity, water saturation and net pay zones for a Kh_45 well, and Kh_56 well as the artificial neural network (ANN) should be applied in further petrophysical studies to avoid the overestimated water saturation in these low resistivity hydrocarbon-bearing sequences. This advanced technique should be run in parallel to another advanced technique to reveal the real mineral composition this type of argillaceous, glauconitic and pyritic reservoirs.

Conclusions
• Processing and interpreting the seismic data enabled figuring out the structural setting of the Cretaceous sequence in Shushan Basin based on a synthetic seismogram which indicating high reflectors strength and reflection coefficient. • Interpreting the seismic profiles delineated the implication of the normal faulting on the Albian-Cenomanian Bahariya reservoirs creating a number of horsts and grabens of slightly tilted blocks and limited throws. • Seismic achievements are supported by the petrophysical properties of the Lower and the Upper Bahariya members that revealed the presence of some prospective zones in both members with a potential net-pay thickness (18-38.0 ft). The porosity (∅e), shale volume (Vsh) and the water saturation (Sw) of the Bahariya sequence are promising (∅e ≥ 16%, Vsh ≤ 32.2%, and Sw ≤ 59.0%). • Finally, delineating the prospectivity of the Albian-Cenomanian Bahariya clastics in Shushan Basin using the proposed workflow is applicable to their subsurface extensions in the other basins in the Western Desert, NE Africa, and the NW parts of the Arabian Peninsula.