Basin and petroleum system analysis in the southeastern Persian Gulf basin: a 2D basin modeling approach

A regional 2D conceptual model was constructed to study the hydrocarbon systems associated with the Middle Sarvak source rock (Cenomanian) in the southeastern Persian Gulf basin. The 2D cross section used for this purpose was 135 km long and encompassed a total of 17 depositional units with the Late Jurassic Surmeh Formation at the base and the Aghajari/Mishan sequences at the top. Compositional type II-S kinetics was considered for the Middle Sarvak source rock, and the model predictions were calibrated against observe data (corrected bottom hole temperatures and vitrinite reflectance data) coming from 4 wells located on the studied cross section. Our results indicate that hydrocarbons predominantly migrate laterally from east to the west of the study area, thereby sequentially charging potential traps en-route. The interaction between facies patterns and geodynamic evolution controls both the distribution and quality of hydrocarbon accumulations in the study area. The calibrated 2D model reliably predicts most of the present-day hydrocarbon occurrences in the study area and explains the present-day variations in their bulk properties. Our findings provide new insights about the unconventional plays associate with the Middle Sarvak source rock especially in the central parts of the southeastern Persian Gulf basin.


Introduction
The Persian Gulf basin is renowned for huge volumes of gas accumulations in the Paleozoic series and considerable amounts of oil reserves within the Jurassic through Cretaceous sediments (Ghazban 2007). Geochemical studies on oils and source rock samples from this basin have testified the presence of several important petroleum systems in the area Hosseiny et al. 2017;Rabbani et al. 2014).
Recent seismic studies in the Iranian sectors of the prolific Persian Gulf region (OEOC, 2014) has resulted in numerous oil and gas discoveries within the Cretaceous Series. The Middle Sarvak Formation (Cenomanian) is believed to be an important source rock within the Cretaceous petroleum system of the eastern Persian Gulf basin (Alipour, 2017;Hosseiny et al. 2016). This source rock is suggested to charge the overlying Mishrif and Ilam carbonates in some oilfields of the study area Alizadeh et al. 2017;Hosseiny et al. 2017). Stratigraphic equivalents of the Middle Sarvak source rock have been geologically and geochemically assessed in the United Arab Emirates (the Shilaif Formation) (Al-Zaabi et al. 2010;Azzam and Taher, 1993;Loutfi and El Bishlawy, 1986) and in the Iraq (the upper parts of the Qmachuqa Formation) (Ameen and Gharib, 2014). The existing organic geochemical data and basin modeling results from the United Arab Emirates collectively suggest that hydrocarbons accumulated within the Mishrif reservoir are sourced from the Cenomanian Shilaif source rock.
Notwithstanding the large volume of geochemical research, the dynamic evolution of the Cretaceous petroleum system has rarely been documented in the eastern Persian Gulf basin and neighboring countries. Previously, modeling techniques were applied to highlight the geohistory evolution of local kitchens in the study area (Alipour, 2017;Alipour et al. 2017). These studies suggested that the timing and extent of hydrocarbon generation from the Middle Sarvak source rock were predominantly controlled by variations in the depth of burial. Previous modeling studies have also suggested that local stratigraphy can critically control the distribution of accumulated hydrocarbons (Alipour et al. 2016a). These studies suggested that biodegradation processes can influence generated hydrocarbons at structurally shallower regions.
Multi-dimensional modeling can be particularly useful to study the generation of hydrocarbons and their migration/ accumulation histories. In addition, these techniques provide new insights about remaining exploration potential of an area by identifying new stratigraphic or unconventional plays (Alipour et al. 2021;Baur et al. 2017;Hantschel and Kauerauf, 2009;İnan et al. 2017). The aim of this study was to gain a deeper insight about the evolution of the Cretaceous petroleum system in the eastern Persian Gulf basin. For this end, a regional 2D model was constructed across the depositional basin of this source rock (Fig. 1). Such cross sections are believed to provide critical information about the elements and processes of the petroleum system (Magoon and Dow, 1994;Peters et al. 2017). Our main emphasize was on analyzing the generation and migration of hydrocarbons from the Middle Sarvak source rock. In addition, we use modeling techniques to study the present-day variations in the bulk properties of accumulations and explain the heavy bitumen occurrences in some fields. Finally, the potential for unconventional oil resources within the Middle Sarvak source rock is assessed based on our modeling results.

Background geology
The Mesozoic stratigraphy of the Persian Gulf basin is characterized by episodic development of intra-shelf basins over the northeastern margin of the Arabian plate (Sharland et al. 2001;Van Buchem et al. 2011). During Aptian and Cenomanian, these basins resulted in deposition of organicrich sediments which adjoined porous rim carbonates on  (Vahrenkamp et al. 2015). The Cenomanian intra-shelf basin was developed following the deposition of the lower parts of the Sarvak Formation, which are regionally continuous and dominated by carbonates (Razin et al. 2010). These carbonates are overlain in the west and east by shallow-marine carbonates of the Mishrif platform rimming an NE-SW trending intra-shelf basin in the central parts of the study area (Vahrenkamp et al. 2015) (Fig. 1a,b). The organic-rich facies of the Middle Sarvak Formation was deposited within this intra-shelf basin, while the riming carbonate platforms could be exposed at times of relative sealevel fall (Al-Zaabi et al. 2010;Van Buchem et al. 2002).
The infill of this basin was provided by deposition of the upper parts of the Sarvak Formation during the subsequent transgression and increase in carbonate production (Vahrenkamp et al. 2015). A major sea level fall during the Turonian resulted in subaerial exposure of the uppermost parts of the Sarvak Formation in most areas across the Persian Gulf basin (Hajikazemi et al. 2010;Motiei, 1993) (Fig. 1c). Following this regional Turonian unconformity, exceptional subsidence (due to ophiolite obduction) led to the deposition of the Laffan Formation (Lawa et al. 2013;Omar et al. 2015;Sharland et al. 2001). This formation is composed of black shales in the study area and acts as a potential seal rock for the underlying Mishrif reservoirs (Alipour, 2017;Ghajar et al. 2011). However, these shales exhibit a general trend of decreasing thickness eastwards in the Persian Gulf basin, which can cause occasional escape of hydrocarbons into the Campanian Ilam carbonates above (Al-Zaabi et al. 2010;Alipour et al. 2016a). In such cases, the thick shales of the Gurpi and Pabdeh formations can act as the ultimate seal rock for these migrating hydrocarbons (Fig. 1c).

Input data
For the purpose of the present study, a regional 2D cross section was extracted from a 3D geological model of the studied area. The section was oriented to encompass the Cenomanian intra-shelf basin (Figs. 1a, 2). We followed a two-step workflow for defining the thickness and geographic extension of the source rock interval in our model. In the first step, well logs (including gamma ray, sonic and resistivity), lithology data (carbonate/shale ratio), and organic geochemistry data (TOC and HI) from a total of 18 offshore wells were used to delineate the extent of the Cenomanian intra-shelf basin in the study area (Figs. 1, 2). The organic-rich intervals are often characterized with high gamma ray readings (e.g., 80-110°API), low sonic transit time, and high resistivity. Similarly, these zones exhibited lower carbonate/shale ratios with higher TOC and HI readings (Fig. 3a). During the second step, the same information were used to define the thickness of the organic-rich interval in our model. The organic-rich zone is generally thicker in the center of the intra-shelf basin, exceeding 40 m in the vicinity of Well-II (Fig. 4). However, its thickness decreases to less than 30 m on the eastern and western margins of the intra-shelf basin until it is eventually replaced by shallow-water, organic-lean carbonates of the Mishrif Formation (Figs. 2, 4).
The main input parameters for constructing the 2D conceptual model, including the geological (thickness, ages, paleo-bathymetries and lithofacies information) and geochemical (Rock-Eval pyrolysis and organic petrography) information, were also obtained from the same 18 wells (Table 1 and Fig. 2).
Geochemical characteristics of the Cenomanian source rock were obtained from Rock-Eval pyrolysis data of a total of 35 samples from different oilfields (see also Hosseiny et al. 2016) (Figs. 2, 3a). The type of organic matter was defined based on previous studies as reactive type II with relatively high sulfur content (Alipour, 2017;Alipour et al. 2016bAlipour et al. , 2019Hosseiny et al. 2016). Due to the lack of published information regarding hydrocarbon generation kinetics from this source rock, a compositional kinetic scheme corresponding to type II-S (Behar et al. 1997) was selected from the default library of PetroMod® 2012 software (Fig. 3b).

Model building
The boundary conditions of the model were set based on our experience from previous modeling studies in the southeastern Persian Gulf basin (Alipour et al. , 2016a(Alipour et al. , 2014. In brief, considering the tectonic quiescence of the study area during the Cretaceous, a constant heat flow of 50 mW/ m 2 was considered through time using McKenzie Model of the PetroMod® 2012 software (Table 1). In addition, a constant thickness of eroded sediments was considered for the upper parts of the Sarvak Formation following the Turonian unconformity (Table 1).
The obduction of ophiolites during the Late Cretaceous is suggested to have caused reactivation of old basement features and initiation of halokinetic movements across the Arabian Plate (Sharland et al. 2001). However, these processes are not well constrained in the study area. Therefore, we did not consider them in the present study because their insertion into our models would incur large uncertainties upon the predictions. The temperature at the sediment-water interface (SWIT) was automatically calculated using the PetroMod® 2012 functionality based on the pre-assigned paleo-water depth (PWD) information (Wygrala, 1989).
Modeling predictions were calibrated against measured values of present-day temperature (i.e., corrected bottomhole temperature) and vitrinite reflectance data (following 1 3 the method of Sweeney and Burnham, 1990) from four wells located on the modeled cross section ( Fig. 2 and Table 2).

Model calibration
After running a simulation on the constructed 2D model, the modeling results should be checked against observation data. A satisfactory calibration means that a reasonable agreement exists between modeling predictions, on the one hand, and the natural conditions, on the other (e.g., measured temperature and vitrinite reflectance data from existing wells). Calibration profiles (Fig. 5) indicate a satisfactory match between model predictions and the observe data, meaning that thermal evolution of the studied source rock can be reliably evaluated using our modeling approach.

Timing and extent of hydrocarbon generation
Thermal maturity profiles for the studied 2D section indicate relatively higher maturity and elevated transformation ratio in the eastern parts (Fig. 6). Generally, the onset of hydrocarbon generation from the Middle Sarvak source rock corresponds to the deposition of the Gurpi and Pabdeh formations (i.e., from c.a. 60 Ma). These formations have a relatively higher thickness in the eastern parts of the studied section (Figs. 4, 7), and therefore, the earliest phases of hydrocarbon generation would correspond to these deeply buried areas (Fig. 7a).  Azzam and Taher (1993). The distribution of offshore wells used for constructing and calibrating the 2D conceptual model along the studied transect is also shown When such early hydrocarbon generation phases occur within a basin, two critical factors should be carefully examined: (a) timing of trap formation, and (b) seal efficiency.
If these requirements are not satisfied in the studied basin, small chances of preservation will exist for the generated hydrocarbons. In our study area, structural deformation  . 4 Snapshot of the present-day stratigraphy of the 2D conceptual model constructed in the eastern Persian Gulf basin (see Fig. 2 for section location). Data from four wells located on the cross section were used for calibration  associated with salt tectonics was active since Late Paleozoic (Motamedi and Gharabeigli, 2019) or Jurassic times (Edgell, 1996). These episodic salt movements resulted in the erosion and karstification of the upper parts of the Sarvak Formation in the crest of many structures in the Persian Gulf basin (Ghazban, 2007;Hajikazemi et al. 2010;Videtich et al. 1988). This observation suggests that traps in our study area were available prior to the beginning of hydrocarbon generation from the source rock. Similarly, interbedded shales of the Laffan Formation or the thick shales of the Gurpi and Pabdeh formations could act as efficient seal rocks for the Upper Sarvak and Ilam reservoirs in the study area. Another phase of pronounced burial and enhanced hydrocarbon generation is evidenced during the deposition of the Middle Miocene Aghajari/Mishan sequences (i.e., from c.a. 13 Ma). Again, the relatively higher thickness of these sediments on the eastern parts could result in relatively higher thermal stress for the source rock compared to the western parts (Fig. 7). These observations suggest that the geodynamic evolution of the study area is largely responsible for the observed variations in the hydrocarbon generation behavior of the Middle Sarvak source rock.
Westward, the level of thermal maturity decreases for the Middle Sarvak source rock due to shallower burial until the source rock lies in the early stages of oil window near Well-II (Figs. 6, 7). This regional trend of maturity evolution has direct implications for defining the geographic extension of the pod of active source rock in the Cretaceous petroleum system. In addition, it suggests that lateral migration would be the prevalent form of fluid migration from the east toward the west of the study area (i.e., from kitchen areas in the east toward accumulations sites on the west) (Fig. 8).

Migration and accumulation
Fluid flow migration was performed on the 2D model using the hybrid method that considered Darcy-flow and flowpath analyses to predict both the magnitude and direction of migration vectors (Hantschel and Kauerauf, 2009). Based on the migration modeling results (Fig. 8), generated hydrocarbons experience upward vertical migration into the porous carbonates of the Upper Sarvak Formation. These carbonates are composed of highly porous facies, especially in the eastern parts of the studied section. Therefore, they can act as potential avenues for upward lateral migration of hydrocarbons (Fig. 8).
The impermeable Laffan Formation provides a relatively good top-seal for these carrier rocks. Therefore, any geometric closure (structural and/or stratigraphic) within the Upper Sarvak Formation would serve as a potential accumulation site for migrating hydrocarbons. Accordingly, three prospects were predicted by our model, which correlated well with the present-day distribution of producing fields in the   (Fig. 8). The measured °API gravity data reported from these oilfields further support the westward migration of hydrocarbons in the study area. The oilfield located in the vicinity of Prospect-I contains an oil with 30°API gravity in the Upper Sarvak reservoir, while other oilfields located near Prospect-II and Prospect-III produce oils with 27°API and 26°API gravity, respectively. This observation suggests that Prospect-I is closest to the kitchen area and that hydrocarbons generally migrate westward in the study area.
Our modeling results also provide evidence to support that the existing °API gravities are not a result of secondary alteration processes such as biodegradation. For this, we analyzed the trap charge dynamics in the three identified prospects. A synthesis of the Upper Sarvak cell histories ( Fig. 9) indicates that the major phase of hydrocarbon emplacement has occurred at c.a. 13 Ma in all three prospects. Interestingly, the temperature of the reservoir cells had already reached pasteurization temperatures around this time (Fig. 9). Head et al. (2003) reported that biological activities tend to be considerably reduced at temperatures above 80 °C. Therefore, the Upper Sarvak reservoir conditions would be favorable for preservation of early charges from degradation by micro-organisms.
However, consistent with the general geology of the study area, modeling results suggest that the risk of finding biodegraded reserves will increase moving upward in the studied section (e.g., within the low-temperature Ilam reservoir). This is consistent with previous modeling studies suggesting moderate levels of biodegradation in the Ilam reservoirs of the eastern Persian Gulf basin (Alipour, 2017;Alipour et al. 2017Alipour et al. , 2016a. In addition, the chance for biological alteration of hydrocarbons progressively increases westward within the studied section (i.e., shallower, cooler regions) (Head et al. 2010(Head et al. , 2003. In other words, migrating hydrocarbons within the Upper Sarvak carbonates would progressively rise up to shallower depths  Table 2) Fig. 6 Evolution of the thermal maturity along the studied section highlighting the transformation ratio (%) of the Middle Sarvak source rock in the vicinity of Well-I and Well-II Fig. 7 Comparison of burial history plots from Well-I (a) and Well-IV (b), illustrating the difference in burial depth and maturation level from east to west of the studied section during updip migration. This is consistent with geological observations regarding the occurrence of heavy bitumen in the Upper Sarvak carbonates at several fields lying to the west of the study area. Therefore, due to the influence of extensive biodegradation, the majority of these hydrocarbons will appear in the form of heavy oil (bitumen) impregnations.
Another important implication that can be drawn from migration modeling results is that a new unconventional play can be associated with the Middle Sarvak source rock, especially in the central parts of the basin (Fig. 10). Several factors contribute to the association of promising shale oil resource potential with this source rock in mentioned areas: (a) the lower thermal maturity of the source rock equivalent to the early oil window, (b) the relatively higher thickness of the source rock interval (i.e., lower expulsion efficiency), (c) the lower permeability of the Upper Sarvak carbonates immediately lying above the source rock which prohibits effective drainage of generated products, and (d) the presence of a thick layer of the impermeable Laffan shales which acts as an ultimate sealing rock (Fig. 10). In short, the Middle Sarvak source rock can be a self-contained unconventional resource system in the central parts of the study area. Similar findings are reported from the United Arab Emirates where the organic-rich Shilaif source rock is in the early stages of oil window and overlain by the tight facies of the Tuwayil Formation (Taher, 2010). Nevertheless, we believe that detailed 3D modeling studies are required to fully capture existing facies variations and to better assess the hydrocarbon migration processes in this system.

Conclusions
A calibrated 2D conceptual model was constructed to investigate thermal evolution of the Middle Sarvak source rock (Cenomanian) and analyze the associated petroleum system. The results indicate a strong lateral migration occurring from east to the west of the study area. In addition, the modeling results suggest that the present-day variations in bulk properties of the Upper Sarvak oils are largely controlled by the initial charging mechanisms rather than biodegradation. However, occurrence of biodegradation cannot be neglected in two cases: (a) in the eastern parts where the Laffan shales are thin and hydrocarbons escape vertically into the shallower Ilam reservoir, (b) in the western areas where some hydrocarbons laterally migrate updip inside the carrier rocks. Eventually, our modeling results provide insights about the new unconventional oil paly associated with the Middle Sarvak source rock in the central parts of the eastern Persian Gulf basin.  Acknowledgements The authors would like to thank department of "petroleum geology and sedimentary basins" of Shahid Chamran University of Ahvaz for supporting this research. Iranian Offshore Oil Company (IOOC) is warmly acknowledged for providing the data and permission to publish. The corresponding author is grateful to the Drilling Research Center (DRC) of Shahid Chamran University of Ahvaz for their support. Comments from anonymous reviewers on earlier versions of this manuscript are acknowledged with thanks.
Funding No funding was provided for doing this research.

Conflict of interest The authors declare no conflict of interest.
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:// creat iveco mmons. org/ licen ses/ by/4. 0/.