The thermal history of Kuwait was reconstructed with Shell’s “Cauldron” basin modeling software and calibrated with maturity indicators from four wells. Maturity measurements were carried out on source rock samples from wells in the Ash-Shaham, Minagish, Raudhatain and Riqua (offshore) fields. The analytical data indicate that ‘top oil window’ is located at a depth of around 2,500 meters (8,250 feet). Due to the low structural dips of the sequences in Kuwait, thermal modeling in one-dimension was found to be adequate. The Cretaceous source rocks subsided without major anomalies in burial rate or heatflow. Modeling results indicate that the Cretaceous Sulaiy, Minagish and Zubair formations entered the oil window during the Late Cretaceous and Early Tertiary, whereas oil expulsion occurred throughout Tertiary time. Quantification of expelled oil volumes suggests that the Sulaiy Formation is the most productive source rock in Kuwait, whereas the Upper Cretaceous Burgan and Maudud formations are presently at the top of the oil window.


The assessment of the maturity of organic matter is one of the main and most difficult tasks in source rock evaluation. There are many different maturity measurements, based on elemental analysis, microscopy, pyrolysis and solvent extraction (Hunt, 1995). Among the most widely used maturity indicators are vitrinite reflectance (VR), thermal alteration index (TAI), spore color index (SCI) and a whole spectrum of molecular parameters based on gas chromatography (GC) and gas chromatography-mass spectrometry (GCMS) (Peters and Moldowan, 1993). All maturity indicators, however, have shortcomings and only through cross checking of several parameters can accurate maturity levels be established.

Maturity indicators are useful in evaluating the present-day depth of the oil window. They do not, however, provide information on the timing of oil generation. Timing of oil generation, and particularly expulsion, is important in relation to the age of reservoirs and traps. Information on timing and duration of oil expulsion is obtained from charge modeling, which is based on the temperature history of a source rock and on the assumption that experimentally determined activation energies can be extrapolated to geological time scales (Connan, 1974; Hood et al., 1975; Waples, 1994; Tissot and Welte, 1984). Charge modeling software is currently undergoing a rapid development. The latest version of Shell’s “Cauldron”, for instance, allows probabilistic sampling of ill-defined input parameters (i.e. thickness and quality of undrilled source rock). In addition to timing of generation, results of charge prediction now typically include volumes and quality (API gravity) of oil expelled from the source rocks.

This study follows our earlier 1996 paper and its objective is to determine the level of thermal maturation and depth of the oil window of the Cretaceous source rocks of Kuwait. In this study, modeled maturity levels are correlated to various maturity indicators in order to reconstruct the thermal history of the area. The software was run in a Monte Carlo mode resulting in different burial histories to fit the measured maturity data.


Kuwait is located in the northwestern corner of the Arabian Gulf Basin (Figure 1). Generally, the beds in Kuwait dip gently to the northeast and the average sedimentary thickness is 6,780 meters (m) or 22,240 feet (ft). The sequence ranges from Pleistocene to Permian (Figure 2; Kuwait Report, 1989). The Mesozoic rocks in Kuwait are a thick sequence of carbonates and clastics which average 3,800 m (12,600 ft) in thickness (Figure 2, modified after Carman, 1996; Yousif and Nouman, 1997).

High tectonic activity characterized Cretaceous time and this is reflected by the three regional unconformities. These occur at the base of the Barremian Zubair, Albian Burgan and above the Turonian-Cenomanian Mishrif Formation. The regional tectonic interpretation of these unconformities is shown in Figure 2 (Carman, 1996).

The main structural features in Kuwait are a series of north-south trending anticlines associated with some minor faults (Carman, 1996). The anomalously thin formations above these structures indicates growth since Middle Cretaceous, or even earlier, as a result of epeirogenic movement. The Ahmadi anticline is the exception as it resulted from regional horizontal shortening (Fox, 1959; Carman, 1996).

Much of Kuwait’s oil is trapped in doubly plunging anticlines (Figure 1). The Kuwait Arch forms the largest structural trap in Kuwait, with culminations in Burgan, Magwa and Ahmadi fields. The Greater Burgan field covers an area of about 770 square kilometers (sq km) and has oil reserves of approximately 60 billion barrels (Adasani, 1965).

The main Cretaceous reservoirs are the Zubair, Burgan and Wara sandstones and the Minagish, Mauddud, Mishrif and Tayarat limestones (Figure 2).


Most of the Kuwait’s oil fields produce from Cretaceous reservoirs (Figure 1 and Table 1). Tertiary and Jurassic reservoirs account for some minor production. The Permian Khuff Formation, a major gas producer in most of the Arabian Gulf region, is not developed in Kuwait (Figure 2). Below the Khuff Formation the deepest well in Burgan field penetrated steeply-dipping beds of probable Precambrian age (Khan, 1989).

Murris (1980) suggested that the Cretaceous oil in Kuwait’s largest field, Burgan (Figure 1), may have migrated laterally from the Luristan Basin in Iran. The potential source rock is the Middle Cretaceous Kazhdumi Formation (Ala, 1979; Beydoun, 1993).

More recently however, Abdullah (1993) and Abdullah and Kinghorn (1996) concluded that Kuwait’s Cretaceous oil was generated from the Minagish and Sulaiy formations (Figure 2). The Sulaiy Formation is equivalent to the Makhul Formation in the Kuwait Oil Company Stratigraphic Column (Lababidi and Hamdan, 1985).

In Kuwait the deeper Jurassic petroleum system is probably separated from the Cretaceous one (Figure 2). Using wire logs, Hussain (1987) suggested that the oil accumulations in the Jurassic Najmah, Sargelu and Marrat reservoirs are generated by these same formations (Beydoun, 1993). Yousif and Nouman (1997) also highlight the good source rock signature, consisting of a high gamma ray and low density log response, in the lower section of the Najmah Formation.


Well logs from four locations, Ash-Shaham, Minagish and Raudhatain onshore, and Riqua offshore (Figure 1) were correlated and formation tops and unconformities were identified. The age of the formations and time gaps were obtained from the stratigraphic section of Kuwait (Figure 2) and the geologic time scale of Harland et al. (1989).

The four locations share a similar geologic setting and therefore only the deepest well from the Raudhatain field was used in the model. Table 2 shows the data from the Raudhatain field.

Maturity measurements were carried out on samples (cores and cuttings) taken from the same four locations. A combination of methods was used to limit the systematic errors of the individual maturity indicators and improve the accuracy of the outcome (Nuccio, 1990, 1991).

The maturity indicators used in this study include two parameters based on Rock Eval Pyrolysis: Tmax and Production Index (Peters, 1986; Espitalie and Joubert, 1987), one based on elemental analysis: (H/C vs. O/C) and the three microscopic indicators vitrinite reflectance, spore color index and thermal alteration index (Staplin, 1969; Tissot and Welte, 1984).

Most maturity indicators are dependent on the type of organic matter in the source rock. The kerogen in the Cretaceous source rocks of Kuwait were characterized as Type II and II/III (Abdullah and Kinghorn, 1996) and the values of the maturity measurements were interpreted accordingly.


Shell’s software “Cauldron” was used to reconstruct the thermal history in the area. “Cauldron” is a sophisticated basin modeling package which calculates subsurface temperatures based on the transient heatflow from the mantle and crust, the burial and uplift history, and an extensive thermal conductivity database. The hydrocarbon generation model uses the kinetics of 16 reaction species. Hydrocarbon expulsion is based on a thermally activated diffusion model (Stainforth and Reinders, 1990). The heatflow is vertical only (1D), an acceptable assumption given the predominantly shallow-dipping strata.

The following constraints and assumptions were made when applying the calculations:

  • (1) The present-day reservoir temperature at 3,000 m (92°C) is taken from a study by Watkin (1979). The geothermal gradient however, is slightly non-linear along the stratigraphic section because thermal conductivity varies due to changes in lithology. The key constraint on the model is that it is consistent with this (present-day) temperature.

  • (2) The average annual sea bottom temperature in the Arabian Gulf at present-day is 20°C (Emery, 1956; Purser and Seibold, 1973). This temperature, with variations of up to ±5°C, was taken as surface temperature from Early Cretaceous to the present. The “surface temperature” refers to a mean near-surface temperature averaged over many thousands of years. Elevation changes are negligible as the area was predominantly a stable shelf.

  • (3) The compaction and surface porosities of the various rock formations were calculated by the program, based on rock compositions given in the stratigraphic description by Lababidi and Hamdan (1985). The compaction algorithm used is based on effective stress.

  • (4) The temperatures measured in the subsurface were used to constrain the modeled present-day heatflow. This produced a present-day surface heatflow of 40 mW/m2, and mantle-derived heatflow of 25 mW/m2, both of which are consistent with the tectonic setting. The mantle-derived heatflow of the area is assumed not to have been affected by tectonic activity, and to have decayed linearly from 30 mW/m2 in the Early Jurassic to the 25 mW/m2 of present time. The high temperatures created around the Zagros collision zone (Perrodon and Masse, 1984; Fowler, 1990) do not affect the temperatures in the modeled wells, as the distance between Kuwait and the suture zone in the Zagros region (400-500 km) is too large.

  • (5) The activation energy ranges of the source rocks were estimated a-priori from the kerogen type (Type II and II/III, Abdullah and Kinghorn, 1996) and the sulfur content (moderate).

  • (6) The estimates of missing overburden are consistent with the rate of deposition before erosion. The main unconformity during the middle Tertiary was assumed to have eroded between 400 to 900 m.

The program requires geochemical parameters such as source rock TOC, hydrogen index and vitrinite reflectance. Measured TAI values were converted into an estimated vitrinite reflectance (VR/E). The best fit with the experimental data was chosen to represent the actual burial and thermal history of the area (Figure 4). The program was run in a Monte Carlo mode, so that many burial histories were created by probabilistic sampling of key parameters for the prediction of oil generation and expulsion timing. The parameters which were allowed to vary were: amount of missing overburden; heat production from shales; mantle heatflow; surface porosity of sediments; compaction coefficient; surface temperature; source rock TOC%, HI, OI, thickness and activation energies. The largest uncertainty in the timing and depth of oil generation is caused by the possible variation in the source rock kinetics.


The measured maturity data of the samples is shown in Table 3 and Figure 3. A comparison of the different parameters, taking into consideration the kerogen type, indicates that oil generation (VR 0.6%) starts at a depth around 2,500 m ± 300 m (8,250 ft ± 1,000 ft). The base of the oil window (onset of gas generation) cannot be established from this sample set.

The Sulaiy Formation has kerogen type II and although it has reached peak oil generation stage it still has an average TOC weight value of 1.5%. The Minagish Formation has an amorphous Type II and II/III kerogen. The Zubair Formation has the lowest kerogen quality, Type III, but with an average TOC weight value of 2%. Using all the input data and taking into consideration the geologic condition of the area, the model effectively concludes that some of the vitrinite reflectance data has measured maturity too low while some of the TAI data is somewhat overestimating maturity.

Applying the charge model to the corrected maturity data on the above source rocks reveals that most oil generation took place during the Late Cretaceous and Early Tertiary time, while most of the oil expulsion occurred during the entire Tertiary time (Figures 4 and 5). The Sulaiy Formation reached its peak oil generation in Late Cretaceous (75 Ma) when it was at a temperature of 100°C and a depth of 3,000 m (Figures 4 and 5).

The Sulaiy Formation reached its peak expulsion during the Paleocene (62 Ma) when the formation was at a depth of 3,750 m (Figures 4 and 5). The two peaks in Figure 5 are related to the highest temperature levels reached during the latest two subsidences which were followed by two uplifts (Figure 4). These occurred during the Paleocene and Oligocene, and are the result of the Zagros mountain building. Note that the Sulaiy has only expelled 25% of its generated oil.

Although the structural traps in the area started growing at different times (Table 4), the Zagros movements rejuvenated most, if not all, of them. This created new traps and possible migration paths from the deep Sulaiy Formation to the shallow Cretaceous reservoirs.

Calculations of the cumulative oil expelled from the three sources rocks are shown in Table 5 and Figure 6. The volumes are based on the assumption of a 150 m average thickness for the Minagish and Sulaiy formations. The 150 m average thickness for the Zubair Formation represents only the shaly facies. The average cumulative oil expelled from the Sulaiy Formation reached 780 million m3/km2, from the Minagish 156 million m3/km2 and from the Zubair 84 million m3/km2. Therefore the Sulaiy Formation is the most productive source rock in the studied geologic section in Kuwait.


The authors are grateful to the Kuwait Oil Company and Kuwait Oil Ministry for providing the samples and logs used in this study. Grateful acknowledgment is made to the University of Kuwait for supporting this research financially. Many thanks are given to Dr. J. Connan for reviewing the paper. The graphics for this paper were prepared by Gulf PetroLink.


Fowzia H.A. Abdullahis an Assistant Professor in the Geology Department of Kuwait University since 1993. She worked as a teaching assistant for three years at Kuwait University. Fowzia joined the Environmental Protection Department in Kuwait as a researcher for two years. She has a BSc in Geology and a MSc in Petroleum Geology from Kuwait University. She obtained her PhD and DIC from Imperial College, London University in 1993 in Organic Geochemistry in Petroleum Exploration.

Mark P. Ormerodis a Shell International Explorationist currently on assignment with Petroleum Development Oman, where his main areas of responsibility include prospect evaluation and basin modeling. He previously spent 5 years in the regional studies team at Nederlands Aardolie Maatschappij (NAM) in the Netherlands. Mark graduated in 1989 from Auckland University (New Zealand) with an MSc (Hons) in Geology & Geophysics. His main area of interest is the quantification of geology through the use of computer simulation and geostatistics.

Peter J.R. Nederlofis Team Leader, Regional Studies in Petroleum Development Oman. Before joining PDO in 1992, he was employed as a Petroleum Geochemist by Shell Canada in Calgary and Shell Research in The Netherlands. Peter’s main interests lay in the fields of reservoir geochemistry and hydrocarbon habitat studies. Peter obtained his MSc and PhD in Chemistry in 1973 and 1978, respectively, from the University of Amsterdam. He did a Post Doctorate at Stanford University in 1979. Peter is an active member of the AAPG and the American Chemical Society and has been a member of the board of the European Association of Organic Geochemists since 1991. Peter is a member of the Editorial Advisory Board of GeoArabia.

Robert R.F. Kinghornis currently a Lecturer in the Geology Department of Imperial College, University of London. He has over 25 years experience in teaching and research, principally in organic geochemistry, the formation of oil and gas and basin modeling. He is the author of a book and several papers on these subjects and has undertaken overseas lecture tours. He has a BSc in Chemistry from Edinburgh University, a BA in Geology from the Open University and PhD from London University.