Introduction

Southwestern Iran is among the richest hydrocarbon production areas in the world. Extensive orogenic activity during the Miocene–Pliocene had formed the Zagros Mountains with two uplifts, one in the northern Lorestan and another one in the southern Fars (James and Wynd 1965).

The study area was a structural subarea in the North Dezful just along the Zagros Fold-Thrust Belt. Extended over an area of 60,000 km2, North Dezful hosts many large oilfields (Bordenave and Hegre 2005; Jackson and Rise 1981; James and Wynd 1965). More than 90% of the crude in this area is produced from Asmari limestones (Oligo-Miocene) and Bangestan Group (Sarvak and Ilam, dated back to Cenomanian–Turonian and Santonian, respectively) (Furst 1990). The limestones of Asmari exhibit very good reservoir quality in most parts of the Dezful Embayment. This favorable reservoir quality is because of the presence of extensive systems of seams and cracks, which are even more extensive near the top of the anticlines. Asmari Formation is capped by the overlying evaporites of Gachsaran Formation (Jackson and Rise 1981). The anticlines within the Dezful Embayment are known to be developed during the Pliocene-Oligocene (Baniasad et al. 2016). Composed dominantly of limestone, Sarvak Formation is overlaid by thick layers of Gurpi and Pabdeh marls (James and Wynd 1965).

In the southwestern Iran, the Jurassic succession has driven mass production of oil and gas. It is believed that most of Iran's oil and gas is generated during the Jurassic (e.g. Sargelu Formation) and Cretaceous (Garu and Kazhdumi formations), with the generated hydrocarbons trapped in the reservoirs across the Dezful Embayment and Abadan Plain during the Cretaceous and Cenozoic (Kobraei et al. 2019). Kazhdumi Formation is believed to have served the Dezful Embayment as the main source of hydrocarbon generation. This source rock has been dated back to Albian and hosts type II kerogen (Bordenave 2002; Versfelt 2001). The Middle Jurassic sequence hosts some of the most significant source rocks in in the Zagros Basin, southwestern Iran and eastern Iraq, including the Sargelu Formation with its high total organic carbon (TOC) (Jassin and Goff 2006). Geochemical studies have revealed that the Pabdeh, Kazhdumi, Garu, and Sargelu Formations are the most important source rocks in the Zagros and Persian Gulf basins, feeding the reservoirs in different areas. For example, Pabdeh Formation has produced oil only in the north of Dezful Embayment (North Dezful) while Kazhdumi Formation serves as source in the southwestern Zagros (Dezful Embayment) (Bordenave and Hegre 2005; Karimi et al. 2016a, b). Quantitative geochemical studies and modeling have been carried out across border areas between Dezful Embayment and Abadan Plain (e.g. Ab-Teymour, Mansuri, Susangerd and Band-e-Karkheh fields). On the border between the Abadan Plain and Dezful Embayment, there is much ambiguity in the process of hydrocarbon migration, leaving it yet to be clear that whether the oils of Ab-Teymour, Band-e-Karkheh and Susangerd fields have originated from a source rock within the Dezful Embayment or came from Abadan Plain.

This study aims at investigating the characteristics of rock units from Jurassic (Sargelu Formation), Cretaceous (Garu and Kazhdumi formations) and Paleocene–Oligocene (Pabdeh Formation) in the study area. To describe the quantity, quality and maturity of the organic matter (OM) in the studied formations and estimate the maturation stage of the source rock and hydrocarbon production in the Dezful Embayment, a one-dimensional (1D) model of burial and thermal history was built. The required data was collected from four wells in the Dezful Embayment (Ab-Teymour, Band-e-Karkheh, Mushtaq-East and Ramshir). Subsequently, two-dimensional (2D) modeling was done across a selected section (northeast–southwest) to study hydrocarbon generation, migration, and storage in the study area. Here, we focus on structural features rather than considering possible lateral changes of facies and lithology, which is usually studied for identifying stratigraphic traps.

Geological setting and stratigraphy

The Zagros Fold-Thrust Belt is part of the Alpine-Himalayan Belt, which is the result of the opening and then closure of the Neo-Tethys Ocean between the Arabian Plate and Iran (Alavi 1994) in the Late Cretaceous (Berberian and King 1981). Extended along a northwest-southeast trend, this belt is about 2000 km long and starts from the Taurus Mountains in Turkey and covers all the north-northeast of Iraq and southwest of Iran before it terminates in the north of the Strait of Hormuz. In terms of geological history, the belt is known to have experienced long periods of subsidence and sedimentation accompanied by short periods of uplift and hiatus (Alavi 2004). The study area is within the Zagros Fold-Thrust Belt in the south of the so-called Dezful Embayment and North Dezful basin (Karimi et al. 2016a, b). The general structure of Dezful Embayment and the surrounding areas (e.g. upward bending, mountain front bending, Kazerun Fault Zone) as well as the north–south and east–west trends observed can be explained by the faults in the bed rock. The oldest evidence of movement in these linear structures goes back to the Late Cretaceous, although the structures remained active until the Jurassic, Triassic and perhaps even earlier in the Middle Oligocene or Miocene. Dezful Embayment is known to experience a subsidence of 300 to 600 m though it has been tectonically more stable and less exposed to folding than the neighboring areas (Motiee 1995). Figure 1 shows a location map of the study area and regional hydrocarbon fields. The action of Kazerun (right) and Balaroud faults (left) has greatly contributed to the formation of this embayment.

Fig. 1
figure 1

Fields and sections studied in the North Dezful Area. On-dimensional modeling was performed on data from four wells (Ab Teymur, Band-e- Karkheh, Mushtaq-East and Ramshir) while the data from two wells were used for 2D studies (Ahvaz and Mansuri)

In terms of geological structure, the Dezful Embayment is composed of a northern part and a southern part separated by the Izeh Fault and the highs of Hindijan. The Dezful Embayment is known to be younger than the Late-Oligocene and Early-Miocene carbonate sediments (Asmari Formation), being a subsidence basin developed since the Middle Miocene to the present (Berberian 1995). In this area, Miocene evaporites provide a good cap rock to the carbonate sediments. Sudden thickness changes in the Gachsaran Formation in this area are associated with folding, faulting, and diapiric activities occurring during and after the deposition. Erosion of the surrounding highlands, especially the Izeh Area in the northeast, has provided more sediment, developing the thickest sedimentary sequence in the Dezful Embayment. Carbonates of Asmari Formation and Bangestan Group are the most important hydrocarbon reservoirs hosted by the embayment. North Dezful is a major zone of hydrocarbon exploration and production in Iran. There are some 60 oilfields in the Dezful Embayment, some of which are so giant, such as Ahvaz, Gachsaran and Bibi Hakimeh fields. Major oilfields in the study area are Tango and Rag Sefid fields (Bordenave 2002; Sherkati and Letouzey 2004).

In the Dezful Embayment, the Kazhdumi Formation has been referred to as the source of the accumulated oils in Bangestan Group. It is the most important source rock for the reservoirs accumulated in Asmari Formation. The organic-rich calcareous shales pertaining to the Kazhdumi Formation, which were deposited in a sealed, intrashelf depression correspond almost exactly to the present-day Dezful Embayment region (Bordenave and Hegre 2005; Karimi et al. 2016a, b). Dated back to Paleocene to Oligocene, Pabdeh Formation exhibits good source rock properties in the Dezful Embayment, especially in the north. Lithologically, the Eocene–Oligocene Pabdeh Formation consists mainly of thin-bedded gray and greenish blue deep-marine pelagic calcareous shale, marl and limestone rich in pelagic micro fauna (Karimi et al. 2016a, b; Alavi 2004; Rabbani and Bagheri Tirtashi 2010). Sargelu Formation underlies Najmeh Formation (below a discontinuity surface) and overlies Alan Formation, with its age estimated as Middle Jurassic. The lithology of this formation is dominantly bituminous gray shale and clayey brown lime with occasional appearance of brown and calcareous dolomites. Dominantly reductive environmental conditions at the time of deposition of Sargelu Formation have made it a source rock of high hydrocarbon generation potential in southwestern Iran (Hosseini and Beldaji 2020; Motiei 1995). Being too deep to drill within the Abadan Plain and Dezful Area, Sargelu Formation is yet to be comprehensively studied in terms of geochemistry. As another effective source rock, Garu Formation is made up of shale and dark brown marl (Arsalan et al. 2015).

In this study, in order to evaluate the maturity of source rocks and hydrocarbon generation and outflow, Ahwaz and Mansuri oilfields were studied at a particular section. Ahwaz Oilfield is located in southwestern Iran. Being in the form of an anticline, it is 67 km long and 6 km wide and produces crude oil and gas (i.e. associated gas) at current average rates of 800,000 barrels per day and more than 13 million cubic meters per day, respectively. Mansuri Oilfield is located in the southwestern part of the Zagros Fold-Thrust Belt in the Dezful Embayment, with its crude oil production capacity estimated at 100,000 barrels per day. Figure 2 shows the generic stratigraphic column of the southwestern Iran.

Fig. 2
figure 2

source rocks

Stratigraphic column of the studied area in the southwest of Iran. The main sources are indicated by two stars while a single star is used to mark the marginal

Materials and methods

Rock–Eval pyrolysis has been widely used in exploratory activities in sedimentary basins. Using this method, it is possible to determine the quality and maturity of organic matter in the source rock and hence estimate hydrocarbon generation potential across the entire petroleum system (Taylor et al. 1998). It works by programmed heating of source rock samples in an oxygen-free neutral environment (helium or nitrogen). The parameters measured in Rock Evil pyrolysis include Tmax(°C), TOC (wt%), S1, S2 and S3 (mg HC/g rock). Using these parameters, the hydrogen index (HI = 100 × S2/TOC), oxygen index (OI = 100 × S3/TOC), production index (PI = S1/(S1 + S2)) and hydrocarbon generation potential (PP = S1 + S2) can be calculated. In this study, 27 samples of Pabdeh and Kazhdumi formations were subjected to Rock–Eval pyrolysis, with the results shown in Table1.

Table 1 Input geochemical data into the model

Basin modeling and hydrocarbon system assessment try to formulate a dynamic predictive model of the geological processes occurred during the evolution of a sedimentary basin (Baniasad et al. 2016; Hantschel and Kauerauf 2009). Modeling tools are widely used in oil exploration work and can help explorationists attain a better understanding of the petroleum systems and the interrelationships among hydrocarbon reserves (Hakimi et al. 2017). Such models have been developed to measure and simulate various processes occurred during the formation and evolution of a sedimentary basin (Baniasad et al. 2016), where the entire package of strata from the oldest to the youngest are simulated (Baur et al. 2010; Beha et al. 2008; Noeth et al. 2002).

Basin modeling can be performed one-dimensionally (1D) or two-dimensionally (2D). For this purpose, we use the geological data acquired along one or more wells, per seismic surveys, or upon structural studies along different sections to simulate the burial history and maturity of the studied section. The first step in basin modeling is to create a hypothetical model of the geological history of the well or area under study in which a sequence of geological events or layers such as sedimentation, interruption, and erosion are defined over a period of time (millions of years ago). To simulate a sedimentary basin, one should define the physical and thermal characteristics of different layers including the original and current thicknesses, lithology, current porosity, cementation, and fractures and faults; add to this the physical and thermal characteristics along the boundaries, including temperature, water–sediment boundary at the time of deposition, and the current and past heat flows. In this study, 1D modeling was performed using the data from four wells, namely Ab-Teymour, Band-e-Karkheh, Mushtaq-East and Ramshir, drilled in the Dezful Embayment. Then, 2D modeling was further performed to determine the migration paths and the location of possible hydrocarbon accumulations based on the source rocks and migration paths across a selected section (northeast–southwest, shown in Fig. 1) in the Dezful Embayment. OpenFlow™ Software was utilized to model and interpret the history of generation, migration and accumulation of hydrocarbons. At the northeast–southwest section, we considered Tertiary (with Pabdeh as source rock), Cretaceous (with Kazhdumi and Garu as source rocks) and Jurassic petroleum systems (with Sargelu as source rock). As far as the hydrocarbon sourcing is concerned, these are the most interesting formations in the Dezful Embayment (Karimi et al. 2016a, b). The well selection was done to have the maximum drilling depth and availability of hydrocarbon maturation and production data.

For 2D modeling across the northeast–southwest section, we began with entering the geological information of the section (e.g. different rock layers, age information, and fault information) into the OpenFlow™ Software. We then introduced into the model, through the Age Assignment and Facies Definition module, the stratigraphic data (e.g. layer names, top depths and ages), lithologic information, erosion and hiatus information, and the position of each layer in the petroleum system, not to mention the geochemical information related to the source formations. This was followed by entering the initial hydrogen index (HI, based on the built-in geochemistry library of the software), total organic carbon (TOC), and the kerogen type. The lithology information was determined by compiling the well data, existing lithological maps, and the stratigraphic column of the current sediments down to the bed rock upon applying the effect of the layer ages.

An important objective of thermal modeling is to measure the accuracy of the initial geological model and hence improve it by adjusting various parameters. The data required for 2D modeling can be extracted from seismic and structural sections, outcrop studies and geological data. Petrological information and chronology of the layers were retrieved from well reports (prepared by Exploration Directorate of National Iranian Oil Company, NIOC Exp.). Geochemical data on the source rock (e.g. TOC, HI, and kerogen type) were extracted from the reports published by the Geochemistry Department of the NIOC Exp. In order to achieve more accurate results, we focused on the hydrocarbon generation from Kazhdumi, Pabdeh, Garu and Sargelu formations only. For 2D modeling, a northeast–southwest cross section was designed in the OpenFlow™.

Results and discussion

In this study, 1D modeling of wells located in Ab-Teymour, Band-e-Karkheh, East-Mushtaq and Ramshir fields in the north of Dezful Embayment was performed. Calibration establishes an optimal fit between the measured parameters and the model and plays an important role in effective modeling of the studied phenomenon. The best calibration results were obtained when a constant heat flux (HF) in the range of 44 to 50 mW/m2 was considered along different wells. Measured values of vitrinite temperature and reflectance were used to check the modeling accuracy.

Rock–Eval pyrolysis

S1 (mg HC/g rock) measures the amount of free hydrocarbons that can be volatilized out of the rock without cracking the kerogen. S2 shows the hydrocarbon yield upon cracking of kerogen, representing the remaining hydrocarbon generation potential of the source rock (Magoon and Dow 1994). S1 values of the Kazhdumi Formation in the north of Dezful Embayment ranged from 0.24 to 142.5 mg HC/g rock, with an average value of 17.69 mg HC/g rock. Pabdeh Formation exhibited S1 values in the range of 0.31 to 20.43 mg HC/g rock, with an average value of 3.91 mg HC/g rock (Table 1). As an output of Rock–Eval pyrolysis, Tmax provides great insights into thermal maturity of the analyzed sample. In this work. Samples from Kazhdumi and Pabdeh Formations exhibited Tmax values in the ranges of 433–445 °C and 296–439°C, respectively (Table 1). Tmax vs. HI plot reveals that most samples of the Kazhdumi Formation were thermally mature while those of the Pabdeh Formations were either immature or in early maturation stages (Fig. 3).

Fig. 3
figure 3

Plot of Tmax vs. HI showing the types of kerogen and maturity of the studied samples in the north of Dezful Embayment

The findings show that the Kazhdumi Formation in the study area (north of Dezful Embayment) is within the oil window, making it a major source rock in the region, while Pabdeh Formation was not buried deep enough to reach the required maturity for producing hydrocarbons although it has produced very small amounts of hydrocarbon in some areas. Based on the Tmax vs. HI plot, it was determined that the studied formations were dominated by type II and a mixed species of type II and type III kerogens (Fig. 3).

Studies on Garu and Sargelu Formations in the North of Dezful Embayment area indicated that these formations serve as source rocks with moderate to high hydrocarbon generation potentials. According to these studies, the Garu Formation was dominated by type II and a mixed species of type II and type III kerogens, while Sargelu Formation contained, for most part, a mixture of type II and type III kerogens. The Garu Formation was found to possess more gas generation potential rather than oil production potential, and Sargelu Formation was found to be capable of gas generation only. According to the results of geochemical studies, Pabdeh Formation was found to be less mature than other source rocks and hence was not the main source rock in this area (Asadi et al. 2010; Hosseini and Beldaji 2020).

1D modeling

In the studied area, the Ab-Teymour well has been drilled into a maximum depth of 5040 m, reaching the Lower Cretaceous deposits of the Fahlian Formation. To achieve more accurate modeling results, the well was extrapolated by defining pseudo-layers to reach the top of the bed rock. For this purpose, information from adjacent wells was used to estimate the tops of Garu and Sargelu formations. Used to calibrate the model, the measured temperature and vitrinite reflectance data in the Ab-Teymour well are shown in Table 2a. Based on this table and upon testing different models, the corresponding HF, which is most consistent with the current bottom-hole temperature (BHT), was taken at 48 mW/m2. Figure 4a shows the maturity diagram along the Ab-Teymour well in terms of vitrinite reflectance. At this well, some 160 m of Aghajari Formation was eroded out, as per seismic data. For all fields, erosion was estimated based on seismic data. This was done by taking the maximum layer thickness at a syncline as reference and measuring the thickness at the respective anticline to identify the lost part, unless evidence of hiatus was observed on the respective anticline/syncline. According to the results of the model, Sargelu and Garu Formations exhibited vitrinite reflectance at 1.4%, making them placed in the wet gas production window. Kazhdumi Formation, with a vitrinite reflectance of 0.85%, was found to be in the oil window, and Pabdeh Formation, with a vitrinite reflectance of 0.58%, was just in the beginning of the oil window.

Table 2 Temperature and vitrinite reflectance data for model calibration
Fig. 4
figure 4

Burial history charts with thermal models (left) and temperature-depth and vitrinite reflectance-depth plots to fit the model results to the measured data (right)

The Band-e-Karkheh well was drilled into a maximum depth of 4285 m, reaching the Upper Cretaceous Sarvak Formation. With no vitrinite reflectance data from this well, we only used well temperature data to calibrate the model (Table 2b). The corresponding maturity diagram indicates that the Aghajari Formation has not been eroded at this well. At the Band-e-Karkheh well, the Sargelu and Garu formations showed vitrinite reflectance of 1.4% and 1.3%, respectively, situating them in the gas window. Finally, vitrinite reflectance of Kazhdumi and Pabdeh Formations was about 0.87% and 0.63%, respectively, placing them in the middle and early oil window, respectively (Fig. 4b).

Drilled into a maximum depth of 4994 m, Mushtaq-East well penetrated through the Upper Cretaceous Sarvak Formation. In the Mushtaq-East Field, which is relatively far from the other studied fields, Aghajari Formation was found to have some 200 m of its original thickness eroded away. The measured values of temperature and vitrinite reflectance were used to calibrate the model according to Table 2c, ending up with a very good agreement between the measured data and the modeling results. The maturity diagram of this well (Fig. 4c) shows that the Sargelu Formation was located at the end of the oil window with a vitrinite reflectance of 1.25%., while Garu, Kazhdumi and Pabdeh Formations were in the middle of oil window with vitrinite reflectance values of 1.15, 0.9 and 0.65%, respectively.

In the Ramshir Field, the Aghajari Formation was found to have some 1100 m of its original thickness eroded away. With a maximum drilling depth of 4660 m, the field was studied down to Sarvak Formation. Given the unavailability of vitrinite reflectance data, only temperature data was used for calibrating the model in this field, leading to a good agreement (Table 2d). Figure 4d shows the maturity of the studied formations at the studied well in Ramshir Field. As shown in the figure, Sargelu and Garu formations exhibited vitrinite reflectance values of 1.8% and 1.7%, respectively, categorizing them in the wet gas production window. Moreover, Kazhdumi Formation, with a vitrinite reflectance of 1.17%, was in the beginning of the gas window. Finally, Pabdeh Formation showed a vitrinite reflectance of 0.74%, which showed that the formation was in the oil window. According to the results of 1D modeling, the studied source rocks were most matured in Ramshir Field, possibly because of their deeper burial depth in this field.

2D modeling

As mentioned earlier, 2D modeling was performed on a selected section (with a northeast–southwest trend) in the study area. The software compiled the seismic, stratigraphic and geological interpretation data to carry out multidimensional simulations of HF, fluid flow and oil migration in the sedimentary basin. As outputs, the 2D modeling gives estimations of the oil/gas generation from the source rock, migration through the carrier rock, and accumulation in the reservoir rock. It further returns the formation pressure history. As shown in Fig. 1, the selected section cut through Ahwaz and Mansuri fields, where the well data were available to evaluate the modeling accuracy. Validation of the modeling results at well location against the well data indicated a good agreement (Table 3). Simultaneous calibration of temperature and vitrinite reflectance and the uniformly increasing trend of vitrinite reflectance with depth highlighted the mild tectonic history of the region (Fig. 5).

Table 3 Temperature and vitrinite reflectance data for calibrating the 2D models of Ahwaz and Mansuri fields
Fig. 5
figure 5

Vitrinite reflectance-depth and temperature-depth diagrams for fitting the models to the measured data

The drilling of Ahwaz well has been completed after reaching a depth of 5000 m in Fahlian Formation. To achieve more accurate results in modeling, the layers were extrapolated to the basement. According to Fig. 6a, the deepest sediments are related to the Lower Fahlian Formation, i.e. Lower Cretaceous, some 145 million years ago. According to Fig. 6a, the Middle-Jurassic Sargelu Formation exhibited a vitrinite reflectance of 1.50% and was found to be within gas production window. In the Middle Cretaceous, the same formation (Sargelu) was found to enter the oil window, with the largest amounts of oil produced from this formation during the Upper Cretaceous. Further at this well, the Garu Formation showed a vitrinite reflectance of 1.40%, placing it is gas production window. But previously, this formation was found to have entered the oil window in the Upper Cretaceous. With a vitrinite reflectance of 0.98%, the Kazhdumi Formation was found to be in the oil window since after the Eocene. The Pabdeh Formation showed a vitrinite reflection of as low as 0.63%, situating it in early stages of oil production where the generation rate is very low.

Fig. 6
figure 6

Burial history charts of Ahwaz a and Mansuri b wells

Completed at a depth of 5000 m, Mansuri well penetrated down to the Lower Fahlian Formation dated back to some 145 million years ago (Fig. 6b). According to the modeling results, Middle-Jurassic Sargelu Formation at Mansuri Well is currently located at a depth of 5773 to 5743 m. Based on Fig. 6b, the Sargelu Formation had reached a vitrinite reflectance of 0.55% some 100 million years ago from when it began to produce hydrocarbons. Presently, this formation has a vitrinite reflectance of 1.25%, implying its position at the end of the oil production window. At this well, the Garu Formation occurred at a depth of 5523 to 5563 m, entering the oil production window in the Upper Cretaceous. At present, this formation exhibits a vitrinite reflectance of 1.15%, keeping it in the oil window. Presently occurred at a depth of 4146–4352 m at Mansuri well, Kazhdumi Formation was dated back to the Middle Cretaceous. It has been found to enter the oil window in the Middle Miocene with a vitrinite reflectance of 0.55%, and is currently at the peak of oil production window with a vitrinite reflectance of 0.37%. Focusing on Pabdeh Formation, Fig. 6b indicates that the vitrinite reflectance had reached 0.55% some 4 million years ago and is now at 0.58%, locating the formation in the beginning of the oil production window. Like the Pabdeh Formation at Ahwaz well, hydrocarbon production It has been insignificant.

Regarding its significant role in the maturity of the source rock and hydrocarbon production and migration, temperature is an important parameter in any geochemical assessment. Basin modeling provides a unique opportunity to study temperature regimes in space and time. Accurate determination of the temperature history is possible only if the heat transfer process is examined in terms of real physical dynamics under the given boundary conditions in the sedimentary basin and its evolution over time. The temperature is also important as it affects the physical properties of sediments and hydrocarbon fluids. As a result, many processes occurred during the evolution of a basin, such as compaction and fluid flow, are directly controlled by temperature. Therefore, a realistic reconstruction of the temperature history is critical for modeling basin evolution and understanding the corresponding highly complex processes and their interactions during basin development (Welte et al. 1997).

In the basin modeling, spatial modeling of temperature variations is of paramount importance as it reveals thermal properties of different rock units and their changes across the underground space. One of the most difficult and sensitive tasks in the basin modeling is to assign reliable values to various geological units comprising the modeling domain.

The modeling results showed that the maximum temperature recorded in the Pabdeh Formation was about 90 °C, which was less than enough, given the shallow burial depth of the formation, to make the rock unit so mature that hydrocarbon can be produced. Kazhdumi and Garu Formations were found to have experienced temperatures in the ranges of 90–105 °C and 160–180 °C, respectively, with Sargelu Formation underwent a maximum temperature of 173 °C (Fig. 7).

Fig. 7
figure 7

Modeling results in terms of temperature profile along the northeast–southwest section

Figure 8 shows the maturity of different formations in the study area in terms of vitrinite reflectance. The figure suggests that the formations had their maturity increased with depth in the synclines. According to the results, the eastern part of the study area had just stepped in the oil window, where the Kazhdumi Formation showed a vitrinite reflectance of about 0.8% on the ridge of the anticline. However, with increasing the depth (synchronization), this value reached as high as 1.3%, making the respective rock units beyond the oil window. The Garu and Sargelu Formations were found to be more mature in synclines rather than anticlines, although the TOC of the Pabdeh Formation was not mature enough to produce hydrocarbons, possibly due to the shallow burial depth. Results of the present maturity model showed a generally decreasing trend of maturity in the source rocks as one moved from the Dezful Embayment to the Abadan Plain.

Fig. 8
figure 8

Plot of the vitrinite reflectance versus time along the studied section

More intensive heat fluxes were observed wherever the bed rock occurred at a shallower depth. The increase in the thickness of sediments in the east of the model increased the depth of burial compared to the east, further contributing to the maturity of the source rocks on this side of the model. Also, based on the results of the 1D modeling of the four wells located in the Ab-Teymour, Band-e-Karkheh, Mushtaq-East and Ramshir fields (in the vicinity of the studied fields, namely Ahwaz and Mansuri) and also those of the 2D modeling, it could be stipulated that Ramshir Field in the east of the section was the most mature field in the region, with the vitrinite reflectance values of the studied source rocks decreased as one moved to the west toward the Ab-Teymour, Band-e-Karkheh and Mushtaq-East fields (Fig. 1). The measured values of the vitrinite reflectance at wells indicated a decrease in the maturity of the source formations as one went from east to the west. In fact, the increased heat flux and the depth of burial at the present time contribute to higher maturity of the source rocks in the east.

An important outcome of basin modeling is the identification (or prediction) of hydrocarbon migration paths and accumulations within the petroleum system. Figure 9 presents the identified migration paths through which the hydrocarbons migrated from source rock to reservoir rock where they were accumulated. The hydrocarbon generation and migration processes can be explained by Darcy's law and Arrhenius kinetic law, respectively (Schneider 2003; Ungerer et al. 1990). Our migration model showed a possibly large accumulation of hydrocarbon beneath the Kazhdumi Formation. At depths below the Kazhdumi Formation, such hydrocarbon accumulation was predicted to be located in between the studied fields, while shallower reservoirs (i.e. above Kazhdumi Formation) contained such accumulations within them. The results of 2D modeling indicated that the petroleum systems developed during the Early and Middle Cretaceous were occasionally connected via faults and special geometry features, with the two separated by the Kazhdumi Formation, which not only served as the source rock for the overlying layers but also sealed the underlying strata as cap rock. This formation was found to have reached the peak of oil production window in the depression area on the east of Ahwaz Field and also in the inter-field areas (Karimi et al. 2016a, b). According to Fig. 9, hydrocarbon migration from the Dezful Embayment to the Abadan Plain occurred at a relatively high intensity. The same was seen in Sarvak Formation in a similar direction but at a lower intensity, while the Ilam Formation had nothing to do with the hydrocarbon migration. Due to the intensity of lateral migration in sub-Kazhdumi reservoirs as well as some of the Sarvak reservoirs, the formations can have their porosity and permeability altered in the inter-field areas, increasing the chances of developing stratigraphic traps, calling for further studies. As shown in Fig. 9, the lighter the weight of the hydrocarbons produced from Kazhdumi Formation, as compared to other fluids (i.e. water), reservoir pressure, presence of vertical migration paths, and relatively high porosity of the reservoir horizons have made the produced hydrocarbon migrate vertically up into Sarvak Formation or horizontally/obliquely into Ilam and Asmari formations on top of the anticline (Zeinalzadeh and Rezaei 2010). The sealing effect of Kazhdumi Formation, however, has prevented the hydrocarbons from leaving the anticline. Ilam and Sarvak formations are dominated by highly fractured carbonates that provide for higher permeability than Pabdeh and Gurpi formations. This situation prevents the evasion of hydrocarbons across the ridges toward Pabdeh and/or Gurpi Formation. Indeed, post-migration internal movement of hydrocarbons within Ilam and Sarvak formations increased the hydrocarbon saturation along the ridge of the anticline. The increase in the saturation creates more thrusts, which together with the cracks and fissures caused on top of the anticline upon folding provide a suitable path for the vertical migration of hydrocarbons into the Asmari Formation. The absence of cracks and high saturation of the areas along the ridges limited the migration of hydrocarbons into Pabdeh and Gurpi formations (adapted from in-house reports of NIOC Exp. and the results of its modeling project).

Fig. 9
figure 9

An output of the hydrocarbon migration and accumulation model along the studied section

The pressure data at the wells indicate a high-pressure zone in the reservoirs underlying the Kazhdumi Formation, especially Fahlian. In the shallower reservoirs (e.g. Asmari, Ilam and Sarvak), no sudden change was observed in the increasing trend of pressure. The Kazhdumi Formation plays two roles in the studied petroleum system, namely as source rock and cap rock, separating the former Cretaceous petroleum system from the Middle Cretaceous-Miocene one. The strength of this cap rock is so high that it can tolerate sudden changes in the pressure of the underlying reservoirs of the Cretaceous petroleum system. Figure 10 indicates the effective roles of Kazhdumi and Gutnia Formations as cap rock. Here, the Gachsaran cap rock and Pabdeh and Gurpi Formations were absent. The most important point regarding the northwest-southeast section is the distribution of the HF across the North Dezful Embayment. As revealed by the results of 1D modeling and measured data at wells (e.g. temperature and vitrinite reflectance), optimal HF rate was between 44 and 50 mW/m2, where the best fit between the measured data and modeling results was obtained. The HF was maximal on the anticlines and generally increased as one moved toward the Abadan Plain (Fig. 11).

Fig. 10
figure 10

High-pressure zones across the studied cross section

Fig. 11
figure 11

Variations of HF across the cross section

Conclusions

According to the geochemical studies, it was found that, in the North Dezful Area, in contrast to the Kazhdumi, Garu and Sargelu Formations, the Pabdeh Formation is yet to be mature enough to produce hydrocarbon. In this study, 1D basin modeling was performed based on the data from four exploratory wells drilled in different fields across North Dezful, including Ab-Teymour, Band-e-Karkheh, Mushtaq-East and Ramshir fields. For this purpose, the maturation and hydrocarbon production processes were investigated in the North Dezful. Well selection was done to ensure maximum drilling depth and hence enable calculations for deeper source rocks at an acceptable level of accuracy. 2D modeling was performed by defining a cross section (northeast—southwest) and then evaluating the maturity levels of kerogen (s) of different sedimentary strata.

The best fit between the measured data and modeling results was obtained with HF values in the range of 44 to 50 mW/m2. Simultaneous calibration of temperature and vitrinite reflectance and the uniform increase in the vitrinite reflectance with depth indicated mild tectonics of the region. According to the results obtained from the 1D modeling, the Sargelu and Garu Formations were found to be in the gas window, and the Kazhdumi Formation was in the oil window, making it possible to effectively charge the regional fields, considering the timing of the Zagros Orogeny. Based on the modeling results, the OM of the Pabdeh Formation was immature, producing no hydrocarbons. The Pabdeh Formation showed an average temperature of 90 °C and vitrinite reflectance of 0.6%. Dated back to Albian, the Kazhdumi and Garu Formations were found to have experienced temperatures in the ranges of 90–105 °C and 160–180 °C, respectively. Sargelu Formation was found to underwent a maximum temperature of 173 °C.

The results of the migration modeling across the studied section showed that the Early Cretaceous and Middle Cretaceous hydrocarbon systems were separated by the Kazhdumi Formation, inhibiting the hydrocarbon migration into deeper layers of Kazhdumi Formation. This often leads to lateral migration of the hydrocarbons into the Abadan Plain. In the Asmari Formation, the charging was found to have occurred within the field due to the presence of faults in the fields. Our findings highlighted the changes in the temperature and maturity of the OM across different parts of the anticline. The Early and Middle Cretaceous petroleum systems were found to be connected through particular geometry features or possibly faults, with the two systems separated by Kazhdumi Formation.