ALBERT

Description: Physical and geochemical characteristics of produced petroleum from the central region of the Llanos basin, Colombia, were analyzed to understand the petroleum charge history and alteration processes. Petroleum properties in the study area are the result of the complex charge history of the oil fields. The amount of gas in fluids is controlled by the migration distance from the late or, possibly, the current generation kitchen located beneath the foothill deformation zone. Gas influx decreases toward the foreland domain, as indicated by lower values of the gas–oil ratio and saturation pressure. The API gravity of the oil samples is mainly controlled by the intensity of biodegradation. Marine-sourced oils accumulated in shallow reservoirs of the foreland prior to the onset of Andean deformation. Those fluids were subjected to different levels of biodegradation, depending on the time they remained at reservoir temperatures lower than 80°C (176°F) and before being buried to their maximum depth. Geochemical data suggest multiple charge pulses from different source kitchens of two main types of source rocks, as well as different biodegradation levels. The proposed petroleum charge and alteration model allows prediction of the temperature history of a reservoir and the most likely physical properties of the petroleum at a specific location. The model can be used as an exploration tool to assess the risk of charge prior to drilling in unexplored areas of the basin.

Description: In this paper, we demonstrate a workflow for constructing kinematic restorations in complex foothill areas devoid of growth strata and other indicators for the chronology of deformation. Our initial reconstructions utilize thermochronometric data, a well-documented structural geometry, and a first-order conversion of exhumation rates into tectonic rates. We then utilize models obtained from the new in-house–developed software FetKin to build a first version of the thermokinematic restoration. The FetKin approach is geared primarily toward testing and further calibration and refinement of the kinematic restoration, based on the extent to which the model result agrees with thermochronometric data from the study area in the form of both discrete ages and inverse-modeled time–temperature envelopes. This analysis also provides rates of shortening and time–temperature paths throughout the model space that can be used to make first-order predictions of when different source rocks entered the oil window. These capabilities are demonstrated in a pilot case study along a cross section in the Colombian Eastern Cordillera. The improved confidence in the reconstruction that this technique provides allows us to show increasing shortening rates in this part of the Andes during the Neogene reaching up to 5 mm/yr (0.20 in./yr) by the Pliocene, and constrain the timing of generation from the most important oil kitchens for the Eastern Cordillera-Llanos basin petroleum system. This approach, therefore, proves to be a useful method for creating high-resolution and high-fidelity kinematic restorations.

Description: Distinguishing axial and lateral sedimentary systems in rift basins is crucial for predicting reservoir distribution and quality, particularly where synrift strata are interrupted by mass transport complexes (MTCs). Upper Jurassic deep-marine synrift successions in the central North Sea have been studied to assess the temporal and spatial relationships of sediments and controls on reservoir quality. In the Late Jurassic, the central graben experienced erosion at rift margins, whereas adjacent grabens were starved and underfilled with marine sediments, supplied by axial and transverse systems. This study focused on sediments adjacent to a major intrabasinal high, the Josephine ridge. Data included seismic, wireline logs from 16 wells, and biostratigraphic and sedimentological analysis of 144 m (472 ft) of core. Synrift strata are dominated by mudstones but include MTCs interbedded with coarse sandstones at the rift margin and fine-grained turbidite sandstones in basinal depocenters. Petrographic and heavy mineral data indicate different provenance between MTCs and basinal turbidites. Turbidites correlate with periods of lowered relative sea level, during the initial rift phase, and record axial sediment supply. The composition of the MTCs corresponds to in situ strata on the adjacent Jade and Judy horsts. The distribution of MTCs implies formation by crestal collapse horsts during the rift climax and represents a transverse system, with no genetic relationship to axial turbidites. In starved deep-marine basins, fine-grained, well-sorted axial systems may provide the most extensive reservoirs. Transverse systems derived from isolated horsts are typically coarse-grained, poorly sorted, and spatially restricted, being unlikely to provide significant reservoir material.

Description: Significant amounts ( $$ 〉 150\hbox{ \hspace{0.17em} }\hbox{ \hspace{0.17em} }{\mathrm{m}}^{3}/\mathrm{day}/\mathrm{well}$$ ) of water are currently being extracted from coalbed methane (CBM) wells in Permian–Carboniferous coal in the Liulin area of the eastern Ordos basin, China. Waters coproduced with CBM have common chemical characteristics that can be an important exploration tool because they relate to the coal depositional environment and hydrodynamic maturation of groundwater and can be used to guide CBM development strategies. The CBM production targets of the No. 3 and 4 coal seams from sandstone in the Shanxi Formation and No. 8, 9, and 10 coal seams in the karst of the Taiyuan Formation were deposited in fluvial-deltaic and epicontinental-sea environments, respectively. This paper combines CBM geology, hydrogeology, CBM recovery, and laboratory data to define mechanisms of CBM preservation including the important influence of groundwater. Relevant indices include fluid inclusions as an indicator of the hydraulic connection between the coal seam reservoir and the overlaying strata and the ensemble characteristics of total dissolved solids (TDS) contents of water, water production rates, and reservoir temperatures as an indication of the current hydraulic connection. The TDS contents of waters from the No. 3 and 4 and No. 9 and 10 coal seams are double those from the subjacent karst No. 8 coal seam, indicating the important control of fast flow in karst. Low-salinity fluid inclusions from the roof of the subjacent-karst No. 8 coal seam also indicate an enduring hydraulic connection with overlaying strata during its burial history. Relatively low current temperatures in the No. 8 (subjacent-karst) coal seam also infer a strong hydraulic connection and active flow regime. Deuterium concentrations are elevated in the mudstone-bounded No. 9 and 10 coal seams, further confirming low rates of fluid transmission. The gas contents of coal seams from the Taiyuan Formation are higher than those from the sandstone-bounded coal seams in Shanxi Formation, also correlating with low rates of water transmission and low permeability. Conceptual models for these fluvial-deltaic and epicontinental-sea environments that are consistent with geology, gas content, and gas and water production rate histories are of gas-pressure sealing for the Shanxi Formation and hydrostatic-pressure sealing for the Taiyuan Formation. These results confirm the important controls of hydrogeological conditions on the preservation of CBM and the utility of hydrogeological indicators in prospecting for CBM.

Description: Using a seismic database from the Qiongdongnan Basin in the South China Sea, this study demonstrates that shelf-edge trajectories and stratal stacking patterns are reliable, but understated, predictors of deep-water sedimentation styles and volumes of deep-water sand deposits, assisting greatly in locating sand-rich environments and in developing a more predictive and dynamic stratigraphy. Three main types of shelf-edge trajectories and their associated stratal stacking patterns were recognized: (1) flat to slightly falling trajectories with negative trajectory angles ( $${T}_{\mathrm{se}}$$ ) (–2° to 0°) and negative shelf-edge aggradation to progradation ratios ( $$\mathrm{d}y/\mathrm{d}x$$ ) (–0.04 to 0) and associated progradational and downstepping stacking patterns with low clinoform relief ( $${R}_{\mathrm{c}}$$ ) (150–550 m [492–1804 ft]) and negative differential sedimentation on the shelf and basin ( $${A}_{\mathrm{s}}/{A}_{\mathrm{b}}$$ ) (–0.6 to 0); (2) slightly rising trajectories with moderate $${T}_{\mathrm{se}}$$ (0°–2°) and medium $$\mathrm{d}y/\mathrm{d}x$$ (0–0.04), and associated progradational and aggradational stacking patterns with intermediate $${R}_{\mathrm{c}}$$ (250–400 m [820–1312 ft]) and intermediate $${A}_{\mathrm{s}}/{A}_{\mathrm{b}}$$ (0–0.6); and (3) steeply rising trajectories with high $${T}_{\mathrm{se}}$$ (2°–6°) and high $$\mathrm{d}y/\mathrm{d}x$$ (0.04–0.10) and associated dominantly aggradational stacking patterns with high $${R}_{\mathrm{c}}$$ (350–650 m [1148–2132 ft]) and high $${A}_{\mathrm{s}}/{A}_{\mathrm{b}}$$ (1–2). Each trajectory regime represents a specific stratal stacking patterns, providing new tools to define a model-independent methodology for sequence stratigraphy. Flat to slightly falling shelf-edge trajectories and progradational and downstepping stacking patterns are empirically related to large-scale, sand-rich gravity flows and associated bigger and thicker sand-rich submarine fan systems. Slightly rising shelf-edge trajectories and progradational and aggradational stacking patterns are associated with mixed sand/mud gravity flows and moderate-scale slope-sand deposits. Steeply rising shelf-edge trajectories and dominantly aggradational stacking patterns are fronted by large-scale mass-wasting processes and associated areally extensive mass-transport systems. Therefore, given a constant sediment supply, then $${T}_{\mathrm{se}}$$ , $$\mathrm{d}y/\mathrm{d}x$$ , $${R}_{\mathrm{c}}$$ , and $${A}_{\mathrm{s}}/{A}_{\mathrm{b}}$$ are all proportional to intensity of mass-wasting processes and to amounts of mass-transport deposits, and are inversely proportional to the intensity of sand-rich gravity flows and to amounts of deep-water sandstone. These relationships can be employed to relate quantitative characteristics of shelf-edge trajectories and stratal stacking patterns to deep-water sedimentation styles.

Description: This multidisciplinary study evaluates the structural and hydrogeologic evolution of Cretaceous-age reservoirs in the Putumayo basin, Colombia. We focused on the Eastern Cordillera fold-thrust belt along the southern Garzón Massif. Many important hydrocarbon accumulations occurred regionally along the proximal foreland basin and frontal fold-thrust belt defining the eastern margin of the northern Andes. To understand why recent Putumayo basin and adjacent thrust belt exploration has resulted in a wide range of oil quality and limited economic discoveries, we reconstructed the structural evolution, timing of oil migration, and timing of groundwater infiltration by (1) assessing regional trends in formation water, oil, and reservoir properties; (2) quantifying the timing of hydrocarbon generation and migration relative to trap formation using (a) two-dimensional (2-D) and three-dimensional seismic data to define and constrain a restorable balanced cross section from the Upper Magdalena Valley to the Putumayo foreland and (b) coupled one-dimensional thermal basin modeling; (3) evaluating the potential roles of Mesozoic extensional faulting and Paleogene shortening in the generation and preservation of structural traps; and (4) assessing groundwater influx from the modern foothills into the reservoir using a 2-D numerical groundwater flow model. We suggest that four-way closure is limited in the study area, where most foreland-verging structures create three-way fault closures that do not effectively trap light hydrocarbons. In addition, east-dipping structures and a relatively large reservoir outcrop area provide water infiltration pathways. Groundwater modeling suggests reservoirs were water washed by 2–200 million pore volumes since Andean uplift. Finally, average reservoir temperatures are 〈80°C (〈176°F), which further facilitated biodegradation.

Description: Characterizing natural fracture systems involves understanding fracture types (faults, joints, and veins), patterns (orientations, sets, and spacing within sets), size distributions (penetration across layering, aperture, and trace length), and timing relationships. Traditionally, observation-based relationships to lithology, mechanical stratigraphy, bed thickness, structural position, failure mode, and stress history have been proposed for predicting fracture spacing along with the relative abundance of opening-mode fracture versus faults in fractured rocks. Developing a conceptual fracture model from these relationships can be a useful process to help predict deformation in a fractured reservoir or other fractured rock systems. A major pitfall when developing these models is using assumptions based on general relationships that are often site specific rather than universal. In this paper, we examine a mixed carbonate-shale sequence that is cut by a seismic-scale normal fault where fracture attributes do not follow commonly reported fracture relationships. Specifically, we find (1) no clear relationship between frequency (or spacing) of opening-mode fractures (joints and veins) and proximity to the main fault trace and (2) no detectable relationship between fracture spacing and bed thickness. However, we did find that (1) the frequency of small-displacement faults is strongly and positively correlated with proximity to the main fault trace, (2) fracture networks change pattern and failure mode (extension versus shear fracture) from pavement to pavement through the mechanically layered stratigraphic section, and (3) faults are more abundant than opening-mode fractures in many areas within the fracture network. We interpret that the major fracturing initiated near maximum burial under relatively high-differential stress conditions where shear failure dominated and that mode-1 extension fracturing occurred later under lower differential stress conditions, filling in between earlier formed shear fractures. We conclude that whenever possible, site-specific observations need to be carefully analyzed prior to developing fracture models and perhaps a different set of fracture network rules apply in rocks where shear failure dominates and mechanical stratigraphy influences deformation.

Description: Tidal heterolithic sandstones are commonly characterized by millimeter- to centimeter-scale intercalations of mudstone and sandstone. Consequently, their effective flow properties are poorly predicted by (1) data that do not sample a representative volume or (2) models that fail to capture the complex three-dimensional architecture of sandstone and mudstone layers. We present a modeling approach in which surfaces are used to represent all geologic heterogeneities that control the spatial distribution of reservoir rock properties (surface-based modeling). The workflow uses template surfaces to represent heterogeneities classified by geometry instead of length scale. The topology of the template surfaces is described mathematically by a small number of geometric input parameters, and models are constructed stochastically. The methodology has been applied to generate generic, three-dimensional minimodels (9 m 3 [~318 ft 3 ] volume) of cross-bedded heterolithic sandstones representing trough and tabular cross bedding with differing proportions of sandstone and mudstone, using conditioning data from two outcrop analogs from a tide-dominated deltaic deposit. The minimodels capture the cross-stratified architectures observed in outcrop and are suitable for flow simulation, allowing computation of effective permeability values for use in larger-scale models. We show that mudstone drapes in cross-bedded heterolithic sandstones significantly reduce effective permeability and also impart permeability anisotropy in the horizontal as well as vertical flow directions. The workflow can be used with subsurface data, supplemented by outcrop analog observations, to generate effective permeability values to be derived for use in larger-scale reservoir models. The methodology could be applied to the characterization and modeling of heterogeneities in other types of sandstone reservoirs.

Description: Tidal heterolithic sandstone reservoirs are heterogeneous at the submeter scale because of the ubiquitous presence of intercalated sandstone and mudstone laminae. Core-plug permeability measurements fail to sample a representative volume of this heterogeneity. Here, we investigate the impact of mudstone drape distribution on the effective permeability of heterolithic, cross-bedded tidal sandstones using three-dimensional, surface-based "minimodels" that capture the geometry of cross beds at an appropriate scale. The impact of seven geometric parameters has been determined: (1) mudstone fraction, (2) sandstone laminae thickness, (3) mudstone drape continuity, (4) toeset dip, (5) climb angle of foreset–toeset surfaces, (6) proportion of foresets to toesets, and (7) trough or tabular geometry of the cross beds. We begin by identifying a representative elementary volume of 1 m 3 (~35 ft 3 ), confirming that the model volume of 9 m 3 (~318 ft 3 ) yields representative permeability values. Effective permeability decreases as the mudstone fraction increases, and it is highly anisotropic: vertical permeability falls to approximately 0.5% of the sandstone permeability at a mudstone fraction of 25%, whereas the horizontal permeability falls to approximately 5% and approximately 50% of the sandstone value in the dip (across mudstone drapes) and strike (parallel to mudstone drapes) directions, respectively. Considerable spread exists around these values, because each parameter investigated can significantly impact effective permeability, with the impact depending upon the flow direction and mudstone fraction. The results yield improved estimates of effective permeability in heterolithic, cross-bedded sandstones, which can be used to populate reservoir-scale model grid blocks using estimates of mudstone fraction and geometrical parameters obtained from core and outcrop-analog data.

Description: The Vaca Muerta–Quintuco system (uppermost lower Tithonian–lower Valanginian) is a thick shallowing-upward sedimentary cycle consisting of dark bituminous shales, marlstones, limestones, and sandstones, cropping out in the Neuquén Basin, west–central Argentina. This paper analyzes three outcrop sections in Chos Malal area, northem Neuquén province. Detailed facies analysis allows us to differentiate six facies associations, representing basinal to proximal outer ramp facies of a homoclinal carbonate ramp system (Vaca Muerta Formation) and basinal to shoreface facies of a mixed carbonate–siliciclastic shelf system (Quintuco Formation), prograding westward from the eastern margin of the basin. Two sequence hierarchies were recognized: 5 composite depositional sequences (third order) and 15 high-frequency sequences (fourth order). Fluctuations in organic matter content within the Vaca Muerta Formation suggest a relationship with depositional sequences, finding the highest values associated with transgressive systems tract, whereas the transition to the Quintuco Formation shows a strong decrease in total organic carbon. The x-ray diffraction studies show an increase of clay minerals and quartz in the transgressive systems tract of the Vaca Muerta Formation and an increase in the content of calcite in highstand systems tracts. This pattern is reversed in the Quintuco Formation. Our sequence stratigraphic approach contributes to the understanding of the relationship between organic matter, clay minerals, facies, stacking pattern, and relative sea level changes in this exceptional shale oil and shale gas unconventional reservoir. This study may be helpful for a better postulate of petrophysical and geomechanical models for unconventional exploration.