SciELO - Scientific Electronic Library Online

 
vol.43 issue2Mapping of indirect indices of hydrocarbon manifestations in the Gulf of Guacanayabo, CubaGeochemistry and mineralogy of the La Pava mine, Muzo-Quípama: implications for emerald exploration in Colombia author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Boletín de Geología

Print version ISSN 0120-0283On-line version ISSN 2145-8553

Bol. geol. vol.43 no.2 Bucaramanga May/Aug. 2021  Epub Mar 15, 2022

https://doi.org/10.18273/revbol.v43n2-2021006 

Artículos científicos

Geochemistry and petrographic analysis of Los Cuervos and Molino formations, Cesar Sub-Basin (Colombia): Implications in the petroleum system evolution

Análisis geoquímico y petrográfico de las formaciones Los Cuervos y Molino, Subcuenca Cesar (Colombia): implicaciones en la evolución del sistema petrolero

Gladys Marcela Avendaño-Sánchez1  * 

Mario García-González2 

Luis Enrique Cruz-Guevara2 

Luis Felipe Cruz-Ceballos¹ 

1Grupo de Investigación en Geología Básica y Aplicada, Universidad Industrial de Santander, Bucaramanga, Colombia.

2Energy and Other Non-Renewable Resources Research Group, Universidad Industrial de Santander, Bucaramanga, Colombia. mgarciag@uis.edu.co; lecruz@uis.edu.co


Abstract

A geochemical characterization of Los Cuervos and Molino formations in the Cesar Sub-Basin was carried out using core samples obtained from the ANH-La Loma-1 Well. A total of 113 Rock-Eval pyrolysis analysis, total organic carbon (TOC), total sulphur content (TS) analysis, 13 vitrinite reflectance analysis (%Ro) and 30 thin-section petrographic analysis were performed. Based on these new data, it was possible to classify the quality of organic matter and the current thermal maturation of Los Cuervos and Molino formations. Additionally, a petrographic characterization of 30 samples allowed the correlation of the lithology with the geochemical results. Also, one-dimensional geochemical modelling was implemented in order to contribute to the knowledge of the evolution of the oil system in the Cesar Sub-Basin. The spatial distribution of the formations used in the modelling was obtained from 2 seismic lines two-way time. The results obtained indicate that Los Cuervos Formation presents TOC values from 0.29 to 66.55%, TS values from 0.02 to 11.29%, their organic matter consisted of type III kerogen which is consistent with an immature thermal maturation stage. In contrast, the Molino Formation presents TOC values from 0.23 to 2.28%, TS values from 0.001 to 1.39%, their organic matter consisted of type II/III kerogen this suggests an early entry to the oil window with a maximum pyrolysis temperature (Tmax) value of 442°C. The geochemical modelling tunes better with measured data from palaeo-geothermometers (%Ro and Tmax). The geochemical modelling shows that, between 60 - 40 million years ago (mya), the Cretaceous formations entered in the oil generation window and it is expected that, between 40 - 30 mya, the Lagunitas, Aguas Blancas, and La Luna formations will be at their peak of hydrocarbon generation.

Keywords: Los Cuervos Formation; Molino Formation; Geochemical modelling; Source rock; Cesar Sub-Basin; Cesar Ranchería Basin.

Resumen

Las muestras de núcleo obtenidas del Pozo ANH-La Loma-1 localizado en la Subcuenca Cesar permitieron la caracterización geoquímica de las formaciones Los Cuervos y Molino. Se realizaron un total de 113 análisis de pirólisis Rock-Eval, de carbono orgánico total (COT) y de contenido de azufre total (TS), 13 análisis de reflectancia de vitrinita (%Ro) y 30 análisis petrográficos de secciones delgadas. Basados en estos nuevos datos se pudo establecer la calidad de la materia orgánica y la condición actual de la maduración termal de las formaciones. Adicionalmente, la caracterización petrográfica de 30 muestras permitió la correlación de la litología con los resultados geoquímicos. También, se realizó un modelamiento geoquímico 1D con la finalidad de aportar en el conocimiento de la evolución del sistema petrolífero de la Subcuenca Cesar. La distribución espacial de las formaciones usada en el modelamiento fue obtenida de 2 líneas sísmicas 2D en tiempo. Los resultados obtenidos indican que la Formación Los Cuervos presenta valores de COT de 0,29-66,55%, valores de TS de 0,02-11,29%, su materia orgánica es de kerógeno tipo III el cual es consistente con un estado de madurez termal inmaduro. En contraste, la Formación Molino presenta valores de COT de 0,23-2,28%, valores de TS de 0,001-1,39%, su materia orgánica es de kerógeno tipo II/III, esto sugiere una temprana entrada en ventana de generación con un valor de Tmáx de 442°C. El modelo geoquímico calibra bien con los datos medidos de paleo-geotermómetros (%Ro y Tmáx) y muestra que entre 60 - 40 Ma las formaciones Cretácicas entraron en ventana de generación de hidrocarburos y se espera que entre 40-30 Ma las formaciones Lagunitas, Aguas Blancas y La Luna se encuentren en su pico de generación de hidrocarburos.

Palabras clave: Formación Los Cuervos; Formación Molino; Modelamiento geoquímico; Roca fuente; Subcuenca Cesar; Cuenca Cesar Ranchería.

Introduction

The study area lies at the south-eastern of the Cesar Sub-Basin, which is part of the Cesar-Ranchería Basin, located at the northern side of Colombia (Figure 1). The Cesar-Ranchería Basin borders to the east with the Perijá Range, to the west with the Sierra Nevada de Santa Marta, to the north with the main lineation of the Oca Fault and to the south with the trace of the Bucaramanga Fault (Barrero et al., 2007). The right side of Figure 1, the generalized column taken from Mesa and Rengifo (2011), shows the ages and thickness reported for each of the existing formations in the Cesar Sub-Basin. The stratigraphical/structural architecture, the main faults, truncated terminations, and onlap strata termination of the Cretaceous, Paleocene, and Neogene age formations are observed in the seismic lines, two-way time (Figure 2) as well as the stratigraphic position of ANH-La Loma-1 Well.

Figure 1 Location map, ANH-La Loma-1 Well, seismic lines used CV-1979-08, CV-1979-03, and wells with geochemical information. To the right, generalized stratigraphic column taken from Mesa and Rengifo (2011)

Multiple studies of potential hydrocarbon establish that this basin contains promising source rocks despite not having commercial hydrocarbon production (García et al., 2008; Hernández and Jaramillo, 2009; Mesa and Rengifo, 2011; Sánchez and Mann, 2015). Particularly, Martínez de Vivas et al. (2012) reported three maturation bands in the source rocks that increase towards the south-eastern of the Cesar-Rancheria Basin with potential for unconventional Liquid Gas type hydrocarbons. This study also highlights the subsurface “hydrocarbon shows” reported in wells within the clayey and calcareous formations of the Cretaceous, such as wet gas in the Cesar F-1 Well, light hydrocarbons in the Paso-3 and Paso-2 wells and commercial gas in the Compae Field Well.

Figure 2 Interpreted 2-D seismic lines in Two-Way Time (TWT). A. CV-1979-08. B. CV-1979-03, with the location of ANH-La Loma-1 Well. 

The geochemical properties of the Cesar Sub-Basin have been described in previous reports (Mesa and Rengifo, 2011) which indicate that the TOC values of the source rock vary from 0 to 5% in La Luna and Aguas Blancas formations, while Los Cuervos Formation presents TOC values of 16% in some mudstones. Aguas Blancas and La Luna formations have organic matter consisted of type II kerogen with capacity to generate liquid hydrocarbons, while Los Cuervos Formation has organic matter consisted of type III kerogen. Thermal values in La Luna, Aguas Blancas, and Lagunitas formations correspond from mature to postmature rocks (Mesa and Rengifo, 2011), with Tmax values from 430 to 475°C and 1.4 of %Ro. Likewise, Mora et al. (2007) processed 181 surface rock samples for TOC analysis and 81 surface samples for Rock-Eval pyrolysis analysis, concluding that La Luna and Ánimas formations must have mainly generated liquid hydrocarbons with gaseous in a lower proportion, while Los Cuervos Formation has favourable conditions for gas generation mainly because its organic matter consisted predominantly of type III kerogen. Correspondingly, Aguilera et al. (2010) reports that Lagunitas, Aguas Blancas, La Luna, and Molino formations have organic matter consisted of type II oil-prone kerogen, while Los Cuervos and Molino formations have organic matter consisted of type III kerogen. Furthermore, Aguas Blancas, La Luna, Molino, and Los Cuervos formations have reached early maturity to overmature conditions in the basin. Additionally, Hernández and Jaramillo (2009) processed 52 core samples for %Ro analysis in the Compae-1 Well for Aguas Blancas and Rio Negro formations and in the Diablito-1E Well for Molino Formation. Although only 7 %Ro measurements were conclusive, 4 of that data belong to the Molino Formation show values from 0.32 to 0.51%, 2 data from the Rio Negro Formation show values of 0.73 and 0.74%, and 1 data from the Aguas Blancas Formation show a value of 0.52%. These results indicate that the organic matter of the deeper samples would be in the oil window. Guo et al. (2018) made a geochemical characterization of sub-bituminous coals from Los Cuervos Formation, specifically from La Jagua and Calenturitas mines (Cesar-Ranchería Basin) and showed the results of %Ro measurements with values from 0.42 to 0.6%, TOC values from 68.5 to 86.0% and Rock-Eval pyrolysis analysis with Tmax values ranging from 419 to 435ºC, which allowed them to conclude that organic matter consisted predominantly of type III kerogen mixed with type II kerogen, the %Ro suggested that the studied coals are sub- bituminous from ranks B to A, which is confirmed by the Tmax values obtained.

This study presents new geochemical and petrographic analysis performed in core samples of source rocks from the ANH-La Loma-1 Well. The samples analysed belong to the Paleocene Los Cuervos Formation and the Upper Cretaceous Molino Formation. Geochemical data include the Rock-Eval pyrolysis of 113 samples, the TOC of 113 samples, and the vitrinite reflectance (%Ro) of 13 samples. Additionally, a petrographic characterization of 30 samples allowed the correlation of the lithology with the geochemical results. This study aims to establish the current maturity state and the potential for hydrocarbon generation (quality and quantity) of the rocks in the Los Cuervos and Molino formations. Furthermore, we would like to reproduce the burial history scenario through numerical modelling with the PetroMod® software (version 2012). The modelling is based on geochemical data, the stratigraphic/structural architecture, as interpreted by 2D seismic reflection data, and the different exhumation events of the Perijá Range (PR) reported by Kellogg (1984), Shagam et al. (1984) and Lugo and Mann (1995), that affected the burial history. It also provides insights that may help to understand the current maturity state, as well as the organic transformation ratios (TR) of other cretaceous formations underneath the Molino Formation. This is fundamental for any understanding of the organic geochemical characteristics of the Los Cuervos and Molino formations, as well as how the sub-basin was filled with sediments. It has direct implications regarding the maturation of organic matter of the units as well as the generation of hydrocarbons.

Geological setting and stratigraphy

During the early Cretaceous a basal conglomerate with gravels of metamorphic, igneous, and sedimentary composition of the Rio Negro Formation were deposited in paleo-depressions (Miller, 1962). This allows the configuration of a smooth topography that favours the deposition of carbonates towards the top of the formation (Gil and Martínez, 1990). Sarmiento-Rojas (2019) argue that the Rio Negro Formation was deposited in a fluvial to the transitional environment in an extensional rift setting, along the underlying unconformable Jurassic rocks from La Quinta Formation (Sánchez and Mann, 2015). Hernández (2003) indicated that the subsidence during the Late Aptian was continuous and that Lagunitas Formation was deposited upon an inner carbonate shelf. Lagunitas Formation corresponds to thick- bedded fossiliferous limestones and very thin-bedded calcareous shales (Sarmiento-Rojas, 2019). During the Albian to Cenomanian, Aguas Blancas Formation (Cesar-Ranchería Basin) and Maraca Formation (Maracaibo Basin) were deposited on an inner to middle marine shelf, controlled by relative sea-level changes (Sarmiento-Rojas, 2019). This transgression is evidenced by bio-micritic limestones and fossiliferous organic-rich shales, glauconitic sandstones, mudstones and sandy limestones (Mesa and Rengifo, 2011).

La Luna Formation was deposited during the Cenomanian to Coniacian on a passive margin setting (Sánchez and Mann, 2015). It corresponds to bituminous biomicrite (pelagic limestones), organic- rich calcareous mudstones, chert, concretions, and abundant content of planktonic and benthic foraminifera such as ammonites, and fish remains. It presumably occurred on an external part of the middle shelf (Mesa and Rengifo, 2011). During the Campanian and Maastrichtian, the Molino Formation was deposited in a deep marine environment that continued until the Middle Campanian (Hernández, 2003). Villamil (1999), determined that in the Maracaibo and Merida basins (Western Venezuelan) the Tres Esquinas Formation was deposited with the presence of glauconite and phosphates, which origin was from the Guyana Craton. Overlying the Tres Esquinas Member, in the area of Maracaibo and Merida basins (Western Venezuela), the Mito Juan and Colón formations were deposited upon an inner to the middle shelf (Sarmiento-Rojas, 2019). In addition, The Molino Formation has been correlated to the Colón - Mito Juan Formations of the Maracaibo Basin (Mesa and Rengifo, 2011).

During the Paleocene, the Barco and Los Cuervos formations were deposited in a continental environment as a product of a slow regression (Hernández, 2003). Barco Formation corresponds to quartz arenites with thin interlayered mudstones and coal (Mesa and Rengifo, 2011). During the Eocene, due to the Pre-Andean Orogeny, the accommodation space for the north of South America decreased significantly, accelerating the exhumation rate of the Central Mountain Range. Between the Eocene and Oligocene, the drainages of the Magdalena and Orinoco rivers took their current direction (Hoorn et al., 1995; Villamil, 1999). Throughout the Late Miocene-Pliocene, during the Andean Orogeny, the main uplift of the Serranía del Perijá was reported (Hernández, 2003). Finally, during the Pliocene-Pleistocene, fans, terraces and alluvial plains were formed covering the Cesar-Ranchería Valley (Hernández, 2003).

Data and methods

Geochemical analyses were performed to report the type of organic matter, current maturation stages, and potential to generate hydrocarbons in the Cesar Sub- Basin. Geochemical data were obtained using 113 core samples from Los Cuervos and Molino formations of ANH-La Loma-1 Well. The geochemical analyses include Rock-Eval pyrolysis analysis, TOC, and TS. Additionally, 13 vitrinite reflectance analysis (%Ro) were obtained through organic petrography. Also, 30 thin sections of samples rocks were described. Vertical sampling was not performed with symmetrical length values (Figure 3). The geochemical analysis and petrographic studies were carried out at the UIS Petroleum Geochemistry Laboratory (LGP) and UIS Petrographic Laboratory, respectively.

The petrographic studies revealed proportions of various minerals present and their spatial distributions in the rock that allowed the correlation of the petrography conditions with the geochemical results. The petrographic description method uses 300 equally spaced points across the thin section. The TOC analysis was carried out with the intention of quantifying the amount of organic matter present in the sample. To determine the percentage of TOC, the content of carbonates present in the sample was eliminated using HCl to dissolve inorganic CaCO3. After that, an automated analyzer Leco SC-144DR was used to establish the amount of TOC in each of the 113 rock samples which included both, kerogen, and bitumen. Together with the TOC, the mineral carbon content (MinC), the volume of free hydrocarbon generated before pyrolysis (S1), the amount of hydrocarbon generated during pyrolysis (S2), and the amount of emitted carbon dioxide (S3) were determined. It also allowed us to calculate the TS as well as the carbonate content. Likewise, to evaluate the quality of this organic matter, a Rock-Eval pyrolysis analysis was performed in each of the samples with Rock-Eval 6 technology equipment. This technique, presented by Barker (1974), Espitalie et al. (1977, 1985) and Behar et al. (2001), uses programmed heating steps of a small amount of rock (70 mg) or coal (30-50 mg) in an inert atmosphere (helium or nitrogen) in order to determine the quantity of free hydrocarbons present in the sample (S1 peak), as well as those that potentially can be released after maturation (S2 peak), and the amount of emitted carbon dioxide (S3). This procedure allows us to establish the Tmax and the current state of thermal maturation. Derivations of these parameters such as hydrogen index (HI), oxygen index (OI), production yield (PY) and production index (PI) were also used and calculated according to the following formulas (Law, 1999): HI = S2 / TOC × 100, OI = S3 / TOC × 100, PY = S1 + S2, and PI = S1 / S1 + S2. HI represents the amount of hydrogen relative to the amount of organic carbon present in a sample. OI represents the amount of oxygen relative to the amount of organic carbon present in a sample. The vitrinite reflectance (%Ro) was measured using Zeiss imager microscopy according to the method proposed by the ASTM D7708-14. The vitrinite reflectance analysis was carried out on core samples to evaluate the diagenetic changes of the organic matter (maturity level). The method starts after calibrating the equipment at an average of 42.61 %Ro, then readings were taken for each polished section. In the readings, the vitrinite macerals were identified on the polished surface, considering that the standard deviation of the 20 readings of %Ro per polished section must be less than 0.003. Geochemical information was plotted to visualize its distribution on profiles using free software PanPlot2 as described by Sieger and Grobe (2013). Six cross-plots, using geochemical data obtained, were made to illustrate and to interpret the results of the geochemical analysis. Interpretation of two 2D seismic lines (CV-1979-08, and CV-1979-03 provided by ANH) in two-way time (TWT) (Figure 2) were used in this research with the purpose of defining the cutting relationships between lithological units present in the basin. Additionally, they unveiled the geographic and stratigraphic position of the ANH-La Loma-1 Well used for this study.

Finally, a one-dimensional basin model was carried out with PetroMod software in order to reconstruct the burial history and thermal maturation of the geological units present in the well. This modelling used the vitrinite thermal maturation model (EASY %Ro) proposed by Sweeney and Burnham (1990), based on the Arrhenius first-order parallel-reaction approach with a distribution of activation energies. EASY %Ro was calibrated to a more rigorous model of vitrinite maturation based on the chemical properties of coal vitrinite. Previous studies, such as Mesa and Rengifo (2011), were also considered. Such studies provided the deposition ages of each unit, and the maximum and minimum thicknesses of the Rio Negro, Lagunitas, Aguas Blancas and lower part of the Molino formations. For Barco, Los Cuervos, and upper Molino formations we used the data obtained from ANH-La Loma-1 Well. The eroded thickness input for the modelling was inferred by the relation to the maximum thickness reported plus those found in the wells. The heat flow evolution was established by PetroMod software using the Mackenzie method, based on the sub-basin evolution proposed by Mesa and Rengifo (2011). It was established that the heat flow from the syn-rift stage goes from 200 to 116 mya and the post-rift stage from 116 to 20 mya, according to Mesa and Rengifo (2011). The Paleo-water depths (PWD) were estimated with biostratigraphic studies (e.g., Patarroyo et al., 2017). The PWD for La Luna and Molino formations, was estimated to be 0 m deep at 120 mya at the Rift stage according to Mesa and Rengifo (2011), which then increased to reach 100 m between 116-103 mya, then it increased to 200 m at 92 mya. Consequently, it was gradually decreased going through 100 m at 79 mya, after that 15m around 70 mya and remained around 0m to the present (Table 1). Sediment Water Interface Temperature (SWIT) was estimated using PetroMod software (e.g., Wygrala, 1989), which assumed that the current position of the study area is 11 degrees north latitude in the northernmost part of South America (Montes et al., 2019). To determine the evolution of the source rock units, the multiple kinetics based on lithologies were used (e.g., Di Primio and Horsfield, 2006). The IES_

TII_brown_limestone kinetic parameter was selected for Lagunitas, Aguas Blancas, La Luna, and Lower Molino formations, considering similarities between facies of these formations and the Brown Limestone Formation. For the middle Molino Formation, the kinetics parameter IES_TII_toarcian_shale was used since both formations present shale facies. For the Los Cuervos Formation, two different kinetic parameters were used, Handil_MahakamDelta_TIII kinetic for shales and IES_TIII_Tertiary_Coal_2C kinetic for coal facies.

Table 1 Paleo-batimetry conditions (PWD) vs time used to generate the geochemical modelling. 0 m deep during the early rift stage (120 mya) according to Mesa and Rengifo (2011). 

Age (mya) PWD (m)
0 0
50 0
70 15
79 100
92 200
103 100
116 100
120 0

Results

In this section we summarize the results obtained, corresponding to petrographic and geochemical analysis for ANH-La Loma-1 Well. Figure 3 shows a stratigraphic column with samples collected. Figures 4 and 5 show the geochemical results.

Petrographic results

The Molino Formation, in samples from LL1-1910.3-P to LL1-437.8-P (Figure 3), shows the presence of gastropods, trochospiral planktonic foraminifera, ostracods, phosphate fragments (bones), echinoderms with coating phosphate material, bryozoan shells and bivalve shells, bivalve fragments with lamellar shells, biserial benthic foraminifera shells, fragments and spicules of echinoderms, and phosphate fragments of teeth and carbonate recrystallization, with exception of some samples presenting angular quartz sandstone with straight extinction, phosphatic skeletal fragments, phosphatic intraclasts, reticulated algae, laminated ostracods fragments, bryozoan fragments, serpulid tubes, micro limestone matrix -granular and limestone cement; additionally, other samples present siltstone of angular and subangular quartz, finely laminated with flat-parallel structure, interstratified with organic material.

Figure 3 Lithologies, type of analysis, petrographic sample names, and main components presents in 30 thin sections from ANH-La Loma-1 Well. 

From Barco Formation, in samples from LL1-344.6-P to LL1-203.5-P (Figure 3), there was presence of siltstone of angular and subangular quartz with straight extinction and feldspar with polysynthetic and Carlsbad twinning, amorphous silica, phyllosilicates, some clastic muscovite, finely laminated interstratified with organic material and sandstone of angular to sub-angular quartz, some quartz were banded, indicating that they are of metamorphic origin, there was also presence of feldspars and chert fragments, and fragments of phosphate material.

The sample LL1-109.5-P from Los Cuervos Formation (Figure 3) is composed of 40% angular and subangular quartz with straight extinction, 3% of altered kaolin feldspars, chert fragments and clastic muscovite with less than 1%, clayey siliceous matrix, phosphate fragments, 7% of chert, clastic muscovite, some fragments of phosphate bone, some zircons and glauconite.

The petrographic analysis from the Diablito-1E and Los Cerezos-1E wells, extracted from the technical report of the Drummond LTD. (2007), describes that, for the top of the Molino Formation there is presence of biomicrite, fine-grained glauconitic quartz arenites and sublitharenite, for the lower member of the Barco Formation, there is presence of mudstones with siderite and fine-grained litharenites with calcareous cement, for the upper member of the Barco Formation there is medium to fine-grained porous litharenites and sublitharenites, and for the lower member of the Los Cuervos Formation there is fine-grained micaceous and argillaceous litharenites. In the exploratory drilling carried out for the technical report of the Drummond LTD. (2007), six samples of the Verdesia- 2E well core were subjected to palynological and micropaleontological analysis; from these samples, the samples taken from 254.3 to 254.4 ft and from 384.6 to 384.7 ft presented a good recovery of palynomorphs, which indicates marine deposit environment with a high continental contribution, located within the Echimonocolpites protofranciscoi zone from Maastrichtian age, identified as the Molino Formation. In the Diablito-1E Well that crossed the Barco-Los Cuervos succession, samples were taken every 5ft for biostratigraphy analysis. The results indicate that the succession formed by the Barco and Los Cuervos formations is of Paleocene age. In addition, for the La Luna-3E Well, this document reported that the Molino Formation was drilled and traverses a monotonous range of dark gray, laminar shales, with algae fragments and echinoderms. Then, there is a sequence of dark biomicrite (packstones), more concentrated towards the top, and towards the base wackestones and mudstones, with the presence of uniserial planktonic and benthic foraminifera; additionally, bivalve shells and pyritization of foraminifera were reported. Ayala-Calvo (2009) states that the Molino Formation has siliceous laminar mudstone facies, packstones/ grainstones facies of bioclasts with glauconite and phosphates. The Molino Formation presents a predominance of siliceous mudstones and decreasing grain bioturbated calcareous cycles and packstones with glauconite and phosphates.

Lithologies and petrographic analysis from the ANH- La Loma-1 Well, compared against the existing bibliography for the Cesar Sub-Basin, indicate that the well crosses a sequence that goes from the Molino Formation to the base from Late Cretaceous age, followed by the Barco Formation from Early Paleocene age, ending up at the top with the Los Cuervos Formation from the middle to late Paleocene age.

Geochemical data

The geochemical values obtained from the Los Cuervos and Barco Formations in ANH- La Loma-1 Well samples were: S1 peak, indicating the volume of free hydrocarbon generated before pyrolysis, which varied from 0.03 to 10.71 mg HC/g rock. S2 values, indicating the amount of hydrocarbon generated during pyrolysis, showing variations from 0.15 to 256.14 mg HC/g rock, and S3 values, indicating the amount of carbon dioxide emitted, showing variations from 0.16 to 2.71 mgCO2/g rock.

The organic richness of the Los Cuervos and Barco Formations ranged from poor to excellent, according to Lewan and Allison (1996). The Paleocene rocks (12.4 to 127m depth) shows TOC values that range from 0.29 to 66.65%. The OI and HI were the most important parameters, of which the correlation indicated the kerogen type of source rocks. The OI value, as calculated, varied from 23.28 to 384.3 mg CO2/g TOC. The HI value, as calculated, varied from 3.88 to 73.9 mg HC/g TOC.

The TS is an important consideration in coal utilization. The TS obtained varied from 0.022 to 11.291 % w/w. The PY was the total amount of hydrocarbon generated/ liberated under the pyrolysis S1 and S2 curves (Ghori, 1998). This index helped to categorize the target horizon in terms of total petroleum-generation potential (Hazra et al., 2019). The PY values, as calculated, varied from 0.21 to 266.85 mg HC/g.

The PI value was taken from a ratio between S1 and S2, which indicated the relationship between the number of free hydrocarbons before pyrolysis with respect to the amount of hydrocarbon generated during pyrolysis, similar to Rock-Eval and Tmax. PI was used as a maturity proxy (Hazra et al., 2019). The PI values, as calculated in ANH- La Loma-1 Well, varied from 0.01 to 0.29 TOC. According to Peters (1986), source rocks with PI values lower than 0.1 TOC are considered thermally immature, while those from 0.1 and 0.4 TOC are considered thermally mature. According to Nuñez-Betelu and Baceta (1994), there is no depletion in the values of these two pyrolysis peaks since the samples were taken from the core and had not been exposed to weathering.

The Tmax values correspond to the maximum hydrocarbon generation point during the pyrolysis (S2), which is widely used as a thermal maturity indicator (Hazra et al., 2019). According to Peters and Cassa (1994) any value of Tmax lower than 435ºC is considered thermally immature (Figure 4). Tmax values obtained ranged from 416 to 436ºC, which indicated that Paleocene rocks are thermally immature, bordering on the limit of 435ºC.

The geochemical values obtained from the Cretaceous rocks in ANH- La Loma-1 Well samples were: S1 peak, as obtained, which varied from 0.03 to 0.85 mg HC/g rock. S2 values, which showed variations from 0.14 to 2.68 mg HC/g rock, and S3 values, which ranged from 0.13 to 1.42 mgCO2/g rock.

The organic richness of the Cretaceous rocks ranged from poor to very good, according to Lewan and Allison (1996). The Cretaceous rocks (131 to 625.8m depth) showed TOC values that ranged from 0.23 to 2.28%. The OI value, as calculated, ranged from 17.48 to 239.05 mg CO2/g TOC. The HI value, as calculated, varied from 16.45 to 238.42 mg HC/g TOC. The TS values obtained ranged from 0.001 to 1.391% w/w. The PY values, as calculated, ranged from 0.18 to 2.8 mg HC/g. The PI ranged from 0.04 to 0.46 %TOC. The Tmax values ranged from 413 to 442ºC which indicates that some Cretaceous rocks are in the inmature oil window, and some Cretaceous rocks are generating oil, according to Peters and Cassa (1994), both of them bordering on the limit of 435ºC (Figure 4).

Figure 4 Graphic distribution of TOC, TS, and Rock-Eval pyrolysis from ANH-La Loma-1 Well core samples at the Cesar Sub-Basin. 

The vitrinite reflectance analysis of 13 core samples in ANH- La Loma-1 Well are presented in Figure 5. The %Ro values, as obtained, varied from 0.49 to 0.59, with a standard deviation (sd) ranging from 0.01 to 0.02. The TS values, as obtained, varied from 1.61 to 11.29%.

The %Ro values correction was made due to suppression of vitrinite generated by presence greater than 1% of TS. The data correlation of TS vs %Ro was performed to determine the reflectance of suppressed vitrinite, according to Mukhopadhyay (1994), who established that the suppression occurs between values of 0.2 to 1.4 %Ro due to excess sulphur of more than 1% w/w. According to Carr (2000), who established that in some situations, vitrinite reflectance values may be either ‘suppressed’ or ‘enhanced’.

Figure 5 Graphic distribution of %Ro, TOC, and TS from 13 core samples of the ANH-La Loma-1 Well at the Cesar Sub-Basin. 

Data analysis and discussion

The geochemical data obtained indicate that Cretaceous rocks are mostly thermally mature, according to Peters (1986). The low vitrinite reflectance of vitrinite suppression may be caused by perhydrous vitrinite, impregnation of vitrinite with hydrogen-rich amorphous material, and in coals with high sulphur content (Carr, 2000). The data correlation of %Ro vs Tmax is presented in Figure 6 to show the change in thermal maturity of Los Cuervos Formation in the La Loma-1Well. Figure 6 shows a good correlation between data. The data of thermal maturation (%Ro and Tmax) reaching values from 415 to 425ºC, show a lithostatic charge of at least 3 kilometres of quartz sandstones, after the deposition of Los Cuervos Formation. Therefore, it is suggested that there was a continuous accumulation of sediments in the sub-basin. This accumulation of sediments ended in the Middle Eocene, when the uplift of Serranía de Perijá began, according to García et al. (2008)

Figure 6 Data correlation of maximum pyrolysis temperature (Tmax) vs vitrinite reflectance (%Ro), to indicate source rock thermal maturity in the ANH-La Loma-1 Well rocks at the Cesar Sub-Basin. 

Quality of organic matter

The cross-plot of S1 vs TOC (Figure 7) was used to evaluate the origin of the hydrocarbon in the Cesar Sub-Basin, using examples of the ANH-La Loma- 1 Well, based on the methodology of Fakhri et al. (2013). It was established that the Cretaceous and Paleocene rocks were no contaminated with migrated hydrocarbons. The hydrocarbon generation potential was proven using the cross-plot of PY vs TOC (Figure 8), based on the proposition by Ghori and Haines (2007). It concludes that the Paleocene mudstones had from good to excellent potential, the Paleocene coals had excellent potential, the Cretaceous limestones had from poor to fair potential, the Cretaceous sandstones had poor potential, and the Cretaceous mudstones had from poor to fair potential.

Figure 7 Cross-plot of S1 vs TOC for ANH-La Loma-1 Well in the Cesar Sub-Basin. 

Figure 8 Cross-plot of Production Yield (PY) vs TOC for ANH-La Loma-1 Well in the Cesar Sub-Basin. 

To assess the kerogen type, the data obtained from ANH-La Loma-1 Well were plotted on HI vs OI diagram, and the modified Van Krevelen diagram (Espitalie et al., 1985) (Figure 9). It concludes that the Cretaceous mudstones and limestones have the type II (oil/gas prone) and III (gas prone) kerogen, and the Paleocene rocks predominantly have type III kerogen (gas prone).

To assess the maturity and kerogen type, the data obtained were plotted on HI vs Tmax diagram (e.g., Peters and Cassa, 1994) (Figure 10). This showed the source rocks with more evolved organic matter, higher Tmax values, more thermally mature stage (Barker, 1974; Espitalie, 1986). It was concluded that, in the Cesar Sub-Basin rocks, the maturation stages increase with age. Both, Paleocene mudstones and coals with values of Tmax from 400 to 425ºC, have not reached the oil generation window, but they would have the capacity to generate mainly gas-type hydrocarbons due to its organic matter type III kerogen. Both, Cretaceous- aged mudstones and limestones with values of Tmax from 425 to 450ºC, are in the gas generation due to its organic matter type III kerogen (gas prone).

These findings are in concordance with Mesa and Rengifo (2011), who affirmed that the Cretaceous sequence, Lagunitas, Aguas Blancas, La Luna, and Molino formations, show very good characteristics as source rocks due to its organic matter type II/ III kerogen, and TOC values higher than 1%. These findings agree with the results reported by Mora et al. (2007), who stated that Los Cuervos Formation has favourable conditions for the generation of gas due to its organic matter, mainly type III kerogen, and it can generate some liquid hydrocarbons due to its organic matter type II kerogen.

The cross plot of Tmax (ºC) vs PI as proposed by Katz (1995) was used to assess the maturation condition of the formations in the Cesar Sub-Basin (Figure 11). It can be seen that the Paleocene mudstone, sandstones and coal samples (with the exception of one case) occupy the immature area, which reinforce the immaturity of the Paleocene rocks of the sub- basin. Most of the Cretaceous mudstone, sandstones and limestone samples occupy the mature area. Tmax (ºC) data vs depth (Figure 12) as done by Peters and Cassa (1994), showed that there was no repetition of the sedimentary sequences in the ANH-La Loma-1 Well, and that there was no presence of displacement due to geological faults; a situation confirmed by the information obtained from the seismic interpretation, as well as from the drilling report.

Figure 9 Van Krevelen diagram, cross-plot of hydrogen index vs oxygen index of ANH-La Loma-1 Well in the Cesar Sub-Basin. 

Figure 10 Cross-plot between Tmax (ºC) vs the hydrogen index (HI) of ANH-La Loma-1 Well in the Cesar Sub-Basin. 

Figure 11 Cross-plot between production index (PI) vs pyrolysis Tmax (ºC) to identify source rock maturity for the ANH-La Loma-1 Well (Cesar Sub-Basin). 

Figure 12 Cross-plot between the depths (meters) vs Tmax (ºC) to demonstrate the distribution of the parameter Tmax with increasing depth in the Los Cuervos Formation of the ANH-La Loma-1 Well (Cesar Sub-Basin). 

Geochemical Modelling

A one-dimensional model of burial history was performed in the location of the ANH-La Loma-1 Well using the information obtained from the interpretation of the seismic lines in time, the boundary conditions heat flow, the paleo-water depth, the SWIT (Figure 13) and the bottom hole temperature (BHT) of the wells (Figure 15B).

The geological input parameters used in the modelling of burial history are presented in Table 2. Initially, the Lower Cretaceous Rio Negro Formation has 220 m of continental conglomerates and mudstones intercalation deposited from 120 to 116 mya. It is overlain by 371 m of limestones and marine/littoral mudstones intercalation from Lagunitas Formation from 116 to 103 mya. Then, 198 m of limestones and marine/ littoral mudstones and sandstones from Aguas Blancas Formation deposited from 103 to 92 mya. From 92 to 79 mya, 340 m intercalation of limestones, cherts, and marine/littoral mudstones with some calcareous concretions from La Luna Formation were deposited. Then, 521 m of marine/littoral mudstones, limestones and sandstones from Molino Formation was deposited from 79 to 65.5 mya. Around 65.5 mya, a large erosive event occurred that eroded 100 m from the Molino Formation. From 63 to 55 mya, 998 m of the non- marine sandstones, claystones and siltstones of Barco and Los Cuervos formations were deposited, with 78 and 920 m thickness, respectively. Then, during the erosional event of the Middle Eocene, 900m of the Los Cuervos Formation were eroded. From 48.6 to 37.2 mya, 3200 m of the Eocene Sedimentites Unit were deposited. Then, from 37 to 21 mya, in the location of ANH-La Loma-1 Well, an erosive event occurred related to the uplift of the Serranía de Perijá eroding the Eocene Sedimentites Unit. From 16 to 12 mya, 200m of the Cuesta Formation was deposited, and finally, from 12 mya to the present the whole record of this formation was eroded.

Table 2 Geological input parameters used for one-dimensional modelling of ANH-La Loma-1 Well. 

Layer Thick. (m) Eroded (m) Depo. from (mya) Depo. to (mya) Eroded from (mya) Eroded to (mya) Lithology (mya)
Cuesta Fm 0 200 16 12 12 0 0 Conglomerate (typical)
Eocene sedimentites 0 2600 48.6 37.2 37.2 21.79 Sandstone (typical)
Los Cuervos Fm. 7 900 61.53 55.8 21.79 16 Shale (typical)
Los Cuervos Fm. [Paleosol] 8 - 61.58 61.53 - - Sandstone (typical)
Los Cuervos Fm. [Coal] 1 1 - 61.59 61.58 - - Coal (pure)
Los Cuervos Fm. [Mudstone] 13 - 61.66 61.59 - - Shale (organic rich, 20% TOC)
Los Cuervos Fm. 6 - 61.71 61.66 - - Siltstone (organic lean)
Barco Fm. [Sandstone] 23 - 62.09 61.71 - - Sandstone (typical)
Barco Fm. [Siltstone] 6 - 62.19 62.09 - - Siltstone (organic lean)
Barco Fm. [Mudstone] 49 - 63 62.19 - - Shale (organic rich, 1% TOC)
Upper Molino Fm. 0 100 69.3 65.5 65.5 63 Shale (typical)
Middle Molino Fm. 17 - 69.6 69.3 - - Sandstone (typical)
Middle Molino Fm. [Mudstone] 44 - 70.5 69.6 - - Shale (organic rich 1-2%, typical)
Middle Molino Fm. 134 - 72.8 70.5 - - Siltstone (organic lean)
Middle Molino Fm. [Limestone] 7 - 72.9 72.8 - - Limestone (organic rich - typical)
Lower Molino Fm. [Limestone] 30 - 73.5 72.9 - - Limestone (organic rich 1-2% TOC)
Lower Molino Fm. 159 - 76.5 73.5 - - Siltstone (organic lean)
Lower Molino Fm. 130 - 79 76.5 - - Shale (organic rich, typical)
Upper La Luna Fm. 57 - 82 79 - - Chert
Middle La Luna Fm. 185 - 88 82 - - Shale (organic rich, 8% TOC)
Lower La Luna Fm. 98 - 92 88 - - Limestone (Chalk, typical)
Upper Aguas Blancas Fm. 66 66 - 95 92 - - Limestone (Chalk, typical)
Middle Aguas Blancas Fm. 66 - 100 95 - - Sandstone (typical)
Lower Aguas Blancas Fm. 66 - 103 100 - - Shale (organic rich, 2% TOC)
Upper Lagunitas Fm. 166 - 107 103 - - Limestone (Chalk, typical)
Middle Lagunitas Fm. 70 - 110 107 - - Limestone (micrite)
Lower Lagunitas Fm. 135 - 116 110 - - Limestone (Chalk, typical)
Upper Rio Negro Fm 80 - 117 116 - - Conglomerate (typical)
Middle Rio Negro Fm 60 - 119 117 - - Siltstone (organic lean)
Lower Rio Negro Fm 80 - 120 119 - - Conglomerate (typical)
120

The Heat Flow (HF) values were calculated using the Mckenzie crystal model option (McKenzie, 1978), considering that the heat flow for the Cesar Sub-Basin is conditioned by local tectonic dynamics. The syn-rift phase goes from 200 to 116 mya (Mesa and Rengifo, 2011) and ends with the deposition of the Rio Negro Formation; the greatest heat flow occurred during this phase, reaching 57.88 mW⁄m2 (milliWatt per square meter) during the Lower Cretaceous. After this, the basin entered the post- rift stage and gradually cooled down until it reached the current heat flow of approximately 34.65 mW⁄m2 (Figure 13). The EASY %Ro was calibrated with paleo- geothermometer data (%Ro and Tmax) using samples extracted from the Los Cuervos Formation (Figure 14A and 14B). Figure 14A shows an excellent fit of the values modelled and measured for %Ro and Tmax.

Figure 13 Boundary conditions: A. Heat flow. B. Paleo water depth and SWIT, used for 1D modelling of burial history of the lithological units in the ANH-La Loma-1 Well. 

Figure 14 A. Calibration of measured %Ro and Tmax with depth in ANH-La Loma-1 Well. B. Calibration of measured %Ro and Tmax in-depth projection. 

The model of burial history (Figure 15), shows that, around 60 mya, the Lagunitas Formation entered in the early oil window. Towards 45 mya, it approached the peak of the increase in values of %Ro. Around 60 mya, the Aguas Blancas Formation entered in the early oil window, and around 40 mya, it approached the peak of hydrocarbon generation, according to the %Ro. Around 40 mya, La Luna Formation approached the oil window, and around 30 mya, it approached the hydrocarbon generation peak, according to the %Ro. Finally, from 40 mya to 0, the Molino Formation entered into the early oil window. According to the vitrinite reflectance diagram (Figure 14B and 15) and, it is expected that the Molino Formation and upper part of the La Luna Formation would be in the early oil window since before the Upper Eocene uplift event. While Aguas Blancas and Lagunitas formations would be at the oil window since the Late Paleocene. After this time, it is expected that its thermal maturation did not progress any further. Additionally, the Los Cuervos Formation is immature in the well. The higher maturation levels in the ANH-La Loma-1 Well might have been reached for the other source rock of Cretaceous formations.

Figure 15 A. Burial history with respect to the increase in %Ro in the time expected for ANH-La Loma-1 Well. This figure is an in-depth projection of ANH-La Loma-1 Well. B. BHT modelled of ANH-La Loma-1 and Cesar A-1X wells. 

The thermal stress can be measured using the Tmax and %Ro values, which is necessary to evaluate the amount of organic matter transformed into hydrocarbon (Pepper and Corvi, 1995). This has could be to explain the scatter in the transformation rates (TR) observed in Figure 16. As a result, the TR were not expected in the study area for Los Cuervos Formation. The TR estimated for the Cretaceous Lagunillas, Aguas Blanca, La Luna, and Molino formations in the ANR-La Loma-1 Well had values of 94, 89, 68, and 12%, respectively. The Lower Cretaceous formations have a greater advance in their TR value at the end of the Paleocene, while the Upper Cretaceous units had almost all their TR increase during the Eocene. This establishes two pulses of hydrocarbon generation in the Paleocene and the Eocene, respectively. This shows that the study area in the Los Cuervos Formation is not in the oil window.

Figure 16 Hydrocarbon transformation ratios from ANH-La Loma-1 Well. 

Conclusions

The mudstone of the Los Cuervos Formation, with 15 samples analysed, showed 20.72% TOC value on average and the presence of organic matter consisting of mixed type III kerogen. It supported a good source rock potential to produce gas type hydrocarbons.

The coals of the Los Cuervos Formation, with 2 sample analysed, presented an average of 60.16% TOC value and the presence of organic matter consisting of type III kerogen. The intervals indicated evidence of excellent gas-generating hydrocarbon.

The mudstones of the Barco Formation, with 8 sample analysed, presented an average of 1.61% TOC value and the presence of organic matter consisting of type III kerogen. It supported a good source rock potential to produce gas type hydrocarbons.

The limestone of the Molino Formation, with 14 sample analysed, presented an average of 0.84% TOC value and the presence of organic matter consisting of type II/III kerogen. Evidence of poor to fair source rock potential to produce oil/gas type hydrocarbons.

The mudstone of the Molino Formation, with 70 samples analysed, presented an average of 0.57% TOC value and the presence of organic matter type II/III kerogen. It supported a poor source rock potential to produce oil/gas type hydrocarbons.

The Los Cuervos Formation in the ANH-La Loma- 1 Well is in the immature stage, but close to the oil window. The Molino Formation in the ANH-La Loma- 1 Well is in the mature stage, but close to the main oil window.

According to the modelling carried out to the ANH- La Loma-1 Well, the Cretaceous formations (Rio Negro, Lagunitas, Aguas Blancas, La Luna, Molino) entered the hydrocarbon generation phase between 60 and 40 mya. The Lagunitas, Aguas Blancas, and La Luna formations approached the peak of hydrocarbon generation between 40 and 30 mya.

Similarly, through the modelling and with %Ro and Tmax data, it is concluded that La Luna formation is in the generation phase and has transformation ratios of 68%, while Lagunitas and Aguas Blancas formations are in the generation phase and have transformation rates of 89% and 94%, respectively.

The petrographic analysis indicates that the Molino Formation was deposited in a marine environment under oxidizing conditions, indicated by the presence of mainly benthic microfossils. Although, the environment was marine, the potential for preservation of organic matter is low because there were no anoxic conditions at the bottom. Petrographic analysis indicate that the rocks from the Molino Formation are fine-grained rocks with hydrocarbon potential, but geochemical data indicate that the Molino Formation is not proven to be a good source rock.

Acknowledgments

The Agencia Nacional de Hidrocarburos (ANH) from Colombia and the Universidad Industrial de Santander (UIS) have granted access and approved the use of all data for this study. Special thanks for providing the resources and facilities for the development of this research. Thanks to the Petroleum Geochemistry Laboratory (LGP) at the UIS for performing the geochemical analysis for this work. Thanks also to the style-writer corrector, German Morales.

References

Aguilera, R.; Sotelo, V.; Burgos, C.; Arce, C.; Gómez, C.; Mojica, J.; Castillo, H.; Jiménez, D.; Osorno, J. (2010). Organic Geochemistry Atlas of Colombia. 2nd Edition. Vol. 14. Special Edition. Earth Sciences Research Journal. [ Links ]

ASTM D7708-14. Standard Test Method for Microscopical Determination of the Reflectance of Vitrinite Dispersed in Sedimentary Rocks. ASTM International, West Conshohocken, PA, 2014. https://doi.org/10.1520/D7708-14Links ]

Ayala-Calvo, R.C. (2009).Análisis tectonoestratigráfico y de procedencia en la Subcuenca de Cesar: relación con los sistemas petroleros. MSc. Tesis, Universidad Simón Bolívar, Sartenejas, Venezuela. [ Links ]

Barker, C. (1974). Pyrolysis techniques for source- rock evaluation. AAPG Bulletin, 58(11), 2349-2361. https://doi.org/10.1306/83D91BAF-16C7-11D7-8645000102C1865DLinks ]

Barrero, D.; Pardo, A.; Vargas, C.A.; Martínez,J. (2007). Colombian sedimentary basins: nomenclature, boundaries and petroleum geology, a new proposal. ANH. [ Links ]

Behar, F.; Beaumont, V.; De Penteado, B. (2001). Rock-Eval 6 technology: performances and developments. Oil & Gas Science and Technology, 56(2), 111-134. https://doi.org/10.2516/ogst:2001013Links ]

Carr, A.D. (2000). Suppression and retardation of vitrinite reflectance, part 1. Formation and significance for hydrocarbon generation. Journal of Petroleum Geology, 23(3), 313-343. https://doi.org/10.1111/j.1747-5457.2000.tb01022.xLinks ]

Di Primio, R.; Horsfield, B. (2006). From petroleum- type organofacies to hydrocarbon phase prediction. AAPG Bulletin, 90(7), 1031-1058. https://doi.org/10.1306/02140605129Links ]

Drummond LTD. (2007). Resumen de algunos aspectos estructurales y estratigráficos del Bloque La Loma. Informe técnico. [ Links ]

Espitalie, J. (1986). Use of Tmax as a maturation index for different types of organic matter. Comparison with vitrinite reflectance. In: J. Burrus (ed). Thermal modelling in sedimentary basins (pp. 475-496). Editions Technip. [ Links ]

Espitalie, J.; Deroo, G.; Marquis, F. (1985). La pyrolyse Rock-Eval et ses applications. Première partie. Revue Institut Francais Petrole, 40(5), 563-579. https://doi.org/10.2516/ogst:1985035Links ]

Espitalie, J.; Laporte, J.L.; Madec, M.; Marquis, F.; Leplat, P.; Paulet, J.; Boutefeu, A. (1977). Méthode rapide de caractérisation des roches mères, de leur potentiel pétrolier et de leur degré d’évolution. Revue Institut Francais Petrole, 32(1), 23-45. https://doi.org/10.2516/ogst:1977002Links ]

Fakhri, M.; Tabatabaei, H.; Amiri, A. (2013). Comparing the potential of hydrocarbon generation of Kazhdomi and Pabdeh formations in Bangestan Anticline (Zagros Basin) according to Rock-Eval pyrolysis data. Journal Earth Science Climate Change, 4(5), 157. https://doi.org/10.4172/2157-7617.1000157Links ]

García, M.; Mier, R.; Arias, A.; Cortes, Y.; Moreno, M.; Salazar, O.; Jiménez, M. (2008). Prospectividad de la Cuenca Cesar-Ranchería. Grupo de Investigación en Geología de Hidrocarburos y Carbones, Universidad Industrial de Santander- ANH, Colombia. [ Links ]

Gil, A.; Martínez, H. (1990). Análisis micro facial del Grupo Cogollo, y Formación La Luna, cuenca Cesar Ranchería, Guajira, Colombia. Tesis, Universidad Nacional de Colombia, Bogotá. [ Links ]

Ghori, K.A.R. (1998). Petroleum generating potential and thermal history of the Palaeozoic, Carnarvon Basin, Western Australia. In: P.G. Purcell, R.R. Purcell (eds). The Sedimentary Basins of Western Australia (pp. 553-568). Vol. 2. Petroleum Exploration Society of Australia Limited. [ Links ]

Ghori, K.A.R.; Haines, P.W. (2007). Paleozoic petroleum systems of the Canning Basin, Western Australia: A review. Search and Discovery, 5-8. [ Links ]

Guo, Q.; Littke, R.; Zieger, L. (2018). Petrographical and geochemical characterization of sub- bituminous coals from mines in the Cesar- Ranchería Basin, Colombia. International Journal of Coal Geology, 191, 66-79. https://doi.org/10.1016/j.coal.2018.03.008Links ]

Hazra, B.; Wood, D.A.; Mani, D.; Singh, P.K.; Singh, A.K. (2019). Evaluation of Shale Source Rocks and Reservoirs. Springer. [ Links ]

Hernández, M. (2003). Memoria explicativa geología, plancha 48, Jagua de Ibirico. Escala 1:100.000, INGEOMINAS. Bogotá, Colombia. [ Links ]

Hernández, O.; Jaramillo, J.M. (2009). Reconstrucción de la historia termal en los sectores de Luruaco y Cerro Cansona-cuenca del Sinú-San Jacinto y en el piedemonte occidental de la Serranía del Perijá entre Codazzi y la Jagua de Ibirico- Cuenca de Cesar-Ranchería. Agencia Nacional de Hidrocarburos, Bogotá, Colombia. [ Links ]

Hoorn, C.; Guerrero, J.; Sarmiento, G.A.; Lorente, M.A. (1995). Andean tectonics as a cause for changing drainage patterns in Miocene northern South America. Geology, 23(3), 237-240. https://doi.org/10.1130/0091-7613(1995)023<0237:ATA ACF>2.3.CO;2Links ]

Katz, B.J. (1995). The Green River Shale: An Eocene carbonate lacustrine source rock. In: B.J. Katz (ed.). Petroleum Source Rocks (pp. 309-324). Springer. [ Links ]

Kellogg, J.N. (1984). Cenozoic tectonic history of the Sierra de Perijá, Venezuela-Colombia, and adjacent basin. In: W.E. Bonini; R.B. Hargraves; R. Shagam (eds.). The Caribbean-South American Plate Boundary and Regional Tectonics (pp. 239- 261). Geological Survey of America, Memoir 162. [ Links ]

Law, C.A. (1999). Evaluating source rocks. In: E.A. Beaumont; N.H. Foster (eds.). Treatise of Petroleum Geology / Handbook of Petroleum Geology: Exploring for Oil and Gas Traps (pp. 6-41). Chapter 6. AAPG Special Volumes. [ Links ]

Lewan, M.D.; Allison A.H. (1996). Gas: oil ratios for source rocks containing Type-I, II, IIS, and III kerogens as determined by hydrous pyrolysis. In: S.T. Dyman; Kuuskraa (eds.). Geologic Studies of Deep Natural Gas Resources (pp. 1-9). Chapter E. U.S. Geological Survey. [ Links ]

Lugo, J.; Mann, P. (1995). Jurassic-Eocene tectonic evolution of Maracaibo Basin, Venezuela. In: A.J. Tankard; R. Suárez S.; H.J. Welsink (eds.). Petroleum basins of South America (pp. 699- 725). AAPG Memoir 62. [ Links ]

Martínez de Vivas, M.; Calderón, W.; Zamora, W.; Rodríguez, I. (2012). Modelamientos numéricos 3D de sistemas petrolíferos en la Cuenca Cesar - Ranchería: nuevas ideas acerca de su potencial petrolífero. 11th Simposio Bolivariano - Exploración Petrolera en las Cuencas Subandinas. Cartagena, Colombia. [ Links ]

McKenzie, D. (1978). Some remarks on the formation of sedimentary basins. Earth Planetary Science Letters, 40, 25-32. [ Links ]

Mesa, A.M.; Rengifo, S. (2011). Geology and hydrocarbon potential regional geology of Colombia. Cesar Ranchería Basin. Vol. 6. ANH- University EAFIT, Medellín, Colombia. [ Links ]

Miller, J.B. (1962). Tectonic trends in Sierra de Perijá and adjacent parts of Venezuela and Colombia. AAPG Bulletin, 46(9), 1565-1595. https://doi.org/10.1306/BC7438D3-16BE-11D7- 8645000102C1865DLinks ]

Montes, C.; Rodríguez-Corcho, A.F.; Bayona, G.; Hoyos, N.; Zapata, S.; Cardona, A. (2019). Continental margin response to multiple arc- continent collisions: The northern Andes Caribbean margin. Earth-Science Reviews, 198. https://doi.org/10.1016/j.earscirev.2019.102903Links ]

Mora, C.; Parra, P.; Navas, G. (2007). Caracterización geoquímica de rocas y crudos en las cuencas de Cesar-Ranchería, Sinú-San Jacinto, Chocó y área de Soápaga (Cuenca Cordillera Oriental). Informe Final. GEMS-ANH. [ Links ]

Mukhopadhyay, P.K. (1994). Vitrinite reflectance as maturity parameter, petrographic and molecular characterization and its applications to basin modelling. In: P.K. Mukhopadhyay; W.G. Dow (eds). Vitrinite Reflectance as a Maturity Parameter (pp. 1-24). Chapter 1. ACS symposium series no. 570. https://doi.org/10.1021/bk-1994- 0570.ch001Links ]

Nuñez-Betelu, L.; Baceta, J.I. (1994). Basics and application of Rock-Eval/TOC pyrolysis: An example from the uppermost Paleocene/ lowermost Eocene in the Basque Basin, Western Pyrenees. Munibe. Ciencias Naturales, 46, 43-62. [ Links ]

Patarroyo, G.D.; Torres, G.A.; Rincón, D.A.; Cárdenas, C.P.; Márquez, R.E. (2017). Bioestratigrafía e inferencias paleoambiéntales de las asociaciones de foraminíferos en las formaciones Cretácicas La Luna-Colón (Cuenca del Catatumbo, Colombia). Boletín de Geología, 39(3), 25-40. https://doi.org/10.18273/revbol.v39n3-2017002Links ]

Patiño, A.M.; Parra, M.; Ramírez, J.C.; Sobel, E.R.; Glodny, J.; Almendral, A.; Echeverri, S. (2019). Thermochronological constraints on Cenozoic exhumation along the southern Caribbean: The Santa Marta range, northern Colombia. In: B. Horton; A. Folguera (eds). Andean Tectonics (pp. 103-132). Chapter 5. Elsevier Inc. https://doi. org/10.1016/B978-0-12-816009-1.00007-1Links ]

Pepper, A.S.; Corvi, P.J. (1995). Simple kinetic models of petroleum formation. Part 1: oil and gas generation from kerogen. Marine and Petroleum Geology, 12(3), 291-319. https://doi.org/10.1016/0264-8172(95)98381-ELinks ]

Peters, K.E. (1986). Guidelines for evaluating petroleum source rock using programmed pyrolysis. AAPG Bulletin, 70(3), 318-329. https://doi.org/10.1306/94885688-1704-11D7- 8645000102C1865DLinks ]

Peters, K.E.; Cassa, M.R. (1994). Applied source-rock geochemistry. In: L.B. Magoon; W.G. Dow (eds.). The petroleum system from source to trap (pp. 93- 120). AAPG Memoir 60. [ Links ]

Sánchez, J.; Mann, P. (2015). Integrated structural and basinal analysis of the Cesar-Ranchería Basin, Colombia: Implications for its tectonic history and petroleum systems. In: J. Sánchez; P. Mann (eds). Petroleum geology and potential of the Colombian Caribbean margin (pp. 431-470). Chapter 16. AAPG Memoir 108. https://doi. org/10.1306/13531945M1083648Links ]

Sarmiento-Rojas, L.F. (2019). Cretaceous stratigraphy and paleo-facies maps of northwestern South America. In: F. Cediel; R.P. Shaw (eds.). Geology and Tectonics of Northwestern South America (pp. 673-747). Chapter 10. Springer. https://doi.org/10.1007/978-3-319-76132-9_10Links ]

Shagam, R.; Kohn, B.P.; Banks, P.O.; Dasch, L.E.; Vargas, R.; Rodriguez, G.I.; Pimentel. N. (1984). Tectonic implications of Cretaceous-Pliocene fission-track ages from rocks of the circum- Maracaibo Basin region of western Venezuela and eastern Colombia. In: W.E. Bonini; R.B. Hargraves; R. Shagam (eds.). The Caribbean- South American Plate Boundary and Regional Tectonics (pp. 385-412). Geological Survey of America, Memoir 162. https://doi.org/10.1130/MEM162-p385Links ]

Sieger, R.; Grobe, H. (2013). PanPlot 2 - software to visualize profiles and time series. Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research: Bremerhaven, Germany. Data Publisher for Earth & Environmental Science. https://doi.org/10.1594/PANGAEA.816201Links ]

Sweeney, J.J.; Burnham, A.K. (1990). Evaluation of a simple model of vitrinite reflectance based on chemical kinetics. AAPG Bulletin, 74(10), 1559-1570. [ Links ]

Villamil, T. (1999). Campanian-Miocene tectonostratigraphy, depocenter evolution and basin development of Colombia and western Venezuela. Palaeogeography, Palaeoclimatology, Palaeoecology, 153(1-4), 239-275. https://doi.org/10.1016/S0031-0182(99)00075-9Links ]

Wygrala, B.P. (1989). Integrated study of an oil field in the southern Po basin, Northern Italy. Thesis, Berichte Kernforschungsanlage Jülich, Germany. [ Links ]

Received: May 04, 2020; Accepted: March 08, 2021

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License