Next Article in Journal
The Surface-to-Atmosphere GHG Fluxes in Rewetted and Permanently Flooded Former Peat Extraction Areas Compared to Pristine Peatland in Hemiboreal Latvia
Next Article in Special Issue
Geogenic Sources of Arsenic and Fluoride in Groundwater: Examples from the Zagros Basin, the Kurdistan Region of Iraq
Previous Article in Journal
Changes in Climate Factors and Catches of Fisheries in the Republic of Korea over the Three Decades
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Meteoric Water Incursion, Crude Oil Degradation and Calcite Cementation of an Upper Cretaceous Reservoir in the Zagros Foreland Basin (Kurdistan Region of Iraq)

1
General Directorate of Scientific Research Center, Salahaddin University-Erbil, Erbil 44002, Kurdistan Region, Iraq
2
Department of Geology, Faculty of Science, Palacký University Olomouc, 17. listopadu 1192/12, 771 46 Olomouc, Czech Republic
3
Department of Geology, College of Science, University of Sulaimani, Sulaymaniyah 46001, Kurdistan Region, Iraq
4
Petroleum Engineering Department, College of Engineering, Al-Kitab University, Kirkuk 36001, Iraq
5
Department of Petroleum Engineering, Soran University, Soran 44008, Kurdistan Region, Iraq
6
Department of Earth Sciences, Khalifa University of Science and Technology, Abu Dhabi P.O. Box 127788, United Arab Emirates
7
School of the Environment, University of Windsor, Windsor, ON N9B 3P4, Canada
8
Departament de Mineralogia, Petrologia i Geologia Aplicada, Facultat de Ciències de la Terra, Universitat de Barcelona, 08028 Barcelona, Spain
*
Author to whom correspondence should be addressed.
Water 2023, 15(10), 1953; https://doi.org/10.3390/w15101953
Submission received: 8 April 2023 / Revised: 14 May 2023 / Accepted: 17 May 2023 / Published: 21 May 2023

Abstract

:
Field observations, together with the results of gas chromatography–mass spectrometry (GC-MS) and stable carbon isotope analysis of bitumen, coupled with fluid inclusion microthermometry and stable isotope analyses of closely associated vug- and fracture-filling columnar calcite in the Upper Cretaceous Bekhme Formation, Kurdistan Region of Iraq, suggest that the degradation of crude oil was caused by the regional incursion of meteoric waters. This incursion, which is interpreted to have occurred during tectonic uplift during the Zagros Orogeny, is evidenced by: (i) the depletion of n-alkanes and acyclic isoprenoid alkanes (pristane and phytane) in the bitumen; (ii) low δ13CVPDB values (−8.5‰ to −3.9‰) and δ18OVPDB values (−22.9‰ to −15.0‰), with more radiogenic Sr isotopic ratios (0.70771–0.70772) compared to Cretaceous seawater; and (iii) low salinity and low temperatures (20 to 40 °C) in fluid inclusions of the columnar calcite. This study demonstrates that regional meteoric water incursion into sedimentary basins can be linked to crude oil degradation accompanied by calcite cementation events in carbonate reservoirs.

1. Introduction

Heavy oil and solid bitumen in hydrocarbon reservoirs create barriers to fluid flow, causing porosity and permeability deterioration and, thus, serious production problems [1]. Undersaturated meteoric waters preferentially remove water-soluble compounds from oil, such as light aromatics (benzene and toluene) [2], which results in an increase of bitumen and heavy oil residues [3,4]. This process is known as water washing [5]. Moreover, aerobic bacteria introduced by meteoric waters have also been suggested to cause the biodegradation of crude oils [6,7]. Water washing and biodegradation may act independently or together [8], but current analytical techniques are not able to precisely recognize their separate roles in the degradation of crude oil [6,7,9]. Apart from the role of meteoric waters, degassing of methane may also result in an enrichment of the asphaltene content in crude oil [10].
Heavy oils and solid bitumens with low API gravity (<27°) have frequently been encountered in recent discoveries and drilling activities in the Kurdistan Region of Iraq, but their origin is poorly understood [11]. It is not clear yet whether crude oil degradation and the formation of heavy oil is local or distributed across the entire Zagros Foreland Basin [7,11,12]. Heavy oil in several recent discoveries in the Bekhme Formation poses a serious challenge to the development of these oil fields in the region [11,13].
This study focuses on identifying the timing and origin of bitumen formation and the closely associated calcite cementation in the Upper Cretaceous Bekhme Formation at an outcrop in the Bekhme anticline in Iraqi Kurdistan. The Formation consists of dolostones containing vugs and fractures filled with coarse crystalline dolomite, calcite, and bitumen. This current study attempts to shed light on the degradation of hydrocarbons, generation of heavy oil, and the associated calcite cementation in the Upper Cretaceous Bekhme Formation and other equivalent reservoirs in the Zagros Foreland Basin.

2. Geological Setting

The study area in the Kurdistan Region of Iraq lies in the Bekhme anticline, which consists of various outcropping units ranging in age from Jurassic–Cretaceous to Miocene (Figure 1 and Figure 2A, [14,15]). The Bekhme Formation belongs to the Upper Campanian–Maastrichtian cycle. This depositional cycle began with a widespread transgression covering all of what is now Iraq, which took place after the termination of the tectonic activities caused by the Middle Cretaceous Orogeny. The cycle also resulted in several formations that show considerable lithological changes both vertically and horizontally, representing shallow to deep marine deposits [16,17].
The tectonic history of the region started by rifting in the Permo-Triassic, which evolved into a passive margin in the Middle Jurassic period and into an active margin in the Late Cretaceous period, owing to the onset of Arabia–Eurasia convergence and formation of the Zagros and Taurus orogenic belts [18,19,20]. A foreland basin developed during the Campanian–Maastrichtian period in northern Iraq along the Zagros margin, with a clastic depocenter represented by the Tanjero Formation, whereas the Bekhme Formation represents an isolated carbonate platform developed within a clastic starved part of the basin [21]. The Bekhme Formation lies within the tectonostratigraphic mega sequence AP9 K170 [22], is relatively short-lived (29 my), and contains an important phase of plate boundary development [17]. The onset of plate boundary development is marked by the major mid-Cretaceous uplift, regression, and associated unconformity, which is recognized across the entire Arabian Plate [22]. The Bekhme Formation was deposited during the Maastrichtian sea-level rise (K 180 MFS of [22]) in the Campanian platform.
The studied Bekhme Formation outcrop is about 315 m in height and is located in the type locality of the formation at the Bekhme Gorge (latitude 36°41′45″; longitude 44°16′ 57″). The Bekhme Formation in the studied area comprises three units: upper, middle, and lower. The lower unit (10 m thick) is composed of limestones and breccias. The middle unit (94 m thick) comprises limestones containing rudest fragments that alternate with fore-reef shoal limestones (Figure 2 and Figure 3). The upper unit is 211 m thick and is composed of bituminous dolostones and limestones [23].

3. Samples and Methods

Bitumen samples were collected from an outcrop in the Bekhme Gorge near the type locality of the Bekhme Formation (Figure 1 and Figure 2B, C). Bitumen occurs in vugs, where it engulfs saddle dolomite and equant calcite, and in large veins, where it engulfs columnar calcite. (Figure 3A–C). The bitumen samples were dissolved in dichloromethane and then deasphalted with n-hexane using 0.45 μm PFTE filters. The dissolved sample was fractioned into saturates, aromatics, resins, and asphaltenes using the SARA method [4].
Carbon isotopic composition of whole bitumen and SARA fractions was performed using a Thermo Finnigan1112 elemental analyzer coupled to a Finnigan Mat Delta-C mass spectrometer. The reference materials were USGS-40 L-glutamic acid, IAEA-CH6 sucrose, IAEA-CH7 polyethylene, IAEA 600 caffeine, and acetanilide. The 13C/12C ratio is reported in “δ” notation, and δ13C relative to VPDB (Vienna Pee Dee Belemnite) (Table 1). Precision of analysis was better than 0.1 per mil based on internal and international standards [24]. Saturated and aromatic hydrocarbons were analyzed by gas chromatography–mass spectrometry (GC-MS). The GC-MS analyses were performed on a Thermo Scientific Trace GC Ultra gas chromatograph with a DB-5 Agilent Technologies column (60 m × 0.25 mm i.d. × 0.1 mm film) coupled to an ITQ900 mass spectrometer. Samples (1 μL) were injected using a splitless injector heated at 280 °C, and helium was used as a carrier gas with a constant flow rate of 1 mL/min. The oven temperature was programmed to run from an initial temperature of 40 °C (held for 1 min) to a final temperature of 300 °C at 2 °C/min and then held at 300 °C for 60 min. The mass spectrometer was operated in selective ion monitoring (SIM) electron impact mode (electron input energy, 40 eV; source temperature, 200 °C).
The following ions were monitored: m/z 57 for n-alkanes, m/z 191 for tricyclic terpanes and hopanes, m/z 217 for steranes, and m/z 231 for triaromatic steroids [10,25,26]. MSD chemstation software was used for the identification of the biomarkers. Individual components were identified by comparison of their retention times and mass spectra with published data [10,24,25]. Relative abundances of triterpanes and steranes were calculated by measuring peak heights in the m/z 191 and m/z 217 fragmentograms, respectively. The extraction procedure and GC-MS analyses were carried out at the Department of Mineralogy, Petrology, and Applied Geology, University of Barcelona.
Twenty-three samples of columnar calcite were subjected to carbon and oxygen isotope analyses at the School of Environment, University of Windsor. The samples were reacted with pure H3PO4 at 25 °C for 1 h, and isotope analyses were then carried out on a Thermo Finnigan Delta Plus mass spectrometer, the analytical precision of which was ± 0.05‰. Isotopic values are given in δ (delta)-notation and reported in per mil (‰) relative to the Vienna Pee Dee Belemnite (VPDB) [27].
The 87Sr/86Sr ratios were obtained from columnar calcite powders by leaching the samples in 0.2 N HCL for 24 h at room temperature. The isotopic measurements were performed using an automated Finnigan 261 mass spectrometer equipped with nine Faraday collectors. Correction for isotopic fractionation during the analyses was made by normalization to 86Sr/88Sr = 0.1194. The mean standard error of mass spectrometer performance was ± 0.00003 for standard NBS-987. Fluid inclusion microthermometry of calcite was performed on two double-polished wafers. A petrographic examination of fluid inclusions was conducted in order to determine the origin (primary, secondary, pseudo-secondary, or uncertain), composition (“aqueous” or “hydrocarbon”), and estimated liquid to vapor ratios.

4. Results

Coarse crystalline columnar calcite lines large fractures filled by bitumen (Figure 3A–C).

4.1. Composition of Bitumen

SARA results show small quantities of saturated and aromatic hydrocarbons (5% and 11%, respectively) and substantial amounts of heavy compounds (i.e., resins and asphaltenes, 7.9% and 66.6%, respectively, and a loss of 9.3%). The mass chromatograms of TIC (total ion current) (Figure 4A) and ion m/z 57 show an UCM (unresolved complex mixture) hump (Figure 4B) and the m/z 191 chromatogram (Figure 5) shows the distributions of terpanes (tricyclics and pentacyclics). The presence of C24 tetracyclic terpane is noticeable. The quantity of trisnorneohopane (Ts) is extremely low, compared with trisnorhopane ™, and consequently, the calculated Ts/(Ts + Tm) ratio is low too (0.15) (Figure 5A). Hopanes are dominated by C29H. C29 norhopane is, overall, more abundant than C30 hopane (Figure 5B), (C29H/C30H ratio = 1.62) and the C35 homohopane shows a similar abundance as C34 homohopane. The steroids (m/z 217 ion) chromatograms (Figure 6) reveal that steranes and diasteranes are present in relatively high quantities in the saturated hydrocarbon fractions of the studied sample. The C27 to C29 diasteranes are also present in low concentrations. Low-molecular-weight steranes (e.g., C21 pregnane and C22 homopregnane) are present in high amounts and pregnane is more prominent than the homopregnane (Figure 6).
The bitumen contains higher amounts of C28-TA (20R, 20S) than C27-TA (20R) and C26-TA (20S) counterparts and the value of C28-TA (20S + 20R)/(C26-TA 20S + C27-TA 20R) is equal 1.49 (Figure 7). The δ13CVPDB values of the whole fraction of bitumen vary between −27.9‰ and −27.8‰.

4.2. δ18. O, 87Sr/86Sr Isotopes and Fluid Inclusions of Columnar Calcite

The vugs (mostly 5–15 cm across) and breccia fractures are filled by coarse crystalline saddle dolomite, equant calcite, and, less frequently, anhydrite. The thrust planes are lined with large crystals of columnar calcite (ca. 3–5 cm long and 1 cm across). These crystals are oriented perpendicular to the thrust plane surfaces.
Columnar calcite has low δ18OVPDB (−22.9‰ to −15.0‰), δ13CVPDB (−9.8‰ and −2.5‰), and 87Sr/86Sr ratios of 0.70771 to 0.70772 (Figure 8 and Figure 9, [28,29]). Mansurbeg et al. [30] reported that: (i) saddle dolomite has narrow ranges of δ18OVPDB (−14.5‰ to −7.4‰), δ13CVPDB (−3.6‰ to +2.6‰), and Sr isotopic ratios (0.70774–0.70780), and (ii) vug/fracture-filling equant calcite cement, which engulfs saddle dolomite, displays wide ranges of δ13CVPDB (−10.5‰ to +2.1‰), δ18OVPDB (−19.2‰ to −4.8‰), and Sr isotopic ratios (0.70773–0.70796).
Fluid inclusions in the columnar calcite crystals are mostly single-phase, liquid type and range in size from about 7 to 25 μm. Final ice melting temperatures (Tmice) (cf. [31]) suggest average salinity of 0.5 wt.% NaCl eq. Furthermore, Mansurbeg et al. [30] reported that the homogenization temperatures (Th) in vug-filling saddle dolomite are 80–120 °C and first ice melting temperatures (Te) are −38 °C to −42 °C. Measured Tmice points to the salinity of 17–22 wt.% NaCl eq. The Th in equant calcite ranges between about 60 °C and 100 °C and the Te between about −50 to −52 °C. Salinity derived from Tmice ranges from 14 to 18 wt.% NaCl eq.
UV light observations made on hydrothermal equant calcite in the Bekhme Formation showed the presence of oil-filled fluid inclusions with pale- to dull-yellow or orange fluorescence [30]. No oil-filled fluid inclusions were reported in any other cements, except for equant calcite [28], which is typical of oil with an API gravity of about 15–25° (cf. [32]).

5. Discussion

5.1. Bulk Bitumen Composition: Clues to Its Origin

The low contents of saturate and aromatic hydrocarbons (5% and 11% respectively) and a large percentage of heavy compounds (resins and asphaltenes) indicate that bitumen consists of heavy or degraded oil [33]. The carbon isotopic values for saturate and aromatic hydrocarbons (−27.8 ‰) suggest a marine origin [34]. The low concentration of C23 to C29 tricyclic terpanes suggests a limited input of terrigenous organic material [35,36]. The presence of C24 tetracyclic terpane is expected for carbonate/evaporitic depositional environments (e.g., [10,37]).
The UCM hump revealed by mass chromatograms of TIC and ion m/z 57 (Figure 4) provides clues to the biodegradation processes. The lack of n-alkanes and isoprenoids, which prevents the determination of the carbon preference index (CPI), as well as ratios of pristane Pr/n-C17 and phytane Ph/n-C18, is also interpreted be the result of biodegradation (e.g., [33]). The prominent peaks present in the TIC and ion m/z 57 are interpreted to represent terpane and sterane biomarkers, which are relatively resistant to biodegradation. The near-absence of tricyclic terpanes (represented only by the C23/3) on the m/z 191 chromatogram (Figure 5) cannot be attributed to biodegradation because these compounds are very resistant to biodegradation [38,39].
The extremely small amounts of trisnorneohopane (Ts) compared to trisnorhopane ™ and the Ts/(Ts + Tm) ratio (0.15) indicate low thermal maturity due to the lower thermal stability of C2717α(H)-22,29,30-trisnorhop™(Tm) compared to that of C2718α(H)-22,29,30-trisnorneohopane (Ts). Nevertheless, the Tm and Ts contents are influenced by the degree of maturation, type of organic matter, and lithology (e.g., [40,41]. According to Peters et al. [10], oils from carbonate source rocks show much lower (< 0.25) ratios compared to those in shales deposited in anoxic environments (0.26–0.35).
This ratio is an indicator for the type of carbonate rock. Carbonates have less than 0.25, and shale more than 0.26. The ratio in the studied sample is 0.15, so the source rock was a carbonate rock.
The presence of hopanes, particularly C29H and the relative abundance of C29 norhopane compared to C30 hopane, (C29H/C30H ratio of 1.62) are typical for carbonate-rich source rocks [42,43]. The presence of gammacerane (Figure 5) suggests the restricted environment with water column stratification during source rock deposition [10,41,44]. This interpretation is further supported by the domination of homohopanes by the C31 homohopane and its decrease with increasing carbon number (Figure 5). The similar amounts of C35 homohopane and C34 homohopane and the predominance of C29/C30 and the C35/C34 hopanes are related to carbonate (hypersaline) environments (e.g., [41,45].
There are relatively high amounts of steranes and diasteranes (Figure 6) in the saturated hydrocarbon fractions of the studied sample. The observed higher proportion of C29 compared with C27 and C28 regular steranes suggests a terrigenous organic matter contribution to the source rock [46,47]. However, partial biodegradation of C27 and C28 steranes cannot be ruled out and, consequently, the ratio C29 over C27 cannot be taken into account. However, the presence of C30 sterane is a characteristic of marine origin. The presence of both C30 sterane and low Ts/Tm ratios (<1) further suggest that these oils were probably derived from a carbonate source rock deposited in a reducing environment [48,49]. The low concentrations of C27 to C29 diasteranes compared to low-molecular-weight steranes, such as C21 pregnane and C22 homopregnane, and the more prominent pregnane than the homopregnane (Figure 6) further support sediments deposited under hypersaline conditions [50].
The small concentrations of aromatic hydrocarbons are attributed to biodegradation and probably water washing. Conversely, the presence of triaromatic hydrocarbons (TA) is due to high resistivity to biodegradation under most depositional and diagenetic conditions [51,52,53]. The presence of C20 and C21 components points to low-to-medium biodegradation [54]. Due to the co-elution of C26-TA (R) with C27-TA 20(S), four other isomers were used to differentiate homologue variations. There are higher amounts of C28-TA (20R, 20S) than C27-TA (20R) and C26-TA (20S) counterparts and the value of C28-TA (20S + 20R)/(C26-TA 20S + C27-TA 20R) is equal to 1.49 (Figure 7).

5.2. Origin of Columnar Calcite

The overall trend of positive correlation between δ18OVPDB and δ13CVPDB values in the columnar calcite cement (Figure 8) indicates the input of 12C upon a decrease in δ18OVSMOW of the pore fluids. The latter isotopic trend might suggest a trend of increased incursion rates of meteoric waters (cf. [30,55]. The higher Sr isotopic ratios of the columnar calcite than the range for Cretaceous seawater (0.7073–0.7075; [29]; Figure 9) corroborates input of 87Sr by the incursion of meteoric waters. The source of this elevated radiogenic Sr can be related to the input of 87Sr from the circulation of meteoric waters in siliciclastic rocks [30].
Integrating field observations, burial history curve, stable isotope analyses, and fluid inclusion microthermometry allows us to decipher the origin of calcite that is closely associated with bitumen in the Upper Cretaceous Bekhme Formation. Evidence shows that bitumen postdates the crystallization of saddle dolomite [30], whereas the columnar calcite covers bitumen and field observations indicate that the bitumen is intermixed with equant calcite and fills the center of vugs and fractures.
The negative δ18OVPDB values (−22.9‰ to −15.0‰) and inferred temperatures (≤50 °C), coupled with low salinity (average 0.5 wt.% NaCl eq.) of the columnar calcite suggest precipitation from meteoric waters [30]. Using these oxygen isotopic values and near-surface temperatures (25–40 °C) (Figure 10), and the oxygen isotope fractionation equation between calcite and water [56], the δ18OVSMOW of the meteoric water is inferred to range from −12.1‰ to −22.4‰. However, for poorly understood reasons, these values are too low for the reported meteoric waters at the paleo-latitude of the basin, which are similar to present-day latitudes of δ18OVSMOW of around −5‰, e.g., [57]. The δ18OVPDB of calcite is instead −7.2‰. However, local factors, such as relief and seasonal variations, may affect δ18O values (cf. [30,57]).
The wide range of δ13CVPDB values (−9.8‰ and −2.5‰) of the columnar calcite cement is attributed to variations in the input of 12C-enriched dissolved carbon derived from the degradation of crude oil δ13CVPDB values <−20‰; [58,59,60]), as well as input of 13C from dissolution of Cretaceous marine carbonates. Therefore, calcite with higher δ13CVPDB values was likely formed during early stages of crude oil degradation, unlike calcite with lower δ13CVPDB values.
The estimated API of the trapped hydrocarbons within the equant calcite (15–25° API; [30], which is paragenetically earlier than the columnar calcite (Figure 11), is higher than the measured API (10–17°; [61]) in the present-day oil, which further supports the role of meteoric water biodegradation.
The burial history curve for the study area (Figure 12; [11,62]) shows a significant uplift/erosion of the Bekhme Formation in the final stages of Zagros Orogeny, bringing the successions close to the basin margin (i.e., with increased hydraulic head that facilitated meteoric water incursion). The presence of the columnar calcite (Figure 3A–C) along the thrust walls suggests that cool (<50 °C), low-salinity (<0.5 wt.% NaCl eq.) meteoric-water flux, coupled with oil degradation and columnar calcite cementation, occurred simultaneously. Oil degradation is a process that generally increases carbonate alkalinity and, hence, promotes carbonate precipitation, as evidenced by the isotopically light δ13C isotopic composition of this calcite [30,60].

6. Conclusions

Integration of gas chromatography, fluid inclusion petrography and microthermometry, carbon, oxygen, and strontium isotope data has provided important insights into the role of meteoric water incursion in crude oil degradation in Cretaceous carbonate reservoirs in the Kurdistan Region of Iraq. The main conclusions derived from this study include:
  • Alteration of hydrocarbons took place via biodegradation and water washing that affected aromatics and steranes. The organic matter precursor was deposited in a restricted carbonatic/evaporitic environment.
  • The depletion of n-alkanes and acyclic isoprenoid (pristane and phytane) confirm that biodegradation affected the hydrocarbons. The prominent peaks are mostly terpane and sterane biomarkers, which are more resistant to the biodegradation. The presence of triaromatic steroid (TAS) hydrocarbons in the analyzed sample is evidence of moderate degradation. Moreover, the presence of C2628 homologues and the C20 and C21 components of TAS are consistent with the interpretation of limited biodegradation of the hydrocarbons.
  • High quantities of steranes and diasteranes in the saturated hydrocarbon fraction and a higher proportion of C29 compared with C27 and C28 regular steranes are attributed to terrigenous organic matter contribution.
  • The negative oxygen and carbon stable isotopes, non-radiogenic Sr isotopic ratios, low salinity, and low temperatures in fluid inclusions of the columnar calcite suggest that calcite cementation occurred simultaneously with oil degradation, which acted as the source of 12C in dissolved carbon needed for calcite precipitation.

Author Contributions

Methodology, I.M. and A.P.; Investigation, N.S., M.A., I.A.-A. and R.M.; Resources, S.S.; Writing—original draft, H.M.; Writing—review & editing, I.M., I.A.-A., A.P. and A.A.-J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

(ISA) would like to express its gratitude to the Natural Science and Engineering Research Council of Canada (NSERC) (812005) for its support. We would also like to extend our sincere appreciation to the General Directorate of Scientific Research Center, Salahaddin University-Erbil, in the Kurdistan Region of Iraq. We appreciate the insightful comments and suggestions from two anonymous reviewers, which have significantly enhanced the quality of the manuscript. Furthermore, we would like to acknowledge the invaluable assistance of editors at MDPI.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Veil, J.; Quinn, J. Water Issues Associated with Heavy Oil Production; Environmental Science Division, Argonne National Laboratory: Lemont, IL, USA, 2008. [CrossRef]
  2. Pytlak, A.; Leis, A.; Prochaska, W.; Sachsenhofer, R.F.; Gross, D.; Linzer, H.G. Light Hydrocarbon Geochemistry of Oils in the Alpine Foreland Basin: Impact of Geothermal Fluids on the Petroleum System. Geofluids 2017, 2017, 7182959. [Google Scholar] [CrossRef]
  3. Mason, P.C.; Burwood, R.; Mycke, B. The reservoir geochemistry and petroleum charging histories of Palaeogene-reservoired fields in the Outer Witch Ground Graben. Geol. Soc. Lond. Spec. Publ. 1995, 86, 281–301. [Google Scholar] [CrossRef]
  4. Bennett, B.; Larter, S.R. Quantitative separation of aliphatic and aromatic hydrocarbons using silver ion-silica solid-phase Extraction. Anal. Chem. 2000, 72, 1039–1044. [Google Scholar] [CrossRef] [PubMed]
  5. Mode, A.W.; Anyiam, O.A.; Amobi, J.O.; Nweke, S.U. Gas chromatographic analysis of whole oil samples: Implications for biodegradation in the Niger Delta. J. Petrol. Explor. Prod. Technol. 2017, 7, 389–398. [Google Scholar] [CrossRef]
  6. Blanc, P.H.; Connan, J. Crude oils in reservoirs: The factors influencing their composition. In Applied Petroleum Geochemistry; Bordenave, M.L., Ed.; Editions Technips: Paris, France, 1993; pp. 151–173. [Google Scholar]
  7. Shlimon, A.G.; Mansurbeg, H.; Othman, R.S.; Gittel, A.; Aitken, C.M.; Head, I.M.; Finster, K.W.; Kjeldsen, K.U. Microbial community composition in crude oils and asphalts from the Kurdistan Region of Iraq. Geomicrobiol. J. 2020, 37, 635–652. [Google Scholar] [CrossRef]
  8. Jobson, A.M.; Cook, F.D.; Westlake, D.W.S. Microbial Utilization of Crude Oil. Appl. Microbiol. 1972, 23, 1082–1089. [Google Scholar] [CrossRef] [PubMed]
  9. Shlimon, A.; Mansurbeg, H.; Othman, R.; Head, I.; Kjeldsen, K.U.; Finster, K. Identity and hydrocarbon degradation activity of enriched microorganisms from natural oil and asphalt seeps in the Kurdistan Region of Iraq (KRI). Biodegradation 2021, 32, 251–271. [Google Scholar] [CrossRef] [PubMed]
  10. Peters, K.E.; Walters, C.C.; Moldowan, J.M. The Biomarker Guide, Vol. 2: Biomarkers and Isotopes in Petroleum Systems and Earth History, 2nd ed.; Cambridge University Press: Cambridge, UK, 2005; p. 1132. [Google Scholar]
  11. English, J.; Davies, L. Regional trends in Jurassic and Triassic thermal maturity in northern Iraq. In Proceedings of the Geological Society of London Conference on Hydrocarbon Exploration in the Zagros Mountains of Iraqi Kurdistan and Iran, London, UK, 23–25 January 2013. [Google Scholar]
  12. Mackertich, D.S.; Samarrai, A.I. History of hydrocarbon exploration in the Kurdistan Region of Iraq. GeoAarabia 2015, 20, 181–220. [Google Scholar] [CrossRef]
  13. EL Guindy, M.; Kelly, S.P. September. Logging While Drilling in Kurdistan’s First Horizontal Well Enhances Production Rates. In Second EAGE Workshop on Iraq; European Association of Geoscientists & Engineers: Utrecht, The Netherlands, 2013; p. cp-355. [Google Scholar]
  14. Csontos, L.; Sasvári, Á.; Pocsai, T.; Kósa, L.; Salae, A.T.; Ali, A. Structural evolution of the northwestern Zagros, Kurdistan Region, Iraq: Implications on oil migration. GeoArabia 2012, 17, 81–116. [Google Scholar]
  15. Mansurbeg, H.; Alsuwaidi, M.; Salih, N.; Shahrokhi, S.; Morad, S. Integration of stable isotopes, radiometric dating and microthermometry of saddle dolomite and host dolostones (Cretaceous carbonates, Kurdistan, Iraq): New insights into hydrothermal dolomitization. Mar. Pet. Geol. 2021, 127, 104989. [Google Scholar] [CrossRef]
  16. Ali, D.O. Sedimentology and Stratigraphy of Bekhme Formation (Upper Cretaceous) in Selected Sections in Kurdistan Region-Iraq. Master’s Thesis, Salahaddin University, Erbil, Iraq, 2005. [Google Scholar]
  17. Jassim, S.Z.; Goff, J.C. Geology of Iraq; Dolin, Prague and Moravian Museum: Brno, Czech Republic, 2006; p. 341. [Google Scholar]
  18. Stöcklin, J. Salt deposits of the Middle East. In Saline Deposits; Mattox, R.B., Ed.; Geological Society of America Special Paper 88; Geological Society of America: Boulder, CO, USA, 1968; pp. 157–181. [Google Scholar]
  19. Stoneley, R. A review of petroleum source rocks in parts of the Middle East. Geol. Soc. Lond. Spec. Publ. 1987, 26, 263–269. [Google Scholar] [CrossRef]
  20. Jackson, J.; Mckenzie, D. Active tectonics of the Alpine-Himalayan Belt between western Turkey and Pakistan. Geophys. J. Int. 1984, 77, 185–264. [Google Scholar] [CrossRef]
  21. Aqrawi, A.A.; Goff, J.C.; Horbury, A.D.; Sadooni, F.N. The Petroleum Geology of Iraq; Scientific Press: Beaconsfield, UK, 2010; p. 424. [Google Scholar]
  22. Sharland, P.R.; Archer, R.; Casey, D.M.; Davies, R.B.; Hall, S.; Heward, A.P.; Horbury, A.D.; Simmons, M.D. Arabian Plate Sequence Stratigraphy; Gulf Petrolink: Manama, Bahrain, 2001; p. 371. [Google Scholar]
  23. Wetzel, R. Stratigraphy Section on Bekhme Gorge; MPC Report, IR/Rw/4 No. 19; INOC Library: Albany, NY, USA, 1946. [Google Scholar]
  24. Permanyer, A.; Marfil, R.; Martí-Martín, J.D.; Estupiñan, J.; Márquez, G.; Arroyo, X. Diagenetic evolution of the Upper Cretaceous carbonate and sandstone exhumed reservoirs (Western Basque-Cantabrian Basin, North Spain). J. Pet. Geol. 2015, 66, 673–694. [Google Scholar] [CrossRef]
  25. Philp, R.P. Biological markers in fossil fuel production. Mass Spectrom. Rev. 1985, 4, 1–54. [Google Scholar] [CrossRef]
  26. Mohialdeen, I.M.J. Source Rock Appraisal and Oil-Source Correlation for the Chia Gara Formation, Kurdistan-North Iraq Unpublished. Ph.D. Thesis, College of Science, University of Sulaimani, Sulaymaniyah, Iraq, 2008. [Google Scholar]
  27. Al-Aasm, I.S.; Taylor, B.E. South, Stable isotope analysis of multiple carbonate samples using selective aciextraction. Chem. Geol. 1990, 80, 119–125. [Google Scholar]
  28. Veizer, J.; Ala, D.; Azmy, K.; Brucksch, P.; Buhl, D.; Bruhn, F.; Carden, G.A.F.; Diener, A.; Ebneth, S.; Godde’ris, Y.; et al. 87Sr/86Sr, d13C and d18O evolution of Phanerozoic seawater. Chem. Geol. 1999, 161, 59–88. [Google Scholar] [CrossRef]
  29. McArthur, J.M.; Donovan, D.T.; Thirlwall, M.F.; Fouke, B.W.; Mattey, D. Strontium isotope profile of the early Toarcian (Jurassic) oceanic anoxic event, the duration of ammonite biozones, and belemnite palaeotemperatures. Earth Planet. Sci. Lett. 2000, 179, 269–285. [Google Scholar] [CrossRef]
  30. Mansurbeg, H.; Morad, D.; Othman, R.; Morad, S.; Ceriani, A.; Al-Aasm, I.; Kolo, K.; Spirov, P.; Proust, N.; Preat, A.; et al. Hydrothermal dolomitization of the Bekhme Formation (Upper Cretaceous), Zagros Basin, Kurdistan Region of Iraq: Record of oil migration. Sediemntary Geol. 2016, 341, 147–162. [Google Scholar] [CrossRef]
  31. Goldstein, H.; Reynolds, J. Systematics of Fluid Inclusions in Diagenetic Minerals; Short Course 31; SEPM Society for Sedimentary Geology: Tulsa, OK, USA, 1994. [Google Scholar]
  32. Mclimans, R.K. The application of fluid inclusions to migration of oil and diagenesis in petroleum reservoirs. Appl. Geochem. 1987, 2, 585–603. [Google Scholar] [CrossRef]
  33. Tissot, B.P.; Welte, D.H. Petroleum Formation and Occurrence; Springer: Berlin, Germany, 1984; p. 538. [Google Scholar]
  34. Sofer, Z. Stable Carbon Isotope composition of crude oils: Application to source depositional environments and petroleum alteration. AAPG Bull. 1984, 68, 31–49. [Google Scholar]
  35. Hunt, J.M. Petroleum Geochemistry and Geology; W.H. Freeman and Company: San Francisco, CA, USA, 1996; 743p. [Google Scholar]
  36. Leila, M.; Moscariello, A. Organic geochemistryof oil and natural gasin the west Dikirnis and El-Tamad Fields, onshore Nile delta, Egypt: Interpretation of potential source rocks. J. Pet. Geol. 2017, 40, 37–58. [Google Scholar] [CrossRef]
  37. Alkhafaji, M.W.; Connan, J.; Engel, M.H.; Al-Jubouri, S.W. Origin, biodegradation, and water washing of bitumen from the Mishraq Sulfur Mine, northern Iraq. Mar. Pet. Geol. 2020, 124, 104786. [Google Scholar] [CrossRef]
  38. Behar, F.H.; Albrecht, P. Correlations between carboxylic acids and hydrocarbons in several crude oils. Alteration by biodegradation. Org. Geochem. 1984, 6, 597–604. [Google Scholar] [CrossRef]
  39. Yang, F.; Wang, T.; Li, M. The Distribution of Triaromatic Steroids and Oil Group Classification of Ordovician Petroleum Systems in the Cratonic Region of the Tarim Basin, NW China. Pet. Sci. Technol. 2015, 33, 1794–1800. [Google Scholar] [CrossRef]
  40. Seifert, W.K.; Moldowan, J.M. The effect of biodegradation on steranes and terpanes in crude oils. Geochim. Cosmochim. Acta 1979, 43, 111–126. [Google Scholar] [CrossRef]
  41. Moldowan, J.M.; Seifert, W.K.; Gallegos, E.J. Relationship between petroleum Compositions and depositional environment of Petroleum source rocks. AAPG Bull. 1985, 69, 1255–1268. [Google Scholar]
  42. Waples, D.W.; Machihara, T. Biomarkers for geologists: A practical guide to theI application of steranes and triterpanes in petroleum geology. AAPG Meth. Exp. 1991, 2, 5–10. [Google Scholar]
  43. Mohialdeen, I.M.J.; Hakimi, M.H. Geochemical characterization of Tithonian- Berriasian Chia Gara organic-rich rocks in northern Iraq with an emphasis on organic matter enrichment and the relationship to the bioproductivity and anoxia conditions. J. Asian Earth Sci. 2016, 116, 181–197. [Google Scholar] [CrossRef]
  44. Damsté, J.S.S.; Kenig, F.; Koopmans, M.P.; Köster, J.; Schouten, S.; Hayes, J.M.; de Leeuw, J.W. Evidence for gammacerane as an indicator of water column stratification. Geochim. Cosmochim. Acta 1995, 59, 1895–1900. [Google Scholar] [CrossRef]
  45. Grantham, P.J. The occurence of unusual C27 and C29 sterane predominances in two types of Oman crude oil. Org. Geochem. 1986, 9, 1–10. [Google Scholar] [CrossRef]
  46. Huang, W.Y.; Meinschein, W.G. Sterols as ecological indicators. Geochim. Cosmochim. Acta 1979, 43, 739–745. [Google Scholar] [CrossRef]
  47. Peters, K.E.; Moldowan, J.M. The Biomarker Guide: Interpreting Molecular Fossils in Petroleum and Ancient Sediments, 1st ed.; Prentice Hall: Hoboken, NJ, USA, 1993; p. 363. [Google Scholar]
  48. Mckirdy, D.M.; Kantsler, A.J.; Emmett, J.K.; Aldridge, A.K. Hydrocarbon genesis and organic facies in Cambrian carbonates of the Eastern Officer Basin, South Australia. In Petroleum Geochemistry and Source Rock Potential of Carbonate Rocks; Palacas, J.G., Ed.; AAPG: Tulsa, OK, USA, 1984; pp. 13–31. [Google Scholar]
  49. Rullkötter, J.; Spiro, B.; Nissenbaum, A. Biological marker characteristics of oils and asphalts from carbonate source rocks in a rapidly subsiding graben, Dead Sea, Israel. Geochim. Cosmochim. Acta 1985, 49, 1357–1370. [Google Scholar] [CrossRef]
  50. Ten Haven, H.L.; De Leeuw, J.W.; Rullkötter, J.; Sinninghe Damsté, J.S.; Schenck, P.A.; Palmer, S.E.; Zumberge, J.E.; Fleet, A.J.; Kelts, K.; Talbot, M.R. Application of biological markers in the recognition of palaeohypersaline environments. Lacustrine Petroleum Source Rocks. Geol. Soc. 1988, 40, 123–130. [Google Scholar] [CrossRef]
  51. Lin, L.H.; Michael, G.E.; Kovachev, G.; Zhu, H.; Philp, R.P.; Lewis, C.A. Biodegradation of tar-sands bitumens from the Ardmore and Anadarko Basins, Carter Country, Oklahoma. Org. Geochem. 1989, 14, 511–523. [Google Scholar] [CrossRef]
  52. Yang, C.; Wang, Z.; Liu, N.; Yang, Z.; Li, Y.; Shah, K.; Zhang, G.; Landriault, M.; Hollebone, B.; Brown, C.; et al. Aromatic Steroids in Crude Oils and Petroleum Products and Their Applications in Forensic Oil Spill Identification. Environ. Forensics 2013, 14, 278–293. [Google Scholar] [CrossRef]
  53. Li, F.J.; Qu, X.L.; Du, L.C.; Jing, X.G.; Cheng, X.Y. Carbonate facies and sedimentary evolution of Upper Devonian Tuqiaozi formation in Ganxi section of Longmen mountain. Lithol. Reserv. 2015, 27, 6–12. [Google Scholar]
  54. Volkman, J.K.; Alexander, R.; Kagi, R.I.; Woodhouse, G.W. Demethylated hopanes in crude oils and their applications in petroleum geochemistry. Geochim. Cosmochim. Acta 1983, 47, 783–794. [Google Scholar] [CrossRef]
  55. Swart, P.K. The geochemistry of carbonate diagenesis: The past, present and future. Sedimentology 2015, 62, 1233–1304. [Google Scholar] [CrossRef]
  56. Friedman, I.; O’neil, J.R. Compilation of stable isotope fractionation factors of geochemical interest. In Data of Geochemistry; Fleisher, M., Ed.; U.S. Geological Survey Professional Paper 440-KK; U.S. Geological Survey: Reston, VR, USA, 1977; p. 12. [Google Scholar]
  57. Clark, I.D.; Fritz, P. Environmental Isotopes in Hydrogeology; Lewis Publishers: Boca Raton, FL, USA, 1997; p. 342. [Google Scholar]
  58. Gluyas, J.; Swarbrick, R. Petroleum Geoscience; Blackwell Publishing: Oxford, UK, 2004; 349p. [Google Scholar]
  59. Waples, D.W. Geochemistry in Petroleum Exploration; International Human Resources Development Corporation: Boston, MA, USA, 1985; p. 232. [Google Scholar]
  60. Schulz-Rojahn, J.; Ryan-Grigor, S.; Anderson, A. Structural Controls on Seismic-Scale Carbonate Cementation in Hydrocarbon-Bearing Jurassic Fluvial and Marine Sandstones from Australia: A Comparison. In Carbonate Cementation in Sandstones; Morad, S., Ed.; Wiley: Hoboken, NJ, USA, 1998; pp. 327–362. [Google Scholar] [CrossRef]
  61. English, J.; Lunn, G.A.; Ferreira, L.; Yacu, G. Geologic evolution of the Iraqi Zagros, and its influence on the distribution of hydrocarbons in the Kurdistan region. AAPG Bull. 2015, 99, 231–272. [Google Scholar] [CrossRef]
  62. Gordon, G. Innovation and Technology Dinarta Block, Kurdistan Region of Iraq. In Proceedings of the Kurdistan-Iraq Oil and Gas Conference, London, UK, 16–18 December 2014; p. 35. [Google Scholar]
Figure 1. Geological map of northern Iraq (A) showing the distribution of the Bekhme Formation (green) and Bekhme Anticline (red line) in the Kurdistan Region (B) (after [14]).
Figure 1. Geological map of northern Iraq (A) showing the distribution of the Bekhme Formation (green) and Bekhme Anticline (red line) in the Kurdistan Region (B) (after [14]).
Water 15 01953 g001
Figure 2. (A) Chronostratigraphic column for the study area (after [15]); (B) sampling location in the Bekhme Gorge, Kurdistan Region of northern Iraq; (C) outcrop showing the thick bitumen seeping out of the Bekhme Formation.
Figure 2. (A) Chronostratigraphic column for the study area (after [15]); (B) sampling location in the Bekhme Gorge, Kurdistan Region of northern Iraq; (C) outcrop showing the thick bitumen seeping out of the Bekhme Formation.
Water 15 01953 g002
Figure 3. (AC) Field images showing bitumen accumulation and vertically arranged columnar calcite crystals (red arrows) occurring along a fracture and thrust surfaces in Bekhme Formation.
Figure 3. (AC) Field images showing bitumen accumulation and vertically arranged columnar calcite crystals (red arrows) occurring along a fracture and thrust surfaces in Bekhme Formation.
Water 15 01953 g003
Figure 4. (A) Total ion current and (B) m/z 57 chromatograms of saturated hydrocarbons showing a strong biodegradation of n-alkanes and an important unresolved complex mixture (UCM).
Figure 4. (A) Total ion current and (B) m/z 57 chromatograms of saturated hydrocarbons showing a strong biodegradation of n-alkanes and an important unresolved complex mixture (UCM).
Water 15 01953 g004
Figure 5. Partially reconstructed chromatogram (m/z 191) displaying the tricyclic terpanes (A) and homohopanes (B) distributions in the extracted sample.
Figure 5. Partially reconstructed chromatogram (m/z 191) displaying the tricyclic terpanes (A) and homohopanes (B) distributions in the extracted sample.
Water 15 01953 g005
Figure 6. Partially reconstructed chromatogram (m/z 217) showing the low-molecule pregnanes, regular steranes, and diasterane distributions in the extracted sample.
Figure 6. Partially reconstructed chromatogram (m/z 217) showing the low-molecule pregnanes, regular steranes, and diasterane distributions in the extracted sample.
Water 15 01953 g006
Figure 7. Partially reconstructed chromatogram (m/z 231) displaying the triaromatic steroids distribution in the extracted sample.
Figure 7. Partially reconstructed chromatogram (m/z 231) displaying the triaromatic steroids distribution in the extracted sample.
Water 15 01953 g007
Figure 8. Cross-plot of the carbon and oxygen isotopic values of columnar calcite. These values are lower than what would be expected for calcite that was precipitated in equilibrium with Cretaceous seawater, as indicated by the pink box [28]. Note the positive correlation between carbon and oxygen isotopic values of the columnar calcite if the three outliers with unusually low δ18O values are removed.
Figure 8. Cross-plot of the carbon and oxygen isotopic values of columnar calcite. These values are lower than what would be expected for calcite that was precipitated in equilibrium with Cretaceous seawater, as indicated by the pink box [28]. Note the positive correlation between carbon and oxygen isotopic values of the columnar calcite if the three outliers with unusually low δ18O values are removed.
Water 15 01953 g008
Figure 9. Cross-plot of strontium and oxygen isotopes [29] of columnar calcite showing that Sr isotopic ratios of the columnar calcite are higher than the postulated range for Cretaceous seawater.
Figure 9. Cross-plot of strontium and oxygen isotopes [29] of columnar calcite showing that Sr isotopic ratios of the columnar calcite are higher than the postulated range for Cretaceous seawater.
Water 15 01953 g009
Figure 10. Cross-plot of homogenization temperatures versus oxygen isotopic values of calcite and dolomite. Note that saddle dolomite, equant calcite, and host dolostones have overlapping values, presumably owing to the genetic link between these carbonates, which is the flux of hydrothermal brines [30]. The negative oxygen isotopic values combined with the very low salinity of fluids responsible for the precipitation of columnar calcite suggest a meteoric origin.
Figure 10. Cross-plot of homogenization temperatures versus oxygen isotopic values of calcite and dolomite. Note that saddle dolomite, equant calcite, and host dolostones have overlapping values, presumably owing to the genetic link between these carbonates, which is the flux of hydrothermal brines [30]. The negative oxygen isotopic values combined with the very low salinity of fluids responsible for the precipitation of columnar calcite suggest a meteoric origin.
Water 15 01953 g010
Figure 11. General paragenetic sequence reconstructed for the Bekhme Formation in relation to burial and orogenic events.
Figure 11. General paragenetic sequence reconstructed for the Bekhme Formation in relation to burial and orogenic events.
Water 15 01953 g011
Figure 12. Burial history curve of the study area shows a significant uplift/erosion of the Bekhme Formation in the Upper Cretaceous period (modified after [11,15,62]).
Figure 12. Burial history curve of the study area shows a significant uplift/erosion of the Bekhme Formation in the Upper Cretaceous period (modified after [11,15,62]).
Water 15 01953 g012
Table 1. Isotopic values for whole bitumen and all SARA fractions.
Table 1. Isotopic values for whole bitumen and all SARA fractions.
Samplesδ13CVPDB (‰)Fraction
SM-1 BIT−27.8Whole bitumen
SM-1 ASP−27.9Asphaltenes
SM-1 RES−27.8Resins
SM-1 ARO−27.8Aromatics
SM-1 SAT−27.9Saturates
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Mansurbeg, H.; Mohialdeen, I.; Al-Juboury, A.; Salih, N.; Alsuwaidi, M.; Shahrokhi, S.; Al-Aasm, I.; Mahmmud, R.; Permanyer, A. Meteoric Water Incursion, Crude Oil Degradation and Calcite Cementation of an Upper Cretaceous Reservoir in the Zagros Foreland Basin (Kurdistan Region of Iraq). Water 2023, 15, 1953. https://doi.org/10.3390/w15101953

AMA Style

Mansurbeg H, Mohialdeen I, Al-Juboury A, Salih N, Alsuwaidi M, Shahrokhi S, Al-Aasm I, Mahmmud R, Permanyer A. Meteoric Water Incursion, Crude Oil Degradation and Calcite Cementation of an Upper Cretaceous Reservoir in the Zagros Foreland Basin (Kurdistan Region of Iraq). Water. 2023; 15(10):1953. https://doi.org/10.3390/w15101953

Chicago/Turabian Style

Mansurbeg, Howri, Ibrahim Mohialdeen, Ali Al-Juboury, Namam Salih, Mohammad Alsuwaidi, Salahadin Shahrokhi, Ihsan Al-Aasm, Rebar Mahmmud, and Albert Permanyer. 2023. "Meteoric Water Incursion, Crude Oil Degradation and Calcite Cementation of an Upper Cretaceous Reservoir in the Zagros Foreland Basin (Kurdistan Region of Iraq)" Water 15, no. 10: 1953. https://doi.org/10.3390/w15101953

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop